×

一种用于弹性波问题三维建模的对比源积分方程方法。 (英语) Zbl 1360.74088号

小结:本文数值实现了最近提出的求解三维弹性波散射问题的积分方程公式。该方法由应力张量和质点速度矢量表示,其中二阶对称张量被分解为其全向分量和偏差分量。随后,利用该积分方程得到了对比源型积分方程。为了求解这些积分方程,我们采用了共轭梯度快速傅里叶变换(CG-FFT)方案,该方案基于提供(二阶)精确近似的求积公式,同时保留了相关积分的卷积性质,使其易于通过快速傅里尔变换进行有效评估。作为线性解算器,我们采用了共轭梯度法残差(CGNR)格式,它总是单调收敛的,但收敛速度较慢,而双共轭梯度稳定(BiCGSTAB)方案效率更高,但稳定性较差。通过使用简单的对角预条件,进一步提高了迭代格式的收敛速度。数值结果表明了所提方法的准确性和效率。

MSC公司:

74J20型 固体力学中的波散射
74B05型 经典线性弹性
45B05型 弗雷德霍姆积分方程
第74页第30页 固体力学中的其他数值方法(MSC2010)
65兰特 积分方程的数值方法
65T50型 离散和快速傅里叶变换的数值方法
PDF格式BibTeX公司 XML格式引用
全文: DOI程序 链接

参考文献:

[1] Aki,K。;Richards,P.G.,《定量地震学:理论和方法》(1998),弗里曼:弗里曼旧金山
[2] 维利厄,J。;Operto,S.,《勘探地球物理学中的全波形反演概述》,地球物理学,74,6,WCC1-WCC26(2009)
[3] 尤因,W.M。;医学硕士Jardetzky。;Press,F.,《分层介质中的弹性波》(1957),McGraw-Hill图书公司:纽约McGraw-Hill图书公司(第1章)·Zbl 0083.23705号
[4] 布罗西耶,R。;艾蒂安,V。;Operto,S。;Virieux,J.,用有限差分和有限元间断Galerkin方法对粘声波进行频域数值模拟,(Dissanayake,D.W.,声波(2010),SCIYO:SCIYO克罗地亚),125-158,(第7章)
[5] Robertson,J.O。;贝德纳尔,B。;布兰奇,J。;科斯托夫,C。;van Manen,D.J.,《地震建模与采集、处理和解释应用补充介绍》,《地球物理学》,72,5,SM1-SM4(2007)
[6] Chironi,C。;摩根大通。;Warner,M.R.,《利用全波场地震层析成像技术对玄武岩和底盐结构进行成像》,J.Geophys。Res.,111,B05313(2006),18页
[7] Abe,S。;拉瑟姆,S。;Mora,P.,粗糙平面断层三维颗粒模拟中的动态破裂,Pure Appl。地球物理学。,163, 9, 1881-1892 (2006)
[8] Carcione,J.M.(卡尔乔内,J.M.)。;赫尔曼,G.C。;ten Kroode,A.P.E.,地震建模,地球物理学,67,1304-1325(2002)
[9] 图黑,T。;Kiuchi,T。;Iwasaki,K.,弹性波在半空间中传播的快速体积积分方程法,国际期刊系统。科学。,48, 3194-3208 (2011)
[10] de Hoop,A.T.,《波的辐射和散射手册》(1995),学术出版社:圣地亚哥学术出版社,(第10章)
[11] 杨,J。;阿布巴卡尔,A。;范登伯格,P.M。;哈巴西,T。;Reitich,F.,解决三维弹性散射问题的应力-速度公式的CG-FFT方法,J.Comput。物理。,227, 10018-10039 (2008) ·Zbl 1154.74051号
[12] Eymé-Bellegarda,E.J。;Habashy,T.M.,《非均匀弹性背景中的多维超声散射》,J.Acoust。《美国社会杂志》,91,6,3104-3115(1992)
[13] 图黑,T。;Kiuchi,T。;Iwasaki,K.,弹性波散射现象中正/反问题的快速体积积分方程方法,国际期刊系统。科学。,46, 3860-3872 (2009) ·Zbl 1176.74099号
[14] de Hoop,A.T。;阿布巴卡尔,A。;Habashy,T.M.,《三维时域弹性动力散射问题的对比源应力-速度积分方程公式:使用张量划分的结构化方法》,J.Acoust。《美国社会》,126,3,1095-1100(2009)
[15] Richmond,J.H.,《任意横截面形状的介质圆柱的散射》,IEEE Trans。天线与传播,13334-341(1965)
[16] Vainikko,G.,多维弱奇异积分方程(1993),Springer-Verlag:Springer-Verlag纽约·兹比尔0789.65097
[17] Zwamborn,P。;van den Berg,P.M.,求解散射问题的共轭梯度FFT方法的三维弱形式,IEEE Trans。微型。《理论技术》,40,1757-1766(1992)
[18] 巴雷特·R。;贝里,M。;Chan,T。;德梅尔,J。;J.多纳托。;Dongarra,J。;埃伊霍特,V。;波佐,R。;罗明,C。;van der Vorst,H.,《线性系统解的模板:迭代方法的构建块》(1994),SIAM:SIAM Philadelphia
[19] Greenbaum,A.,求解线性系统的迭代方法(1997),SIAM:SIAM Philadelphia,(第8章)·Zbl 0883.65022号
[20] van der Vorst,H.A.,BiCGSTAB:非对称线性系统解的BICG的一种快速平滑收敛变体,SIAM J.Sci。统计计算。,13, 2, 631-644 (1992) ·Zbl 0761.65023号
[21] Gutknecht,M.H.,《复谱矩阵的BiCGSTAB变体》,SIAM J.Sci。计算。,14, 5, 1020-1033 (1993) ·兹比尔083765031
[22] 阿布巴卡尔,A。;van den Berg,P.M.,《应用于井间感应传感器的三维逆散射》,IEEE Trans。地质科学。遥感,38,4,1669-1681(2000)
[23] 阿布巴卡尔,A。;van den Berg,P.M.,基于三维电磁物体的区域积分方程的迭代正演和逆演算法,J.Compute。物理。,195, 236-262 (2004) ·Zbl 1086.78006号
[24] 范登伯格,P.M。;Abubakar,A.,基于对比源积分表示的逆散射算法,逆概率。科学。,10, 6, 559-576 (2002)
[25] 佩莱卡诺斯,G。;阿布巴卡尔,A。;van den Berg,P.M.,《弹性动力学中的对比源反演方法》,J.Acoust。《美国社会杂志》,114,52825-2834(2003)
[26] 范登伯格,P.M。;Kleinman,R.E.,对比源反演方法,反演问题,13,6,1607-1620(1997)·Zbl 0945.35090号
[27] Eymé-Bellegarda,E.J。;Habashy,T.M.,埋在多层弹性背景中多维物体的超声逆散射,IEEE Trans。超声波。铁电极。频率控制,39,1,11-18(1992)
[28] Eymé-Bellegarda,E.J。;哈巴西,T.M。;Bellegarda,J.R.,层状弹性背景介质中埋藏多维物体的自动识别,Opt。工程师,31,12,2572-2579(1992)
[29] Martin,G.S.,Marmousi2:Marmousi的弹性升级,The Leading Edge,25,2,156-166(2006)
[30] L.Vaillant,SEG-EAGE盐数据,技术代表,斯坦福勘探项目,2002年。;L.Vaillant,SEG-EAGE盐数据,技术代表,斯坦福勘探项目,2002年。
[31] 出版社,W.H。;Teukolsky,S.A。;韦特林,W.T。;弗兰纳里,B.P.,《FORTRAN中的数字公式——科学计算的艺术》(1992),剑桥大学出版社:剑桥大学出版社·Zbl 0778.65002号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。