何家歡

(1.中國(guó)石油西南油氣田分公司，四川 成都 610051;2.四川省頁(yè)巖氣評(píng)價(jià)與開(kāi)采重點(diǎn)實(shí)驗(yàn)室，四川 成都 610051;3.西南石油大學(xué)，四川 成都 615000)

0 引言

頁(yè)巖氣作為具有良好潛力的非常規(guī)天然氣資源，在水平井及水力壓裂等生產(chǎn)技術(shù)的不斷推動(dòng)下，近年來(lái)取得了快速發(fā)展，中國(guó)頁(yè)巖氣勘探開(kāi)發(fā)總體上仍處于起步階段。國(guó)內(nèi)針對(duì)頁(yè)巖儲(chǔ)層物性的相關(guān)實(shí)驗(yàn)主要依靠引進(jìn)國(guó)外設(shè)備，但對(duì)于測(cè)試頁(yè)巖基質(zhì)滲透率的原理、理論模型研究卻不多。

國(guó)外針對(duì)頁(yè)巖滲透率的技術(shù)研究起步較早，早在1968年，Brace[1]針對(duì)常規(guī)穩(wěn)態(tài)法測(cè)低滲樣品流量太小不便計(jì)量的困難，建立壓力脈沖衰減法。1988年，Dicker[2]解決了傳統(tǒng)的壓力脈沖衰減法實(shí)驗(yàn)時(shí)間過(guò)長(zhǎng)的問(wèn)題，建立改良的壓力脈沖衰減法。1992年，Xiuxu Ning詳細(xì)闡述了壓力脈沖衰減法的基本原理，建立針對(duì)均質(zhì)柱塞樣、含有裂縫柱塞樣的數(shù)學(xué)模型。1993年，Luffel列舉3種當(dāng)時(shí)在美國(guó)常用測(cè)試頁(yè)巖基質(zhì)滲透率的方法，其中包括壓力脈沖衰減法、粉碎法(GRI方法)和脫氣測(cè)試法。2011年，李智鋒[3]對(duì)原地測(cè)試法、改良的壓力脈沖衰減法和解吸流動(dòng)法3種針對(duì)頁(yè)巖儲(chǔ)層滲透率測(cè)試方法作了對(duì)比研究。2012年于榮澤[4－5]也提到Luffel列舉的3種方法，只是將 GRI方法翻譯成“巖屑脈沖衰減法”。2013年，何家歡[6]根據(jù)粉碎法測(cè)頁(yè)巖基質(zhì)滲透率的實(shí)驗(yàn)流程，考慮不同頁(yè)巖顆粒形態(tài)、直徑大小，建立相應(yīng)的數(shù)學(xué)模型，通過(guò)實(shí)測(cè)數(shù)據(jù)與理論計(jì)算曲線相擬合的方式，測(cè)試了頁(yè)巖的基質(zhì)滲透率。

1 問(wèn)題的提出

一些學(xué)者的研究[7－18]表明，頁(yè)巖主要由有機(jī)物與無(wú)機(jī)物混合組成，其中有機(jī)物(干酪根)的孔隙主要由微孔隙(孔徑小于2 nm)和中孔隙(孔徑為2～50 nm)組成，平均孔徑小于5 nm。由于有機(jī)物孔隙占總孔隙比例大，且孔隙非常小，烴類中甲烷分子直徑(約0.38 nm)雖然遠(yuǎn)小于頁(yè)巖小孔隙的平均直徑，但由于小孔隙擁有較大的比表面積，因此大量甲烷及烴類分子吸附在表面上(常規(guī)天然氣由于吸附量少，可不予考慮)，而且孔隙空間狹小，分子平均自由程變短。在不發(fā)育裂縫的情況下，大多數(shù)頁(yè)巖的滲透能力非常差，氣體的運(yùn)移方式主要依靠擴(kuò)散。李曉強(qiáng)[19]在2011年通過(guò)數(shù)學(xué)計(jì)算認(rèn)為，忽略基質(zhì)中的擴(kuò)散流動(dòng)將嚴(yán)重低估氣井的產(chǎn)量，且滲透率越低，估算的誤差越大。

