地震数据处理-反褶积课件

第三章反褶积3.1、反褶积及褶积模型3.2、反滤波3.3、最佳维纳滤波与最小平方反褶积3.4、脉冲反褶积3.5、预测反褶积3.6、子波整形反褶积3.7、同态反褶积3.8、地表一致性反褶积第三章反褶积3.1、反褶积及褶积模型1

反褶积又称反滤波。为了消除大地滤波及接收系统滤波对地震数据的影响而作出的滤波处理。反滤波本质上是一种频率滤波。从数学上看，它是一种褶积运算，故称反褶积。3.1、反褶积及褶积模型一、反褶积的概念反褶积又称反滤波。为了消除大地滤波及接收系统滤波2反褶积处理：是常用处理方法之一。可以用于叠前和叠后，也可以多次使用。作用：压缩地震子波，提高分辨率。可以压制多次波和短周期鸣震等干扰，提高地震资料信噪比。反褶积处理：3震源爆炸使地下介质形成三个区域：

地层对震源脉冲的改造作用，相当于对它进行了一次低通滤波，此滤波器常称为大地滤波器。

震源爆炸产生尖脉冲传播到弹性区起始边界时，已经变成了有一定延续时间的稳定波形——地震子波。震源爆炸使地下介质形成三个区域：地层对震源脉冲的改造4

假设震源脉冲在地下介质中传播未受大地改造，脉冲信号入射到分界面、反射信号返回地面，被检波器接收、传输到仪器被记录下来。如果接收系统未对震源脉冲进行改造，则地震记录为反射系数序列：反褶积就是要获得未经系统作用的地震波形。假设震源脉冲在地下介质中传播未受大地改造，脉冲信号5

地震子波同震源子波，其概念是有区别的，它与许多因素有关。根据地震波传播过程中影响因素的不同，地震子波可描述为：式中而地震子波同震源子波，其概念是有区6干扰波是由非激发干扰（次生）、背景噪声

及规则（或称相干）（由激发产生）干扰

叠加而成:规则干扰分两类：

一类与地质构造有关，包括多次波、转换波、绕射波、伴随波、折射波、瑞利波、勒夫波和斯通利波等,这类波在特定的条件下可转化为有效波;另一类与地质构造无关,如水中震鸣、气泡效应、地表及海面散射等(也包括地下震鸣、薄层微曲多次波)。实际处理时,要根据不同的勘探情况,分别对待。

干扰波是由非激发干扰（次生）、背景噪声规则干扰7反褶积的关键

是如何设计一个反滤波器去抵消另一个滤波器的作用。设计反滤波器的方法：

反褶积的关键8

由已知地震子波计算反褶积算子，称确定性反褶积，主要用于去除记录系统的响应、海上震源子波反褶积等方面；通过统计方法求取最佳反褶积算子，如脉冲反褶积、预测反褶积等。由已知地震子波计算反褶积算子，称确定性反褶积，主9二、褶积模型理想模型：加噪模型：其中：（3-1）（3-2）二、褶积模型理想模型：加噪模型：其中：（3-1）（3-2）10反褶积的假设条件：（1）地下地层是水平层状介质；（2）地震波是垂直入射反射的平面波；（3）地震子波在传播过程中保持波形不变；（4）地震记录中无噪声；（5）地震子波已知；（6）反射系数序列为白噪序列；（7）地震子波是最小相位的。若假设条件与实际不吻合，势必会造成褶积模型与实际地震记录存在一定差异。反褶积的假设条件：11改进模型：①海上“特征反褶积”模型其中：改进模型：①海上“特征反褶积”模型其中：12②沙漠地区可控震源地震记录模型其中：②沙漠地区可控震源地震记录模型其中：13③仿真褶积模型其中：③仿真褶积模型其中：14

