Semi-analytic modeling of transient multi-layer flow with TTim

TTim is a free code for the semi-analytic simulation of transient flow in multi-layer systems consisting of an arbitrary number of layers. No grid or time-stepping is required, nor does a closed model boundary need to be specified in any of the layers. Currently, TTim includes multi-layer wells and line-sinks, which may be used to simulate transient flow to a variety of hydrogeologic features, including wells with a skin and wellbore storage, incompletely sealed abandoned wells, streams with leaky beds, vertical faults, and horizontal wells; transient forcing needs to be represented by a step function. Other features that may be simulated include vertical anisotropy and the delayed response of the water table. Behind the scenes of TTim, the Laplace-transform analytic element method is applied. TTim is written in Python, with Python scripts used as input files. TTim has many practical applications, including the design of riverbank filtration systems, analysis of aquifer tests near surface-water bodies, design and evaluation of recirculation wells, and modeling of the transient pressure response of proposed carbon geologic sequestration projects. In addition, the short and simple input files and the one-to-one link between analytic elements and hydrogeologic features make TTim well suited for education.RésuméTTim est un code libre de droit pour la simultation semi-analytique d’écoulement transitoire de systèmes multi-couches composés d’un nombre arbitraire de couches. Aucune grille ou pas de temps n’est requise, pas plus que la limite fermée du modèle doit être spécifiée pour l’ensemble des couches. Actuellement, TTim comprend des puits multi-couche et des points d’engouffrement, qui peuvent être utilisés pour simuler l’écoulement transitoire pour une variété de configurations hydrogéologiques, telles que des puits avec un effet de peau et un effet capacitif au puits, des puits abandonnés à étanchéité incomplète, des cours d’eau à lit infiltrant, des failles verticales, et des puits horizontaux; le forçage transitoire nécessite d’être représenté par une fonction en échelon. D’autres configurations pouvant être simulées comprennent l’anisotropie verticale et la réponse retardée du niveau piézométrique. Derrière le montage de TTim, la méthode analytique élémentaire de la transformée de Laplace est appliquée. TTim est écrit en langage Python, avec des scripts utilisés comme fichiers d’entrée. TTim possède plusieurs applications pratiques, telles que la conception de systyèmes d’infiltration au niveau des rives de cours d’eau, l’analyse des essais d’aquifère à proximité de cours d’eau, la conception et l’évaluation de puits permettant la recirculation de l’eau et la modélisation des réponses en régime transitoire de la pression dans le cadre de projets de séquestration géologique du CO2. De plus, les fichiers d’entrée sont courts et simples et la relation entre les éléments analytiques et les configurations hydrogéologiques rendent TTim approprié pour l’enseignement.ResumenTTim es un programa libre para la simulación semianalítica de flujo transitorio en sistemas multicapas consistentes en un número arbitrario de capas. No se requiere una grilla ni un paso temporal, tampoco es necesario especificar un modelo de bordes cerrados para ninguna de las capas. Actualmente, TTim incluye pozos multicapas y sumideros en líneas, que puede ser usado para simular el flujo transitorio para una variedad de aspectos hidrogeológicos, incluyendo pozos con un almacenamiento pelicular y en el pozo, pozos abandonados incompletamente sellados, corrientes con capas filtrantes, fallas verticales y pozos horizontales; las necesidades transitorias obligan a estar representadas por una función escalón. Otras características que pueden ser simuladas incluyendo la anisotropía vertical y la respuesta retardada de los niveles freáticos. Detrás de las presentaciones de TTim, se aplicó un método de elementos analíticos de la transformada de Laplace. TTim está escrito en Python, con escritura de Python usada como archivos de entrada. TTim tiene muchas aplicaciones prácticas, inclusive el diseño de sistema de filtración de las márgenes del río, análisis de ensayos de acuíferos próximos a los cuerpos de agua superficial, diseño y evaluación de pozos de recirculación, y modelado de respuesta transitoria de la presión en proyectos de secuestro geológico de carbón. Además, los cortos y simples archivos de entradas y el enlace uno a uno entre los elementos analíticos y los aspectos hidrogeológicos hacen que el TTim sea muy apropiado para la educación.ResumoO TTim é um software livre para simulação semi-analítica do escoamento transitório em sistemas multi-camada constituídos por um número arbitrário de camadas. O software não requer nenhuma discretização espacial ou passos de tempo, nem é necessário especificar um limite de modelo fechado para qualquer das camadas. Presentemente, o TTim inclui poços e sumidouros lineares multi-camada que podem ser usados para simular o escoamento transitório para uma variedade de entidades hidrogeológicas, incluindo poços com efeito de pele e com armazenamento, poços abandonados com selagem incompleta, cursos de água com leito drenante, falhas verticais e poços horizontais; o caráter transitório tem que ser representado por uma função em degraus. Outras funcionalidades que podem ser simuladas incluem a anisotropia vertical e a resposta diferida do nível freático. Por trás do TTrim é aplicado o método de elementos analíticos da transformada de Laplace. O TTim é escrito em linguagem Python, sendo usados ficheiros Python como entrada de dados. O TTim tem muitas aplicações práticas, incluindo o projeto de sistemas de filtração em margem de cursos de água, a análise de ensaios em aquíferos próximos de corpos de água superficial, o projeto e a avaliação de poços de recirculação e a modelação da pressão transitória resultante dos propostos projetos de sequestro geológico de carbono. Para além disso, os ficheiros de entrada, curtos e simples, e a conexão direta entre os elementos analíticos e as caraterísticas hidrogeológicas tornam o TTim bem adaptado ao ensino.

