CN115480296A - 天然气水合物稳定带底界预测方法及装置 - Google Patents
天然气水合物稳定带底界预测方法及装置 Download PDFInfo
- Publication number
- CN115480296A CN115480296A CN202110605588.3A CN202110605588A CN115480296A CN 115480296 A CN115480296 A CN 115480296A CN 202110605588 A CN202110605588 A CN 202110605588A CN 115480296 A CN115480296 A CN 115480296A
- Authority
- CN
- China
- Prior art keywords
- work area
- depth
- seabed
- gas hydrate
- natural gas
- 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
- 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 252
- 238000000034 method Methods 0.000 title claims abstract description 50
- 230000015572 biosynthetic process Effects 0.000 claims description 39
- 238000005755 formation reaction Methods 0.000 claims description 39
- 238000004590 computer program Methods 0.000 claims description 14
- 238000012545 processing Methods 0.000 claims description 12
- 239000013535 sea water Substances 0.000 claims description 12
- 238000003860 storage Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 5
- 239000013049 sediment Substances 0.000 claims description 5
- 238000010586 diagram Methods 0.000 description 20
- 230000006870 function Effects 0.000 description 14
- 238000004891 communication Methods 0.000 description 12
- 238000011160 research Methods 0.000 description 11
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Natural products C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 9
- 238000005259 measurement Methods 0.000 description 7
- 239000000872 buffer Substances 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 239000003345 natural gas Substances 0.000 description 3
- -1 natural gas hydrates Chemical class 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 150000004677 hydrates Chemical class 0.000 description 2
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000003381 stabilizer Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
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
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
本发明公开了一种天然气水合物稳定带底界预测方法及装置，该方法包括：确定工区的区域地层深度与天然气水合物稳定带底界的预测温度间的关联关系；所述区域地层深度为海底深度和海底以下地层深度的组合；根据工区的不同海底以下地层深度的实测温度，建立工区的变地温梯度场；针对工区的每一位置，对该位置处不同的海底以下地层深度进行逐一确定，当海底以下地层深度满足：在该海底以下地层深度下，所述关联关系中的天然气水合物稳定带底界的预测温度，与变地温梯度场所表征的实测温度相同，确定该海底以下地层深度为该位置处的天然气水合物稳定带底界的预测深度。本发明可提高天然气水合物稳定带底界预测的适用性和预测精度。
Description
技术领域
本发明涉及天然气水合物勘探技术领域，尤其涉及天然气水合物稳定带底界预测方法及装置。
背景技术
本部分旨在为权利要求书中陈述的本发明实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。
天然气水合物广泛分布于全球大多数大陆边缘的深水地区以及多年冻土地区，它在资源、环境和灾害等方面具有不可忽视的作用。
通常天然气水合物是通过地震剖面上BSR(Barrier Synchronization Register，海底模拟反射层，也称作似海底反射)的识别来预测，但随着海洋天然气水合物的勘探开发的深入，发现水合物有时也会在无BSR的地层中出现，但传统的天然气水合物稳定带底界预测方法，无法在工区无似海底反射的地层中，进行天然气水合物稳定带底界的预测，导致目前天然气水合物稳定带底界预测方法的适用性较低。
同时，常规的天然气水合物稳定带底界预测主要是通过单点地温梯度场或常地温梯度场进行，但用该方法进行大面积三维空间的水合物稳定带底界的预测会出现误差，这主要是因为海底地温梯度在空间上变化的，并不是常数，因此会导致天然气水合物稳定带底界的预测精度降低。
发明内容
本发明实施例提供一种天然气水合物稳定带底界预测方法，用以提高天然气水合物稳定带底界预测的适用性和预测精度，该方法包括：
根据工区的天然气水合物稳定带底界的预测温度，确定工区的区域地层深度与天然气水合物稳定带底界的预测温度间的关联关系；所述天然气水合物稳定带底界的预测温度对应不同的区域地层深度；所述区域地层深度为海底深度与海底以下地层深度之和；
根据工区的不同海底以下地层深度的实测温度，建立工区的变地温梯度场；所述变地温梯度场用于表征在工区的不同位置处，不同的海底以下地层深度的连续地温梯度数据；
针对工区的每一位置，对该位置处不同的海底以下地层深度进行逐一确定，当海底以下地层深度满足：在该海底以下地层深度下，所述关联关系中的天然气水合物稳定带底界的预测温度，与变地温梯度场所表征的实测温度相同，确定该海底以下地层深度为该位置处的天然气水合物稳定带底界的预测深度。
本发明实施例还提供一种天然气水合物稳定带底界预测装置，用以提高天然气水合物稳定带底界预测的适用性和预测精度，该装置包括：
关联关系确定模块，用于根据工区的天然气水合物稳定带底界的预测温度，确定工区的区域地层深度与天然气水合物稳定带底界的预测温度间的关联关系；所述天然气水合物稳定带底界的预测温度对应不同的区域地层深度；所述区域地层深度为海底深度与海底以下地层深度之和；
变地温梯度场建立模块，用于根据工区的不同海底以下地层深度的实测温度，建立工区的变地温梯度场；所述变地温梯度场用于表征在工区的不同位置处，不同的海底以下地层深度的连续地温梯度数据；
预测深度模块，用于针对工区的每一位置，对该位置处不同的海底以下地层深度进行逐一确定，当海底以下地层深度满足：在该海底以下地层深度下，所述关联关系中的天然气水合物稳定带底界的预测温度，与变地温梯度场所表征的实测温度相同，确定该海底以下地层深度为该位置处的天然气水合物稳定带底界的预测深度。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述天然气水合物稳定带底界预测方法。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述天然气水合物稳定带底界预测方法的计算机程序。
本发明实施例中，根据工区的天然气水合物稳定带底界的预测温度，确定工区的区域地层深度与天然气水合物稳定带底界的预测温度间的关联关系；所述天然气水合物稳定带底界的预测温度对应不同的区域地层深度；所述区域地层深度为海底深度和海底以下地层深度的组合；根据工区的不同海底以下地层深度的实测温度，建立工区的变地温梯度场；所述变地温梯度场用于表征在工区的不同位置处，不同的海底以下地层深度的连续地温梯度数据；针对工区的每一位置，对该位置处不同的海底以下地层深度进行逐一确定，当海底以下地层深度满足：在该海底以下地层深度下，所述关联关系中的天然气水合物稳定带底界的预测温度，与变地温梯度场所表征的实测温度相同，确定该海底以下地层深度为该位置处的天然气水合物稳定带底界的预测深度，从而可通过工区的天然气水合物稳定带底界的预测温度和工区的不同海底以下地层深度的实测温度，确定出工区每一位置处的天然气水合物稳定带底界的预测深度，直接实现天然气水合物稳定带底界的预测，不再需要借助似海底反射进行预测，避免了现有技术中需借助似海底反射，导致在工区无似海底反射的地层中无法进行天然气水合物稳定带底界预测的问题，提高了天然气水合物稳定带底界预测的适用性；同时，通过建立工区的变地温梯度场，可准确描述不同海底以下地层深度下的连续地温梯度数据，避免了现有技术下仅能借助单点地温梯度场或常地温梯度场，导致无法准确实现对海底下天然气水合物稳定带底界预测的问题，提高了天然气水合物稳定带底界预测的预测精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本发明实施例中一种天然气水合物稳定带底界预测方法的具体示例图；
图2为本发明实施例中一种天然气水合物稳定带底界预测图的具体示例图；
图3为本发明实施例中一种变地温梯度场的具体示例图；
图4为本发明实施例中一种海底层位时间拾取时长和天然气水合物稳定带底界的地震双程旅行时间预测的示意图；
图5为本发明实施例中一种天然气水合物稳定带底界预测方法的流程示意图；
图6为本发明实施例中一种天然气水合物稳定带底界预测方法的具体示例图；
图7为本发明实施例中一种天然气水合物稳定带底界预测方法的具体示例图；
图8为本发明实施例中一种天然气水合物稳定带底界预测装置的结构示意图；
图9为本发明实施例中一种天然气水合物稳定带底界预测装置的具体示例图；
图10为本发明实施例中用于天然气水合物稳定带底界预测的计算机设备示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
目前，天然气水合物广泛分布于全球大多数大陆边缘的深水地区以及多年冻土地区，它在资源、环境和灾害等方面具有不可忽视的作用。通常天然气水合物是通过地震剖面上BSR的识别来预测，但随着海洋天然气水合物勘探开发的深入，发现水合物有时也会在无BSR的地层中出现，所以亟需解决无BSR时水合物储层空间分布的预测问题。
常规的水合物稳定带底界预测主要是通过单点或常地温梯度场进行，但用该方法进行大面积三维空间的水合物稳定带底界预测会出现误差，这主要是因为海底地温梯度在空间上变化的，并不是常数。
随着海洋天然气水合物勘探开发的深入，发现水合物有时既会在地震剖面中有BSR也会在无BSR的地层中出现。对于有BSR的地层，我们可以用地震数据解释的BSR进行天然气水合物的预测，而对于无BSR地层中的天然气水合物预测问题，亟待需要新的解决方法。
为了提高水合物稳定带底界的计算精度，并解决无BSR时水合物储层空间分布的预测问题，本发明实施例提供一种天然气水合物稳定带底界预测方法，可用以提高天然气水合物稳定带底界预测的适用性和预测精度，可针对海底地温梯度场的可变特点，对实测的离散地温梯度数据进行空间插值，形成地温梯度的平面分布场，进而用于水合物稳定带底界的预测，提高水合物预测精度。
本发明实施例提供的一种天然气水合物稳定带底界预测方法，用以提高天然气水合物稳定带底界预测的适用性和预测精度，如图5所示，该方法包括：
步骤501：根据工区的天然气水合物稳定带底界的预测温度，确定工区的区域地层深度与天然气水合物稳定带底界的预测温度间的关联关系；上述天然气水合物稳定带底界的预测温度对应不同的区域地层深度；上述区域地层深度为海底深度与海底以下地层深度之和；
步骤502：根据工区的不同海底以下地层深度的实测温度，建立工区的变地温梯度场；上述变地温梯度场用于表征在工区的不同位置处，不同的海底以下地层深度的连续地温梯度数据；
步骤503：针对工区的每一位置，对该位置处不同的海底以下地层深度进行逐一确定，当海底以下地层深度满足：在该海底以下地层深度下，上述关联关系中的天然气水合物稳定带底界的预测温度，与变地温梯度场所表征的实测温度相同，确定该海底以下地层深度为该位置处的天然气水合物稳定带底界的预测深度。
本发明实施例中，根据工区的天然气水合物稳定带底界的预测温度，确定工区的区域地层深度与天然气水合物稳定带底界的预测温度间的关联关系；上述天然气水合物稳定带底界的预测温度对应不同的区域地层深度；上述区域地层深度为海底深度和海底以下地层深度的组合；根据工区的不同海底以下地层深度的实测温度，建立工区的变地温梯度场；上述变地温梯度场用于表征在工区的不同位置处，不同的海底以下地层深度的连续地温梯度数据；针对工区的每一位置，对该位置处不同的海底以下地层深度进行逐一确定，当海底以下地层深度满足：在该海底以下地层深度下，上述关联关系中的天然气水合物稳定带底界的预测温度，与变地温梯度场所表征的实测温度相同，确定该海底以下地层深度为该位置处的天然气水合物稳定带底界的预测深度，从而可通过工区的天然气水合物稳定带底界的预测温度和工区的不同海底以下地层深度的实测温度，确定出工区每一位置处的天然气水合物稳定带底界的预测深度，直接实现天然气水合物稳定带底界的预测，不再需要借助似海底反射进行预测，避免了现有技术中需借助似海底反射，导致在工区无似海底反射的地层中无法进行天然气水合物稳定带底界预测的问题，提高了天然气水合物稳定带底界预测的适用性；同时，通过建立工区的变地温梯度场，可准确描述不同海底以下地层深度下的连续地温梯度数据，避免了现有技术下仅能借助单点地温梯度场或常地温梯度场，导致无法准确实现对海底下天然气水合物稳定带底界预测的问题，提高了天然气水合物稳定带底界预测的预测精度。
具体实施时，首先根据工区的天然气水合物稳定带底界的预测温度，确定工区的区域地层深度与天然气水合物稳定带底界的预测温度间的关联关系；上述天然气水合物稳定带底界的预测温度对应不同的区域地层深度；上述区域地层深度为海底深度与海底以下地层深度之和。
实施例中，根据工区的天然气水合物稳定带底界的预测温度，确定工区的区域地层深度与天然气水合物稳定带底界的预测温度间的关联关系，可以包括：
对工区的天然气水合物稳定带底界的预测温度，与不同的区域地层深度，进行指数拟合，得到工区的区域地层深度与天然气水合物稳定带底界的预测温度间的关联关系。
在一个实施例中，按如下公式进行指数拟合，得到工区的区域地层深度与天然气水合物稳定带底界的预测温度间的关联关系：
其中，a和b为拟合参数，无量纲；T1为工区的海底以下地层深度为Dmbsf处的天然气水合物稳定带底界的预测温度，单位为摄氏度；Dmbsf为工区的海底以下地层深度，单位为米；Dsf为工区的海底深度，单位为米。
在上述实施例中，工区的区域地层深度与天然气水合物稳定带底界的预测温度间的关联关系，用于表述区域地层深度与甲烷形成水合物温度之间关系，是根据研究区的地层深度和甲烷形成水合物温度的交会，并通过指数拟合得到的两者间关系式。其中，区域地层深可为海水深度和海底以下地层深度之和，其中，海水深度可称作海底深度。
在上述实施例中，可确定工区的区域地层深度与天然气水合物稳定带底界的预测温度间的关联关系，有助于后续步骤中确定出工区每一位置处的天然气水合物稳定带底界的预测深度，直接实现天然气水合物稳定带底界的预测。
具体实施时，在根据工区的天然气水合物稳定带底界的预测温度，确定工区的区域地层深度与天然气水合物稳定带底界的预测温度间的关联关系后，根据工区的不同海底以下地层深度的实测温度，建立工区的变地温梯度场；上述变地温梯度场用于表征在工区的不同位置处，不同的海底以下地层深度的连续地温梯度数据。
实施例中，根据工区的不同海底以下地层深度的实测温度，建立工区的变地温梯度场，如图6所示，可以包括：
步骤601：针对工区的每一位置，对工区的该位置处不同的海底以下地层深度的实测温度，与不同的海底以下地层深度进行线性拟合，确定对应工区该位置处的不同的海底以下地层深度的实测温度，与海底以下地层深度的线性拟合关系；
步骤602：根据工区每一位置处的上述线性拟合关系，得到工区的离散地温梯度场；上述离散地温梯度场用于表征在工区的不同位置处，不同的海底以下地层深度的离散的地温梯度数据；
步骤603：对工区的离散地温梯度场进行差值处理，建立工区的变地温梯度场。
在一个实施例中，按如下公式确定对应工区该位置处的不同的海底以下地层深度的实测温度，与海底以下地层深度的线性拟合关系：
T2＝Tsf+G×Dmbsf
其中，Tsf为工区该位置处的海底实测温度，单位为摄氏度；G为工区该位置处的地温梯度，单位为摄氏度每米；Dmbsf为工区该位置处的海底以下地层深度，单位为米；T2为工区该位置处的海底以下地层深度为Dmbsf处的实测温度，单位为摄氏度。
在上述实施例中，工区的不同海底以下地层深度的实测温度，可为工区的离散地温梯度数据。上述离散地温梯度数据为工区的海上实际测量数据；上述变地温梯度场为离散地温梯度场经过数据插值后形成的连续的地温梯度场数据。通过建立工区的变地温梯度场，可准确描述不同海底以下地层深度下的连续地温梯度数据，避免了现有技术下仅能借助单点地温梯度场或常地温梯度场，导致无法准确实现对海底下天然气水合物稳定带底界预测的问题，提高了天然气水合物稳定带底界预测的预测精度。
具体实施时，在根据工区的不同海底以下地层深度的实测温度，建立工区的变地温梯度场后，针对工区的每一位置，对该位置处不同的海底以下地层深度进行逐一确定，当海底以下地层深度满足：在该海底以下地层深度下，上述关联关系中的天然气水合物稳定带底界的预测温度，与变地温梯度场所表征的实测温度相同，确定该海底以下地层深度为该位置处的天然气水合物稳定带底界的预测深度。
实施例中，针对工区的每一位置，对该位置处不同的海底以下地层深度进行逐一确定，当海底以下地层深度满足：在该海底以下地层深度下，上述关联关系中的天然气水合物稳定带底界的预测温度，与变地温梯度场所表征的实测温度相同，确定该海底以下地层深度为该位置处的天然气水合物稳定带底界的预测深度，可以包括：
针对工区的每一位置，通过逐渐增加海底以下地层深度，分别计算上述变地温梯度场所表征的实测温度和上述关联关系，分别在该深度的温度，当两者温度相等时的深度即为甲烷水合物稳定带底界深度。
在上述实施例中，通过对该位置处不同的海底以下地层深度进行逐一确定，当海底以下地层深度满足：在该海底以下地层深度下，上述关联关系中的天然气水合物稳定带底界的预测温度，与变地温梯度场所表征的实测温度相同，确定该海底以下地层深度为该位置处的天然气水合物稳定带底界的预测深度，可确定出工区每一位置处的天然气水合物稳定带底界的预测深度，直接实现天然气水合物稳定带底界的预测，不再需要借助似海底反射进行预测，避免了现有技术中需借助似海底反射，导致在工区无似海底反射的地层中无法进行天然气水合物稳定带底界预测的问题，提高了天然气水合物稳定带底界预测的适用性；同时，通过建立工区的变地温梯度场，可准确描述不同海底以下地层深度下的连续地温梯度数据，避免了现有技术下仅能借助单点地温梯度场或常地温梯度场，导致无法准确实现对海底下天然气水合物稳定带底界预测的问题，提高了天然气水合物稳定带底界预测的预测精度。
具体实施时，本发明实施例提供的一种天然气水合物稳定带底界预测方法，如图7所示，还可以包括：
步骤701：对工区的地震数据进行层位拾取，确定工区的海底深度；
步骤702：根据工区的海底深度，确定工区的天然气水合物稳定带底界的预测温度。
实施例中，可按如下公式对工区的地震数据进行层位拾取，确定工区的海底深度：
其中，Dsf为工区的海底深度，单位为米；Tsf为对工区的地震数据进行层位拾取而得到的海底层位时长，单位为秒；Vw为海水速度，单位为米每秒。
实施例中，对工区的地震数据进行层位拾取，可得到的海底层位时长，可以包括工区的地震数据中海底双程旅行时。
在一个实施例中，根据工区的海底深度，确定工区的天然气水合物稳定带底界的预测温度，可以包括：
基于CSMHyd软件，根据工区的海底深度，确定工区的天然气水合物稳定带底界的预测温度。
在上述实施例中，可利用科罗拉多矿业学院开发的开源程序CSMHyd，结合研究区实际水深，计算研究区甲烷水合物稳定带底界的温度。
在上述实施例中，通过对工区的地震数据进行层位拾取，确定工区的海底深度；进而根据工区的海底深度，可确定出工区的天然气水合物稳定带底界的预测温度，可实现根据工区的天然气水合物稳定带底界的预测温度，确定工区的区域地层深度与天然气水合物稳定带底界的预测温度间的关联关系的目的。
具体实施时，本发明实施例提供的一种天然气水合物稳定带底界预测方法，还可以包括：针对工区的每一位置，根据确定的工区该位置处的天然气水合物稳定带底界的预测深度，计算工区该位置处的天然气水合物稳定带底界的地震旅行时长。
实施例中，按如下公式根据确定的工区该位置处的天然气水合物稳定带底界的预测深度，计算工区该位置处的天然气水合物稳定带底界的地震旅行时长：
其中，Tbghsz为工区该位置处的天然气水合物稳定带底界的地震旅行时长，单位为秒；Dsf为工区的海底深度，单位为米；Dbghsz为工区该位置处的天然气水合物稳定带底界的预测深度，单位为米；Vw为海水速度，单位为米每秒；Vs为海底以下沉积物的平均速度，单位为米每秒。
在上述实施例中，通过计算工区每一位置处的天然气水合物稳定带底界的地震旅行时长，可进一步实现天然气水合物稳定带底界的预测，有助于工作人员根据工区每一位置处的天然气水合物稳定带底界的地震旅行时长，确定天然气水合物稳定带底界。
下面给出一个具体实施例，来说明本发明的方法的具体应用，该实施例中，如图1所示，可以包括如下步骤：
1、首先地震数据层位拾取。
通过对工区的地震数据进行层位拾取，可确定工区的海底深度；
通过工区的海底深度，可确定工区的天然气水合物稳定带底界的预测温度。
2、计算海水深度和似海底反射的深度。根据步骤1拾取的海底和似海底反射的地震双程旅行时间，可通过将时间转换成海底和似海底反射的深度。
3、甲烷水合物稳定带底界温度随深度变化曲线，即工区的区域地层深度与天然气水合物稳定带底界的预测温度间的关联关系。
可利用科罗拉多矿业学院开发的开源程序CSMHyd，结合研究区实际水深，计算研究区甲烷水合物稳定带底界的温度，绘制温度随地层深度变化的点线，通过指数关系式拟合，获取研究区的区域地层深度与甲烷水合物稳定带底界温度的指数关系式，如下所示：
其中，a和b为拟合参数，无量纲；T1为工区的海底以下地层深度为Dmbsf处的天然气水合物稳定带底界的预测温度，单位为摄氏度；Dmbsf为工区的海底以下地层深度，单位为米；Dsf为工区的海底深度，单位为米。
4、实际地层温度随深度变化曲线。根据实际测量的不同深度的地层温度，绘制实际地层温度随深度变化曲线，即针对工区的每一位置，对工区的该位置处不同的海底以下地层深度的实测温度，与不同的海底以下地层深度进行线性拟合，确定对应工区该位置处的不同的海底以下地层深度的实测温度，与海底以下地层深度的线性拟合关系。
5、地温梯度。对步骤4中的曲线进行线性拟合，获取实际地层深度与温度的线性关系式，如下所示：
T2＝Tsf+G×Dmbsf
其中，Tsf为工区该位置处的海底实测温度，单位为摄氏度；G为工区该位置处的地温梯度，单位为摄氏度每米；Dmbsf为工区该位置处的海底以下地层深度，单位为米；T2为工区该位置处的海底以下地层深度为Dmbsf处的实测温度，单位为摄氏度。
6、离散地温梯度场。利用步骤5里的公式对研究区不同位置的实测数据进行线性拟合，获取研究区的离散地温梯度场，即根据工区每一位置处的上述线性拟合关系，得到工区的离散地温梯度场；上述离散地温梯度场用于表征在工区的不同位置处，不同的海底以下地层深度的离散的地温梯度数据。
7、连续变化的变地温梯度场，通过数据插值对步骤6中的离散地温梯度场进行插值，获取研究区的连续变化地温梯度场。
图3是一种变地温梯度场的示意图，图3显示了一种研究区连续变化的地温梯度场，其横坐标为主测线号，其纵坐标为联络测线号，图中G1到G5为实际测量温度的离散位置，通过插值形成如图所示的变地温梯度场。
8、研究区甲烷水合物稳定带底界的温度场预测，根据步骤2的似海底反射深度，通过逐渐增加海底以下地层的深度，利用步骤3里的指数公式预测研究区甲烷水合物稳定带底界的温度场。
9、研究区甲烷水合物稳定带底界的实测温度场，根据步骤7连续变化的地温梯度场，通过逐渐增加海底以下地层的深度，利用步骤5里的线性公式计算研究区甲烷水合物稳定带底界的实测温度场。
10、通过比较步骤8和步骤9里的温度场，在某一海底以下地层深度下，上述关联关系中的天然气水合物稳定带底界的预测温度，与变地温梯度场所表征的实测温度相同的深度，即为甲烷水合物稳定带底界预测深度，即，即针对工区的每一位置，对该位置处不同的海底以下地层深度进行逐一确定，当海底以下地层深度满足：在该海底以下地层深度下，上述关联关系中的天然气水合物稳定带底界的预测温度，与变地温梯度场所表征的实测温度相同，确定该海底以下地层深度为该位置处的天然气水合物稳定带底界的预测深度。
如图1所示，图1中-Tm为甲烷水合物稳定带底界温度，单位℃；Tcal为计算的甲烷水合物稳定带底界温度，单位℃；Tm与Tcal与差的绝对值小于1E-2为10-2。
如图2所示，图2是一种甲烷水合物稳定带底界预测的示意图。图2显示了地层温度随深度的变化。图中线条1周围的点(即图中颜色较深的点)，为实测地层温度随深度变化数据；线条1为图中线条1周围的点，经线性公式拟合得到；图中线条2周围的点(即图中颜色较浅的点)为在研究区内利用CSMHyd软件预测的地层温度随深度变化数据；线条2为对线条2周围的点经指数公式拟合得到。
通过逐渐增加海底以下地层深度，可分别计算线性和指数公式在该深度的温度，当两者温度相等时的深度即为甲烷水合物稳定带底界深度。
11、最后将预测深度通过转换为地震双程旅行时间：
按如下公式计算工区该位置处的天然气水合物稳定带底界的地震旅行时长：
其中，Tbghsz为工区该位置处的天然气水合物稳定带底界的地震旅行时长，单位为秒；Dsf为工区的海底深度，单位为米；Dbghsz为工区该位置处的天然气水合物稳定带底界的预测深度，单位为米；Vw为海水速度，单位为米每秒；Vs为海底以下沉积物的平均速度，单位为米每秒。
如图4所示，图4是海底层位时间拾取图和甲烷水合物稳定带底界的地震双程旅行时间预测图(线条1是拾取的海底层位时间，点线2和点线3分别用常地温梯度场和变地温梯度场预测的甲烷水合物稳定带底界，通过比较可以发现变地温梯度场预测的甲烷水合物稳定带底界与地震剖面中BSR位置较为吻合，其预测精度要明显高于常地温梯度预测深度)。
本发明实施例中，根据工区的天然气水合物稳定带底界的预测温度，确定工区的区域地层深度与天然气水合物稳定带底界的预测温度间的关联关系；上述天然气水合物稳定带底界的预测温度对应不同的区域地层深度；上述区域地层深度为海底深度和海底以下地层深度的组合；根据工区的不同海底以下地层深度的实测温度，建立工区的变地温梯度场；上述变地温梯度场用于表征在工区的不同位置处，不同的海底以下地层深度的连续地温梯度数据；针对工区的每一位置，对该位置处不同的海底以下地层深度进行逐一确定，当海底以下地层深度满足：在该海底以下地层深度下，上述关联关系中的天然气水合物稳定带底界的预测温度，与变地温梯度场所表征的实测温度相同，确定该海底以下地层深度为该位置处的天然气水合物稳定带底界的预测深度，从而可通过工区的天然气水合物稳定带底界的预测温度和工区的不同海底以下地层深度的实测温度，确定出工区每一位置处的天然气水合物稳定带底界的预测深度，直接实现天然气水合物稳定带底界的预测，不再需要借助似海底反射进行预测，避免了现有技术中需借助似海底反射，导致在工区无似海底反射的地层中无法进行天然气水合物稳定带底界预测的问题，提高了天然气水合物稳定带底界预测的适用性；同时，通过建立工区的变地温梯度场，可准确描述不同海底以下地层深度下的连续地温梯度数据，避免了现有技术下仅能借助单点地温梯度场或常地温梯度场，导致无法准确实现对海底下天然气水合物稳定带底界预测的问题，提高了天然气水合物稳定带底界预测的预测精度。
如上述，本发明实施例可针对大区域常地温梯度条件下预测水合物稳定带底界误差较大的问题，通过建立连续的变地温梯度场，构建海底温度、地层温度与甲烷形成水合物温度的指数关系，从而确定天然气水合物稳定带底界，并与地震资料相互验证。本发明能够解决常地温梯度场地震数据预测甲烷水合物稳定带底界误差大的问题，提高了甲烷水合物稳定带的预测精度，也提供了无BSR反射的甲烷水合物稳定带预测新方法。该发明有效提高工作效率，在实际应用中发挥了作用，显示出较好应用潜力。本发明实施例可能够解决地震数据预测甲烷水合物稳定带底界误差大的问题，能够提高甲烷水合物稳定带底界的预测精度，能够用于大面积二维或三维地震数据无BSR处的甲烷水合物稳定带底界预测，有效提高工作效率。
本发明实施例中还提供了一种天然气水合物稳定带底界预测装置，如下面的实施例上述。由于该装置解决问题的原理与天然气水合物稳定带底界预测方法相似，因此该装置的实施可以参见天然气水合物稳定带底界预测方法的实施，重复之处不再赘述。
本发明实施例还提供一种天然气水合物稳定带底界预测装置，用以提高天然气水合物稳定带底界预测的适用性和预测精度，如图8所示，该装置包括：
关联关系确定模块01，用于根据工区的天然气水合物稳定带底界的预测温度，确定工区的区域地层深度与天然气水合物稳定带底界的预测温度间的关联关系；上述天然气水合物稳定带底界的预测温度对应不同的区域地层深度；上述区域地层深度为海底深度与海底以下地层深度之和；
变地温梯度场建立模块02，用于根据工区的不同海底以下地层深度的实测温度，建立工区的变地温梯度场；上述变地温梯度场用于表征在工区的不同位置处，不同的海底以下地层深度的连续地温梯度数据；
预测深度模块03，用于针对工区的每一位置，对该位置处不同的海底以下地层深度进行逐一确定，当海底以下地层深度满足：在该海底以下地层深度下，上述关联关系中的天然气水合物稳定带底界的预测温度，与变地温梯度场所表征的实测温度相同，确定该海底以下地层深度为该位置处的天然气水合物稳定带底界的预测深度。
在一个实施例中，本发明实施例提供的一种天然气水合物稳定带底界预测装置，如图9所示，还可以包括：
预测温度确定模块04，用于：
对工区的地震数据进行层位拾取，确定工区的海底深度；
根据工区的海底深度，确定工区的天然气水合物稳定带底界的预测温度。
在一个实施例中，预测温度确定模块，具体用于：按如下公式对工区的地震数据进行层位拾取，确定工区的海底深度：
其中，Dsf为工区的海底深度，单位为米；Tsf为对工区的地震数据进行层位拾取而得到的海底层位时长，单位为秒；Vw为海水速度，单位为米每秒。
在一个实施例中，关联关系确定模块，具体用于：
对工区的天然气水合物稳定带底界的预测温度，与不同的区域地层深度，进行指数拟合，得到工区的区域地层深度与天然气水合物稳定带底界的预测温度间的关联关系。
在一个实施例中，关联关系确定模块，具体用于：按如下公式进行指数拟合，得到工区的区域地层深度与天然气水合物稳定带底界的预测温度间的关联关系：
其中，a和b为拟合参数，无量纲；T1为工区的海底以下地层深度为Dmbsf处的天然气水合物稳定带底界的预测温度，单位为摄氏度；Dmbsf为工区的海底以下地层深度，单位为米；Dsf为工区的海底深度，单位为米。
在一个实施例中，变地温梯度场建立模块，具体用于：
针对工区的每一位置，对工区的该位置处不同的海底以下地层深度的实测温度，与不同的海底以下地层深度进行线性拟合，确定对应工区该位置处的不同的海底以下地层深度的实测温度，与海底以下地层深度的线性拟合关系；
根据工区每一位置处的上述线性拟合关系，得到工区的离散地温梯度场；上述离散地温梯度场用于表征在工区的不同位置处，不同的海底以下地层深度的离散的地温梯度数据；
对工区的离散地温梯度场进行差值处理，建立工区的变地温梯度场。
在一个实施例中，地温梯度场建立模块，具体用于：按如下公式确定对应工区该位置处的不同的海底以下地层深度的实测温度，与海底以下地层深度的线性拟合关系：
T2＝Tsf+G×Dmbsf
其中，Tsf为工区该位置处的海底实测温度，单位为摄氏度；G为工区该位置处的地温梯度，单位为摄氏度每米；Dmbsf为工区该位置处的海底以下地层深度，单位为米；T2为工区该位置处的海底以下地层深度为Dmbsf处的实测温度，单位为摄氏度。
在一个实施例中，本发明实施例提供的一种天然气水合物稳定带底界预测装置，还可以包括：
地震旅行时长计算模块，用于：
针对工区的每一位置，根据确定的工区该位置处的天然气水合物稳定带底界的预测深度，计算工区该位置处的天然气水合物稳定带底界的地震旅行时长。
在一个实施例中，地震旅行时长计算模块，具体用于：
按如下公式根据确定的工区该位置处的天然气水合物稳定带底界的预测深度，计算工区该位置处的天然气水合物稳定带底界的地震旅行时长：
其中，Tbghsz为工区该位置处的天然气水合物稳定带底界的地震旅行时长，单位为秒；Dsf为工区的海底深度，单位为米；Dbghsz为工区该位置处的天然气水合物稳定带底界的预测深度，单位为米；Vw为海水速度，单位为米每秒；Vs为海底以下沉积物的平均速度，单位为米每秒。
本发明实施例提供一种用于实现上述天然气水合物稳定带底界预测方法中的全部或部分内容的计算机设备的实施例上述计算机设备具体包含有如下内容：
处理器(processor)、存储器(memory)、通信接口(Communications Interface)和总线；其中，上述处理器、存储器、通信接口通过上述总线完成相互间的通信；上述通信接口用于实现相关设备之间的信息传输；该计算机设备可以是台式计算机、平板电脑及移动终端等，本实施例不限于此。在本实施例中，该计算机设备可以参照实施例用于实现天然气水合物稳定带底界预测方法的实施例及用于实现天然气水合物稳定带底界预测装置的实施例进行实施，其内容被合并于此，重复之处不再赘述。
图10为本申请实施例的计算机设备1000的系统构成的示意框图。如图10所示，该计算机设备1000可以包括中央处理器1001和存储器1002；存储器1002耦合到中央处理器1001。值得注意的是，该图10是示例性的；还可以使用其他类型的结构，来补充或代替该结构，以实现电信功能或其他功能。
一实施例中，天然气水合物稳定带底界预测功能可以被集成到中央处理器1001中。其中，中央处理器1001可以被配置为进行如下控制：
根据工区的天然气水合物稳定带底界的预测温度，确定工区的区域地层深度与天然气水合物稳定带底界的预测温度间的关联关系；上述天然气水合物稳定带底界的预测温度对应不同的区域地层深度；上述区域地层深度为海底深度与海底以下地层深度之和；
根据工区的不同海底以下地层深度的实测温度，建立工区的变地温梯度场；上述变地温梯度场用于表征在工区的不同位置处，不同的海底以下地层深度的连续地温梯度数据；
针对工区的每一位置，对该位置处不同的海底以下地层深度进行逐一确定，当海底以下地层深度满足：在该海底以下地层深度下，上述关联关系中的天然气水合物稳定带底界的预测温度，与变地温梯度场所表征的实测温度相同，确定该海底以下地层深度为该位置处的天然气水合物稳定带底界的预测深度。
在另一个实施方式中，天然气水合物稳定带底界预测装置可以与中央处理器1001分开配置，例如可以将天然气水合物稳定带底界预测装置配置为与中央处理器1001连接的芯片，通过中央处理器的控制来实现天然气水合物稳定带底界预测功能。
如图10所示，该计算机设备1000还可以包括：通信模块1003、输入单元1004、音频处理器1005、显示器1006、电源1007。值得注意的是，计算机设备1000也并不是必须要包括图10中所示的所有部件；此外，计算机设备1000还可以包括图10中没有示出的部件，可以参考现有技术。
如图10所示，中央处理器1001有时也称为控制器或操作控件，可以包括微处理器或其他处理器装置和/或逻辑装置，该中央处理器1001接收输入并控制计算机设备1000的各个部件的操作。
其中，存储器1002，例如可以是缓存器、闪存、硬驱、可移动介质、易失性存储器、非易失性存储器或其它合适装置中的一种或更多种。可储存上述与失败有关的信息，此外还可存储执行有关信息的程序。并且中央处理器1001可执行该存储器1002存储的该程序，以实现信息存储或处理等。
输入单元1004向中央处理器1001提供输入。该输入单元1004例如为按键或触摸输入装置。电源1007用于向计算机设备1000提供电力。显示器1006用于进行图像和文字等显示对象的显示。该显示器例如可为LCD显示器，但并不限于此。
该存储器1002可以是固态存储器，例如，只读存储器(ROM)、随机存取存储器(RAM)、SIM卡等。还可以是这样的存储器，其即使在断电时也保存信息，可被选择性地擦除且设有更多数据，该存储器的示例有时被称为EPROM等。存储器1002还可以是某种其它类型的装置。存储器1002包括缓冲存储器1021(有时被称为缓冲器)。存储器1002可以包括应用/功能存储部1022，该应用/功能存储部1022用于存储应用程序和功能程序或用于通过中央处理器1001执行计算机设备1000的操作的流程。
存储器1002还可以包括数据存储部1023，该数据存储部1023用于存储数据，例如联系人、数字数据、图片、声音和/或任何其他由计算机设备使用的数据。存储器1002的驱动程序存储部1024可以包括计算机设备的用于通信功能和/或用于执行计算机设备的其他功能(如消息传送应用、通讯录应用等)的各种驱动程序。
通信模块1003即为经由天线1008发送和接收信号的发送机/接收机1003。通信模块(发送机/接收机)1003耦合到中央处理器1001，以提供输入信号和接收输出信号，这可以和常规移动通信终端的情况相同。
基于不同的通信技术，在同一计算机设备中，可以设置有多个通信模块1003，如蜂窝网络模块、蓝牙模块和/或无线局域网模块等。通信模块(发送机/接收机)1003还经由音频处理器1005耦合到扬声器1009和麦克风1010，以经由扬声器1009提供音频输出，并接收来自麦克风1010的音频输入，从而实现通常的电信功能。音频处理器1005可以包括任何合适的缓冲器、解码器、放大器等。另外，音频处理器1005还耦合到中央处理器1001，从而使得可以通过麦克风1010能够在本机上录音，且使得可以通过扬声器1009来播放本机上存储的声音。
本发明实施例还提供一种计算机可读存储介质，存储有执行上述天然气水合物稳定带底界预测方法的计算机程序。
本发明实施例中，根据工区的天然气水合物稳定带底界的预测温度，确定工区的区域地层深度与天然气水合物稳定带底界的预测温度间的关联关系；上述天然气水合物稳定带底界的预测温度对应不同的区域地层深度；上述区域地层深度为海底深度和海底以下地层深度的组合；根据工区的不同海底以下地层深度的实测温度，建立工区的变地温梯度场；上述变地温梯度场用于表征在工区的不同位置处，不同的海底以下地层深度的连续地温梯度数据；针对工区的每一位置，对该位置处不同的海底以下地层深度进行逐一确定，当海底以下地层深度满足：在该海底以下地层深度下，上述关联关系中的天然气水合物稳定带底界的预测温度，与变地温梯度场所表征的实测温度相同，确定该海底以下地层深度为该位置处的天然气水合物稳定带底界的预测深度，从而可通过工区的天然气水合物稳定带底界的预测温度和工区的不同海底以下地层深度的实测温度，确定出工区每一位置处的天然气水合物稳定带底界的预测深度，直接实现天然气水合物稳定带底界的预测，不再需要借助似海底反射进行预测，避免了现有技术中需借助似海底反射，导致在工区无似海底反射的地层中无法进行天然气水合物稳定带底界预测的问题，提高了天然气水合物稳定带底界预测的适用性；同时，通过建立工区的变地温梯度场，可准确描述不同海底以下地层深度下的连续地温梯度数据，避免了现有技术下仅能借助单点地温梯度场或常地温梯度场，导致无法准确实现对海底下天然气水合物稳定带底界预测的问题，提高了天然气水合物稳定带底界预测的预测精度。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (20)
1.一种天然气水合物稳定带底界预测方法，其特征在于，包括：
根据工区的天然气水合物稳定带底界的预测温度，确定工区的区域地层深度与天然气水合物稳定带底界的预测温度间的关联关系；所述天然气水合物稳定带底界的预测温度对应不同的区域地层深度；所述区域地层深度为海底深度与海底以下地层深度之和；
根据工区的不同海底以下地层深度的实测温度，建立工区的变地温梯度场；所述变地温梯度场用于表征在工区的不同位置处，不同的海底以下地层深度的连续地温梯度数据；
针对工区的每一位置，对该位置处不同的海底以下地层深度进行逐一确定，当海底以下地层深度满足：在该海底以下地层深度下，所述关联关系中的天然气水合物稳定带底界的预测温度，与变地温梯度场所表征的实测温度相同，确定该海底以下地层深度为该位置处的天然气水合物稳定带底界的预测深度。
2.如权利要求1所述的方法，其特征在于，还包括：
对工区的地震数据进行层位拾取，确定工区的海底深度；
根据工区的海底深度，确定工区的天然气水合物稳定带底界的预测温度。
3.如权利要求2所述的方法，其特征在于，按如下公式对工区的地震数据进行层位拾取，确定工区的海底深度：
其中，Dsf为工区的海底深度，单位为米；Tsf为对工区的地震数据进行层位拾取而得到的海底层位时长，单位为秒；Vw为海水速度，单位为米每秒。
4.如权利要求1所述的方法，其特征在于，根据工区的天然气水合物稳定带底界的预测温度，确定工区的区域地层深度与天然气水合物稳定带底界的预测温度间的关联关系，包括：
对工区的天然气水合物稳定带底界的预测温度，与不同的区域地层深度，进行指数拟合，得到工区的区域地层深度与天然气水合物稳定带底界的预测温度间的关联关系。
5.如权利要求4所述的方法，其特征在于，按如下公式进行指数拟合，得到工区的区域地层深度与天然气水合物稳定带底界的预测温度间的关联关系：
其中，a和b为拟合参数，无量纲；T1为工区的海底以下地层深度为Dmbsf处的天然气水合物稳定带底界的预测温度，单位为摄氏度；Dmbsf为工区的海底以下地层深度，单位为米；Dsf为工区的海底深度，单位为米。
6.如权利要求1所述的方法，其特征在于，根据工区的不同海底以下地层深度的实测温度，建立工区的变地温梯度场，包括：
针对工区的每一位置，对工区的该位置处不同的海底以下地层深度的实测温度，与不同的海底以下地层深度进行线性拟合，确定对应工区该位置处的不同的海底以下地层深度的实测温度，与海底以下地层深度的线性拟合关系；
根据工区每一位置处的所述线性拟合关系，得到工区的离散地温梯度场；所述离散地温梯度场用于表征在工区的不同位置处，不同的海底以下地层深度的离散的地温梯度数据；
对工区的离散地温梯度场进行差值处理，建立工区的变地温梯度场。
7.如权利要求6所述的方法，其特征在于，按如下公式确定对应工区该位置处的不同的海底以下地层深度的实测温度，与海底以下地层深度的线性拟合关系：
T2＝Tsf+G×Dmbsf
其中，Tsf为工区该位置处的海底实测温度，单位为摄氏度；G为工区该位置处的地温梯度，单位为摄氏度每米-；Dmbsf为工区该位置处的海底以下地层深度，单位为米；T2为工区该位置处的海底以下地层深度为Dmbsf处的实测温度，单位为摄氏度。
8.如权利要求1所述的方法，其特征在于，还包括：
针对工区的每一位置，根据确定的工区该位置处的天然气水合物稳定带底界的预测深度，计算工区该位置处的天然气水合物稳定带底界的地震旅行时长。
9.如权利要求8所述的方法，其特征在于，按如下公式根据确定的工区该位置处的天然气水合物稳定带底界的预测深度，计算工区该位置处的天然气水合物稳定带底界的地震旅行时长：
其中，Tbghsz为工区该位置处的天然气水合物稳定带底界的地震旅行时长，单位为秒；Dsf为工区的海底深度，单位为米；Dbghsz为工区该位置处的天然气水合物稳定带底界的预测深度，单位为米；Vw为海水速度，单位为米每秒；Vs为海底以下沉积物的平均速度，单位为米每秒。
10.一种天然气水合物稳定带底界预测装置，其特征在于，包括：
关联关系确定模块，用于根据工区的天然气水合物稳定带底界的预测温度，确定工区的区域地层深度与天然气水合物稳定带底界的预测温度间的关联关系；所述天然气水合物稳定带底界的预测温度对应不同的区域地层深度；所述区域地层深度为海底深度与海底以下地层深度之和；
变地温梯度场建立模块，用于根据工区的不同海底以下地层深度的实测温度，建立工区的变地温梯度场；所述变地温梯度场用于表征在工区的不同位置处，不同的海底以下地层深度的连续地温梯度数据；
预测深度模块，用于针对工区的每一位置，对该位置处不同的海底以下地层深度进行逐一确定，当海底以下地层深度满足：在该海底以下地层深度下，所述关联关系中的天然气水合物稳定带底界的预测温度，与变地温梯度场所表征的实测温度相同，确定该海底以下地层深度为该位置处的天然气水合物稳定带底界的预测深度。
11.如权利要求10所述的装置，其特征在于，还包括：
预测温度确定模块，用于：
对工区的地震数据进行层位拾取，确定工区的海底深度；
根据工区的海底深度，确定工区的天然气水合物稳定带底界的预测温度。
12.如权利要求11所述的装置，其特征在于，预测温度确定模块，具体用于：按如下公式对工区的地震数据进行层位拾取，确定工区的海底深度：
其中，Dsf为工区的海底深度，单位为米；Tsf为对工区的地震数据进行层位拾取而得到的海底层位时长，单位为秒；Vw为海水速度，单位为米每秒。
13.如权利要求10所述的装置，其特征在于，关联关系确定模块，具体用于：
对工区的天然气水合物稳定带底界的预测温度，与不同的区域地层深度，进行指数拟合，得到工区的区域地层深度与天然气水合物稳定带底界的预测温度间的关联关系。
14.如权利要求13所述的装置，其特征在于，关联关系确定模块，具体用于：按如下公式进行指数拟合，得到工区的区域地层深度与天然气水合物稳定带底界的预测温度间的关联关系：
其中，a和b为拟合参数，无量纲；T1为工区的海底以下地层深度为Dmbsf处的天然气水合物稳定带底界的预测温度，单位为摄氏度；Dmbsf为工区的海底以下地层深度，单位为米；Dsf为工区的海底深度，单位为米。
15.如权利要求10所述的装置，其特征在于，变地温梯度场建立模块，具体用于：
针对工区的每一位置，对工区的该位置处不同的海底以下地层深度的实测温度，与不同的海底以下地层深度进行线性拟合，确定对应工区该位置处的不同的海底以下地层深度的实测温度，与海底以下地层深度的线性拟合关系；
根据工区每一位置处的所述线性拟合关系，得到工区的离散地温梯度场；所述离散地温梯度场用于表征在工区的不同位置处，不同的海底以下地层深度的离散的地温梯度数据；
对工区的离散地温梯度场进行差值处理，建立工区的变地温梯度场。
16.如权利要求15所述的装置，其特征在于，地温梯度场建立模块，具体用于：按如下公式确定对应工区该位置处的不同的海底以下地层深度的实测温度，与海底以下地层深度的线性拟合关系：
T2＝Tsf+G×Dmbsf
其中，Tsf为工区该位置处的海底实测温度，单位为摄氏度；G为工区该位置处的地温梯度，单位为摄氏度每米；Dmbsf为工区该位置处的海底以下地层深度，单位为米；T2为工区该位置处的海底以下地层深度为Dmbsf处的实测温度，单位为摄氏度。
17.如权利要求10所述的装置，其特征在于，还包括：
地震旅行时长计算模块，用于：
针对工区的每一位置，根据确定的工区该位置处的天然气水合物稳定带底界的预测深度，计算工区该位置处的天然气水合物稳定带底界的地震旅行时长。
18.如权利要求17所述的装置，其特征在于，地震旅行时长计算模块，具体用于：
按如下公式根据确定的工区该位置处的天然气水合物稳定带底界的预测深度，计算工区该位置处的天然气水合物稳定带底界的地震旅行时长：
其中，Tbghsz为工区该位置处的天然气水合物稳定带底界的地震旅行时长，单位为秒；Dsf为工区的海底深度，单位为米；Dbghsz为工区该位置处的天然气水合物稳定带底界的预测深度，单位为米；Vw为海水速度，单位为米每秒；Vs为海底以下沉积物的平均速度，单位为米每秒。
19.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至9任一所述方法。
20.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至9任一所述方法的计算机程序。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110605588.3A CN115480296A (zh) | 2021-05-31 | 2021-05-31 | 天然气水合物稳定带底界预测方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110605588.3A CN115480296A (zh) | 2021-05-31 | 2021-05-31 | 天然气水合物稳定带底界预测方法及装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115480296A true CN115480296A (zh) | 2022-12-16 |
Family
ID=84419611
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110605588.3A Pending CN115480296A (zh) | 2021-05-31 | 2021-05-31 | 天然气水合物稳定带底界预测方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115480296A (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117408186A (zh) * | 2023-12-14 | 2024-01-16 | 青岛海洋地质研究所 | 地震剖面中天然气水合物似海底反射层的判断方法 |
| CN117784239A (zh) * | 2023-12-14 | 2024-03-29 | 青岛海洋地质研究所 | 过去甲烷水合物稳定带底界面位置的判断方法 |
| CN117784239B (zh) * | 2023-12-14 | 2024-06-04 | 青岛海洋地质研究所 | 过去甲烷水合物稳定带底界面位置的判断方法 |
-
2021
- 2021-05-31 CN CN202110605588.3A patent/CN115480296A/zh active Pending
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117408186A (zh) * | 2023-12-14 | 2024-01-16 | 青岛海洋地质研究所 | 地震剖面中天然气水合物似海底反射层的判断方法 |
| CN117408186B (zh) * | 2023-12-14 | 2024-03-12 | 青岛海洋地质研究所 | 地震剖面中天然气水合物似海底反射层的判断方法 |
| CN117784239A (zh) * | 2023-12-14 | 2024-03-29 | 青岛海洋地质研究所 | 过去甲烷水合物稳定带底界面位置的判断方法 |
| CN117784239B (zh) * | 2023-12-14 | 2024-06-04 | 青岛海洋地质研究所 | 过去甲烷水合物稳定带底界面位置的判断方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Bourne et al. | Predictive modelling of naturally fractured reservoirs using geomechanics and flow simulation | |
| EP2718742B1 (en) | Mapping geologic features | |
| RU2590265C2 (ru) | Системы и способы для оценки моментов прорыва флюида в местонахождениях добывающих скважин | |
| CN109188520A (zh) | 薄储层厚度预测方法及装置 | |
| US10060228B2 (en) | Pseudo phase production simulation: a signal processing approach to assess quasi-multiphase flow production via successive analogous step-function relative permeability controlled models in reservoir flow simulation in order to rank multiple petro-physical realizations | |
| CN115480296A (zh) | 天然气水合物稳定带底界预测方法及装置 | |
| US20190310392A1 (en) | Global surface paleo-temperature modeling tool | |
| CN113433587B (zh) | 岩石物理反演方法及系统 | |
| FR3040809A1 (zh) | ||
| CN114384584A (zh) | 裂缝建模方法及装置 | |
| CN112147688A (zh) | 一种储层油气识别方法和装置 | |
| CN114578428A (zh) | 一种速度场高效构建方法及装置 | |
| CN111751879A (zh) | 储层孔隙度预测方法及装置 | |
| AU2013387229A1 (en) | Improvements in determining sub-surface temperature | |
| Craig et al. | Using Maps of Microseismic Events to Define Reservoir Discontinuities | |
| Ballin et al. | New reservoir dynamic connectivity measurement for efficient well placement strategy analysis under depletion | |
| CN112462423A (zh) | 生物礁厚度预测方法及装置 | |
| CN112255679A (zh) | 地震资料绕射深度偏移处理方法及装置 | |
| CN114578445A (zh) | 基于重力资料确定断裂位置的方法及装置 | |
| CN113933901A (zh) | 深度域地震数据拼接方法及装置 | |
| CN113075747B (zh) | 储层裂缝发育区域的预测方法及装置 | |
| Absalyamov et al. | Coupled static and dynamic high resolution modelling on high performance computer clusters | |
| Ainsworth et al. | Predicting Connected Hydrocarbon Volume Uncertainty Ranges: From Static Fault Model to Dynamic Simulator | |
| CN113156496A (zh) | 扇三角洲有利储层的预测方法及装置 | |
| CN114622902A (zh) | 全剖面水平主应力测量方法及装置 |
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 |