CN107255835B - 一种用于压制机械因素谐波干扰的方法 - Google Patents
一种用于压制机械因素谐波干扰的方法 Download PDFInfo
- Publication number
- CN107255835B CN107255835B CN201710448762.1A CN201710448762A CN107255835B CN 107255835 B CN107255835 B CN 107255835B CN 201710448762 A CN201710448762 A CN 201710448762A CN 107255835 B CN107255835 B CN 107255835B
- Authority
- CN
- China
- Prior art keywords
- big gun
- harmonic wave
- wave interference
- harmonic
- ground force
- 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
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. analysis, for interpretation, for correction
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
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)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明涉及地震资料保幅去噪技术领域,特别是涉及一种用于压制机械因素谐波干扰的方法,该压制方法可靠性高、易于实现,可以有效的压制机械因素产生的谐波干扰,且不会对与谐波干扰频率相近的有效信号产生损害。该压制方法包括如下步骤:获取可控震源原始单炮记录、单炮地面力信号以及声波测井曲线;计算第i炮中各次谐波分量、第i炮的总谐波、地下反射系数;计算谐波畸变信号产生的谐波干扰;计算得到第i炮产生的总谐波干扰;确定并消除第i炮对其它临近炮的谐波干扰。
Description
技术领域
本发明涉及地震资料保幅去噪技术领域,特别是涉及一种用于压制机械因素谐波干扰的方法。
背景技术
在可控震源勘探过程中,由于种种因素震源车向地下传输扫描信号时不可避免的会产生谐波干扰。其中产生谐波干扰最重要的原因之一就是机械因素,由于可控震源的机械装置、震动装置以及液压伺服系统的非线性,导致从可控震源输出的信号就存在谐波干扰。这种谐波干扰与地表结构无关,也就是说,无论在什么样的地表结构下作业,都会产生这种谐波干扰。这种谐波信号可以通过可控震源震板上的输出信号记录下来,该谐波干扰存在两个明显特征:①谐波分量的频率为有效波频率的整数倍。如二次谐波的频率范围为有效波频率范围两倍,且它们的频率随有效波频率的变化而变化;②谐波干扰主要有本炮干扰和临炮干扰两种类型,线性升频信号本炮干扰出现在负时间轴,对临炮干扰,采用滑动扫描采集时,后一炮产生在负时间轴上的谐波与前一炮混在一起,对前一炮会造成污染。
为衰减或消除上述机械因素产生的谐波干扰,现有技术通常采用压制机械系统谐波的方法例如:纯相移滤波法、分数Fourier方法等。
纯相移滤波法处理比较简单,容易操作,只需要利用K次谐波的相位谱设计一个纯相移滤波器,就可以进行滤波处理。但仅适用于当前炮的谐波压制,不适合高效采集资料,并且很容易损伤有效信号和谐波频率相近的成分,造成资料失真,保真度不好。分数Fourier是一种特殊的时频分析工具,它在对非平稳信号的分析和处理中具有许多优良的特性。可以对可控震源记录信号进行分数Fourier变换,对可控震源扫描信号在某个分数阶上实现能量聚焦,表现为一个脉冲函数,而噪声的能量却分布在整个时频域,在任何分数阶上都不会实现能量聚焦。这样我们就可以利用分数Fourier的这一特性,压制噪音,该方法尽管可以一定程度上压制谐波干扰,但针对性不强。上述这些方法处理比较简单,但容易损伤有效信号和谐波相近的频率成分。因此设计一种专门针对机械因素谐波的衰减方法是非常重要的。
发明内容
本发明提供了一种用于压制机械因素谐波干扰的方法,该压制方法可靠性高、易于实现,可以有效的压制机械因素产生的谐波干扰,且不会对与谐波干扰频率相近的有效信号产生损害。
为解决上述技术问题,本发明采用了如下技术方案:
一种用于压制机械因素谐波干扰的方法,所述方法包括如下步骤:
步骤1:获取可控震源原始单炮记录、单炮地面力信号以及声波测井曲线;
步骤2:根据第i炮的地面力信号,计算第i炮中各次谐波分量,并计算得到第i炮的总谐波,其中i为大于等于2的正整数;
根据第i炮的地面力信号以及从声波测井曲线中提取的滑动扫描第i炮振动记录,计算得到地下反射系数;
步骤3:将步骤2所得第i炮的总谐波与从声波测井曲线中提取的滑动扫描第i炮振动记录做互相关处理,得到谐波畸变信号产生的谐波干扰HARMi;
步骤4:根据步骤2所得地下反射系数以及步骤3所得谐波畸变信号产生的谐波干扰HARMi,计算得到第i炮产生的总谐波干扰;
步骤5:从声波测井曲线中提取的滑动扫描第i炮振动记录中确定出扫描信号长度T、记录长度tl以及滑动时间th,确定并消除第i炮对其它临近炮的谐波干扰。
进一步的,所述方法还包括:
步骤6:重复循环步骤2-步骤5,分别逐项确定并消除第i-1炮、第i-2炮、…、第1炮对其它临近炮的谐波干扰。
优选的,所述步骤2中根据第i炮的地面力信号,计算第i炮中各次谐波分量的具体过程可描述为:
将第i炮的地面力信号转换为频率域表示,得到频率域第i炮的地面力信号gi(ω),其中,频率域第i炮的地面力信号gi(ω)满足:
gi(ω)满足:
gi(ω)为频率域第i炮的地面力信号,Fi(ω)为频率域第i炮的基波信号,hik(ω)为频率域第i炮的k次谐波,i≧2,k≧2且均为正整数;
将得到的频率域第i炮的k次谐波与频率做互相关,即可得到第i炮的各次谐波分量。
本发明提供了一种用于压制机械因素谐波干扰的方法,该压制方法中包括有获取可控震源原始单炮记录、单炮地面力信号以及声波测井曲线;计算第i炮中各次谐波分量、第i炮的总谐波、地下反射系数;计算谐波畸变信号产生的谐波干扰HARMi;计算得到第i炮产生的总谐波干扰;确定并消除第i炮对其它临近炮的谐波干扰;进一步分别逐项确定并消除第i-1炮、第i-2炮、…、第1炮对其它临近炮的谐波干扰等步骤。具有上述步骤的本方法用于专门压制机械系统因素所产生的谐波干扰,其具体压制方式设置简单,操作易于实现运算速度快,适合应用在对三维可控震源地震资料的处理过程中。
附图说明
图1为本发明提供的一种用于压制机械因素谐波干扰的方法的流程示意图;
图2为畸变信号相关记录;
图3为声波时差曲线的示意图;
图4为消除第i炮对某一炮的谐波干扰前后的频率域显示效果;
图5为消除第i炮对某一炮的谐波干扰前后的单炮效果显示。
具体实施方式
本发明提供了一种用于压制机械因素谐波干扰的方法,该压制方法可靠性高、易于实现,可以有效的压制机械因素产生的谐波干扰,且不会对与谐波干扰频率相近的有效信号产生损害。
下面结合下述附图对本发明实施例做详细描述。
本发明提供了一种用于压制机械因素谐波干扰的方法,如图1所示,该压制方法具体包括有:
步骤1:获取可控震源原始单炮记录、单炮地面力信号以及声波测井曲线;
为验证本发明提供的压制方法确实有效,以XX油田XX区块的实际资料作为靶区进行了试验。在对靶区进行测井的过程中采用的是声波升频测井扫描方式,基波信号初始频率选取为10HZ,终止频率选取为72HZ,扫描时间为12s,采样间隔为0.002s,滑动时间为9.5s,接收时间为7s。测井过程中分别收集并获取得到可控震源原始单炮记录、单炮地面力信号以及声波测井曲线。
步骤2:根据第i炮的地面力信号,计算第i炮中各次谐波分量,并计算得到第i炮的总谐波,其中i为大于等于2的正整数;根据第i炮的地面力信号以及从声波测井曲线中提取的滑动扫描第i炮振动记录,计算得到地下反射系数;
进一步计算第i炮中各次谐波分量以及第i炮的总谐波。其中,i为大于等于2的正整数。值得注意的是,在本发明确定并消除谐波干扰的过程中,首先消除第i炮(优选将第i炮选定为最后一炮)产生的谐波干扰,而后顺次消除第i-1炮产生的谐波干扰、第i-2炮产生的谐波干扰、…、直至第一炮产生的谐波干扰。
具体计算第i炮中各次谐波分量的过程可描述为:首先利用频率域第i炮的地面力信号gi(ω)与频率域第i炮的基波信号Fi(ω)之间的关联,确定频率域第i炮的k次谐波hik(ω)。
上述三者满足:(i≧2,k≧2且均为正整数;)
而得到频率域第i炮的k次谐波hik(ω)之后,将频率域第i炮的k次谐波与频率做互相关,即可得到第i炮的各次谐波分量;进一步将第i炮的各次谐波分量累加即可得到第i炮的的总谐波。
而地下反射系数满足:Di'=giRi式(2);
其中:Ri为第i炮所对应的地下反射系数;gi为第i炮的地面力信号;Di'为从声波测井曲线中提取的滑动扫描第i炮振动记录。而进一步还可对式(2)做频率域表达处理,则可将式(2)变为:
Di'(ω)=gi(ω)Ri(ω)式(3);
其中,Ri(ω)表示第i炮与接收点中点位置处反射系数的频率域表示。
步骤3:将步骤2所得第i炮的总谐波与从声波测井曲线中提取的滑动扫描第i炮振动记录做互相关处理,得到谐波畸变信号产生的谐波干扰HARMi;
步骤4:根据步骤2所得地下反射系数以及步骤3所得谐波畸变信号产生的谐波干扰HARMi,计算得到第i炮产生的总谐波干扰;
在完成步骤2的基础上,进一步计算第i炮产生的总谐波干扰。其中第i炮产生的总谐波干扰hi满足:hi=HARMi*Ri式(4);
其中,谐波畸变信号产生的谐波干扰HARMi是由第i炮的总谐波与从声波测井曲线中提取的滑动扫描第i炮振动记录互相关处理得到的;Ri则通过步骤2的公式(2)即可得到。具体以测井扫描为例,产生的谐波畸变信号的相关记录可参考如图2。
值得注意的是,为提高反射系数的精确度,可进一步利用相邻层的密度以及声波时差曲线计算所得速度对反射系数进行修正。具体的,反射系数的修正满足:
其中,a代表某一地层,a-1代表a的下一地层,ρ为密度,v是相同位置处的声波时差曲线计算所得速度;声波时差曲线可由声波测井曲线中得出,具体如图3所示。
步骤5:从声波测井曲线中提取的滑动扫描第i炮振动记录中确定出扫描信号长度T、记录长度tl以及滑动时间th,确定并消除第i炮对其它临近炮的谐波干扰。
在完成步骤3、步骤4的基础上,进一步确定第i炮对其它临近炮的谐波干扰。首先,从声波测井曲线中提取的滑动扫描第i炮振动记录中确定出扫描信号长度T、记录长度ti以及滑动时间th。具体的,第i炮影响邻近的炮的影响长度为T+ti。在确定第i炮对其它炮的谐波干扰后,在其它炮振动记录中将第i炮产生的谐波干扰去掉即可。具体的可参照图4、图5,其中图4为消除第i炮对某一炮的谐波干扰前后的频率域显示效果;图5为消除第i炮对某一炮的谐波干扰前后的单炮显示效果。
如图1所示,进一步的,本发明压制方法还包括有:
步骤6:重复循环步骤2-步骤5,分别逐项确定并消除第i-1炮、第i-2炮、…、第1炮对其它临近炮的谐波干扰。
在完成步骤5的基础上,进一步循环步骤2-步骤5,依次确定并消除第i-1炮、第i-2炮、…,直至消除第1炮对其它临近炮的谐波干扰。也就是说,通过循环步骤2-步骤5可将所有炮对其它炮的谐波干扰(在第i炮为最后一炮的前提下)。
本发明提供了一种用于压制机械因素谐波干扰的方法,该压制方法中包括有获取可控震源原始单炮记录、单炮地面力信号以及声波测井曲线;计算第i炮中各次谐波分量、第i炮的总谐波、地下反射系数;计算谐波畸变信号产生的谐波干扰HARMi;计算得到第i炮产生的总谐波干扰;确定并消除第i炮对其它临近炮的谐波干扰;进一步分别逐项确定并消除第i-1炮、第i-2炮、…、第1炮对其它临近炮的谐波干扰等步骤。具有上述步骤的本方法用于专门压制机械系统因素所产生的谐波干扰,其具体压制方式设置简单,操作易于实现运算速度快,适合应用在对三维可控震源地震资料的处理过程中。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。
Claims (2)
1.一种用于压制机械因素谐波干扰的方法,其特征在于,所述方法包括如下步骤:
步骤1:获取可控震源原始单炮记录、单炮地面力信号以及声波测井曲线;
步骤2:根据第i炮的地面力信号,计算第i炮中各次谐波分量,并计算得到第i炮的总谐波,其中i为大于等于2的正整数;
所述根据第i炮的地面力信号,计算第i炮中各次谐波分量的具体过程可描述为:
将第i炮的地面力信号转换为频率域表示,得到频率域第i炮的地面力信号gi(ω),其中,频率域第i炮的地面力信号gi(ω)满足:
gi(ω)满足:
gi(ω)为频率域第i炮的地面力信号,Fi(ω)为频率域第i炮的基波信号,hik(ω)为频率域第i炮的k次谐波,i≧2,k≧2且均为正整数;其中,i为炮的次序,k为第i炮中的总谐波数,m为第i炮中谐波的次序;
将得到的频率域第i炮的k次谐波与频率做互相关,即可得到第i炮的各次谐波分量;
根据第i炮的地面力信号以及从声波测井曲线中提取的滑动扫描第i炮振动记录,计算得到地下反射系数;
步骤3:将步骤2所得第i炮的总谐波与从声波测井曲线中提取的滑动扫描第i炮振动记录做互相关处理,得到谐波畸变信号产生的谐波干扰HARMi;
步骤4:根据步骤2所得地下反射系数以及步骤3所得谐波畸变信号产生的谐波干扰HARMi,计算得到第i炮产生的总谐波干扰;
步骤5:从声波测井曲线中提取的滑动扫描第i炮振动记录中确定出扫描信号长度T、记录长度tl以及滑动时间th,确定并消除第i炮对其它临近炮的谐波干扰。
2.根据权利要求1所述的一种压制机械因素谐波干扰的方法,其特征在于,所述方法还包括:
步骤6:重复循环步骤2-步骤5,分别逐项确定并消除第i-1炮、第i-2炮、…、第1炮对其它临近炮的谐波干扰。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710448762.1A CN107255835B (zh) | 2017-06-14 | 2017-06-14 | 一种用于压制机械因素谐波干扰的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710448762.1A CN107255835B (zh) | 2017-06-14 | 2017-06-14 | 一种用于压制机械因素谐波干扰的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107255835A CN107255835A (zh) | 2017-10-17 |
CN107255835B true CN107255835B (zh) | 2019-01-01 |
Family
ID=60023721
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710448762.1A Active CN107255835B (zh) | 2017-06-14 | 2017-06-14 | 一种用于压制机械因素谐波干扰的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107255835B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6161076A (en) * | 1997-11-14 | 2000-12-12 | Baker Hughes Incorporated | Seismic data acquisition and processing using non-linear distortion in a vibratory output signal |
CN102692649A (zh) * | 2011-03-23 | 2012-09-26 | 中国石油天然气集团公司 | 一种可控震源地震勘探中消除谐波干扰的方法 |
CN102798894A (zh) * | 2011-05-27 | 2012-11-28 | 中国石油天然气集团公司 | 一种基于稀疏反演的滑动扫描谐波压制方法 |
CN103885085A (zh) * | 2012-12-21 | 2014-06-25 | 中国石油集团东方地球物理勘探有限责任公司 | 一种压制可控震源谐波干扰的方法 |
CN104111477A (zh) * | 2014-07-04 | 2014-10-22 | 中国石油天然气集团公司 | 一种用于压制谐波的滑动时间获取方法及系统 |
CN105277987A (zh) * | 2015-05-27 | 2016-01-27 | 中石化石油工程地球物理有限公司胜利分公司 | 基于预测滤波法和纯相移法的可控震源谐波压制方法 |
-
2017
- 2017-06-14 CN CN201710448762.1A patent/CN107255835B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6161076A (en) * | 1997-11-14 | 2000-12-12 | Baker Hughes Incorporated | Seismic data acquisition and processing using non-linear distortion in a vibratory output signal |
CN102692649A (zh) * | 2011-03-23 | 2012-09-26 | 中国石油天然气集团公司 | 一种可控震源地震勘探中消除谐波干扰的方法 |
CN102798894A (zh) * | 2011-05-27 | 2012-11-28 | 中国石油天然气集团公司 | 一种基于稀疏反演的滑动扫描谐波压制方法 |
CN103885085A (zh) * | 2012-12-21 | 2014-06-25 | 中国石油集团东方地球物理勘探有限责任公司 | 一种压制可控震源谐波干扰的方法 |
CN104111477A (zh) * | 2014-07-04 | 2014-10-22 | 中国石油天然气集团公司 | 一种用于压制谐波的滑动时间获取方法及系统 |
CN105277987A (zh) * | 2015-05-27 | 2016-01-27 | 中石化石油工程地球物理有限公司胜利分公司 | 基于预测滤波法和纯相移法的可控震源谐波压制方法 |
Non-Patent Citations (1)
Title |
---|
可控震源两种谐波产生机理与特征研究;李振春 等;《石油物探》;20160331;第55卷(第2期);第159-172页 |
Also Published As
Publication number | Publication date |
---|---|
CN107255835A (zh) | 2017-10-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101551465B (zh) | 一种自适应识别和消除地震勘探单频干扰的方法 | |
Atkinson et al. | Attenuation and source characteristics of the 23 June 2010 M 5.0 Val-des-Bois, Quebec, earthquake | |
CN102692649A (zh) | 一种可控震源地震勘探中消除谐波干扰的方法 | |
CN108549100A (zh) | 基于非线性高次拓频的时间域多尺度全波形反演方法 | |
CN104950326B (zh) | 基于目的层频谱的可控震源非线性扫描信号的设计方法 | |
RU2011135740A (ru) | Способ обнаружения или мониторинга структуры размером с углеводородный пласт-коллектор | |
CN112327358B (zh) | 一种粘滞性介质中声波地震数据正演模拟方法 | |
Pavlenko | Nonlinear seismic effects in soils: numerical simulation and study | |
CN102706555A (zh) | 一种复解析最优小波解调法 | |
US2808577A (en) | Method of and apparatus for determining the travel time of signals | |
CN107229075B (zh) | 深度域地震子波的确定方法和装置 | |
CN108957551B (zh) | 基于重构地面力信号的可控震源谐波压制方法 | |
CN110879410A (zh) | 一种多分量地震面波勘探方法 | |
CN107255835B (zh) | 一种用于压制机械因素谐波干扰的方法 | |
CN111427086B (zh) | 利用纵横波振幅频率分量预测地层最大衰减参数的方法 | |
JP3467208B2 (ja) | 表面波位相速度検出システム、及び表面波の位相速度検出方法 | |
CN105277987A (zh) | 基于预测滤波法和纯相移法的可控震源谐波压制方法 | |
CN108318919B (zh) | 一种动态参数的可控震源非线性扫描信号设计系统及方法 | |
CN109212609B (zh) | 基于波动方程延拓的近地表噪音压制方法 | |
CN114779328B (zh) | 一种基于自相关噪音成像的地球内部结构探测方法 | |
CN113031070B (zh) | 一种深度域合成地震记录的制作方法 | |
CN109143345A (zh) | 基于模拟退火的品质因子q非线性反演方法及系统 | |
CN106646599A (zh) | 针对地表响应因素产生谐波的自动识别与衰减方法 | |
Zhang et al. | Extraction Method of Vibroseis Phase Signal Based on Time-Varying Narrowband Filtering | |
CN112526611A (zh) | 表层地震波品质因子的提取方法及装置 |
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 |