張延旭，姜 晶，王 濤，王 凱

1.中海石油(中國)有限公司蓬勃作業(yè)公司，天津 300459；2.中海石油(中國)有限公司天津分公司，天津 300459；3. 中海油田服務(wù)股份有限公司，天津 300450；4. 中海油能源發(fā)展股份有限公司工程技術(shù)分公司，天津 300452

國內(nèi)外已有的研究和應(yīng)用成果表明，油氣藏是封閉條件良好的地下儲氣庫，可實現(xiàn)CO2的長期埋存[1]。影響油藏CO2埋存效果的因素很多，且各影響因素與埋存效果之間是一種非線性、不確定的復(fù)雜關(guān)系，致使常規(guī)的數(shù)模、實驗方法費時費力，且多維回歸方法也存在模型建立困難、計算量大、精度較低的缺點[2-3]。而神經(jīng)網(wǎng)絡(luò)法具有較強的自學(xué)習、自適應(yīng)、自調(diào)整能力，其強大的非線性逼近能力使其很好的逼近被測對象的實際值[4-5]。因此，在綜合考慮各影響因素的基礎(chǔ)上，筆者引入反映油藏埋存CO2效果的5個無因次變量，嘗試采用BP神經(jīng)網(wǎng)絡(luò)對油藏埋存CO2效果進行評價預(yù)測，發(fā)現(xiàn)人工神經(jīng)網(wǎng)絡(luò)法預(yù)測精度較高，具有一定的應(yīng)用價值。

1 神經(jīng)網(wǎng)絡(luò)預(yù)測方法的基本原理

神經(jīng)網(wǎng)絡(luò)技術(shù)實際上是“人工神經(jīng)網(wǎng)絡(luò)”或“模擬神經(jīng)網(wǎng)絡(luò)”的簡稱，是將人腦功能加以模擬化的一種計算機信息處理系統(tǒng)。到目前為止，已研制出了多種人工神經(jīng)網(wǎng)絡(luò)模型和相應(yīng)的學(xué)習算法[6]。

通過人工神經(jīng)網(wǎng)絡(luò)法求出的公式是如下“隱含”表達式：

y=ANN(x1,x2,…,xm)

式中，ANN是一個非線性函數(shù)。這個函數(shù)不能用通常的數(shù)學(xué)公式表示，稱為“知識庫”。筆者采用目前應(yīng)用最廣泛的誤差反向傳播神經(jīng)網(wǎng)絡(luò)，簡稱BP網(wǎng)絡(luò)。這種神經(jīng)網(wǎng)絡(luò)具有1個由相互獨立的許多神經(jīng)元構(gòu)成的輸入層，用于輸入相關(guān)數(shù)據(jù)信息；1個或多個隱層，用于分析、模擬過程中復(fù)雜的中間計算；1個輸出層，輸出運算、模擬結(jié)果。

輸入層N個節(jié)點分別對應(yīng)N個輸入分量構(gòu)成的輸入向量X=(x1,x2,…,xN)，輸出層L個節(jié)點對應(yīng)L個輸出分量構(gòu)成的輸出向量Y=(y1,y2,…,yL)。隱含節(jié)點Zj與輸出節(jié)點Yk分別為：

BP算法的學(xué)習過程由正向傳播和反向傳播過程組成。在正向傳播過程中，輸入信息從輸入層經(jīng)隱含層處理，并傳向輸出層，每一層神經(jīng)元的狀態(tài)只影響下一層神經(jīng)元的。如輸出層得不到期望的輸出，則轉(zhuǎn)入反向傳播，將誤差信號沿原連接通道返回，通過修改各層神經(jīng)元的權(quán)值，使得誤差信號最小。

2 油藏埋存CO2影響參數(shù)的無因次化及評價指標的建立

為對油藏埋存CO2的效果進行評價和預(yù)測，引入CO2的埋存系數(shù)St的概念。CO2的埋存量為累積注入量減去累積采出量：

QsCO2=NiCO2-NpCO2

埋存系數(shù)，其數(shù)學(xué)表達式為埋存量與累積注入量的比值：

式中：QsCO2為CO2氣體的埋存量，104m3；NiCO2為CO2的累積注入量，104m3；NpCO2為CO2的累積采出量，104m3；St為CO2埋存系數(shù)。

3 CO2埋存系數(shù)神經(jīng)網(wǎng)絡(luò)預(yù)測模型的建立

3.1 BP神經(jīng)網(wǎng)絡(luò)模型的設(shè)計

3.2 學(xué)習訓(xùn)練樣本的建立

用訓(xùn)練樣本集訓(xùn)練網(wǎng)絡(luò)，即讓網(wǎng)絡(luò)自動尋找該函數(shù)的內(nèi)部表達形式，從而實現(xiàn)映射。由于學(xué)習樣本的離散性，學(xué)習樣本的選擇將直接影響埋存系數(shù)預(yù)測的結(jié)果。本文將提出的5個無因次變量賦值后進行組合數(shù)值模擬，結(jié)果見表1。將數(shù)值模擬結(jié)果作為BP神經(jīng)網(wǎng)絡(luò)的學(xué)習樣本集，結(jié)果見表2。用誤差反向傳播算法自動找出5個無因次變量和埋存系數(shù)之間的內(nèi)部表達形式，據(jù)此可以預(yù)測油藏埋存CO2的效果。

表1 試驗因素及取值水平

表2 神經(jīng)網(wǎng)絡(luò)學(xué)習訓(xùn)練樣本數(shù)據(jù)

續(xù)表2

3.3 神經(jīng)網(wǎng)絡(luò)學(xué)習及預(yù)測效果評價

訓(xùn)練誤差曲線見圖1。

圖1 訓(xùn)練誤差曲線

從圖1看出，由訓(xùn)練將學(xué)習樣本輸入后，經(jīng)過數(shù)10 000次誤差反向傳導(dǎo)迭代運算，達到了很好的訓(xùn)練效果，誤差控制在0.001內(nèi)。

通過網(wǎng)絡(luò)訓(xùn)練得到網(wǎng)絡(luò)權(quán)值和閾值后，再輸入預(yù)測樣本進行網(wǎng)絡(luò)預(yù)測，可檢驗本文所建模型的正確性。預(yù)測樣本的選取為5個無因次變量重新組合后數(shù)模計算的不同于學(xué)習樣本的15組數(shù)據(jù)，結(jié)果見表3?？傮w上來看，多數(shù)實際輸出值與期望輸出值的誤差率在8%的范圍內(nèi)，數(shù)值模擬計算埋存系數(shù)與神經(jīng)網(wǎng)絡(luò)預(yù)測埋存系數(shù)的相關(guān)系數(shù)為0.94，這對于油藏埋存CO2的效果預(yù)測是可允許且較精確的誤差范圍。因此，文中建立神經(jīng)網(wǎng)絡(luò)預(yù)測模型在預(yù)測油藏埋存CO2埋存系數(shù)時是切實可行的。

表3 神經(jīng)網(wǎng)絡(luò)預(yù)測樣本數(shù)據(jù)

4 結(jié)論

1)油藏埋存CO2效果影響因素與埋存系數(shù)之間是一種非線性、不確定的復(fù)雜關(guān)系， 采用常規(guī)的實驗及多元線性回歸方法不能滿足實際預(yù)測的需要。

2)BP網(wǎng)絡(luò)預(yù)測結(jié)果表明，預(yù)測值與實際值較為吻合，誤差相對較小，預(yù)測模型可靠，訓(xùn)練好的網(wǎng)絡(luò)具有良好的簡單性和通用性。該神經(jīng)網(wǎng)絡(luò)提高了油藏埋存CO2埋存效果預(yù)測的準確性和效率，提供了一種新的方法。

[1] 沈平平,楊永智.溫室氣體在石油開采中資源化利用的科學(xué)問題[J].中國基礎(chǔ)科學(xué),2006,8(3):23-31.

[2] 何俊,陳小凡,樂平,等.線性回歸方法在油氣產(chǎn)量遞減分析中的應(yīng)用[J].巖性油氣藏,2009,21(2):103-105.

[3] 李先鵬.膠結(jié)指數(shù)的控制因素及評價方法[J].巖性油氣藏,2008,20(4):105-108.

[4] 孫勤華,劉曉梅,劉建新,等.利用波形分析技術(shù)半定量預(yù)測塔中碳酸鹽巖儲層[J].巖性油氣藏,2010,22(1):101-103.

[5] 李虎,蒲春生,吳飛鵬.基于廣義回歸神經(jīng)網(wǎng)絡(luò)的CO2驅(qū)最小混相壓力預(yù)測[J].巖性油氣藏,2012,24(1):108-111.

[6] 楊建,楊程博,張巖,等.基于改進神經(jīng)網(wǎng)絡(luò)的滲透率預(yù)測方法[J].巖性油氣藏,2011,23(1):98-102.

[7] 王濤,姚約東,朱黎明,等.油藏埋存二氧化碳影響因素Box-Behnken法優(yōu)化研究[J].中國石油大學(xué)勝利學(xué)院學(xué)報,2009,23(4):5-8.

[8] 戴文戰(zhàn).基于三層BP網(wǎng)絡(luò)的多指標綜合評估方法及應(yīng)用[J].系統(tǒng)工程理論與實踐,1999,19(53):29-34.