An inverse problem methodology to identify flow channels in fractured media using synthetic steady-state head and geometrical data

Abstract We present a methodology for identifying highly-localized flow channels embedded in a significantly less permeable medium using steady-state head and geometrical data. This situation is typical of fractured media where flows are often strongly channeled at the scales of interest (10 m–1 km). The objective is to identify both geometrical and hydraulic characteristics of the conducting structures. Channels are identified in decreasing order of importance by successive optimizations of an objective function. The identification strategy takes advantage of the hierarchical flow organization to restrict the dimension of the solution space of each individual optimization step. The characteristics of the secondary channels are strongly determined by the main flow channels. The latter are slightly modified by the secondary channels through the addition of a regularization term to the main channel characteristics in the objective function. As the objective function is strongly non-convex with numerous local minima, inversion is performed using a stochastic algorithm (simulated annealing). We assess the possibilities of the hierarchical identification strategy on simple synthetic steady-state flow configurations where hydraulic data are made up of 25 regularly spaced heads and of the boundary conditions. Those flow structures that are dominated by at most two simple channels can be identified with these head data only. Configurations comprising up to three complex and interconnected channels can still be identified with additional geometrical information including the distances of piezometers to their closest channel. The capabilities of the hierarchical identification strategy are limited to flow structures dominated by at most three equivalent flow channels. We finally discuss the perspectives of application of the method to transient-state data obtained on a more restricted number of piezometers.

[1]  Hakim Benabderrahmane,et al.  Evaluation of retention properties of a semi-synthetic fractured block from modelling at performance assessment time scales (Äspö Hard Rock Laboratory, Sweden) , 2009 .

[2]  A. Tarantola.  Inverse problem theory : methods for data fitting and model parameter estimation , 1987 .

[3]  Chin-Fu Tsang,et al.  Flow channeling in a single fracture as a two‐dimensional strongly heterogeneous permeable medium , 1989 .

[4]  David E. Dougherty,et al.  Optimal groundwater management: 2. Application of simulated annealing to a field-scale contamination site , 1993 .

[5]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.

[6]  S. P. Neuman,et al.  Three‐dimensional numerical inversion of pneumatic cross‐hole tests in unsaturated fractured tuff: 2. Equivalent parameters, high‐resolution stochastic imaging and scale effects , 2001 .

[7]  Antoine Mensch,et al.  An inverse technique for developing models for fluid flow in fracture systems using simulated annealing , 1993 .

[8]  Klaus Mosegaard,et al.  MONTE CARLO METHODS IN GEOPHYSICAL INVERSE PROBLEMS , 2002 .

[9]  A Blocking Markov Chain Monte Carlo Method for Inverse Stochastic Hydrogeological Modeling , 2009 .

[10]  Mark Grasmueck,et al.  3-D ground‐penetrating radar applied to fracture imaging in gneiss , 1996 .

[11]  S. P. Neuman,et al.  Estimation of Aquifer Parameters Under Transient and Steady State Conditions: 1. Maximum Likelihood Method Incorporating Prior Information , 1986 .

[12]  D. McLaughlin,et al.  A Reassessment of the Groundwater Inverse Problem , 1996 .

[13]  H. Yeh,et al.  Determination of the parameter pattern and values for a one-dimensional multi-zone unconfined aquifer , 2008 .

[14]  Stephen E. Silliman,et al.  An interpretation of the difference between aperture estimates derived from hydraulic and tracer tests in a single fracture , 1989 .

[15]  Nikolaus Hansen,et al.  Adapting arbitrary normal mutation distributions in evolution strategies: the covariance matrix adaptation , 1996, Proceedings of IEEE International Conference on Evolutionary Computation.

[16]  W. Yeh,et al.  Optimal observation network design for parameter structure identification in groundwater modeling , 2005 .

[17]  G. de Marsily,et al.  Four decades of inverse problems in hydrogeology , 2000 .

[18]  Frederick D. Day-Lewis,et al.  Identifying fracture‐zone geometry using simulated annealing and hydraulic‐connection data , 2000 .

[19]  K. Karasaki,et al.  Hydraulic well testing inversion for modeling fluid flow in fractured rocks using simulated annealing: a case study at Raymond field site, California , 2000 .

[20]  Jean-Raynald de Dreuzy,et al.  Hydraulic properties of two‐dimensional random fracture networks following a power law length distribution: 1. Effective connectivity , 2001 .

[21]  Graham E. Fogg,et al.  Groundwater Flow and Sand Body Interconnectedness in a Thick, Multiple-Aquifer System , 1986 .

