储层大孔径孔隙度计算方法、装置、电子设备及介质与流程

1.本发明涉及石油地球物理勘探技术领域，更具体地，涉及一种储层大孔径孔隙度计算方法、装置、电子设备及介质。


背景技术：

2.随着致密砂岩等孔隙结构较复杂的储层成为各油公司增储上产的主要领域，孔隙结构评价成为勘探开发技术中的一项重要内容。利用测井资料开展孔隙结构评价，前人已完成了一些工作，例如：
3.将储层孔隙简化为一根弯曲毛细管，推导了孔喉平均半径等计算公式，建立自然伽马曲线与公式参数的拟合关系，实现孔隙结构参数的计算；
4.用含油饱和度消除储层含烃的影响，然后横向驰豫时间计算孔隙结构参数，实现储层孔隙结构连续表征；
5.将储层孔隙简化为一组毛细管，推导孔隙结构参数与电阻率参数的关系，从而形成了由电阻率资料计算孔隙结构参数的方法；
6.通过沉积相、成岩相和裂缝相等特征，总结岩石物理相与储层孔隙结构关系，结合测井资料划分岩石物理相来实现孔隙结构的分类评价；
7.根据压汞曲线对储层孔隙结构分类，提取测井参数利用自组织神经网络技术实现孔隙结构分类判别。
8.综合分析前人的工作，要么将储层孔隙简化为毛细管模型，然后推导孔隙结构参数计算公式，由于储层孔隙模型过于理想，与实际地层情况差异大，适用范围有限；要么将储层孔隙结构依据成岩相、毛管压力曲线特征等信息分类，将这种分类与一些测井参数建立关系，从而利用图版或神经网络等方法实现孔隙结构的定性评价；要么利用核磁共振测井，核磁共振测井能连续定量评价储层孔隙结构，但属于特殊测井，资料采集成本较高，大部分井无核磁资料，限制了该方法的使用。
9.因此，有必要开发一种储层大孔径孔隙度计算方法、装置、电子设备及介质。
10.公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。


技术实现要素：

