[发明专利]微米尺度平面孔隙网络内气液驱替图像的饱和度获取方法

权利要求书

1.微米尺度平面孔隙网络内气液驱替图像的饱和度获取方法，其特征在于，包括以下步骤：

步骤一，获取平面孔隙网络模型图像，对该图像进行高斯滤波除去噪声后转化为灰度图像，对灰度图像进行坐标划分，并利用阈值分割方法将灰度图像转化像素尺寸为宽L_x和高L_y的黑白二值图像；

步骤二，裁剪黑白二值图像中入口与出口区域，保留中间孔隙网络区域，生成校验图像；

步骤三，统计校验图像中白色区域像素数，获得白色区域总面积，计算流动区域总面积A_total；

步骤四，获取平面孔隙网络气相驱替形态图像，裁剪平面孔隙网络气相驱替形态图像入口与出口区域，并进行灰度阈值分割生成黑白二值待测图像；

步骤五，统计黑白二值待测图像中白色区域像素数，获得白色区域总面积，计算气相和液相流动区域面积；

步骤六，对比校验图像与黑白二值待测图像，获得气液两相流动区域驱替形态图像，并求取气相饱和度s_g和液相饱和度s_l以及驱替流量q_g和q_l。

2.根据权利要求1所述的微米尺度平面孔隙网络内气液驱替图像的饱和度获取方法，其特征在于，步骤一中对灰度图像进行坐标划分：以平面孔隙网络模型灰度图像的左上角为坐标原点，建立直角坐标系，将图像分割为若干区块。

3.根据权利要求2所述的微米尺度平面孔隙网络内气液驱替图像的饱和度获取方法，其特征在于，求取每个区块的灰度阈值，将其转化为黑白二值图像，黑白二值图像中任意坐标(x,y)处的灰度值为I(x,y)，黑色区域的灰度I(x,y)＝255，白色区域的灰度I(x,y)＝0，坐标(x,y)处的像素值记为p(x,y)；

4.根据权利要求1所述的微米尺度平面孔隙网络内气液驱替图像的饱和度获取方法，其特征在于，裁剪黑白二值图像入口与出口区域的具体过程如下：

(1)将黑白二值图像的左边界上所有坐标点的像素值p(0,y_n)与其相邻的坐标点的像素值p(0,y_n-1)和p(0,y_n+1)分别进行“与”运算，搜索运算结果为0的坐标点A_i(0,y_A)和B_i(0,y_B)即为二值图像左边界上像素灰度变化点；

(2)将平面孔隙网络右边界上所有坐标点的像素值p(L_x,y_m)与其相邻的坐标点p(L_x,y_m-1)和p(L_x,y_m+1)的像素点进行“与”运算，搜索运算结果为0的坐标点C_o(L_x,y_C)、D_o(L_x,y_D)即为二值图像右边界上像素灰度变化点；

(3)分别沿直线y＝y_A和直线y＝y_B向x轴正方向按步骤(1)的方法搜索获取左边界上像素灰度变化点A(x_A,y_A)、B(x_B,y_B)；

(4)分别沿直线y＝y_C和直线y＝y_D向x轴负方向按步骤(2)的方法搜索获取右边界上像素灰度变化点C(x_C,y_C)、D(x_D,y_D)；

(5)以直线A-B和直线C-D为界限，裁剪黑白二值图像，去除平面孔隙网络入口和出口区域，只保留中间孔隙网络区域，获得像素尺寸为宽(x_D-x_A)和高L_y或宽(x_C-x_B)和高L_y的校验图像。