金忠康，方全堂，王　磊，趙　龍（.中國石化江蘇油田分公司勘探開發(fā)研究院，江蘇揚(yáng)州5009；.西南石油大學(xué)石油與天然氣工程學(xué)院，成都60500；.中國石油新疆油田分公司勘探開發(fā)研究院，新疆克拉瑪依84000）

考慮儲(chǔ)集層參數(shù)時(shí)變效應(yīng)的數(shù)值模擬方法與應(yīng)用

金忠康1，方全堂2，王磊1，趙龍3
（1.中國石化江蘇油田分公司勘探開發(fā)研究院，江蘇揚(yáng)州225009；2.西南石油大學(xué)石油與天然氣工程學(xué)院，成都610500；3.中國石油新疆油田分公司勘探開發(fā)研究院，新疆克拉瑪依834000）

受長(zhǎng)期注入水的沖刷作用，注水開發(fā)油田的儲(chǔ)集層物性及滲流特征參數(shù)發(fā)生較大變化，從而引起巖石相對(duì)滲透率變化范圍加大，變化后的相對(duì)滲透率曲線形態(tài)更為復(fù)雜。常規(guī)數(shù)值模擬軟件未能考慮儲(chǔ)集層物性和相對(duì)滲透率參數(shù)的變化，導(dǎo)致模擬結(jié)果與實(shí)際情況不完全符合。通過建立蘇北盆地東部低滲透油藏沙20斷塊的滲流機(jī)理模型，評(píng)價(jià)了長(zhǎng)期注水開發(fā)過程中油藏典型特征參數(shù)的變化規(guī)律，實(shí)現(xiàn)對(duì)儲(chǔ)集層物性和相對(duì)滲透率參數(shù)時(shí)變數(shù)值模擬方法的表征。結(jié)果表明，針對(duì)中低滲透油藏，長(zhǎng)期注水沖刷作用導(dǎo)致儲(chǔ)集層物性變差，平均剩余油飽和度計(jì)算值有所增大，尤其是在注水井附近，但考慮時(shí)變效應(yīng)前后剩余油分布模式基本一致。

注水開發(fā)；儲(chǔ)集層；時(shí)變效應(yīng)；數(shù)值模擬；剩余油分布

經(jīng)過長(zhǎng)期注入水的沖刷作用，高含水期油藏儲(chǔ)集層物性會(huì)發(fā)生變化［1-6］，但通常的油藏?cái)?shù)值模擬研究是從油藏投入開發(fā)開始，儲(chǔ)集層物性參數(shù)均被視為定值，從而導(dǎo)致模擬結(jié)果與實(shí)際情形有一定出入。目前，儲(chǔ)集層物性參數(shù)的時(shí)變數(shù)值模擬研究主要是通過修改傳統(tǒng)黑油模型，編制相應(yīng)程序，實(shí)現(xiàn)物性參數(shù)時(shí)變模擬［7-12］，這種方法能夠靈活地考慮參數(shù)時(shí)變效應(yīng)，但商用化難度大，難以滿足現(xiàn)場(chǎng)快速應(yīng)用的要求。筆者采用相應(yīng)的數(shù)學(xué)處理手段，將前人建立的滲透率變化倍數(shù)與注水沖刷倍數(shù)的關(guān)系進(jìn)行了修改，推導(dǎo)出滲透率變化倍數(shù)與網(wǎng)格節(jié)點(diǎn)含水飽和度的變化關(guān)系，以等效實(shí)現(xiàn)注水開發(fā)過程中儲(chǔ)集層滲透率的時(shí)變數(shù)值模擬，并將這一成果應(yīng)用于蘇北盆地東部低滲透油藏沙20斷塊，分析對(duì)比考慮儲(chǔ)集層滲透率時(shí)變效應(yīng)前后的井底流壓、剩余油分布以及時(shí)變效應(yīng)對(duì)儲(chǔ)集層的影響，結(jié)果表明，這一方法可用于評(píng)價(jià)長(zhǎng)期注水開發(fā)過程中油藏儲(chǔ)集層物性及滲流特征參數(shù)。

1　相對(duì)滲透率時(shí)變規(guī)律

文獻(xiàn)［13］根據(jù)大量礦場(chǎng)不同開發(fā)階段的取巖心井物性分析資料與室內(nèi)水沖刷實(shí)驗(yàn)成果，研究油田注水開發(fā)過程中滲透率變化規(guī)律，得出了勝利油區(qū)主力砂巖油藏相關(guān)數(shù)學(xué)關(guān)系。

（1）滲透率變化倍數(shù)與注水沖刷倍數(shù)的關(guān)系

式中MK——滲透率變化倍數(shù)，注水沖刷后滲透率與注水沖刷前滲透率之比；

R——注水沖刷倍數(shù)。根據(jù)注水沖刷程度的不同，可重新建立與之相應(yīng)的滲透率參數(shù)場(chǎng)。

（2）相對(duì)滲透率曲線端點(diǎn)變化量與滲透率變化倍數(shù)的關(guān)系

式中

ΔSwc——束縛水飽和度變化量，無因次；

ΔSor——?dú)堄嘤惋柡投茸兓?，無因次。

在實(shí)際油藏開發(fā)過程中，平均注水沖刷倍數(shù)遠(yuǎn)小于（1）式中的注入倍數(shù)臨界值1 210，但由于油井投產(chǎn)時(shí)間的不同及油水滲流場(chǎng)的變化，在主流線方向上，注水沖刷程度明顯強(qiáng)于其他區(qū)域，即油藏存在注水沖刷程度的非均質(zhì)性。為了解決注水沖刷程度非均質(zhì)性描述的難題，本文借助于典型油藏機(jī)理模型，建立注水沖刷倍數(shù)與含水飽和度變化倍數(shù)的關(guān)系式，以實(shí)現(xiàn)油藏?cái)?shù)值模擬中儲(chǔ)集層時(shí)變參數(shù)的表征。

2　時(shí)變模擬方法的建立

2.1機(jī)理模型

利用商用油藏?cái)?shù)值模擬器開展參數(shù)時(shí)變研究，在無法直接獲取單一網(wǎng)格節(jié)點(diǎn)注水沖刷量的情況下，可以通過建立以典型油藏參數(shù)為基準(zhǔn)的一維流動(dòng)模型，模擬不同注水量的含水飽和度變化，建立滲透率變化倍數(shù)與網(wǎng)格節(jié)點(diǎn)含水飽和度變化之間的關(guān)系式，以滿足在油藏?cái)?shù)值模擬器中時(shí)變效應(yīng)的表征。

在典型油藏模型中抽取一維流動(dòng)網(wǎng)格（圖1），兩端分別部署注水井（沙20-13井）和采油井（沙20-10井），設(shè)定初期采油速度不超過5%，且保持注采平衡。以油藏綜合含水率接近100%為模擬終止條件，即可獲得不同沖刷強(qiáng)度下的油藏滲透率和相對(duì)滲透率特征參數(shù)。

選定一維流動(dòng)模型某一網(wǎng)格節(jié)點(diǎn)為研究對(duì)象，假設(shè)注水沖刷使網(wǎng)格孔隙體積全部充滿，統(tǒng)計(jì)不同模擬時(shí)間下此網(wǎng)格的含水飽和度以及流經(jīng)此網(wǎng)格節(jié)點(diǎn)的總注水體積和前緣孔隙存水量，其中，水驅(qū)前緣網(wǎng)格節(jié)點(diǎn)的注水沖刷量為總注水體積與前緣孔隙存水量之差。以水驅(qū)前緣的含水飽和度變化和通過網(wǎng)格節(jié)點(diǎn)的注水沖刷量為描述對(duì)象，建立前緣含水飽和度變化量與注水沖刷倍數(shù)的關(guān)系（圖2）。

圖1　一維流動(dòng)機(jī)理模型

