微米尺度平面孔隙网络内气液驱替图像的饱和度获取方法技术

本发明专利技术公开了微米尺度平面孔隙网络内气液驱替图像的饱和度获取方法，获取平面孔隙网络模型图像，对该图像进行高斯滤波除去噪声、坐标划分得到黑白二值图像，然后对黑白二值图像中入口与出口区域裁剪统计校验图像中白色区域像素数，获得白色区域总面积，同时平面孔隙网络气相驱替形态图像入口与出口区域，并进行灰度阈值分割生成黑白二值待测图像；统计黑白二值待测图像中白色区域像素数获得白色区域总面积，对比校验图像与黑白二值待测图像，获得气液两相流动区域驱替形态图像从而计算得到气相饱和度和液相饱和度以及驱替流量，本发明专利技术测量结果稳定，鲁棒性强，本方法可高效、精确地处理大量图像数据，特别是对低驱替速度下的饱和度分析和驱替流量计算，更加便捷、有效。

Saturation acquisition method of gas-liquid displacement image in micron scale planar pore network

【技术实现步骤摘要】
微米尺度平面孔隙网络内气液驱替图像的饱和度获取方法
本专利技术属于油气藏强化开采领域，适用于二氧化碳地质封存领域的微尺度平面孔隙多孔介质中气液两相流动实验，具体涉及一种微米尺度平面孔隙网络内气液两相驱替图像的饱和度获取方法。
技术介绍
平面孔隙网络中的气液两相驱替研究是获取天然多孔介质中气液两相驱替机理的主要手段。天然多孔介质中的气液两相驱替过程发生在众多的工业过程中，例如致密油气藏的气驱油、水驱气、二氧化碳的地质封存等过程。毛管力与流体饱和度的关系曲线，即毛管力曲线，是油气储层中不混溶两相渗流的基本特征。流体相对渗透率与其饱和度的关系曲线，即相对渗透率曲线，是研究不混溶流体渗流的基础。不混溶的流体在多孔介质中驱替渗流时，每种流体会占据一定的孔隙体积，多孔介质孔隙中任意一点P处m相流体的饱和度定义为P点周围的特征体V*内的m相流体所占据的孔隙体积的百分比，即其中是特征体V*内的孔隙体积，是特征体V*内的m相流体的体积，特征体内各相流体饱和度之和为1。若取特征体为整个多孔介质，则多孔介质内m流体的饱和度sm为多孔介质孔隙Vp内m相流体所占据的孔隙体积的百分比，即其中Vm是整个多孔介质孔隙内m相流体的体积，多孔介质内各相流体饱和度之和为1。毛管力曲线和相对渗透率曲线广泛应用于油气田开发参数计算、动态分析以及油藏数值模拟等，是油气田开发中的两个重要关系曲线。获得毛管力曲线和相对渗透率曲线的关键参数是确定各相流体的饱和度。在微尺度平面孔隙网络中进行两相不混溶流体驱替实验是研究多孔介质渗流机理的常用方法。微米尺度平面孔隙网络内的渗流是小流量驱替工况，特别是进行准静态驱替实验时，此时难以通过实时测量各相流体的流量来确定孔隙内各相流体的体积，进一步地难以确定各相流体饱和度。油藏岩心样品中进行两相不混溶流体驱替实验时，往往先要测量多孔介质的孔隙度，再通过测量流经多孔介质的入口和出口流量之差来确定多孔介质内的流体的饱和度，该方法对多孔介质入口与出口处的流量传感器精度要求较高。以上方法不适用微米尺度的平面孔隙网络内的流体饱和度的测量。
技术实现思路
本专利技术的目的在于提供微米尺度平面孔隙网络内气液驱替图像的饱和度获取方法，以克服现有技术的不足，本专利技术能够从平面孔隙网络气液两相驱替图像中识别出孔隙网络的固体区域、各流体流动区域、相界面，同时测量出各相流体流动区域的面积，各相流体的饱和度以及驱替流量。为达到上述目的，本专利技术采用如下技术方案：微米尺度平面孔隙网络内气液驱替图像的饱和度获取方法，包括以下步骤：步骤一，获取平面孔隙网络模型图像，对该图像进行高斯滤波除去噪声后转化为灰度图像，对灰度图像进行坐标划分，并利用阈值分割方法将灰度图像转化像素尺寸为宽Lx和高Ly的黑白二值图像；步骤二，裁剪黑白二值图像中入口与出口区域，保留中间孔隙网络区域，生成校验图像；步骤三，统计校验图像中白色区域像素数，获得白色区域总面积，计算流动区域总面积Atotal；步骤四，获取平面孔隙网络气相驱替形态图像，裁剪平面孔隙网络气相驱替形态图像入口与出口区域，并进行灰度阈值分割生成黑白二值待测图像；步骤五，统计黑白二值待测图像中白色区域像素数，获得白色区域总面积，计算气相和液相流动区域面积；步骤六，对比校验图像与黑白二值待测图像，获得气液两相流动区域驱替形态图像，并求取气相饱和度sg和液相饱和度sl以及驱替流量qg和ql。进一步的，步骤一中对灰度图像进行坐标划分：以平面孔隙网络模型灰度图像的左上角为坐标原点，建立直角坐标系，将图像分割为若干区块。进一步的，求取每个区块的灰度阈值，将其转化为黑白二值图像，黑白二值图像中任意坐标(x,y)处的灰度值为I(x,y)，黑色区域的灰度I(x,y)＝255，白色区域的灰度I(x,y)＝0，坐标(x,y)处的像素值记为p(x,y)；进一步的，裁剪黑白二值图像入口与出口区域的具体过程如下：(1)将黑白二值图像的左边界上所有坐标点的像素值p(0,yn)与其相邻的坐标点的像素值p(0,yn-1)和p(0,yn+1)分别进行“与”运算，搜索运算结果为0的坐标点Ai(0,yA)和Bi(0,yB)即为二值图像左边界上像素灰度变化点；(2)将平面孔隙网络右边界上所有坐标点的像素值p(Lx,ym)与其相邻的坐标点p(Lx,ym-1)和p(Lx,ym+1)的像素点进行“与”运算，搜索运算结果为0的坐标点Co(Lx,yC)、Do(Lx,yD)即为二值图像右边界上像素灰度变化点；(3)分别沿直线y＝yA和直线y＝yB向x轴正方向按步骤(1)的方法搜索获取左边界上像素灰度变化点A(xA,yA)、B(xB,yB)；(4)分别沿直线y＝yC和直线y＝yD向x轴负方向按步骤(2)的方法搜索获取右边界上像素灰度变化点C(xC,yC)、D(xD,yD)；(5)以直线A-B和直线C-D为界限，裁剪黑白二值图像，去除平面孔隙网络入口和出口区域，只保留中间孔隙网络区域，获得像素尺寸为宽(xD-xA)和高Ly或宽(xC-xB)和高Ly的校验图像。进一步的，步骤三中求取流动区域总面积的具体过程如下：(1)从步骤二中生成的校验图像的左上角开始计数，搜索整个校验图像区域，白色区域的像素数为整个校验图像区域所有的像素值之和；(2)按照校验图像像素尺寸与长度尺寸的比例，将像素数换算为白色区域总面积，即为平面孔隙网络中流动区域总面积Atotal。进一步的，步骤四中裁剪气相驱替形态图像：(1)以平面孔隙网络中气相驱替形态图像的左上角为坐标原点，建立直角坐标系；以步骤二中获得的左边界上像素灰度变化点A(xA,yA)、B(xB,yB)建立直线A-B，以右边界上像素灰度变化点C(xC,yC)、D(xD,yD)建立直线C-D；(2)以直线A-B和直线C-D为界限，裁剪二值图像，去除气相驱替形态图像的入口和出口区域，只保留中间区域，获得像素尺寸为宽(xD-xA)和高Ly或宽(xC-xB)和高Ly的驱替形态图像。进一步的，将裁剪完成的驱替形态图像分割为若干区块，取每个区块的灰度阈值，将其转化为黑白二值图像，生成驱替形态待测图像，黑白二值图像中任意坐标(x,y)处的灰度值为I(x,y)，黑色区域的灰度I(x,y)＝255，白色区域的灰度I(x,y)＝0，坐标(x,y)处的像素值记为p(x,y)；进一步的，步骤五中获取气相和液相流动区域面积的具体过程如下：(1)从步骤四生成的黑白二值待测图像的左上角开始计数，搜索整个黑白二值待测图像区域，白色区域的像素数为整个黑白二值待测图像区域所有的像素值之和；(2)按照黑白二值待测图像像素尺寸与长度尺寸的比例，将白色区域像素数换算为白色区域总面积，即为气液两相驱替中气相流动面积Ag，流动区域总面积Atotal减去气相流动面积即得液相流动面积Al：Al＝Atotal-Ag(5)。进一步的，步骤六中获得气液本文档来自技高网...

