Zonificación Del Gradiente Geotérmico en La Cuenca Oriente de Ecuador A Partir de Temperatura de Fondo de Pozos

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 21

vol. 45, n.

° 1, enero-abril de 2023
ISSN impreso: 0120-0283 ISSN en línea: 2145-8553

Zonificación del gradiente geotérmico en la cuenca


Oriente de Ecuador a partir de temperatura de
fondo de pozos
Elena Angulo-Romero1* ; Bernardo Beate1 ; Christian Romero-Cóndor2
1
Ingeniería en Geología, Escuela Politécnica Nacional, Quito, Ecuador. (*) [email protected];
[email protected]
2
Instituto de Investigación Geológico y Energético, Quito, Ecuador. [email protected]

Resumen
La cuenca Oriente (Cretácico-Cenozoico tardío) es una de las más complejas y atractivas cuencas
sedimentarias de antepaís de los Andes septentrionales. Está constituida por diversos ciclos tectónicos
y sedimentarios, desarrollada sobre un basamento Precámbrico. Con el fin de evaluar el potencial
geotérmico teórico y los usos potenciales del calor en la cuenca, datos de temperaturas de fondo de
pozos (BHT) de 1683 pozos petroleros profundos, distribuidos en 202 campos individuales, han sido
ploteados en un mapa georreferenciado, cubriendo un área de 57.400 km2. El tratamiento estadístico
de los datos para el cálculo del gradiente geotérmico incluyó correcciones de las BHT, de acuerdo con
el método propuesto por Willet y Chapman. A partir de los datos trazados, se hace evidente que a lo
largo del Corredor Sacha-Shushufindi, muestra valores bajos de gradiente (7,38 a 26,1ºC/km). Por el
contrario, los valores altos de gradiente (26,1 a 69,01ºC/km) están presentes en el Sistema Subandino
y en el Sistema Capirón-Tiputini. La distribución de anomalías de temperatura parece estar relacionada
con el control estructural NNE-SSO de la cuenca, cuyo patrón de permeabilidad es bien conocido.
La extracción de calor para usos comerciales se ve altamente facilitada por la existencia de pozos
profundos que penetran en zonas de alto gradiente, lo que ahorra los costos generalmente prohibitivos
de la perforación profunda.

Palabras clave: Recurso geotérmico; Exploración y utilización de energía geotérmica; Energía


geotérmica en pozos petroleros.

Zoning of the geothermal gradient in the Oriente Basin of Ecuador from


Bottom Hole Temperatures
Abstract
The Oriente Basin (Late Cretaceous-Cenozoic) is one of the most complex and attractive foreland
sedimentary basins of the northern Andes. It is made up of various tectonic and sedimentary cycles
developed on a Precambrian basement. To evaluate the theoretical geothermal potential and the potential
uses of heat in the basin, data on Bottom Hole Temperatures (BHT) of 1683 deep oil wells, distributed in
202 individual fields, have been plotted on a geo-referenced map, covering an area of 57.400 km2. The
statistical treatment of the data for the calculation of the geothermal gradient included corrections of
the BHT according to the method proposed by Willet and Chapman. From the plotted data, it becomes
clear that the Sacha-Shushufindi Corridor shows low gradient values (7.38 to 26.1ºC/ km). Instead,
high gradient values (26.1 to 69.01ºC/km) are present in the Sub-Andean System and Capirón-Tiputini
System. The distribution of temperature anomalies seems to be related to the NNE-SSO structural
control of the basin, whose permeability pattern is well known. The extraction of heat for commercial
uses is highly facilitated by the existence of deep wells that penetrate high-gradient zones, saving the
generally prohibitive costs of deep drilling.

Keywords: Geothermal resource; Exploration and utilization of geothermal energy; Geothermal


energy in oil wells.

Forma de citar: Angulo-Romero, E.; Beate, B.; Romero-Cóndor, C. (2023). Zonificación del gradiente
geotérmico en la cuenca Oriente de Ecuador a partir de temperatura de fondo de pozos. Boletín de Geología,
45(1), 119-139. https://doi.org/10.18273/revbol.v45n1-2023006

119
Zonificación del gradiente geotérmico en la cuenca Oriente de Ecuador a partir de temperatura de fondo de pozos

Introducción niveles someros de la corteza (arco volcánico y punto


caliente activos, y una cuenca sedimentaria profunda),
La energía geotérmica es el calor contenido en el interior lo cual implica la existencia de importantes recursos
de la Tierra, disponible para su extracción y explotación. geotérmicos en su territorio con anomalías de calor
Los recursos geotérmicos se encuentran en todo el considerables, para el aprovechamiento de energía
mundo, pero los sistemas geotérmicos explotables se geotérmica. La exploración de recursos geotérmicos
encuentran principalmente en regiones con gradientes en el país se inició en 1979, cuando se descubrieron
geotérmicos normales o anómalos. Aunque el potencial y caracterizaron varios prospectos (Chachimbiro,
geotérmico está altamente concentrado en regiones Chalpatán, Chacana-Jamanco, Chalupas, Guapán,
volcánicas, también se puede encontrar como agua Chacana-Cachiyacu, Tufiño, Chimborazo, Chacana-
subterránea caliente en formaciones sedimentarias y Oyacachi, Baños de Cuenca y Alcedo), definidos como
con flujo de calor lo suficientemente alto como para de alta y media temperatura, y que presentaban un
ser utilizado (Saemundsson et al., 2009). mayor potencial al concentrarse en regiones volcánicas
(Beate y Urquizo, 2015). En algunas cuencas se estudió
Las cuencas sedimentarias y los recursos geotérmicos el gradiente geotérmico, el antearco (Manabí, Golfo de
están estrechamente relacionados, ya que están Guayaquil y la Provincia Geológica Progreso-Tumbes),
directamente influenciados por el ajuste tectónico. que presentan un gradiente geotérmico promedio bajo
Estos sistemas ocurren en capas sedimentarias (17 a 24°C/km), debido a la presencia de un basamento
permeables a gran profundidad (>30°C/km), debido oceánico frío de edad Cretácico medio (Barba, 2017).
a las fuentes de calor radiogénico en la corteza En el caso de la cuenca Oriente, se reportan gradientes
profunda o al levantamiento tectónico en la región geotérmicos promedios de 22°C/km (Burgos et al.,
(Saemundsson et al., 2009). Estos sistemas son de 2004), 21,4°C/km para un intervalo de 1,5 a 3,5 km de
naturaleza conductiva y albergan sistemas de acuíferos profundidad (Hamza et al., 2005) y 23,9°C/km (Barba
principales donde se puede producir y utilizar el agua et al., 2021).
termal, donde las fracturas y las fallas desempeñan un
papel importante (Moeck, 2014). Durante los últimos 50 años, la exploración petrolera
en la cuenca Oriente ha dado como resultado
Actualmente, la energía geotérmica almacenada en los la acumulación de miles de registros de pozos
yacimientos de hidrocarburos también presenta un gran (exploración, producción y reinyección) en los que
potencial, no solo porque existe energía geotérmica se han medido datos de Temperatura de Fondo de
masiva en los yacimientos de petróleo, sino también Pozos (BHT). Esta investigación evaluó más de 2500
porque los campos petroleros tienen enormes ventajas registros de BHT de la cuenca Oriente, almacenados
para desarrollar esa energía (Wang et al., 2018). en la Biblioteca de Información Petrolera del Ecuador
(BIPE), de los cuales 1683 datos fueron utilizados,
Ante la disminución de las reservas, el aumento del estos están distribuidos en 202 campos petrolíferos
costo de operación, los precios volátiles del petróleo individuales, que abarcan 57.400 km2, es decir, 2/3 de
y la tendencia de la energía verde, las compañías de la cuenca (Figura 1). El tratamiento de los datos de
petróleo y gas comenzaron a explorar y utilizar la temperaturas para el cálculo del gradiente geotérmico
energía geotérmica de los yacimientos de petróleo, aparente incluyó los valores de temperaturas media
buscando soluciones para reducir el costo de operación, anual superficial, así como las correcciones de BHT,
extender la vida económica de los campos antiguos y de acuerdo con el método dispuesto por Willet y
lograr beneficios medioambientales y sociales (Wang Chapman (1987).
et al., 2016). Como un tipo de recurso que coexiste con
hidrocarburos en las cuencas sedimentarias, el recurso El análisis de este conjunto de datos tiene como
geotérmico de los campos petroleros se caracteriza en objetivo estudiar la distribución de temperaturas en la
función de su temperatura, la cual es intermedia a baja, cuenca Oriente, de tal manera que se pueda determinar
dado que las temperaturas de los fluidos producidos el gradiente geotérmico y establecer su zonificación,
oscilan entre 65 y 150°C (Liu et al., 2018). además de caracterizar su potencial para un futuro
aprovechamiento.
Ecuador posee una situación geodinámica favorable
para la presencia de anomalías térmicas positivas en

120 Boletín de Geología - vol. 45, n.° 1


Elena Angulo-Romero; Bernardo Beate; Christian Romero-Cóndor

Figura 1. Ubicación del área de estudio: se incluyen los datos de temperaturas de fondo de pozos disponibles, así como los
principales campos petrolíferos y los dominios tectónicos, modificado de Baby et al. (2013).

Área de estudio Baby et al. (2013) proponen tres dominios tectónicos:


Subandino, Sacha-Shushufindi y Capirón-Tiputini.
La cuenca Oriente forma parte del sistema de cuencas
de trasarco de la margen noroccidental sudamericana, es La estratigrafía de la cuenca Oriente (Figura 2), se
contigua a las cuencas Putumayo en el sur de Colombia compone de diversos ciclos tectónicos y sedimentarios
y Santiago-Marañón en el norte de Perú (Higley, 2001). desarrollados sobre un basamento Precámbrico (Díaz
La cuenca es asimétrica, se superpone gradualmente al et al., 2004). En el basamento se reportaron granulitas
escudo guyanés al este y está limitada por el cinturón en el pozo Tapi, gneis micáceos en el pozo Tiputini
plegado y sobrecorrido Subandino en el oeste (Ruiz 1 y esquistos micáceos en el pozo Dureno 1 (Vera,
et al., 2007). Respecto a la estructura de la cuenca, 2016). En el límite sur de la cuenca Putumayo, en el

Boletín de Geología - vol. 45, n.° 1 121


Zonificación del gradiente geotérmico en la cuenca Oriente de Ecuador a partir de temperatura de fondo de pozos

