Assimilation of hydrogeophysical data for the characterization of subsurface heterogeneity using Ensemble Kalman Filter (EnKF)
KANG Xueyuan1, SHI Xiaoqing1, DENG Yaping1, LIAO Kaihua2, WU Jichun1
1. Key Laboratory of Surficial Geochemistry, School of Earth Science and Engineering, Nanjing University, Nanjing 210023, China; 2. Key Laboratory of Watershed Geographic Sciences, Nanjing Institute of Geography and Limnology, Chinese Academy of Sciences, Nanjing 210008, China
Abstract:Characterization of spatial variability of hydrogeologic properties is the key to simulate and predict the fate and transport of contaminants in the subsurface. In this study, we present a sequential data assimilation framework to estimate the heterogeneous saturated hydraulic conductivity fields through the assimilation of Electrical Resistivity Tomography (ERT)-monitored data and groundwater flow/transport observation data. This framework is integrated Ensemble Kalman Filter (EnKF), groundwater flow/transport models and effective medium resistivity model. To test the performance of the framework, synthetic cases of contaminant transport are reconstructed. We compare the performance of the coupled and uncoupled methods. The factors to control the performance of coupled and uncoupled methods are also discussed in a number of different scenarios. Results showed that both methods can effectively estimate the spatial distribution of hydraulic conductivity via time-lapse ERT-monitored data. The coupled method performs better than the uncoupled one when the prior statistics are close to real field. Meanwhile, the uncoupled method is more robust when the prior statistics is biased. The accuracy of estimated heterogeneous parameter field could be improved when integrating of multiple type observations including ERT-monitored data and a few observations of groundwater flow/transport model (i. e., concentration). As the uncoupled method requires a small computational effort compared to the coupled one, it is suggested to use the uncoupled method as a preliminary inversion before refining the results with a fully coupled method. We conclude that integrating multiple types of observations is recommended to improve the ability to delineate subsurface heterogeneity.
康学远, 施小清, 邓亚平, 廖凯华, 吴吉春. 基于EnKF融合地球物理数据刻画含水层非均质性[J]. 水科学进展, 2018, 29(1): 40-49.
KANG Xueyuan, SHI Xiaoqing, DENG Yaping, LIAO Kaihua, WU Jichun. Assimilation of hydrogeophysical data for the characterization of subsurface heterogeneity using Ensemble Kalman Filter (EnKF). Advances in Water Science, 2018, 29(1): 40-49.
郭高轩, 吴吉春. 应用GPR获取水文地质参数研究初探[J]. 水文地质工程地质, 2005, 32(1):89-93. (GUO G X, WU J C. Approach to the estimation of hydrogeological parameters using GPR[J]. Hydrogeology and Engineering Geology, 2005, 32(1):89-93. (in Chinese))
[2]
HUISMAN J A, RINGS J, VRUGT J A, et al. Hydraulic properties of a model dike from coupled Bayesian and multi-criteria hydrogeophysical inversion[J]. Journal of Hydrology, 2010, 380(1):62-73.
[3]
KEMNA A, KULESSA B, VEREECKEN H. Imaging and characterisation of subsurface solute transport using electrical resistivity tomography (ERT) and equivalent transport models[J]. Journal of Hydrology, 2002, 267(3):125-146.
[4]
IRVING J, SINGHA K. Stochastic inversion of tracer test and electrical geophysical data to estimate hydraulic conductivities[J]. Water Resources Research, 2010, 46(11):2387-2392.
[5]
TRAN A P, VANCLOOSTER M, ZUPANSKI M, et al. Joint estimation of soil moisture profile and hydraulic parameters by ground-penetrating radar data assimilation with maximum likelihood ensemble filter[J]. Water Resources Research, 2014, 50(4):3131-3146.
[6]
EVENSEN G. The ensemble Kalman filter for combined state and parameter estimation[J]. IEEE Control Systems Magazine, 2009, 29(3):83-104.
[7]
史良胜, 张秋汝, 宋雪航, 等. 地下水位在非饱和水流数据同化中的应用[J]. 水科学进展, 2015, 26(3):404-412. (SHI L S, ZHANG Q R, SONG X H, et al. Application of groundwater level data to data assimilation for unsaturated flow[J]. Advances in Water Science, 2015, 26(3):404-412. (in Chinese))
[8]
南统超, 吴吉春. 集合卡尔曼滤波估计水文地质参数的局域化修正[J]. 水科学进展, 2010, 21(5):613-621. (NAN T C, WU J C. Localization corrections for the estimation of hydrogeological parameters using ensemble Kalman filter[J]. Advances in Water Science, 2010, 21(5):613-621. (in Chinese))
[9]
LI L, ZHOU H, HENDRICKS H J, et al. Groundwater flow inverse modeling in non-MultiGaussian media:performance assessment of the normal-score Ensemble Kalman Filter[J]. Hydrology and Earth System Sciences, 2012, 16(2):573-590.
[10]
CAMPORESE M, CASSIANI G, DEIANA R, et al. Coupled and uncoupled hydrogeophysical inversions using ensemble Kalman filter assimilation of ERT-monitored tracer test data[J]. Water Resources Research, 2015, 51(5):3277-3291.
[11]
JUHOJUNTTI N, KAMM J. Joint inversion of seismic refraction and resistivity data using layered models:applications to groundwater investigation[J]. Geophysics, 2015, 80(1):N43-N55.
[12]
SINGHA K, DAY-LEWIS F D, JOHNSON T, et al. Advances in interpretation of subsurface processes with time-lapse electrical imaging[J]. Hydrological Processes, 2015, 29(6):1549-1576.
[13]
ARCHIE G E. The electrical resistivity log as an aid in determining some reservoir characteristics[J]. Transactions of the American Institute of Mining and Metallurgical Engineers, 1942, 146:54-61.
[14]
BERG C. An effective medium algorithm for calculating water saturation at any frequency[J]. Geophysics, 2007, 72(2):E59-E67.
[15]
ZHENG C, WANG P P. MT3DMS:A modular three-dimensional multispecies transport model for simulation of advection, dispersion, and chemical reactions of contaminants in groundwater systems; documentation and user's guide[R]. Vicksburg:U S Army Engineer Research and Development Center Contract Report SERDP-99-1, 1999.
[16]
朱建友, 邓亚平, 施小清, 等. 高密度电阻率法探测DNAPLs污染的适宜性探讨[J]. 水文地质工程地质, 2017, 44(1):144-151. (ZHU J Y, DENG Y P, SHI X Q, et al. Discussion on the suitability of electrical resistivity tomography method to detect DNAPLs contamination[J]. Hydrogeology and Engineering Geology, 2017, 44(1):144-151. (in Chinese))
[17]
ATEKWANA E A, ATEKWANA E, LEGALL F D, et al. Biodegradation and mineral weathering controls on bulk electrical conductivity in a shallow hydrocarbon contaminated aquifer[J]. Journal of Contaminant Hydrology, 2005, 80(3/4):149-167.
[18]
SEN P N. Influence of temperature on electrical conductivity on shaly sands[J]. Geophysics, 1992, 57(1):89-96.
[19]
褚楠, 黄春林, 杜培军. 基于集合卡尔曼平滑算法的土壤水分同化[J]. 水科学进展, 2015, 26(2):243-249. (CHU N, HUANG C L, DU P J. Ensemble Kalman smoother for soil moisture data assimilation[J]. Advances in Water Science, 2015, 26(2):243-249. (in Chinese))
[20]
施小清, 吴吉春, 袁永生, 等. 含水介质各向异性对渗透系数空间变异性统计的影响[J]. 水科学进展, 2005, 16(5):679-684. (SHI X Q, WU J C, YUAN Y S, et al. Effect of the anisotropy in porous media on the spatial variability of the hydraulic conductivity[J]. Advances in Water Science, 2005, 16(5):679-684. (in Chinese))
[21]
DEUTSCH C V, JOURNEL A G. Geostatistical software library and user's guide[M]. New York:Oxford University Press, 1998:170-174.
[22]
SINGHA K, GORELICK S M. Saline tracer visualized with three-dimensional electrical resistivity tomography:field-scale spatial moment analysis[J]. Water Resources Research, 2005, 41(5):305-316.
[23]
YEH T J, SIMU。NEK J. Stochastic fusion of information for characterizing and monitoring the vadose zone[J]. Vadose Zone Journal, 2002, 1(2):207-221.
[24]
HERCKENRATH D, FIANDACA G, AUKEN E, et al. Sequential and joint hydrogeophysical inversion using a field-scale groundwater model with ERT and TDEM data[J]. Hydrology and Earth System Sciences, 2013, 17(10):4043-4060.
[25]
MICHAEL T C, ZHA Y, YEH T J, et al. The relative importance of head, flux, and prior information in hydraulic tomography analysis[J]. Water Resources Research, 2016, 52(1):3-20.