CN109655893A - 一种波形自适应匹配的可控震源谐波消除方法及系统 - Google Patents
一种波形自适应匹配的可控震源谐波消除方法及系统 Download PDFInfo
- Publication number
- CN109655893A CN109655893A CN201710947734.4A CN201710947734A CN109655893A CN 109655893 A CN109655893 A CN 109655893A CN 201710947734 A CN201710947734 A CN 201710947734A CN 109655893 A CN109655893 A CN 109655893A
- Authority
- CN
- China
- Prior art keywords
- signal
- force signal
- amplitude
- reference scan
- controlled source
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 230000003044 adaptive effect Effects 0.000 title claims abstract description 29
- 230000008030 elimination Effects 0.000 title claims abstract description 14
- 238000003379 elimination reaction Methods 0.000 title claims abstract description 14
- 238000001914 filtration Methods 0.000 claims abstract description 49
- 238000001228 spectrum Methods 0.000 claims description 36
- 230000009466 transformation Effects 0.000 claims description 12
- 239000000284 extract Substances 0.000 claims description 11
- 230000004044 response Effects 0.000 claims description 10
- 230000001133 acceleration Effects 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 5
- 230000000717 retained effect Effects 0.000 abstract description 5
- 230000006870 function Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 241000252206 Cypriniformes Species 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提出了一种波形自适应匹配的可控震源谐波消除方法及系统，该方法包括：提取可控震源特征信号中的参考扫描信号和出力信号；将参考扫描信号作为波形匹配道，采用自适应匹配方法，将出力信号与参考扫描信号道进行匹配滤波；求取出力信号滤波前后的振幅变化系数；将振幅变化系数作用于滤波后的出力信号。本发明将震源特征信号中的参考扫描信号作为匹配滤波道，设计自适应滤波算子，作用于产生谐波畸变的震源力信号，消除力信号中的谐波干扰。求取原始力信号与滤波后的力信号振幅变化系数，将该系数对滤波后的振幅进行处理，有效保留力信号中有效信号振幅信息，最终达到消除可控震源力信号中谐波干扰的目的。
Description
技术领域
本发明属于地震资料处理领域，通过提取可控震源出力信号中有效信号、消除谐波干扰信号，提供一种波形自适应匹配的可控震源谐波消除方法及系统。
背景技术
在可控震源勘探过程中，由于震源机械装置和震动装置的非线性，以及震板与大地的耦合问题，导致可控震源输出的信号产生畸变。震源平板向地下传输扫描信号的同时，也向地下传输参考信号整数倍的谐波干扰。根据其频率的倍数，将这种谐波干扰依次分为二次谐波、三次谐波及n次谐波。
国内外学者提出了许多消除地震数据中高阶谐波，提取有效波信号的方法。Li等(1995)针对线性扫描方式提出了纯相移滤波方法，可有效地压制近场和远场的谐波畸变，它是针对传统的单炮点采集提出的，不适用于可控震源滑动扫描采集方式。针对滑动扫描数据，Fleure(2002)提出了一种在频率域分离谐波并对其压制的技术，Meunier和Bianchi等(2002)提出了一种通过预测谐波，再在每道记录中依次减去谐波的方法。Moerig(2007)提出通过对相关后数据中负时间域的谐波噪声进行统计，来估计地震记录中谐波能量，进而消除谐波畸变的方法。Abd El-Aal(2010)和Claudio Bagaini(2010)分别在时间域和频率域推导谐波预测算子来预测相关后地震记录中的谐波噪声，但它需要用到相关前数据，无法直接处理相关后数据。Charles Sicking(2009)利用地面力信号和扫描起始时间设计滤波器，将滤波器与基波能量褶积来估计该炮对上一炮记录产生的谐波噪声，该方法可直接用来处理相关后记录。
因此，本领域急需一种在有效保留原始地震道振幅变化信息的同时，消除震源信号中的谐波干扰信号的方法。
发明内容
可控震源激发过程中，震源平板向地下激发有效信号的同时，有效信号产生畸变生成整数倍频率的谐波干扰。为了解决现有技术中的问题，本发明提供一种自适应的波形匹配滤波方法，并对滤波后的振幅进行处理，在有效保留原始地震道振幅变化信息的同时，消除震源信号中的谐波干扰信号。
根据本发明的一个方面，提供一种波形自适应匹配的可控震源谐波消除方法，该方法包括：
提取可控震源特征信号中的参考扫描信号和出力信号；
将参考扫描信号作为波形匹配道，采用自适应匹配方法，将出力信号与参考扫描信号道进行匹配滤波；
求取出力信号滤波前后的振幅变化系数；
将振幅变化系数作用于滤波后的出力信号。
进一步地，提取可控震源特征信号包括：
提取可控震源辅助道中参考扫描信号、出力信号、平板加速度信号、重锤加速度信号；
分别将参考扫描信号和出力信号的两道地震信号数据提取出来。
进一步地，将出力信号与参考扫描信号进行匹配滤波包括：
设计自适应滤波器；
将自适应滤波器作用于出力信号，得到滤波后的出力信号道。
进一步地，设计自适应滤波器包括：
采用最小均方LMS算法作为自适应滤波器；
对自适应滤波器进行初始化，给定n时刻的抽头输入向量；
计算自适应滤波器输出对输入信号的响应；
比较输出结果与期望响应产生的估计误差；
根据估计误差自动调整滤波器参数。
进一步地，求取出力信号滤波后的振幅变化系数包括以下步骤：
对参考扫描信号作傅里叶变换，获得参考扫描信号的振幅谱；
对滤波前出力信号作傅里叶变换，获得出力信号的振幅谱；
计算每一个频率点处参考扫描信号振幅值与滤波前出力信号振幅值的比值，将该比值作为出力信号滤波后在此频率点处的振幅变化系数。
进一步地，将振幅变化系数作用于滤波后的出力信号道包括：
对滤波后的出力信号作傅里叶变换，获得滤波后的出力信号的振幅谱和相位谱；
将振幅变化系数和滤波后出力信号振幅谱作乘积，作为新的振幅谱，相位谱保持不变；
对新的振幅谱和相位谱作傅里叶逆变换，得到振幅变化作用后的出力信号。
根据本发明的另一方面，提供一种波形自适应匹配的可控震源谐波消除系统，该系统包括：
存储器，存储有计算机可执行指令；
处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：
提取可控震源特征信号中的参考扫描信号和出力信号；
将参考扫描信号作为波形匹配道，采用自适应匹配方法，将出力信号与参考扫描信号道进行匹配滤波；
求取出力信号滤波前后的振幅变化系数；
将振幅变化系数作用于滤波后的出力信号。
本发明利用震源特征信号中的参考扫描信号，消除出力信号中谐波干扰，从而获得激发的实际震源有效信号。
附图说明
通过结合附图对本公开示例性实施方式进行更详细的描述，本公开的上述以及其它目的、特征和优势将变得更加明显，其中，在本公开示例性实施方式中，相同的参考标号通常代表相同部件。
图1示出了本发明实施例的参考扫描信号的示意图。
图2示出了本发明实施例的含谐波干扰的出力信号。
图3示出了本发明实施例的去除谐波的出力信号和参考扫描信号对比图。
图4示出了本发明实施例的滤波后的出力信号与参考扫描信号的误差值。
图5示出了本发明实施例的实际出力信号。
图6示出了本发明实施例的出力信号的时频分析图。
图7示出了本发明实施例的滤波后的出力信号。
图8示出了本发明实施例的滤波后出力信号的时频图。
图9示出了本发明方法的流程图。
具体实施方式
下面将参照附图更详细地描述本公开的优选实施方式。虽然附图中显示了本公开的优选实施方式，然而应该理解，可以以各种形式实现本公开而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了使本公开更加透彻和完整，并且能够将本公开的范围完整地传达给本领域的技术人员。
本发明提供了一种基于波形匹配的滤波方法，将震源特征信号中的参考扫描信号作为匹配滤波道，设计自适应滤波算子，作用于产生谐波畸变的震源出力信号，消除出力信号中的谐波干扰。求取原始出力信号与滤波后的出力信号振幅变化系数，将该系数对滤波后的振幅进行处理，有效保留出力信号中有效信号振幅信息，最终达到消除可控震源出力信号中谐波干扰的目的。
根据本发明的一个实施方式，如图9所示，提供一种波形自适应匹配的可控震源谐波消除方法，该方法包括：
提取可控震源特征信号中的参考扫描信号和出力信号；
将参考扫描信号作为波形匹配道，采用自适应匹配方法，将出力信号与参考扫描信号道进行匹配滤波；
求取出力信号滤波前后的振幅变化系数；
将振幅变化系数作用于滤波后的出力信号。
优选地，提取可控震源特征信号包括：
提取可控震源辅助道中参考扫描信号、出力信号、平板加速度信号、重锤加速度信号；
分别将参考扫描信号和出力信号的两道地震信号数据提取出来。
优选地，将出力信号与参考扫描信号进行匹配滤波包括：
设计自适应滤波器；
具体地，采用最小均方LMS算法作为自适应滤波器，即将均方误差替换为瞬时误差的平方，以瞬时误差信号平方的梯度作为均方误差函数梯度的估计；
对自适应滤波器进行初始化，给定n时刻的抽头输入向量；
计算自适应滤波器输出对输入信号的响应；
比较输出结果与期望响应产生的估计误差；
根据估计误差自动调整滤波器参数；
最后将自适应滤波器作用于出力信号，得到滤波后的出力信号道。
优选地，求取出力信号滤波后的振幅变化系数包括以下步骤：
对参考扫描信号作傅里叶变换，获得参考扫描信号的振幅谱；
对滤波前出力信号作傅里叶变换，获得出力信号的振幅谱；
计算每一个频率点处参考扫描信号振幅值与滤波前出力信号振幅值的比值，将该比值作为出力信号滤波后在此频率点处的振幅变化系数。
最后，将振幅变化系数作用于滤波后的出力信号道包括：
对滤波后的出力信号作傅里叶变换，获得滤波后的出力信号的振幅谱和相位谱；
将振幅变化系数和滤波后出力信号振幅谱作乘积，作为新的振幅谱，相位谱保持不变；
对新的振幅谱和相位谱作傅里叶逆变换，得到振幅变化作用后的出力信号。
根据本发明的另一实施方式，提供一种波形自适应匹配的可控震源谐波消除系统，该系统包括：
存储器，存储有计算机可执行指令；
处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：
提取可控震源特征信号中的参考扫描信号和出力信号；
将参考扫描信号作为波形匹配道，采用自适应匹配方法，将出力信号与参考扫描信号道进行匹配滤波；
求取出力信号滤波前后的振幅变化系数；
将振幅变化系数作用于滤波后的出力信号。
本发明提供了一种基于波形匹配的滤波方法，将震源特征信号中的参考扫描信号作为匹配滤波道，设计自适应滤波算子，作用于产生谐波畸变的震源力信号，消除力信号中的谐波干扰。求取原始力信号与滤波后的力信号振幅变化系数，将该系数对滤波后的振幅进行处理，有效保留力信号中有效信号振幅信息。
为便于理解本发明实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
为了验证这种基于自适应波形匹配的可控震源谐波消除方法的方法可行性，对模拟和实际资料进行测试。设计频率范围为4～84hz线性升频扫描信号，采样率为0.5ms，扫描时间2s，如图1所示。
参考扫描信号发生畸变，生成频率整数倍的谐波信号，二阶谐波能量为基波信号的0.8倍，三阶谐波能量为基波信号的0.5倍。因此，含谐波干扰的出力信号，其中基波信号能量与参考扫描信号能量一致，如图2所示。
采用最小均方(LMS)算法来进行自适应滤波，即将均方误差替换为瞬时误差的平方，以瞬时误差信号平方的梯度作为均方误差函数梯度的估计。
首先，初始化线性自适应滤波器；
给定n时刻的抽头输入向量：
计算线性滤波器输出对输入信号的响应；
比较输出结果与期望响应产生的估计误差；
根据估计误差自动调整滤波器参数。
采用上述迭代过程，对图2中的谐波干扰进行滤波，得到不含谐波的出力信号如图3所示，图4为去除谐波后的出力信号与参考扫描信号误差。可以看到，除信号初始部分由于未达到完全收敛误差值较大外，随着迭代收敛误差值迅速减小，控制在2％以内。
当出力信号与参考扫描信号之间存在能量差异时，需要求取出力信号滤波前后振幅变化系数，并将该振幅变化系数作用于滤波后出力信号。
对某工区实际出力信号进行自适应滤波消除高阶谐波信号(如图5所示)。从出力信号的时频分析图中可以看到(如图6所示)，信号中存在能量较强的高阶谐波干扰。
采用滤波器阶数为200，步长因子为2.0*10-11的自适应滤波器，以参考扫描信号为目标对出力信号进行滤波处理，图7为滤除谐波干扰后的信号，图8为滤波后信号与参考扫描信号误差，图9为滤波后信号的时频图。可以看到谐波干扰很好得到消除。
以上已经描述了本公开的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施例的原理、实际应用或对市场中的技术改进，或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。
Claims (10)
1.一种波形自适应匹配的可控震源谐波消除方法，其特征在于，该方法包括：
提取可控震源特征信号中的参考扫描信号和出力信号；
将参考扫描信号作为波形匹配道，采用自适应匹配方法，将出力信号与参考扫描信号道进行匹配滤波；
求取出力信号滤波后的振幅变化系数；
将振幅变化系数作用于滤波后的出力信号。
2.根据权利要求1所述的波形自适应匹配的可控震源谐波消除方法，其特征在于，提取可控震源特征信号包括：
提取可控震源辅助道中参考扫描信号、出力信号、平板加速度信号、重锤加速度信号；
分别将参考扫描信号和出力信号的两道地震信号数据提取出来。
3.根据权利要求1所述的波形自适应匹配的可控震源谐波消除方法，其特征在于，将出力信号与参考扫描信号进行匹配滤波包括：
设计自适应滤波器；
将自适应滤波器作用于出力信号，得到滤波后的出力信号道。
4.根据权利要求3所述的波形自适应匹配的可控震源谐波消除方法，其特征在于，设计自适应滤波器包括：
采用最小均方LMS算法作为自适应滤波器；
对自适应滤波器进行初始化，给定n时刻的抽头输入向量；
计算自适应滤波器输出对输入信号的响应；
比较输出结果与期望响应产生的估计误差；
根据估计误差自动调整滤波器参数。
5.根据权利要求1所述的波形自适应匹配的可控震源谐波消除方法，其特征在于，求取出力信号滤波后的振幅变化系数包括以下步骤：
对参考扫描信号作傅里叶变换，获得参考扫描信号的振幅谱；
对滤波前出力信号作傅里叶变换，获得出力信号的振幅谱；
计算每一个频率点处参考扫描信号振幅值与滤波前出力信号振幅值的比值，将该比值作为出力信号滤波后在此频率点处的振幅变化系数。
6.根据权利要求1所述的波形自适应匹配的可控震源谐波消除方法，其特征在于，将振幅变化系数作用于滤波后的出力信号道包括：
对滤波后的出力信号作傅里叶变换，获得滤波后的出力信号的振幅谱和相位谱；
将振幅变化系数和滤波后出力信号振幅谱作乘积，作为新的振幅谱，相位谱保持不变；
对新的振幅谱和相位谱作傅里叶逆变换，得到振幅变化作用后的出力信号。
7.一种波形自适应匹配的可控震源谐波消除系统，其特征在于，该系统包括：
存储器，存储有计算机可执行指令；
处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：
提取可控震源特征信号中的参考扫描信号和出力信号；
将参考扫描信号作为波形匹配道，采用自适应匹配方法，将出力信号与参考扫描信号道进行匹配滤波；
求取出力信号滤波前后的振幅变化系数；
将振幅变化系数作用于滤波后的出力信号。
8.根据权利要求7所述的波形自适应匹配的可控震源谐波消除系统，其特征在于，将出力信号与参考扫描信号进行匹配滤波包括：
设计自适应滤波器；
将自适应滤波器作用于出力信号，得到滤波后的出力信号道。
其中，设计自适应滤波器包括：
采用最小均方LMS算法作为自适应滤波器；
对自适应滤波器进行初始化，给定n时刻的抽头输入向量；
计算自适应滤波器输出对输入信号的响应；
比较输出结果与期望响应产生的估计误差；
根据估计误差自动调整滤波器参数。
9.根据权利要求7所述的波形自适应匹配的可控震源谐波消除系统，其特征在于，求取出力信号滤波后的振幅变化系数包括以下步骤：
对参考扫描信号作傅里叶变换，获得参考扫描信号的振幅谱；
对滤波前出力信号作傅里叶变换，获得出力信号的振幅谱；
计算每一个频率点处参考扫描信号振幅值与滤波前出力信号振幅值的比值，将该比值作为出力信号滤波后在此频率点处的振幅变化系数。
10.根据权利要求9所述的波形自适应匹配的可控震源谐波消除系统，其特征在于，将振幅变化系数作用于滤波后的出力信号道包括：
对滤波后的出力信号作傅里叶变换，获得滤波后的出力信号的振幅谱和相位谱；
将振幅变化系数和滤波后出力信号振幅谱作乘积，作为新的振幅谱，相位谱保持不变；
对新的振幅谱和相位谱作傅里叶逆变换，得到振幅变化作用后的出力信号。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710947734.4A CN109655893B (zh) | 2017-10-12 | 2017-10-12 | 一种波形自适应匹配的可控震源谐波消除方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710947734.4A CN109655893B (zh) | 2017-10-12 | 2017-10-12 | 一种波形自适应匹配的可控震源谐波消除方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109655893A true CN109655893A (zh) | 2019-04-19 |
| CN109655893B CN109655893B (zh) | 2020-06-19 |
Family
ID=66109451
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710947734.4A Active CN109655893B (zh) | 2017-10-12 | 2017-10-12 | 一种波形自适应匹配的可控震源谐波消除方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109655893B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112666597A (zh) * | 2019-10-16 | 2021-04-16 | 中国石油天然气集团有限公司 | 可控震源扫描方法和装置 |
| CN113050166A (zh) * | 2019-12-27 | 2021-06-29 | 中国石油天然气集团有限公司 | 地震激发信号修正方法及装置 |
| CN113127799A (zh) * | 2020-01-15 | 2021-07-16 | 中国石油天然气集团有限公司 | 一种扫描信号校正方法及装置 |
| CN113138417A (zh) * | 2020-01-20 | 2021-07-20 | 中国石油天然气集团有限公司 | 一种基于力信号反馈的扫描信号自动校正方法及装置 |
| CN113608018A (zh) * | 2021-06-30 | 2021-11-05 | 中冶南方都市环保工程技术股份有限公司 | 提高谐波检测精度的自适应vmd检测方法、装置及存储介质 |
| CN116068613A (zh) * | 2021-10-29 | 2023-05-05 | 中国石油化工股份有限公司 | 一种基于谐波倍频特征的可控震源高效扫描方法及系统 |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100085837A1 (en) * | 2008-10-03 | 2010-04-08 | Kambiz Iranpour | Harmonic attenuation using multiple sweep rates |
| US20100142320A1 (en) * | 2007-03-09 | 2010-06-10 | Claudio Bagaini | Vibroseis acquisition method |
| US20110273966A1 (en) * | 2010-05-05 | 2011-11-10 | Emmanuel Coste | Reducing Noise in Seismic Surveying Environment |
| WO2011144215A2 (fr) * | 2010-05-20 | 2011-11-24 | Entreprise Nationale De Geophysique - Enageo- Filiale Du Groupe Sonatrach | Méthode d'atténuation du bruit harmonique en vibrosismique par filtrage temps-variant avec référence |
| CN102262243A (zh) * | 2010-05-31 | 2011-11-30 | 中国石油天然气集团公司 | 一种滤波法可控震源地震数据谐波干扰压制方法 |
| CN102478671A (zh) * | 2010-11-23 | 2012-05-30 | 中国石油天然气集团公司 | 一种压制可控震源谐波干扰的方法 |
| CN102692649A (zh) * | 2011-03-23 | 2012-09-26 | 中国石油天然气集团公司 | 一种可控震源地震勘探中消除谐波干扰的方法 |
| CN104216010A (zh) * | 2013-06-05 | 2014-12-17 | 上海青凤致远地球物理地质勘探科技有限公司 | 利用可控震源谐波提高地震数据质量的方法 |
| US20150316672A1 (en) * | 2014-05-01 | 2015-11-05 | Conocophillips Company | Deterministic phase correction and application |
| CN105277983A (zh) * | 2014-06-26 | 2016-01-27 | 中石化石油工程地球物理有限公司胜利分公司 | 新型压制可控震源滑动扫描谐波干扰的方法 |
| CN105277986A (zh) * | 2015-05-27 | 2016-01-27 | 中石化石油工程地球物理有限公司胜利分公司 | 基于自适应匹配滤波算子的可控震源谐波压制方法 |
| CN106094024A (zh) * | 2016-06-21 | 2016-11-09 | 中国石油天然气集团公司 | 一种可控震源滑动扫描数据采集处理方法及装置 |
-
2017
- 2017-10-12 CN CN201710947734.4A patent/CN109655893B/zh active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100142320A1 (en) * | 2007-03-09 | 2010-06-10 | Claudio Bagaini | Vibroseis acquisition method |
| US20100085837A1 (en) * | 2008-10-03 | 2010-04-08 | Kambiz Iranpour | Harmonic attenuation using multiple sweep rates |
| US20110273966A1 (en) * | 2010-05-05 | 2011-11-10 | Emmanuel Coste | Reducing Noise in Seismic Surveying Environment |
| WO2011144215A2 (fr) * | 2010-05-20 | 2011-11-24 | Entreprise Nationale De Geophysique - Enageo- Filiale Du Groupe Sonatrach | Méthode d'atténuation du bruit harmonique en vibrosismique par filtrage temps-variant avec référence |
| CN102262243A (zh) * | 2010-05-31 | 2011-11-30 | 中国石油天然气集团公司 | 一种滤波法可控震源地震数据谐波干扰压制方法 |
| CN102478671A (zh) * | 2010-11-23 | 2012-05-30 | 中国石油天然气集团公司 | 一种压制可控震源谐波干扰的方法 |
| CN102692649A (zh) * | 2011-03-23 | 2012-09-26 | 中国石油天然气集团公司 | 一种可控震源地震勘探中消除谐波干扰的方法 |
| CN104216010A (zh) * | 2013-06-05 | 2014-12-17 | 上海青凤致远地球物理地质勘探科技有限公司 | 利用可控震源谐波提高地震数据质量的方法 |
| US20150316672A1 (en) * | 2014-05-01 | 2015-11-05 | Conocophillips Company | Deterministic phase correction and application |
| CN105277983A (zh) * | 2014-06-26 | 2016-01-27 | 中石化石油工程地球物理有限公司胜利分公司 | 新型压制可控震源滑动扫描谐波干扰的方法 |
| CN105277986A (zh) * | 2015-05-27 | 2016-01-27 | 中石化石油工程地球物理有限公司胜利分公司 | 基于自适应匹配滤波算子的可控震源谐波压制方法 |
| CN106094024A (zh) * | 2016-06-21 | 2016-11-09 | 中国石油天然气集团公司 | 一种可控震源滑动扫描数据采集处理方法及装置 |
Non-Patent Citations (4)
| Title |
|---|
| YAO JIANJUN，等: "Harmonic cancellation for electro-hydraulic servo shaking table based on LMS adaptive algorithm", 《JOURNAL OF VIBRATION AND CONTROL》 * |
| 冯玉苹，等: "可控震源白适应谐波压制技术集成及应用", 《2017年物探技术研讨会》 * |
| 曲英铭，等: "自适应匹配预测滤波压制可控震源谐波", 《石油地球物理勘探》 * |
| 李璐璐，等: "自适应匹配预测滤波可控震源谐波压制方法研究", 《中国地球科学联合学术年会 2015》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112666597A (zh) * | 2019-10-16 | 2021-04-16 | 中国石油天然气集团有限公司 | 可控震源扫描方法和装置 |
| CN112666597B (zh) * | 2019-10-16 | 2024-03-26 | 中国石油天然气集团有限公司 | 可控震源扫描方法和装置 |
| CN113050166A (zh) * | 2019-12-27 | 2021-06-29 | 中国石油天然气集团有限公司 | 地震激发信号修正方法及装置 |
| CN113127799A (zh) * | 2020-01-15 | 2021-07-16 | 中国石油天然气集团有限公司 | 一种扫描信号校正方法及装置 |
| CN113127799B (zh) * | 2020-01-15 | 2023-12-01 | 中国石油天然气集团有限公司 | 一种扫描信号校正方法及装置 |
| CN113138417A (zh) * | 2020-01-20 | 2021-07-20 | 中国石油天然气集团有限公司 | 一种基于力信号反馈的扫描信号自动校正方法及装置 |
| CN113608018A (zh) * | 2021-06-30 | 2021-11-05 | 中冶南方都市环保工程技术股份有限公司 | 提高谐波检测精度的自适应vmd检测方法、装置及存储介质 |
| CN116068613A (zh) * | 2021-10-29 | 2023-05-05 | 中国石油化工股份有限公司 | 一种基于谐波倍频特征的可控震源高效扫描方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109655893B (zh) | 2020-06-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109655893A (zh) | 一种波形自适应匹配的可控震源谐波消除方法及系统 | |
| Xiao et al. | Adaptive MOMEDA based on improved advance-retreat algorithm for fault features extraction of axial piston pump | |
| CN109884192B (zh) | 用于焊缝特征导波缺陷回波信号特性提取的稀疏表示方法 | |
| CN101627322A (zh) | 可控震源采集方法 | |
| CN102692649A (zh) | 一种可控震源地震勘探中消除谐波干扰的方法 | |
| KR20140021584A (ko) | 스펙트럼 형상을 이용하는 전 파동장 역산의 수렴 레이트 | |
| Qu et al. | An effective method for reducing harmonic distortion in correlated vibroseis data | |
| CN105786038B (zh) | 一种时频域混合自适应滤波的结构振动主动控制方法 | |
| CN102116868A (zh) | 一种地震波分解方法 | |
| CN101545981A (zh) | 可控震源地震数据零相位子波最小相位化方法 | |
| CN106680874A (zh) | 基于波形形态特征稀疏化建模的谐波噪声压制方法 | |
| CN102419972B (zh) | 一种声音信号检测和识别的方法 | |
| CN104614767A (zh) | 基于分段延拓的时变地震子波相位校正方法 | |
| CN105954803B (zh) | 叠后地震反演方法及装置 | |
| CN105277983B (zh) | 压制可控震源滑动扫描谐波干扰的方法 | |
| CN114355446A (zh) | 一种基于高分辨时频变换的可控震源地震资料黑三角区噪音压制方法 | |
| CN103728661B (zh) | 一种高精度反q滤波地震资料处理方法 | |
| CN102735759A (zh) | 一种基于脊的兰姆波信号去噪方法 | |
| CN117849192A (zh) | 一种基于广义频散模态分解方法的超声导波模式分离方法 | |
| CN107390261B (zh) | 基于线性Bregman算法的表面多次波和子波估计方法及系统 | |
| CN120948633A (zh) | 基于非线性最小二乘法引导的s-g滤波超声信号降噪方法 | |
| CN111650655B (zh) | 一种非负矩阵分解的有监督瞬变电磁信号降噪方法 | |
| CN116047589B (zh) | 基于牛顿迭代的表层多次波压制方法、系统及识别仪 | |
| CN114428346A (zh) | 一种谐波压制方法及系统 | |
| Hu et al. | Multistep full-waveform inversion based on waveform-mode decomposition |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |