致密變形介質(zhì)油藏五點(diǎn)面積井網(wǎng)水驅(qū)動(dòng)態(tài)分析

曹耐，董平川*，雷剛

1 中國(guó)石油大學(xué)(北京)石油工程學(xué)院，北京 102249

2 College Petroleum of Engineering & Geoscience, King Fahd University of Petroleum and Minerals, Dhahran 31261, Kingdom of Saudi Arabia

0 引言

致密變形介質(zhì)油藏的滲流特征與常規(guī)油藏的滲流特征顯著不同[1-2]，有著其特有的滲流規(guī)律，主要表現(xiàn)為2個(gè)方面：①孔喉細(xì)小、孔喉比大、毛管壓力的影響顯著，非均質(zhì)性強(qiáng)，致使流體滲流偏離Darcy定律，出現(xiàn)非線性-擬線性的非Darcy滲流規(guī)律[3-4]；②實(shí)際生產(chǎn)中，儲(chǔ)層骨架產(chǎn)生部分或全部的不可逆變形，巖石滲透率和孔隙度降低[5-6]，對(duì)油田的動(dòng)態(tài)特征有明顯的影響[7]。而且在實(shí)際開(kāi)發(fā)過(guò)程中也發(fā)現(xiàn)了很多問(wèn)題：難以形成有效的驅(qū)動(dòng)壓力系統(tǒng)，地層壓力下降快，產(chǎn)量迅速遞減，原油有效動(dòng)用率低[8-10]。

國(guó)內(nèi)外學(xué)者對(duì)介質(zhì)變形理論做了大量研究：Samaniego等考慮巖石應(yīng)力敏感和流體性質(zhì)，對(duì)儲(chǔ)層流體非穩(wěn)態(tài)滲流特征進(jìn)行了分析[11]。熊建和李凌峰對(duì)巖心驅(qū)替實(shí)驗(yàn)結(jié)果進(jìn)行統(tǒng)計(jì)學(xué)分析，發(fā)現(xiàn)低滲透氣藏巖心絕對(duì)滲透率隨有效應(yīng)力增加指數(shù)下降，并基于此得到了變形介質(zhì)低滲透氣藏產(chǎn)量預(yù)測(cè)公式[12]。Samimi和Ghorbani等采用無(wú)單元數(shù)值模擬方法對(duì)變形多孔介質(zhì)中多相流體流動(dòng)規(guī)律進(jìn)行分析[13-14]。王美楠等結(jié)合二次梯度項(xiàng)及動(dòng)邊界理論，對(duì)低滲透變形介質(zhì)油藏滲流規(guī)律進(jìn)行了探討和分析[15]。郭平和張楠等建立了啟動(dòng)壓力梯度和應(yīng)力敏感效應(yīng)共同影響下的單井產(chǎn)能模型，并分析了非達(dá)西滲流效應(yīng)對(duì)低滲透儲(chǔ)層單井產(chǎn)能的影響[16-17]。雷剛等基于致密砂巖氣藏中氣體非達(dá)西滲流特征，考慮儲(chǔ)層絕對(duì)滲透率應(yīng)力敏感效應(yīng)影響，推導(dǎo)了致密砂巖在擬穩(wěn)態(tài)流動(dòng)階段的垂直裂縫井產(chǎn)能方程[18]。Pesavento等對(duì)變形介質(zhì)中多相流動(dòng)研究進(jìn)行了總結(jié)和分析[19]。

然而，上述研究工作基本上都以巖心實(shí)驗(yàn)、介質(zhì)變形系數(shù)對(duì)單井產(chǎn)能和滲流規(guī)律影響為研究目標(biāo)；針對(duì)致密油藏面積井網(wǎng)，考慮介質(zhì)變形影響的研究卻不多[20-22]。齊亞?wèn)|等建立了特低滲透斷塊油藏不規(guī)則三角形井網(wǎng)有效動(dòng)用系數(shù)計(jì)算模型，并以某特低滲透斷塊油藏為例，對(duì)不規(guī)則三角形井網(wǎng)有效動(dòng)用系數(shù)進(jìn)行了計(jì)算分析[23]。呂棟梁等建立了低滲透油田反九點(diǎn)井網(wǎng)面積波及效率數(shù)學(xué)模型[24]。郭粉轉(zhuǎn)等考慮啟動(dòng)壓力梯度影響，建立數(shù)學(xué)模型對(duì)低滲透油田四點(diǎn)井網(wǎng)和五點(diǎn)井網(wǎng)面積波及效率進(jìn)行了計(jì)算和分析[25-26]。上述研究工作對(duì)于分析低滲透油藏面積井網(wǎng)流動(dòng)規(guī)律具有重大意義，但由于上述文獻(xiàn)并未考慮儲(chǔ)層的介質(zhì)變形，建立的理論模型不適于變形介質(zhì)油藏面積井網(wǎng)生產(chǎn)動(dòng)態(tài)分析，需增加對(duì)變形介質(zhì)油藏面積井網(wǎng)的研究及討論。因此，從滲流理論上揭示變形介質(zhì)油藏井網(wǎng)的流動(dòng)特征是提高油藏原油有效動(dòng)用率亟待解決的首要問(wèn)題之一。

筆者充分考慮致密油藏油水兩相滲流的非線性特征及多孔介質(zhì)變形特征，推導(dǎo)了考慮介質(zhì)變形面積井網(wǎng)的地層壓力、勢(shì)函數(shù)、流函數(shù)和平面速度計(jì)算表達(dá)式，根據(jù)物質(zhì)平衡理論及特征線方法，得到考慮介質(zhì)變形系數(shù)面積井網(wǎng)單元流場(chǎng)分布、油藏波及系數(shù)、見(jiàn)水時(shí)間和剩余油飽和度分布等生產(chǎn)指標(biāo)，并研究了不同介質(zhì)變形系數(shù)對(duì)油藏生產(chǎn)影響。

1 致密介質(zhì)變形模型驅(qū)替理論

建立考慮介質(zhì)變形油水兩相滲流數(shù)學(xué)模型的假設(shè)條件為：①水平均質(zhì)等厚的無(wú)限大地層布有無(wú)限延伸的五點(diǎn)井網(wǎng)；②地層只有油水兩相，流體在多孔介質(zhì)中流動(dòng)均為二維流動(dòng)；③不考慮毛管壓力影響；④油井以定液量生產(chǎn)，注入井注入量與油井產(chǎn)液量相等；⑤多孔介質(zhì)發(fā)生變形，儲(chǔ)層滲透率隨壓力發(fā)生變化，而孔隙度隨壓力變化較小不予考慮。

根據(jù)基本假設(shè)條件，可以得到多孔介質(zhì)中油水兩相流體控制方程為：

油水兩相速度方程分別為：

式中k(p)為巖石絕對(duì)滲透率。研究結(jié)果表明k(p)為巖石孔隙結(jié)構(gòu)、彈性模量和泊松比等參數(shù)的函數(shù)。目前，國(guó)內(nèi)外學(xué)者主要采用指數(shù)函數(shù)來(lái)表征巖石絕對(duì)滲透率隨有效應(yīng)力的變化規(guī)律[27-28]，關(guān)系表達(dá)式為：

將方程(1～4)聯(lián)立并整理(qw+qo=0)，可以得到：

基于不變流線假設(shè)，計(jì)算地層壓力場(chǎng)時(shí)，不考慮流度在空間上變化，方程(5)可以變形為：

其中

2 模型計(jì)算

2.1 壓力計(jì)算

方程(6)為調(diào)和方程，采用疊加原理，計(jì)算無(wú)限大地層矩形5點(diǎn)井網(wǎng)穩(wěn)態(tài)壓力分布。如圖1所示，根據(jù)所建立直角坐標(biāo)系可以得到，油井坐標(biāo)為(2nd+md,ma)，注入井坐標(biāo)(2nd+md+d,ma)，其中(m=0，±1，±2，...；n=0，±1，±2，…)。對(duì)于地層中任意點(diǎn)M(x,y)，第m排井排中第n口油井在M點(diǎn)的擬勢(shì)函數(shù)為：

其中

第m排井排中第n口注入井在M點(diǎn)的擬勢(shì)函數(shù)為：

