Groundwater contaminant source identification based on iterative local update ensemble smoother

Identification of the location and intensity of groundwater pollution source contributes to the effect of pollution remediation,and is called groundwater contaminant source identification.This is a kind of typical groundwater inverse problem,and the solution is usually ill-posed.Especially considering the spatial variability of hydraulic conductivity field,the identification process is more challenging.In this paper,the solution framework of groundwater contaminant source identification is composed with groundwater pollutant transport model(MT3DMS)and a data assimilation method(Iterative local update ensemble smoother,ILUES).In addition,Karhunen-Loève expansion technique is adopted as a PCA method to realize dimension reduction.In practical problems,the geostatistical method is usually used to characterize the hydraulic conductivity field,and only the contaminant source information is inversely calculated in the identification process.In this study,the identification of contaminant source information under Kriging K-field is compared with simultaneous identification of source information and K-field.The results indicate that it is necessary to carry out simultaneous identification under heterogeneous site,and ILUES has good performance in solving high-dimensional parameter inversion problems.

Groundwater contaminant source identification based on QS-ILUES

When groundwater pollution occurs,to come up with an efficient remediation plan,it is particularly important to collect information of contaminant source(location and source strength)and hydraulic conductivity field of the site accurately and quickly.However,the information can not be obtained by direct observation,and can only be derived from limited measurement data.Data assimilation of observations such as head and concentration is often used to estimate parameters of contaminant source.As for hydraulic conductivity field,especially for complex non-Gaussian field,it can be directly estimated by geostatistics method based on limited hard data,while the accuracy is often not high.Better estimation of hydraulic conductivity can be achieved by solving inverse groundwater problem.Therefore,in this study,the multi-point geostatistics method Quick Sampling(QS)is proposed and introduced for the first time and combined with the iterative local updating ensemble smoother(ILUES)to develop a new data assimilation framework QS-ILUES.It helps to solve the contaminant source parameters and non-Gaussian hydraulic conductivity field simultaneously by assimilating hydraulic head and pollutant concentration data.While the pilot points are utilized to reduce the dimension of hydraulic conductivity field,the influence of pilot points’layout and the ensemble size of ILUES algorithm on the inverse simulation results are further explored.