涂彬 李杰

中國石油大學(xué)（北京）油氣資源與探測國家重點實驗室

0 引言

探井與開發(fā)井儲層含油飽和度分析是評價油藏潛力和開發(fā)效果、制定開發(fā)方案與調(diào)整方案等的重要依據(jù)和基礎(chǔ)［1］。求取地層含油飽和度的方法很多：李寶同應(yīng)用測井方法求取地層油水飽和度［2］；Botset［3］（1939）和 Rathmell［4］（1973）等利用常規(guī)巖心方法分析水驅(qū)油藏的剩余油飽和度；賈振堯（1995）等分析比較了常規(guī)密閉巖心與保壓密閉取心分析方法［5］，和李萬順［6］、馬力寧［7］等人得到了比較一致的認(rèn)識，都認(rèn)為常規(guī)密閉取心法相對準(zhǔn)確、直觀與經(jīng)濟(jì)，在國內(nèi)應(yīng)用很廣泛，為油田勘探開發(fā)研究試驗提供了可靠的油層物性數(shù)據(jù)，但他們也發(fā)現(xiàn)利用常規(guī)密閉取心方法測量的油水飽和度之和一般都小于100%，有的甚至低至70%以下，說明常規(guī)密閉取心法存在比較明顯的誤差，需要進(jìn)行校正，才能得到真實的地層流體飽和度。

針對常規(guī)密閉取心存在的精度不足，目前有4種類型的飽和度校正方法，即校正系數(shù)法、室內(nèi)實驗法、數(shù)理統(tǒng)計回歸法和數(shù)學(xué)模型校正法。Rathmell在文獻(xiàn)［3］中首先提出了校正系數(shù)法，Kazemi［8］（1977）和 Emmanuel［9］（1997）等人在此基礎(chǔ)上提出了不同的校正系數(shù)值。這類方法比較簡單，考慮的因素較少，校正得到的飽和度與實際相比誤差較大。楊勝來［10］（2004）、張亮［11］（2009）和劉麗［12］（2009）等人采用室內(nèi)實驗法對巖心飽和度進(jìn)行校正，葉青（2016）進(jìn)行了長巖心水驅(qū)油標(biāo)定飽和度損失［13］。這些實驗的過程比較類似，都是通過高壓下配制帶有溶解氣的實驗用油，然后降壓脫氣，測量脫氣前后水飽和度并回歸得到二者的關(guān)系式。該方法應(yīng)用時不用描繪降壓脫氣過程中復(fù)雜的物理化學(xué)變化過程，但復(fù)配地層油氣水組成及含量存在誤差。馬名臣等（1993）統(tǒng)計江蘇油田密閉取心數(shù)據(jù)，發(fā)現(xiàn)油水飽和度間存在線性關(guān)系，提出了數(shù)理統(tǒng)計校正方法［14］，王藝景等（2000）引入曲線平移和旋轉(zhuǎn)對該方法進(jìn)行了修正［15］；孫佩等（2012）在應(yīng)用該方法時先對不同孔隙結(jié)構(gòu)的巖心進(jìn)行分類，然后再進(jìn)行校正［16］。這些操作本質(zhì)上是對一定量的測試飽和度數(shù)據(jù)進(jìn)行統(tǒng)計回歸分析，因此歸類為數(shù)理統(tǒng)計回歸法。該方法比較簡單、容易理解，但未考慮巖心和流體性質(zhì)的差異性和取心過程中狀態(tài)變化，統(tǒng)計得到的損失率受不同巖石性質(zhì)的影響，樣本點的歸類標(biāo)準(zhǔn)不客觀，存在較大的誤差。辛治國（2012）基于流體分流率原理提出了密閉取心校正數(shù)學(xué)模型［17］，程會明（2012）建立了類似的模型［18］，王玉環(huán)（2014）采用指數(shù)函數(shù)處理辛治國和程會明模型中的水油相滲比［19］，總體上這3種模型的形式和求解方法基本不變。馬勇新（2016）引入了應(yīng)力敏感效應(yīng)推導(dǎo)了一種新的密閉取心飽和度校正模型，相對前面3種模型有了改進(jìn)，但其本質(zhì)仍然沒有改變，模型僅考慮了油和水的分流率卻沒有考慮溶解氣的影響［20］?？傮w來看，現(xiàn)有的數(shù)學(xué)模型方法考慮了從地下到地面后巖心孔隙體積和流體體積的改變，并基于分流理論來確定取心過程中流體的變化，但存在2個問題：一是僅僅考慮了從地下到地面這2個狀態(tài)下巖心孔隙體積變化和流體變化，沒有考慮整個壓降過程；二是沒有考慮取心時油氣溶解與分離的整個過程，只考慮了地面條件下巖心飽和度的分布結(jié)果，沒有體現(xiàn)出壓力下降過程中各種參數(shù)的變化過程。

基于目前研究方法及存在的問題，本文將建立一種動態(tài)校正模型，同時考慮巖心上提過程中因壓力不斷下降導(dǎo)致的孔隙形變和油氣分離逸散，以獲得較精確的密閉取心飽和度動態(tài)校正結(jié)果。

1 密閉取心飽和度動態(tài)校正模型

1.1 密閉取心飽和度影響因素

文獻(xiàn)［17-19］中，辛治國和王玉環(huán)等人綜合大量的研究成果指出，密閉取心測量的油水飽和度受到5個方面的影響：（1）巖石孔隙體積變化和流體體積變化；（2）巖心內(nèi)的原油降壓脫氣，出現(xiàn)溶解氣驅(qū)；（3）巖心至地面后轉(zhuǎn)運(yùn)過程中油水的揮發(fā)；（4）鉆井取心時鉆井液入侵；（5）實驗誤差。

文獻(xiàn)［17］中辛治國指出這5個因素中，實驗誤差無所不在，鉆井液侵入和巖心轉(zhuǎn)運(yùn)過程也不是很重要的影響因素，因此不做重點考慮。綜上所述，密閉取心巖心飽和度的主要影響因素是巖心上提過程中巖石孔隙體積變化、流體體積變化和溶解氣驅(qū)替，因此對這3種影響因素的校正也是眾多研究關(guān)注的重點問題。一般來說，這3個因素受溫度和壓力影響較大，而密閉取心在上提至地面時，由于時間較短，巖心筒內(nèi)的溫度變化幅度較小，壓力變化比較明顯，因此壓力的改變是導(dǎo)致巖心內(nèi)流體飽和度變化的主要原因，所以文中重點是分析壓力變化導(dǎo)致的飽和度變化。

1.2 取心過程油水飽和度動態(tài)變化

在上提過程中，壓力持續(xù)不斷下降，巖心內(nèi)的流體飽和度也是持續(xù)不斷變化的，這是一個動態(tài)的連續(xù)過程。可以基于離散原理將一個連續(xù)的物理過程用有限個階段過程來近似表示，即可以將巖心上提過程中連續(xù)的壓力變化用階段壓力來替代，每降低Δp，巖心內(nèi)的流體會出現(xiàn)如下的變化。

（1）壓力變化Δp后，如果巖心筒內(nèi)的壓力高于飽和壓力，則巖心內(nèi)只有孔隙體積和液體體積發(fā)生變化。在這二者的綜合作用下，巖心內(nèi)的液體被排擠出來，從而導(dǎo)致巖心內(nèi)流體飽和度的變化，此時排擠出來的油水體積可認(rèn)為是按照油水分流量分配。

（2）如果壓力變化Δp后，巖心筒內(nèi)的壓力低于飽和壓力，此時會出現(xiàn)溶解氣分離現(xiàn)象，因此除了液體膨脹、孔隙體積縮小之外，溶解氣也會膨脹且體積變化幅度更大，從而被排擠出來的流體體積更多，巖心內(nèi)流體體積的變化幅度也更大一些。

隨著巖心筒的繼續(xù)上行，筒內(nèi)壓力進(jìn)一步下降。巖心筒內(nèi)的飽和度變化會重復(fù)上述（1）或（2）的過程，當(dāng)壓力比飽和壓力低得多時，氣體析出量增加，成為連續(xù)相，因此氣體更多的是逸出巖心，且攜帶液體流出巖心，表現(xiàn)出比較明顯的溶解氣驅(qū)的特征。