圖1 矩形五點(diǎn)面積井網(wǎng)示意圖Fig. 1 The sketch map of rectangularfive-spot well pattern

其中

根據(jù)勢(shì)疊加原理，無(wú)窮井網(wǎng)在M點(diǎn)擬勢(shì)函數(shù)為：

其中

對(duì)式(9)進(jìn)行處理，可以得到：

由于ψ=λtU(x,y)，可以得到：

根據(jù)無(wú)因次壓力定義，可以得到地層壓力：

其中

若注入井注入壓力為piwf，生產(chǎn)井井底壓力為pwf，則根據(jù)壓力計(jì)算公式，可得到無(wú)因次注采壓差為：

2.2 流函數(shù)計(jì)算

令vt=vo+vw，由式(2)、(3)、(4)可得：

通過(guò)方程(10)可以得到五點(diǎn)面積井網(wǎng)任意位置在x，y方向滲流速度：

由于dy=adyD，通過(guò)對(duì)速度積分可以得到流函數(shù)為：

式中：

2.3 有效動(dòng)用系數(shù)計(jì)算

對(duì)于五點(diǎn)井網(wǎng)，x，y軸上主流線均為直線。主流線任意點(diǎn)在x方向上的滲流速度為：

由此可確定x方向上見(jiàn)水時(shí)間：

同理可確定y方向上見(jiàn)水時(shí)間：

所以可得到見(jiàn)水時(shí)間：

根據(jù)(12)，并利用物質(zhì)平衡方法，可得到井網(wǎng)單元有效動(dòng)用系數(shù)為：

由于地層均質(zhì)，井網(wǎng)單元有效動(dòng)用系數(shù)即為無(wú)窮大面積井網(wǎng)有效動(dòng)用系數(shù)。

2.4 含水飽和度計(jì)算

由方程(1)、(10)可得到：

微分方程(14)特征微分方程為：

沿特征線有：

式(16)說(shuō)明沿特征線含水飽和度不變，即特征線就為等飽和度線。通過(guò)式(15)得到等飽和度面運(yùn)動(dòng)特征方程：

求解特征方程(17)可以得到任意時(shí)間地層含水飽和度分布。

3 計(jì)算實(shí)例

為研究?jī)?chǔ)層有效動(dòng)用情況，以長(zhǎng)慶油田某低滲透砂巖實(shí)際儲(chǔ)層特性數(shù)據(jù)及生產(chǎn)數(shù)據(jù)為例，對(duì)推導(dǎo)的數(shù)學(xué)模型進(jìn)行了計(jì)算。通過(guò)室內(nèi)巖心滲透率應(yīng)力敏感實(shí)驗(yàn)，得到介質(zhì)變形系數(shù)為0.014 MPa-1。計(jì)算分析所采用的基礎(chǔ)數(shù)據(jù)如表1所示。

3.1 典型壓力場(chǎng)、流場(chǎng)分布

通過(guò)計(jì)算可以得到地層壓力分布如圖2、圖3所示。圖2為地層某區(qū)塊壓力分布，圖3為地層某五點(diǎn)面積井網(wǎng)單元壓力分布。地層壓力場(chǎng)分布圖展現(xiàn)了面積井網(wǎng)生產(chǎn)特征，注入井附近區(qū)域壓力值較高，生產(chǎn)井附近區(qū)域地層壓力值較低，由于面積井網(wǎng)對(duì)稱性及地層均質(zhì)性，地層壓力表現(xiàn)為對(duì)稱分布。

圖4為地層某五點(diǎn)面積井網(wǎng)單元流場(chǎng)分布，地層流場(chǎng)分布圖展現(xiàn)了面積井網(wǎng)生產(chǎn)特征，近井附近區(qū)域流線密，遠(yuǎn)井附近區(qū)域流線疏，主要原因在于近井附近區(qū)域生產(chǎn)壓差大，從而造成流線密集分布。

3.2 介質(zhì)變形系數(shù)對(duì)流場(chǎng)分布的影響

對(duì)于不同介質(zhì)變形系數(shù)，在流函數(shù)值相等情況下，可得到地層流線分布示意圖如圖5、圖6、圖7所示。

表1 模型基礎(chǔ)數(shù)據(jù)Table 1 The reservoir parameters of the model

圖2 區(qū)塊壓力分布(α=0.014 MPa-1)Fig. 2 Pressure field distribution of formation area block(α=0.014 MPa-1)

由圖5、圖6和圖7可見(jiàn)，介質(zhì)變形系數(shù)越大，流線波及面積越小。主要原因在于隨著介質(zhì)變形系數(shù)增大，地層滲流阻力增大，滲流速度減小，從而導(dǎo)致流線波及面積減小。

3.3 介質(zhì)變形系數(shù)對(duì)見(jiàn)水時(shí)間、動(dòng)用程度的影響

根據(jù)得到壓力計(jì)算公式，對(duì)于不同介質(zhì)變形系數(shù)，產(chǎn)量與生產(chǎn)壓差關(guān)系曲線圖如圖8所示。

從圖8可以看出，隨著介質(zhì)變形系數(shù)增大，產(chǎn)量增幅隨著生產(chǎn)壓差增大而減小。主要原因在于隨著生產(chǎn)壓差增大，地層滲透率降低，生產(chǎn)阻力變大，從而使得產(chǎn)量與生產(chǎn)壓差關(guān)系不再滿足達(dá)西定律。

圖3 五點(diǎn)面積井網(wǎng)單元壓力分布(α=0.014 MPa-1)Fig. 3 Pressure field distribution of five-spot well pattern unit(α=0.014 MPa-1)

圖4 五點(diǎn)井網(wǎng)單元流場(chǎng)分布(α=0.014 MPa-1)Fig. 4 Streamline distribution of 5-spot well pattern unit(α=0.014 MPa-1)

通過(guò)式(11)、(12)可計(jì)算見(jiàn)水時(shí)間隨介質(zhì)變形系數(shù)變化規(guī)律如圖9所示。

從圖9可以看出，見(jiàn)水時(shí)間隨介質(zhì)變形系數(shù)增加而增大，說(shuō)明多孔介質(zhì)變形增大了滲流阻力，降低了流體滲流速度。

對(duì)于本文所給實(shí)例，當(dāng)介質(zhì)變形系數(shù)為0.014 MPa-1時(shí)，油井見(jiàn)水時(shí)間為149.6 d，現(xiàn)場(chǎng)分析結(jié)果顯示油井見(jiàn)水時(shí)間為166 d，本文計(jì)算結(jié)果相對(duì)誤差為9.9 %，符合性較好，說(shuō)明本文計(jì)算模型正確。根據(jù)式(13)可計(jì)算地層有效動(dòng)用系數(shù)為11.87 %。

圖5 流場(chǎng)分布(α=0.013 MPa-1)Fig. 5 Streamline distribution (α=0.013 MPa-1)

圖6 流場(chǎng)分布(α=0.017 MPa-1)Fig. 6 Streamline distribution (α=0.017 MPa-1)

圖7 流場(chǎng)分布(α=0.028 MPa-1)Fig. 7 Streamline distribution (α=0.028 MPa-1)

圖8 產(chǎn)量與生產(chǎn)壓差關(guān)系圖Fig. 8 The graph of yield and production pressure drop

圖9 見(jiàn)水時(shí)間與介質(zhì)變形系數(shù)關(guān)系圖Fig. 9 The breakthrough time and media deformation factor

3.4 介質(zhì)變形系數(shù)對(duì)剩余油飽和度的影響

在不同介質(zhì)變形系數(shù)下，通過(guò)式(17)計(jì)算生產(chǎn)50 d后地層剩余油飽和度分布，結(jié)果如圖10和圖11所示。由圖10、圖11可見(jiàn)，隨著介質(zhì)變形系數(shù)增大，地層剩余油飽和度增大，說(shuō)明應(yīng)力敏感效應(yīng)不利于油藏生產(chǎn)，且應(yīng)力敏感效應(yīng)越強(qiáng)烈，油藏采出程度越低。

3.5 介質(zhì)變形系數(shù)對(duì)水驅(qū)前緣的影響

在不同介質(zhì)變形系數(shù)下，通過(guò)式(17)計(jì)算生產(chǎn)130 d后地層含水飽和度分布，結(jié)果如圖12所示。由圖12可見(jiàn)，隨著介質(zhì)變形系數(shù)增大，水驅(qū)前緣移動(dòng)距離減小。在同一位置，含水飽和度隨變形系數(shù)的增大而減小。說(shuō)明應(yīng)力介質(zhì)變形對(duì)水驅(qū)效率有不利影響，巖石應(yīng)力效應(yīng)越強(qiáng)烈，水驅(qū)油效果越差。

圖10 剩余油飽和度分布(α=0.028 MPa-1)Fig. 10 Remaining oil saturation distribution (α=0.028 MPa-1)

圖11 剩余油飽和度分布(α=0.014 MPa-1)Fig. 11 Remaining oil saturation distribution (α=0.01 MPa-1)

圖12 含水飽和度與介質(zhì)變形系數(shù)關(guān)系圖Fig. 12 The water saturation with media deformation factor

4 結(jié)論

(1)基于油水兩相滲流方程，結(jié)合介質(zhì)變形方程，推導(dǎo)了考慮介質(zhì)變形的五點(diǎn)面積井網(wǎng)壓力、流函數(shù)解析公式。

(2)根據(jù)地層流場(chǎng)分布，得到油井見(jiàn)水時(shí)間，結(jié)合物質(zhì)平衡理論，推導(dǎo)出無(wú)因次井網(wǎng)單元有效動(dòng)用系數(shù)計(jì)算公式，進(jìn)而可以定量計(jì)算低滲透介質(zhì)變形油田五點(diǎn)面積井網(wǎng)的有效動(dòng)用系數(shù)。理論計(jì)算見(jiàn)水時(shí)間與實(shí)際生產(chǎn)情況基本一致，相對(duì)誤差為9.9 %，符合性較好，說(shuō)明本文計(jì)算模型正確。根據(jù)油水兩相滲流方程，得到了等飽和度面運(yùn)動(dòng)特征方程，進(jìn)而可以得到地層剩余油飽和度分布。

(3)介質(zhì)變形系數(shù)造成面積井網(wǎng)動(dòng)用程度降低，油藏采出程度降低。隨著介質(zhì)變形系數(shù)增大，地層滲流阻力增大，滲流速度減小，從而導(dǎo)致流線波及面積減小、油井見(jiàn)水時(shí)間增長(zhǎng)；同時(shí)，由于介質(zhì)變形系數(shù)增大，生產(chǎn)阻力隨生產(chǎn)壓差的增大而增大，水驅(qū)前緣移動(dòng)距離減小，進(jìn)而對(duì)儲(chǔ)層滲透率產(chǎn)生重大影響，使產(chǎn)量與生產(chǎn)壓差的關(guān)系不再遵循達(dá)西定律，產(chǎn)量增幅變小。

符號(hào)注釋：

式中：v—為流體滲流速度，m/d；k0—巖石初始滲透率，μm2；k(p)—地層壓力為p時(shí)地層絕對(duì)滲透率，μm2；p—地層壓力，MPa；p0—初始地層壓力，MPa；kr—相對(duì)滲透率；μ—流體黏度，mPa·s；α—介質(zhì)變形系數(shù)，1/MPa；q—源匯項(xiàng)，1/d；φ—地層孔隙度，%；S—流體飽和度，%；t—時(shí)間，d；U—無(wú)因次壓力；λt—混合流體流度，μm2/(mPa·s)；▽—Laplace算子，表示內(nèi)積；Q—油井產(chǎn)量，m3/d；C,C1,C2—常數(shù)；r—距離，m；h—油層厚度，m；ψ—擬勢(shì)函數(shù)，m2/(d·MPa)；vt—混合流體速度，m/d；T—見(jiàn)水時(shí)間，d；A—井網(wǎng)單元面積；ω—積分變量；fw—Leverrett分流函數(shù)；下標(biāo)：o—油相，w—水相；m—井排排數(shù)，n—油井或注入井井?dāng)?shù)，M—地層中任意點(diǎn)。