×

地震层析成像高维反演问题的贝叶斯线性模型。 (英语) Zbl 1288.62045号

摘要:我们将线性贝叶斯模型应用于地震层析成像,这是地球物理中的一个高维反演问题。其目的是根据地表测量数据估计地球内部的三维结构。由于这通常涉及估计数千个或更多未知值,因此一直被视为线性(化)优化问题。
我们提出了一个贝叶斯层次模型来估计地球结构和震源参数的联合分布。椭圆空间先验允许适应地幔的分层性质。利用我们的高效算法,我们可以对大规模线性反问题的后验分布进行采样,并根据给定数据的参数分布和可信区间提供精确的不确定性量化。我们将该方法应用于一个成熟的层析成像问题,即北美西部上部构造的反演,涉及11000多个参数。在对模拟数据和实际数据的研究中,我们表明,我们的方法检索了地球内部的主要结构以及经典最小二乘最小化,同时还提供了不确定性评估。

MSC公司:

2015年1月62日 贝叶斯推断
第62页,第35页 统计学在物理学中的应用
86A32型 地理统计学
62小时12分 多元分析中的估计
65C60个 统计中的计算问题(MSC2010)
62J05型 线性回归;混合模型
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] Aki,K.和Lee,W.(1976年)。利用局部地震的初至P时间确定地震阵列下的三维速度异常:1。同质初始模型。地球物理研究杂志81 4381-4399。
[2] Amestoy,P.R.、Davis,T.A.和Duff,I.S.(1996)。一种近似最小度排序算法。SIAM J.矩阵分析。申请。17 886-905. ·Zbl 0861.65021号 ·doi:10.1137/S0895479894278952
[3] Andrieu,C.和Thoms,J.(2008)。自适应MCMC教程。统计计算。18 343-373. ·doi:10.1007/s11222-008-9110-y
[4] Aster,R.、Borchers,B.和Thurber,C.(2005)。参数估计和反演问题(国际地球物理),Har/cdr编辑,圣地亚哥学术出版社·Zbl 1088.35081号
[5] Banerjee,S.、Gelfand,A.E.和Carlin,B.P.(2003)。空间数据的层次建模与分析。查普曼和霍尔/CRC,佛罗里达州博卡拉顿·Zbl 1053.62105号
[6] Bardsley,J.M.(2012)。基于MCMC的不确定性量化图像重建。SIAM J.科学。计算。34 A1316-A1332·Zbl 1246.15022号 ·doi:10.1137/11085760X文件
[7] Bodin,T.和Sambridge,M.(2009年)。利用可逆跳跃算法进行地震层析成像。国际地球物理杂志178 1411-1436。
[8] Bodin,T.、Sambridge,M.、Tkalć,H.、Arroucau,P.、Gallagher,K.和Rawlinson,N.(2012年a)。接收函数和面波色散的二维反演。地球物理研究杂志117 B02301。
[9] Bodin,T.、Sambridge,M.、Rawlinson,N.和Arroucau,P.(2012b)。未知数据噪声的三维层析成像。《国际地球物理杂志》189 1536-1556。
[10] Brooks,S.、Gelman,A.、Jones,G.L.和Meng,X.-L.编辑(2011年)。马尔可夫链蒙特卡罗手册。佛罗里达州博卡拉顿CRC出版社·Zbl 1218.65001号 ·doi:10.1201/b10905
[11] Bui-Thanh,T.、Ghattas,O.和Higdon,D.(2011年)。基于自适应Hessian的非平稳高斯过程响应面方法,用于概率密度近似,并应用于大规模反问题的贝叶斯解。技术报告,德克萨斯大学计算工程与科学研究所,奥斯汀·Zbl 1257.62035号 ·数字对象标识代码:10.1137/10851419
[12] Burdick,S.、Li,C.、Martynov,V.、Cox,T.、Eakins,J.、Mulder,T.,Astiz,L.、Vernon,F.L.、Pavlis,G.L.和van der Hilst,R.D.(2008)。利用全球和美国阵列可运输阵列数据进行旅行时间层析成像,得出北美下方的上地幔非均质性。地震研究快报79 389-392。
[13] Celeux,G.、Forbes,F.、Robert,C.P.和Titterington,D.M.(2006年)。缺失数据模型的偏差信息标准。贝叶斯分析。1 651-673(电子版)·Zbl 1331.62329号 ·doi:10.1214/06-BA122
[14] Congdon,P.(2003)。应用贝叶斯模型。奇切斯特·威利·Zbl 1023.62026号 ·doi:10.1002/0470867159
[15] Crosson,R.S.(1976年)。地震数据的地壳结构建模1。震源和速度参数的同时最小二乘估计。地球物理研究杂志81 3036-3046。
[16] Cui,T.、Fox,C.和O'Sullivan,M.J.(2011)。采用新的自适应延迟验收Metropolis-Hastings算法对大型地热储层模型进行贝叶斯校准。水资源研究47 26页。
[17] Dahlen,F.A.、Hung,S.H.和Nolet,G.(2000年)。有限频率旅行时间的Fréhet核——I.理论。《国际地球物理杂志》141 157-174。
[18] De \ogonek bski,W.(2010)。蒙特卡罗采样地震层析成像。《纯粹和应用地球物理学》167 131-152。
[19] Dziewonski,A.M.和Anderson,D.L.(1981)。初步参考地球模型。《地球和行星内部物理学》25 297-356。
[20] Fichtner,A.、Kennett,B.L.N.、Igel,H.和Bunge,H.P.(2009)。使用伴随方法对澳大利亚地区上地幔结构进行全波形层析成像。《国际地球物理杂志》179 1703-1725。
[21] Flath,H.P.、Wilcox,L.C.、Akçelik,V.、Hill,J.、van Bloemen Waanders,B.和Ghattas,O.(2011年)。基于低秩部分Hessian近似的大规模线性反问题中贝叶斯不确定性量化的快速算法。SIAM J.科学。计算。33 407-432. ·Zbl 1229.65174号 ·数字对象标识代码:10.1137/090780717
[22] Frühwirth-Schnatter,S.和Pyne,S.(2010年)。单变量和多元偏正态分布和斜(t)分布有限混合的贝叶斯推断。生物统计学11 317-336。
[23] Haario,H.、Saksman,E.和Tamminen,J.(2001)。自适应大都会算法。伯努利7 223-242·Zbl 0989.65004号 ·doi:10.2307/3318737
[24] Haario,H.、Laine,M.、Mira,A.和Saksman,E.(2006年)。DRAM:高效自适应MCMC。统计计算。16 339-354. ·doi:10.1007/s11222-006-9438-0
[25] Hämäläinen,m.S.和Ilmoniemi,R.J.(1994)。解释大脑的磁场:最小范数估计。医学和生物工程与计算32 35-42。
[26] Herbei,R.、McKeague,I.W.和Speer,K.G.(2008)。环流和喷流:海洋环流结构的示踪数据反演。物理海洋学杂志38 1180-1202。
[27] Hoerl,E.和Kennard,R.W.(1970年)。岭回归:非正交问题的有偏估计。技术计量12 55-67·Zbl 0202.17205号 ·数字对象标识代码:10.2307/1267351
[28] Kaipio,J.和Somersalo,E.(2007年)。统计反问题:离散化、模型简化和反犯罪。J.计算。申请。数学。198 493-504. ·Zbl 1101.65008号 ·doi:10.1016/j.cam.2005.09.027
[29] Kass,R.E.、Carlin,B.P.、Gelman,A.和Neal,R.M.(1998)。实践中的马尔可夫链蒙特卡罗:圆桌讨论。阿默尔。统计师。52 93-100.
[30] Kennett,B.L.N.和Engdahl,E.R.(1991)。全球地震位置和相位识别的行程时间。国际地球物理杂志105 429-465。
[31] Khan,A.、Zunino,A.和Deschamps,F.(2011年)。根据表面波相速度的贝叶斯反演,北美大陆下方的热化学和物理结构。地球物理研究杂志116 23页。
[32] Kolehmainen,V.、Vanne,A.、Siltanen,S.、Järvenpää,S.,Kaipio,J.P.、Lassas,M.和Kalke,M.(2007年)。三维牙科X射线成像的贝叶斯反演方法。e&i Elektrotechnik und Informationstechnik 124 248-253。
[33] Lieberman,C.、Willcox,K.和Ghattas,O.(2010年)。大规模统计反问题的参数和状态模型约简。SIAM J.科学。计算。32 2523-2542. ·Zbl 1217.65123号 ·doi:10.1137/090775622
[34] Liu,J.S.(2008)。科学计算中的蒙特卡罗策略。纽约州施普林格·Zbl 1132.65003号
[35] Martin,J.、Wilcox,L.C.、Burstedde,C.和Ghattas,O.(2012年)。大规模统计反演问题的随机牛顿MCMC方法及其在地震反演中的应用。SIAM J.科学。计算。34 A1460-A1487·兹比尔1250.65011 ·数字对象标识代码:10.1137/10845598
[36] Mercerat,D.和Nolet,G.(2013年)。关于有限频率层析成像中互相关延迟时间的线性。国际地球物理杂志192 681-687。
[37] Mohammad-Djafari,A.(2012年)。具有先验模型的贝叶斯方法,在信号和图像处理中增强稀疏性。2012年《EURASIP信号处理进展杂志》52。
[38] Montelli,R.、Nolet,G.、Dahlen,F.A.、Masters,G.,Engdahl,E.R.和Hung,S.H.(2004)。有限频率层析成像揭示了地幔中的各种羽流。科学303 338-343。
[39] Mosca,I.、Cobden,L.、Deuss,A.、Ritsema,J.和Trampert,J.(2012)。概率层析成像下地幔的地震和矿物学结构。地球物理研究杂志117 26页。
[40] Mosegaard,K.和Tarantola,A.(1995年)。反问题解的蒙特卡罗抽样。地球物理研究杂志100 12431-12447。
[41] Mosegaard,K.和Tarantola,A.(2002年)。《国际地震与工程地震学手册:反演问题的概率方法》237-265。圣地亚哥学术出版社。
[42] Neal,R.M.(2010)。MCMC使用哈密顿动力学。《马尔可夫链蒙特卡罗手册》(S.Brooks、A.Gelman、G.Jones和X.L.Meng编辑)。查普曼和霍尔/CRC,佛罗里达州博卡拉顿·Zbl 1229.65018号
[43] Nolet,G.(2008)。地震层析成像简介:地球和太阳内部成像。剑桥大学出版社,剑桥·Zbl 1166.86001号
[44] Paige,C.C.和Saunders,M.A.(1982年)。LSQR:稀疏线性方程组和稀疏最小二乘的算法。ACM事务处理。数学。软件8 43-71·Zbl 0478.65016号 ·数字对象标识代码:10.1145/355984.355989
[45] Pettitt,A.N.、Weir,I.S.和Hart,A.G.(2002年)。不规则间隔多元数据的条件自回归高斯过程,用于建模大型二进制数据集。统计计算。12 353-367. ·doi:10.1023/A:1020792130229
[46] Rue,H.和Held,L.(2005)。高斯马尔可夫随机场:理论与应用。统计学和应用概率专著104。查普曼和霍尔/CRC,佛罗里达州博卡拉顿·邮编1093.60003
[47] Sahu,S.K.、Dey,D.K.和Branco,M.D.(2003)。一类新的多元偏斜分布及其在贝叶斯回归模型中的应用。加拿大。J.统计。31 129-150. ·Zbl 1039.62047号 ·doi:10.2307/3316064
[48] Sambridge,M.和Mosegaard,K.(2002年)。地球物理反演问题中的蒙特卡罗方法。地球物理学评论40 1-29·Zbl 1015.65001号 ·doi:10.1088/0266-5611/18/3/201
[49] Schmandt,B.和Humphreys,E.(2010年)。美国西部上地幔体波层析成像揭示了复杂的俯冲和小规模对流。《地球与行星科学快报》297 435-445。
[50] Sieminski,A.、Liu,Q.、Trampert,J.和Tromp,J.(2007)。基于伴随方法的体波对各向异性的有限频率敏感性。《国际地球物理杂志》171 368-389。
[51] Sigloch,K.(2008)。多频体波层析成像。新泽西州普林斯顿大学博士论文。
[52] Sigloch,K.(2011)。来自北美地幔省的多频P波层析成像。地球化学、地球物理、地球系统12 Q02W08。
[53] Sigloch,K.、McQuarrie,N.和Nolet,G.(2008)。根据有限频率层析成像推断北美洲下两阶段俯冲历史。自然地理1 458-462。
[54] Sigloch,K.和Nolet,G.(2006年)。测量有限频率体波振幅和旅行时间。国际地球物理杂志167 271-287。
[55] Simpson,D.P.、Turner,I.W.和Pettitt,A.N.(2008)。使用Krylov子空间方法从高斯马尔可夫随机场快速采样。澳大利亚布里斯班昆士兰大学QUT ePrints ID 14376。
[56] Spiegelhalter,D.J.、Best,N.G.、Carlin,B.P.和van der Linde,A.(2002)。模型复杂度和拟合度的贝叶斯度量。J.R.统计社会服务。B统计方法。64 583-639. ·Zbl 1067.62010年 ·数字对象标识代码:10.1111/1467-9868.00353
[57] Stramer,O.和Tweedie,R.L.(1999)。Langevin型模型。二、。MCMC算法的自我定位候选。Methodol公司。计算。申请。普罗巴伯。1 307-328. ·兹标0946.60063 ·doi:10.1023/A:101090512027
[58] Swindel,B.F.(1976年)。基于先验信息的良好岭估计。通信统计。理论方法A5 1065-1075·Zbl 0342.62035号 ·doi:10.1080/03610927608827423
[59] Tarantola,A.(2004)。模型参数估计的反问题理论和方法,第1版,SIAM,费城·Zbl 1111.62025号
[60] Ter Braak,C.J.F.(2006)。遗传算法差分进化的马尔可夫链蒙特卡罗版本:实际参数空间的简易贝叶斯计算。统计计算。16 239-249. ·doi:10.1007/s11222-006-8769-1
[61] Tian,Y.、Sigloch,K.和Nolet,G.(2009)。美国西部上地幔的多频SH波层析成像。国际地球物理杂志178 1384-1402。
[62] Tian,Y.、Montelli,R.、Nolet,G.和Dahlen,F.A.(2007年)。在有限频率层析成像中计算旅行时和振幅灵敏度核。J.计算。物理学。226 2271-2288. ·Zbl 1147.86300号 ·doi:10.1016/j.jcp.2007.07.004
[63] Tromp,J.、Tape,C.和Liu,Q.(2005)。地震层析成像,伴随方法,时间反转和巴纳花生核。国际地球物理杂志160 195-216。
[64] Uutela,K.、Häämäläinen,m.和Somersalo,E.(1999年)。使用最小电流估计值可视化脑磁图数据。神经影像10 173-180。
[65] Vrugt,J.A.、ter Braak,C.J.F.、Diks,C.G.H.、Robinson,B.A.、Hyman,J.M.和Higdon,D.(2009)。通过自适应随机子空间采样的差分进化加速马尔可夫链蒙特卡罗模拟。国际非线性科学与数值模拟杂志10 273-290·Zbl 06942400号
[66] Watzenig,D.和Fox,C.(2009年)。电容层析成像的统计建模和推断综述。测量科学与技术20 22页。
[67] Wilkinson,D.J.和Yeung,S.K.H.(2002年)。高结构高斯系统的条件模拟,应用于块MCMC,用于超大线性模型的贝叶斯分析。统计计算。12 287-300. ·doi:10.1023/A:1020711129064
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。