申请/专利权人：中海石油(中国)有限公司;中海石油(中国)有限公司北京研究中心

申请日：2021-06-11

公开（公告）日：2024-03-08

公开（公告）号：CN113362312B

主分类号：G06F17/00

分类号：G06F17/00;G06T7/00

优先权：

专利状态码：有效-授权

法律状态：2024.03.08#授权;2021.09.24#实质审查的生效;2021.09.07#公开

摘要：本发明公开了一种储层不连续界限方位及宽度的检测方法，包括：计算储层不连续界限的二值化界限；以所述储层不连续界限的二值化界限作为掩模，通过四连通标记法获得储层不连续界限骨架，构建全方位等间隔方位算子及全方位等间隔宽度算子；计算储层不连续界限的宽度属性和方位属性。本发明公开的储层不连续界限方位及宽度的检测方法，能够通过检测储层不连续界限的宽度信息来评价其等级，检测储层不连续界限的方位信息来判断其与砂体物源之间的关系，从而为储层不连续界限阻渗能力的综合评价提供依据。

主权项：1.一种储层不连续界限方位及宽度的检测方法，其特征在于，包括以下步骤：步骤S1：输入储层不连续界限平面属性，计算全局最优二值化阈值，根据全局最优二值化阈值计算储层不连续界限的二值化界限，具体包括以下步骤：步骤S11：对储层不连续界限平面属性进行预处理，将其初始化为256级灰度数据，并统计归一化直方图；步骤S12：以最大类间方差法计算全局最优二值化阈值T，具体包括以下步骤：步骤S121：计算类间方差Tm，其表达式为：Tm＝p0×p1×u0-u1×u0-u11式中，m表示前景与背景的分割灰度，m∈[0,255]；p0表示前景像素点数所占比例；p1表示背景像素点数所占比例；u0表示前景像素点的平均灰度；u1表示背景像素点的平均灰度；步骤S122：计算全局最优二值化阈值T，其表达式为：T＝maxTm2式中，T为全局最优二值化阈值；Tm为类间方差；步骤S13：根据全局最优二值化阈值T计算表征储层不连续界限的二值化界限，具体包括以下步骤：步骤S131：建立储层不连续界限的直角坐标系，在储层不连续界限内选定某一点A，其横坐标为x，纵坐标为y，A点的属性值Bx,y的计算公式为： 其中，fx,y表示储层不连续界限平面属性中某点坐标x,y的属性值，x,y表示直角坐标系的坐标；步骤S132：重复所述步骤S131，将储层不连续界限的所有点遍历，建成二值化界限的数据体，得到所有点的储层不连续界限的二值化界限；步骤S2：以所述储层不连续界限的二值化界限作为掩模，通过四连通标记法获得储层不连续界限骨架，具体包括以下步骤：步骤S21：建立四连通标记算子S，具体如下： 其中，四连通标记算子S中0表示不参加运算的数据点；四连通标记算子S中1表示参加运算的数据点；四连通标记算子S中2表示标记的数据点；步骤S22：利用四连通标记算子S对所述储层不连续界限的二值化界限进行卷积运算，求取所述储层不连续界限的二值化界限中每一点的标记值Qx,y，其计算表达式为：Qx,y＝Bx,y*S4其中，“*”表示卷积运算；步骤S23：根据所述储层不连续界限的二值化界限中每一点的标记值Qx,y计算储层不连续界限骨架中每一点的属性值Kx,y，其计算表达式为： 步骤S23：将储层不连续界限骨架中每一点的属性值Kx,y组建成数据体，形成储层不连续界限骨架；步骤S3：构建全方位等间隔方位算子di及全方位等间隔宽度算子wj，具体包括以下步骤：步骤S31：构建全方位等间隔方位算子di，其中，全方位等间隔方位算子di表示在0～360°中等间隔方向的方位算子，i表示全方位等间隔方位算子di的方向，i∈[1,8]，所述全方位等间隔方位算子di具体包括d1、d2、d3、d4、d5、d6、d7和d8，具体形式分别表示如下： 步骤S32：构建全方位等间隔宽度算子wj，且该全方位等间隔宽度算子wj与全方位等间隔方位算子di正交，其中，j表示全方位等间隔宽度算子wj的方向，j∈[1,8]，所述全方位等间隔的宽度算子wj具体包括w1、w2、w3、w4、w5、w6、w7和w8，具体形式分别表示如下： 其中，w1、w2、w3、w4、w5、w6、w7和w8中的2、3和4表示参加运算数据点的不同权值；步骤S4：计算储层不连续界限的宽度属性Wx,y和方位属性Dx,y，具体包括以下步骤：步骤S41：计算储层不连续界限的方位属性D，具体包括以下步骤：步骤S411：在储层不连续界限骨架中任选一点P，其横坐标为x，纵坐标为y，P点的属性值Kx,y，则该点处的方位属性Dx,y计算方式如下：Dx,y＝maxKx,y*di7此时，为了计算储层不连续界限的宽度属性W，设置方位属性标记矩阵N，则P点的方位属性标记值Nx,y的表达式为：Nx,y＝i'8其中，i'表示P点时Dx,y对应算子的方向；步骤S412：重复所述步骤S411，将储层不连续界限骨架的所有点遍历一次，建成储层不连续界限的方位属性的数据体形成所有点的储层不连续界限的方位属性D，以及建成储层不连续界限的方位属性标记矩阵的数据体得到所有点的储层不连续界限的方位属性标记矩阵N；步骤S42：计算储层不连续界限的宽度属性W，具体包括以下步骤：步骤S421：在储层不连续界限的方位属性标记矩阵N中选取一点O，其横坐标为x，纵坐标为y，O点的标记值为方位属性标记值Nx,y，此时O点的全方位等间隔宽度算子wj的算子方向j’的计算表达式为：j'＝Nx,y9步骤S422：设O点处的全方位等间隔宽度算子为wj'，根据O点处的全方位等间隔宽度算子和储层不连续界限的二值化界限计算O点处储层不连续界限的宽度属性Wx,y，其计算表达式为：Wx,y＝Bx,y*wj'10步骤S423：重复上述步骤S421和S422，将储层不连续界限的方位属性标记矩阵N中每个点遍历一次，得到每个点的储层不连续界限宽度属性以形成储层不连续界限宽度属性W。

全文数据：

权利要求：

百度查询： 中海石油(中国)有限公司;中海石油(中国)有限公司北京研究中心 一种储层不连续界限方位及宽度的检测方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD