Sensitivity of groundwater flow with respect to the drain–aquifer leakage coefficient

Mitigation measures may be used to prevent soil and water pollution from waste disposal, landfill sites, septic or chemical storage tanks. Among them, drains and impervious barriers may be set up. The efficiency of this technique can be evaluated by means of groundwater modeling tools. The groundwater flow and the leakage drain–aquifer interactions are implemented in a conforming finite element method (FEM) and a mixed hybrid FEM (MHFEM) in a horizontal two-dimensional domain modeling regional aquifer below chemical storage tanks. Considering the influence of uncertainties in the drain–aquifer exchange rate parameter and using an automatic differentiation (AD) tool, the aim of this paper is to carry out a sensitivity analysis with respect to the leakage coefficient for the piezometric head, velocity field, and streamlines to provide a new insight into groundwater waterbody exchanges. Computations are performed with both an ideal homogeneous hydraulic conductivity and a realistic heterogeneous one. The tangent linear codes are validated using Taylor tests performed on the head and the velocity field. The streamlines computed using AD are well approximated in comparison with the nondifferentiated codes. Piezometric head computed by the MHFEM is the more sensitive, particularly near to the drain, than the FEM one.

[1]  Ilio Galligani,et al.  Mathematical Aspects of Finite Element Methods , 1977 .

[2]  Vincent Guinot,et al.  Sensitivity analysis of two-dimensional steady-state aquifer flow equations. Implications for groundwater flow model calibration and validation , 2010 .

[3]  Andreas Griewank,et al.  ADIFOR - Generating Derivative Codes form Fortran Programs , 1992, Sci. Program..

[4]  Thierry A. Mara,et al.  Bayesian sparse polynomial chaos expansion for global sensitivity analysis , 2017 .

[5]  H. Franssen,et al.  Field evidence of a dynamic leakage coefficient for modelling river–aquifer interactions , 2007 .

[6]  K. R. Rushton,et al.  Possible mechanisms for leakage between aquifers and rivers , 1979 .

[7]  G. Chavent,et al.  Application of the mixed hybrid finite element approximation in a groundwater flow model: Luxury or necessity? , 1994 .

[8]  J. Bear.  Hydraulics of Groundwater , 1979 .

[9]  Laurent Hascoët,et al.  The Tapenade automatic differentiation tool: Principles, model, and specification , 2013, TOMS.

[10]  J. Morgenstern,et al.  How to compute fast a function and all its derivatives: a variation on the theorem of Baur-strassen , 1985, SIGA.

[11]  Andreas Griewank,et al.  Evaluating derivatives - principles and techniques of algorithmic differentiation, Second Edition , 2000, Frontiers in applied mathematics.

[12]  Jean Utke,et al.  Fast higher-order derivative tensors with Rapsodia , 2009, Optim. Methods Softw..

[13]  Michael Bruen,et al.  Sensitivity of stream–aquifer seepage to spatial variability of the saturated hydraulic conductivity of the aquifer , 2004 .

[14]  Graham E Fogg,et al.  River‐Aquifer Interactions, Geologic Heterogeneity, and Low‐Flow Management , 2006, Ground water.

[15]  David Elizondo,et al.  Automatic versus manual model differentiation to compute sensitivities and solve non-linear inverse problems , 2002 .

[16]  Jean E. Roberts,et al.  A unified physical presentation of mixed, mixed-hybrid finite elements and standard finite difference approximations for the determination of velocities in waterflow problems , 1991 .

[17]  Udo Meissner,et al.  A mixed finite element model for use in potential flow problems , 1973 .

[18]  Marwan Fahs,et al.  An easy and efficient combination of the Mixed Finite Element Method and the Method of Lines for the resolution of Richards' Equation , 2009, Environ. Model. Softw..

[19]  W. Kinzelbach,et al.  Continuous Groundwater Velocity Fields and Path Lines in Linear, Bilinear, and Trilinear Finite Elements , 1992 .

[20]  Uwe Naumann,et al.  The Art of Differentiating Computer Programs - An Introduction to Algorithmic Differentiation , 2012, Software, environments, tools.

[21]  P. Ackerer,et al.  Mixed finite elements for solving 2‐D diffusion‐type equations , 2010 .

[22]  Jon C. Helton,et al.  Sensitivity Analysis of Model Output: SAMO 2004 , 2006, Reliab. Eng. Syst. Saf..

[23]  I. Charpentier,et al.  Local Sensitivity Analysis of a Numerical Model of Volcanic Plinian Columns through Automatic Differentiation , 2005 .

[24]  Isabelle Charpentier.  The MesODiF package for gradient computations with the atmospheric model Meso-NH , 2000, Environ. Model. Softw..

[25]  Andreas Griewank,et al.  Algorithm 755: ADOL-C: a package for the automatic differentiation of algorithms written in C/C++ , 1996, TOMS.