CN110895349B - 基于中值滤波处理需求的独立扫描信号的设计方法 - Google Patents
基于中值滤波处理需求的独立扫描信号的设计方法 Download PDFInfo
- Publication number
- CN110895349B CN110895349B CN201811062296.4A CN201811062296A CN110895349B CN 110895349 B CN110895349 B CN 110895349B CN 201811062296 A CN201811062296 A CN 201811062296A CN 110895349 B CN110895349 B CN 110895349B
- Authority
- CN
- China
- Prior art keywords
- aliasing
- design method
- length
- scanning
- scan
- 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
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000001914 filtration Methods 0.000 title claims abstract description 17
- 230000001629 suppression Effects 0.000 claims abstract description 13
- 238000001514 detection method Methods 0.000 claims abstract description 8
- 230000002452 interceptive effect Effects 0.000 claims description 4
- 238000003384 imaging method Methods 0.000 abstract description 3
- 238000005516 engineering process Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Images
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
- G01V1/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
-
- 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
- G01V2210/324—Filtering
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
本发明提出了一种基于中值滤波处理需求的独立扫描信号的设计方法，其中，所述设计方法包括以下步骤：S101：在实施可控震源独立扫描区域内建立二维地质模型，同时激发变化扫描长度的干扰炮和固定扫描长度的有效炮获得多个混叠波场；S102：采用固定扫描信号分别与所获得的多个混叠波场进行相关以获取多个混叠单炮；S103：将所述多个混叠单炮分别抽成共检波点道集；S104：采用中值滤波分别对共检波点道集进行噪音压制；S105：根据所述噪音压制的效果，选择合理的独立扫描信号。本发明的设计方法可以实现提高采集效率的同时，降低混叠波场对目的层成像的影响。
Description
技术领域
本发明涉及油气地球物理勘探领域，特别涉及一种基于中值滤波处理需求的独立扫描信号的方法。
背景技术
可控震源技术是目前西部地震勘探采集的重要技术，应用广泛，独立扫描技术的核心及难点是如何分离波场，这分为在采集过程中如何设计独立扫描信号，采用相关性较差的独立扫描信号，以便有利于后期的波场分离。而现阶段还缺少这方面的技术。
发明内容
为解决现有技术的缺陷，本发明的目的在于针对以上问题展开的研究，采用不同扫描长度信号进行混叠波场的模拟，在获取不同影响大小的混叠波场后采用中值滤波方式分别进行压制分析，从压制效果可以直观的分析判断适合于独立扫描信号，提高采集效率的同时，降低混叠波场对目的层成像的影响。为实现该目的，本发明提出了一种基于中值滤波处理需求的独立扫描信号的设计方法，其中，所述设计方法包括以下步骤：
S101：在实施可控震源独立扫描区域内建立二维地质模型，同时激发变化扫描长度的干扰炮和固定扫描长度的有效炮获得多个混叠波场；
S102：采用固定扫描信号分别与所获得的多个混叠波场进行相关以获取多个混叠单炮；
S103：将所述多个混叠单炮分别抽成共检波点道集；
S104：采用中值滤波分别对共检波点道集进行噪音压制；
S105：根据所述噪音压制的效果，选择合理的独立扫描信号。
如上所述的设计方法，其中，在所述步骤S101中，所述有效炮的扫描长度为20s。
如上所述的设计方法，其中，在所述步骤S101中，所述干扰炮的扫描长度逐渐依次改变。
如上所述的设计方法，其中，所述干扰炮的扫描长度从10s逐渐增加到20s，其中，增加的步长为1s。
如上所述的设计方法，其中，在所述步骤S102中，所述固定扫描信号的扫描长度为20s。
如上所述的设计方法，其中，在所述步骤S103中，以所述有效波的sps为主，将所述多个混叠单炮分别抽成共检波点道集。
如上所述的设计方法，其中，在所述共检波点道集中，所述有效炮的波场是连续的，所述干扰炮的同相轴被打散，所述有效炮的扫描长度与所述干扰炮的扫描长度不同，所获得的随机干扰能量大小不同。
如上所述的设计方法，其中，选择所述有效炮的扫描长度与所述干扰炮的扫描长度之间的差值为4s。
本发明提供一种基于中值滤波处理需求的独立扫描信号设计方法。该设计方法包括建立二维地质模型，采用不同扫描长度信号进行混叠波场的模拟，在获取不同影响大小的混叠波场后采用中值滤波方式分别进行压制分析，从压制效果可以直观的分析判断适合于独立扫描信号，提高采集效率的同时，降低混叠波场对目的层成像的影响。
附图说明
在此描述的附图仅用于解释目的，而不意图以任何方式来限制本发明公开的范围。另外，图中的各部件的形状和比例尺寸等仅为示意性的，用于帮助对本发明的理解，并不是具体限定本发明各部件的形状和比例尺寸。本领域的技术人员在本发明的教导下，可以根据具体情况选择各种可能的形状和比例尺寸来实施本发明。
图1为本发明的基于中值滤波处理需求的独立扫描信号的设计方法的流程图；
图2为本发明的不同扫描长度压制前后对比图。
具体实施方式
结合附图和本发明具体实施方式的描述，能够更加清楚地了解本发明的细节。但是，在此描述的本发明的具体实施方式，仅用于解释本发明的目的，而不能以任何方式理解成是对本发明的限制。在本发明的教导下，技术人员可以构想基于本发明的任意可能的变形，这些都应被视为属于本发明的范围，下面将结合附图对本发明作进一步说明。
为实现本发明的上述目的，本发明的设计方法可通过如下技术措施来实现：
步骤1，根据所要实施可控震源独立扫描区域的地质构造，建立二维地质模型，在该模型上进行可控震源扫描信号模拟，有效炮固定扫描长度，而干扰炮扫描长度逐渐依次改变；
步骤2，在获取混叠波场之后，采用固定扫描信号分别与混叠波场进行相关获取混叠单炮；
步骤3，将混叠单炮抽成共检波点道集；
步骤4，采用中值滤波分别对这些共检波点道集进行噪音压制；
步骤5，从噪音压制效果选择合理的独立扫描信号。
为使本发明的目的、特征和优点能更明显易懂，结合图1和图2，现进行如下详细说明：
如图1所示，图1为一种本发明的基于中值滤波处理需求的独立扫描信号的设计方法的流程图。
在步骤101，建立二维地质模型，在该模型上进行可控震源扫描信号模拟，要求同时激发两炮，一炮作为有效炮，另外一炮作为干扰炮，有效炮固定扫描长度，而干扰炮扫描长度逐渐依次改变，在一具体实施例中，有效炮扫描长度为20s，干扰炮扫描长度从10s逐渐增大到20s，步长为1s，总共获得11个模拟记录，每一个模拟记录中都混叠存在2炮母记录。流程进入到步骤102。
在一实施例中，步骤102为，将这11个模拟记录分别于20s扫描信号进行相关，获得11个混叠波场波场的单炮，有效炮与干扰炮混叠在一起。进入到步骤103。
在步骤103，将这些混叠单炮，以有效炮的sps为主，抽成共检波点道集，在共检波点道集中，有效炮仍然是波场连续记录，而干扰炮同相轴则被打散，表现为随机干扰的特征，不同扫描长度所获得的随机干扰能量大小不同。流程进入到步骤104。
在步骤104，采用中值滤波方式压制这些随机干扰，不同能量大小的随机干扰被压制的效果不同。流程进入到步骤105。
在步骤105，从这些压制后的效果分析采用多大差异的扫描长度信号进行独立扫描合适，更有利于独立扫描后混叠波场的分离。流程结束。
如图2所示，图2是通过不同扫描长度信号作为干扰炮，从20s逐渐变化到10s，采用固定的20s扫描长度信号作为有效波，通过正演模拟后与20s扫描信号进行相关，获取单炮后抽检波点道集，干扰炮表现为随机干扰，采用中值滤波对这些混叠波场单炮进行压制，从压制效果来看，采用扫描长度差异为4s的信号可以满足处理的需要。
Claims (7)
1.基于中值滤波处理需求的独立扫描信号的设计方法，其特征在于，所述设计方法包括以下步骤：
S101：在实施可控震源独立扫描区域内建立二维地质模型，同时激发变化扫描长度的干扰炮和固定扫描长度的有效炮获得多个混叠波场；
S102：采用固定扫描信号分别与所获得的多个混叠波场进行相关以获取多个混叠单炮；
S103：将所述多个混叠单炮分别抽成共检波点道集；
S104：采用中值滤波分别对共检波点道集进行噪音压制；
S105：根据所述噪音压制的效果，选择合理的独立扫描信号；
其中，在所述步骤S101中，所述干扰炮的扫描长度以设定的步长逐渐依次增加。
2.根据权利要求1所述的设计方法，其特征在于，
在所述步骤S101中，所述有效炮的扫描长度为20s。
3.根据权利要求1所述的设计方法，其特征在于，
所述干扰炮的扫描长度从10s逐渐增加到20s，其中，增加的步长为1s。
4.根据权利要求1所述的设计方法，其特征在于，
在所述步骤S102中，所述固定扫描信号的扫描长度为20s。
5.根据权利要求1至4中任意一项所述的设计方法，其特征在于，
在所述步骤S103中，以所述有效波的sps为主，将所述多个混叠单炮分别抽成共检波点道集。
6.根据权利要求5所述的设计方法，其特征在于，
在所述共检波点道集中，所述有效炮的波场是连续的，所述干扰炮的同相轴被打散，所述有效炮的扫描长度与所述干扰炮的扫描长度不同，所获得的随机干扰能量大小不同。
7.根据权利要求3所述的设计方法，其特征在于，选择所述有效炮的扫描长度与所述干扰炮的扫描长度之间的差值为4s的信号作为独立扫描信号。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811062296.4A CN110895349B (zh) | 2018-09-12 | 2018-09-12 | 基于中值滤波处理需求的独立扫描信号的设计方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811062296.4A CN110895349B (zh) | 2018-09-12 | 2018-09-12 | 基于中值滤波处理需求的独立扫描信号的设计方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110895349A CN110895349A (zh) | 2020-03-20 |
| CN110895349B true CN110895349B (zh) | 2021-11-05 |
Family
ID=69784879
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811062296.4A Active CN110895349B (zh) | 2018-09-12 | 2018-09-12 | 基于中值滤波处理需求的独立扫描信号的设计方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110895349B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112327356A (zh) * | 2020-10-30 | 2021-02-05 | 中国海洋大学 | 基于同相轴迭代追踪提取的混叠记录分离方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7881160B2 (en) * | 2008-04-04 | 2011-02-01 | Ion Geophysical Corporation | Seismic vibrator array and methods of operation |
| US9304216B2 (en) * | 2009-02-05 | 2016-04-05 | Westerngeco L.L.C. | Seismic acquisition system and technique |
| US10371839B2 (en) * | 2013-06-05 | 2019-08-06 | Saudi Arabian Oil Company | Blended land seismic data acquisition employing dispersed source arrays with variable sweep length |
| CN104216010B (zh) * | 2013-06-05 | 2017-03-01 | 上海青凤致远地球物理地质勘探科技有限公司 | 利用可控震源谐波提高地震数据质量的方法 |
| CN108181648B (zh) * | 2017-11-23 | 2019-07-09 | 中国石油天然气集团公司 | 连续记录混叠地震采集数据的质控方法和装置 |
| CN108205155B (zh) * | 2017-12-27 | 2019-11-08 | 中国石油天然气集团公司 | 一种可控震源交替独立同步激发方法、装置及系统 |
-
2018
- 2018-09-12 CN CN201811062296.4A patent/CN110895349B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN110895349A (zh) | 2020-03-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3132288B1 (en) | Near-continuous time-based marine seismic data acquisition and processing | |
| CN110907995B (zh) | 井中vsp地震数据的逆时偏移方法及装置 | |
| CN111158049A (zh) | 一种基于散射积分法的地震逆时偏移成像方法 | |
| CA2859286A1 (en) | Method and device for separating seismic signals from seismic sources | |
| CN105242309A (zh) | 一种地震数据中规则干扰的压制方法及装置 | |
| CN102798894A (zh) | 一种基于稀疏反演的滑动扫描谐波压制方法 | |
| CN110895349B (zh) | 基于中值滤波处理需求的独立扫描信号的设计方法 | |
| US10788596B2 (en) | Method and system for generating geophysical data | |
| WO2015001058A2 (en) | Method and device for de-blending seismic data using source signature | |
| GB2530126A (en) | Seismic Data Processing | |
| CN102103215B (zh) | 三维高密度地震勘探记录叠前面波压制方法 | |
| CN113534259B (zh) | 一种可控震源高效采集实时叠前时间偏移成像方法 | |
| CN104181600A (zh) | 地震数据线性噪声衰减方法及装置 | |
| CN105277983A (zh) | 新型压制可控震源滑动扫描谐波干扰的方法 | |
| CN105093282A (zh) | 基于频率约束的能量置换面波压制方法 | |
| CN114791623A (zh) | 一种微动采集方法 | |
| US11422277B2 (en) | Seismic data filtering based on distances between seismic sources | |
| CN110673210B (zh) | 一种地震原始数据信噪比定量分析评价方法 | |
| CN108363739B (zh) | 一种基于稀疏采集的地震资料高低频拓展方法 | |
| CN109270576B (zh) | 多通道三维地震物理模拟数据采集方法及装置 | |
| Yang et al. | Improved FK migration based on interpolation method for GPR imaging | |
| CN111665563A (zh) | 基于聚焦分析的叠前偏移垂向分辨率评价方法 | |
| CN104865601B (zh) | 确定地震勘探中采集的地震记录的类型的方法 | |
| CN113721297B (zh) | 利用速度特征压制沙丘鸣震的方法 | |
| CN103454675A (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 |