Gricelcy Perez.pdf - Departamento de Ciencias de la Tierra ...

UNIVERSIDAD SIMÓN BOLÍVAR
DECANATO DE ESTUDIOS PROFESIONALES
COORDINACIÓN DE INGENIERÍA GEOFÍSICA
MODELADO 3D DE VELOCIDADES INTERVALICAS DE ONDA P
A PARTIR DE REGISTROS DE POZO
CAMPO SANTA BARBARA
Por
Gricelcy Nairam Pérez Flores
INFORME FINAL DE CURSOS EN COOPERACIÓN
Presentado ante la Ilustre Universidad Simón Bolívar
como requisito parcial para optar al título de
Ingeniero Geofísico
Sartenejas, Noviembre 2008

UNIVERSIDAD SIMÓN BOLÍVAR
DECANATO DE ESTUDIOS PROFESIONALES
COORDINACIÓN DE INGENIERÍA GEOFÍSICA
MODELADO 3D DE VELOCIDADES INTERVALICAS DE ONDA P
A PARTIR DE REGISTROS DE POZO
CAMPO SANTA BARBARA
Por
Gricelcy Nairam Pérez Flores
Realizado con la Asesoría de:
Tutora Académica: Dra. Milagrosa Aldana
Tutor Industrial: Ing. Manuel Bolívar
INFORME FINAL DE CURSOS EN COOPERACIÓN
Presentado ante la Ilustre Universidad Simón Bolívar
como requisito parcial para optar al título de
Ingeniero Geofísico
Sartenejas, Noviembre 2008

MODELADO 3D DE VELOCIDADES INTERVALICAS DE ONDA P
A PARTIR DE REGISTROS DE POZO
CAMPO SANTA BARBARA
Por
Gricelcy Nairam Pérez Flores
RESÚMEN
La Cuenca Oriental de Venezuela, es la segunda cuenca petrolífera más importante de
Venezuela. Además es una de las provincias geológicas más complejas; ésta complejidad y
la poca calidad de la imagen sísmica, ha obstaculizado la interpretación sísmica de la zona,
y los modelos estructurales existentes han honrado medianamente la geología de la zona.
Por esto, el objetivo de este trabajo es obtener un modelo de velocidades interválicas del
Campo Santa Bárbara, para mejorar la calidad de la imagen sísmica mediante el futuro
reprocesamiento de los datos sísmicos del área. Se presenta una metodología que muestra el
control y caracterización de los registros sónicos y “Check Shot” existentes en el área,
identificando las tendencias de las velocidades por estrato geológico en cada pozo, para
obtener así ecuaciones de velocidad en función de la profundidad para cada región definida
las cuales son: Somero, Bloque Alóctono, Carapita y Yacimiento. Cabe destacar que la
construcción de un buen modelo de velocidades interválicas ayudará con la obtención de
buenas imágenes sísmicas, paso esencial para la interpretación y la generación de modelos
estructurales.

iv
DEDICATORIA
A Dios
A Francy, mi mamá
A Julio, mi papá
Y a mis hermanos, Amanda y Octavio
Esto es de ustedes tanto como mío.

AGRADECIMIENTOS
Agradezco a mis padres, Francy y Julio por demostrarme su amor día a día y por ser un
apoyo incondicional, sin importar cuán lejos estemos sé que siempre me tienen en su mente
y corazón. Quien soy hoy es gracias a ustedes. Los amo.
Gracias a Amanda y a Octavio por su ayuda. Sin ustedes no hubiese terminado a tiempo,
son muy importantes para mí, su amor y ayuda me fortalece. Recuerden siempre que la
batalla de la vida no se gana por ser más rápido o más fuerte, sino que tarde o temprano,
quien gana es aquel que siempre pensó en ganar. Los adoro.
A mi tutor Manuel Bolívar, gracias por abrirme las puertas de PDVSA, por brindarme sus
conocimientos y por guiarme en la elaboración de este trabajo. Gracias a Mario Díaz y
Ariannella Levanti por haberme enseñando a usar GOCAD.
Agradezco la profesora Milagrosa Aldana, por acceder a ser mi tutora académica y
ayudarme en todo lo necesario con respecto a la elaboración del informe final.
Gracias a todos las personas que conocí durante mi estadía en las instalaciones de
PDVSA Exploración y Producción Oriente. A los ingenieros de la Gerencia de Operaciones
Geofísicas y a los del Centro de Procesamiento Geofísico.
A mis amigos Karla, Ender, Cynthia y Eugenio por los buenos ratos que pasamos durante
los almuerzos y los últimos meses de mi estadía en Puerto La Cruz.
Gracias a mis amigos Luis, Miguel, Iraivel, Picardo, Mafe, Adriana, a todos los
muchachos de la Proveeduría y a todos los compañeros Geofísicos que conocí y con los
cuales compartí durante toda la carrera.
Stephane, gracias por acompañarme de corazón en todo esto. Por estar presente siempre
que te necesite. Gracias por tu amor, comprensión y amistad. ¡Te Amo!
Finalmente, a todas las personas que de una u otra manera contribuyeron a que llegara
hasta aquí. ¡GRACIAS!
v

INDICE GENERAL
RESÚMEN ............................................................................................................................ iii
DEDICATORIA .................................................................................................................... iv
AGRADECIMIENTOS ......................................................................................................... v
INDICE GENERAL ............................................................................................................. vi
INDICE DE FIGURAS ...................................................................................................... viii
INDICE DE TABLAS .......................................................................................................... xi
INTRODUCCIÓN ................................................................................................................. 1
CAPÍTULO 1 ......................................................................................................................... 4
MARCO GEOGRÁFICO Y GEOLOGICO ...................................................................... 4
1.1. Ubicación del área de estudio ............................................................................. 4
1.2. Geología Regional .............................................................................................. 5
1.2.1. Cuenca Oriental de Venezuela ....................................................................... 5
1.3. Geología Local ................................................................................................... 7
1.3.1. Subcuenca de Maturín .................................................................................... 7
1.3.2. Evolución Téctono-Estratigráfica de la Sub-Cuenca de Maturín ................. 10
1.3.3. Campo Santa-Bárbara/Pirital ........................................................................ 16
1.3.4. Campo Tácata ............................................................................................... 19
CAPÍTULO 2 ....................................................................................................................... 22
MARCO TEÓRICO ........................................................................................................... 22
2.1. Métodos Sísmicos ............................................................................................. 22
2.2. Velocidades Sísmicas ....................................................................................... 22
2.2.1. Velocidad interválica Vint ............................................................................. 23
2.2.2. Velocidad promedio (V) ............................................................................... 23
2.2.3. Velocidad promedio (Vrsm) ......................................................................... 24
2.3. Migración sísmica ............................................................................................. 24
2.4. Modelo de velocidades ..................................................................................... 25
vi

2.5. La herramienta sónica ....................................................................................... 25
2.6. La herramienta sónica compensada .................................................................. 27
2.7. Disparos de verificación sísmica ...................................................................... 29
CAPÍTULO 3 ....................................................................................................................... 32
METODOLOGIA ............................................................................................................... 32
3.1. Recopilación de datos de registros de pozos .................................................... 32
3.2. Control de Calidad de Registros de Pozos ....................................................... 33
3.4. Modelo de velocidades ..................................................................................... 44
3.5. Construcción de Voxet o malla 3D .................................................................. 47
3.6. Carga y modelado de superficies ...................................................................... 48
3.7. Propagación de velocidades ............................................................................. 49
CAPÍTULO 4 ....................................................................................................................... 50
RESULTADOS Y ANALISIS DE RESULTADOS ......................................................... 50
4.1 Estudio y Control de calidad de Registros ........................................................... 50
4.2 Región Somera ..................................................................................................... 54
4.3 Región Bloque Alóctono ...................................................................................... 62
4.4 Región Carapita .................................................................................................... 69
4.5 Región Yacimiento ............................................................................................... 76
4.7 Modelo de Velocidades ........................................................................................ 83
4.8 Histogramas de Velocidades (Modelo Final) ..................................................... 104
CONCLUSIONES ............................................................................................................. 106
RECOMENDACIONES .................................................................................................... 108
vii

INDICE DE FIGURAS
Figura 1-1 Ubicación geográfica del área de estudio
viii
............................................................ 4
Figura 1-2 Cuenca Oriental de Venezuela ............................................................................. 5
Figura 1-3 Marco tectónico y geográfico de mayores campos de la Cuenca Oriental de
Venezuela
........................................................................................................... 7
Figura 1-4 Corte geológico conceptual noreste-sureste desde la plataforma Margarita-Los
Testigos hasta el río Orinoco . ............................................................................. 8
Figura 1-5 Mapa Estructural de la Serranía del Interior y el frente de corrimientos . ............ 9
Figura 1-6 Columna estratigráfica de la Cuenca Oriental de Venezuela. ............................ 14
Figura 1-7 Corte geológico a partir de la interpretación sísmica, Bloque Este .................... 16
Figura 1-8 Cortes geológicos a partir de la interpretación sísmica ...................................... 17
Figura 1-9 Desarrollo de la sísmica, campo Santa Bárbara entre los años 1986 a 1999 ...... 18
Figura 1-10 Mapa de ubicación de la información sísmica ................................................. 18
Figura 1-11 Rasgos estructurales principales definidos en la estructura de Tácata
............ 20
Figura 1-12 Levantamientos sísmicos 3D del Volumen Integrado Tacat_2000. ................ 21
Figura 2-1 Herramienta sónica
............................................................................................ 26
Figura 2-2 Diagrama de viaje de las ondas .......................................................................... 27
Figura 2-3 Problemas con las herramientas de un solo transmisor. ..................................... 28
Figura 2-4 Secuencia de las mediciones sónico compensado
............................................. 29
Figura 3-1 Ubicación geográfica pozos con registros sónicos en el área ............................. 33
Figura 3-2 Registro Check Shot pozo PIC-21 con valor anómalo ....................................... 35
Figura 3-3 Registro Sónico de pozo PIC-8 ........................................................................... 37
Figura 3-4 Registro Sónico de pozo PIC-08 con medianas de 100, 200 y 300 datos ........... 39
Figura 3-5 Registro Sónico de pozo PIC-28 con mediana de 300 ........................................ 40
Figura 3-6 Registro Sónico y “Check Shot” de pozo PIC-28 ............................................. 42
Figura 3-7 Registro Sónico y “Check Shot” pozo SBC-15 .................................................. 43
Figura 3-8 Imagen de pozo desviado en GOCAD ................................................................ 45
Figura 3-9 Imagen con propiedades de velocidad ................................................................ 46
Figura 3-10 Imagen del Voxet Campo Santa Bárbara/Pirital y Tácata con “celdas” ........... 47

Figura 4-1 Paquetes litológicos identificados en el área norte (a) y área sur (b) del campo
Santa Bárbara/Pirital. ......................................................................................... 51
Figura 4-2 Regiones reconocidas en registro sónico y columna estratigráfica asociada para
el área sur del campo Santa Bárbara/Pirital ....................................................... 52
Figura 4-3 Regiones reconocidas en registro sónico y columna estratigráfica asociada para
el área sur del campo Santa Bárbara/Pirital ....................................................... 53
Figura 4-4 Registros Sónicos región somera. Notar como no existe correspondencia entre
todas las curvas de velocidad ............................................................................. 54
Figura 4-5 Velocidad interválica y curva de ajuste lineal para la Serie I ............................. 56
Figura 4-6 Velocidad interválica y curva de ajuste lineal para la Serie II ............................ 56
Figura 4-7 Velocidad interválica y curva de ajuste lineal para la Serie III .......................... 57
Figura 4-8 Velocidad interválica y curva de ajuste lineal para la Serie IV .......................... 57
Figura 4-9 Velocidad interválica y curva de ajuste lineal para la Serie V ........................... 58
Figura 4-10. Velocidad interválica y curva de ajuste lineal para la Serie VI ....................... 58
Figura 4-11. Comparación de curvas de registro Sónico (en azul) y Check Shot (puntos
rojos) en región Somera ..................................................................................... 60
Figura 4-12 Velocidad interválica pozo PIC-12 con curva de ajuste lineal ......................... 62
Figura 4-13 Pozos con “Check Shot” Bloque Alóctono ....................................................... 63
Figura 4-14 Registro Sónico Bloque Alóctono mostrando picos abruptos y tendencia
general ............................................................................................................... 64
Figura 4-15 Curvas de velocidad interválica con ajuste lineal para la series I, II, III y IV de
la región Bloque Alóctono. ................................................................................ 65
Figura 4-16 Velocidades interválicas con ajuste lineal para región Bloque Alóctono. ........ 67
Figura 4-17. Comparación registros sónicos y “Check Shot” región Blque Alóctono. ...... 68
Figura 4-18 Pozos con Check Shot Tendencia Carapita ...................................................... 69
Figura 4-19 Curvas de velocidad interválica mostrando tendencia general y alta correlación
entre los datos .................................................................................................... 70
Figura 4-20 Velocidad interválica pozo PIC-12 con curva de ajuste lineal ......................... 71
Figura 4-21 Velocidad interválica y curva de ajuste lineal para la Serie I ........................... 72
Figura 4-22. Velocidad interválica y curva de ajuste lineal para la Serie II ......................... 73
Figura 4-23. Velocidad interválica y curva de ajuste lineal para la Serie III ....................... 73
ix

Figura 4-24. Velocidad interválica y curva de ajuste lineal para la Serie IV ....................... 74
Figura 4-25. Velocidad interválica y curva de ajuste lineal para la Serie V ........................ 74
Figura 4-26. Registros Sónicos con tendencia en región Carapita ...................................... 75
Figura 4-27. Curvas de velocidad tipo de la región Yacimiento .......................................... 76
Figura 4-28. Velocidades Región Yacimiento .................................................................... 77
Figura 4-29. Velocidades interválicas a partir de registro Check Shot del pozo PIC-25,
región Yacimiento de color verde. ..................................................................... 78
Figura 4-30. Ajuste lineal para las velocidades derivadas del registro Check Shot Región
Yacimiento ......................................................................................................... 79
Figura 4-31. Histograma de velocidades interválicas de la región Somera .......................... 80
Figura 4-32. Histograma de velocidades interválicas de la región Bloque Alóctono .......... 81
Figura 4-33. Histograma de velocidades interválicas de la región Carapita ........................ 81
Figura 4-34. Histograma de velocidades interválicas de la región Yacimiento ................... 82
Figura 4-35. Pozos del área de estudio visto en ambiente GOCAD . ................................... 83
Figura 4-36. Voxet de trabajo área extensa, definido por los pozos del campo Santa
Bárbara/Pirital y Tácata ..................................................................................... 85
Figura 4-37. Vista en planta de Voxet “área reducida” con la ubicación de pozos
pertenecientes al volumen de interés. ................................................................ 86
Figura 4-38. Comparación de velocidades derivadas de registros sónicos del pozo SBC-35,
curva típica de campo Santa Bárbara, y del pozo TAG-11, curva típica de
campo Tácata. .................................................................................................... 87
Figura 4-39. Curva de velocidades interválicas con marcadores del pozo TAG-11 ............ 88
Figura 4-40 Superficies generadas en GOCAD dentro del VOXET “área reducida” .......... 89
Figura 4-41 Regiones definidas dentro del VOXET “área reducida” .................................. 90
Figura 4-42 Modelo de velocidades en VOXET Extenso. Cortado por las superficies ....... 91
Figura 4-43. Modelos 3D de Velocidades Interválicas Campo Santa Bárbara .................... 95
Figura 4-44. Modelo 3D de Velocidades Interválica Campo Santa Bárbara. Velocidades
para la Región Somera ....................................................................................... 96
Figura 4-45. Modelo 3D de Velocidades Interválica Campo Santa Bárbara. Velocidades
para la Región Bloque Alóctono. ...................................................................... 97
x

Figura 4-46. Modelo 3Dde Velocidades Interválica Campo Santa Bárbara. Velocidades
para la Región Carapita. .................................................................................... 98
Figura 4-47. Modelo 3D de Velocidades Interválica Campo Santa Bárbara. Velocidades
para la Región Yacimiento. ............................................................................... 99
Figura 4-48. Comparación del modelo de velocidades con la interpretación sísmica I.
Presencia de impedancias importantes entre la base de la región Bloque
Alóctono y el tope de la región Carapita ........................................................ 101
Figura 4-49. Comparación del modelo de velocidades con la interpretación sísmica II.
Presencia de impedancias importantes entre la base de la región Bloque
Alóctono y el tope de la región Carapita ........................................................ 102
Figura 4-50. Influencia del alóctono de Pirital sobre la calidad de la sísmica ................... 103
Figura 4-51. Histograma de velocidades del modelo final. (a) Histograma región Somero,
(b) Histograma región Bloque Alóctono, (c) Histograma región Carapita. .... 104
Figura 4-52 Histograma de velocidades para el modelo final de velocidades interválicas
generado en el Voxet “área reducida” ............................................................. 105
INDICE DE TABLAS
Tabla 1-1. Descripción breve de formaciones geológicas de interés. .................................. 15
Tabla 3-1. Check Shot SBC12 .............................................................................................. 34
Tabla 4-1. Pozos con tendencia Somero ............................................................................... 55
Tabla 4-2 Pozos con registro Check Shot región Somero .................................................... 59
Tabla 4-3. Registro “Check Shot” región Bloque Alóctono ............................................... 63
Tabla 4-4.Pozos con tendencia en región Bloque Alóctono ................................................. 65
Tabla 4-5. Pozos con tendencia en la región Carapita .......................................................... 72
Tabla 4-6. Dimensiones de VOXET extenso ....................................................................... 84
Tabla 4-7. Dimensiones VOXET área reducida ................................................................... 85
xi

INTRODUCCIÓN
Día a día se hace más difícil para la industria petrolera descubrir nuevos yacimientos
petrolíferos, pues a medida que pasan los años se descubren yacimientos de estructura
compleja y se convierte en un reto constante el estudio integrado de las geociencias
(Bastidas, 2008). En la última década ha habido un constante desarrollo tecnológico que ha
acarreado la optimización de métodos de caracterización del subsuelo.
El modelado 3D para procesar y visualizar la información geocientífica es un hecho. Sin
embargo, se enfrenta a una falta de información debido a que el subsuelo es investigado por
lo general de forma irregular mediante registros de pozo (puntos de muestreo),
afloramientos, secciones sísmicas, etc. Esta incertidumbre en cuanto a la forma real implica
que las aplicaciones clásicas de diseño y visualización geocientífica podrían ser
irrelevantes. Y más aun cuando la situación puede ser compleja por la presencia de
discontinuidades tales fallas, fracturas, corrimientos, etc. Es entonces cuando entra en
juego el modelado de objetos naturales en 3D. Mediante el uso de un software de modelado
con herramientas geoestadísticas las propiedades físicas del subsuelo son estimadas a lo
largo y ancho de un volumen y hace posible reducir la incertidumbre que los métodos de
prospección geocientífica aportan (Royer, 2004)
El campo Santa Bárbara desde el punto de vista de producción petrolera es un campo muy
maduro que actualmente está en desarrollo con reservas significativas, probadas y
remanentes de gas y condensado. La complejidad estructural del área y las variaciones en la
calidad de los datos están asociadas a un cierto nivel de incertidumbre, que genera la
necesidad de actualizar continuamente los modelos disponibles y procurar la integración de
los datos sísmicos, petrofísicos, geoquímicos, sedimentológicos y estratigráficos, en pro de
una mejor y actualizada definición del sistema petrolífero.
La unidad estructural que resalta del campo Santa Bárbara/Pirital es el bloque alóctono de
Pirital que se encuentra encima del anticlinal de rampa de El Furrial/Carito/Santa-Bárbara.
Está constituido por un anticlinal de rampa, con un desplazamiento horizontal hacia el Sur
que varía entre 20 y 30 Km y con un salto vertical de 4 a 5 Km. Este corrimiento tiene un
1

umbo oblicuo al rumbo general de los demás corrimientos. Su existencia provoca efectos
sobre la calidad de los datos sísmicos. Debido a su gran espesor, la energía de las ondas
sísmicas pierde intensidad, bien sea por absorción o por dispersión. También, el contraste
importante de impedancia entre las formaciones del bloque alóctono y las lutitas de la
Formación Carapita que se encuentran por debajo, hace que las ondas sísmicas se reflejen
sobre el plano de falla constituido por el corrimiento de Pirital, produciendo reflexiones
múltiples que perturban la señal en los niveles más profundos.
Además, el campo no cuenta con un levantamiento sísmico único que cubra toda su
extensión, siendo la información actual un resultado de mezclar varios levantamientos sin
llegar aun a procesar un volumen de datos “mezclados” de manera ajustada a los
requerimientos de la caracterización, por lo que se tienen planes de hacer un
reprocesamiento de la información sísmica. Es por esto que en el presente trabajo se
pretende elaborar un modelo de velocidades que se ajuste a las características, tanto
estratigráficas como estructurales del campo con el fin de contar con un buen
entendimiento de las velocidades de propagación de ondas compresionales, para la
aplicación en el procesamiento de migración pre apilamiento. Para ello se integraron los
datos provenientes de perfiles sísmicos “Check Shot”, información geológica, así como
información estructural del área proporcionada por el modelo estructural vigente del campo
Santa Bárbara/Pirital para poder obtener un modelo cónsono con la estratigrafía y estructura
del volumen en estudio.
Para lograr la construcción del modelo 3D de velocidades interválicas de onda
compresional, a partir de los registros de pozo, sónico y “Check Shot”, para el campo Santa
Bárbara, al norte del estado Monagas se procederá a:
1. Validar y compilar información de los registros sónicos y “Check Shot” de los
pozos disponibles en el área de estudio.
2. Calcular y construir graficas de velocidades interválicas de propagación de onda
compresional a partir de registros de pozo sónico y “Check Shot”
3. Identificar las tendencias verticales de los registros de pozo sónico y “Check Shot”.
2

4. Modelar un marco estructural de capas, horizontes y topes de formación para crear
regiones de interés en el volumen estudiado.
5. Efectuar la propagación de las velocidades dentro de las regiones previamente
identificadas en el estudio haciendo uso de las estrategias y herramientas de
modelado geoestadístico.
Actualmente el campo Santa Bárbara/Pirital está siendo explotado por la empresa
Petróleos de Venezuela S.A. que es la corporación estatal de la República Bolivariana de
Venezuela que se encarga de la exploración, producción, manufactura, transporte y
mercadeo de los hidrocarburos. Este trabajo fue realizado en el edificio sede de PDVSA –
División Oriente en la Dirección Ejecutiva de Exploración, específicamente en la Gerencia
de Operaciones Corporativas de Geofísica y Geodesia en colaboración a la Gerencia de
Estudios Integrados de Yacimiento en el marco de la pasantía larga Junio-Noviembre 2008.
3

1.1. Ubicación del área de estudio
CAPÍTULO 1
MARCO GEOGRÁFICO Y GEOLOGICO
El área de estudio se encuentra ubicada al noreste del estado Monagas, dentro de la
Cuenca Oriental de Venezuela, en la Sub Cuenca de Maturín. Desde el punto de vista
operacional incluye los campos Santa Bárbara/Pirital y Tácata (Figura 1-1)
Figura 1-1 Ubicación geográfica del área de estudio (Modificado de WEC, 1997)

1.2. Geología Regional
1.2.1. Cuenca Oriental de Venezuela
La Cuenca Oriental es la segunda cuenca petrolífera más importante de Venezuela.
Consiste en una depresión estructural ubicada en la región Centro-Este del país.
Específicamente se encuentra entre los 8° y 11° de latitud Norte y los 61° y 66° de longitud
Oeste (Figura 1-2). La depresión topográfica y estructural que forma la Cuenca Oriental de
Venezuela tiene una longitud aproximada de 800 Km en sentido este-oeste, una anchura
promedio de 200 Km de norte a sur y un área total aproximada de 160.000 km2 que cubre
los Estados Guarico, Anzoátegui, Monagas, Delta Amacuro y una extensión menor del
Estado Sucre y se prolonga por el golfo de Paria hasta Trinidad. (Léxico estratigráfico de
Venezuela, 1997)
Debido a sus características sedimentológicas, tectónicas y estratigráficas la cuenca ha
sido subdividida operacionalmente en dos subcuencas: la de Guárico y la de Maturín
(González de Juana et al., 1980).
Figura 1-2 Cuenca Oriental de Venezuela (Modificado del WEC, 1997)
5

Los límites de la Cuenca Oriental de Venezuela son:
• Al norte, la Cordillera de la Costa/Villa de Cura y los cinturones ígneos y metamórficos
de Araya y Paria.
• Al sur, el escudo de Guayana que está compuesto por diferentes litologías de rocas de
edad Pre-cámbrica. El escudo de Guayana ha sido la mayor fuente de sedimentos de la
historia geológica de la cuenca.
• Al este, el límite no está bien definido pero se considera que lo constituye la corteza
oceánica del Atlántico Ecuatorial.
• Al oeste, el Arco del Baúl, compuesto por Rocas ígneas de edad Paleozoico
Debido a los esfuerzos en dirección aproximada norte-sur que actuaron durante el
Terciario Inferior el flanco norte de la Cuenca Oriental de Venezuela está caracterizado por
plagamiento y fallamiento de tipo compresivo. Presenta un sistema de fallas de corrimiento
buzando al norte con una dirección general N60°E y N90°E como lo son los corrimientos
de Anaco y Pirital, el Furrial y las fallas de Onado, Guanaco y San Juan; fallas
trasncurrentes dextrales de rumbo general N80°O que desplazan a las primeras y están
representadas por las fallas de Urica, San Francisco, Bohordal y los Bajos (Figura 1-3).
En General la cuenca tiene un mas de 8 Km de sedimentos terciarios. El flanco Sur buza
hacia el norte y el flanco norte está limitado por los cinturones de cabalgamientos y
plegamientos de la Serranía del Interior, el Golfo de Paria y Trinidad. La continuación
hacia el mar de esta cuenca es la plataforma costa afuera del Orinoco (Plataforma Deltana),
la cual se une hacia el Sur con el Margen pasivo Atlántico de Guayana. (Erlich y Barrett,
1992).
6

Figura 1-3 Marco tectónico y geográfico de mayores campos de la Cuenca Oriental de
1.3. Geología Local
1.3.1. Subcuenca de Maturín
Venezuela (Modificado de Parnaud, 1995)
La Subcuenca de Maturín está ubicada al Este de la Cuenca Oriental y constituye la
principal unidad petrolífera de la Cuenca Oriental, con reservas en el orden de 25
MMMBls. Ha sido objeto de la exploración y producción de hidrocarburos desde 1930
aproximadamente (Yoris y Ostos, 1997). Además podría afirmarse que la deformación
estructural y los acuñamientos de las unidades estratigráficas hacia el sur definen dos
dominios operacionales: uno al norte del Corrimiento de Pirital y otro al sur.
Es una Cuenca alargada de dirección N50°E paralela a la Serranía del interior, asimétrica,
que posee un flanco Sur pasivo apoyado en el escudo de Guayana y un flanco norte activo.
Además la región ha recibido sedimentos casi de forma ininterrumpida desde el Cretácico
7

Inferior al Reciente. La Sub-cuenca se formo a partir del Eoceno, encima de una cuenca de
margen pasivo, en ella se presentan cambios de facies y discordancias. La Cuenca es de
carácter marino dominante con algunas turbiditas en el flanco norte, depósitos
sedimentarios sintectónicos y continentales (Figura 1-4)
Figura 1-4 Corte geológico conceptual noreste-sureste desde la plataforma Margarita-Los
Testigos hasta el río Orinoco (Tomado de Yoris y Ostos, 1997).
La transcolisión de la placa Caribe con la Sur Americana produjo tres familias de fallas
principales en la Cuenca Oriental de Venezuela: fallas transcurrentes dextrales de rumbo
Este-Oeste (falla de El Pilar), fallas transcurrentes dextrales de rumbo N60°O-N80ºO
(como las fallas principales de Urica y de San Francisco) y corrimientos de rumbo N60°E-
N90°E (como los corrimientos de Furrial o de Pirital). Las fallas transcurrentes, de
direcciones sub paralelas al esfuerzo principal constituyen los rieles y rampas laterales de
los corrimientos que tienen direcciones sub perpendiculares al esfuerzo principal. (Bosset J.
y Gou Y. 2004)
En el bloque de Bergantín, al Sur de la Serranía del Interior, y soterrado por los depósitos
del Plio – Pleistoceno, se ubica el corrimiento de Pirital, cuyo rumbo es en general oblicuo
con respecto a los otros corrimientos de la cordillera. Esta falla presenta un desplazamiento
8

mayor hacia el Oeste, cerca de la falla de Urica (Figura 1-5). El corrimiento de Pirital
involucra una sección sumamente espesa que parece incluir unidades Pre-Cretácicas,
desconocidas en la Serranía (Parnaud et al., 1991)
BLOQUE DE
SANTA ROSA
FALLA FALLA FALLA FALLA DE DE DE DE URICA URICA URICA URICA
C. C. C. C. DE DE DE DE ANACO ANACO ANACO ANACO
BLOQUE DE
BERGANTÍN
Cretácico Inferior
Cretácico Superior
Paleoceno - Mioceno
FALLA FALLA FALLA FALLA EL EL EL EL PILAR PILAR PILAR PILAR
FALLA FALLA FALLA FALLA
FRANCISCO
FRANCISCO
FRANCISCO
FRANCISCO
SAN SAN SAN SAN
Rocas
ígneo-metamórficas
C. C. C. C. DE DE DE DE PIRITAL PIRITAL PIRITAL PIRITAL
BLOQUE DE
CARIPE
C. C. FURRIAL FURRIAL
C. C. OROCUAL
OROCUAL
FRENTE FRENTE FRENTE FRENTE DE DE DE DE DEFORMACIÓN
DEFORMACIÓN
DEFORMACIÓN
DEFORMACIÓN
Retrocorrimientos
Corrimientos Fallas Transcurrentes
Figura 1-5 Mapa Estructural de la Serranía del Interior y el frente de corrimientos
(Modificado de Bosset y Gou, 2004).
Los corrimientos de Pirital, Orocual y Furrial, representan el límite sur del dominio para-
alóctono. Esta cuenca contiene entre 5 y 8 kilómetros de sedimentos neógenos y se extiende
hacia el Sur hasta el Río Orinoco, lo que marca el límite entre los afloramientos de la
Serranía y los afloramientos del basamento del Escudo de Guayana. (Roure et al., 1994)
El Rasgo estructural característico de esta provincia es el corrimiento de Pirital, el cual
corre paralelo a la Serranía del Interior. En esta provincia de transición se observa también
un alineamiento de volcanes y Diapiros de Barro, paralelos a los anticlinales presentes
(Subieta et al., 1987).
9

1.3.2. Evolución Téctono-Estratigráfica de la Sub-Cuenca de Maturín
La evolución geodinámica de la Cuenca, a partir del Paleozoico, se puede dividir en 4
grandes episodios (Eva, Burke, Mann y Wadge, 1989 citado en Parnaud et al., 1995):
1. Fase de Pre Ruptura en el paleozoico
2. Fase de Ruptura y Deriva durante el Jurásico y el Cretácico temprano, caracterizada
por la generación grabenes.
3. Período de Margen Pasivo o Deriva Continental durante el Cretácico y Paleógeno
4. Fase final durante la Colisión oblicua, en el Neógeno y Cuaternario, que resultó en
la formación de la Serranía del Interior y en la transformación de cuenca de margen
pasivo en una cuenca de antepaís. Esta colisión de la placa del Caribe con la placa
de Sur América fue diacrónica con el desplazamiento Oeste-Este
Fase de Pre Ruptura (Pre- Rifting)
El área donde se formó la Cuenca Oriental era parte del súper continente Pangea. Se ha
sugerido que esta área fue un escenario Paleozoico continental a marino marginal antes de
la fase de ruptura, basándose en los ambientes y en las tasas de sedimentación de las rocas
que se depositaron aquí, correspondiendo de esta forma a una megasecuencia de edad
Paleozoico. (González de Juana et al., 1980)
Megasecuencia de Rift
Se desarrolló durante el Jurásico Tardío – Cretácico Temprano y ha sido descrito como la
Formación La Quinta, al oeste del área de estudio en el graben de Espino. Esta formación,
que se depositó en un ambiente continental, está compuesta por lutitas rojas con sills
basálticos. Esta Megasecuencia alcanza un espesor de 0,6 segundos o de aproximadamente
3600 metros.
Megasecuencia de Margen Pasivo
Esta Megasecuencia ocurre del Cretácico al Paleógeno y está caracterizada por tres fases
transgresivas principales que se desarrollaron de norte a sur y culminaron durante el
10

