Un modelo 3D de la densidad del subsuelo en el graben de Espino: un aulacógeno Jurásico en el oriente de Venezuela

 

A 3D model of the subsurface density in the Espino graben: a Jurassic aulacogen in eastern Venezuela

 

Juan García-Abdeslem¹*, Jimmy José Regalado-Sosa² y Hugo Rafael Cerquone-Ravelo²

 

¹ Centro de Investigación Científica y de Educación Superior de Ensenada (CICESE), División de Ciencias de la Tierra, Departamento de Geofísica Aplicada, Carretera Ensenada-Tijuana 3918, Zona Playitas, 22860, Ensenada, Baja California, México. *jgarcia@cicese.mx

² Petróleos de Venezuela S.A., Exploración y Producción, Gerencia Corporativa de Geofísica, Puerto La Cruz, Anzoátegui, Venezuela.

 

Manuscrito recibido: Septiembre 19, 2012
Manuscrito corregido recibido: Diciembre 10, 2012
Manuscrito aceptado: Diciembre 17, 2012

 

RESUMEN

Como parte de un programa de generación de oportunidades de exploración petrolera en el oriente de Venezuela, hemos obtenido un modelo 3D de la densidad del subsuelo en la región del graben de Espino mediante la inversión lineal de anomalías de gravedad de Bouguer residual isostática. En el modelo obtenido se interpreta que la superficie cuyo contraste en densidad es cero corresponde al basamento económico y sugiere un relleno sedimentario de ~6500 m de espesor en la región central del graben. Sobre esa superficie se observa una región cóncava hacia arriba, de ~120 km de largo por ~80 km de ancho, con paredes inclinadas de manera consistente con las fallas que delimitan el graben. Se interpreta que la región cóncava, donde el contraste en densidad es negativo, contiene el mayor espesor sedimentario sintectónico del Jurásico (~3500 m). La porción más somera del modelo sugiere un relleno con secuencias sedimentarias poco deformadas de ~3000 m de espesor que interpretamos fueron depositadas después del periodo de extensión cortical.

Palabras clave: anomalía gravimétrica, inversión lineal 3D, graben de Espino, Venezuela.

 

ABSTRACT

As apart of aprogram to generate opportunities for oil exploration in eastern Venezuela, we have obtained a 3D subsurface density model over the Espino graben region, by the linear inversion of isostatic residual Bouguer gravity anomalies. The obtained model suggests that the sedimentary column is at least 6500 m thick in the central portion of the graben. We interpreted the model zero density contrast surface as the economic basement, and on that surface we observe an upward concave region with north and south edges inclined consistently with the border faults of the graben. It is interpreted that within the concave region, where the density contrast is negative, the graben accommodates some 3500-m thick sedimentary sequences that were deposited during an episode of crustal extension in Jurassic times. The shallow portion of the model suggests 3000-m thick, sub horizontal strata, interpreted as sedimentary sequences deposited after the episode of crustal extension.

Key words: gravity anomalies, 3D lineal inversión, Espino graben, Venezuela.

 

INTRODUCCIÓN

Se ha sugerido que el graben de Espino en el oriente de Venezuela (Figura 1) es un aulacógeno que se originó como un rift abortado de una unión triple RRR durante la apertura del Atlántico central en el Jurásico, con posibilidad de albergar un sistema petrolero en rocas pre-cretácicas (Burke, 1977; Moticska, 1985; Salazar-Tomey, 2006). Diversos estudios de reflexión sísmica y de estratigrafía sísmica por secuencias, aunque no alcanzan a iluminar claramente la base del graben, han mostrado que el graben de Espino se extiende en dirección NE con una longitud de unos 120 km, y está limitado hacia el noreste por las fallas Machete, Ruiz-Sabán y Guama-Sabán, que son fallas normales con deslizamiento al SE; hacia el sureste, el graben está limitado por la falla normal Altamira con deslizamiento al NW (Figura 1).

