04-偏移算法的实现和结果分析.docVIP

第4章 偏移算法的实现和结果分析 通过前面内容的介绍，我们了解了波动方程法偏移的基本原理和复杂地表条件下的叠前深度偏移方法的基本原理。本章主要通过理论模型来检验的正确性和有效性，比较各种方法的成像效果，并对其成像结果进行分析。这里用到的是一个二维均匀介质点绕射的例子。模型如图4-1所示。模型设计的目的是检验输出剖面的脉冲响应。模型大小为200*200m2,绕射点位于中心处的（100，100）。 图4-1 均匀介质单点绕射模型 介质速度为2000m/s，网格大小为10m*10m,正演所用的子波为雷克子波，子波的主频为35Hz。 图4-2 地面接收的点绕射双曲线 图3-4 点绕射逆时偏移结果 图4-2为其地面接收的绕射双曲线，图4-3为逆时偏移得到的成像结果。可见绕射曲线收敛很好,能量集中在绕射点附近。 虽然点绕射体似乎在实际地震勘探中不如反射界面那么常见，但是由于可以把反射界面或绕射体当作许多这样的点的集合，所以点绕射体的偏移具有一般性的意义。它的成像证明了这种逆时计算的合理性和可行性。 4.1.2 逆冲断层成像 这里选用了一个逆冲断层的例子。速度模型如图4-4所示，大小为800*500，网格大小为2*2m。速度参数大小为由上至下：1500m/s，1800m/s, 2200m/s, 3000m/s, 3500m/s, 4000m/s。模型中有断层，背斜，向斜等地质构造。 图4-4 逆冲断层模型 在合成的叠加记录（图4-5）上可以看到断点绕射波，斜界面的同相轴，以及多次波的存在，在逆时偏移剖面（图4-6）中可以看到同相轴位置与速度模型一致，多次波影响消除，但由于是隔一道采样做的偏移处理，由于网隔的加大，存在频散，导致分辨率降低。从上述数值模拟例子可以知道此方法没有倾角限制，对陡倾角和复杂构造同样可以清晰地成像。 图4-5 逆冲断层的合成叠加剖面 图4-6 任意差分精细积分结果(隔一道采集做的偏移) 4.1.3 复杂构造成像 这里以如图4-7所示的模型分别用Kirchhoff深度偏移、相移法和有限差分深度偏移方法与逆时偏移对同一叠加剖面（图4-8）成像，比较了这些偏移方法的效果。 图4-7 模型示意 模型由一个断层和一个背斜构造构成，大小为1200*400m2，网格大小为1*1m，速度由上至下分别为1700m/s，2500m/s，3000m/s。图4-8为其对应的叠加剖面。 图4-8 模型对应的叠加剖面（共1200道） 图4-9为Kirchhoff深度偏移的结果，图4-10为相移法偏移结果转换到深度域的剖面，图4-11为显式有限差分深度偏移的结果。图4-12为逆时偏移的结果。 图4-9 Kirchhoff深度偏移剖面 在图4-9上可以看到主要的同相轴都成像很好，在速度发生横向变化的地方陡倾角断面以及下伏的背斜平界面都能清晰准确的成像。由于Kirchhoff自身方法的影响，在偏移剖面中可以看到有背景噪音存在。 图4-10 相移法偏移剖面 图4-10是相移法偏移结果经时深转换变换到深度域的结果，以便于和其他方法对比。由于相移法属于时间偏移范围，不能很好地处理速度的横向变化问题。因此可以看到在上图中由于上伏断层引起的速度横向变化的影响使得下面的平界面不能很好得成像。与Kirchhoff深度偏移相比成像效果要差很多。与图4-6所示的显式有限差分深度偏移和任意差分精细积分逆时偏移相比也差很多。 图4-11 显式有限差分深度偏移剖面 在图4-11显式有限差分深度偏移的偏移剖面上可以看到偏移成像效果很好，只是边界有一些残余能量。陡倾角的断面成像很好，下伏地层的背斜和平界面也没有发生畸变。 图4-12 逆时偏移剖面 图4-12是逆时偏移的结果，可以看到同相轴归位很好，波形也没有发生畸变，无论是陡倾角还是下伏地层均成像很好。通过与上述三种偏移方法对应剖面的对比可以看出逆时偏移效果很好。在模型速度准确已知的情况下成像最准确，70°显式有限差分偏移的效果与之最接近，Kirchhoff深度偏移次之。相移法最差，值得一提的是逆时偏移计算算法复杂，计算量大所以计算也最费时。 4.2 marmousi模型 marmousi模型，它是一个比较复杂的2-D地质模型。模型上部主要有三大断裂组成；中部有一盐丘；深部两侧有高速体，中间在2400m以下有一低速目的层。其速度场如下图所示： 图4-13 marmousi模型速度场 所用炮记录来自于SEG/EAGE提供的声波模拟记录。该模型记录有240炮，每炮96道，750个采样点，4ms采样，最小偏移距为200m，道间距为25m；速度场为497*750，CDP间距为12.5m，深度采样为4m。单炮记录如图4-所示。