• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 28
    • 下载费用:30 金币  

    重庆时时彩后三杀一码技巧: 岩石三维图像裂缝提取方法.pdf

    关 键 词:
    岩石 三维 图像 裂缝 提取 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201710030343.6

    申请日:

    2017.01.17

    公开号:

    CN106803253A

    公开日:

    2017.06.06

    当前法律状态:

    实审

    有效性:

    审中

    法律详情: 实质审查的生效IPC(主分类):G06T 7/00申请日:20170117|||公开
    IPC分类号: G06T7/00(2017.01)I; G06T7/10(2017.01)I; G06T15/00(2011.01)I; G06K9/62 主分类号: G06T7/00
    申请人: 四川大学
    发明人: 滕奇志; 夏晨木; 何小海; 张耀; 吴晓红; 卿粼波
    地址: 610065 四川省成都市武侯区一环路南一段24号
    优先权:
    专利代理机构: 代理人:
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201710030343.6

    授权公告号:

    |||

    法律状态公告日:

    2017.06.30|||2017.06.06

    法律状态类型:

    实质审查的生效|||公开

    摘要

    本发明公开了一种岩石三维图像裂缝提取方法。该方法对三维岩石孔隙模型每个连通分量依次做表面重建、拉普拉斯网格平滑及网格简化操作;然后基于三角网格单元三角面积和单位法向量特征,利用K??均值聚类计算聚类中心;根据聚类中心重新将三角网格单元分类;对每个三角网格类中的三角网格单元构成的三维面结构执行孔洞修补,以获得封闭面结构集合;依据形状因子、面积比等参数判定当前封闭面结构集合是否为裂缝;对封闭面结构集合表面及内部体素点集执行形态学膨胀处理并与原始三维岩石孔隙模型连通分量的体素点集进行逻辑与操作,与操作得到的体素点集合即为裂缝。与传统方法相比,本发明可以正确识别并准确提取岩石裂缝。

    权利要求书

    1.岩石三维图像裂缝提取方法,其特征在于经过图像预处理后,算法包括以下步骤:
    步骤1:计算预处理得到的三角网格模型Aq中每个三角网格单元的单位法向量与三角面
    积,并调整单位法向量方向,其中q=1,2,3…k,k为三角网格模型个数;接着利用三角网格
    单元的单位法向量及三角面积特征对Aq的三角网格单元进行K-均值聚类;
    步骤2:根据步骤1得到的聚类中心将Aq的三角网格单元重新分类,分类得到多个三角网
    格类;
    步骤3:提取出步骤2中每个三角网格类的三角网格连通分量;
    步骤4:剔除步骤3得到的三角网格连通分量中的噪声成分,重新计算三角网格类的聚
    类中心,并将重新计算得到的聚类中心相邻的所有三角网格类合并为一类;接着对每个三
    角网格类分别独立执行步骤5到步骤10的操作;
    步骤5:步骤4得到的每个三角网格类中的所有三角网格单元组成一个三维面结构,对
    这个三维面结构执行三维模型孔洞修补操作,操作后可得到一个或多个封闭的三维面结
    构,这些三维面结构组成一个封闭面结构集合;
    步骤6:计算步骤5获得的封闭面结构集合的最小外接长方体;
    步骤7:计算步骤5获得的封闭面结构集合的形状因子;计算步骤6得到的最小外接长方
    体的形状因子;计算面积比参数;
    步骤8:利用步骤7得到的两个形状因子及面积比等参数,判定当前封闭面结构集合是
    否为裂缝;若判定结果为真,则继续执行后续操作;若判定结果为假,则当前封闭面结构集
    合不是裂缝;
    步骤9:提取步骤8中判定结果为真的封闭面结构集合表面及内部的体素点,对提取的
    体素点集作三维图像形态学膨胀处理;
    步骤10:将步骤9中已膨胀的体素点集与原始三维岩石孔隙模型连通分量的体素点集
    进行逻辑与操作,操作后得到的体素点集即为裂缝。
    2.根据权利要求1所述的方法,其特征在于步骤1中所述基于网格单位法向方向及网格
    三角面积特征的K-均值聚类,计算三角网格单元聚类中心,采用下述计算方法:
    K-均值聚类算法流程如下:
    ①确定初始聚类中心;
    ②将所有三角网格单元分配到距其最近的聚类中心所属三角网格类;
    ③更新聚类中心;
    ④反复执行流程②,③直到满足终止条件;
    K-均值聚类算法中相关计算公式如下所示:
    流程①中初始聚类中心Vi的计算公式如下所示,其中i=1,2,3:
    <mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>V</mi> <mn>1</mn> </msub> <mo>=</mo> <mi>&sigma;</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>V</mi> <mn>2</mn> </msub> <mo>=</mo> <mi>&sigma;</mi> <mo>&times;</mo> <mi>&xi;</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>V</mi> <mn>3</mn> </msub> <mo>=</mo> <mi>&sigma;</mi> <mo>&times;</mo> <mrow> <mo>(</mo> <mi>&sigma;</mi> <mo>&times;</mo> <mi>&xi;</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
    其中σ为Aq中三角面积最大的三角网格单元的单位法向量;ξ为Aq中三角面积最小的三
    角网格单元的单位法向量;
    令Δj为Aq中第j个三角网格单元,λj为Δj单位法向量,其中j∈[1,m],m为Aq含有的三角
    网格单元个数;为Aq中第i个三角网格类,其聚类中心为Vi,其中i=1,2,3;为第i个三
    角网格类中第j个三角网格单元,为单位法向量,为三角面积,其中j∈[1,n],n
    为第i个三角网格类中三角网格单元个数;
    流程②中三角网格单元Δj与各聚类中心Vi距离θij的计算公式如下:
    <mrow> <msub> <mi>&theta;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>&pi;</mi> <mo>-</mo> <msub> <mi>&phi;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>&phi;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>&gt;</mo> <mi>&pi;</mi> <mo>/</mo> <mn>2</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&phi;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> </mtd> <mtd> <mrow> <msub> <mi>&phi;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>&le;</mo> <mi>&pi;</mi> <mo>/</mo> <mn>2</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
    其中
    <mrow> <msub> <mi>cos&phi;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>&lambda;</mi> <mi>j</mi> </msub> <msub> <mi>V</mi> <mi>i</mi> </msub> </mrow> <mrow> <mo>|</mo> <msub> <mi>&lambda;</mi> <mi>j</mi> </msub> <mo>|</mo> <mo>|</mo> <msub> <mi>V</mi> <mi>i</mi> </msub> <mo>|</mo> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
    流程③中聚类中心计算公式如下:
    <mrow> <msub> <mi>V</mi> <mi>i</mi> </msub> <mo>=</mo> <mfrac> <msubsup> <mi>&beta;</mi> <mi>n</mi> <mi>i</mi> </msubsup> <mrow> <mo>|</mo> <msubsup> <mi>&beta;</mi> <mi>n</mi> <mi>i</mi> </msubsup> <mo>|</mo> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
    其中
    且Add定义为:
    <mrow> <mi>A</mi> <mi>d</mi> <mi>d</mi> <mrow> <mo>(</mo> <msub> <mi>J</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>J</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>J</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>J</mi> <mn>2</mn> </msub> </mrow> </mtd> <mtd> <mrow> <mo>|</mo> <msub> <mi>J</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>J</mi> <mn>2</mn> </msub> <mo>|</mo> <mo>&GreaterEqual;</mo> <mo>|</mo> <msub> <mi>J</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>J</mi> <mn>2</mn> </msub> <mo>|</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>J</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>J</mi> <mn>2</mn> </msub> </mrow> </mtd> <mtd> <mrow> <mo>|</mo> <msub> <mi>J</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>J</mi> <mn>2</mn> </msub> <mo>|</mo> <mo>&lt;</mo> <mo>|</mo> <msub> <mi>J</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>J</mi> <mn>2</mn> </msub> <mo>|</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
    <mrow> <msubsup> <mi>&omega;</mi> <mi>j</mi> <mi>i</mi> </msubsup> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mn>5</mn> </mtd> <mtd> <mrow> <msubsup> <mi>S</mi> <mi>j</mi> <mi>i</mi> </msubsup> <mo>&gt;</mo> <msubsup> <mi>S</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> <mi>i</mi> </msubsup> <mo>/</mo> <mn>2</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mrow> <msubsup> <mi>S</mi> <mi>j</mi> <mi>i</mi> </msubsup> <mo>&le;</mo> <msubsup> <mi>S</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> <mi>i</mi> </msubsup> <mo>/</mo> <mn>2</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
    <mrow> <msubsup> <mi>S</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> <mi>i</mi> </msubsup> <mo>=</mo> <mi>m</mi> <mi>a</mi> <mi>x</mi> <mo>{</mo> <msubsup> <mi>S</mi> <mi>j</mi> <mi>i</mi> </msubsup> <mo>|</mo> <mi>j</mi> <mo>&Element;</mo> <mo>&lsqb;</mo> <mn>1</mn> <mo>,</mo> <mi>n</mi> <mo>&rsqb;</mo> <mo>}</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>
    其中n为第i个三角网格类中三角网格单元个数;
    流程④中迭代终止条件为:迭代次数超过50次或者聚类中心变化幅度小于0.001;第K
    次迭代聚类中心变化幅度值的计算公式如下:
    <mrow> <msup> <mo>&part;</mo> <mi>k</mi> </msup> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>3</mn> </munderover> <mrow> <mo>(</mo> <mo>|</mo> <msub> <mi>x</mi> <msubsup> <mi>V</mi> <mi>i</mi> <mi>k</mi> </msubsup> </msub> <mo>-</mo> <msub> <mi>x</mi> <msubsup> <mi>V</mi> <mi>i</mi> <mrow> <mi>k</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> </msub> <mo>|</mo> <mo>+</mo> <mo>|</mo> <msub> <mi>y</mi> <msubsup> <mi>V</mi> <mi>i</mi> <mi>k</mi> </msubsup> </msub> <mo>-</mo> <msub> <mi>y</mi> <msubsup> <mi>V</mi> <mi>i</mi> <mrow> <mi>k</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> </msub> <mo>|</mo> <mo>+</mo> <mo>|</mo> <msub> <mi>z</mi> <msubsup> <mi>V</mi> <mi>i</mi> <mi>k</mi> </msubsup> </msub> <mo>-</mo> <msub> <mi>z</mi> <msubsup> <mi>V</mi> <mi>i</mi> <mrow> <mi>k</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> </msub> <mo>|</mo> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>
    其中为第k次迭代结束时聚类中心的三个坐标值。
    3.根据权利要求1所述的方法,其特征在于步骤2所述三角网格单元分类方法,采用下
    述计算方法:
    三角网格单元Δj与聚类中心Vi的距离为θij,且
    θkj=min{θij|i∈[1,2,3]} (10)
    若θkj<0.59359877则Δj属于第k三角网格类否则Δj不属于任何一类。
    4.根据权利要求1所述的方法,其特征在于步骤4所述三角网格连通分量噪声剔除及合
    并相近三角网格类,采用下述计算方法:
    三角网格类的三角网格单元三角面积和为三角网格连通分量的三角网格单元
    三角面积和为其中t∈[1,k],k是三角网格类的三角网格连通分量个数;计算公式如
    下:

    <mrow> <msub> <mi>S</mi> <msubsup> <mi>&eta;</mi> <mi>t</mi> <mi>i</mi> </msubsup> </msub> <mo>=</mo> <munder> <mo>&Sigma;</mo> <mrow> <msubsup> <mi>&Delta;</mi> <mi>j</mi> <mi>i</mi> </msubsup> <mo>&Element;</mo> <msubsup> <mi>&eta;</mi> <mi>t</mi> <mi>i</mi> </msubsup> </mrow> </munder> <msubsup> <mi>S</mi> <mi>j</mi> <mi>i</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow>
    其中为三角网格单元三角面积;
    对中三角网格连通分量若

    则将三角网格连通分量视为噪声干扰,并从三角网格类中剔除;
    设聚类中心Vi与聚类中心Vj的距离为χij,其计算公式如下:
    <mrow> <msub> <mi>&chi;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>&pi;</mi> <mo>-</mo> <msub> <mi>&phi;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>&phi;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>&gt;</mo> <mi>&pi;</mi> <mo>/</mo> <mn>2</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&phi;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> </mtd> <mtd> <mrow> <msub> <mi>&phi;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>&le;</mo> <mi>&pi;</mi> <mo>/</mo> <mn>2</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>14</mn> <mo>)</mo> </mrow> </mrow>
    其中
    <mrow> <msub> <mi>cos&phi;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>V</mi> <mi>i</mi> </msub> <msub> <mi>V</mi> <mi>j</mi> </msub> </mrow> <mrow> <mo>|</mo> <msub> <mi>V</mi> <mi>i</mi> </msub> <mo>|</mo> <mo>|</mo> <msub> <mi>V</mi> <mi>j</mi> </msub> <mo>|</mo> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>15</mn> <mo>)</mo> </mrow> </mrow>
    若χij<0.7853981,则视两个聚类中心所属三角网格类相邻,将三个三角网格类中所有
    相邻的类合并为一个类;合并后的三角网格类仍记为其中1≤i≤3。
    5.根据权利要求1所述的方法,其特征在于步骤8所述裂缝的判定,采用以下判定标准:
    当下列条件均成立时:
    <mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>&Phi;</mi> <mo>&Pi;</mo> </msub> <mo>&le;</mo> <mi>M</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>&Lambda;</mi> <mo>&GreaterEqual;</mo> <mi>E</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>&Phi;</mi> <mi>&psi;</mi> </msub> <mo>&le;</mo> <mn>0.05</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>16</mn> <mo>)</mo> </mrow> </mrow>
    或当以下条件均成立时:
    <mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>&Phi;</mi> <mo>&Pi;</mo> </msub> <mo>&le;</mo> <mn>0.04</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>&Lambda;</mi> <mo>&GreaterEqual;</mo> <mn>0.17</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>&Phi;</mi> <mi>&psi;</mi> </msub> <mo>&le;</mo> <mn>0.05</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>17</mn> <mo>)</mo> </mrow> </mrow>
    其中
    <mrow> <mi>M</mi> <mo>=</mo> <mn>0.6</mn> <mo>*</mo> <mfrac> <mrow> <msub> <mi>log</mi> <mn>2</mn> </msub> <msub> <mi>V</mi> <mo>&Pi;</mo> </msub> </mrow> <mn>20</mn> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>18</mn> <mo>)</mo> </mrow> </mrow>
    <mrow> <mi>E</mi> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mn>0.22</mn> </mtd> <mtd> <mrow> <mi>M</mi> <mo>&lt;</mo> <mn>0.45</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>0.17</mn> </mtd> <mtd> <mrow> <mi>M</mi> <mo>&GreaterEqual;</mo> <mn>0.45</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>19</mn> <mo>)</mo> </mrow> </mrow>
    <mrow> <msub> <mi>&Phi;</mi> <mi>&psi;</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mn>36</mn> <msup> <msub> <mi>&pi;V</mi> <mi>&psi;</mi> </msub> <mn>2</mn> </msup> </mrow> <mrow> <msup> <msub> <mi>S</mi> <mi>&psi;</mi> </msub> <mn>3</mn> </msup> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>20</mn> <mo>)</mo> </mrow> </mrow>
    <mrow> <msub> <mi>&Phi;</mi> <mo>&Pi;</mo> </msub> <mo>=</mo> <mfrac> <mrow> <mn>36</mn> <msup> <msub> <mi>&pi;V</mi> <mo>&Pi;</mo> </msub> <mn>2</mn> </msup> </mrow> <mrow> <msup> <msub> <mi>S</mi> <mo>&Pi;</mo> </msub> <mn>3</mn> </msup> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>21</mn> <mo>)</mo> </mrow> </mrow>
    <mrow> <mi>&Lambda;</mi> <mo>=</mo> <mfrac> <msub> <mi>S</mi> <mi>&psi;</mi> </msub> <msub> <mi>S</mi> <mo>&Pi;</mo> </msub> </mfrac> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>22</mn> <mo>)</mo> </mrow> </mrow>
    Vψ为封闭面结构集合ψ的体积,Sψ为ψ孔洞修补前所有三角网格单元三角面积之和,VΠ为
    封闭面结构集合ψ最小外接长方体Π的体积,SΠ为Π的表面积;
    则判定当前封闭面结构集合ψ表面及内部包含的体素点集δ属于裂缝的一部分;否则判
    定为非裂缝。
    6.根据权利要求1所述的方法,其特征在于步骤10所述裂缝的准确提取,使用以下提取
    方法:
    裂缝体素点集ζ计算公式为:
    ζ={p|p∈Hq∩δ'} (23)
    其中Hq为Aq对应的原始三维岩石孔隙模型连通分量的体素点集,p为任意体素点;δ'为
    步骤9中形态学膨胀处理后得到的体素点集。

    说明书

    岩石三维图像裂缝提取方法

    技术领域

    本发明涉及一种图像分割与识别技术,尤其涉及一种岩石三维图像裂缝提取方
    法,属于图像分割及识别技术领域。

    技术背景

    岩石在应力作用下产生机械性破坏,无明显位移的断裂构造叫裂缝。裂缝是一种
    特殊的孔隙。在石油地质领域,裂缝识别是油藏和储存非均质研究的重要内容,对油气开发
    和原油采收率的提高具有重要意义。

    本发明针对的是岩石CT序列图像的裂缝识别。本发明需先对岩石CT序列图中的每
    幅图像进行二值化,再利用二值化的序列图构建出三维岩石孔隙模型,并在三维空间识别
    出岩石裂缝。经由CT获取的岩石图像是一组连续的二维序列图,利用传统二维图像分割算
    法,可依次对序列图中每一幅CT图像作裂缝提取,进而便可得到裂缝的三维图像。但上述方
    法仅利用了二维图像信息,忽略了岩石裂缝在多张连续岩石CT图像之间的关联信息。这将
    导致CT图像呈狭长线条状而三维空间呈非裂缝形态的普通孔隙易被误判为裂缝,或CT图像
    呈非狭长线条状而三维空间呈裂缝形态的裂缝易被误判为非裂缝。

    目前已有的基于三维空间的裂缝识别方法是依据每个孤立三维孔隙的最小外接
    球半径与等效球半径的比值、形状因子、近似最小外接长方体最长边与最短边的比值等参
    数来识别岩石裂缝。但由于裂缝通常不是独立存在的,裂缝与普通孔隙之间关联紧密,故此
    方法对裂缝的误识率较高,如将与普通孔隙相连的裂缝判为非裂缝,或将与裂缝相连的普
    通孔隙判为裂缝。

    近年来,研究人员在不断地研究与探索三维岩石图像裂缝识别问题,但目前还未
    提出一种可对三维岩石图像裂缝进行高效准确识别的方法。本发明正是基于这种研究现
    状,提出一种可高效准确识别三维岩石图像裂缝的方法。

    发明内容

    本发明的目的在于克服现有技术存在的缺陷和不足,提供一种三维岩石图像裂缝
    的识别方法。

    本发明提供的三维岩石图像裂缝识别方法基于下述思想:裂缝是一种特殊的孔
    隙,常与普通孔隙相连接,故不能直接判定某个孔隙是否为裂缝。本发明对三维岩石孔隙模
    型的每个连通分量执行表面重建、拉普拉斯网格平滑、网格简化、网格单元单位法向量计算
    等一系列操作。由于裂缝存在面状特征,位于裂缝表面的三角网格单元被简化的概率高。三
    角网格单元被简化后网格三角面积将变大,且网格单位法向量方向趋于统一。根据网格三
    角面积和网格单位法向量方向特征,对三角网格单元执行K-均值聚类,根据聚类后得到的
    聚类中心将三角网格类重新划分,位于同一裂缝表面的三角网格单元将聚于同一类中。进
    而通过形状因子等参数即可判定每个三角网格类构成的三维空间结构是否具有裂缝特征。
    为了提高裂缝识别的精确度,对具有裂缝特征的三维空间结构执行三维图像形态学膨胀操
    作,之后将膨胀得到的三维空间结构的体素点集与原始三维岩石孔隙模型连通分量的体素
    点集进行逻辑与操作,与操作的结果即为岩石裂缝。

    本发明提供的三维岩石图像裂缝的识别方法,包括以下操作步骤:

    预处理:对一组待识别的岩石二值序列图构建其相应的三维岩石孔隙模型,并提
    取模型所有连通分量;剔除其中体素点过少的连通分量;使用移动立方体等值面提取
    (Marching Cube)算法对每个连通分量进行表面重建,得到每个连通分量相应的三角网格
    模型Aq,其中q=1,2,3…k,k为三角网格模型个数;对得到的三角网格模型作拉普拉斯网格
    平滑处理,并使用网格简化算法对三角网格模型进行网格简化;接着对每一个Aq独立执行
    下述步骤1到步骤10的操作;

    步骤1:计算预处理得到的三角网格模型Aq中每个三角网格单元的单位法向量与
    三角面积,并调整单位法向量方向;接着利用三角网格单元的单位法向量及三角面积特征
    对Aq的三角网格单元进行K-均值聚类;

    步骤2:根据步骤1得到的聚类中心将Aq的三角网格单元重新分类,分类得到多个
    三角网格类;

    步骤3:提取出步骤2中每个三角网格类的三角网格连通分量;

    步骤4:剔除步骤3得到的三角网格连通分量中的噪声成分,重新计算三角网格类
    的聚类中心,并将重新计算得到的聚类中心相邻的所有三角网格类合并为一类;接着对每
    个三角网格类分别独立执行步骤5到步骤10的操作;

    步骤5:步骤4得到的每个三角网格类中的所有三角网格单元组成一个三维面结
    构,对这个三维面结构执行三维模型孔洞修补操作,操作后可得到一个或多个封闭的三维
    面结构,这些三维面结构组成一个封闭面结构集合;

    步骤6:计算步骤5获得的封闭面结构集合的最小外接长方体;

    步骤7:计算步骤5获得的封闭面结构集合的形状因子;计算步骤6得到的最小外接
    长方体的形状因子;计算面积比参数;

    步骤8:利用步骤7得到的两个形状因子及面积比等参数,判定当前封闭面结构集
    合是否为裂缝。若判定结果为真,则继续执行后续操作;若判定结果为假,则当前封闭面结
    构集合不是裂缝;

    步骤9:提取步骤8中判定结果为真的封闭面结构集合表面及内部的体素点,对提
    取的体素点集作三维图像形态学膨胀处理;

    步骤10:将步骤9中已膨胀的体素点集与原始三维岩石孔隙模型连通分量的体素
    点集进行逻辑与操作,操作后得到的体素点集即为裂缝。

    上述技术方案中,预处理中所述剔除体素点过少的连通分量,优先采用的体素点
    个数阈值为100。

    上述技术方案中,预处理中所述对三角网格模型执行的拉普拉斯网格平滑处理,
    优先采用以下参数设置:拉普拉斯网格平滑指数为0.2,平滑次数为100。

    上述技术方案中,预处理中所述三角网格简化处理,优先采用网格抽取
    (Decimation of Triangle Meshes)算法,网格抽取比例控制在0.9,即移除90%的三角网
    格单元。

    上述技术方案中,步骤1中所述对三角网格单元单位法向量方向的调整,优先采用
    以下调整方式:

    设计算得到单位法向量为(x,y,z),若z<0时,则调整单位法向量为(-x,-y,-z)。

    上述技术方案中,步骤1中所述对三角网格单元的K-均值聚类,计算聚类中心,优
    先采用下述方式进行计算:

    K-均值聚类算法流程如下:

    (1)确定初始聚类中心;

    (2)将所有三角网格单元分配到距其最近的聚类中心所属三角网格类;

    (3)更新聚类中心;

    (4)反复执行(2)(3)步直到满足终止条件。

    K-均值聚类算法中相关计算公式如下所示:

    流程(1)中初始聚类中心Vi的计算公式如下所示,其中i=1,2,3:


    其中σ为Aq中三角面积最大的三角网格单元的单位法向量;ξ为Aq中三角面积最小
    的三角网格单元的单位法向量。

    令Δj为Aq中第j个三角网格单元,λj为Δj单位法向量,其中j∈[1,m],m为Aq含有的
    三角网格单元个数;为Aq中第i个三角网格类,其聚类中心为Vi,其中i=1,2,3;为第i
    个三角网格类中第j个三角网格单元,为单位法向量,为三角面积,其中j∈[1,
    n],n为第i个三角网格类中三角网格单元个数。

    流程(2)中三角网格单元Δj与各聚类中心Vi距离θij的计算公式如下:


    其中


    流程(3)中聚类中心计算公式如下:


    其中


    Add定义为:




    其中n为第i个三角网格类中三角网格单元个数。

    流程(4)中迭代终止条件为:迭代次数超过50次或聚类中心变化幅度小于0.001;

    第K次迭代聚类中心变化幅度值的计算公式如下:


    其中为第k次迭代结束时聚类中心的三个坐标值。

    上述技术方案中,步骤2中所述三角网格单元的重新分类,优先采用下述方法进行
    分类:

    三角网格单元Δj与聚类中心Vi的距离为θij,且

    θkj=min{θij|i∈[1,2,3]} (10)

    若θkj<0.59359877,则Δj属于第k三角网格类否则Δj不属于任何一类。

    上述技术方案中,步骤3中所述提取三角网格类的三角网格连通分量,优先采用下
    述方法:

    若三角网格单元Δi的三角顶点与三角网格单元Δj的三角顶点有至少一个顶点相
    同,则Δi与Δj空间相邻。每个三角网格类中,按空间相邻关系提取三角网格连通分量
    其中t∈[1,k],k是三角网格类的三角网格连通分量个数。三角网格类与三角网格
    连通分量的关系如下:


    上述技术方案中,步骤4中所述对三角网格连通分量噪声剔除的方法,优先采用下
    述方法:

    三角网格类的三角网格单元三角面积之和为三角网格连通分量的三角
    网格单元三角面积之和为其中t∈[1,k],k是三角网格类的三角网格连通分量个数。
    计算公式如下:



    其中为三角网格单元三角面积。

    对中所有三角网格连通分量若


    则将三角网格连通分量视为噪声干扰,并从三角网格类中剔除。

    上述技术方案中,步骤4中所述合并相近三角网格类的方法,优先采用下述方法进
    行合并:

    设聚类中心Vi与聚类中心Vj的距离为χij,其计算公式如下:


    其中


    若χij<0.7853981,则视两个聚类中心所属三角网格类相邻。将三个三角网格类中
    所有相邻的类合并为一个类。合并后的三角网格类仍记为其中1≤i≤3。

    上述技术方案中,步骤6中所述封闭面结构集合的最小外接长方体计算,优先采用
    以下文献中所述的方法计算:三维岩心图像裂缝自动识别[期刊论文]邓知秋,滕奇志-《计
    算机与数字工程》2013年1期。

    上述技术方案中,步骤7中所述形状因子、面积比的计算,优先采用下述方法:

    ①封闭面结构集合ψ的形状因子Φψ计算公式如下:


    ②最小外接长方体形Π的形状因子ΦΠ计算公式如下:


    ③面积比Λ计算公式如下:


    其中Vψ为封闭面结构集合ψ的体积;Sψ为ψ孔洞修补前所有三角网格单元三角面积
    之和;VΠ为封闭面结构集合最小外接长方体Π的体积;SΠ为Π的表面积。

    上述技术方案中,步骤8中所述裂缝的判定标准,优先使用以下判定标准:

    当下列条件均成立时:


    或当以下条件均成立时:


    其中



    则判定当前封闭面结构集合ψ表面及内部的体素点集δ属于裂缝的一部分,并执行
    后续步骤;否则判定为非裂缝。

    上述技术方案中,步骤9中所述三维图像形态学膨胀操作,优先采用以下方法:


    其中B为中心位于(0,0,0)、尺寸为3*3的立方体;z为空间任意点(z1,z2,z3)

    (B)z={c|c=b+z,b∈B} (25)

    上述技术方案中,步骤10中所述逻辑与操作,优先采用以下方法:

    裂缝体素点集ζ计算公式为:

    ζ={p|p∈Hq∩δ'} (26)

    其中Hq为Aq对应的原始三维岩石孔隙模型连通分量的体素点集,p为任意体素点。
    δ'为步骤9中形态学膨胀处理后得到的体素点集。

    本发明与现有技术相比,所具有的优点及有益的技术效果如下:

    本发明所述的三维岩石图像裂缝识别方法,以实际的岩石二值序列图构建三维岩
    石孔隙模型,使用本发明所述的方法对其裂缝进行识别,并对识别结果进行观察,从而验证
    了本发明所用方法的可靠性与实用性。本发明解决了传统方法中裂缝识别不正确以及裂缝
    提取不准确的问题,本发明可以正确识别并精确提取岩石裂缝。

    本发明受国家自然科学基金“岩石微观非均质结构三维图像重建及分辨率提升技
    术研究(61372174)”资助。

    附图说明

    图1-1是本发明实施例中岩石二值序列图之一;

    图1-2是本发明实施例中由岩石二值序列图构造的三维岩石孔隙模型;

    图2-1是本发明实施例中对图1-2中模型去噪后的三维岩石孔隙模型;

    图2-2是本发明实施例中对图2-1中模型进行表面重建的结果;

    图3-1是本发明实施例中对图2-2中模型执行拉普拉斯平滑后的局部放大图片;

    图3-2是本发明实施例中对图2-2中模型执行拉普拉斯平滑后的局部放大图片;

    图4-1是本发明实施例中对图3-1中模型执行网格抽取后的局部放大图片;

    图4-2是本发明实施例中对图3-1中模型执行网格抽取后的局部放大图片;

    图5是本发明实施例中对图4-1中模型计算网格单位法向量的局部示意图;

    图6-1是本发明实施例中对图5中模型的三角网格单元K-均值聚类后重新分类的
    结果图;

    图6-2是本发明实施例中对图5中模型的三角网格单元K-均值聚类后重新分类的
    结果图;

    图7-1是本发明实施例中对图6-1中分类后的三角网格单元执行连通分量标记的
    结果图;

    图7-2是本发明实施例中对图6-1中分类后的三角网格单元执行连通分量标记的
    结果图;

    图7-3是本发明实施例中对图6-1中分类后的三角网格单元执行连通分量标记的
    结果图;

    图8-1是本发明实施例中三维岩石孔隙模型裂缝识别的效果图;

    图8-2是本发明实施例中图8-1中第一条裂缝的最小外接长方体示意图;

    图8-3是本发明实施例中图8-1中第二条裂缝的最小外接长方体示意图;

    图8-4是本发明实施例中图8-1中第三条裂缝的最小外接长方体示意图;

    图8-5是本发明实施例中三维岩石孔隙模型裂缝识别的效果图;

    图8-6是本发明实施例中图8-5中裂缝与原始孔隙同时显示的效果图;

    图8-7是本发明实施例中三维岩石孔隙模型裂缝识别的效果图;

    图8-8是本发明实施例中图8-7中裂缝与原始孔隙同时显示的效果图;

    图8-9是本发明实施例中三维岩石孔隙模型裂缝识别的效果图;

    图8-10是本发明实施例中图8-9中裂缝与原始孔隙同时显示的效果图;

    图9是本发明所述的方法流程图。

    具体实施方式

    下面通过具体实施例并结合附图对本发明作进一步详细说明,有必要在此指出的
    是所述实施例只是用于对本发明的进一步描述,但并不意味着是对本发明?;し段У娜魏?br />限定。

    实施例:为了使本发明所述识别方法更加便于理解和接近于真实应用,下面从岩
    石二值序列图像开始到三维岩石孔隙模型裂缝识别完成为止,对整个操作流程进行说明。
    具体操作步骤如下:

    (1)准备一组待识别的岩石二值序列图,图片数量为256,格式为BMP,尺寸256×
    256,其中代表如图1-1所示;由岩石二值序列图构建的三维岩石孔隙模型,如图1-2所示,图
    中灰色部分为岩石孔隙。

    (2)提取三维岩石孔隙模型的所有连通分量,剔除其中体素点个数小于100的连通
    分量,去噪结果如图2-1所示。用移动立方体等值面提取(Marching Cube)算法对每个连通
    分量进行表面重建,得到每个连通分量相应的三角网格模型Aq,其中q=1,2,3……k,k为三
    角网格模型个数。表面重建结果如图2-2所示。

    (3)对每个表面重建结果Aq执行平滑指数0.2,平滑次数100的拉普拉斯网格平滑
    处理;网格平滑效果如模型局部放大图3-1和3-2所示。

    (4)使用三角网格抽取算法对平滑后的三角网格模型执行网格简化处理,网格抽
    取比例控制在0.9。网格简化效果如模型局部放大图4-1和4-2所示,图中直线段为三角网格
    的边界。

    (5)计算所有三角网格单元的单位法向量,设计算得到单位法向量为(x,y,z),若z
    <0时,则调整单位法向量为(-x,-y,-z)。单位法向量计算结果如图5所示,图中三角网格单
    元上的箭头为其单位法向量。

    (6)利用三角网格单元的单位法向量及三角网格单元三角面积对Aq中的三角网格
    单元进行K-均值聚类,得到三个聚类中心V1、V2、V3。

    K-均值聚类算法流程如下:

    1)求Aq中三角面积最大的三角网格单元的单位法向量σ及Aq中三角面积最小的三
    角网格单元的单位法向量ξ。则令σ、σ×ξ、σ×(σ×ξ)分别为初始的三个聚类中心;

    2)计算三角网格单元Δj与各聚类中心Vi距离θij,计算公式为:


    其中


    其中λj为三角网格Δj的单位法向量。

    将三角网格单元分配到距离其最近的聚类中心所属三角网格类。

    3)令Δj为Aq中第j个三角网格单元,λj为Δj单位法向量,其中j∈[1,m],m为Aq含有
    的三角网格单元个数;为Aq中第i个三角网格类,其聚类中心为Vi,其中i=1,2,3;为第
    i个三角网格类中第j个三角网格单元,为单位法向量,为三角面积,其中j∈
    [1,n],n为第i个三角网格类中三角网格单元个数;

    重新计算聚类中心Vi,计算公式如下:


    其中


    Add定义为:




    其中n为第i个三角网格类中三角网格单元个数。

    4)计算第k次迭代结束时,聚类中心变化幅度值


    其中为第k次迭代结束时聚类中心的三个坐标值。

    若迭代次数k超过50或者聚类中心变化幅度小于0.001则结束迭代,否则返回流
    程2)。

    (7)计算每个三角网格单元Δj与聚类中心Vi的距离θij,且

    θkj=min{θij|i∈[1,2,3]} (9)

    若θkj<0.59359877则将Δj归为第k三角网格类否则Δj不属于任何一类。

    (8)提取三角网格类的三角网格连通分量:若三角网格单元Δi的三角顶点与
    Δj的三角顶点至少有一个顶点相同,则Δi与Δj空间相邻。每个三角网格类中,按空间相
    邻关系提取三角网格连通分量其中t∈[1,k],k是三角网格类的三角网格连通分量
    个数。三角网格连通分量提取结果如图7-1、7-2和7-3所示。

    (9)计算每个三角网格类中所有三角网格单元的三角面积和计算三角网格
    连通分量中所有三角网格单元的三角面积和其中t∈[1,k],k是三角网格类的三
    角网格连通分量个数。计算公式分别为:



    其中为三角网格单元的三角面积。

    对中所有三角网格连通分量若


    则将此三角网格连通分量视为噪声干扰,并从三角网格类中剔除。

    剔除噪声后,重新计算每个类的聚类中心,并合并相近的类。计算聚类中心Vi与聚
    类中心Vj的距离χij,计算公式如下:


    其中


    若χij<0.7853981,则视两个聚类中心所属三角网格类相邻。将三个三角网格类中
    所有相邻的类合并为一个类。合并后的三角网格类仍记为其中1≤i≤3。相邻类的合并
    结果如图6-1、6-2所示,图中不同灰度分别表示一个类,白色三角网格不属于任何一类。

    (10)三维岩石孔隙连通分量Aq中的每个三角网格类中的所有三角网格单元组成
    一个三维面结构。对这个三维面结构执行三维模型孔洞修补操作,操作后得到一个或多个
    封闭的三维面结构,这些三维面结构组成一个封闭面结构集合。

    (11)按如下文献计算封闭面结构集合的最小外接长方体:三维图像裂缝自动识别
    [期刊论文]邓知秋,滕奇志-《计算机与数字工程》2013年1期。

    (12)计算封闭面结构集合ψ的形状因子Φψ:


    计算最小外接长方体形Π的形状因子ΦΠ:


    计算面积比Λ:


    其中Vψ为封闭面结构集合ψ的体积;Sψ为ψ孔洞修补前的所有三角网格单元三角面
    积之和;VΠ为封闭面结构集合最小外接长方体Π的体积;SΠ为Π的表面积。

    (13)当下列条件均成立时:


    或者当以下条件均成立时:


    其中



    则判定当前封闭面结构集合ψ表面及内部包含的体素点集属于裂缝的一部分,并
    执行下一步骤;否则判定为非裂缝。

    (14)提取(13)中判定为真的封闭面结构集合的表面以及内部的所有体素点,对提
    取得到的体素点集进行三维图像形态学膨胀处理。

    (15):将(14)中膨胀得到的体素点集与原始三维岩石孔隙模型连通分量的体素点
    集进行逻辑与操作,操作后得到的体素点集即为裂缝。与操作结果如图8-1至8-10所示:

    图8-1为提取的裂缝示意图,图中共有3个裂缝;

    图8-2至8-4为每个裂缝的最小外接长方体示意图;

    图8-5至8-10为裂缝不同方位显示图,以及将其与原始孔隙同时显示的效果图。

    本实施例中,以实际的岩石二值序列图构建三维模型并使用本发明的方法对其进
    行裂缝识别。对识别结果进行观察,从而验证了本发明识别方法的可靠性与准确性。

    本发明上述实施例只是本发明的优选实施例,并不是对本发明所述技术方案的任
    何限制,只要是不经过创造性劳动即在上述实施例的基础上实现的技术方案,均应视为落
    入本发明所?;さ姆段?。

    关于本文
    本文标题:岩石三维图像裂缝提取方法.pdf
    链接地址://www.4mum.com.cn/p-6001245.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
  • 幸运飞艇计划软件下载手机版 怎样买3d稳赚不赔 北京pk10 6码2期计划 重庆欢乐生肖免费计划 pk10定位独胆技巧 买福彩3d稳赚不赔的绝密方法 3d彩票六码复式是多少钱 倍投资长龙 江苏时时诈骗案 双色球开奖各省数据 免费计划软件手机版下载 扎金花技巧顺口溜 重庆时时五星彩走势图360 pk10免费永久计划软件 北京pk赛车5码稳赚技巧 3分pk拾是真还是假的