杨鹏良；罗曼·布罗西耶；Métiver，卢多维奇；让·维利厄；周伟（Zhou，Wei） 粘声多参数全波形反演的时域预处理截断牛顿方法。 （英文） Zbl 1394.65045号 SIAM J.科学。计算。 40，编号4，B1101-B1130（2018）. 摘要：基于二阶伴随态公式，发展了一种用于粘声介质中时域全波形反演的截断牛顿（TRN）方法。时域梯度估计和Hessian矢量积是通过重新计算入射波场和伴随波场来管理的，在没有数值不稳定性的情况下，同时对入射波场进行重新计算，以缓解记忆问题。由于链式规则，提出了通用算法工作流来在参数化之间切换。为了提高内共轭梯度迭代的收敛速度，提出了一种适用于多参数结构的高效预条件器。预处理器中引入了额外的用户定义缩放，以缓解特定参数对波形变化的弱敏感性。通过一个玩具示例验证了逆Hessian对缓解参数间权衡的重要性。通过一个基于北海实际数据应用的真实二维合成案例，在（V_p，\rho，Q^{-1}）参数化下的令人鼓舞的数值结果表明，考虑Hessian影响显著改善了多参数重建，减轻了参数之间的耦合，与标准准牛顿优化策略相比，计算成本有了合理的增加。 引用于10文件 MSC公司： 65K10码 数值优化和变分技术 35兰特 PDE的反问题 86甲15 地震学（包括海啸建模）、地震 关键词：地震衰减；地震层析成像；逆理论；波传播 软件：L-BFGS公司；旋转；SEISCOPE公司 PDF格式BibTeX公司 XML格式引用 \textit{P.Yang}等人，SIAM J.Sci。计算。40，第4号，B1101--B1130（2018；Zbl 1394.65045） 全文： 内政部 参考文献： [1] T.Alkhalifah和R.Plessix，各向异性介质中实用全波形反演方法：解析参数分辨率研究《地球物理学》，79（2014），第R91–R101页。 [2] J.E.Anderson、L.Tan和D.Wang，RTM和FWI的时间覆盖检查点方法《地球物理学》，77（2012），第S93–S103页。 [3] A.Askan、V.Akcelik、J.Bielak和O.Ghattas，非均质结构中地震速度和滞弹性损失的全波形反演，公牛。地震。《美国社会杂志》，97（2007），第1990-2008页。 [4] O.Barkved、P.Heavey、J.H.Kommedal、J.-P.van Gestel、R.S.ynnöve、H.Pettersen、C.Kent和U.Albertin，Valhall全波形反演的业务影响，《SEG技术计划扩展摘要》，2010年，第925–929页。 [5] C.Boehm、M.Hanzich、J.de la Puente和A.Fichtner，全波形反演中伴随方法的波场压缩《地球物理学》，81（2016），第R385–R397页。 [6] C.Boehm和M.Ulbrich，地震层析成像中约束参数识别的半光滑Newton-CG方法，SIAM J.科学。计算。，37（2015），第S334–S364页·Zbl 1325.86001号 [7] Y.Capdeville、L.Guillot和J.-J.Marigo，P-SV波在高尺度弹性介质中的二维非周期均匀化，地球物理。《国际期刊》，182（2010），第903–922页。 [8] J.Carcione、D.Kosloff和R.Kosloff，线性粘声介质中的波传播模拟，地球物理。《国际期刊》，93（1988），第393-401页·Zbl 0635.73037号 [9] J.Carcione、D.Kosloff和R.Kosloff，弹性各向异性（横观各向同性）固体中的波传播模拟，夸脱。J.机械。申请。数学。，41（1988），第319-345页·Zbl 0656.73020号 [10] E.Causse、R.Mittet和B.Ursin，粘声波介质中全波形反演的预处理《地球物理学》，64（1999），第130–145页。 [11] C.Cerjan、D.Kosloff、R.Kosloff和M.Reshef，离散声波和弹性波方程的无反射边界条件《地球物理学》，50（1985），第2117–2131页。 [12] G.Chavent，反问题的非线性最小二乘法，施普林格，多德雷赫特，海德堡，伦敦，纽约，2009年·Zbl 1191.65062号 [13] P.Chen、T.Jordan和L.Zhao，全三维层析成像：散射积分和伴随波场方法的比较，地球物理。《国际期刊》，170（2007），第175–181页。 [14] Y.Choi和C.Shin，基于新伪海西矩阵的频域弹性全波形反演：弹性Marmousi（2）合成数据的经验，公牛。地震波。《美国社会杂志》，98（2008），第2402–2415页。 [15] S.M.Day和J.B.Minster，用Pade近似方法数值模拟衰减波场，地球物理。《国际期刊》，78（1984），第105–118页·Zbl 0587.73150号 [16] E.Duveneck和P.M.Bakker，倾斜TI介质中反向偏移的稳定P波模拟《地球物理学》，76（2011），第S65–S75页。 [17] S.C.Eisenstat和H.F.Walker，用不精确牛顿法选择强迫项，SIAM J.科学。计算。，17（1996），第16-32页·Zbl 0845.65021号 [18] H.Emmerich和M.Korn，在地震波场时域计算中引入衰减《地球物理学》，52（1987），第1252-1264页。 [19] I.Epanomeritakis、V.Akçelik、O.Ghattas和J.Bielak，大规模三维弹性全波形地震反演的Newton-CG方法《反问题》，24（2008），第1-26页·Zbl 1142.65052号 [20] A.Fichtner、B.L.N.Kennett、H.Igel和H.P.Bunge，径向各向异性结构的全波形层析成像：对澳大利亚上地幔现在和过去状态的新见解，地球行星。科学。莱特。，290（2010），第270-280页。 [21] A.Fichtner和J.Trampert，基于伴随技术的地震数据泛函Hessian核，地球物理。《国际期刊》，185（2011），第775-798页。 [22] G.H.F.Gardner、L.W.Gardner和A.R.Gregory，地层速度和密度——地层圈闭的诊断基础《地球物理学》，39（1974），第770-780页。 [23] S.Greenhalgh、B.Zhou和A.Green，局部最小搜索地球物理反演的解决方案、算法和相互关系《地球物理学杂志》。Eng.，3（2006），第101-113页。 [24] A.Griewank，在反向自动微分中实现时间和空间复杂性的对数增长，最佳。方法软。，1（1992年），第35-54页。 [25] A.Griewank和A.Walther，算法799:Revolve：计算微分的反向或伴随模式的检查点实现，ACM变速器。数学。《软件》，26（2000），第19-45页·Zbl 1137.65330号 [26] B.Hak和W.A.Mulder，具有因果关系的地震衰减成像，地球物理。《国际期刊》，184（2011），第439–451页。 [27] M.Hinze、R.Pinnau、M.Ulbrich和S.Ulbich，基于PDE约束的优化，数学。模型理论应用。23，施普林格，多德雷赫特，2009年·Zbl 1167.49001号 [28] K.A.Innanen，临界前反射全波形反演中的地震AVO和逆Hessian，地球物理。《国际期刊》，199（2014），第717-734页。 [29] A.Kalogeropoulos、J.van der Kruk、J.Hugenschmidt、S.Busch和K.Merz，用探地雷达全波形反演法评价混凝土中的氯化物和水分，近地表地球物理。，9（2011年），第277-285页。 [30] R.Kamei和R.G.Pratt，粘声波形反演的反演策略，地球物理。《国际期刊》，194（2013），第859-894页。 [31] A.Klotzsche、J.van der Kruk、G.A.Meles和H.Vereecken，含水层系统中作为优先流动路径的波导的井间探地雷达全波形反演《地球物理学》，77（2012），第H57–H62页。 [32] N.Korta、A.Fichtner和V.Sallares，用于全波形反演预处理的块对角近似Hessian，《第75届EAGE会议和展览的扩展摘要》，包括SPE EUROPEC 2013，伦敦，2013。 [33] P.Lailly，作为叠前偏移序列的地震反演问题《逆散射、理论与应用会议》，R.Bednar和Weglein主编，SIAM，费城，1983年，第206-220页。 [34] A.R.Levander，四阶有限差分P-SV地震图，《地球物理学》，53（1988），第1425-1436页。 [35] D.C.Liu和J.Nocedal，大规模优化的有限内存BFGS方法，数学。编程，45（1989），第503-528页·Zbl 0696.90048号 [36] H.-P.Liu、D.L.Anderson和H.Kanamori，滞弹性引起的速度色散；地震学和地幔成分的意义，地球物理。《国际期刊》，47（1976），第41-58页。 [37] J.Martin、L.C.Wilcox、C.Burstede和O.Ghattas，大规模统计反演问题的随机Newton-MMC方法及其在地震反演中的应用，SIAM J.科学。计算。，34（2012），第A1460–A1487页·Zbl 1250.65011号 [38] L.Meítiver、F.Bretaudeau、R.Brossier、S.Operto和J.Virieux，全波形反演和截断牛顿法：复杂地下结构的定量成像，地球物理。前景。，62（2014），第1353-1375页。 [39] L.Meítiver和R.Brossier，SEISCOPE优化工具箱：基于反向通信的大规模非线性优化库《地球物理学》，81（2016），第F11–F25页。 [40] L.Meítiver、R.Brossier、S.Operto和J.Virieux，用于重建纵波速度、密度和衰减的声学多参数FWI：预处理截断牛顿法，见SEG技术计划扩展摘要，2015年，第1198–1203页。 [41] L.Meítiver、R.Brossier和J.Virieux，结合渐近线性化反演和全波形反演，地球物理。《国际期刊》，201（2015），第1682-1703页。 [42] L.Meítiver、R.Brossier、J.Virieux和S.Operto，全波形反演和截断牛顿法，SIAM J.科学。计算。，35（2013），第B401–B437页·Zbl 1266.86002号 [43] J.Minet、S.Lambot、E.C.Slob和M.Vanclooster，薄层存在下全波形探地雷达信号反演估算土壤表层含水量，IEEE传输。地质科学。《遥感》，48（2010），第1138–1150页。 [44] P.Moczo和J.Kristek，地震波传播时域方法中的流变模型，地球物理。Res.Lett.公司。，32（2005），L01306。 [45] J.L.Morales和J.Nocedal，有限记忆准牛顿更新的自动预处理、SIAM J.Optim.、。，10（2000），第1079–1096页·Zbl 1020.65019号 [46] W.A.Mulder和B.Hak，衰减散射成像中的模糊性，地球物理。《国际期刊》，178（2009），第1614-1624页。 [47] S.G.Nash，截断牛顿方法综述，J.计算。申请。数学。，124（2000），第45-59页·Zbl 0969.65054号 [48] J.Nocedal，用有限存储更新拟Newton矩阵，数学。公司。，35（1980），第773-782页·Zbl 0464.65037号 [49] J.Nocedal和S.J.Wright，数值优化第二版，施普林格出版社，纽约，2006年·Zbl 1104.65059号 [50] S.Operto、R.Brossier、Y.Ghoma、L.Meítiver、V.Prieux、A.Ribodetti和J.Virieux，多分量数据多参数全波形反演导览：从理论到实践《前沿》，32（2013），第1040–1054页。 [51] S.Operto、A.Miniussi、R.Brossier、L.Combe、L.Meítiver、V.Monteiller、A.Ribodetti和J.Virieux，海底电缆数据的高效三维频域单参数全波形反演：应用于Valhall的粘声垂直横向各向同性近似，地球物理。《国际期刊》，202（2015），第1362-1391页。 [52] S.Operto、J.Virieux、J.X.Dessa和G.Pascal，基于多倍海底地震仪数据的频域全波形层析成像地壳地震成像：在南开海槽东部的应用《地球物理学杂志》。第111号决议（2006年）。 [53] R.E.普莱西斯，计算泛函梯度的伴随状态方法及其地球物理应用综述，地球物理。《国际期刊》，167（2006），第495-503页。 [54] R.E.Plessix和C.Perkins，深海海底地震仪数据集的全波形反演《第一次突破》，第28页（2010年），第71–78页。 [55] R.G.Pratt、C.Shin和G.J.Hicks，频率空间地震波形反演中的Gauss-Newton和full Newton方法，地球物理学。《国际期刊》，133（1998），第341-362页。 [56] V.Prieux、R.Brossier、Y.Ghoma、S.Operto、J.Virieux、O.Barkved和J.Kommedal，各向同性全波形反演中各向异性的足迹：Valhall案例研究，地球物理。《国际期刊》，187（2011），第1495-1515页。 [57] V.Prieux、R.Brossier、S.Operto和J.Virieux，Valhall多分量OBC数据的多参数全波形反演。第1部分：成像压缩波速、密度和衰减，地球物理。《国际期刊》，194（2013），第1640–1664页。 [58] V.Prieux、R.Brossier、S.Operto和J.Virieux，Valhall多分量OBC数据的多参数全波形反演。第2部分：成像纵波和横波速度，地球物理。《国际期刊》，194（2013），第1665-1681页。 [59] J.O.A.Robertson，地形影响下弹性/粘弹性有限差分建模的数值自由曲面条件《地球物理学》，61（1996），第1921-1934页。 [60] F.Santosa和W.W.Symes，反射地震学反问题最小二乘解的Hessian计算《反问题》，4（1988），第211–233页·Zbl 0642.65052号 [61] C.Shin、S.Jang和D.J.Min，利用逆散射理论改进叠前深度偏移的振幅保持，地球物理。前景。，49（2001），第592-606页。 [62] L.Sirgue、O.I.Barkved、J.Dellinger、J.Etgen、U.Albertin和J.H.Kommedal，全波形反演：瓦尔霍尔成像技术的下一个飞跃《第一次突破》，第28页（2010年），第65-70页。 [63] 孙文华、傅立义，减少逆向时间迁移中数据存储的两种有效方法，计算。和地质科学。，56（2013年），第69-75页。 [64] W.W.Symes，具有最佳检查点的反向时间迁移《地球物理学》，72（2007），第SM213–SM221页。 [65] C.Tape、Q.Liu、A.Maggi和J.Tromp，基于谱元和伴随方法的南加州地壳地震层析成像，地球物理。《国际期刊》，180（2010），第433-462页。 [66] A.塔兰托拉，声学近似中地震反射数据的反演《地球物理学》，49（1984），第1259–1266页。 [67] L.A.Thomsen，弱弹性各向异性《地球物理学》，51（1986），第1954-1966页。 [68] D.Vigh、K.Jiao、D.Watts和D.Sun，地震数据采集多分量测量弹性全波形反演应用《地球物理学》，79（2014），第R63–R77页。 [69] J.Virieux和S.Operto，勘探地球物理全波形反演综述《地球物理学》，74（2009），第WCC1–WCC26页。 [70] Y.Wang、L.Dong、Y.Liu和J.Yang，基于块对角伪海西近似的二维频域弹性全波形反演《地球物理学》，81（2016），第R247–R259页。 [71] M.Warner、A.Ratcliffe、T.Nangoo、J.Morgan、A.Umpleby、N.Shah、V.Vinje、I.Stekl、L.Guasch、C.Win、G.Conroy和A.Bertrand，各向异性\（3\）D全波形反演《地球物理学》，第78页（2013年），第R59–R80页。 [72] P.Yang、R.Brossier、L.Meítiver和J.Virieux，粘弹性介质中三维多参数全波形反演系统公式综述，地球物理。《国际期刊》，207（2016），第129-149页。 [73] P.Yang、R.Brossier、L.Meítiver和J.Virieux，衰减介质中的波场重建：一种检查点辅助的反向模拟方法《地球物理学》，81（2016），第R349–R362页。 [74] P.Yang、R.Brossier和J.Virieux，反向传播波场重建的降采样加插值方法，第78届EAGE会议和展览扩展摘要，2016年，SBT5 08。 [75] P.Yang、R.Brossier和J.Virieux，从大量减少的边界重建波场《地球物理学》，80（2016），第T197–T209页。 [76] W.Zhou、R.Brossier、S.Operto和J.Virieux，基于尺度分离的阻抗反演用于速度建模的潜水波和反射波全波形反演，地球物理。《国际期刊》，202（2015），第1535-1554页。 [77] H.Zhu、E.Bozdağ、D.Peter和J.Tromp，伴随层析成像揭示的欧洲上地幔结构，国家地质科学。，5（2012），第493-498页。 [78] H.Zhu、E.Bozdağ和J.Tromp，基于伴随层析成像的欧洲上地幔地震结构，地球物理。《国际期刊》，201（2015），第18-52页。 此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配，可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。