[发明专利]一种页岩油储层有机孔隙和无机孔隙的划分与表征方法

权利要求书

1.一种页岩油储层有机孔隙和无机孔隙的划分与表征方法，其特征在于，包括以下步骤：

步骤一、岩心图像序列采集：对目标页岩油储层段进行取心，钻取直径为1英寸的岩样，表面用砂纸及氩离子束抛光，利用聚焦离子束-扫描电镜仪器对岩样进行连续切割扫描，采集不少于700张SEM图像；

步骤二、图像对齐：通过图像之间的全局平移和旋转进行图像对齐处理，二维图像中，坐标点(x,y)经刚体变换到点(x',y')的变换公式为：

其中为旋转角度，Δx、Δy为平移距离；

步骤三、几何校正：通过对SEM图像乘以校正因子将SEM图像尺寸变换为样品截面的实际尺寸；

步骤四、图像分割：根据孔隙灰度分布特征，采用双阈值对岩心数据进行两次分割，通过逻辑运算与代数运算实现有机孔隙和无机孔隙的精细划分；

步骤五、孔隙结构定量表征：分析得有机孔隙和无机孔隙孔径分布，从而实现页岩油储层有机孔隙和无机孔隙的定量表征；

所述步骤四的具体步骤为：

S1、基于有机孔隙和无机孔隙的分布特征，确定孔隙灰度分布特征，由于景深的影响，有机孔隙灰度值最小，有机质次之，无机质灰度值最大，无机孔隙较为复杂，一部分无机孔隙灰度值接近有机孔隙或者有机质，一部分无机孔隙灰度值介于有机质和无机质之间；对三维岩心数据二值分割后，图像变成黑白两部分，即0和1，对同一块岩心数据进行两次二值分割后，对数据体进行图像逻辑运算与代数运算可以得到有机孔隙和无机孔隙空间分布；

S2、第一次二值化分割，以有机孔隙最大灰度值为分割阈值，将岩心分割成两部分，分别为：有机孔隙、部分无机孔隙，记为A，有机质、无机质、部分无机孔隙，记为B；

S3、第二次二值化分割，以无机孔隙最大灰度值为分割阈值，将岩心分割成两部分，分别为：有机孔隙、有机质、无机孔隙，记为C，无机质记为D；

S4、获取无机孔隙数据体，有机孔隙分布在有机质中，无机孔隙分布在无机质中，采用MATLAB中imfill函数对数据体D中孔隙相进行填充可得到所有无机孔隙和无机质数据体E，利用数据体E与数据体D做减法运算，得到无机孔隙数据体INOP；

S5、对数据体E进行逻辑取反操作，得到所有有机孔隙和有机质数据体F，有机孔隙分布在有机质中，对数据体A和数据体F进行逻辑与运算，可得到有机孔隙数据体OP。