可控震源

可控震源

-

2023年2月27日发(作者：邮票知识)

可控震源零相位转换问题探讨

杨全斌;刘凤智;吴永国;欧荣生;张立军

【摘要】Thetransformationofzerophaseintosmallphaseinvibroseis

seismicdataprocessingisaproblemofcontroversy:Isitanecessary

process?whenwillisbeconducted?Whatisitsinfluenceonseismic

wavelet?Theauthorsanalyzedsmallphasetransforma-tionsofthe

autocorrelationwaveletsweepsignalofvibroseis,thesimulatedearth

hors

havefoundthatphasedriftphenomenonisinevitableinthephase

transformationprocess,andproposethemodificationmethod,whichcanto

horsaffirmed

thenecessityofphasetransformation,andputforwardtheoptimalorder

toconducttheprocessing.%当前对于可控震源零相位资料小相位转换存在较大

争议,是否为必需流程?何时进行转换?转换后对地震子波的影响如何?本文通过对可

控震源的扫描信号自相关子波、模拟大地吸收后可控震源的子波和实际可控震源采

集数据做小相位转换处理分析,发现相位转换过程中的相位漂移现象不可避免,并明

确给出修改方法,一定程度上解决了相位漂移的问题;文章肯定了相位转换的必要性,

同时给出处理次序.

【期刊名称】《物探与化探》

【年(卷),期】2018(042)003

【总页数】5页(P555-559)

【关键词】可控震源;零相位;小相位;相位漂移

【作者】杨全斌;刘凤智;吴永国;欧荣生;张立军

【作者单位】东方地球物理勘探有限责任公司青海物探处,甘肃敦煌736202;东方

地球物理勘探有限责任公司青海物探处,甘肃敦煌736202;东方地球物理勘探有限

责任公司青海物探处,甘肃敦煌736202;东方地球物理勘探有限责任公司青海物

探处,甘肃敦煌736202;东方地球物理勘探有限责任公司青海物探处,甘肃敦煌

736202

【正文语种】中文

【中图分类】P631.4

0引言

随着高效采集的发展及人们对环保的要求越来越高，可控震源作为一种高效、环保、

相对经济的勘探工具应用越来越广泛，零相位资料的小相位转换处理是可控震源资

料处理中必不可少的部分。在可控震源资料的处理过程中，需要经常考虑到可控震

源资料小相位转换的问题。在一些资料书籍和处理人员的认识中，对于零相位资料

小相位转换存在较大争议，比如：用设计扫描信号对经过大地吸收后的地震记录做

互相关，地震资料的相位还是不是零相位；针对设计扫描信号的小相位转换因子，

对经过大地吸收后的地震记录做小相位转换有没有影响；在资料处理环节中什么时

候应用小相位转换最合理等［1-2］。文中以柴达木盆地平台地区为例，应用

GeoEast地震资料处理系统，针对以上这些问题进行了分析，提出了解决方案，

获得了较好的效果。

1可控震源记录的相位分析

通常，在不考虑大地吸收的情况下，人们认为可控震源相关后记录是零相位的，但

由于受大地吸收的影响，设计的扫描信号和接收到的经过大地吸收后的扫描信号互

相关子波还是不是零相位子波？为了验证这个问题，选取了一个可控震源相关前的

单炮记录数据及其扫描信号，通过子波模拟的方式来分析大地吸收对可控震源资料

子波的影响：图1a是设计的可控震源扫描信号，图1c是该扫描信号的频谱，图

1d是使用该扫描信号激发的相关前单炮的频谱；那么图1d频谱的包络线可以近

似的认为是频率域的大地吸收滤波曲线；在本文实例中，应用俞氏10～50Hz带

通滤波器曲线模拟它。用图1d模拟的大地吸收滤波曲线对扫描信号图1a滤波，

可以得到经过大地吸收后的扫描信号图1b。

图2显示的是大地吸收滤波对相关子波相位的影响。图2a是图1a可控震源扫描

信号的自相关子波，图2c是图1b模拟的可控震源扫描信号经大地吸收后的自相

关子波，毫无疑问，它们是一个以时间0轴对称的零相位子波；图2b是图1a与

图1b的互相关子波，它同样是一个以时间0轴对称的零相位子波。因此，我们认

为对可控震源相关前记录用设计扫描信号相关后，其相位仍然是零相位的。

图1模拟大地吸收对可控震源扫描信号的影响

图2模拟大地吸收对可控震源相关子波的影响

2可控震源资料的小相位转换

2．1零相位资料转换成小相位的必要性

根据前面的分析，可控震源相关后记录是零相位，因此存在以下问题：

1）炸药、可控震源资料联合处理时，两种激发方式获得的地震资料相位不统一，

影响叠加效果，降低资料信噪比。

2）最小平方脉冲、预测反褶积算法要求输入地震子波是最小相位的，而最小平方

预测反褶积基本是目前常规处理中不可或缺的处理方法。

2．2零相位资料的小相位转换

目前对可控震源资料小相位转换的常用办法为：先根据理论扫描信号求出它的自相

关子波；然后通过多项式求根法、最小平方法、希尔伯特变换法等求出自相关子波

的小相位转换因子（即小相位转换因子是根据设计信号计算获得的），再对可控震

源相关后记录与小相位转换因子褶积，得到小相位转换后的可控震源记录。

图3不同频宽子波小相位转换结果

从图3可以看出：转换后的小相位子波主瓣不如零相位子波突出，能量衰减缓慢，

分辨率变低；在实际资料处理中，我们是以理论的扫描信号计算小相位转换因子，

由于大地的滤波作用，造成实际接收到的信号振幅谱变窄，应用理论扫描信号求得

的转换因子进行小相位转换后，信号起跳点向前漂移，吸收越多，向前的时移越大

（图2f、图2d相对于图2b向前漂移），因此，在实际资料处理中，浅层由于受

大地吸收影响小，地震资料经过小相位转换后向前漂移少，越到深层，子波受大地

吸收影响越大，子波向前漂移现象越明显。而理论上，零相位子波经过正确的小相

位转换后起跳点位于0时间处。

图4是柴达木盆地平台地区可控震源资料小相位转换前后叠加剖面实例：对于

0．9～1s之间的基底低频反射同相轴来看，小相位转换后同相轴向前时移；

0．8～0．9s处的同相轴小相位转换后变弱，这是由于小相位转换后子波主瓣没

有零相位子波主瓣突出造成的；在不考虑后续处理的情况下，在分辨率方面，零相

位可控震源资料小相位转换后，分辨率并未提高。

图4小相位转换前后叠加剖面对比

3小相位转换中问题讨论

3．1小相位转换对去噪的影响

在对可控震源资料的处理过程中，大家认为可控震源资料是零相位的，所以很多人

一开始资料处理就对其进行小相位转换处理［2］，但从模拟的零相位子波小相位

转换的过程中可以清晰看到，零相位子波主瓣比其小相位转换后的子波主瓣更加突

出（见图3）。因此，在叠前去噪处理中零相位子波比其转换后的小相位子波更有

优势，例如在线性干扰压制方面零相位资料处理效果更佳（见图5）。

图5小相位转换与线性噪声压制先后顺序对比单炮

3．2小相位转换对反褶积的影响

尽管人们认为零相位的资料不满足最小平方预测反褶积要求（输入时子波是最小相

位的假设），如果对零相位资料直接进行反褶积，造成反褶积后输出的子波不是我

们期望的最小相位子波。但很多时候对零相位资料处理时仍然使用预测反褶积，对

零相位资料使用预测反褶积后在压缩子波方面效果仍然很好（见图6），反褶积后

子波主瓣更尖锐，表现在地震剖面上就是地震反射同相轴明显变细。

图6零相位子波小相位转换、反褶积效果对比

3．3小相位转换与反褶积的应用次序问题

在小相位转换与反褶积使用次序方面，一般认为预测反褶积是在最小相位的假设下

才正确，对于零相位资料应当先小相位化，然后反褶积，但通过图6c、图6e的对

比，我们认识到先做小相位化，再做反褶积和先做反褶积再应用小相位化，其结果

一致。这并不是偶然现象，对于小相位转换与反褶积的应用不分先后次序这一特点，

可以通过公式证明。最小平方法预测反褶积的数学模如下：

输入信号：x（t）；设预测滤波因子：c（t）＝［c（0），c（1），…，c

（m）］；期望输出：x（t＋τ）（τ＞0）。通过最小平方法，我们能推出预测反

