范利軍，劉佳佳，李憲騰，趙東亞，李兆敏，鹿 騰，楊建平

1.中國(guó)石油大學(xué)（華東）化學(xué)工程學(xué)院，山東青島 266580

2.中國(guó)石油大學(xué)（華東）石油工程學(xué)院，山東青島 266580

3.國(guó)家能源稠（重）油開(kāi)采研發(fā)中心遼河油田公司，遼寧盤(pán)錦 124000

煙道氣驅(qū)分為混相驅(qū)和非混相驅(qū)，本文考慮的是混相驅(qū)情況。煙道氣、溶劑與稠油作用后，原油黏度減小，表面張力降低，邦德數(shù)增大，從而降低采出井中非油相組分的含量，進(jìn)一步提高了原油采收率。

對(duì)于油藏注采的優(yōu)化，通常采用指標(biāo)對(duì)比法。即通?；谟筒?cái)?shù)值模擬軟件（CMG、VIP和Eclipse），在建立常規(guī)油藏驅(qū)油模型后，對(duì)分別給定的不同開(kāi)采方案進(jìn)行油藏?cái)?shù)值模擬，然后根據(jù)定義的指標(biāo)選擇相對(duì)較好的開(kāi)發(fā)方案。但是所有油藏?cái)?shù)值模擬軟件的本質(zhì)是數(shù)學(xué)模型的集成，是大型商業(yè)公司封裝好的商業(yè)產(chǎn)品，只能單方面進(jìn)行模擬，無(wú)法結(jié)合一般的優(yōu)化方法進(jìn)行開(kāi)源優(yōu)化。而混相煙道氣驅(qū)的數(shù)學(xué)模型涉及一系列的耦合偏微分方程和許多物理化學(xué)代數(shù)方程，問(wèn)題復(fù)雜度高，計(jì)算復(fù)雜，難以操作。因此，本文基于CMG軟件的采樣數(shù)據(jù)，采用模糊神經(jīng)網(wǎng)絡(luò)的建模方法對(duì)混相煙道氣驅(qū)進(jìn)行數(shù)據(jù)建模，以便基于該模型進(jìn)行開(kāi)采方案優(yōu)化，給出具體的最優(yōu)開(kāi)采策略。

1 模糊神經(jīng)網(wǎng)絡(luò)的建模方法

神經(jīng)網(wǎng)絡(luò)技術(shù)是一種被廣泛用于非線(xiàn)性系統(tǒng)建模中的基于數(shù)據(jù)的學(xué)習(xí)訓(xùn)練方法。然而，在神經(jīng)網(wǎng)絡(luò)系統(tǒng)建模中，建模的精度和模型的泛化能力是一對(duì)廣泛存在的矛盾體[1]。出現(xiàn)上述問(wèn)題的一個(gè)重要原因是神經(jīng)網(wǎng)絡(luò)是僅僅基于數(shù)據(jù)的，缺乏經(jīng)驗(yàn)知識(shí)的結(jié)合。而在模糊系統(tǒng)中，模糊集、隸屬度函數(shù)、模糊規(guī)則的設(shè)計(jì)是建立在經(jīng)驗(yàn)知識(shí)基礎(chǔ)上的，能夠避免純客觀(guān)因素帶來(lái)的泛化性不足的問(wèn)題[2-7]。大量文獻(xiàn)表明，Takagi-Sugeno（T-S）模型具有諸多優(yōu)異性能，比如計(jì)算簡(jiǎn)單，易于優(yōu)化等[8-12]。

本文基于T-S模型的模糊神經(jīng)網(wǎng)絡(luò)，采用混合螺旋優(yōu)化算法（HSO）訓(xùn)練網(wǎng)絡(luò)權(quán)值以及隸屬度函數(shù)的相關(guān)參數(shù)。

1.1 模糊系統(tǒng)的T-S模型

由于多輸入多輸出（MIMO）的模糊規(guī)則可分解為多個(gè)多輸入單輸出（MISO）模糊規(guī)則，因此下面只討論MISO模糊系統(tǒng)的模型[13]。

設(shè)輸入向量，x=[x1，x2，…，xn]T，式中每個(gè)分量 xi均為模糊語(yǔ)言變量[6，13-16]。Takagi-Sugeno 所提出的模糊規(guī)則后件是輸入變量的線(xiàn)性組合，模糊系統(tǒng)的輸出量為每條規(guī)則的輸出量的加權(quán)平均，即[14-16]：

式中：y為模糊系統(tǒng)的輸出量；αi為權(quán)重；yj為各規(guī)則后件的輸出量；為權(quán)重平均值；m為整數(shù)。

1.2 系統(tǒng)結(jié)構(gòu)

該網(wǎng)絡(luò)是由前件網(wǎng)絡(luò)和后件網(wǎng)絡(luò)兩部分組成的，前件網(wǎng)絡(luò)用來(lái)匹配模糊規(guī)則，后件網(wǎng)絡(luò)用來(lái)產(chǎn)生模糊規(guī)則[15-19]。

（1）前件網(wǎng)絡(luò)。前件網(wǎng)絡(luò)由四層組成。第一層為輸入層，它的每一個(gè)節(jié)點(diǎn)直接與輸入向量的各分量xi連接，其作用為將輸入值x=[x1，x2，…，xn]T傳送到下次層[13-19]；第二層的每個(gè)節(jié)點(diǎn)代表一個(gè)語(yǔ)言變量值，如NM、PS等，其作用是計(jì)算各輸入分量屬于各語(yǔ)言變量值模糊集合的隸屬函數(shù)[13-19]；第三層的每個(gè)節(jié)點(diǎn)代表一條模糊規(guī)則，其作用是匹配模糊規(guī)則的前件，并計(jì)算出每一條規(guī)則的適應(yīng)度；第四層所實(shí)現(xiàn)的是歸一化計(jì)算。

（2）后件網(wǎng)絡(luò)。后件網(wǎng)絡(luò)由各結(jié)構(gòu)相同的并列子網(wǎng)絡(luò)組成，每個(gè)子網(wǎng)絡(luò)產(chǎn)生一個(gè)輸出量[18-19]。第一層為子網(wǎng)絡(luò)的輸入層，它將輸入變量傳送到第二層，其作用是提供模糊規(guī)則后件中的常數(shù)項(xiàng)[14-19]；子網(wǎng)絡(luò)的第二層共有m個(gè)節(jié)點(diǎn)，每個(gè)節(jié)點(diǎn)代表一條規(guī)則，該層的作用是計(jì)算每一條規(guī)則的后件；子網(wǎng)絡(luò)的第三層是計(jì)算系統(tǒng)的輸出。

1.3 參數(shù)學(xué)習(xí)

上述過(guò)程即為T(mén)-S模糊系統(tǒng)模型。本文采用混合螺旋優(yōu)化算法訓(xùn)練模糊神經(jīng)網(wǎng)絡(luò)的權(quán)值，以及隸屬度函數(shù)的參數(shù)[18-20]。

假設(shè)T-S模型中前件第二層的隸屬度函數(shù)為采用高斯函數(shù)表示的鈴型函數(shù)，即：

式中：cij和σij分別表示隸屬函數(shù)的中心和寬度[15-19]。

對(duì)于本文中的基于T-S模型的模糊神經(jīng)網(wǎng)絡(luò)，需要學(xué)習(xí)的參數(shù)包括前件網(wǎng)絡(luò)第二層中高斯函數(shù)的中心cij和寬度σij，后件網(wǎng)絡(luò)第二層中的權(quán)重。

2 基于模糊神經(jīng)網(wǎng)絡(luò)的煙道氣驅(qū)建模

2.1 油藏描述

油藏采用四注一采方式，即模型中包含4口注入井和1口生產(chǎn)井，其井位分布情況如圖1所示。

圖1四注一采井位示意

油藏厚5 m，長(zhǎng)421.02 m，寬443.8 m，按x、y、z三個(gè)方向?qū)τ筒貏澐志W(wǎng)格，x方向網(wǎng)格邊長(zhǎng)46.78 m，y方向網(wǎng)格邊長(zhǎng)44.38 m，總供劃分為9×10網(wǎng)格。油藏上表面距離地表的深度為1 200 m，上層孔隙度為0.31，下層孔隙度為0.32，孔隙體積為5.9 m3。油藏初始含油飽和度為0.65，初始含水飽和度0.35，原油密度為0.95 g/cm3，地質(zhì)中原油儲(chǔ)量為3.634×105t，油藏原始?jí)毫?2 MPa，非均勻滲透。滲透率如圖2所示。

