[發(fā)明專利]面向碳酸鹽巖復(fù)雜儲(chǔ)層數(shù)字巖心的各組分提取方法

在審

• 申請(qǐng)?zhí)枺?/td>: 201510875673.6
• 申請(qǐng)日: 2015-12-03
• 公開（公告）號(hào): CN105551018A
• 公開（公告）日: 2016-05-04
• 發(fā)明（設(shè)計(jì)）人: 張翔;肖小玲;屈樂(lè);杜環(huán)虹;余春昊
• 申請(qǐng)（專利權(quán)）人: 長(zhǎng)江大學(xué)
• 主分類號(hào): G06T7/00
• 分類號(hào): G06T7/00
• 代理公司: 武漢河山金堂專利事務(wù)所(普通合伙) 42212
• 代理人: 胡清堂

1.一種面向碳酸鹽巖復(fù)雜儲(chǔ)層數(shù)字巖心的各組分提取方法，其特征在于，所述面向碳 酸鹽巖復(fù)雜儲(chǔ)層數(shù)字巖心的各組分提取方法包括以下步驟：

S1、加載電成像測(cè)井資料，將電成像測(cè)井資料進(jìn)行預(yù)處理，得到電成像測(cè)井靜態(tài)圖像；

S2、采用圖像修復(fù)方法生成全井眼電成像圖像；

S3、將全井眼電成像圖像轉(zhuǎn)換為灰度圖像，在灰度圖像上采用雙類間最大方差法得到 分組閾值，通過(guò)分組閾值將碳酸鹽巖儲(chǔ)層的組分進(jìn)行自動(dòng)提取；

S4、基于Visualc++6.0開發(fā)平臺(tái)，設(shè)置閾值調(diào)節(jié)滑動(dòng)條，采用人機(jī)交互方式手動(dòng)局部 調(diào)節(jié)分割雙閾值，對(duì)碳酸鹽巖三類儲(chǔ)層各組分交互精細(xì)提取。

2.根據(jù)權(quán)利要求1所述的面向碳酸鹽巖復(fù)雜儲(chǔ)層數(shù)字巖心的各組分提取方法，其特征 在于，所述步驟S3包括以下子步驟：

S31、使用類間最大方差法，根據(jù)前景圖象和背景圖象的方差最大時(shí)得到碳酸鹽巖復(fù)雜 儲(chǔ)層數(shù)字巖心分類的最佳閾值，通過(guò)最佳閾值將圖像分割成前景圖像和背景圖像兩類；

S32、在分割后的兩類中，再次使用最大類間方差法分別計(jì)算前景圖像和背景圖像的子 類中的最佳分類，并求得兩個(gè)最佳分類的類間最大方差；

S33、分別選取兩個(gè)類間最大方差最大的方差對(duì)應(yīng)閾值作為碳酸鹽巖三類儲(chǔ)層的分割 閾值。

3.根據(jù)權(quán)利要求2所述的面向碳酸鹽巖復(fù)雜儲(chǔ)層數(shù)字巖心的各組分提取方法，其特征 在于，

所述圖像的總平均灰度值：

u＝w₀*u₀+w₁*u₁(1)

所述前景圖象和背景圖象的方差：

g＝w₀*(u₀-u)*(u₀-u)+w₁(u₁-u)*(u₁-u)(2)

＝w₀*w₁*(u₀-u₁)*(u₀-u₁)

其中：t為前景圖像與背景圖像的分割閾值，前景圖像點(diǎn)數(shù)占圖像比例為w₀，平均灰度為 u₀，背景圖像點(diǎn)數(shù)占圖像比例為w₁，平均灰度為u₁。

4.根據(jù)權(quán)利要求3所述的面向碳酸鹽巖復(fù)雜儲(chǔ)層數(shù)字巖心的各組分提取方法，其特征 在于，

當(dāng)前景和背景類間方差最大時(shí)，所對(duì)應(yīng)的灰度t是兩類分類的最佳閾值。

5.根據(jù)權(quán)利要求4所述的面向碳酸鹽巖復(fù)雜儲(chǔ)層數(shù)字巖心的各組分提取方法，其特征 在于，所述采用雙類間最大方差法進(jìn)行碳酸鹽巖儲(chǔ)層的組分自動(dòng)提取包括以下分步驟：

S31、根據(jù)最佳閾值t將碳酸鹽巖電成像靜態(tài)灰度圖像p分割成p₁和p₂兩部分；

S32、在分割后的兩部分中，對(duì)p₁部分采用類間最大方差法確定閾值t₂，對(duì)p₂部分采用類 間最大方差法確定閾值t₃，并根據(jù)閾值t₂將p₁分割為p₁₁和p₁₂，根據(jù)閾值t₃將p₂分割為p₂₁和 p₂₂；

S33、從t₂和t₃兩個(gè)閾值中選取兩者中類間方差最大的對(duì)應(yīng)閾值作為最佳閾值t₁,由t₁和 t將碳酸鹽巖電成像靜態(tài)灰度圖像分割成三類，即碳酸鹽巖儲(chǔ)層組分。