鄧 鵬，林 亮，鄂廣杏，朱蘇陽(yáng)，彭小龍，楊宇光

(1.油氣藏地質(zhì)與開(kāi)發(fā)工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，西南石油大學(xué)，成都 610500；2.中聯(lián)煤層氣有限責(zé)任公司，北京 100016； 3.中石化普光分公司采氣廠，四川普光 636150)

隨著非常規(guī)天然氣藏的開(kāi)發(fā)，氣體在致密多孔介質(zhì)中的滑脫效應(yīng)越來(lái)越受到研究人員的關(guān)注[1-5]。煤巖的滲透率一般采用氣測(cè)滲透率實(shí)驗(yàn)，然后利用Klinkenberg提出的校正方法，延長(zhǎng)低壓條件下實(shí)驗(yàn)數(shù)據(jù)的擬合線，得煤巖的絕對(duì)滲透率[6-10]。然而，根據(jù)流體在微納通道中的流動(dòng)狀態(tài)，當(dāng)流體壓力達(dá)到一定數(shù)值后，流體在致密多孔介質(zhì)中的連續(xù)性變強(qiáng)，滑脫效應(yīng)消失[10-12]，因此，理論上氣測(cè)滲透率并非沿著低壓段的擬合線變化。因此，本文通過(guò)分析滑脫效應(yīng)實(shí)驗(yàn)中滲透率變化，提出了一種滲透率校正的改進(jìn)方法。

1 克氏校正

氣測(cè)滲透率實(shí)驗(yàn)中，低滲巖石的滲透率測(cè)量結(jié)果出現(xiàn)了較液測(cè)結(jié)果偏高且隨壓力變化的現(xiàn)象。該現(xiàn)象于1941年由L.J.Klinkenberg引入致密多孔介質(zhì)的滲透率測(cè)量實(shí)驗(yàn)，并將其解釋為氣體分子在巖石孔隙壁上的滑脫行為，即所謂的滑脫效應(yīng)或Klinkenberg效應(yīng)。氣測(cè)滲透率實(shí)驗(yàn)時(shí)，一般采用氣體滲流Darcy公式[13]

(1)

式中：qg為氣體在壓力p0下的流量cm3/s；p0為流量計(jì)算的基準(zhǔn)壓力(通常為大氣壓，MPa；p1為巖心進(jìn)口端壓力，MPa；p2為巖心出口端壓力，MPa；g為氣體在平均壓力下的黏度，mPa·s；kg為巖心的氣測(cè)滲透率，D；A為巖心的橫截面積，cm2；ΔL為巖心的橫截面積，cm；

氣體的物性取測(cè)試過(guò)程中氣體平均壓力時(shí)的黏度物性，而氣體的平均壓力為巖心入口端和出口端壓力的平均值，計(jì)算公式為

(2)

氣測(cè)滲透率實(shí)驗(yàn)時(shí)，氣體流動(dòng)偏離達(dá)西流動(dòng)，Klinkenberg對(duì)滲透率進(jìn)行校正，并通常采用式(3)對(duì)其進(jìn)行校正[14]

(3)

其中，b為滑脫因子，MPa。由公式(4)計(jì)算[15]

(4)

式中：c為近似于1的比例常數(shù)；λ為分子的平均自由程，m；r為平均孔隙半徑，m。

式(3)中氣測(cè)滲透率隨平均壓力增大而減小，當(dāng)平均壓力的倒數(shù)接近于0時(shí)，測(cè)試壓力接近于無(wú)窮大時(shí)，氣測(cè)滲透率接近液測(cè)滲透率(圖1)。因此，k被稱作等效液體滲透率或克氏滲透率[16]。

圖1 克氏校正中滲透率與平均壓力倒數(shù)的關(guān)系Figure 1 Relationship between permeability and reciprocal of average pressure in Klinkenberg modification

圖1中，由于實(shí)驗(yàn)條件限制，實(shí)測(cè)數(shù)據(jù)一般集中在低壓階段通過(guò)擬合低壓階段的直線關(guān)系并延長(zhǎng)得到克氏滲透率。然而，當(dāng)提高測(cè)試系統(tǒng)中的平均壓力時(shí)，克氏滲透率校正曲線并不沿直線變化。對(duì)此，本文開(kāi)展了實(shí)驗(yàn)研究并分析了曲線變化的原因。

2 實(shí)驗(yàn)測(cè)試

煤巖的克氏滲透率校正實(shí)驗(yàn)可以采用自研的孔滲聯(lián)測(cè)儀完成(圖2)，設(shè)備的進(jìn)氣口可以直接接入高壓氮?dú)馄俊?biāo)準(zhǔn)規(guī)格的氮?dú)馄繅毫υ?3MPa左右，因此實(shí)驗(yàn)過(guò)程中僅需要計(jì)量進(jìn)口壓力以及流量，并不需要驅(qū)替泵增加進(jìn)口端壓力。同時(shí)，為了更好的模擬平均氣體壓力對(duì)克氏滲透率的影響，不增加高速非達(dá)西效應(yīng)對(duì)滲透率測(cè)試過(guò)程的影響，實(shí)驗(yàn)不采用增加壓差提高系統(tǒng)平均壓力的方法，而是在設(shè)備中增加了回壓閥和壓力傳感器，通過(guò)保持穩(wěn)定的流動(dòng)壓差，提高回壓以提高測(cè)試系統(tǒng)中的平均壓力。

