CN105372705B - 一种基于多波资料的地层切片方法 - Google Patents
一种基于多波资料的地层切片方法 Download PDFInfo
- Publication number
- CN105372705B CN105372705B CN201510707117.8A CN201510707117A CN105372705B CN 105372705 B CN105372705 B CN 105372705B CN 201510707117 A CN201510707117 A CN 201510707117A CN 105372705 B CN105372705 B CN 105372705B
- Authority
- CN
- China
- Prior art keywords
- data
- wavelength
- offset data
- divisions
- wave
- 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 38
- 230000008569 process Effects 0.000 claims abstract description 7
- 238000001615 p wave Methods 0.000 claims description 13
- 241001269238 Data Species 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 7
- 230000008021 deposition Effects 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 4
- 239000003550 marker Substances 0.000 abstract 1
- 238000004321 preservation Methods 0.000 abstract 1
- 230000010363 phase shift Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 208000035126 Facies Diseases 0.000 description 3
- 230000007704 transition Effects 0.000 description 3
- 230000007812 deficiency Effects 0.000 description 2
- 238000000280 densification Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 238000006467 substitution reaction 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/30—Analysis
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
本发明提出了一种基于多波资料的地层切片方法，包括：步骤1，对包含PP波及PS波的地震数据进行保幅处理，获取叠前道集的PP波分偏移距数据和PS波分偏移距数据；步骤2，联合全波列测井资料，根据PP波分偏移距数据和PS波分偏移距数据联合同时反演方法获得至少包括纵波阻抗、横波阻抗、泊松比在内的弹性参数体；步骤3，在等时性分析的控制下对PP波分偏移距数据追踪获取等时标志层；步骤4，将弹性参数体进行转换获得地质年代体；步骤5，根据地质年代体在等时标志层之间以等比例内插的方式对地层进行切片处理，获得具有近似年代地层意义的地层切片，通过多张地层切片的平面相解释，获得反映砂组尺度的沉积特征，恢复沉积演化过程。
Description
技术领域
本发明涉及石油地球物理勘探领域，尤指一种基于多波资料的地层切片方法。
背景技术
目前，地震沉积学是继地震地层学和层序地层学之后出现的一个新的学科。近年来国内外学者开展了广泛的研究。概括地讲，地震沉积学是一门主要利用地震信息和技术研究有关沉积岩及其形成过程的学科。具体来讲，它侧重于研究地震岩石学、地震地貌学、沉积结构、沉积史、沉积体系和沉积相平面展布(曾洪流，朱筱敏,朱如凯等，陆相坳陷型盆地地震沉积学研究规范，石油勘探与开发，2012，39(3)：275-284；董春梅，张宪国，林承焰，有关地震沉积学若干问题的探讨，石油地球物理勘探，2006，41(4)：405-409；刘化清，卫平生，李相博等，对地震沉积学理论创新的认识，岩性油气藏，2012，24(1)：7-11)。其主要的研究方法为：90度相位转换、分频解释和地层切片技术。90度相位转换的目的是相位转换后同相轴具有了岩性意义，从而使地震道与岩性测井更加具有一致性。基于地震数据的频率控制了同相轴倾角和内部反射结构这一理念，分频解释可以适当调整地震数据的频率确保解释的等时性。利用地层切片技术进行沉积时间界面的地震地貌学研究。
现行的所有的地层切片方法都是基于90度相移数据及叠后反演数据。之所以选择这两种数据，就是因为这两种数据可能具备区分岩性的功能，但实际上在复杂岩性区、致密油气区等这两种数据并不能有效区分岩性，那么基于这两种数据的地层切片进行地震沉积学研究就失去了意义。因此，现有的地层切片方法还存在一定的不足。
发明内容
为解决基于90度相移数据及叠后反演数据的地层切片不适用于复杂岩性区、致密油气区等地区的问题，本发明提供一种基于多波资料的地层切片方法，首次将多波资料应用到地震沉积学研究中，可以有效克服现有技术的不足，提高地层切片在不同地区的适用性。
本发明提出的基于多波资料的地层切片方法包括：步骤1，对包含PP波及PS波的地震数据进行保幅处理，获取叠前道集的PP波分偏移距数据和PS波分偏移距数据；步骤2，联合全波列测井资料，根据所述PP波分偏移距数据和所述PS波分偏移距数据联合同时反演方法获得至少包括纵波阻抗、横波阻抗、泊松比在内的弹性参数体；步骤3，在等时性分析的控制下，对所述PP波分偏移距数据追踪获取等时标志层；步骤4，将所述弹性参数体进行转换获得地质年代体；步骤5，根据所述地质年代体，在所述等时标志层之间以等比例内插的方式对地层进行切片处理，获得具有近似年代地层意义的地层切片，通过多张所述地层切片的平面相解释，获得反映砂组尺度的沉积特征，恢复沉积演化过程。
进一步的，在步骤1中，至少获取3个以上的叠前道集的PP波分偏移距数据和PS波分偏移距数据。
进一步的，步骤2，联合全波列测井资料，根据所述PP波分偏移距数据和所述PS波分偏移距数据联合同时反演方法获得至少包括纵波阻抗、横波阻抗、泊松比在内的弹性参数体，包括：利用所述全波列测井资料分别对所述PP波分偏移距数据和所述PS波分偏移距数据进行标定，将所述PS波分偏移距数据对齐到所述PP波分偏移距数据的时间域，根据所述PP波分偏移距数据和所述PS波分偏移距数据联合同时反演方法获得至少包括纵波阻抗、横波阻抗、泊松比在内的弹性参数体。
进一步的，步骤3，在等时性分析的控制下，对所述PP波分偏移距数据追踪获取等时标志层，包括：对所述PP波分偏移距数据进行追踪解释，获取大套标志层；对所述大套标志层进行小波分频，对小波分频后的数据进行地震倾角估算，获得高、低频的倾角；当高、低频的倾角差值低于阈值时，判定地震反射轴是等时的，在等时性分析的控制下，对所述PP波分偏移距数据追踪获取等时标志层。
进一步的，步骤4，将所述弹性参数体进行转换获得地质年代体，包括：将所述弹性参数体进行Wheeler域变换，获得地质年代体。
通过本发明提出的基于多波资料的地层切片方法，可以在地震PP+PS分偏移距资料的基础上，进行叠前PP+PS联合同时反演方法求取高分辨率高精度的纵波阻抗、横波阻抗、泊松比等弹性参数体，创新性地在弹性参数体上进行地层切片研究，可以解决复杂岩性区、致密油气区等地区利用90度相移数据或者叠后反演数据的地层切片不能区分岩性等问题，并将多波资料获得的纵波阻抗、横波阻抗、泊松比等弹性参数体应用到地震沉积学研究中，可以提高地震沉积学技术在复杂岩性区、致密油气区等不同地区的适用性。
附图说明
此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中：
图1为本发明一实施例的基于多波资料的地层切片方法流程图。
图2为某三维工区Line753多波叠前反演获得的弹性参数剖面。
图3为横波阻抗弹性参数体典型地层切片。
图4为泊松比弹性参数体典型地层切片。
具体实施方式
以下配合图示及本发明的较佳实施例，进一步阐述本发明为达成预定发明目的所采取的技术手段。
图1为本发明一实施例的基于多波资料的地层切片方法流程图。如图1所示，该方法包括：
步骤1，对包含PP波及PS波的地震数据进行保幅处理，获取叠前道集的PP(纵波)波分偏移距数据和PS(转换波)波分偏移距数据；
步骤2，联合全波列测井资料，根据所述PP波分偏移距数据和所述PS波分偏移距数据联合同时反演方法获得至少包括纵波阻抗、横波阻抗、泊松比在内的弹性参数体；
步骤3，在等时性分析的控制下，对所述PP波分偏移距数据追踪获取等时标志层；
步骤4，将所述弹性参数体进行转换获得地质年代体；
步骤5，根据所述地质年代体，在所述等时标志层之间以等比例内插的方式对地层进行切片处理，获得具有近似年代地层意义的地层切片，通过多张所述地层切片的平面相解释，获得反映砂组尺度的沉积特征，恢复沉积演化过程。
具体而言，在步骤1中，获取保幅处理的n(n≥3)个叠前道集的PP波分偏移距数据和n(n≥3)个PS波分偏移距数据。
在步骤2中，利用所述全波列测井资料分别对所述PP波分偏移距数据和所述PS波分偏移距数据进行标定，将所述PS波分偏移距数据对齐到所述PP波分偏移距数据的时间域，根据所述PP波分偏移距数据和所述PS波分偏移距数据联合同时反演方法获得至少包括纵波阻抗、横波阻抗、泊松比在内的弹性参数体。
可参考图2所示，图2为某三维工区Line753多波叠前反演获得的弹性参数剖面，相比于常规地震沉积学主要技术手段(90度相移)，多波叠前反演获得的弹性参数体具有高精度、高分辨率的特点，可以有效区分岩性，并可以识别有效储层。
在步骤3中，首先，对所述PP波分偏移距数据进行追踪解释，获取大套标志层。
然后，对所述大套标志层进行小波分频，对小波分频后的数据进行地震倾角估算，获得高、低频的倾角。
当高、低频的倾角差值低于阈值时，判定地震反射轴是等时的，在等时性分析的控制下，对所述PP波分偏移距数据追踪获取等时标志层。如果高、低频的倾角差高于阈值，则需要修改解释方案，直到满足上述条件，使得在等时性分析的控制下追踪出大套的等时标志层层位。
步骤4的技术方案是撇开常规的90度相移分频然后切片的地震沉积学研究思路，直接在叠前PP+PS波分偏移距数据联合同时反演获得的纵波阻抗、横波阻抗、泊松比等弹性参数体上进行Wheeler域变换，将其转换为地质年代体，然后可以进行地层切片研究。
例如，结合图3及图4所示，图3为横波阻抗弹性参数体典型地层切片示意图，图4为泊松比弹性参数体典型地层切片示意图。图3图4中数字表示为油气井，其中横波阻抗弹性数据体(图3)可以反映主砂体的展布，泊松比弹性参数体(图4)可以反映有效储层的分布范围，统计工区内72口油气井上的砂体厚度及含油气性并标定图3图4，图3符合率达到85.8％，图4符合率达到75.7％。由此可知，基于弹性参数数据体的地层切片可以有效反映研究区储层的平面展布及空间演化过程。
最后，步骤5可以进行地震沉积相分析，进行综合评价。具体而言，根据所述地质年代体，在所述等时标志层之间以等比例内插的方式对地层进行切片处理，获得具有近似年代地层意义的地层切片，通过多张所述地层切片的平面相解释，获得反映砂组尺度的沉积特征，恢复沉积演化过程。
通过本发明提出的基于多波资料的地层切片方法，可以在地震PP+PS分偏移距资料的基础上，进行叠前PP+PS联合同时反演方法求取高分辨率高精度的纵波阻抗、横波阻抗、泊松比等弹性参数体，创新性地在弹性参数体上进行地层切片研究，可以解决复杂岩性区、致密油气区等地区利用90度相移数据或者叠后反演数据的地层切片不能区分岩性等问题，并将多波资料获得的纵波阻抗、横波阻抗、泊松比等弹性参数体应用到地震沉积学研究中，可以提高地震沉积学技术在复杂岩性区、致密油气区等不同地区的适用性。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (5)
1.一种基于多波资料的地层切片方法，其特征在于，该方法包括：
步骤1，对包含PP波及PS波地震数据进行保幅处理，获取叠前道集的PP波分偏移距数据和PS波分偏移距数据；
步骤2，联合全波列测井资料，根据所述PP波分偏移距数据和所述PS波分偏移距数据联合同时反演方法获得至少包括纵波阻抗、横波阻抗、泊松比在内的弹性参数体；
步骤3，在等时性分析的控制下，对所述PP波分偏移距叠加数据追踪获取等时标志层；
步骤4，将所述弹性参数体进行转换获得地质年代体；
步骤5，根据所述地质年代体，在所述等时标志层之间以等比例内插的方式对地层进行切片处理，获得具有近似年代地层意义的地层切片，通过多张所述地层切片的平面相解释，获得反映砂组尺度的沉积特征，恢复沉积演化过程。
2.根据权利要求1所述的基于多波资料的地层切片方法，其特征在于，在步骤1中，至少获取3个以上的叠前道集的PP波分偏移距数据和PS波分偏移距数据。
3.根据权利要求1所述的基于多波资料的地层切片方法，其特征在于，步骤2，联合全波列测井资料，根据所述PP波分偏移距数据和所述PS波分偏移距数据联合同时反演方法获得至少包括纵波阻抗、横波阻抗、泊松比在内的弹性参数体，包括：
利用所述全波列测井资料分别对所述PP波分偏移距数据和所述PS波分偏移距数据进行标定，将所述PS波分偏移距数据对齐到所述PP波分偏移距数据的时间域，根据所述PP波分偏移距数据和所述PS波分偏移距数据联合同时反演方法获得至少包括纵波阻抗、横波阻抗、泊松比在内的弹性参数体。
4.根据权利要求1所述的基于多波资料的地层切片方法，其特征在于，步骤3，在等时性分析的控制下，对所述PP波分偏移距数据追踪获取等时标志层，包括：
对所述PP波分偏移距数据进行追踪解释，获取大套标志层；
对所述大套标志层进行小波分频，对小波分频后的数据进行地震倾角估算，获得高、低频的倾角；
当高、低频的倾角差值低于阈值时，判定地震反射轴是等时的，在等时性分析的控制下，对所述PP波分偏移距数据追踪获取等时标志层。
5.根据权利要求1所述的基于多波资料的地层切片方法，其特征在于，步骤4，将所述弹性参数体进行转换获得地质年代体，包括：
将所述弹性参数体进行Wheeler域变换，获得地质年代体。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510707117.8A CN105372705B (zh) | 2015-10-27 | 2015-10-27 | 一种基于多波资料的地层切片方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510707117.8A CN105372705B (zh) | 2015-10-27 | 2015-10-27 | 一种基于多波资料的地层切片方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105372705A CN105372705A (zh) | 2016-03-02 |
| CN105372705B true CN105372705B (zh) | 2017-10-27 |
Family
ID=55375053
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510707117.8A Active CN105372705B (zh) | 2015-10-27 | 2015-10-27 | 一种基于多波资料的地层切片方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105372705B (zh) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105717542B (zh) * | 2016-04-11 | 2019-02-22 | 北京科胜伟达石油科技股份有限公司 | 基于三维地震的Wheeler域切片砂体雕刻方法 |
| US11604909B2 (en) | 2019-05-28 | 2023-03-14 | Chevron U.S.A. Inc. | System and method for accelerated computation of subsurface representations |
| US11249220B2 (en) | 2019-08-14 | 2022-02-15 | Chevron U.S.A. Inc. | Correlation matrix for simultaneously correlating multiple wells |
| US11187826B2 (en) | 2019-12-06 | 2021-11-30 | Chevron U.S.A. Inc. | Characterization of subsurface regions using moving-window based analysis of unsegmented continuous data |
| US11010969B1 (en) | 2019-12-06 | 2021-05-18 | Chevron U.S.A. Inc. | Generation of subsurface representations using layer-space |
| US10984590B1 (en) | 2019-12-06 | 2021-04-20 | Chevron U.S.A. Inc. | Generation of subsurface representations using layer-space |
| US11320566B2 (en) | 2020-01-16 | 2022-05-03 | Chevron U.S.A. Inc. | Multiple well matching within subsurface representation |
| US11263362B2 (en) | 2020-01-16 | 2022-03-01 | Chevron U.S.A. Inc. | Correlation of multiple wells using subsurface representation |
| US11397279B2 (en) | 2020-03-27 | 2022-07-26 | Chevron U.S.A. Inc. | Comparison of wells using a dissimilarity matrix |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102866426A (zh) * | 2012-09-20 | 2013-01-09 | 成都晶石石油科技有限公司 | 一种利用avo大角度道集分析岩体油气信息的方法 |
| CN103454686A (zh) * | 2013-08-12 | 2013-12-18 | 中国石油天然气股份有限公司 | 基于地层切片的小尺度沉积相进行储层预测的方法及系统 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2965066B1 (fr) * | 2010-09-20 | 2012-10-26 | Total Sa | Procede d'estimation de parametres elastiques |
| CN102830423B (zh) * | 2012-08-30 | 2015-04-08 | 中国石油天然气股份有限公司 | 储层流体检测方法和储层流体检测装置 |
| CN103439739B (zh) * | 2013-04-08 | 2016-08-17 | 中国石油集团东方地球物理勘探有限责任公司 | 地球物理勘探用纵横波匹配方法及匹配装置 |
| CN103454678B (zh) * | 2013-08-12 | 2016-02-10 | 中国石油天然气股份有限公司 | 一种地震切片等时性的确定方法及系统 |
-
2015
- 2015-10-27 CN CN201510707117.8A patent/CN105372705B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102866426A (zh) * | 2012-09-20 | 2013-01-09 | 成都晶石石油科技有限公司 | 一种利用avo大角度道集分析岩体油气信息的方法 |
| CN103454686A (zh) * | 2013-08-12 | 2013-12-18 | 中国石油天然气股份有限公司 | 基于地层切片的小尺度沉积相进行储层预测的方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105372705A (zh) | 2016-03-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105372705B (zh) | 一种基于多波资料的地层切片方法 | |
| CN101907727B (zh) | 一种面波多分量转换波静校正方法 | |
| Kandilarov et al. | Crustal structure of the ultra-slow spreading Knipovich Ridge, North Atlantic, along a presumed ridge segment center | |
| WO2018044960A1 (en) | Attenuation of multiple reflections | |
| Maresh et al. | Seismic attenuation of Atlantic margin basalts: Observations and modeling | |
| Lü et al. | The role of rifting in the development of the continental margins of the southwest subbasin, South China Sea: Insights from an OBS experiment | |
| CN104316966A (zh) | 一种流体识别方法及系统 | |
| Godfrey et al. | Shear-wave velocity structure of the Tongariro Volcanic Centre, New Zealand: fast Rayleigh and slow Love waves indicate strong shallow anisotropy | |
| Maries et al. | Cross-profile seismic data acquisition, imaging, and modeling of iron-oxide deposits: A case study from Blötberget, south-central Sweden | |
| Carr et al. | Shear-wave studies in glacial till | |
| CN104133244A (zh) | 一种震源横波勘探的近地表静校正处理方法 | |
| CN105527648A (zh) | 用于各向异性参数反演的敏感度矩阵的计算方法及系统 | |
| Ding et al. | Reverse-time ray-tracing method for microseismic source localization | |
| Wang et al. | Research on the anisotropy of gas hydrate reservoirs in South China Sea | |
| Bown | Legacy seismic investigations of karst surfaces: Implications for heavy oil extraction from the Devonian Grosmont formation, northeastern Alberta, Canada | |
| Ghalayini | Geophysical Investigation of Carrizo Formation by Using Two-Dimensional Seismic Surveys in the Tullos-Urania Oilfield in LaSalle Parish, LA | |
| CN105372708B (zh) | 一种基于等地质时间的致密薄储层预测方法 | |
| Snavely et al. | Stratigraphic and structural interpretations using diffraction seismograms | |
| Hajnal et al. | The earth's crust under the Williston Basin in eastern Saskatchewan and western Manitoba | |
| Campbell | Layer stripping the response from sedimentary basins in teleseismic data using transfer functions | |
| Sihoyiya | Cost-effective and novel seismic methods for mineral and coal exploration: Examples from Witwatersrand goldfields and Bushveld Complex | |
| Alali et al. | Data Conditioning of a Modern Midcontintent Data Volume using POCS Five-Dimensional Interpolation | |
| Velasquez-Espejo | 3D multicomponent seismic characterization of a clastic reservoir in the Middle Magdalena Valley Basin, Colombia | |
| Györfi et al. | High-resolution 2D and 3D seismic imaging of structurally complex hardrock environments hosting high-grade uranium ore, Athabasca Basin, Canada | |
| Boniface | Feasibility study of basin structure mapping using passive seismic exploration technique case study: Ahero area, Nyanza basin, Kisumu county, Kenya |
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 |