Turonience, el Paleoceno – Eoceno Temprano y el Oligoceno, respectivamente. La base de
esta megasecuencia no se reconoce ni en pozos ni en afloramientos.
La primera fase transgresiva comienza con la depositación de areniscas basales de la
Formación Barranquín. El máximo avance transgresivo está marcado por la depositación de
la plataforma carbonática, desarrollándose ambientes marinos hacia el norte y más
continentales hacia el sur. Esta transgresión se define como la secuencia Cretácica
intermedia.
En el área de El Pilar, hacia el Norte, se depositan calizas y areniscas de la Formación
Barranquín en el Barremiense y en el área más meridional de los afloramientos se depositan
intercalaciones de areniscas, calizas y lutitas correspondientes a la Formación El Cantil en
el Albiense. Más al norte y noroeste se deposita la formación Tigre con sus miembros
Guavinita, Calizas de Infante y La Cruz.
En el área de El Furrial, se halló otra facie intermedia entre las formaciones Tigre y San
Antonio, la cual está caracterizada por areniscas ricas en materia orgánica, limonitas
dolomíticas con contenido fósil (fragmentos de equinodermos) y presencia de glauconita.
Son comunes capas fosfáticas parcialmente silificadas asociadas con limolitas dolomíticas.
La principal roca madre se depositó durante esta fase transgresiva. Son depósitos de
lodolitas ricas en materia orgánica marina depositados en un ambiente batial entre el alto de
Pirital y la falla de El Pilar, y constituyen una facie orgánica mixta que se desarrolló en un
ambiente de plataforma hacia el Sur, entre el alto de Pirital y el frente de deformación.
Esta roca madre corresponde a las formaciones Querecual y San Antonio de edad
Cenomaniense - Campanéense y son contemporáneos con la depositación de la plataforma
de carbonatos hacia el Sur.
Durante el Paleoceno – Eoceno, la siguiente transgresión fue seguida por la regresión del
Mastristiense; esta está representada por la Formación San Juan, la cual en afloramientos
muestra facies regresivas de abanicos submarinos y en el área de El Furrial muestra facies
deltáicas - estuarinas.
11

La transgresión final que se desarrolló durante el Oligoceno comenzó con la depositación
de la arenisca basal de la Formación Merecure. Esta arenisca es el principal reservorio del
área de El Furrial. La Formación Merecure se depositó en un ambiente continental en la
parte sur de la cuenca (área de Cerro Negro) y en un ambiente interior de plataforma en la
parte norte donde los afloramientos de la Serranía del Interior están compuestos
principalmente por alternancias de areniscas (de grano fino a grueso) con lutitas; la fuente
de sedimentos proviene inicialmente el Escudo Guayanés.
Megasecuencia de Colisión Oblicua
La magasecuencia de margen pasivo finaliza durante el Oligoceno cuando se da la
colisión de la placa del Caribe contra la placa de Sur América y la cuenca se transforma en
una cuenta de antepaís. La colisión oblicua migró progresivamente hacia el Este durante el
Oligoceno Tardío al Mioceno Temprano, dividiendo la cuenca de antepaís en tres áreas: (1)
el área sur (desde Cerro Negro hasta Oritupano) que corresponde con la zona de
plataforma, (2) la zona central (desde Acema-Casma hasta Pirital) que corresponde a la
antefosa, y (3) el área norte (al norte de la falla de Pirital) que corresponde al área de
sobrecorrimiento.
En la parte oriental se depositaron los sedimentos fluvio–deltáicos y de plataforma de las
Formaciones Merecure y Oficina integrados por intercalaciones de areniscas, lutitas y
lignitos. La Formación Oficina alcanza espesores máximos de 7500 pies en Guárico
Occidental y presenta una disminución constante hacia el este en el Edo. Anzoátegui. Las
Formaciones Chaguaramas, Merecure y Oficina son las principales productoras de petróleo
en la Cuenca Oriental.
Durante el Plioceno-Pleistoceno en gran parte de la Cuenca Oriental se establecieron
ambientes fluvio–deltáicos y neríticos costeros, consecuencia del continuo retroceso del
mar hacia el este. En estos ambientes se sedimenta la Formación Las Piedras, caracterizada
por arcillas, limolitas, lignitos, arenas y en algunas partes conglomerados.
Durante el Mioceno se producen cambios fundamentales en la forma de la cuenca y en la
distribución de los sedimentos. En la parte norte de Guárico la sedimentación se hace
12

predominantemente continental. El aparato deltaico de Oficina se extiende
considerablemente hacia el sur y sureste hasta una línea cercana al curso actual del río
Orinoco. Hacia el este la profundización de ambientes continúa determinando la
sedimentación de Carapita. Durante este período el tectonismo adquiere especial
importancia y se manifiesta con mayor intensidad en el flanco norte de la Cuenca.
(González de Juana et al, 1980)
La culminación del proceso sedimentario de la Cuenca Oriental está representada por la
Formación Mesa, de ambiente continental. Se extiende sobre los llanos orientales de
Guárico, Anzoátegui y Monagas. Esta formación es considerada de edad pleistocena debido
a su posición discordante sobre la Formación Las Piedras del Plioceno.
Finalmente después de estos procesos de desarrollo y sedimentación de la Cuenca
Oriental de Venezuela se tiene la columna estratigráfica resultante para la zona costera y
continental del oriente venezolano (Figura 1-6)
En la Tabla 1.1 se presenta además una breve descripción tomada del Código Geológico
de Venezuela 2007, de las unidades estratigráficas que intervienen en el modelado de
velocidades del campo Santa Bárbara.
13

Figura 1-6 Columna estratigráfica de la Cuenca Oriental de Venezuela.
(Tomado del Léxico Estratigráfico de Venezuela)
14

Tabla 1-1. Descripción breve de formaciones geológicas de interés.
FORMACIÓN CRONOLOGÍA LITOLOGÍA DOMINANTE
Fm. La Mesa Pleistoceno Depositos horizontales, deltaicos y palustres
Fm. Las Piedras Mioceno Superior a
Plioceno
15
Sedimentos finos mal consolidados, incluyen
areniscas y limolitas carbonáceas, lutitas
arcillosas
Fm. La Pica Mioceno Superior Alternancia de lutitas y arenas
interlaminadas con arenas y limos finos
Fm. Carapita Oligoceno al Mioceno Lutitas marinas
Fm. El Cantil Aptiense al
Cenomaniense
Fm. Querecual Cenomaniense al
Santoniense
Fm. San Antonio Santoniense al
Maestrichtiense
Calizas fosilíferas macizas a cantidades
apreciables de arenisca lutitas y caliza
Carbonatos pelágicos en alternancia con
calizas lutiticas
Calizas y lutitas negras con numerosas capas
de areniscas duras de color gris claro.
Fm. San Juan Maestrichtiense Alternancia monótona de arenicas duras de
Fm. Vidoño Maestrichtiense Superior
al Eoceno Inferior
grano fino, intercaladas con capas
centimetricas de lutitas negras localmente
calcaras y limolitas negras
Lutitas negras ricas en foraminíferos, capas
menores de areniscas y limolitas calcareos
Fm. Caratas Eoceno Medio Dos secuencias litológicas: pelitico arenosa
que es la formación Caratas y carbonática,
bioclástica llamada Miembro Tinajitas

1.3.3. Campo Santa-Bárbara/Pirital
El campo consiste en un corrimiento buzando hacia el norte con un pliegue asociado de
rumbo N75°E y se subdivide en dos bloques. El Bloque Este y el Bloque Oeste se
encuentran separados por una zona de fallas transcurrentes dextrales. El Bloque Oeste está
caracterizado por rampas laterales y retrocorrimientos y el Bloque Este está definido
igualmente por el corrimiento y por fallas transcurrentes de rumbo N30°W. Ya hacia el
oeste se presentan las rampas laterales de Urica y hacia el este el anticlinal de Carito lo cual
suma complejidad a la zona en ambos extremos (Figura 1-7 y Figura 1-8).
Las arenas más prospectivas en el Campo Santa Bárbara se ubican en los sedimentos de
edad Crétacico – Paleoceno (Formación San Antonio y San Juan) y del Terciario
(Formación Naricual y Carapita Basal), con un espesor promedio de 1000 a 4000 pies. La
profundidad promedio de los pozos es de 16800 pies (Figura 1-7).
Figura 1-7 Corte geológico a partir de la interpretación sísmica, Bloque Este
(Tomado de Bosset y Gou, 2004)
16

Figura 1-8 Cortes geológicos a partir de la interpretación sísmica
(Tomado de Bosset y Gou, 2004)
Los datos sísmicos 3D del campo Santa Bárbara son el producto de la unión de varios
levantamientos independientes adquiridos en el Norte de Monagas (Figura 1-9). El cubo
compilado (merge) resultante, utilizado en la interpretación y modelado estructural vigente
del área, presenta un tamaño de malla (bin) de 20 x 20 metros, y una frecuencia
predominante de 20 Hertz (Figura 1-10). La calidad de los datos en tiempo es en general de
buena a regular, y disminuye progresivamente en sentido Este – Oeste, de tal manera que el
área correspondiente al campo Carito presenta mejor calidad que la correspondiente a los
Campos Santa Bárbara, Pirital y Bosque (Parra, 2006).
17

Figura 1-9 Desarrollo de la sísmica en el campo Santa Bárbara entre los años 1986 a 1999.
1100000
1080000
1060000
1040000
380000 380000
21
1
ESTADO ESTADO ANZOATEGUI
ANZOATEGUI
ESTADO ESTADO MONAGAS
MONAGAS
400000 400000
FF
2
EE
FIGURA FIGURA 8.16 8.16
DD
5
25 28
51
22
37 36
49
53
18
(Tomado de U.E Pirital, PDVSA)
420000 420000
CC
43
FIGURA FIGURA 8.14 8.14
6
26
24
42
31
AA
BB
34
19
46
50
52 52
35
47
54
55
8
15
FIGURA FIGURA 8.15 8.15
FIGURA FIGURA 8.8 8.8
14
9
4
SBC-4 30 45 48
440000
0m 10000m 20000m
20
MATURIN
Figura 1-10 Mapa de ubicación de la información sísmica (Tomado de Parra M. 2006)
3
FIGURA FIGURA 8.7 8.7
460000
NORTE
480000
PAG-3X
Lìmite entre Estados
Lìneas Sìsmicas 2D
Cubo Sìsmico
(tiempo y profundidad)
Pozos Claves
Campos Petroleros
18

1.3.4. Campo Tácata
Al Oeste del Campo Santa Bárbara se encuentra el campo Tácata. El campo está
prácticamente dividido por dos Sistemas de Fallas Principales o Mayores y un
retrocorrimiento: El Sistema de Fallas de Urica – Pirital al Este; El Sistema de Fallas de
Tala al Norte; y el retrocorrimiento al Sur.
La estructura Tácata-Tacat es tectónicamente compleja, con rumbo E-O y está ubicada
dentro de una zona triangular limitada hacia el norte por los corrimientos mayores de Tala y
Pirital, y hacia el sur por el frente de deformación definiendo el techo de la zona triangular,
lo que se puede distinguir en la Figura 1.5. Asociados a estas estructuras principales se
encuentran una serie de corrimientos y retrocorrimientos menores, que dan origen a las
diferentes repeticiones de las formaciones Capaya y Carapita, donde se encuentran los
yacimientos de interés.
La estructura de Tácata se define como un anticlinal fallado dentro de una zona triangular
que involucra unidades que comprenden edades desde el Mioceno temprano hasta el
Mioceno Medio. La zona triangular se encuentra definida por tres corrimientos (T0,T1 y
T2) y tres retrocorrimientos (BT1, BT2, BT3), donde los retrocorrimientos se forman en la
etapa tardía del Mioceno Temprano y los corrimientos en la etapa temprana del Mioceno
Medio, involucrando las secuencias depositadas durante el Mioceno Temprano hasta el
Mioceno Medio. (Figura 1-11)
19

Figura 1-11 Rasgos estructurales principales definidos en la estructura de Tácata (Tomado
de Informe técnico área Tacata-Tacat, PDVSA 2006)
El cubo sísmico empleado para el estudio del Campo Tácata-Tacat abarca un área
aproximada de 750 km 2 , donde se ha realizado la unión de levantamientos sísmicos 3D
adyacentes, generando lo que se conoce como Volumen Integrado Tacat_2000, el cual
comprende la integración de los volúmenes sísmicos: Bosque 98, Casupal 91, Casupal 92,
Pato 98 y ASTRA 98 (Figura 1-12).
20

10
70
00
0.
0
10
60
00
0.
0
10
50
00
0.
0
380000.0 390000.0 400000.0 410000.0 420000.0
CONVENIO
QUIAMARE - LA CEIBA
ASTRA 3D/98
AN
ZO
AT
EG
MO
NA
GA
S
PATO 3D/98
CASUPAL NORTE 3D/92
UI BOSQUE 3D/98
CASUPAL SUR 3D/91
0 10
380000.0 390000.0 400000.0 410000.0 420000.0
Figura 1-12 Levantamientos sísmicos 3D que conforman el Volumen Integrado Tacat_2000
(Tomado de Informe técnico área Tacata-Tacat, PDVSA 2006).
Km.
21
10
70
00
0.
0
10
60
00
0.
0
10
50
00
0.
0

2.1. Métodos Sísmicos
CAPÍTULO 2
MARCO TEÓRICO
Los métodos sísmicos representan una herramienta muy importante dentro de la
prospección geofísica. La sísmica de reflexión es usada como método geofísico de
exploración para la búsqueda que hidrocarburos debido al gran poder que tiene de
resolución y penetración no invasiva.
Su uso se enfoca principalmente en la caracterización del suelo mediante la grabación de
los tiempos de viaje o llegada, variaciones de amplitud, fase y frecuencia de las ondas que
se introducen en el subsuelo, todo esto con el objetivo de reconstruir la disposición de los
estratos y rocas en el subsuelo así como también obtener sus características físicas.
La medición del tiempo de viaje de las ondas entre las capas con diferentes propiedades
físicas es la base de los métodos sísmicos; este fenómeno esencialmente sucede por
contraste entre impedancias acústicas. La energía sísmica se genera mediante un impacto o
explosión controlada en superficie. En un medio homogéneo las ondas se propagan
esféricamente desde la fuente a través del subsuelo interaccionando con los distintos
estratos, de manera que una parte de la energía que transportan se refleja y permanece en el
mismo medio que la energía incidente, y el resto se transmite al otro medio con un fuerte
cambio de la dirección de propagación. (Telford, 1990)
2.2. Velocidades Sísmicas
Las ondas se propagan en el subsuelo constantemente, pero esta propagación depende
totalmente de las propiedades elásticas de las rocas, por eso es importante discutir algunos
conceptos de elasticidad y medios elásticos.

La elasticidad es la propiedad mecánica de los materiales de sufrir deformaciones
reversibles cuando se encuentran sujetos a la acción de fuerzas exteriores y de recuperar la
forma original si estas fuerzas exteriores se eliminan, también se traduce como la capacidad
de los materiales a resistir estos cambios (Telford, 1990).
La velocidad de onda es la velocidad a la cual se transmite la onda de energía en un
medio. Ésta depende de las propiedades características del medio, en el caso de ondas
sísmicas depende de los parámetros y la densidad. La velocidad de onda se relaciona
siempre con la frecuencia (f) y la longitud de onda (λ) por la ecuación:
2.2.1. Velocidad interválica Vint
23
v = f λ (2.1)
Es la velocidad sísmica sobre un intervalo de profundidad “Z”. Si el tipo de roca es
uniforme en ese intervalo, entonces la velocidad de intervalo es igual a la velocidad de la
formación. Si el intervalo de profundidad abarca un número de estratos de roca de
diferentes litologías, entonces la velocidad de intervalo es igual a la velocidad promedio
(V) calculada sobre la distancia “Z”. Si Z i es el espesor del i-ésimo intervalo y t i es el
tiempo de recorrido de un camino a través de él, se tiene:
2.2.2. Velocidad promedio (V)
Vint= zi / ti (2.2)
Es el promedio de todas las velocidades interválicas desde la superficie hasta un horizonte
en particular (Cantos, 1974). Estas velocidades son medidas durante la corrida de un
registro de velocidades en el pozo (“check shot”), en donde un geófono es introducido en
este a diferentes profundidades, y una fuente sísmica es activada en la superficie.
La velocidad es obtenida calculando los tiempos de viaje Δt, de las ondas sísmicas, desde
la posición de la fuente hasta el geófono. Con los datos obtenidos, la velocidad promedio
desde la posición de la fuente hasta la profundidad donde es colocado el geófono, es
calculada mediante:

V = ∑ Vi . Δti (2.3)
∑ Δti
donde Vi es la velocidad interválica y Δt los tiempos de viaje. Esta medida se utiliza para la
construcción de curvas tiempo – profundidad (T-Z), las cuales permiten realizar la
conversión a profundidad de datos sísmicos en tiempo.
2.2.3. Velocidad promedio (Vrsm)
Es definida como la raíz cuadrada de la velocidad promedio ponderada en tiempo, cuya
expresión matemática es:
2.3. Migración sísmica
Vrms² = ∑ Vi² . Δti (2.4)
∑ Δti
Los programas óptimos de exploración y explotación requieren la mejor “imagen” posible
del subsuelo. La migración sísmica es una de las fases básicas del proceso de datos que
optimiza y presenta la mejor “imagen” del subsuelo. La migración sísmica corrige el
desplazamiento geométrico de los datos de un reflector buzante o los cambios laterales de
velocidad, y coloca los reflectores en su verdadera posición espacial en lugar del CMP
entre la fuente y el receptor. Esto se logra a través de diversas soluciones de la ecuación de
las ondas elásticas a través de las rocas.
La migración puede realizarse en dos dominios, el tiempo o la profundidad, y antes o
después del apilamiento. Ciertos problemas de generación de imágenes se pueden resolver
con la migración en domino del tiempo, pero los más complejos requieren la migración en
el dominio de la profundidad. La migración en profundidad es un procesamiento muy
sensible al modelo de velocidad. Si no existe un buen control de las velocidades con la
información de pozo, se pueden presentar inconsistencias entre las reflexiones asociadas a
niveles estratigráficos y las correlaciones geológicas en los pozos.
24

2.4. Modelo de velocidades
La migración en profundidad antes del apilamiento, es capaz de crear una imagen de los
reflectores de superficie que se encuentren bajo la presencia de gradientes laterales severos
de velocidad (Kim et al, 1996); además de que se puede aplicar a las estructuras y los
campos de velocidad más complejos.
La construcción de los modelos de velocidades, típicamente es iterativo y requiere
muchas corridas computacionales intensivas de migración antes del apilamiento. La
migración iterativa en profundidad antes del apilamiento es utilizada para generar un
modelo de velocidades exacto.
La combinación de la migración de rayo con la migración en profundidad antes del
apilamiento, permite la construcción iterativa de modelos de velocidades sin la necesidad
de generar migración en profundidad de los CMP “gathers” en cada posición de superficie.
El análisis convencional de velocidades, es una de las más laboriosas tareas en el
procesamiento sísmico, pero puede ser evadido actualizando las funciones de velocidad de
capas basadas sobre la inspección visual de los “gathers” migrados en profundidad.
El objetivo del procesamiento sísmico es conseguir que los datos de reflexiones del
subsuelo se asemejen lo más posible a una sección geológica de la tierra a lo largo del
plano del arreglo de los geófonos. La presencia de registros de pozos en el área facilita y
optimiza la creación del modelo de velocidades, como es nuestro caso, la presencia de
registro sónico y de “check shot”, nos da información de las tendencias verticales de las
velocidades en el área, y por lo tanto, con esta información podemos diferenciar paquetes
que sigan los mismos patrones.
2.5. La herramienta sónica
La herramienta sónica básicamente mide la velocidad del sonido en las formaciones
penetradas por un pozo. La herramienta presenta un transmisor y dos receptores colocados
en el cuerpo de la sonda, con un espaciamiento típico entre el transmisor y el receptor más
cercano de 3 pies, y de 2 pies entre los receptores, como podemos ver en la figura 2-1.
25

Figura 2-1 Herramienta sónica. Tomado de Dewan (1984).
El principio físico de la herramienta consiste en aplicar un pulso de corriente o voltaje al
transmisor, por lo que éste generará un pulso de presión corto y oscilatorio a una frecuencia
de alrededor de 25 kHz, en el lodo. Este proceso iniciará seis diferentes ondas viajando
hacia arriba y hacia abajo por el hoyo (figura 2-2). Estas ondas son la compresional y la de
cizalla las cuales se refractan a través de la formación, dos ondas directas que viajan a
través de la herramienta y del lodo, y finalmente dos ondas superficiales a lo largo de la
pared del hoyo que son la pseudo-Rayleigh y Stoneley. Las velocidades de estas ondas
varían entre los 4000 pies/s y los 25000 pies/s.
Al poco tiempo que el transmisor ha pulsado, el receptor cercano detecta la llegada de
varios frentes de onda, los cuales, un poco después, llegan hasta el detector más lejano. La
secuencia normal de las llegadas de las ondas se caracteriza por la llegada de la onda
compresional de primero, seguida inmediatamente por la onda de cizalla que viaja como
una onda compresional en el lodo.
El interés principal está en las ondas compresionales y de cizalla. Por definición, una
onda compresional es aquella en la cual las partículas del medio vibran en la misma
dirección en el cual se propaga la energía, en este caso una dirección paralela al pozo. En
onda de cizalla por otra parte, las partículas vibran en una dirección perpendicular a la
dirección de propagación, en este caso perpendicular al eje del pozo.
26

Figura 2-2 Diagrama de viaje de las ondas. Tomado de Dewan (1984)
La herramienta sónica estándar mide solo el tiempo de viaje de las ondas compresionales.
Para lograr esto, el transmisor es pulsado una vez y un circuito electrónico mide el tiempo
hasta la primera incursión negativa de la llegada compresional en el receptor más cercano.
El transmisor se vuelve a pulsar y el circuito mide el tiempo hasta el receptor lejano. La
diferencia de tiempos es calculada y dividida por la distancia entre los receptores. El
resultado es presentado en un registro como tiempo de tránsito en la formación, en
microsegundos/pie (μs/pie). La precisión de la medida es bastante buena, aproximadamente
de más o menos 0.25 μs/pie.
2.6. La herramienta sónica compensada
La primera generación de herramientas sónicas, las de un solo transmisor y dos
receptores, sufría el defecto de que cuando existía un hoyo agrandado o cuando la
herramienta presentaba una inclinación (figura 2-3) con respecto al eje del pozo, el tiempo
de tránsito medido era mayor.
27

Figura 2-3 Problemas con las herramientas de un solo transmisor.
Tomado de Western Atlas (1992).
Se realiza una secuencia de 4 mediciones de tiempo (figura 2-4): tiempo desde el
transmisor 1 hasta el receptor 1 o T11, desde el transmisor 1 hasta el receptor 2 o T12,
desde el transmisor 2 hasta el receptor 1 o T21 y desde el transmisor 2 hasta el receptor 2 o
T22. El tiempo de tránsito (lentitud) correcto es obtenido al promediar las cuatro lecturas de
tiempo en μs, y dividir el resultado entre la separación de los receptores, en pies,
en donde Δt estará dando μs/pie.
Δt = ( ( T12 - T11 ) – ( T22 – T21 ) / 2 ) / X (2.5)
28

Figura 2-4 Secuencia de las mediciones sónico compensado. Tomado de Dewan (1984)
Esta técnica elimina los efectos de hoyo agrandado y de inclinación de la herramienta y
con ella, la calidad de los registros sónicos en hoyos muy derrumbados fue mejorada
grandemente. Normalmente, la herramienta sónica es corrida centrada en el hoyo para que
las contribuciones a la señal en el receptor de diferentes partes del hoyo estén en fase (si el
hoyo es redondo) y sea maximizada la relación señal/ruido.
2.7. Disparos de verificación sísmica
Es un tipo de medición sísmica realizada en los pozos, diseñada para medir el tiempo de
viaje sísmico a partir de la superficie hasta una profundidad conocida. La velocidad de las
ondas P en las formaciones encontradas en la pared de las perforaciones puede ser medida
directamente mediante un geófono que es bajado a cada formación de interés, desde la
superficie de la tierra se genera y se envía una fuente de energía, y este geófono recibe las
perturbaciones u ondas creadas por la fuente, las cuales viajan por las formaciones hasta el
receptor donde se graban las mediciones de los tiempos de viaje. Las mediciones pueden
ser correlacionadas con la data sísmica de superficie haciendo uso de las correcciones del
registro sónico y generando un sismograma sintético para confirmar o modificar la
interpretación sísmica (Schlumberger Oilfield Services, 2006).
29

Los disparos de verificación sísmica o “check shots” se llevan a cabo empleando una
fuente de ondas compresionales (tal como un cañón de aire o un camión vibrador), el cual
se encuentra ubicado en la superficie. Además de la fuente de energía, se utiliza en
superficie un sistema de control, un sistema registrador, un compresor de aire (sólo cuando
se utiliza el cañón de aire), un hidrófono y geófonos de superficie.
De la misma manera el instrumental de pozos consiste en un dispositivo de anclaje, un
sistema hidráulico de alta presión, un geófono y un amplificador. La herramienta de pozo
de uso más común es la herramienta de sísmica de pozos (Well Seismic Tool o WST) con
un arreglo de cuatro geófonos uniaxiales que han sido diseñados para operaciones en
ambientes de alta temperatura.
El producto final son todos los datos grabados digitalmente y en un procesador
(computador) con el resultado del apilamiento de los trenes de ondas por cada profundidad
de disparo. En estos registros digitales se muestra el tiempo medido entre la señal de
referencia del hidrófono y la señal del geófono dentro del pozo. Los registros sónicos que
han sido calibrados a través de “check shots” son la referencia básica para llevar una
sección sísmica en tiempo al dominio de la profundidad, proceso conocido como
conversión tiempo – profundidad.
Los tiempos de tránsito grabados de las llegadas directas son usados para ajustar el
registro sónico por deriva, que es la diferencia entre los tiempos sísmicos y los tiempos
integrados del sónico; la cual es causada principalmente por dispersión y efectos de pozo.
Los disparos de verificación sísmica difieren del perfil sísmico vertical (VSP), en el
número y densidad de las profundidades recibidas y registradas; las posiciones de los
geófonos pueden ser ampliamente e irregularmente localizadas en el fondo del pozo; sin
embargo el perfil sísmico vertical, usualmente, tiene numerosos geófonos posicionados en
intervalos cercanos y regularmente espaciados en el fondo del pozo
30

El registro sónico y los disparos de verificación sísmica se diferencian en la geometría de
la adquisición, en el tipo y ubicación de la fuente utilizada. Los disparos de verificación se
utilizan principalmente para generar curvas de tiempo-profundidad, con la cual se calibra el
registro sónico y se generan los sismogramas sintéticos. La velocidad interválica se puede
derivar de los disparos de verificación mediante la fórmula de Dix que relaciona la
velocidad RMS y los tiempos de viaje de los reflectores (Schlumberger, Glosario en línea)
La fórmula de Dix
Vint = {(t2 * Vrms2 2 – t1 * Vrms1 2 )/ (t 2 – t1)} 1/2
(2.6)
Los registros sónicos se utilizan con la finalidad de medir con mayor detalle las
velocidades de los paquetes o formaciones geológicas de interés. Es una potente
herramienta para obtener información sobre fracturación y cambios litológicos.
31

CAPÍTULO 3
METODOLOGIA
3.1. Recopilación de datos de registros de pozos
Para la construcción del modelo 3D de velocidades interválicas a partir de registros
sónicos y Check Shot, para el Campo Santa Bárbara/Pirital al norte del estado Monagas se
siguió la metodología que se describe a continuación.
En principio se realizó una revisión de la información geológica disponible de la Cuenca
Oriental de Venezuela para obtener de esta forma una ubicación geográfica, espacial y
teórica de la complejidad estructural y geológica presente en la región. También se estudió
la geología local, la evolución tectónica, la estratigrafía y la sedimentología de la Sub
Cuenca de Maturín.
Cabe destacar que el proyecto de obtención del modelo de velocidades 3D para el Campo
Santa Bárbara fue iniciado con anterioridad por el Ing. Manuel Bolívar por lo que la
obtención y selección de datos de registros de pozos del área ya había sido realizada y se
compilaron una serie de datos de pozo. La selección de pozos básicamente estuvo
condicionada por la existencia de registros sísmicos (Check Shot) y registros sónicos.
También se realizó un compendio de trabajos previos de Campo Santa Bárbara, entre los
cuales se haya informes y tesis con metodologías distintas para el modelo de velocidades
interválicas. Además se compilo información del modelo estructural vigente del campo,
tesis e informes de interpretaciones estructurales y estratigáficas.

3.2. Control de Calidad de Registros de Pozos
Dentro de la exploración y prospección petrolera es de gran importancia la aplicación de
métodos petrofísicos. La importancia de éstos radica en que se obtienen mediciones que
abarcan desde la parte somera hasta la mayor profundidad posible dentro de un pozo, es
decir, profundidad de yacimiento, la cual varía dependiendo del área del estudio. En el caso
del Campo Santa Bárbara, los pozos seleccionados poseen registros desde profundidades
someras (entre 0 y 1000 pies), hasta profundidades de yacimiento (14000 pies en adelante).
Otra de las cosas que se tomó en cuenta para la discriminación de los registros de pozos es
que al menos en los datos se observase de manera marcada las tendencias de cambios de
velocidad en intervalos que representen o correspondan a cambios geológicos.
Al tener los datos de registros de pozos, se comenzó por hacer una revisión detallada de
cada uno. En principio se verificó la ubicación geográfica de los pozos con registros
sónicos y Check Shot, lo que en total suma 43 pozos con sónico y 56 Check Shot. En la
figura 3-1 se observa la ubicación y distribución de estos pozos.
Norte (UTM)
1076000
1073000
1070000
1067000
1064000
1061000
TAG-13
TAG-12
TAG-17
TAG-15
TAG-18
TAG-16
TAG-21
TAG-11E
TAG-20
TAG-19
Pozos con Registro Sónico
Campo Santa Bárbara
TAC-02X
TAC-03
TAG-04
PIC-01E
PIC-17
SBC-33
SBC-50
Figura 3-1 Ubicación geográfica pozos con registros sónicos en el área
PIC-25
TAC-01E
SBC-37E
PIC-26
PIC-06E
SBC-15
SBC-49
SBC-39
SBC-130E
PIC-12
MGC-1E
PIC-08
SBC-35
SBC-11
SBC-136
PIC-28
SBC-32
SBC-04
SBC-31
SBC-142
SBC-143
SBC-22E
SBC-131
MUC-23
SBC-12
1058000
385000 395000 405000 415000 425000 435000 445000
Este (UTM)
SBC-10E

El modelo de velocidades interválicas puede obtenerse haciendo solamente uso de los
disparos de verificación, pero para tener un modelo confiable, sólido y exacto es necesaria
la integración de éstos con los registros sónicos.
El disparo de verificación es un registro sísmico que reporta el tiempo de las primeras
llegadas de ondas compresionales a profundidades conocidas y con una serie de
correcciones y cálculos se obtiene la velocidad interválica de dos niveles del Check Shot
(Dewan, 1984). La velocidad de intervalo procede de la fórmula de DIX (ecuación 2.6) o de
la profundidad vertical medida desde el SRD (Seismic Referente Datum) o nivel del mar y
el tiempo de viaje vertical desde el SRD.
Tabla 3-1. Check Shot SBC12
MD TVDSS TWT Vint Vint
(ft) (ft) (ms) (ft/s) (ft/s)
642 0 0
1600 -958 367 5224 5220,71
2600 -1958 640 7324 7326,01
3400 -2758 847 7708 7729,47
4062 -3420 1018 7750 7742,69
4915 -4273 1225 8255 8241,55
5900 -5256 1461 8330 8330,51
7170 -6522 1747 8862 8853,15
8000 -7349 1916 9749 9786,98
9100 -8446 2162 8933 8918,70
11000 -10340 2573 9219 9216,55
12610 -11943 2878 10522 10511,48
13700 -13029 3077 10902 10914,57
14500 -13827 3240 9754 9791,41
15550 -14876 3438 10603 10595,96
16664 -15989 3658 10116 10118,18
17000 -16325 3704 14600 14608,70
18000 -17358 3827 16888 16796,75
En este caso la velocidad interválica se calcula de acuerdo a la expresión siguiente:
Vint = (Profvert 2– Profvert 1)/ (Tviaje 2 – Tviaje 1)
34

donde: Vint es la profundidad del intervalo
Profvert 1 es la profundidad vertical menor
Profvert 2 es la profundidad vertical mayor
Tviaje 1 es el tiempo de viaje menor hasta Profvert 1
Tviaje 2 es el tiempo de viaje mayor hasta Profvert 2
En el estudio de velocidades con los disparos de verificación (Check Shot), el problema
de la dispersión es mínimo ya que la densidad de datos es menor en comparación con el
registro sónico (Díaz, 2006). El inconveniente se presenta cuando en los reportes de
adquisición la información no está completa, o cuando existe algún punto con valores que
escapan de forma drástica y observable de la tendencia general de los demás datos. Al tener
anomalías en un registro, tal como se muestra en la figura 3-2 para el pozo PIC-21,en el
cual se tienen valores muy altos de velocidad para un punto poco profundo, simplemente no
se toma en cuenta este valor anómalo y se hace referencia a la omisión para no pasarlo por
alto al hacer el cálculo de tendencias lineales.
Vint (pie/s)
15000
13000
11000
9000
7000
Tendencia Somera en Registro Check Shot
Pozo PIC-21
5000
0 1000 2000 3000
TVDSS (pie)
4000 5000 6000
Figura 3-2 Registro Check Shot pozo PIC-21 con valor anómalo
35

La ventaja de los Registros Sónicos es su alta resolución. Poseen una gran densidad de
mediciones; estas se hacen por lo general cada ½ pie en profundidad. En vista de la
cantidad de datos presentes en cada uno de los registros en este estudio, se tomó la decisión
de hacer un muestreo que abarcara solo valores correspondientes a mediciones cada pie;
por ejemplo registros que tenían una densidad de datos de al menos 30 mil valores cada ½
pie pasaron a contener solo 15 mil valores, para tener menor volumen de datos al momento
del análisis y manipulación de los mismos.
A pesar de la reducción de datos, existen variaciones significativas para los valores de
profundidad mínima, lo que resulta de la densidad de mediciones que posee el registro
sónico. Esto se puede observar en la figura 3-3 generada con Microsoft Excel 2003
La herramienta Sónica proporciona el tiempo de tránsito (µ/pie) de una onda dentro de la
formación; por esto a partir de los registros sónicos se hizo en primer lugar el cálculo de las
velocidades con la simple aplicación de una conversión de medidas para obtener los valores
de velocidad (pie/s). La conversión es la siguiente:
36
Vint = 10 6 * (∆t) -1 (Ec. 3.1)
donde Vint es la velocidad interválica y ∆t es el tiempo de tránsito.

Figura 3-3 Registro Sónico de pozo PIC-8
37

Debido a la dispersión de los datos presentes en el registro sónico se aplicó un método de
“suavizado” de la curva. En la figura 3-4 se pueden ver las medianas que fueron calculadas
para el conjunto de valores de lentitud del pozo PIC-28. Fueron deducidas medianas de
100, 200 y 300. El mejor ajuste lo da un cálculo de mediana de 300 datos por que
representa las variaciones de las curvas sin hacer un suavizado muy fuerte y/o reproducir la
dispersión de datos mientras que para 100 y 200 se mantienen picos abruptos, (Díaz, 2006).
Además se observa que con una mediana de 300 se mantiene relativamente la variación real
de la curva, sin perder tendencias o detalle importante y logrando eliminar los picos o
puntos anómalos dentro del registro. Este proceso se aplicó a todos los registros sónicos
con el fin de mejorar el muestreo y eliminar los problemas de dispersión y las variaciones
exageradas a escasos pies de profundidad.
La mediana opera de modo simple, ya que un conjunto finito de valores los ordena de
manera creciente y toma aquel valor que divide al conjunto en dos partes iguales, de forma
que el número de datos mayor o igual a la mediana es igual al número de valores menores o
igual a éstos. Si el número de valores muestreados es impar, la mediana es el valor medio,
el cual corresponde al dato Xn/2. Cuando el número de valores en el conjunto es par, no
existe un solo valor medio, sino que existen dos valores medios; en tal caso, la mediana es
el promedio de los valores, es decir, la mediana es numéricamente igual a Xn/2 + Xn/2+1 / 2
(Spiegel et. al, 2003).
38

Figura 3-4 Registro Sónico de pozo PIC-08 con medianas de 100, 200 y 300 datos
39

Figura 3-5 Registro Sónico de pozo PIC-28 con mediana de 300
40

En la figura 3-5 se puede observar como en el registro se disminuyó la dispersión y los
puntos aislados que provienen de malas mediciones y ruido. Despues del control de calidad
de los disparos de verificación y los registros Sónicos, lo que se hizo fue acoplarlos. Este
procedimiento permite observar la manera en que los disparos de verificación honran los
valores de registro sónico y viceversa; a su vez permite verificar errores de medición en el
registro. Es decir, permite hacer un control de calidad previo para identificar posibles fallas
de la herramienta.
Nótese que los puntos de Check Shot son bastante fieles a la curva de Vint derivada del
registro sónico del pozo PIC-28 en la figura 3-6, se puede decir que definen de igual forma
para la mayoría de los casos tendencias generales de velocidad y sus cambios. De este
análisis, para algunos registros de pozos, se pudo observar que la herramienta sónica fallaba
y no obtenía lecturas, mientras que el registro sísmico si contenía información de ese
intervalo, lo cual es de ayuda importante y los valores son muy probablemente confiables
(figura 3-7).
41

Figura 3-6 Registro Sónico y “Check Shot” de pozo PIC-28
42

Figura 3-7 Registro Sónico y “Check Shot” pozo SBC-15
43

3.3. Identificación de Tendencias en registros de velocidad interválica
Luego de disminuir los problemas de dispersión en los datos y de adecuar para el análisis
un gran cúmulo de registro, se comenzó la observación e identificación de los cambios de
velocidad en las curvas. Estas inflexiones son las que permitieron reconocer los distintos
dominios en los datos, los cuales describen una curva de ajuste que relaciona la velocidad
con la profundidad. Después a cada paquete de lecturas agrupadas por tendencias se le
aplicó una regresión por mínimos cuadrados para estimar la ecuación de la tendencia.
Es importante señalar que en este tipo de operación es posible generar un factor
denominado coeficiente de Pearson el cual es un índice estadístico que mide la relación
lineal entre dos o más variables cuantitativas, es decir, indica la relación existente entre la
desviación típica de una o varias muestras y su media (Spiegel et. al, 2003). El valor del
índice de correlación varía en el intervalo [-1, +1] a medida que R 2 es mayor o igual a 1 se
entiende que existe una correlación que va de positiva a perfecta, que es lo que se espera en
buena medida para nuestro estudio.
Estadísticamente un coeficiente de correlación por encima de 0.7 para un ajuste lineal, es
una regresión aceptable; por este motivo podemos asignar la ecuación de velocidad en
función de la profundidad al sector definido por la tendencia que estamos estudiando,
cuando tengamos un coeficiente R 2
mayor o igual a 0.7. (Díaz, 2006)
3.4. Modelo de velocidades
El modelo de velocidades interválicas se obtuvo en “GOCAD”, un paquete de la
compañía “PARADIGM”, cuyo objetivo es desarrollar nuevos métodos de modelado
interactivo de la geometría y de propiedades de objetos geológicos complejos, para
aplicaciones en geología, geofísica e ingeniería de depósitos (GOCAD User Guide 2007).
A continuación se explica de manera breve la metodología seguida para la obtención del
modelo de velocidades.
44

En primer lugar se cargaron los datos de los pozos tales como desviación, velocidad
interválica, calculada a partir del registro sónico y velocidad interválica derivada de los
registros check shot.
El buzamiento de los estratos, la alternancia de litologías y de dureza de los estratos,
ocasiona cierta desviación del pozo, es decir, abandona una trayectoria vertical. Para pozos
muy profundos, en ocasiones el cambio es imperceptible pero, en su mayoría, los pozos
presentan cierta variación, la cual se mide como la diferencia entre las coordenadas Este
(X), Oeste (Y) de una serie de puntos en profundidad con las coordenadas del pozo en
superficie (figura 3-8).
Figura 3-8 Imagen de pozo desviado en GOCAD
45

Las velocidades interválicas calculadas a partir del registro sónico y check shot fueron
cargadas en GOCAD como dos propiedades distintas “Vint_sonico” y “Vint_CS”. Estas
propiedades luego fueron combinadas con la aplicación de un algoritmo en una tercera
propiedad, “Vint_t”. El script da prioridad a los valores de velocidad provenientes del
registro sónico en cada pozo; si no existe registro en el pozo o en ciertas profundidades del
pozo, entonces “Vint_t” tendrá como valores de velocidad los derivados del registro “check
shot”. Se toma como primordial el registro sónico por la densidad de muestreo del registro
(figura 3-9).
Figura 3-9 Imagen con propiedades de velocidad
46

3.5. Construcción de Voxet o malla 3D
El VOXET es un volumen tridimensional que limita el espacio en una región, la cual se
define de tal manera, que dentro estén ubicados los pozos con los cuales se realiza el
estudio. Además, sus dimensiones vienen dadas en coordenadas UTM. La dimensión de las
celdas que se definen en cada eje del VOXET es esencial para la obtención del modelo de
velocidades (figura 3-10). Debe tenerse en cuenta la densidad de pozos presentes en el área
de estudio para así poder diseñar un tamaño ideal, cuyo principio está en ubicar al menos
un pozo por cada una de las celdas definidas. Dentro del VOXET se definen igualmente
propiedades derivadas de los registros de pozos previamente cargados, regiones
provenientes de horizontes interpretados de la sísmica, marcadores de estudios
estratigráficos o datos de pozo.
Figura 3-10 Imagen del Voxet Campo Santa Bárbara/Pirital y Tácata con “celdas”
Dentro de las opciones que ofrece GOCAD para definir el volumen de trabajo, existe el
llamado “SGRID”. Su construcción es similar a la del VOXET, con la diferencia de que el
“SGRID” permite modificar cada celda añadiendo a las superficies de éstas inclinaciones
derivadas de la estratigrafía, lo que constituye una herramienta importante al momento de
respetar capas definidas en la sedimentología de la zona. Su uso se restringe en este estudio
47

ya que la velocidad es una propiedad que no depende de la sedimentología sino que es
inversamente proporcional al cambio en profundidad (Díaz. 2006).
3.6. Carga y modelado de superficies
El control y estudio preliminar de los datos de registros sísmicos “check shot” y registros
sónicos en las graficas realizadas de velocidad interválica, permitió establecer cuatro
horizontes de importancia dentro del área de estudio. La presencia de los horizontes se hizo
evidente en los registros como cambios de tendencia en las curvas de velocidad calculadas.
Estos cambios pueden ser asociados directamente con los cambios litológicos presentes en
las formaciones del área. El tope y base de cada cambio fue cargado en GOCAD dentro de
la propiedad marcadores y estos definen para cada pozo cuatro regiones de importancia
derivadas de los marcadores: Somero, Carapita, Bloque Alóctono y Yacimiento (ver figura
del Capítulo 4 número 1).
Las superficies pueden ser generadas a partir de los marcadores litológicos. Ésta se
ajustaría de manera exacta a la tendencia local de cada tope o base, según sea el caso. Pero,
además, es posible agregar a la base del proyecto en GOCAD, las superficies interpretadas
provenientes del modelo estructural del área. Estas superficies son cargadas en GOCAD
como una malla de puntos de densidad variable. Al ser muy densa la malla con respecto a
su extensión espacial, se aplica un filtro para disminuir la densidad pero sin perder detalle
en la definición de la superficie. La superficie se genera dentro del área de interés mediante
una interpolación DSI. Y además puede ser ajustada a los marcadores presentes en cada
pozo. De esta forma la superficie honra los marcadores definidos y los horizontes
interpretados de la sísmica.
Al estar definidas dentro del Voxet las superficies, deben ser creadas las regiones. Cada
región se construye a partir de un tope y una base, es decir a partir de dos superficies. Es
dentro de estas regiones donde se realizará más adelante las interpolaciones a fin de obtener
el modelo de velocidades. El rango de valores que interviene en el proceso de propagación
se construye igualmente a partir de las regiones construidas a partir de las superficies que
honran los marcadores y los horizontes interpretados de la sísmica.
48

3.7. Propagación de velocidades
La propagación de los datos de velocidad constituye el último paso en la construcción del
modelo de velocidades. Por medio del uso de distintos métodos de interpolación, como DSI
(Discrete Smooth Interpolation) y Kriging se logra poblar el volumen de interés, asignando
así un valor de velocidad a cada celda del VOXET. Los valores de velocidad interválica
provenientes de los registros son usados como variable principal en la aplicación de los
métodos de interpolación usados y como segunda variable de control son usadas la o las
ecuaciones que describen la tendencia lineal identificada por región para el universo de
pozos con registros “Sónicos” y “Check Shot”.
Los métodos de interpolación como DSI y Kriging con Deriva Externa respetan las
variables estadísticas de los datos. Además generan una aproximación de la variación de
velocidad, la cual es una expresión suavizada que se logra extrapolando los valores de
velocidad.
49

CAPÍTULO 4
RESULTADOS Y ANALISIS DE RESULTADOS
4.1 Estudio y Control de calidad de Registros
Para el estudio y realización del modelado 3D fueron cargados en la plataforma GOCAD
un total de pozos, de los cuales, 43 contaban con registro Sónico, 56 con registro “Check
Shot”. Además de los pozos con registros, también fueron cargados al modelo varios pozos
que sólo contenían marcadores correspondientes a las regiones de interés identificadas
previamente en los registros.
Los datos de registros de pozo no sólo correspondían al Campo Santa Bárbara. También
fueron incluidos 15 pozos del Campo Tácata con sus respectivos marcadores geológicos.
Dentro de las curvas de velocidad es posible identificar cambios de tendencia asociados a
los cambios litológicos. En el presente estudio se identificaron tres y cuatro tendencias
principales dependiendo de la ubicación geográfica del pozo.
En las curvas de velocidad pertenecientes a pozos ubicados al norte del área de estudio
fueron definidos cuatro paquetes litológicos o regiones de interés (Figura 4-1 (a)): Somero
(Fm. Morichito), Alóctono, Carapita y Yacimiento. Para pozos ubicados al sur del área de
interés se establecieron como paquetes litológicos o regiones principales: Somero (Fm. La
Pica), Carapita y Yacimiento (Figura 4-1 (b)).

Figura 4-1 Paquetes litológicos identificados en el área norte (a) y área sur (b) del campo
Santa Bárbara/Pirital.
Las tendencias pueden ser comparadas con los modelos estratigráficos del campo. La
región denominada Somero, ubicada por encima del primer cambio de velocidades,
corresponde a las Formaciones Mesa, las Piedras, La Pica o Morichito cuyas profundidades
promedio se ubican entre los 0 y 8700 pies. Este cambio se asocia directamente con el tope
de Carapita (al sur) y tope del Bloque Alóctono (al norte). En las figuras 4-2 y 4-3 se dan
ejemplos de la región somera para dos pozos al sur y al norte del campo.
El siguiente cambio de tendencia, para los pozos al sur del área de estudio, se asocia al
tope de la formación Naricual que en promedio se ubica a partir de los 15000 pies de
profundidad. La región donde se incluye Naricual fue denominada Yacimiento y
comprende a partir del tope de esta formación en adelante. Para los pozos al norte sucede
que el segundo cambio de tendencia en la curvas de velocidad corresponde al tope de
51

Carapita que en esta zona se encuentra entre los 14,15 y 18 mil pies hasta 21500 pies (ver
Figuras 4-2 y 4-3)
En estos mismos pozos se reconoció un tercer y último cambio de tendencia
correspondiente al tope de la región Yacimiento la cual se define a partir del tope de la
formación Naricual.
Figura 4-2 Regiones reconocidas en registro sónico y columna estratigráfica asociada para
el área sur del campo Santa Bárbara/Pirital
52

Figura 4-3 Regiones reconocidas en registro sónico y columna estratigráfica asociada para
el área sur del campo Santa Bárbara/Pirital
53

4.2 Región Somera
La primera región definida es Somera; cuenta aproximadamente con un espesor que va
desde los 0 pies hasta los 8700 pies. Se asocia directamente a la formación Mesa, las
Piedras, La Pica y/o Morichito, las dos últimas al sur y al norte del área, respectivamente.
De manera general, para la región Somera se esperaba obtener distintas tendencias por la
presencia hacia el sur de la formación la “La Pica” y al norte la formación Morichito. La
disparidad en los valores de velocidad se hace evidente al observar las curvas de
velocidades interválicas para los pozos en donde se identificó esta región. Esto se ve en la
Figura 4-4, donde las curvas de velocidad no tienen correspondencia unas con otras.
Figura 4-4 Registros Sónicos región somera. Notar como no existe correspondencia entre
todas las curvas de velocidad
La variación de la velocidad posiblemente se debe a la separación impuesta al campo por
la presencia del Bloque Alóctono de Pirital. En la zona norte se deposito Morichito y en la
sur La Pica. Principalmente se diferencian porque la Formación La Pica se depositó en una
cuenca con influencia marina y la Formación Morichito muestra un ambiente puramente
54

continental, por tanto difieren litológicamente teniendo en La Pica alternancia de lutitas y
arenas interlaminadas mientras que Morichito se compone de conglomerados heterogéneos
con grandes guijarros y bloques derivados de las formaciones sedimentarias Cretácicas de
la Serranía del Interior.
Con los registros agrupados en serie es posible calcular, por medio de una regresión
lineal, la ecuación que mejor se ajuste y describa la tendencia general de los datos de cada
serie. Las series deben estar compuestas por pozos cuyas velocidades sean similares para
los mismos intervalos. Luego la combinación de cada una de las líneas de tendencia
calculadas para cada una de las series resulta como una recta de ajuste muy general para las
tendencias de todos los pozos En las figura 4-5, 4-6, 4-7, 4-8, 4-9 y 4-10 se pueden ver las
series I, II, III, IV, V y VI respectivamente y los pozos que componen las series se reflejan
en la Tabla 4.3.
Tabla 4-1. Pozos con tendencia Somero
Pozos con tendencia en región Somero
SBC-04 SBC-35
Serie 1 SBC-10E Serie 4 SBC-39
SBC-11 SBC-49
SBC-22E SBC-131
Serie 2 PIC-25 Serie 5 SBC-136
PIC-8 SBC-138
PIC-28 SBC-37E
Serie 3 PIC-26 Serie 6
SBC-130
PIC-6E SBC-142
55

