CN105093326A - 一种地震勘探资料中波长剩余静校正处理方法 - Google Patents
一种地震勘探资料中波长剩余静校正处理方法 Download PDFInfo
- Publication number
- CN105093326A CN105093326A CN201410221656.6A CN201410221656A CN105093326A CN 105093326 A CN105093326 A CN 105093326A CN 201410221656 A CN201410221656 A CN 201410221656A CN 105093326 A CN105093326 A CN 105093326A
- Authority
- CN
- China
- Prior art keywords
- static correction
- medium wavelength
- gather data
- shot
- data
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种地震勘探资料中波长剩余静校正处理方法,属于地球物理勘探地震资料处理领域。本方法包括:(1)输入构造形态控制的高程静校正数据和叠加速度,然后抽取共检波点道集数据、共炮点道集数据;(2)速度扫描获取线性动校正的速度;(3)优选参与叠加的炮检距范围及叠加;(4)检波点和炮点中波长静校正量提取及应用;(5)多次迭代求取准确的地表一致性剩余静校正量并应用。
Description
技术领域
本发明属于地球物理勘探地震资料处理领域,具体涉及一种地震勘探资料中波长剩余静校正处理方法,可应用于黄土塬地区地震资料静校正处理以及类似条件下复杂静校正的地震资料处理。
背景技术
由于黄土塬地区受表层结构特征(高程、风化层的速度、厚度和潜水面埋深)在区域上的不均匀性等诸多因素的影响,除了存在严重的短波长静校正问题外,还存在一定的中、长波长静校正问题,静校正问题严重且极其复杂是黄土塬地区资料的主要特点,它严重影响到均方根速度的准确求取和构造的落实。如何确保静校正的真实可靠应是黄土塬地区资料处理中的最大难点。当前黄土塬地震资料处理中,静校正处理技术还不是很成熟,静校正处理问题仍然有待研究解决。在采用基于标志层构造形态控制的高程静校正方法处理后,还存在较大的中短波长的静校正问题。
发明内容
本发明的目的在于解决上述现有技术中存在的难题,提供一种地震勘探资料中波长剩余静校正处理方法,基于共炮点、共检波点道集数据叠加提取中波长剩余静校正处理方法,它与高程静校正方法、地表一致性剩余静校正方法联合应用,可较好地解决黄土塬地区严重的静校正问题,可提高处理成果构造的可靠性,为解释提供可靠的基础资料。因此,该发明在黄土塬地区地震勘探的处理中具有很好的应用前景。
本发明是通过以下技术方案实现的:
一种地震勘探资料中波长剩余静校正处理方法,包括:
(1)输入构造形态控制的高程静校正数据和叠加速度,然后抽取共检波点道集数据、共炮点道集数据;
(2)速度扫描获取线性动校正的速度;
(3)优选参与叠加的炮检距范围及叠加;
(4)检波点和炮点中波长静校正量提取及应用;
(5)多次迭代求取准确的地表一致性剩余静校正量并应用。
所述步骤(2)是这样实现的:
利用速度扫描方式对步骤(1)得到的共检波点道集数据或共炮点道集数据进行线性动校正叠加,根据初至波校平的情况,选取线性动校正的速度。
所述步骤(3)是这样实现的:
对步骤(2)得到线性动校正的速度进行线性动校正后,避开近炮检距的直达波部分,保留稳定的折射波部分,并对线性动校正后的共检波点道集数据或线性动校正后的共炮点道集数据进行叠加。
所述步骤(4)是这样实现的:
基于共检波点道集数据叠加的初至波趋势拾取用于提取中波长剩余静校正的零时间线,基于共检波点道集数据叠加的初至波拾取用于提取中波长剩余静校正的校正时间线,将拾取的校正时间线减去零时间线,即可获得检波点的中波长静校正量;
在可控震源激发无井深,且炮点与检波点的站号相同的情况下,只提取检波点的中波长静校正量,然后将其映射并应用到炮点;
如果炮点与检波点的站号不同,或在炸药激发存在井深的情况下,需分别提取炮点和检波点的中波长静校正量;炮点的中波长静校正量提取方法与检波点的中波长静校正量提取方法相同。
所述步骤(5)是这样实现的:
抽取共深度点(CDP)道集数据,在应用中波长静校正的基础上,利用前面高程静校正处理后求取的速度谱、通过速度拾取获得准确叠加速度(步骤(1)中提供),采用多次迭代地表一致性剩余静校正的处理方法,便可逐步获得更为精确的剩余静校正量。
与现有技术相比,本发明的有益效果是:本发明方法适应于应用高程静校正方法后,且存在大值的中波长静校正量的情况,通过将它与高程静校正方法、地表一致性剩余静校正方法联合应用,克服了黄土塬地区地震资料的静校正技术难题,它能使求取的静校正量更加准确。通过该发明应用,可以提高资料处理成果的品质和构造的可靠性。
附图说明
图1是本发明方法的步骤框图。
图2是线性动校正后某一共检波点道集数据(全炮检距范围)。
图3是共检波点道集的叠加剖面(采用1600米/秒的线性动校正速度,检波点站号160附近初至波不连续)。
图4是共检波点道集的叠加剖面(采用1650米/秒的线性动校正速度,总体初至波连续)。
图5是共检波点道集的叠加剖面(采用1700米/秒的线性动校正速度,检波点站号580附近初至波不连续)。
图6是线性动校正后某一共检波点道集数据(优选部分炮检距范围,与图2对应)。
图7是共炮点道集数据叠加剖面。
图8是基于共检波点道集数据叠加的初至波趋势拾取用于提取剩余静校正的“零时间线”示意图。
图9是基于共检波点道集数据叠加的初至波拾取用于提取剩余静校正的“校正时间线”的示意图。
图10是提取的炮点中波长剩余静校正量。
图11是提取的检波点中波长剩余静校正量。
图12是中波长剩余静校正应用前的炮记录。
图13是中波长剩余静校正应用后的炮记录(与图12对应)。
图14是高程静校正后的叠加剖面(部分地段信噪比低,部分地段信噪比高,但构造格局基础好)。
图15是采用本发明的基于初至波叠加提取的中波长剩余静校正的叠加剖面。
图16是应用本发明静校正量的基础上进行多次迭代剩余静校正后的叠加剖面(与图17对比,构造可靠)。
图17是数值模型剖面。
图18是应用层析静校正的叠加剖面(与图16、图17对比,构造严重失真)。
具体实施方式
下面结合附图对本发明作进一步详细描述:
本发明可在一定程度上解决中波长的静校正难题,为进一步短波长的剩余静校正处理奠定良好的基础,最终获得较为真实可靠的剩余静校正量。
如图1所示,本方法包括:
(1)输入构造形态控制的高程静校正数据(通过执行高程静校正,质量控制标志层的构造形态可靠即可)和叠加速度,抽取共检波点道集数据、共炮点道集数据;
为了获得中波长的剩余静校正量,在高程静校正处理基础上,首先需要将数据抽取成共检波点道集数据(见图2)、共炮点道集数据。
(2)速度扫描获取线性动校正的速度
利用速度扫描方式对步骤(1)得到的共检波点道集数据或共炮点道集数据进行线性动校正叠加,根据初至波校平的情况,选取线性动校正的速度(看初至波的波组是否连续,如图3、图4、图5所示,速度合适时初至波连续,如图4),以便叠加可以获得可靠清晰的初至波。
(3)优选参与叠加的炮检距范围及叠加
共检波点道集数据进行线性动校正后(利用现有的线性动校正处理模块对步骤(2)是得到线性动校正的速度进行动校正),避开近炮检距(指避开小炮检距的直达波对折射初至波的影响)的直达波部分(见图6),保留稳定的折射波部分,并进行叠加,以便提取精确的中波长静校正量。
(4)检波点和炮点中波长静校正量提取及应用
在炮点与检波点的站号不同,或炸药激发存在井深的情况下,炮点和检波点的中波长静校正量需要分别提取。故需要分别进行共检波点道集数据叠加、共炮点道集数据叠加(见图7)。应用是把提取的静校正量执行掉,即对数据进行时移处理。
在可控震源激发无井深,且炮点与检波点的站号相同的情况下,炮点和检波点的中波长静校正量是相同的,只需提取检波点的静校正量,然后将其映射并应用到炮点即可(通过相同的站号的校正量可以拷贝,如果炮点与检波点的站号不同,炮点的静校正量提取方法与检波点的静校正量提取方法相同)。
图8给出了基于共检波点道集数据叠加的初至波趋势拾取用于提取中波长剩余静校正的“零时间线”,图9给出了基于共检波点道集数据叠加的初至波拾取用于提取中波长剩余静校正的“校正时间线”,将拾取的“校正时间线”减去“零时间线”,即可获得检波点的中波长静校正量(见图10)。同样的方法用在共炮点道集数据叠加剖面上,可以获得炮点的中波长静校正量(见图11)。
图12、图13分别给出了中波长剩余静校正应用前后的炮记录,可见应用后的记录上反射波双曲线规律得到改善。
图14是基于标志层构造形态控制的高程静校正叠加剖面,可见其构造格局基础较好,部分地段信噪比高,但部分地段信噪比低。
图15是采用本发明的基于初至波叠加提取的中波长剩余静校正的叠加剖面,可见剖面中波组的连续性得到了一定改善。
(5)多次迭代求取准确的静校正应用
在应用中波长静校正的基础上,利用前面处理获得的准确叠加速度(在高程静校正后进行速度谱求取,在信噪比高的地段进行速度分析。),采用多次迭代地表一致性剩余静校正的处理方法(由于低速带的速度和厚度在横向上的变化,使野外表层参数测量不准,或无法测量,故使野外静校正后,爆炸点和接收点的静校正量还残存着或正或负的误差,这个误差称为“剩余静校正量”。静校正采用地表一致性模型,即炮点和检波点在地面同一位置时,它们具有相同的静校正量,称为“地表一致性剩余静校正”。),便可逐步获得更为精确的剩余静校正量(将求取的静校正量应用后,再次进行剩余静校正处理。),调用处理系统的叠加模块,获得高品质的叠加剖面(叠加剖面是检验静校正量质量的有效手段(见图16)。
图16为应用本发明静校正量的基础上进行多次地表一致性剩余静校正后的叠加剖面,图17为数值模型剖面,图18为应用层析静校正的叠加剖面。对比图16、图17和图18可见:图16与图17的界面形态相似度高,图18与图17的界面形态相似度低,说明该发明方法优于层析静校正方法,可在一定程度上解决复杂的静校正问题,获得信噪比较高且构造形态可靠的处理成果。
上述技术方案只是本发明的一种实施方式,对于本领域内的技术人员而言,在本发明公开了应用方法和原理的基础上,很容易做出各种类型的改进或变形,而不仅限于本发明上述具体实施方式所描述的方法,因此前面描述的方式只是优选的,而并不具有限制性的意义。
Claims (5)
1.一种地震勘探资料中波长剩余静校正处理方法,其特征在于:所述方法包括:
(1)输入构造形态控制的高程静校正数据和叠加速度,然后抽取共检波点道集数据、共炮点道集数据;
(2)速度扫描获取线性动校正的速度;
(3)优选参与叠加的炮检距范围及叠加;
(4)检波点和炮点中波长静校正量提取及应用;
(5)多次迭代求取准确的地表一致性剩余静校正量并应用。
2.根据权利要求1所述的地震勘探资料中波长剩余静校正处理方法,其特征在于:所述步骤(2)是这样实现的:
利用速度扫描方式对步骤(1)得到的共检波点道集数据或共炮点道集数据进行线性动校正叠加,根据初至波校平的情况,选取线性动校正的速度。
3.根据权利要求2所述的地震勘探资料中波长剩余静校正处理方法,其特征在于:所述步骤(3)是这样实现的:
对步骤(2)得到线性动校正的速度进行线性动校正后,避开近炮检距的直达波部分,保留稳定的折射波部分,并对线性动校正后的共检波点道集数据或线性动校正后的共炮点道集数据进行叠加。
4.根据权利要求3所述的地震勘探资料中波长剩余静校正处理方法,其特征在于:所述步骤(4)是这样实现的:
基于共检波点道集数据叠加的初至波趋势拾取用于提取中波长剩余静校正的零时间线,基于共检波点道集数据叠加的初至波拾取用于提取中波长剩余静校正的校正时间线,将拾取的校正时间线减去零时间线,即可获得检波点的中波长静校正量;
在可控震源激发无井深,且炮点与检波点的站号相同的情况下,只提取检波点的中波长静校正量,然后将其映射并应用到炮点;
如果炮点与检波点的站号不同,或在炸药激发存在井深的情况下,需分别提取炮点和检波点的中波长静校正量;炮点的中波长静校正量提取方法与检波点的中波长静校正量提取方法相同。
5.根据权利要求4所述的地震勘探资料中波长剩余静校正处理方法,其特征在于:所述步骤(5)是这样实现的:
抽取共深度点道集数据,在应用中波长静校正的基础上,利用前面高程静校正处理后求取的速度谱、通过速度拾取获得准确叠加速度,采用多次迭代地表一致性剩余静校正的处理方法,便可逐步获得更为精确的剩余静校正量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410221656.6A CN105093326B (zh) | 2014-05-23 | 2014-05-23 | 一种地震勘探资料中波长剩余静校正处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410221656.6A CN105093326B (zh) | 2014-05-23 | 2014-05-23 | 一种地震勘探资料中波长剩余静校正处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105093326A true CN105093326A (zh) | 2015-11-25 |
CN105093326B CN105093326B (zh) | 2017-09-19 |
Family
ID=54574180
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410221656.6A Active CN105093326B (zh) | 2014-05-23 | 2014-05-23 | 一种地震勘探资料中波长剩余静校正处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105093326B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105866839A (zh) * | 2016-06-28 | 2016-08-17 | 中国矿业大学(北京) | 基于共姿态道集的静校正方法及装置 |
CN109188526A (zh) * | 2018-11-02 | 2019-01-11 | 中国石油天然气股份有限公司大港油田分公司 | 一种改进型地表一致性剩余静校正方法 |
CN110749923A (zh) * | 2018-07-24 | 2020-02-04 | 中国石油化工股份有限公司 | 一种基于范数方程提高分辨率的反褶积方法 |
CN111624655A (zh) * | 2019-02-27 | 2020-09-04 | 中国石油天然气集团有限公司 | 初至波剩余静校正量的确定方法及装置 |
CN111983685A (zh) * | 2020-07-21 | 2020-11-24 | 中国海洋大学 | τ-p域地表非一致性长波长静校正方法 |
CN112444879A (zh) * | 2019-08-30 | 2021-03-05 | 中国石油化工股份有限公司 | 地震勘探数据的处理方法及存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020075759A1 (en) * | 2001-01-21 | 2002-06-20 | Xixiang Zhou | Static correction method for exploration seismic data using first arrivals of seismic waves |
CN102998704A (zh) * | 2012-12-13 | 2013-03-27 | 石颖 | 一种地球物理勘探地震资料处理方法 |
CN103728665A (zh) * | 2012-10-12 | 2014-04-16 | 中国石油化工股份有限公司 | 一种精细剩余静校正方法 |
-
2014
- 2014-05-23 CN CN201410221656.6A patent/CN105093326B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020075759A1 (en) * | 2001-01-21 | 2002-06-20 | Xixiang Zhou | Static correction method for exploration seismic data using first arrivals of seismic waves |
CN103728665A (zh) * | 2012-10-12 | 2014-04-16 | 中国石油化工股份有限公司 | 一种精细剩余静校正方法 |
CN102998704A (zh) * | 2012-12-13 | 2013-03-27 | 石颖 | 一种地球物理勘探地震资料处理方法 |
Non-Patent Citations (3)
Title |
---|
王琪: "多域叠加交叉时折射静校正方法研究", 《中国优秀硕士学位论文全文数据库•信息科学辑》 * |
王西文 等: "三湖地区地震资料处理中的静校正方法研究", 《石油地球物理勘探》 * |
马义忠 等: "复杂地表低信噪比地震数据处理研究", 《地球物理学进展》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105866839A (zh) * | 2016-06-28 | 2016-08-17 | 中国矿业大学(北京) | 基于共姿态道集的静校正方法及装置 |
CN110749923A (zh) * | 2018-07-24 | 2020-02-04 | 中国石油化工股份有限公司 | 一种基于范数方程提高分辨率的反褶积方法 |
CN109188526A (zh) * | 2018-11-02 | 2019-01-11 | 中国石油天然气股份有限公司大港油田分公司 | 一种改进型地表一致性剩余静校正方法 |
CN111624655A (zh) * | 2019-02-27 | 2020-09-04 | 中国石油天然气集团有限公司 | 初至波剩余静校正量的确定方法及装置 |
CN112444879A (zh) * | 2019-08-30 | 2021-03-05 | 中国石油化工股份有限公司 | 地震勘探数据的处理方法及存储介质 |
CN112444879B (zh) * | 2019-08-30 | 2022-10-14 | 中国石油化工股份有限公司 | 地震勘探数据的处理方法及存储介质 |
CN111983685A (zh) * | 2020-07-21 | 2020-11-24 | 中国海洋大学 | τ-p域地表非一致性长波长静校正方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105093326B (zh) | 2017-09-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105093326A (zh) | 一种地震勘探资料中波长剩余静校正处理方法 | |
CN102778693B (zh) | 一种基于反射波层拉平提取并消除的绕射波分离处理方法 | |
CN102221709B (zh) | 基于地层参数信息的速度分析与动校正方法 | |
CN102841379B (zh) | 一种基于共散射点道集的叠前时间偏移与速度分析方法 | |
CN104459794B (zh) | 共反射点道集时变时间差值的校正方法及装置 | |
CN102901985B (zh) | 一种适用于起伏地表的深度域层速度修正方法 | |
CN104483705A (zh) | 一种三维剩余静校正方法 | |
Amoroso et al. | S‐wave identification by polarization filtering and waveform coherence analyses | |
CN102073064B (zh) | 一种利用相位信息提高速度谱分辨率的方法 | |
CN104199103B (zh) | 转换波静校正方法及装置 | |
CN103454681B (zh) | 评价三维地震观测系统成像效果的方法和设备 | |
CN103645503A (zh) | 一种三维时间域照明分析及振幅补偿方法 | |
CN104280775A (zh) | 一种基于全波形矢量偏移叠加的微地震监测定位方法 | |
CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
CN105093319A (zh) | 基于三维地震数据的地面微地震静校正方法 | |
CN102053260B (zh) | 获得地震纵波的方位速度的方法及处理地震数据的方法 | |
CN102338885B (zh) | 三分量vsp资料初至时间自动拾取方法 | |
CN102478666B (zh) | 一种提高可控震源地震记录初至拾取精度的方法 | |
CN103777242A (zh) | 一种深度聚焦和道集同相轴拉平联合的速度判别方法 | |
CN104133244A (zh) | 一种震源横波勘探的近地表静校正处理方法 | |
CN108375789A (zh) | 联合采集地震数据的同步匹配方法 | |
CN105425300B (zh) | 一种剩余静校正方法 | |
CN106338760B (zh) | 误差补偿的起伏地表偏移方法 | |
CN104076395A (zh) | 基于滤波组合的镜面能量提取与成像方法 | |
CN104536048A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |