CN110703328B - 一种超压界面识别方法、装置、设备及存储介质 - Google Patents
一种超压界面识别方法、装置、设备及存储介质 Download PDFInfo
- Publication number
- CN110703328B CN110703328B CN201910973426.8A CN201910973426A CN110703328B CN 110703328 B CN110703328 B CN 110703328B CN 201910973426 A CN201910973426 A CN 201910973426A CN 110703328 B CN110703328 B CN 110703328B
- Authority
- CN
- China
- Prior art keywords
- overpressure
- seismic
- current
- phase
- data
- 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 47
- 239000003245 coal Substances 0.000 claims description 129
- 230000015572 biosynthetic process Effects 0.000 claims description 36
- 238000005553 drilling Methods 0.000 claims description 28
- 208000035126 Facies Diseases 0.000 claims description 26
- 239000011435 rock Substances 0.000 claims description 10
- 238000004590 computer program Methods 0.000 claims description 8
- 238000011161 development Methods 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 7
- 239000011148 porous material Substances 0.000 description 7
- 238000012545 processing Methods 0.000 description 7
- 230000003287 optical effect Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 239000003086 colorant Substances 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000005056 compaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000008707 rearrangement Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000007704 transition Effects 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
-
- 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
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
技术领域
本发明实施例涉及地层孔隙压力预测技术领域，尤其涉及一种超压界面识别方法、装置、设备及存储介质。
背景技术
世界上的许多盆地均存在超压现象。常压是指地层的孔隙压力梯度范围大于等于1.0g/cm3而小于1.2g/cm3，超压是指地层的孔隙压力梯度大于1.2g/cm3。通常地层孔隙压力梯度等于1.2g/cm3的深度位置称之为超压界面。
近年来，常规超压界面预测方法是根据超压地层随着地层加深速度降低的特征，采用基于欠压实成因的Eaton法以及基于流体膨胀成因的Bowers法，这些方法均是基于超压地层具有地层速度降低的特征，即在钻前获得地层速度，通过地层速度定量预测地层孔隙压力。
但是，随着中深层地震资料品质降低，地震速度较难求准，因此超压界面的定量预测存在不确定性。
发明内容
本发明实施例的目的是提供一种超压界面识别方法、装置、设备及存储介质，以提升超压界面预测准确性。
第一方面，本发明实施例提供了一种超压界面识别方法，其中，包括：
确定当前区域的地震剖面图中的煤系地层数据，并确定所述煤系地层数据的煤系地层顶面、煤系地层底面以及位于所述煤系地层顶面和所述煤系地层底面之间的至少一个同相轴；
获取各所述同相轴的振幅数据，根据所述振幅数据确定所述当前区域地震剖面图的超压带地震相；
根据所述超压带地震相确定所述当前区域地震剖面图的超压界面。
第二方面，本发明实施例还提供了一种超压界面识别装置，其中，包括：
第一确定模块，用于确定当前区域的地震剖面图中的煤系地层数据，并确定所述煤系地层数据的煤系地层顶面、煤系地层底面以及位于所述煤系地层顶面和所述煤系地层底面之间的至少一个同相轴；
获取模块，用于获取各所述同相轴的振幅数据，根据所述振幅数据确定所述当前区域地震剖面图的超压带地震相；
第二确定模块，用于根据所述超压带地震相确定所述当前区域地震剖面图的超压界面。
第三方面，本发明实施例还提供了一种超压界面识别设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其中，所述处理器执行所述计算机程序时实现如第一方面中任一项所述的超压界面识别方法。
第四方面，本发明实施例还提供了一种包含计算机可执行指令的存储介质，其中，所述计算机可执行指令在由计算机处理器执行时实现如第一方面中任一项所述的超压界面识别方法。
本发明实施例的技术方案，通过确定当前区域的地震剖面图中的煤系地层数据，并确定煤系地层数据的煤系地层顶面、煤系地层底面以及位于煤系地层顶面和所述煤系地层底面之间的至少一个同相轴，获取各同相轴的振幅数据，根据振幅数据确定当前区域地震剖面图的超压带地震相，根据超压带地震相确定当前区域地震剖面图的超压界面，解决了现有技术中随着深层地震资料品质降低，地震速度较难求准的问题，达到通过确定煤系地层区域并对煤系地层区域的地震数据的反射特征确定超压界面的目的，实现提升超压界面预测准确性的效果。
附图说明
图1为本发明实施例一提供的一种超压界面识别方法的流程示意图；
图2为本发明实施例二提供的一种超压界面识别方法的流程示意图；
图3为本发明实施例二提供的一种地震数据的界面图；
图4为本发明实施例三提供的一种超压界面识别方法的流程示意图；
图5为本发明实施例三提供的一种超压界面识别方法的界面图；
图6为本发明实施例四提供的一种超压界面识别装置的结构示意图；
图7为本发明实施例五提供的一种超压界面识别设备的结构示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释本发明，而非对本发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与本发明相关的部分而非全部结构。
实施例一
图1为本发明实施例一提供的一种超压界面识别方法的流程示意图，本实施例可适用于通过地震数据识别超级界面的情况，该方法可以由超压界面识别装置来执行，其中该装置可由软件和/或硬件实现，并一般集成在终端中。具体参见图1所示，该方法可以包括如下步骤：
S110，确定当前区域的地震剖面图中的煤系地层数据，并确定煤系地层数据的煤系地层顶面、煤系地层底面以及位于煤系地层顶面和煤系地层底面之间的至少一个同相轴。
其中，当前区域可以为海洋、陆地以及其他地区。一般情况下，对于煤系地层区域，在煤层的生排烃作用下，会促进超压界面形成，因而，从地震相数据中确定了煤系地层区域就可以进一步确定超压界面。
可选地，当前区域的地震剖面图的煤系地层数据可以包括以下至少一个特征：分布方位为沉积盆地的边缘区域、地层类型为碎屑岩地层、煤层发育情况为发育煤层以及岩性类别为砂岩、泥岩和煤岩中的至少一种。这样，用户可以通过判断当前区域的是否位于大型沉积盆地的边缘、地层是否为碎屑岩地层、岩性是否以砂岩、泥岩为主以及地层中是否发育煤层，确定当前区域的地震数据是否为煤系地层数据。
进一步地，确定了煤系地层数据后，将煤系地层数据输入到超压界面识别设备，超压界面识别设备将不同振幅强度的数据显示为不同颜色，每一种颜色对应一个振幅强度，超压界面识别设备确定出煤系地层顶面、煤系地层底面，并通过分析颜色以及振幅强度确定至少一个同相轴。其中，同相轴为相位相同的振幅数据的极值的连线，用来记录地震数据的规律性波动。
S120，获取各同相轴的振幅数据，根据振幅数据确定当前区域地震剖面图的超压带地震相。
在煤系地层区域处，当煤层频繁出现后，煤层的生排作用可以促进超压的形成，这样该区域的地层具有煤层发育与超压出现的双重特征，由于煤系地层区域的波阻抗强度比砂岩区域以及泥岩区域小，因此，在煤系地层区域的地震数据的显示界面上，会出现振幅较强且具有负反射特征的地震相数据，这样，超压界面识别设备根据地震数据的振幅强度以及反射特征确定超压带地震相。其中，波阻抗为地震应力在岩石区域、煤层出现区域以及砂岩区域中传播时，运动着的岩石区域、煤层出现区域以及砂岩区域的质点数据产生单位速度所需的扰动力，它可以表明岩石区域、煤层出现区域以及砂岩区域对扰动力的抵抗能力。
S130，根据超压带地震相确定当前区域地震剖面图的超压界面。
可以理解的是，超压界面一般出现在煤系地层区域的超压带地震相范围内，因而，通过识别超压带地震相可以进一步确定超压界面。大量研究得出，超压界面通常位于超压带地震相的顶部，也就是说，超压界面一般对应超压带地震相的第一根负相位的同相轴，因此，超压界面识别设备确定了超压带地震相之后，根据各同相轴的位置就可以确定超压界面，达到定量预测超压界面的目的。
本发明实施例提供了一种超压界面识别方法，通过确定当前区域的地震剖面图中的煤系地层数据，并确定煤系地层数据的煤系地层顶面、煤系地层底面以及位于煤系地层顶面和所述煤系地层底面之间的至少一个同相轴，获取各同相轴的振幅数据，根据振幅数据确定当前区域地震剖面图的超压带地震相，根据超压带地震相确定当前区域地震剖面图的超压界面，解决了现有技术中随着深层地震资料品质降低，地震速度较难求准的问题，达到通过确定煤系地层区域并对煤系地层区域的地震数据的反射特征确定超压界面的目的，实现提升超压界面预测准确性的效果。
实施例二
图2为本发明实施例二提供的一种超压界面识别方法的流程示意图。本实施例的技术方案在上述实施例的基础上进行了细化，可选地，所述根据所述振幅数据确定所述当前区域地震剖面图的超压带地震相，包括：对于当前同相轴，将所述当前同相轴的振幅数据与第一预设振幅阈值进行比较，以及将朝向所述煤系地层顶面方向的所述当前同相轴相邻的三个同相轴的振幅数据分别与第二预设振幅阈值进行比较；如果所述当前同相轴的振幅数据小于所述第一预设振幅阈值，且所述当前同相轴相邻的三个同相轴的振幅数据的绝对值均小于所述第二预设振幅阈值，则将所述当前同相轴至与所述煤系地层底面相邻的同相轴之间的区域确定为所述当前区域地震剖面图的超压带地震相。具体参见图2所示，本实施例的方法可以包括如下步骤：
S210，确定当前区域的地震剖面图中的煤系地层数据，并确定煤系地层数据的煤系地层顶面、煤系地层底面以及位于煤系地层顶面和煤系地层底面之间的至少一个同相轴。
S220，在获取各同相轴的振幅数据后，对于当前同相轴，将当前同相轴的振幅数据与第一预设振幅阈值进行比较，以及将朝向煤系地层顶面方向的当前同相轴相邻的三个同相轴的振幅数据分别与第二预设振幅阈值进行比较。
其中，当前同相轴可以为煤系地层顶面、煤系地层底面之间的任意一根同相轴。在上述实施例的基础上，在超压带地震相处，由于煤层数据波阻抗较小，超压带地震相会出现纵向连续的强振幅波组特征，在常压带地震相处，由于煤层数据波阻抗较大，常压带地震相会出现中弱振幅波组的地震相特征，这样，可以通过比较地震数据的振幅强度与预设的振幅强度的大小确定出超压带地震相。
一般情况下，在常压带地震相处，纵向连续的若干个同相轴的振幅数据小于第二预设振幅阈值，在超压带地震相的临界区域，同相轴的振幅数据发生变化，这样，通过将当前同相轴的振幅数据与第一预设振幅阈值进行比较，以及将朝向煤系地层顶面方向的当前同相轴相邻的三个同相轴的振幅数据分别与第二预设振幅阈值进行比较，就可以确定出超压地震相以及常压地震相。
S230，如果当前同相轴的振幅数据小于第一预设振幅阈值，且当前同相轴相邻的三个同相轴的振幅数据的绝对值均小于第二预设振幅阈值，则将当前同相轴至与煤系地层底面相邻的同相轴之间的区域确定为当前区域地震剖面图的超压带地震相。
可选地，第一预设振幅阈值为常压带地震相与超压带地震相过渡区域的振幅强度，第二预设振幅阈值为超压带地震相的最小振幅强度，并且超压带地震相的任一同相轴的振幅数据的绝对值超过第二预设振幅阈值。其中，第一预设振幅阈值，一般取值为-0.67M，第二预设振幅阈值一般取值为0.67M。
示例性地，参见图3，在地震数据的显示界面上，超压界面识别设备确定了煤系地层顶面1和煤系地层底面2以及在获取煤系地层顶面1和煤系地层地面2范围内的各同相轴的振幅数据，如果同相轴5处的振幅数据小于-0.67M，并且当前同相轴5相邻的三个同相轴的振幅数据的绝对值均小于0.67M，则将当前同相轴5至与煤系地层底面2相邻的同相轴之间的区域确定为当前区域地震剖面图的超压带地震相，即图3中的4区域为超压带地震相。
可选地，如图3所示，还可以将与煤系地层顶面1相邻的同相轴与当前同相轴5之间的区域确定为当前区域地震剖面图的常压带地震相，即图3中的3区域为常压带地震相。其中，常压带地震相的岩性类别主要包括砂岩和泥岩，在特定深度处，常压带地震相的砂岩的波阻抗和泥岩的波阻抗的差值小于误差阈值。
S240，根据超压带地震相确定当前区域地震剖面图的超压界面。
在上述实施例的基础上，超压界面一般对应超压带地震相的第一根负相位的同相轴，由图3可以看出，当前同相轴5为第一根负相位的同相轴，可以将当前同相轴5对应的界面确定为超压界面。
本发明实施例提供了一种超压界面识别方法，超压界面识别设备通过将当前同相轴的振幅数据与第一预设振幅阈值进行比较，以及将朝向煤系地层顶面方向的当前同相轴相邻的三个同相轴的振幅数据分别与第二预设振幅阈值进行比较，根据比较结果确定超压带地震相以及超压界面，解决了现有技术中随着深层地震资料品质降低，地震速度较难求准的问题，达到通过确定煤系地层区域并对煤系地层区域的地震数据的反射特征确定超压界面的目的，实现提升超压界面预测准确性的效果。
实施例三
图4为本发明实施例三提供的一种超压界面识别方法的流程示意图。本实施例的技术方案在上述实施例的基础上增加了新的步骤，可选地，还包括：获取所述当前区域的多组钻井资料，根据每组所述钻井资料分别确定所述当前区域的压力系数曲线；其中，所述压力系数曲线包括所述当前区域的钻井深度和钻井压力值；将处于相同深度区间的多个所述压力系数曲线的多个压力值分别与所述超压界面进行匹配；获取匹配数量与设定阈值的比较结果，根据所述比较结果调整所述深度区间的深度范围，以更新处于调整后的深度区间所述超压界面的面积。具体参见图4所示，本实施例的方法可以包括如下步骤：
S410，确定当前区域的地震剖面图中的煤系地层数据，并确定煤系地层数据的煤系地层顶面、煤系地层底面以及位于煤系地层顶面和煤系地层底面之间的至少一个同相轴。
S420，获取各同相轴的振幅数据，根据振幅数据确定当前区域地震剖面图的超压带地震相。
S430，根据超压带地震相确定当前区域地震剖面图的超压界面。
S440，获取当前区域的多组钻井资料，根据每组钻井资料分别确定当前区域的压力系数曲线。
S450，将处于相同深度区间的多个压力系数曲线的多个压力值分别与超压界面进行匹配。
S460，获取匹配数量与设定阈值的比较结果，根据比较结果调整深度区间的深度范围，以更新处于调整后的深度区间超压界面的面积。
可以理解的是，通过确定煤系地层区域并根据煤系地层区域的地震数据的反射特征确定超压界面可能会出现误差。为了提高超压界面确定的准确性，用户可以对当前区域进行多次钻井，并将得到的多组钻井资料输入到超压界面识别设备中，超压界面识别设备分别根据每组钻井资料确定当前区域的每组钻井的压力系数曲线，其中，该压力系数曲线用于描述钻井深度与地层孔隙压力梯度的对应关系。一般情况下，当地层孔隙的压力值超过1.2g/cm3时，将超过1.2g/cm3对应深度的平面作为超压界面。
如图5所示，伽马曲线为每组钻井曲线，压力系数曲线对应每组钻井资料，超压面为通过分析煤系地层区域的地震数据的反射特征确定的超压界面，如果超压面与每组钻井资料发生重叠，则说明该组钻井资料与超压面匹配。可理解的是，超压界面识别设备确定了每组钻井的压力系数曲线以及地震剖面图的超压界面之后，将处于相同深度区间的每组压力系数曲线的多个压力值分别与超压界面进行匹配，如果匹配数量超过了超过设定阈值，则说明通过分析煤系地层区域的地震数据的反射特征确定超压界面较为准确，如果匹配数量小于设定阈值，则可以在误差范围内调整深度区间的深度范围，更新处于调整后的深度区间超压界面的面积。
示例性地，如果超压界面识别设备通过分析煤系地层区域的地震数据的反射特征确定超压界面位于[2800，3500]米处，并得到位于[2800，3500]米处的10组钻井资料，分别将该10组钻井资料的压力值与超压界面进行匹配，如果匹配数量超过6组，则说明通过分析煤系地层区域的地震数据的反射特征确定超压界面准确，如果匹配数量小于6组，在深度误差范围内调整深度区间，例如将深度区间调整为[3000，3600]米，并重新获得调整后的深度区间的超压界面，直到调整后的深度区间的超压界面与通过钻井资料确定的超压界面的匹配数量不小于6组，并确定该深度区间的超压界面的面积。
本发明实施例提供了一种超压界面识别方法，通过将多组钻井资料确定的超压界面与根据煤系地层区域的地震数据的反射特征确定超压界面进行匹配，进一步验证确定的超压界面是否合理，以及在深度误差范围内调整深度区间，重新确定超压界面，提高了超压界面确定的准确性。
实施例四
图6为本发明实施例四提供的一种超压界面识别装置的结构示意图。参见图6所示，该装置包括：第一确定模块61、获取模块62以及第二确定模块63。
其中，第一确定模块61，用于确定当前区域的地震剖面图中的煤系地层数据，并确定煤系地层数据的煤系地层顶面、煤系地层底面以及位于煤系地层顶面和煤系地层底面之间的至少一个同相轴；获取模块62，用于获取各同相轴的振幅数据，根据振幅数据确定当前区域地震剖面图的超压带地震相；第二确定模块63，用于根据超压带地震相确定当前区域地震剖面图的超压界面。
在上述各技术方案的基础上，当前区域的地震剖面图的煤系地层数据包括以下至少一个特征：分布方位为沉积盆地的边缘区域、地层类型为碎屑岩地层、煤层发育情况为发育煤层以及岩性类别为砂岩、泥岩和煤岩中的至少一种。
在上述各技术方案的基础上，获取模块62还用于，对于当前同相轴，将当前同相轴的振幅数据与第一预设振幅阈值进行比较，以及将朝向煤系地层顶面方向的当前同相轴相邻的三个同相轴的振幅数据分别与第二预设振幅阈值进行比较；如果当前同相轴的振幅数据小于第一预设振幅阈值，且当前同相轴相邻的三个同相轴的振幅数据的绝对值均小于第二预设振幅阈值，则将当前同相轴至与煤系地层底面相邻的同相轴之间的区域确定为当前区域地震剖面图的超压带地震相。
在上述各技术方案的基础上，第二确定模块63还用于，将当前同相轴对应的界面确定为超压界面。
在上述各技术方案的基础上，还包括：第三确定模块。第三确定模块，用于将与煤系地层顶面相邻的同相轴与当前同相轴之间的区域确定为当前区域地震剖面图的常压带地震相，其中，常压带地震相的岩性类别包括砂岩和泥岩，在特定深度处，常压带地震相的砂岩的波阻抗和泥岩的波阻抗的差值小于误差阈值。
在上述各技术方案的基础上，超压带地震相的任一同相轴的振幅数据的绝对值超过第二预设振幅阈值。
在上述各技术方案的基础上，还包括：调整模块。调整模块，用于获取当前区域的多组钻井资料，根据每组钻井资料分别确定当前区域的压力系数曲线；其中，压力系数曲线包括当前区域的钻井深度和钻井压力值；将处于相同深度区间的多个压力系数曲线的多个压力值分别与超压界面进行匹配；获取匹配数量与设定阈值的比较结果，根据比较结果调整深度区间的深度范围，以更新处于调整后的深度区间超压界面的面积。
本实施例提供了一种超压界面识别装置，通过确定当前区域的地震剖面图中的煤系地层数据，并确定煤系地层数据的煤系地层顶面、煤系地层底面以及位于煤系地层顶面和煤系地层底面之间的至少一个同相轴，获取各同相轴的振幅数据，根据振幅数据确定当前区域地震剖面图的超压带地震相，根据超压带地震相确定当前区域地震剖面图的超压界面，解决了现有技术中随着深层地震资料品质降低，地震速度较难求准的问题，达到通过确定煤系地层区域并对煤系地层区域的地震数据的反射特征确定超压界面的目的，实现提升超压界面预测准确性的效果。
实施例五
图7为本发明实施例五提供的一种超压界面识别设备的结构示意图。图7示出了适于用来实现本发明实施方式的示例性超压界面识别设备12的框图。图7显示的超压界面识别设备12仅仅是一个示例，不应对本发明实施例的功能和使用范围带来任何限制。
如图7所示，超压界面识别设备12以通用计算设备的形式表现。超压界面识别设备12的组件可以包括但不限于：一个或者多个处理器或者处理单元16，系统存储器28，连接不同系统组件(包括系统存储器28和处理单元16)的总线18。
总线18表示几类总线结构中的一种或多种，包括存储器总线或者存储器控制器，外围总线，图形加速端口，处理器或者使用多种总线结构中的任意总线结构的局域总线。举例来说，这些体系结构包括但不限于工业标准体系结构(ISA)总线，微通道体系结构(MAC)总线，增强型ISA总线、视频电子标准协会(VESA)局域总线以及外围组件互连(PCI)总线。
超压界面识别设备12典型地包括多种计算机系统可读介质。这些介质可以是任何能够被超压界面识别设备12访问的可用介质，包括易失性和非易失性介质，可移动的和不可移动的介质。
系统存储器28可以包括易失性存储器形式的计算机系统可读介质，例如随机存取存储器(RAM)30和/或高速缓存存储器32。超压界面识别设备12可以进一步包括其它可移动/不可移动的、易失性/非易失性计算机系统存储介质。仅作为举例，存储系统34可以用于读写不可移动的、非易失性磁介质(图7未显示，通常称为“硬盘驱动器”)。尽管图7中未示出，可以提供用于对可移动非易失性磁盘(例如“软盘”)读写的磁盘驱动器，以及对可移动非易失性光盘(例如CD-ROM,DVD-ROM或者其它光介质)读写的光盘驱动器。在这些情况下，每个驱动器可以通过一个或者多个数据介质接口与总线18相连。存储器28可以包括至少一个程序产品，该程序产品具有一组(例如超压界面识别装置的第一确定模块61、获取模块62以及第二确定模块63)程序模块，这些程序模块被配置以执行本发明各实施例的功能。
具有一组(例如第一确定模块61、获取模块62以及第二确定模块63)程序模块46的程序/实用工具44，可以存储在例如存储器28中，这样的程序模块46包括但不限于操作系统、一个或者多个应用程序、其它程序模块以及程序数据，这些示例中的每一个或某种组合中可能包括网络环境的实现。程序模块46通常执行本发明所描述的实施例中的功能和/或方法。
超压界面识别设备12也可以与一个或多个外部设备14(例如键盘、指向设备、显示器24等)通信，还可与一个或者多个使得用户能与该超压界面识别设备12交互的设备通信，和/或与使得该超压界面识别设备12能与一个或多个其它计算设备进行通信的任何设备(例如网卡，调制解调器等等)通信。这种通信可以通过输入/输出(I/O)接口22进行。并且，超压界面识别设备12还可以通过网络适配器20与一个或者多个网络(例如局域网(LAN)，广域网(WAN)和/或公共网络，例如因特网)通信。如图所示，网络适配器20通过总线18与超压界面识别设备12的其它模块通信。应当明白，尽管图中未示出，可以结合超压界面识别设备12使用其它硬件和/或软件模块，包括但不限于：微代码、设备驱动器、冗余处理单元、外部磁盘驱动阵列、RAID系统、磁带驱动器以及数据备份存储系统等。
处理单元16通过运行存储在系统存储器28中的程序，从而执行各种功能应用以及数据处理，例如实现本发明实施例所提供的一种超压界面识别方法，该方法包括：
确定当前区域的地震剖面图中的煤系地层数据，并确定煤系地层数据的煤系地层顶面、煤系地层底面以及位于煤系地层顶面和煤系地层底面之间的至少一个同相轴；
获取各同相轴的振幅数据，根据振幅数据确定当前区域地震剖面图的超压带地震相；
根据超压带地震相确定当前区域地震剖面图的超压界面。处理单元16通过运行存储在系统存储器28中的程序，从而执行各种功能应用以及数据处理，例如实现本发明实施例所提供的一种超压界面识别方法。
当然，本领域技术人员可以理解，处理器还可以实现本发明任意实施例所提供的一种超压界面识别方法的技术方案。
实施例六
本发明实施例六还提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现如本发明实施例所提供的一种超压界面识别方法，该方法包括：
确定当前区域的地震剖面图中的煤系地层数据，并确定煤系地层数据的煤系地层顶面、煤系地层底面以及位于煤系地层顶面和煤系地层底面之间的至少一个同相轴；
获取各同相轴的振幅数据，根据振幅数据确定当前区域地震剖面图的超压带地震相；
根据超压带地震相确定当前区域地震剖面图的超压界面。当然，本发明实施例所提供的一种计算机可读存储介质，其上存储的计算机程序不限于如上的方法操作，还可以执行本发明任意实施例所提供的一种超压界面识别方法中的相关操作。
本发明实施例的计算机存储介质，可以采用一个或多个计算机可读的介质的任意组合。计算机可读介质可以是计算机可读信号介质或者计算机可读存储介质。计算机可读存储介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的系统、系统或器件，或者任意以上的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括：具有一个或多个导线的电连接、便携式计算机磁盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑磁盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本文件中，计算机可读存储介质可以是任何包含或存储程序的有形介质，该程序可以被指令执行系统、系统或者器件使用或者与其结合使用。
计算机可读的信号介质可以包括在煤系地层数据、煤系地层顶面、煤系地层底面以及位于所述煤系地层顶面和所述煤系地层底面之间的至少一个同相轴等，其中承载了计算机可读的程序代码。这种传播的煤系地层数据、煤系地层顶面、煤系地层底面以及位于所述煤系地层顶面和所述煤系地层底面之间的至少一个同相轴等形式。计算机可读的信号介质还可以是计算机可读存储介质以外的任何计算机可读介质，该计算机可读介质可以发送、传播或者传输用于由指令执行系统、系统或者器件使用或者与其结合使用的程序。
计算机可读介质上包含的程序代码可以用任何适当的介质传输，包括——但不限于无线、电线、光缆、RF等等，或者上述的任意合适的组合。
可以以一种或多种程序设计语言或其组合来编写用于执行本发明操作的计算机程序代码，程序设计语言包括面向对象的程序设计语言—诸如Java、Smalltalk、C++，还包括常规的过程式程序设计语言—诸如”C”语言或类似的程序设计语言。程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络——包括局域网(LAN)或广域网(WAN)—连接到用户计算机，或者，可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。
值得注意的是，上述超压界面识别装置的实施例中，所包括的各个模块只是按照功能逻辑进行划分的，但并不局限于上述的划分，只要能够实现相应的功能即可；另外，各功能单元的具体名称也只是为了便于相互区分，并不用于限制本发明的保护范围。
注意，上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解，本发明不限于这里所述的特定实施例，对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此，虽然通过以上实施例对本发明进行了较为详细的说明，但是本发明不仅仅限于以上实施例，在不脱离本发明构思的情况下，还可以包括更多其他等效实施例，而本发明的范围由所附的权利要求范围决定。
Claims (9)
1.一种超压界面识别方法，其特征在于，包括：
确定当前区域的地震剖面图中的煤系地层数据，并确定所述煤系地层数据的煤系地层顶面、煤系地层底面以及位于所述煤系地层顶面和所述煤系地层底面之间的至少一个同相轴；
获取各所述同相轴的振幅数据，根据所述振幅数据确定所述当前区域地震剖面图的超压带地震相；
根据所述超压带地震相确定所述当前区域地震剖面图的超压界面；
所述根据所述振幅数据确定所述当前区域地震剖面图的超压带地震相，包括：
对于当前同相轴，将所述当前同相轴的振幅数据与第一预设振幅阈值进行比较，以及将朝向所述煤系地层顶面方向的所述当前同相轴相邻的三个同相轴的振幅数据分别与第二预设振幅阈值进行比较；
如果所述当前同相轴的振幅数据小于所述第一预设振幅阈值，且所述当前同相轴相邻的三个同相轴的振幅数据的绝对值均小于所述第二预设振幅阈值，则将所述当前同相轴至与所述煤系地层底面相邻的同相轴之间的区域确定为所述当前区域地震剖面图的超压带地震相。
2.根据权利要求1所述的方法，其特征在于，所述当前区域的地震剖面图的煤系地层数据包括以下至少一个特征：分布方位为沉积盆地的边缘区域、地层类型为碎屑岩地层、煤层发育情况为发育煤层以及岩性类别为砂岩、泥岩和煤岩中的至少一种。
3.根据权利要求1所述的方法，其特征在于，所述根据所述超压带地震相确定所述当前区域地震剖面图的超压界面，包括：
将所述当前同相轴对应的界面确定为所述超压界面。
4.根据权利要求1所述的方法，其特征在于，所述方法还包括：
将与所述煤系地层顶面相邻的同相轴与所述当前同相轴之间的区域确定为所述当前区域地震剖面图的常压带地震相，其中，所述常压带地震相的岩性类别包括砂岩和泥岩，在特定深度处，所述常压带地震相的砂岩的波阻抗和泥岩的波阻抗的差值小于误差阈值。
5.根据权利要求1所述的方法，其特征在于，所述超压带地震相的任一同相轴的振幅数据的绝对值超过所述第二预设振幅阈值。
6.根据权利要求1所述的方法，其特征在于，还包括：
获取所述当前区域的多组钻井资料，根据每组所述钻井资料分别确定所述当前区域的压力系数曲线；其中，所述压力系数曲线包括所述当前区域的钻井深度和钻井压力值；
将处于相同深度区间的多个所述压力系数曲线的多个压力值分别与所述超压界面进行匹配；
获取匹配数量与设定阈值的比较结果，根据所述比较结果调整所述深度区间的深度范围，以更新处于调整后的深度区间所述超压界面的面积。
7.一种超压界面识别装置，其特征在于，包括：
第一确定模块，用于确定当前区域的地震剖面图中的煤系地层数据，并确定所述煤系地层数据的煤系地层顶面、煤系地层底面以及位于所述煤系地层顶面和所述煤系地层底面之间的至少一个同相轴；
获取模块，用于获取各所述同相轴的振幅数据，根据所述振幅数据确定所述当前区域地震剖面图的超压带地震相；
第二确定模块，用于根据所述超压带地震相确定所述当前区域地震剖面图的超压界面；
获取模块，还用于对于当前同相轴，将当前同相轴的振幅数据与第一预设振幅阈值进行比较，以及将朝向煤系地层顶面方向的当前同相轴相邻的三个同相轴的振幅数据分别与第二预设振幅阈值进行比较；如果当前同相轴的振幅数据小于第一预设振幅阈值，且当前同相轴相邻的三个同相轴的振幅数据的绝对值均小于第二预设振幅阈值，则将当前同相轴至与煤系地层底面相邻的同相轴之间的区域确定为当前区域地震剖面图的超压带地震相。
8.一种超压界面识别设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现如权利要求1-6中任一项所述的超压界面识别方法。
9.一种包含计算机可执行指令的存储介质，其特征在于，所述计算机可执行指令在由计算机处理器执行时实现如权利要求1-6中任一项所述的超压界面识别方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910973426.8A CN110703328B (zh) | 2019-10-14 | 2019-10-14 | 一种超压界面识别方法、装置、设备及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910973426.8A CN110703328B (zh) | 2019-10-14 | 2019-10-14 | 一种超压界面识别方法、装置、设备及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110703328A CN110703328A (zh) | 2020-01-17 |
| CN110703328B true CN110703328B (zh) | 2021-09-24 |
Family
ID=69199318
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910973426.8A Active CN110703328B (zh) | 2019-10-14 | 2019-10-14 | 一种超压界面识别方法、装置、设备及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110703328B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115680638B (zh) * | 2021-07-26 | 2024-07-23 | 中国石油化工股份有限公司 | 一种利用压力衰减梯度识别超压顶封层的方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6374186B1 (en) * | 2000-03-24 | 2002-04-16 | Petrophysical Consulting, Inc. | Method for overpressure detection from compressional-and- shear-wave data |
| CN105093274A (zh) * | 2014-05-07 | 2015-11-25 | 中国石油化工股份有限公司 | 一种水力压裂裂缝震源机制的反演方法及系统 |
| CN106481337A (zh) * | 2015-08-27 | 2017-03-08 | 中国石油化工股份有限公司 | 超压顶界面的预测方法 |
| CN109061729A (zh) * | 2018-08-22 | 2018-12-21 | 西安石油大学 | 一种高温高压气藏含气敏感曲线重构方法 |
| CN109143353A (zh) * | 2018-08-20 | 2019-01-04 | 电子科技大学 | 一种基于深度卷积生成对抗网络的叠前地震波形分类方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7349807B2 (en) * | 2004-03-08 | 2008-03-25 | Geomechanics International, Inc. | Quantitative risk assessment applied to pore pressure prediction |
-
2019
- 2019-10-14 CN CN201910973426.8A patent/CN110703328B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6374186B1 (en) * | 2000-03-24 | 2002-04-16 | Petrophysical Consulting, Inc. | Method for overpressure detection from compressional-and- shear-wave data |
| CN105093274A (zh) * | 2014-05-07 | 2015-11-25 | 中国石油化工股份有限公司 | 一种水力压裂裂缝震源机制的反演方法及系统 |
| CN106481337A (zh) * | 2015-08-27 | 2017-03-08 | 中国石油化工股份有限公司 | 超压顶界面的预测方法 |
| CN109143353A (zh) * | 2018-08-20 | 2019-01-04 | 电子科技大学 | 一种基于深度卷积生成对抗网络的叠前地震波形分类方法 |
| CN109061729A (zh) * | 2018-08-22 | 2018-12-21 | 西安石油大学 | 一种高温高压气藏含气敏感曲线重构方法 |
Non-Patent Citations (1)
| Title |
|---|
| "莺歌海盆地深层高温超压地层弱振幅下储层预测技术研究与应用";李芳 等;《中国石油学会2019 年物探技术研讨会》;20190930;第1064-1066页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110703328A (zh) | 2020-01-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7983845B2 (en) | Method and system for analyzing a laminated sand/shale formation | |
| CN106324676B (zh) | 一种确定断层封闭性的方法及装置 | |
| US12117580B2 (en) | Evaluation of rock physical properties from drill sounds through minimizing the effect of the drill bit rotation | |
| CN106896406A (zh) | 基于阻抗域属性预测优质储层的方法和装置 | |
| CN110703328B (zh) | 一种超压界面识别方法、装置、设备及存储介质 | |
| CN109143375B (zh) | 一种确定钻井方位的方法及装置 | |
| CN103717833A (zh) | 将测量测井与预测测井相关的方法和系统 | |
| CN110850502A (zh) | 页岩油气储层页理缝分布规律的预测方法、设备及系统 | |
| CN111983671A (zh) | 基于微古地貌恢复的浅水湖盆储层预测方法和装置 | |
| CN111723471B (zh) | 剥蚀厚度确定方法、装置、设备及存储介质 | |
| US9874094B2 (en) | Identifying formation layer boundaries on well log measurements | |
| US11397278B2 (en) | Determining a numerical age for geological events within a scheme | |
| Silversides et al. | BIF-hosted deposit unit differentiation using multivariate Gaussian processes on measure while drilling data | |
| US20230313665A1 (en) | Detection of potential geological control points with logging-while-drilling sensors | |
| US11709287B2 (en) | Well logging to identify low resistivity pay zones in a subterranean formation using elastic attributes | |
| CN114114390B (zh) | 油气勘探方法、装置、设备以及存储介质 | |
| CN107765318A (zh) | 一种确定地质层位的方法及装置 | |
| Guo et al. | Advancement of Drilling Hazard Prevention Practices in HPHT Exploration Well–A Case Study in South China Sea | |
| Lerche et al. | Estimating subsurface fracture density in core: effects resulting from variable fracture spacing | |
| CN114428306A (zh) | 变系数的强反射分离方法、装置、电子设备及介质 | |
| CN118226549A (zh) | 一种砂岩储层流体解释方法、装置、介质及电子设备 | |
| CN118295018A (zh) | 一种分区信息融合速度建模方法、装置、设备和存储介质 | |
| CN112213772A (zh) | 一种井震结合识别地下沉积体形态的方法、设备及系统 | |
| CN118050821A (zh) | 一种基于机器学习的页岩toc含量预测方法和装置 | |
| CN118859304A (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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |