TFM Tol Car PDF
TFM Tol Car PDF
TFM Tol Car PDF
CARTAGENA
A todos los compañeros que han colaborado y ayudado para que pueda realizar
este Máster sin ser Ingeniero, sólo un maestro y licenciado en Matemáticas.
A mi familia Rosa, Esther, Laura, Pedro Manuel y Rosi que ha hecho posible mi
dedicación durante muchos días a este Trabajo Fin de Máster.
2
ÍNDICE
3
CAPÍTULO 6. CAMPO DE TEMPERATURAS Y FLUJOS EN POZOS EN LAS CERCANIAS DE FOCOS
TÉRMICOS. IMPLEMENTACIÓN EN FAHET ......................................................................................... 65
6.1 Desarrollo del trabajo ....................................................................................................................... 65
6.2 Características de los dominios de simulación, condiciones de frontera y datos de entrada.......... 67
6.3 Resultados ......................................................................................................................................... 82
CAPÍTULO 7. CONCLUSIONES.......................................................................................................... 250
CAPÍTULO 8. BIBLIOGRAFIA ............................................................................................................ 253
4
CAPITULO 1. INTRODUCCIÓN
Sin duda la presencia del recurso geotérmico y la magnitud del depósito es la primera
motivación de desarrollo de las tecnologías y de los métodos de aprovechamiento de energía
geotérmica. En este caso, el estado del arte resulta diferentemente caracterizado. A escala mundial, si
hablamos de geotermia de alta entalpía, los principales países de referencia en el sector geotérmico son
USA (por una capacidad eléctrica producida de 3440 MW sobre 12000 MW), Filipinas, Nueva Zelanda,
México y países Europeos como Islandia, Italia; para sistemas de media alta entalpía, Turquía, Filipinas,
Alemania y Francia; y Noruega y Suecia para sistemas de baja entalpía [11] y [12].
5
1.2 Objetivos
Los objetivos específicos del presente Trabajo de Fin de Máster son estudiar numéricamente,
con ayuda de un software específicamente desarrollado para el estudio de problemas relacionados con
los procesos de transmisión de calor y masa en medios porosos, el tema de energía geotérmica
propuesto. En particular, se abordarán las siguientes cuestiones, que se ha relacionado con las distintas
actividades reflejadas más adelante en la planificación y el diagrama Gantt:
1. Revisión bibliográfica en bases de datos especializadas y estado del arte en el campo del cálculo
numérico en Geotermia. (Actividades 1 y 2)
2. Descripción de los fundamentos teóricos que rigen los procesos de flujo de fluidos con
transporte de calor asociado. (Actividad 3)
3. Estudio y diseño mediante FAHET de los modelos en red correspondientes al problema que se
va a abordar. (Actividad 4)
4. Realización de experimentos con el programa FAHET, así como con el resto de programas en los
que se apoya, para hallar la solución numérica del modelo en red empleado para el cálculo y
presentación de los resultados obtenidos. Con las simulaciones se encontrarán las respuestas
del pozo y del dominio al variar de las características geométricas y físico-químicas del depósito
geotérmico. (Actividades 5 y 6)
5. Discusión de los resultados obtenidos y presentación de conclusiones. (Actividades 7 y 8)
6
1.3 Planificación del proyecto. Diagrama “Gantt”
En la tabla siguiente se presenta la planificación seguida para la realización del Trabajo de Fin de
Máster, y se incluye la planificación en diagrama de “Gantt”.
7
1.4 Estructura de la memoria
La memoria está divida en 8 capítulos:
1. Introducción: en esta parte se introducen los objetivos principales del trabajo, la planificación
de proyecto, con esquema en tablas y en forma de diagrama “Gantt”, donde se resumen todas
las actividades, las fechas de inicio y finalización de cada una de ellas, y el trabajo en su
totalidad.
2. Energía geotérmica: se estudia y explica las principales características de los recursos
geotérmicos, el potencial geotérmico y se hace una caracterización y definición de los
principales factores que intervienen en el problema de cálculo numérico que se aborda y que
son necesarios para definir las principales variables empleadas en la implementación de
modelos en FAHET.
3. Fundamentos teóricos que rigen los procesos de flujo de fluidos en el transporte de calor: se
estudian y abordan los principales procesos físicos sobre la transmisión de calor introduciendo
las principales propiedades y ecuaciones que gobiernan el problema.
4. Métodos de simulación de problemas en Geotermia: en este capítulo se explican los métodos
de simulación numérica más empleados y se presenta un resumen del estado del arte,
definiendo y abordando en particular el método de simulación por redes eléctricas, que es la
base de desarrollo del programa FAHET, y el método a las diferencias finitas.
5. El programa FAHET: se explica el funcionamiento del programa y las técnicas de
implementación de dominios geotérmicos y los modos para determinar las variables.
6. Campo de temperaturas y flujos en pozos en las cercanías de focos térmicos. Implementación
en FAHET: se solucionan y explican los casos que forman objeto de estudio de este trabajo y se
presentan y discuten los resultados obtenidos. Se abordan distintos casos de configuración del
sistema.
7. Conclusiones: se presenta un resumen con los resultados más destacables obtenidos tras la
realización de los experimentos.
8. Bibliografía: se incluye toda la bibliografía usada para la realización de este trabajo.
8
CAPÍTULO 2 ENERGÍA GEOTÉRMICA
2.1 Introducción
En este capítulo del TFM se aborda el estado del arte en el campo de la energía geotérmica
enumerando sus principales características y utilización. Se describe el potencial geotérmico distribuido
en el mundo, las variables que influyen sobre la respuesta térmica de cada sitio, y el estado del arte en
el campo de la resolución de problemas numéricos, distinguiendo y explicando algunos de los diversos
programas de cálculo numérico comerciales y de investigación existentes.
Hay otros sistemas y aplicaciones que aprovechan potenciales geotérmicos de baja entalpía o
pequeño gradiente térmico, que suelen ser para profundidades del orden 5-20 m, y en los que el suelo
se emplea como medio de aislamiento, con la intención de mantener una temperatura constante
durante todo el año. Sin entrar en argumentaciones que no son objeto de este trabajo, podemos
9
describir la variación de la temperatura con la profundidad 𝑥(m), y el tiempo t (días), representada en la
figura 2.2 y 2.4, a través de la formulación descrita por Van Wijk y de Vries [51] y [53]:
𝑥 𝑧 𝑥 𝑥
𝑇 𝑥, 𝑡 = 𝑇𝑎 + 𝐴𝑒𝑥𝑝 − sin 𝜔𝑡 − + 𝜙 + 𝐴𝑑 exp − sin 𝜔𝑡 − + 𝜙𝑑
𝐷 𝐷 𝐷𝑑 𝐷
Como puede apreciarse, la respuesta del suelo depende de la temperatura media ambiental
durante un año (𝑇𝑎 ) en ºC, su variación a lo largo del año (A =𝑇max- 𝑇𝑎 ), de características térmicas del
suelo como “Damping” (D) calculada como D= (2 𝑘/ 𝜔) ½, la difusividad térmica (Dd), y la frecuencia
radial 𝜔[rad/s] (2𝜋/365dias), y 𝜙 que es la fase 𝑘 ∗ 𝜔 ∗ 𝑡 y 𝜙𝑑 es la fase calculada en función de la
profundidad.
Figura 2.2 Variación de la temperatura en función del tiempo a lo largo de un año para diferentes profundidades en la
ciudad de Ottawa, Canadá (de Geothermal heat pump systems:Status review and comparison with other heating options,
Stuart J. Self, Bale V. Reddy, Marc A. Rosen).
Para cada potencial geotérmico o entalpía (alta- media o baja) hay asociada una tipología de
sistema de producción/generación de energía eléctrica o térmica. La elección de un sistema u otro
depende de la temperatura, figura 2.3, y el gradiente geotérmico.
Figura 2.3 Recursos y tecnologías geotérmicas en función de la temperatura. Figura 2.4 Variación de la tempera-
(Utilization and tourism valorisation of geothermal waters in Croatia tura en función de la profundidad
Stasa Borovic, Izidora Markovi). (Unión geotérmica Italia).
10
Los principales sistemas de aprovechamiento de la energía geotérmica se dividen en sistemas
con disposición vertical y horizontal, que pueden ser de circuito cerrado o de circuito abierto, y con
pozos de extracción de fluido y/o pozos de aporte.
Los sistemas con disposición vertical son los más usados, porque se adaptan a todos los tipos de
instalaciones (alta, media y baja entalpía). Esta tipología puede llegar a profundidades del orden de
kilómetros (geotermia de alta profundidad) y ser de tipo cerrado y abierto, figura 2.5. Si hablamos de
alta y media entalpía para la producción de electricidad, los sistemas vienen asociados a una planta de
generación de energía con ciclo Rankine, compuestos por uno o más intercambiadores de calor, turbinas
y generadores de electricidad. En algunas aplicaciones, se aprovecha la temperatura del fluido para
alimentar sistemas mixtos de generación de electricidad y calor “district heating”.
Los sistemas con disposición horizontal trabajan, en la mayoría de los casos, con circuitos
cerrados, intercambiadores y bombas de calor, y se usan en geotermia de baja entalpía.
Figura 2.5 Sistemas con disposición vertical cerrada y abierta (Geothermal heat pump systems: Status review and
comparison with other heating options Stuart J. Self, Bale V. Reddy, Marc A. Rosen).
Q = Q* + DA
𝑇 = 𝑇surface + DΓ
donde el flujo de calor Q y la conductividad térmica de las rocas K están relacionadas por la expresión
Q=𝐾Γ
11
Los principales lugares en el mundo con gradiente geotérmico más elevado son los lugares con
espesores pequeños del manto y/o donde existen discontinuidades en el mismo, como sitios con
presencia de fenómenos tectónicos (movimientos de fallas) y de discontinuidades. Así, conociendo la
temperatura necesaria para producir electricidad a partir del recurso geotérmico (aproximativamente
150oC), podemos localizar áreas del mundo donde será más fácil y rentable aprovechar este recurso. Por
tanto, es de esperar que se pueda obtener una mayor producción de electricidad o de calor en aquellos
lugares donde el potencial sea más elevado, figura 2.6.
Figura 2.6 Localización de placas oceánicas y sitios de fácil aprovechamiento de energía geotérmica para la
producciónde electricidad (Chandrasekharam D., Bundschuh J., 2008. Low-enthalpy Geothermal Resources
for Power Generation. Balkema/CRC Press/Francis & Taylor, Boca Raton and London)
El país con mayor potencial instalado del mundo es los Estados Unidos, con una capacidad de
3440 MW sobre unos 12000 MW de producción mundial en 2015. Los principales países donde
encontramos capacidades de generación de energía a partir de recursos geotérmicos son Indonesia,
Italia (conocida por disponer de la primera instalación del mundo para generación eléctrica, Larderello,
1901, Florencia), Nueva Zelanda e Islandia, este último es el país donde se encuentra la instalación más
grande de Europa (Hellisheidarvirkjun power station, con capacidad para producir electricidad con una
potencia de unos 330 MW y, al mismo tiempo, agua caliente sanitaria con una capacidad térmica de
unos 404 MWth, con 50 pozos de 2200 metros de profundidad), figura 2.7.
Figura 2.7 Distribución mundial de recurso geotérmico. Capacidad instalada en el mundo entre 2007 y 2012.(Jennejohn D,
Hines B, Gawell K, Blodgett L. Geothermal: 2012 international market overview report. Washington, D.C. USA: Geothermal
Energy Association; 2012).
12
2.3. Principales factores a considerar para el diseño de una planta geotérmica
Cuando se aborda cualquier problema científico y de diseño de obras, después haber definido el
problema y las principales variables en juego, se intenta predecir y predeterminar la respuesta del
sistema, usando herramientas de cálculo que faciliten y permitan determinar la solución más precisa
posible y de fácil interpretación. En el campo del cálculo numérico geotérmico, hay muchos programas
de cálculo convencionales y de investigación. En este apartado se explicará el problema del cálculo, no
solo desde un punto de vista de la respuesta térmica del suelo, sino también desde el punto de vista de
la definición de variables para el diseño final de la instalación. Para el diseño de la planta geotérmica
puede ser útil dividir el problema en cinco partes:
Figura 2.8 Estudio y evaluación del potencial geotérmico en Jigiri y Yokoke. (Integration on of Geothermal Exploration Data and
Numerical Simulation Data Using GIS in a Hot Spring Area, Satoru Higuchi, Jun Nishijima and Yasuhiro Fujimitsu).
13
2.4 Principales factores a considerar en el problema del cálculo numérico en
geotermia
No vamos a entrar en detalles en este trabajo acerca de la implementación de programas GIS,
que son el primer paso a realizar para elegir el tipo de instalación. En segundo lugar, la caracterización
del terreno se hace mediante la implementación de ecuaciones de conservación de la masa, del
momento y de la energía, dependiendo de la tipología de los depósitos. Uno de los factores más
relevantes en la definición del problema es la composición químico-física del subsuelo de análisis. En
este sentido, podemos dividir el problema en dos tipos.
𝜕 (𝜀 𝜌)
+ ∇ . 𝜌𝑢 = 0
𝜕𝑡
𝑢
𝜕 (𝜌 ) 𝑢 𝑢 𝑢 𝜇
𝜕𝑡
𝜀
+ ∇ . 𝜌 𝜀 . 𝜀 = − ∇𝑃 . ∇ . 𝜇 𝑒𝑞𝑣 ∇ 𝜀 - 𝜇 𝐾 𝑢 + ρg =0
𝜕 (𝜀 𝜌𝐶𝑝 f T f )
𝜕𝑡
+ 𝑢 . ∇ . (ρCp )f = ∇. Tf + K eff
f + ha (Tf + Ts)
𝜕 (𝜀 (1 − 𝜀 𝜌 𝐶𝑝 )𝑠 𝑇𝑠 ) 𝑒𝑓𝑓
= ∇. 𝐾𝑠 . ∇𝑇𝑠 + 𝑎 𝑇𝑓 + 𝑇𝑠
𝜕𝑡
𝜌 = 𝜌0 − Δ 𝜌Θ
con densidad lineal en un intervalo de temperatura dado por T0 = 20 °C y T1 = 250 °C, y con temperatura
normalizada,
𝑇− 𝑇0
Θ=
𝑇1 − 𝑇0
14
La variación de densidad puede depender de factores como la presión y la viscosidad dinámica.
Por ejemplo, Knabner P y Frolkovic P [32], o Intraval y Graf [37], proporcionan las siguientes expresiones
en las que la presión influye en la densidad por debajo del 3% hasta una profundidad de 600 m,
𝜌 1+ 𝜌 0
P = P0 – g (𝜌0 𝑧 + 2
) 𝑧2
𝜌 = 𝜌𝑜 e (4,5.10exp-10 (p – p0 ))
Otra característica que puede influir sobre la repuesta del dominio en términos de variación de
densidad son la viscosidad y la conductividad hidráulica, figura 2.9. La variación de la viscosidad
dinámica con la temperatura, de acuerdo con Holzbecher y Molson [36], viene dada por,
Figura 2.9 Variación de la viscosidad dinámica [mPa] en función de la temperatura (de Fórmulas Modeling
Density-Driven Flow in Porous Media, Ekkehard O. Holzbecher) y variación de la viscosidad dinámica con la
densidad (de Simulation of Geothermal Flow in Deep sedimentary basins in Alberta).
En los próximos capítulos se definirán los procesos básicos y las formulaciones numéricas para
resolver el problema de cálculo.
15
2.5 Estado del arte del cálculo numérico en geotermia
El cálculo numérico para la predicción de la respuesta de un depósito geotérmico es una de las
principales operaciones para determinar las variables útiles para el dimensionamiento de instalaciones
geotérmicas. Existen muchos programas de cálculo numérico, comerciales y de investigación útiles para
determinar la respuesta del dominio geotérmico, que se diferencian por las propiedades de cada
sistema y de las variables a determinar. Sin embargo, podemos definir dos problemas básicos, ya
descritos previamente: el primero corresponde a la trasmisión de calor en un sistema impermeable o
con baja permeabilidad, en el cual el fenómeno dominante es la conducción, y el segundo a la
transmisión de calor en un medio poroso, que viene regulado por el principio de la diferencia de
densidad (Darcy) y variación de presión con transporte de fluido [38].
16
2.5.1 Método de los elementos finitos (o diferencias finitas)
El método de los elementos finitos puede resumirse en cinco pasos básicos [38]:
En presencia de agua subterránea, la ecuación de balance de flujo de calor de una masa sólida
puede ser escrita en la forma [39],
𝜕𝑇
𝜌𝑐 = ∇ 𝜆∇𝑇 − ∇ 𝜌𝑤 𝑐𝑤 𝑣𝑇 + 𝐻 (𝑥, 𝑦, 𝑧)
𝜕𝑡
El flujo de agua de densidad constante en tres dimensiones en un medio poroso que cumple la
ley de Darcy, puede ser descrito por la ecuación en derivadas parciales,
𝜕 𝜕 𝜕 𝜕 𝜕 𝜕 𝜕
𝐾𝑥𝑥 + 𝐾𝑦𝑦 + 𝐾𝑧𝑧 − 𝑤 = 𝑆𝑠
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑧 𝜕𝑧 𝜕𝑡
donde x, y, z son las coordenadas cartesianas ortonormales alineadas a lo largo de las direcciones
principales del tensor de conductividad hidráulica Kij, h el potencial hidráulico por unidad de peso
(z + p/ϒ + v2/2g expresado en metros), W es el término de fuente-volumen de agua por unidad de
tiempo inyectado/extraído por unidad de volumen de acuífero en el punto (x, y, z)-, Ss el coeficiente de
almacenamiento específico del medio, y t el tiempo.
HYDROGEOSPHERE (FEM)
Es un software numérico de elementos finitos con control de volumen 3D, desarrollado por la
Universidad de Waterloo y de Laval, en grado de solucionar sistemas en los regímenes de flujo
superficial y subterráneo, en condiciones saturadas y no-saturadas, de estratigrafías diferentes y
teniendo en cuenta problemas de transporte de materiales usando la técnica de formulaciones
de la densidad variable, resolviendo las ecuaciones de flujo, transporte de energía y sólidos
simultáneamente, proveyendo de un balance completo.
17
FEFLOW (FEM)
Es un paquete de Software que usa el método de elementos finitos para el análisis y resolución
de diversas aplicaciones físicas y de ingeniería, desarrollado por la Universidad Tecnológica Real
de Estocolmo, Suecia. Cuenta con diversos módulos específicos, entre ellos se tienen el de flujo
subsuperficial para modelación de acuíferos. Aplica las ecuaciones de Richard, la ley de Darcy y
las ecuaciones de Navier-Stokes para el análisis de flujos subsuperficiales así como para la
modelación de transporte, reacción de sólidos y transporte de calor en medios porosos.
Este software simula el flujo subterráneo en un medio poroso en una, dos o tres dimensiones,
fue desarrollado por el Servicio Geológico de los Estados Unidos (U.S.G.S.), es considerado como
un estándar internacional para la simulación y predicción de las condiciones del agua
subterránea así como las interacciones del agua superficial y subterránea, utiliza el método
numérico de diferencias finitas para resolver mediante diversas interacciones la ecuación de
flujo del agua subterránea, por medio de la interfaz gráfica Visual MODFLOW se proporciona el
ambiente más adecuado para modelar de forma más completa y fácil pues se pueden visualizar
las simulaciones bidimensionales y tridimensionales del flujo de agua subterráneo y el
transporte de contaminantes en cualquier momento durante el desarrollo del modelo.
18
2.5.2 Métodos de los volúmenes finitos
El “método de los volúmenes finitos” es una técnica numérica utilizada para resolver ecuaciones
diferenciales en derivadas parciales, principalmente las relacionadas con los fenómenos de transporte.
Una de las diferencias importantes del método de los volúmenes finitos comparado con el método de
las diferencias finitas, es que puede ser implementado en una malla estructurada o no estructurada.
También, a diferencia del método de los elementos finitos, las condiciones de contorno en el método de
los volúmenes finitos son no invasivas, dando más estabilidad al proceso numérico [40].
TOUGH 2 (FEV)
TOUGH 2 (Pruess, 1991), basado en versión anterior MULKOM (Pruess, 1983), es un simulador
de yacimientos de propósito general capaz de simular flujos no isotérmicos de múltiples
componentes, fluidos de fases múltiples en medios porosos y fracturados. Su formulación de
volúmenes finitos le da la flexibilidad para modelar en 2D o 3D los sistemas geotérmicos [12]
teniendo en cuenta, con el modulo EWASG (Equacion of state for Wáter, Salt and Gas), la
presencia simultánea de muchas componentes.
19
2.5.3 Método de simulación por redes eléctricas MESIR
Otra metodología para resolver el problema de depósitos geotérmicos y cálculo de las variables
puede ser la utilización de programas que utilizan la metodología MESIR. En este trabajo todas las
simulaciones de depósitos geotérmicos serán hechas con el programa FAHET. El método de simulación
por redes es una técnica para el estudio de cualquier proceso físico que pueda definirse mediante un
modelo matemático o conjunto completo de ecuaciones, generalmente en derivadas parciales
acopladas o no. Partiendo de estas el procedimiento consiste en dos etapas bien diferenciadas: en
primer lugar, elaborar un “modelo en red” o circuito eléctrico equivalente al proceso, y en segundo
lugar, simular dicho proceso obteniendo la solución del modelo en red mediante un programa adecuado
de resolución de circuitos eléctricos [41].
FAHET (MESIR)
FAHET hace uso de la analogía o equivalencia entre el transporte eléctrico y el de fluido con
difusión de calor en un medio poroso y se presenta al usuario a través de un entorno de
comunicación ameno, tipo ventanas, que dirige paso a paso las acciones y opciones posibles
tales como selección y definición del problema, entrada de datos, creación y manipulación de
archivos de modelos, opciones de simulación, simulación avanzada, presentación de resultados,
etc. Los archivos de modelos en red se ejecutan en PSpice (1994), y los resultados de simulación
se ofrecen directamente en el entorno, gráfico o tabulado de salida de PSpice o bien, mediante
manipulaciones adecuadas, en el entorno gráfico del propio programa y, en mayor detalle, en el
entorno gráfico del software MATLAB (1997), merced a rutinas auxiliares incorporadas al
programa. Asimismo, FAHET permite presentar soluciones animadas de los isolíneas de
concentración y flujo (función de corriente) en problemas transitorios [41].
20
CAPÍTULO 3. FUNDAMENTOS TEÓRICOS QUE RIGEN LOS PROCESOS DE
FLUJO DE FLUIDOS EN EL TRANSPORTE DE CALOR
3.1 Introducción al problema
Si hay un gradiente térmico en un sistema, o si dos o más sistemas con temperaturas diferentes
están puestos a contacto, hay transmisión de energía. Este proceso, mediante el cual hay intercambio
de energía entre sistemas, se llama transmisión de calor. En éste, el calor no puede ser medido, pero
pueden ser medidos los efectos que se van a producir sobre el sistema, como variación de su energía
interna. Los principios básicos, que explican el calor y la energía interna de un sistema, están definidos
en termodinámica mediante el primer principio, que habla de la conservación y transformación de la
energía, y con el segundo principio, que implica que no existe ninguna máquina que trabaje por ciclos
pueda producir trabajo si solo se toma calor desde una región a temperatura mayor y no se cede parte
del mismo a otra a temperatura menor. Cuando hablamos de estos procesos, no podemos pararnos a
estudiar solo la Termodinámica, porque en todos los casos, la transformación real se producirá en
ausencia de equilibrio. Por esta razón, el estudio debe basarse en otras ramas de la ciencia, que van a
traer en cuenta procesos como transporte de materia, variación de características químico-físicas de los
materiales, etc. [42].
𝜕
S 𝜕𝑡 + ∇ . 𝑞 = 0
𝜌0 𝑔 ρ + ρ0
𝑞= − 𝑘 ∇h+ ∇z
𝜇 ρ0
donde
La ley que describe el movimiento de flujo y transporte de calor se ajusta bien con la ley de
Darcy (para el flujo en cada región). El flujo regional o local de agua es un flujo de tipo advectivo o
convectivo (arrastre), incluido por las condiciones de contorno en el cual hay una altura piezométrica
constante, flujo constante, recarga o extracción (con bombeo) y algunas veces es considerado como de
convección forzada según los autores Holzbecher y Yusa (1995). En la formulación de Darcy la variable
21
de la cual depende la variación de densidad del agua es la temperatura. Este fenómeno es el primero del
cual depende el flujo de fluido: si la densidad del fluido de la parte superior es mayor que la dela parte
inferior, tenderá a descender, creando lugares o porciones del medio, donde tendríamos una
recirculación continua de fluido. La formulación introducida antes, puede ser escrita o, en términos de
altura piezométrica h (Holzbecher [16]) o a través de la formulación tratada por Diersch [35], en
términos de presión,
p = 𝜌0 gh
𝑘 𝜌
𝑞= ∇ 𝑝 − 𝜌𝑔 = −𝑘 ∇ − 𝑔
𝜇 𝜌0
siendo:
- k el tensor de permeabilidad,
- μ la viscosidad dinámica,
- k la conductividad hidráulica efectiva.
𝑘𝜌0 𝑔
𝐾=
𝜇
La ley que gobierna el fenómeno de la trasmisión de calor en medios porosos es la ecuación de
conducción definida como:
𝜕𝑇 𝜕𝑇
∇𝑗. = − 𝜌 ∙ 𝑐𝑒 ∙ = − 𝜌𝑚 ∙ 𝑐𝑒,𝑚 ∙ = 𝑘𝑚 ∇2 𝑇
𝜕𝑡 𝜕𝑡
donde:
- j es la densidad de flujo de calor.
- 𝜌 𝑦 𝜌𝑚 son la densidad del agua y la densidad referida a un medio poroso formado por dos
componentes, respectivamente.
- ce y ce,m son el calor específico del agua y el correspondiente al medio poroso, respectivamente.
- 𝑇es la temperatura.
- t es el tiempo.
- km es la conductividad media del medio.
Las principales ecuaciones que gobiernan el fenómeno de transporte de flujo de fluido y calor
son el teorema de transporte de Reynolds y la ecuación del flujo de fluido combinada con la ley de
Darcy.
𝑑 𝑑𝜏
𝜏𝑑𝑅 = 𝑑𝑅 + 𝜏 𝑉𝜏 𝑛𝑑𝑆
𝑑𝑡 𝑑𝑡
𝑄 𝑄 𝑆𝑢
22
𝑑𝜏
𝑑𝑅 + ∇ τVt 𝑑𝑅 = 0
𝑑𝑡
𝑅 𝑆
𝜕𝜏
+ ∇ τVt = 0
𝜕𝑡
𝜕𝜏
+ ∇ 𝜏 𝑉𝜏 + ∇ 𝜏 𝑉𝜏 − 𝑉 =0
𝜕𝑡
𝜕 𝜌𝜑
+ ∇ 𝜌𝜑𝑉 = 0
𝜕𝑡
𝜕 𝜌𝜑 𝑘
= ∇ 𝜌 (∇ 𝜌 − 𝜌𝑔)
𝜕𝑡 𝜇
Si no hay generación interna de calor y disipación térmica por efecto viscoso, el balance de calor
se escribe como:
𝜕𝑇
𝜌𝑓 𝑐𝑒,𝑓 = ∇ 𝜍 + 𝑉 ∇𝑇 = 𝑘𝑚 ∇2 𝑇
𝜕𝑡
donde f es el subíndice relativo al fluido, y m a las propiedades globales del medio. El término 𝜍, y km
(combinación de propiedad del fluido y del medio), permiten estudiar la relación del calor específico por
unidad de masa entre el medio y el fluido como
𝜑 𝐶𝑒 + 1 − 𝜑 𝜌𝑆 𝐶𝑒,𝑠 𝜌𝑚 𝐶𝑒,𝑚
𝜍= =
𝜌𝑓 𝐶𝑒,𝑓 𝜌𝑓 𝐶𝑒,𝑓
𝑘𝑚 = 𝜑𝑘𝑓 1 − 𝜑 𝑘𝑠
En la mayoría de los casos, donde se estudian problemas geotérmicos, la transmisión del calor
viene descrita a través de la formulación de la función de corriente (𝜓, en m2/s) para flujo 2D, en la cual,
usando la aproximación de Boussinesq,
𝜕𝜓 𝜕𝜓
= 𝑞𝑧 = − 𝑞𝑥
𝜕𝑥 𝜕𝑧
23
La ecuación que describe el flujo de fluido, junto con la ley de Darcy para el caso más general,
conduce a la ecuación de flujo:
𝜕 𝜇 𝜕𝜓 𝜕 𝜇 𝜕𝜓 𝜕 𝜇 𝜕 𝜇 𝜕𝜌
+ = 𝑞𝑧 + 𝑞𝑥 = −𝑔
𝜕𝑧 𝑘𝑧 𝜕𝑥 𝜕𝑥 𝑘𝑥 𝜕𝑧 𝜕𝑥 𝑘𝑧 𝜕𝑧 𝑘𝑥 𝜕𝑥
En los problemas que hemos abordado hasta ahora, se han encontrado formulaciones con
variables dimensionales y adimensionales. Es conveniente uniformizar las unidades de medida y tratar el
problema de manera adimensional. Para hacer esto, podemos adimensionalizar todas las variables de
cálculo siguiendo el método propuesto por Henry o el propuesto por Holzbencher [16], Tabla 3.1 y 3.2:
x x’ 𝑦 𝑦
Anchura 𝑙𝑑𝑜 𝐻
𝑇 𝑇′ 𝑇 𝑇 − 𝑇𝑚𝑖𝑛
Temperatura 𝑇0 𝑇𝑚𝑎𝑥 − 𝑇𝑚𝑖𝑛
q q’ 𝐷 𝐷
𝑞 𝑞
Velocidad Darcy 𝑄 𝑄
t t’ 𝐷
t 𝐷
Tiempo 𝑑2 𝑡
𝐻2
V V’ 𝐻
𝑉
Velocidad flujo 𝐷
𝜓 𝜓′ 𝜓
Función de corriente 𝐷
De este modo, las fórmulas que describen el fenómeno y las funciones de corriente se convierten en
Tipo de adimensionalización
Formulación
𝜕 2 𝜓′ 𝜕 2 𝜓′ 1𝜕𝑇′
Ec. De flujo + =
𝜕𝑥 ′ 2 𝜕𝑥 ′ 2 𝑎 𝜕𝑥′
Henry 𝜕 2 𝑇′ 𝜕 2 𝑇′ 1 𝜕𝜓′ 𝜕𝑇′ 𝜕𝜓′ 𝜕𝑇′ 𝜕𝑇′
Ec. De transporte ′ 2
+ ′ 2
− − =
𝜕𝑥 𝜕𝑥 𝑏 𝜕𝑧 𝜕𝑥 𝜕𝑥′ 𝜕𝑧′ 𝜕𝑡′
𝜕 2 𝜓′ 𝜕 2 𝜓′ 𝜕𝑐′
Ec. De flujo + = 𝑅𝑎
𝜕𝑥 ′ 2 𝜕𝑥 ′ 2 𝜕𝑥′
Holzbencher
𝜕 2 𝑇′ 𝜕 2 𝑇′ 1 𝜕𝜓′ 𝜕𝑇′ 𝜕𝜓′ 𝜕𝑇′ 𝜕𝑇′
Ec. De transporte ′ 2
+ ′ 2
− − =
𝜕𝑥 𝜕𝑥 𝑏 𝜕𝑧 𝜕𝑥 𝜕𝑥′ 𝜕𝑧′ 𝜕𝑡′
Tabla 3.2. Formas adimensionales de las ecuaciones de flujo y de transporte
24
Ra es el número de Rayleigh, siendo
𝑘Δ𝜌𝑔𝐻
𝑅𝑎 =
𝐷𝜇
Para terminar se estudian las condiciones de frontera, de manera que cada variable
dependiente afín al problema, pueda ser determinada matemáticamente y tenga una sola solución. En
la Tabla 3.3 se especifican las condiciones de frontera de primera y segunda clase:
Condiciones de contorno
Tipo Variable Flujo Variable Temperatura
De primera clase Se especifica el valor de la función de Se especifica el valor de la
corriente temperatura
(condición de Dirichlet)
De segunda clase Se especifican los valores de las Se especifica el valor del flujo de
velocidades paralela y perpendicular calor en la frontera
(condición de Neumann) a la frontera
Tabla 3.3 Condiciones de contorno y su significado físico
25
En las Tablas 3.4.a y 3.4.b tenemos las expresiones de los tipos de condiciones de contorno más
comunes.
26
CAPÍTULO 4. MÉTODOS DE SIMULACIÓN NUMÉRICA EN PROBLEMAS DE
GEOTERMIA
4.1 Introducción
En este capítulo, vienen descritos los principales métodos para solucionar los problemas en el
campo del cálculo numérico en geotermia. Aquí serán descritos los métodos de simulación por redes
eléctricas, (utilizado en este trabajo, ver capítulo 6) y a continuación el método de resolución de las
diferencias finitas. Una descripción detallada del método, con numerosas aplicaciones a los campos de
transporte a través de membranas, transferencia de calor, sistemas de reacción química, transferencia
de masa en soluciones electrolíticas y no electrolíticas y fenómenos electrocinéticos en suspensiones
coloidales, puede encontrarse en el libro de González- Fernández (2002). En los últimos años el método
ha sido aplicado con éxito en otros campos de investigación, tales como transporte de calor en fluidos,
sistemas caóticos, vibraciones mecánicas, elasticidad, problemas inversos, etc., incluyendo el campo que
nos ocupa de flujo asociado a densidad variable con transporte de soluto, así como explican Moreno y
col. (2007), Alhama y col. (2003 a-b, 2004), Zueco y Alhama (2006, 2007), Zueco y col. (2004, 2005, 2006
a-b) y Soto y col. (2007 b-d), así como en el diseño de programas educativos Lopera y col. (2009) y
Alhama y Del Cerro (2010 a-b).
27
4.2 Método MESIR
En este primer apartado viene estudiada la técnica utilizada, para solucionar los problemas
geotérmicos usando el software PSpice y el programa FAHET, a través del método MESIR. El método de
Simulación por redes eléctricas permite el estudio de problemas físicos, en este caso de transmisión de
calor en medios porosos, mediante un modelo matemático o conjunto completo de ecuaciones,
generalmente en derivadas parciales acopladas o no. Partiendo de estas el procedimiento consiste en
dos etapas bien diferenciadas: en primer lugar, se elabora un “modelo en red” o circuito eléctrico
equivalente (en FAHET), y en segundo lugar, se simula dicho proceso, obteniendo la solución del modelo
en red mediante un programa (PSpice), adecuado para la resolución de circuitos eléctricos. La
equivalencia formal, entre el modelo en red y el proceso físico en estudio, reside en que ambos se rigen
por las mismas ecuaciones discretizadas en el espacio, es decir por las mismas ecuaciones diferenciales
en diferencias finitas, referidas a un elemento de volumen o celda y las mismas ecuaciones discretizadas
para las condiciones de contorno [41].
En esta técnica, viene establecida una correspondencia entre variables flujo por un lado
(densidad de corriente eléctrica con flujo de calor, flujo de masa…) y variables tipo potencial por otro
(potencial eléctrico con temperatura, concentración…) pero es posible establecer otras analogías aún en
procesos físicos que describan el transporte de una determinada magnitud. Por ello, el que los procesos
de flujo debido a densidad y transporte de soluto admitan redes eléctricas equivalentes, supone no
sólo la equivalencia matemática sino, también, la equivalencia física entre las variables características de
unos y otros procesos. El MESIR es, por otro lado, un método de simulación en tanto que incluye la
resolución numérica del modelo en red. Así, las variables flujo y fuerza características del mismo deben
satisfacer las leyes de Kirchhoff, y sus relaciones constitutivas determinarán los elementos de circuito
correspondientes.
Ahora bien, en cada proceso concreto y una vez elegidas las variables conjugadas, la información
de qué elementos de circuito intervienen en el modelo en red y cómo se conectan entre sí, se obtiene
del modelo matemático y no de consideraciones de tipo físico acerca del papel que juegan estas
variables.
28
En síntesis, en la teoría de redes, la viabilidad de un modelo en red supone:
Las relaciones entre flujo y fuerza asociados a una rama y sus (dos) nudos límite, que pueden
incluir o no variaciones temporales de estas variables que se dicen conjugadas, definen los elementos
concretos del circuito equivalente a esa rama. La relación causa-efecto entre las variables conjugadas es
completamente arbitraria con tal que sea consistente con ii) y iii).
MONOPUERTAS BÁSICAS
A la red se le asocia un conjunto de flujos que obedecen a una ley de balance local y un conjunto
de fuerzas que satisfacen la condición de unicidad. Tales requisitos dan cuenta de la topología de la red
relativa al proceso. Las propiedades topológicas dependen únicamente de la asignación de conexiones
entre los diferentes puntos o de las posibles combinaciones de trayectorias que unen un nudo dado con
otros nudos. Son independientes de las medidas y, desde un punto de vista topológico, dos grafos son
iguales o isomorfos si las asignaciones de vértices y ramas son las mismas.
Las leyes de Kirchhoff establecen relaciones entre flujos y fuerzas por separado, pero no
expresan ningún tipo de relación entre flujos y fuerzas entre sí. Las relaciones entre el par conjugado
flujo-fuerza se conocen como ecuaciones constitutivas o fenomenológicas y definen los elementos de
circuito que expresan características específicas de cada proceso. Se dice que dos grafos son
geométricamente iguales si los potenciales y flujos de cada par de puntos y su rama correspondiente
son iguales para cualquier conjunto de valores que puedan ser elegidos para los flujos o las fuerzas. Las
propiedades geométricas de la red, es decir, sus características métricas, se siguen de las relaciones
constitutivas, que se pueden establecer entre las variables de un par flujo-fuerza, en cuyo caso se habla
de monopuerta.
La construcción del modelo eléctrico equivalente viene hecho a través de elementos como
monopuertas pasivas (resistiva, capacitiva, inercial o inductiva) y activas (constantes, dependiente del
tiempo y controladas). La implementación de cada una de estas permite construir el sistema eléctrico
equivalente y solucionar el conjunto de ecuaciones. Las leyes de Kirchhoff establecen relaciones entre
flujos y fuerzas por separado, pero no expresan ningún tipo de relación entre flujos y fuerzas entre sí.
Una primera clasificación está relacionada con lo que en electricidad se conoce como elementos
pasivos y activos. Los elementos pasivos no generan potencia sino que la disipan (transformación
energética), tienen la capacidad de almacenarla y/o entregarla a la red y constituyen lo que llamamos
monopuertas pasivas.
Las fuentes de tensión y corriente son elementos activos, generan potencia de acuerdo a una
determinada ley y son las monopuertas activas o fuentes (es posible, no obstante, que una relación
constitutiva correspondiente a una monopuerta pasiva pueda ser representada mediante una
monopuerta activa donde la función de control es una constante).
29
Monopuertas pasivas. En función de la relación expresa existente entre las variables LCK y LVK,
las monopuertas pasivas tienen nombres específicos:
-La monopuerta resistiva es un elemento de circuito asociado a una relación entre las derivadas
temporales de las variables flujo y fuerza de una misma rama, mediante una función independiente del
tiempo que llamaremos resistencia R, que puede depender o no del flujo o de la fuerza:
A partir de esta expresión es posible relacionar las variables en forma finita y escribir
𝑋 𝑡 = 𝐹𝑅 (𝐽)
Una monopuerta resistiva es lineal cuando la relación entre las variables X(t) y J(t) lo es, es decir
X(t)= RJ(t); naturalmente R es una constante en este caso. Su acción es instantánea, no importa cual sea
su estado anterior: en este sentido carecen de memoria. En su analogía física representan efectos
disipativos, fricciones, efectos viscosos, energías de reacción, etc., y desde el punto de vista
termodinámico son elementos generadores de entropía.
Las monopuertas resistivas no lineales se definen a través de las funciones que las
caracterizan,𝐽 𝑡 = 𝐹𝑅−1 (𝑋)o𝑋 𝑡 = 𝐹𝑅 𝐽 . Constituyen, en definitiva fuentes controladas de corriente
o tensión, respectivamente. En la figura 4.1 se muestra la representación simbólica de la monopuerta
resistiva lineal y no lineal.
Monopuerta Resistiva
Tipo Lineal No lineal No lineal
Componente
eléctrico
usado
Relaciones 𝑋 𝑡 = 𝑅𝐽(𝑡) 𝑗 𝑡 = 𝐹𝑅 (𝐽) 𝑉 𝑡 = 𝐹𝑅−1 (𝑗0)
Figura 4.1 Representación simbólica de una monopuerta resistiva lineal y no lineal.
-La monopuerta capacitiva es un elemento de circuito asociado a una relación entre la variable
flujo y la derivada temporal de la variable fuerza, de una misma rama, mediante una función no
dependiente del tiempo que designaremos como capacidad C. En ésta se produce almacenamiento sin
pérdidas (no hay disipación energética), y su estado, que no cambia instantáneamente, tiene en cuenta
todas las operaciones llevadas a cabo en el pasado (tiene memoria). En su analogía, representa procesos
físicos en los que se produce algún tipo de almacenamiento como condensadores, tanques, etc.
𝑑𝑋(𝑡)
𝐽 𝑡 =𝐶
𝑑𝑡
𝑑𝑞 𝐹𝐶 (𝑋)
𝐶= =𝑑
𝑑𝑋 𝑑𝑋
30
Así como para el caso anterior cuando la dependencia 𝑞 = 𝐹𝐶 (𝑋) es lineal, 𝐶 = 𝑞/𝑋. Las
dependencias 𝑞 = 𝐹𝐶 (𝑋)no lineales, ejemplos de las cuales se presentarán en esta memoria, deben
estudiarse en cada caso. En la figura 4.2 se muestra la representación simbólica de la monopuerta
capacitiva lineal. La traducción al modelo en red es un condensador eléctrico de valor C faradios.
𝑑𝑗 (𝑡)
𝑋 𝑡 =𝐿
𝑑𝑡
que equivale a la relación, no dependiente del tiempo, entre las variables flujo y momento.
𝑝 = 𝐹𝐿 (𝐽)
Ésta es constante cuando la dependencia 𝑝= FL(J) es lineal, L= 𝑝/J. Al igual que en la monopuerta
capacitiva, las dependencias 𝑝=FL(J) no lineales deben estudiarse particularmente en cada caso.
El que cada elemento pueda ser caracterizado por un par de variables conjugadas con única
ecuación constitutiva entre ellas es una hipótesis básica en el MESIR que deriva de la teoría de redes.
Físicamente equivale a decir que es posible elegir un elemento de volumen lo suficientemente pequeño
como para que su tiempo de relajación interna sea mucho menor que el del sistema global, pero
suficientemente grande como para que las fluctuaciones de las variables que describe el sistema en él
sean despreciables.
31
Monopuertas activas. En éstas se produce una aportación o extracción de energía al sistema,
podemos distinguir:
-Fuentes dependientes del tiempo: la relación constitutiva entre las variables tiene la misma
forma de las fuentes constantes; además, 𝑋 = 𝑋 (𝑡) y 𝐽 = 𝐽 (𝑡) según se trate de fuentes de fuerza o de
flujo. En la figura 4.4.b se muestra la representación simbólica de fuentes dependientes del tiempo.
𝑋 = 𝐹𝑋 𝑋𝑐
𝑋 = 𝐹𝐽 𝐽𝑐
𝐽 = 𝐹𝐽 𝐽𝑐
𝐽 = 𝐹𝑋 𝑋𝑐
según se trate de i) fuentes de tensión controladas por tensión, ii) de tensión controladas por corriente,
iii) de corriente controladas por corriente y iv) de corriente controladas por tensión, respectivamente; F
designa una función arbitraria de la variable o variables de control.
La acción de control puede ser ejercida por más de una variable y las funciones de control
pueden ser complejas. Aunque la monopuerta puede especificarse arbitrariamente, su implementación
como elemento de circuito puede no ser posible en tanto que no esté contenida en las librerías del
software elegido. Sin embargo, la teoría de circuitos permite, mediante circuitos auxiliares, resolver
prácticamente todos los casos de diseño de la red eléctrica que se necesiten para cualquier tipo
complejo de fuente controlada. En la figura 4.4.c se muestra la representación simbólica de fuentes
controladas por una variable y en la 4.4.d fuentes controladas por varias variables.
El potencial de estas monopuertas activas para establecer los modelos en red de sistemas
fuertemente no lineales es inmenso ya que su uso permite imponer a la monopuerta el valor de una
variable (en función de variables de otras monopuertas) sin influir en la otra variable, cuyo valor, se
ajusta a la topología y geometría del modelo en red.
32
Figura 4.4 Representación simbólica de monopuertas activas. a) fuentes constantes, b) fuentes dependientes del tiempo, c)
fuentes controladas por una variable, d) fuente controlada por varias variables.
Elementos lineales:
Resistencias (R)
Condensadores (C)
Elementos no lineales:
33
4.3 Método de las diferencias finitas
Entre los métodos matemáticos más utilizados para solucionar problemas de geotermia, se
encuentra el método de las diferencias finitas. Recordando los resultados presentados en el capítulo
anterior, donde se incluye la ecuación estacionaria del movimiento de un fluido de densidad 𝜌𝑓 , y
permeabilidad k, podemos escribir la ecuación general en función de la altura piezométrica [45],
𝜌𝑓 𝑔 𝑘 𝜕 𝜕 𝜕 𝜕
∇= ∇0 + 𝜌𝑟 ∇𝑧 =0 → 𝑘𝑥 + 𝑘𝑧 ] =0
𝜇 𝜕𝑥 𝜕𝑥 𝜕𝑧 𝜕𝑧 + 𝜌𝑟
𝑣2 𝑃
= + 𝑧+
2𝑔 𝑔𝜌
𝑘 𝑑𝑝
𝑣=
𝜇 𝑑𝑥
La ecuación estacionaria del fluido puede escribirse como la suma de dos partes: un término
asociado a la difusividad y otro asociado con el transporte de fluido. La ecuación estacionaria del calor,
para un suelo con la conductividad 𝜆, puede escribirse como:
𝜕 𝜕𝑇 𝜕 𝜕𝑇
∇ 𝜆∇𝑇 − 𝜌𝑓 𝑐𝑓 𝑣𝑇 = 0 → 𝜆 − 𝜌𝑓 𝑐𝑓 𝑣𝑥 𝑇 + 𝜆 − 𝜌𝑓 𝑐𝑓 𝑣𝑧 𝑇 = 0
𝜕𝑥 𝜕𝑥 𝜕𝑧 𝜕𝑧
La conductividad térmica del suelo, puede ser calculada como suma de la del medio sólido y la
del fluido, teniendo en cuenta la porosidad de la tierra,
𝜆 𝜙, 𝑇, 𝑃 = 1 − 𝜙 𝜆𝑟 𝑇 + 𝜙𝜆𝑓 𝑇, 𝑃
donde la conductividad térmica de las rocas 𝜆𝑟 y su calor específico cr son funciones de la temperatura,
y se pueden calcular como:
770
𝜆𝑟 𝑇 > 𝑇𝐶 = + 0,7
350 + 𝑇
El precedente sistema de ecuaciones acopladas no puede ser resuelto con métodos analíticos,
por lo que es necesario discretizarlo mediante diferencias finitas. El método consiste en dividir el
dominio en una malla de cálculo, como representamos en la figura 4.5, solucionando cada celda, en la
cual aproximamos las derivadas parciales con diferencias finitas, obteniendo así un sistema de
ecuaciones lineales que pueden ser resueltas numéricamente.
34
Figura 4.5 Ejemplo de implementación de la malla de cálculo y los nodos de las celdas.
Una vez desarrollado el sistema de cálculo (celdas y nodos), se escriben las ecuaciones
dinámicas de la altura piezométrica y las ecuaciones térmicas como sistemas de ecuaciones
discretizadas,
Para solucionar la ecuación de arriba, las constantes se definen del siguiente modo:
𝑘𝑥 𝑘𝑥 𝑘𝑧 𝑘𝑧 𝑘𝑥 𝑘𝑧
𝐴= ; 𝐵= ;𝐶 = ;𝐷 = 𝜌𝑟 ; 𝐸 = ; 𝐹= 𝜌𝑟 ;
∆𝑥 2 ∆𝑥 2 ∆𝑧 2 ∆𝑧 2 ∆𝑧 2 ∆𝑧 2
𝐺 =𝐴+𝐵+𝐶+𝐸
𝛼 = 0; 𝑠𝑖 𝑣 > 0;
𝑇𝑖 +,𝑘 = 1 − 𝛼 𝑇𝑖,𝑘 + 𝛼 𝑇𝑖+1,𝑘 ;
𝛼 = 1; 𝑠𝑖 𝑣 < 0;
35
CAPÍTULO 5. EL PROGRAMA FAHET
FAHET hace uso de la analogía o equivalencia entre el transporte eléctrico y el de fluido con
difusión de calor en un medio poroso y se presenta al usuario a través de un entorno de comunicación
ameno, tipo ventanas, que dirige paso a paso las acciones y opciones posibles tales como selección y
definición del problema, entrada de datos, creación y manipulación de archivos de modelos, opciones
de simulación, simulación avanzada, presentación de resultados, etc. Los archivos de modelos en red se
ejecutan en PSpice (1994) y los resultados de simulación se ofrecen directamente en el entorno, gráfico
o tabulado de salida de PSpice o bien, mediante manipulaciones adecuadas en el entorno gráfico del
propio programa y (en mayor detalle) en el entorno gráfico del software MATLAB (1997), merced a
rutinas auxiliares incorporadas al programa. Asimismo, FAHET permite presentar soluciones animadas
de las isolíneas de concentración y flujo (función de corriente) en problemas transitorios. La figura 5.2
muestra un esquema del funcionamiento básico del programa. Su puesta en marcha da acceso,
directamente, a la entrada de datos: geometría de la reticulación, características físicas del medio
poroso, fluido y soluto, condiciones de contorno etc. Una vez completada la especificación del problema
se puede crear un archivo de texto básico del modelo que permite su manipulación directa y su
modificación.
36
Figura 5.2 Diagrama de funcionamiento de FAHET.
La asignación de nombres a los elementos del modelo en red así como la asignación de nodos
sigue una regla sencilla e intuitiva y lógica, permitiendo al usuario localizar inmediatamente el elemento,
sección o punto del modelo del que se desea obtener información acerca del valor de las variables
dependientes, flujos, potenciales y temperaturas tanto en el circuito de flujo de fluido como en el de
temperaturas. FAHET incorpora un gráfico directo en el que se muestra la disposición de las celdas o
volúmenes de control mostrando la leyenda de los nudos centrales de cada celda, que sirve también
para identificar todos los componentes de la celda (resistencias y condensadores), sus propiedades y las
condiciones de contorno, si existen. Como opción añadida pueden modificarse y actualizarse datos o
información sobre la propia retícula, tanto para una celda como para un conjunto seleccionado de ellas.
Por último, un archivo de ayuda accesible desde cualquier paso del programa da información al
usuario de cómo resolver e interpretar las dificultades que surgen en la explotación del mismo.
37
5.2 Creación de archivos de modelos
La carpeta FAHET contiene las subcarpetas mostradas en la figura 5.3 entre las que se encuentra
la del ejecutable “FAHET.exe” que da acceso a la pantalla de entrada del programa, “Designer FAHET-A”
figura 5.4. Esta pantalla inicial permite cargar un modelo ya existe a través de la ruta adecuada, opción
“Load” figura 5.5, o iniciar la creación de un nuevo modelo, opción “New” figura 5.6.
Figura 5.3 Subcarpetas de la carpeta “FAHET” con acceso al ejecutable “FAHET.exe” (Manual FAHET).
Figura 5.5 Pantalla de ruta para cargar modelos ya creados (Manual FAHET)
38
Figura 5.6 Pantalla de creación de nuevo modelo (Manual FAHET).
Los primeros datos para la creación de un modelo nuevo son el número de celdas horizontales y
verticales, su tamaño y las propiedades del fluido y del medio, viscosidad, gravedad e intervalo de
cambio de la densidad del fluido en unidades internacionales. La opción “aceptar” da acceso a la
retícula, figura 5.7 que contiene la numeración de celdas y nombres de las capas, siguiendo una regla
simple: los dos primeros dígitos corresponden a la posición horizontal (empezando por la celda 01 a la
izquierda) y los dos últimos a la posición vertical (empezando a contar de abajo a arriba).
Al hacer clic con el botón izquierdo del ratón, con el cursor ubicado en una celda cualquiera,
aparece información de la misma en un cuadro ampliado: número de celda y condición de contorno,
figura 5.8 para una celda del interior y figura 5.9 para una celda de contorno.
39
Figura 5.8 Detalle de información de una celda del interior (Manual FAHET).
Figura 5.9 Detalle de la información de una celda del contorno (Manual FAHET).
Al hacer clic con el botón derecho del ratón, en una celda o un conjunto de ellas, figuras 5.10 y
5.11, respectivamente, se presenta un cuadro de diálogo que da acceso a la edición o reedición de las
propiedades de la celda o de sus condiciones de contorno. Al pulsar el botón de “Editar” tras seleccionar
un grupo de celdas, se presenta la pantalla de la figura 5.12 en la que se muestra el cuadro de
parámetros a asignar a dichas celdas. Si éstas son del contorno se pueden editar tanto las propiedades
como las condiciones de contorno, figura 5.13, en una pantalla auxiliar se puede seleccionar el tipo de
condición (temperatura o función de corriente), la posición en donde se implementa la condición (parte
superior, inferior, izquierda o derecha de la celda) y tipo de condición (constante o adiabática), figura
5.14. Además, se pueden seleccionar celdas o grupos de celdas para deshabilitarlas, esto es, para crear
huecos en el dominio. Basta para ello marcar el botón “Enable cells” una vez se ha seleccionado la celda
o celdas, figura 5.15.
40
Figura 5.10 Pantalla de edición o reedición de la información de una celda (Manual FAHET).
Figura 5.11 Pantalla de selección simultánea de celdas para editar (Manual FAHET).
Figura 5.12 Subpantalla para editar las celdas seleccionadas (Manual FAHET).
41
Figura 5.13 Subpantalla para editar condiciones de contorno de las celdas seleccionadas (Manual FAHET).
42
Figura 5.15 Botón para deshabilitación de celdas, creación de huecos, en el dominio (Manual FAHET).
Figura 5.16 Condiciones de contorno. Aviso de error en entrada de datos (Manual FAHET).
43
El botón “circuits” del cuadro de diálogo principal da acceso a la pantalla mostrada en la figura
5.17. Dicha pantalla permite mostrar el código (archivo de texto del modelo) correspondiente al circuito
de la capa previamente seleccionada (véase que el circuito conserva el mismo informe de la capa). Esta
opción permite modificar sobre la pantalla, si se quiere, los parámetros y componentes del circuito o
incluso modificar por completo el subcircuito del conjunto de celdas que componen cada una de las
capas ya definidas.
El último botón del menú principal “Other simulation data” presenta al usuario la pantalla de la
figura 5.18. En ésta se introducen dos tipos de datos. Por un lado, el intervalo de tiempos de la
simulación (tiempo inicial, final e intervalo de tiempo con el que se imprimen los datos tabulados de la
solución). Estos datos se incluyen en la sentencia “.TRAN” que define el transitorio. En este mismo
grupo de datos se introduce la tolerancia relativa de cálculo “RELTOL” y el número de dígitos con el que
se presentan los datos de archivo de salida. Por otro, los datos de simulación de un parámetro.
Figura 5.18 Cuadro de diálogo de introducción de tiempos de simulación y simulación de parámetro (Manual FAHET).
44
El programa permite ejecutar simultáneamente un conjunto de archivos cuya diferencia está en
el valor de uno de los parámetros que definen el problema. Para ello se introduce el valor inicial del
parámetro, el valor final y el intervalo que define el conjunto de valores discretos del parámetro. Esta
opción permite presentar en el entorno gráfico de salida de PSpice, simultáneamente, las soluciones de
una variables para el conjunto de valores definidos de un parámetro. Una vez completado el diseño del
modelo, la opción “Generar CIR” o “Generar y procesar”, figura 5.19, da paso a la pantalla de la figura
5.20 que presenta el archivo de texto del modelo sobre el que se pueden hacer comprobaciones y/o
manipulaciones. Terminadas éstas puede guardarse, opción “Guardar CIR” (siguiendo la ruta indicada),
o ejecutarse, opción “Simular”. En este último caso FAHET abre PSpice y procede a la simulación.
Figura 5.19 Botones “Generar CIR” o “Generar y procesar” para crear el archivo de texto del modelo (Manual FAHET).
Figura 5.20 Archivo.cir del modelo con opciones de “Guardar CIR” o “Simular” (Manual FAHET).
45
5.3 Criterios para la numeración de celdas, nodos y elementos del modelo
Como en FATSIM-A, FAHET genera automáticamente la numeración de las celdas siguiendo un
criterio lógico consistente en atribuir a cada una un conjunto de 4 dígitos, de los cuales los dos primeros
indican la posición horizontal de la celda y los dos siguientes la posición vertical. Los nodos y celdas
asociados al circuito de concentración añaden una “c” a la numeración mientras que los
correspondientes a la función de corriente añaden una “f”. Esta misma nomenclatura se sigue para
nombrar los subcircuitos. Se llama subcircuito al conjunto de componentes eléctricos del modelo que
implementan las ecuaciones en diferencias finitas aplicadas a la celda o elemento de volumen. En
general, basta definir en el modelo dos subcircuitos (uno para la concentración y otro para la función de
corriente) para cada capa o conjunto de celdas idénticas y luego conectar eléctricamente los
subcircuitos de temperatura y función de corriente a las celdas contiguas y a las condiciones de
contorno correspondientes. Cada subcircuito tiene una denominación propia (asociada a la celda que
implementa) para distinguirlo de los demás.
Con la numeración de la figura anterior es fácil identificar la posición relativa de cada punto del
mallado a partir del número de celdas que contiene y solicitar los datos de las variables dependientes
temperatura y función de corriente en los puntos requeridos una vez realizada la simulación. Esta
definición de nudos es muy útil cuando se trata de buscar los errores o fallos de programación del
archivo usando directamente los resultados de la simulación mostrados en el entorno de salida gráfico
de PSpice.
46
En relación con la denominación de elementos del modelo, estos se definen con una letra inicial
que los identifica (R para la resistencia; C para el condensador; V para un generador de tensión
constante o pila; I para un generador de corriente constante; E para un generador de tensión controlado
por tensión, etc.) seguida de los números correspondientes a la celda a la que pertenecen. Debido al
diseño simétrico de la celda, se añade “izq” o “der” y “sup” o “inf” a las resistencias para identificar su
posición relativa en la celda. Con esta identificación intuitiva el usuario puede encontrar fácilmente el
elemento correspondiente a la posición deseada para solicitar la información que contiene (flujos y
potenciales).
Finalmente, en relación con los elementos de contorno se sigue una regla lógica para
identificarlos. La condición adiabática se implementa mediante resistencias conectadas entre el nudo
correspondiente del contorno y masa mientras que las condiciones de concentración o función de
corriente constante y las de convección se implementan mediante generadores entre los mismos nudos.
47
5.4 Estructura de los archivos de texto de modelos
La estructura de estos archivos está dividida en los bloques siguientes, por este orden:
o Nombre del archivo ( con la opción de incluir, a continuación del nombre, una
descripción general del problema).
o Parámetros físicos, geométricos y de reticulación.
o Descripción de los subcircuitos (isoconcentración y función de corriente) de las celdas
correspondientes a cada capa donde se especifican los componentes de los mismos y la
denominación de sus nudos internos.
o Listado de interconexión de subcircuitos, especificando el nombre del subcircuito y la
numeración de nudos externos.
o Listado de elementos que implementan las condiciones de contorno de cada variable,
indicando el tipo de elemento, su valor y los nudos de conexión.
o Listado de variables a imprimir.
o Sentencias de opciones de simulación.
El archivo, pues, está encabezado por el nombre que lo identifica y, opcionalmente, una
descripción no limitada del problema a que se refiere. Ejemplo de encabezamiento:
La sección siguiente del archivo de texto está formada por un listado de las variables que usa.
Éstas se refieren a los parámetros físicos (conductividades hidráulicas, permeabilidades, viscosidades,
conductividades térmicas, calores específicos y porosidad), parámetros geométricos del problema
(longitudes del medio), parámetros asociados a las condiciones de contorno (temperaturas de
referencia y valores de la función de corriente) y tamaño de las celdas de cada capa (ancho y alto). La
denominación de estas variables, que toman el valor dado en la especificación del problema o lo
deducen de los datos de entrada si están definidas mediante operaciones matemáticas, es una
abreviatura de su nombre completo con objeto de identificarlas fácilmente. El siguiente cuadro muestra
un ejemplo de listado de variables (las líneas que comienzan con asterisco “*” son comentarios de
aclaración para el usuario):
48
El tercer bloque define los subcircuitos. Cada par de ellos, asociados a la misma posición, se
corresponden con celdas de una misma capa: diferentes capas (de diferentes propiedades) tienen
distintos subcircuitos. La primera línea define el nombre del subcircuito seguido de la numeración de
nodos internos del mismo. En problemas 2-D, como los que nos ocupan, los nudos son 6 (izquierdo,
derecho, inferior, superior, central y masa, por este orden). Los componentes que contiene para el caso
general suelen ser cuatro resistencias correspondientes a los términos lineales de las EDP, dispuestas
simétricamente en la celda, un condensador conectado entre el centro de la celda y el nudo de
referencia para implementar el término de almacenamiento de las EDP, más los generadores
controlados de corriente para implementar los términos no lineales y/o acoplados. Los valores se
escriben directamente mediante números o mediante expresiones de las funciones de los parámetros
de los que dependen; para el caso de generadores controlados estas expresiones vienen encerradas
entre llaves. La sentencia “.ENDS nombre de la capa” cierra la especificación del subcircuito. Ejemplo de
descripción de subcircuitos del problema de Henry en el que un mismo subcircuito con 11 nudos
externos define tanto los componentes de la ecuación de transporte como los de la ecuación de flujo.
49
El siguiente bloque es el correspondiente a las condiciones de contorno de las celdas sometidas
a esta condición. En modelos 2-d los contornos son los bordes izquierdo, derecho, superior e inferior del
modelo más los bordes correspondientes a los huecos dentro del medio, caso de que existan. Los
componentes que implementan estas condiciones, siguiendo una numeración ordenada, se especifican
uno a uno. Dado que las líneas de programa tienen un número limitado de dígitos, es frecuente que el
listado se separa en bloques de variables ocupando un gran número de líneas. El siguiente ejemplo
corresponde a las condiciones del problema de Henry:
** **
* Superior * Inferior
R Superior 0c 0126xc 0 1E + 15 R Inferior 0c 0101xc 0 1E + 15
R Superior 1c 0226xc 0 1E + 15 R Inferior 1c 0201xc 0 1E + 15
R Superior 2c 0326xc 0 1E + 15 R Inferior 2c 0301xc 0 1E + 15
R Superior 3c 0426xc 0 1E + 15 ……
…… ……
R Superior 47c 4826xc 0 1E + 15 ……
R Superior 48c 4926xc 0 1E + 15 R Inferior 48c 4901xc 0 1E + 15
R Superior 49c 5026xc 0 1E + 15 R Inferior 49c 5001xc 0 1E + 15
** **
* Izquierda * Derecha
V Izquierda 0c 0101yc 00 V Derecha 0c 5101yc 01
V Izquierda 1c 0102yc 00 V Derecha 1c 5102yc 01
V Izquierda 2c 0103yc 00 V Derecha 2c 5103yc 01
…… ……
V Izquierda 22c 0123yc 00 ……
V Izquierda 23c 0124yc 00 V Derecha 23c 5124yc 01
V Izquierda 24c 0125yc 00 V Derecha 24c 5125yc 01
** **
* Superior * Inferior
V Superior 0s 0126xs 0 6.6000000515487E-05 V Inferior 0s 0101xs 0 0
V Superior 1s 0226xs 0 6.6000000515487E-05 V Inferior 1s 0201xs 0 0
V Superior 2s 0326xs 0 6.6000000515487E-05 V Inferior 2s 0301xs 0 0
…… ……
…… V Inferior 47s 4801xs 0 0
V Superior 48s 4926xs 0 6.6000000515487E-05 V Inferior 48s 4901xs 0 0
V Superior 49s 5026xs 0 6.6000000515487E-05 V Inferior 49s 5001xs 0 0
** **
* Izquierda * Derecha
R Izquierda 0s 0101ys 0 1E + 15 R Derecha 0s 5101ys 0 1E + 15
R Izquierda 1s 0102ys 0 1E + 15 R Derecha 1s 5102ys 0 1E + 15
…… ……
R Izquierda 23s 0124ys 0 1E + 15 R Derecha 23s 5124ys 0 1E + 15
R Izquierda 24s 0125ys 0 1E + 15 R Derecha 24s 5125ys 0 1E + 15
**
50
El bloque en el que se listan las variables, cuyos resultados de simulación se desean obtener en
forma tabulada en el archivo “.out”, constituye la siguiente sección del archivo de texto del modelo. Por
defecto, siempre se solicita la impresión de la tensión, tanto del circuito de concentraciones como del
circuito de funciones de corriente, en todos los centros de las celdas durante el transitorio, de acuerdo
con el intervalo de tiempo solicitado para la impresión. Cualquier otro valor que desee ser tabulado
debe solicitarse añadiendo al archivo las sentencias adecuadas. Para un modelo de 50x25 celdas el
listado que aparece por defecto es el siguiente:
El último bloque contiene un grupo de sentencias fijas. La que define la ventana del transitorio
de tiempos de la simulación, ya introducidos al diseñar el modelo, sentencia “.TRAN”, la asociada a la
precisión requerida en los cálculos, sentencia “.OPTIONS RELTOL”, la que define el número de dígitos
con que se presentan los resultados tabulados, sentencia “.OPTIONS NUMDIG”, la que activa el entorno
gráfico de PSpice, sentencia “.PROBE” y, finalmente, la sentencia de cierre del archivo modelo, “.END”.
Un ejemplo de estas sentencias es:
**
. TRAN 500 40000 0 UIC
.OPTIONS RELTOL 1E20
.OPTIONS NUMDGT 5
.PROBE
.END
**
El primero de los valores de la sentencia .TRAN define el intervalo de tiempo de impresión de los
datos tabulados de salida, el segundo valor se refiere al tiempo total solicitado en la simulación y el
tercero al valor a partir del cual se imprimen los resultados.
51
5.5 Pantallas de presentación de resultados
Como mencionamos en 5.2, una vez creado el archivo de modelo su ejecución es inmediata
pulsando el botón “Simulate” de la pantalla “Simulation”, figura 5.20, que aparece al pulsar el botón
“Generate CIR” o “Generate and process” de la pantalla principal. Esta acción arranca PSpice y ejecuta
el modelo presentando la pantalla que muestra la figura 5.22.
Figura 5.22 Pantalla de PSpice con información del proceso de simulación (Manual FAHET).
Es importante recordar que FAHET trabaja los decimales con punto y no con coma. En panel de
control, en configuración regional y de idioma debes cambiar esta opción si al simular con PSpice
aparece un mensaje de error en la pantalla de PSpice. El entorno de PSpice, figura 5.22, contiene tres
zonas separadas: la superior, donde se muestra el archivo de texto del modelo en ejecución; la inferior
izquierda que proporciona información de los posibles errores del modelo (caso de no convergencia o
errores en el circuito) y la inferior derecha donde se muestra el intervalo temporal de la simulación y el
tiempo de paso de la misma, así como el porcentaje de tiempo simulado. El programa PSpice varía
continuamente el tiempo de paso de simulación, de acuerdo con la tendencia uniforme o cambiante de
los resultados actuales, para reducir al máximo (sin merma de la precisión especificada por el
programador para los resultados numéricos) los tiempos de computación totales. Si se ha incluido la
opción de presentación de resultados en el entorno gráfico de PSpice, sentencia “.PROBE”, una vez
finalizada la simulación FAHET muestra directamente este entorno, figura 5.23, que consiste en una
cuadrícula vacía cuyo eje horizontal muestra una escala de tiempos cuya extensión es la del transitorio y
cuyo eje vertical, de momento sin escala, contendrá los valores numéricos de la variable solicitada
(concentraciones, funciones de corriente, flujos, etc.).
Esta misma pantalla puede abrirse durante el transitorio para observar los resultados parciales y
continuos de alguna variable mediante la opción “VIEW” + “SIMULATION RESULTS”, figura 5.24,
seguida de la opción “TRACE” + “ADD TRACE”, figura 5.25. Al pulsar este último botón se accede a un
listado de todas las variables, figura 5.26, en donde se seleccionan aquellas que queremos presentar
52
durante la simulación o al final de la misma. La opción “VIEW” + “SIMULATION RESULTS” es una
herramienta interesante en tanto que si tenemos una idea de los resultados o su tendencia podemos
permitir la continuación del cálculo o abortarlo para ahorrarnos la simulación de un modelo mal
diseñado. La escala vertical tiene siempre unidades de voltios V o amperios A, que hemos de traducir
directamente por concentración o función de corriente y por flujos de calor y de fluido según el
subcircuito del que se trate. El final de la simulación se muestra, en el propio entorno de PSpice, en la
pantalla de la figura 5.27.
Figura 5.24 Opción “VIEW”+”SIMULATION RESULTS” para visualizar progresivamente los resultados numéricos (M. FAHET).
53
Figura 5.25 Opción “TRACE”+”ADD TRACE” para seleccionar las variables que queremos representar (Manual FAHET).
Figura 5.26 Listado de todas las variables del problema (Manual FAHET).
54
Figura 5.27 Pantalla del entorno de PSpice indicando el final de la simulación (Manual FAHET).
Las posibilidades del entorno gráfico de PSpice son enormes y pueden consultarse en el manual
del programa. Así, por ejemplo, permiten representar simultáneamente o por separado los entornos
gráficos de distintos modelos y ajustar el tamaño de las pantallas para poder comparar sus resultados;
permite implementar un puntero que, al desplazarse muestra en un cuadro auxiliar el valor numérico
cambiante de la variable respecto a una referencia dada.
No obstante lo anterior, la forma en que PSpice presenta los resultados en su entorno gráfico
está limitada en parte, pues solo es capaz de presentar el transitorio de una variable o de un conjunto
de variables. No existe la posibilidad de mostrar representaciones instantáneas de un grupo de variables
correspondientes a un mismo modelo (por ejemplo perfiles de concentración o flujos de calor) ni otro
tipo de representaciones 2-D, opción que sí integra FAHET como veremos a continuación.
La figura 5.28 muestra una representación de variables típica del entorno de PSpice. El color
negro de fondo desaparece al transportar el gráfico al procesador de texto mediante las rutas de las
figuras 5.29 y 5.30 que permiten elegir nuevos colores para la pantalla, opción “cambio de colores”. La
figura 5.31 muestra el entorno PSpice una vez transportada al procesador Word.
55
Figura 5.28 Representación de variables en el entorno gráfico de PSpice (Manual FAHET).
Figura 5.29 Ruta para transportar la información gráfica de PSpice. Figura 5.30 Opciones de cambio de colores.
(Manual FAHET)
56
Figura 5.31 Pantalla de PSpice con fondo claro (Manual FAHET).
Los resultados tabulados pueden obtenerse, desde PSpice o bien desde Word, en el archivo de
salida de extensión “.out” que se crea tras la simulación. Este archivo, en el que también aparece al
principio del mismo todos los fallos de diseño del archivo de texto del modelo (cuando no se puede
simular éste) contiene ordenadamente un listado de todas las variables de las que se ha solicitado
información tabulada mediante la sentencia “.PRINT TRAN”. En columnas se muestra el valor de la
variable en cada instante, en los intervalos de tiempo especificados por la sentencia “.TRAN”. Estos
resultados pueden transportarse fácilmente (copiar y pegar) a una hoja de cálculo para su manipulación
y con ellos elaborar nuevas representaciones gráficas de perfiles, curvas de concentración constante,
etc. Un ejemplo de este listado para la variable concentración en los nodos (0101) y (0204) se muestra
en la figura 5.32 donde el intervalo de tiempo de paso es de 0.01 s.
Figura 5.32 Representación de los resultados tabulados en el archivo de texto de salida “.out” (Manual FAHET).
El entorno gráfico de FAHET, aunque limitado, permite acceder a una representación gráfica de
regiones (coloreadas) de concentración y función de corriente, obtenidas por interpolación de los
resultados tabulados anteriores, figura 5.33. El acceso a estas representaciones es directo una vez
simulado el modelo.
57
Figura 5.33 Regiones de isotemperatura del entorno gráfico de FAHET (Manual FAHET).
Por otro lado, FAHET incorpora la presentación de gráficos usando Matlab. Así, una vez
procesado el modelo por PSpice, aparece en la pantalla (en la parte inferior) el nuevo botón “Simulate”,
figura 5.20, que, al pulsarlo, nos muestra la pantalla auxiliar que arranca Matlab para estas nuevas
representaciones gráficas, mucho más precisas. Esta pantalla auxiliar, figura 5.34, tiene dos
subpantallas, “Display” y “Export”. En “Display” se accede en primer lugar a los botones “Display axis”
y “Current plane” que proporcionan diferentes vistas de la representación, figura 5.35, (estas opciones
han sido incluidas en FAHET con vistas a representaciones gráficas en 3-D en futuras versiones del
programa); se debe mantener siempre el eje Z para una mejor representación. En segundo lugar, el
botón “Subcircuit” permite seleccionar la variable a representar (temperatura o función de corriente),
mientras que el botón “time” se refiere al instante de tiempo solicitado dentro del transitorio, figura
5.36 y 5.37. Por último, el botón “General” presenta un desplegable con nuevos botones (“Shoe cell
borders”, “Show contours”,“Show circuit name”, “Border color”,“Shor colormap”, “Star color” y “End
final”) que permiten ajustar detalles de las representaciones, figuras 5.38 y 5.39.
Figura 5.34 Pantalla de procesado de malla para la manipulación con MATLAB (Manual FAHET).
58
Figura 5.35 Opciones de visualización (Manual FAHET).
Figura 5.37 Pantalla de selección del tiempo a representar dentro del transitorio (Manual FAHET).
59
Figura 5.39 Opciones de color de representación gráfica (Manual FAHET).
Ejecutada la simulación lo más conveniente es guardar el archivo de datos “.out” creado por
PSpice, que contiene toda la información de entrada del procesado de malla. El botón “Save out” de la
pantalla “Simulation results → Export” permite realizar esta operación, figura 5.40. Por defecto, este
archivo se guarda en la misma carpeta donde se encuentra el modelo. El archivo “Save DAT” sólo es útil
para representar gráficas en el entorno PSpice y ocupa una gran cantidad de memoria por lo que no
suele guardarse.
Figura 5.40 Pantalla de FAHET una vez simulado el modelo para el procesado de la malla (Manual FAHET).
60
El botón “Generate.m” presenta la ruta de acceso para guardar el gráfico Matlab que se va a
crear. Existen diversas opciones (botón “Opciones”) de representación, figura 5.41: “Fill contours”, para
presentar las regiones entre isolíneas con distinta coloración; “Grayscale image” para eliminar los
colores; “Show legend”, para incluir la leyenda de las líneas y tiempos representados, y “Define contour
lines”, para especificar las líneas requeridas (por defecto, se presentan un número típico de líneas entre
los valores máximo y mínimo definidos en las condiciones de contorno), figura 5.42. Para ver la figura
generada se pulsa el botón “Generar gráfico Matlab”. Esta acción arranca Matlab y procede a la
generación del gráfico, que se presenta automáticamente al cabo de unos pocos segundos en el formato
seleccionado por “Opciones”.
Figura 5.41 Opciones de representación gráfica de isolineas con MATLAB (Manual FAHET).
Figura 5.42 Selección de los valores de las isolineas MATLAB (Manual FAHET).
La figura 5.43 a-d muestra cuatro gráficos típicos (en blanco y negro y color), con la leyenda de
isolíneas de temperatura y funcion de corriente.
61
Figura 5.43 a-b-c-d Gráficos típicos de MATLAB (Manual FAHET).
FAHET contiene una subrutina que genera un archivo Excel con los datos de salida del archivo
“.out” de PSpice (concentración o función de corriente según se haya seleccionado previamente), para
su tratamiento y/o representación. Esta acción se realiza mediante el botón “Exportar datos Excel” de la
pantalla “Simulation result → Export Data Excel”, figura 5.44 que muestra la ruta para guardar este
archivo.
62
La opción “Generate Animation” de la pantalla “Export” da acesso a la pantalla de la figura 5.45.
FAHET puede generar animaciones en distintos formatos, figura 5.46. Se especifica el intervalo de
tiempo de la animación, el retardo entre imágenes y el tamaño. El número de animaciones es el
indicado en la sentencia “.TRAN” del archivo de modelo “.cir” de PSpice. Una vez seleccionadas las
opciones de animación, se especifica la carpeta donde se guardará la animación, figura 5.47. La opción
“Separate Matlab Figures” permite generar la animación y representar gráficamente las imágenes
generadas en la ruta especificada.
63
Figura 5.46 Formato de las animaciones (Manual FAHET)
Figura 5.47 Ruta de acceso para guardar las figuras de la representación animada (Manual FAHET).
Las figuras y animaciones Matlab ya guardadas, al abrirlas dan acceso a la pantalla Matlab de la
figura 5.41. Este archivo puede manipularse si el usuario conoce el código para generar la figura
guardada u otras cuyo número de isolíneas, color, etc., pueden especificarse en esta pantalla definiendo
las líneas de contorno.
64
CAPITULO 6. CAMPO DE TEMPERATURAS Y FLUJOS EN POZOS EN LAS
CERCANÍAS DE FOCOS TÉRMICOS. IMPLEMENTACIÓN EN FAHET
En la primera fila se indican las tres configuraciones consideradas para el pozo de extracción:
Caso I: Pozo de extracción con función de flujo de valor 4.9x10 -6 a la izquierda del pozo,
3.9x10-6 en el fondo del pozo y 3.0x10-6 a la derecha del pozo.
Caso II: Pozo de extracción con función de flujo de valor 4.9x10 -6 a la izquierda del pozo,
3.0x10-6 en el fondo del pozo y 1.0x10-6 a la derecha del pozo.
Caso III : Pozo de extracción con función de flujo de valor 4.9x10 -6 a la izquierda del pozo,
5.9x10-6 en el fondo del pozo y 6.7x10-6 a la derecha del pozo.
En las columnas las variables calculadas para las influencias indicadas, con un intervalo de
tiempo de 5.0x1011s. y un incremento de tiempo de 5.0x109 s.
65
CONFIGURACIÓN GEOMÉTRICA DEL POZO DE EXTRACCIÓN
CASO I CASO II CASO III
-6 -6
Ψizq=4.9X10 Ψizq=4.9X10 Ψizq=4.9X10-6
Tabla. 6.1 Esquema de las distintas configuraciones de análisis con las variables calculadas en cada caso.
66
6.2 Características de los dominios de simulación, condiciones de frontera y
datos de entrada
CASO I
El dominio, de altura 1000 m y anchura 6250 m, será modelado a través una malla regular,
constituida por 20 x 50 celdas (20 celdas sobre la vertical z y 50 sobre la horizontal x), de dimensiones
∆𝑥 = 125 𝑚. ∆𝑦 = 50 𝑚. Las principales características físico-químicas a implementar son la
permeabilidad del subsuelo (constante en las dos direcciones x y z), su viscosidad dinámica y su
difusividad térmica. Se considera una velocidad del fluido igual a 9.8 x 10- 8 m/s en sentido izquierda-
derecha, dos paredes laterales adiabáticas y aisladas (con valores de flujo y temperatura adiabáticas), un
foco caliente en la parte baja izquierda que actúa sobre 5, 10, 15 y 20 celdas (con valor de temperatura
adimensional constante igual a 1), una función de flujo con valor igual acero en las cincuenta celdas
inferiores, y un valor de temperatura adiabático para las restantes celdas de la parte inferior del
dominio. La parte superior consta de un valor adimensional de temperatura igual a cero y una función
de flujo de valor 4.9x10 -6 a la izquierda del pozo, 3.9x10-6 en la parte inferior del pozo y 3.0x10-6a la
derecha del pozo. Se presentan en la figura 6.2-I.1 las principales características del dominio.
𝑇=0 𝑇=0
Ψ = 4.9x10−6 Ψ = 3.0x10−6
1250 m
1000 m
𝑇 = Adiabática
𝑇 = Adiabática
Ψ = Adiabática
Ψ = Adiabática 600m
625 m
1250 m
1875 m 𝑇=1 Ψ= 0 𝑇 = Adiabática Ψ=0
2500 m
6250 m
Figura 6.2-I.1 Esquema del modelo físico con las condiciones de contorno.
67
En la tabla 6.2-I.1 se muestran las características geométricas y físico-químicas del dominio de
análisis.
Cuando introducimos un pozo, en este caso con una profundidad de 600 m., hay que tener en
cuenta las condiciones de contorno adicionales que caracterizan las discontinuidades producidas (en
términos de transmisión de flujo térmico y de flujo de fluido) en el dominio de análisis: en la parte
superior del pozo (coincidente con la del dominio), en los bordes derecho e izquierdo del pozo y en la
celda inferior del pozo. Para solucionar este problema, se modelan las celdas en las direcciones 𝑥 𝑦 𝑧,
asignando condiciones de frontera en función de las componentes del vector de velocidad del flujo de
agua 𝑣 , a partir de las cuales se puede obtener el valor del flujo correspondiente Ψ en cada cara de la
celda.
Tomando una velocidad de fluido igual a 9.8 x 10-8 m/s, podemos definir para cada caso la
condición de frontera según se indica en la siguiente tabla 6.2-I.2 Por tanto, en la parte superior
izquierda del dominio y en las celdas del borde izquierdo del pozo la condición de contorno de la función
de flujo (Ψ𝐼𝑍 ) será igual a 4.9 x 10-6 m2/s, en la celda inferior del pozo (Ψ𝐼𝑁 ) será igual a 3.9 x 10-6 m2/s,
en las celdas del borde derecho y en la parte superior derecha del modelo (Ψ𝐷𝐸 ) será igual a 3.0 x 10-6
m2/s. Las restantes condiciones de frontera, en términos de temperatura, serán de paredes adiabáticas,
según se indica en la figura 6.2-I.1
El caudal total de agua en una celda del dominio se puede calcular a partir de la velocidad y las
dimensiones de las celdas (dx y dz),
𝑄𝑡𝑜𝑡 = 𝑣 × 𝑑𝑥 × 𝑑𝑧
68
CONFIGURACIÓN 2: INFLUENCIA DE LA PROFUNDIDAD DEL POZO
Con la segunda configuración se han estudiado las cuatro situaciones reflejadas en la figura
6.2-I.2, donde se indican las distintas profundidades del pozo estudiadas. Se ha elegido para la posición
del pozo a 1250 m del borde izquierdo y el tamaño del foco térmico es de 15 celdas que equivalen a
1875 m.
𝑇=0 𝑇=0
Ψ = 4.9x10−6 Ψ = 3.0x10−6
1250 m
200m
1000 m
400m 𝑇 = Adiabática
𝑇 = Adiabática
Ψ = Adiabática
Ψ = Adiabática
600m
800m
6250 m
Figura 6.2-I.2 Esquema del modelo físico con las condiciones de contorno.
Como antes, el pozo implica imponer diferentes condiciones de frontera para las celdas de
referencia elegidas en términos de temperatura (adiabática si la celda de referencia colinda con el pozo)
y función de flujo.
Las condiciones de contorno para el dominio de análisis se calculan igual que para la
configuración 1 del caso I.
69
CONFIGURACIÓN 3: INFLUENCIA DEL FACTOR DE PERMEABILIDAD K
Con la tercera configuración se han estudiado las cuatro situaciones reflejadas en la figura
6.2-I.3, donde se indican las distintos valores de permeabilidad K estudiados. Se ha elegido para la
posición del pozo a 1250 m del borde izquierdo, con una profundidad de 600 m. y el tamaño del foco
térmico es de 15 celdas que equivalen a 1875 m.
𝑇=0 𝑇=0
Ψ = 4.9x10−6 Ψ = 3.0x10−6
1250 m
1000 m
𝑇 = Adiabática
𝑇 = Adiabática
Ψ = Adiabática
Ψ = Adiabática 600m
6250 m
Figura 6.2-I.3 Esquema del modelo físico con las condiciones de contorno.
K1 = 1x10-13 m2
K2 = 5x10-14 m2
K3 = 1x10-14 m2
K4 = 5x10-15 m2
Las condiciones de contorno para el dominio de análisis se calculan igual que para la
configuración 1 del caso I.
70
CONFIGURACIÓN 4: INFLUENCIA DEL FACTOR DE DIFUSIVIDAD D
Con la cuarta configuración se han estudiado las cuatro situaciones reflejadas en la figura
6.2-I.4, donde se indican los distintos valores de difusividad D estudiados. Se ha elegido para la posición
del pozo a 1250 m del borde izquierdo, con una profundidad de 600 m. y el tamaño del foco térmico es
de 15 celdas que equivalen a 1875 m.
𝑇=0 𝑇=0
Ψ = 4.9x10−6 Ψ = 3.0x10−6
1250 m
1000 m
𝑇 = Adiabática
𝑇 = Adiabática
Ψ = Adiabática
Ψ = Adiabática 600m
6250 m
Figura 6.2-I.4 Esquema del modelo físico con las condiciones de contorno.
D1 = 1x10-5 m2/s
D2 = 5x10-6 m2/s
D3 = 1x10-6 m2/s
D4 = 5x10-7 m2/s
Las condiciones de contorno para el dominio de análisis se calculan igual que para la
configuración 1 del caso I.
71
CASO II
El dominio, de altura 1000 m y anchura 6250 m, será modelado a través una malla regular,
constituida por 20 x 50 celdas (20 celdas sobre la vertical z y 50 sobre la horizontal x), de dimensiones
∆𝑥 = 125 𝑚. ∆𝑦 = 50 𝑚. Las principales características físico-químicas a implementar son la
permeabilidad del subsuelo (constante en las dos direcciones x y z), su viscosidad dinámica y su
difusividad térmica. Se considera una velocidad del fluido igual a 9.8 x 10- 8 m/s en sentido izquierda-
derecha, dos paredes laterales adiabáticas y aisladas (con valores de flujo y temperatura adiabáticas), un
foco caliente en la parte baja izquierda que actúa sobre 5, 10, 15 y 20 celdas (con valor de temperatura
adimensional constante igual a 1), una función de flujo con valor igual a cero en las cincuenta celdas
inferiores, y un valor de temperatura adiabático para las restantes celdas de la parte inferior del
dominio. La parte superior consta de un valor adimensional de temperatura igual a cero y una función
de flujo de valor 4.9x10 -6 a la izquierda del pozo, 3.0x10-6 en la parte inferior del pozo y 1.0x10-6 a la
derecha del pozo. Se presentan en la figura 6.2-II.1 las principales características del dominio.
𝑇=0 𝑇=0
Ψ = 4.9x10−6 Ψ = 1.0x10−6
1250 m
1000 m
𝑇 = Adiabática
𝑇 = Adiabática
Ψ = Adiabática
Ψ = Adiabática 600m
625 m
1250 m
1875 m 𝑇=1 Ψ= 0 𝑇 = Adiabática Ψ=0
2500 m
6250 m
Figura 6.2-II.1 Esquema del modelo físico con las condiciones de contorno.
72
En la tabla 6.2-II.1 se muestran las características geométricas y físico-químicas del dominio de
análisis.
Tomando una velocidad de fluido igual a 9.8 x 10-8 m/s, podemos definir para cada caso la condición de
frontera según se indica en la siguiente tabla 6.2-II.2 Por tanto, en la parte superior izquierda del
dominio y en las celdas del borde izquierdo del pozo la condición de contorno de la función de flujo
(Ψ𝐼𝑍 ) será igual a 4.9 x 10-6 m2/s, en la celda inferior del pozo (Ψ𝐼𝑁 ) será igual a 3.0 x 10-6 m2/s, en las
celdas del borde derecho y en la parte superior derecha del modelo (Ψ𝐷𝐸 ) será igual a 1.0 x 10-6 m2/s.
Las restantes condiciones de frontera, en términos de temperatura, serán de paredes adiabáticas, según
se indica en la figura 6.2-II.1.
El caudal total de agua en una celda del dominio se puede calcular a partir de la velocidad y las
dimensiones de las celdas (dx y dz),
𝑄𝑡𝑜𝑡 = 𝑣 × 𝑑𝑥 × 𝑑𝑧
73
CONFIGURACIÓN 2: INFLUENCIA DE LA PROFUNDIDAD DEL POZO
Con la segunda configuración se han estudiado las cuatro situaciones reflejadas en la figura
6.2-II.2, donde se indican las distintas profundidades del pozo estudiadas. Se ha elegido para la posición
del pozo a 1250 m del borde izquierdo y el tamaño del foco térmico es de 15 celdas que equivalen a
1875 m.
𝑇=0 𝑇=0
Ψ = 4.9x10−6 Ψ = 1.0x10−6
1250 m
200m
1000 m
400m 𝑇 = Adiabática
𝑇 = Adiabática
Ψ = Adiabática
Ψ = Adiabática
600m
800m
6250 m
Figura 6.2-II.2 Esquema del modelo físico con las condiciones de contorno.
Como antes, el pozo implica imponer diferentes condiciones de frontera para las celdas de
referencia elegidas en términos de temperatura (adiabática si la celda de referencia colinda con el pozo)
y función de flujo.
Las condiciones de contorno para el dominio de análisis se calculan igual que para la
configuración 1 del caso II.
74
CONFIGURACIÓN 3: INFLUENCIA DEL FACTOR DE PERMEABILIDAD K
Con la tercera configuración se han estudiado las cuatro situaciones reflejadas en la figura
6.2-II.3, donde se indican las distintos valores de permeabilidad K estudiados. Se ha elegido para la
posición del pozo a 1250 m del borde izquierdo, con una profundidad de 600 m. y el tamaño del foco
térmico es de 15 celdas que equivalen a 1875 m.
𝑇=0 𝑇=0
Ψ = 4.9x10−6 Ψ = 1.0x10−6
1250 m
1000 m
𝑇 = Adiabática
𝑇 = Adiabática
Ψ = Adiabática
Ψ = Adiabática 600m
6250 m
Figura 6.2-II.3 Esquema del modelo físico con las condiciones de contorno.
K1 = 1x10-13 m2
K2 = 5x10-14 m2
K3 = 1x10-14 m2
K4 = 5x10-15 m2
Las condiciones de contorno para el dominio de análisis se calculan igual que para la
configuración 1 del caso II.
75
CONFIGURACIÓN 4: INFLUENCIA DEL FACTOR DE DIFUSIVIDAD D
Con la cuarta configuración se han estudiado las cuatro situaciones reflejadas en la figura
6.2-II.4, donde se indican los distintos valores de difusividad D estudiados. Se ha elegido para la posición
del pozo a 1250 m del borde izquierdo, con una profundidad de 600 m. y el tamaño del foco térmico es
de 15 celdas que equivalen a 1875 m.
𝑇=0 𝑇=0
Ψ = 4.9x10−6 Ψ = 1.0x10−6
1250 m
1000 m
𝑇 = Adiabática
𝑇 = Adiabática
Ψ = Adiabática
Ψ = Adiabática 600m
6250 m
Figura 6.2-II.4 Esquema del modelo físico con las condiciones de contorno.
D1 = 1x10-5 m2/s
D2 = 5x10-6 m2/s
D3 = 1x10-6 m2/s
D4 = 5x10-7 m2/s
Las condiciones de contorno para el dominio de análisis se calculan igual que para la
configuración 1 del caso II.
76
CASO III
El dominio, de altura 1000 m y anchura 6250 m, será modelado a través una malla regular,
constituida por 20 x 50 celdas (20 celdas sobre la vertical z y 50 sobre la horizontal x), de dimensiones
∆𝑥 = 125 𝑚. ∆𝑦 = 50 𝑚. Las principales características físico-químicas a implementar son la
permeabilidad del subsuelo (constante en las dos direcciones x y z), su viscosidad dinámica y su
difusividad térmica. Se considera una velocidad del fluido igual a 9.8 x 10- 8 m/s en sentido izquierda-
derecha, dos paredes laterales adiabáticas y aisladas (con valores de flujo y temperatura adiabáticas), un
foco caliente en la parte baja izquierda que actúa sobre 5, 10, 15 y 20 celdas (con valor de temperatura
adimensional constante igual a 1), una función de flujo con valor igual a cero en las cincuenta celdas
inferiores, y un valor de temperatura adiabático para las restantes celdas de la parte inferior del
dominio. La parte superior consta de un valor adimensional de temperatura igual a cero y una función
de flujo de valor 4.9x10 -6 a la izquierda del pozo, 5.9x10-6 en la parte inferior del pozo y 6.7x10-6 a la
derecha del pozo. Se presentan en la figura 6.2-III.1 las principales características del dominio.
𝑇=0 𝑇=0
Ψ = 4.9x10−6 Ψ = 6.7x10−6
1250 m
1000 m
𝑇 = Adiabática
𝑇 = Adiabática
Ψ = Adiabática
Ψ = Adiabática 600m
625 m
1250 m
1875 m 𝑇=1 Ψ= 0 𝑇 = Adiabática Ψ=0
2500 m
6250 m
Figura 6.2-III.1 Esquema del modelo físico con las condiciones de contorno.
77
En la tabla 6.2-III.1 se muestran las características geométricas y físico-químicas del dominio de
análisis.
Tomando una velocidad de fluido igual a 9.8 x 10-8 m/s, podemos definir para cada caso la
condición de frontera según se indica en la siguiente tabla 6.2-III.2 Por tanto, en la parte superior
izquierda del dominio y en las celdas del borde izquierdo del pozo la condición de contorno de la función
de flujo (Ψ𝐼𝑍 ) será igual a 4.9 x 10-6 m2/s, en la celda inferior del pozo (Ψ𝐼𝑁 ) será igual a 5.9 x 10-6 m2/s,
en las celdas del borde derecho y en la parte superior derecha del modelo (Ψ𝐷𝐸 ) será igual a 6.7 x 10-6
m2/s. Las restantes condiciones de frontera, en términos de temperatura, serán de paredes adiabáticas,
según se indica en la figura 6.2-III.1
El caudal total de agua en una celda del dominio se puede calcular a partir de la velocidad y las
dimensiones de las celdas (dx y dz),
𝑄𝑡𝑜𝑡 = 𝑣 × 𝑑𝑥 × 𝑑𝑧
78
CONFIGURACIÓN 2: INFLUENCIA DE LA PROFUNDIDAD DEL POZO
Con la segunda configuración se han estudiado las cuatro situaciones reflejadas en la figura
6.2-III.2, donde se indican las distintas profundidades del pozo estudiadas. Se ha elegido para la posición
del pozo a 1250 m del borde izquierdo y el tamaño del foco térmico es de 15 celdas que equivalen a
1875 m.
𝑇=0 𝑇=0
Ψ = 4.9x10−6 Ψ = 6.7x10−6
1250 m
200m
1000 m
400m 𝑇 = Adiabática
𝑇 = Adiabática
Ψ = Adiabática
Ψ = Adiabática
600m
800m
6250 m
Figura 6.2-III.2 Esquema del modelo físico con las condiciones de contorno.
Como antes, el pozo implica imponer diferentes condiciones de frontera para las celdas de
referencia elegidas en términos de temperatura (adiabática si la celda de referencia colinda con el pozo)
y función de flujo.
Las condiciones de contorno para el dominio de análisis se calculan igual que para la
configuración 1 del caso III.
79
CONFIGURACIÓN 3: INFLUENCIA DEL FACTOR DE PERMEABILIDAD K
Con la tercera configuración se han estudiado las cuatro situaciones reflejadas en la figura
6.2-III.3, donde se indican las distintos valores de permeabilidad K estudiados. Se ha elegido para la
posición del pozo a 1250 m del borde izquierdo, con una profundidad de 600 m. y el tamaño del foco
térmico es de 15 celdas que equivalen a 1875 m.
𝑇=0 𝑇=0
Ψ = 4.9x10−6 Ψ = 6.7x10−6
1250 m
1000 m
𝑇 = Adiabática
𝑇 = Adiabática
Ψ = Adiabática
Ψ = Adiabática 600m
6250 m
Figura 6.2-III.3 Esquema del modelo físico con las condiciones de contorno.
K1 = 1x10-13 m2
K2 = 5x10-14 m2
K3 = 1x10-14 m2
K4 = 5x10-15 m2
Las condiciones de contorno para el dominio de análisis se calculan igual que para la
configuración 1 del caso III.
80
CONFIGURACIÓN 4: INFLUENCIA DEL FACTOR DE DIFUSIVIDAD D
Con la cuarta configuración se han estudiado las cuatro situaciones reflejadas en la figura
6.2-III.4, donde se indican los distintos valores de difusividad D estudiados. Se ha elegido para la
posición del pozo a 1250 m del borde izquierdo, con una profundidad de 600 m. y el tamaño del foco
térmico es de 15 celdas que equivalen a 1875 m.
𝑇=0 𝑇=0
Ψ = 4.9x10−6 Ψ = 6.7x10−6
1250 m
1000 m
𝑇 = Adiabática
𝑇 = Adiabática
Ψ = Adiabática
Ψ = Adiabática 600m
6250 m
Figura 6.2-III.4 Esquema del modelo físico con las condiciones de contorno.
D1 = 1x10-5 m2/s
D2 = 5x10-6 m2/s
D3 = 1x10-6 m2/s
D4 = 5x10-7 m2/s
Las condiciones de contorno para el dominio de análisis se calculan igual que para la
configuración 1 del caso III.
Se han calculado T(𝑥, 𝑧) y Ψ 𝑥, 𝑧 en todo el dominio y su derivada con respecto a 𝑧 para los
distintos intervalos de tiempo indicados.
81
6.3 Resultados
A continuación se presentan los resultados obtenidos para las tres CASOS y cuatro
CONFIGURACIONES del dominio de estudio y se resumen los principales resultados de análisis obtenidos
para las diferentes opciones.
CASO I
En la figura 6.3-I.1-𝑇.1 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 1.0 x 1010 s. y los distintos tamaños del foco térmico considerados (5, 10, 15
y 20 celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
Podemos concluir que para el intervalo de tiempo considerado, el tamaño del foco térmico sólo
influye en la distancia que las isolíneas ocupan en la dirección horizontal (eje 𝑥).
10
Figura 6.3-I.1-T.1 Isolíneas de la temperatura adimensional para un intervalo de tiempo 1.0x10 s en el dominio de estudio.
82
En la figura 6.3-I.1-𝑇.2 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 3.0 x 1010 s. y los distintos tamaños del foco térmico considerados (5, 10, 15
y 20 celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
Se puede apreciar la evolución de las isolíneas, con cierta curvatura, hacia la superficie del
dominio y el tamaño del foco térmico sólo influye en la distancia que las isolíneas ocupan en la dirección
horizontal (eje 𝑥).
10
Figura 6.3-I.1-T.2 Isolíneas de la temperatura adimensional para un intervalo de tiempo 3.0x10 s en el dominio de estudio.
83
En la figura 6.3-I.1-𝑇.3 se presentan las gráficas con las isolíneas de temperatura adimensional
para unintervalo de tiempo 9.0 x 1010 s. y los distintos tamaños del foco térmico considerados (5, 10, 15
y 20 celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
Las isolíneas siguen evolucionando hacia la superficie del dominio y tienen una mayor curvatura
con valores máximos. Cuando el foco térmico (15 y 20 celdas) alcanza la vertical del pozo, las isolíneas
de la temperatura a cada lado del pozo tienen distintos valores.
10
Figura 6.3-I.1-T.3 Isolíneas de la temperatura adimensional para un intervalo de tiempo 9.0x10 s en el dominio de estudio.
84
En la figura 6.3-I.1-𝑇.4 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 1.2 x 1011 s. y los distintos tamaños del foco térmico considerados (5, 10, 15
y 20 celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
11
Figura 6.3-I.1-T.4 Isolíneas de la temperatura adimensional para un intervalo de tiempo 1.2x10 s en el dominio de estudio.
85
En la figura 6.3-I.1-𝑇.5 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 3.2 x 1011 s. y los distintos tamaños del foco térmico considerados (5, 10, 15
y 20 celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
11
Figura 6.3-I.1-T.5 Isolíneas de la temperatura adimensional para un intervalo de tiempo 3.2x10 s en el dominio de estudio.
86
En la figura 6.3-I.1-𝑇.6 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 5.0 x 1011 s. y los distintos tamaños del foco térmico considerados (5, 10, 15
y 20 celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
11
Figura 6.3-I.1-T.6 Isolíneas de la temperatura adimensional para un intervalo de tiempo 5.0x10 s en el dominio de estudio.
87
CASO I.1-𝑻: Influencia del tamaño del foco térmico
88
En segundo lugar, vamos a estudiar la distribución de flujo de corriente Ψ en el dominio
geotérmico.
En la figura 6.3-I.1-Ψ.1 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 1.0 x 1010 s. y los distintos tamaños del foco térmico considerados (5, 10, 15 y 20
celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
Se puede apreciar la influencia del foco al desplazarse hacia la derecha los vértices inferiores de
las isolíneas con la ampliación del tamaño del foco térmico hacia la derecha.
10
Figura 6.3-I.1-𝚿.1 Isolíneas de la funciónde corriente para un intervalo de tiempo 1.0x10 s en el dominio de estudio.
89
En la figura 6.3-I.1-Ψ.2 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 3.0 x 1010 s. y los distintos tamaños del foco térmico considerados (5, 10, 15 y 20
celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
10
Figura 6.3-I.1-𝚿.2 Isolíneas de la función de corriente para un intervalo de tiempo 3.0x10 s en el dominio de estudio.
90
En la figura 6.3-I.1-Ψ.3 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 9.0 x 1010 s. y los distintos tamaños del foco térmico considerados (5, 10, 15 y 20
celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
Se desarrolla la formación de isolíneas ovoidales en la zona situada sobre la vertical del foco
térmico, que se va desplazando hacia la derecha del dominio.
10
Figura 6.3-I.1-𝚿.3 Isolíneas de la función de corriente para un intervalo de tiempo 9.0x10 s en el dominio de estudio.
91
En la figura 6.3-I.1-Ψ.4 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 1.2 x 1011 s. y los distintos tamaños del foco térmico considerados (5, 10, 15 y 20
celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
11
Figura 6.3-I.1-𝚿.4 Isolíneas de la función de corriente para un intervalo de tiempo 1.2x10 s en el dominio de estudio.
92
En la figura 6.3-I.1-Ψ.5 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 3.2 x 1011 s. y los distintos tamaños del foco térmico considerados (5, 10, 15 y 20
celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
Se desarrolla la formación de isolíneas ovoidales como en el caso anterior. En todos los casos
aumenta el valor de las isolíneas más interiores.
11
Figura 6.3-I.1-𝚿.5 Isolíneas de la función de corriente para un intervalo de tiempo 3.2x10 s en el dominio de estudio.
93
En la figura 6.3-I.1-Ψ.6 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 5.0 x 1011 s. y los distintos tamaños del foco térmico considerados (5, 10, 15 y 20
celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
Las formas y valores de las isolíneas parecen haber llegado a su máximo desarrollo, pues no hay
grandes cambios con respecto al intervalo del caso anterior.
11
Figura 6.3-I.1-𝚿.6 Isolíneas de la función de corriente para un intervalo de tiempo 5.0x10 s en el dominio de estudio.
94
CASO I.1-𝚿: Influencia del tamaño del foco térmico
95
6.3_I.2 Configuración de análisis 2: Influencia de la profundidad del pozo.
En la figura 6.3-I.2-𝑇.1 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 1.0 x 1010 s. y las distintas profundidades del pozo consideradas (200 m,
400 m, 600 m y 800 m) en todo el dominio de estudio.
Podemos concluir que para el intervalo de tiempo considerado, la profundidad del pozo no
influye sobre las isolíneas, sólo en el pozo de 800 m las isolíneas pierden su forma horizontal.
10
Figura 6.3-I.2-T.1 Isolíneas de la temperatura adimensional para un intervalo de tiempo 1.0x10 s en el dominio de estudio.
96
En la figura 6.3-I.2-𝑇.2 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 3.0 x 1010 s. y las distintas profundidades del pozo consideradas (200 m,
400 m, 600 m y 800 m) en todo el dominio de estudio.
Las isolíneas evolucionan hacia la superficie del dominio, por tanto la presencia del pozo
modifica la forma horizontal de las isolíneas.
10
Figura 6.3-I.2-T.2 Isolíneas de la temperatura adimensional para un intervalo de tiempo 3.0x10 s en el dominio de estudio.
97
En la figura 6.3-I.2-𝑇.3 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 9.0 x 1010 s. y las distintas profundidades del pozo consideradas (200 m,
400 m, 600 m y 800 m) en todo el dominio de estudio.
Las isolíneas siguen evolucionando hacia la superficie del dominio y se inicia la forma ovoidal.
10
Figura 6.3-I.2-T.3 Isolíneas de la temperatura adimensional para un intervalo de tiempo 9.0x10 s en el dominio de estudio.
98
En la figura 6.3-I.2-𝑇.4 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 1.2 x 1011 s. y las distintas profundidades del pozo consideradas (200 m,
400 m, 600 m y 800 m) en todo el dominio de estudio.
11
Figura 6.3-I.2-T.4 Isolíneas de la temperatura adimensional para un intervalo de tiempo 1.2x10 s en el dominio de estudio.
99
En la figura 6.3-I.2-𝑇.5 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 3.2 x 1011 s. y las distintas profundidades del pozo consideradas (200 m,
400 m, 600 m y 800 m) en todo el dominio de estudio.
11
Figura 6.3-I.2-T.5 Isolíneas de la temperatura adimensional para un intervalo de tiempo 3.2x10 s en el dominio de estudio.
100
En la figura 6.3-I.2-𝑇.6 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 5.0 x 1011 s. y las distintas profundidades del pozo consideradas (200 m,
400 m, 600 m y 800 m) en todo el dominio de estudio.
11
Figura 6.3-I.2-T.6 Isolíneas de la temperatura adimensional para un intervalo de tiempo 5.0x10 s en el dominio de estudio.
101
CASO I.2-𝑻: Influencia de la profundidad del pozo
102
En segundo lugar, vamos a estudiar la distribución de flujo de corriente Ψen el dominio
geotérmico.
En la figura 6.3-I.2-Ψ.1 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 1.0 x 1010 s. y las distintas profundidades del pozo consideradas (200 m, 400 m,
600 m y 800 m) en todo el dominio de estudio.
Se puede apreciar la influencia del pozo al encontrarse isolíneas de valores mayores hacia la
parte superior del dominio y a la izquierda del pozo. Encontramos dos distorsiones de la forma de las
curvas, una sobre la vertical del foco térmico y otra en la posición donde se encuentra el pozo en cada
caso. Se puede apreciar que cuando más profundo es el pozo, las isolíneas de corriente son menos
horizontales.
10
Figura 6.3-I.2-𝚿.1 Isolíneas de la función de corriente para un intervalo de tiempo 1.0x10 s en el dominio de estudio.
103
En la figura 6.3-I.2-Ψ.2 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 3.0 x 1010 s. y las distintas profundidades del pozo consideradas (200 m, 400 m,
600 m y 800 m) en todo el dominio de estudio.
10
Figura 6.3-I.2-𝚿.2 Isolíneas de la función de corriente para un intervalo de tiempo 3.0x10 s en el dominio de estudio.
104
En la figura 6.3-I.2-Ψ.3 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 9.0 x 1010 s. y las distintas profundidades del pozo consideradas (200 m, 400 m,
600 m y 800 m) en todo el dominio de estudio.
La forma ovoidal de las isolíneas se desarrolla en todos los casos. La diferencia está en la zona
alrededor de la parte inferior del pozo.
10
Figura 6.3-I.2-𝚿.3 Isolíneas de la función de corriente para un intervalo de tiempo 9.0x10 s en el dominio de estudio.
105
En la figura 6.3-I.2-Ψ.4 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 1.2 x 1011 s. y las distintas profundidades del pozo consideradas (200 m, 400 m,
600 m y 800 m) en todo el dominio de estudio.
11
Figura 6.3-I.2-𝚿.4 Isolíneas de la función de corriente para un intervalo de tiempo 1.2x10 s en el dominio de estudio.
106
En la figura 6.3-I.2-Ψ.5 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 3.2 x 1011 s. y las distintas profundidades del pozo consideradas (200 m, 400 m,
600 m y 800 m) en todo el dominio de estudio.
11
Figura 6.3-I.2-𝚿.5 Isolíneas de la función de corriente para un intervalo de tiempo 3.2x10 s en el dominio de estudio.
107
En la figura 6.3-I.2-Ψ.6 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 5.0 x 1011 s. y las distintas profundidades del pozo consideradas (200 m, 400 m,
600 m y 800 m) en todo el dominio de estudio.
Las formas y valores de las isolíneas parecen haber llegado a su máximo desarrollo, pues no hay
grandes cambios con respecto al intervalo de tiempo anterior. Podemos concluir que la profundidad del
pozo apenas afecta sobre la respuesta del dominio geotérmico y lo hace en la zona alrededor de la parte
inferior del pozo.
11
Figura 6.3-I.2-𝚿.6 Isolíneas de la función de corriente para un intervalo de tiempo 5.0x10 s en el dominio de estudio.
108
CASO I.2-𝚿: Influencia de la profundidad del pozo
109
6.3_I.3 Configuración de análisis 3: Influencia del factor de permeabilidad K
En la figura 6.3-I.3-𝑇.1 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 1.0 x 1010 s. y los distintos valores del factor de permeabilidad K
(K1 = 1x10-13 m2, K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
Se aprecia la evolución de las isolíneas hacia la parte superior del dominio y la forma ovoidal
para el mayor valor del factor de permeabilidad K.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
10
Figura 6.3-I.3-𝑻.1 Isolíneas de la temperatura adimensional para un intervalo de tiempo 1.0x10 s en el dominio de estudio.
110
En la figura 6.3-I.3-𝑇.2 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 3.0 x 1010 s. y los distintos valores del factor de permeabilidad K
(K1 = 1x10-13 m2, K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
Aumenta la evolución de las isolíneas hacia la parte superior del dominio y la forma ovoidal para
los dos mayores valores del factor de permeabilidad K.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
10
Figura 6.3-I.3-𝑻.2 Isolíneas de la temperatura adimensional para un intervalo de tiempo 3.0x10 s en el dominio de estudio.
111
En la figura 6.3-I.3-𝑇.3 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 9.0 x 1010 s. y los distintos valores del factor de permeabilidad K
(K1 = 1x10-13 m2, K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
10
Figura 6.3-I.3-𝑻.3 Isolíneas de la temperatura adimensional para un intervalo de tiempo 9.0x10 s en el dominio de estudio.
112
En la figura 6.3-I.3-𝑇.4 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 1.2 x 1011 s. y los distintos valores del factor de permeabilidad K
(K1 = 1x10-13 m2, K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
11
Figura 6.3-I.3-𝑻.4 Isolíneas de la temperatura adimensional para un intervalo de tiempo 1.2x10 s en el dominio de estudio.
113
En la figura 6.3-I.3-𝑇.5 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 3.2 x 1011 s. y los distintos valores del factor de permeabilidad K
(K1 = 1x10-13 m2, K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
11
Figura 6.3-I.3-𝑻.5 Isolíneas de la temperatura adimensional para un intervalo de tiempo 3.2x10 s en el dominio de estudio.
114
En la figura 6.3-I.3-𝑇.6 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 5.0 x 1011 s. y los distintos valores del factor de permeabilidad K
(K1 = 1x10-13 m2, K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
11
Figura 6.3-I.3-𝑻.6 Isolíneas de la temperatura adimensional para un intervalo de tiempo 5.0x10 s en el dominio de estudio.
115
CASO I.3-𝑻: Influencia del factor de permeabilidad K
116
En segundo lugar, vamos a estudiar la distribución de flujo de corriente Ψ en el dominio
geotérmico.
En la figura 6.3-I.3-Ψ.1 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 1.0 x 1010 s. y los distintos valores del factor de permeabilidad K (K1 = 1x10-13 m2,
K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
En el caso del mayor valor de K, la función de corriente es igual a cero excepto en la zona situada
junto al pozo y a la vertical del foco térmico, donde las isolíneas tienen forma ovoidal. En el segundo
mayor valor de K se inicia la forma ovoidal de las isolíneas. En los otros casos, las isolíneas de valores
mayores están en la parte superior izquierda del dominio.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
10
Figura 6.3-I.3-𝚿.1 Isolíneas de la función de corriente para un intervalo de tiempo 1.0x10 s en el dominio de estudio.
117
En la figura 6.3-I.3-Ψ.2 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 3.0 x 1010 s. y los distintos valores del factor de permeabilidad K (K1 = 1x10-13 m2,
K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
En los dos casos de los mayores valores de K, la función de corriente es igual a cero excepto en
la zona situada junto al pozo y a la vertical del foco térmico, donde las isolíneas tienen dos distorsiones
con forma ovoidal. En los otros casos, las isolíneas de valores mayores están en la parte superior
izquierda del dominio.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
10
Figura 6.3-I.3-𝚿.2 Isolíneas de la función de corriente para un intervalo de tiempo 3.0x10 s en el dominio de estudio.
118
En la figura 6.3-I.3-Ψ.3 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 9.0 x 1010 s. y los distintos valores del factor de permeabilidad K (K1 = 1x10-13 m2,
K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
En los dos casos de los mayores valores de K, la función de corriente es igual a cero excepto en
la zona situada junto a la vertical del foco térmico y junto al pozo, donde las isolíneas tienen dos formas
ovoidales con distintos valores en cada forma ovoidal. En los otros casos, las isolíneas de valores
mayores están en la parte superior izquierda del dominio y en el tercer valor de K se inicia una forma
ovoidal de las isolíneas.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
10
Figura 6.3-I.3-𝚿.3 Isolíneas de la función de corriente para un intervalo de tiempo 9.0x10 s en el dominio de estudio.
119
En la figura 6.3-I.3-Ψ.4 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 1.2 x 1011 s. y los distintos valores del factor de permeabilidad K (K1 = 1x10-13 m2,
K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
11
Figura 6.3-I.3-𝚿.4 Isolíneas de la función de corriente para un intervalo de tiempo 1.2x10 s en el dominio de estudio.
120
En la figura 6.3-I.3-Ψ.5 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 3.2 x 1011 s. y los distintos valores del factor de permeabilidad K (K1 = 1x10-13 m2,
K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
Todos los casos siguen el mismo comportamiento del caso anterior, iniciándose en el cuarto
valor de K la forma ovoidal de las isolíneas.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
11
Figura 6.3-I.3-𝚿.5 Isolíneas de la función de corriente para un intervalo de tiempo 3.2x10 s en el dominio de estudio.
121
En la figura 6.3-I.3-Ψ.6 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 5.0 x 1011 s. y los distintos valores del factor de permeabilidad K (K1 = 1x10-13 m2,
K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
11
Figura 6.3-I.3-𝚿.6 Isolíneas de la función de corriente para un intervalo de tiempo 5.0x10 s en el dominio de estudio.
122
CASO I.3-𝚿: Influencia del factor de permeabilidad K
123
6.3_I.4 Configuración de análisis 4: Influencia del factor de difusividad D
En la figura 6.3-I.4-𝑇.1 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 1.0 x 1010 s. y los distintos valores del factor de difusividad D
(D1 = 1x10-5 m2/s, D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
Se aprecia la evolución de las isolíneas hacia la superficie del dominio, en función de cuanto
mayor es el factor de difusividad D, más cerca de la superficie están las isolíneas.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
10
Figura 6.3-I.4-𝑻.1 Isolíneas de la temperatura adimensional para un intervalo de tiempo 1.0x10 s en el dominio de estudio.
124
En la figura 6.3-I.4-𝑇.2 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 3.0 x 1010 s. y los distintos valores del factor de difusividad D
(D1 = 1x10-5 m2/s, D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
10
Figura 6.3-I.4-𝑻.2 Isolíneas de la temperatura adimensional para un intervalo de tiempo 3.0x10 s en el dominio de estudio.
125
En la figura 6.3-I.4-𝑇.3 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 9.0 x 1010 s. y los distintos valores del factor de difusividad D
(D1 = 1x10-5 m2/s, D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
La evolución de las isolíneas sigue el mismo comportamiento del caso anterior, perdiéndose la
forma horizontal de las isolíneas.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
10
Figura 6.3-I.4-𝑻.3 Isolíneas de la temperatura adimensional para un intervalo de tiempo 9.0x10 s en el dominio de estudio.
126
En la figura 6.3-I.4-𝑇.4 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 1.2 x 1011 s. y los distintos valores del factor de difusividad D
(D1 = 1x10-5 m2/s, D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
La evolución de las isolíneas sigue el mismo comportamiento del caso anterior. Se aprecia el
inicio de la curvatura de las isolíneas.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
11
Figura 6.3-I.4-𝑻.4 Isolíneas de la temperatura adimensional para un intervalo de tiempo 1.2x10 s en el dominio de estudio.
127
En la figura 6.3-I.4-𝑇.5 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 3.2 x 1011 s. y los distintos valores del factor de difusividad D
(D1 = 1x10-5 m2/s, D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
La evolución de las isolíneas sigue el mismo comportamiento del caso anterior. Se aprecia, en el
último caso, la forma curvilínea de las isolíneas.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
11
Figura 6.3-I.4-𝑻.5 Isolíneas de la temperatura adimensional para un intervalo de tiempo 3.2x10 s en el dominio de estudio.
128
En la figura 6.3-I.4-𝑇.6 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 5.0 x 1011 s. y los distintos valores del factor de difusividad D
(D1 = 1x10-5 m2/s, D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
11
Figura 6.3-I.4-𝑻.6 Isolíneas de la temperatura adimensional para un intervalo de tiempo 5.0x10 s en el dominio de estudio.
129
CASO I.4-𝑻: Influencia del factor de difusividad D
130
En segundo lugar, vamos a estudiar la distribución de flujo de corriente Ψ en el dominio
geotérmico.
En la figura 6.3-I.4-Ψ.1 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 1.0 x 1010 s. y los distintos valores del factor de difusividad D (D1 = 1x10-5 m2/s,
D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
Para los valores mayores del factor de difusividad D, se aprecia la formación de isolíneas
ovoidales junto a la vertical del foco térmico y a la derecha del pozo. También se aprecian isolíneas de
valores mayores en la parte superior izquierda del dominio.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
10
Figura 6.3-I.4-𝚿.1 Isolíneas de la función de corriente para un intervalo de tiempo 1.0x10 s en el dominio de estudio.
131
En la figura 6.3-I.4-Ψ.2 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 3.0 x 1010 s. y los distintos valores del factor de difusividad D (D1 = 1x10-5 m2/s,
D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
Todos los casos siguen el mismo comportamiento anterior. Se van desarrollando las isolíneas
ovoidales.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
10
Figura 6.3-I.4-𝚿.2 Isolíneas de la función de corriente para un intervalo de tiempo 3.0x10 s en el dominio de estudio.
132
En la figura 6.3-I.4-Ψ.3 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 9.0 x 1010 s. y los distintos valores del factor de difusividad D (D1 = 1x10-5 m2/s,
D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
Todos los casos siguen el comportamiento de los casos anteriores, con las isolíneas ovoidales en
todos los casos.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
10
Figura 6.3-I.4-𝚿.3 Isolíneas de la función de corriente para un intervalo de tiempo 9.0x10 s en el dominio de estudio.
133
En la figura 6.3-I.4-Ψ.4 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 1.2 x 1011 s. y los distintos valores del factor de difusividad D (D1 = 1x10-5 m2/s,
D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
11
Figura 6.3-I.4-𝚿.4 Isolíneas de la función de corriente para un intervalo de tiempo 1.2x10 s en el dominio de estudio.
134
En la figura 6.3-I.4-Ψ.5 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 3.2 x 1011 s. y los distintos valores del factor de difusividad D (D1 = 1x10-5 m2/s,
D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
11
Figura 6.3-I.4-𝚿.5 Isolíneas de la función de corriente para un intervalo de tiempo 3.2x10 s en el dominio de estudio.
135
En la figura 6.3-I.4-Ψ.6 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 5.0 x 1011 s. y los distintos valores del factor de difusividad D (D1 = 1x10-5 m2/s,
D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
11
Figura 6.3-I.4-𝚿.6 Isolíneas de la función de corriente para un intervalo de tiempo 5.0x10 s en el dominio de estudio.
136
CASO I.4-𝚿: Influencia del factor de difusividad D
137
CASO II
En la figura 6.3-II.1-𝑇.1 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 1.0 x 1010 s. y los distintos tamaños del foco térmico considerados (5, 10, 15
y 20 celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
Podemos concluir que para el intervalo de tiempo considerado, el tamaño del foco térmico sólo
influye en la distancia que las isolíneas ocupan en la dirección horizontal (eje 𝑥).
10
Figura 6.3-II.1-T.1 Isolíneas de la temperatura adimensional para un intervalo de tiempo 1.0x10 s en el dominio de estudio.
138
En la figura 6.3-II.1-𝑇.2 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 3.0 x 1010 s. y los distintos tamaños del foco térmico considerados (5, 10, 15
y 20 celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
Se puede apreciar la evolución de las isolíneas, con cierta curvatura, hacia la superficie del
dominio y el tamaño del foco térmico sólo influye en la distancia que las isolíneas ocupan en la dirección
horizontal (eje 𝑥).
10
Figura 6.3-II.1-T.2 Isolíneas de la temperatura adimensional para un intervalo de tiempo 3.0x10 s en el dominio de estudio.
139
En la figura 6.3-II.1-𝑇.3 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 9.0 x 1010 s. y los distintos tamaños del foco térmico considerados (5, 10, 15
y 20 celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
Las isolíneas siguen evolucionando hacia la superficie del dominio y tienen una mayor curvatura
con valores máximos. Cuando el foco térmico (15 y 20 celdas) alcanza la vertical del pozo, las isolíneas
de la temperatura a cada lado del pozo tienen distintos valores.
10
Figura 6.3-II.1-T.3 Isolíneas de la temperatura adimensional para un intervalo de tiempo 9.0x10 s en el dominio de estudio.
140
En la figura 6.3-II.1-𝑇.4 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 1.2 x 1011 s. y los distintos tamaños del foco térmico considerados (5, 10, 15
y 20 celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
11
Figura 6.3-II.1-T.4 Isolíneas de la temperatura adimensional para un intervalo de tiempo 1.2x10 s en el dominio de estudio.
141
En la figura 6.3-II.1-𝑇.5 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 3.2 x 1011 s. y los distintos tamaños del foco térmico considerados (5, 10, 15
y 20 celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
11
Figura 6.3-II.1-T.5 Isolíneas de la temperatura adimensional para un intervalo de tiempo 3.2x10 s en el dominio de estudio.
142
En la figura 6.3-II.1-𝑇.6 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 5.0 x 1011 s. y los distintos tamaños del foco térmico considerados (5, 10, 15
y 20 celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
11
Figura 6.3-II.1-T.6 Isolíneas de la temperatura adimensional para un intervalo de tiempo 5.0x10 s en el dominio de estudio.
143
CASO II.1-𝑻: Influencia del tamaño del foco térmico
144
En segundo lugar, vamos a estudiar la distribución de flujo de corriente Ψ en el dominio
geotérmico.
En la figura 6.3-II.1-Ψ.1 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 1.0 x 1010 s. y los distintos tamaños del foco térmico considerados (5, 10, 15 y 20
celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
Se puede apreciar la influencia del foco al desplazarse hacia la derecha los vértices inferiores de
las isolíneas con la ampliación del tamaño del foco térmico hacia la derecha. Se puede apreciar, como
diferencia con el caso I, que las isolíneas a la derecha del pozo tienen valores menores que en el
6.3-I.1-Ψ.1
10
Figura 6.3-II.1-𝚿.1 Isolíneas de la función de corriente para un intervalo de tiempo 1.0x10 s en el dominio de estudio.
145
En la figura 6.3-II.1-Ψ.2 se presentan las gráficas con las isolíneas de la función de corriente
para un intervalo de tiempo 3.0 x 1010 s. y los distintos tamaños del foco térmico considerados (5, 10, 15
y 20 celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
10
Figura 6.3-II.1-𝚿.2 Isolíneas de la función de corriente para un intervalo de tiempo 3.0x10 s en el dominio de estudio.
146
En la figura 6.3-II.1-Ψ.3 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 9.0 x 1010 s. y los distintos tamaños del foco térmico considerados (5, 10, 15 y 20
celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
Se desarrolla la formación de isolíneas ovoidales en la zona situada sobre la vertical del foco
térmico, que se va desplazando hacia la derecha del dominio. Se mantiene la diferencia indicada en los
casos anteriores con el caso I.
10
Figura 6.3-II.1-𝚿.3 Isolíneas de la función de corriente para un intervalo de tiempo 9.0x10 s en el dominio de estudio.
147
En la figura 6.3-II.1-Ψ.4 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 1.2 x 1011 s. y los distintos tamaños del foco térmico considerados (5, 10, 15 y 20
celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
11
Figura 6.3-II.1-𝚿.4 Isolíneas de la función de corriente para un intervalo de tiempo 1.2x10 s en el dominio de estudio.
148
En la figura 6.3-II.1-Ψ.5 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 3.2 x 1011 s. y los distintos tamaños del foco térmico considerados (5, 10, 15 y 20
celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
Se desarrolla la formación de isolíneas ovoidales muy parecido al caso anterior. En todos los
casos aumenta el valor de las isolíneas más interiores. Se mantiene la diferencia indicada en los casos
anteriores con el caso I.
11
Figura 6.3-II.1-𝚿.5 Isolíneas de la función de corriente para un intervalo de tiempo 3.2x10 s en el dominio de estudio.
149
En la figura 6.3-II.1-Ψ.6 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 5.0 x 1011 s. y los distintos tamaños del foco térmico considerados (5, 10, 15 y 20
celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
Las formas y valores de las isolíneas parecen haber llegado a su máximo desarrollo, no hay
grandes cambios con respecto al intervalo del caso anterior. Se mantiene la diferencia indicada en los
casos anteriores con el caso I.
11
Figura 6.3-II.1-𝚿.6 Isolíneas de la función de corriente para un intervalo de tiempo 5.0x10 s en el dominio de estudio.
150
CASO II.1-𝚿: Influencia del tamaño del foco térmico
151
6.3_II.2 Configuración de análisis 2: Influencia de la profundidad del pozo.
En la figura 6.3-II.2-𝑇.1 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 1.0 x 1010 s. y las distintas profundidades del pozo consideradas (200 m,
400 m, 600 m y 800 m) en todo el dominio de estudio.
Podemos concluir que para el intervalo de tiempo considerado, la profundidad del pozo no
influye sobre las isolíneas, sólo en el pozo de 800 m las isolíneas pierden su forma horizontal.
10
Figura 6.3-II.2-T.1 Isolíneas de la temperatura adimensional para un intervalo de tiempo 1.0x10 s en el dominio de estudio.
152
En la figura 6.3-II.2-𝑇.2 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 3.0 x 1010 s. y las distintas profundidades del pozo consideradas (200 m,
400 m, 600 m y 800 m) en todo el dominio de estudio.
Las isolíneas evolucionan hacia la superficie del dominio, por tanto la presencia del pozo
modifica la forma horizontal de las isolíneas.
10
Figura 6.3-II.2-T.2 Isolíneas de la temperatura adimensional para un intervalo de tiempo 3.0x10 s en el dominio de estudio.
153
En la figura 6.3-II.2-𝑇.3 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 9.0 x 1010 s. y las distintas profundidades del pozo consideradas (200 m,
400 m, 600 m y 800 m) en todo el dominio de estudio.
Para el intervalo de tiempo considerado, las isolíneas siguen evolucionando hacia la superficie
del dominio y se inicia la forma curvilínea.
10
Figura 6.3-II.2-T.3 Isolíneas de la temperatura adimensional para un intervalo de tiempo 9.0x10 s en el dominio de estudio.
154
En la figura 6.3-II.2-𝑇.4 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 1.2 x 1011 s. y las distintas profundidades del pozo consideradas (200 m,
400 m, 600 m y 800 m) en todo el dominio de estudio.
11
Figura 6.3-II.2-T.4 Isolíneas de la temperatura adimensional para un intervalo de tiempo 1.2x10 s en el dominio de estudio.
155
En la figura 6.3-II.2-𝑇.5 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 3.2 x 1011 s. y las distintas profundidades del pozo consideradas (200 m,
400 m, 600 m y 800 m) en todo el dominio de estudio.
11
Figura 6.3-II.2-T.5 Isolíneas de la temperatura adimensional para un intervalo de tiempo 3.2x10 s en el dominio de estudio.
156
En la figura 6.3-II.2-𝑇.6 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 5.0 x 1011 s. y las distintas profundidades del pozo consideradas (200 m,
400 m, 600 m y 800 m) en todo el dominio de estudio.
11
Figura 6.3-II.2-T.6 Isolíneas de la temperatura adimensional para un intervalo de tiempo 5.0x10 s en el dominio de estudio.
157
CASO II.2-𝑻: Influencia de la profundidad del pozo
158
En segundo lugar, vamos a estudiar la distribución de flujo de corriente Ψ en el dominio
geotérmico.
En la figura 6.3-II.2-Ψ.1 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 1.0 x 1010 s. y las distintas profundidades del pozo consideradas (200 m, 400 m,
600 m y 800 m) en todo el dominio de estudio.
Se puede apreciar la influencia del pozo al encontrase isolíneas de valores mayores hacia la
parte superior del dominio y a la izquierda del pozo. Encontramos dos distorsiones de la forma de las
curvas, una sobre la vertical del foco térmico y otra en la posición donde se encuentra el pozo en cada
caso. Cuando más profundo es el pozo, las isolíneas de corriente son menos horizontales. Como
diferencia con el caso I, las isolíneas a la derecha del pozo tienen valores menores que en el 6.3-I.2-Ψ.1.
10
Figura 6.3-II.2-𝚿.1 Isolíneas de la función de corriente para un intervalo de tiempo 1.0x10 s en el dominio de estudio.
159
En la figura 6.3-II.2-Ψ.2 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 3.0 x 1010 s. y las distintas profundidades del pozo consideradas (200 m, 400 m,
600 m y 800 m) en todo el dominio de estudio.
10
Figura 6.3-II.2-𝚿.2 Isolíneas de la función de corriente para un intervalo de tiempo 3.0x10 s en el dominio de estudio.
160
En la figura 6.3-II.2-Ψ.3 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 9.0 x 1010 s. y las distintas profundidades del pozo consideradas (200 m, 400 m,
600 m y 800 m) en todo el dominio de estudio.
La forma ovoidal de las isolíneas se desarrolla en todos los casos. La diferencia está en la zona
alrededor de la parte inferior del pozo.
10
Figura 6.3-II.2-𝚿.3 Isolíneas de la función de corriente para un intervalo de tiempo 9.0x10 s en el dominio de estudio.
161
En la figura 6.3-II.2-Ψ.4 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 1.2 x 1011 s. y las distintas profundidades del pozo consideradas (200 m, 400 m,
600 m y 800 m) en todo el dominio de estudio.
11
Figura 6.3-II.2-𝚿.4 Isolíneas de la función de corriente para un intervalo de tiempo 1.2x10 s en el dominio de estudio.
162
En la figura 6.3-II.2-Ψ.5 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 3.2 x 1011 s. y las distintas profundidades del pozo consideradas (200 m, 400 m,
600 m y 800 m) en todo el dominio de estudio.
11
Figura 6.3-II.2-𝚿.5 Isolíneas de la función de corriente para un intervalo de tiempo 3.2x10 s en el dominio de estudio.
163
En la figura 6.3-II.2-Ψ.6 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 5.0 x 1011 s. y las distintas profundidades del pozo consideradas (200 m, 400 m,
600 m y 800 m) en todo el dominio de estudio.
Las formas y valores de las isolíneas parecen haber llegado a su máximo desarrollo, no hay
grandes cambios con respecto al intervalo de tiempo anterior.
11
Figura 6.3-II.2-𝚿.6 Isolíneas de la función de corriente para un intervalo de tiempo 5.0x10 s en el dominio de estudio.
164
CASO II.2-𝚿: Influencia de la profundidad del pozo
165
6.3_II.3 Configuración de análisis 3: Influencia del factor de permeabilidad K
En la figura 6.3-II.3-𝑇.1 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 1.0 x 1010 s. y los distintos valores del factor de permeabilidad K
(K1 = 1x10-13 m2, K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
Se aprecia la evolución de las isolíneas hacia la parte superior del dominio y la forma ovoidal
para el mayor valor del factor de permeabilidad K.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
10
Figura 6.3-II.3-𝑻.1 Isolíneas de la temperatura adimensional para un intervalo de tiempo 1.0x10 s en el dominio de estudio.
166
En la figura 6.3-II.3-𝑇.2 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 3.0 x 1010 s. y los distintos valores del factor de permeabilidad K
(K1 = 1x10-13 m2, K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
Aumenta la evolución de las isolíneas hacia la parte superior del dominio y la forma ovoidal para
los dos mayores valores del factor de permeabilidad K.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
10
Figura 6.3-II.3-𝑻.2 Isolíneas de la temperatura adimensional para un intervalo de tiempo 3.0x10 s en el dominio de estudio.
167
En la figura 6.3-II.3-𝑇.3 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 9.0 x 1010 s. y los distintos valores del factor de permeabilidad K
(K1 = 1x10-13 m2, K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
10
Figura 6.3-II.3-𝑻.3 Isolíneas de la temperatura adimensional para un intervalo de tiempo 9.0x10 s en el dominio de estudio.
168
En la figura 6.3-II.3-𝑇.4 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 1.2 x 1011 s. y los distintos valores del factor de permeabilidad K
(K1 = 1x10-13 m2, K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
11
Figura 6.3-II.3-𝑻.4 Isolíneas de la temperatura adimensional para un intervalo de tiempo 1.2x10 s en el dominio de estudio.
169
En la figura 6.3-II.3-𝑇.5 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 3.2 x 1011 s. y los distintos valores del factor de permeabilidad K
(K1 = 1x10-13 m2, K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
11
Figura 6.3-II.3-𝑻.5 Isolíneas de la temperatura adimensional para un intervalo de tiempo 3.2x10 s en el dominio de estudio.
170
En la figura 6.3-II.3-𝑇.6 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 5.0 x 1011 s. y los distintos valores del factor de permeabilidad K
(K1 = 1x10-13 m2, K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
11
Figura 6.3-II.3-𝑻.6 Isolíneas de la temperatura adimensional para un intervalo de tiempo 5.0x10 s en el dominio de estudio.
171
CASO II.3-𝑻: Influencia del factor de permeabilidad K
172
En segundo lugar, vamos a estudiar la distribución de flujo de corriente Ψ en el dominio
geotérmico.
En la figura 6.3-II.3-Ψ.1 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 1.0 x 1010 s. y los distintos valores del factor de permeabilidad K (K1 = 1x10-13 m2,
K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
En el caso del mayor valor de K, la función de corriente es igual a cero excepto en la zona situada
junto al pozo y sobre la vertical del foco térmico, donde las isolíneas tienen forma ovoidal. En los otros
casos, las isolíneas de valores mayores están en la parte superior izquierda del dominio. Como diferencia
con el caso I, las isolíneas a la derecha del pozo tienen valores menores que en el 6.3-I.3-Ψ.1.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
10
Figura 6.3-II.3-𝚿.1 Isolíneas de la función de corriente para un intervalo de tiempo 1.0x10 s en el dominio de estudio.
173
En la figura 6.3-II.3-Ψ.2 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 3.0 x 1010 s. y los distintos valores del factor de permeabilidad K (K1 = 1x10-13 m2,
K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
En los dos casos de los mayores valores de K, la función de corriente es igual a cero excepto en
la zona situada junto al pozo y sobre la vertical del foco térmico, donde las isolíneas tienen dos
distorsiones con forma ovoidal. En los otros casos, las isolíneas de valores mayores están en la parte
superior izquierda del dominio.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
10
Figura 6.3-II.3-𝚿.2 Isolíneas de la función de corriente para un intervalo de tiempo 3.0x10 s en el dominio de estudio.
174
En la figura 6.3-II.3-Ψ.3 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 9.0 x 1010 s. y los distintos valores del factor de permeabilidad K (K1 = 1x10-13 m2,
K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
En los dos casos de los mayores valores de K, la función de corriente es igual a cero excepto en
la zona situada sobre la vertical del foco térmico y junto al pozo, donde las isolíneas tienen dos formas
ovoidales con distintos valores en cada forma ovoidal. En los otros casos, las isolíneas de valores
mayores están en la parte superior izquierda del dominio y en el tercer valor de K se inicia una forma
ovoidal.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
10
Figura 6.3-II.3-𝚿.3 Isolíneas de la función de corriente para un intervalo de tiempo 9.0x10 s en el dominio de estudio.
175
En la figura 6.3-II.3-Ψ.4 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 1.2 x 1011 s. y los distintos valores del factor de permeabilidad K (K1 = 1x10-13 m2,
K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
11
Figura 6.3-II.3-𝚿.4 Isolíneas de la función de corriente para un intervalo de tiempo 1.2x10 s en el dominio de estudio.
176
En la figura 6.3-II.3-Ψ.5 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 3.2 x 1011 s. y los distintos valores del factor de permeabilidad K (K1 = 1x10-13 m2,
K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
Todos los casos siguen el mismo comportamiento del caso anterior, iniciándose en el cuarto
valor de K las isolíneas de forma ovoidal.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
11
Figura 6.3-II.3-𝚿.5 Isolíneas de la función de corriente para un intervalo de tiempo 3.2x10 s en el dominio de estudio.
177
En la figura 6.3-II.3-Ψ.6 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 5.0 x 1011 s. y los distintos valores del factor de permeabilidad K (K1 = 1x10-13 m2,
K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
11
Figura 6.3-II.3-𝚿.6 Isolíneas de la función de corriente para un intervalo de tiempo 5.0x10 s en el dominio de estudio.
178
CASO II.3-𝚿: Influencia del factor de permeabilidad K
179
6.3_II.4 Configuración de análisis 4: Influencia del factor de difusividad D
En la figura 6.3-II.4-𝑇.1 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 1.0 x 1010 s. y los distintos valores del factor de difusividad D
(D1 = 1x10-5 m2/s, D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
Se aprecia la evolución de las isolíneas hacia la superficie del dominio, en función de cuanto
mayor es el factor de difusividad D, más cerca de la superficie están las isolíneas.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
10
Figura 6.3-II.4-𝑻.1 Isolíneas de la temperatura adimensional para un intervalo de tiempo 1.0x10 s en el dominio de estudio.
180
En la figura 6.3-II.4-𝑇.2 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 3.0 x 1010 s. y los distintos valores del factor de difusividad D
(D1 = 1x10-5 m2/s, D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
10
Figura 6.3-II.4-𝑻.2 Isolíneas de la temperatura adimensional para un intervalo de tiempo 3.0x10 s en el dominio de estudio.
181
En la figura 6.3-II.4-𝑇.3 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 9.0 x 1010 s. y los distintos valores del factor de difusividad D
(D1 = 1x10-5 m2/s, D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
La evolución de las isolíneas sigue el mismo comportamiento del caso anterior, perdiéndose la
forma horizontal de las isolíneas.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
10
Figura 6.3-II.4-𝑻.3 Isolíneas de la temperatura adimensional para un intervalo de tiempo 9.0x10 s en el dominio de estudio.
182
En la figura 6.3-II.4-𝑇.4 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 1.2 x 1011 s. y los distintos valores del factor de difusividad D
(D1 = 1x10-5 m2/s, D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
La evolución de las isolíneas sigue el mismo comportamiento del caso anterior. Se aprecia el
inicio de la curvatura de las isolineas.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
11
Figura 6.3-II.4-𝑻.4 Isolíneas de la temperatura adimensional para un intervalo de tiempo 1.2x10 s en el dominio de estudio.
183
En la figura 6.3-II.4-𝑇.5 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 3.2 x 1011 s. y los distintos valores del factor de difusividad D
(D1 = 1x10-5 m2/s, D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
La evolución de las isolíneas sigue el mismo comportamiento del caso anterior. Se aprecia, en el
último caso, la forma ovoidal de las isolíneas.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
11
Figura 6.3-II.4-𝑻.5 Isolíneas de la temperatura adimensional para un intervalo de tiempo 3.2x10 s en el dominio de estudio.
184
En la figura 6.3-II.4-𝑇.6 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 5.0 x 1011 s. y los distintos valores del factor de difusividad D
(D1 = 1x10-5 m2/s, D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
11
Figura 6.3-II.4-𝑻.6 Isolíneas de la temperatura adimensional para un intervalo de tiempo 5.0x10 s en el dominio de estudio.
185
CASO II.4-𝑻: Influencia del factor de difusividad D
186
En segundo lugar, vamos a estudiar la distribución de flujo de corriente Ψ en el dominio
geotérmico.
En la figura 6.3-II.4-Ψ.1 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 1.0 x 1010 s. y los distintos valores del factor de difusividad D (D1 = 1x10-5 m2/s,
D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
Para los valores mayores del factor de difusividad D, se aprecia la formación de isolíneas
ovoidales junto a la vertical del foco térmico y a la derecha del pozo. También se aprecian isolíneas de
valores mayores en la parte superior izquierda del dominio. Como diferencia con el caso I, las isolíneas a
la derecha del pozo tienen valores menores que en el 6.3-I.4-Ψ.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
10
Figura 6.3-II.4-𝚿.1 Isolíneas de la función de corriente para un intervalo de tiempo 1.0x10 s en el dominio de estudio.
187
En la figura 6.3-II.4-Ψ.2 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 3.0 x 1010 s. y los distintos valores del factor de difusividad D (D1 = 1x10-5 m2/s,
D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
Todos los casos siguen el mismo comportamiento anterior. Se van desarrollando las isolíneas
ovoidales.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
10
Figura 6.3-II.4-𝚿.2 Isolíneas de la función de corriente para un intervalo de tiempo 3.0x10 s en el dominio de estudio.
188
En la figura 6.3-II.4-Ψ.3 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 9.0 x 1010 s. y los distintos valores del factor de difusividad D (D1 = 1x10-5 m2/s,
D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
Todos los casos siguen el mismo comportamiento anterior, con las isolíneas ovoidales en todos
los casos.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
10
Figura 6.3-II.4-𝚿.3 Isolíneas de la función de corriente para un intervalo de tiempo 9.0x10 s en el dominio de estudio.
189
En la figura 6.3-II.4-Ψ.4 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 1.2 x 1011 s. y los distintos valores del factor de difusividad D (D1 = 1x10-5 m2/s,
D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
11
Figura 6.3-II.4-𝚿.4 Isolíneas de la función de corriente para un intervalo de tiempo 1.2x10 s en el dominio de estudio.
190
En la figura 6.3-II.4-Ψ.5 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 3.2 x 1011 s. y los distintos valores del factor de difusividad D (D1 = 1x10-5 m2/s,
D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
11
Figura 6.3-II.4-𝚿.5 Isolíneas de la función de corriente para un intervalo de tiempo 3.2x10 s en el dominio de estudio.
191
En la figura 6.3-II.4-Ψ.6 se presentan las gráficas con las isolíneas de la función de corriente para
un intervalo de tiempo 5.0 x 1011 s. y los distintos valores del factor de difusividad D (D1 = 1x10-5 m2/s,
D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
11
Figura 6.3-II.4-𝚿.6 Isolíneas de la función de corriente para un intervalo de tiempo 5.0x10 s en el dominio de estudio.
192
CASO II.4-𝚿: Influencia del factor de difusividad D
193
CASO III
En la figura 6.3-III.1-𝑇.1 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 1.0 x 1010 s. y los distintos tamaños del foco térmico considerados (5, 10, 15
y 20 celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
Para el intervalo de tiempo considerado, el tamaño del foco térmico sólo influye en la distancia
que las isolíneas ocupan en la dirección horizontal (eje χ).
10
Figura 6.3-III.1-T.1 Isolíneas de la temperatura adimensional para un intervalo de tiempo 1.0x10 s en el dominio de estudio.
194
En la figura 6.3-III.1-𝑇.2 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 3.0 x 1010 s. y los distintos tamaños del foco térmico considerados (5, 10, 15
y 20 celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
Se puede apreciar la evolución de las isolíneas, con cierta curvatura, hacia la superficie del
dominio y el tamaño del foco térmico sólo influye en la distancia que las isolíneas ocupan en la dirección
horizontal (eje 𝑥).
10
Figura 6.3-III.1-T.2 Isolíneas de la temperatura adimensional para un intervalo de tiempo 3.0x10 s en el dominio de estudio.
195
En la figura 6.3-III.1-𝑇.3 se presentan las gráficas con las isolíneas de temperatura adimensional
para unintervalo de tiempo 9.0 x 1010 s. y los distintos tamaños del foco térmico considerados (5, 10, 15
y 20 celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
Las isolíneas siguen evolucionando hacia la superficie del dominio y tienen una mayor curvatura
con valores máximos. Cuando el foco térmico (15 y 20 celdas) alcanza la vertical del pozo, las isolíneas
de la temperatura a cada lado del pozo tienen distintos valores.
10
Figura 6.3-III.1-T.3 Isolíneas de la temperatura adimensional para un intervalo de tiempo 9.0x10 s en el dominio de estudio.
196
En la figura 6.3-III.1-𝑇.4 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 1.2 x 1011 s. y los distintos tamaños del foco térmico considerados (5, 10, 15
y 20 celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
11
Figura 6.3-III.1-T.4 Isolíneas de la temperatura adimensional para un intervalo de tiempo 1.2x10 s en el dominio de estudio.
197
En la figura 6.3-III.1-𝑇.5 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 3.2 x 1011 s.y los distintos tamaños del foco térmico considerados (5, 10, 15
y 20 celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
11
Figura 6.3-III.1-T.5 Isolíneas de la temperatura adimensional para un intervalo de tiempo 3.2x10 s en el dominio de estudio.
198
En la figura 6.3-III.1-𝑇.6 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 5.0 x 1011 s. y los distintos tamaños del foco térmico considerados (5, 10, 15
y 20 celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
11
Figura 6.3-III.1-T.6 Isolíneas de la temperatura adimensional para un intervalo de tiempo 5.0x10 s en el dominio de estudio.
199
CASO III.1-𝑻: Influencia del tamaño del foco térmico
200
En segundo lugar, vamos a estudiar la distribución de flujo de corriente Ψ en el dominio
geotérmico.
En la figura 6.3-III.1-Ψ.1 se presentan las gráficas con las isolíneas de la función de corriente
para un intervalo de tiempo 1.0 x 1010 s. y los distintos tamaños del foco térmico considerados (5, 10, 15
y 20 celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
Se puede apreciar la influencia del foco al desplazarse hacia la derecha los vértices inferiores de
las isolíneas con la ampliación del tamaño del foco térmico hacia la derecha. También como diferencia
con el CASO I y II, en este CASO III encontramos isolíneas de valores mayores hacia la parte superior
derecha del dominio y a la derecha del pozo ya que la función de corriente en esta zona es Ψ = 6.7x10-6
10
Figura 6.3-III.1-𝚿.1 Isolíneas de la función de corriente para un intervalo de tiempo 1.0x10 s en el dominio de estudio.
201
En la figura 6.3-III.1-Ψ.2 se presentan las gráficas con las isolíneas de la función de corriente
para un intervalo de tiempo 3.0 x 1010 s. y los distintos tamaños del foco térmico considerados (5, 10, 15
y 20 celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
10
Figura 6.3-III.1-𝚿.2 Isolíneas de la función de corriente para un intervalo de tiempo 3.0x10 s en el dominio de estudio.
202
En la figura 6.3-III.1-Ψ.3 se presentan las gráficas con las isolíneas de la función de corriente
para un intervalo de tiempo 9.0 x 1010 s. y los distintos tamaños del foco térmico considerados (5, 10, 15
y 20 celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
Se desarrolla la formación de isolíneas ovoidales en la zona situada sobre la vertical del foco
térmico, que se va desplazando hacia la derecha del dominio con la ampliación del tamaño del foco
térmico hacia la derecha.
10
Figura 6.3-III.1-𝚿.3 Isolíneas de la función de corriente para un intervalo de tiempo 9.0x10 s en el dominio de estudio.
203
En la figura 6.3-III.1-Ψ.4 se presentan las gráficas con las isolíneas de la función de corriente
para un intervalo de tiempo 1.2 x 1011 s. y los distintos tamaños del foco térmico considerados (5, 10, 15
y 20 celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
11
Figura 6.3-III.1-𝚿.4 Isolíneas de la función de corriente para un intervalo de tiempo 1.2x10 s en el dominio de estudio.
204
En la figura 6.3-III.1-Ψ.5 se presentan las gráficas con las isolíneas de la función de corriente
para un intervalo de tiempo 3.2 x 1011 s. y los distintos tamaños del foco térmico considerados (5, 10, 15
y 20 celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
Se desarrolla la formación de isolíneas ovoidales muy parecido al caso anterior. En todos los
casos aumenta el valor de las isolíneas más interiores.
11
Figura 6.3-III.1-𝚿.5 Isolíneas de la función de corriente para un intervalo de tiempo 3.2x10 s en el dominio de estudio.
205
En la figura 6.3-III.1-Ψ.6 se presentan las gráficas con las isolíneas de la función de corriente
para un intervalo de tiempo 5.0 x 1011 s. y los distintos tamaños del foco térmico considerados (5, 10, 15
y 20 celdas de tamaño del foco en la parte inferior izquierda) en todo el dominio de estudio.
Las formas y valores de las isolíneas parecen haber llegado a su máximo desarrollo, no hay
grandes cambios con respecto al intervalo del caso anterior.
11
Figura 6.3-III.1-𝚿.6 Isolíneas de la función de corriente para un intervalo de tiempo 5.0x10 s en el dominio de estudio.
206
CASO III.1-𝚿: Influencia del tamaño del foco térmico
207
6.3_III.2 Configuración de análisis 2: Influencia delaprofundidad del pozo.
En la figura 6.3-III.2-𝑇.1 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 1.0 x 1010 s. y las distintas profundidades del pozo consideradas (200 m,
400 m, 600 m y 800 m) en todo el dominio de estudio.
Podemos concluir que para el intervalo de tiempo considerado, la profundidad del pozo no
influye sobre las isolíneas, sólo en el pozo de 800 m las isolíneas pierden su forma horizontal.
10
Figura 6.3-III.2-T.1 Isolíneas de la temperatura adimensional para un intervalo de tiempo 1.0x10 s en el dominio de estudio.
208
En la figura 6.3-III.2-𝑇.2 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 3.0 x 1010 s. y las distintas profundidades del pozo consideradas (200 m,
400 m, 600 m y 800 m) en todo el dominio de estudio.
Las isolíneas evolucionan hacia la superficie del dominio, por tanto la presencia del pozo
modifica la forma horizontal de las isolíneas.
10
Figura 6.3-III.2-T.2 Isolíneas de la temperatura adimensional para un intervalo de tiempo 3.0x10 s en el dominio de estudio.
209
En la figura 6.3-III.2-𝑇.3 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 9.0 x 1010 s. y las distintas profundidades del pozo consideradas (200 m,
400 m, 600 m y 800 m) en todo el dominio de estudio.
Las isolíneas siguen evolucionando hacia la superficie del dominio y se inicia la forma curvilínea.
10
Figura 6.3-III.2-T.3 Isolíneas de la temperatura adimensional para un intervalo de tiempo 9.0x10 s en el dominio de estudio.
210
En la figura 6.3-III.2-𝑇.4 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 1.2 x 1011 s. y las distintas profundidades del pozo consideradas (200 m,
400 m, 600 m y 800 m) en todo el dominio de estudio.
11
Figura 6.3-III.2-T.4 Isolíneas de la temperatura adimensional para un intervalo de tiempo 1.2x10 s en el dominio de estudio.
211
En la figura 6.3-III.2-𝑇.5 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 3.2 x 1011 s. y las distintas profundidades del pozo consideradas (200 m,
400 m, 600 m y 800 m) en todo el dominio de estudio.
11
Figura 6.3-III.2-T.5 Isolíneas de la temperatura adimensional para un intervalo de tiempo 3.2x10 s en el dominio de estudio.
212
En la figura 6.3-III.2-𝑇.6 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 5.0 x 1011 s. y las distintas profundidades del pozo consideradas (200 m,
400 m, 600 m y 800 m) en todo el dominio de estudio.
11
Figura 6.3-III.2-T.6 Isolíneas de la temperatura adimensional para un intervalo de tiempo 5.0x10 s en el dominio de estudio.
213
CASO III.2-𝑻: Influencia de la profundidad del pozo
214
En segundo lugar, vamos a estudiar la distribución de flujo de corriente Ψ en el dominio
geotérmico.
En la figura 6.3-III.2-Ψ.1 se presentan las gráficas con las isolíneas de la función de corriente
para un intervalo de tiempo 1.0 x 1010 s. y las distintas profundidades del pozo consideradas (200 m,
400 m, 600 m y 800 m) en todo el dominio de estudio.
Encontramos dos distorsiones de la forma de las curvas, una sobre la vertical del foco térmico y
otra en la posición del pozo en cada caso. Se puede apreciar que cuando más profundo es el pozo, las
isolíneas de corriente son menos horizontales. También como diferencia con el CASO I y II, en este CASO
III encontramos isolíneas de valores mayores hacia la parte superior derecha del dominio y a la derecha
del pozo ya que la función de corriente en esta zona es Ψ = 6.7x10-6.
10
Figura 6.3-III.2-𝚿.1 Isolíneas de la función de corriente para un intervalo de tiempo 1.0x10 s en el dominio de estudio.
215
En la figura 6.3-III.2-Ψ.2 se presentan las gráficas con las isolíneas de la función de corriente
para un intervalo de tiempo 3.0 x 1010 s. y las distintas profundidades del pozo consideradas (200 m,
400 m, 600 m y 800 m) en todo el dominio de estudio.
10
Figura 6.3-III.2-𝚿.2 Isolíneas de la función de corriente para un intervalo de tiempo 3.0x10 s en el dominio de estudio.
216
En la figura 6.3-III.2-Ψ.3 se presentan las gráficas con las isolíneas de la función de corriente
para un intervalo de tiempo 9.0 x 1010 s. y las distintas profundidades del pozo consideradas (200 m,
400 m, 600 m y 800 m) en todo el dominio de estudio.
10
Figura 6.3-III.2-𝚿.3 Isolíneas de la función de corriente para un intervalo de tiempo 9.0x10 s en el dominio de estudio.
217
En la figura 6.3-III.2-Ψ.4 se presentan las gráficas con las isolíneas de la función de corriente
para un intervalo de tiempo 1.2 x 1011 s. y las distintas profundidades del pozo consideradas (200 m,
400 m, 600 m y 800 m) en todo el dominio de estudio.
11
Figura 6.3-III.2-𝚿.4 Isolíneas de la función de corriente para un intervalo de tiempo 1.2x10 s en el dominio de estudio.
218
En la figura 6.3-III.2-Ψ.5 se presentan las gráficas con las isolíneas de la función de corriente
para un intervalo de tiempo 3.2 x 1011 s. y las distintas profundidades del pozo consideradas (200 m,
400 m, 600 m y 800 m) en todo el dominio de estudio.
11
Figura 6.3-III.2-𝚿.5 Isolíneas de la función de corriente para un intervalo de tiempo 3.2x10 s en el dominio de estudio.
219
En la figura 6.3-III.2-Ψ.6 se presentan las gráficas con las isolíneas de la función de corriente
para un intervalo de tiempo 5.0 x 1011 s. y las distintas profundidades del pozo consideradas (200 m,
400 m, 600 m y 800 m) en todo el dominio de estudio.
Las formas y valores de las isolíneas parecen haber llegado a su máximo desarrollo, pues no hay
grandes cambios respecto al intervalo anterior.
11
Figura 6.3-III.2-𝚿.6 Isolíneas de la función de corriente para un intervalo de tiempo 5.0x10 s en el dominio de estudio.
220
CASO III.2-𝚿: Influencia de la profundidad del pozo
221
6.3_III.3 Configuración de análisis 3: Influencia del factor de permeabilidad K
En la figura 6.3-III.3-𝑇.1 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 1.0 x 1010 s. y los distintos valores del factor de permeabilidad K
(K1 = 1x10-13 m2, K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
Se aprecia la evolución de las isolíneas hacia la parte superior del dominio y la forma ovoidal
para el mayor valor del factor de permeabilidad K.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
10
Figura 6.3-III.3-𝑻.1 Isolíneas de la temperatura adimensional para un intervalo de tiempo 1.0x10 s en el dominio de estudio.
222
En la figura 6.3-III.3-𝑇.2 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 3.0 x 1010 s. y los distintos valores del factor de permeabilidad K
(K1 = 1x10-13 m2, K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
Aumenta la evolución de las isolíneas hacia la parte superior del dominio y la forma ovoidal para
los dos mayores valores del factor de permeabilidad K.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
10
Figura 6.3-III.3-𝑻.2 Isolíneas de la temperatura adimensional para un intervalo de tiempo 3.0x10 s en el dominio de estudio.
223
En la figura 6.3-III.3-𝑇.3 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 9.0 x 1010 s. y los distintos valores del factor de permeabilidad K
(K1 = 1x10-13 m2, K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
10
Figura 6.3-III.3-𝑻.3 Isolíneas de la temperatura adimensional para un intervalo de tiempo 9.0x10 s en el dominio de estudio.
224
En la figura 6.3-III.3-𝑇.4 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 1.2 x 1011 s. y los distintos valores del factor de permeabilidad K
(K1 = 1x10-13 m2, K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
11
Figura 6.3-III.3-𝑻.4 Isolíneas de la temperatura adimensional para un intervalo de tiempo 1.2x10 s en el dominio de estudio.
225
En la figura 6.3-III.3-𝑇.5 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 3.2 x 1011 s. y los distintos valores del factor de permeabilidad K
(K1 = 1x10-13 m2, K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
11
Figura 6.3-III.3-𝑻.5 Isolíneas de la temperatura adimensional para un intervalo de tiempo 3.2x10 s en el dominio de estudio.
226
En la figura 6.3-III.3-𝑇.6 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 5.0 x 1011 s. y los distintos valores del factor de permeabilidad K
(K1 = 1x10-13 m2, K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
11
Figura 6.3-III.3-𝑻.6 Isolíneas de la temperatura adimensional para un intervalo de tiempo 5.0x10 s en el dominio de estudio.
227
CASO III.3-𝑻: Influencia del factor de permeabilidad K
228
En segundo lugar, vamos a estudiar la distribución de flujo de corriente Ψ en el dominio
geotérmico.
En la figura 6.3-III.3-Ψ.1 se presentan las gráficas con las isolíneas de la función de corriente
para un intervalo de tiempo 1.0 x 1010 s. y los distintos valores del factor de permeabilidad K
(K1 = 1x10-13 m2, K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
En los dos valores mayores valor de K se desarrolla la forma ovoidal de las isolíneas.
También como diferencia con el CASO I y II, en este CASO III encontramos isolíneas de valores mayores
hacia la parte superior derecha del dominio y a la derecha del pozo ya que la función de corriente en
esta zona es Ψ = 6.7x10-6
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
10
Figura 6.3-III.3-𝚿.1 Isolíneas de la función de corriente para un intervalo de tiempo 1.0x10 s en el dominio de estudio.
229
En la figura 6.3-III.3-Ψ.2 se presentan las gráficas con las isolíneas de la función de corriente
para un intervalo de tiempo 3.0 x 1010 s. y los distintos valores del factor de permeabilidad K
(K1 = 1x10-13 m2, K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
En el caso del mayor valor de K la función de corriente es igual a cero excepto en la zona situada
junto al pozo y a la vertical del foco térmico. Además en los dos casos de los mayores valores de K las
isolíneas tienen dos distorsiones con forma ovoidal. También como diferencia con el CASO I y II, en este
CASO III encontramos isolíneas de valores mayores hacia la parte superior derecha del dominio y a la
derecha del pozo ya que la función de corriente en esta zona es Ψ = 6.7x10-6
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
10
Figura 6.3-III.3-𝚿.2 Isolíneas de la función de corriente para un intervalo de tiempo 3.0x10 s en el dominio de estudio.
230
En la figura 6.3-III.3-Ψ.3 se presentan las gráficas con las isolíneas de la función de corriente
para un intervalo de tiempo 9.0 x 1010 s. y los distintos valores del factor de permeabilidad K
(K1 = 1x10-13 m2, K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
En los dos casos de los valores mayores de K las isolíneas tienen tres distorsiones con forma
ovoidal. También como diferencia con el CASO I y II, en este CASO III encontramos isolíneas de valores
mayores hacia la parte superior derecha del dominio y a la derecha del pozo ya que la función de
corriente en esta zona es Ψ = 6.7x10-6
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
10
Figura 6.3-III.3-𝚿.3 Isolíneas de la función de corriente para un intervalo de tiempo 9.0x10 s en el dominio de estudio.
231
En la figura 6.3-III.3-Ψ.4 se presentan las gráficas con las isolíneas de la función de corriente
para un intervalo de tiempo 1.2 x 1011 s. y los distintos valores del factor de permeabilidad K
(K1 = 1x10-13 m2, K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
Todos los casos siguen el mismo comportamiento del caso anterior, iniciándose en el tercer
valor de K la forma ovoidal de las isolíneas.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
11
Figura 6.3-III.3-𝚿.4 Isolíneas de la función de corriente para un intervalo de tiempo 1.2x10 s en el dominio de estudio.
232
En la figura 6.3-III.3-Ψ.5 se presentan las gráficas con las isolíneas de la función de corriente
para un intervalo de tiempo 3.2 x 1011 s. y los distintos valores del factor de permeabilidad K
(K1 = 1x10-13 m2, K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
11
Figura 6.3-III.3-𝚿.5 Isolíneas de la función de corriente para un intervalo de tiempo 3.2x10 s en el dominio de estudio.
233
En la figura 6.3-III.3-Ψ.6 se presentan las gráficas con las isolíneas de la función de corriente
para un intervalo de tiempo 5.0 x 1011 s. y los distintos valores del factor de permeabilidad K
(K1 = 1x10-13 m2, K2 = 5x10-14 m2, K3 = 1x10-14 m2, K4 = 5x10-15 m2) en todo el dominio de estudio.
-13 2 -14 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
-14 2 -15 2
Factor de permeabilidad K=1x10 m Factor de permeabilidad K=5x10 m
11
Figura 6.3-III.3-𝚿.6 Isolíneas de la función de corriente para un intervalo de tiempo 5.0x10 s en el dominio de estudio.
234
CASO III.3-𝚿: Influencia del factor de permeabilidad K
235
6.3_III.4 Configuración de análisis 4: Influencia del factor de difusividad D
En la figura 6.3-III.4-𝑇.1 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 1.0 x 1010 s. y los distintos valores del factor de difusividad D
(D1 = 1x10-5 m2/s, D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
Se aprecia la evolución de las isolíneas hacia la superficie del dominio, en función de cuánto
mayor es el factor de difusividad D, más cerca de la superficie están las isolíneas.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
10
Figura 6.3-III.4-𝑻.1 Isolíneas de la temperatura adimensional para un intervalo de tiempo 1.0x10 s en el dominio de estudio.
236
En la figura 6.3-III.4-𝑇.2 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 3.0 x 1010 s. y los distintos valores del factor de difusividad D
(D1 = 1x10-5 m2/s, D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
10
Figura 6.3-III.4-𝑻.2 Isolíneas de la temperatura adimensional para un intervalo de tiempo 3.0x10 s en el dominio de estudio.
237
En la figura 6.3-III.4-𝑇.3 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 9.0 x 1010 s. y los distintos valores del factor de difusividad D
(D1 = 1x10-5 m2/s, D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
La evolución de las isolíneas sigue el mismo comportamiento del caso anterior, perdiéndose la
forma horizontal de las isolíneas.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
10
Figura 6.3-III.4-𝑻.3 Isolíneas de la temperatura adimensional para un intervalo de tiempo 9.0x10 s en el dominio de estudio.
238
En la figura 6.3-III.4-𝑇.4 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 1.2 x 1011 s. y los distintos valores del factor de difusividad D
(D1 = 1x10-5 m2/s, D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
La evolución de las isolíneas sigue el mismo comportamiento del caso anterior. Se aprecia el
inicio de la curvatura de las isolíneas.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
11
Figura 6.3-III.4-𝑻.4 Isolíneas de la temperatura adimensional para un intervalo de tiempo 1.2x10 s en el dominio de estudio.
239
En la figura 6.3-III.4-𝑇.5 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 3.2 x 1011 s. y los distintos valores del factor de difusividad D
(D1 = 1x10-5 m2/s, D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
La evolución de las isolíneas sigue el mismo comportamiento del caso anterior. Se aprecia, en el
último caso, la forma curvilínea de las isolíneas.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
11
Figura 6.3-III.4-𝑻.5 Isolíneas de la temperatura adimensional para un intervalo de tiempo 3.2x10 s en el dominio de estudio.
240
En la figura 6.3-III.4-𝑇.6 se presentan las gráficas con las isolíneas de temperatura adimensional
para un intervalo de tiempo 5.0 x 1011 s. y los distintos valores del factor de difusividad D
(D1 = 1x10-5 m2/s, D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
11
Figura 6.3-III.4-𝑻.6 Isolíneas de la temperatura adimensional para un intervalo de tiempo 5.0x10 s en el dominio de estudio.
241
CASO III.4-𝑻: Influencia del factor de difusividad D
242
En segundo lugar, vamos a estudiar la distribución de flujo de corriente Ψ en el dominio
geotérmico.
En la figura 6.3-III.4-Ψ.1 se presentan las gráficas con las isolíneas de la función de corriente
para un intervalo de tiempo 1.0 x 1010 s. y los distintos valores del factor de difusividad D
(D1 = 1x10-5 m2/s, D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
Se aprecian alteraciones de la forma de las isolíneas sobre la vertical del foco térmico y junto al
pozo. También como diferencia con el CASO I y II, en este CASO III encontramos isolíneas de valores
mayores hacia la parte superior derecha del dominio y a la derecha del pozo ya que la función de
corriente en esta zona es Ψ = 6.7x10-6.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
10
Figura 6.3-III.4-𝚿.1 Isolíneas de la función de corriente para un intervalo de tiempo 1.0x10 s en el dominio de estudio.
243
En la figura 6.3-III.4-Ψ.2 se presentan las gráficas con las isolíneas de la función de corriente
para un intervalo de tiempo 3.0 x 1010 s. y los distintos valores del factor de difusividad D
(D1 = 1x10-5 m2/s, D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
Para los dos mayores valores del factor de difusividad D, se aprecia la formación de isolíneas
ovoidales.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
10
Figura 6.3-III.4-𝚿.2 Isolíneas de la función de corriente para un intervalo de tiempo 3.0x10 s en el dominio de estudio.
244
En la figura 6.3-III.4-Ψ.3 se presentan las gráficas con las isolíneas de la función de corriente
para un intervalo de tiempo 9.0 x 1010 s. y los distintos valores del factor de difusividad D
(D1 = 1x10-5 m2/s, D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
10
Figura 6.3-III.4-𝚿.3 Isolíneas de la función de corriente para un intervalo de tiempo 9.0x10 s en el dominio de estudio.
245
En la figura 6.3-III.4-Ψ.4 se presentan las gráficas con las isolíneas de la función de corriente
para un intervalo de tiempo 1.2 x 1011 s. y los distintos valores del factor de difusividad D
(D1 = 1x10-5 m2/s, D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
Todos los casos siguen el mismo comportamiento del caso anterior. Ampliando la formación de
isolíneas ovoidales en el tercer caso.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
11
Figura 6.3-III.4-𝚿.4 Isolíneas de la función de corriente para un intervalo de tiempo 1.2x10 s en el dominio de estudio.
246
En la figura 6.3-III.4-Ψ.5 se presentan las gráficas con las isolíneas de la función de corriente
para un intervalo de tiempo 3.2 x 1011 s. y los distintos valores del factor de difusividad D
(D1 = 1x10-5 m2/s, D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
11
Figura 6.3-III.4-𝚿.5 Isolíneas de la función de corriente para un intervalo de tiempo 3.2x10 s en el dominio de estudio.
247
En la figura 6.3-III.4-Ψ.6 se presentan las gráficas con las isolíneas de la función de corriente
para un intervalo de tiempo 5.0 x 1011 s. y los distintos valores del factor de difusividad D
(D1 = 1x10-5 m2/s, D2 = 5x10-6 m2/s, D3 = 1x10-6 m2/s, D4 = 5x10-7 m2/s) en todo el dominio de estudio.
Todos los casos siguen el mismo comportamiento del caso anterior. Apreciando la formación de
isolíneas ovoidales en todos los casos.
-5 2 -6 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
-6 2 -7 2
Factor de difusividad D=1x10 m /s Factor de difusividad D=5x10 m /s
11
Figura 6.3-III.4-𝚿.6 Isolíneas de la función de corriente para un intervalo de tiempo 5.0x10 s en el dominio de estudio.
248
CASO III.4-𝚿: Influencia del factor de difusividad D
249
CAPÍTULO 7. CONCLUSIONES
En este Trabajo de Fin de Máster de título “Caracterización de escenarios geotérmicos con un
pozo de extracción, flujo regional y foco de calor en el fondo del dominio” se ha desarrollado el
estudio de la respuesta de modelos geotérmicos y las variaciones de temperatura y de flujo con el
tiempo y la profundidad en el dominio de estudio, esto es, la respuesta de la reserva geotermal para las
condiciones iniciales planteadas y diferentes intervalos de tiempo.
A continuación resumimos las principales conclusiones del trabajo para los CASOS I, II y III:
CASO I
Ψizq=4.9X10-6
Ψinf=3.9X10-6
Ψdcha=3.0X10-6
250
CASO II
Ψizq=4.9X10-6
Ψinf=3.0X10-6
Ψdcha=1.0X10-6
Para la función de corriente los CASOS I y II son iguales excepto en los valores de las isolíneas a la
derecha del foco térmico y del pozo, donde son inferiores en el CASO II.
251
CASO III
Ψizq=4.9X10-6
Ψinf=5.9X10-6
Ψdcha=6.7X10-6
Para la temperatura los CASOS I, II y III son muy similares en sus cuatro configuraciones.
Para la función de corriente en el CASO III encontramos isolíneas de valores mayores hacia la parte
superior derecha del dominio y a la derecha del pozo, ya que el valor de la función de corriente en
esta zona es más alta 𝚿 = 6.7x10-6, a diferencia de los CASOS I y II que están en la parte superior
izquierda del dominio y a la izquierda del pozo.
252
CAPÍTULO 8. BIBLIOGRAFIA
[1] Stuart J. Self, V. Reddy, Marc A. Rosen, (2013), Geothermal heat pump systems: Status review and
comparison with other heating options, Applied Energy 101 341-348, ELSEVIER.
[2] Tauqeer Abbas, Aqeel Ahmed Bazmi, Abdul Waheed Bhutto, Gholamreza Zahedi, (2014), Greener
energy: Issues and challenges for Pakistan- geothermal energy prospective, Renewable and Sustainable
Energy Reviews 31 258-269, ELSEVIER.
[3] Alhama I., Soto Meca A. y Alhama F. (2010 a). Programa FAHET (Flow and Heat Transport Simulator).
UPCT (Universidad Politécnica de Cartagena), Cartagena.
[4] Gómez- Lopera S.A., Mariñoso- Pascual J.M, Fernández- Bernal J.A., (2013-2104) “NANOTIO 2”
Copyright 2013-2014.
[5] Salvador A. Gómez- Lopera, Juan M. Mariñoso- Pascual, José A. Fernández- Bernal, (2102-2104),-
UPCT. Gómez - Lopera S.A., “DOPCALC”. Copyright 2012- 2014, Salvador A. Gómez – Lopera- UPCT.
[6] Gómez – Lopera S.A., Alhama F., del Cerro Velázquez F., Prodasim 1, (2005-2014) Copyright
2005-2014, UPCT.
[7] González-Fernández C.F. (2002). Applications of the network simulation method to transport
processes, in Network Simulation Method. Ed. J. Horno, Research, Signpost, Trivandrum, India.
[8] F. Del Cerro Velázquez, S.A. Gómez- Lopera, F. Alhama (2008). A powerful and versatile educational
software to simulate transient heat transfer processes in simple fins. Computer Applications in
Engineering Education, 16 (1), 72-82.
[9] Darcy H. (1856). Les fontaines publiques de la ville de Dijon. Libraire des Corps impériaux des ponts
et chaussées et des mines. Paris.
[10] Awkward K. W. y Lewis R. D. (1991). Models for determining thermal profiles of bond wires using
pspice analysis. Procc. 7° Annual Semiconductor Thermal> Measurement and Management Symposium,
AZ, USA.
[11] Bejan A. (1987). Convective heat transfer in a porous media. In: Kakac S., Shah R. K., y Aung W.
Eds. Handbook of single- phase convective heat transfer. Wiley, New York.
[12] Bénard H. (1900). Les tourbullons cellulaires dans nappe liqiude transportant de la chaleur par
convections en regime permanent. Rev. Gen. Sci. Pures Appl Bull Assoc 11, 1261-1271, 1309-1328
Chapman. A. J. (1984). Heat Transfer, 4ª ed. Macmillan Publishing Company. New York.
[13] Elder J. W. (1967). Transient convention in a porous médium, J. Fluid Mech., 27 (39) 69-96.
[14] Harrington D. A. y Van Der Driessche, P., (2004). Equivalent circuits for some surface
electrochemical mechanism. J. of Electroanlytical Chemistry, 5667, 153-166.
[15] Holzbecher E. y Yusa Y. (1995). Numerical experiments on free and forced convection in porous
media. Int. J. Heat and Mass Transfer 38, 2109-2115.
[16] Holzbecher E. (1998). Modelling Density-driven Flow in Porous Media. Springer, Berlín.
253
[17] Lapwood E. (1948). Convection of a fluid in a porous medium. Proc. Camb. Phil. Soc. 44, 508-521.
[18] Narayan K.A. y Amstrong D. (1995). Simulation of groundwater interception at Lake Ranfurly.
Victoria, incorporating variable density flow and solute transport. J. Hydrol. 165, 161-184.
[20] Patankar S. V. (1980). Numerical heat transfer and fluid flow. Hemisphere, New York.
[21] Peusner L. (1987). The principles of network thermodynamics: Theory and biophysical applications.
Entropy, Lincoln, Massachusetts.
[22] PSPICE 6.0 (1994), Microsim Corporation Fairbanks, Irvine, California 92718.
[23] Senger R.K. y Fogg G.E. (1990). Stream functions and equivalent freshwater heads for modelling
regional flow of variable-density groundwater, I, Review of theory and verification, Water Res. 26 (9),
2089-2096.
[24] Yusa Y. y Oishi I. (1989). Theoretical study of two-phase flow through porous medium (II), J. of the
Geothermal Research Society of Japan. 11 (3) 217-237 Tauqeer Abbas, Aqeel AhmedBazmi a, Abdul
WaheedBhutto c.
[25] Gholamreza Zahedi, Greener energy: Issues and challenges for Pakistan-geothermal energy
prospective.
[26] S.A. Gómez- Lopera, V. Martínez Martínez, (2014) Estudio numérico del problema geotérmico de
Elder mediante FAHET, influencia del Reyleight, UPCT (Universidad Politécnica de Cartagena), Cartagena,
Trabajo de Fin de Master.
[27] T. T. Chow, H. Long, H.Y. Mok, K. W. Li. Estimation of soil temperature profile in Hong Kong from
climatic variables.
[28] Stuart J. Self, V. Reddy, Marc A. Rosen, (2000) Publ., Geothermal heat pump systems: Status review
and comparison with other heating options, Faculty of Engineering and Applied Science, University of
Ontario Institute of Technology, Simcoe Street North, Oshawa, Ontario, Canada.
[29] Satoru Higuchi, Jun Nishijima and Yasuhiro Fujimitsu, Integration of Geothermal Exploration Data
Numerical Simulation Data Using Gis in a Hot Spring Area, Department of Earth Resources Engineering,
Graduate School of Engineering, Kyushu University, Fukuoka.
[30] Jiliang Chen, Fangming Jiang, Designing multi-well layout for enhanced geothermal system to
better exploit hot dry rock geothermal energy, Laboratory of Advanced Energy Systems, CAS Key
Laboratory of Renewable Energy, Guangzhou Institute of Energy Conversion, Chinese Academy of
Sciences.
[31] Benamar Bouhacina, Rachid Saim, Hakan F. Oztop, Numerical investigation of a novel tube design
for the geothermal borehole heat exchanger, Energetic and Applied Thermal Laboratory (ETAP), Faculty
of Technology, Abou Bakr Belkaid University, BP 230, 13000 Tlemcen, Algeria, Department of
Mechanical Engineering, Technology Faculty, Firat University, 23119 Elazig, Turkey.
254
[32] Knaber P, Frolkovic P., (1996), Consistent velocity approximation for finite volume or element
discretizations of density driven flow in porous media. In Aldama AA et al., editors. Computacional
methodsinwaterresources XI, vol. 1: Computational method sin sub-surface flow and transport
problems. Southampton: Computational Mechanics Publications; 1996. p. 93-100. (referencia 8).
[33] E. Holzbecher and Y. Yusa, (1983), Numerical experiments of free and forced convection in porous
media, Kyoto University, Beppu Geophysical Research Laboratory, Oita 874, Japan.
[34] Y. Yusa, Numerical experiment of groundwater motion under geothermal conditions, J. Geotherm.
Res. Sot. Jap. 5, 23-38.
[35] H.-J.G. Diersch, O. Kolditz Wasy, (Received 8 October 2001), Variable-density flow and transport in
porous media approaches and challenges Institute for Water Resources Planning and Systems Research,
Watersdorfer Str. 105, Berlin, Germany Center of Applied Geosciences, University of Tuubingen,
Germany.
[36] Ekkehard O. Holzbecher, (1998), Modeling Density-driven Flow in Porous Media, Springer Science &
Business Media, pag 16-cap 2.1.3.
[37] T. Graf, Simulation of Geothermal Flow in Deep sedimentary basins in Alberta, Energy resource
conservacion bord, Alberta geological Survey, Canada.
[38] Flores Fernández Giovanni Carlo, Sánchez Quispe, Sonia Tatiana, Correa González Alejandra y Ávila
Olivera Jorge Alejandro, (Octubre 2014), AMH, XXIII Congreso Nacional de hidráulica, Puerto Vallarta,
Jalisco, México, Modelos matemáticos de agua subterránea y de su importancia en México.
[39] M. Nabi, R. Al-Khoury, An efficient finite volume model for shallow geothermal systems. Part I:
Model formulation, Faculty of Civil Engineering and Geosciences, Delft University of Technology, The
Netherlands.
[40] Adrian E. Croucher, Michael J. O’Sullivan, Application of the computer code TOUGH2 to the
simulation of supercritical conditions in geothermal systems, Department of Engineering Science,
University of Auckland, Private Bag 92019, Auckland, New Zealand.
[42] Frank Kreit, (1973) Principi di trasmissione del calore, Liguori editore, de Principles of Heat Transfert
(Third Edicion), Intext Educational Publisher, New York.
[43] Marwan Fahs, Anis Younes, Thierry Alex Mar, A new benchmark semi-analvtical solution for
density-driven flow in porous media, a Laboratorie d’Hydrologie de Geochimie de Strasbourg, University
of Strasbourg, CNRS, UMR 7517, 1 rue Blessig, 67000 Strasbourg, France Department de Physique,
Piment, University of La Réunion, Moufia, La Réunion, France.
255
[44] H-J.G Diersch, O. Kolditz, Wasy, (2001), Variable-density flow en tranport in pourous media:
approches and challenges Institute for Water Resources Planning and Systems Research, Watersdorfer
Str. 105, Berlin, Germany, Center of Applied Geosciences, University of Tubingen, Germany,.
[45] M. Filippucci, A Tallarico, (2012) La modellizzazione físico matemática per l’ottimizzazione dello
sfruttamento delle risorse geotermiche, Università degli Studi di Bari, incontro dibattito, energía
sostenibile per tutti.
[46] Gómez- Lopera S.A., Mariñoso- Pascual J.M., Fernéndez- Bernal J.A., (2013-2014), “NANOTIO2”.
“2013-2014, Salvador A. Gómez- Lopera, Juan M. Mariñoso- Pascual, J. A. Fernández- Bernal – UPCT
[47] Alhama I., Soto Meca A. y Alhama F. (2010). Programa FAHET (Flow and Heat Transport Simulator).
“UPCT (Universidad Politécnica de Cartagena), Cartagena Parameter identification and uncertainty
analysis for heat transfer at the KTB drill site using a 2-D inverse method Holger Lehmann, Kelin Wang,
Christoph Clauser.
[48] Christian Vogt, Norbert Klitzscha, Volker Rath On self-potential data for estimating permeability in
enhanced geothermal systems,.
[49] O. Kolditz, C. Clauser, Numerical simulation of flow and heat transfer in fractured.
[50] H. Lehmann, K. Wang, C. Clauser, Parameter identification and uncertainty analysis for heat
transfer at the KTB drill site using a 2-D inverse method.
[51] W. R. Van Wijk, D.A De Vries, (1963), Periodic temperature variation in a homogeneous soil, Physics
of plant environment, Amsterdam, North Holland.
[52] Joaquin Cruces de Abia, (2006-2007), Contaminación ambiental, modelo MODFLOW, Master en
ingeniería y gestión del agua.
[53] J.L. Cenis, (1989) Temperature evaluación in solarized soil by Fourier analisis. Phytopatology 79:
509-510.
[54] Assesment of geothermal Resource of United Stated (1975), Geological Survey Circular 726, pag.
84.
256