頁巖氣藏滲流數值模擬及井底壓力動態分析

尹 虎，王新海，姜 永，陳光喜，王婭妮(油氣資源與勘探技術教育部重點實驗室(長江大學)，湖北 荊州 434023)

汪金如 (華東石油局試采大隊，江蘇 泰州 225300)

頁巖氣藏滲流數值模擬及井底壓力動態分析

尹 虎，王新海，姜 永，陳光喜，王婭妮(油氣資源與勘探技術教育部重點實驗室(長江大學)，湖北 荊州 434023)

汪金如 (華東石油局試采大隊，江蘇 泰州 225300)

根據Langmuir等溫吸附方程和Fick擴散定律，結合滲流理論建立了頁巖氣井滲流數學模型，并對模型進行了數值求解，繪制了井底壓力變化曲線，分析了吸附氣解吸擴散、表皮系數、Langmuir壓力、Langmuir體積對井底壓力的影響。結果表明，考慮吸附氣解吸擴散，氣井井底壓力降低減慢；表皮系數增大，井筒附近污染嚴重，壓力降低加快；Langmuir壓力接近原始地層壓力，井底壓力降低最慢；Langmuir體積增大，地層中吸附氣解吸量增加，井底壓力降低緩慢。

頁巖氣；吸附解吸；擴散；數值模擬； 井底壓力

頁巖氣是以吸附和游離2種狀態儲存在于具有自身生氣能力的高碳泥頁巖或暗色泥頁巖中，也可存在粉砂巖和泥質粉砂巖中[1]。與常規氣藏相比，頁巖層既是生氣層，也是運移和保存天然氣的蓋層和儲層，表現為原地成藏的模式，具有飽含氣性、隱蔽聚集、含氣面積大、多種巖性封閉以及相對短距離運移等特點，其儲層孔隙度、滲透率較低、氣流阻力比常規天然氣大[2]。

頁巖氣以吸附狀態賦存于基質微孔隙內表面，這是與常規氣藏最主要的區別。頁巖氣與煤層氣相比，煤層氣中吸附氣占總量的70%～95%，裂縫中游離氣占總量的10%～20%[3]，而頁巖氣中吸附氣占總量的20%～85%，裂縫中游離氣占總量的50%左右。在開采過程中，煤層氣藏需要排水降壓開采，而頁巖氣藏在開采過程中就可使吸附氣降壓解吸，實現持續開采。在解吸過程中，煤層氣的解吸率主要受地層壓力、煤階和原始含氣量影響，頁巖氣的解吸率主要受泥質含量和頁巖頁理發育程度影響[4]。

美國是頁巖氣勘探與開發技術相對較高的國家，2009年美國的頁巖氣產量占同年美國的天然氣總產量的13%～15%。我國對頁巖氣的研究還處于探索階段,勘探資料表明我國寒武系、志留系、二疊系的頁巖地層十分發育，四川盆地、鄂爾多斯盆地、塔里木盆地、準噶爾盆地等具備頁巖氣成藏的基本條件，初步估計我國頁巖氣資源儲量約為(15～30)×1012m3，具有廣闊的開發前景[5-7]。

頁巖氣藏中吸附氣存在吸附解吸動態平衡，游離氣具有濃度擴散和達西滲流特征，張金川[8]全面分析頁巖氣成藏機理和氣體相態轉換；段永剛[9]根據頁巖氣吸附解吸動態平衡，建立了頁巖氣雙重介質壓裂井滲流數學模型以及無限導流壓裂井評價模型；Gregory R.King[10]推導出泥盆系頁巖氣藏物質平衡方程。這些文獻大多考慮了頁巖氣吸附解吸動態平衡，而未考慮解吸氣擴散特征。下面筆者結合滲流理論，應用Langmuir等溫吸附方程、Fick擴散定律等建立了頁巖氣在地層中的滲流數學模型。

1 吸附氣解吸擴散方程

