基于EnKF融合地球物理数据刻画含水层非均质性

含水层非均质性的刻画是模拟地下水中污染物运移的关键。以渗透系数为研究对象,构建了综合集合卡尔曼滤波方法、有效电阻率模型与地下水运移模型的同化框架,通过融合地球物理观测数据与污染物浓度观测数据来推估渗透系数的空间分布。基于理想算例,验证了该同化框架刻画含水层非均质渗透系数场的有效性,并针对不同初始参数信息与观测类型对比了耦合与非耦合水文地球物理方法的适用性。研究结果表明:基于集合卡尔曼滤波方法同化多种类型的观测数据,可有效地推估非均质参数空间分布。当初始信息较准确时,耦合方法的参数推估精度更高;初始信息存在偏差时,非耦合方法有更好的同化效果。由于非耦合方法计算成本较低且对初始信息缺失时适用性更强,在实际应用中可先基于非耦合方法初步估计参数,再利用耦合方法进一步提高参数推估精度。融合多种类型观测数据可有效提高参数推估效果。     

基于EnKF综合水头和浓度观测数据推估地下水流模型参数

地下水反应运移模型具有参数个数众多,观测数据类型多样的特点。为了探究不同类型观测数据在反应运移模拟数据同化中的数据价值,构建了三氯乙烯降解反应运移模型的理想算例,基于水头和浓度两种类型观测数据,采用集合卡尔曼滤波方法推估渗透系数和贮水系数的非均质空间分布,讨论了影响同化结果的因素。结果表明:与仅同化水头数据的结果相比,联合同化水头和浓度观测数据推估渗透系数场和贮水系数场时具有更高的精度,在观测数据拟合和模型预测方面也有更好的表现。与目前溶质运移模型、非饱和流模型等地下水模型中的研究结果相似,数据同化结果受样本数量,观测井的数量和位置的影响,合理优化布置监测井和选择样本数量可有效改善数据同化效果并提高计算效率。     

基于集合卡尔曼滤波的多相流模型参数估计——以室内二维砂箱中重质非水相污染物入渗为例

重质非水相有机污染物(DNAPL)泄漏到地下后,其运移与分布特征受渗透率非均质性影响显著。为刻画DNAPL污染源区结构特征,需进行参数估计以描述水文地质参数的非均质性。本研究构建了基于集合卡尔曼滤波方法(EnKF)与多相流运移模型的同化方案,通过融合DNAPL饱和度观测数据推估非均质介质渗透率空间分布。通过二维砂箱实际与理想算例,验证了同化方法的推估效果,并探讨了不同因素对同化的影响。研究结果表明:基于EnKF方法同化饱和度观测资料可有效地推估非均质渗透率场;参数推估精度随观测时空密度的增大而提高;观测点位置分布对同化效果有所影响,布置在污染集中区域的观测数据对于参数估计具有较高的数据价值。     

地下水污染强度及渗透系数场的反演识别研究

在制定地下水污染修复方案时,污染源参数和渗透系数场是最重要的地下水数值模型参数,但前人研究多集中于单一类型参数的识别。文章中采用地下水污染物运移模型（MT3DMS）和数据同化方法（迭代局部更新集合平滑器,ILUES）构成地下水污染源识别的求解框架,并利用Karhunen-Loève展开技术实现渗透系数场的参数降维,最后通过同化水头与浓度数据实现地下水污染源强和渗透系数场的联合反演。结果表明：（1）ILUES算法能精确识别污染源参数和渗透系数场,并且具有很高的普适性;（2）精确表征渗透系数在空间上呈现出的非均质性,是预测污染物迁移路径、反演污染强度的关键;（3）ILUES算法参数影响着反演效果,综合考虑计算效率和计算精度等,可以得到算例的最佳样本集合大小（Ne=4 000）和ILUES算法最佳参数组合（局部临近样本集合占比α=0.4,相对权重b=4）。但在实际工程案例中,如果对精度的要求不是过高,经验组合（α=0.1,b=1）更值得推荐。研究结果对于区域地下水资源调查、评价和管理等工作具有较强的实践意义,并可为后期地下水污染预测及地下水监测井网优化提供技术支撑。     

基于局域化EnKF同时同化地下水水流和溶质运移数据的渗透系数场识别研究

