CN105093328B - 一种滑动扫描谐波压制方法及装置 - Google Patents
一种滑动扫描谐波压制方法及装置 Download PDFInfo
- Publication number
- CN105093328B CN105093328B CN201510504580.2A CN201510504580A CN105093328B CN 105093328 B CN105093328 B CN 105093328B CN 201510504580 A CN201510504580 A CN 201510504580A CN 105093328 B CN105093328 B CN 105093328B
- Authority
- CN
- China
- Prior art keywords
- mrow
- correlation
- track data
- harmonic wave
- imf components
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及石油勘探技术领域，尤其涉及一种滑动扫描谐波压制方法及装置，该方法包括以下步骤：获取滑动扫描相关后道数据，并拾取所述相关后道数据的初至时间；将所述相关后道数据进行互补总体平均的经验模态分解CEEMD，得到一组本征模函数IMF分量；分析谐波在所述IMF分量中的分布情况，确定含有谐波的IMF分量和不含有谐波的IMF分量；将所述含有谐波的IMF分量从所述初至时间开始进行谐波衰减处理；利用谐波衰减处理后的IMF分量以及所述不含有谐波的IMF分量进行重构得到去谐波后的相关后道数据。本发明实现了在无需震源力信号的情况下对滑动扫描谐波的有效压制，提高了地震数据的信噪比。
Description
技术领域
本发明涉及石油勘探技术领域，尤其涉及一种滑动扫描谐波压制方法及装置。
背景技术
目前，在石油地震勘探中，滑动扫描技术是国际上可控震源高效采集技术中成熟、常用的技术之一。在该技术中，谐波干扰是客观存在的，并且一般表现为后一炮中的谐波畸变对前一炮产生影响，从而降低了地震资料的品质，因此在处理时需要专门的压制技术。
然而，现有可控震源谐波压制方法(如模型法、滤波法等)都需要震源力信号来计算一个滤波因子，从而估计谐波分量，再从原始记录减去谐波分量，达到谐波衰减的目的。震源力信号是震源实际传到地下的信号，通常震源力信号的记录需要额外的设备和手段，在高效采集施工中，常常会有某些震次的震源力信号没记录上的情况。因此，如何在没有震源力信号或震源力信号不完整的情况下实现对可控震源滑动扫描谐波的压制，以提高地震数据信噪比是目前石油勘探技术领域亟需解决的技术问题。
发明内容
本申请实施例的目的在于提供一种滑动扫描谐波压制方法及装置，通过基于互补总体平均的经验模态分解(Complementary Ensemble Empirical Mode Decomposition,CEEMD)的分解手段，在无需震源力信号的情况下实现对滑动扫描谐波的压制，提高地震数据信噪比。
为达到上述目的，一方面，本申请实施例提供了一种滑动扫描谐波压制方法，包括：
获取滑动扫描相关后道数据，并拾取所述相关后道数据的初至时间；
将所述相关后道数据进行CEEMD分解，得到一组本征模函数(Intrinsic ModeFunction,IMF)分量；
分析谐波在所述IMF分量中的分布情况，确定含有谐波的IMF分量和不含有谐波的IMF分量；
将所述含有谐波的IMF分量从所述初至时间开始进行谐波衰减处理；
利用谐波衰减处理后的IMF分量以及所述不含有谐波的IMF分量进行重构得到去谐波后的相关后道数据。
另一方面，本申请实施例还提供一种滑动扫描谐波压制装置，包括：
获取模块，用于获取滑动扫描相关后道数据，并拾取所述相关后道数据的初至时间；
CEEMD分解模块，用于将所述相关后道数据进行CEEMD分解，得到一组IMF分量；
谐波分析模块，用于分析谐波在所述IMF分量中的分布情况，确定含有谐波的IMF分量和不含有谐波的IMF分量；
谐波衰减模块，用于将所述含有谐波的IMF分量从所述初至时间开始进行谐波衰减处理；
重构模块，用于利用谐波衰减处理后的IMF分量以及所述不含有谐波的IMF分量进行重构得到去谐波后的相关后道数据。
本申请实施例是基于CEEMD分解的一种滑动扫描谐波压制方法，通过对相关后道数据进行CEEMD分解得到IMF分量，随后分析谐波在IMF分量重的分布情况并对含有谐波的IMF分量进行衰减处理，再利用衰减后的含谐波的IMF分量和不含谐波的IMF分量重构了去谐波后的道数据，从而实现了在无需震源力信号的情况下对滑动扫描谐波的有效压制，提高了地震数据的信噪比。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1是本申请实施例的滑动扫描谐波压制方法流程图；
图2是本申请实施例的滑动扫描谐波压制装置的结构框图；
图3是本申请一实施例的滑动扫描相关后的原始单道数据示意图；
图4是本申请一实施例的原始单道数据对数变换前后对比示意图；
图5是本申请一实施例的原始单道数据谐波压制前后的对比示意图
图6是本申请一实施例的谐波压制前的单炮记录示意图；
图7是本申请一实施例的谐波压制后的单炮记录示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
下面结合附图，对本申请实施例的具体实施方式作进一步的详细说明。
参考图1所示，本申请实施例的滑动扫描谐波压制方法，包括以下步骤：
S1、获取滑动扫描相关后道数据，并拾取所述相关后道数据初至时间。
其中，本申请实施例的相关后道数据指地震仪器接收到的数据与可控震源参考信号互相关后的共炮点道集地震数据；初至时间指的是滑动扫描相关后道数据初至波达到的时间。如图3所示为获取的原始单道数据X(t)，其中X(t)有N个样点值，N＞1024，在该图中，纵坐标表示信号强度，横坐标表示为时间；由图可以看出中深层的谐波干扰非常严重。
S2、将所述相关后道数据进行CEEMD分解，得到一组IMF分量。具体的，CEEMD分解包括：
向所述相关后道数据X^(t)中加入Nit组辅助白噪声，辅助白噪声是以正、负对的方式加入，从而生成两套集合IMF：
其中，X^(t)为相关后道数据；Noise(t)为辅助噪声，其均值为0，标准差为1；X^1(t)，X^2(t)分别为加入正负成对的噪声后的信号，从而得到集合信号的个数为2Nit。
将集合中的每一个信号进行经验模态分解(Empirical Mode Decomposition,EMD)，每个信号得到一组IMF分量；
通过多组分量组合的方式得到分解结果：
其中，Cij表示EMD分解后第i个信号的第j个IMF分量，Cj表示CEEMD分解最终得到的第j个IMF分量。
S3、分析谐波在所述IMF分量中的分布情况，确定含有谐波的IMF分量和不含有谐波的IMF分量。本申请实施例中，根据IMF的频率从高到低的分布特征及可控震源勘探中谐波的频带范围，第一个IMF分量中的主要成分是在CEEMD分解过程中产生的随机噪声及部分高频谐波成分，其余谐波主要分布在第2、3、4阶IMF中，即C2，C3，C4中，大于9阶的IMF分量为CEEMD分解过程中的低频假象。
S4、将所述含有谐波的IMF分量从所述初至时间开始进行谐波衰减处理。本申请实施例中谐波衰减处理目的是保护浅层高频信息，衰减中深层的谐波。具体的，可以将含有谐波的IMF分量从初至时间开始乘以布莱克曼窗函数，具体公式为：
其中：
Winblackman(t)为布莱克曼窗函数，长度L由IMF分量中谐波开始出现的位置确定，一般为500-1500个样点值；Cm(t)为含有谐波的IMF分量，m＝2，3，4；C^m(t)为谐波衰减处理后的IMF分量。
S5、利用谐波衰减处理后的IMF分量以及所述不含有谐波的IMF分量进行重构得到去谐波后的相关后道数据。本申请实施例中根据下述公式进行重构得到去谐波后的相关后道数据：
C^m(t)为谐波衰减处理后的IMF分量；Ci(t)为不含有谐波的IMF分量，其中n＝log2(N)，N为样点值个数；X^denoise(t)为去谐波后的相关后道数据。
在步骤S2之前，本申请实施例还可以包括将所述相关后道数据进行标准化预处理，对应的，步骤S2中所述相关后道数据为标准化预处理后的相关后道数据。本申请实施例中所述的标准化预处理具体包括：
将所述相关后道数据进行对数变换；
计算对数变换后的相关后道数据的标准差；
将所述对数变换后的相关后道数据除以所述对数变换后的相关后道数据的标准差获得标准化预处理后的相关后道数据。
其中，本申请实施例中对数变换目的在于放大弱信号而不增加人为的假象，有利于互CEEMD分解过程中保护弱信号；同时标准化预处理后的相关后道数据的标准差和CEEMD过程中加入的辅助白噪声标准差都为1，从而有利于辅助白噪声的加入。如图4所示为原始单道数据对数变换前后的数据对比，从图中可以看出对数变换增强了弱信号。所述对数变换计算公式为：
Xlog(t)＝sign[slog(p+|X(t)|)]
其中：
其中，X(t)表示相关后道数据，sign和p表示条件系数，Xlog(t)表示对数变换后的相关后道数据，M表示允许输出的最大值，I表示输入信号X(t)中的最大值。
在步骤S5之后，本申请实施例还可以包括将所述得到的去谐波后的道数据乘以所述对数变换后的相关后道数据的标准差，从而将其转换到原始数据的能量等级，并得到最终的去谐波后的相关后道数据。如图5所述为最终的谐波压制后的道数据与原始单道数据的对比，可以看出本申请实施例对谐波进行了有效压制。图6和图7为谐波压制前后的炮集数据图，可以看到经过谐波压制处理后的单炮记录信噪比比原始单炮记录显著提高，本申请实施例的谐波压制效果是非常显著的。
本申请实施例中通过对相关后道数据标准化预处理手段，有效保护了弱信号而不增加人为的假象，然后再利用CEEMD分解对含有谐波的IMF分量进行衰减处理，再利用衰减后的含谐波的IMF分量和不含谐波的IMF分量重构去谐波后的相关后道数据，从而实现了在无需震源力信号的情况下对滑动扫描谐波的有效压制，提高了地震数据的信噪比。
参考图2所示，与上述滑动扫描谐波压制方法对应，本申请实施例的滑动扫描谐波压制装置，包括：
获取模块21，用于获取滑动扫描相关后道数据，并拾取所述相关后道数据初至时间。
其中，本申请实施例的相关后道数据指地震仪器接收到的数据与可控震源参考信号互相关后的共炮点道集地震数据；初至时间指的是滑动扫描相关后道数据初至波达到的时间。
CEEMD分解模块22，用于将所述相关后道数据进行CEEMD分解，得到一组IMF分量。具体的，CEEMD分解模块包括：
白噪声加入子模块，用于向标准化预处理后的相关后道数据X^(t)中加入Nit组辅助白噪声，辅助白噪声是以正、负对的方式加入，从而生成两套集合IMF：
其中，X^(t)为标准化处理后的相关后道数据；Noise(t)为辅助噪声，其均值为0，标准差为1；X^1(t)，X^2(t)分别为加入正负成对的噪声后的信号，从而得到集合信号的个数为2Nit。
EMD分解子模块，用于将集合中的每一个信号进行EMD分解，每个信号得到一组IMF分量；
多组分量组合子模块，用于将EMD分解的IMF分量通过多组分量组合的方式得到分解结果：
其中，Cij表示EMD分解后第i个信号的第j个IMF分量，Cj表示CEEMD分解最终得到的第j个IMF分量。
谐波分析模块23，用于分析谐波在所述IMF分量中的分布情况，确定含有谐波的IMF分量和不含有谐波的IMF分量。本申请实施例中，根据IMF的频率从高到低的分布特征及可控震源勘探中谐波的频带范围，第一个IMF分量中的主要成分是在CEEMD分解过程中产生的随机噪声及部分高频谐波成分，其余谐波主要分布在第2、3、4阶IMF中，即C2，C3，C4中，大于9阶的IMF分量为CEEMD分解过程中的低频假象。
谐波衰减模块24，用于将所述含有谐波的IMF分量从所述初至时间开始进行谐波衰减处理。本申请实施例中谐波衰减处理目的是保护浅层高频信息，衰减中深层的谐波。具体的，可以将含有谐波的IMF分量从所述初至时间开始乘以布莱克曼窗函数，具体公式为：
其中：
Winblackman(t)为布莱克曼窗函数，长度L由IMF分量中谐波开始出现的位置确定，一般为500-1500个样点值；Cm(t)为含有谐波的IMF分量，m＝2，3，4；C^m(t)为谐波衰减处理后的IMF分量。
重构模块25，用于利用谐波衰减处理后的IMF分量以及所述不含有谐波的IMF分量进行重构得到去谐波后的相关后道数据。本申请实施例中根据下述公式进行重构得到去谐波后的相关后道数据：
C^m(t)为谐波衰减处理后的IMF分量；Ci(t)为不含有谐波的IMF分量，其中n＝log2(N)，N为样点值个数；X^denoise(t)为去谐波后的相关后道数据。
本申请实施例还可以包括标准化预处理模块，用于在将所述相关后道数据进行CEEMD分解之前，对所述相关后道数据进行标准化预处理；对应的，CEEMD分解模块中的相关后道数据为标准化预处理后的相关后道数据。本申请实施例中，所述标准化预处理模块具体包括：
对数变换子模块，用于将所述相关后道数据进行对数变换；
标准差子模块，用于计算对数变换后的相关后道数据的标准差；
标准化子模块，用于将所述对数变换后的相关后道数据除以所述对数变换后的相关后道数据的标准差获得标准化预处理后的相关后道数据。
其中，本申请实施例中对数变换目的在于放大弱信号而不增加人为的假象，有利于互CEEMD分解过程中保护弱信号；同时标准化预处理后的相关后道数据的标准差和CEEMD过程中加入的辅助白噪声标准差都为1，从而有利于辅助白噪声的加入。对数变换具体公式为：
Xlog(t)＝sign[slog(p+|X(t)|)]
其中：
其中，X(t)表示相关后道数据，sign和p表示条件系数，Xlog(t)表示对数变换后的相关后道数据，M表示允许输出的最大值，I表示输入信号X(t)中的最大值。
本申请实施例中通过对相关后道数据标准化预处理手段，有效保护了弱信号而不增加人为的假象，然后再利用CEEMD分解对含有谐波的IMF分量进行衰减处理，再利用衰减后的含谐波的IMF分量和不含谐波的IMF分量重构去谐波后的相关后道数据，从而实现了在无需震源力信号的情况下对滑动扫描谐波的有效压制，提高了地震数据的信噪比。
本领域技术人员还可以了解到本申请实施例列出的各种说明性逻辑块、单元和步骤可以通过硬件、软件或两者的结合来实现。至于是通过硬件还是软件来实现取决于特定的应用和整个系统的设计要求。本领域技术人员可以对于每种特定的应用，可以使用各种方法实现所述的功能，但这种实现不应被理解为超出本申请实施例保护的范围。
本申请实施例中所描述的各种说明性的逻辑块，或单元都可以通过通用处理器，数字信号处理器，专用集成电路(ASIC)，现场可编程门阵列或其它可编程逻辑装置，离散门或晶体管逻辑，离散硬件部件，或上述任何组合的设计来实现或操作所描述的功能。通用处理器可以为微处理器，可选地，该通用处理器也可以为任何传统的处理器、控制器、微控制器或状态机。处理器也可以通过计算装置的组合来实现，例如数字信号处理器和微处理器，多个微处理器，一个或多个微处理器联合一个数字信号处理器核，或任何其它类似的配置来实现。
本申请实施例中所描述的方法或算法的步骤可以直接嵌入硬件、处理器执行的软件模块、或者这两者的结合。软件模块可以存储于RAM存储器、闪存、ROM存储器、EPROM存储器、EEPROM存储器、寄存器、硬盘、可移动磁盘、CD-ROM或本领域中其它任意形式的存储媒介中。示例性地，存储媒介可以与处理器连接，以使得处理器可以从存储媒介中读取信息，并可以向存储媒介存写信息。可选地，存储媒介还可以集成到处理器中。处理器和存储媒介可以设置于ASIC中，ASIC可以设置于用户终端中。可选地，处理器和存储媒介也可以设置于用户终端中的不同的部件中。
在一个或多个示例性的设计中，本申请实施例所描述的上述功能可以在硬件、软件、固件或这三者的任意组合来实现。如果在软件中实现，这些功能可以存储与电脑可读的媒介上，或以一个或多个指令或代码形式传输于电脑可读的媒介上。电脑可读媒介包括电脑存储媒介和便于使得让电脑程序从一个地方转移到其它地方的通信媒介。存储媒介可以是任何通用或特殊电脑可以接入访问的可用媒体。例如，这样的电脑可读媒体可以包括但不限于RAM、ROM、EEPROM、CD-ROM或其它光盘存储、磁盘存储或其它磁性存储装置，或其它任何可以用于承载或存储以指令或数据结构和其它可被通用或特殊电脑、或通用或特殊处理器读取形式的程序代码的媒介。此外，任何连接都可以被适当地定义为电脑可读媒介，例如，如果软件是从一个网站站点、服务器或其它远程资源通过一个同轴电缆、光纤电缆、双绞线、数字用户线(DSL)或以例如红外、无线和微波等无线方式传输的也被包含在所定义的电脑可读媒介中。所述的碟片(disk)和磁盘(disc)包括压缩磁盘、镭射盘、光盘、DVD、软盘和蓝光光盘，磁盘通常以磁性复制数据，而碟片通常以激光进行光学复制数据。上述的组合也可以包含在电脑可读媒介中。
以上所述的具体实施例，对本申请的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本申请实施例的具体实施例而已，并不用于限定本申请的保护范围，凡在本申请的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。
Claims (8)
1.一种滑动扫描谐波压制方法，其特征在于，包括以下步骤：
获取滑动扫描相关后道数据，并拾取所述相关后道数据的初至时间；
将所述相关后道数据进行互补总体平均的经验模态分解CEEMD，得到一组本征模函数IMF分量；
分析谐波在所述IMF分量中的分布情况，确定含有谐波的IMF分量和不含有谐波的IMF分量；
将所述含有谐波的IMF分量从所述初至时间开始进行谐波衰减处理；
利用谐波衰减处理后的IMF分量以及所述不含有谐波的IMF分量进行重构得到去谐波后的相关后道数据；其中，
在所述将相关后道数据进行互补总体平均的经验模态分解CEEMD，得到一组本征模函数IMF分量之前，还包括：
将所述相关后道数据进行标准化预处理；
对应的，所述将相关后道数据进行互补总体平均的经验模态分解CEEMD中，所述相关后道数据为标准化预处理后的相关后道数据；
所述将相关后道数据进行标准化预处理，具体包括：
将所述相关后道数据进行对数变换；
计算对数变换后的相关后道数据的标准差；
将所述对数变换后的相关后道数据除以所述对数变换后的相关后道数据的标准差获得标准化预处理后的相关后道数据。
2.根据权利要求1所述的滑动扫描谐波压制方法，其特征在于，所述将含有谐波的IMF分量从所述初至时间开始进行谐波衰减处理，具体包括：
将所述含有谐波的IMF分量从所述初至时间开始乘以布莱克曼窗函数。
3.根据权利要求1所述的滑动扫描谐波压制方法，其特征在于，所述利用谐波衰减处理后的IMF分量以及所述不含有谐波的IMF分量进行重构得到去谐波后的相关后道数据的计算公式为：
<mrow>
<msub>
<msup>
<mi>X</mi>
<mo>^</mo>
</msup>
<mrow>
<mi>d</mi>
<mi>e</mi>
<mi>n</mi>
<mi>o</mi>
<mi>i</mi>
<mi>s</mi>
<mi>e</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>2</mn>
</mrow>
<mn>4</mn>
</munderover>
<msubsup>
<mi>C</mi>
<mi>m</mi>
<mo>^</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>5</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>C</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
其中，C^m(t)为谐波衰减处理后的IMF分量；Ci(t)为不含有谐波的IMF分量；X^denoise(t)为去谐波后的相关后道数据，n＝log2(N)，N为样点值个数。
4.根据权利要求1所述的滑动扫描谐波压制方法，其特征在于，所述将相关后道数据进行对数变换的具体计算公式为：
Xlog(t)＝sign[slog(p+|X(t)|)]
其中：
<mrow>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
<mo>...</mo>
<mo>...</mo>
<mi>X</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>></mo>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
<mo>...</mo>
<mo>...</mo>
<mi>X</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>=</mo>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>-</mo>
<mn>1</mn>
<mo>...</mo>
<mn>..</mn>
<mi>X</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo><</mo>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<mi>p</mi>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
<mo>...</mo>
<mo>...</mo>
<mo>|</mo>
<mi>X</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>|</mo>
<mo><</mo>
<mn>1</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
<mo>...</mo>
<mo>...</mo>
<mo>|</mo>
<mi>X</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>|</mo>
<mo>></mo>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<mi>s</mi>
<mo>=</mo>
<mfrac>
<mi>M</mi>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>p</mi>
<mo>+</mo>
<mi>I</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
X(t)表示相关后道数据，sign和p表示条件系数，Xlog(t)表示对数变换后的相关后道数据，M表示允许输出的最大值，I表示X(t)中的最大值。
5.一种滑动扫描谐波压制装置，其特征在于，包括：
获取模块，用于获取滑动扫描相关后道数据，并拾取所述相关后道数据的初至时间；
CEEMD分解模块，用于将所述相关后道数据进行CEEMD分解，得到一组IMF分量；
谐波分析模块，用于分析谐波在所述IMF分量中的分布情况，确定含有谐波的IMF分量和不含有谐波的IMF分量；
谐波衰减模块，用于将所述含有谐波的IMF分量从所述初至时间开始进行谐波衰减处理；
重构模块，用于利用谐波衰减处理后的IMF分量以及所述不含有谐波的IMF分量进行重构得到去谐波后的相关后道数据；
标准化预处理模块，用于在将所述相关后道数据进行CEEMD分解之前，将所述相关后道数据进行标准化预处理；
对应的，所述将相关后道数据进行CEEMD分解中，所述相关后道数据为标准化预处理后的相关后道数据；其中，
所述标准化预处理模块，具体包括：
对数变换子模块，用于将所述相关后道数据进行对数变换；
标准差子模块，用于计算对数变换后的相关后道数据的标准差；
标准化子模块，用于将所述对数变换后的相关后道数据除以所述对数变换后的相关后道数据的标准差获得标准化预处理后的相关后道数据。
6.根据权利要求5所述的滑动扫描谐波压制装置，其特征在于，所述将含有谐波的IMF分量从所述初至时间开始进行谐波衰减处理，具体包括：
将所述含有谐波的IMF分量从所述初至时间开始乘以布莱克曼窗函数。
7.根据权利要求5所述的滑动扫描谐波压制装置，其特征在于，所述利用谐波衰减处理后的IMF分量以及所述不含有谐波的IMF分量进行重构得到去谐波后的相关后道数据的计算公式为：
<mrow>
<msub>
<msup>
<mi>X</mi>
<mo>^</mo>
</msup>
<mrow>
<mi>d</mi>
<mi>e</mi>
<mi>n</mi>
<mi>o</mi>
<mi>i</mi>
<mi>s</mi>
<mi>e</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>2</mn>
</mrow>
<mn>4</mn>
</munderover>
<msubsup>
<mi>C</mi>
<mi>m</mi>
<mo>^</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>5</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>C</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
其中，C^m(t)为谐波衰减处理后的IMF分量；Ci(t)为不含有谐波的IMF分量；X^denoise(t)为去谐波后的相关后道数据，n＝log2(N)，N为样点值个数。
8.根据权利要求5所述的滑动扫描谐波压制装置，其特征在于，所述将相关后道数据进行对数变换的具体计算公式为：
Xlog(t)＝sign[slog(p+|X(t)|)]
其中：
<mrow>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
<mo>...</mo>
<mo>...</mo>
<mi>X</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>></mo>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
<mo>...</mo>
<mo>...</mo>
<mi>X</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>=</mo>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>-</mo>
<mn>1</mn>
<mo>...</mo>
<mn>..</mn>
<mi>X</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo><</mo>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<mi>p</mi>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
<mo>...</mo>
<mo>...</mo>
<mo>|</mo>
<mi>X</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>|</mo>
<mo><</mo>
<mn>1</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
<mo>...</mo>
<mo>...</mo>
<mo>|</mo>
<mi>X</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>|</mo>
<mo>></mo>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<mi>s</mi>
<mo>=</mo>
<mfrac>
<mi>M</mi>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>p</mi>
<mo>+</mo>
<mi>I</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
X(t)表示已拾取初至时间的相关后道数据，sign和p表示条件系数，Xlog(t)表示对数变换后的道数据，M表示允许输出的最大值，I表示X(t)中的最大值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510504580.2A CN105093328B (zh) | 2015-08-17 | 2015-08-17 | 一种滑动扫描谐波压制方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510504580.2A CN105093328B (zh) | 2015-08-17 | 2015-08-17 | 一种滑动扫描谐波压制方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105093328A CN105093328A (zh) | 2015-11-25 |
| CN105093328B true CN105093328B (zh) | 2018-02-02 |
Family
ID=54574182
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510504580.2A Active CN105093328B (zh) | 2015-08-17 | 2015-08-17 | 一种滑动扫描谐波压制方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105093328B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106842323A (zh) * | 2015-12-04 | 2017-06-13 | 中国石油化工股份有限公司 | 一种基于分频处理的滑动扫描谐波干扰压制方法 |
| CN105700020A (zh) * | 2016-03-23 | 2016-06-22 | 中国石油天然气集团公司 | 一种地震数据随机噪声压制方法及装置 |
| CN105911590B (zh) * | 2016-05-09 | 2019-08-09 | 刘宏杰 | 一种谐波噪声的压制处理方法及装置 |
| CN106556647B (zh) * | 2016-11-18 | 2019-02-22 | 湘潭大学 | 一种冲击回波数据处理方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6950371B2 (en) * | 2003-05-20 | 2005-09-27 | Chevron U.S.A. Inc. | Method for signal-to-noise ratio enhancement of seismic data using frequency dependent true relative amplitude noise attenuation |
| CN102073066B (zh) * | 2009-11-25 | 2012-07-18 | 中国石油天然气集团公司 | 一种消除地震数据谐波干扰的方法 |
| CN102262243B (zh) * | 2010-05-31 | 2013-04-24 | 中国石油天然气集团公司 | 一种滤波法可控震源地震数据谐波干扰压制方法 |
| CN102944895B (zh) * | 2012-10-26 | 2015-04-15 | 中国石油大学(北京) | 一种地震记录面波干扰压制方法及装置 |
-
2015
- 2015-08-17 CN CN201510504580.2A patent/CN105093328B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105093328A (zh) | 2015-11-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wang et al. | Residual learning of deep convolutional neural network for seismic random noise attenuation | |
| CN105093328B (zh) | 一种滑动扫描谐波压制方法及装置 | |
| CN104395779B (zh) | 用于估计和减弱地震数据中的噪声的系统和方法 | |
| CA2599958A1 (en) | Removal of noise from seismic data using radon transformations | |
| CN102945548A (zh) | 一种基于方向金字塔滤波的图像处理方法及装置 | |
| CN107179550B (zh) | 一种数据驱动的地震信号零相位反褶积方法 | |
| CN104849757A (zh) | 消除地震信号中随机噪声系统及方法 | |
| CN105319593A (zh) | 基于曲波变换和奇异值分解的联合去噪方法 | |
| CN104133248A (zh) | 一种高保真声波干扰压制方法 | |
| You et al. | Attenuation of seismic swell noise using convolutional neural networks in frequency domain and transfer learning | |
| Li et al. | Distributed acoustic sensing vertical seismic profile data denoising based on multistage denoising network | |
| Yuan et al. | Attenuation of linear noise based on denoising convolutional neural network with asymmetric convolution blocks | |
| Zhang et al. | Investigation of microseismic signal denoising using an improved wavelet adaptive thresholding method | |
| Li et al. | Removing abnormal environmental noise in nodal land seismic data using deep learning | |
| CN117031539A (zh) | 一种自监督深度学习地震数据低频重建方法及系统 | |
| Cheng et al. | Multiscale recurrent-guided denoising network for distributed acoustic sensing-vertical seismic profile background noise attenuation | |
| CN115797318A (zh) | 一种光谱数据预处理方法、装置、计算机设备及存储介质 | |
| Lin et al. | Removal of a series of spikes from magnetic resonance sounding signal by combining empirical mode decomposition and wavelet thresholding | |
| CN109212609A (zh) | 基于波动方程延拓的近地表噪音压制方法 | |
| CN103954994A (zh) | 基于连续小波变换的地震信号增强方法及装置 | |
| CN116009076A (zh) | 地震信号处理方法、介质及设备 | |
| Sun et al. | Application of adaptive iterative low-rank algorithm based on transform domain in desert seismic signal analysis | |
| CN104422961A (zh) | 一种地震随机噪声衰减方法 | |
| Oboué et al. | An advanced median filter for improving the signal-to-noise ratio of seismological datasets | |
| CN113158797B (zh) | 一种微震资料去噪方法、系统、电子设备及存储介质 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |