CN102854532A - 三维叠前炮检域随机噪声压制方法 - Google Patents
三维叠前炮检域随机噪声压制方法 Download PDFInfo
- Publication number
- CN102854532A CN102854532A CN2011101805108A CN201110180510A CN102854532A CN 102854532 A CN102854532 A CN 102854532A CN 2011101805108 A CN2011101805108 A CN 2011101805108A CN 201110180510 A CN201110180510 A CN 201110180510A CN 102854532 A CN102854532 A CN 102854532A
- Authority
- CN
- China
- Prior art keywords
- data
- dimensional
- big gun
- territory
- rearrangement
- 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
本发明是地球物理勘探地震数据处理三维叠前炮检域随机噪声压制方法,将动校正之后的数据在线方向进行线号重排,在沿测线方向按炮序重排,使重排之后每条CMP线相当于一条二维测线,修改道头数据,顺序重排,进行三维叠前炮检域随机噪声处理,得到全部模型道,将全部模型道与动校正之后数据混波,完成噪声压制。本方法通过对观测系统的变形重排技术,可以在不同域中实现对随机噪声压制,去噪效果更加显著。
Description
技术领域
本发明涉及地球物理勘探地震数据处理方法,属于压制噪声的技术范畴,具体一种三维叠前炮检域随机噪声压制方法。
背景技术
在陆上地震数据采集中,由于地表复杂及地下地质结构等原因,采集的地震数据信噪比极低,特别是沙漠地区,随机噪声比较严重,常规做法是在CMP(沿测线方向)域中处理,这种方法很难满足需求。
在CMP域中的实现方法是首先将数据抽道CMP域中,一个CMP内的道数比较少,受实际覆盖次数的影响,预测方法使用的空间道数就少,预测效果就受影响,另外在CMP域预测有一个缺陷,使信噪比高的地方改进很大,但信噪比低的地方改进不大,而且后续再进行其它域中去噪效果不好。
发明内容
本发明目的是提供一种通过对观测系统的变形重排技术,去噪效果更加显著的三维叠前炮检域随机噪声压制方法。
本发明通过以下步骤实现:
1)采集地震数据,输入观测系统定义以后的数据,经过处理得到动校正之后的CMP道集数据,并保留;
步骤1)所述的处理是做常规叠加处理;所述的保留是将对动校正之后的数据记录保留,准备做三维叠前炮检域随机噪声压制。
2)将动校正之后的数据在线方向进行线号重排,在沿测线方向(CMP)按炮序重排,使重排之后每条CMP线相当于一条二维测线;
3)修改道头数据,顺序重排;
步骤3)所述的数据重排是在三维施工时,一炮一般对应多个接收排列,每个接收排列与每个平行炮线相当于一条二维线,把每个二维线按地表连续位置重排,每个排列对应一炮。
4)进行三维叠前炮检域随机噪声处理,得到全部模型道;
步骤4)所述的三维叠前炮检域随机噪声压制是首先对四维数据体按照空间上的参数为线数,炮数,道数,时间上给定一个时窗,这样就形成一个小的三维数据体,把这个三维数据体变换到相应的F-XYZ域,在F-XYZ域求取预测算子,然后对这个小数据体进行预测滤波,得到预测结果;
然后在空间上重复道数,在时间上重复样点数,得到下一个小的三维数据体,然后求取下一个三维数据体的预测结果,反复在空间和时间上滑动求取,完成对整个三维叠前数据的处理。
5)将全部模型道与步骤1)数据混波,完成噪声压制。
步骤5)所述的混波是首先将步骤4)全部模型道做叠加,与步骤1)步保留的常规叠加做对比,根据两者剖面噪比情况,确定混波比例。
所述的确定混波比例是信噪比低的部位,混入模型道比例加大。
本方法通过对观测系统的变形重排技术,可以在不同域中实现对随机噪声压制,去噪效果更加显著。
附图说明
图1是数据重排前的炮检位置图;
图2是数据重排后的炮检位置图;
图3是随机噪声压制的时窗滑动示意图;
图4是随机噪声压制前道集;
图5是随机噪声压制后道集;
图6是随机噪声压制前叠加;
图7是随机噪声压制后叠加。
具体实施方式
本发明是随机噪声衰减方法,用于三维地震采集的数据。本发明基于三维频率空间的F-XYZ域预测去噪技术。它假设地震记录中的有效波在F-XYZ域具有可预测性,而随机噪声无此特性。利用多道复数最小平方原理求取三维预测算子,并用该预测算子对该频率成分的四维地震数据体进行预测滤波,达到衰减随机噪声的目的。
本发明具体步骤如下:
第一步:采集地震数据,输入观测系统定义之后的数据,经过一系列处理,得到动校正之后的CMP道集数据。一系列处理包括:静校正,振幅保持处理,规则噪声压制处理,反褶积(可选),动校正,剩余静校正等常规处理。最后得到动校正之后的CMP道集数据。
第二步:将动校正之后的数据在线方向进行线号重排,在沿测线方向(CMP)按炮序重排,使重排之后每条CMP线相当于一条二维测线;见附图1和2,重排前后的炮检位置图。从图中看到,数据重排之后,位置图更加规则。如果按三维显示的化,数据的最后一维是一个排列的炮集,所以叫炮检域RNA技术。这一步也叫数据准备。
第三步:进行三维叠前炮检域随机噪声压制处理,得到全部模型道,这一步也叫三维叠前炮检域随机噪声压制。
三维叠前炮检域随机噪声压制是首先对四维数据体按照空间上的参数为线数,炮数,道数,时间上给定一个时窗,这样就形成一个小的三维数据体,把这个三维数据体变换到相应的F-XYZ域,在F-XYZ域求取预测算子,然后对这个小数据体进行预测滤波,得到预测结果;
然后在空间上重复道数,在时间上重复样点数,得到下一个小的三维数据体,然后求取下一个三维数据体的预测结果,反复在空间和时间上滑动求取,见图3,完成对整个三维叠前数据的处理。这步的输出通常叫模型道。
第四步将全部模型道与第一步数据叠前混波,完成噪声压制。
这里所述的混波是首先将第三步全部模型道做叠加,与第一步保留的常规叠加做对比,根据两者剖面噪比情况,确定混波比例,然后用到叠前数据。
所述的确定混波比例是信噪比低的部位,混入模型道比例加大。
本方法通过对观测系统的变形重排技术,可以在不同域中实现对随机噪声压制,去噪效果更加显著。图4与图5噪声压制前后的道集数据,图6与图7是噪声压制前后的叠加数据,效果明显。
Claims (5)
1.一种三维叠前炮检域随机噪声压制方法,特点是通过以下步骤实现:
1)地震采集,将观测系统定义后数据处理,得到动校正之后的CMP道集数据,并保留;
2)将动校正之后的数据在线方向进行线号重排,在沿测线方向CMP按炮序重排,使重排之后每条CMP线相当于一条二维测线;
3)修改道头数据,顺序重排;
4)进行三维叠前炮检域随机噪声处理,得到全部模型道;
5)将全部模型道与步骤1)数据混波,完成噪声压制。
2.根据权利要求1所述的方法,特点是步骤1)所述的处理是做常规叠加处理;所述的保留是将对动校正之后的数据记录保留,准备做三维叠前炮检域随机噪声压制。
3.根据权利要求1所述的方法,特点是步骤3)所述的数据重排是在三维施工时,一炮一般对应多个接收排列,每个接收排列与每个平行炮线相当于一条二维线,把每个二维线按地表连续位置重排,每个排列对应一炮。
4.根据权利要求1所述的方法,特点是步骤4)所述的三维叠前炮检域随机噪声压制是首先对四维数据体按照空间上的参数为线数,炮数,道数,时间上给定一个时窗,这样就形成一个小的三维数据体,把这个三维数据体变换到相应的F-XYZ域,在F-XYZ域求取预测算子,然后对这个小数据体进行预测滤波,得到预测结果;
然后在空间上重复道数,在时间上重复样点数,得到下一个小的三维数据体,然后求取下一个三维数据体的预测结果,反复在空间和时间上滑动求取,完成对整个三维叠前数据的处理。
5.根据权利要求1所述的方法,特点是步骤5)所述的混波是首先将步骤4)全部模型道做叠加,与步骤1)步保留的常规叠加做对比,根据两者剖面噪比情况,确定混波比例。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110180510.8A CN102854532B (zh) | 2011-06-30 | 2011-06-30 | 三维叠前炮检域随机噪声压制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110180510.8A CN102854532B (zh) | 2011-06-30 | 2011-06-30 | 三维叠前炮检域随机噪声压制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102854532A true CN102854532A (zh) | 2013-01-02 |
CN102854532B CN102854532B (zh) | 2015-07-01 |
Family
ID=47401278
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110180510.8A Active CN102854532B (zh) | 2011-06-30 | 2011-06-30 | 三维叠前炮检域随机噪声压制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102854532B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104020492A (zh) * | 2013-07-01 | 2014-09-03 | 西安交通大学 | 一种三维地震资料的保边滤波方法 |
CN106338769A (zh) * | 2015-07-07 | 2017-01-18 | 中国石油化工股份有限公司 | 地震数据去噪方法及系统 |
CN108646296A (zh) * | 2018-05-16 | 2018-10-12 | 吉林大学 | 基于自适应谱峭度滤波器的沙漠地震信号噪声消减方法 |
CN109031410A (zh) * | 2018-07-16 | 2018-12-18 | 中国石油大学(华东) | 共偏移距域结果约束的炮检域双边束合成方法及系统 |
CN112782766A (zh) * | 2019-11-11 | 2021-05-11 | 中国石油天然气股份有限公司 | 去除地震资料侧源干扰的方法及装置 |
CN112799132A (zh) * | 2019-11-13 | 2021-05-14 | 中国石油天然气股份有限公司 | 微局部线性噪声压制方法及装置 |
CN112882101A (zh) * | 2019-11-29 | 2021-06-01 | 中国石油天然气集团有限公司 | 叠前地震数据的随机噪音衰减方法及装置 |
CN113126162A (zh) * | 2019-12-30 | 2021-07-16 | 中国石油天然气集团有限公司 | 随机噪声衰减计算方法及装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1797040A (zh) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | 压制低信噪比地震记录中随机噪声的方法 |
CN101419293A (zh) * | 2007-10-25 | 2009-04-29 | 中国石油天然气集团公司 | 一种提高地震数据信噪比的方法 |
-
2011
- 2011-06-30 CN CN201110180510.8A patent/CN102854532B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1797040A (zh) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | 压制低信噪比地震记录中随机噪声的方法 |
CN101419293A (zh) * | 2007-10-25 | 2009-04-29 | 中国石油天然气集团公司 | 一种提高地震数据信噪比的方法 |
Non-Patent Citations (1)
Title |
---|
苏贵士等: "频率空间(三维)F-XYZ域预测去噪技术", 《石油地球物理勘探》, vol. 33, no. 1, 28 February 1998 (1998-02-28), pages 95 - 103 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104020492A (zh) * | 2013-07-01 | 2014-09-03 | 西安交通大学 | 一种三维地震资料的保边滤波方法 |
CN104020492B (zh) * | 2013-07-01 | 2015-10-28 | 西安交通大学 | 一种三维地震资料的保边滤波方法 |
CN106338769A (zh) * | 2015-07-07 | 2017-01-18 | 中国石油化工股份有限公司 | 地震数据去噪方法及系统 |
CN108646296A (zh) * | 2018-05-16 | 2018-10-12 | 吉林大学 | 基于自适应谱峭度滤波器的沙漠地震信号噪声消减方法 |
CN109031410A (zh) * | 2018-07-16 | 2018-12-18 | 中国石油大学(华东) | 共偏移距域结果约束的炮检域双边束合成方法及系统 |
CN109031410B (zh) * | 2018-07-16 | 2019-07-02 | 中国石油大学(华东) | 共偏移距域结果约束的炮检域双边束合成方法及系统 |
CN112782766A (zh) * | 2019-11-11 | 2021-05-11 | 中国石油天然气股份有限公司 | 去除地震资料侧源干扰的方法及装置 |
CN112799132A (zh) * | 2019-11-13 | 2021-05-14 | 中国石油天然气股份有限公司 | 微局部线性噪声压制方法及装置 |
CN112799132B (zh) * | 2019-11-13 | 2023-08-22 | 中国石油天然气股份有限公司 | 微局部线性噪声压制方法及装置 |
CN112882101A (zh) * | 2019-11-29 | 2021-06-01 | 中国石油天然气集团有限公司 | 叠前地震数据的随机噪音衰减方法及装置 |
CN112882101B (zh) * | 2019-11-29 | 2024-04-30 | 中国石油天然气集团有限公司 | 叠前地震数据的随机噪音衰减方法及装置 |
CN113126162A (zh) * | 2019-12-30 | 2021-07-16 | 中国石油天然气集团有限公司 | 随机噪声衰减计算方法及装置 |
CN113126162B (zh) * | 2019-12-30 | 2024-05-28 | 中国石油天然气集团有限公司 | 随机噪声衰减计算方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN102854532B (zh) | 2015-07-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102854532B (zh) | 三维叠前炮检域随机噪声压制方法 | |
EP2707757B1 (en) | Couple time-distance dependent swept frequency source acquisition design and data de-noising | |
CN101487899B (zh) | 三维双传感器拖缆数据中的波场分离方法 | |
US12007515B2 (en) | Optimal survey design | |
CN104062683A (zh) | 一种基于曲波变换和全变差的联合衰减随机噪声处理方法 | |
CN103954992B (zh) | 一种反褶积方法及装置 | |
US10634802B2 (en) | Prime number survey design | |
CN103364826A (zh) | 基于独立分量分析的地震盲源反褶积方法 | |
CN105467460A (zh) | 电磁勘探方法和装置 | |
CN103901473B (zh) | 一种基于非高斯性最大化的双检信号上下行波场分离方法 | |
CN114185095A (zh) | 一种三维平面波域地震数据多次波压制的方法 | |
Shabrawi et al. | How single-sensor seismic improved image of Kuwait’s Minagish Field | |
Retailleau et al. | Advanced 3D Land Internal Mutiple Modeling and Subtraction, a WAZ Oman Case Study | |
El-Emam et al. | Maximizing the value of Single-Sensor measurements, Kuwait experience | |
Al-Shuhail et al. | Processing of Seismic Reflection Data Using MATLAB | |
CN104880732B (zh) | 一种十字子集的构建方法及装置 | |
CN113835120B (zh) | 不同激发能量大小的正演模拟分析方法 | |
CN113391351B (zh) | 一种基于被动源地震波场分析提取矿集区结构的方法 | |
CN112711067B (zh) | 薄储层预测方法及装置 | |
CN104990774B (zh) | 一种地震数据插值方法 | |
Urosevic et al. | Acquisition of the Nested 3D seismic survey at Harvey | |
Wu et al. | Least-squares reverse time migration of simultaneous sources with deep-learning-based denoising | |
CN103760599A (zh) | 一种微小断层检测方法及断层检测装置 | |
CN117724149A (zh) | 一种五维数据规则化处理方法及装置 | |
Mojica et al. | Using SVD filters for velocity analysis and ground-roll attenuation |
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 |