目前國(guó)內(nèi)外主要測(cè)試頁(yè)巖滲透率的方法均未考慮擴(kuò)散的影響。在美國(guó)，頁(yè)巖基質(zhì)滲透率為1.0×10－8～1.0 ×10－4μm2時(shí)使用壓力脈沖衰減法，基質(zhì)滲透率小于1.0×10－8μm2時(shí)采用粉碎法測(cè)試更為合適。WeatherFord、TerraTek等實(shí)驗(yàn)室能夠利用粉碎法測(cè)頁(yè)巖基質(zhì)滲透率，最初提出該方法的研究機(jī)構(gòu)命名為GRI(Gas Research Institute，天然氣研究院)方法，后來(lái)又結(jié)合頁(yè)巖儲(chǔ)層的其他物性分析項(xiàng)目，統(tǒng)稱為SRP(Shale Rock Property，頁(yè)巖物性分析技術(shù))方法。不同的公司根據(jù)GRI方法研制了不同的實(shí)驗(yàn)系統(tǒng)，分別命名不同的名稱，如:TerraTek實(shí)驗(yàn)室將其稱為TRA(Tight Rock Analysis，致密巖分析技術(shù))方法，美國(guó)巖心公司則命名為SMP(Shale Matrix Permeability)方法。雖然名稱各異，但這些實(shí)驗(yàn)系統(tǒng)的流程和原理類似。粉碎法測(cè)頁(yè)巖基質(zhì)滲透率技術(shù)目前在國(guó)內(nèi)外的頁(yè)巖滲透率測(cè)試中得到了廣泛使用，但該方法模型和壓力脈沖衰減法的模型一樣，并未考慮擴(kuò)散的影響，而是在滲流方程中沿用了達(dá)西定律的一般形式。

隨著頁(yè)巖滲流機(jī)理研究的深入，出現(xiàn)了一個(gè)新的問(wèn)題，不考慮擴(kuò)散的影響會(huì)影響測(cè)試頁(yè)巖的滲透率嗎?要回答上述問(wèn)題，必須要明白“滲透率”和“擴(kuò)散”2個(gè)專業(yè)名詞的意義。

2 滲透率概念的適用性

所謂滲透率[20]，是指在一定壓差下，巖石允許流體通過(guò)的能力。柱塞樣氣測(cè)滲透率測(cè)試具體做法是:用加壓氣體(用氮?dú)馄炕驂猴L(fēng)機(jī))方法在巖樣兩端建立壓力差，測(cè)量進(jìn)、出口壓力及出口流量，按式(1)可計(jì)算巖石氣測(cè)滲透率。

式中:Ka為氣測(cè)滲透率，μm2;Q0為在壓差Δp下通過(guò)巖心的流量，cm3/s;A為巖心截面積，cm2;L為巖心長(zhǎng)度，cm;μ為通過(guò)巖心的氣體黏度，mPa·s;p1、p2分別為巖心進(jìn)口端、出口端的壓力，MPa;p0為穩(wěn)定壓力，MPa。

參照定義，滲透率主要表征巖石在一定壓差下允許流體通過(guò)的能力，式(1)計(jì)算出來(lái)的值是對(duì)在某個(gè)壓差狀態(tài)下其滲透率最直觀的表示，至于流體是否滿足達(dá)西定律與滲透率測(cè)試沒(méi)有任何關(guān)系。例如，實(shí)測(cè)低滲樣品[21]的滲流指示曲線(圖1)包括高速直線段和低速曲線段，兩部分流體通過(guò)巖石的能力并不相同，因此同一巖心樣品在不同的壓差環(huán)境下具有2個(gè)(甚至更多)流體通過(guò)能力，即滲透率。顯然，遇到圖1中的情況，滲透率測(cè)試的實(shí)驗(yàn)結(jié)果出現(xiàn)2個(gè)值，并不意味著滲透率測(cè)試不準(zhǔn)確，相反，這更為直接地反映了巖石在不同條件下允許流體通過(guò)的能力。

圖1 實(shí)測(cè)滲流指示曲線(李傳亮，2013年)

對(duì)于低滲樣品(滲透率為1.0×10－8～1.0×10－4μm2)面臨的氣體難以準(zhǔn)確計(jì)量的問(wèn)題，美國(guó)通常采用壓力脈沖衰減法。該方法的數(shù)學(xué)模型中，滲流微分方程采用以下方式表述:

從式(2)的形式看，該方程滿足達(dá)西定律，由此引發(fā)爭(zhēng)議。不少學(xué)者認(rèn)為，式(2)沿用到低滲樣品時(shí)，并未在滲流方程中增加啟動(dòng)壓力梯度或者其他影響項(xiàng)。但研究后認(rèn)為，這并不會(huì)影響其滲透率值計(jì)算的準(zhǔn)確性。原因如下:壓力脈沖衰減法的精髓在于，巖心在穩(wěn)定壓力p0狀態(tài)時(shí)加一個(gè)很小的壓力脈沖Δp，記錄從壓力脈沖施加到壓力重新恢復(fù)平衡的壓力變化與時(shí)間的關(guān)系，經(jīng)過(guò)與理論模型計(jì)算曲線進(jìn)行對(duì)比擬合，得到相應(yīng)的滲透率。從整個(gè)實(shí)驗(yàn)過(guò)程來(lái)看，實(shí)驗(yàn)中壓力變化并不大，基于上述情況，該方程假設(shè)滲透率在壓力微小變化的情況下是一個(gè)常數(shù)，因此，該方法計(jì)算出的滲透率可看成是壓力從p0+Δp到p0時(shí)的一個(gè)平均滲透率，由于Δp相對(duì)于p0而言很小，壓力脈沖衰減法計(jì)算出的滲透率近似看成是p0時(shí)的滲透率。特別要指出的是，式(2)能夠適用，是因?yàn)樵摲匠滩⑽聪薅B透率在所有情況下均為常數(shù)，只要從p0+Δp到p0滿足常數(shù)的要求，式(2)均成立。

上述方法是美國(guó)針對(duì)低滲樣品測(cè)試的主要手段，從工程應(yīng)用角度來(lái)看，該方法較好地反映了某一壓力時(shí)低滲樣品允許流體通過(guò)的能力。國(guó)內(nèi)有學(xué)者提出疑問(wèn)，認(rèn)為應(yīng)該在式(2)中增加啟動(dòng)壓力梯度項(xiàng)，這其實(shí)犯了一個(gè)邏輯上的錯(cuò)誤，所謂的啟動(dòng)壓力梯度本來(lái)是穩(wěn)態(tài)法實(shí)驗(yàn)測(cè)試結(jié)果理論分析后得到的概念，壓力脈沖衰減法實(shí)驗(yàn)測(cè)試結(jié)果反映的是實(shí)際壓力情況下樣品的滲透率值，這是對(duì)客觀事實(shí)最為直接的表征，其測(cè)試結(jié)果相當(dāng)于實(shí)驗(yàn)中的“原始數(shù)據(jù)”，實(shí)驗(yàn)的原始數(shù)據(jù)如果還未采集就遭到理論上的“修正”，則已失去了本身存在的意義。

美國(guó)針對(duì)頁(yè)巖常用的GRI方法以及何家歡提出的粉碎法測(cè)頁(yè)巖基質(zhì)滲透率均是在壓力脈沖衰減法思路的基礎(chǔ)上，針對(duì)不同形狀的實(shí)驗(yàn)樣品建立的相應(yīng)數(shù)學(xué)模型并求解，其測(cè)試出的滲透率值也是直接客觀反映某一壓力時(shí)低滲樣品宏觀上允許流體通過(guò)的能力。至于流體是以哪種具體形式(包括究竟是不是擴(kuò)散)通過(guò)低滲樣品，并不在滲透率測(cè)試討論的范圍內(nèi)。

從目前頁(yè)巖氣勘探開(kāi)發(fā)的情況來(lái)看，其基質(zhì)滲透率測(cè)試結(jié)果的意義主要體現(xiàn)在:判斷頁(yè)巖氣藏是否具有開(kāi)發(fā)價(jià)值;壓裂施工中確定有利層段;為氣藏?cái)?shù)值模擬提供數(shù)據(jù)支撐。以滲透率作為特征參數(shù)，較使用擴(kuò)散系數(shù)或者其他參數(shù)更能達(dá)到評(píng)價(jià)頁(yè)巖儲(chǔ)層物性的目的，更容易實(shí)現(xiàn)上述目標(biāo)。國(guó)外學(xué)者認(rèn)為，當(dāng)滲透率小于1.0×10－12μm2時(shí)，頁(yè)巖的基質(zhì)滲透率太低，沒(méi)有任何經(jīng)濟(jì)價(jià)值;當(dāng)滲透率為1.0×10－12～1.0×10－9μm2時(shí)，基質(zhì)滲透率是頁(yè)巖氣井產(chǎn)量的一個(gè)重要控制因素;當(dāng)滲透率大于1.0×10－9μm2時(shí)，基質(zhì)滲透率與產(chǎn)量關(guān)系并不大，其產(chǎn)量主要受裂縫的性質(zhì)控制。這里的滲透率均是代表巖石允許流體通過(guò)的能力，GRI方法和粉碎法測(cè)頁(yè)巖基質(zhì)滲透率均能直接滿足頁(yè)巖氣勘探開(kāi)發(fā)最基本的需要，具有明確的指示性和較好的實(shí)用性。

3 擴(kuò)散概念的適用性

擴(kuò)散(Diffusion)是指分子通過(guò)隨機(jī)分子運(yùn)動(dòng)從高濃度區(qū)域向低濃度區(qū)域的網(wǎng)狀傳播，擴(kuò)散的結(jié)果是緩慢地將物質(zhì)混合起來(lái)，在溫度恒定的空間中，忽略外部分子的相互作用力，擴(kuò)散的過(guò)程是完全混合或達(dá)到一種平衡狀態(tài)。

從擴(kuò)散的定義和解釋來(lái)看，頁(yè)巖氣開(kāi)發(fā)主要依靠擴(kuò)散并不合理:①頁(yè)巖氣在開(kāi)發(fā)過(guò)程中，即便是在基質(zhì)中，只要宏觀上有氣體定向運(yùn)移，就意味著該系統(tǒng)不是一個(gè)封閉系統(tǒng)，認(rèn)為其是擴(kuò)散，并不滿足擴(kuò)散概念上的要求，套用Fick定律、Knudsen定律均不符合物理事實(shí);②擴(kuò)散微觀上隨機(jī)無(wú)序，不可能形成定向的氣流，決定氣體流動(dòng)的情況是頁(yè)巖基質(zhì)中壓力大于裂縫中的壓力，實(shí)質(zhì)上還是滲流的特點(diǎn)。

擴(kuò)散過(guò)程符合熱力學(xué)第二定律，參照熱力學(xué)第二定律的推論——熵增加原理，“對(duì)于一個(gè)孤立系統(tǒng)，其內(nèi)部自發(fā)進(jìn)行的與熱相關(guān)的過(guò)程必然向熵增的方向進(jìn)行，而孤立系統(tǒng)不受外界任何影響，且系統(tǒng)最終處于平衡態(tài)”。頁(yè)巖氣開(kāi)發(fā)過(guò)程中，氣體通過(guò)裂縫滲流到井眼，而裂縫與基質(zhì)連通，基質(zhì)與裂縫接觸的面與基質(zhì)內(nèi)部必然存在一個(gè)壓力差，并不符合絕熱孤立系統(tǒng)的條件，因此，頁(yè)巖基質(zhì)中的氣體顯然不滿足這樣的條件，即便很多學(xué)者認(rèn)為，氣體分子的平均自由程與有機(jī)孔直徑在同一數(shù)量級(jí)，將宏觀上的運(yùn)移認(rèn)為是擴(kuò)散引發(fā)也并不合適。

綜上所述，頁(yè)巖氣開(kāi)發(fā)過(guò)程中考慮基質(zhì)中氣體的運(yùn)移機(jī)理是否是擴(kuò)散有待商榷，而認(rèn)為測(cè)試頁(yè)巖基質(zhì)滲透率的實(shí)驗(yàn)該考慮擴(kuò)散的影響同樣是無(wú)從說(shuō)起。

4 頁(yè)巖樣品的特殊性

頁(yè)巖儲(chǔ)層脆性礦物含量高，有利于儲(chǔ)層壓裂改造，卻也給柱塞巖樣的制備帶來(lái)了巨大的難度。從地下條件到地表?xiàng)l件存在一個(gè)卸壓過(guò)程，由于脆性礦物含量高，巖石會(huì)因?yàn)閼?yīng)力的變化而誘導(dǎo)新的裂縫產(chǎn)生，帶來(lái)更多的不可控因素。以N3井龍馬溪組頁(yè)巖儲(chǔ)層為例，通過(guò)常規(guī)滲透率測(cè)試方法測(cè)得滲透率大于1.0×10－4μm2的樣品數(shù)占到54.9%，如果此測(cè)試結(jié)果能夠反映真實(shí)地層情況，那么該井不壓裂也應(yīng)具備一定的滲流能力，但實(shí)際的情況是頁(yè)巖氣井不壓裂根本無(wú)法生產(chǎn)。一些學(xué)者往往把具有裂縫、微裂縫的柱塞樣品進(jìn)行擴(kuò)散系數(shù)的測(cè)試，顯然不再適用，此時(shí)頁(yè)巖中氣體運(yùn)移的主要通道是裂縫、微裂縫，不再是那些平均孔徑小于5 nm的孔隙喉道。

5 結(jié)論

(2)滲透率是體現(xiàn)巖石允許流體通過(guò)能力大小的一個(gè)宏觀表征量，與流體究竟以哪種方式通過(guò)沒(méi)有直接關(guān)系，氣體是否以擴(kuò)散的方式通過(guò)頁(yè)巖，并不影響頁(yè)巖滲透率的測(cè)定。

(3)頁(yè)巖氣開(kāi)發(fā)過(guò)程中，氣體在基質(zhì)中的運(yùn)移并不符合物理化學(xué)上“擴(kuò)散”的基本概念，用Fick定律或者Knudsen定律描述氣體在頁(yè)巖中的運(yùn)移時(shí)須謹(jǐn)慎。