pozo Caima 3, edades U/Pb de circones obtenidos La Formación Santiago es una sucesión marina del
de granulitas atribuidas al basamento indican Jurásico, expuesta a lo largo del río Santiago, en el
edades 1,4 Ga y 0,95-0,98 Ga (Ibañez-Mejia et al., sur de la Cordillera Cutucú, al sureste de Ecuador
2011). Secuencias monótonas de cuarcitas y filitas (Goldschmid, 1941). En la base de la Formación
intensamente deformadas de la Formación Pumbuiza Santiago se observan monótonas secuencias de calizas
fueron asignadas al Devónico con base en asociación negras a grises intercaladas con areniscas y lutitas con
de acritarcos (Ordóñez et al., 2006). La Formación materia orgánica, mientras que al tope se reportaron
Pumbuiza es sobreyacida en discordancia por calizas, intercalaciones de rocas volcánicas y volcanoclásticas
lutitas fosilíferas, areniscas calcáreas y esporádicos (Tschopp, 1953; Romeuf et al., 1995). Gaibor et al.
niveles volcánicos de la Formación Macuma (Tschopp, (2008) discriminaron a la Formación Santiago en tres
1953). La abundante fauna reportada en la Formación miembros: Río Santiago, Yuquianza y Patuca, con base
Macuma indica que estas rocas posiblemente fueron en análisis palinológico y en el registro de estructuras
depositadas en ambientes marino-costeros en el sedimentarias.
Carbonífero (Christophoul, 1998; Ordóñez et al., 2006).

Figura 2. Relaciones litoestratigráficas de las secuencias sedimentarias de la cuenca Oriente, modificado de Gutiérrez et al.
(2019) y Vallejo et al. (2021).

122 Boletín de Geología - vol. 45, n.° 1


Elena Angulo-Romero; Bernardo Beate; Christian Romero-Cóndor

La amplia gama de amonites y bivalvos de la Formación et al., 2000). Secuencias de conglomerados, brechas,
Santiago, recolectados en la Cordillera Cutucú, indican litoareniscas y limolitas vinculadas a secuencias de
una edad Hettangiano a Pliensbachiano (Aspden y relleno de valle aluvial peneplanizan el relieve sobre
Ivimey-Cook, 1992; Gaibor et al., 2008; Romero et al., el cual se depositó la Formación Hollín (White et al.,
2019). Vallejo et al. (2021) sugieren que la base de la 1995; Vallejo et al., 2021). La Formación Hollín puede
Formación Santiago fue depositada en una plataforma ser dividida en dos miembros: el miembro Inferior y el
carbonatada somera en una cuenca de trasarco. Mientras miembro Superior (Canfield et al., 1982). El miembro
que, hacia el tope de la Formación Santiago series Inferior se compone de secuencias estrato-decrecientes
de deltas progradaron hacia el este y se depositaron de conglomerados clasto y matriz-soportados de
de manera coetánea al arco volcánico Piuntza, que cuarzo, cuarzoareniscas con estratificación cruzada,
representaría la expresión superficial del magmatismo gradación normal y arenisca finas masivas (Macellari,
del Jurásico inferior. La Formación Sacha, definida al 1998; Romero et al., 2019; Vallejo et al., 2021).
norte de la zona Subandina (Díaz et al., 2004), podría El miembro Inferior se vincula a la sedimentación
representar un equivalente clástico no marino de la continental de cinturones de ríos trenzados (White
Formación Santiago (Vallejo et al., 2021). et al., 1995; Shanmugam et al., 2000; Romero et al.,
2019; Vallejo et al., 2021). En onlap sobre el miembro
La Formación Chapiza cubre en discordancia angular Inferior se deposita el miembro Superior, que contiene
a las formaciones Santiago y Sacha, agrupa a limolitas múltiples intercalaciones rítmicas de areniscas y
multicolores con yeso, litoareniscas, areniscas arcillolitas con laminación lenticular, ondulitica,
volcánicas, conglomerados volcánicos depositados flaser, crinkle, wispy, superficies de reactivación y
en planicies de inundación, abanicos aluviales con ripples bidireccionales que expresan ambiente costero
canales fluviales trenzados y depósitos menores de con influencia mareal (White et al., 1995; Jaillard et
debris flows de edad Jurásico superior a Cretácico al., 1997; Romero et al., 2019; Vallejo et al., 2021).
inferior (Tschopp, 1953; Romero et al., 2021; Vallejo De esta manera, la litoestratigrafía de la Formación
et al., 2021). La Formación Chapiza fue depositada Hollín registra la transición desde ambientes fluviales
de manera sincrónica a un periodo de extensión a costeros influenciados por mareas, ocasionado por la
regional en una cuenca de trasarco, donde se formaron subida del nivel del mar en el Albiano y documentada
estructuras de semigraben que fueron parcialmente ampliamente en la margen noroccidental sudamericana
rellenadas por depósitos clásticos de la Formación (Bulot et al., 2005).
Chapiza (Vallejo et al., 2021). Al oeste de la cuenca,
la Formación Misahuallí incluye andesitas, andesitas La Formación Napo cubre en concordancia a la
basálticas, riolitas y depósitos piroclásticos vinculados Formación Hollín (Figura 2). La Formación Napo
a un arco volcánico continental de afinidad calco- puede ser analizada en función de varios ciclos de
alcalina (Romeuf et al., 1995), que es considerado sedimentación controlados por variaciones eustáticas
la expresión extrusiva del magmatismo del Jurásico en el intervalo Albiano a Maastrichtiano (Brookfield et
medio a superior (Litherland et al., 1994). al., 2009; Vallejo et al., 2021). La Formación Napo tiene
un espesor que varía entre 240 a 900 m (Dashwood y
Un evento compresivo en el Cretácico inferior afectó Abbotts, 1990; Mello et al., 1995). La Formación Napo
parte del Ecuador, ocasionando rápida exhumación y incluye a doce miembros informales (Jaillard et al.,
enfriamiento de rocas precretácicas en la Cordillera 1997; Vallejo et al., 2021).
Real y Zona Subandina (Litherland et al., 1994;
Ruiz et al., 2007). En la cuenca Oriente, este evento La secuencia del Albiano superior es representada
corresponde a un hiato sedimentario y una discordancia por lutitas negras laminadas agrupadas en el miembro
angular entre las formaciones Hollín y Chapiza (Vallejo Lutitas Napo Basal (White et al., 1995; Vallejo et al.,
et al., 2021). 2021). Un intervalo transgresivo es reportado en la
secuencia del Albiano superior al Cenomaniano inferior
En el Cretácico tardío, la sedimentación en la cuenca (Jaillard et al., 1997; Vallejo et al., 2021), el cual
Oriente fue controlada por variaciones del nivel incluye los miembros Arenisca T, Caliza B y Lutitas
eustático y la generación de espacio de acomodación Napo Medio (Figura 2). El miembro Areniscas T se
(Brookfield et al., 2009). La Formación Hollín se vincula a depósitos de estuario que incluyen ambientes
depositó en discordancia angular sobre un basamento intermareales a submareales. El miembro Caliza B
heterogéneo intensamente deformado (Shanmugam fue depositado sobre la Arenisca T y es una monótona

Boletín de Geología - vol. 45, n.° 1 123


Zonificación del gradiente geotérmico en la cuenca Oriente de Ecuador a partir de temperatura de fondo de pozos

secuencia de lutitas negras. El miembro Lutitas Napo Caribe contra la margen noroccidental sudamericana
Medio incluye lutitas negras marinas intercaladas con registrada en el intervalo 70 a 73 Ma (Luzieux et al.,
estratos de lodolitas calcáreas (Vallejo et al., 2021). 2006; Vallejo et al., 2021).

En el intervalo Cenomaniano superior a Turoniano En el Cenozoico, la sedimentación en la cuenca Oriente


superior fueron incluidos los miembros Arenisca U y está controlada por periodos de exhumación y erosión
Caliza A. El miembro Arenisca U se vincula a facies de los Andes. Sobre la Formación Tena en discordancia
de estuario que evolucionaron a plataforma marina angular se deposita la Formación Tiyuyacu, que consiste
(Shanmugam et al., 2000). El miembro Caliza A en gruesos paquetes de conglomerados y areniscas
cubre al miembro Arenisca U (Jaillard et al., 1997). El fluviales del Eoceno inferior a medio (Christophoul et
miembro Caliza A representa una secuencia calcárea, al., 2002; Romero et al., 2021). En el centro norte de
que agrupa a lodolitas negras a grises y calizas tipo la cuenca, la Formación Chalcana agrupa a limolitas y
wackestone abundantemente fosilíferas (Vallejo et al., areniscas con abundante materia orgánica de carácter
2021). continental (Christophoul et al., 2002), mientras que
al norte de la cuenca la Formación Orteguaza incluye
La secuencia Turoniano a Santoniano incluye a los a intercalaciones de areniscas y limolitas costeras
miembros: Arenisca M2, Caliza M2, Caliza M1 y (Rivadeneira y Baby, 1999). Tanto las formaciones
Lutita Napo Superior. El miembro Arenisca M2 está Chalcana como Orteguaza son vinculadas al Eoceno
restringido al noreste de la cuenca Oriente, e incluye medio-Oligoceno temprano (Christophoul et al., 2002).
secuencias de areniscas glauconitas intercaladas El relleno Neógeno comprende depósitos gruesos
con calizas y lutitas calcáreas (Jaillard et al., 1997). no marinos de las formaciones Arajuno y Chambira
Vallejo et al. (2021) afirman que en el resto de la lateralmente continuos con depósitos marinos y
cuenca el equivalente lateral del miembro Arenisca M2 lacustres poco profundos de la Formación Curaray al
corresponde a un intervalo de lutitas negras marinas este de la cuenca (Christophoul et al., 2002). Estas
entre los miembros Caliza A y Caliza M2. Secuencias secuencias están cubiertas por el mega abanico aluvial
de lutitas negras con calizas tipo wackestone con altos de Pastaza, que se extiende desde la depresión Pastaza
valores de TOC fueron agrupadas en el miembro Caliza hasta la cuenca Marañón (Bès de Berc et al., 2005).
M2 (Jaillard et al., 1997; Vallejo et al., 2002). Un
marcador estratigráfico regional dentro del miembro
Caliza M1 es el marcador L, que fue interpretado Metodología
como una superficie de máxima inundación del ciclo
La información utilizada fue obtenida de las
Napo y es utilizado como dato cronoestratigráfico para
tablas de datos suministradas por el Banco de
correlaciones de pozos (Brookfield et al., 2009; Vallejo
Información Petrolera de Ecuador, BIPE (Secretaría
et al., 2017). En el Campaniano, el miembro Arenisca
de Hidrocarburos). Las temperaturas de fondo de
M1 es vinculado a un ambiente deltaico debido a la
pozos (BHT) utilizadas para este análisis se tomaron
presencia de depósitos gravitacionales y tendencia
de los encabezados de los registros de pozos; se
estrato-creciente.
consideraron valores de TD, MD y TVD. La selección
se hizo teniendo en cuenta criterios de calidad de la
La Formación Tena cubre en discordancia angular a
información e incluyó la eliminación de registros de
la Formación Napo (Jaillard et al., 1997; Gutiérrez et
pozo sin información, eliminación de redundancia de
al., 2019). La Formación Tena fue depositada en el
datos, eliminación de registros inconsistentes. De una
intervalo Maastrichtiano a Paleoceno y representa a
base inicial de 2500 registros se utilizaron 1683 parejas
secuencias de sedimentos continentales que marcan el
de datos de temperatura-profundidad para este estudio.
final de la sedimentación marina dentro de la cuenca
Oriente, es decir, el comienzo de la sedimentación
El tratamiento de la información incluyó tres etapas:
continental en la cuenca. Este cambio probablemente
1) organización, selección de la información y
ocurrió en respuesta al acelerado levantamiento de los
conversión de unidades; 2) corrección de temperaturas
Andes y marcó el cambio a un sistema de cuenca de
de fondo de pozos y estimación de gradientes
antepaís (Martin-Gombojav y Winkler, 2008; Spikings
geotérmicos aparentes; y 3) aplicación de un método
et al., 2010; Gutiérrez et al., 2019; Vallejo et al., 2021).
de interpolación y generación del mapa.
Este cambio coincide con la colisión oblicua del Plateu

