井间地震弹性波逆时偏移成像的数值模拟.docxVIP

井间地震弹性波逆时偏移成像的数值模拟 0 井间地震资料处理 简言之，井间地震是在两口井内配置震源（激发信号）和探测器（接收信号）的方法，并通过探测器获得的信号进行分析和处理，在两口井之间的地质结构中获得地震勘探方法。井间地震的资料处理主要分为两类：一类是层析成像，主要利用接收信号中直达波的到时来得到井间的速度模型；另一类是反射波成像，它是在层析成像所得的速度场基础上，将接收信号中的反射波归位到地层的正确位置。与层析成像相比，反射波成像的垂直向分辨率更高、勘探区域更广，是井间地震资料处理中不可缺少的环节[1-2]。 井间地震反射波成像的算法可以分为两类。一类是利用波的几何射线理论，选用不同的道集并根据其特点进行成像的方法。其代表性的算法有XSP-CDP转化法 1 井间地震波的模拟处理 井间地震弹性波的波场模拟是使用计算效率高的时域有限差分法(Finite-difference timedomain,FDTD) 在图1模型中，将模型的左上角设为坐标原点，向右为x轴的正方向，表示水平的距离；向下为z轴的正方向，表示所处的深度。模型大小为300 m×600 m，网格大小为?x=?z=0.5 m，时间间隔约为57μs，最大的记录时间是0.4 s。震源坐标在(0 m,300 m)处，如图1中左侧星号所对应的位置。震源使用200 Hz的雷克子波，加载到V 本文对数值模拟获得的人工合成信号进行处理和反射波偏移。首先是对接收到的信号进行波场分离。由于直达波的幅度较大，对界面成像起干扰作用，需要被去除。常用的方法是中值滤波，先将需要去除掉的直达波信号进行排齐，再对其进行中值滤波。用排齐后的信号减去中值滤波的信号后再解排，就可得到去除直达波后的波列，如图3。 获得处理过的井间地震信号后，即可进行逆时偏移。逆时偏移是利用了波的时间一致性原理。当波传播到地层界面时会发生反射，在反射界面同时同点处产生入射波(上行波)和反射波(下行波)。逆时偏移就是根据这样的原理来对地层界面进行成像的。它的实现通常分为三步：一是震源波场的正向计算，来模拟地震当中的下行波场；二是检波器波场的逆向延拓，来模拟地震当中的上行波场；三是选用合适的成像条件，将上行和下行波场两者结合起来，来成像地层界面。其具体的实现方法可以参考文献 逆时偏移最常用到互相关成像条件，它的计算公式如下： 式(1)中，S 2 层模型试验结果 逆时偏移的成像条件有多种，如激发时间成像条件、反褶积成像条件等。不同的成像条件，它的实现方式和成像效果是不同的 可以看到，归一化的成像条件就是在常规互相关成像的基础上用震源波场或检波器波场来对其进行归一化。为了防止分母为零的情况，也可以在分母处添加一个小数。由于互相关成像条件需要对波场进行积分，就使得其得到的成像值物理意义变得模糊，没有合适的量纲与之匹配。而使用不同的波场对互相关条件成像值进行归一化后，得到的结果就具有了量纲。归一化成像条件的量纲接近反射系数，从而使其结果具有了保幅性。这里对这3种成像条件在井间地震中的成像效果进行探讨。 仍旧使用图1的地层模型和参数进行试验。所不同的是图4是震源在300 m深度处的单炮偏移结果，而下面的偏移结果都是多炮叠加后的成像。在模型的左侧均匀放置震源。模型的垂直长度是600 m，每炮的间隔是10 m，共激发61炮。检波器相邻两个的间距为2 m。每激发一炮共接收301道数据并对所得的信号进行偏移成像，将61炮偏移成像的结果进行叠加，就得到了多炮逆时偏移的成像。将图1模型分别对不同的成像条件进行偏移，就得到了不同成像条件下的偏移结果，显示在图5。 图5分别是检波器归一化成像条件、震源归一化成像条件和常规互相关成像条件的结果。从图5(a)中可以看到靠近震源一侧的界面可以被成像出来，同时在震源一测产生了噪声。图5(b)和图5(c)相比，图5(b)中靠近检波器一侧的界面成像被保留和增强，而靠近震源附近的界面成像被减弱。这里对产生这种现象的原因进行解释：归一化的成像条件其实就是对互相关成像条件后的结果进行加权，只是加权的系数由震源或检波器波场进行提供。例如式(2)的分母，由检波器波场计算得到，可以将式(2)的分母看作是一幅在检波器一侧的像素值较大、而在震源附近幅度值较小的图像。用这样的一幅图与互相关成像的结果进行相除，可以使互相关成像中震源一侧的像被增强，而检波器一端的像被减弱，这就是造成在图5(a)中只成像出了靠近震源这侧的界面，而靠近检波器端界面的像几乎消失的原因。由于这种加权作用不仅对界面成像有效，同时也对互相关结果中单个炮点震源附近的波场进行增强，导致在图5(a)震源位置附近产生了噪声，这种噪声可以通过使用更多的炮点叠加来进行降低。同理，震源归一化成像条件会导致震源附近的像被减弱，检波器一侧的像被增强，如图5(b)所示。 通过对图