Theory of transient streaming potentials associated with axial-symmetric flow in unconfined aquifers

SUMMARY We present a semi-analytical solution for the transient streaming potential response of an unconfined aquifer to continuous constant rate pumping. We assume that flow occurs without leakage from the unit below a transverse anisotropic aquifer and neglect flow in the unsaturated zone by treating the water-table as a moving material boundary. In the development of the solution to the streaming potential problem, we impose insulating boundary conditions at land surface and the lower boundary of the lower confining unit. We solve the problem exactly in the double Laplace–Hankel transform space and obtain the inverse transforms numerically. The solution is used to analyse transient streaming potential data collected during dipole hydraulic tests conducted at the Boise Hydrogeophysical Research Site in 2007 June. This analysis yields estimates of aquifer hydraulic parameters. The estimated hydraulic parameters, namely, hydraulic conductivity, transverse hydraulic anisotropy, specific storage and specific yield, compare well to published estimates obtained by inverting drawdown data collected at the field site.

[1]  W. Barrash,et al.  Significance of porosity for stratigraphy and textural composition in subsurface, coarse fluvial deposits: Boise Hydrogeophysical Research Site , 2004 .

[2]  Daniel M. Tartakovsky,et al.  Type curve interpretation of late‐time pumping test data in randomly heterogeneous aquifers , 2007 .

[3]  Thomas Wieder,et al.  Algorithm 794: numerical Hankel transform by the Fortran program HANKEL , 1999, TOMS.

[4]  Kristopher L. Kuhlman,et al.  A semi-analytical solution for transient streaming potentialsassociated with confined aquifer pumping tests , 2009 .

[5]  G. Petiau Second Generation of Lead-lead Chloride Electrodes for Geophysical Applications , 2000 .

[6]  Neven Kresic,et al.  Hydrogeology and groundwater modeling , 2006 .

[7]  André Revil,et al.  Streaming potentials of granular media: Influence of the Dukhin and Reynolds numbers , 2007 .

[8]  Salvatore Straface,et al.  Self‐potential signals associated with pumping tests experiments , 2004 .

[9]  André Revil,et al.  A Sandbox Experiment of Self‐Potential Signals Associated with a Pumping Test , 2004 .

[10]  Christophe Tournassat,et al.  Modeling the composition of the pore water in a clay-rock geological formation (Callovo-Oxfordian, France) , 2007 .

[11]  Y. Ilyin,et al.  Electrokinetic spontaneous polarization in porous media: petrophysics and numerical modelling , 2002 .

[12]  Reprint from the Proceedings of the 32nd Symposium on Engineering Geology and Geo- technical Engineering, March 26-28, 1997, Boise, ID. LITHOLOGIC, HYDROLOGIC AND PETROPHYSICAL CHARACTERIZATION OF AN UNCONSOLIDATED COBBLE-AND-SAND AQUIFER CAPITAL STATION SITE, BOISE, IDAHO , 1997 .

[13]  Jörg Renner,et al.  Self‐potential signals induced by periodic pumping tests , 2008 .

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

[15]  Xiaoye S. Li,et al.  A Comparison of Three High-Precision Quadrature Schemes , 2003, Exp. Math..

[16]  W. R. Sill,et al.  Self-potential modeling from primary flows , 1983 .

[17]  A Revil,et al.  A Potential‐Based Inversion of Unconfined Steady‐State Hydraulic Tomography , 2009, Ground water.

[18]  Tom Clemo,et al.  Field, laboratory, and modeling investigation of the skin effect at wells with slotted casing, Boise Hydrogeophysical Research Site , 2006 .

[19]  A. Ogilvy,et al.  DEFORMATIONS OF NATURAL ELECTRIC FIELDS NEAR DRAINAGE STRUCTURES , 1973 .

[20]  Pascal Sailhac,et al.  Analytic potentials for the forward and inverse modeling of SP anomalies caused by subsurface fluid flow , 2001 .

[21]  Masatake Mori,et al.  Double Exponential Formulas for Numerical Integration , 1973 .

[22]  André Revil,et al.  Principles of electrography applied to self‐potential electrokinetic sources and hydrogeological applications , 2003 .

[23]  Pascal Sailhac,et al.  Estimating aquifer hydraulic properties from the inversion of surface Streaming Potential (SP) anomalies , 2003 .

[24]  S. P. Neuman,et al.  Theory of flow in aquicludes adjacent to slightly leaky aquifers , 1968 .

[25]  L. Rosenhead Conduction of Heat in Solids , 1947, Nature.

[26]  Richard Chandler,et al.  Transient streaming potential measurements on fluid‐saturated porous structures: An experimental verification of Biot’s slow wave in the quasi‐static limit , 1980 .

[27]  Daniel M. Tartakovsky,et al.  Ergodicity of pumping tests , 2007 .

[28]  M. Voltz,et al.  Monitoring of an infiltration experiment using the self‐potential method , 2006 .

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

[30]  Robert W. Gillham,et al.  Unsaturated and Saturated Flow in Response to Pumping of an Unconfined Aquifer: Field Evidence of Delayed Drainage , 1992 .

[31]  R. Gillham,et al.  A Comparative Study of Specific Yield Determinations for a Shallow Sand Aquifer , 1984 .

[32]  P. Donaldson,et al.  Geophysical Surveys across the Boise Hydrogeophysical Research Site To Determine Geophysical Parameters of a Shallow, Alluvial Aquifer (1999) , 1999 .

[33]  S. Troisi,et al.  Numerical modelling of self-potential signals associated with a pumping test experiment , 2005 .

[34]  Allen F. Moench,et al.  Specific Yield as Determined by Type‐Curve Analysis of Aquifer‐Test Data , 1994 .

[35]  W. Barrash,et al.  Hierarchical geostatistics and multifacies systems: Boise Hydrogeophysical Research Site, Boise, Idaho , 2002 .

[36]  S. P. Neuman,et al.  Three‐dimensional saturated‐unsaturated flow with axial symmetry to a partially penetrating well in a compressible unconfined aquifer , 2007 .

[37]  Cass T. Miller,et al.  Reconstruction of the Water Table from Self‐Potential Data: A Bayesian Approach , 2009, Ground water.