CN114518309A - 一种页岩油储层中夹层孔隙有效三维表征的方法及装置 - Google Patents
一种页岩油储层中夹层孔隙有效三维表征的方法及装置 Download PDFInfo
- Publication number
- CN114518309A CN114518309A CN202011314324.4A CN202011314324A CN114518309A CN 114518309 A CN114518309 A CN 114518309A CN 202011314324 A CN202011314324 A CN 202011314324A CN 114518309 A CN114518309 A CN 114518309A
- Authority
- CN
- China
- Prior art keywords
- porosity
- rock sample
- pores
- threshold value
- obtaining
- 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
- 239000011148 porous material Substances 0.000 title claims abstract description 196
- 238000000034 method Methods 0.000 title claims abstract description 45
- 239000003079 shale oil Substances 0.000 title claims abstract description 37
- 239000011229 interlayer Substances 0.000 title claims abstract description 26
- 239000011435 rock Substances 0.000 claims abstract description 141
- 238000012545 processing Methods 0.000 claims abstract description 42
- 238000009826 distribution Methods 0.000 claims abstract description 34
- 238000001914 filtration Methods 0.000 claims abstract description 22
- 230000011218 segmentation Effects 0.000 claims abstract description 13
- 238000012512 characterization method Methods 0.000 claims description 25
- 239000001307 helium Substances 0.000 claims description 22
- 229910052734 helium Inorganic materials 0.000 claims description 22
- SWQJXJOGLNCZEY-UHFFFAOYSA-N helium atom Chemical compound [He] SWQJXJOGLNCZEY-UHFFFAOYSA-N 0.000 claims description 22
- 238000005259 measurement Methods 0.000 claims description 14
- 238000004626 scanning electron microscopy Methods 0.000 claims 2
- 238000004364 calculation method Methods 0.000 abstract description 5
- 238000005516 engineering process Methods 0.000 description 6
- 238000011161 development Methods 0.000 description 4
- 230000018109 developmental process Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000012216 screening Methods 0.000 description 4
- 239000007789 gas Substances 0.000 description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- XKRFYHLGVUSROY-UHFFFAOYSA-N argon Substances [Ar] XKRFYHLGVUSROY-UHFFFAOYSA-N 0.000 description 2
- 229910052786 argon Inorganic materials 0.000 description 2
- -1 argon ion Chemical class 0.000 description 2
- 239000010459 dolomite Substances 0.000 description 2
- 229910000514 dolomite Inorganic materials 0.000 description 2
- 238000000265 homogenisation Methods 0.000 description 2
- 238000005498 polishing Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 235000019738 Limestone Nutrition 0.000 description 1
- 241000425037 Toona sinensis Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000006028 limestone Substances 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 1
- 229910052753 mercury Inorganic materials 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
Images
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
- 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
-
- 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
- G01N2015/0846—Investigating permeability, pore-volume, or surface area of porous materials by use of radiation, e.g. transmitted or reflected light
Landscapes
- Chemical & Material Sciences (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明公开了一种页岩油储层中夹层孔隙有效三维表征的方法及装置，该方法包括：获取岩石样品的实测孔隙度；将岩石样品制备成与CT设备对应的岩石样品；获取与CT设备对应的岩石样品的灰度图像；设置图像处理的灰度阈值，对滤波降噪处理后的灰度图像进行分割和二值化处理，初步识别岩石样品的孔隙和骨架；基于初步识别的孔隙和骨架，获得岩石样品的计算孔隙度；以实测孔隙度值为标定，调整灰度阈值，获得岩石样品的孔隙度和孔径分布特征。本发明以岩石样品的孔隙度实测值为标定，在图像处理中调整灰度阈值，进而使得计算获得的孔隙度近似等于孔隙度实测值，实现对页岩油储层非泥岩夹层孔隙的有效真实反应，有效地避免CT应用技术局限性。
Description
技术领域
本发明属于油气勘探与开发技术领域，具体涉及一种页岩油储层中夹层孔隙有效三维表征的方法及装置。
背景技术
近年来，伴随石油工业的发展，非常规页岩油气成为重要的勘探开发项目。伴随着国内外学者对北美典型页岩油藏的剖析，页岩油储层中非泥岩夹层(砂岩、碳酸盐岩)类似于常规油藏，常作为优势岩相，具有物性好、连通性优、含油性佳的特点。由此，非泥岩夹层孔隙的有效表征及研究对页岩油气的勘探和开发具有重要的意义。目前对识别出的孔隙孔隙结构的全息表征从常规的岩心观察描述、薄片鉴定及氩离子抛光扫描等定性描述发展到低温液氮吸附—高压压汞联合、微纳米CT及FIB-SEM三维构建等定量分析。其中CT技术在对页岩层系储层(包括泥/页岩、砂岩夹层、白云岩夹层和灰岩夹层)进行三维表征方面具有巨大潜力。同时结合图像处理技术(Avizo等)可以对识别出的孔隙连通性、各向异性及孔隙度等进行定量评价。
苟启洋等(2018)利用纳米CT及三维重构评价了川东南焦石坝地区JY1井下志留统龙马溪组页岩的孔隙特征及连通性，其中基于纳米CT扫描所计算孔隙约为2.71％，小于前人的研究成果(3.90％和4.40％)。朱如凯等(2016)采用纳米CT(研究尺度100-6500nm)对四川盆地志留系龙马溪组页岩进行三维构建及定量分析，同时将纳米CT计算孔隙度与氦气法孔隙进行比较发现，计算孔隙度2.2％小于氦气法孔隙度5.6％。而造成这种情况的原因，主要为：(1)页岩样品的非均质化；(2)灰度图像处理中孔隙灰度值范围的设定不合理；(3)CT测量范围受制于制备样品尺寸，太过局限性。在此背景下，这就使得所表征的孔隙结构具有局限性，其有效性值得商榷。而对于非泥岩夹层而言，类似于常规油藏储层，非均质化的影响要小得多，其影响更多的是，灰度图像处理中灰度值的调整以及制备样品尺寸。
为此，亟需探索一种新方法以弥补以上CT技术表征非泥岩夹层孔隙结构所存在的不足，以使得非泥岩夹层孔隙得到有效表征。
发明内容
本发明的目的是提出一种使得非泥岩夹层孔隙得到有效表征的页岩油储层中夹层孔隙有效三维表征的方法及装置。
第一方面，本发明提供一种页岩油储层中夹层孔隙有效三维表征的方法，包括：获取甄选后的岩石样品的实测孔隙度；将所述甄选后的岩石样品制备成与CT设备对应的岩石样品；获取与所述CT设备对应的岩石样品的灰度图像；对所述灰度图像进行滤波降噪处理；设置图像处理的灰度阈值，对所述滤波降噪处理后的灰度图像进行分割和二值化处理，初步识别所述岩石样品的孔隙和骨架；基于初步识别的孔隙和骨架，获得所述岩石样品的计算孔隙度；基于所述计算孔隙度与所述实测孔隙度值的差值，以所述实测孔隙度值为标定，调整灰度阈值，获取最终识别的孔隙和骨架；基于所述最终识别的孔隙和骨架，获得所述岩石样品的孔隙度和孔径分布特征。
所述基于所述计算孔隙度与所述实测孔隙度值的差值，以所述实测孔隙度值为标定，调整灰度阈值，获取最终识别的孔隙和骨架包括：计算所述计算孔隙度和实测孔隙度的差值，基于所述差值，以所述实测孔隙度值为标定，调整灰度阈值，获得新的计算孔隙度；重复此过程直至所述差值小于预设阈值时；在最后一次调整的灰度阈值下，获取最终识别的孔隙和骨架。
可选的，所述基于所述最终识别的孔隙和骨架，获得所述岩石样品的孔隙度和孔径分布特征包括：基于所述最终识别的孔隙和骨架，建立所述岩石样品的三维数字模型；基于所述三维数字模型，获得所述岩石样品的孔隙度和孔径分布特征。
可选的，将通过薄片鉴定或扫描电镜选取的均一性岩石样品作为甄选后的岩石样品。
可选的，通过氦气法孔隙度测试仪获取甄选后的岩石样品的实测孔隙度。
可选的，所述将所述甄选后的岩石样品制备成与CT设备对应的岩石样品包括：根据所述CT设备的测量范围，将所述甄选后的岩石样品制备成与所述CT设备对应的尺寸。
可选的，所述CT设备为微米CT设备。
第二方面，本发明还提供一种页岩油储层中夹层孔隙有效三维表征的装置，包括：氦气法孔隙度测试仪，获取甄选后的岩石样品的氦气法实测孔隙度；CT设备，获取与所述CT设备对应的岩石样品的灰度图像；处理器，所述处理器与所述CT设备连接，所述处理器执行以下步骤：对所述灰度图像进行滤波降噪处理；设置图像处理的灰度阈值，对所述滤波降噪处理后的灰度图像进行分割和二值化处理，初步识别所述岩石样品的孔隙和骨架；基于初步识别的孔隙和骨架，获得所述岩石样品的计算孔隙度；基于所述计算孔隙度与所述实测孔隙度值的差值，以所述实测孔隙度值为标定，调整灰度阈值，获取最终识别的孔隙和骨架；基于所述最终识别的孔隙和骨架，获得所述岩石样品的孔隙度和孔径分布特征。
所述基于所述计算孔隙度与所述实测孔隙度值的差值，以所述实测孔隙度值为标定，调整灰度阈值，获取最终识别的孔隙和骨架包括：计算所述计算孔隙度和实测孔隙度的差值，基于所述差值，以所述实测孔隙度值为标定，调整灰度阈值，获得新的计算孔隙度；重复此过程直至所述差值小于预设阈值时；在最后一次调整的灰度阈值下，获取最终识别的孔隙和骨架。
可选的，所述基于所述最终识别的孔隙和骨架，获得所述岩石样品的孔隙度和孔径分布特征包括：基于所述最终识别的孔隙和骨架，建立所述岩石样品的三维数字模型；基于所述三维数字模型，获得所述岩石样品的孔隙度和孔径分布特征。
可选的，将通过薄片鉴定或扫描电镜选取的均一性岩石样品作为甄选后的岩石样品。
可选的，根据所述CT设备的测量范围，将所述甄选后的岩石样品制备成与所述CT设备对应的尺寸。
可选的，所述CT设备为微米CT设备。
可选的，所述氦气法孔隙度测试仪执行GB/T34533-2017标准获取甄选后的岩石样品的氦气法实测孔隙。
本发明的有益效果在于：本发明的页岩油储层中夹层孔隙有效三维表征的方法通过遴选均一化夹层岩石样品，以岩石样品的孔隙度实测值为标定，在图像处理中调整灰度阈值，进而使得计算获得的孔隙度近似等于孔隙度实测值，实现对页岩油储层非泥岩夹层孔隙的有效真实反应，有效地避免前人在应用CT技术时的局限性。
本发明具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施例中将是显而易见的，或者将在并入本文中的附图和随后的具体实施例中进行详细陈述，这些附图和具体实施例共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施方式进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施方式中，相同的参考标号通常代表相同部件。
图1示出了根据本发明的一个实施例的一种页岩油储层中夹层孔隙有效三维表征的方法的流程图。
图2示出了根据本发明的一个实施例的一种页岩油储层中夹层孔隙有效三维表征的方法的孔隙结构三维特征。
图3示出了根据本发明的一个实施例的一种页岩油储层中夹层孔隙有效三维表征的方法的孔隙分布特征图。
图4示出了根据本发明的一个实施例的一种页岩油储层中夹层孔隙有效三维表征的装置的框图。
102、氦气法孔隙度测试仪；104、CT设备；106、处理器。
具体实施方式
下面将更详细地描述本发明的优选实施方式。虽然以下描述了本发明的优选实施方式，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。
本发明提供一种页岩油储层中夹层孔隙有效三维表征的方法，包括：获取甄选后的岩石样品的实测孔隙度；将甄选后的岩石样品制备成与CT设备对应的岩石样品；对灰度图像进行滤波降噪处理；设置图像处理的灰度阈值，对滤波降噪处理后的灰度图像进行分割和二值化处理，初步识别岩石样品的孔隙和骨架；基于初步识别的孔隙和骨架，获得岩石样品的计算孔隙度；基于计算孔隙度与实测孔隙度值的差值，以实测孔隙度值为标定，调整灰度阈值，获取最终识别的孔隙和骨架；基于最终识别的孔隙和骨架，获得岩石样品的孔隙度和孔径分布特征。
具体的，通过氦气法测定甄选后的岩石样品的孔隙度实测值，将甄选后的岩石样品制备成与CT设备对应的样品尺寸，CT设备获取对应岩石样品的灰度图像，基于Avizo软件对灰度图像进行滤波降噪处理后，设置图像处理的灰度阈值，以使得计算孔隙度近似等于实测孔隙度，然后对所述灰度图像进行阈值分割和二值化，识别出孔隙和骨架；基于孔隙和骨架，对识别出的孔隙进行数据统计，获得岩石样品的计算孔隙度，以实测孔隙度值作为标定，调整灰度阈值，重新识别孔隙和骨架，获得计算孔隙度，当计算孔隙度近似于实测孔隙度时，此时即确定识别出的孔隙、骨架为最终识别的结果。基于最终识别出的孔隙、骨架，建立岩石样品的三维数字模型，根据三维数字模型，获得岩石样品的孔隙度和孔径分布特征。
根据示例性的实施方式，页岩油储层中夹层孔隙有效三维表征的方法通过遴选均一化夹层岩石样品，以岩石样品的孔隙度实测值为标定，在图像处理中调整灰度阈值，进而使得计算获得的孔隙度近似等于孔隙度实测值，实现对页岩油储层非泥岩夹层孔隙的有效真实反应，有效地避免前人在应用CT技术时的局限性。
作为可选方案，基于计算孔隙度与实测孔隙度值的差值，以实测孔隙度值为标定，调整灰度阈值，获取最终识别的孔隙和骨架包括：计算计算孔隙度和实测孔隙度的差值，基于差值，以实测孔隙度值为标定，调整灰度阈值，获得新的计算孔隙度；重复此过程直至差值小于预设阈值时；在最后一次调整的灰度阈值下，获取最终识别的孔隙和骨架。
具体的，先对灰度图像进行滤波降噪处理，然后设置图像处理的灰度阈值，对图像进行阈值分割和二值化，识别出孔隙、骨架；对识别出的孔隙进行数据统计，计算输出岩石样品的孔隙度，以实测孔隙度值作为标定，调整灰度阈值，当计算孔隙度近似于实测孔隙度，此时即确定识别出的孔隙、骨架。
作为可选方案，基于最终识别的孔隙和骨架，获得岩石样品的孔隙度和孔径分布特征包括：基于最终识别的孔隙和骨架，建立岩石样品的三维数字模型；基于三维数字模型，获得岩石样品的孔隙度和孔径分布特征。
具体的，通过选用最终识别的孔隙和骨架的数百张图像，将这些二维CT图像进行堆叠排列并重构，形成三维数字模型，通过三维数字模型输出孔隙结构的三维特征；基于分水岭算法识别孔隙边界并对各孔隙进行标记，然后对各孔隙进行单独分析，通过统计各孔隙中的体素个数及边界上的像素个数来表征孔隙的体积，并按孔隙大小(即等效直径)分别对孔隙数量、体积的分布进行统计，计算输出孔隙度、孔径分布特征。
作为可选方案，将通过薄片鉴定或扫描电镜选取的均一性岩石样品作为甄选后的岩石样品。
具体的，中通过薄片鉴定、扫描电镜尽量遴选均一性岩石样品。
作为可选方案，通过氦气法孔隙度测试仪获取甄选后的岩石样品的实测孔隙度。
具体的，氦气法孔隙度测定执行GB/T34533-2017标准。
作为可选方案，将甄选后的岩石样品制备成与CT设备对应的岩石样品包括：根据CT设备的测量范围，将甄选后的岩石样品制备成与CT设备对应的尺寸。
具体的，根据不同测量范围的CT设备，制备与CT设备相应的样品尺寸。
作为可选方案，CT设备为微米CT设备。
具体的，岩油储层的孔隙尺寸适用范围为微米级以上，因此需选择对应的微米CT。
第二方面，本发明还提供一种页岩油储层中夹层孔隙有效三维表征的装置，包括：氦气法孔隙度测试仪，获取甄选后的岩石样品的氦气法实测孔隙度；CT设备，获取与CT设备对应的岩石样品的灰度图像；处理器，处理器与CT设备连接，处理器执行以下步骤：对灰度图像进行滤波降噪处理；设置图像处理的灰度阈值，对滤波降噪处理后的灰度图像进行分割和二值化处理，初步识别岩石样品的孔隙和骨架；基于初步识别的孔隙和骨架，获得岩石样品的计算孔隙度；基于计算孔隙度与实测孔隙度值的差值，以实测孔隙度值为标定，调整灰度阈值，获取最终识别的孔隙和骨架；基于最终识别的孔隙和骨架，获得岩石样品的孔隙度和孔径分布特征。
具体的，通过氦气法测定甄选后的岩石样品的孔隙度实测值，将甄选后的岩石样品制备成与CT设备对应的样品尺寸，CT设备获取对应岩石样品的灰度图像，基于Avizo软件对灰度图像进行滤波降噪处理后，设置图像处理的灰度阈值，以使得计算孔隙度近似等于实测孔隙度，然后对所述灰度图像进行阈值分割和二值化，识别出孔隙和骨架；基于孔隙和骨架，对识别出的孔隙进行数据统计，获得岩石样品的计算孔隙度，以实测孔隙度值作为标定，调整灰度阈值，重新识别孔隙和骨架，获得计算孔隙度，当计算孔隙度近似于实测孔隙度时，此时即确定识别出的孔隙、骨架为最终识别的结果。基于最终识别出的孔隙、骨架，建立岩石样品的三维数字模型，根据三维数字模型，获得岩石样品的孔隙度和孔径分布特征。
根据示例性的实施方式，页岩油储层中夹层孔隙有效三维表征的方法通过遴选均一化夹层岩石样品，以岩石样品的孔隙度实测值为标定，在图像处理中调整灰度阈值，进而使得计算获得的孔隙度近似等于孔隙度实测值，实现对页岩油储层非泥岩夹层孔隙的有效真实反应，有效地避免前人在应用CT技术时的局限性。
作为可选方案，基于计算孔隙度与实测孔隙度值的差值，以实测孔隙度值为标定，调整灰度阈值，获取最终识别的孔隙和骨架包括：计算计算孔隙度和实测孔隙度的差值，基于差值，以实测孔隙度值为标定，调整灰度阈值，获得新的计算孔隙度；重复此过程直至差值小于预设阈值时；在最后一次调整的灰度阈值下，获取最终识别的孔隙和骨架。
具体的，先对灰度图像进行滤波降噪处理，然后设置图像处理的灰度阈值，对图像进行阈值分割和二值化，识别出孔隙、骨架；对识别出的孔隙进行数据统计，计算输出岩石样品的孔隙度，以实测孔隙度值作为标定，调整灰度阈值，当计算孔隙度近似于实测孔隙度，此时即确定识别出的孔隙、骨架。
作为可选方案，基于最终识别的孔隙和骨架，获得岩石样品的孔隙度和孔径分布特征包括：基于最终识别的孔隙和骨架，建立岩石样品的三维数字模型；基于三维数字模型，获得岩石样品的孔隙度和孔径分布特征。
具体的，通过选用最终识别的孔隙和骨架的数百张图像，将这些二维CT图像进行堆叠排列并重构，形成三维数字模型，通过三维数字模型输出孔隙结构的三维特征；基于分水岭算法识别孔隙边界并对各孔隙进行标记，然后对各孔隙进行单独分析，通过统计各孔隙中的体素个数及边界上的像素个数来表征孔隙的体积，并按孔隙大小(即等效直径)分别对孔隙数量、体积的分布进行统计，计算输出孔隙度、孔径分布特征。
作为可选方案，将通过薄片鉴定或扫描电镜选取的均一性岩石样品作为甄选后的岩石样品。
具体的，中通过薄片鉴定、扫描电镜尽量遴选均一性岩石样品。
作为可选方案，根据CT设备的测量范围，将甄选后的岩石样品制备成与CT设备对应的尺寸。
具体的，根据不同测量范围的CT设备，制备与CT设备相应的样品尺寸。
作为可选方案，CT设备为微米CT设备。
具体的，岩油储层的孔隙尺寸适用范围为微米级以上，因此需选择对应的微米CT。
作为可选方案，氦气法孔隙度测试仪执行GB/T34533-2017标准获取甄选后的岩石样品的氦气法实测孔隙。
实施例一
图1示出了根据本发明的一个实施例的一种页岩油储层中夹层孔隙有效三维表征的方法的流程图。图2示出了根据本发明的一个实施例的一种页岩油储层中夹层孔隙有效三维表征的方法的孔隙结构三维特征。图3示出了根据本发明的一个实施例的一种页岩油储层中夹层孔隙有效三维表征的方法的孔径分布特征图。
结合图1、图2和图3所示，该页岩油储层中夹层孔隙有效三维表征的方法，包括：
步骤1：获取甄选后的岩石样品的实测孔隙度；
其中，将通过薄片鉴定或扫描电镜选取的均一性岩石样品作为甄选后的岩石样品。
其中，通过氦气法孔隙度测试仪获取甄选后的岩石样品的实测孔隙度。
步骤2：将甄选后的岩石样品制备成与CT设备对应的岩石样品；
其中，将甄选后的岩石样品制备成与CT设备对应的岩石样品包括：根据CT设备的测量范围，将甄选后的岩石样品制备成与CT设备对应的尺寸。
其中，CT设备为微米CT设备。
步骤3：获取与CT设备对应的岩石样品的灰度图像；
步骤4：对灰度图像进行滤波降噪处理；
步骤5：设置图像处理的灰度阈值，对滤波降噪处理后的灰度图像进行分割和二值化处理，初步识别岩石样品的孔隙和骨架；
步骤6：基于初步识别的孔隙和骨架，获得岩石样品的计算孔隙度；
步骤7：基于计算孔隙度与实测孔隙度值的差值，以实测孔隙度值为标定，调整灰度阈值，获取最终识别的孔隙和骨架；
其中，基于计算孔隙度与实测孔隙度值的差值，以实测孔隙度值为标定，调整灰度阈值，获取最终识别的孔隙和骨架包括：计算计算孔隙度和实测孔隙度的差值，基于差值，以实测孔隙度值为标定，调整灰度阈值，获得新的计算孔隙度；重复此过程直至差值小于预设阈值时；在最后一次调整的灰度阈值下，获取最终识别的孔隙和骨架。
步骤8：基于最终识别的孔隙和骨架，获得岩石样品的孔隙度和孔径分布特征。
其中，基于最终识别的孔隙和骨架，获得岩石样品的孔隙度和孔径分布特征包括：基于最终识别的孔隙和骨架，建立岩石样品的三维数字模型；基于三维数字模型，获得岩石样品的孔隙度和孔径分布特征。
具体的，通过选用最终识别的孔隙和骨架的数百张图像，将这些二维CT图像进行堆叠排列并重构，形成三维数字模型，通过三维数字模型输出孔隙结构的三维特征；基于分水岭算法识别孔隙边界并对各孔隙进行标记，然后对各孔隙进行单独分析，通过统计各孔隙中的体素个数及边界上的像素个数来表征孔隙的体积，并按孔隙大小(即等效直径)分别对孔隙数量、体积的分布进行统计，计算输出孔隙度、孔径分布特征。
以W99-A岩石样品(含泥白云岩)为例，孔隙空间为白云石晶间孔均匀分布。CT设备选用微米CT，该设备最大管电压为180kV，样品为1mm左右时最小的体元像素尺寸(分辨率)为0.5μm，能够分析的样品直径最大约为2.5cm，氦气法孔隙度测试仪的压力传感精度不低于±0.01psi，样品室体积20cm3-100cm3)、采用Avizo软件进行图像处理。
通过薄片、氩离子抛光将甄选后的岩石样品制备成直径2.5cm小岩心柱和直径2mm长4mm岩心柱。通过氦气法测得的孔隙度实测值为12.96％。
微米CT共计扫描1000张16bit灰度图片，总计492516bit灰度图片。图片为BMP格式，像素分辨率为1μm。Avizo软件对灰度图像滤波降噪处理后，调整灰度阈值，对灰度图像阈值分割和二值化，三维数字模型重构，输出孔隙结构的三维特征；对识别出的孔隙数据统计，从而计算输出岩石样品的孔隙度、主体孔径，计算孔隙度值近似于实测值(12.96％)，具体计算结果如表1所示。
表1W99-A孔隙灰度值调整前后孔径分布特征比较
实施例二
图4示出了根据本发明的一个实施例的一种页岩油储层中夹层孔隙有效三维表征的装置的框图。
如图4所示，该页岩油储层中夹层孔隙有效三维表征的装置，包括：
氦气法孔隙度测试仪102，获取甄选后的岩石样品的氦气法实测孔隙度；
CT设备104，获取与CT设备对应的岩石样品的灰度图像；
处理器106，处理器106与CT设备104连接，处理器106执行以下步骤：对灰度图像进行滤波降噪处理；设置图像处理的灰度阈值，对滤波降噪处理后的灰度图像进行分割和二值化处理，初步识别岩石样品的孔隙和骨架；基于初步识别的孔隙和骨架，获得岩石样品的计算孔隙度；基于计算孔隙度与实测孔隙度值的差值，以实测孔隙度值为标定，调整灰度阈值，获取最终识别的孔隙和骨架；基于最终识别的孔隙和骨架，获得岩石样品的孔隙度和孔径分布特征。
其中，基于计算孔隙度与实测孔隙度值的差值，以实测孔隙度值为标定，调整灰度阈值，获取最终识别的孔隙和骨架包括：计算计算孔隙度和实测孔隙度的差值，基于差值，以实测孔隙度值为标定，调整灰度阈值，获得新的计算孔隙度；重复此过程直至差值小于预设阈值时；在最后一次调整的灰度阈值下，获取最终识别的孔隙和骨架。
其中，基于最终识别的孔隙和骨架，获得岩石样品的孔隙度和孔径分布特征包括：基于最终识别的孔隙和骨架，建立岩石样品的三维数字模型；基于三维数字模型，获得岩石样品的孔隙度和孔径分布特征。
其中，将通过薄片鉴定或扫描电镜选取的均一性岩石样品作为甄选后的岩石样品。
其中，根据CT设备104的测量范围，将甄选后的岩石样品制备成与CT设备对应的尺寸。
其中，CT设备104为微米CT设备。
其中，氦气法孔隙度102测试仪执行GB/T34533-2017标准获取甄选后的岩石样品的氦气法实测孔隙。
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
Claims (10)
1.一种页岩油储层中夹层孔隙有效三维表征的方法，其特征在于，包括：
获取甄选后的岩石样品的实测孔隙度；
将所述甄选后的岩石样品制备成与CT设备对应的岩石样品；
获取与所述CT设备对应的岩石样品的灰度图像；
对所述灰度图像进行滤波降噪处理；
设置图像处理的灰度阈值，对所述滤波降噪处理后的灰度图像进行分割和二值化处理，初步识别所述岩石样品的孔隙和骨架；
基于初步识别的孔隙和骨架，获得所述岩石样品的计算孔隙度；
基于所述计算孔隙度与所述实测孔隙度值的差值，以所述实测孔隙度值为标定，调整灰度阈值，获取最终识别的孔隙和骨架；
基于所述最终识别的孔隙和骨架，获得所述岩石样品的孔隙度和孔径分布特征。
2.根据权利要求1所述的页岩油储层中夹层孔隙有效三维表征的方法，其特征在于，所述基于所述计算孔隙度与所述实测孔隙度值的差值，以所述实测孔隙度值为标定，调整灰度阈值，获取最终识别的孔隙和骨架包括：
计算所述计算孔隙度和实测孔隙度的差值，基于所述差值，以所述实测孔隙度值为标定，调整灰度阈值，获得新的计算孔隙度；重复此过程直至所述差值小于预设阈值时；
在最后一次调整的灰度阈值下，获取最终识别的孔隙和骨架。
3.根据权利要求2所述的页岩油储层中夹层孔隙有效三维表征的方法，其特征在于，所述基于所述最终识别的孔隙和骨架，获得所述岩石样品的孔隙度和孔径分布特征包括：
基于所述最终识别的孔隙和骨架，建立所述岩石样品的三维数字模型；
基于所述三维数字模型，获得所述岩石样品的孔隙度和孔径分布特征。
4.根据权利要求1所述的页岩油储层中夹层孔隙有效三维表征的方法，其特征在于，将通过薄片鉴定或扫描电镜选取的均一性岩石样品作为甄选后的岩石样品。
5.根据权利要求1所述的页岩油储层中夹层孔隙有效三维表征的方法，其特征在于，通过氦气法孔隙度测试仪获取甄选后的岩石样品的实测孔隙度。
6.根据权利要求1所述的页岩油储层中夹层孔隙有效三维表征的方法，其特征在于，所述将所述甄选后的岩石样品制备成与CT设备对应的岩石样品包括：
根据所述CT设备的测量范围，将所述甄选后的岩石样品制备成与所述CT设备对应的尺寸。
7.根据权利要求1所述的页岩油储层中夹层孔隙有效三维表征的方法，其特征在于，所述CT设备为微米CT设备。
8.一种页岩油储层中夹层孔隙有效三维表征的装置，其特征在于，包括：
氦气法孔隙度测试仪，获取甄选后的岩石样品的氦气法实测孔隙度；
CT设备，获取与所述CT设备对应的岩石样品的灰度图像；
处理器，所述处理器与所述CT设备连接，所述处理器执行以下步骤：
对所述灰度图像进行滤波降噪处理；
设置图像处理的灰度阈值，对所述滤波降噪处理后的灰度图像进行分割和二值化处理，初步识别所述岩石样品的孔隙和骨架；
基于初步识别的孔隙和骨架，获得所述岩石样品的计算孔隙度；
基于所述计算孔隙度与所述实测孔隙度值的差值，以所述实测孔隙度值为标定，调整灰度阈值，获取最终识别的孔隙和骨架；
基于所述最终识别的孔隙和骨架，获得所述岩石样品的孔隙度和孔径分布特征。
9.根据权利要求8所述的页岩油储层中夹层孔隙有效三维表征的装置，其特征在于，所述基于所述计算孔隙度与所述实测孔隙度值的差值，以所述实测孔隙度值为标定，调整灰度阈值，获取最终识别的孔隙和骨架包括：
计算所述计算孔隙度和实测孔隙度的差值，基于所述差值，以所述实测孔隙度值为标定，调整灰度阈值，获得新的计算孔隙度；
计算所述新的计算孔隙度与所述实测孔隙度值的差值，直至所述新的计算孔隙度与所述实测孔隙度值的差值小于预设阈值时，在最后一次调整的灰度阈值下，获取最终识别的孔隙和骨架。
10.根据权利要求8所述的页岩油储层中夹层孔隙有效三维表征的装置，其特征在于，将通过薄片鉴定或扫描电镜选取的均一性岩石样品作为甄选后的岩石样品。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011314324.4A CN114518309A (zh) | 2020-11-20 | 2020-11-20 | 一种页岩油储层中夹层孔隙有效三维表征的方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011314324.4A CN114518309A (zh) | 2020-11-20 | 2020-11-20 | 一种页岩油储层中夹层孔隙有效三维表征的方法及装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114518309A true CN114518309A (zh) | 2022-05-20 |
Family
ID=81594213
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011314324.4A Pending CN114518309A (zh) | 2020-11-20 | 2020-11-20 | 一种页岩油储层中夹层孔隙有效三维表征的方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114518309A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2024108689A1 (zh) * | 2022-11-21 | 2024-05-30 | 中国石油天然气股份有限公司 | 一种花斑状碳酸盐岩的储层孔隙度和/或储层渗透率确定方法 |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103513285A (zh) * | 2013-09-27 | 2014-01-15 | 中国石油天然气股份有限公司 | 一种确定横向表面弛豫速率的方法及装置 |
| CN105205296A (zh) * | 2014-06-09 | 2015-12-30 | 中国石油化工股份有限公司 | 一种求取页岩气储层孔隙度的方法 |
| CN105372168A (zh) * | 2015-12-08 | 2016-03-02 | 清华大学 | 一种混凝土孔隙三维结构的构建方法 |
| CN106546521A (zh) * | 2016-10-12 | 2017-03-29 | 北京师范大学 | 一种基于ct扫描技术量化土壤大孔隙空间网络结构的方法 |
| CN108444881A (zh) * | 2016-11-03 | 2018-08-24 | 中国石油化工股份有限公司 | 一种适用于陆相泥页岩微纳米尺度储集空间的表征方法 |
| CN109283116A (zh) * | 2018-11-14 | 2019-01-29 | 中国地质大学（武汉） | 一种页岩储层孔隙-裂隙表征方法和装置 |
| CN109285222A (zh) * | 2018-09-06 | 2019-01-29 | 中国地质大学（北京） | 有机页岩高分辨率数字岩心构建与分析方法 |
| CN109697752A (zh) * | 2018-11-20 | 2019-04-30 | 中国石油天然气集团有限公司 | 基于岩心ct图像孔隙信息提取定量表征岩心非均质性方法 |
| CN111175207A (zh) * | 2018-11-12 | 2020-05-19 | 中石化石油工程技术服务有限公司 | 一种岩石有效孔隙度的检测方法 |
-
2020
- 2020-11-20 CN CN202011314324.4A patent/CN114518309A/zh active Pending
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103513285A (zh) * | 2013-09-27 | 2014-01-15 | 中国石油天然气股份有限公司 | 一种确定横向表面弛豫速率的方法及装置 |
| CN105205296A (zh) * | 2014-06-09 | 2015-12-30 | 中国石油化工股份有限公司 | 一种求取页岩气储层孔隙度的方法 |
| CN105372168A (zh) * | 2015-12-08 | 2016-03-02 | 清华大学 | 一种混凝土孔隙三维结构的构建方法 |
| CN106546521A (zh) * | 2016-10-12 | 2017-03-29 | 北京师范大学 | 一种基于ct扫描技术量化土壤大孔隙空间网络结构的方法 |
| CN108444881A (zh) * | 2016-11-03 | 2018-08-24 | 中国石油化工股份有限公司 | 一种适用于陆相泥页岩微纳米尺度储集空间的表征方法 |
| CN109285222A (zh) * | 2018-09-06 | 2019-01-29 | 中国地质大学（北京） | 有机页岩高分辨率数字岩心构建与分析方法 |
| CN111175207A (zh) * | 2018-11-12 | 2020-05-19 | 中石化石油工程技术服务有限公司 | 一种岩石有效孔隙度的检测方法 |
| CN109283116A (zh) * | 2018-11-14 | 2019-01-29 | 中国地质大学（武汉） | 一种页岩储层孔隙-裂隙表征方法和装置 |
| CN109697752A (zh) * | 2018-11-20 | 2019-04-30 | 中国石油天然气集团有限公司 | 基于岩心ct图像孔隙信息提取定量表征岩心非均质性方法 |
Non-Patent Citations (2)
| Title |
|---|
| 张文涛，等: "《基于流体吸入实验的页岩纳米孔隙连通性分析方法》", 《石油实验地质》, vol. 42, no. 3 * |
| 赵秀才;姚军;房克荣;: "合理分割岩心微观结构图像的新方法", 中国石油大学学报(自然科学版), no. 01 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2024108689A1 (zh) * | 2022-11-21 | 2024-05-30 | 中国石油天然气股份有限公司 | 一种花斑状碳酸盐岩的储层孔隙度和/或储层渗透率确定方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| RU2543698C1 (ru) | Анализ петрографических изображений для определения капиллярного давления в пористых средах | |
| CN105115874B (zh) | 多源信息融合的多组分三维数字岩心构建方法 | |
| Liu et al. | Pore-scale characterization of tight sandstone in Yanchang Formation Ordos Basin China using micro-CT and SEM imaging from nm-to cm-scale | |
| Ramandi et al. | Micro-CT image calibration to improve fracture aperture measurement | |
| CN110415167B (zh) | 一种基于数字图像技术的粗糙面裂隙生成方法及试验系统 | |
| Cnudde et al. | 3D characterization of sandstone by means of X-ray computed tomography | |
| CN103513285B (zh) | 一种确定横向表面弛豫速率的方法及装置 | |
| Ketcham et al. | Accurate measurement of small features in X‐ray CT data volumes, demonstrated using gold grains | |
| Jasti et al. | Microscopic imaging of porous media with X-ray computer tomography | |
| CN109490165A (zh) | 表征碳酸盐岩非组构选择性储集空间的方法 | |
| CN108335310B (zh) | 一种便携式粒形粒度检测方法及系统 | |
| Fogden* et al. | Micro-CT analysis of pores and organics in unconventionals using novel contrast strategies | |
| CN111624146A (zh) | 一种定量表征致密储层全尺度孔径分布特征的方法 | |
| CN112686917A (zh) | 提高岩心非均质性表征精度的数字岩心建模方法及装置 | |
| CN115147340A (zh) | 岩石裂缝变形信息获取方法、装置、设备和存储介质 | |
| CN114518309A (zh) | 一种页岩油储层中夹层孔隙有效三维表征的方法及装置 | |
| CN113129275B (zh) | 一种基于岩土体材料数字图像三维结构表征方法 | |
| CN112147680A (zh) | 岩石内部微裂缝提取和刻画方法及系统 | |
| CN110672496A (zh) | 孔隙测量方法、装置、设备及存储介质 | |
| CN112435288B (zh) | 一种基于图像的孔隙特征计算方法 | |
| Fogden et al. | 3D mapping of pore and organic matter distributions in unconventional reservoirs utilizing a digital rocks approach | |
| Sheppard et al. | Quantitative properties of complex porous materials calculated from x-ray μCT images | |
| RU2621371C1 (ru) | Способ исследования фильтрационно-емкостных свойств горных пород | |
| CN116413183A (zh) | 一种孔隙连通性评价方法和装置 | |
| CN116168073A (zh) | 一种数字岩心渗流特性分析的rev确定方法和装置 |
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 |