440C為美國(guó)牌號(hào)，近似于我國(guó)的102Cr17Mo，是一種高碳高鉻的馬氏體不銹鋼，其碳含量（質(zhì)量分?jǐn)?shù)）達(dá)到1%左右。該鋼適用于制造在腐蝕環(huán)境和無(wú)潤(rùn)滑強(qiáng)氧化氣氛中工作的軸承、軸套，燃油附件中的活門(mén)襯套和其他既耐磨又耐蝕的零件。由于其高碳高鉻的特性，熱處理后碳化物數(shù)量多，耐磨性好，耐蝕性好。在大氣、水、海水中以及某些酸類(lèi)和鹽類(lèi)的水溶液中有優(yōu)良的耐蝕性。該鋼由于含碳量高，生產(chǎn)時(shí)易出現(xiàn)脫碳、淬裂和殘留奧氏體等問(wèn)題，是較難進(jìn)行熱處理的鋼種。

圖2 絕對(duì)滲透率分布/mD

當(dāng)采出井液中非油驅(qū)替相質(zhì)量分?jǐn)?shù)fa達(dá)到97%時(shí)，開(kāi)始從4口注入井中注入煙道氣，初始含水飽和度如圖3所示，初始?jí)毫θ鐖D4所示。油藏中煙道氣的初始體積分?jǐn)?shù)為0。油藏?cái)?shù)據(jù)及流體數(shù)據(jù)分別見(jiàn)表1和表2。注入方式采用三段塞，4口注入井的注入量Qin均為15 m3/d，生產(chǎn)井產(chǎn)液量Qout為60 m3/d[20]。采用CMG軟件（2010版本）進(jìn)行數(shù)值模擬，得到樣本數(shù)據(jù)。

圖3 初始含水飽和度分布

2.2 段塞注入策略

對(duì)煙道氣驅(qū)模型，系統(tǒng)輸入為注入井的煙道氣注入速度，這里只考慮注入速度變化，假設(shè)注入壓力、段塞固定，驅(qū)替時(shí)間固定；系統(tǒng)的輸出為采出井的非油相組分質(zhì)量分?jǐn)?shù)fa，為了方便，由非油相組分求得采油量。煙道氣的注入采用段塞形式，注入速度可按下式表示：

圖4 初始?jí)毫Ψ植?MPa

表1 油藏參數(shù)

表2 流體數(shù)據(jù)

式中：v（t）為注入速度，m3/d；v（wk）為第w口注入井在第k個(gè)段塞長(zhǎng)度內(nèi)注入驅(qū)替液中煙道氣的注入速度，m3/d；t為段塞時(shí)長(zhǎng)，月；k為段塞數(shù)；P為最大段塞數(shù)；tp、tf分別表示注入速度為零的段塞上、下限。

為了建立煙道氣驅(qū)的辨識(shí)模型，基于油藏?cái)?shù)值模擬軟件CMG進(jìn)行了油藏?cái)?shù)值模擬，得到注入井煙道氣注入速度v與采出井非油相組分質(zhì)量分?jǐn)?shù)fa間的曲線(xiàn)關(guān)系樣本數(shù)據(jù)，而后用于基于HSO和T-S模型的模糊神經(jīng)網(wǎng)絡(luò)的建模。

2.3 建模過(guò)程

采用三段塞注入策略，注采方案為先水驅(qū)后煙道氣驅(qū)。在2005年12月31日至2013年12月31日時(shí)間段內(nèi)，按月采樣數(shù)據(jù)，共取97個(gè)時(shí)間點(diǎn)，煙道氣和溶劑的注入時(shí)間為2006年10月1日至2010年9月1日。注入煙道氣段塞的時(shí)間節(jié)點(diǎn)分別為10～21月、22～50月、51～56月，具體注入方式如圖5所示。

圖5 三段塞劃分

為了使建模更具說(shuō)服力，在煙道氣注入速度v∈[5，30]（m3/d） 范圍內(nèi)隨機(jī)生成100組注入策略，分別運(yùn)用到CMG建模的上述的油藏區(qū)塊進(jìn)行數(shù)值模擬，得到與100組注入策略相應(yīng)的采出井非油相組分質(zhì)量分?jǐn)?shù)fa，即得到了100組輸入、輸出樣本數(shù)據(jù)。