局域化改进集合卡尔曼滤波（EnKF）可以克服EnKF方法在使用小集合时,对参数识别精度较低的缺陷,其能同化 地下水位观测数据有效识别渗透系数场。实际工作中,溶质运移数据也较容易获得。崔凯鹏（2013）尝试增加溶质运移 数据以改进只同化水流数据对渗透系数的估计结果,但是精度提高有限。本文在其基础上修改模型,进一步增加溶质注 入井,探究同时同化水头和溶质运移数据,对渗透系数场识别效果,之后对比了局域化EnKF与非局域化EnKF参数识别结 果,并分析了溶质影响范围与参数识别的关系。结果表明：同时同化溶质运移和水头资料,比同化单一种类观测数据识别 的渗透系数精度更高;相同实现数目下,局域化EnKF比EnKF对渗透系数场的估计结果与真实场更为接近;仅考虑溶质影 响范围内的渗透系数,同化水头数据在最后时刻参数识别结果好于同化溶质运移数据参数识别结果,但差别不大。     

基于遗传算法-支持向量机的滑坡渗透系数反演

求解库岸边坡岩土体的渗透系数是研究滑坡渗流场及多场演化的基础,一般通过原位试验和室内试验求得,但试验成本较高且试验位置具有一定的随机性。本文以三峡库区马家沟滑坡为例,提出一种利用地下水位动态观测资料反演滑坡岩土层渗透系数的方法。具体步骤为:（1）依据滑坡的勘察资料和水位观测数据,构建滑坡数值模型;（2）利用SPSS生成不同渗透系数正交试验组合,并将渗透系数代入数值模型中计算监测井的水位,得到不同渗透系数及其对应的模拟水位数据;（3）应用遗传算法优化的支持向量机构建坡体模拟水位与渗透系数的非线性映射关系,再通过代入实际动态监测水位值求得滑坡岩土层的渗透系数;（4）将求得的渗透系数代入数值模型,用计算的模拟水位与实际观测水位进行对比验证。研究结果表明:遗传算法优化的支持向量机具有良好的学习预测效果,能准确预测渗透系数与水位的关系。该反演方法具有高效、准确的优点,反演结果的精度满足实际应用需要。     

二维平面非恒定流数学模型的遥感水位数据同化

为了在水流计算中定量化利用遥感水位数据,基于偏微分方程最优化控制理论建立了变分模型来融合二维平面非恒定流的数学模型和数据。根据遥感数据空间信息密集的特点,提出遥感水位数据同化的新算法。采用人工合成数据考察了面域遥感数据对于定常参数和时变两类参数的反演效果。试验结果表明,遥感数据提供的空间分布式信息有利于空间分布式参数的反演识别,而且通过引入考虑水面空间变化信息的附加项,可以改善观测信息的同化,更好地辨识时变参数(流量过程)。以Moselle河的RADARSAT卫星遥感水位数据检验了模型的实用性。     

基于不同地质统计方法的渗透系数场对污染物运移的影响

渗透系数场的空间变异性是影响污染物运移结果的决定因素,而地质统计方法是解决渗透系数空间变异性的主要技术手段。本文利用野外场地实测数据,采用普通克里格法和指示克里格法、顺序高斯模拟法和顺序指示模拟法四种地质统计方法,插值估测和模拟再现随机渗透系数场,进而对比研究四种渗透系数场对大尺度污染物运移的影响。研究结果表明,污染羽的质心位置（一阶矩）主要由渗透系数的平均值来决定;污染羽在空间上的展布范围（二阶矩）主要受渗透系数空间变异方差的影响;条件模拟克服了估计法的平滑效果,较好地再现真实曲线的波动性,渗透系数（ lnK）估计方差与污染羽空间二阶矩随着条件模拟次数的增加而减小,并且顺序指示模拟程度更加明显。     

基于集合卡尔曼平滑算法的土壤水分同化