【技术保护点】
1.微米尺度平面孔隙网络内气液驱替图像的饱和度获取方法，其特征在于，包括以下步骤：/n步骤一，获取平面孔隙网络模型图像，对该图像进行高斯滤波除去噪声后转化为灰度图像，对灰度图像进行坐标划分，并利用阈值分割方法将灰度图像转化像素尺寸为宽L

【技术特征摘要】
1.微米尺度平面孔隙网络内气液驱替图像的饱和度获取方法，其特征在于，包括以下步骤：
步骤一，获取平面孔隙网络模型图像，对该图像进行高斯滤波除去噪声后转化为灰度图像，对灰度图像进行坐标划分，并利用阈值分割方法将灰度图像转化像素尺寸为宽Lx和高Ly的黑白二值图像；
步骤二，裁剪黑白二值图像中入口与出口区域，保留中间孔隙网络区域，生成校验图像；
步骤三，统计校验图像中白色区域像素数，获得白色区域总面积，计算流动区域总面积Atotal；
步骤四，获取平面孔隙网络气相驱替形态图像，裁剪平面孔隙网络气相驱替形态图像入口与出口区域，并进行灰度阈值分割生成黑白二值待测图像；
步骤五，统计黑白二值待测图像中白色区域像素数，获得白色区域总面积，计算气相和液相流动区域面积；
步骤六，对比校验图像与黑白二值待测图像，获得气液两相流动区域驱替形态图像，并求取气相饱和度sg和液相饱和度sl以及驱替流量qg和ql。


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,yn)与其相邻的坐标点的像素值p(0,yn-1)和p(0,yn+1)分别进行“与”运算，搜索运算结果为0的坐标点Ai(0,yA)和Bi(0,yB)即为二值图像左边界上像素灰度变化点；
(2)将平面孔隙网络右边界上所有坐标点的像素值p(Lx,ym)与其相邻的坐标点p(Lx,ym-1)和p(Lx,ym+1)的像素点进行“与”运算，搜索运算结果为0的坐标点Co(Lx,yC)、Do(Lx,yD)即为二值图像右边界上像素灰度变化点；
(3)分别沿直线y＝yA和直线y＝yB向x轴正方向按步骤(1)的方法搜索获取左边界上像素灰度变化点A(xA,yA)、B(xB,yB)；
(4)分别沿直线y＝yC和直线y＝yD向x轴负方向按步骤(2)的方法搜索获取右边界上像素灰度变化点C(xC,yC)、D(xD,yD)；
(5)以直线A-B和直线C-D为界限，裁剪黑白二值图像，去除平面孔隙网络入口和出口区域，只保留中间孔隙网络区域，获得像素尺寸为宽(xD-xA)和高Ly或宽(xC-xB)和高Ly的校验图像。


5.根据权利要求1所述的微米尺度平面孔隙网络内气液驱替图像的饱和度...

【专利技术属性】
技术研发人员：文孟刚，李云，雷达，宋岩，
申请(专利权)人：西安交通大学，
类型：发明
国别省市：陕西;61