Cubos de Propiedades Petrofísicas y Petrosísmicas, Utilizando Los Pozos Bauxite-B5, Emerald-A9 y Quartz-A2, Del Campo Cloudspin
Cubos de Propiedades Petrofísicas y Petrosísmicas, Utilizando Los Pozos Bauxite-B5, Emerald-A9 y Quartz-A2, Del Campo Cloudspin
Cubos de Propiedades Petrofísicas y Petrosísmicas, Utilizando Los Pozos Bauxite-B5, Emerald-A9 y Quartz-A2, Del Campo Cloudspin
SEMINARIO DE PETROFÍSICA
TRABAJO FINAL
INGENIERO GEÓLOGO
P R E S E N T A N:
ABSTRACT 14
INTRODUCCIÓN 15
OBJETIVOS 16
UBICACIÓN 17
CAPÍTULO 2: GEOFÍSICA 30
7
2.3.3 SATURACIÓN (ЅF) 42
2.4 GEOESTADÍSTICA 43
2.4.1 VARIOGRAMA 44
3.1 METODOLOGÍA 47
3.2 HORIZONTES 49
3.3 EXTRACCIÓN DE ATRIBUTOS 54
3.3.1 ATRIBUTOS SÍSMICOS 55
3.3.2 STRUCTURAL SMOOTHING 55
3.3.3 VARIANZA 56
3.3.4 CHAOS 56
3.3.5 FRECUENCIA INSTANTÁNEA 57
3.3.6 ATENUACIÓN. 58
3.4 INTERPRETACIÓN DE FALLAS GEOLÓGICAS 59
3.5 MODELADO 3D DE LAS FALLAS. 61
3.6 INTERPRETACIÓN DE ATRIBUTOS SÍSMICOS 62
3.6.1 CHAOS 63
3.6.2 ATENUACIÓN 64
3.6.3 FRECUENCIA INSTANTÁNEA 65
3.7 COMPARACIÓN ATRIBUTOS DE FLUIDO 66
3.8 ZONAS DE PAGA 67
3.9 COMPARACIÓN CUBOS SGEMS VS PETREL®. 68
4.1 METODOLOGÍA 70
4.1.1 CONTROL DE CALIDAD IP® 71
4.1.2 CIMAS Y GRADIENTE GEOTÉRMICO DE LOS POZOS 72
4.1.3 DETERMINACIÓN DEL VOLUMEN DE ARCILLA 73
4.1.4 DETERMINACIÓN DE SATURACIÓN DE AGUA 74
4.1.5 DETERMINACIÓN DE LA POROSIDAD Y PERMEABILIDAD 74
4.1.6 DETERMINACIÓN DE PAYFLAG 75
4.2 RESULTADOS IP® 76
4.2.1 POZO BAUXITE-B5 76
4.2.2 POZO EMERALD-A9 78
4.2.3 POZO QUARTZ-A2 80
8
5.1 METODOLOGÍA 82
5.1.1 IMPORTAR DATOS 82
5.1.2 MALLA CARTESIANA O GRID 83
5.1.3 CREACIÓN DE VARIOGRAMAS 84
5.2 RESULTADOS SGEMS® 86
5.2.1 CUBO POROSIDAD EFECTIVA 86
5.2.2 CUBO SATURACIÓN DE AGUA 87
5.2.3 CUBO VOLUMEN DE ARCILLA 88
5.2.4 CUBO ZONAS DE PAGA 89
CONCLUSIONES 90
BIBLIOGRAFÍA: 91
9
ÍNDICE DE ILUSTRACIONES
Ilustración 1Zona de Estudio Abanico del Mississippi, Google Earth (2019). __________________________ 17
Ilustración 1.1.1.1 A) Posición inicial previo a la ruptura; B) Rotación, extensión de la corteza continental y
expansión del fondo marino; C) Posición actual. (Bird, 2005) _____________________________________ 18
Ilustración 1.1.1.2 Distribución de los principales depocentros en la cuenca al Norte del Golfo de México.
(modificado después de Martin, 1977; y Galloway, 2011). _______________________________________ 20
Ilustración 1.1.1.3 Mapa tectónico del norte del Golfo de México (modificado después de Diegel, 1995).
Tenga en cuenta que la mayor parte sal ha sido evacuada del norte por la ruptura de la plataforma producto
de la carga de sedimentos. ________________________________________________________________ 20
Ilustración 1.1.1.4 Tres ríos importantes que proporcionan sedimentos al área de estudio. _____________ 21
Ilustración 1.1.1.5 Sección (N-S) a través de la cuenca norte del Golfo de México. (mapa modificado y
simplificado después de McBride, 1998) ______________________________________________________ 21
Ilustración 1.1.2.1. Sección sísmica que muestra fallas extensionales desarrolladas en la parte superior de
cuerpos de sal. (Kim, 2012) ________________________________________________________________ 22
Ilustración 1.1.2.2 Sección sísmica que muestra las fallas de rotura de empuje desarrolladas por el cuerpo
salino. (Kim, 2012) _______________________________________________________________________ 22
Ilustración 1.1.2.3 Sección sísmica de tendencia NE-SW aplanada por un límite de secuencia de 0.7 Ma. (Kim,
2012) _________________________________________________________________________________ 23
Ilustración 1.1.2.4 Sección sísmica de tendencia NW-SE que muestra fallas extensionales en la parte superior
de la sal (amarilla) y fallas de empuje que se extinguen en la superficie soldada con sal (roja) (Kim, 2012) 23
Ilustración 1.2.1.1 Abanico del Mississippi, Golfo de México (O´Connell, 1986) _______________________ 24
Tabla 1.2.6.1 Formaciones documentadas para la región Sur de Luisiana (Eguiluz,2001) _______________ 29
Ilustración 2.1.1. Esquema básico de la emisión y recepción y recepción de los rayos reflejados en las distintas
capas reflectoras (CBT, tratamiento sísmico, propiedad de CGG) __________________________________ 31
Ilustración 2.2.1. Clasificación de los registros geofísicos (Martínez M.,2017) ________________________ 34
Ilustración 2.2.1.1. Representación esquemática del desarrollo del potencial espontáneo en un pozo.
Adaptado de Dewan; (Ellis y Singer; 2017 ____________________________________________________ 35
Ilustración 2.2.1.2, Resumen esquemático del comportamiento de la curva SP en una variedad de diferentes
circunstancias de registro comúnmente encontradas. De asquith. (Ellis y Singer,2017). ________________ 35
Ilustración 2.2.4.1. Ejemplo de un registro registrado por la herramienta RAB. El aumento de la resistividad
de bits en A indica que hay una gran cantidad de servidores y. Se detuvo la perforación para colocar la
carcasa. Adaptado de Bonner et al .(Ellis y Singer, 2017) ________________________________________ 37
Ilustración 2.3.1.1 Representación conceptual de la porosidad (Dandekar, 2006) _____________________ 38
Ecuación 1 Porosidad. ____________________________________________________________________ 38
Ecuación 2 Porosidad Total ________________________________________________________________ 39
Ecuación 3 Porosidad Efectiva ______________________________________________________________ 39
______________________________________________________________________________________ 39
Ilustración 2.3.1.1.1 Clasificación de los tipos de porosidad con respecto a su origen (Murray, 1960) _____ 39
Ecuación 4 Permeabilidad _________________________________________________________________ 40
Ecuación 5 Permeabilidad Absoluta _________________________________________________________ 41
Ecuación 6 Permeabilidad Efectiva __________________________________________________________ 41
Ecuación 7 Permeabilidad Relativa __________________________________________________________ 41
Ecuación 8 Saturación ____________________________________________________________________ 42
10
Ecuación 9 Saturación de Aceite ____________________________________________________________ 42
Ecuación 10 Saturación de Gas _____________________________________________________________ 42
Ecuación 11 Saturación de Agua ____________________________________________________________ 42
Ecuación 12 Relaciones Matemáticas de Saturaciones de Fluidos __________________________________ 43
Ecuación 13 Variograma __________________________________________________________________ 44
Ilustración 2.4.1.2.1, Modelos teóricos de variograma, (Miranda M. y Condal A., 2003). _______________ 45
Ilustración 2.4.1.3.., Parámetros del variograma (Ganan J., 2014). ________________________________ 46
Ilustración 3.1.1, Distribución de los pozos en la ventana Petrel ® _________________________________ 47
Ilustración 3.1.2 Metodología ______________________________________________________________ 48
Ilustración 3.2.1, Horizonte Kobe en Petrel ® __________________________________________________ 49
Ilustración 3.2.2 Configuración del Horizonte Caracas en Petrel ® _________________________________ 50
Ilustración 3.2.3 Configuración del Horizonte Dallas en Petrel® ___________________________________ 51
Ilustración 3.2.4, Configuración del Horizonte Houston en Petrel ® _________________________________ 52
Ilustración 3.2.5 Configuración del Horizonte París en Petrel ® ____________________________________ 53
Ilustración 3.2.6, Configuración del Horizonte Salt en Petrel ® ____________________________________ 54
Ilustración 3.3.2.1, Visualización del cubo sísmico con el atributo Structural Smoothing, (Petrel®). _______ 55
Ilustración 3.3.3.1, Visualización del cubo sísmico con el atributo de Varianza, (Petrel®). _______________ 56
Ilustración 3.3.4.1, Cubo sísmico con el atributo de Chaos (Petrel®). _______________________________ 56
Ilustración 3.3.5.1, Cubo sísmico con el atributo de Frecuencia Instantánea (Petrel®). _________________ 57
Ilustración 3.3.6.1, Cubo sísmico con el atributo de Atenuación, (Petrel®) ___________________________ 58
Ilustración 3.3.6.2, Atributo de Atenuación vista en planta, (Petrel®). ______________________________ 58
Ilustración 3.4.1, Fallas macadas (Falla color azul turquesa; Falla 1, Falla color amarillo; Falla 2 y Falla color
azul rey; Falla 3) (Petrel®). _________________________________________________________________ 59
Ilustración 3.4.2, Atributo de Chaos aplicado para el trazado de fallas (Petrel®). _____________________ 60
Ilustración 3.4.3, Atributo de Varianza aplicado para el trazado de fallas (Petrel®). ___________________ 60
Ilustración 3.5.1, Horizontes Caracas, Dallas y Kobe de arriba hacia abajo consecutivamente; principales
afectados por las fallas 1, 2 y 3 (Petrel®). _____________________________________________________ 61
Ilustración 3.5.2, Horizontes: Houston, París y Salt. _____________________________________________ 61
Ilustración 3.6.1 Distribución de los pozos en el cubo sísmico, (Petrel®) _____________________________ 62
Ilustración 3.6.1.1, Time Slice a la altura del Horizonte París, con el atributo sísmico Chaos (Petrel®). _____ 63
Ilustración 3.6.2.1, Time Slice a la altura del Horizonte París, con el atributo sísmico Atenuación (Petrel®). 64
Ilustración 3.6.3.1, Time Slice a la altura del horizonte parís, con el atributo sísmico Frecuencia Instantánea
(Petrel®). _______________________________________________________________________________ 65
Ilustración 3.7.1 Comparacion de los tres atributos (Atenuacion, Chaos y Frecuencia Instantánea) ______ 66
Ilustración 3.8.1 Vista 3D, zonas de Paga importadas de IP® _____________________________________ 67
Ilustración 3.8.2, Vista de costado, zonas de Paga importadas de IP® ______________________________ 67
Ilustración 2, Metodología utilizada en el software IP® __________________________________________ 70
______________________________________________________________________________________ 71
Tabla 4.1.1.1, Control de calidad de los registros geofísicos de los pozos, Bauxite-B5, Emerald-A9 y Quartz-A2
IP® ____________________________________________________________________________________ 71
______________________________________________________________________________________ 72
Tabla 4.1.2.1 Cimas y bases de los pozos Bauxite-B5, Emerald,-A9, Quartz-A2, IP® ____________________ 72
Ilustración 4.1.2.2, cálculo de gradiente geotérmico. ___________________________________________ 72
Ilustración 4.1.3.1, Curvas volumen de arcilla total calculada en el software IP® ______________________ 73
Ilustración 4.1.4.1, Porosity and Water Saturation del software IP®., sirve para calcular saturación de agua 74
Ilustración 4.1.5.1, curvas porosidad y permeabilidad del pozo Bauxite- B5, IP®. _____________________ 74
11
Ilustración 4.1.6.1, parámetros necesarios para la determinación de zonas de paga el software IP® ______ 75
Ilustración 4.1.6.1 Payflag y Resvoir del pozo Bauxite, resultados, IP® ______________________________ 75
Ilustración 4.2.1. , Pozo Bauxite-B5, IP® ______________________________________________________ 76
Ilustración 4.2.2, Pozo Emerald-A9, IP®. _____________________________________________________ 78
Ilustración 4.2.3,Pozo Quartz-A2, IP®. _______________________________________________________ 80
Ilustración 5.1 Metodología utilizada en S-GeMS _______________________________________________ 82
Ilustración 3.1.1 Formato (.txt) para cargar en el Software SGeMs ® _______________________________ 83
Tabla 5.1.2 Valores utilizados en la creación de la malla cartesiana o grid __________________________ 84
Ilustración 5.1.2 Visualización de los pozos y la malla cartesiana, SGeMs®. __________________________ 84
Ilustración 5.2.1, Cubo de porosidad efectiva (PHIE,), SGeMs. _____________________________________ 86
Ilustración 5.2.2, Cubo saturación de agua (Sw), SGeMs. ________________________________________ 87
Ilustración 5.2.3 Cubo volumen de arcilla (Vcl), SGeMs.__________________________________________ 88
Ilustración 5.2.4 Cubo zonas de paga (PayFlag), SGeMs. _________________________________________ 89
12
RESUMEN
Se realizó la evaluación de tres pozos: Bauxite-B5, Emerald A-9 y Quartz A-2; las
formaciones que se interpretaron fueron: Kobe, Caracas, Dallas, Houston, París y Salt.
13
ABSTRACT
The Cloudspin Field is located in the northeast of the Gulf of Mexico and in the Mississippi
fan, it is known as an area with hydrocarbon potential.
The purpose of this study is to carry out the petrophysical and seismic evaluation of this
field with the help of geology, in order to facilitate the interpretation of lithologies and
structures. These characteristics allow us to have an overview about the characteristics of
the area.
The wells analyzed were: Bauxite-B5, Emerald A-9 and Quartz A-2; the formations were:
Kobe, Caracas, Dallas, Houston, Paris and Salt.
As a support material, two softwares specialized in seismic and petrophysics were used.
Petrel® was used for the interpretation of seismic sections, fault modeling, generation of
horizons and application of seismic attributes to facilitate interpretation, on the other hand,
Interactive Petrophysics software (IP®) was used for Well-log data, estimation of
lithologies, facies and payment zones, the latest was imported to Petrel®, where 3D
observation was facilitated.
Finally, integrating the data and results obtained in the interpretation of seismic,
petrophysical and geological data, the evaluation of each formation was carried out, in
order to define those with more potential.
14
INTRODUCCIÓN
Los registros geofísicos de pozo, son técnicas utilizadas para determinar las propiedades y
parámetros físicos de las rocas, para realizar esta evaluación y posteriormente una
interpretación se utilizó el software INTERACTIVE PETROPHYSICS®
La Geología en conjunto con la Geofísica nos permiten determinar las propiedades de las
rocas, es por ello, que en este trabajo se determinaron las características del Yacimiento en
conjunto, es decir sustentando los resultados, de manera que los resultados obtenidos de
cada software e información geológica se pudieran correlacionar entre sí.
15
OBJETIVOS
Objetivo General
Objetivos Específicos
Estimar la litología de cada uno de los pozos con Registros Geofísicos, así como
evaluar las zonas de interés.
16
UBICACIÓN
El Campo Cloudspin o Abanico del Mississippi se localiza en la parte este del Golfo de
México (Ilustración1), el cual es una amplia acumulación arqueada de sedimentos del
Pleistoceno con una cubierta delgada del Holoceno. Esta estructura se presenta en forma de
cuenca desde el Cañón Viejo de Mississippi y causa un leve abultamiento hacia los
contornos batimétricos.
El abanico está delimitado al este por el Escarpe de Florida, al noreste por el Escarpe de
Sigsbee y al sur por el Escarpe de Campeche. Al sureste y al oeste, el abanico termina en la
zona sureste del Estado Florida y en las llanuras abisales de Sigsbee, respectivamente, en
profundidades marinas de aproximadamente 3,200 a 3,400 m.
17
Capítulo 1: Marco Geológico
1.1 Geología Regional
1.1.1 Evolución tectónica del Golfo de México
La etapa del rifting en la cuenca duro alrededor de 70 Ma. y resultó en una rotación en
sentido antihorario del bloque de Yucatán. Se cree que el momento de inicio de esta
rotación fue durante el Calloviano-Oxfordiense, esta rotación y movimiento hacia el sur del
bloque de Yucatán causó una extensión regional e inició la etapa de margen pasivo junto
con la Formación de una corteza oceánica durante el Kimmeridgiano (Salvador et. al,
1991). La ruptura provocó la división de los depósitos de sal de Louann en dos partes: una
al Norte que limita con los Estados Unidos de América y la otra al Sur que limita con
México (Buffler y Thomas, 1994) (Ilustración 1.1.1.1)
18
que limitan a los principales depocentros. En la parte Norte del Golfo la inclinación de las
evaporitas del Jurásico Superior hacia el final del Cretácico Tardío dio lugar al movimiento
descendente a gran escala de la sal por gravedad y la Formación de los primeros diapiros
salinos. El movimiento de la sal durante el Eoceno Medio y Tardío causó la Formación de
una gran capa de sal cubierta por una espesa secuencia de sedimentos terrígenos. Durante el
Mioceno, los depocentros de sedimentos terrígenos se trasladaron al sector de Luisiana del
margen continental y este cambio provocó la Formación de una nueva capa de sal que
cubría una gran área de la plataforma y el talud continental. Una cantidad creciente de
acumulación de sedimentos dio lugar a la creación de diversas estructuras de sal en el límite
Mioceno-Plioceno alrededor de 5 Ma. (Konyukhov, 2008).
19
Ilustración 1.1.1.2 Distribución de los principales depocentros en la cuenca al Norte del Golfo de
México. (modificado después de Martin, 1977; y Galloway, 2011).
Ilustración 1.1.1.3 Mapa tectónico del norte del Golfo de México (modificado después de Diegel,
1995). Tenga en cuenta que la mayor parte sal ha sido evacuada del norte por la ruptura de la
plataforma producto de la carga de sedimentos.
20
La historia de depósito del pleistoceno en el área de estudio está dominada por tres ríos:
Red, Mississippi y Tennessee (Ilustración 1.1.1.4).
Ilustración 1.1.1.4 Tres ríos importantes que proporcionan sedimentos al área de estudio.
El área de estudio está situada en el norte del Golfo de México y se extiende a ambos lados
de la plataforma que separa el talud continental y la Plataforma de Luisiana. El ancho de la
plataforma de Texas Luisiana oscila entre 32 y 90 km. La batimetría de la pendiente está
muy dominada por la actividad salina y la actividad de la sal da lugar a la Formación de
fallas complejas en la pendiente (Ewing, 1991). Según Diegel, (1995) de acuerdo a su
clasificación de provincias tectonoestratigráficas, el área de estudio se encuentra en el
límite del desprendimiento de Plioceno-Pleistoceno en la provincia Tabular Salt-Minibasin.
(Ilustración 1.1.1.5)
Ilustración 1.1.1.5 Sección (N-S) a través de la cuenca norte del Golfo de México. (mapa
modificado y simplificado después de McBride, 1998)
21
1.1.2 Interpretación estructural al noroeste del Golfo de México
El área de estudio se puede caracterizar por patrones estructurales altamente complejos que
consisten en cuerpos de sal alóctonos, fallas extensionales, pliegues y fallas de contracción.
La mayoría de las fallas están relacionadas con la deformación de los cuerpos salinos
alóctonos en el norte del Golfo de México, provocada por la carga de sedimentos (Diegel,
1995; Rowan, 1998). En el área de estudio, las fallas normales comúnmente se
desarrollaron en los límites norte, este y oeste de la mini cuenca (Ilustración 1.1.2.1) y en
general en las direcciones NNE-SSW.
Ilustración 1.1.2.1. Sección sísmica que muestra fallas extensionales desarrolladas en la parte
superior de cuerpos de sal. (Kim, 2012)
Hacia el sur se desarrollaron dos tipos de fallas inversas que se dirigen de este a oeste, las
cuales se interpretan como fallas de rotura y empuje (Ilustración 1.1.2.2) (después de
Rowan et al., 1998), que se hunden tanto hacia el interior como hacia el fondo de la cuenca,
cortando el pliegue salino (Rowan, 1999).
Ilustración 1.1.2.2 Sección sísmica que muestra las fallas de rotura de empuje desarrolladas por
el cuerpo salino. (Kim, 2012)
22
Las fallas que se observan en la parte sureste de la cuenca disminuyen significativamente
hacia el oeste (Ilustración 1.1.2.3). Se cree que estas fallas se activaron en la última etapa
del plegado y continuaron hasta aproximadamente 0.6 Ma. Después de 0.6 Ma, las
actividades de la falla disminuyeron debido a la evacuación de sal desde la parte central de
la cuenca. La otra falla, vista en planta muestra un patrón similar a la falla de rotura. Esta
falla también es una falla inversa, que desciende hacia el interior y se desarrolla
principalmente por traslación de sedimentos depositados en el cuerpo de sal alóctona
(Ilustración 1.1.2.4), (Rowan et al., 1999).
Ilustración 1.1.2.3 Sección sísmica de tendencia NE-SW aplanada por un límite de secuencia de
0.7 Ma. (Kim, 2012)
Ilustración 1.1.2.4 Sección sísmica de tendencia NW-SE que muestra fallas extensionales en la
parte superior de la sal (amarilla) y fallas de empuje que se extinguen en la superficie soldada
con sal (roja) (Kim, 2012)
23
1.2 Geología local
1.2.1 Abanico del Mississippi
El Abanico del Mississippi se encuentra al noreste del Golfo de México y se extiende hacia
el mar de la plataforma continental de Luisiana a profundidades abisales (Ilustración,
1.2.1.1). Su forma de abanico ha sido resultado de las enormes cantidades de turbiditas que
se han depositado desde el Cuaternario y se ha extendido hasta una distancia de más de 600
km del delta del Río Mississippi, alcanzando las Planicies Abisales de Sigsbee y de Florida.
Cubre un área de más de 300,000 km2 y tiene un volumen de aproximadamente
290,000km3. En su parte central, el espesor de sedimentos depositados es de
aproximadamente 4.5km con un tirante de agua de aproximadamente 2500m en la base del
talud continental (Stuart and Caughey, 1976; Moore, et al, 1979).
24
Este abanico es un gran océano abierto de aguas profunda depositado en un margen pasivo.
La mayor parte del Abanico del Mississippi fue depositado durante el Pleistoceno y su
acumulación ocurrió en varias etapas. Al final del periodo de glaciación, en el Wisconsin,
el curso del Río Mississippi, se movió de la región este a la posición que actualmente se
encuentra. Este cambio cortó el sistema del Cañón de Mississippi, como un conducto para
el transporte de sedimento del abanico (Bouma et al, 1989). Su crecimiento no fue
influenciado por las dimensiones de la cuenca. Sin embargo, se depositaron grandes
cantidades de material de origen, como arena y arcilla debido a cambios en el nivel del mar
(Stuart y Caughey, 1976).
La forma de Abanico del Mississippi consiste en varias salientes en forma de lóbulos; cada
uno tiene una forma alargada y son conectados al cañón submarino que atraviesa la
plataforma y el talud continental, llevando grandes flujos masivos de sedimentos que en su
final son depositados en secuencias de turbiditas. El complejo del abanico lobular del cañón
del Mississippi básicamente ha sido formado por los depósitos derramados sobre los bancos
de canales (O’Connell, 1986).
1.2.2. Fisiografía
El Abanico de Mississippi se compone de varios lóbulos del abanico, cada uno tiene una
forma alargada. Durante el Pleistoceno, una migración general de la deposición del lóbulo
del abanico desde el oeste hasta el este y hacia aguas más profundas se produjo un lóbulo
de abanico es básicamente un complejo de canales sobre bancos que se puede dividir en
cuatro áreas (Coleman, 1996). La fisiografía del abanico está caracterizada por el relieve
topográfico de un cañón submarino que atraviesa, las tres provincias fisiográficas del
25
abanico. Basados en datos de perfiles de reflexión sísmica, (Moore 1978), divide el abanico
en tres unidades estratigráficas, abanico superior, medio e inferior.
2. El abanico medio tiene forma de convexo en secciones cruzadas, con canales de forma
senoidal que atraviesan a lo largo del talud; su estructura principal, es un depósito lobular
amplio con una superficie lisa moderadamente y se compone de un complejo masivo de
canales llenos (hasta sus niveles naturales), en donde el material se derrama sobre sus
bancos. El canal central que se encuentra en el talud mide aproximadamente 3 km a lo
ancho, muestra en varios segmentos su fondo del lecho marino. El canal único en el
abanico central es de naturaleza migratoria y de carácter agrario y muestra una similitud
general con los sistemas fluviales serpenteantes. La parte basal parte del relleno del canal
generalmente está cubierto por guijarros barro y arena. El transporte de sedimentos gruesos
debe haber sido principalmente depositado dentro del canal. Los volúmenes de sedimentos
deben haber sido muy grande y la densidad fluye más gruesa que la efectiva.
3. El abanico inferior, se caracteriza por varios canales semiparalelos, es una zona que
tiene varios canales: uno activo y otros abandonados. Su superficie es casi plana, con una
suave inclinación en su pendiente que emerge hasta la planicie abisal hacia el Sureste y
Suroeste. Es caracterizado por tener numerosos canales pequeños de poca profundidad
agrupados en la base del abanico, los cuales han cambiado a través del tiempo. Estas
distribuciones son el resultado de corrientes de turbiditas que se mueven sobre los canales.
Se caracteriza por una depositación de arena en forma de lámina. La litología general de los
sedimentos de la cuenca es de grano muy fino. El área superior del abanico inferior
contiene una canal y un número de canales más o menos paralelos, rastros de depósitos de
tipo dique, y sedimentos con importantes cantidades de arena, derivadas de una fuente con
una baja relación arena / arcilla, se movieron a través del abanico a su área distal, una
distancia de transporte superior a 600 km.
26
1.2.3. Abanicos submarinos
Los abanicos submarinos representan los sedimentos de aguas profundas depositados por
corrientes de turbidez que se originaron en aguas marinas poco profundas. (Normark 1970)
definió tres morfológicas de abanicos de aguas profundas:
27
1.2.5 Geología estructural
1.2.6 Estratigrafía
Wilcox Paleoceno Constituida litológicamente por una secuencia de lutita gris oscuro,
Medio-Superior semidura, arenosa y calcárea, con alternancia de areniscas gris
claro y gris oscuro, de granos finos de cuarzo, regularmente
clasificados, semicompacta en matriz arcillosa y cementante
calcáreo, (Eguiluz, 2001)
Claiborne Eoceno Inferior Litológicamente constituida por una secuencia de lutitas gris claro
y gris oscuro, suave a semidura y ligeramente calcárea. Alternando
con cuerpos de arena, de color gris claro, de granos finos de
cuarzo, subredondeados, regularmente clasificados,
semicompactas, en matriz arcillosa y cementante calcáreo.
(Eguiluz, 2001)
28
Jackson Eoceno Medio- Litológicamente es una secuencia alternada de lutitas y areniscas,
Superior las primeras de color gris claro; las areniscas, son de color gris
claro, de granos finos a medios de cuarzo, subredondeados,
regularmente clasificados, compacta, en matriz arcillosa y
cementante calcáreo, (Eguiluz, 2001)
29
Capítulo 2: Geofísica
2.1 Prospección sísmica
El método sísmico de reflexión se basa en las reflexiones del frente de ondas sísmico sobra
las distintas interfaces del subsuelo. Estas interfaces (reflectores) responden a contrastes de
impedancia que posteriormente se relacionaran con las distintas capas geológicas. Las
reflexiones son detectadas por los receptores (hidrófonos) que se ubican en superficie y que
están alineados con la fuente emisora. Dado que las distancias entre la fuente y los
hidrófonos son pequeñas respecto a la profundidad de penetración que se alcanza
(Ilustración, 2.1.1), el dispositivo experimental soporta que se esté operando en “Corto
Ángulo”, asegurando así la obtención de reflexiones y, distinguiéndose de la sísmica de
refracción o de “Gran Ángulo”.
30
Ilustración 2.1.1. Esquema básico de la emisión y recepción y recepción de los rayos reflejados
en las distintas capas reflectoras (CBT, tratamiento sísmico, propiedad de CGG)
El primer tipo de anisotropía es la Isotropía Transversal Vertical (VTI por sus siglas en
inglés), y la mayor causa para la VTI es la delgada estratificación de las arcillas en el
subsuelo. El segundo tipo es la Isotropía Transversal Horizontal (HTI), y su principal causa
es la presencia de fracturas alineadas verticalmente.( Moreno Hernández, 2014).
31
4.-Efectos relacionados con la interfaz: Cuando una onda encuentra un cambio abrupto en
las propiedades elásticas, alguna parte de la energía se refleja y otra parte se refracta.(Lopéz
J., 2018)
En un diseño ortogonal 3D, es el área que se encuentra delimitada por dos líneas de disparo
adyacentes y dos líneas receptoras adyacentes.
Bin: Se refiere a una pequeña área rectangular que usualmente agrupa dos puntos de
disparo consecutivos y dos puntos receptores consecutivos; se asume que todos los
puntos medios que caen dentro de esta área (bin), corresponden al mismo punto
medio común (CMP), por lo tanto todos los trazos que estén dentro de un mismo
bin serán apiladas en un CMP y contribuirán a la cobertura o fold de ese bin
Dirección Cross – Line: Es la dirección perpendicular a la dirección de las líneas
receptoras.
Dirección In – Line: Es la dirección paralela a las líneas receptoras.
Diseño Ortogonal: Es una herramienta de ingeniería que simplifica los esfuerzos del
diseño estadístico, que permite examinar simultáneamente varias fuentes de
información, gracias al arreglo lineal del diseño, que generan un conjunto especial
de cuadrados.
Líneas Receptoras: Son las líneas o “caminos”, a lo largo del cual se plantan los
receptores o geófonos a intervalos regulares (intervalo de receptores o RI) que es
usualmente el doble de la dimensión del bin en la dirección inline.
Source Interval – (SI): Distancia entre disparos, que es usualmente el doble de la
distancia de la dimensión del bin en la dirección Cross-line.
Líneas de disparo o tiro: Es una línea a lo largo de la cual se colocan los puntos de
tiro (explosivo) a intervalos regulares (intervalo de fuente o SI) es usualmente el
doble de la distancia de la dimensión del bin en la dirección Cross-line
Source Line Interval – (SLI): Distancia entre líneas de disparo sucesivas
Shot Density – (SD): Densidad de puntos de tiro, que se define como el número de
puntos de tiro por Km2.
Templete: Es la suma de un Patch más un Salvo; existen dos tIP®os de
movimientos del templete, que deben ser considerados en un levantamiento sísmico
3D: desplazamiento in-line y desplazamiento Cross-line.
Cobertura o Fold: Corresponde al número de trazos que son apilados dentro de un
punto medio común (CMP). Aunque normalmente se habla de un número promedio
de fold, este varía de bin a bin y para diferentes offsets. La cobertura nominal (full
fold), es la cobertura para el máximo offset y la mayoría de los bin tienen cobertura
nominal.
32
Punto Medio: Es el punto localizado exactamente en la mitad de la distancia entre
una fuente y un receptor.(León D., 2018)
Los Registros Geofísicos de pozos son mediciones de las propiedades físicas de los
yacimientos tanto las rocas y los fluidos, siendo obtenidos a través de métodos geofísicos
indirectos aplicados mediante instrumentos o sondas a lo largo del pozo (Ilustración, 2.2.1).
Los registros se obtienen al hacer pasar los sensores de la sonda enfrente de la Formación
moviendo la herramienta lentamente con el cable.
Los registros son una representación de los datos capturados por la(s) herramienta(s), en
función de la profundidad en algunos casos del tiempo. El registro comienza con las
herramientas posadas en la parte inferior del pozo, para obtener al momento del despegue
que definirá la profundidad total de pozo cuando se recoja el cable, este ascenso se realiza a
una velocidad constante que va a depender de la herramienta. (Martínez M., 2017)
33
Ilustración 2.2.1. Clasificación de los registros geofísicos (Martínez M.,2017)
Medición: registra un fenómeno físico en el pozo, producido por la interacción del agua de
la Formación, el fluido de perforación y las lutitas o arcillas.
Escala: escala logarítmica, de 0.2 a 2000 MM/ m (o línea de 0-20 ohm) medido en mv
(multivolts). La escala de SP no tiene un valor de cero absoluto. Escala típica -20 mV a +80
mV. (Martínez M., 2017)
Aplicaciones:
Indicador de litologías
Indicador de la arcilla de la capa
Identificación de litología
Determina la resistividad del agua de Formación RW
Determinación del diámetro de invasión
34
Ilustración 2.2.1.1. Representación esquemática del desarrollo del potencial espontáneo en un pozo. Adaptado de Dewan;
(Ellis y Singer; 2017)
2.2.2 Sónico
Principio: El sonido emitido desde el transmisor choca con las paredes del agujero,
generando ondas de compresión y cizallamiento dentro de la Formación, ondas de
superficie a lo largo de la pared del agujero y ondas dirigidas en la columna del fluido.
35
Medición: mide el tiempo de retardo que tarda en viajar el pulso de sonido entre transmisor
y el receptor. (Martínez M.,2017)
Aplicaciones:
Principio: los rayos gamma chocan con el detector ocasionando un centello que excita un
fotocátodo, cada fotón de luz libera electrones, estos son acelerados en un campo eléctrico
y producen una corriente proporcional a la energía incidente.
Medición: correlación entre el contenido de arena y la actividad Gamma y mide los niveles
de energía de los Gamma Ray emitidos por la Formación y determina concentraciones de
potasio (k) Uranio (u) y Torio (Th). (Martínez M., 2017)
Aplicación:
36
2.2.4 Registro de resistividad (LLD, MSFL, LLS)
Aplicación:
37
2.3 Propiedades petrofísicas
𝒗𝒐𝒍𝒖𝒎𝒆𝒏 𝒅𝒆 𝒑𝒐𝒓𝒐𝒔
𝛟=
𝒗𝒐𝒍𝒖𝒎𝒆𝒏 𝒕𝒐𝒕𝒂𝒍 𝒅𝒆 𝒍𝒂 𝒓𝒐𝒄𝒂
Ecuación 1 Porosidad.
38
2.3.1.1 Clasificación de la porosidad
La porosidad según la comunicación de sus poros puede ser total o efectiva. La porosidad
total (ϕT) (Ecuación 2): Es la relación de volumen de poros y el volumen total de la roca;
mientras que la porosidad efectiva (ϕe) (Ecuación 3): Es la fracción de volumen poroso que
corresponde a los poros que están conectados entre sí, es afectada por la presencia de arcilla
y esta porosidad es de interés para determinar presencia de hidrocarburos (Pirson, 1965).
Matemáticamente se expresan de la siguiente manera:
Ilustración 2.3.1.1.1 Clasificación de los tipos de porosidad con respecto a su origen (Murray,
1960)
39
2.3.2 Permeabilidad (ĸ)
Se define como la capacidad que tiene una roca de permitir el flujo de fluidos a través de
sus poros interconectados, su unidad de medida es en milidarcy y se simboliza mD.
𝒒∗𝝁∗𝑳
𝚱=
𝑨 ∗ ∆𝑷
Ecuación 4 Permeabilidad
Dónde:
ĸ= Permeabilidad (Darcy).
A= Área (cm2).
40
2.3.2.1 Clasificación de la permeabilidad
Permeabilidad absoluta (𝒌𝒂 ): Es aquella donde hay un solo tipo de fluido fluyendo a
través del espacio poroso de la roca. (Ecuación 5)
𝒒𝒇 ∗ 𝝁𝒇 ∗ 𝑳
𝚱𝒂 =
𝑨 ∗ ∆𝑷
Ecuación 5 Permeabilidad Absoluta
Dónde:
A= Área (cm2).
Permeabilidad efectiva (𝒌𝒆𝒇 ): Es aquella donde el medio poroso de la roca está saturado
por 2 o más fluidos. Considerando que un yacimiento está constituido de tres fluidos (agua,
gas y aceite). (Ecuación 6)
𝒒𝒇 ∗ 𝝁𝒇 ∗ 𝑳
𝚱 𝒆𝒇 =
𝑨 ∗ ∆𝑷
Ecuación 6 Permeabilidad Efectiva
𝚱 𝒆𝒇
𝚱𝒓 =
𝚱 𝒂𝒃𝒔
Ecuación 7 Permeabilidad Relativa
41
2.3.3 Saturación (Ѕf)
La saturación es la relación que existe entre el volumen de los fluidos presentes en una roca
y el volumen de poros total de la roca, (Panesso, 1997). Esta propiedad se expresa
matemáticamente por la siguiente relación (Ecuación 8):
𝑽𝒇
𝑺𝒇 =
𝑽𝒑
Ecuación 8 Saturación
Dónde:
Dónde:
42
Por definición, la suma de las saturaciones de todos los fluidos que se encuentran presentes
en el espacio poroso de una roca, debe ser igual a 1, (Barberii, 1998).
𝑺𝒐 + 𝑺𝒘 = 𝟏
𝑺𝒘 + 𝑺𝒈 = 𝟏
𝑺𝒐 + 𝑺𝒘 + 𝑺𝒈 = 𝟏
Saturación residual (Sxr): Corresponde a la saturación de dicha fase (agua, aceite o gas)
que queda en el yacimiento en la zona barrida, después de un proceso de desplazamiento,
(Ferrer et. al, 1984).
Saturación crítica (Sxc): Corresponde a la mínima saturación requerida para que una fase
pueda moverse en el yacimiento, es decir, corresponde a la máxima saturación a la cual la
permeabilidad relativa de dicha fase es cero, (Ferrer et. al, 1984).
2.4 Geoestadística
La geoestadística es una rama de la estadística que trata fenómenos espaciales (Journel &
Huijbregts, 1978). Su interés primordial es la estimación, predicción y simulación de dichos
fenómenos (Myers, 1987).
La geoestadística está asociada con un conjunto de técnicas, para analizar y predecir los
valores de una variable distribuida en el espacio. Se asumen tales valores implícitamente,
para ser supuestos en correlación entre sí, y el estudio de semejante correlación
normalmente se le conoce como análisis estructural o variograma.
43
2.4.1 Variograma
El variograma se describe como la media aritmética de todas las diferencias entre pares de
valores experimentales separados por una distancia h, (Journel A., 1989), y viene dada por
la siguiente expresión matemática (Ecuación 13):
𝑵𝒑
(𝒁(𝒙𝒊 ) − 𝒁(𝒙𝒊 − 𝒉))𝟐
𝒚(𝒉) = ∑ =
𝟐𝑵𝒑
𝒊=𝟏
Ecuación 13 Variograma
Dónde:
Sill (Meseta): Es la cota superior del variograma. También puede definirse como el límite
del variograma cuando la distancia h tiende a infinito. La meseta puede ser o no finita. Los
variograma que tienen meseta finita cumplen con la hipótesis de estacionariedad fuerte;
mientras que cuando ocurre lo contrario, el variograma define un fenómeno natural que
cumple sólo con la hipótesis intrínseca.
44
2.4.1.2 Modelos matemáticos de variograma
Para interpretar los distintos comportamientos que pueden tener la distribución espacial de
una propiedad no basta con un variograma experimental debido a que por lo general los
resultados de este son imperfectos y presentan gran incertidumbre al momento de
identificar características claves del variograma como la meseta o el rango de influencia, es
por ello por lo que se han desarrollado modelos de variograma teóricos generados a partir
de una ecuación matemática, (Miranda M. y Condal A., 2003).
Los modelos teóricos más usados en geoestadística son los modelos Esférico, Exponencial,
Gaussiano y Lineal. En la Ilustración 2.4.1.2.1 se puede visualizar la gráfica
correspondiente a cada modelo junto con su respectiva función matemática.
45
Ilustración 2.4.1.3.., Parámetros del variograma (Ganan J., 2014).
Dónde:
Lag: El valor de “lag” es la cantidad de veces que se repite la distancia que indica la
separación con la que se realiza el cálculo.
Es una técnica de interpolación basada en la media móvil ponderada, dicho de otra manera,
esta herramienta genera valores en lugares no referenciados a partir de datos conocidos, su
creador Daniel Gerhardus Krige (1960) lo desarrolló para evaluar la distribución de
minerales en minas de Sudáfrica, el principio de esta técnica es hacer que la suma de
errores tienda a cero y que el cuadrado de la desviación sea mínimo.
46
Capítulo 3: Interpretación Sísmica 3D
El área de interés se delimitó por la ubicación de 3 pozos (Ilustración 3.1.1) a partir de los
cuales se obtuvo la información necesaria para su correlación.
3.1 METODOLOGÍA
47
3. Se cargó la información al software, ubicando los pozos en el plano y se importó a
cada uno de los pozos sus desviaciones, es importante cargar correctamente las
desviaciones a los pozos correspondientes ya que errar en este paso podría ocasionar
un mal picado de los horizontes productores.
4. Como tercer paso, se procedió a la extracción de atributos de volumen, esto nos
auxiliará en el siguiente paso que es la interpretación; en el proyecto se realizaron
los siguientes atributos; Atenuación, Chaos, Frecuencia Instantánea y Ant tracking;
para posteriormente realizar la extracción de los valores de los atributos a las
superficies interpretadas.
5. Se realizó la interpretación del cubo sísmico mediante sus propiedades petrofísicas,
se puede realizar al mismo tiempo la interpretación de horizontes y de fallas o bien,
realizar ambas interpretaciones por separado.
6. Una vez que interpretados los horizontes en las secciones sísmicas, se realiza la
configuración del horizonte, esto dará como resultado una superficie “con relieve”
que está acotada en milisegundos y podría tomarse como una representación
“topográfica” del relieve del horizonte interpretado.
7. Posteriormente, se procede a extraer atributos a estas superficies; nos auxiliamos de
nuestros atributos de volumen previamente realizados para generar mapas de
valores de cada superficie, esto nos auxilia en la correlación de la información
petrofísica interpretada en IP®.
8. Una vez finalizados los pasos anteriores el resultado es una interpretación del cubo
sísmico, mediante la aplicación de los atributos en el campo petrolero Cloudspin.
Importar zonas
de paga a Petrel®
y Conlcusiones
48
3.2 HORIZONTES
Este Horizonte está representado por una variación de colores que van de naranja (menor
tiempo), a morado (mayor tiempo).
49
Ilustración 3.2.2 Configuración del Horizonte Caracas en Petrel ®
Este Horizonte está representado por una variación de colores que van de amarillo (menor
tiempo), a azul fuerte (mayor tiempo).
50
Ilustración 3.2.3 Configuración del Horizonte Dallas en Petrel®
Este Horizonte está representado por una variación de colores que van de naranja (menor
tiempo), a azul claro (mayor tiempo).
51
Ilustración 3.2.4, Configuración del Horizonte Houston en Petrel ®
Este Horizonte está representado por una variación de colores que van de naranja (menor
tiempo), a azul (mayor tiempo).
52
Ilustración 3.2.5 Configuración del Horizonte París en Petrel ®
Este Horizonte está representado por una variación de colores que van de rojo (menor
tiempo), a morado (mayor tiempo).
53
Ilustración 3.2.6, Configuración del Horizonte Salt en Petrel ®
Este Horizonte está representado por una variación de colores que van de rojo (menor
tiempo), a morado (mayor tiempo).
Dentro del análisis de los datos sísmicos, se utiliza el término de atributo sísmico para
definir algunos de sus parámetros obtenidos mediante la aplicación de la traza compleja
(Taner y Sheriff; 1979), estas son diferentes opciones de análisis y forma de desplegar las
reflexiones sísmicas con la aplicación de diferentes algoritmos matemáticos.
54
una herramienta auxiliar de interpretación, por lo que se pueden utilizar una vez calibrados
con los registros geofísicos de pozo, para realizar una interpretación cualitativa y en
algunos casos, cuantitativa (Chopra y Marfut; 2007).
El software permite generar volúmenes de atributos a partir del cubo sísmico original, o a
partir de otro cubo de atributos. De igual forma se pueden extraer mapas de atributos con
opción para calcularlos entre dos valores de tiempo, dos superficies o en una misma.
Los atributos sísmicos los utilizamos para una mejor interpretación del área de estudio
principalmente en la estimación de propiedades, posibles aéreas de interés y delimitar
estructuras tales como horizontes, fallas y anomalías referentes a los fluidos relacionados al
campo de estudio.
Ilustración 3.3.2.1, Visualización del cubo sísmico con el atributo Structural Smoothing,
(Petrel®).
55
3.3.3 VARIANZA
En este atributo es la estimación de la varianza local en la señal podemos tener una mejor
definición de los horizontes reflectores e identificarse de las zonas de fallas.
Ilustración 3.3.3.1, Visualización del cubo sísmico con el atributo de Varianza, (Petrel®).
3.3.4 CHAOS
Con este atributo marcamos las fallas de nuestro campo ya que este atributo ilumina las
fallas y discontinuidades relacionadas a características geológicas locales y puede ser
afectado por la migración de gas, intrusión de cuerpos salinos, Ilustración 3.3.4.1.
56
3.3.5 FRECUENCIA INSTANTÁNEA
Este atributo (Ilustración 3.3.5.1) nos ayuda para la visualizar cambios graduales de espesor
y litologías, variaciones como interfaces (hidrocarburo-agua). Las frecuencias más bajas se
observan en reflectores debajo de arenas con gas condesados y yacimientos de aceite, las
bajas frecuencias se asocian solamente en reflectores que se encuentran inmediatamente
debajo de la zona petrolífera.
57
3.3.6 ATENUACIÓN.
58
3.4 INTERPRETACIÓN DE FALLAS GEOLÓGICAS
Para identificar las fallas a primera vista se utilizaron los diferentes atributos, una vez
marcado los horizontes se observaron las discontinuidades y los desplazamientos de estos.
Ilustración 3.4.1, Fallas macadas (Falla color azul turquesa; Falla 1, Falla color amarillo; Falla
2 y Falla color azul rey; Falla 3) (Petrel®).
Teniendo una falla más grande con una tendencia a E-W, esto ocasionado por el intrusivo
salino que se encuentra en la parte SW del bloque.
59
Para marcar fallas se utilizaron los atributos Chaos (Ilustración 3.4.2) y Varianza
(Ilustración 3.4.3); esto para definir de una manera más precisa el comportamiento de las
fallas.
Las zonas que presentan colores brillantes son las que, de acuerdo a las características de
estos atributos, nos pueden indicar zonas de falla.
En estos dos atributos se puede constatar la presencia de tres fallas, tienen tendencia de
fallas normales, las cuales son características de zonas sometidas a regímenes
extensionales.
60
3.5 MODELADO 3D DE LAS FALLAS.
Ilustración 3.5.1, Horizontes Caracas, Dallas y Kobe de arriba hacia abajo consecutivamente;
principales afectados por las fallas 1, 2 y 3 (Petrel®).
Se pueden observar los distintos colores que corresponden a una escala determinada de
tiempo (milisegundos), lo que permite identificar altos y bajos estructurales.
61
Las fallas denominadas 1,2 y 3 inician a una profundidad aproximada de -2000ms en la
escala de tiempo y atraviesan los horizontes Caracas, Dallas y Kobe. Son fallas de tipo
normal, resaltando que la falla número 3 es la más pequeña.
62
3.6.1 CHAOS
Ilustración 3.6.1.1, Time Slice a la altura del Horizonte París, con el atributo sísmico Chaos
(Petrel®).
La más importante dentro de los pozos evaluados en Bauxite B-5, es que se presenta una
anomalía, lo que podría indicar migración de gas.
63
3.6.2 ATENUACIÓN
Ilustración 3.6.2.1, Time Slice a la altura del Horizonte París, con el atributo sísmico Atenuación
(Petrel®).
El Time Slice del atributo de Atenuación (Ilustración 3.6.2.1), nos permite identificar zonas
de fallamiento, ya que se basa en la densidad de fractura, lo cual nos puede indicar un
movimiento del fluido que puede ir migrando a través de fallas hasta llegar a una roca
almacenadora.
Los pozos Bauxite-B5 y Quartz-A2 se encuentran en áreas con de alta atenuación, mientras
que Emerald-A9 muestra muy baja atenuación; lo cual puede indicar la presencia de
movimiento de fluidos en los primeros dos pozos.
64
3.6.3 FRECUENCIA INSTANTÁNEA
Ilustración 3.6.3.1, Time Slice a la altura del horizonte parís, con el atributo sísmico Frecuencia
Instantánea (Petrel®).
Las muy bajas frecuencias nos indican reflectores debajo de arenas con gas condesados y
yacimientos de aceite, las bajas frecuencias se asocian solamente en reflectores que se
encuentran inmediatamente debajo de la zona petrolífera.
65
3.7 COMPARACIÓN ATRIBUTOS DE FLUIDO
Los atributos de fluidos (Ilustración 3.7.1) son un excelente indicador, sin embargo, no
podemos conocer qué tipo de fluido es del que se trata. Utilizamos tres atributos para poder
definir de una manera más certera su presencia.
De acuerdo a las anomalías se puede concluir que en el pozo Bauxite B-5, tenemos más
probabilidad de presencia de fluidos, de tal modo, que puede resultar de interés.
66
3.8 ZONAS DE PAGA
Las zonas de paga fueron importadas del software IP®, y observadas en PETREL® en 3D
(Ilustración 3.8.1), se visualiza Time Slice de la Formación París por ser de interés, con el
atributo Frecuencia Instantánea.
Se puede observar la importancia de los pozos Bauxite B-5 y Quartz A-2, ya que las zonas
de paga se muestran con más abundancia (Ilustración 3.8.2).
67
3.9 COMPARACIÓN CUBOS SGeMs VS PETREL®.
Los atributos utilizados, son un indicador de presencia de fluidos, cada pozo presenta
diferentes peculiaridades, sin embargo, no podemos concluir que tipo de fluido existe ahí
debido a que no se utilizan propiedades petrofísicas en este proceso, de modo que la
sísmica solo es una representación parcial de lo que existe en el subsuelo. Al realizar la
comparación con un software que permite obtener una interpretación mediante modelos
matemáticos y estadísticos, logramos tener un buen panorama de lo que puede haber en la
zona de estudio.
Haciendo una comparación con los atributos Chaos, Atenuación y Frecuencia instantánea,
se detectaron anomalías similares, marcadas con cambios de coloración distintivos de cada
uno como lo muestran las ilustraciones de la Ilustración 3.9.1.
Por otra parte, los Pozos Emerald A-9 y Quartz A-2, dentro de la sección de profundidad
realizada en el software SGeMs, están ubicados en una zona ligeramente azul, ambos pozos
68
se ubican cerca, pero mantienen una desviación evidente en los Time Slice de Petrel®, la
zona en la que se encuentran ubicados, posiblemente sea una zona de paga, pero debido a la
coloración es difícil concluir la presencia de fluido en ese lugar, siendo que en los atributos
del software Petrel® el pozo Emerald A-9 presenta una anomalía únicamente Chaos, y
Quartz A-2 es consistente con una en Atenuación.
Podemos concluir, que en los softwares Petrel® y SGeMs, nos arrojan datos similares,
aplicando diferentes métodos de interpretación, (uno por atributos y el otro por métodos
estadísticos) en el pozo Bauxite-B5, siendo este el que tiene mayores anomalías de fluidos,
detectadas por los tres atributos aplicados (Chaos, Atenuación y Frecuencia instantánea) de
Petrel®, coincidentes con las zonas de paga marcadas en la sección de profundidad de
SGeMs.
69
Capítulo 4: Software IP®
4.1 Metodología
Cimas y
Determinación
Control de gradiente
del volumen de
calidad geotérmico de los
arcilla
pozos
Determinación de Determinación de
Determinación de
la porosidad y saturación de
Payflag
permeabilidad agua
Interpretación y
conclusiones
70
4.1.1 Control de calidad IP®
Para realizar el control de calidad se verifico los registros geofísicos de los pozos Bauxite-
B5, Emerald-A9 y Quartz-A2 en software IP®. Se desplegaron las curvas en los carriles
para verificar los registros con lo que se contaban y la profundidad en las que se tenían
estos registros, al igual que estas fueran de buena calidad. Lo anterior se plasmó en la
siguiente Tabla 4.1.1.1 comparativa. En la cual se observa la ausencia de registros como
rayos gamma, densidad y neutrón, los cuales son de importancia para el cálculo de
propiedades como volumen de minerales entre otras.
Tabla 4.1.1.1, Control de calidad de los registros geofísicos de los pozos, Bauxite-B5, Emerald-
A9 y Quartz-A2 IP®
71
4.1.2 Cimas y gradiente geotérmico de los pozos
Después del control de calidad se realiza una zonificación empleada para determinar las
posibles zonas litológicas presentes en el pozo, estos espesores se determinaron con el
software Petrel® correlacionándose con la geología de la zona. Los espesores se plasmaron
en la tabla 4.1.2.1, en la que podemos observar a la Formación “Caracas” al principio de
todos los pozos evaluados y a la formación “Salt” como la formación más profunda en los
pozos de Bauxite-B5 y Emerald-A9, en el pozo Quartz-A2 se tiene a la Formación Paris
como la más profunda.
Tabla 4.1.2.1 Cimas y bases de los pozos Bauxite-B5, Emerald,-A9, Quartz-A2, IP®
El cálculo del gradiente geotérmico se realizó con 1 grado cada 1000ft, y la temperatura
usada de superficie fue de 35 grados la cual se verifico con el promedio de la temperatura
que se tiene en la zona (Ilustración 4.1.2.2).
72
4.1.3 Determinación del volumen de arcilla
73
4.1.4 Determinación de saturación de agua
Esta propiedad se calculó con el registro sónico (DT), resistividad (ILD), el gradiente
geotérmico y el volumen de arcilla calculado anteriormente. Para posteriormente usarse en
las propiedades de permeabilidad y Payflag (Ilustración 4.1.4.1 ).
Ilustración 4.1.4.1, Porosity and Water Saturation del software IP®., sirve para calcular
saturación de agua
Ilustración 4.1.5.1, curvas porosidad y permeabilidad del pozo Bauxite- B5, IP®.
74
4.1.6 Determinación de Payflag
Payflag son el producto final de la evaluación realizada en los pozos, se le conoce también
como zonas de paga debido a que son zonas donde se puede encontrar petróleo y se puede
extraer con un beneficio monetario. La determinación de Payflag y Reservoroirflag se
realiza con los parámetros de saturación de agua, volumen de arcilla y porosidad total
(Ilustración 4.1.6.1).
75
4.2 Resultados IP®
4.2.1 Pozo Bauxite-B5
76
El pozo Bauxite-B5 se encuentra entre las profundidades 2040.5 y 9439.5ft. En el cual se
van a evaluar las formaciones Caracas, Dallas, Houston, Houston Base, Kobe, Paris y Salt
77
4.2.2 Pozo Emerald-A9
78
El pozo Emerald-A9 se encuentra entre las profundidades 2043 y 8999.5ft. En el cual se
van a evaluar las formaciones Caracas, Dallas, Houston, Houston Base, Kobe, Paris y Salt.
79
4.2.3 Pozo Quartz-A2
80
El pozo Quartz-A2 se encuentra entre las profundidades 2802 y 10366ft. En el cual se van a
evaluar las formaciones Caracas, Dallas, Houston, Houston Base, Kobe, Paris y Salt.
81
Capítulo 5: Software SGeMs®
5.1 Metodología
Realizar la
Cubo de propiedades simulación a partir
del Kriging Ordinario
Antes de importar los datos en el sistema SGeMs® es necesario crear un archivo de entrada
(.txt), el cual consta de un encabezado, que comprende el nombre del campo, el número de
columnas y el nombre de las variables. Específicamente el cuerpo de los datos corresponde
a la coordenada x, coordenada y, profundidad, y los parámetros petrofísicos previamente
calculados en IP®. (Ilustración, 5.1.1)
82
Ilustración 3.1.1 Formato (.txt) para cargar en el Software SGeMs ®
Una vez concluido con el archivo de entrada se procede a cargarlo, mediante la opción
Object en el menú superior y posteriormente Load Object. En la ventana que se despliega
se selecciona Point Set y se establece en que columna se encuentra X, Y, y Z.
Para realizar un modelo 3D se tienen que elaborar una malla cartesiana o grid, para ello se
despliega la pestaña Objects y se selecciona New Cartesian Grid, posteriormente se
selecciona la pestaña Info donde se muestra las coordenadas máximas y mínimas de X, Y, y
Z, estas coordenadas nos servirán para delimitar nuestra área donde se realizara la
simulación de cada propiedad.
En la tabla 5.1.2, se muestran los valores utilizados para crear la malla cartesiana o grid y
en la ilustración 5.1.2. se visualiza el resultado.
83
Distancia entre
Mínima Número de
Coordenadas Máxima (Vmax.) puntos
(Vmin.) cubos
(Vmax-Vmin)
Con los pozos y la malla cargada lo siguiente es crear los variogramas para cada propiedad,
los cuales servirán para llevar a cabo la simulación de nuestras propiedades. Para ello en el
menú Data Analysis se escoge la opción variograma, donde se selecciona la propiedad.
Posteriormente es introducido el número de lags, la separación y tolerancia de los lags, el
azimut, DIP, tolerancia y el ancho de la banda. Una vez generado el variograma se coloca
el efecto Nugget, Sill y también se selecciona el tipo de curva que mejor se ajusta a los
datos (exponencial, esférica y gaussiana).
84
Después de crear el variograma se procede a realizar la simulación. En la pestaña
Simulation se selecciona la opción Sequential Gaussian simulation, donde se carga el
variograma, la malla donde se generará la simulación y se selecciona el tipo de
interpolación (kriging), para este trabajo se empleó Ordinary Kriging.
85
5.2 RESULTADOS SGeMs®
5.2.1 Cubo porosidad efectiva
Teniendo mayor porosidad efectiva en la parte inferior del cubo donde se estaría ubicando a
las formaciones Kobe y París. Por otro lado, tenemos una intercalación de porosidad alta y
baja en la parte media entre las formaciones de Houston y Houston Base. Finalmente, en la
parte superior se tiene una intercalación de porosidad alta y baja, predominando la
porosidad baja en las formaciones de Dallas y Caracas.
86
5.2.2 Cubo Saturación de agua
Con respecto a la Ilustración 5.2.2, es posible observar áreas donde hay mayor
concentración de agua (color rojo), esto nos sugiere que podría estar relacionado con la
presencia de materiales arcillosos. También se puede observar zonas donde la saturación de
agua es baja (color azul claro), consideradas como zonas no invadidas.
Por lo tanto, se puede observar que en la parte superior dentro de las formaciones Caracas y
Dallas, la saturación de agua es elevada. En la parte media dentro de las formaciones
Houston y Houston Base sigue predominando la saturación de agua elevada; mientras que,
en la parte inferior donde se ubica a las formaciones Kobe y París disminuye
considerablemente la saturación de agua respecto a las formaciones anteriores.
87
5.2.3 Cubo Volumen de arcilla
Teniendo una mayor concentración de arcilla en la parte superior entre las formaciones
Caracas y Dallas; posteriormente en la parte media, entre las formaciones Houston y
Houston Base se tiene bajo volumen de arcilla en la mayoría de las zonas; por último, en la
parte inferior entre París y Kobe predominan los valores bajos de arcillosidad.
88
5.2.4 Cubo zonas de paga
En la Ilustración 5.2.4, se puede observar que las zonas de paga son reducidas (color azul),
y la mayoría de ellas se encuentran en la parte sur del cubo, por lo tanto, se puede inferir
que en dichas zonas se cumple con los parámetros de corte. Estos parámetros de corte
consisten en lo siguiente:
89
CONCLUSIONES
En los softwares IP® y SGeMs® podemos correlacionar a la Formación París con respecto
a las propiedades evaluadas como son Porosidad efectiva alta, Saturación de agua baja y
Volumen de arcilla bajo, lo cual quiere decir que posiblemente sea la roca productora
principal de nuestros tres pozos. Por otro lado, tenemos el Horizonte Kobe con propiedades
parecidas al Horizonte París por lo que sería nuestra segunda roca productora.
Al igual en Houston podemos inferir de acuerdo a lo observado en los cubos y pozos que
tiene condiciones para poder ser la roca almacén y en Caracas sería la roca sello
En Petrel®, se pueden constatar las afirmaciones anteriores, los intervalos más importantes
de cada pozo se ubican en el Horizonte París, al realizar las correlaciones con los softwares
anteriores, se determinó que aquí se encuentra el Horizonte con mayor importancia.
90
BIBLIOGRAFÍA:
91
Ellis, D. Y Singer, J. (2008). Well Logging For Earth Scientists.
Http://Old.Geology.Lnu.Edu.Ua/Phis_Geo/Fourman/Library-
Earth/Well%20Logging%20for%20Earth%20Scientists.Pdf, Consultado el 14 de
Abril de (2019).
Ewing, T. E. (1991). Structural Framework. The Gulf Of México Basin: Geological
Society of America, The Geology Of North America, V. J, 31-52.
Ferrer, M. P. (1984). Fundamentos De Ingeniería De Yacimientos.
Galloway, W. E., Whiteaker, T. L., and Ganey-Curry, P. (2011). History Of
Cenozoic North American Drainage Basin Evolution, Sediment Yield, And
Accumulation In The Gulf Of México Basin. Geosphere, 7(4), 938-973.
Ganan J., (2014). “Caracterización de un reservorio aplicando técnicas
geoestadísticas para un campo del distrito Amazónico”. Tesis de Grado. Escuela
Politécnica Nacional. Quito-Ecuador. Págs. 224.
Google Earth, Abanico Del Mississippi Marzo (2019).
Harding, T. P., and J. D. Lowell, (1979), Structural styles, their plate-tectonic
habitats and hydrocarbon traps in petroleum provinces: AAPG Bulletin, v. 63, p.
1016-1058.
Journel A., (1989). “Fundamentals Of Geoestatistics In Five Lessons”. American
Geophysical Union. Volumen 8. 1-57.
Journel, A. G. And Huijbregts. (1978). Mining Geostatistics. London: Academic
Press.
Kim, B. (2012). Seismic Sequence Stratigraphy Of Pliocene-Pleistocene Turbidite
Systems, Ship® Shoal South Addition, Northwestern Gulf Of México. Master
Degree. Texas A&M University.
Konyukhov, A. I. (2008). Geological Structure, Evolution Stages, And Petroliferous
Complexes Of The Gulf Of México Basin. Lithology And Mineral Resources,
43(4), 380-393.
Lehner, P. (1969). Salt tectonics and Pleistocene stratigraphy on continental slope of
northern Gulf of Mexico. American Association of Petroleum Geologists Bulletin
53, 2431–79.
López, J., Ramírez, O. Y León, D. (2018). Análisis Petrosísmico 3-D De Las
Formaciones Caracas, Houston, Kobe, Sal, Bloque 3 Del Campo Cloudspin (Tesis
De Licenciatura). Instituto Politécnico Nacional, CDMX, México.
Martin, R. G. (1984). Northern and Eastern Gulf of Mexico continental margin:
stratigraphic and structural framework. In Framework, Facies, and Oil-trapping
Characteristics of the Upper Continental, pp. 21–42. American Association of
Petroleum Geologists, Studies in Geology no. 7
Martin, R. G. Jr. (1978) Northern And Eastern Gulf Of Mexico Continental Margin:
Stratigraphic And Structural Framework.In: Framework, Facies, And Oil-Trapping
92
Characteristics Of The Upper Continental Margin (Eds A. H. Bouma, G. T. Moore,
And J. M. Coleman), Am. Assoc, Petrol Geol Stud. Geol No. 7, Pp. 21-42
Martin, R. G., And Bouma, A. H., (1977). Physiography Of Gulf Of Mexico. In
Bouma, A. H., Moore, G. T, And Coleman, J. M. (Eds.), Framework, Fades, And
Oil-Trapping Characteristics Of The Upper Continental Margin. Am. Assoc. Pet.
Geol. Stud. Geol., 7: 3-19.
Martínez, M. (2017). Estudio Petrosismico 3-D De La Formación JSK En Un
Campo Del Suereste De México (Tesis De Licenciatura). Instituto Politécnico
Nacional, CDMX, México.
Mcbride, B. C. (1998). The Evolution Of Allochthonous Salt Along A Mega-
Regional Profile Across The Northern Gulf Of México Basin. AAPG Bulletin,
82(5B), 1037-1054.
Miranda M. Y Condal A., (2003). Importancia Del Análisis Estadístico Exploratorio
En El Proceso De Interpolación Espacial: Caso De Estudio Reserva Forestal
Valdivia. Bosque, Vol. 24 N° 2, Pp. 29-42.
Moore, G. T, Starke, G. W., Bonham, L. C , And Woodbury, H. O., (1978).
Mississippi Fan-Gulf of Mexico—Physiography, Stratigraphy, And Sedimentational
Patterns. In Bouma, A. H., Moore, G. T., and Coleman, J. M. (Eds.), Framework,
Fades, And Oil-Trapping Characteristics Of The Upper Continental Margin. Am.
Assoc. Pet. Geol. Stud. Geol., 7:155-191.
Moore, G. T., Woodbury, H. O., Worzel, J. L., Watkins, J. S., And Starke, G. W.,
(1979). Investigations Of The Mississippi Fan. Gulf Of Mexico, Am. Assoc. Pet.
Geol. Mem., 29:393-402.
Morales, S. (2018). Análisis Petrosísmico 3-D De Las Formaciones Dallas, Houston
Base Y París, Bloque 8 Del Campo Cloudspin (Tesis de Licenciatura). Instituto
Politécnico Nacional, CDMX, México.
Moreno, R. Y Moreno, J. (2014). Secuencia De Procesamiento Y Migración "Pre-
STM" De Datos Sísmicos Marinos (Tesis de Licenciatura). Instituto Politécnico
Nacional, CDMX, México.
Murray, R. C. (1960). Origin of Porosity in Carbonate Rocks. J. Sediment. Petrol
30, 59-84.
Myers, D. E. (1987). Optimization of Sampling Locations For Variogram
Calculations. Water Resources Research. 23(3): 283(93).
Normark, W. R., 1970. Growth Pattern Of Deep-Sea Fans. Am. Assoc. Pet. Geol.
Bull., 54:2170-2195.
O'Connell, S., (1986). Anatomy of Modern Submarine Depositional And
Distributary Systems. Ph.D. Dissertation, Columbia University, New York, N.Y.,
302 Pp. (Unpublished)
Panesso, R. D. (1997). Saturación De Fluidos: Conceptos Y Métodos De Cálculo
Utilizando Los Registros Eléctricos. Ciudad De Bolívar.
93
Pirson, S. J. (1965). Ingeniería de Yacimientos Petrolíferos. Barcelona: Ediciones
Omega.
Randen, PCT Patent Application Nopct/IB99/01040 AUTOMATED
STRATIGRAPHIC AND FAULT INTERPRETATION (2002).
Rowan, G. M., And P. Weimer, (1999), Salt-Sediment Interaction, Northern Green
Canyon and Ewing Bank (Offshore Louisiana), Northern Gulf Of México: AAPG
Bulletin, V. 82, No. 5B, P. 1055-1082.
Rowan, M. G., M. P. A. Jackson, And B. D. Trudgill, (1998), Salt-Related Fault
Families and Fault Welds in the Northern Gulf Of México: AAPG Bulletin, V. 83,
No. 9, P. 1454-1484.
Salvador, A. (1987). Late Triassic-Jurassic Paleogeography And Origin Of Gulf Of
México Basin. Am. Assoc. Pet. Geol., Bull.; (United States), 71(4).
Salvador, A. (1991). Origin And Development Of The Gulf Of México Basin. The
Gulf Of México Basin, 389-444.
Schlumberger Wireline & Testing, “Log Interpretation Principles/Applications”,
Seventh Printing, Schlumberger (2009), Chapter 3, Pg. 7
Stuart, C. J., And Caughey, C. A., (1976). Form and Composition Of The
Mississippi Fan. Trans. Gulf. Coast Assoc. Geol. Soc, 26:333-343.
Taner, M.T, F. Koehler And R.E Sheriff (1979), Complex Seismic Trace Analysis
Geophysics, V.
94