CN113077074A - 基于储层预测因子的储层预测方法及装置 - Google Patents
基于储层预测因子的储层预测方法及装置 Download PDFInfo
- Publication number
- CN113077074A CN113077074A CN202110243211.8A CN202110243211A CN113077074A CN 113077074 A CN113077074 A CN 113077074A CN 202110243211 A CN202110243211 A CN 202110243211A CN 113077074 A CN113077074 A CN 113077074A
- Authority
- CN
- China
- Prior art keywords
- reservoir
- reservoir prediction
- prediction factor
- curve
- factor
- 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
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000010219 correlation analysis Methods 0.000 claims abstract description 25
- 238000011160 research Methods 0.000 claims abstract description 25
- 238000009826 distribution Methods 0.000 claims description 29
- 239000004215 Carbon black (E152) Substances 0.000 claims description 28
- 229930195733 hydrocarbon Natural products 0.000 claims description 28
- 238000004364 calculation method Methods 0.000 claims description 18
- 125000001183 hydrocarbyl group Chemical group 0.000 claims description 17
- 238000004590 computer program Methods 0.000 claims description 16
- 238000004458 analytical method Methods 0.000 claims description 12
- 238000003860 storage Methods 0.000 claims description 12
- 150000002430 hydrocarbons Chemical class 0.000 claims description 11
- 238000010276 construction Methods 0.000 claims description 6
- 238000010586 diagram Methods 0.000 description 15
- 230000006870 function Effects 0.000 description 13
- 238000004891 communication Methods 0.000 description 12
- 238000012545 processing Methods 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 239000000872 buffer Substances 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 239000012530 fluid Substances 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 239000000284 extract Substances 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
- 239000003129 oil well Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Mining
Abstract
本发明提供一种基于储层预测因子的储层预测方法及装置，方法包括：获取测井曲线数据、测井解释成果和地震数据体；根据所述测井曲线数据构建储层预测因子；将所述储层预测因子与所述测井曲线数据中的GR测井曲线进行双参数相关分析，得到所述储层预测因子中的待定系数值；基于所述测井曲线数据以及所述待定系数值计算储层预测因子曲线；根据所述储层预测因子曲线和所述测井解释成果确定储层预测因子阈值；基于所述地震数据体和所述待定系数值计算储层预测因子地震数据体；根据所述储层预测因子阈值和所述储层预测因子地震数据体进行研究区含油气储层地震预测。本发明能够有效提高含油气储层预测的精度，降低储层预测的风险。
Description
技术领域
本发明涉及石油勘探技术领域，具体涉及一种基于储层预测因子的储层预测方法及装置。
背景技术
随着石油勘探开发程度的深入，油气勘探的目标不断复杂化、隐蔽化，勘探手段不断多样化、精细化。近些年随着储层预测技术的进步，新的理论与技术手段层出不穷，在油气勘探目标不断演化的同时地球物理勘探技术也得到了长足的发展，尤其是在基于地震数据的基于储层预测因子的储层预测方面。20世纪80年代AVO技术开始发展，利用振幅变化预测砂岩储层的流体分布，随着岩石物理技术以及电子计算机技术的进步，叠前AVO反演实现了直接从地震数据中获取纵波速度、横波速度及密度数据，为储层预测因子构建提供了更多可能选项，大大提高了储层预测的精度。叠前AVO反演利用地震反射振幅随偏移距变化的规律，结合测井信息从地震反射波数据体中提取出纵波速度、横波速度及密度信息，在保留AVO信息的同时将反射振幅信息转化到速度、密度等层属性信息，为后续构建弹性参数为基础的储层预测提供了基础。
近年来，基于岩石弹性信息构建储层预测因子已经成为储层预测领域研究的重点，常规的属性因子包括纵横波速度比、泊松阻抗、Gassmann流体因子等，这些属性因子最多只含有一个待定系数，有的不含待定系数，在实际应用过程中由于各地区地质背景具有差异性与复杂性导致固定的属性因子不能很好的适应工区条件，限制了储层预测因子的实用性，参数缺乏可变性还导致了储层预测的精度降低。在实际勘探过程中，发现纵波速度、横波速度以及密度三者对于含油气储层响应最为敏感，而常规的阻抗域属性因子由于固有的阶数限制没有充分发挥出速度以及密度的优势，导致地震储层预测阻抗域同阶属性因子敏感性受限问题。
发明内容
针对现有技术中的问题，本发明提供一种基于储层预测因子的储层预测方法及装置，能够有效提高含油气储层预测的精度，降低储层预测的风险。
为解决上述技术问题，本发明提供以下技术方案：
第一方面，本发明提供一种基于储层预测因子的储层预测方法，包括：
获取测井曲线数据、测井解释成果和地震数据体；
根据所述测井曲线数据构建储层预测因子；
将所述储层预测因子与所述测井曲线数据中的GR测井曲线进行双参数相关分析，得到所述储层预测因子中的待定系数值；
基于所述测井曲线数据以及所述待定系数值计算储层预测因子曲线；
根据所述储层预测因子曲线和所述测井解释成果确定储层预测因子阈值；
基于所述地震数据体和所述待定系数值计算储层预测因子地震数据体；
根据所述储层预测因子阈值和所述储层预测因子地震数据体进行研究区含油气储层地震预测。
其中，所述根据所述测井曲线数据构建储层预测因子，包括：
根据所述测井曲线数据中的纵波速度曲线、横波速度曲线和密度曲线构建储层预测因子。
其中，所述将所述储层预测因子与所述测井曲线数据中的GR测井曲线进行双参数相关分析，得到所述储层预测因子中的待定系数值，包括：
利用双参数扫描方法，对储层预测因子与所述GR测井曲线进行相关性分析，得到相关系数最大时的待定系数的取值。
其中，所述根据所述储层预测因子曲线和所述测井解释成果确定储层预测因子阈值，包括：
对所述储层预测因子曲线进行直方图分析，根据直方图分布确定含油气储层对应的储层预测因子阈值；
其中，直方图中横坐标为储层预测因子数值，纵坐标为百分比数值。
其中，所述根据所述储层预测因子阈值和所述储层预测因子地震数据体进行研究区含油气储层地震预测，包括：
根据所述储层预测因子地震数据体确定含油气储层分布范围；
确定所述含油气储层分布范围上的储层预测因子小于所述储层预测因子阈值的区域为含油气储层。
第二方面，本发明提供一种基于储层预测因子的储层预测装置，包括：
数据获取单元，用于获取测井曲线数据、测井解释成果和地震数据体；
储层预测因子构建单元，用于根据所述测井曲线数据构建储层预测因子；
待定系数计算单元，用于将所述储层预测因子与所述测井曲线数据中的GR测井曲线进行双参数相关分析，得到所述储层预测因子中的待定系数值；
储层预测因子测井曲线计算单元，用于基于所述测井曲线数据以及所述待定系数值计算储层预测因子曲线；
阈值确定单元，用于根据所述储层预测因子曲线和所述测井解释成果确定储层预测因子阈值；
储层预测因子地震数据计算单元，用于基于所述地震数据体和所述待定系数值计算储层预测因子地震数据体；
储层预测单元，用于根据所述储层预测因子阈值和所述储层预测因子地震数据体进行研究区含油气储层地震预测。
其中，所述储层预测因子构建单元包括：
储层预测因子构建子单元，用于根据所述测井曲线数据中的纵波速度曲线、横波速度曲线和密度曲线构建储层预测因子。
其中，所述待定系数计算单元包括：
待定系数计算子单元，用于利用双参数扫描方法，对储层预测因子与所述GR测井曲线进行相关性分析，得到相关系数最大时的待定系数的取值。
其中，所述阈值确定单元包括：
阈值确定子单元，用于对所述储层预测因子曲线进行直方图分析，根据直方图分布确定含油气储层对应的储层预测因子阈值；
其中，直方图中横坐标为储层预测因子数值，纵坐标为百分比数值。
其中，所述储层预测单元包括：
储层分布范围子单元，用于根据所述储层预测因子地震数据体确定含油气储层分布范围；
储层预测子单元，用于确定所述含油气储层分布范围上的储层预测因子小于所述储层预测因子阈值的区域为含油气储层。
第三方面，本发明提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现所述的基于储层预测因子的储层预测方法的步骤。
第四方面，本发明提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现所述的基于储层预测因子的储层预测方法的步骤。
由上述技术方案可知，本发明提供一种基于储层预测因子的储层预测方法及装置，通过获取测井曲线数据、测井解释成果和地震数据体；根据所述测井曲线数据构建储层预测因子；将所述储层预测因子与所述测井曲线数据中的GR测井曲线进行双参数相关分析，得到所述储层预测因子中的待定系数值；基于所述测井曲线数据以及所述待定系数值计算储层预测因子曲线；根据所述储层预测因子曲线和所述测井解释成果确定储层预测因子阈值；基于所述地震数据体和所述待定系数值计算储层预测因子地震数据体；根据所述储层预测因子阈值和所述储层预测因子地震数据体进行研究区含油气储层地震预测，能够有效提高含油气储层预测的精度，降低储层预测的风险。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例中的基于储层预测因子的储层预测方法的流程示意图；
图2为本发明实施例中储层预测因子岩石物理分析示意图；
图3为本发明实施例中储层预测因子曲线与GR测井曲线相关性分析图；
图4为本发明实施例中储层预测因子直方图分析图；
图5为本发明实施例中基于储层预测因子识别的含油气储层平面分布图；
图6为本发明实施例中的基于储层预测因子的储层预测装置的结构示意图。
图7为本发明实施例中的电子设备的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整的描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明提供一种基于储层预测因子的储层预测方法的实施例，参见图1，所述基于储层预测因子的储层预测方法具体包含有如下内容：
S101：获取测井曲线数据、测井解释成果和地震数据体；
在本步骤中，测井曲线数据包括：纵波速度曲线、横波速度曲线、密度曲线以及自然伽马测井曲线(GR测井曲线)。测井解释成果为测井岩性及油气水解释结论。地震数据体包括：基于常规叠前AVO反演得到的纵波速度、横波速度以及密度地震数据体。
S102：根据所述测井曲线数据构建储层预测因子；
在本步骤中，是根据测井曲线数据中的纵波速度曲线、横波速度曲线和密度曲线构建储层预测因子。其中，储层预测因子包含两项待定系数。
构建储层预测因子的公式如下：
其中，F为储层预测因子，Vp为纵波速度曲线，Vs为横波速度曲线，c1和c2为待定系数，ρ为密度曲线。参见图2所示储层预测因子岩石物理分析示意图，图2中横坐标为纵波阻抗(Vp·ρ)，纵坐标为纵横波速度比(Vp/Vs)，其中，随着纵横波速度比的提高，储层预测因子的值也变得更高。根据岩石物理理论含油气储层位于左下角椭圆区域，非储层位于其他区域。从图2中可以看出储层预测因子低值完全对应了含油气储层区域，高值则对应了非储层区域，表明储层预测因子能够有效的区分含油气储层与非储层。
S103：将所述储层预测因子与所述测井曲线数据中的GR测井曲线进行双参数相关分析，得到所述储层预测因子中的待定系数值；
在本步骤中，利用双参数扫描方法，对储层预测因子与所述GR测井曲线进行相关性分析，得到相关系数最大时的待定系数的取值。
在实施过程中，利用双参数扫描方法，对储层预测因子F与GR测井曲线进行相关性分析，得到相关系数最大时的待定系数(c1和c2)的取值，从而确定研究工区储层预测因子的计算公式，其中相关性分析公式如下所示：
其中，Cov(F,GR)为储层预测因子F与GR测井曲线的协方差，Var[F]为储层预测因子F的方差，Var[GR]为GR测井曲线的方差，R为相关系数，c1和c2为待定系数。当相关系数R达到最大值时，此时对应的c1和c2值即为该研究区对应的储层预测因子待定系数值。图3所示储层预测因子曲线与GR测井曲线相关性分析图，图3中横坐标为c1，纵坐标为c2。从图3中可以看出当c1＝1.8，c2＝2.25时(黑色圆圈位置)储层预测因子F与GR测井曲线的相关性最好，此时的相关系数R＝0.94，从而确定该研究工区内储层预测因子公式中的待定系数值，即c1＝1.8，c2＝2.25。
S104：基于所述测井曲线数据以及所述待定系数值计算储层预测因子曲线；
在实施过程中，将纵波速度曲线中的纵波速度、横波速度曲线中的横波速度、密度曲线中的密度和不收S103计算的待定系数值(c1＝1.8，c2＝2.25)代入到储层预测因子公式中，获得目标井对应的储层预测因子曲线，储层预测因子曲线计算公式为：
S105：根据所述储层预测因子曲线和所述测井解释成果确定储层预测因子阈值；
在本步骤中，对所述储层预测因子曲线进行直方图分析，根据直方图分布确定含油气储层对应的储层预测因子阈值；其中，直方图中横坐标为储层预测因子数值，纵坐标为百分比数值。
在实施过程中，对目标井的储层预测因子曲线进行直方图分析。根据直方图分布规律确定含油气储层对应的储层预测因子阈值，为后续含油气储层预测提供支持。图4展示了储层预测因子曲线直方图分析图，从图4中可以看出储层预测因子能够有效的区分含油气储层与非储层(黑色虚线位置)，且在储层预测因子域含油气储层与非储层的重叠部分占比非常小，当F<1.55时，黑色虚线左侧对应了含油气储层，当F≥1.55时，黑色虚线右侧基本对应了非储层区域，可以确定该研究工区储层预测因子的阈值为F＝1.55。
S106：基于所述地震数据体和所述待定系数值计算储层预测因子地震数据体；
在实施过程中，将纵波速度曲线的纵波速度、横波速度曲线的横波速度以及密度曲线的密度代入到步骤S104中的储层预测因子计算公式，得到研究工区对应的储层预测因子地震数据体，为后续全工区含油气储层分布规律预测提供数据支持。
S107：根据所述储层预测因子阈值和所述储层预测因子地震数据体进行研究区含油气储层地震预测。
在本步骤中，根据所述储层预测因子地震数据体确定含油气储层分布范围；确定所述含油气储层分布范围上的储层预测因子小于所述储层预测因子阈值的区域为含油气储层。
在实施过程中，根据直方图分析确定的储层预测因子阈值以及计算得到的储层预测因子地震数据体确定研究区含油气储层分布范围，即储层预测因子数值小于1.55的区域对应了含油气储层，大于1.55的区域为非储层。图5展示了研究区基于储层预测因子识别的含油气储层平面分布规律图，在研究区范围有3口已钻井，其中W-2和W-3井均为商业油井，日产油分别达到300桶/天和160桶/天，其测井解释结论为高孔高渗透性砂岩储层，而W-1井为失利井，测井解释结论为水井，后续的地质研究发现W-1井位于油水分界面以下。在储层预测因子识别的含油气储层平面分布规律图中，W-1井位于非储层分布范围内，W-2井和W-3井位于含油气储层分布范围内，表明基于储层预测因子得到的预测结果与钻井信息吻合，证实本方案的有效性。
从上述描述可知，本发明实施例提供的基于储层预测因子的储层预测方法，通过获取测井曲线数据、测井解释成果和地震数据体；根据所述测井曲线数据构建储层预测因子；将所述储层预测因子与所述测井曲线数据中的GR测井曲线进行双参数相关分析，得到所述储层预测因子中的待定系数值；基于所述测井曲线数据以及所述待定系数值计算储层预测因子曲线；根据所述储层预测因子曲线和所述测井解释成果确定储层预测因子阈值；基于所述地震数据体和所述待定系数值计算储层预测因子地震数据体；根据所述储层预测因子阈值和所述储层预测因子地震数据体进行研究区含油气储层地震预测，能够有效提高含油气储层预测的精度，降低储层预测的风险。
本发明实施例提供一种能够实现所述基于储层预测因子的储层预测方法中全部内容的基于储层预测因子的储层预测装置的具体实施方式，参见图6，所述基于储层预测因子的储层预测装置具体包括如下内容：
数据获取单元10，用于获取测井曲线数据、测井解释成果和地震数据体；
储层预测因子构建单元20，用于根据所述测井曲线数据构建储层预测因子；
待定系数计算单元30，用于将所述储层预测因子与所述测井曲线数据中的GR测井曲线进行双参数相关分析，得到所述储层预测因子中的待定系数值；
储层预测因子测井曲线计算单元40，用于基于所述测井曲线数据以及所述待定系数值计算储层预测因子曲线；
阈值确定单元50，用于根据所述储层预测因子曲线和所述测井解释成果确定储层预测因子阈值；
储层预测因子地震数据计算单元60，用于基于所述地震数据体和所述待定系数值计算储层预测因子地震数据体；
储层预测单元70，用于根据所述储层预测因子阈值和所述储层预测因子地震数据体进行研究区含油气储层地震预测。
其中，所述储层预测因子构建单元包括：
储层预测因子构建子单元，用于根据所述测井曲线数据中的纵波速度曲线、横波速度曲线和密度曲线构建储层预测因子。
其中，所述待定系数计算单元包括：
待定系数计算子单元，用于利用双参数扫描方法，对储层预测因子与所述GR测井曲线进行相关性分析，得到相关系数最大时的待定系数的取值。
其中，所述阈值确定单元包括：
阈值确定子单元，用于对所述储层预测因子曲线进行直方图分析，根据直方图分布确定含油气储层对应的储层预测因子阈值；
其中，直方图中横坐标为储层预测因子数值，纵坐标为百分比数值。
其中，所述储层预测单元包括：
储层分布范围子单元，用于根据所述储层预测因子地震数据体确定含油气储层分布范围；
储层预测子单元，用于确定所述含油气储层分布范围上的储层预测因子小于所述储层预测因子阈值的区域为含油气储层。
本发明提供的基于储层预测因子的储层预测装置的实施例具体可以用于执行上述实施例中的基于储层预测因子的储层预测方法的实施例的处理流程，其功能在此不再赘述，可以参照上述方法实施例的详细描述。
从上述描述可知，本发明实施例提供的基于储层预测因子的储层预测装置，通过获取测井曲线数据、测井解释成果和地震数据体；根据所述测井曲线数据构建储层预测因子；将所述储层预测因子与所述测井曲线数据中的GR测井曲线进行双参数相关分析，得到所述储层预测因子中的待定系数值；基于所述测井曲线数据以及所述待定系数值计算储层预测因子曲线；根据所述储层预测因子曲线和所述测井解释成果确定储层预测因子阈值；基于所述地震数据体和所述待定系数值计算储层预测因子地震数据体；根据所述储层预测因子阈值和所述储层预测因子地震数据体进行研究区含油气储层地震预测，能够有效提高含油气储层预测的精度，降低储层预测的风险。
本申请提供一种用于实现所述基于储层预测因子的储层预测方法中的全部或部分内容的电子设备的实施例所述电子设备具体包含有如下内容：
处理器(processor)、存储器(memory)、通信接口(Communications Interface)和总线；其中，所述处理器、存储器、通信接口通过所述总线完成相互间的通信；所述通信接口用于实现相关设备之间的信息传输；该电子设备可以是台式计算机、平板电脑及移动终端等，本实施例不限于此。在本实施例中，该电子设备可以参照实施例用于实现所述基于储层预测因子的储层预测方法的实施例及用于实现所述基于储层预测因子的储层预测装置的实施例进行实施，其内容被合并于此，重复之处不再赘述。
图7为本申请实施例的电子设备9600的系统构成的示意框图。如图7所示，该电子设备9600可以包括中央处理器9100和存储器9140；存储器9140耦合到中央处理器9100。值得注意的是，该图7是示例性的；还可以使用其他类型的结构，来补充或代替该结构，以实现电信功能或其他功能。
一实施例中，基于储层预测因子的储层预测功能可以被集成到中央处理器9100中。其中，中央处理器9100可以被配置为进行如下控制：
获取测井曲线数据、测井解释成果和地震数据体；根据所述测井曲线数据构建储层预测因子；将所述储层预测因子与所述测井曲线数据中的GR测井曲线进行双参数相关分析，得到所述储层预测因子中的待定系数值；基于所述测井曲线数据以及所述待定系数值计算储层预测因子曲线；根据所述储层预测因子曲线和所述测井解释成果确定储层预测因子阈值；基于所述地震数据体和所述待定系数值计算储层预测因子地震数据体；根据所述储层预测因子阈值和所述储层预测因子地震数据体进行研究区含油气储层地震预测。
从上述描述可知，本申请的实施例提供的电子设备，通过获取测井曲线数据、测井解释成果和地震数据体；根据所述测井曲线数据构建储层预测因子；将所述储层预测因子与所述测井曲线数据中的GR测井曲线进行双参数相关分析，得到所述储层预测因子中的待定系数值；基于所述测井曲线数据以及所述待定系数值计算储层预测因子曲线；根据所述储层预测因子曲线和所述测井解释成果确定储层预测因子阈值；基于所述地震数据体和所述待定系数值计算储层预测因子地震数据体；根据所述储层预测因子阈值和所述储层预测因子地震数据体进行研究区含油气储层地震预测，能够有效提高含油气储层预测的精度，降低储层预测的风险。
在另一个实施方式中，基于储层预测因子的储层预测装置可以与中央处理器9100分开配置，例如可以将基于储层预测因子的储层预测配置为与中央处理器9100连接的芯片，通过中央处理器的控制来实现基于储层预测因子的储层预测功能。
如图7所示，该电子设备9600还可以包括：通信模块9110、输入单元9120、音频处理器9130、显示器9160、电源9170。值得注意的是，电子设备9600也并不是必须要包括图7中所示的所有部件；此外，电子设备9600还可以包括图7中没有示出的部件，可以参考现有技术。
如图7所示，中央处理器9100有时也称为控制器或操作控件，可以包括微处理器或其他处理器装置和/或逻辑装置，该中央处理器9100接收输入并控制电子设备9600的各个部件的操作。
其中，存储器9140，例如可以是缓存器、闪存、硬驱、可移动介质、易失性存储器、非易失性存储器或其它合适装置中的一种或更多种。可储存上述与失败有关的信息，此外还可存储执行有关信息的程序。并且中央处理器9100可执行该存储器9140存储的该程序，以实现信息存储或处理等。
输入单元9120向中央处理器9100提供输入。该输入单元9120例如为按键或触摸输入装置。电源9170用于向电子设备9600提供电力。显示器9160用于进行图像和文字等显示对象的显示。该显示器例如可为LCD显示器，但并不限于此。
该存储器9140可以是固态存储器，例如，只读存储器(ROM)、随机存取存储器(RAM)、SIM卡等。还可以是这样的存储器，其即使在断电时也保存信息，可被选择性地擦除且设有更多数据，该存储器的示例有时被称为EPROM等。存储器9140还可以是某种其它类型的装置。存储器9140包括缓冲存储器9141(有时被称为缓冲器)。存储器9140可以包括应用/功能存储部9142，该应用/功能存储部9142用于存储应用程序和功能程序或用于通过中央处理器9100执行电子设备9600的操作的流程。
存储器9140还可以包括数据存储部9143，该数据存储部9143用于存储数据，例如联系人、数字数据、图片、声音和/或任何其他由电子设备使用的数据。存储器9140的驱动程序存储部9144可以包括电子设备的用于通信功能和/或用于执行电子设备的其他功能(如消息传送应用、通讯录应用等)的各种驱动程序。
通信模块9110即为经由天线9111发送和接收信号的发送机/接收机9110。通信模块(发送机/接收机)9110耦合到中央处理器9100，以提供输入信号和接收输出信号，这可以和常规移动通信终端的情况相同。
基于不同的通信技术，在同一电子设备中，可以设置有多个通信模块9110，如蜂窝网络模块、蓝牙模块和/或无线局域网模块等。通信模块(发送机/接收机)9110还经由音频处理器9130耦合到扬声器9131和麦克风9132，以经由扬声器9131提供音频输出，并接收来自麦克风9132的音频输入，从而实现通常的电信功能。音频处理器9130可以包括任何合适的缓冲器、解码器、放大器等。另外，音频处理器9130还耦合到中央处理器9100，从而使得可以通过麦克风9132能够在本机上录音，且使得可以通过扬声器9131来播放本机上存储的声音。
本发明的实施例还提供能够实现上述实施例中的基于储层预测因子的储层预测方法中全部步骤的一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述实施例中的基于储层预测因子的储层预测方法的全部步骤，例如，所述处理器执行所述计算机程序时实现下述步骤：
获取测井曲线数据、测井解释成果和地震数据体；根据所述测井曲线数据构建储层预测因子；将所述储层预测因子与所述测井曲线数据中的GR测井曲线进行双参数相关分析，得到所述储层预测因子中的待定系数值；基于所述测井曲线数据以及所述待定系数值计算储层预测因子曲线；根据所述储层预测因子曲线和所述测井解释成果确定储层预测因子阈值；基于所述地震数据体和所述待定系数值计算储层预测因子地震数据体；根据所述储层预测因子阈值和所述储层预测因子地震数据体进行研究区含油气储层地震预测。
从上述描述可知，本发明实施例提供的计算机可读存储介质，通过获取测井曲线数据、测井解释成果和地震数据体；根据所述测井曲线数据构建储层预测因子；将所述储层预测因子与所述测井曲线数据中的GR测井曲线进行双参数相关分析，得到所述储层预测因子中的待定系数值；基于所述测井曲线数据以及所述待定系数值计算储层预测因子曲线；根据所述储层预测因子曲线和所述测井解释成果确定储层预测因子阈值；基于所述地震数据体和所述待定系数值计算储层预测因子地震数据体；根据所述储层预测因子阈值和所述储层预测因子地震数据体进行研究区含油气储层地震预测，能够有效提高含油气储层预测的精度，降低储层预测的风险。
虽然本发明提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或客户端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。
本领域技术人员应明白，本说明书的实施例可提供为方法、装置(系统)或计算机程序产品。因此，本说明书实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、装置(系统)和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。本发明并不局限于任何单一的方面，也不局限于任何单一的实施例，也不局限于这些方面和/或实施例的任意组合和/或置换。而且，可以单独使用本发明的每个方面和/或实施例或者与一个或更多其他方面和/或其实施例结合使用。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围，其均应涵盖在本发明的权利要求和说明书的范围当中。
Claims (12)
1.一种基于储层预测因子的储层预测方法，其特征在于，包括：
获取测井曲线数据、测井解释成果和地震数据体；
根据所述测井曲线数据构建储层预测因子；
将所述储层预测因子与所述测井曲线数据中的GR测井曲线进行双参数相关分析，得到所述储层预测因子中的待定系数值；
基于所述测井曲线数据以及所述待定系数值计算储层预测因子曲线；
根据所述储层预测因子曲线和所述测井解释成果确定储层预测因子阈值；
基于所述地震数据体和所述待定系数值计算储层预测因子地震数据体；
根据所述储层预测因子阈值和所述储层预测因子地震数据体进行研究区含油气储层地震预测。
2.根据权利要求1所述的基于储层预测因子的储层预测方法，其特征在于，所述根据所述测井曲线数据构建储层预测因子，包括：
根据所述测井曲线数据中的纵波速度曲线、横波速度曲线和密度曲线构建储层预测因子。
3.根据权利要求1所述的基于储层预测因子的储层预测方法，其特征在于，所述将所述储层预测因子与所述测井曲线数据中的GR测井曲线进行双参数相关分析，得到所述储层预测因子中的待定系数值，包括：
利用双参数扫描方法，对储层预测因子与所述GR测井曲线进行相关性分析，得到相关系数最大时的待定系数的取值。
4.根据权利要求1所述的基于储层预测因子的储层预测方法，其特征在于，所述根据所述储层预测因子曲线和所述测井解释成果确定储层预测因子阈值，包括：
对所述储层预测因子曲线进行直方图分析，根据直方图分布确定含油气储层对应的储层预测因子阈值；
其中，直方图中横坐标为储层预测因子数值，纵坐标为百分比数值。
5.根据权利要求1所述的基于储层预测因子的储层预测方法，其特征在于，所述根据所述储层预测因子阈值和所述储层预测因子地震数据体进行研究区含油气储层地震预测，包括：
根据所述储层预测因子地震数据体确定含油气储层分布范围；
确定所述含油气储层分布范围上的储层预测因子小于所述储层预测因子阈值的区域为含油气储层。
6.一种基于储层预测因子的储层预测装置，其特征在于，包括：
数据获取单元，用于获取测井曲线数据、测井解释成果和地震数据体；
储层预测因子构建单元，用于根据所述测井曲线数据构建储层预测因子；
待定系数计算单元，用于将所述储层预测因子与所述测井曲线数据中的GR测井曲线进行双参数相关分析，得到所述储层预测因子中的待定系数值；
储层预测因子测井曲线计算单元，用于基于所述测井曲线数据以及所述待定系数值计算储层预测因子曲线；
阈值确定单元，用于根据所述储层预测因子曲线和所述测井解释成果确定储层预测因子阈值；
储层预测因子地震数据计算单元，用于基于所述地震数据体和所述待定系数值计算储层预测因子地震数据体；
储层预测单元，用于根据所述储层预测因子阈值和所述储层预测因子地震数据体进行研究区含油气储层地震预测。
7.根据权利要求6所述的基于储层预测因子的储层预测装置，其特征在于，所述储层预测因子构建单元包括：
储层预测因子构建子单元，用于根据所述测井曲线数据中的纵波速度曲线、横波速度曲线和密度曲线构建储层预测因子。
8.根据权利要求6所述的基于储层预测因子的储层预测装置，其特征在于，所述待定系数计算单元包括：
待定系数计算子单元，用于利用双参数扫描方法，对储层预测因子与所述GR测井曲线进行相关性分析，得到相关系数最大时的待定系数的取值。
9.根据权利要求6所述的基于储层预测因子的储层预测装置，其特征在于，所述阈值确定单元包括：
阈值确定子单元，用于对所述储层预测因子曲线进行直方图分析，根据直方图分布确定含油气储层对应的储层预测因子阈值；
其中，直方图中横坐标为储层预测因子数值，纵坐标为百分比数值。
10.根据权利要求6所述的基于储层预测因子的储层预测装置，其特征在于，所述储层预测单元包括：
储层分布范围子单元，用于根据所述储层预测因子地震数据体确定含油气储层分布范围；
储层预测子单元，用于确定所述含油气储层分布范围上的储层预测因子小于所述储层预测因子阈值的区域为含油气储层。
11.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现权利要求1至5任一项所述的基于储层预测因子的储层预测方法的步骤。
12.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现权利要求1至5任一项所述的基于储层预测因子的储层预测方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110243211.8A CN113077074B (zh) | 2021-03-05 | 2021-03-05 | 基于储层预测因子的储层预测方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110243211.8A CN113077074B (zh) | 2021-03-05 | 2021-03-05 | 基于储层预测因子的储层预测方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113077074A true CN113077074A (zh) | 2021-07-06 |
| CN113077074B CN113077074B (zh) | 2024-03-05 |
Family
ID=76610035
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110243211.8A Active CN113077074B (zh) | 2021-03-05 | 2021-03-05 | 基于储层预测因子的储层预测方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113077074B (zh) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101187961A (zh) * | 2007-11-20 | 2008-05-28 | 中国海洋石油总公司 | 射孔完井产能预测系统 |
| CN101630013A (zh) * | 2008-07-16 | 2010-01-20 | 中国石油天然气集团公司 | 一种叠前地震数据泊松比参数反演的方法 |
| CN102759748A (zh) * | 2012-07-04 | 2012-10-31 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于目标分析的交汇分析地震储层预测方法 |
| CN103091709A (zh) * | 2012-12-25 | 2013-05-08 | 中国石油天然气集团公司 | 获得纵波、转换波地震数据时间匹配关系的方法和装置 |
| CN111830562A (zh) * | 2019-04-16 | 2020-10-27 | 中国石油天然气股份有限公司 | 一种油气储层渗透率预测方法及装置 |
-
2021
- 2021-03-05 CN CN202110243211.8A patent/CN113077074B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101187961A (zh) * | 2007-11-20 | 2008-05-28 | 中国海洋石油总公司 | 射孔完井产能预测系统 |
| CN101630013A (zh) * | 2008-07-16 | 2010-01-20 | 中国石油天然气集团公司 | 一种叠前地震数据泊松比参数反演的方法 |
| CN102759748A (zh) * | 2012-07-04 | 2012-10-31 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于目标分析的交汇分析地震储层预测方法 |
| CN103091709A (zh) * | 2012-12-25 | 2013-05-08 | 中国石油天然气集团公司 | 获得纵波、转换波地震数据时间匹配关系的方法和装置 |
| CN111830562A (zh) * | 2019-04-16 | 2020-10-27 | 中国石油天然气股份有限公司 | 一种油气储层渗透率预测方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113077074B (zh) | 2024-03-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111025384B (zh) | 基于波形分类交汇融合的储层预测方法和装置 | |
| CN105425289A (zh) | 确定低频波阻抗的方法和装置 | |
| CN103869362A (zh) | 体曲率获取方法和设备 | |
| CN112147688B (zh) | 一种储层油气识别方法和装置 | |
| Clark | The prediction of hydraulic fracture azimuth through geological, core, and analytical studies | |
| Qobi et al. | Permeability determination from Stoneley waves in the Ara group carbonates, Oman | |
| CN113077074A (zh) | 基于储层预测因子的储层预测方法及装置 | |
| CN112684498A (zh) | 一种基于宽方位地震数据的储层裂缝预测方法及系统 | |
| CN112817042B (zh) | 沉积储层中油气识别方法及装置 | |
| CN111751879A (zh) | 储层孔隙度预测方法及装置 | |
| CN113933901B (zh) | 深度域地震数据拼接方法及装置 | |
| CN113075747B (zh) | 储层裂缝发育区域的预测方法及装置 | |
| CN112817059B (zh) | 峰峰组内幕储层确定方法及装置 | |
| CN106772596A (zh) | 一种确定叠前时间偏移速度场的方法及装置 | |
| Douglas et al. | Geostatistical model for the Arab-D Reservoir, North’Ain Dar Pilot, Ghawar Field, Saudi Arabia: An improved reservoir simulation model | |
| CN114859407A (zh) | 火山岩储层声学特征参数的确定方法及装置 | |
| CN112130209A (zh) | 岩溶储层预测方法及装置 | |
| CN113138412A (zh) | 深层页岩气孔隙度地震预测方法及装置 | |
| CN113945973B (zh) | 一种储层特征分析方法、存储介质及电子设备 | |
| CN112523740A (zh) | 用于超深海井下油气勘探的岩石类型识别方法及装置 | |
| CN113325476B (zh) | 地震资料静校正方法及装置 | |
| CN110632650B (zh) | 全波形反演数据的分时窗提取方法及装置 | |
| CN112859171B (zh) | 利用有效裂缝进行储层连通性分析的方法和装置 | |
| CN112444857B (zh) | 一种基于叠后地震资料的膏盐地球物理预测方法 | |
| US11867857B2 (en) | Method and system for updating a seismic velocity model |
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 |