專利名稱：油藏井間動(dòng)態(tài)連通性反演方法
技術(shù)領(lǐng)域：
本發(fā)明涉及一種油藏連通性的評(píng)價(jià)方法，尤其是采用油藏注入和采出動(dòng)態(tài)數(shù)據(jù)反
演得到油藏井間動(dòng)態(tài)連通性的方法。
背景技術(shù)：
在油田開(kāi)發(fā)中，油藏連通性是油藏評(píng)價(jià)的重要內(nèi)容，也是影響油田開(kāi)發(fā)動(dòng)態(tài)的重 要參數(shù)。油藏連通性的認(rèn)識(shí)對(duì)提高油藏原油采收率具有重要意義。隨著油藏的不斷開(kāi)發(fā)， 油藏參數(shù)已經(jīng)發(fā)生了很大變化，靜態(tài)連通性已經(jīng)不能準(zhǔn)確反映地層性質(zhì)和儲(chǔ)層流體連通狀 況。了解井間動(dòng)態(tài)連通性、弄清注水井水流方向，不僅能夠指導(dǎo)高含水期調(diào)剖堵水等工藝措 施的優(yōu)化實(shí)施，而且對(duì)于剩余油分布的定量描述具有重要意義。 目前，在油田主要通過(guò)示蹤劑測(cè)試、試井、油藏?cái)?shù)值模擬等方法來(lái)研究油藏井間動(dòng) 態(tài)連通性。但示蹤劑測(cè)試、試井在實(shí)施時(shí)會(huì)影響油田的正常生產(chǎn)，且費(fèi)用昂貴；數(shù)值模擬方 法需要掌握油層的靜態(tài)資料和動(dòng)態(tài)資料，如油層的孔隙度、滲透率、流體飽和度、油層厚度 等參數(shù)場(chǎng)，不同時(shí)刻油水井的注采數(shù)據(jù)等，這些參數(shù)要準(zhǔn)備齊全且要真正符合油藏實(shí)際是 比較困難的。