圖2　注水沖刷倍數(shù)與網(wǎng)格節(jié)點(diǎn)含水飽和度變化量關(guān)系

2.2基于含水飽和度場(chǎng)變化的參數(shù)時(shí)變模型

2.2.1滲透率的時(shí)變

進(jìn)一步將圖2中回歸的公式轉(zhuǎn)換為以含水飽和度變化量為自變量的公式，即將原有的以注水沖刷倍數(shù)為界限的公式變?yōu)橐院柡投茸兓喀膕w為劃分依據(jù)的公式:

故可將（1）式變換為

2.2.2時(shí)變數(shù)值模擬的實(shí)現(xiàn)

通過建立適用于目標(biāo)區(qū)塊的注水沖刷倍數(shù)與含水飽和度變化量回歸關(guān)系式，將難以描述的注水沖刷倍數(shù)轉(zhuǎn)換為不同節(jié)點(diǎn)處的含水飽和度變化量，進(jìn)而通過在油藏?cái)?shù)值模擬過程中，計(jì)算出各點(diǎn)的含水飽和度變化量，結(jié)合對(duì)應(yīng)的滲透率以及相對(duì)滲透率曲線的時(shí)變公式，即可在常規(guī)的油藏?cái)?shù)值模擬過程中描述時(shí)變現(xiàn)象，其中滲透率時(shí)變和相對(duì)滲透率曲線的時(shí)變?cè)谀壳俺Ｓ玫腅CLISPE軟件中主要采用關(guān)鍵字修改的方式來實(shí)現(xiàn)。

（1）滲透率時(shí)變的實(shí)現(xiàn)研究某一時(shí)間節(jié)點(diǎn)發(fā)生時(shí)變的滲透率場(chǎng)時(shí)，通過計(jì)算此時(shí)間節(jié)點(diǎn)的含水飽和度變化量，根據(jù)（6）式計(jì)算出不同含水飽和度變化量下的滲透率變化倍數(shù)，并確定滲透率變化倍數(shù)分布模型，將其與初始滲透率參數(shù)場(chǎng)在PETREL軟件中做疊加計(jì)算，即可得到對(duì)應(yīng)時(shí)間節(jié)點(diǎn)發(fā)生時(shí)變后新的滲透率場(chǎng)。為進(jìn)一步研究此時(shí)間節(jié)點(diǎn)滲透率變化，需要通過PETREL軟件將滲透率變化倍數(shù)場(chǎng)導(dǎo)出，將其轉(zhuǎn)化為ECLIPSE軟件中的傳導(dǎo)率倍乘因子，并通過在ECLIPSE軟件Schedule文件中發(fā)生滲透變化的時(shí)間節(jié)點(diǎn)處添加關(guān)鍵字MULTX-和MULTY-來修改網(wǎng)格傳導(dǎo)率，從而實(shí)現(xiàn)ECLIPSE軟件中對(duì)滲透率的修改。

（2）相對(duì)滲透率時(shí)變的實(shí)現(xiàn)首先要通過（6）式計(jì)算出不同含水飽和度變化量下的滲透率變化倍數(shù)。然后根據(jù)（2）式和（3）式計(jì)算出不同滲透率變化倍數(shù)下的相對(duì)滲透率端點(diǎn)值的變化，生成發(fā)生時(shí)變后多條不同的相對(duì)滲透率曲線。最后按照滲透率變化倍數(shù)與變化后的相對(duì)滲透率曲線編號(hào)對(duì)應(yīng)關(guān)系，利用PETREL軟件分別計(jì)算出每一個(gè)網(wǎng)格節(jié)點(diǎn)所對(duì)應(yīng)的相對(duì)滲透率曲線編號(hào)，將此文件導(dǎo)出后轉(zhuǎn)化為ECLIPSE軟件中的Satnum文件，即可實(shí)現(xiàn)相對(duì)滲透率的時(shí)變。

2.3參數(shù)時(shí)變效應(yīng)分析