[1]  M. Bakker Simulating groundwater flow to surface water features with leaky beds using analytic elements , 2007 .

[2]  Mark Bakker,et al.  Multilayer analytic element modeling of radial collector wells. , 2005, Ground water.

[3]  H. M. Haitjema Analytic Element Modeling , 1995 .

[4]  Henk M. Haitjema,et al.  Analytic Element Modeling of Groundwater Flow , 1995 .

[5]  A. N. Stokes,et al.  An Improved Method for Numerical Inversion of Laplace Transforms , 1982 .

[6]  Fernando Pérez,et al.  Python: An Ecosystem for Scientific Computing , 2011, Computing in Science & Engineering.

[7]  S. P. Neuman,et al.  Laplace-transform analytic element solution of transient flow in porous media , 2003 .

[8]  M. Bakker,et al.  Groundwater flow through anisotropic fault zones in multiaquifer systems , 2008 .

[9]  Hongbin Zhan,et al.  Groundwater flow to a horizontal or slanted well in an unconfined aquifer , 2002 .

[10]  G. Kruseman,et al.  Analysis and Evaluation of Pumping Test Data , 1983 .

[11]  F. Carlson,et al.  MLU: A Windows Application for the Analysis of Aquifer Tests and the Design of Well Fields in Layered Systems , 2012 .

[12]  Kristopher L. Kuhlman,et al.  Laplace-transform analytic-element method for transient porous-media flow , 2008 .

[13]  Mark Bakker,et al.  Simulation of axi-symmetric flow towards wells: A finite-difference approach , 2012, Comput. Geosci..

[14]  Travis E. Oliphant,et al.  Python for Scientific Computing , 2007, Computing in Science & Engineering.

[15]  Kristopher L. Kuhlman,et al.  Computational issues and applications of line-elements to model subsurface flow governed by the modified Helmholtz equation , 2011 .

[16]  C. J. Hemker,et al.  Unsteady flow to wells in layered and fissured aquifer systems , 1987 .

[17]  Jens Birkholzer,et al.  Analytical solutions for pressure perturbation and fluid leakage through aquitards and wells in multilayered‐aquifer systems , 2011 .

[18]  Otto D Strack,et al.  Three‐Dimensional Streamlines in Dupuit‐Forchheimer Models , 1984 .

[19]  Otto D Strack,et al.  Groundwater mechanics , 1989 .

[20]  Robert J. Moorhead,et al.  Robust Loop Detection for Interactively Placing Evenly Placed Streamlines , 2007, Computing in Science & Engineering.

[21]  Charles McLane AnAqSim: Analytic Element Modeling Software for Multi‐Aquifer, Transient Flow , 2012 .

[22]  A. W. Harbaugh MODFLOW-2005 : the U.S. Geological Survey modular ground-water model--the ground-water flow process , 2005 .

[23]  M. Bakker Analytic Modeling of Transient Multilayer Flow , 2013 .

[24]  C. R. Fitts Modeling aquifer systems with analytic elements and subdomains , 2010 .

[25]  Otto D Strack,et al.  Theory and applications of the Analytic Element Method , 2003 .

[26]  Randall J Hunt,et al.  Ground Water Modeling Applications Using the Analytic Element Method , 2006, Ground water.

[27]  Arnold Verruijt,et al.  Theory of groundwater flow , 1970 .

[28]  John D. Hunter,et al.  Matplotlib: A 2D Graphics Environment , 2007, Computing in Science & Engineering.

[29]  S. P. Neuman Theory of flow in unconfined aquifers considering delayed response of the water table , 1972 .

[30]  Randal J. Barnes,et al.  High-order line elements in modeling two-dimensional groundwater flow , 1999 .

[31]  Mark Bakker,et al.  Writing Analytic Element Programs in Python , 2009, Ground water.