基于快速推进迎风双线性插值法的三维地震波走时计算
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
三维地震波走时计算技术是三维地震反演、层析成像、偏移成像等诸多地震数据处理技术中非常重要的正演计算工具.为了获得精度高且兼顾效率的三维走时计算方法:首先,在常规双线性插值公式推导过程中,充分利用平面波双线性假设的结论,获得了二元极小值超越方程的解析解,进而推导出了准确的局部走时计算公式,同时构造性地证明了该计算公式满足地震波的传播规律和Eikonal方程;其次,引入迎风差分的基本思想,提出迎风双线性插值的局部走时计算策略,该计算策略能简化算法、提高效率且保证无条件稳定性;然后,将上述计算公式和迎风双线性插值策略与常规快速推进法中的窄带技术结合,获得了一种新的基于快速推进迎风双线性插值法的三维地震波走时计算方法;最后,通过精度和效率分析检验了新算法的精度、效率和正确性,并通过计算实例验证了算法在面对复杂介质时的稳定性和有效性.
As a very important forward computational tool,3Dtraveltime computational scheme has been widely used in many 3Dseismic data processing techniques such as traveltime inversion,tomography,migration,etc.In the literatures,finite-difference scheme and conventional bilinearinterpolation method are frequently used for calculating the 3Dtraveltime.The finite-difference scheme has high-efficiency,but is not accurate enough.The conventional bilinear interpolation method has high-accuracy,but is not efficient enough.To build a method for computing the traveltime in 3Dspace with high-efficiency and high-accuracy,we firstly propose the analytic solutions of the binary minimum transcendental equation by making full use of the results of the plane wave bilinear assumption.Using the above analytic solutions and the deducing process of the conventional bilinear interpolation formulas, we derive some new local computational formulas of 3Dtraveltime.Secondly,we present a local computational strategy of 3Dtraveltime by introducing the upwind idea that was proposed in the fast marching method.The above new local computational formulas and strategy make up a new local computational method of 3D traveltime named upwind bilinear interpolation scheme.The new local computational scheme of3 Dtraveltime is more accuracy and efficient than the conventional bilinear interpolation scheme.We also prove that the new local scheme complies with the propagation law of seismic wave and the eikonal equation.Then,to calculate the traveltime distributed in the whole 3D model,we adapt the above local computational formulas and upwind bilinear interpolation strategy to the general narrow band technique of the fast marching method and build a global algorithm.This global algorithm is named fast marching upwind bilinear interpolation method.Finally,we analyze the accuracy and efficiency of the new method in the 3D homogeneous media and 3D layered media,and verify the stability and effectiveness of the new method with a modified 3D Marmousi model.Based on the above analysis and verification,we can make the following conclusions:1A constructive proof indicates that the new method proposed in this paper abides strictly by the Fermat principle and eikonal equation;2 The new method presented here is accurate and efficient than the 1st-order finite-difference method;3For introducing the upwind idea into the conventional bilinear interpolation formulas,the new method proposed in this paper is concise and unconditionally stable;4The new method presented here is stable and flexible in3 Dcomplex media.
引文
Asakawa E,Kawanaka T.1993.Seismic ray tracing using lineartraveltime interpolation.Geophysical Prospecting,41(1):99-111.
    Gao L,Li Y M,Chen X R,et al.2000.An attempt to seismic raytracing with symplectic algorithm.Chinese J.Geophys.(inChinese),43(3):402-410,doi:10.3321/j.issn:0001-5733.2000.03.014.
    Jing X L,Yang C C,Wang S Q.2007.A improved seismic reflectiontomographic method.Chinese J.Geophys.(in Chinese),50(6):1831-1836.
    Julian B R,Gubbins D.1977.Three-dimensional seismic raytracing.J.Geophys.,43:95-113.
    Li P M,Mei S Q,Ma Q B.2013.An improved bilinearinterpolation travel-time ray tracing method.Oil GeophysicalProspecting(in Chinese),48(4):553-558.
    Liu F,Zhang D,Yang Y,et al.2012.A fast numerical solution tominimum equation in 3-D seismic LTI ray tracing.Journal ofWuhan University(Natural Science Edition)(in Chinese),58(5):395-400.
    Mei S Q,Deng F,Zhong B S,et al.2010.The 3Dray tracingmethod base on the improved bilinear traveltime interpolation.Computing Techniques for Geophysical and GeochemicalExploration(in Chinese),32(2):152-157.
    Moser T J.1991.Shortest path calculation of seismic rays.Geophysics,56(1):59-67.
    Qin M Z,Chen J B.2000.Maslov asymptotic theory and symplecticalgorithm.Chinese J.Geophys.(in Chinese),43(4):522-533,doi:10.3321/j.issn:0001-5733.2000.04.013.
    Qu C,Zhou H L,Zhao D P.2007.Deep structure beneath the westmargin of Philippine Sea Plate and South China Sea from P andS wave travel time tomography.Chinese J.Geophys.(in Chinese),50(6):1757-1768.
    Sethian J A,Popovici A M.1999.3-D traveltime computation usingthe fast marching method.Geophysics,64(2):516-523.
    Sun J G.1998.On the limited aperture migration in twodimensions.Geophysics,63(3):984-994.
    Sun J G.1999.On the aperture effect in 3D Kirchhoff-typemigration.Geophysical Prospecting,47(6):1045-1076.
    Sun J G.2000.Limited-aperture migration.Geophysics,65(2):584-595.
    Sun J G.2004.True-amplitude weight functions in 3Dlimitedaperture migration revisited.Geophysics,69(4):1025-1036.
    Sun J G,Sun Z Q,Han F X.2011.A finite difference scheme forsolving the eikonal equation including surface topography.Geophysics,76(4):T53-T63.
    Vidale J E.1988.Finite-difference calculation of traveltimes.Bull.Seism.Soc.Am.,78(6):2062-2076.
    Vinje V,Iversen E,Gjystdal H.1993.Traveltime and amplitudeestimation using wavefront construction.Geophysics,58(8):1157-1166.
    Wang H Z,Ma Z T.1999.Traveltime calculation in 3Dmedia witharbitrary velocity distribution.∥69th Annual InternationalMeeting,SEG,Expanded Abstracts,1778-1781.
    Zhang D,Fu X R,Yang Y,et al.2009.3Dseismic ray tracingalgorithm based on LTI and partition of grid interface.ChineseJ.Geophys.(in Chinese),52(9):2370-2376,doi:10.3969/j.issn.0001-5733.2009.09.023.
    Zhang D,Zhang T T,Qiao Y F,et al.2013.A 3-D ray tracingmethod based on B-spline traveltime interpolation.Oil GeophysicalProspecting(in Chinese),48(4):559-566.
    Zhang M G,Jia Y G,Wang M Y,et al.2006.A global minimum
    traveltime ray tracing algorithm of wavefront expanding withinterface points as secondary sources.Chinese J.Geophys.(inChinese),49(4):1169-1175.
    高亮,李幼铭,陈旭荣等.2000.地震射线辛几何算法初探,地球物理学报,43(3):402-410,doi:10.3321/j.issn:0001-5733.2000.03.014.
    井西利,杨长春,王世清.2007.一种改进的地震反射层析成像方法.地球物理学报,50(6):1831-1836.
    李培明,梅胜全,马青坡.2013.一种改进的双线性插值射线追踪方法.石油地球物理勘探,48(4):553-558.
    刘锋,张东,杨艳等.2012.三维LTI射线追踪极小值方程的快速数值解法.武汉大学学报(理学版),58(5):395-400.
    梅胜全,邓飞,钟本善等.2010.基于改进的双线性走时插值的三维射线追踪.物探化探计算技术,32(2):152-157.
    秦孟兆,陈景波.2000.Maslov渐近理论与辛几何算法.地球物理学报,43(4):522-533,doi:10.3321/j.issn:0001-5733.2000.04.013.
    瞿辰,周蕙兰,赵大鹏.2007.使用纵波和横波走时层析成像研究菲律宾海板块西边缘带和南海地区的深部结构.地球物理学报,50(6):1757-1768.
    张东,傅相如,杨艳等.2009.基于LTI和网格界面剖分的三维地震射线追踪算法.地球物理学报,52(9):2370-2376,doi:10.3969/j.issn.0001-5733.2009.09.023.
    张东,张婷婷,乔友锋等.2013.三维走时场B样条插值射线追踪方法.石油地球物理勘探,48(4):559-566.
    张美根,贾豫葛,王妙月等.2006.界面二次源波前扩展法全局最小走时射线追踪技术.地球物理学报,49(4):1169-1175.

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