在頁巖氣藏中，吸附氣的含量占總含量的20%～85%，具有吸附解吸和擴散特征，是頁巖氣與常規天然氣的最主要的區別。Langmuir在1916年從動力學角度推導出單分子層等溫吸附狀態方程。Langmuir等溫吸附曲線能夠表示吸附氣含量與游離氣含量處于平衡狀態的關系[11]：

(1)

式中，Ve為等溫吸附量，m3/m3；VL為飽和吸附氣含量，m3/m3；PL為吸附氣含量達到最大吸附量的50%壓力，MPa；Pg為氣體壓力，MPa。

吸附氣體發生解吸作用后，基巖內部和表面氣體存在濃度差，內部氣體以擴散方式向外運移。根據Fick擴散定律，吸附氣平均含量變化率與吸附氣平均含量和表面吸附氣含量之差成正比，基質巖塊擴散量與基質中吸附氣平均含量變化率成正比[12]：

(2)

式中，Vm為吸附氣平均含量，m3/m3；Dm為擴散系數；Fs為形狀因子；τ為吸附時間常數；qm為擴散量，m3/(m3·d)；FG為幾何因子。

2 頁巖氣滲流數學模型

頁巖氣由儲層流入井筒包括解吸、擴散、滲流3個過程[13-15]，基本假設條件如下：①游離氣以層流形式流向井筒，服從Darcy定律；吸附氣解吸，服從Langmuir等溫吸附方程；解吸氣擴散運移，服從Fick定律；②巖石、流體均可壓縮；③考慮氣藏為均質地層，各向同性；④忽略重力作用；⑤不考慮開采過程中頁巖氣藏溫度變化，滲流為等溫過程；⑥氣體粘度為常數。

圖1 平面徑向流模型

平面徑向流對應的物理模型如圖1所示，地層為水平圓盤狀，等厚，其厚度為h，dr為變化量，m；r為滲流半徑，m。

(3)

頁巖氣藏的孔隙度和滲透率極低，氣體在孔隙中滲流時，由于滲流阻力的作用，必然伴隨著能量的損失。法國工程師達西經過大量的試驗研究總結出的流量與壓力差、巖心截面、滲透率、粘度之間的相互關系即為達西定律。根據達西定律可以得出氣體運動方程：

(4)

式中，ρg為氣體密度，kg/m3；ρgsc為天然氣密度，kg/m3；vg為氣體滲流速度，m/s；kg為氣體滲透率，μm2；μg為氣體粘度，mPa·s。

滲流是一個運動過程，氣體狀態和巖石孔隙也隨之變化，頁巖氣不是理想氣體，不服從理想氣體定律，在工程上引入天然氣壓縮因子。隨著開采持續，地層壓力發生變化，孔隙度也隨之改變，用壓縮系數來表示：

(5)

式中，Cg為氣體彈性壓縮系數，MPa-1；Vg為氣體絕對體積，m3；z為壓縮因子；φ為孔隙度；Cf為孔隙壓縮系數，MPa-1；Ct為氣相綜合壓縮系數，MPa-1。

將上述推導的解吸方程(1)、擴散方程(2)、運動方程(4)和狀態方程(5)帶入方程(3)可得頁巖氣滲流數學模型：

內邊界條件：

外邊界條件：

式中，S為表皮系數；C為井筒儲集系數，m3/MPa。

3 井底壓力變化曲線

應用有限差分法對數學模型進行數值求解，繪制出井底壓力變化曲線圖。影響井底壓力變化的因素很多，筆者著重討論了吸附氣解吸擴散、表皮系數、Langmuir壓力、Langmuir體積對井底壓力變化的影響。圖2反映了考慮吸附氣解吸擴散與不考慮解吸擴散對井底壓力變化的影響，開采初期，壓力波在井口附近，吸附氣的解吸作用表現得不明顯，隨著開采時間增加，氣藏整體壓力降低，吸附氣解吸量增加，部分氣藏壓力損失得到補償，考慮解吸擴散的井底壓力導數曲線在徑向流前期偏離直線段，且較晚探測到地層邊界。

圖3反映了不同表皮系數對井底壓力變化的影響。隨著表皮系數的增大，井底壓力差變大，井底壓力降低較快，壓力導數“拱”的弧度增大，可見要保證氣井長期穩產，必須降低表皮系數的影響。

圖2 解吸擴散對井底壓力的影響 圖3 表皮系數對井底壓力的影響

圖4反映了Langmuir壓力(PL)對井底壓力變化的影響。由圖4可以看出，當Langmuir壓力接近原始地層壓力時，井底壓力降低最慢；當Langmuir壓力遠離地層壓力時，井底壓力降低加快。

圖5反映了Langmuir體積(VL)對井底壓力變化的影響。由圖5可以看出，隨著Langmuir體積的不斷增大，氣藏中的吸附氣的解吸量增大，補償井底壓力損失明顯，井底壓力降低緩慢。

圖4 Langmuir壓力對井底壓力的影響 圖5 Langmuir體積對井底壓力的影響

4 結 語

頁巖氣是一種重要的非常規油氣資源，高效開發頁巖氣藏將會對天然氣市場產生重大影響。當前，對頁巖氣藏的數值模擬多為單相模擬，對流體在孔隙中的滲流機理及表征方法多采用煤層氣的研究成果。因此，對頁巖氣藏中的吸附解吸機制和滲流機理需要加大物理模擬實驗，對氣、水兩相的頁巖氣藏數值模擬需要進一步研究和改善。

[1]張金川，薛會，王艷芳，等.中國非常規天然氣勘探雛議[J].天然氣工業，2006,26(12)：53-56.

[2]蒲泊伶，包書景，王毅，等.頁巖氣成藏條件分析——以美國頁巖氣盆地為例[J].石油地質與工程，2008,22(3)：33-36.

[3]江山，王新海，鄭愛玲.煤層氣開采的二維非平衡吸附模型[J].天然氣工業，2005,25(6)：81-82.

[4]張利萍，潘仁芳.頁巖氣的主要成藏要素與氣儲改造[J].中國石油勘探，2009,3：20-23.

[5]胡文瑞，翟光明，李景明.中國非常規油氣的潛力和發展[J].中國工程科學，2010,12(5)：25-29.

[6]張金川，徐波，聶海寬，等.中國頁巖氣資源勘探潛力[J].天然氣工業，2008,28(6)：136-140.

[7]劉德華，肖佳林，關富佳.頁巖氣開發技術現狀及研究方向[J].石油天然氣學報(江漢石油學院學報)，2011,33(1)：119-123.

[8]張金川，金之鈞，袁明生.頁巖氣成藏機理和分布[J].天然氣工業，2004,24(7)：15-18.

[9]段永剛，魏明強，李建秋，等.頁巖氣藏滲流機理及壓裂井產能評價[J].重慶大學學報，2011,34(4)：62-65.

[10]King G R.Material Balance Techniques for Coal seem and Devonian Shale Gas Reservoirs[J].SPE 20730,1990.

[11]張群.煤層氣儲層數值模擬模型及應用的研究[D].西安：煤炭科學研究總院西安分院，2003.

[12]胡玲，徐蕓蕓，吳瑾.Fick第二定律的應用研究現狀與展望[J].河北建筑科技學院學報，2005,22(3)：50-53.

[13]Sawyer W K，Paul G W.Development and application of a 3D coal-bed simulator[J].SPE 90119,1990.

[14]段永剛，李建秋.頁巖氣無限導流壓裂井壓力動態分析[J].天然氣工業，2010,30(10)：1-4.

[15]KIing G R，Turgay E,Fred C S.Numerical simulation of the transient behavior of coal seem degasification wells[J].SPE Formation Evaluation,1986,1(2)：165-183.

10.3969/j.issn.1673-1409(N).2012.08.022

TE331

A

1673-1409(2012)08-N068-04

2012-03-23

國家油氣重大專項(2011ZX05015)。

尹虎(1987-)，男，2010年大學畢業，碩士生，現主要從事試井、油藏工程方面的研究工作。

[編輯] 洪云飛