为了把地震子波压缩成尖脉冲（必需去掉大地滤波器的作用），使地震记录变为反射系数序列，出现了各种反褶积方法，而实际处理结果往往不如人愿。其原因有三：①地震记录已知，地震子波未知，求反射系数序列，必须有若干假设条件限定解的唯一性，否者是多解的；假设条件与实际情况越接近，反褶积效果越好。为了把地震子波压缩成尖脉冲（必需去掉大地滤波器的作15②反褶积方法依赖地震记录的褶积模型，模型中地震子波是大地滤波器的脉冲响应，而大地滤波的作用复杂，模型不太可靠。只有先彻底解决正演问题，才能使反褶积得到发展。③反褶积方法可能会提高噪声水平，有必要同时发展提高分辨率及信噪比的方法。反褶积方法很多，有些（如最大熵、卡尔曼、时变Q等）未能在常规处理中获得一席之位。②反褶积方法依赖地震记录的褶积模型，模型中地震子波是大地滤波16反射系数剖面地震剖面反射系数剖面地震剖面17地震数据处理-反褶积课件18地震数据处理-反褶积课件19第三章反褶积3.1、反褶积及褶积模型3.2、反滤波3.3、最佳维纳滤波与最小平方反褶积3.4、脉冲反褶积3.5、预测反褶积3.6、子波整形反褶积3.7、同态反褶积3.8、地表一致性反褶积第三章反褶积3.1、反褶积及褶积模型203.2反滤波一、反滤波的概念1、概念3.2反滤波21假定地震记录不含干扰，即2、反子波对应的频域形式则可得到假定地震记录不含干扰，即2、反子波对应的频域形式则可得到22写成时域形式为：反子波与子波褶积为：

由子波和反射系数求地震记录，是一褶积过程（正演）；已知反子波和地震记录求反射系数，称为反褶积或反滤波。写成时域形式为：反子波与子波褶积为：由子波和反射系数23二、地震子波的求取

确定性反褶积，需已知子波。故先讨论子波求取方法，有5种方法：直接观测法（适用于海上）；自相关法；多项式求根法；利用测井资料求子波；对数分解法。确定性反褶积处理步骤：先提子波，再求反子波，然后进行反褶积。二、地震子波的求取确定性反褶积，需已知子波。故先讨242、自相关法选择一段质量较高的地震记录，时窗长度为T：其Z变换为假设反射系数是白噪声序列，其z变换为则的自相关的z变换：2、自相关法选择一段质量较高的地震记录，时窗长度为T：其Z变25地震记录的z变换为地震记录自相关的z变换为将代入，有：由于都是实数序列，所以有：地震记录的z变换为地震记录自相关的z变换为将26因此有：也有：未知，现在来确定它（3-19）（3-18）假如地震子波是最小相位的物理可实现序列，则其z变换为：因此有：也有：未知，现在来确定它（3-19）（3-18）假如27由物理可实现性知：当时，,对下式（3-19）两端取对数，有：令因而得到（3-20）由物理可实现性知：当时，28根据复变函数理论其中C为常量，和因而得到：因为是实数，由（3-20）式知，于是可得C根据复变函数理论其中C为常量，和因而得到：因为29（3-20）对上求希尔伯特变换，即可求出相位谱：令（3-19）求出相位谱后，对下式求付立叶反变换，得最小相位子波（3-23）复杂内容简单化：（3-20）对上求希尔伯特变换，即可求出相位谱：令（3-1930假如地震子波是零相位的，下式（3-19）中的相位谱为零，即，因此有零相位地震子波：（3-32）假如地震子波是零相位的，下式（3-19）中的相位谱为零，即313、多项式法选择一段质量较高的地震记录

设反射系数为白噪声序列，则记录自相关与子波自相关等价，即褶积模型：令，则（3-34）3、多项式法选择一段质量较高的地震记录设反射系数为白噪声序32将上式两端乘以，则有：由于显然，应有2M个根。鉴于系数均为实数，所以2M个根是M对互为倒数的，即若则另一根为：将上式两端乘以，则有：由于显然，应33

根据这M对根在单位圆内、外的位置，可以组成2M个不同相位的地震子波，其中必有一个是最小相位，一是最大相位的。根据“最小相位序列z域零点在单位圆内”这一特点，选出模小于1的根，便可组成最小相位子波，其z变换为：令z=0,得由此得最小相位子波根据这M对根在单位圆内、外的位置，可以组成2M个不34例如，已知记录的自相关其z变换为：两端乘以z2得：求的根，有选出模小于1的根，求得最小相位子波为：例如，已知记录的自相关其z变换为：两端乘以z2得：求354、用测井资料求取子波

这种方法要求有较好的声波测井和密度测井资料，并在井旁有质量较高的地震记录．

首先将声波时差转化为声波速度设声波单位为，则然后，进行深时转换，是双层旅行时4、用测井资料求取子波首先将声波时差36

然后，计算反射系数对井旁地震记录和反射系数进行傅氏变换后，可得到子波的频谱对傅氏反变换就得到地震子波,即：然后，计算反射系数375、对数分解法时域模型：此法不需假设反射系数是白噪声，不需假设子波是最小相位。频域模型：对频域模型两端取对数，则将子波与反射系数分离开来称为对数谱。5、对数分解法时域模型：此法不需假设反射系数是白噪声，不需假38用付立叶反变换对数谱的时间信号：由于分布在时间轴原点附近；分布离原点较远区域。若二者分离较好，则可用低通滤波将分离出来，便可求出子波。即先付氏正变换，再取指数，然后进行付氏反变换用付立叶反变换对数谱的时间信号：由于分布在时39

由于很难确定对数谱在时间轴上的分布区域，故用N道地震记录来求取。求平均值：由于，所以因此可获得较准确的子波。由于很难确定对数谱在时间轴上的分布40三、反滤波的实现求出地震子波后，可用付氏变换或z变换或最小平方法求反子波，然后对地震记录做褶积运算得。如三、反滤波的实现求出地震子波后，可用付氏变换或z41第三章反褶积3.1、反褶积及褶积模型3.2、反滤波3.3、最佳维纳滤波与最小平方反褶积3.4、脉冲反褶积3.5、预测反褶积3.6、子波整形反褶积3.7、同态反褶积3.8、地表一致性反褶积第三章反褶积3.1、反褶积及褶积模型423.3最佳维纳滤波及最小平方反褶积一、最佳维纳滤波

维纳滤波即最小平方滤波,它是使实际输出与期望输出，在误差能量最小条件下，求滤波因子的方法，所代表滤波器是最佳滤波器。求解关系

输入信号滤波因子实际输出

期望输出误差能量3.3最佳维纳滤波及最小平方反褶积一、最佳维纳滤波43最小平方最小最小平方最小44滤波方程写成矩阵形式当反射系数为白噪声时，记录和子波自相关等价，也就不用求子波了。滤波方程写成矩阵形式当反射系数为白噪声时，记录和子波自相关等45地震数据处理-反褶积课件46表3-1常用的期望输出根据滤波目的设定期望输出表3-1常用的期望输出根据滤波目的设定期望输出47二、最小平方反褶积

将最佳维纳滤波原理应用于反褶积问题，就是最小平方反褶积方法。地震记录为：地震记录含随机噪声后,褶积模型变为：设计滤波因子为：二、最小平方反褶积将最佳维纳滤波原理应用于反褶积问题48

输出为：误差的平方和为：输出为：误差的平方和为：49式中：式中：50表明：地震记录的自相关等于子波自相关与噪音自相关的和表明：地震记录的自相关等于子波自相关与噪音自相关的和51地震数据处理-反褶积课件52表明：地震记录与期望输出的互相关等于期望脉冲因子与子波的自相关表明：地震记录与期望输出的互相关等于期望脉冲因子与子波的自相53写成矩阵形式（3-76）（3-75）写成矩阵形式（3-76）（3-75）54

地震记录含噪声后,求得的反子波与不含噪时会有差别(托布里兹矩阵主对角线元素上加了噪声的自相关值e),但这一点并不影响反褶积的效果,反而可增加方程组求解的稳定性。有时为了实际需要,要人为加进一些噪声,这就是预白化问题。地震记录含噪声后,求得的反子波与不含噪时会有差别(托55地震数据处理-反褶积课件56例1:（1）求实际输出序列（输入：最小相位子波，期望输出：零延迟脉冲）按原理求例1:（1）求实际输出序列（输入：最小相位子波，期望输出：零57（2）求滤波因子（2）求滤波因子58地震数据处理-反褶积课件59例2:（输入：最大相位子波，期望输出：任意延迟脉冲）按方程求例2:（输入：最大相位子波，期望输出：任意延迟脉冲）按方程求60地震数据处理-反褶积课件61地震数据处理-反褶积课件62三、最小平方反褶积中的预白化处理1.问题的提出对于反滤波方程可以得到反子波的频谱为此时,如果是带限信号,那么在频带外计算出的将会很大,甚至趋于无穷.从信号分析的角度讲,这样的滤波器是不稳定的.

实际情况是肯定是带限信号,这时如何解决带限问题呢?（3-43）（不讲）三、最小平方反褶积中的预白化处理1.问题的提出（3-43）（63

对（3-43）式描述的是频谱形式,对其改造为

（3-43’）为了解决带限问题,在地震信号的功率谱P(ω)中,从低频到高频统一加一白噪声。2.预白化处理（3-44）对（3-43）式描述的是频谱形式,对其改造为2.64当反射系数为白噪声时，子波的自相关=地震记录的自相关。当反射系数为白噪声时，子波的自相关=地震记录的自相关。65地震数据处理-反褶积课件66从以上方程可以看出：

从以上方程可以看出：67第三章反褶积3.1、反褶积及褶积模型3.2、反滤波3.3、最佳维纳滤波与最小平方反褶积3.4、脉冲反褶积3.5、预测反褶积3.6、子波整形反褶积3.7、同态反褶积3.8、地表一致性反褶积第三章反褶积3.1、反褶积及褶积模型683.4脉冲反褶积一、脉冲反褶积原理

当期望输出为时，除外，其余全为零。方程（矩阵形式）：3.4脉冲反褶积一、脉冲反褶积原理当期望输出为69（3-80）（3-79）脉冲反褶积（对应3-80式）要求输入子波为最小相位的，实际中往往是非最小相位的，因此要作某种优化处理。（3-80）（3-79）脉冲反褶积（对应3-80式）要求输入70

最小平方反滤波，对子波相位已无要求，可以任意。

反滤波因子在理论上应为无穷项，实际只能取有限项，并且其主要部分反滤波因子的形状和位置要根据地震子波的形状来决定。最小平方反滤波，对子波相位已无要求，可以任意。71详情见书p78最小延迟，m0=0，m最大延迟，m0=-n-m混合延迟，m0=-m1详情见书p78最小延迟，m0=0，m最大延72二、参数选择

（效果好的标准是什么？子波是否被压缩）（1）反滤波因子长度m，通过试验来选择；（2）相关时窗长度m+n，常为2m~记录长度；（3）稳定常数，一般取的百分数。干扰小，取0.005~0.01；干扰大，取0.02~0.05。二、参数选择（1）反滤波因子长度m，通过试验来选择；（2）相73下图是延迟脉冲反褶积效果分析图

（误差能量）延迟时间输入为非最小相位子波，期望输出下图是延迟脉冲反褶积效果分析图（误差能量）74第三章反褶积3.1、反褶积及褶积模型3.2、反滤波3.3、最佳维纳滤波与最小平方反褶积3.4、脉冲反褶积3.5、预测反褶积3.6、子波整形反褶积3.7、同态反褶积3.8、地表一致性反褶积第三章反褶积3.1、反褶积及褶积模型753.5预测反褶积一、预测滤波原理

预测问题是已知某个物理量的过去值和现在值,通过对已知信息加工处理来获得未来某个时刻的预测值。其数学描述如下:

设为现在值，为过去值,如果定义预测步长为，用现在值和过去值来预测将来时刻的预测值，即（3-83）3.5预测反褶积一、预测滤波原理预测问题是已知某个76使预测值与实际未来值的误差—预测误差最小，按最小平方原理，得如下方程：（3-86）矩阵形式为：（3-87）使预测值与实际未来值的误差—预测误差最小，按最小平方原理，得77二、预测反褶积原理预测反褶积是要从含有多次波的地震记录中，根据多次波具有周期性特点，预测出多次波来，用前者减去后者得含一次波的地震记录（预测误差）。即：地震记录=一次波（不可预测）+多次波（可预测）地震记录-多次波（可预测）

=一次波（不可预测）1、预测反褶积原理二、预测反褶积原理预测反褶积是要从含有多次波的地震记录中，根78

设地震子波满足最小相位条件，反射系数为白噪声,褶积模型为则时刻的输出值为：含未来时的信息含过去、现在时的信息设地震子波满足最小相位条件，反射系数为白噪声,褶积模79比较如下两式：可以看出:Ⅱ是预测值(即鸣震干扰),Ⅰ是预测误差(一次反射波).

设即相当于取了子波的前部分,则①式变为Ⅰ的意义：①②比较如下两式：可以看出:设即80②式表明：将一个子波的前部分与反射系数的褶积就得到了预测误差（一次反射波）。

这种方法也压缩了子波的长度，从而提高地震资料的分辨率。两步法实现过程：设输入序列为5点序列：，预测步长，期望输出：②式表明：两步法实现过程：设输入序列为5点序列：81①、求预测因子②、求预测值③、求预测误差即先做预测反滤波，再求预测误差，故称两步法。能否直接求预测误差？必须找出预测反滤波因子和预测误差反滤波因子的关系①、求预测因子②、求预测值③、求预测误差即先做预测反滤波，82预测反褶积数学描述：其中为预测误差反滤波因子，它为：③预测反褶积使子波长度变为，故也称其为子波切除反褶积。当时，预测反褶积变为脉冲反褶积。预测反褶积数学描述：其中为预测误差反滤波因子832、预测因子的求取即：输入信号预测因子实际输出期望输出2、预测因子的求取即：输入信号843、计算举例已知地震记录存在虚反射,其时间序列为预测反褶积来消除虚反射。解：设预测因子长，预测步长。3、计算举例已知地震记录存在虚反射,其时间序列为预测反褶积来85(2)解方程得预测因子：(3)求虚反射,即预测值:(2)解方程得预测因子：(3)求虚反射,即预测值:86(4)预测误差,消除虚反射后的那一段地震数据(5)压制虚反射后的地震记录原地震记录周期为5因周期为5，故预测步长取5(4)预测误差,消除虚反射后的那一段地震数据(5)压制虚反射87或MATLAB程序：或MATLAB程序：88地震数据处理-反褶积课件89地震数据处理-反褶积课件90三、预测反褶积压缩反射脉冲（自学）结论：三、预测反褶积压缩反射脉冲（自学）91四、参数选择（1）预测步长，关键参数，效果与此有关；（2）预测因子长度，要适宜，小，波形尾部波动；大，增加运算量，相邻反射及噪声影响；（3）预白化量，增加求解的稳定性，需根据噪声水平来定。后有三图，以助说明。四、参数选择（1）预测步长，关键参数，效果与此有关；92预测因子长度128ms，预测步长2ms，效果最好，其余有延续时间（非脉冲信号）。（a）为脉冲序列，（b）为合成记录预测因子长度128ms，预测步长2ms，效果最好，其余有延续93预测步长2ms,预测因子长度128ms，效果好。小长度因子，尾部波动明显预测步长2ms,预测因子长度128ms，效果好。小长度因子，94预白化量0.1，效果较好；预白化量增大，尾部波动增大，效果变差。预白化量0.1，效果较好；预白化量增大，尾部波动增大，效果变95第三章反褶积3.1、反褶积及褶积模型3.2、反滤波3.3、最佳维纳滤波与最小平方反褶积3.4、脉冲反褶积3.5、预测反褶积3.6、子波整形反褶积3.7、同态反褶积3.8、地表一致性反褶积第三章反褶积3.1、反褶积及褶积模型963.6子波整形反褶积１．子波的相位与分辨率子波的相位通常有三种，即最小相位，混合相位和最大相位，这些子波是单边的物理可实现信号．海上勘探和陆上勘探爆炸震源产生的地震子波都接近最小相位，所以，经常假设和讨论最小相位．还有一种相位的子波是零相位子波，虽说是一种非因果信号，是物理不可实现的，但在数字滤波、反褶积和反演中经常用到。3.6子波整形反褶积97为什么要用到零相位子波，基于以下两个考虑：①陆上可控震源子波因是通过自相关处理得到的，所以是零相位的；②零相位子波，比其他相位子波的分辨率高．为什么要用到零相位子波，基于以下两个考虑：98

２．子波与反子波的时域分布的特点对于褶积公式

假设子波为有限长度，

从理论上来说，我们可以得到无限长的子波利用Ｚ变换，我们可以推导出有限长子波与所求反子波的关系２．子波与反子波的时域分布的特点99３.相位对反褶积精度的影响

在前面的讨论中，为了讨论问题方便，都假设了子波为最小相位．实际上子波有不同的相位．在反褶积中，其他条件相同，相位不同，反褶积的精度是不同的．对于最小平方反褶积，滤波方程为３.相位对反褶积精度的影响在前面的讨论中，为了讨论问题100

我们分两种相位的期望输出来讨论反褶积问题．１）期望输出为（0，０，1）这时的方程为

解得反子波为，反褶积输出为与期望输出（0，0，1）的误差能量为．我们分两种相位的期望输出来讨论反褶积问题．101

从中我们也可以看到，此时的误差比最小相位是１／２１大多了．２）期望输出为（0,1,0）这是的滤波方程为得到的反子波为，反褶积输出为这与期望输出(0,1,0)的误差能量为．式中我们可以看到此时的反褶积比上面的精度要高

从中我们也可以看到，此时的误差比最小相位是１／２１大多了．102

由此我们可以得到一个重要的结论：①子波振幅相同时，最小相位子波对期望输出为零延迟的反褶积，误差最小；②在子波为混合相位和最大相位时，期望输出的相位应与子波的相位相匹配，有一个最佳延迟，只有这样才能得到合适的反褶积结果．由此我们可以得到一个重要的结论：103最小相位子波，期望输出波形为零延迟尖脉冲；最大相位子波，期望输出波形为最大延迟尖脉冲；混合相位子波，期望输出波形为零延迟尖脉冲与最大延迟尖脉冲之间的非零延迟尖脉冲。最小相位子波，期望输出波形为零延迟尖脉冲；104地震数据处理-反褶积课件105地震数据处理-反褶积课件106第三章反褶积3.1、反褶积及褶积模型3.2、反滤波3.3、最佳维纳滤波与最小平方反褶积3.4、脉冲反褶积3.5、预测反褶积3.6、子波整形反褶积3.7、同态反褶积3.8、地表一致性反褶积第三章反褶积3.1、反褶积及褶积模型107第三章反褶积3.1、反褶积及褶积模型3.2、反滤波3.3、最佳维纳滤波与最小平方反褶积3.4、脉冲反褶积3.5、预测反褶积3.6、子波整形反褶积3.7、同态反褶积3.8、地表一致性反褶积第三章反褶积3.1、反褶积及褶积模型1083.8地表一致性反褶积目的：消除由于近地表变化对地震子波波形的影响。

地表一致性谱分解（Taner1981）认为任意道的频谱是对应的炮点、接收点、共中心点和炮检距的频谱响应的乘积，对其取对数，可分解炮点和接收点的振幅响应和相位响应。振幅响应被用于振幅一致性补偿，相位响应被用于时差校正。

地表一致性谱分解的基础是褶积模型和预测原理，并扩展为多道平均，把等效子波分解成多个分量。3.8地表一致性反褶积目的：消除由于近地表变化对地震子波109地震数据处理-反褶积课件110地表一致性模型：其中：地表一致性模型：其中：111求频谱：其振幅谱和相位谱分别为：对振幅谱取对数:求112

设是实际地震道的振幅谱，根据地表一致性的假设，一个地震道由相应分量组成，而每一分量对相关的地震记录的贡献是一样的，则实际振幅谱与模型估计振幅谱之差的最小平方和为：使误差能量达到最小，得到关于4个分量的方程组。用迭代法求解取每个分量，在最小相位条件下，估算它们的反算子，进行反褶积。设是实际地震道的振113

与单道反褶积的区别：炮点（或接收点）分量里仅含与观测点有关因素，对相关的所有道是相同的。与单道反褶积的区别：炮点（或接收点）分量里仅含与观测114一．方法原理１．褶积模型对于地震道有式中w(t)为一中和子波，其他符号意义同前．考虑到地表四个因素：

共炮点共接收点共中心点共炮检距点（3-118）一．方法原理共炮点共接收点共中心点共炮检距点（3-118）115式中式中116上式在频域式分别写出它们的振幅谱和相位谱的形式上式在频域式分别写出它们的振幅谱和相位谱的形式1172.对数谱分析方法以上子波包括4个振幅谱和4个相位谱，假设w(t)为最小相位子波，这样只须讨论子波的振幅谱。为了计算方便，对（3-120）式两边取对数，得到于是，褶积关系变成相加关系，下标i代表第i炮，下标j代表第j道。简记对数振幅谱为，并误差函数为

（3-123）2.对数谱分析方法118

这样在每个i、j上，可由（3-122）得到一个，它与实际地震道有误差。令E最小，用Gauss-Seidel方法，可以求出4个对数振幅谱分量，再用反对数变换得到他们的振幅谱。

３．反褶积因子的计算有了子波的振幅谱后，再假定子波为最小相位，可求出４个反子波，对数据道相继采用这些反子波，就可以完成地表一次性反在褶积．这样在每个i、j上，可由（3-122）119二、地表一致性反褶积的应用

１．反褶积之前的有关处理

在反褶积之前，必须先进行去噪声处理，去除面波以及规则的干扰；同时要做好各种振幅恢复补偿处理．

２．时窗的选取

为了得到各个频率分量的振幅谱和自相关函数，要控制剖面上地质构造的形态，划分时窗，每一道时窗按半个时窗互相重叠．得到的各个分时窗的自相关函数后，再做统计平均，作为托布里兹矩阵的自相关函数．二、地表一致性反褶积的应用120

３．反褶积因子的选取

用以上４个分量实现的反褶积，计算量十分巨大，是普通反褶积的８０倍．在同一ＣＭＰ道集中，共中心位置分量来自相同地下界面，计算时可以忽略．因此，通常情况下，只采用３个分量对道进行反褶积．但计算量还是很大的，实际使用时，不是对每一分量求反褶积因子，而是把３个分量的振幅谱求和，求一个综合的反褶积因子．另外，还可以在时间域，用迭代的方法求解．３．反褶积因子的选取121

如果忽略检波点的因素，只考虑炮点和接收点的情况，地表一致性反褶积实际变成了两步法统计子波反褶积．再退一步，如果只做炮点的反褶积，就变成一步反褶积，此时已不再称为一致性反褶积了．

４．效果分析

图3-19是爆炸源的叠加地震剖面分别做普通反褶积和地表一致性反褶积的结果．比较可看出，地表一致性反褶积信噪比要比普通反褶积高，并且振幅变化均匀，特别是在０.２～０.４之间的浅层．如果忽略检波点的因素，只考虑炮点和接收点的情况，地表一致122

地表一致性反褶积的主要目标是校正子波的振幅谱，所以，这种方法适合于地表条件变化大的地区．但是真正方法并不是着重展宽频谱，所以，分辨率并不能提高很多．地表一致性反褶积的主要目标是校正子波的振幅谱，所以，这种1233.9反Q滤波及谱白化

由于地层的吸收作用,地震波经地层传播后,能量被衰减损耗,频率变低,特别是深层,分辨率大大下降.因此,为恢复地震波原能量,必须做吸收补偿,即Q补偿或滤波.3.9反Q滤波及谱白化由于地层的吸收作用,地震波1241、Q因子的物理含义波在介质中传播，产生的弹性能逐渐被介质吸收，最后转换为热能，其过程称为吸收。一个波长内，原地震波能量与传输所损耗的能量之比为Q。1、Q因子的物理含义波在介质中传播，产生的弹性能逐渐被介质吸1252、Q值的估算(1)振幅包络法；(2)李氏经验公式：3、反Q滤波的方法实现设反Q滤波器为最小相位，振幅谱为求出滤波因子，便可滤波。2、Q值的估算3、反Q滤波的方法实现设反Q滤波器为最小相位，126地震数据处理-反褶积课件127•SetDeconvolutionTypetoMinimumPhasePredictive•UseanOperatorLength(ms)of250•SetthePredictionDistance(ms)to96•LetWhiteNoise(percent)be.1,thedefault•UseonetimegatewithaStartTime(ms)of700andanInterval(ms)of400.PredictiveDeconvolutionProcessingParameterOptionsandTimeGatedesignforDeconvolution•SetDeconvolutionTypetoMi128Line127With96msPredictiveDeconvolutionLine127WithoutPost-StackProcessingAppliedLine127With96msPredictive129Line127With12msPredictiveDeconvolutionLine127With24msPredictiveDeconvolutionLine127With12msPredictive130PowerSpectrumfromline127withnodecon.Powerspectrumofline127afterminimumphasepredictivedeconvolutionwitha36msgapPowerSpectrumfromline127w131PowerSpectrumfromline127withnodecon.Powerspectrumofline127afterminimumphasepredictivedeconvolutionwitha24msgapPowerSpectrumfromline127w132PowerSpectrumfromline127withnodecon.Powerspectrumofline127afterpre