李戈?yáng)|，馬猛，王建立，陳銘帥，杜波，韓明剛

(1.中海油能源發(fā)展股份有限公司 工程技術(shù)分公司，天津 300452；2.中海石油(中國(guó))有限公司 天津分公司，天津 300459)

FLAIR實(shí)時(shí)地層流體分析錄井(簡(jiǎn)稱FLAIR錄井)能實(shí)時(shí)檢測(cè)鉆井液中C1～C8烴類組分以及CO2和H2S等非烴類組分含量[1]。本文針對(duì)研究區(qū)的基本地質(zhì)情況和FLAIR實(shí)時(shí)地層流體錄井特點(diǎn)，利用統(tǒng)計(jì)分析的ReliefF算法對(duì)氣測(cè)敏感參數(shù)進(jìn)行計(jì)算優(yōu)選[2]，選出對(duì)含油性敏感度最高的兩個(gè)參數(shù)，根據(jù)優(yōu)選后的參數(shù)建立圖板，并編寫算法，由軟件自動(dòng)劃分油層與非油層區(qū)域，免除了人為劃分區(qū)域的主觀性，并對(duì)已建立的圖板解釋效果進(jìn)行驗(yàn)證。

1 解釋參數(shù)優(yōu)選

與常規(guī)氣測(cè)錄井技術(shù)相比，F(xiàn)LAIR錄井具有脫氣效率高、分辨率高(C1/C2比高達(dá)8000)、測(cè)量范圍更廣(烴類可測(cè)C1～C8)、數(shù)據(jù)質(zhì)量更高的特點(diǎn)[3]。

1.1 相關(guān)參數(shù)介紹

氣測(cè)原始數(shù)據(jù)所受影響因素較多，因此在建立氣測(cè)交會(huì)圖版過程中多采用其派生參數(shù)，共篩選出以下25個(gè)氣測(cè)相關(guān)派生參數(shù)[4-9]:

(1)CR1:100*C1/∑C；(2)CR6+:100*(nC6+nC7+nC8+C7H14)/∑C；(3)(C4+C5)/(C1+C2)；(4)(C1+C2)/C3；(5)(C1+C2)/(C3+C4) ；(6)(C2+C3+C4)/C1；(7)(C3+C4)/C4；(8)C1/C3；(9)RL-M：10*C1/(C2+C3)2；(10)RL-H：100*(C1+C2)/(iC4+nC4+iC5+nC5)3；(11)RH-M：(iC4+nC4+iC5+nC5)2/C3；(12)氣測(cè)幅度指數(shù)：Tg/Tj/Ht；(13)重?zé)N指數(shù)：C1*Tj/Tg* (iC4+nC4+iC5+nC5)/(0.0001+C3)；(14)豐度系數(shù)：C1/(C1-C2-C3)*∑C*(Tg-Tj)/Tg；(15)氣測(cè)含油豐度指數(shù)：[C1/(C1-C2)]2*∑C*(Tg-Tj)/Tg；(16)油性指數(shù)：Hq/Ht*(Tg-Tj)/Tj*[(nC4+nC5)/(iC4+iC5)]2*(iC4+nC4+iC5+nC5)2/C3；(17)(C1-C2)/C1；(18)C1/(C1-C2-C3)；(19)nC4/iC4；(20)nC7+nC8+C6H6+C7H8+C7H14；(21)iC5+nC5+nC6；(22)(nC7+nC8+C6H6+C7H8+C7H14)/(iC5+nC5+nC6)；(23)iC5+nC5+nC6+nC7+nC8；(24)C6H6+C7H8+C7H14；(25)nC8+C6H6+C7H8+C7H14。式中：C1、C2、C3、C4、iC4、nC4、iC5、nC5、nC6、C6H6、nC7、C7H8、C7H14、nC8為FLAIR錄井在油氣層段的氣體組分檢測(cè)值，%；∑C為FLAIR錄井氣體各組分含量平均值的加和，%；Tg為氣測(cè)異常顯示段全烴平均值，%；Tj為氣測(cè)異常顯示段上部非顯示段氣測(cè)基值，%；Ht為氣測(cè)異常顯示段對(duì)應(yīng)的低鉆時(shí)厚度，以半幅點(diǎn)為界，m；Hq為低鉆時(shí)段對(duì)應(yīng)的氣測(cè)異常顯示厚度，以半幅點(diǎn)為界，m。

