CN112946737B - 一种利用纵横波速度增量交汇图识别天然气水合物的方法 - Google Patents
一种利用纵横波速度增量交汇图识别天然气水合物的方法 Download PDFInfo
- Publication number
- CN112946737B CN112946737B CN202110073859.5A CN202110073859A CN112946737B CN 112946737 B CN112946737 B CN 112946737B CN 202110073859 A CN202110073859 A CN 202110073859A CN 112946737 B CN112946737 B CN 112946737B
- Authority
- CN
- China
- Prior art keywords
- wave velocity
- natural gas
- longitudinal
- transverse wave
- gas hydrate
- 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
- NMJORVOYSJLJGU-UHFFFAOYSA-N methane clathrate Chemical compound C.C.C.C.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O NMJORVOYSJLJGU-UHFFFAOYSA-N 0.000 title claims abstract description 134
- 238000000034 method Methods 0.000 title claims abstract description 44
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 56
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Natural products C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims abstract description 52
- 239000003345 natural gas Substances 0.000 claims abstract description 37
- -1 natural gas hydrates Chemical class 0.000 claims abstract description 22
- 229920006395 saturated elastomer Polymers 0.000 claims abstract description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 15
- 238000010586 diagram Methods 0.000 claims description 13
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 13
- 239000011707 mineral Substances 0.000 claims description 13
- 239000011148 porous material Substances 0.000 claims description 8
- 239000004927 clay Substances 0.000 claims description 5
- 238000005553 drilling Methods 0.000 abstract description 15
- 238000005755 formation reaction Methods 0.000 description 49
- 239000007789 gas Substances 0.000 description 5
- 230000002159 abnormal effect Effects 0.000 description 4
- 239000007787 solid Substances 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000001747 exhibiting effect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000011800 void material Substances 0.000 description 2
- 241001251371 Betula chinensis Species 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012822 chemical development Methods 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000002803 fossil fuel Substances 0.000 description 1
- 150000004677 hydrates Chemical class 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- NIFIFKQPDTWWGU-UHFFFAOYSA-N pyrite Chemical compound [Fe+2].[S-][S-] NIFIFKQPDTWWGU-UHFFFAOYSA-N 0.000 description 1
- 229910052683 pyrite Inorganic materials 0.000 description 1
- 239000011028 pyrite Substances 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
- 239000000126 substance Substances 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
- 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
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- 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/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- 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
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (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
本发明属于海洋能源勘探技术领域，涉及一种利用纵横波速度增量交汇图识别天然气水合物的方法，S1、获得实际的纵波速度Vp和横波速度Vs；S2、计算含水饱和地层的纵波速度Vp0和横波速度Vs0；S3、计算纵波速度增量ΔVp和横波速度增量ΔVs；S4、确定阈值曲线；S5、画出纵横波速度增量交汇图；S6、判断：是否含天然气水合物地层；S7、确定天然气水合物的孔隙度理论线；S8、根据天然气水合物的孔隙度理论线估算待测天然气水合物的饱和度范围。本发明所述方法提高天然气水合物识别的可信度，准确勘探识别天然气水合物的位置，减少干井率，节省钻井费用；保障海上钻井安全。
Description
技术领域
本发明属于海洋能源勘探技术领域，涉及一种利用纵横波速度增量交汇图识别天然气水合物的方法，尤其涉及一种利用纵横波速度增量交汇图定性识别天然气水合物和定量估计天然气水合物饱和度范围的方法。
背景技术
天然气水合物是由水分子和天然气分子在特定的温度和压强等条件下生成的类似冰状的固体化合物，具有可燃性，燃烧能量高，又称为可燃冰。它主要分布于深海海底以下或永久冻土带之下。由于天然气水合物具有分布范围广泛以及储量巨大等优势，被视为未来化石燃料的替代能源，受到国内外的广泛关注。目前主要以地球物理和地球化学手段探测天然气水合物，随着海洋科技的发展，水下成像技术逐渐被应用，但其成本较高。物探法是目前普遍使用的天然气水合物的勘探方法，有着较高的可信度，尤其是地震勘探技术应用最为广泛(樊栓狮，1999，参见参考文献[1]；李小森，2008，参见参考文献[2]；刘培栋等，2016，参见参考文献[3]；栾锡武等，2008，参见参考文献[5])。
与本发明相关的现有技术一的技术方案简述如下：
利用地震剖面上的似海底反射(BSR)识别天然气水合物：
在天然气水合物存在的地震剖面上，常常会伴随出现一条近似平行于海底的强反射轴，称之为BSR，通常被用来作为天然气水合物识别的标志，BSR具有如下特点：(1)近似平行于海底，与沉积地层斜交；(2)强反射轴，极性与海底反射轴相反；(3)BSR上方呈弱振幅或振幅空白带；(4)BSR上方具有高速特征。假如在地震剖面上出现有上述特征的反射轴，则可以推断其上有天然气水合物的存在(马在田等，2002，参见参考文献[6])。
以上所述现有技术一具有以下缺点：
由于BSR与天然气水合物的存在关系并不明确，因此该方法具有较低的可信度。并非所有的天然气含水合物地层都会出现BSR，对于平缓的海底，即使存在天然气水合物，也可能并不形成BSR，这和天然气水合物与游离气的饱和度以及储存方式等因素相关。BSR通常出现在海域环境的斜坡及海底沉积地层。在没有BSR出现的情况下，如何识别出天然气水合物是天然气水合物勘探面临的一个难点问题。此外，值得注意的是，并非所有BSR出现的地方都会有天然气水合物存在。另外，与常规不同的是，尽管绝大部分天然气水合物层都在BSR上方，但深海钻探资料显示，并非所有的天然气水合物层都位于BSR之上，也可能在其下方，因此，用BSR识别水合物具有一定的难度，且可信度较低(宋海斌等，2001，参见参考文献[7])。
与本发明相关的现有技术二的技术方案简述如下：
利用纵横波速度异常识别天然气水合物：
研究表明，当地层中充填天然气水合物时，地层的纵横波速度升高。因此，有学者提出用纵横波速度升高的特征识别天然气水合物。换言之，如果地层的纵横波速度比其他地层背景速度高，则推断该处存在天然气水合物。
以上所述现有技术二具有以下缺点：
由于当地层中充填其它矿物时，如，黄铁矿地层的纵横波速度也会升高，因此非天然气水合物因素同样会导致纵横波速度异常升高，从而无法指示是否是由于天然气水合物而引起的。因此，利用纵横波速度异常升高识别天然气水合物的方法也存在缺陷。
与本发明相关的现有技术三的技术方案简述如下：
一种利用纵横波波阻抗增量比识别天然气水合物的方法(刘学伟，CN200910089934.6，即参考文献[4])：
该发明针对现有技术二中，地层的纵横波高速度异常识别天然气水合物可信度低的问题，提出用纵横波波阻抗增量比的方法，提高天然气水合物识别的准确性。给出了含天然气水合物地层与非含天然气水合物纵横波波阻抗增量比之间的阈值，当纵横波波阻抗增量比大于阈值，则表示为含天然气水合物地层，反之为非含天然气水合物地层。
以上所述现有技术三具有以下缺点：
由于该方法需要反演地震横波信息，横波信息的准确性对后续计算纵横波波阻抗增量比至关重要，横波波阻抗作为纵横波波阻抗增量比比值中的分母，如果横波信息出现误差，那么会导致计算的比值出现较大误差，影响最终天然气水合物识别的结果。
本发明利用纵横波速度增量交汇图识别判断地层的高速特征是由含天然气水合物地层还是非含天然气水合物地层而引起的，从而提高天然气水合物识别的可信度，并估算天然气水合物饱和度的范围。
关键术语定义列表
(1)纵波速度增量：ΔVp＝Vp-Vp0；
其中，Vp：地层的实测纵波速度。如果由于某种地质因素(如孔隙空间中充填矿物或天然气水合物)使得地层纵波速度增加，则VP表现为高速特征。
Vp0：计算的含水饱和地层(孔隙空间充填水时的地层)的纵波速度。即假设引起地层速度增加的因素不存在，且孔隙空间完全充填水时，地层的纵波速度。
(2)横波速度增量：ΔVs＝Vs-Vs0；
其中，Vs：地层的实测横波速度。如果由于某种地质因素(如孔隙空间中充填矿物或天然气水合物)使得地层横波速度增加，则Vs表现为高速特征。
Vs0：计算的含水饱和地层(孔隙空间充填水时的地层)的横波速度。即假设引起地层速度增加的因素不存在，且孔隙空间完全充填水时，地层的横波速度。
发明内容
针对上述现有技术存在的缺陷，本申请提供了一种能够提高天然气水合物识别可信度的利用纵横波速度增量交汇图识别天然气水合物的方法。本发明解决了根据纵横波速度异常升高识别天然气水合物可信度低的问题，明确了纵横波速度增量与天然气水合物之间的内在关系，从而为准确识别天然气水合物提供有力支撑，本发明采取的技术方案是：
一种利用纵横波速度增量交汇图识别天然气水合物的方法，包括以下步骤：
S1、利用地震数据反演地层实际的纵波速度Vp和横波速度Vs；
S2、计算出含水饱和地层的纵波速度Vp0和横波速度Vs0；
S3、利用公式(1)计算出纵波速度增量ΔVp，利用公式(2)计算出横波速度增量ΔVs，
ΔVp＝Vp-Vp0 (1)
ΔVs＝Vs-Vs0 (2)
S4、经研究发现，当地层孔隙度为30％时，以粘土的弹性参数为参考，设置临界矿物的体积模量为20.9Gpa，剪切模量为6.5Gpa，密度为2.5g/cm3，计算临界矿物饱和度从0到50％变化的纵横波速度增量交点曲线，可以作为阈值曲线；
S5、画出地层纵横波速度增量的交汇图以及阈值曲线；
S6、进行以下判断：
当纵横波速度增量交点在阈值曲线的左侧，则表示为含天然气水合物地层；当纵横波速度增量交点在阈值曲线的右侧，则表示为非含天然气水合物地层；
S7、在地层孔隙度分别为30％、50％和70％时，计算天然气水合物饱和度从0到50％变化的纵横波速度增量交点曲线，称为孔隙度理论线；
S8、在所述地层纵横波速度增量的交汇图上，将待测天然气水合物地层纵横波速度增量交点的位置与孔隙度理论线上的天然气水合物饱和度位置范围进行对比，估算待测天然气水合物的饱和度范围。
在上述技术方案的基础上，所述纵波速度Vp0为：假设引起地层速度增加的因素不存在，且孔隙空间完全充填水时，地层的纵波速度；所述横波速度Vs0为：假设引起地层速度增加的因素不存在，且孔隙空间完全充填水时，地层的横波速度。
在上述技术方案的基础上，所述纵波速度Vp0和横波速度Vs0通过岩石物理模型计算。
在上述技术方案的基础上，步骤S5所述地层纵横波速度增量的交汇图首先是以横波速度增量为横坐标，以纵波速度增量为纵坐标，形成二维坐标平面；
在所述交汇图上，将步骤S3计算的地层的横波速度增量ΔVs作为纵横波速度增量交点的横坐标；
将步骤S3计算的地层的相应的(即横波速度增量ΔVs与纵波速度增量ΔVp为同一点)纵波速度增量ΔVp作为纵横波速度增量交点的纵坐标；最终形成纵横波速度增量交点，所述交点的横坐标为ΔVs，所述交点的纵坐标为ΔVp。
在上述技术方案的基础上，步骤S4-S6所述阈值曲线也位于所述地层纵横波速度增量的交汇图上。
在上述技术方案的基础上，步骤S8的具体步骤为：在所述孔隙度理论线上标出天然气水合物饱和度值对应的纵横波速度增量交点，根据待测天然气水合物地层纵横波速度增量交点所处的位置，将其位置所对应的天然气水合物饱和度值的范围作为估算待测天然气水合物的饱和度范围。
本发明的有益技术效果如下：
本发明针对当前天然气水合物识别方法可信度低的问题，提出了应用纵横波速度增量交汇图的方法区分天然气水合物和非天然气水合物因素引起的地层高速异常，从而提高天然气水合物识别的可信度，简述如下：
(1)海上钻探费用高昂，准确勘探识别天然气水合物的位置，提高天然气水合物识别可信度，减少干井率，可以节省钻井费用；
(2)天然气水合物一般赋存在海底以下较浅位置的地层中，常规深水油气钻探可能会穿过含有天然气水合物的地层；假设地层中含有天然气水合物，如果勘探不当就会破坏天然气水合物稳定的温压(温度和压力)条件，天然气水合物分解，引发海底塌陷；严重情况下，还会导致钻井平台倾斜或倒塌。为了防范海上钻井事故的发生，应该提前勘探海底天然气水合物的分布情况。提高天然气水合物的识别可信度，保障海上钻井安全。
附图说明
本发明有如下附图：
图1(a)是1245E井位计算的与实测的纵波速度示意图；
图1(b)是1245E井位计算的与实测的横波速度示意图；
图2是1245E井位纵横波速度增量交点识别天然气水合物示意图
图3是1245E井位纵横波速度增量交点估算天然气水合物饱和度的示意图；
图4(a)是1247B井位计算的与实测的纵波速度示意图；
图4(b)是1247B井位计算的与实测的横波速度示意图；
图5是1247B井位纵横波速度增量交点识别天然气水合物示意图；
图6是1247B井位纵横波速度增量交点估算天然气水合物饱和度的示意图；
图7(a)是1250F井位计算的与实测的纵波速度示意图；
图7(b)是1250F井位计算的与实测的横波速度示意图；
图8是1250F井位纵横波速度增量交点识别天然气水合物示意图；
图9是1250F井位纵横波速度增量交点估算天然气水合物饱和度的示意图。
具体实施方式
下面结合附图和具体实施例对本发明做进一步说明，但不作为对本发明的限定。
实施例1
大洋钻探204航次在俄勒冈水合物海岭1245E井位钻井，发现有天然气水合物存在。本例利用1245E井位已知的纵横波速度测井资料，计算该井位位置处的含水饱和地层的纵横波速度，求出纵横波速度增量，根据已确定的地层参数计算阈值曲线，判断地层是否含有天然气水合物，与实际资料天然气水合物位置进行对比。
利用地层纵横波速度增量交汇图识别天然气水合物以及估算天然气水合物饱和度范围的方法，包括如下步骤：
(1)用岩石物理模型计算的含水饱和地层纵横波速度(Vp0，Vs0)，简称计算速度，和实测的地层纵横波速度(Vp，Vs)，简称实测速度，曲线如图1(a)和图1(b)所示；
(2)根据公式(1)计算出纵波速度增量ΔVp，根据公式(2)计算出横波速度增量ΔVs，
ΔVp＝Vp-Vp0 (1)
ΔVs＝Vs-Vs0 (2)
(3)当地层孔隙度为30％时，以粘土的弹性参数为参考，设置临界矿物的体积模量为20.9Gpa，剪切模量为6.5Gpa，密度为2.5g/cm3，计算临界矿物饱和度从0到50％变化的纵横波速度增量交点曲线作为阈值曲线；
(4)如图2所示，以横波速度增量为横坐标，以纵波速度增量为纵坐标，画出地层的纵横波速度增量交点以及阈值曲线(在图中用实线表示，简称阈值)；根据实际资料，在所述纵横波速度增量交点中，地层中含天然气水合物的点(在图中简称水合物点)用方形表示，非含天然气水合物的点(在图中简称非水合物点)用圆形表示，阈值曲线用实线表示；
(5)判断：
当地层的纵横波速度增量交点位于阈值曲线左侧时，表示为含天然气水合物点；否则为非含天然气水合物点。从图2中可见，在阈值曲线左侧主要分布着含天然气水合物点，右侧为非含天然气水合物点，与实际情况基本保持一致。
(6)在地层孔隙度分别为30％、50％和70％时，计算天然气水合物饱和度从0到50％变化的纵横波速度增量交点曲线，称为孔隙度理论线，如图3所示，分别用实线、虚线和点划线表示地层孔隙度为30％、50％和70％，其中，星号、加号和叉号分别表示天然气水合物饱和度为10％、20％和30％的点。图3仍为地层纵横波速度增量交点图，其中含天然气水合物点用方形表示，非含天然气水合物点用圆形表示。从图3中可见，大部分的含天然气水合物点分布在饱和度为0-10％范围内，与实际钻探资料中的结果基本一致。
实施例2
大洋钻探204航次在俄勒冈水合物海岭1247B井位钻井，发现有天然气水合物存在。本例利用1247B井位已知的纵横波速度测井资料，计算该井位位置处的含天然气水饱和地层的纵横波速度，求出纵横波速度增量，根据已确定的地层参数计算阈值曲线，判断地层是否含有天然气水合物，与实际资料天然气水合物位置进行对比。
利用地层纵横波速度增量交汇图识别天然气水合物以及估算天然气水合物饱和度范围的方法，包括如下步骤：
(1)用岩石物理模型计算的含水饱和地层纵横波速度(Vp0，Vs0)和实测的地层纵横波速度(Vp，Vs)曲线如图4(a)和图4(b)所示；
(2)根据公式(1)计算出纵波速度增量ΔVp，根据公式(2)计算出横波速度增量ΔVs，
ΔVp＝Vp-Vp0 (1)
ΔVs＝Vs-Vs0 (2)
(3)当地层孔隙度为30％时，以粘土的弹性参数为参考，设置临界矿物的体积模量为20.9Gpa，剪切模量为6.5Gpa，密度为2.5g/cm3，计算临界矿物饱和度从0到50％变化的纵横波速度增量交点曲线作为阈值曲线；
(4)如图5所示，以横波速度增量为横坐标，以纵波速度增量为纵坐标，画出地层的纵横波速度增量交点以及阈值曲线(在图中用实线表示，简称阈值)；根据实际资料，在所述纵横波速度增量交点中，地层中含天然气水合物的点(在图中简称水合物点)用方形表示，非含天然气水合物的点(在图中简称非水合物点)用圆形表示，阈值曲线用实线表示；
(5)判断：
当地层的纵横波速度增量交点位于阈值曲线左侧时，表示为含天然气水合物点；否则为非含天然气水合物点。从图5中可见，在阈值曲线左侧主要分布着含天然气水合物点，右侧为非含天然气水合物点，与实际情况基本保持一致。
(6)在地层孔隙度分别为30％、50％和70％时，计算天然气水合物饱和度从0到50％变化的纵横波速度增量交点曲线，称为孔隙度理论线，如图6所示，分别用实线、虚线和点划线表示地层孔隙度为30％、50％和70％，其中，星号、加号和叉号分别表示天然气水合物饱和度为10％、20％和30％的点。图6仍为地层纵横波速度增量交点图，其中含天然气水合物点用方形表示，非含天然气水合物点用圆形表示。从图6中可见，大部分的含天然气水合物点分布在饱和度为0-20％范围内，与实际钻探资料中的结果基本一致。
实施例3
大洋钻探204航次在俄勒冈水合物海岭1250F井位钻井，发现有天然气水合物存在。本例利用1250F井位已知的纵横波速度测井资料，计算该井位位置处的含天然气水饱和地层的纵横波速度，求出纵横波速度增量，根据已确定的地层参数计算阈值曲线，判断地层是否含有天然气水合物，与实际资料天然气水合物位置进行对比。
利用地层纵横波速度增量交汇图识别天然气水合物以及估算天然气水合物饱和度范围的方法，包括如下步骤：
(1)用岩石物理模型计算的含水饱和地层纵横波速度(Vp0，Vs0)和实测的地层纵横波速度(Vp，Vs)曲线如图7(a)和图7(b)所示；
(2)根据公式(1)计算出纵波速度增量ΔVp，根据公式(2)计算出横波速度增量ΔVs，
ΔVp＝Vp-Vp0 (1)
ΔVs＝Vs-Vs0 (2)
(3)当地层孔隙度为30％时，以粘土的弹性参数为参考，设置临界矿物的体积模量为20.9Gpa，剪切模量为6.5Gpa，密度为2.5g/cm3，计算临界矿物饱和度从0到50％变化的纵横波速度增量交点曲线作为阈值曲线；
(4)如图8所示，以横波速度增量为横坐标，以纵波速度增量为纵坐标，画出地层的纵横波速度增量交点以及阈值曲线(在图中用实线表示，简称阈值)；根据实际资料，在所述纵横波速度增量交点中，地层中含天然气水合物的点(在图中简称水合物点)用方形表示，非含天然气水合物的点(在图中简称非水合物点)用圆形表示，阈值曲线用实线表示；
(5)判断：
当地层的纵横波速度增量交点位于阈值曲线左侧时，表示为含天然气水合物点；否则为非含天然气水合物点。从图8中可见，在阈值曲线左侧主要分布着含天然气水合物点，右侧为非含天然气水合物点，与实际情况基本保持一致。
(6)在地层孔隙度分别为30％、50％和70％时，计算天然气水合物饱和度从0到50％变化的纵横波速度增量交点曲线，称为孔隙度理论线，如图9所示，分别用实线、虚线和点划线表示地层孔隙度为30％、50％和70％，其中，星号、加号和叉号分别表示天然气水合物饱和度为10％、20％和30％的点。图9仍为地层纵横波速度增量交点图，其中含天然气水合物点用方形表示，非含天然气水合物点用圆形表示。从图9中可见，大部分的含天然气水合物点分布在饱和度为0-20％范围内，与实际钻探资料中的结果基本一致。
本发明的技术关键点和欲保护点如下所述：
技术关键点如下：
1、发现了含天然气水合物地层的纵波速度增量与横波速度增量与非含天然气水合物地层的纵横波速度增量之间的差异。基于这种差异创建了识别天然气水合物的构架；
2、创建了区分含天然气水合物地层与非含天然气水合物地层的纵横波速度增量的阈值方法；
3、发明了基于纵波速度增量与横波速度增量识别天然气水合物的交汇图法；
4、发明了估计天然气水合物饱和度分布范围的方法。
欲保护点如下：
1.利用含天然气水合物地层与非含天然气水合物地层的纵横波速度增量差异识别天然气水合物；
2.提出利用地层纵横波速度增量交汇图方法识别天然气水合物；
3.提出区分含天然气水合物地层与不含天然气水合物地层的纵横波速度增量阈值方法；
4.利用含天然气水合物地层纵横波速度增量估计天然气水合物饱和度的分布范围。
参考文献(如专利/论文/标准)如下所列：
[1]樊栓狮.笼型水合物研究进展[J].化工进展,1999,18(1):5-7.
[2]李小森.天然气水合物能源的勘探与开发[J].现代化工，2008，28(6)：1-13.
[3]刘培栋，王锡文，袁洪仁，邓国振，王勇，薛超.厄瓜多尔西海岸天然气水合物地震识别方法及其分布规律[J].地球物理学报，2016，31(4),1633-1638.
[4]刘学伟.一种利用纵横波波阻抗增量比识别天然气水合物方法.中国，CN200910089934.6
[5]栾锡武，赵克斌，孙冬胜，等.海域天然气水合物勘测的地球物理方法[J].地球物理学进展，2008，23(1)：210-219.
[6]马在田，耿建华，董良国，宋海斌.海洋天然气水合物的地震识别方法研究[J].海洋地质与第四纪地质，2002，22(I):1-8
[7]宋海斌，松林修，杨胜雄，江为为.海洋天然气水合物的地球物理研究(Ⅱ)地震方法[J].地球物理学进展，2001,16(3):110-118.
显然，本发明的上述实施例仅仅是为清楚地说明本发明所做的举例，而并非是对本发明实施方式的限定，对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动，这里无法对所有的实施方式予以穷举，凡是属于本发明的技术方案所引申出的显而易见的变化或变动仍处于本发明的保护范围之列。
本说明书中未做详细描述的内容属于本领域专业技术人员公知的现有技术。
Claims (6)
1.一种利用纵横波速度增量交汇图识别天然气水合物的方法，其特征在于，包括以下步骤：
S1、利用地震数据反演地层实际的纵波速度Vp和横波速度Vs；
S2、计算出含水饱和地层的纵波速度Vp0和横波速度Vs0；
S3、利用公式(1)计算出纵波速度增量ΔVp，利用公式(2)计算出横波速度增量ΔVs，
ΔVp＝Vp-Vp0 (1)
ΔVs＝Vs-Vs0 (2)
S4、当地层孔隙度为30％时，以粘土的弹性参数为参考，设置临界矿物的体积模量为20.9Gpa，剪切模量为6.5Gpa，密度为2.5g/cm3，计算临界矿物饱和度从0到50％变化的纵横波速度增量交点曲线，作为阈值曲线；
S5、画出地层纵横波速度增量的交汇图以及阈值曲线；
S6、进行以下判断：
当纵横波速度增量交点在阈值曲线的左侧，则表示为含天然气水合物地层；当纵横波速度增量交点在阈值曲线的右侧，则表示为非含天然气水合物地层；
S7、在地层孔隙度分别为30％、50％和70％时，计算天然气水合物饱和度从0到50％变化的纵横波速度增量交点曲线，称为孔隙度理论线；
S8、在所述地层纵横波速度增量的交汇图上，将待测天然气水合物地层纵横波速度增量交点的位置与孔隙度理论线上的天然气水合物饱和度位置范围进行对比，估算待测天然气水合物的饱和度范围。
2.如权利要求1所述的利用纵横波速度增量交汇图识别天然气水合物的方法，其特征在于：
所述纵波速度Vp0为：假设引起地层速度增加的因素不存在，且孔隙空间完全充填水时，地层的纵波速度；所述横波速度Vs0为：假设引起地层速度增加的因素不存在，且孔隙空间完全充填水时，地层的横波速度。
3.如权利要求1所述的利用纵横波速度增量交汇图识别天然气水合物的方法，其特征在于：
所述纵波速度Vp0和横波速度Vs0通过岩石物理模型计算。
4.如权利要求1所述的利用纵横波速度增量交汇图识别天然气水合物的方法，其特征在于：
步骤S5所述地层纵横波速度增量的交汇图首先是以横波速度增量为横坐标，以纵波速度增量为纵坐标，形成二维坐标平面；
在所述交汇图上，将步骤S3计算的地层的横波速度增量ΔVs作为纵横波速度增量交点的横坐标；
将步骤S3计算的地层的相应的纵波速度增量ΔVp作为纵横波速度增量交点的纵坐标；最终形成纵横波速度增量交点，所述交点的横坐标为ΔVs，所述交点的纵坐标为ΔVp。
5.如权利要求1所述的利用纵横波速度增量交汇图识别天然气水合物的方法，其特征在于：
步骤S4-S6所述阈值曲线也位于所述地层纵横波速度增量的交汇图上。
6.如权利要求1所述的利用纵横波速度增量交汇图识别天然气水合物的方法，其特征在于：
步骤S8的具体步骤为：在所述孔隙度理论线上标出天然气水合物饱和度值对应的纵横波速度增量交点，根据待测天然气水合物地层纵横波速度增量交点所处的位置，将其位置所对应的天然气水合物饱和度值的范围作为估算待测天然气水合物的饱和度范围。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110073859.5A CN112946737B (zh) | 2021-01-20 | 2021-01-20 | 一种利用纵横波速度增量交汇图识别天然气水合物的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110073859.5A CN112946737B (zh) | 2021-01-20 | 2021-01-20 | 一种利用纵横波速度增量交汇图识别天然气水合物的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112946737A CN112946737A (zh) | 2021-06-11 |
| CN112946737B true CN112946737B (zh) | 2023-10-31 |
Family
ID=76235730
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110073859.5A Active CN112946737B (zh) | 2021-01-20 | 2021-01-20 | 一种利用纵横波速度增量交汇图识别天然气水合物的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112946737B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113960664B (zh) * | 2021-10-25 | 2023-07-28 | 中国石油天然气股份有限公司 | 一种纵横波反射率模版估算水合物和游离气饱和度的方法 |
| CN117031545B (zh) * | 2023-07-28 | 2024-02-13 | 中国海洋大学 | 基于敏感参数识别砂质储层水合物与游离气共存层的方法 |
Citations (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101644781A (zh) * | 2009-07-28 | 2010-02-10 | 刘学伟 | 一种利用纵横波波阻抗增量比识别天然气水合物的方法 |
| JP2010530976A (ja) * | 2007-06-21 | 2010-09-16 | シュルムバーガー ホールディングス リミッテッド | マルチアトリビュート地震探査データによるガスハイドレートの特徴付け |
| RU2490676C1 (ru) * | 2012-04-11 | 2013-08-20 | Александр Валентинович Воробьев | Способ обнаружения подводных залежей газовых гидратов |
| WO2014071385A1 (en) * | 2012-11-05 | 2014-05-08 | Fugro Geoconsulting, Inc. | Method of indicating the presence of gas hydrate and shallow gas in deepwater environment |
| CN105911587A (zh) * | 2016-04-22 | 2016-08-31 | 中国地质大学（北京） | 一种利用单程波算子的双程波叠前深度偏移方法 |
| CN109212162A (zh) * | 2018-08-17 | 2019-01-15 | 西南石油大学 | 一种估算成岩天然气水合物储层力学参数的方法 |
| CN109298463A (zh) * | 2018-09-27 | 2019-02-01 | 北京珠玛阳光科技有限公司 | 一种联合纵、横波速度和密度预测天然气水合物方法 |
| CN110579802A (zh) * | 2019-10-09 | 2019-12-17 | 中国科学院海洋研究所 | 一种天然气水合物储层物性参数的高精度反演方法 |
| CN111221039A (zh) * | 2020-01-21 | 2020-06-02 | 中国石油大学(北京) | 天然气水合物的弹性波波速确定方法、装置和设备 |
| CN111722282A (zh) * | 2020-06-18 | 2020-09-29 | 中国科学院海洋研究所 | Avo预测天然气水合物储层顶部水合物饱和度的方法 |
| CN111812711A (zh) * | 2020-07-20 | 2020-10-23 | 中国石油大学(北京) | 基于反射波场的水合物地层地震物理模拟测试设备及方法 |
| CN111856560A (zh) * | 2020-07-08 | 2020-10-30 | 中国海洋大学 | 一种天然气水合物储层信息评估方法及其应用 |
| CN112133377A (zh) * | 2020-08-28 | 2020-12-25 | 中国石油天然气集团有限公司 | 一种海域天然气水合物赋存状态的判别方法及系统 |
| CN112149282A (zh) * | 2020-08-28 | 2020-12-29 | 中国石油天然气集团有限公司 | 一种井中天然气水合物饱和度岩石物理计算方法及系统 |
| CN112182840A (zh) * | 2020-08-28 | 2021-01-05 | 中国石油天然气集团有限公司 | 一种天然气水合物地层衰减建模方法与装置 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2014082018A1 (en) * | 2012-11-23 | 2014-05-30 | Fugro Geoconsulting, Inc. | Method and system for identification of gas hydrates and fee gas in geologic beds |
-
2021
- 2021-01-20 CN CN202110073859.5A patent/CN112946737B/zh active Active
Patent Citations (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2010530976A (ja) * | 2007-06-21 | 2010-09-16 | シュルムバーガー ホールディングス リミッテッド | マルチアトリビュート地震探査データによるガスハイドレートの特徴付け |
| CN101644781A (zh) * | 2009-07-28 | 2010-02-10 | 刘学伟 | 一种利用纵横波波阻抗增量比识别天然气水合物的方法 |
| RU2490676C1 (ru) * | 2012-04-11 | 2013-08-20 | Александр Валентинович Воробьев | Способ обнаружения подводных залежей газовых гидратов |
| WO2014071385A1 (en) * | 2012-11-05 | 2014-05-08 | Fugro Geoconsulting, Inc. | Method of indicating the presence of gas hydrate and shallow gas in deepwater environment |
| CN105911587A (zh) * | 2016-04-22 | 2016-08-31 | 中国地质大学（北京） | 一种利用单程波算子的双程波叠前深度偏移方法 |
| CN109212162A (zh) * | 2018-08-17 | 2019-01-15 | 西南石油大学 | 一种估算成岩天然气水合物储层力学参数的方法 |
| CN109298463A (zh) * | 2018-09-27 | 2019-02-01 | 北京珠玛阳光科技有限公司 | 一种联合纵、横波速度和密度预测天然气水合物方法 |
| CN110579802A (zh) * | 2019-10-09 | 2019-12-17 | 中国科学院海洋研究所 | 一种天然气水合物储层物性参数的高精度反演方法 |
| CN111221039A (zh) * | 2020-01-21 | 2020-06-02 | 中国石油大学(北京) | 天然气水合物的弹性波波速确定方法、装置和设备 |
| CN111722282A (zh) * | 2020-06-18 | 2020-09-29 | 中国科学院海洋研究所 | Avo预测天然气水合物储层顶部水合物饱和度的方法 |
| CN111856560A (zh) * | 2020-07-08 | 2020-10-30 | 中国海洋大学 | 一种天然气水合物储层信息评估方法及其应用 |
| CN111812711A (zh) * | 2020-07-20 | 2020-10-23 | 中国石油大学(北京) | 基于反射波场的水合物地层地震物理模拟测试设备及方法 |
| CN112133377A (zh) * | 2020-08-28 | 2020-12-25 | 中国石油天然气集团有限公司 | 一种海域天然气水合物赋存状态的判别方法及系统 |
| CN112149282A (zh) * | 2020-08-28 | 2020-12-29 | 中国石油天然气集团有限公司 | 一种井中天然气水合物饱和度岩石物理计算方法及系统 |
| CN112182840A (zh) * | 2020-08-28 | 2021-01-05 | 中国石油天然气集团有限公司 | 一种天然气水合物地层衰减建模方法与装置 |
Non-Patent Citations (8)
| Title |
|---|
| A new approach for the identification of gas hydrate in marine sediments;Tian, D., & Liu, X;《Marine Geophysical Research》;第41卷(第3期);1-4 * |
| Gas hydrate identified in sand-rich inferred sedimentary section using downhole logging and seismic data in Shenhu area, South China Sea;Wang, X., Lee, M., Collett, T., Yang, S., Guo, Y., & Wu, S;《Marine and Petroleum Geology》(第51期);298-306 * |
| Joint analysis of P-wave velocity and resistivity for morphology identification and quantification of gas hydrate;Liu, T., Liu, X., & Zhu, T;《Marine and Petroleum Geology》(第112期);1-5 * |
| Natural gas hydrates - a review;COLLETT T S;《Browse Collections》(第89期);146-219 * |
| 天然气水合物和游离气层的弹性参数异常特征及其在饱和度估算中的应用;冯凯;李敏锋;刘竣川;刘学伟;;《石油地球物理勘探》;第43卷(第04期);431-434 * |
| 天然气水合物和游离气饱和度估算的影响因素;王秀娟;吴时国;刘学伟;;《地球物理学报》(第02期);204-211 * |
| 天然气水合物横波速度等效介质模型预测方法;孟大江;文鹏飞;张如伟;赵斌;李延;;《石油地球物理勘探》(第01期);130-138 * |
| 海洋天然气水合物地震检测技术及其应用;张光学;张明;杨胜雄;雷新华;徐华宁;刘学伟;梁金强;沙志彬;;《海洋地质与第四纪地质》;第31卷(第04期);52-55 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112946737A (zh) | 2021-06-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101644781B (zh) | 一种利用纵横波波阻抗增量比识别天然气水合物的方法 | |
| CN112946737B (zh) | 一种利用纵横波速度增量交汇图识别天然气水合物的方法 | |
| Boehm et al. | Fluidized sandstone intrusions as an indicator of paleostress orientation, Santa Cruz, California | |
| CN111856560B (zh) | 一种天然气水合物储层信息评估方法及其应用 | |
| CN105468886A (zh) | 基于岩石物理学参数计算地层压力的方法 | |
| CN103293553A (zh) | 一种复杂海底上下缆地震采集数据边界元延拓校正方法 | |
| Suter et al. | The Shannon Sandstone and isolated linear sand bodies: interpretations and realizations | |
| CN111487681B (zh) | 天然气水合物及下伏游离气储集体地震响应特征分析方法 | |
| CN111190228B (zh) | 一种海底扇水道储层构型识别方法 | |
| Song et al. | Controls on shallow gas distribution, migration, and associated geohazards in the Yangtze subaqueous delta and the Hangzhou Bay | |
| CN101609166B (zh) | 一种水域区近地表结构的测量方法 | |
| CN113219531A (zh) | 致密砂岩气水分布的识别方法及装置 | |
| CN113835121A (zh) | 一种利用avo属性分析识别孔隙型天然气水合物方法 | |
| CN112965106B (zh) | 利用纵横波速度的奇异性识别天然气水合物赋存类型方法 | |
| Wardhana et al. | Exploring the overlooked seurula potential in Arun field; a preliminary study from outcrops and subsurface data in North Sumatra Basin | |
| Nainggolan et al. | Depositional sequence interpretation using seismic and well data of offshore Central Sumatra Basin | |
| CN113740911B (zh) | 一种基于坐标旋转波阻抗反演提高储层预测精度的方法 | |
| CN117408186B (zh) | 地震剖面中天然气水合物似海底反射层的判断方法 | |
| Guo et al. | Quantitative engraving technology and application of ordovician carbonate reservoir in Tarim Basin | |
| CN103135130B (zh) | 一种海洋阶梯形电缆地震接收方法 | |
| Zhou et al. | Evaluation of CO 2 hydrate storage potential in the Qiongdongnan Basin via combining the phase equilibrium mechanism and the volumetric method. | |
| Zhao et al. | Application of seismic multi-attribute technique in fracture prediction of carbonate reservoir | |
| Meckel et al. | Offshore CCS in the Gulf of Mexico, with Emphasis on the Inner-Shelf Cenozoic Stratigraphy of Texas, USA | |
| Yu et al. | Fine Depiction of the Single Sand Body and Connectivity Unit of a Deltaic Front Underwater Distributary Channel: Taking the Third Member of the Dongying Formation in the Cha71 Fault Block of the Chaheji Oilfield as an Example | |
| Koehn et al. | Developing a representative model for offshore CO2 storage in the Gulf of Mexico continental shelf |
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 |