124 Boletín de Geología - vol. 45, n.° 1


Elena Angulo-Romero; Bernardo Beate; Christian Romero-Cóndor

Corrección de las BHT y cálculo del gradiente Willet y Chapman (1987), proponen una corrección
geotérmico aparente empírica aplicada a mediciones de pozos para fines
Durante la excavación de un pozo petrolero, se supone geotérmicos, comparando las temperaturas de fondo
que la temperatura máxima está en el fondo del pozo, con las temperaturas medidas en equilibrio. El ajuste de
por lo que da una medida directa de la temperatura la diferencia entre estas dos temperaturas, por mínimos
del fluido de perforación a la profundidad actual del cuadrados, dio lugar a una ecuación polinomial de
pozo, o una “temperatura del fondo del pozo” (BHT). tercer orden mediante la cual se obtiene la corrección
La relación entre BHT y la “temperatura de la roca de la temperatura de fondo de pozo como una función
virgen” (VRT) es compleja y hay casi tantas técnicas de la profundidad. La ecuación para la corrección se
para corregir las BHT como publicaciones sobre el describe a continuación:
tema (Beardsmore y Cull, 2001).

V RT = BHT + Tc (1)

Tc = 6, 93z − 1, 67z 2 + 0, 101z 3 + 0, 0026z 4 (2)

Donde Tc es la corrección de temperatura en °C. El gradiente es aparente debido a que no tiene en cuenta
polinomio de tercer orden empleado para la corrección la variación vertical asociada a los cambios litológicos
de BHT corresponde a la ecuación 2, donde z es la y a las propiedades de transferencia de calor de la roca
profundidad en km (Beardsmore y Cull, 2001). (Alfaro et al., 2015). La validez de este método se
fundamenta en que las mediciones de temperatura son
Se debe considerar, que cualquier corrección empírica las menos perturbadas durante el período de medición
produce un valor aproximado de la magnitud de la después de terminar las actividades de perforación
corrección promedio del BHT y tiene en cuenta que las (Deming, 1989).
temperaturas reales de formación son mayores a las de
fondo de pozo (alrededor de 10 a 15°C en promedio); Suponiendo una relación lineal de temperatura con
que la cantidad de corrección incrementa con la la profundidad, el gradiente geotérmico se calcula
profundidad pero no de forma lineal y que puede aún dividiendo la diferencia entre la temperatura corregida
disminuir a grandes profundidades; que la cantidad de y la temperatura media anual superficial, por la
corrección disminuye con el tiempo de cierre y que profundidad, de acuerdo con la ecuación empírica:
en todo caso, la aplicación de cualquier corrección
a los datos de BHT incorpora un error desconocido BHT − T SU P
(Deming, 1989). GGA = (3)
ZBHT
El gradiente geotérmico es la variación de la
temperatura con la profundidad. En la corteza, Donde BHT es la temperatura del pozo corregida en
esta variación es más irregular que en el manto y el °C a una profundidad z en km, TSUP es la temperatura
núcleo, y registra grandes cambios laterales, lo cual media anual superficial en °C.
depende en gran medida del contenido de isótopos
radiactivos en las rocas, de la tectónica y del régimen Análisis y representación de los datos
hidrológico (Clauser, 2009). En la corteza, el gradiente Para la generación del mapa de gradientes geotérmicos
de temperatura es típicamente 30°C/km, pero puede se utilizó el método de interpolación de mínima
ser tan alto como 150°C/km en áreas geotérmicas curvatura con tamaño de celdas de 1673 m (tanto
calientes (ESMAP, 2012). en X como en Y) y distancia de blanqueo de 5 km,
en el software Geosoft-Target de ArcGIS. Todas
El cálculo del gradiente geotérmico a partir de las representaciones se generaron a partir del mejor
temperaturas de fondo de pozo, corresponde al gradiente variograma teórico que se ajustó al variograma
estimado de la medición de la temperatura del fondo experimental de los datos. El mapa fue editado en el
a la profundidad máxima respecto a la temperatura software ArcMap.
medida en superficie (Carvalho y Vacquier, 1977). Este

Boletín de Geología - vol. 45, n.° 1 125


Zonificación del gradiente geotérmico en la cuenca Oriente de Ecuador a partir de temperatura de fondo de pozos

Resultados En todos los campos de la cuenca Oriente se registraron


mediciones de temperatura a diferentes profundidades,
Los gradientes geotérmicos fueron calculados a partir los gradientes registran variación con la profundidad
de 1683 pares de datos de temperatura y profundidad en en un mismo campo. Esta variación presumiblemente
pozos petroleros. Las temperaturas registradas varían está relacionada con factores como cambios en
entre 29,53 y 139,35°C a una profundidad entre 0,27 litología y textura, en las porosidades de las rocas, en
y 4,9 km. La dispersión de las temperaturas corregidas la conductividad térmica, en los fluidos circulantes,
en función de la profundidad para la cuenca Oriente se entre otros factores, los cuales no fueron evaluados y
ilustra en la Figura 3; se observa que mantienen una no forman parte del alcance de este trabajo.
tendencia directamente proporcional entre el aumento
de la profundidad con el aumento de la temperatura. Se En la cuenca Oriente se han podido identificar dos
evidencia una clara acumulación de las temperaturas zonas con gradientes geotérmicos anómalos, mayores
intermedias a altas en profundidades entre los 2,2 km a 40°C/km, y por lo general se encuentran asociados
y los 3,8 km. a profundidades someras (<2 km de profundidad
de los reservorios). Dichas zonas anómalas (Figura
Los ajustes de mínimos cuadrados al conjunto de datos 4) se encuentran en los campos individuales de: 1)
de BHT aplicando la corrección propuesta por Willet Pungarayacu-San Carlos, cuyas temperaturas de fondo
y Chapman (1987), indican un gradiente geotérmico oscilan entre los 30 y 110°C a profundidades entre
de 26,1°C/km para un intervalo de profundidad entre los 0,27 y 4,07 km, y 2) Cowi cuya única temperatura
0,2 y 4,9 km. Esto es ligeramente mayor al calculado de fondo es 130°C a 0,385  km de profundidad.
por Burgos et al. (2004) y Hamza et al. (2005). A Adicionalmente, hay dos zonas con anomalías entre
partir de los datos trazados, se hace evidente que, la 30 y 40°C/km asociados a profundidades <2,5 km,
parte central de la cuenca, a lo largo del Corredor cuyos campos individuales son: 1) Bermejo-El Rayo,
Sacha-Shushufindi, muestra los valores más bajos de cuyas temperaturas de fondo oscilan entre los 55 y
gradiente, que van desde 7,38 a 26,1ºC/km. Por el 110°C a profundidades entre los 1,3 y 1,7 km; y 2) ITT
contrario, los valores más altos de gradiente, que van cuyas temperaturas de fondo oscilan entre 70 y 95°C a
desde 26,1 a 69,01ºC/km, están presentes en la parte profundidades entre los 1,6 y 2,3 km.
occidental (Levantamiento Napo y Subandino) y en
la parte oriental (Corredor Capirón-Tiputini), como se
ilustra en la Figura 4.

Figura 3. Variación de BHT corregidas con la ecuación de Willet y Chapman con la profundidad. La temperatura máxima de
139,25°C, se registra a 4,4 km y representa al pozo Cangaime 001.

126 Boletín de Geología - vol. 45, n.° 1


Elena Angulo-Romero; Bernardo Beate; Christian Romero-Cóndor

Figura 4. Mapa de gradientes geotérmicos de la cuenca Oriente. Se observa que las zonas anómalas con gradientes superiores
a los 40°C/km se encuentran en el Levantamiento Napo y el Corredor Capirón-Tiputini. Los colores cálidos indican anomalía
positiva, y los colores fríos, anomalías negativas, con un rango de variación desde 7,38°C/km a 69,01°C/km.

Anomalías Pungarayacu-San Carlos y Bermejo-El positiva creando condiciones favorables para la


