地震波传播速度的偏移分析.docxVIP

地震波传播速度的偏移分析 地震波的传播速度是地震数据处理中最基本、最重要的参数之一。地球物理工作者正在寻找最合适的速度模型。近年来,油气勘探目标已转向了地表和地下地质构造条件均十分复杂的地区,基于波动方程理论的叠前时间(或深度)偏移成像技术得到了快速的发展和应用,与此相应的以成像质量最优作为速度正确判断标准的偏移速度分析方法,如深度聚焦分析(DFA)、剩余曲率分析(RCA)和基于共聚焦(CFP)的速度分析方法等,已成为当前速度分析研究领域里的重点和热点,一些学者对此已经作了不同程度的归纳总结。本文主要对针对提高计算效率的叠前爆炸反射面模拟偏移速度分析方法、针对提高速度分析精度的相似谱判断标准非线性反演的全波动方程的速度模型建立、基于射线理论关系速度模型的时深转转以及针对消除浅层速度异常体影响等几种新方法技术的方法原理和应用效果进行了介绍。 1 叠前时间偏移 波场延拓偏移成像时,人们更愿意采用迭代波场延拓偏移速度分析的结果而不是传统的迭代Kirchhoff偏移结果。同样,若不考虑计算代价,人们倾向于采用波场层析算子进行速度场更新,而不是采用传统省时的射线层析方法,但是基于波动方程偏移与正演的各类方法,其计算代价严重制约了这些方法的实际推广和应用。 Biondi提出了一种叠前爆炸反射面模拟方法,大大缩减了基于波动方程偏移速度分析的计算时间。该方法既可应用于传统射线层析的偏移速度分析,也可用于波场延拓算子的偏移速度分析。它用叠前成像结果作为初始条件来正演单个地下偏移距域的共成像点道集(SODCIG)。虽然道集数目少于原始炮集数目,但却能包含速度分析所需的所有信息。 当地下反射界面水平时,采用如下成像条件 ΡSxiξ(t=0,x=xξ-hξ,z=zξ)=Ι(zξ,xiξ,hξ)ΡGxiξ(t=0,x=xξ+hξ,z=zξ)=Ι(zξ,xiξ,hξ)(1)PSxiξ(t=0,x=xξ?hξ,z=zξ)=I(zξ,xiξ,hξ)PGxiξ(t=0,x=xξ+hξ,z=zξ)=I(zξ,xiξ,hξ)(1) 由炮点波场(PSxiξ)和检波点波场(PGxiξ)来模拟虚拟炮数据(DSxiξ)和虚拟检波点数据(DGxiξ)。式中:I(zξ,xξ,hξ)为叠前成像结果;zξ为深度;xξ为水平坐标;hξ为地下水平偏移距;xiξiξ为SODCIG中相应的水平坐标。检波点波场沿时间t方向传播,而炮点波场沿时间反方向传播。 用下式 DSxiξ(t,x)=ΡSxiξ(t,x,z=0)DGxiξ(t,x)=ΡGxiξ(t,x,z=0)(2)DSxiξ(t,x)=PSxiξ(t,x,z=0)DGxiξ(t,x)=PGxiξ(t,x,z=0)(2) 提取所有时间和空间位置处在地表的波场。当地下反射界面非水平时,就得使震源函数沿着反射界面的法向方向,此时解释的偏移距是相对于地层倾角的,而非水平方向的。 图1是某合成数据由SODCIG转换来的角度域共成像点道集(ADCIG)。图1a是对合成数据用低10%的速度炮检点偏移的结果,图1b至图1e是用提出的方法从偏移数据集中获得的结果,其中Δxiξiξ分别为640,320,160,80m。除了图1d和图1e中假频稍严重之外,图1a和图1b都提供了非常有用的速度信息。 为避免速度迭代分析与拾取耗费大量的时间和精力,Douglas等提出利用叠前时间偏移道集中反射波的剩余正常时差直接校正原始均方根偏移速度函数的方法。Douglas等利用 e2+2e-v2(t2corrh-t20)/h2=0(3)e2+2e?v2(t2corrh?t20)/h2=0(3) 求出速度比例e后,再用 v=ve/(1+e)(4)v=ve/(1+e)(4) 直接求出校正后的速度。式中:ve为用于偏移或动校正的待修正的速度;h为偏移距;t0为中点处的双程旅行时;tcorrh为叠前时间偏移道集中与h对应的反射波至的旅行时间。 2 速度场的建立和应用 Soubaras等提出一种基于相似谱判断标准非线性反演的全波动方程速度模型建立方法。首先在粗网格上建立一个初始速度模型,然后调节炮偏移(modulated-shot migration)形成P道集,最后通过使相似谱 S(θ)=|∑xn(θ)|2Ν∑|xn(θ)|2(5)S(θ)=∣∣∑xn(θ)∣∣2N∑|xn(θ)|2(5) 最大化来线性反演速度场。式中:xn(θ)为与实际地震数据期望匹配的模型数据;N为参与计算的道数。 Sigsbee模型盐下速度估计难度较大。Soubaras等运用该模型对所研究的方法进行了测试,主要考察40 000～45 000 ft 和20 000～30 000 ft(1 ft≈0.304 8 m)深度范围内断层的成像效果,断层正好位于该模型顶部的复杂盐体之下。图2是用初始速度模型、Soub