采用基于HSO和T-S模型的模糊神經(jīng)網(wǎng)絡(luò)進(jìn)行訓(xùn)練建模，從而建立了煙道氣驅(qū)油藏的辨識(shí)模型，得到注入井煙道氣注入速度和采出井非油組分質(zhì)量分?jǐn)?shù)之間的關(guān)系，其數(shù)學(xué)描述如下：

式中：Ψ（v）為復(fù)雜的非線(xiàn)性函數(shù)，代表模糊神經(jīng)網(wǎng)絡(luò)訓(xùn)練的結(jié)果。

由于把100組建模數(shù)據(jù)畫(huà)在一個(gè)圖中無(wú)法清晰顯示，因此，從100組非油相組分質(zhì)量分?jǐn)?shù)fa中選取9組，對(duì)建模的效果進(jìn)行測(cè)試。選取的9組輸入（注入井的煙道氣注入速度v）如表3所示。

注采策略的建模輸出（模糊神經(jīng)網(wǎng)絡(luò)模型的計(jì)算結(jié)果）與實(shí)際輸出（CMG模型的模擬計(jì)算結(jié)果）都十分接近，對(duì)于整個(gè)驅(qū)油時(shí)間段，建模誤差均小于2.2%，且絕大部分誤差都小于1%，見(jiàn)圖6～7。因此，本文提出的基于HSO和T-S模型的模糊神經(jīng)網(wǎng)絡(luò)，可建立準(zhǔn)確的煙道氣驅(qū)模型。

為了測(cè)試模型的泛化能力，在[5，30]（m3/d）范圍內(nèi)隨機(jī)選取一組注入策略，分別代入CMG軟件和建立的模型中，得到模型輸出。若選取4口注入井的煙道氣注入速度為{（30，20，12），（25，15，10），（25，15，10），（20，10，6），（10，12，19）}m3/d，則其相應(yīng)的測(cè)試輸出曲線(xiàn)如圖8所示，測(cè)試輸出誤差曲線(xiàn)如圖9所示。

從圖8～9中可以發(fā)現(xiàn)，模型的輸出基本和實(shí)際值相吻合。定義如下指標(biāo)為模型輸出與實(shí)際值的相對(duì)誤差的均值st：

表3 煙道氣注入策略/（m3·d-1）

圖6 9種策略的模糊神經(jīng)網(wǎng)絡(luò)建模輸出結(jié)果和實(shí)際值

圖7 9種策略的建模計(jì)算誤差

式中：fai表示模型輸出采出井非油組分質(zhì)量分?jǐn)?shù)，%；表示實(shí)際采出井非油組分質(zhì)量分?jǐn)?shù)，%；i表示模型計(jì)算次數(shù)；N表示模型計(jì)算最大次數(shù)。

通過(guò)計(jì)算可知，平均相對(duì)誤差st=0.011 3%，其中含水率的單位為%，相對(duì)誤差絕對(duì)值的最大值為1.14%，模型輸出與真實(shí)值的差值的方差為0.125。從上述數(shù)據(jù)可知，該模型的泛化能力較好，與實(shí)際值十分接近。

綜上所述，基于HSO和T-S模型的模糊神經(jīng)網(wǎng)絡(luò)建模方法可有效地用于煙道氣驅(qū)油藏的建模，它具有較高的精度和可靠性。

3 結(jié)束語(yǔ)

圖8 模型測(cè)試輸出結(jié)果

圖9 模型的測(cè)試誤差

本文結(jié)合神經(jīng)網(wǎng)絡(luò)、模糊系統(tǒng)和智能算法提出了一種基于HSO和T-S模型的模糊神經(jīng)網(wǎng)絡(luò)建模方法。本方法繞過(guò)了煙道氣驅(qū)的機(jī)理模型，避免了復(fù)雜的滲流方程和能量方程的計(jì)算，避免了帶耦合的高階偏微分方程組的計(jì)算，簡(jiǎn)化了煙道氣驅(qū)的模型，減小了對(duì)于先驗(yàn)知識(shí)的依賴(lài)性。數(shù)值仿真計(jì)算驗(yàn)證了該建模方法在建立煙道氣驅(qū)模型上的準(zhǔn)確性和可靠性。

