Evaluation of the benefits of improved permeability estimation on high-resolution characterization of DNAPL distribution in aquifers with low-permeability lenses

• Published in: Journal of hydrology (Amsterdam) 2021-12, Vol.603, p.126955, Article 126955
• Main Authors: Guo, Qiongze, Shi, Xiaoqing, Kang, Xueyuan, Hao, Shefeng, Liu, Ling, Wu, Jichun
• Format: Article
• Language: English
• Subjects: Clay content; Combined ESMDA-DS and ESMDA method; DNAPL source zone characterization; Electrical resistivity tomography
• ISSN: 0022-1694; 1879-2707
• DOI: 10.1016/j.jhydrol.2021.126955

•The improved permeability characterization using the ESMDA-DS method can improve DNAPL mapping.•Integrating multiple data (especially geophysical data) helps characterize the DNAPL source zone architecture.•Iterative inversion framework is used to incorporate the clay content information for DNAPL source zone characterization. High-resolution characterization of dense non-aqueous phase liquid (DNAPL) source zone architecture (SZA) is crucial for developing remediation strategies. Inverse modeling approaches are widely used to characterize DNAPL SZA. However, unknown parameters (permeability and DNAPL saturation) are usually assumed as Gaussian distribution in traditional geostatistical inverse methods (e.g., ensemble smoother with multiple data assimilation, ESMDA). This assumption may not be true in channels of fluvial deposits. The presence of clay lenses will bring significant challenges in reconstructing the DNAPL SZA using geostatistical inverse methods based on the Gaussian assumption. Moreover, the DNAPL characterization also suffers from a limited number of borehole data in the field. In this study, we use ESMDA's variant, ensemble smoother with multiple data assimilation-direct sampling (ESMDA-DS), to improve the estimation of non-Gaussian permeability field and evaluate the benefits of improved permeability estimation on the reconstruction of DNAPL source zone. The inversion technique is based on the iterative inversion framework to incorporate the clay content information into the DNAPL characterization. The performance of our proposed framework is examined and discussed through a numerical synthetic case with challenging pool-dominated SZA. To overcome data scarcity, we integrated hydrogeological and geophysical (electrical resistance tomography) datasets. Results show that the permeability is characterized well using ESMDA-DS method compared with ESMDA method. With the improved permeability estimation, the DNAPL source zone can be characterized in higher resolution. Compared with ESMDA method, the root mean square error (RMSE) of permeability is reduced by 38.5% using the ESMDA-DS method. Then the estimation RMSE of DNAPL saturation is reduced by 24.9% by improved permeability estimation. We also demonstrated the significant benefits of improved permeability estimation on DNAPL depletion behavior and plume evolution. Numerical simulations demonstrate with the improved permeability estimation, the DNAPL depletion prediction is highly improved,