CN112983394A - 基于测井资料的曲线构建方法、装置和存储介质 - Google Patents
基于测井资料的曲线构建方法、装置和存储介质 Download PDFInfo
- Publication number
- CN112983394A CN112983394A CN202110168757.1A CN202110168757A CN112983394A CN 112983394 A CN112983394 A CN 112983394A CN 202110168757 A CN202110168757 A CN 202110168757A CN 112983394 A CN112983394 A CN 112983394A
- Authority
- CN
- China
- Prior art keywords
- reservoir
- frequency spectrum
- capillary pressure
- curve
- porosity
- 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
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/06—Measuring temperature or pressure
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本申请提供了一种基于测井资料的曲线构建方法、装置和存储介质。该方法包括:通过实际测量的微电阻率成像测井资料计算得到孔隙度频率谱;对岩心压汞毛管压力曲线进行储层分类,得到每类储层的平均岩心压汞毛管压力曲线;基于岩心压汞毛管压力曲线、孔隙度频率谱以及储层分类,得到每类储层的平均孔隙度频率谱反向累积曲线;建立每类储层的平均岩心压汞毛管压力曲线和平均孔隙度频率谱反向累积曲线之间的刻度函数;结合孔隙度频率谱和每类储层的刻度函数分类构建深度连续的毛管压力曲线。实现了基于微电阻率成像测井资料定量评价储层孔隙结构的目的,对储层孔隙结构进行了更加精准的评价。
Description
技术领域
本申请涉及油气地质勘探与开发技术领域,特别涉及一种基于测井资料的曲线构建方法、装置和存储介质。
背景技术
缝洞型碳酸盐岩储层是油气勘探中的一类重要储层,该类型储层具有孔、洞、缝发育的特征,孔、缝、洞的发育是油气获得高产的关键。针对该类储层孔隙结构的研究是储层有效性评价的关键。
相关技术中,储层毛管压力曲线是用于定量计算孔喉半径、排驱压力等孔隙结构评价参数、进行储层孔隙结构评价的重要资料。但是由于取心数量有限,单井储层毛管压力曲线资料十分有限,只能利用离散的数据评价储层孔隙结构。此外,相关技术中还应用核磁共振测井资料评价孔隙结构,其中,核磁共振测井资料中的核磁T2谱的形态及其获取的参数,只能定性评价储层孔隙结构;又因为核磁易受油气影响,对储层孔、缝、洞响应的敏感度较差,基于核磁共振测井资料中其他内容对储层孔隙结构进行定量评价的方法具有一定局限性。
所以相关技术中评价该类型储层孔隙结构的方法得到的结果存在缺陷,对储层孔隙结构进行定量评价的结果的精确度低。
发明内容
本申请提供了一种基于测井资料的曲线构建方法、装置和存储介质,提高了对缝洞型碳酸盐岩储层孔隙结构评价的准确性。
第一方面,提供了一种基于测井资料的曲线构建方法,该方法包括:
基于微电阻率成像测井资料计算孔隙度频率谱;对岩心压汞毛管压力曲线进行储层分类,基于分类后的岩心压汞毛管压力曲线,求取每类储层的平均岩心压汞毛管压力曲线;基于孔隙度频率谱和岩心压汞毛管压力曲线,结合储层分类,求取每类储层的平均孔隙度频率谱反向累积曲线;建立每类储层的平均岩心压汞毛管压力曲线和平均孔隙度频率谱反向累积曲线之间的刻度函数;结合孔隙度频率谱和每类储层的刻度函数分类构建深度连续的毛管压力曲线。
微电阻率成像测井资料中包括微电流成像,可以通过图像连续地展示储层的孔、缝、洞分布情况;此外,微电阻率成像测井资料显示的是更微观的测井资料,从而最终得到的深度连续的毛管压力曲线所反映出的信息也是更微观的,更有助于精准地评价储层的孔隙结构。
在一种可能的实现方式中,对岩心压汞毛管压力曲线进行储层分类,包括:结合地质分类标准,将岩心压汞毛管压力曲线根据对应的岩心孔隙度进行分类。
在一种可能的实现方式中,基于孔隙度频率谱和岩心压汞毛管压力曲线,结合储层分类,求取每类储层的平均孔隙度频率谱反向累积曲线,包括:从孔隙度频率谱中读取岩心压汞毛管压力曲线对应深度的孔隙度频率谱数据;对读取的孔隙度频率谱数据进行反向累积,结合储层分类,求取各类储层的平均孔隙度频率谱反向累积曲线。
在一种可能的实现方式中,基于孔隙度频率谱和岩心压汞毛管压力曲线,结合储层分类,求取每类储层的平均孔隙度频率谱反向累积曲线,包括:基于每类储层,从孔隙度频率谱中估选孔隙度频率谱数据;对估选的孔隙度频率谱数据进行反向累积,求取该每类储层的平均孔隙度频率谱反向累积曲线。
在一种可能的实现方式中,建立每类储层的平均岩心压汞毛管压力曲线和平均孔隙度频率谱反向累积曲线之间的刻度函数,包括:采用分段非线性转换刻度方法,建立每类储层的平均岩心压汞毛管压力曲线和平均孔隙度频率谱反向累积曲线之间的刻度函数。
在一种可能的实现方式中,结合所述孔隙度频率谱和所述每类储层的刻度函数分类构建深度连续的毛管压力曲线之后,还包括:根据构建的深度连续的毛管压力曲线,计算得到最大孔喉半径图像和孔喉半径分布图像。
在一种可能的实现方式中,该方法还包括:展示曲线构建的结果图像,该结果图像包括微电阻率成像资料中的微电阻率成像、孔隙度频率谱、构建的深度连续的毛管压力曲线和基础测井信息。
第二方面,提供了一种基于测井资料的曲线构建装置,该装置包括:第一计算模块,用于基于微电阻率成像测井资料计算孔隙度频率谱;第一求取模块,用于对岩心压汞毛管压力曲线进行储层分类,基于分类后的岩心压汞毛管压力曲线,求取每类储层的平均岩心压汞毛管压力曲线;第二求取模块,用于基于孔隙度频率谱和岩心压汞毛管压力曲线,结合储层分类,求取每类储层的平均孔隙度频率谱反向累积曲线;第一构建模块,用于建立每类储层的平均岩心压汞毛管压力曲线和平均孔隙度频率谱反向累积曲线之间的刻度函数;第二构建模块,用于结合所述孔隙度频率谱和所述每类储层的刻度函数分类构建深度连续的毛管压力曲线。
在一种可能的实现方式中,第一求取模块包括:分类单元,用于结合地质分类标准,将岩心压汞毛管压力曲线根据对应的岩心孔隙度进行分类。
在一种可能的实现方式中,第二求取模块包括:第一读取单元,用于从孔隙度频率谱中读取岩心压汞毛管压力曲线对应深度的孔隙度频率谱数据;第一求取单元,对读取的孔隙度频率谱数据进行反向累积,结合储层分类,求取各类储层的平均孔隙度频率谱反向累积曲线。
在一种可能的实现方式中,第二求取模块包括:第二读取单元,用于基于每类储层,从孔隙度频率谱中估选孔隙度频率谱数据;第二求取单元,用于对估选的孔隙度频率谱数据进行反向累积,结合储层分类,求取各类储层的平均孔隙度频率谱反向累积曲线。
在一种可能的实现方式中,第一构建模块用于采用分段非线性转换刻度方法,建立每类储层的平均岩心压汞毛管压力曲线和平均孔隙度频率谱反向累积曲线之间的刻度函数。
在一种可能的实现方式中,该装置还包括:第二计算模块,用于根据构建的深度连续的毛管压力曲线,计算得到最大孔喉半径图像和孔喉半径分布图像。
在一种可能的实现方式中,该装置还包括:显示模块,用于展示曲线构建的结果图像,该结果图像包括微电阻率成像资料中的微电阻率成像、孔隙度频率谱、构建的深度连续的毛管压力曲线和基础储层信息。
第三方面,提供一种计算机可读存储介质,该计算机可读存储介质中存储有至少一条程序代码,程序代码由处理器加载并执行以使计算机实现如第一方面中任一项的基于测井资料的曲线构建方法。
第四方面,提供了一种计算机程序或计算机程序产品,该计算机程序或计算机程序产品中存储有至少一条计算机指令,该至少一条计算机指令由处理器加载并执行,以使计算机实现上述第一方面中任一项基于测井资料的曲线构建方法。
附图说明
图1是本申请实施例提供的一种基于测井资料的曲线构建方法的流程图;
图2是本申请实施例提供的一种岩心压汞毛管压力曲线;
图3是本申请实施例提供的一种平均岩心压汞毛管压力曲线;
图4是本申请实施例提供的一种平均孔隙度频率谱反向累积曲线;
图5是本申请实施例提供的一种毛管压力Pc和孔隙度频率谱中各孔隙度值分量的倒数1/POR值之间的转换关系图;
图6是本申请实施例提供的一种微电阻率成像构建深度连续的毛管压力曲线结果;
图7是本申请实施例提供的一种基于测井资料的曲线构建装置的示意图。
具体实施方式
对于储层勘探工作,准确评价储层的孔隙结构、孔洞分布的情况是对储层进行分类、进行后续油气开采工作的基础,相关技术中的曲线构造方法主要基于核磁共振测井资料计算毛管压力曲线,由于核磁共振测井对碳酸盐储层的适应性较差,且核磁共振测井资料只能反映宽度较大的缝洞,导致最终得到的孔隙结构的评价结果的精准性较差。对此,本申请实施例提供了基于微电阻率成像测井资料的毛管压力曲线的构建方法。参见图1,该方法包括但不限于以下步骤101-步骤105:
101、基于微电阻率成像测井资料计算孔隙度频率谱。
该微电阻率成像测井资料通过专业的测井技术部门实际测量得到,微电阻率测井可以利用相应的电极结构,在不同深度对井中情况进行探测。微电阻率成像测井资料中包括微电阻率成像,可以图像形式展示各深度对应的储层中孔、缝、洞的分布情况,该微电阻率成像的图像是深度连续的,图像中的像素点对应于储层中的实际位置的电阻率数据。
基于该实际测量的微电阻率成像测井资料的数据,计算得到孔隙度频率谱,计算方式可以包括对阿尔奇公式的应用等多种计算方法,本申请实施例对此不做限定。另外,计算得到对应深度的孔隙度曲线,孔隙度曲线反映各深度的储层对应的孔隙度(Porosity,POR)的大小,将孔隙度曲线作为对孔隙度频率谱进行储层分类的依据;还将该孔隙度曲线作为基础储层信息,对储层孔隙度的数值进行基本展示。其中,孔隙度曲线和孔隙度频率谱中的曲线是深度连续的。
102、对岩心压汞毛管压力曲线进行储层分类,基于分类后的岩心压汞毛管压力曲线,求取每类储层的平均岩心压汞毛管压力曲线。
在储层勘探的过程中,会通过取心作业采集不同深度的多个岩心样本,对岩心样本进行压汞实验,获得每个岩心样本对应的岩心压汞毛管压力曲线;此外,对每个岩心样本进行其他的实验测算,获取到岩心样本对应的储层的深度、岩心样本的孔隙度、岩心样本的排驱压力等储集空间标量参数。因而,每条岩心压汞毛管压力曲线也对应确定的深度和孔隙度。
依据被勘探的储层相对应的孔隙度和地质分类标准,对岩心压汞毛管压力曲线进行储层分类,该地质分类标准是依据储层的孔隙度将储层分为不同的类型的分类标准。依据分类结果,计算每一类储层的平均岩心压汞毛管压力曲线。该平均岩心压汞毛管压力曲线反映每类岩心样本的进汞饱和度(Shg)和毛管压力(Pc)之间的函数关系,计算得到的平均岩心压汞毛管压力曲线反映每类储层的进汞饱和度(Shg)和毛管压力(Pc)之间的函数关系Shg=f(PC)。
103、基于孔隙度频率谱和岩心压汞毛管压力曲线,结合储层分类,求取前述每类储层的平均孔隙度频率谱反向累积曲线。
由于岩心压汞毛管压力曲线具有对应的孔隙度和深度,所以基于结合102中实验得到的岩心压汞毛管压力曲线,在101中计算得到的孔隙度频率谱中能够得到与岩心压汞毛管压力曲线对应的孔隙度频率谱数据,并计算其反向累积饱和度;根据前述的进行储层分类的地质分类标准,进一步计算求取前述步骤102中的每类储层的平均孔隙度谱反向累积曲线。该平均孔隙度频率谱反向累积曲线反映前述每一类储层的反向累积饱和度(Sw)和孔隙度频率谱中各孔隙度值分量的倒数1/POR值之间的函数关系Sw=f(1/POR)。反向累积饱和度可以反映储层的含水饱和度。其中,求取前述每类储层的平均孔隙度谱反向累积曲线有以下两种场景。
场景一:岩心压汞毛管压力曲线对应深度具有微电阻率成像测井资料。
微电阻率成像测井资料是通过图像反映储层的孔隙结构,会对应不同深度的储层。在该种可能的实现场景中,进行压汞实验的岩心样本对应的深度同时进行了微电阻率成像技术的测量,则该岩心压汞毛管压力曲线对应深度具有微电阻率成像测井资料。
在该种可能的实现方式中,读取每条岩心压汞毛管压力曲线对应深度点的孔隙度频率谱数据,对其进行反向累积,计算得到反向累积饱和度,用以得到多条孔隙度频率谱反向累积曲线,孔隙度频率谱反向累积曲线反映反向累积饱和度(Sw)和孔隙度值分量的倒数1/POR的函数关系;根据102中的进行储层分类的地质分类标准,对多条孔隙度频率谱反向累积曲线进行储层分类,得到与102中储层分类对应的每类储层的平均孔隙度频率谱反向累积曲线。
场景二:岩心压汞毛管压力曲线对应深度不具有微电阻率成像测井资料。
在该种可能的实现方式中,因为测井资料的缺少,无法获取每一条岩心压汞毛管压力曲线对应深度点的孔隙度频率谱数据,则基于每一类储层分别估选孔隙度频率谱数据将其进行反向累积,进而求取每一类储层的平均孔隙度频率谱反向累积曲线。其中,估选孔隙度频率谱数据的方式可以是,根据岩心压汞毛管压力曲线对应的孔隙度,在孔隙度频率谱中选取体现的孔隙度与该岩心压汞毛管压力曲线的孔隙度相近的孔隙度频率谱数据。其中,孔隙度相近是指岩心压汞毛管压力曲线的孔隙度与孔隙度频率谱中的某一数据体现的孔隙度相同或者偏差程度小于一定阈值。该阈值可由技术人员预先设定或者根据实际环境情况进行调整。
104、建立每类储层的平均岩心压汞毛管压力曲线和平均孔隙度频率谱反向累积曲线之间的刻度函数。
依据102中的平均岩心压汞毛管压力曲线和103中的平均孔隙度频率谱反向累积曲线,采用分段非线性转换刻度方法,建立每类储层的平均岩心压汞毛管压力曲线和平均孔隙度频率谱反向累积曲线之间的刻度函数。该步骤包括但不限于1041和1042两个分步骤:
1041、根据饱和度相等的原则,针对每一类储层,在含水饱和度与进汞饱和度相同时,确定相同饱和度条件下的孔隙度值分量的倒数1/POR值和毛管压力Pc的值,建立1/POR和Pc之间的关系图像。
1042、根据1041中得到的图像显示的每一类储层毛管压力Pc和1/POR之间的转换关系特征,分别选取其拐点,区分每一类关系图像中的大孔段和小孔段,然后分段进行拟合,分别得到每一类储层中大小孔的毛管压力Pc和1/POR之间的刻度关系,即刻度函数。
105、结合孔隙度频率谱和每类储层的刻度函数分类构建深度连续的毛管压力曲线。
在本步骤中,最终毛管压力曲线的构建过程包括但不限于以下步骤1051~1053:
1051、依据前述地质分类标准,基于101中计算得到的孔隙度曲线,对101中计算得到的孔隙度频率谱进行储层分类,对每一类储层,应用104中已建立的该储层的毛管压力Pc和孔隙度值分量的倒数1/POR之间的刻度函数,求取毛管压力Pc的值,得到每类储层对应的多个毛管压力Pc的值;
1052、根据在101中计算得到的孔隙度频率谱数据,进行反向累积,求取得到反向累积饱和度Sw;在一种可能的实现方式中,该步骤中的Sw可以直接应用103中计算得到的反向累积饱和度。
1053、将求取的毛管压力Pc的值与反向累积饱和度Sw结合,即可构建得到深度连续的毛管压力曲线。
在本步骤中,由于步骤1051中基于深度连续的孔隙度频率谱进行储层分类和计算,且步骤1051得到的毛管压力Pc的值对应多类储层,所以将得到每类储层深度连续的毛管压力曲线。
在一种可能的实现方式中,将得到的每类储层深度连续的毛管压力曲线构建在同一图像中进行展示。
进一步的,根据构建得到的深度连续的毛管压力曲线可以计算得到最大孔喉半径和孔喉半径分布的连续图像。
在一种可能的实现方式中,展示曲线构建的结果图像,该结果图像包括微电阻率成像资料中的微电阻率成像、孔隙度频率谱、构建的深度连续的毛管压力曲线和基础储层信息。即最终将微电阻率成像测井资料、构建的深度连续的毛管压力曲线以及计算得到的最大孔喉半径和孔喉半径分布的连续图像等图像以及基础储层信息进行组合显示。基础储层信息包括从测井资料得到的井径、声波时差等信息,本申请实施例对此不作限定。
本申请实施例提供了一种基于微电阻成像测井资料的毛管压力曲线构建方法。通过实际测量的微电阻率成像测井资料计算得到孔隙度频率谱,再结合岩心压汞毛管压力曲线和储层的分类标准,建立平均岩心压汞毛管压力曲线和平均孔隙度频率谱反向累积曲线之间的刻度函数,建立起孔隙度频率谱与毛管压力曲线之间的关系,将孔隙度频率谱转换为毛管压力曲线,以实现微电阻率成像测井资料定量评价储层孔隙结构的目的,更精准地评价了储层的孔隙结构。
下面将结合示例性实施例介绍本申请提供的方法:
201、基于微电阻率成像测井资料计算得到孔隙度频率谱。
从专业的测井技术部门获取到实际测量的微电阻率成像测井资料,该实际测量的微电阻率成像测井资料包括以图像形式展示的微电阻率成像,直观展示储层中孔、缝、洞的分布情况。示例性地,该微电阻率成像中像素点的亮度反映出储层中对应深度、对应位置的电阻率,像素点呈暗色,则该点对应的电阻率较低,则实际储层中对应位置的孔、缝、洞发育情况较好,利于油气的存储,也就意味着该储层为有利储层的概率较高,开发价值较高。
基于实际测量获得的微电阻率成像测井资料的内容,利用图像中像素点对应的电阻率进行计算,得到孔隙度频率谱,还计算得到对应深度的孔隙度曲线。其中,孔隙度频率谱反映微电阻率成像对应深度的孔隙度的分布情况,孔隙度曲线反映对应深度的储层的孔隙度的大小;该步骤中计算时可以应用阿尔齐公式。孔隙度曲线则作为对孔隙度频率谱进行储层分类的依据,以及用于对储层孔隙度情况的基本展示。其中,孔隙度曲线和孔隙度频率谱中的曲线是深度连续的。
202、对岩心压汞毛管压力曲线进行储层分类。
在储层勘探过程中,不只获取数据、图像等测井资料,也会通过取心作业获取一定数量的岩心样本。对采集到的样本进行实验分析,获取岩心样本的储集空间标量参数。示例性地,获取每个岩心样本对应的实际深度、测量计算每个样本的孔隙度等。
在本申请实施例中,还对每个岩心样本进行压汞实验,得到每个岩心样本的岩心压汞毛管压力曲线,用以反映每个岩心样本的进汞饱和度(Shg)和毛管压力(Pc)之间的关系。每条岩心压汞毛管压力曲线不仅对应于一份岩心样本,而且对应于特定的深度和孔隙度。应用孔隙度将得到的岩心压汞毛管压力曲线分为三类储层。针对不同储层的地质情况,如地理位置、海拔、储层类型等,有多种依据岩心孔隙度的地质分类标准,示例性地,地质分类标准可以是孔隙度为2.0~6.0%时为Ⅲ类储层;孔隙度为6.0~12.0%时为Ⅱ类以及孔隙度>12.0%时的Ⅰ类储层。
在本申请中,针对示例性实施例中获取到的碳酸盐储层的多个岩心样本,测量对应的孔隙度,然后进行压汞实验,得到多条岩心压汞毛管压力曲线,并结合对应的地质分类标准(参见表1),对得到的岩心压汞毛管压力曲线按照地质分类标准加以分类,得到如2所示的岩心压汞毛管压力曲线图,图中相同颜色的曲线对应的岩心样本属于相同类型的储层。
表1地质分类标准
储层类型 | 分类标准(孔隙度,%) |
Ⅰ类 | >7.0% |
Ⅱ类 | 4.0~7.0% |
Ⅲ类 | 2.0~4.0% |
203、基于储层分类,求取每类储层的平均岩心压汞毛管压力曲线。
根据对岩心压汞毛管压力曲线分类后得到的图2,对每一类岩心压汞毛管压力曲线进行平均,分别得到Ⅰ类、Ⅱ类、Ⅲ类储层的平均岩心压汞毛管压力曲线,如图3所示,反映每一类储层的进汞饱和度(Shg)和毛管压力(Pc)之间的函数关系Shg=f(PC)。
204、基于孔隙度频率谱和岩心压汞毛管压力曲线,结合储层分类,求取前述每类储层的平均孔隙度频率谱反向累积曲线。
应用201中得到的孔隙度频率谱和202中得到的岩心压汞毛管压力曲线,结合202中的储层分类情况,求取与步骤202中储层对应的每类储层的平均孔隙度频率谱反向累积曲线。该曲线反映每一类储层的反向累积饱和度(Sw)和孔隙度频率谱中各孔隙度值分量的倒数1/POR值之间的函数关系,该函数关系表示为Sw=f(1/POR)。其中,反向累积饱和度(Sw)可以表示储层中的含水饱和度。
在一种可能的实现方式中,岩心压汞毛管压力曲线图中各曲线的对应深度具有微电阻率成像测井资料,则从201中得到的孔隙度频率谱中读取对应深度点的孔隙度频率谱数据,对获取的孔隙度频率谱数据进行反向累积,得到与孔隙度频率谱数据对应的反向累积饱和度,用以得到孔隙度频率谱反向累积曲线;并根据202中的储层分类标准,对得到孔隙度频率谱反向累积曲线进行分类;对每类孔隙度频率谱反向累积曲线求平均,得到平均孔隙度频率谱反向累积曲线。
在一种可能的实现方式中,岩心压汞毛管压力曲线图中各曲线对应深度未进行微电阻率成像测井时,针对每一类储层分别估选孔隙度频率谱数据,将其进行反向累积,得到孔隙度频率谱反向累积曲线,进而求取每类储层的平均孔隙度频率谱反向累积曲线。其中,估选孔隙度频率谱数据的一种可选方式为:根据岩心压汞毛管压力曲线图中各曲线对应的孔隙度,在孔隙度频率谱中选取对应的孔隙度与各曲线对应的孔隙度相近的孔隙度频率谱数据。示例性地,当岩心压汞毛管压力曲线对应的孔隙度的大小和孔隙度频率谱中某一孔隙度频率谱数据体现的孔隙度大小相同或者偏差在一定阈值内,则选取该孔隙度频率谱数据。其中,阈值的数值可以有技术人员预先设定,也可以根据储层的地质情况进行调整。
在本申请示例性实施例中,岩心压汞毛管压力曲线图中各曲线对应深度具有微电阻率成像测井资料,采用相应的计算方法,得到了如图4所示的平均孔隙度频率谱反向累积曲线。
205、建立每类储层的平均岩心压汞毛管压力曲线和平均孔隙度频率谱反向累积曲线之间的刻度函数。
本步骤中,采用分段非线性转换刻度方法,建立每类储层的岩心压汞平均毛管压力曲线和平均孔隙度频率谱反向累积曲线之间的刻度函数。具体步骤包括但不限于如下2051~2052:
2051、依据204中获得的平均孔隙度频率谱反向累积曲线和203中获得的平均岩心压汞毛管压力曲线,根据饱和度相等的原则,针对每一类储层,在含水饱和度与进汞饱和度相同时,确定相同饱和度条件下的孔隙度值分量的倒数1/POR的值和毛管压力Pc的值,建立1/POR和Pc之间的关系图像。
2052、根据2051中得到的图像显示的每一类储层毛管压力Pc的值和孔隙度值分量的倒数1/POR的值之间的关系图像,分别选取其拐点,然后分段进行拟合,得到每一类储层中Pc和1/POR之间的刻度关系,即刻度函数。
在大孔段,该刻度函数的数学形式为:
式中:Pc表示压汞毛管压力,单位为MPa;1/POR表示孔隙度频率谱中各孔隙度值分量的倒数,单位1/%;m1为待定参数,由非线性拟合得到;n1为待定参数,由非线性拟合得到;
在小孔段,该刻度函数的数学形式为:
式中:Pc表示压汞毛管压力,单位为MPa;1/POR表示孔隙度频率谱中各孔隙度值分量的倒数,单位1/%;m2为待定参数,由非线性拟合得到;n2为待定参数,由非线性拟合得到;
其中,大孔段和小孔段的划分依据人员经验和1/POR和Pc之间的关系图像表现出的关系特征,如拐点对应的孔喉半径,进行确定。
在本申请示例性实施例中,针对缝洞型碳酸盐储层获得的1/POR和Pc之间的关系图像如图5所示。由图5可知,每类储层对应的1/POR和Pc的曲线具有明显的拐点位置。根据图像信息进行分段拟合,得到图5中显示的y和x之间的函数关系,即每类储层中的大孔段和小孔段的刻度函数。其中,y代表毛管压力Pc,x代表孔隙度频率谱中各孔隙度值分量的倒数1/POR,R为相关系数,R2越接近1,刻度函数的精准度越高。
206、结合孔隙度频率谱和每类储层的刻度函数分类构建深度连续的毛管压力曲线。
该步骤可以依据以下2061~2063步骤进行。
2061、结合前述的地质分类标准,依据201中计算得到的孔隙度曲线,对201中计算得到的孔隙度频率谱进行储层分类。对分类得到的每一类储层,应用205中已建立的毛管压力Pc和1/POR之间的刻度函数,求取毛管压力Pc的值;205中针对每一类储层建立了刻度函数,所以,在对微电阻率成像资料进行计算的过程中,需要根据孔隙度频率谱、依据上述内容中的分类标准进行储层分类;然后应用相应的刻度函数进行计算,得到每一类储层对应的多个毛管压力Pc。
2062、根据在201中计算得到的孔隙度频率谱数据,进行反向累积,求取得到反向累积饱和度Sw;该步骤中的Sw可以直接应用204中计算得到的反向累积饱和度。
2063、将步骤2061中求取的每类储层对应的多个毛管压力Pc,与相对应的反向累积饱和度Sw结合,即可构建得到每一类储层的毛管压力曲线。由于毛管压力Pc和反向累积饱和度Sw是基于深度连续的孔隙度频率谱得到的,构建的每一类储层的毛管压力曲线是深度连续的。
在示例性实施例中,依据地质分类标准将储层分为三类,得到图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”)或“响应于确定”或“响应于检测到”。类似地,根据上下文,短语“若确定...”或“若检测到[所陈述的条件或事件]”可被解释为意指“在确定...时”或“响应于确定...”或“在检测到[所陈述的条件或事件]时”或“响应于检测到[所陈述的条件或事件]”。
还应理解,说明书通篇中提到的“一个实施例”、“一实施例”、“一种可能的实现方式”意味着与实施例或实现方式有关的特定特征、结构或特性包括在本申请的至少一个实施例中。因此,在整个说明书各处出现的“在一个实施例中”或“在一实施例中”、“一种可能的实现方式”未必一定指相同的实施例。此外,这些特定的特征、结构或特性可以任意适合的方式结合在一个或多个实施例中。
以上所述仅为本申请的可选实施例,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本申请的保护范围之内。
Claims (10)
1.一种基于测井资料的曲线构建方法,其特征在于,所述方法包括:
基于微电阻率成像测井资料计算孔隙度频率谱;
对岩心压汞毛管压力曲线进行储层分类,基于分类后的岩心压汞毛管压力曲线,求取每类储层的平均岩心压汞毛管压力曲线;
基于所述孔隙度频率谱和所述岩心压汞毛管压力曲线,结合所述储层分类,求取所述每类储层的平均孔隙度频率谱反向累积曲线;
建立所述每类储层的平均岩心压汞毛管压力曲线和平均孔隙度频率谱反向累积曲线之间的刻度函数;
结合所述孔隙度频率谱和所述每类储层的刻度函数分类构建深度连续的毛管压力曲线。
2.根据权利要求1所述方法,其特征在于,所述对岩心压汞毛管压力曲线进行储层分类,包括:
结合地质分类标准,将所述岩心压汞毛管压力曲线根据对应的岩心孔隙度进行分类。
3.根据权利要求1所述方法,其特征在于,所述基于所述孔隙度频率谱和所述岩心压汞毛管压力曲线,结合所述储层分类,求取所述每类储层的平均孔隙度频率谱反向累积曲线,包括:
从所述孔隙度频率谱中读取所述岩心压汞毛管压力曲线对应深度的孔隙度频率谱数据;
对读取的所述孔隙度频率谱数据进行反向累积,结合所述储层分类,求取各类储层的平均孔隙度频率谱反向累积曲线。
4.根据权利要求1所述方法,其特征在于,所述基于所述孔隙度频率谱和所述岩心压汞毛管压力曲线,结合所述储层分类,求取所述每类储层的平均孔隙度频率谱反向累积曲线,包括:
基于所述每类储层,从所述孔隙度频率谱中估选孔隙度频率谱数据;
对估选的所述孔隙度频率谱数据进行反向累积,结合所述储层分类,求取所述各类储层的平均孔隙度频率谱反向累积曲线。
5.根据权利要求1所述方法,其特征在于,所述建立所述每类储层的平均岩心压汞毛管压力曲线和平均孔隙度频率谱反向累积曲线之间的刻度函数,包括:
采用分段非线性转换刻度方法,建立所述每类储层的平均岩心压汞毛管压力曲线和平均孔隙度频率谱反向累积曲线之间的刻度函数。
6.根据权利要求1所述方法,其特征在于,所述结合所述孔隙度频率谱和所述每类储层的刻度函数分类构建深度连续的毛管压力曲线之后,还包括:
根据构建的深度连续的毛管压力曲线,计算得到最大孔喉半径图像和孔喉半径分布图像。
7.根据权利要求1-6任一所述方法,其特征在于,所述方法还包括:
展示曲线构建的结果图像,所述结果图像包括所述微电阻率成像资料中的微电阻率成像、所述孔隙度频率谱、所述构建的深度连续的毛管压力曲线和基础储层信息。
8.一种基于测井资料的曲线构建装置,其特征在于,所述装置包括:
第一计算模块,用于基于微电阻率成像测井资料计算孔隙度频率谱;
第一求取模块,用于对岩心压汞毛管压力曲线进行储层分类,基于分类后的岩心压汞毛管压力曲线,求取每类储层的平均岩心压汞毛管压力曲线;
第二求取模块,用于基于所述孔隙度频率谱和所述岩心压汞毛管压力曲线,结合所述储层分类,求取所述每类储层的平均孔隙度频率谱反向累积曲线;
第一构建模块,用于建立所述每类储层的平均岩心压汞毛管压力曲线和平均孔隙度频率谱反向累积曲线之间的刻度函数;
第二构建模块,结合所述孔隙度频率谱和所述每类储层的刻度函数分类构建深度连续的毛管压力曲线。
9.根据权利要求8所述的装置,其特征在于,所述装置还包括:
第二计算模块,用于根据构建的深度连续的毛管压力曲线,计算得到最大孔喉半径图像和孔喉半径分布图像。
10.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质中存储有至少一条程序代码,所述程序代码由处理器加载并执行以实现如权利要求1至权利要求7中任一项所述的基于测井资料的曲线构建方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110168757.1A CN112983394A (zh) | 2021-02-07 | 2021-02-07 | 基于测井资料的曲线构建方法、装置和存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110168757.1A CN112983394A (zh) | 2021-02-07 | 2021-02-07 | 基于测井资料的曲线构建方法、装置和存储介质 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112983394A true CN112983394A (zh) | 2021-06-18 |
Family
ID=76348927
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110168757.1A Pending CN112983394A (zh) | 2021-02-07 | 2021-02-07 | 基于测井资料的曲线构建方法、装置和存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112983394A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113818867A (zh) * | 2021-09-17 | 2021-12-21 | 成都理工大学 | 一种构建伪毛管压力曲线方法、系统、介质、设备及应用 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NO930575D0 (no) * | 1992-03-16 | 1993-02-18 | Texaco Development Corp | Fremgangsmaate for aa vurdere permeabiliteten til formasjoner ut fra logger som anvender nervenettverk |
CN103592690A (zh) * | 2013-10-24 | 2014-02-19 | 长江大学 | 基于电成像测井孔隙度谱信息自动识别储层裂缝的方法 |
CN105386753A (zh) * | 2015-10-28 | 2016-03-09 | 中国地质大学(北京) | 一种利用核磁共振测井构造伪毛管压力曲线方法 |
CN110618158A (zh) * | 2019-10-28 | 2019-12-27 | 中石化石油工程技术服务有限公司 | 利用核磁共振信息构建岩心毛管压力曲线的方法 |
CN112098291A (zh) * | 2019-06-18 | 2020-12-18 | 中国石油化工股份有限公司 | 基于电成像孔隙度谱的渗透率计算方法及系统 |
-
2021
- 2021-02-07 CN CN202110168757.1A patent/CN112983394A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NO930575D0 (no) * | 1992-03-16 | 1993-02-18 | Texaco Development Corp | Fremgangsmaate for aa vurdere permeabiliteten til formasjoner ut fra logger som anvender nervenettverk |
CN103592690A (zh) * | 2013-10-24 | 2014-02-19 | 长江大学 | 基于电成像测井孔隙度谱信息自动识别储层裂缝的方法 |
CN105386753A (zh) * | 2015-10-28 | 2016-03-09 | 中国地质大学(北京) | 一种利用核磁共振测井构造伪毛管压力曲线方法 |
CN112098291A (zh) * | 2019-06-18 | 2020-12-18 | 中国石油化工股份有限公司 | 基于电成像孔隙度谱的渗透率计算方法及系统 |
CN110618158A (zh) * | 2019-10-28 | 2019-12-27 | 中石化石油工程技术服务有限公司 | 利用核磁共振信息构建岩心毛管压力曲线的方法 |
Non-Patent Citations (1)
Title |
---|
肖亮: "利用井眼电成像测井评价裂缝性致密砂岩储层孔隙结构的方法", AAPG BULLETIN, vol. 104, no. 1, pages 205 - 226 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113818867A (zh) * | 2021-09-17 | 2021-12-21 | 成都理工大学 | 一种构建伪毛管压力曲线方法、系统、介质、设备及应用 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20090260883A1 (en) | Continuous measurement of heterogeneity of geomaterials | |
Sanderson et al. | Line sampling of fracture swarms and corridors | |
US8290711B2 (en) | Method of estimating the fracture density in a rock medium | |
CA2149647C (fr) | Methode d'analyse des traces sismiques utilisant une technique de calibrage statistique pour en deduire des proprietes geologiques | |
US20110292764A1 (en) | Bayesian DHI for Seismic Data | |
CN103026202A (zh) | 获取多孔介质的相容和综合物理性质的方法 | |
US9977996B2 (en) | Characterizing porosity distribution from a borehole image | |
US20150234069A1 (en) | System and Method for Quantifying Vug Porosity | |
CN110320139A (zh) | 缝洞型储层裂缝孔隙度定量评价方法及系统 | |
US10928548B2 (en) | Rock type based free water level inversion | |
US20180342073A1 (en) | Methods and systems for automatically analyzing an image representative of a formation | |
US20140297186A1 (en) | Rock Classification Based on Texture and Composition | |
AU2014290779B2 (en) | System and method for estimating porosity distribution in subterranean reservoirs | |
CN111144225B (zh) | 一种地质勘探岩芯的识别方法、装置和存储介质 | |
CN106777649B (zh) | 一种裂缝型储集层孔隙结构定量评价方法 | |
Sastre Jurado et al. | Probabilistic 3D modelling of shallow soil spatial variability using dynamic cone penetrometer results and a geostatistical method | |
CN112983394A (zh) | 基于测井资料的曲线构建方法、装置和存储介质 | |
CN103899305A (zh) | 一种碳酸盐岩储层流体性质的识别方法及设备 | |
CN108035710B (zh) | 基于数据挖掘划分深层岩石地质相的方法 | |
CN114060015B (zh) | 一种致密砂岩含气性的评价方法及装置 | |
Berre et al. | Experience on controlled-source electromagnetic performance for exploration in Norway | |
CN111963162A (zh) | 基于岩性与水性的双参数流体性质识别方法 | |
CN115749755A (zh) | 岩屑地层流体性质评价方法、装置、设备及存储介质 | |
CN114428365A (zh) | 一种基于含气特征的致密砂岩气层测井识别方法及装置 | |
US11905809B2 (en) | Determining reservoir heterogeneity for optimized drilling location |
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 |