Con base en pozos perforados en la región, se sabe que el basamento cristalino al sur de la Falla Altamira es del Precámbrico y que sobre ese basamento yacen rocas paleozoicas de las Formaciones Hato Viejo y Carrizal del Cámbrico temprano y medio, respectivamente, y se infiere que la base del graben está ocupada por rocas afines a estas formaciones.

La Formación Hato Viejo (Hedberg, 1942, en PDVSA-Intevep, 2012) está caracterizada principalmente por arenisca de grano fino a grueso y conglomerado con clastos de cuarcita a intervalos irregulares; yace discordantemente sobre las rocas precámbricas del basamento del cratón de Guayana y subyace discordantemente a la Formación Carrizal. Su espesor es de unos 90 m y desaparece al este del Estado Anzoátegui, donde el Grupo Temblador es discordante sobre rocas precámbricas del cratón de Guayana. El ambiente de sedimentación es continental y sus sedimentos representan el relleno de la cuenca en facies fluvial y/o de piedemonte.

La Formación Carrizal (Hedberg, 1942, en PDVSA-Intevep, 2012) está constituida por una secuencia de arcilita, roja, dura y compacta. Contiene algunas capas de limolita y arenisca, generalmente bioturbadas, y hacia la base de la formación hay estratos de arenisca. El tope de la Formación Carrizal es siempre erosional y está cubierto por sedimentos pre-cretácicos del Grupo Temblador (Hedberg, 1942, en PDVSA-Intevep, 2012) o de la Formación Oficina (Hedberg et al., 1947, en PDVSA-Intevep, 2012), y sobreyace a la Formación Hato Viejo.

De acuerdo con la interpretación estratigráfica de los pozos profundos NZZ-88 y TM-45 (Figura 1), el graben de Espino acomoda rocas sedimentarias clásticas del Paleozoico al Jurásico, deformadas y dislocadas durante el proceso de rifting del Jurásico que inició la apertura del Océano Atlántico, cuyo final está marcado por una discordancia de extensión regional sobre la cual yacen rocas sedimentarias clásticas cretácicas tardías a neógenas. Éstas incluyen al Grupo Temblador del Cretácico, caracterizado por arenisca con intercalaciones delgadas de lutita, marga y caliza, y a la Formación Merecure (Funkhouser et al., 1948, en PDVSA-Intevep, 2012) del Oligoceno y a la Formación Oficina del Mioceno. Estas dos últimas caracterizadas por una alternancia de lutita, arenisca y limolita.

En la región del graben de Espino, el Basalto Altamira es la única evidencia de actividad volcánica en el Jurásico, con edad absoluta de 160 ± 5 Ma (Moticska, 1985). Este basalto fue descubierto a ~2000 m de profundidad al perforar el pozo NZZ-88, justo al norte de la Falla Altamira (Figura 1). Está constituido por tres flujos de lava que en conjunto alcanzan un espesor de ~50 m. Estos basaltos representan un episodio volcánico aislado en la fase terminal del rifting, y de acuerdo con sus relaciones de contacto, se interpreta que fueron emitidos en condiciones subaéreas sobre estratos rojos (MARAVEN, 1983; Moticska, 1985).

La estructura regional y la profundidad al límite corteza-manto en Venezuela ha sido inferida con base a la interpretación de datos derivados de experimentos de refracción sísmica profunda, de reflexión sísmica de gran ángulo, y de eventos telesísmicos analizados mediante la técnica de funciones de receptor (Schmitz et al., 2002 y 2008; Bezada et al., 2008; Fenglin et al., 2007). Dichos estudios han mostrado que en la región de la Cuenca Oriental, donde el graben de Espino está ubicado (Figura 1), la corteza continental alcanza un espesor de ~40 km, lo cual sugiere que en la región del graben de Espino la corteza es al menos 5 km más gruesa que el grosor cortical promedio de la corteza continental y unos 5 km más delgada que el grosor de corteza estimado para el cratón de Guayana, al sur del Río Orinoco (Figura 1), donde se ha estimado un grosor cortical entre 42 y 46 km (Schmitz et al., 2002). Con base en esas diferencias de espesor cortical, se puede interpretar que el proceso de extensión que originó el graben de Espino durante el Jurásico, afectó una corteza continental gruesa y térmicamente madura (v.g., fría), equivalente a la corteza del cratón de Guayana, con una rigidez flexural suficientemente grande como para soportar el desbalance de masa debido al proceso de extensión, con adelgazamiento cortical marginal y escaso volcanismo subaéreo. Esta interpretación concuerda con Salazar-Tomey (2006), quien infiere un estiramieto cortical (β) de 1.13 a 1.16 en la región del graben de Espino.