Rayo migración de fluidos.
Estas zonas anómalas por lo general se encuentran
asociadas a profundidades someras (<2 km) y con Los campos Pungarayacu y Bermejo están siendo
gradientes >40°C/km. La parte occidental de la cuenca afectados por la última reactivación tectónica iniciada
Oriente se caracteriza por una tectónica profunda, que en el Plioceno (Rivadeneira y Baby, 2004). Debido
expone las secuencias pre-Cretácicas como resultado a esta reactivación tectónica, es posible que los
del levantamiento Napo-Galeras y la Cordillera fluidos más profundos estén alcanzando niveles más
del Cutucú. En esta parte se ubican las estructuras someros siguiendo caminos a lo largo de fracturas y
Bermejo-El Rayo y Pungarayacu-San Carlos. De fallas (Beate et al., 2020), lo que explicaría que estos
acuerdo con Rivadeneira y Baby (2004), el Campo campos presenten anomalías positivas de gradientes
Bermejo está asociado al sistema de transpresión geotérmicos. La presencia de estas anomalías positivas
dextral del frente de empuje Subandino; por su parte, de calor perceptibles está fuertemente influenciada
el Campo Pungarayacu se ubica en el periclinal sur del por la cercanía del basamento debido a procesos
gran anticlinal Napo-Galeras, el mismo que a finales tectónicos (levantamiento de las formaciones y erosión
del Cretácico empezó a actuar como una estructura de orógenos), además, los campos Pungarayacu-San

Boletín de Geología - vol. 45, n.° 1 127


Zonificación del gradiente geotérmico en la cuenca Oriente de Ecuador a partir de temperatura de fondo de pozos

Carlos y Bermejo-El Rayo se encuentran cerca de distribución de permeabilidad y porosidad en las rocas
los volcanes activos Sumaco (alcalino) y Reventador reservorios; la distribución de sistemas de fallas y
(calco-alcalino), por lo que pueden considerarse como fracturas principales en la cuenca y fuentes/tasas de
las fuentes de calor adicional que abastecen al sistema recarga de fluidos.
geotérmico.
Basamento y fuente de calor
El valor de gradiente geotérmico más alto obtenido en Para la cuenca Oriente se ha considerado que la
la cuenca Oriente corresponde al pozo Pungarayacu IP- fuente de calor podría estar asociada y depender de
5B, con un valor de 69,01°C/km a una profundidad de la geometría del basamento y, en algunas partes, del
328 m llegando a niveles de Hollín Basal e ingresando volcanismo reciente; por ejemplo, en el Levantamiento
a la Formación Chapiza, lo cual aparenta ser una Napo, la cercanía a los volcanes Sumano y Reventador,
anomalía local, ya que entrando al basamento más influiría en el gradiente geotérmico. Por otro lado, el
profundo es más frío, como se evidencia en el pozo espesor de las secuencias sedimentarias que cubren al
IP-17 el cual alcanzó una profundidad de 4,14 km y se basamento representa una profundidad adecuada (1,5-
obtuvo un gradiente de 21,53°C/km. 2,5 km) para que formaciones geológicas permeables,
capaces de contener y dejar circular fluidos, extraigan
Anomalías Cowi e ITT el calor del basamento. También hay que considerar
Estas dos zonas anómalas presentan gradientes que una de las mayores fuentes de calor terrestre es la
geotérmicos entre 30 y 40°C/km asociados a desintegración de isótopos radiactivos presentes en las
profundidades someras (<2,5 km). En la parte oriental rocas, ese calor varía con la composición química y
de la cuenca Oriente se han desarrollado dos grupos con la edad de las rocas (Moeck, 2014).
de fallas de deslizamiento de rumbo inverso de alto
ángulo (reverse strike-slip), respectivamente, con En el Ecuador se ha descrito al basamento de la
orientaciones NNE y NNO (Lin et al., 2015). La cuenca como rocas ígneas y metamórficas (Díaz et al.,
distribución del sistema de fallas generalmente se 2004); de acuerdo con Litherland et al. (1994), rocas
considera como resultado de la inversión tectónica del cratón Amazónico fueron datadas por Rb-Sr en
de fallas activas y la ocurrencia de un movimiento ~1,6 Ga en muestras del pozo Tapi. Litológicamente
de deslizamiento de rumbo inverso (Lin et al., 2015). se ha descrito como granulitas en el pozo Tapi, gneis
Esta distribución de fallas activas permite inferir que micáceos en el pozo Tiputini 1 y esquistos micáceos
el mecanismo de conducción de calor es a través del en el pozo Dureno 1 (Vera, 2016). En el pozo Caima
basamento cristalino. Como se observa en el mapa 3, localizado en la cuenca Putumayo, Ibañez-Mejia et
expuesto en la Figura 4, se registra un incremento en al. (2011) dataron rocas atribuidas al basamento en
el gradiente geotérmico hacia el oriente de la cuenca 1,4 Ga, 0,98 Ga y 0,95 Ga con el método de U-Pb.
con una tendencia suroriental. El corredor Capirón- De acuerdo con Gutiérrez et al. (2019), el cratón
Tiputini muestra una distancia vertical más cercana amazónico (3,2-2,6 Ga) comprende un núcleo de rocas
al basamento caliente, lo que pudiera significar la ígneas y metamórficas Arqueanas con cinturones de
presencia de estas anomalías positivas de calor y acreción adyacentes de tendencia NE.
puede estar asociado a concentraciones de elementos
radiactivos en el basamento. Si bien el calor se mueve mayormente vía convección
siempre y cuando haya zonas con alta permeabilidad,
Discusión generalmente la conducción es el mecanismo
dominante cuando el transporte de calor se estudia
Las cuencas sedimentarias poseen sistemas dentro de la corteza terrestre. En el caso de las
geotérmicos de intermedia a baja temperatura, el mapa cuencas sedimentarias el transporte de calor por
de gradientes geotérmicos define varias anomalías advección desempeña un papel importante, donde el
positivas localizadas en la cuenca. En ese contexto, a dominio de alta permeabilidad y las fallas enraizadas
continuación, se describen los componentes del sistema permiten la circulación profunda de aguas meteóricas
geotérmico hipotético en la cuenca Oriente. El modelo (Moeck, 2014). Estos sistemas geotérmicos rara vez
presentado es teórico considerando la información se asocian con reservorios hidrotermales, sino que
disponible y publicada. Las características del sistema son el resultado de sistemas de circulación profunda
geotérmico hipotético en la cuenca consideran la asociados con fallas complejas de escamas corticales
naturaleza y la profundidad de la fuente de calor; el mayores. Estas áreas típicamente contienen flujo de
mecanismo de transferencia de calor dominante; la calor bajo a moderado (Allen et al., 2006).

128 Boletín de Geología - vol. 45, n.° 1


Elena Angulo-Romero; Bernardo Beate; Christian Romero-Cóndor

Principales reservorios sedimentarios en la cuenca (Mancilla et al., 2008). En el centro norte de la cuenca
Oriente su espesor alcanza 150 m (Brookfield et al., 2009) y
Los principales reservorios en la cuenca Oriente son las posee amplia distribución en toda la cuenca acuñándose
areniscas de las formaciones cretácicas Hollín y Napo, al este (Vallejo et al., 2021).
los cuales proponemos para este modelo geotérmico
hipotético. Los valores de porosidad y permeabilidad para
la Formación Hollín se resumen en la Tabla 1.
Areniscas de la Formación Hollín: la Formación Romero (2018) afirma que en la Formación Hollín
Hollín es el principal reservorio de la cuenca, agrupa la permeabilidad es inversamente proporcional al
a múltiples sets de areniscas cuarzosas de origen contenido de láminas continuas de lodo, debido a que
fluvial (miembro Inferior) que evolucionan a finas estas actúan como un sello e impiden el paso de fluidos,
intercalaciones rítmicas de areniscas y arcillolitas mientras que, la porosidad es directamente proporcional
vinculadas a ambientes costeros con influencia mareal al tamaño de grano y grado de sorteo, disminuyendo en
que gradualmente son cubiertos por depósitos de las rocas con tamaño de grano más fino.
plataforma marina epicontinental (miembro Superior)

Tabla 1. Rango de porosidad y permeabilidad en la Formación Hollín.


Porosidad (%) Permeabilidad (mD) Referencia
12,1-20,4 66,2-350,67 Romero (2018)
Hollín, miembro
12-25 20-2000 Dashwood y Abbotts (1990)
Inferior
16 400-800 Canfield et al. (1982)
7-19,7 375,2-732,1 Romero (2018)
Hollín, miembro
20 658,61 Zhongzhen et al. (2017)
Superior
12-25 20-2000 Yinfu et al. (2010)

Arenisca de la Formación Napo: dentro de la lodolitas (U Superior), representando facies de estuario


Formación Napo se consideran reservorios las areniscas (Dashwood y Abbotts, 1990; Vallejo et al., 2021). El
“T”, “U” y “M1”. Estas secuencias progradan desde el espesor de las Areniscas U fue reportado en 70 m al este
este de la cuenca, su geometría indica que se acuñan de la cuenca (Estupiñán et al., 2006); U Inferior tiene
hacia el centro-oeste de la cuenca (Shanmugam et al., las mejores características petrofísicas, con valores de
2000). porosidad que varían entre 15 y 20%, mientras que para
U Superior el rango va de 10 a 18% (Rivadeneira y
Las Areniscas T (Albiano tardío a Cenomaniano Almeida, 2014).
temprano) agrupan areniscas de grano medio a fino (T
Inferior) e intercalaciones de areniscas glauconíticas La Arenisca M1 (Campaniano Medio) agrupa a
y lodolitas (T Superior) de ambiente intermareal a brechas y areniscas finas deltaicas (Vallejo et al.,
marino. El espesor de T Inferior al este de la cuenca se 2017) distribuidas al este de la cuenca, su espesor no
mantiene constante próximo a los 50 m (Smith, 1989). es homogéneo y varía entre 10 a 80 m, con valores de
Zhongzhen et al. (2017) indican que T Inferior posee porosidad de 18-32%, y una permeabilidad promedio
buena permeabilidad y alta porosidad. El espesor de de 1013,25 mD (Zhongzhen et al., 2017).
T Superior al centro-este de la cuenca alcanza los 40
m (Smith, 1989). Dashwood y Abbotts (1990) afirman Principales capas sellos de la cuenca Oriente
que T Superior contiene láminas de arcilla y alta
Dentro de la secuencia estratigráfica cretácica las
proporción de glauconita, estas características la hacen
capas sello son definidas como zonas impermeables
menos permeable. De acuerdo con Dashwood y Abbotts
sin porosidad que impiden la movilidad de fluidos;
(1990), las condiciones de porosidad y permeabilidad
según la calidad de las capas estas pueden conservar
típicos son 10-25% y 100-6000 mD para la Arenisca T.
características físicas del fluido como son la temperatura
(Sass y Götz, 2012). En el miembro Superior de la
Las Areniscas U (Cenomaniano tardío a Turoniano
Formación Hollín, estratos métricos de arcillolitas
temprano) presentan morfología de barras mareales,
ricas en materia orgánica vinculadas a planicies de
compuestas principalmente por areniscas calcáreas
inundación lodosas en ambientes intermareales y
(U Inferior) e intercalaciones rítmicas de areniscas y
Boletín de Geología - vol. 45, n.° 1 129
Zonificación del gradiente geotérmico en la cuenca Oriente de Ecuador a partir de temperatura de fondo de pozos