[1]Brace，Walsh，F(xiàn)rangos.Permeability of granite under high pressure[J].J.Geophys.Res，1968，73(6):25 －36.

[2]Dicker，Smits.A practical approach for determining permeability from laboratory pressure－pulse decay measurements[C].SPE17578，1988:285 －292.

[3]李智鋒，李治平，王楊，等.頁(yè)巖氣儲(chǔ)層滲透率測(cè)試方法對(duì)比分析[J].斷塊油氣田，2011，18(6):761－764.

[4]于榮澤，卞亞南，張曉偉，等.頁(yè)巖儲(chǔ)層非穩(wěn)態(tài)滲透率測(cè)試方法綜述[J].科學(xué)技術(shù)與工程，2012，12(9):7019－7025.

[5]于榮澤，卞亞南，張曉偉，等.頁(yè)巖儲(chǔ)層流動(dòng)機(jī)制綜述[J]. 科技導(dǎo)報(bào)，2012，30(24):75 －79.

[6]何家歡，郭靜姝，王蘭生，等.一種頁(yè)巖氣體滲透率測(cè)定儀:中國(guó)，ZL201320142392.6[P].2013－10－16.

[7]閻存章，等.中國(guó)頁(yè)巖氣地質(zhì)研究進(jìn)展[M].北京:石油工業(yè)出版社，2011:61－63.

[8]秦佳，張威，劉晶，等.美國(guó)頁(yè)巖氣開(kāi)發(fā)狀況分析[J].大慶石油地質(zhì)與開(kāi)發(fā)，2014，33(4):170－174.

[9]李建青，高玉巧，花彩霞，等.北美頁(yè)巖氣勘探經(jīng)驗(yàn)對(duì)建立中國(guó)南方海相頁(yè)巖氣選區(qū)評(píng)價(jià)體系的啟示[J].油氣地質(zhì)與采收率，2014，21(4):23 －27，32.

[10]葛忠偉，樊莉.頁(yè)巖氣研究中應(yīng)注意的問(wèn)題[J].油氣地質(zhì)與采收率，2013，20(6):19－22.

[11]陳志明，汪偉英，蔡雨桐，等.致密砂巖和頁(yè)巖滲透率實(shí)驗(yàn)研究[J].斷塊油氣田，2013，20(1):80－84.

[12]曲秀英.擴(kuò)散系數(shù)實(shí)驗(yàn)技術(shù)研究[J].大慶石油地質(zhì)與開(kāi)發(fā)，2012，31(4):46 －49.

[13]郭平，王德龍，汪周華，等.頁(yè)巖氣藏儲(chǔ)層特征及開(kāi)發(fā)機(jī)理綜述[J].地質(zhì)科技情報(bào)，2012，31(6):118－123.

[14]王瑞，張寧生，劉曉娟，等.頁(yè)巖氣擴(kuò)散系數(shù)和視滲透率的計(jì)算[J].西北大學(xué)學(xué)報(bào):自然科學(xué)版，2013，43(1):75－88.

[15]糜利棟，姜漢橋，李俊鍵，等.頁(yè)巖儲(chǔ)層氣體擴(kuò)散機(jī)理[J]. 大慶石油地質(zhì)與開(kāi)發(fā)，2014，33(1):154－159.

[16]張磊，石軍太，楊先倫，等.頁(yè)巖裂縫和孔隙中氣體運(yùn)移方式評(píng)價(jià)及確定[J].大慶石油地質(zhì)與開(kāi)發(fā)，2014，33(4):166－169.

[17]孫海，姚軍，孫致學(xué)，等.頁(yè)巖氣數(shù)值模擬技術(shù)進(jìn)展及展望[J].油氣地質(zhì)與采收率，2012，19(1):46－49.

[18]郭瑞忠，汪洋，劉海成，等.頁(yè)巖氣生產(chǎn)機(jī)理及影響因素分析[J].特種油氣藏，2014，21(2):84－87.

[19]李曉強(qiáng)，周志宇，馮光，等.頁(yè)巖基質(zhì)擴(kuò)散流動(dòng)對(duì)頁(yè)巖氣井產(chǎn)能的影響[J].油氣藏評(píng)價(jià)與開(kāi)發(fā)，2011，5(1):67－70.

[20]何更生.油層物理[M].北京:石油工業(yè)出版社，2004:30－48.

[21]李傳亮，朱蘇陽(yáng).再談啟動(dòng)壓力梯度[J].巖性油氣藏，2013，25(4):1 －5.