基于T-P模型的高低溫環(huán)境型煤對(duì)甲烷吸附性能預(yù)測(cè)*

岳基偉，王兆豐,2

(1.河南理工大學(xué) 安全科學(xué)與工程學(xué)院，河南 焦作 454000;2.煤礦災(zāi)害預(yù)防與搶險(xiǎn)救災(zāi)教育部工程研究中心，河南 焦作 454000)

0　引言

吸附等溫線是研究煤儲(chǔ)存甲烷性能的重要參數(shù)之一，其對(duì)于計(jì)算煤層瓦斯含量至關(guān)重要[1]。測(cè)定煤層瓦斯含量的方法有直接法和間接法2種，直接法一般應(yīng)用于取樣深度小于1 000 m以下的煤層，對(duì)于取樣深度大于1 000 m以上的煤層只能以間接法獲得，大多采用Langmuir模型預(yù)測(cè)煤對(duì)甲烷吸附等溫線，但只考慮了壓力的影響，未考慮溫度的影響[2]，Langmuir模型只能獲得該溫度下的等溫吸附曲線，而不能獲得其他溫度下的等溫吸附曲線[3]。

在煤礦瓦斯治理中，隨著開采深度、地溫及煤層瓦斯壓力的增加，導(dǎo)致深部煤樣較難取得，其吸附等溫線也較難測(cè)得[4]，若能以少量等溫吸附試驗(yàn)數(shù)據(jù)，精確預(yù)測(cè)其他溫度和壓力下的吸附量，則對(duì)減少實(shí)驗(yàn)工作量大有益處[5]。姜偉等[6]、蘇現(xiàn)波等[7]、藺亞兵等[8]、柴琳等[9]、馬向攀等[10]以某幾個(gè)溫度下等溫吸附數(shù)據(jù)預(yù)測(cè)其他溫度下的等溫吸附數(shù)據(jù)，其模型中的參數(shù)為經(jīng)驗(yàn)參數(shù)，且采用的顆粒煤，取芯管取出的煤樣一般是塊狀煤，更符合煤的原始賦存狀態(tài)，采用顆粒煤不能較好的反映煤層瓦斯的解吸規(guī)律[11]。由于取芯管取樣，鉆桿與孔壁之間的摩擦，造成煤樣溫度升高、瓦斯解吸加快[12-14]，王兆豐等[15]提出低溫取芯的方法，低溫可以抑制瓦斯解吸，降低取芯過程中的損失量，但其對(duì)于50 ℃以下的低溫環(huán)境較難達(dá)到；岳高偉[4]等通過采用吸附熱理論對(duì)顆粒煤高低溫環(huán)境下的吸附等溫線進(jìn)行了預(yù)測(cè)，但未對(duì)型煤的吸附性能進(jìn)行研究，且公式中參數(shù)的選取只是經(jīng)驗(yàn)參數(shù)。

因此，本文基于T-P模型，采用型煤分別在293.15，273.15，253.15 K環(huán)境下進(jìn)行等溫吸附實(shí)驗(yàn)，并對(duì)公式中的參數(shù)進(jìn)行了合理的求解，同時(shí)也驗(yàn)證了T-P模型預(yù)測(cè)高低溫環(huán)境下吸附等溫線的準(zhǔn)確性。

1　煤樣的制備及吸附特征曲線

1.1　煤樣的制備

試驗(yàn)煤樣取自山西興無煤礦，其煤質(zhì)為焦煤，將取得的新鮮煤樣制成0.17～0.25 mm煤樣。型煤制作流程如下：

1)將制備好的0.17～0.25 mm的煤樣均勻混合，添加20%的蒸餾水，充分?jǐn)嚢琛?/p>

2)把充分?jǐn)嚢韬蟮臐駶櫭簶蛹尤胄兔耗＞咧校捎肊HC-3100型微機(jī)控制電液伺服萬能機(jī)進(jìn)行型煤的壓制，設(shè)置壓制負(fù)荷為60 kN，并穩(wěn)壓20 min。

3)安全緩慢的退出型煤。

1.2　試驗(yàn)系統(tǒng)

實(shí)驗(yàn)采用河南理工大學(xué)瓦斯裝備與技術(shù)研究所自制的高低溫試驗(yàn)裝置，如圖1所示。裝置有高低溫變頻控制系統(tǒng)、吸附平衡單元、真空脫氣單元、解吸測(cè)量單元、數(shù)據(jù)實(shí)時(shí)采集系統(tǒng)組成。溫度測(cè)量范圍：-50～100℃；溫度偏差：≤±1℃；溫度波動(dòng)度：±0.5℃。

圖1　高/低溫吸附裝置示意Fig.1　Diagram of high/low temperature adsorption device

1.3　試驗(yàn)方法

1)將型煤煤樣烘干；2)對(duì)實(shí)驗(yàn)裝置進(jìn)行抽真空，抽至20 Pa以下；3)將煤樣罐放置在高低溫變頻控制系統(tǒng)中，根據(jù)實(shí)驗(yàn)需要將高低溫變頻控制系統(tǒng)分別設(shè)置為253.15，273.15，293.15 K，采用吸附平衡單元分別設(shè)定不同的吸附平衡壓力，并測(cè)試吸附量。

1.4　吸附試驗(yàn)結(jié)果及分析

分別測(cè)試293.15 ，273.15 ，253.15 K煤的瓦斯吸附量，其吸附等溫線，如圖2所示。

圖2　不同溫度下的煤的瓦斯吸附等溫線Fig.2　Adsorption isotherm at different temperature

由圖2可知，同一吸附平衡壓力下，溫度越低，煤的瓦斯吸附量越大。因?yàn)闇囟仍降?，甲烷分子的能量越低，其運(yùn)動(dòng)越慢，游離的甲烷分析從而更容易被吸附。已經(jīng)處于吸附狀態(tài)的甲烷分子更難獲得能量克服煤與甲烷分子之間的作用力變成游離態(tài)。

1.5　煤對(duì)甲烷的ε-ω吸附特征曲線

根據(jù)Polanyi的吸附勢(shì)理論ε：

ε=RTln(ps/p)

(1)

式中：ε為吸附勢(shì)，kJ·mol-1；ps為飽和蒸氣壓，MPa；p為吸附平衡壓力；R為氣體常數(shù)，J/(mol·K)。

飽和蒸氣壓ps可由Amankwah[16]改進(jìn)的計(jì)算氣體飽和蒸氣壓的公式得出：

(2)

式中：pc為甲烷的臨界壓力，4.62 MPa；Tc為甲烷臨界溫度，190.6 K；m為與吸附體系有關(guān)的參數(shù)。

吸附空間ω：

(3)

ρa(bǔ)d=ρbexp[-0.002 5(T-Tb)]

(4)

式中：ρb為沸點(diǎn)下甲烷的吸附相密度，0.424 g/cm3; T為實(shí)驗(yàn)溫度，K；Tb為甲烷在沸點(diǎn)時(shí)的溫度，111.5 K。

以m=3為例，采用293.15 ，273.15，253.15 K的吸附數(shù)據(jù)，繪制ε-ω吸附特征曲線，如圖3所示。

圖3　ε-ω吸附特征曲線Fig.3　ε-ω adsorption characteristic curve

由圖中可以看出，293.15 ，273.15，253.15 K的ε-ω吸附特征曲線基本上在一條曲線上，因此ε-ω吸附特征曲線與溫度無關(guān)。吸附勢(shì)能ε和吸附空間ω滿足ε=-alnω+b的關(guān)系，ε-ω吸附特征曲線呈現(xiàn)對(duì)數(shù)的形式。

2　合理m值的確定

m為與吸附體系有關(guān)的參數(shù)，為使m的值較為準(zhǔn)確，分別繪出m值為1.5，1.7，2，2.2，2.3，2.4，2.5，2.6，2.7，2.8，2.9，3.5，3.7，4，4.2，4.5，5，5.2，5.3，5.4，5.5及6時(shí)型煤的ε-ω吸附特征曲線，并對(duì)曲線進(jìn)行擬合分析，m與擬合度(R2)的曲線圖，如圖4所示。

圖4　m與擬合度(R2)的曲線Fig.4　The curve of m and the degree of fitting (R2)

由圖4可知，擬合度R2與m滿足拋物線關(guān)系，即R2=cm2+dm+e(c<0)，其擬合效果較好。因?yàn)閏<0，R2=cm2+dm+e有最大值。其最大值為R2(-d/2c)=(4ce-d2)/4c，因此，可求出擬合度R2的最大值及取最大值時(shí)所對(duì)應(yīng)的m值，R2(3.59)max=0.968。

3　T-P(溫度-壓力)綜合模型的推導(dǎo)及吸附等溫線的預(yù)測(cè)

3.1　T-P(溫度-壓力)綜合模型的推導(dǎo)

由式ε=-alnω+b可知，ω與吸附量Q滿足線性關(guān)系，因此，式ε=-alnω+b可變換為ε=-flnQ+g。

可得：lnQ=(g-ε)/f，Q=eg/f/eε/f，令eg/f=A，Q=A/eRTln(pc(T/Tc)m/p)/f，對(duì)其取對(duì)數(shù)lnQ=lnA-RTln(pc(T/Tc)m/p)/f，lnQ=lnA+RTlnp/f-RT/fln(pc(T/Tc)m)，令M=RT/fN=lnA-RT/fln(pc(T/Tc)m)，則lnQ=Mlnp+N。

采用溫度為283.15 K吸附等溫線計(jì)算特征常數(shù)M和N，如圖5所示。

圖5　lnQ與lnp的關(guān)系曲線Fig.5　The relation curve between lnQ and lnp

由圖5可知：M=0.360 74，N=2.477 96。因此可得：A=34.557 3，f=6.525 8，因此，吸附量Q關(guān)于T-P綜合模型Q=34.557 3/eRTln(pc(T/Tc)3.59/p)/6.525 8

3.2　吸附等溫線的預(yù)測(cè)

根據(jù)吸附量Q關(guān)于T-P綜合模型Q=34.557 3/eRTln(pc(T/Tc)3.59/p)/6.525 8，預(yù)測(cè)溫度為263.15 和303.15 K的吸附等溫線，并把預(yù)測(cè)值與實(shí)測(cè)值進(jìn)行對(duì)比，如圖6所示。

圖6　吸附等溫線的實(shí)測(cè)值與預(yù)測(cè)值對(duì)比Fig.6　Compassion measured value of isothermal adsorption curve to predicted value

由圖6可知，實(shí)測(cè)值與預(yù)測(cè)值差別不大，是一致的。其相對(duì)誤差不超過5%。采用T-P綜合模型Q=34.557 3/eRTln(pc(T/Tc)3.59/p)/6.525 8預(yù)測(cè)其他溫度下的吸附等溫線。

采用T-P模型預(yù)測(cè)溫度為333.15，323.15，313.15，293.15，283.15，273.15，253.15，243.15，233.15，223.15，213.15 K的吸附等溫線，如圖7所示。

圖7　型煤吸附甲烷等溫線預(yù)測(cè)Fig.7　Prediction of isothermal adsorption curve of methane in coal briquette

4　結(jié)論

1)同一吸附平衡壓力下，溫度越低，煤的瓦斯吸附量越大。

2)ε-ω吸附特征曲線與溫度無關(guān)，其呈現(xiàn)對(duì)數(shù)曲線的形式，擬合效果較好。

3)與吸附體系有關(guān)的參數(shù)m和擬合度R2滿足拋物線的關(guān)系，可得到擬合效果最好時(shí)的參數(shù)m值。

4)采用T-P模型預(yù)測(cè)得到的吸附等溫線與實(shí)測(cè)的吸附等溫線，差別不大，其相對(duì)誤差不超過5%。

[1]蘇現(xiàn)波，林萌，林曉英，等. 吸附勢(shì)理論在煤層甲烷吸附中的應(yīng)用[J]. 中國煤層氣，2006，3(2): 28-30.

SU Xianbo, LIN Meng, LIN Xiaoying, et al. Study of methane adsorption on coal with adsorption potential theory[J]. China Coal Bed Methane, 2006, 3(2): 28-30.

[2]鐘玲文，鄭玉柱，員爭榮，等. 煤在溫度和壓力綜合影響下的吸附性能及氣含量預(yù)測(cè)[J]. 煤炭學(xué)報(bào)，2002，27(6):581-585.

ZHONG Lingwen, ZHENG Yuzhu, YUAN Zhengrong, et al. The adsorption capability of coal under integrated influence of temperature and pressure and predicted of content quantity of coal bed gas [J]. Journal of China Coal Society, 2002, 27(6):581-585.

[3]趙振國.吸附作用應(yīng)用原理[M].北京:化學(xué)工業(yè)出版社, 2005.

[4]岳高偉，王兆豐，康博，等. 基于吸附熱理論的煤—甲烷高低溫等溫吸附線預(yù)測(cè)[J]. 天然氣地球科學(xué)，2015，26(1)：148-153.

YUE Gaowei, WANG Zhaofeng, KANG Bo, et al. Prediction for isothermal adsorption curve of coal/ CH4based on adsorption heat theory[J]. Natural Gas Geoscience, 2015, 26(1):148-153.

[5]陳紹杰，陳學(xué)習(xí)，柏松，等. 基于吸附勢(shì)理論的煤—甲烷吸附等溫線預(yù)測(cè)[J]. 華北科技學(xué)院學(xué)報(bào)，2009，6(4)：30-32.

CHEN Shaojie, CHEN Xuexi, BAI Song, et al. Prediction for isothermal adsorption curve of coal/CH4based on adsorption potential theory[J]. Journal of North China Institute of Science and Technology, 2009, 6(4): 30-32.

[6]姜偉，吳財(cái)芳，姜瑋，等. 吸附勢(shì)理論在煤層氣吸附解吸研究中的應(yīng)用[J]. 煤炭科學(xué)技術(shù)，2011，39(5)：102-104.

JIANG Wei, WU Caifang, JIANG Wei, et al. Application of adsorption potential theory to study on adsorption-desorption of coalbed methane[J]. Coal Science and Technology, 2011, 39(5): 102-104.

[7]蘇現(xiàn)波，林曉英，趙孟軍，等. 儲(chǔ)層條件下煤吸附甲烷能力預(yù)測(cè)[J]. 天然氣工業(yè)，2006，26(8):34-36.

SU Xianbo, LIN Xiaoying, ZHAO Mengjun, et al. Prediction on coal adsorption capacity under reservoir conditions[J]. Natural Gas Industry, 2006, 26(8): 34-36.

[8]藺亞兵，馬東民，劉鈺輝，等. 溫度對(duì)煤吸附甲烷的影響實(shí)驗(yàn)[J]. 煤田地質(zhì)與勘探，2012，40(6)：24-28.

LIN Yabing, MA Dongmin, LIU Yuhui, et al. Experiment of the influence of temperature on coal bed methane adsorption[J]. Coal Geology and Exploration, 2012, 40(6): 24-28.

[9]柴琳，吳世躍，牛熠. 基于吸附勢(shì)理論的煤吸附超臨界甲烷研究[J]. 煤礦安全，2017，48(1)：21-23，27.

CHAI Lin，WU Shiyue，NIU Yu. Study on coal adsorbing supercritical methane based on adsorption potential theory[J]. Safety in Coal Mines, 2017, 48(1):21-23, 27.

[10]馬向攀，王兆豐，任浩洋. 基于吸附勢(shì)理論的低溫環(huán)境煤甲烷吸附等溫線預(yù)測(cè)[J]. 煤礦安全，2017，48(5)：22-25.

MA Xiangpan, WANG Zhaofeng, REN Haoyang. Prediction for CH4isothermal adsorption curve of coal at low temperature based on adsorption potential theory[J]. Safety in Coal Mines, 2017, 48(5):22-25.

[11]胡千庭，鄒銀輝，文光才，等．瓦斯含量法預(yù)測(cè)突出危險(xiǎn)新技術(shù)[J]．煤炭學(xué)報(bào)，2007，32(3):276-280.

HU Qianting, ZOU Yinhui, WEN Guangcai, et al. New technology of outburst danger prediction by gas content[J]. Journal of China Coal Society, 2007, 32(3):276-280.

[12]TANG Xu, LI Zhiqiang, RIPEPI NINO, et al. Temperature-dependent diffusion process of methane through dry crushed coal[J]. Journal of Natural Gas Science and Engineering, 2015, 22:609-617.

[13]李志強(qiáng)，成強(qiáng)，段正鵬，等. 不同溫壓柱狀煤芯瓦斯吸附飽和度和進(jìn)擴(kuò)散時(shí)間的確定[J].中國安全生產(chǎn)科學(xué)技術(shù)，2017，13(7)：92-99.

LI Zhiqiang, CHENG Qiang, DUAN Zhengpeng, et al. Determination of adsorption saturation and diffusion time of gas in cylindrical coal core under different temperatures and pressures[J]. Journal of Safety Science and Technology, 2017, 13(7):92-99.

[14]聶百勝，楊濤，李祥春，等. 煤粒瓦斯解吸擴(kuò)散規(guī)律實(shí)驗(yàn)[J]. 中國礦業(yè)大學(xué)學(xué)報(bào)，2013，42(6):975-980.

NIE Baisheng, YANG Tao, LI Xiangchun, et al. Reaearch on diffusion of methane in coal particles[J]. Journal of China University of Mining and Technology, 2013, 42(6):975-980.

[15]王兆豐，岳高偉，康博，等．低溫環(huán)境對(duì)煤的瓦斯解吸抑制效應(yīng)試驗(yàn)[J]．重慶大學(xué)學(xué)報(bào)，2014，37(9):106-112，143．

WANG Zhaofeng, YUE Gaowei, KANG Bo, et al. Gas desorption inhibitory effect of coal in low temperature environment[J]. Journal of Chongqing University, 2014, 37(9): 106-112, 143．

[16]AMANKWAH KAG, SCHWAKZ J A. A modified approach for estimating pseudo-vapor pressures in the application of the Dubinin-Astakhov equation[J]. Carbon, 1995, 33(9):1313-1319.