面向碳酸鹽巖復雜儲層數字巖心的各組分提取方法
【技術領域】
[0001]本發明涉及測井技術數據處理領域，具體涉及一種面向碳酸鹽巖復雜儲層數字巖心的各組分提取方法。
【背景技術】
[0002]電成像測井資料縱向分辨率高，能夠清晰地反映井眼所穿過地層的各種儲層特征，通過研究基于成像測井資料的儲層各組分提取，結合巖心鑄體薄片資料，進行三維數字巖心構建，使數字巖心模型最大程度地與實際的巖石相同，為精確地預測儲層的巖石物理特性與分析不同的巖石物理響應之間的內在關系提供研究基礎，碳酸鹽巖儲層各組分(孔隙、連通部分和骨架)提取是三維數字巖心構建的基礎，對三維數字巖心的建模精度影響很大。目前，主要利用電成像測井資料，采用圖像分割方法提取儲層各組分。其得到的結果精度較低，數字巖心模型與實際的巖石具有較大差異，因此，急待研究一種精度高，且得到的數字巖心模型最大程度地與實際的巖石相同的碳酸鹽巖儲層各組分提取方法。
【發明內容】
[0003]有鑒于此，有必要提供一種精度高，且得到的數字巖心模型最大程度地與實際的巖石相同的碳酸鹽巖儲層各組分提取方法。
[0004]—種面向碳酸鹽巖復雜儲層數字巖心的各組分提取方法，其包括以下步驟:
[0005]S1、加載電成像測井資料，將電成像測井資料進行預處理，得到電成像測井靜態圖像；
[0006]S2、采用圖像修復方法生成全井眼電成像圖像；
[0007]S3、將全井眼電成像圖像轉換為灰度圖像，在灰度圖像上采用雙類間最大方差法進行碳酸鹽巖儲層的組分自動提取；
[0008]S4、基于Visual c++6.0開發平臺，設置閾值調節滑動條，采用人機交互方式手動局部調節分割雙閾值，對碳酸鹽巖三類儲層各組分交互精細提取。
[0009]本發明所述面向碳酸鹽巖復雜儲層數字巖心的各組分提取方法，其利用電成像測井資料，采用雙類間最大方差法進行碳酸鹽巖儲層的組分自動提取，實現準確對儲層的孔隙部分、連通部分和骨架部分進行提取。該方法具有精度高，可靠性強的特點。解決了碳酸鹽巖三維數字巖心構建，使數字巖心模型最大程度地與實際的巖石相同。
【附圖說明】
[0010]圖1為本發明所述面向碳酸鹽巖復雜儲層數字巖心的各組分提取方法的流程框圖；
[0011]圖2為圖1中步驟S3的子流程框圖；
[0012]圖3為圖1中步驟S3的另一子流程框圖；
[0013]圖4為全井眼電成像圖像與碳酸鹽巖三類儲層各組分自動提取結果圖；其中，左邊為原始電成像圖像，右邊為雙類間最大方差法碳酸鹽巖三類儲層組分自動提取結果；
[0014]圖5為碳酸鹽巖三類儲層各組分交互精細提取結果圖；其中，左邊為雙閾值交互界面，右邊為碳酸鹽巖三類儲層組分交互精細提取結果。
【具體實施方式】
[0015]為了使本發明的目的、技術方案及優點更加清楚明白，以下結合附圖及實施例，對本發明進行進一步詳細說明，應當理解，此處所描述的具體實施例僅僅用以解釋本發明，并不用于限定本發明。
[0016]本發明提供一種面向碳酸鹽巖復雜儲層數字巖心的各組分提取方法，如圖1所示，其包括以下步驟:
[0017]S1、加載電成像測井資料，將電成像測井資料進行預處理，得到電成像測井靜態圖像，具體的，所述預處理包括極板對齊、均衡化、產生動態與靜態圖像；
[0018]S2、采用圖像修復方法生成全井眼電成像圖像；
[0019]S3、將全井眼電成像圖像轉換為灰度圖像，在灰度圖像上采用雙類間最大方差法進行碳酸鹽巖儲層的組分自動提取；
[0020]S4、基于Visual c++6.0開發平臺，設置閾值調節滑動條，采用人機交互方式手動局部調節分割雙閾值，對碳酸鹽巖三類儲層各組分交互精細提取。
[0021]其中，如圖2所示，所述步驟S3包括以下子步驟:
[0022]S31、使用類間最大方差法，根據前景和背景的之間方差最大得到碳酸鹽巖復雜儲層數字巖心分類的最佳閾值，通過最佳閾值將圖像分割成前景圖像和背景圖像兩類；
[0023]S32、在分割后的兩類中，再次使用類間最大方差法分別計算前景圖像和背景圖像的子類中的最佳分類，并求得兩個最佳分類的類間最大方差；
[0024]S33、分別選取兩個類間最大方差的對應閾值作為碳酸鹽巖三類儲層的分割閾值。
[0025]所述碳酸鹽巖三類儲層包括孔隙部分、連通部分和骨架部分，對碳酸鹽巖儲層各組分等三類分類問題，需要兩個閾值進行分類，因此，提出雙類間最大方差法，所述雙類間最大方差法為對灰度圖像進行兩次最大類間方差法分隔，屬于基于類間最大方差法的一種面向三類問題分割算法。
[0026]具體的，在雙類間最大方差法中，設t為前景圖像與背景圖像的分割閾值，前景圖像點數占圖像比例為wo，平均灰度為UQ，背景圖像點數占圖像比例為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]其中，當前景和背景類間方差最大時，前景圖像和背景圖像差異最大，所對應的灰度t是兩類分類的最佳閾值。
[0033]如圖3所示，所述采用雙類間最大方差法進行碳酸鹽巖儲層的組分自動提取包括以下分步驟:
[0034]S31、根據最佳閾值t將碳酸鹽巖電成像靜態灰度圖像P分割成pi和P2兩部分；
[0035]在分割后的兩部分中，對P1部分采用類間最大方差法確定閾值^，對?2部分采用類間最大方差法確定閾值t3，并根據閾值t2將Pl分割為Pll和P12，根據閾值t3將P2分割為P21和
P22；
[0036]從t2和t3兩個閾值中選取兩者中類間方差最大的對應閾值作為最佳閾值，由和