面向碳酸鹽巖復(fù)雜儲(chǔ)層數(shù)字巖心的各組分提取方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及測井技術(shù)數(shù)據(jù)處理領(lǐng)域，具體涉及一種面向碳酸鹽巖復(fù)雜儲(chǔ)層數(shù)字巖心的各組分提取方法。
【背景技術(shù)】
[0002]電成像測井資料縱向分辨率高，能夠清晰地反映井眼所穿過地層的各種儲(chǔ)層特征，通過研究基于成像測井資料的儲(chǔ)層各組分提取，結(jié)合巖心鑄體薄片資料，進(jìn)行三維數(shù)字巖心構(gòu)建，使數(shù)字巖心模型最大程度地與實(shí)際的巖石相同，為精確地預(yù)測儲(chǔ)層的巖石物理特性與分析不同的巖石物理響應(yīng)之間的內(nèi)在關(guān)系提供研究基礎(chǔ)，碳酸鹽巖儲(chǔ)層各組分(孔隙、連通部分和骨架)提取是三維數(shù)字巖心構(gòu)建的基礎(chǔ)，對(duì)三維數(shù)字巖心的建模精度影響很大。目前，主要利用電成像測井資料，采用圖像分割方法提取儲(chǔ)層各組分。其得到的結(jié)果精度較低，數(shù)字巖心模型與實(shí)際的巖石具有較大差異，因此，急待研究一種精度高，且得到的數(shù)字巖心模型最大程度地與實(shí)際的巖石相同的碳酸鹽巖儲(chǔ)層各組分提取方法。

【發(fā)明內(nèi)容】

[0003]有鑒于此，有必要提供一種精度高，且得到的數(shù)字巖心模型最大程度地與實(shí)際的巖石相同的碳酸鹽巖儲(chǔ)層各組分提取方法。
[0004]—種面向碳酸鹽巖復(fù)雜儲(chǔ)層數(shù)字巖心的各組分提取方法，其包括以下步驟:
[0005]S1、加載電成像測井資料，將電成像測井資料進(jìn)行預(yù)處理，得到電成像測井靜態(tài)圖像；
[0006]S2、采用圖像修復(fù)方法生成全井眼電成像圖像；
[0007]S3、將全井眼電成像圖像轉(zhuǎn)換為灰度圖像，在灰度圖像上采用雙類間最大方差法進(jìn)行碳酸鹽巖儲(chǔ)層的組分自動(dòng)提??；
[0008]S4、基于Visual c++6.0開發(fā)平臺(tái)，設(shè)置閾值調(diào)節(jié)滑動(dòng)條，采用人機(jī)交互方式手動(dòng)局部調(diào)節(jié)分割雙閾值，對(duì)碳酸鹽巖三類儲(chǔ)層各組分交互精細(xì)提取。
[0009]本發(fā)明所述面向碳酸鹽巖復(fù)雜儲(chǔ)層數(shù)字巖心的各組分提取方法，其利用電成像測井資料，采用雙類間最大方差法進(jìn)行碳酸鹽巖儲(chǔ)層的組分自動(dòng)提取，實(shí)現(xiàn)準(zhǔn)確對(duì)儲(chǔ)層的孔隙部分、連通部分和骨架部分進(jìn)行提取。該方法具有精度高，可靠性強(qiáng)的特點(diǎn)。解決了碳酸鹽巖三維數(shù)字巖心構(gòu)建，使數(shù)字巖心模型最大程度地與實(shí)際的巖石相同。
【附圖說明】
[0010]圖1為本發(fā)明所述面向碳酸鹽巖復(fù)雜儲(chǔ)層數(shù)字巖心的各組分提取方法的流程框圖；
[0011]圖2為圖1中步驟S3的子流程框圖；
[0012]圖3為圖1中步驟S3的另一子流程框圖；
[0013]圖4為全井眼電成像圖像與碳酸鹽巖三類儲(chǔ)層各組分自動(dòng)提取結(jié)果圖；其中，左邊為原始電成像圖像，右邊為雙類間最大方差法碳酸鹽巖三類儲(chǔ)層組分自動(dòng)提取結(jié)果；
[0014]圖5為碳酸鹽巖三類儲(chǔ)層各組分交互精細(xì)提取結(jié)果圖；其中，左邊為雙閾值交互界面，右邊為碳酸鹽巖三類儲(chǔ)層組分交互精細(xì)提取結(jié)果。
【具體實(shí)施方式】
[0015]為了使本發(fā)明的目的、技術(shù)方案及優(yōu)點(diǎn)更加清楚明白，以下結(jié)合附圖及實(shí)施例，對(duì)本發(fā)明進(jìn)行進(jìn)一步詳細(xì)說明，應(yīng)當(dāng)理解，此處所描述的具體實(shí)施例僅僅用以解釋本發(fā)明，并不用于限定本發(fā)明。
[0016]本發(fā)明提供一種面向碳酸鹽巖復(fù)雜儲(chǔ)層數(shù)字巖心的各組分提取方法，如圖1所示，其包括以下步驟:
[0017]S1、加載電成像測井資料，將電成像測井資料進(jìn)行預(yù)處理，得到電成像測井靜態(tài)圖像，具體的，所述預(yù)處理包括極板對(duì)齊、均衡化、產(chǎn)生動(dòng)態(tài)與靜態(tài)圖像；
[0018]S2、采用圖像修復(fù)方法生成全井眼電成像圖像；
[0019]S3、將全井眼電成像圖像轉(zhuǎn)換為灰度圖像，在灰度圖像上采用雙類間最大方差法進(jìn)行碳酸鹽巖儲(chǔ)層的組分自動(dòng)提?。?br>[0020]S4、基于Visual c++6.0開發(fā)平臺(tái)，設(shè)置閾值調(diào)節(jié)滑動(dòng)條，采用人機(jī)交互方式手動(dòng)局部調(diào)節(jié)分割雙閾值，對(duì)碳酸鹽巖三類儲(chǔ)層各組分交互精細(xì)提取。
[0021]其中，如圖2所示，所述步驟S3包括以下子步驟:
[0022]S31、使用類間最大方差法，根據(jù)前景和背景的之間方差最大得到碳酸鹽巖復(fù)雜儲(chǔ)層數(shù)字巖心分類的最佳閾值，通過最佳閾值將圖像分割成前景圖像和背景圖像兩類；
[0023]S32、在分割后的兩類中，再次使用類間最大方差法分別計(jì)算前景圖像和背景圖像的子類中的最佳分類，并求得兩個(gè)最佳分類的類間最大方差；
[0024]S33、分別選取兩個(gè)類間最大方差的對(duì)應(yīng)閾值作為碳酸鹽巖三類儲(chǔ)層的分割閾值。
[0025]所述碳酸鹽巖三類儲(chǔ)層包括孔隙部分、連通部分和骨架部分，對(duì)碳酸鹽巖儲(chǔ)層各組分等三類分類問題，需要兩個(gè)閾值進(jìn)行分類，因此，提出雙類間最大方差法，所述雙類間最大方差法為對(duì)灰度圖像進(jìn)行兩次最大類間方差法分隔，屬于基于類間最大方差法的一種面向三類問題分割算法。
[0026]具體的，在雙類間最大方差法中，設(shè)t為前景圖像與背景圖像的分割閾值，前景圖像點(diǎn)數(shù)占圖像比例為wo，平均灰度為UQ，背景圖像點(diǎn)數(shù)占圖像比例為W1，平均灰度為m;
[0027]則所述圖像的總平均灰度值:
[0028]u=wo*uo+wi*ui (I)
[0029]所述前景圖象和背景圖象的方差:
[0030]g=wo* (Uo-U )*( Uo-U)+wi (u1-u )*( u1-u) (2)
[0031 ] =wo*wi*(uo-ui)*(uo-ui)
[0032]其中，當(dāng)前景和背景類間方差最大時(shí)，前景圖像和背景圖像差異最大，所對(duì)應(yīng)的灰度t是兩類分類的最佳閾值。
[0033]如圖3所示，所述采用雙類間最大方差法進(jìn)行碳酸鹽巖儲(chǔ)層的組分自動(dòng)提取包括以下分步驟:
[0034]S31、根據(jù)最佳閾值t將碳酸鹽巖電成像靜態(tài)灰度圖像P分割成pi和P2兩部分；
[0035]在分割后的兩部分中，對(duì)P1部分采用類間最大方差法確定閾值^，對(duì)?2部分采用類間最大方差法確定閾值t3，并根據(jù)閾值t2將Pl分割為Pll和P12，根據(jù)閾值t3將P2分割為P21和
P22；
[0036]從t2和t3兩個(gè)閾值中選取兩者中類間方差最大的對(duì)應(yīng)閾值作為最佳閾值，由和