Post on 27-Jun-2022
1
Facultad de Ciencias
Departamento de Geociencias
Interpretación de Geoformas Paleozoicas en una Localidad de los
Llanos Orientales en Colombia aplicando Atributos Sísmicos.
Proyecto de grado presentado por:
Laura Carolina Esquivel Arenales
Directora:
Indira Molina PhD.
Bogotá DC, Colombia
Mayo 2020
____________________ ____________________
Firma Estudiante Firma Directora
2
AGRADECIMIENTOS
Este trabajo fue posible con la ayuda y apoyo de varias personas e instituciones. A mi
directora Ph.D Claudia Indira Molina por el apoyo a lo largo del trabajo, al SGC por permitir
el uso de información necesaria para realizar este proyecto muy especialmente a la geóloga
Fanny Johana Villamizar, quien con su ayuda y enseñanza me guio en cada una de las etapas
de este trabajo. Agradezco a mi familia por el apoyo incondicional, a mi padre por brindarme
siempre ayuda con su conocimiento, a mi madre por su motivación a lo largo de este camino
y a mi hermano por siempre estar a mi lado. Por último, a grandes amigos que me
acompañaron por mi recorrido en la universidad, William Carantón, Daniel Castañeda,
Paula Cabrejo, Gabriela Rojas y Nicole Salazar. Sin ellos, nada de esto hubiese sido posible.
3
RESUMEN
El presente estudio se enfoca en producir mapas de geoformas paleozoicas en una región de
la cuenca de los Llanos Orientales. Este estudio es de gran importancia ya que aporta
resultados innovadores sobre el Paleozoico para el área de estudio, estos resultados lograron
determinar un ambiente deposicional Marino Profundo en el Cámbrico y Marino Somero en
el Ordovícico. Esto se obtuvo partiendo del análisis de secuencias sismoestratigráficas,
electrofacies y de estructuras geológicas de un cubo sísmico 3D, 4 líneas sísmicas 2D y 4
pozos. Se aplicaron los atributos de Coherencia, Varianza, Amplitud RMS y Descomposición
Espectral con el fin de generar mapas geomorfológicos correspondientes a las secuencias
Ordovicico Medio, Ordovicico Inferior y Cambrico.
Palabras Clave: Cuenca de los llanos Orientales, Símica, Registros de Pozo, atributos
sísmicos, Geoformas Paleozoicas.
ABSTRACT
The present study focuses on producing maps of Paleozoic geoforms in a region of the Llanos
Orientales Basin. The importance of this study relays in the Paleozoic innovative results
about the depositional environment which suggest deep marine depositional environment in
the Cambrian and a shallow marine environment in the Ordovician. This was done by
carrying out seismic stratigraphy, electrofacies and geological analysis of a 3D seismic cube,
four seismic lines and four wells. The attributes of Coherence, Variance, RMS amplitude and
Spectral Decomposition were applied in order to generate geomorphological maps
corresponding to the Middle Ordovician, Lower Ordovician and Cambrian sequences.
Key Words : Llanos Orientales Basin, Seismic, Well Logs seismic attributes, Paleozoic
Geoforms.
4
Tabla de Contenido
RESUMEN ................................................................................................................................... 3
ABSTRACT ................................................................................................................................ 3
1. INTRODUCCIÓN ............................................................................................................... 7
2. OBJETIVOS ........................................................................................................................ 8
2.1. Objetivo General .......................................................................................................... 8
2.2. Objetivos Específicos.................................................................................................... 8
3. MARCO GEOLÓGICO ...................................................................................................... 8
3.1. Evolución Tectónica y Estratigráfica ......................................................................... 8
4. DATOS ................................................................................................................................. 9
5. MARCO TEÓRICO .......................................................................................................... 10
5.1. Sísmica de Reflexión .................................................................................................. 10
5.2. Sísmica 3D .................................................................................................................. 10
............................................................................................................................................ 11
5.3. Estratigrafía de Secuencias ........................................................................................ 11
5.4. Sismoestratigrafía ....................................................................................................... 12
5.5. Registros de Pozo ....................................................................................................... 13
5.5.1. Gamma Ray ........................................................................................................ 13
5.5.2. Resistividad ......................................................................................................... 13
5.5.3. Densidad .............................................................................................................. 13
5.5.4. Sónico .................................................................................................................. 14
5.5.5. Densidad Neutrón ............................................................................................... 14
5.6. Análisis de Electrofacies con registros de Pozo. ...................................................... 14
5.7. Atributos Sísmicos ..................................................................................................... 15
5.7.1. Atributo de Coherencia ...................................................................................... 15
5.7.2. Descomposición Espectral .................................................................................. 15
5.7.3. Atributo de Amplitud RMS .............................................................................. 16
5.7.4. Atributo de Varianza ......................................................................................... 16
6. METODOLOGÍA .............................................................................................................. 16
6.1. Resolución Sísmica Vertical....................................................................................... 17
5
6.2. Correlación de Pozos con la Sísmica 2D y 3D ......................................................... 17
6.2.1. Sismograma Sintético ......................................................................................... 17
6.2.2. Correlación de pozos con la sísmica 2D ............................................................ 18
6.2.3. Interpretación de Horizontes Sísmicos .............................................................. 18
6.3. Interpretación de Fallas ............................................................................................. 18
6.4. Análisis de Electrofacies con registros de Pozo........................................................ 18
6.5. Análisis Secuencial ..................................................................................................... 18
6.6. Sismoestratigrafía ....................................................................................................... 19
6.7. Mapas estructurales ................................................................................................... 19
6.8. Aplicación de atributos sísmicos ............................................................................... 19
7. RESULTADOS Y ANALISIS........................................................................................... 20
7.1. Sismograma Sintético ................................................................................................. 20
7.2. Correlación de pozos con la sísmica 2D .................................................................... 21
7.3. Interpretación de fallas .............................................................................................. 22
7.4. Análisis de Electrofacies con registros de Pozo........................................................ 23
7.4.1. Pozo-1 .................................................................................................................. 24
7.4.2. Pozo-3 .................................................................................................................. 24
7.5. Análisis Secuencial ..................................................................................................... 25
7.6. Sismoestratigrafía ....................................................................................................... 27
7.6.1. Secuencia Sísmica C ........................................................................................... 28
7.6.2. Secuencia Sísmica D ........................................................................................... 29
7.6.3. Secuencia Sísmica E ........................................................................................... 30
7.7. Mapas Estructurales .................................................................................................. 31
7.8. Atributos Sísmicos ..................................................................................................... 31
8. INTERPRETACIÓN DE RESULTADOS ...................................................................... 37
9. CONCLUSIONES .............................................................................................................. 40
10. BIBLOGRAFÍA ............................................................................................................. 42
11. Anexo A .......................................................................................................................... 45
13.1. Marco Geológico ..................................................................................................... 45
13.2. Datos ....................................................................................................................... 48
6
13.3. Procesamiento Pre-Stack Time Migration (PSTM) ............................................ 48
13.5. Atributo de Coherencia.......................................................................................... 50
14. Anexo B: Complementos Metodología ......................................................................... 50
14.1. Resolución Sísmica ................................................................................................. 51
14.2. Sismograma Sintéticos ........................................................................................... 51
14.3. Correlación de la Información de Pozos con la Sísmica 2D y 3D ....................... 53
15. Anexo C: Resultados ..................................................................................................... 53
15.2. Dato Palinológico ................................................................................................... 56
................................................................................................................................................ 56
15.3. Análisis de Estratigrafía sísmica ........................................................................... 56
15.3.1. Secuencia A ......................................................................................................... 56
15.3.2. Secuencia B ......................................................................................................... 57
15.3.3. Secuencia F ......................................................................................................... 30
15.3.4. Mapas Estructurales ....................................................................................... 58
15.3.5. Atributos Sísmicos .............................................................................................. 59
15.3.6. Atributos Sísmicos a diferentes tiempos (cada 170ms) ................................... 61
7
1. INTRODUCCIÓN
La exploración sísmica es un método geofísico que consiste en el uso de una fuente, la cual
produce ondas acústicas que se propagan en un medio, esto con el fin de determinar la forma
y distribución de las diferentes unidades litológicas en profundidad (López et al., 2008). A lo
largo del tiempo, este método se ha convertido en una herramienta indispensable para
obtener información del subsuelo. En Colombia, se han realizado varios estudios con esta
finalidad, especialmente en la región oriental del país para las eras Cenozoica y Mesozoica
debido a el potencial hidrocarburífero existente en la región, Suarez y Solano, 2012.
Por otro lado, la era Paleozoica en Colombia ha sido poco estudiada, entre los estudios más
relevantes se encuentra el presentado por Martínez 2008, en el cual se destacan los diferentes
eventos geológicos de las cuencas subandinas en el Paleozoico. El presente proyecto pretende
contribuir a estos estudios con diferentes análisis partiendo de la información geofísica
proporcionada (sísmica 2D-3D y registro de pozos) en una región de los Llanos Orientales.
Esto fue posible debido a que los datos utilizados para este proyecto, entre ellos la sísmica,
fue especialmente procesada para resaltar elementos de formaciones Paleozoicas; sin
embargo, es importante mencionar que existieron limitantes como la resolución sísmica, la
cual dificultó la interpretación de diferentes elementos geológicos en profundidad.
Este estudio se desarrolló por medio de la interpretación de estructuras geológicas de un
cubo sísmico 3D; además se realizó el análisis de secuencias sismoestratigráficas, secuencial
y de electrofacies paleozoicas de la información geofísica proporcionada. Finalmente, se
aplicaron los atributos sísmicos de Coherencia, Varianza y Amplitud RMS, los cuales se
utilizan en la interpretación sísmica de estructuras, estratigrafías y litologías (Chopra y
Marfurt, 2007). La evaluación, integración y análisis de información, permitió generar mapas
geomorfológicos con información geológica sobre los ambientes de sedimentación en el área
de estudio, esto se determinó para las secuencias paleozoicas correspondientes a los periodos:
Ordovicico Medio, Ordovícico Inferior y Cámbrico. El área de estudio de este proyecto se
encuentra ubicada al sureste de Cuenca de los Llanos Orientales como se muestra en la figura
1.
Figura 1.Ubicación de la Cuenca de los Llanos Orientales en Colombia. En rojo, la localización del
área de estudio. Modificado de Sarmiento (2011) y ANH (2010)
8
2. OBJETIVOS
2.1. Objetivo General
Realizar un análisis sismo-estratigráfico y estructural, aplicando atributos símicos, con el fin
de determinar geoformas de cuerpos geológicos paleozoicos en una región de los Llanos
Orientales de Colombia.
2.2. Objetivos Específicos
Determinar los atributos sísmicos que caractericen mejor la información
estratigráfica y estructural de la zona de estudio.
Analizar secuencias estratigráficas y estructurales a partir de los datos de pozo y la
información sísmica proporcionada por el SGC.
Analizar las posibles geoformas presentes en una localidad de la cuenca de los llanos
orientales en el Paleozoico con base a la información proporcionada (líneas sísmicas,
registros de pozo y literatura).
Producir mapas de las geoformas paleozoicas presentes en la cuenca de los Llanos
Orientales a partir de la compilación de la información analizada.
3. MARCO GEOLÓGICO
La Cuenca de los Llanos Orientales se distingue por ser una región con bajo relieve al Este
de la Cordillera Oriental. Esta cuenca se caracteriza por ser de tipo antepaís debido a su
geometría y se encuentra limitada por el sistema de fallas de Guaicaramo al oeste y por el
Escudo de Guyana al Este.
3.1. Evolución Tectónica y Estratigráfica
La Formación de la Cuenca de los Llanos Orientales tiene como origen diferentes procesos
tectónicos y sedimentarios a lo largo del tiempo geológico. A inicios del Paleozoico, en el
periodo Cambrico, se presentó un margen pasivo divergente, figura 2, en el cual se forma
una cuenca epicontinental, a causa de las fallas normales que generó este evento las cuales
formaron un semigraben en el cual se depositaron principalmente sucesiones clásticas y
silicoclasticas debido a un aumento del nivel del mar (Ecopetrol & Beicip, 1995, Cediel,
2019).
Posteriormente, a inicios del Ordovícico, se produjo la orogenia Caparonensis debido a una
acreción de retro-arco entre Laurasia y Gondwana. Durante este evento ocurrieron altas
tasas de subsidencia, acreciones de terrenos e inversiones de fallas en el semigraben
previamente rellenado. Estos eventos causaron una extensa deposición de secuencias marinas
transgresivas, conformadas por intercalaciones de lutitas y areniscas figura 2 (Cediel,2019)
(Martínez- Gonzales,2008). Durante el Silúrico continua un régimen extensional, en donde
comienzan a cerrarse los océanos para formar Prepangea (Martinez,2008).
9
A finales del Devónico y e inicios del Carbonífero ocurre una transgresión marina que
sedimenta arcillolitas, areniscas y calizas, figura 2. Posteriormente un evento compresional,
es causado por la orogenia Hercinica, la cual levanta las Cuencas subandinas y las somete a
eventos erosivos. A finales del Carbonífero, esta misma orogenia ensambla Pangea debido a
la colisión de Gondwana con Laurasia (Ecopetrol & Beicip, 1995).
4. DATOS
La información sísmica y geofísica utilizada para este proyecto fue proporcionada por el
Servicio Geológico Colombiano (SGC). La distribución de estos datos se muestra en la figura
3 y consistió en:
Un bloque sísmico 3D de 297km2 con un procesamiento sísmico PSTM (Pre-Stack
Time Migration). El procesamiento sísmico PSTM, se enfoca en realizar la migración
de las trazas sísmicas después de las correcciones estáticas, atenuación de ruido y
deconvolucion para obtener una mejor visualización de la sísmica (ANH,2010).
4 Pozos símicos con registros Gamma Ray, Densidad, Neutrón, Resistividad y Sónico
(únicamente para los pozos 1-2-3).
4 Líneas sísmicas de 35.5 km, 77 km ,60 km y 33 km.
Figura 2.Síntesis estructural de la Cuenca de los llanos orientales en el Paleozoico. Modificado de Martínez, 2008.
10
Figura 3. Distribución de información sísmica 2D (líneas sísmicas), información sísmica 3D (el cubo) y pozos.
Una tabla con la información geofísica proporcionada con más detalle se encuentra en el
anexo A.
5. MARCO TEÓRICO
5.1. Sísmica de Reflexión
La sísmica de reflexión es un método geofísico que proporciona información geológica del
subsuelo; ésta consiste en la detonación de un explosivo en la superficie de un terreno que se
transforma en energía vibratoria a través del medio elástico hacia la profundidad (ANH,
2010). Estas reflexiones son definidas por el coeficiente de reflexión obtenido a partir de la
impedancia acústica (Z); resultante del producto de la densidad (ρ) y la velocidad (v).
5.2. Sísmica 3D
La sísmica 3D consiste en cubo de reflexiones tridimensionales que permiten observar la
distribución de las diferentes estructuras y litologías presentes en el área de estudio, la
dirección de adquisición sísmica y elementos de la sísmica 3D se muestran en la figura 4
(Inlines, Crosslines y Time Slice). Para este estudio se interpretaron las geomorfologías
presentes dentro de un cubo sísmico 3D con procesamiento PSTM.
11
5.3. Estratigrafía de Secuencias
La estratigrafía de secuencias se define como la relación genética de diferentes formaciones
rocosas en un tiempo geológico definido (Posamentier et al., 1998). Así mismo, una secuencia
estratigráfica se define como una unidad limitada por dos discontinuidades o superficies de
erosión. El análisis de estas secuencias permite determinar ciclos estratigráficos de una zona
específica (Vail et al., 1987).
Las secuencias estratigráficas se componen de Tratos de sistema o System Tracks en inglés,
los cuales son sucesiones de sistemas deposicionales representadas en el registro rocoso. La
interpretación de estos registros rocosos indica los cambios eustaticos (a nivel del mar) que
dieron origen al depósito. Se reconocen los siguientes tipos de Tratos de sistema:
TMB, Trato de Mar Bajo en inglés, LST: Lowstand System Track
Consiste en depósitos acumulados durante la caída del nivel del mar.
TMT, Trato Transgresivo, en inglés TST: Transgressive System Track
Demarca depósitos acumulados al iniciar la transgresión marina.
TMA, Trato de Mar Alto, en inglés HST: Highland System Track
Ocurre en eventos de máxima inundación, en donde cambia el régimen de deposición y
disposición de los sedimentos.
Estos cambios eustaticos definen ciclos estratigráficos, consistentes en la relación entre el
espacio de acomodación y suministro de sedimentos (A/S) dentro de dos superficies (Van
Wagoner et al., 1988). Estos ciclos definen superficies de máxima inundación (MSI),
secuencias Regresivas-Transgresivas (TR) y superficies de máxima regresión (SMR).
Figura 4.Elementos de la sísmica 3D. (Tomado de Tellford,1990).
12
5.4. Sismoestratigrafía
La Estratigrafía Sísmica o Sismoestratigrafía, es el estudio de características estratigráficas
y facies sedimentarias en secciones sísmicas de reflexión. A partir de este análisis se pueden
identificar diferentes rasgos geológicos dependiendo de la geometría de los reflectores. En la
figura 5 se muestran las posibles terminaciones de una sección sísmica. En las terminaciones
Downlap, es visible la progradación en la base, las terminaciones Onlap, se presentan con un
ángulo bajo y las truncaciones de los e stratos se presentan a lo largo de una superficie de
disconformidad debido a los efectos erosiónales o estructurales post-deposición.
Figura 5.Terminaciones de los reflectores sísmicos. Tomado del glosario de Schlumlberger,2020.
Según Sangree & Widmier (1977), a partir de la interpretación de los elementos geológicos
mencionados, es posible realizar mapas A-B/C con el fin de reconocer facies sísmicas, figura
6. Para ello se analizan las propiedades de los reflectores teniendo en cuenta criterios como
continuidad, amplitud y frecuencia, mostradas en la figura 7, además del patrón interno de
las reflexiones y la forma dominante de la secuencia como se muestra en la figura 8.
Figura 6.Sistema de códigos para la elaboración de mapas A-B/C y los tipos de configuraciones internas para
definirla (Tomado de Mitchum et al.,1997).
13
5.5. Registros de Pozo
5.5.1. Gamma Ray
Los registros de rayos gamma (GR) son una medición de radioactividad natural de las
formaciones dada en unidades API. Los valores de rayos gamma se basan en la medida de
la emisión de elementos radiactivos como U, Th y K a diferentes profundidades del pozo.
Para este estudio se utilizaron los valores que proporcionan información sobre la litología,
según Stuart et al., (2014), los registros de Gamma Ray muestran valores mayores a 110
API para Arcillolitas y Lodolitas; Limolitas arenosas poseen valores entre 90-110 API y
Arenas corresponden a valores hasta de 90 API.
5.5.2. Resistividad
Los registros de resistividad se caracterizan por medir la propiedad de un material a oponerse
al paso de corriente eléctrica. Este se mide en unidades de ohm, y resulta de gran importancia
debido a que permite diferenciar los componentes de los fluidos en una formación. La
interpretación de estos registros se puede relacionar inversamente con la porosidad, en donde
se asocian altos valores resistivos a bajas porosidades y bajas resistividades con altas
porosidades
5.5.3. Densidad
El registro de densidad (RHOB) en unidades 𝑔
𝑐𝑚3 de formación mide la densidad de electrones
en unidades geológicas por medio del flujo de rayos gamma.
Figura 7. Altas y bajas amplitudes, frecuencias y
continuidades. (Tomado de Sangree &
Widmier,1977).
Figura 8.Patrones de las configuraciones internas de
las reflexiones. (Tomado de Mitchum et al.,1997).
14
5.5.4. Sónico
El registro sónico (DT) es una medición de tiempo en unidades 𝜇𝑠
𝑓𝑡, en el cual una onda
compresional acústica viaja desde la fuente hasta la formación y regresa al receptor, es útil
para identificar litologías y calibrar datos sísmicos.
5.5.5. Densidad Neutrón
El registro de densidad neutrón es una medida de porosidad que se basa en la concentración
de iones de hidrogeno en una formación. Estos neutrones son emitidos por una fuente química
y colisionan con el material de la formación, perdiendo energía. Las unidades de los registros
densidad neutrón se basan en unidades volumétricas(𝑣
𝑣) y se interpretan como un porcentaje
de porosidad.
5.6. Análisis de Electrofacies con registros de Pozo.
Los análisis de electrofacies se realizan con el objetivo de encontrar unidades que responden
de forma similar a diferentes registros eléctricos. A partir de la interpretación de estos
registros se puede obtener información de ambientes deposicionales en un intervalo geológico
de interés. Para este estudio se utilizó la interpretación de los registros proporcionados para
tener una mejor aproximación a la composición del área de estudio.
En la figura 9, se muestran los principales patrones de clasificación de electrofacies asociados
con eventos deposicionales de registro GR. Adicionalmente, en la tabla 1, se muestran los
valores de densidad asociados a diferentes composiciones.
Figura 9.Clasificación de electrofacies por medio de respuestas de registros GR.(Modificado de Rider, 2002 y
Serra,1984).
15
Tabla 1. Valores de registro de densidad obtenidos de Schlumberger, 1981.
5.7. Atributos Sísmicos
Los atributos sísmicos se definen, según Chopra y Marfurt, 2007, como “una medición
cuantitativa de una característica de interés, extraída directamente de los datos sísmicos”
estas características pueden ser geométricas, cinemáticas, dinámicas o estadísticas (Chen and
Sidney, 1997). El objetivo de estos atributos es proveer información detallada al intérprete
de los parámetros estructurales, estratigráficos y litológicos de la sísmica (Taner, 2001).
Adicionalmente, los atributos pueden ser medidos en un instante de tiempo o sobre una
ventana de tiempo (varios tiempos); tanto en una sola traza como en un conjunto de ellas.
De forma general, los atributos derivados del tiempo proporcionan información estructural;
los derivados de la amplitud y de frecuencia, proporcionan información estratigráfica.
5.7.1. Atributo de Coherencia
La coherencia es un atributo geométrico instantáneo, el cual mide la similaridad y semejanza
entre las trazas (Neidell y Tanner, 1971); este atributo se basa en la variación de la traza
sísmica en términos de amplitud, frecuencia, y fase en función de los contrastes de impedancia
acústica y el espesor de las capas. Así pues, este atributo permite extraer información
estratigráfica y estructural como fallas y canales. De esta forma, el atributo presentará
valores de coeficientes cercanos a uno (similitud máxima) al encontrar zonas homogéneas y
al presentarse regiones con discontinuidades estratigráficas, la similitud disminuirá y su valor
se acercará a cero. Este atributo permite visualizar la evolución de rasgos estratigráficos y
estructurales los cuales favorecen la comprensión de depósitos sedimentarios, plegamientos
tectónicos o fallas (Chopra y Marfurt, 2007). La información sobre el cálculo de este atributo
se muestra en el anexo A.
5.7.2. Descomposición Espectral
La descomposición espectral es un atributo general, el cual consiste en la separación de
señales sísmicas en componentes de frecuencia por medio de transformaciones matemáticas
como la Transformada de Fourier. La respuesta de contenido espectral se encuentra asociada
16
a la respuesta sísmica de los reflectores, los cuales dependen de las propiedades y los estratos
de rocas (Neidell y Tanner, 1971). La descomposición espectral realiza un promedio de estas
frecuencias y resulta de gran utilidad para interpretar rasgos deposicionales y elementos
estructurales. (Chopra y Marfurt, 2007).
5.7.3. Atributo de Amplitud RMS
La amplitud RMS (Root Mean Signal) es un atributo general y estadístico de amplitud, el
cual proporciona una estimación de la envolvente de la traza. En general, este atributo
resalta el contenido de energía de una traza sísmica y se utiliza para distinguir canales y
litologías; por ejemplo, altos valores de amplitud RMS se relacionan con composiciones de
altas porosidades como arenas. (Pereira,et al.,2009).
5.7.4. Atributo de Varianza
El atributo de varianza, de forma general hace una correlación de trazas cercanas. Esta se
estima como una medida de dispersión sobre una ventana de tiempo. Este atributo encuentra
el promedio de todas las muestras de ventana, sustrae el promedio de cada muestra y
posteriormente calcula el cuadrado de la diferencia, ecuación (1) (Van Bemmel,P . et al.,
2000); el resultado consiste en el promedio de estos valores. Este atributo se utiliza para
resaltar discontinuidades relacionadas a terminaciones estratigráficas y lineamientos
estructurales (Landmark, 2004).Este atributo difiere del atributo de Coherencia porque
calcula directamente la diferencia del valor de la traza con el valor medio dentro de la
ventana definida.
𝜎2 = ∑ (𝑥𝑖 − �̅�)2 𝑙𝑡=1 (1)
6. METODOLOGÍA
La metodología que se llevó a cabo para alcanzar los objetivos propuestos en el desarrollo
del trabajo de investigación, contempló ocho fases especificadas en la figura 10.
Figura 10.Metodología realizada para el estudio, en amarillo cada fase realizada en este proyecto.
17
Para evaluar la información proporcionada, se utilizó el software de integración geológica y
geofísica Petrel versión 2016, en donde se cargaron los datos sísmicos 2D, 3D y pozos ( de
allí se obtuvo la distribución mostrada en la figura 3).
6.1. Resolución Sísmica Vertical
La resolución vertical se define como el espesor mínimo que puede resolver la sísmica de
reflexión. Esta se obtuvo analizando las frecuencias sobre una ventana de tiempo en el
intervalo de interés. Utilizando la herramienta “Análisis Espectral” de Petrel, se obtuvo un
valor de 45 Hz como frecuencia máxima utilizable a -10 dB (anexo B1). De acuerdo al campo
de velocidad PSTM suministrado, se determinó una velocidad RMS (Root Mean Square)
4400 m/s para la zona de interés entre 1800 ms y 3000 ms.
Finalmente, se usó el criterio de Raleigh, el cual establece que el límite de resolución es de
¼ de la longitud de onda dominante (Sheriff,1997). Límite de resolución sísmica (λ/4),
obteniendo una resolución de 24.4m. Este valor corresponde a una buena resolución sísmica
(para la profundidad del objetivo) muy útil para reconocer elementos estructurales, sin
embargo, dificulta el análisis de elementos estratigráficos ya que no permite reconocer objetos
menores a 24m.
6.2. Correlación de Pozos con la Sísmica 2D y 3D
Luego de cargar el pozo al software con el fin de obtener información geológica lo más
completa posible, se procedió correlacionar los datos sísmicos con los datos de pozo. Esto se
realizó por medio de un sismograma sintético.
6.2.1. Sismograma Sintético
El sismograma sintético busca reproducir la imagen de una traza sintética, a partir de datos
geológicos y parámetros geofísicos tomados en la dirección del pozo (Brown, 1996). Este
sismograma se genera a partir del registro sónico (DT) y registro de densidad (RHOB); con
estos valores se calcula la impedancia acústica y los coeficientes de reflexión (ecuación 1),
para la profundidad total del pozo como se muestra en la figura 11. La convolución estos
coeficientes de reflectividad con la ondicula predeterminada derivan en las trazas de los datos
sísmicos; las uniones de estas trazas conforman un sismograma sintético.
Figura 11.Proceso de obtención de un sismograma sintético (Martínez Cabañas, 1996).
18
Para realizar el sismograma sintético, se utilizó la herramienta de Petrel “Well tie”, la cual
relaciona tiempo con profundidad. Así se logró amarrar las secciones sísmicas (líneas sísmicas
2D 1 y 3, e Inline 650 del cubo 3D), con los registros de los pozos, pozo-1, pozo-2 y pozo-
3. Para todos los sismogramas se utilizó una ondicula sintética de Ormsby de fase cero que
consiste en un filtro trapezoidal que incluye el espectro de las frecuencias que se ajusten
mejor a la sísmica, como se muestra en el anexo B (figura B1). Al sismograma sintético
generado se le debe realizar un shift en tiempo para que las trazas coincidan con las de la
sísmica; esto se realiza moviendo el sismograma levemente en función de tiempo para
asegurar el correcto “amarre” (esto implica un ajuste en el modelo de velocidades).
6.2.2. Correlación de pozos con la sísmica 2D
A partir del amarre de los pozos con la sísmica, se procedió a ubicar los topes de las
formaciones presentes. Esto se realizó subiendo un archivo de texto a la sección “Well tops”,
en donde se encontraban las profundidades de las formaciones proporcionadas por informes
de perforación. Debido a que los pozos perforaron hasta el Ordovícico, la interpretación del
Cámbrico es netamente interpretación.
6.2.3. Interpretación de Horizontes Sísmicos
A continuación, se procedió a interpretar los horizontes correspondientes a los topes de cada
formación. Lo anterior se realizó utilizando la herramienta de demarcación predeterminada
para horizontes de Petrel. Esto consistió en seguir un reflector horizontalmente por las líneas
sísmicas 2D y el volumen 3D. Los horizontes se trazaron tomando como referencia los topes
mostrados en los pozos. Partiendo de estos horizontes se logró definir 6 horizontes que
corresponden a 6 secuencias sísmicas (Secuencia A-B-C-D-E y F).
6.3. Interpretación de Fallas
La interpretación se realizó con el fin de conocer los elementos estructurales presentes en la
zona de estudio con sus respectivas direcciones y tener una mejor visualización del área de
trabajo. Lo anterior se efectuó para todo el bloque 3D y fue posible utilizando de la
herramienta predeterminada para resaltar fallas de Petrel. Adicionalmente para comprobar
la posición de las fallas se realizó un procesamiento de Amplitud ABS (Absoluta) para la
Inline 650. Este procesamiento consiste en el cálculo del valor absoluto de la traza y se
emplea para resaltar los eventos estructurales en el dominio de la Inline.
6.4. Análisis de Electrofacies con registros de Pozo
El análisis de electrofacies se realizó para el pozos-1 y pozo-3, ya que en estos pozos se
tiene información en el Paleozoico, periodo de interés para este estudio. La interpretación
de estos registros se realizó utilizando la figura 8 y la tabla 2.
6.5. Análisis Secuencial
El análisis secuencial se realizó con el propósito de relacionar los registros de pozo como el
Gamma Ray con eventos deposicionales según lo planteado por Sloss el at., (1949). Para
esto, se identificaron los límites de variaciones en los registros GR y Resistivos. Estas
19
variaciones se relacionaron con las litológicas y los límites de formaciones para obtener las
superficies de máxima y mínima inundación y finalmente determinar las secuencias de pozos.
6.6. Sismoestratigrafía
Luego de determinar las secuencias sismoestratigráficas, determinadas en la interpretación
de los horizontes, se procedió a analizarlas en términos de amplitudes y frecuencias. En donde
las amplitudes están dadas por el contraste de densidad y velocidad de las ondas acústicas a
través de las diferentes litologías. La frecuencia es un indicador de la distancia entre
reflectores y depende del espesor entre los estratos (Mitchum et al., 1977). Por último, la
continuidad corresponde a la extensión lateral de los reflectores. Esto con el fin de realizar
los mapas A-B/C. (figura 5,6 y 7).
6.7. Mapas estructurales
Los mapas estructurales se realizaron para el tope de cada secuencia sísmica incluyendo las
fallas para facilitar la comprensión y visualización de las estructuras geológicas en la zona
de estudio. En estos se ubicaron las fallas más relevantes encontradas en cada una de las
Inlines y Crosslines del bloque sísmico 3D.
6.8. Aplicación de atributos sísmicos
Los atributos sísmicos se aplicaron con el objetivo obtener una información más completa de
las geoformas presentes en cada una de las secuencias sísmicas. Se aplicaron los atributos de
Coherencia, Descomposición Espectral, Amplitud RMS y Varianza. Para este estudio se
corrieron varios atributos adicionales como: el atributo de suavizado, frecuencia instantánea,
cuadratura y curvatura; sin embargo, se escogieron los cuatro primeros ya que fueron los que
brindaron mejor información. Cada uno de estos atributos proporcionó información geológica
relevante para reconstruir geomorfologías paleozoicas en la zona de estudio. Esto se fue
posible aplicando la herramienta de Petrel “atributos de superficie” a los topes de las
secuencias que corresponden a los horizontes Ordovícico Medio, Ordovícico Inferior y
Cambrico del bloque sísmico.
En el documento se muestran y analizan los resultados para las secuencias C, D y E (definidas
previamente por los horizontes) ya que pertenecen al intervalo de interés; sin embargo, para
la secuencia D y E, la resolución vertical dificulto la interpretación de geomorfologías. Es
importante mencionar que se aplicaron atributos sísmicos para todas las secuencias con el
objetivo de mostrar en detalle los cambios morfológicos del área de estudio (Secuencia A y
B en el Anexo C). Para la secuencia F, la profundidad es muy alta, la resolución sísmica
muy baja o la posible presencia de un evento de metamorfismo ligero no permitió visualizar
información en esta secuencia.
Para el análisis de los atributos sísmicos se utilizó la tabla 2, la cual resume las propiedades
sísmicas que analiza cada atributo y su interpretación geológica encontrados en la literatura.
20
Atributo Símico
Propiedad sísmica que analiza
Interpretación
Coherencia Similitud entre trazas Fallas y Canales
Descomposición
espectral
Frecuencias
Rasgos deposicionales y elementos
estructurales
Amplitud RMS Energía de la traza (amplitud) Porosidades, Litologías.
Varianza Correlación de trazas cercanas Discontinuidades Estratigráficas y
Estructurales
Tabla 2.Atributos sísmicos utilizados para este estudio con las propiedades sísmicas que analizan y su
interpretación geológica.( Compilado de lo encontrado en Landmark 2004, Pereira, et al.,2009, Chopra y
Marfurt, 2007).
7. RESULTADOS Y ANALISIS
7.1. Sismograma Sintético
El sismograma sintético del pozo-3 con la Inline 560 del bloque sísmico, se muestra en la
figura 12. Allí se pueden observar que los reflectores de la sísmica coinciden con los generados
por en sismograma sintético; esto significa que se realizó un amarre correcto. El mismo
proceso se realizó para el pozo 2 con la línea sísmica 3 y el pozo 3 con el Inline 560 del cubo
sísmico; sus sismogramas sintéticos se presentan en el anexo B.
Figura 12.Sismograma sintético 1, para la Inline 560 del bloque sísmico 3D con el pozo-3. A la izquierda
(a), sismograma generado por Petrel. A la derecha (b), sismograma ajustado a la sísmica.
21
7.2. Correlación de pozos con la sísmica 2D
En la figura 13, se muestran los horizontes de línea sísmica 1 con el amarre y respectivos
topes del pozo-1. El amarre de los pozos proporcionó informacion hasta el horizonte
Palozoico relacionado a el Ordovícico Medio. Para el Ordovícico inferior y Cámbrico, se
utilizo el analisis de estratigrafia de secuencias; es decir, que para cada horizonte se encontro
una disconformidad que se relaciono con un reflector y se interpreto para las lineas sismicas
2D y el volumen 3D
En la figura 14, se muestran los horizontes encontrados en las líneas sísmicas que se
extrapolaron junto con los topes del pozo-3 y pozo-4. Lo anterior amarró de tal forma, que
se sabe con certeza la ubicación de los topes de las formaciones León, Carbonera y Paleozoico
en el bloque sísmico.
Figura 13.Interpretación de la Crossline 659 amarrada a los pozo-3 y pozo-4.
.
Figura 14.Interpretación de horizontes de la línea sísmica 1 a partir del amarre del pozo-1
22
7.3. Interpretación de fallas
La interpretación de fallas se realizó para todo el bloque sísmico. En la figura 15, la Inline
650 del bloque sísmico sin interpretar, en la figura 16, la misma Inline interpretada.
Adicionalmente en la figura 17, se muestra la Inline 650 con un procesamiento ABS (Anexo
B) y en la figura 18 su interpretación. En estas interpretaciones se puede observar, un sistema
de fallas inversas, predominantes en el Paleozoico y algunas fallas que alcanzan formaciones
Cenozoicas. Las fallas presentes tienen buzamientos en dirección NW y SE. Estas se
prolongan aparentemente hasta edades más antiguas al Cambrico (En el anexo C se
presentan las imágenes con mayor tamaño).
Figura 15. Inline 650 sin interpretar.
Figura 16. Inline 650 Interpretada.
Figura 18.Inline 650 con procesamiento ABS simple interpretada
.
Figura 17. Inline 650 con procesamiento ABS
simple, sin interpretar.
23
La presencia de estas fallas inversas a lo largo de todo el bloque sísmico, coinciden con la
información encontrada en la literatura de un antiguo semigraben (fallas normales) formado
por esfuerzos extensionales que fue rellenado y posteriormente se reactivó en fallas inversas
debido a eventos compresionales y que se prolongaron a superficie por ser zonas de debilidad
estructural. En la figura 18, se muestran las fallas presentes en el basamento prolongadas en
superficie a formaciones más recientes como León y Carbonera; a pesar de que esta
interpretación no permite observar la presencia de fallas normales, se muestran fallas inversas
que corresponderían a una inversión que convirtió estas fallas normales en inversas.
7.4. Análisis de Electrofacies con registros de Pozo
De los 4 pozos proporcionados para este estudio mostrados en la figura 19, únicamente los
pozos-1 y pozo-3 poseen registros correspondientes al Paleozoico como se muestra en la figura
20. Por tal motivo, el análisis de resultados para el paleozoico se trabajó con estos dos pozos.
Los pozos-2 y pozo-4, contribuyeron únicamente con información acerca de los topes de las
formaciones superiores.
Figura 19.Registros Gamma Ray para los pozos con
límites de formaciones.
En este análisis de pozos se utilizaron los registros eléctricos de Gamma Ray, Resistividad,
Densidad y Densidad Neutrón. Para determinar las litologías presentes, se usó la metodología
planteada por Stuart et al., (2014) usando las asociaciones de valores Gamma ray y
utilizando los patrones de respuesta de registro Gamma Ray. Se lograron identificar los
siguientes elementos arquitecturales: depósitos sedimentarios en ambientes marino somero y
marino profundo debido a la regresión marina e incremento en la presencia de arenas
encontradas para el Ordovicico Medio.
Figura 20Registros Gamma Ray de los pozos
en el Paleozoico.
24
7.4.1. Pozo-1
Los registros del pozo 1, figura 21, muestran valores de Gamma Ray, para un tiempo de
3634 ms menores a 90 API asociados a arenas. Valores mayores a 110 API correspondientes
a arcillolitas y lodolitas, estos se observan para una profundidad de 3628 a 3600 ft con un
pico en 3633 ft; lo anterior indica un registro con tendencia granodecreciente. La respuesta
del perfil para esta electrofacie es de forma de embudo aserrado. Debido a esto y según lo
reportado en la literatura esta electrofacie se clasifica como un depósito sedimentario de
ambiente marino somero.
Los valores de densidad están en un rango ente 2.4 y 2.7 g/ cm3; asociados con la presencia
de cuarzo y calcita respectivamente; estos cambios en los valores corresponden a su vez, un
cambio de litología. Además, los registros de neutrón muestran una disminución de la
porosidad en profundidad con valores entre 0.16 (16%), asociados a porosidades bajas
correspondientes a lodos y 0.32 (32%), porosidades medianamente bajas que se podrían
asociar preliminarmente a arenas. Por último, las resistividades muestran valores entre 19.5
y 97 Ohm, una baja resistividad asociada con altas porosidades y buenas conductividades
asociadas a arenas y una alta resistividad asociada a baja porosidad y conductividad, que a
su vez se relacionan con composiciones lodosas, según lo planteado por Schlumberger (1981).
7.4.2. Pozo-3
Los registros del pozo-3, figura 22, muestran valores de Gamma Ray mayores a 110 API;
correspondientes a altos contenidos de arcillolitas y lodolitas con tendencia granocreciente.
La respuesta del perfil para esta electrofacie es de forma cilíndrica lisa descrito en la figura
8, en la parte inferior y campana aserrada en la superior. Debido a lo anterior y comprobando
con lo reportado en la literatura, esta electrofacie se podría clasificar como un depósito de
ambiente marino profundo.
Se observa que los registros de densidad corresponden a valores entre 2.6 a 2.9 g/ cm3
asociados a composiciones de Cuarzo, Anhidrita y Dolomita; estas variaciones en los valores
son asociados a cambios de litología. Además, registros neutrónicos se muestran
relativamente constantes para toda la electrofacie con valores de 0.26 (26%), asociados a
Figura 21.Esquema de elementos litológicos y arquitecturales para el pozo-1.
25
porosidades bajas, los cuales que reafirman la presencia de arcillolitas y lodolitas en la sección
analizada. Por último, los registros de resistividad se muestran también constantes a lo largo
de toda la electrofacie con valores de 35.5 Ohm, lo que indica una alta resistividad asociadas
a baja porosidad y baja conductividad, que a su vez se relacionan con composiciones lodosas,
según lo planteado por Schlumberger,1981.
Para los datos bioestratigraficos, se revisó el informe final de perforación del pozo-3, en donde
se encuentran estudios palinológicos utilizados para determinar las edades geológicas. En el
anexo C, se muestran los datos palinológicos obtenidos para todo el pozo-3. Estos indican
que para el Paleozoico a una profundidad: 4380’ a 4410’ probablemente Ordovícico, se
encontró la presencia de acritarcos tales como Baltisphaeridium sp., Multiplicisphaeridium
sp., Veyhachium sp., que sugieren la presencia de sedimentitas del Paleozoicos y condiciones
de sedimentación de un ambiente marino. Este dato proporciona información acerca del
ambiente en el cual vivía el organismo, correspondiente a un ambiente marino.
7.5. Análisis Secuencial.
Para el análisis secuencial se utilizaron los registros GR y resistivo del pozo-3, debido a que
este pozo se encuentra dentro del bloque sísmico. A partir de los datos de registros de Gamma
Ray del pozo-3 (figura 19), se lograron identificar tres superficies de máxima inundación
SMI, obteniendo como resultado la figura 23, la cual se construyó en base a los datos
interpretados. (sin incluir la encontrada en la formación Guayabo, que no se tuvo en cuenta
para este estudio). Sin embargo, debido a que el intervalo de interés se encuentra en el
Paleozoico, se realizó la interpretación correspondiente a ese marcador, presentado en la
figura 24.
Figura 22.Esquema de elementos litológicos y arquitecturales para el pozo-3.
26
Este marcador corresponde al tope del Paleozoico, secuencia de pozo C, con un límite de 443
Ma datado biostratigraficamente como Ordovícico. Bajo esta superficie se identifica un
evento evidenciado en la disminución de relación A/S, limitado hacia la base con una
superficie de máxima inundación y hacia el tope con una superficie de inconformidad.
En esta sección se observan patrones progradantes correspondientes a un trato de nivel alto
(HST) y según el análisis electrofacial, correspondería principalmente a un ambiente marino
somero. Por otro lado, sobre este marcador se identifica un evento correspondiente a un
trato de nivel bajo (LST), en donde se muestran patrones retrogradantes y un aumento de
la relación A/S. El tope se encuentra limitado por una superficie de máxima inundación,
correspondiente a un trato de nivel transgresivo (TST).
Figura 24. Análisis estratigráfico basado en cambios de nivel base en el Paleozoico para el pozo-3.
Figura 23.Análisis estratigráfico basado en cambios de nivel base para el pozo-3.
27
La figura 25, muestra la relación estratigráfica entre el pozo-3 y pozo-1 para la secuencia de
pozo C, en donde se muestran los cambios de tratos de nivel y sus límites con las secuencias.
De forma general, se observa un comportamiento similar en ambos registros GR, que indican
que esta secuencia se encuentra progradando; esto debido a que para ambos registros se
evidencia principalmente una superficie de máxima inundación que cambia a un trato de
nivel transgresivo.
Figura 25.Correlacion estratigráfica basada en los cambios de nivel base para el pozo-3.
7.6. Sismoestratigrafía
Para el área de estudio se identificaron seis horizontes correspondientes a seis secuencias
sísmicas, como se muestra en la figura 26, entre los cuales para el Paleozoico se encontraron
tres secuencias sísmicas correspondientes al Ordovícico medio, Ordovícico inferior y
Cámbrico. De forma general se observan secuencias deposicionales concordantes y paralelas
hacia el Sureste y patrones discordantes sin continuidad en la parte Norte del bloque sísmico.
Para las correlaciones de edades en la sísmica se utilizaron los datos palinológicos reportados
en los informes finales del Pozo-3.
Figura 26.Secuencias establecidas para la zona de estudio Inline 650.
28
A continuación, se muestran las interpretaciones, según la metodología planteada por
Sangree & Widmier (1977) de las secuencias sísmicas ubicadas en el intervalo del Paleozoico;
es decir, las secuencias sísmicas C, D, E y F, las cuales se presentan de la más joven a la
más antigua. Asi mismo, las interpretaciones de las secuencias sísmicas A y B se muestran
en el anexo C.
7.6.1. Secuencia Sísmica C
La secuencia sísmica C, presentada en la tabla 3, corresponde a formaciones paleozoicas cuyo
horizonte base se encuentra en 471.8 Ma y su tope en 443 Ma; es decir, a unidades
Ordovícicas medias a tardías según dataciones bioestratigrafías. Esta secuencia se caracteriza
por conservar espesores a lo largo de la zona de estudio, amplitudes medias a fuertes,
frecuencias altas y terminaciones concordantes al sur y dowlap al norte. Las reflexiones
internas se presentan con patrones subparalelos y continuos a lo largo de la zona Sur; por
otro lado, se observan configuraciones un poco interrumpidas y contorsionadas en el Centro-
Norte de la zona de estudio.
En la sísmica, la sismofacie C-C/P a subparalela contiene reflexiones de amplitudes
moderadamente fuertes, especialmente en la parte media de la sección; esto corresponde a
cuerpos arenosos con intercalaciones de cuerpos y reflexiones de menor amplitud asociadas
a lodolitas, corroborados con los análisis de electrofacies presentados en la figura 21 y 22. La
sismofacie C-D/W, contiene reflexiones medias a fuertes en la parte inferior e
interdigitaciones de cuerpos finos arenosos.
Secuencia Sismofacies Configuración Interna Ejemplo
SECUENCIA C
C-C/
P a
subparalela
Los reflectores son
paralelos, continuos y
presentan amplitudes y
frecuencias
moderadamente altas.
Sismofacie correspondiente al sur de la zona de
estudio. (Crossline- 155).
C-D/ W
Los reflectores son
ondulosos,
interrumpidos con
amplitudes medias y
frecuencias altas.
Sismofacie correspondiente al norte de la zona de
estudio. (Crossline- 1255)
Tabla 3.Resultados de análisis de facies sísmicas de la secuencia C
Tabla 4.
29
7.6.2. Secuencia Sísmica D
La secuencia sísmica D, presentada en la tabla 4, corresponde a formaciones paleozoicas cuyo
horizonte base se encuentra en 488 Ma y su tope en 471 Ma; es decir, a unidades Ordovícicas
tempranas. Esta secuencia se caracteriza por conservar espesores a lo largo de la zona de
estudio, amplitudes y frecuencias moderadamente fuertes y terminaciones concordantes hacia
el centro-sur y terminaciones onlap al norte del área de estudio. Las reflexiones internas se
presentan con dos comportamientos diferentes en la parte superior. Se presentan patrones
subparalelos y continuos a lo largo de la zona Suroeste, además de configuraciones un poco
interrumpidas y contorsionadas hacia el Centro y Noroeste de la zona de estudio.
La sismofacie C-C/P a subparalelos, presenta amplitudes fuertes con reflectores uniformes
en la parte superior de la secuencia y amplitudes medias con reflectores subparalelos y un
poco contorsionados en la parte inferior; esto corresponde a cuerpos de arenas finas con
intercalaciones que presentan reflexiones de menor amplitud, preliminarmente asociados
lodolitas y arcillolitas relacionados a ambientes marinos. La sismofacie Top-C/W se
caracteriza por reflexiones de amplitudes medias y reflectores interrumpidos asociados
corrientes de turbidez relacionadas con depósitos marinos profundos
Secuencia Sismofacies Configuración Interna Ejemplo
SECUENCIA D
C-C/P a
subparalelo
Los reflectores son
continuos, n continuos,
presentan amplitudes
moderadamente fuertes y
altas frecuencias.
Sismofacie correspondiente al suroeste de zona de
estudio.(Crossline- 185).
Top-C/ W
Los reflectores tienen
amplitudes y frecuencias
moderadamente altas, con
ondulaciones, caóticas y
contorsionadas.
Sismofacie correspondiente al noroeste de la zona de
estudio.( Crossline- 1177)
Tabla 4. Resultados de análisis de facies sísmicas de la secuencia C
30
7.6.3. Secuencia Sísmica E
La secuencia sismica E, presentada en la tabla 5, corresponde a formaciones paleozoicas cuyo
horizonte base se encuentra en 542 Ma y su tope en 488.3 Ma; es decir, a unidades Cámbricas.
Esta secuencia presenta reflectores con configuraciones subparalela, ondulosa, interrumpida,
contorsionada y caótica en ciertas zonas. Adicionalmente presenta terminaciones
concordantes con amplitudes moderadamente bajas a lo largo de toda la zona de estudio
Esta sismofacie presenta reflexiones con amplitudes bajas y algunos con fuertes amplitudes
hacia el medio de la sección. Es importante mencionar que la resolución sísmica no permite
interpretar con claridad la configuración de la sismofacie debido a que son unidades muy
profundas o posiblemente se podría tratar de la presencia de una unidad con un
metamorfismo ligero.
13.1.1. Secuencia F
La secuencia F, presentada en la tabla 6, corresponde a formaciones cuyo horizonte tope está
en 542 Ma y su base edades <542 Ma; es decir unidades precámbricas asociadas al
basamento. Esta secuencia presenta reflectores de amplitudes bajas con configuración
caótica, sin continuidad y con poca coherencia. Sus terminaciones en la parte inferior son
concordantes y la superior se asumen son truncaciones.
En la sísmica, esta secuencia presenta algunas reflexiones con amplitudes variadas de tope a
base. No es evidente una continuidad y los reflectores se comportan de forma caótica; debido
a esto, se infiere que esta secuencia corresponde a un basamento cristalino o un basamento
que ha sufrido algún tipo de metamorfismo debido a que la sísmica impide observar elementos
geológicos en detalle.
Secuencia Sismofacies Configuración
Interna
Ejemplo
SECUENCIA E
C-C/W
Los reflectores se
presentan de forma
subparalela,
contorsionada y
caótica con
amplitudes
moderadamente
bajas y frecuencias
medias.
Sismofacie correspondiente al centro de zona de estudio.
Tabla 5.Resultados de análisis de facies sísmicas de la secuencia E.
31
7.7. Mapas Estructurales
La interpretación sísmica del volumen 3D permitió determinar el nivel de las secuencias
sismoestratigráficas C, D y E. Estas secuencias se muestran de forma paralela a lo largo de
todo el bloque símico y definieron un sistema de fallas con dirección predominante SW-NE,
(tendencia de las fallas), evidenciadas en la figura 27.
7.8. Atributos Sísmicos
Los atributos sísmicos se aplicaron con el objetivo de conocer la distribución de los elementos
estratigráficos y estructurales para el área de estudio. De forma general, permitieron
identificar elementos como litologías, canales y fallas (tabla 2).
Estos atributos se calcularon sobre los horizontes sísmicos interpretados, los cuales
corresponden a los topes de las secuencias sísmicas C, D y F y pueden ser visualizados como
mapas. Los resultados se muestran en las figuras 28, figura 29 y figura 30. Las secuencias A
y B se presentan en el anexo C.
Los atributos de Coherencia y Amplitud RMS para este estudio permitieron identificar
principalmente elementos relacionados con litologías y canales, elementos estratigráficos. Por
otro lado, los atributos de Descomposición Espectral y Varianza facilitaron la interpretación
de algunos canales, sinuosidades y fallas únicamente para la secuencia C, pues para
secuencias más profundas no se identificaron claramente elementos estructurales. De esta
forma, se identificó una diferencia de comportamientos entre la zona NW y SE del bloque
sísmico para las secuencias C y D. Debido a esto se realizó una distinción para estas
secuencias en NW y SE. Los resultados de estas interpretaciones se presentan en la tabla 7.
Secuencia Sismofacies Configuración Interna Ejemplo
SECUENCIA F
T-
C/Caótico
Los reflectores se presentan
de forma subparalela,
caótica y sin coherencia
Sismofacie correspondiente al centro de zona
de estudio.
Tabla 6. Resultados de análisis de facies sísmicas de la secuencia F.
32
Figura 27. Mapas estructurales de las secuencias C, D y E. Las fallas interpretadas en el bloque se encuentran
resaltadas en negro, adicionalmente se muestra la ubicación de los pozos (Pozo 3- cruz roja, Pozo 4 – Cruz
Azul).
Mapas Estructurales
33
Para la secuencia C (figura 28) se lograron identificar depósitos de litologías arcillosas o
lodosas debido a al comportamiento de sus propiedades sísmicas (comprobados por los pozos
1 y 3 mostrados en la figura 28, con más detalle en análisis estratigráfico de pozos). Por otro
lado, esta secuencia presenta fallas hacia el Noroeste y predominancia de canales hacia el
Sureste. Estos canales presentan patrones ondulosos evidenciados en los mapas generados
por los atributos y se asocian a canales intermareales, debido a los reportes de literatura que
indican ambientes marinos.
ATRIBUTO SISMICO
SECUENCIA
Coherencia
Descomposición
Espectral
Amplitud
RMS
Varianza
Interpretación
Secuencia C
NW
Alta
similitud
Bajas frecuencias
con presencia de
segmentos con
frecuencias media a
altos
Bajas energías
Porosidades y
amplitudes
Bajas
correlaciones
con segmentos
de correlaciones
medias
DEPOSITOS de
litologías arcillosas y
lodosas
FALLAS
SE
Medias a
bajas para
morfologías
ondulosas
Altas frecuencias en
la morfología
ondulosa
Altas energías
Porosidades y
amplitudes en la
morfología
ondulosa
Bajas
correlaciones
con presencia de
altas
correlaciones
que demarcan
los cuerpos
ondulosos
Sedimentación de
litologías lodosas con
intercalaciones de
arenas
CANALES
Secuencia D
NW
Alta
similitud
Bajas frecuencias
con presencia de
segmentos con
frecuencias medias.
Bajas energías,
porosidades y
amplitudes
Altas
correlaciones
DEPOSITOS de
litologías arcillosas y
lodosas
SE
Valores
medios de
similitud
Frecuencias medias
a bajas para los
patrones ondulosos
Medias a bajas
energías,
porosidades y
amplitudes
Correlaciones
medias a altas
SINUOSIDADES de
litologías arcillosas y
lodosas
Secuencia E
Altas
similitudes
Bajas frecuencias
Bajas energías
Porosidades y
amplitudes
Altas
correlaciones
DEPOSITOS de
litologías arcillosas y
lodosas
Tabla 7.Resultados e interpretación de los atributos sísmicos aplicados
34
SECUENCIA C
Figura 28. Atributos de a) Coherencia, b) Descomposición Espectral, c) Amplitud RMS y d) Varianza para el tope de la secuencia C.
(Horizonte Ordovícico Medio). En la parte inferior el amarre del Pozo1 con la Línea 1 y Pozo-3 con la Inline 650.
35
Para la secuencia D (figura 29), los atributos sísmicos permitieron identificar depósitos de
litologías arcillosas con sinuosidades arcillosas y lodosas. Estructuralmente no fue posible
observar ninguna falla ya que esta secuencia tiene una profundidad alta, en donde la
resolución sísmica impidió la visualización de estos elementos.
Figura 29. Atributos de a) Coherencia, b) Descomposición Espectral, c) Amplitud RMS y d) Varianza
aplicados para el tope de la secuencia D (Horizonte Ordovícico Inferior).
SECUENCIA D
36
Para la secuencia E (figura 30) la interpretación se dificulto debido a que la resolución sísmica
no permitió observar elementos en gran detalle. Sin embargo, se lograron reconocer
comportamientos similares para todo el horizonte asociados con depósitos de litologías
lodosas y arcillosas.
Figura 30.Atributos de a) Coherencia, b) Descomposición Espectral, c) Amplitud RMS y d) Varianza
,aplicados para el tope de la secuencia E (Horizonte Cambrico).
SECUENCIA E
37
8. INTERPRETACIÓN DE RESULTADOS
Con el fin de analizar los mapas de geomorfologías (generados por los atributos)
detalladamente y establecer un ambiente deposicional para la zona de estudio durante el
Paleozoico, se realizó un análisis integral de los resultados obtenidos por secuencias. A
continuación, se presentan los resultados interpretados de la secuencia más joven a la más
antigua para este estudio:
Secuencia C
Esta secuencia se relaciona al periodo Ordovicico Medio debido a datos palinológicos
reportados. Presenta una disminución de relación A/S, limitados hacia la base con una
superficie de máxima inundación que se relaciona a un trato de nivel alto (HTS) y hacia el
tope con una superficie de inconformidad limitada por un trato de nivel bajo a transgresivo
(LST+TST).
Los registros del pozo-3 indican presencia de minerales como cuarzo, calcita y dolomita con
composiciones litológicas dominantemente lodosas. Así mismo, los registros del pozo-1
muestran composiciones lodosas en la base de la sección y presentan interdigitaciones de
cuerpos finos de arenas hacia el tope.
De lo anterior se tiene que unos tratos de nivel alto asociados a superficies de máxima
inundación, depositan composiciones lodosas en la base de la secuencia con presencia de
minerales como dolomita y calcita. Mientras que, para la parte superior de la secuencia se
presenta una disminución en la relación A/S representada en las intercalaciones finas de
cuerpos arenosos. Esto se relaciona con un ambiente deposicional marino somero, figura 31.
Por otro lado, esta secuencia presenta hacia el Sur reflectores paralelos, continuos con
terminaciones concordantes y hacia el Norte los reflectores son ondulosos, interrumpidos con
terminaciones concordantes en el tope y downlaps la base y que definen una discordancia
observada en la tabla 3. Esto se puede relacionar con los mapas estructurales, en donde se
evidencian fallas con dirección SW-NE las cuales predominan en el Noroeste de área de
estudio para esta sección y las cuales se relacionan con eventos compresionales que causaron
la inversión de fallas en la zona de estudio.
Por último, los atributos símicos permitieron observar la ubicación de los depósitos de
sedimentos lodosos y los canales intermareales de intraplaya, que realizaron un leve aporte
de arenas comprobados por el análisis secuencial y de electrofacies. Además de elementos
estructurales predominantes al NW de la sección corroborados con los mapas estructurales.
La formación de estos canales se podría asociar al comportamiento compresivo que afecto la
zona de estudio durante el Paleozoico y que se puede relacionar con la presencia de fallas
inversas para el área.
38
Secuencia D
Para la secuencia D, asociada al Ordovícico Inferior, los reflectores hacia el Sureste de la
zona de estudio son continuos, con terminaciones concordantes. Por otro lado, hacia el
Noroeste de la secuencia se presentan amplitudes y frecuencias moderadamente fuertes, con
ondulaciones, terminaciones Toplap a concordantes. Esto se relaciona con los mapas
estructurales, en donde se muestran fallas con dirección SW-NE predominantes en la zona
centro y sur de la sección, las cuales pudieron dar origen a estas terminaciones estructurales.
Por otro lado, el análisis de sismofacies junto con los atributos sísmicos permitieron dar
información acerca de la litología para este horizonte. Estos análisis indican grandes depósitos
de cuerpos lodosos asociados a depósitos marinos someros o profundos (figura 31). Para esta
secuencia la interpretación de los resultados fue compleja debido a la profundidad de la
secuencia, en donde la resolución sísmica o la presencia de un metamorfismo ligero (debido
a eventos compresivos) impidió la visualización de diferentes elementos estructurales y
estratigráficos.
Secuencia E
En la secuencia E asociada al Cambrico, los reflectores se presentan de forma subparalelos,
contorsionada y caótica con amplitudes moderadamente bajas y frecuencias medias se
asocian a sedimentaciones lodosas y arcillosas. El comportamiento de esta secuencia es
homogéneo para toda el área de estudio, debido a que estructuralmente se presenta un tren
de fallas paralelas con la misma tendencia y dirección SW-NE.
Adicionalmente, a pesar de que la resolución sísmica no permitió una correcta interpretación
de la sección debido a su profundidad o un posible evento metamórfico ligero, los atributos
sugieren que para existen altas similitudes litológicas, lo cual se relaciona con un ambiente
deposicional lodoso muy pasivo asociado a un ambiente marino profundo, figura 31,
corroborando lo interpretado para la estratigrafía de secuencias de la sección.
Es importante tener en cuenta que las secuencias D y E no se han perforado, lo que significa
que no se tienen registros de pozo. Debido a esto, la interpretación geofísica de estas secciones
se realizó a partir de la información de sismofacie, atributos sísmicos y fallas. Por otro lado,
la secuencia F proporciono información acerca del basamento, en donde se identificó un
basamento cristalino o con un metamorfismo debido a distribución de los reflectores.
De la anterior integración de resultados se realizaron los siguientes mapas geomorfológicos,
presentados en la figura 32 en donde se evidencian las litologías, canales intermareales y
fallas.
39
Figura 31. Elementos geomorfológicos de Ambiente a) Marino somero y b) Marino Profundo. Modificado de
Serra,1984.
La integración y análisis de resultados permitieron determinar que, para el área de estudio
durante el Paleozoico, específicamente para el Ordovícico Medio, Ordovícico Inferior y
Cambrico ocurrió un evento progradacional que termina en una discordancia regional, esto
Figura 32. Mapas de Geomorfológicos de las secuencias C (1192ms a 1872ms), D (1808ma a 2196ms) y E (1808ms a 2796ms).
40
debido a la erosión de la secuencia paleozoica superior y mesozoica a causa de las orogenias.
Este evento progradacional se interpreta como una transición de ambientes marino profundo
en el Cambrico a marino somero en el Ordovícico Medio.
9. CONCLUSIONES
El análisis de estratigrafía secuencial utilizando la información de pozos y sísmica 2D y 3D,
facilitó la comprensión de los elementos estratigráficos y estructurales del área de estudio
para el intervalo del Paleozoico. Debido a esto se logró determinar un ambiente deposición
marino profundo para el Cambrico a marino somero en el Ordovícico Inferior, resultados que
aportan información geológica innovadora para esta era en el área de estudio.
El análisis secuencial permitió determinar tres secuencias estratigráficas A, B y C. De estas
se realizó el análisis únicamente para la secuencia C, ya que es la única perteneciente al
Paleozoico. Para esta secuencia se presentan composiciones de material arcilloso-lodoso con
intercalaciones finas de arenas correspondientes a un trato de nivel alto (HST) hacia la base,
relacionado con una superficie de inundación máxima (SMI) y posteriormente a una
disminución en la relación A/S y posterior trato de nivel transgresivo (LST+TST) hacia el
tope.
Se identificaron dos electrofacies, una para el pozo-1 y otra para el pozo-3 correspondientes
a depósitos de ambientes marinos someros y marinos profundos respectivamente. Para los
pozo-1 y pozo-4 no fue posible realizar este análisis ya que los registros alcanzaron
únicamente al tope de la formación paleozoica.
El análisis de sismofacies se realizó para las secuencias sismoestratigráficas A, B, C, D, E y
F. Este análisis se permitió determinar diferencias en los comportamientos en las facies
sísmicas entre la zona Sur y Noroeste de las secuencias C, D y E.
Las fallas inversas con dirección NE-SW, son predominantes en la zona norte del bloque
sísmico para la secuencia C y en todo el bloque para las secuencias D y E, se relacionan con
los eventos reportados en la literatura. Estos indican un rifting a inicios del Cámbrico que
produjo un semigraben de fallas normales que se reactivaron en el Paleozoico debido a un
fallamiento de tipo inverso, estas fallas se prologaron a secuencias más superficiales debido
a que son zonas de debilidad estructural.
Los atributos sísmicos aplicados a la información sísmica 3D (Coherencia, Amplitud RMS,
Descomposición Espectral y Varianza) resultaron ser una herramienta de gran utilidad que
permitió encontrar la distribución de litologías, fallas y asociarlas a eventos. En el Cambrico
predominan lodolitas y arcillolitas relacionadas a ambientes marinos profundos y en el
Ordovícico Inferior a Medio, se identifica una transición en la cual se presentaron canales
41
intermareales de litologías lodosas y arcillosas con presencia de cuerpos arenosos finos hacia
el Suroeste, relacionados a ambientes marinos someros.
A partir de los atributos símicos, se obtuvieron los mapas de geoformas que facilitaron el
entendimiento deposicional de la zona de estudio. Sin embargo, la interpretación de los
atributos no fue tan detallada, para las secuencias sismoestratigráficas D y F, debido a que
la resolución sísmica a estas profundidades no permitió identificar más elementos
estructurales y estratigráficos. De lo anterior se infiere que la calidad de los atributos sísmicos
y el máximo aprovechamiento de los mismos, se encuentra estrechamente relacionada con la
calidad de los datos sísmicos.
Se sugieren registros de pozos que alcancen mayores profundidades, los cuales permitan el
análisis secuencial y de electrofacies y así disminuir la incertidumbre sobre la genética de las
secuencias sismoestratigráficas D y E.
Este estudio proporciona información geofísica innovadora para la exploración en
hidrocarburos; ya que este se encuentra enfocado a la era Paleozoica, en donde se encontraron
geoformas que no habían sido identificadas previamente.
42
10. BIBLOGRAFÍA
Agencia Nacional de Hidrocarburos (ANH) (2010), Universidad Nacional de Colombia. Plan
Cubrimiento Sísmico – Tomo .
ANH, (2010) Manual para la adquisición y procesamiento de sísmica terrestre y su aplicación
en Colombia.
Amorocho, P. R., & Badillo, J. (2012). Influencia de la composición mineral de rocas
arcillosas en la estabilidad de pozos petroleros. Boletín de Geología, 34(1).
Brown, A. (1996). Seismic attributes and their classification. The Leading Edge, 15(10),1090.
Cabañas, W. M. (1996). Interpretación geológica de sismogramas sintéticos. Un ejemplo
aplicado al sondeo Río Segura G-1. Geogaceta, 20(1), 153-156.
Cediel, F. (2019). Phanerozoic Orogens of Northwestern South America: Cordilleran-Type
Orogens. Taphrogenic Tectonics. The Maracaibo Orogenic Float. The Chocó-Panamá
Indenter. In Geology and Tectonics of Northwestern South America (pp. 3-95).
Springer, Cham.
Chen, Q., & Sidney, S. (1997). Seismic attribute technology for reservoir forecasting and
monitoring. The Leading Edge, 16(5), 445-448.
Chopra, S. y K. J. Marfurt, 2007b. SEG, Tulsa, Seismic Attribute for prospect identification
and reservoir characterization. Pp 45-72
Ecopetrol, & Beicip. (1995). Cuenca de los Llanos Orientales, Estudio geológico regional:
santafé de Bogotá (p. 198).
Glosario de yacimientos petrolíferos. Schlumberger Oilfield Glossary (tomado en Abril 2020).
https://www.glossary.oilfield.slb.com/en/Terms/d/downlap.aspx
Landmark (2004). Post Stack Family Reference Manual. Houston: Landmark Graphics
Corporation.
López, J., et al., 2008. Aplicación del método de refracción sísmica para la determinación de
velocidades de ondas P. Tesis Doctoral. Universidad de El Salvador.
Lozano, E., & Zamora, N (2014). Compilación de la cuenca de los Llanos Orientales.
Martínez-González, J. A. (2008). Síntesis Estructural del Paleozoico en las Cuencas
Subandinas.
43
Mitchum, R., Vail, P. & Thompson, S. (1977). The depositional sequence as a basic unit for
stratigraphic analysis: Memoir 26 Seismic Stratigraphy – Applications to
hydrocarbon exploration. America Association of Petroleum Geologists Bulletin, 1,
53-62.
Neidell, N. and Tanner, M. (1971) Semblance and other coherency measures for multichannel
data. Geophysics, 36(3),482–497
Pereira, L. Azevedo, R. Pinheiro, L. Abbassi, H. (2009). Tesis de grado: Seismic Attributes
in Hydrocarbon Reservoirs Characterization. Oporto, Portugal Universidad de Aveiro
Geosciences Pág. 67, 94
Posamentier, H.W., Jervey, M.T. & Vail, P.R. (1988). Eustatic controls on clastic deposition
conceptual framework en: Wilgus, C.K., Hastings, B.S., Kendall, C.G.St.C.,
Posamentier, H.W., Ross, C.A., Van Wagoner, J.C. Sea Level Changes - An
Integrated Approach. Special Publication Society of Economic Paleontologists and
Mineralogists (SEPM), 42, 110–124.
Rider, M. (2002). The geological interpretation of well logs. (2ª ed.). Scotland: Rider– French
Consulting Ltd.
Sangree, J. & Widmier, J. (1977). Seismic interpretation of clastic depositional facies:
Memoir 26 Seismic Stratigraphy – Applications to hydrocarbon exploration. America
Association of Petroleum Geologists Bulletin, 165-184.
Sarmiento, L. F. (2011). Llanos Basins Vol. 9 (p. 184). Medellín, Colombia: ANH-University
EAFIT. Department of Geology.
Schlumberger. (1981). Principios/aplicaciones de la interpretación de registros de pozo
Schlumberger. USA
Serra, O., & Desparmet, J. R. (1984). Análisis de ambientes sedimentarios mediante perfiles
de pozo. Schlumberger
Sheriff, R. E. (1997). Seismic resolution a key element. AAPG Explorer, 18(10), 44-51.
Stuart, J., Mountney, N., Mccaffrey, W., Lang, S. & Collinson, J. (2014). Prediction of
channel connectivity and fluvial style in the flood-basin successions of the Upper
Permian Rangal coal measures (Queensland). America Association of Petroleum
Geologists Bulletin, 98(2), 191–212.
Sloss, L. L. (1949). Integrated Facies Analysis In Sedimentary facies in geologic history:
Conference at meeting of the Geological Society of America held in New York, New
York, November 11, 1948 (Vol. 39, p. 91). Geological Society of America.
44
Suárez, H., Solano Y. 2012. El paleozoico en los llanos Orientales de Colombia: una ventana
en la búsqueda de fuentes de Hidrocarburos. Revista GEO petróleo, 14, 22.
Taner, M. T. 2001. Seismic attributes, CSEG Recorder, Houston, U.S.A., 48-56. US Energy
Information Administration.
Telford, W. M., Telford, W. M., Geldart, L. P., Sheriff, R. E., & Sheriff, R. E.
(1990). Applied geophysics. Cambridge university press.
Vail, P. R., Mitchum Jr, R. M., & Thompson III, S. (1977). Seismic stratigraphy and global
changes of sea level: Part 3. Relative changes of sea level from Coastal Onlap: section
2. Application of seismic reflection Configuration to Stratigrapic Interpretation.
Van Bemmel, P and Pepper R.E.F (2000) Seismic signal processing method and apparatus
for generating a cube of Variance values. U.S. Patent Number6,151,555.
Van Wagoner, J.C., Posamentier, H.W., Mitchum, R.M., Vail, P.R., Sarg, J.F., Loutit, T.S.
& Hardenbol, J. (1988). An overview of sequence stratigraphy and key definitions
en: Wilgus, C.K., Hastings, B.S., Kendall, C.G.St.C., Posamentier, H.W., Ross, C.A.,
Van Wagoner, J.C. Sea Level Changes — An Integrated Approach. Special
Publication Society of Economic Paleontologists and Mineralogists (SEPM), 42, 39–
45
45
11. Anexo A
13.2. Marco Geológico
Precámbrico
A finales del Precámbrico, durante el periodo Ediacariano, ocurre un rompimiento final del
supercontiente Rodinia, formando Laurentia, Siberia, Báltica y Proto-Godwana. Durante
este periodo de rifting, se presentan depósitos intracratonicos y pericratonicos, en donde se
depositaron principalmente areniscas, arcillolitas y calizas (Martínez, 2008).
Mesozoico
En el Mesozoico inferior esfuerzos extensionales, figura A1b, continúan con la apertura de la
cuenca, originando fosas tectónicas en donde se depositaron sedimentos de ambientes
fluviales del Jurásico–Triásico (Lozano-Zamora,2014), como se muestra en la figura A2. En
el Cretácico tardío, ocurre la deposición de sedimentos marinos (figura A1c) debido a un
régimen transgresivo causado por la subducción de la corteza oceánica con el continente,
denominada la orogenia Andina, esta subducción genera el levantamiento de la Cordillera
Central (Cooper et al 1995).
Cenozoico
Los esfuerzos compresionales con el levantamiento de la Cordillera Central dieron origen a
una megacuenca ante-arco, conformada por las cuencas de Valle Medio del Magdalena,
Cordillera Oriental y Llanos; debido al movimiento de la orogenia se inicia la fragmentación
de la megacuenca ante-arco (figura A1d y A1e), ocasionando un decaimiento flexural en la
cuenca de los Llanos Orientales donde se depositaron sedimentos continentales sobre los
marinos. Hacia el final del Paleógeno se levantó́ la Cordillera Oriental (figura A1f) originando
la cuenca ante-arco de los Llanos en donde continuaron depositándose sedimentos
continentales durante el Neógeno que tenían como fuente esta Cordillera (Lozano-Zamora,
2014).
46
Figura A 1.Evolución tectónica de la cuenca de los llanos orientales. Modificado de Sarmiento 2011.
47
Figura A 2. Columna estratigrafica generalizada de la Cuenca de los LLanos orientales. Modificado de Amorocho
y Badillo, 2012 , Lozano-Zamora,2011 y Martínez, 2008.
48
13.3. Datos
En la siguiente tabla se muestran los datos proporcionados por el servicio geológico
utilizados para este proyecto.
13.4. Procesamiento Pre-Stack Time Migration (PSTM)
La PSTM utilizada en este proyecto, es un tipo de procesamiento sísmico cuya principal
característica es la de realizar la migración de las trazas sísmicas después de las correcciones
estáticas, atenuación de ruido y deconvolucion. Estas migraciones se realizan en el dominio
de los planos del offset y generan un volumen de trazas que se apilan posteriormente. Lo
anterior se realiza con el objetivo de ubicar los reflectores en su verdadera localización y
colapsar las difracciones permitiendo resolver mucho mejor las fallas (ANH,2010).
13.5. Clasificación de los Atributos Sísmicos
Taner (2001), divide los atributos en dos categorías generales: geométricos y físicos. Los
atributos geométricos mejoran la visibilidad de la geometría de las características de los datos
sísmicos como rumbo, buzamiento y continuidad. Los atributos físicos se refieren a la física
de los parámetros del subsuelo como la litología, amplitudes, fases y frecuencias.
Datos
Información
Sísmica
3D
Un cubo sísmico 3D de 297 km2 con procesamiento PSTM (Pre-Stack
Time Migration, ver anexo A).
Sísmica
2D
Cuatro líneas con procesamiento PSTM con las siguientes longitudes
Línea-1 35.5 km
Línea-2 77 km
Línea-3 60 km
Línea-4 33 km
Pozos
Cuatro Pozos con los siguientes registros geofísicos.
Pozo Gamma
Ray
Densidad Neutrón Resistividad Sónico
Pozo-1 x X x x x
Pozo-2 x X x x x
Pozo-3 x X x x x
Pozo-4 x X x x
Tabla A 1 Datos sísmicos y geofísicos proporcionados por SGC para el desarrollo del Proyecto.
49
Por ejemplo, los atributos instantáneos son los derivados a partir de una traza sísmica
compleja como la amplitud, fase y frecuencia instantánea. Por otro lado, los atributos
generales se basan en transformadas matemáticas partiendo de los datos sísmicos.
Otra clasificación comúnmente utilizada es de pre-apilado y post-apilado:
Los atributos preapilados: Se calculan antes de apilar las trazas de un registro. Uno
de los más conocidos es la técnica de AVO, la cual consiste en estudiar y caracterizar
las variaciones de amplitudes sísmicas con respecto al offset o ángulo de incidencia
para una familia CDP.
Los atributos postapilado: Son los atributos que se aplican después del apilado de
migración. Estos pueden ser divididos en “atributos de tipo instantáneo” y “atributos
de intervalo” o como “atributos de horizonte” y “atributos de ventana” (Taner et al.,
1994). Los atributos instantáneos, producen un valor del atributo para cada muestra
de la traza y los atributos de intervalo como su nombre lo indica, dan un valor de
atributo para una ventana de trazas seleccionadas (Brown, 1996)
Brown (1996b) clasifica los atributos usando una estructura de árbol en tiempo, amplitud,
frecuencia y atenuación que se dividen a su vez en pre apilado y post apilado (figura A3).
Figura A 3.Clasificación de atributos sísmicos según Brown (1996)
50
13.6. Atributo de Coherencia
En una sección sísmica, los cambios de respuesta corresponden a el comportamiento de la
ondicula sísmica en profundidad debido a su convolución con la geología del suelo. Estos
cambios dependen en gran medida del contraste de la impedancia, los cuales son afectados
por la litología, porosidad y densidad de las capas del subsuelo. Debido a esto, los datos
sísmicos en una sección presentan fuertes variaciones laterales causados por los contrastes de
impedancia como se muestra en la figura A4. El cálculo de la coherencia se realiza mediante
la búsqueda de discontinuidades de la forma de la ondicula.
Figura A 4.Ejemplo de la variación lateral en las formas de las ondas sísmicas: (a) una forma de onda plana,
lateralmente invariante, (b) una forma de onda de geometría sinclinal, lateralmente invariante, (c) una forma de
onda variable lateralmente, indicando cambios laterales de impedancia o espesor, y (d) una forma de onda muy
variante lateralmente asociada a tres canales (Chopra y Marfurt, 2007b).
Su interpretación facilita en gran medida el entendimiento de deposiciones sedimentarias y
eventuales fallamiento y plegamientos. En estos sistemas deposicionales generalmente se
resaltan los canales con geometrías meandricas. Las fallas y canales no son fáciles de
identificar en secciones verticales ya que están relacionadas con la calidad de los datos y el
objetivo de la interpretación(Chopra y Marfurt, 2007).
14. Anexo B: Complementos Metodología
51
14.1. Resolución Sísmica
En la figura B1 se muestra el espectro de frecuencia utilizado para el intervalo de interés
(Paleozoico) en la sísmica proporcionada. Este espectro se utilizó para obtener la resolución
sísmica vertical.
Figura B 1. Espectro de frecuencia obtenida para el intervalo de la Inline 650.
14.2. Sismograma Sintéticos
En la figura B2 se muestra la ondicula de Omsby aplicada al pozo 3 para realizar el
“amarre” con la sísmica.
Figura B 2. Ondicula utilizada para el sismograma sintético del pozo-3
Las figuras B3 y B4 corresponden a los sismogramas sintéticos obtenidos para los pozos 1 y
2 correspondientemente.
52
Figura B 3 .Sismograma sintético para la línea 1 con el pozo-1..
Figura B 4.Sismograma sintético del Pozo-2 con la línea sísmica 3
53
14.3. Correlación de la Información de Pozos con la Sísmica 2D y 3D
En las figuras B5 y se muestra el amarre del pozo 2 con la línea sísmica 3. Por otro lado,
en la figura B6 se muestran los horizontes obtenidos para la línea sísmica 2.
Figura B 5.Amarre de la Línea sísmica 3 con el Pozo 2.
Figura B 6. Interpretación de la Línea sísmica 2.
15. Anexo C: Resultados
15.1. Interpretación de Fallas
54
A continuación, se muestran las figuras ampliadas de la interpretación de las fallas
Figura C 1..Inline 650 interpretada.
Figura C 2nline 650 sin interpretar.
55
Figura C 3.Inline 650 con procesamiento ABS simple, sin interpretar.
Figura C 4.nline 650 con procesamiento ABS simple, interpretada.
56
15.2. Dato Palinológico
Figura C 5.Reporte Palinológico Pozo-3.
15.3. Análisis de Estratigrafía sísmica
15.3.1. Secuencia A
La secuencia A, figura C3, corresponde a la formación León cuyo horizonte base está en el
límite de los 23 Ma y tope a 15.9 Ma. Esta secuencia se caracteriza por mantener espesores
constantes y terminaciones concordantes a lo largo de la zona de estudio. En cuanto al
análisis de las facies sísmicas, las reflexiones internas de amplitudes y frecuencias
moderadamente bajas. Los patrones de reflexión se presentan con una configuración
subparalela e uniforma y en la parte inferior son de forma subparalela.
Las configuraciones de sismofacie C-C/P presentan amplitudes bajas las cuales se asocian a
cuerpos finos de arenas de ambientes marinos someros. Las reflexiones de amplitudes
moderadamente bajas corresponden a ambientes de baja energía. (Sangree, 1978)
Secuencia Sismofacies Configuración Interna Ejemplo
57
15.3.2. Secuencia B
La secuencia B, figura C3, corresponde a la formación Carbonera cuyo horizonte base está
en el límite de los 37.9 Ma y tope a aproximadamente 23 Ma. Esta secuencia se caracteriza
por mantener espesores constantes y terminaciones concordantes. Al Suroriente y centro se
observan reflectores con una configuración subparalela y continua; hacia el Noroeste, se
presentan reflexiones onduladas, contorsionadas e interrumpidas. Las reflexiones internas de
la secuencia presentan amplitudes moderadamente altas.
En la sísmica, la sismofacie C-C/P a subparalela presenta reflexiones de altas amplitudes
asociadas con cuerpos arenosos; las intercalaciones con cuerpos más continuos y bajas
SECUENCIA A
C-C/ P
Los reflectores se
encuentran uniformes y
subparalelos con
pequeñas interrupciones
y bajas amplitudes.
Figura C 6.Resultados de análisis de facies sísmicas de la secuencia A.
Secuencia Sismofacies Configuración Interna Ejemplo
SECUENCIA B
C-C/P a
subparalela
Los reflectores son
subparalelos,
continuos, un poco
ondulados y presentan
fuertes amplitudes.
Sismofacie correspondiente al centro de la
zona de estudio.
C-C/ Ondulada
(W)
Los reflectores tienen
un comportamiento
contorsionado,
ondulosos y poca
continuidad con
amplitudes
moderadamente altas
Sismofacie correspondiente al Noroeste de
la zona de estudio.
Figura C 7.Resultados de análisis de facies sísmicas de la secuencia B.
58
amplitudes correspondientes a llanuras de inundación. La sismofacie C-C/ondulada, contiene
amplitudes moderadamente altas con reflectores interrumpidos corresponden a ambientes
marinos someros
15.3.3. Mapas Estructurales
En la figura C10 se presentan los mapas estructurales de las secuencias A y B
Mapas Estructurales
Figura C 8 Mapas Estructurales de las Secuencias A y B, en rojo y azul los pozos 3 y 4 correspondientemente..
59
15.3.4. Atributos Sísmicos
15.3.4.1. Secuencia A
SECUENCIA A
Figura C11. Resultados de atributos para la secuencia A.
60
15.3.4.2. Secuencia B
SECUENCIA B
Figura C 12. Resultados de atributos para la secuencia B.
61
15.3.5. Atributos Sísmicos a diferentes tiempos (cada 170ms)
SECUENCIA C (-1192 a -1872): Esta secuencia se subdividió en tiempos de 170ms para
visualizar mejor el cambio.
Atributo de
Coherencia
-1192 ms
-1362 ms
-1532 ms
-1702 ms
Ventana de tiempo
t(m/s)
Atributo Aplicado
62
Ventana de tiempo
t(m/s)
Atributo Aplicado
Atributo de Descomposición
Espectral Generalizada
-1192 ms
-1362 ms
-1532 ms
-1702 ms
63
Atributo de
Amplitud RMS
-1192 ms
-1362 ms
-1532 ms
-1702 ms
Ventana de tiempo
t(m/s)
Atributo Aplicado
64
-1192 ms
-1362 ms
-1532 ms
-1702 ms
Ventana de tiempo
t(m/s)
Atributo Aplicado
Atributo de
Varianza
65
SECUENCIA D ( -1808-a - 2196): Esta secuencia se subdividió en tiempos de 150ms para
visualizar mejor el cambio.
Ventana de tiempo
t(m/s)
Atributo Aplicado
-1872 ms
Atributo de
Coherencia
-2022 ms
66
Ventana de tiempo
t(m/s)
Atributo Aplicado
-1872 ms Atributo de Descomposición
Espectral Generalizada
-2022 ms
67
Ventana de tiempo
t(m/s)
Atributo Aplicado
-1872 ms
Atributo de
Amplitud RMS
-2022 ms
68
Ventana de tiempo
t(m/s)
Atributo Aplicado
Atributo de
Varianza
-1872 ms
-2022 ms
69
SECUENCIA E ( -1808 a tiempos mayores de - 2196) : calculado únicamente para -
2040ms, pues a tiempos mayores no la resolución del bloque sísmico no permitió nada.
Ventana de tiempo
t(m/s)
Atributo Aplicado
Atributo de
Coherencia -2042 ms
Atributo de Descomposición
Espectral Generalizada -2042 ms
Atributo de
Amplitud RMS
-2042 ms
70
Atributo de
Varianza
-2042 ms