[1]劉毅娟，雷鳴，何旸，等.人工神經(jīng)網(wǎng)絡(luò)概述[J].電子測(cè)試，2015（6）：74-75.

[2]RAO S S，CHEN L.Numericalsolution of fuzzy linear equations in engineering analysis[J].International Journal for Numerical Methods in Engineering，2015，43（3）：391-408.

[3]TALAEI M，MOGHADDAM B F，PISHVAEE M S，et al.A robust fuzzy optimization model for carbon-efficient closed-loop supply chain network design problem:a numerical illustration in electronics industry [J].Journal of Cleaner Production，2015，113：662-673.

[4]YAGHOOBI A，BAKHSHI-JOOYBARI M，GORJI A，et al.Application of adaptive neuro fuzzy inference system and genetic algorithm for pressure path optimization in sheet hydroforming process[J].International Journal of Advanced Manufacturing Technology，2016：1-11.

[5]DENG Y，REN Z，KONG Y，et al.Ahierarchicalfused fuzzy deep neural network for data classification[J].IEEE Transactions on Fuzzy Systems，2016：1-1.

[6]張芳.基于改進(jìn)GA的模糊神經(jīng)網(wǎng)絡(luò)參數(shù)學(xué)習(xí)算法及其應(yīng)用的研究[D].長(zhǎng)沙：中南大學(xué)，2007.

[7]唐海燕.基于模糊神經(jīng)網(wǎng)絡(luò)的電梯群優(yōu)化控制研究[D].哈爾濱：哈爾濱工業(yè)大學(xué)，2006.

[8]劉駿，殷曉明，顧幸生.一種改進(jìn)的T-S模糊模型建模及優(yōu)化方法[J].華東理工大學(xué)學(xué)報(bào)（自然科學(xué)版），2016，42（2）：233-239.

[9]張健.模糊神經(jīng)網(wǎng)絡(luò)模型算法研究與應(yīng)用[D].大慶：大慶石油學(xué)院，2002.

[10]孫海蓉.模糊神經(jīng)網(wǎng)絡(luò)的研究及其應(yīng)用[D].保定：華北電力大學(xué)（河北），2006.

[11]何春梅.模糊神經(jīng)網(wǎng)絡(luò)的性能及其學(xué)習(xí)算法研究[D].南京：南京理工大學(xué)，2010.

[12]周琛暉.模糊神經(jīng)網(wǎng)絡(luò)穩(wěn)定性分析及其應(yīng)用[D].成都：電子科技大學(xué)，2013.

[13]張妍妍.連續(xù)攪拌反應(yīng)釜的智能控制器設(shè)計(jì)[D].鄭州：鄭州大學(xué)，2007.

[14]李國(guó)勇.智能預(yù)測(cè)控制及其MATLAB實(shí)現(xiàn) [M].第2版.北京：電子工業(yè)出版社，2010.

[15]劉坤.電液伺服系統(tǒng)的智能控制研究 [D].秦皇島：燕山大學(xué)，2003.

[16]楊勇.EMD和模糊神經(jīng)網(wǎng)絡(luò)在滾動(dòng)軸承故障診斷中的研究與應(yīng)用[D].太原：太原理工大學(xué)，2008.

[17]秦曉娜.基于T-S模糊神經(jīng)網(wǎng)絡(luò)的CFB-FGD過(guò)程脫硫效率的預(yù)報(bào)研究[D].沈陽(yáng)：東北大學(xué)，2009.

[18]畢作文.模糊神經(jīng)網(wǎng)絡(luò)在CASS工藝污水處理系統(tǒng)中的研究與應(yīng)用[D].蘭州：蘭州理工大學(xué)，2011.

[19]郝敏.基于模糊聚類(lèi)算法的自適應(yīng)模糊神經(jīng)網(wǎng)絡(luò)研究[D].哈爾濱：哈爾濱理工大學(xué)，2007.

[20]雷陽(yáng).高溫高鹽油藏聚合物驅(qū)最優(yōu)控制方法研究[D].青島：中國(guó)石油大學(xué)（華東），2013.