楊　龍，梅海燕，張茂林，袁　恩

考慮干酪根中溶解氣的頁巖氣藏儲量計算方法

楊龍，梅海燕，張茂林，袁恩

（西南石油大學(xué)石油與天然氣工程學(xué)院，成都610500）

前人在計算頁巖氣儲量時未考慮干酪根中溶解的頁巖氣，但已有學(xué)者證明，干酪根中也溶解大量的頁巖氣，在進行儲量計算時不能忽略。不考慮干酪根中溶解氣的儲量計算方法，不能真實、準確地計算頁巖氣藏的儲量。在前人關(guān)于頁巖氣儲量計算方法的基礎(chǔ)上，建立了同時考慮吸附相密度、基質(zhì)和裂縫孔隙體積隨壓力變化、吸附相體積隨壓力變化以及儲存在干酪根中溶解氣的物質(zhì)平衡方程。通過實例計算可知，溶解氣占據(jù)總儲量的4.69%左右，不考慮溶解氣計算出來的自由氣儲量、吸附氣儲量和總儲量都偏大。因此，為了確保儲量計算的準確性，在進行儲量計算時必須將溶解氣考慮在內(nèi)。

干酪根；頁巖氣；基質(zhì)；裂縫；溶解氣；孔隙度；物質(zhì)平衡方程；儲量計算

頁巖氣是一種連續(xù)生成的生物化學(xué)作用氣、熱裂解作用氣或兩者混合而成的天然氣，具有運移距離短、多種封閉機理、成藏隱蔽、地層飽和氣等成藏特點［1］。頁巖氣賦存方式有3種：第一種是以游離狀態(tài)存在于頁巖孔隙和天然裂縫中；第二種是以吸附狀態(tài)存在干酪根和黏土顆粒表面；第三種是以溶解狀態(tài)存在干酪根和瀝青質(zhì)中。前人研究認為，溶解氣含量較少，進行儲量計算時可以不考慮溶解氣的影響［2-4］。而文獻［5］通過實驗發(fā)現(xiàn)，干酪根中的溶解氣也是氣體的重要儲存方式，忽略溶解氣會導(dǎo)致儲量計算結(jié)果不符合實際。

在頁巖氣資源勘探開發(fā)過程中，儲量評價十分重要［6］，物質(zhì)平衡法是目前評價頁巖氣藏儲量的常用方法［7］。文獻［7］考慮吸附氣的影響建立了頁巖氣藏的物質(zhì)平衡方程，提出Z*計算方法；文獻［8］重新定義了Z因子，簡化了頁巖氣藏物質(zhì)平衡方程，但是卻沒有考慮裂縫系統(tǒng)；文獻［2］同時考慮基質(zhì)孔隙體積和裂縫孔隙體積隨地層壓力變化，建立了封閉性頁巖氣藏物質(zhì)平衡方程；文獻［3］建立了同時考慮有效應(yīng)力和基質(zhì)收縮對孔隙度影響的物質(zhì)平衡方程；文獻［4］在計算游離氣體積時，考慮了吸附相所占有的孔隙體積；文獻［9］建立了考慮吸附相密度變化的物質(zhì)平衡方程；文獻［10］建立了同時考慮基質(zhì)孔隙體積和裂縫孔隙體積隨地層壓力變化，以及吸附氣所占體積的物質(zhì)平衡方程。但是，上述方法均沒有考慮干酪根中的溶解氣對儲量計算的影響，本文在前人的研究基礎(chǔ)上，提出了考慮干酪根中溶解氣的頁巖氣藏物質(zhì)平衡方程。

1　頁巖氣藏物質(zhì)平衡方程的建立

由于頁巖氣和煤層氣的吸附機理相似，因此采用Langmuir吸附等溫式作為頁巖氣藏中吸附氣量的表達式。假設(shè)頁巖氣藏開發(fā)處于等溫狀態(tài)，則頁巖氣等溫吸附式［11］為

基質(zhì)總體積為

地面條件下吸附氣體積為

由質(zhì)量守恒定律可得，ρ1V1=ρ2V2.因此，地層條件下吸附相的體積為

吸附相的體積變化量為

吸附氣的解吸量為

隨著地層壓力的下降，巖石顆粒變形和束縛水發(fā)生彈性膨脹，則基質(zhì)和裂縫中孔隙體積的變化值為

瀝青質(zhì)中甲烷溶解度的經(jīng)驗公式可由（9）式來計算，為了避免在壓力較低時溶解度出現(xiàn)負值，在地層溫度為25～100℃、壓力為0時，溶解度也為0。

式中b1=-0.018 931；b2=-0.850 480；

b3=827.260；b4=-635.260.

甲烷在干酪根中的溶解度和在瀝青質(zhì)中一樣，則溶解氣體積的計算公式為［13-14］

又因為總干酪根占基質(zhì)的體積分數(shù)Vtker可以表示為［15］

因此，干酪根的總體積

Vtker的計算式為［16］

則溶解氣的儲量為

當(dāng)?shù)貙訅毫ο陆档絧時，在地層條件下，溶解氣的擴散量可由下式計算：

根據(jù)體積守恒定律可得到

累計產(chǎn)氣量=基質(zhì)內(nèi)自由氣膨脹體積+顆粒巖石變形和束縛水的膨脹量+解吸的吸附氣體積+擴散的溶解氣體積-吸附相變化體積+裂縫內(nèi)自由氣的膨脹體積+裂縫內(nèi)孔隙體積減小和束縛水膨脹體積

因此可以得到：

令，代入（18）式，可以簡化得

2　模型求解

令Y=VPBg，M=Bg-Bgi+CcfBgiΔp，

整理（19）式可以得到

兩邊同時除以M可以得到

3　實例分析

某頁巖氣藏的基本參數(shù)如下：pi=24.13 MPa，Bgi=4.82×10-3m3/m3，Smwi=0.25，Cm=4.35×10-4MPa-1，Cf= 2.0×10-2MPa-1，Cw=4.35×10-4MPa-1，Sfwi=0，φmt=0.112，VL= 11.32m3/t，pL=2.41MPa，ρsc=7.7×10-4g/cm3，ρs=0.34g/cm3，ρb=2.65 g/cm3，T=366.49 K，CTOC=5%，φorg=0.012，φads= 0.002 3，ρko=1.325 g/cm3，生產(chǎn)數(shù)據(jù)如表1.

表1　某頁巖氣藏生產(chǎn)數(shù)據(jù)

利用上述數(shù)據(jù)，使用本文中推導(dǎo)的物質(zhì)平衡方程來計算頁巖氣藏儲量，結(jié)果見圖1和表2.

圖1　本文推導(dǎo)物質(zhì)平衡方程的儲量回歸

表2　不同方法計算儲量結(jié)果對比108m3

由計算結(jié)果可知，考慮溶解氣計算的基質(zhì)中自由氣的儲量比不考慮溶解氣的計算結(jié)果要少6.36%，裂縫中自由氣的儲量少9.26%，吸附氣儲量少6.36%，總儲量少2.30%，考慮溶解氣計算的溶解氣儲量為0.984 7×108m3，占據(jù)總儲量的4.69%.

忽略溶解氣的影響會導(dǎo)致累計產(chǎn)氣量中自由氣的貢獻率偏大，導(dǎo)致計算出來的自由氣、吸附氣儲量偏大，進而導(dǎo)致計算出來的總儲量偏大。由此可見，估算頁巖氣藏儲量的時候，考慮溶解氣會更加精確。

4　結(jié)論

本文建立同時考慮吸附相密度、基質(zhì)和裂縫中孔隙體積隨壓力變化、吸附相體積隨壓力變化和儲存在干酪根中溶解氣等多因素影響的物質(zhì)平衡方程。

通過實例計算發(fā)現(xiàn)，忽略溶解氣的影響會導(dǎo)致累計產(chǎn)氣量中自由氣的貢獻率偏大，導(dǎo)致計算出來的自由氣、吸附氣儲量偏大，進而導(dǎo)致計算出來的總儲量偏大。而且，溶解氣的儲量占據(jù)總儲量的4.69%，在進行儲量計算時不可忽略溶解氣的影響。

符號注釋

Bg——當(dāng)前地層壓力下氣體的體積系數(shù)，m3/m3；

Bgi——頁巖氣的體積系數(shù)，m3/m3；

Cf——裂縫壓縮系數(shù)，MPa-1；

Cm——基質(zhì)巖石的壓縮系數(shù)，MPa-1；

Cw——地層水壓縮系數(shù)，MPa-1；

CTOC——總有機碳含量，%；

c（p）——干酪根中甲烷的溶解度，m3/m3；

Va——地層條件下吸附相體積，108m3；

Vasc——地面條件下吸附氣體積，108m3；

Vdsc（p）——在地面條件下溶解氣的體積，10m；

Vdsk（p）——干酪根中溶解氣體積，108m3；

Vf——裂縫中游離氣儲量，108m3；

Vm——基質(zhì)中自由氣的體積，108m3；

ΔVa——解吸氣在地層條件下的體積，108m3；

ΔVd——溶解氣的擴散量在地層條件下的體積，108m3；

ΔVf——裂縫中孔隙體積的變化量，108m3；

ΔVm——基質(zhì)中孔隙體積的變化量，108m3；

p——氣藏壓力，MPa；

pi——氣藏原始壓力，MPa；

pL——蘭氏壓力，MPa；

Sfwi——裂縫中束縛水飽和度，f；

Smwi——基質(zhì)中束縛水飽和度，f；

T——氣藏溫度，K；

Vdiff——固體干酪根占基質(zhì)總體積的體積分數(shù)，f；

VE——地層壓力p下的等溫吸附量，m3/t；

VL——蘭氏體積，m3/t；

Vrock——基質(zhì)總體積，108m3；

Vsk——干酪根總體積，108m3；

Vtker——總干酪根占氣藏總體積的體積分數(shù)，f；

ΔVa——吸附相體積的變化量，108m3；

φads——吸附相的孔隙度，f；

φorg——有機質(zhì)中自由氣所占據(jù)的孔隙度，f；

φmt——基質(zhì)孔隙度，f；

ρb——巖石密度，g/cm3；

ρko——干酪根密度，g/cm3；

ρr——干酪根的相對密度，無量綱；

ρs——地層條件下的吸附相密度，g/cm3；

ρsc——地面條件下的頁巖氣密度，g/cm3.

［1］張金川，薛會，張德明，等.頁巖氣及其成藏機理［J］.現(xiàn)代地質(zhì)，2003，17（4）：466. ZHANG Jinchuan，XUE Hui，ZHANG Deming，et al.The accumulation mechanism for shale gas reservoir［J］.Geoscience，2003，17（4）：466.

［2］劉鐵成，唐海，劉鵬超，等.裂縫性封閉頁巖氣藏物質(zhì)平衡方程及儲量計算方法［J］.天然氣勘探與開發(fā)，2011，34（2）：28-30. LIU Tiecheng，TANG Hai，LIU Pengchao，et al.Material balance equation and reserve calculation method of fractured and closed shale-gas reservoir［J］.Natural Gas Exploration and Development，2011，34（2）：28-30.

［3］劉波濤，尹虎，王新海，等.修正巖石壓縮系數(shù)的頁巖氣藏物質(zhì)平衡方程及儲量計算［J］.石油與天然氣地質(zhì)，2013，34（4）：471-474. LIU Botao，YIN Hu，WANG Xinhai，et al.Material balance equation with revised rock compressibility for shale gas reserve calculation［J］.Oil&Gas Geology，2013，34（4）：471-474.

［4］楊浩瓏，戚志林，李龍，等.頁巖氣儲量計算的新物質(zhì)平衡方程［J］.油氣田地面工程，2013，32（8）：1-3. YANG Haolong，QI Zhilin，LI Long，et al.New material balance equation for shale gas reserves calculation［J］.Oil-Gas Field Surface Engineering，2013，32（8）：1-3.

［5］AGUILERA R，LOPEZ B.Evaluation of quintuple porosity in shale petroleum reservoirs［R］.SPE 165681，2013.

［6］曲占慶，林珊珊，張杰，等.多組分和吸附對頁巖氣儲量計算的影響［J］.特種油氣藏，2012，19（3）：114-116. QU Zhanqing，LIN Shanshan，ZHANG Jie，et al.Impacts of multicomponent and adsorption on shale gas reserve estimation［J］.Special Oil&Gas Reservoirs，2012，19（3）：114-116.

［7］KING G R.Material-balance techniques for coal-seam and Devonian shale gas reservoirs with limited water influx［J］.SPE Reservoir Engineering，1993，8（1）：67-72.

［8］MOGHADAM S，JEJE O，MATTAR L.Advanced gas material balance in simplified format［J］.Journal of Canadian Petroleum Technology，2011，50（1）：90-98.

［9］張烈輝，陳果，趙玉龍，等.改進的頁巖氣藏物質(zhì)平衡方程及儲量計算方法［J］.天然氣工業(yè)，2013，33（12）：66-70. ZHANG Liehui，CHEN Guo，ZHAO Yulong，et al.A modified material equation for shale gas reservoirs and a calculation method of shale gas reserves［J］.Natural Gas Industry，2013，33（12）：66-70.

［10］熊鈺，熊萬里，劉啟國，等.考慮吸附相體積的頁巖氣儲量計算方法［J］.地質(zhì)科技情報，2015，34（4）：139-143. XIONG Yu，XIONG Wanli，LIU Qiguo，et al.Calculation of shale gas reserves based on the volume of adsorbed gas［J］.Geological Science and Technology Information，2015，34（4）：139-143.

［11］傅獻彩，沈文霞，姚天揚，等.物理化學(xué)［M］.北京：高等教育出版社，2006：360-375. FU Xiancai，SHEN Wenxia，YAO Tianyang，et al.Physical Chemistry［M］.Beijing：Higher Education Press，2006：360-375.

［12］MEHROTRA A K，SVRCEK W Y.Correlations for properties of bitumen saturated with CO2，CH4and N2，and experiments with combustion gas mixtures［J］.Journal of Canadian Petroleum Technology，1982，21（6）：95-104.

［13］SWAMI V，SETTARI A T，JAVADPOUR F.A numerical model for multi-mechanism flow in shale gas reservoirs with application to laboratory scale testing［R］.SPE 164840，2013.

［14］SWAMI V，SETTARI A.A pore scale gas flow model for shale gas reservoir［R］.SPE 155756，2012.

［15］LOPEZ B，AGUILERA R.Petrophysical quantification of multiple porosities in shale petroleum reservoirs［R］.SPE 171638，2014.

［16］WU P，AGUILERA R.Investigation of gas shales at nanoscale using scan electron microscopy，transmission electron microscopy and atomic force microscopy［R］.SPE 159887，2012.

（編輯曹元婷）

A New Method to Calculate Shale Gas Reserves Considering Dissolved Gas in Kerogen

YANG Long，MEI Haiyan，ZHANG Maolin，YUAN En
（School of Petroleum and Natural Gas Engineering，Southwest Petroleum University，Chengdu，Sichuan 610500，China）

Predecessors did not consider the dissolved shale gas in kerogen when calculating the shale gas reserves.But it has been proved by some researchers that large amount of shale gas dissolved in kerogen can not be ignored in the calculation.The reserves calculation method ignoring the dissolved gas in kerogen can’t be used to truly and accurately calculate the reserves in shale gas reservoirs.Based on the previous calculation methods for shale gas reserves，a new material balance equation is established considering adsorbed phase density，volumetric variations of matrix and fracture pores with pressure，changes of adsorbed phase volume with pressure and dissolved gas in kerogen.Case study shows that dissolved gas accounts for about 4.69%of the total reserves.Reserves of free gas，adsorbed gas and total gas will be overestimated if dissolved gas is not taken into consideration.Therefore，in order to ensure the accuracy of reserves calculation，dissolved gas must be taken into account in the reserves calculation.

kerogen；shale gas；matrix；fracture；dissolved gas；porosity；material balance equation；reserves calculation

TE15

A

1001-3873（2016）05-0602-04DOI：10.7657/XJPG20160520

2016-03-08

2016-06-23

國家科技重大專項（2016ZX05027004，2016ZX05060006）

楊龍（1992-），男，湖北襄陽人，碩士研究生，油氣藏數(shù)值模擬與動態(tài)分析，（Tel）18190837618（E-mail）1137171129@ qq.com