考慮啟動(dòng)壓力梯度的低滲透滲流模型改造

劉 剛 孫建博 尹錦濤 徐 杰

陜西延長石油(集團(tuán))有限責(zé)任公司研究院,陜西 西安 710075

?

考慮啟動(dòng)壓力梯度的低滲透滲流模型改造

劉 剛 孫建博 尹錦濤 徐 杰

陜西延長石油(集團(tuán))有限責(zé)任公司研究院,陜西 西安 710075

如何充分反映啟動(dòng)壓力梯度的影響一直是低滲透油藏?cái)?shù)值模擬研究領(lǐng)域的難點(diǎn),以低滲透油藏滲流機(jī)理為研究起點(diǎn),通過對低滲透油藏滲流模型的改造,得到了一種啟動(dòng)壓力梯度在數(shù)值模擬方程中新的表現(xiàn)形式,并研發(fā)出相應(yīng)的低滲透油藏?cái)?shù)值模擬編輯器,通過與專業(yè)石油軟件進(jìn)行對比驗(yàn)證了該模型的高精確度,進(jìn)一步結(jié)合實(shí)際油田資料完成了該模型的實(shí)際生產(chǎn)應(yīng)用,結(jié)果表明改造后的模型很大程度上加強(qiáng)了啟動(dòng)壓力梯度對低滲透數(shù)值模擬過程的影響程度,并且能夠更加準(zhǔn)確地反映出低滲透油藏的實(shí)際生產(chǎn)動(dòng)態(tài)。

低滲透油藏;啟動(dòng)壓力梯度;模型改造;油藏?cái)?shù)值模擬

0 前言

作為低滲透油藏中最重要的滲流特征,啟動(dòng)壓力梯度一直是低滲透油藏的研究重點(diǎn)。從20世紀(jì)90年代開始程林松、吳文瑞等廣大學(xué)者通過理論分析和實(shí)驗(yàn)研究對啟動(dòng)壓力梯度進(jìn)行了深入研究,分析出啟動(dòng)壓力梯度與儲(chǔ)層、流體特征之間存在一定的相關(guān)關(guān)系,并總結(jié)出關(guān)于啟動(dòng)壓力梯度的多種表達(dá)公式[1-8]。然而我國低滲透油藏儲(chǔ)層非均質(zhì)極強(qiáng),若精確反映整個(gè)油藏的啟動(dòng)壓力梯度,對于網(wǎng)格數(shù)量級(jí)達(dá)到數(shù)十萬級(jí)別的油藏?cái)?shù)值模擬耗時(shí)繁瑣,實(shí)效性較差。為求簡便高效,應(yīng)用擬啟動(dòng)壓力梯度[9-10]參與數(shù)值模擬成為廣大學(xué)者的首選,但筆者通過研究發(fā)現(xiàn)該方法同樣存在缺陷:處在微分項(xiàng)內(nèi)的擬啟動(dòng)壓力梯度值在前后網(wǎng)格數(shù)值差分過程中被相互抵消,極大地削弱了啟動(dòng)壓力梯度的影響程度。

本文通過對低滲透油藏儲(chǔ)層流體性質(zhì)和滲流機(jī)理進(jìn)行深入研究[11-13],實(shí)現(xiàn)對低滲油藏滲流模型的改造,力求使啟動(dòng)壓力梯度對低滲透油藏的影響在數(shù)值模擬過程中得到較大程度的體現(xiàn)。

1 低滲透滲流模型建立

以低滲透油藏的質(zhì)量守恒方程和流體運(yùn)動(dòng)方程為基礎(chǔ),建立了考慮啟動(dòng)壓力梯度的低滲透油藏的滲流模型[14-15]:

(1)

2 低滲透滲流模型改造

低滲透油藏微可壓縮流體的密度一般表征為壓力的函數(shù),因此可將式(1)左端的第一項(xiàng)做微分展開:

ρm

(2)

其中密度關(guān)于壓力的函數(shù)通過流體壓縮系數(shù)Cm的定義推導(dǎo)而得[16]：

(3)

式中:V為體積,cm3。

微可壓縮流體的Cm作常量處理,由式(3)得到：

(4)

從而得到流體密度關(guān)于壓力的函數(shù)表達(dá)式如下:

(5)

將式(5)代入到微分展開式(2)中,整理可得:

(6)

ρmp+qm

(7)

3 模型精確度驗(yàn)證

在改造后的低滲透油藏滲流模型原理基礎(chǔ)上,利用Visual Basic軟件開發(fā)出了相應(yīng)的低滲透油藏滲流數(shù)值模擬器[16-22],該模擬器可用于實(shí)際低滲透油藏?cái)?shù)值模擬研究，見圖1。

