Improved Characterization of DNAPL Source Zones via Sequential Hydrogeophysical Inversion of Hydraulic‐Head, Self‐Potential and Partitioning Tracer Data

High‐resolution characterization of hydraulic properties and dense nonaqueous phase liquid (DNAPL) contaminant source is crucial to develop efficient remediation strategies. However, DNAPL characterization suffers from a limited number of borehole data in the field, resulting in a low‐resolution estimation. Moreover, high‐resolution DNAPL characterization requires a large number of unknowns to be estimated, presenting a computational bottleneck. In this paper, a low‐cost geophysical approach, the self‐potential method, is used as additional information for hydraulic properties characterization. Joint inversion of hydraulic head and self‐potential measurements is proposed to improve hydraulic conductivity estimation, which is then used to characterize the DNAPL saturation distribution by inverting partitioning tracer measurements. The computational barrier is overcome by (a) solving the inversion by the principal component geostatistical approach, in which the covariance matrix is replaced by a low‐rank approximation, thus reducing the number of forward model runs; (b) using temporal moments of concentrations instead of individual concentration data points for faster forward simulations. To assess the ability of the proposed approach, numerical experiments are conducted in a 3‐D aquifer with 104 unknown hydraulic conductivities and DNAPL saturations. Results show that with realistic DNAPL sources and a limited number of hydraulic heads, the traditional hydraulic/partitioning tracer tomography roughly reconstructs subsurface heterogeneity but fails to resolve the DNAPL distribution. By adding self‐potential data, the error is reduced by 24% in hydraulic conductivity estimation and 68% in DNAPL saturation characterization. The proposed sequential inversion framework utilizes the complementary information from multi‐source hydrogeophysical data sets, and can provide high‐resolution characterizations for realistic DNAPL sources.

[1]  R. H. Brooks,et al.  Hydraulic properties of porous media , 1963 .

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

[3]  Jack C. Parker,et al.  A model for hysteretic constitutive relations governing multiphase flow: 3. Refinements and numerical simulations , 1989 .

[4]  Peter K. Kitanidis,et al.  Orthonormal residuals in geostatistics: Model criticism and parameter estimation , 1991 .

[5]  E. Eric Adams,et al.  Field study of dispersion in a heterogeneous aquifer: 1. Overview and site description , 1992 .

[6]  C. R. Dietrich,et al.  A fast and exact method for multidimensional gaussian stochastic simulations , 1993 .

[7]  Charles F. Harvey,et al.  Temporal Moment‐Generating Equations: Modeling Transport and Mass Transfer in Heterogeneous Aquifers , 1995 .

[8]  Minghui Jin,et al.  AN ITERATIVE STOCHASTIC INVERSE METHOD: CONDITIONAL EFFECTIVE TRANSMISSIVITY AND HYDRAULIC HEAD FIELDS , 1995 .

[9]  Kamy Sepehrnoori,et al.  Partitioning Tracer Test for Detection, Estimation, and Remediation Performance Assessment of Subsurface Nonaqueous Phase Liquids , 1995 .

[10]  P. Kitanidis.  Quasi‐Linear Geostatistical Theory for Inversing , 1995 .

[11]  Mark L. Brusseau,et al.  Field Study of the Partitioning Tracer Method for Detection of Dense Nonaqueous Phase Liquid in a Trichloroethene- Contaminated Aquifer , 1996 .

[12]  John A. Cherry,et al.  Dense Chlorinated Solvents and other DNAPLs in Groundwater , 1996 .

[13]  Wendy D. Graham,et al.  Optimal estimation of residual non–aqueous phase liquid saturations using partitioning tracer concentration data , 1997 .

[14]  Brian Berkowitz,et al.  A generalized growth model for simulating initial migration of dense non‐aqueous phase liquids , 1998 .

[15]  M. Grismer,et al.  Kinetics of volatile organic compound sorption/desorption on clay minerals , 1999 .

[16]  W. Graham,et al.  Estimation of spatially variable residual nonaqueous phase liquid saturations in nonuniform flow fields using partitioning tracer data , 2000 .

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

[18]  Linda M. Abriola,et al.  The influence of field-scale heterogeneity on the infiltration and entrapment of dense nonaqueous phase liquids in saturated formations , 2000 .

[19]  Thomas C. Harmon,et al.  Inverse modeling for locating dense nonaqueous pools in groundwater under steady flow conditions , 2000 .

[20]  Peter K. Kitanidis,et al.  Sensitivity of temporal moments calculated by the adjoint-state method and joint inversing of head and tracer data , 2000 .

[21]  L. Yarrington,et al.  Gravity‐destabilized nonwetting phase invasion in macroheterogeneous porous media: Near‐pore‐scale macro modified invasion percolation simulation of experiments , 2001 .

[22]  A. Binley,et al.  Vadose zone flow model parameterisation using cross-borehole radar and resistivity imaging , 2001 .

[23]  Michelle Sneed.  Hydraulic and mechanical properties affecting ground-water flow and aquifer-system compaction, San Joaquin Valley, California , 2001 .

[24]  Spatial characterization of a hydrogeochemically heterogeneous aquifer using partitioning tracers: Optimal estimation of aquifer parameters , 2001 .

[25]  G. Pope,et al.  Inverse modeling of partitioning interwell tracer tests: A streamline approach , 2002 .

[26]  David P. Lesmes,et al.  Electrical‐hydraulic relationships observed for unconsolidated sediments , 2002 .

[27]  André Revil,et al.  Relationship between self‐potential (SP) signals and redox conditions in contaminated groundwater , 2003 .

[28]  Tian-Chyi J. Yeh,et al.  Characterization of aquifer heterogeneity using transient hydraulic tomography , 2004 .

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

[30]  P. Kitanidis,et al.  Estimation of historical groundwater contaminant distribution using the adjoint state method applied to geostatistical inverse modeling , 2004 .

[31]  E. Sudicky,et al.  Transition Probability/Markov Chain Analyses of DNAPL Source Zones and Plumes , 2006, Ground water.

[32]  T. Illangasekare,et al.  Influence of dense nonaqueous phase liquid pool morphology on the performance of partitioning tracer tests: Evaluation of the equilibrium assumption , 2006 .

[33]  Quantifying PCE and TCE in DNAPL Source Zones: Effects of Sampling Methods Used for Intact Cores at Varied Contaminant Levels and Media Temperatures , 2006 .

[34]  Charles J. Newell,et al.  Multiyear Temporal Changes in Chlorinated Solvent Concentrations at 23 Monitored Natural Attenuation Sites , 2006 .

[35]  Kurt D. Pennell,et al.  Estimating mass discharge from dense nonaqueous phase liquid source zones using upscaled mass transfer coefficients: An evaluation using multiphase numerical simulations , 2006 .

[36]  W. Nowak,et al.  Geostatistical inference of hydraulic conductivity and dispersivities from hydraulic heads and tracer data , 2006 .

[37]  Burke J. Minsley,et al.  Three‐dimensional self‐potential inversion for subsurface DNAPL contaminant detection at the Savannah River Site, South Carolina , 2007 .

[38]  André Revil,et al.  Tomography of the Darcy velocity from self‐potential measurements , 2007 .

[39]  Tian-Chyi J. Yeh,et al.  Hydraulic/partitioning tracer tomography for characterization of dense nonaqueous phase liquid source zones , 2007 .

[40]  A. Hördt,et al.  Hydraulic conductivity estimation from induced polarisation data at the field scale — the Krauthausen case history , 2007 .

[41]  S. Finsterle,et al.  Electrokinetic coupling in unsaturated porous media. , 2007, Journal of colloid and interface science.

[42]  B. Parker,et al.  Plume persistence caused by back diffusion from thin clay layers in a sand aquifer following TCE source-zone hydraulic isolation. , 2008, Journal of contaminant hydrology.

[43]  A. Revil,et al.  Stochastic joint inversion of temperature and self-potential data to determine the permeability of geological units in sedimentary basins , 2009 .

[44]  Cheng-Haw Lee,et al.  River stage tomography: A new approach for characterizing groundwater basins , 2009 .

[45]  Hiromitsu Saegusa,et al.  Hydraulic tomography in fractured granite: Mizunami Underground Research site, Japan , 2009 .

[46]  A. Villa,et al.  Spectral induced polarization for the characterization of free-phase hydrocarbon contamination of sediments with low clay content , 2009 .

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

[48]  Junfeng Zhu,et al.  Analysis of tracer tomography using temporal moments of tracer breakthrough curves , 2009 .

[49]  Roelof Versteeg,et al.  Improved hydrogeophysical characterization and monitoring through parallel modeling and inversion of time-domain resistivity andinduced-polarization data , 2010 .

[50]  James J. Butler,et al.  Inherent Limitations of Hydraulic Tomography , 2010, Ground water.

[51]  Olaf A. Cirpka,et al.  Fully coupled hydrogeophysical inversion of synthetic salt tracer experiments , 2010 .

[52]  André Revil,et al.  Influence of oil saturation upon spectral induced polarization of oil-bearing sands , 2010 .

