CN104751473B - 碳酸岩多尺度孔洞特征的确定方法及装置 - Google Patents
碳酸岩多尺度孔洞特征的确定方法及装置 Download PDFInfo
- Publication number
- CN104751473B CN104751473B CN201510169106.9A CN201510169106A CN104751473B CN 104751473 B CN104751473 B CN 104751473B CN 201510169106 A CN201510169106 A CN 201510169106A CN 104751473 B CN104751473 B CN 104751473B
- Authority
- CN
- China
- Prior art keywords
- hole
- scan image
- sample
- carbonatite
- 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
Links
- 239000010430 carbonatite Substances 0.000 title claims abstract description 107
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000006243 chemical reaction Methods 0.000 claims abstract description 45
- 238000002591 computed tomography Methods 0.000 claims abstract description 33
- 238000010586 diagram Methods 0.000 claims description 61
- 239000011435 rock Substances 0.000 claims description 17
- 239000002245 particle Substances 0.000 claims description 13
- 238000012545 processing Methods 0.000 claims description 11
- 230000015572 biosynthetic process Effects 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 6
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 4
- 229910052799 carbon Inorganic materials 0.000 claims description 4
- 230000004323 axial length Effects 0.000 claims 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 13
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 11
- 229910052753 mercury Inorganic materials 0.000 description 11
- 239000011148 porous material Substances 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 7
- 238000005266 casting Methods 0.000 description 7
- 238000002347 injection Methods 0.000 description 6
- 239000007924 injection Substances 0.000 description 6
- BVKZGUZCCUSVTD-UHFFFAOYSA-N carbonic acid Chemical compound OC(O)=O BVKZGUZCCUSVTD-UHFFFAOYSA-N 0.000 description 5
- 230000035699 permeability Effects 0.000 description 4
- 238000006073 displacement reaction Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 239000004744 fabric Substances 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 229920005479 Lucite® Polymers 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 239000002734 clay mineral Substances 0.000 description 2
- 230000006378 damage Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 239000003822 epoxy resin Substances 0.000 description 2
- 238000002493 microarray Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 229920000647 polyepoxide Polymers 0.000 description 2
- 239000004926 polymethyl methacrylate Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- KYKQHSMYWLWROM-UHFFFAOYSA-N ac1l4yjn Chemical compound [Hg].[Hg] KYKQHSMYWLWROM-UHFFFAOYSA-N 0.000 description 1
- 239000002253 acid Substances 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002389 environmental scanning electron microscopy Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000007747 plating Methods 0.000 description 1
- 231100000614 poison Toxicity 0.000 description 1
- 230000007096 poisonous effect Effects 0.000 description 1
- 230000037452 priming Effects 0.000 description 1
- 238000012797 qualification Methods 0.000 description 1
- 238000004626 scanning electron microscopy Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000005406 washing Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明实施例提供了一种碳酸岩多尺度孔洞特征的确定方法及装置,其中,该方法包括:对直径大于全直径的碳酸岩样品进行毫米CT扫描,获得第一扫描图像;在碳酸岩样品上取直径为微米级的孔隙或孔洞分布密集的部分作为第一柱塞样品,对该第一柱塞样品进行微米CT扫描,获得第二扫描图像;在第一柱塞样品上取直径为纳米级的孔隙或孔洞分布密集的部分作为第二柱塞样品,对第二柱塞样品进行纳米CT扫描,获得第三扫描图像;根据二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像,确定碳酸岩样品的孔隙度以及不同直径的孔隙和孔洞在所述碳酸岩样品上的分布方式。该方案可以实现对碳酸岩多尺度孔洞特征进行定量、精细、可靠地分析。
Description
技术领域
本发明涉及碳酸岩储层微观结构研究技术领域,特别涉及一种碳酸岩多尺度孔洞特征的确定方法及装置。
背景技术
实践表明碳酸盐岩储层的孔洞、裂缝存在多尺度特性,大的洞穴可到几厘米,甚至更大,小的孔隙可到几个微米,甚至更小,表现出很强的非均质性。目前,对碳酸盐岩储层的孔、洞、缝微观结构特征研究的方法和认识较少,主要采用常规的铸体薄片、扫描电镜和压汞等手段,但是,这些手段都存在一些缺陷,使得不能定量、精细地对碳酸盐岩储层的孔、洞、缝微观结构进行可靠分析。
例如,基于铸体薄片和扫描电镜技术手段的储层研究中的微观分析技术。孔隙铸体薄片是将地下含油水岩石经洗油烘干,放到铸体仪中将带色的有机玻璃或环氧树脂通过真空灌注到岩石的孔隙裂缝中,再经高温高压固化、酸化去除岩石颗粒后,只保留孔隙骨架,对孔隙骨架表面镀模处理后放在显微镜下观察,可以直接观察岩心薄片中的面孔率、孔隙、喉道、孔喉配位数以及碎屑组分等微观形貌信息。此方法简单易行,成本低廉,是目前储层孔隙结构研究中的常用方法。但是,该分析方法存在以下缺陷:
1)利用铸体薄片图像求取储层孔隙度,需要专业技术人士去判读,对操作人员的技术水平经验依赖性大,铸体薄片的厚薄、平整度的差别会直接影响鉴定结果的好坏;而实际工作中,因操作人员、摄像参数设置等的不同,分析结果可能存在差别;
2)识别图像代表二维局部孔隙,代表性差。由于铸体薄片尺度太小,放大后的视域很小,而岩石的非均质性很强,因此,基于很小的视域很难对碳酸盐岩的微观结构做出评价;
3)铸体薄片在向岩石的孔隙裂缝灌注带色的有机玻璃或环氧树脂的过程中,对原来的孔隙及其内部的粘土矿物等会产生破坏,产生人工诱导孔隙和裂缝,使得形成的孔隙骨架反映不出真实的碳酸盐岩孔隙结构特征。
例如,基于压汞法的手段的储层研究中的微观分析技术。常规压汞法以毛细管束模型为基础,假设多孔介质由直径大小不同的毛细管束组成。将非润湿相的汞注入被抽空的岩样空间中去时,一定要克服岩石孔隙***对汞的毛细管阻力。汞作为非润湿相,汞的注入过程中汞首先进入较大的孔隙喉道中;随着注入压力不断增加,汞即不断进入较小的孔隙。岩样的毛细管压力与孔径分布、孔隙体积分布、孔喉连通关系、孔隙度、渗透率及流体饱和度等都有关。因此,利用毛管曲线可以分析不同渗透率岩心的孔喉分布及渗透率贡献率分布特征,探讨渗透率、平均喉道半径、分选系数等孔隙结构参数。但是其也存在以下缺点:
1)压汞实验过程中,使用的汞有毒,会对人体和环境造成危害;
2)压汞实验时,汞的注入过程中会对原来的孔隙及其内部的粘土矿物等会产生破坏,使得反映不出真实的岩石孔隙结构特征,而且注入压力越高,这种影响越大。
发明内容
本发明实施例提供了一种碳酸岩多尺度孔洞特征的确定方法,以实现对碳酸岩多尺度孔洞特征进行定量、精细、可靠地分析。该方法包括:对直径大于全直径的碳酸岩样品进行毫米CT扫描,获得第一扫描图像;在所述碳酸岩样品上取直径为微米级的孔隙或孔洞分布密集的部分作为第一柱塞样品,对该第一柱塞样品进行微米CT扫描,获得第二扫描图像;在所述第一柱塞样品上取直径为纳米级的孔隙或孔洞分布密集的部分作为第二柱塞样品,对该第二柱塞样品进行纳米CT扫描,获得第三扫描图像;对所述第一扫描图像、所述第二扫描图像和所述第三扫描图像进行二值化处理,根据二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像,确定所述碳酸岩样品的孔隙度以及不同直径的孔隙和孔洞在所述碳酸岩样品上的分布方式,其中,在二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像中,孔隙和孔洞的像素数相同,孔隙和孔洞的像素数与所述碳酸岩样品的颗粒的像素数不同。
在一个实施例中,在所述第一扫描图像中,孔隙或孔洞的直径大于等于0.5毫米,在所述第二扫描图像中,孔隙或孔洞的直径大于等于0.2微米小于0.5毫米,在所述第三扫描图像中,孔隙或孔洞的直径大于等于50纳米小于0.2微米。
在一个实施例中,根据二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像,确定所述碳酸岩样品的孔隙度,包括:根据孔隙的像素数分别计算二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像中孔隙的面积,将二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像中孔隙的面积相加得到所有孔隙的总面积;根据二值化处理后的第一扫描图像中孔隙的像素数和颗粒的像素数计算所述碳酸岩样品的CT图像面积;将所述所有孔隙的总面积与碳酸岩样品的CT图像面积的比值,确定为所述碳酸岩样品的孔隙度。
在一个实施例中,根据二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像,确定不同直径的孔隙和孔洞在所述碳酸岩样品上的分布方式,包括:将二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像中孔洞和孔隙的图形处理为轮廓图,其中,轮廓图是孔洞或孔隙的图形的轮廓线形成的封闭的图形;根据轮廓图的分布方式确定不同直径的孔隙和孔洞在所述碳酸盐样品上的分布方式。
在一个实施例中,根据轮廓图的分布方式确定不同直径的孔隙和孔洞在所述碳酸岩样品上的分布方式,包括:根据孔洞和孔隙的轮廓图,获取孔洞和孔隙的骨架图,骨架图是由线段组成的表示孔洞或孔隙的拓扑结构的图形;根据骨架图在碳酸岩样品的CT图像上的坐标,确定孔洞或孔隙在所述碳酸岩样品上的位置。
在一个实施例中,还包括:对于每个轮廓图,将轮廓图划分为多个图块,在每个图块中建立与图块边缘相切的椭圆,多个图块中椭圆的短轴方向相同;将每个图块短轴的长度与该图块面积相乘,得到每个图块的乘积,将多个图块的乘积相加,得到的和为与该轮廓图边缘相切的椭圆的短轴长度,确定为孔洞或孔隙的短轴长度。
本发明实施例还提供了一种碳酸岩多尺度孔洞特征的确定装置,以实现对碳酸岩多尺度孔洞特征进行定量、精细、可靠地分析。该装置包括:第一图像扫描模块,用于对直径大于全直径的碳酸岩样品进行毫米CT扫描,获得第一扫描图像;第二图像扫描模块,用于在所述碳酸岩样品上取直径为微米级的孔隙或孔洞分布密集的部分作为第一柱塞样品,对该第一柱塞样品进行微米CT扫描,获得第二扫描图像;第三图像扫描模块,用于在所述第一柱塞样品上取直径为纳米级的孔隙或孔洞分布密集的部分作为第二柱塞样品,对该第二柱塞样品进行纳米CT扫描,获得第三扫描图像;确定模块,用于对所述第一扫描图像、所述第二扫描图像和所述第三扫描图像进行二值化处理,根据二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像,确定所述碳酸岩样品的孔隙度以及不同直径的孔隙和孔洞在所述碳酸岩样品上的分布方式,其中,在二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像中,孔隙和孔洞的像素数相同,孔隙和孔洞的像素数与所述碳酸岩样品的颗粒的像素数不同。
在一个实施例中,在所述第一扫描图像中,孔隙或孔洞的直径大于等于0.5毫米,在所述第二扫描图像中,孔隙或孔洞的直径大于等于0.2微米小于0.5毫米,在所述第三扫描图像中,孔隙或孔洞的直径大于等于50纳米小于0.2微米。
在一个实施例中,所述确定模块,包括:第一面积计算单元,用于根据孔隙的像素数分别计算二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像中孔隙的面积,将二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像中孔隙的面积相加得到所有孔隙的总面积;第二面积计算单元,用于根据二值化处理后的第一扫描图像中孔隙的像素数和颗粒的像素数计算所述碳酸岩样品的CT图像面积;孔隙度确定单元,用于将所述所有孔隙的总面积与碳酸岩样品的CT图像面积的比值,确定为所述碳酸岩样品的孔隙度。
在一个实施例中,所述确定模块,还包括:轮廓图处理单元,用于将二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像中孔洞和孔隙的图形处理为轮廓图,其中,轮廓图是孔洞或孔隙的图形的轮廓线形成的封闭的图形;分布确定单元,用于根据轮廓图的分布方式确定不同直径的孔隙和孔洞在所述碳酸盐样品上的分布方式。
在一个实施例中,所述确定模块,还包括:骨架图处理单元,用于根据孔洞和孔隙的轮廓图,获取孔洞和孔隙的骨架图,骨架图是由线段组成的表示孔洞或孔隙的拓扑结构的图形;位置确定单元,用于根据骨架图在碳酸岩样品的CT图像上的坐标,确定孔洞或孔隙在所述碳酸岩样品上的位置。
在一个实施例中,还包括:图块划分模块,用于对于每个轮廓图,将轮廓图划分为多个图块,在每个图块中建立与图块边缘相切的椭圆,多个图块中椭圆的短轴方向相同;轴长确定模块,用于将每个图块短轴的长度与该图块面积相乘,得到每个图块的乘积,将多个图块的乘积相加,得到的和为与该轮廓图边缘相切的椭圆的短轴长度,确定为孔洞或孔隙的短轴长度。
在本发明实施例中,通过对碳酸岩样品进行毫米CT扫描,在获得的第一扫描图像中可以显示直径为毫米级的孔隙和孔洞的分布情况,在进一步对直径为微米级的孔隙或孔洞分布密集的样品部分进行微米CT扫描,在获得的第二扫描图像中可以显示直径为微米级的孔隙和孔洞的分布情况,为了更精细地确定孔隙和孔洞的分布情况,对直径为纳米级的孔隙或孔洞分布密集的样品部分进行纳米CT扫描,在获得的第三扫描图像中可以显示直径为纳米级的孔隙和孔洞的分布情况,最终将获得的第一扫描图像、第二扫描图像和第三扫描图像综合起来,可以精确、定量得分析不同直径尺度的孔隙和孔洞在碳酸岩样品上的分布,进而可以实现对碳酸岩多尺度孔洞特征进行定量、精细、可靠地分析,使得有助于为碳酸盐岩储层的储层评价、地质建模,油藏的有效驱替和调堵、压裂酸化改造等提供明确的各种尺度的孔洞分布特征及其连通关系和改造效果的评价,为油田开发与增产提供有力的基础和清晰地认识。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,并不构成对本发明的限定。在附图中:
图1是本发明实施例提供的一种碳酸岩多尺度孔洞特征的确定方法的流程图;
图2是本发明实施例提供的一种孔洞和孔隙的轮廓图的示意图;
图3是本发明实施例提供的一种孔洞和孔隙的骨架图的示意图;
图4是本发明实施例提供的一种孔洞和孔隙的轮廓图和骨架图重叠的示意图;
图5是本发明实施例提供的一种毫米级孔洞和孔隙分布情况的示意图;
图6是本发明实施例提供的一种碳酸岩多尺度孔洞特征的确定装置的结构框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施方式和附图,对本发明做进一步详细说明。在此,本发明的示意性实施方式及其说明用于解释本发明,但并不作为对本发明的限定。
在本发明实施例中,提供了一种碳酸岩多尺度孔洞特征的确定方法,如图1所示,该方法包括:
步骤101:对直径大于全直径的碳酸岩样品进行毫米CT扫描,获得第一扫描图像;
步骤102:在所述碳酸岩样品上取直径为微米级的孔隙或孔洞分布密集的部分作为第一柱塞样品,对该第一柱塞样品进行微米CT扫描,获得第二扫描图像;
步骤103:在所述第一柱塞样品上取直径为纳米级的孔隙或孔洞分布密集的部分作为第二柱塞样品,对该第二柱塞样品进行纳米CT扫描,获得第三扫描图像;
步骤104:对所述第一扫描图像、所述第二扫描图像和所述第三扫描图像进行二值化处理,根据二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像,确定所述碳酸岩样品的孔隙度以及不同直径的孔隙和孔洞在所述碳酸岩样品上的分布方式,其中,在二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像中,孔隙和孔洞的像素数相同,孔隙和孔洞的像素数与所述碳酸岩样品的颗粒的像素数不同。
由图1所示的流程可知,在本发明实施例中,通过对碳酸岩样品进行毫米CT扫描,在获得的第一扫描图像中可以显示直径为毫米级的孔隙和孔洞的分布情况,在进一步对直径为微米级的孔隙或孔洞分布密集的样品部分进行微米CT扫描,在获得的第二扫描图像中可以显示直径为微米级的孔隙和孔洞的分布情况,为了更精细地确定孔隙和孔洞的分布情况,对直径为纳米级的孔隙或孔洞分布密集的样品部分进行纳米CT扫描,在获得的第三扫描图像中可以显示直径为纳米级的孔隙和孔洞的分布情况,最终将获得的第一扫描图像、第二扫描图像和第三扫描图像综合起来,可以精确、定量得分析不同直径尺度的孔隙和孔洞在碳酸岩样品上的分布,进而可以实现对碳酸岩多尺度孔洞特征进行定量、精细、可靠地分析,使得有助于为碳酸盐岩储层的储层评价、地质建模,油藏的有效驱替和调堵、压裂酸化改造等提供明确的各种尺度的孔洞分布特征及其连通关系和改造效果的评价,为油田开发与增产提供有力的基础和清晰地认识。
具体实施时,上述碳酸岩样品可以是碳酸盐岩样品、碳酸砂岩样品或者是碳酸页岩样品。
具体实施时,为了可以精细地分析不同直径的孔隙和孔洞的分布情况,在本实施例中,在对碳酸岩样品进行毫米CT扫描时,可以采用高能毫米CT仪器,其最小分辨率为0.05毫米,这样在获得的第一扫描图像中可以识别毫米级的孔洞和孔隙,例如,在第一扫描图像中可以识别直径大于等于0.5毫米的孔隙和孔洞;为了进一步观察微米级直径的孔隙和孔洞,在碳酸岩样品上取直径为微米级的孔隙或孔洞分布密集的部分作为第一柱塞样品,该第一柱塞样品的直径可以为φ25毫米,对第一柱塞样品进行微米CT扫描时,以采用微米CT仪器,其最小分辨率为0.5微米,这样在获得的第二扫描图像中可以识别微米级的孔洞和孔隙,例如,在第二扫描图像中可以识别直径大于等于0.2微米小于0.5毫米的孔隙和孔洞;类似的,在第一柱塞样品上取直径为纳米级的孔隙或孔洞分布密集的部分作为第二柱塞样品,进行微米CT扫描,可以采用纳米CT仪器进行扫描,在获得的第三扫描图像中可以识别直径大于等于50纳米小于0.2微米的孔隙和孔洞。
具体实施时,可以基于二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像,精确地确定所述碳酸岩样品的孔隙度,例如,根据孔隙的像素数分别计算二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像中孔隙的面积,将二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像中孔隙的面积相加得到所有孔隙的总面积;然后,根据二值化处理后的第一扫描图像中孔隙的像素数和颗粒的像素数计算所述碳酸岩样品的CT图像面积,最后,将所有孔隙的总面积与碳酸岩样品的CT图像面积的比值,确定为所述碳酸岩样品的孔隙度。
具体实施时,还可以基于二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像,确定不同直径的孔隙和孔洞在所述碳酸岩样品上的分布方式。例如,将二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像中孔洞和孔隙的图形处理为轮廓图,如图2所示(如图2中所示的不规则图形即为孔洞或孔隙的轮廓图),轮廓图是孔洞或孔隙的图形的轮廓线形成的封闭的图形;根据轮廓图的分布方式确定不同直径的孔隙和孔洞在所述碳酸岩样品上的分布方式。
具体实施时,为了确定孔洞和孔隙在碳酸岩样品上的具***置,在本实施例中,根据孔洞和孔隙的轮廓图,获取孔洞和孔隙的骨架图(如图3、4所示),骨架图是由线段组成的表示孔洞或孔隙的拓扑结构的图形;根据骨架图在碳酸岩样品的CT图像上的坐标,确定孔洞或孔隙在碳酸岩样品上的位置。具体的,在扫描碳酸岩样品获得CT图像时,建立碳酸岩样品上的位置与CT图像坐标的对应关系,根据骨架图在碳酸岩样品的CT图像上的坐标,确定孔洞或孔隙在碳酸岩样品上的位置。
具体实施时,为了定量地确定不同尺度的孔洞和孔隙的大小,在本实施例中,对于每个轮廓图,将轮廓图划分为多个图块,在每个图块中建立与图块边缘相切的椭圆形,多个图块中椭圆的短轴方向相同;将每个图块短轴的长度与该图块面积相乘,得到每个图块的乘积,将多个图块的乘积相加,得到的和为与该轮廓图边缘相切的椭圆的短轴长度,确定为孔洞或孔隙的短轴长度。具体的,确定了孔洞或孔隙的短轴长度后,即可以确定孔洞或孔隙的直径大小,进而可以统计出第一扫描图像、第二扫描图像以及第三扫描图像中不同直径孔洞或孔隙的分布率,例如,如图5所示的在第一扫描图像中毫米级的不同直径孔洞或孔隙的分布率。
具体的,还可以定量地确定孔洞和孔隙的形态,在本实施例中,将每个图块长轴的长度与该图块面积相乘,得到每个图块的乘积,将多个图块的乘积相加,得到的和为与该轮廓图边缘相切的椭圆的长轴长度;根据轮廓图的长轴长度与短轴长度的大小关系,确定孔洞或孔隙的图形形态。例如,当轮廓图的长轴长度大于短轴长度,则轮廓图的图形为椭圆,轮廓图对应的孔洞或孔隙的形态为椭圆;当轮廓图的长轴长度等于短轴长度,则轮廓图的图形为圆形,轮廓图对应的孔洞或孔隙的形态为圆形。
基于同一发明构思,本发明实施例中还提供了一种碳酸岩多尺度孔洞特征的确定装置,如下面的实施例所述。由于碳酸岩多尺度孔洞特征的确定装置解决问题的原理与碳酸岩多尺度孔洞特征的确定方法相似,因此碳酸岩多尺度孔洞特征的确定装置的实施可以参见碳酸岩多尺度孔洞特征的确定方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
图6是本发明实施例的碳酸岩多尺度孔洞特征的确定装置的一种结构框图,如图6所示,包括:第一图像扫描模块601、第二图像扫描模块602、第三图像扫描模块603和确定模块604,下面对该结构进行说明。
第一图像扫描模块601,用于对直径大于全直径的碳酸岩样品进行毫米CT扫描,获得第一扫描图像;
第二图像扫描模块602,与第一图像扫描模块601连接,用于在所述碳酸岩样品上取直径为微米级的孔隙或孔洞分布密集的部分作为第一柱塞样品,对该第一柱塞样品进行微米CT扫描,获得第二扫描图像;
第三图像扫描模块603,与第二图像扫描模块602连接,用于在所述第一柱塞样品上取直径为纳米级的孔隙或孔洞分布密集的部分作为第二柱塞样品,对该第二柱塞样品进行纳米CT扫描,获得第三扫描图像;
确定模块604,与第三图像扫描模块603连接,用于对所述第一扫描图像、所述第二扫描图像和所述第三扫描图像进行二值化处理,根据二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像,确定所述碳酸岩样品的孔隙度以及不同直径的孔隙和孔洞在所述碳酸岩样品上的分布方式,其中,在二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像中,孔隙和孔洞的像素数相同,孔隙和孔洞的像素数与所述碳酸岩样品的颗粒的像素数不同。
在一个实施例中,在所述第一扫描图像中,孔隙或孔洞的直径大于等于0.5毫米,在所述第二扫描图像中,孔隙或孔洞的直径大于等于0.2微米小于0.5毫米,在所述第三扫描图像中,孔隙或孔洞的直径大于等于50纳米小于0.2微米。
在一个实施例中,所述确定模块,包括:第一面积计算单元,用于根据孔隙的像素数分别计算二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像中孔隙的面积,将二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像中孔隙的面积相加得到所有孔隙的总面积;第二面积计算单元,用于根据二值化处理后的第一扫描图像中孔隙的像素数和颗粒的像素数计算所述碳酸岩样品的CT图像面积;孔隙度确定单元,与第二面积计算单元和第一面积计算单元连接,用于将所述所有孔隙的总面积与碳酸岩样品的CT图像面积的比值,确定为所述碳酸岩样品的孔隙度。
在一个实施例中,所述确定模块,还包括:轮廓图处理单元,用于将二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像中孔洞和孔隙的图形处理为轮廓图,其中,轮廓图是孔洞或孔隙的图形的轮廓线形成的封闭的图形;分布确定单元,与轮廓图处理单元连接,用于根据轮廓图的分布方式确定不同直径的孔隙和孔洞在所述碳酸盐样品上的分布方式。
在一个实施例中,所述确定模块,还包括:骨架图处理单元,用于根据孔洞和孔隙的轮廓图,获取孔洞和孔隙的骨架图,骨架图是由线段组成的表示孔洞或孔隙的拓扑结构的图形;位置确定单元,与骨架图处理单元连接,用于根据骨架图在碳酸岩样品的CT图像上的坐标,确定孔洞或孔隙在所述碳酸岩样品上的位置。
在一个实施例中,还包括:图块划分模块,用于对于每个轮廓图,将轮廓图划分为多个图块,在每个图块中建立与图块边缘相切的椭圆,多个图块中椭圆的短轴方向相同;轴长确定模块,与图块划分模块连接,用于将每个图块短轴的长度与该图块面积相乘,得到每个图块的乘积,将多个图块的乘积相加,得到的和为与该轮廓图边缘相切的椭圆的短轴长度,确定为孔洞或孔隙的短轴长度。
在本发明实施例中,通过对碳酸岩样品进行毫米CT扫描,在获得的第一扫描图像中可以显示直径为毫米级的孔隙和孔洞的分布情况,在进一步对直径为微米级的孔隙或孔洞分布密集的样品部分进行微米CT扫描,在获得的第二扫描图像中可以显示直径为微米级的孔隙和孔洞的分布情况,为了更精细地确定孔隙和孔洞的分布情况,对直径为纳米级的孔隙或孔洞分布密集的样品部分进行纳米CT扫描,在获得的第三扫描图像中可以显示直径为纳米级的孔隙和孔洞的分布情况,最终将获得的第一扫描图像、第二扫描图像和第三扫描图像综合起来,可以精确、定量得分析不同直径尺度的孔隙和孔洞在碳酸岩样品上的分布,进而可以实现对碳酸岩多尺度孔洞特征进行定量、精细、可靠地分析,使得有助于为碳酸盐岩储层的储层评价、地质建模,油藏的有效驱替和调堵、压裂酸化改造等提供明确的各种尺度的孔洞分布特征及其连通关系和改造效果的评价,为油田开发与增产提供有力的基础和清晰地认识。
显然,本领域的技术人员应该明白,上述的本发明实施例的各模块或各步骤可以用通用的计算装置来实现,它们可以集中在单个的计算装置上,或者分布在多个计算装置所组成的网络上,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,并且在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样,本发明实施例不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种碳酸岩多尺度孔洞特征的确定方法,其特征在于,包括:
对直径大于全直径的碳酸岩样品进行毫米CT扫描,获得第一扫描图像;
在所述碳酸岩样品上取直径为微米级的孔隙或孔洞分布密集的部分作为第一柱塞样品,对该第一柱塞样品进行微米CT扫描,获得第二扫描图像;
在所述第一柱塞样品上取直径为纳米级的孔隙或孔洞分布密集的部分作为第二柱塞样品,对该第二柱塞样品进行纳米CT扫描,获得第三扫描图像;
对所述第一扫描图像、所述第二扫描图像和所述第三扫描图像进行二值化处理,根据二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像,确定所述碳酸岩样品的孔隙度以及不同直径的孔隙和孔洞在所述碳酸岩样品上的分布方式,其中,在二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像中,孔隙和孔洞的像素数相同,孔隙和孔洞的像素数与所述碳酸岩样品的颗粒的像素数不同;
在所述第一扫描图像中,孔隙或孔洞的直径大于等于0.5毫米,在所述第二扫描图像中,孔隙或孔洞的直径大于等于0.2微米小于0.5毫米,在所述第三扫描图像中,孔隙或孔洞的直径大于等于50纳米小于0.2微米。
2.如权利要求1所述的方法,其特征在于,根据二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像,确定所述碳酸岩样品的孔隙度,包括:
根据孔隙的像素数分别计算二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像中孔隙的面积,将二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像中孔隙的面积相加得到所有孔隙的总面积;
根据二值化处理后的第一扫描图像中孔隙的像素数和颗粒的像素数计算所述碳酸岩样品的CT图像面积;
将所述所有孔隙的总面积与碳酸岩样品的CT图像面积的比值,确定为所述碳酸岩样品的孔隙度。
3.如权利要求1所述的方法,其特征在于,根据二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像,确定不同直径的孔隙和孔洞在所述碳酸岩样品上的分布方式,包括:
将二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像中孔洞和孔隙的 图形处理为轮廓图,其中,轮廓图是孔洞或孔隙的图形的轮廓线形成的封闭的图形;
根据轮廓图的分布方式确定不同直径的孔隙和孔洞在所述碳酸岩样品上的分布方式。
4.如权利要求3所述的方法,其特征在于,根据轮廓图的分布方式确定不同直径的孔隙和孔洞在所述碳酸岩样品上的分布方式,包括:
根据孔洞和孔隙的轮廓图,获取孔洞和孔隙的骨架图,骨架图是由线段组成的表示孔洞或孔隙的拓扑结构的图形;
根据骨架图在碳酸岩样品的CT图像上的坐标,确定孔洞或孔隙在所述碳酸岩样品上的位置。
5.如权利要求3所述的方法,其特征在于,还包括:
对于每个轮廓图,将轮廓图划分为多个图块,在每个图块中建立与图块边缘相切的椭圆,多个图块中椭圆的短轴方向相同;
将每个图块短轴的长度与该图块面积相乘,得到每个图块的乘积,将多个图块的乘积相加,得到的和为与该轮廓图边缘相切的椭圆的短轴长度,确定为孔洞或孔隙的短轴长度。
6.一种碳酸岩多尺度孔洞特征的确定装置,其特征在于,包括:
第一图像扫描模块,用于对直径大于全直径的碳酸岩样品进行毫米CT扫描,获得第一扫描图像;
第二图像扫描模块,用于在所述碳酸岩样品上取直径为微米级的孔隙或孔洞分布密集的部分作为第一柱塞样品,对该第一柱塞样品进行微米CT扫描,获得第二扫描图像;
第三图像扫描模块,用于在所述第一柱塞样品上取直径为纳米级的孔隙或孔洞分布密集的部分作为第二柱塞样品,对该第二柱塞样品进行纳米CT扫描,获得第三扫描图像;
确定模块,用于对所述第一扫描图像、所述第二扫描图像和所述第三扫描图像进行二值化处理,根据二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像,确定所述碳酸岩样品的孔隙度以及不同直径的孔隙和孔洞在所述碳酸岩样品上的分布方式,其中,在二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像中,孔隙和孔洞的像素数相同,孔隙和孔洞的像素数与所述碳酸岩样品的颗粒的像素数不 同
在所述第一扫描图像中,孔隙或孔洞的直径大于等于0.5毫米,在所述第二扫描图像中,孔隙或孔洞的直径大于等于0.2微米小于0.5毫米,在所述第三扫描图像中,孔隙或孔洞的直径大于等于50纳米小于0.2微米。
7.如权利要求6所述的装置,其特征在于,所述确定模块,包括:
第一面积计算单元,用于根据孔隙的像素数分别计算二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像中孔隙的面积,将二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像中孔隙的面积相加得到所有孔隙的总面积;
第二面积计算单元,用于根据二值化处理后的第一扫描图像中孔隙的像素数和颗粒的像素数计算所述碳酸岩样品的CT图像面积;
孔隙度确定单元,用于将所述所有孔隙的总面积与碳酸岩样品的CT图像面积的比值,确定为所述碳酸岩样品的孔隙度。
8.如权利要求6所述的装置,其特征在于,所述确定模块,还包括:
轮廓图处理单元,用于将二值化处理后的第一扫描图像、第二扫描图像和第三扫描图像中孔洞和孔隙的图形处理为轮廓图,其中,轮廓图是孔洞或孔隙的图形的轮廓线形成的封闭的图形;
分布确定单元,用于根据轮廓图的分布方式确定不同直径的孔隙和孔洞在所述碳酸岩样品上的分布方式。
9.如权利要求8所述的装置,其特征在于,所述确定模块,还包括:
骨架图处理单元,用于根据孔洞和孔隙的轮廓图,获取孔洞和孔隙的骨架图,骨架图是由线段组成的表示孔洞或孔隙的拓扑结构的图形;
位置确定单元,用于根据骨架图在碳酸岩样品的CT图像上的坐标,确定孔洞或孔隙在所述碳酸岩样品上的位置。
10.如权利要求9所述的装置,其特征在于,还包括:
图块划分模块,用于对于每个轮廓图,将轮廓图划分为多个图块,在每个图块中建立与图块边缘相切的椭圆,多个图块中椭圆的短轴方向相同;
轴长确定模块,用于将每个图块短轴的长度与该图块面积相乘,得到每个图块的乘积,将多个图块的乘积相加,得到的和为与该轮廓图边缘相切的椭圆的短轴长度,确定为孔洞或孔隙的短轴长度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510169106.9A CN104751473B (zh) | 2015-04-10 | 2015-04-10 | 碳酸岩多尺度孔洞特征的确定方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510169106.9A CN104751473B (zh) | 2015-04-10 | 2015-04-10 | 碳酸岩多尺度孔洞特征的确定方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104751473A CN104751473A (zh) | 2015-07-01 |
CN104751473B true CN104751473B (zh) | 2017-10-27 |
Family
ID=53591092
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510169106.9A Active CN104751473B (zh) | 2015-04-10 | 2015-04-10 | 碳酸岩多尺度孔洞特征的确定方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104751473B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105649615B (zh) * | 2015-12-28 | 2019-01-18 | 中国石油天然气股份有限公司 | Ct定量、三维可视化测试储层致密油赋存状态的方法 |
CN106918607B (zh) * | 2017-03-06 | 2019-04-05 | 西南石油大学 | 一种孔隙结构获取方法及装置 |
CN109187305A (zh) * | 2018-08-24 | 2019-01-11 | 中国地质大学(武汉) | 一种页岩结构表征方法和装置 |
CN111257188B (zh) * | 2018-11-30 | 2022-10-04 | 中国石油天然气股份有限公司 | 碳酸盐岩孔隙***的多尺度表征方法及装置 |
CN109736794A (zh) * | 2019-01-24 | 2019-05-10 | 中国石油天然气股份有限公司 | 裂缝性孔隙型致密砂岩储层特征确定方法、装置及*** |
CN110441220A (zh) * | 2019-09-02 | 2019-11-12 | 西南石油大学 | 微生物碳酸盐岩的孔隙表征方法、装置及服务器 |
CN112986087B (zh) * | 2019-12-02 | 2023-04-07 | 中国石油天然气股份有限公司 | 表征多尺度碳酸盐岩中孔、洞和缝的方法和装置 |
CN111366521B (zh) * | 2020-04-14 | 2023-08-11 | 清能艾科(深圳)能源技术有限公司 | 多尺度确定孔隙度的方法及相关设备 |
CN112084884A (zh) * | 2020-08-19 | 2020-12-15 | 厦门理工学院 | 一种扫描电镜图像孔隙识别方法、终端设备及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101556703A (zh) * | 2009-05-16 | 2009-10-14 | 中国石油大学(华东) | 基于连续切片图像的网络模型建立方法 |
CN102222359A (zh) * | 2011-05-24 | 2011-10-19 | 中国石油天然气股份有限公司 | 一种岩心三维孔隙结构重构方法 |
CN103575751A (zh) * | 2013-09-10 | 2014-02-12 | 济南大学 | 一种基于ct技术的多孔混凝土孔结构表征方法 |
CN103698803A (zh) * | 2012-09-27 | 2014-04-02 | 中国石油天然气股份有限公司 | 一种岩石孔隙结构表征方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8583410B2 (en) * | 2010-05-28 | 2013-11-12 | Ingrain, Inc. | Method for obtaining consistent and integrated physical properties of porous media |
-
2015
- 2015-04-10 CN CN201510169106.9A patent/CN104751473B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101556703A (zh) * | 2009-05-16 | 2009-10-14 | 中国石油大学(华东) | 基于连续切片图像的网络模型建立方法 |
CN102222359A (zh) * | 2011-05-24 | 2011-10-19 | 中国石油天然气股份有限公司 | 一种岩心三维孔隙结构重构方法 |
CN103698803A (zh) * | 2012-09-27 | 2014-04-02 | 中国石油天然气股份有限公司 | 一种岩石孔隙结构表征方法及装置 |
CN103575751A (zh) * | 2013-09-10 | 2014-02-12 | 济南大学 | 一种基于ct技术的多孔混凝土孔结构表征方法 |
Non-Patent Citations (2)
Title |
---|
Pore Scale Characterization of Carbonates at Multiple Scales:Integration of MicroCT,BSEM and FIBSEM;Robert M.Sok 等;《Petrophysics》;20101231;第51卷(第6期);379-387 * |
利用多尺度CT成像表征致密砂岩微观孔喉结构;白斌 等;《石油勘探与开发》;20130630;第40卷(第3期);329-333 * |
Also Published As
Publication number | Publication date |
---|---|
CN104751473A (zh) | 2015-07-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104751473B (zh) | 碳酸岩多尺度孔洞特征的确定方法及装置 | |
CN104535475B (zh) | 碳酸盐岩微观结构的确定方法及装置 | |
CN107817199B (zh) | 一种致密砂岩多尺度孔隙模型的构建方法和应用 | |
CN104237103B (zh) | 一种孔隙连通性定量表征方法及装置 | |
CN105487121B (zh) | 基于ct扫描图像与电成像图像融合构建多尺度数字岩心方法 | |
CN105115874B (zh) | 多源信息融合的多组分三维数字岩心构建方法 | |
CN104619952B (zh) | 具可靠多相渗透性决定数字岩石分析***及方法 | |
CN113609696B (zh) | 基于图像融合的多尺度多组分数字岩心构建方法及*** | |
CN104729972B (zh) | 确定碳酸盐岩气水相对渗透率的方法及装置 | |
CN108802073A (zh) | 一种基于数字岩心的岩石电性参数获取方法及装置 | |
Moreno et al. | Tillage and soil type effects on soil surface roughness at semiarid climatic conditions | |
CN107368670A (zh) | 基于大数据深度学习的胃癌病理诊断支持***和方法 | |
CN106355166A (zh) | 一种基于监控视频与遥感图像的灰霾扩散路径绘制与来源确定方法 | |
CN108876901A (zh) | 一种基于二维图像和多点统计学的数字岩心重建方法 | |
CN109490165A (zh) | 表征碳酸盐岩非组构选择性储集空间的方法 | |
CN108956416A (zh) | 一种基于Matlab分析致密砂岩储层孔隙表征的方法 | |
CN110320137A (zh) | 一种基于数字岩心的多尺度融合方法 | |
CN107132171B (zh) | 一种基于压汞-氮吸附联测数据确定致密储层孔径分布的方法 | |
CN103698803A (zh) | 一种岩石孔隙结构表征方法及装置 | |
JP2011214983A (ja) | シミュレーション装置及びそのコンピュータプログラム | |
CN103971395B (zh) | 一种纤维过滤介质结构的拟态化重建及其性能计算方法 | |
CN108061697B (zh) | 土体三维孔隙率计算方法 | |
CN110441209A (zh) | 一种基于致密储层数字岩心计算岩石渗透率的方法 | |
CN109387468A (zh) | 页岩储层纳米孔隙结构特征参数测试分析方法及*** | |
CN109635629A (zh) | 一种基于深度学习的公交站台人群密度检测方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C41 | Transfer of patent application or patent right or utility model | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20151013 Address after: 100007 Beijing, Dongzhimen, North Street, No. 9, No. Applicant after: PetroChina Company Limited Address before: 102249 Beijing city Changping District Road No. 18 Applicant before: China University of Petroleum (Beijing) |
|
GR01 | Patent grant | ||
GR01 | Patent grant |