Efficient multiscale imaging of subsurface resistivity with uncertainty quantification using ensemble Kalman inversion
Tóm tắt
Electrical resistivity tomography (ERT) is widely used to image the Earth’s subsurface and has proven to be an extremely useful tool in application to hydrological problems. Conventional smoothness-constrained inversion of ERT data is efficient and robust, and consequently very popular. However, it does not resolve well sharp interfaces of a resistivity field and tends to reduce and smooth resistivity variations. These issues can be problematic in a range of hydrological or near-surface studies, for example mapping regolith-bedrock interfaces. While fully Bayesian approaches, such as those using Markov chain Monte Carlo sampling, can address the above issues, their very high computation cost makes them impractical for many applications. Ensemble Kalman inversion (EKI) offers a computationally efficient alternative by approximating the Bayesian posterior distribution in a derivative-free manner, which means only a relatively small number of ‘black-box’ model runs are required. Although common limitations for ensemble Kalman filter-type methods apply to EKI, it is both efficient and generally captures uncertainty patterns correctly. We propose the use of a new EKI-based framework for ERT which estimates a resistivity model and its uncertainty at a modest computational cost. Our EKI framework uses a level-set parametrization of the unknown resistivity to allow efficient estimation of discontinuous resistivity fields. Instead of estimating level-set parameters directly, we introduce a second step to characterize the spatial variability of the resistivity field and infer length scale hyperparameters directly. We demonstrate these features by applying the method to a series of synthetic and field examples. We also benchmark our results by comparing them to those obtained from standard smoothness-constrained inversion. Resultant resistivity images from EKI successfully capture arbitrarily shaped interfaces between resistivity zones and the inverted resistivities are close to the true values in synthetic cases. We highlight its readiness and applicability to similar problems in geophysics.
Từ khóa
Tài liệu tham khảo
Aalstad, 2018, Ensemble-based assimilation of fractional snow-covered area satellite retrievals to estimate the snow distribution at Arctic sites, Cryosphere, 12, 247, 10.5194/tc-12-247-2018
Aghasi, 2011, Parametric level set methods for inverse problems, SIAM J. Imag. Sci., 4, 618, 10.1137/100800208
Andersen, 2003, Bayesian inversion of geoelectrical resistivity data, J. R. Stat. Soc., B, 65, 619, 10.1111/1467-9868.00406
Arridge, 2019, Solving inverse problems using data-driven models, Acta Numer., 28, 1, 10.1017/S0962492919000059
Asunbo, 2007, Improving the hydrological characterization of peatlands using geophysics, Master’s thesis
Auken, 2004, Layered and laterally constrained 2D inversion of resistivity data, Geophysics, 69, 752, 10.1190/1.1759461
Baù, 2015, Ensemble smoothing of land subsidence measurements for reservoir geomechanical characterization, Int. J. Numer. Anal. Methods Geomech., 39, 207, 10.1002/nag.2309
Binley, 2015, Tools and techniques: electrical methods, Treatise on Geophysics, 2nd end, 233, 10.1016/B978-0-444-53802-4.00192-5
Binley, 2020, Resistivity and Induced Polarization: Theory and Applications to the Near-Surface Earth, 10.1017/9781108685955
Binley, 2002, Vadose zone flow model parameterisation using cross-borehole radar and resistivity imaging, J. Hydrol., 267, 147, 10.1016/S0022-1694(02)00146-4
Binley, 2004, Characterization of heterogeneity in unsaturated sandstone using borehole logs and cross-borehole tomography, Aquifer Characterization, 129, 10.2110/pec.04.80.0129
Binley, 2005, Relationship between spectral induced polarization and hydraulic properties of saturated and unsaturated sandstone, Water Resour. Res., 41, W12417, 10.1029/2005WR004202
Blanchy, 2020, ResIPy, an intuitive open source software for complex geoelectrical inversion/modeling, Comp. Geosci., 137, 104423, 10.1016/j.cageo.2020.104423
Blanchy, 2020, Timeintensive geoelectrical monitoring under winter wheat, Near Surf. Geophys., 18, 413, 10.1002/nsg.12107
Bouchedda, 2012, 2-D joint structural inversion of cross-hole electrical resistance and ground penetrating radar data, J. appl. Geophys., 78, 52, 10.1016/j.jappgeo.2011.10.009
Bouchedda, 2017, Constrained electrical resistivity tomography Bayesian inversion using inverse Matérn covariance matrix, Geophysics, 82, E129, 10.1190/geo2015-0673.1
Calvetti, 2018, Inverse problems: From regularization to Bayesian inference, WIREs Comput. Stat., 10, 1, 10.1002/wics.1427
Camporese, 2015, Coupled and uncoupled hydrogeophysical inversions using ensemble Kalman filter assimilation of ERT-monitored tracer test data, Water Resour. Res., 51, 3277, 10.1002/2014WR016017
Canchumuni, 2019, Towards a robust parameterization for conditioning facies models using deep variational autoencoders and ensemble smoother, Comp. Geosci., 128, 87, 10.1016/j.cageo.2019.04.006
Cardiff, 2009, Bayesian inversion for facies detection: an extensible level set framework, Water Resour. Res., 45, 1, 10.1029/2008WR007675
Cassiani, 2005, Modeling unsaturated flow in a layered formation under quasi-steady state conditions using geophysical data constraints, Adv. Water Resour., 28, 467, 10.1016/j.advwatres.2004.12.007
Chada, 2018, Parameterizations for ensemble Kalman inversion, Inverse Prob., 34, 055009, 10.1088/1361-6420/aab6d9
Chambers, 2014, Derivation of lowland riparian wetland deposit architecture using geophysical image analysis and interface detection, Water Resour. Res., 50, 5886, 10.1002/2014WR015643
Chan, 2004, Level set and total variation regularization for elliptic inverse problems with discontinuous coefficients, J. Comput. Phys., 193, 40, 10.1016/j.jcp.2003.08.003
Chen, 2015, Using the level set method to study the effects of heterogeneity and anisotropy on hyporheic exchange, Water Resour. Res., 51, 3617, 10.1002/2014WR016444
Chen, 2013, Application of ensemble-based data assimilation techniques for aquifer characterization using tracer data at Hanford 300 area, Water Resour. Res., 49, 7064, 10.1002/2012WR013285
Cheng, 2019, Characterization of karst structures using quasi-3D electrical resistivity tomography, Environ. Earth Sci., 78, 10.1007/s12665-019-8284-2
Chou, 2016, Estimation of saturated hydraulic conductivity during infiltration test with the aid of ERT and level-set method, Vadose Zone J., 15, vzj2015.05.0082, 10.2136/vzj2015.05.0082
Chung, 2005, Electrical impedance tomography using level set representation and total variational regularization, J. Comput. Phys., 205, 357, 10.1016/j.jcp.2004.11.022
Claes, 2020, Parameterization of a hydrologic model with geophysical data to simulate observed subsurface return flow paths, Vadose Zone J., 19, e20024, 10.1002/vzj2.20024
Comas, 2004, Low-frequency electrical properties of peat, Water Resour. Res., 40, 10.1029/2004WR003534
Comas, 2004, Geophysical evidence for peat basin morphology and stratigraphic controls on vegetation observed in a Northern Peatland, J. Hydrol., 295, 173, 10.1016/j.jhydrol.2004.03.008
de Pasquale, 2019, Probabilistic inference of subsurface heterogeneity and interface geometry using geophysical data, Geophys. J. Int., 217, 816
Dorn, 2008, History matching of petroleum reservoirs using a level set technique, Inverse Prob., 24, 035015, 10.1088/0266-5611/24/3/035015
Dunlop, 2017, Hierarchical Bayesian level set inversion, Stat. Comput., 27, 1555, 10.1007/s11222-016-9704-8
Emerick, 2016, Analysis of the performance of ensemble-based assimilation of production and seismic data, J. Petrol. Sci. Eng., 139, 219, 10.1016/j.petrol.2016.01.029
Emerick, 2013, Ensemble smoother with multiple data assimilation, Comput. Geosci., 55, 3, 10.1016/j.cageo.2012.03.011
Englert, 2016, Comparison of smoothness-constrained and geostatistically based cross-borehole electrical resistivity tomography for characterization of solute tracer plumes, Water Sci. Eng., 9, 274, 10.1016/j.wse.2017.01.002
Fernández-Álvarez, 2017, Combined GPR and ERT exploratory geophysical survey of the Medieval Village of Pancorbo Castle (Burgos, Spain), J. Appl. Geophys., 144, 86, 10.1016/j.jappgeo.2017.07.002
Galetti, 2018, Transdimensional electrical resistivity tomography, J. geophys. Res., 123, 6347, 10.1029/2017JB015418
Hermans, 2017, Facies discrimination with electrical resistivity tomography using a probabilistic methodology: effect of sensitivity and regularization, Near Surf. Geophys., 15, 13, 10.3997/1873-0604.2016047
Hermans, 2020, Combined inversion of land and marine electrical resistivity tomography for submarine groundwater discharge and saltwater intrusion characterization, Geophys. Res. Lett., 47, 10.1029/2019GL085877
Hermans, 2012, Imaging artificial salt water infiltration using electrical resistivity tomography constrained by geostatistical data, J. Hydrol., 438–439, 168, 10.1016/j.jhydrol.2012.03.021
Holmes, 2020, Four-dimensional electrical resistivity tomography for continuous, near-real-time monitoring of a landslide affecting transport infrastructure in British Columbia, Canada, Near Surf. Geophys., 18, 337, 10.1002/nsg.12102
Iglesias, 2016, A Bayesian level set method for geometric inverse problems, Interfaces Free Bound., 18, 181, 10.4171/IFB/362
Iglesias, 2018, Bayesian inversion in resin transfer molding, Inverse Prob., 34, 105002, 10.1088/1361-6420/aad1cc
Iglesias, 2016, A regularizing iterative ensemble Kalman method for PDE-constrained inverse problems, Inverse Prob., 32, 025002, 10.1088/0266-5611/32/2/025002
Iglesias, 2021, Adaptive regularisation for ensemble Kalman inversion, Inverse Prob., 37, 025008, 10.1088/1361-6420/abd29b
Iglesias, 2014, Well-posed Bayesian geometric inverse problems arising in subsurface flow, Inverse Prob., 30, 114001, 10.1088/0266-5611/30/11/114001
Irving, 2010, Stochastic inversion of tracer test and electrical geophysical data to estimate hydraulic conductivities, Water Resour. Res., 46, W11514, 10.1029/2009WR008340
Johnson, 2010, Improved hydrogeophysical characterization and monitoring through parallel modeling and inversion of time-domain resistivity and induced-polarization data, Geophysics, 75, WA27, 10.1190/1.3475513
Johnson, 2012, Characterization of a contaminated wellfield using 3D electrical resistivity tomography implemented with geostatistical, discontinuous boundary, and known conductivity constraints, Geophysics, 77, EN85, 10.1190/geo2012-0121.1
Ju, 2018, An adaptive Gaussian process-based iterative ensemble smoother for data assimilation, Adv. Water Resour., 115, 125, 10.1016/j.advwatres.2018.03.010
Kang, 2018, 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, J. Hydrol., 567, 149, 10.1016/j.jhydrol.2018.10.019
Kang, 2019, Coupled hydrogeophysical inversion to identify non-Gaussian hydraulic conductivity field by jointly assimilating geochemical and time-lapse geophysical data, J. Hydrol., 578, 10.1016/j.jhydrol.2019.124092
Kitanidis, 2015, Three-dimensional ERT imaging by the geostatistical approach, Proceedings of the American Geophysical Union Fall Meeting
Köpke, 2019, Hydrogeophysical parameter estimation using iterative ensemble smoothing and approximate forward solvers, Front. Environ. Sci., 7, 34, 10.3389/fenvs.2019.00034
Kuras, 2016, Geoelectrical monitoring of simulated subsurface leakage to support high-hazard nuclear decommissioning at the Sellafield Site, UK, Sci. Total Environ., 566-567, 350, 10.1016/j.scitotenv.2016.04.212
Labrecque, 2001, Difference inversion of ERT data: a fast inversion method for 3-D in situ monitoring, J. Environ. Eng. Geophys., 6, 83, 10.4133/JEEG6.2.83
Lan, 2018, Joint inversion of physical and geochemical parameters in groundwater models by sequential ensemble-based optimal design, Stoch. Environ. Res. Risk Assess., 32
Le, 2016, An adaptive ensemble smoother with multiple data assimilation for assisted history matching, SPE J., 21, 2195, 10.2118/173214-PA
Linde, 2006, Improved hydrogeophysical characterization using joint inversion of cross-hole electrical resistance and ground-penetrating radar traveltime data, Water Resour. Res., 42, WR005131, 10.1029/2006WR005131
Linde, 2017, On uncertainty quantification in hydrogeology and hydrogeophysics, Adv. Water Resour., 110, 166, 10.1016/j.advwatres.2017.10.014
Lindgren, 2011, An explicit link between Gaussian fields and Gaussian Markov random fields: the stochastic partial differential equation approach, J. R. Stat. Soc., B, 73, 423, 10.1111/j.1467-9868.2011.00777.x
Liu, 2018, Stochastic nonlinear inversion of seismic data for the estimation of petroelastic properties using the ensemble smoother and data reparameterization, Geophysics, 83, M25, 10.1190/geo2017-0713.1
Loke, 2003, A comparison of smooth and blocky inversion methods in 2D electrical imaging surveys, Explor. Geophys., 34, 182, 10.1071/EG03182
McLachlan, 2020, Electrical resistivity monitoring of river-groundwater interactions in a chalk river and neighbouring riparian zone, Near Surf. Geophys., 18, 385, 10.1002/nsg.12114
Muir, 2020, Geometric and level set tomography using ensemble Kalman inversion, Geophys. J. Int., 220, 967, 10.1093/gji/ggz472
Nguyen, 2016, Data-driven selection of the minimum-gradient support parameter in time-lapse focused electric imaging, Geophysics, 81, A1, 10.1190/geo2015-0226.1
Nussbaumer, 2019, Simulation of fine-scale electrical conductivity fields using resolution-limited tomograms and area-to-point kriging, Geophys. J. Int., 218, 1322, 10.1093/gji/ggz185
Ramirez, 2005, Stochastic inversion of electrical resistivity changes using a Markov Chain Monte Carlo approach, J. geophys. Res., 110, B02101, 10.1029/2004JB003449
Roininen, 2014, Sparse approximations of fractional Matérn fields, 45, 194
Ronczka, 2020, Monitoring freshwater–saltwater interfaces with SAMOS–installation effects on data and inversion, Near Surf. Geophys., 18, 369, 10.1002/nsg.12115
Sassen, 2012, Reactive facies: an approach for parameterizing field-scale reactive transport models using geophysical methods, Water Resour. Res., 48, 1, 10.1029/2011WR011047
Scholer, 2011, Estimating vadose zone hydraulic properties using ground penetrating radar: the impact of prior information, Water Resour. Res., 47, 1, 10.1029/2011WR010409
Sebacher, 2017, An adaptive Pluri-Gaussian simulation model for geological uncertainty quantification, J. Petrol. Sci. Eng., 158, 494, 10.1016/j.petrol.2017.08.038
Si, 2015, TetGen, a Delaunay-based quality tetrahedral mesh generator, ACM Trans. Math. Software, 41, 1, 10.1145/2629697
Slater, 2003, Evaluation of permeable reactive barrier (PRB) integrity using electrical imaging methods, Geophysics, 68, 911, 10.1190/1.1581043
Slater, 2002, Investigating peatland stratigraphy and hydrogeology using integrated electrical geophysics, Geophysics, 67, 365, 10.1190/1.1468597
Song, 2019, Delineating facies spatial distribution by integrating ensemble data assimilation and indicator geostatistics with levelset transformation, Water Resour. Res., 55, 2652, 10.1029/2018WR023262
Tso, 2019, Enhancing the information content of geophysical data for nuclear site characterisation, PhD thesis
Tso, 2017, Improved characterisation and modelling of measurement errors in electrical resistivity tomography (ERT) surveys, J. appl. Geophys., 146, 103, 10.1016/j.jappgeo.2017.09.009
Tso, 2021, Example scripts for ensemble Kalman inversion of geophysical data
Tso, 2019, On the field estimation of moisture content using electrical geophysics: the impact of petrophysical model uncertainty, Water Resour. Res., 55, 7196, 10.1029/2019WR024964
Tso, 2020, Integrated hydrogeophysical modelling and data assimilation for geoelectrical leak detection, J. Contam. Hydrol., 234, 10.1016/j.jconhyd.2020.103679
Uhlemann, 2017, Four-dimensional imaging of moisture dynamics during landslide reactivation, J. geophys. Res., 122, 398, 10.1002/2016JF003983
Ullrich, 2008, Electrical resistivity tomography methods for archaeological prospection, Layers of Perception. Proceedings of the 35th International Conference on Computer Applications and Quantitative Methods in Archaeology (CAA)
Wainwright, 2014, Bayesian hierarchical approach and geophysical data sets for estimation of reactive facies over plume scales, Water Resour. Res., 50, 4564, 10.1002/2013WR013842
Ward, 2014, Distribution-based fuzzy clustering of electrical resistivity tomography images for interface detection, Geophys. J. Int., 197, 310, 10.1093/gji/ggu006
Whalley, 2017, Methods to estimate changes in soil water for phenotyping root activity in the field, Plant Soil, 415, 407, 10.1007/s11104-016-3161-1
Wu, 2019, Accurate imaging of hydraulic fractures using templated electrical resistivity tomography, Geothermics, 81, 74, 10.1016/j.geothermics.2019.04.004
Yang, 2014, Uncertainty quantification of CO2 saturation estimated from electrical resistance tomography data at the Cranfield site, Int. J. Greenhouse Gas Contl., 27, 59, 10.1016/j.ijggc.2014.05.006
Yeh, 2002, A geostatistically based inverse model for electrical resistivity surveys and its applications to vadose zone hydrology, Water Resour. Res., 38, WR001204, 10.1029/2001WR001204
Zahner, 2016, Image synthesis with graph cuts: a fast model proposal mechanism in probabilistic inversion, Geophys. J. Int., 204, 1179, 10.1093/gji/ggv517
Zheglova, 2013, 2-D reconstruction of boundaries with level set inversion of traveltimes, Geophys. J. Int., 192, 688, 10.1093/gji/ggs035
Zheglova, 2018, Multiple level-set joint inversion of traveltime and gravity data with application to ore delineation: a synthetic study, 83