CN106353796A - 一种提高地面地震资料分辨率方法 - Google Patents
一种提高地面地震资料分辨率方法 Download PDFInfo
- Publication number
- CN106353796A CN106353796A CN201510423011.5A CN201510423011A CN106353796A CN 106353796 A CN106353796 A CN 106353796A CN 201510423011 A CN201510423011 A CN 201510423011A CN 106353796 A CN106353796 A CN 106353796A
- Authority
- CN
- China
- Prior art keywords
- frequency
- section
- seismic data
- seismic
- resolution ratio
- 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.)
- Pending
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及勘探地球物理资料处理技术领域,具体是一种提高地面地震资料分辨率方法。本发明首先通过小波包分频方法将原始的地震剖面分解成高频剖面和低频剖面，进一步根据基于系统辨识提高地震资料分辨率方法利用测井资料作为约束提高高频剖面分辨率。最后将拓频后的高频剖面与原始的低频剖面叠加重构获得最终的拓频地震剖面。常规基于系统辨识理论拓频方法会造成拓频后的地震资料低频信息有所损失，改变了原始地震资料中所包含的低频成分，而本发明方法获得的地震剖面更好保持了原始剖面中低频能量成分，拓频结果更加真实可靠。
Description
技术领域
本发明涉及勘探地球物理资料处理技术领域,具体是一种提高地面地震资料分辨率方法。
背景技术
测井、井间地震、VSP和地面地震等不同类别的地球物理技术可以对地下同一目标地质体进行不同尺度性质的反映，通过它们之间的联合作用提高地震资料的分辨率是一个新的研究方向。声波测井资料在井位置处具有最高的分辨率和最大程度的分辨可靠性，由其得到的反射系数序列可以认为是地震子波为单位脉冲响应条件下的地震记录，而且测井资料几乎在所有的区域都大量存在，另外由于地层介质的滞弹性和非均质性，地震波在地层中的传播要经历吸收作用，这种吸收作用与信号的频率有关，即地层对高频信号的吸收作用要远远大于对低频信号的作用。基于此刘浩杰等提出了“基于系统辨识提高地震资料分辨率”(参考文献：《油气地球物理》2010年4月第8卷第2期)。该方法基于系统辨识理论利用测井资料提高地震资料分辨率，但是直接利用该方法进行地震资料拓频，虽然能够提高地震资料分辨率，但在多数情况下会造成拓频后的地震资料低频信息有所损失，改变了原始地震资料中所包含的低频成分，造成拓频后的地震剖面失真。
发明内容
本发明的目的就是针对现有技术存在的上述问题，提供一种采用小波包分频与“基于系统辨识提高地震资料分辨率”方法相结合的提高地面地震资料分辨率方法。
本发明的提高地面地震资料分辨率方法，包括以下步骤:
步骤1：选择“sym4”小波基对原始地震剖面进行二级小波包分解，获得从低到高四个地震剖面P1,P2,P3,P4；
步骤2：将地震剖面P1作为低频剖面PL不做处理，同时将剖面P2,P3,P4进行叠加获得高频剖面PH，其计算公式如下
PH＝P2+P3+P4 公式1
步骤3：利用测井资料依据基于系统辨识拓频方法对高频剖面PH进行提高分辨率处理获得拓频高频剖面PHH；
步骤4：将步骤2中获得的低频剖面PL与步骤3中获得的拓频高剖面PHH进行叠加获得最总的拓频剖面PF，其计算公式如下
PF＝PHH+PL
步骤5：重复步骤1-4完成所有地震剖面的处理，获得最终的三维拓频数据体。
该方法将原始地震数据进行小波包分频，分频后的高频部分利用基于系统辨识拓频方法进行拓频，最后将拓频后的高频部分与原始低频部分叠加重构，最终获得保真度高的高分辨率地面地震剖面。
附图说明
图1本发明的技术流程图；
图2原始地震剖面截图；
图3小波包分解第一剖面P1截图；
图4小波包分解第二剖面P2截图；
图5小波包分解第三剖面P3截图；
图6小波包分解第四剖面P4截图；
图7叠加高频剖面PH截图；
图8拓频高频剖面PHH截图；
图9小波包分频拓频剖面截图；
图10常规拓频剖面截图。
具体实施方式
下面结合附图对本发明作进一步详细描述：
参照图1，本发明的具体流程为：
步骤1：选择”sym4”小波基对原始地震剖面(如图2所示)进行二级小波包分解，获得从低到高四个地震剖面P1,P2,P3,P4(如图3-6所示)。
步骤2：将地震剖面P1作为低频剖面PL(如图3所示)不做处理，同时将剖面P2,P3,P4进行叠加获得高频剖面PH(如图7所示)。
步骤3：利用高频测井资料依据基于系统辨识理论拓频方法对高频剖面PH进行提高分辨率处理获得拓频高频剖面PHH(如图8所示)。
步骤4：将步骤2中获得的低频剖面PL(如图3所示)与步骤3中获得的拓频高剖面PHH进行叠加获得最终的拓频剖面PF(如图9所示)
步骤5：重复步骤1-4完成所有地震剖面的处理，获得最终的三维拓频地震数据体。
利用常规基于系统辨识理论拓频方法获得了拓频后的地震剖面(如图10)和三维地震数据体。
从图3-6可以看出经过小波包分频后的4个剖面具有不同的主频，它们所凸显的信息存在着明显差异，通过小波包分解达到了对地震剖面进行分频的目的。进一步进行高低频剖面重构后可以看出，低频剖面PL(如图3)突出显示了地震剖面的强反射层，包含更多的地层构造信息，高频剖面PH(如图7)中包含更多的是地层微幅构造信息。
对比分析图1、9、10可以看出两种方法得到的拓频资料都能有效区分复合波、加强弱反射，使拓频地震剖面细节更加丰富，拓频后的地震资料更有利于研究岩性的空间变化；对比图9和10可以看出：常规基于系统辨识理论拓频方法会造成拓频后的地震资料低频信息有所损失，改变了原始地震资料中所包含的低频成分，而本发明方法获得的地震剖面更好保持了原始剖面中低频能量成分，拓频结果更加真实可靠。
Claims (1)
1.一种提高地面地震资料分辨率方法，其特征在于包括以下步骤:
步骤1：选择“sym4”小波基对原始地震剖面进行二级小波包分解，获得从低到高四个地震剖面P1,P2,P3,P4；
步骤2：将地震剖面P1作为低频剖面PL不做处理，同时将剖面P2,P3,P4进行叠加获得高频剖面PH，其计算公式如下
PH＝P2+P3+P4 公式1
步骤3：利用测井资料依据基于系统辨识拓频方法对高频剖面PH进行提高分辨率处理获得拓频高频剖面PHH；
步骤4：将步骤2中获得的低频剖面PL与步骤3中获得的拓频高剖面PHH进行叠加获得最总的拓频剖面PF，其计算公式如下
PF＝PHH+PL
步骤5：重复步骤1-4完成所有地震剖面的处理，获得最终的三维拓频数据体。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510423011.5A CN106353796A (zh) | 2015-07-17 | 2015-07-17 | 一种提高地面地震资料分辨率方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510423011.5A CN106353796A (zh) | 2015-07-17 | 2015-07-17 | 一种提高地面地震资料分辨率方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106353796A true CN106353796A (zh) | 2017-01-25 |
Family
ID=57842381
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510423011.5A Pending CN106353796A (zh) | 2015-07-17 | 2015-07-17 | 一种提高地面地震资料分辨率方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106353796A (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110646842A (zh) * | 2018-06-27 | 2020-01-03 | 中国石油化工股份有限公司 | 一种地震资料叠前偏移的面波压制的方法 |
| CN112835095A (zh) * | 2019-11-25 | 2021-05-25 | 中国石油天然气股份有限公司 | 地震资料低幅度构造成图方法及装置 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102323616A (zh) * | 2011-06-08 | 2012-01-18 | 浙江大学 | 提高灰岩出露区地震数据分辨率的分频匹配方法 |
| CN104216014A (zh) * | 2014-09-25 | 2014-12-17 | 郭平 | 一种地震信号分频处理方法 |
-
2015
- 2015-07-17 CN CN201510423011.5A patent/CN106353796A/zh active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102323616A (zh) * | 2011-06-08 | 2012-01-18 | 浙江大学 | 提高灰岩出露区地震数据分辨率的分频匹配方法 |
| CN104216014A (zh) * | 2014-09-25 | 2014-12-17 | 郭平 | 一种地震信号分频处理方法 |
Non-Patent Citations (6)
| Title |
|---|
| HAOJIE LIU ET AL.: "High Frequency Restoration of Surface Seismic Data Based on System Identification", 《SEG LAS VEGAS 2008 ANNUAL MEETING》 * |
| LIU SHUCHONG ET AL.: "Seismic signals wavelet packet de-noising method based on improved threshold function and adaptive threshold", 《COMPUTER MODELLING & NEW TECHNOLOGIES》 * |
| 刘浩杰等: "基于系统辨识提高地震资料分辨率", 《油气地球物理》 * |
| 刘浩杰等: "基于系统辨识提高地震资料分辨率研究", 《地球物理学进展》 * |
| 刘霞等: "基于小波阈值的地震信号去噪处理", 《科学技术与工程》 * |
| 杨培辉: "小波包技术研究及储层预测评价", 《中国优秀博硕士学位论文全文数据库（硕士） 基础科学辑》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110646842A (zh) * | 2018-06-27 | 2020-01-03 | 中国石油化工股份有限公司 | 一种地震资料叠前偏移的面波压制的方法 |
| CN112835095A (zh) * | 2019-11-25 | 2021-05-25 | 中国石油天然气股份有限公司 | 地震资料低幅度构造成图方法及装置 |
| CN112835095B (zh) * | 2019-11-25 | 2023-08-22 | 中国石油天然气股份有限公司 | 地震资料低幅度构造成图方法及装置 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Regone et al. | Geologic model building in SEAM Phase II—Land seismic challenges | |
| CN104570125B (zh) | 一种利用井数据提高成像速度模型精度的方法 | |
| CN102736107B (zh) | 一种能量约束非均质储层厚度识别系统 | |
| CN101158724B (zh) | 基于偶极小波的储层厚度预测方法 | |
| CN103675897B (zh) | 一种地震绕射波分离成像方法 | |
| CN101334483B (zh) | 一种在地震数据处理中衰减瑞雷波散射噪声的方法 | |
| CN104237945B (zh) | 一种地震资料自适应高分辨处理方法 | |
| CN103527184A (zh) | 一种白云岩储层的预测方法和系统 | |
| CN104516018A (zh) | 一种地球物理勘探中岩性约束下的孔隙度反演方法 | |
| Koesoemadinata et al. | Seismic reservoir characterization in Marcellus shale | |
| CN107390270B (zh) | 一种基于弹性波逆时偏移ADCIGs的AVA分析方法 | |
| Yang et al. | Integrated application of 3D seismic and microseismic data in the development of tight gas reservoirs | |
| CN104570067A (zh) | 一种地球物理勘探中相控地震反演方法 | |
| CN107255831A (zh) | 一种叠前频散属性的提取方法 | |
| CN103792573A (zh) | 一种基于频谱融合的地震波阻抗反演方法 | |
| CN105116449B (zh) | 一种弱反射储层的识别方法 | |
| CN103926619A (zh) | 一种三维vsp数据的逆时偏移方法 | |
| CN109738951A (zh) | 一种基于地震同相轴子波谱的时变反褶积方法 | |
| CN106405637A (zh) | 一种用于砂泥岩薄互层砂体厚度预测的方法及系统 | |
| CN104297800A (zh) | 一种自相控叠前反演方法 | |
| CN105092343B (zh) | 去除薄层调谐效应的方法和识别预测薄储层及气层的方法 | |
| CN102385066B (zh) | 一种叠前地震定量成像方法 | |
| CN104422960A (zh) | 基于信号低频强异常自适应提取的地震资料流体识别方法 | |
| CN106353796A (zh) | 一种提高地面地震资料分辨率方法 | |
| Ningkai et al. | Stepped and detailed seismic prediction of shallow-thin reservoirs in Chunfeng oilfield of Junggar Basin, NW China |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170125 |