首页> 中国专利> 基于测井资料的曲线构建方法、装置和存储介质

基于测井资料的曲线构建方法、装置和存储介质

摘要

本申请提供了一种基于测井资料的曲线构建方法、装置和存储介质。该方法包括:通过实际测量的微电阻率成像测井资料计算得到孔隙度频率谱;对岩心压汞毛管压力曲线进行储层分类,得到每类储层的平均岩心压汞毛管压力曲线;基于岩心压汞毛管压力曲线、孔隙度频率谱以及储层分类,得到每类储层的平均孔隙度频率谱反向累积曲线;建立每类储层的平均岩心压汞毛管压力曲线和平均孔隙度频率谱反向累积曲线之间的刻度函数;结合孔隙度频率谱和每类储层的刻度函数分类构建深度连续的毛管压力曲线。实现了基于微电阻率成像测井资料定量评价储层孔隙结构的目的,对储层孔隙结构进行了更加精准的评价。

著录项

  • 公开/公告号CN112983394A

    专利类型发明专利

  • 公开/公告日2021-06-18

    原文格式PDF

  • 申请/专利权人 中国石油天然气股份有限公司;

    申请/专利号CN202110168757.1

  • 申请日2021-02-07

  • 分类号E21B47/00(20120101);E21B47/06(20120101);E21B49/00(20060101);

  • 代理机构11138 北京三高永信知识产权代理有限责任公司;

  • 代理人侯晓雅

  • 地址 100007 北京市东城区东直门北大街9号中国石油大厦

  • 入库时间 2023-06-19 11:29:13

说明书

技术领域

本申请涉及油气地质勘探与开发技术领域,特别涉及一种基于测井资料的曲线构建方法、装置和存储介质。

背景技术

缝洞型碳酸盐岩储层是油气勘探中的一类重要储层,该类型储层具有孔、洞、缝发育的特征,孔、缝、洞的发育是油气获得高产的关键。针对该类储层孔隙结构的研究是储层有效性评价的关键。

相关技术中,储层毛管压力曲线是用于定量计算孔喉半径、排驱压力等孔隙结构评价参数、进行储层孔隙结构评价的重要资料。但是由于取心数量有限,单井储层毛管压力曲线资料十分有限,只能利用离散的数据评价储层孔隙结构。此外,相关技术中还应用核磁共振测井资料评价孔隙结构,其中,核磁共振测井资料中的核磁T

所以相关技术中评价该类型储层孔隙结构的方法得到的结果存在缺陷,对储层孔隙结构进行定量评价的结果的精确度低。

发明内容

本申请提供了一种基于测井资料的曲线构建方法、装置和存储介质,提高了对缝洞型碳酸盐岩储层孔隙结构评价的准确性。

第一方面,提供了一种基于测井资料的曲线构建方法,该方法包括:

基于微电阻率成像测井资料计算孔隙度频率谱;对岩心压汞毛管压力曲线进行储层分类,基于分类后的岩心压汞毛管压力曲线,求取每类储层的平均岩心压汞毛管压力曲线;基于孔隙度频率谱和岩心压汞毛管压力曲线,结合储层分类,求取每类储层的平均孔隙度频率谱反向累积曲线;建立每类储层的平均岩心压汞毛管压力曲线和平均孔隙度频率谱反向累积曲线之间的刻度函数;结合孔隙度频率谱和每类储层的刻度函数分类构建深度连续的毛管压力曲线。

微电阻率成像测井资料中包括微电流成像,可以通过图像连续地展示储层的孔、缝、洞分布情况;此外,微电阻率成像测井资料显示的是更微观的测井资料,从而最终得到的深度连续的毛管压力曲线所反映出的信息也是更微观的,更有助于精准地评价储层的孔隙结构。

在一种可能的实现方式中,对岩心压汞毛管压力曲线进行储层分类,包括:结合地质分类标准,将岩心压汞毛管压力曲线根据对应的岩心孔隙度进行分类。

在一种可能的实现方式中,基于孔隙度频率谱和岩心压汞毛管压力曲线,结合储层分类,求取每类储层的平均孔隙度频率谱反向累积曲线,包括:从孔隙度频率谱中读取岩心压汞毛管压力曲线对应深度的孔隙度频率谱数据;对读取的孔隙度频率谱数据进行反向累积,结合储层分类,求取各类储层的平均孔隙度频率谱反向累积曲线。

在一种可能的实现方式中,基于孔隙度频率谱和岩心压汞毛管压力曲线,结合储层分类,求取每类储层的平均孔隙度频率谱反向累积曲线,包括:基于每类储层,从孔隙度频率谱中估选孔隙度频率谱数据;对估选的孔隙度频率谱数据进行反向累积,求取该每类储层的平均孔隙度频率谱反向累积曲线。

在一种可能的实现方式中,建立每类储层的平均岩心压汞毛管压力曲线和平均孔隙度频率谱反向累积曲线之间的刻度函数,包括:采用分段非线性转换刻度方法,建立每类储层的平均岩心压汞毛管压力曲线和平均孔隙度频率谱反向累积曲线之间的刻度函数。

在一种可能的实现方式中,结合所述孔隙度频率谱和所述每类储层的刻度函数分类构建深度连续的毛管压力曲线之后,还包括:根据构建的深度连续的毛管压力曲线,计算得到最大孔喉半径图像和孔喉半径分布图像。

在一种可能的实现方式中,该方法还包括:展示曲线构建的结果图像,该结果图像包括微电阻率成像资料中的微电阻率成像、孔隙度频率谱、构建的深度连续的毛管压力曲线和基础测井信息。

第二方面,提供了一种基于测井资料的曲线构建装置,该装置包括:第一计算模块,用于基于微电阻率成像测井资料计算孔隙度频率谱;第一求取模块,用于对岩心压汞毛管压力曲线进行储层分类,基于分类后的岩心压汞毛管压力曲线,求取每类储层的平均岩心压汞毛管压力曲线;第二求取模块,用于基于孔隙度频率谱和岩心压汞毛管压力曲线,结合储层分类,求取每类储层的平均孔隙度频率谱反向累积曲线;第一构建模块,用于建立每类储层的平均岩心压汞毛管压力曲线和平均孔隙度频率谱反向累积曲线之间的刻度函数;第二构建模块,用于结合所述孔隙度频率谱和所述每类储层的刻度函数分类构建深度连续的毛管压力曲线。

在一种可能的实现方式中,第一求取模块包括:分类单元,用于结合地质分类标准,将岩心压汞毛管压力曲线根据对应的岩心孔隙度进行分类。

在一种可能的实现方式中,第二求取模块包括:第一读取单元,用于从孔隙度频率谱中读取岩心压汞毛管压力曲线对应深度的孔隙度频率谱数据;第一求取单元,对读取的孔隙度频率谱数据进行反向累积,结合储层分类,求取各类储层的平均孔隙度频率谱反向累积曲线。

在一种可能的实现方式中,第二求取模块包括:第二读取单元,用于基于每类储层,从孔隙度频率谱中估选孔隙度频率谱数据;第二求取单元,用于对估选的孔隙度频率谱数据进行反向累积,结合储层分类,求取各类储层的平均孔隙度频率谱反向累积曲线。

在一种可能的实现方式中,第一构建模块用于采用分段非线性转换刻度方法,建立每类储层的平均岩心压汞毛管压力曲线和平均孔隙度频率谱反向累积曲线之间的刻度函数。

在一种可能的实现方式中,该装置还包括:第二计算模块,用于根据构建的深度连续的毛管压力曲线,计算得到最大孔喉半径图像和孔喉半径分布图像。

在一种可能的实现方式中,该装置还包括:显示模块,用于展示曲线构建的结果图像,该结果图像包括微电阻率成像资料中的微电阻率成像、孔隙度频率谱、构建的深度连续的毛管压力曲线和基础储层信息。