En este trabajo presentamos la metodología utilizada y el resultado de la modelación de anomalías de gravedad residual isostática de la región del graben de Espino. La modelación de los datos gravimétricos se realizó utilizando un algoritmo de inversión lineal, cuyo diseño se basa en métodos descritos por Jackson (1979) y Tarantola (2005), para inferir el contraste en densidad en un cuboide constituido por un conjunto de prismas rectangulares, de forma tal que su efecto gravimétrico explique la gravedad observada. El resultado de la modelación inversa sugiere que en la región del graben de Espino hay una columna sedimentaria de al menos 6500 m de espesor.

 

DATOS GRAVIMÉTRICOS

Para estudiar la región del graben de Espino utilizamos datos de anomalías de gravedad de aire libre de una compilación regional que incluye una selección de los datos de mejor calidad disponibles en tierra y datos derivados de altimetría satelital (Andersen et al., 2010; Technical University of Denmark, 2012). Calculamos la anomalía de Bouguer completa utilizando datos de topografía de la Shuttle Radar TopographicMission (<http://strm.csi.cgiar.org/>) y batimetría de la General Bathymetric Chart of the Oceans (<http:///www.gebco.net>). Este cálculo se extendió a una distancia de 167 km alrededor del área de estudio, utilizando una densidad de 2000 kg/m³, para los sedimentos que cubren las tierras bajas (< 500 m SNM) de la Cuenca Oriental y de la Cuenca Barinas-Apure, un valor consistente con los registros de densidad disponibles, y 2670 kg/m³ para la Serranía del Interior, la Cordillera de la Costa y el cratón de Guayana (Figura 1). A la anomalía de Bouguer calculada se le restó un efecto regional, calculado a partir de la raíz cortical derivada de un modelo isostático de tipo regional, cuyo cálculo se describe en el Apéndice A. La anomalía residual isostática obtenida se muestra en la Figura 2.

Vale la pena mencionar que la secuencia del procesamiento de los datos gravimétricos implica que la anomalía de Bouguer residual isostática debe ser interpretada como debida al contraste de densidad entre las rocas objeto del estudio y una densidad de referencia que típicamente se designa como 2670 kg/m³. Una narrativa de aspectos históricos sobre el valor asignado típicamente a la densidad de referencia se puede consultar en Hinze (2003).

 

EL PROBLEMA DIRECTO

La geometría del modelo es un cuboide que consiste de M_e estratos. Cada estrato tiene un espesor constante y está constituido por una malla regular de M_r × M_c prismas, y en cada prisma la densidad es constante.

En un punto de coordenadas x_i, y_i, z_i, localizado fuera del cuboide, la componente vertical del campo de gravedad está dada por

en donde ρ_j es la densidad de j-ésimo prisma, M=M_e×M_c×M_r es el número de prismas en el cuboide, y

donde Γ es la constante gravitacional, R=(X²+Y²+Z²)^1/2, X=x_j-x_i,Y=y_j-y_i, y Z=z_j-z_i en donde x_j, y_j, z_j son coordenadas de la fuente, y el j-ésimo prisma está limitado entre X¹_j < x < X²_j, Y¹_j< y < Y²_j y Z¹_j < z < Z²_j donde x y y representan distancia horizontal y z es profundidad, definida positiva hacia abajo.

 

EL PROBLEMA INVERSO

El objetivo del problema inverso es inferir el contraste de densidad del subsuelo en el dominio de un cuboide rectangular, cuyo efecto gravitacional explique la anomalía de gravedad observada. Se busca que el cuboide represente la estructura del graben y que en general sea consistente con los registros de densidad obtenidos en pozos que PDVSA ha perforado en el graben de Espino.

Los parámetros del modelo, v.g., la densidad de los prismas en el cuboide, se describe como m ∈ E^M, un vector en el espacio Euclidiano de dimensión M, donde M corresponde al número de prismas en el cuboide. Similarmente, el vector g ∈ E^N denota la solución al problema directo, donde N es el número de datos, y el vector d ∈ E^N representa la gravedad observada.

Suponiendo que asociado a cada dato tenemos una estimación de su error σ_i (i =1,..., N) y que los errores son estadísticamente independientes, incluimos la incertidumbre en los datos en la diagonal principal de su matriz de covarianza, C_d = σ_iI. Proponemos a priori un modelo inicial, m_prior, que consiste en una estimación de la densidad en cada uno de los estratos que constituyen el cuboide, y una estimación de su incertidumbre σ_j (j =1, M), la cual se incluye en la matriz de covarianza del modelo, C_m = σ_jI. Adicionalmente, buscamos una solución cuya variación espacial sea suave, incluyendo un operador de primera derivada, D, que actúa sobre el modelo en las direcciones x, y, z. De acuerdo con Jackson (1979) y Tarantola (2005), la solución al problema inverso lineal está dada por

en donde α es un escalar que asigna importancia al acotamiento en suavidad y A representa la matriz de derivadas del modelo directo con respecto a la densidad, es decir, la solución al problema directo para densidad unitaria.

A fin de optimizar los recursos computacionales disponibles para resolver el problema inverso, las coordenadas de los datos de gravedad y de topografía fueron rotadas 14 grados en sentido horario, y se seleccionó la región rectangular de 160 por 110 km indicada en la Figura 2. El problema inverso incluye 17871 datos de gravedad distribuidos regularmente en una malla con nodos espaciados 1 km, y postulamos que la incertidumbre en los datos es de 1 mGal. La geometría del modelo del subsuelo es un cuboide que incluye 12240 prismas repartidos en 15 estratos (M_e = 15) de espesor constante; cada estrato está constituido por un arreglo de M_r = 24 por M_c = 34 prismas, y la extensión horizontal de cada prisma se fijó en 5×5 km². En el modelo a priori (Figura 3, Tabla 1), el contraste en densidad del estrato más somero se asignó considerando que en la región de estudio afloran sedimentos pobremente consolidados; a los estratos más profundos se les asignó un contraste en densidad que disminuye con la profundidad y que, a grosso modo, sigue la variación observada en los registros de densidad de los pozos perforados por PDVSA en la región.

Los resultados de la inversión se presentan en las Figuras 4 y 5. En la Figura 4 se compara la anomalía de gravedad observada con la anomalía calculada mediante isolíneas, que en general se traslapan pues el desajuste es menor a ±1 mGal. En la Figura 5, la profundidad (z) es negativa y se muestra en hectómetros (1 hm = 100 m). La Figura 5a muestra cortes transversales del modelo en las direcciones x y z, donde se aprecia la variación lateral del contraste en densidad a profundidad, sugiriendo que en algunos segmentos el graben pudiera extenderse hasta ~11 km de profundidad. La Figura 5b muestra la superficie de cero contraste en densidad, la cual se interpreta como el basamento económico. Por debajo de esa superficie se incrementa la densidad, implicando que disminuyen la porosidad y la permeabilidad, y no esperamos encontrar las condiciones para el desarrollo de un sistema petrolero.

 

VALIDACIÓN DEL MODELO

Para validar la solución al problema inverso le sumamos 2670 kg/m³ al modelo 3D de contraste en densidad para compararlo con los registros de densidad de diez pozos. El modelo 3D de densidad fue muestreado en la vecindad de los pozos para obtener perfiles verticales de densidad comparables con los registros de densidad que se observan en la Figura 6. Esta comparación muestra que la densidad del modelo 3D es comparable en magnitud y tendencia en profundidad con los registros de densidad.

 

CONCLUSIONES

Mediante la inversión lineal de anomalías de gravedad residual isostática obtuvimos un modelo 3D del contraste de densidad en el subsuelo en la región del graben de Espino. El modelo explica adecuadamente la anomalía de Bouguer residual isostática y en general es consistente con los registros de densidad de diez pozos, de lo cual se concluye que el modelo describe de manera razonable la densidad en el subsuelo.

En el modelo se interpreta la superficie de cero contraste de densidad como el basamento económico; esta superficie describe una región cóncava hacia arriba, en cuyos bordes se aprecian paredes inclinadas de manera consistente con las fallas que limitan el graben. Se interpreta que el graben acomoda ~3500 m de rocas sedimentarias, formadas en el Jurásico durante el episodio de extensión cortical, sobre las que yacen ~3000 m de estratos sedimentarios sub horizontales generados con posterioridad al episodio de extensión cortical.

 

APENDICE A: LA ANOMALÍA REGIONAL ISOSTÁTICA

La anomalía de Bouguer residual isostática se obtiene al sustraer de la anomalía de Bouguer una anomalía de gravedad regional causada por la raíz cortical calculada con un modelo de compensación isostática. En este apéndice mostramos el cálculo del modelo de compensación isostático regional utilizado y la anomalía de gravedad regional que se deriva.

Los diferentes tipos de modelos de compensación isostática se clasifican como modelos de compensación local o regional. Para cualesquiera de estos dos modelos, la masa compensadora localizada por debajo de cadenas montañosas se conoce como la raíz, en tanto que por debajo de las cuencas oceánicas se conoce como la anti raíz. En ambos casos suponemos que la geometría de la raíz cortical se aproxima al límite corteza-manto.

En los modelos locales, la compensación isostática se puede obtener por cambios en el grosor de un estrato cortical cuya densidad es constante (Airy, 1855). En los modelos de compensación del tipo regional (Venning-Meinesz, 1948, en Watts, 2001) la masa compensadora se distribuye lateralmente dependiendo de la rigidez flexural de una placa elástica delgada, continua y homogénea que representa a la corteza.

En el modelo tipo Airy la flexión debida a una carga superficial está dada por

en donde h(x,y) representa la topografía (h>0) o batimetría (h<0) y w(x,y) representa la flexión producida, r_m es la densidad del manto, r₀ es la densidad de la topografía y r_f representa la densidad del fluido que la topografía desplaza: aire en regiones continentales y agua en regiones oceánicas. Sumando a la flexión el grosor cortical al nivel del mar (Z) obtenemos la profundidad al límite corteza-manto

El modelo de compensación regional toma en cuenta la rigidez flexural (D) de la placa elástica:

caracterizada por su módulo de Young (E), su relación de Poisson (v) y grosor elástico efectivo (T_e).

Para estimar la geometría y profundidad al límite corteza-manto en el norte de Suramérica (Figura A1), hemos utilizado un modelo de compensación isostática de tipo regional, para obtener la flexión de una placa elástica delgada, continua y homogénea debido a una carga distribuida en su superficie, mediante la solución derivada por Banks et al. (1977) en el dominio del número de onda:

en donde H(k_x, k_y) representa la topografía (y/o batimetría) en el dominio del número de onda y W(k_x, k_y) representa la flexión producida; k_x = 2π / λ_x y k_y = 2π / λ,_y son números de onda en las direcciones x y y, donde λ_x y λ_y representan longitudes de onda. La flexión obtenida en el dominio del número de onda se transforma al dominio del espacio mediante la transformada inversa de Fourier y se obtiene una estimación de la profundidad al límite corteza-manto utilizando la Ecuación A2.

La topografía y batimetría del norte de Suramérica se muestra en la Figura A1. Los principales rasgos fisiográficos en el occidente de esta región son los Andes de Colombia y Venezuela que alcanzan alturas de más de 5000 m sobre el nivel del mar (SNM), el cratón de Guayana cuya altura se incrementa de oriente a occidente desde ~500 a 2000 m s.n.m., y en el norte de Venezuela, la Cordillera de la Costa y la Serranía del Interior que alcanzan alturas de hasta ~1500 m s.n.m. La región de los llanos de Colombia y Venezuela está caracterizada por tierras bajas cuya elevación es en general menor a 200 m s.n.m.. Fuera de costa los rasgos más sobresalientes se encuentran en la región del Mar Caribe y en el Océano Atlántico, donde la batimetría alcanza profundidades abisales de más de 5000 m.

En la Tabla A1 se enlistan las densidades y parámetros elásticos utilizados en el cálculo de la profundidad a la frontera corteza-manto mediante la Ecuación A4. Para acomodar en este cálculo las variaciones laterales de densidad en regiones continentales, se supuso que la densidad de los sedimentos en las tierras bajas de la región de los llanos de Colombia y Venezuela, con alturas (z) menores a 500 m s.n.m. es de 2000 kg/m³ y que para el resto de las tierras elevadas la densidad de las rocas es de 2670 kg/m³. Utilizando estas densidades se calculó una pseudotopografía (h), cuya densidad es de 2670 kg/m³, utilizando la siguiente relación:

En la Figura A2 se muestra la pseudotopografía y batimetría seleccionada para este trabajo.

Tomando en cuenta las inferencias de la profundidad al límite corteza-manto en Venezuela reportadas por Schmitz et al. (2002, 2008), para la región de la cuenca Oriental se asignó un espesor cortical al nivel del mar (Z) de 40 km. Se procedió a estimar la profundidad al límite corteza-manto para la región, utilizando diversos valores de grosor elástico efectivo en la Ecuación A3; los resultados obtenidos utilizando T_e = 5, 10 y 15 km se muestran en la Figura A3.

En el norte de Suramérica existen inferencias recientes del grosor elástico efectivo: para los Andes de Colombia, Galán y Casallas (2010) infieren T_e = 20 km mediante la admitancia entre la anomalía de aire libre y la topografía; Arnaiz-Rodríguez et al. (2011) infieren T_e = 16 km para el occidente de Venezuela, incluyendo los Andes de Merida, la región del Lago Maracaibo y la cuenca Barinas-Apure.

Utilizando el resultado obtenido con un grosor elástico efectivo T_e = 15 km, debido a su parecido con la profundidad a la base de la corteza inferida por Schmitz et al. (2002, 2008) se calculó la anomalía regional isostática mostrada en la Figura A4, utilizando un contraste de densidad Δρ = 350 kg/m³ (positivo para el efecto de la anti raíz en regiones oceánicas y negativo para la raíz en regiones continentales), la cual se restó de la anomalía de Bouguer para obtener la anomalía de Bouguer residual isostática mostrada en la Figura 2.

 

AGRADECIMIENTOS

Deseamos agradecer la ayuda generosa de Arturo Martín Barajas, por revisar la primer versión del presente manuscrito, y de los revisores asignados por la RMCG: René Chávez Segura y Claudia B. Prezzi: sus valiosos comentarios y sugerencias han ayudado a mejorar substancialmente este trabajo.

 

REFERENCIAS

Airy, G.B., 1855, On the computations of the effect of the attraction of mountain-masses, as disturbing the apparent astronomical latitude of stations of geodetic surveys: Philosophical Transactions of the Royal Society of London, 145, 101-104.         [ Links ]

Andersen, O.B., Knudsen, P., Berry, P., 2010, The DNSC08GRA global marine gravity field from double re-tracked satellite altimetry: Journal of Geodesy, 84, 3, doi: 10.1007/s00190-009-0355-9.         [ Links ]

Arnaiz-Rodríguez, M.S., Rodríguez-Millán, I., Audemard, F., 2011, Análisis gravimétrico y flexural del occidente de Venezuela: Revista Mexicana de Ciencias Geológicas, 28 (3) 420-438.         [ Links ]

Banks, R.J., Parker, R.L., Huestis, S.P., 1977, Isostatic compensation on a continental scale: local versus regional mechanisms: Geophysical Journal of the Royal Astronomical Society, 51 (2), 431,452.         [ Links ]

Bezada, M., Schmitz, M., Jácome, M.I., Rodríguez, J., Audemard, F., Izarra, C., and the BOLIVAR Active Seismic Working Group, 2008, Crustal structure in the Falcón Basin area, northwestern Venezuela, from seismic and gravimetric evidence: Journal of Geodynamics, 45 (4-5), 191-200.         [ Links ]

Burke, K., 1977, Aulacogenes and continental breakup: Annual Reviewof Earth and Planetary Sciences, 5, 371-396.         [ Links ]

Fenglin, N., Bravo, T., Pavlis, G., Vernon, F., Rendón, H., Bezada, M., Levander, A., 2007, Receiver function study of the crustal structure of the southeastern Caribbean plate boundary and Venezuela: Journal of Geophysical Research, 112, B11308, doi:10.1029/2006JB004802.         [ Links ]

Galán, R.A., Casallas, I.F., 2010, Determination of effective elastic thickness of the Colombian Andes using satellite derived gravity data: Earth Sciences Research Journal, 14 (1), 7-16.         [ Links ]

Hinze, W.J., 2003, Bouguer reduction density, why 2.67?: Geophysics, 68, (5), 1559-1560.         [ Links ]

Jackson, D.D., 1979, The use of a priori data to resolve non uniqueness in linear inversion: Geophysical Journal of the Royal Astronomical Society, 57 (1), 137-157.         [ Links ]

MARAVEN, 1983, Área Zuata, en Evaluación Exploratoria de la Faja petrolífera del Orinoco, v. VII: Petróleos de Venezuela, S.A., Informe técnico MARAVEN FP-0206.         [ Links ]

Moticska, P., 1985, Volcanismo Mesozóico en el subsuelo de la Faja Petrolífera del Orinoco, Estado Guárico, Venezuela en VI Congreso Geológico Venezolano, v. 3, 1929-1943.         [ Links ]

PDVSA-Intevep, 2012, Código Estratigráfico de las Cuencas Petroleras de Venezuela: <www.pdv.com/lexico/lexicoh.htm>         [ Links ].

Salazar-Tomey, B.A., 2006, Evolución estructural e implicaciones tectónicas del graben de Espino: Universidad Simón Bolívar, Decanato de Estudios de Posgrado, Tesis de maestría, 197 pp.         [ Links ]

Schmitz, M., Chalbaud, D., Castillo, J., Izarra, C., 2002, The crustal structure of the Guayana Shield, Venezuela, from seismic refraction and gravity data: Tectonophysics, 345(1), 103-118.         [ Links ]

Schmitz, M., Ávila, J., Bezada, M., Vieira, E., Yáñez, M., Levander, A., Zelt, C. A., Jácome, M.I., Magnani, M.B., and the BOLIVAR active seismic working group, 2008, Crustal variations in Venezuela from deep seismic observations: Tectonophysics, 459(1-4), 14-26.         [ Links ]

Tarantola, A., 2005, Inverse Problem Theory and Methods for Model Parameter Estimation: Philadelphia, Society for Industrial and Applied Mathematics, 342 pp.         [ Links ]

Technical University of Denmark, 2012, Global gravity field model: <http://www.space.dtu.dk/>         [ Links ]

Venning-Meinesz, F.A., 1948, Gravity expeditions at sea, 1923-1938, v. IV, Complete results with isostatic reduction interpretation of the results: Delf, The Netherlands, Geodetic Commission, 24 pp.         [ Links ]

Watts, A.B., 2001, Isostasy and flexure of the lithosphere: Cambridge University Press, 458 pp.         [ Links ]