地震勘探原理之反褶积

基础信息

地震资料数据常规处理流程

• 反褶积：提高地震资料分辨率，压制部分干扰波
• 水平叠加：提高地震资料信噪比
• 偏移成像：提高地震资料保真度

所谓褶积的概念相对简单，其实褶积与卷积在数学上是等价的。地球物理当中的褶积主要是指代地震波在地下介质当中传播的过程。换句话说，就是震源经过褶积的过程，得到了检波器记录到的地震记录。
而反褶积是一种处理方法，由于这样子直接得到的地震记录有很多缺点，反褶积可以有效处理这个问题。

震源产生尖脉冲在地下介质传播，经过反射界面返回地面。在理想情况下，如果是一系列尖脉冲，每一个尖脉冲对应一个反射界面，就可以确定地层界面。可以用公式表示：

$X(t)=W(t)r(t)$

即地震记录$X(t)$等于震源脉冲$W(t)$乘以反射系数$r(t)$。
经过如下一系列的处理，震源发出的尖脉冲变成了地震子波。

由上图可知，地震记录是许多反射波叠加的结果。地震记录$W(t)$可以看作地震子波$X(t)$与反射系数$r(t)$褶积的结果。
下图就是“褶积”的具体运算公式。

如果考虑干扰波的影响，那么实际的地震记录就满足如下图的公式。

引出问题

在地震记录上，一个反射界面的反射波，对应着一个延续时长为几十毫秒的波形。由于地下反射界面一般是由相距几米或者几十米的底层交错组成，他们发出的反射波到达时间一般相差几毫秒或几十毫秒，所以它们在地震记录上发生干涉，难以分开。
所以为了提高地震分辨率，希望得到的地震记录中每个反射波都表现为一个窄的脉冲。

脉冲的性质

• 脉冲的强弱与界面反射系数大小成正比
• 脉冲的极性对应着界面反射系数的正负

反褶积的过程

在Robinson反褶积模型中，最重要的是求出地震反子波（又名反滤波因子或反褶积因子）。方法介绍如下：

然而，在反褶积处理过程中，地震子波的求取精度对于反褶积效果影响较大，目前求取地震子波的方法如下：

反褶积方法如下：

常用反褶积介绍

脉冲反褶积

（不够具体）
一句话概括：脉冲反褶积是在地震数据处理时将实际地震记录期望改造为与反射界面相对应的尖脉冲的处理过程。
在地震勘探中，地震波在地层内传播时会遇到不同的介质分界面，发生反射和折射等现象。这些反射波会记录在地震勘探的记录中。由于地震波在地下传播时存在散射和衰减等效应，使得地震记录中的反射波具有较大的宽度，降低了分辨率。因此，通过这种办法，将有宽度的子波尽可能压缩成脉冲，从而将地下连续的反射系数序列转换成脉冲序列。这样可以提高地震记录的分辨率，更准确地估计地下反射界面的反射系数。
脉冲反褶积是在最小平方误差准则下，利用地震子波和反射系数伪随机过程的假设条件来求解反褶积因子。

反射系数序列是一个随机过程，意味着地震记录具有地震子波的特征，也就是地震记录和地震子波自相关和振幅谱相似。这就允许使用地震记录的自相关可以代替未知的地震子波的自相关。
脉冲反褶积的目标：

• 将输出的振幅谱展平。即将地震记录中的反射系数序列转换成脉冲序列，以便更准确地分析地下地质结构等信息。
• 使输出的波形逼近零相位的尖脉冲。（一种理想化的脉冲信号，其相位延迟为零，即在时间上没有相位延迟。）
  具体效果可以参考这个图
  
  在下图例子中，共炮点地震记录经过反褶积后分辨率获得极大提升。
  

预测反褶积

设计一个算子，对已知物理量的过去值和现在值进行处理，获得未来某个时刻的值，这个过程就是预测反褶积。使用的过去值和现在值的时差叫做预测步长。
预测反褶积是采用预测误差滤波器对地震道进行滤波的过程。由于可以适当的选择预测步长，有针对性的压制和消除在该尺度上的重复振动。所以经常用于消除海上鸣震混响，消除多次波等可预测的干扰波。

由于预测反褶积希望输出是最小相位值（最小相位延迟），因此它也能在地震数据处理中压缩地震子波，提高地震分辨率。