ambientes supra-mareales son impermeables y no constituyen un reservorio de agua salobre justificado


tienen porosidad; dependiendo de sus espesores estos por la presencia de glauconita y niveles conchíferos
pueden ser considerados capas sello (Smith, 1989; intercalados. Los registros de pozos y pruebas
Romero, 2018). Sin embargo, no son regionales y en testigos de perforación documentan un patrón
pueden variar de espesor a corta distancia. consistente de aumento de la salinidad de este a oeste
dentro de las areniscas individuales y verticalmente de
En la Formación Napo las capa sello, que van variando areniscas más profundas a menos profundas (Smith,
de espesor hacia el este de la cuenca, son: las Lutitas 1989). Además, el modelo de flujo hidrodinámico
Napo Inferior (Albiano medio), que agrupa a niveles descrito por Smith (1989), explica que el agua ingresa
de lutitas calcáreas fosilíferas (Jaillard et al., 1997), (zona de recarga) a la cuenca a través de las secciones
en el centro de la cuenca su espesor alcanza los 100 elevadas y expuestas hacia el oeste de la cuenca, y se
m. Lutitas Napo Medio (Cenomaniano), que incluye descarga a lo largo del borde poco profundo oriental de
a arcillolitas calcáreas y lutitas fosilíferas (Jaillard la cuenca (Figura 5).
et al., 1997), en el noroeste de la cuenca presentan
espesor próximo a los 50 m (Rivadeneira, 1986). Las mediciones de salinidad derivadas de algoritmos
Lutitas Napo Superior (Santoniano) que se compone en el estudio de Smith (1989) revelan que: 1) el agua de
de secuencias monótonas de lutitas negras calcáreas, formación de la arenisca Hollín en la base de la sección
en el centro-este de la cuenca su espesor supera los del Cretáceo es uniformemente fresca en gran parte de
90 m (Rivadeneira, 1986). A nivel regional, capas la cuenca central; 2) las areniscas de Hollín tienden a
de limolitas multicolores, con materia orgánica de engrosarse hacia el oeste dentro del subsuelo y en la
la Formación Tena (Maastrichtiano-Paleoceno), que zona Subandina; 3) las areniscas individuales de Napo
corresponden a capas sello de la secuencia cretácica; tienen su desarrollo máximo en el área central y un
su espesor no es uniforme, varían de ~100 m en el adelgazamiento regional hacia el oeste, en la región
Subandino a 500 m en el centro-este de la cuenca. Subandina está representada solo por lutitas marinas
distales y capas delgadas de areniscas glauconíticas; 4)
Características de los fluidos en la cuenca Oriente al este, las secciones de Hollín se van acuñando y las
Las concentraciones de salinidad de las aguas de secciones de Napo son menos profundas y finalmente
formación en los reservorios Cretácicos varían entre se fusionan para convertirse en una capa de areniscas
500 y 130.000 ppm de NaCl equivalente (Smith, indiferenciadas antes de entrar en contacto con las
1989). De acuerdo con Rivadeneira y Almeida (2014), rocas precámbricas expuestas del escudo Guyanés; y
en Hollín Inferior el agua de formación tiene muy 5) las aguas de formación de las areniscas cretácicas
baja salinidad, alrededor de 2000 a 3000 ppm de se caracterizan por gradientes regionales de salinidad
NaCl, mientras que las areniscas de Hollín Superior dentro y entre areniscas individuales.

Figura 5. Sección estratigráfica cretácica en la parte nororiental de la cuenca Oriente. Las flechas representan el movimiento del
agua hacia el este y hacia arriba en areniscas más jóvenes. Tomado de Smith (1989).

130 Boletín de Geología - vol. 45, n.° 1


Elena Angulo-Romero; Bernardo Beate; Christian Romero-Cóndor

Trampas y fallas principales gas en todo el mundo. Por ejemplo, en China, se ha


En la cuenca Oriente principalmente se desarrollaron descubierto que las principales cuencas petrolíferas
trampas estructurales de bajo relieve en la parte son ricas en recursos geotérmicos, como el campo
occidental y central, y trampas estratigráficas en la petrolífero Daqing, Liaohe y Huabei, donde las
parte oriental (Zhongzhen et al., 2017). La cuenca reservas totales fueron de hasta 424 EJ (1 EJ=1018 J)
desarrolló dos grupos de fallas con tendencias NO-SE de recurso geotérmico recuperable (Wang et al., 2016).
y NE-SO, por separado (Yinfu et al., 2010). A partir Estados Unidos es uno de los principales países en
de secciones geológicas regionales, la cuenca Oriente investigación y utilización de la generación de energía
puede dividirse en tres cinturones estructurales de oeste geotermoeléctrica de intermedia a baja temperatura.
a este: el cinturón de empuje del piedemonte occidental
(Cinturón Plegado Sobrecorrido Subandino), el sistema La primera estación de energía geotermoeléctrica del
de rift invertido central a lo largo del Corredor Sacha- mundo se instaló en 2006 en el campo petrolero Teapot
Shushufindi y el Sistema Invertido Capirón-Tiputini Dome en el norte de Wyoming. La potencia instalada
o región de deslizamiento de rumbo (strike-slip) de esta estación con el principio de circulación Organic
conjugado oriental (desde la falla inversa Fanny-Dorine Rankine Cycle (ORC) fue de 0,25 MW (Zhang et al.,
y más al este la falla inversa del ITT) (Baby et al., 2013; 2016). Esta estación genera 0,18 MW de capacidad
Lin et al., 2015). La orientación de los tres cinturones instalada usando 40.000 barriles H2O/día (73,71 L/s) a
tectónicos es casi N-S, donde la parte norte de la cuenca 90,6-98,9°C (Wang et al., 2018).
está más desarrollada estructuralmente que la parte sur.
La segunda estación de energía geotermoeléctrica
De acuerdo con Mancilla et al. (2008), las principales asociada a un campo petrolero se construyó en China,
trampas están asociadas a estructuras anticlinales en el 2011, en el yacimiento petrolífero de Liubei en el
originadas a partir de la inversión tectónica de los Campo Huabei; se instaló una planta que genera 0,40
semigrábenes, las cuales se desarrollan al menos en MW de capacidad instalada usando 2880 m3H2O/día
dos etapas principales de deformación. Durante la (33,33 L/s) a 110°C (Wang et al., 2018). La capacidad
primera etapa (Turoniano superior-Maastrichtiano), de generación de energía esperada de la estación de
coincidente con la depositación de las formaciones Liubei es de 1,5 millones de kWh/año, con un aumento
Napo Superior y Tena, se origina un incipiente anual promedio de 120.000 toneladas de petróleo
estado de deformación compresiva que causa la (Zhang et al., 2016).
inversión tectónica de los antiguos depocentros. Con
la depositación de la Formación Napo, la cuenca Por su parte, en Colombia, se ha construido la primera
tuvo un cambio en la fuente, de marina a no marina. planta geotérmica en el campo Las Maracas, en
Esta etapa representa un acortamiento en la cuenca y Casanare. Tiene una capacidad instalada de 100 kW,
una reactivación de fallas extensionales heredadas e capaz de generar hasta 72.000 kWh/mes (Richter,
inversión de la cuenca, posiblemente vinculadas a la 2021).
acreción del Caribbean Large Ingeous Province en
~73Ma (Gutiérrez et al., 2019). La segunda etapa se La Tabla 2 resume los proyectos que aplican
desarrolla durante el Eoceno inferior a medio, donde la geotermia como usos directos, tanto de pozos
se produce el segundo episodio compresivo y de productores como abandonados en todo el mundo.
inversión en la cuenca (Mancilla et al., 2008). Esto El uso directo geotérmico de campos petroleros
puede relacionarse con el avance del cinturón plegado específicos se lleva a cabo principalmente en campos
sobrecorrido debido a la exhumación de los Andes, que producen petróleo altamente viscoso (Wang
consistente con la evidencia termocronológica durante et al., 2018). Por otro lado, la Tabla 3 resume los
el Eoceno (Spikings et al., 2010). proyectos implementados de generación de energía
geotermoeléctrica en campos petroleros en una planta
de energía binaria que utiliza agua de formación.
Posible aprovechamiento geotérmico en la cuenca
Oriente
Numerosos estudios han demostrado el abundante
recurso geotérmico en los yacimientos de petróleo y

Boletín de Geología - vol. 45, n.° 1 131


Zonificación del gradiente geotérmico en la cuenca Oriente de Ecuador a partir de temperatura de fondo de pozos

Tabla 2. Resumen de proyectos que emplean la geotermia como uso directo de energía en yacimientos petrolíferos. Modificado
de Wang et al. (2018).

Ubicación del Recurso Temperatura y


Tipo de utilización Observaciones
campo petrolero geotérmico caudal

Pozos Calefacción de Capacidad instalada de 4.497 MWt y un consumo


Albania 65,5°C
abandonados invernaderos anual de energía para calefacción de 23,26 TJ/año
Aguas geotérmicas utilizadas en los balnearios
Pozos
Estiria, Australia Complejos de spa 47°C Loipersdorf Waltersdorf. Capacidad instalada de 10
abandonados
MWt y un consumo anual de 247 TJ/año
Generación de agua
Algyő, Hungría
potable 25 L/s a 87°C y Agua geotérmica usada para calentar la tubería de
Productor de
Recuperación una profundidad recuperación en Sávoly. Capacidad instalada de 42
agua
Sávoly, Hungría de petróleo por de 1,8 km MWt y un consumo de 695 TJ/año
trazado térmico
Calentador
Productor de de espacios. Para un área de 120.000 km2 se estima un recurso de
Daqing, China 50-110°C
agua Transportación de 2,905 × 108 J de energía geotermoeléctrica
crudo pesado
Productor de Calentador de Para un área de 25.000 km2 se estima un recurso de
Liaohe, China 50-110°C
agua espacios 1,008 × 108 J de energía geotermoeléctrica
Recuperación
Pozos 600 m3/día a Para un área de 32.000 km2 se estima un recurso de
Huabei, China de petróleo por
abandonados 100-110°C 7,099 × 108 J de energía geotermoeléctrica
trazado térmico

