魏玉梅，劉文鈺，劉俊東，劉永河，竇如勝，王順利，孟姝婳

(1.中國石油集團測井有限公司 天津分公司，天津 300280；2.中國石油天然氣股份有限公司大港油田分公司 第四采油廠，天津 300280；3.中國石油集團渤海鉆探工程有限公司 第一錄井分公司，天津 300280)

1 引 言

儲層巖性、物性、電性及含油性可以用測井資料進行表征，而電性對含油性最為敏感[1,2]，同樣，電性又受巖性、物性的影響較大，也有可能會受地層水礦化度[3]及泥巖電阻率的影響。由于電阻率的影響因素較多[4]，在進行儲層流體性質評價中會出現(xiàn)多解性[5]。針對這一難題，通過對測井資料在油水層的表征進行分析，厘清油水層敏感曲線[6]，利用數(shù)理統(tǒng)計技術，將能夠綜合反映儲層流體性質特征的多種測井資料結合，建立儲層分類評價的多參數(shù)圖版[7]和定量判別函數(shù)，對未知類型的儲層進行分類。該方法能夠降低異常數(shù)據對分類結果的不良影響[8]，綜合考慮多項測井參數(shù)，使分類結果更加可靠真實。

2 數(shù)理統(tǒng)計技術在儲層流體性質識別中的應用

2.1 數(shù)理統(tǒng)計方法原理

本次方法研究，主要是利用數(shù)理統(tǒng)計方法中的均值聚類和Fisher判別分析法[9]，通過對一些樣本的諸多已知性質進行分析提取、建立樣本分類規(guī)則來實現(xiàn)對未分類樣本的歸類。它的幾何意義是通過選擇適當?shù)耐队胺较?，把由幾個變量構成的高維空間的點群向低維空間投影，并進行適當?shù)淖儞Q，把幾個原始變量組合成一個新的變量R，使得不同類樣本的R值盡可能分離開來，達到將樣本歸類的目的[10,11]。

根據數(shù)理統(tǒng)計評價原理，將能夠反映儲層流體性質的測井參數(shù)歸類，達到對油層、油水同層、水層的定性和定量評價。

2.2 儲層流體性質識別方法研究

2.2.1 選取儲層流體識別參數(shù)

選取大港油田板中東斷塊2019～2020年已試采井的測井資料，遵循參數(shù)間應盡可能獨立的原則，從測井響應參數(shù)中選取最能反映油水層特點的參數(shù)[12]：反映儲層含油性的深電阻率(Rt)；反映儲層鄰近純泥巖的電阻率(Rsh)；反映儲層巖性的自然伽馬(GR)；反映儲層鄰近純砂巖和純泥巖的自然伽馬(GRmn、GRmx)；反映儲層物性特征參數(shù)[13]的聲波(AC)，將以上參數(shù)作為判別分析的觀察值[14]。選取的特征參數(shù)，基本能消除由于巖性、物性等因素對電阻率的影響，更準確的反映真實儲層流體性質。

選取有代表性的試采井14口井93個數(shù)據點(油層64個數(shù)據點、油水同層16個數(shù)據點、水層13個數(shù)據點)作為已知油層、油水同層和水層3類儲層的已知樣本。

為了驗證選擇數(shù)據的可靠性，將93個數(shù)據點進行了聚類分析，從聚類譜系圖(圖1)中可以看出，儲層可以分為3個類別，正好符合研究中油層、油水同層、水層三類儲層的劃分，也證明資料對儲層流體性質的劃分是可靠的。

圖1 聚類譜系圖Fig.1 Cluster pedigree

2.2.2 數(shù)理統(tǒng)計技術在儲層流體性質識別中的應用

1)均值聚類分析建立多參數(shù)油水解釋圖版

傳統(tǒng)的雙參數(shù)——聲波與深電阻率油水解釋圖版(圖2)，由于存在高電阻率水層，點子雜亂，同一區(qū)域，油層、油水同層、水層有重合的現(xiàn)象，無法區(qū)分。

圖2 雙參數(shù)油水解釋圖版Fig.2 Oil-water interpretation chart of two-parameters

通過均值聚類分析，把每個樣品作為一類，然后把最靠近的樣品(即距離最小的樣品)首先聚為小類，再將已聚合的小類按其類間距離再合并，不斷繼續(xù)下去，最后把一切子類都聚合到一個大類[15,16]。將試采井的93個數(shù)據點進行歸類，建立了多參數(shù)油水解釋圖版，從圖中看，油層、油水同層、水層能夠明顯的區(qū)分開(圖3)[17]，與雙參數(shù)圖版對比，優(yōu)勢明顯。

圖3 多參數(shù)油水解釋圖版Fig.3 Oil-water interpretation chart of multi-parameters

同時，在圖版建立之時，根據判別分析技術，對各特征參數(shù)求取判別向量，即第 1、第 2判別向量，進而計算出綜合指標函數(shù)1和函數(shù)2，即公式(1)和公式(2)。

公式中：y1、y2為圖3中函數(shù)1、函數(shù)2，單位是無量綱；GR為自然伽馬，單位是API；GRmx為純泥巖自然伽馬，單位是API；GRmn為純砂巖自然伽馬，單位是API；Rsh為純泥巖電阻率，單位是Ω·m；Rt為地層深電阻率，單位是Ω·m；AC為聲波，單位是μs/m。

利用已知數(shù)據對圖版的準確性進行驗證，能夠正確地對94.6 %的原數(shù)據進行分類(表1)，準確率較高。

表1 分類結果Table 1 Classification results

2)Fisher判別分析法建立儲層分類函數(shù)模型

另外，根據已知樣品類型，運用判別分析技術，建立了油層、油水同層、水層3類層流體性質識別函數(shù)模型[18]，詳見公式(3)至公式(5),實現(xiàn)了油水層定量評價標準[19]。

油層:

(3)

油水同層:

(4)

水層:

式中：GR′為自然伽馬相對值，見公式(6)，無量綱。

最終計算的3個函數(shù)值，哪個值最大，該點就屬于相應的哪一類。該方法達到對流體性質進行定量識別的目的。

3 應用效果分析

應用數(shù)理統(tǒng)計技術，將能夠反映儲層特征的多種測井資料進行融合，建立了多參數(shù)油水解釋圖版和定量判別技術，對板中東斷塊的17口井進行重新評價，解釋符合率由之前的81.6 %提高到了92.1 %(表2)，取得了較好的應用效果[20]。

表2 解釋符合率統(tǒng)計Table 2 The coincidence rate of interpretation

板A井位于板中東斷塊，106號層頂部(3 550～3 554.6 m)與109號層，電阻率值較高，均在35 Ω·m左右，巖性純、物性好，兩個層的差異不明顯，利用常規(guī)解釋標準很難將兩個層區(qū)分開。將本次的研究成果應用在板A井中，兩個層的差異明顯就體現(xiàn)出來，在定性解釋圖版中(圖4)，106號層頂部落在油水同層區(qū)，109號層落在油層區(qū)；定量解釋結果(圖5)與定性結果一致，106號層頂部油水同層，109號層為油層。

圖4 板A井106、109號層圖版位置Fig.4 Plate location of layers 106 and 109 in well ban A

圖5 板A井判別分析解釋成果Fig.5 Diagram of discriminant analysis interpretation of well ban A

106(3 550.5～3 554 m)號層和109(3 593.8～3 596.2 m)號層合采，日產油1.58噸，氣3 469方，水34.2方，含水95.6 %，試采結果與解釋結果一致，驗證了研究方法適用性。

4 結 論

利用數(shù)理統(tǒng)計技術，建立了函數(shù)判別圖版和定量評價函數(shù)模型，能夠將綜合反映儲層特征的眾多信息相結合，在油水層識別方面具有很強的優(yōu)勢。該方法應用在板中東斷塊油水層識別中見到較好的效果。

1)數(shù)理統(tǒng)計方法，能夠對已知性質進行提取、歸類，適用在儲層流體性質識別中；

2)應用數(shù)理統(tǒng)計方法，選取目標區(qū)塊的試采井資料，將對油水層反映敏感的測井資料進行融合、歸類，建立的多參數(shù)流體性質識別圖版與傳統(tǒng)的雙參數(shù)圖版對比，效果顯著；

3)建立了判別分析函數(shù)模型，達到了流體性質的定量判別，在油水層識別中具有智能、快速的優(yōu)勢；

4)在目標區(qū)塊中應用效果顯著，引用多參數(shù)，實現(xiàn)了流體識別技術的創(chuàng)新。