[1]张双圣,强静,刘汉湖,等.基于拉丁超立方抽样的改进型DRAM算法求解地下水污染反问题[J].郑州大学学报(工学版),2020,41(03):72-78.[doi:10.13705/j.issn.1671-6833.2019.02.016]
 Zhang Shuangsheng,strong static,Liu Han Lake,et al.Improved MCMC method based on latin hypercube sampling for inverse problems of underground water pollution[J].Journal of Zhengzhou University (Engineering Science),2020,41(03):72-78.[doi:10.13705/j.issn.1671-6833.2019.02.016]
点击复制

基于拉丁超立方抽样的改进型DRAM算法求解地下水污染反问题()
分享到:

《郑州大学学报(工学版)》[ISSN:1671-6833/CN:41-1339/T]

卷:
41卷
期数:
2020年03期
页码:
72-78
栏目:
出版日期:
2020-07-29

文章信息/Info

Title:
Improved MCMC method based on latin hypercube sampling for inverse problems of underground water pollution
作者:
张双圣强静刘汉湖刘喜坤孙韶华
1. 中国矿业大学环境与测绘学院;2. 徐州市城区水资源管理处;  3. 中国矿业大学数学学院;4. 山东省城市供排水水质监测中心 
Author(s):
Zhang Shuangsheng12strong static 3Liu Han Lake 1Liu Xikun 2Sun Shaohua 4
1. School of Environment and Surveying, China University of Mining and Technology; 2. Xuzhou Urban Water Resources Management Office;  3. School of Mathematics, China University of Mining and Technology; 4. Shandong Provincial Water Supply and Drainage Water Quality Monitoring Center 
关键词:
二维水质模型贝叶斯-马尔科夫链蒙特卡洛法拉丁超立方抽样延迟拒绝自适应Metropolis算法污染源识别
Keywords:
2D water quality modelBayesian-Markov chain Monte Carlo methodLatin hypercube samplingDelay rejection adaptive Metropolis algorithmpollution source identification
DOI:
10.13705/j.issn.1671-6833.2019.02.016
文献标志码:
A
摘要:
针对运用贝叶斯统计方法求解地下水污染反问题时,经典MCMC算法求解结果受样本初始点影响的问题,提出了一种基于拉丁超立方抽样方法的改进型MCMC算法.将贝叶斯统计方法与二维水质对流-扩散方程相耦合,建立地下水污染源识别模型.构建一个污染物在地下水含水层中瞬时排放的算例,分别进行经典MCMC算法与改进型MCMC算法对污染源信息(污染源强度、排放位置(x,y)和排放时刻)的反求.算例研究表明,经典MCMC算法受样本初始点影响,容易出现反演结果局部最优,或者反演结果难以收敛的问题;改进型MCMC算法可实现反演结果的全局最优,而且反演结果不受样本初始点的影响,反演结果与真值基本一致,准确性与稳定性显著提高.
Abstract:
Aiming at the problem of the result affected by samples’ initial values with classical MCMC method, when the inverse problems of underground water pollution were solved by Bayesian statistical methods, an improved MCMC method based on latin hypercube sampling was presented. An underground water pollution source identification model was built by coupling Bayesian statistical methods to two-dimensional water quality convection-diffusion equation. An example of a pollutant in the underground aquifer discharged instantly was put forward, and the pollution source information including source’’s position, intensity and discharging time was solved by classical MCMC method and improved MCMC method separately. The example showed that the inversion results affected by initial values with classical MCMC method were locally optimal or difficult to convergence; on the contrary, the improved MCMC method could achieve a global optimization. The computed values solved by the improved MCMC method were basically the same as the truth values. The accuracy and stability were improved significantly.
更新日期/Last Update: 2020-07-29