CN104570121B - 定向地震波畸变信号消除方法 - Google Patents
定向地震波畸变信号消除方法 Download PDFInfo
- Publication number
- CN104570121B CN104570121B CN201410832029.6A CN201410832029A CN104570121B CN 104570121 B CN104570121 B CN 104570121B CN 201410832029 A CN201410832029 A CN 201410832029A CN 104570121 B CN104570121 B CN 104570121B
- Authority
- CN
- China
- Prior art keywords
- seismic
- data
- road
- record
- filtering
- 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.)
- Expired - Fee Related
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种定向地震波畸变信号消除方法,利用地震波束定向后的记录与原始中心记录的局部相关系数,对地震波束定向后的记录进行滤波,即可消除主波束方向外的畸变信号,得到无畸变的地震波束定向记录。本发明与现有技术相比能够在原始数据信噪比低的情况下,在保证定向地震波数据在勘探目标上形成高质量数据情况下,同时有效消除了主波束方向外的畸变信号。该方法对强背景噪声干扰条件下,提高地震勘探数据质量具有重要的应用价值。
Description
技术领域:
本发明涉及一种地震数据畸变信号的消除方法,尤其是针对低信噪比条件下的定向地震波数据主波束方向外畸变信号的消除方法。
背景技术:
波束定向的思想最早源于相控雷达领域。由于波束定向的方法能够有效加强目标体上的有效信号,很快被引进地震勘探领域。地震波束形成最初在震源端实现,之后提出了基于接收端的地震波束合成方法,涉及的相关专利主要有:
CN101320095公开了一种《单震源定向照明地震信号合成方法》;CN101339252《基于单震源的定向照明地震勘探方法》;CN102590861A《基于接收阵列的地震波束形成方法》;CN103984019A《局部相关加权地震波束合成方法》。
地震波束形成方法使无方向的地震波场变成具有方向性的波场,主波束方向上的有效信号能量得到加强的同时,主波束方向外的波场相互混淆,引起信号畸变。这种畸变会使后期地震处理和解释出现较大的偏差。因此,在通过定向地震波技术提高勘探目标信号信噪比的同时,有效地消除地震波非主波束的信号畸变,成为地震波束定向技术应用于实际地震勘探工作的关键。
CN103984019A《局部相关加权地震波束合成方法》提出了一种消除定向地震波畸变的方法,该方法工作原理是在合成定向地震波过程中同时抑制畸变,但该方法对原始数据信噪比要求较高,计算量大,在原始数据信噪比低的条件下存在较大的误差。
发明内容:
本发明的目的就在于针对上述地震数据信噪比较低的情况,在合成定向地震波地震波束数据后,对主波束方向外为信号产生的畸变问题,提供一种定向地震波数据畸变信号变消除方法。
本发明的主要思想是:定向地震波束具有在主波束方向上信号被加强的特点,但是主波束方向外由于信号不是同相叠加,因此会产生畸变形成干扰。主波束方向上的信号与原始波束形成的中心记录(延时为零)上的信号具有很强的局部相关性,局部相关系数大,而主波束方向外的干扰与原始中心记录上信号的局部相关性较弱,局部相关系数较小。因而利用地震波束定向后的记录与原始中心记录的局部相关系数,对地震波束定向后的记录进行滤波,即可消除主波束方向外的畸变信号,得到无畸变的地震波束定向记录。
本发明的目的是通过以下技术方案实现的:
定向地震波畸变信号消除方法,包括以下步骤:
a、输入n炮炮点连续的地震记录,记为U1(t,x)、U2(t,x)、…、Um(t,x)、…、Un-1(t,x)、Un(t,x),t为时间,x为道号,n取为大于2的奇数,其中
b、保持第m炮记录Um(t,x)不变,其他炮记录利用延时参数τ1和道分选参数λ,依次作等间隔延时和等道数分选,延时和分选结果为U1(t-(m-1)τ1,x+(m-1)λ)、U2(t-(m-2)τ1,x+(m-2)λ)、…、Um(t,x)、…、Un-1(t+(m-2)τ1,x-(m-2)λ)、Un(t+(m-1)τ1,x-(m-1)λ),其中d1为炮间距,d2为道间距,k取正整数;
c、对延时后的所有记录作线性叠加合成地震波束记录。波束合成后的数据记为D(t)。则
d、以Um(t,x)为基准,取数据D(t)。记Um(t,x)的第一道数据为u1,D(t)的第一道数据为s1;
e、以时窗T选取u1和s1的第一个时窗数据u1,1(t)和s1,1(t),T不小于一个有效地震子波信号的周期;
f、将u1,1(t)和s1,1(t)作互相关,计算得到R1,1(τ)=∫u1,1(t)s1,1(t-τ)dt,最大互相关系数c(1,1)=max[R1,1(τ)]=R1,1[τc(1,1)],τc为最大互相关系数对应的延时量;
g、对u1和s1,同时沿第一道移动时窗T,重复步骤d~f,可得第一道数据的最大互相关系数曲线c1=(c(1,1),c(1,2),…,c(1,v))T及其对应的延时曲线τc(1)=(τc(1,1),τc(1,2),…,τc(1,v))T,v为一道数据总的采样点数,(·)T表示转置;
h、对Um(t,x)和D(t)对应的其他道数据,重复步骤d~g,可得最大互相关系数谱C=(c1,c2,c3,…cσ)及其对应的延时谱τc=(τc(1),τc(2),τc(3),…τc(σ)),σ为总道数;
i、对C进行延时滤波,延时滤波采用如下公式: 为延时滤波后的互相关系数谱,f为地震子波主频;
j、对K进行单道阈值滤波,单道阈值滤波采用如下公式: 为阈值,δ为滤波因子取值在1~8,i=1,2,3,…,σ,j=1,2,3,…,v;对K的所有道进行单道阈值滤波,就得到单道阈值滤波后的结果
k、以∏中第q道和第p个采样点为中心,对∏进行l道(l为奇数)横向中值滤波,l道横向中值滤波的输出为ρ(p,q)=median[π(i,j)],i=p;j=q-(l-1)/2,…,q+(l-1)/2;i、p为采样点号,j为道号,q为中心道号,median(·)为取中值;对∏中的每个采样点都进行l道横向中值滤波,即可得到中值滤波后的结果
1、最后对P进行编码,得到最终的主波束信号提取系数W,编码采用如下公式:
m、消除D(t)中畸变信号后保留主波束方向上的有效信号De(t)=WD(t)。有益效果:经试验,本发明公开的定向地震波数据的畸变信号消除方法,能够在原始数据信噪比低的情况下,在保证定向地震波数据在勘探目标上形成高质量数据情况下,同时有效消除了主波束方向外的畸变信号。该方法对强背景噪声干扰条件下,提高地震勘探数据质量具有重要的应用价值。
附图说明:
图1地震波束定向后的地震记录
图2消除畸变信号后的地震记录
具体实施方式:
下面结合附图和实施例做进一步的详细说明:
在本实施例中以n=9,τ1=1ms,λ=1,信号主频f=60Hz为例对9炮地震记录进行波束定向和畸变信号的消除,但该方法不受实例中给出的参数的限制。
定向地震波畸变信号消除方法,包括以下步骤:
a、输入n=9炮炮点连续的地震记录,记为U1(t,x)、U2(t,x)、…、U5(t,x)、…、U8(t,x)、U9(t,x);
b、保持第m=5炮记录U5(t,x)不变(延时为0),其他炮记录利用延时参数τ1=1ms和道分选参数λ=1(d1=d2k取1),依次作等间隔延时和等道数分选,延时和分选结果为U1(t-4τ1,x+4λ)、U2(t-3τ1,x+3λ)、…、U5(t,x)、…、U8(t+3τ1,x-3λ)、U9(t+4τ1,x-4λ);
c、对延时后的所有记录作线性叠加合成地震波束记录。波束合成后的数据记为D(t),则 如图1所示;
d、以U5(t,x)为基准,取数据D(t),记U5(t,x)的第一道数据为u1,D(t)的第一道数据为s1;
e、以时窗T=40ms选取u1和s1的第一个时窗数据u1,1(t)和s1,1(t);
f、将u1,1(t)和s1,1(t)作互相关,计算得到R1,1(τ)=∫u1,1(t)s1,1(t-τ)dt,最大互相关系数c(1,1)=max[R1,1(τ)]=R1,1[τc(1,1)],τc为最大互相关系数对应的延时量;
g、对u1和s1,同时沿第一道移动时窗T,重复步骤d~f,可得第一道数据的最大互相关系数曲线c1=(c(1,1),c(1,2),…,c(1,100))T及其对应的延时曲线τc(1)=(τc(1,1),τc(1,2),…,τc(1,100))T;
h、对U5(t,x)和D(t)的其他道数据,重复步骤d~g,可得最大互相关系数谱C及其对应的延时谱τc;
i、对C进行延时滤波,延时滤波采用如下公式: 为延时滤波后的互相关系数谱;
j、对K进行单道阈值滤波,单道阈值滤波采用如下公式: 为阈值,δ=6,i为道号;对K的所有道进行单道阈值滤波,就得到单道阈值滤波后的结果
k、对∏中所有采样点进行7道横向中值滤波。以第q道和第p个采样点为中心的7道横向中值滤波的输出为ρ(p,q)=median[π(i,j)],i=p;j=q-3,…,q+3;i、p为采样点号,j为道号,q为中心道号;对∏中的每个采样点都进行7道横向中值滤波,得到中值滤波后的结果
1、最后对P进行编码,得到最终的主波束信号提取系数W,编码采用如下公式:
m、消除D(t)中畸变信号后保留主波束方向上的有效信号De(t)=WD(t),如图2所示。
Claims (1)
1.一种定向地震波畸变信号消除方法,包括以下步骤:
a、输入n炮炮点连续的地震记录,记为U1(t,x)、U2(t,x)、…、Um(t,x)、…、Un-1(t,x)、Un(t,x),t为时间,x为道号,n取为大于2的奇数,其中
b、保持第m炮记录Um(t,x)不变,其他炮记录利用延时参数τ1和道分选参数λ,依次作等间隔延时和等道数分选,延时和分选结果为
U1(t-(m-1)τ1,x+(m-1)λ)、U2(t-(m-2)τ1,x+(m-2)λ)、…、Um(t,x)、…、Un-1(t+(m-2)τ1,x-(m-2)λ)、Un(t+(m-1)τ1,x-(m-1)λ),其中d1为炮间距,d2为道间距,k取正整数;
c、对延时后的所有记录作线性叠加合成地震波束记录,波束合成后的数据记为D(t),则
d、以Um(t,x)为基准,取数据D(t),记Um(t,x)的第一道数据为u1,D(t)的第一道数据为s1;
e、以时窗T选取u1和s1的第一个时窗数据u1,1(t)和s1,1(t),T不小于一个有效地震子波信号的周期;
f、将u1,1(t)和s1,1(t)作互相关,计算得到R1,1(τ)=∫u1,1(t)s1,1(t-τ)dt,最大互相关系数c(1,1)=max[R1,1(τ)]=R1,1[τc(1,1)],τc为最大互相关系数对应的延时量;
g、对u1和s1,同时沿第一道移动时窗T,重复步骤d~f,得第一道数据的最大互相关系数曲线c1=(c(1,1),c(1,2),…,c(1,v))T及其对应的延时曲线τc(1)=(τc(1,1),τc(1,2),…,τc(1,v))T,v为一道数据总的采样点数,(·)T表示转置;
h、对Um(t,x)和D(t)对应的其他道数据,重复步骤d~g,可得最大互相关系数谱C=(c1,c2,c3,…cσ)及其对应的延时谱τc=(τc(1),τc(2),τc(3),…τc(σ)),σ为总道数;
i、对C进行延时滤波,延时滤波采用如下公式:
i=1,2,3,…,σ,j=1,2,3,…,v,为延时滤波后的互相关系数谱,f为地震子波主频;
j、对Κ进行单道阈值滤波,单道阈值滤波采用如下公式: 为阈值,δ为滤波因子取值在1~8,i=1,2,3,…,σ,j=1,2,3,…,v;对Κ的所有道进行单道阈值滤波,就得到单道阈值滤波后的结果
k、以Π中第q道和第p个采样点为中心,对Π进行l道(l为奇数)横向中值滤波,l道横向中值滤波的输出为ρ(p,q)=median[π(i,j)],i=p;j=q-(l-1)/2,…,q+(l-1)/2;i、p为采样点号,j为道号,q为中心道号,median(·)为取中值;对Π中的每个采样点都进行l道横向中值滤波,即可得到中值滤波后的结果
l、最后对Ρ进行编码,得到最终的主波束信号提取系数W,编码采用如下公式:i=1,2,3,…,σ,j=1,2,3,…,v,
m、消除D(t)中畸变信号后保留主波束方向上的有效信号De(t)=WD(t)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410832029.6A CN104570121B (zh) | 2014-12-26 | 2014-12-26 | 定向地震波畸变信号消除方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410832029.6A CN104570121B (zh) | 2014-12-26 | 2014-12-26 | 定向地震波畸变信号消除方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104570121A CN104570121A (zh) | 2015-04-29 |
CN104570121B true CN104570121B (zh) | 2017-01-04 |
Family
ID=53086635
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410832029.6A Expired - Fee Related CN104570121B (zh) | 2014-12-26 | 2014-12-26 | 定向地震波畸变信号消除方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104570121B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0289670A1 (en) * | 1987-05-05 | 1988-11-09 | Atlantic Richfield Company | Method of suppressing vibration seismic signal correlation noise |
CN101320095A (zh) * | 2008-07-22 | 2008-12-10 | 吉林大学 | 单震源定向照明地震信号合成方法 |
CN101339252A (zh) * | 2008-08-08 | 2009-01-07 | 吉林大学 | 基于单震源的定向照明地震勘探方法 |
CN102590861A (zh) * | 2012-01-17 | 2012-07-18 | 吉林大学 | 基于接收阵列的地震波束形成方法 |
CN103984019A (zh) * | 2014-06-07 | 2014-08-13 | 吉林大学 | 局部相关加权地震波束合成方法 |
-
2014
- 2014-12-26 CN CN201410832029.6A patent/CN104570121B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0289670A1 (en) * | 1987-05-05 | 1988-11-09 | Atlantic Richfield Company | Method of suppressing vibration seismic signal correlation noise |
CN101320095A (zh) * | 2008-07-22 | 2008-12-10 | 吉林大学 | 单震源定向照明地震信号合成方法 |
CN101339252A (zh) * | 2008-08-08 | 2009-01-07 | 吉林大学 | 基于单震源的定向照明地震勘探方法 |
CN102590861A (zh) * | 2012-01-17 | 2012-07-18 | 吉林大学 | 基于接收阵列的地震波束形成方法 |
CN103984019A (zh) * | 2014-06-07 | 2014-08-13 | 吉林大学 | 局部相关加权地震波束合成方法 |
Non-Patent Citations (3)
Title |
---|
余弦振幅加权定向地震波束形成方法;宋健等;《仪器仪表学报》;20140831;第35卷(第8期);第1729-1737页 * |
相控震源与组合震源信号检测结果对比;姜弢等;《控制与决策》;20041231;第19卷(第12期);第1345-1348、1353页 * |
相控震源对地震信号信噪比的改善研究;姜弢等;《地球物理学报》;20061130;第49卷(第6期);第1819-1825页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104570121A (zh) | 2015-04-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103901416B (zh) | 一种基于稳健主成分分析法的多通道杂波抑制方法 | |
CN102819043B (zh) | 阵列信号随机噪声自适应模型去噪方法 | |
CN104199001B (zh) | 一种认知雷达抗速度欺骗干扰的相位编码方法 | |
CN103543469B (zh) | 一种基于小波变换的小尺度阈值去噪方法 | |
CN101598812B (zh) | 去除数字检波器单点接收地震记录中的异常噪声方法 | |
CN103675899A (zh) | 一种基于子波压缩拓展叠后地震数据频带的方法 | |
CN103091714B (zh) | 一种自适应面波衰减方法 | |
CN111045077B (zh) | 一种陆地地震数据的全波形反演方法 | |
Jazayeri et al. | Sparse blind deconvolution of ground penetrating radar data | |
CN103984019B (zh) | 局部相关加权地震波束合成方法 | |
AU2016229330A1 (en) | Wideband channel equalization for signals propagated in lossy transmission media | |
CN105403918A (zh) | 一种三分量微地震数据有效事件识别方法及系统 | |
CN103630936A (zh) | 基于波束定向原理的地震单炮记录随机噪声压制方法 | |
CN102854532A (zh) | 三维叠前炮检域随机噪声压制方法 | |
CN103913770A (zh) | 基于vsp资料对地震数据进行处理的方法 | |
CN109270574B (zh) | 一种基于多种震源采集数据的联合反褶积方法 | |
CN104570121B (zh) | 定向地震波畸变信号消除方法 | |
CN104765069B (zh) | 一种压制同步激发采集邻炮干扰的方法 | |
CN110780341B (zh) | 一种各向异性地震成像方法 | |
EP3260887A1 (en) | Methods and data processing apparatus for seismic signal separation | |
CN105093327A (zh) | 多分量地震数据的矢量均值滤波方法 | |
Nguyen et al. | A comprehensive performance comparison of RFI mitigation techniques for UWB radar signals | |
CN105785322A (zh) | 一种波束形成方法 | |
CN104122583A (zh) | 一种拓宽地震数据频谱的方法和装置 | |
CN112327356A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170104 Termination date: 20211226 |