复杂介质多波地震正演模拟的动力学射线追踪法.docxVIP

复杂介质多波地震正演模拟的动力学射线追踪法 1 地震射线追踪 辐射法是一种与光的性质和行为研究密切相关的几何地震方法，在地震学的发展过程中发挥着很好的作用。早在19世纪末,一些学者在地震学研究中,就开始用射线法计算各种球向对称结构介质中的地震波理论走时,以便将计算结果同观测资料进行比较,从而根据地震波的走时曲线研究地球的内部结构。经过一个多世纪的发展,射线追踪法目前已成为许多地球物理问题研究的基础,在诸如地震记录正演模拟、地震叠前偏移和地震层析成像等地球物理正、反演工作中皆得到广泛应用,并发挥着极其重要的作用。特别是在地震叠前偏移处理中,射线追踪法仍是目前计算速度最快和最有效的正演方法。 目前,动力学射线追踪在射线法中起着重要作用,为此,许多与此相关的方法得到了发展,而动力学射线追踪最重要的应用莫过于旁轴射线近似和高斯射线束的计算 2 zeopprits方程的基本原理 弹性波在传播过程中若遇到弹性突变的地层分界面,就象光在非均匀介质中传播的情形一样,能量在分界面上会重新分配,并产生反射和透射。而且,在反射和透射波中,还会出现不同于入射波类型的波,称之为转换波。图1展示了纵波入射于各向同性介质1与介质2构成的地层分界面上的情况,其中,除了产生反射纵波和透射纵波这两种同类波(P波)之外,还会产生反射横波和透射横波这两种转换波(P-SV波)。为计算反射波和透射波的振幅值,可如图1选取坐标轴,Z轴垂直向下,XOY平面与分界面重合。假定所研究的弹性波为平面谐波,并假设纵波位移沿射线传播方向为正,横波位移在射线传播方向的右侧为正(见图1射线上黑点处的箭头方向),那么,根据分界面上应力和位移连续的边界条件,即可推出计算反射纵波、透射纵波、反射横波和透射横波位移振幅系数的Zoeppritz方程为 式中:A、B、C、D分别为反射P波、反射P-SV波、透射P波、透射P-SV波的位移振幅系数,亦即各波分量振幅与入射波振幅的比值;θ 入射角和各波反射角及透射角的正弦与波传播速度之间的关系则由斯奈尔定律确定。 同样,当在垂直平面内偏振的横波入射到由两种介质构成的弹性分界面时,除了产生反射横波和透射横波这两种同类波(SV波)之外,还会产生反射纵波和透射纵波这两种转换波(SV-P波),如图2所示。采取同上一样的处理措施,可推得SV波入射时反射横波、透射横波、反射纵波和透射纵波位移振幅系数的计算公式为 式中:E、F、G、H分别为反射SV-P波、反射SV波、透射SV-P波、透射SV波的位移振幅系数;φ 根据(1)和(2)式,即可计算波在各向同性介质中传播、并以任意角度入射到弹性分界面上时的反射波和透射波能量的分配情况。不过,当入射角大于临界角时,反射波的性质会发生变化,这样就不能直接应用(1)、(2)式计算各反射波的振幅系数。在多波地震勘探中,由于反射转换波是以纵波速度下行传播,而以横波速度上行传播,因此,在多波正演动力学射线追踪过程中,可只考虑纵波超临界角入射情况。当纵波超临界角入射于分界面上时,根据斯奈尔定律可知cosθ 3 界面编码的基础 在应用斯奈尔定律进行射线追踪时,为了使射线能按正确的传播路径传播,就必须根据实际介质结构对界面进行正确的编码,即将各界面进行合理的排序,否则就会导致某些错误的追踪结果。因此,正确的界面编码是应用斯奈尔定律进行射线追踪的基础,也是保证获得正确追踪结果的前提。 图3为一结构比较简单的层状介质模型,对这样的结构,其界面编码不会存在困难,因为射线在旅行过程中所穿过的介质层序下行时依次递增,而上行时则依次递减,因此只要按常规方法从浅至深依次递增排列介质层序即可,图中R 此外,对于回转型和存在逆断层的界面,在界面编码中应把单一的物理界面划分成若干个序号递增、界面相连的局部界面(图7、8)。图7a为一回转型界面,图中G点能接收来自界面AB段和BC段两部分的反射波,为使射线追踪能处理这样的情况,就必须将一种介质人为地分成性质相同的两种介质,如图7b所示,将界面段AB编码为R 在射线追踪界面编码处理中,只要根据本文阐述的编码原则,运用上述几个具有代表性的模型界面编码规则,就可以处理复杂的介质结构。 4 连续介质介质层的简化 一般情况下,实际地下介质常可近似成均匀的层状介质,每层介质的速度都可看成在垂向或横向上是不变的。然而,在沉积旋回比较明显的地区,垂向上常由许多层速度渐变的薄层组成,另外,不少地区由于沉积物在横向上逐渐变化,导致介质层速度出现横向渐变。因此,有时把一些实际介质视为连续介质将可能更具有实际意义。那么,在射线追踪中,对于垂向速度连续变化的介质层,可以将其简化为若干个等厚的、速度均一的水平层进行追踪(图9a);而对于横向速度连续变化的介质层,则可以将其简化为若干等宽的、速度均一的垂直层进行追踪(图9b)。图9中的ΔZ、ΔX