| 一种岩石内部流体饱和度空间分布的定量分析装置及方法 |
| 江文滨; 林缅 ; 姬莉莉; 曹高辉; 徐志朋; 郑思平
|
| 2021-10-19
|
| Rights Holder | 中国科学院力学研究所
|
| Abstract | 本发明属于石油科研技术领域。针对当前岩石内流体饱和分布的定量分析并未达到孔隙尺度的现状，本发明提供了一种岩石内部流体饱和度空间分布的定量分析装置及方法，覆压状态下无液体、注入流体岩心三维成像D1、D2，D1和D2三维配准；由D2采样得到与D1匹配的D3；基于覆压孔隙度确定阈值T1，D1转换为空隙、基质二值矩阵；基于总流体饱和度确定阈值T2，根据D3中空隙像素的灰度，修改D1中相应像素取值；基于D1的空隙像素提取孔隙网络，获得孔隙、孔喉与对应像素列表；统计孔隙、孔喉中两种流体所占像素数及总像素数，得到孔隙流体饱和度；统计不同流体饱和度的孔隙数量、空间分布及与形状因子关系。本发明能够无损、定量获得岩心内部孔隙尺度的流体饱和度的空间分布。 |
| Application Date | 2020-12-08
|
| Application Number | CN202011441351.8
|
| Patent Number | CN112577979B
|
| Claim | 1.一种岩石内部流体饱和度空间分布的定量分析方法，采用一种岩石内部流体饱和度空间分布的定量分析装置进行数据采集，再进行数据分析，其特征在于，具体包括如下步骤：
（1）首先在夹持器中放入经过烘干后、孔隙中不含液体的岩石样品，施加一定的环压和轴压，待稳定后，进行CT扫描，获得第1个三维数据体，D1，像素灰度与对X射线吸收能力正相关；
（2）以一定的流速或一定的压力向岩石样品中注入某种流体或两种流体，待稳定后，进行CT扫描，获得第2个三维数据体，D2，像素灰度与对X射线吸收能力正相关；
（3）对于D1和D2数据体进行三维配准，采用图像配准算法进行D1和D2数据体的三维配准，通过重采样算法，获得D1数据体中各像素点对应的D2值，该数据体命名为D3，即，D1和D3数据体中每个像素对应的是岩石中三维空间中的同一个点；
（4）对D1数据体，以相同覆压状态下氦气测孔隙度为约束，确定孔隙和岩石骨架像素的划分阈值T1，使得划分后的孔隙像素数占岩石总像素数的比例与氦气测孔隙度一致，将D1数据转换为三维二值矩阵，其中孔隙为1，骨架为0；对于D1数据中骨架像素，将D3中相应像素的值设为0；
（5）对于D3数据体，以入口端注入两相流体量和出口端计量的两相流体量之差确定的样品内的流体饱和度为约束，确定两相流体的划分阈值T2，对于D1数据体中取值为1的所有像素，如果D3中对应像素的灰度大于阈值T2，则D1中像素值改为2；
（6）对D1数据，以取值为1或2的像素为孔隙像素，采用中轴-最大球孔隙网络抽提算法，进行孔隙网络抽提，得到若干的孔隙和孔喉，获得孔隙和孔喉的参数和孔隙与孔喉的拓扑连接关系，建立孔隙和孔喉与D1数据体中对应像素的映射关系；
（7）遍历所有的孔隙和孔喉，基于映射关系确定D1数据体中的对应像素，统计其中取值为1的像素个数及总的像素个数，两者之比即为该孔隙/孔喉中对X射线吸收较少的流体的饱和度S1，像素个数乘以像素尺寸的三次方即为该流体占据的体积；1-S1即为该孔隙/孔喉中另一流体的饱和度S2，像素个数乘以像素尺寸的三次方即为另一流体占据的体积；
（8）统计S2分别为从0%到100%、以10%为间隔的孔隙的数量分布，统计沿流动方向和垂直流动方向上S2大于一定值的孔隙数量与总的孔隙数量的比值变化，统计不同形状因子孔隙的S2平均值；
所述一种岩石内部流体饱和度空间分布的定量分析装置，其特征在于，包括X射线三维显微成像系统，
所述X射线三维显微成像系统包括射线源、夹持器、样品台和接收器，所述射线源和接收器相互对应安装；岩石样品设置为圆柱体形状，其安装于X射线穿透的夹持器中，夹持器固定安装在样品台上，通过样品台转动控制夹持器转动；成像时，夹持器随着样品台转动，接收器记录不同角度的图像，通过三维重建算法建立三维数据体。 |
| Language | 中文
|
| Classification | 发明授权
|
| Status | 有效
|
| Note | 授权
|
| Country | 中国
|
| Agency | 北京和信华成知识产权代理事务所(普通合伙)
|
| Document Type | 专利
|
| Identifier | http://dspace.imech.ac.cn/handle/311007/88107
|
| Collection | 流固耦合系统力学重点实验室
|
Recommended Citation
GB/T 7714 |
江文滨,林缅,姬莉莉,等. 一种岩石内部流体饱和度空间分布的定量分析装置及方法. CN112577979B[P]. 2021-10-19.
|
|
File name:
|
20211019_0C_CN_0 (2).pdf
|
|
Format:
|
Adobe PDF
|
This file does not support browsing at this time
Items in the repository are protected by copyright, with all rights reserved, unless otherwise indicated.
Edit Comment