针对复杂地形的三种地震波走时算法及对比
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
复杂地形条件下地震波走时算法对于研究复杂地形地区的成像问题有着重要的意义.为了得到精度高且适应于复杂地形的走时算法,首先提出阶梯网格迎风差分法.然后将该方法与不等距网格有限差分法和混合网格线性插值法进行对比研究,得出如下结论:混合网格线性插值法的计算精度最高,但其计算效率最低;阶梯网格迎风差分法的计算精度最低,但其计算效率最高;不等距网格有限差分法的计算精度和计算效率均居中;而究竟选取哪种算法作为给定复杂地形模型的地震波走时算法,应该综合考虑地形的特点、所研究问题对计算精度及计算效率的要求等因素.最后通过一个计算实例验证了三种算法在面对复杂地形、近地表及地下复杂介质等复杂地质条件时均有很好的适应性和稳定性.
The traveltime computation scheme in complex topographical conditions has a great significance to seismic imaging in complex topographical regions.To obtain a scheme that can treat the complex topography with high precision and good flexibility,we firstly present a upwind finite difference scheme with stepwise grid.Then we compare this scheme with the following two schemes that we have published: finite difference scheme with nonuniform grid spacing and linear interpolation scheme with hybrid grid,and we can make the following conclusions: the linear interpolation scheme with hybrid grid has the best computational accuracy,but the worst computational efficiency;the upwind finite difference scheme with setpwise grid has the worst computational accuracy,but the best computational efficiency;the finite difference scheme with nonuniform grid spacing has the moderate computational accuracy and efficiency.We should choose the scheme by synthetically considering the characteristic of the topography,the computational accuracy and efficiency requirements of our object,and the other factors.At last,a numerical test shows that the three schemes all can treat the complex topographical region problem with satisfactory stability and flexibility.
引文
[1]孙建国.复杂地表条件下地球物理场数值模拟方法评述.世界地质,2007,26(3):345-362.Sun J G.Methods for numerical modeling of geophysicalfields under complex topographical conditions:a criticalreview.Global Geology(in Chinese),2007,26(3):345-362.
    [2]Sun J G.True-amplitude weight functions in 3Dlimited-aperture migration revisited.Geophysics,2004,69(4):1025-1036.
    [3]Sun J G.Limited-aperture migration.Geophysics,2000,65(2):584-595.
    [4]Sun J G.On the aperture effect in 3D Kirchhoff-typemigration.Geophysical Prospecting,1999,47(6):1045-1076.
    [5]Sun J G.On the limited aperture migration in twodimensions.Geophysics,1998,63(3):984-994.
    [6]井西利,杨长春,王世清.一种改进的地震反射层析成像方法.地球物理学报,2007,50(6):1831-1836.Jing X L,Yang C C,Wang S Q.A improved seismicreflection tomographic method.Chinese J.Geophys.(inChinese),2007,50(6):1831-1836.
    [7]瞿辰,周蕙兰,赵大鹏.使用纵波和横波走时层析成像研究菲律宾海板块西边缘带和南海地区的深部结构.地球物理学报,2007,50(6):1757-1768.Qu C,Zhou H L,Zhao D P.Deep structure beneath the westmargin of Philippine Sea Plate and South China Sea from Pand S wave travel time tomography.Chinese J.Geophys.(inChinese),2007,50(6):1757-1768.
    [8]Julian B R,Gubbins D.Three-dimensional seismic raytracing.J.Geophys.,1977,43(1-2):95-113.
    [9]徐涛,徐果明,高尔根等.三维复杂介质的块状建模和试射射线追踪.地球物理学报,2004,47(6):1118-1126.Xu T,Xu G M,Gao E G,et al.Block modeling and shootingray tracing in complex 3-D media.Chinese J.Geophys.(inChinese),2004,47(6):1118-1126.
    [10]Vinje V,Iversen E,Gjystdal H.Traveltime and amplitudeestimation using wavefront construction.Geophysics,1993,58(8):1157-1166.
    [11]Han F X,Sun J G,Sun Z Q,et al.Positioning of grid pointsin wave front construction.Applied Geophysics,2009,6(3):248-258.
    [12]Vidale J E.Finite-difference calculation of travel times.Bull.Seism.Soc.Am.,1988,78(6):2062-2076.
    [13]Sethian J A,Popovici A M.3-D travel time computationusing the fast marching method.Geophysics,1999,64(2):516-523.
    [14]Rawlinson N,Sambridge M.Wave front evolution instrongly heterogeneous layered media using the fast marchingmethod.Geophys.J.Int.,2004,156(3):631-647.
    [15]Sun J G,Yang H,Han F X.A finite-difference scheme forsolving the eikonal equation with varying grid spacing.77thAnn.Internat.Mtg.,Soc.Expl.Geophys.,ExpandedAbstracts,2007:2120-2124.
    [16]Schneider W A Jr,Ranzinger K A,Balch A H,et al.Adynamic programming approach to first arrival traveltimecomputation in media with arbitrarily distributed velocities.Geophysics,1992,57(1):39-50.
    [17]Asakawa E,Kawanaka T.Seismic ray tracing using lineartraveltime interpolation.Geophysical Prospecting,1993,41(1):99-111.
    [18]Wang H Z,Ma Z T,Hubral P.Traveltime calculation in 3Dmedia with arbitrary velocity distribution.69th Ann.Internat.Mtg.,Soc.Expl.Geophys.,Expanded Abstracts1999:1778-1781.
    [19]Moser T J.Shortest path calculation of seismic rays.Geophysics,1991,56(1):59-67.
    [20]Cao S H,Greenhalgh S.Calculation of the seismic first-breaktime field and its ray path distribution using a minimumtraveltime tree algorithm.Geophys.J.Int.,1993,114(3):593-600.
    [21]张美根,贾豫葛,王妙月等.界面二次源波前扩展法全局最小走时射线追踪技术.地球物理学报,2006,49(4):1169-1175.Zhang M G,Jia Y G,Wang M Y,et al.A global minimumtraveltime raytracing algorithm of wavefront expanding withinterface points as secondary sources.Chinese J.Geophys.(in Chinese),2006,49(4):1169-1175.
    [22]李满树,方伍宝,周腾等.初至波走时信息在复杂地区近地表速度反演中的应用.石油物探,2004,43(1):72-75.Li M S,Fang W B,Zhou T,et al.Inversion of near-surfacevelocity through first-arrival traveltime in complex region.Oil Geophysical Prospecting(in Chinese),2004,43(1):72-75.
    [23]刘少勇,王华忠,张兵.起伏地表Kirchhoff积分法叠前深度偏移方法研究与应用.2010年国际石油地球物理技术交流会会议专刊,科学出版社,2010:49-54.Liu S Y,Wang H Z,Zhang B.Kirchhoff integral PSDM forrugged topography:technology and application(in Chinese).Conference on Exploration Geophysics in the West of China,Science Press,2010:49-54.
    [24]刘国峰,刘洪,李博等.起伏地表直接叠前时间偏移.石油地球物理勘探,2010,45(2):196-200.Liu G F,Liu H,Li B,et al.Direct pre-stack time migrationon rugged topography.Oil Geophysical Prospecting(inChinese),2010,45(2):196-200.
    [25]张建中,陈世军,徐初伟.动态网络最短路径射线追踪.地球物理学报,2004,47(5):899-904.Zhang J Z,Chen S J,Xu C W.A method of shortest pathraytracing with dynamic networks.Chinese J.Geophys.(inChinese),2004,47(5):899-904.
    [26]张东,谢宝莲,杨艳等.一种改进的线性走时插值射线追踪算法.地球物理学报,2009,52(1):200-205.Zhang D,Xie B L,Yang Y,et al.A ray tracing methodbased on improved linear traveltime interpolation.Chinese J.Geophys.(in Chinese),2009,52(1):200-205.
    [27]蒋先艺.基于二维与三维复杂结构模型正演的地震数据采集设计方法研究[博士论文].成都:成都理工大学地球物理学院,2004.Jiang X Y.The Study of Seismic Acquisition Method Basedon Modeling of 2Dand 3DComplex Geologic Structure[Ph.D.thesis](in Chinese).Chengdu:College of Geophysics ofChengdu University of Technology,2004.
    [28]赵瑞,白超英.复杂层状模型中多次波快速追踪———一种基于非规则网格的最短路径算法.地震学报,2010,32(4):433-444.Zhao R,Bai C Y.Fast multiple ray tracing within complexlayered media:the shortest path method based on irregulargrid cells.Acta Seismologica Sinica(in Chinese),2010,32(4):433-444.
    [29]黄国娇,白超英.二维复杂层状介质中地震多波走时联合反演成像.地球物理学报,2010,53(12):2972-2981.Huang G J,Bai C Y.Simultaneous inversion with multipletraveltimes within 2-D complex layered media.Chinese J.Geophys.(in Chinese),2010,53(12):2972-2981.
    [30]白超英,黄国娇,李忠生.三维复杂层状介质中多震相走时联合反演成像.地球物理学报,2011,54(1):182-192.Bai C Y,Huang G J,Li Z S.Simultaneous inversioncombining multiple-phase travel times within 3D complexlayered media.Chinese J.Geophys.(in Chinese),2011,54(1):182-192.
    [31]Sun J G,Sun Z Q,Han F X.A finite difference scheme forsolving the eikonal equation including surface topography.Geophysics,2011,76(4):T53-T63
    [32]Sun J G,Sun Z Q,Han F X.A finite difference scheme forsolving the eikonal equation including surface topography.79th Ann.Internat.Mtg.,Soc.Expl.Geophys.,ExpandedAbstracts,2009:2667-2671.
    [33]Sun J G,Sun Z Q,Han F X.A method for computing the 2Dseismic traveltime including the surface topography bycombining the linear interpolation and the narrow bandtechnique.CPS/SEG 2009International Geophysical Conference&Exposition,2009:ID-1187.
    [34]孙章庆,孙建国,韩复兴.复杂地表条件下基于线性插值和窄带技术的地震波走时计算.地球物理学报,2009,52(11):2846-2853.Sun Z Q,Sun J G,Han F X.Traveltimes computation usinglinear interpolation and narrow band technique under complextopographical conditions.Chinese J.Geophys.(in Chinese),2009,52(11):2846-2853.

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