吳克柳，李相方，樊兆琪，李武廣，李元生，羊新州

(中國(guó)石油大學(xué)石油工程教育部重點(diǎn)實(shí)驗(yàn)室，北京102249)

油水相對(duì)滲透率的準(zhǔn)確獲取對(duì)油藏開(kāi)發(fā)方案制定、開(kāi)發(fā)效果預(yù)測(cè)及油藏管理具有重要意義［1-9］。低滲-特低滲儲(chǔ)層由于孔喉細(xì)小、孔隙結(jié)構(gòu)復(fù)雜及孔隙內(nèi)表面積大，具有很大的固液表面分子作用力和毛管力，其流體滲流顯著差異于常規(guī)中高滲儲(chǔ)層中的流體，表現(xiàn)出很強(qiáng)的非達(dá)西滲流特征［10-15］。與常規(guī)中高滲油藏油水相對(duì)滲透率相比，啟動(dòng)壓力梯度及毛管力對(duì)低滲-特低滲油藏油水滲流影響顯著，且當(dāng)儲(chǔ)層傾斜度較大時(shí)，重力效應(yīng)的影響也不可忽略［16-18］。筆者結(jié)合低滲-特低滲油藏油水滲流特征，考慮啟動(dòng)壓力梯度、重力及毛管力的影響，建立低滲-特低滲油藏非穩(wěn)態(tài)油水相對(duì)滲透率計(jì)算模型，為該類油藏的有效開(kāi)發(fā)提供理論支撐。

1 模型建立

假設(shè):儲(chǔ)層為均勻多孔介質(zhì);驅(qū)動(dòng)力不變，為水驅(qū);油水性質(zhì)保持不變;油水不發(fā)生反應(yīng)，且無(wú)相間傳質(zhì)現(xiàn)象;忽略儲(chǔ)層及流體的壓縮性。

1.1 運(yùn)動(dòng)方程

考慮低滲-特低滲儲(chǔ)層多孔介質(zhì)流體滲流特性，則油、水相非達(dá)西流運(yùn)動(dòng)方程分別為

式中，vo和vw分別為油相和水相滲流速度，m/s;k為多孔介質(zhì)滲透率，m2;Kro和Krw分別為油相和水相相對(duì)滲透率;μo和μw分別為油相和水相黏度，Pa·s;po和pw分別為油相和水相壓力，MPa;x為流動(dòng)距離，m;λo和λw分別為油相和水相啟動(dòng)壓力梯度，MPa/m;ρo和 ρw分別為油相和水相密度，g/cm3;g為重力加速度，m/s2;a為流體在儲(chǔ)層多孔介質(zhì)中流動(dòng)流線與水平面的夾角，(°)。

根據(jù)毛管力的定義有

式中，pc為毛管壓力，是含水飽和度 Sw的函數(shù)，MPa。

結(jié)合式(3)則式(1)可變形為

儲(chǔ)層多孔介質(zhì)中流體總滲流速度為

則油、水分流量可分別表示為

式中，v為總滲流速度，m/s;fo和fw分別為油相和水相分流量。

由式(1)和式(2)變形得

式(8)減式(7)，聯(lián)立式(6)得

物質(zhì)平衡關(guān)系有fo=1-fw，則式(9)變?yōu)?/p>

由式(10)變形得

式中，t為流動(dòng)時(shí)間，s。

1.2 Buckley-Leverett方程

忽略油、水壓縮性，一維均質(zhì)地層水驅(qū)油過(guò)程中的油、水相連續(xù)性方程分別為

式中，φ為多孔介質(zhì)孔隙度;So為含油飽和度。

結(jié)合式(6)，式(13)可變形為

由式(14)變形整理，可得等含水飽和度面在多孔介質(zhì)中的移動(dòng)速度為

這就是Buckley-Leverett方程，油水兩相流體不可壓縮時(shí)才適用。其中含水率fw用式(11)代入，則式(15)即為考慮了啟動(dòng)壓力梯度、重力及毛管力影響的Buckley-Leverett方程。

1.3 油水相對(duì)滲透率計(jì)算方程

巖心兩端壓差Δp與相對(duì)滲透率的關(guān)系通過(guò)式(2)變形可表達(dá)為

將式(6)代入式(16)得

巖石多孔介質(zhì)多為水濕，故巖心兩端壓差用含水相的參數(shù)表達(dá)較為合理，為

式中，L為巖心長(zhǎng)度，m。

依據(jù)等含水飽和度面推進(jìn)速度，由式(15)可推導(dǎo)出

其中

式中，f′w2為在巖心末端的分流量對(duì)含水飽和度的導(dǎo)數(shù);f′w為分流量對(duì)含水飽和度的導(dǎo)數(shù);為累積注入孔隙體積倍數(shù);QIw(t)為累積注入水量，m3;A為巖心橫截面積，m2。

將式(17)和(19)代入式(18)得

將式(19)代入式(20)并兩端求導(dǎo)，整理得水相相對(duì)滲透率為

聯(lián)立式(21)和(11)得油相相對(duì)滲透率表達(dá)式為

式中，下標(biāo)“2”代表巖心末端處。

從水油兩相相對(duì)滲透率式(21)和(22)可知:①如果不考慮啟動(dòng)壓力梯度、重力和毛管壓力，式(21)和(22)退化為常規(guī)的JBN方法油水相對(duì)滲透率計(jì)算公式;②毛管力僅對(duì)油相相對(duì)滲透率有影響，即對(duì)潤(rùn)濕相(水相)無(wú)作用，這一現(xiàn)象與現(xiàn)場(chǎng)情況相符;③油水相對(duì)滲透率均受啟動(dòng)壓力梯度和重力的影響。

1.4 含水飽和度及梯度計(jì)算

由式(21)和(22)可知，油水相對(duì)滲透率計(jì)算首先必須求取巖心末端的含水飽和度及其梯度。由物質(zhì)平衡原理有

式中，Swa為平均含水飽和度;Swc為束縛水飽和度;∑Qo為累積產(chǎn)油量，m3。

巖心末端的含水飽和度可表示為

式中，Sw2為巖心末端處含水飽和度;fo2為巖心末端處油相分流量。

由偏導(dǎo)數(shù)基本性質(zhì)可得

結(jié)合式(15)，式(25)變形為

式中，Q為累積總產(chǎn)液量，其與累積注入水量QIw(t)幾乎相等，m3。

根據(jù)物質(zhì)平衡關(guān)系，巖心中含水飽和度變化的比例等于累積產(chǎn)液量變化的比例［16］，即

將式(27)變形整理得

式中，Sw(x，t)為t時(shí)刻巖心x處含水飽和度(x，t)為t時(shí)刻巖心［0，x］區(qū)間內(nèi)的平均含水飽和度。

對(duì)式(27)兩邊求導(dǎo)可得

巖心中水驅(qū)油時(shí)，有

對(duì)式(30)連續(xù)兩次求導(dǎo)并代入式(29)得

將式(31)代入式(26)，則巖心末端處含水飽和度梯度為

聯(lián)立式(21)、(22)、(24)和(32)就可計(jì)算出低滲-特低滲儲(chǔ)層多孔介質(zhì)中油水滲流時(shí)考慮啟動(dòng)壓力梯度、重力及毛管力影響的油水相對(duì)滲透率曲線。

2 非穩(wěn)態(tài)試驗(yàn)

為了計(jì)算低滲-特低滲油藏油水相對(duì)滲透率，分析啟動(dòng)壓力梯度、重力及毛管力對(duì)相對(duì)滲透率的影響，需通過(guò)試驗(yàn)獲取計(jì)算相對(duì)滲透率的基礎(chǔ)參數(shù)。試驗(yàn)裝置見(jiàn)文獻(xiàn)［20］。

2.1 試驗(yàn)樣品

原油:密度為0.809 g/cm3，黏度為0.66 mPa·s，體積系數(shù)為1.045，原始?xì)庥捅葹?.23 m3/t。

地層水:根據(jù)地層水礦化度及離子含量配置試驗(yàn)用地層水，總礦化度為158.842 g/L，pH值為5.64，水型為 CaCl2型，密度為 1.01 g/cm3，黏度為0.93 mPa·s。

巖心:試驗(yàn)所用巖心取自延長(zhǎng)油田長(zhǎng)6儲(chǔ)層，巖心直徑為2.37 cm，長(zhǎng)度為4.79 cm，氣測(cè)孔隙度為10.8%，氣測(cè)滲透率為2.35×10-3μm2。

2.2 試驗(yàn)條件

采用非穩(wěn)態(tài)法測(cè)定油水相對(duì)滲透率曲線。試驗(yàn)溫度為地層溫度(32℃)，試驗(yàn)壓力模擬原始地層條件下有效覆壓 7.5 MPa［19］，驅(qū)替速度為 0.6 mL/min。

2.3 試驗(yàn)步驟

(1)對(duì)巖心進(jìn)行清洗、烘干，測(cè)定基本參數(shù);

(2)巖心抽真空飽和地層水，裝入巖心夾持器;

(3)調(diào)整出口處的油、水體積計(jì)量裝置;

(4)測(cè)定巖心的含油飽和度和束縛水飽和度，用油驅(qū)水的方法建立束縛水，直至不出水為止，記錄驅(qū)出的水量，計(jì)算巖心的含油飽和度和束縛水飽和度;

(5)在試驗(yàn)條件下進(jìn)行水驅(qū)油，驅(qū)替至含水率達(dá)到或接近99%，記錄各個(gè)時(shí)刻的驅(qū)替壓力、產(chǎn)油量和產(chǎn)水量。

2.4 試驗(yàn)結(jié)果分析

2.4.1 油水相對(duì)滲透率曲線

把試驗(yàn)測(cè)得基礎(chǔ)數(shù)據(jù)等代入式(21)、(22)、(24)和(32)，其中水相啟動(dòng)壓力梯度為0.013 MPa/m，油相啟動(dòng)壓力梯度為0.016 MPa/m，巖心傾斜度為10°，計(jì)算巖心油水兩相滲流時(shí)考慮啟動(dòng)壓力梯度、重力及毛管力影響的油水相對(duì)滲透率曲線，結(jié)果見(jiàn)圖1。

圖1 非穩(wěn)態(tài)油水相對(duì)滲透率曲線Fig.1 Unsteady-state oil-water relative permeability curves

由圖1可以看出:考慮啟動(dòng)壓力梯度時(shí)，油水相對(duì)滲透率均下降，但油相對(duì)滲透率下降幅度較小，等滲點(diǎn)向右下方移動(dòng);考慮重力時(shí)，水相對(duì)滲透率下降，油相對(duì)滲透率升高，油水相對(duì)滲透率等滲點(diǎn)幾乎不變;考慮毛管力時(shí)，油相對(duì)滲透率升高，水相對(duì)滲透率不變，油水相對(duì)滲透率等滲點(diǎn)向左上方移動(dòng);隨著含水飽和度的增加，啟動(dòng)壓力梯度、重力及毛管力對(duì)油水相對(duì)滲透率的影響逐漸增強(qiáng)。

啟動(dòng)壓力梯度阻礙油水相的滲流，水相對(duì)滲透率下降。油相啟動(dòng)壓力梯度一般大于水相啟動(dòng)壓力梯度，油相滲流能力比水相下降更明顯，所以油相對(duì)滲透率也下降;對(duì)于傾斜儲(chǔ)層，生產(chǎn)井在高部位時(shí)重力作為油水相滲流阻力，水相對(duì)滲透率下降，油相對(duì)滲透率提高，等滲點(diǎn)幾乎不變;巖心水濕，毛管力對(duì)水相無(wú)影響，對(duì)油相而言，毛管力作為動(dòng)力提高了油相滲流能力，油相對(duì)滲透率增大。

2.4.2 敏感性分析

為了分析毛管力對(duì)油水相對(duì)滲透率影響程度，分別選取了4條毛管力曲線(圖2)，計(jì)算不同毛管力下的油水相對(duì)滲透率曲線。

圖2 毛管力曲線Fig.2 Capillary pressure curves

分析不同啟動(dòng)壓力梯度下重力和毛管力對(duì)油水相對(duì)滲透率的影響，結(jié)果如圖3所示。

圖3表明:啟動(dòng)壓力梯度越大，油水相對(duì)滲透率下降越明顯;隨著含水飽和度的增加，啟動(dòng)壓力梯度對(duì)油水相對(duì)滲透率的影響逐漸增強(qiáng);儲(chǔ)層傾斜度越大，油水相對(duì)滲透率變化越明顯;隨著含油飽和度的增加，油相對(duì)滲透率也增加，且受毛管力的影響越大，即在低含水飽和度時(shí)毛管力大，作為動(dòng)力更有助于油相滲流。

