The facies distribution of a reservoir is one of the biggest concerns for geologists,geophysicists,reservoir modelers,and reservoir engineers due to its high importance in the setting of any reliable decisionmaking/op...The facies distribution of a reservoir is one of the biggest concerns for geologists,geophysicists,reservoir modelers,and reservoir engineers due to its high importance in the setting of any reliable decisionmaking/optimization of field development planning.The approach for parameterizing the facies distribution as a random variable comes naturally through using the probability fields.Since the prior probability fields of facies come either from a seismic inversion or from other sources of geologic information,they are not conditioned to the data observed from the cores extracted from the wells.This paper presents a regularized element-free Galerkin(R-EFG)method for conditioning facies probability fields to facies observation.The conditioned probability fields respect all the conditions of the probability theory(i.e.all the values are between 0 and 1,and the sum of all fields is a uniform field of 1).This property achieves by an optimization procedure under equality and inequality constraints with the gradient projection method.The conditioned probability fields are further used as the input in the adaptive pluri-Gaussian simulation(APS)methodology and coupled with the ensemble smoother with multiple data assimilation(ES-MDA)for estimation and uncertainty quantification of the facies distribution.The history-matching of the facies models shows a good estimation and uncertainty quantification of facies distribution,a good data match and prediction capabilities.展开更多
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 gro...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.展开更多
针对未知的污染场地,为了准确估计污染物运移模型的参数,提出一种基于多重数据同化集合平滑器(ensemble smoother with multiple data assimilation,ES-MDA)算法的地下水模型参数反演方法,通过融合由高密度电阻率(electrical resistance...针对未知的污染场地,为了准确估计污染物运移模型的参数,提出一种基于多重数据同化集合平滑器(ensemble smoother with multiple data assimilation,ES-MDA)算法的地下水模型参数反演方法,通过融合由高密度电阻率(electrical resistance tomography,ERT)法采集的ERT观测数据,实现对污染源源强和渗透系数场的联合反演。以此为基础设计3组数值算例,比较不同类型观测数据对反演精度的影响。研究结果表明:融合ERT数据的ES-MDA算法对模型参数的反演精度更高,并且将ERT数据和传统的质量浓度与水头观测数据相结合,能进一步优化反演结果。展开更多
In this paper, a new state-parameter estimation approach is presented based on the dual ensemble Kalman smoother(DEn KS) and simple biosphere model(Si B2) to sequentially estimate both the soil properties and soil moi...In this paper, a new state-parameter estimation approach is presented based on the dual ensemble Kalman smoother(DEn KS) and simple biosphere model(Si B2) to sequentially estimate both the soil properties and soil moisture profile by assimilating surface soil moisture observations. The Arou observation station, located in the upper reaches of the Heihe River in northwestern China, was selected to test the proposed method. Three numeric experiments were designed and performed to analyze the influence of uncertainties in model parameters, atmospheric forcing, and the model's physical mechanics on soil moisture estimates. Several assimilation schemes based on the ensemble Kalman filter(En KF), ensemble Kalman smoother(En KS), and dual En KF(DEn KF) were also compared in this study. The results demonstrate that soil moisture and soil properties can be simultaneously estimated by state-parameter estimation methods, which can provide more accurate estimation of soil moisture than traditional filter methods such as En KF and En KS. The estimation accuracy of the model parameters decreased with increasing error sources. DEn KS outperformed DEn KF in estimating soil moisture in most cases, especially where few observations were available. This study demonstrates that the DEn KS approach is a useful and practical way to improve soil moisture estimation.展开更多
构建了基于通用陆面模型(CoLM,Common Land Model)、微波辐射传输模型L-MEB(Lband Microwave Emission of the Biosphere)和集合平滑算法(EnKS,Ensemble Kalman Smoother)的土壤水分数据同化框架,用于联合同化MODIS地表温度和机载L波段...构建了基于通用陆面模型(CoLM,Common Land Model)、微波辐射传输模型L-MEB(Lband Microwave Emission of the Biosphere)和集合平滑算法(EnKS,Ensemble Kalman Smoother)的土壤水分数据同化框架,用于联合同化MODIS地表温度和机载L波段被动微波亮温数据。以2012年HiWATER试验期间中游大满超级站为实验站点,分析了3种LAI数据产品对土壤温度模拟结果的影响,进而分析了联合同化地表温度和微波亮度温度对土壤水分估计结果的影响。研究结果表明:3种LAI数据对土壤温度模拟结果的影响显著,MODIS LAI产品在该研究区显著低估,导致土壤温度模拟结果高估4~6K;同化亮度温度、同化地表温度以及联合同化两者均可以改进土壤水分的估计精度,联合同化地表温度和亮度温度对于土壤水分的改进最为显著,土壤水分同化结果的RMSE减少31%~53%。展开更多
文摘The facies distribution of a reservoir is one of the biggest concerns for geologists,geophysicists,reservoir modelers,and reservoir engineers due to its high importance in the setting of any reliable decisionmaking/optimization of field development planning.The approach for parameterizing the facies distribution as a random variable comes naturally through using the probability fields.Since the prior probability fields of facies come either from a seismic inversion or from other sources of geologic information,they are not conditioned to the data observed from the cores extracted from the wells.This paper presents a regularized element-free Galerkin(R-EFG)method for conditioning facies probability fields to facies observation.The conditioned probability fields respect all the conditions of the probability theory(i.e.all the values are between 0 and 1,and the sum of all fields is a uniform field of 1).This property achieves by an optimization procedure under equality and inequality constraints with the gradient projection method.The conditioned probability fields are further used as the input in the adaptive pluri-Gaussian simulation(APS)methodology and coupled with the ensemble smoother with multiple data assimilation(ES-MDA)for estimation and uncertainty quantification of the facies distribution.The history-matching of the facies models shows a good estimation and uncertainty quantification of facies distribution,a good data match and prediction capabilities.
基金supported by the Fundamental Research Funds for the Central Universities(No.22120190013)National Natural Science Foundation of China(No.41807187)
文摘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.
文摘针对未知的污染场地,为了准确估计污染物运移模型的参数,提出一种基于多重数据同化集合平滑器(ensemble smoother with multiple data assimilation,ES-MDA)算法的地下水模型参数反演方法,通过融合由高密度电阻率(electrical resistance tomography,ERT)法采集的ERT观测数据,实现对污染源源强和渗透系数场的联合反演。以此为基础设计3组数值算例,比较不同类型观测数据对反演精度的影响。研究结果表明:融合ERT数据的ES-MDA算法对模型参数的反演精度更高,并且将ERT数据和传统的质量浓度与水头观测数据相结合,能进一步优化反演结果。
基金supported by the Natural National Science Foundation of China(Grant Nos.91325106&41271358)the Hundred Talent Program of the Chinese Academy of Sciences(Grant No.29Y127D01)+1 种基金the Cross-disciplinary Collaborative Teams Program for ScienceTechnology and Innovation of the Chinese Academy of Sciences
文摘In this paper, a new state-parameter estimation approach is presented based on the dual ensemble Kalman smoother(DEn KS) and simple biosphere model(Si B2) to sequentially estimate both the soil properties and soil moisture profile by assimilating surface soil moisture observations. The Arou observation station, located in the upper reaches of the Heihe River in northwestern China, was selected to test the proposed method. Three numeric experiments were designed and performed to analyze the influence of uncertainties in model parameters, atmospheric forcing, and the model's physical mechanics on soil moisture estimates. Several assimilation schemes based on the ensemble Kalman filter(En KF), ensemble Kalman smoother(En KS), and dual En KF(DEn KF) were also compared in this study. The results demonstrate that soil moisture and soil properties can be simultaneously estimated by state-parameter estimation methods, which can provide more accurate estimation of soil moisture than traditional filter methods such as En KF and En KS. The estimation accuracy of the model parameters decreased with increasing error sources. DEn KS outperformed DEn KF in estimating soil moisture in most cases, especially where few observations were available. This study demonstrates that the DEn KS approach is a useful and practical way to improve soil moisture estimation.
文摘构建了基于通用陆面模型(CoLM,Common Land Model)、微波辐射传输模型L-MEB(Lband Microwave Emission of the Biosphere)和集合平滑算法(EnKS,Ensemble Kalman Smoother)的土壤水分数据同化框架,用于联合同化MODIS地表温度和机载L波段被动微波亮温数据。以2012年HiWATER试验期间中游大满超级站为实验站点,分析了3种LAI数据产品对土壤温度模拟结果的影响,进而分析了联合同化地表温度和微波亮度温度对土壤水分估计结果的影响。研究结果表明:3种LAI数据对土壤温度模拟结果的影响显著,MODIS LAI产品在该研究区显著低估,导致土壤温度模拟结果高估4~6K;同化亮度温度、同化地表温度以及联合同化两者均可以改进土壤水分的估计精度,联合同化地表温度和亮度温度对于土壤水分的改进最为显著,土壤水分同化结果的RMSE减少31%~53%。