1.2 解釋參數(shù)優(yōu)選方法

由于研究區(qū)域油藏埋深跨度較大，且各層位的油質(zhì)不同，導(dǎo)致淺層和中深層的氣測(cè)錄井特征具有明顯差異，因此需要將樣本數(shù)據(jù)分為淺層樣本數(shù)據(jù)和中深層樣本數(shù)據(jù)，并針對(duì)淺層和中深層分別優(yōu)選敏感參數(shù)[10]。對(duì)于淺層與中深層的劃分，這里根據(jù)樣本層中油層的氣測(cè)相對(duì)組分含量隨深度的變化情況確定其界限，相對(duì)組分含量采用公式Cn/∑C計(jì)算獲得，其中Cn為氣測(cè)異常段各組分含量平均值，∑C為各組分含量平均值的加和[11]。

圖1反映了研究區(qū)域不同深度油層氣測(cè)相對(duì)組分含量的變化情況。萊州灣凹陷在1600 m之上，C1的比例占到了90%以上，C2及其以后各組分占了極小一部分比例，甚至在1600 m之下，C1所占比例逐漸降低，其他各組分比例均逐漸升高，組分齊全。根據(jù)這一特征，將1600 m作為萊州灣凹陷淺層與中深層的分界，分類情況見表1。

表1 研究區(qū)域樣本分類情況研究區(qū)域縱向分界樣本層數(shù)萊州灣凹陷區(qū)域(墾利)1600 m以上1600 m以下88112圖1 萊州灣凹陷區(qū)域FLAIR在不同深度下氣測(cè)比例變化

由于參數(shù)較多，隨機(jī)抽選參數(shù)建立圖板比較復(fù)雜，因此使用ReliefF算法對(duì)研究區(qū)域的各派生參數(shù)進(jìn)行特征權(quán)重計(jì)算。Relief算法最早由Kira提出，最初局限于兩類數(shù)據(jù)的分類問題。Relief算法是一種特征權(quán)重算法(Feature Weighting Algorithms)，根據(jù)各個(gè)特征和類別的相關(guān)性賦予特征不同的權(quán)重，權(quán)重小于某個(gè)閾值的特征將被移除。Relief算法比較簡(jiǎn)單，但運(yùn)行效率高，并且結(jié)果也比較令人滿意，因此得到廣泛應(yīng)用，但是其局限性在于只能處理兩類別數(shù)據(jù)，因此1994年Kononeill對(duì)其進(jìn)行了擴(kuò)展，得到了ReliefF算法，可以處理多類別問題。ReliefF算法在處理多類問題時(shí)，每次從訓(xùn)練樣本集中隨機(jī)取出一個(gè)樣本R，然后從和R同類的樣本集中找出R的k個(gè)近鄰樣本(near Hits)，從每個(gè)R的不同類的樣本集中均找出k個(gè)近鄰樣本(near Misses)，然后更新每個(gè)特征的權(quán)重[12]。

本文采用ReliefF算法來計(jì)算1.1節(jié)中提到的25個(gè)參數(shù)的各個(gè)特征的權(quán)重，選取權(quán)重最大的兩個(gè)值建立解釋圖版。由于算法在運(yùn)行過程中，會(huì)選擇隨機(jī)樣本R，隨機(jī)數(shù)的不同將導(dǎo)致結(jié)果權(quán)重有一定的出入，因此本文采取平均的方法，將主程序運(yùn)行20次，避免單次運(yùn)行出現(xiàn)誤差。如圖2、圖3所示，縱坐標(biāo)為計(jì)算得到的特征權(quán)重，橫坐標(biāo)為各氣測(cè)相關(guān)派生參數(shù)(對(duì)應(yīng)1.1中提到的(1)～(25))。