[22]  C. McDermott,et al.  New experimental techniques for pneumatic tomographical determination of the flow and transport parameters of highly fractured porous rock samples , 2003 .

[23]  Mrinal K. Sen,et al.  Vertical fracture detection by exploiting the polarization properties of ground-penetrating radar signals , 2004 .

[24]  Jin-Ping Gwo,et al.  In search of preferential flow paths in structured porous media using a simple genetic algorithm , 2001 .

[25]  Philippe Davy,et al.  Statistical characteristics of flow as indicators of channeling in heterogeneous porous and fractured media , 2010 .

[26]  F. Mazzarini,et al.  A structural and geophysical approach to the study of fractured aquifers in the Scansano-Magliano in Toscana Ridge, southern Tuscany, Italy , 2009 .

[27]  F. Tsai,et al.  A Combinatorial Optimization Scheme for Parameter Structure Identification in Ground Water Modeling , 2003, Ground water.

[28]  J. Bruyelle,et al.  Automated characterization of fracture conductivities from well tests inversion , 2009 .

[29]  Jean-Raynald de Dreuzy,et al.  Hydraulic properties of two‐dimensional random fracture networks following a power law length distribution: 2. Permeability of networks based on lognormal distribution of apertures , 2001 .

[30]  P. Patrick Wang,et al.  Parameter structure identification using tabu search and simulated annealing , 1996 .

[31]  Harald Stehfest,et al.  Remark on algorithm 368: Numerical inversion of Laplace transforms , 1970, CACM.

[32]  S. P. Neuman,et al.  Theory, modeling, and field investigation in hydrogeology : a special volume in honor of Shlomo P. Neuman's 60th birthday , 2000 .

[33]  S. Benson,et al.  Hydrologic characterization of heterogeneous geologic media with an inverse method based on iterated function systems , 1994 .

[34]  Ghislain de Marsily,et al.  Three‐dimensional interference test interpretation in a fractured aquifer using the Pilot Point Inverse Method , 2001 .

[35]  D. A. Zimmerman,et al.  A comparison of seven geostatistically based inverse approaches to estimate transmissivities for modeling advective transport by groundwater flow , 1998 .

[36]  G. Marsily,et al.  Modeling fracture flow with a stochastic discrete fracture network: calibration and validation: 1. The flow model , 1990 .

[37]  C. Tsang,et al.  Flow channeling in heterogeneous fractured rocks , 1998 .

[38]  L. Ingber.  Adaptive Simulated Annealing (ASA) , 1993 .

[39]  Mary C Hill,et al.  The Practical Use of Simplicity in Developing Ground Water Models , 2006, Ground water.

[40]  Timothy A. Davis,et al.  Algorithm 832: UMFPACK V4.3---an unsymmetric-pattern multifrontal method , 2004, TOMS.

[41]  Peter K. Kitanidis,et al.  An interactive Bayesian geostatistical inverse protocol for hydraulic tomography , 2006 .

[42]  Peter Dietrich,et al.  A travel time based hydraulic tomographic approach , 2003 .

[43]  J. Hanor.  Effective hydraulic conductivity of fractured clay beds at a hazardous waste landfill, Louisiana Gulf Coast , 1993 .

[44]  G Bruckner,et al.  An inverse problem from 2D ground-water modelling , 1998 .

[45]  Dominique Bruel,et al.  Using the migration of the induced seismicity as a constraint for fractured Hot Dry Rock reservoir modelling , 2007 .

[46]  Andres Alcolea,et al.  Inverse problem in hydrogeology , 2005 .

[47]  Emmanuel Ledoux,et al.  Modeling fracture flow with a stochastic discrete fracture network: Calibration and validation: 2. The transport model , 1990 .

[48]  Arnaud Lange,et al.  An extended evolution strategy for the characterization of fracture conductivities from well tests , 2009, GECCO.

[49]  Yonghong Hao,et al.  Hydraulic Tomography for Detecting Fracture Zone Connectivity , 2008, Ground water.

[50]  William W.-G. Yeh,et al.  Aquifer parameter identification with optimum dimension in parameterization , 1981 .

[51]  Jesús Carrera,et al.  An analysis of hydraulic conductivity scale effects in granite (Full‐scale Engineered Barrier Experiment (FEBEX), Grimsel, Switzerland) , 2005 .

[52]  P. Witherspoon,et al.  Porous media equivalents for networks of discontinuous fractures , 1982 .

[53]  Lester Ingber,et al.  Simulated annealing: Practice versus theory , 1993 .

[54]  C. Renshaw.  Estimation of fracture zone geometry from steady-state hydraulic head data using iterative sequential cokriging , 1996 .