第三方面,提供一种计算机可读存储介质,该计算机可读存储介质中存储有至少一条程序代码,程序代码由处理器加载并执行以使计算机实现如第一方面中任一项的基于测井资料的曲线构建方法。

第四方面,提供了一种计算机程序或计算机程序产品,该计算机程序或计算机程序产品中存储有至少一条计算机指令,该至少一条计算机指令由处理器加载并执行,以使计算机实现上述第一方面中任一项基于测井资料的曲线构建方法。

附图说明

图1是本申请实施例提供的一种基于测井资料的曲线构建方法的流程图;

图2是本申请实施例提供的一种岩心压汞毛管压力曲线;

图3是本申请实施例提供的一种平均岩心压汞毛管压力曲线;

图4是本申请实施例提供的一种平均孔隙度频率谱反向累积曲线;

图5是本申请实施例提供的一种毛管压力P

图6是本申请实施例提供的一种微电阻率成像构建深度连续的毛管压力曲线结果;

图7是本申请实施例提供的一种基于测井资料的曲线构建装置的示意图。

具体实施方式

对于储层勘探工作,准确评价储层的孔隙结构、孔洞分布的情况是对储层进行分类、进行后续油气开采工作的基础,相关技术中的曲线构造方法主要基于核磁共振测井资料计算毛管压力曲线,由于核磁共振测井对碳酸盐储层的适应性较差,且核磁共振测井资料只能反映宽度较大的缝洞,导致最终得到的孔隙结构的评价结果的精准性较差。对此,本申请实施例提供了基于微电阻率成像测井资料的毛管压力曲线的构建方法。参见图1,该方法包括但不限于以下步骤101-步骤105:

101、基于微电阻率成像测井资料计算孔隙度频率谱。

该微电阻率成像测井资料通过专业的测井技术部门实际测量得到,微电阻率测井可以利用相应的电极结构,在不同深度对井中情况进行探测。微电阻率成像测井资料中包括微电阻率成像,可以图像形式展示各深度对应的储层中孔、缝、洞的分布情况,该微电阻率成像的图像是深度连续的,图像中的像素点对应于储层中的实际位置的电阻率数据。

基于该实际测量的微电阻率成像测井资料的数据,计算得到孔隙度频率谱,计算方式可以包括对阿尔奇公式的应用等多种计算方法,本申请实施例对此不做限定。另外,计算得到对应深度的孔隙度曲线,孔隙度曲线反映各深度的储层对应的孔隙度(Porosity,POR)的大小,将孔隙度曲线作为对孔隙度频率谱进行储层分类的依据;还将该孔隙度曲线作为基础储层信息,对储层孔隙度的数值进行基本展示。其中,孔隙度曲线和孔隙度频率谱中的曲线是深度连续的。

102、对岩心压汞毛管压力曲线进行储层分类,基于分类后的岩心压汞毛管压力曲线,求取每类储层的平均岩心压汞毛管压力曲线。

在储层勘探的过程中,会通过取心作业采集不同深度的多个岩心样本,对岩心样本进行压汞实验,获得每个岩心样本对应的岩心压汞毛管压力曲线;此外,对每个岩心样本进行其他的实验测算,获取到岩心样本对应的储层的深度、岩心样本的孔隙度、岩心样本的排驱压力等储集空间标量参数。因而,每条岩心压汞毛管压力曲线也对应确定的深度和孔隙度。

