孫召勃,賈曉飛,楊馥榕
(1.中國海洋石油國際有限公司,北京 100028;2.中海石油(中國)有限公司天津分公司,天津 300459)
渤海油田稠油資源具有儲量規(guī)模大、油藏埋藏深、黏度范圍廣等特點,水平井是開發(fā)稠油油藏的有效手段之一[1-3]。目前已有學者提出了稠油油藏水平井產(chǎn)能計算公式,但其中只有極少數(shù)人考慮到井筒摩阻損失,且相關(guān)的井筒摩阻計算均采用了牛頓流體管流模型,不符合非牛頓流體的管流特征,夸大了摩阻的影響[4-13]。客觀準確計算非牛頓流體特性和摩阻損失的影響,對優(yōu)化水平井長度、準確預測產(chǎn)能意義重大[14-20]。
該文同時考慮稠油油藏的非牛頓特性和水平井段的摩阻損失,利用滲流阻力法建立了冪律型稠油油藏水平井產(chǎn)能計算模型,分析了水平井長度、冪律指數(shù)和稠度系數(shù)對水平井產(chǎn)能的影響,為評價冪律型稠油油藏水平井產(chǎn)能提供理論指導。
冪律型流體的流變性表征方程為:
(1)
式中:τ為剪切應力,MPa;n為冪律指數(shù);H為稠度系數(shù),mPa·sn;r為徑向距離,m;u為速度,m/s。
冪律型稠油油藏中流體的運動方程為:
(2)
其中:
μb=μeffvn-1
(3)
(4)
式中:v為滲流速度,m/s;k為滲透率,mD;p為壓力,MPa;μb為冪律流體的視黏度,mPa·s;μeff為等效黏度,mPa·s;φ為孔隙度。
將水平井的滲流場簡化為2個相互連接的滲流區(qū)域,第1流動區(qū)域為平面上從驅(qū)動邊界流向水平井方向的滲流區(qū),該區(qū)域受儲層厚度控制;第2流動區(qū)域為垂直平面內(nèi)圍繞水平井段向水平井的滲流區(qū),該區(qū)域滲流受水平井長度控制。
1.1.1 外部平面內(nèi)徑向流
外部平面滲流場可近似為穩(wěn)定的平面徑向流,相當于泄油半徑為re、厚度為h的圓柱體泄油區(qū)向井筒半徑為L/4的“直井”滲流。
由式(2)和式(3)可得外部穩(wěn)定滲流時徑向流方程為:
(5)
式中:qh為外部平面內(nèi)徑向流的流量,m3/d;Bo為原油體積系數(shù);h為儲層厚度,m。
對式(5)進行分離變量法求解,可得外部平面內(nèi)徑向流流量為:
(6)
式中:Δp1為外部平面內(nèi)徑向流區(qū)域的壓差,MPa;re為供給半徑,m;L為水平井的水平井段長度,m。
1.1.2 內(nèi)部垂直面內(nèi)徑向流
在垂直平面上,水平井滲流場受水平井段L控制,垂直徑向流外部邊界半徑為h/(2π),內(nèi)部邊界半徑為井半徑rw。由式(2)和式(3)可得穩(wěn)定滲流時徑向流方程為:
(7)
式中:qv為內(nèi)部垂直面內(nèi)徑向流的流量,m3/d。
對式(7)進行分離變量法求解,可得內(nèi)部垂直面內(nèi)徑向流流量為:
(8)
式中:Δp2為內(nèi)部垂直面內(nèi)徑向流的壓差,MPa。
由于原油黏度較大,流量相對不高,在水平井段內(nèi)通常出現(xiàn)層流,即流體只沿著軸向流動而無橫向運動。
水平井段內(nèi)層流的流量為:
(9)
式中:qf為水平井段內(nèi)層流的流量,m3/d;D為水平段井筒的直徑,m。
將式(1)帶入式(9)并分部積分可得:
(10)
(11)
式中:Δp3為井筒附加壓降,MPa。
式(11)變形可得:
(12)
(13)
(14)
式中:q為水平井產(chǎn)量,m3/d;Δp為總的生產(chǎn)壓差,MPa。
由于n趨于或等于1時,存在如下關(guān)系式:
(15)
(16)
因此,當流體逐漸由冪律型流體趨于牛頓流體時,若不考慮摩阻損失,式(14)可逐漸接近直至變?yōu)锽orisov水平井產(chǎn)能公式[14]:
(17)
而對于層流摩阻而言,在n=1時,水平井段井筒內(nèi)摩阻損失造成的滲流阻力退化為常見的牛頓流體層流造成的滲流阻力:
(18)
式中:R為水平井段內(nèi)的滲流阻力,MPa·d·m-3。
綜合以上分析,n=1時,式(13)可以簡化為考慮摩阻損失的常規(guī)油藏水平井產(chǎn)能公式:
(19)
由上述分析可知,式(13)是適用范圍更廣的考慮摩阻損失的冪律型稠油油藏水平井產(chǎn)能公式,Borisov水平井產(chǎn)能公式和牛頓流體摩阻損失公式都是其特例。
分別利用上述產(chǎn)能公式對水平井產(chǎn)量進行實例計算,對比分析各公式的適應性。
已知油藏基本參數(shù)為:re=1 000 m,rw=0.12 m,L=500 m,k=3 000 mD,h=30 m,φ=0.3,Bo=1.08,n=0.9,H=100 mPa·sn。水平井生產(chǎn)壓差Δp=2 MPa,實際產(chǎn)量為38.50 m3/d。采用不同的方法計算水平井產(chǎn)量,結(jié)果見表1。
表1 各種方法計算的水平井產(chǎn)量結(jié)果Table 1 Production results of horizontal wells calculated by various methods
由表1可以看出,式(17)和式(19)計算結(jié)果遠大于實際值,說明非牛頓特性對產(chǎn)能影響很大,不可忽略;考慮非牛頓特性的后3種計算方法所得結(jié)果誤差明顯降低,其中考慮非牛頓流體管流摩阻的式(13)計算誤差最小,僅為-5.77%,其他2個公式相對誤差的絕對值大于10%,尤其是地層滲流阻力采用非牛頓流體而管流摩阻損失采用牛頓流體的最后一種方法,放大了摩阻損失的影響,誤差較大,低估了水平井產(chǎn)能。綜上對比分析,該文所建立的綜合形式得到的式(13)考慮因素最為全面、準確,誤差最小,滿足工程要求,可以用于礦場實踐。
當其他參數(shù)一定時,計算不同水平井長度下的產(chǎn)能,如圖1所示。圖1中藍線為考慮了摩阻損失的冪律型稠油油藏水平井產(chǎn)能公式,紅線為對應公式在不考慮摩阻損失情況下的形式。由圖1可以看出,水平井越長,兩公式求得的水平井產(chǎn)量越高。其中,考慮了摩阻損失的公式計算出的產(chǎn)能相對更低,且隨著水平井長度的不斷增加,產(chǎn)能的增加幅度也逐漸變小。當水平井長度超過600 m后,由于摩阻損失造成的產(chǎn)能降低超過了20%。根據(jù)Dikken研究成果可知,此地層流體參數(shù)下水平井最優(yōu)長度為600 m,因此對于稠油油藏而言,水平井不宜過長。
圖1 水平井長度對產(chǎn)能的影響Fig.1 Influence of horizontal well length on productivity
當其他參數(shù)一定時,計算不同冪律指數(shù)對水平井產(chǎn)能的影響,如圖2所示。由圖2可以看出,隨著冪律指數(shù)的降低,考慮了摩阻損失的冪律型稠油油藏水平井產(chǎn)能公式和不考慮摩阻損失情況下的產(chǎn)能公式計算出的產(chǎn)量均明顯降低,且當冪律指數(shù)下降至0.9時,兩者的產(chǎn)能降低都超過了80%。冪律指數(shù)的下降表示原油的非牛頓特性增強,從而影響了水平井的產(chǎn)能。因此,在礦場實際生產(chǎn)中,對于強非牛頓特性的稠油,通過物理化學等手段降低非牛頓特性可以顯著提高水平井產(chǎn)能,是重要的增產(chǎn)措施之一。
圖2 冪律指數(shù)對產(chǎn)能的影響Fig.2 Influence of power-law index on productivity
當其他參數(shù)一定時,計算不同稠度系數(shù)對水平井產(chǎn)能的影響,如圖3所示。由圖3可以看出,隨著稠度系數(shù)的增大,兩公式計算出的水平井產(chǎn)量均不斷降低,且降低幅度逐漸減小。由于稠度系數(shù)反應了非牛頓流體的黏性程度,隨著稠度系數(shù)增大,相應的剪切黏度和有效黏度也會增大,滲流阻力增加,因此導致水平井產(chǎn)能降低。
圖3 稠度系數(shù)對產(chǎn)能的影響Fig.3 Influence of consistency coefficient on productivity
S油田為某典型稠油油藏,部署4口水平井進行局部細分層系調(diào)整。已知基本參數(shù)為:rw=0.12 m,φ=0.31,Bo=1.08,n=0.95,H=186 mPa·sn,其他參數(shù)及產(chǎn)量見表2。采用新方法計算的產(chǎn)能結(jié)果相對誤差小于7%,滿足礦場要求。
表2 4口井其他基本參數(shù)
1)利用冪律型稠油油藏水平井產(chǎn)能計算公式,研究了水平井長度對稠油油藏產(chǎn)能的影響。結(jié)果表明,水平井越長,產(chǎn)能增加幅度越小,水平井長度超過600 m后,摩阻損失造成的產(chǎn)能降低將超過20%。
2)稠油的非牛頓流體特性對水平井產(chǎn)能影響顯著,原油的非牛頓特性越強,水平井產(chǎn)能下降越明顯,當冪律指數(shù)下降至0.9時,產(chǎn)能降低超過80%。
3)礦場實例表明,綜合考慮了非牛頓流體特性和摩阻損失后建立的冪律型稠油油藏水平井產(chǎn)能公式計算誤差最小,僅為-5.77%,較常用公式適應性更強,計算結(jié)果更準確,更能滿足工程要求。