圖2 萊州灣凹陷區(qū)域1600 m以上FLAIR氣測(cè)派生參數(shù)特征權(quán)重 圖3 萊州灣凹陷區(qū)域1600 m以下FLAIR氣測(cè)派生參數(shù)特征權(quán)重

針對(duì)不同區(qū)域分別對(duì)以上的25個(gè)參數(shù)進(jìn)行特征權(quán)重計(jì)算，各選出兩個(gè)特征權(quán)重最高的參數(shù)用以建立圖版。因此，萊州灣凹陷1600 m以上區(qū)域采用(C3+C4)/C4與(C1-C2)/C1兩個(gè)派生參數(shù)，萊州灣凹陷1600 m以下區(qū)域采用CR6+與RH-M兩個(gè)派生參數(shù)。

2 解釋圖版建立

通過以上分析，可以得出研究區(qū)對(duì)儲(chǔ)集層含油性最敏感的兩個(gè)參數(shù)，然后使用這兩個(gè)參數(shù)分別建立圖板(見圖4、圖5)。將優(yōu)選出的氣測(cè)派生參數(shù)導(dǎo)入軟件后，再導(dǎo)入編寫好的分區(qū)公式，軟件生成的交會(huì)圖版將直接自動(dòng)畫線區(qū)分油層與非油層區(qū)域，避免了人為劃區(qū)域的主觀性。

圖4 萊州灣凹陷區(qū)域FLAIR1600 m以上(C3+C4)/C4—(C1-C2)/C1圖版 圖5 萊州灣凹陷區(qū)域FLAIR1600 m以下CR6+—RH-M圖版

經(jīng)過計(jì)算后建立圖版，可得研究區(qū)不同深度圖版識(shí)別率，見表2。

表2 FLAIR錄井圖版識(shí)別率統(tǒng)計(jì)

3 應(yīng)用效果分析

針對(duì)兩個(gè)研究區(qū)，選取了5口新鉆井的10個(gè)已確定流體類型的儲(chǔ)集層用于FLAIR錄井解釋圖版應(yīng)用驗(yàn)證，識(shí)別結(jié)果見表3。

表3 FLAIR圖版應(yīng)用驗(yàn)證統(tǒng)計(jì)

與測(cè)試結(jié)論對(duì)比，F(xiàn)LAIR解釋圖版對(duì)油層與非油層的判別準(zhǔn)確率可以達(dá)到90.0%，證實(shí)了所建圖板的有效性。不符合層是因其均位于油水界面處，緊鄰其上方即為油層，氣測(cè)數(shù)據(jù)受到了油層的影響，所以表現(xiàn)出油層的特征。

4 結(jié)論

(1)利用ReliefF算法對(duì)氣測(cè)派生參數(shù)進(jìn)行篩選，根據(jù)不同區(qū)域不同深度，從25個(gè)參數(shù)中分別優(yōu)選出兩個(gè)參數(shù)(萊州灣凹陷1600 m以上區(qū)域：(C3+C4)/C4與(C1-C2)/C1；萊州灣凹陷1600 m以下區(qū)域：CR6+與RH-M)并導(dǎo)入軟件、利用公式進(jìn)行了自動(dòng)分區(qū)，建立了油層與非油層判別圖版。

(2)在兩個(gè)研究區(qū)選取了5口新鉆井的10個(gè)已確定流體類型的儲(chǔ)集層用于FLAIR圖版應(yīng)用驗(yàn)證，準(zhǔn)確率達(dá)到90.0%，圖版應(yīng)用效果良好。該解釋圖版建立方法簡(jiǎn)便、快捷，具有良好的推廣應(yīng)用價(jià)值。