×

时间分数阶扩散波方程反问题的基于多尺度模型降阶的两级集合卡尔曼滤波器。 (英语) Zbl 1416.65163号

综述:集合卡尔曼滤波器(EnKF)已广泛应用于动态系统的状态估计和参数估计,其中观测数据是按时间顺序获得的。为了更新大量EnKF系综样本,需要对正向问题进行非常繁重的模拟。这将降低EnKF对大型和高维模型的分析效率。为了减少不确定性和加速后验推理,提出了一种两阶段集成卡尔曼滤波器来改进EnKF的序列分析。众所周知,最终的后整体可能集中在初始前整体支撑的一小部分。如果我们首先通过部分观测建立一个新的先验,并且只在新先验的有效区域上构造一个代理,那么效率会更高。为此,我们使用广义多尺度有限元方法(GMsFEM)构建了一个非常粗糙的模型,并在第一阶段生成了一个新的先验集合。GMsFEM提供了一组在粗块中支持的分层多尺度基函数。这为选择自由度来构建简化模型提供了灵活性和适应性。在第二阶段,我们使用GMsFEM和基于稀疏广义多项式混沌(gPC)的随机配置方法,基于新先验建立初始代理模型。为了改进初始代理模型,我们动态更新代理模型,该模型适用于数据的顺序可用性和更新的分析。两阶段EnKF可以实现比标准EnKF更好的估计,并显著提高集合分析(后验探索)的更新效率。为了提高贝叶斯反问题的适用性和灵活性,我们将两阶段EnKF扩展到非高斯模型和层次模型。本文主要研究多孔介质中时间分数阶扩散波模型,并利用所提出的两阶段EnKF研究其Bayes反问题。通过几个数值算例,验证了考虑渗透率场和源函数参数和结构反演的两阶段EnKF方法的性能。

MSC公司:

65J22型 抽象空间中反问题的数值解法
2015年1月62日 贝叶斯推断
62M20型 随机过程推断和预测

软件:

贝叶斯DA
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] Bondarenko,A.N。;Ivaschenko,D.S.,求解变系数时间分数阶扩散方程反问题的数值方法,J.逆病态问题。,17, 419-440 (2009) ·兹比尔1205.26009
[2] Duane,S。;A.D.肯尼迪。;彭德尔顿,B.J。;Roweth,D.,混合蒙特卡罗,物理。莱特。B、 195216-222(1987)
[3] 伊芬迪耶夫,Y。;Datta-Gupta,A。;金廷,V。;马,X。;Mallick,B.,《动态数据集成的有效两阶段马尔可夫链蒙特卡罗方法》,《水资源》。决议,41(2005)
[4] 伊芬迪耶夫,Y。;加尔维斯,J。;Hou,T.Y.,广义多尺度有限元方法(GMsFEM),J.Compute。物理。,251, 116-135 (2013) ·Zbl 1349.65617号
[5] 埃默里克,A.A。;Reynolds,A.C.,利用多重数据同化实现更平滑的集成,计算。地质科学。,55, 3-15 (2013)
[6] Ernst,O.G。;斯普伦克,B。;Starkloff,H.J.,贝叶斯反问题和卡尔曼滤波器,(从复杂系统中提取可量化信息。从复杂系统提取可量化的信息,Lect.Notes Compute.Sci.Eng.(2014),Springer),133-159·Zbl 1328.93260号
[7] Ernst,O.G。;斯普伦克,B。;Starkloff,H.J.,贝叶斯反问题中的集合和多项式混沌卡尔曼滤波器分析,SIAM/ASA J.不确定性。量化。,3823-851(2015年)·Zbl 1339.60041号
[8] Evensen,G.,用非线性准营养模型Issing-Monte-Carlo方法预测误差统计的序贯数据同化,J.Geophys。海洋研究,99,10143-10162(1994)
[9] Gelman,A。;Carlin,J.B。;斯特恩,H.S。;Dunson,D.B.,贝叶斯数据分析(2013),查普曼和霍尔/CRC
[10] 吉尔克斯,W.R。;理查森,S。;Spiegelhalter,D.J.,《马尔可夫链蒙特卡罗实践》(1996),查普曼和霍尔:查普曼与霍尔伦敦·Zbl 0832.00018号
[11] Hansen,P.C.,利用L曲线分析离散不适定问题,SIAM Rev.,34561-580(1992)·Zbl 0770.65026号
[12] Hou,T.Y。;Wu,X.H.,复合材料和多孔介质椭圆问题的多尺度有限元方法,J.Compute。物理。,134, 169-189 (1997) ·Zbl 0880.73065号
[13] Iglesias,医学硕士。;法律,K。;Stuart,A.M.,反问题的集合卡尔曼方法,反问题。,29,第045001条pp.(2013)·兹比尔1311.65064
[14] 江,L。;Li,Q.,基于简化混合GMsFE基方法的模型稀疏表示,J.Compute。物理。,338, 285-312 (2017) ·Zbl 1415.65259号
[15] 江,L。;Ou,N.,地下水流贝叶斯反问题的多尺度模型简化方法,J.Compute。申请。数学。,319, 188-209 (2017) ·兹比尔1524.65478
[16] 江,L。;Ou,N.,基于时间分数扩散方程的粗多尺度模型的中间分布贝叶斯推断,多尺度模型。模拟。,16, 327-355 (2018) ·Zbl 1393.65032号
[17] 凯皮奥,J。;Somersalo,E.,《统计计算反问题》(2006),Springer Science&Business Media·Zbl 0927.35134号
[18] 李,C。;Zeng,F.,分数微积分的数值方法(2015),Chapman&Hall/CRC·Zbl 1326.65033号
[19] 李,C。;赵,Z。;陈永清,带次扩散和超扩散的非线性分数阶微分方程的数值逼近,计算。数学。申请。,62, 855-875 (2011) ·Zbl 1228.65190号
[20] 李·G。;顾伟(Gu,W.)。;Jia,X.,时间分数阶扩散方程中空间相关扩散系数的数值反演,J.逆病态问题。,20, 339-366 (2012) ·Zbl 1279.35091号
[21] 李,J。;Marzouk,Y.M.,逆向问题贝叶斯解代用品的自适应构造,SIAM J.Sci。计算。,36,A1163-A1186(2014)·Zbl 1415.65009号
[22] 李,J。;Xiu,D.,一种基于广义多项式混沌的高精度集合卡尔曼滤波器,J.Compute。物理。,228, 5454-5469 (2009) ·Zbl 1280.93084号
[23] Liu,J.S。;Chen,R.,动态系统的序贯蒙特卡罗方法,美国统计协会,93,1032-1044(1998)·Zbl 1064.65500号
[24] 罗,X。;Moroz,I.M.,《带无味变换的集成卡尔曼滤波器》,Physica D,238549-562(2009)·兹比尔1159.93358
[25] Mallat,S.,《信号处理的小波之旅:稀疏方法》(2008),学术出版社
[26] 马丁·J。;Wilcox,L.C.公司。;伯斯特德,C。;Ghattas,O.,《大规模统计反演问题的随机牛顿MCMC方法及其在地震反演中的应用》,SIAM J.Sci。计算。,34,A1460-A1487(2012)·Zbl 1250.65011号
[27] 米勒,L。;Yamamoto,M.,分数阶扩散方程的系数反问题,逆问题。,29,第075013条pp.(2013)·Zbl 1278.35268号
[28] Miller,K.S。;Ross,B.,分数微积分和分数微分方程导论(1993),Wiley Intersci。出版物·Zbl 0789.26002号
[29] Oldham,K。;Spanier,J.,《分数阶微积分理论及其在任意阶微分和积分中的应用》(1974),Elsevier Science·Zbl 0292.26011号
[30] Restrepo,P。;拉科维奇,O。;Noh,S.J.,《在业务水文预测中推进数据同化:进展、挑战和新机遇》,《水文》。地球系统。科学。,2012年6月16日
[31] 萨夫达里·瓦伊加尼,A。;Heryudono,A。;Larsson,E.,金融应用中产生的对流扩散方程的径向基函数单位配点法,J.Sci。计算。,64, 341-367 (2015) ·Zbl 1325.65139号
[32] 席林斯,C。;Stuart,A.M.,《反问题的集合卡尔曼滤波器分析》,SIAM J.Numer。分析。,55, 1264-1290 (2017) ·Zbl 1366.65101号
[33] 塞勒,A。;埃文森,G。;Skjervheim,J.A.,使用EnKF进行复杂油藏模型的历史拟合和不确定性量化,(大规模反演问题的计算方法和不确定性的量化(2011)),247-271
[34] Stuart,A.M.,《逆向问题:贝叶斯观点》,《数值学报》。,19, 451-559 (2010) ·Zbl 1242.65142号
[35] Vogel,C.R.,《反问题的计算方法》(2002),SIAM:SIAM Philadelphia,PA·Zbl 1008.65103号
[36] Wyss,W.,分数扩散方程,J.Math。物理。,27, 2782-2785 (1986) ·Zbl 0632.35031号
[37] 秀,D。;Karniadakis,G.,随机微分方程的Wiener-Askey多项式混沌,SIAM J.Sci。计算。,24, 619-644 (2002) ·Zbl 1014.65004号
[38] Yan,L。;郭,L.,反问题贝叶斯解的(L_1)最小化随机配置算法,SIAM J.Sci。计算。,37,A1410-A1435(2015)·Zbl 1328.65200号
[39] 周,H。;戈梅斯·赫南德斯,J.J。;Franssen,H.J.H.,处理集合卡尔曼滤波中参数和状态变量的非高斯性的方法,《水资源研究进展》,34,844-864(2011)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。