[發明專利]油藏井間動態連通性反演方法

無效

• 申請號: 200910256524.6
• 申請日: 2009-12-15
• 公開（公告）號: CN101725346A
• 公開（公告）日: 2010-06-09
• 發明（設計）人: 侯健;杜慶軍;劉業俊;于波;張言輝;高達;王傳飛
• 申請（專利權）人: 中國石油大學(華東)
• 主分類號: E21B49/00
• 分類號: E21B49/00;G06N3/12
• 代理公司: 青島高曉專利事務所 37104
• 代理人: 楊大興

1.一種油藏井間動態連通性反演方法，將油藏中的注水井、生產井以及井間孔道看作一個完整的系統，根據該系統中生產井的產液量和注水井的注入量動態數據反演得到油藏井間動態連通性，其特征是，在注入信號的時滯性和衰減性規律研究的基礎上，應用信號處理技術，構建包含對注入信號進行預處理的井間動態連通性反演的多元回歸模型，將遺傳算法引入反演模型的優化求解中，反演得到表征油藏井間動態連通性的連通系數以及表征注入信號時滯性和衰減性的時間常數，具體按以下步驟實現：

(1)注采動態數據的選取：從油田月度生產動態數據報表中直接獲取注入井的月度注水量、生產井的月度產液量數據，為了具有更好的反演效果，注采動態數據的選取考慮動態數據要有幅度為±20％的波動，各單井注采量保持連續性，選取數據時間段應大于1年；

(2)生成初始群體：根據時間常數與導壓系數、井距以及孔隙度的回歸關系式估計時間常數初始值，式中，τ為時間常數，month；φ為孔隙度，小數；d為井距，m；為導壓系數，cm²/s，以此為基礎，產生一組實數型的時間常數值構成一個個體，個體實數個數即為時間常數待定參數個數，然后在時間常數初始值的50％～150％區間內，按均勻分布隨機產生N_l個這種個體構成初始群體；

(3)求解井間動態連通系數：結合時間常數群體數據應用褶積器對注入量信號進行預處理，得到修正后的注入數據，然后應用最小二乘法求得表征井間動態連通程度的多元線性回歸的權重系數，權重系數的最小二乘估計滿足下式：

∂∂βij[Σk=1Nk(Qj(k)-Q^j(k))2]=0]]>

式中，Q_j為第j口生產井產液量實際值，m³/d；為第j口生產井產液量估計值，m³/d；β_ij為第j口生產井和第i口注水井間的多元線性回歸權重，定義為井間動態連通系數；N_k為數據點總數；

(4)適應度函數定義：目標函數F_j可定義為

Fj(τp,τ1j,τ2j,······,τNIj)=Σk=1Nk|(Q^j(k)-Qj(k))/Qj(k)|,(j=1,···,NP)]]>

式中，τ_p為產液量初值影響的時間常數，month；τ_ij為對應于第j口生產井和第i口注水井間的時間常數，month；N_I口注水井井數；N_P口生產井井數。該目標函數值反映的是第j口生產井產液量估計值與實際值之間的累積相對誤差，以產液量估計值與實際值的相對誤差定義適應度函數，應保證產液量估計值與實際值之間的誤差越大，適應度越小，取適應度函數f_j＝1/F_j；

(5)適應性評估：以適應度函數來確定個體適應度值以評估其優劣，并以此作為遺傳操作的依據，計算時間常數群體中各個體的適應度f_jl(l＝1，…，N_l)；

(6)遺傳變異操作：進行遺傳變異操作，操作的目的是根據進化原則從當前群體中選出優良的個體，從而形成下一代，首先進行停止規則判斷，若發現占群體80％～90％以上的個體在相對誤差為設定值范圍內是同一個體，或者算法迭代步數超出設定的閾值，就終止算法迭代，則當前最優個體組成即為所求時間常數的待定值，進而再重新計算井間動態連通系數，輸出結果，否則進行選擇、交叉和變異操作，產生子代，重復步驟(3)～(6)，繼續進化；

(7)輸出反演結果：將得到的井間動態連通系數繪制動態連通圖。