Tes Is Julian Guzman

CENTRO DE INVESTIGACIN Y DE ESTUDIOS AVANZADOS
DEL INSTITUTO POLITCNICO NACIONAL

DEPARTAMENTO DE COMPUTACIN
Modelo de Vibracin y Simulacin de Comportamiento

Dinmico de Riser, con Interaccin Agua, Estructura y
Suelo
Tesis que presenta:
Julin Edgardo Guzmn Corts
Para obtener el grado de
Maestro en Ciencias en Computacin
Director de la Tesis: Dr. Sergio Vctor Chapa Vergara
Codirector de la Tesis: Dr. Amilcar Meneses Viveros
Mxico, D. F.
Febrero, 2010
II
III
Agradecimientos
A lo largo de este trabajo de investigacin y durante mis estudios de maestra, he sido

afortunado de contar con el apoyo de diversas personas, a las cuales quisiera agradecer
por impulsarme a cumplir una ms de mis metas.
En primer lugar, quisiera agradecer a Dios por permitirme llegar a este momento
tan importante de mi vida profesional.
A mis padres que siempre me han apoyado incondicionalmente en mis decisiones y
me alientan para seguir adelante. Gracias por su amor, comprensin y apoyo que me
brindaron durante este tiempo que estuve lejos de casa.
A mis abuelitos que me brindan su cario y siempre estn al pendiente de mis pasos.
A mis tos que siempre me escuchan y estn presentes cuando necesito su ayuda.
Gracias por sus consejos que me han ayudado a tomar las mejores decisiones.
A todos mis viejos y nuevos amigos que conoc en el CIVESTAV. Gracias por compartir tantas aventuras y desveladas de estudio durante este tiempo.
A Pamela Snchez por apoyarme y tolerarme en los momentos difciles. Gracias
por permitirme compartir contigo tantos momentos y por hacer que la estancia en la
escuela, lejos de mi familia, fuera menos difcil.
A mi asesor, el Dr. Sergio V. Chapa Vergara le agradezco la oportunidad que me
ofreci para trabajar con l. Gracias por todo el apoyo y su ayuda para poder realizar
este trabajo de tesis. He aprendido mucho bajo su asesora.
Al Dr. Amilcar Meses Viveros le agradezco su ayuda y el tiempo que me brindo
durante el desarrollo de esta tesis. Gracias por la paciencia que tuvo para explicarme
los temas que me eran desconocidos, pero sobretodo, gracias por su amistad.
A mis sinodales: Dr. Jos Matas Alvarado y Dr. Oliver Schtze gracias por su
consejos, comentarios y sugerencias que me ayudaron a mejorar mi trabajo de investigacin.
Quisiera agradecer a todos los doctores del departamento de computacin, por todas
sus enseanza, consejos y por compartir sus conocimientos conmigo.
A Sofa Reza por brindarme su amistad y su ayuda para realizar todos los tramites
que se requieren durante estancia de la maestra. Gracias por su trato tan gentil y por
hacernos sentir como en casa.
Finalmente quisiera agradecer al CINVESTAV por permitirme realizar mis estudios
de maestra en la institucin y al Consejo Nacional de Ciencia y Tecnologa (CONACyT)
por el apoyo econmico que me permiti solventar los gastos durante mis estudios.
IV
Contenido
ndice de Figuras
IX
ndice de Tablas
XII
ndice de Algoritmos
XIV
Resumen
Abstract
Introduccin
Antecedentes
2.1 El petrleo en aguas profundas . . . . . . . . . . . . . . . . . . . . . . . .
2.1.1 Importancia del petrleo . . . . . . . . . . . . . . . . . . . . . . . .
2.1.2 El petrleo en Mxico . . . . . . . . . . . . . . . . . . . . . . . . .
2.1.3 Exploracin y explotacin de petrleo en aguas profundas . . . .
2.1.4 Retos y dificultades en la exploracin y explotacin en aguas profundas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.2 Desarrollo de campos petroleros en alta mar . . . . . . . . . . . . . . . .
2.2.1 Sistema submarino . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.2.2 Lneas de flujo, tuberas y sistemas riser . . . . . . . . . . . . . . .
2.2.3 Plataformas fijas y flotantes . . . . . . . . . . . . . . . . . . . . . .
2.2.4 Sistema de procesamiento superior . . . . . . . . . . . . . . . . . .
2.3 Modelos oceanogrficos . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.3.1 Base de datos oceanogrfica . . . . . . . . . . . . . . . . . . . . . .
2.3.2 Mtodo geostrfico . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.4 Simulacin . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.4.1 El proceso de simulacin por computadoras . . . . . . . . . . . .
2.5 Vibracin . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.5.1 Ejemplos de estructuras afectadas por el fenmeno de vibracin .
2.6 Trabajos relacionados . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.7 Resumen . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
V
9
9
9
12
13
16
17
17
19
21
22
22
23
23
25
26
27
28
30
31
VI
3
CONTENIDO
Modelo del sistema
3.1 Modelo fsico del sistema . . . . . . . . . . . . . . . . . . . . . .
3.2 Movimiento Armnico Simple (MAS) . . . . . . . . . . . . . . .
3.2.1 Sistema de un grado de libertad . . . . . . . . . . . . . .
3.2.2 Sistema de varios grados de libertad . . . . . . . . . . . .
3.3 Modelo matemtico . . . . . . . . . . . . . . . . . . . . . . . . . .
3.3.1 Modelo matemtico con interaccin a primeros vecinos .
3.3.2 Modelo matemtico con interaccin a segundos vecinos
3.4 Resumen . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Diseo del sistema
4.1 Diseo Conceptual . . . . . . . . . . . . . . .
4.2 Modelo del riser (cadenas de vibracin) . . .
4.2.1 Problema de valores propios . . . . .
4.2.2 Anlisis de Fourier . . . . . . . . . . .
4.3 Modelo de corrientes basado en Monte Carlo
4.3.1 Mtodo Monte Carlo . . . . . . . . . .
4.3.2 Distribucin normal . . . . . . . . . .
4.4 Visualizacin . . . . . . . . . . . . . . . . . . .
4.5 Resumen . . . . . . . . . . . . . . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
Desarrollo e implementacin del sistema

5.1 Arquitectura general . . . . . . . . . . . . . . . . . . . . . . . . .
5.1.1 Plataforma del subsistema de cmputo de datos . . . . .
5.1.2 Plataforma del subsistema de visualizacin . . . . . . . .
5.2 Desarrollo del subsistema de cmputo de datos . . . . . . . . . .
5.2.1 Proceso para obtener las matrices de movimiento . . . .
5.2.2 Proceso para obtener los valores y vectores propios . . .
5.2.3 Proceso para obtener modos normales de vibracin . . .
5.2.4 Proceso para dividir el riser en segmentos . . . . . . . . .
5.2.5 Proceso para generar frecuencias aleatorias . . . . . . . .
5.2.6 Proceso para calcular los modos aleatrios de vibracin .
5.3 Desarrollo del subsistema para la interfaz grfica . . . . . . . . .
5.3.1 Mdulo de grficas de resultados . . . . . . . . . . . . . .
5.3.2 Mdulo de animacin de partculas 2D . . . . . . . . . .
5.3.3 Mdulo de visualizacin interactiva 3D . . . . . . . . . .
5.4 Resumen . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Validacin y pruebas del simulador
6.1 Descripcin general del uso del simulador .
6.1.1 Aplicacin de cmputo de datos . .
6.1.2 Aplicacin de visualizacin . . . . .
6.2 Diseo de la pruebas . . . . . . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
33
34
35
37
37
38
39
42
45
.
.
.
.
.
.
.
.
.
47
48
49
50
51
52
53
55
56
57
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
59
59
60
60
61
62
63
64
65
67
70
72
74
75
77
78
.
.
.
.
79
79
80
83
87
CONTENIDO
6.3
6.4
6.5
7
VII
Resultados de las pruebas . . . . . . . . . . . .

6.3.1 Prueba con 8 segmentos . . . . . . . . .
6.3.2 Prueba con 16 segmentos . . . . . . . .
6.3.3 Prueba con 32 segmentos . . . . . . . .
Comparacin de resultados con otros modelos
Resumen . . . . . . . . . . . . . . . . . . . . . .
Resultados, discusin, y conclusiones

7.1 Resultados . . . . . . . . . . . . .
7.2 Discusin . . . . . . . . . . . . . .
7.3 Conclusiones . . . . . . . . . . . .
7.4 Trabajo futuro . . . . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
88
89
91
94
96
99
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
101
101
102
104
105
Bibliografa
106
A Instalacin del sistema

A.1 Requerimientos de Instalacin . . . . . . . . . . . . . . . . . . . . .
A.1.1 Requerimientos para el subsistema de Cmputo de Datos
A.1.2 Requerimientos para el subsistema de Visualizacin . . . .
A.2 Tutorial de compilacin . . . . . . . . . . . . . . . . . . . . . . . .
A.2.1 Subsistema de Cmputo de Datos . . . . . . . . . . . . . .
A.2.2 Subsistema de Visualizacin . . . . . . . . . . . . . . . . .
111
111
111
112
112
112
113
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
VIII
CONTENIDO
Lista de Figuras
2.1
2.2
2.3
2.4
2.5
2.6
2.7
2.8
2.9
2.10
2.11
2.12
2.13
2.14
2.15
Proceso de formacin del petrleo . . . . . . . . . . . . . . .

Pozo de alquitrn . . . . . . . . . . . . . . . . . . . . . . . .
Principales usos del petrleo . . . . . . . . . . . . . . . . . .
Regiones de aguas profundas . . . . . . . . . . . . . . . . . .
Superficie de las Aguas Profundas en el Golfo de Mxico . . .
Procesos para desarrollar un proyecto en Aguas Profundas . .
Elementos para el desarrollo de campos en alta mar . . . . . .
Clasificacin de sistemas submarinos . . . . . . . . . . . . .
Lnea de flujo, tuberas y sistemas riser . . . . . . . . . . . . .
Tipos de sistema de riser para aguas profundas . . . . . . . .
Tipos de plataformas que se utilizan en los campos petroleros
Puente Tacoma . . . . . . . . . . . . . . . . . . . . . . . . .
Puente del Milenio, Londres. . . . . . . . . . . . . . . . . . .
Hotel Burj Al Arab, Dubai. . . . . . . . . . . . . . . . . . . .
Tipos de riser utilizados en los campos de aguas profundas. . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
10
10
11
14
15
16
18
19
20
21
22
28
28
29
30
3.1
Representacin del modelo fsico. Tubera sumergida en agua dividida en n

partculas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
3.2 Movimiento del riser provocado por las corrientes marinas . . . . . . . . . . .
3.3 Grfica de Movimiento Armnico Simple . . . . . . . . . . . . . . . . . . . .
3.4 Sistema de un grado de libertad . . . . . . . . . . . . . . . . . . . . . . . . .
3.5 Transformacin del modelo fsico al modelo matemtico . . . . . . . . . . . .
3.6 Cadena de vibracin con interaccin a primeros vecinos . . . . . . . . . . . .
3.7 Cadena de vibracin con interaccin a primeros vecinos con partculas imaginarias en los extremos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
3.8 Ecuaciones de movimiento con interaccin a primeros vecinos . . . . . . . . .
3.9 Cadena de vibracin con interaccin a segundos vecinos . . . . . . . . . . . .
3.10 Cadena de vibracin con interaccin a segundos vecinos con partculas imaginarias en los extremos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
4.1
4.2
4.3
Divisin del simulador en dos subsistemas . . . . . . . . . . . . . . . . . . .

Mdulos que componen el subsistema de cmputo de datos . . . . . . . . . .
Mdulos que componen el subsistema de Visualizacin . . . . . . . . . . . . .
IX
34
35
36
37
38
39
40
41
43
44
48
48
49
LISTA DE FIGURAS
4.4
4.5
4.6
49
53
4.7
4.8
Representacin general del sistema . . . . . . . . . . . . . . . . . . . . . . . .

Divisin del riser en partculas y segmentos . . . . . . . . . . . . . . . . . . .
Generacin de frecuencias aleatrias para cada partcula, con diferenes condiciones en cada segmento . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Campa de la distribucin normal. . . . . . . . . . . . . . . . . . . . . . . . .
Distribucin normal con rangos de 1 a 1. . . . . . . . . . . . . . . . . . . . .
5.1
5.2
5.3
5.4
5.5
5.6
5.7
5.8
5.9
5.10
5.11
5.12
Arquitectura general del sistema . . . . . . . . . . . . . . . . . .

Diagrama de flujo de datos del Subsistema de Cmputo de Datos
Asignacin de partculas por segmento . . . . . . . . . . . . . .
Clculo de la superposicin de frecuencias aletorias . . . . . . . .
Diagrama de clases del sistema de visualizacin . . . . . . . . . .
Procesos del mdulo de grficas de resultados . . . . . . . . . . .
Simulacin del mdulo de grficas . . . . . . . . . . . . . . . . .
Interfaz para generar las grficas de las partculas . . . . . . . . .
Procesos del mdulo de animacin de partculas . . . . . . . . .
Animacin del partculas . . . . . . . . . . . . . . . . . . . . . .
Procesos del mdulo de visualizacin interactiva . . . . . . . . .
Aimacin del sistema en 3D . . . . . . . . . . . . . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
60
61
67
70
73
74
75
75
76
77
77
78
6.1
6.2
6.3
6.4
6.5
6.6
6.7
6.8
6.9
6.10
6.11
6.12
6.13
6.14
6.15
6.16
6.17
6.18
6.19
6.20
6.21
6.22
6.23
Pantalla de inicio de la aplicacin de cmputo de datos . . . . . . . . . . . . .

Pantalla de error en el nmero de segmentos . . . . . . . . . . . . . . . . . .
Pantalla de la creacin del nombre del directorio . . . . . . . . . . . . . . . .
Pantalla de matrices de movimiento con sus valores propios y vector propio . .
Pantalla del clculo de modos normales de vibracin . . . . . . . . . . . . . .
Pantalla de los datos de media y varianza de los segmentos . . . . . . . . . . .
Pantalla de inicio del subsistema de visualizacin . . . . . . . . . . . . . . . .
Mens del subsistema de visualizacin . . . . . . . . . . . . . . . . . . . . .
Barra de herramientas del subsistema de visualizacin . . . . . . . . . . . . .
Pantalla para abrir archivo. . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Pantalla de inicio de simulacin . . . . . . . . . . . . . . . . . . . . . . . . .
Pantalla del proceso de simulacin. . . . . . . . . . . . . . . . . . . . . . . .
Pantalla del sistema que muestra la grfica de los datos. . . . . . . . . . . . .
Distribucin de partculas en el riser para las pruebas . . . . . . . . . . . . . .
Archivos que se generaron de la prueba de 8 segmentos . . . . . . . . . . . . .
Resultado de la visualizacin de los modos aleatorios para primeros vecinos . .
Resultado de la visualizacin de los modos aleatorios para segundos vecinos .
Resultado de la visualizacin de la superposicin para primeros vecinos . . . .
Resultado de la visualizacin de la superposicin para segundos vecinos . . .
Archivos que se generaron de la prueba de 16 segmentos . . . . . . . . . . . .
Visualizacin de los modos aleatorios con 16 segmento para primeros vecinos .
Visualizacin de los modos aleatorios con 16 segmento para segundos vecinos
Visualizacin de la superposicin con 16 segmentos para primeros vecinos . .
80
81
81
82
82
83
84
84
85
85
86
86
87
88
89
89
90
90
91
92
92
92
93
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
54
55
56
LISTA DE FIGURAS
6.24
6.25
6.26
6.27
6.28
6.29
6.30
Visualizacin de la superposicin con 16 segmentos para segundos vecinos . .

Archivos que se generaron de la prueba de 32 segmentos . . . . . . . . . . . .
Visualizacin de los modos aleatorios con 32 segmento para primeros vecinos .
Visualizacin de los modos aleatorios con 32 segmentos para segundos vecinos
Visualizacin de la superposicin con 32 segmentos para primeros vecinos . .
Visualizacin de la superposicin con 32 segmentos para segundos vecinos . .
Grfica de comparacin de frecuencias naturales del riser . . . . . . . . . . . .
XI
93
94
95
95
96
96
99
XII
LISTA DE FIGURAS
Lista de Tablas
2.1
Principales pases productores de petrleo en aguas profundas . . . . . . . . .
14
6.1
6.2
6.3
6.4
6.5
6.6
6.7
6.8
6.9
6.10
6.11
6.12
Datos de entrada para la ejecucin de la prueba de 8 segmentos . . . . . . . .

Resultados de los modos aleatorios de vibracin. Riser con 8 segmentos. . . . .
Resultados de la superposicin de los movimientos. Riser con 8 segmentos . .
Resultados de los modos aleatorios de vibracin para la prueba de 16 segmentos
Resultados de la superposicin de los movimientos. Riser con 16 segmentos . .
Resultados de los modos aleatorios de vibracin. Riser con 32 segmentos . . . .
Resultados de la superposicin de movimientos. Riser con 32 segmentos . . . .
Parmetros del riser de ejemplo . . . . . . . . . . . . . . . . . . . . . . . . . .
Parmetros del ejemplo para nuestro modelo . . . . . . . . . . . . . . . . . .
Comparacin de resultados con otros modelos . . . . . . . . . . . . . . . . . .
89
90
91
91
93
94
94
95
96
97
98
98
XIII
XIV
LISTA DE TABLAS
Lista de Algoritmos
1
2
3
4
5
Genera Matriz de Movimiento . . . . .

Modos Normales de Vibracin . . . . .
Divisin de Segmentos . . . . . . . . . .
Frecuencias aleatorias con Monte Carlo
Movimientos Aleatorios de Vibracin .
XV
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
63
65
66
68
71
XVI
LISTA DE ALGORITMOS
Resumen
ctualmente existen diversas formas y fuentes proveedoras de energa, y el petrleo

es la ms importante. Este energtico porporciona el mayor porcentaje de energa que
se consume en el mundo y se utiliza como materia prima para fabricar una gran variedad de productos de consumo cotidiano. Histricamente el crecimiento econmico de
los pases industrializados depende directamente del abastecimiento de petrleo que
tengan. Por esta razn, los principales pases productores de petrleo se han dado a
la tarea de realizar la exploracin y explotacin de petrleo en aguas profundas. Esto
representa una gran oportunidad para nuestro pas, debido al potencial de petrleo
que podra existir en aguas profundas del Golfo de Mxico.
Llevar a cabo el proceso de exploracin y explotacin de petrleo en aguas profundas no es sencillo, ya que se enfrentan retos y dificultades diferentes a los que se
presentan al trabajar en aguas someras. Uno de los problemas ms importantes, es que
las estructuras que se utilizan para la extraccin del petrleo sean capaces de resistir la
interaccin de las corrientes marinas.
En esta investigacin se presenta el diseo y desarrollo de un simulador que ayuda
a realizar un anlisis de las vibraciones que se presentan en los riser, provocadas por
la fuerza de las corrientes marinas. Este simulador permite construir un modelo representativo del sistema para conocer o predecir el comportamiento dinmico de las
estructuras en diversas situaciones a partir de bases tericas. Los resultados obtenidos
nos ayudan a conocer los puntos crticos de la estructura y permiten tomar desiciones
al momento de su construccin para hacerlas ms flexibles y resistentes.
El desarrollo del simulador esta basado en un modelo matemtico de cadenas de vibracin. Adems, cuenta con un modelo estocstico de corrientes basado en el mtodo
Monte Carlo que nos permite analizar situaciones complejas y que son difciles de reproducir fsicamente. El sistema se compone de dos subsistemas: el subsistema de
cmputo de datos que se encarga de utilizar los modelos desarrollados para obtener
los resultados del comportamiento del sistema. Por su parte, el subsistema de visualizacin utiliza los datos obtenidos para generar la simulacin grfica que nos permita
interpretar ms facilmente los resultados.
Finalmente el simulador fue sometido a diversas pruebas variando las condiciones
iniciales del sistema. De esta manera, se verific que tiene una presicin adecuada y se
garantiza que los resultados estan apegados al modelo matemtico.
Abstract
owadays there are various sources and forms of energy suppliers, and oil is the
most important one. Oil provides the highest percentage of energy consumed in the
world and is used as a raw material to manufacture a wide variety of daily consumed
pruducts. Historically the economical growth of industrialized countries depends directly on its oil supply. For this reason, the major oil producers countries have been
dealing with the task of performing the exploration and exploitation of oil in deepwater. This represents a great opportunity for our country, based on the posibility of oil
existence in deepwater of Gulf of Mexico.
The process of exploration and exploitation of oil in deepwater is not simple, mostly
because different conditions founded in contrast to the usual one when working in
shallow water. One of the most important problems is, the characteristics of resistance
related to the structures used for oil extraction against the interaction of ocean currents.
This research introduces the desing and development of a simulator that helps to
analyze the vibrations in the riser, caused by the forces of ocean current. This simulator
allows the construction of a representative model of the system in order to know or
predic the dynamic behavior of the structures in diverse situations from theoretical
bases. The results help us to know the critical points of the structure allowing us to
take decisions at their construction time to make them more flexible and resistant
The development of the simulator is based on a mathematical model of strings vibration. Moreover, it has a stochastic model of currents based on the Monte Calo method,
whole allowing us to analyze complex situations which are difficult to reproduce physically. The system consist of two subsystems: The computing subsystem which is responsible for the use of the models developed to get the system behavior results. Moreover, The display subsystem uses the data to generate graphical simulation allowing
us more easily interpret the results.
Finally, the simulator was subject to several test varying the initial conditions of the
system. Thus, we verified an adequate precision of the simulator ensuring that the
results are related to the mathematical model.
Captulo 1
Introduccin
No existe en el mundo nada ms
poderoso que una idea a la que le ha
llegado su tiempo
Victor Hugo
n nuestros das, el petrleo se ha convertido en el principal energtico de la humanidad. Actualmente su uso se ha extendido y para satisfacer la demanda, su produccin va en aumento. Debido a que este energtico es un recurso no renovable,
algunos pases han buscado alternativas para tener una mayor produccin de petrleo. Se ha descubierto que las mayores reservas petrolferas se encuentran en zonas
conocidas como aguas profundas.
Las aguas profundas representan una gran oportunidad para nuestro pas. En Mxico las aguas profundas se encuentran ubicadas en las aguas territoriales del Golfo de
Mxico. En esta zona a partir de informacin geolgica y geofsica, PEMEX ha identificado la existencia de estructuras geolgicas capaces de contener hidrocarburos.
Sin embargo, el reto de explorar y explotar yacimientos de petrleo en aguas profundas no es sencillo, ya que es totalmente diferente a realizarlo en tierra o cerca de las
costas (aguas someras) como actualmente se hace en nuestro pas. Este es un proceso
largo, complicado y muy costoso; ya que se tiene que contar con el financiamiento,
la infraestructura, y los conocimientos necesarios para llevarlo a cabo. Trabajar en estas condiciones lleva a las empresas exploradoras a enfrentarse a nuevos retos como
son: las corrientes marinas, los cambios de temperatura entre la superficie y el lecho
marino, las presiones anormales al momento de la perforacin, la geometra de los yacientos, la operacin remota a travs de robots y desarrollar nueva tecnologa que a
pases pioneros en el proceso de explotacin de yacimientos en aguas profundas (Estados Unidos, Brasil, Noruega, etc.), les ha tomado varios aos de desarrollo.
En nuestro pas, el lapso del tiempo de estudio que toma decidir excavar un nuevo
pozo es de entre 12 a 18 meses aproximadamente, lo que repercute en los costos de
operacin del proyecto [1]. Durante este periodo se hacen diversos estudios, pruebas y
5
6
simulaciones sobre cada una de las partes que conforman el proyecto de explorar y explotar yacimientos de crudo. El tener un control o una prediccin sobre las situaciones
de riesgo, agilizar este proceso y disminuir el tiempo para la puesta en produccin de
los campos petroleros, permite aumentar considerablemente las utilidades.
En esta investigacin abordamos el problema de poder predecir el comportamiento
dinmico de los riser (estructuras mecnicas que se utilizan para la explotacin de petrleo en aguas profundas), cuyas estructuras se ven afectadas directamente por las
corrientes marinas.
Las diversas condiciones de fuerzas que las corrientes marinas ejercen sobre las estructuras producen un movimiento continuo, lo que ocasiona que los riser presenten
vibraciones a lo largo de su cuerpo de metal. Poder predecir esta clase de comportamiento es de gran importancia, ya que si la frecuencia del movimiento provocado
por las corrientes se encuentra cerca de la frecuencia natural de la estructura, se puede
caer en Resonancia [2]. Este fenmeno produce desgastes y deformaciones que provocan daos en la estructura, que en el peor de los casos pueden ocasionar accidentes.
Cuando se trabaja con sistemas en los cuales se tienen condiciones que pueden llegar
a ser impredecibles (como es el caso de las corrientes marinas); se tiene un alto indice
de riesgos. Actualmente se cuenta con herramientas como simuladores y modelos
matemticos de vibracin que nos pueden auxiliar en el estudio del comportamiento
de las estructuras mecnicas. Sin embargo, este tipo de tecnologas necesarias para la
exploracin y explotacin de yacimientos en aguas profundas actualmente son muy
costosas, debido a que son escasas y su disponibilidad en los mercados es limitada.
El desarrollo de simuladores estocsticos que utilicen mtodos numricos y nuevos
algoritmos de computacin, son de las herramientas principales que pueden ayudar a
solucionar diversos problemas en el proceso de exploracin y explotacin de petrleo
en aguas profundas, especficamente predecir el comportamiento dinmico de los riser
en operacin, provocado por las fuertes corrientes marinas.
El objetivo principal de este trabajo de investigacin es presentar el diseo y desarrollo de un simulador de vibraciones con un modelo estocstico de corrientes, que
nos ayude a predecir el comportamiento dinmico de los riser, tomando en consideracin la interaccin con su ambiente (Agua, Estructura y Suelo). Para su desarrollo se
siguieron los siguientes pasos:
Modelo fsico del sistema: El modelo fsico del sistema se presenta como una tubera de metal sumergida en el mar, sobre la cual la fuerza de las corrientes marinas
producen vibraciones. Esta tubera se divide en partculas de masa equidistantes sobre
las cuales se analiza la fuerza que ejercen las corrientes marinas. El sistema se representa como una cadena de vibracin unidimensional de n grados de libertad formada
por una serie de partculas y unidas por medio de resortes.
Modelo Matemtico: El modelo fsico se transform a un modelo matemtico de
valores propios, basado en Movimiento Armnico Simple homogneo MAS [3]. Este
modelo permite obtener las matrices de movimiento para la interaccin de partculas
de masa del riser con interaccin a primeros y segundos vecinos. MAS sirve como base
para nuestro simulador y nos proporciona confiabilidad en nuestros resultados [4]. La
CAPTULO 1. INTRODUCCIN
solucin del modelo matemtico se obtiene de forma numrica. Para ello, se usan bibliotecas computacionales que nos ofrecen soluciones con diversos mtodos como el
de Kramer, el mtodo de Jacobi, algoritmo QR, etc. Particularmente se utiliz la biblioteca IMSL (FORTRAN Numerical Library) la cual, es una biblioteca computacional con
algortmos matemticos y estadsticos para lenguaje FORTRAN. La biblioteca recibe
como entrada la matriz de movimiento de nuestro modelo y nos ofrece como resultado
las frecuencias (valores propios) y las Amplitudes (vectores propios) correspondientes.
Generacin de Frecuencias Aleatorias: Para la reproduccin de fenmenos naturales (como las corrientes marinas), es necesario la simulacin de nmeros aleatorios.
De esta manera, el simulador cuenta con la implementacin de un mtodo numrico
no determnistico conocido como Monte Carlo (basado en una distribucin normal),
el cual nos ayuda a generar las condiciones de fuerzas de excitacin y frecuencias que
afectan a las estructuras mecnicas. Gracias a esto, vamos a poder analizar situaciones
que son muy costosas y difciles de producir fsicamente.
Anlisis de coeficientes de Fourier: Para obtener el comportamiento general de
nuestro sistema, calculamos la superposicin de las ondas que se obtuvieron como
resultado de nuestro modelo matemtico en cada uno de los nodos que conforman el
riser.
Desarrollo del sistema: El sistema se encarga de realizar el clculo de las vibraciones en la estructura para posteriormente desplegarlos en una simulacin grfica. El
desarrollo del sistema se encuentra dividido en 2 mdulos principales:
1. Mdulo de cmputo de datos: El mdulo principal se encarga de recibir los
valores de entrada para nuestro simulador (longitud del riser, masa, elasticidad, nmero de partculas, etc). Con estos valores se generan las matrices de
movimiento correspondientes a los modelos de primeros y segundos vecinos.
Adems, el mdulo se encarga de generar las frecuencias aleatorias (por medio
del mtodo de Monte Carlo) y realizar los clculos correspondientes para obtener
el comportamiento de nuestro sistema. Los resultados se almacenan en archivos
que posteriormente son utilizados en el mdulo de la Interfaz Grfica para desplegar la simulacin. Esta aplicacin fu desarrollada en lenguaje FORTRAN
(F90) y funciona en plataformas tipo Unix.
2. Interfaz Grfica: Este mdulo es una aplicacin totalmente independiente a la del
mdulo principal, es decir, que la representacin grfica de los resultados no se
realiza en tiempo real. Esta aplicacin hace uso de los datos generados por el mdulo principal y despliega los resultados en base a los datos contenidos en ellos.
La interfaz grfica es amigable para el usuario, adems permite que los resultados se puedan interpretar fcil y rpidamente. La interfaz grfica se encuentra
dividida en tres secciones: La primer seccin es una simulacin que nos indica
el comportamiento individual de cada uno de los nodos que conforman el sistema. La segunda seccin nos genera las grficas correspondientes a cada uno de
los nodos por separado. Y la tercer secccin nos ofrece una simulacin en 3D de
nuestro sistema. El desarrollo de la interfaz grfica est escrita en lenguaje C++
8
con ayuda de la biblioteca QT para el desarrollo de la interfaz grfica, adems se
usa OpenGL para la simulacin 3D.
Para verificar los resultados y la buena presicin del simulador, se realizaron diversas pruebas variando las condiciones iniciales del sistema. Estos resultados nos
permitieron determinar que nuestro simulador se encuentra apegado a su modelo
matemtico y que su margen de error es despreciable. Sin embargo, para comprobar totalmente nuestro simulador es necesario llevar a cabo una prueba fsca, la cual es
conocida como Prueba de Tanque; desafortunadamente este tipo de pruebas requieren
de equipo especializado y son demasiado costosas, por lo que esta prueba queda fuera
del alcance del trabajo de investigacin.
El presente documento se encuentra organizado de tal forma que el lector pueda
seguir cada uno de los pasos que fueron necesarios para el diseo y desarrollo de este
simulador. La informacin se encuentra organizada en 7 captulo que se describen a
continuacin.
En el Captulo 2: Antecedentes se ofrece un panorama general sobre la situacin
actual de la exploracin y explotacin de petrleo en aguas profundas, as como las
tcnicas y tecnologas que se usan actualmente para realizar ese proceso. Adems, se
describe la importancia de la simulacin en este tipo de sistemas y se presentan algunos
trabajos relacionados.
El Captulo 3: Modelo del sistema describe la respresentacin del modelo fsico
del sistema. Se construye el modelo matemtico con interaccin a primeros y segundos
vecinos para finalmente obtener nuestras matrices de movimiento.
En el Captulo 4: Diseo del simulador se describen los principales mdulos que
forman el sistema y las metodologas que se utilizaron para desarrollar cada uno de
ellos.
En el Captulo 5: Desarrollo e implementacin se muestra la arquitectura general
de nuestro simulador y se detallan los mdulos de programacin que se realizaron
durante su desarrollo.
En el Captulo 6: Validacin y pruebas se define la selecin de pruebas que se
realizaron para validar los resultados de nuestro simulador y se hace una comparacin
de resultados con otros modelos.
En el Captulo 7: Resultados, discusin y conclusiones se describen los resultados
obtenidos durante este trabajo de investigacin y se hace una discusin final haciendo
referencia a los resultados que se obtuvieron durante la etapa de pruebas. Adems, se
describen las conclusiones a las que se llegaron y se presenta una lista sobre los trabajos
futuros para mejorar el simulador.
Finalmente, el documento presenta un anexo en donde se detallan los requirimientos
necesarios para instalar y compilar los subsistemas que forman el simulador.
Captulo 2
Antecedentes
Vale ms saber alguna cosa de todo,
que saberlo todo de una sola cosa.
Blaise Pascal
ara desarrollar este trabajo de tesis, se realiz una investigacin sobre los
antecedentes, y el proceso general de exploracin y explotacin de petrleo.
En este captulo se abordan temas referentes a la importancia del petrleo en el mundo y especialmente en Mxico. Adems, se analizan las tcnicas actuales que se utilizan para explorar y explotar yacimientos de petrleo
en aguas profundas y los problemas a los que se enfrentan las empresas al
realizar este proceso.
Por otra parte, en este catulo se describe el concepto de simulacin por
computadora y la forma en como nos apoya para enfrentar problemas como
las vibraciones que se presentan en los riser. Finalmente se mencionan algunos desarrollos de software que se utilizan en la industria petrolera para
administrar y monitorear los sistemas riser.
2.1
2.1.1
El petrleo en aguas profundas

Importancia del petrleo
El petrleo es una forma de energa qumica potencial [5]. Es lquido de origen natural
compuesto por una mezcla de hidrocarburos. La formacin del petrleo se debe al desarrollo de rocas sedimentarias depositadas en ambientes marinos o prximos al mar,
y es el resultado de un proceso de descomposicin de organismos de origen vegetal y
animal, que en tiempos remotos quedaron incorporados en esos depsitos (Fig. 2.1 1 )
[6].
1
Imagen tomada de sitio http://www.telefonica.net/web2/marcillaot/
10
2.1. EL PETRLEO EN AGUAS PROFUNDAS
Figura 2.1: Proceso de formacin del petrleo

El nombre de petrleo proviene del trmino griego que significa Aceite de Roca
(Petro: Roca y Oleum: aceite) y se utiliz por primera vez en el tratado De Natura
Fossilium por el mineralogista Geor Bauer (conocido tambin como Georgius Agricola)[7].
No se sabe exactamente cundo fue utilizado por primera vez el petrleo por la
humanidad, sin embargo, las primeras civilizaciones supieron reconocer y aprovechar
la utilidad de ste energtico. Se tiene conocimiento que se usaba el asfalto (una forma
muy viscosa de petrleo) para la construccin de barcos en los aos 6000 a. de C.
y en la construccin de las pirmides de Egipto alrededor del ao 3000 a. de C. [8].
Inclusive, en algunos relatos bblicos se hace mencin al uso de brea o asfalto durante
la construccin de la torre de babel y por No para la construccin del arca.
En un principio las nicas fuentes de petrleo que se conocan eran los pozos de
alquitrn (Fig. 2.2 2 ), sin embargo estos no eran muy productivos. Debido al aumento de la demanda del petrleo, se tuvieron que buscar nuevas fuentes proveedoras
de dicho combustible. En el ao de 1859 el coronel Edwin L. Drake realiz el descubrimiento y la perforacon del primer pozo para extraer petrleo cerca de Titusville,
Pensilvania[5]; siendo el inicio de la industria petrolera en el mundo.
Figura 2.2: Pozo de alquitrn

A partir de los descubrimientos de la dcada de 1850, se comenzaron a construir numerosas refineras para convertir el petrleo en queroseno y aceites lubricantes que se
2
Imagen tomada de http://blog.webosaurs.com/2009/10/page/2/
CAPTULO 2. ANTECEDENTES
11
utilizaban en lmparas y en las mquinas de la revolucin industrial. No obstante, la

importancia del petrleo para la humanidad comenz cuando esta fuente de energa
reemplaz al carbn como principal combustible para las mquinas que se utilizaban
durante esa poca, esto debido principalmente a que el petrleo es una fuente de energa ms verstil y potente que el carbn.
Con la aparicin de los primeros automviles, gracias a la invensin del motor de
gasolina (1885) y el motor de diesel (1892) [8], aument la demanda de este tipo de
combustibles y propici que esta industria fuera avanzando cada vez ms, hasta convertirse en un elemento primordial para el desarrollo econmico e industrial.
Actualmente, el petrleo es la fuente de energa ms importante y se ha convertido
en una de las materias primas de mayor importancia para la humanidad. Es el recurso
que proporciona el mayor porcentaje del total de la energa que se consume en todo
el mundo. Alrededor del 40% del total del suministro de energa a nivel mundial es
suministrada por el petrleo [9]. Por su parte el transporte depende de un 90% del
petrleo (ya sea transporte terrestre, martimo o areo). Histricamente el consumo
mundial de petrleo ha seguido muy de cerca la trayectoria del crecimiento del producto interno bruto de muchos pases. Por lo que el crecimiento econmico de los
pases industrializados va acompaado por un consumo creciente de petrleo. Se estima que la demanda de petrleo seguir creciendo ms rpidamente en los pases con
mayor desarrollo. Por lo tanto, la economa y el estilo de vida de las sociedades industrializadas, dependen en gran medida del suministro suficiente de petrleo que se
tenga as como de su precio [10].
Aunque el petrleo es la principal fuente de energa, tambin se utiliza como materia prima en la fabricacin una gran variedad de productos, entre los que encontramos:
plsticos, pinturas, fertilizantes, detergentes, cosmticos, insecticidas e inclusive algunas grasas y complementos alimenticios (Fig. 2.3 3 ).
Figura 2.3: Principales usos del petrleo

3
Imagen tomada del sitio http://atinachile.bligoo.com/content/view/22576
12
Debido a que este energtico es un recurso no renovable y a la dificultad de sustituir

el petrelo por otras fuentes fsiles de energa, algunos pases han buscado alternativas para tener una mayor produccin de petrleo. Se ha descubierto que las mayores
reservas petrolferas se encuentran en zonas conocidas como aguas profundas, hasta las
que hay que descender para poder extraer el petrleo.[11]
2.1.2
El petrleo en Mxico
En Mxico se tiene conocimiento del petrleo desde la poca prehispnica en la que los
indgenas le daban diversos usos a los depsitos de chapopotli en los cuales el petrleo acumulado era utilizado con fines medicinales y para la elaboracin de figurillas
de arcilla, impermeabilizante, material de construccin, etc. [12].
A pesar de los abundantes depsitos de betn negro (conocidos como chapopoteras)
en la zona del Golfo de Mxico, durante el periodo de la poca colonial los espaoles
no le dieron un uso prctico ya que cuando se presentaba en la superficie de la tierra
era visto como un elemento nocivo para las actividades agrcolas y ganaderas [13].
Por otra parte, la historia comercial comienza hasta 1863 cuando el cura Manuel Gil
y Senz descubri un pozo superficial llamado Mina de Petrleo de San Fernando, en
el estado de Tabasco. Durante esa poca el Emperador Maximiliano intento promover
las actividades petroleras mediante el otorgamiento de concesiones para la explotacin
de depsitos naturales.
Sin embargo, la etapa de produccin de petrleo a gran escala en Mxico, comienza
a partir de 1901 con el inicio de las operaciones de la empresa Mexican Petroleum [12].
En se mismo ao el Ing. Ezequiel Ordez descubre un yacimiento petrolero en el
campo de el bano en San Luis Potos que fue llamado El Pez. De esta manera inicia
de forma organizada la exploracin petrolera en Mxico [14]. En un principio no se
tuvo mucho xito en la perforacin de varios pozos. Fue hasta el ao de 1904 cuando
la empresa Mexican Petroleum tuvo el primer gran pozo productivo.
Para 1920 haba alrededor de 80 compaias productoras de petrleo en el pas y 17
compaias exportadoras, en las cuales el mayor porcentaje de su capital era angloamericano. Para fomentar la inversin de capital nacional en la industria del petrleo,
el gobierno cre Petrleos de Mxico. El 18 de marzo, el presidente Lzaro Crdenas
decreta la disponibilidad de Mxico para indemnizar a las compaias petroleras sus
inversiones. De esta manera se crea Petrleos Mexicanos (PEMEX) como organismo
encargado de la adminsitracin y explotacin de los hidrocarburos [1].
El petrleo ha sido un factor determinante en Mxico ya que ha influido directamente en el desarrollo de su economa, las finanzas pblicas, la industria, el desarrollo
tecnolgico, la balanza comercial y en las relaciones exteriores que mantiene con otros
pases, principalmente con Estados Unidos [14].
Actualmente la actividad de exploracin en Mxico se ha localizado en dos regiones
[15]:
1. Exploracin en tierra: Esta exploracin se realizaba a principios del siglo pasado
en lugares conocidos como chapopoteras.
13
2. Exploracin en aguas someras: Esta etapa de exploracin se hizo a finales del

siglo pasado en la regin de Campeche, y en las costas de Veracruz y Tamaulipas.
Actualmente PEMEX ha descubierto la oportunidad de realizar la exploracin y explotacin de yacimientos de petrleo en aguas profundas. Las aguas profundas representan una gran oportunidad para el pas, principalmente por el volumen de petrleo
que se estima podra existir.
En el golfo de Mxico se han descubierto la existencia de estructuras geolgicas capaces de contener hidrocarburos. Los trabajos de exploracin realizados por Petrleos
Mexicanos han permitido obtener informacin geolgica y geofsica diversa para estimar el potencial e identificar las principales cuencas petroleras que pueden existir en
esa regin.
2.1.3
Exploracin y explotacin de petrleo en aguas profundas
El concepto de ingeniera offshore o mar adentro, surge en el siglo pasado (en la dcada
de lo 60s) en el continente europeo. Sin embargo, por los elevados costos que se requeran para generarla, no tuvo un gran desarrollo. En la dcada de los 70s debido al
conflicto blico rabe-israel, queda al descubierto la importancia que tiene el mercado
del petrleo en el desarrollo econmico en el mundo [16]. Esta situacin provoc el aumento de precio del petrleo a nivel mundial, con lo cual se pudieron sufragar los altos
costos de la tecnologa requerida para explorar y explotar nuevos yacimientos ubicados en el fondo del mar. De esta forma la ingeniera offshore se desarroll exitosamente,
principalmente por la participacin y las capacidades de los centros de investigacin y
la asociacin entre universidades y empresas.
A mediados de la dcada de los 80s, la actividad de exploracin de campos de desarrollo se desplaz a las aguas profundas [17]. Principalmente porque en esta poca
se descubrieron grandes cantidades de petrleo en el Golfo de Mxico; lo que gener
la necesidad de nuevas estrategias de exploracin y aplicar nuevas tecnologas e infraestructura ms especializada.
Segn el Minerals Management Service (MMS) de Estados Unidos, las aguas profundas significan distancias mayores a 305 mts. (aproximadamente 1000 pies) hacia el
fondo del mar[18]. En nuestro pas, la exploracin y explotacin de petrleo en aguas
profundas se refiere a regiones ubicadas en tirantes de agua mayores a 500 metros (distancia entre la superficie y el lecho marino, figura 2.4 4 )[1].
Imagen tomada del sitio http://www.pemex.com [1]
14
Figura 2.4: Regiones de aguas profundas

Hoy en da, el desarrollo de yacimientos en aguas profundas est dejando atrs el
desarrollo de yacimientos en tierra y en aguas poco profundas como se vena haciendo
en muchos pases hace algunos aos. Las principales razones de esto son [18]:
Las fuentes de petrleo en tierra o en aguas someras son muy limitadas.
Las reservas que se tienen en aguas profundas son mucho ms productivas (Aproximadamente 20 veces ms en petrleo y 8 veces ms en gas)
Es necesario realizar una mayor inversin en investigacin, pero se tienen ms
retribuciones.
Se mejora y se desarrolla nueva tecnologa.
Debido al potencial de petrleo que se estima puede haber en aguas profundas, los
pases se han dado a la tarea de buscar nuevas fuentes proveedoras en esos campos.
La tabla 2.1 muestra los principales pases productores de petrleo en aguas profundas
[1].
Principales Pases
1.- Brasil
2.- Nigeria
3.- Estados Unidos
4.- Angola
5.- Egipto
6.- Gran Bretaa
7.- Noruega
8.- India
9.- Filipinas
10.- Guinea Ecuatorial
Tabla 2.1: Principales pases productores de petrleo en aguas profundas
15
En Mxico las aguas profundas se encuentran en una extensin de 575 mil kilmetros cuadrados ubicada en las aguas territoriales del Golfo de Mxico (fig 2.5 5 ). En
esta zona, PEMEX ha identificado la existencia de estructuras geolgicas que pudiesen
contener hidrocarburos. Esto representa una gran oportunidad para el pas, debido
a la gran cantidad de petrleo que puede existir en esta zona de aguas profundas.
Sin embargo, algunos de estos yacimietnos son transfronterizos, es decir, que geogrficamente se localizan en las fronteras de dos o ms pases. Varios de estos campos
se encuentran localizados en la frontera con Estados Unidos en donde ya se tiene un
tiempo de exploracin de 15 aos y que en los ltimos 10 aos perforan alrededor de
100 pozos por ao [1]. Por esta razn, es urgente tomar medidas para la exploracin
de dichos campos ya que de no hacerlo se producir un efecto popote y se veran
afectados los campos que estn de lado mexicano.
Figura 2.5: Superficie de las Aguas Profundas en el Golfo de Mxico
Para el MMS de Estados Unidos, el Golfo de Mxico se ha convertido en su mayor

proveedor de gas y petrleo. El desarrollo de recursos ha convertido al Golfo de
Mxico en el principal punto de exploracin y produccin de gas y petrleo a nivel
mundial. Con la obtencin de petrleo en aguas profundas del Golfo de Mxico, el
gobierno americano espera que la produccin se duplique para la prxima dcada, es
decir, que para el 2011 se obtenga casi el 80% de la produccin total [19].
Para el gobierno de los Estados Unidos, las aguas profundas representan un papel
fundamental en el futuro de la energa de esa nacin. Actualmente cuentan con un gran
inventario de contratos de arrendamiento disponibles para la industria de exploracin.
Esto genera que la tecnologa siga avanzando, a fin de facilitar las actividades en las
situaciones extremas de las aguas profundas. Todos estos factores, aseguran que las
aguas profundas del Golfo de Mxico sigan siendo uno de los principales campos de
desarrollo en el mundo [20].
5
16
2.1.4
Retos y dificultades en la exploracin y explotacin en aguas

profundas
El reto de explorar y explotar yacimientos de petrleo en aguas profundas no es sencillo, ya que es totalmente diferente a realizarlo en tierra o cerca de las costas como
actualmente se hace en Mxico. Estudios realizados por PEMEX han determinado que
el tiempo que llevara realizar un proyecto en aguas profundas es de aproximadamente
8 a 10 aos, ya que se tiene que pasar por diversos procesos (Fig. 2.6 6 ) [1].
Figura 2.6: Procesos para desarrollar un proyecto en Aguas Profundas

Cuando se explotan yacimientos en aguas profundas las empresas se enfrentan a
riesgos geolgicos, de operacin y financieros. Estos riegos producen grandes retos y
dificultades que las empresas tienen que solucionar [1]:
1. Corrientes marinas: Estas fuerzas externas provocan el movimiento de las estructuras que producen vibraciones en las tuberas y desgastan los equipos de
perforacin.
2. Cambios de temperatura: Se presentan diferentes cambios de temperatura entre
la superficie, el lecho marino y las zonas perforadas; que pueden complicar el
bombeo del fluido. Adems, las bajas temperaturas que se pueden presentar
implican el uso de diversas tcnicas y metodologas para que los fluidos puedan
circular por las tuberas [17].
3. Problemas de perforacin: Cuando se atraviesan formaciones someras, se pueden
presentar flujos de agua de alta presin, flujos de gas y presiones anormales.
4. Instalaciones submarina: Despus de la perforacin, el pozo debe ser controlado
por un sistema de produccin. As que se instala una estructura submarina en el
fondo mar para su control. Esto requiere la aplicacin de nuevas tecnologas, as
como un alto grado de automatizacin.
5. Operacin remota: Se requiere de la implemetacin de tecnologa robtica que
permita operar la instalacin submarina desde el fondo del mar.
6
17
6. Altos costos: Los campos son desarrollados con un menor nmero de pozos que
en aguas someras, debido a que la inversin que se requiere en cada uno de ellos
es muy elevada.
7. Geometra de yacimientos: En aguas profundas los yacimientos prodran ser
diferente a los de aguas someras, lo que dificulta su exploracin y explotacin.
8. Infraestructura: Es necesario renovar toda la infraestructura con la que se cuenta
actualmente y utilizar nuevas metodologas y tecnologas para llevar a cabo el
desarrollo del proyecto.
Para poder trabajar bajo estas condiciones y enfrentar los retos que se presentan, es
necesario utilizar nuevas tecnologas, las cuales a pases como Estados Unidos, Brasil
y Noruega les ha tomado alrededor de 15 aos desarrollar. Adicionalmente tambin
se requiere contar con los conocimientos y las habilidades necesarias para poder seleccionar correctamente las tecnologas, utilizarlas adecuadamente y resolver cualquier
contratiempo que se pueda presentar.
2.2
Desarrollo de campos petroleros en alta mar
Como se ha mencionado anteriormente, el desarrollo de campos petrolferos es un

proceso complicado y muy costoso que requiere de la implementacin de tecnologa
de vanguardia. Actualmente el desarrollo de yacimientos en alta mar esta formado por
cuatro elementos principales [18](Fig. 2.7 7 )
1. Sistema Submarino.
2. Lneas de flujo, tuberas y riser.
3. Plataformas fijas o flotantes.
4. Sistema de procesamiento superior.
2.2.1
Sistema submarino
El sistema submarino es necesario para controlar el sistema de produccin una vez que
se ha completado el proceso de perforacin. En campos desarrollados en tierra o en
aguas someras se implementa un sistema llamado dry completion (finalizacin seca),
que va desde la cabeza del pozo y se extiende hasta la parte superior de la plataforma.
Sin embargo, en aguas profundas no es tan sencillo instalar este tipo de sistemas, por
lo que se instala una estructura submarina que controla la cabeza del pozo en el fondo
marino llamda wet completion (Finalizacin humeda) [17].
El sistema submarino se puede dividir en tres partes principales [18]:
7
Imagen tomada del libro Introduction to Offshore Pipelines and Risers [18]
18
2.2. DESARROLLO DE CAMPOS PETROLEROS EN ALTA MAR
Figura 2.7: Elementos para el desarrollo de campos en alta mar

1. Estructura de la cabeza del pozo (rbol de navidad): Es una estructura de la
parte superior del armazn de perforacin, que se utiliza para montar un panel
de control de vlvulas, la cual tiene una forma de pino.
2. Sistema de control: Es un sistema para el control de bobinas, vlvulas y monitorear la presin, temperatura y la posicin de los sensores que se montan en
sistema submarino.
3. Sistema de conexin: Se utiliza para realizar puentes y terminacin de tuberas.
Los sistemas submarinos se clasifican por el nmero de pozos y los mtodos de
conexin de la siguiente manera [17] (Fig. 2.8 8 ):
Single Well Tieback: Es un pozo submarino que est conectado a la plataforma
flotante de produccin por una lnea de flujo simple o dual.
Daisy Chain Tieback: Pozos submarinos individuales que son unidos por una
lnea de flujo simple que forman un ciclo desde el pozo principal a los dems
pozos.
Cluster Well Manifold: Un colector es empleado en varios pozos submarinos individuales agrupados, generalmente en un centro de perforacin comn.
Multi-well Template: Proporciona una estructura que facilita la operacin de perforacin desde un solo punto. Los puntos de produccin se alcanzan mediante
perforacin direccional y todo el rbol de pozos est localizado sobre la estructura lo que permite un fcil acceso.
8
Imagen tomada del artculo Design and Installation of Deepwater Petroleum Pipelines [17]
19
Figura 2.8: Clasificacin de sistemas submarinos
2.2.2
Lneas de flujo, tuberas y sistemas riser
Las lneas de flujo tambin se denominan lnea de produccin o de importacin y se

encargan de transportar fluidos sin procesar (petrleo crudo o gas). Generalmente la
transmisn de fluido puede llevar algunos slidos como parafina y arena (Fig. 2.9 9 )
[18].
Por su parte, el sistema de tuberas transporta petrleo y gas procesado, por lo que
primero se realiza un proceso para separar el gas, el agua y algunos slidos que vienen
mezclados con el petrleo al momento de su extraccin. En general el tamao de las
tuberas es mayor a las lneas de flujo, adems presentan bajas temperaturas y baja
presin. [18] (Fig. 2.9).
Un sistema riser generalmente es un tubo el cual conecta la tubera submarina al sistema de la plataforma flotante. En el desarrollo del campo de aguas profundas el riser
puede estar hecho por una tubera rgida o flexible. El riser permite el acceso desde el
fondo del mar a las instalaciones de la plataforma para el desarrollo de estaciones submarinas, tuberas para importar y exportar productos, y una plataforma de servicios
(lneas elctricas y de control) [21] (Fig. 2.9).
9
20
2.2. DESARROLLO DE CAMPOS PETROLEROS EN ALTA MAR
Figura 2.9: Lnea de flujo, tuberas y sistemas riser
A continuacin se describen los sistemas de riser que se tienen actualmente para