11.本发明提出了一种储层大孔径孔隙度计算方法、装置、电子设备及介质，其能够通过渗透率与大孔径孔隙度良好的相关性，由渗透率计算大孔径孔隙度，进而与测井参数建立计算关系，简单实用，可操作性强，在实际应用中效果较好。
12.第一方面，本公开实施例提供了一种储层大孔径孔隙度计算方法，包括：
13.针对每块岩心样品计算样品大孔径孔隙度；
14.建立所述样品大孔径孔隙度与所述渗透率的第一拟合关系；
15.根据所述第一拟合关系与目标区域的渗透率，计算所述目标区域的大孔径孔隙度；
16.计算孔隙度差，建立所述目标区域的大孔径孔隙度与所述孔隙度差的第二拟合关系；
17.根据目标井的孔隙度差与所述第二拟合关系，计算所述目标井的大孔径孔隙度。
18.优选地，针对每块岩心样品计算样品大孔径孔隙度包括：
19.计算毛管力，确定所述毛管力对应的进汞饱和度；
20.通过所述进汞饱和度与岩心孔隙度，计算所述样品大孔径孔隙度。
21.优选地，通过公式(1)计算毛管力：
[0022][0023]
通过公式(2)计算所述样品大孔径孔隙度：
[0024][0025]
其中，pc为毛管力，rc为计算孔径，sc为毛管力为pc对应的进汞饱和度，为岩心孔隙度。
[0026]
优选地，计算孔隙度差包括：
[0027]
在测井曲线上提取测井声波时差值与密度值；
[0028]
计算声波时差孔隙度与密度孔隙度；
[0029]
通过所述声波时差孔隙度与所述密度孔隙度，计算所述孔隙度差。
[0030]
优选地，通过公式(3)计算声波时差孔隙度：
[0031][0032]
其中，pora为声波时差孔隙度，tm为岩石骨架声波时差，ac为声波时差值。
[0033]
优选地，通过公式(4)计算密度孔隙度：
[0034][0035]
其中，pord为密度孔隙度，dg为岩石密度骨架，df为地层水密度，den为密度值。
[0036]
优选地，通过公式(5)计算所述孔隙度差：
[0037]
δpor＝pord-pora
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)
[0038]
其中，
△
por为孔隙度差。
[0039]
作为本公开实施例的一种具体实现方式，
[0040]
第二方面，本公开实施例还提供了一种储层大孔径孔隙度计算装置，包括：
[0041]
样品计算模块，针对每块岩心样品计算样品大孔径孔隙度；
[0042]
第一拟合模块，建立所述样品大孔径孔隙度与所述渗透率的第一拟合关系；
[0043]
目标区域计算模块，根据所述第一拟合关系与目标区域的渗透率，计算所述目标区域的大孔径孔隙度；
[0044]
第二拟合模块，计算孔隙度差，建立所述目标区域的大孔径孔隙度与所述孔隙度差的第二拟合关系；
[0045]
目标井计算模块，根据目标井的孔隙度差与所述第二拟合关系，计算所述目标井的大孔径孔隙度。
[0046]
优选地，针对每块岩心样品计算样品大孔径孔隙度包括：
[0047]
计算毛管力，确定所述毛管力对应的进汞饱和度；
[0048]
通过所述进汞饱和度与岩心孔隙度，计算所述样品大孔径孔隙度。
[0049]
优选地，通过公式(1)计算毛管力：
[0050][0051]
通过公式(2)计算所述样品大孔径孔隙度：
[0052][0053]
其中，pc为毛管力，rc为计算孔径，sc为毛管力为pc对应的进汞饱和度，为岩心孔隙度。
[0054]
优选地，计算孔隙度差包括：
[0055]
在测井曲线上提取测井声波时差值与密度值；
[0056]
计算声波时差孔隙度与密度孔隙度；
[0057]
通过所述声波时差孔隙度与所述密度孔隙度，计算所述孔隙度差。
[0058]
优选地，通过公式(3)计算声波时差孔隙度：
[0059][0060]
其中，pora为声波时差孔隙度，tm为岩石骨架声波时差，ac为声波时差值。
[0061]
优选地，通过公式(4)计算密度孔隙度：
[0062][0063]
其中，pord为密度孔隙度，dg为岩石密度骨架，df为地层水密度，den为密度值。
[0064]
优选地，通过公式(5)计算所述孔隙度差：
[0065]
δpor＝pord-pora
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)
[0066]
其中，
△
por为孔隙度差。
[0067]
第三方面，本公开实施例还提供了一种电子设备，该电子设备包括：
[0068]
存储器，存储有可执行指令；
[0069]
处理器，所述处理器运行所述存储器中的所述可执行指令，以实现所述的储层大孔径孔隙度计算方法。
[0070]
第四方面，本公开实施例还提供了一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现所述的储层大孔径孔隙度计算方法。
[0071]
本发明的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
[0072]
通过结合附图对本发明示例性实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施例中，相同的参考标号通常代表相同部件。
[0073]
图1示出了根据本发明的一个实施例的储层大孔径孔隙度计算方法的步骤的流程图。
[0074]
图2示出了根据本发明的一个实施例的储层孔径分布的示意图。
[0075]
图3示出了根据本发明的一个实施例的大孔径孔隙度与渗透率的拟合关系的示意图。
[0076]
图4示出了根据本发明的一个实施例的大孔径孔隙度与
△
por的拟合关系的示意图。
[0077]
图5示出了根据本发明的一个实施例的xx26井测井计算大孔径孔隙度与压汞实验得到的大孔径孔隙度的对比示意图。
[0078]
图6示出了根据本发明的一个实施例的一种储层大孔径孔隙度计算装置的框图。
[0079]
附图标记说明：
[0080]
201、样品计算模块；202、第一拟合模块；203、目标区域计算模块；204、第二拟合模块；205、目标井计算模块。
具体实施方式
[0081]
下面将更详细地描述本发明的优选实施方式。虽然以下描述了本发明的优选实施方式，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。
[0082]
本发明提供一种储层大孔径孔隙度计算方法，包括：
[0083]
针对每块岩心样品计算样品大孔径孔隙度；在一个示例中，针对每块岩心样品计算样品大孔径孔隙度包括：计算毛管力，确定毛管力对应的进汞饱和度；通过进汞饱和度与岩心孔隙度，计算样品大孔径孔隙度。
[0084]
在一个示例中，通过公式(1)计算毛管力：
[0085][0086]
通过公式(2)计算样品大孔径孔隙度：
[0087][0088]
其中，pc为毛管力，rc为计算孔径，sc为毛管力为pc对应的进汞饱和度，为岩心孔隙度。
[0089]
具体地，分析目标油气田压汞实验数据，获取岩心样品孔喉分布情况，确定大孔径孔隙度对应的孔径下限。一般而言，致密砂岩等复杂储层中来自优质储层的样品孔径分布一般呈多峰分布，取最大孔径峰和次大孔径峰之间波谷对应的孔径为大孔径孔隙下限，记为rc。通过公式(1)计算毛管力，确定毛管力对应的进汞饱和度；通过进汞饱和度与岩心孔隙度，通过公式(2)计算样品大孔径孔隙度。
[0090]
建立样品大孔径孔隙度与渗透率的第一拟合关系；由压汞实验得到的大孔径孔隙
度与物性实验得到的渗透率k建立第一拟合关系为：
[0091][0092]
其中，a、b为拟合系数。
[0093]
根据第一拟合关系与目标区域的渗透率，通过公式(6)计算目标区域的大孔径孔隙度。
[0094]
计算孔隙度差，建立目标区域的大孔径孔隙度与孔隙度差的第二拟合关系；在一个示例中，计算孔隙度差包括：在测井曲线上提取测井声波时差值与密度值；计算声波时差孔隙度与密度孔隙度；通过声波时差孔隙度与密度孔隙度，计算孔隙度差。
[0095]
在一个示例中，通过公式(3)计算声波时差孔隙度：
[0096][0097]
其中，pora为声波时差孔隙度，tm为岩石骨架声波时差，ac为声波时差值。
[0098]
在一个示例中，通过公式(4)计算密度孔隙度：
[0099][0100]
其中，pord为密度孔隙度，dg为岩石密度骨架，df为地层水密度，den为密度值。
[0101]
在一个示例中，通过公式(5)计算孔隙度差：
[0102]
δpor＝pord-pora
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)
[0103]
其中，
△
por为孔隙度差。
[0104]
具体地，根据目标区域的大孔径孔隙度，在对应深度上提取测井声波时差值ac和密度值den，通过公式(3)计算声波时差孔隙度pora，通过公式(4)计算密度孔隙度pord，进而通过公式(5)计算孔隙度差
△
por，由孔隙度差与目标区域的大孔径孔隙度进行拟合，获得第二拟合关系为：
[0105][0106]
其中，c、d为拟合系数。
[0107]
根据公式(4)-(6)计算目标井的孔隙度差，将目标井的孔隙度差代入第二拟合关系，计算获得目标井的大孔径孔隙度。
[0108]
本发明还提供一种储层大孔径孔隙度计算装置，包括：。
[0109]
样品计算模块，针对每块岩心样品计算样品大孔径孔隙度；在一个示例中，针对每块岩心样品计算样品大孔径孔隙度包括：计算毛管力，确定毛管力对应的进汞饱和度；通过进汞饱和度与岩心孔隙度，计算样品大孔径孔隙度。
[0110]
在一个示例中，通过公式(1)计算毛管力：
[0111][0112]
通过公式(2)计算样品大孔径孔隙度：
[0113][0114]
其中，pc为毛管力，rc为计算孔径，sc为毛管力为pc对应的进汞饱和度，为岩心孔
隙度。
[0115]
具体地，分析目标油气田压汞实验数据，获取岩心样品孔喉分布情况，确定大孔径孔隙度对应的孔径下限。一般而言，致密砂岩等复杂储层中来自优质储层的样品孔径分布一般呈多峰分布，取最大孔径峰和次大孔径峰之间波谷对应的孔径为大孔径孔隙下限，记为rc。通过公式(1)计算毛管力，确定毛管力对应的进汞饱和度；通过进汞饱和度与岩心孔隙度，通过公式(2)计算样品大孔径孔隙度。
[0116]
第一拟合模块，建立样品大孔径孔隙度与渗透率的第一拟合关系；由压汞实验得到的大孔径孔隙度与物性实验得到的渗透率k建立第一拟合关系为公式(6)。
[0117]
目标区域计算模块，根据第一拟合关系与目标区域的渗透率，通过公式(6)计算目标区域的大孔径孔隙度。
[0118]
第二拟合模块，计算孔隙度差，建立目标区域的大孔径孔隙度与孔隙度差的第二拟合关系；在一个示例中，计算孔隙度差包括：在测井曲线上提取测井声波时差值与密度值；计算声波时差孔隙度与密度孔隙度；通过声波时差孔隙度与密度孔隙度，计算孔隙度差。
[0119]
在一个示例中，通过公式(3)计算声波时差孔隙度：
[0120][0121]
其中，pora为声波时差孔隙度，tm为岩石骨架声波时差，ac为声波时差值。
[0122]
在一个示例中，通过公式(4)计算密度孔隙度：
[0123][0124]
其中，pord为密度孔隙度，dg为岩石密度骨架，df为地层水密度，den为密度值。
[0125]
在一个示例中，通过公式(5)计算孔隙度差：
[0126]
δpor＝pord-pora
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)
[0127]
其中，
△
por为孔隙度差。
[0128]
具体地，根据目标区域的大孔径孔隙度，在对应深度上提取测井声波时差值ac和密度值den，通过公式(3)计算声波时差孔隙度pora，通过公式(4)计算密度孔隙度pord，进而通过公式(5)计算孔隙度差
△
por，由孔隙度差与目标区域的大孔径孔隙度进行拟合，获得第二拟合关系为公式(7)。
[0129]
目标井计算模块，根据公式(3)-(5)计算目标井的孔隙度差，将目标井的孔隙度差代入第二拟合关系，计算获得目标井的大孔径孔隙度。
[0130]
本发明还提供一种电子设备，电子设备包括：存储器，存储有可执行指令；处理器，处理器运行存储器中的可执行指令，以实现上述的储层大孔径孔隙度计算方法。
[0131]
本发明还提供一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现上述的储层大孔径孔隙度计算方法。
[0132]
为便于理解本发明实施例的方案及其效果，以下给出四个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
[0133]
实施例1
[0134]
图1示出了根据本发明的一个实施例的储层大孔径孔隙度计算方法的步骤的流程图。
[0135]
如图1所示，该储层大孔径孔隙度计算方法包括：步骤101，针对每块岩心样品计算样品大孔径孔隙度；步骤102，建立样品大孔径孔隙度与渗透率的第一拟合关系；步骤103，根据第一拟合关系与目标区域的渗透率，计算目标区域的大孔径孔隙度；步骤104，计算孔隙度差，建立目标区域的大孔径孔隙度与孔隙度差的第二拟合关系；步骤105，根据目标井的孔隙度差与第二拟合关系，计算目标井的大孔径孔隙度。
[0136]
图2示出了根据本发明的一个实施例的储层孔径分布的示意图。
[0137]
对某井区利用常规测井资料获取大孔径孔隙度，如图2所示，由来自优质储层的样品压汞实验得到的储层孔径分布中可知，样品孔隙呈三峰分布，从右至左第一个峰与第二个峰之间的波谷对应的孔隙约为0.1μm，因此确定大孔径孔隙的孔径下限为0.1μm。
[0138]
以hh1057-3岩心深度为2230.7m的样品为例计算大孔径孔隙度，大孔径孔隙的孔径下限为0.1μm，由公式(1)得到对应的毛管力为7.35mpa，由该样品的压汞实验数据可知，7.35mpa时的进汞饱和度为28.73％，该样品岩心物性实验得到孔隙度为7.5％，由公式(2)得到该样品大孔径孔隙度为2.15％，由此可得到目标区开展压汞实验的19个岩心样品的大孔径孔隙度。
[0139]
图3示出了根据本发明的一个实施例的大孔径孔隙度与渗透率的拟合关系的示意图。
[0140]
将压汞实验样品得到的大孔径孔隙度与物性实验的渗透率建立拟合关系，如图3所示，可知公式(6)中a＝9.8922，b＝0.5729。
[0141]
由公式(6)计算目标区184块开展了物性实验(渗透率测量)的岩心样品的大孔径孔隙度。
[0142]
图4示出了根据本发明的一个实施例的大孔径孔隙度与
△
por的拟合关系的示意图。
[0143]
将184块进行了物性实验的岩心样品，在对应深度上提取声波时差测井值、密度测井值，由公式(3)-(5)计算密度孔隙度与声波孔隙度差值
△
por，将由184块样品由渗透率计算得到的大孔径孔隙度与
△
por建立拟合关系，如图4所示，可知公式(7)c＝53.185，d＝2.8516。
[0144]
图5示出了根据本发明的一个实施例的xx26井测井计算大孔径孔隙度与压汞实验得到的大孔径孔隙度的对比示意图。
[0145]
对需要进行孔隙结构的井(以xx26井为例)按公式(3)-(5)计算密度孔隙度与声波孔隙度差值
△
por，按公式(7)计算大孔径孔隙度，如图5第五道，phix为测井计算的大孔径孔隙度，cphix(杆状图显示)为压汞实验得到的大孔径孔隙度，两者吻合较好，显示了发明在研究区有较好的应用效果。
[0146]
实施例2
[0147]
图3示出了根据本发明的一个实施例的一种储层大孔径孔隙度计算装置的框图。
[0148]
如图3所示，该储层大孔径孔隙度计算装置，包括：
[0149]
样品计算模块201，针对每块岩心样品计算样品大孔径孔隙度；
[0150]
第一拟合模块202，建立样品大孔径孔隙度与渗透率的第一拟合关系；
[0151]
目标区域计算模块203，根据第一拟合关系与目标区域的渗透率，计算目标区域的大孔径孔隙度；
[0152]
第二拟合模块204，计算孔隙度差，建立目标区域的大孔径孔隙度与孔隙度差的第二拟合关系；
[0153]
目标井计算模块205，根据目标井的孔隙度差与第二拟合关系，计算目标井的大孔径孔隙度。
[0154]
作为可选方案，针对每块岩心样品计算样品大孔径孔隙度包括：
[0155]
计算毛管力，确定毛管力对应的进汞饱和度；
[0156]
通过进汞饱和度与岩心孔隙度，计算样品大孔径孔隙度。
[0157]
作为可选方案，通过公式(1)计算毛管力：
[0158][0159]
通过公式(2)计算样品大孔径孔隙度：
[0160][0161]
其中，pc为毛管力，rc为计算孔径，sc为毛管力为pc对应的进汞饱和度，为岩心孔隙度。
[0162]
作为可选方案，计算孔隙度差包括：
[0163]
在测井曲线上提取测井声波时差值与密度值；
[0164]
计算声波时差孔隙度与密度孔隙度；
[0165]
通过声波时差孔隙度与密度孔隙度，计算孔隙度差。
[0166]
作为可选方案，通过公式(3)计算声波时差孔隙度：
[0167][0168]
其中，pora为声波时差孔隙度，tm为岩石骨架声波时差，ac为声波时差值。
[0169]
作为可选方案，通过公式(4)计算密度孔隙度：
[0170][0171]
其中，pord为密度孔隙度，dg为岩石密度骨架，df为地层水密度，den为密度值。
[0172]
作为可选方案，通过公式(5)计算孔隙度差：
[0173]
δpor＝pord-pora
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)
[0174]
其中，
△
por为孔隙度差。
[0175]
实施例3
[0176]
本公开提供一种电子设备包括，该电子设备包括：存储器，存储有可执行指令；处理器，处理器运行存储器中的可执行指令，以实现上述储层大孔径孔隙度计算方法。
[0177]
根据本公开实施例的电子设备包括存储器和处理器。
[0178]
该存储器用于存储非暂时性计算机可读指令。具体地，存储器可以包括一个或多个计算机程序产品，该计算机程序产品可以包括各种形式的计算机可读存储介质，例如易失性存储器和/或非易失性存储器。该易失性存储器例如可以包括随机存取存储器(ram)
和/或高速缓冲存储器(cache)等。该非易失性存储器例如可以包括只读存储器(rom)、硬盘、闪存等。
[0179]
该处理器可以是中央处理单元(cpu)或者具有数据处理能力和/或指令执行能力的其它形式的处理单元，并且可以控制电子设备中的其它组件以执行期望的功能。在本公开的一个实施例中，该处理器用于运行该存储器中存储的该计算机可读指令。
[0180]
本领域技术人员应能理解，为了解决如何获得良好用户体验效果的技术问题，本实施例中也可以包括诸如通信总线、接口等公知的结构，这些公知的结构也应包含在本公开的保护范围之内。
[0181]
有关本实施例的详细说明可以参考前述各实施例中的相应说明，在此不再赘述。
[0182]
实施例4
[0183]
本公开实施例提供一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现所述的储层大孔径孔隙度计算方法。
[0184]
根据本公开实施例的计算机可读存储介质，其上存储有非暂时性计算机可读指令。当该非暂时性计算机可读指令由处理器运行时，执行前述的本公开各实施例方法的全部或部分步骤。
[0185]
上述计算机可读存储介质包括但不限于：光存储介质(例如：cd－rom和dvd)、磁光存储介质(例如：mo)、磁存储介质(例如：磁带或移动硬盘)、具有内置的可重写非易失性存储器的媒体(例如：存储卡)和具有内置rom的媒体(例如：rom盒)。
[0186]
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
[0187]
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。