地震逆信号求取方法与研究.pdfVIP

·848 · 中国地球科学联合学术年会 2014 地震逆信号的求取方法研究 1* 2 2 2 1 黄广谭 梁鸿贤 葛大明 孙朋朋 张军华 1.中国石油大学（华东）地学院 青岛 266580; 2.胜利油田物探研究院 东营 257022 1.引言 在地震数据处理中，很多情况下需要计算反射地震道的逆信号。通常情况下我们采用的方法是对整 个地震道进行最小二乘法求逆，但由于地震道样点数往往较大，利用最小二乘法求解时计算量会很大。 因此，我们在这里给出地震道逆信号的另一种解法，以提高计算效率和计算精度。 2.方法原理 逆地震信号就是将整个地震道的信号进行求逆，通常我们采用的方法是对整个地震道信号进行最小 二乘法求取，当地震数据采样点较小时，我们可以很快地获得该逆信号，但随着地震采样点数的逐渐增 多，求逆所需要的时间将消耗很多的时间。为此，我们提出基于褶积模型和自相关函数的方法。 假设 为地震道，f 为反射系数序列，s 为震源子波， 的褶积模型表达式为， p f ∗s (1) 假设 f 为白噪序列，其自相关信号有以下性质， Corr (f ,f ) δ ⋅ [∑ f (t )f (t )] (2) 0 t 其中， 表示零时间点脉冲信号。假设 表示信号的倒序，则地震道 p 的自相关有， * ( ) * ( ) ( ) * * ( * ) ( * ) Corr p ,p p ∗ p f ∗s ∗ f ∗s f * s * f * s f * f * s * s          (3)  ( ) ( ) [ () ()] ( ) Corr f ,f ∗Corr s,s ∑f t f t δ ∗Corr s,s t 0 即有， −1 * ( )−1 [ ( ) ( )] (4) p p * Corr s,s / ∑tf t f t 舍去常数项，得到地震道的逆 −1 * ( )−1 p p * Corr s,s (5) 3.模型测试及效果分析 为验证本文方法的效果，我们对该方法求取逆信号与将整道进行最小二乘法求逆的计算效果进行了 比较。经比较发现，基于褶积模型和自相关函数法计算效率比最小二乘求逆的计算效率要高，而且随着 采样点数的逐渐增加，本文方法求逆所需时间增长较为稳定，并且增长幅度非常小，而对整道进行最小 二乘法求逆所用的时间增加得很快，当样点值增长到一定数量时，计算时间会相差多个数量级。  4.认识与结论  在地震勘探的处理及反演过程当中，经常会出现对地震逆信号的求取，例如基于 SRME 方法压制 层间多次波的方法中，就需要获得地震道的逆信号。本文提出的方法可以应用于整个地震道求逆的运算 中，式中子波自相关逆的计算可以采用最小二乘法实现。经比较分析，本文提到的方法在计算量上远小 于对整道直接进