CN102426386A - 脉冲编码可控震源的多维匹配冲击方法 - Google Patents
脉冲编码可控震源的多维匹配冲击方法 Download PDFInfo
- Publication number
- CN102426386A CN102426386A CN2011103298964A CN201110329896A CN102426386A CN 102426386 A CN102426386 A CN 102426386A CN 2011103298964 A CN2011103298964 A CN 2011103298964A CN 201110329896 A CN201110329896 A CN 201110329896A CN 102426386 A CN102426386 A CN 102426386A
- Authority
- CN
- China
- Prior art keywords
- signal
- multidimensional
- coupling
- vibroseis
- seismic
- 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
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种脉冲编码可控震源的多维匹配冲击方法。依据多维匹配冲击的信号流程首先产生4个相互关联的匹配冲击信号,分别作为脉冲编码可控震源的激发信号于同等采集条件下依次激发得到4个匹配的原始地震记录,而后依据多维匹配冲击的信号流程对4个匹配的原始地震记录进行解码叠加运算得到脉冲编码可控震源多维匹配冲击的综合解码地震剖面。本发明与脉冲编码可控震源的常规线性扫描冲击方法相比,能够进一步提高地震信号的信噪比,压制地震剖面中的相关噪声和随机噪声干扰,是一种有效的脉冲编码可控震源编码方案。
Description
技术领域
本发明涉及一种地球物理勘探方法,尤其是适用于以脉冲编码可控震源多维匹配冲击作为震源进行激发的高分辨率地震勘探。
背景技术
震源是地震勘探中信号的源头,其信号质量对于地震记录的分辨率和信噪比具有决定性作用。炸药震源是一种陆地地震勘探中的有效震源形式,然而其对于环境的破坏和污染等弊端已无法适应城市及工程等领域中对于非破坏性地震勘探的要求。可控震源作为一种非破坏性震源形式在陆地地震勘探中得到了广泛的应用,其中具有代表性的是基于Vibroseis技术的连续振动方式可控震源。Vibroseis可控震源以Chirp信号作为震源扫描信号以线性扫频的形式向地下引入连续的地震波,其原始地震记录经相关解码等运算后可压缩为单脉冲形式的地震剖面。然而由于Chirp信号的自相关函数旁瓣能量较强,使得Vibroseis可控震源Chirp扫描的地震剖面中存在明显的旁瓣效应,降低了地震信号的分辨率。
为了压制可控震源Chirp线性扫描地震剖面中的旁瓣效应,数字通信领域中的伪随机编码方案被应用到连续振动方式可控震源的信号设计中。其中二元m-序列伪随机扫描方法可以有效压制Chirp扫描所产生的旁瓣效应,然而由于相关运算的作用而使地震剖面中出现了能量较强的相关噪声,降低了地震信号的信噪比。基于连续振动方式可控震源的匹配扫描方法依次激发两个相互匹配的伪随机扫描信号,通过将采集到的两个匹配的原始地震记录进行解码和叠加运算,能够有效压制地震响应剖面中的旁瓣效应和相关噪声干扰,同时提高了地震信号的分辨率和信噪比,是一种有效的连续振动方式可控震源编码方式。Vibroseis可控震源除了旁瓣效应及相关噪声等干扰之外,由于可控震源基板与地面耦合的频率响应的非线性以及耦合条件的差异等因素,容易使震源信号产生畸变,在解码地震剖面中出现的谐波干扰以常规数据处理方法难以消除。
基于编码冲击方式的脉冲编码可控震源不同于连续振动方式的Vibroseis可控震源,其结合了Vibroseis技术线性扫频的思想和Mini-SOSIE技术随机冲击的思想,按照一定的时间间隔规律向地下激发扫描冲击序列,各次冲击之间的时间间隔按线性规律变化。扫描冲击方式的脉冲编码可控震源不受连续振动方式Vibroseis可控震源耦合条件的限制和影响,所产生的冲击序列信号具有良好的稳定性和一致性,近年来在浅层高分辨率地震勘探等领域中得到了广泛的应用和不断的完善。目前适用于脉冲编码可控震源的编码方案为线性扫描冲击编码方案。
发明内容
本发明的目的就在于针对上述现有技术的不足,提共一种适用于冲击式脉冲编码可控震源的多维匹配冲击方法。
该方法通过依次激发4个相关联的多维匹配冲击信号,并对采集得到的4个相互关联的你原始地震记录进行解码叠加等运算来进一步提高地震信号的信噪比,压制解码地震剖面中出现的相关噪声和随机噪声等干扰信号。
本发明的目的是通过以下技术方案实现的:
①首先确定用于产生多维匹配冲击信号的匹配伪随机序列偶(ak}和{bk}的长度N。
②依次构造由“1”、“-1”所构成且长度为N的两个相关联的匹配伪随机序列偶(ak)和{bk)。
③以正弦载波信号对两个匹配伪随机序列偶(ak}和{bk}进行调制产生连续振动方式可控震源的匹配扫描信号pa(t)和pb(t)。
④将pa(t)和pb(t)中以0相位起始的各正弦波保留,以p相位起始的各正弦波置零得到两个子信号pa1(t)和pb1(t);将pa(t)和pb(t)中以0相位起始的各正弦波置零,以p相位起始的各正弦波保留并做反相处理得到两个子信号pa2(t)和pb2(t)。
⑤将4个子信号pa1(t)、pb1(t)、pa2(t)、pb2(t)各正弦波周期中的单个正弦波对应一个元素“1”,零信号对应一个元素“0”,得到4个由“0”、“1”构成的多维匹配伪随机子序列{ak1}、{bk1}、{ak2}、{bk2}。
⑥以地震子波信号对多维匹配伪随机子序列{ak1}、{bk1}、{ak2}、{bk2}进行调制得到可由脉冲编码可控震源激发的4个多维匹配冲击信号qa1(t)、qb1(t)、qa2(t)、qb2(t)。
⑦以多维匹配冲击信号qa1(t)、qb1(t)、qa2(t)、qb2(t)分别作为脉冲编码可控震源的激发信号进行同等条件下4次激发,采集得到4个匹配的原始地震记录Ra1(xi,t)、Rb1(xi,t)、Ra2(xi,t)、Rb2(xi,t),(其中i=1,2,3,...,M,M为地震道数)。
⑧将(Ra1(xi,t)-Ra2(xi,t))和(Rb1(xi,t)-Rb2(xi,t))分别进行解码运算得到两个匹配的解码地震剖面Ca(xi,t)和Cb(xi,t)。
⑨最后对两个匹配的解码地震剖面Ca(xi,t)和Cb(xi,t)进行叠加运算得到脉冲编码可控震源多维匹配冲击的综合解码地震剖面C(xi,t)。
有益效果:本发明为了能够对比脉冲编码可控震源的多维匹配冲击方法相对于常规线性扫描冲击方法所具有的优势,在相同的地震波采集系统及采集条件下先后进行了大锤震源激发和脉冲编码可控震源的线性扫描冲击以及多维匹配冲击。在图6(a)所示的大锤震源激发的地震剖面中,远偏移距地震道中存在严重的随机噪声干扰,淹没了有效的地震信号,信噪比较低;在图6(b)所示的脉冲编码可控震源线性扫描冲击的解码地震剖面中,远偏移距地震道中的随机噪声干扰相比大锤震源激发情形得到了明显的衰减,然而其剖面下部存在线状相关噪声干扰;在图6(c)所示的脉冲编码可控震源多维匹配冲击的综合解码地震剖面中,同样能够有效压制远偏移距地震道中的随机噪声干扰,同时相比线性扫描冲击情形,相关噪声的能量要微弱得多,使得地震剖面的信噪比得到了进一步的提高。本发明的多维匹配冲击方法是一种有效的脉冲编码可控震源编码方案,对于现有的线性扫描冲击可控震源,仅需更改震源控制器中的编码信号,即可实现脉冲编码可控震源的多维匹配冲击过程。
附图说明
图1:脉冲编码可控震源多维匹配冲击信号流程图。
图2:连续振动方式可控震源匹配扫描信号pa(t)和pb(t)信号波形图。
图3:匹配扫描子信号pa1(t)、pb1(t)、pa2(t)、pb2(t)信号波形图。
图4:多维匹配冲击信号qa1(t)、qb1(t)、qa2(t)、qb2(t)信号波形图。
图5:多维匹配冲击原始地震记录(截取5s~15s部分)。
图6:(a)大锤震源激发的原始地震剖面;
(b)脉冲编码可控震源线性扫描冲击的解码地震剖面;
(c)脉冲编码可控震多维匹配冲击的综合解码地震剖面C(xi,t)。
具体实施方式
下面结合附图和实施例作进一步的详细说明:
按照图1中的多维匹配冲击信号流程图对脉冲编码可控震源多维匹配冲击的完整过程做进一步详细说明。
①首先确定匹配伪随机序列偶{ak}和{bk}的长度。为叙述方便,取N=10。
②依次构造由“1”、“-1”构成且长度N=10的两个相关联的匹配伪随机序列偶{ak}和{bk}为
{ak}={1,-1,-1,1,-1,1,-1,-1,-1,1},
{bk}={1,-1,-1,-1,-1,-1,-1,1,1,-1}。
③以频率f=50Hz的正弦载波信号对两个匹配伪随机序列偶{ak}和{bk}进行调制产生连续振动方式可控震源的匹配扫描信号pa(t)和pb(t),其信号波形图如图2所示。
④将信号pa(t)和pb(t)中以0相位起始的各正弦波保留,以p相位起始的各正弦波置零得到两个子信号pa1(t)和pb1(t);将pa(t)和pb(t)中以0相位起始的各正弦波置零,以p相位起始的各正弦波保留并做反相处理得到两个子信号pa2(t)和pb2(t),4个子信号pa1(t)、pb1(t)、pa2(t)、pb2(t)的波形图如图3所示。
⑤将4个子信号pa1(t)、pb1(t)、pa2(t)、pb2(t)中各正弦波周期中单个正弦波对应一个元素“1”,零信号对应一个元素“0”,得到4个由“0”、“1”构成的多维匹配伪随机子序列{ak1}、{bk1}、{ak2}、{bk2}为
{ak1}={1,0,0,1,0,1,0,0,0,1},
{bk1}={1,0,0,0,0,0,0,1,1,0},
{ak2}={0,1,1,0,1,0,1,1,1,0},
{bk2}={0,1,1,1,1,1,1,0,0,1}。
⑥以主频为90Hz的Ricker子波作为地震子波信号对多维匹配伪随机子序列{ak1}、{bk1}、{ak2}、{bk2}进行调制得到可由脉冲编码可控震源进行激发的4个多维匹配冲击信号qa1(t)、qb1(t)、qa2(t)、qb2(t),信号波形图如图4所示。
⑦以多维匹配冲击信号qa1(t)、qb1(t)、qa2(t)、qb2(t)分别作为脉冲编码可控震源的激发信号进行同等条件下4次激发采集得到4个匹配的原始地震记录Ra1(xi,t)、Rb1(xi,t)、Ra2(x,t)、Rb2(xi,t),(其中1#i 28,i代表各地震道)。图5显示了Ra1(xi,t)、Rb1(xi,t)、Ra2(xi,t)、Rb2(xi,t)的5s~15s之间的原始地震记录。
⑧将(Ra1(xi,t)-Ra2(xi,t))和(Rb1(xi,t)-Rb2(xi,t))分别进行解码运算得到两个匹配的解码地震剖面Ca(xi,t)和Cb(xi,t)。
⑨最后对两个匹配的解码地震剖面Ca(xi,t)和Cb(xi,t)进行叠加运算得到脉冲编码可控震源多维匹配冲击的综合解码地震剖面C(xi,t),其地震剖面如图6(c)所示。
Claims (1)
1.一种脉冲编码可控震源的多维匹配冲击方法,其特征在于:通过使脉冲编码可控震源在相同采集条件下依次激发4个相互关联的匹配冲击信号来完成多维匹配冲击过程,多维匹配冲击方法完整的信号流程是通过以下9个步骤来顺序实现:
a、首先确定用于产生多维匹配冲击信号的匹配伪随机序列偶{ak}和{bk}的长度N;
b、依次构造由“1”、“-1”所构成且长度为N的两个相关联的匹配伪随机序列偶{ak}和{bk};
c、以正弦载波信号对两个匹配伪随机序列偶{ak}和{bk}进行调制产生连续振动方式可控震源的匹配扫描信号pa(t)和pb(t);
d、将pa(t)和pb(t)中以0相位起始的各正弦波保留,以p相位起始的各正弦波置零得到两个子信号pa1(t)和pb1(t);将pa(t)和pb(t)中以0相位起始的各正弦波置零,以p相位起始的各正弦波保留并做反相处理得到两个子信号pa2(t)和pb2(t);
e、将4个子信号pa1(t)、pb1(t)、pa2(t)、pb2(t)各正弦波周期中的单个正弦波对应一个元素“1”,零信号对应一个元素“0”,得到4个由“0”、“1”构成的多维匹配伪随机子序列{ak1}、{bk1}、{ak2}、{bk2};
f、以地震子波信号对多维匹配伪随机子序列{ak1}、{bk1}、{ak2}、{bk2}进行调制得到由脉冲编码可控震源激发的4个多维匹配冲击信号qa1(t)、qb1(t)、qa2(t)、qb2(t);
g、以多维匹配冲击信号qa1(t)、qb1(t)、qa2(t)、qb2(t)分别作为脉冲编码可控震源的激发信号进行同等条件下4次激发,采集得到4个匹配的原始地震记录Ra1(xi,t)、Rb1(xi,t)、Ra2(xi,t)、Rb2(xi,t),(其中i=1,2,3,...,M,M为地震道数)
h、将(Ra1(xi,t)-Ra2(xi,t))和(Rb1(xi,t)-Rb2(xi,t))分别进行解码运算得到两个匹配的解码地震剖面Ca(xi,t)和Cb(xi,t);
i、最后对两个匹配的解码地震剖面Ca(xi,t)和Cb(xi,t)进行叠加运算得到脉冲编码可控震源多维匹配冲击的综合解码地震剖面C(xi,t)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110329896 CN102426386B (zh) | 2011-10-27 | 2011-10-27 | 脉冲编码可控震源的多维匹配冲击方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110329896 CN102426386B (zh) | 2011-10-27 | 2011-10-27 | 脉冲编码可控震源的多维匹配冲击方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102426386A true CN102426386A (zh) | 2012-04-25 |
CN102426386B CN102426386B (zh) | 2013-06-26 |
Family
ID=45960388
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110329896 Expired - Fee Related CN102426386B (zh) | 2011-10-27 | 2011-10-27 | 脉冲编码可控震源的多维匹配冲击方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102426386B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102879825A (zh) * | 2012-09-18 | 2013-01-16 | 吉林大学 | 可控震源地震数据的强脉冲噪声检测及压制方法 |
CN104035128A (zh) * | 2013-03-06 | 2014-09-10 | 中国石油集团东方地球物理勘探有限责任公司 | 可控震源伪随机扫描信号生成方法 |
CN111856552A (zh) * | 2019-04-30 | 2020-10-30 | 中国石油天然气集团有限公司 | 可控震源扫描信号的生成方法及装置 |
CN112505747A (zh) * | 2020-12-22 | 2021-03-16 | 吉林大学 | 基于多信号发生器协同可控震源振动畸变抑制系统及方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101251602A (zh) * | 2008-04-15 | 2008-08-27 | 吉林大学 | 海洋组合可控震源短时分段扫描方法 |
CN101545981A (zh) * | 2008-03-28 | 2009-09-30 | 中国石油天然气集团公司 | 可控震源地震数据零相位子波最小相位化方法 |
-
2011
- 2011-10-27 CN CN 201110329896 patent/CN102426386B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101545981A (zh) * | 2008-03-28 | 2009-09-30 | 中国石油天然气集团公司 | 可控震源地震数据零相位子波最小相位化方法 |
CN101251602A (zh) * | 2008-04-15 | 2008-08-27 | 吉林大学 | 海洋组合可控震源短时分段扫描方法 |
Non-Patent Citations (1)
Title |
---|
林君等: "用于浅层地震勘查的高频可控震源系统", 《地质装备》, vol. 2, no. 3, 30 September 2001 (2001-09-30), pages 20 - 26 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102879825A (zh) * | 2012-09-18 | 2013-01-16 | 吉林大学 | 可控震源地震数据的强脉冲噪声检测及压制方法 |
CN104035128A (zh) * | 2013-03-06 | 2014-09-10 | 中国石油集团东方地球物理勘探有限责任公司 | 可控震源伪随机扫描信号生成方法 |
CN104035128B (zh) * | 2013-03-06 | 2016-08-03 | 中国石油集团东方地球物理勘探有限责任公司 | 可控震源伪随机扫描信号生成方法 |
CN111856552A (zh) * | 2019-04-30 | 2020-10-30 | 中国石油天然气集团有限公司 | 可控震源扫描信号的生成方法及装置 |
CN111856552B (zh) * | 2019-04-30 | 2023-07-25 | 中国石油天然气集团有限公司 | 可控震源扫描信号的生成方法及装置 |
CN112505747A (zh) * | 2020-12-22 | 2021-03-16 | 吉林大学 | 基于多信号发生器协同可控震源振动畸变抑制系统及方法 |
CN112505747B (zh) * | 2020-12-22 | 2021-10-01 | 吉林大学 | 基于多信号发生器协同可控震源振动畸变抑制系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102426386B (zh) | 2013-06-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU2010213635B2 (en) | Multiple seismic signal inversion | |
Quiros et al. | Seismic interferometry of railroad induced ground motions: Body and surface wave imaging | |
US11774618B2 (en) | Marine vibrator directive source survey | |
US9551798B2 (en) | Seismic vibrator to produce a continuous signal | |
Feng et al. | Effectiveness of multi-mode surface wave inversion in shallow engineering site investigations | |
AU2017248407B2 (en) | Coded signals for marine vibrators | |
CN102426386B (zh) | 脉冲编码可控震源的多维匹配冲击方法 | |
CA2682087A1 (en) | Method for optimizing energy output from a seismic vibrator array | |
CN106291684B (zh) | 一种盲源地震波场的地震响应恢复与虚源道集构建方法 | |
CN104216010A (zh) | 利用可控震源谐波提高地震数据质量的方法 | |
CN101251602B (zh) | 海洋组合可控震源短时分段扫描方法 | |
Cheng et al. | Testing of a permanent orbital surface source and distributed acoustic sensing for monitoring of unconventional reservoirs: Preliminary results from the Eagle Ford Shale | |
CN102508292B (zh) | 可控震源匹配扫描方法 | |
Herbst et al. | Seismic sources for shallow investigations: A field comparison from Northern Germany | |
Baradello et al. | Vibroseis deconvolution: A comparison of pre and post correlation vibroseis deconvolution data in real noisy data | |
CN105319583A (zh) | 基于分频动态编码的可控震源混叠数据成像方法 | |
Pugin et al. | Update on recent observations in multi-component seismic reflection profiling | |
Wang et al. | Notch effect and frequency compensation of dual-sensor OBC data in shallow water | |
CN103018772B (zh) | 震电勘探信号的去噪处理方法 | |
Miller et al. | Enhancing high-frequency vibroseis images of near surface: Uncorrelated is key | |
Li et al. | Diffraction Separation for the Ground Penetrating Radar Data by Masking Filters | |
Pecholcs et al. | Distance Separated Custom Slip-Sweep-A New High-Productivity Method | |
Naskar et al. | A Superior Seismic Data Preprocessing Technique to Improve the Resolution of Surface wave Velocity Spectrums | |
RU2388020C1 (ru) | Способ получения полевых сейсмограмм, свободных от многократных волн | |
Hanafy et al. | Locating trapped miners using super-stacking and super-resolution properties of time reversal mirrors |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130626 Termination date: 20131027 |