aguas profundas, as como las caractersticas de cada uno de ellos [17] (Fig. 2.10 10 ):
Flexible Riser: Puede tener gran longitud, es robusto con respecto a la dinmica y
tiene un bajo costo de instalacin.
Steel Catenary Riser (SCR): Es un riser de acero que pude tener grandes dimetros
y resiste altas temperaturas.
Top Tension Riser (TTR): Riser de alta tensin que puede contener largas trayectorias.
Tension Leg Riser (TLR): Este riser es la combinacin de SCR y un sistema de riser
flexible. Es considerado como una de las mejores opciones en cuanto al costo
efectivo. Este tipo de sistemas pueden tener un buen aislamiento y menor carga
en los buques.
10
Imagen tomada del artculo Design and Installation of Deepwater Petroleum Pipelines [17]
21
Figura 2.10: Tipos de sistema de riser para aguas profundas

La seleccin del sistema de riser adecuado depende de diversos factores como: la
profundidad del agua, el medio ambiente, el mantenimiento de la estacin, la produccin de fluidos, la instalacin del sistema, entre otros.
2.2.3
Plataformas fijas y flotantes
Para realizar el proceso de extraccin de petrleo, es necesario contar con una estructura en la superficie del agua que se encargue de recibir los fluidos que se obtengan
del yacimiento. La seleccin de este tipo de superficies depende de la profundidad del
agua, el ambiente, el mantenimiento de la estacin, el proceso de fluido, el equipo de
procesamiento, ect. [18]
Actualmente existen dos tipos de plataformas (Fig. 2.11 11 ):
11
22
2.3. MODELOS OCEANOGRFICOS

Plataformas Fijas: Este tipo de estructuras estn fijas en el fondo del mar y se
pueden utilizar cuando la profundidad es poca.
Plataformas Flotantes: Este tipo de estructuras se encuentran sujetas por cables
o cadenas en el fondo del mar y son utilizadas cuando la profundidad es muy
grande.
Figura 2.11: Tipos de plataformas que se utilizan en los campos petroleros

Los sistemas de produccin flotantes son considerados como la solucin ms practica y comercialmente ms atractiva para el desarrollo de campos en aguas profundas
[21].
2.2.4
Sistema de procesamiento superior
Antes de enviar el crudo a las refineras, es tratado por las instalaciones que se encuentran en las plataformas que se instalan sobre la superficie marina (fijas y flotantes). El
diseo de este tipo de instalaciones tiene que ser compacta debido al lmite de espacio
y peso que se tiene en las plataformas [18].
Los requisitos del diseo de los sistemas de preprocesamiento superior, dependern
directamente de las condiciones en las que se est desarrollando el pozo, y los planes
para su desarrollo futuro.
2.3
Modelos oceanogrficos
Los modelos oceanogrficos permiten entender la estructura de los movimientos del

oceano para ser capaces de predecir su comportamiento en un estado futuro. A travs
de ciencias bsicas como la fsica, la qumica, la biologa, la geologa y las matemticas;
la oceanografa se encarga del estudio de los procesos que ocurren en los ocanos y se
23
divide en cuatro ramas: biolgica, fsica, qumica y geolgica. Actualmente la investigacin oceanogrfica ha incorporado diversos modelos numricos para describir la
hidrografa de las masas de agua y el comportamiento del movimiento del ocano.
Para el desarrollo de modelos oceanogrficos se utiliza la oceanografa fsica, la cual
se encarga de estudiar la distribucin de propiedades del agua tales como: la temperatura, salinidad, densidad, oxgeno, ect. Adems, se encarga de estudiar el movimiento del ocano en respuesta a las fuerzas que actan sobre l (corrientes, mareas
y oleaje). La oceanografa fsica est dividida a su vez en dos subramas: oceanografa
descriptiva y oceanografa dinmica.
La obtencin de lo datos oceanogrficos para la construccin de los modelos, provienen de diferentes fuentes como: buques oceanogrficos, aviones, satlites y sensores posicionados estratgicamente para recolectar informacin y analizarla posteriormente. Las propiedades del agua ms importantes para el estudio de dinmica de
fluidos son la presin, la temperatura y la salinidad.
El estudio del ocano es un fenmeno fsico complejo que involucra el espacio y
el tiempo, que manifiesta caractersticas de mesoescala como giros y frentes. El estudio de la evolucin, deformacin e interaccin de estas caractersticas dinmicas y su
modelado es la esencia de la investigacin en la circulacin ocenica.
2.3.1
Base de datos oceanogrfica
Este tipo de base de datos se encuentra conformada con la informacin recabada durante el desarrollo de campaas oceanogrfica. La informacin se obtiene a travs de
las mediciones realizadas por las diversas estaciones hidrogrficas ubicadas en el Golfo
de Mxico, principalmente con dos medios:
1. A travs de una sonda hidrogrfica CTD (Conductivity, Temperature y Depth) que
registra datos continuos fsicos y qumicos desde la superficie del mar hasta el
suelo del ocano en diferentes posiciones geogrficas del Golfo de Mxico.
2. Recolectando muestras de agua de mar a difrentes profundidades, las cuales son
analizadas posteriormente a bordo del buque oceanogrfico para determinar sus
concentraciones de oxgeno, salinidad y micronutrientes disueltos. Los resultados obtenidos se utilizan para calibrar los datos del CTD.
Todos los datos crudos obtenidos y toda la informacin que resulta del procesamiento
numrico y grfico que se aplica de los datos originales, conforman una base de datos
oceanogrfica, que por sus caractersticas en nuestro pas seria de gran importancia y
utilidad.
2.3.2
Mtodo geostrfico
El mtodo geostrfico es uno de los ms utilizados en oceanografa dinmica para el

clculo de corrientes y flujos ocenicos. En este mtodo los efectos de friccin, junto
24
2.3. MODELOS OCEANOGRFICOS
con la fuerza de gravedad del sol y la luna son despreciables, por lo que el balance se
da solamente entre las fuerzas horizontales de presin y la fuerza de Coriolis.
El mtodo geostrfico requiere la utilizacin de tareas de procesamiento cientfico
para realizar el clculo de las propiedades que se describen a continuacin.
Topografa dinmica: Para el clculo de la altura dinmica de cada estacin se utilizan
datos de conductividad, temperatura y presin registrados pro el CTD, tomando como
niveles de referencia de movimiento cero profundidades de 100, 250, 500, 1000, 1500,
2000 y 3000 mts.
Campo de velocidad geostrfico: El clculo de velocidad geostrfico de planta, relativo
a superficies de movimiento cero, se calcula con base en las anomalas geopotenciales
corresponientes a cada estacin hidrogrfica. A partir de esta distribucin irregular de
alturas dinmicas, se calculan nuevos valores mediante un algoritmo de interpolacin
bidimensional. Al realizar esta interpolacin dentro de una malla regular definida,
se obtiene una distribucin uniforme de alturas dinmicas dentro del rea de estudio.
Con base a esta distribucin uniforme de anomalas geopotenciales, para cada nodo de
la malla se calculan los componentes norte-sur y este-oeste de la velocidad geostrfica
utilizando la ecuacin 2.1:
V =
Da Db
2wsen
(2.1)
donde V es la velocidad geostrfica, Da Db es la distancia geopotencial, w la velocidad

angular de la tierra y la latitud media entre nodos adyacentes.
Secciones verticales de velocidad: Se obtienen de las alturas dinmicas de cada uno de
los trasectos latitudinales y meridionales. Para cada una de las estaciones hidrogrficas se calculan los perfiles de velocidad de la superficie a diferentes profundidades
utilizando la ecuacin 2.1. Posteriormente se ajustan los datos mediante un producto
tensorial de interpolacin bidimensional para graficarlos y representar el campo de
velocidad de los transectos latitudinales y meridionales, en donde las velocidades positivas tienen direccin hacia el norte y el este, y las velocidades negativas hacia el sur
y el oeste
Transporte geostrfico horizontal para todos los transectos (latitudinales y meridionales),
se estima a partir de las velocidades geostrficas calculadas. El transporte geostrfico
se define en la ecuacin 2.2
Q=V A
(2.2)
donde V es la velocidad y A representa el rea perpendicular o normal al flujo.

Para determinar el transporte, se toma la malla definida para las interpolaciones
de velocidad y se calcula el transporte para cada rectngulo de la misma mediante la
ecuacin 2.3
X X vi,j + vi,j+1 + vi+1,j + vi+1,j+1
xy
(2.3)
Tl =
4
i
j
25
donde Tl es el transporte horizontal del transecto 1, vi,j es la velocidad horizontal del

nodo (i, j), x es la longitud horizontal de un rectngulo de la malla y y la longitud
vertical de un rectngulo de la malla.
Transporte vertical: las velocidades verticales se estiman del balance entre la divergencia horizontal y el gradiente vertical. La distribucin se calcula utillizando superficies
separadas por intervalos de profundiad. Cada una de las superficies se representa en
una malla compuesta de rectngulos. En cada uno de los nodos se hace una interpolacin de la velocidad vertical y se utilizan los valores de velocidad de los cuatro
vertices para obtener una velocidad vertical promedio, la cual se multiplica por el rea
del rectngulo para obtener el transporte vertical asociado al rectngulo. Finalmente,
se calcula el transporte total para cada superficie sumando los transportes de todos los
rectangulos.
El trasnporte positivo indica flujo hacia arriba y el negativo flujo hacia

abajo.
2.4
Simulacin
La simulacin se puede ver como una tcnica que nos permite construir un modelo
representativo de la realidad, con el cual se pueden realizar experimentos para analizar
su comportamiento en diversas situaciones. Por su parte, la simulacin por computadora se pude definir como una tcnica numrica para conducir experimentos en una
computadora digital, los cuales requieren ciertos tipos de modelos lgicos y matemticos, que describen el comportamiento de un sistema (o algn componente de este) en
periodos extensos de tiempo real [22].
Una simulacin tiene los siguientes propsitos:
Conocer el funcionamiento de un sistema. La observacin detallada de lo que
se est simulando permite un mejor conocimiento, proporcionando informacin
para mejorarlo.
Ayuda a identificar problemas.
Permite experimentar con nuevas situaciones y generar polticas de decisin.
Experimentar en un tiempo reducido. La simulacin permite, en algunos segundos, analizar procesos que llevara meses o aos desarrollarlo en tiempo real.
Reduccin de los procesos analticos de investigacin y desarrollo convencional.
Permite estudiar situaciones complejas.
Para educacin y entrenamiento.
Para tener una visin global del sistema.
26
2.4. SIMULACIN
Por su parte, la simulacin por computadora requiere de un modelo que pueda ser
desarrollado por medio de un programa en computadora. La simulacin por computadora se encuentra dentro del campo de la Computacin Cientfica, y nos ofrece la
oportunidad de estudiar y entender el comportamiento de fenmenos naturales [4].
En la actualidad se han desarrollado un gran nmero de simuladores, los cuales
tienen diversos objetivos (simuladores de vuelo, juegos militares, simuladores de economa,
fenmenos fsicos, ect.). El fundamento principal para usar herramientas de simulacin es adquirir conocimientos que permitan la prediccin del futuro.
2.4.1
El proceso de simulacin por computadoras
Para llevar a cabo el desarrollo de experimetos de simulacin en computadoras, es

necesario utilizar diversas tcnicas como: estadstica matemtica, anlisis numrico,
programacin, ect. El proceso de experimentos de simulacin en computadora est
formado por las siguientes etapas [22]:
1. Formulacin del problema: Al igual que en otras reas de la computacin cientfica,
se tiene que comenzar con la formulacin de un problema y definir claramente
los objetivos de la investigacin que se va a desarrollar.
2. Recoleccin y procesamiento de datos: Esta etapa consiste en captar de la realidad,
los hechos y la informacin disponible que pueda ser procesada posteriormente.
En muchas ocaciones, es necesario realizar una conversin de la informacin
obtenida, a datos especficos que puedan ser interpretados y procesados de manera ms eficiente.
3. Formulacin del modelo matemtico: Esta etapa consiste en construir un modelo
para representar el fenmeno que se va a simular. El modelo nos va a proporcionar descripciones o predicciones sobre comportamiento del sistema dado.
4. Estimacin de parmetros a partir de los datos reales: Una vez que se ha formado
el modelo matemtico, es necesario estimar los valores de los parmetros del
modelo y posteriormente probar su significado estadstico.
5. Evaluacin del modelo: En este paso se va a probar la exactitud de nuestro modelo
para poder tener un juicio inicial de nuestro sistema.
6. Desarrollo del programa para computadora: Para llevar a cabo esta etapa, se debe de
representar el modelo mediante un algoritmo y posteriormente elegir el lenguaje
en el cual se va a desarrollar la aplicacin. El desarrollo del programa generalmente debe permitir especificar las condiciones iniciales del sistema, procesar los
datos y generar reportes de salida grficos para presentarlos a los usuario.
7. Validacin: Para poder validar la aplicacin, es necesario hacer un proceso de
comparacin de resultados. El problema de validar modelos de simulacin es
27
uno de los procesos ms complicados ya que implica diferentes complejidades

(terica, prctica, estadstica, etc.).
8. Diseo de experimentos: Despus de estar satisfechos con la validacin de nuestro modelo y de la aplicacin desarrollada, podemos dirigir de manera eficaz los
experimentos. En esta fase vamos a poder seleccionar las combinaciones de factores que intervienen en el proceso (condiciones iniciales, datos aleatorios, etc.) y
seleccionar el orden en el que se deben llevar a cabo las simulaciones.
9. Anlisis de los datos simulados: Finalmente se requiere hacer un anlisis de los datos
generados por la computadora y darles una interpretacin.
El proceso de simulacin por computadora nos ayuda a estudiar fenmenos naturales que resultan muy difciles o inclusive imposibles de reproducir en la realidad.
Este medio nos permite resolver conjeturas y tener una prediccin sobre el comportamiento que pudiera tener el sistema en condiciones especiales [4].
2.5
Vibracin
En general los sistemas de ingeniera que poseen masa y elasticidad pueden presentar
movimientos relativos. Si este movimiento se repite despus de un determinado intervalo de tiempo el movimiento se conoce como vibracin [3]. La vibracin se puede ver
como una forma de energa disipada. La vibracin puede presentarse en dos formas:
Vibracin Libre: Es el movimiento cuando el sistema se desplaza de su posicin
de equilibrio. En este movimiento interviene la elasticidad de la estructura, la
friccin y el amortiguamiento.
Vibracin Forzada: Se presenta cuando actan fuerzas externas sobre el sistema
durante su movimiento. Las fuerzas son generalmente de tipo F (t) = F0 sin t o
F (t) = F0 cos t.
En la actualidad una gran variedad de aplicaciones derivadas de estructuras mecnicas, sistemas acsticos, circuitos elctrico, fluidos mecnicos, etc. requieren un anlisis
de vibraciones en sus estructuras. Cuando una estructura es excitada por fuerzas externas en las cuales las frecuencias se encuentran cerca de sus frecuencias naturales, la
vibracin en la estructura se amplifica y el sistema se vuelve inestable. Este fenmeno
se conoce como Resonancia [2]. Dependiendo de los sistemas o de la aplicacin, el
fenmeno de resonancia puede ser deseado o no deseado. En las estructuras mecnicas
las vibraciones pueden provocar desgastes que pueden producir catstrofes. Por otra
parte, en muchas aplicaciones de ingeniera elctrica como al sintonizar una estacin
de radio, el inters esta en encontrar una amplitud de resonancia tan grande como sea
posible.
28
2.5. VIBRACIN
2.5.1
Ejemplos de estructuras afectadas por el fenmeno de vibracin
En la actualidad el fenmeno de vibracin y resonancia es muy estudiado sobre todo

en la construccin de grandes estructuras. Una de las estructuras donde comnmente
se presenta el fenmeno de la vibracin son los puentes, en los cuales si no se hace un
clculo adecuado de la frecuencia natural pueden llegar a colapsarse. Uno de los casos
ms importantes fue el del puente de Tacoma (Fig. 2.12 12 ) que es un puente colgante
de 1600 mts. de longitud. En el ao de 1940 la estructura tuvo un dramtico colapso
debido al viento.
Figura 2.12: Puente Tacoma

Uno de los casos ms recientes en este tipo de estructuras fue el que se di en Junio
del 2000 en el Puente del Milenio, el cual tiene una longitud de 320 mts. y se encuentra
localizado sobre el ro Tmesis en Londres (Fig. 2.13 13 ). Esta estructura en el da de
su inaguracin comenz a tambalearse de forma alarmante debido al paso de cientos
de personas que haban asistido a dicha inaguracin. Dos das despus el puente fue
cerrado y posteriormente fue equipado con amortiguadores para evitar que se tambaleara nuevamente.
Figura 2.13: Puente del Milenio, Londres.

Otro tipo de estructuras que estn expuestas a este tipo de fenmenos de vibracin
son los grandes edificios. En este tipo de estructuras tambin es necesario implementar
12
13
Imagen tomada del sitio http://puentetacoma.blogspot.com/

Imagen tomada del artculo The Quadratic Eigenvalue Problem [2]
29
mtodos numricos complejos para evitar o amortiguar las excitaciones a las que estn
expuestas, principalmente al viento. Uno de los casos ms recientes se dio en el Hotel
Burj Al Arab en Dubai, considerado el hotel ms alto del mundo (Fig. 2.14 14 ). Debido a que sta construccin fue inspirada en una embarcacin de vela y a su cercana
con la costa, al momento de hacer pruebas de viento simulado se present un fenmeno el cual, bajo ciertas condiciones de viento soplando sobre los vrtices afilados
de la estructura, se podan crear minitornados que producan peligrosas vibraciones
que sacudan el edificio y podan llegar a destruirlo. Para solucionar este problema
recurrieron a un peso colgante llamado amortiguador de masa sincronizado que hace
que cuando el viento sople, se mezan 5 toneladas de peso en lugar de la estructura,
amortiguando las vibraciones [23].
Figura 2.14: Hotel Burj Al Arab, Dubai.
Durante el proceso de exploracin y explotacin de petrleo en aguas profundas,

tambin se presenta el fenmeno de vibracin. Este fenmeno se produce principalmente por la fuerza que ejercen las corrientes marinas sobre las estructuras que se
utilizan en la extraccin del petrleo (Fig. 2.15 15 ). Para evitar accidentes debido a las
vibraciones, es necesario tener un modelo que nos permita predecir el comportamiento
de los riser y determinar la resistencia que necesita para soportar la interaccin de las
corrientes.
14
15
Imagen tomada del artculo Burj Al Arab. Una superestructura de 7 estrellas [23]
Imagen tomada del Deepwater Reference Book [21]
30
2.6. TRABAJOS RELACIONADOS
Figura 2.15: Tipos de riser utilizados en los campos de aguas profundas.
2.6
Trabajos relacionados
En la industria petrolera no se cuenta con mucha informacin referente al software

que se utliza para la exploracin y explotacin de petrleo en aguas profundas, ya que
este es desarrollado por empresas privadas las cuales lo venden a costos muy altos
y no proporcionan al pblico mucha informacin respecto a ellos. A continuacin se
describen algunos desarrollos de software utlizados para administrar y monitorear los
sistemas de riser:
Flexcom 3D: Desarrollado por MCS Internacional. Este es un sistema dinmico no
lineal que utiliza el mtodo del elemento finito para anlizar el comportamiento
de una amplia gama de estructuras martimas. Su capacidad de anlisis incluye
la instalacin y la operacin de las estructuras en condiciones extremas. Esta
herramienta es utlizada por las empresas petroleras y fabricantes de equipo para
el anlisis de tuberas flexibles [21] [24].
Visual Orcaflex: Desarrollado por Orcina. Es un software para el anlisis dinmico
de sistemas marinos de alta mar. Esta aplicacin nos ofrece un modelado en
3D, no lineal. Permite realizar un anlisis de fatiga, de tensin, flexin y torsin. Adems permite realizar un anlsis de la maniobra de buques (velocidad
de avance y giro). [21] [25].
Optima: Desarrollado por MCS Internacional. Es un sistema para administracin
de herramientas, planeacin y monitoreo de perforaciones o finalizacin y redireccionamiento de riser [26].
Shear 7: Desarrollado por el MIT. El Shear 7 es una herramienta de modelado para
la prediccin de vibraciones inducida por vrtices (VIV). El programa se encarga
de la prediccin de vibraciones en estructuras martimas que son sometidas a
la variacin de corrientes. El programa permite calcular las respuestas de las
31
estructuras, las medidas de fuerza de arrastre y la prediccin de daos por fatiga

[21] [27].
BOAST: Es una herramienta de simulacin tridimendional trifsica para yacimientos del tipo Black Oil (petrleo negro). Este simulador es de cdigo abierto, se
encuentra programado en FORTRAN y fue desarrollado por el Laboratorio Nacional de Tecnologa Energtica (NETL), el cual se encuentra asociado al Departamento de Energa de los Estados Unidos (DOE). Este simulador puede resultar un poco deficiente, sobretodo a nivel de solucin de sistemas de ecuaciones.
Posteriormente a la versin inicial, sali una versin mejorada BOAST-NFR. Las
mejoras fundamentales de esta versin es que puede ser ejecutada bajo Windows,
los datos pueden ser cargados en hojas de clculo de excel y el cdigo se ejecuta
en subrutinas programadas en Visual Basic.
2.7
Resumen
En este captulo se di un panorama general sobre el proceso y los conocimientos tcnicos que se requieren para la comprensin de este proyecto de investigacin. En la
primera seccin analizamos los antecentes del petrleo y las nuevas fuentes proveedoras que se han descubierto. Adems, se describieron los principales retos y problemas
a los que se enfrenta las empresas durante el proceso de exploracin y explotacin de
petrleo en aguas profundas.
En la segunda seccin se analizaron los cuatro elementos principales que se requieren para el desarrollo de un campo de petrleo en alta mar. Se describireron las
funciones de cada uno de estos elementos, la metodologa y tcnicas requeridas para
su desarrollo e implementacin.
En la tercer seccin se describe en que consiste un modelo oceanogrfico y la importancia que tiene en el anlisis del comportamiento del movimiento de los ocanos.
En la cuarta seccin describimos el concepto de simulacin y su aplicacin en el
desarrollo de aplicaciones computacionales para el estudio de distintos fenmenos.
Se mencionaron las ventajas de la implementacin de simuladores. Y finalmente se
analizaron las partes principales y el proceso general para la construccin de un simulador.
En la quinta seccin del captulo se vi el fenmeno de vibracin y las dos formas en
las que se puede presentar. Se analiz el concepto de resonancia y la forma en como
afecta a las estructras mecnicas. Tambin se vieron algunos ejemplos de estructuras
donde se presentan estos fenmenos y la forma en como se solucionaron.
Finalmente en la ltima seccin se presenta una investigacin sobre simuladores que
existen actualmente en la industria petrolera y el diseo de algunos softwares especficos para los sistemas riser.
32
2.7. RESUMEN
Captulo 3
Modelo del sistema
Todo aquello que un hombre puede
imaginar, otros lo podran realizar.
Julio Verne
l modelo de un sistema, es una estructura que permite conocer o deducir su comportamiento a partir de bases tericas, y que tiene una gran
semejanza al sistema real que se est analizando [28].
Como se coment en el Captulo 1, la eleccin de un modelo adecuado
para representar el sistema, es uno de los pasos ms importantes en el proceso de la simulacin. La eleccin del sistema debe estar en funcin de los
objetivos que se platearon en el inicio de la investigacin y de la presicin
que se requiera, para que los resultados obtenidos que se obtengan estn
apegados a la realidad.
Para formular un modelo, es necesario realizar un proceso de observacin del fenmeno que se est estudiando que nos permita hacer una
representacin fsica del sistema y formular una o ms hiptesis relativas
a su funcionamiento. La abstraccin de esta representacin fsica nos va
a permitir la formulacin de un modelo matemtico que describa su comportamiento. Adems nos va a permitir desarrollar modelos a escala para
poderlo someter a pruebas.
El modelo matemtico nos va a ofrecer un soporte confiable para nuestra
simulacin [4], ya que la validacin de los resultados obtenidos dependen
directamente del modelo matemtico que se formul.
Una de las consideraciones ms importantes que se deben de tomar en
cuenta para la formulacin de un modelo matemtico es su complejidad.
De esto depende la dificultad de programacin y el tiempo de cmputo que
requieran para resolverlos [22].
33
34
3.1. MODELO FSICO DEL SISTEMA

En este captulo describiremos la representacin fsica de nuestro sistema y los pasos que se realizaron para la construccin del modelo matemtico
que nos va a permitir predecir su comportamiento.
3.1
Modelo fsico del sistema
Nuestro modelo fsico representa un riser, que es una tubera de metal que se encuentra sumergida en el ocano y que est conectada a la plataforma flotante y una estacin
submarina (ver seccin 2.2.2). Sobre esta tubera, las corrientes marinas generan vibraciones que producen fatigas y desgastes sobre la estructura. El propsito de nuestro
modelo es predecir el comportamiento de estas estructuras en condiciones extremas,
producidas principalmente por las corrientes marinas.
Por ello, el sistema se puede representar como una tubera de metal sumergida en
agua, sobre la cual se presentan fuerzas externas que producen vibraciones en la estructura metlica. Para analizar las vibraciones a lo largo de la estructura, el riser se
dividi en varios puntos (partculas), sobre los cuales se analiza la fuerza que ejercen
las corrientes marinas y se determina que tanto afectan estas fuerzas a la estructura en
un punto determinado. Estas corrientes marinas se pueden presentar con una fuerza
y/o direcccin diferente en cada una de las partculas (Fig 3.1). Consideramos que la
divisin de las partculas de masa en la que se dividir el sistema ser homognea, de
tal forma que cada uno de los puntos se encontrarn equidistantes uno de otro.
Figura 3.1: Representacin del modelo fsico. Tubera sumergida en agua dividida en n
partculas
Este sistema se puede representar como una cadena de vibracin monoatmica unidimensional que esta formado por un conjunto de partculas (con masas iguales m)
que son equidistantes, y unidas por medio de resortes (con una constante elstica k)
(Fig. 3.2). Cuando se hace vibrar una cadena conformada de esta forma, los diversos
movimientos que pueda tener respecto a un sistema de coordenadas cartesianas (en
CAPTULO 3. MODELO DEL SISTEMA
35
este caso cada una de nuestras partculas), definen los grados de libertad de nuestro
sistema [28].
Figura 3.2: Movimiento del riser provocado por las corrientes marinas
De esta manera obtuvimos un sistema con n grados de libertad y en consecuencia
tendremos n modos normales de vibracin, en el cual las masas son iguales para todas
sus partculas m0 = m1 = m2 = mn1 = mn = m y de igual forma los resortes
tambin tienen el mismo valor de elasticidad k0 = k1 = k2 = kn1 = kn = k.
En cada uno de los modos normales de vibracin, todas las partculas de la cadena
oscilarn con la misma frecuencia y fase. Por lo que cualquier vibracin en nuestro
sistema ser el resultado de una superposicin de sus modos normales de vibracin,
es decir, el comportamiento del sistema estar en funcin de sus modos de vibracin.
3.2
Movimiento Armnico Simple (MAS)
El Movimiento Armnico Simple (MAS) es la forma ms sencilla de un movimiento peridico. Una partcula presenta un MAS cuando tiene un movimiento rectilneo en el
cual su aceleracin es siempre proporcional a la distancia a un punto fijo de la trayectoria [3], es decir, cuando su posicin esta en funcin del tiempo de una sinusoide. Si
se presenta un movimiento peridico de vibracin, sea simple o complejo (como el que
se representan en nuestro modelo fsico), puede considerarse formado por un MAS
un nmero de MAS de diferentes amplitudes y frecuencias por medio de una serie de
Fourier.
La ecuacin diferencial que representa al MAS esta dada por la ecuacin 3.1
m
x + cx + kx = F (t)
(3.1)
36
3.2. MOVIMIENTO ARMNICO SIMPLE (MAS)
donde m representa la masa, c el amortiguamiento y k es el coeficiente de rigidez. La

solucin general de este tipo de sistemas para cuando c = 0 esta dada por la siguiente
ecuacin:
x = Asen (t + )
(3.2)
donde A representa la amplitud del movimiento, es la frecuencia angular, es la fase
inicial y t representa el tiempo. Debido a que los valores de la funcin seno van de +1
a -1 el movimiento se realiza entre una regin que va desde +A hasta A (fig. 3.3) [29].
Figura 3.3: Grfica de Movimiento Armnico Simple

Si derivamos al ecuacin 3.2 con respecto al tiempo (t), podemos obtener la velocidad
a la que se desplaza la partcula [30].
v=
dx
= A cos (t + )
dt
(3.3)
Si derivamos nuevamente la ecuacin 3.3 con respecto al tiempo (la doble derivada
de 3.2), obtenemos la aceleracin de la partcula [30].
a=
dv
= 2 Asen (t + )
dt
(3.4)
donde:
r
k
k
2 =
= =
(3.5)
m
m
Puesto que la funcin seno es peridica, el movimiento se repite despus de un
tiempo 2/, lo que representa el periodo T del movimiento [30]
T =
2
2
=p
k/m
(3.6)
A partir de la ecuacin 3.6 podemos concluir que a mayor masa, ms lenta ser la
oscilacin (mayor periodo). Y si el resorte es ms blando (k grande), tambin se
tendr una oscilacin ms lenta.
Por el contrario, pequeas masas y resortes duros (menor K) darn como resultados oscilaciones rpidas (de alta frecuencia).
3.2.1
37
Sistema de un grado de libertad
Este tipo de sistemas estn restringidos de tal forma que solo pueden vibrar de una
manera, es decir, que solo se necesita una coordenada independiente para determinar
por completo la localizacin geomtrica de la masa del sistema en el espacio [3].
Si se tuviera un sistema de masa resorte (como el que se muestra en la figura 3.4) en
cual estuviese restringido para moverse horizontalmente, se dice que el sistema solo
posee un grado de libertad.
Figura 3.4: Sistema de un grado de libertad
3.2.2
Sistema de varios grados de libertad
Cuando en un sistema (como el de nuestro modelo) se necesitan n coordenadas independientes para determinar las posiciones de las masas de un sistema, el sistema se
denomina como de n grados de libertad. El movimiento de un sistema con n grados de
libertad se puede representar por n ecuaciones diferenciales [3].
Las ecuaciones diferenciales de un sistema de n masas se pueden expresar de manera
general en forma de matrices de la siguiente manera:
m11 x1 + m12 x2 + + m1n xn + k11 x1 + k12 x2 + + k1n xn = 0
m21 x1 + m22 x2 + + m2n xn + k21 x1 + k22 x2 + + k2n xn = 0
..
.
mn1 x1 + mn2 x2 + + mnn xn + kn1 x1 + kn2 x2 + + knn xn = 0
utilizando una notacin matricial la podemos representar de la siguiente manera:

m11 m11 m1n
x1
k11 k11 k1n
x1
0
m21 m21 m2n x2 k21 k21 k2n x2 0

..
..
.. .. + ..
.. . .
.. .. = ..
.
.
.
.
.
.
. . .
.
. . .
mn1 mn2 mnn
xn
kn1 kn2 knn
xn
0
[M ] {
x} + [K] {x} = 0
(3.7)
38
3.3. MODELO MATEMTICO
donde [M ] es la matriz de la masa de las partculas y [K] es la matriz de rigidez de los

resortes.
El uso de matrices en el anlisis de este tipo de sistemas de vibracin, nos ayuda
a simplificar el trabajo y a comprender el procedimiento usado en la solucin. Para
este tipo de sistemas debido al gran nmero de posibilidades que existen, encontrar
las frecuencias naturales y evaluar las respuestas, generalmente requiere mucho ms
trabajo matemtico [3].
3.3
Modelo matemtico
Una vez que se a descrito la representacin de nuestro modelo fsico, el siguiente paso
consiste en transformarlo a un modelo matemtico. Como se coment anteriormente,
el modelo del sistema se puede considerar como una cadena monoatmica lineal, formada con n partculas ( con n grados de libertad) homogneas equidistantes, con sus
extremos fijos (Fig. 3.2).
Este tipo de sistemas presentan un movimiento peridico y se componen de movimientos armnicos de diferentes amplitudes y frecuencias que se puden representar
en una matriz diagonal (Matriz ) (Fig. 3.5).
En el desarrollo de este trabajo de tesis se propone el desarrollo de dos modelos:
1. Modelo con Interaccin a primeros vecinos
2. Modelo con Interaccin a segundos vecinos
los cuales se describirn en las siguientes secciones.
Figura 3.5: Transformacin del modelo fsico al modelo matemtico
3.3.1
39
Modelo matemtico con interaccin a primeros vecinos
En una cadena de vibracin como la de nuestro caso de estudio, las partculas que
componen el sistema y que estan unidas por resortes elsticos, se ven afectadas por
el movimiento de sus partculas vecinas superior e inferior; es decir, que la partcula
i-sima solo tiene una interaccin con la partcula i 1 (partcula superior) y con la
partcula i + 1 (partcula inferior), lo que se conoce como Interaccin a Primeros Vecinos
(Fig. 3.6) [31].
Figura 3.6: Cadena de vibracin con interaccin a primeros vecinos

La fuerza resultante de una partcula de la cadena se puede obtener mediante la
ecuacin de la ley de Hooke (ecuacin 3.8)
fi =
n
X
kij (xi+j + xij 2xi )
(3.8)
j=1
donde fi es la fuerza resultante de la partcula i, xl es el desplazamiento que tiene la

partcula l y kij es el resorte que liga a las partculas i y j [31].
Tomando la ecuacin 3.8 y si consideramos que solo se tiene una interaccin a primeros
vecinos obtenemos la ecuacin 3.9
fi = k (xi+j + xij 2xi ) .
(3.9)
40
Usando la segunda ley de Newton [32], encontramos que la fuerza actuando sobre
la partcula i-sima del sistema esta dada por:
n
X
fi = mi xi .
(3.10)
i=1
Al igualar la ecuacin 3.10 con la ecuacin 3.9 obtenemos:

mi xi = k (xi+1 + xi1 2xi )
(3.11)
Sin embargo, debido a que en nuestro sistema tenemos fronteras fijas, la ecuacin
3.11 tiene restricciones para las partculas de los extremos superior e inferior. La partcula 1 del sistema no cuenta con un vecino superior y la partcula n no cuenta con un
vecino inferior, por lo que se tomaron dos partculas imaginarias 0 y n + 1 las cuales
tendrn un desplazamiento con valor de 0 (Fig. 3.7).
Figura 3.7: Cadena de vibracin con interaccin a primeros vecinos con partculas imaginarias
en los extremos
Debido a lo anterior se obtuvieron tres ecuaciones de movimiento diferentes, una

para la partcula x1 que representa la frontera superior, posteriormente la partcula
41
xi que representa todas las partculas del centro y finalmente la partcula xn que representa la frontera inferior del sistema (Fig. 3.8). De esta manera las ecuaciones de
movimiento que obtuvimos quedaron de la siguiente manera:
m1 x1 = 2kx1 + kx2
m2 x2 = kx1 2kx2 + kx3
m3 x3 = kx2 2kx3 + kx4
mi xi = kxi1 2kxi + kxi+1
mn1 xn1 = kxn2 2kxn1 + kxn
mn xn = kxn1 2kxn
Figura 3.8: Ecuaciones de movimiento con interaccin a primeros vecinos

El sistema de ecuaciones obtenido, representan el nmero de MAS por los que est
compuesto el sistema, entonces fue necesario encontrar una solucin del tipo:
xi = Ai sen (t)
(3.12)
en donde al obtener la segunda derivada de la ecuacin 3.12 con respecto al tiempo (t)
tenemos:
xi = 2 Ai sen (t) .
(3.13)
Las ecuaciones de movimiento producen un sistema de ecuaciones diferenciales.
Este sistema se puede representar en forma matricial reemplazando las ecuaciones 3.12
y 3.13 en las ecuaciones de movimiento de nuestro sistema. Si se simplifica obtenemos
el siguiente sistema de ecuaciones homogneo:
42
2k
m
k
m
..
.
0
..
.
0
k
m
2k
m
..
.
0
..
.
0
0 0
k
m
0
..
..
..
.
.
.
k
0 m
..
..
..
.
.
.
0 0
0
0
..
.
2k
m
..
.
0
0
0
..
.
..
.
0
0
..
.
k
m
0
..
..
..
.
.
.
k
0 m

0
A1
A1
A2
0 A2
.

..
.
.
. ..
.
= 2
0 Ai
Ai
.
.
..
..
. ..
2k
An
An
m
(3.14)
De esta forma, si por ejemplo se tuviera un sistema con 4 grados de libertad obtendramos el siguiente modelo matemtico con interaccin a primeros vecinos:
2k
m
k
m
0
0
k
m
2k
m
k
m
0
k
m
2k
m
k
m

0
A1
A1
0 A2
2 A2
=
k
A3
A3
m
2k
A4
A4
m
Este sistema se puede expresar de forma
M X = 2X
(3.15)
Una forma de resolver este sistema de ecuaciones diferenciales de segundo orden, es

obtener los valores propios diagonalizando la matriz M [28]. Donde 2 son los valores
propios de la matriz cuadrada y X el vector propio de los coeficientes A1 , A2 . . . An que
representan las amplitudes de los modos normales de vibracin.
3.3.2
Modelo matemtico con interaccin a segundos vecinos
Se propuso un segundo modelo en el cual, el movimiento de las partculas del sistema

son afectadas por sus dos partculas vecinas superiores y sus dos partculas vecinas inferiores; es decir, que la partcula i-sima del sistema tiene interaccin con las partculas
i2, i1 superior y con i+1, i+2 inferior, lo que se conoce como Interaccin a Segundos
Vecinos (Fig. 3.9) [28].
43
Figura 3.9: Cadena de vibracin con interaccin a segundos vecinos
Al igual que en el modelo anterior, en este modelo tambin se aplica la ley de Hooke,
por lo que tomando en consideracin las ecuaciones 3.8 y 3.10 obtenemos la ecuacin
general de nuestro sistema para una interaccin a segundos vecinos.
fi = mi xi = k (xi+1 + xi1 2xi ) + k 0 (xi+2 + xi2 2xi )
(3.16)
Si tomamos en consideracin que la constante elastica k es igual para todos los resortes que ligan a las partculas del sistema k = k 0 . De la ecuacin 3.17 obtenemos:
mi xi = kxi2 + kxii 4kxi + kxi1 + kxi2
(3.17)
La ecuacin (3.17) presenta restricciones para las partculas 1, 2, n1 y n, ya que estas

no cuentan con todas su partculas vecinas (2 vecinas superiores y 2 inferiores). Para
afrontar estas restricciones, igual que en el modelo de primeros vecinos se propusieron
agregar partculas imaginarias en los extremos con un valor de desplazamiento 0 (Fig.
3.10).
44
Figura 3.10: Cadena de vibracin con interaccin a segundos vecinos con partculas imaginarias en los extremos
Tomando en consideracin lo anterior, se obtuvieron 5 ecuaciones distintas para este
modelo. La primer ecuacin es para la partcula x1 que representa la frontera superior y
que no cuenta con sus dos vecinas superiores. La segunda ecuacin es para la partcula
x2 la cual slo cuenta con una partcula vecina superior. La tercera ecuacin representa
al movimiento de todas las partculas del centro del sistema xi y que cuenta con sus
dos vecinos superiores e inferiores. La cuarta ecuacin es para la partcula xn1 que
es la penltima del sistema y que slo cuenta con una partcula vecina inferior. Por
ltimo la ecuacin para la partcula xn que representa la frontera inferior del sistema y
que no cuenta con sus dos vecinas inferiores. As las ecuaciones que se obtienen son:
m1 x1 = 4kx1 + kx2 + kx3
m2 x2 = kx1 4kx2 + kx3 + kx4
mi xi = kxi2 + kxii 4kxi + kxi1 + kxi2
mn1 xn1 = kxn3 + kxn2 4kxn1 + kxn
mn xn = kxn2 + kxn1 4kxn
45
De la misma forma que en el modelo de primeros vecinos, estas ecuaciones de movimiento tambin generan un sistema de ecuaciones diferenciales. Al aplicar la solucin general 3.12 y 3.13 obtenemos el siguiente sistema de ecuaciones representado en
forma matricial:
4k
m
k
m
k
m
k
k
m
m
0
0
4k
k
k
0
m
m
m
k
4k
k
k
m m m m
..
..
..
..
.
.
.
.
0
0
0
0
0
0
0
0
...
0
0
0
..
.
0
0
0
..
.
0
0
0
..
.
0
0
0
..
.
Ai
A2
A3
..
.
Ai
A2
A3
..
.
..
= 2
k
k
4k
k
0
An1
m m m m An1
k
k
4k
An
An
0
0 m
m
m
(3.18)
Un ejemplo de este modelo para un sistema con 5 partculas quedara representado
de la siguiente manera:

4k
k
k
m
m
0
0
A1
A1
m
4k
k
k
k
m m 0 A2
A2
m
m
4k
k
k
= 2 A3
k k
A
3
m
m
m
m

m
4k
k
A4
0 k k
A4
m
m
m
m
k
k
4k
A5
A5
0
0 m m m
Este sistema de ecuaciones diferenciales correspondiente al modelo matemtico con
interaccin a segundos vecinos, tambin puede resolverse obtiendo los valores propios
y el vector propio correspondientes que representan las frecuencias y las amplitudes
de los modos normales de vibracin respectivamente.
3.4
Resumen
En este captulo se detallan los pasos para la construccin del modelo del sistema. En
la primer seccin se describi la forma en como fue representado el modelo fsico del
sistema. El modelo que se obtuvo se trat como un sistema de modos normales de
vibracin con n grados de libertad.
En la segunda seccin se describe la forma ms sencilla del movimiento peridico, el
Movimietno Armnico Simple (MAS). Se analiz la relacin del MAS con nuestro modelo fsico y se explic la ecuacin diferencial que representa a este tipo de movimiento.
Adems se present la solucin general para los sistemas que estn conformados por
un MAS o un nmero de MAS.
Finalmente en la tecer seccin se present la transformacin del modelo fsico del
sistema al modelo matemtico. Se desarrollaron los modelos matemticos correspondientes a la interaccin de primeros y segundos vecinos. Con estos modelos se obtuvo
un sistema homogneo de ecuaciones diferenciales para cada uno de ellos.
46
3.4. RESUMEN
Captulo 4
Diseo del sistema
La duda es la madre de la invencin.
Galileo Galilei
l obtetivo de este trabajo de investigacin es el diseo y construccin de

un sistema que nos permita predecir el comportamiento de las estructuras
que se ultizan para la extraccin de petrleo en aguas profundas. En los
captulos anteriores describimos el planteamiento del problema, los objetivos planteados, los antecedentes del sistema y finalmente el desarrollo de
nuestro modelo fsico y matemtico.
En este captulo planteamos el diseo de nuestro sistema, el cual nos
proporciona una idea general de la funcionalidad de ste a travs de sus
componentes [33].
Generalmente los sistemas de simulacin se descomponen en subsistemas y estos a su vez en mdulos, los cuales proporcionan algn conjunto
de servicios. El proceso de diseo que se lleva a cabo en este captulo, nos
ayudar a identificar los subsistemas y mdulos del sistema, as como la
comunicacin que existe entre ellos.
El objetivo de este captulo es dar a conocer el diseo de los diferentes
subsistemas y mdulos por los que est compuesto nuestro sistema. En
cada uno de estos mdulos nos enfrentamos a diversos problemas y para
resolverlos fue necesario usar diversas herramientas y mtodos numricos,
los cuales se describen en las siguientes secciones.
Un diseo modular como el que se usa, reduce la complejidad del sistema, facilita los cambios y la deteccin de errores en los mdulos. Adems,
el mantenimiento de los mdulos es ms sencillo ya que atacan problemas especficos. Este diseo en mdulos da como resultado una implementacin ms sencilla [34].
47
48
4.1. DISEO CONCEPTUAL
4.1
Diseo Conceptual
Al analizar los requerimientos del sistema descritos en nuestro planteamiento del problema y los objetivos de la investigacin, se propuso dividir el simulador en dos subsistemas (Fig. 4.1)
Figura 4.1: Divisin del simulador en dos subsistemas

Subsistema para el cmputo de datos: El primer subsistema es el encargado de
resolver el modelo matemtico del riser propuesto en el captulo anterior. Adems,
genera las condiciones de corrientes aleatorias a las que ser sometido a prueba el sistema. Finalmente crea archivos con los resultados de las pruebas, que nos permiten
determinar el comportamiento dinmico del riser en las condiciones de corrientes generadas. El subsistema de cmputo de datos est formado por los siguientes mdulos
(Fig. 4.2):
1. Mdulo del modelo matemtico del riser.
2. Mdulo del modelo hidrodinmico de corrientes
Figura 4.2: Mdulos que componen el subsistema de cmputo de datos

Subsistema de Visualizacin de resultados: Es el encargado de tomar los resultados generados por el primer subsistema para convertirlos en imgenes y animaciones,
CAPTULO 4. DISEO DEL SISTEMA
49
que permitan representar de una forma sencilla el comportamiento del sistema. Este
subsistema est conformado por los siguientes mdulos (Fig. 4.3):
1. Grfica de Resultados.
2. Animacin de partculas 2D.
3. Visualizacin Interactiva 3D.
Figura 4.3: Mdulos que componen el subsistema de Visualizacin

En resumen el diseo general del sistema est estructurado como se muestra en la
figura 4.4.
Figura 4.4: Representacin general del sistema
4.2
Modelo del riser (cadenas de vibracin)
Como vimos en el Captulo 3, el modelo de nuestro sistema es tratado como una cadena
de vibracin monoatmica, formada por partculas de masa equidistantes unidas por
4.2. MODELO DEL RISER (CADENAS DE VIBRACIN)
50
resortes. De este modelo matemtico se obtuvieron dos matrices de movimiento, una

para las partculas con interaccin a primeros vecinos (Ecuacin 3.14) y la segunda
para partculas con interaccin a segundos vecinos (Ecuacin 3.18).
Las matrices de movimiento para los dos modelos matemticos que se desarrollaron,
nos generan un sistema de ecuaciones diferenciales de segundo orden. Para resolver
este sistema de ecuaciones, se obtienen el vector y los valores propios correspondientes
al sistema, que representan las frecuencias y las amplitudes de las partculas respectivamente.
4.2.1
Problema de valores propios
El problema de calcular valores propios () se puede presentar en diversos campos

de ingeniera y las matemticas, como son: ecuaciones diferenciales, estabilidad de
sistemas lineales, sistemas elctricos, ect. Los valores y vectores propios se definen
formalmente de la siguiente manera:
Si A es una matriz de nxn. Un vector v distinto de cero es un vector propio de A si
para cierto escalar de ,
Av = v
(4.1)
donde es un valor propio de A asociado con el vector v.
Como se defini en el captulo 3, el sistema de ecuaciones homogneo que se obtuvo
para los modelos de primeros y segundos vecinos se puede representar de la forma:
M X = X
(4.2)
La ecuacin 4.2 se puede representar como una Ecuacin diferencial para los Modos
Normales de Vibracin [31]
= X
X
(4.3)
donde representa las frecuencias naturales o de resonancia del sistema.

Para resolver el sistema (Ecuacion 4.2 y 4.3) es necesario determinar las frecuencias
naturales o de resonancia y los modos normales de vibracin. Esto implica encontrar
los valores propios () y los vectores propios (X) de la matriz de movimiento (M ).
Debido que nuestro sistema est representado por una cadena de vibracin homognea, es posible la resolucin del problema por medios analticos. De otra forma,
sera necesario utilizar mtodos numricos para obtener espectros de frecuencia y de
esta manera obtener los modos normales de vibracin.[31].
Las soluciones obtenidas de las matrices de movimiento (3.14 y 3.18), son llamadas
vibracin libre o frecuencias de resonancia del sistema (valores propios 1 , 2 , . . . , n )
y dan la variacin del movimiento de las partculas del sistema por unidad de tiempo.
Por su parte los elementos del vector propio (X) se denominan coordenadas normales
del sistema. Cada una de estas coordenadas corresponde a una vibracin del sistema
con una sola frecuencia. Por lo tanto al conjunto de los componentes de la oscilacin
51
que conforman el vector propio se denominan Modos Normales de Vibracin[32]. Particularmente en nuestro sistema todas las partculas oscilan con la misma frecuencia y
fase, por lo que la vibracin total de nuestro sistema ser el resultado de la suma de los
modos normales (superposicin de los modos de vibracin).
4.2.2
Anlisis de Fourier
En diversos problemas de anlisis de ingeniera, es necesario trabajar con funciones

peridicas, es decir, funciones que se repiten cada cierto intervalo de tiempo y donde
este intervalo es conocido como periodo. Las funciones trigonomtricas son funciones
peridicas se utilizan comnmente y son relativamente fciles de trabajar ya que tienen
la propiedad de ortogonalidad. Debido a estas ventajas, es recomendable expandir la
funcin peridica en funciones trigonomtricas, donde esa expansin se conoce como
serie de Fourier [35]. El anlisis de Fourier es la clave para el anlisis de la seales de
vibracin [36].
Debido a que las ecuaciones obtenidas para nuestro sistema son ecuaciones lineales,
podemos aplicar el princpio de superposicin que nos dice que la suma de las soluciones
obtenidas es tambin una solucin. Sin embargo debido a que nuestro sistema esta
compuesto de n modos vibracin y n frecuencias naturales, se debe de considerar la
suma de n soluciones armnicas. Esto se obtiene mediante una solucin de la Serie de
Fourier [37].
X
1
[an cos (n t) + bn sin (n t)]
f (x) = a0 +
2
n=1
(4.4)
donde se conoce como la frecuencia, an (n = 1, 2, . . .) y bn (n = 1, 2, . . .) son conocidos

como los coeficientes de Fourier de la ecuacin f (x)
De esta manera, por medio de una serie de Fourier, se puede determinar que un movimiento peridico de vibracin esta formado por un nmero de MAS de diferentes
amplitudes y frecuencias.
Bajo ciertas circunstancias las series de Fourier se pueden simplificar. Esto se da
cuando se tiene una funcin peridica par (Ecuacin 4.5) y otra impar (Ecuacin 4.6)
que satisfagan ciertas relaciones de simetra, entonces la serie de Fourier puede no
contener cualquiera de los dos trminos [35].
f (t) = f (t)
(4.5)
f (t) = f (t)
(4.6)
Tomando en consideracin lo anterior para nuestro modelo, las funciones peridicas

las podemos expandir a una serie de senos de Fourier (Ecuacin 4.7).
X
1
f (x) = a0 +
[bn sin (n t)]
2
n=1
(4.7)
52
4.3. MODELO DE CORRIENTES BASADO EN MONTE CARLO
De esta manera, la ecuacin 4.7 se va aplicar a nuestro modelo para obtener el comportamiento general del sistema con sus frecuencias naturales de resonancia y con las
frecuencias aleatorias generadas con nuestro modelo de corrientes basada en el mtodo
de Monte Carlo (el cual se describir ms adelante).
4.3
Modelo de corrientes basado en Monte Carlo
Debido a que el sistema se ve afectado por un fenmeno natural (en nuestro caso las
corrientes marinas), es necesario el desarrollo de un modelo estocstico de corrientes,
el cual es independiente al modelo del riser. Este modelo genera condiciones de fuerza
y frecuencias aleatorias de corrientes que afectan al sistema, y se utiliza para probar
que el modelo matemtico que se desarroll funciona correctamente en diversas situaciones.
Para poder desarrollar este tipo de modelos estocsticos, es necesario la generacin
de nmeros aleatorios que permitan la representacin de fenmenos como el de las
corrientes marinas, ya que para este tipo de sistemas es poco prctico poder examinar
todos los casos o condiciones posibles.
Antes de comezar a generar las frecuencias aleatorias, fue necesario acondicionar
nuestro modelo fsico de tal forma, que nos permitiera producir diversas condiciones
de corrientes en diferentes segmentos del riser.
Para llevar a cabo lo anterior, adems de la divisin en partculas, el riser se dividi
en varios segmentos. Esta divisin nos va a permitir determinar condiciones especficas para la generacin de frecuencias en un segmento determinado (Fig. 4.5). Lo
anterior nos va a permitir simular condiciones de corrientes en dos sentidos diferentes.
53
Figura 4.5: Divisin del riser en partculas y segmentos

Una vez que se hicieron las modificaciones a nuestro modelo fsico, se dise el
proceso encargado de generar las frecuencias aleatorias, a las cuales ser sometido
nuestro sistema para analizar su comportimiento.
Para la generacin de los nmeros aleatorios utilizamos un mtodo numrico no
determinstico conocido como Mtodo Monte Carlo, con el cual logramos simular las
variables aleatorias que representan la fuerza y frecuencias de exitacin del sistema.
4.3.1
Mtodo Monte Carlo
El mtodo Monte Carlo es utilizado en diferente reas de la ciencia (fsica, matemticas, biologa, etc.), sin embargo en un inicio este mtodo no tuvo mucho auge debido
a la gran cantidad de clculos que se necesitan realizar para para resolver diversos
problemas. Fue hasta despus de la invencin de la computadora que tuvo una gran
impacto, ya que estas herramientas permiten realizar un gran nmero de clculos en
un experimento numrico y se tiene una gran fiabilidad en sus resultados. Con ayuda
de las computadoras la solucin de un problema o experimento numrico en el que se
implemente el mtodo Monte Carlo, es mucho ms rpido [38].
La funcin principal del mtodo Monte Carlo consiste en la generacin de nmeros
aleatorios. Para generarlos, este mtodo requiere utilizar una distribucin de probabilidad, con la cual se pueden crear muestras con un comportamiento determinado
[38].
La implementacin de este mtodo, fue necesario para generar el desalojo de cargas de frecuencia para cada una de las partculas que componen nuestra cadena de
vibracin.
54
4.3. MODELO DE CORRIENTES BASADO EN MONTE CARLO
Para cada una de las partculas de masa del riser es necesario generar n veces un
mismo vector que contenga m nmeros aleatorios, es decir, que para cada partcula se
necesitan m nmeros aleatorios (frecuencias) y cada uno de estos nmeros se deben
generar n veces (Fig. 4.6). Cada uno de estos nmeros del vector representan las frecuencias que afectan a cada una de las partculas en una unidad de tiempo.
Figura 4.6: Generacin de frecuencias aleatrias para cada partcula, con diferenes condiciones
en cada segmento
La generacin de cada uno de estos nmeros aleatorios se realiza en base a una
Distribucin Normal, en la cual los valores de media () y varianza ( 2 ) determinan
las condiciones para generar cada uno de estos nmeros aleatorios (esto se describir
ms a detalle en la siguiente seccin). Es importante mencionar que cada uno de los
segmentos en los cuales se dividi nuestro sistema, deber tener diferentes valores de
media y varianza con lo cual representaremos diferentes condiciones de corrientes.
Cabe mencionar que nuestro modelo hidrodinmico de corrientes debe estar basado
en los modos normales de vibracin de nuestro sistema. Por esta razn, las frecuencias
aleatorias para cada una de nuestras partculas se obtiene del producto del vector de
nmeros aleatorios por la frecuencia normal de vibracin de cada partcula; es decir, el
vector de cada partcula se multiplica por su valor propio correspondiente (ecuacin
4.8).

1
n1
2 n2
. .
. .
. .
=
i ni
. .
.. ..
m
Am
55
(4.8)
La aproximacin de resultados reales en el mtodo Monte Carlo, depende de diversos

factores como lo es la exactitud de los valores, el nmero de experimentos realizados
(en ese caso, el tamao del vector de nmeros aleatrios) y las varianzas de la muestra
[38].
4.3.2
Distribucin normal
La distribucin Normal o Gaussiana, es la distribucin de probablidad que ms se utiliza

para modelar fenmenos naturales. Esta distribucin cuenta con una grfica en forma
de campana, la cual se denomina curva normal y es simtrica respecto a un parmetro
[39] (Fig. 4.7).
Figura 4.7: Campa de la distribucin normal.

Una variable aleatoria normal tiene una ecuacin matemtica que depende de dos
parmetros: media () y desviacin estandar (). n (x : , ) y su funcin de densidad
esta dada por:
2
1
n (x : , ) =
e(1/2)[(x)]
(4.9)
2
Estos dos parmetros determinan la forma como se comporta la grfica de la curva.
Esta distribucin de probabilidad nos permite describir una gran cantidad de fenmenos naturales, industriales, econmicos y de investigacin. Cuando este tipo de
fenmenos son estudiados, las muestras que se toman de ellos tienden hacia una distribucin normal, conforme se aumenta el tamao de las muestras [38].
Con esta distribucin de probabilidad vamos a generar los nmeros aleatorios que
se requieren para implementar el mtodo Monte Carlo. Como se coment en la seccin
56
4.4. VISUALIZACIN
anterior, estos nmeros aleatorios se deben genenar para cada una de las partculas que
conforman el sistema, en donde los valores de media y varianza deben ser diferentes
para cada uno de los segmentos en los que se dividi el sistema (Fig. 4.6).
Para la generar estos nmeros aleatorios, primero es necesario definir el rango en
el cual deben estar los valores generados. Para el caso de nuestro sistema, los valores
debern estar en un rango de 1 < x < 1 (Fig. 4.8). Esto ser una restriccin al
momento de solicitar los valores de y .
Figura 4.8: Distribucin normal con rangos de 1 a 1.
4.4
Visualizacin
Una vez que se tuvo el diseo de los mdulos del primer subsistema, se diseo el
segundo subsistema, el cual permite la visualizacin de los resultados obtenidos por
los mdulos anteriores.
La visualizacin cientfica nos permite transformar la informacin y los datos generados por los mdulos anteriores en imgenes y animaciones, que permitan expresar
y comunicar de una manera ms sencilla el comportamiento de nuestro sistema. La
simulacin del fenmeno que estamos desarrollando nos va a producir una gran cantidad de informacin numrica, por lo que resulta muy complicado observar e interpretar los resultados. Afortunadamente, en la actualidad se cuenta con un gran poder de
cmputo y se han desarrollado tcnicas que nos facilitan la interpretacin visual de la
informacin. Estas tcnicas se conocen como visualizacin cientfica [40].
Debido a que el modelo del sistema depende del tiempo, es importante el desarrollo
de animaciones y visualizaciones interactivas que permitan realizar un anlisis ms
completo de la informacin. Una buena visualizacin nos permite explicar a personas
que no son expertas en el tema lo que esta sucediendo en el fenmeno simulado. Sin
embargo, es importante mencionar que las imgenes y animaciones que se generen,
dependen directamente de los algoritmos de extraccin de datos, por esto es necesario
estar completamente seguros de su correcto funcionamiento.
Para poder interpretar los resultados de nuestro sistema se disearon tres mdulos
de visualizacin:
1. Grfica de Resultados: El primer mdulo de visualizacin consiste en el uso
57
de grficas en dos dimensiones (x, y), en el cual se van a representar los datos
generados por cada partcula del sistema. Este mdulo nos permite comparar el
comportamiento y tendencias de todas las partculas.
2. Animacin de Partculas 2D: Este mdulo consiste en una animacin en 2D en el
cual se podr observar el comportamiento de cada partcula en cada unidad de
tiempo.
3. Visualizacin Interactiva 3D: El tercer mdulo consiste una animacin 3D en el
cual vamos a poder observar el comportamiento de todo el sistema en conjunto.
Esta visualizacin nos permite cambiar el punto de vista del sistema mediante
rotaciones (de los ejes x, y y z) y acercamientos.
4.5
Resumen
En este captulo se present el proceso general de nuestro simulador. En la primera seccin se describieron los subsistemas y mdulos por los cuales esta conformado nuestro
simulador. Se analiz la interaccin entre los procesos y se desarroll un diagrama
general.
En la segunda seccin se analiz el modelo del riser y los fundamentos tericos que
se utlizaron para resolver cada una de las partes del modelo.
En la tercera seccin se describe el desarollo del modelo hidrodinmico de corrientes
en donde se describi como se lleva a cabo la implementacin del mtodo Monte Carlo
para la generacin de las frecuencias aleatorias de las corrientes marinas.
Finalmente, en la cuarta seccin se describe el subsistema de visualizacin y los mdulos por los que est conformado (Gfica de Resultados, Animacin de Partculas 2D
y Visualizacin Interactiva 3D).
58
4.5. RESUMEN
Captulo 5
Desarrollo e implementacin del
sistema
Actuar es fcil, pensar es difcil;
actuar segn se piensa es an ms
difcil.
Johann Wolfgang Goethe
n el captulo anterior tratamos el diseo de nuestro sistema, en el cual

se describen los mdulos que lo conforman y la interaccin que hay entre
ellos. Este diseo es la base para el desarrollo de nuestro sistema. Seguir
estas especificaciones de diseo reduce el riesgo de errores durante la etapa
de desarrollo.
Como se decribe en el Captulo 4, el sistema que se desarroll est formado de dos subsistemas que se encargan de realizar los clculos para
obtener los modos de vibracin del sistema y posteriormente desplegarlos
en una animacin grfica.
En este captulo se presenta la arquitectura general de nuestro sistema.
Adems, se describen los algoritmos de programacin que se implementaron para desarrollar cada uno de los mdulos que conforman el sistema.
Cada uno de los mdulos que se desarrollarn utilizan diversas herramientas de hardware y software la cuales tambin se describen en esta captulo.
5.1
Arquitectura general
Una vez que se tuvo el diseo de nuestro sistema y que se tuvieron definidos los subsistemas y cada uno de los mdulos que los conformaban, el siguiente paso consisti
en definir la arquitectura para nuestro sistema. Para esto se realiz un anlisis para
59
60
5.1. ARQUITECTURA GENERAL
definir la plataforma, el lenguaje y las bibliotecas necesarias para poder llevar a cabo
la programacin de los subsistemas de nuestro simulador. La arquitectura general de
nuestro simulador esta definida como se muestra en la figura 5.1.
Figura 5.1: Arquitectura general del sistema

En la figura 5.1 se presentan los componentes de la arquitectura para cada subsistema. Estos componentes de la arquitectura se describen a continuacin.
5.1.1
Plataforma del subsistema de cmputo de datos
Debido a que la solucin del modelo matemtico y el modelo de corrientes requieren

un gran de poder de cmputo, el subsitema de cmputo de datos se desarroll en
un servidor el cual cuenta con un procesador Quad Core de 8 ncleos. El servidor
trabaja con un sistema operativo Mac OS Server 3.1. Este sistema operativo es usado en
servidores de la familia Macintosh y es comercializado por Apple.
Por su parte el desarrollo de esta aplicaccin se realiz en lenguaje Fortan F90. Este
lenguaje nos permite tener una interaccin con la biblioteca IMSL, la cual porporciona
rutinas que permitan resolver el modelo matemtico e implementar el mtodo Monte
Carlo para generar las frecuencias aleatorias de corrientes.
5.1.2
Plataforma del subsistema de visualizacin
El desarrollo de este subsistema se llev a cabo en una computadora con menor poder
de cmputo. La computadora que se utiliz para su desarrollo cuenta con un procesador AMD con 1GB en RAM, con sistema operatvo Linux.
Para el desarrollo de esta aplicacin de visualizacin, se utilizaron diversas herramientas de software libre. En primer lugar el sistema se diseo para trabajar en
plataformas basadas en Linux. El lenguaje base para el desarrollo de esta aplicacin fue
CAPTULO 5. DESARROLLO E IMPLEMENTACIN DEL SISTEMA
61
C++ con ayuda de la biblioteca QT versin 4.3.4, la cual nos permite crear aplicaciones
con una intefaz grfica. Una de la ventajas de trabajar con este lenguaje y esta biblioteca, es su disponibilidad para varias plataformas como Mac OS, Windows, sistemas
Unix, Linux, inclusive para algunos sistemas operativos para dispositivos mviles.
Otra de las herramientas que se utilizaron para desarrollar esta aplicacin, fue la
implementacion de OpenGL (Open Graphics Library), la cual es una interfaz de software
que nos permite hacer uso eficiente de las tarjetas aceleradoras de gficos. Esta interfaz
se compone de ms de 250 funciones que permiten programar grficos en dos o tres
dimensiones, los cuales se deben de crear a partir de primitivas geomtricas como
puntos, lnales y polgonos.
5.2
Desarrollo del subsistema de cmputo de datos
Para llevar a cabo el desarrollo del subsistema de Cmputo de Datos, se elabor un

diagrama de flujo de datos (DFD) para identificar los procesos que conforman a los
dos mdulos generales de este subsistema y la relacin que existe entre ellos (Fig. 5.2).
Figura 5.2: Diagrama de flujo de datos del Subsistema de Cmputo de Datos

En la figura 5.2 podemos observar que para comenzar la ejecucin del sistema se
requieren de algunos datos inicales, los cuales sern proporcionados por el usuario.
Los datos iniciales determinan las condiciones generales para ejecutar la simulacin.
Los datos de entrada que se requieren son:
62
5.2. DESARROLLO DEL SUBSISTEMA DE CMPUTO DE DATOS

Nmero de partculas (n): Define el nmero de partculas en las que se va a dividir
el riser y sobre las cuales se realizar el anlisis de vibraciones.
Masa de las partculas (m): Especifica la masa que va a tener cada una de las
partculas del sistema. Debido a que nuestro sistema es homogneo, solo es necesario definir la masa de una de las partculas.
Constante de elasticidad (k): Define la constante de elasticidad de los resortes que
unen a las partculas de masa de nuestra cadena de vibracin.
Longitud del riser (Long) : Con este dato vamos a establecer la longitud que tendr
nuestro riser y realizaremos la distribucin de las partculas.
Numero de Segmentos (N umSeg): Especifica el nmero de segmentos en los que
vamos a dividir el riser, en donde cada uno de estos tendrn diferentes condiciones para generar las frecuencias de corrientes aleatorias.
Con estos datos iniciales podremos comenzar a ejecutar los procesos del subsistema,
los cuales generan una serie de archivos que contienen los resultados de la simulacin.
En los procesos que se muestran en el DFD de la figura 5.2, se contemplan las rutinas para generar los clculos de los dos modelos matemticos que se desarrollaron
(primeros y segundos vecinos). El proceso para el desarrollo de cada uno de estos
mdulos se describen en las siguientes secciones.
5.2.1
Proceso para obtener las matrices de movimiento
Las matrices de movimiento son el resultado de las ecuaciones obtenidas por los modelos matemticos (primeros y segundos vecinos) que se desarrollaron en el Captulo
3. A travs de ellas podemos obtener los modos normales de vibracin del sistema con
las condiciones iniciales que se definen previamente.
Para obtener estas matrices de movimiento, es necesario contar con diversos datos
iniciales como son: el nmero de partculas, la masa de las partculas y constante de
elasticidad de los resortes. Con estos datos se desarroll el proceso que se encarga de
generar las matrices de movimiento para los modelos de primeros y segundos vecinos
(Ecuaciones 3.14 y 3.18).
Para almacenar los datos de nuestras matrices de movimiento, se necesit de dos
matrices cuadradas cuyo tamao est definido por el nmero de partculas en las que
se divide el riser. Posteriormente, se desarroll un algoritmo para llenar stas matrices con los datos correspondientes a los modelos de primeros y segundos vecinos
obtenidos en el Captulo 3.
Este algoritmo se encarga de recorrer todos los elementos de la matriz de movimiento. Cuando el algoritmo identifica que est en un elemento de la diagonal principal asigna el valor correspondiente (dependiendo si el modelo es de primeros o segundos vecinos). Adems, se encarga de asignar a los elementos vecinos de la diagonal principal el valor correspondiente a un vecino. Finalmente a los elementos que
63
no pertenecen a la diagonal principal y que no son vecinos de estos elementos se les

asigna un valor de 0.
Algoritmo 1 Genera Matriz de Movimiento
Require: Nmero de partculas n, Masa m, Elasticidad k, Matriz M atrix[n][n]
1: for i = 0 to n do
2:
{Se recorren las filas de la matriz}
3:
for j = 0 to n do
4:
{Se recorren las columnas de la matriz}
5:
if i = j then
6:
{Si es un elemento de la diagonal principal}
7:
M atrix[i][j] (2 k/m){Para matriz de Primeros Vecinos}
M atrix[i][j] (4 k/m) {Para matriz de Segundos Vecinos}
8:
else if j = (i 1) j = (i + 1) j = (i 2) j = (i + 2) then
9:
{Si es un elemento vecino}
10:
M atrix[i][j] (k/m)
11:
else
12:
{Si es cualquier otro elemento}
13:
M atrix[i][j] 0
14:
end if
15:
end for
16: end for
5.2.2
Proceso para obtener los valores y vectores propios
Una vez que se tuvieron las matrices de movimiento correspondientes a los dos modelos matemticos que se desarrollaron, el siguiente paso fue resolver el sistema de
ecuaciones. Como se describi en el Captulo 4, este proceso se realiza obteniendo
los valores propios y el vector propio correspondientes a cada uno de lo sistemas de
ecuaciones desarrollados.
Para obtener los valores y vectores propios correspondientes se uso una funcin
matemtica de la biblioteca IMSL:
CALL EVCRG(N, MATIRX, LDA, EVAL, EVEC, LDVEC)
Esta funcin de la biblioteca IMSL se encarga de calcular todos los valores propios y
el vector propio de una matriz con datos de tipo real.
Los parmetros que recibe la funcin EVCRG son:
N: Nos indica el orden de la matriz (Dato de entrada).
64

MATRIX: Es un arreglo bidimensional de punto flotante que contiene los datos
de la matriz (Dato de entrada).
LDA: Valor de tipo entero que contiene la dimensin de la matriz tal como se hizo
en su declaracin (Dato de entrada).
EVAL:Arreglo de tamao N que contienen los valores propios obtenidos de la
matriz (Datos de salida).
EVEC: Contiene la matriz de vectores propios obtenidos (Datos de Salida).
LDVEC: Valor de tipo entero que contiene la dimensin de la matriz de vectores
propios EVEC (Datos de entrada).
La funcin se encarga en primer lugar de balancear la matriz utilizando una transformacin ortogonal. Para el clculo de los valores y vectores propios se utiliza el
algoritmo QR, que consiste en obtener el producto de una matriz ortogonal por una
matriz triangular superior (ecuacin 5.1) [41].
M AT RIZ = QR
(5.1)
donde Q representa la matriz ortogonal y R es la matriz triangular superior.
5.2.3
Proceso para obtener modos normales de vibracin
Despus de realizar el clculo de los valores y vectores propios correspondientes a

las matrices de movimiento de los modelos que se desarrollaron, el siguiente proceso
consisti en obtener los modos normales de vibracin y la superposicin de los movimientos del sistema por cada unidad de tiempo.
Para esto fue necesario implementar un algoritmo, el cual ejecuta un ciclo que representa el tiempo. En cada iteracin se recorre el arreglo de valores y vectores propios
para calcular los modos normales de vibracin del sistema aplicando la ecuacin 5.2.
xi = Ai sin (t)
(5.2)
En este mismo ciclo se calcula la superposicin de estos modos normales para cada
partcula. Para realizar este clculo se implemento la ecuacin 5.3.
xi =
[Ai sin (n t)]
(5.3)
n=1
Para almacenar los resultados de los clculos, el algoritmo requiere de una matriz
para los modos normales de vibracin y un vector para almacenar los resultados de
la sumatoria. Sin embargo, al final de nuestro algoritmo estos resultados son almacenados en archivos los cuales se utilizan en los mdulos de visualizacin para hacer la
representacin grafica. Este proceso queda representado en el algoritmo 2.
65
Algoritmo 2 Modos Normales de Vibracin

Require: Matriz de vectores propios para primeros y segundos vecinos EVECM1,
EVECM2. Vector de valores propios para primeros y segundos vecinos EVALM1,
EVALM2. Nmero de partculas n. Lmite de Tiempo lim.
1: {Se declaran las matrices para almacenar los modos normales de vibracin}
2: real ResultM 1[n][n], ResultM 2[n][n]
3: {Se declaran los vectores para almacenar la superposicin de las vibraciones}
4: real SumaM 1[n], SumaM 2[n]
5: for t = 0 to lim do
6:
{Se recorren todos los valores propios}
7:
for i = 0 to n do
8:
{Se inicializa en cero el valor de la sumatoria para la superposicin}
9:
SumaM 1[i] 0
10:
SumaM 2[i] 0
11:
{Se recorre todos los vectores propios}
12:
for j = 0 to n do
13:
{Se calculan los modos normales de vibracion de cada partcula}
14:
ResultM 1[i][j] EV ECM 1[i][j] sin(EV ALM 1[i] t)
15:
ResultM 2[i][j] EV ECM 2[i][j] sin(EV ALM 2[i] t)
16:
{Se calcula la superposicin de los modos normales}
17:
SumaM 1[i] SumaM 1[i] + EV ECM 1[i][j] sin(EV ALM 1[i] t)
18:
SumaM 2[i] SumaM 2[i] + EV ECM 2[i][j] sin(EV ALM 2[i] t)
19:
end for
20:
{Se almacenan los resultados en archivos}
21:
M odN ormV ec1.txt ResultM 1
22:
M odN ormV ec2.txt ResultM 2
23:
{Se almacenan los resultados de la superposicin}
24:
SumaN ormV ec1.txt SumaM 1
25:
SumaN ormV ec2.txt SumaM 2
26:
end for
27: end for
5.2.4
Proceso para dividir el riser en segmentos
Como se describi en los Captulos 3 y 4, para realizar el clculo de los modos aleatoreos
de vibracin, primeramente se requera dividir el riser en varios segmentos, los cuales
van a determinar las condiciones para generar las frecuencias aleatorias para cada una
de las partculas de dichos segmentos.
Para llevar a cabo la divisin de segmentos se implement un algoritmo en el cual
obtiene una relacin del nmero de partculas por segmento de tal forma que cada
segmento tenga un nmero de partculas equivalente o aproximado. El proceso se
describe en el algoritmo 3.
66
Algoritmo 3 Divisin de Segmentos

Require: Vector para asignar segmentos AsignaSeg. Nmero de partculas n. Nmero
de segmentos N umSeg
1: {Variables para almacenar el nmero de partculas por segmento y restantes}
2: integer N odosSeg, N odosRestan, iN odo = 1
3: {Se obtiene el nmero de nodos por segmento}
4: N odosSeg n/N umSeg
5: {Obtenemos los nodos que restan de la divisin}
6: N odosRestan n MODULO N unSeg
7: {Se recorren los segmentos para asignarles sus nodos correspondientes}
8: for iSeg = 1 to N umSeg do
9:
{Asignamos los nodos que le tocan al segmento}
10:
for j = 1 to N odosSeg do
11:
AsignaSeg[iN odo] iSeg
12:
iN odo iN odo + 1
13:
end for
14:
{Si el segmento es menor al nmero de nodos restantes}
15:
if iSeg <= N odosRestan then
16:
{Asignamos un nodo mas a este segmento}
17:
AsignaSeg[iN odo] iSeg
18:
iN odo iN odo + 1
19:
end if
20: end for
El algoritmo requiere como datos de entrada el nmero de partculas en las que se

dividi el riser, el nmero de segmentos que se requiere tener y por ltimo un vector donde se tendr una relacin que indica a que segmento pertenece cada una de
las partculas. El proceso consiste en dividir el nmero de partculas entre el nmero
de segmentos, para asignar un nmero de parculas equivalentes. Posteriormente, las
partculas restantes de la divisin son asignadas entre los primeros segmentos, esto
implica que en algunos casos no se tenga el mismo nmero de partculas en cada segmento. La asignacin de las partculas se lleva a cabo en un arreglo de datos de tipo
entero, en donde el ndice del vector indica el nmero de la partcula y el contenido
hace referencia al segmento al que pertenece. La asignacin se puede observar en la
figura 5.3
67
Figura 5.3: Asignacin de partculas por segmento
5.2.5
Proceso para generar frecuencias aleatorias
Una vez que se tuvo la divisin del sistema en segmentos y que se realiz la asignacin
de las partculas a cada uno de ellos, el siguiente proceso consisti en implementar el
Mtodo Montecarlo para generar las frecuencias de corriente aleatorias para cada una de
las partculas.
Como se describi en el captulo 4, para generar las frecuencias aleatorias mediante
el Mtodo Montecarlo, fue necesario utilizar la distribucin Normal. Esta distribucin de
probabilidad nos va a permitir determinar las condiciones (media y varianza ) para
generar los nmeros aleatorios en cada segmento.
Para llevar a cabo este proceso se desarroll un algoritmo el cual hace uso de las
funciones de la Biblioteca IMSL, lo que nos permite generar nmeros aleatrios basados
en una distribucin normal con valores de media y varianza distintos.
El algoritmo que se implement (algoritmo 4), requiere como entradas el nmero de
segmentos, el nmero de partculas, el vector de asignacin de partculas con segmentos y los vectores de valores propios correspondientes a los dos modelos (primeros y
segundos vecinos). Lo primero que se tiene que definir es el tamao del vector de frecuencias para cada una de las partculas y el nmero de repeticiones que se tiene hacer
para generarlo.
El proceso que realiza el algoritmo se encarga de recorrer los segmentos en los que
se dividi el riser. Para cada uno de los segmentos, se producen o se solicitan (segn se
seleccione) los datos de y que sern utilizados para generar las frecuencias aleatorias. Posteriormente, se recorren cada una de las partculas que fueron asignadas al
segmento para generar su vector de frecuencias. Cada vector de frecuencias es generado un nmero de veces determinado (en este caso 1000 veces) antes de ser asignado
a la partcula. Finalmente, las frecuencias generadas deben estar basadas en los modos
normales de vibracin del sistema, por lo que el vector de nmeros aleatorios generado
68
se multiplica por el valor propio obtenido para la partcula y se almacena en una matriz
de frecuencias aleatorias.
Algoritmo 4 Frecuencias aleatorias con Monte Carlo
Require: Matriz de frecuencias aleatorias M atrixRandom1, M atrixRandom2.
Nmero de segmentos N umSeg. Vector de asignacin de segmento AsignaSeg.
Nmero de partculas n. Vector de valores propios EVALM1, EVALM2
1: {Variables para definir el nmero de frecuencias y repeticiones}
2: integer N umF rec = 100, N umRep = 1000
3: {Variables para la media, varianza y el vector para almacenar las frecuencias}
4: real xm, xstd, V ecRandom[N umF rec]
5: CALL RNSET(0) {Fijamos la semilla (funcin IMSL)}
6: {Recorremos todos los segmentos para generar los nmeros aleatrios}
7: for i = 1 to N umSeg do
8:
if Genera == 1 then
9:
CALL GeneraMediaVarianza(xm, xstd) {Genera media y varianza}
10:
else
11:
xm Leer valor, xstd Leer valor
12:
end if
13:
{Mientras el nodos pertenezca al segmento se generan las frecuencias}
14:
while AsignaSeg[iN odo] == i do
15:
{Generamos el vector de frecuencias varias veces para esa nodo}
16:
for j = 1 to N umRep do
17:
{Generamos el vector aletorio con la distribucin normal (funcin IMSL)}
18:
CALL RNNOR(NumFrec, VecRandom)
19:
{Establecemos la desviacin estndar (funcin IMSL)}
20:
CALL SSCAL(NumFrec, xstd, VecRandom, 1)
21:
{Establecemos la media de la distribucin (funcin IMSL)}
22:
CALL SADD(NumFrec, xm, VecRandom, 1)
23:
end for
24:
{Multiplicamos el vector por el valor propio correspondiente al nodo}
25:
for j = 1 to N umF rec do
26:
M atrixRandom1[iN odo][j] V ecRandom[j] EV ALM 1[i]
27:
M atrixRandom2[iN odo][j] V ecRandom[j] EV ALM 2[i]
28:
end for
29:
iN odo iN odo + 1 {Incrementamos el ndice de los nodos}
30:
end while
31: end for
Como podemos observar en el algoritmo 4 se implementaron diversas funciones de
la biblioteca IMSL, las cuales nos ayudaron principalmente a que la generacin de los
nmeros aleatrios estuvieran basados en un distribucin Normal. Las funciones que se
utilizaron en este algoritmo se describen a continuacin:
Funcin RNSET(0)
Inicializa la semilla para generar los nmeros aleatorios. Recibe como parmetro [42]:
ISEED=0: Indica que se utiliza el reloj del sistema para generar la semilla (Dato
de entrada).
Funcin RNNOR(NumFrec, VecRandom)

