CN106569281A - 消除面波的椎体相干方法 - Google Patents
消除面波的椎体相干方法 Download PDFInfo
- Publication number
- CN106569281A CN106569281A CN201610983042.0A CN201610983042A CN106569281A CN 106569281 A CN106569281 A CN 106569281A CN 201610983042 A CN201610983042 A CN 201610983042A CN 106569281 A CN106569281 A CN 106569281A
- Authority
- CN
- China
- Prior art keywords
- data
- face
- ripple
- wave noise
- apparent velocity
- 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
- 238000000034 method Methods 0.000 title claims abstract description 23
- 230000001427 coherent effect Effects 0.000 claims abstract description 17
- 230000002452 interceptive effect Effects 0.000 claims abstract description 15
- 238000006073 displacement reaction Methods 0.000 claims abstract description 10
- 238000005070 sampling Methods 0.000 claims description 18
- 230000008030 elimination Effects 0.000 claims description 15
- 238000003379 elimination reaction Methods 0.000 claims description 15
- 238000013459 approach Methods 0.000 claims description 13
- 241001269238 Data Species 0.000 claims description 8
- 230000000694 effects Effects 0.000 abstract description 5
- 230000000007 visual effect Effects 0.000 abstract 4
- 230000001629 suppression Effects 0.000 abstract 2
- 238000001914 filtration Methods 0.000 description 4
- 235000012149 noodles Nutrition 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 238000005755 formation reaction Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic 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
本发明提出了一种消除面波的椎体相干方法，包括：采用交互速度拾取工具计算面波的视速度范围，对视速度范围进行射线分组，模拟面波在炮集记录上形成的多个视速度椎体；对每个视速度椎体进行处理得到当前面波的开始位置；对当前道面波开始样点所对应的相邻多道数据进行相关运算，得到每道数据与当前道数据最大相关时间位移，得到当前道最大线性相似的相干面波噪音；根据原始数据自相关、面波噪音自相关和原始数据与面波噪音互相关系数比值计算出道中每个样点处面波噪音压制系数；利用每个样点处面波噪音压制系数，在原始数据中减去面波噪音得到更新后的原始数据。本发明有效避免低频有效波数据受到损伤，去噪方法简单，去噪效果和振幅保真性好。
Description
技术领域
本发明涉及地震勘探领域，特别涉及一种消除面波的椎体相干方法。
背景技术
为了确定地下地层是否蕴含石油天然气资源，需要对地层进行地震勘探。目前主要采取反射波勘探技术，在油气勘探的区域布置二维和三维测线，使用炸药或者可控震源进行激发地震波，地震波经过地下岩层时发生反射传播到地面布置的检波器，同时也有部分地震波沿地表传播到检波器被记录下来，形成了面波。
面波是陆上地震勘探中广泛存在的一种规则强能量干扰波，沿地表传播，具有低频、低速、衰减较慢、频散的特点。目前在资料处理过程中采用滤波方法压制面波干扰，根据面波出现在地震数据上的时间和空间采用时间与空间约束的高通滤波方式进行频率域滤波。但是面波频率和有效信号频率分布存在交叉范围，这样滤波会使有效波成分受到一定程度的损伤。通过高通滤波方法能够在一定程度上压制面波，但是也有可能使有效信号受到损伤。
现代地震勘探处理要求在压制面波噪音的同时，要保护有效波的低频成分，即提高地震资料信噪比的同时，对资料的保真性提出了更高要求。
发明内容
本发明的目的旨在至少解决所述技术缺陷之一。
为此，本发明的目的在于提出一种消除面波的椎体相干方法，可以有效避免低频有效波数据受到损伤，去噪方法简单，去噪效果好，振幅保真性好。
为了实现上述目的，本发明的实施例提供一种消除面波的椎体相干方法，包括如下步骤：
步骤S1，对地震数据上的面波，采用交互速度拾取工具计算所述面波的视速度范围，对所述视速度范围进行射线分组，模拟面波在炮集记录上形成的多个视速度椎体，针对每个划分的视速度执行如下步骤S2至步骤S5；
步骤S2，对每个所述视速度椎体进行处理，得到视速度投影在道数据上的样点号，作为当前面波的开始位置；
步骤S3，对当前道面波开始样点所对应的相邻多道数据进行相关运算，得到每道数据与当前道数据最大相关时间位移，位移后数据进行叠加，得到当前道最大线性相似的相干面波噪音；
步骤S4，根据原始数据自相关、面波噪音自相关和原始数据与面波噪音互相关系数比值计算出道中每个样点处面波噪音压制系数；
步骤S5，利用所述每个样点处面波噪音压制系数，在原始数据中减去面波噪音得到更新后的原始数据。
进一步，在所述步骤S1中，所述采用交互速度拾取工具计算所述面波的视速度范围，包括：采用交互速度拾取工具沿面波干扰设置一条连线，由所述交互速度拾取工具自动计算当前面波对应的视速度。
进一步，在所述步骤S2中，当视速度投影在实际道有效数据之外时，按照有效数据的开始为起点进行计算，直到本道数据的结束样点。
进一步，在所述步骤S3中，所述对当前道面波开始样点所对应的相邻多道数据进行相关运算，包括：对当前道的左右各5道在预设位移范围内做相关运算，得到各道与本道数据最大近似的位移量。
进一步，在所述步骤S4中，所述原始数据指第一次刚进入处理的原始数据或者已经除去视速度面波噪音的更新后数据。
根据本发明实施例的消除面波的椎体相干方法，通过交互速度拾取工具确定面波噪音的视速度范围，在视速度椎体范围内按照速度增量进行分段，对相邻的多道数据作相干运算，统计出相干面波噪音系数并予以压制。本发明采用视速度射线椎体表示面波所在范围，获取压制面波区域准确度高，并且采用相干系数方法评估信噪比，能够准确检测出面波噪声并压制。本发明可以有效避免低频有效波数据受到损伤，去噪方法简单，去噪效果好，振幅保真性好。
本发明附加的方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
图1为根据本发明实施例的消除面波的椎体相干方法的流程图；
图2为根据本发明实施例的面波噪音视速度测量的示意图；
图3为根据本发明实施例的椎体相干技术去除面波噪音后的示意图。
具体实施方式
下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。
本发明针对地震勘探中反射地震资料处理，提出了一种消除面波的椎体相干方法，对面波进行噪音衰减。
如图1所示，本发明实施例的消除面波的椎体相干方法，包括如下步骤：
步骤S1，对地震数据上的面波，采用交互速度拾取工具计算面波的视速度范围，对视速度范围进行射线分组，模拟面波在炮集记录上形成的多个视速度椎体，针对每个划分的视速度执行如下步骤S2至步骤S5。
在本步骤中，由于面波沿地表传播，在地震数据上表现为线性低频干扰，根据线性干扰在相邻几道上的空间距离和时间差值，即可算出其大致的传播速度。采用交互速度拾取工具计算面波的视速度范围，包括：采用交互速度拾取工具沿面波干扰设置一条连线，由交互速度拾取工具自动计算当前面波对应的视速度。
为了测试椎体相干技术去除面波噪音效果，本步骤可以选取图2所示的炮道集地震数据，并进行视速度测量，测得的视速度为100-350米每秒。
步骤S2，对每个视速度椎体进行处理，得到视速度投影在道数据上的样点号，作为当前面波的开始位置。
在本发明的一个实施例中，当视速度投影在实际道有效数据之外时，按照有效数据的开始为起点进行计算，直到本道数据的结束样点。
步骤S3，对当前道面波开始样点所对应的相邻多道数据进行相关运算，得到每道数据与当前道数据最大相关时间位移，位移后数据进行叠加，得到当前道最大线性相似的相干面波噪音。
具体地，对当前道面波开始样点所对应的相邻多道数据进行相关运算，包括：对当前道的左右各5道在预设位移范围内做相关运算，得到各道与本道数据最大近似的位移量，然后进行线性叠加。由于这里是按照面波噪音线性方向进行的叠加，可以得到有效信号在此叠加道中不具有线性关系而被减弱。
步骤S4，根据原始数据自相关、面波噪音自相关和原始数据与面波噪音互相关系数比值计算出道中每个样点处面波噪音压制系数。在本步骤中，原始数据指第一次刚进入处理的原始数据或者已经除去视速度面波噪音的更新后数据。
步骤S5，利用每个样点处面波噪音压制系数，在原始数据中减去面波噪音得到更新后的原始数据。图3为对图2中地震数据，应用本发明的椎体相干技术去除面波噪音后的示意图。
根据本发明实施例的消除面波的椎体相干方法，通过交互速度拾取工具确定面波噪音的视速度范围，在视速度椎体范围内按照速度增量进行分段，对相邻的多道数据作相干运算，统计出相干面波噪音系数并予以压制。本发明采用视速度射线椎体表示面波所在范围，获取压制面波区域准确度高，并且采用相干系数方法评估信噪比，能够准确检测出面波噪声并压制。本发明可以有效避免低频有效波数据受到损伤，去噪方法简单，去噪效果好，振幅保真性好。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。本发明的范围由所附权利要求及其等同限定。
Claims (5)
1.一种消除面波的椎体相干方法，其特征在于，包括如下步骤：
步骤S1，对地震数据上的面波，采用交互速度拾取工具计算所述面波的视速度范围，对所述视速度范围进行射线分组，模拟面波在炮集记录上形成的多个视速度椎体，针对每个划分的视速度执行如下步骤S2至步骤S5；
步骤S2，对每个所述视速度椎体进行处理，得到视速度投影在道数据上的样点号，作为当前面波的开始位置；
步骤S3，对当前道面波开始样点所对应的相邻多道数据进行相关运算，得到每道数据与当前道数据最大相关时间位移，位移后数据进行叠加，得到当前道最大线性相似的相干面波噪音；
步骤S4，根据原始数据自相关、面波噪音自相关和原始数据与面波噪音互相关系数比值计算出道中每个样点处面波噪音压制系数；
步骤S5，利用所述每个样点处面波噪音压制系数，在原始数据中减去面波噪音得到更新后的原始数据。
2.如权利要求1所述的消除面波的椎体相干方法，其特征在于，在所述步骤S1中，所述采用交互速度拾取工具计算所述面波的视速度范围，包括：采用交互速度拾取工具沿面波干扰设置一条连线，由所述交互速度拾取工具自动计算当前面波对应的视速度。
3.如权利要求1所述的消除面波的椎体相干方法，其特征在于，在所述步骤S2中，当视速度投影在实际道有效数据之外时，按照有效数据的开始为起点进行计算，直到本道数据的结束样点。
4.如权利要求1所述的消除面波的椎体相干方法，其特征在于，在所述步骤S3中，所述对当前道面波开始样点所对应的相邻多道数据进行相关运算，包括：对当前道的左右各5道在预设位移范围内做相关运算，得到各道与本道数据最大近似的位移量。
5.如权利要求1所述的消除面波的椎体相干方法，其特征在于，在所述步骤S4中，所述原始数据指第一次刚进入处理的原始数据或者已经除去视速度面波噪音的更新后数据。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610983042.0A CN106569281B (zh) | 2016-11-09 | 2016-11-09 | 消除面波的锥体相干方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610983042.0A CN106569281B (zh) | 2016-11-09 | 2016-11-09 | 消除面波的锥体相干方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106569281A true CN106569281A (zh) | 2017-04-19 |
| CN106569281B CN106569281B (zh) | 2019-01-04 |
Family
ID=58540619
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610983042.0A Active CN106569281B (zh) | 2016-11-09 | 2016-11-09 | 消除面波的锥体相干方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106569281B (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2011144623A1 (en) * | 2010-05-19 | 2011-11-24 | Cggveritas Services Sa | Passive monitoring method for seismic events |
| CN103091714A (zh) * | 2011-10-28 | 2013-05-08 | 中国石油化工股份有限公司 | 一种自适应面波衰减方法 |
| US20140095079A1 (en) * | 2012-10-01 | 2014-04-03 | Cgg Services Sa | System and method for three-dimensional iterative filtering of scattered waves in cross spread seismic systems |
| CN105093282A (zh) * | 2014-05-16 | 2015-11-25 | 中国石油化工股份有限公司 | 基于频率约束的能量置换面波压制方法 |
-
2016
- 2016-11-09 CN CN201610983042.0A patent/CN106569281B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2011144623A1 (en) * | 2010-05-19 | 2011-11-24 | Cggveritas Services Sa | Passive monitoring method for seismic events |
| CN103091714A (zh) * | 2011-10-28 | 2013-05-08 | 中国石油化工股份有限公司 | 一种自适应面波衰减方法 |
| US20140095079A1 (en) * | 2012-10-01 | 2014-04-03 | Cgg Services Sa | System and method for three-dimensional iterative filtering of scattered waves in cross spread seismic systems |
| CN105093282A (zh) * | 2014-05-16 | 2015-11-25 | 中国石油化工股份有限公司 | 基于频率约束的能量置换面波压制方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106569281B (zh) | 2019-01-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108345031B (zh) | 弹性介质主动源和被动源混采地震数据全波形反演方法 | |
| CN103913776B (zh) | 盾构施工中球状孤石的探测方法 | |
| CN101334483B (zh) | 一种在地震数据处理中衰减瑞雷波散射噪声的方法 | |
| CN101551465B (zh) | 一种自适应识别和消除地震勘探单频干扰的方法 | |
| CN109669212B (zh) | 地震数据处理方法、地层品质因子估算方法与装置 | |
| CN106526677A (zh) | 一种海上自适应压制鬼波的宽频逆时偏移成像方法 | |
| CN104360388B (zh) | 一种三维地震观测系统评价方法 | |
| CN103424777B (zh) | 一种提高地震成像分辨率的方法 | |
| CN107065013B (zh) | 一种地震尺度下的层速度确定方法及装置 | |
| CN102099705A (zh) | 处理地震数据以供解读的系统和方法 | |
| CN109188520B (zh) | 薄储层厚度预测方法及装置 | |
| US10598807B2 (en) | Correction of sea surface state | |
| CN104330826A (zh) | 一种去除复杂地表条件下多种噪音的方法 | |
| CN105158801B (zh) | 光缆耦合噪声的压制处理方法和装置 | |
| CN101021568A (zh) | 基于最大能量旅行时计算的三维积分叠前深度偏移方法 | |
| CN109765615A (zh) | 一种地层品质因子反演方法及装置 | |
| CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
| CN104199088B (zh) | 一种提取入射角道集的方法及系统 | |
| CN106199699B (zh) | 利用传输矩阵法清除鬼波的方法 | |
| WO2021155754A1 (zh) | 光纤声波传感地震数据中井筒波干扰的去除方法及装置 | |
| CN112230274B (zh) | 面向随钻导向的声波方程频率域逆时偏移快速成像方法 | |
| CN109975873B (zh) | 一种逆时偏移成像去除低频噪音的方法及系统 | |
| CN106569281B (zh) | 消除面波的锥体相干方法 | |
| CN105093275B (zh) | 基于速度模型的有效信号提取iss波场分离方法 | |
| CN112198547A (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 |