CN105372705A - 一种基于多波资料的地层切片方法 - Google Patents

一种基于多波资料的地层切片方法 Download PDF

Info

Publication number
CN105372705A
CN105372705A CN201510707117.8A CN201510707117A CN105372705A CN 105372705 A CN105372705 A CN 105372705A CN 201510707117 A CN201510707117 A CN 201510707117A CN 105372705 A CN105372705 A CN 105372705A
Authority
CN
China
Prior art keywords
data
wavelength
offset data
division offset
division
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
Application number
CN201510707117.8A
Other languages
English (en)
Other versions
CN105372705B (zh
Inventor
窦玉坛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Natural Gas Co Ltd
Original Assignee
China Petroleum and Natural Gas Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Petroleum and Natural Gas Co Ltd filed Critical China Petroleum and Natural Gas Co Ltd
Priority to CN201510707117.8A priority Critical patent/CN105372705B/zh
Publication of CN105372705A publication Critical patent/CN105372705A/zh
Application granted granted Critical
Publication of CN105372705B publication Critical patent/CN105372705B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis

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域变换,获得地质年代体。
CN201510707117.8A 2015-10-27 2015-10-27 一种基于多波资料的地层切片方法 Active CN105372705B (zh)

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 true CN105372705A (zh) 2016-03-02
CN105372705B 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)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105717542A (zh) * 2016-04-11 2016-06-29 北京科胜伟达石油科技股份有限公司 基于三维地震的Wheeler域切片砂体雕刻方法
US10984590B1 (en) 2019-12-06 2021-04-20 Chevron U.S.A. Inc. Generation of subsurface representations using layer-space
US11010969B1 (en) 2019-12-06 2021-05-18 Chevron U.S.A. Inc. Generation of subsurface representations using layer-space
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
US11249220B2 (en) 2019-08-14 2022-02-15 Chevron U.S.A. Inc. Correlation matrix for simultaneously correlating multiple wells
US11263362B2 (en) 2020-01-16 2022-03-01 Chevron U.S.A. Inc. Correlation of multiple wells using subsurface representation
US11320566B2 (en) 2020-01-16 2022-05-03 Chevron U.S.A. Inc. Multiple well matching within subsurface representation
US11397279B2 (en) 2020-03-27 2022-07-26 Chevron U.S.A. Inc. Comparison of wells using a dissimilarity matrix
US11604909B2 (en) 2019-05-28 2023-03-14 Chevron U.S.A. Inc. System and method for accelerated computation of subsurface representations

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102830423A (zh) * 2012-08-30 2012-12-19 中国石油天然气股份有限公司 储层流体检测方法和储层流体检测装置
CN102866426A (zh) * 2012-09-20 2013-01-09 成都晶石石油科技有限公司 一种利用avo大角度道集分析岩体油气信息的方法
US20130176822A1 (en) * 2010-09-20 2013-07-11 Total Sa Method for estimating elastic parameters by inverting 4d seismic measurements
CN103439739A (zh) * 2013-04-08 2013-12-11 中国石油集团东方地球物理勘探有限责任公司 地球物理勘探用纵横波匹配方法及匹配装置
CN103454686A (zh) * 2013-08-12 2013-12-18 中国石油天然气股份有限公司 基于地层切片的小尺度沉积相进行储层预测的方法及系统
CN103454678A (zh) * 2013-08-12 2013-12-18 中国石油天然气股份有限公司 一种地震切片等时性的确定方法及系统

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130176822A1 (en) * 2010-09-20 2013-07-11 Total Sa Method for estimating elastic parameters by inverting 4d seismic measurements
CN102830423A (zh) * 2012-08-30 2012-12-19 中国石油天然气股份有限公司 储层流体检测方法和储层流体检测装置
CN102866426A (zh) * 2012-09-20 2013-01-09 成都晶石石油科技有限公司 一种利用avo大角度道集分析岩体油气信息的方法
CN103439739A (zh) * 2013-04-08 2013-12-11 中国石油集团东方地球物理勘探有限责任公司 地球物理勘探用纵横波匹配方法及匹配装置
CN103454686A (zh) * 2013-08-12 2013-12-18 中国石油天然气股份有限公司 基于地层切片的小尺度沉积相进行储层预测的方法及系统
CN103454678A (zh) * 2013-08-12 2013-12-18 中国石油天然气股份有限公司 一种地震切片等时性的确定方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
窦玉坛,等: "鄂尔多斯盆地SX地区PP+PS波地震流体检测技术", 《天然气地球科学》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105717542A (zh) * 2016-04-11 2016-06-29 北京科胜伟达石油科技股份有限公司 基于三维地震的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
US10984590B1 (en) 2019-12-06 2021-04-20 Chevron U.S.A. Inc. Generation of subsurface representations using layer-space
US11010969B1 (en) 2019-12-06 2021-05-18 Chevron U.S.A. Inc. Generation of subsurface representations using layer-space
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
US11263362B2 (en) 2020-01-16 2022-03-01 Chevron U.S.A. Inc. Correlation of multiple wells using subsurface representation
US11320566B2 (en) 2020-01-16 2022-05-03 Chevron U.S.A. Inc. Multiple well matching within subsurface representation
US11397279B2 (en) 2020-03-27 2022-07-26 Chevron U.S.A. Inc. Comparison of wells using a dissimilarity matrix

Also Published As

Publication number Publication date
CN105372705B (zh) 2017-10-27

Similar Documents

Publication Publication Date Title
CN105372705A (zh) 一种基于多波资料的地层切片方法
EP2247964B1 (en) Method for determining geological information related to a subsurface volume of interest
Behm et al. Application of stacking and inversion techniques to three-dimensional wide-angle reflection and refraction seismic data of the Eastern Alps
CN104678434B (zh) 一种预测储层裂缝发育参数的方法
CN103454685A (zh) 利用测井约束波阻抗反演预测砂体厚度的方法和装置
Kandilarov et al. Crustal structure of the ultra-slow spreading Knipovich Ridge, North Atlantic, along a presumed ridge segment center
CN102313900A (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
Rointan et al. Improvement of seismic velocity model by selective removal of irrelevant velocity variations
Martorana et al. Assessment of a geological model by surface wave analyses
Kammann et al. Deep onshore reflection seismic imaging of the chalk group strata using a 45 kg accelerated weight-drop and combined recording systems with dense receiver spacing
Hanafy et al. New insights on structure and stratigraphic interpretation for assessing the hydrocarbon potentiality of the offshore Nile Delta basin, Egypt
Lin et al. Marine 3D seismic volumes from 2D seismic survey with large streamer feathering
Chandler et al. Utility of the horizontal-to-vertical spectral ratio passive seismic method for estimating thickness of Quaternary sediments in Minnesota and adjacent parts of Wisconsin
Kaiser et al. Detailed images of the shallow Alpine Fault Zone, New Zealand, determined from narrow-azimuth 3D seismic reflection data
Grobys et al. Extensional and magmatic nature of the Campbell Plateau and Great South Basin from deep crustal studies
Toosy et al. Depositional architecture of Eocene and Cretaceous reservoirs (blending seismic and well data) in Mubarik Area, Central Indus Basin, Pakistan
Victor et al. Joint Inversion of High‐Frequency Receiver Functions and Surface‐Wave Dispersion: Case Study in the Parnaíba Basin of Northeast Brazil
Frifita et al. A renewed view of basement structural geometry beneath the Southern Atlas Front in Tunisia inferred from gravity, seismic reflection and earthquake data
Fischer et al. Development of an exploration concept for the Shuaiba Formation using seismic sequence and facies analysis with forward modeling
Colombo et al. Seismic-airborne TEM joint inversion and surface consistent refraction analysis: New technologies for complex near surface corrections
Bown Legacy seismic investigations of karst surfaces: Implications for heavy oil extraction from the Devonian Grosmont formation, northeastern Alberta, Canada
Majdański et al. Near-surface structure of the Carpathian Foredeep marginal zone in the Roztocze Hills area
Jang et al. Application of Reverse Time Migration to Faults Imaging in Rakhine Basin, Myanmar
McBride et al. Geophysical characterization of volcanic layering

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