根據(jù)上述方法得出考慮滲透率時(shí)變效應(yīng)前后的井底流壓對(duì)比圖（圖3）以及剩余油飽和度分布對(duì)比圖（圖4），這里分別選擇1口采油井和1口注水井作為典型特例進(jìn)行分析。

圖3　采油井沙20-10井和注水井沙20-13井井底流壓對(duì)比

由圖3可以看出，隨著注水時(shí)間的增加，考慮儲(chǔ)集層滲透率時(shí)變效應(yīng)后的注水井井底流壓高于未考慮儲(chǔ)集層滲透率時(shí)變效應(yīng)的井底流壓，而采油井井底流壓低于未考慮滲透率時(shí)變效應(yīng)的井底流壓。說明針對(duì)低滲透油藏，長(zhǎng)期注水開發(fā)，可能使得儲(chǔ)集層物性變差，導(dǎo)致采油井生產(chǎn)能力及注水井注入能力下降，這與沙20斷塊巖心室內(nèi)水沖刷實(shí)驗(yàn)后巖心滲透率有所降低的結(jié)論一致。從剩余油飽和度分布分析，長(zhǎng)期注水開發(fā)引起的儲(chǔ)集層滲透率時(shí)變效應(yīng)后，全區(qū)平均剩余油飽和度計(jì)算值有所增加（圖4），儲(chǔ)集層滲透率時(shí)變效應(yīng)前、后計(jì)算剩余飽和度平均值分別為0.371 和0.379.但無論是否考慮滲透率時(shí)變效應(yīng)，其剩余油分布模式基本一致，剩余油主要分布在注采井間非優(yōu)勢(shì)通道區(qū)域，優(yōu)勢(shì)通道區(qū)域剩余油飽和度較低。

圖4　考慮滲透率時(shí)變效應(yīng)前（a）后（b）剩余油飽和度分布

注水開發(fā)過程中，由于注采井間優(yōu)勢(shì)通道區(qū)域注入水沖刷作用大，滲透率變化大；而非優(yōu)勢(shì)通道區(qū)域注入水沖刷作用小，滲透率變化小。儲(chǔ)集層滲透率的增大，一方面增加了流體的流動(dòng)能力，增大了水驅(qū)驅(qū)油效率；另一方面，加劇了儲(chǔ)集層平面非均質(zhì)性，降低了水驅(qū)波及系數(shù)。兩個(gè)方面共同作用、相互矛盾，但由于每口井的控制面積相對(duì)固定，當(dāng)生產(chǎn)至極限含水率，波及系數(shù)對(duì)水驅(qū)開發(fā)效果的影響較大［14-16］，最終導(dǎo)致水驅(qū)開發(fā)效果變差。

3　結(jié)論

（1）以典型油藏參數(shù)為基礎(chǔ)建立一維水驅(qū)機(jī)理模型，分析了注水沖刷過程中含水飽和度與注水沖刷倍數(shù)的相關(guān)性，并以水驅(qū)前緣含水飽和度的變化為描述對(duì)象，建立注水沖刷倍數(shù)與含水飽和度變化的函數(shù)關(guān)系是考慮儲(chǔ)集層參數(shù)時(shí)變效應(yīng)的基礎(chǔ)。

（2）綜合分析典型油藏注水開發(fā)過程中注水沖刷倍數(shù)非均質(zhì)性特征，將注水沖刷倍數(shù)與參數(shù)時(shí)變函數(shù)關(guān)系轉(zhuǎn)換為含水飽和度變化量與參數(shù)時(shí)變的函數(shù)關(guān)系，可以實(shí)現(xiàn)在商用數(shù)值模擬器內(nèi)分時(shí)段進(jìn)行油藏參數(shù)時(shí)變的模擬研究。

（3）針對(duì)長(zhǎng)期注水開發(fā)的低滲透油田，在不考慮油藏參數(shù)時(shí)變情況下的油藏?cái)?shù)值模擬得到的平均剩余油飽和度低于油藏實(shí)際值，影響較為明顯的是注水井附近，但其剩余油分布模式與是否考慮時(shí)變效應(yīng)基本無關(guān)。