为研究观测资料稀少情况下土壤质地及有机质对土壤水分同化的影响,发展了集合卡尔曼平滑(Ensemble Kalman Smooth, EnKS)的土壤水分同化方案。利用黑河上游阿柔冻融观测站2008年6月1日至10月29日的观测数据,使用EnKS算法将表层土壤水分观测数据同化到简单生物圈模型(Simple Biosphere Model 2, SiB2)中,分析不同方案对土壤水分估计的影响,并与集合卡尔曼滤波算法(EnKF)的结果进行比较。研究结果表明,土壤质地和有机质对表层土壤水分模拟结果影响最大而对深层的影响相对较小;利用EnKF和EnKS算法同化表层土壤水分观测数据,均能够显著提高表层和根区土壤水分估计的精度,EnKS算法的精度略高于EnKF且所受土壤质地和有机质的影响小于EnKF;当观测数据稀少时,EnKS算法仍然可以得到较高精度的土壤水分估计。     

土壤水分同化系统的敏感性试验研究

利用1998年7月6日至8月9日青藏高原GAME-Tibet试验区MS3608站点的4cm、20cm和100cm的土壤水分观测数据同化SiB2模型输出的表层、根区和深层土壤水分,探讨了一个基于集合卡尔曼滤波和简单生物圈模型的单点土壤水分同化方案。分析和评价了集合大小、同化周期、模型误差、背景场误差以及观测误差对同化系统性能的影响。结果表明:①增加集合数目可以减小土壤水分同化系统的误差,但同时又降低了运行效率;②对于集合卡尔曼滤波,初始场的估计是否准确对同化系统性能影响不大;③模型误差和观测误差的准确估计可以提高土壤水分的估计精度;④利用数据同化的方法对土壤水分的估计有显著提高。     

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.     

Pattern Recognition in a Bimodal Aquifer Using the Normal-Score Ensemble Kalman Filter

The ensemble Kalman filter (EnKF) is now widely used in diverse disciplines to estimate model parameters and update model states by integrating observed data. The EnKF is known to perform optimally only for multi-Gaussian distributed states and parameters. A new approach, the normal-score EnKF (NS-EnKF), has been recently proposed to handle complex aquifers with non-Gaussian distributed parameters. In this work, we aim at investigating the capacity of the NS-EnKF to identify patterns in the spatial distribution of the model parameters (hydraulic conductivities) by assimilating dynamic observations in the absence of direct measurements of the parameters themselves. In some situations, hydraulic conductivity measurements (hard data) may not be available, which requires the estimation of conductivities from indirect observations, such as piezometric heads. We show how the NS-EnKF is capable of retrieving the bimodal nature of a synthetic aquifer solely from piezometric head data. By comparison with a more standard implementation of the EnKF, the NS-EnKF gives better results with regard to histogram preservation, uncertainty assessment, and transport predictions.     

A Simple Approach to Account for Radial Flow and Boundary Conditions When Kriging Hydraulic Head Fields for Confined Aquifers

The estimation and mapping of realistic hydraulic head fields, hence of flow paths, is a major goal of many hydrogeological studies. The most widely used method to obtain reliable head fields is the inverse approach. This approach relies on the numerical approximation of the flow equation and requires specifying boundary conditions and the transmissivity of each grid element. Boundary conditions are often unknown or poorly known, yet they impose a strong signature on the head fields obtained by inverse analysis. A simpler alternative to the inverse approach is the direct kriging of the head field using the measurements obtained at observation wells. The kriging must be modified to incorporate the available information. Use of the dual kriging formalism enables simultaneously estimating the head field, the aquifer mean transmissivity, and the regional hydraulic gradient from head data in steady or transient state conditions. In transient state conditions, an estimate of the storage coefficient can be obtained. We test the approach on simple analytical cases, on synthetic cases with solutions obtained numerically using a finite element flow simulator, and on a real aquifer. For homogeneous aquifers, infinite or bounded, the kriging estimate retrieves the exact solution of the head field, the exact hydrogeological parameters and the flow net. With heterogeneous aquifers, kriging accurately estimates the head field with prediction errors of the same magnitude as typical head measurement errors. The transmissivities are also accurately estimated by kriging. Moreover, if inversion is required, the kriged head along boundaries can be used as realistic boundary conditions for flow simulation.     

Review of self-potential methods in hydrogeophysics