在儲(chǔ)層、流體條件及工作制度相同的情況下,應(yīng)用商業(yè)軟件Eclipse與所開發(fā)的低滲透油藏?cái)?shù)值模擬器進(jìn)行單井開發(fā)單元數(shù)值模擬對比分析,完成模型器精確度驗(yàn)證,具體的儲(chǔ)層和流體數(shù)據(jù)見表1～2。生產(chǎn)制度規(guī)定如下:定產(chǎn)油量生產(chǎn),日產(chǎn)油量為10 t/d,當(dāng)井底壓力降到2 MPa時(shí),改變生產(chǎn)制度,轉(zhuǎn)為定井底壓力生產(chǎn)。

圖1 低滲透油藏?cái)?shù)值模擬器

表1 研究區(qū)儲(chǔ)層數(shù)據(jù)

儲(chǔ)層厚度/m孔隙度/(%)原始地層壓/MPa初始含水飽和度/(%)平均滲透率/μm2152614.7300.005

表2 研究區(qū)流體數(shù)據(jù)

油密度/(g·cm-3)水密度/(g·cm-3)油黏度/(mPa·s)水黏度/(mPa·s)油壓縮系數(shù)/MPa-1水壓縮系數(shù)/MPa-10.7812.60.69.4×10-42.3×10-4

圖2 日產(chǎn)油量模擬對比圖

圖3 井底壓力預(yù)測對比圖

對比模擬結(jié)果表明,應(yīng)用低滲透模擬器模擬出的單井井底壓力、日產(chǎn)油量變化趨勢和Eclipse軟件模擬出的結(jié)果吻合度很高,表明應(yīng)用該模擬器進(jìn)行油藏?cái)?shù)值模擬具有較高的精確度，見圖2～3。

4 模型應(yīng)用

應(yīng)用該數(shù)值模擬器對新疆吐哈溫米油田米登區(qū)塊年產(chǎn)油量進(jìn)行歷史擬合,米登區(qū)塊在1993年對中侏羅統(tǒng)三間房組開始投產(chǎn),儲(chǔ)層物性以低孔低滲為主,目前尚未采取大規(guī)模注水措施,屬于低滲透天然能量衰竭式開發(fā)模式。通過收集區(qū)塊內(nèi)幾十口生產(chǎn)井的儲(chǔ)層、流體參數(shù),綜合分析得出整個(gè)區(qū)塊的擬啟動(dòng)壓力梯度約為0.053 MPa/m。應(yīng)用研發(fā)的低滲數(shù)值模擬器,采用外邊界封閉,內(nèi)邊界定壓的生產(chǎn)制度,完成溫米區(qū)塊1995年到2010年年產(chǎn)量的歷史擬合。分析歷史擬合結(jié)果表明,考慮啟動(dòng)壓力梯度影響的擬合結(jié)果與實(shí)際生產(chǎn)歷史吻合程度更高,能夠更加準(zhǔn)確地反映油田的實(shí)際生產(chǎn)情況,特別是油田生產(chǎn)中后期,誤差達(dá)到萬噸以上，見圖4。這主要是由于隨著流體的不斷采出,油田生產(chǎn)后期地層能量供給不足,對原油驅(qū)替能力下降,克服啟動(dòng)壓力的難度更大,此外油田生產(chǎn)中后期孔隙水鎖現(xiàn)象嚴(yán)重,也加大了原油驅(qū)替的難度,最終導(dǎo)致產(chǎn)量下降更加明顯,因此油田后期生產(chǎn)階段要重視啟動(dòng)壓力梯度對油田生產(chǎn)的影響,及時(shí)調(diào)整油田生產(chǎn)工作制度,加大注水開采力度,補(bǔ)充地層能量,盡可能地消除啟動(dòng)壓力梯度對油田生產(chǎn)的影響。

圖4 米登區(qū)塊年產(chǎn)量歷史擬合

5 結(jié)論

1)本文分析了低滲透油藏?cái)?shù)值模擬中啟動(dòng)壓力梯度的研究難點(diǎn),以低滲透滲流機(jī)理為研究基礎(chǔ),實(shí)現(xiàn)了對低滲透油藏滲流模型的改造,使啟動(dòng)壓力梯度的影響在數(shù)值模擬過程中得到更加充分的體現(xiàn),在此基礎(chǔ)上開發(fā)出相應(yīng)的低滲透數(shù)值模擬器。

2)通過對數(shù)值模擬器進(jìn)行精確度驗(yàn)證和實(shí)際生產(chǎn)應(yīng)用,表明啟動(dòng)壓力梯度對低滲透油藏生產(chǎn)動(dòng)態(tài)影響非常顯著,在油田生產(chǎn)中后期尤為明顯。因此在進(jìn)行生產(chǎn)動(dòng)態(tài)方案預(yù)測時(shí)需充分考慮啟動(dòng)壓力梯度的影響,并及時(shí)調(diào)整油田生產(chǎn)工作制度,加大注水開采力度,補(bǔ)充地層能量,盡可能地消除啟動(dòng)壓力梯度對油田生產(chǎn)的影響。