［1］劉天鶴，郭莉，左毅.高含水油田注水后油藏參數(shù)變化規(guī)律研究［J］.特種油氣藏，2005，12（2）:42-44.

LIU Tianhe，GUO Li，ZUO Yi.The research on the parameters change rule of high water reservoir after water injection［J］.Special Oil and Gas Reservoirs，2005，12（2）:42-44.

［2］王洪光，蔣明，張繼春，等.高含水期油藏儲(chǔ)集層物性變化特征模擬研究［J］.石油學(xué)報(bào)，2004，25（6）:53-58.

WANG Hongguang，JIANG Ming，ZHANG Jichun，et al.Simulation on variation of physical properties in high water-cut reservoir［J］. Acta Petrolei Sinica，2004，25（6）:53-58.

［3］胡治華，馬奎前，劉宗賓，等.海上S油田注水開發(fā)后期儲(chǔ)層物性變化規(guī)律及應(yīng)用［J］.科學(xué)技術(shù)與工程，2014，14（15）:164-168.

HU Zhihua，MA Kuiqian，LIU Zongbin，et al.Research and application of the properties variation regularity of water injection reservoirs in Soffshore oilfield［J］.Science Technology and Engineering，2014，14（15）:164-168.

［4］陳小凡，劉峰，王博，等.長(zhǎng)期沖刷條件下的滲透率變化和微觀模型研究［J］.特種油氣藏，2012，19（1）:66-69.

CHEN Xiaofan，LIU Feng，WANG Bo，et al.Research on permeability changes and microscopic model under the condition of long-term washou［tJ］.Special Oil and Gas Reservoirs，2012，19（1）:66-69.

［5］吳素英.長(zhǎng)期注水沖刷儲(chǔ)層參數(shù)變化規(guī)律及對(duì)開發(fā)效果的影響［J］.大慶石油地質(zhì)與開發(fā)，2006，25（4）:35-37.

WU Suying.Variation rule of oil layer parameters washed by longterm injected water and its impact on development effect［J］.Petroleum Geology&Oilfield Development in Daqing，2006，25（4）: 35-37.

［6］ZHANGHongxing，LIUQingnan，LIFengqin.Variationsofpetrophysical parameters after sandstone reservoirs watered out in Daqing oilfield［R］.SPE 30884，1995.

［7］HAUGEN A，F(xiàn)ERNO M，GRAUE A.Numerical sinulation and sensityvity analysis of in-situ fluid flow in MRI laboratory waterfloods of fracturde carbonate rocks at different wettabilities［R］.SPE 116145，2008.

［8］姚風(fēng)英，陳曉軍，宗習(xí)武，等.勝坨油田注水開發(fā)過程中的潤(rùn)濕性變化［J］.油氣地質(zhì)與采收率，2002，9（4）:58-60.

YAO Fengying，CHEN Xiaojun，ZONGXiwu，etal.Wettability change during the waterflooding process in Shengtuo oilfield［J］.Petroleum Geology and Recovery Efficiency，2002，9（4）:58-60.

［9］趙壽元.注水開發(fā)儲(chǔ)層物性參數(shù)變化數(shù)值模擬研究［D］.青島:中國石油大學(xué)，2009.

ZHAO Shouyuan.Research of numerical simulation on the change of reservoir parameters during flood development［D］.Qingdao:China University of Petroleum，2009.

［10］李曉燕.注水開發(fā)后儲(chǔ)層物理特征參數(shù)變化的數(shù)模研究［D］.北京:中國石油大學(xué)，2007.

LI Xiaoyan.Research of numerical simulation on the change of reservoir parameters after flood development［D］.Beijing:China University of Petroleum，2007.

［11］崔傳智，趙曉燕.考慮油藏儲(chǔ)層參數(shù)變化的數(shù)值模擬研究［J］.水動(dòng)力學(xué)研究與進(jìn)展，2004，19（增刊）:912-915.

