CN106680872A - 一种高效的海洋天然气水合物地震研究方法 - Google Patents

一种高效的海洋天然气水合物地震研究方法 Download PDF

Info

Publication number
CN106680872A
CN106680872A CN201710066004.3A CN201710066004A CN106680872A CN 106680872 A CN106680872 A CN 106680872A CN 201710066004 A CN201710066004 A CN 201710066004A CN 106680872 A CN106680872 A CN 106680872A
Authority
CN
China
Prior art keywords
bsr
hydrate
polarity
seismic
amplitude
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
CN201710066004.3A
Other languages
English (en)
Other versions
CN106680872B (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 University of Petroleum East China
Original Assignee
China University of Petroleum East China
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 University of Petroleum East China filed Critical China University of Petroleum East China
Priority to CN201710066004.3A priority Critical patent/CN106680872B/zh
Publication of CN106680872A publication Critical patent/CN106680872A/zh
Application granted granted Critical
Publication of CN106680872B publication Critical patent/CN106680872B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
    • 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
    • G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Oceanography (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种高效的海洋天然气水合物地震研究方法，属于天然气水合物地震勘探技术领域。该方法可以精确判断似海底反射层及下伏强振幅反射区是否与水合物及下伏游离气区相关，并提供了一种高效率刻画水合物底界BSR在不同深度分布特征的方法，并通过提取BSR到下部合适深度时窗的属性剖面，得到下伏游离气区的分布特征。本方法能够用于准确识别与水合物相关的BSR，并得到BSR与下伏游离气区的平面分布范围，对于认识水合物体系具有十分重要的作用。

Description

一种高效的海洋天然气水合物地震研究方法
技术领域
本发明属于天然气水合物勘探技术领域，涉及一种高效的海洋天然气水合物地震研究方法，具体地说，涉及一种高效的利用地震数据研究海洋天然气水合物体系分布特征的方法。
背景技术
天然气水合物俗称‘可燃冰’，是洁净的未来新能源。水合物的稳定性对环境条件有很高要求，一般存在于高压(>10MPa)和低温(0-10℃)条件下。当温压条件发生改变时，水合物会快速分解，释放出大量气体和水。除蕴含巨大的资源量外，受其稳定性影响，水合物还可能引起引发地质灾害如海底滑坡、毁坏海上作业机器甚至引发海啸。因此，有必要对水合物进行深入研究。
地震勘探是最常用的天然气水合物勘探方法，用地震反射数据进行研究是水合物先期研究的重点工作，主要是利用似海底反射层BSR作为识别水合物的地震标志。但不是所有的BSR都与水合物有关。与成岩作用有关的BSR，也具有类似的地震反射特征，强反射、斜切等时地层，与海底大致平行。目前对BSR的地震研究一般很少重视除振幅外的其他属性。
另外，目前对BSR的地震研究一般是通过解释二维地震剖面和研究与BSR相交的等时地层的平面图来分析BSR的分布特征，这种方法不能直观地反应BSR在三维空间的分布特征，且作图任务繁重，需要解释大量与BSR相交的等时地层反射层，制作不同反射层的平面图来得到BSR在不同深度的分布范围。
发明内容
本发明的目的在于克服上述技术存在的缺陷，提供一种高效的海洋天然气水合物地震研究方法，是一种确定BSR及下伏强振幅带是否与水合物及下伏游离气有关、高效率精细刻画水合物底界BSR及下伏游离气区FGZs的空间展布特征的方法。
其具体技术方案为：
一种高效的海洋水合物地震研究方法，包括以下步骤：
步骤①：在工区主测线和联络测线方向分别选取具有BSR反射特征的地震剖面各两条；
步骤②：提取这四条地震剖面对应的视极性属性体；
步骤③：判断BSR是否与水合物有关：若BSR为正极性则说明其与水合物无关，结束水合物研究；若BSR为负极性，则说明其与水合物有关，可以进行水合物研究，继续进行下面④-⑦步骤；
步骤④：提取地震数据体的视极性属性体；
步骤⑤：BSR代表水合物底界，为水合物和下伏游离气的相界面，在分辨率较高的地震数据上，BSR会表现出不连续的阶梯状、叠瓦状特征，大网格地震解释后进行空间插值会造成解释结果的偏差。因此，需要在工区内对BSR反射层进行小网格2*2精细地震解释、追踪，在解释网格的基础上进行线性插值，制作BSR的时间构造平面图；线性插值即曲线y＝f(x)上两点P1(x1,y1)和P2(x2,y2)之间用直线y＝ax+b代替。
步骤⑥：在BSR时间构造平面图的基础上制作包含信息量远大于常规等时地层平面图的BSR的均方根(Root-Mean-Square：RMS)振幅属性平面图，时窗范围选择BSR上下各5ms，高效率刻画水合物底界在不同深度的分布特征；
RMS振幅是指在某一点将指定时窗内的振幅值的平方的平均值再开平方，对较大的振幅非常敏感，可用于识别亮点、暗点，指示烃类等。
式中，ai为瞬时振幅。
步骤⑦：利用视极性属性体结合地震剖面判断BSR下伏强振幅带是否为游离气FGZs，负极性说明为游离气FGZs，统计BSR反射层下部的游离气FGZs分布的最大时窗范围Tmax；
步骤⑧：通过提取BSR反射层到下部Tmax时窗范围内的RMS振幅属性平面图，得到下伏游离气区(Free Gas Zones:FGZs)的分布特征。
进一步，步骤①中所述BSR反射特征包括斜切等时地层、强振幅、近似平行于海底。
进一步，所述步骤②③通过提取地震数据的视极性属性，来辅助判断剖面上的BSR是否与水合物有关。成岩作用也可能类似BSR的反射轴发生，在地震剖面上，成岩作用导致的BSR也表现为斜切等时地层、强振幅、与海底大致平行等特征。通过提取不同地震属性进行对比分析，发现视极性属性能够有效判断BSR类型。视极性属性是实地震道的符号位，假设零相位子波、视极性与反射系数的极性相同，通常用于振幅异常的品质分析。水合物相关的BSR反射层具有负视极性，而成岩作用相关的BSR反射层具有正视极性。
进一步，步骤⑧所制作的反映FGZs分布特征的RMS振幅属性图中，并不是所有强振幅反射区都是FGZs的分布范围，需要结合视极性属性体进一步判断，RMS振幅图上具有负视极性的强振幅区才是FGZs的分布范围。
与现有技术相比，本发明的有益效果为：
(1)由于不是所有的BSR及下伏强振幅带都与水合物和下伏游离气区FGZs有关。与成岩作用有关的BSR，也具有类似的地震反射特征。通过优选视极性属性剖面，结合地震剖面特征进行判断，若BSR和下伏强振幅带表现为强的负视极性，则能够确定BSR与水合物相关，且下部强振幅带为FGZs；
(2)制作BSR的地震属性平面图来代替大量的等时地层反射层平面图，提供了一种高效率刻画水合物底界在不同深度分布特征的方法，能够反映BSR反射层的空间展布特征。而常规等时地震层位的属性平面图只能反映BSR在某一深度的特征，因此需要追踪、解释大量的地震等时层位来了解BSR在不同深度的平面分布特征；
(3)通过选取BSR下部合适时窗范围，制作反映FGZs分布范围的RMS振幅属性平面图。
附图说明
图1为本发明高效的海洋水合物地震研究方法流程图。
图2为BSR的剖面特征，其中图2(a)地震剖面；图2(b)视极性属性剖面。
图3为线性插值示意图。
图4(a)BSR的RMS属性平面图，图4(b)地震剖面显示BSR及下伏强振幅FGZ。
图5为BSR至下部50ms时窗内的RMS属性图，显示FGZs分布特征。
具体实施方式
为了使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解，下面结合附图和具体实例进一步阐述本发明。
如图1、图2、图3、图4、图5所示，一种高效的海洋天然气水合物地震研究方法。
此种高效的海洋天然气水合物地震研究方法步骤如下：
步骤①：在工区主测线和联络测线方向分别选取具有BSR反射特征(斜切等时地层、强振幅、与海底大致平行)的地震剖面各两条；
步骤②：提取这四条地震剖面对应的视极性属性体；
步骤③：判断BSR是否与水合物有关：
若BSR为正极性则说明其与水合物无关，结束水合物研究；
若BSR为负极性，则说明其与水合物有关，可以进行水合物研究，继续进行下面④-⑦步骤；
步骤④：提取地震数据体的视极性属性体；
步骤⑤：在工区内对BSR反射层进行小网格2*2精细地震解释、追踪，在解释网格的基础上进行线性插值，制作BSR的时间构造平面图；
步骤⑥：在BSR时间构造平面图的基础上制作包含信息量远大于常规等时地层平面图的BSR的均方根(Root-Mean-Square：RMS)振幅属性平面图，时窗范围选择BSR上下各5ms，高效率刻画水合物底界在不同深度的分布特征；
步骤⑦：利用视极性属性体结合地震剖面判断BSR下伏强振幅带是否为游离气FGZs，负极性说明为游离气FGZs，统计BSR反射层下部的游离气FGZs分布的最大时窗范围Tmax；
步骤⑧：通过提取BSR反射层到下部Tmax时窗范围内的RMS振幅属性平面图，得到下伏游离气区(Free Gas Zones:FGZs)的分布特征。
按照以上方法得到了西非Mauritania海域的水合物体系特征。图1显示该地区BSR及下伏强振幅带具有负视极性，表明其与水合物及下伏游离气FGZs有关；通过对BSR进行全区追踪、解释，得到了BSR的RMS属性平面图，显示了BSR在底辟构造顶部的条带状分布特征；通过在BSR至下伏50ms的时窗范围内提取RMS属性平面图，得到了BSR下伏游离气FGZs的分布范围，与BSR分布特征相似，FGZ也在底辟构造顶部分布，显示受底辟构造的控制。
以上所述，仅为本发明最佳实施方式，任何熟悉本技术领域的技术人员在本发明披露的技术范围内，可显而易见地得到的技术方案的简单变化或等效替换均落入本发明的保护范围内。

Claims (4)

1.一种高效的海洋水合物地震研究方法，其特征在于，包括以下步骤：
步骤①：在工区主测线和联络测线方向分别选取具有BSR反射特征的地震剖面各两条；
步骤②：提取这四条地震剖面对应的视极性属性体；
步骤③：判断BSR是否与水合物有关：若BSR为正极性则说明其与水合物无关，结束水合物研究；若BSR为负极性，则说明其与水合物有关，则进行水合物研究，继续进行下面④-⑦步骤；
步骤④：提取地震数据体的视极性属性体；
步骤⑤：BSR代表水合物底界，为水合物和下伏游离气的相界面，在分辨率较高的地震数据上，BSR会表现出不连续的阶梯状、叠瓦状特征，大网格地震解释后进行空间插值会造成解释结果的偏差，需要在工区内对BSR反射层进行小网格2*2精细地震解释、追踪，在解释网格的基础上进行线性插值，制作BSR的时间构造平面图；
线性插值即曲线y＝f(x)上两点P1(x1,y1)和P2(x2,y2)之间用直线y＝ax+b代替；
步骤⑥：在BSR时间构造平面图的基础上制作包含信息量远大于常规等时地层平面图的BSR的均方根振幅属性平面图，时窗范围选择BSR上下各5ms，高效率刻画水合物底界在不同深度的分布特征；
RMS振幅是指在某一点将指定时窗内的振幅值的平方的平均值再开平方，对较大的振幅非常敏感，能用于识别亮点、暗点，指示烃类；
式中，ai为瞬时振幅；
步骤⑦：利用视极性属性体结合地震剖面判断BSR下伏强振幅带是否为游离气FGZs，负极性说明为游离气FGZs，统计BSR反射层下部的游离气FGZs分布的最大时窗范围Tmax；
步骤⑧：通过提取BSR反射层到下部Tmax时窗范围内的RMS振幅属性平面图，得到下伏游离气区的分布特征。
2.根据权利要求1所述的高效的海洋水合物地震研究方法，其特征在于，步骤①中所述BSR反射特征包括斜切等时地层、强振幅、近似平行于海底。
3.根据权利要求1所述的高效的海洋水合物地震研究方法，其特征在于，所述步骤②③通过提取地震数据的视极性属性，来辅助判断剖面上的BSR是否与水合物有关；成岩作用也可能类似BSR的反射轴发生，在地震剖面上，成岩作用导致的BSR也表现为斜切等时地层、强振幅、与海底大致平行的特征，通过提取不同地震属性进行对比分析，发现视极性属性能够有效判断BSR类型，视极性属性是实地震道的符号位，假设零相位子波、视极性与反射系数的极性相同，通常用于振幅异常的品质分析，水合物相关的BSR反射层具有负视极性，而成岩作用相关的BSR反射层具有正视极性。
4.根据权利要求1所述的高效的海洋水合物地震研究方法，其特征在于，步骤⑧所制作的反映FGZs分布特征的RMS振幅属性图中，并不是所有强振幅反射区都是FGZs的分布范围，需要结合视极性属性体进一步判断，RMS振幅图上具有负视极性的强振幅区才是FGZs的分布范围。
CN201710066004.3A 2017-02-06 2017-02-06 一种高效的海洋天然气水合物地震研究方法 Active CN106680872B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710066004.3A CN106680872B (zh) 2017-02-06 2017-02-06 一种高效的海洋天然气水合物地震研究方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710066004.3A CN106680872B (zh) 2017-02-06 2017-02-06 一种高效的海洋天然气水合物地震研究方法

Publications (2)

Publication Number Publication Date
CN106680872A true CN106680872A (zh) 2017-05-17
CN106680872B CN106680872B (zh) 2018-11-06

Family

ID=58859352

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710066004.3A Active CN106680872B (zh) 2017-02-06 2017-02-06 一种高效的海洋天然气水合物地震研究方法

Country Status (1)

Country Link
CN (1) CN106680872B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107817522A (zh) * 2017-10-31 2018-03-20 西安恒泰艾普能源发展有限公司 一种海上定性识别水合物顶底界面的确定性方法
CN107861158A (zh) * 2017-11-07 2018-03-30 中国石油大学(华东) 一种确定天然气水合物bsr界面与海底相交位置的研究方法
CN110456024A (zh) * 2019-08-23 2019-11-15 中国石油大学(华东) 一种分析天然气水合物稳定区边界处碳循环过程的方法及系统
CN111781640A (zh) * 2020-06-30 2020-10-16 广州海洋地质调查局 一种基于bsr特征快速判别水合物含量的方法及处理终端
CN112255678A (zh) * 2020-10-21 2021-01-22 广州海洋地质调查局 一种提高天然气水合物矿体分辨率方法及处理终端
CN113341464A (zh) * 2021-06-04 2021-09-03 中国石油大学(北京) 天然气水合物储层的识别方法、装置、设备和存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101644781A (zh) * 2009-07-28 2010-02-10 刘学伟 一种利用纵横波波阻抗增量比识别天然气水合物的方法
CN102200588A (zh) * 2011-03-22 2011-09-28 成都理工大学 地震资料的波形相似性体曲率分析方法
CN103278852A (zh) * 2013-03-07 2013-09-04 林昌荣 利用地震数据体波形结构特征模型预测油气的方法
CN105224704A (zh) * 2014-06-25 2016-01-06 中国科学院海洋研究所 一种海底不稳定性的评价方法
CN106383360A (zh) * 2016-08-16 2017-02-08 广东石油化工学院 一种基于地震数据的海底沉积物温压场的识别方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101644781A (zh) * 2009-07-28 2010-02-10 刘学伟 一种利用纵横波波阻抗增量比识别天然气水合物的方法
CN102200588A (zh) * 2011-03-22 2011-09-28 成都理工大学 地震资料的波形相似性体曲率分析方法
CN103278852A (zh) * 2013-03-07 2013-09-04 林昌荣 利用地震数据体波形结构特征模型预测油气的方法
CN105224704A (zh) * 2014-06-25 2016-01-06 中国科学院海洋研究所 一种海底不稳定性的评价方法
CN106383360A (zh) * 2016-08-16 2017-02-08 广东石油化工学院 一种基于地震数据的海底沉积物温压场的识别方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
刘培栋 等: "厄瓜多尔西海岸天然气水合物地震识别方法及其分布规律", 《地球物理学进展》 *
唐志远 等: "海上Z区天然气水合物地震研究与成藏远景预测", 《西安石油大学学报（自然科学版）》 *
杨金秀 等: "BSR及其下伏游离气区的分布特征与控制因素", 《石油与天然气地质》 *
杨金秀 等: "流体运移对水合物及下伏游离气分布的控制作用", 《中国地球科学联合学术年会2014》 *
焦春波: "珠江口盆地东部天然气水合物地震识别依据", 《辽宁化工》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107817522A (zh) * 2017-10-31 2018-03-20 西安恒泰艾普能源发展有限公司 一种海上定性识别水合物顶底界面的确定性方法
CN107861158A (zh) * 2017-11-07 2018-03-30 中国石油大学(华东) 一种确定天然气水合物bsr界面与海底相交位置的研究方法
CN110456024A (zh) * 2019-08-23 2019-11-15 中国石油大学(华东) 一种分析天然气水合物稳定区边界处碳循环过程的方法及系统
CN111781640A (zh) * 2020-06-30 2020-10-16 广州海洋地质调查局 一种基于bsr特征快速判别水合物含量的方法及处理终端
CN112255678A (zh) * 2020-10-21 2021-01-22 广州海洋地质调查局 一种提高天然气水合物矿体分辨率方法及处理终端
CN112255678B (zh) * 2020-10-21 2022-02-11 广州海洋地质调查局 一种提高天然气水合物矿体分辨率方法及处理终端
CN113341464A (zh) * 2021-06-04 2021-09-03 中国石油大学(北京) 天然气水合物储层的识别方法、装置、设备和存储介质
CN113341464B (zh) * 2021-06-04 2024-01-26 中国石油大学(北京) 天然气水合物储层的识别方法、装置、设备和存储介质

Also Published As

Publication number Publication date
CN106680872B (zh) 2018-11-06

Similar Documents

Publication Publication Date Title
CN106680872A (zh) 一种高效的海洋天然气水合物地震研究方法
Hansen et al. Snøhvit: The history of injecting and storing 1 Mt CO2 in the fluvial Tubåen Fm
Awdal et al. The characteristics of fracture networks in the Shiranish formation of the Bina Bawi Anticline; comparison with the Taq Taq field, zagros, Kurdistan, NE Iraq
Palermo et al. Insights into a new super-giant gas field-sedimentology and reservoir modeling of the Coral Reservoir Complex, Offshore Northern Mozambique
Robertson Reservoir management using 3D seismic data
Demyttenaere et al. Brunei deep water exploration: from sea floor images and shallow seismic analogues to depositional models in a slope turbidite setting
Kim et al. Site characterization and geotechnical aspects on geological storage of CO 2 in Korea
Hollinsworth et al. Geological evaluation of suprasalt carbon storage opportunities in the Silverpit Basin, United Kingdom Southern North Sea
CN114706125A (zh) 一种基于广角反射信息的潜山裂缝储层预测方法及系统
Willems et al. Geology of the Upper Jurassic to Lower Cretaceous geothermal aquifers in the West Netherlands Basin–an overview
Bashari Khuff formation Permian-Triassic carbonate in the Qatar-South Fars arch hydrocarbon province of the Persian Gulf
Meneley Oil and gas fields in the East Coast and Arctic basins of Canada
Toublanc et al. Ekofisk Field: fracture permeability evaluation and implementation in the flow model
Ozkaya Modeling finite-fracture networks in a partially fractured reservoir in the Middle East
Ciftci et al. Outcrop-based three-dimensional modeling of the Tensleep sandstone at Alkali Creek, Bighorn Basin, Wyoming
Garland et al. Appraisal and development of the Taq Taq field, Kurdistan region, Iraq
Wilson et al. Fracture model of the Upper Freeport coal: Marshall County West Virginia pilot ECBMR and CO2 sequestration site
Seggie et al. Awakening of a sleeping giant: Sunrise-Troubadour gas-condensate Field
Smith Modeling natural fracture networks: Establishing the groundwork for flow simulation at Teapot Dome, Wyoming
Pashin et al. Mesozoic Structure and Petroleum Systems in the De Soto Canyon Salt Basin in the Mobile, Pensacola, Destin Dome, and Viosca Knoll Areas of the MAFLA Shelf
Hamilton Approaches to identifying reservoir heterogeneity in barrier/strandplain reservoirs and the opportunities for increased oil recovery: an example from the prolific oil-producing Jackson-Yegua trend, south Texas
Wach et al. Well placement, cost reduction, and increased production using reservoir models based on outcrop, core, well-log, seismic data, and modern analogs: Onshore and offshore Western Trinidad
Pointe et al. 3-D reservoir and stochastic fracture network modeling for enhanced oil recovery, Circle Ridge Phosphoria/Tensleep reservoir, Wind River Reservation, Arapaho and Shoshone Tribes, Wyoming
Duvbiama et al. 3D Static Modelling of an offshore field in the Niger-delta
Lothe et al. Updated estimate of storage capacity and evaluation of Seal for selected Aquifers (D26)

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant