CONDITIONING DISCRETE FRACTURE NETWORK MODELS OF GROUNDWATER FLOW

Many geological formations consist of crystalline rock that have very low matrix permeability but allow flow through an interconnected net- work of fractures. Understanding the flow of groundwater through such rocks is important in considering disposal of radioactive waste in underground reposi- tories. A specific area of interest is the conditioning of fracture transmissivities on measured values of pressure in these formations. While there are exist- ing methods to condition transmissivity fields on transmissivity, pressure and flow measurements for a continuous porous medium, considerably less work has been devoted to conditioning discrete fracture networks. This article presents two new methods for conditioning fracture transmissivities on measured pres- sures in a discrete fracture network. The first approach adopts a linear ap- proximation when fracture transmissivities are mildly heterogeneous, while the minimisation of a suitable objective function is undertaken when fracture trans- missivities are highly heterogeneous. The second conditioning algorithm is a Bayesian method that finds a maximum a posteriori (MAP) estimator which maximises the posterior distribution defined by Bayes' theorem using informa- tion from the prior distribution of fracture transmissivities and observations in the form of measured pressures. The conditioning methods are tested on two separate, large scale test cases that model a potential site for radioactive waste disposal. Results from these test cases are shown and comparisons between the two conditioning methods are made.

[1]  Andrew M. Stuart,et al.  Inverse problems: A Bayesian perspective , 2010, Acta Numerica.

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

[3]  S. P. Neuman,et al.  A statistical approach to the inverse problem of aquifer hydrology: 2. Case study , 1980 .

[4]  S. P. Neuman,et al.  Estimation of aquifer parameters under transient and steady-state conditions: 2 , 1986 .

[5]  S. P. Neuman,et al.  A statistical approach to the inverse problem of aquifer hydrology: 1. Theory , 1979 .

[6]  J. Bear Dynamics of Fluids in Porous Media , 1975 .

[7]  D. T. Snow,et al.  The frequency and apertures of fractures in rock , 1970 .

[8]  Claire R. Tiedeman,et al.  Characterizing a Ground Water Basin in a New England Mountain and Valley Terrain , 1998 .

[9]  Conditioning stochastic groundwater flow models on head data , 1994 .

[10]  William H. Press,et al.  Numerical recipes , 1990 .

[11]  S. P. Neuman Calibration of distributed parameter groundwater flow models viewed as a multiple‐objective decision process under uncertainty , 1973 .

[12]  P. Kitanidis On the geostatistical approach to the inverse problem , 1996 .

[13]  C. Tiedeman,et al.  Effective Groundwater Model Calibration , 2007 .

[14]  Lanru Jing,et al.  A review of techniques, advances and outstanding issues in numerical modelling for rock mechanics and rock engineering , 2003 .

[15]  Gudmundur S. Bodvarsson,et al.  Hydraulic conductivity of rock fractures , 1996 .

[16]  R. Marrett Aggregate properties of fracture populations , 1996 .

[17]  Vladimir Cvetkovic,et al.  Inference of field‐scale fracture transmissivities in crystalline rock using flow log measurements , 2010 .