李傳亮 朱蘇陽(yáng)

西南石油大學(xué)石油與天然氣工程學(xué)院

0 引言

巖石壓縮系數(shù)是油氣藏工程的一個(gè)基本參數(shù)，在物質(zhì)平衡分析、試井解釋和數(shù)值模擬計(jì)算中有著重要的用途。傳統(tǒng)的巖石壓縮系數(shù)測(cè)量一直采用體積法，即通過(guò)測(cè)量巖石孔隙體積隨孔隙壓力的變化來(lái)確定巖石的壓縮系數(shù)[1-4]。由于存在表皮效應(yīng)，測(cè)量結(jié)果出現(xiàn)了邏輯反轉(zhuǎn)的現(xiàn)象，即孔隙度越高或巖石越疏松，巖石壓縮系數(shù)卻越小，意味著巖石越難壓縮；孔隙度越低或巖石越致密，巖石壓縮系數(shù)反而越大，意味著巖石越容易壓縮[5-6]。而且，實(shí)驗(yàn)測(cè)量值普遍大于地層水和油的壓縮系數(shù)，極端條件下甚至高于氣體。很顯然，這樣的實(shí)驗(yàn)結(jié)果對(duì)生產(chǎn)實(shí)踐沒(méi)有指導(dǎo)意義，反而會(huì)誤導(dǎo)生產(chǎn)實(shí)踐。

為了克服巖石壓縮系數(shù)測(cè)量方法上存在的缺陷，筆者于2003年在雙重有效應(yīng)力理論框架下提出了巖石壓縮系數(shù)的理論公式[6]，于2008年又提出了巖石壓縮系數(shù)的測(cè)量新方法—彈性模量法[7]，這才使得巖石壓縮系數(shù)的測(cè)量結(jié)果出現(xiàn)了正常的邏輯關(guān)系，即巖石孔隙度越高，巖石壓縮系數(shù)也越大，意味著巖石越容易壓縮；巖石孔隙度越低，巖石壓縮系數(shù)也越小，意味著巖石越難壓縮。而且，巖石的壓縮系數(shù)都普遍小于了流體的壓縮系數(shù)。這樣的實(shí)驗(yàn)測(cè)量結(jié)果不僅符合物理常識(shí)，而且對(duì)生產(chǎn)實(shí)踐真正具有了指導(dǎo)意義。

然而，很多人對(duì)新方法不理解，先是羅睿蘭提出了質(zhì)疑[8-9]，后是竇宏恩[10-12]、王厲強(qiáng)[13-14]和高有瑞[15]提出了質(zhì)疑，筆者都一一進(jìn)行了回復(fù)。最近，李子豐也撰文對(duì)其提出質(zhì)疑[16]，為了讓大家更好地理解巖石的壓縮性質(zhì)，筆者在此與李子豐作如下商榷。

1 關(guān)于巖石壓縮系數(shù)的概念問(wèn)題

與普通材料不同，巖石有骨架體積、孔隙體積和外觀體積共3個(gè)體積。巖石的物質(zhì)和體積構(gòu)成如圖1所示。

圖1 巖石物質(zhì)和體積構(gòu)成示意圖

巖石的3個(gè)體積滿足：

式中Vs表示骨架體積，m3；Vp表示孔隙體積，m3；Vb表示外觀體積，m3。

孔隙體積與外觀體積的比值，定義為巖石的孔隙度，即

式中φ表示孔隙度。

巖石通常同時(shí)受到外應(yīng)力（外壓，圍壓）和內(nèi)應(yīng)力（內(nèi)壓，孔隙壓力）的共同作用。在巖石的內(nèi)部，還存在第3個(gè)應(yīng)力，即骨架應(yīng)力。因此，巖石共有骨架應(yīng)力、內(nèi)應(yīng)力、外應(yīng)力3個(gè)應(yīng)力。巖石的受力和應(yīng)力構(gòu)成如圖2所示。

圖2 巖石的受力和應(yīng)力構(gòu)成示意圖

巖石的3個(gè)應(yīng)力滿足下面的應(yīng)力關(guān)系方程[17]：

式中σ表示外應(yīng)力，MPa；p表示內(nèi)應(yīng)力，MPa；σs表示骨架應(yīng)力，MPa。

當(dāng)巖石的任一應(yīng)力發(fā)生變化時(shí)，巖石的3個(gè)體積都將發(fā)生變化，因此，巖石有9個(gè)壓縮系數(shù)，其數(shù)學(xué)表達(dá)式分別為：

式中cbb表示外觀體積對(duì)外應(yīng)力的壓縮系數(shù)，MPa－1，簡(jiǎn)寫(xiě)成cb；cbs表示外觀體積對(duì)骨架應(yīng)力的壓縮系數(shù)，MPa－1；cbp表示外觀體積對(duì)內(nèi)應(yīng)力的壓縮系數(shù)，MPa－1；csb表示骨架體積對(duì)外應(yīng)力的壓縮系數(shù)，MPa－1；css表示骨架體積對(duì)骨架應(yīng)力的壓縮系數(shù)，MPa－1，簡(jiǎn)寫(xiě)成cs；csp表示骨架體積對(duì)內(nèi)應(yīng)力的壓縮系數(shù)，MPa－1；cpb表示孔隙體積對(duì)外應(yīng)力的壓縮系數(shù)，MPa－1；cps表示孔隙體積對(duì)骨架應(yīng)力的壓縮系數(shù)，MPa－1；cpp表示孔隙體積對(duì)內(nèi)應(yīng)力的壓縮系數(shù)，MPa－1，簡(jiǎn)寫(xiě)成 cp。

由于式（1）和式（3）的約束關(guān)系，3個(gè)體積和3個(gè)應(yīng)力中都分別只有2個(gè)是獨(dú)立的，因此，9個(gè)壓縮系數(shù)中只有4個(gè)是獨(dú)立的。

油氣藏工程只關(guān)心孔隙體積的變化，因?yàn)橹挥锌紫扼w積的變化對(duì)開(kāi)采油氣才有意義，其他體積的變化對(duì)開(kāi)采油氣沒(méi)有直接意義。由于巖石的外應(yīng)力（如上覆壓力）在油氣藏開(kāi)采過(guò)程中不發(fā)生變化，孔隙體積的變化都是因?yàn)榭紫秹毫Φ淖兓?，因此，在外?yīng)力為常數(shù)的情況下，巖石孔隙體積隨孔隙壓力的變化情況可以表征為[18-19]：

式（13）定義的壓縮系數(shù)實(shí)際上就是巖石孔隙體積對(duì)孔隙壓力的壓縮系數(shù)，但油氣藏工程習(xí)慣將其簡(jiǎn)稱為巖石壓縮系數(shù)，這樣做也是為了簡(jiǎn)便和與流體壓縮系數(shù)相對(duì)應(yīng)，因?yàn)橛蜌獠毓こ讨皇褂眠@一個(gè)巖石壓縮系數(shù)和3個(gè)流體壓縮系數(shù)（油、氣、水），在這種情況下使用簡(jiǎn)稱，并不會(huì)產(chǎn)生混亂。

實(shí)際上，巖石壓縮系數(shù)的測(cè)量工作屬于巖石物理學(xué)的研究范疇[1,19]，在該領(lǐng)域可以把名稱叫的詳細(xì)一些，而油氣藏工程是巖石壓縮系數(shù)的應(yīng)用領(lǐng)域，沒(méi)必要使用全稱，使用全稱反而顯得冗長(zhǎng)和繁瑣。就像發(fā)電領(lǐng)域需指明火電、水電與核電，而用戶則完全不關(guān)心電的來(lái)源。

可是，羅瑞蘭、竇宏恩和李子豐都覺(jué)得巖石壓縮系數(shù)的概念十分混亂，主張使用“巖石孔隙壓縮系數(shù)”而不是“巖石壓縮系數(shù)”。其實(shí)，巖石有3個(gè)孔隙壓縮系數(shù)（孔隙體積對(duì)外應(yīng)力、孔隙體積對(duì)內(nèi)應(yīng)力、孔隙體積對(duì)骨架應(yīng)力），即使采用巖石孔隙壓縮系數(shù)，仍然不清楚是指哪一個(gè)，并沒(méi)有提高嚴(yán)謹(jǐn)性。

2 關(guān)于一些工程認(rèn)識(shí)問(wèn)題

李子豐批評(píng)大家對(duì)一些工程問(wèn)題認(rèn)識(shí)不清，總共有4條[16]，在此一并回復(fù)。

“（1）巖石顆粒不可壓縮的假設(shè)不合理”[16]。這是土力學(xué)的一個(gè)基本假設(shè)，近似適用于疏松的土介質(zhì)，后來(lái)被巖石力學(xué)沿用，石油領(lǐng)域的孔隙度覆壓實(shí)驗(yàn)也一直采用這個(gè)假設(shè)。實(shí)際上，巖石為致密介質(zhì)，并不適合采用這個(gè)假設(shè)，筆者在巖石壓縮系數(shù)測(cè)量新方法中廢棄了這個(gè)假設(shè)[7]。土介質(zhì)屬于疏松介質(zhì)，受到應(yīng)力作用時(shí)會(huì)產(chǎn)生結(jié)構(gòu)變形，而且變形量特別大。所謂結(jié)構(gòu)變形，是指因骨架顆粒排列方式的改變而導(dǎo)致的巖石整體變形[17]。在結(jié)構(gòu)變形過(guò)程中，骨架顆粒自身的體積并不發(fā)生變化（因變形量小而忽略不計(jì)），變化的是骨架顆粒的排列方式（圖3）。圖4所示為結(jié)構(gòu)變形過(guò)程中的體積變化。很顯然，結(jié)構(gòu)變形過(guò)程中孔隙度發(fā)生了巨大變化。結(jié)構(gòu)變形實(shí)際上就是介質(zhì)的壓實(shí)變形（compaction），是一個(gè)不可逆過(guò)程。

圖3 多孔介質(zhì)的結(jié)構(gòu)變形示意圖

圖4 結(jié)構(gòu)變形的體積變化示意圖

“（2）巖石骨架與孔隙等比例脹縮的假設(shè)不合理。（3）孔隙度是常數(shù)的假設(shè)不合理”[15]。這兩條實(shí)際上說(shuō)的是同一個(gè)事情。與疏松的土介質(zhì)不同，巖石為致密介質(zhì)，經(jīng)過(guò)壓實(shí)和膠結(jié)作用之后，顆粒之間的排列方式便固定下來(lái)了，受到應(yīng)力作用之后不會(huì)產(chǎn)生結(jié)構(gòu)變形，即不會(huì)產(chǎn)生顆粒的重排，只會(huì)產(chǎn)生本體變形。所謂本體變形，是指因骨架顆粒自身的變形而導(dǎo)致的介質(zhì)整體變形[17]。在本體變形過(guò)程中，骨架顆粒的排列方式并不發(fā)生變化，變化的是骨架顆粒自身的體積（圖5）。圖6所示為本體變形過(guò)程中的體積變化。很顯然，本體變形過(guò)程中孔隙度為常數(shù)。本體變形實(shí)際上就是介質(zhì)的壓縮變形（compression），是一個(gè)可逆過(guò)程。

圖5 多孔介質(zhì)的本體變形示意圖

圖6 本體變形的體積變化示意圖

巖石孔隙度是粒度分布與排列方式的函數(shù)。對(duì)于特定的油氣藏巖石，粒度組成是確定的，若排列方式?jīng)]發(fā)生變化，孔隙度就不會(huì)發(fā)生變化，如圖5的本體變形；若排列方式發(fā)生了變化，孔隙度也跟著發(fā)生變化，如圖3所示的結(jié)構(gòu)變形。

巖石就像一個(gè)房間，若墻壁（骨架）不發(fā)生變化，房間的容積（孔隙）則不會(huì)發(fā)生變化。若墻壁變矮了，房間容積也跟著變小，房間容積是跟著墻壁等比例變化的。若房間垮塌了，即出現(xiàn)了結(jié)構(gòu)變形，房間容積則不隨墻壁等比例變化，而是變化量遠(yuǎn)大于墻壁。

實(shí)際上，巖石孔隙度隨埋藏深度（D）的變化一般如圖7所示，φo為地表松散沉積物的孔隙度，數(shù)值極高。在淺層（a點(diǎn)—b點(diǎn)）經(jīng)歷了壓實(shí)（結(jié)構(gòu)變形）階段，孔隙度大幅度減小。到了成巖深度（b點(diǎn)）以下，孔隙度則不再發(fā)生變化，而是一個(gè)常數(shù)，即進(jìn)入了壓縮（本體變形）階段。因此，巖石的孔隙度并不都是常數(shù)，只是在壓縮階段（b點(diǎn)—c點(diǎn)）才是常數(shù)。把地下巖石取到地面進(jìn)行實(shí)驗(yàn)，也只有壓縮變形，孔隙度并不會(huì)發(fā)生變化。

圖7 巖石孔隙度變化趨勢(shì)示意圖

“（4）巖石的平均壓應(yīng)力為常數(shù)的假設(shè)，在油氣藏內(nèi)一般不嚴(yán)格成立”[16]。不知道李子豐先生從哪里得到的這個(gè)認(rèn)識(shí)。實(shí)際上，油氣藏內(nèi)每一點(diǎn)的壓力都不相等，而是存在一個(gè)壓力分布，且隨時(shí)間發(fā)生變化。若用數(shù)值方法進(jìn)行研究，可以通過(guò)不同的網(wǎng)格壓力來(lái)模擬整個(gè)油氣藏。若用物質(zhì)平衡方法進(jìn)行研究，則采用油氣藏的平均壓力，相對(duì)較粗，但不影響巖石的壓縮性研究，因?yàn)橛闷骄鶋航涤?jì)算的膨脹量為用各點(diǎn)壓降計(jì)算的膨脹量之和。巖石壓縮系數(shù)通常不隨壓力變化，只是壓縮量或膨脹量隨壓力變化。

3 關(guān)于巖石壓縮系數(shù)與孔隙度的關(guān)系

巖石是由骨架顆粒組成的，顆粒之間為孔隙?？紫恫皇俏镔|(zhì)，只是多孔介質(zhì)的存在形式。人們無(wú)法給孔隙單獨(dú)施加壓力去壓縮孔隙，孔隙的壓縮是因?yàn)楣羌艿膲嚎s所致，骨架壓縮帶動(dòng)孔隙一起壓縮。骨架越硬，孔隙就越難壓縮；若骨架為剛體不可壓縮，則孔隙的壓縮系數(shù)為0。因此，孔隙的壓縮系數(shù)反映了骨架的壓縮系數(shù)。把骨架壓縮系數(shù)按孔隙度折算到孔隙上，即得孔隙壓縮系數(shù)。因此，孔隙的壓縮系數(shù)還反映了巖石的孔隙度。根據(jù)式（2），可以得出孔隙體積與骨架體積之間的關(guān)系：

根據(jù)式（3），在外應(yīng)力為常數(shù)的情況下，孔隙壓力與骨架應(yīng)力變化量之間的關(guān)系為：

把式（14）和式（15）代入式（5），得孔隙壓縮系數(shù)[6]：

式（16）就是折算后得到的孔隙壓縮系數(shù)，它是骨架壓縮系數(shù)和孔隙度的綜合反映，既反映了巖石的軟硬程度（通過(guò)骨架壓縮系數(shù)），又反映了巖石的疏松程度（通過(guò)孔隙度）。很顯然，巖石越硬，骨架壓縮系數(shù)越小，孔隙壓縮系數(shù)也就越小，巖石就越難壓縮。巖石越疏松，孔隙度越高，孔隙壓縮系數(shù)就越大，巖石就越容易壓縮。如果孔隙體積小于骨架體積（φ＜0.5），孔隙壓縮系數(shù)則小于骨架壓縮系數(shù)；如果孔隙體積大于骨架體積（φ＞0.5），孔隙壓縮系數(shù)則大于骨架壓縮系數(shù)。

李子豐認(rèn)為式（16）是錯(cuò)誤的，因?yàn)楫?dāng)φ＜0.5時(shí)，cp＜cs，并認(rèn)為孔隙不可能比骨架更難壓縮[15]。實(shí)際上，孔隙不是物質(zhì)，孔隙本身不能單獨(dú)被壓縮，不能把孔隙當(dāng)作物質(zhì)與骨架進(jìn)行對(duì)比，孔隙壓縮系數(shù)反映的是骨架壓縮系數(shù)和孔隙度的性質(zhì)。壓縮系數(shù)是一個(gè)強(qiáng)度指標(biāo)，是單位壓力下體積的變化率。根據(jù)式（15），若孔隙度為0.2（＜0.5），孔隙壓力增大1 MPa，骨架應(yīng)力才減小0.25 MPa，骨架應(yīng)力的變化量小于孔隙壓力，即骨架體積變化很小，孔隙體積當(dāng)然也跟著變化很小。當(dāng)孔隙度趨于0，孔隙壓力增大1 MPa，骨架應(yīng)力的減小量也趨于0，即骨架沒(méi)啥變化，孔隙體積當(dāng)然也沒(méi)啥變化，孔隙壓縮系數(shù)趨于0。骨架壓縮系數(shù)為常數(shù)，孔隙壓縮系數(shù)卻不是常數(shù)，而是一個(gè)變量?？紫抖雀撸▓D8-a），骨架少，支撐力度就弱，孔隙壓縮系數(shù)就大；孔隙度低（圖8-b），骨架多，支撐力度就強(qiáng)，孔隙壓縮系數(shù)就小。圖8-a的孔隙壓縮系數(shù)明顯大于圖8-b。式（16）恰好說(shuō)明了孔隙壓縮系數(shù)的這種變化規(guī)律，用其繪制的巖石壓縮系數(shù)曲線如圖9所示[6]。

圖8 高孔隙度與低孔隙度巖石示意圖

圖9 巖石壓縮系數(shù)曲線

4 關(guān)于鼓脹系數(shù)

李子豐先生提出了巖石孔隙鼓脹系數(shù)的概念，又重新定義了巖石孔隙壓縮系數(shù)的概念。

巖石孔隙鼓脹系數(shù)是指在圍壓為常數(shù)時(shí)，孔隙壓力增加一個(gè)單位壓力時(shí)，孔隙體積的變化率[15]。該系數(shù)是指因孔隙內(nèi)壓力增加，導(dǎo)致孔隙空間增加，不會(huì)產(chǎn)生歧義，可以用來(lái)計(jì)算孔隙體積變化量。巖石孔隙鼓脹系數(shù)的定義式，即原文式（13）[16]為：

式中cpp表示巖石孔隙鼓脹系數(shù)，MPa－1。

實(shí)際上，這個(gè)巖石孔隙鼓脹系數(shù)cpp就是大家常用的孔隙壓縮系數(shù)cp，已包括在9個(gè)壓縮系數(shù)之中了，只不過(guò)李子豐僅用它表征升壓過(guò)程，把降壓過(guò)程排除在外了。油氣藏生產(chǎn)是一個(gè)動(dòng)態(tài)過(guò)程，注入時(shí)升壓，采出時(shí)降壓，地層壓力一直處于變化之中，既然鼓脹系數(shù)只能表征升壓過(guò)程，降壓過(guò)程又該如何表征呢？在地層壓力波動(dòng)過(guò)程中它還有實(shí)際用途嗎？

李子豐重新定義的巖石孔隙壓縮系數(shù)是指在孔隙壓力為常數(shù)時(shí)，圍壓增加一個(gè)單位壓力時(shí)，孔隙體積的變化率[16]。該系數(shù)是指因巖石圍壓的增加，導(dǎo)致孔隙空間減小，不會(huì)產(chǎn)生歧義，可以用來(lái)計(jì)算孔隙體積變化量。巖石孔隙壓縮系數(shù)的定義式，即原文式（16）[16]為：

式中cps表示巖石孔隙壓縮系數(shù)，MPa－1。

實(shí)際上，這個(gè)巖石孔隙壓縮系數(shù)cps就是孔隙體積對(duì)外應(yīng)力的壓縮系數(shù)cpb，已包括在9個(gè)壓縮系數(shù)之中了，只不過(guò)李子豐又把它定義為了巖石孔隙壓縮系數(shù)，這與人們通常使用的孔隙壓縮系數(shù)完全不是一個(gè)概念，人們通常使用的巖石孔隙壓縮系數(shù)是孔隙體積對(duì)孔隙壓力的壓縮系數(shù)。由于油氣藏生產(chǎn)過(guò)程中外應(yīng)力（如上覆壓力）通常為常數(shù)，而且孔隙壓力也不可能為常數(shù)，因此，李子豐重新定義的這個(gè)巖石孔隙壓縮系數(shù)，根本找不到用處。

其實(shí)，任何物質(zhì)升壓都會(huì)壓縮，降壓都會(huì)膨脹。壓縮與膨脹是一個(gè)可逆過(guò)程，受壓力的變化方向控制。不管是壓縮，還是膨脹，都可以用壓縮系數(shù)來(lái)表征，物理學(xué)也沒(méi)有定義壓縮系數(shù)和膨脹系數(shù)來(lái)分別表征壓縮過(guò)程和膨脹過(guò)程。壓縮系數(shù)定義為單位壓力的體積變化率。壓力可以升高，也可以降低，由于壓力升高會(huì)導(dǎo)致體積壓縮和壓力降低會(huì)導(dǎo)致體積膨脹，為了讓壓縮系數(shù)取正值，定義式中加了一個(gè)負(fù)號(hào)，即

式中c表示物質(zhì)的壓縮系數(shù)，MPa－1；V表示物質(zhì)的體積，m3。

由式（19）可以計(jì)算體積的變化量：

壓力變化可正可負(fù)，體積變化也可正可負(fù)。若式（20）的dp為正，則為升壓過(guò)程；為負(fù)，則為降壓過(guò)程。若式（20）的dV為正，則為膨脹過(guò)程；為負(fù)，則為壓縮過(guò)程。式（20）既可以表征升壓壓縮過(guò)程，也可以表征降壓膨脹過(guò)程，一個(gè)公式完美地表征了兩個(gè)過(guò)程。數(shù)學(xué)上的正負(fù)，不就是為了表征兩個(gè)相反方向的過(guò)程而提出的嗎？還有必要自找麻煩分別表征嗎？

雖然式（19）和式（20）中的c稱為壓縮系數(shù)，但不一定必須用來(lái)表征壓縮過(guò)程，也同樣可以用來(lái)表征膨脹過(guò)程，因?yàn)榕蛎涍^(guò)程就是壓縮過(guò)程的逆過(guò)程，二者在數(shù)學(xué)上是統(tǒng)一的。物理學(xué)還定義了熱膨脹系數(shù)的概念用來(lái)表征物體遇熱膨脹、遇冷收縮的性質(zhì)，但并沒(méi)有定義膨脹系數(shù)和收縮系數(shù)來(lái)分別表征膨脹和收縮兩個(gè)過(guò)程。由此可見(jiàn)，李子豐定義巖石孔隙鼓脹系數(shù)來(lái)表征膨脹過(guò)程、定義巖石孔隙壓縮系數(shù)來(lái)表征壓縮過(guò)程，完全沒(méi)有必要，反而會(huì)使研究過(guò)程更加混亂。

5 結(jié)論

1）巖石有9個(gè)壓縮系數(shù)，油氣藏工程只使用孔隙體積對(duì)孔隙壓力的壓縮系數(shù)，并將其簡(jiǎn)稱為巖石壓縮系數(shù)，這樣做是為了簡(jiǎn)便和與流體壓縮系數(shù)相對(duì)應(yīng)，并不會(huì)產(chǎn)生混亂。

2）巖石的孔隙度在壓實(shí)過(guò)程中是不斷變化的，僅在壓縮過(guò)程中為常數(shù)。

3）巖石孔隙壓縮系數(shù)是骨架壓縮系數(shù)和孔隙度的綜合反映，可以大于骨架壓縮系數(shù)，也可以小于骨架壓縮系數(shù)。

4）巖石的膨脹和壓縮可以用同一個(gè)壓縮系數(shù)參數(shù)進(jìn)行表征，沒(méi)必要另外定義鼓脹系數(shù)來(lái)表征膨脹過(guò)程、定義壓縮系數(shù)來(lái)表征壓縮過(guò)程，應(yīng)用起來(lái)很不方便。