褶积的法方程系（推导过程见文献［2-5］）：

解上述方程组，获得反褶积因子c（t），将方程写成矩阵形式：

以上方程的系数矩阵左端由x（t）的自相关函数组成，中间的系数矩阵为预测反

褶积因子，右边的系数矩阵是x（t）与期望输出子波的互相关，因此上述方程组

左端及右端不管是零相位子波还是与其对应的小相位子波，根据相关运算的原理，

它们的值是相同的；因此，一个零相位子波与其对应的小相位子波求出的反褶积因

子c（t）是一样的，所以，对一个零相位子波直接进行预测反褶积，其结果与零

相位子波小相位转换后做预测反褶积，只差一个反褶积前零相位到小相位的转换因

子，对零相位子波反褶积后的结果应用其反褶积前对应的小相位转换因子进行相位

校正，结果是相同的。图7是实际资料的处理效果对比。

图7小相位转换与预测反褶积应用先后顺序对比

对比零相位与小相位转换后做预测反褶积的结果（图6中的c和d的对比），零

相位子波小相位转换后反褶积的子波更突出，旁瓣更小，衰减速度更快，因此分辨

率更高；预测反褶积对零相位子波小相位转换中造成的主瓣不够突出得到改善，转

换过程中产生的强尾波由于预测反褶积对鸣震压制效果而减弱，而零相位子波直接

进行反褶积后只是子波主瓣变的尖锐，而分辨率不如对其小相位转换后再反褶积的

效果好。

4结论

通过前面的分析，对于可控震源资料小相位转换过程中存在的问题得出以下结论：

1）可控震源相关后记录是零相位的，但对可控震源记录使用由扫描信号求取的小

相位转换因子进行小相位转换存在问题，由于大地的吸收作用造成转换后的小相位

子波向前漂移；可以尝试在小相位转换之前应用反Q补偿等拓频手段，减轻大地

吸收对子波的影响，使之与设计子波频宽接近后再转换，以减少子波向前漂移，在

有井资料的情况下需要进一步的井控处理，以提高成像的精度。

2）小相位转换在最小平方预测反褶积前后做效果一致，因此不要在资料处理一开

始就做零相位资料小相位转换，可以先在零相位资料上做去噪、补偿等工作，然后

再做零相位子波小相位转换处理。

3）零相位资料小相位转换后，小相位子波不如零相位子波好，但零相位资料和转

换后的小相位资料串用同样的预测反褶积之后，小相位资料要比零相位资料好，因

此，尽管用扫描信号求出的小相位转换因子对可控震源资料进行小相位转换时会产

生相位漂移现象，但对可控震源资料进行小相位转换处理仍是必要的。

参考文献：

［1］李庆忠．走向精确勘探的道路［M］．北京：石油工业出版社，1993．

［2］熊翥．地震数据处理应用技术［M］．北京：石油工业出版社，2008．

［3］渥·伊尔马兹．地震资料分析［M］．刘怀山等，译．北京：石油工业出版

社，2006．

［4］高少武，魏庚雨，赵波，等．可控震源地震数据子波最小相位化方法

［J］．石油地球物理勘探，2009，44（6）：685-689．

［5］皮科克·K·L，特雷特尔·斯文．预测反褶积的理论与实践［J］．石油物探，

1974，2（1）：87-104．

［6］OsvaldodeOliveiraDuarte．含噪声数据的预测反褶积［C］／／美国勘

探地球物理学家学会第62届年会论文集，1992．

［7］闵成花．零相位校正和子波反褶积［J］．油气地球物理，2004，2（2）：

13-16．

［8］李国发，李观寿，张立勤．零相位改进型最小熵反褶积方法［J］．石油地

球物理勘探，1996，31（3）：359-366．

［9］宋锡熊．一种零相位反褶积方法［J］．石油地球物理勘探，1984，6（3）：

266-271．

［10］李强，何晓松，王立．一种实现井约束零相位化的处理方案［J］．石油地

球物理勘探，2009，44（S1）：76-78．

［11］吕景英，蒋进勇．地震子波零相位化改造的应用［J］．勘探地球物理进展，

2003，26（2）：121-123．

［12］舍恩伯格·M．最小相位和零相位信号的分辩率比较［J］．石油地球物理勘

探，1976，11（S1）：34．