1.3 密閉取心飽和度多相流體的校正方程

油氣兩相滲流或者是油氣水三相滲流的過程是不穩(wěn)定的，但在總過程中每一瞬間可近似看成穩(wěn)定狀態(tài)，這樣，總過程的不穩(wěn)定狀態(tài)就可以看成是無數(shù)個穩(wěn)定狀態(tài)的疊加，這種方法稱為穩(wěn)定狀態(tài)逐次替換法（引自葛家理［21］（1982））。當(dāng)時間間隔或者壓力間隔取得很小時，用此方法求得的結(jié)果將基本符合實際情況。注水開發(fā)油藏壓力一般大于其飽和壓力，巖心上提時壓力下降可分為泡點壓力之上和之下2個階段，因此在進(jìn)行飽和度校正時應(yīng)對這2個階段分別進(jìn)行考慮。

1.3.1 巖心壓力在泡點之上的校正方程 當(dāng)壓力下降，如果巖心壓力大于泡點壓力，溶解氣不分離，此時液體膨脹和巖石孔隙體積壓縮導(dǎo)致流體被排擠出來的流體總量為

按照物質(zhì)平衡原理，巖心降壓后剩余油量為降壓前油量與被排擠的油量之差，且認(rèn)為被排擠出來的油水的量按照分流量方程分配，因此可以得到

將式（1）中的總液量ΔNe帶入式（2）并化簡可以得到

其中，fw是含水率，按照分流量方程有

式（4）中kro和krw關(guān)于飽和度的函數(shù)，此時對應(yīng)的飽和度應(yīng)為降壓之前和之后的平均飽和度：Sw=（Sw+Sw-Δp）/2?？紤]到此時只有油水兩相，Sw+So=1，因此式（4）是一個關(guān)于飽和度Sw-Δp的隱函數(shù)方程，可以采用迭代法的方式求解得到Sw-Δp。

1.3.2 巖心壓力在泡點之下的校正方程 當(dāng)壓力下降后，如果巖心壓力小于泡點壓力，出現(xiàn)溶解氣分離，此時氣體膨脹、液體膨脹和巖石孔隙體積壓縮導(dǎo)致流體被排擠出來的流體總量ΔNe為

降壓比較小的時候，近似認(rèn)為在巖心內(nèi)的滲流是穩(wěn)定滲流，因此被排擠出巖心的油氣水三相可近似按照分流量fw、fo、fg來分配。

油組分的物質(zhì)平衡原理為：巖心降壓后剩余油量=降壓前油量-排出油量，由此可知油飽和度的校正方程為

同理可知水飽和度校正方程為

氣組分的物質(zhì)平衡原理為：降壓后游離氣量=降壓前游離氣量+降壓前溶解氣量-降壓后溶解氣量-排出氣量，因此可以得到氣體飽和度校正方程為

在上面的3個校正方程中，油水氣三者的分流量fw、fo、fg的表達(dá)式分別為

1.4 巖心測試飽和度校正的計算流程

根據(jù)巖心上提過程中油水飽和度變化的過程分析，可以采用圖1的流程來擬合實際巖心測試的飽和度。

圖1 巖心測試飽和度校正的計算流程圖Fig. 1 Calculation flow diagram of measured core saturation correction

（1）設(shè)定初始的油藏飽和度，如果注水開發(fā)油藏地層壓力高于飽和壓力，則可以得到Sg=0；

（2）降低壓力，判斷此時巖心筒內(nèi)的壓力是否低于泡點壓力，如果高于泡點壓力，則只需要根據(jù)油水兩相流來校正巖心飽和度的變化；如果低于泡點壓力，則巖心內(nèi)流體脫氣，出現(xiàn)溶解氣驅(qū)，此時需要用三相流來校正巖心飽和度的變化。

