反褶积处理方法

辨率的主要手段。在Robinson褶积模型中，均假定地震子波是最小相 位，并且是时不变的，反射系数是白噪的。实际地震记录一般不满足这 些假设，因此直接进行常规的反褶积处理就很难得到预期的效果。另一 方面，在反褶积处理中，期望输出的选择也是影响反褶积效果的重要因 素。在流行的反褶积软件中，常用的期望输出有脉冲、带通子波和 Ricker子波。这些子波在反褶积处理中虽然见到了一定的效果，但也存 在难以克服的缺陷。针对这些问题，研究出了地表一致性俞氏子波反褶 积处理技术，并在CONVEX-SPP大型并行机上开发了地表一致性俞氏子波 反褶积地震处理软件，取得了满意的处理效果。
本文讨论两种实现地表一致性反褶积方法，它基本上代表两种不同 类型的实现方式，并且在生产实践中均已经得到应用。一种是在频率域 内，谱分解法；另一种是建立反褶积输出判别准则的时域迭代法。 1、谱分解法（复谱法）
利用谱分解法来实现地表一致性反褶积，主要有三大步骤，即 谱分析、谱分解和反滤波因子的应用。图2是地表一致性反褶积处理流 程，图2中的B框进一步划分为图3。
域反褶积及频率域预测反褶积。 （7）为改进可控震源反褶积的效果，提出了不相关可控震源反
褶积方法。 （8）直接利用子波整形来解决问题的有子波均衡法。 （9）一种全新的发展是用人工神经网络作最小方差反褶积。 二、反褶积处理方法原理 一个记录地震道的最简单的模型是：
S(t)=W(t)*r(t)+n(t) （1）
如果我们认为，地震数据道是由地震子波和反射系数褶积而成，即
（3）
那么，为了消除对分辨率的影响，我们可以设计一个反滤波因子，
使其
（4）
把作用在上，就可得到反射系数序列：
（5）
只能得到的一个近似解而不能得到精确的解。用表示近似解所产生
的误差，用表示总的误差能量。那么有
（6）
和
（7）
据最小平方的设计思想，是要使误差能量取最小值。因此我们求对
（1）结果与目的背离 反褶积的目的是由地震数据恢复反射系数。作用主要是压缩地震
反射脉冲的长度，提高反射地震记录的分辨能力，并进一步估计地下反 射界面的反射系数。一般的确定性反褶积子波估计并不准确，达不到准 确估计反射系数的效果。而估计性反褶积就只能是拓宽频带，提高一些
分辨率而已。这就出现了一个有趣的现象反褶积背离了它的初衷，将提 取反射系数的任务退化为只是提高分辨率了。
（15） 求取，使达到最小值。要满足：
所以
（16）
上式就是最小平方预测滤波方程组，解方程组可得预测滤波因子，
与褶积可得。实际使用时，我们同样使用记录的自相关，来代替子波的
自相关，从而引入相应的假设条件。到此为止，我们完成了预测反褶积
的第一步：预测滤波。
第二步实现预测反褶积。求出以后，有两种途径来实现预测反褶
最普通的反褶积方法是反滤波。它具有这一性质，即反褶积后的反 射系数其频带是有限的。这种带限的后果是不能清晰地分辨反射界面， 这就会给详细解释带来严重问题。为了克服带限的不良响应，提出了各 种各样的反滤波方法。 （一）脉冲/预测反褶积
该模块用维纳-莱文森最小二乘法来设计并应用脉冲/预测反褶积算 子，作为一个选件，可以纯相位或纯振幅模式来应用它。当预测步长为 1个样点时，预测反褶积就是脉冲反褶积，因此在某些应用系统中，这 两个数学模型被合并在一个模块之中。 1、脉冲反褶积
相关等于子波的自相关。由（3）式
把（9）和（10）式代入方程组（8）就得到了脉冲反褶积的法方程
系。考虑到自相关函数的对称性，即，则有：
（11）
或者写成：
（12）
在推导（12）式时，引入了两个假设条件，即子波是最小相位的，
反射系数序列是白噪的。当分析脉冲反褶积的效果时，应时刻不要忘记
这两个应用的基本前提假设。由（12）式求出反褶积因子与地震数据进
正文
一、反褶积 （一）研究目的和意义
1、研究目的 （1）弄清各种反褶积处理方法的原理。 （2）弄清反褶积处理模块的参数意义。 （3）掌握地震资料数字处理的基本流程及处理方法。 （4）完善反褶积方法，提高地震资料处理的分辨率，保持信噪
比，振幅均匀化。 2、研究意义
反褶积是地震资料数字处理流程中最关键的一环，也是提高地震勘 探分辨率最有效的方法。一个处理流程包括许多处理步骤。而每一个处 理步骤又要涉及到好几个处理模块。一个处理流程通常由预处理、叠前 处理和叠后处理三部分组成。其中反褶积是最重要的一个部分，如图1 所示。
行褶积，就完成了脉冲反褶积的处理。
2、预测反褶积
预测反褶积在某种意义上可以说是一种更广义的最小平方反褶积，
它能包括脉冲反褶积，并能用于研究一般的反褶积问题。预测反褶积的
起由，是为了消除长短周期的多次波及混响对一次波的影响。做法上， 是根据地震记录一次反射和干扰的信息预测出纯干扰部分，再由包括一 次波和干扰的地震记录中减去纯干扰部分，得到消除干扰后的一次反射 信号，以消除一次反射后面的海上鸣震等多次波干扰。
预测反褶积是分作两步来实现的。第一步是预测时间以后的波形信 号，第二步是从信号中减去预测出的信号。不难看出，第一步是要求一 个预测滤波因子，建立最小平方预测滤波数学模型。设输入信号为，要 求的滤波因子为，即,期望输出为，称为预测时间或预测步长，实际输 出是：
（13） 输出误差为：
（14） 总的误差能量E为：
脉冲反褶积一般属于叠前处理手段，是生产实践中最常用的也是最 基本的一种反褶积方法，它的期望输出是一个尖脉冲，即
（2） 在常规处理流程中，原始记录输入以后，做一个高通滤波，压制很 强的面波干扰，紧接着就是脉冲反褶积。因此，它还是一个处理初期使 用的基本模块。一般情况下，程序设计给用户安排反褶积因子长度、相 关时窗长度、白噪系数这几个处理参数，这些参数直接控制着反褶积的 处理效果，通常要通过实验来进行选择，而且要根据剖面的具体情况进 行时变和空变。
式中，W(t)是地震子波，r(t)是反射系数，n(t)是附加噪声，而 S(t)是地震信号。它是双程旅行时t的函数。实践中，此信号是经过采 样获得的。信号样本可分别地用向量S,W,r和n表示。
反褶积的目的是从给定的一个观测结果S中恢复反射系数r。为了做 到这点，要先估算出W。在某些方法中它是与r同时估算的，而在另外一 些方法中却假设它为已知的。
同一时间上，与数据道的振幅有同样的水平。时窗不能太小，整道有3
至5个时窗即可。
(2)直接求反褶积输出道
从预测的观点出发，反褶积的输出是预测误差，根据（14）式我们
可以直接写出反滤波因子的具体形式：
（16）式可写成Toeplitz矩阵形式： （17）
（二） 地表一致性反褶积 以Robinson褶积模型为基础的反褶积处理仍是目前提高地震资料分
在地表一致性假设中，地震道被看成是震源算子，接受算子、反射 系数算子（按CDP排列）和与炮检距有关的算子的褶积。在对数-傅立叶 域中，褶积变成求和。对于给定的频率，其信号振幅的对数是震源、接 收、构造和与炮检距有关项的和。由于我们认为，对地表同一位置，滤 波作用与地震波的入射角无关，无论是浅、中、深层反射，其滤波作用 均相同。因此，我们把实现这种反滤波功能的方法，称为“地表一致性 反褶积”。
地震子波数学褶积模型： （18）
式中 —炮点坐标；—检波点坐标；—炮道的有效波； —第号震源脉冲；—该道接收点响应； —与偏移距有关的地层响应，该道半偏移距为； —该道的地层脉冲响应，即反射系数函数； 对上式求复谱
（19）
分解为振幅谱及相位谱：
（波是最小相位的，则只需考虑振幅谱。对振幅谱求对数：
当前地震资料处理解释已经基本实现了数据化、自动化，我国各大 解释公司、研究所、高等院校都已有了较为先进数字化处理软件，在处 理数字化的地震数据时表现出了很好的速度性和准确性。反褶积可分为 确定性反褶积和估计性反褶积两种。目前常用的反褶积有最小平方反褶 积、预测反褶积、同态反褶积、地表一致性反褶积、最大熵反褶积、变 模反褶积、Q反褶积等等；特殊的反褶积有Noah反褶积、最小信息反褶 积等。
输入（CMP道集或记录道集） 高通滤波（消除面波） 地形高程校正 叠加
振幅补偿（可选） 叠后去噪处理（可选）
反褶积 振幅均衡（可选） 去噪处理（可选）
偏移 叠加速度分析
滤波 动校正 显示动平衡 剩余静校正 输出 图1 地震资料数字处理流程图
（二）反褶积的局限性、国内外现状和发展趋势 1、反褶积的局限性
反褶积的目的就是为了分离子波和反射系数序列。子波就像无线电 中的载波，反射系数序列就像无线电中的声波。只有消除高频载波才能 提取声波。子波在地层中传播，携带着反射系数序列这种有用的地质信
息返回地面，只有消除子波才能恢复反射系数序列的本来面目。反射系 数序列中有波阻抗随时间变化的信息，这就提供了速度和密度随时间变 化的信息,随之就可得到地层、岩性及构造在地下中间分布的信息。在 有利条件下还可得到岩石孔隙率、渗透率、孔隙流体性质（油、气、 水）乃至地层压力的信息。反褶积提高了分辨率，拓展了频带，保持了 信噪比。
的偏导数，并令其等于零，就可得到最小平方反褶积求解反因子的基本
方程组：
（8）
式中，为子波的自相关函数；为期望输出与子波的互相关函数。
考虑脉冲反褶积的期望输出形式，方程组（8）右端的具体形式。
如果认为是物理可实现的，即当时，，就可以将表示成。于是就可以得
到：
（9）
如果认为反射系数序列是一个白噪序列，则可以得到：数据道的自
（2）频带的限制 地震道是限带的，而期望输出反射系数序列是宽带的，要从限
带输入中得出宽带输出是—个极大的矛盾，就是说要得到反射系数脉冲 序列几乎是不可能的。
（3）提高了分辨率降低了信噪比 实际上地震记录中不可能没有噪音。经反褶积后，只有在靠近反褶 积算子振幅负峰的很小频段内，噪音振幅谱才会降低，而在其两侧噪音 迅速提高，信噪比降低更多，情况恶化了。为了滤去信噪比恶化的频 段，要专门设计一个滤波算子。 （4）反褶积早被反演理论所超越 2、国内外现状及发展趋势 目前常用的反褶积有最小平方反褶积、预测反褶积、同态反褶积、 地表一致性反褶积、最大熵反褶积、变模反褶积、Q反褶积等等；特殊 的反褶积有Noah反褶积、最小信息反褶积、Kalman反褶积等。 反褶积方法有许多局限性和缺陷，这一方面限制了反褶积的发展， 另一方面也促使人们来弥补，使反褶积方法仍不断有所发展。发展的方 向大致有如下五个：一是对最小平方方法及变模法的前提进行修改，突 破已有的限制；二是改进具体的算法，使之更稳定，更唯一，更快速； 三是发展子波整形技术，直接为提高信噪比，保真度、一致性及分辨率 服务；四是提高子波估计水平及质量，为子波反褶积提供更确定的前 提；五是发展全新的概念，建立新的目标函数和判别准则。 基于上面反褶积发展的五个方向，反褶积的新进展有： （1）为摆脱相位限制，提出了混合相位最小平方反褶积。 （2）为摆脱反射系数序列是白噪的假设，提出了非白反褶积，约 瑟夫反褶积及蓝色补偿反褶积。 （3）为弥补地震道的限带性质，提出了限带脉冲反褶积及广义正 向迭代反褶积。 （4）为了提高估计反褶积算子中所需的自相关数据质量，提出了 约束迭代谱反褶积。 （5）为了改善最小熵反褶积的效果，除前述的一些方法外还有： 最佳滞后最小熵反褶积及包络敏感反褶积。 （6）为了提高预测反褶积的效果，提出了波动方程反褶积，速度