子波零相位化处理流程.pdf

子波零相位化处理流程

由于海上气枪震源产生的地震子波是混合相位子波，而反褶积的要求是最小相位子波

（零相位子波的效果最好），如果是混合相位子波达不到十分理想的效果，所以理论上应该

先把地震子波做最小相位化处理（零相位化处理的效果更好）再做反褶积。

首先打开Waveval或者Matcal，选择不建工区，进去之后的界面如下：

选择上图红圈所示的Load选件,按下图所示选择参数，选定远场子波0000.dag，OK即可。

选择主菜单的“×”号，把起跳设置成负值。按下图所示选择参数。

下图就是子波翻转之后的对比情况，起跳由原来的正值变成了负值。

由上图翻转后的子波可以看出，第一个负的非零采样点位于236ms处，要把它移动到0ms

处，选择上图红圈所示的Time-Shift选件，按下图的参数填写。

下图就是时移后的子波，该子波由于只包含震源鬼波，所以还要把检波器鬼波加进去。

选择下图红圈所示的WaveletModling选件，再选择下图中的GhostWavelet选件。

按下图所示的参数模拟接收器鬼波。

下图就是模拟出来的接收器鬼波，是个零相位子波。

把模拟出来的接收器鬼波用褶积的方式加到翻转时移之后的远场子波中，得到的就是实

际地震资料的近场子波。选择上图红圈所示的convolution选件把翻转时移后的子波和鬼波

进行褶积计算。参数如下图所示：

下图就是加接收器鬼波前的子波波形和加子波之后的子波波形。

结下来要把加入接收器鬼波之后的子波（该子波既包含了震源鬼波由包含了接收器鬼波）

进行去气泡（de-bulling）处理。选择上图红圈所示的WaveletModling选件，然后选择下图

的PredictiveDeconvolution选件。

之后按下图所示的参数填上即可。

上图子波长度必须包含最小的旗袍周期，gap值不大于100，最小也要大于第二个零点

位置。下图就是预测反褶积压制鬼波之后的效果。

同时这一步也产生了一个de-buling因子，波形如下：

下图是将求取的去鬼波因子作用到原始子波上进行QC的情况。红色是去鬼波后的子波。

接下来就是求零相位化子波了。选择上图红圈所示的OperatorCalc选件，然后选择

OperatorDesion选件，按下图所示的参数填写即可。

下图就是产生的零相位化因子的波形，把这个零相位化子波和之前产生的那个去气泡子

波（de-buling）进行褶积，得到的就是最终要求的子波了。

下图是将求得的零相位化子波作用到原始的子波上进行QC的情况。图上红色是零相位化的

子波，黑色的是原始子波。

下面是将两个子波进行褶积求取最终子波的褶积计算。

下图就把最终得到的子波以滤波器的形式保存。

到此为止，要求的子波已经完成，下面是一些效果图展示。

最终因子QC

最后因子褶积

最后因子QC

保存最终因子

效果：

前两步：乘-1，时移

加接收器鬼波：

去气泡：

零相位化因子

最后QC：

翻转时移