各向异性TI介质qP反射波走时层析成像
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
地震走时层析成像是反演地层各向异性参数分布的有效方法,但是关于地震各向异性介质走时层析成像的研究并不多,其技术远远没有达到成熟的阶段.在野外数据采集时,地表反射波观测方式相对井间和垂直地震剖面观测方式的成本更低,利用qP反射波走时反演各向异性参数具有更加广泛的实用价值.本文实现的TI介质地震走时层析成像方法结合了TI介质反射波射线追踪算法、走时扰动方程和非线性共轭梯度算法,它可以对任意强度的TI介质模型进行反演,文中尝试利用qP反射波走时重建TI介质模型的参数图像.利用qP反射波对层状介质模型和块状异常体模型进行走时反演,由于qP波相速度对弹性模量参数和Thomsen参数的偏微分不同,所以可以分别反演弹性模量参数和Thomsen参数.数值模拟结果表明:利用qP反射波可以反演出TI介质模型的弹性模量参数与Thomsen参数,不同模型的走时迭代反演达到了较好的收敛效果,与各向同性介质走时反演结果相比较,各向异性介质走时反演结果具有较好的识别能力.
Seismic traveltime tomography is an effective method to invert underground anisotropic parameters.Research about the anisotropic traveltime tomographic method remains insufficient.Most studies on this issue are based on weak anisotropy assumption which was proposed by Leon Thomsen according to previous research.But the degree of anisotropic property of strata and rocks is uncertain in real applications.Therefore,a kind of traveltime tomography by use of qP reflection waves is proposed to invert arbitrary anisotropic TI media.First arrival traveltime tomography is often used for near-surface and crosswell anisotropic parameters inversion.In fact,traveltime tomography of qP reflection waves has a much moreextensive practical value,because the surface observation system is more economic and general than the systems of crossholes and vertical seismic profiles.Here,we present a new nonlinear traveltime inversion method for the surface observation system which combines several important features:(1)A robust reflected wave ray tracing method is used for arbitrary anisotropic TI media;(2)The first-order traveltime perturbation equation is not the eigenvector form of Cerveny′s linearized formula,which suffers from a singularity problem for the two quasi-shear waves;(3)An effective computation of the Jacobian matrix is employed for an arbitrary anisotropic TI media;and(4)It adopts a fast,local minimization search style of nonlinear inversion.These features make the traveltime tomography of qP reflection waves can be used to invert anisotropic TI media with arbitrary anisotropic degree.In numerical simulation section,this seismic reflection tomography method was used to invert for a layered model and a blocky abnormal body model,respectively.The elastic moduli parameter and Thomsen parameter of these anisotropic models were obtained,respectively.Firstly,we used the traveltimes of qP reflection waves to invert the elastic moduli parameters of the layered model and blocky abnormal body model.(1)According to the tomograms of the layered model,the velocity image is not reconstructed very well although the ray paths can cover the whole model perfectly.The velocity values of the first and third layers differ from the true velocity value greatly.And the bottom interface of the second layer is not recovered very well.Fortunately,the elastic moduli parameters are all reconstructed very well.We can see these images are all very close to sections of the true elastic moduli parameters.(2)After inversion of the blocky abnormal body model,the velocity image is not reconstructed very well.Because the velocity values of the blocky abnormal body differ from the true velocity values greatly.On the contrary,the elastic moduli parameters are very close to the profiles of true elastic moduli parameters.But the c11 parameter is not as good as the other three parameters.Secondly,we used the traveltimes of qP reflection waves to invert the Thomsen parameters of the layered model and blocky abnormal body model.(1)According to the tomograms of the layered model,the velocity image is still not reconstructed very well,because the velocity values of the second layer differ from the true velocity values greatly.The Thomsen parameters are reconstructed very well,except the image ofβ0is slightly not so good.We can see the inverted and true images ofα0,εandδ*are very close to each other.(2)After inversion of the blocky abnormal body model,the velocity image is not reconstructed correctly,because its numerical value and shape differ from the true velocity model greatly.However,the profiles ofα0,β0andδ*are reconstructed successfully,except theεparameter.The images of these three inverted parameters are very close to that of the true parameters.Besides,these four numerical experiments all reach satisfactory convergence levels for different anisotropic models. Hence,the simulation results show that traveltime of qP reflection waves can invert anisotropic parameters correctly.The traveltimes of qP reflection waves have been used to invert anisotropic parameters for the layered model and blocky abnormal body model successfully.The images of the reconstructed results are very close to the images of the true anisotropic parameters.The tomographic method proposed here has superior capability in recovering theoretical models compared with the isotropic traveltime tomographic method.The more information is used,the less uncertainty of the inversion results according to the inverse theory.In other words,the inversion results can be improved greatly if the traveltimes of qP,qSV and qSH waves are used in the traveltime inversion.Therefore,joint traveltime inversion of qP,qSV and qSH waves should worth further research.
引文
Cerveny V.2001.Seismic Ray Theory.Cambridge:CambridgeUniversity Press.
    Chapman C H,Pratt R G.1992.Traveltime tomography inanisotropic media—I Theory.Geophysical Journal International,109(1):1-19.
    Crampin S.1984.Effective anisotropic elastic constants for wavepropagation through cracked solids.Geophysical Journal of theRoyal Astronomical Society,76(1):135-145.
    Helbig K.1981.Systematic classification of layered-inducedtransverse isotropy.Geophysical Prospecting,29(4):550-577.
    Jech J.1988.Three-dimensional inversion problem for inhomogeneoustransversely isotropic media.Studia Geophysica et Geodetica,32(2):136-143.
    Kreyszig E.1993.Advanced Engineering Mathematics.New York:John Wiley&Sons Inc,790.
    Li F,Cao S Y,Yao J.2012.Calculation of phase and groupvelocities in an arbitrary anisotropic medium.Chinese Journalof Geophysics(in Chinese),55(10):3420-3426.
    Liu Y Z,Wang G Y,Dong L G,et al.2014.Joint inversion of VTIparameters using nonlinear traveltime tomography.ChineseJournal of Geophysics(in Chinese),57(10):3402-3410.
    Liu Y Z,Yang J Z.2014.Offset-weighted seismic tomography aimedat using the first arrival efficiently.Geophysical Prospecting forPetroleum(in Chinese),53(1):99-115.
    Lu M H,Tang J H,Yang H Z,et al.2005.P-wave traveltimeanalysis and Thomsen parameters inversion in orthorhombicmedia.Chinese Journal of Geophysics(in Chinese),48(5):1167-1171.
    Pratt R G,Chapman C H.1992.Traveltime tomography in anisotropicmedia—II Application.GeophysicalJ ournal International,109(1):20-37.
    Tryggvason A,Schmelzbach C,Juhlin C.2009.Traveltimetomographic inversion with simultaneous static corrections—well worth the effort.Geophysics,74(6):WCB25-WCB33.
    Wang X X,Tsvankin I.2013.Ray-based gridded tomography fortilted transversely isotropic media.Geophysics,78(1):C11-C23.
    Wang Y H.2014.Seismic ray tracing in anisotropic media:Amodified Newton algorithm for solving highly nonlinear systems.Geophysics,79(1):T1-T7.
    Zelt C A,Azaria A,Levander A.2006.3D seismic refractiontraveltime tomography at a groundwater contamination site.Geophysics,71(5):H67-H78.
    Zhang W S,He Q D.1999.Statistical tomography in anisotropicmedia by using crosshole traveltimes of qP wave and qSV wave.Chinese Journal of Geophysics(in Chinese),42(S1):195-203.
    Zhou B,Greenhalgh S A,Sinadinovski C.1992.Iterativealgorithms for the damped minimum norm,least squares andconstraining problem in seismic tomography.ExplorationGeophysics,23:497-505.
    Zhou B,Greenhalgh S A.2005.Analytic expressions for thevelocity sensitivity to the elastic moduli for the most generalanisotropic media.Geophysical Prospecting,53(4):619-641.
    Zhou B,Greenhalgh S A.2006.Raypath and traveltime computationsfor 2Dtransversely isotropic media with dipping symmetry axes.Exploration Geophysics,37(2):150-159.
    Zhou B,Greenhalgh S A,Green A.2008.Nonlinear traveltimeinversion scheme for crosshole seismic tomography in tiltedtransversely isotropic media.Geophysics,73(4):D17-D33.
    Zhou H,He Q D.1995.Seismic traveltime inversion for transversely
    isotropic media.Geophysical Prospecting for Petroleum(inChinese),34(3):63-68.
    Zhou H W.2006.Multiscale deformable-layer tomography.Geophysics,71(3):R11-R19.
    李芳,曹思远,姚健.2012.任意各向异性介质相(群)速度的计算.地球物理学报,55(10):3420-3426,doi:10.6038/j.issn.0001-5733.2012.10.025.
    刘玉柱,王光银,董良国等.2014.VTI介质多参数联合走时层析成像方法.地球物理学报,57(10):3402-3410.
    刘玉柱,杨积忠.2014.有效利用初至信息的偏移距加权地震层析成像方法.石油物探,53(1):99-115.
    卢明辉,唐建侯,杨慧珠等.2005.正交各向异性介质P波走时分析及Thomsen参数反演.地球物理学报,48(5):1167-1171,doi:10.3321/j.issn:0001-5733.2005.05.026.
    张文生,何樵登.1999.各向异性介质中利用跨孔qP波和qSV波走时的统计法成像.地球物理学报,42(S1):195-203.
    周辉,何樵登.1995.横向各向同性介质走时反演.石油物探,34(3):63-68.

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