基于主成分分析的測井曲線識別不整合面方法

姜虹

摘 要：本文介紹了利用測井曲線識別不整合面的方法。利用歸一化和主成分分析方法，對測井曲線進行主成分提取，將主成分變量進行最優分割計算，構建一條綜合分層曲線，通過分析綜合分層曲線數值與不整合面的對應關系，實現了對不整合面的識別。

關鍵詞：歸一化；主成分分析；綜合分層曲線；不整合面

0 引言

利用測井曲線識別不整合面需要優選多條測井曲線，選取的原則：能夠綜合反映地層不整合的信息；且在不整合面附近的響應特征出現明顯異常的變化。此方法優點是信息量大且全面，但在對其進行運算過程中，變量多，反映的地層信息有所重復。因此，通過對多條測井曲線進行主成分分析，利用降維處理，消除或減少重復信息，將多條測井曲線復合成少數幾個綜合變量，降維的同時基本不損失原始曲線所反映的地層信息。在此基礎上，對利用主成分分析提取出的少數幾個變量進行最優分割方法的計算，構建的綜合分層曲線可以對不整合面進行識別。

1 測井曲線識別不整合面方法

基于測井曲線及其對應的地質意義，在利用測井曲線識別不整合面時，需要選用以下常規曲線綜合分析來進行劃分：AC、DEN、CNL、COND、RLML、RNML、GR、SP。由于各測井參數的量綱不同，不能直接將它們放在一起計算。因此，為了消除數量級差異給測井參數融合帶來的影響，需對測井參數進行歸一化處理。

假設測井信號n個深度點中第i個深度點的第j條測井曲線數值為Xij，對選取的測井數據進行歸一化處理后，得出第i個深度點的第j條測井曲線數據歸一化結果為Xij[1]。

對以上8組數據進行歸一化處理之后，以主成分累積方差百分比大于80%為準，進行主成分提取，見表1，選用變量F1，F2，F3，這三個變量就能表征原來8個變量85.529%的信息，降維的同時，消除或減少重復的信息，提取了有效信息。

提取的主成分F1，F2，F3的表達式如式（1-3，4，5），式中的變量Xi*都為標準化之后的變量，標準化規則為，其中μi為樣本均值，Sij為樣本方差。

由自變量解釋能力數值可以看出，主成分F1主要反映了地層的孔滲以及巖性特征，主成分F2反映了地層的巖性及流體性質特征，主成分F3則反映了一些補充信息。因此，選取F1，F2，F3這三個主成分變量進一步構建綜合分層曲線。

將主成分分析得到的F1，F2，F3這三個主成分變量，進行最優分割計算，構建得到一條綜合分層曲線。曲線的極小值點，為最優分割點，最小值位置對應的深度，其地質意義為地層變化最劇烈的位置，即為不整合面深度。

2 測井曲線識別不整合面的應用

圖1中的綜合分層曲線是將某A井的8條測井曲線進行主成分提取，得到F1，F2，F3這三個主成分變量，運用最優分割計算對這三個主成分進行處理得到一條綜合分層曲線，可以準確識別出Es2/Es3不整合界面的位置在曲線的極小值點1197m處。結合地震及巖心分析等資料，劃分的不整合界面位置準確。

3 結論

基于測井曲線及其對應的地質意義，優選8條常規測井曲線，利用歸一化和主成分分析方法對其進行處理后，有效提取了三個主成分，降維的同時提取了有效信息。將得到的主成分進行最優分割計算，構建了一條綜合分層曲線，其極小值對應深度為不整合面位置，從而實現了對地層不整合面的識別。

參考文獻：

[1]張蕾.濟陽凹陷下第三系不整合結構測井評價[D].山東：中國石油大學（華東），2008：1-78.