[發(fā)明專利]一種含一條直線封閉邊界油藏曲線擬合的處理方法

有效

• 申請?zhí)枺?/td>: 201410758423.X
• 申請日: 2014-12-10
• 公開（公告）號: CN104392068B
• 公開（公告）日: 2017-09-29
• 發(fā)明（設(shè)計）人: 王磊;何桃;程翔;趙錦波;胡麗琴
• 申請（專利權(quán)）人: 中國科學(xué)院合肥物質(zhì)科學(xué)研究院
• 主分類號: G06F17/50
• 分類號: G06F17/50
• 代理公司: 北京科迪生專利代理有限責(zé)任公司11251
• 代理人: 成金玉,孟卜娟

技術(shù)領(lǐng)域

本發(fā)明涉及計算機技術(shù)和石油工程領(lǐng)域，尤其涉及一種含一條直線封閉邊界油藏曲線擬合的處理方法。

背景技術(shù)

現(xiàn)有的試井分析方法主要包括：直線段分析方法、曲線擬合法及非線性回歸方法。

其中，現(xiàn)有曲線擬合法的本質(zhì)是通過試湊法實現(xiàn)曲線擬合；即，不斷的更改參數(shù)，每更改一次后，查看擬合結(jié)果，直至曲線的外形與弧度近似；再移動曲線使其重合。該方法的缺陷主要為：1)它是試湊方法，對解釋工程師水平有很高要求；2)需要多次反復(fù)的更改參數(shù)，并進行相應(yīng)計算，工作效率較低；3)當(dāng)參數(shù)范圍較大時更改的次數(shù)非常多，無法快速的進行曲線擬合。

發(fā)明內(nèi)容

本發(fā)明的技術(shù)解決問題：克服現(xiàn)有技術(shù)的不足，提供一種含一條直線封閉邊界油藏曲線擬合的處理方法，提高工作效率。

本發(fā)明的目的是通過以下技術(shù)方案實現(xiàn)的：

一種曲線擬合的處理方法，該方法包括：

將實測曲線及對應(yīng)的實測導(dǎo)數(shù)曲線繪制在雙對數(shù)圖版上，使得所述實測曲線及實測導(dǎo)數(shù)曲線的初始部分，與經(jīng)過原點且斜率為預(yù)定值s的斜線重合；該實測導(dǎo)數(shù)曲線的第一個水平段，與預(yù)定值為m的水平線重合；該實測導(dǎo)數(shù)曲線的結(jié)尾部分與數(shù)值為預(yù)定值n的水平線重合；

繪制一理論曲線及對應(yīng)的理論導(dǎo)數(shù)曲線，輸入無量綱組合參數(shù)C_De^2S的值，使理論導(dǎo)數(shù)曲線與實測導(dǎo)數(shù)曲線峰值重合；

接收光標(biāo)控制指令，根據(jù)該指令控制并移動光標(biāo)，當(dāng)光標(biāo)移動至該實測導(dǎo)數(shù)曲線第一個水平段的終點時，計算該位置對應(yīng)的無量綱組合參數(shù)L_De^S；

根據(jù)計算出的L_De^S參數(shù)計算并繪制一理論曲線及對應(yīng)的理論導(dǎo)數(shù)曲線，完成曲線擬合。

所述實測曲線包括：

含一條直線封閉邊界地層的井底壓力數(shù)據(jù)曲線，其無量綱井底壓力表達式為：

其中，J₀(u)與J₁(u)為0階與1階第一類Bessel函數(shù)(貝塞爾函數(shù))；Y₀(u)與Y₁(u)為0階與1階第二類Bessel函數(shù)；為無量綱壓力表達式；為無量綱時間；為無量綱井筒存儲常數(shù)；S為表皮系數(shù)；為無量綱直線封閉邊界距離；p_wf為井底流動壓力(MPa)；p_i為原始地層壓力(MPa)；k為地層滲透率(D)；h為地層有效厚度(m)；q為井產(chǎn)量(m³/d)；μ為流體粘度(mPa.s)；t為流動時間(h)；r_w為井筒半徑(m)；φ為地層孔隙度；C_t為綜合壓縮系數(shù)(1/Mpa)；C為井筒存儲參數(shù)(m³/MPa)；L為直線封閉邊界距離(m)；Ei為指數(shù)積分函數(shù)。

無量綱參數(shù)L_De^S與實測導(dǎo)數(shù)曲線第一個水平段終點對應(yīng)的t_D/C_D的關(guān)系為：

t_D/C_D＝0.3248(L_De^S)²/(C_De^2S)。

所述當(dāng)光標(biāo)移動至該實測導(dǎo)數(shù)曲線第一個水平段的終點時，計算該位置對應(yīng)的無量綱組合參數(shù)L_De^S包括：

當(dāng)光標(biāo)移動至該實測導(dǎo)數(shù)曲線第一個水平段的終點處時，根據(jù)光標(biāo)在屏幕中的位置計算所述實測導(dǎo)數(shù)曲線第一個水平段終點對應(yīng)的無量綱參數(shù)t_D/C_D及所述實測導(dǎo)數(shù)曲線對應(yīng)的L_De^S參數(shù)；

所述根據(jù)光標(biāo)在屏幕中的位置計算所述實測導(dǎo)數(shù)曲線第一個水平段終點對應(yīng)的無量綱參數(shù)t_D/C_D及所述實測導(dǎo)數(shù)曲線對應(yīng)的L_De^S參數(shù)包括：

根據(jù)光標(biāo)在屏幕中的位置計算所述實測導(dǎo)數(shù)曲線該實測導(dǎo)數(shù)曲線第一個水平段的終點對應(yīng)的t_D/C_D，其公式為：其中，x為光標(biāo)在屏幕中的橫坐標(biāo)，x₀為雙對數(shù)圖的坐標(biāo)原點在屏幕中的橫坐標(biāo)，l_x為圖在屏幕上的寬度；