Sevillano Lucascantinelli M
Sevillano Lucascantinelli M
Sevillano Lucascantinelli M
CAMPINAS
2013
i
ii
UNIVER
RSIDADE ESTADUAL DE CAMPINAS
S
FACULA
ADADE DE ENGENHARIA MECÂNICA
A
E INSTITUTO DE GEOCIÊNCIAS
LU
UCAS CANTINELLI SEVILLANO
ESTUDO D
DE MOVIMENTAÇÃO DE UMA
U
PLATAFORMA DDE PERFURAÇÃO FLUTUA ANTE COM
BOP
P SUSPENSO PELO RISER
_______________________________________
Orientador
CAMPINAS
2013
iii
Ficha catalográfica
Universidade Estadual de Campinas
Biblioteca da Área de Engenharia e Arquitetura
Rose Meire da Silva - CRB 8/5974
Título em outro idioma: Study of the movement of a floating drilling platform with the BOP
suspended by the riser
Palavras-chave em inglês:
Maritime structures - Dynamics
Petroleum engineering
Drilling riser
Área de concentração: Explotação
Titulação: Mestre em Ciências e Engenharia de Petróleo
Banca examinadora:
Celso Kazuyuki Morooka [Orientador]
Paulo Roberto Ribeiro
Heitor Rodrigues de Paula Lima
Data de defesa: 29-08-2013
Programa de Pós-Graduação: Ciências e Engenharia de Petróleo
iv
UNIVERSIDADE ESTADUAL DE CAMPINAS
FACULADADE DE ENGENHARIA MECÂNICA
E INSTITUTO DE GEOCIÊNCIAS
______________________________________________
Prof. Dr. Celso Kazuyuki Morooka, Presidente
DEP/FEM/UNICAMP
______________________________________________
Prof. Dr. Paulo Roberto Ribeiro
DEP/FEM/UNICAMP
______________________________________________
Dr. Heitor Rodrigues de Paula Lima
PETROBRAS
v
vi
DEDICATÓRIA
Dedico este trabalho à minha família. Aos meus pais, Mathias e Lenita e ao meu irmão
Daniel, pelo amor, incentivo e motivação ao longo de toda minha vida.
vii
viii
AGRADECIMENTOS
Agradeço aos meus pais e a todos os membros da minha família pelo apoio dado durante
toda a minha formação.
Agradeço ao Prof. Dr. Celso K. Morooka pela oportunidade oferecida de trabalhar sob sua
tutela. Agradeço ao Prof. Dr. José R. P. Mendes pela coorientação. Agradeço ao Dr. Kazuo Miura
da Petrobras pelo apoio e interesse demonstrado neste trabalho.
Agradeço a todos os amigos que fiz durante meu tempo na UNICAMP e que me
acompanharam em grande parte do período na universidade. Em especial, aos meus colegas de
laboratório, Bruno, Raphael, Denis, Maurício, Michele, Renato, Luis, Danilo, Marcus, Natália,
Maiara e Pedro.
ix
x
“Perseverança, segredo de todos os
triunfos”
Victor Hugo
xi
xii
RESUMO
xiii
xiv
ABSTRACT
The location of newly found oilfields under ultra deepwater makes the running and
retrieving of submarine equipment lengthy operations. These operations can be executed with a
marine drilling riser. Regarding the particular case of a BOP, its installation at the wellhead is
fundamental for drilling the well. When this phase ends, the BOP is brought up to the surface and
the drilling rig leaves to drill a new well elsewhere. This study proposes to consider not retrieving
the BOP and, instead, letting it hanging of the drilling riser, during the change of location to drill
a new well. Thus, the present work has for objectives to first analyze riser behavior during this
operation, through computational simulations, while it faces several different cases of
environmental loading, and then propose a procedure that determines operational limits that must
be fulfilled, in order to preserve the riser’s physical integrity.
xv
xvi
SUMÁRIO
1. INTRODUÇÃO ................................................................................................. 1
2. DESENVOLVIMENTO ................................................................................... 7
3. METODOLOGIA ........................................................................................... 25
4. ESTUDOS DE CASO...................................................................................... 41
5. RESULTADOS ................................................................................................ 53
xvii
5.3. Mapas Operacionais ....................................................................................... 59
6. CONCLUSÕES ............................................................................................... 89
xviii
LISTA DE FIGURAS
Figura 1-1 – Tendências na exploração de petróleo offshore (Mahoney et al, 2012) ..................... 2
Figura 1-2 Linha de regressão para tempo de descida e recolhimento de BOP/LMRP em função
da profundidade de lâmina d’água (Holand, 2001) ......................................................................... 3
Figura 2-1- Riser Marítimo de Perfuração e componentes (adaptado de API RP 16Q, 1993) ....... 9
Figura 2-2 – Vista expandida de Aranha e Gimbal Fonte: Aker Solutions, 2010 ......................... 10
Figura 2-3- Discretização de riser (A), Elemento de viga 2D (B) ................................................. 11
Figura 2-4- Elemento diferencial de riser ..................................................................................... 12
Figura 2-5-Elemento de viga 3D ................................................................................................... 13
Figura 2-6-Carregamentos atuantes sobre o riser.......................................................................... 16
Figura 2-7 – Espectro de mar irregular.......................................................................................... 18
Figura 2-8 – Graus de liberdade de uma embarcação ................................................................... 19
Figura 2-9 – Transformação de ondas em resposta ....................................................................... 21
Figura 2-10 Perfil de Corrente Marítima. ...................................................................................... 23
Figura 3-1 Determinação do comprimento do riser. Adaptado de API RP 16Q (1993). .............. 26
Figura 3-2 – Tensões locais em elemento da parede do riser ....................................................... 30
Figura 3-3 – Esquema do contato do riser com estruturas abaixo do deck de perguração ........... 33
Figura 3-4 – Exemplo de histórico de tensão de von Mises. ......................................................... 33
Figura 3-5 - Tensão de von Mises máxima no riser em função de parâmetros de ondas ............. 35
Figura 3-6 – Curvas de nível de tensão de von Mises máxima no riser em função de parâmetros
de onda........................................................................................................................................... 35
Figura 3-7 – Mapa operacional gerado a partir do critério de projeto de tensão em função de
parâmetros de onda ........................................................................................................................ 37
Figura 3-8 – Curvas de nível de dupla amplitude significativa de heave da embarcação em função
de parâmetros de onda ................................................................................................................... 37
Figura 3-9 – Curvas de nível de tensão de von Mises máxima no riser em função de parâmetros
de heave ......................................................................................................................................... 39
Figura 3-10 – Mapa operacional gerado a partir do critério de projeto de tensão em função de
parâmetros de heave ...................................................................................................................... 39
xix
Figura 4-1 – Configuração do riser, em operação (A) e em navegação (B) ................................. 42
Figura 4-2 – Funções de transferências da embarcação semissubmersível: Surge, Sway e Heave 47
Figura 4-3 – Funções de transferências da embarcação semissubmersível: Roll, Pitch e Yaw ..... 48
Figura 4-4 – Carregamentos ambientais e direções de avanço da embarcação ............................. 51
Figura 5-1 – Amostras de resultados de simulações numéricas .................................................... 53
Figura 5-2 – Comparação de envelopes de tração efetiva ............................................................. 54
Figura 5-3 – Curvas de envoltórias de deslocamento horizontal (A) e deslocamento vertical (B).
kθ = 686 kNm/grau, v = 0,5 nó, θ = 180°, H1/3= 5 m e Tz = 13 s.................................................. 56
Figura 5-4 – Curvas de envoltórias de rotação no plano XZ (A) e de força axial efetiva (B). kθ =
686 kNm/grau, v = 0,5 nó, θ = 180°, H1/3= 5 m e Tz = 13 s .......................................................... 57
Figura 5-5 – Curvas de envoltórias de momento fletor no plano XZ (A) e de tensão de von Mises
(B). kθ = 686 kNm/grau, v = 0,5 nó, θ = 180°, H1/3= 5 m e Tz = 13 s ......................................... 58
Figura 5-6 – Curvas de nível de tensão de von Mises máxima no riser; kθ = 686 kNm/grau,
v = 0,5 nó, θ = 180° ....................................................................................................................... 60
Figura 5-7 - Curvas de nível de inclinação máxima no topo do riser; kθ = 686 kNm/grau,
v = 0,5 nó, θ = 180° ....................................................................................................................... 61
Figura 5-8 – Curvas de nível de tensão de von Mises máxima no riser em função dos parâmetros
de heave; kθ = 686 kNm/grau, v = 0,5 nó, θ = 180° ..................................................................... 62
Figura 5-9 – Curvas de nível de rotação no topo do riser em função dos parâmetros de heave; kθ
= 686 kNm/grau, v = 0,5 nó, θ = 180° .......................................................................................... 63
Figura 5-10 – Mapa operacional para direção de avanço da embarcação de 180 graus, gimbal de
rigidez rotacional de 676 kNm/grau. Parâmetros de onda............................................................. 65
Figura 5-11 – Mapa operacional para direção de avanço da embarcação de 180 graus, gimbal de
rigidez rotacional de 676 kNm/grau. Parâmetros de heave. .......................................................... 66
Figura 5-12 – Mapa operacional para direção de avanço da embarcação de 180 graus, gimbal de
rigidez rotacional de 338 kNm/grau. Parâmetros de onda............................................................. 67
Figura 5-13 – Mapa operacional para direção de avanço da embarcação de 180 graus, gimbal de
rigidez rotacional de 338 kNm/grau. Parâmetros de heave. .......................................................... 67
Figura 5-14 – Mapa operacional para direção de avanço da embarcação de 180 graus, gimbal de
rigidez rotacional de 169 kNm/grau. Parâmetros de onda............................................................. 68
xx
Figura 5-15 – Mapa operacional para direção de avanço da embarcação de 180 graus, gimbal de
rigidez rotacional de 169 kNm/grau. Parâmetros de heave. .......................................................... 68
Figura 5-16 – Mapa operacional para direção de avanço da embarcação de 135 graus, gimbal de
rigidez rotacional de 676 kNm/grau. Parâmetros de onda............................................................. 69
Figura 5-17 – Mapa operacional para direção de avanço da embarcação de 135 graus, gimbal de
rigidez rotacional de 676 kNm/grau. Parâmetros de heave. .......................................................... 69
Figura 5-18 – Mapa operacional para direção de avanço da embarcação de 135 graus, gimbal de
rigidez rotacional de 338 kNm/grau. Parâmetros de onda............................................................. 70
Figura 5-19 – Mapa operacional para direção de avanço da embarcação de 135 graus, gimbal de
rigidez rotacional de 338 kNm/grau. Parâmetros de heave. .......................................................... 70
Figura 5-20 – Mapa operacional para direção de avanço da embarcação de 135 graus, gimbal de
rigidez rotacional de 169 kNm/grau. Parâmetros de onda............................................................. 71
Figura 5-21 – Mapa operacional para direção de avanço da embarcação de 135 graus, gimbal de
rigidez rotacional de 169 kNm/grau. Parâmetros de heave. .......................................................... 71
Figura 5-22 – Mapa operacional para direção de avanço da embarcação de 90 graus, gimbal de
rigidez rotacional de 676 kNm/grau. Parâmetros de onda............................................................. 72
Figura 5-23 – Mapa operacional para direção de avanço da embarcação de 90 graus, gimbal de
rigidez rotacional de 676 kNm/grau. Parâmetros de heave. .......................................................... 72
Figura 5-24 – Mapa operacional para direção de avanço da embarcação de 90 graus, gimbal de
rigidez rotacional de 338 kNm/grau. Parâmetros de onda............................................................. 73
Figura 5-25 – Mapa operacional para direção de avanço da embarcação de 90 graus, gimbal de
rigidez rotacional de 338 kNm/grau. Parâmetros de heave. .......................................................... 73
Figura 5-26 – Mapa operacional para direção de avanço da embarcação de 90 graus, gimbal de
rigidez rotacional de 169 kNm/grau. Parâmetros de onda............................................................. 74
Figura 5-27 – Mapa operacional para direção de avanço da embarcação de 90 graus, gimbal de
rigidez rotacional de 169 kNm/grau. Parâmetros de heave. .......................................................... 74
Figura 5-28 – Mapa operacional para direção de avanço da embarcação de 45 graus, gimbal de
rigidez rotacional de 676 kNm/grau. Parâmetros de onda............................................................. 75
Figura 5-29 – Mapa operacional para direção de avanço da embarcação de 45 graus, gimbal de
rigidez rotacional de 676 kNm/grau. Parâmetros de heave. .......................................................... 75
xxi
Figura 5-30 – Mapa operacional para direção de avanço da embarcação de 45 graus, gimbal de
rigidez rotacional de 338 kNm/grau. Parâmetros de onda............................................................. 76
Figura 5-31 – Mapa operacional para direção de avanço da embarcação de 45 graus, gimbal de
rigidez rotacional de 338 kNm/grau. Parâmetros de heave. .......................................................... 76
Figura 5-32 – Mapa operacional para direção de avanço da embarcação de 45 graus, gimbal de
rigidez rotacional de 169 kNm/grau. Parâmetros de onda............................................................. 77
Figura 5-33 – Mapa operacional para direção de avanço da embarcação de 45 graus, gimbal de
rigidez rotacional de 169 kNm/grau. Parâmetros de heave. .......................................................... 77
Figura 5-34 – Mapa operacional para direção de avanço da embarcação de 0 grau, gimbal de
rigidez rotacional de 676 kNm/grau. Parâmetros de onda............................................................. 78
Figura 5-35 – Mapa operacional para direção de avanço da embarcação de 0 grau, gimbal de
rigidez rotacional de 676 kNm/grau. Parâmetros de heave. .......................................................... 78
Figura 5-36 – Mapa operacional para direção de avanço da embarcação de 0 grau, gimbal de
rigidez rotacional de 338 kNm/grau. Parâmetros de onda............................................................. 79
Figura 5-37 – Mapa operacional para direção de avanço da embarcação de 0 grau, gimbal de
rigidez rotacional de 338 kNm/grau. Parâmetros de heave. .......................................................... 79
Figura 5-38 – Mapa operacional para direção de avanço da embarcação de 0 grau, gimbal de
rigidez rotacional de 169 kNm/grau. Parâmetros de onda............................................................. 80
Figura 5-39 – Mapa operacional para direção de avanço da embarcação de 0 grau, gimbal de
rigidez rotacional de 169 kNm/grau. Parâmetros de heave. .......................................................... 80
Figura 5-40 – Curvas de nível de força axial efetiva máximo no riser. ........................................ 86
Figura 5-41 – Curvas de nível de momento fletor máximo no riser. ............................................ 87
Figura A-1 - Diagrama de corpo livre de um segmento tubular ................................................. 135
Figura B 1 – Comparação entre efeitos de carregamentos regulares no envelope de deslocamento
na direção inline........................................................................................................................... 105
Figura B 2 – Comparação entre efeitos de carregamentos regulares no envelope de deslocamento
na direção cross-flow ................................................................................................................... 105
Figura B 3 – Comparação entre efeitos de carregamentos regulares no envelope de deslocamento
de direção axial ............................................................................................................................ 106
Figura B 4 – Comparação entre efeitos de carregamentos regulares no envelope de rotação no
plano definido pelas direções vertical e inline............................................................................. 106
xxii
Figura B 5 – Comparação entre efeitos de carregamentos regulares no envelope de rotação no
plano definido pelas direções vertical cross-flow........................................................................ 107
Figura B 6 – Comparação entre efeitos de carregamentos regulares no envelope de força na
direção inline ............................................................................................................................... 107
Figura B 7 – Comparação entre efeitos de carregamentos regulares no envelope de força na
direção cross-flow........................................................................................................................ 108
Figura B 8 – Comparação entre efeitos de carregamentos regulares no envelope de força axial 108
Figura B 9 – Comparação entre efeitos de carregamentos regulares no envelope de momento
fletor no plano definido pelas direções vertical e inline .............................................................. 109
Figura B 10 – Comparação entre efeitos de carregamentos regulares no envelope de momento
fletor no plano definido pelas direções vertical e cross-flow ...................................................... 109
Figura B 11 – Comparação entre efeitos de carregamentos regulares na tensão devido ao
momento fletor ............................................................................................................................ 110
Figura B 12 – Comparação entre efeitos de carregamentos regulares na tensão devido à força
axial ............................................................................................................................................. 110
Figura B 13 – Comparação entre efeitos de carregamentos regulares no envelope de tensão
cisalhante ..................................................................................................................................... 111
Figura B 14 – Comparação entre efeitos de carregamentos regulares no envelope de tensão de
von Mises..................................................................................................................................... 111
Figura B 15 – Comparação do efeito do acréscimo de VIV no envelope de deslocamento na
direção inline ............................................................................................................................... 112
Figura B 16 – Comparação do efeito do acréscimo de VIV no envelope de deslocamento na
direção cross-flow........................................................................................................................ 112
Figura B 17 – Comparação do efeito do acréscimo de VIV no envelope de deslocamento na
direção axial................................................................................................................................. 113
Figura B 18 – Comparação do efeito do acréscimo de VIV no envelope de rotação no plano
definido pelas direções vertical inline ......................................................................................... 113
Figura B 19 – Comparação do efeito do acréscimo de VIV no envelope de rotação no plano
definido pelas direções vertical e cross-flow ............................................................................... 114
Figura B 20 – Comparação do efeito do acréscimo de VIV no envelope de força na direção inline
..................................................................................................................................................... 114
xxiii
Figura B 21 – Comparação do efeito do acréscimo de VIV no envelope de força na direção cross-
flow .............................................................................................................................................. 115
Figura B 22 – Comparação do efeito do acréscimo de VIV no envelope de força axial ............ 115
Figura B 23 – Comparação do efeito do acréscimo de VIV no envelope de momento fletor no
plano definido pelas direções vertical e inline............................................................................. 116
Figura B 24 – Comparação do efeito do acréscimo de VIV no envelope de momento fletor no
plano definido pelas direções vertical e cross-flow ..................................................................... 116
Figura B 25 – Comparação do efeito do acréscimo de VIV na tensão devido ao momento fletor
..................................................................................................................................................... 117
Figura B 26 – Comparação do efeito do acréscimo de VIV no envelope de tensão devido à força
axial ............................................................................................................................................. 117
Figura B 27 – Comparação do efeito do acréscimo de VIV no envelope de tensão de
cisalhamento ................................................................................................................................ 118
Figura B 28 – Comparação do efeito do acréscimo de VIV no envelope de tensão de von Mises
..................................................................................................................................................... 118
Figura C 1- Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 686 kNm/grau, velocidade de avanço de 0,5 nó, direção de 135 graus ..................... 119
Figura C 2- Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 686 kNm/grau, velocidade de avanço de 0,5 nó, direção de 90 graus ....................... 119
Figura C 3 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 686 kNm/grau, velocidade de avanço de 0,5 nó, direção de 45 graus ....................... 120
Figura C 4 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 686 kNm/grau, velocidade de avanço de 0,5 nó, direção de 0 graus ......................... 120
Figura C 5 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 686 kNm/grau, velocidade de avanço de 1,0 nó, direção de 180 graus ..................... 120
Figura C 6 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 686 kNm/grau, velocidade de avanço de 1,0 nó, direção de 135 graus ..................... 121
Figura C 7 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 686 kNm/grau, velocidade de avanço de 1,0 nó, direção de 90 graus ....................... 121
Figura C 8 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 686 kNm/grau, velocidade de avanço de 1,0 nó, direção de 45 graus ....................... 121
xxiv
Figura C 9 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 686 kNm/grau, velocidade de avanço de 1,0 nó, direção de 0 graus ......................... 122
Figura C 10 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 686 kNm/grau, velocidade de avanço de 1,5 nó, direção de 45 graus ...................... 122
Figura C 11 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 338 kNm/grau, velocidade de avanço de 0,5 nó, direção de 180 graus ..................... 122
Figura C 12 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 338 kNm/grau, velocidade de avanço de 0,5 nó, direção de 135 graus ..................... 123
Figura C 13 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 338 kNm/grau, velocidade de avanço de 0,5 nó, direção de 90 graus ....................... 123
Figura C 14 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 338 kNm/grau, velocidade de avanço de 0,5 nó, direção de 45 graus ....................... 123
Figura C 15 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 338 kNm/grau, velocidade de avanço de 0,5 nó, direção de 0 graus ......................... 124
Figura C 16 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 338 kNm/grau, velocidade de avanço de 1,0 nó, direção de 180 graus ..................... 124
Figura C 17 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 338 kNm/grau, velocidade de avanço de 1,0 nó, direção de 135 graus ..................... 124
Figura C 18 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 338 kNm/grau, velocidade de avanço de 1,0 nó, direção de 90 graus ...................... 125
Figura C 19 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 338 kNm/grau, velocidade de avanço de 1,0 nó, direção de 45 graus ....................... 125
Figura C 20 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 338 kNm/grau, velocidade de avanço de 1,0 nó, direção de 0 graus ......................... 125
Figura C 21 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 338 kNm/grau, velocidade de avanço de 1,5 nó, direção de 180 graus ..................... 126
Figura C 22 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 338 kNm/grau, velocidade de avanço de 1,5 nó, direção de 135 graus ..................... 126
Figura C 23- Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 338 kNm/grau, velocidade de avanço de 1,5 nó, direção de 90 graus ....................... 126
xxv
Figura C 24 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 338 kNm/grau, velocidade de avanço de 1,5 nó, direção de 45 graus ....................... 127
Figura C 25 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 338 kNm/grau, velocidade de avanço de 1,5 nó, direção de 0 graus ......................... 127
Figura C 26 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 169 kNm/grau, velocidade de avanço de 0,5 nó, direção de 180 graus ..................... 127
Figura C 27 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 169 kNm/grau, velocidade de avanço de 0,5 nó, direção de 135 graus .................... 128
Figura C 28 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 169 kNm/grau, velocidade de avanço de 0,5 nó, direção de 90 graus ....................... 128
Figura C 29 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 169 kNm/grau, velocidade de avanço de 0,5 nó, direção de 45 graus ....................... 128
Figura C 30 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 169 kNm/grau, velocidade de avanço de 0,5 nó, direção de 0 graus ......................... 129
Figura C 31 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 169 kNm/grau, velocidade de avanço de 1,0 nó, direção de 180 graus ..................... 129
Figura C 32 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 169 kNm/grau, velocidade de avanço de 1,0 nó, direção de 135 graus ..................... 129
Figura C 33- Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 169 kNm/grau, velocidade de avanço de 1,0 nó, direção de 90 graus ...................... 130
Figura C 34 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 169 kNm/grau, velocidade de avanço de 1,0 nó, direção de 45 graus ....................... 130
Figura C 35 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 169 kNm/grau, velocidade de avanço de 1,0 nó, direção de 0 graus ......................... 130
Figura C 36 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 169 kNm/grau, velocidade de avanço de 1,5 nó, direção de 180 graus ..................... 131
Figura C 37- Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 169 kNm/grau, velocidade de avanço de 1,5 nó, direção de 135 graus ..................... 131
Figura C 38 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 169 kNm/grau, velocidade de avanço de 1,5 nó, direção de 90 graus ....................... 131
xxvi
Figura C 39 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 169 kNm/grau, velocidade de avanço de 1,5 nó, direção de 45 graus ....................... 132
Figura C 40 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do
riser; kθ = 169 kNm/grau, velocidade de avanço de 1,5 nó, direção de 0 graus ......................... 132
Figura D-1 – Mapa Operacional – Parâmetros de Onda. Altura significativa de onda máxima em
função de comprimento de riser e de velocidade de navegação. ................................................ 133
Figura D-2 – Mapa Operacional – Parâmetros de Heave. Dupla amplitude significativa de heave
máxima em função de comprimento de riser e de velocidade de navegação.............................. 134
xxvii
xxviii
LISTA DE TABELAS
xxix
xxx
LISTA DE NOMENCLATURAS
Letras latinas
c Velocidade de onda
{d} Vetor deslocamento dos nós do riser discretizado
{ } Vetor velocidade dos nós do riser discretizado
{ } Vetor aceleração dos nós do riser discretizado
dm Peso específico do fluido de perfuração
dw Peso específico da água
fbt Fator de perda de flutuabilidade dos flutuadores
fwt Fator de tolerância aplicado ao peso submerso do riser
k Número de onda
[k] Matriz de rigidez do elemento
[ke] Matriz de rigidez elástica do elemento
[kg] Matriz de rigidez geométrica do elemento
kθ Rigidez rotacional
mn Momento espectral de ordem n
n Número de tracionadores que podem falhar de maneira súbita
ri raio interno do tubo
ro raio externo do tubo
t Tempo de simulação
u Velocidade (aceleração) da partícula de água
v Velocidade de avanço da embarcação
Velocidade (aceleração) do riser
z Heave instantâneo
z1/3 Dupla amplitude significativa de heave
A Área de seção do riser
Ai Área definida pelo raio interno do riser
Ao Área definida pelo raio externo do riser
xxxi
[B] Matriz de amortecimento global do sistema
Bn Força de empuxo devido aos módulos de flutuabilidade
CA Coeficiente de massa adicionada
CD Coeficiente de arrasto
D0 Diâmetro externo do riser
E Módulo de Elasticidade
{F} Vetor que representa as forças que atuam sobre os nós do riser
F Força hidrodinâmica por unidade de comprimento
Fc Força de arrasto devido à corrente marítima
Fn Fator normalizante do espectro de onda de Jonswap
FX Força axial sobre o riser
Fx0 Força horizontal devido à pressão externa
Fxi Força horizontal devido à pressão interna
Fy0 Força vertical devido à pressão externa
Fyi Força vertical devido à pressão interna
H1/3 Altura significativa de onda
Hm Altura da lâmina d’água
Hw Altura da coluna de fluido de perfuração
I Momento de inércia
IXY Momento de inércia ao redor dos eixos X e Y
[K] Matriz de rigidez global do sistema
[M] Matriz de massa global do sistema
MX Momento fletor ao redor do eixo X
MY Momento fletor ao redor do eixo Y
N Número de tracionadores em operação
P Pressão
Rf Fator de eficiência do sistema de tracionamento
Sζ Espectro de onda
SJ Espectro de onda de Jonswap
SPM Espectro de onda de Pierson-Moskowitz
Sz Espectro da resposta de heave
xxxii
T Tração
Tmin Tração mínima no topo do riser
TSRmin Tração mínima no anel de tracionamento
Tp Período de pico
Tz Período de zero ascendente
Uc Velocidade de corrente
V Força de cisalhamento
Ws Peso submerso do riser
Letras Gregas
xxxiii
Subscritos
0 Externo
i Interno
s Relativo à parede do riser
Siglas
xxxiv
1. INTRODUÇÃO
1.1. Motivação
Em diversas partes do globo, a busca por petróleo em alto-mar tem levado à exploração de
reservatórios de hidrocarbonetos em lâminas d’água de grande profundidade, como os campos do
pré-sal descobertos na bacia de Santos. Neste contexto, dois dos equipamentos essenciais para a
construção dos poços que darão acesso à acumulação de petróleo na subsuperfície são o riser e o
BOP.
O riser é a tubulação que interliga a cabeça de poço à embarcação. De acordo com seu
propósito ele pode ser classificado como riser de produção, de perfuração, de completação ou de
injeção.
Risers de produção servem como um meio para escoar os hidrocarbonetos produzidos pelo
reservatório, enquanto risers de injeção escoam fluidos que são injetados para o interior do
reservatório. Tanto um quanto o outro pode ser do tipo rígido ou flexível. Eles podem adotar uma
configuração vertical ou em catenária.
1
instalação, ou trazido de volta à superfície quando a etapa de perfuração termina e o poço é
abandonado temporariamente.
A Figura 1-1 apresenta a progressão ao longo dos anos da profundidade das lâminas d’água
na exploração offshore. A classificação adotada para profundidade de lâminas d’água é a
seguinte: até 400 metros (aproximadamente 1300 pés) têm-se águas rasas. Profundidades entre
400 e 1500 metros (aproximadamente entre 1300 e 5000 pés) são chamadas de águas profundas.
A partir de 1500 metros (ou aproximadamente 5000 pés) passa-se a operar em águas
ultraprofundas. Observa-se que o limite de águas rasas foi superado no final da década de 1960,
enquanto que o limite de águas profundas foi superado em meados da década de 1980. O recorde
observado em 2012 foi de perfurar-se sob uma lâmina d’água de 3107 m.
2
posteriormente recolhe-lo. Para uma lâmina d’água de 6500 pés, ou aproximadamente 2000 m, o
tempo total de operação seria de 100 horas, ou quatro dias.
120
100
80
60
40
20
0
0
400
800
1200
1600
2000
2400
2800
3200
3600
4000
4400
4800
5200
5600
6000
6400
6800
7200
Profundidade [pés]
Figura 1-2 Linha de regressão para tempo de descida e recolhimento de BOP/LMRP em função
da profundidade de lâmina d’água (Holand, 2001)
3
poucos os estudos que tratam dela, e não há norma ou boletim emitido por organismo ou
instituição da indústria de petróleo que façam algum tipo de recomendação à sua execução
segura.
O riser durante esta operação está sujeito a esforços gerados pela ação contínua de
correnteza, ondas e movimentos da sonda de perfuração. Estudos publicados a respeito de risers
suspensos tratam principalmente do comportamento apresentado pelo riser para casos extremos
de carregamentos ambientais, como tempestades, que exigem a desconexão do riser da cabeça do
poço. Steddum (2003) discorre a respeito dos fatores que comprometem a resposta do riser e
afetam sua previsão. Ambrose (2001) apresenta um estudo para determinar o efeito que o modo
empregado para ancorar o riser tem na sua capacidade de sobreviver a eventos extremos.
1.2. Objetivos
Este estudo teve como um de seus objetivos analisar a resposta de um riser desconectado à
cabeça de poço, portanto suspenso, mas conectado a um BOP, durante o deslocamento da sonda
de perfuração para um novo local de perfuração. Para tanto, foram feitas análises do
comportamento do riser frente a uma gama variada de esforços ambientais, de forma a reproduzir
uma quantidade apropriada de cenários de carregamentos que podem ser encontrados durante esta
operação.
4
O segundo objetivo deste trabalho tratou, justamente, da delimitação das condições de
carregamentos ambientais, de modo que a execução do deslocamento do riser suspenso pela
sonda respeite uma série de parâmetros de projetos, garantindo a integridade física do riser. O
conhecimento claro dos limites a que está sujeito o riser é essencial para que uma operação
pouco comum como o transporte do BOP suspenso pelo riser possa se tornar rotineira.
O Capítulo 3 apresenta, por sua vez, a metodologia empregada neste estudo. São abordados
os procedimentos adotados para definição da geometria do riser e dos carregamentos ambientais,
considerando-se a inserção destes dados no simulador numérico. É definida a abordagem adotada
para se analisar os resultados obtidos, e como partir destes resultados chega-se aos mapas
operacionais propostos.
5
6
2. DESENVOLVIMENTO
O riser marítimo de perfuração, ou simplesmente riser, é uma tubulação longa e esbelta que
se estende desde a cabeça de poço até a sonda de perfuração, na superfície do mar. Seu corpo
principal consiste de segmentos, ou juntas de riser. Estas são tubos de aço de grande diâmetro, de
comprimento entre 50 e 75 pés, unidas no deck de perfuração e então descidas à água. Esta
operação procede até que o riser atinja seu comprimento total.
Externamente à junta de riser estão acopladas linhas auxiliares com diferentes funções.
Tem-se a linha de kill, que bombeia fluidos de alta densidade para o interior do poço, como forma
de impedir que os fluidos no seu interior cheguem à superfície, como quando ocorre um influxo
de fluidos da formação rochosa sendo perfurada para o interior do poço, ocorrência denominada
kick. A linha de choke conduz o kick, de maneira segura, até a sonda de perfuração. Pode-se ter
uma linha de circulação de fluido de perfuração, chamada booster, para auxiliar no transporte do
cascalho gerado pela operação da broca e também há uma linha por onde são bombeados fluidos
hidráulicos usados para atuar o BOP.
Outros elementos da configuração do riser são a junta telescópica, a junta flexível e a pup
joint.
7
Pup joints são juntas de riser de comprimento menor que o padrão. São utilizadas na
montagem do riser para este alcançar o comprimento projetado.
A junta flexível é um elemento de tubulação instalado tanto no topo do riser, entre ele e a
sonda, quanto na sua extremidade inferior, entre ele e o BOP. Esta junta permite a inclinação
entre o riser e a sonda, ou BOP, diminuindo a intensidade dos esforços sobre o riser devido ao
momento fletor a que o riser está sujeito.
Risers operando em lâminas d’água de até 200 pés não necessitam de outros meios de
suporte além da conexão com a embarcação (Maclachlam, 1987). Em lâminas d’água maiores, e
consequentemente para riser maiores, o risco de falha do riser devido ao seu peso próprio e ao
efeito das forças ambientais fazem necessário um meio de se sustentar o riser. O sistema de
tracionamento de sondas de perfuração supre essa necessidade mantendo o riser sob uma tração
quase constante. Além de sustentar o riser, a tração que atua sobre o riser o torna menos
suscetível a realizar grandes deslocamentos laterais devido à ação das forças ambientais.
A força de tração é exercida por cabos de aço ligados ao anel de tracionamento. Os cabos
são distribuídos ao longo do anel e uma mesma tração uniforme é aplicada. A tração nos cabos é
exercida por cilindros hidráulicos dispostos diagonalmente opostos com suas câmaras ligadas a
uma bateria de garrafas de ar comprimido de alta pressão e grande volume, de modo que o
deslocamento relativo não provoque alteração no valor da tração. Os cabos são ligados aos
cilindros por sistemas de polias, compatibilizando o curso admissível da junta telescópica com o
curso dos pistões.
8
spider repousa sobre outro eqquipamento, chamado gimbal. O gimbal redist
stribui as cargas que
atuam no topo do riser e tambbém atenua as rotações da embarcação que são passados
p ao riser.
9
Spider
Gimbal
Figura 2-2 – Vista expandida de Aranha e Gimbal Fonte: Aker Solutions, 2010
Além disso, o riser também é usado para levar o equipamento BOP até o leito marinho.
Antes do início da montagem do riser, o BOP é conectado à primeira junta que compõe o riser
por intermédio do Lower Marine Riser Package, LMRP. E na medida em que mais juntas são
conectadas, o BOP continua a sua descida.
A
Figura 2-3- Discretização de riser (A), Elemento de viga 2D (B)
O riser pode ser modelado como um elemento de viga carregado axialmente e
transversalmente, este último devido a carregamentos ambientais, e também sofrendo os efeitos
de pressão devido a fluidos internos e externos (Morooka et al, 2005). Patel e Witz (1991)
apresentam um breve histórico dos métodos e técnicas empregados na análise de risers até
chegar-se a equação diferencial que governa o comportamento de risers, para uma análise
realizada no plano.
y y
− + − − γ +γ −γ = (2-1)
11
onde EI é a rigidez à flexão da viga, P0 é a pressão externa, Pi é a pressão interna, A0 é a área
definida pelo diâmetro externo do riser, Ai é a área definida pelo diâmetro interno do riser e As,
por sua vez, é a área definida pela parede do riser. γ0, γi e γs são os pesos específicos do fluido
externo ao riser, do fluido interno ao riser e do material do riser, respectivamente. T representa a
tração sobre o riser e q o carregamento por unidade de comprimento agindo sobre o riser.
O primeiro termo do lado esquerdo da Equação (2-1) representa a rigidez a flexão do riser
devido às propriedades do material que compõem o riser. O segundo termo do lado esquerdo da
equação representa a rigidez à flexão devido aos efeitos da força axial sobre o riser, de acordo
com o conceito de tensão efetiva sob o riser.
O conceito de tensão efetiva foi elaborado como uma forma de agrupar parâmetros de
efeitos similares. A pressão hidrostática externa sob o riser tem um efeito similar ao de uma força
de tração, enquanto que a pressão interna influencia o comportamento do riser como uma força
de compressão.
O terceiro termo do lado esquerdo da equação representa o efeito devido à massa do riser e
o empuxo resultante.
ds
12
Na Figura 2-4, T representa a tração sobre o riser, V representa a força de cisalhamento
sobre suas paredes, Fx0 e Fxi representam as forças horizontais devido às pressões externa e
interna, respectivamente, Fy0 e Fyi representam as forças verticais devido às pressões externa e
interna, respectivamente, e Fc representa a força de arrasto devido à corrente marítima. D0 e Di
são, respectivamente, os diâmetros externo e interno do riser.
13
2.2.2. Análise Estática
K d = F → d = K F (2-2)
onde [K] é a matriz de rigidez global do sistema, {d} é o vetor dos deslocamentos dos nós do
riser discretizado e {F} representa a força que atua sobre o riser.
Adotando que o riser está sujeito a grandes deslocamentos, quando sujeito à ação de
grandes correntes, por exemplo, o problema em mãos se torna não linear, devido à dependência
da matriz de rigidez em relação aos deslocamentos. A solução da equação matricial do equilíbrio
estático é encontrada a partir da aplicação incremental das cargas que atuam sobre o riser.
Devido à natureza oscilatória dos carregamentos a que o riser está sujeito quando em
operação, é necessária a realização de uma análise dinâmica para se determinar como a estrutura
se comporta. A Equação (2-3) apresenta a equação matricial de equilíbrio dinâmico do riser.
M d + B d + K d = F t (2-3)
14
onde M é a matriz de massa global do sistema, B é a matriz de amortecimento global do
sistema, d é o vetor velocidade dos nós que compõem o riser e d é o vetor aceleração dos nós
do riser.
Risers em operação estão sujeitos à ação direta de ondas incidentes e da corrente marítima.
Ondas também atuam indiretamente sobre o riser. Elas causam os movimentos oscilatórios da
embarcação flutuante, os quais são passados ao topo do riser, através de seu ponto de conexão
com a embarcação. No caso da análise realizada neste estudo, se soma ao efeito do carregamento
devido à corrente o efeito devido à velocidade de avanço da embarcação, podendo o perfil de
15
corrente composto ser mais severo ou mais brando que o perfil de corrente original. A Figura 2-6
apresenta os carregamentos atuantes considerados neste estudo.
Movimentos
da sonda
Corrente
ζ t = ζ cos − + (2-4)
16
onde, ζn, ωn, kn e εn representam, respectivamente, a amplitude de onda, frequência, número de
onda e ângulo de fase aleatório da componente de onda n. Enquanto que t e x representam o valor
do tempo para um instante da simulação e a posição da onda ao longo da sua direção de
propagação, respectivamente.
As ondas de mar irregulares de um dado estado de mar podem ser caracterizadas por um
espectro de onda, S(ω), que é a função de densidade de energia espectral do deslocamento
vertical da superfície do mar. A amplitude de cada componente de onda, na, relaciona-se com
sua respectiva frequência a partir do espectro de onda como mostra a Equação (2-5):
1
= Δω
2
(2-5)
17
S(ω)
Altura
Soma Soma
Tempo
Soma
Na estatística de curto prazo, estes parâmetros são considerados constantes e cada conjunto
deles caracteriza um dado estado de mar. A escolha do espectro de onda e de seus parâmetros
característicos é baseada em medições realizadas na posição geográfica de interesse.
18
H / 1 16π 1
= 4π ∙ ∙ exp − ∙
T ω T ω
(2-6)
=F ∙ ∙ ∙ (2-7)
19
Estes movimentos de natureza oscilatória são comumente assumidos gerados unicamente
pela ação das ondas, o vento por sua vez é considerado um carregamento estático, tendo efeito
sobre a embarcação semelhante ao da corrente marítima.
O espectro de energia dos movimentos da embarcação pode ser escrito de maneira análoga
ao espectro de energia de onda, Sζ(ω). Portanto, o espectro de energia para o heave, z(ω,t), pode
ser definido por:
1
Δω =
(2-8)
2
z
= ζ ω S
(2-9)
Para um espectro qualquer S(ω) podem ser calculadas grandezas chamadas momentos
espectrais. A formulação genérica de momento espectral de ordem n, mn, é apresentada pela
Equação (2-10):
20
Onda RAO Resposta
Sζ(ω) Sz(ω)
za/ζa
Componentes Componentes
ω de onda ω regulares da ω
regulares resposta
ζ(t) z(t)
tempo tempo
onda irregular resposta irregular
̅ /
=2 (2-11)
O período de zero ascendente, , por sua vez, é calculado a partir da relação entre o
momento espectral de ordem zero e o momento espectral de ordem dois, como mostra a Equação
(2-12).
=2 ∙ (2-12)
2.3.3. Correntes
A soma vetorial de cada corrente gerada por um diferente fator resulta na corrente total,
cuja velocidade e direção para profundidades especificadas são representadas por meio de um
perfil de corrente, Figura 2-10. Na maioria dos casos a corrente é turbulenta, mas adota-se um
valor correspondente à média observada. No design de estruturas offshore é usual considerar a
corrente como invariante no tempo, geralmente adotando-se o valor máximo para períodos de 10
ou 100 anos, obtido a partir de medições locais.
F= + | | + − + −
= = | |= + −
(2-13)
22
onde F é força hidrodinâmica por unidade de comprimento atuando sobre o riser, ρ é a densidade
do fluido externo ao riser, D0 é o diâmetro externo do riser, CD é o coeficiente de arrasto, CA é o
coeficiente de massa adicionada, u ( é a velocidade (aceleração) da partícula de água, Uc é a
velocidade de corrente e ( ) é a velocidade (aceleração) do riser.
Outro efeito importante devido à ação das correntes sobre o riser é o fenômeno conhecido
por Vibrações Induzidas por Vórtices, ou VIV.
Fluidos escoando ao redor de corpos esbeltos, como risers, geram na esteira deste
escoamento turbilhões que induzem forças cíclicas e desbalanceadas no riser. De tal forma que o
riser passa a oscilar ou vibrar na direção perpendicular à do fluxo observado.
VIV observadas em risers podem levar ao aumento do tempo não produtivo, visto a
necessidade de inspeções e reparos, além de poder causar a falha do equipamento e eventual
perda do poço. Visto sua natureza cíclica, um dos maiores efeitos deletérios decorrentes das VIV
está no seu efeito cumulativo no dano à fadiga que o riser sofre.
Perfil de Corrente
2000
Altura Acima do Leito
1500
Marinho [m]
1000
500
0
00 0,1
0.1 0,2
0.2 0,3
0.3 0,4
0.4
Velocidade [m/s]
Figura 2-10 Perfil de Corrente Marítima.
23
24
3. METODOLOGIA
O objetivo deste capítulo é apresentar os passos adotados na idealização do riser que foi
alvo deste estudo. Além disso, são definidos os limites operacionais adotados para a resposta do
riser e como eles foram elaborados Por fim, apresenta-se o conceito de mapas operacionais, que
delimitam as condições ambientais, sob as quais a mudança de locação do BOP suspenso pelo
riser pode ser realizada respeitando-se aos critérios de operação estabelecidos.
O riser enquanto conectado à cabeça de poço, seja em operação ou não, ou mesmo quando
desconectado devido a eventos meteorológicos extremos, apresenta uma configuração com todos
os elementos e equipamentos descritos no Capítulo 2.
De qualquer forma, para se definir a configuração do riser suspenso, é preciso conhecer sua
configuração enquanto conectado ao poço. O comprimento do riser durante esta etapa pode ser
determinado a partir de procedimento descrito pela norma API RP 16Q (1993). A Figura 3-1
apresenta as dimensões que devem ser levadas em consideração neste cálculo.
25
A dimensão A corresponde à altura da cabeça do poço, medida a partir do leito marinho. A
dimensão B corresponde à altura do conjunto LMRP e BOP. C é o comprimento do riser, o qual
se deseja calcular. D é o comprimento da junta telescópica, quando esta está numa posição de
meio curso. E é a distância do topo do diverter até a mesa rotativa. F é a distância da mesa
rotativa até o leito marinho, enquanto G é o comprimento total da coluna. O comprimento do
riser pode ser determinado tanto pela subtração de F pelas dimensões A, B, D e E, como da
subtração de G pelas dimensões B, D e E.
Mesa Rotativa
E
Junta telescópica D
a meio curso
Base da junta
Altura média da lâmina Comprimento
telescópica
F d’água mais distância do G total da coluna
nível do mar até mesa
rotativa
Comprimento do C
riser
Topo da junta
flexível
Altura do BOP B
LMRP
BOP
Cabeça do poço
Mudline Altura da cabeça A
de poço
Figura 3-1 Determinação do comprimento do riser. Adaptado de API RP 16Q (1993).
26
Juntas de riser revestidas por módulos de flutuabilidade são utilizadas na montagem do
riser de modo a diminuir a capacidade exigida do sistema de tracionamento. De acordo com a
API RP 16Q (1993), a determinação da tração mínima no topo do riser, Tmin, segue o
procedimento de cálculo descrito pelas equações (3-1) e (3-2):
=T ∙
∙ −
(3-1)
T =W ∙ −B ∙ + d ∙ −d ∙ (3-2)
e Ws é o peso submerso do riser, fwt é um fator de tolerância aplicado ao peso submerso do riser,
usualmente igual a 1,05, Bn corresponde à forca de empuxo devido aos módulos de
flutuabilidade, fbt é um fator de perda de flutuabilidade devido ao uso contínuo dos flutuadores,
usualmente igual a 0,96, Ai é a área de seção interna do riser e de todas as linhas auxiliares, dm é
peso específico do fluido de perfuração, Hm corresponde à altura da coluna de fluido de
perfuração, dw corresponde ao peso específico da água, Hw é a altura da lâmina d’água, N é o
número de tracionadores atuando, n corresponde ao número de tracionadores que podem falhar
de maneira súbita e Rf corresponde a um fator de eficiência que relaciona a tração efetiva no anel
de tracionamento com a tração exercida pelos tracionadores, valor usualmente adotado entre 0,90
e 0,95.
Esta metodologia também pode ser empregada para calcular-se o comprimento do trecho de
riser revestido por flutuadores para um valor estabelecido de tração de topo inferior à capacidade
máxima de tracionamento do sistema da embarcação.
Durante a mudança de locação da embarcação se considera que o riser está suspenso pela
spider, a qual está apoiada sobre gimbal, no deck de perfuração da sonda de perfuração, o que
caracteriza a condição de hard hangoff. Os movimentos translacionais da embarcação são
passados integralmente ao topo do riser, mas as rotações são mais ou menos amortecidas a
depender das propriedades dos componentes do gimbal. Este efeito do gimbal sobre o topo do
27
riser pode ser modelado adotando-se que o gimbal é um elemento escalar com valores definidos
de rigidez rotacional nas três direções de rotação.
Para o caso de análise de riser realizada neste estudo não existem critérios de projetos
consolidados na indústria. De tal forma que foram adotados critérios que se julgaram pertinentes
para julgar a viabilidade da operação de transferência do BOP submerso.
Este estudo focou na análise de extremos do riser. A partir do registro da resposta do riser
para cada cálculo realizado, verificou-se se os valores de certas variáveis ficaram abaixo de um
valor estabelecido. Dano à fadiga que o riser incorre durante a mudança de locação do BOP não
foi adotado como critério de projeto.
Limites operacionais de parâmetros chaves para o projeto de risers, como inclinação das
juntas flexíveis inferior e superior, tensões médias e extremas, e fator de segurança apropriado
28
para o limite de tensão dinâmica (DTL) são selecionados com base em princípios de engenharia
sólidos e práticas operacionais bem sucedidas.
A Tabela 3-1 apresenta diretrizes de projeto para três modos operacionais, riser conectado
enquanto ocorre a perfuração do poço, riser conectado sem que se esteja perfurando o poço e
riser desconectado da cabeça de poço.
O terceiro modo é o que mais se assemelha à situação proposta pelo estudo, para ele foi
definido um critério de inclinação de flexjoint e um critério de tensão admissível. O critério de
tensão admissível será adotado como critério de projeto 1 da metodologia aplicada a os estudos
de casos desta dissertação. Os dois métodos propostos para se calcular a tensão admissível,
método A aplicável a água rasas e método B aplicável a águas profundas, são iguais. A tensão
admissível máxima sobre o riser é igual a 67% da tensão de escoamento, σy,do material do riser.
29
Por tensão admissível a norma API RP 16Q (1993) se refere à tensão equivalente de von
Mises ao longo do riser. A tensão equivalente de von Mises é a composição das diversas tensões
que atuam em diferentes direções num dado ponto do riser. Mais formalmente, ela é definida
como a tensão de tração uniaxial que cria uma energia de distorção equivalente à energia de
distorção gerada pelas tensões aplicadas (Hibbeler, 2004).
Os valores de tensão de von Mises, σvm, podem ser calculados como mostra a Equação
(3-3), conforme a norma API RP 16Q (1993), onde σ1, σ2 e σ3 correspondem às tensões
principais atuantes:
1
σ = σ −σ + σ −σ + σ −σ
2
(3-3)
De acordo com a norma API RP 2RD (1998), quando o cisalhamento e a torção são
negligenciáveis, σ1, σ2 e σ3 passam a ser, respectivamente, σpr, σpz e σpθ, onde r, z e θ referem a
tensão radial, axial e circunferencial. A Figura 3-2 apresenta estas tensões atuando sobre um
elemento da parede do riser, onde ri é o raio interno do riser, ro é o raio externo, r é o raio de um
ponto da parede do riser e θ é o ângulo entre o ponto da parede do riser e o eixo Y do referencial
adotado.
r
ri
σpz
ro
Z σpr
θ Y
X σpθ
30
Para um tubo de parede espessa, a tensão circunferencial e a tensão radial podem ser
calculadas como mostrado pelas Equações (3-4) e (3-5), respectivamente, onde C1 e C2 são
definidos pelas Equações (3-6) e (3-7), sendo Pi e Po pressão interna e pressão externa,
respectivamente:
C
σ =C +
r
(3-4)
C
σ =C −
r
(3-5)
Pr − P r
C =
r −r
(3-6)
P −P r r
C =
r −r
(3-7)
A tensão axial tem dois componentes, um componente, σM, devido aos momentos fletores
sobre o riser, Mx e My,e outro componente, σFx, devido à força axial real, Fx, e a força axial
induzida pelas pressões externa e interna. σFx é calculada como mostra a Equação (3-8), onde Ai e
Ao são respectivamente área da seção com raio ri e área da seção com raio ro:
pA −p A F
σ = +
A −A A −A
(3-8)
A contribuição devido ao momento fletor, σM, varia ao longo da seção transversal do riser.
Para um ponto específico da seção, σM é calculado decompondo-se os valores de Mx e My,
respectivamente, momento fletor na direção do eixo X e momento fletor na direção do eixo Y, de
acordo com o referencial da Figura 3-2. Para o ângulo θ da Figura 3-2 a Equação(3-9) indica
como se calcula σM, sendo que Ixy é o momento de inércia do tubo ao redor dos eixos X e Y:
M senθ + M cosθ
σ =
I
(3-9)
31
3.2.2. Critério de Projeto 2
A junta flexível inferior continua instalada, no entanto para simplificar a análise realizada
por este estudo sua inclinação foi desconsiderada em prol de outro critério de inclinação máxima,
pertinente ao gimbal e à geometria da sonda.
Enquanto o riser está apoiado no conjunto spider-gimbal ele vai se inclinar, devido à ação
de correnteza, em relação a este ponto de conexão, que se encontra no nível do deck de
perfuração. O ângulo de topo limite para se evitar colisões depende do espaço existente entre o
diâmetro externo da junta de riser e o diâmetro interno do diverter housing, como mostra o
esquema da Figura 3-3. Tipicamente o ângulo que limita o contato entre riser e diverter housing é
próximo de 3 graus (Chakrabarti, 2005), contato com outras estruturas abaixo do deck de
perfuração ocorreriam para inclinações maiores.
Quando comparado com a configuração de quando o riser está suspenso pelo gancho, o
limite de inclinação para o riser ancorado na spider é maior porque o riser gira ao redor de um
ponto que está bem mais abaixo; tipicamente o diverter housing está localizado quinze pés abaixo
do deck de perfuração.
Assim, o segundo critério de projeto adotado neste estudo, referido por critério de projeto
dois, determina que para se evitar o contato entre o riser e a parede do diverter housing, o ângulo
no topo do riser não pode ultrapassar o valor de três graus.
32
Inclinação do
topo do riser
Ancoragem do
riser na spider
Figura 3-3 – Esquema do contato do riser com estruturas abaixo do deck de perfuração
O riser durante a mudança de locação da sonda está sujeito à ação de corrente marítima,
ondas e, no ponto de conexão com a sonda, aos movimentos da sonda induzidos por ondas. O
perfil de corrente pode ser considerado constante durante a mudança de locação, enquanto que as
ondas se alteram continuamente ao longo do tempo, assim como os movimentos da plataforma.
A resposta do riser devido a estes dois últimos carregamentos oscila, portanto, ao longo do
tempo, ao redor de um valor médio determinado pelos efeitos do perfil de corrente, trazendo o
riser para mais próximo ou para mais distante dos limites de operação estabelecidos. Esta
resposta calculada longo do tempo pode ser registrada para pontos de interesse da estrutura,
criando-se assim gráficos de históricos como o mostrado na Figura 3-4. Nesta figura é possível
observar a variação de tensão de von Mises no topo do riser ao longo de um certo período de
tempo, para uma dada condição de carregamento.
350
Histórico de Tensão de von Mises
Tensão de von
Mises [MPa]
330
310
290
270
250
7000 7050 7100 7150 7200 7250 7300 7350 7400
Tempo de Simulação [s]
Figura 3-4 – Exemplo de histórico de tensão de von Mises.
33
Do histórico extrai-se o valor máximo da resposta do riser, neste caso tensão de von Mises,
e determina-se se o respectivo critério, de tensão admissível no riser, foi respeitado ou não
durante o período de tempo que a resposta do riser foi analisada. O resultado máximo obtido
pode ser extrapolado e ser atribuído não apenas ao período de tempo simulado, mas de maneira
geral ao estado de mar que caracteriza as condições de carregamento impostas durante a
simulação numérica.
Os estados de mar utilizados nas simulações são descritos por um espectro de Pierson-
Moskowitz modificado, que é um espectro de dois parâmetros: altura significativa de onda, H1/3,
e período de zero ascendente, TZ. Cada par H1/3 e TZ tem um valor de tensão de von Mises e
inclinação máxima no topo atribuído a ele, de forma que pode-se apresentar os resultados dos
cálculos realizados compilados na forma de um gráfico de superfície tridimensional, como mostra
a Figura 3-5.
Nos eixos do plano horizontal estão os parâmetros de onda: altura significativa e período de
zero ascendente. No eixo vertical estão os valores de tensão de von Mises máxima no riser. Mas,
a partir do gráfico de superfície é possível criar gráficos de linhas de contorno, como mostrado na
Figura 3-6. Este gráfico apresenta a vantagem de ser mais intuitivo para se trabalhar, visto a
maior facilidade em se determinar o valor atribuído para um par H1/3 e TZ.
34
Superfície de Tensão de von Mises Máxima no Riser
Figura 3-5 - Tensão de von Mises máxima no riser em função de parâmetros de ondas
Figura 3-6 – Curvas de nível de tensão de von Mises máxima no riser em função de parâmetros
de onda
35
Assim, numa situação real de mudança de locação da sonda de perfuração com o BOP
suspenso pelo riser, para as mesmas condições de configuração do riser e de navegação da
embarcação para as quais se criou o gráfico de linhas de contorno da Figura 3-6 se esperaria que
para um estado de mar descrito por um espectro de Pierson-Moskowitz com altura significativa
de onda de 3 metros e período de zero ascendente de 12 segundos, por exemplo, a tensão máxima
de von Mises no riser fosse de 360 MPa, enquanto que para H1/3 = 5 m e TZ = 14 s a tensão de
von Mises seria 420 MPa.
A Figura 3-7 apresenta o mapa operacional criado a partir das curvas de nível de tensão de
von Mises da Figura 3-6 e do critério de projeto 1, da Tabela 3-2, adotando-se σy = 550 MPa.
Analisando-se o mapa operacional, observa-se que enquanto o estado de mar durante a mudança
de locação do BOP suspenso puder ser localizado no interior da região azul do mapa, ou de
maneira mais geral enquanto a altura significativa de onda for inferior a 3,0 metros, o critério de
projeto é respeitado e a integridade do riser é preservada.
A criação de gráficos de curvas de nível pode ser realizada para outros resultados obtidos
pelos cálculos computacionais. A Figura 3-9 apresenta as curvas de nível de dupla amplitude
significativa de heave, HZ1/3, da embarcação. Observa-se como o valor de HZ1/3 tende a aumentar
na medida em que o estado de mar se desloca para o canto superior direito do gráfico, onde tanto
altura significativa de onda e o período característico são maiores.
36
Mapa Operacional – Parâmetros de Onda
Figura 3-7 – Mapa operacional gerado a partir do critério de projeto de tensão em função de
parâmetros de onda
Dupla Amplitude Significativa de Heave [m]
Figura 3-8 – Curvas de nível de dupla amplitude significativa de heave da embarcação em função
de parâmetros de onda
37
Os valores de tensão de von Mises e dupla amplitude significativa de heave referentes a
um mesmo par H1/3 e TZ podem ser relacionados ao valor correspondente calculado de período de
zero ascendente de heave. Dessa forma é possível apresentar as curvas de nível tensão de von
Mises em função dos parâmetros de heave: dupla amplitude e período de zero ascendente, como
mostra a Figura 3-9.
A partir da Figura 3-9 e da adoção do valor para tensão admissível no riser, 366 MPa,
mesmo valor adotado na criação do mapa operacional da Figura 3-7, é possível delimitar a
resposta de heave que não representa riscos à integridade do riser durante a mudança de locação
do BOP.
Observa-se que a zona segura de operação se estende do ponto de dupla amplitude igual a
0,5 metros e período 12 segundos e alcançando o ponto de dupla amplitude de 1.5 metros e
período 18,5 segundos. Isto torna evidente a necessidade de monitoramento não só da amplitude
da resposta como de seu período também.
38
Curvas de Nível de Tensão de von Mises [MPa]
Figura 3-9 – Curvas de nível de tensão de von Mises máxima no riser em função de parâmetros
de heave
Mapa Operacional – Parâmetros de Heave
Figura 3-10 – Mapa operacional gerado a partir do critério de projeto de tensão em função de
parâmetros de heave
39
40
4. ESTUDOS DE CASO
Este capítulo visa apresentar os estudos de casos gerados para se empregar a metodologia
descrita para a obtenção dos mapas operacionais da mudança de locação da sonda de perfuração
com o BOP suspenso pelo riser.
Neste estudo adotou-se que a operação de perfuração foi inicialmente realizada por uma
plataforma semissubmersível sob uma lâmina d’água de 2150 metros. O comprimento do riser
utilizado nesta operação foi determinado a partir de procedimento descrito pela norma API RP
16Q e apresentado no capítulo anterior.
Para o riser operando na lâmina d’água de 2150 m e adotando-se que a mesa rotativa esteja
26 m acima do nível do mar, a dimensão F da Figura 3-1 é igual a 2176 m. Considerando que a
altura da cabeça do poço e altura do conjunto LMRP mais BOP sejam, respectivamente, 1,5 m e
12 m, que o comprimento da junta telescópica a meio curso seja de 26,3 m e também que a
distância da base da junta flexível superior até a mesa rotativa seja de 4,4 m, então o
comprimento do riser é igual a 2131,8 m.
A Tabela 4-1 resume os valores das propriedades da junta de riser utilizada para os cálculos
deste estudo.
41
A Tabela 4-2 abaixo apresenta os valores assumidos para o cálculo do comprimento do
trecho do riser revestido por flutuadores.
A força devido ao empuxo dos flutuadores, Bn, para esta configuração de riser e sistema de
tracionamento é de 861956 lbf, ou 3834,17 kN. O comprimento de flutuadores, com as
características apresentadas na Tabela 4-1, que corresponde a este valor de força é de 2927,2 pés,
ou 892,8 m. A Figura 4-1 apresenta a configuração do riser enquanto em operação de perfuração
e também durante a navegação da embarcação. Na Figura 4-1b está representado o sistema de
referência adotado para este estudo, com origem no leito marinho e o eixo vertical ascendente.
Para a navegação assumiu-se que o riser foi recolhido de um trecho de 150 m.
Juntas de deriser
Junta riser ee
Juntas de deriser
Junta riser ee flutuadores
flutuadores
2132 m 1982 m
flutuadores
flutuadores
2150 m 2150 m Juntas
Junta dede riser
riser
Juntas
Junta de riser
de riser
Junta flexível
42
4.1. Carregamentos Ambientais
Um único perfil de corrente, descrito na Tabela 4-3, foi utilizado ao longo do estudo. Isto se
deve porque o estudo focou em delimitar as condições de navegação da embarcação com o riser e
BOP suspensos em função da resposta da embarcação, que por sua vez é função das
características da onda incidente. Os valores apresentados na Tabela 4-3 são valores médios
calculados a partir de medições realizadas em campo.
As ondas utilizadas neste estudo são ondas irregulares criadas a partir de um espectro de
Pierson-Moskowitz. Os parâmetros que definem um espectro de Pierson-Moskowitz, a altura
significativa e o período de zero ascendente, foram sistematicamente variados, dentro das faixas 1
e 10 m e 8 e 16 s, respectivamente, de forma que foram criados noventa casos de carregamento,
43
abrangendo assim um alto número de estados de mar que a embarcação poderia encontrar
enquanto navegando.
Os movimentos da embarcação devido à ação das ondas são calculados a partir de uma
função de transferência para cada um dos seis graus de liberdade. A Figura 4-2 e a Figura 4-3
apresentam a função RAO para os movimentos de heave, surge, sway, roll, pitch e yaw da
plataforma semissubmersível adotada para este estudo.
Cada tipo de embarcação tem a sua resposta característica representada por suas respectivas
funções RAO. As amplitudes de deslocamento e de rotação são dadas em função do período da
onda incidente de altura unitária. Diferentes ângulos de incidência de onda resultam em respostas
distintas da embarcação, isso pode ser observado para os ângulos de incidência de 0, 45º, 90º,
135° e 180°.
Como este estudo foca na análise de extremos, a exclusão do efeito das VIV, a princípio,
não compromete significativamente os resultados obtidos, visto sua maior influência no estudo da
fadiga a que o riser estaria sujeito, o que está fora do escopo deste estudo. O Apêndice B
apresenta, a título de ilustração, alguns cálculos que contemplam as VIV, além de mostrar como
este e outros carregamentos influencia a resposta global do riser.
44
A partir da interface gráfica do simulador, o usuário pode criar um modelo que será
utilizado na análise da estrutura de interesse. Esta tarefa consiste da geração da geometria do
problema e a especificação dos carregamentos a que a estrutura está sujeita.
Para modelar-se o riser foi necessário criar diversos segmentos distintos, correspondentes a
cada um dos diferentes equipamentos que compõem o riser concebido para este estudo. Assim
têm-se os segmentos correspondentes ao BOP, à junta flexível entre o BOP e o riser, à stress
joint entre o BOP e o riser, às juntas de riser sem módulo de flutuabilidade e às juntas de riser
com módulo de flutuabilidade, cada uma com suas características próprias.
A transição entre segmentos com elementos com comprimentos distintos é feita por
segmentos com elementos de comprimentos variáveis. Definidos os comprimentos do elemento
inicial e do elemento final, o simulador calcula por meio de um algoritmo o número de elementos
45
finitos necessários para compor este elemento, além de também recalcular os comprimentos dos
elementos inicial e final, de acordo com seu algoritmo.
O conjunto BOP e LMRP foi modelado como uma junta de riser com valores muito
maiores de massa linear e rigidez. Foi adotado para o conjunto BOP-LMRP um valor de massa
total de 430 toneladas no ar e um valor de massa aparente na água de 373 toneladas. Estes valores
foram retirados de Ambrose et al (2001).
De acordo com a conexão estabelecida com a sonda, este ponto material está engastado à
sonda, sujeito a seus movimentos induzidos pela ação de ondas, a partir da função de
transferência adotada. As três translações da plataforma são passadas integralmente a este ponto
do riser, enquanto que as três rotações são atenuadas de acordo com as características do gimbal
empregado.
2 Junta Flexível 0 1 0
46
1 200
150
Amplitude de Surge Xa /ζa [-]
0.8
150
Fase de Sway [ graus]
0.8
100
0.6 50
0.4 0
-50
0.2
-100
0 -150
0 10 20 30 0 10 20 30
Período [s] Período [s]
2.5 200
150
Amplitude de Heave Za /ζa [-]
2
Fse de Heave [graus]
100
1.5 50
0
1 -50
-100
0.5
-150
0 -200
0 10 20 30 0 10 20 30
Período [s] Período [s]
0.14 0
-50
0.07
-100
0 -150
0 10 20 30 0 10 20 30
Período [s] Período [s]
0.35 200
Amplitude de Pitch θa /kζa [-]
150
0.28
Fase de Pitch [°] 100
0.21 50
0.14 0
-50
0.07
-100
0 -150
0 10 20 30 0 10 20 30
Período [s] Período [s]
0.0035 200
150
Amplitude de Yaw ψa /kζa [-]
0.0028
100
Fase de Yaw [°]
0.0021 50
0.0014 0
-50
0.0007
-100
0 -150
0 10 20 30 0 10 20 30
Período [s] Período [s]
A Tabela 4-5 apresenta os perfis de corrente alterados para a embarcação navegando com
uma direção de avanço que faz um ângulo θ = 180 grau em relação à direção da corrente
marítima na superfície do mar. Os perfis de corrente alterada para as outras direções de avanço da
embarcação são apresentados no anexo A.
49
4.3. Cenários de Simulação
Foi adotado neste estudo que o riser está em seu topo ancorado na spider, que por sua vez
repousa sobre o gimbal A rigidez rotacional, kθ, do conjunto gimbal-spider determina o quanto
das rotações da embarcação é transferido à extremidade superior do riser. Este valor de rigidez
rotacional é desconhecido, portanto será uma das variáveis de interesse desse estudo.
50
A direção de avanço, θ, que a embarcação pode tomar durante a navegação é medida em
relação à direção de onda e corrente. Foram consideradas neste estudo cinco direções de avanço,
a 0, 45°, 90°, 135° e 180°, correspondentes às direções de a favor do perfil de corrente, à direção
de contracorrente, e três direções intermediárias.
As direções 225°, 270° e 315° puderam ser omitidas devido à simetria. Para o módulo da
velocidade de avanço considerou-se inicialmente três valores: 0,5 nó, 1 nó, e 1,5 nós,
respectivamente 0,257 m/s, 0,514 m/s e 0,781 m/s.
Direção de avanço
X
180°
Z
Y X
Direção de avanço
135°
Z
Onda
X
Direção de avanço
Y
90°
Z
Corrente
X
Direção de avanço
45°
Z
Direção de avanço
X Z
0°
Y
Figura 4-4 – Carregamentos ambientais e direções de avanço da embarcação
51
4.3.3. Quanto aos Carregamentos Ambientais
O foco deste estudo foi a delimitação dos estados de mar sob os quais o riser de perfuração
não falhe durante o transporte do BOP, portanto também variou-se de uma simulação numérica
para outra o espectro de onda de Pierson-Moskowitz adotado, de forma a se cobrir vários
cenários distintos.
O período de zero ascendente dos estados de mar varia dentro da faixa de 8 a 16 segundos,
enquanto que sua altura significativa varia entre 1 e 10 metros.
52
5. RESULTADOS
T [s]
Zh-MIN
Oscilação em Z [m]
Zh-MIN Zh-MAX
Mapas operacionais como apresentados neste estudo representam um conceito novo dentro
da linha de pesquisa de risers. É possível encontrar na literatura científica trabalhos que tratam da
53
analise acoplada da resposta do conjunto embarcação, riser e equipamento submarino, como
apresentados por Nishimoto et al (1989) e Morooka et al (1990). Outros trabalhos publicados
analisam a resposta do riser suspenso. Steddum (2003) discute a influência de vários fatores na
resposta do riser, mas não apresenta resultados.
Portanto, a Figura 5-2 apresenta a comparação do envelope de tração efetiva obtido por
Ambrose et al (2001) com o envelope de tração efetiva obtido pelo simulador numérico
empregado neste estudo.
2500
3750
5000
6250
7500
Ambrose et al (2001)
8750
Sevillano
10000
-1000 -500 0 500 1000 1500 2000 2500 3000
Tração Efetiva [kips]
Figura 5-2 – Comparação de envelopes de tração efetiva
54
É possível observar a grande semelhança entre as curvas referentes aos dois cálculos
computacionais realizados e como elas se sobrepõem para trechos extensos do riser. Levando em
consideração as diferenças inerentes aos programas utilizados, a convergência de resultados
observada pode ser tomada como indicadora da precisão dos resultados obtidos ao longo do
estudo.
A análise a partir dos resultados das simulações computacionais se baseia nos valores
extremos das grandezas medidas. Dentro dos diferentes cenários de velocidade de navegação da
sonda, cada simulação numérica pode ser identificada a partir das características da onda
irregular a que o riser está sujeito.
Assim, o maior valor calculado para deslocamentos, rotações, forças, momentos ou tensões
ao longo do riser, não importando a posição do ponto onde ele ocorre, é atribuído a um par altura
significativa de onda e período de zero ascendente de onda. Independente do tempo a que a
embarcação, num caso real, esteja sujeita a um dado estado de mar, é assumido que a resposta do
riser seja igual à resposta extrema registrada na simulação.
Devido a restrições de espaço não é possível apresentar os resultados obtidos para todas as
condições de carregamento previstas. A Figura 5-3, a Figura 5-4, e a Figura 5-5apresentam os
resultados referentes à simulação numérica na qual o gimbal tem valor de rigidez rotacional de
676 kNm/grau, a embarcação se desloca a 0,5 nó com uma direção de avanço de 180 graus, o
estado de mar é caracterizado pelo espectro de Pierson-Moskowitz com altura significativa de 5
metros e período de zero ascendente de 13 segundos.
55
O Apêndice B desta dissertação apresenta um estudo de sensibilidade do comportamento
do riser frente aos carregamentos a que ele está sujeito, quando estes são impostos
individualmente. No Apêndice B é possível observar de maneira detalhada como cada
carregamento atuante sobre o riser afeta sua reposta, e também apontar, a partir dos critérios
adotados, quais são os carregamentos mais críticos para a operação de trânsito do BOP suspenso.
2000 2000
1800 1800
1600 1600
1400 1400
1200 1200
1000 1000
800 800
600 600
400 400
200 200
0 0
-5 -2.5 0 2.5 5 7.5 10 12.5 15 -4 -3 -2 -1 0 1 2 3
Deslocamento Horizontal [m] Deslocamento Vertical [m]
A B
Valor Máximo Calculado pela Análise Dinâmica
Legenda Valor Mínimo Calculado pela Análise Dinâmica
Valor Calculado pela Análise Estática
Figura 5-3 – Curvas de envoltórias de deslocamento horizontal (A) e deslocamento vertical (B).
kθ = 686 kNm/grau, v = 0,5 nó, θ = 180°, H1/3= 5 m e Tz = 13 s
56
A Figura 5-3 (B) apresenta as curvas de envoltórias de deslocamento vertical. A curva
referente à análise estática indica que o riser se alonga aproximadamente 1,5 m devido ao seu
peso próprio. Para a condição de carregamento particular apresentada, a amplitude máxima de
deslocamento vertical ao longo do riser foi de aproximadamente dois metros. Resultados
similares foram observados para outras condições de carregamento de onda. Isto é, para o sistema
considerado, o BOP, assim como o riser ao longo de seu comprimento, acompanha o movimento
vertical de heave da plataforma flutuante de perfuração.
A Figura 5-4 (A) apresenta as curvas de envoltórias de rotação no plano XZ. O valor
máximo de rotação ao longo do riser foi de aproximadamente 2,1 graus, para um ponto do riser
20 metros abaixo de seu topo. A inclinação no topo do riser não ultrapassou 1,4 graus, portanto
respeitando o critério de inclinação máxima de 3,0 graus.
2000 2000
1800 1800
1600 1600
1400 1400
1200 1200
1000 1000
800 800
600 600
400 400
200 200
0 0
-0.5 0 0.5 1 1.5 2 2.5 3000 4000 5000 6000 7000 8000
Rotação no Plano XZ [graus] Força Axial [kN]
A B
Valor Máximo Calculado pela Análise Dinâmica
Legenda Valor Mínimo Calculado pela Análise Dinâmica
Valor Calculado pela Análise Estática
Figura 5-4 – Curvas de envoltórias de rotação no plano XZ (A) e de força axial efetiva (B). kθ =
686 kNm/grau, v = 0,5 nó, θ = 180°, H1/3= 5 m e Tz = 13 s
57
A Figura 5-4 (B) apresenta as curvas de envoltórias de força axial efetiva. É possível
identificar três trechos distintos nas curvas de força axial efetiva, referentes aos trechos com
diferentes pesos lineares. A partir da extremidade inferior do riser tem-se um trecho composto
por juntas de riser sem flutuador, um trecho composto por juntas de riser com flutuador e,
finalmente, outro trecho composto por juntas de riser sem flutuador. As propriedades das juntas
de riser empregadas no estudo são apresentadas na Tabela 4-1.
A Figura 5-5 (A) apresenta as curvas de envoltórias de momento fletor no plano XZ. O
valor do momento fletor ao longo da maior parte do riser é relativamente baixo, se comparado
com os valores de momento fletor na região próxima ao topo do riser, onde ocorre um salto
significativo no valor de momento fletor.
58
A Figura 5-5Erro! Fonte de referência não encontrada. (B) apresenta as curvas de
envoltória de tensão de von Mises. Também está representada na figura uma linha vertical que
representa o valor de tensão admissível no riser. Observa-se que para esta condição particular de
configuração do riser e carregamentos ambientais, os níveis de tensão ao longo do riser
ultrapassam o valor limite de 366 MPa, definido pelo critério de projeto 1, tornado não
recomendável a navegação sob ação deste estado de mar.
O efeito do salto no valor do momento fletor observado na Figura 5-5 (A) pode ser
observado nas envoltórias de tensão de von Mises, Figura 5-5 (B). O valor de tensão também
aumenta abruptamente próximo à região do topo do riser. Para o restante do riser, a curva de
tensão de von Mises acompanha o formato da curva de força axial efetiva. Resultados obtidos a
partir dos cálculos computacionais realizados indicam que o máximo valor de tensão de von
Mises ao longo do riser ocorre sempre na região próxima ao topo.
Nos eixos das abscissas e das ordenadas deste gráfico se tem os parâmetros de onda período
de zero ascendente e altura significativa de onda, respectivamente. Vale lembrar que estes são os
parâmetros necessários para o cálculo do espectro de onda de Pierson-Moskowitz. As curvas de
nível representam o valor máximo esperado de tensão de von Mises ao longo do riser para os
estados de mar simulados.
59
Curvas de Nível de Tensão de von Mises [MPa]
Figura 5-6 – Curvas de nível de tensão de von Mises máxima no riser; kθ = 686 kNm/grau,
v = 0,5 nó, θ = 180°
É possível observar que o nível de tensão de von Mises varia entre 350 MPA e 650 MPA.
O critério de projeto 1 prevê que a tensão de von Mises não pode ultrapassar dois terços da tensão
de escoamento do material. Este limite está demarcado no mapa operacional pela curva de nível
de valor 366 MPa.
Portanto o limite de operação de acordo com esse critério de projeto é um estado de mar
descrito por um espectro de mar de Pierson-Moskowitz com altura significativa de onda entre 3,0
e 3,5 metros, para toda a faixa de valores de períodos de zero ascendente.
A Figura 5-7 apresenta o mapa operacional gerado para o Critério de Projeto 2, para o
mesmo cenário do mapa operacional da Figura 5-6. A mesma descrição quanto aos eixos do
gráfico se aplica a esta figura, o que a difere da anterior é que agora as curvas de níveis
representam valores de inclinação máxima no topo do riser, na conexão entre este e a
embarcação.
60
O Critério de Projeto 2 prevê uma inclinação para o riser neste ponto, de forma a garantir
que não haja contato entre o corpo do riser e a embarcação. O ângulo no topo do riser varia entre
1,0 graus e 2,5 graus, portanto abaixo do limite imposto de 3,0 graus. Para qualquer estado de
mar representado no mapa operacional, a presente combinação de rigidez rotacional do gimbal e
velocidade de avanço da embarcação respeita o critério de inclinação máxima de 3,0 graus.
Figura 5-7 - Curvas de nível de inclinação máxima no topo do riser; kθ = 686 kNm/grau,
v = 0,5 nó, θ = 180°
Figura 5-8 – Curvas de nível de tensão de von Mises máxima no riser em função dos parâmetros
de heave; kθ = 686 kNm/grau, v = 0,5 nó, θ = 180°
Figura 5-9 – Curvas de nível de rotação no topo do riser em função dos parâmetros de heave;
kθ = 686 kNm/grau, v = 0,5 nó, θ = 180°
O mapa operacional de tensão de von Mises pode também ser apresentado de uma outra
maneira, de forma a tornar explícito no mapa o critério para o valor de tensão de von Mises
adotado para sua criação, ou permitir que um usuário do mapa adote algum outro critério que
achar mais apropriado para o valor de tensão de von Mises. Isso é feito dividindo o valor das
curvas de níveis pelo valor adotado de tensão de escoamento do material, tornando as curvas de
níveis adimensionais.
63
De forma que, segundo o critério adotado neste estudo, a curva de nível referente ao valor
0,67 passa a demarcar a separação entre os estados de mar que respeitam o valor de tensão
admissível daqueles estados de mar que o ultrapassam.
Uma vez que um critério foi estabelecido, apenas a curva referente a este critério basta para
usar o mapa operacional, assim é possível plotar no mesmo mapa operacional as curvas referentes
aos dois critérios de projetos, sem perda de legibilidade. Em favor da segurança, a curva que
delimite uma condição de operação mais restritiva deve ser usada para delimitar os estados de
mar que a embarcação pode enfrentar durante a mudança de locação. Caso estas duas curvas se
cruzem, é possível mesclá-las.
A região em azul demarca no mapa operacional para quais estados de mar a mudança de
locação do BOP suspenso pelo riser pode ser executada respeitando-se ambos os critérios
estabelecidos, para uma velocidade de avanço de 0,5 nó. Portanto, o mapa operacional da Figura
5-13 indica que para esta combinação de rigidez rotacional do gimbal e direção de avanço da
embarcação, o transporte do riser e do BOP suspensos poderia ocorrer a uma velocidade a 0,5 nó,
desde que a altura significativa de onda não ultrapasse os 3,0 metros.
64
Mapa Operacional – Parâmetros de Onda
Figura 5-10 – Mapa operacional para direção de avanço da embarcação de 180 graus, gimbal de
rigidez rotacional de 676 kNm/grau. Parâmetros de onda.
As figuras de número 5-15 a 5-42, a seguir, apresentam os mapas operacionais gerados para
as cinco direções de navegação e os três valores adotado de rigidez rotacional do gimbal, tanto
em função dos parâmetros de onda como em função dos parâmetros de respostra da embarcação.
No Apêndice C se encontram os grafícos de curvas de nível a partir dos quais criaram-se os
mapas operacionais apresentados.
65
Mapa Operacional – Parâmetros de Heave
Figura 5-11 – Mapa operacional para direção de avanço da embarcação de 180 graus, gimbal de
rigidez rotacional de 676 kNm/grau. Parâmetros de heave.
66
Mapa Operacional – Parâmetros de Onda
Figura 5-12 – Mapa operacional para direção de avanço da embarcação de 180 graus, gimbal de
rigidez rotacional de 338 kNm/grau. Parâmetros de onda.
Mapas Operacionais – Parâmetros de Heave
Figura 5-13 – Mapa operacional para direção de avanço da embarcação de 180 graus, gimbal de
rigidez rotacional de 338 kNm/grau. Parâmetros de heave.
67
Mapa Operacional – Parâmetros de Onda
Figura 5-14 – Mapa operacional para direção de avanço da embarcação de 180 graus, gimbal de
rigidez rotacional de 169 kNm/grau. Parâmetros de onda.
Mapas Operacionais – Parâmetros de Heave
Figura 5-15 – Mapa operacional para direção de avanço da embarcação de 180 graus, gimbal de
rigidez rotacional de 169 kNm/grau. Parâmetros de heave.
68
Mapa Operacional – Parâmetros de Onda
Figura 5-16 – Mapa operacional para direção de avanço da embarcação de 135 graus, gimbal de
rigidez rotacional de 676 kNm/grau. Parâmetros de onda.
Mapas Operacionais – Parâmetros de Heave
Figura 5-17 – Mapa operacional para direção de avanço da embarcação de 135 graus, gimbal de
rigidez rotacional de 676 kNm/grau. Parâmetros de heave.
69
Mapaa Operacional – Parâmetros de Onda
70
Mapaa Operacional – Parâmetros de Onda
71
Mapa Operacional – Parâmetros de Onda
Figura 5-22 – Mapa operacional para direção de avanço da embarcação de 90 graus, gimbal de
rigidez rotacional de 676 kNm/grau. Parâmetros de onda.
Mapas Operacionais – Parâmetros de Heave
Figura 5-23 – Mapa operacional para direção de avanço da embarcação de 90 graus, gimbal de
rigidez rotacional de 676 kNm/grau. Parâmetros de heave.
72
Mapaa Operacional – Parâmetros de Onda
73
Mapaa Operacional – Parâmetros de Onda
74
Mapaa Operacional – Parâmetros de Onda
75
Mapaa Operacional – Parâmetros de Onda
76
Mapa Operacional – Parâmetros de Onda
Figura 5-32 – Mapa operacional para direção de avanço da embarcação de 45 graus, gimbal de
rigidez rotacional de 169 kNm/grau. Parâmetros de onda.
Figura 5-33 – Mapa operacional para direção de avanço da embarcação de 45 graus, gimbal de
rigidez rotacional de 169 kNm/grau. Parâmetros de heave.
77
Mapaa Operacional – Parâmetros de Onda
Figura 5-34 – Mapa operaccional para direção de avanço da embarcação de 0 grau, gimbal de
rigidez rootacional de 676 kNm/grau. Parâmetros de onda..
Figura 5-35 – Mapa operaccional para direção de avanço da embarcação de 0 grau, gimbal de
rigidez rottacional de 676 kNm/grau. Parâmetros de heave.
78
Mapaa Operacional – Parâmetros de Onda
Figura 5-36 – Mapa operaccional para direção de avanço da embarcação de 0 grau, gimbal de
rigidez rootacional de 338 kNm/grau. Parâmetros de onda..
Figura 5-37 – Mapa operaccional para direção de avanço da embarcação de 0 grau, gimbal de
rigidez rottacional de 338 kNm/grau. Parâmetros de heave.
79
Mapa Operacional – Parâmetros de Onda
Figura 5-38 – Mapa operacional para direção de avanço da embarcação de 0 grau, gimbal de
rigidez rotacional de 169 kNm/grau. Parâmetros de onda.
Mapas Operacionais – Parâmetros de Heave
Figura 5-39 – Mapa operacional para direção de avanço da embarcação de 0 grau, gimbal de
rigidez rotacional de 169 kNm/grau. Parâmetros de heave.
80
A navegação a 1,5 nós não obedece aos critérios de projetos, para qualquer combinação de
rigidez rotacional do gimbal e direção de avanço, por isso não é apresentada nos mapas
operacionais.
É possível observar que a condição mais restritiva de navegação é para a direção de avanço
de 135 graus e para o riser ancorado em um gimbal com rigidez rotacional de 676 kNm/grau,
retratada na Figura 5-16 e na Figura 5-17, visto que a região branca destes mapas operacionais é a
maior de todas.
Nestas condições, os dois critérios de projeto são respeitados apenas para uma velocidade
de avanço de 0,5 nó. Para esta velocidade, a altura significativa de onda limite varia entre 4,2 e
2,3 metros.
A amplitude significativa de heave pode ter o valor máximo de 0,5 metro para um período
de zero ascendente de aproximadamente 13,3 segundos, enquanto que para um período de 20
segundos a amplitude significativa de heave não pode ultrapassar 1,3 metros.
A condição menos severa de navegação, por sua vez, ocorre para a direção de avanço de 45
graus e para o riser ancorado em um gimbal com rigidez rotacional de 169 kNm/grau, retratada
na Figura 5-35 e na Figura 5-36, visto que estes mapas estão completamente coloridos de verde.
Nestas condições, a navegação à velocidade de 0,5 nó pode ser realizada para qualquer
estado de mar dentre os representados no mapa operacional. Em relação à navegação à
velocidade de 1,0 nó, esta pode ser efetuada para qualquer estado de mar cuja altura significativa
não ultrapasse 8,0 metros, independente do valor de período de zero ascendente.
O mapa operacional criado em função dos parâmetros de heave apresenta uma região verde
que o cobre por completo. Isto indica que qualquer que seja a resposta de heave da embarcação
dentre os limites estabelecidos no mapa, o transporte do BOP submerso pode ocorrer a uma
velocidade de navegação de até 1,0 nó.
O efeito da rigidez rotacional adotada nas simulações numéricas é observado no que diz
respeito a qual critério de projeto acaba por ser restringir a menor área de operação dentro dos
mapas operacionais.
81
Observa-se que nas simulações numéricas em que o valor de rigidez rotacional do gimbal
foi de 686 kNm/grau ou de 338 kNm/grau, a tensão de von Mises máxima no riser foi o fator
determinante na hora de delimitar os estados de mar que respeitam os critérios de projeto.
Já nas simulações numéricas em que foi adotado para o gimbal o valor de rigidez rotacional
de 169 kNm/grau, a inclinação máxima no topo do riser foi, na maioria dos casos, o fator
determinante na hora de delimitar os estados de mar.
A Tabela 5-1 resume os resultados das simulações. Nela são apresentadas para cada uma
das diversas combinações de velocidade e de rigidez rotacional se os critérios adotados são
satisfeitos ou não.
A variação nos valores de tensão de von Mises observada para os diversos cálculos
computacionais executados é explicada pela variação no valor do momento fletor frente a
diferentes carregamentos e condições de contorno no topo do riser.
A força axial efetiva mantém-se praticamente constante ao longo das simulações, como
mostram as Figura 5-40 (A) a (D), que comparam as curvas de nível de força axial efetiva no
riser geradas a partir dos resultados de simulações computacionais executadas para quatro
cenários distintos, em que se alteraram os valores de velocidade de deslocamento da embarcação,
sua direção de avanço e/ou o valor de rigidez rotacional do gimbal.
O momento fletor máximo no riser, por sua vez, apresenta valores significativamente
maiores para a configuração de riser com gimbal de rigidez rotacional de 686 kNm//grau, como
mostra as Figura 5-41 (A) e (B). Para simulações numéricas com condições de contorno de topo
idênticas, tanto a direção de avanço como a velocidade de avanço afetaram o valor do momento
fletor, como pode ser observado nas Figura 5-41 (C) e (D).
82
Estes resultados mostram a importância do valor do momento fletor a que o riser está
submetido. A análise das curvas de envoltórias, como as representadas na Figura 5-3, na Figura
5-4 e na Figura 5-5, mostram que o valor máximo de momento ocorre sempre no topo do riser.
Isso implica não só na importância de se modelar corretamente esta condição de contorno, mas de
que procedimentos especiais podem ser adotados para se mitigar a tensão no topo do riser, como,
por exemplo, o uso de uma stress joint, ao invés de uma junta de riser comum, para se ancorar o
riser.
Em relação à inclinação no topo do riser, ela passa a ser restritiva para as operações de
deslocamento do riser com BOP suspenso para valores de velocidade de avanço superiores a 1,0
nó. Isso porque foi adotado o limite de três graus de inclinação, obviamente o valor máximo de
inclinação a que o riser pode estar sujeito vai depender do arranjo espacial da embarcação que
realiza o transporte, o que pode vir a afrouxar ou enrijecer este critério, tornando possível o
deslocamento a velocidades maiores.
A análise de extremos realizada neste trabalho contempla dois critérios. Outros critérios
podem ser acrescentados em análises futuras, ou até substituir os critérios originais, quanto à
delimitação de estados de mar.
Como este trabalho realizou uma análise de extremos, fadiga do riser não foi um tema
abordado. Justamente devido ao seu maior impacto no dano à fadiga no riser, o efeito das VIV
não foi incorporado aos cálculos realizados, de forma a acelerar e simplificar a implementação do
trabalho computacional requerido para este estudo.
83
Por mais que se tentou criar uma gama variada de cenários de carregamentos ambientais, na
medida em que o perfil de corrente onde a operação ocorre mude, os resultados passam a ser
menos confiáveis, chegando ao ponto de que novos cálculos tenham que ser realizados.
84
Tabela 5-1 – Resumo dos resultados obtidos no estudo
Rigidez Rotacional do Gimbal
169 338 676
[kNm/grau]
85
Curvas de Nível de Força Axial Efetiva Máxima no Riser [kN]
A B
C D
Figura 5-40 – Curvas de nível de força axial efetiva máximo no riser.
(A) v = 0,5 nó, θ = 180°, kθ = 686 kNm/grau. (B) v = 0,5 nó, θ = 180°, kθ = 169 kNm/grau.
(C) v = 0,5 nó, θ = 90°, kθ = 686 kNm/grau. (D) v = 1,0 nó, θ = 180°, kθ = 686 kNm/grau.
86
Curvas de Nível de Momento Fletor Máximo no Riser [kNm]
A B
C D
Figura 5-41 – Curvas de nível de momento fletor máximo no riser.
(A) v = 0,5 nó, θ = 180°, kθ = 686 kNm/grau. (B) v = 0,5 nó, θ = 180°, kθ = 169 kNm/grau.
(C) v = 0,5 nó, θ = 90°, kθ = 686 kNm/grau. (D) v = 1,0 nó, θ = 180°, kθ = 686 kNm/grau.
87
88
6. CONCLUSÕES
Portanto, a tensão de von Mises ao longo do riser não poderia ultrapassar um valor de
tensão admissível, de forma a garantir que o riser não escoe, e a inclinação no topo do riser
não poderia ultrapassar um valor máximo de inclinação, de forma a se respeitar restrições
geométricas da sonda, garantindo que o riser não entre em contato com estruturas abaixo
do deck de perfuração.
• O estudo do dano à fadiga que o riser incorre durante esta operação e qual seu
impacto na vida útil do equipamento e como incluir a fadiga aos critérios de
projeto já estabelecidos.
90
• A análise de práticas, recomendações e estudos a respeito de quanto tempo BOP e
riser podem ser mantidos submersos, antes que seja necessário trazê-los à
superfície para inspecionar suas condições.
• O estudo de como se realizar a aquisição dos dados ambientais dos carregamentos
de ondas e corrente que atuam sobre o riser.
91
92
7. REFERÊNCIAS BIBLIOGRÁFICAS
AKER SOLUTIONS, 2010. Deep water CLIP RiserTM for marine drilling operations.
Disponível em: <http://www.akersolutions.com/Documents/Drilling%20Technologies/
Drilling%20risers/Updated-brochure-with-fold-out-low-res.pdf>. Acesso em 21 mar. 2013.
AMBROSE, B. D., GREALISH, F., WHOOLEY, K. Soft Hangoff Method for Drilling
Risers in Ultra Deepwater. In: Offshore Technology Conference, 2001, Houston, Texas.
BAI, Y. Pipeline and Risers. Oxford: Elsevier Science, 2001. 526p. Vol. 3.
BAI, Y, BAI Q. Subsea Pipelines and Risers. Oxford: Elsevier Science, 2005. 841p.
HIBBELER, R. C. Resistência dos Materiais. 5 ed. São Paulo: Prentice Hall, 2004.
93
HOLAND, P. Reliability of Deepwaters Subsea Blowout Preventers. SPE Drilling and
Completion, vol 16, no. 1, p. 12-18, 2001.
HULL, J. P., BOP-deployed move saves time, money. Offshore, volume 62, número 6, p.
36, 2002.
MAHONEY, C. N., SUPAN, C. 2012 Deepwaer Solutions and Records for Concept
Selection. Offshore, 2012. Disponível em: <www.offshore-mag.com/content/dam/
offshore/print-articles/ Volume%2072/may/0512OS_DeepwaterPoster040512Ads.pdf>.
Acesso em 21 mar. 2013.
94
24th International Conference on Offshore Mechanics and Arctic Engineering (OMAE),
Halkidiki, Grécia, 2005.
MOROOKA, C. K., COELHO, F. M., MATT, C. G. C., FRANCIS, R., A Time Domain
Simulation Tool for Rigid Riser’s Design. In: Rio Oil & Gas Expo and Conference, 2006,
Rio de Janeiro: Instituto Brasileiro de Petróleo e Gás, 2006, IBP.
95
APÊNDICE A – CARREGAMENTOS AMBIENTAIS
97
Tabela A2 – Perfil de corrente alterado para a direção de avanço de 45 graus
Avanço a 0,5 nós Avanço a 1,0 nós Avanço a 1,5 nós
Profundidade [m] Velocidade Azimute Velocidade Azimute Velocidade Azimute
[m/s] [°] [m/s] [°] [m/s] [°]
0 0.47 203 0.72 211 0.97 214
10 0.47 203 0.72 211 0.97 214
25 0.48 194 0.71 205 0.95 210
50 0.59 208 0.84 213 1.09 216
75 0.56 203 0.80 210 1.05 214
100 0.56 210 0.81 215 1.06 217
200 0.46 202 0.70 210 0.95 214
300 0.41 206 0.66 213 0.91 216
400 0.36 216 0.61 220 0.87 222
500 0.34 214 0.60 219 0.85 221
600 0.38 210 0.63 216 0.88 218
700 0.29 233 0.55 229 0.81 228
800 0.33 218 0.58 221 0.84 222
900 0.35 221 0.60 222 0.86 223
1000 0.25 239 0.50 232 0.76 230
1250 0.30 234 0.56 230 0.81 228
1500 0.34 239 0.59 233 0.85 230
1750 0.32 235 0.57 230 0.83 229
2000 0.29 234 0.55 230 0.80 228
2150 0.26 225 0.51 225 0.77 225
98
Tabela A3 – Perfil de corrente alterado para a direção de avanço de 90 graus
Avanço a 0,5 nós Avanço a 1,0 nós Avanço a 1,5 nós
Profundidade [m] Velocidade Azimute Velocidade Azimute Velocidade Azimute
[m/s] [°] [m/s] [°] [m/s] [°]
0 0.36 226 0.57 244 0.62 298
10 0.36 226 0.57 244 0.62 298
25 0.34 215 0.53 238 0.55 299
50 0.49 226 0.70 241 0.67 288
75 0.44 222 0.64 239 0.62 290
100 0.47 230 0.68 244 0.69 291
200 0.35 225 0.56 244 0.61 299
300 0.32 233 0.55 250 0.65 303
400 0.31 250 0.56 259 0.73 307
500 0.29 249 0.54 259 0.71 308
600 0.30 241 0.54 254 0.68 306
700 0.31 271 0.57 271 0.82 313
800 0.29 254 0.54 262 0.73 310
900 0.31 255 0.56 262 0.75 308
1000 0.29 281 0.55 276 0.83 316
1250 0.32 271 0.58 271 0.82 312
1500 0.36 271 0.62 271 0.85 310
1750 0.33 270 0.59 270 0.83 311
2000 0.31 272 0.57 271 0.82 313
2150 0.26 270 0.51 270 0.77 315
99
Tabela A4 – Perfil de corrente alterado para a direção de avanço de 135 graus
Avanço a 0,5 nós Avanço a 1,0 nós Avanço a 1,5 nós
Profundidade [m] Velocidade Azimute Velocidade Azimute Velocidade Azimute
[m/s] [°] [m/s] [°] [m/s] [°]
0 0.19 249 0.38 287 0.87 303
10 0.19 249 0.38 287 0.87 303
25 0.15 230 0.31 286 0.80 304
50 0.32 240 0.46 273 0.91 295
75 0.26 236 0.40 275 0.86 297
100 0.30 247 0.47 278 0.93 297
200 0.18 250 0.37 289 0.86 304
300 0.18 267 0.40 296 0.90 307
400 0.23 290 0.47 303 0.98 309
500 0.21 292 0.46 305 0.97 310
600 0.19 281 0.43 301 0.93 309
700 0.30 308 0.56 311 1.07 313
800 0.23 297 0.48 307 0.99 311
900 0.25 294 0.50 305 1.01 310
1000 0.32 318 0.58 317 1.09 316
1250 0.31 308 0.56 311 1.08 313
1500 0.34 303 0.60 308 1.11 311
1750 0.32 305 0.57 310 1.08 312
2000 0.30 309 0.56 312 1.07 313
2150 0.26 315 0.51 315 1.03 315
100
Tabela A5 – Perfil de corrente alterado para a direção de avanço de 180 graus
Avanço a 0,5 nós Avanço a 1,0 nós Avanço a 1,5 nós
Profundidade [m] Velocidade Azimute Velocidade Azimute Velocidade Azimute
[m/s] [°] [m/s] [°] [m/s] [°]
0 0.01 360 0.26 360 0.52 360
10 0.01 360 0.26 360 0.52 360
25 0.07 108 0.24 15 0.50 7
50 0.13 228 0.20 332 0.44 348
75 0.08 208 0.19 348 0.44 355
100 0.11 246 0.24 335 0.48 348
200 0.02 44 0.27 3 0.53 1
300 0.07 2 0.32 0 0.58 0
400 0.15 349 0.41 356 0.67 357
500 0.15 356 0.41 358 0.67 359
600 0.11 358 0.37 359 0.63 360
700 0.27 348 0.52 354 0.78 356
800 0.18 354 0.44 358 0.69 358
900 0.18 347 0.44 354 0.69 356
1000 0.32 355 0.57 357 0.83 358
1250 0.27 347 0.52 353 0.78 355
1500 0.28 338 0.53 349 0.79 352
1750 0.27 343 0.52 351 0.78 354
2000 0.27 349 0.53 354 0.78 356
2150 0.26 360 0.51 360 0.77 360
101
102
APÊNDICE B – EFEITOS DE CARREGAMENTOS
Pela própria natureza da corrente, que é assumida constante, e também de como são
realizados os cálculos computacionais, o riser assume uma configuração inicial decorrente
da ação da corrente e passa a oscilar ao redor desta configuração devido aos carregamentos
dinâmicos.
O caso 10 difere do anterior pelo acréscimo das Vibrações Induzidas por Vórtices
(VIV). Apesar de ser dito no capítulo 4 desta dissertação que a versão do simulador
numérico utilizado não contemplar as VIV, o grupo de pesquisa da UNICAMP conseguiu
implementar numa versão anterior deste programa, de código aberto, o cálculo das VIV. As
figuras B1 a B28 apresentam as curvas de envoltórias geradas a partir dos resultados
obtidos dos cálculos realizados.
O efeito no valor extremo de tensão de von Mises devido ao momento fletor gerado a
partir das VIV está na ordem de kilo Pascais. Por este motivo, as VIV foram
desconsideradas no estudo apresentado por esta dissertação, que foca nos valores extremos
da resposta do riser.
104
Comparação - Deslocamento Direção Inline
2200
2000 Cor
Cor + OR
1800 Cor + He
Cor + Su
1600 Cor + Sw
Posição no R iser [m ]
Cor + Ro
1400 Cor + Pi
Cor + OR + Mov
1200
1000
800
600
400
200
0
-2 0 2 4 6 8 10 12
Deslocamento Direção Cross-Flow [m]
Figura B 1 – Comparação entre efeitos de carregamentos regulares no envelope de
deslocamento na direção inline
2000
1800
1600
Posição no R iser [m ]
1400
1200
1000
Cor
800 Cor + OR
Cor + He
600 Cor + Su
Cor + Sw
400 Cor + Ro
Cor + Pi
200 Cor + OR + Mov
0
-0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8
Deslocamento Direção Cross-Flow [m]
Figura B 2 – Comparação entre efeitos de carregamentos regulares no envelope de
deslocamento na direção cross-flow
105
Comparação - Deslocamento Direção Axial
2200
2000
1800
1600
Posição no R iser [m ]
1400
1200
1000
Cor
800 Cor + OR
Cor + He
600 Cor + Su
Cor + Sw
400 Cor + Ro
Cor + Pi
200 Cor + OR + Mov
0
-2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2
Deslocamento Direção Cross-Flow [m]
Figura B 3 – Comparação entre efeitos de carregamentos regulares no envelope de
deslocamento de direção axial
2000
1800
1600
Posição no R iser [m ]
1400
1200
1000
Cor
800 Cor + OR
Cor + He
600 Cor + Su
Cor + Sw
400 Cor + Ro
Cor + Pi
200 Cor + OR + Mov
0
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1
Rotação [graus]
Figura B 4 – Comparação entre efeitos de carregamentos regulares no envelope de
rotação no plano definido pelas direções vertical e inline
106
Comparação - Rotação no Plano Definido Pelas Direções Vertical e Cross-Flow
2200
2000
1800
1600
Posição no R iser [m ]
1400
1200
1000
Cor
800 Cor + OR
Cor + He
600 Cor + Su
Cor + Sw
400 Cor + Ro
Cor + Pi
200 Cor + OR + Mov
0
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1
Rotação [graus]
Figura B 5 – Comparação entre efeitos de carregamentos regulares no envelope de
rotação no plano definido pelas direções vertical cross-flow
2000
1800
1600
Posição no R iser [m ]
1400
1200
1000
Cor
800 Cor + OR
Cor + He
600 Cor + Su
Cor + Sw
400 Cor + Ro
Cor + Pi
200 Cor + OR + Mov
0
-250 -200 -150 -100 -50 0 50 100
Força [kN]
Figura B 6 – Comparação entre efeitos de carregamentos regulares no envelope de
força na direção inline
107
Comparação - Força Direção Cross-Flow
2200
2000
1800
1600
Posição no R iser [m ]
1400
1200
1000
Cor
800 Cor + OR
Cor + He
600 Cor + Su
Cor + Sw
400 Cor + Ro
Cor + Pi
200 Cor + OR + Mov
0
-150 -100 -50 0 50 100 150
Força [kN]
Figura B 7 – Comparação entre efeitos de carregamentos regulares no envelope de
força na direção cross-flow
2000 Cor
Cor + OR
1800 Cor + He
Cor + Su
1600 Cor + Sw
Posição no R iser [m ]
Cor + Ro
1400 Cor + Pi
Cor + OR + Mov
1200
1000
800
600
400
200
0
0 1000 2000 3000 4000 5000 6000 7000 8000
Força [kN]
Figura B 8 – Comparação entre efeitos de carregamentos regulares no envelope de
força axial
108
Comparação - Momento Fletor no Plano Definido pelas Direções Vertical e Inline
2200
2000
1800
1600
Posição no R iser [m ]
1400
1200
1000
Cor
800 Cor + OR
Cor + He
600 Cor + Su
Cor + Sw
400 Cor + Ro
Cor + Pi
200 Cor + OR + Mov
0
-1400 -1200 -1000 -800 -600 -400 -200 0 200 400
Momento [kN*m]
Figura B 9 – Comparação entre efeitos de carregamentos regulares no envelope de
momento fletor no plano definido pelas direções vertical e inline
2000
1800
1600
Posição no R iser [m ]
1400
1200
1000
Cor
800 Cor + OR
Cor + He
600 Cor + Su
Cor + Sw
400 Cor + Ro
Cor + Pi
200 Cor + OR + Mov
0
-800 -600 -400 -200 0 200 400 600 800
Momento [kN*m]
Figura B 10 – Comparação entre efeitos de carregamentos regulares no envelope de
momento fletor no plano definido pelas direções vertical e cross-flow
109
Comparação - Tensão devido ao Momento Fletor
2200
2000
1800
1600
Posição no R iser [m ]
1400
1200
1000
Cor
800 Cor + OR
Cor + He
600 Cor + Su
Cor + Sw
400 Cor + Ro
Cor + Pi
200 Cor + OR + Mov
0
0 50 100 150 200 250 300 350
Tensão [MPa]
Figura B 11 – Comparação entre efeitos de carregamentos regulares na tensão devido
ao momento fletor
2000 Cor
Cor + OR
1800 Cor + He
Cor + Su
1600 Cor + Sw
Posição no R iser [m ]
Cor + Ro
1400 Cor + Pi
Cor + OR + Mov
1200
1000
800
600
400
200
0
-50 0 50 100 150 200 250
Tensão [MPa]
Figura B 12 – Comparação entre efeitos de carregamentos regulares na tensão devido à
força axial
110
Comparação - Tensão de Cisalhamento
2200
2000
1800
1600
Posição no R iser [m ]
1400
1200
1000
Cor
800 Cor + OR
Cor + He
600 Cor + Su
Cor + Sw
400 Cor + Ro
Cor + Pi
200 Cor + OR + Mov
0
-20 -15 -10 -5 0 5
Tensão [MPa]
Figura B 13 – Comparação entre efeitos de carregamentos regulares no envelope de
tensão cisalhante
2000
1800
1600
Posição no R iser [m ]
1400
1200
1000
Cor
800 Cor + OR
Cor + He
600 Cor + Su
Cor + Sw
400 Cor + Ro
Cor + Pi
200 Cor + OR + Mov
0
0 100 200 300 400 500 600
Tensão [MPa]
Figura B 14 – Comparação entre efeitos de carregamentos regulares no envelope de tensão
de von Mises
111
Comparação - Deslocamento Direção Inline
2000
1800
1600
P o siç ão no R iser [m ]
1400
1200
1000
800
600
H1/3=3 m Tz=8 s sem VIV
H1/3=3 m Tz=8 s com VIV
400 H1/3=5 m Tz=13 s sem VIV
H1/3=5 m Tz=13 s com VIV
200 H1/3=3 m Tz=16 s sem VIV
H1/3=3 m Tz=16 s com VIV
0
-2 0 2 4 6 8 10 12 14
Deslocamento Direção Cross-Flow [m]
Figura B 15 – Comparação do efeito do acréscimo de VIV no envelope de deslocamento na
direção inline
1800
1600
Po siçã o n o R ise r [m ]
1400
1200
1000
800
H1/3=3 m Tz=8 s sem VIV
H1/3=3 m Tz=8 s com VIV
600
H1/3=5 m Tz=13 s sem VIV
H1/3=5 m Tz=13 s com VIV
400 H1/3=3 m Tz=16 s sem VIV
H1/3=3 m Tz=16 s com VIV
200
0
-0.25 -0.2 -0.15 -0.1 -0.05 0 0.05 0.1 0.15 0.2
Deslocamento Direção Cross-Flow [m]
Figura B 16 – Comparação do efeito do acréscimo de VIV no envelope de deslocamento na
direção cross-flow
112
Comparação - Deslocamento Direção Axial
2000
1800
1600
P o siç ão no R iser [m ]
1400
1200
1000
800
600
H1/3=3 m Tz=8 s sem VIV
H1/3=3 m Tz=8 s com VIV
400 H1/3=5 m Tz=13 s sem VIV
H1/3=5 m Tz=13 s com VIV
200 H1/3=3 m Tz=16 s sem VIV
H1/3=3 m Tz=16 s com VIV
0
-2 -1.5 -1 -0.5 0 0.5 1 1.5
Deslocamento Direção Cross-Flow [m]
Figura B 17 – Comparação do efeito do acréscimo de VIV no envelope de deslocamento na
direção axial
1800
1600
Po siçã o n o R ise r [m ]
1400
1200
1000
800
600
H1/3=3 m Tz=8 s sem VIV
H1/3=3 m Tz=8 s com VIV
400 H1/3=5 m Tz=13 s sem VIV
H1/3=5 m Tz=13 s com VIV
200 H1/3=3 m Tz=16 s sem VIV
H1/3=3 m Tz=16 s com VIV
0
-0.5 0 0.5 1 1.5 2 2.5
Rotação [graus]
Figura B 18 – Comparação do efeito do acréscimo de VIV no envelope de rotação no plano
definido pelas direções vertical inline
113
Comparação - Rotação no Plano Definido pelas Direções Vertical e Transversal
2000
1800
1600
P o siç ão no R iser [m ]
1400
600
400
200
0
-0.08 -0.06 -0.04 -0.02 0 0.02 0.04 0.06 0.08
Rotação [graus]
Figura B 19 – Comparação do efeito do acréscimo de VIV no envelope de rotação no plano
definido pelas direções vertical e cross-flow
1800
1600
Po siçã o n o R ise r [m ]
1400
600
400
200
0
-60 -50 -40 -30 -20 -10 0 10 20 30 40
Força [kN]
Figura B 20 – Comparação do efeito do acréscimo de VIV no envelope de força na direção
inline
114
Comparação - Força Direção Cross-Flow
2000
1800
1600
P o siç ão no R iser [m ]
1400
600
400
200
0
-2 -1.5 -1 -0.5 0 0.5 1 1.5 2 2.5 3
Força [kN]
Figura B 21 – Comparação do efeito do acréscimo de VIV no envelope de força na direção
cross-flow
1400
1200
1000
800
600
400
200
0
-9000 -8000 -7000 -6000 -5000 -4000 -3000 -2000 -1000 0
Força [kN]
Figura B 22 – Comparação do efeito do acréscimo de VIV no envelope de força axial
115
Comparação - Momento Fletor no Plano Definido pelas Direções Vertical e Inline
2200
2000
1800
1600
P o siç ão no R iser [m ]
1400
1200
1000
800
1800
1600
Po siçã o n o R ise r [m ]
1400
600
400
200
0
-20 -15 -10 -5 0 5 10 15
Momento [kN*m]
Figura B 24 – Comparação do efeito do acréscimo de VIV no envelope de momento fletor
no plano definido pelas direções vertical e cross-flow
116
Comparação - Tensão devido ao Momento Fletor
2200
2000
1800
1600
P o siç ão no R iser [m ]
1400
1200
1000
800
1400
1200
1000
800
600
400
200
0
-50 0 50 100 150 200 250 300
Tensão [MPa]
Figura B 26 – Comparação do efeito do acréscimo de VIV no envelope de tensão devido à
força axial
117
Comparação - Tensão de Cisalhamento
2200
2000
1800
1600
P o siç ão no R iser [m ]
1400
1200
1000
800
2000
1800
1600
Po siçã o n o R ise r [m ]
1400
H1/3=3 m Tz=8 s sem VIV
H1/3=3 m Tz=8 s com VIV
1200
H1/3=5 m Tz=13 s sem VIV
H1/3=5 m Tz=13 s com VIV
1000
H1/3=3 m Tz=16 s sem VIV
H1/3=3 m Tz=16 s com VIV
800
600
400
200
0
0 100 200 300 400 500 600 700
Tensão [MPa]
Figura B 28 – Comparação do efeito do acréscimo de VIV no envelope de tensão de von
Mises
118
APÊNDICE C – CURVAS DE NÍVEL
Figura C 2- Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo
do riser; kθ = 686 kNm/grau, velocidade de avanço de 0,5 nó, direção de 90 graus
119
Figura C 3 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 686 kNm/grau, velocidade de avanço de 0,5 nó, direção de 45 graus
Figura C 4 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 686 kNm/grau, velocidade de avanço de 0,5 nó, direção de 0 graus
Figura C 5 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 686 kNm/grau, velocidade de avanço de 1,0 nó, direção de 180 graus
120
Figura C 6 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 686 kNm/grau, velocidade de avanço de 1,0 nó, direção de 135 graus
Figura C 7 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 686 kNm/grau, velocidade de avanço de 1,0 nó, direção de 90 graus
Figura C 8 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 686 kNm/grau, velocidade de avanço de 1,0 nó, direção de 45 graus
121
Figura C 9 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 686 kNm/grau, velocidade de avanço de 1,0 nó, direção de 0 graus
Figura C 10 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 686 kNm/grau, velocidade de avanço de 1,5 nó, direção de 45 graus
Figura C 11 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 338 kNm/grau, velocidade de avanço de 0,5 nó, direção de 180 graus
122
Figura C 12 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 338 kNm/grau, velocidade de avanço de 0,5 nó, direção de 135 graus
Figura C 13 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 338 kNm/grau, velocidade de avanço de 0,5 nó, direção de 90 graus
Figura C 14 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 338 kNm/grau, velocidade de avanço de 0,5 nó, direção de 45 graus
123
Figura C 15 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 338 kNm/grau, velocidade de avanço de 0,5 nó, direção de 0 graus
Figura C 16 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 338 kNm/grau, velocidade de avanço de 1,0 nó, direção de 180 graus
Figura C 17 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 338 kNm/grau, velocidade de avanço de 1,0 nó, direção de 135 graus
124
Figura C 18 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 338 kNm/grau, velocidade de avanço de 1,0 nó, direção de 90 graus
Figura C 19 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 338 kNm/grau, velocidade de avanço de 1,0 nó, direção de 45 graus
Figura C 20 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 338 kNm/grau, velocidade de avanço de 1,0 nó, direção de 0 graus
125
Figura C 21 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 338 kNm/grau, velocidade de avanço de 1,5 nó, direção de 180 graus
Figura C 22 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 338 kNm/grau, velocidade de avanço de 1,5 nó, direção de 135 graus
Figura C 23- Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 338 kNm/grau, velocidade de avanço de 1,5 nó, direção de 90 graus
126
Figura C 24 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser;
kθ = 338 kNm/grau, velocidade de avanço de 1,5 nó, direção de 45 graus
Figura C 25 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 338 kNm/grau, velocidade de avanço de 1,5 nó, direção de 0 graus
Figura C 26 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 169 kNm/grau, velocidade de avanço de 0,5 nó, direção de 180 graus
127
Figura C 27 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 169 kNm/grau, velocidade de avanço de 0,5 nó, direção de 135 graus
Figura C 28 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 169 kNm/grau, velocidade de avanço de 0,5 nó, direção de 90 graus
Figura C 29 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 169 kNm/grau, velocidade de avanço de 0,5 nó, direção de 45 graus
128
Figura C 30 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 169 kNm/grau, velocidade de avanço de 0,5 nó, direção de 0 graus
Figura C 31 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 169 kNm/grau, velocidade de avanço de 1,0 nó, direção de 180 graus
Figura C 32 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 169 kNm/grau, velocidade de avanço de 1,0 nó, direção de 135 graus
129
Figura C 33- Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 169 kNm/grau, velocidade de avanço de 1,0 nó, direção de 90 graus
Figura C 34 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 169 kNm/grau, velocidade de avanço de 1,0 nó, direção de 45 graus
Figura C 35 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 169 kNm/grau, velocidade de avanço de 1,0 nó, direção de 0 graus
130
Figura C 36 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 169 kNm/grau, velocidade de avanço de 1,5 nó, direção de 180 graus
Figura C 37- Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 169 kNm/grau, velocidade de avanço de 1,5 nó, direção de 135 graus
Figura C 38 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 169 kNm/grau, velocidade de avanço de 1,5 nó, direção de 90 graus
131
Figura C 39 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 169 kNm/grau, velocidade de avanço de 1,5 nó, direção de 45 graus
Figura C 40 - Curvas de nível de tensão de von Mises máxima no riser e inclinação no topo do riser; kθ
= 169 kNm/grau, velocidade de avanço de 1,5 nó, direção de 0 graus
132
APÊNDICE D – EFEITO
O DO COMPRIMENTO DO RISE
ER NOS
MAPAS OPERACIONAIIS
A mudança de locação de B
BOP suspenso por riser de perfuração obviamen
nte não está restrita a
uma configuração de riser. No quee diz respeito ao comprimento do riser, de maneira
m geral, risers
menores possuem janelas operacionnais maiores do que risers mais extensos (Sevillaano et al, 2013).
Os mapas operacionais abaixxo foram criados para o seguinte cenário:
133
Figura D-2 – Mapa Operacionnal – Parâmetros de Heave. Dupla amplitude siignificativa de heave
máxima em função de comprimentoo de riser e de velocidade de navegação.
Observa-se que a operação dde transferência do BOP pode ser realizada a vellocidades maiores na
medida em que o comprimento do rriser diminui, ou para condições mais severas dee carregamento.
Feitas as devidas consideraçções entre tempo gasto recolhendo-se o riser e tempo ganho com a
viagem a uma velocidade maior, com este tipo de mapa operacional seria po
ossível determinar o
comprimento ótimo para realizar a m
mudança de locação do BOP com o riser suspen
nso.
134
TAMENTO ESTÁTICO DO RISER
ANEXO A – COMPORT R
ds
As forças identificadas aginndo sobre o riser são: força axial ou tração de topo (T), força
cisalhante (V), força horizontal reesultante das pressões interna e externa (Fxo+F
+ xi), força vertical
resultante das pressões interna e exxterna (Fyo+Fyi), peso próprio do elemento (W
W) e força de arrasto
devido à correnteza (N), que é deccomposta em dois componentes, um normal e outro tangencial ao
elemento, sendo considerado apenass o efeito da componente normal.
As equações de equilíbrioo do elemento de riser são escritas a partirr do somatório dos
componentes das forças nas direçõees x e y:
135
Direção x: F =0
Direção y: F =0
Então se tem:
cos dθ = 1 sen dθ = dθ
F =0
F =0
136
Tcos θ + Vsen θ dθ + d Tsen θ − dVcos θ + F +F − Nrdθ cos θ − W = 0 (A6)
As forças Fx0, Fxi, Fy0 e Fyi originam-se da pressão hidrostática (Patel, Witz, 1991) e podem ser
expressas por:
W = γ A rdθ (A10)
onde γs é o peso especifico do material do riser e A área da seção transversal da parede do riser.
Para efeitos de simplificação, a pressão dinâmica exercida pelo escoamento do fluido interno será
desprezada. As forças hidrostáticas não atuam na extremidade inferior do riser por esta estar engastada
e, portanto, sem realizar deslocamentos.
θ ≠ dθ ds ≠ dx
Então:
137
dy
dx dy dy dx
cos θ = = 1+ sen θ = =
ds dx ds
dy
1+
dx
dy d y
dθ 1 d arctan dx dx dx
= = =
ds r dx ds
dy
1+
dx
dV dV dx dV dy
= = 1+ (A12)
ds dx ds dx dx
d y dy dV dy
T + P A − PA 1+ − + γ A −γ A −γ A −N 1+ =0
dx dx dx dx
(A13)
dV d d x
= EI
dy dy dy
(A14)
d d x d x dx
EI − T + P A − PA =N+ γ A +γ A −γA =0
ds dy dy dy
(A15)
A equação (A15) representa a equação diferencial que governa a deflexão lateral estática do riser.
138
ANEXO B – SOLUÇÃO DA EQUAÇÃO DINÂMICA DO RISER
Métodos de integração no tempo têm como característica fundamental aproximar o valor das
derivadas que surgem nos sistema de equações de movimento e gerar uma solução passo a passo com
intervalo de tempo ∆t. A solução obtida para os deslocamentos no final de cada intervalo de tempo é
utilizada nos cálculos efetuados para o intervalo seguinte.
∆ + ∆ + 1+ R ∆ − F ∆ + R − F =0 (B1)
F ∆ = K ∆ ∆x + F (B2)
139
A partir deste ponto, as aproximações de Newmark são introduzidas na equação B4, assim como
seu procedimento iterativo. O método Newton Raphson também é utilizado.
∆ = + 1− ∆ x + ∆ ∆ (B5)
1
= +∆ +∆ − +∆
∆
2 ∆ (B6)
1 1 1
x = − − − − x
∆
β∆ ∆
β∆ 2 (B7)
γ γ 1 ∆
x = + 1− ∆ x + − + −γ − x
β∆ β 2 β
(B8)
∆ ∆
1 1 1
− − −− x
β∆ ∆
β∆ 2
γ γ
(B9)
+ + 1− + 1− x∆ + =
β∆ β 2
Onde
140
1 γ
a = ; a =
β∆ β∆
1 γ
a = ; a = −1
β∆ 2 (B10)
γ γ
a = −1 ; a = −1 ∆
2
a + a + =
1+α R − F − −a −a x − −a −a x −α R − F
(B11)
∆ ∆
∆ = R ∆ (B12)
A solução desta equação corresponde a solução de um problema não linear estático e que pode
ser resolvida através de método incremental-iterativo. Após a convergência do processo não linear, e de
posse do vetor ∆x retorna-se aos operadores de integração para calcular as velocidades e acelerações no
intervalo de tempo t + ∆t
x ∆ = a ∆ −a −a x (B13)
x ∆ = +a x +a x ∆ (B14)
Sendo, a = 1 − γ ∆t e a = γ∆t
141