Datos Junin - Pregrado
Datos Junin - Pregrado
Datos Junin - Pregrado
=
=
n
i
R i R
y y SS
1
2
) (
=
=
n
i
AB i T
y y SS
1
2
) (
33
Donde:
i: corresponde a cada estacin;
n = nmero total de estaciones.
i
y
(SS
R
) = Valor de la Anomala Regional de grado k.
R
y
= Media aritmtica de la Anomala Regional de grado k.
i
y
(SS
T
) = Valor de la Anomala de Bouguer.
AB
y
= Media aritmtica de la Anomala de Bouguer.
Las anomalas regionales representan la configuracin estructural de la corteza
terrestre, caracterizadas por una longitud de onda larga. En cambio, las anomalas
residuales se deben a efectos de carcter local, es decir a eventos del subsuelo
relativamente prximos a superficie. Los mapas regionales y residuales entonces
proporcionarn informacin que permite interpretar la configuracin estructural del
basamento, as como de diversas estructuras ms someras.
A partir de la interpretacin cualitativa del mapa de anomala de Bouguer se
eligen los posibles perfiles a modelar, contemplando que stos cumplan con los
objetivos del estudio. De dichos perfiles se obtendr la informacin de la gravedad
observada, que permitir realizar el modelado gravimtrico, ajustndose a valores de
densidad de los cuerpos de roca propios de la zona (obtenidos por ncleos, registros
de pozo, etc.) y a la profundidad y geometra de los mismos (mediante modelos
geolgicos previos, ssmica 2D y 3D, etc.).
3.2 Mtodo Magntico
El mtodo magntico, quiz la ms antigua de las tcnicas de exploracin
geofsica, inici su auge despus de la llegada de los levantamientos aeromagnticos
en la Segunda Guerra Mundial. Con las mejoras en la instrumentacin, filtrado de
34
datos y visualizacin, hoy da es posible cartografiar superficies magnetizadas
completas a diferentes escalas, desde el basamento ms intensamente magntico a
escala regional, hasta contactos sedimentarios dbilmente magnetizados a escala local
(Nabighian et al., 2005).
Antiguas civilizaciones, entre ellos griegos, chinos, rabes y escandinavos, ya
usaban los principios que rigen los mtodos magnticos, desarrollndose los primeros
estudios del magnetismo terrestre en el siglo XVI. En la actualidad el mtodo
magntico es la herramienta de exploracin primaria en la bsqueda de minerales
(hierro, molibdeno, diamantes y otras piedras preciosas) (Nabighian et al., 2005).
En cuanto a la exploracin petrolera, los mtodos magnticos han sido
utilizados desde principios del siglo pasado, pero slo para investigar zonas
sumamente falladas y para cartografiar el basamento gneo-metamrfico. Sin
embargo, recientes avances indican que, bajo condiciones favorables (y
especialmente en combinacin con otros mtodos geofsicos y geoqumicos), las
tcnicas magnticas pueden jugar un rol ms importante en la localizacin de
hidrocarburos (Eventov, 1997).
Nabighian et al. (2005) exponen que los nuevos avances tecnolgicos, antes
descritos, han impulsado el uso del mtodo magntico en ramas como la arqueologa,
geotecnia, entre otras, as como la aplicacin en un rango ms amplio en la
prospeccin petrolera, para fines como la localizacin de fallas intrasedimentarias,
definicin de contactos litolgicos poco perceptibles antes y cartografiar domos
salinos en sedimentos dbilmente magnetizados.
Por su lado, Eventov (1997) indica que el uso del mtodo magntico para
detectar directa o indirectamente localizaciones de hidrocarburo est basado en la
deteccin de magnetita diagentica, causada por la filtracin de hidrocarburo. La
35
presencia de cuerpos magnticos sobre acumulaciones de gas o petrleo ha sido
establecida para varias reas productoras del mundo, en yacimientos someros.
3.2.1 Campo Magntico Terrestre
La componente del campo magntico de la Tierra (entre el 80% y 90%) es
originada por la conveccin del hierro lquido en el ncleo externo del planeta
(Campbell, 1997), el cual es monitoreado por una red mundial de observatorios
magnticos y varios satlites (Langel y Hinze, 1997). Una aproximacin a este campo
dipolar sera de unos 50.000 nanoTesla (nT).
El campo magntico terrestre cambia lentamente con el tiempo,
constituyndose este fenmeno como la fuente de estudio del Paleomagnetismo.
Adems existe una influencia fuerte de campos externos, que influyen sobre el campo
terrestre, como es el caso de las tormentas solares (Nabighian et al., 2005).
3.2.2 Aplicaciones de las Medidas Magnticas
Existen diferentes formas de adquirir las mediciones magnticas: desde tierra,
aire, en los ocanos, en el espacio y dentro de los mismos pozos petroleros; cubriendo
un amplio rango de escalas y propsitos. Este tipo de levantamientos excepto el de los
pozos, son realizados con la premisa de identificar las anomalas en el campo
magntico producido por variaciones laterales en la magnetizacin de la corteza
terrestre; las medidas correspondientes a los pozos se enfocan en las variaciones
verticales de la magnetizacin en las vecindades del pozo en cuestin (Nabighian et
al., 2005).
36
De acuerdo con el Teorema de Poisson, el campo potencial de un cuerpo es
proporcional a la componente gravitatoria en la direccin de magnetizacin. Esta
relacin fundamental entre las propiedades magnticas y gravitacionales de un cuerpo
ha ayudado a la interpretacin de las anomalas magnticas (Cook, 1948).
Cook en 1948 expone que decenas de frmulas han sido desarrolladas con la
intencin de expresar la intensidad magntica sobre un cuerpo en funcin de la
susceptibilidad magntica, forma, dimensiones, disposicin del cuerpo, y la fuerza y
direccin del campo magntico terrestre, todas variables que afectan la respuesta
magntica de los cuerpos en el subsuelo.
3.2.3 Separacin Regional - Residual de Mapas Magnticos
La separacin regional- residual es un paso crucial en la interpretacin de los
datos magnticos en la deteccin de minerales a nivel somero, sin embargo para
aplicaciones petroleras su importancia decrece, ya que el rango en profundidad de la
exploracin de hidrocarburos se extiende a lo largo de una columna sedimentaria
(Nabighian et al., 2005). En todo caso, la respuesta regional de las anomalas
magnticas se atribuye al basamento gneo-metamrfico, mientras que las tendencias
locales corresponden a cuerpos magnetizados, ms prximos a superficie.
3.3 Mtodo Ssmico
Los diversos mtodos ssmicos de exploracin geofsica son actualmente las
herramientas ms utilizadas en la bsqueda indirecta de los hidrocarburos. El
fundamento de estos mtodos radica en que para litologas distintas, las ondas
corpreas se propagan con velocidades diferentes.
37
El mtodo de reflexin ssmica suministra informacin acerca de las
diferentes propiedades elsticas de las rocas y su objetivo principal es obtener
estimaciones de las posiciones de los estratos o estructuras del subsuelo al interpretar
las ondas ssmicas que se reflejan en ellas. Estas respuestas se miden en trminos de
las propiedades fsicas de amplitud y tiempo.
Por su parte Schlumberger (1980) define una seccin ssmica como un plano
vertical en el que aparecen representados los diferentes reflectores, que segn
Gonzlez (2007), no son ms que una serie de superficies de diferente impedancia
acstica, simulando capas sedimentarias.
3.3.1 Interpretacin Ssmica - Estructural
Segn Sheriff et al. (1990) la interpretacin ssmica se basa en dos aspectos
fundamentales: la continuidad de los reflectores y la correlacin. El primer aspecto
radica en la posibilidad de reconocer un pulso ssmico en trazas sucesivas, creando
alineaciones que dan la caracterstica de continuidad a lo largo de las secciones
ssmicas. Y el segundo permite hacer un reconocimiento de patrones como la
longitud, amplitud o forma de las ondas. Con ello pueden fcilmente correlacionarse
las propiedades de la traza ssmica con propiedades petrofsicas y geolgicas.
3.3.2 Calibracin Ssmica Pozo
La unidad vertical en la que se presenta una seccin ssmica, originalmente, es
tiempo doble de viaje, no profundidad. Es por ello, que los horizontes y fallas de
inters pueden ser interpretados en tiempo, posterior a una calibracin entre los datos
ssmicos y los registros de pozo. Este procedimiento se realiza mediante la
construccin del sismograma sinttico.
38
Trorey (1962) definen al sismograma sinttico como la reconstruccin de una
traza ssmica en la direccin del pozo. Para manejar matemticamente su significado es
necesario conocer la serie de reflectividad de las formaciones que son penetradas por
el pozo, y una ondcula que puede ser terica o extrada de la ssmica.
La serie de reflectividad o de coeficientes de reflexin representa los cambios
de impedancia acstica entre formaciones litolgicas cuya densidad y velocidad
varan con la profundidad. Por ello su expresin viene dada por la siguiente ecuacin:
i i
i i
i i i i
i i i i
i
Z Z
Z Z
V V
V V
CR
+
=
+
=
+
+
+ +
+ +
1
1
1 1
1 1
; Ec. 3.9
Donde:
CR = coeficiente de reflexin
= densidad de la formacin
V = velocidad intervlica
Z = impedancia acstica y
La capa i suprayace sobre la capa i+1
Segn esto, Schlumberger (1980) recomienda el empleo de los registros snico y
de densidad para determinar los coeficientes de reflexin (CR) correspondientes para
cada profundidad dentro del pozo. Para realizar un sismograma sinttico es preciso
calibrar las velocidades intervlicas, medidas con el registro snico, con las
velocidades promedio estimadas con los tiros de verificacin ssmica (TVS) o check
shots.
De acuerdo con Schlumberger (1980) para representar el sismograma sinttico
es necesaria la convolucin de un pulso terico o extrado de la ssmica con la serie
de reflectividad propia del pozo en cuestin (Figura 3.3).
39
Figura 3.3. Modelo convolucional del sismograma sinttico.
Las aplicaciones del sismograma sinttico, segn Trorey (1962), radican en los
siguientes aspectos:
Conversin exacta y continua de tiempo ssmico en profundidad.
Identificacin y correlacin de eventos ssmicos tanto de mltiples como
primarios.
Correlacin de la informacin ssmica con los registros de pozos.
3.3.3 Perfil ssmico vertical
Un perfil ssmico vertical, conocido por sus siglas en ingls VSP (vertical
seismic profile), consiste en inducir una seal ssmica desde la superficie, que es
grabada en tiempo por medio de sensores en un pozo espaciados a varios pies de
profundidad. El VSP es una evolucin de los TVS donde se analiza no slo el arribo
directo sino la onda emitida completamente. Esta tcnica guarda informacin acerca
de litologa, estructuras y comportamiento de las ondas en el subsuelo.
40
El campo de ondas presente en un VSP puede dividirse en dos grandes grupos
(figura 3.5): los eventos ssmicos que llegan al gefono desde la superficie (ondas
"D") y los que lo hacen desde abajo del receptor, tales como las reflexiones (ondas
"U").
Existen actualmente variantes que se podran clasificar como tipos de VSP
(Figura 3.4), las cuales son:
Zero-offset VSP: La fuente est ubicada en superficie a una distancia muy
corta del pozo, mientras que los receptores quedan espaciados, al igual que en
todos los VSPs, dentro del pozo.
Offset VSP: la fuente se encuentra en superficie a una distancia especfica del
pozo.
Walk-away VSP: la fuente se encuentra en superficie pero variando las
distancias entre la fuente y el pozo.
Azimuthal VSP: consiste en variar radialmente la posicin de la fuente con
respecto al pozo.
Tambin hay otros tipos de configuraciones en la adquisicin, pero no son
ms que combinaciones entre los diferentes arreglos de VSP descritos anteriormente.
41
Figura 3.4. Adquisicin de VSP, tomado de www.slb.com.
El registro final de VSP puede cumplir mltiples funciones, entre las cuales se
puede destacar: la calibracin de ssmica pozo, la identificacin de reflectores
dentro y por debajo del pozo y la conversin tiempo profundidad de alta precisin.
3.3.4 Conversin Tiempo Profundidad
Puesto que los horizontes ssmicos, fallas y otros rasgos estructurales son
interpretados comnmente en tiempo, es necesaria la conversin tiempo
profundidad (TZ) que permita manejar la informacin interpretada en magnitudes de
metros, pies, etc., de mayor relevancia en la industria de hidrocarburos.
Dicha conversin puede realizarse mediante varios mecanismos, que permitan
establecer una relacin directa entre los valores de tiempo y sus correspondientes
valores de profundidad. Para ello se pueden utilizar las velocidades provenientes del
procesamiento ssmico, las calibraciones de los sismogramas sintticos de la zona y
los tiros de verificacin ssmica.
42
Los tiros de verificacin ssmica se basan en la generacin de ondas ssmicas
en la boca de un pozo, el cual posee una serie de receptores ubicados a distintas
profundidades. stos registrarn los tiempos de viaje de las ondas (Figura 3.5),
permitiendo as la estimacin de las velocidades ssmicas medias (Dobrin, 1961).
Como resultado se obtiene una tabla de valores de tiempo-profundidad (curvas
TZ), que una vez verticalizados permitir establecer una relacin entre estas dos variables
(Schlumberger, 1980).
Finalmente, para obtener la conversin tiempo profundidad, es necesaria la
elaboracin de un modelo de velocidades que caracterice a la zona en estudio. ste
ser ms preciso mientras mayor sea la cantidad de datos empleados (curvas TZ,
velocidades RMS, sismogramas, etc.), siempre y cuando sea comprobada
previamente la calidad de los mismos.
Figura 3.5. Esquema de adquisicin de los registros de VSP. Tomado de
http://www.energia3.mecon.gov.ar.
3.4 Registros de Pozo
Los registros de pozo son utilizados en las diferentes etapas de vida de los
yacimientos, desde la exploracin hasta la produccin de los hidrocarburos. Existen
TVS: ONDA D (Directa)
43
actualmente una amplia gama de registros petrofsicos que contribuyen a la
interpretacin del subsuelo y a la extraccin del crudo, entre los cuales se pueden
resaltar: registro snico, registro de densidad, registro Gamma Ray, entre otros.
3.4.1 Registro Snico
El registro snico mide el tiempo de trnsito y la atenuacin de una onda
acstica que se propaga a travs de la formacin (Schlumberger, 1980). La
herramienta de adquisicin contiene dos trasmisores y cuatro receptores, tal como se
muestra en la Figura 3.6. La adquisicin del registro snico consiste en medir el
tiempo de refraccin de las ondas generadas en las distintas formaciones en
profundidad.
Figura 3.6. Esquema de herramienta de registro
snico. Tomado y modificado de Casas (2006).
44
Detector
de espaciamiento Grande
Detector
de espaciamiento corto
Enjarre
Fuente
Formacin
3.4.2 Registro de Densidad
Segn Casas (2006) los registros de densidad se obtienen por medio de una
sonda que se acopla a la pared del pozo. La herramienta de medicion consta de una
fuente emisora de rayos gamma colocada en la parte inferior de la sonda y dos
detectores de rayos gamma colocados a distancias diferentes de la fuente, tal como se
indica en la Figura 3.7.
Este tipo de registro tiene amplias aplicaciones, asociadas con la deteccin de:
litologas complejas, minerales pesados, formaciones gasferas, formaciones de baja
porosidad, fracturas y mineraloga de arcillas (Schlumberger, 1980).
Figura 3.7. Esquema de la herramienta de registro de
densidad. Tomado y modificado de Casas (2006).
3.4.3 Registro de Rayos Gamma
Esta herramienta define las concentraciones relativas de diferentes elementos
en la formacin, con el objeto de poder hacer una evaluacin completa de la misma
45
en trminos de saturacin, porosidad, arcillosidad, entre otras. La determinacin de
las concentraciones relativas de los elementos se basa en el anlisis de la energa de
los rayos gamma, que son inducidos en la formacin al bombardearla con neutrones
de alta energa (Schlumberger, 1980).
3.5 Anlisis Geoestadstico
La evaluacin de reservas minerales tiles, fue la actividad fundamental que
motiv a partir de la dcada de 1950, la aplicacin de la teora de funciones aleatorias
al reconocimiento y estimacin de fenmenos naturales. As surge la Geoestadstica,
trmino concebido por George Matheron a partir de trabajos previos de H. Sichel,
D.G. Krige, y B. Matern, que se desarrolla y consolida en los ltimos 30 aos como
ciencia aplicada y que da respuesta a necesidades prcticas y concretas dentro de la
caracterizacin de los yacimientos (Cuador, 2000).
La Geoestadstica, segn Cuador (2000), estudia las variables distribuidas
espacialmente, partiendo de una muestra representativa del fenmeno en estudio,
utilizando como elemento fundamental el anlisis de la distribucin espacial de la
informacin disponible. Se pretende entonces, minimizar la varianza del error de
estimacin, obtenindose el mejor estimador lineal insesgado, el Kriging, que es
actualmente utilizado en la mayora de las empresas mineras y petroleras alrededor
del mundo. Las herramientas ms comunes de la estadstica descriptiva para la
evaluacin de datos geofsicos son: las medidas de tendencia central, medidas de
dispersin, las distribuciones de frecuencia y las pruebas paramtricas y no-
paramtricas.
Medidas de dispersin: son conocidas como medidas de variacin e indican el
grado en que los valores numricos de un conjunto de datos tienden a
expandirse alrededor de un valor promedio (Spiegel y Stephens, 2002).
46
3,00 4,00 5,00 6,00 7,00 8,00 9,00 10,00
Anomala de Bouger (mGal)
0
1
2
3
4
5
6
7
F
r
e
c
u
e
n
c
i
a
Medidas de tendencia central: segn Spiegel y Stephens (2002) estas medidas
representan valores tpicos o representativos de un conjunto de datos, que
suelen situarse hacia el centro de su distribucin de probabilidad ordenados
por magnitud; entre los ms utilizados estn la moda, la mediana y la media
aritmtica.
Histograma de frecuencia: es una representacin grfica de una variable en
forma de barras, donde la superficie de cada barra es proporcional a la
frecuencia de los valores representados (Figura 3.8). En el eje vertical se
representan las frecuencias, y en el eje horizontal los valores de las variables,
normalmente sealando las marcas de clase; permitiendo tener una idea
grfica de cmo se distribuyen los datos.
Figura 3.8. Ejemplo de histograma de frecuencia. Datos
reales de Anomala de Bouguer para este proyecto.
Diagrama de caja: es una representacin grfica de la distribucin de un total
de datos (Figura 3.9), donde es posible reconocer los siguientes elementos: el
mnimo (medida ms pequea), el cuartil menor o primer cuartil (corta por
debajo del 25% de los datos), la mediana (dentro de la caja paralela al eje X),
47
el cuartil superior o tercer cuartil (que corta por encima del 75% de los datos)
y el mximo u observacin ms grande en magnitud.
Figura 3.9. Ejemplo de diagrama de caja. Datos reales de
Anomala de Bouguer para este proyecto.
En lneas generales, Gonzlez y Tabare (2007) describen que el diagrama de
caja divide a los datos en cuatro reas de igual frecuencia (Figura 3.9). Las lneas
verticales que parten de la caja son conocidos como bigotes, en donde el inferior tiene
un extremo en el primer cuartil Q
1
y el otro extremo en el valor dado por el primer
cuartil menos 0.5 veces la distancia intercuartil, esto es ( )
1 3 1
5 . 1 Q Q Q .
El bigote superior tiene un extremo en el tercer cuartil Q
3
y el otro en el valor
dado por el tercer cuartil ms 1.5 veces la distancia intercuartil (Figura 3.9), es decir
( )
1 3 3
5 . 1 Q Q Q + . Un valor se considera atpico si est a ms de 1.5 veces la distancia
intercuartil de alguno de los cuartiles Q
1
o Q
3
y aquellos a ms de 3 veces la distancia
intercuartil se les llama atpicos extremos (Gonzlez y Tabare, 2007).
Q
3
Q
1
Mediana
Anomala de Bouguer
CAPTULO 4:
MARCO METODOLGICO
Flujograma de Trabajo
Mtodo Gravimtrico
Mtodo Magntico
Mtodo Ssmicos
Registros de Pozo
Modelado Estructural
49
4.1 Flujograma de Trabajo
El desarrollo de este proyecto se inici con una primera etapa de recopilacin
bibliogrfica en la literatura especializada acerca de los rasgos estratigrficos,
tectnicos, estructurales, petrofsicos y ssmicos del rea, tanto a nivel regional como
local. As mismo, se investigaron los trabajos previos realizados en la zona en estudio
y los aspectos tericos que soportan el proyecto.
Posteriormente fue necesaria la realizacin de un inventario, a manera de
agrupar ordenadamente los datos gravimtricos, magnticos, registros de pozo, as
como la informacin ssmica 2D y 3D disponible de la zona en estudio. Ello permiti
el manejo de los datos y la interaccin entre ellos. En la figura 4.1 se muestra el
flujograma de trabajo seguido en el proyecto.
En el caso de los datos gravimtricos y magnticos se hizo especial nfasis en
el anlisis geoestadstico, para conocer la distribucin de sus poblaciones, as como
distinguir las mediciones anmalas que pudieron interferir con la calidad de los
resultados a ser mostrados en el prximo captulo.
A continuacin se presenta la metodologa utilizada para cada uno de los
mtodos geofsicos de prospeccin envueltos en la interpretacin integrada del campo
Sincor.
5
0
Flujograma de trabajo
Mtodo
Ssmico
Registro
de pozos
Mtodo
Magntico
Calibracin
ssmica-pozo
Elaboracin de mapas
estructurales en tiempo
Elaboracin de
modelo de velocidad
Interpretacin
estructural en tiempo
Elaboracin de
mapas estructurales en
profundidad
Mtodo
Gravimtrico
Elaboracin de
sismograma-sinttico
Estimacin de las
densidades de las
formaciones
Empleo de los
marcadores geolgicos en
tiempo y profundidad
Validacin y
correccin de datos
Generacin del
mapa de anomalas
de Bouguer
Separacin
regional-residual
Seleccin del
Perfil gravimtrico-
magntico
Validacin y
correccin de datos
Generacin del
mapa de IMT
Anlisis
Espectral
Modelo estructural del
Perfil A-A
Figura 4.1. Flujograma de trabajo.
51
4.2 Mtodo Gravimtrico
En este proyecto de investigacin se integr informacin proveniente de datos
gravimtricos obtenidos en diversos levantamientos, llevados a cabo durante la
dcada de 1980, dentro del plan de PDVSA de explorar las cuencas petrolferas de
Venezuela, que a travs de su filial Maraven evalu el potencial del bloque Zuata
(actual bloque Junn) de la Faja Petrolfera del Orinoco. Estos datos fueron sometidos
a un control de calidad que permiti corregirlos y validarlos.
El nmero total de estaciones gravimtricas en el rea es de 92, distribuidas de
manera heterognea en la zona (Figura 4.2), en su mayora a lo largo de la carretera
que conduce a la poblacin de San Diego de Cabrutica (direccin NO - SE). La
Gerencia Tcnica de Exploracin de INTEVEP (Filial de PDVSA) fue quien
proporcion la base de datos (Apndice 1) con la ubicacin (coordenadas
geogrficas), cota, los valores de gravedad observada, anomala de aire libre y
anomala de Bouguer para cada estacin (Figura 4.2).
Para validar los datos gravimtricos se procedi a cotejar los valores de cota
de cada estacin medidos durante las campaas de adquisicin de 1980, con sus
homlogos correspondientes a los levantamientos ssmicos y los mapas topogrficos
de la zona.
Empleando las ecuaciones 3.1, 3.2 y 3.3 se calcularon los valores de gravedad
terica, correccin de Bouguer y correccin de aire libre respectivamente, empleando
el programa Microsoft Excel
,
generada por el Departamento de Geofsica de la UCV. Luego, se procedi a
recalcular los valores de anomala de Bouguer usando la ecuacin 3.5 con una
densidad de Bouguer de 2,1 g/cm
3
; valor obtenido mediante un perfil de Nettleton en
la zona elaborado por la compaa Suelopetrol (Colmenares, 2007).
Posteriormente, se cartografi un primer esbozo del mapa de anomala de
Bouguer en el rea, mediante el uso del programa Surfer
Versin
12.0, en donde se efectuaron histogramas de frecuencia, diagramas de caja, pruebas
paramtricas P-P y Q-Q, as como los clculos de los valores de tendencia central y
de dispersin.
El histograma de frecuencia correspondiente a una muestra de 92 datos de
anomala de Bouguer es mostrado en la figura 4.3, y se presenta con un
comportamiento, que en lneas generales, podra categorizarse como unimodal. En
cuanto a la distribucin de los valores de anomala de Bouguer, puede interpretarse
que no presentan normalidad, ya que los valores de moda, mediana y media presentan
magnitudes distantes entre s.
54
Figura 4.3. Histograma de frecuencia del total de datos de anomala de Bouguer.
En la Figura 4.3 se presenta la distribucin de la totalidad de datos de
anomala de Bouguer (92). Los datos constituyen una sola poblacin, en donde se
reconoce a un valor aislado del resto (enmarcado en rojo), correspondiente al valor
mnimo de la muestra, que es mostrado al igual que las medidas de tendencia central
y dispersin para el total de datos de anomala de Bouguer en la tabla 4.1.
Tabla 4.1. Medidas de tendencia central y dispersin para los valores de
anomala de Bouguer. Muestra correspondiente a 92 estaciones
gravimtricas.
Estadstico mGal
Media 2,2374
Mediana 1,7418
Moda 2,0189
Desviacin tpica 0,2497
Varianza 8,273
Mnimo -6,58
Mximo 9,99
El diagrama de caja para el total de datos gravimtricos es mostrado en la
Figura 4.4, donde se aprecia que por arriba del tercer cuartil se encuentran varios
10,000000 5,000000 0,000000 -5,000000 -10,000000
Anomala de Bouguer (mGal)
30
25
20
15
10
5
0
F
r
e
c
u
e
n
c
i
a
55
datos fuera de la mayor distribucin, presentados con un crculo y su respectiva
etiqueta de cada estacin, significando que son mediciones atpicas.
_
Figura 4.4. Diagrama de caja para los datos de anomala de Bouguer. Encerrado en
el crculo rojo queda reflejado lo excluida que est la estacin 84 del resto de la
poblacin de datos.
La estacin 84, resaltada en la Figura 4.4 en rojo, se presenta como el nico
dato atpico extremo negativo de toda la muestra que se encuentra totalmente
relegado de la mayor distribucin de datos. Por esta razn se procedi a realizarse un
estudio de mayor detalle evaluando esta estacin. Para este fin se solaparon las
estaciones gravimtricas sobre los mapas geolgico (Figura 4.5) y topogrfico de la
zona, con el objetivo de indagar el posible efecto de la geologa superficial y el
relieve de la zona en estudio sobre la estacin 84. Los mapas mencionados son los
siguientes.
Anomala de Bouguer
10,000000
5,000000
0,000000
-5,000000
-10,000000
20
55
84
56
265000 270000 275000 280000 285000 290000 295000 300000 305000 310000 315000
895000
900000
905000
910000
915000
920000
925000
930000
935000
940000
San Diego de Cabrutica
Escala 1:300.000
m
Mapa de Geologa Superficial
con Ubicacin de las
Estaciones Gravimtricas
UNIVERSIDAD CENTRAL
DE VENEZUELA
FACULTAD DE INGENIERA
ESCUELA DE GEOLOGA,
MINAS Y GEOFSICA
DEPARTAMENTO DE
GEOFSICA
DAVID CONTRERAS
HERMAN ROJAS
2007
ELABORADO POR:
A'
A
Perfil A - A'
Leyenda
Estacin 84
Poblado
Campo Sincor
0 5000 10000 15000
Fm. Mesa
Fm. La Pica
Fm. Las Piedras
Venezuela. Direccin de Cartografa Nacional. (1983) San Diego de
Cabrutica. Hoja 7241. Caracas. Escala 1:100.000. Mapa geolgico.
Venezuela. U.S Geological Survey, FUNVISIS y UCV. (2005) Venezuela.
Hoja Norte. Escala 1:750.000. Mapa Geolgico del Relieve de Venezuela.
Figura 4.5. Mapa de geologa superficial con ubicacin de las estaciones gravimtricas. La estacin
84 se representa con verde para diferenciarla del resto de estaciones mostradas en rojo. Proyeccin:
UTM, Huso 20; Datum Geodsico: La Canoa.
La estacin 84 fue adquirida en una zona donde aflora el Cuaternario reciente
(Figura 4.5), representado por las arenas de la Formacin Mesa. A nivel de geologa
de superficie, no hay ningn rasgo fisiogrfico al que pudiera otorgrsele la causa por
la que el valor de la anomala de Bouguer sea tan elevado en magnitud, cercano a los
-6,58 mGal (Apndice 1), respecto al resto de la muestra.
57
10,000000 7,500000 5,000000 2,500000 0,000000 -2,500000
Anomala de Bouguer
20
15
10
5
0
F
r
e
c
u
e
n
c
i
a
Mientras tanto, al comparar los valores reportados de cota durante la
adquisicin de esta estacin con la cota asociada en el mapa topogrfico, se presenta
una diferencia de aproximadamente 40 m, razn por la cual, se imposibilita que la
correccin topogrfica sea efectiva.
Una vez excluida la estacin 84 se procedi a graficar nuevamente el
histograma de frecuencia (Figura 4.6) y el diagrama de caja (Figura 4.7) para la
poblacin de datos de anomala de Bouguer. El comportamiento de normalidad en la
distribucin de los datos fue evaluado por medio de pruebas paramtricas P-P y Q-Q
(Figura 4.8).
Figura 4.6. Histograma de frecuencia para los datos de anomala de Bouguer excluida la estacin 84.
58
Anomala de Bouguer
10,00
7,50
5,00
2,50
0,00
-2,50
20
43
Figura 4.7. Diagrama de caja para los datos de anomala de Bouguer excluida la estacin 84.
En la Figura 4.6 se muestra el histograma correspondiente a los valores de
anomala de Bouguer luego de haber excluido la estacin 84, donde se puede notar
una mejor distribucin de los datos, con tendencia a la normalidad. Respecto a la
Figura 4.7 se debe decir que se presentan ciertos datos lejos de la mayor densidad de
valores, que a pesar de ser atpicos, no hay justificacin geolgica para su exclusin.
A continuacin se muestra la tabla 4.2 con la estadstica descriptiva de los
datos de anomala de Bouguer, para un total de 91 estaciones, luego de haber sido
excluida la estacin 84, lo cual genera un nuevo mnimo para el total de datos.
59
Tabla 4.2. Medidas de tendencia central y dispersin para los valores de
anomala de Bouguer luego de la exclusin de estacin 84. Muestra
correspondiente a 91 estaciones gravimtricas.
Estadstico mGal
Media 2,3343
Mediana 1,8097
Moda 2,1283
Desviacin tpica 2,82959
Varianza 7,491
Mnimo -2,11
Mximo 9,99
Las pruebas paramtricas P-P y Q-Q permitieron hacer el anlisis del
comportamiento normal del total de datos de anomala de Bouguer, antes de la
exclusin de la estacin 84 (Figuras 4.8a y 4.8b) y despus de haber excluido la
estacion 84 (Figuras 4.8c y 4.8d). Hay que denotar que los datos siguen presentando
distribucin normal, al acoplarse con la recta de 45 en cada uno de los casos, donde
la exclusin de la estacin 84 permite ver grficamente dicha distribucin para una
muestra menor.
6
0
Figura 4.8. Anlisis de la distribucin de los datos gravimtricos. a) Prueba paramtrica P-P para el total de datos de
anomala de Bouguer. b) Prueba paramtrica Q-Q para el total de datos de anomala de Bouguer. c) Prueba paramtrica
P-P para los datos de anomala de Bouguer sin estacin 84. d) Prueba paramtrica Q-Q para los datos de anomala de
Bouguer sin la estacin 84.
61
4.2.2 Generacin del Mapa de Isoanomalas de Bouguer
Una vez concluido el anlisis estadstico, se gener el mapa de anomala de
Bouguer definitivo, utilizando el programa Surfer
Versin 12.0.
En la Figura 4.13 se presenta el histograma de frecuencia de los datos de
Intensidad Magntica Total (IMT) del rea, donde se evidencia una sola poblacin
que presenta distribucin parcialmente sesgada hacia los valores de mayor frecuencia.
Figura 4.13. Histograma de frecuencia del total de datos magnticos.
A continuacin se presenta la tabla 4.3 donde son mostrados los valores de
tendencia central y dispersin para el total de datos de IMT. Acorde con estos datos
fue generado el diagrama de caja y bigotes mostrado en la Figura 4.15.
67
Tabla 4.3. Medidas de tendencia central y dispersin. Muestra
correspondiente a 6473 valores de Intensidad Magntica Total.
Estadstico gammas
Media 9795,9669
Mediana 9802,54
Moda 9842,53
Desviacin tpica 58,2691
Varianza 3395,288
Mnimo 9668,74
Mximo 9968,02
El diagrama de caja correspondiente al total de datos de IMT (Figura 4.14)
comprueba la existencia de una sola poblacin de datos, donde el nico punto atpico
es la estacin 2474 con 9968,02 gammas, que corresponde al valor mximo de la
muestra. Se decide entonces excluir a esta estacin como dato inicial para la
elaboracin del mapa de IMT, que ser presentado y analizado en el prximo
captulo. La mayor parte de la muestra est en el rango de 9755 y 9840 gammas.
Figura 4.14. Diagrama de caja para el total de datos magnticos.
La normalidad de los datos de IMT queda comprobada al apreciar se ajustan a
la recta de 45 en el grfico que representa la prueba paramtrica P-P, que es
68
0,0 0,2 0,4 0,6 0,8 1,0
Probabilidad Acumulativa Observada
0,0
0,2
0,4
0,6
0,8
1,0
P
r
o
b
a
b
i
l
i
d
a
d
A
c
u
m
u
l
a
t
i
v
a
E
s
p
e
r
a
d
a
9.400 9.600 9.800 10.000 10.200
Valor Observado
9.400
9.600
9.800
10.000
10.200
V
a
l
o
r
N
o
r
m
a
l
E
s
p
e
r
a
d
o
mostrada en la Figura 4.15. Lo mismo sucede con la prueba paramtrica Q-Q (Figura
4.16).
Figura 4.15. Prueba paramtica P-P para el total de datos magnticos.
Figura 4.16. Prueba paramtrica Q-Q para el total de datos magnticos.
69
9650
9700
9750
9800
9850
9900
9950
0 5 10 15 20 25 30 35 40
Distancia (km)
I
M
T
(
g
a
m
m
a
s
)
S
A A
N
N10O
4.3.2 Generacin del Mapa de Isoanomalas de IMT
Finalizado el anlisis estadstico de los datos magnticos se procedi a
elaborar el mapa de isoanomalas de IMT (Figura 5.5), que ser presentado y
discutido en el siguiente captulo. Para este fin se implement el programa Surfer
Versin 8.0. El perfil magntico a modelar es el mismo A-A descrito para el mtodo
gravimtrico y se puede observar en la Figura 4.17, como en lneas generales los
valores de IMT disminuyen hacia el norte, mostrando estabilidad para los 9825
gammas entre los kilmetros 15 y 25 del perfil.
Figura 4.17. Intensidad Magntica Total para el Perfil A A.
Puesto que la respuesta magntica dentro del campo Sincor es producida
presumiblemente en su totalidad por el basamento gneo-metamrfico, se considera
que la cobertura sedimentaria tiene un efecto magntico despreciable. Es por esta
razn que no tiene sentido geolgico-geofsico realizar la separacin regional
residual de la seal magntica, dado que todo el efecto tiene carcter regional. Sin
embargo, se realiz el clculo de la bondad de ajuste de las superficies polinmicas
(Figura 4.18), para corroborar la anterior afirmacin.
70
50
55
60
65
70
75
80
85
90
95
0 1 2 3 4 5 6
Grado del Polinomio
B
o
n
d
a
d
d
e
A
j
u
s
t
e
(
%
)
Figura 4.18. Grfico de Bondad de ajuste para datos de IMT.
En la Figura 4.18 se evidencia efectivamente, que ninguna superficie
polinmica logra modelar matemticamente la respuesta magntica de la zona, ya que
el ajuste de la respuesta de IMT aumenta con el grado del polinomio utilizado. Por
esta razn no se realiz la separacin regional residual correspondiente a los datos
aeromagnticos.
4.3.3 Anlisis Espectral de los Valores de IMT
El anlisis espectral aplicado para los valores de IMT fueron realizados con la
misma metodologa llevada a cabo en el mtodo gravimtrico, utilizando nuevamente
el programa Surfer
y Gocad
de la plataforma OpenWorks
.
Posteriormente, se estim la frecuencia mxima y mnima considerando una
atenuacin de 12 decibeles (dB), lo que implica una frecuencia dominante a de la
amplitud mxima. La frecuencia dominante fue estimada por medio de la siguiente
ecuacin:
+
=
2
fmin fmax
fdom , Ec. 5.1
Donde:
fdom= frecuencia dominante
fmax = frecuencia mxima y
fmin= frecuencia mnima
75
Figura 4.20. Espectro de amplitud ssmica en la cercana del pozo P6.
4.4.2 Conversin Tiempo Profundidad
La conversin a profundidad de los horizontes ssmicos interpretados en
tiempo, se realiz mediante la construccin de un modelo de velocidades, generado a
partir de los tiros de verificacin ssmica en los 61 pozos mostrados en la Tabla 4.6.
Tabla 4.6. Inventario de pozos con tiros de verificacin ssmica
Nombre de Pozo Check Shot Nombre de Pozo Check Shot
P1 X P46 X
P2 X P47 X
P3 X P48 X
P4 X P49 X
P5 X P50 X
P6 X P51 X
P8 X P52 X
P9 X P53 X
P10 X P54 X
P15 X P55 X
P26 X P56 X
fmin fmax
76
Continuacin Tabla 4.6
Nombre de Pozo Check Shot Nombre de Pozo Check Shot
P27 X P57 X
P28 X P58 X
P29 X P59 X
P30 X P60 X
P31 X P61 X
P32 X P62 X
P33 X P63 X
P34 X P64 X
P35 X P65 X
P36 X P66 X
P37 X P67 X
P38 X P68 X
P39 X P69 X
P40 X P70 X
P41 X P71 X
P42 X P72 X
P43 X P73 X
P44 X P74 X
P45 X P75 X
Se emple el mtodo DSI (discrete smooth interpolation o interpolacin suave
discreta) de la plataforma Gocad
de la plataforma OpenWorks
.
La Tabla 4.7 muestra la totalidad de los pozos utilizados para la calibracin
ssmica pozo (con excepcin de los pozos P24 y P25 usados como control de
calidad del modelo). Como requisito cada uno de ellos deba que tener los registros
RHOB y DTC.
Para las calibraciones fueron empleadas las curvas T-Z derivadas de los pozos
que contenan check shots. Para los pozos que solo tenan registros VSP, las curvas T-
Z fueron extradas del campo descendente de dichos registros; mientras que para los
pozos que no contienen ninguno de estos dos registros se utiliz la curva T-Z del
pozo ms cercano. En la Figura 4.24 se muestra las pruebas de calibracin de los
sismogramas sintticos con las diferentes ondculas.
81
Figura 4.24. Prueba de calibracin ssmica del Pozo P6. A la derecha se observa el
sismograma sinttico a partir de la convolucin con la ondcula extrada de la ssmica con un
coeficiente de correlacin de 57%. A la izquierda el sismograma sinttico a partir de la
convolucin con ondcula Ricker con un coeficiente de correlacin de 65%.
En la figura 4.24 se aprecian las calibraciones de los sismogramas (para cada
ondcula utilizada) con la ssmica, seleccionndose finalmente aquellos obtenidos
mediante la convolucin de la ondcula Ricker por presentar mejores coeficientes de
correlacin. Posteriormente se aplicaron a los sismogramas sintticos elegidos los
filtros Ricker y trapezoidal, ajustando su contenido frecuencial con las frecuencias de
la ssmica en la zona de inters.
Los filtros Ricker fueron construidos al centrar las ondculas fase cero en las
frecuencias dominantes de la ssmica, mientras que para elaborar los filtros
82
trapezoidales se analizaron los espectros de poder en donde se definieron las
frecuencias de paso y las frecuencias de corte para -12dB.
El filtro Ricker fue el elegido para realizar las calibraciones definitivas, dado
que los coeficientes de correlacin entre los sismogramas sintticos y la ssmica
presentaron mayor magnitud para una ventana seleccionada (figura 4.25).
Figura 4.25. Pruebas de filtros Ricker en los sismogramas sinttico y VSP del pozo P6.
Para realizar las calibraciones de los sismogramas sintticos con la ssmica
(Figura 4.26), se correlacion desde el tope de la Formacin Oficina hasta el ltimo
marcador en profundidad, considerando un rango de 10 trazas ssmicas alrededor del
pozo. En general, se busc que existiera una buena similitud entre los eventos de los
sismogramas sintticos y las reflexiones ssmicas, variando empricamente el tiempo
de reajuste (time shift) hasta conseguir los valores ptimos para la calibracin.
8
3
Figura 4.26. Ejemplo de hoja de calibracin del pozo P6
84
4.5.2 Calibracin entre los VSP y los Datos Ssmicos
Para cargar los datos de VSP en la plataforma OpenWorks
, fue necesaria la
conversin del formato SEGY a un formato ASCII, que permiti su carga como
registros de pozo. Esta conversin no la efecta OpenWorks
directamente, por lo
que fue necesario programar un cdigo en Matlab
de la plataforma OpenWorks
para
la consecucin de esta tarea.
4.6 Modelado Estructural
A travs de la construccin del modelo estructural de la zona en estudio,
correspondiente al perfil A-A, se pretende estimar la profundidad y geometra de los
diferentes estratos sedimentarios presentes en el subsuelo. Para ello se procedi en
primer lugar a cargar en el programa GM_SYS