CN114077784A - 覆压孔隙度的确定方法、装置、存储介质和计算机设备 - Google Patents

覆压孔隙度的确定方法、装置、存储介质和计算机设备 Download PDF

Info

Publication number
CN114077784A
CN114077784A CN202010848459.2A CN202010848459A CN114077784A CN 114077784 A CN114077784 A CN 114077784A CN 202010848459 A CN202010848459 A CN 202010848459A CN 114077784 A CN114077784 A CN 114077784A
Authority
CN
China
Prior art keywords
porosity
overburden
pressure
model
target reservoir
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
Application number
CN202010848459.2A
Other languages
English (en)
Inventor
杨振恒
鲍云杰
韩志艳
申宝剑
卢龙飞
张培先
党伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Sinopec Exploration and Production Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Exploration and Production Research Institute
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Sinopec Exploration and Production Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN202010848459.2A priority Critical patent/CN114077784A/zh
Publication of CN114077784A publication Critical patent/CN114077784A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供一种覆压孔隙度的确定方法、装置、存储介质和计算机设备。其中,覆压孔隙度的确定方法,包括以下步骤:构建目标储层的覆压孔隙度模型,所述目标储层的覆压孔隙度模型为目标储层的覆压孔隙度与上覆压力之间的对数关系模型;基于目标储层的覆压孔隙度模型,根据目标储层的上覆压力和常压孔隙度,确定目标储层的覆压孔隙度。本发明提供一种覆压孔隙度的确定方法,能够快速准确地计算不同压力下目标储层的覆压孔隙度,为页岩气藏评价提供数据基础,对于评价页岩储层及页岩油气储量有重要的意义。并可以用于我国南方页岩气勘探开发评价中确定覆压条件下的孔隙度。

Description

覆压孔隙度的确定方法、装置、存储介质和计算机设备
技术领域
本发明涉及石油、地质、矿业勘探和开发技术领域,尤其涉及一种覆压孔隙度的确定方法、装置、存储介质和计算机设备。
背景技术
在页岩气勘探开发中,页岩的孔隙度是研究评价目标储层的一个重要参数。通常情况下,反映储层物性特征的岩心孔隙度是在常压下通过实验分析取得的。实际上,页岩深埋地下,地层压力等影响因素导致孔隙度与地面常规实验条件下通过岩心分析所确定的孔隙度的值不一致。因此,需通过覆压校正才能准确的还原地下页岩储层孔隙度的真实特征。
通过文献调研发现,现有的储层物性压力敏感性研究均认为地下储层孔隙度与上覆压力有着密切的关系,不同学者针对致密覆压下岩石的孔隙度进行了研究。有的研究发现了孔隙度随上覆岩层压力变大而变小的趋势,但缺乏更深入的研究,没有建立相对应的数学模型。有的研究建立了关于目标储层的孔隙度的数学模型,但是用该模型计算目标储层的准确度较差,无法实际使用。
因此,亟需一种更加准确的确定目标储层的孔隙度的方法。
发明内容
本发明的主要目的是提供一种覆压孔隙度的确定方法、装置、存储介质和计算机设备,以准确确定目标储层的覆压孔隙度。
第一方面,本申请提供一种覆压孔隙度的确定方法,包括以下步骤:构建目标储层的覆压孔隙度模型,所述目标储层的覆压孔隙度模型为目标储层的覆压孔隙度与上覆压力之间的对数关系模型;基于目标储层的覆压孔隙度模型,根据目标储层的上覆压力和常压孔隙度,确定目标储层的覆压孔隙度。
在一个实施例中,通过以下步骤构建覆压孔隙度模型:获取若干份目标储层的岩心样品,对于每份岩心样品,依次施加常压以及高于常压的多个不同上覆压力,测量其在常压下的常压孔隙度以及在高于常压的各个不同上覆压力下的覆压孔隙度;建立目标储层的覆压孔隙度初始模型,所述覆压孔隙度初始模型包括待确定的压力项系数和模型常数;对于每份岩心样品,基于覆压孔隙度初始模型,根据其在常压下的常压孔隙度以及在高于常压的各个不同上覆压力下的覆压孔隙度,利用曲线拟合的方法,确定与该份岩心样品对应的覆压孔隙度模型,其中,与该份岩心样品对应的覆压孔隙度模型包括与该份岩心样品对应的压力项系数和模型常数;根据与每份岩心样品对应的压力项系数确定覆压孔隙度初始模型中待确定的压力项系数,根据每份岩心样品的常压孔隙度和与每份岩心样品对应的模型常数确定覆压孔隙度初始模型中的模型常数;将确定了压力项系数和模型常数的覆压孔隙度初始模型作为目标储层的覆压孔隙度模型。
在一个实施例中,根据与每份岩心样品对应的压力项系数确定覆压孔隙度初始模型中待确定的压力项系数,包括:覆压孔隙度初始模型中待确定的压力项系数等于每份岩心样品对应的压力项系数的均值。
在一个实施例中,根据每份岩心样品的常压孔隙度和与每份岩心样品对应的模型常数确定覆压孔隙度初始模型中的模型常数,包括:根据每份岩心样品的常压孔隙度和与每份岩心样品对应的模型常数,利用线性拟合的方法,建立覆压孔隙度初始模型中的模型常数与常压孔隙度之间的对应关系;基于所述对应关系,根据目标储层的常压孔隙度确定目标储层的覆压孔隙度初始模型中的模型常数。
在一个实施例中,所述覆压孔隙度初始模型为:
ΦP=-a ln(P)+b
其中,ΦP表示目标储层在上覆压力P下的覆压孔隙度,P表示目标储层的上覆压力,a表示待确定的压力项系数,b表示模型常数。
在一个实施例中,覆压孔隙度初始模型中的模型常数与常压孔隙度之间的关系为以下线性关系:
b=cΦ0+d
其中,b表示模型常数,c表示斜率,Φ0表示目标储层的常压孔隙度,d表示常数。
在一个实施例中,所获取的目标储层的岩心样品的份数至少为12份,对于每份岩心样品,所施加的高于常压的不同上覆压力至少为6个。
第二方面,本申请提供一种覆压孔隙度的确定装置,包括:模型构建模块,用于构建目标储层的覆压孔隙度模型,所述目标储层的覆压孔隙度模型为目标储层的覆压孔隙度与上覆压力之间的对数关系模型;孔隙度确定模块,用于基于目标储层的覆压孔隙度模型,根据目标储层的上覆压力和常压孔隙度,确定目标储层的覆压孔隙度。
第三方面,本申请提供一种存储介质,存储有计算机程序,其特征在于,所述计算机程序被处理器执行时,实现如上文所述的覆压孔隙度的确定方法的步骤。
第四方面,本申请提供一种计算机设备,包括处理器和存储有程序代码的存储介质,所述程序代码被所述处理器执行时,实现如上文所述的覆压孔隙度的确定方法的步骤。
本发明提供一种覆压孔隙度的确定方法,能够快速准确地计算不同压力下目标储层的覆压孔隙度,为页岩气藏评价提供数据基础,对于评价页岩储层及页岩油气储量有重要的意义。并可以用于我国南方页岩气勘探开发评价中确定覆压条件下的孔隙度。
附图说明
构成本申请的一部分的说明书附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定,在附图中:
图1为根据本申请一示例性实施方式的覆压孔隙度的确定方法的流程图;
图2为根据本申请一示例性实施方式的目标储层的覆压孔隙度模型的构建方法的流程图;
图3为根据本申请一具体实施例的覆压孔隙度的确定方法的流程图;
图4为根据本申请一具体实施例的岩心样品在2MPa上覆压力下的实测孔隙度与通过覆压孔隙度模型计算所得覆压孔隙度的对比图;
图5为根据本申请一具体实施例的岩心样品在41.6MPa上覆压力下的实测孔隙度与通过覆压孔隙度模型计算所得覆压孔隙度的对比图。
具体实施方式
需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
实施例一
图1为根据本申请一示例性实施方式的覆压孔隙度的确定方法的流程图。如图1所示,本实施例提供一种覆压孔隙度的确定方法,可以包括以下步骤:
S100:构建目标储层的覆压孔隙度模型,所述目标储层的覆压孔隙度模型为目标储层的覆压孔隙度与上覆压力之间的对数关系模型。
S200:基于目标储层的覆压孔隙度模型,根据目标储层的上覆压力和常压孔隙度,确定目标储层的覆压孔隙度。
本发明提供的覆压孔隙度的确定方法,能够快速准确地计算不同压力下目标储层的覆压孔隙度,为页岩气藏评价提供数据基础,对于评价页岩储层及页岩油气储量有重要的意义。
实施例二
图2为根据本申请一示例性实施方式的目标储层的覆压孔隙度模型的构建方法的流程图。如图2所示,可以通过以下步骤构建覆压孔隙度模型:
S110:获取若干份目标储层的岩心样品,对于每份岩心样品,依次施加常压以及高于常压的多个不同上覆压力,测量其在常压下的常压孔隙度以及在高于常压的各个不同上覆压力下的覆压孔隙度。
S120:建立目标储层的覆压孔隙度初始模型,所述覆压孔隙度初始模型包括待确定的压力项系数和模型常数。
目标储层的覆压孔隙度初始模型可以为:
ΦP=-a ln(P)+b
其中,ΦP表示目标储层在上覆压力P下的覆压孔隙度,P表示目标储层的上覆压力,a表示待确定的压力项系数,b表示模型常数。
S130:对于每份岩心样品,基于覆压孔隙度初始模型,根据其在常压下的常压孔隙度以及在高于常压的各个不同上覆压力下的覆压孔隙度,利用曲线拟合的方法,确定与该份岩心样品对应的覆压孔隙度模型,其中,与该份岩心样品对应的覆压孔隙度模型包括与该份岩心样品对应的压力项系数和模型常数。
其中,所获取的目标储层的岩心样品的份数至少为12份,对于每份岩心样品,所施加的高于常压的不同上覆压力至少为6个。
作为一个示例,在施加高于常压的多个不同上覆压力时,在较低的上覆压力范围可以对岩心样品施加相互之间差值较小的多个不同的上覆压力,在较高的上覆压力范围可以对岩心样品施加相互之间差值较大的多个不同的上覆压力范围。
S140:根据与每份岩心样品对应的压力项系数确定覆压孔隙度初始模型中待确定的压力项系数,根据每份岩心样品的常压孔隙度和与每份岩心样品对应的模型常数确定覆压孔隙度初始模型中的模型常数。
其中,根据与每份岩心样品对应的压力项系数确定覆压孔隙度初始模型中待确定的压力项系数,可以包括:覆压孔隙度初始模型中待确定的压力项系数等于每份岩心样品对应的压力项系数的均值。
其中,根据每份岩心样品的常压孔隙度和与每份岩心样品对应的模型常数确定覆压孔隙度初始模型中的模型常数,可以包括:根据每份岩心样品的常压孔隙度和与每份岩心样品对应的模型常数,利用线性拟合的方法,建立覆压孔隙度初始模型中的模型常数与常压孔隙度之间的对应关系;基于所述对应关系,根据目标储层的常压孔隙度确定目标储层的覆压孔隙度初始模型中的模型常数。
覆压孔隙度初始模型中的模型常数与常压孔隙度之间的关系可以为以下线性关系:
b=cΦ0+d
其中,b表示模型常数,c表示斜率,Φ0表示目标储层的常压孔隙度,d表示常数。
S150:将确定了压力项系数和模型常数的覆压孔隙度初始模型作为目标储层的覆压孔隙度模型。
本实施例通过研究目标储层的岩心物性与上覆压力之间的关系,构建了的目标储层的覆压孔隙度模型,建立了目标储层的覆压孔隙度与上覆压力之间的对数关系,相比于现有技术,本实施例的目标储层的覆压孔隙度模型能够更加准确地通过目标储层的上覆压力计算目标储层的覆压孔隙度。
实施例三
图3为根据本申请一具体实施例的覆压孔隙度的确定方法的流程图。
本实施例提供的覆压孔隙度的确定方法可以包括以下步骤:
第一步,选取南方五峰组-龙马溪组页岩,取若干份芯井样品,对岩心样品开展常压实验,并获取实验结果,其中,该实验结果包括取芯样品的常压孔隙度。在本实施例中,以取12份岩心样品为例进行说明。常压压力为0.10133MPa。
第二步,对经过了常压实验的取芯样品开展覆压实验,并获取实验结果,其中,该实验结果包括取芯样品的覆压孔隙度。在本实施例中,利用6个不同的上覆压力对每份取芯样品进行覆压实验,并获取每份取芯样品在各个不同的上覆压力下的覆压孔隙度。
其中,所获取的实验结果如表1所示。在表1中,TOC(Total Organic Carbon)表示总有机碳,Rb表示沥青反射率。
表1
Figure BDA0002643896050000061
第三步,建立目标储层的覆压孔隙度初始模型,所述覆压孔隙度初始模型包括待确定的压力项系数和模型常数。
目标储层的覆压孔隙度初始模型可以为:
ΦP=-a ln(P)+b
其中,ΦP表示目标储层在上覆压力P下的覆压孔隙度,P表示目标储层的上覆压力,a表示待确定的压力项系数,b表示模型常数。
第四步,确定覆压孔隙度初始模型中的a值和b值。
首先,对于每份岩心样品,基于所建立的覆压孔隙度初始模型,根据每份岩心样品在常压下的常压孔隙度以及在高于常压的6个不同上覆压力下的覆压孔隙度,利用曲线拟合的方法,确定与该份岩心样品对应的覆压孔隙度模型。其中,与每份岩心样品对应的覆压孔隙度模型包括与该份岩心样品对应的压力项系数和模型常数(分别如表1中最右侧两列的A值和B值)。
在本实施例中,我们得到的与每份岩心样品对应的覆压孔隙度模型如下:
1号样品模型:ΦP=-0.511ln(P)+4.3404,R2=0.8508;
2号样品模型:ΦP=-0.379ln(P)+3.5884,R2=0.9902;
3号样品模型:ΦP=-0.329ln(P)+3.3392,R2=0.9359;
4号样品模型:ΦP=-0.322ln(P)+3.5331,R2=0.9141;
5号样品模型:ΦP=-0.69ln(P)+4.4064,R2=0.9964;
6号样品模型:ΦP=-0.483ln(P)+6.6533,R2=0.9901;
7号样品模型:ΦP=-0.754ln(P)+7.7577,R2=0.9752;
8号样品模型:ΦP=-0.504ln(P)+5.3064,R2=0.9978;
9号样品模型:ΦP=-0.322ln(P)+6.5015,R2=0.9986;
10号样品模型:ΦP=-0.172ln(P)+6.7673,R2=0.7613;
11号样品模型:ΦP=-0.325ln(P)+6.6521,R2=0.9906;
12号样品模型:ΦP=-0.625ln(P)+5.9842,R2=0.9488。
其中,R2表示相关系数。
接下来,根据与每份岩心样品对应的压力项系数(即表1中的A值)确定覆压孔隙度初始模型中待确定的压力项系数(即a值),根据每份岩心样品的常压孔隙度(Φ0)和与每份岩心样品对应的模型常数(即表1中的B值)确定覆压孔隙度初始模型中的模型常数(即b值)。
其中,根据与每份岩心样品对应的压力项系数(即表1中的A值)确定覆压孔隙度初始模型中待确定的压力项系数(即a值),可以包括:覆压孔隙度初始模型中待确定的压力项系数等于每份岩心样品对应的压力项系数的均值。
具体的,12个样品的均值模型:ΦP=-0.451ln(P)+5.4025,R2=0.9978,则a=0.451。
其中,根据每份岩心样品的常压孔隙度和与每份岩心样品对应的模型常数确定覆压孔隙度初始模型中的模型常数,可以包括:根据每份岩心样品的常压孔隙度和与每份岩心样品对应的模型常数,利用线性拟合的方法,建立覆压孔隙度初始模型中的模型常数与常压孔隙度之间的对应关系;基于所述对应关系,根据目标储层的常压孔隙度确定目标储层的覆压孔隙度初始模型中的模型常数。
对每份岩心样品的常压孔隙度和其对应的模型常数进行线性拟合,可以得到:b=0.8157Φ0+0.1184,R2=0.8766。
最后,将确定了压力项系数和模型常数的覆压孔隙度初始模型作为目标储层的覆压孔隙度模型。也即,目标储层的覆压孔隙度模型为:
ΦP=-0.451ln(P)+0.8157Φ0+0.1184
该模型即为目标储层的覆压孔隙度与上覆压力之间的对数关系模型。
当测得目标储层的上覆压力和常压孔隙度时,可以基于所构建的目标储层的覆压孔隙度模型,根据目标储层的上覆压力和常压孔隙度,确定目标储层的覆压孔隙度。
图4为根据本申请一具体实施例的岩心样品在2MPa上覆压力下的实测孔隙度与通过覆压孔隙度模型计算所得覆压孔隙度的对比图。图5为根据本申请一具体实施例的岩心样品在41.6MPa上覆压力下的实测孔隙度与通过覆压孔隙度模型计算所得覆压孔隙度的对比图。
通过图4和图5可以看出,对目标储层实测得到的孔隙度与通过所构建的覆压孔隙度模型所得到的的覆压孔隙度高度吻合。
本发明提供一种覆压孔隙度的确定方法,能够快速准确地计算不同压力下目标储层的覆压孔隙度,为页岩气藏评价提供数据基础,对于评价页岩储层及页岩油气储量有重要的意义。并可以用于我国南方页岩气勘探开发评价中确定覆压条件下的孔隙度。
实施例四
本实施例提供一种覆压孔隙度的确定装置,包括:模型构建模块,用于构建目标储层的覆压孔隙度模型,所述目标储层的覆压孔隙度模型为目标储层的覆压孔隙度与上覆压力之间的对数关系模型;孔隙度确定模块,用于基于目标储层的覆压孔隙度模型,根据目标储层的上覆压力和常压孔隙度,确定目标储层的覆压孔隙度。
在本实施例中,覆压孔隙度的确定装置还可以包括:处理器和存储器,其中所述处理器用于执行存储在存储器中的以下程序模块:模型构建模块和孔隙度确定模块,以目标储层的上覆压力和常压孔隙度,确定目标储层的覆压孔隙度。
实施例五
本实施例提供一种存储介质,存储有计算机程序,其特征在于,所述计算机程序被处理器执行时,实现如上文所述的覆压孔隙度的确定方法的步骤:
构建目标储层的覆压孔隙度模型,所述目标储层的覆压孔隙度模型为目标储层的覆压孔隙度与上覆压力之间的对数关系模型;
基于目标储层的覆压孔隙度模型,根据目标储层的上覆压力和常压孔隙度,确定目标储层的覆压孔隙度。
本领域内的技术人员应明白,本发明的实施例可提供为方法或计算机程序产品。因此,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法和计算机程序产品的流程图来描述的。应理解可由计算机程序指令实现流程图中的每一流程以及流程图中的流程的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程中指定的功能的步骤。
其中,存储介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。
存储介质的例子包括但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。
实施例六
本实施例提供一种计算机设备,包括处理器和存储有程序代码的存储介质,所述程序代码被所述处理器执行时,实现如上文所述的覆压孔隙度的确定方法的步骤:
构建目标储层的覆压孔隙度模型,所述目标储层的覆压孔隙度模型为目标储层的覆压孔隙度与上覆压力之间的对数关系模型;
基于目标储层的覆压孔隙度模型,根据目标储层的上覆压力和常压孔隙度,确定目标储层的覆压孔隙度。
在一个示例中,计算机设备可以包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash FLASH RAM)。内存是计算机可读介质的示例。
需要注意的是,这里所使用的的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
需要说明的是,本申请的说明书和权利要求书及附图中的术语“第一”、“第二”等是用于区别类似的对象,而不用于描述特定的顺序或先后次序。应该理解这样使用的术语在适当情况下可以互换,以便这里描述的本申请的实施方式例如能够除了在这里图示或描述的那些以外的顺序实施。
应当理解的是,本说明书中的示例性实施方式可以由多种不同的形式来实施,并且不应当被解释为只限于这里所阐述的实施方式。提供这些实施方式是为了使得本申请的公开彻底且完整,并且将这些示例性实施方式的构思充分传达给本领域普通技术人员,而不应当理解为对本发明的限制。

Claims (10)

1.一种覆压孔隙度的确定方法,其特征在于,包括以下步骤:
构建目标储层的覆压孔隙度模型,所述目标储层的覆压孔隙度模型为目标储层的覆压孔隙度与上覆压力之间的对数关系模型;
基于目标储层的覆压孔隙度模型,根据目标储层的上覆压力和常压孔隙度,确定目标储层的覆压孔隙度。
2.根据权利要求1所述的覆压孔隙度的确定方法,其特征在于,通过以下步骤构建覆压孔隙度模型:
获取若干份目标储层的岩心样品,对于每份岩心样品,依次施加常压以及高于常压的多个不同上覆压力,测量其在常压下的常压孔隙度以及在高于常压的各个不同上覆压力下的覆压孔隙度;
建立目标储层的覆压孔隙度初始模型,所述覆压孔隙度初始模型包括待确定的压力项系数和模型常数;
对于每份岩心样品,基于覆压孔隙度初始模型,根据其在常压下的常压孔隙度以及在高于常压的各个不同上覆压力下的覆压孔隙度,利用曲线拟合的方法,确定与该份岩心样品对应的覆压孔隙度模型,其中,与该份岩心样品对应的覆压孔隙度模型包括与该份岩心样品对应的压力项系数和模型常数;
根据与每份岩心样品对应的压力项系数确定覆压孔隙度初始模型中待确定的压力项系数,根据每份岩心样品的常压孔隙度和与每份岩心样品对应的模型常数确定覆压孔隙度初始模型中的模型常数;
将确定了压力项系数和模型常数的覆压孔隙度初始模型作为目标储层的覆压孔隙度模型。
3.根据权利要求2所述的覆压孔隙度的确定方法,其特征在于,根据与每份岩心样品对应的压力项系数确定覆压孔隙度初始模型中待确定的压力项系数,包括:
覆压孔隙度初始模型中待确定的压力项系数等于每份岩心样品对应的压力项系数的均值。
4.根据权利要求2所述的覆压孔隙度的确定方法,其特征在于,根据每份岩心样品的常压孔隙度和与每份岩心样品对应的模型常数确定覆压孔隙度初始模型中的模型常数,包括:
根据每份岩心样品的常压孔隙度和与每份岩心样品对应的模型常数,利用线性拟合的方法,建立覆压孔隙度初始模型中的模型常数与常压孔隙度之间的对应关系;
基于所述对应关系,根据目标储层的常压孔隙度确定目标储层的覆压孔隙度初始模型中的模型常数。
5.根据权利要求2所述的覆压孔隙度的确定方法,其特征在于,所述覆压孔隙度初始模型为:
ΦP=-aln(P)+b
其中,ΦP表示目标储层在上覆压力P下的覆压孔隙度,P表示目标储层的上覆压力,a表示待确定的压力项系数,b表示模型常数。
6.根据权利要求4所述的覆压孔隙度的确定方法,其特征在于,覆压孔隙度初始模型中的模型常数与常压孔隙度之间的关系为以下线性关系:
b=cΦ0+d
其中,b表示模型常数,c表示斜率,Φ0表示目标储层的常压孔隙度,d表示常数。
7.根据权利要求2所述的覆压孔隙度的确定方法,其特征在于,所获取的目标储层的岩心样品的份数至少为12份,对于每份岩心样品,所施加的高于常压的不同上覆压力至少为6个。
8.一种覆压孔隙度的确定装置,其特征在于,包括:
模型构建模块,用于构建目标储层的覆压孔隙度模型,所述目标储层的覆压孔隙度模型为目标储层的覆压孔隙度与上覆压力之间的对数关系模型;
孔隙度确定模块,用于基于目标储层的覆压孔隙度模型,根据目标储层的上覆压力和常压孔隙度,确定目标储层的覆压孔隙度。
9.一种存储介质,存储有计算机程序,其特征在于,所述计算机程序被处理器执行时,实现如权利要求1至7中任一项所述的覆压孔隙度的确定方法的步骤。
10.一种计算机设备,包括处理器和存储有程序代码的存储介质,所述程序代码被所述处理器执行时,实现如权利要求1至7中任一项所述的覆压孔隙度的确定方法的步骤。
CN202010848459.2A 2020-08-21 2020-08-21 覆压孔隙度的确定方法、装置、存储介质和计算机设备 Pending CN114077784A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010848459.2A CN114077784A (zh) 2020-08-21 2020-08-21 覆压孔隙度的确定方法、装置、存储介质和计算机设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010848459.2A CN114077784A (zh) 2020-08-21 2020-08-21 覆压孔隙度的确定方法、装置、存储介质和计算机设备

Publications (1)

Publication Number Publication Date
CN114077784A true CN114077784A (zh) 2022-02-22

Family

ID=80282200

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010848459.2A Pending CN114077784A (zh) 2020-08-21 2020-08-21 覆压孔隙度的确定方法、装置、存储介质和计算机设备

Country Status (1)

Country Link
CN (1) CN114077784A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114624163A (zh) * 2022-04-18 2022-06-14 成都理工大学 基于覆压孔渗测试的富浊沸石碎屑岩储层物性校正方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114624163A (zh) * 2022-04-18 2022-06-14 成都理工大学 基于覆压孔渗测试的富浊沸石碎屑岩储层物性校正方法
CN114624163B (zh) * 2022-04-18 2023-05-05 成都理工大学 基于覆压孔渗测试的富浊沸石碎屑岩储层物性校正方法

Similar Documents

Publication Publication Date Title
WO2009034253A1 (fr) Procede, programme et systeme informatique de mise a l'echelle de donnees de modele de reservoir d'hydrocarbure
WO2009027598A1 (fr) Procede, programme et systeme informatique de conciliation de donnees de modele de reservoir d'hydrocarbure
CN104272140A (zh) 用于校准在储藏层建模中使用的渗透性的系统和方法
MX2014004885A (es) Metodo para determinar en tiempo real la porosidad y la saturacion de agua de una formacion subterranea usando datos de registro de gas y perforacion.
CN107956465B (zh) 基于关联井的全区多井测井曲线标准化方法及装置
Yang et al. Reservoir fluid data acquisition using advanced mud logging gas in shale reservoirs
Askari et al. A fully integrated method for dynamic rock type characterization development in one of Iranian off-shore oil reservoir
Yang et al. Predicting reservoir fluid properties from advanced mud gas data
CN110470581B (zh) 确定储层应力敏感程度的方法、装置及存储介质
WO2017003828A1 (en) Oilfield reservior saturation and permeability modeling
CN114077784A (zh) 覆压孔隙度的确定方法、装置、存储介质和计算机设备
Moret‐Fernández et al. A modified multiple tension upward infiltration method to estimate the soil hydraulic properties
CN113552036B (zh) 中低孔渗储层油水相对渗透率的确定方法及装置
CN111241652B (zh) 一种确定地层原油粘度的方法及装置
Yang et al. An integrated approach to uncertainty assessment for coalbed methane model
Albrecht et al. Laboratory measurements of fluid transport properties on tight gas sandstones and applications
CN114486620B (zh) 一种黏土吸附气含量的确定方法、存储介质和计算机设备
CN112711076B (zh) 提取石油钻井中泥浆侵入地层深度的方法和装置
Amanipoor Providing a subsurface reservoir quality maps in oil fields by geostatistical methods
Dastan et al. A New Look at Nonlinear Regression in Well Test Interpretation
Worthington Quality-assured evaluation of effective porosity using fit-for-purpose estimates of clay-mineral volume fraction
CN108647417A (zh) 一种确定页岩气储层含气饱和度的简易方法
CN112328953B (zh) 气井的水侵识别方法、装置、设备及可读存储介质
Friedel et al. Optimum parameterization for saturation modeling in heterogeneous hydrocarbon reservoirs using well log inversion
Dubost et al. Automated hydraulic units, fluid types, and free-fluid levels resolved using new algorithms applicable to distributed pressure measurements

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