圖3 油水相對(duì)滲透率敏感性分析Fig.3 Sensitivity analysis on oil-water relative permeability

3 結(jié)論

(1)考慮低滲-特低滲儲(chǔ)層多孔介質(zhì)中油水滲流特征建立并求解了低滲-特低滲油藏非穩(wěn)態(tài)油水相對(duì)滲透率計(jì)算數(shù)學(xué)模型。

(2)油水啟動(dòng)壓力梯度越大，油水相對(duì)滲透率下降越明顯;隨著含水飽和度的增加，啟動(dòng)壓力梯度對(duì)油水相對(duì)滲透率的影響逐漸增強(qiáng);儲(chǔ)層傾斜度越大，油水相對(duì)滲透率變化越明顯;隨著含水飽和度的增加，油相對(duì)滲透率增大且受毛管力的影響越大，即在低含水飽和度時(shí)毛管力大，作為動(dòng)力更有助于油相滲流。

(3)啟動(dòng)壓力梯度對(duì)油水相對(duì)滲透率影響最為顯著，其次是重力，而毛管力僅對(duì)油相對(duì)滲透率有影響，對(duì)水相對(duì)滲透率無(wú)影響。

［1］ 張繼成，宋考平.相對(duì)滲透率特征曲線及其應(yīng)用［J］.石油學(xué)報(bào)，2007，28(4):104-107.

ZHANG Ji-cheng，SONG Kao-ping.Eigen curve of relative permeability and its application ［J］.Acta Petrolei Sinica，2007，28(4):104-107.

［2］ 王曙光，趙國(guó)忠，余碧君.大慶油田油水相對(duì)滲透率統(tǒng)計(jì)規(guī)律及其應(yīng)用［J］.石油學(xué)報(bào)，2005，26(3):78-81.

WANG Shu-guang，ZHAO Guo-zhong，YU Bi-jun.Statistical regularity of oil-water relative permeability in Daqing Oilfield［J］.Acta Petrolei Sinica，2005，26(3):78-81.

［3］ 趙培華，譚廷棟，宋社民，等.利用生產(chǎn)測(cè)井資料確定水驅(qū)油藏產(chǎn)層剩余油飽和度方法［J］.石油學(xué)報(bào)，1997，18(2):54-60.

ZHAO Pei-hua，TAN Ting-dong，SONG She-min，et al.Determination of remaining oil saturation of water-drive pool from production logging data［J］.Acta Petrolei Sinica，1997，18(2):54-60.

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

WANG Hong-guang，JIANG Ming，ZHANG Ji-chun，et al.Simulation on variation of physical properties in high water-cut reservoir［J］.Acta Petrolei Sinica，2004，25(6):53-58.

［5］ JONES S C，ROSZELLE W O.Graphical techniques for determining relative permeability from displacement experiments［J］.JPT，1978，30(5):807-817.

［6］ SIGMUND P M，MCCAFFERY F G.An improved unsteady-state procedure for determining the relative permeability characteristics of heterogeneous porous media［J］.SPEJ，1979，19(1):593-608.

［7］ 桓冠仁，沈平平.一種非穩(wěn)態(tài)油水相對(duì)滲透率曲線計(jì)算方法［J］.石油勘探與開(kāi)發(fā)，1982，9(2):52-58.

HUAN Guan-ren，SHEN Ping-ping.A new unsteadystate method for calculating relative permeability［J］.Petroleum Exploration and Development，1982，9(2):52-58.

［8］ 楊小平.精確計(jì)算相對(duì)滲透率的方法［J］.石油勘探與開(kāi)發(fā)，1998，25(6):63-66.

YANG Xiao-ping.A method to calculate the relative permeability accurately［J］.Petroleum Exploration and Development，1998，25(6):63-66.

［9］ 李克文，戴志堅(jiān)，江義容.穩(wěn)定流方法與非穩(wěn)定流方法的對(duì)比研究［J］.石油勘探與開(kāi)發(fā)，1995，22(1):47-51.

LI Ke-wen，DAI Zhi-jian，JIANG yi-rong.A comparison of water-oil relative permeability by steady-state flow and unsteady-state flow techniques［J］.Petroleum Exploration and Development，1995，22(1):47-51.

［10］ 李道品.低滲透油田高效開(kāi)發(fā)決策論［M］.北京:石油工業(yè)出版社，2003:8-9.

［11］ 黃延章.低滲透油層滲流機(jī)理［M］.北京:石油工業(yè)出版社，1998:12.

［12］ 楊正明，張英芝，郝明強(qiáng)，等.低滲透油田儲(chǔ)層綜合評(píng)價(jià)方法［J］. 石油學(xué)報(bào)，2006，27(2):64-67.

YANG Zheng-ming， ZHANG Ying-zhi， HAO Mingqiang，et al.Comprehensive evaluation of reservoir in low-permeability oilfields ［J］.Acta Petrolei Sinica，2006，27(2):64-67.

［13］ 楊正明，李治碩，王學(xué)武，等.特低滲透油田相對(duì)滲透率曲線測(cè)試新方法［J］.石油學(xué)報(bào)，2010，31(4):629-632.

YANG Zheng-ming，LI Zhi-shuo，WANG Xue-wu，et al.A new method for testing relative permeability of ultra-low permeability reservoirs［J］.Acta Petrolei Sinica，2010，31(4):629-632.

［14］ 林玉保，劉春林，衛(wèi)秀芬，等.特低滲透儲(chǔ)層油水滲流特征研究［J］.大慶石油地質(zhì)與開(kāi)發(fā)，2005，24(6):42-44.

LIN Yu-bao，LIU Chun-lin，WEI Xiu-fen，et al.Oil and water flow characteristics in ultra-low permeable reservoirs［J］.Petroleum Geology＆ Oilfield Development in Daqing，2005，24(6):42-44.

［15］ 高慧梅，姜漢橋，陳民鋒.儲(chǔ)層孔隙結(jié)構(gòu)對(duì)油水兩相相對(duì)滲透率影響微觀模擬研究［J］.西安石油大學(xué)學(xué)報(bào):自然科學(xué)版，2007，22(2):56-60.

GAO Hui-mei， JIANG Han-qiao， CHEN Min-feng.Simulation study on the effect of the microscopic parameters of reservoir pore structure on oil-water relative permeability［J］.Journal of Xi'an Shiyou University(Natural Science Edition)，2007，22(2):56-60.

［16］ 李克文，沈平平，秦同洛.考慮毛管壓力時(shí)計(jì)算油水相對(duì)滲透率的新方法［J］.力學(xué)與實(shí)踐，1994，16(2):46-48.

LI Ke-wen，SHEN Ping-ping，QIN Tong-luo.A new approach for calculating oil-water relative permeabilities considering the capillary pressure［J］.Mechanics in Engineering，1994，16(2):46-48.

［17］ 鄧英爾，王允誠(chéng)，劉慈群，等.低滲非達(dá)西滲流相對(duì)滲透率計(jì)算方法及特征［J］.西南石油學(xué)院學(xué)報(bào)，2000，22(3):34-36.

DENG Ying-er，WANG Yun-cheng，LIU Ci-qun，et al.Calculating relative permeability of non-Darcy flow in low permeability media［J］.Journal of Southwest Petroleum Institute，2000，22(3):34-36.

［18］ CIVAN F，DONALDSON E C.Relative permeability from unsteady-state displacements with capillary pressure included ［J］.Journal of SPE Formation Evaluation，1989，4(2):189-193.

［19］ 呂成遠(yuǎn).油藏條件下油水相對(duì)滲透率實(shí)驗(yàn)研究［J］.石油勘探與開(kāi)發(fā)，2003，30(4):102-105.

［20］ 李武廣，楊勝來(lái)，李芳，等.冷凍與常規(guī)巖心油水相對(duì)滲透率實(shí)驗(yàn)對(duì)比［J］.中國(guó)海上油氣，2012，24(3):41-43.

LI Wu-guang，YANG Sheng-lai，LI Fang，et al.An experimental comparison of oil-water relative permeability between frozen and conventional cores［J］.China Offshore Oil and Gas，2012，24(3):41-43.