Tabla 3. Proyectos de generación de energía geotermoeléctrica en campos petroleros implementados notablemente. Modificado
de Wang et al. (2018).

Caudal Temperatura del agua Potencia


Ubicación Observaciones
(barriles/día) (°C) instalada (kW)
Produjo más de 2210 MWh en los primeros
Wyoming, USA 40.000 (73,6 L/s) 90,6-98,9 132
3,5 años
Primer proyecto comercial de generación
Dakota del Norte,
30.000 (55,2 L/s) 98 250 de energía geotermoeléctrica en un pozo de
USA
petróleo y gas en EE. UU

Huabei, China 18.114 (33,3 L/s) 110 310 Generó ~31 × 10⁴ kWh/año a finales de 2011

De acuerdo con Moeck (2014), las partes someras de apliquen tecnologías para mejorar la productividad. Por
las cuencas sedimentarias (<3 km) con anomalías de su parte, Liu et al. (2018) consideran que la temperatura
calor superiores al gradiente geotérmico promedio de fondo debe encontrarse entre 65-150°C para que un
pueden ser aprovechadas directamente. En el caso de campo petrolero sea viable para posibles aplicaciones
la cuenca Oriente, se han podido identificar dos zonas geotérmicas. Se estima que las tecnologías actualmente
con gradientes geotérmicos anómalos, por lo general, disponibles requieren al menos 15.000 barriles/día
asociados a profundidades someras, que podrían (27,6 L/s) de agua de formación con una temperatura
tener un aprovechamiento directo para sistemas de mínima de 101,7°C para una generación de energía
calefacción y refrigeración urbana, sistemas de clima eléctrica eficiente y económica en campos petroleros
controlado (como en los invernaderos), deshidratación para aplicaciones geotérmicas (Liu et al., 2018).
de productos agrícolas y varios procesos industriales
adicionales. Dichas zonas anómalas se encuentran en En un informe presentado por PETROECUADOR
los campos individuales de Pungarayacu-San Carlos y (2012), la producción de aguas de formación en las
Cowi (Tabla 4). áreas donde opera (Auca, Shushufindi, Lago Agrio,
Libertador y Cuyabeno) fue de 288.841 barriles/día
De igual manera, Moeck (2014) considera que las partes (∼531,55 L/s), superando el caudal que se considera
más profundas de las cuencas (>3 km) con rocas de baja viable para la producción de energía eléctrica en
porosidad/permeabilidad pueden ser adecuadas para la campos petroleros.
producción de energía eléctrica siempre y cuando se

132 Boletín de Geología - vol. 45, n.° 1


Tabla 4. Factibilidad de uso de recursos geotérmicos en la cuenca Oriente-Ecuador.
Boletín de Geología - vol. 45, n.° 1

Ubicación Temperatura,
Recurso Tipo, configuración tectónica y ambiente
del campo gradiente Tipo de aprovechamiento
geotérmico geológico
petrolero geotérmico y caudal
Tipo conductivo
BHT: 30-110°C
Pungarayacu- Pozos petroleros asociado a un Acuíferos someros de baja (0,270-4 km); GGA: Sistemas de calefacción/refrigeración urbana,
San Carlos exploratorios cinturón plegado salinidad en la arenisca 17,02-69,01°C/km sistemas de calefacción con clima controlado
sobrecorrido porosa y permeable de la
(invernaderos), deshidratación de productos
Formación cretácica Hollín
Tipo conductivo BHT: 130°C a 0,385 agrícolas y varios procesos industriales
(<2 km) controlado por
Cowi Pozo abandonado asociado a una km; GGA: 66,3°C/ adicionales

Elena Angulo-Romero; Bernardo Beate; Christian Romero-Cóndor


fracturas y fallas
cuenca de antepaís km
BHT: 97,57-
118,68°C (3-3,5
Tipo conductivo
Pozos productores km); GGA: 23,49-
Lago Agrio asociado a una
de petróleo 30,29°C/km; Q:
cuenca de antepaís
531,55 L/s, en total
1’595.558 bles/año
BHT: 74,09-
113,14°C (1,6-3,5
Tipo conductivo
Pozos productores km); GGA: 19,37-
Libertador asociado a una
de petróleo 31,06°C/km; Q:
cuenca de antepaís
531,55 L/s, en total
26’054.091 bles/año
Acuíferos profundos de BHT: 86,32-
Producción de energía eléctrica. Se estima que
baja salinidad en la arenisca 115,23°C (2,5-2,7
Tipo conductivo la tecnología actualmente disponible requiere al
Pozos productores porosa y permeable de la km); GGA: 20,72-
Cuyabeno asociado a una menos un caudal de 15.000 bles/día (~27,6 L/s) de
de petróleo Formación cretácica Hollín 36,69°C/km; Q:
cuenca de antepaís agua de formación con una temperatura mínima
(>3 km) controlado por 531,55 L/s, en total
de 101,7°C
fracturas y fallas 26’724.123 bles/año
BHT: 74,09-
113,14°C (1,6-3,5
Tipo conductivo
Pozos productores km); GGA: 19,37-
Shushufindi asociado a una
de petróleo 31,06°C/km; Q:
cuenca de antepaís
531,55 L/s, en total
34’487.127 bles/año
BHT: 83,91-
127,14°C (2,6-3,7
Tipo conductivo
Pozos productores km); GGA: 15,04-
Auca asociado a una
de petróleo 32,96°C/km; Q:
cuenca de antepaís
531,55 L/s, en total
16’864.860 bles/año
133
Zonificación del gradiente geotérmico en la cuenca Oriente de Ecuador a partir de temperatura de fondo de pozos

Si se consideran únicamente las temperaturas de fondo Corredor Capirón-Tiputini, pueden estar asociadas a
de pozos y la producción de aguas de formación, concentraciones anómalas de elementos radiactivos en
estás se encuentran dentro de los rangos establecidos el basamento.
por Liu et al. (2018). Como solo se pudo obtener
información de producción de agua de formación en En la cuenca Oriente se ha considerado a la Formación
los campos Auca, Shushufindi, Lago Agrio, Libertador Hollín como el principal reservorio geotérmico debido a
y Cuyabeno, se considera que ellos podrían tener sus condiciones favorables de porosidad/permeabilidad
un aprovechamiento para la generación de energía y baja salinidad. Por su parte, las areniscas de Napo
eléctrica y abastecer a las redes energéticas principales se consideran como unos reservorios locales, debido
o dotar de energía eléctrica a las zonas de influencia a su restringida distribución en la cuenca; la Arenisca
directa donde operan. Para este estudio no se ha “T” se extiende en la parte central y oriental de la
considerado la infraestructura mínima requerida para cuenca, por su parte las Areniscas “U” y “M1” se
implementar una aplicación geotérmica en campos restringen solo a la parte más oriental de la cuenca.
petroleros ecuatorianos, es un aporte inicial que puede El agua ingresa (zona de recarga) a la cuenca a través
ser tomado en cuenta para futuras evaluaciones. de las secciones elevadas, y se descarga a lo largo del
borde poco profundo oriental de la cuenca. Las lutitas
Conclusiones de Tena Basal y Napo son el sello regional principal.
La deformación tectónica que sufrió la cuenca afectó a
La ley de variación de la temperatura con la profundidad las rocas de cobertura y basamento, y está relacionada
para la cuenca Oriente es TEQ=20,14z + 41,87. Esta ley con la inversión de los sistemas de fallas extensionales
es la que consigue la mejor correlación lineal de los precretácicas, durante tres períodos de inversión
pares de datos temperatura/profundidad de la cuenca. (Cretácico tardío-Paleoceno, Eoceno temprano
Esta ecuación permite modelar la temperatura a una y Mioceno), desarrollando trampas estructurales
profundidad constante y es un promedio para toda la de bajo relieve en la parte occidental y central, y
cuenca, asumiendo un régimen conductivo y regional algunas trampas estratigráficas en la parte oriental.
para la transferencia de calor. Adicionalmente, Hollín Superior es estratigráfico al
oeste de la cuenca, “T” y “U” son estratigráficos en la
En el mapa de gradientes geotérmicos, se hace evidente parte central de la cuenca y “M1” es estratigráfico al
que el Corredor Sacha-Shushufindi muestra los valores este de la cuenca. La existencia de estas fallas puede
más bajos de gradiente geotérmico (7,38-26,1°C/km). aumentar o disminuir la permeabilidad, y de alguna
Por el contrario, los valores más altos de gradiente manera puede influir en la compartimentación del
geotérmico (26,1-69,01°C/km) están presentes al reservorio, generando accesos al recurso geotérmico
occidente (Levantamiento Napo) y al oriente (Corredor ubicados en niveles permeables.
Capirón-Tiputini).
Las zonas anómalas de la cuenca Oriente podrían
En la cuenca Oriente se identificaron con claridad tener un aprovechamiento directo para sistemas
cuatro zonas anómalas, las cuales están localizadas en de calefacción/refrigeración urbana, sistemas de
los campos individuales de Pungarayacu-San Carlos y clima controlado en invernaderos, deshidratación de
Bermejo-El Rayo ubicados en el Levantamiento Napo productos agrícolas y varios procesos industriales
hacia el borde occidental de la cuenca, Cowi (borde en las zonas anómalas de Pungarayacu-San Carlos y
centro-occidental del Corredor Capirón-Tiputini) e Cowi, debido a la existencia de acuíferos someros en la
Ishpingo-Tambococha-Tiputini (borde más oriental arenisca porosa de Hollín a profundidades <2 km con
del Corredor Capirón-Tiputini). La presencia de estas temperaturas de fondo que oscilan entre los 30-90°C.
anomalías positivas en el Levantamiento Napo está Incluso, pueden ser adecuados para la producción
fuertemente influenciada por la cercanía del basamento, de energía eléctrica, siempre y cuando el caudal y la
debido a procesos tectónicos (levantamiento de las temperatura sean eficientes (>27,6 L/s y >101,7°C) y
formaciones y erosión de orógenos), además, los se apliquen tecnologías para mejorar la productividad
campos Pungarayacu-San Carlos y Bermejo-El Rayo del campo, ya que al aumentar la profundidad,
se encuentran cerca de los volcanes activos Sumaco disminuye la permeabilidad de las rocas. De acuerdo
(alcalino) y Reventador (calco-alcalino) por lo que con la información disponible se puede considerar
pueden considerarse como probables fuentes de calor como prospectos a las áreas de Lago Agrio, Libertador,
adicional que abastecen al sistema geotérmico. Por Cuyabeno, Shushufindi y Auca.
su parte, las anomalías de calor perceptibles en el

134 Boletín de Geología - vol. 45, n.° 1


Elena Angulo-Romero; Bernardo Beate; Christian Romero-Cóndor

Agradecimientos Beate, B.; Urquizo, M.; Lloret, A. (2020). Geothermal


Country Update for Ecuador, 2015–2020. World
Al Banco de Información Petrolera del Ecuador (BIPE) Geothermal Congress, Reykjavik, Iceland.
y la Secretaría de Hidrocarburos, en especial al Ing.
César Andrade, por otorgar acceso a la base de datos. Bès de Berc, S.; Soula, J.C.; Baby, P.; Souris, M.;
Christophoul, F.; Rosero, J. (2005). Geomorphic
Referencias evidence of active deformation and uplift in
Alfaro, C.; Alvarado, I.; Manrique, A. (2015). Heat a modern continental wedge-top–foredeep
Flow Evaluation at Eastern Llanos Sedimentary transition: example of the eastern Ecuadorian
Basin, Colombia. World Geothermal Congress, Andes. Tectonophysics, 399(1-4), 351-380.
Melbourne, Australia. https://doi.org/10.1016/j.tecto.2004.12.030

Allen, D.M.; Grasby, S.E.; Voormeij, D.A. (2006). Brookfield, M.E.; Hemmings, D.P.; Van Straaten,
Determining the circulation depth of thermal P. (2009). Paleoenvironments and origin of the
springs in the southern Rocky Mountain Trench, sedimentary phosphorites of the Napo Formation
south-eastern British Columbia, Canada using (Late Cretaceous, Oriente Basin, Ecuador).
geothermometry and borehole temperature logs. Journal of South American Earth Sciences,
Hydrogeology Journal, 14(1-2), 159-172. https:// 28(2), 180-192. https://doi.org/10.1016/j.
doi.org/10.1007/s10040-004-0428-z jsames.2009.02.004

Aspden, J.A.; Ivimey-Cook, H.C. (1992). Nuevos datos Bulot, L.G.; Kennedy, W.J.; Jaillard, E.; Robert, E.
paleontológicos del Centro y Sureste del Ecuador. (2005). Late Middle–early Late Albian ammonites
Boletín Geológico Ecuatoriano, 3, 33-42. from Ecuador. Cretaceous Research, 26(3), 450-
459. https://doi.org/10.1016/j.cretres.2005.01.008
Baby, P.; Rivadeneira, M.; Barragán, R.; Christophoul,
F. (2013). Thick-skinned tectonics in the Oriente Burgos, J.; Baby, P.; Christophoul, F.; Soula, J.;
foreland basin of Ecuador. Geological Society, Rochat, P. (2004). Cuantificación de las erosiones
London, Special Publications, 377, 59-76. https:// Terciarias y Plio-Cuaternarias en la parte sur de
doi.org/10.1144/SP377.1 la Cuenca Oriente. En: P. Baby, M. Rivadeneira,
R. Barragán (eds.). La Cuenca Oriente:
Barba, D. (2017). Gradiente geotermal de las cuencas Geología y Petróleo (pp. 115-130). IFEA-IRD-
de ante-arco del Ecuador. VIII Jornadas en Petroamazonas EP.
Ciencias de la Tierra, Quito, Ecuador.
Canfield, R.W.; Bonilla, G.; Robbins, R.K. (1982). Sacha
Barba, D.; Barragán, R.; Gallardo, J.; Ormasa, A.; oil field of Ecuadorian Oriente. AAPG Bulletin, 66(8),
Salguero, A. (2021). Geothermal Gradients in the 1076-1090. https://doi.org/10.1306/03B5A642-
Upper Amazon Basin derived from BHT data. 16D1-11D7-8645000102C1865D
International Journal of Terrestrial Heat Flow
and Applied Geothermics, 4(1), 85-94. https://doi. Carvalho, H.D.S.; Vacquier, V. (1977). Method
org/10.31214/ijthfa.v4i1.58 for determining terrestrial heat flow in oil
fields. Geophysics, 42(3), 584-593. https://doi.
Beardsmore, G.R.; Cull, J.P. (2001). Thermal Gradient. org/10.1190/1.1440729
En: Crustal Heat Flow – A Guide to Measurement
and Modelling (pp. 47-87). Cambridge University Christophoul, F. (1998). Comisión zona Gualaquiza,
Press. carretera Méndez-Patuca-Santiago-Morona y
Cerro Macuma. Orstom-Petroproducción, Quito.
Beate, B.; Urquizo, M. (2015). Geothermal country
update for Ecuador: 2010–2015. World Christophoul, F.; Baby, P.; Dávila, C. (2002).
Geothermal Congress, Melbourne, Australia. Stratigraphic responses to a major tectonic event
in a foreland basin: the Ecuadorian Oriente from

Boletín de Geología - vol. 45, n.° 1 135


Zonificación del gradiente geotérmico en la cuenca Oriente de Ecuador a partir de temperatura de fondo de pozos

Eocene to Oligocene times. Tectonophysics, Tectonics (pp. 237-268). Elsevier. https://doi.


345(1-4), 281-298. https://doi.org/10.1016/ org/10.1016/B978-0-12-816009-1.00011-3
S0040-1951(01)00217-7
Hamza, V.; Silva, F.; Gomes, A.; Delgadilho, Z.
Clauser, C. (2009). Heat Transport processes in the (2005). Numerical and functional representations
Earth’s Crust. Surveys in Geophysics, 30(3), 163- of regional heat flow in South America. Physics
191. https://doi.org/10.1007/s10712-009-9058-2 of the Earth and Planetary Interiors, 152(4), 223-
256. https://doi.org/10.1016/j.pepi.2005.04.009
Dashwood, M.F.; Abbotts, I.L (1990). Aspects of
the petroleum geology of the Oriente Basin, Higley, D. (2001). The Putumayo-Oriente-Maranon
Ecuador. Geological Society, London, Special Province of Colombia, Ecuador, and Peru-
Publications, 50, 89-117. https://doi.org/10.1144/ Mesozoic-Cenozoic and Paleozoic Petroleum
GSL.SP.1990.050.01.06 Systems. U.S. Geological Survey Digital Data.

Deming, D. (1989). Application of bottom-hole Ibañez-Mejia, M.; Ruiz, J.; Valencia, V.; Cardona,
temperature corrections in geothermal studies. A.; Gehrels, G.; Mora, A. (2011). The Putumayo
Geothermics, 18(5-6), 775-786. https://doi. Orogen of Amazonia and its implications
org/10.1016/0375-6505(89)90106-5 for Rodinia reconstructions: New U-Pb
geochronological insights into the Proterozoic
Díaz, M.; Baby, P.; Rivadeneira, M.; Christophoul, tectonic evolution of northwestern South America.
F. (2004). El pre-Aptense de la cuenca Oriente Precambrian Research, 191(1-2), 58-77. https://
ecuatoriana. En: P. Baby, M. Rivadeneira, R. doi.org/10.1016/j.precamres.2011.09.005
Barragán (eds.). La Cuenca Oriente, Geología y
Petróleo (pp. 23-44). IFEA-IRD-Petroamazonas Jaillard, E.; Caron, M.; Dhondt, A.; Ordóñez, M.;
EP. Andrade, R.; Bengtson, P.; Bulot, L.; Cappeta,
H.; Dávila, C.; Díaz, R.; Huacho, J.; Huamán,
ESMAP (2012). Manual de Geotermia: cómo planificar C.; Jiménez, D.;Jiménez, N.; Montenegro, J.;
y financiar la generación de electricidad. Néraudeau, D.; Rivadenerira, M.; Toro, J.;
Programa de Asistencia para la Gestión del Sector Villagómez, R.; Zambrano, Í. (1997). Síntesis
Energético (ESMAP) y Grupo del Banco Mundial. estratigráfica y sedimentológica del Cretáceo y
Paleógeno de la cuenca oriental del Ecuador.
Estupiñán, J.; Marfil, R.; Permanyer, A. (2006). Orstom-Petroproducción.
Diagénesis y estratigrafía secuencial de la
arenisca “U” de la Fm Napo del Oriente, Ecuador. Lin, J.; Zhang, X.; HU, Q.; Zhao, X.; Good, D.;
Geogaceta, 40, 283-286. Tong, X.; Orozco, R.; Zhang, Z; Xie, Y. (2015).
Model experiments to simulate compressional-
Gaibor, J.; Hochuli, J.P.A.; Winkler, W.; Toro, J. (2008). strike fault pattern in Oriente Basin, Ecuador.
Hydrocarbon source potential of the Santiago SPE Latin American and Caribbean Petroleum
Formation, Oriente Basin, SE of Ecuador. Journal Engineering Conference, Quito, Ecuador. https://
of South American Earth Sciences, 25(2), 145-156. doi.org/10.2118/177175-MS
https://doi.org/10.1016/j.jsames.2007.07.002
Litherland, M.; Aspden, J.A.; Jemielita, R.A. (1994).
Goldschmidt, K. (1941). Geological Compilation on Metamorphic belts of Ecuador. Overseas Geology
the Cutucu Macuma-Cangaime Region. SHELL- and Mineral Resources. British Geological
Petroproducción. Survey.

Gutiérrez, E.; Horton, B.; Vallejo, C.; Jackson, Liu, X.; Falcone, G.; Alimonti, C. (2018). A
L.; George, S. (2019). Provenance and systematic study of harnessing low-temperature
geochronological insights into Late Cretaceous- geothermal energy from oil and gas reservoirs.
Cenozoic foreland basin development in the Energy, 142, 346-355. https://doi.org/10.1016/j.
Subandean Zone and Oriente Basin of Ecuador. energy.2017.10.058
En: B. Horton, A. Folguera (eds.). Andean

136 Boletín de Geología - vol. 45, n.° 1


Elena Angulo-Romero; Bernardo Beate; Christian Romero-Cóndor

Luzieux, L.D.A.; Heller, F.; Spikings, R.; Vallejo, produced fluids from oil production. https://www.
C.F.; Winkler, W. (2006). Origin and Cretaceous thinkgeoenergy.com/first-geothermal-power-
tectonic history of the coastal Ecuadorian forearc plant-inaugurated-in-colombia/
between 1°N and 3°S: Paleomagnetic, radiometric
and fossil evidence. Earth and Planetary Rivadeneira, M. (1986). Evaluación geoquímica
Science Letters, 249(3-4), 400-414. https://doi. de rocas madres de la Cuenca Amazónica
org/10.1016/j.epsl.2006.07.008 Ecuatoriana. IV Congreso Ecuatoriano de
Geología, Minas y Petróleo, Quito, Ecuador.
Macellari, C.E. (1988). Cretaceous paleogeography
and depositional cycles of western South America. Rivadeneira, M.; Baby, P. (1999). La Cuenca
Journal of South American Earth Sciences, Oriente: estilo tectónico, etapas de deformación
1(4), 373-418. https://doi.org/10.1016/0895- y características de los principales campos
9811(88)90024-7 de Petroproducción. Informe técnico. IRD-
Petroamazonas.
Mancilla, O.; Albariño, L.; Meissinger, V.; Rivadeneira,
M.; Sciamanna, S. (2008). Sistemas petroleros Rivadeneira, M.; Baby, P. (2004). Características
de la cuenca Oriente, Ecuador. VII Congreso de geológicas generales de los principales campos
Exploración y Desarrollo de Hidrocarburos, Mar petroleros de Petroproducción. En: P. Baby,
del Plata, Argentina. M. Rivadeneira, R. Barragán (eds.). La Cuenca
Oriente, Geología y Petróleo (pp. 229-295).
Martin-Gombojav, N.; Winkler, W. (2008). Recycling IFEA-IRD-Petroamazonas EP.
of Proterozoic crust in the Andean Amazon
foreland of Ecuador: implications for orogenic Rivadeneira, M.; Almeida, P. (2014). Características de
development of the Northern Andes. Terra Nova, los reservorios Cretácicos de la cuenca Oriente.
20(1), 22-31. https://doi.org/10.1111/j.1365- En: P. Baby, M. Rivadeneira, R. Barragán (eds.).
3121.2007.00782.x La Cuenca Oriente, Geología y Petróleo (pp.
279-326). IFEA-IRD-Petroamazonas EP.
Mello, M.R.; Koutsoukos, E.; Erazo, W. (1995).
The Napo Formation, Oriente Basin, Romero, C. (2018). Identificación y caracterización de
Ecuador: Hydrocarbon source potential and Facies de la Formación Hollín en Centro Shaime:
paleoenvironmental assessment. En: Katz, B.J. El registro de una transición fluvio-marina en la
(eds.). Petroleum Source Rocks (pp. 167-181). región Sur Oriental del Ecuador. Tesis, Escuela
https://doi.org/10.1007/978-3-642-78911-3_10 Politécnica Nacional, Quito, Ecuador.

Moeck, I. (2014). Catalog of geothermal play types Romero, C.; Rivadeneira, M.; Calderon, E.; Naranjo,
based on geologic controls. Renewable and M.; Meneses, M.; Gramal, A.; Toainga, S. (2019).
Sustainable Energy Reviews, 37, 867-882. https:// Reconstruyendo el pasado del planeta: El registro
doi.org/10.1016/j.rser.2014.05.032 estratigráfico y sedimentológico de la Formación
Tena en la carretera Tiwintza–Puerto Morona.
Ordóñez, M.; Jiménez, N.; Suárez, J. (2006). GeoLatitud, 2(2), 2-13.
Micropaleontologia ecuatoriana: Datos
bioestratigráficos y paleoecológicos de las Romero, C.; Escobar, V.; Calderón, D.; Menéndez,
cuencas: Graben de Jambelí. Progreso. Manabí, B.; Gallardo, A.; Gramal, A.; Vélez, T.; Condoy,
Esmeraldas y Oriente; del levantamiento de la D. (2021). La evolución de abanicos aluviales
Península de Santa Elena, y de las cordilleras documentada en el registro estratigráfico de la
Colonche, costera y occidental. Petroproducción. Formación Tiyuyacu. GeoLatitud, 4(1), 15.

PETROECUADOR EP. (2012). Informe estadístico de Romeuf, N.; Aguirre, L.; Soler, P.; Feraud, G.; Jaillard,
la industria hidrocarburífera Ecuatoriana 1972- E.; Ruffet, G. (1995). Middle Jurassic volcanism
2012. in the Northern and Central Andes. Andean
Geology, 22(2), 245-259.
Richter, A. (2021). The first geothermal power unit
has started operations in Colombia, utilizing co-

Boletín de Geología - vol. 45, n.° 1 137


Zonificación del gradiente geotérmico en la cuenca Oriente de Ecuador a partir de temperatura de fondo de pozos

Ruiz, G.M.H.; Seward, D.; Winkler W. (2007). M.; Leal, J.; Cuenca, D. (2017). Geology of
Evolution of the Amazon Basin in Ecuador with the Campanian M1 sandstone oil reservoir of
special reference to hinterland tectonics: data from eastern Ecuador: A delta system sourced from the
zircon fission-track and heavy mineral analysis. Amazon Craton. Marine and Petroleum Geology,
Developments in Sedimentology, 58, 907-934. 86, 1207-1223. https://doi.org/10.1016/j.
https://doi.org/10.1016/S0070-4571(07)58036-2 marpetgeo.2017.07.022

Saemundsson, K.; Axelsson, G.; Straingrimsson, B. Vallejo, C.; Romero, C.; Horton, B.K.; Spikings, R.A.;
(2009). Geothermal system in global perspective. Gaibor, J.; Winkler, W.; Esteban, J.J.; Thomsen,
Short Course on Geothermal Drilling, Resource T.B.; Mariño, E. (2021). Jurassic to Early
Development and Power Plants, El Salvador. Paleogene sedimentation in the Amazon region
of Ecuador: Implications for the paleogeographic
Sass, I.; Götz, A.E. (2012). Geothermal reservoir evolution of northwestern South America.
characterization: a thermofacies concept. Terra Global and Planetary Change, 204. https://doi.
Nova, 24(2), 142-147. https://doi.org/10.1111/ org/10.1016/j.gloplacha.2021.103555
j.1365-3121.2011.01048.x
Vera, R. (2016). Geology of Ecuador: An introduction
Shanmugam, G.; Poffenberger, M.; Toro, J. (2000). to the unique geology of Ecuador. 2da edición.
Tide-dominated estuarine facies in the Hollin and Elsevier.
Napo (“T” and “U”) Formations (Cretaceous),
Sacha field, Oriente basin, Ecuador. AAPG Wang, S.; Yan, J.; Li, F.; Hu, J.; Li, K. (2016).
Bulletin, 84(5), 652-682. https://doi.org/10.1306/ Exploitation and utilization of oilfield geothermal
C9EBCE7D-1735-11D7-8645000102C1865D resources in China. Energies, 9(10). https://doi.
org/10.3390/en9100798
Smith, L. (1989). Regional Variations in Formation
Water Salinity, Hollin and Napo Formations Wang, K.; Yuan, B.; Ji, G.; Wu, X. (2018). A
(Cretaceous), Oriente Basin, Ecuador. comprehensive review of geothermal energy
AAPG Bulletin, 73(6), 757-776. https:// extraction and utilization in oilfields. Journal of
d o i . o rg / 1 0 . 1 3 0 6 / 4 4 B 4 A 2 5 8 - 1 7 0 A - 11 D 7 - Petroleum Science and Engineering, 168, 465-
8645000102C1865D 477. https://doi.org/10.1016/j.petrol.2018.05.012

Spikings, R.A.; Crowhurst, P.V.; Winkler, W.; White, H.J.; Skopec, R.; Ramírez, F.; Rodas, J.;
Villagomez, D. (2010). Syn- and post- Bonilla, G. (1995). Reservoir characterizations
accretionary cooling history of the Ecuadorian of the Hollin and Napo formations, western
Andes constrained by their in-situ and detrital Oriente basin, Ecuador. En: A.J. Tankard, R.
thermochronometric record. Journal of South Suárez, H.J. Welsink (eds.). Petroleum basins of
American Earth Sciences, 30(3-4), 121-133. South America (pp. 573-596). AAPG. https://doi.
https://doi.org/10.1016/j.jsames.2010.04.002 org/10.1306/M62593C30

Tschopp, H.J. (1953). Oil explorations in the Oriente Willet, S.; Chapman, D. (1987). Analysis of
of Ecuador, 1938-1950. AAPG Bulletin, 37(10), Temperatures and Thermal Processes in the
2303-2347. https://doi.org/10.1306/5CEADD94- Uinta Basin. En: C. Beaumont y A.J. Tankard
16BB-11D7-8645000102C1865D (eds.). Sedimentary Basins and Basin Forming
Mechanisms (pp. 447-461). AAPG.
Vallejo, C.; Hochuli, P.; Winkler, W.; von Salis, K.
(2002). Palynological and sequence stratigraphic Yinfu, X.; Jihancheng; Yongdi, S.; Ying, H. (2010).
analysis of the Napo Group in the Pungarayacu Petroleum geology and exploration potential of
30 well, Sub-Andean Zone, Ecuador. Cretaceous Oriente-Maranon Basin. Petroleum Exploration
Research, 23(6), 845-859. https://doi.org/10.1006/ and Development, 37(1), 51-56. https://doi.
cres.2002.1028 org/10.1016/S1876-3804(10)60014-6

Vallejo, C.; Tapia, D.; Gaibor, J.; Steel, R.; Cardenas, Zhang, J.; Cho, H.; Knizley, A. (2016). Evaluation of
M.; Winkler, W.; Valdez, A.; Esteban, J.; Figuera, financial incentives for combined heat and power

138 Boletín de Geología - vol. 45, n.° 1


Elena Angulo-Romero; Bernardo Beate; Christian Romero-Cóndor

(CHP) systems in U.S. regions. Renewable and hydrocarbon plays in Putumayo-Oriente-Marañon


Sustainable Energy Reviews, 59, 738-762. https:// Basin, South America. Petroleum Exploration
doi.org/10.1016/j.rser.2016.01.012 and Development, 44(2), 247-256. https://doi.
org/10.1016/S1876-3804(17)30027-7
Zhongzhen, M.A.; Chen, H.; Yinfu, X.; Zhiwei, Z.;
Yaming, L.; Xiaofa, Y.; Yubing, Z.; Dandan, W. Fecha de recibido: 13 de septiembre de 2021
(2017). Division and resources evaluation of Fecha de aprobado: 15 de noviembre de 2022

Boletín de Geología - vol. 45, n.° 1 139

También podría gustarte