基于LTI和网格界面剖分的三维地震射线追踪算法
|
摘要
将二维线性走时插值射线追踪算法(LTI)推广应用至三维模型,并结合网格界面剖分方式,提出了一种三维射线追踪算法.该算法既可获得高精度的全局最小走时和射线路径,又具有快速稳定的特点.三维模型计算结果表明,在模型参数包括网格密度完全相同情况下,本文算法较传统的三维最短路径算法在计算效率、走时和射线的计算精度上均有明显改进.
In the paper a new three-dimensional seismic ray tracing method is put forward, which is derived from the improvement of two-dimensional seismic ray tracing algorithm based on linear travel-time interpolation and partition of grid interface. This method inherits the advantage of not only the high accuracy of the minimum travel-times and ray paths but also the high calculation speed associated with partition of grid interface. The calculating results on some threedimensional models demonstrate that, compared to the traditional three-dimensional shortest path ray tracing algorithm, the speed of this new three-dimensional seismic ray tracing calculation and the accuracy of the minimum travel-times and ray paths are greatly improved.
In the paper a new three-dimensional seismic ray tracing method is put forward, which is derived from the improvement of two-dimensional seismic ray tracing algorithm based on linear travel-time interpolation and partition of grid interface. This method inherits the advantage of not only the high accuracy of the minimum travel-times and ray paths but also the high calculation speed associated with partition of grid interface. The calculating results on some threedimensional models demonstrate that, compared to the traditional three-dimensional shortest path ray tracing algorithm, the speed of this new three-dimensional seismic ray tracing calculation and the accuracy of the minimum travel-times and ray paths are greatly improved.
引文
[1] 黄中玉.从2007年度EAGE年会看地球物理技术的发展.勘探地球物理进展,2007,30(6) :473-485 Huang Z Y. Developing trends of geophysical technologies:Highlights of 2007 EAGE conference. Progress in Exploration Geophysics (in Chinese), 2007, 30(6) : 473-485
[2] 王辉,常旭.基于图形结构的三维射线追踪方法.地球物理学报,2000,43(4) :534-541 Wang H, Chang X. 3-D ray tracing method based on graphic structure. Chinese J . Geophys. (in Chinese) , 2000 , 43(4) :534-541
[3] Vidale J E. Finite-difference traveltime calculation. Bull. ,Seis. Sot. Am. , 1988, 78(6) : 2062-2076
[4] Vidale J E. Finite-difference calculation of traveltimes in three dimensions. Geophysics, 1990, 55(5) : 521-526
[5] Schneider W A Jr. Robust and efficient upwind finite-difference traveltime calculations in three dimensions.Geophysics, 1995, 60(4) : 1108-1117
[6] Vinje V, Iversen E, Gjoystdal H. Traveltime and amplitude estimation using wavefront construction. Geophysics, 1993,58(8) : 1157-1166
[7] Asakawa Eiichi, Kawanaka Taku. Seismic ray tracing using linear traveltime interpolation. Geophysical Prospecting, 1993, 41(1) : 99-112
[8] Nakanishi I, Yamaguchi K. A numerical experiment on nonlinear image reconstruction from first-arrival times for two-dimensional island arc structure. J. Phys. Earth, 1986, 34(2) : 195-201
[9] Moser T J. Shortest path calculation of seismic rays.Geophysics, 1991, 56(1) : 59-67
[10] Shunhua Cao, Stewart Greenhalgh. Calculation of the seismic first-break time field and its ray path distribution using aminimum traveltime tree algorithm. Geophys. J. Int, 1993,114(3) :593-600
[11] Ludek Klimes and Michal Kvasnicka. 3-D network ray tracing. Geophys. J. Int, 1994, 116(3) :726-738
[12] Bube K P, Washbourne J K. Wavetracing: Ray tracing for the propagation of band-limited signals: Part 1-theory.Expanded Abstract of the 71st Annual International SEG Meeting, 2001． 1131-1134
[13] 张东,谢宝莲,杨艳等.一种改进的线性走时插值射线追踪算法.地球物理学报,2009,52(1) :200-205 Zhang D, Xie B L, Yang Y, et al. A ray tracing method based on improved linear traveltime interpolation. Chinese J,Geophys. (in Chinese) , 2009, 52(1) :200-205
[14] Bai C Y, Greenhalgh S, Zhou B. 3D ray tracing using a modified shortest-path method. Geophysics, 2007,72(4) :T27-T36
[15] 张美根,贾豫葛,王妙月等.界面二次源波前扩展法全局最小走时射线追踪技术.地球物理学报,2006,49(4) :1169-1175 Zhang M G, Jia Y G, Wang M Y, et al. A Global Minimum Traveltime Raytracing Algorithm of Wavefront Expanding with Interface Points as Secondary Sources. Chinese J.Geophys. (in Chinese) , 2006, 49(4) : 1169-1175
[16] 张美根,程冰洁,李小凡等.一种最短路径射线追踪的快速算法.地球物理学报,2006,49(5) :1467-1474 Zhang M G, Jia Y G, Wang M Y, et al. A fast algorithm of shortest path ray tracing. Chinese J. Geophys. (in Chinese),2006, 49(5) : 1467-1474
[17] 张建中,陈世军,余大洋.最短路径射线追踪方法及其改进.地球物理学进展,2003,18(1) :146-150 Zhang J Z, Chen S J , YU D Y. Improvement of shortest path ray tracing method. Progress in Geophysics (in Chinese),2003, 18(1) : 146-150
[18] William H Press, Brian P Flannery, Saul A Teukolsky, et al. Numerical Recipes in C:The Art of Scientific Computing (Second Edition). New York: Syndicate of the University of Cambridge, 1992
[2] 王辉,常旭.基于图形结构的三维射线追踪方法.地球物理学报,2000,43(4) :534-541 Wang H, Chang X. 3-D ray tracing method based on graphic structure. Chinese J . Geophys. (in Chinese) , 2000 , 43(4) :534-541
[3] Vidale J E. Finite-difference traveltime calculation. Bull. ,Seis. Sot. Am. , 1988, 78(6) : 2062-2076
[4] Vidale J E. Finite-difference calculation of traveltimes in three dimensions. Geophysics, 1990, 55(5) : 521-526
[5] Schneider W A Jr. Robust and efficient upwind finite-difference traveltime calculations in three dimensions.Geophysics, 1995, 60(4) : 1108-1117
[6] Vinje V, Iversen E, Gjoystdal H. Traveltime and amplitude estimation using wavefront construction. Geophysics, 1993,58(8) : 1157-1166
[7] Asakawa Eiichi, Kawanaka Taku. Seismic ray tracing using linear traveltime interpolation. Geophysical Prospecting, 1993, 41(1) : 99-112
[8] Nakanishi I, Yamaguchi K. A numerical experiment on nonlinear image reconstruction from first-arrival times for two-dimensional island arc structure. J. Phys. Earth, 1986, 34(2) : 195-201
[9] Moser T J. Shortest path calculation of seismic rays.Geophysics, 1991, 56(1) : 59-67
[10] Shunhua Cao, Stewart Greenhalgh. Calculation of the seismic first-break time field and its ray path distribution using aminimum traveltime tree algorithm. Geophys. J. Int, 1993,114(3) :593-600
[11] Ludek Klimes and Michal Kvasnicka. 3-D network ray tracing. Geophys. J. Int, 1994, 116(3) :726-738
[12] Bube K P, Washbourne J K. Wavetracing: Ray tracing for the propagation of band-limited signals: Part 1-theory.Expanded Abstract of the 71st Annual International SEG Meeting, 2001． 1131-1134
[13] 张东,谢宝莲,杨艳等.一种改进的线性走时插值射线追踪算法.地球物理学报,2009,52(1) :200-205 Zhang D, Xie B L, Yang Y, et al. A ray tracing method based on improved linear traveltime interpolation. Chinese J,Geophys. (in Chinese) , 2009, 52(1) :200-205
[14] Bai C Y, Greenhalgh S, Zhou B. 3D ray tracing using a modified shortest-path method. Geophysics, 2007,72(4) :T27-T36
[15] 张美根,贾豫葛,王妙月等.界面二次源波前扩展法全局最小走时射线追踪技术.地球物理学报,2006,49(4) :1169-1175 Zhang M G, Jia Y G, Wang M Y, et al. A Global Minimum Traveltime Raytracing Algorithm of Wavefront Expanding with Interface Points as Secondary Sources. Chinese J.Geophys. (in Chinese) , 2006, 49(4) : 1169-1175
[16] 张美根,程冰洁,李小凡等.一种最短路径射线追踪的快速算法.地球物理学报,2006,49(5) :1467-1474 Zhang M G, Jia Y G, Wang M Y, et al. A fast algorithm of shortest path ray tracing. Chinese J. Geophys. (in Chinese),2006, 49(5) : 1467-1474
[17] 张建中,陈世军,余大洋.最短路径射线追踪方法及其改进.地球物理学进展,2003,18(1) :146-150 Zhang J Z, Chen S J , YU D Y. Improvement of shortest path ray tracing method. Progress in Geophysics (in Chinese),2003, 18(1) : 146-150
[18] William H Press, Brian P Flannery, Saul A Teukolsky, et al. Numerical Recipes in C:The Art of Scientific Computing (Second Edition). New York: Syndicate of the University of Cambridge, 1992
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心