地震数据处理第四章：动校正及叠加

第四章动校正及叠加（NMOandStack）本章内容第一节动校正（normalmoveoutcorrection）第二节水平叠加第三节剩余时差及叠加特性第四节水平叠加存在的问题第一节动校正反射波传播时间是地层深度、倾角、介质速度及炮检距的函数。动校正目的：消除炮检距（与观测方式有关）对传播时间的影响。动校正基本假设：均匀介质或等效均匀介质。正常时差：由炮检距引起的时差。反射波时距曲线：（1）共炮点：)sin2(22222xddxvt当倾角=0时，为水平地层；当倾角0时，时间最小点向上倾方向偏移，其横向距离为sin2dxm共炮点反射波时距曲线是以炮点位置的法向深度d为参数的双曲线。（2）共中心点：2/sinxddsm炮点法向深度ds与中心点法向深度dm之间的关系：222,02222,02224)(cos4)(cos2NMOMMvxtvxtxdvtcosvvNMO其中：vdtNMO2)(是动校正速度；它表示反射波的横向视速度，界面倾斜时，它大于地层速度；界面水平时，它等于介质速度。共中心点反射波时距曲线是以中点位置的法向深度d为参数的双曲线。是中心点M的法向双程旅行时，当倾角为0时，则退化为中心点M的垂向双程时。一、动校正的概念2222)0()(vxtxt（4-1）式中x—炮检距；v—界面上覆介质速度；t(x)—x处的地震波旅行时；t(0)—0处的地震波旅行时。反射波时距曲线只有当界面为水平时，共深度点D在地面的投影与中心点M重合。在多次覆盖勘探中，多个炮检距上都接收到来自共深度点的反射波，但反射波在不炮检距的到达时间不同，由于零炮检距反射波与与地下构造有直接对应关系，因此需要将非零炮检距上的反射波旅行时校正到零炮检距旅行时，其时差为：（4-2）)0()0()0()()(22tvxttxtxt这种时差称为正常时差或动态时差（在同一道上，正常时差是反射时间的函数，不同的反射时间，正常时差也不相同，动校正就是把所有的反射波校正到各自的法向时间t0上，所用的叠加速度与地层倾角、深度和上部地层介质有关，相当于等效视速度。）将不同炮检距的反射时间校正到零炮检距反射时的过程称为动校正。“动”的概念体现在同一地震道上不同反射时间的动校正量不同。动校正量是反射时间、炮检距和地层速度的函数，它通常随炮检距递增，随深度和速度递减。CDP道集速度正确同相轴校直速度偏大校正不足速度偏小校正过量单层水平反射界面，动校正速度=上覆介质的层速度；水平层状介质和更复杂的上覆构造，动校正速度与上覆地层速度之间的关系较复杂，反射波时距曲线的轨迹也不再是标准的双曲线。二、水平层状介质的动校正左图为N层水平层状介质模型vi表示第i层的层速度；Δti表示地震波在第i层的垂直双程旅行时间；地震波由震源S点出发，到达反射点D后返回接收点G；地面中心点M与反射点D在同一铅垂线上；炮检距为x。水平层状介质的反射时间t不能表示为炮检距x的显函数，可近似展开为：63422222)0()(xcxcvxtxtrms式中c2、c3—与地层厚度和速度有关函数；vrms——均方根速度。Niiirmstvtv122)0(1（4-3）（4-4）Niitt1)0(，上式表明，水平层状介质情况下，共深度点反射波时距曲线，不再是标准的双曲线。当排列较短、炮检距小于反射点深度时，截断（4-3）式中的高次项，时距曲线可近似为：2222)0()(rmsvxtxt（4-5）短排列近似条件下，水平层状介质的动校正速度近似为地层的均方根速度。三、单一倾斜层的动校正地震波从震源点S出发，传播至反射点D，再返回接收点G。中心点M不再是反射点D沿界面法方向在地表的投影。CMP（共中心点）道集不严格等价于CDP（共深度点）和CRP（共反射点）道集。CMP道集中的反射波不再来自地下同一反射点。单一倾斜层反射波时距曲线：222cos/)0()(vxtxt（4-6）cosvvNMO动校正速度为：（4-7）该方程是双曲线方程地层倾斜，动校正速度大于上覆地层的层速度，产生剩余时差。三维情况，设地层倾角为φ，炮点和检波点连线方向与地层倾向的夹角，即方位角θ，地层在炮检方向的视倾角φ’，有：cossinsincosvvNMO动校正速度为：222222cossin-1)0()(vxtxt三维倾斜地层，反射波时距曲线方程为：（4-8）（4-9）（4-10）四、任意倾斜层状介质动校正地震波从震源点S出发，传播至反射点D，再返回接收点G的旅行时t与炮检距x的关系。中心点M自激自收的反射点D’与炮检距为x的反射点D不再是同一反射点。CMP道集中的反射波来自于一段界面的反射，而不是同一点的反射。Hubral和Krey（1980）导出了下面的反射波时距方程：)()0()(2NMO222xTvxtxtH（4-11）其中112212022coscoscos)0(1ikkkNiiiNMOtvtv（4-12）高阶阶项—)(个界面的双程旅行时i零偏移距地震波在第——t角；个界面的入射角和透射i传播时在第DM零偏移距地震波沿—β,αiiixTH；112212022coscoscos)0(1ikkkNiiiNMOtvtv（4-12）0cosvvNMO为单倾斜层时，式（4-12）变为为水平层状介质时，式（4-12）变为Niiirmstvtv122)0(1当地层倾角不大、排列较短，任意倾斜层状介质反射波时距曲线仍可近似为双曲线2NMO222)0()(vxtxt（4-13）五、数字动校正方法和动校正拉伸设Δt为采样间隔，炮检距为x的地震记录为yx(iΔt)(i=1,…,N)，动校正速度v(iΔt)(i=1,…,N)，动校正后的地震记录为yo(iΔt)见下图，第k个点yo(kΔt)应为)()(0xytky（4-14）其中)()222tkvxtk（（4-15）并未落到的离散样点上，而是在kx和kx+1之间，即无时刻的样值，需要利用相邻几点样值内插或用采样定理来恢复。)(tyx)(tyx*1、数字动校正方法采样定理恢复公式：Nixxtittittiyy1)()(sin)()(（4-16）2、离散动校正对地震记录波形的影响地震记录上的子波由若干离散点组成，由于各离散点的动校正量不同，动校正后的子波将不再保持原来的形态，通常会发生相对畸变。它们分别为：了时间，是校正后子波起始和终了时间；是校正前子波起始和终2121,,222111--（4-17）式中动校正时差。点的和分别是、212121112212--TT由于浅层的动校正时差大于深层的动校正时差，所以TT（4-18）动校正后，子波的波形被拉伸了，称其为动校正拉伸。浅层、大炮检距的拉伸最为严重。为了定量地表示动校正拉伸，引入了拉伸系数的概念，其定义为：TTTTT-动校正前后子波主频的相对变化以及拉伸系数与动校正量和反射时间的关系：子波起始时刻的时距曲线为子波终了时刻的时距曲线为222121vx）（222121)(vxTTT）（（4-20）（4-21）展开(4-21)式后与(4-20)式相减，得到212122）（）（TTTTTT整理后21112-2TTTT）（）（（4-22）（4-23）舍去右端的高阶项，并注意,（4-23）变为11-NMOtTTTtNMO）（1（4-24）当远大于T时，则到子波拉伸与动校正量的关系11NMOtTT（4-25）以子波的主频表示子波的延续时间，有mfmfT1（4-26）于是，得到以子波波形、主频和反射时间相对变化的拉伸系数0ttffTTNMOmm（4-28）由此式可知，反射深度越浅，炮检距越大（越大），动校正拉伸越严重，子波主频向低频转移越严重。动校正拉伸引起波形畸变，破坏了动校正后CMP同相轴的相关性，降低了叠加质量，对纵向分辨率尤为有害，故在叠加之前采用动校正拉伸切除。NMOt动校正误差来自四个方面：（1）地层、构造或岩性等因素破坏前提假设条件；（2）速度误差引起动校正误差；（3）动校正拉伸量随t0的减小而严重，浅层和远炮检距的拉伸比较大；（4）离散采样。第二节水平叠加一、水平叠加的原理设共中心点道集，其中M为样点数，N为道集中的道数，各道已经进行了正常时差校正，要确定一个标准道，使得标准道与各记录道的差别最小，现讨论如何确定这个标准道。利用最小平方原理，计算任意地震道与标准道的误差能量),...,2,1,,...,2,1)((NjMiixj),...,2,1)((Miiy)(ixj)(iy),...,2,1(,)]()([21NjiyixQMijj（4-29）总误差能量NjMijNjjiyixQQ1121)]()([（4-30）),...,2,1(0)(MiiyQ要求（4-31）})]()([{)()(112NjMijkykxiyiyQNjjNjjixNiyiNyix11)(1)(0)()(（4-32）标准道就是N道叠加的平均，这正是水平叠加的理论基础，叠加次数N为有效叠加次数（不含死道、切除道）。二、自适应水平叠加在式（4-32）中，参加叠加各道的加权系数是相等的，而且各道的加权系数不随时间变化，加权系数都为1。实际上，参加叠加的各道的质量是有差别的，当CMP道集中各道的品质差异较大时，等权叠加不会取得理想的叠加效果。1、基本原理地震记录道的质量在空间和时间上都会有差异，可以根据它们在空间和时间上质量的差异来控制它们参与叠加的成分，使用不同的权系数，用最小平方法来确定加权系数。NjjixNiy1)(1)(（1）标准道的形成能较好地反映叠加剖面特征的地震道，可以是相邻几个CMP道集的叠加，也可以是在地震剖面上进行了信号增强处理后的自适应叠加道。（2）求加权系数第j个地震道xj(t)乘上加权系数wj(t)后，应该与标准道y(t)在最小平方意义下最接近。某一段地震记录，其中心时间t0,时窗长度为T，该时窗内加权后的地震记录xj(t)wj(t)与标准道y(t)的误差能量为：2/2/200)]()()([TtTttjjjtytwtxQ（4-33）])2/,2/[(0)(000TtTtttwQjj要求：由此解得2/02/02/02/00)()()()()(TtTttjjTtTttjjtxtxtxtytw（4-35）计算步骤：（1）计算标准道（2）计算权系数按（4-35）计算，并异常值编辑、平滑滤波。（3）地震道加权（4）加权道叠加)()()(twtxtrjjjNjjtrNts1)(1)(N为叠加次数第三节剩余时差及叠加特性1、剩余时差各种干扰都有自已的正常时差；一次反射波的正常时差随水平界面和倾斜界面不同而不同；用一次波正常时差校正公式校正后，除一次波外，其他类型的波仍存在一定量的时差，将动校正后残留的时差叫剩余时差。2、全程多次波剩余时差全程多次波时距曲线方程为22202ddvxtt（4-38）假设一次反射和多次反射的t0相等，按一次反射作动校正后，多次波的剩余时差为22202)11(2qxvvtxtdd（4-39）其中q为剩余时差系数)11(21220vvtqd（4-40）由（4-39）式看出，剩余时差与炮检距的平方成正比，其轨迹近似抛物线形。各叠加道的剩余时差不同，叠后多次反射被削弱。由此知，水平叠加技术是利用各波剩余时差的不同来压制干扰波的一种方法，它能较好的压制多次反射波。3、叠加特性设地面共中心点的一次反射