CUI Chuanzhi，ZHAO Xiaoyan.The reservoir numerical simulation study with the variety of reservoir parameters［J］.Journal of Hydrodynamics，2004，19（Supp.）:912-915.

［12］蓋英杰，呂德靈，郭元靈.高含水期油藏分階段數(shù)值模擬［J］.油氣采收率技術(shù)，2000，7（1）:54-56.

GAI Yingjie，LV Deling，GUO Yuanling.Numerical simulation by stages about the reservoir at high water cut period［J］.Oil & Gas Recovery Technology，2000，7（1）:54-56.

［13］劉顯太.中高滲透砂巖油藏儲(chǔ)層物性時(shí)變數(shù)值模擬技術(shù)［J］.油氣地質(zhì)與采收率，2012，19（5）:58-62.

LIU Xiantai.Study on numerical simulation technology based on time varying physical properties in mid-high permeability sandstone reservoirs［J］.Petroleum Geology and Recovery Efficiency，2012，19（5）:58-62.

［14］黃炳光，劉蜀知.實(shí)用油藏工程與動(dòng)態(tài)分析方法［M］.北京:石油工業(yè)出版社，1997:77-79.

HUANG Bingguang，LIU Shuzhi.Practical and dynamic analysis method of reservoir engineering［M］.Beijing:Petroleum Industry Press，1997:77-79.

［15］陳元千，王小林，姚尚林，等.加密井提高注水開發(fā)油田采收率的評(píng)價(jià)方法［J］.新疆石油地質(zhì)，2009，30（6）:706-709.

CHEN Yuanqian，WANG Xiaolin，YAO Shanglin，et al.A method for evaluation of EOR of waterflooding oilfield by infill well process ［J］.Xinjiang Petroleum Geology，2009，30（6）:706-709.

［16］王文濤.江蘇油田復(fù)雜斷塊油藏水驅(qū)開發(fā)效果評(píng)價(jià)［D］.成都:西南石油大學(xué)，2012.

WANG Wentao.Water flooding effect evaluation of the Jiangsu oilfield complex faulted block reservoirs［D］.Chengdu:Southwest Petroleum University，2012.

（編輯顧新元）

Numerical Simulation Method and Its Application when Considering Time-Dependent Effect of Reservoir Parameters

JIN Zhongkang1，F(xiàn)ANG Quantang2，WANG Lei1，ZHAO Long3
（1.Research Institute of Exploration and Development，Jiangsu Oilfield Company，Sinopec，Yangzhou，Zhejiang 225009，China；2.School of Petroleum Engineering，Southwest Petroleum University，Chengdu，Sichuan 610500，China；
3.Research Institute of Exploration and Development，Xinjiang Oilfield Company，Petrochina，Karamay，Xinjiang834000，China）

Due to the influence of long-term scouring of injected-water，the physical properties and characteristic percolation parameters of reservoirs tend to change greatly in water-injection oilfields，which will result in large variations of relative permeability of rocks and more complexity of relative permeability curves after the changes.Conventional numerical modeling software never considers the variations of physical properties and relative permeabilities of reservoirs，so that the modeling results can't match with the actual situation completely. A percolation mechanism model has been established for low-permeability Sha-20 fault block in eastern Subei basin to evaluate the changes of typical characteristic reservoir parameters during long-term water injection，by which the characterization of time-dependent numerical simulation method can be realized for physical properties and relative permeability of reservoirs.The result shows that long-term scouring of injected water will result in poorer physical properties and enlargement of the calculated average remaining oil saturation regarding medium and low permeability reservoirs，especially in the surroundings of injectors.But the remaining oil distribution patterns are basically the same before and after considering the time-depending effect.

water flooding development；reservoir；time-dependent effect；numerical simulation；remaining oil distribution

TE311

A

1001-3873（2016）03-0342-04

10.7657/XJPG20160319

2015-12-02

2016-02-25

金忠康（1970-），男，江蘇吳江人，高級(jí)工程師，油藏工程，（Tel）15895708695（E-mail）jinzk.jsyt@sinopec.com