（3）計算降壓后的巖心飽和度，如果巖心筒內(nèi)的壓力仍然高于大氣壓或者至地面給定參考壓力，則繼續(xù)降低壓力，用此時計算的巖心飽和度替代初始飽和度，重復(fù)（2）的過程直至巖心筒內(nèi)的壓力降低至給定參考壓力。

（4）將計算得到的地面壓力條件下的巖心飽和度與實驗測試的巖心飽和度進(jìn)行比較，判斷這兩個數(shù)值之間的差值是否滿足誤差要求，如果滿足則說明初始給定的飽和度就是要尋求的油藏飽和度，反之則重新給定一個初始飽和度進(jìn)行計算。重復(fù)（1）～（4）的過程，直到滿足誤差為止。

2 巖心飽和度隨壓力變化及參數(shù)敏感性分析

2.1 巖心相關(guān)參數(shù)

密閉取心測量飽和度之后，根據(jù)建立的動態(tài)模型進(jìn)行校正，校正過程中，所需相關(guān)參數(shù)：地層水體積系數(shù)Bw為1.020，地層水黏度為0.304 mPa·s。地層油和天然氣的高壓物性數(shù)據(jù)見表1，包括原油黏度、天然氣黏度、原油體積系數(shù)、天然氣體積系數(shù)及溶解氣油比。

表1 油藏區(qū)塊PVT數(shù)據(jù)Table 1 PVT data of oil reservoir

圖2為原油和水的兩相相滲曲線，圖3為天然氣和液體的兩相相滲曲線。三相相滲可根據(jù)這2組相滲曲線應(yīng)用Stone I模型［22］計算得到（此處略）。

圖2 原油與水的相對滲透率曲線Fig. 2 Oil/water relative permeability curve

2.2 取心過程油氣水三相飽和度變化特征

利用油田取得的油藏性質(zhì)、巖石及流體物性資料，結(jié)合圖2、圖3的相對滲透率數(shù)據(jù)對巖心取心過程中的飽和度變化進(jìn)行計算。地層壓力為13.5MPa，飽和壓力為11.5 MPa，地層原油黏度為2.7 mPa·s。假設(shè)初始油水飽和度分別為50%和50%，采用本文模型（見1.4計算過程①～③）計算的結(jié)果可再現(xiàn)油氣水三相飽和度損失的過程，結(jié)果如圖4所示。

圖3 天然氣與液體的相對滲透率曲線Fig. 3 Gas/liquid relative permeability curve

圖4 半解析飽和度校正計算結(jié)果Fig. 4 Calculation result of semianalytical saturation correction

從圖4可以看出，隨著壓力的降低，巖心氣體飽和度迅速增加，而油和水的飽和度相應(yīng)降低。這說明氣體在膨脹過程中，會將一部分原油和水從巖心中排擠出去。當(dāng)巖心壓力降低至地面大氣壓時，氣體飽和度約為22%左右，原油的損失量約為10%，水的損失量約為11%。

2.3 密閉取心飽和度影響因素敏感性評價

根據(jù)實際油藏的參數(shù)，對影響密閉取心飽和度的因素進(jìn)行了敏感性分析，方案設(shè)計為：以基本參數(shù)為基準(zhǔn)值，通過曲線平移或拉伸等方法，使每個影響因素分別增加或減少30%、20%、10%，這樣每個影響因素共計算7個方案。計算結(jié)果見圖5和圖6，圖中橫坐標(biāo)為各影響參數(shù)的變化幅度，縱坐標(biāo)分別為剩余在巖心中油水飽和度。

從圖5和圖6可以看出每個因素的影響程度，并得到以下6個方面的認(rèn)識。

（1）氣體黏度對巖心油水飽和度的影響最大。泡點時氣體黏度從比基準(zhǔn)值小30%增加至比基準(zhǔn)值大30%時，巖心油飽和度減少3.3%，巖心水飽和度減少2.3%，巖心綜合油水飽和度改變率為9.3%。說明氣體黏度越高，溶解氣運(yùn)移攜帶油和水的能力越強(qiáng)。其中，巖心綜合油水飽和度改變率=（巖心油飽和度改變值+巖心水飽和度改變值）÷對應(yīng)參數(shù)變化幅度百分比×100%。