Figura 4-5 Velocidad interválica y curva de ajuste lineal para la Serie I
Figura 4-6 Velocidad interválica y curva de ajuste lineal para la Serie II
56

Figura 4-7 Velocidad interválica y curva de ajuste lineal para la Serie III
Figura 4-8 Velocidad interválica y curva de ajuste lineal para la Serie IV
57

Figura 4-9 Velocidad interválica y curva de ajuste lineal para la Serie V
Figura 4-10. Velocidad interválica y curva de ajuste lineal para la Serie VI
58

Además se hizo una comparación entre los valores de velocidad provenientes de los
registros sónicos y los derivados de los registros “Check Shot”. Con los registros “Check
Shot” resulta menos exacto el cálculo de una tendencia lineal debido a la poca densidad de
datos que se poseen, pero de todas formas con la nube de datos de velocidades del registro
Check Shot igualmente se realiza un ajuste lineal con el fin de comparar globalmente la
tendencia de los datos. En la grafica 4-11 se observa cómo difiere la recta de ajuste para las
curvas de registros sónicos (línea azul) de la recta de ajuste para los registros “Check Shot”
(línea roja). Esta diferencia es posible porque están juntos los registros del lado norte
(Morichito) y sur del campo (La Pica). Incluso existen registros, tanto sónico como Check
Shot, con velocidades más altas de las esperadas para la región Somera. Esto se debe a que
en algunos pozos donde existía dificultad para identificar los cambios de tendencia en la
curva se hizo uso de los topes y bases formacionales provenientes de la interpretación
estructural vigente del campo y estos definen las regiones igualmente pero incluyendo
velocidades altas que se apartan de la tendencia general esperada para la región.
PIC-8
PIC-10E
PIC-25
SBC-1E
SBC-2E
Tabla 4-2 Pozos con registro Check Shot región Somero
SBC-4
SBC-6E
SBC-7
SBC-8
SBC-9
Pozos con registro Check Shot de la región Somero
SBC-11
SBC-12
SBC-13
SBC-14
SBC-15
SBC-16
SBC-22E
SBC-25
SBC-26
SBC-29
SBC-34
SBC-36
SBC-37E
SBC-40
SBC-51
SBC-55
SBC-61
SBC-64
SBC-65
SBC-74
SBC-83
SBC-84
SBC-110
SBC-130
SBC-136
59

Figura 4-11. Comparación de curvas de registro Sónico (en azul) y Check Shot (puntos
rojos) en región Somera
Por esto, se tomó la decisión de trabajar con las ecuaciones calculadas para cada serie. De
esta forma se consigue un buen ajuste lineal, a pesar de que un par de series no alcanzan un
factor de regresión > 0.7, pero se acercan bastante. Las ecuaciones de ajuste lineal se
muestran a continuación con el coeficiente de Pearson correspondiente:
Serie I Vint = 0,5255*Z + 6127,8 R 2 = 0,9072 (Ec. 4.1)
Serie II Vint = 0,6343*Z + 7734,5 R 2 =0,6535 (Ec. 4.2)
Serie III Vint = 0,8151*Z + 8216,4 R 2 =0,7016 (Ec. 4.3)
Serie IV Vint = 0,4295*Z + 6596,7 R 2 =0,9011 (Ec. 4.4)
Serie V Vint = 0,62 * Z + 6364,3 R 2 =0,6771 (Ec. 4.5)
Serie VI Vint = 1,0076*Z + 6212,9 R 2 =0,8716 (Ec. 4.6)
60

donde V es la velocidad interválica en pies/s y Z es la profundidad en pies. La Serie II y
Serie V poseen valores por debajo del minino establecido, sin embargo, para el coeficiente
menor ésta diferencia es del 7,11% con respecto a R 2 = 0,7.
Por tanto, se puede decir que, las ecuaciones derivadas del ajuste lineal muestran buena
correlación con los valores de velocidad. Las mismas, son válidas dentro de los intervalos
definidos por la región Somero y serán usadas más adelante como deriva externa para
modelar las velocidades en la capa Somera.
61

4.3 Región Bloque Alóctono
El bloque Alóctono de Pirital es una estructura que confiere complejidad en el modelo
estructural del Campo Santa Bárbara/Pirital (Parra ,2006). El campo está cubierto en mayor
parte por el bloque Alóctono de Pirital que corresponde al pliegue con rumbo N75E,
asociado al corrimiento de Pirital que buza hacia el norte. La región Bloque Alóctono
aumenta su espesor hacia el norte, alcanzando un máximo de 20000 pies dentro del
volumen de estudio.
A partir de las curvas de velocidad derivadas de los registros sónicos se obtiene la
tendencia regional de la velocidad. Es decir, es posible modelar la variación de velocidades
interválicas para esta región aplicando una ecuación lineal para relacionar de esta forma la
velocidad con la profundidad. Esto se puede comprobar con los datos de velocidad
interválica del pozo PIC-12 presentados en la Figura 4-12.
Figura 4-12 Velocidad interválica pozo PIC-12 con curva de ajuste lineal
62

En total se contó con 17 registros sónicos y 29 registros Check Shot en los cuales se
identificó la región Bloque Alóctono. Los registros sónicos proporcionan una mejor
aproximación debido a la cantidad de datos disponibles.
Por otro lado, en relación a los valores de velocidad calculados a partir de los registros
Check Shot, a pesar de que se tiene menor cantidad de muestras, éstas se ajustan a una
tendencia lineal con una correlación > 0.7 (Figura 4-13). El resultado de esta regresión
lineal para los valores de velocidad de tiros de verificación para la región Bloque Alóctono
es:
Vint = 0.493 * Z + 9512.7 (Ec. 4.7)
donde V es la velocidad interválica en pies/s y Z es la profundidad en pies.
PIC-1E
PIC-2E
PIC-3E
PIC-4
PIC-8
Tabla 4-3. Registro “Check Shot” región Bloque Alóctono
Pozos con registro "Check Shot" para la región Bloque Aloctono
PIC-9
PIC-10E
PIC-11
PIC-14
PIC-15
PIC-16
PIC-20
PIC-21
PIC-22
PIC-25
PIC-27
PIC-28
SBC-6E
SBC-17
SBC-18
SBC-22E
SBC-33
SBC-37E
SBC-50
SBC-61
Figura 4-13 Pozos con “Check Shot” Bloque Alóctono
SBC-62
SBC-95
SBC-130
SBC-136
63

En la grafica 4-14 se observan las curvas de velocidad derivadas del registro sónico. Se
distingue en éstas variaciones de los valores de velocidad, picos y bajos, en una misma
curva. Al ver con detalle se evidencia patrones similares en todas las curvas por lo que no
se sugiere errores en la toma de registros.
Figura 4-14 Registro Sónico Bloque Alóctono mostrando picos abruptos y tendencia
general
Para el cálculo de la ecuación de ajuste, se agruparon las curvas en serie, de la misma
forma como se realizó con las curvas de la región Somera. Resultaron cuatro series, a las
cuales se le ajusto una recta de tendencia lineal. Luego estas ecuaciones lineales se juntan
en una serie, a la cual se le calcula una relación lineal correspondiente a la tendencia
general de las curvas de velocidad, derivadas de los registros sónicos.
En las figura 4-15 se muestran las series formadas por los registros de los pozos
presentados en la tabla 6.
64

Tabla 4-4.Pozos con tendencia en región Bloque Alóctono
Figura 4-15 Curvas de velocidad interválica con ajuste lineal para la series I, II, III y IV de
la región Bloque Alóctono.
65

Las ecuaciones que resultan de la regresión lineal hecha a las series de la tabla 6 son:
Serie I Vint = 0.5634 * Z + 9604.7 R 2 = 0.8185 (Ec. 4.8)
Serie II Vint = 0.3297 * Z + 11136 R 2 = 0.7338 (Ec. 4.9)
Serie III Vint = 0.4377 * Z + 10062 R 2 = 0.6405 (Ec. 4.10)
Serie IV Vint = 0.4887 * Z + 8862.1 R 2 = 0.6635 (Ec. 4.11)
donde V es la velocidad interválica en pies/s y Z es la profundidad en pies.
La Serie III y Serie IV poseen valores por debajo del minino establecido, sin embargo,
para el coeficiente menor ésta diferencia es del 9.28% con respecto a R 2 = 0,7. Combinando
las cuatro series se obtiene como función lineal de ajuste para las curvas de velocidad
derivada de los registros sónicos la siguiente ecuación:
Vint = 0.4549 * Z + 9916 (Ec. 4.12)
donde V es la velocidad interválica en pies/s y Z es la profundidad en pies (Figura 4-16).
Esta ecuación es válida dentro de la región Bloque Alóctono, ya que representa un ajuste
lineal promedio de las tendencias presentes en las curvas de velocidad.
66

Figura 4-16 Velocidades interválicas con ajuste lineal para región Bloque Alóctono.
Además de calcular la ecuación de tendencia lineal, se hizo una comparación entre las
tendencias mostradas en los registros sónicos y las tendencias expuestas por los registros
check shot de cada pozo. En la figura 4-17 se observa que la aproximación lineal calculada
para ambos casos es muy similar, la diferencia de pendientes de ubica en un 9,5%. De
forma general, esto sugiere que la ecuación 4-12 representa un buen ajuste lineal para la
distribución de velocidades, por lo que será usada para describir la variación en
profundidad y modelar las velocidades dentro del volumen de interés del campo Santa
Bárbara/Pirital.
67

Figura 4-17. Comparación registros sónicos y “Check Shot” región Blque Alóctono.
68

4.4 Región Carapita
La región Carapita está presente en todo el campo Santa Bárbara. Se ubica entre los 2000
y 19000 pies de profundidad alcanzando 21000 pies en la parte norte del campo. En primer
lugar se realizó un análisis de las tendencias presentes en los disparos de verificación de
esta región. En la figura 4-18 se puede ver cómo la mayoría de los valores de velocidad
interválica se agrupan entre los 8 y 12 mil pies. La ecuación lineal que describe la tendencia
de los valores de velocidad es:
Vint = 0.0943*Z + 8605.9 (Ec. 4.13)
donde V es la velocidad interválica en pies/s y Z es la profundidad en pies,
Figura 4-18 Pozos con Check Shot Tendencia Carapita
La ecuación lineal describe una tendencia de los datos bastante general, ya que los valores
no se ajustan del todo a la recta. A su vez, existe dispersión en los datos, lo que se pone en
evidencia al apreciar el coeficiente de correlación el cual tiene un valor de 0,5896 para la
aproximación lineal aplicada a los valores de velocidad del registro “Check Shot”. Se
69

podría decir que esta ecuación es válida entre los 4000 y 18000 pies de profundidad, pero el
factor de correlación no supera el límite de lo esperado como valido en el estudio.
Se intentó ajustar la tendencia con ecuaciones polinomiales de segundo y hasta tercer
grado, pero la correlación de los datos con la curva no mejoraba; por tanto se mantuvo la
aproximación lineal de primer grado como la más cercana y descriptiva de la tendencia
general.
También se analizó el comportamiento de las curvas de velocidad calculadas a partir del
registro sónico. Al estudiar el comportamiento de todas las curvas se observa una mayor
correlación de los datos con respecto a una relación lineal (Figura 4-19). Este fenómeno se
esperaba para la región ya que litológicamente la formación Carapita consiste casi
exclusivamente de lutitas y las areniscas pero de manera general son escasas. En algunos
pozos existen variaciones abruptas de la velocidad o comportamientos irregulares como
picos y cambios bruscos de velocidad en intervalos cortos de profundidad.
Figura 4-19 Curvas de velocidad interválica mostrando tendencia general y alta correlación
entre los datos
70

En la figura 4-20 se tiene una curva de velocidad de la región Carapita para un pozo del
área de estudio. El comportamiento de la velocidad en este caso es lineal, menos en la parte
más profunda de la capa donde se observan sinuosidades que están dentro del rango de
velocidades para Carapita pero es un comportamiento anómalo que se presenta en la
mayoría de los registros de esta capa. Es posible que estas anomalías del registro necesiten
un estudio a detalle a fin de obtener una prospección más exacta de lo que realmente sucede
con la formación Carapita en el Campo Santa Bárbara. La curva de velocidad del pozo
SBC-33 (Figura 4-20) muestra para la zona de Carapita un ajuste lineal con un R 2
=0.7212;
por tanto la correlación de los datos con la ecuación de primer grado es buena.
Figura 4-20 Velocidad interválica pozo PIC-12 con curva de ajuste lineal
71

Luego se procedió a agrupar los registros de pozo que contenían datos pertenecientes a la
zona de Carapita. Se dividió en series los registros con tendencias más similares y se
calcularon aproximaciones lineales para las series que al igual que el registro del pozo
SBC-33 describían tendencias lineales con coeficientes de correlación alto. Resultaron
cinco series de los 16 pozos en donde se identificó la región Carapita. (Tabla 7).
Tabla 4-5. Pozos con tendencia en la región Carapita
Luego de hacer la regresión lineal a las cinco series, se obtuvo una ecuación general
promediando las ecuaciones de cada serie. En las figuras 4-21, 4-22, 4-23, 4.24 y 4-25 se
pueden ver las series formadas con las curvas de velocidad y las aproximaciones aplicadas.
Figura 4-21 Velocidad interválica y curva de ajuste lineal para la Serie I
72

Figura 4-22. Velocidad interválica y curva de ajuste lineal para la Serie II
Figura 4-23. Velocidad interválica y curva de ajuste lineal para la Serie III
73

Figura 4-24. Velocidad interválica y curva de ajuste lineal para la Serie IV
Figura 4-25. Velocidad interválica y curva de ajuste lineal para la Serie V
74

Mediante un promedio de todas las ecuaciones, finalmente se obtuvo una ecuación que
representa la tendencia general de velocidades en la región (Figura 4-26). Esta se usará
como deriva externa para el modelado de velocidades en la región Carapita. La ecuación
obtenida para Carapita fue:
Vint = 0.0948 * Z + 8783.9 (Ec. 4.14)
donde V es la velocidad interválica en pies/s y Z es la profundidad en pies.
Figura 4-26. Registros Sónicos con tendencia en región Carapita
Esta ecuación es válida y puede ser aplicada en la región Carapita, que va desde 4000
pies hasta los 18000 pies de profundidad. Es posible aplicar la ecuación como deriva
externa. Además, en la figura anterior se nota cómo la tendencia lineal de las velocidades
en la región Carapita se ajusta de manera muy cercana a las curvas de velocidad. Esto se
debe a la similitud de las tendencias en los distintos pozos donde se reconoció la capa.
La ecuación 4.14 será usada para modelar las velocidades en la región Carapita del
modelo de velocidades realizado en Gocad.
75

4.5 Región Yacimiento
Se obtuvieron velocidades para la región Yacimiento a partir de 21 registros sónicos y 52
registros Check Shot. Las velocidades que presenta esta capa superan los 10000 pies/s. En
la figura XX se observa la curva de velocidades derivada del registro sónico para el pozo
PIC-08, típica de la región Yacimiento.
Figura 4-27. Curvas de velocidad tipo de la región Yacimiento
76

El problema que se presentó con la región Yacimiento fue que asignar una tendencia
lineal a las curvas de registro sónico de esta capa fue una tarea difícil, pues existe un alto
nivel de dispersión de los datos y no sólo de forma individual, sino que al agrupar las
curvas de velocidad es posible observar como difieren altamente unas de otra. Este
fenómeno impide la asignación de tendencias lineales o de primer grado, por lo que no se
propone un ajuste para esta capa (Figura 4-28).
Figura 4-28. Velocidades Región Yacimiento
Con las curvas de velocidad derivadas de los registros Check Shot sucede que al ser
menor la densidad de datos es posible aproximar a una tendencia lineal la variación real de
los valores de velocidad de la región. Esto se puede ver en la figura 4-29 del registro Check
Shot del pozo PIC-25 donde la región Yacimiento (en verde) fácilmente se le es aplicada
una aproximación por mínimos cuadrados de donde resulta una recta que ajusta una
tendencia para la región Yacimiento de este pozo.
77

Figura 4-29. Velocidades interválicas a partir de registro Check Shot del pozo PIC-25,
región Yacimiento de color verde.
Las velocidades de los registros Check Shot de manera individual, es decir, para cada
pozo, se ajustan mejor a una aproximación lineal por lo que son más confiables si se desea
hacer una propagación de los datos usando esta tendencia lineal como deriva externa.
Globalmente la relación entre las velocidades de los pozos en la región es baja. A diferencia
de la región Carapita, Yacimiento cuenta con menos cantidad de puntos de velocidad por
eso la dispersión en este caso aumenta (Figura 4-30)
78

