CN112037306A - 砂体厚度图的绘制方法、装置、计算机设备及存储介质 - Google Patents
砂体厚度图的绘制方法、装置、计算机设备及存储介质 Download PDFInfo
- Publication number
- CN112037306A CN112037306A CN202010672420.XA CN202010672420A CN112037306A CN 112037306 A CN112037306 A CN 112037306A CN 202010672420 A CN202010672420 A CN 202010672420A CN 112037306 A CN112037306 A CN 112037306A
- Authority
- CN
- China
- Prior art keywords
- sand
- well point
- thickness
- depth domain
- sand body
- 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
- 239000004576 sand Substances 0.000 title claims abstract description 330
- 238000000034 method Methods 0.000 title claims abstract description 56
- 238000003860 storage Methods 0.000 title claims abstract description 16
- 230000035945 sensitivity Effects 0.000 claims abstract description 46
- 238000004590 computer program Methods 0.000 claims description 9
- 238000006243 chemical reaction Methods 0.000 claims description 7
- 238000013507 mapping Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000002349 favourable effect Effects 0.000 abstract 3
- 230000008021 deposition Effects 0.000 description 15
- 238000010586 diagram Methods 0.000 description 14
- 238000012937 correction Methods 0.000 description 8
- 238000009826 distribution Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 4
- 238000012795 verification Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 3
- 230000007480 spreading Effects 0.000 description 3
- 238000003892 spreading Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000012417 linear regression Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 208000035126 Facies Diseases 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000000875 corresponding effect Effects 0.000 description 1
- 239000010432 diamond Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/20—Drawing from basic elements, e.g. lines or circles
- G06T11/206—Drawing of charts or graphs
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Abstract
本发明实施例提供了一种砂体厚度图的绘制方法、装置、计算机设备及存储介质，其中，该方法包括：获取井点地震敏感属性或井点反演时间域砂体厚度，获取井点测井解释深度域砂体厚度；根据所述井点测井解释深度域砂体厚度，将所述井点地震敏感属性或所述井点反演时间域砂体厚度转化为井点反演深度域砂体厚度，得到反演深度域砂体厚度面；计算所述井点测井解释深度域砂体厚度与所述井点反演深度域砂体厚度的差值，基于所述差值获取差值趋势面；基于所述差值趋势面校正所述反演深度域砂体厚度面，得到砂体厚度图。该方案，有利于确保井点砂体厚度的一致性，符合沉积规律，有利于提高砂体厚度值的准确性，有利于使得无井区砂体展布更合理可靠。
Description
技术领域
本发明涉及石油地质技术领域，特别涉及一种砂体厚度图的绘制方法、装置、计算机设备及存储介质。
背景技术
砂体厚度图是油田储量的确定以及井位部署，优化的基础数据，如何准确地计算和绘制砂体厚度图是油田各项工作开展的关键。目前绘制砂体厚度图的方法有多种，大体可以分为两类：第一类为在井数较多的情况下，根据已钻井点的砂厚值，结合地质研究，在了解物源方向，沉积相展布、储层分布情况下，直接手工绘制或者使用软件插值成图；第二类是地震平面趋势约束(敏感属性或者反演结果)，提取能够反映砂体厚度的属性，通过井点校正成图。
第一类方法，由于只是基于宏观地质认识，依据井点砂体数据，用手工直接绘制或者软件插值成图，井间砂体分布缺乏资料信息支撑，人为因素影响较大，结果的合理性很难把握。
第二类基于地震平面趋势约束的井点校正成图方法，地震敏感属性(其中反演可定量反映砂体展布)可以较好地定性反映砂体展布，在一定程度上能够弥补第一类方法的不足，但由于井点数据是硬数据，直接用井点数据校正往往会导致井点零误差而较大的改变井间地震敏感属性的趋势。
目前第二类方法在实际应用中，主要有两种方式：一种是通过井点砂厚与属性(反演)时间域厚度建立关系，通过关系式，直接将属性(或反演)时间厚度转换为深度厚度；反演时间砂厚与井点砂体厚度具有一定的相关性，但相关系数基本只保持在0.5～0.8左右，时间域砂厚直接转换为深度厚度成图，虽然符合沉积规律，但在部分井点和井间的砂体厚度相差较大；第二种是通所有井点砂厚插值，属性(反演)时间域砂厚约束成图。所有井点插值时间厚度约束成图的方法，虽然井点厚度一致，但存在部分井点砂厚与反演时间域砂厚没有相关性的情况，所有井点插值，反演时间域砂厚约束成图的方法会出现局部与沉积规律不一致的砂体厚度预测情况。
经调研，以上方法均存在无法准确地绘制砂体厚度图的缺陷。
发明内容
本发明实施例提供了一种，以解决现有技术中砂体厚度图的绘制存在准确性低的技术问题。该方法包括：
获取井点地震敏感属性或井点反演时间域砂体厚度，获取井点测井解释深度域砂体厚度；
根据所述井点测井解释深度域砂体厚度，将所述井点地震敏感属性或所述井点反演时间域砂体厚度转化为井点反演深度域砂体厚度，得到反演深度域砂体厚度面；
计算所述井点测井解释深度域砂体厚度与所述井点反演深度域砂体厚度的差值，基于所述差值获取差值趋势面；
基于所述差值趋势面校正所述反演深度域砂体厚度面，得到砂体厚度图。
本发明实施例还提供了一种砂体厚度图的绘制装置，以解决现有技术中砂体厚度图的绘制存在准确性低的技术问题。该装置包括：
数据获取模块，用于获取井点地震敏感属性或井点反演时间域砂体厚度，获取井点测井解释深度域砂体厚度；
数据转换模块，用于将所述井点地震敏感属性或所述井点反演时间域砂体厚度转化为井点反演深度域砂体厚度，得到反演深度域砂体厚度面；
差值计算模块，用于计算所述井点测井解释深度域砂体厚度与所述井点反演深度域砂体厚度的差值，基于所述差值获取差值趋势面；
砂体厚度图绘制模块，用于基于所述差值趋势面校正所述反演深度域砂体厚度面，得到砂体厚度图。
本发明实施例还提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述任意的砂体厚度图的绘制方法，以解决现有技术中砂体厚度图的绘制存在准确性低的技术问题。
本发明实施例还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述任意的砂体厚度图的绘制方法的计算机程序，以解决现有技术中砂体厚度图的绘制存在准确性低的技术问题。
在本发明实施例中，提出了根据井点测井解释深度域砂体厚度，将井点地震敏感属性或井点反演时间域砂体厚度转化为井点反演深度域砂体厚度，并得到反演深度域砂体厚度面；进而计算井点测井解释深度域砂体厚度与所述井点反演深度域砂体厚度的差值，并基于所述差值获取差值趋势面；最后，基于差值趋势面校正所述反演深度域砂体厚度面，即可得到最终的砂体厚度图。由于将井点地震敏感属性或井点反演时间域砂体厚度转化为域砂体厚度符合沉积规律，井点反演深度域砂体厚度与井点测井解释深度域砂体厚度存井点反演深度在相关性，使得基于井点测井解释深度域砂体厚度与所述井点反演深度域砂体厚度的差值获取的差值趋势面也符合沉积规律，进而使用差值趋势面来校正所述反演深度域砂体厚度面以得到砂体厚度图，与现有技术相比，可以避免人工的影响因素，可以避免井间砂体分布缺乏资料信息支撑、井点和井间的砂体厚度相差较大等问题，有利于确保井点砂体厚度的一致性，符合沉积规律，有利于提高砂体厚度值的准确性，有利于使得无井区砂体展布更合理可靠，进而可以更好地为油田井位的部署提供指导数据。
附图说明
此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中：
图1是本发明实施例提供的一种砂体厚度图的绘制方法的流程图；
图2是本发明实施例提供的一种区块构造图及反演时间域砂体厚度图的对比示意图；
图3是本发明实施例提供的一种井点反演时间域砂体厚度与井点测井解释深度域砂体厚度的拟合相关关系示意图；
图4是本发明实施例提供的一种应用上述砂体厚度图的绘制方法得到的砂体厚度图的示意图；
图5是本发明实施例提供的一种实施上述砂体厚度图的绘制方法的示意流程图；
图6是本发明实施例提供的一种应用上述砂体厚度图的绘制方法得到的砂体厚度图的对比图；
图7是本发明实施例提供的一种砂体厚度误差图；
图8是本发明实施例提供的一种计算机设备的结构框图；
图9是本发明实施例提供的一种砂体厚度图的绘制装置的结构框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施方式和附图，对本发明做进一步详细说明。在此，本发明的示意性实施方式及其说明用于解释本发明，但并不作为对本发明的限定。
在本发明实施例中，提供了一种砂体厚度图的绘制方法，如图1所示，该方法包括：
步骤102：获取井点地震敏感属性或井点反演时间域砂体厚度，获取井点测井解释深度域砂体厚度；
步骤104：根据所述井点测井解释深度域砂体厚度，将所述井点地震敏感属性或所述井点反演时间域砂体厚度转化为井点反演深度域砂体厚度，得到反演深度域砂体厚度面；
步骤106：计算所述井点测井解释深度域砂体厚度与所述井点反演深度域砂体厚度的差值，基于所述差值获取差值趋势面；
步骤108：基于所述差值趋势面校正所述反演深度域砂体厚度面，得到砂体厚度图。
由图1所示的流程可知，在本发明实施例中，提出了根据井点测井解释深度域砂体厚度，将井点地震敏感属性或井点反演时间域砂体厚度转化为井点反演深度域砂体厚度，并得到反演深度域砂体厚度面；进而计算井点测井解释深度域砂体厚度与所述井点反演深度域砂体厚度的差值，并基于所述差值获取差值趋势面；最后，基于差值趋势面校正所述反演深度域砂体厚度面，即可得到最终的砂体厚度图。由于将井点地震敏感属性或井点反演时间域砂体厚度转化为域砂体厚度符合沉积规律，井点反演深度域砂体厚度与井点测井解释深度域砂体厚度存井点反演深度在相关性，使得基于井点测井解释深度域砂体厚度与所述井点反演深度域砂体厚度的差值获取的差值趋势面也符合沉积规律，进而使用差值趋势面来校正所述反演深度域砂体厚度面以得到砂体厚度图，与现有技术相比，可以避免人工的影响因素，可以避免井间砂体分布缺乏资料信息支撑、井点和井间的砂体厚度相差较大等问题，有利于确保井点砂体厚度的一致性，符合沉积规律，有利于提高砂体厚度值的准确性，有利于使得无井区砂体展布更合理可靠，进而可以更好地为油田井位的部署提供指导数据。
具体实施时，针对采用地震敏感属性或者地震反演进行了地震储层预测的工区，采用上述砂体厚度图的绘制方法可以准确地绘制砂体厚度的情况。
具体实施时，地震敏感属性或者反演时间域砂体厚度能够反映沉积规律、刻画砂体平面展布。
具体实施时，可以先采集区块的反演时间域砂体厚度及测井解释深度域砂体厚度数据，进而再提取出井点反演时间域砂体厚度及井点测井解释深度域砂体厚度。具体的，区块的反演时间域砂体厚度如图2所示，图2中左侧为区块的构造图，图2中右侧为区块的反演时间域砂体厚度图。
具体实施时，为了进一步提高砂体厚度图的准确性，在本实施例中，通过以下步骤使得井点测井解释深度域砂体厚度与井点反演深度域砂体厚度之间存在相关性，例如，根据所述井点测井解释深度域砂体厚度，将所述井点地震敏感属性或所述井点反演时间域砂体厚度转化为井点反演深度域砂体厚度，包括：
根据所述井点地震敏感属性和所述井点测井解释深度域砂体厚度建立第一相关关系，或者，根据所述井点反演时间域砂体厚度和所述井点测井解释深度域砂体厚度建立第二相关关系；
根据所述第一相关关系，将所述井点地震敏感属性转化为所述井点反演深度域砂体厚度；或者，根据所述第二相关关系，将所述井点反演时间域砂体厚度转化为所述井点反演深度域砂体厚度。
具体的，通过对井点地震敏感属性和井点测井解释深度域砂体厚度进行拟合，可以建立二者之间的第一相关关系，通过对井点反演时间域砂体厚度和井点测井解释深度域砂体厚度进行拟合，可以建立二者之间的第二相关关系，例如，如图3所示，以井点反演时间域砂体厚度和井点测井解释深度域砂体厚度进行拟合得到第二相关关系为例，图3的横坐标为井点反演时间域砂体厚度，图3的纵坐标为井点测井解释深度域砂体厚度，图3中的菱形表示井点砂体厚度，既包括井点反演时间域砂体厚度也包括井点测井解释深度域砂体厚度，直线表示二者拟合得到的第二相关关系。
具体的，如图3所示，井点反演时间域砂体厚度和井点测井解释深度域砂体厚度进行拟合得到的第二相关关系可以通过以下公式来表示：
y＝2.3756x-21.591
其中，y为井点测井解释深度域砂体厚度；x为井点反演时间域砂体厚度。
具体实施时，对井点地震敏感属性(或井点反演时间域砂体厚度)和井点测井解释深度域砂体厚度进行线性回归拟合建立相关关系的过程中，还可以计算相关系数，相关系数值越大说明相关性越好。例如，以井点反演时间域砂体厚度和井点测井解释深度域砂体厚度进行线性回归拟合建立相关关系为例，如图3所示，得到的相关系数R2为0.8072。
具体实施时，在建立相关关系时，可以将地震敏感属性(或者反演时间域砂体厚度)与井点测井解释深度域砂体厚度差异大的井挑选出来，进行测井解释核实，核实后参与相关关系的拟合。
具体实施时，得到相关关系后，即可将井点地震敏感属性或井点反演时间域砂体厚度转化为井点反演深度域砂体厚度，具体的，根据上述第一相关关系，即可将井点地震敏感属性转化为井点反演深度域砂体厚度，根据上述第二相关关系，即可将井点反演时间域砂体厚度转化为井点反演深度域砂体厚度；得到井点反演深度域砂体厚度后，进而根据移动平均法进行插值，可以得到反演深度域砂体厚度面。
具体实施时，得到井点反演深度域砂体厚度后，即可计算井点测井解释深度域砂体厚度与井点反演深度域砂体厚度的差值，进而基于差值采用移动平均法进行插值，即可得到差值趋势面。
具体实施时，得到差值趋势面后，即可采用差值趋势面来校正反演深度域砂体厚度面，进而即可得到最终的砂体厚度平面图。具体的，通过上述砂体厚度图的绘制方法得到的砂体厚度图如图4所示，由图4和图2可知，得到的砂体厚度图与采集的反演时间域砂体厚度图沉积规律一致。
具体实施时，以外蒙塔19区塔19-56区块为例，以下结合示例来描述上述砂体厚度图的绘制方法的实施过程，如图5所示，该过程包括以下步骤：
1)收集外蒙塔19区塔19-56区块地震敏感属性(或反演时间域砂体厚度)及测井解释深度域砂体厚度数据；
收集地震敏感属性(或反演时间域砂体厚度，如图2所示，外蒙塔19区塔19-56区块构造及反演时间域砂体厚度图，反演时间域砂体厚度图反映了研究区的沉积特征，研究区整体为东断西超、中央凹四周高的构造特征，受沉积和构造影响，砂体呈现在东部构造高点处薄，东部断层处砂体厚，并向西部变薄的特征)和33口井的测井解释深度域砂体厚度数据，进而提取井点的地震敏感属性(或井点反演时间域砂体厚度)，对井点地震敏感属性(或井点反演时间域砂体厚度)与井点测井解释深度域砂体厚度进行拟合，得到二者的相关关系。具体的，井点反演时间域砂体厚度与井点测井解释深度域砂体厚度进行拟合，得到二者的相关关系如图3所示。在建立相关关系时，可以将反演时间域砂体厚度与井点测井解释深度域砂体厚度差异大的井挑选出来，进行测井解释核实，核实后再参与相关关系式的拟合。
2)根据上一步骤求取的相关关系式，将井点地震敏感属性(或井点反演时间域砂体厚度)转换为井点反演深度域砂体厚度，进而得到反演深度域砂体厚度面；
3)提取井点反演深度域砂体厚度，求取井点测井解释深度域砂体厚度与井点反演深度域砂体厚度的差值；
4)基于差值根据移动平均法进行插值，求取井点测井解释深度域砂体厚度与井点反演深度域砂体厚度的差值趋势面；
5)通过差值趋势面对反演深度域砂体厚度面进行校正，求取最终的砂体厚度平面图，如图4所示，由图4和图2可知，得到的砂体厚度图与采集的反演时间域砂体厚度图沉积规律一致。
具体实施时，为了进一步验证上述砂体厚度图的绘制方法的准确性、可靠性，在本实施例中，在以西部两口井T19-81、T19-82参与校正与不参与校正的两种情况下实施上述砂体厚度图的绘制方法，西部两口井T19-81、T19-82参与校正的情况，即实施上述砂体厚度图的绘制方法，基于该两口井T19-81、T19-82的井点反演深度域砂体厚度和井点测井解释深度域砂体厚度计算差值以及差值趋势面，并基于该差值趋势面校正反演深度域砂体厚度面，得到砂体厚度图，如图6中的左图所示；西部两口井T19-81、T19-82不参与校正的情况，即实施上述砂体厚度图的绘制方法的过程中，不基于该两口井T19-81、T19-82的井点反演深度域砂体厚度和井点测井解释深度域砂体厚度计算差值，而是基于区块内其它井的反演深度域砂体厚度和井点测井解释深度域砂体厚度计算差值以及差值趋势面，区块内其它井不包括两口井T19-81、T19-82的砂体厚度的差值如下表1所示，形成的砂体厚度误差图如图7所示，进而采用求得的差值趋势面对得到的反演深度域砂体厚度面进行校正，得到砂体厚度图，如图6中的右图所示。通过对比图6中的左右两张图可知，两张图的砂体展布特征基本一致，图中黑线框中为T19-81、T19-82两口井对应的位置，可见上述砂体厚度图的绘制方法的实施并不会因为采用不同井点的数据而影响砂体厚度图的准确性，也可见上述砂体厚度图的绘制方法与现有技术中基于井点插值的方式不同，通过差值趋势面来校正反演深度域砂体厚度面，可以使得砂体厚度图符合沉积规律，井点、井间及无井区砂体展布更加合理准确；其中，不参与校正的两口井T19-81、T19-82处的砂体厚度误差百分比分别为1.88％、3.73％，从图6、图7及表1可以看出，采用上述砂体厚度图的绘制方法绘制的砂体厚度图对井点吻合度高，计算的无井区砂体厚度准确度高，绘制的砂体厚度图符合沉积规律，井点、井间及无井区砂体展布更加合理准确，可以更好地为井位部署与调整提供参考数据。
表1
在本实施例中，提供了一种计算机设备，如图8所示，包括存储器802、处理器804及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述任意的砂体厚度图的绘制方法。
具体的，该计算机设备可以是计算机终端、服务器或者类似的运算装置。
在本实施例中，提供了一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述任意的砂体厚度图的绘制方法的计算机程序。
具体的，计算机可读存储介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机可读存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读存储介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
基于同一发明构思，本发明实施例中还提供了一种砂体厚度图的绘制装置，如下面的实施例所述。由于砂体厚度图的绘制装置解决问题的原理与砂体厚度图的绘制方法相似，因此砂体厚度图的绘制装置的实施可以参见砂体厚度图的绘制方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。
图9是本发明实施例的砂体厚度图的绘制装置的一种结构框图，如图9所示，该装置包括：
数据获取模块902，用于获取井点地震敏感属性或井点反演时间域砂体厚度，获取井点测井解释深度域砂体厚度；
数据转换模块904，用于将所述井点地震敏感属性或所述井点反演时间域砂体厚度转化为井点反演深度域砂体厚度，得到反演深度域砂体厚度面；
差值计算模块906，用于计算所述井点测井解释深度域砂体厚度与所述井点反演深度域砂体厚度的差值，基于所述差值获取差值趋势面；
砂体厚度图绘制模块908，用于基于所述差值趋势面校正所述反演深度域砂体厚度面，得到砂体厚度图。
在一个实施例中，所述数据转换模块，包括：
相关关系建立单元，用于根据所述井点地震敏感属性和所述井点测井解释深度域砂体厚度建立第一相关关系，或者，根据所述井点反演时间域砂体厚度和所述井点测井解释深度域砂体厚度建立第二相关关系；
数据转化单元，用于根据所述第一相关关系，将所述井点地震敏感属性转化为所述井点反演深度域砂体厚度；或者，根据所述第二相关关系，将所述井点反演时间域砂体厚度转化为所述井点反演深度域砂体厚度。
在一个实施例中，所述第二相关关系通过以下公式表示：
y＝2.3756x-21.591
其中，y为井点测井解释深度域砂体厚度；x为井点反演时间域砂体厚度。
在一个实施例中，所述差值计算模块，具体用于基于所述差值采用移动平均法进行插值，得到所述差值趋势面。
本发明实施例实现了如下技术效果：提出了根据井点测井解释深度域砂体厚度，将井点地震敏感属性或井点反演时间域砂体厚度转化为井点反演深度域砂体厚度，并得到反演深度域砂体厚度面；进而计算井点测井解释深度域砂体厚度与所述井点反演深度域砂体厚度的差值，并基于所述差值获取差值趋势面；最后，基于差值趋势面校正所述反演深度域砂体厚度面，即可得到最终的砂体厚度图。由于将井点地震敏感属性或井点反演时间域砂体厚度转化为域砂体厚度符合沉积规律，井点反演深度域砂体厚度与井点测井解释深度域砂体厚度存井点反演深度在相关性，使得基于井点测井解释深度域砂体厚度与所述井点反演深度域砂体厚度的差值获取的差值趋势面也符合沉积规律，进而使用差值趋势面来校正所述反演深度域砂体厚度面以得到砂体厚度图，与现有技术相比，可以避免人工的影响因素，可以避免井间砂体分布缺乏资料信息支撑、井点和井间的砂体厚度相差较大等问题，有利于确保井点砂体厚度的一致性，符合沉积规律，有利于提高砂体厚度值的准确性，有利于使得无井区砂体展布更合理可靠，进而可以更好地为油田井位的部署提供指导数据。
显然，本领域的技术人员应该明白，上述的本发明实施例的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行，并且在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样，本发明实施例不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种砂体厚度图的绘制方法，其特征在于，包括：
获取井点地震敏感属性或井点反演时间域砂体厚度，获取井点测井解释深度域砂体厚度；
根据所述井点测井解释深度域砂体厚度，将所述井点地震敏感属性或所述井点反演时间域砂体厚度转化为井点反演深度域砂体厚度，得到反演深度域砂体厚度面；
计算所述井点测井解释深度域砂体厚度与所述井点反演深度域砂体厚度的差值，基于所述差值获取差值趋势面；
基于所述差值趋势面校正所述反演深度域砂体厚度面，得到砂体厚度图。
2.如权利要求1所述的砂体厚度图的绘制方法，其特征在于，根据所述井点测井解释深度域砂体厚度，将井点地震敏感属性或井点反演时间域砂体厚度转化为井点反演深度域砂体厚度，包括：
根据所述井点地震敏感属性和所述井点测井解释深度域砂体厚度建立第一相关关系，或者，根据所述井点反演时间域砂体厚度和所述井点测井解释深度域砂体厚度建立第二相关关系；
根据所述第一相关关系，将所述井点地震敏感属性转化为所述井点反演深度域砂体厚度；或者，根据所述第二相关关系，将所述井点反演时间域砂体厚度转化为所述井点反演深度域砂体厚度。
3.如权利要求2所述的砂体厚度图的绘制方法，其特征在于，所述第二相关关系通过以下公式表示：
y＝2.3756x-21.591
其中，y为井点测井解释深度域砂体厚度；x为井点反演时间域砂体厚度。
4.如权利要求1至3中任一项所述的砂体厚度图的绘制方法，其特征在于，基于所述差值获取差值趋势面，包括：
基于所述差值采用移动平均法进行插值，得到所述差值趋势面。
5.一种砂体厚度图的绘制装置，其特征在于，包括：
数据获取模块，用于获取井点地震敏感属性或井点反演时间域砂体厚度，获取井点测井解释深度域砂体厚度；
数据转换模块，用于将所述井点地震敏感属性或所述井点反演时间域砂体厚度转化为井点反演深度域砂体厚度，得到反演深度域砂体厚度面；
差值计算模块，用于计算所述井点测井解释深度域砂体厚度与所述井点反演深度域砂体厚度的差值，基于所述差值获取差值趋势面；
砂体厚度图绘制模块，用于基于所述差值趋势面校正所述反演深度域砂体厚度面，得到砂体厚度图。
6.如权利要求5所述的砂体厚度图的绘制装置，其特征在于，所述数据转换模块，包括：
相关关系建立单元，用于根据所述井点地震敏感属性和所述井点测井解释深度域砂体厚度建立第一相关关系，或者，根据所述井点反演时间域砂体厚度和所述井点测井解释深度域砂体厚度建立第二相关关系；
数据转化单元，用于根据所述第一相关关系，将所述井点地震敏感属性转化为所述井点反演深度域砂体厚度；或者，根据所述第二相关关系，将所述井点反演时间域砂体厚度转化为所述井点反演深度域砂体厚度。
7.如权利要求6所述的砂体厚度图的绘制装置，其特征在于，所述第二相关关系通过以下公式表示：
y＝2.3756x-21.591
其中，y为井点测井解释深度域砂体厚度；x为井点反演时间域砂体厚度。
8.如权利要求5至7中任一项所述的砂体厚度图的绘制装置，其特征在于，所述差值计算模块，具体用于基于所述差值采用移动平均法进行插值，得到所述差值趋势面。
9.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至4中任一项所述的砂体厚度图的绘制方法。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至4中任一项所述的砂体厚度图的绘制方法的计算机程序。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010672420.XA CN112037306A (zh) | 2020-07-14 | 2020-07-14 | 砂体厚度图的绘制方法、装置、计算机设备及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010672420.XA CN112037306A (zh) | 2020-07-14 | 2020-07-14 | 砂体厚度图的绘制方法、装置、计算机设备及存储介质 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112037306A true CN112037306A (zh) | 2020-12-04 |
Family
ID=73579549
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010672420.XA Pending CN112037306A (zh) | 2020-07-14 | 2020-07-14 | 砂体厚度图的绘制方法、装置、计算机设备及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112037306A (zh) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103454685A (zh) * | 2013-08-09 | 2013-12-18 | 中国石油天然气股份有限公司 | 利用测井约束波阻抗反演预测砂体厚度的方法和装置 |
| CN107329171A (zh) * | 2017-06-07 | 2017-11-07 | 中国石油天然气股份有限公司 | 深度域储层地震反演方法及装置 |
| CN108051860A (zh) * | 2017-11-28 | 2018-05-18 | 北京金双狐油气技术有限公司 | 一种利用趋势面形态生成误差校正面的方法 |
-
2020
- 2020-07-14 CN CN202010672420.XA patent/CN112037306A/zh active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103454685A (zh) * | 2013-08-09 | 2013-12-18 | 中国石油天然气股份有限公司 | 利用测井约束波阻抗反演预测砂体厚度的方法和装置 |
| CN107329171A (zh) * | 2017-06-07 | 2017-11-07 | 中国石油天然气股份有限公司 | 深度域储层地震反演方法及装置 |
| CN108051860A (zh) * | 2017-11-28 | 2018-05-18 | 北京金双狐油气技术有限公司 | 一种利用趋势面形态生成误差校正面的方法 |
Non-Patent Citations (9)
| Title |
|---|
| 于建国, 路慎强, 王金铎, 于正军, 谷玉田, 崔世凌: "测井约束反演技术在史南浊积砂体描述中的应用", 石油物探, vol. 40, no. 02, 30 July 2001 (2001-07-30), pages 102 - 107 * |
| 刘宝国;: "适用于水平井轨迹设计的深度域随机模拟反演技术", 石油物探, vol. 53, no. 01, 25 January 2014 (2014-01-25), pages 93 - 98 * |
| 张兆辉;李胜军;张忠义;蔡刚;苏明军;: "目标砂体厚度精细预测技术在歧南地区的应用", 地球物理学进展, vol. 29, no. 05, 15 October 2014 (2014-10-15), pages 2315 - 2318 * |
| 李姊桐;张庆国;邢小林;: "基于测井系列校正的趋势面分析法标准化研究", 地质找矿论丛, vol. 33, no. 02, 15 June 2018 (2018-06-15), pages 272 - 277 * |
| 李宗田, 刘伟: "非线性约束储层反演技术在卫城油田的应用", 西安石油大学学报(自然科学版), vol. 20, no. 05, 25 September 2005 (2005-09-25), pages 20 * |
| 洪有密 编: "《石油大学测井技术论文集 第二辑》", vol. 1, 30 April 1994, 石油大学出版社, pages: 205 - 207 * |
| 熊冉;高亮;杨姣;赵继龙;: "曲线重构反演在储层横向预测中的应用", 西南石油大学学报(自然科学版), vol. 34, no. 1, 8 November 2011 (2011-11-08), pages 83 - 89 * |
| 边立恩;于茜;张京思;王国强;刘腾;: "容积法储量计算中砂体有效厚度高精度成图方法", 大庆石油地质与开发, vol. 35, no. 03, 1 June 2016 (2016-06-01), pages 43 - 48 * |
| 邢兴;赵勇刚;赵永峰;: "地震反演地震地质储层参数转换应用――以歧口地区为例", 重庆科技学院学报(自然科学版), vol. 19, no. 03, 15 June 2017 (2017-06-15), pages 66 - 70 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106405651B (zh) | 一种基于测井匹配的全波形反演初始速度模型构建方法 | |
| CN110955982B (zh) | 一种变质岩储层渗透率计算方法、装置及计算机存储介质 | |
| FR2842321A1 (fr) | Methode pour contraindre un champ de permeabilite heterogene representant un reservoir souterrain par des donnees dynamiques | |
| NO20101734A1 (no) | Hastighetsmodeller for en enkeltbronner og for et sett av bronner | |
| CN111008482B (zh) | 一种变质岩储层含水饱和度计算方法及装置 | |
| CN112037306A (zh) | 砂体厚度图的绘制方法、装置、计算机设备及存储介质 | |
| CN110309166B (zh) | 一种可溯源的地理高程数据补齐方法 | |
| CN104832161A (zh) | 基于双尺度相关对比的自动深度校正方法及系统 | |
| US11927713B2 (en) | Method, device and medium for acquiring logging parameters | |
| CN112666615B (zh) | 基于Tilt法获取精确航空磁测数据的方法、装置 | |
| CN112630840B (zh) | 基于统计特征参数的随机反演方法及处理器 | |
| CN115113268A (zh) | 一种基于米氏旋回的三级层序划分方法 | |
| CN112363243A (zh) | 预测页岩油优质储层的方法、装置、计算机设备及介质 | |
| CN113496070B (zh) | 地层俘获截面曲线的处理方法、装置、设备及介质 | |
| CN113495296B (zh) | 层析静校正量的确定方法、装置、设备及可读存储介质 | |
| CN111812716A (zh) | 页岩气储层总有机碳含量叠前定量预测方法、装置及设备 | |
| CN112415580B (zh) | 消除速度模型突变界面的方法、叠前深度偏移处理方法 | |
| CN113219534B (zh) | 一种叠前深度偏移速度质控方法、装置、介质及电子设备 | |
| CN112269212A (zh) | 测井小分层的地震解释层位确定方法、装置、设备及介质 | |
| CN114063165B (zh) | 一种三维地震数据的拼接方法及装置 | |
| CN111898091B (zh) | 油水界面的确定方法及装置 | |
| CN109581498B (zh) | 储层厚度分布确定方法、系统、设备、可读介质 | |
| CN110824565B (zh) | 一种空变子波体的构建方法及系统 | |
| CN111856585A (zh) | 时深转换速度场的构建方法、装置、设备及可读存储介质 | |
| CN113820745A (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 |