[發(fā)明專(zhuān)利]一種油藏地質(zhì)剖面圖矢量化方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201710444737.6
• 申請(qǐng)日: 2017-06-13
• 公開(kāi)（公告）號(hào): CN107273608B
• 公開(kāi)（公告）日: 2020-06-30
• 發(fā)明（設(shè)計(jì)）人: 孫歧峰;段友祥;李洪強(qiáng);李寧寧
• 申請(qǐng)（專(zhuān)利權(quán)）人: 中國(guó)石油大學(xué)（華東）
• 主分類(lèi)號(hào): G06F30/20
• 分類(lèi)號(hào): G06F30/20;G06K9/62;G06T7/13;G06T7/136
• 代理公司: 暫無(wú)信息
• 代理人: 暫無(wú)信息

1.一種油藏地質(zhì)剖面圖矢量化方法，其特征在于包括以下步驟：

步驟S1：運(yùn)用MATLAB軟件，通過(guò)將Direct數(shù)字化油藏表征軟件系統(tǒng)輸出的地質(zhì)剖面圖作為輸入，將該油藏地質(zhì)剖面圖根據(jù)本文的矢量化方法獲得矢量化結(jié)果；

步驟S2：對(duì)步驟S1中的輸入進(jìn)行預(yù)處理，采用高斯濾波對(duì)圖像進(jìn)行平滑濾波，然后在平滑的基礎(chǔ)上采用拉普拉斯銳化方法來(lái)增強(qiáng)圖像；

步驟S3：對(duì)預(yù)處理后的油藏地質(zhì)剖面圖由RGB顏色空間模型轉(zhuǎn)為L(zhǎng)ab顏色空間模型，再進(jìn)行顏色聚類(lèi)分割，得出初始油藏地質(zhì)剖面圖；

步驟S4：采用改進(jìn)的自適應(yīng)Canny邊緣檢測(cè)算法對(duì)所述初始油藏地質(zhì)剖面圖進(jìn)行邊緣檢測(cè)，獲取所述初始油藏地質(zhì)剖面圖的邊緣輪廓；

步驟S5：對(duì)所述初始油藏地質(zhì)剖面圖的邊緣輪廓進(jìn)行矢量化，得到目標(biāo)油藏地質(zhì)剖面圖；

步驟S6：輸出矢量化結(jié)果；

其中，步驟S4包括以下步驟：

步驟S41：首先利用一維高斯函數(shù)，對(duì)圖像進(jìn)行低通平滑濾波；

步驟S42：用Sobel算子對(duì)圖像進(jìn)行卷積求像素點(diǎn)的梯度，采用X軸、45°和135°方向上的一階梯度模板，模板如下所示：

；

三個(gè)方向上的一階梯度分量分別為G_X(i,j)、G_45°(i,j)、G_135°(i,j)，由三個(gè)對(duì)應(yīng)的一階梯度模板與圖像進(jìn)行卷積得到，其中(i,j)表示圖像中i行j列的像素點(diǎn)，每個(gè)像素點(diǎn)的梯度幅值為：梯度方向?yàn)椋篏_Y(i,j)為Y軸方向的一階梯度分量；

步驟S43：依據(jù)大小為n行m列的梯度圖像M，計(jì)算以(i,j)為起點(diǎn)，大小為(i+a)行(j+a)列的窗口的梯度的累加和sum,avg為點(diǎn)(i,j)在該窗口的梯度累加和的平均值，則其中i,j分別表示圖像i行j列的像素點(diǎn),i∈[0，n-1-a],j∈[0，m-1-a]，其中a為取某一像素點(diǎn)的鄰域窗口的大小，鄰域窗口的移動(dòng)步長(zhǎng)設(shè)為step；

步驟S44：矩陣Matrix為梯度圖像M中每一點(diǎn)的梯度平均值累加和矩陣，其中每一點(diǎn)存放了該點(diǎn)所經(jīng)過(guò)的窗口的梯度平均值的累加和，矩陣Matrix1為每一點(diǎn)的梯度平均值的累加次數(shù)，計(jì)算Matrix中每一點(diǎn)的平均值為矩陣Matrix2；

步驟S45：判斷如果M中某一點(diǎn)的梯度值小于Matrix2中對(duì)應(yīng)點(diǎn)梯度值的b倍，b為倍數(shù)因子，則該點(diǎn)為梯度值為0，否則該點(diǎn)的梯度值保持原值不變；

步驟S46：對(duì)梯度幅值進(jìn)行非極大值抑制；

步驟S47：用雙閾值算法檢測(cè)和連接邊緣，得到最終的邊緣圖像。

2.根據(jù)權(quán)利要求1所述的一種油藏地質(zhì)剖面圖矢量化方法，其特征在于步驟S5包括以下步驟：

步驟S51：對(duì)邊緣檢測(cè)的結(jié)果進(jìn)行連通區(qū)域標(biāo)記；

步驟S52：對(duì)標(biāo)記的每個(gè)區(qū)域采用Douglas–Peucker矢量數(shù)據(jù)壓縮算法提取特征點(diǎn)；

步驟S53：由三次貝塞爾曲線(xiàn)對(duì)提取的特征點(diǎn)進(jìn)行擬合，完成油藏地質(zhì)剖面圖的矢量化。