依据被勘探的储层相对应的孔隙度和地质分类标准,对岩心压汞毛管压力曲线进行储层分类,该地质分类标准是依据储层的孔隙度将储层分为不同的类型的分类标准。依据分类结果,计算每一类储层的平均岩心压汞毛管压力曲线。该平均岩心压汞毛管压力曲线反映每类岩心样本的进汞饱和度(Sh

103、基于孔隙度频率谱和岩心压汞毛管压力曲线,结合储层分类,求取前述每类储层的平均孔隙度频率谱反向累积曲线。

由于岩心压汞毛管压力曲线具有对应的孔隙度和深度,所以基于结合102中实验得到的岩心压汞毛管压力曲线,在101中计算得到的孔隙度频率谱中能够得到与岩心压汞毛管压力曲线对应的孔隙度频率谱数据,并计算其反向累积饱和度;根据前述的进行储层分类的地质分类标准,进一步计算求取前述步骤102中的每类储层的平均孔隙度谱反向累积曲线。该平均孔隙度频率谱反向累积曲线反映前述每一类储层的反向累积饱和度(S

场景一:岩心压汞毛管压力曲线对应深度具有微电阻率成像测井资料。

微电阻率成像测井资料是通过图像反映储层的孔隙结构,会对应不同深度的储层。在该种可能的实现场景中,进行压汞实验的岩心样本对应的深度同时进行了微电阻率成像技术的测量,则该岩心压汞毛管压力曲线对应深度具有微电阻率成像测井资料。

在该种可能的实现方式中,读取每条岩心压汞毛管压力曲线对应深度点的孔隙度频率谱数据,对其进行反向累积,计算得到反向累积饱和度,用以得到多条孔隙度频率谱反向累积曲线,孔隙度频率谱反向累积曲线反映反向累积饱和度(S

场景二:岩心压汞毛管压力曲线对应深度不具有微电阻率成像测井资料。

在该种可能的实现方式中,因为测井资料的缺少,无法获取每一条岩心压汞毛管压力曲线对应深度点的孔隙度频率谱数据,则基于每一类储层分别估选孔隙度频率谱数据将其进行反向累积,进而求取每一类储层的平均孔隙度频率谱反向累积曲线。其中,估选孔隙度频率谱数据的方式可以是,根据岩心压汞毛管压力曲线对应的孔隙度,在孔隙度频率谱中选取体现的孔隙度与该岩心压汞毛管压力曲线的孔隙度相近的孔隙度频率谱数据。其中,孔隙度相近是指岩心压汞毛管压力曲线的孔隙度与孔隙度频率谱中的某一数据体现的孔隙度相同或者偏差程度小于一定阈值。该阈值可由技术人员预先设定或者根据实际环境情况进行调整。

104、建立每类储层的平均岩心压汞毛管压力曲线和平均孔隙度频率谱反向累积曲线之间的刻度函数。

依据102中的平均岩心压汞毛管压力曲线和103中的平均孔隙度频率谱反向累积曲线,采用分段非线性转换刻度方法,建立每类储层的平均岩心压汞毛管压力曲线和平均孔隙度频率谱反向累积曲线之间的刻度函数。该步骤包括但不限于1041和1042两个分步骤:

1041、根据饱和度相等的原则,针对每一类储层,在含水饱和度与进汞饱和度相同时,确定相同饱和度条件下的孔隙度值分量的倒数1/POR值和毛管压力P

1042、根据1041中得到的图像显示的每一类储层毛管压力P

105、结合孔隙度频率谱和每类储层的刻度函数分类构建深度连续的毛管压力曲线。

在本步骤中,最终毛管压力曲线的构建过程包括但不限于以下步骤1051~1053:

1051、依据前述地质分类标准,基于101中计算得到的孔隙度曲线,对101中计算得到的孔隙度频率谱进行储层分类,对每一类储层,应用104中已建立的该储层的毛管压力P

1052、根据在101中计算得到的孔隙度频率谱数据,进行反向累积,求取得到反向累积饱和度S

1053、将求取的毛管压力P

在本步骤中,由于步骤1051中基于深度连续的孔隙度频率谱进行储层分类和计算,且步骤1051得到的毛管压力P

在一种可能的实现方式中,将得到的每类储层深度连续的毛管压力曲线构建在同一图像中进行展示。

进一步的,根据构建得到的深度连续的毛管压力曲线可以计算得到最大孔喉半径和孔喉半径分布的连续图像。

在一种可能的实现方式中,展示曲线构建的结果图像,该结果图像包括微电阻率成像资料中的微电阻率成像、孔隙度频率谱、构建的深度连续的毛管压力曲线和基础储层信息。即最终将微电阻率成像测井资料、构建的深度连续的毛管压力曲线以及计算得到的最大孔喉半径和孔喉半径分布的连续图像等图像以及基础储层信息进行组合显示。基础储层信息包括从测井资料得到的井径、声波时差等信息,本申请实施例对此不作限定。

本申请实施例提供了一种基于微电阻成像测井资料的毛管压力曲线构建方法。通过实际测量的微电阻率成像测井资料计算得到孔隙度频率谱,再结合岩心压汞毛管压力曲线和储层的分类标准,建立平均岩心压汞毛管压力曲线和平均孔隙度频率谱反向累积曲线之间的刻度函数,建立起孔隙度频率谱与毛管压力曲线之间的关系,将孔隙度频率谱转换为毛管压力曲线,以实现微电阻率成像测井资料定量评价储层孔隙结构的目的,更精准地评价了储层的孔隙结构。

下面将结合示例性实施例介绍本申请提供的方法:

201、基于微电阻率成像测井资料计算得到孔隙度频率谱。

从专业的测井技术部门获取到实际测量的微电阻率成像测井资料,该实际测量的微电阻率成像测井资料包括以图像形式展示的微电阻率成像,直观展示储层中孔、缝、洞的分布情况。示例性地,该微电阻率成像中像素点的亮度反映出储层中对应深度、对应位置的电阻率,像素点呈暗色,则该点对应的电阻率较低,则实际储层中对应位置的孔、缝、洞发育情况较好,利于油气的存储,也就意味着该储层为有利储层的概率较高,开发价值较高。

基于实际测量获得的微电阻率成像测井资料的内容,利用图像中像素点对应的电阻率进行计算,得到孔隙度频率谱,还计算得到对应深度的孔隙度曲线。其中,孔隙度频率谱反映微电阻率成像对应深度的孔隙度的分布情况,孔隙度曲线反映对应深度的储层的孔隙度的大小;该步骤中计算时可以应用阿尔齐公式。孔隙度曲线则作为对孔隙度频率谱进行储层分类的依据,以及用于对储层孔隙度情况的基本展示。其中,孔隙度曲线和孔隙度频率谱中的曲线是深度连续的。

202、对岩心压汞毛管压力曲线进行储层分类。

在储层勘探过程中,不只获取数据、图像等测井资料,也会通过取心作业获取一定数量的岩心样本。对采集到的样本进行实验分析,获取岩心样本的储集空间标量参数。示例性地,获取每个岩心样本对应的实际深度、测量计算每个样本的孔隙度等。

在本申请实施例中,还对每个岩心样本进行压汞实验,得到每个岩心样本的岩心压汞毛管压力曲线,用以反映每个岩心样本的进汞饱和度(Sh

在本申请中,针对示例性实施例中获取到的碳酸盐储层的多个岩心样本,测量对应的孔隙度,然后进行压汞实验,得到多条岩心压汞毛管压力曲线,并结合对应的地质分类标准(参见表1),对得到的岩心压汞毛管压力曲线按照地质分类标准加以分类,得到如2所示的岩心压汞毛管压力曲线图,图中相同颜色的曲线对应的岩心样本属于相同类型的储层。

表1地质分类标准

203、基于储层分类,求取每类储层的平均岩心压汞毛管压力曲线。

根据对岩心压汞毛管压力曲线分类后得到的图2,对每一类岩心压汞毛管压力曲线进行平均,分别得到Ⅰ类、Ⅱ类、Ⅲ类储层的平均岩心压汞毛管压力曲线,如图3所示,反映每一类储层的进汞饱和度(Sh

204、基于孔隙度频率谱和岩心压汞毛管压力曲线,结合储层分类,求取前述每类储层的平均孔隙度频率谱反向累积曲线。

应用201中得到的孔隙度频率谱和202中得到的岩心压汞毛管压力曲线,结合202中的储层分类情况,求取与步骤202中储层对应的每类储层的平均孔隙度频率谱反向累积曲线。该曲线反映每一类储层的反向累积饱和度(S

在一种可能的实现方式中,岩心压汞毛管压力曲线图中各曲线的对应深度具有微电阻率成像测井资料,则从201中得到的孔隙度频率谱中读取对应深度点的孔隙度频率谱数据,对获取的孔隙度频率谱数据进行反向累积,得到与孔隙度频率谱数据对应的反向累积饱和度,用以得到孔隙度频率谱反向累积曲线;并根据202中的储层分类标准,对得到孔隙度频率谱反向累积曲线进行分类;对每类孔隙度频率谱反向累积曲线求平均,得到平均孔隙度频率谱反向累积曲线。

在一种可能的实现方式中,岩心压汞毛管压力曲线图中各曲线对应深度未进行微电阻率成像测井时,针对每一类储层分别估选孔隙度频率谱数据,将其进行反向累积,得到孔隙度频率谱反向累积曲线,进而求取每类储层的平均孔隙度频率谱反向累积曲线。其中,估选孔隙度频率谱数据的一种可选方式为:根据岩心压汞毛管压力曲线图中各曲线对应的孔隙度,在孔隙度频率谱中选取对应的孔隙度与各曲线对应的孔隙度相近的孔隙度频率谱数据。示例性地,当岩心压汞毛管压力曲线对应的孔隙度的大小和孔隙度频率谱中某一孔隙度频率谱数据体现的孔隙度大小相同或者偏差在一定阈值内,则选取该孔隙度频率谱数据。其中,阈值的数值可以有技术人员预先设定,也可以根据储层的地质情况进行调整。

在本申请示例性实施例中,岩心压汞毛管压力曲线图中各曲线对应深度具有微电阻率成像测井资料,采用相应的计算方法,得到了如图4所示的平均孔隙度频率谱反向累积曲线。

205、建立每类储层的平均岩心压汞毛管压力曲线和平均孔隙度频率谱反向累积曲线之间的刻度函数。

本步骤中,采用分段非线性转换刻度方法,建立每类储层的岩心压汞平均毛管压力曲线和平均孔隙度频率谱反向累积曲线之间的刻度函数。具体步骤包括但不限于如下2051~2052:

2051、依据204中获得的平均孔隙度频率谱反向累积曲线和203中获得的平均岩心压汞毛管压力曲线,根据饱和度相等的原则,针对每一类储层,在含水饱和度与进汞饱和度相同时,确定相同饱和度条件下的孔隙度值分量的倒数1/POR的值和毛管压力P

2052、根据2051中得到的图像显示的每一类储层毛管压力P

在大孔段,该刻度函数的数学形式为:

式中:P

在小孔段,该刻度函数的数学形式为:

式中:P

其中,大孔段和小孔段的划分依据人员经验和1/POR和P

在本申请示例性实施例中,针对缝洞型碳酸盐储层获得的1/POR和P

206、结合孔隙度频率谱和每类储层的刻度函数分类构建深度连续的毛管压力曲线。

该步骤可以依据以下2061~2063步骤进行。

2061、结合前述的地质分类标准,依据201中计算得到的孔隙度曲线,对201中计算得到的孔隙度频率谱进行储层分类。对分类得到的每一类储层,应用205中已建立的毛管压力P

2062、根据在201中计算得到的孔隙度频率谱数据,进行反向累积,求取得到反向累积饱和度S

2063、将步骤2061中求取的每类储层对应的多个毛管压力P

在示例性实施例中,依据地质分类标准将储层分为三类,得到图5表示的三类储层的刻度函数,且每类储层中的大、小孔存在各自相应的刻度函数。在将孔隙度频率谱进行储层分类后,依据图5计算得到每类储层对应的多个毛管压力,结合相应的反向累积饱和度分别构建每类储层的深度连续的毛管压力曲线。

在一种可能的实现方式中,将每类储层的深度连续的毛管压力曲线在同一图像中显示。

进一步的,根据构建得到的深度连续的毛管压力曲线可以计算得到最大孔喉半径和孔喉半径分布的连续图像。

进一步的,展示曲线构建的结果图像,该结果图像包括微电阻率成像资料中的微电阻率成像、孔隙度频率谱、构建的深度连续的毛管压力曲线和基础储层信息。最终将微电阻率成像测井资料、构建的深度连续的毛管压力曲线以及计算分析得到的最大孔喉半径和孔喉半径分布的连续图像等图像和基础储层信息进行组合显示。

在本申请示例性实施例中,基于微电阻率成像测井资料构建得到的深度连续的毛管压力曲线的结果图像如图6所示。图6中第1-4道为基础储层信息展示图像,依次包括第一道的井径曲线(borehole diameter,CAL)、自然伽马曲线(natural gamma ray,GR),第二道的深度(DEPTH),第三道的声波时差曲线(acoustic,AC)、中子曲线(neutron,CNL)、密度曲线(density,DEN),第四道的地层真电阻率(true formation resistivity,RT)、冲洗带电阻率曲线(flushed zone formation resistivity,Rxo);第5道为储层分类,其中,RL1、RL2、RL3分别对应Ⅰ类储层、Ⅱ类储层、Ⅲ类储层,第6道为孔隙度(Porosity,POR)曲线,第7道为微电阻率动态成像(IMAGEDYNAMIC),第8道为孔隙度频率谱(PHIDIST),第9道为应用该方法构建得到的深度连续的毛管压力曲线(PCDIST),第10道为应用构造的深度连续的毛管压力曲线计算得到最大孔喉半径(RCMAX),第11道为孔喉半径分布(RCDIST)。

本申请具体实施例提供了一种基于微电阻成像的缝洞型碳酸盐岩储层毛管压力曲线构建方法。通过实际测量的微电阻率成像测井资料计算得到孔隙度频率谱,通过建立平均岩心压汞毛管压力曲线和平均孔隙度频率谱反向累积曲线之间的刻度函数,建立起孔隙度频率谱与毛管压力曲线之间的关系,将孔隙度频率谱转换为毛管压力曲线,以实现应用微电阻率成像测井资料定量评价缝洞型碳酸盐储层孔隙结构的目的,更精准地了解储层的孔隙结构。

参见图7,本申请实施例提供了一种基于测井资料的曲线构建装置,该装置包括:

第一计算模块701,用于基于微电阻率成像测井资料计算孔隙度频率谱;第一求取模块702,用于对岩心压汞毛管压力曲线进行储层分类,基于分类后的岩心压汞毛管压力曲线,求取每类储层的平均岩心压汞毛管压力曲线;第二求取模块703,用于基于孔隙度频率谱和岩心压汞毛管压力曲线,结合储层分类,求取每类储层的平均孔隙度频率谱反向累积曲线;第一构建模块704,用于建立每类储层的平均岩心压汞毛管压力曲线和平均孔隙度频率谱反向累积曲线之间的刻度函数;第二构建模块705,用于结合孔隙度频率谱和每类储层的刻度函数分类构建深度连续的毛管压力曲线。

在一种可能的实现方式中,第一求取模块702包括:分类单元,用于结合地质分类标准,将岩心压汞毛管压力曲线根据对应的岩心孔隙度进行分类。

在一种可能的实现方式中,第二求取模块703包括:第一读取单元,用于从孔隙度频率谱中读取岩心压汞毛管压力曲线对应深度的孔隙度频率谱数据;第一求取单元,对读取的孔隙度频率谱数据进行反向累积,结合储层分类,求取各类储层的平均孔隙度频率谱反向累积曲线。

在一种可能的实现方式中,第二求取模块703包括:第二读取单元,用于基于每类储层,从孔隙度频率谱中估选孔隙度频率谱数据;第二求取单元,用于对估选的孔隙度频率谱数据进行反向累积,结合储层分类,求取各类储层的平均孔隙度频率谱反向累积曲线。

在一种可能的实现方式中,第一构建模块704,用于采用分段非线性转换刻度方法,建立每类储层的平均岩心压汞毛管压力曲线和平均孔隙度频率谱反向累积曲线之间的刻度函数。

在一种可能的实现方式中,该装置还包括:第二计算模块,用于根据构建的深度连续的毛管压力曲线,计算得到最大孔喉半径图像和孔喉半径分布图像。

在一种可能的实现方式中,该装置还包括:显示模块,用于展示曲线构建的结果图像,该结果图像包括微电阻率成像资料中的微电阻率成像、孔隙度频率谱、构建的深度连续的毛管压力曲线和基础储层信息。

应理解的是,上述图7提供的装置在实现其功能时,仅以上述各功能模块的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的功能模块完成,即将设备的内部结构划分成不同的功能模块,以完成以上描述的全部或者部分功能。另外,上述实施例提供的装置与方法实施例属于同一构思,其具体实现过程详见方法实施例,这里不再赘述。

在示例性实施例中,还提供了一种计算机可读存储介质,该存储介质中存储有至少一条程序代码,该至少一条程序代码由处理器加载并执行,以使计算机实现上述任一种基于测井资料的曲线构建方法。

可选地,上述计算机可读存储介质可以是只读存储器(Read-Only Memory,ROM)、随机存取存储器(Random Access Memory,RAM)、只读光盘(Compact Disc Read-OnlyMemory,CD-ROM)、磁带、软盘和光数据存储设备等。

在示例性实施例中,还提供了一种计算机程序或计算机程序产品,该计算机程序或计算机程序产品中存储有至少一条计算机指令,该至少一条计算机指令由处理器加载并执行,以使计算机实现上述任一种基于测井资料的曲线构建方法。

应理解,在本申请的各个实施例中,各个过程的序号的大小并不意味着执行顺序的先后,各过程的执行顺序应以其功能和内在逻辑确定,而不应对本申请实施例的实施过程构成任何限定。

本申请中术语“至少一个”的含义是指一个或多个,本申请中术语“多个”的含义是指两个或两个以上。

应理解,在本文中对各种所述示例的描述中所使用的术语只是为了描述特定示例,而并非旨在进行限制。如在对各种所述示例的描述和所附权利要求书中所使用的那样,单数形式“一个(“a”,“an”)”和“该”旨在也包括复数形式,除非上下文另外明确地指示。

还应理解,术语“包括”(也称“includes”、“including”、“comprises”和/或“comprising”)当在本说明书中使用时指定存在所陈述的特征、整数、步骤、操作、元素、和/或部件,但是并不排除存在或添加一个或多个其他特征、整数、步骤、操作、元素、部件、和/或其分组。

应理解,根据A确定B并不意味着仅仅根据A确定B,还可以根据A和/或其它信息确定B。

还应理解,术语“若”和“如果”可被解释为意指“当..时”(“when”或“upon”)或“响应于确定”或“响应于检测到”。类似地,根据上下文,短语“若确定...”或“若检测到[所陈述的条件或事件]”可被解释为意指“在确定...时”或“响应于确定...”或“在检测到[所陈述的条件或事件]时”或“响应于检测到[所陈述的条件或事件]”。

还应理解,说明书通篇中提到的“一个实施例”、“一实施例”、“一种可能的实现方式”意味着与实施例或实现方式有关的特定特征、结构或特性包括在本申请的至少一个实施例中。因此,在整个说明书各处出现的“在一个实施例中”或“在一实施例中”、“一种可能的实现方式”未必一定指相同的实施例。此外,这些特定的特征、结构或特性可以任意适合的方式结合在一个或多个实施例中。

以上所述仅为本申请的可选实施例,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本申请的保护范围之内。

去获取专利,查看全文>

相似文献

  • 专利
  • 中文文献
  • 外文文献
获取专利

客服邮箱:kefu@zhangqiaokeyan.com

京公网安备:11010802029741号 ICP备案号:京ICP备15016152号-6 六维联合信息科技 (北京) 有限公司©版权所有
  • 客服微信

  • 服务号