Esta funcin se encarga de generar nmeros pseudoaleatorios con una distribucin
normal estndar. Lo parmetros que recibe son [42]:
NumFrec: Cantidad de nmeros aleatorios que se desea generar (Dato de entrada).
VecRandom: Vector de longitud Numfrec que contiene los nmeros generados
(Datos de Salida).
Funcin SSCAL(NumFrec, XSTD, VecRandom, 1)

Esta funcin se encarga de establecer la desviacin estandar que se utiliza en la distribucin. Los parmetros que recibe son [42]:
NumFrec: Cantidad de nmeros aleatorios que se van a generar.
XSTD: Desviacin que se utiliza en la distribucin
VecRandom: Vector de longitud Numfrec que contiene los nmeros generados.
Funcin SADD(NumFrec, XM, VecRandom, 1)

Esta funcin se encarga de establecer la desviacin estndar que se utiliza en la distribucin. Los parmetros que recibe son [42]:
NumFrec: Cantidad de nmeros aleatorios que se van a generar.
XM: Valor de la media que se va a utilizar
VecRandom: Vector de longitud Numfrec que contiene los nmeros generados.
69
70
5.2.6
Proceso para calcular los modos aleatrios de vibracin
Despus de que se tuvieron generadas las frecuencias aletorias para cada una de las
partculas, el siguiente paso fue desarrollar un proceso que se encargara de calcular
los movimientos aleatorios de vibracin de cada partcula y posteriormente obtener la
superposicin de estos movimientos del sistema por unidad de tiempo.
Para calcular los movimientos aleatorios de vibracin de cada partcula, fue necesario aplicar la ecuacin de movimiento 5.2 para cada uno de los valores del vector de
frecuencias generado y posteriormente obtener el promedio de los resultados (ecuacin
5.4).
PnF rec
Ai sin (k t)
(5.4)
nF rec
donde nF rec representa el tamao del vector de frecuencias para cada partcula, Ai
es el vector propio correspondiente a la partcula, k es la frecuencia aleatoria que se
esta tratando y t representa el tiempo.
Por su parte, para calcular la superposicin de estos movimientos fue necesario desarrollar un proceso mucho ms complicado. Este clculo consiste en obtener la sumatoria de los promedios de aplicar la ecuacin de movimiento 5.2, con todos los vectores
de las partculas que pertenecen a su mismo segmento (ecuacin 5.5). Lo anterior se
puede observar ms claramente en la figura 5.4.
xi =
k=1
nP arSeg
xi =
X PnF rec Ai sin (k t)

k=1
l=1
nF rec
(5.5)
Figura 5.4: Clculo de la superposicin de frecuencias aletorias

donde nP arSeg representa el nmero de partculas en el segmento, N umF rec es
el tamao del vector de frecuencias, Ai es el vector propio de la partcula, k es la
frecuencia aleatoria actual y t el tiempo.
71
Algoritmo 5 Movimientos Aleatorios de Vibracin

Require: Vector de valores propios EVALM1, EVALM2. Matriz de vectores propios
EVECM1, EVECM2. Nmero de partculas n. Lmite de Tiempo lim. Nmero de
Segmentos N umSeg
1: {Vectores para almacenar los resultados del movimiento y superposicin}
2: real ResultM 1[n], ResultM 2[n], SumaM 1[n], SumaM 2[n]
3: {Hacemos la divisin de segmentos y generamos las frecuencias aleatorias}
4: CALL AsigaSegmentos(AsignaSeg, n, NumSeg)
5: CALL GeneraFrecuenciasMC(NumFrec, Matrix1, Matrix2, NumSeg, AsignaSeg, n,
EVALM1, EVALM2)
6: for t = 0 to lim do
7:
{Se recorren todas las parculas del sistema}
8:
for i = 1 to n do
9:
ResultM 1[i] 0, ResultM 2[i] 0
10:
{Se recorren las frecuencias generadas para esa partcula}
11:
12:
{Obtenemos el movimiento de la partcula con la frecuencia aleatoria}
13:
ResultM 1[i] ResultM 1[i] + EV ECM 1[i][j] sin(M atrix1[i][j] t)
14:
ResultM 2[i] ResultM 2[i] + EV ECM 2[i][j] sin(M atrix2[i][j] t)
15:
end for
16:
{Obtenemos el promedio de los movimientos de la partcula}
17:
ResultM 1[i] ResultM 1[i]/N umf rec, ResultM 2[i] ResultM 2[i]/N umf rec
18:
{Proceso para calcular la superposicin por segmento}
19:
SumaM 1[i] 0, SumaM 2[i] 0, iAsig 1
20:
{Recorremos el vector de asignaciones.}
21:
while AsignaSeg[i] <= AsignaSeg[iAsig] iAsug <= n do
22:
{Si la particin es igual a la particin de la partcula actual}
23:
if AsignaSeg[i] <= AsignaSeg[iAsig] then
24:
SumaT 1 0, SumaT 2 0 {Variables temporales para promedios}
25:
{Se recorren las frecuencias generadas}
26:
27:
{Obtenemos el movimiento de la partcula con la frecuencia aleatoria}
28:
SumaT 1 SumaT 1 + EV ECM 1[i][j] sin(M atrix1[iAsig][j] t)
29:
SumaT 2 SumaT 2 + EV ECM 2[i][j] sin(M atrix2[iAsig][j] t)
30:
end for
31:
{Obtenemos la suma de los promedios}
32:
SumaM 1[i] SumaM 1[i] + (SumaT 1 N umF rec)
33:
SumaM 2[i] SumaM 2[i] + (SumaT 2 N umF rec)
34:
{Se almacenan los resultados en archivos}
35:
M ovAleaV ec1.txt ResultM 1, M ovAleaV ec2.txt ResultM 2
36:
SumaAleaV ec1.txt SumaM 1, SumaAleaV ec2.txt SumaM 2
37:
end if
38:
end while
39:
end for
40: end for
72
5.3. DESARROLLO DEL SUBSISTEMA PARA LA INTERFAZ GRFICA
Como podemos observar el algoritmo 5, se encarga de realizar los clculos necesarios para obtener los movimientos aleatorios y la superposicin de estos movimiento
para cada ua de las partculas del sistema. Los resultados de los clculos obtenidos al
ejecutar este algoritmo son almacenados en archivos de texto, en los cuales los archivos
MovAleaVec1.txt y MovAleaVec2.txt corresponden a los resultados de los movimientos
aleatorios de cada partcula para el modelo de primeros y segundos vecinos respectivamente. Por otra parte los archivos SumaAleaVec1.txt y SumaAleaVec2.txt corresponden
al clculo de la superposicin de los movimientos para el modelo de primeros y segundos vecinos.
Los archivos de resultados que genera este algoritmo se utilizan para generar las grficas y animaciones del subsistema de visualizacin que se describe en las siguientes
secciones.
5.3
Desarrollo del subsistema para la interfaz grfica
Como se observa en el diagrama general del sistema (Fig 4.4), se necesit desarrollar
un subsistema independiente al de cmputo de datos, el cual permitiera hacer una
visualizacin grfica en forma de imgenes y animaciones de los resultados obtenidos
por el subsistema anterior.
Este subsistema utiliza los archivos generados por el subsistema de cmputo de datos
para poder desplegar en pantalla los resultados grficamente, de tal forma que el
usuario pueda entender e interpretar el comportamiento del sistema de una manera
ms sencilla.
El subsistema de visualizacin consiste en una aplicacin con una interfaz grfica dividida en tres segmentos, que estn desarrollados en el lenguaje de programacin C++
con ayuda de QT. Los segmentos en los que se divide la interfaz grfica representan
cada uno de los mdulo que conforman el subsistema de visualizacin descritos en el
Captulo 4.
El desarrollo de este subsistema se llev a cabo en un esquema orientado a objetos,
por lo que los mdulos que lo conforman estan compuestos por clases, las cuales se
encargan de generar las imgenes y animaciones que se despliegan en pantalla.
Para describir mejor este subsistema, se desarroll un diagrama general de las clases
que conforman cada uno de los mdulos(Fig. 5.5).
73
Figura 5.5: Diagrama de clases del sistema de visualizacin

La clase principal denominanda mainwindow se encargada de generar la interfaz
principal de nuestro sistema y la divide en segmentos para la visualizacin de cada
uno de los mdulos. Adems, se encarga de establecer los diferentes eventos que se
pueden ejecutar en la simulacin.
Abrir archivo: Permite seleccionar un archivo de resultados y generar la simulacin con sus datos.
Salir: Finaliza la aplicacin.
Ejecutar: Esta opcin se habilita cuando se abre un archivo y nos permite comenzar la simulacin.
Pausar: Cuando una animacin esta en proceso, este evento permite detenerla
temporalmente para poder continuar posteriormente con su ejecucin.
Siguiente: Esta opcin nos permite ver el comportamiento del sistema paso a
paso, siempre y cuando la simulacin se encuentre en pausa.
Detener: Detiene la simulacin.
Reinicia: Permite volver a comenzar la simulacin del archivo que se abri en un
inicio.
Graficar: Genera una grfica con los datos del archivo que se tiene abierto.
Cuando se abre un archivo, esta clase se encarga de obtener los datos de la cabecera
del archivo, con los cuales verifica que ste sea un archivo vlido para la aplicacin.
Posteriormente, obtiene la informacin adicional que proporciona la cabecera del archivo
74
como el nmero de partculas, la longitud del riser, el nmero de segmentos, el lmite

de tiempo, la masa y la constante de elasticidad para las partculas y resortes.
Una vez que se pudo abrir el archivo de resultados y que se obtuvieron los datos
de la cabecera, la clase mainwindow se encarga de dividir la interfaz principal en tres
regiones y genera instancias de las clases grafica, graphwidget y glwidget para inicializar
la animacin en cada uno de los segmentos.
Cuando se selecciona el evento para ejecutar la simulacin, se crea un timer en el
cual cada determinado intervalo de tiempo se lee una lnea de datos del archivo. Estos datos son enviados a las instancias de las clases que se crearon para que puedan
ser desplegados en cada uno de los mdulos de visualizacin. El proceso para los diversos mdulos de visualizacin y las clases que se implementaron se describen en la
siguientes secciones.
5.3.1
Mdulo de grficas de resultados
Este mdulo de visualizacin se encarga de ir graficando los resultados de cada una

de las partculas del sistema que son ledos del archivo por los clase mainwindow.
Para el desarrollo de este mdulo se implement una clase denominada grafica, la
cual cuenta con tres procesos principales (Fig. 5.6):
Figura 5.6: Procesos del mdulo de grficas de resultados

1. Inicializar Escena: Este proceso se ejecuta en el momento de abrir un archivo de
resultados y generar una instancia de la clase grafica. En este proceso la clase se
encarga de construir toda la escena donde se van generar las grficas. A partir
del nmero de nodos, la clase asigna un espacio para cada una de las partculas
del sistema en donde se van poder graficar los resultados ledos del archivo.
2. Graficar resultado: Cada que se cumple el intervalo de tiempo del timer en la
clase mainwindow, se realiza una lectura al archivo de datos. Los datos ledos son
almacenados en una cadena y enviados a una funcin de la clase grafica. Esta clase
se encarga de tomar la cadena de resultados y dividirla de manera que a cada
partcula corresponda un resultado. Cada uno de estos resultados representa un
punto en la grfica que se genera para cada partcula, por lo que despus de hacer
la divisin, se manda pintar un punto en el espacio asignado para cada grfica
(Fig. 5.7).
75
3. Aumentar y disminuir tamao: Debido a la cantidad de grficas que se tienen

que generar y al lmite del espacio que se tiene asignado para cada ellas, la clase
grafica cuenta con rutinas que permiten aumentar o disminuir el tamao de la
escena en general. Para esto se implement una funcin que captura eventos del
teclado, de tal forma que cuando se presiona la tecla + aumenta de tamao los
componentes de la escena. Por el contrario cuando se selecciona la tecla - los
elementos disminuyen su tamao.
Figura 5.7: Simulacin del mdulo de grficas

Adems el mdulo del grficas anterior, la aplicacin cuenta con un mdulo adicional que nos permite generar las grficas de todas las partculas de forma inmediata.
Para desarrollar este proceso se necesit impementar dos clases adicionales plotter y
plotsetings. Esta opcin se ejecuta cuando se selecciona el evento graficar de la interfaz principal. El proceso se encarga de realizar una lectura al archivo de resultados
independiente a la lectura de la simulacin para graficar todos los resultados de las
partculas (Fig 5.8).
Figura 5.8: Interfaz para generar las grficas de las partculas
5.3.2
Mdulo de animacin de partculas 2D
Adems de las grficas, se requeri un mdulo que nos permitiera observar el comportamiento independiente de cada una de las partculas del sistema. Para cumplir
con este requerimiento, se desarroll un mdulo que genera una animacin en 2D y
76
nos permite analizar el movimiento independiente de las partculas basndonos en el

archivo de resultados ledo.
Este mdulo consta de dos clases principales graphwidget y nodo. Al igual que en el
mdulo anterior, este tambin cuenta con tres procesos principales (Fig. 5.9), los cuales
se describen a continuacin:
Figura 5.9: Procesos del mdulo de animacin de partculas
1. Inicializar Escena: Este proceso se ejecuta cuando se abre un archivo de resultados y se genera una instancia de la clase graphwidget. El proceso se encarga de
construir la escena y crear los elementos necesario para realizar la animacin de
las partculas. Tomando como base el nmero de partculas del sistema, la clase
crea un arreglo de nodos. Cada uno de estos nodos representan una instancia de
la clase nodo, lo que nos va a permitir darles movimiento de manera independiente.
2. Mover Nodos: Al igual que en mdulo anterior, cuando se cumple el intervalo de

tiempo del timer, se envan los datos ledos a una funcin de la clase graphwidget,
la cual se encarga de separar los resultados por partculas. Con estos datos se
actualiza la posicin de los nodos en la escena. Esto se lleva a cabo en intervalos
de tiempo definidos por el timer, lo que da una perspectiva de movimiento en las
partculas (Fig. 5.10).
3. Aumentar y disminuir tamao: En este mdulo tambin fue necesario implementar una rutina que permita interactuar con el tamao de los elemento de la
escena. Para esto se implemet una funcin que captura los eventos del teclado
y que nos permite aumentar (tecla +) y disminuir (tecla -) el tamao de los componente de la escena.
77
Figura 5.10: Animacin del partculas
5.3.3
Mdulo de visualizacin interactiva 3D
Este mdulo nos va permitir analizar el funcionamiento general del sistema. Para llevar a cabo desarrollo de este mdulo se implementaron rutinas del lenguaje OpenGL.
Estas rutinas nos va permitir generar un modelo en 3D de nuestro sistema y actualizarlo con los datos del archivo de resultados.
El mdulo esta conformado por una clase denominada glwidget, encargada de construir el modelo del sistema y generar la animacin. Los procesos que conforman este
mdulo se muestran en la figura 5.11:
Figura 5.11: Procesos del mdulo de visualizacin interactiva

1. Generar modelo: Este proceso inicia cuando la clase mainwindow crea una instancia de la clase glwidget. El proceso se encarga de generar la escena y construir un
modelo que representa nuestro sistema. El modelo que se construye esta dividido
en partes, las cuales van a representar las partculas del sistema.
2. Animacin del Modelo: Al igual que en los mdulos anteriores, la animacin de
nuestro modelo se realiza en cada intervalo de tiempo del timer. Los resultados
obtenidos del archivo para cada partcula, representan la posicin de los nodos
por los que se conforma el sistema. En cada intervalo de tiempo estos nodos van
actualizando sus posiciones, lo que produce el movimiento de la estructura (Fig.
5.12).
3. Rotar Vista del Modelo: Este proceso nos permite ver el comportamiento del
sistema desde distintos ngulos. Para llevar a cabo esto se implementaron fun-
78
5.4. RESUMEN
ciones que nos permiten mover el modelo del sistema en sus distintos ejes (x, y
y z). Para llevar a cabo esta rotacin, se implement una funcin que detecta los
eventos del mouse y conforme a su movimiento hace rotar el modelo. Adems,
se implementaron tres barras de rotacin que tambin nos permiten realizar este
tipo de movimientos en los tres ejes.
Figura 5.12: Aimacin del sistema en 3D
5.4
Resumen
En este captulo se describi el proceso de como se llev a cabo el desarrollo del simulador. En la primer seccin se describi la arquitectura general de nuestro sistema. Se
detallarn los lenguajes de programacin que se utilizaron durante el desarrollo del
sistema y se describieron las libreras adicionales que se utilizaron.
En la segunda seccin se describieron los mdulos y los diversos algoritmos que se
implementaron para desarrollar el subsistema de cmputo de datos.
Finalmente en la ltima seccin, se describi el proceso y las clases que se implementaron para llevar a cabo el desarrollo del subsistema del Visualizacin.
Captulo 6
Validacin y pruebas del simulador
Todos somos muy ignorantes. Lo que
ocurre es que no todos ignoramos las
mismas cosas.
Albert Einstein
a validacin y pruebas de un proceso de simulacin por computadora,

nos permite comprobar los resultados obtenidos y analizar que tan exactas
son las predicciones del comportamiento del sistema.
El objetivo de este captulo es presentar una descripcin general de las
pruebas que se realizaron del sistema y comprobar su correcta funcionalidad.
En esta captulo se describe la forma de operar los dos subsistemas que
conforman el simulador y de esta forma obtener los resultados correctos.
Adems, se presenta el diseo de las pruebas que se llevaron a cabo y los
resultados obtenidos.
6.1
Descripcin general del uso del simulador
Como se describi en los Captulos 3 y 4 (diseo y desarrollo), nuestro sistema est formado por dos subsistemas. El primero se encarga realizar el clculo de las vibraciones
y almacena los resultados en archivos de texto. Por su parte, el segundo subsistema se
encarga de visualizar los resultados generados por medio de grficas y animaciones.
Estas dos aplicaciones son indepedientes ya que se pueden estar ejecutandos en estaciones de trabajo distintas e inclusive pueden estar realizando pruebas totalmente diferentes. Sin embargo, la aplicacin encargada de la visualizacin depende de los resultados que genera la aplicacin de cmputo de datos. Por lo cual, para llevar a cabo las
pruebas se requiere ejecutar primero la aplicacin de cmputo de datos y posteriormente
79
80
6.1. DESCRIPCIN GENERAL DEL USO DEL SIMULADOR
la aplicacin de visualizacin. La forma de operar estos dos subsistema se describe en

las siguientes secciones.
6.1.1
Aplicacin de cmputo de datos
Para comenzar a utilizar esta aplicacin es necesario contar con el archivo binario ejecutable correspondiente a la aplicacin (SistComputo). Este programa puede ejecutarse
desde lnea de comandos mediante la siguiente instruccin:
./SistComputo
Una vez que se haya ejecutado este comando, en la pantalla de la consola aparecer un mensaje (Fig. 6.1) el cual nos indica que la aplicacin inici correctamente.
Posteriormente, la aplicacin solicitar los datos de las condiciones iniciales para la simulacin. Los datos que se solicitan son el nmero de nodos, la masa de los nodos (n),
la constante de elasticidad (k), la longitud (Long) y el nmero de segmentos (N umSeg).
Figura 6.1: Pantalla de inicio de la aplicacin de cmputo de datos
Si el nmero de segmentos que se proporciona es mayor al nmero de partculas

del sistema, la aplicacin desplegar un mensaje de error como el que se muestra en
la figura 6.2 y volver a pedir el nmero de segmentos hasta que se proporcione un
nmero correcto.
CAPTULO 6. VALIDACIN Y PRUEBAS DEL SIMULADOR
81
Figura 6.2: Pantalla de error en el nmero de segmentos
Despus de proporcionar los datos iniciales, se genera un directorio en el cual se

almacenarn los archivos de resultados. El nombre de este directorio esta conformado
por la fecha y la hora en que se ejecut la aplicacin (Fig. 6.3).
Figura 6.3: Pantalla de la creacin del nombre del directorio
Posteriormente la aplicacin se encarga de generar las matrices de movimiento para

los dos modelos y obtiene los valores propios y el vector propio correspondientes a
estas matrices (Fig. 6.4).
82
Figura 6.4: Pantalla de matrices de movimiento con sus valores propios y vector propio
Con estos datos el sistema comienza a realizar el clculo de los modos normales
de vibracin de las partculas y la superposicin de los movimientos (Fig. 6.5). Los
resultados obtenidos son almacenados en cuatro archivos:
1. ModNormVec1.txt: Contiene los resultados de los modos normales de vibracin
del modelo de primeros vecinos.
2. ModNormVec2.txt: Contiene los resultados de los modos normales de vibracin
del modelo de segundos vecinos.
3. SumaNormVec1.txt: Contiene los resultados de la superposicin de los modos normales de vibracin del modelo de primeros vecinos
4. SumaNormVec2.txt: Contiene los resultados de la superposicin de los modos
normales de vibracin del modelo de segundos vecinos
Figura 6.5: Pantalla del clculo de modos normales de vibracin
83
Una vez obtenidos los resultados de los modos normales de vibracin, el sistema
realiza la divisin de segmento del sistema y genera la matriz de frecuencias aleatorias.
Para esto el sistema pregunta si los valores de media y varianza para cada uno de
los segmentos se van a generar o se van a proporcionar por el usuario(Fig. 6.6).
Figura 6.6: Pantalla de los datos de media y varianza de los segmentos

Finalmente la aplicacin realiza el clculo de los modos aleatorios de vibracin de
las partculas y la superposicin de estos movimientos. Los resultados obtenidos se
almacenan en otros cuatro archivos:
1. MovAleaVec1.txt: Contiene los resultados de los modos aleatorios de vibracin
del modelo de primeros vecinos.
2. ModAleaVec2.txt: Contiene los resultados de los modos aleatorios de vibracin
del modelo de segundos vecinos.
3. SumaAleaVec1.txt: Contiene los resultados de la superposicin de los modos aleatorios de vibracin del modelos de primeros vecinos
4. SumaAleaVec2.txt: Contiene los resultados de la superposicin de los modos
aleatorios de vibracin del modelos de segundos vecinos
Para evitar que el usuario este proporcionando los datos necesarios para ejecutar la
aplicacin, se pude crear un achivo (DatosEntrada.txt) que contenga todos estos datos
de entrada y ejecutar la aplicacin de la siguente manera:
./SistComputo < DatosEntrada.txt
6.1.2
Aplicacin de visualizacin
Para comenzar a utilizar la aplicacin de visualizacin, se require del archivo ejecutable

SistVisualizacion. Este programa puede ejecutarse desde la interfaz grfica haciendo
84
doble click sobre el icono del programa o bien desde lnea de comandos mediante la
instruccin:
./SistVisualizacion
Al inicio se desplegar en pantalla una interfaz que contienen un men principal
y una barra de herramientas con diversas opciones para realizar la simulacin de un
archivo de resultados (Fig. 6.7).
Figura 6.7: Pantalla de inicio del subsistema de visualizacin

El men principal est compuesto por cuatro submens (Fig. 6.8) que contienen las
siguientes opciones:
Figura 6.8: Mens del subsistema de visualizacin

1. Archivo
(a) Abrir: Abre un archivo para visualizar sus resultados.
(b) Salir: Cierra la aplicacin.
2. Herramientas
(a) Graficar: Grafica los datos del archivo en una ventana independiente a la
interfaz principal.
85
3. Ejecutar
(a) Reiniciar: Reinicia la simulacin.
(b) Ejecutar: Comienza la simulacin.
(c) Pausa: Detiene momentneamente la simulacin
(d) Detener: Detiene totalmente la simulacin.
(e) Siguiente: Genera paso a paso la simulacin
4. Ayuda
(a) Acerca de: Contiene informacin sobre la versin del simulador.
Por su parte la barra de herramientas se compone de las siguientes secciones (Fig.
6.9):
Figura 6.9: Barra de herramientas del subsistema de visualizacin

1. Seccin Archivo: Contiene las mismas opciones que el men archivo.
2. Seccin Ejecutar: Contiene las opciones del men ejecutar
3. Seccin Herramientas: Contiene las opciones del men herramientas.
Para comenzar con la simulacin, es necesario seleccionar la opcin Abrir archivo. Al
seleccionar esta opcin, se despliega una ventana (Fig. 6.10) en la que se puede elegir
alguno de los archivos generados por el susbsistema de cmputo de datos.
Figura 6.10: Pantalla para abrir archivo.

Al abrir un archivo de datos, la intefaz se divide en tres segmentos en los que se
generan las escenas correspondientes para hacer la visualizacin de los datos (Fig. 6.11)
y se habilitan las diversas opciones del men para manejar la simulacin.
86
Figura 6.11: Pantalla de inicio de simulacin
Posteriormente se puede seleccionar la opcin ejecutar para comenzar a generar las

grficas y animaciones con los datos del archivo (Fig. 6.12). Si se desea ejecutar la
simulacin paso a paso, es necesario seleccionar la opcion pausar para detener momentneamente la simulacin y despus seleccionar la opcin siguiente.
Figura 6.12: Pantalla del proceso de simulacin.
Para obtener la grfica completa de todos los datos del archivos, se debe seleccionar
la opcin graficar, la cual despliega una ventana con la grfica generada (Fig. 6.13).
Este proceso se puede ejecutar en cualquier momento sin necesidad de detener la simulacin principal.
87
Figura 6.13: Pantalla del sistema que muestra la grfica de los datos.
Finalmente, para salir de la aplicacin se debe seleccionar la opcion salir o bien seleccionar el botn cerrar ventana de la interfaz.
6.2
Diseo de la pruebas
Para validar el funcionamiento de nuestro simulador, se disearon diversa pruebas con

diferentes condiciones. Durante estas pruebas, se verific que se ejecuten cada uno de
los procesos que conforman el sistema y que se generen resultados de acuerdo a las
condiciones establecidas.
Los datos principales que fueron considerados para el diseo de las pruebas fueron:
la longitud del riser, el nmero de partculas y el nmero de segmentos en los que se
divide el sistema.
Para determinar la longitud del sistema, se tom como referencia la profundidad
mxima que pudiera tener un riser en aguas profundas, que como se mencion en el
Captulo 1 es de 1500 mts (Fig. 2.4). Tomando como referencia esta longitud y debido a
que se quera realizar un anlisis de las vibraciones en distancias relativamente cortas,
se estableci tener una cantidad de 256 partculas, las cuales distribudas entre los 1500
metros de longitud se tendr una partcula cada 5.85 mts (Fig. 6.14).
88
6.3. RESULTADOS DE LAS PRUEBAS
Figura 6.14: Distribucin de partculas en el riser para las pruebas

Con estas condiciones iniciales se disearon tres pruebas para analizar el comportamiento del sistema en diversas condiciones. Para esto el sistema se divid en segmentos de la siguiente manera:
Sistema con 32 segmentos y 8 partculas en cada segmento.
Para estas pruebas la masa de las partculas y la constante de elasticidad se definieron
con un valor de 1. Por su parte los valores de media y varianza para generar las frecuencias de cada segmento, se calcularon de forma aleatoria con la finalidad de generar
distintas condiciones de corriente en cada uno de los segmentos.
En las pruebas realizadas, el clculo de los datos para cada partcula se hizo en un
ciclo de 1 a 3600 seg., en donde cada iteracin representa una unidad de tiempo. Esto
con la finalidad de analizar si el comportamiento de las partculas se repite despus de
un intervalo de tiempo determinado.
Con estos valores se ejecut el sistema de cmputo de datos para generar los archivos
de resultados correspondientes y posteriormente desplegar estos resultados en el sistema de visualizacin.
6.3
Resultados de las pruebas
En esta seccin se mostrarn los resultados que se obtuvieron al realizar las pruebas
que se disearon en la seccin anterior. Cada una de estas pruebas se realizaron por
separado y generaron distintos archivos de resultados.
Para agilizar la ejecucin de las pruebas, se crearon archivos con los datos de entrada
para evitar estar proporcionando los datos a travs del teclado.
6.3.1
89
Prueba con 8 segmentos
Para comenzar a realizar los clculos para esta prueba, se cre un archivo que contiene
los datos de entrada para el sistema de computo

de datos. Los datos para la ejecucin
de esta prueba se muestran en la tabla 6.1.
Datos de entrada
Nmero de partculas n
Masa de las partculas m
Constante de elasticidad k
Longitud del riser Long
Nmero de Segmentos N umSeg
256
1
1
1500
8
Tabla 6.1: Datos de entrada para la ejecucin de la prueba de 8 segmentos

Posteriormente se ejecut el sistema con los datos de entrada definidos y se generaron los siguientes archivos (Fig. 6.15).
Figura 6.15: Archivos que se generaron de la prueba de 8 segmentos

Una vez que se finaliz el cmputo de los datos, los archivos que contenan los resultados de las frecuencias aleatorias se sometieron al sistema de visualizacin. El primer
archivo que se someti al sistema de visualizacin, es el que contiene los resultados de
los modos aleatorios de vibracin para el modelo de primeros vecinos MovAleaVec1.txt.
La simulacin y la grfica que se gener en el sistema se muestra en la figura 6.16.
Figura 6.16: Resultado de la visualizacin de los modos aleatorios para primeros vecinos
90
El siguiente archivo que se ejecut en el sistema de visualizacin, es el que contiene

los modos aleatorios de vibracin para el modelo de segundos vecinos MovAleaVec2.txt.
Los resultados obtenidos en el sistema de visualizacin se muestran en la figura 6.17.
Figura 6.17: Resultado de la visualizacin de los modos aleatorios para segundos vecinos
En la tabla 6.2 se muestra una comparacin de los resultados obtenidos de algunas
de las partculas en los tiempos t100 t103 para los modelos de primeros y segundos
vecinos.
n
10
20
30
40
50
60
70
80
90
100
t100 (1er Vecinos, 2do Vecinos)

1.907847E-05 -1.992678E-04
-1.813446E-03 -1.124881E-03
-1.246394E-04 -3.265392E-03
-3.350988E-02 1.990384E-02
-4.408304E-02 -9.005854E-03
-5.551505E-02 -5.186288E-02
-7.175901E-03 1.292170E-02
-7.305894E-04 -5.214892E-03
-2.965253E-03 -4.334845E-03
2.615707E-03 2.836561E-04
Modos aleatorios de vibracin. Riser con 8 segmentos

t101 (1er Vecinos, 2do Vecinos) t102 (1er Vecinos, 2do Vecinos)
-3.929516E-05 -1.699451E-04 -9.593633E-05 -1.479500E-04
-1.809459E-03 -1.344840E-03 -1.801484E-03 -1.517096E-03
-6.290861E-05 -3.695666E-03 9.108428E-07 -4.033606E-03
-3.241201E-02 1.880272E-02
-3.122393E-02 1.757011E-02
-4.307781E-02 -8.612059E-03 -4.195805E-02 -8.158773E-03
-5.500680E-02 -5.079360E-02 -5.436418E-02 -4.938839E-02
-7.595513E-03 1.436066E-02
-7.868181E-03 1.554352E-02
-5.016610E-04 -5.079770E-03 -3.422102E-04 -4.961209E-03
-2.882498E-03 -3.916809E-03 -2.827634E-03 -3.518403E-03
2.443388E-03 3.041067E-04
2.267575E-03 3.126351E-04

-1.504193E-04 -1.325698E-04
-1.789948E-03 -1.638451E-03
6.807496E-05 -4.272896E-03
-2.994924E-02 1.621585E-02
-4.072702E-02 -7.649585E-03
-5.358897E-02 -4.765817E-02
-7.990567E-03 1.637010E-02
-2.547934E-04 -4.832680E-03
-2.796132E-03 -3.159899E-03
2.088836E-03 3.094971E-04
Tabla 6.2: Resultados de los modos aleatorios de vibracin. Riser con 8 segmentos.
Despus de visualizar los datos de los modos aleatorios de vibracin de los modelos
de primeros y segundos vecinos, se procedi a desplegar los resultados de los archivos
que contienen la superposicin de los movimientos de ambos modelos.
La visualizacin que se gener para el archivo de superposicin del modelo de
primeros vecinos (SumaAleaVec1.txt) se muestra en la figura 6.18.
Figura 6.18: Resultado de la visualizacin de la superposicin para primeros vecinos
91
Finalmente se visualiz el archivo que contiene los resultados de la superposicin

del modelo de segundos vecinos SumaAleaVec2.txt (Fig. 6.19).
Figura 6.19: Resultado de la visualizacin de la superposicin para segundos vecinos

La comparacin de los resultados correspondientes la superposicin de los movimientos para esta prueba se muestran en la tabla 6.3.
n
10
20
30
40
50
60
70
80
90
100

-7.063910E-03 -7.742336E-03
-1.410334E-02 -1.584999E-02
-2.099716E-02 -2.325837E-02
-1.14011
-1.11607
-1.39728
-1.40756
-1.63629
-1.68989
-2.510596E-02 -2.021624E-02
-2.767916E-02 -2.272061E-02
-2.983056E-02 -2.689450E-02
-1.368259E-02 -1.443152E-02
Superposicin de los movimientos. Riser con 8 segmentos

-7.019826E-03 -7.797924E-03 -6.907654E-03 -7.791796E-03
-1.401532E-02 -1.601859E-02 -1.379137E-02 -1.606765E-02
-2.086613E-02 -2.345954E-02 -2.053269E-02 -2.347956E-02
-1.11454
-1.09076
-1.08604
-1.06267
-1.36595
-1.37611
-1.33101
-1.34100
-1.59959
-1.65262
-1.55868
-1.61079
-2.362785E-02 -1.984045E-02 -2.202464E-02 -1.929876E-02
-2.604954E-02 -2.220884E-02 -2.428201E-02 -2.151776E-02
-2.807429E-02 -2.580012E-02 -2.616937E-02 -2.453260E-02
-1.591373E-02 -1.661902E-02 -1.826016E-02 -1.890768E-02

-6.726176E-03 -7.721555E-03
-1.342904E-02 -1.599171E-02
-1.999326E-02 -2.331088E-02
-1.05468
-1.03187
-1.29258
-1.30233
-1.51367
-1.56452
-2.030496E-02 -1.853681E-02
-2.238609E-02 -2.059305E-02
-2.412609E-02 -2.306439E-02
-2.073010E-02 -2.130510E-02
Tabla 6.3: Resultados de la superposicin de los movimientos. Riser con 8 segmentos
6.3.2
Al igual que en la prueba anterior, tambin se cre un archivo con los datos de entrada,
los cuales se muestran la tabla 6.4.
Datos de entrada
256
1
1
1500
16

Con estos datos de entrada se ejecut el proceso de cmputo de datos y se generaron
los archivos que se muestran en la figura 6.20.
92

Una vez que se generaron los archivos de resultados, se hizo la visualizacin de los
archivos que contienen los modos aleatorios de vibracin para los modelos de primeros
y segundos vecinos. La visualizacin correspondiente al archivo MovAleaVec1.txt se
muestra en la figura 6.21.
Figura 6.21: Visualizacin de los modos aleatorios con 16 segmento para primeros vecinos
Por su parte la visualizacin que se gener para el archivo MovAleaVec2.txt se muestra en la figura 6.22.
Figura 6.22: Visualizacin de los modos aleatorios con 16 segmento para segundos vecinos
Por su parte, la tabla 6.5 mestra los resultados obtenidos de los modos aleatorios de
vibracin que obtuvieron con el sistema para esta prueba.

n
10
20
30
40
50
60
70
80
90
100

-3.678883E-04 9.624943E-04
-6.906331E-04 4.767009E-03
5.285075E-03 5.038357E-04
2.230595E-04 4.435430E-03
3.338236E-02 3.968179E-03
4.311035E-02 2.621725E-02
-6.986971E-03 9.441948E-03
2.673260E-03 -2.097178E-03
-8.699710E-03 5.818764E-03
8.079447E-05 -9.115489E-04

-2.971953E-04 9.215801E-04
-2.171444E-04 8.696534E-04
-7.981448E-04 5.123731E-03
-9.059251E-04 5.431342E-03
5.379449E-03 2.732549E-04
5.458830E-03 5.111218E-05
7.361718E-05 4.480458E-03
-8.221404E-05 4.500663E-03
3.284801E-02 3.895328E-03
3.224033E-02 3.801157E-03
4.253648E-02 2.529160E-02
4.186080E-02 2.421246E-02
-6.866159E-03 9.452593E-03
-6.754513E-03 9.450705E-03
2.977755E-03 -1.881491E-03 3.272270E-03 -1.685780E-03
-8.751728E-03 5.919788E-03
-8.774925E-03 5.904634E-03
-2.119206E-04 -9.772019E-04 -4.824622E-04 -1.046817E-03
93
-1.312430E-04 8.034877E-04
-1.011415E-03 5.681378E-03
5.521549E-03 -1.565643E-04
-2.413793E-04 4.494882E-03
3.156222E-02 3.687246E-03
4.108621E-02 2.299237E-02
-6.651716E-03 9.434445E-03
3.556486E-03 -1.510610E-03
-8.770801E-03 5.769524E-03
-7.311178E-04 -1.119169E-03
Tabla 6.5: Resultados de los modos aleatorios de vibracin para la prueba de 16 segmentos
Posteriormente se realiz la visualizacin de los archivos que contienen los resultados de la superposicin de los movimientos. La visualizacin correspondiente al
archivo SumaAleaVec1.txt se muestra en la figura 6.23.
Figura 6.23: Visualizacin de la superposicin con 16 segmentos para primeros vecinos

Por ltimo se realiz la visualizacin del archivo que contiene los resultados de la
superposin correspondiente al modelo de segundos vecinos SumaAleaVec2.txt (Fig.
6.24)
Figura 6.24: Visualizacin de la superposicin con 16 segmentos para segundos vecinos

En la tabla 6.6 se compran los resultados correspondientes a la superposicin de los
movimientos aleatorios para la prueba con 16 segmentos.
94
n
10
20
30
40
50
60
70
80
90
100

-2.904481E-03 -2.318145E-03
-6.729351E-03 -9.387511E-03
-1.001872E-02 -1.345996E-02
3.765200E-03 7.211054E-04
0.559681
0.564643
0.655415
0.681291
-4.979879E-02 -5.186686E-02
-5.490282E-02 -5.700000E-02
-1.402523E-02 -1.247462E-02
1.246617E-02 1.139727E-02

-2.837077E-03 -2.466894E-03 -2.774331E-03 -2.628331E-03
-6.583652E-03 -8.958607E-03 -6.437876E-03 -8.563327E-03
-9.801791E-03 -1.287640E-02 -9.584764E-03 -1.233636E-02
4.458917E-03 1.388127E-03
5.127800E-03 2.014992E-03
0.549596
0.554466
0.538216
0.542966
0.643604
0.669007
0.630279
0.655050
-4.840570E-02 -5.026108E-02 -4.707154E-02 -4.873393E-02
-5.336696E-02 -5.524845E-02 -5.189607E-02 -5.358185E-02
-1.408295E-02 -1.162544E-02 -1.404520E-02 -1.073559E-02
1.037795E-02 9.235727E-03
8.338261E-03 7.125000E-03

-2.715826E-03 -2.795956E-03
-6.298355E-03 -8.200939E-03
-9.377045E-03 -1.184012E-02
5.766595E-03 2.595129E-03
0.525593
0.530193
0.615496
0.639489
-4.579417E-02 -4.728381E-02
-5.048776E-02 -5.199837E-02
-1.391439E-02 -9.815786E-03
6.348027E-03 5.066688E-03
Tabla 6.6: Resultados de la superposicin de los movimientos. Riser con 16 segmentos
6.3.3
La ltima de las pruebas que se realiz, consiste en un sistema con 32 segmentos y 8

partculas en cada uno de ellos. Los datos de entrada para esta prueba se muestran en
la tabla 6.7.
Datos de entrada
256
1
1
1500
32

Los archivos de resultados que fueron generados por el sistema de cmputo de datos
se muestran en la figura 6.25.

Al igual que en las pruebas anteriores, tambin se realizaron las visualizaciones de
los resultados correspondientes a los modos aleatorios de vibracin. La simulacin y
la grfica que se gener para el archivo MovAleaVec1.txt se muestra en la figura 6.26.
95
Figura 6.26: Visualizacin de los modos aleatorios con 32 segmento para primeros vecinos
Posteriormente se hizo la visualizacin para el archivo que contiene los resultados
para los modos aleatorios de vibracin del modelo de segundos vecinos MovAleaVec2.txt
(Fig. 6.27).
Figura 6.27: Visualizacin de los modos aleatorios con 32 segmentos para segundos vecinos
La tabla 6.8 muestra una comparacin de los movimientos aletorios de vibracin
obtenidos para los modelos de primeros y segundos vecinos.
n
10
20
30
40
50
60
70
80
90
100

6.876415E-04 -8.850906E-04
-4.036903E-04 4.386333E-04
-1.502593E-03 4.304219E-03
-1.446174E-03 -3.393252E-04
4.417054E-03 9.598906E-05
-5.210994E-04 9.901992E-03
-1.827724E-04 -4.078187E-03
-1.033875E-02 1.173887E-02
-4.143782E-04 2.221534E-03
-7.358078E-03 -1.783243E-04

6.828061E-04 -8.441567E-04 6.738878E-04 -7.997585E-04
-3.535228E-04 4.923437E-04
-3.069670E-04 5.486439E-04
-1.427628E-03 4.331979E-03
-1.330414E-03 4.311639E-03
-1.742488E-03 -1.285509E-04 -1.987391E-03 1.928559E-05
4.218118E-03 8.527814E-05
3.988178E-03 7.397072E-05
-5.576548E-04 9.453511E-03
-5.840030E-04 8.937095E-03
-5.563378E-05 -4.705289E-03 1.087742E-04 -5.362171E-03
-1.041135E-02 1.198066E-02
-1.043335E-02 1.211056E-02
-3.586058E-04 2.288717E-03
-2.719905E-04 2.301937E-03
-6.713003E-03 -1.901416E-04 -6.090183E-03 -2.019299E-04

6.608858E-04 -7.524574E-04
-2.639329E-04 6.071637E-04
-1.213468E-03 4.243405E-03
-2.176788E-03 9.841991E-05
3.730056E-03 6.492515E-05
-6.016193E-04 8.358353E-03
3.090075E-04 -6.031497E-03
-1.040436E-02 1.211814E-02
-1.560925E-04 2.268540E-03
-5.489950E-03 -2.138223E-04
Tabla 6.8: Resultados de los modos aleatorios de vibracin. Riser con 32 segmentos
Una vez que se tuvieron las visualizaciones de los modos aleatorios de vibracin
para los modelos de primeros y segundos vecinos, se procedi a generar la visualizacin de la superposicin de los movimientos para ambos modelos. El primer archivo
que se ejecut en el simulador fue el correspondiente a la superposicin de los movimientos del modelo de primeros vecinos SumaAleaVec1.txt(Fig. 6.28)
96
6.4. COMPARACIN DE RESULTADOS CON OTROS MODELOS
Figura 6.28: Visualizacin de la superposicin con 32 segmentos para primeros vecinos

Finalmente se gener la simulacin correspondiente al archivo que contiene los resultados de la superposicin de los movimientos correspondientes al modelo de segundos vecinos (Fig. 6.29).
Figura 6.29: Visualizacin de la superposicin con 32 segmentos para segundos vecinos

En la tabla 6.9 se muestran los resultados correspondientes a la superposicin de los
movimientos obtenidos para la prueba con 32 segmentos.
n
10
20
30
40
50
60
70
80
90
100

-3.002102E-04 -6.024964E-04
-2.570995E-03 -2.361001E-03
-5.866234E-03 -5.350759E-03
-9.734591E-03 -1.007392E-02
-8.916042E-04 -4.443963E-04
3.514747E-03 -1.472781E-03
4.039155E-04 1.159850E-03
7.860347E-03 1.959921E-02
1.746936E-03 -4.621875E-03
-9.515047E-02 -9.466318E-02

-2.033571E-04 -4.529762E-04 -1.048360E-04 -2.958368E-04
-2.249602E-03 -2.197622E-03 -1.936035E-03 -2.040577E-03
-5.801875E-03 -4.665426E-03 -5.672919E-03 -3.997867E-03
-9.905836E-03 -1.003439E-02 -1.004572E-02 -1.002643E-02
-1.279588E-03 -9.231197E-04 -1.746433E-03 -1.464311E-03
3.480793E-03 -1.428056E-03 3.444052E-03 -1.339760E-03
5.015741E-04 1.073712E-03
6.165056E-04 1.055959E-03
6.460517E-03 1.844117E-02
4.917956E-03 1.702852E-02
1.785556E-03 -4.204112E-03 1.833161E-03 -3.694265E-03
-9.029475E-02 -8.982857E-02 -8.556519E-02 -8.512270E-02

-5.790993E-06 -1.334611E-04
-1.629300E-03 -1.887964E-03
-5.483500E-03 -3.364217E-03
-1.014799E-02 -1.004957E-02
-2.301221E-03 -2.075825E-03
3.403790E-03 -1.207379E-03
7.455526E-04 1.110428E-03
3.268740E-03 1.538688E-02
1.887605E-03 -3.104652E-03
-8.096351E-02 -8.054730E-02
Tabla 6.9: Resultados de la superposicin de movimientos. Riser con 32 segmentos
6.4
Comparacin de resultados con otros modelos
En la literatura actual, existen investigaciones que utilizan diversos mtodos con los
cuales se puede hacer una anlisis y monitoreo de las vibraciones que sufren los riser
97
por la condiciones marinas. Estos mtodos implementan diversos modelos hidrodinmicos de corrientes, los cuales son diferentes y ms complejos al que se implement en
este trabajo de investigacin.
Debido a las diferencias ente los modelos hidrodinmicos que se implementan para
el monitoreo de los riser, es muy complicado hacer una comparacin directa de los resultados de nuestro modelo con los otros resultados. Otro de los factores que complica
la comparacin de resultados entre los modelos, es que utilizan diferentes parmetros
de entrada para realizar los experimentos y en muchas ocaciones no se cuenta con los
datos necesarios para reproducir el mismo experimento con otro modelo. Por su parte,
los simuladores que existen actualmente para el monitoreo de los riser, son desarrollados por empresas privadas y obtener una licencia de estos para realizar la comparacin
de resultados sera muy costoso.
Es posible hacer una comparacin de las frecuencias normales del sistema de nuestro
modelo con algunos otros modelos implementados actualmente. Para esto, se necesit
investigar algunas publicaciones que ofrecieran resultados numricos de sus pruebas
y los parmetros que utilizaron para poder reproducir el experimento.
En los artculos presentados por Senjanovic [43] y Sparks [44], se realiza un anlisis de las vibraciones naturales de riser con diferente tensin, a travs de un mtodo
de segmentacin. El mtodo que utilizan tiene caractersticas similares al nuestro ya
que tambin hacen una segmentacin del sistema en nodos, adems proporcionan los
parmetros que utilizaron para realizar sus pruebas. Esto nos permite reproducir las
pruebas que ellos realizaron, principalmente para obtener las frecuencias naturales del
riser.
Los parmetros que se utilizan [44] y [43] para las pruebas que realizaron se muestran en la tabla 6.10:
Caractersticas del riser
Longitud del riser L
Resistencia de flexin E (Mdulo de Young)
Peso Lineal w
Tensin superior Tt
Tensin de fondo Tb
Masa del riser m
Modos Normales n
2000m
318.6 106 N m2
3433.5N/m
7.5537 106 N
0.6867 106 N
1200Kg/m
50
Tabla 6.10: Parmetros del riser de ejemplo

Para reproducir estas pruebas, de la tabla 6.10 se identificaron los datos que necesita
nuestro modelo. Los parmetros que se requieren son: la longitud del riser (L), el
nmero de partculas(n), la constante de elasticidad (k) y la masa de las partculas
(m). Para la masa de las partculas, en la tabla de datos se proporciona la cantida de
masa por metro, por lo que fue necesario hacer el clculo para obtener la masa de cada
partcula. Los valores de los parmetros que se utilizaron para reproducir el ejemplo
con nuestro modelo se muestran en la tabla 6.11
98
6.4. COMPARACIN DE RESULTADOS CON OTROS MODELOS

Parmetros utilizados para nuestro modelo
Longitud del riser L
2000m
Constante de elasticidad k (Mdulo de Young) 318.6 106 N m2
Masa del riser m
2400000Kg
48000Kg
Modos Normales (Partculas) n
50
Tabla 6.11: Parmetros del ejemplo para nuestro modelo

Una vez que se reprodujo la prueba, fue necesario hacer la conversin de los resultados, ya que las frecuencias obtenidas por nuestro modelo estan en grados y los
resultados reportados esta en radianes.
La obtensin de las frecuencias en nuestro modelo se hace por medio de los valores
propios del sistema, el cual se representa con la ecuacin 6.1
M X = 2X
(6.1)
Para obtener la frecuencia natural del riser es necesario obtener el valor de y hacer
la conversin a radianes mediante la ecuacin 6.2

2
(6.2)
=
180
De esta forma se obtuvieron los resultados para compararlos con los resultados reportados en el artculo [43]. La comparacin de los resultados se muestra en la tabla
6.12
n
1
2
3
4
5
10
20
30
40
50
Exacta
0.07973
0.16140
0.24273
0.32395
0.40511
0.81072
1.62170
2.43263
3.24353
4.05443
Frecuencias naturales del riser n

Resultados Literatura
Nuestros Modelos
Asinttica Segmentacin FEM
1er Vecinos 2do Vecinos Promedio
0.08108
0.08115
0.07983 0.08761
0.19357
0.10050
0.16217
0.16238
0.16176 0.17506
0.38670
0.20158
0.24326
0.24399
0.24370 0.26238
0.57847
0.30242
0.32435
0.32610
0.32602 0.34946
0.76840
0.40305
0.40544
0.40886
0.40891 0.43623
0.95586
0.50340
0.81089
0.83553
0.83580 0.86213
1.83570
0.99846
1.62179
1.77296
1.77331 1.64311
2.87707
1.88499
2.43268
2.84559
2.84630 2.26944
3.15289
2.66326
3.24357
4.07467
4.07600 2.68221
3.42449
3.45741
4.05447
5.47942
5.48210 2.84252
3.55103
4.24400
Tabla 6.12: Comparacin de resultados con otros modelos

Como se observa en la grfica 6.30, los resultados obtenidos por nuestros modelos
tiene un comportamiento muy similar a las frecuencias obtenidas con los otros mtodos
y la diferencia entre ellos es muy pequea. Sin embargo, se observa que el modelo de
99
segundos vecinos es el que tienen un comportamiento un poco diferente a los dems,

lo que nos puede indicar que la interaccin de ms partculas vecinas puede influir en
el comportamiento del sistema, principalmente en su rigidez y flexibilidad.
Figura 6.30: Grfica de comparacin de frecuencias naturales del riser
6.5
Resumen
En este captulo se describi la forma en como fue validado nuestro simulador y las
pruebas que se hicieron para corroborar su funcionamiento. En la primer seccin se
describen los pasos a seguir para hacer un uso adecuado de nuestro simulador.
En la segunda seccin se describe el diseo de las pruebas que se realizaron con el
simulador. Se mostraron los diferentes valores de entrada (masa, elasticidad, longitud
de la estructura, etc.) con los que se realizaron las pruebas.
En la tercera seccin, se mostraron los resultados obtenidos de cada una de las pruebas realizadas para los modelos de primeros y segundos vecinos.
Finalmente en la cuarta seccin se hizo el diseo de una prueba adicional para hacer
una comparacin de los resultados de nuestro modelo con otros modelos, tomando
como base los parmetros de entrada reportados por los otros modelos.
100
6.5. RESUMEN
Captulo 7
Resultados, discusin, y conclusiones
La aventura podr ser loca, pero el
aventurero ha de ser cuerdo.
Gilbert Keith Chesterton
7.1
Resultados
En este trabajo de tesis se abord el problema de las vibraciones que afectan a los riser
durante el proceso de exploracin y explotacin de petrleo en aguas profundas. Para
ayudar a contrarrestar este problema, se propuso la construccin de un simulador que
permita analizar el comportamiento de estas estructuras en diversas condiciones.
Para desarrollar este trabajo de tesis se realiz una investigacin para conocer algunos antecedentes del petrleo, as como la importancia y la influencia que tiene ste
energtico en la economa en mundial y especialmente en nuestro pas. Adems, se
di a conocer el proceso de exploracin y explotacin de petrleo en aguas profundas
y se describieron los diferentes retos que se presentan al realizar este proceso. Especficamente se analiz el problema de las vibraciones que generan las corrientes marinas
sobre las tuberas (conocidas como riser), que se utlizan para la extracin de petrleo en
aguas profundas. Durante esta investigacin, tambin se describeron algunos desarrollos de software que se utilizan en la industria petrolera para controlar y monitorear
los sistemas de riser.
Con la informacin anterior, se realiz la contruccin del modelo matemtico sobre
el cual se basa el simulador. Para el desarrollo de este modelo, se hizo una representacin fsica que nos permiti determinar que el sistema se poda tratar como una
cadena de vibracin monoatomica unidimensional formada por partculas unidas por
resortes.
Una vez que se defini la forma en como se poda tratar el sistema, se desarrollaron
dos modelos matemticos: interaccin a primeros y segundos vecinos. A partir de
estos modelos se obtuvieron dos matrices de movimiento, que nos permiten obtener
los modos de vibracin del sistema y predecir su comportamiento.
101
102
7.2. DISCUSIN
Posteriormente, se desarroll un modelo estocstico de corrientes basado en el mtodo

