Modelación geoestadística de medios geológicos naturalmente fracturados

Abstract

La modelación del flujo en medios geológicos naturalmente fracturados depende principalmente de la composición estructural de la relación matriz-fractura y de la forma en que se aproxime el intercambio de fluido entre ambas fases. Las dos características mencionadas (composición estructural de la doble porosidad y la función de transferencia de flujo) incorporan el mayor grado de incertidumbre en la modelación de los estados transitorios generados por bombeo. Para poder mejorar la interpretación del ensayo de bombeo y con el objetivo último de obtener parámetros característicos confiables del acuífero (tales como la conductividad y el coeficiente de almacenamiento para ambas fases), la metodología propuesta en esta tesis aprovecha las herramientas que brinda la geoestadística para simular realizaciones de campos de fisuras cuyo ensamble permite a su vez estimar los parámetros de forma de la función de transferencia matriz-fisura (FTMF). En efecto, los métodos clásicos disponibles para interpretar resultados de ensayos de bombeo en acuíferos heterogéneos, en especial aquellos que se presentan naturalmente fracturados y con una compleja relación entre las dos componentes porosas, no contienen una vinculación directa entre la estructura interna de un acuífero con doble porosidad y los parámetros hidrogeológicos que definen el flujo a través de ella. Esto sólo se ha logrado con geometrías muy simples de la estructura matriz-fisura o con realizaciones estocásticas de la red de fisuras que luego requieren un gran esfuerzo computacional al simularlos con modelos hidrodinámicos. Lo que se presenta en este trabajo de tesis es una metodología sencilla y práctica para lograr mejorar la estimación de los parámetros de un ensayo de bombeo en este tipo de acuíferos confinados naturalmente fracturados, incorporando por primera vez la estructura y la proporción de las familias de fisuras en las ecuaciones de pozo a través de la estimación directa del factor de forma de la función de transferencia, incidiendo principalmente en las etapas tempranas del régimen transitorio del bombeo. Para verificar su funcionamiento, la metodología desarrollada fue puesta a prueba en algunos casos sintéticos orientados a la interpretación de ensayos de bombeo en régimen variable. Además, se aplicó dicha metodología a la interpretación de las características hidrodinámicas de un acuífero cretácico confinado del Distrito uranífero Cerro Solo, en la provincia de Chubut, Argentina. A partir de las aplicaciones mencionadas, se discuten los resultados obtenidos por comparación con los ejemplos tomados de la literatura presentando el contraste con soluciones analíticas desarrolladas para casos particulares de acuíferos con doble porosidad y con la re-interpretación de la evolución temporal de las depresiones observadas en el acuífero ensayado. En las conclusiones se destaca la importancia de la consideración de la doble porosidad en las ecuaciones de flujo, alentando a su aplicación concreta en la práctica de interpretación de los ensayos de bombeo con la metodología presentada en esta tesis y a proseguir su evolución en etapas de desarrollo vinculadas al transporte de solutos en este tipo de acuíferos complejos.