[發(fā)明專利]一種考慮孔喉末端的孔隙喉道識(shí)別方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201410527749.1
• 申請(qǐng)日: 2014-10-09
• 公開(kāi)（公告）號(hào): CN104778678A
• 公開(kāi)（公告）日: 2015-07-15
• 發(fā)明（設(shè)計(jì)）人: 侯健;韋貝;于波;杜慶軍;劉永革;周康;陸努;夏志增
• 申請(qǐng)（專利權(quán)）人: 中國(guó)石油大學(xué)（華東）
• 主分類號(hào): G06T7/00
• 分類號(hào): G06T7/00;G06K9/00;G06K9/46;G01N15/08
• 代理公司: 暫無(wú)信息
• 代理人: 暫無(wú)信息

1.一種考慮孔喉末端的孔隙喉道識(shí)別方法，其特征在于以下步驟：

(1)對(duì)圖像進(jìn)行二值化處理，區(qū)分巖石顆粒與孔喉空間；

(2)以巖石顆粒為前景像素對(duì)孔喉空間進(jìn)行歐式長(zhǎng)度變換，變換計(jì)算公式為：

Di=min(D(p(xi,yi),q(xj,yj)))=min((xi-xj)2+(yi-yj)2)]]>

其中p為孔喉空間內(nèi)一點(diǎn)，q為p點(diǎn)周圍任意的巖石顆粒點(diǎn)；

(3)由基于長(zhǎng)度變換的轉(zhuǎn)折點(diǎn)算法提取圖像中線，此時(shí)的中線未考慮末端；

(4)基于轉(zhuǎn)折點(diǎn)中線，進(jìn)行不考慮末端時(shí)的孔喉識(shí)別；

(5)利用燃燒算法提取中線，查找燃燒中線尾端端點(diǎn)，確定末端所在區(qū)域；

(6)借助燃燒中線尾端端點(diǎn)，利用轉(zhuǎn)折點(diǎn)算法提取末端中線；

(7)將步驟(4)中的轉(zhuǎn)折點(diǎn)中線與步驟(6)中的末端轉(zhuǎn)折點(diǎn)中線合并，對(duì)末端區(qū)域進(jìn)行孔喉識(shí)別，然后將末端識(shí)別結(jié)果與步驟(4)中未考慮末端區(qū)域的孔喉識(shí)別結(jié)果合并，得到考慮末端的孔喉識(shí)別結(jié)果。

2.如權(quán)利要求1所述的一種可以識(shí)別孔喉末端的孔隙喉道識(shí)別方法，其特征在于：

其所述步驟(4)中，孔喉識(shí)別的步驟如下：

①刪除位于孔隙內(nèi)部的轉(zhuǎn)折點(diǎn)中線路徑，將剩余中線路徑所對(duì)應(yīng)像素按順序編號(hào)，并以每個(gè)像素的長(zhǎng)度變換結(jié)果作為此像素對(duì)應(yīng)的孔隙半徑；

②以像素編號(hào)為橫坐標(biāo)，對(duì)應(yīng)的孔隙半徑為縱坐標(biāo)，做出每條路徑的沿程曲線；

③將沿程曲線看作信號(hào)，選擇離散Dmey波為小波基函數(shù)，對(duì)信號(hào)進(jìn)行二尺度分解，選擇合理閾值截取小波系數(shù)，然后重構(gòu)信號(hào)完成小波消噪，閾值選擇公式為：

其中σ為噪聲標(biāo)準(zhǔn)方差，N是信號(hào)長(zhǎng)度；

④在去噪后的沿程曲線上查找最靠近兩端的兩個(gè)局部極小值，將其對(duì)應(yīng)的像素作為孔隙和喉道的分界點(diǎn)即喉道的瓶頸，在瓶頸對(duì)應(yīng)的位置畫(huà)出孔喉識(shí)別線；

⑤孔喉識(shí)別線和巖石邊界將非基質(zhì)空間分成了若干封閉區(qū)域，將包含中線結(jié)點(diǎn)的區(qū)域作為孔隙，其余部分作為喉道。

3.如權(quán)利要求1所述的一種可以識(shí)別孔喉末端的孔隙喉道識(shí)別方法，其特征在于：

其所述步驟(6)中，提取末端中線步驟如下：

①?gòu)牟襟E(4)識(shí)別后的孔隙或喉道中，提取包含燃燒中線端點(diǎn)像素的喉道或孔隙區(qū)域即為末端區(qū)域；

②以燃燒中線的末端像素為原點(diǎn)，沿正交方向?qū)⑻崛〉哪┒藚^(qū)域切分為四個(gè)子圖，應(yīng)用轉(zhuǎn)折點(diǎn)算法分別提取四個(gè)子圖的中線；

③將四個(gè)子圖的中線合并在一起并清理分支，所得結(jié)果即為轉(zhuǎn)折點(diǎn)中線在末端處的分支。