發(fā)明內(nèi)容
本發(fā)明的目的是為了克服常規(guī)油藏井間連通性研究手段成本高、資料準(zhǔn)備困難的 缺點(diǎn)，提供一種油藏井間動(dòng)態(tài)連通性反演方法，利用油藏注入和采出動(dòng)態(tài)數(shù)據(jù)反演得到油 藏井間動(dòng)態(tài)連通性。 本發(fā)明采取的技術(shù)路線是首先建立油藏井間動(dòng)態(tài)連通性反演的多元回歸模型， 根據(jù)注入量信號(hào)在油藏中的傳播特征構(gòu)建一階線性系統(tǒng)褶積器，應(yīng)用褶積器對(duì)注入量信號(hào) 進(jìn)行預(yù)處理，然后將遺傳算法引入反演模型的優(yōu)化求解中，通過(guò)對(duì)待優(yōu)化參數(shù)的系統(tǒng)訓(xùn)練， 反演得到表征油藏井間動(dòng)態(tài)連通性的連通系數(shù)以及表征注入信號(hào)時(shí)滯性和衰減性的時(shí)間 常數(shù)。 由于油藏中的注水井、采油井以及井間孔道是一個(gè)完整的系統(tǒng)，注水井注入量是 該系統(tǒng)的輸入信號(hào)，油井產(chǎn)液量則是該系統(tǒng)的輸出信號(hào)。注入量信號(hào)在油藏介質(zhì)中擴(kuò)散傳 播時(shí)存在一定的衰減、時(shí)滯，如附圖1所示。 圖2所示為系統(tǒng)單位階躍響應(yīng)，q(①)是響應(yīng)的終值，與激勵(lì)值接近，0.5q(①) 則是響應(yīng)曲線終值的一半。其動(dòng)態(tài)性能指標(biāo)中的延遲時(shí)間Td定義為響應(yīng)曲線第一次達(dá)到 其終值一半所需的時(shí)間。基于階躍響應(yīng)動(dòng)態(tài)性能指標(biāo)，定義了表征注入信號(hào)時(shí)滯性和衰減 性的指標(biāo)t，定義為 T = i (1)
y 式中，t為表征注入信號(hào)時(shí)滯性和衰減性的時(shí)間常數(shù)，month ;Td為延遲時(shí)間， month ; y為延遲時(shí)間與時(shí)間常數(shù)的折算系數(shù)，其經(jīng)驗(yàn)取值為3. 0。
結(jié)合油藏?cái)?shù)值模擬進(jìn)行注入信號(hào)時(shí)滯性和衰減性影響因素研究，得到了時(shí)間常數(shù)
與導(dǎo)壓系數(shù)、井距以及孔隙度的關(guān)系式
<formula>formula see original document page 5</formula>
式中，小為孔隙度，小數(shù)；d為井距，m;化為導(dǎo)壓系數(shù)，cm7s。 在油田生產(chǎn)中，每一口井的產(chǎn)量變化都與周圍與之連通的所有注水井的共同作用 相關(guān)聯(lián)。根據(jù)多元線性回歸思想，對(duì)于一個(gè)由N工口注水井和Np 口生產(chǎn)井組成的注采系統(tǒng)， 生產(chǎn)井j的產(chǎn)液量可以由相鄰的注入井注入量表示為
<formula>formula see original document page 5</formula>
式中，Q.為第j 口生產(chǎn)井產(chǎn)液量估計(jì)值，mVd ; 13 。為表征注采不平衡的常數(shù)項(xiàng)，mV d ; 13 p產(chǎn)液量初值影響的權(quán)重系數(shù)；Qj(k。)為第j 口生產(chǎn)井產(chǎn)液量初值，m3/d ; t p為產(chǎn)液量 初值影響的時(shí)間常數(shù)，month ;k為采樣時(shí)間，month ;k。為初始采樣時(shí)間，month ;Ii(m)為第 m月第i 口注水井的注水量，m3/d ; t ij為對(duì)應(yīng)于第j 口生產(chǎn)井和第i 口注水井間的時(shí)間常 數(shù)，month ; 13 ij為第j 口生產(chǎn)井和第i 口注水井間的多元線性回歸權(quán)重，定義為井間動(dòng)態(tài)連 通系數(shù)。該式右端項(xiàng)包含三部分第一部分為表征注采不平衡的常數(shù)項(xiàng)；第二部分為產(chǎn)液 量初始值的影響；第三部分為注入信號(hào)預(yù)處理后修正值。 一般情況下，產(chǎn)液量初始值的影響 較小。當(dāng)油藏注采平衡且產(chǎn)液量初值為零時(shí)，模型右端項(xiàng)僅含第三部分。
本發(fā)明具體按以下步驟進(jìn)行 (1)選取注采動(dòng)態(tài)數(shù)據(jù)用于油藏井間動(dòng)態(tài)連通性反演的注入井的月度注水量、
生產(chǎn)井的月度產(chǎn)液量數(shù)據(jù)可從油田月度生產(chǎn)動(dòng)態(tài)數(shù)據(jù)報(bào)表中直接獲取。為了具有更好的反
演效果，注采動(dòng)態(tài)數(shù)據(jù)的選取一般盡可能考慮動(dòng)態(tài)數(shù)據(jù)有一定的波動(dòng)，波動(dòng)幅度一般大于
±20%，各單井注采量保持一定的連續(xù)性，選取數(shù)據(jù)時(shí)間段一般應(yīng)大于1年。
(2)生成初始群體根據(jù)時(shí)間常數(shù)與導(dǎo)壓系數(shù)、井距以及孔隙度的回歸關(guān)系式<formula>formula see original document page 5</formula>
，估計(jì)出時(shí)間常數(shù)初始值。以此為基礎(chǔ)，產(chǎn)生一組實(shí)數(shù)型的時(shí)間常
數(shù)值構(gòu)成一個(gè)個(gè)體，個(gè)體實(shí)數(shù)個(gè)數(shù)即為時(shí)間常數(shù)待定參數(shù)個(gè)數(shù)。然后在時(shí)間常數(shù)初始值的 50 % 150%區(qū)間內(nèi)，按均勻分布隨機(jī)產(chǎn)生&個(gè)這種個(gè)體構(gòu)成初始群體。
(3)求解井間動(dòng)態(tài)連通系數(shù)結(jié)合時(shí)間常數(shù)群體數(shù)據(jù)應(yīng)用褶積器對(duì)注入量信號(hào)進(jìn) 行預(yù)處理，預(yù)處理的目的是引入時(shí)間常數(shù)修正注入信號(hào)的時(shí)滯性和衰減性對(duì)油藏井間動(dòng)態(tài) 連通性反演結(jié)果的影響。得到修正后的注入數(shù)據(jù)后，應(yīng)用最小二乘法求得表征井間動(dòng)態(tài)連 通程度的多元線性回歸的權(quán)重系數(shù)，權(quán)重系數(shù)的最小二乘估計(jì)滿足下式
<formula>formula see original document page 5</formula>
式中，Qj為第j 口生產(chǎn)井產(chǎn)液量實(shí)際值，m3/d ;Nk為數(shù)據(jù)點(diǎn)總數(shù)。
(4)適應(yīng)度函數(shù)定義對(duì)基于信號(hào)處理的井間動(dòng)態(tài)連通性反演模型的求解可看作
連續(xù)參數(shù)優(yōu)化問(wèn)題，即優(yōu)化時(shí)間常數(shù)使模型計(jì)算的產(chǎn)液量曲線與實(shí)際觀察的產(chǎn)液量曲線最為接近，達(dá)到minFj。目標(biāo)函數(shù)Fj可定義為巧(、，、，、,......, )=￡ (込W - Q W) / Q (y. = i,…,W) (5) 該目標(biāo)函數(shù)值反映的是第j 口生產(chǎn)井產(chǎn)液量估計(jì)值與實(shí)際值之間的累積相對(duì)誤 差。以產(chǎn)液量估計(jì)值與實(shí)際值的相對(duì)誤差定義適應(yīng)度函數(shù)，應(yīng)保證產(chǎn)液量估計(jì)值與實(shí)際值 之間的誤差越大，適應(yīng)度越小，取適應(yīng)度函數(shù)fj = 1/Fj。
(5)適應(yīng)性評(píng)估以適應(yīng)度函數(shù)來(lái)確定個(gè)體適應(yīng)度值以評(píng)估其優(yōu)劣，并以此作為
遺傳操作的依據(jù)。計(jì)算時(shí)間常數(shù)群體中各個(gè)體的適應(yīng)度&(1 = 1，…，N》。
(6)遺傳變異操作操作的目的是根據(jù)進(jìn)化原則從當(dāng)前群體中選出優(yōu)良的個(gè)體，
從而形成下一代。首先進(jìn)行停止規(guī)則判斷，若發(fā)現(xiàn)占群體80% 90%以上的個(gè)體在相對(duì)誤
差為設(shè)定值范圍內(nèi)是同一個(gè)體，或者算法迭代步數(shù)超出設(shè)定的閾值，則終止算法迭代，當(dāng)前
最優(yōu)個(gè)體組成即為所求時(shí)間常數(shù)的待定值，進(jìn)而再重新計(jì)算井間動(dòng)態(tài)連通系數(shù)，輸出結(jié)果；
否則進(jìn)行選擇、交叉和變異等操作，產(chǎn)生子代，重復(fù)步驟(3) (6)，繼續(xù)進(jìn)化。
(7)輸出反演結(jié)果將得到的井間動(dòng)態(tài)連通系數(shù)繪制動(dòng)態(tài)連通圖。
本發(fā)明的有益效果是通過(guò)油藏注入和采出動(dòng)態(tài)數(shù)據(jù)反演得到油藏井間動(dòng)態(tài)連通
性，引入了信號(hào)處理技術(shù)，消除了注入信號(hào)時(shí)滯性和衰減性對(duì)井間動(dòng)態(tài)連通性反演效果的
影響，實(shí)現(xiàn)了僅利用易于獲取的注采數(shù)據(jù)就可得到較準(zhǔn)確的動(dòng)態(tài)連通性反演結(jié)果。具有以
下優(yōu)點(diǎn)①反演所需要的資料較少，僅限于注采動(dòng)態(tài)數(shù)據(jù)，容易準(zhǔn)備齊全；②生產(chǎn)動(dòng)態(tài)數(shù)據(jù)
資料容易獲取，并且不影響油田的正常生產(chǎn)；③所需資料獲取不用額外增加測(cè)量和測(cè)試儀
器，反演過(guò)程簡(jiǎn)單便捷。因此，利用開(kāi)發(fā)動(dòng)態(tài)數(shù)據(jù)來(lái)反演油藏井間動(dòng)態(tài)連通性能得到較準(zhǔn)確
有效的結(jié)果，為進(jìn)一步指導(dǎo)油田開(kāi)發(fā)方案的編制和生產(chǎn)優(yōu)化提供了理論依據(jù)。


圖1為注入信號(hào)的時(shí)滯性與衰減性示意圖。
圖2為系統(tǒng)的階躍響應(yīng)圖。 圖3為依據(jù)本發(fā)明所提出的油藏井間動(dòng)態(tài)連通性反演求解計(jì)算機(jī)程序框圖。
圖4為均質(zhì)概念模型井間動(dòng)態(tài)連通圖。
圖5為陸9井區(qū)J2x4油藏井間動(dòng)態(tài)連通圖。
具體實(shí)施例方式
下面結(jié)合實(shí)施例和附圖來(lái)進(jìn)一步描述本發(fā)明。 基于信號(hào)處理的油藏井間動(dòng)態(tài)連通性反演求解步驟如附圖3所示，以本發(fā)明在陸 9井區(qū)J2x4油藏的現(xiàn)場(chǎng)應(yīng)用為例，來(lái)進(jìn)一步說(shuō)明本發(fā)明。 陸梁油田位于準(zhǔn)格爾盆地腹部古爾班通古特沙漠腹地，陸9井區(qū)位于陸梁隆起三 個(gè)泉凸起西北部陸9井一帶點(diǎn)，工區(qū)內(nèi)共發(fā)育10多條斷層，集中分布于陸9井區(qū)背斜構(gòu)造 的東傾末端的南北兩翼。陸9井區(qū)西山窯組J2x4油藏為邊水活躍的低飽和、底水塊狀、低 幅度背斜砂巖油藏。原油密度低、粘度低，含蠟量高，凝固點(diǎn)高。該區(qū)地層壓力為20.6MPa， 地層溫度為65. 6t:，飽和壓力為4. 19MPa，屬于常溫常壓系統(tǒng)。該油藏于2000年10月投入 開(kāi)發(fā)，2001年11月投入注水開(kāi)發(fā)。截至目前，陸9井區(qū)J2x4油藏共投產(chǎn)油井78 口 ，水井30口。油井日產(chǎn)液水平1158. 9t/d，綜合含水59. 08%，水井日注水平1356m 月注采比1. 10。
采用本發(fā)明對(duì)陸9井區(qū)J2x4油藏進(jìn)行井間動(dòng)態(tài)連通性研究，具體的實(shí)施過(guò)程為
(1)注采動(dòng)態(tài)數(shù)據(jù)的選取 選取2001年12月 2007年7月為評(píng)價(jià)時(shí)間段，從油田月度生產(chǎn)動(dòng)態(tài)數(shù)據(jù)報(bào)表中 直接獲取得到這一時(shí)間段內(nèi)涉及的30 口注入井的月度注水量，以及78 口生產(chǎn)井的月度產(chǎn) 液量數(shù)據(jù)。整體來(lái)看，動(dòng)態(tài)數(shù)據(jù)有一定的波動(dòng)，各單井產(chǎn)量也保持了一定的連續(xù)性，符合井 間動(dòng)態(tài)連通性反演方法應(yīng)用條件。
(2)初始群體的生成 根據(jù)時(shí)間常數(shù)與導(dǎo)壓系數(shù)、井距以及孔隙度的回歸關(guān)系式(即式(2))估計(jì)時(shí)間常
數(shù)初始值。對(duì)于不同井組，導(dǎo)壓系數(shù)、井距以及孔隙度有所不同，估計(jì)的時(shí)間常數(shù)初始值也
不同。產(chǎn)生的一組實(shí)數(shù)型的時(shí)間常數(shù)值構(gòu)成一個(gè)個(gè)體，個(gè)體實(shí)數(shù)個(gè)數(shù)為時(shí)間常數(shù)待定參數(shù)
個(gè)數(shù)。然后在時(shí)間常數(shù)初始值的50% 150%區(qū)間內(nèi)按均勻分布隨機(jī)產(chǎn)生^ = 1000個(gè)這
種個(gè)體構(gòu)成初始群體。 (3)井間動(dòng)態(tài)連通系數(shù)的求解 結(jié)合時(shí)間常數(shù)群體數(shù)據(jù)應(yīng)用褶積器對(duì)注入量信號(hào)進(jìn)行預(yù)處理，得到修正后的注入 數(shù)據(jù)，然后應(yīng)用最小二乘法求得表征井間動(dòng)態(tài)連通程度的多元線性回歸的權(quán)重系數(shù)。權(quán)重 系數(shù)的最小二乘估計(jì)滿足式(4)。
(4)適應(yīng)度函數(shù)定義 對(duì)基于信號(hào)處理的井間動(dòng)態(tài)連通性反演模型的求解可看作是一個(gè)連續(xù)參數(shù)優(yōu)化 問(wèn)題，即優(yōu)化時(shí)間常數(shù)使模型計(jì)算的產(chǎn)液量曲線與實(shí)際觀察的產(chǎn)液量曲線最為接近，達(dá)到 minFj，按照式(5)定義目標(biāo)函數(shù)Fj，并取適應(yīng)度函數(shù)f j = 1/Fj。
(5)適應(yīng)性評(píng)估 以適應(yīng)度函數(shù)來(lái)確定個(gè)體適應(yīng)度值以評(píng)估其優(yōu)劣，并以此作為遺傳操作的依據(jù)。 計(jì)算時(shí)間常數(shù)群體中各個(gè)體的適應(yīng)度= 1，…，N》。
(6)遺傳變異操作 操作的目的是根據(jù)進(jìn)化原則從當(dāng)前群體中選出優(yōu)良的個(gè)體，從而形成下一代。首 先進(jìn)行停止規(guī)則判斷，若發(fā)現(xiàn)占群體90%以上的個(gè)體在相對(duì)誤差為1%范圍內(nèi)是同一個(gè) 體，或者算法迭代步數(shù)超出300計(jì)算時(shí)步時(shí)，終止算法迭代，則當(dāng)前最優(yōu)個(gè)體組成即為所求 時(shí)間常數(shù)的待定值，進(jìn)而再重新計(jì)算井間動(dòng)態(tài)連通系數(shù)，輸出結(jié)果；否則進(jìn)行選擇、交叉和 變異等操作，產(chǎn)生子代，重復(fù)步驟(3) (6)，繼續(xù)進(jìn)化。
(7)輸出反演結(jié)果 將得到的井間動(dòng)態(tài)連通系數(shù)繪制動(dòng)態(tài)連通圖，如附圖4所示，圖中箭頭大小反映
生產(chǎn)井和注入井的連通關(guān)系，箭頭越大，表示兩口井間的連通程度就越好。 通過(guò)上述流程，本發(fā)明能較容易地得到全局最優(yōu)解，最終利用注采動(dòng)態(tài)數(shù)據(jù)簡(jiǎn)捷
地反演出油藏井間動(dòng)態(tài)連通性。 陸9井區(qū)J2x4油藏井間動(dòng)態(tài)連通性計(jì)算結(jié)果如附圖5所示，該圖直觀地反映了各 井組注水井與周圍生產(chǎn)井之間相對(duì)連通性較好的方向，評(píng)價(jià)結(jié)果主要反映的是注入水的影 響?？梢钥闯觯捎趦?chǔ)層非均質(zhì)特征，各井組在不同方向連通情況存在一定差異。整體來(lái)看， 西區(qū)油水主要流動(dòng)方向?yàn)楸睎|_南西方向，中部區(qū)域以南_北方向?yàn)橹鳎瑬|區(qū)以西北_東南方向?yàn)橹?，這與該區(qū)儲(chǔ)層沉積特征基本一致。該區(qū)內(nèi)部斷層均為開(kāi)啟的，因此評(píng)價(jià)結(jié)果顯示 斷層幾乎未對(duì)兩側(cè)注采井之間的動(dòng)態(tài)連通性產(chǎn)生影響，符合油藏井間實(shí)際連通情況。
權(quán)利要求
一種油藏井間動(dòng)態(tài)連通性反演方法，將油藏中的注水井、生產(chǎn)井以及井間孔道看作一個(gè)完整的系統(tǒng)，根據(jù)該系統(tǒng)中生產(chǎn)井的產(chǎn)液量和注水井的注入量動(dòng)態(tài)數(shù)據(jù)反演得到油藏井間動(dòng)態(tài)連通性，其特征是，在注入信號(hào)的時(shí)滯性和衰減性規(guī)律研究的基礎(chǔ)上，應(yīng)用信號(hào)處理技術(shù)，構(gòu)建包含對(duì)注入信號(hào)進(jìn)行預(yù)處理的井間動(dòng)態(tài)連通性反演的多元回歸模型，將遺傳算法引入反演模型的優(yōu)化求解中，反演得到表征油藏井間動(dòng)態(tài)連通性的連通系數(shù)以及表征注入信號(hào)時(shí)滯性和衰減性的時(shí)間常數(shù)，具體按以下步驟實(shí)現(xiàn)(1)注采動(dòng)態(tài)數(shù)據(jù)的選取從油田月度生產(chǎn)動(dòng)態(tài)數(shù)據(jù)報(bào)表中直接獲取注入井的月度注水量、生產(chǎn)井的月度產(chǎn)液量數(shù)據(jù)，為了具有更好的反演效果，注采動(dòng)態(tài)數(shù)據(jù)的選取考慮動(dòng)態(tài)數(shù)據(jù)要有幅度為±20％的波動(dòng)，各單井注采量保持連續(xù)性，選取數(shù)據(jù)時(shí)間段應(yīng)大于1年；(2)生成初始群體根據(jù)時(shí)間常數(shù)與導(dǎo)壓系數(shù)、井距以及孔隙度的回歸關(guān)系式估計(jì)時(shí)間常數(shù)初始值，式中，τ為時(shí)間常數(shù)，month；φ為孔隙度，小數(shù)；d為井距，m；為導(dǎo)壓系數(shù)，cm2/s，以此為基礎(chǔ)，產(chǎn)生一組實(shí)數(shù)型的時(shí)間常數(shù)值構(gòu)成一個(gè)個(gè)體，個(gè)體實(shí)數(shù)個(gè)數(shù)即為時(shí)間常數(shù)待定參數(shù)個(gè)數(shù)，然后在時(shí)間常數(shù)初始值的50％～150％區(qū)間內(nèi)，按均勻分布隨機(jī)產(chǎn)生Nl個(gè)這種個(gè)體構(gòu)成初始群體；(3)求解井間動(dòng)態(tài)連通系數(shù)結(jié)合時(shí)間常數(shù)群體數(shù)據(jù)應(yīng)用褶積器對(duì)注入量信號(hào)進(jìn)行預(yù)處理，得到修正后的注入數(shù)據(jù)，然后應(yīng)用最小二乘法求得表征井間動(dòng)態(tài)連通程度的多元線性回歸的權(quán)重系數(shù)，權(quán)重系數(shù)的最小二乘估計(jì)滿足下式 <mrow><mfrac> <mo>&PartialD;</mo> <mrow><mo>&PartialD;</mo><msub> <mi>&beta;</mi> <mi>ij</mi></msub> </mrow></mfrac><mo>[</mo><munderover> <mi>&Sigma;</mi> <mrow><mi>k</mi><mo>=</mo><mn>1</mn> </mrow> <msub><mi>N</mi><mi>k</mi> </msub></munderover><msup> <mrow><mo>(</mo><msub> <mi>Q</mi> <mi>j</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>-</mo><msub> <mover><mi>Q</mi><mo>^</mo> </mover> <mi>j</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>)</mo> </mrow> <mn>2</mn></msup><mo>]</mo><mo>=</mo><mn>0</mn> </mrow>式中，Qj為第j口生產(chǎn)井產(chǎn)液量實(shí)際值，m3/d；為第j口生產(chǎn)井產(chǎn)液量估計(jì)值，m3/d；βij為第j口生產(chǎn)井和第i口注水井間的多元線性回歸權(quán)重，定義為井間動(dòng)態(tài)連通系數(shù)；Nk為數(shù)據(jù)點(diǎn)總數(shù)；(4)適應(yīng)度函數(shù)定義目標(biāo)函數(shù)Fj可定義為 <mrow><msub> <mi>F</mi> <mi>j</mi></msub><mrow> <mo>(</mo> <msub><mi>&tau;</mi><mi>p</mi> </msub> <mo>,</mo> <msub><mi>&tau;</mi><mrow> <mn>1</mn> <mi>j</mi></mrow> </msub> <mo>,</mo> <msub><mi>&tau;</mi><mrow> <mn>2</mn> <mi>j</mi></mrow> </msub> <mo>,</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>,</mo> <msub><mi>&tau;</mi><mrow> <msub><mi>N</mi><mi>I</mi> </msub> <mi>j</mi></mrow> </msub> <mo>)</mo></mrow><mo>=</mo><munderover> <mi>&Sigma;</mi> <mrow><mi>k</mi><mo>=</mo><mn>1</mn> </mrow> <msub><mi>N</mi><mi>k</mi> </msub></munderover><mo>|</mo><mrow> <mo>(</mo> <msub><mover> <mi>Q</mi> <mo>^</mo></mover><mi>j</mi> </msub> <mrow><mo>(</mo><mi>k</mi><mo>)</mo> </mrow> <mo>-</mo> <msub><mi>Q</mi><mi>j</mi> </msub> <mrow><mo>(</mo><mi>k</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><mo>/</mo><msub> <mi>Q</mi> <mi>j</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>|</mo><mo>,</mo><mrow> <mo>(</mo> <mi>j</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>,</mo> <msub><mi>N</mi><mi>P</mi> </msub> <mo>)</mo></mrow> </mrow>式中，τp為產(chǎn)液量初值影響的時(shí)間常數(shù)，month；τij為對(duì)應(yīng)于第j口生產(chǎn)井和第i口注水井間的時(shí)間常數(shù)，month；NI口注水井井?dāng)?shù)；NP口生產(chǎn)井井?dāng)?shù)。該目標(biāo)函數(shù)值反映的是第j口生產(chǎn)井產(chǎn)液量估計(jì)值與實(shí)際值之間的累積相對(duì)誤差，以產(chǎn)液量估計(jì)值與實(shí)際值的相對(duì)誤差定義適應(yīng)度函數(shù)，應(yīng)保證產(chǎn)液量估計(jì)值與實(shí)際值之間的誤差越大，適應(yīng)度越小，取適應(yīng)度函數(shù)fj＝1/Fj；(5)適應(yīng)性評(píng)估以適應(yīng)度函數(shù)來(lái)確定個(gè)體適應(yīng)度值以評(píng)估其優(yōu)劣，并以此作為遺傳操作的依據(jù)，計(jì)算時(shí)間常數(shù)群體中各個(gè)體的適應(yīng)度f(wàn)jl(l＝1，…，Nl)；(6)遺傳變異操作進(jìn)行遺傳變異操作，操作的目的是根據(jù)進(jìn)化原則從當(dāng)前群體中選出優(yōu)良的個(gè)體，從而形成下一代，首先進(jìn)行停止規(guī)則判斷，若發(fā)現(xiàn)占群體80％～90％以上的個(gè)體在相對(duì)誤差為設(shè)定值范圍內(nèi)是同一個(gè)體，或者算法迭代步數(shù)超出設(shè)定的閾值，就終止算法迭代，則當(dāng)前最優(yōu)個(gè)體組成即為所求時(shí)間常數(shù)的待定值，進(jìn)而再重新計(jì)算井間動(dòng)態(tài)連通系數(shù)，輸出結(jié)果，否則進(jìn)行選擇、交叉和變異操作，產(chǎn)生子代，重復(fù)步驟(3)～(6)，繼續(xù)進(jìn)化；(7)輸出反演結(jié)果將得到的井間動(dòng)態(tài)連通系數(shù)繪制動(dòng)態(tài)連通圖。F2009102565246C00011.tif,F2009102565246C00012.tif,F2009102565246C00014.tif
全文摘要
本發(fā)明公開(kāi)了一種油藏井間動(dòng)態(tài)連通性反演方法，它是基于信號(hào)處理技術(shù)并采用油藏注入和采出動(dòng)態(tài)數(shù)據(jù)反演得到油藏井間動(dòng)態(tài)連通性的方法。首先建立油藏井間動(dòng)態(tài)連通性反演的多元回歸模型，并根據(jù)注入量信號(hào)在油藏中的傳播特征構(gòu)建一階線性系統(tǒng)褶積器，應(yīng)用褶積器對(duì)注入量信號(hào)進(jìn)行了預(yù)處理，然后將遺傳算法引入反演模型的優(yōu)化求解中，通過(guò)對(duì)待優(yōu)化參數(shù)的系統(tǒng)訓(xùn)練，反演得到表征油藏井間動(dòng)態(tài)連通性的連通系數(shù)以及表征注入信號(hào)時(shí)滯性和衰減性的時(shí)間常數(shù)。本發(fā)明實(shí)用性強(qiáng)，通用性高，主要用于評(píng)價(jià)油藏井間動(dòng)態(tài)連通性，進(jìn)一步指導(dǎo)油田開(kāi)發(fā)方案的編制和生產(chǎn)優(yōu)化。
文檔編號(hào)E21B49/00GK101725346SQ20091025652
公開(kāi)日2010年6月9日 申請(qǐng)日期2009年12月15日 優(yōu)先權(quán)日2009年12月15日
發(fā)明者于波, 侯健, 劉業(yè)俊, 張言輝, 杜慶軍, 王傳飛, 高達(dá) 申請(qǐng)人:中國(guó)石油大學(xué)(華東)