[53]  Walter A Illman,et al.  Hydraulic/partitioning tracer tomography for DNAPL source zone characterization: small-scale sandbox experiments. , 2010, Environmental science & technology.

[54]  Johan Alexander Huisman,et al.  Hydraulic properties of a model dike from coupled Bayesian and multi-criteria hydrogeophysical inversion , 2010 .

[55]  J. Chambers,et al.  Hydrogeophysical imaging of deposit heterogeneity and groundwater chemistry changes during DNAPL source zone bioremediation. , 2010, Journal of contaminant hydrology.

[56]  M. Riva,et al.  Joint Inversion of Steady-State Hydrologic and Self-Potential Data for 3D Hydraulic Conductivity Distribution at the Boise Hydrogeophysical Research Site , 2011 .

[57]  M. Cardiff,et al.  3‐D transient hydraulic tomography in unconfined aquifers with fast drainage response , 2011 .

[58]  Tian-Chyi J. Yeh,et al.  Robustness of joint interpretation of sequential pumping tests: Numerical and field experiments , 2011 .

[59]  Peter K. Kitanidis,et al.  A field proof‐of‐concept of aquifer imaging using 3‐D transient hydraulic tomography with modular, temporarily‐emplaced equipment , 2012 .

[60]  Olaf A. Cirpka,et al.  Fully coupled hydrogeophysical inversion of a laboratory salt tracer experiment monitored by electrical resistivity tomography , 2012 .

[61]  H. Vereecken,et al.  Coupled hydrogeophysical inversion of electrical resistances and inflow measurements for topsoil hydraulic properties under constant head infiltration , 2012 .

[62]  Michael A. Mooney,et al.  Saline pulse test monitoring with the self‐potential method to nonintrusively determine the velocity of the pore water in leaking areas of earth dams and embankments , 2012 .

[63]  Niklas Linde,et al.  Derivation of Soil‐Specific Streaming Potential Electrical Parameters from Hydrodynamic Characteristics of Partially Saturated Soils , 2012, 1209.6335.

[64]  M. Zamora,et al.  Laboratory study of self-potential signals during releasing of CO2 and N2 plumes , 2012 .

[65]  A. Revil,et al.  Coupled hydromechanical and electromagnetic disturbances in unsaturated porous materials , 2013, Water resources research.

[66]  E. Miller,et al.  A geometric approach to joint inversion with applications to contaminant source zone characterization , 2013, 1303.5109.

[67]  Yonghong Hao,et al.  A temporal sampling strategy for hydraulic tomography analysis , 2013 .

[68]  A. Revil,et al.  Stochastic joint inversion of hydrogeophysical data for salt tracer test monitoring and hydraulic conductivity imaging , 2013 .

[69]  A. Revil,et al.  Effective conductivity and permittivity of unsaturated porous materials in the frequency range 1 mHz–1GHz , 2013, Water resources research.

[70]  P. Kitanidis,et al.  Aquifer heterogeneity characterization with oscillatory pumping: Sensitivity analysis and imaging potential , 2013 .

[71]  Ivan Lunati,et al.  Challenges in modeling unstable two‐phase flow experiments in porous micromodels , 2014 .

[72]  André Revil,et al.  Hydraulic conductivity field characterization from the joint inversion of hydraulic heads and self‐potential data , 2014 .

[73]  Peter K. Kitanidis,et al.  Large‐scale hydraulic tomography and joint inversion of head and tracer data using the Principal Component Geostatistical Approach (PCGA) , 2014 .

[74]  T. Goto,et al.  Self-potential Inversion for the Estimation of Permeability Structure , 2013 .

[75]  C. H. Ward,et al.  Chlorinated Solvent Source Zone Remediation , 2014 .

[76]  P. Kitanidis,et al.  Principal Component Geostatistical Approach for large-dimensional inverse problems , 2014, Water resources research.

[77]  J. Gerhard,et al.  Improved time-lapse electrical resistivity tomography monitoring of dense non-aqueous phase liquids with surface-to-horizontal borehole arrays , 2015 .

[78]  Jichun Wu,et al.  Influence of flow velocity and spatial heterogeneity on DNAPL migration in porous media: insights from laboratory experiments and numerical modelling , 2015, Hydrogeology Journal.

[79]  Niklas Linde,et al.  Monitoring of saline tracer movement with vertically distributed self-potential measurements at the HOBE agricultural test site, Voulund, Denmark , 2015, 1501.02277.

[80]  Wolfgang Nowak,et al.  Predicting DNAPL mass discharge and contaminated site longevity probabilities: Conceptual model and high‐resolution stochastic simulation , 2015 .

[81]  G. Fiandaca,et al.  Influence of non-aqueous phase liquid configuration on induced polarization parameters: Conceptual models applied to a time-domain field case study , 2015 .

[82]  Bernd Kulessa,et al.  Resistivity and self-potential tomography applied to groundwater remediation and contaminant plumes: Sandbox and field experiments , 2015 .

[83]  Lucy S. Oxby,et al.  Geoelectrical monitoring of simulated subsurface leakage to support high-hazard nuclear decommissioning at the Sellafield Site, UK. , 2016, The Science of the total environment.

[84]  André Revil,et al.  Joint inversion of hydraulic head and self‐potential data associated with harmonic pumping tests , 2016 .

[85]  Peter K. Kitanidis,et al.  Scalable subsurface inverse modeling of huge data sets with an application to tracer concentration breakthrough data from magnetic resonance imaging , 2016 .

[86]  W. Nowak,et al.  Identification of contaminant source architectures—A statistical inversion that emulates multiphase physics in a computationally practicable manner , 2016 .

[87]  André Revil,et al.  Specific storage and hydraulic conductivity tomography through the joint inversion of hydraulic heads and self-potential data , 2016 .

[88]  Peter K. Kitanidis,et al.  Imaging geochemical heterogeneities using inverse reactive transport modeling: An example relevant for characterizing arsenic mobilization and distribution , 2016 .

[89]  Jichun Wu,et al.  Quantitative assessment of electrical resistivity tomography for monitoring DNAPLs migration – Comparison with high-resolution light transmission visualization in laboratory sandbox , 2017 .

[90]  A. Revil,et al.  Redox potential distribution of an organic-rich contaminated site obtained by the inversion of self-potential data , 2017 .

[91]  G. Fiandaca,et al.  In fl uence of non-aqueous phase liquid con fi guration on induced polarization parameters : Conceptual models applied to a time-domain fi eld case study , 2017 .

[92]  A. Revil.  Transport of water and ions in partially water-saturated porous media. Part 1. Constitutive equations , 2017 .

[93]  Multidimensional Investigation of Bedrock Heterogeneity/Unconformities at a DNAPL‐Impacted Site , 2017, Ground water.

[94]  P. Kitanidis,et al.  Fast Large-Scale Joint Inversion for Deep Aquifer Characterization Using Pressure and Heat Tracer Measurements , 2018, Transport in Porous Media.

[95]  Ruben Juanes,et al.  Improved characterization of heterogeneous permeability in saline aquifers from transient pressure data during freshwater injection , 2017 .

[96]  Eric Darve,et al.  Riverine Bathymetry Imaging With Indirect Observations , 2017 .

[97]  James Irving,et al.  On uncertainty quantification in hydrogeology and hydrogeophysics , 2017 .

[98]  M. Cardiff,et al.  Oscillatory hydraulic testing as a strategy for NAPL source zone monitoring: Laboratory experiments. , 2017, Journal of contaminant hydrology.

[99]  M. Jackson,et al.  Self‐Potential as a Predictor of Seawater Intrusion in Coastal Groundwater Boreholes , 2018, Water Resources Research.

[100]  Jichun Wu,et al.  Coupled hydrogeophysical inversion of DNAPL source zone architecture and permeability field in a 3D heterogeneous sandbox by assimilation time-lapse cross-borehole electrical resistivity data via ensemble Kalman filtering , 2018, Journal of Hydrology.

[101]  N. Linde,et al.  Impact of petrophysical uncertainty on Bayesian hydrogeophysical inversion and model selection , 2018 .

[102]  A. Jardani,et al.  3D geostatistical inversion of induced polarization data and its application to coal seam fires , 2018 .

[103]  M. Jackson,et al.  Remote Detection of Saline Intrusion in a Coastal Aquifer Using Borehole Measurements of Self‐Potential , 2018 .

[104]  André Revil,et al.  The Self-Potential Method: Theory and Applications in Environmental Geosciences , 2018 .

[105]  Y. Zha,et al.  Exploitation of pump-and-treat remediation systems for characterization of hydraulic heterogeneity , 2019, Journal of Hydrology.

[106]  A. Revil,et al.  Induced polarization of volcanic rocks. 3. Imaging clay cap properties in geothermal fields , 2019, Geophysical Journal International.

[107]  Jichun Wu,et al.  Coupled hydrogeophysical inversion to identify non-Gaussian hydraulic conductivity field by jointly assimilating geochemical and time-lapse geophysical data , 2019, Journal of Hydrology.

[108]  V. Guihard,et al.  Induced polarization as a tool to non-intrusively characterize embankment hydraulic properties , 2020, Engineering Geology.