CN117074274B - 一种基于全直径岩心样品的孔隙度计算方法、装置及介质 - Google Patents

一种基于全直径岩心样品的孔隙度计算方法、装置及介质 Download PDF

Info

Publication number
CN117074274B
CN117074274B CN202311104089.1A CN202311104089A CN117074274B CN 117074274 B CN117074274 B CN 117074274B CN 202311104089 A CN202311104089 A CN 202311104089A CN 117074274 B CN117074274 B CN 117074274B
Authority
CN
China
Prior art keywords
porosity
core
full
energy
diameter
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.)
Active
Application number
CN202311104089.1A
Other languages
English (en)
Other versions
CN117074274A (zh
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.)
Oil & Gas Survey Cgs
Original Assignee
Oil & Gas Survey Cgs
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 Oil & Gas Survey Cgs filed Critical Oil & Gas Survey Cgs
Priority to CN202311104089.1A priority Critical patent/CN117074274B/zh
Publication of CN117074274A publication Critical patent/CN117074274A/zh
Application granted granted Critical
Publication of CN117074274B publication Critical patent/CN117074274B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
    • G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
    • G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]

Landscapes

  • Health & Medical Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Physics & Mathematics (AREA)
  • Dispersion Chemistry (AREA)
  • Engineering & Computer Science (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pulmonology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Theoretical Computer Science (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Length-Measuring Devices Using Wave Or Particle Radiation (AREA)

Abstract

本发明涉及了一种基于全直径岩心样品的孔隙度计算方法、装置及介质，一种基于全直径岩心样品的孔隙度计算方法，包括：将待扫描的多个全直径岩心样品置于全直径岩心CT扫描仪托盘上，底部放置一玻璃棒作为标定样品；在第一能量下进行全直径岩心CT扫描，获取岩心的第一孔隙度；在第二能量下进行全直径岩心CT扫描，获取岩心的第二孔隙度，其中第一能量小于第二能量；基于两种能量下的CT数据体计算出密度数据体和有效原子序数数据体，计算得到岩心第三孔隙度和岩心第四孔隙度；根据观察判断结果，选取匹配的孔隙度曲线作为孔隙度计算结果。本发明可直接获取全直径岩心的孔隙度结果，数据的结果相对目前已有手段精度更高，代表性更强。

Description

一种基于全直径岩心样品的孔隙度计算方法、装置及介质
技术领域
本发明涉及岩心的孔隙度计算技术领域。尤其涉及一种基于全直径岩心样品的孔隙度计算方法、装置及介质。
背景技术
石油和天然气储存于地下岩石内部，地下岩石储集空间的大小是计算油气资源量的重要参数。地下岩石的储集空间在油气行业内以孔隙度作为标准参数，即地下岩样中所有孔隙空间体积之和与该岩样体积的比值。目前，孔隙度的计算方法主要有地球物理测井法和岩心分析法。
地球物理测井是通过将探测工具，置入油气钻井的井岩中，通过岩石声学、电学和放射性特征，计算岩石的孔隙度。该方法的优势在于能实现大段地层的连续探测，缺点在于为间接探测手段，解释成果的控制因素较多(仪器状态、地下温度、压力等)，精度较低(采样点间距多为0.125米)且成本高。岩心作为地下钻取的实物资料，能直接反映地下储集层的真实信息。但现有手段主要是将地下钻取的原始岩心(又称全直径岩心，直径多为6-10cm直径的圆柱形岩心)上钻取小样(又称柱塞样岩心，直径2.5-3.5cm的圆柱形岩心)置于特定的仪器内测试，结果较为准确，但由于柱塞岩心的样本身体积小，数据代表性有限。
现有技术可以通过基于油气钻井内的地球物理测井获取地下储集体的连续数字，随后基于小型柱塞岩石样品(直径多为2.5cm，长度3-5cm)测试结果标定，从而获取地下不同深度岩石的孔隙度。现有手段中地球物理测井技术对于地下岩石的裂缝表征精度较低，且受仪器性能影响大，成本高。岩心分析则主要基于岩石表面或者小样品分析，通常每米取3块样品开展，对于页岩等较为致密的油气储集体，代表性不够。
发明内容
本发明提供一种基于全直径岩心样品的孔隙度计算方法、装置及介质。以解决现有手段中地球物理测井技术表征精度较低，准确性差，常规的柱塞样岩心测试结果代表性不够的问题。
为实现上述目的，本发明提供一种基于全直径岩心样品的孔隙度计算方法，包括以下步骤：
将待扫描的多个全直径岩心样品置于全直径岩心CT扫描仪托盘上，按所述全直径岩心样品的采样深度由深至浅依次由下至上摆放，在所述全直径岩心样品的底部放置一玻璃棒作为标定样品；
在第一能量下进行全直径岩心CT扫描，获取岩心的第一孔隙度；
在第二能量下进行全直径岩心CT扫描，获取岩心的第二孔隙度，其中所述第一能量小于所述第二能量；
基于两种能量下的CT数据体计算出密度数据体和有效原子序数数据体，基于所述密度数据体计算得到岩心第三孔隙度，基于所述有效原子序数数据体计算得到岩心第四孔隙度；
根据观察判断结果，选取所述第一孔隙度、第二孔隙度、第三孔隙度和第四孔隙度中匹配度高于预设条件的孔隙度曲线作为孔隙度计算结果。
优选的，所述在第一能量下进行全直径岩心CT扫描，具体包括：首先设置第一能量扫描参数，参数为扫描电压80-120kV，电子束流80-100mA，基在第一能量下扫描的CT图像可提取出毫米级的孔隙和裂缝特征，计算出第一孔隙度，所述第一孔隙度为毫米级的孔隙比例。
优选的，所述在第二能量下进行全直径岩心CT扫描，获取岩心的第二孔隙度，具体包括：重置CT扫描仪参数，重置后的参数为扫描电压120-140kV，电子束流100-120mA；基于在第二能量下扫描的CT图像的CT值，计算第二孔隙度P2。
优选的，所述第二孔隙度P2的计算公式为：P2＝a*(CTmax-CTi)/(CTmax-CTlow)；
式中CTmax岩心在所有扫描的深度范围内最大的CT值，单位为Hu；CTi：岩心在某一深度上的CT值，单位为Hu；CTlow为岩在所有扫描的深度范围内最低的CT值，单位为Hu；a为修正因子，根据玻璃棒的CT数调整；不同深度上的CT值为岩心垂长轴方向上切片所有体素点的算术平均值。
优选的，所述密度数据体和所述有效原子序数数据体的计算公式为：
ρe＝A×CThigh+B×CTlow+C
Zen×ρe＝D×(CTlow-CThigh)+E
其中，ρe为样品电子密度，单位g/cm3；②CThigh为第二能量CT扫描图像的CT数，单位为Hounsfield，缩写Hu；③CTlow为第一能量CT扫描图像的CT数，单位为Hu；Ze有效原子序数；n为有效原子序数指数，取值为3.6；A,B,C,D,E均为统计模型系数，其中A、B和D单位为cm3/g，C和E为常数。
优选的，所述第三孔隙度P3计算公式为：P3＝a*(ρmax-ρi)/(ρmax-ρlow)。
式中：ρmax：岩心在所有扫描的深度范围内最大的密度值；ρi：岩心在某一深度上的密度值，单位g/cm3；ρlow：岩心在所有扫描的深度范围内最低的密度值，单位g/cm3；a：修正因子，根据玻璃棒的密度(取值2.65g/cm3)进行调整，无量纲。不同深度上的密度值为岩心垂长轴方向上切片所有体素点的算术平均值；
所述第四孔隙度P4计算公式为：P4＝a*(Zemax-Zei)/(Zemax-Zelow)；式中，Zemax为岩心在所有扫描的深度范围内最大的有效原子序数；Zei为岩心在某一深度上的有效原子序数；Zelow为岩心在所有扫描的深度范围内最低的有效原子序数。a为修正因子，根据玻璃棒的有效原子序数进行调整，取值11.78。
优选的，所述选取所述第一孔隙度、第二孔隙度、第三孔隙度和第四孔隙度中匹配度最高的孔隙度曲线作为孔隙度计算结果之前，还包括：对全直径岩心样品后续钻取了柱塞样岩心并在实验室测试了孔隙度，则基于柱塞样岩心的孔隙度值对所述第一孔隙度、第二孔隙度、第三孔隙度和第四孔隙度进行校正。
优选的，所述根据观察判断结果，选取所述第一孔隙度、第二孔隙度、第三孔隙度和第四孔隙度中匹配度高于预设条件的孔隙度曲线作为孔隙度计算结果，
具体为：基于地球物理测井计算所述多个全直径岩心样品的第五孔隙隙度；
所述第一孔隙度与第二孔隙度的的数值和岩石的流动能力成正比，通过将所述的第一孔隙度、所述第二孔隙度、所述第三孔隙度和所述第四孔隙度分别与所述第五孔隙度对比，选取其中与所述第五孔隙度匹配度最高的作为孔隙度计算结果；
或：根据岩心观察结果，选取所述第一孔隙度、所述第二孔隙度、第三孔隙度和所述第四孔隙度中与所述岩心观察结果匹配较好的孔隙度曲线作为孔隙度计算结果。
第二方面，本发明还涉及一种基于全直径岩心样品的孔隙度计算装置，应用于多个全直径岩心样品进行全直径岩心CT扫描仪扫描的过程中，其中将待扫描的多个全直径岩心样品置于全直径岩心CT扫描仪托盘上，按深度由深至浅依次由下至上摆放，在所述全直径岩心样品的底部放置一玻璃棒作为标定样品，该装置包括：第一能量孔隙度获取模块，用于基于第一能量下进行全直径岩心CT扫描的结果，获取岩心的第一孔隙度；
第二能量孔隙度获取模块，基于第二能量下进行全直径岩心CT扫描的结果，获取岩心的第二孔隙度，其中所述第一能量小于所述第二能量；
密度数据体和有效原子序数数据体获取模块：用于基于两种能量下的CT数据体计算出密度数据体和有效原子序数数据体；
第三第四空隙度获取模块：用于基于所述密度数据体计算得到岩心第三孔隙度，基于所述有效原子序数数据体计算得到岩心第四孔隙度；以及
获取孔隙度结果模块，用于根据观察判断结果，选取所述第一孔隙度、第二孔隙度、第三孔隙度和第四孔隙度中匹配度高于预设条件的孔隙度曲线作为孔隙度计算结果。
第三方面，本发明还涉及一种计算机可读存储介质，存储介质中存储有指令，所述指令运行时执行上述中任一项所述的一种基于全直径岩心样品的孔隙度计算方法。
本发明涉及的一种基于全直径岩心样品的孔隙度计算方法、装置及介质，相比于现有技术，有以下有益效果：
可直接获取全直径岩心的孔隙度结果，数据的结果相对目前已有手段精度更高，代表性更强。
孔隙度值的间距为数百微米(不同CT设备略有不同，通常在65微米-500微米)，高分辨率的表征结果可准确、有效地确定页岩等致密储层的有利开发层段。
附图说明
图1为本发明实施例一中一种基于全直径岩心样品的孔隙度计算方法的方法流程图一；
图2是本发明实施例一中第一孔隙度、第二孔隙度、第三空隙度以及第四孔隙度的曲线示意图；
图3是本发明实施例二中一种基于全直径岩心样品的孔隙度计算装置的结构示意图；
图4为本发明实施例一种关于柱塞样岩心在实验室测试孔隙度的折线图；
图5为本发明实施例一或二中全直径岩心CT扫描仪托盘的结构示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释本发明，而非对本发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与本发明相关的部分而非全部结构。
实施例一
一种基于全直径岩心样品的孔隙度计算方法，请参阅图1-2、4-5，包括以下步骤：
S10:将待扫描的多个全直径岩心样品置于全直径岩心CT扫描仪托盘上，按全直径岩心样品的采样深度由深至浅依次由下至上摆放，在全直径岩心样品的底部放置一玻璃棒作为标定样品。其中，全直径岩心的取样是用特定工具在地下连续取样，短则连续或不连续数十米，长则连续取数百米。
具体的，对全直径岩心样品进行整理，将待扫描的全直径岩心样品整理好后，置于全直径岩心CT扫描仪托盘上，按全直径岩心样品的深度由深至浅，依次由下至上摆放。在全直径岩心样品样品的底部放置一玻璃棒，作为标定样品。岩心样品摆放完毕后，开始启动CT扫描仪。
S20:在第一能量下进行全直径岩心CT扫描，获取岩心的第一孔隙度。
具体包括：首先设置第一能量扫描参数，本实施例中第二能量参数为扫描电压80-120kV，电子束流80-100mA，基于在第一能量下扫描的CT图像可提取出毫米级的孔隙和裂缝特征，计算出第一孔隙度P1，第一孔隙度为毫米级的孔隙比例。
首先设置低能量扫描参数，对岩心进行扫描，建议的参数为扫描电压80-120kV，电子束流80-100mA。本次扫描将获取第一能量下的三维岩石图像，基于CT图像可提取出毫米级的孔隙和裂缝特征，计算出毫米级的孔隙比例，即第一孔隙度P1。
S30:在第二能量下进行全直径岩心CT扫描，获取岩心的第二孔隙度，其中,第一能量小于第二能量。
具体的：重置CT扫描仪参数，重置后的第二能量参数为扫描电压120-140kV，电子束流100-120mA；基于在第二能量下扫描的CT图像的CT值，计算第二孔隙度P2。
本实施例中，在完成低能量岩心扫描后，重置CT扫描仪参数，开展高能量下的岩心CT扫描。建议的参数为，扫描电压120-140kV，电子束流100-120mA。本次扫描将获取另一种CT图像，根据CT图像的CT值，计算第二孔隙度P2。CT值是表征X射线穿透不同物质后的衰减程度的数值，不同密度的物质CT值不同。通常以水和空气作为标准，其中，水的CT值为0，空气的CT值为-1000。
第二孔隙度P2的计算公式为：P2＝a*(CTmax-CTi)/(CTmax-CTlow)；
式中CTmax岩心在所有扫描的深度范围内最大的CT值，单位为Hu；CTi：岩心在某一深度上的CT值，单位为Hu；CTlow为岩在所有扫描的深度范围内最低的CT值，单位为Hu；a为修正因子，根据玻璃棒的CT数调整,本申请中取值1500，无量纲；不同深度上的CT值为岩心垂长轴方向上切片所有体素点的算术平均值。
S40:基于两种能量下的CT数据体计算出密度数据体和有效原子序数数据体，基于密度数据体计算得到岩心第三孔隙度P3，基于有效原子序数数据体计算得到岩心第四孔隙度P4。
本实施例中，密度数据体和有效原子序数数据体的计算公式为：
ρe＝A×CThigh+B×CTlow+C
Zen×ρe＝D×(CTlow-CThigh)+E
其中，ρe为样品电子密度，单位g/cm3；②CThigh为第二能量CT扫描图像的CT数，单位为Hounsfield，缩写Hu；③CTlow为第一能量CT扫描图像的CT数，单位为Hu；Ze有效原子序数；n为有效原子序数指数，取值为3.6；A,B,C,D,E均为统计模型系数，其中A、B和D单位为cm3/g，C和E为常数。
第三孔隙度P3计算公式为：P3＝a*(ρmax-ρi)/(ρmax-ρlow)。
式中：ρmax：岩心在所有扫描的深度范围内最大的密度值；ρi：岩心在某一深度上的密度值，单位g/cm3；ρlow：岩心在所有扫描的深度范围内最低的密度值，单位g/cm3；a：修正因子，根据玻璃棒的密度进行调整，本实施例中取值2.65g/cm3，无量纲。不同深度上的密度值为岩心垂长轴方向上切片所有体素点的算术平均值。
第四孔隙度P4计算公式为：P4＝a*(Zemax-Zei)/(Zemax-Zelow)；式中，Zemax为岩心在所有扫描的深度范围内最大的有效原子序数；Zei为岩心在某一深度上的有效原子序数；Zelow为岩心在所有扫描的深度范围内最低的有效原子序数。a为修正因子，根据玻璃棒的有效原子序数进行调整，取值11.78，无量纲。不同深度上的有效原子序数为岩心垂长轴方向上切片所有体素点的算术平均值。
S50:根据观察判断结果，选取第一孔隙度、第二孔隙度、第三孔隙度和第四孔隙度中匹配度高于预设条件的孔隙度曲线作为孔隙度计算结果。
一些实施例中，观察结果有两种情况：1.无标定数据时，仅凭四种孔隙度的曲线形态(即不同深度的数值变化)，选择匹配形态较好的2条以上曲线，将两条或两条以上曲线形态取算术平均值。2.有实验室获取的不同深度岩心数据点孔隙度值或地球物理测井曲线时，选取四条孔隙度值中和岩心数据点孔隙度值或地球物理测井曲线匹配较好的曲线。
S50具体为：基于地球物理测井计算多个全直径岩心样品的第五孔隙度；
第一孔隙度与第二孔隙度的的数值和岩石的流动能力成正比，通过将第一孔隙度、第二孔隙度、第三孔隙度和第四孔隙度分别与第五孔隙度对比，选取其中与第五孔隙度匹配度最高的作为孔隙度计算结果。
或：S50具体为，根据岩心观察结果，选取第一孔隙度、第二孔隙度、第三孔隙度和第四孔隙度中与岩心观察结果匹配较好的孔隙度曲线作为孔隙度计算结果。
在一些实施例中，在S40之后，S50之前还包括：S60:对全直径岩心样品后续钻取了柱塞样岩心并在实验室测试了孔隙度，则基于柱塞样岩心的孔隙度值对第一孔隙度、第二孔隙度、第三孔隙度和第四孔隙度进行校正。
例如，如图4所示：任一孔隙度曲线校正方法：
假设孔隙度曲线中不同深度的数值为Pi，令校正后的孔隙度为Yi＝aPi+b，将所有岩心数据Zi点统计出来，计算不同深度点测井值与岩心值的差值平方，计算出使N值最小时的a值和b值，完成孔隙度校正。
实施例二
如图2-3、5所示，一种基于全直径岩心样品的孔隙度计算装置，本实施例的基于全直径岩心样品的孔隙度计算装置应用于多个全直径岩心样品进行全直径岩心CT扫描仪扫描的过程中，其中，将待扫描的多个全直径岩心样品置于全直径岩心CT扫描仪托盘上，按深度由深至浅依次由下至上摆放，在全直径岩心样品的底部放置一玻璃棒作为标定样品，该装置可为带有中央处理器的控制中心，如个人电脑、上位机、服务器等电子设备实现。
实施例中，如图3所示，基于全直径岩心样品的孔隙度计算装置还包括：第一能量孔隙度获取模块71、第二能量孔隙度获取模块72、密度数据体和有效原子序数数据体获取模块73、第三第四空隙度获取模块74以及获取孔隙度结果模块75。
第一能量孔隙度获取模块71，用于基于第一能量下进行全直径岩心CT扫描的结果，获取岩心的第一孔隙度；
第二能量孔隙度获取模块72，基于第二能量下进行全直径岩心CT扫描的结果，获取岩心的第二孔隙度，其中第一能量小于第二能量；
密度数据体和有效原子序数数据体获取模块73：用于基于两种能量下的CT数据体计算出密度数据体和有效原子序数数据体；
第三第四空隙度获取模块74：用于基于密度数据体计算得到岩心第三孔隙度，基于有效原子序数数据体计算得到岩心第四孔隙度；
获取孔隙度结果模块75，用于根据观察判断结果，选取第一孔隙度、第二孔隙度、第三孔隙度和第四孔隙度中匹配度高于预设条件的孔隙度曲线作为孔隙度计算结果。
本实施例的一种基于全直径岩心样品的孔隙度计算装置，其实现过程方法和效果均与实施例一描述的一种基于全直径岩心样品的孔隙度计算方法相同，在此不再赘述。
实施例三
本发明涉及一种计算机可读存储介质，存储介质中存储有指令，指令运行时执行实施一的一种基于全直径岩心样品的孔隙度计算方法，其运行时实现过程方法和效果均与实施例一描述的一种基于全直径岩心样品的孔隙度计算方法相同，在此不再赘述。
需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者装置不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者装置所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括该要素的过程、方法、物品或者装置中还存在另外的相同要素。
以上仅为本发明的优选实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其他相关的技术领域，均同理包括在本发明的专利保护范围内。

Claims (9)

1.一种基于全直径岩心样品的岩石孔隙度计算方法，其特征在于，包括以下步骤：
将待扫描的多个全直径岩心样品置于全直径岩心CT扫描仪托盘上，按所述全直径岩心样品的采样深度由深至浅依次由下至上摆放，在所述全直径岩心样品的底部放置一玻璃棒作为标定样品；
在第一能量下进行全直径岩心CT扫描，获取岩心的第一孔隙度；
在第二能量下进行全直径岩心CT扫描，获取岩心的第二孔隙度，其中所述第一能量小于所述第二能量；
基于两种能量下的CT数据体计算出密度数据体和有效原子序数数据体，基于所述密度数据体计算得到岩心第三孔隙度，基于所述有效原子序数数据体计算得到岩心第四孔隙度；
根据观察判断结果，选取所述第一孔隙度、所述第二孔隙度、所述第三孔隙度和所述第四孔隙度中匹配度高于预设条件的孔隙度曲线作为孔隙度计算结果；
所述第三孔隙度P3计算公式为：P3＝a*(ρmax-ρi)/(ρmax-ρlow)；
式中：ρmax：岩心在所有扫描的深度范围内最大的密度值；ρi：岩心在某一深度上的密度值，单位g/cm3；ρlow：岩心在所有扫描的深度范围内最低的密度值，单位g/cm3；a：修正因子，根据玻璃棒的密度进行调整，无量纲；不同深度上的密度值为岩心垂长轴方向上切片所有体素点的算术平均值；
所述第四孔隙度P4计算公式为：P4＝a*(Zemax-Zei)/(Zemax-Zelow)；
式中，Zemax为岩心在所有扫描的深度范围内最大的有效原子序数；Zei为岩心在某一深度上的有效原子序数；Zelow为岩心在所有扫描的深度范围内最低的有效原子序数，a为修正因子，根据玻璃棒的有效原子序数进行调整，其中，所述玻璃棒的有效原子序数为11.78。
2.根据权利要求1所述的一种基于全直径岩心样品的孔隙度计算方法，其特征在于，所述在第一能量下进行全直径岩心CT扫描，具体包括：首先设置第一能量扫描参数，参数为扫描电压80-120kV，电子束流80-100mA，基在第一能量下扫描的CT图像可提取出毫米级的孔隙和裂缝特征，计算出第一孔隙度，所述第一孔隙度为毫米级的孔隙比例。
3.根据权利要求2所述的一种基于全直径岩心样品的孔隙度计算方法，其特征在于，所述在第二能量下进行全直径岩心CT扫描，获取岩心的第二孔隙度，具体包括：重置CT扫描仪参数，重置后的参数为扫描电压120-140kV，电子束流100-120mA；基于在第二能量下扫描的CT图像的CT值，计算第二孔隙度P2。
4.根据权利要求3所述的一种基于全直径岩心样品的孔隙度计算方法，其特征在于，所述第二孔隙度P2的计算公式为：P2＝a*(CTmax-CTi)/(CTmax-CTlow)；
式中CTmax为岩心在所有扫描的深度范围内最大的CT值，单位为Hu；CTi：岩心在某一深度上的CT值，单位为Hu；CTlow为岩在所有扫描的深度范围内最低的CT值，单位为Hu；a为修正因子，根据玻璃棒的CT数调整；不同深度上的CT值为岩心垂长轴方向上切片所有体素点的算术平均值。
5.根据权利要求1所述的一种基于全直径岩心样品的孔隙度计算方法，其特征在于，所述密度数据体和所述有效原子序数数据体的计算公式为：
Pe＝A×CThigh+B×CTlow+C
Zen×ρe＝D×(CTlow-CThigh)+E
其中，ρe为样品电子密度，单位g/cm3；②CThigh为第二能量CT扫描图像的CT数，单位为Hounsfield，缩写Hu；③CTlow为第一能量CT扫描图像的CT数，单位为Hu；Ze有效原子序数；n为有效原子序数指数，取值为3.6；A,B,C,D,E均为统计模型系数，其中A、B和D单位为cm3/g，C和E为常数。
6.根据权利要求1所述的一种基于全直径岩心样品的孔隙度计算方法，其特征在于，所述选取所述第一孔隙度、第二孔隙度、第三孔隙度和第四孔隙度中匹配度最高的孔隙度曲线作为孔隙度计算结果之前，还包括：对全直径岩心样品后续钻取了柱塞样岩心并在实验室测试了孔隙度，则基于柱塞样岩心的孔隙度值对所述第一孔隙度、第二孔隙度、第三孔隙度和第四孔隙度进行校正。
7.根据权利要求1所述的一种基于全直径岩心样品的孔隙度计算方法，其特征在于：所述根据观察判断结果，选取所述第一孔隙度、第二孔隙度、第三孔隙度和第四孔隙度中匹配度高于预设条件的孔隙度曲线作为孔隙度计算结果，
具体为：基于地球物理测井计算所述多个全直径岩心样品的第五孔隙度；
所述第一孔隙度与第二孔隙度的的数值和岩石的流动能力成正比，通过将所述第一孔隙度、所述第二孔隙度、所述第三孔隙度和所述第四孔隙度分别与所述第五孔隙度对比，选取其中与所述第五孔隙度匹配度最高的作为孔隙度计算结果；
或：根据岩心观察结果，选取所述第一孔隙度、所述第二孔隙度、所述第三孔隙度和所述第四孔隙度中与所述岩心观察结果匹配较好的孔隙度曲线作为孔隙度计算结果。
8.一种基于全直径岩心样品的孔隙度计算装置，其特征在于：应用于多个全直径岩心样品进行全直径岩心CT扫描仪扫描的过程中，其中将待扫描的多个全直径岩心样品置于全直径岩心CT扫描仪托盘上，按深度由深至浅依次由下至上摆放，在所述全直径岩心样品的底部放置一玻璃棒作为标定样品，该装置包括：第一能量孔隙度获取模块，用于基于第一能量下进行全直径岩心CT扫描的结果，获取岩心的第一孔隙度；
第二能量孔隙度获取模块，基于第二能量下进行全直径岩心CT扫描的结果，获取岩心的第二孔隙度，其中所述第一能量小于所述第二能量；
密度数据体和有效原子序数数据体获取模块：用于基于两种能量下的CT数据体计算出密度数据体和有效原子序数数据体；
第三第四空隙度获取模块：用于基于所述密度数据体计算得到岩心第三孔隙度，基于所述有效原子序数数据体计算得到岩心第四孔隙度；所述第三孔隙度P3计算公式为：P3＝a*(ρmax-ρi)/(ρmax-ρlow)；
式中：ρmax：岩心在所有扫描的深度范围内最大的密度值；ρi：岩心在某一深度上的密度值，单位g/cm3；ρlow：岩心在所有扫描的深度范围内最低的密度值，单位g/cm3；a：修正因子，根据玻璃棒的密度进行调整，无量纲；不同深度上的密度值为岩心垂长轴方向上切片所有体素点的算术平均值；
所述第四孔隙度P4计算公式为：P4＝a*(Zemax-Zei)/(Zemax-Zelow)；
式中，Zemax为岩心在所有扫描的深度范围内最大的有效原子序数；Zei为岩心在某一深度上的有效原子序数；Zelow为岩心在所有扫描的深度范围内最低的有效原子序数，a为修正因子，根据玻璃棒的有效原子序数进行调整；以及
获取孔隙度结果模块，用于根据观察判断结果，选取所述第一孔隙度、第二孔隙度、第三孔隙度和第四孔隙度中匹配度高于预设条件的孔隙度曲线作为孔隙度计算结果。
9.一种计算机可读存储介质，其特征在于：存储介质中存储有指令，所述指令运行时执行如权利要求1-7中任一项所述的一种基于全直径岩心样品的孔隙度计算方法。
CN202311104089.1A 2023-08-30 2023-08-30 一种基于全直径岩心样品的孔隙度计算方法、装置及介质 Active CN117074274B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311104089.1A CN117074274B (zh) 2023-08-30 2023-08-30 一种基于全直径岩心样品的孔隙度计算方法、装置及介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311104089.1A CN117074274B (zh) 2023-08-30 2023-08-30 一种基于全直径岩心样品的孔隙度计算方法、装置及介质

Publications (2)

Publication Number Publication Date
CN117074274A CN117074274A (zh) 2023-11-17
CN117074274B true CN117074274B (zh) 2024-05-14

Family

ID=88705908

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311104089.1A Active CN117074274B (zh) 2023-08-30 2023-08-30 一种基于全直径岩心样品的孔隙度计算方法、装置及介质

Country Status (1)

Country Link
CN (1) CN117074274B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4982086A (en) * 1988-07-14 1991-01-01 Atlantic Richfield Company Method of porosity determination in porous media by x-ray computed tomography
CN103718016A (zh) * 2011-07-26 2014-04-09 因格瑞恩股份有限公司 用双能x射线计算机断面显像技术估算岩样的有效原子序数和体积密度的方法
CN104169714A (zh) * 2012-01-13 2014-11-26 领英股份有限公司 用多能量x射线成像确定储层性质和品质的方法
WO2015034472A1 (en) * 2013-09-04 2015-03-12 Ingrain, Inc. Combining multiple energy x-ray imaging and well data to obtain high-resolution rock, mechanical, and elastic property profiles
CN112576244A (zh) * 2020-11-17 2021-03-30 中海油田服务股份有限公司 一种孔隙度计算方法和装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4982086A (en) * 1988-07-14 1991-01-01 Atlantic Richfield Company Method of porosity determination in porous media by x-ray computed tomography
CN103718016A (zh) * 2011-07-26 2014-04-09 因格瑞恩股份有限公司 用双能x射线计算机断面显像技术估算岩样的有效原子序数和体积密度的方法
CN104169714A (zh) * 2012-01-13 2014-11-26 领英股份有限公司 用多能量x射线成像确定储层性质和品质的方法
WO2015034472A1 (en) * 2013-09-04 2015-03-12 Ingrain, Inc. Combining multiple energy x-ray imaging and well data to obtain high-resolution rock, mechanical, and elastic property profiles
CN112576244A (zh) * 2020-11-17 2021-03-30 中海油田服务股份有限公司 一种孔隙度计算方法和装置

Also Published As

Publication number Publication date
CN117074274A (zh) 2023-11-17

Similar Documents

Publication Publication Date Title
CN102628354B (zh) 孔隙微米级油水分布识别量化方法
US5164590A (en) Method for evaluating core samples from x-ray energy attenuation measurements
US4982086A (en) Method of porosity determination in porous media by x-ray computed tomography
Muszynski et al. Particle shape estimates of uniform sands: visual and automated methods comparison
CN110927035A (zh) 一种低渗致密砂岩束缚水饱和度计算方法
CN107941670B (zh) 一种岩屑孔隙度测定方法
CN103353462A (zh) 一种基于核磁共振成像的岩石非均质性定量评价方法
CN103091226A (zh) 一种检测饱和土孔隙率装置及方法
CN110320137A (zh) 一种基于数字岩心的多尺度融合方法
CN117074274B (zh) 一种基于全直径岩心样品的孔隙度计算方法、装置及介质
CN109444190B (zh) 利用数字岩心分析提高复杂油气层含水饱和度评价精度的方法
CN114235862A (zh) 一种基于双能ct对烃源岩的评估方法
CN114460119A (zh) 力学参数获取方法及装置、电子设备及计算机存储介质
WO2022204297A1 (en) The combined discrete gaussian analysis of micp and nmr t2 distributions of multi-modal carbonate rocks
Lundy et al. Measurement of snow density and microstructure using computed tomography
Freitag et al. Petrophysical characterization, BIB-SEM imaging, and permeability models of tight carbonates from the Upper Jurassic (Malm ß), SE Germany
CN114518309A (zh) 一种页岩油储层中夹层孔隙有效三维表征的方法及装置
Le Roux et al. MicroCT-based bulk density measurement method for soils
CN114428049B (zh) 一种计算古老碳酸盐岩储层沥青含量的方法
Al-Bulushi et al. Capillary Pressure Corrections, Quality Control and Curve Fitting Workflow
CN111693427A (zh) 油气藏流体可动性的分析方法
CN113405962B (zh) 碳酸盐岩储层的储集空间评价方法和装置
CN112160307B (zh) 基于多元采样属性分析的岩土取样布局和数据分析***
CN118150407A (zh) 一种基于三维数字岩心定量表征砂岩润湿性的方法
Mahmoud et al. Evaluating Flow Characteristics of Reservoir Rock Using Digital Rock Physics (DRP)

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