CN112698391A - 砂体刻画方法、装置、设备及存储介质 - Google Patents
砂体刻画方法、装置、设备及存储介质 Download PDFInfo
- Publication number
- CN112698391A CN112698391A CN202011276426.1A CN202011276426A CN112698391A CN 112698391 A CN112698391 A CN 112698391A CN 202011276426 A CN202011276426 A CN 202011276426A CN 112698391 A CN112698391 A CN 112698391A
- Authority
- CN
- China
- Prior art keywords
- target
- sand body
- sand
- seismic data
- depth
- 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
Links
- 239000004576 sand Substances 0.000 title claims abstract description 320
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000003860 storage Methods 0.000 title claims abstract description 14
- 238000005553 drilling Methods 0.000 claims abstract description 17
- 230000008569 process Effects 0.000 claims abstract description 15
- 238000001514 detection method Methods 0.000 claims abstract description 14
- 238000009826 distribution Methods 0.000 claims description 21
- 238000012512 characterization method Methods 0.000 claims description 19
- 230000015654 memory Effects 0.000 claims description 17
- 238000000605 extraction Methods 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 4
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 abstract 2
- 239000003345 natural gas Substances 0.000 abstract 1
- 239000003209 petroleum derivative Substances 0.000 abstract 1
- 238000012545 processing Methods 0.000 description 21
- 230000006870 function Effects 0.000 description 12
- 230000002093 peripheral effect Effects 0.000 description 10
- 230000001133 acceleration Effects 0.000 description 9
- 230000000694 effects Effects 0.000 description 9
- 238000010586 diagram Methods 0.000 description 7
- 238000004891 communication Methods 0.000 description 6
- 238000013507 mapping Methods 0.000 description 6
- 230000003287 optical effect Effects 0.000 description 5
- 230000035945 sensitivity Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000004590 computer program Methods 0.000 description 3
- 238000013473 artificial intelligence Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 239000000919 ceramic Substances 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 208000035126 Facies Diseases 0.000 description 1
- 238000013499 data model Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000005055 memory storage Effects 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
- 239000003208 petroleum Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- 238000013519 translation 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/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/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
技术领域
本申请涉及石油天然气开采领域，特别涉及一种砂体刻画方法、装置、设备及存储介质。
背景技术
砂体刻画是指预测砂体在地层中的分布情况，例如，预测砂体在地层中的位置以及砂体的形状等。随着对石油资源的勘探以及开发程度的不断深入，对砂体刻画精度的要求越来越高。
相关技术中，在钻井过程中检测到砂体后，会根据检测到砂体的位置获取该位置附近区域的地震资料，根据地震资料，将该附近区域中振幅达到预设条件的区域直接作为该砂体所在的区域。
然而，由于对地震资料的分析过于简单，导致砂体刻画的精度较低。
发明内容
本申请实施例提供了一种砂体刻画方法、装置、设备及存储介质，能够提高砂体刻画的精度。所述技术方案如下：
一方面，提供了一种砂体刻画方法，所述方法包括：
根据钻井过程中检测到的砂体的厚度，确定所述砂体对应的敏感频率，所述敏感频率为对所述砂体检测能力最强的地震波的频率；
根据所述敏感频率所在的敏感频率范围，从三维地震数据体中提取多个目标频率对应的分频地震数据体，所述敏感频率范围包括所述多个目标频率；
从多个分频地震数据体中提取所述砂体所在位置处的振幅；
从所述多个分频地震数据体中选择目标地震数据体，所述目标地震数据体对应的振幅与所述砂体的厚度之间的相关度最大；
根据所述砂体顶部距离地面的第一深度，以及所述砂体底部距离所述地面的第二深度，从所述目标地震数据体中截取所述砂体所在的目标区域；
显示所述目标区域的画面，所述目标区域包含所述砂体。
在一种可能的实现方式中，所述显示所述目标区域的画面，包括：
过滤掉所述目标区域中的第一区域，所述第一区域中的任意位置对应的振幅不大于参考振幅；
显示过滤后的目标区域的画面，所述过滤后的目标区域包含所述砂体。
在另一种可能的实现方式中，所述根据所述砂体顶部距离地面的第一深度，以及所述砂体底部距离所述地面的第二深度，从所述目标地震数据体中截取所述砂体所在的目标区域，包括：
将所述第一深度转换成第一时长，所述第一时长表示地震波由地面传播至所述砂体顶部的所用时长；
将所述第二深度转换成第二时长，所述第二时长表示地震波由地面传播至所述砂体底部的所用时长；
根据所述第一时长和所述第二时长，在所述目标地震数据体中标注出所述砂体所在的位置；
从标注后的目标地震数据体中截取所述砂体所在的目标区域。
在另一种可能的实现方式中，所述根据所述砂体顶部距离地面的第一深度，以及所述砂体底部距离所述地面的第二深度，从所述目标地震数据体中截取所述砂体所在的目标区域，包括：
在所述第一深度的基础上减少参考深度，得到第三深度；
在所述第二深度的基础上增加所述参考深度，得到第四深度；
根据所述第三深度和所述第四深度，从所述目标地震数据体中截取所述砂体所在的目标区域。
在另一种可能的实现方式中，所述显示所述目标区域的画面，包括：
对所述目标区域进行种子点追踪，得到目标点集，所述目标点集所在的空间区域表示所述砂体的分布区域；
显示所述目标点集所在的空间区域的画面。
在另一种可能的实现方式中，所述显示所述目标点集所在的空间区域的画面，包括：
确定所述目标点集中的多个目标点对应的时间长度，所述时间长度表示地震波由地面传播至目标点的所用时长；
根据所述多个目标点对应的时间长度，以及地震波由地面传播至所述多个目标点的平均速度，确定所述多个目标点距离所述地面的深度；
根据所述多个目标点距离所述地面的深度，显示所述目标点集所在的空间区域的画面。
另一方面，提供了一种砂体刻画装置，所述装置包括：
频率确定模块，被配置为根据钻井过程中检测到的砂体的厚度，确定所述砂体对应的敏感频率，所述敏感频率为对所述砂体检测能力最强的地震波的频率；
数据体提取模块，被配置为根据所述敏感频率所在的敏感频率范围，从三维地震数据体中提取多个目标频率对应的分频地震数据体，所述敏感频率范围包括所述多个目标频率；
振幅提取模块，被配置为从多个分频地震数据体中提取所述砂体所在位置处的振幅；
数据体选择模块，被配置为从所述多个分频地震数据体中选择目标地震数据体，所述目标地震数据体对应的振幅与所述砂体的厚度之间的相关度最大；
区域截取模块，被配置为根据所述砂体顶部距离地面的第一深度，以及所述砂体底部距离所述地面的第二深度，从所述目标地震数据体中截取所述砂体所在的目标区域；
砂体显示模块，被配置为显示所述目标区域的画面，所述目标区域包含所述砂体。
在一种可能的实现方式中，所述砂体显示模块，被配置为过滤掉所述目标区域中的第一区域，所述第一区域中的任意位置对应的振幅不大于参考振幅；显示过滤后的目标区域的画面，所述过滤后的目标区域包含所述砂体。
在另一种可能的实现方式中，所述区域截取模块，被配置为将所述第一深度转换成第一时长，所述第一时长表示地震波由地面传播至所述砂体顶部的所用时长；将所述第二深度转换成第二时长，所述第二时长表示地震波由地面传播至所述砂体底部的所用时长；根据所述第一时长和所述第二时长，在所述目标地震数据体中标注出所述砂体所在的位置；从标注后的目标地震数据体中截取所述砂体所在的目标区域。
在另一种可能的实现方式中，所述区域截取模块，被配置为在所述第一深度的基础上减少参考深度，得到第三深度；在所述第二深度的基础上增加所述参考深度，得到第四深度；根据所述第三深度和所述第四深度，从所述目标地震数据体中截取所述砂体所在的目标区域。
在另一种可能的实现方式中，所述砂体显示模块，包括：
种子点追踪单元，被配置为对所述目标区域进行种子点追踪，得到目标点集，所述目标点集所在的空间区域表示所述砂体的分布区域；
画面显示单元，被配置为显示所述目标点集所在的空间区域的画面。
在另一种可能的实现方式中，所述画面显示单元，被配置为确定所述目标点集中的多个目标点对应的时间长度，所述时间长度表示地震波由地面传播至目标点的所用时长；根据所述多个目标点对应的时间长度，以及地震波由地面传播至所述多个目标点的平均速度，确定所述多个目标点距离所述地面的深度；根据所述多个目标点距离所述地面的深度，显示所述目标点集所在的空间区域的画面。
另一方面，提供了一种计算机设备，所述计算机设备包括处理器和存储器，所述存储器中存储有至少一条程序代码，所述程序代码由所述处理器加载并执行以实现上述任一种可能实现方式中的砂体刻画方法中执行的操作。
另一方面，提供了一种计算机可读存储介质，所述计算机可读存储介质中存储有至少一条程序代码，所述程序代码由处理器加载并执行以实现上述任一种可能实现方式中的砂体刻画方法中执行的操作。
再一方面，提供了一种计算机程序产品，所述计算机程序产品中包括至少一条程序代码，所述程序代码由处理器加载并执行以实现上述任一种可能实现方式中的砂体刻画方法中执行的操作。
本申请实施例提供的技术方案带来的有益效果至少包括：
本申请实施例所提供的技术方案，考虑到不同频率对不同厚度砂体的敏感程度不同，属于敏感频率的地震波对砂体的检测能力最强，对砂体分布的表现效果更好，因此利用钻井过程中检测到的砂体的厚度，获取砂体对应的敏感频率，根据敏感频率提取分频地震数据体，则提取的分频地震数据体能够很好地表现该砂体的分布。并且，考虑到获取的敏感频率可能具有误差，因此在根据敏感频率提取分频地震数据体时，根据敏感频率所在的敏感频率范围，从三维地震数据体中提取敏感频率范围中的多个目标频率对应的分频地震数据体，以避免敏感频率的误差导致获取的分频地震数据体不准确。并且，考虑到振幅与厚度砂体的相关度越大的三维地震数据体对砂体的识别能力越强，从提取的多个分频地震数据体中选择振幅与砂体厚度的相关度最大的目标地震数据体，进一步保证了目标地震数据体对砂体分布的表现效果，后续从该目标地震数据体中截取砂体所在的目标区域，显示目标区域的画面，能够保证画面中砂体位置和形状的准确性，提高砂体的刻画精度。
附图说明
为了更清楚地说明本申请实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本申请实施例提供的一种砂体刻画方法的流程图；
图2是本申请实施例提供的一种砂体刻画方法的流程图；
图3是本申请实施例提供的一种标注砂体位置的目标地震数据体的示意图；
图4是本申请实施例提供的一种刻画出的砂体的示意图；
图5是本申请实施例提供的一种砂体刻画装置的框图；
图6是本申请实施例提供的一种终端的结构示意图；
图7是本申请实施例提供的一种服务器的结构示意图。
具体实施方式
为使本申请的目的、技术方案和优点更加清楚，下面将结合附图对本申请实施方式作进一步地详细描述。
本申请所使用的术语“第一”、“第二”、“第三”、“第四”等可在本文中用于描述各种概念，但除非特别说明，这些概念不受这些术语限制。这些术语仅用于将一个概念与另一个概念区分。举例来说，在不脱离本申请的范围的情况下，可以将第一深度称为深度，且类似地，可将第二深度称为第一深度。
本申请所使用的术语“至少一个”、“多个”、“每个”、“任一”，至少一个包括一个、两个或两个以上，多个包括两个或两个以上，而每个是指对应的多个中的每一个，任一是指多个中的任意一个。举例来说，多个目标频率包括3个目标频率，而每个是指这3个目标频率中的每一个目标频率，任一是指这3个目标频率中的任意一个，可以是第一个，可以是第二个、也可以是第三个。
图1是本申请实施例提供的一种砂体刻画方法的流程图。执行主体为计算机设备，参见图1，该实施例包括：
101：根据钻井过程中检测到的砂体的厚度，确定砂体对应的敏感频率，敏感频率为对砂体检测能力最强的地震波的频率。
102：根据敏感频率所在的敏感频率范围，从三维地震数据体中提取多个目标频率对应的分频地震数据体，敏感频率范围包括多个目标频率。
103：从多个分频地震数据体中提取砂体所在位置处的振幅。
104：从多个分频地震数据体中选择目标地震数据体，目标地震数据体对应的振幅与砂体的厚度之间的相关度最大。
105：根据砂体顶部距离地面的第一深度，以及砂体底部距离地面的第二深度，从目标地震数据体中截取砂体所在的目标区域。
106：显示目标区域的画面，目标区域包含砂体。
本申请实施例所提供的技术方案，考虑到不同频率对不同厚度砂体的敏感程度不同，属于敏感频率的地震波对砂体的检测能力最强，对砂体分布的表现效果更好，因此利用钻井过程中检测到的砂体的厚度，获取砂体对应的敏感频率，根据敏感频率提取分频地震数据体，则提取的分频地震数据体能够很好地表现该砂体的分布。并且，考虑到获取的敏感频率可能具有误差，因此在根据敏感频率提取分频地震数据体时，根据敏感频率所在的敏感频率范围，从三维地震数据体中提取敏感频率范围中的多个目标频率对应的分频地震数据体，以避免敏感频率的误差导致获取的分频地震数据体不准确。并且，考虑到振幅与厚度砂体的相关度越大的三维地震数据体对砂体的识别能力越强，从提取的多个分频地震数据体中选择振幅与砂体厚度的相关度最大的目标地震数据体，进一步保证了目标地震数据体对砂体分布的表现效果，后续从该目标地震数据体中截取砂体所在的目标区域，显示目标区域的画面，能够保证画面中砂体位置和形状的准确性，提高砂体的刻画精度。
在一种可能的实现方式中，显示目标区域的画面，包括：
过滤掉目标区域中的第一区域，第一区域中的任意位置对应的振幅不大于参考振幅；
显示过滤后的目标区域的画面，过滤后的目标区域包含砂体。
在另一种可能的实现方式中，根据砂体顶部距离地面的第一深度，以及砂体底部距离地面的第二深度，从目标地震数据体中截取砂体所在的目标区域，包括：
将第一深度转换成第一时长，第一时长表示地震波由地面传播至砂体顶部的所用时长；
将第二深度转换成第二时长，第二时长表示地震波由地面传播至砂体底部的所用时长；
根据第一时长和第二时长，在目标地震数据体中标注出砂体所在的位置；
从标注后的目标地震数据体中截取砂体所在的目标区域。
在另一种可能的实现方式中，根据砂体顶部距离地面的第一深度，以及砂体底部距离地面的第二深度，从目标地震数据体中截取砂体所在的目标区域，包括：
在第一深度的基础上减少参考深度，得到第三深度；
在第二深度的基础上增加参考深度，得到第四深度；
根据第三深度和第四深度，从目标地震数据体中截取砂体所在的目标区域。
在另一种可能的实现方式中，显示目标区域的画面，包括：
对目标区域进行种子点追踪，得到目标点集，目标点集所在的空间区域表示砂体的分布区域；
显示目标点集所在的空间区域的画面。
在另一种可能的实现方式中，显示目标点集所在的空间区域的画面，包括：
确定目标点集中的多个目标点对应的时间长度，时间长度表示地震波由地面传播至目标点的所用时长；
根据多个目标点对应的时间长度，以及地震波由地面传播至多个目标点的平均速度，确定多个目标点距离地面的深度；
根据多个目标点距离地面的深度，显示目标点集所在的空间区域的画面。
上述所有可选技术方案，可以采用任意结合形成本申请的可选实施例，在此不再一一赘述。
图2是本申请实施例提供的一种砂体刻画方法的流程图。参见图2，该实施例包括：
201：计算机设备根据钻井过程中检测到的砂体的厚度，确定砂体对应的敏感频率，敏感频率为对砂体检测能力最强的地震波的频率。
地震波具有不同的频率，不同频率的地震波对不同厚度的砂体的检测能力不同。对于特定厚度的砂体来说，对该砂体检测能力最强的地震波的频率即为该砂体对应的敏感频率。
钻井过程中在钻遇砂体时，能够直接检测该砂体的厚度，计算机设备则能够根据该砂体厚度来确定砂体对应的敏感频率。
在一种可能的实现方式中，计算机设备存储砂体厚度和敏感频率之间的映射关系，根据检测的砂体的厚度和该映射关系，确定该砂体对应的敏感频率。
可选地，该映射关系为如下述公式(1)：
其中，f表示敏感频率，v表示砂体所在地层的层速度，h表示砂体的厚度。
例如，砂体所在地层的层速度为2000米/秒，砂岩厚度为10米时，则根据该公式
确定砂体的敏感频率为50Hz。
需要说明的一点是，若通过多个钻井检测到该砂体，且多个钻井检测的砂体的厚度不同，能够取多个厚度的平均值作为该砂体的厚度，或者将多个厚度中最大的厚度作为砂体的厚度。
202：计算机设备根据敏感频率所在的敏感频率范围，从三维地震数据体中提取多个目标频率对应的分频地震数据体，敏感频率范围包括多个目标频率。
三维地震数据体是一个三维的地震数据模型，其中具有多个高密度的信息点，每个信息点具有丰富的地震信息，能够反映地下的地质构造情况。
示例性的，通过人工制造地震波方法，采集地震信息，将该地震信息输入到计算机设备，按不同要求用一系列功能不同的程序进行处理运算，把地震信息进行归类编排，突出有效的地震信息，除去无效和干扰的信息，最后把经过各种处理的地震信息进行叠加和偏移，生成三维地震数据体。
分频地震数据体是对三维地震数据体进行地震分频处理后得到的地震数据体。对三维地震数据体进行地震分频处理，能够得到对应不同频率的分频地震数据体。
其中，地震分频处理是利用不同的地层厚度对应不同频率的地震成像解释方法，能够整体揭示地层的纵向变化规律、沉积相带的空间演变模式。经过地震分频处理后的分频地震数据体，波形特征十分清晰，对砂体的空间变化极度敏感。
计算机设备确定敏感频率后，还要确定敏感频率所在的敏感频率范围。在一种可能的实现方式中，计算机设备将敏感频率增加参考频率得到第一频率，将敏感频率减少参考频率得到第二频率，则第一频率、第二频率以及第一频率和第二频率之间的频率则构成敏感频率范围。
从三维地震数据体中提取多个目标频率对应的分频地震数据体之前，计算机设备先要从敏感频率范围中确定目标频率。可选地，计算机设备将敏感频率范围中的每个频率作为目标频率，或者，计算机设备按照参考频率间隔从敏感频率范围中提取多个目标频率。如此，能够减少提取的分频地震数据体的数据量，降低计算成本。例如，参考频率范围为【40Hz-60Hz】，参考频率间隔为5Hz，则提取的目标频率分别为40Hz、45Hz、50Hz、55Hz、60Hz。
可选地，从三维地震数据体中提取多个目标频率对应的分频地震数据体的实现方式包括：通过分频处理函数对三维地震数据体进行分频处理后，选择目标频率对应的分频地震数据体。
可选地，分频处理函数如下述公式(2)：
其中，a为尺度因子，b为平移因子，f(t)为实信号，t为时间，ψ(t)为小波变换的母小波，上面的短横“—”表示取负数共轭，(wψf)(a，b)为小波变换的系数。
需要说明的一点是，上述分频处理函数采用的是小波变换，由于小波变换在时间域和频率域中都具有表征信号局部变化的能力，对信号的动态瞬时分析十分有利。因此，得到的分频地震数据体更加准确。当然，也能够采用离散傅里叶变换、S变换等其他分频处理函数，本申请实施例对此不做限制。
203：计算机设备从多个分频地震数据体中提取砂体所在位置处的振幅。
可选地，计算机设备对每个分频地震数据体进行振幅属性分析，得到每个分频地震数据体中砂体所在位置处的振幅。
204：计算机设备从多个分频地震数据体中选择目标地震数据体，目标地震数据体对应的振幅与砂体的厚度之间的相关度最大。
计算机设备从多个分频地震数据体中选择目标地震数据体之前，先要确定每个分频地震数据体与砂体的厚度之间的相关度。
可选地，计算机设备存储分频地震数据体与砂体厚度的相关度与砂体的厚度之间的映射关系，根据该映射关系和砂体的厚度确定每个分频地震数据体对应的相关度。
可选地，该映射关系为下述公式(3)。
其中，r为相关系数，n为已钻井数，x为砂体厚度，y为砂体所在位置处的振幅值。
在本申请实施例中，考虑到振幅与厚度砂体的相关度越大的分频地震数据体对砂体的识别能力越强，从提取的多个分频地震数据体中选择振幅与砂体厚度的相关度最大的目标地震数据体，进一步保证了目标地震数据体对砂体分布的表现效果。
205：计算机设备根据砂体顶部距离地面的第一深度，以及砂体底部距离地面的第二深度，从目标地震数据体中截取砂体所在的目标区域。
可选地，钻井过程中在钻遇砂体时，直接检测该砂体顶部距离地面的第一深度，以及砂体底部距离地面的第二深度。
在本申请实施例中，根据砂体顶部距离地面的第一深度，以及砂体底部距离地面的第二深度，从目标地震数据体中截取砂体所在的目标区域，能够保证后续显示砂体所在目标区域的画面时，显示的画面中砂体所占区域的比例较大，从而能够提高砂体刻画的精度。
需要说明的一点是，由于目标地震数据体是属于时间域的数据体，在目标地震数据体中，对于任意一个信息点，通过地震波从地面到达该信息点的时长以及地震波从地面到达该信息点所在地层的平均速度来表示该信息点距离地面的深度。因此，在从目标地震数据体中截取砂体所在的目标区域前，要先将深度转换为时长，根据时长从目标地震数据体中截取砂体所在的目标区域。
相应的，计算机设备根据砂体顶部距离地面的第一深度，以及砂体底部距离地面的第二深度，从目标地震数据体中截取砂体所在的目标区域，包括：计算机设备将第一深度转换成第一时长，第一时长表示地震波由地面传播至砂体顶部的所用时长；将第二深度转换成第二时长，第二时长表示地震波由地面传播至砂体底部的所用时长；计算机设备根据第一时长和第二时长，在目标地震数据体中标注出砂体所在的位置；从标注后的目标地震数据体中截取砂体所在的目标区域。
在一种可能的实现方式中，计算机设备根据砂体顶部距离地面的第一深度，以及砂体底部距离地面的第二深度，从目标地震数据体中截取砂体所在的目标区域，包括：计算机设备在第一深度的基础上减少参考深度，得到第三深度；在第二深度的基础上增加参考深度，得到第四深度；根据第三深度和第四深度，从目标地震数据体中截取砂体所在的目标区域。可选地，参考深度为任意数值，本申请实施例对此不做限制。可选地，该参考深度根据该砂体所在地层上部的泥层的厚度来确定。例如，将该参考深度设置为小于该泥层厚度的任意数值。
在本申请实施例中，考虑到砂体的各部分的厚度可能不同，直接根据第一深度和第二深度来截取目标区域，可能不能将砂体所在目标区域完整地截取出来，因此，在第一深度的基础上减少参考深度，得到第三深度，在第二深度的基础上增加参考深度，得到第四深度，根据第三深度和第四深度进行截取，以保证能够将砂体所在区域完整地截取出来，从而保证砂体刻画的准确性。
需要说明的另一点是，若通过多个钻井检测到该砂体，则根据每个钻井检测的砂体顶部的深度和砂体底部的深度，从目标地震数据体中截取砂体所在的目标区域。实现方式为：计算机设备根据多个钻井对应的砂体顶部的深度和砂体底部的深度，在目标地震数据体中标注砂体所在的多个位置，然后从标注后的目标地震数据体中截取砂体所在的目标区域。参考图3，为标注砂体位置后的目标地震数据体的示意图。其中，水平方向上的具有一定宽度的线条表示砂体所在地层的走向。
需要说明的再一点是，在从标注后的目标地震数据体中截取砂体所在的目标区域时，要沿着砂体所在地层的走向来截取。如此，能够保证截取的砂体所在区域的完整性，从而保证砂体刻画的精度。
206：计算机设备过滤掉目标区域中的第一区域，第一区域中的任意位置对应的振幅不大于参考振幅。
其中，参考振幅为砂体和泥岩的振幅分界线，不同地区对应的振幅分界线可能不同，可选地，根据砂体所在地区确定对应的振幅分界线作为参考振幅。例如，参考振幅为5000。由于低于参考振幅的第一区域为泥岩所在区域，因此，过滤掉目标区域中的第一区域，能够避免泥岩对砂体刻画的影响，提高砂体刻画的精度。
207：计算机设备显示过滤后的目标区域的画面，过滤后的目标区域包含砂体。
在一种可能的实现方式中，计算机设备显示目标区域的画面，包括：计算机设备对目标区域进行种子点追踪，得到目标点集，目标点集所在的空间区域表示砂体的分布区域；计算机设备显示目标点集所在的空间区域的画面。
可选地，对目标区域进行种子点追踪的实现方式为：计算机设备从目标区域中选取信息点作为种子点，根据种子点周围信息点的地震信息与种子点的地震信息的相似度，从种子点周围的信息点中选取新的种子点，然后再根据新选取的种子点继续选取种子点。其中，可选地，根据种子点周围信息点的地震信息与种子点的地震信息的相似度，从种子点周围的信息点中选取新的种子点的实现方式为：将种子点周围信息点中，地震信息与种子点的地震信息的相似度小于参考相似度的信息点作为新的种子点。其中，参考相似度根据需要设置为任意数值，本申请实施例对此不做限制。
计算机设备对目标区域进行种子点追踪之后，将追踪到的所有种子点构成目标点集，其中，目标点集中的种子点即为目标点。
在本申请实施例中，对目标区域进行种子点追踪，得到目标点集，显示目标点集所在的空间区域的画面，在点的层面上去刻画砂体的位置和形态，因此画面中砂体的位置和形状更加准确，提高了砂体刻画的精度。
在一种可能的实现方式中，计算机设备显示目标点集所在的空间区域的画面，包括：计算机设备确定目标点集中的多个目标点对应的时间长度，时间长度表示地震波由地面传播至目标点的所用时长；根据多个目标点对应的时间长度，以及地震波由地面传播至多个目标点的平均速度，确定多个目标点距离地面的深度；根据多个目标点距离地面的深度，显示目标点集所在的空间区域的画面。
其中，根据多个目标点对应的时间长度，以及地震波由地面传播至多个目标点的平均速度，确定多个目标点距离地面的深度的实现方式为：将每个目标点对应的时间长度与地震波由地面传播至该目标点的平均速度的乘积，确定为该目标点距离地面的深度。由于目标点集中的目标点是通过地震波从地面传播至该目标点的时间长度来表示的，因此，将时间长度先转换成深度，然后再显示目标点集所在的空间区域的画面，能够更加直观地体现砂体的深度和厚度，方便用户查看，从而提高用户粘性。参考图4，为刻画出的砂体的示意图。其中，各个数值表示砂体各部分对应的深度。
需要说明的一点是，上述206是可选地，也即是，能够不对目标区域进行过滤处理，而直接显示目标区域的画面，如此能够减少计算机设备的计算量，提高砂体刻画的效率。
本申请实施例所提供的技术方案，考虑到不同频率对不同厚度砂体的敏感程度不同，属于敏感频率的地震波对砂体的检测能力最强，对砂体分布的表现效果更好，因此利用钻井过程中检测到的砂体的厚度，获取砂体对应的敏感频率，根据敏感频率提取分频地震数据体，则提取的分频地震数据体能够很好地表现该砂体的分布。并且，考虑到获取的敏感频率可能具有误差，因此在根据敏感频率提取分频地震数据体时，根据敏感频率所在的敏感频率范围，从三维地震数据体中提取敏感频率范围中的多个目标频率对应的分频地震数据体，以避免敏感频率的误差导致获取的分频地震数据体不准确。并且，考虑到振幅与厚度砂体的相关度越大的三维地震数据体对砂体的识别能力越强，从提取的多个分频地震数据体中选择振幅与砂体厚度的相关度最大的目标地震数据体，进一步保证了目标地震数据体对砂体分布的表现效果，后续从该目标地震数据体中截取砂体所在的目标区域，显示目标区域的画面，能够保证画面中砂体位置和形状的准确性，提高砂体的刻画精度。
图5是本申请实施例提供的一种砂体刻画方法的流程图。参见图5，该装置包括：
频率确定模块501，被配置为根据钻井过程中检测到的砂体的厚度，确定砂体对应的敏感频率，敏感频率为对砂体检测能力最强的地震波的频率；
数据体提取模块502，被配置为根据敏感频率所在的敏感频率范围，从三维地震数据体中提取多个目标频率对应的分频地震数据体，敏感频率范围包括多个目标频率；
振幅提取模块503，被配置为从多个分频地震数据体中提取砂体所在位置处的振幅；
数据体选择模块504，被配置为从多个分频地震数据体中选择目标地震数据体，目标地震数据体对应的振幅与砂体的厚度之间的相关度最大；
区域截取模块505，被配置为根据砂体顶部距离地面的第一深度，以及砂体底部距离地面的第二深度，从目标地震数据体中截取砂体所在的目标区域；
砂体显示模块506，被配置为显示目标区域的画面，目标区域包含砂体。
在一种可能的实现方式中，砂体显示模块506，被配置为过滤掉目标区域中的第一区域，第一区域中的任意位置对应的振幅不大于参考振幅；显示过滤后的目标区域的画面，过滤后的目标区域包含砂体。
在另一种可能的实现方式中，区域截取模块505，被配置为将第一深度转换成第一时长，第一时长表示地震波由地面传播至砂体顶部的所用时长；将第二深度转换成第二时长，第二时长表示地震波由地面传播至砂体底部的所用时长；根据第一时长和第二时长，在目标地震数据体中标注出砂体所在的位置；从标注后的目标地震数据体中截取砂体所在的目标区域。
在另一种可能的实现方式中，区域截取模块505，被配置为在第一深度的基础上减少参考深度，得到第三深度；在第二深度的基础上增加参考深度，得到第四深度；根据第三深度和第四深度，从目标地震数据体中截取砂体所在的目标区域。
在另一种可能的实现方式中，砂体显示模块506，包括：
种子点追踪单元，被配置为对目标区域进行种子点追踪，得到目标点集，目标点集所在的空间区域表示砂体的分布区域；
画面显示单元，被配置为显示目标点集所在的空间区域的画面。
在另一种可能的实现方式中，画面显示单元，被配置为确定目标点集中的多个目标点对应的时间长度，时间长度表示地震波由地面传播至目标点的所用时长；根据多个目标点对应的时间长度，以及地震波由地面传播至多个目标点的平均速度，确定多个目标点距离地面的深度；根据多个目标点距离地面的深度，显示目标点集所在的空间区域的画面。
本申请实施例所提供的技术方案，考虑到不同频率对不同厚度砂体的敏感程度不同，属于敏感频率的地震波对砂体的检测能力最强，对砂体分布的表现效果更好，因此利用钻井过程中检测到的砂体的厚度，获取砂体对应的敏感频率，根据敏感频率提取分频地震数据体，则提取的分频地震数据体能够很好地表现该砂体的分布。并且，考虑到获取的敏感频率可能具有误差，因此在根据敏感频率提取分频地震数据体时，根据敏感频率所在的敏感频率范围，从三维地震数据体中提取敏感频率范围中的多个目标频率对应的分频地震数据体，以避免敏感频率的误差导致获取的分频地震数据体不准确。并且，考虑到振幅与厚度砂体的相关度越大的三维地震数据体对砂体的识别能力越强，从提取的多个分频地震数据体中选择振幅与砂体厚度的相关度最大的目标地震数据体，进一步保证了目标地震数据体对砂体分布的表现效果，后续从该目标地震数据体中截取砂体所在的目标区域，显示目标区域的画面，能够保证画面中砂体位置和形状的准确性，提高砂体的刻画精度。
需要说明的是：上述实施例提供的砂体刻画装置在进行砂体刻画时，仅以上述各功能模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能模块完成，即将计算机设备的内部结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。另外，上述实施例提供的砂体刻画装置与砂体刻画方法实施例属于同一构思，其具体实现过程详见方法实施例，这里不再赘述。
本申请实施例还提供了一种计算机设备，该计算机设备包括处理器和存储器，存储器中存储有至少一条程序代码，该至少一条程序代码由处理器加载并执行，以实现上述实施例的砂体刻画方法中执行的操作。
可选地，该计算机设备提供为终端。图6示出了本申请一个示例性实施例提供的终端600的结构框图。该终端600可以是：智能手机、平板电脑、MP3播放器(Moving PictureExperts Group Audio Layer III，动态影像专家压缩标准音频层面3)、MP4(MovingPicture Experts Group Audio Layer IV，动态影像专家压缩标准音频层面4)播放器、笔记本电脑或台式电脑。终端600还可能被称为用户设备、便携式终端、膝上型终端、台式终端等其他名称。
终端600包括有：处理器601和存储器602。
处理器601可以包括一个或多个处理核心，比如4核心处理器、8核心处理器等。处理器601可以采用DSP(Digital Signal Processing，数字信号处理)、FPGA(Field－Programmable Gate Array，现场可编程门阵列)、PLA(Programmable Logic Array，可编程逻辑阵列)中的至少一种硬件形式来实现。处理器601也可以包括主处理器和协处理器，主处理器是用于对在唤醒状态下的数据进行处理的处理器，也称CPU(Central ProcessingUnit，中央处理器)；协处理器是用于对在待机状态下的数据进行处理的低功耗处理器。在一些实施例中，处理器601可以集成有GPU(Graphics Processing Unit，图像处理器)，GPU用于负责显示屏所需要显示的内容的渲染和绘制。一些实施例中，处理器601还可以包括AI(Artificial Intelligence，人工智能)处理器，该AI处理器用于处理有关机器学习的计算操作。
存储器602可以包括一个或多个计算机可读存储介质，该计算机可读存储介质可以是非暂态的。存储器602还可包括高速随机存取存储器，以及非易失性存储器，比如一个或多个磁盘存储设备、闪存存储设备。在一些实施例中，存储器602中的非暂态的计算机可读存储介质用于存储至少一个程序代码，该至少一个程序代码用于被处理器601所执行以实现本申请中方法实施例提供的砂体刻画方法。
在一些实施例中，终端600还可选包括有：外围设备接口603和至少一个外围设备。处理器601、存储器602和外围设备接口603之间可以通过总线或信号线相连。各个外围设备可以通过总线、信号线或电路板与外围设备接口603相连。具体地，外围设备包括：射频电路604、显示屏605、摄像头组件606、音频电路607、定位组件608和电源609中的至少一种。
外围设备接口603可被用于将I/O(Input/Output，输入/输出)相关的至少一个外围设备连接到处理器601和存储器602。在一些实施例中，处理器601、存储器602和外围设备接口603被集成在同一芯片或电路板上；在一些其他实施例中，处理器601、存储器602和外围设备接口603中的任意一个或两个可以在单独的芯片或电路板上实现，本实施例对此不加以限定。
射频电路604用于接收和发射RF(Radio Frequency，射频)信号，也称电磁信号。射频电路604通过电磁信号与通信网络以及其他通信设备进行通信。射频电路604将电信号转换为电磁信号进行发送，或者，将接收到的电磁信号转换为电信号。可选地，射频电路604包括：天线系统、RF收发器、一个或多个放大器、调谐器、振荡器、数字信号处理器、编解码芯片组、用户身份模块卡等等。射频电路604可以通过至少一种无线通信协议来与其它终端进行通信。该无线通信协议包括但不限于：城域网、各代移动通信网络(2G、3G、4G及5G)、无线局域网和/或WiFi(Wireless Fidelity，无线保真)网络。在一些实施例中，射频电路604还可以包括NFC(Near Field Communication，近距离无线通信)有关的电路，本申请对此不加以限定。
显示屏605用于显示UI(User Interface，用户界面)。该UI可以包括图形、文本、图标、视频及其它们的任意组合。当显示屏605是触摸显示屏时，显示屏605还具有采集在显示屏605的表面或表面上方的触摸信号的能力。该触摸信号可以作为控制信号输入至处理器601进行处理。此时，显示屏605还可以用于提供虚拟按钮和/或虚拟键盘，也称软按钮和/或软键盘。在一些实施例中，显示屏605可以为一个，设置终端600的前面板；在另一些实施例中，显示屏605可以为至少两个，分别设置在终端600的不同表面或呈折叠设计；在另一些实施例中，显示屏605可以是柔性显示屏，设置在终端600的弯曲表面上或折叠面上。甚至，显示屏605还可以设置成非矩形的不规则图形，也即异形屏。显示屏605可以采用LCD(LiquidCrystal Display，液晶显示屏)、OLED(Organic Light-Emitting Diode,有机发光二极管)等材质制备。
摄像头组件606用于采集图像或视频。可选地，摄像头组件606包括前置摄像头和后置摄像头。前置摄像头设置在终端的前面板，后置摄像头设置在终端的背面。在一些实施例中，后置摄像头为至少两个，分别为主摄像头、景深摄像头、广角摄像头、长焦摄像头中的任意一种，以实现主摄像头和景深摄像头融合实现背景虚化功能、主摄像头和广角摄像头融合实现全景拍摄以及VR(Virtual Reality，虚拟现实)拍摄功能或者其它融合拍摄功能。在一些实施例中，摄像头组件606还可以包括闪光灯。闪光灯可以是单色温闪光灯，也可以是双色温闪光灯。双色温闪光灯是指暖光闪光灯和冷光闪光灯的组合，可以用于不同色温下的光线补偿。
音频电路607可以包括麦克风和扬声器。麦克风用于采集用户及环境的声波，并将声波转换为电信号输入至处理器601进行处理，或者输入至射频电路604以实现语音通信。出于立体声采集或降噪的目的，麦克风可以为多个，分别设置在终端600的不同部位。麦克风还可以是阵列麦克风或全向采集型麦克风。扬声器则用于将来自处理器601或射频电路604的电信号转换为声波。扬声器可以是传统的薄膜扬声器，也可以是压电陶瓷扬声器。当扬声器是压电陶瓷扬声器时，不仅可以将电信号转换为人类可听见的声波，也可以将电信号转换为人类听不见的声波以进行测距等用途。在一些实施例中，音频电路607还可以包括耳机插孔。
定位组件608用于定位终端600的当前地理位置，以实现导航或LBS(LocationBased Service，基于位置的服务)。定位组件608可以是基于美国的GPS(GlobalPositioning System，全球定位系统)、中国的北斗系统、俄罗斯的格雷纳斯系统或欧盟的伽利略系统的定位组件。
电源609用于为终端600中的各个组件进行供电。电源609可以是交流电、直流电、一次性电池或可充电电池。当电源609包括可充电电池时，该可充电电池可以支持有线充电或无线充电。该可充电电池还可以用于支持快充技术。
在一些实施例中，终端600还包括有一个或多个传感器610。该一个或多个传感器610包括但不限于：加速度传感器611、陀螺仪传感器612、压力传感器613、指纹传感器614、光学传感器615以及接近传感器616。
加速度传感器611可以检测以终端600建立的坐标系的三个坐标轴上的加速度大小。比如，加速度传感器611可以用于检测重力加速度在三个坐标轴上的分量。处理器601可以根据加速度传感器611采集的重力加速度信号，控制显示屏605以横向视图或纵向视图进行用户界面的显示。加速度传感器611还可以用于游戏或者用户的运动数据的采集。
陀螺仪传感器612可以检测终端600的机体方向及转动角度，陀螺仪传感器612可以与加速度传感器611协同采集用户对终端600的3D动作。处理器601根据陀螺仪传感器612采集的数据，可以实现如下功能：动作感应(比如根据用户的倾斜操作来改变UI)、拍摄时的图像稳定、游戏控制以及惯性导航。
压力传感器613可以设置在终端600的侧边框和/或显示屏605的下层。当压力传感器613设置在终端600的侧边框时，可以检测用户对终端600的握持信号，由处理器601根据压力传感器613采集的握持信号进行左右手识别或快捷操作。当压力传感器613设置在显示屏605的下层时，由处理器601根据用户对显示屏605的压力操作，实现对UI界面上的可操作性控件进行控制。可操作性控件包括按钮控件、滚动条控件、图标控件、菜单控件中的至少一种。
指纹传感器614用于采集用户的指纹，由处理器601根据指纹传感器614采集到的指纹识别用户的身份，或者，由指纹传感器614根据采集到的指纹识别用户的身份。在识别出用户的身份为可信身份时，由处理器601授权该用户执行相关的敏感操作，该敏感操作包括解锁屏幕、查看加密信息、下载软件、支付及更改设置等。指纹传感器614可以被设置终端600的正面、背面或侧面。当终端600上设置有物理按键或厂商Logo时，指纹传感器614可以与物理按键或厂商Logo集成在一起。
光学传感器615用于采集环境光强度。在一个实施例中，处理器601可以根据光学传感器615采集的环境光强度，控制显示屏605的显示亮度。具体地，当环境光强度较高时，调高显示屏605的显示亮度；当环境光强度较低时，调低显示屏605的显示亮度。在另一个实施例中，处理器601还可以根据光学传感器615采集的环境光强度，动态调整摄像头组件606的拍摄参数。
接近传感器616，也称距离传感器，设置在终端600的前面板。接近传感器616用于采集用户与终端600的正面之间的距离。在一个实施例中，当接近传感器616检测到用户与终端600的正面之间的距离逐渐变小时，由处理器601控制显示屏605从亮屏状态切换为息屏状态；当接近传感器616检测到用户与终端600的正面之间的距离逐渐变大时，由处理器601控制显示屏605从息屏状态切换为亮屏状态。
本领域技术人员可以理解，图6中示出的结构并不构成对终端600的限定，可以包括比图示更多或更少的组件，或者组合某些组件，或者采用不同的组件布置。
可选地，该计算机设备提供为服务器。图7是本申请实施例提供的一种服务器的结构示意图，该服务器700可因配置或性能不同而产生比较大的差异，可以包括一个或一个以上处理器(central processing units，CPU)701和一个或一个以上的存储器702，其中，所述存储器702中存储有至少一条程序代码，所述至少一条程序代码由所述处理器701加载并执行以实现上述各个方法实施例提供的砂体刻画方法。当然，该服务器还可以具有有线或无线网络接口、键盘以及输入输出接口等部件，以便进行输入输出，该服务器还可以包括其他用于实现设备功能的部件，在此不做赘述。
本申请实施例还提供了一种计算机可读存储介质，该计算机可读存储介质中存储有至少一条程序代码，该至少一条程序代码由处理器加载并执行，以实现上述实施例的砂体刻画方法中执行的操作。
本申请实施例还提供了一种计算机程序，该计算机程序中存储有至少一条程序代码，该至少一条程序代码由处理器加载并执行，以实现上述实施例的砂体刻画方法中执行的操作。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成，也可以通过程序来指令相关的硬件完成，所述的程序可以存储于一种计算机可读存储介质中，上述提到的存储介质可以是只读存储器，磁盘或光盘等。
以上所述仅为本申请的可选实施例，并不用以限制本申请，凡在本申请的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。
Claims (10)
1.一种砂体刻画方法，其特征在于，所述方法包括：
根据钻井过程中检测到的砂体的厚度，确定所述砂体对应的敏感频率，所述敏感频率为对所述砂体检测能力最强的地震波的频率；
根据所述敏感频率所在的敏感频率范围，从三维地震数据体中提取多个目标频率对应的分频地震数据体，所述敏感频率范围包括所述多个目标频率；
从多个分频地震数据体中提取所述砂体所在位置处的振幅；
从所述多个分频地震数据体中选择目标地震数据体，所述目标地震数据体对应的振幅与所述砂体的厚度之间的相关度最大；
根据所述砂体顶部距离地面的第一深度，以及所述砂体底部距离所述地面的第二深度，从所述目标地震数据体中截取所述砂体所在的目标区域；
显示所述目标区域的画面，所述目标区域包含所述砂体。
2.根据权利要求1所述的方法，其特征在于，所述显示所述目标区域的画面，包括：
过滤掉所述目标区域中的第一区域，所述第一区域中的任意位置对应的振幅不大于参考振幅；
显示过滤后的目标区域的画面，所述过滤后的目标区域包含所述砂体。
3.根据权利要求1所述的方法，其特征在于，所述根据所述砂体顶部距离地面的第一深度，以及所述砂体底部距离所述地面的第二深度，从所述目标地震数据体中截取所述砂体所在的目标区域，包括：
将所述第一深度转换成第一时长，所述第一时长表示地震波由地面传播至所述砂体顶部的所用时长；
将所述第二深度转换成第二时长，所述第二时长表示地震波由地面传播至所述砂体底部的所用时长；
根据所述第一时长和所述第二时长，在所述目标地震数据体中标注出所述砂体所在的位置；
从标注后的目标地震数据体中截取所述砂体所在的目标区域。
4.根据权利要求1所述的方法，其特征在于，所述根据所述砂体顶部距离地面的第一深度，以及所述砂体底部距离所述地面的第二深度，从所述目标地震数据体中截取所述砂体所在的目标区域，包括：
在所述第一深度的基础上减少参考深度，得到第三深度；
在所述第二深度的基础上增加所述参考深度，得到第四深度；
根据所述第三深度和所述第四深度，从所述目标地震数据体中截取所述砂体所在的目标区域。
5.根据权利要求1所述的方法，其特征在于，所述显示所述目标区域的画面，包括：
对所述目标区域进行种子点追踪，得到目标点集，所述目标点集所在的空间区域表示所述砂体的分布区域；
显示所述目标点集所在的空间区域的画面。
6.根据权利要求5所述的方法，其特征在于，所述显示所述目标点集所在的空间区域的画面，包括：
确定所述目标点集中的多个目标点对应的时间长度，所述时间长度表示地震波由地面传播至目标点的所用时长；
根据所述多个目标点对应的时间长度，以及地震波由地面传播至所述多个目标点的平均速度，确定所述多个目标点距离所述地面的深度；
根据所述多个目标点距离所述地面的深度，显示所述目标点集所在的空间区域的画面。
7.一种砂体刻画装置，其特征在于，所述装置包括：
频率确定模块，被配置为根据钻井过程中检测到的砂体的厚度，确定所述砂体对应的敏感频率，所述敏感频率为对所述砂体检测能力最强的地震波的频率；
数据体提取模块，被配置为根据所述敏感频率所在的敏感频率范围，从三维地震数据体中提取多个目标频率对应的分频地震数据体，所述敏感频率范围包括所述多个目标频率；
振幅提取模块，被配置为从多个分频地震数据体中提取所述砂体所在位置处的振幅；
数据体选择模块，被配置为从所述多个分频地震数据体中选择目标地震数据体，所述目标地震数据体对应的振幅与所述砂体的厚度之间的相关度最大；
区域截取模块，被配置为根据所述砂体顶部距离地面的第一深度，以及所述砂体底部距离所述地面的第二深度，从所述目标地震数据体中截取所述砂体所在的目标区域；
砂体显示模块，被配置为显示所述目标区域的画面，所述目标区域包含所述砂体。
8.根据权利要求7所述的装置，其特征在于，
所述砂体显示模块，被配置为过滤掉所述目标区域中的第一区域，所述第一区域中的任意位置对应的振幅不大于参考振幅；显示过滤后的目标区域的画面，所述过滤后的目标区域包含所述砂体。
9.一种计算机设备，其特征在于，所述计算机设备包括处理器和存储器，所述存储器中存储有至少一条程序代码，所述程序代码由所述处理器加载并执行以实现如权利要求1至权利要求6任一项所述的砂体刻画方法所执行的操作。
10.一种计算机可读存储介质，其特征在于，所述存储介质中存储有至少一条程序代码，所述程序代码由处理器加载并执行以实现如权利要求1至权利要求6任一项所述的砂体刻画方法所执行的操作。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011276426.1A CN112698391B (zh) | 2020-11-16 | 2020-11-16 | 砂体刻画方法、装置、设备及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011276426.1A CN112698391B (zh) | 2020-11-16 | 2020-11-16 | 砂体刻画方法、装置、设备及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112698391A true CN112698391A (zh) | 2021-04-23 |
| CN112698391B CN112698391B (zh) | 2023-10-31 |
Family
ID=75505888
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011276426.1A Active CN112698391B (zh) | 2020-11-16 | 2020-11-16 | 砂体刻画方法、装置、设备及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112698391B (zh) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106526670A (zh) * | 2016-09-21 | 2017-03-22 | 中石化石油工程技术服务有限公司 | 一种碎屑岩储层中地震属性砂体空间分布描述及评价的方法 |
| CN109709612A (zh) * | 2018-11-07 | 2019-05-03 | 中国石油大学(北京) | 砂体分布预测方法、装置、设备及可读存储介质 |
-
2020
- 2020-11-16 CN CN202011276426.1A patent/CN112698391B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106526670A (zh) * | 2016-09-21 | 2017-03-22 | 中石化石油工程技术服务有限公司 | 一种碎屑岩储层中地震属性砂体空间分布描述及评价的方法 |
| CN109709612A (zh) * | 2018-11-07 | 2019-05-03 | 中国石油大学(北京) | 砂体分布预测方法、装置、设备及可读存储介质 |
Non-Patent Citations (1)
| Title |
|---|
| 安鹏 等: "敏感频率地震属性在薄层砂体预测中的应用――以松辽盆地肇源地区为例", 物探与化探, vol. 44, no. 02, pages 321 - 328 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112698391B (zh) | 2023-10-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109696706B (zh) | 盐构造检测方法、装置及存储介质 | |
| CN112922589B (zh) | 尖灭线确定方法、装置、终端及存储介质 | |
| CN110174699B (zh) | 确定断层封闭性的方法、装置及存储介质 | |
| CN110389380B (zh) | 地震剖面的同相轴自动拾取方法、装置及存储介质 | |
| CN111005715A (zh) | 确定气井产量的方法、装置及存储介质 | |
| CN111077570B (zh) | 确定油藏中辉绿岩的侵入位置的方法、装置和存储介质 | |
| CN111427082A (zh) | 地震响应模板的确定方法、装置及存储介质 | |
| CN110673214B (zh) | 预测水平井的入口靶点和末端点的深度的方法和装置 | |
| CN109281648B (zh) | 确定油藏的合理井网密度的方法和装置 | |
| CN109116424B (zh) | 地震波数据的低波数噪音分离方法、装置及存储介质 | |
| CN111736215B (zh) | 断层断距确定方法及装置 | |
| CN112698391B (zh) | 砂体刻画方法、装置、设备及存储介质 | |
| CN109343113B (zh) | 一种预测油气藏位置的方法和装置 | |
| CN110058306B (zh) | 确定三维速度体的方法、装置及计算机存储介质 | |
| CN112096373B (zh) | 选定页岩气区域的方法、装置及计算机存储介质 | |
| CN110630254B (zh) | 油水层的区分方法、装置及计算机可读存储介质 | |
| CN112987105B (zh) | 定量预测地下岩盐层分布的方法、装置、终端和存储介质 | |
| CN115263244B (zh) | 气水层的排水控制方法、装置及计算机存储介质 | |
| CN113586043B (zh) | 束缚水饱和参数的确定方法、装置和计算机设备 | |
| CN111812744B (zh) | 地层可压性确定方法、装置及计算机存储介质 | |
| CN115708001A (zh) | 预测油气储层的方法、装置、设备及存储介质 | |
| CN115576009A (zh) | 断层的确定方法、装置、计算机设备和存储介质 | |
| CN109116415B (zh) | 地震波数据分离方法、装置及存储介质 | |
| CN114527505A (zh) | 薄层识别方法、装置、设备及存储介质 | |
| CN115711123A (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 |