Study of groundwater flow from temperature data using inverse methods. Applications to hydrogeological scenarios around mar menor lagoon

Resumen

Resumen de la tesis: Los objetivos principales de este trabajo son la caracterización adimensional de escenarios, la obtención de soluciones universales, así como su representación gráfica, y la propuesta de protocolos de problema inverso para la estimación tanto del caudal de agua subterránea como de las propiedades térmicas y geométricas del medio poroso. La aplicación de los resultados obtenidos al escenario del Mar Menor quedó como un objetivo añadido debido a la gran cantidad de escenarios caracterizados. Los distintos casos estudiados han sido: i) escenario de flujo vertical unidimensional (ascendente y descendente), ii) flujo horizontal en acuíferos saturados, iii) flujo horizontal a través de una capa permeable situada entre dos regiones impermeables, iv) flujo oblicuo de velocidad constante, v) escenarios en los cuales los potenciales hidráulicos son constantes en las fronteras del dominio y, vi) interacción río acuífero. La técnica utilizada para la obtención de soluciones universales ha sido la adimensionalización discriminada y normalizada del modelo matemático del problema. El método de simulación por redes, como herramienta para el diseño de los modelos en red, posibilita el estudio de cualquier proceso que pueda definirse mediante un conjunto de ecuaciones que formen un modelo matemático. Para el escenario de flujo vertical 1D, se ha introducido una longitud característica (magnitud oculta en el enunciado del problema) para la cual los efectos de conducción y advección son comparables. Para la condición de temperatura constante en la superficie se ha demostrado que las soluciones son aplicables tanto para el caso de temperatura armónica en superficie (siempre que se refieran a la temperatura media del ciclo) como para el caso de temperatura constante en la superficie. Para la condición de variación sinusoidal de la temperatura en la superficie del terreno se aporta, por primera vez en la literatura científica, la distinción entre acuíferos profundos y acuíferos superficiales desde el punto de vista térmico. Para ello, es necesario introducir una longitud (o profundidad) característica vertical que define la extensión del acuífero donde la onda armónica de temperatura es aún apreciable. La aplicación más directa de los resultados obtenidos es su uso para la estimación de flujos (velocidades) y/o propiedades térmicas y geométricas del acuífero en la forma de problema inverso. En el segundo escenario estudiado (flujo horizontal en acuíferos saturados), la aparición de una longitud característica horizontal, permite caracterizar los perfiles temperatura-profundidad dentro de la región delimitada por dicha longitud. Una de las contribuciones más relevantes es la referente a la demostración de que los gradientes térmicos horizontales en ningún caso son constantes, como se asume por hipótesis en los trabajos científicos para estos escenarios. Con base a los resultados obtenidos, se propone un protocolo de problema inverso clásico mediante simulaciones numéricas. Respecto al problema de flujo de agua horizontal a través de una capa permeable, la introducción, por un lado, de una longitud característica horizontal, a lo largo de la cual se desarrollan los perfiles verticales de temperatura, y por otro, de un tiempo característico que marca la duración del período transitorio, permiten la estimación, de forma directa, de las dependencias de las formas adimensionales de estas dos incógnitas sobre el resto de grupos del problema. Se establecen dos protocolos de problema inverso para la estimación de la velocidad horizontal del flujo de agua subterránea. Para el escenario de flujo oblicuo de velocidad constante, existe una región característica en la cual los perfiles de temperatura dependen de las componentes vertical y horizontal de la velocidad del flujo. Más allá de esta región, los perfiles verticales de temperatura únicamente dependen de la componente vertical y son independientes de la posición horizontal. Para el escenario de flujo oblicuo de potenciales hidráulicos de valor constante, la velocidad de flujo no es constante en cada punto del medio poroso, sino que se va desarrollando gradualmente dentro de una región más allá de la cual, el flujo horizontal se hace nulo. En escenarios de flujo constante se presenta un protocolo de problema inverso para la estimación de las componentes vertical y horizontal, respectivamente. Finalmente, para el escenario de interacción río-acuífero, el flujo no es constante a lo largo de una región horizontal a lo largo de la cual la componente vertical irá perdiendo importancia frente a la horizontal. Más allá de esta región definida a partir de una longitud característica hidráulica, la velocidad vertical es despreciable. En este escenario, se puede desacoplar el problema hidráulico del problema térmico. Por lo tanto, el caudal adimensional, los potenciales hidráulicos adimensionales y el valor de la longitud característica horizontal hidráulica, sólo dependen de la geometría del escenario, de las componentes de la permeabilidad y del valor de los potenciales hidráulicos en la frontera. Sin embargo, los patrones de temperatura, además de depender de las condiciones térmicas y mecánicas en los contornos, dependerán de las propiedades térmicas e hidráulicas del medio, de la geometría y de la posición. Por último, el protocolo de problema inverso clásico propuesto para el caso de flujo horizontal ha sido aplicado en una zona del entorno del Mar Menor situada en Los Alcázares, estimándose la velocidad horizontal y la temperatura de entrada del agua subterránea. http://repositorio.bib.upct.es/dspace/