×

使用贝叶斯框架进行地下污染物运移预测。 (英语) Zbl 1472.62038号

摘要:在监测地下含水层污染时,我们希望利用具有专业知识和有限数据的地下流体流动模型预测数量——污染物浓度的分数流量曲线。这里考虑的是贝叶斯方法,与仿真研究相关的复杂性是一个持续的实际挑战。在两阶段马尔可夫链蒙特卡罗(MCMC)方法中,我们将渗透率场的Karhunen-Loève展开与GPU计算结合使用。通过运行几个MCMC链可以进一步降低计算成本。我们比较收敛标准来量化预测的不确定性。我们的贡献是双重的:我们首先提出了一个多元潜在标度缩减因子(MPSRF)数据的拟合程序,该程序允许我们估计收敛的迭代次数。然后,我们对分数流量曲线的集合进行了仔细的分析,表明对于当前的问题,通过MPSRF分析进行收敛所需的迭代次数过多。因此,对于实际应用,我们的结果表明,对感兴趣量的后验分布的分析为终止MCMC模拟以量化不确定性提供了可靠的标准。

MSC公司:

2015年1月62日 贝叶斯推断
86A05型 水文学、水文学、海洋学
86A32型 地理统计学

软件:

博阿
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] 希格顿,D。;Lee,H。;Bi,Z.,《利用粗尺度和细尺度信息表征反问题中不确定性的贝叶斯方法》,IEEE Trans。信号处理。,50, 389-399 (2002) ·Zbl 1369.62047号
[2] Lee,H。;希格顿,D。;Bi,Z。;费雷拉,M。;West,M.,多孔介质中流体流动模拟中高维参数的马尔可夫随机场模型,技术报告(2002),Technometrics
[3] 马,X。;Al-Harbi,M。;Datta-Gupta,A。;Efendiev,Y.,《历史匹配地质模型中不确定性量化的有效两阶段抽样方法》,SPE J.,77-87(2008)
[4] 希格顿,D。;Reese,C.S。;莫尔顿,J.D。;Vrugt,J.A。;Fox,C.,计算密集型正向模型的后验探索,(Brooks,S.;Gelman,A.;Jones,G.;Meng,X.-L.,《马尔可夫链蒙特卡罗手册》(2011),Chapman&Hall/CRC),401-418·Zbl 1229.65016号
[5] 金廷,V。;佩雷拉,F。;Rahunanthan,A.,《地下水流的多阶段贝叶斯预测框架》,《国际不确定性杂志》。数量。,499-522年3月6日(2013年)·Zbl 07694415号
[6] 金廷,V。;佩雷拉,F。;Rahunanthan,A.,多孔介质中流量预测的多马尔可夫链蒙特卡罗方法,Procedia Compute。科学。,9, 707-716 (2012)
[7] 金廷,V。;佩雷拉,F。;Rahunanthan,A.,预测多孔介质流动的预取技术,计算。地质科学。,1861-675年5月18日(2014年)·Zbl 1392.76081号
[8] Brockwell,A.,《预取的并行马尔可夫链蒙特卡罗模拟》,J.Comp。图表。统计。,15, 246-261 (2006)
[9] 金廷,V。;佩雷拉,F。;Rahunanthan,A.,地下水流的多物理马尔可夫链蒙特卡罗方法,数学。计算。同时。,118, 224-238 (2015) ·Zbl 07313418号
[10] 马蒙,A。;佩雷拉,F。;Rahunanthan,A.,地下水流问题MCMC方法的收敛分析,(Gervasi,O.;Murgante,B.;Misra,S.;Stankova,E.;Torre,C.M.;Rocha,A.M.A.;Taniar,D.;Apduhan,B.O.;Tarantino,E.,Ryu,Y.,《计算科学及其应用——ICCSA 2018(2018)》,Springer International Publishing,305-317
[11] 金廷,V。;佩雷拉,F。;Rahunanthan,A.,《快速量化油藏渗透率和孔隙度的不确定性以实现预测模拟》,数学。计算。同时。,99, 139-152 (2014) ·Zbl 07312586号
[12] 布鲁克斯,S。;Gelman,A.,《监测迭代模拟收敛性的通用方法》,J.Compute。图表。Stat.,7434-455(1998)
[13] 布鲁克斯,S。;Roberts,G.,马尔可夫链蒙特卡罗算法的收敛性评估,统计计算。,8,319-335(1998年)
[14] Oliver,D.S。;库尼亚,L.B。;Reynolds,A.C.,《根据压力数据调节渗透率场的马尔可夫链蒙特卡罗方法》,数学。地质。,29, 1, 61-91 (1997)
[15] Loève,M.,概率论(1997),施普林格:柏林施普林格出版社
[16] Wong,E.,《信息和动力系统中的随机过程》(1971),纽约麦格劳-希尔出版社·Zbl 0245.60001号
[17] 伊芬迪耶夫,Y。;Datta-Gupta,A。;金廷,V。;马,X。;Mallick,B.,《动态数据集成的有效两阶段马尔可夫链蒙特卡罗方法》,《水资源》。决议,41(12)(2005)
[18] Douglas,C。;伊芬迪耶夫,Y。;尤因,R。;金廷,V。;Lazarov,R.,《随机环境中的动态数据驱动模拟》,《计算》,77(4),321-333(2006)·Zbl 1122.65006号
[19] 伊芬迪耶夫,Y。;Hou,T。;Luo,W.,使用粗尺度模型预处理马尔可夫链蒙特卡罗模拟,SIAM J.Sci。计算。,28(2), 776-803 (2006) ·Zbl 1111.65003号
[20] 金廷,V。;佩雷拉,F。;普雷绍,M。;Wo,S.,两阶段马尔可夫链蒙特卡罗方法在使用代理流模型表征裂缝性储层中的应用,计算。地质科学。,15 (4), 691-707 (2011) ·兹比尔1237.76190
[21] 贝尔·J。;Cheng,A.H.-D.,《地下水流动和污染物运移建模》(2010年),施普林格出版社·Zbl 1195.76002号
[22] 陈,Z。;Huan,G。;Ma,Y.,多孔介质中多相流的计算方法(2006),SIAM:SIAM Philadelphia,PA·Zbl 1092.76001号
[23] 佩雷拉,F。;Rahunanthan,A.,GPU上两相流的数值模拟,第九届计算科学高性能计算国际会议(VECPAR’10),加州伯克利(2010)
[24] 佩雷拉,F。;Rahunanthan,A.,三维两相流近似的半离散中心格式,数学。计算。同时。,81(10), 2296-2306 (2011) ·Zbl 1419.76459号
[25] Christen,J。;Fox,C.,《使用近似的马尔可夫链蒙特卡罗》,J.Compute。图表。统计学,14(4),795-810(2005)
[26] 罗伯特·C。;Casella,G.,Monte Carlo统计方法(2004),纽约斯普林格·Zbl 1096.62003年
[27] Durlowsky,L.,非均质多孔介质等效网格块渗透率张量的数值计算,水资源。决议,27(5),699-708(1991)
[28] M.K.Cowles,B.C.,《马尔可夫链蒙特卡罗收敛诊断:比较综述》,《美国统计协会期刊》,91,883-904(1996)·Zbl 0869.62066号
[29] Mengersen,K。;罗伯特,C.P。;Guihenneuc-Jouyaux,C.,MCMC收敛诊断:综述,(Bernardo,M.;Berger,J.O.;Dawid,A.P.;Smtith,A.F.M.,Bayesian Statistics 6,牛津大学出版社(1999)),415-440·Zbl 0957.62019号
[30] Malyshkina,N.,《马尔可夫转换模型:道路安全应用》,普渡大学博士论文(2008)
[31] 阿克巴拉巴迪,M。;博尔赫斯,M。;简,A。;佩雷拉,F。;Piri,M.,验证地下水流模型的贝叶斯框架:合成实验,计算机。地质科学。,19, 1231-1250 (2015) ·Zbl 1390.86044号
[32] Smith,B.J.,Boa:MCMC输出收敛评估和后验推断的R包,J.Stat.Softw。,21, 1-37 (2007)
[33] 特纳,B.M。;布朗,S.D。;Sederberg,P.B。;Steyvers,M.,《从相关维度分布中有效抽样的方法》,美国心理学杂志。协会,18,368-384(2013)
[34] Neal,R.M.,MCMC使用哈密顿动力学,Handb。马尔可夫链蒙特卡罗,54,113-162(2010)·Zbl 1229.65018号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。