[1] 馮文光.非達(dá)西低速滲流的研究現(xiàn)狀與展望[J].石油勘探與開發(fā),1986,13(4):76-80

Feng Wenguang.Research Status and Prospect of Non-Darcy Low-Velocity Percolation[J].Petroleum Exploration and Development,1986,13(4):76-80

[2] 許建紅,程林松,周 穎,等.一種求解低滲透油藏啟動(dòng)壓力梯度的新方法[J].石油勘探與開發(fā),2007,34(5):594-597.

Xu Jianhong,Cheng Linsong,Zhou Ying,et al.A New Method for Calculating Kickoff Pressure Gradient in Low Permeability Reservoirs[J].Petroleum Exploration and Development,2007,34(5):594-597.

[3] 郝 斐,程林松,李春蘭,等.特低滲透油藏啟動(dòng)壓力梯度研究[J].西南石油學(xué)院學(xué)報(bào),2006,28(6):29-32.

Hao Fei,Cheng Linsong,Li Chunlan,et al.Study on Threshold Pressure Gradient in Ultra Low Permeability Reservoir[J].Journal of Southwest Petroleum Institute,2006,28(6):29-32.

[4] 侯秀林,鄧宏文,谷麗冰.一個(gè)新的描述低滲透油藏油水兩相滲流啟動(dòng)壓力梯度的公式[J].石油天然氣學(xué)報(bào),2009,31(5):342-344.

Hou Xiulin,Deng Hongwen,Guo Libing.A New Description of Low Permeability Reservoir Oil-Water Two Phase Flow of Start-Up Pressure Gradient Formula[J].Journal of and Gas Technology,2009,31(5):342-344.

[5] 彭春洋,歐陽云麗,柯立文,等.低滲透油藏滲流啟動(dòng)壓力梯度研究[J].油氣地球物理,28(6):64-66.

Peng Chunyang,OuYang Yunli,Ke Liwen,et al.The Study of Actuating Pressure Gradient of the Fluid Flow in Low Permeable Reservoir[J].Petroleum Geophysics,28(6):64-66.

[6] 張代燕,王子強(qiáng),王殿生,等.低滲透油藏最小啟動(dòng)壓力梯度實(shí)驗(yàn)研究[J].新疆地質(zhì),2011,29(1):106-109.

Zhang Daiyan,Wang Ziqiang,Wang Diansheng,et al.An Experiment Study on Measuring Minimum Starting Pressure for Low Permeability Reservoir[J].Xinjiang Geology,2011,29(1):106-109.

[7] Li Lezhong,Li Xiangfang,He Dongbo,et al.The Establishment of a New Deliverability Equation Considering Threshold Pressure Gradient[J].Engineering Sciences,2009,7(3):84-87.

[8] 田 敏,李 楠,吳文瑞.考慮啟動(dòng)壓力梯度的油水兩相滲流壓力分布[J].復(fù)雜油氣藏,2012,5(1):51-54.Tian Min,Li Nan,Wu Wenrui.Pressure Distribution of Oil-Water Seepage Considering Start-Up Pressure Gradient[J].Complex Hydrocarbon Reservoir,2012,5(1):51-54.

[9] 熊 偉,雷 群,劉先貴,等.低滲透油藏?cái)M啟動(dòng)壓力梯度[J].石油勘探與開發(fā),2009,36(2):232-235.

Xiong Wei,Lei Qun,Liu Xiangui,et al.Pseudo Threshold Pressure Gradient to Flow for Low Permeability Reservoirs[J].Petroleum Exploration and Development,2009,36(2):232-235.

[10] 徐德敏,黃潤秋,鄧英爾.低滲透軟弱巖非達(dá)西滲流擬啟動(dòng)壓力梯度試驗(yàn)研究[J].水文地質(zhì)工程地質(zhì),2008,35(3):57-60.

Xu Demin,Huang Runqiu,Deng Yinger.Non-Darcy Flow Quasi-Threshold Pressure Gradient Experimental Study for Low Permeability Soft Rock[J].Hydrogeology & Engineering Geology,2008,35(3):57-60.

[11] 李中鋒,何順利.低滲透儲(chǔ)層非達(dá)西滲流機(jī)理探討[J].特種油氣藏,2005,12(2):35-38.

Li Zhongfeng,He Shunli.Non-Darcy Percolation Mechanism in Low Permeability Reservoir[J].Special Oil and Gas Reservoir,2005,12(2):35-38.

[12] 竇宏恩,楊 旸.低滲透油藏流體滲流再認(rèn)識(shí)[J].石油勘探與開發(fā),2012,39(5):633-640.

