本發(fā)明涉及一種識別井間連通性及油水動態(tài)預(yù)測的方法，具體涉及一種結(jié)合地層及流體屬性資料和注采動態(tài)數(shù)據(jù)資料，并借助優(yōu)化算法進(jìn)行相關(guān)參數(shù)反演進(jìn)而識別油水井連通性及油水動態(tài)預(yù)測的新方法，屬于油氣開發(fā)
技術(shù)領(lǐng)域：
。
背景技術(shù)：
：目前，我國大多數(shù)老油田均已進(jìn)入開發(fā)中后期，地下滲流情況復(fù)雜，注水效率不高，經(jīng)常出現(xiàn)注入水竄流等嚴(yán)重問題，從而導(dǎo)致油井開采期間含水率迅速升高，最終采收率難以提高。油藏井間連通性研究是油田精細(xì)注水及開發(fā)模式轉(zhuǎn)換的基礎(chǔ)，也是老油田實(shí)現(xiàn)低成本開發(fā)的靈魂，可以指導(dǎo)老油田注采結(jié)構(gòu)調(diào)整、分析剩余油分布規(guī)律、制定加密調(diào)剖、研究轉(zhuǎn)注時機(jī)等措施方案。此外，由于我國石油勘探開發(fā)逐漸向老油區(qū)外圍及地層復(fù)雜區(qū)域發(fā)展，油藏地質(zhì)構(gòu)造及地層屬性越來越復(fù)雜，地層非均質(zhì)性嚴(yán)重，局部區(qū)域出現(xiàn)斷層、尖滅等復(fù)雜情況，導(dǎo)致進(jìn)一步勘探開發(fā)難度加大。而井間連通性的研究可以幫助指導(dǎo)復(fù)雜區(qū)域的油田勘探開發(fā)，大大減小開發(fā)成本，提高采收率，最終提升低油價下油田經(jīng)濟(jì)效益。因此如何有效利用現(xiàn)有的數(shù)據(jù)，經(jīng)濟(jì)高效地評價油藏井間連通性即顯得尤為重要。目前，油田所采用的井間連通性的評價方式需要耗費(fèi)大量的時間和財力去獲取詳細(xì)的數(shù)據(jù)，有些數(shù)據(jù)的獲取相當(dāng)困難或者難以精確測量，使得井間連通性無法定量評價或者不能準(zhǔn)確評價，加大開發(fā)成本。而井的注采動態(tài)數(shù)據(jù)的獲取卻相當(dāng)簡單，也便于測量，通過簡單的模型便可建立起注采井間的動態(tài)連通關(guān)系，為油田開發(fā)提供有效指導(dǎo)?，F(xiàn)有技術(shù)目前，油田常用的井間連通性評價方法主要包括靜態(tài)分析法和動態(tài)分析法兩大類。而油水動態(tài)指標(biāo)計算則主要采用經(jīng)驗(yàn)公式法，難以滿足各種類型油藏等復(fù)雜情況?，F(xiàn)有技術(shù)種類：1)電測井曲線對比通過油水井的電測井曲線對比可以定性判斷井的連通狀況。由于不同地層的巖石類型、沉積韻律及組合特征各不相同，反映在電測井曲線上也具有不同的特征。通過對比兩口井的電測井曲線，即可推斷出其構(gòu)造深度、沉積特征、巖石類型等是否在同一層位，進(jìn)而定性判斷兩口井之間是否具有連通關(guān)系。2)各井的油藏參數(shù)比較法通過對比各井的原油密度、組分等屬性參數(shù)定性判斷井的連通狀況。例如對油田某一井組相鄰六口井的原油進(jìn)行化學(xué)檢測，可知，這6口井的原油密度為(0.95～0.96)×103kg/m3，原油組分也相似，同時又處于同一水動力學(xué)系統(tǒng)，因此可以推斷這六口井的連通狀況較好。3)試井分析方法試井是認(rèn)識油藏的一種重要方法，也可用于判斷井間連通性，常用的試井方法主要有不穩(wěn)定試井、干擾試井和脈沖試井。干擾試井是通過改變一口井的工作制度，同時在周圍幾口井底裝入高精度壓力計，進(jìn)而觀測井底壓力的變化來判斷井間連通狀況。該方法雖然比較準(zhǔn)確，但是需要改變井的工作制度，影響油田正常經(jīng)營，同時下入壓力計成本也較高。4)示蹤劑測試方法該方法是油田最常用且最直觀的判斷井間連通性的方法。通過在注入流體中加入示蹤劑，同時跟蹤示蹤劑軌跡，在采油井中檢測示蹤劑含量、峰值濃度及變化等參數(shù)，結(jié)合數(shù)值模擬及生產(chǎn)動態(tài)分析判斷儲層特性。5)地球化學(xué)方法應(yīng)用色譜指紋技術(shù)可判別井間連通狀況。通過對原油進(jìn)行全烴色譜檢測得到油藏各層原油的基礎(chǔ)指紋數(shù)據(jù)庫，再對周圍井各小層原油進(jìn)行色譜檢測，判別原油性質(zhì)相似的儲層，若特征相似，則兩口井之間連通。6)數(shù)值模擬方法數(shù)值模擬方法通常與示蹤劑測試結(jié)果相結(jié)合分析井的連通性，首先通過示蹤劑測試得到峰值濃度、突破時間等參數(shù)，然后通過數(shù)值模擬軟件擬合這些參數(shù)，進(jìn)而驗(yàn)證連通性的分析結(jié)果?，F(xiàn)有技術(shù)缺點(diǎn)：以上現(xiàn)有方法實(shí)施起來比較困難，成本大，周期長。對于小層厚度較薄，油層、水層和油水同層分布交錯復(fù)雜，井間的橫向?qū)Ρ纫约皽y井資料等很難確定地層的連通性，對比的正確性也很難保證；壓力測試、示蹤劑測試、試井在實(shí)施時會影響油田生產(chǎn)的正常進(jìn)行，且花費(fèi)昂貴；數(shù)值模擬需要掌握油層的靜態(tài)資料和動態(tài)資料，如油層的孔、滲、飽、厚參數(shù)場，不同時期油井的地層壓力、流動壓力、產(chǎn)水量、產(chǎn)油量和注水量以及油層的潤濕性、油水相對滲透率等，這些參數(shù)要準(zhǔn)備齊全且要真正符合油藏的實(shí)際是非常困難的。因此，提供一種新的方法來研究油藏連通性問題已經(jīng)成為本領(lǐng)域亟需解決的技術(shù)問題。技術(shù)實(shí)現(xiàn)要素：為了解決上述的缺點(diǎn)和不足，本發(fā)明的目的在于提供一種識別井間連通性及油水動態(tài)預(yù)測的方法。為達(dá)到上述目的，本發(fā)明提供一種識別井間連通性及油水動態(tài)預(yù)測的方法，其包括以下步驟：(1)、根據(jù)物質(zhì)平衡方程建立考慮油水兩相的水驅(qū)油藏井間連通性評價模型；(2)、根據(jù)物質(zhì)平衡原理建立水驅(qū)油藏油水兩相的含油飽和度追蹤方程；(3)、通過Koval方法計算單井的有效控制體積；(4)、給定連通系數(shù)λji、干擾系數(shù)δji及時間常數(shù)τj的初值，分別為：λji0、δji0、τj0；再根據(jù)步驟(1)及步驟(2)中的連通性評價模型、含油飽和度追蹤方程，求得第j口油井第一個時間步長的平均含油飽和度及產(chǎn)液量qj(1)；再求解油水前緣的含油飽和度(5)、根據(jù)油水井的注入和產(chǎn)出數(shù)據(jù)，建立最小化目標(biāo)函數(shù)，并采用約束的最優(yōu)化算法對所述目標(biāo)函數(shù)進(jìn)行最小化計算，得到最優(yōu)的λjik、δjik、(6)、根據(jù)步驟(5)中得到的最優(yōu)的λjik判別注采井之間的連通關(guān)系，即：λjik越大，則該井組連通性越好；(7)、將步驟(5)中得到的最優(yōu)的λjik、δjik、帶入步驟(1)-(3)中所述的連通性評價模型、含油飽和度追蹤方程、Koval方法中分別得到更新后的模型，再根據(jù)更新后的模型進(jìn)行動態(tài)指標(biāo)的預(yù)測。根據(jù)本發(fā)明所述的方法，優(yōu)選地，步驟(1)中所述物質(zhì)平衡方程如下式(1)所示：式(1)中，：為綜合壓縮系數(shù)，MPa-1；為生產(chǎn)井j的泄油控制體積，m3；為生產(chǎn)井j泄油區(qū)內(nèi)的平均壓力，MPa；ui(t)為注入井i在第t時間步長的注入量，m3/d；qk(t)為生產(chǎn)井k在第t時間步長的產(chǎn)液量，m3/d；I為注入井?dāng)?shù)，口；P為生產(chǎn)井?dāng)?shù)，口；λji為注入量ui(t)的權(quán)重系數(shù)，即注入井i和生產(chǎn)井j之間的連通系數(shù)；δjk為生產(chǎn)井間的干擾系數(shù)。根據(jù)本發(fā)明所述的方法，優(yōu)選地，步驟(1)中所述水驅(qū)油藏井間連通性評價模型如下式(9)所示：式(9)中，n為某一生產(chǎn)時間點(diǎn)；n0為初始生產(chǎn)時刻；Δn為生產(chǎn)時間的步長；q(n)為某一生產(chǎn)時間點(diǎn)油井的產(chǎn)液量向量，m3/d；Nc為生產(chǎn)井間的干擾系數(shù)向量；Mc為連通系數(shù)向量；τ為時間常數(shù)向量；u(m)為某一時間點(diǎn)注入井注入量向量，m3/d。根據(jù)本發(fā)明所述的方法，優(yōu)選地，步驟(2)中所述物質(zhì)平衡方程如下式(10)所示：式(10)中，Vb為單井控制體積，m3；φ為孔隙度；ρo為原油密度，kg/m3；為平均含油飽和度。根據(jù)本發(fā)明所述的方法，優(yōu)選地，步驟(2)中所述含油飽和度追蹤方程如下式(11)所示：式(11)中，k為某一時間步，Δt為時間步長，Cf為巖石壓縮系數(shù)，1/MPa；Co為原油壓縮系數(shù)，1/MPa。根據(jù)本發(fā)明所述的方法，優(yōu)選地，步驟(3)中Koval方法所用的關(guān)系式如下所示：其中，式中：為井底含水率；kval為Koval系數(shù)，即地層非均質(zhì)性系數(shù)；tD為累計注入水孔隙體積倍數(shù)；Vpj為單井的有效控制體積，m3。根據(jù)本發(fā)明所述的方法，優(yōu)選地，步驟(3)具體包括以下步驟：分別給定Koval系數(shù)及單井的有效控制體積的初值kval0、Vpj0，擬合實(shí)際的單井含水率數(shù)據(jù)，得到最優(yōu)的單井有效控制體積Vpjk。根據(jù)本發(fā)明所述的方法，優(yōu)選地，步驟(4)中根據(jù)油水飽和度追蹤方程求解油水前緣的含油飽和度其中，步驟(4)中所述的油水飽和度追蹤方程(B-L方程)為本領(lǐng)域公知的方程，詳見《油氣層滲流力學(xué)》，張建國等編著，中國石油大學(xué)出版社，2010年，P221。根據(jù)本發(fā)明所述的方法，優(yōu)選地，步驟(4)中根據(jù)以下公式給定連通系數(shù)λji、干擾系數(shù)δji及時間常數(shù)τj的初值；式中：為第i口水井與第j口油井間的平均滲透率，10-3μm2；為第i口水井與第j口油井間的滲流截面積，m2；Lij為第i口水井與第j口油井間的井距，m；φ為孔隙度；η為導(dǎo)壓系數(shù)，μm2/s，根據(jù)本發(fā)明所述的方法，優(yōu)選地，步驟(5)中所述目標(biāo)函數(shù)如下式12所示：約束條件如下：式(12)中：為油井j在第k時間步產(chǎn)液量的模型計算值；為油井j在第k時間步產(chǎn)液量的實(shí)際觀測值；為油井j在第k時間步產(chǎn)油量的模型計算值；為油井j在第k時間步產(chǎn)油量的實(shí)際觀測值；nt為總的時間步長。根據(jù)本發(fā)明所述的方法，優(yōu)選地，步驟(5)中所述動態(tài)指標(biāo)包括產(chǎn)液量、產(chǎn)油量、含油飽和度及含水率。本發(fā)明提供一種識別井間連通性及油水動態(tài)預(yù)測的方法，其具體包括以下步驟：步驟1：根據(jù)物質(zhì)平衡方程建立考慮油水兩相的連通性評價模型，其具體包括以下步驟；假設(shè)研究區(qū)域有I口注水井，P口生產(chǎn)井，建立考慮流體壓縮性，忽略毛管力及重力影響的物質(zhì)平衡方程，如下所示：式(1)中，：為綜合壓縮系數(shù)，MPa-1；為生產(chǎn)井j的泄油控制體積，m3；為生產(chǎn)井j泄油區(qū)內(nèi)的平均壓力，MPa；ui(t)為注入井i在第t時間步長的注入量，m3/d；qk(t)為生產(chǎn)井k在第t時間步長的產(chǎn)液量，m3/d；I為注入井?dāng)?shù)，口；P為生產(chǎn)井?dāng)?shù)，口；λji為注入量ui(t)的權(quán)重系數(shù)，即注入井i和生產(chǎn)井j之間的連通系數(shù)；δjk為生產(chǎn)井間的干擾系數(shù)。生產(chǎn)井產(chǎn)液量與生產(chǎn)井井底壓力又存在以下線性關(guān)系：式(2)中：Jj為產(chǎn)液指數(shù)，m3/(MPa·d)；Pj,wf為生產(chǎn)井井底流壓，MPa；將式(2)帶入式(1)，同時忽略井底壓力的變化，可以得到連通性模型的最終形式為：式(3)中：q(t)＝[q1(t)q2(t)...qP(t)]T(4)；u(t)＝[u1(t)u2(t)...uI(t)]T(5)；τ＝[τ'1τ'2...τ'P]T(6)；Mc＝[λj1,…,λjI](7)；Nc＝[δj1,…,δjp](8)；其中：j＝{1,...,N}定義為時間常數(shù)；為采油指數(shù)，Mj為流體的流度。將式(3)離散化后可得：式(9)中，n為某一生產(chǎn)時間點(diǎn)；n0為初始生產(chǎn)時刻；Δn為生產(chǎn)時間的步長。步驟2：建立油水兩相的飽和度追蹤方程；根據(jù)物質(zhì)平衡原理建立含油飽和度方程如下所示：式(10)中，Vb為單井控制體積，m3；φ為孔隙度；ρo為原油密度，kg/m3；為平均含油飽和度；Cf為巖石壓縮系數(shù)，1/MPa；Co為原油壓縮系數(shù)，1/MPa。通過半解析法求解上述微分方程可得：式(11)中，k為某一時間步；Δt為時間步長。步驟3：通過Koval方法計算單井的有效控制體積給出Koval方法的關(guān)系式如下：其中，式中：為井底含水率；kval為Koval系數(shù)，即地層非均質(zhì)性系數(shù)；tD為累計注入水孔隙體積倍數(shù)；Vpj為單井的有效控制體積，m3。分別給定Koval系數(shù)及單井的有效控制體積的初值kval0、Vpj0，應(yīng)用MATLAB優(yōu)化工具箱等最優(yōu)化軟件，擬合實(shí)際的單井含水率數(shù)據(jù)，得到最優(yōu)的單井有效控制體積Vpjk。步驟4：給定初值，進(jìn)行產(chǎn)液量及含油飽和度初算，其包括以下步驟：近似認(rèn)為同一井組內(nèi)油井流壓相同，則：式中：為第i口水井與第j口油井間的平均滲透率，10-3μm2；為第i口水井與第j口油井間的滲流截面積，m2；Lij為第i口水井與第j口油井間的井距，m；φ為孔隙度；η為導(dǎo)壓系數(shù)，μm2/s，根據(jù)油田現(xiàn)場地質(zhì)資料，由上述公式初步估計連通系數(shù)λji、干擾系數(shù)δji及時間常數(shù)τj的初值，分別為：λji0、δji0、τj0；再根據(jù)步驟(1)及步驟(2)中的連通性評價模型、含油飽和度追蹤方程，求得第j口油井第一個時間步長的平均含油飽和度及產(chǎn)液量qj(1)；再根據(jù)油水飽和度追蹤方程(B-L方程)進(jìn)而求解油水前緣的含油飽和度步驟5：建立目標(biāo)函數(shù)，進(jìn)行目標(biāo)函數(shù)的最小化處理根據(jù)油水井的產(chǎn)出動態(tài)特征，建立最小化目標(biāo)函數(shù)如下：約束條件如下：式(12)中：為油井j在第k時間步產(chǎn)液量的模型計算值；為油井j在第k時間步產(chǎn)液量的實(shí)際觀測值；為油井j在第k時間步產(chǎn)油量的模型計算值；為油井j在第k時間步產(chǎn)油量的實(shí)際觀測值；nt為總的時間步長。應(yīng)用MATLAB最優(yōu)化工具箱等最優(yōu)化軟件，采用約束的最優(yōu)化算法對(12)式進(jìn)行最小化計算，得到最優(yōu)的λjik、δjik、步驟6：井間連通性判別及油水動態(tài)指標(biāo)計算通過最優(yōu)化得到的連通系數(shù)λji，即可判別注采井之間的連通關(guān)系，連通系數(shù)越大，則該井組連通性越好。更新后的模型可以進(jìn)行產(chǎn)液量、產(chǎn)油量、含油飽和度、含水率等動態(tài)指標(biāo)的預(yù)測。油藏可看成一個水動力學(xué)平衡的系統(tǒng)，而注采井則可以看作信號的發(fā)射器和接收器，注入井注入?yún)?shù)的變化就可以引起生產(chǎn)井相關(guān)參數(shù)的改變，生產(chǎn)井相關(guān)參數(shù)的變化受注入井注入?yún)?shù)、地層屬性及其他井的影響。油田現(xiàn)場一般根據(jù)此方法定性判斷井間連通性，但是存在主觀性強(qiáng)、考慮因素不夠多、結(jié)果不夠準(zhǔn)確的問題。許多學(xué)者也利用注采數(shù)據(jù)建立了多種反演模型，如多元線性回歸模型、電容模型、神經(jīng)網(wǎng)絡(luò)、小波分析、流線模型等，但是這些模型目前仍存在以下問題：模型太理想化，反演參數(shù)無明確地質(zhì)意義，不能分層進(jìn)行連通性計算，只能對產(chǎn)液數(shù)據(jù)進(jìn)行預(yù)測擬合，難以表征層間非均質(zhì)性及反演結(jié)果的多解性等問題。因此就需要我們對基于注采數(shù)據(jù)的井間連通性模型進(jìn)行更加準(zhǔn)確的表征。本發(fā)明主要將注入井、生產(chǎn)井及井間孔道看成一個整體，結(jié)合油田靜動態(tài)數(shù)據(jù)，建立水驅(qū)油藏井間連通性定量識別模型，并對模型進(jìn)行求解，進(jìn)而識別井間動態(tài)連通效果。同時將油水兩相的影響因素考慮到模型當(dāng)中，進(jìn)行含油飽和度、含水率、控制體積等相關(guān)指標(biāo)的計算。儲層建模和其他模擬等方法均需要大量的地質(zhì)及地球物理數(shù)據(jù)，評價結(jié)果也不夠準(zhǔn)確，而示蹤劑技術(shù)雖然比較準(zhǔn)確，但需耗費(fèi)大量時間和財力。而本發(fā)明采用基于注采數(shù)據(jù)的反演模型評價井間連通性卻能很好解決這些問題，該方法僅需井的注采數(shù)據(jù)，花費(fèi)少量時間便可以對井間連通性進(jìn)行準(zhǔn)確的評價，同時也可以計算得到相關(guān)油水動態(tài)指標(biāo)。本發(fā)明通過采用油田中最容易獲取的注采動態(tài)數(shù)據(jù)，建立連通性反演模型及飽和度追蹤模型進(jìn)行井組間的連通性評價及油水動態(tài)指標(biāo)的計算。該模型將油藏系統(tǒng)離散成一系列由井間連通系數(shù)、時間常數(shù)及生產(chǎn)井間干擾系數(shù)等參數(shù)表征的井間連通單元，以連通單元為對象建立物質(zhì)平衡方程。以此為基礎(chǔ)，利用MATLAB優(yōu)化工具箱等最優(yōu)化軟件進(jìn)行優(yōu)化，通過對實(shí)際動態(tài)注采參數(shù)擬合，反演模型參數(shù)，最后得到井間連通系數(shù)。同時結(jié)合Koval方法及飽和度追蹤方程進(jìn)行油水兩相的連通性計算及動態(tài)指標(biāo)預(yù)測。結(jié)合含有封閉斷層油藏模型的實(shí)例應(yīng)用，本發(fā)明取得了很好的動態(tài)擬合和預(yù)測效果，反演參數(shù)和實(shí)際油藏具有良好的一致性，驗(yàn)證了該模型的可靠性。相對當(dāng)前的連通性評價及油水動態(tài)指標(biāo)計算方法，本發(fā)明計算速度快，成本低周期短，可同時進(jìn)行多口井的連通性反演及動態(tài)指標(biāo)計算，對指導(dǎo)油田生產(chǎn)意義極大。本發(fā)明所提供的方法可以很好地進(jìn)行井間連通性的定量識別及油水動態(tài)指標(biāo)的準(zhǔn)確預(yù)測。相對傳統(tǒng)的連通性評價方法及經(jīng)驗(yàn)公式計算油水動態(tài)指標(biāo)，該方法所需數(shù)據(jù)量少且容易獲取，不需要改變井的工作制度且可以進(jìn)行大規(guī)模油田的計算，在當(dāng)前低油價形勢下，可以極大地減少成本，縮短周期，有效指導(dǎo)油田的生產(chǎn)運(yùn)營。附圖說明圖1為本發(fā)明實(shí)施例中含有封閉斷層的油藏平面滲透率場圖；圖2a-圖2d分別為本發(fā)明實(shí)施例中含有封閉斷層油藏P1-P4四口生產(chǎn)井產(chǎn)液量擬合圖；圖3為實(shí)施例中含有封閉斷層油藏P1-P4四口生產(chǎn)井的井時間常數(shù)反演圖；圖4為實(shí)施例中含有封閉斷層油藏P1-P4四口生產(chǎn)井的平均含油飽和度變化圖；圖5為實(shí)施例中含有封閉斷層油藏P1-P4四口生產(chǎn)井的點(diǎn)含油飽和度變化圖；圖6為本發(fā)明實(shí)施例中含有封閉斷層油藏井間連通性圖。具體實(shí)施方式為了對本發(fā)明的技術(shù)特征、目的和有益效果有更加清楚的理解，現(xiàn)結(jié)合以下具體實(shí)施例及說明書附圖對本發(fā)明的技術(shù)方案進(jìn)行以下詳細(xì)說明。顯而易見，下面描述中的附圖僅僅是本發(fā)明的具體實(shí)施例，所描述的實(shí)施例僅僅是本發(fā)明一部分實(shí)施例，而不是全部實(shí)施例?；诒景l(fā)明中的實(shí)施例，本領(lǐng)域普通技術(shù)人員在沒有做出創(chuàng)造性勞動前提下所獲得的所有其他實(shí)施例，都屬于本發(fā)明保護(hù)范圍。實(shí)施例1本實(shí)施例提供了一種識別井間連通性及油水動態(tài)預(yù)測的方法，其包括以下步驟：1)具體實(shí)施方法第一步：通過所建立的連通性評價模型，擬合實(shí)際注采動態(tài)數(shù)據(jù)，進(jìn)行連通系數(shù)等相關(guān)參數(shù)的反演。根據(jù)連通系數(shù)的大小，結(jié)合實(shí)際油藏動靜態(tài)資料分析井組的連通狀況。第二步：根據(jù)Koval方法及飽和度追蹤方程進(jìn)行油水兩相動態(tài)指標(biāo)計算，結(jié)合連通性反演模型對含水率、含油飽和度、單井有效控制體積等有效參數(shù)進(jìn)行計算及預(yù)測。2)具體實(shí)施例借助Eclipse軟件，應(yīng)用本發(fā)明所提供的連通性模型對含有斷層的典型油藏模型進(jìn)行井間動態(tài)連通性反演。所建Eclipse模型如圖1所示，該模型共有21×21×1＝441個網(wǎng)格，網(wǎng)格長度為20m，平均有效厚度為5m。其中，黑色部分滲透率為200mD，白色部分滲透率為0，孔隙度為0.25。采用五點(diǎn)法井網(wǎng)(5注4采)，油井定壓生產(chǎn)。5口注水井的月動態(tài)數(shù)據(jù)如表1所示。表1(1)給出連通系數(shù)λji、干擾系數(shù)δji及時間常數(shù)τj初值，如表2所示：表2λjiI1(i＝1)I2(i＝2)I3(i＝3)I4(i＝4)I5(i＝5)τjδjiP1(j＝1)0.250.250.250.250.250.200P2(j＝2)0.250.250.250.250.250.200P3(j＝3)0.250.250.250.250.250.200P4(j＝4)0.250.250.250.250.250.200同時將表1及表2中的數(shù)據(jù)帶入如下式(9)中，進(jìn)而可以計算出每個月的產(chǎn)液量：qj(n)，如表3所示。表3(2)應(yīng)用MATLAB優(yōu)化工具箱，擬合模型計算的產(chǎn)液量和實(shí)際的單井產(chǎn)液量數(shù)據(jù)，擬合效果圖如圖2a-圖2d所示，得到優(yōu)化后的λji、δji、τj，如表4、5所示：表4λjiI1(i＝1)I2(i＝2)I3(i＝3)I4(i＝4)I5(i＝5)τP1(j＝1)00.490000.27P2(j＝2)0.6800.490.490.300.39P3(j＝3)00.480000.15P4(j＝4)0.2900.480.480.670.30表5δjiP1(p＝1)P2(p＝2)P3(p＝3)P4(p＝4)P1(j＝1)0.030-0.030P2(j＝2)00.040-0.04P3(j＝3)-0.0300.040P4(j＝4)0-0.0400.05(3)給定平均含油飽和度初值結(jié)合上一步計算的模型參數(shù)、Koval方法得出的單井有效控制體積Vpj＝33235m3；及式(11)，求解出每一時間步的平均含油飽和度如表6所示，并通過B-L方程計算油水前緣的含油飽和度如表7所示。表6表7通過相滲關(guān)系獲取相應(yīng)的油水相對滲透率，這里采用相關(guān)經(jīng)驗(yàn)公式：式中：krw為水的相對滲透率，mD；kro為油的相對滲透率，mD；Sw為含水飽和度；Swr為束縛水飽和度；Sor為殘余油飽和度。進(jìn)而求得每一步長的流度如表8所示：式中μw＝0.72，為水的粘度，mPa·s；μo＝1.63，為油的粘度，mPa·s。表8求取更新后的時間常數(shù)τj'＝τj/Mj，同時帶入到式(9)當(dāng)中。(4)將上述步驟求解的產(chǎn)液量數(shù)據(jù)(如表9)，Koval方法得到的含水率進(jìn)而求取產(chǎn)油量數(shù)據(jù)(如表10)，帶入式(12)中；約束條件如下：表9表10得到最優(yōu)化的λji'、、τj”如表11所示，δji'如表12所示。表11λjiI1(i＝1)I2(i＝2)I3(i＝3)I4(i＝4)I5(i＝5)τjP1(j＝1)00.490000.20P2(j＝2)0.6800.490.490.300.35P3(j＝3)00.490000.20P4(j＝4)0.3000.490.490.680.35表12δjiP1(p＝1)P2(p＝2)P3(p＝3)P4(p＝4)P1(j＝1)0.050-0.040P2(j＝2)00.040-0.04P3(j＝3)-0.0400.050P4(j＝4)0-0.0400.05優(yōu)化后時間常數(shù)的變化如圖3所示，平均含油飽和度及井點(diǎn)含油飽和度變化分別如圖4及圖5所示，最終連通性大小采用如圖6中的箭頭表示，箭頭越大則連通性越好，否則連通性越小。(5)通過上述優(yōu)化后的參數(shù)代入模型中可進(jìn)行產(chǎn)液量、產(chǎn)油量、前緣含油飽和度等動態(tài)指標(biāo)的計算，計算結(jié)果分別如表13-15所示。表13表14表15當(dāng)前第1頁1 2 3