The self-potential (SP) method is a passive geophysical method based on the natural occurrence of electrical fields on the Earth's surface. Combined with other geophysical methods, SP surveys are especially useful for localizing and quantifying groundwater flows and pollutant plume spreading, and estimating pertinent hydraulic properties of aquifers (water table, hydraulic conductivity). Laboratory experiments have shown that the involved coupling coefficients mainly depend on the fluid chemistry, conductivity and pH, and on the soil or rock properties. The interpretation of SP observations can be done qualitatively, for instance, by correlation of SP gradients with water fluxes (through electrokinetics) or salt fluxes (through electro-diffusion). In recent years, the interpretation has been improved with the help of modelling or/and inversion of the Poisson equation and endeavours to estimate hydraulic parameters by means of the intensity of electric current sources caused by underground flows.     

Empirical Orthogonal Functions Analysis Applied to the Inverse Problem in Hydrogeology: Evaluation of Uncertainty and Simulation of New Solutions

To fulfil the need to generate more realistic solutions, stochastic inverse simulations in hydrogeology are now constrained on both piezometric head and hydraulic conductivity data. These inverse techniques, often based on geostatistics, allow modifications of an initial solution conditioned only on hydraulic conductivity data to arrive at a final solution that also matches observed heads. By repeating the process as many times as necessary with different initial solutions, one generates an ensemble of final solutions thereby addressing the uncertainty of the inverse problem. This requires a method able to handle the whole ensemble and to work on its relevant characteristics. From this standpoint, the analysis by Empirical Orthogonal Functions (EOF) appears promising. The method builds an orthogonal decomposition of the covariance matrix, calculated over the whole set of solutions, and the areas in space where the first functions have a greater influence corresponding to locations of maximum uncertainty in the solutions. These locations depend both on the hydraulic characteristics of the flow problem and on the spatial distribution of available data. The EOF analysis is used on a synthetic problem that mimics a possible behavior of the Culebra aquifer of the Waste Isolation Pilot Plant (WIPP, New Mexico). The method also allows new solutions to be generated at lower computational cost by a random composition of the functions obtained by the EOF analysis. These new solutions keep the main characteristics of the initial ensemble and because they can be conditioned, they return very good results when they are used to solve the direct problem.     

Oscillatory Pumping Test to Estimate Aquifer Hydraulic Parameters in a Bayesian Geostatistical Framework

Comprehensive information about the spatial distribution of the subsurface hydraulic properties is crucial to model groundwater flow, to predict solute transport in aquifers and to design remediation actions. In this work, a Bayesian Geostatistical approach, as implemented in bgaPEST, was adopted to estimate the hydraulic properties of a well field located at the Campus of Science and Technology of the University of Parma (Northern Italy), in a contest of a highly parameterized inversion. Head data, collected by means of multi frequency oscillatory pumping tests, were used to both estimate the hydraulic parameters and validate the results. The groundwater flow processes were modelled by means of MODFLOW 2005 and an adjoint-state formulation of the same software was used to efficiently calculate the sensitivity matrix, required by the inverse procedure. The Bayesian Geostatistical approach estimated the hydraulic conductivity and specific storage fields, handling a large number of parameters. The results of the inversion are consistent with the alluvial nature of the investigated aquifer and the preliminary traditional pumping tests carried out at the site.     

A cost-effective device and methodology to compute aquifer transmissivity and piezometry from free-flowing artesian wells

Artesian aquifers offer interesting opportunities for water supply by providing a low-vulnerability groundwater resource that is easily abstracted without any installation of pumps or power supply costs. However, hydraulic tests are challenging to perform, notably where the piezometric head is above ground level with free-flowing wells not equipped with valves and open for years. This paper describes a low-cost, easy to reproduce and adaptable device, the free-flowing artesian well device (FFAWD), which is mainly designed with a set of PVC tubes equipped with a pressure probe and a valve. This device is used to perform hydraulic tests on free-flowing artesian wells, to measure the piezometric head of the aquifer and to compute its transmissivity. The practical use of the FFAWD is described and a method is proposed to compute the piezometric head and the transmissivity of the aquifer from this data set (free-flowing well discharge and pressure increase measurements) with any adapted analytical solution, using the Houpeurt-Pouchan method. Artefacts such as post-production effects, surge effects, and the impact of a leaky well are identified to avoid any misinterpretation. The FFAWD was applied to the volcano-sedimentary artesian plain of Pasuruan (Indonesia). The advantages and limitations of using the device, along with the interpretation methodology, are also discussed.