一种岩心ct扫描分析方法及装置制造方法
【专利摘要】本申请公开了一种岩心CT扫描分析方法及装置，该方法包括如下步骤：获取岩心样品的CT图像，该CT图像含有岩心样品中所有像素点的CT值；提取目标区域内CT图像上每个像素点的CT值，该目标区域包含有n×m个像素点，其中n,m为正整数；利用所提取的CT值，计算目标区域内CT图像中每个像素点的储层评价参数；构建目标区域内所有像素点的储层评价参数的频率分布曲线；利用所述储层评价参数的频率分布曲线，对所述岩心样品进行分析。该装置相应的包括CT图像获取单元、CT值提取单元、储层评价参数计算单元、频率分布曲线构建单元以及分析单元。通过本申请中的技术方案能实现从微观视角来定量分析岩心及其中流体分布特征的目的。
【专利说明】-种岩心CT扫描分析方法及装置
【技术领域】
[0001] 本申请涉及岩心储层评价分析【技术领域】，尤其涉及一种岩心CT扫描分析方法及 装直。
【背景技术】
[0002] 计算机断层扫描（Computed Tomography，简称CT)技术是一种专门用来进行无损 检测和探伤的技术。近年来，随着计算机工业的发展，CT扫描系统性能不断提高，在石油工 业许多领域得到了广泛的应用。目前，CT扫描技术已成功应用于油气田开发中的储层评价 和渗流规律研究，如岩心非均质性表征、岩心在线饱和度测量以及岩心驱替流动实验评价 等领域。
[0003] CT扫描技术最大的优势就是其能够在无损岩心的前提下获取岩心内部的信息，而 CT扫描成像原理也决定其满足能够在较小尺度范围内开展定量分析的先决条件。因此，CT 扫描技术除了无损岩心特性外，其具备从微观视角来定量分析岩心内部特征的特性。
[0004] 不过就目前而言，现有的CT扫描分析方法普遍还是基于宏观角度，如获取岩心整 体的平均孔隙度和平均饱和度以及孔隙度和饱和度参数的沿程分布特性等，缺少微观定量 的分析手段，微观定量的分析结果也很少，而未来油层物理的发展趋势必定是朝着更微观 更定量的方向发展。因此，从微观角度来定量分析岩心内部特征显得很有必要。
【发明内容】
[0005] 本申请实施例的目的是提供一种岩心CT图像频率分布分析方法及装置，以实现 从微观视角来定量分析岩心及其中流体分布特征。
[0006] 为解决上述技术问题，本申请实施例通过以下技术方案来实现：
[0007] 本申请提供了一种岩心CT扫描分析方法，包括如下步骤：
[0008] 获取岩心样品的计算机断层扫描CT图像，所述CT图像含有岩心样品中所有像素 点的CT值；
[0009] 提取目标区域内CT图像上每个像素点的CT值，所述目标区域包含有nXm个像素 点，其中n，m为正整数；
[0010] 利用所提取的CT值，计算目标区域内CT图像上每个像素点的储层评价参数；
[0011] 构建目标区域内所有像素点的储层评价参数的频率分布曲线；
[0012] 利用所述储层评价参数的频率分布曲线，对所述岩心样品进行分析。
[0013] 优选的，所述计算目标区域内CT图像中每个像素点的储层评价参数包括：
[0014] 将所述nXm个像素点整合成N个像素点，其中# = 4^，N和M均为正整数并且 M M>1 ；
[0015] 计算整合成的N个像素点的CT值；
[0016] 利用计算得到的CT值，计算整合成的N个像素点的储层评价参数。
[0017] 优选的，所述计算目标区域内CT图像上每个像素点的储层评价参数包括：
[0018] 计算目标区域内CT图像上每个像素点的孔隙度和流体饱和度。
[0019] 优选的，所述计算整合成的N个像素点的储层评价参数包括：
[0020] 计算整合成的N个像素点的孔隙度和流体饱和度。
[0021] 优选的，所述计算目标区域内CT图像中每个像素点的孔隙度包括利用下述公式 计算有伤害储层的岩心样品在目标区域内每个像素点的孔隙度小：
[0022]
【权利要求】
1. 一种岩心CT扫描分析方法，其特征在于，包括如下步骤： 获取岩心样品的计算机断层扫描CT图像，所述CT图像含有岩心样品中所有像素点的 CT值； 提取目标区域内CT图像上每个像素点的CT值，所述目标区域包含有nXm个像素点， 其中n，m为正整数； 利用所提取的CT值，计算目标区域内CT图像上每个像素点的储层评价参数； 构建目标区域内所有像素点的储层评价参数的频率分布曲线； 利用所构建的储层评价参数的频率分布曲线，对所述岩心样品进行分析。
2. 根据权利要求1所述的岩心CT扫描分析方法，其特征在于，所述计算目标区域内CT 图像中每个像素点的储层评价参数包括： η X in 将所述nXm个像素点整合成N个像素点，其中W =i^-，N和M均为正整数并且M>1 ; 计算整合成的N个像素点的CT值； 利用计算得到的CT值，计算整合成的N个像素点的储层评价参数。
3. 根据权利要求1所述的岩心CT扫描分析方法，其特征在于，所述计算目标区域内CT 图像中每个像素点的储层评价参数包括： 计算目标区域内CT图像中每个像素点的孔隙度和流体饱和度。
4. 根据权利要求2所述的岩心CT扫描分析方法，其特征在于，所述计算整合成的N个 像素点的储层评价参数包括： 计算整合成的N个像素点的孔隙度和流体饱和度。
5. 根据权利要求3或4所述的岩心CT扫描分析方法，其特征在于，所述计算目标区域 内CT图像中每个像素点的孔隙度包括利用下述公式计算有伤害储层的岩心样品在目标区 域内各个像素点的孔隙度Φ :
上式中CTt为驱替过程中t时刻岩心断层面的CT值；CTgrain为岩石骨架的CT值，S fi为 驱替过程中岩心断层面中所含第i种流体的饱和度，CTi为第i种流体的CT值，j为驱替过 程中岩心断层面中所含流体的种类数。
6. 根据权利要求5所述的岩心CT扫描分析方法，其特征在于，所述计算目标区域内CT 图像中每个像素点的流体饱和度包括利用下述公式计算岩心样品在目标区域内各个像素 点的流体饱和度：
上式中CTtoy为干岩心断层面的CT值；CTfluidwrt为岩心100%饱和流体后断层面的CT 值；CTt为驱替过程中时刻t岩心断层面的CT值；CTfuild为驱替流体的CT值；CT fi为岩心样 品中所含第i种流体的CT值；CT&为空气的CT值，Sfi为驱替过程中岩心断层面中所含第 i种流体的饱和度。
7. -种岩心CT扫描分析装置，其特征在于，包括： CT图像获取单元，所述CT图像获取单元用于获取岩心样品的计算机断层扫描CT图像， 所述CT图像含有岩心样品中所有像素点的CT值； CT值提取单元，所述CT值提取单元用于提取目标区域内CT图像上每个像素点的CT 值，所述目标区域包含有nXm个像素点,其中n, m为正整数； 储层评价参数计算单元，所述储层评价参数计算单元用于利用所提取的CT值，计算目 标区域内CT图像上每个像素点的储层评价参数； 频率分布曲线构建单元，所述频率分布曲线构建单元用于构建目标区域内所有像素点 的储层评价参数的频率分布曲线； 分析单元，所述分析单元用于利用所述储层评价参数的频率分布曲线，对所述岩心样 品进行储层评价分析。
8. 根据权利要求7所述的岩心CT扫描分析装置，其特征在于，储层评价参数计算单元 包括： 第一子单元，所述第一子单元用于将所述nXm个像素点整合成N个像素点，其中
N和M均为正整数； 第二子单元，所述第二子单元用于计算整合成的N个像素点的CT值； 第三子单元，所述第三子单元用于利用计算得到的CT值，计算整合成的N个像素点的 储层评价参数。
9. 根据权利要求7所述的岩心CT扫描分析装置，其特征在于，所述储层参数计算单元 还包括计算子单元，所述计算子单元用于计算目标区域内CT图像中每个像素点的孔隙度 和流体饱和度。
10. 根据权利要求8所述的岩心CT扫描分析装置，其特征在于，第三子单元包括第三计 算子单元，所述第三计算子单元用于计算整合成的N个像素点的孔隙度和流体饱和度。
【文档编号】G01N15/08GK104374682SQ201410645526
【公开日】2015年2月25日 申请日期:2014年11月12日 优先权日:2014年11月12日
【发明者】冷振鹏, 刘庆杰, 吕伟峰, 马德胜, 李彤, 贾宁洪, 杨思玉, 吴康云 申请人:中国石油天然气股份有限公司