CN103499836A - 空变-多时窗融合高精度剩余静校正方法 - Google Patents
空变-多时窗融合高精度剩余静校正方法 Download PDFInfo
- Publication number
- CN103499836A CN103499836A CN201310302572.0A CN201310302572A CN103499836A CN 103499836 A CN103499836 A CN 103499836A CN 201310302572 A CN201310302572 A CN 201310302572A CN 103499836 A CN103499836 A CN 103499836A
- Authority
- CN
- China
- Prior art keywords
- static correction
- residual static
- space
- window
- variant
- 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
技术领域
本发明属于石油地震勘探,对地震资料技术处理并解释的技术领域,具体地说涉及基于变差函数拟合重构的高精度三维静校正的空变—多时窗融合高精度剩余静校正方法。
背景技术
现有技术中,基础静校正以后,一般需要通过剩余静校正处理进一步消除短波长静校正量的影响。目前剩余静校正处理主要有三类方案:一类是基于相关类的地表一致性剩余静校正方法,通过相关等方案拾取地震道的时移量,通过将拾取的时移量进行反演获得炮、检点剩余静校正量,此类方法计算效率相对较高,但是当近地表结构复杂、地下构造复杂(即为“双复杂”区域)、资料信噪比偏低时,地表一致性假设条件不能被很好的满足,此时可能产生“周期跳跃”,此类方法通常会失效。第二类是基于非地表一致性假设的剩余静校正方案,目前常用的是TRIM剩余静校正,但是此类方法如果不能很好的控制剩余静校正参数,很容易产生假构造,特别是针对陌生区域,处理人员对地下构造认知度不高的时候,一般是慎用的。第三类是基于全局寻优思路的剩余静校正方案,此类方法理论上可以解决大周期的剩余静校正问题,但是此类方法计算效率偏低,特别是在低信噪比、复杂地下构造区域应用时,局部收敛过慢,全局搜索能力也会得到一定抑制,且其计算结果很难预知,对处理人员综合素质要求较高,因此针对大规模三维工区,这类方法也很少被应用。
为了避免双复杂及低信噪比区域因动校正速度不准而引起剩余静校正畸变,通常需要选取信噪比较高、同相轴连续性较强的构造区域建立模型道,这样才能降低由动校正速度不准而引起的对剩余静校正畸变,增强相关的稳定性,提高标准模型道的可靠性,从而提高剩余静校正量的计算精度。然而,复杂勘探区地下构造在空间的展布变化比较剧烈,很难保证在一个时窗内计算出高标准的模型道。按照以往的剩余静校正计算思路,我们通常会选取空间上固定的时窗进行剩余静校正量计算,但是如果选取大时窗,那么过大的时窗内,往往会引入一些动校正速度落实不准确的区域,从而在产生模型道过程中降低精度;如果选取的时窗过小,又不具备统计效应,同样会降低模型道的精度。
通过对野外实际资料的测试,采用常规的单一时窗地表一致性剩余静校正方法取得的叠加剖面信噪比虽然得到了一定程度的提高,但是局部地方细节并不清晰。
发明内容
本发明的目的在于为克服双复杂及低信噪比区域因动校正速度不准而引起的剩余静校正畸变,避免像Trim等非地表一致性算法容易产生假构造的弊端,而提供一种空变—多时窗融合高精度剩余静校正方法;该方法不仅可以消除速度误差对剩余静校正量计算的影响,并且在一定程度上可以有效解决非地表一致性静校正问题。
本发明的目的是通过如下方法来实现的:空变—多时窗融合高精度剩余静校正方法;该方法包括以下步骤:
A、获取准确的基础静校正量,彻底解决好长波长问题;
B、做好叠前去噪,提高信噪比工作;
C、对地下构造空间变化趋势进行三维空间的宏观掌控,以选取并确定空变时窗;
D、根据地下构造的纵横向变化趋势,确定空变时窗的个数及空变时窗的时间范围;
E、根据之前选取的空变时窗的时间范围,分别用不同的时窗计算全区的模型及剩余静校正量;
F、分别进行剩余静校正量的分解,依据不同时窗计算静校正量的叠加效果,对不同时窗计算的剩余静校正量进行三维空间的拼接;
G、采用高斯——赛德尔迭代方式统一将剩余静校正量进行分解,将步骤F拼接的剩余静校正量分解成各自炮点、检波点及CMP剩余静校正量。
本发明的有益效果是:
(1)采用本发明的方法得到的叠加剖面不仅提高了剖面的整体信噪比,剖面整体面貌较常规的单一时窗剩余静校正后的结果更加的清晰,局部细节也更加的清楚。
(2)采用本发明的方法不仅可以消除速度误差对剩余静校正量计算的影响,并且在一定程度上可以有效解决非地表一致性静校正问题。
(3)采用本发明的方法在双复杂区及低信噪比区是可行、有效的。
附图说明
图1是本发明在中非乍得Daniela东区块剩余静校正不同计算时窗选取(时窗1:0s—2s;时窗2);
图2是本发明在中非乍得Daniela东区块剩余静校正空变时窗的选取(时窗1:W1:700—2100ms,时窗2:W2:200—1000ms,时窗3:W3:200—1500ms);
图3是本发明三个空变时窗的空间分布范围(时窗1:W1:700—2100ms,时窗2:W2:200—1000ms,时窗3:W3:200—1500ms);
图4是本发明常规的单一时窗剩余静校正叠加剖面;
图5是本发明提供的空变—多时窗融合高精度剩余静校正叠加剖面。
具体实施方式
下面例举对本发明的方法作进一步详细阐述
空变—多时窗融合高精度剩余静校正方法;该方法步骤如下:
1)获取准确的基础静校正量,彻底解决好长波长问题;
2)做好叠前去噪等提高信噪比工作;
3)对地下构造空间变化趋势进行三维空间的宏观掌控,以便确定空变时窗的选取;
4)根据地下构造的纵横向变化趋势,确定空变时窗的个数及空变时窗的时间范围;分别为时窗1:700—2100ms、时窗2:200—1000ms、时窗3:200—1500ms;
5)根据之前选取的空变时窗的时间范围,分别用不同的时窗计算全区的模型及剩余静校正量;
6)分别进行剩余静校正量的分解,依据不同时窗计算静校正量的叠加效果,对不同时窗计算的剩余静校正量进行三维空间的拼接;
7)采用高斯——赛德尔迭代(一种常用的数学分解公式,具有计算效率高,分解稳定的特性)等方式统一将剩余静校正量进行分解,将步骤6拼接的剩余静校正量分解成各自炮点、检波点及CMP剩余静校正量。
本发明以中非乍得Daniela东区块三维实际资料为例进行了反复实验和野外实际资料的测试,根据前期叠加剖面及现场提供的叠前时间偏移剖面,可以看到,地下构造南北向变化剧烈,采用图1中的两种时窗计算(时窗1:0s—2s,时窗2)的剩余静校正量,叠加效果都不理想,不能使所有的同相轴都能得到有效加强,仍然存在残余剩余静校正量。
而依据本发明提出的新的静校正思路,分别采用了时窗1:W1:700—2100ms,时窗2:200—1000ms,时窗3:200—1500ms(如图2所示)三个计算时窗进行了剩余静校正量的计算,然后按照各自窗口进行剩余静校正量的分解,利用三套分离后的剩余静校正量进行叠加,叠加后依据叠加效果将不同时窗计算的剩余静校正量进行分析对比,在空间上找出三个时窗计算的剩余静校正量各自叠加效果较好的区域,然后进行空间范围的确定(如图3所示)。这种方法不仅可以消除速度误差对剩余静校正量计算的影响,并且在一定程度上可以有效解决非地表一致性静校正问题。由于每个时窗内采用的计算方法仍然是地表一致性算法,因此可以避免像Trim等非地表一致性算法容易产生的假构造。
采用常规的单一时窗地表一致性剩余静校正方法取得的叠加剖面信噪比虽然得到了一定程度的提高,但是局部地方细节并不清晰(如图4所示)。
而通过空变—多时窗融合高精度剩余静校正方法得到的叠加剖面不仅提高了剖面的整体信噪比,剖面整体面貌较常规的单一时窗剩余静校正后的结果更加的清晰,局部细节也更加的清楚(如图5所示)。因此,可以说明空变—多时窗融合高精度剩余静校正方法在双复杂区、低信噪比区是可行、有效的;通过上述地区的实验和野外资料测试足可以证明本发明方法设计思路的实现和应用效果。
Claims (1)
1.空变—多时窗融合高精度剩余静校正方法,其特征在于:该方法包括以下步骤:
A、获取准确的基础静校正量,彻底解决好长波长问题;
B、做好叠前去噪,提高信噪比工作;
C、对地下构造空间变化趋势进行三维空间的宏观掌控,以选取并确定空变时窗;
D、根据地下构造的纵横向变化趋势,确定空变时窗的个数及空变时窗的时间范围;
E、根据之前选取的空变时窗的时间范围,分别用不同的时窗计算全区的模型及剩余静校正量;
F、分别进行剩余静校正量的分解,依据不同时窗计算静校正量的叠加效果,对不同时窗计算的剩余静校正量进行三维空间的拼接;
G、采用高斯——赛德尔迭代方式统一将剩余静校正量进行分解,将步骤F拼接的剩余静校正量分解成各自炮点、检波点及CMP剩余静校正量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310302572.0A CN103499836B (zh) | 2013-07-18 | 2013-07-18 | 空变-多时窗融合高精度剩余静校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310302572.0A CN103499836B (zh) | 2013-07-18 | 2013-07-18 | 空变-多时窗融合高精度剩余静校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103499836A true CN103499836A (zh) | 2014-01-08 |
CN103499836B CN103499836B (zh) | 2016-04-06 |
Family
ID=49865059
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310302572.0A Expired - Fee Related CN103499836B (zh) | 2013-07-18 | 2013-07-18 | 空变-多时窗融合高精度剩余静校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103499836B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103837894A (zh) * | 2014-03-10 | 2014-06-04 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 获取剩余静校正量的方法 |
CN104914468A (zh) * | 2015-06-09 | 2015-09-16 | 中南大学 | 一种矿山微震信号p波初至时刻联合拾取方法 |
CN105629306A (zh) * | 2014-10-27 | 2016-06-01 | 中国石油化工股份有限公司 | 一种信噪比模型建立方法 |
CN105974467A (zh) * | 2016-05-05 | 2016-09-28 | 中国海洋大学 | 一种地震记录变时窗自动增益控制方法 |
CN106443829A (zh) * | 2016-08-26 | 2017-02-22 | 郑鸿明 | 一种近地表模型构建方法及装置 |
CN106547024A (zh) * | 2015-09-16 | 2017-03-29 | 中国石油化工股份有限公司 | 针对微地震射孔数据的剩余静校正量估计方法和装置 |
CN109975884A (zh) * | 2017-12-28 | 2019-07-05 | 核工业北京地质研究院 | 一种放射性地球物理测量数据融合方法 |
CN111624655A (zh) * | 2019-02-27 | 2020-09-04 | 中国石油天然气集团有限公司 | 初至波剩余静校正量的确定方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5157638A (en) * | 1992-01-13 | 1992-10-20 | Conoco Inc. | Method of deriving statics corrections from common reflection point gathers |
US6188964B1 (en) * | 1999-09-14 | 2001-02-13 | Ut-Battelle, Llc | Method for using global optimization to the estimation of surface-consistent residual statics |
CN101315427A (zh) * | 2007-05-29 | 2008-12-03 | 中国石油天然气集团公司 | 一种复杂地区地震勘探资料处理的方法及系统 |
CN102692651A (zh) * | 2011-03-24 | 2012-09-26 | 中国石油集团东方地球物理勘探有限责任公司 | 速度空变的初至波剩余静校正方法 |
-
2013
- 2013-07-18 CN CN201310302572.0A patent/CN103499836B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5157638A (en) * | 1992-01-13 | 1992-10-20 | Conoco Inc. | Method of deriving statics corrections from common reflection point gathers |
US6188964B1 (en) * | 1999-09-14 | 2001-02-13 | Ut-Battelle, Llc | Method for using global optimization to the estimation of surface-consistent residual statics |
CN101315427A (zh) * | 2007-05-29 | 2008-12-03 | 中国石油天然气集团公司 | 一种复杂地区地震勘探资料处理的方法及系统 |
CN102692651A (zh) * | 2011-03-24 | 2012-09-26 | 中国石油集团东方地球物理勘探有限责任公司 | 速度空变的初至波剩余静校正方法 |
Non-Patent Citations (2)
Title |
---|
李培明等: "《复杂地表建模及静校正技术》", 《科技成果管理与研究》 * |
潘文勇: "《基于改进模拟退火算法的剩余静校正及程序实现》", 《物探与化探》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103837894B (zh) * | 2014-03-10 | 2016-08-31 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 获取剩余静校正量的方法 |
CN103837894A (zh) * | 2014-03-10 | 2014-06-04 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 获取剩余静校正量的方法 |
CN105629306B (zh) * | 2014-10-27 | 2018-01-16 | 中国石油化工股份有限公司 | 一种信噪比模型建立方法 |
CN105629306A (zh) * | 2014-10-27 | 2016-06-01 | 中国石油化工股份有限公司 | 一种信噪比模型建立方法 |
CN104914468B (zh) * | 2015-06-09 | 2017-06-06 | 中南大学 | 一种矿山微震信号p波初至时刻联合拾取方法 |
CN104914468A (zh) * | 2015-06-09 | 2015-09-16 | 中南大学 | 一种矿山微震信号p波初至时刻联合拾取方法 |
CN106547024A (zh) * | 2015-09-16 | 2017-03-29 | 中国石油化工股份有限公司 | 针对微地震射孔数据的剩余静校正量估计方法和装置 |
CN106547024B (zh) * | 2015-09-16 | 2018-08-07 | 中国石油化工股份有限公司 | 针对微地震射孔数据的剩余静校正量估计方法和装置 |
CN105974467A (zh) * | 2016-05-05 | 2016-09-28 | 中国海洋大学 | 一种地震记录变时窗自动增益控制方法 |
CN105974467B (zh) * | 2016-05-05 | 2017-04-12 | 中国海洋大学 | 一种地震记录变时窗自动增益控制方法 |
CN106443829A (zh) * | 2016-08-26 | 2017-02-22 | 郑鸿明 | 一种近地表模型构建方法及装置 |
CN109975884A (zh) * | 2017-12-28 | 2019-07-05 | 核工业北京地质研究院 | 一种放射性地球物理测量数据融合方法 |
CN109975884B (zh) * | 2017-12-28 | 2020-10-23 | 核工业北京地质研究院 | 一种放射性地球物理测量数据融合方法 |
CN111624655A (zh) * | 2019-02-27 | 2020-09-04 | 中国石油天然气集团有限公司 | 初至波剩余静校正量的确定方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN103499836B (zh) | 2016-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103499836B (zh) | 空变-多时窗融合高精度剩余静校正方法 | |
CN107843925B (zh) | 一种基于修正相位的反射波波形反演方法 | |
US9632192B2 (en) | Method of processing seismic data by providing surface offset common image gathers | |
CN102841375A (zh) | 一种复杂条件下基于角度域共成像点道集的层析速度反演方法 | |
US10670750B2 (en) | Multistage full wavefield inversion process that generates a multiple free data set | |
Boonyasiriwat et al. | Applications of multiscale waveform inversion to marine data using a flooding technique and dynamic early-arrival windows | |
CN102841379B (zh) | 一种基于共散射点道集的叠前时间偏移与速度分析方法 | |
CN102841376A (zh) | 一种基于起伏地表的层析速度反演方法 | |
CN102884447A (zh) | Q层析成像方法 | |
CN102890290A (zh) | 一种起伏地表条件下的叠前深度偏移方法 | |
US11579323B2 (en) | Noise attenuation | |
CN102944896A (zh) | 表层调查数据的模型法静校正方法 | |
CN102937720A (zh) | 井控提高地震资料分辨率的方法 | |
CN110471113A (zh) | 基于非稳态地震资料的反演动校正方法、装置及存储介质 | |
CN104977615B (zh) | 一种基于模型统计拾取的深水obc资料多次波压制方法 | |
CN111060961B (zh) | 基于多信息约束反演的品质因子确定方法、装置及系统 | |
An et al. | Auto-pick first breaks with complex raypaths for undulate surface conditions | |
Montel et al. | RTM and Kirchhoff angle domain common-image gathers for migration velocity analysis | |
Shi et al. | A layer-stripping method for 3D near-surface velocity model building using seismic first-arrival times | |
CN111208568B (zh) | 一种时间域多尺度全波形反演方法及系统 | |
CN107193047A (zh) | 非地表一致性时变动态剩余静校正方法和装置 | |
CN107229074A (zh) | 低勘探程度海相探区的变速成图方法和装置 | |
CN106569278A (zh) | 一种多道相似相干速度谱计算方法 | |
Cheng et al. | Structure enhanced least-squares migration by deep learning based structural preconditioning | |
Gao et al. | Multimode wave-equation-based Scholte wave inversion |
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 |
Granted publication date: 20160406 Termination date: 20180718 |
|
CF01 | Termination of patent right due to non-payment of annual fee |