FMM射线追踪方法在地震学正演和反演中的应用
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
FMM(Fast Marching Method)射线追踪方法是一种求解非线性程函方程数值解的射线追踪方法.本文在FMM射线追踪算法流程的基础上引入一个数据链表,方便求解射线路径,另外还引入一个与解精度有关的方向参数,可以更好的控制解的精度.然后在一些理论模型中对FMM方法的正演和反演计算做了若干验证,结果表明,FMM方法可以较快而且稳定地得到令人满意的结果.
Fast Marching Method (FMM) ray tracing is a ray tracing method using the numerical solutions of Eikonal equation.In this paper a new data array structure and a direction parameter are introduced, which are different from the original FMM method.It is convenient to get the ray path from the new data array structure.The resolution accuracy can be well controlled by the direction parameter.At last some test results are shown based on several velocity models.The resut shows that we can get the content result from the FMM method, and the FMM method is much faster and more stable.
引文
[1] Zhao D, Hasegawa A, Horiuchi S.Tomographic Imaging of P and S Wave Velocity Structure Beneath Northeasetern Japan[J].J.Geophys.Res., 1992, 97(B13) : 19909~19928.
    [2] Zhao D, Hasegawa A, Kanamori H.Deep structure of Japan subduction zone as derived from local, regional, and teleseismic events[J].J.Geophys.Res., 1994, 99(B11) :22313~22329.
    [3] 刘福田,李强,吴华,等.用于速度图象重建的层析成象法[J].地球物理学报,1989,32(1) :46~57. Liu F T, Li Q, Wu H, et al.On the tomographic inverse method used in velocity image reconstruction[J].Chinese J.Geophys.(Acta Geophysia Sinica) (in chinese), 1989, 32(1) : 46~57.
    [4] 雷建设,周蕙兰,赵大鹏.帕米尔及邻区地壳上地幔P波三维 速度结构的研究[J].地球物理学报,2002,45(6) :802~811. Lei J S, Zhou H L, Zhao D P.3-D velocity structure of Pwave in the crust and upper-mantle beneath Pamir and adjacent region[J].Chinese J.Geophys.(in chinese), 2002, 45(6) :802~811.
    [5] 王椿镛,Mooney W D,王溪莉,等.川滇地区地壳上地幔三维速度结构研究[J].地震学报,2002,24(1) :1~16. Wang C Y, Mooney W D, Wang X L, et al.Study on 3-D velocity structure of crust and upper mantle in Sichuan-Yunnan region[J],China.Acta Seismologica Sinica, 2002, 24(1) :1~16.
    [6] 丁志峰,何正勤,孙为国,等.青藏高原东部及其边缘地区的地壳上地幔三维速度结构[J].地球物理学报,1999,42(2) :197~205. Ding Z F, He Z Q, Sui W G, et al.3-D crust and upper mantle velocity structure in eastern Tibetan plateau and its surrounding areas[J].Chinese J.Geophys.(in chinese),1999, 42(2) : 197~205.
    [7] 胥颐,刘建华,郝天珧,等.中国东部海域及邻区岩石层地幔的P波速度结构与构造分析[J].地球物理学报,2006,49(4) :1053~1061.Xu Y, Liu J H, Hao T Y, et al.P wave velocity structure and tectonics analysis of lithospheric mantle in eastern China seas and adjacent regions[J].Chinese J.Geophys.(in Chinese),2006, 49(4) : 1053~1061.
    [8] 徐佩芬,刘福田,王清晨,等.大别-苏鲁碰撞造山带的地震层析成像研究--岩石圈三维速度结构[J].地球物理学报,2000,43(3) :377~386.Xu P F, Liu F T, Wang Q C, et al.Seismic tomography beneath the Dabie-Sulu collision orogen--3-D velocity structures of lithosphere[J].Chinese J.Geophys.( in Chinese), 2000, 43(3) :377~386.
    [9] 成谷,马在田,张宝金,等.地震层析成像中存在的主要问题及应对策略[J].地球物理学进展.2003,18(3) :512~518. Cheng G, Ma Z T, Zhang B J, et al.Primary problems and according strategies in seismic tomography[J].Progress in Geophysics, 2003, 18(3) : 512~518.
    [10] Huang J L, Zhao D P.High-resolution mantle tomography of China and surrounding regions[J].J.Geophys.Res.,2006, 111(B09305) .
    [11] Rawlinson N, Kennett B L N, Heintz M.Insights into the structure of the upper mantle beneath the Murray Basin from 3D teleseismic tomography[J].Australian Journal of Earth Sciences, 2006, 53: 595~604.
    [12] Rawlinson N, Reading A M, Kennett B L N.Lithospheric structure of Tasmania from a novel form of teleseismic tomography[J].J.Geophys.Res.2006, 111(B02301) .
    [13] Julian B R, Gubbins D.Three-dimensional seismic ray tracing[J].J.Geophys., 1977, 43(1) : 95~113.
    [14] Thurber C H, Ellsworth W L.Rapid solution of ray tracing problems in heterogeneous media[J].Bull.Seis.Soc.Am., 1980, 70(4) : 1137~1148.
    [15] 王辉,常旭.基于图形结构的三维射线追踪方法[J].地球物理学报,2000,43(4) :534~541. Wang H, Chang X.3-D ray tracing method based on graphic structure[J].Chinese J.Geophys.(in chinese), 2000, 43(4) : 534~541.
    [16] 张建中,陈世军,徐初伟.动态网格最短路径射线追踪[J].地球物理学报,2004,47(5) :899~904. Zhang J Z, Chen S J, Xu C W.A method of shortest path ray tracing with dynamic networks[J].Chinese J.Geophys.(in Chinese), 2004, 47(5) : 899~904.
    [17] 张美根,程冰洁,李小凡,等.一种最短路径射线追踪的快速算法[J].地球物理学报.2006,49(5) :1467~1474. Zhang M G, Cheng B J, Li X F, et al.A fast algorithm of shortest path ray tracing[J].Chinese J.Geophys.(in chinese), 2006, 49(5) : 1467~1474.
    [18] 赵爱华,张中杰,王光杰,等.非均匀介质中地震波走时与射线路径快速计算技术[J].地震学报,2000,22(2) :151~157.Zhao A H, Zhang Z J, Wang G J, et al.A new scheme for fast calculation of seismic traveltimes and ray paths in heterogeneous media[J].Acta Seismol.Sin.(in chinese),2000, 22(2) : 151~157.
    [19] Zhao A H, Zhang Z J, Teng J W.Minimum travel time tree algorithm for seismic ray tracing: improvement in efficiency[J].J.Geophys.Eng.2004(1) : 245~251.
    [20] Moser T J.Shortest path calculation of seimic rays[J].Geophysics, 1991, 56(1) : 59~67.
    [21] Vidale J.Finite-difference calculation of travel times[J].Bull.Seis.Soc.Am.,1988, 78(6) : 2062~2076.
    [22] 朱金明,王丽燕.地震波走时的有限差分算法[J].地球物理学报,1992,35(1) :86~92. Zhou J M, Wang L Y.Finite difference calculation of seismic travel times[J].Chinese J.Geophys.(Acta Geophysia Sinica) (in chinese), 1992, 35(1) : 86~92.
    [23] 赵改善,郝守玲,杨尔皓,等.基于旅行时线性插值的地震射线追踪算法[J].石油物探,1998,37(2) :14~24. Zhao G S, Hao S L, Yang E H, et al.Seismic ray tracing algorithm based on the linear traveltime interpolation[J].GPP, 1998, 37(2) : 14~24.
    [24] 刘洪,孟凡林,李幼铭.计算最小走时和射线路径的界面网全局方法[J].地球物理学报,1995,38(6) :823~832. Liu H, Meng F L, Li Y M.The interface grid method for seeking global minimum travel-time and the correspondent ray path[J].Chinese J.Geophys.(Acta Geophysia Sinica)(in chinese), 1995, 38(6) : 823~832.
    [25] 鲁彬,周立发,孔省吾,等.迭代优化的网络最短路径射线追踪方法研究[J].地球物理学进展,2009,24(4) :1420~1425. Lu B, Zhou L F, Kong S W, et al.A method of shortest path raytracing by iterative optimization[J].Progress in Geophys, (in chinese), 2009, 24(4) : 1420~1425.
    [26] 张东,谢宝莲,杨艳,等.一种改进的线性走时插值射线追踪算法[J].地球物理学报,2009,52(1) :200~205. Zhang D, Xie B L, Yang Y, et al.A ray tracing method based on improved linear traveltime interpolation[J].Chinese J.Geophys.(in Chinese), 2009, 52(1) : 200~205.
    [27] Sethian J A.A fast marching level set method for monotonically advancing fronts[J].Proc.Natl.Acad.Sci.,1996, 93: 1591~1595.
    [28] Sethian J A, Popovici A M.3-D traveltime computation using the fast marching method[J].Geophysics, 1999, 64 (2) :516~523.
    [29] Popovici A M, Sethian J A.3-D imaging using higher order fast marching traveltimes[J].Geophysics, 2002, 67 (2) :604~609.
    [30] Rawlinson N, Kennett B L N.Teleseismie tomography of the upper mantle beneath the southern Lachlan Orogen,Australia[J].Phys.Earth Planet.Inter., 2008, 167: 84~97.
    [31] Rawlinson N, Sambridge M.Seismic traveltime tomography of the crust and lithosphere[J].Advances in Geophysics,2003, 46: 81~197.
    [32] Paige C C, Saunders M A.LSQR: An algorithm for Sparse Linear Equations and Sparse Least Squares[J].ACM Trans.Math.Softw.,1982, 8(1) : 43~71.
    [33] Paige C C, Saunders M A.LSQR: Sparse Linear Equations and Least Squares Problems[J].ACM Trans.Math.Softw.,1982, 8(2) : 195~209.
    [34] 成谷,张宝金.反射地震走时层析成像中的大型稀疏矩阵压缩存储和求解[J].地球物理学进展,2008,23(3) :674~680. Cheng G, Zhang B J.Compression storage and solution of large sparse matrix in traveltime tomography of reflection seismic data[J].Progress in Geophysics, 2008, 23 (3) :674~680

版权所有:© 2023 中国地质图书馆 中国地质调查局地学文献中心