圖2 滲透率克氏校正實(shí)驗(yàn)設(shè)備Figure 2 Permeability Klinkenberg modification laboratory installations

克氏滲透率校正實(shí)驗(yàn)具體步驟如下：①大致保持流動(dòng)壓差在0.5～0.8MPa內(nèi)不變，以減少高速非達(dá)西效應(yīng)對(duì)滲透率測(cè)試的影響，緩慢增加回壓，使凈圍壓依次為0.2、0.5、1、1.5、2、2.5、3.0、3.5、4.0、4.5 、5.0、5.5MPa。系統(tǒng)中的平均壓力采用巖心的進(jìn)口壓力和出口壓力的平均值；②每一壓力點(diǎn)持續(xù)30min后，等待流動(dòng)穩(wěn)定，計(jì)量進(jìn)口壓力與出口壓力以及氣體流量，通過(guò)公式1計(jì)算煤巖的克氏滲透率滲透率。③繪制煤樣的滲透率與系統(tǒng)平均壓力的倒數(shù)曲線，研究克氏滲透率在高平均壓力下的響應(yīng)。本研究測(cè)試了兩塊來(lái)自貴州六盤水煤礦的礦井獲取的煤樣。

3 結(jié)果與討論

3.1 實(shí)驗(yàn)結(jié)果

根據(jù)克氏滲透率校正曲線可知(圖3)，低壓條件下，煤巖的克氏滲透率與平均壓力的倒數(shù)基本呈現(xiàn)線性變換。然而，在高壓條件下，煤巖的克氏滲透率對(duì)壓力變化并不敏感，基本處于不變的狀態(tài)。

根據(jù)測(cè)試結(jié)果可知，煤樣1屬于致密煤巖，克氏滲透率形狀變化的轉(zhuǎn)折壓力為3.97MPa，回壓校正方法獲得的滲透率為0.122mD，而克氏校正得到的液測(cè)滲透率為0.113mD，誤差在7.96%。而煤樣2 屬于低滲煤巖，克氏滲透率形狀變化的轉(zhuǎn)折壓力為1.87MPa，回壓校正方法獲得的滲透率為1.69mD，而克氏校正得到的液測(cè)滲透率為1.61mD，測(cè)試誤差在4.97%。由此可知，采用克氏滲透率校正方法或過(guò)低的估計(jì)煤巖的絕對(duì)滲透率。

(a)樣品1 致密煤巖 (b)樣品2 低滲煤巖圖3 煤巖心滑脫校正曲線Figure 3 Coal-rock core slippage modification curves

3.2 文獻(xiàn)數(shù)據(jù)驗(yàn)證

中-低滲巖心校正后的滲透率與原方法滲透率有一定的差距。W.Q. An[17]通過(guò)一定回壓下中-低滲巖心的滑脫流動(dòng)實(shí)驗(yàn)發(fā)現(xiàn)(圖4)，通過(guò)回壓校正方法得到的滲透率分別為0.32mD和3.87mD，而通過(guò)常規(guī)校正得到的滲透率為0.29mD和3.82mD，相對(duì)差距分別為6.89%和1.31%。實(shí)驗(yàn)得到的轉(zhuǎn)折壓力分別為3.84MPa與1.87MPa。

而校正后的滲透率差距在致密巖心上更加明顯，Shuliang Li等人[18]對(duì)5塊不同滲透率的巖心進(jìn)行變回壓實(shí)驗(yàn)發(fā)現(xiàn)，隨著回壓的逐漸增加，各巖心的滑脫效應(yīng)逐漸減弱，其中滲透率最低的S5巖心的回壓達(dá)到7.16MPa后，流量隨壓力梯度呈現(xiàn)線性關(guān)系，滑脫效應(yīng)消失(圖5)。

而滑脫校正過(guò)程中(圖5b),當(dāng)回壓大于7.16Mpa后，滑脫滲透率不隨流體壓力變化，計(jì)算得到滲透率為0.053mD，而通過(guò)常規(guī)滑脫校正后得到的滲透率為0.046mD,校正后的滲透率比常規(guī)方法得到的滲透率大15.21%。在利用LBM模擬得到的滲透率動(dòng)態(tài)存在相同的問(wèn)題[19]，文獻(xiàn)中通過(guò)格子氣方法，模擬得到的滲透率動(dòng)態(tài)中，常規(guī)滑脫校正后得到的滲透率為62nD，而經(jīng)過(guò)回壓校正之后滲透率為84nD，較常規(guī)方法大35.48%。因此，采用克氏校正方法會(huì)低估致密多孔介質(zhì)的滲透率。

(a)樣品1 低滲巖心 (b)樣品3 中滲巖心圖4 中-低滲巖心滑脫校正曲線Figure 4 Medium-low permeability core slippage modification curves

(a)流動(dòng)曲線 (b)滑脫校正圖5 S5巖心不同回壓流動(dòng)曲線Figure 5 Core S5 different back-pressure flow curves

3.3 流動(dòng)階段討論

高壓條件下，克氏滲透率曲線發(fā)生變化，這是流體連續(xù)性變強(qiáng)導(dǎo)致的。為研究流體流動(dòng)時(shí)的狀態(tài)，描述流體的連續(xù)性與離散性特征，通常采用流體力學(xué)中的Knudsen數(shù)進(jìn)行描述，Knudsen數(shù)的計(jì)算方法為[20]

(5)

式中：λ為分子平均自由程，m；D為圓筒的特征長(zhǎng)度，m。

氣體的平均自由程可以表示為

(6)

式中：kB為玻爾茲曼常數(shù)，取值為1.380 5×10-23J/K；T為氣體的熱力學(xué)溫度，K；p為氣體壓強(qiáng)，MPa；d為分子的直徑，m。

根據(jù)Knudsen數(shù)對(duì)流體流動(dòng)階段的劃分[17]如圖2所示，當(dāng)氣體分子處于稠密狀態(tài)時(shí)(Kn<10-3)，通常采用連續(xù)介質(zhì)方法研究流體的流動(dòng)。當(dāng)氣體壓力下降，分子稠密程度下降時(shí)，流體流動(dòng)不再符合連續(xù)介質(zhì)假設(shè)，但仍可以采用Navier-Stokes方程與Slip邊界層進(jìn)行描述。而當(dāng)氣體壓力進(jìn)一步下降，氣體分子處于過(guò)渡流狀態(tài)，氣體流動(dòng)沒(méi)有統(tǒng)一的描述方法，一般采用Burmett方程或統(tǒng)計(jì)方法進(jìn)行描述。當(dāng)Knudsen數(shù)大于10，只能采用Boltzmann分子運(yùn)動(dòng)方程或分子模擬的方法對(duì)流體進(jìn)行描述。

而滲透率滑脫效應(yīng)實(shí)驗(yàn)的校正方法，通常采用延長(zhǎng)低壓條件下滑移流動(dòng)擬合線的方法。然而，當(dāng)流體的平均壓力達(dá)到一定值時(shí)，氣體將從滑移狀態(tài)變?yōu)檫B續(xù)狀態(tài)，此時(shí)，滑脫效應(yīng)消失，滲透率將不隨壓力變化。因此，流體壓力達(dá)到一定數(shù)值后，氣測(cè)滲透率應(yīng)為一條像理論曲線一樣的水平直線(圖6)，而不是向長(zhǎng)滑移流動(dòng)直線的延長(zhǎng)線方向變化。然而，延長(zhǎng)滑移流動(dòng)階段的直線是目前克氏校正獲得液相滲透率的方式。綜上，在滑脫效應(yīng)實(shí)驗(yàn)對(duì)滲透率校正的計(jì)算中存在著與滲透率變化規(guī)律不符以及現(xiàn)有實(shí)驗(yàn)與方法矛盾的問(wèn)題。因此，在滑脫實(shí)驗(yàn)中對(duì)滲透率的校正并不妥當(dāng)。

4 結(jié)論

本文分析了分析煤巖在實(shí)驗(yàn)過(guò)程中不同回壓條件下的流動(dòng)階段，研究發(fā)現(xiàn)在滑脫流動(dòng)擬合線延長(zhǎng)過(guò)程中，流動(dòng)從滑移流動(dòng)變化為連續(xù)流動(dòng)的過(guò)程中，氣測(cè)滲透率并不符合滑移流動(dòng)延長(zhǎng)線的變化趨勢(shì)。研究通過(guò)實(shí)驗(yàn)以及文獻(xiàn)中的數(shù)據(jù)驗(yàn)證了提出了滲透率校正方法。提出了一種滲透率校正的改進(jìn)方法，通過(guò)研究得到以下結(jié)論：

(1)滑脫效應(yīng)校正方法采用延長(zhǎng)擬合直線的方法，然而在延長(zhǎng)過(guò)程中，流動(dòng)從滑移流動(dòng)變化為連續(xù)流動(dòng)，其滲透率變化并不符合滑移流動(dòng)延長(zhǎng)線的變化趨勢(shì)；

(2)煤巖以及低滲巖樣的液測(cè)滲透率校正若采用克氏校正方法，會(huì)低估巖石的液測(cè)滲透率，通過(guò)回壓校正方法，得到的巖心滲透率較常規(guī)滑脫校正方法得到的滲透率大，兩者的差距在中等滲透率巖心上并不明顯，滲透率越小，兩者差距越大。

圖6 根據(jù)克努森數(shù)劃分的流動(dòng)階段與滲流率變化情況Figure 6 Flow stages and permeability variation based on Knudsen number