Dou Hongen,Yang Yang.Further Understanding on Fluid Flow Through Multi-Porous Media in Low Permeability Reservoirs[J].Petroleum Exploration and Development,2012,39(5):633-640.

[13] 孔祥言.高等滲流力學(xué)[M].安徽:中國科學(xué)技術(shù)大學(xué)出版社,1999:30-57.

Kong Xiangyan.Advanced Mechanics of Fluids in Porous Media[M].Anhui:Press of University of Science and Technology of China,1999:30-57.

[14] 郭永存,盧德堂.有啟動(dòng)壓力梯度滲流的數(shù)學(xué)模型[J].中國科學(xué)技術(shù)大學(xué)學(xué)報(bào),2005,35(4):492-497.

Guo Yongcun,Lu Detang.Mathematical Model of Fluid Flow in Porous Media Related to Start-Up Pressure Gradients[J].Journal of University of Science and Technology of China,2005,35(4):492-497.

[15] 李松泉,程林松,李秀生,等.特低滲透油藏非線性滲流模型[J].石油勘探與開發(fā),2008,35(2):606-612.

Li Songquan,Cheng Linsong,Li Xiusheng,et al.Non-Linear Seepage Flow Models of Ultra-Low Permeability Reservoirs[J].Petroleum Exploration and Development,2008,35(2):606-612.

[16] 韓大匡,陳欽雷.油藏?cái)?shù)值模擬基礎(chǔ)[M].北京:石油工業(yè)出版社,1993:37-40.

Han Dakuang,Chen Qinlei.Basic of Reservoir Numerical Simulation[M].Beijing:Petroleum Industry Press,1993:37-40.

[17] 楊正明,于榮澤,蘇致新,等.特低滲透油藏非線性滲流數(shù)值模擬[J].石油勘探與開發(fā),2010,37(1):94-98.

Yang Zhengming,Yu Rongze,Su Zhixin,et al.Numerical Simulation of Nonlinear Seepage in Ultra-Low Permeability Reservoirs[J].Petroleum Exploration and Development,2010,37(1):94-98.

[18] 劉 剛,李治平,馮 彬,等.基于預(yù)處理GMRES算法的油藏?cái)?shù)值模擬研究[J].天然氣與石油,2011,29(5):44-46.

Liu Gang,Li Zhiping,Feng Bin,et al.Reservoirs Numerical Simulation Basing on Preprocessing GMRES Algorithm[J].Natural Gas and Oil,2011,29(5):44-46.

[19] 趙國忠.變啟動(dòng)壓力梯度三維三相滲流數(shù)值模擬方法[J].石油學(xué)報(bào),2006,27(增刊1):119-123.

Zhao Guozhong.Numerical Simulation of 3D and Three-Phase Flow with Variable Start-Up Pressure Gradient[J].Acta Petrolei Sinica,2006,27(Suppl 1):119-123.

[20] 張烈輝,向祖平,馮國慶.低滲氣藏考慮啟動(dòng)壓力梯度的單井?dāng)?shù)值模擬[J].天然氣工業(yè),2008,28(1):108-109.

Zhang Liehui,Xiang Zuping,Feng Guoqing.Single Well Numerical Simulation of Low-Permeability GasReservoir with Consideration of Starting Pressure Gradient[J].Natural Gas Industry,2008,28(1):108-109.

[21] 尹 虎,王新海,劉 洪,等.考慮啟動(dòng)壓力梯度的頁巖氣藏?cái)?shù)值模擬[J].天然氣與石油,2012,30(4):43-47.

Yin Hu,Wang Xinhai,Liu Hong,et al.Numerical Simulation of Shale Gas Reservoirs Under Consideration of Start-Up Pressure Gradient[J].Natural Gas and Oil,2012,30(4):43-47.

[22] 許進(jìn)進(jìn),徐旭輝,凡哲元,等.考慮啟動(dòng)壓力梯度下的火山巖氣藏?cái)?shù)值模擬研究[J].石油與天然氣地質(zhì),2010,31(2):240-243.

Xu Jinjin,Xu Xuhui,Fan Zheyuan,et al.A Numerical Simulation of Volcanic Gas Reservoirs Under Consideration of Start-Up Pressure Gradient[J].Oil & Gas Geology,2010,31(2):240-243.

2016-05-05

國家科技重大專項(xiàng)(2011 ZX 05009-006)；國家“863”計(jì)劃項(xiàng)目(SS 2013 AA 061001)

劉 剛(1985-)，男，河北保定人，工程師，碩士，主要從事低滲透油藏?cái)?shù)值模擬和非常規(guī)油氣勘探開發(fā)領(lǐng)域的研究。

10.3969/j.issn.1006-5539.2016.04.010