Figura 4-30. Ajuste lineal para las velocidades derivadas del registro Check Shot Región
Yacimiento
El ajuste lineal calculado para las velocidades derivadas de los registros Check Shot de la
capa Yacimiento viene dado por la ecuación:
Vint = 0,7078 * Z + 2467.2 (Ec. 4.15)
donde V es la velocidad interválica en pies/s y Z es la profundidad en pies, a profundidades
mayores de 15000 pies. El ajuste de correlación de esta ecuación es muy bajo, aunado a
esto las curvas de velocidad derivadas del registro sónico muestran que no existe una
tendencia general que se pueda atribuir a la región y de forma individual sucede lo mismo.
Por este motivo para la región Yacimiento se sugirió un modelado a partir de los mismos
datos velocidad del pozo haciendo una propagación DSI.
79

4.6 Histogramas de velocidades
Adicional a las tendencias se obtuvieron para las regiones histogramas de velocidades
para tener un mejor control de la distribución de frecuencias de velocidades interválicas en
cada una de las regiones identificadas. De esta forma se puede hacer una comparación entre
las velocidades que son introducidas para realizar las interpolaciones y las velocidades que
terminan poblando cada región o volumen.
La distribución de las velocidades en la Región Somera tiene un comportamiento
moderadamente asimétrico o sesgado. El valor promedio de velocidad para Somero es de
8000 pies/s y posee una distribución total de velocidades que van desde los 6000 pies/s
hasta los 13000 pies/s. (Ver figura
Figura 4-31. Histograma de velocidades interválicas de la región Somera
En la región Bloque Alóctono, la frecuencia de distribución de las velocidades tiene un
comportamiento moderadamente asimétrico o sesgado nuevamente. El valor promedio de
velocidad para Bloque Alóctono está en el rango de 8000 y 9000 pies/s y posee una
distribución total de velocidades que van desde los 6000 pies/s hasta los 16000 pies/s. (Ver
figura 4-32)
80

Figura 4-32. Histograma de velocidades interválicas de la región Bloque Alóctono
Figura 4-33. Histograma de velocidades interválicas de la región Carapita
81

Al contrario de los anteriores, en la región Carapita la frecuencia de distribución de las
velocidades tiene un comportamiento moderadamente asimétrico o sesgado pero en este
caso es hacia la izquierda (figura 4-33). El valor promedio de velocidad para Carapita está
en el rango de 11000 pies/s y posee una distribución total de velocidades que van desde los
8000 pies/s hasta los 12000 pies/s.
Figura 4-34. Histograma de velocidades interválicas de la región Yacimiento
Por último tenemos el histograma de la región Yacimiento (Figura 4-34). En éste la
frecuencia de distribución de las velocidades tiene un comportamiento simétrico unimodal.
El valor promedio de velocidad para Yacimiento se ubica en 11000 pies/s y posee una
distribución total de velocidades que van desde los 10000 pies/s hasta los 18000 pies/s.
82

4.7 Modelo de Velocidades
GOCAD es un software que se opera en una estación de trabajo “SUN Microsystems”
bajo el sistema operativo SOLARIS. El modelo de velocidades se generó en GOCAD y este
proceso consta de varias etapas. La primera de ellas la constituye la carga de los datos de
pozo en la plataforma, empezando por las coordenadas de ubicación y desvío de los pozos,
marcadores de cambio de región o tendencia, hasta la inclusión de las propiedades o
registros de velocidad interválica, todos en el dominio de la profundidad.
Una vez cargados todos los datos es posible visualizar los pozos con cada una de sus
propiedades agregadas al proyecto (Figura 4-35).
Figura 4-35. Pozos del área de estudio visto en ambiente GOCAD. En rojo la desviación de
los pozos.
Los marcadores de los pozos se interpretaron a partir de las curvas de velocidad derivadas
de los registros sónicos, porque es el registro muestreado con mayor frecuencia, cada pie;
en cambio los registros Check Shot tienen mediciones cada 500 pies o más, por lo que
definir un cambio de tendencia es complicado con tan pocos datos.
83

Los cambios de tendencia sugieren cambios litológicos (ver Figura 4-1), por tanto en esos
puntos se ubicaban los marcadores. Además se cargaron pozos de los cuales solo se tenía
desviaciones y marcadores de la interpretación estructural. Estos marcadores también
fueron usados de control. Se comparaban los establecidos por cambios de velocidad y los
marcadores estructurales para corregir o mantener los ya establecidos. Finalmente cada
pozo contenía entre tres y cuatro marcadores dependiendo de su ubicación en el campo (ver
Figura 4-3 y 4-4).
Seguidamente fueron cargados al proyecto los horizontes provenientes de la
interpretación sísmica estructural de Bosset J. y Gou Y., 2004. Estos fueron la base del
Corrimiento de Pirital, base del Corrimiento de Tala, base de Morichito, tope de Carapita y
tope de Naricual. Los marcadores, junto a las superficies que fueron cargadas como puntos
en el proyecto, son los que luego permitirán modelar las regiones estructurales.
Luego de haber cargado todos estos datos se crea el VOXET de trabajo. Este es un
volumen que abarca el área y profundidad de interés para el estudio. Fueron definidos dos
VOXET en el proyecto. En principio se creó un volumen de trabajo llamado Voxet
“extenso” que agrupaba los pozos ubicados tanto en el campo Santa Bárbara/Pirital como
los del campo Tácata (Figura 4-36), con un área de 705 km 2
de extensión.
Tabla 4-6. Dimensiones de VOXET extenso
DIMENSIONES DE VOXET EXTENSO
Norte - Sur 1075490 - 1059360
Este - Oeste 432000 - 388300
Los VOXET no son simplemente cubos, sino que a su vez éstos se subdividen en
pequeñas celdas. Para el Voxet “Extenso” se definieron celdas de dimensión 200x200x100
dando un total de 4 millones de celdas (Figura 4-36).
84

Figura 4-36. Voxet de trabajo área extensa, definido por los pozos del campo Santa
Bárbara/Pirital y Tácata
El segundo VOXET definido es llamado Voxet de “área reducida”; cuenta con una
extensión de 384 km 2
(Figura 4-37). En éste se ubican pozos del campo Santa
Bárbara/Pirital. Con esta configuración se reduce la incertidumbre ocasionada por la
ausencia de datos de pozo al Noroeste del VOXET extenso.
Tabla 4-7. Dimensiones VOXET área reducida
DIMENSIONES VOXET ÁREA REDUCIDA
Norte – Sur 1076000 – 1060000
Este - Oeste 432000 - 408000
De igual forma que para el Voxet extenso, para el Voxet “área Reducida” se definió un
tamaño de celda de 150x150x100 haciendo un total de 2250000 celdas. La profundidad
máxima de los dos volúmenes de interés depende directamente del pozo más profundo que
conforme el estudio. Para ambos volúmenes la profundidad es de 21524 pies.
85

Figura 4-37. Vista en planta de Voxet “área reducida” con la ubicación de pozos
pertenecientes al volumen de interés.
No sólo la incertidumbre por falta de pozos al noroeste del área hizo necesaria la
separación de los volúmenes de trabajo. Para empezar, los campos a pesar de ser contiguos
cambian litológica y estructuralmente. Las curvas de velocidad típicas de Santa
Bárbara/Pirital difieren de las de Tácata debido a la existencia de una serie de corrimientos
y retrocorrimientos menores, que dan origen a las diferentes repeticiones de las
formaciones Capaya y Carapita. Incluso la formación Naricual se encuentra a mayor
profundidad en el campo Tácata que en Santa Bárbara. (Bosset y Gou, 2004)
Las curvas de velocidad interválica son bastante diferentes para ambos campos; en Santa
Bárbara las tendencias de velocidad se comportan de manera relativamente constante en
intervalos amplios de profundidad, lo que permite establecer paquetes o regiones litológicas
principales en el campo. Mientras que en Tácata la complejidad estructural añade cambios
bruscos de velocidad y repeticiones (Figura 4-38).
86

Figura 4-38. Comparación de velocidades derivadas de registros sónicos del pozo SBC-35,
curva típica de campo Santa Bárbara, y del pozo TAG-11, curva típica de campo Tácata.
Para los pozos de Tácata también se contaba con marcadores provenientes de la
interpretación estructural (Figura 4-39). Lo único que no se obtuvo para esta zona fueron
los horizontes de la interpretación sísmica; sólo el corrimiento de Tala, una prolongación
del corrimiento de Pirital y el tope de Carapita fue lo que se logró concretar gracias a los
marcadores de los pozos de Tácata.
87

Figura 4-39. Curva de velocidades interválicas con marcadores del pozo TAG-11
Los horizontes cargados como puntos en la plataforma GOCAD permiten junto a los
marcadores de los pozos, crear las superficies de interés con mayor control. Por ejemplo el
tope de la región Yacimiento coincide con el horizonte Naricual Superior. Lo importante
del uso de los horizontes es que contribuyen en la construcción de superficies de manera
más exacta. Igualmente se tomaron otros horizontes para la creación de las superficies.
Las superficies (Ver figura 4-40) fueron modeladas haciendo uso de una herramienta de
Modelado Estructural que pone a disposición GOCAD. La creación de superficies se lleva a
cabo con los horizontes interpretados y luego se realiza una interpolación y se ajusta a los
marcadores de pozo las superficies. Es así como se generan superficies robustas que honran
88

la interpretación sísmica y el dato duro que son los marcadores provenientes de la
interpretación de los datos de pozo.
Figura 4-40 Superficies generadas en GOCAD dentro del VOXET “área reducida”
Es importante aclarar que la construcción del segundo Voxet se llevó a cabo por la
necesidad de proyectar las superficies de manera que abarcaran un volumen completo. A
pesar que las superficies son elementos independientes del “Voxet”, para la creación de las
regiones es necesario que éstas se extiendan en el volumen de interés, es decir que corten el
volumen. Para el “Voxet Extenso” no se pudo lograr esto por la falta de información. Y es
89

entonces cuando se construye el Voxet “área reducida” que justamente ocupa el área del
campo Santa Bárbara/Pirital.
Ahora bien, las superficies creadas definen los límites de cada región y haciendo uso de
ellas es posible generar los volúmenes asociados a las bases y topes geológicos. Las
regiones fueron creadas dentro del Voxet “área reducida” por lo antes mencionado. Para el
estudio se elaboraron cuatro regiones (Ver figura 4-41); la primera llamada Región Somera
que va desde el suelo, profundidad cero, hasta el tope del corrimiento en la zona norte del
volumen y hasta el tope de Carapita en la zona sur del volumen; la Región Bloque Alóctono
que tiene como tope la base de Región Somera y como base el Corrimiento de Pirital; la
Región Carapita cuyos límites son la base de la Región Somera y el tope de Naricual y por
último la Región Yacimiento que va desde el tope de Naricual Superior hasta los 21000
pies de profundidad, que es la máxima profundidad a la que llegan los pozos.
Figura 4-41 Regiones definidas dentro del VOXET “área reducida”
90

Luego de haber definido los parámetros anteriores, comienza el proceso de modelado de
velocidades mediante la propagación de los datos usando distintas técnicas geoestadísticas
que pone a la disposición GOCAD.
Para comenzar fue generado un modelo de velocidades para el Voxet “Extenso”, sin
definir en él regiones ni superficies, por medio de la interpolación DSI (Discrete Smooth
Interpolation). Esta técnica de interpolación se aplica utilizando los datos de pozo, es decir,
tomando en cuenta la propiedad Vint creada para cada pozo. La técnica de interpolación
DSI, como su nombre lo indica, es una técnica de atenuación; el problema para el modelo
generado de esta forma es que en las zonas donde no existen registros de pozos se asigna
igualmente valores muy suavizados para los cuales no se tiene control (Figura 4-42)
De igual manera se construyó un segundo modelo, pero esta vez cortando el VOXET con
las superficies existentes. Con esto se logra visualizar la forma en que el volumen es
cortado por las superficies (Ver figura 4-42). Además se puede ver cómo las regiones
dentro de la región Bloque Alóctono aumentan las velocidades en promedio al resto de las
regiones. Esta interpolación permite observar una distribución de velocidades promedio.
Figura 4-42 Modelo de velocidades en VOXET Extenso. Cortado por las superficies
91

En general el resultado es lo que se esperaba con la aplicación del método DSI, un
modelo suavizado con velocidades que aumentan en función de la profundidad
gradualmente en la región Somera y al entrar en la Región Alóctono las velocidades se
elevan de forma vertiginosa con la profundidad. Para la región Carapita las velocidades
incrementan de manera paulatina, y se observan, al igual que en todas las regiones, zonas
donde se concentran velocidades altas que son producto de las variaciones laterales y de
cambios de velocidad bruscos con respecto a los pozos vecinos. En último lugar para la
Región Yacimiento se observa un incremento de velocidad alcanzando zonas de hasta
18000 pies/s; esto sucede en las zonas más profundas de la Región Yacimiento.
El haber cortado con las superficies el VOXET permite captar los cambios de tendencias
de una manera más sencilla ya que es justo en estos límites donde coinciden los cambios de
velocidad y los cambios litológicos. Se puede entonces decir que haciendo uso del método
DSI con los datos de pozo como valores de control se obtiene un modelo global del cambio
aunque no completamente confiable en las zonas donde no se tienen pozos de control y la
interpolación realiza aproximaciones que podrían no ser las más idóneas.
Para la elaboración del modelo de velocidades por región se utilizó como métodos de
interpolación el Kriging con Deriva Externa, el cual es un método geoestadístico de
estimación de puntos que utiliza un modelo de variograma como función de ponderación
para la obtención de datos. Es un procedimiento simple y eficiente que incorpora una
segunda variable para la estimación de una primaria (Deutsch y Journel, 1998). Es una
extensión del Krigeaje con modelo de tendencia o Krigeaje Universal, el cual honra las
tendencias locales de los datos e integra condiciones de universalidad suplementaria
relativa a una o varias variables externas medidas de forma exhaustiva en todo el dominio
donde se desea estimar la variable de interés.
Los valores de velocidad presentes en cada pozo son la variable de interés primaria, la
cual es escasa por la distribución de pozos en el VOXET. La mayor concentración de datos
proviene justamente de las ecuaciones halladas como tendencias lineales de velocidad para
cada región, la variable secundaria o deriva externa, y de los variogramas que atribuyen el
92

comportamiento espacial de la variable primaria sobre un área, es decir, cambios laterales
de la velocidad en el VOXET.
Se había discutido anteriormente (Capítulo 4., sección 4.2) que para la Región Somera las
tendencias de las curvas de velocidad de los pozos presentaban linealidad, pero ésta no se
correspondía a una sola ecuación lineal aplicable entre los pozos, sino que la variación de
las curvas hacia necesaria la aplicación de regresiones lineales pero por series o
agrupamiento de curvas de rangos de velocidad similares. Así se obtuvieron seis ajustes
buenos que definieron de forma general las tendencias presentes en la región.
Para la región Somero la estrategia de modelado fue ligeramente distinta a las demás
regiones a la que se les aplicó la interpolación Kriging con deriva externa. Se creó una
propiedad llamada “Vint_ec_Somera” en cada uno de los pozos donde se reconoció el
intervalo somero; mediante un “script” (código de programación) se fijo la ecuación a la
propiedad y se crearon rectas de velocidad a partir de la aproximación lineal calculada para
cada serie de pozos. Esta propiedad se propagó en la región con el método DSI propio de
GOCAD y esta misma es la que se usará como variable secundaria en la interpolación por
Kriging con Deriva Externa de la región Somera.
El modelado para las regiones Carapita y Alóctono contó con la creación de una
propiedad a la cual se le aplicó igualmente un “script” para interpolar por la región la
ecuación de tendencia lineal y así usarla como variable secundaria al realizar Kriging con
Deriva Externa.
Finalmente la región Yacimiento se modeló con el método DSI debido a la ausencia de
tendencia lineal general que describiera la variación de la velocidad en esta región.
En la figura 4-43 podemos ver un volumen que integra todas las regiones que conforman
el campo Santa Bárbara en profundidad. Se puede ver cómo las velocidades cambian en
función de la profundidad de forma que para la región Somera (figura 4-44) corresponden
velocidades bajas con variaciones laterales sutiles producto posiblemente de cambios
litológicos, ya que dentro de esta región se ubican las formaciones Mesa-Las Piedras, La
Pica y Morichito respectivamente. Estas formaciones desde un punto de vista general
93

consisten en areniscas, calizas, arcillas, lutitas y conglomerados constituidos por peñones,
peñas y gravas intercaladas con arenas y limolitas.
En la figura 4-45 vemos la distribución de las velocidades en la región Bloque Alóctono.
El alóctono está compuesto de rocas de edad Mesozoico y Cenozoico. Las formaciones que
le componen, en particular la calizas de las formaciones Chimana y El Cantil del Cretáceo y
las rocas más antiguas encontradas en el bloque alóctono, las “capas rojas” de la Formación
La Quinta (edad Jurásico-Triásico), ocasionan el aumento de la velocidad en esta región.
(Bosset y Gow, 2004). Se puede ver cómo las velocidades aumentan con respecto a las
regiones vecinas encontrándose velocidades que van desde 12000 pies/s hasta 19000 pies/s,
aproximadamente.
Las velocidades que se observan en la región Carapita son mayores a las de la región
Somera pero menores a las presentes en la región Bloque Alóctono (Figura 4-46). Las
velocidades aumentan lentamente para esta región y se ubican entre los 8000pies/s y 12000
pies/s. Existen variados contrastes laterales que pueden ser asociados a la presencia
eventual de lenguas de areniscas interestratificadas con las lutitas marinas de Carapita
según González de Juana et al (1980).
Finalmente en la región Yacimiento se puede observar el incremento de la velocidad de
forma rápida. Las velocidades interválicas en esta región se ubican entre los 10000 pies/s y
los 20000 pies/s. Es un rango bastante amplio y las razones por las que se obtuvo este
resultado después de hacer la propagación, es que posiblemente se hayan incluido para
prospectar la región Yacimiento velocidades bajas pertenecientes a la región Carapita
(Figura 4-47). Otra razón es que no son muchos los pozos que tenían registros por debajo
de la región Yacimiento; por lo tanto la estimación de las velocidades en la región
permanecen con alta incertidumbre.
En la Figura 4-43 se muestra el modelo 3D de velocidades obtenido después de realizar
las propagaciones pertinentes en cada región.
94

Figura 4-43. Modelos 3D de Velocidades Interválicas Campo Santa Bárbara

Figura 4-44. Modelo 3D de Velocidades Interválica Campo Santa Bárbara. Velocidades
para la Región Somera

Figura 4-45. Modelo 3Dde Velocidades Interválica Campo Santa Bárbara. Velocidades
para la Región Bloque Alóctono.
97

Figura 4-46. Modelo 3Dde Velocidades Interválica Campo Santa Bárbara. Velocidades
para la Región Carapita.
98

Figura 4-47. Modelo 3D de Velocidades Interválica Campo Santa Bárbara. Velocidades
para la Región Yacimiento.
99

Entre las formaciones del bloque alóctono (cuyas velocidades varían entre 12000 y 19000
pies/s) y las lutitas de la Formación Carapita que se encuentran por debajo (con velocidades
interválicas entre los 8000pies/s y 12000 pies/s), existe una inversión importante de la
velocidad (Figura 4-43).
Esta inversión se traduce como un alto índice de impedancia, lo cual trae como
consecuencia que las ondas sísmicas se reflejen sobre el plano de falla constituido por el
corrimiento de Pirital, produciendo reflexiones múltiples que perturban la señal en los
niveles más profundos (Figura 4-48).
Además, debido al espesor del bloque alóctono (hasta 17000-18000 pies), la energía de
las ondas sísmicas pierde intensidad, bien sea por absorción cruzando todo el bloque
(Figura 4-49); o por dispersión (dentro del bloque alóctono las ondas se reflejan debido a
los contrastes de impedancia fuertes). Por ende, debajo del bloque alóctono unas áreas
pueden aparecer transparente sísmicamente (zona sorda) (Figura 4-50).
100

Figura 4-48. Comparación del modelo de velocidades con la interpretación sísmica I.
Presencia de impedancias importantes entre la base de la región Bloque Alóctono y el tope
de la región Carapita. (Modificada de Bosset J. y Gou Y. 2004)
101

Figura 4-49. Comparación del modelo de velocidades con la interpretación sísmica II.
Presencia de impedancias importantes entre la base de la región Bloque Alóctono y el tope
de la región Carapita. (Modificada de Bosset J. y Gou Y. 2004)
102

Figura 4-50. Influencia del alóctono de Pirital sobre la calidad de la sísmica
(Tomado Bosset J. y Gou Y. 2004)
103

4.8. Histogramas de Velocidades (Modelo Final)
Luego de haber realizado el modelado 3D de velocidades interválicas de onda
compresional en el campo Santa Bárbara se generaron histogramas de las velocidades
asignadas por los algoritmos de propagación en cada región. Esto nos permite garantizar
que los datos de salida corresponden a datos reales de entrada existentes en la zona y que
son empleados en la realización de este trabajo.
Figura 4-51. Histograma de velocidades del modelo final. (a) Histograma región Somero,
(b) Histograma región Bloque Alóctono, (c) Histograma región Carapita.
Para el histgrama de la región Somera (Figura 4-45 (a)) se tiene un comportamiento
simétrico. El valor promedio de velocidad para Somero es de 8340 pies/s y posee una
distribución total de velocidades que van desde los 4700 pies/s hasta los 13500 pies/s
En el histograma para la región Bloque Alóctono (Figura 4-45 (b)) la frecuencia de
distribución de las velocidades tiene un comportamiento simétrico. El valor promedio de
104

velocidad para Somero en el rango de 12000 y 13000 pies/s y posee una distribución total
de velocidades que van desde los 7200 pies/s hasta los 18000 pies/s. Estas velocidades
difieren bastante de las analizadas por el histograma aplicado a las curvas de velocidad de
los pozos de la región Bloque Alóctono. La causa de esto es que para el modelo se añadió
una ecuación lineal, es decir, el bloque alóctono a pesar de ocupar gran parte del campo
posee hacia el norte un vacio de valores de velocidad para control, en esta zona la ecuación
lineal domina la propagación y arroja valores altos de velocidad.
En histograma para la región Carapita (Figura 4-45 (c)) la frecuencia de distribución de
las velocidades tiene un comportamiento moderadamente asimétrico. El valor promedio de
velocidad para Carapita se ubica en 10300 pies/s y posee una distribución total de
velocidades que van desde los 8000 pies/s hasta los 16000 pies/s.
Finalmente en la Figura 4-52 se observa el histograma para el modelo de velocidades
final generado en el Voxet “área reducida”. Se distingue un comportamiento simétrico con
un alto porcentaje de valores entre los 9000 y 12000 pies/s, lo que principalmente se debe a
que este rango de velocidades es el más frecuente en la zona.
Figura 4-52 Histograma de velocidades para el modelo final de velocidades interválicas
generado en el Voxet “área reducida”
105

CONCLUSIONES
Para la realización del modelado 3D de velocidades interválicas de onda compresional a
partir de registros de pozo, en el campo Santa Bárbara se incluyeron múltiples disciplinas
de las geociencias empleadas en los estudios integrados de yacimiento. Las técnicas de
caracterización del subsuelo incorporadas en el estudio fueron los registros de pozo,
sísmica de pozo y métodos geoestadísticos.
Es de gran importancia el control y validación de los datos de pozo para la creación del
modelo de velocidades, de esta forma se garantiza un resultado similar al comportamiento
real de las velocidades en el área de estudio.
Se pudo comprobar que mediante el análisis de los valores velocidad interválica
derivados de los registros sónico y “Check Shot” es posible definir regiones conformadas
por paquetes litológicos, donde la propiedad de velocidad interválica mantiene tendencias
similares. Es posible estimar para estas tendencias curvas de ajuste, por medio de las cuales
se generará el modelado de las velocidades.
A partir de las velocidades interválicas calculadas para cada pozo y la estratigrafía, el
campo Santa Bárbara puede subdividirse a grandes rasgos en cuatro regiones asociadas a la
litología, la primera región denominada Somera está conformada por las formaciones Mesa,
Las Piedras, Morichito (al norte del área) y La Pica (al sur del área). La segunda región
llamada Bloque Alóctono la componen las formaciones Vidoño, San Juan, San Antonio,
Querecual, Chimana El Cantil y formación Barranquín. La tercera región definida es
Carapita constituida por la formación del mismo nombre. Y por último se definió
Yacimiento compuesta por la formación Naricual, principal productora y reserva de
hidrocarburos en los campos profundos del norte de Monagas.
La región Somera identificada en las curvas de velocidad derivadas de los registros
sónicos y “Check Shot” poseen una relación lineal; sin embargo, no existe correlación total
entre las estimaciones obtenidas para cada curva de velocidad. Por esto se agrupó en series
la propiedad de pozos que tuvieran tendencias de velocidad similares y para estas series fue
estimada la ecuación de ajuste lineal.
106

Las tendencias encontradas en las curvas de velocidad derivadas de los registros sónicos y
“Check Shot” para las regiones Carapita y Bloque Alóctono, poseen una relación lineal
cuya curva de ajuste es una ecuación de primer grado, que describe la relación entre las
variables dependiente (la velocidad) y la variable independiente (la profundidad) de las
tendencias.
La región Yacimiento a diferencia de las otras tres regiones no presentó una relación
lineal ni de otro tipo con la cual pudiera ser ajustada la variación de la velocidad interválica
de la región. Por tanto se hará uso exclusivamente de los datos de pozo para la construcción
del modelo
El uso de GOCAD para la elaboración del modelo, significó un reto desde el punto de
vista operacional, a pesar de ser una herramienta bastante amigable al usuario. Son
múltiples las aplicaciones brindadas por el software en la integración de técnicas de
prospección geocientífica.
La región Somera presentó variaciones de velocidad interválica entre los 5000 y 10000
pies/s, las velocidades para la región Bloque Alóctono se ubicaron entre los 12000 y 18000
pies/s, mientras que para Carapita se estipularon velocidades que van de los 8000 a los
12000 pies/s. Por último, para la región Yacimiento se establecieron velocidades entre los
11000 y 18000 pies/s.
La comparación del modelo de velocidades con perfiles sísmicos del área permite
corroborar los efectos provocados por la litología y la complejidad estructural del campo.
Finalmente con el modelado de velocidades interválicas se dan a conocer de manera más
exacta las variaciones horizontales y verticales presentes en el campo, las cuales
contribuirán más adelante en la aplicación de migración pre-apilamiento de los datos
sísmicos llevada a cabo como parte de la reinterpretación pautada para los levantamientos
del área.
107

RECOMENDACIONES
Para futuros trabajos en el área o en trabajos con metodología similar para la creación de
modelo de velocidad se recomienda comparar las anomalías presentes en los registros
sónicos y “Check Shot” con registros que midan otra propiedad con la finalidad de
corregirlos y garantizar que los valores usados sean correctos y no errores de medición.
Además se propone utilizar velocidades de procesamiento para comparar con las
obtenidas en el modelado. No incluir estas velocidades como variables secundarias en el
modelado ya que la idea es construir un modelo independiente a estas para mejorar los
problemas presentados en procesamientos anteriores.
También se exhorta a usar este modelado de velocidades como guía al momento
establecer velocidades de procesamiento para realizar migración pre- apilamiento. Esto con
el fin de verificar la configuración geométrica y ubicación de los principales eventos
sísmicos que sustentan los modelos estructurales del campo.
Asimismo se sugiere incluir nuevos modelos estructurales del campo Santa Bárbara y
Tácata con mayores horizontes interpretados para controlar y modelar los cambios de
velocidad tomando más en cuenta la litología y estructura ya que a pesar de ser un área muy
madura desde el punto de vista exploratorio aun existen estructuras muy atractivas, que
representan una nueva oportunidad exploratoria.
108

BIBLIOGRAFIA
Bastidas, C. (2008). Modelado Geofísico Integrado en el norte de Monagas, aplicando la
metodología de Interpretación Sísmica 2D e Inversión Sísmica. Trabajo Especial de Grado
para optar al título de Ingeniero Geofísico. Departamento de Ciencias de la Tierra,
Coordinación de Ingeniería Geofísica, Universidad Simón Bolívar. Caracas, Venezuela.
Bolívar, M. (2005). Reporte: Análisis de incertidumbre de interpretación sísmica de
Pirital. Exploración y Producción PDVSA. Estudios Integrados Distrito Norte. Puerto La
Cruz, Venezuela.
Bosset J. y Gou Y. (2004). Campo Santa Bárbara – Pirital. Revisión Modelo Estructural.
Exploración y Producción PDVSA, Estudios Integrados Norte de Monagas. Puerto La
Cruz, Venezuela.
Carrillo, M., F. Ruíz, R. Herrera, M. Vargas, H. Oropeza, (2005) Modelo Sedimentológico y
Estratigráfico detallado del Terciario y el Cretácico del Campo Santa Bárbara-Pirital,
Norte de Monagas. INT-10657,2005. Los Teques. PDVSA Intevep.
Castillo, C., R. García, G. Gedler, R. Porjesz, O. Rondon, (2000) Introducción a la
Geoestadística petrolera. Manual, Los Teques, INTEVEP
Cobos S. (2002). Interpretación estructural en el cinturón plegado de Monagas, por medio
de la integración de datos sísmicos, datos de pozo y geología de superficie. Trabajo
Especial de Grado para optar al título de Ingeniero Geofísico. Departamento de Ciencias de
la Tierra, Coordinación de Ingeniería Geofísica, Universidad Simón Bolívar. Caracas,
Venezuela.
Dewam, J. (1984). Modern Open-hole Log Interpretation. Pág 139-142
Di Croce J., Bally A. y Vail P. (1999). Sequence stratigraphy of the eastern venezuelan
basin. En P. Mann (Ed), Caribbean basins. Sedimentary basins of the world (pp. 419 –
476). Amsterdam: Elsevier Science.
109

Díaz, M. (2006). Modelado 3D de velocidades interválicas de onda compresional
utilizando registros de pozos, del Campo Carito, Edo. Monagas, Venezuela. Trabajo
Especial de Grado para optar al título de Ingeniero Geofísico, Escuela de Ingeniería
Geofísica, Universidad Central de Venezuela. Caracas, Venezuela.
Dix, C., (1955). Seismic velocities from surface measurements: Geophysics. 20 (1): 68-86.
Dobrin, M., (1961). Introducción a la prospección geofísica, Ediciones Omega, S.A.,
Barcelona, España. 483 p.
GOCAD Suite 2.5 Guía de Usuario (2008). Paradigm Geophysical Corp.
González de Juana, C; J. Iturralde, y X. Picard, (1980), Geología de Venezuela y de sus
Cuencas Petrolíferas, Tomo II, Ediciones FONINDES, Caracas-Venezuela.
Morales, M., N. Santiago, J. Cordova, T. Maita, H. Rojas, M. Moya. Informe Final Área
Tácata-Tacat (2006). Exploración y Producción PDVSA. Puerto La Cruz, Venezuela.
Parnaud, F., I. Truskowski, Y. Gou, O. Gallango, J. Pascual, F. Roure, J Di Croce, y H.
Passalacqua, (1991) Modelo geológico integrado del transecto Chacopata – Uverito.
Cuenca Oriental de Venezuela. INT-02337,91. Los Teques. PDVSA Intevep.
Parra, M. (2006). Modelado Estructural y Restauración de la Región Noroccidental de la
Subcuenca de Maturín. Trabajo Especial de Grado para optar al título de Magíster en
Ciencias de la Tierra. Departamento de Ciencias de la Tierra, Coordinación de Ingeniería
Geofísica, Universidad Simón Bolívar. Caracas, Venezuela.
Repilloza, H. (2008). Caracterización Geofísica de la formación Naricual en el norte de
Monagas, aplicando la metodología de Inversión Sísmica. Trabajo Especial de Grado para
optar al título de Ingeniero Geofísico. Departamento de Ciencias de la Tierra, Coordinación
de Ingeniería Geofísica, Universidad Simón Bolívar. Caracas, Venezuela.
110

Roure F., E. Hernandez, N. Bordas-Lefloch, C. Rivero, J. Toro, C. Bourg, W. Sassi, N.
Guilhaumoul, S. Lecornec-Lance, (2003). Petroleum Systems and Reservoir Appraisal in
the Sub-Andean Basins (Eastern Venezuela and Eastern Colombian Foothills). Memoria 79
AAPG, p. 750–775.
Royer, J. (2004). 3D Modeling and Visualization. 9 th International CODATA Conference.
Berlin, Alemania.
Summaa L., E. Goodmanb, M. Richardsona, I. Nortonb, A. Greenb, (2003). Hydrocarbon
systems of Northeastern Venezuela: plate through molecular scale-analysis of the genesis
and evolution of the Eastern Venezuela Basin. Marine and Petroleum Geology 20 (2003), p.
323–349. Amsterdam: Elsevier Science.
Western Atlas (1992), Introduction to Wireline Log Analysis. P. 142
Yoris F. y Ostos M. (1997). Geología de Venezuela. En Schlumberger (Ed.), Venezuela
WEC 1997, evaluación de pozos, p. 3-44. Caracas: Editores.
111