圖5 巖心相關(guān)參數(shù)變化與含油飽和度關(guān)系曲線Fig. 5 Relationship of variation of core parameters vs. oil saturation

圖6 巖心相關(guān)參數(shù)變化與含水飽和度關(guān)系曲線Fig. 6 Relationship of variation of core parameters vs. water saturation

（2）溶解氣油比對巖心油水飽和度的影響較大，且對巖心油水飽和度的影響并不是線性增加。當(dāng)泡點時溶解氣油比從比基準(zhǔn)值小30%增加至比基準(zhǔn)值大30%時，巖心油飽和度減少2.6%，巖心水飽和度減少2.0%，巖心綜合油水飽和度改變率為7.6%。

（3）原油體積系數(shù)對巖心油和水的飽和度影響程度也都比較大，原油體積系數(shù)越大，巖心中的油和水飽和度越高。當(dāng)泡點時原油體積系數(shù)從比基準(zhǔn)值小30%增加至比基準(zhǔn)值大30%時，巖心油飽和度增加2.2%，巖心水飽和度增加2.0%，巖心綜合油水飽和度改變率為6.9%。由于此處設(shè)計的情形可能導(dǎo)致原油體積系數(shù)小于1，因此計算結(jié)果僅具參考意義。為確保體積系數(shù)大于1，本文考慮了原油體積系數(shù)從比基準(zhǔn)值小3%增加至比基準(zhǔn)值大3%的情況，計算結(jié)果顯示，此時巖心油飽和度改變0.25%，巖心水飽和度改變0.2%，巖心綜合油水飽和度改變率為7.5%，和6.9%相差不大。

（4）原油黏度對巖心油飽和度的影響比對巖心水飽和度的影響大許多。當(dāng)泡點時原油黏度從比基準(zhǔn)值小30%增加至比基準(zhǔn)值大30%時，巖心油飽和度增加2.5%，而水飽和度減少0.8%，巖心綜合油水飽和度改變率為5.5%。說明油黏度越大，巖心油飽和度越高，被溶解氣攜帶出來的水越多。

（5）相滲對巖心油水飽和度的影響也各不相同，油相滲對巖心油飽和度的影響更大，水相滲對巖心水飽和度的影響更大。當(dāng)油相滲從比基準(zhǔn)值小30%增加至比基準(zhǔn)值大30%時，巖心油飽和度減少2.9%，水飽和度增加0.9%，巖心綜合油水飽和度改變率6.3%；當(dāng)水相滲從基準(zhǔn)值的-30%增加至30%時，巖心油飽和度僅增加0.1%，水飽和度卻減少3.2%，巖心綜合油水飽和度改變率為5.5%。

（6）壓縮系數(shù)對巖心油水飽和度的影響不大。其中巖石壓縮系數(shù)的影響最小，當(dāng)巖石壓縮系數(shù)從比基準(zhǔn)值小30%增加至比基準(zhǔn)值大30%時，巖心油水飽和度的變化值都不超過0.1%；原油壓縮系數(shù)僅對巖心油飽和度有影響，當(dāng)原油壓縮系數(shù)從比基準(zhǔn)值小30%增加至比基準(zhǔn)值大30%時，巖心油飽和度減少1.8%，水飽和度僅減少0.2%，巖心綜合油水飽和度改變率為3.3%。

綜合各參數(shù)的影響程度，可以認(rèn)為取心過程影響巖心飽和度變化的參數(shù)排序為：氣體黏度、溶解氣油比、原油體積系數(shù)、原油黏度、油水相滲和原油壓縮系數(shù)，影響最小的是巖石壓縮系數(shù)。

3 密閉取心飽和度動態(tài)校正模型驗證

3.1 油藏簡況

國內(nèi)西部某油田區(qū)塊，埋藏深度1 380 m左右，油層有效厚度約30 m，原始地層壓力13.7 MPa，飽和壓力11.4 MPa，地飽壓差2.3 MPa，原始含油飽和度平均為60%左右。該油藏以細(xì)砂巖、粉砂巖及中砂巖為主，孔隙度14.9%～12.5%，平均13.9%，滲透率10～27.2 mD，平均 15 mD。

3.2 密閉取心飽和度校正結(jié)果評價

W400217井地層水飽和度來自于生產(chǎn)單位，是基于巖心資料獲取相關(guān)飽和度計算參數(shù)，而后根據(jù)測井評價方法對巖心四性參數(shù)進(jìn)行標(biāo)定求取的流體飽和度值，和文中提供的動態(tài)模型計算結(jié)果進(jìn)行比較，結(jié)果如圖7所示，可以看出，密閉取心飽和度經(jīng)過模型校正后，與標(biāo)定飽和度的相對誤差范圍平均為0.8%，和靜態(tài)模型［17］計算結(jié)果（相對誤差平均為1.8%）相比，與標(biāo)定飽和度更接近，說明文中提供的動態(tài)校正方法具有一定可信度。從動態(tài)模型的建立和計算過程來看，本文模型在2個方面取得了進(jìn)展：其一是考慮了取心過程中壓力的動態(tài)變化過程，更符合飽和度變化的實際情況；其二是同時考慮孔隙體積變化和流體體積變化，更符合取心過程中巖心的相態(tài)變化規(guī)律。

圖7 W400217井含水飽和度測試及校正結(jié)果Fig. 7 Water saturation measured in Well 400217 and its correction result

4 結(jié)論

（1）密閉取心過程中，巖心內(nèi)流體飽和度變化是動態(tài)過程，孔隙和油水體積、原油脫氣等隨著壓力降低持續(xù)不斷變化，這是造成密閉取心測試飽和度出現(xiàn)誤差的根本原因，因此需要針對這些因素進(jìn)行飽和度校正。

（2）基于穩(wěn)定狀態(tài)逐次替換原理和多相流分流量理論，建立了巖心油水飽和度動態(tài)校正數(shù)學(xué)模型。模型既考慮了壓力的動態(tài)變化過程，又將取心過程中孔隙體積和流體體積這二者的狀態(tài)變化過程同時進(jìn)行計算，更加符合取心時巖心中的相態(tài)變化規(guī)律。

（3）參數(shù)敏感性分析結(jié)果表明，溶解氣對巖心油水飽和度的影響都比較大，氣體黏度從標(biāo)準(zhǔn)值的-30%增加至30%時，油水飽和度一共減少5.6%；溶解氣油比對應(yīng)的油水飽和度共減少4.6%。原油體積系數(shù)的影響程度次之，該參數(shù)從-3%增加3%時，油水飽和度共增加0.45%左右。原油黏度和油相滲透率僅對巖心油飽和度影響較大，而水相滲僅對巖心水飽和度影響較大。巖石壓縮系數(shù)對巖心油水飽和度的影響較小，影響程度不超過0.1%。綜合而言，影響程度大小排序依次為：氣體黏度、溶解氣油比、原油黏度、原油體積系數(shù)，油水相滲和原油壓縮系數(shù)，影響最小的是巖石壓縮系數(shù)。

符號說明：

Bo為原油體積系數(shù)；Bo-Δp、Bw-Δp、Bg-Δp分別為壓力下降時，壓力下的油、水和氣的體積系數(shù)；Bo-Δp/2、Bw-Δp/2、Bg-Δp/2分別為壓力下降時，壓力區(qū)間平均壓力下的油、水和天然氣的體積系數(shù)；Cf、Co、Cw分別為巖石、地層原油、地層水體積壓縮系數(shù)，MPa-1；fo、fw、fg分別為含油率、含水率和含氣率，%；kro、krw、krg分別為油、水和氣的相對滲透率；Rs為地層原油溶解氣油比，m3/m3；So、Sw、Sg分別為含油飽和度、含水飽和度和含氣飽和度，%；Swc、Sor分別為束縛水飽和度和殘余油飽和度，%；Vb為巖石外觀體積，m3；μo、μw、μg分別為油、水和氣的黏度，mPa·s；φ為孔隙度，%。