Interpretation and inversion of low‐frequency head observations

[1] An asymptotic expression for transient head variations, valid at low-frequencies, forms the basis for an efficient scheme for estimating hydraulic conductivity. The critical computational step is equivalent to solving the governing equation for steady state head. Thus, model parameter sensitivities, relating changes in head to changes in hydraulic conductivity, of the fully transient problem can be computed with the equivalent of four steady state head computations. A comparison of model parameter sensitivities computed using the low-frequency asymptotic approach and sensitivities computed using a purely numerical approach indicates good agreement. An inversion of synthetic hydraulic tomography data indicates that it is possible to estimate overall permeability variations using the technique. In an actual application to truncated crosswell pressure tests from the Raymond field site, we image two high permeability fracture zones, in agreement with a conceptual model of the region. The location of the two fracture zones correlates with the position of transmissive fractures, as measured by borehole conductivity logs.

[1]  S. Jonathan Chapman,et al.  Ray Theory for High-Péclet-Number Convection-Diffusion , 1999, SIAM Journal on Applied Mathematics.

[2]  James J. Butler,et al.  Pumping tests in networks of multilevel sampling wells: Motivation and methodology , 1999 .

[3]  Michael A. Saunders,et al.  Algorithm 583: LSQR: Sparse Linear Equations and Least Squares Problems , 1982, TOMS.

[4]  J. Virieux,et al.  Asymptotic theory for diffusive electromagnetic imaging , 1994 .

[5]  T. Yeh,et al.  Hydraulic tomography: Development of a new aquifer test method , 2000 .

[6]  A. Gringarten,et al.  Use of Downhole Permanent Pressure Gauge Data to Diagnose Production Problems in a North Sea Horizontal Well , 2003 .

[7]  M. Kline,et al.  Electromagnetic theory and geometrical optics , 1965 .

[8]  S. Finsterle,et al.  Numerical trajectory calculations for the efficient inversion of transient flow and tracer observations , 2004 .

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

[10]  Joseph B. Keller,et al.  ASYMPTOTIC METHODS FOR PARTIAL DIFFERENTIAL EQUATIONS: THE REDUCED WAVE EQUATION AND MAXWELL'S EQUATION , 1995 .

[11]  A. Cohen.  Hydrogeologic characterization of a fractured granitic rock aquifer, Raymond, California , 1993 .

[12]  A. Tikhonov.  Mathematical basis of the theory of electromagnetic soundings , 1965 .

[13]  K. Karasaki,et al.  A multidisciplinary fractured rock characterization study at Raymond field site, Raymond, CA. , 2000 .

[14]  R. Bracewell.  The Fourier Transform and Its Applications , 1966 .

[15]  David L. Freyberg,et al.  A natural gradient experiment on solute transport in a sand aquifer: 2. Spatial moments and the advection and dispersion of nonreactive tracers , 1986 .

[16]  James J. Butler,et al.  Pumping tests in nonuniform aquifers — The radially symmetric case , 1988 .

[17]  A. Datta-Gupta,et al.  Rapid inverse modeling of pressure interference tests using trajectory‐based traveltime and amplitude sensitivities , 2006 .

[18]  Joseph B. Keller,et al.  Asymptotic expansions of solutions of (▽2+k2)u=0 , 1955 .

[19]  Don W. Vasco,et al.  Inversion of pressure observations: an integral formulation , 2001 .

[20]  W. Durand.  Dynamics of Fluids , 1934 .

[21]  R. Parker.  Geophysical Inverse Theory , 1994 .

[22]  S. P. Neuman,et al.  Field Determination of the Three-Dimensional Hydraulic Conductivity Tensor of Anisotropic Media: 2. Methodology and Application to Fractured Rocks , 1985 .

[23]  Frederick L. Paillet,et al.  Using borehole geophysics and cross-borehole flow testing to define hydraulic connections between fracture zones in bedrock aquifers , 1993 .

[24]  Poroelasticity: observations of anomalous near surface tilt induced by ground water pumping , 2003 .

[25]  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 .

[26]  Michael A. Saunders,et al.  LSQR: An Algorithm for Sparse Linear Equations and Sparse Least Squares , 1982, TOMS.

[27]  Jack K. Cohen,et al.  A Ray Method for the Asymptotic Solution of the Diffusion Equation , 1967 .

[28]  K. Karasaki,et al.  Estimation of reservoir properties using transient pressure data: An asymptotic approach , 2000 .

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

[30]  James J. Butler,et al.  Steady shape analysis of tomographic pumping tests for characterization of aquifer heterogeneities , 2002 .

[31]  F. B..  Introduction to Bessel Functions , 1939, Nature.

[32]  M. A. Pérez-Flores,et al.  Forward modelling of direct current and low-frequency electromagnetic fields using integral equations , 2002 .

[33]  Ronald Smith.  The early stages of contaminant dispersion in shear flows , 1981, Journal of Fluid Mechanics.

[34]  D. Vasco,et al.  A simultaneous inversion of seismic traveltimes and amplitudes for velocity and attenuation , 1996 .

[35]  Velimir V. Vesselinov,et al.  Three‐dimensional numerical inversion of pneumatic cross‐hole tests in unsaturated fractured tuff: 1. Methodology and borehole effects , 2001 .

[36]  D. Vasco.  An algebraic formulation of geophysical inverse problems , 2000 .

[37]  D. Vasco.  Intersections, ideals, and inversion , 1999 .

[38]  D. Vasco,et al.  Estimation of flow properties using surface deformation and head data: A trajectory‐based approach , 2004 .

[39]  K. Karasaki,et al.  Sensitivity Study on Hydraulic Well Testing Inversion Using Simulated Annealing , 1997, Ground water.

[40]  R. Dingle.  Asymptotic expansions : their derivation and interpretation , 1975 .

[41]  Dean S. Oliver,et al.  The influence of nonuniform transmissivity and storativity on drawdown , 1993 .