利用ANFIS模型預(yù)測吸水剖面.docx

1、斷塊油氣田FAULT-BLOCKOIL&GASFIELD文章編號：1005-8907(2009)04-072-03利用ANFIS模型預(yù)測吸水剖面單玲李慧莉高秀田（中原油田分公司勘探開發(fā)科學(xué)研究院，河南濮陽457001）基金項目：中國石油天然氣集團公司石油科技中青年創(chuàng)新菽金“高金水期油藏大孔道動態(tài)變化機理與識別模型"（04E7029）資助摘要隨著開發(fā)時間的延長，嚴重非均質(zhì)油藏層間干擾越來越嚴重，需要對開發(fā)屢系進行調(diào)整重組。根據(jù)油藏動、靜態(tài)資料，對胡7南油藏沙三中的嚴重非均質(zhì)儲層的滲透率、級差分布進行合理選取，設(shè)計長巖心水驅(qū)油實茲，通過實驗確定不同滲透率巖心動用規(guī)律；在實膾基礎(chǔ)上

2、利用修正的數(shù)值模擬模型進行多層聯(lián)合水驅(qū)動用條件下的分層動用規(guī)律研霓；提出用綜合影響因子來反映非均步性對低滲透層的動用特征、驅(qū)替效果的定量影響程度；并找出不同滲透屈組合的合理界限，用以指導(dǎo)此類嚴重非均質(zhì)油喊.開發(fā)政策的制定。關(guān)鍵詞非均質(zhì)；開發(fā)層系；疹透率級差；組合界限中圖分類號:TE319文獻標(biāo)識碼：AForecastingiijectionprofilebyANFISmodelShanLingLiHuiliGaoXiutian(Exploration&DevelopmentResearchInstitute,ZhongyuanOilfieldCompany,SINOPEC,Puyang

3、,457001,China)Asthedevelopmenttimegoingon,theinterlayerinterferenceofextremelyheterogeneousreservoirbecomesmoreserious.Thelayerseriesofdevelopmentshouldbeadjusted.Basedonthedynamicandstaticreservoirdata,reasonableselectionsaremadeforthepermeabilityandgradingdistributionoftheextremelyheterogeneousres

4、ervoirinthemiddlememberofthethirdsectionofShahejieFormationinthesouthareaofHu7BlockWaterdisplacingoiloflongcoreexperimentsisdesigned,whichisusedtodetermineregularproducingpatternofcoreswithdifferentpermeabilities,basedontheexperiment,layeredproducingpatternwasstudiedundertheconditionofdevelopmentwit

5、hmultilayersseriesbyusingthemodifiednumericalsimulator.Itisputforwardthatthequantitativeinfluenceofheterogeneitytoproductioncharacteristicsanddisplacementeffectsoflayerswithlowpermeabilitywillbereflectedbyacombinedfactor.Areasonabletechnicallimitforcombinationsofdifferentpermeablelayersisidentifiedi

6、nordertoguidetheformulationofdevelopmentpoliciestosuchextremelyheterogeneousreservoirs.Keywords:heterogeneity,developmentlayers,permeabilitygrading,combinationlimit.1吸水剖面影響因素分析注水井吸水剖面的變化受多種因素影響，根據(jù)現(xiàn)場實際情況,初步分析認為，影響注水井吸水剖面變化的主要因素為注水井各小層滲透率、小層有效厚度、地層原油黏度、注水壓力等。同時,考慮到吸水剖面的變化同樣受注采系統(tǒng)的綜合影響，因此將注水井對應(yīng)見效油井的相關(guān)參數(shù)

7、也作為影響因素，其中包括對應(yīng)油井點的各層滲透率、有效厚度、油井產(chǎn)液雖及油水井井距。各影響因素對注水井各小層相對吸水量:變化的影響程度各不相同，為了定量了解其影響程度，采用層次分析法計算這些因素對吸水剖面影響的權(quán)數(shù)分配，以各自的權(quán)數(shù)作為因素相對重要性的評價指標(biāo)。1.1建立遞階層次結(jié)構(gòu)模型及判斷矩陣應(yīng)用層次分析方法定量評價注水井吸水剖面變化的影響因素，需要建立體現(xiàn)影響吸水剖面因素之間相互關(guān)系的遞階層次結(jié)構(gòu)模型（見圖1）。圖1吸水剖面綜合評價遞階層次結(jié)構(gòu)模型收稿日期：2008-】2-26；改回日期；2009-05-12.作者簡介：單玲，女，1968年生，1988年畢業(yè)于K春地質(zhì)學(xué)院.現(xiàn)從事油氣田開發(fā)

8、研究工作。電話：（0393）4816164。層次結(jié)構(gòu)模型一經(jīng)建立，就確定了上下層因素之間的隸屬關(guān)系，進而可以構(gòu)造體現(xiàn)這種隸屬關(guān)系的判斷矩陣。判斷矩陣是通過兩兩比較來構(gòu)造的。12各因素相對要性定指標(biāo)計算1.2.1計算方法采用和積法的公式計算權(quán)重向缺W=（W|,W2,，WJ。計算步驟如下：1）將判斷矩陣每一列正規(guī)化：牛一（ij=l,2,n）*N12）每一列經(jīng)正規(guī)化后的判斷矩陣按行相加：n祐£插（i=l,2,n）（2）i13）對向最直取,電,，WT正規(guī)化得到的W=叫即為所求特征向量。4）計算矩陣最大特征根：(3)3粉j«I八綺式中：W為權(quán)向髭;W,為權(quán)向量元素；吟為判斷矩陣的正規(guī)

9、化矩陣;K為判斷矩陣;如為判斷矩陣元素為最大特征根；（AW）j為向AW的第£個元素。1.2.2判斷矩陣的一致性檢驗判斷矩陣的一致性檢驗公式為N其中幾一1式中/為一致性比例;N為平均一致性指標(biāo)湖為一致性指標(biāo)。當(dāng)£<0.1時，判斷矩陣具有滿意的一致性，否則，要重新建立判斷矩陣直到達到滿意的一致性為止。1.2.3計算各層因素相對于目標(biāo)層的合成權(quán)重計算各層因素相對于目標(biāo)層的合成權(quán)重，就是通過進行層次總排序計算出各層因素的權(quán)重。計算中需要對整個遞階層次結(jié)構(gòu)在K層及其以上水平的所有判斷矩陣進行一致性檢驗，檢驗公式為ic.乙=號（4）i=J式中:G為因素權(quán)重值。如果1x0.1,則認

10、為遞階層次結(jié)構(gòu)在K層及其以上的所有判斷矩陣具有滿意的一致性，計算出的合成權(quán)重合理A）。1.2.4計算結(jié)果根據(jù)區(qū)塊的實際情況，建立吸水剖面影響因素對于吸水剖面變化（各小層相對吸水量）的層次結(jié)構(gòu)模型，并利用上述和積法的計算結(jié)果及通過層次總排序，可以獲得各層因素的權(quán)重（其中=0.039801<0.1）（見表1）。«1吸水剖面變化影響因素權(quán)項目因素“權(quán)底W厚度G0.223注水井滲透率G0.278參數(shù)注入壓力G0.114原油黏度G0.022厚度G0.103見效油井漆透率G0.112參數(shù)產(chǎn)液量G0.043注采井距G0.105由計算結(jié)果看，注水井小層滲透率權(quán)重值最高，表明注水井的小層滲透率是

11、影響吸水剖面分布的最重要因素；其次是小層有效厚度；水井注入壓力、對應(yīng)見效油井的小層滲透率、有效厚度及注采井距等屬于中等程度影響因素；油井產(chǎn)液量及原油黏度屬于影響程度較低的因素。2吸水剖面變化預(yù)測2.1自適應(yīng)模糊神經(jīng)網(wǎng)絡(luò)ANFIS模型綜合利用了模糊系統(tǒng)規(guī)則的可解釋性和神經(jīng)網(wǎng)絡(luò)自適應(yīng)學(xué)習(xí)的特點，因此可以有效地解決非線性系統(tǒng)的建模與控制問題A）。ANFIS模型的隸屬度函數(shù)和模糊規(guī)則是通過實際測試的吸水剖面數(shù)據(jù)得到的。進行模糊推理需要以下4個步驟:1）劃分輸入空間和確定隸屬函數(shù)。2）利用相應(yīng)的Fuzzy算子（一般為乘法算子或取小算子）求出各條加則的適用度。3）產(chǎn)生每一規(guī)則的后續(xù)輸出，即T-S模型中后面

12、的部分。4）合并步驟3）中產(chǎn)生的所有輸出，形成一個模糊輸出，這一步也稱為“解模糊”。ANFIS模型的結(jié)構(gòu)是一個多層前饋網(wǎng)絡(luò)，各層之間的連接權(quán)值可以通過給定的樣本數(shù)據(jù)來進行自適應(yīng)調(diào)節(jié)，其參數(shù)學(xué)習(xí)采用混合學(xué)習(xí)算法：在前向傳遞時，輸入各組吸水剖面數(shù)據(jù)后，數(shù)據(jù)通過各節(jié)點函數(shù)計算網(wǎng)絡(luò)每層各節(jié)點的輸出，然后通過線性最小二乘估計對模糊規(guī)則后件的結(jié)論參數(shù)進行辨識，計算每組輸入數(shù)據(jù)的輸出誤差;在反向傳遞時，將誤差用最速下降梯度法從輸出節(jié)點反向傳遞到輸入節(jié)點，并對模糊規(guī)則的前件參數(shù)(即隸屬度函數(shù)形狀參數(shù))進行調(diào)整。如此反復(fù)循環(huán)，直到滿足誤差標(biāo)準或達到迭代次數(shù)為止。2.2吸水剖面變化研究結(jié)果.2.2.1指標(biāo)確定及樣

13、本數(shù)據(jù)采集為了借助ANFIS模型開展吸水剖面變化規(guī)律的研究,根據(jù)實際情況，選用注水井及對應(yīng)采油井相關(guān)參數(shù)作為模型的輸入樣本，包括注水井小層滲透率、有效厚度、注入壓力(換算到各小層)、采油井小層滲透率、有效厚度及注采井距等。將測試的注水井各小層相對吸水技作為輸出樣本。2.2.2模型模擬實現(xiàn)借助MATLAB的模糊工具箱提供的輔助自適應(yīng)模糊推力工具箱中的相關(guān)函數(shù)：anils、genfis2、evalfis等，建立并模擬實現(xiàn)基于自適應(yīng)神經(jīng)模糊推理系統(tǒng)的吸水剖面變化預(yù)測模型。計算過程中經(jīng)過對輸入/輸出數(shù)據(jù)以初始步長為0.5的180次迭代，得到經(jīng)過調(diào)整的模糊神經(jīng)網(wǎng)絡(luò)模型。圖2為訓(xùn)練樣本數(shù)據(jù)與ANFIS模型

14、輸出數(shù)據(jù)的比較圖。350從圖2可以看出，訓(xùn)練用的實際數(shù)據(jù)與模型輸出數(shù)據(jù)吻合較好，計算二者均方差為3.925x10-5,滿足擬合要求。選用某注水井未參與訓(xùn)練的吸水剖面數(shù)據(jù)進行檢驗，根據(jù)模型預(yù)測結(jié)果與實際值對比狀況,檢驗樣本的均方誤差為3.3288x10-5,誤差較小，表明利用預(yù)測模型精度較高，可拳性好。3利用吸水剖面進行層系優(yōu)化組合3.1層系優(yōu)化組合思路針對礦場實際悄況，層系優(yōu)化組合的主要思路是油藏的開發(fā)不再按照固定的層系劃分模式，而是從單并點入手，以注水井為中心逐井進行各小層的優(yōu)化重組，以實現(xiàn)改善吸水剖面的目標(biāo)。3.2利用吸水剖面進行層系優(yōu)化組合的主要做法利用吸水剖面進行層系優(yōu)化組合的具體做法

15、是，針對注水井單井，首先劃分大的注水層段，劃分后的層段需要保持一定的儲ht規(guī)模，層段內(nèi)各小層性質(zhì)接近，高滲層不宜過多，好差層間互的頻率不能過高，然后利用ANFIS模型預(yù)測劃分后各個小層的吸水剖面分布情況。如果預(yù)測結(jié)果表明吸水剖面分布較為均勻，則確定該劃分方式；如果預(yù)測結(jié)果表明吸水剖面仍然不均勻，則更新劃分注水層段,最終達到各小層較為均勻吸水的目的。4結(jié)束語量化分析了影響注水井吸水剖面的地質(zhì)開發(fā)因素，提出了將重點影響因素作為輸入?yún)?shù)并利用模糊神經(jīng)網(wǎng)絡(luò)進行吸水剖面預(yù)測的方法。該吸水剖面預(yù)測法可以預(yù)測出不同地質(zhì)開發(fā)條件下注水井的吸水剖面變化，從而可以量化指導(dǎo)注水井的分層注水優(yōu)化組合措施。參考文獻tn張昌民.儲層研究中的層次分析方法J.石油與天然氣地質(zhì)，1992,13(3):344-350.2 王蓮芬，許樹柏.層次分析法引論Mj.北京：中國人民大學(xué)出版社,1990.3 趙煥臣，許樹柏，和金生.層次分析法M.北京：科學(xué)出版社,1986.4 ZhangZhiyong.AreexaminationofmethodsofhierarchiccompositionintheAHPJJournalofSystemsScienceandSystemsEngineerin