Monte Carlo. Con este modelo podemos generar diferentes condiciones de corrientes
para probar que nuestros modelos matemticos funcionan adecuadamente.
Con estos modelos, se realiz el desarrollo del simulador. Para esto, se hizo un
diseo general del sistema en donde se describen las diversas herramientas y mtodos
numricos que necesitaban para su desarrollo.
Como resultado final de este trabajo de tesis, desarrollamos un simulador el cual
esta conformado de dos aplicaciones: cmputo de datos y visualizacin. Las cuales se
encargan de realizar los clculos de los datos y la animacin grfica de los resultados
respectivamente.
El simulador se someti a varias pruebas para comprobar su funcionalidad y analizar
el comportamiento del sistema con los dos modelos desarrollados. Los resultados de
las pruebas realizadas se pueden observar en la seccin 6.3 de este docuento.
7.2
Discusin
Durante este trabajo de tesis nos enfrentamos a diversos problemas y para resolverlos
se utilizaron diversas tcnicas, estrategias y herramientas que fueron descritas en este
documento.
Para elegir una manera adecuada de tratar nuestro problema, se realiz una investigacin para conocer el proceso de exploracin y explotacin de petrleo en aguas
profundas, y analizar cmo afectan las corrientes a las tuberas conocidas como riser
que se utilizan para este proceso. Con esta informacin se hizo una reprentacin fsica
del sistema, la cual ayud a determinar que ste tiene un comportamiento similar al
que se present en el problema del puente del milenio descrito en la seccin 2.4.1 [2].
De esta manera, se decidi tratar el problema como una cadena de vibracin y se opt
por la construccin de los dos modelos matemticos que son la base para de nuestro
simulador.
En un principio de esta investigacin, slo se tena contemplado la construccin del
modelo matemtico, sin embargo, se necesitaba tener un mecanismo adicional para
generar diversas condiciones de corrientes que nos permitiera probar nuestro modelo
matemtico. Por esta razn, se desarroll un modelo estocstico basado en el mtodo
Mnte Carlo, que se encargara de generar estas condiciones de corrientes de manera
aleatoria.
Debido a que nuestro modelo de corrientes esta basado en el mtodo Monte Carlo,
fue necesario elegir una distribucin de probabilidad para poder generar datos aleatorios con un cierto comportamiento. A travs de una investigacin que realizamos,
observamos que la distribucin ms usada para simular fenmenos naturales, como
las corrientes marinas, es la distribucin normal.
Por otra parte, fue necesario determinar el tamao de la muestra para cada partcula
y el nmero de veces que se tenan que genarar estos datos. Respecto al vector de frecuencias para cada partcula se determin que fuera de 100 elementos y que cada uno
CAPTULO 7. RESULTADOS, DISCUSIN, Y CONCLUSIONES
103
de estos elementos se deba generar 1000 veces, ya que al realizar pruebas con el sistema se consider que estas cantidades generan condiciones aleatorias adecuadas para
probar el sistema. Sin embargo, para comprobar que estos datos sean correctos y que
esten apegados a la realidad, es necesario hacer pruebas fsicas y comparar los resultados. Otro aspecto que se tom en consideracin para la generacin de frecuencias, es
que deban estar en funcin de nuestro modelo matemtico.
Para poder generar diversas condiciones de corriente a lo largo de la estructura, se
hizo la divisin del riser en segmentos, la cual consiste en dividir el nmero de partculas del riser en un nmero de segmentos proporcionado por el usuario, de tal forma
que se tenga una distribucin uniforme de partculas en todos los segmentos. Esta
divisin, permite determinar los rangos (valores de y ) en los que la distribucin
de probabilidad va a generar los nmeros aleatorios para cada una de las partculas
del segmento. Con esta divisin de segmentos se van a poder generar condiciones de
corriente en 1 dimensin en 2 direcciones distintas.
Con el simulador que se desarroll se pueden disear un gran nmero de pruebas,
variando las condiciones iniciales del sistema (nmero de partculas, segmentos, longitud de riser, etc.). La pruebas que se eligeron para presentar en este documento,
muestran el funcionamiento del sistema con la longitud mxima que puede tener un
riser en aguas profundas. Adems, se hizo la divisin de las partculas en una distancia
relativamente corta y se prob con 3 nmeros de segmentos diferentes. Tomado los resultados obtenidos por las pruebas realizadas, se hizo un anlisis del comportamiento
del sistema en los diversos escenarios.
Respecto a los resultados obtenidos en las pruebas realizadas, en la prueba con 8
segementos y 256 partculas para los modelos de primeros y segundos vecinos (Fig.
6.16 y 6.17), podemos observar que al calcular los movimientos aleatorios de vibracin
de las partculas, se presentan movimientos muy pequeos a lo largo de la estructura.
Adems, se observa que en el modelo de segundos vecinos se tienen ms oscilaciones
(aunque muy pequeas) y en el modelo de primeros vecinos las partculas mantienen
una desplazamiento estable que se va reduciendo conforme trancurre el tiempo. Sin
embargo, al realizar la superposicin de estos movimientos (Fig. 6.18 y 6.19), la estructura se ve claramente afectada, principalmente al inicio de la simulacin ya que
es cuando partculas del sistema tienen un desplazamiento mayor y al realizar la superposicin se tienen movimientos considrables. Como podemos observar el comportamiendo en los dos modelos muy parecido, esto se debe principalmente a que se utilizaron las mismas condiciones iniciales y los mismos valores de frecuencias aleatorias
para realizar los clculos.
Para el caso de la prueba con 16 segmentos y 256 partculas (Fig. 6.21, 6.22, 6.23 y
6.24), pudimos observar que al realizar el clculo de los modos aleatorios de vibracin
y la superposicin de estos, la amplitud de los movimientos de la estructura son ms
pequeos en comparacin con la prueba anterior. Esto se debe principalmente a que al
tener ms segmentos se generan ms condiciones diferentes de fecuencias, por lo que el
movimiento de las partculas de un segmento se puede ver afectado por el movimiento
de las partculas de otro segemento y se reduce su amplitud. Finalmente pudimos
104
7.3. CONCLUSIONES
observar que el comportamiento del sistema es muy parecido al de la prueba anterior,

donde al realizar la superposicin de los movimientos la estructura se ve afectada al
inicio de la simulacin.
En la prueba con 32 segmento y 256 partculas (Fig. 6.26, 6.27, 6.28 y 6.29), se observ
que presentaba una mayor oscilacin de las partculas a lo largo de la simualcin, lo
que puede ser provocado principalmente a los valores de frecuencia que se generaron
para ejecutar la simulacin.
7.3
Conclusiones
En este proyecto de investigacin se desarroll un sistema que permite predecir y

analizar el comportamiento dinmico de los riser durante el proceso de exploracin de
petrleo en aguas profundas. El sistema esta compuesto por un modelo matemtico,
un modelo de corrientes y un sistema de visualizacin de resultados.
Durante la investigacin, nos dimos cuenta que desarrollar un simulador con estas
caractersticas no es una tarea sencilla. Adems de requerir conocimientos tcnicos y
matemticos, es necesario conocer ampliamente el proceso de exploracin de petrleo
en aguas profundas para poder cubrir todos los requerimientos del sistema.
Los resultados obtenidos por el sistema nos ayudaran a identificar problemas en
las estructuras y determinar puntos crticos que pueden ser suceptibles a accidentes.
Adems, nos van permitir generar polticas de desicin que ayuden en la construccin
de los riser, para lograr que sean ms flexibles, resistentes y puedan soportar la fuerza
de la corrientes marinas.
Uno de los aspectos ms importantes durante el desarrollo del simulador, fue la eleccin y la construccin de nuestro modelo matemtico, el cual esta basado en cadenas
de vibracin. Este modelo es la base de nuestro simulador ya que nos permite conocer
o predecir el comportamiento del sistema a partir de ciertas condiciones iniciales.
Por otra parte, el modelo estocstico de corrientes que se desarroll y la divisin
de segmentos, nos permiten experimentar con el sistema en diversas condiciones que
pueden ser muy difciles de generar fisicamente. A pesar de que este modelo de corrientes es simple, sirve como base para poder generar nuevos modelos ms complicados
y que puedan generan condiciones de corrientes diferentes.
Utilizar la biblioteca IMSL nos permiti desarrollar ms rpida y adecudamente el
sistema, gracias a que cuenta con rutinas especializadas y optimizadas que nos ayudaron a resolver el modelo matemtico y el modelo de corrientes.
Con las pruebas que se realizaron, comprobamos que los resultados obtenidos estn
apegados a nuestro modelo matemtico, por lo que podemos concluir que el sistema
esta listo para utilizar datos oceanogrficos reales. Sin embargo, para poder comprobar totalmente el funcionamiento del simulador es neceario realizar una prueba fsica
conocida como Prueba de Tanque, la cual requiere de equipo especializado y es demasiado costosa.
Debido a la situacin actual de la industria petrolera y al inters por buscar nuevas
CAPTULO 7. RESULTADOS, DISCUSIN, Y CONCLUSIONES
105
fuentes proveedoras, se considera que este proyecto surge en un momento adecuado.

Principalmente por que en nuestro pas, PEMEX est buscando fuentes de energa en
aguas profundas y requieren de informacin, conocimientos y tecnologas como la que
se desarroll en este proyecto.
Por ltimo, las aportaciones que se tuvieron con este proyecto de investigacin son:
El desarrollo de un modelo matemtico basado en cadenas de vibracin para
analizar el comportamiento dinmico de los riser.
El desarrollo de un modelo estocstico de corrientes basado en Monte Carlo.
Diseo y desarrollo de un sistema que implementa los dos modelos creados y
genera los resultados del comportamiento del sistema.
Diseo y desarrollo de un sistema de visulizacin para generar imgenes y animaciones a partir de los resultados obtenidos por el sistema anterior.
7.4
Trabajo futuro
Algunas recomendaciones que se pueden realizar para mejorar nuestro simulador son:
Aunque las rutinas de la biblioteca IMSL estan optimizadas para funcionar en
plataformas multicore, es necesario realizar un anlisis para determinar los algoritmos se pueden paralelizar y reducir el tiempo de procesamiento.
Someter el sistema a pruebas con datos oceanogrficos reales, como los que obtiene CICESE a partir de un estudio que esta realizando en el golfo de Mxico.
Adems, realizar la prueba de tanque que nos permita comprobar totalmente el
funcionamiento del simulador.
Realizar un anlisis de los materiales que se utilizan para la construccin de los
riser y determinar los lmites de confianza del sistema.
Desarrollar un modelo ms complejo que nos permita generar corrientes en 2 o
ms dimensiones y que permita simular situaciones como tener dos o mas riser
trabajando juntos.
Mejorar el sistema de visualizacin agregandole ms herramientas que permitan
tener opciones adicionales como: graficar resultados por segmento, generar solo
un tipo de visualizacin, etc.
Mejorar el diseo de la animacin 3D o generar animaciones con otras herramientas como el Quicktime reality.
106
7.4. TRABAJO FUTURO
Bibliografa
[1] Petrleos Mexicanos.
El petrleo. historia, 2009.
http://www.pemex.com. Citada en septiembre 2009.
Disponible en
[2] Francoise Tisseur and Karl Meerbergen. The quadratic eigenvalue problem. Society for Industrial and Applied Mathematics, 43(2):235286, 2001.
[3] William W. Seto. Vibraciones Mecnicas. 1970.
[4] Amilcar Meneses Viveros. Simulacin de paquetes de onda y teora de Weyl. PhD
thesis, Centro de Invetsigacin y de Estudios Avanzados IPN, Mxico, D.F., 2009.
[5] Shelly Bauman and Rockford MI. Intermediate Energy Infobook. NEED, Manassas,
VA 20108 PO Box 10101, need edition, 2008-2009.
[6] Instituto Mexicano del Petrleo.
El petroleo, 2009.
http://www.imp.mx. Citada en septiembre 2009.
Disponible en
[7] Jean A. Bandy Georgius Agricola, Mark Chance Bandy. De Natura Fossilium.
Dover Publications, 2004.
[8] Alphonsus Fagan. An introduction to the petroleum industry. Government of Newfoundland and Labrador, Department of mines and energy, November 1991.
[9] Kjell Aleklett Bengt Sderbergh. Canadas oil sands resources and its future impact on global oil supply. UPPSALA UNIVERSIT, Box 535 SE-751 21 Sweden,
page 5, 2005.
[10] Energy Watch Group. Crude oil the supply outlook. Report to the Energy Watch
Group, EGW-Series No. 3/2007, page 18, Octubre 2007.
[11] Georg Bauer and Rolf Jakoby. A simple prediction model for line of sight coverage
in urban scenarios. In VTC Fall, pages 15. IEEE, 2008.
[12] Joel lvarez de la Borda. Crnoica del petrleo en Mxico de 1863 hasta nuestros das.
Archivo histrico de petrleos mexicanos, primera edicin edition, 2006.
[13] Joel lvarez de la Borda. Los orgenes de la industria petrolera en Mxico 1900-1925.
Archivo histrico de petrleos mexicanos, primera edicin edition, 2005.
107
108
BIBLIOGRAFA
[14] Gerardo Gil Valdivia and Susana Chacn Domnguez. La crsis del petrleo en Mxico. Tecnolgico de Monterrey, Abril 2008.
[15] Abelardo Cant Chapa. Golfo de Mxico. PetroQuiMex, pages 2430, Marzo-Abril
2008.
[16] Amparo Garca Acevedo. Rumbo a aguas profundas; el prximo reto en educacin martima. Fidenatica, (1):812, Enero-Febrero 2008.
[17] Jaeyoung Lee. Design and installation of deepwater petroleum pipelines. In The
World Congress of Korean and Korean Ethnic Scientists and Engineers, Seoul, Korea,
2002. Marine Pipeline Engineer.
[18] Jaeyoung Lee. Introduction to Offshore Pipelines and Risers. Primera edicin edition,
2009.
[19] Mineral Management Service. Deep water, where the energy is. MMS: Securing
Ocean Energy And Economic Value For America, 2004.
[20] Rita Tuub. Deepwater gulf of mexico: Americas offshore energy future. Pipeline
and Gas Journal, pages 4246, Octubre 2008.
[21] J. C. Berger J. C. Cutler. Deepwater Reference Book, volume Volumen 1, 2, 3. septiembre 2000.
[22] Tomas H. Naylor, Josep L. Balintfy, Donald S. Burdick, and Kong Chu. Tcnicas de
Simulacin en Computadoras. 1980.
[23] Mara Cristina Rojas. Burj Al Arab. una superestructura de 7 estrellas. Metal Actual,
Mayo-Junio 2008.
[24] MCS. Flexcom, 2009. Disponible en http://www.mcs.com/site.php?page=15.
Citada en septiembre 2009.
[25] Orcina. Orcaflex, 2009. Disponible en
http://www.orcina.com/SoftwareProducts/OrcaFlex/index.php. Citada en septiembre 2009.
[26] MCS fugro. OPTIMA. Rise Operations Planning and Management. Technical report,
MCS, 2008.
[27] Massachusetts Institute of Technology (MIT. Vortex-induced vibration, 2009.
Disponible en http://web.mit.edu/shear7/shear7.html. Citada en septiembre
2009.
[28] Erasto Victor Vergara Nava. Vibraciones de una Red Diatmica con Interaccin a Segundos Vecinos, 1974. Escuela Superior de Fsica y Matemticas- IPN.
BIBLIOGRAFA
109
[29] Angel Franco Garca. Fsica con ordenador. Curso Interactivo de fsica en Internet, 2009.
Disponible en http://www.sc.ehu.es/sbweb/fisica/. Citada en septiembre 2009.
[30] David Halliday Robert Resnick. Fsica parte I, volume Volumen I. 1977.
[31] Isidro Romero Medina. Vibraciones en Redes Unidimensionales con Interacciones a
Primero Vecinos, 1971. Escuela Superior de Fsica y Matemticas- IPN.
[32] Jesus Ortega Campos. Estudio de Vibraciones en Cadenas Lineales con Iteraccin a
Segundos Vecinos, 1971. Escuela Superior de Fsica y Matemticas- IPN.
[33] Ian Sommerville. Ingeniera de Software. Septima edicion edition, 2005.
[34] Roger S. Presssman. Ingeniera del Software. Un enfoque Prctico. 5 edicin edition,
2002.
[35] Leonard Meirovitch. Fundamentals of Vibrations. 2001.
[36] Clarence W. de Silva. VIBRATIONS Fundamentals and Practice. 2000.
[37] Mark Cartwright. FOURIER METHODS for Mathematicians, scientists and engieners. 1990.
[38] Fabiola Ortega Robles. Ncleo de Monte Carlo y Camino Aleatorio en Ambientes de
Alto Desempeo. PhD thesis, Centro de Invetsigacin y de Estudios Avanzados
IPN, Mxico, D.F., 2008.
[39] Sharon L. Myers Ronald E. Walpole, Raymond H. Myers. Probabilidad y Estadstica
par Ingenieros. Sexta edicin edition, 1999.
[40] Luis Miguel de la Cruz Salas. Visualizacin cientfica. Revista Digital Universitaria
UNAM, 6(12):19, Diciembre 2005.
[41] Visual Numerics. IMSL Math/Library, 1997.
[42] Visual Numerics. IMSL Stat/Library, 1997.
[43] Sparks C.P. Natural vibration analysis of tensioned risers by segmentation
method. Oil Gas Science and Technology, 61(5):647659, 2006.
[44] Sparks C.P. Transverse modal vibrations of vertical tensioned risers. Oil Gas Science and Technology, 57(1):7186, 2002.
110
BIBLIOGRAFA
Apndice A
Instalacin del sistema
Para obtener xito en el mundo, hay
que parecer loco y ser sabio.
Montesquieu
n este anexo se describen los requerimientos necesarios para poder utilizar el simulador. Adems, se describe el procedimiento de instalacin,
as como los pasos para poder compilar el sistema y obtener los archivos
ejecutables a partir del cdigo fuente.
A.1
Requerimientos de Instalacin
Para ejecutar el subsistema de cmputo de datos y vusializacin es necesario contar

con los archivos ejecutables correspondientes a cada uno de ellos. Sin embargo, si el
usuario requiere compilar nuevamente los cdigos, debern cumplir con ciertos requerimientos para llevar a cabo este proceso.
A.1.1
Requerimientos para el subsistema de Cmputo de Datos
Los requerimientos con los que se debe contar para compilar esta aplicacin son:
1. Compilador Absoft Fortran en su versin 11.5 o superior
2. Tener instalada la bilioteca IMSL Fortran Numerical Library
Esta aplicacin se ha desarrollado y probado en una plataforma Mac Os Server 3.1,
sin embargo el compilador Absoft y la biblioteca IMSL esta disponibles para diversas
plataformas (Windows, Linux, etc.), por lo que la compilacin del cdigo no debe de
implicar ningn cambio.
111
112
A.1.2
A.2. TUTORIAL DE COMPILACIN
Requerimientos para el subsistema de Visualizacin
Lo requerimientos para poder compilar esta aplicacin son:

1. Compilador Ansi C (gcc) version 4.2.4
2. Se deber contar con las libreras de Qt correspondientes a la versin 4.0
3. Finalmente se debe contar con la biblioteca OpenGl (Open Graphics Library)
Esta aplicacin se ha desarrollado y probado en una plataforma Linux Ubuntu 8.10
con kernel 2.6.24
A.2
Tutorial de compilacin
Al igual que en los requerimientos, el proceso de compilacin e instalacin de cada una

de los subsistemas es distinto, por lo que se dar una explicacin por separado para
cada uno de ellos.
A.2.1
Subsistema de Cmputo de Datos
Este subsistema se compone de los siguientes archivos:

SistComputo.f90: Contiene las funciones para generar el clculo de los resultados
SistComputo: Archivo ejecutable
compila: Archivo que contiene el script para compilar el cdigo fuente
Para poder compilar esta aplicacin se tiene que hacer desde lnea de comando mediante la siguiente instruccin:
f90 SistComputo.f90 -o SistComputo -m64 $LINK_f90
donde la instruccin -m64 indica que se esta compilando en un procesador de 64 bits y
$LINK_f90 indica que hacemos uso de las rutinas de la biblioteca IMSL.
Esta instruccin se encuentra en un scritp llamado compila, por lo que para compilar
solo es necesario ejecutar la siguiente instruccin:
./compila
Esta instruccin se encarga de generar nuestro archivo ejecutable el cual tendr el
nombre del SistComputo y para ejecutar este archivo se realiza mediante el siguiente
comando:
./SistComputo
APNDICE A. INSTALACIN DEL SISTEMA
A.2.2
113
Subsistema de Visualizacin
Para poder compilar el sistema de visualizacin, es necesario tener todos los archivos
correspondientes a la aplicacin (cdigo fuente, archivos de configuracin, etc.) en una
carpeta y dentro de sta un directorio con las imgenes que se utlizan en la interfaz
grfica.
Para comenzar a compilar el proyecto, es necesario generar el archivo de configuracin (.pro) y el archivo Makefile que nos permite generar el archivo ejecutable. Para
generar el archivo de configuracin debemos posicionarnos en la carpeta donde tenemos los archivos fuentes e introducir el comando:
qmake-qt4 -project
Esto nos crear un archivo con el nombre de la carpeta y extensin .pro (i.e. Simulador.pro). Con este archivo vamos a generar nuestro Makefile, pero antes de generarlo
debido a que en nuestra aplicacin usamos rutinas de la biblioteca OpenGL, es necesario agregar una lnea al archivo .pro indicndole que vamos hacer uso de la rutinas
de OpenGL. Para esto es necesario abrir el archivo y agregarle la siguiente instruccin
al final:
QT += opengl
Posteriormente podemos generar nuestro archivo Makefile mediante la siguiente instruccin:
qmake-qt4
Finalmete ejecutamos nuestro archivo Makefile y obtendremos como resultado el
archivo ejecutable para poder utilizar el simulador.
./Makefile
Si no se ha realizado ninguna modificacin al sistema y se requiere compilar nuevamente el cdigo, no es necesario repetir todo el proceso, slo se tendr que ejecutar
nuevamente el archivo Makefile que generamos anteriormente.

Tes Is Julian Guzman

Cargado por

Copyright:

Formatos disponibles

Tes Is Julian Guzman

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Tes Is Julian Guzman

Cargado por

Copyright:

Formatos disponibles

CENTRO DE INVESTIGACIN Y DE ESTUDIOS AVANZADOS

DEL INSTITUTO POLITCNICO NACIONAL

Modelo de Vibracin y Simulacin de Comportamiento

Tesis que presenta:

Julin Edgardo Guzmn Corts

Para obtener el grado de

Maestro en Ciencias en Computacin

Director de la Tesis: Dr. Sergio Vctor Chapa Vergara

Codirector de la Tesis: Dr. Amilcar Meneses Viveros

A lo largo de este trabajo de investigacin y durante mis estudios de maestra, he sido

Desarrollo e implementacin del sistema

Resultados de las pruebas . . . . . . . . . . . .

Resultados, discusin, y conclusiones

A Instalacin del sistema

Proceso de formacin del petrleo . . . . . . . . . . . . . . .

Representacin del modelo fsico. Tubera sumergida en agua dividida en n

Divisin del simulador en dos subsistemas . . . . . . . . . . . . . . . . . . .

Representacin general del sistema . . . . . . . . . . . . . . . . . . . . . . . .

Arquitectura general del sistema . . . . . . . . . . . . . . . . . .

Pantalla de inicio de la aplicacin de cmputo de datos . . . . . . . . . . . . .

Visualizacin de la superposicin con 16 segmentos para segundos vecinos . .

Principales pases productores de petrleo en aguas profundas . . . . . . . . .

Datos de entrada para la ejecucin de la prueba de 8 segmentos . . . . . . . .

Genera Matriz de Movimiento . . . . .

ctualmente existen diversas formas y fuentes proveedoras de energa, y el petrleo

El petrleo en aguas profundas

Imagen tomada de sitio http://www.telefonica.net/web2/marcillaot/

2.1. EL PETRLEO EN AGUAS PROFUNDAS

Figura 2.1: Proceso de formacin del petrleo

Figura 2.2: Pozo de alquitrn

Imagen tomada de http://blog.webosaurs.com/2009/10/page/2/

utilizaban en lmparas y en las mquinas de la revolucin industrial. No obstante, la

Figura 2.3: Principales usos del petrleo

Imagen tomada del sitio http://atinachile.bligoo.com/content/view/22576

2.1. EL PETRLEO EN AGUAS PROFUNDAS

Debido a que este energtico es un recurso no renovable y a la dificultad de sustituir

2. Exploracin en aguas someras: Esta etapa de exploracin se hizo a finales del

Exploracin y explotacin de petrleo en aguas profundas

Imagen tomada del sitio http://www.pemex.com [1]

2.1. EL PETRLEO EN AGUAS PROFUNDAS

Figura 2.4: Regiones de aguas profundas

Tabla 2.1: Principales pases productores de petrleo en aguas profundas

Figura 2.5: Superficie de las Aguas Profundas en el Golfo de Mxico

Para el MMS de Estados Unidos, el Golfo de Mxico se ha convertido en su mayor

Imagen tomada del sitio http://www.pemex.com [1]

2.1. EL PETRLEO EN AGUAS PROFUNDAS

Retos y dificultades en la exploracin y explotacin en aguas

Figura 2.6: Procesos para desarrollar un proyecto en Aguas Profundas

Imagen tomada del sitio http://www.pemex.com [1]

Desarrollo de campos petroleros en alta mar

Como se ha mencionado anteriormente, el desarrollo de campos petrolferos es un

2.2. DESARROLLO DE CAMPOS PETROLEROS EN ALTA MAR

Figura 2.7: Elementos para el desarrollo de campos en alta mar

Figura 2.8: Clasificacin de sistemas submarinos

Lneas de flujo, tuberas y sistemas riser

Las lneas de flujo tambin se denominan lnea de produccin o de importacin y se

2.2. DESARROLLO DE CAMPOS PETROLEROS EN ALTA MAR

Figura 2.9: Lnea de flujo, tuberas y sistemas riser