Research on adsorption characteristics and adsorption model of CH4 and CO2 in Leping Formation coal seam of Fengcheng Mining Area
-
摘要:
煤体对气体的吸附特征决定了煤层的地质封存能力。对丰城矿区主力B4煤层进行了不同温度条件下的等温吸附实验,研究了煤层对气体的吸附特征,采用9种不同吸附数学模型对CO2吸附实验数据进行拟合,检验不同模型的拟合效果。结果表明:丰城矿区B4煤层CH4、CO2等温吸附实验规律总体表现为压力正效应、温度负效应;煤层对CO2的吸附能力明显高于对CH4的吸附能力,二者最大吸附量的比值约为2.0;9种模型中,L模型、T模型、LF模型、EL模型的拟合效果最好,最适合表征煤层对CO2的吸附特征。
Abstract:The gas adsorption characteristic of coal body determines the geological storage capacity of coal seam. Isothermal adsorption experiments were conducted in the main B4 coal seam in the mining area under different temperature conditions to study the gas adsorption characteristics of the coal seam. Nine different adsorption mathematical models were used to fit and analyze the CO2 adsorption experimental data, and the fitting effects of different models were tested. The results indicate that: the overall experimental pattern of CH4 and CO2 isothermal adsorption in B4 coal seam of Fengcheng Mining Area shows that the pressure is a positive effect and the temperature is a negative effect; the adsorption capacity of coal seams for CO2 is significantly stronger than that for CH4, and the ratio of the maximum adsorption capacity of the two is about 2.0; the L model, T model, LF model, and EL model have the best fitting effect and are most suitable for characterizing the adsorption characteristics of CO2 in coal seams.
-
地下资源开采后,破坏其原有平衡状态,经过覆岩传递到地表,是1个复杂的动态移动过程,涉及地质、采矿、空间等多方面因素,许多学者在开采沉陷动态变形原理、模型、精度等方面进行了研究,并取得了大量成果。开采沉陷动态预测模型主要有动态力学模型和时间函数模型2类,目前主要采用时间函数构建开采沉陷的动态变形过程;常用的时间函数有Knothe时间函数模型、双曲线时间函数模型、Gompertz时间函数模型、logistic 曲线时间函数模型及 Weibull 曲线时间函数模型等。
王悦汉等[1]将动态移动过程分为顶板初次垮落前、顶板垮落未充满采空区、顶板垮落已充满采空区、表土移动4个阶段,采动岩体分为连续介质、拟连续介质、非连续介质3种介质,分阶段建立了动态力学预测模型;吴侃等[2]提出了基于时间序列分析的动态预测模型;崔希民等[3]建立基于Knothe时间函数的概率积分预计方法,根据临界开采尺寸确定时间影响系数;张兵等[4-5]建立了1种新的优化分段Knothe时间函数,并提出了基于实测数据的反算时间函数对比求参法、基于概率积分参数的直接计算法的2种优化分段Knothe时间函数参数求取方法;刘玉成等[6-7]分析了下沉动态过程的Knothe时间函数模型的不足,提出了幂指数的Knothe时间函数,分析了常见地表沉陷时间函数模型的下沉、下沉速度、下沉加速度与时间关系,认为Weibull 曲线时间函数模型能完整地描述地表沉陷的动态过程;孙闯等[8]建立了考虑松散层下沉特点的双因素Knothe时间函数;王军保等[9]借鉴岩石流变力学中非定常流变模型,构建了变时间影响系数的Knothe的地表动态沉陷预测模型;许国胜等[10]回归分析得到地表动态移动变形参数与地质及开采技术参数之间的关系;李春意等[11-12]进行了基于Logistic时间函数、正态分布时间函数的动态沉陷预测;张永胜等[13]提出了以地表点最大下沉速度时刻为分界点,结合偏差改正、生长函数模型的1种分段Weibull时间函数;刘东海[14]根据常村观测资料,指出Weibull时间函数参数不固定,两参数变化规律为距边界变化类偏态分布;JAROSZ等[15]考虑开采产生的空间收敛,建立了双参数Sroka-Schober时间函数;KOWALSKI[16]在Knothe和Sroka-Schober双参数时间函数基础上,提出了广义时间函数;González Nicieza等[17]在Knothe时间函数基础上,引入正态分布函数,建立了正态分布时间函数。
综上,总体认为地表动态下沉量与时间的关系曲线呈“S”形曲线,下沉速度为0→Vmax→0过程,下沉加速度为0→+amax→-amax→0过程。
1. 地表动态移动过程覆岩变化特征
工作面上方的覆岩及地表在井下资源开采后,经历弯曲、垮落、再次压缩等过程,但在工作面不同位置的岩层完整厚度、位移、受力等经历过程不同,根据工作面回采、覆岩受力与垮落情况,将工作面回采结束后的地表分为3个部分,工作面动态分析分区图如图1。
1)竖向整体受压区(Ⅰ区)。处于工作面外侧煤柱区域,工作面开采边界至开采影响边界之间。此部分区域在整个回采期间,覆岩经历了弯曲,但没有垮落、再次压缩等过程,在竖向上是完整的。在整个开采期间竖向方向整体是受压状态。
2)中间区(Ⅱ区)。处于工作面内测靠近煤柱区域,工作面开采边界至覆岩动态完全区边界(覆岩垮落角确定边界与工作面开采边界之间)。此部分区域在整个回采期间,覆岩经历了弯曲、垮落、非完全再次压缩过程等部分或全部过程,主要是由于岩梁的存在导致垮落非充分及垮落岩石受力非充分,在竖向上部分破裂,且离工作面边界(煤柱)越远,破裂高度越大,即完整覆岩厚度减小。在整个回采期间主要经历如下几个阶段:①受到影响到工作面推到下方前,在竖直方向上整体是受到压缩变形,竖向受力增大;②推过后至覆岩垮落前,在竖直方向上整体是受到拉伸变形、竖向受力减小;③垮落后至移动稳定,在竖直方向上整体是受到压缩变形,竖向受力增大至稳定。垮落岩石部分再次受到压缩变形。
3)动态完全区(Ⅲ区)。处于工作面内测远离煤柱区域,覆岩动态完全区(覆岩垮落角确定边界至采空区中心)。此部分区域在整个回采期间,覆岩经历了弯曲、垮落、完整再次压缩过程,在竖向上部分破裂,破裂高度一致。在整个回采期间主要经历如下几个阶段:①受到影响到工作面推到下方(垮落前),在竖直方向上整体是受到压缩变形,竖向受力增大;②到再次受到压缩前,在竖直方向上整体是受到拉伸变形,竖向受力减小;③到移动稳定,在竖直方向上整体是受到压缩变形,竖向受力增大至稳定。
根据开采沉陷岩梁相关理论,地表下沉可简化为各种类型岩梁弯曲,其弯曲下沉值与刚度(弹性模量与惯性矩之积)成反比,其中惯性矩与岩梁高度的三次方成正比。结合图1,在同种情况下,上覆岩层刚度、惯性矩变化规律为:Ⅰ区(整体一致,最大)>Ⅱ区(与离开切眼距离成反比)>Ⅲ区(整体一致,最小)。
目前以时间函数为基础的各种地表动态沉陷预测模型参数主要包括:最大变形值、时间函数(时间函数不同,其参数不同),各类参数与覆岩、开采等因素有关,体现其差异,但目前时间函数都采用相同参数,即在工作面不同位置(煤柱内外)的时间函数参数都一样,这与处在工作面不同位置的受力过程、变形过程、岩层破坏情况及岩石力学性质等明显不符。
2. Weibull时间函数参数变化特点
Weibull时间函数f(t)数学表达式为:
$$ f\left(t\right)=\left\{\begin{array}{cc}km{t}^{m-1}{{\mathrm{e}}}^{-k{t}^{m}}& t\geqslant 0\\ 0& t < 0\end{array}\right. $$ (1) 式中:k、m为参数,正数;t为时间。
据此建立开采沉陷过程中地表下沉时间函数、下沉速度时间函数、下沉加速度时间函数,开采沉陷动态下沉过程曲线如图2。
下沉时间函数W(t):
$$ W\left(t\right)={W}_{0}\left(1-{{\mathrm{e}}}^{-k{t}^{m}}\right) $$ (2) 式中:W0为该点下沉的最终值。
下沉速度时间函数v(t):
$$ v\left(t\right)={W}_{0}km{t}^{m-1}{{\mathrm{e}}}^{-k{t}^{m}} $$ (3) 下沉加速度时间函数a(t):
$$ a\left(t\right)={W}_{0}km\left(m-1-km{t}^{m}\right){t}^{m-2}{{\mathrm{e}}}^{-k{t}^{m}} $$ (4) 根据理论分析及Weibull时间函数上述公式可知:① 当地面达到最大下沉速度时,下沉加速度为0,此时时间tvmax满足${t}_{vmax}=\sqrt[\leftroot{-1}\uproot{8}{{m}}]{{(m-1)}/{km}} $;② 由下沉加速度变化规律可知$m > 1 $;③ k反映了加速时间、加速度大小,k越小,加速度越小,到达最大下沉速度的时间越大;④ m越小,到达最大下沉速度的时间越小,最大下沉速度越小。
根据前面分析可知:工作面不同位置在回采期间,其覆岩移动、破坏、受力、岩梁刚度等不同,则其动态时间函数参数应不同,m、k决定了Weibull时间函数的形态变化过程,因此m、k在工作面不同位置其存在差异,以垮落区边界为原点、指向煤柱方向为正,根据前述岩梁弯曲特征简述及相关岩石力学理论,建立其与地面点距离垮落区域边界距离S的函数,工作面不同位置m、k变化情况图如图3。
在煤柱侧区域(Ⅰ区)m、k为:
$$ m={a}_1 {{\mathrm{e}}}^{-b_1\frac{s}{H}}+c_1 $$ (5) $$ k={a}_2 \frac{{{S}}_{0}}{H}+b_2 $$ (6) 在垮落区外侧至工作面边界(Ⅱ区)m、k为:
$$ {m}={a}_1 {{\mathrm{e}}}^{-b_1\frac{S}{H}}+c_1 $$ (7) $$ {k}={a}_2 \frac{{s}}{{H}}+{b}_2 $$ (8) 在垮落区内侧(Ⅲ区)m、k为:
$$ {m}={a}_1+c_1 $$ (9) $$ {k}={b}_2 $$ (10) 式中: a1、b1、c1、a2、b2为系数;S0为煤柱边界距离;H为工作面采高。
从m、k计算公式及图3中可以看出:
1)在垮落区外侧为负指数关系,以工作面边界为界,此区域分为2部分;Ⅰ区(煤柱范围)m变化较小,基本为常数;Ⅱ区(中间区)m变化较大,在垮落区内侧,m为一常数。
2)以工作面边界、垮落区边界为界,此区域分为3部分;Ⅰ区(煤柱范围)k为常数;Ⅱ区(中间区)m为线性变化;在垮落区内侧,k为一常数。
由于不同位置上覆岩层经历的受力、破坏、变形过程不同,根据开采沉陷相关研究,可采用如下方法确定垮落范围d。① 方法1:采用覆岩力学理论,根据岩梁的相关破裂理论确定垮落边界;② 方法2:采用概率积分法拐点偏移距;③ 方法3:采用顶板覆岩垮落角φ确定,$d={H}/{\mathrm{tan}\;\;\varphi } $。
3. 案例分析
某工作面采用长壁一次性采全高采煤法回采。走向为774 m,倾斜宽为220.7 m。工作面开采深度为470~515 m。根据走向观测数据,采用上面的计算模型,获得的反演模型系数见表1。
表 1 反演模型系数Table 1. Coefficients of inversion model点号 k m R2 Z17 0.271 3 0.294 3 0.469 7 Z18 0.163 1 0.386 6 0.629 3 Z19 0.621 7 0.430 6 0.882 6 Z20 0.131 8 0.401 2 0.945 1 Z21 0.293 1 0.322 1 0.899 5 Z22 0.271 7 0.365 3 0.880 4 Z23 0.135 8 0.530 0 0.925 4 Z24 0.009 9 1.133 0 0.924 2 Z25 0.002 2 1.494 0 0.928 5 Z26 0.001 4 1.545 0 0.952 1 Z27 0.001 7 1.449 0 0.941 8 Z28 0.002 5 1.301 0 0.927 0 Z29 0.002 7 1.261 0 0.939 5 Z30 0.001 0 1.414 0 0.941 2 根据工作面开采情况,以工作面内侧90 m为垮落边界,以此确定各点距离边界的S值,根据表1反演动态时间函数$\varphi \left(t\right) $如下:
$$ \varphi \left(t\right)=1-{{\mathrm{e}}}^{-k{t}^{m}} $$ (11) 在采空区垮落内侧m、k为:
$$ m=1.419\;4, k=0.001\;908 $$ (12) 在采空区垮落外侧m、k为:
$$ m\left(S\right)=1.099\cdot {{\mathrm{e}}}^{-12.97\times \frac{S}{495}}+0.320\;4 \;\;\;(R^{2}=0.921\;6) $$ (13) $$ k=1.517\times \frac{{S}}{495}+0.001\;908 \;\;\;(R^{2}=0.919\text{,}S<0.2H) $$ $$ k=0.305 \qquad\qquad (S>0.2H)$$ (14) 根据工作面开采情况及参数确定情况,分别以开采结束后实测结果作为最终值、概率积分法预测值作为最终值,进行开采第523 d的动态预测,实测结果、概率积分法预测值作为最终值的动态预测与实测对比如图4(图中距离以工作面边界为原点,采空区方向为正)。
以实测结果作为最终值:动态预测中误差为110 mm,为最大下沉值的3.5%;以概率积分法预测结果作为最终值:动态预测中误差为130 mm,为最大下沉值的4.2%。
4. 结 语
1)分析了工作面上方不同位置的覆岩在开采过程中在竖向上的受力、变形、破坏等变化过程,根据煤柱、垮落情况将覆岩与地表划分为竖向整体受压区、中间区、动态完全区3个区域。
2)以Weibull时间函数动态预测模型为例,分析了m、k 2个参数的空间分布规律及特地,构建了2个参数的分区数学模型及选取准则。
3)分别以实测结果、概率积分法预测结果作为工作面开采结束后的最终值,依据本文动态预测模型及参数确定方法,动态预测中误差分别为110、130 mm,为最大下沉值的3.5%、4.2%,证明本模型的可靠性。
-
表 1 丰城矿区B4煤层煤质特征
Table 1 Coal quality characteristics of B4 coal seam in Fengcheng Mining Area
井田 水分/% 灰分/% 挥发分/% 硫分/% 煤类 曲江 1.39 22.69 21.82 2.09 焦煤 尚庄 1.45 25.35 17.21 1.36 瘦煤 表 2 不同温度条件下曲江煤样(QJG01)吸附模型的拟合度
Table 2 Fitting degrees of adsorption models for Qujiang coal samples (QJG01) under different temperature conditions
温度/℃ 拟合度R2 L模型 F模型 T模型 LF模型 EL模型 DR模型 二次模型 三次模型 对数模型 20 0.998 0.996 0.999 0.999 0.999 0.952 0.974 0.994 0.998 30 0.999 0.996 1.000 1.000 1.000 0.980 0.983 0.996 0.998 40 0.997 0.998 0.999 0.998 0.999 0.991 0.984 0.993 0.994 50 0.996 0.999 0.999 0.999 0.999 0.993 0.984 0.993 0.991 表 3 9种模型拟合平均相对误差和平均标准方差
Table 3 Mean relative error and mean standard deviation of fitting for nine models
模型参数 L模型 F模型 T模型 LF模型 EL模型 DR模型 二次模型 三次模型 对数模型 平均相对误差 0.024 0.017 0.012 0.013 0.012 0.047 0.046 0.028 0.036 平均标准方差 0.559 0.641 0.330 0.336 0.331 1.751 1.265 0.712 0.728 -
[1] 袁志刚. 碳达峰碳中和国家战略行动路线图[M]. 北京:中国经济出版社,2021. [2] 撒占友,吴静波,陆卫东,等. 注二氧化碳强化煤层气开采研究进展[J]. 煤矿安全,2022,53(10):32−43. SA Zhanyou, WU Jingbo, LU Weidong, et al. Research progress of enhanced coalbed methane exploitation by carbon dioxide injection[J]. Safety in Coal Mines, 2022, 53(10): 32−43.
[3] 桑树勋,刘世奇,陆诗建,等. 工程化CCUS全流程技术及其进展[J]. 油气藏评价与开发,2022,12(5):711−725. SANG Shuxun, LIU Shiqi, LU Shijian, et al. Engineered full flowsheet technology of CCUS and its research progress[J]. Petroleum Reservoir Evaluation and Development, 2022, 12(5): 711−725.
[4] 米剑锋,马晓芳. 中国CCUS技术发展趋势分析[J]. 中国电机工程学报,2019,39(9):2537−2543. MI Jianfeng, MA Xiaofang. Development trend analysis of carbon capture, utilization and storage technology in China[J]. Proceedings of the CSEE, 2019, 39(9): 2537−2543.
[5] 张洪涛,文冬光,李义连,等. 中国CO2地质埋存条件分析及有关建议[J]. 地质通报,2005,24(12):1107−1110. doi: 10.3969/j.issn.1671-2552.2005.12.004 ZHANG Hongtao, WEN Dongguang, LI Yilian, et al. Conditions for CO2 geological sequestration in China and some suggestions[J]. Geological Bulletin of China, 2005, 24(12): 1107−1110. doi: 10.3969/j.issn.1671-2552.2005.12.004
[6] 刘廷,马鑫,刁玉杰,等. 国内外CO2地质封存潜力评价方法研究现状[J]. 中国地质调查,2021,8(4):101−108. LIU Ting, MA Xin, DIAO Yujie, et al. Research status of CO2 geological storage potential evaluation methods at home and abroad[J]. Geological Survey of China, 2021, 8(4): 101−108.
[7] 李甫成,张杨,张晓娟,等. 深部咸水层CO2地质储存适宜性评价方法研究[J]. 冰川冻土,2014,36(3):649−660. LI Fucheng, ZHANG Yang, ZHANG Xiaojuan, et al. Suitability evaluation method of CO2 geological sequestration in deep saline aquifers[J]. Journal of Glaciology and Geocryology, 2014, 36(3): 649−660.
[8] SHUKLA R, RANJITH P, HAQUE A, et al. A review of studies on CO2 sequestration and caprock integrity[J]. Fuel, 2010, 89(10): 2651−2664. doi: 10.1016/j.fuel.2010.05.012
[9] 桑树勋,王冉,周效志,等. 论煤地质学与碳中和[J]. 煤田地质与勘探,2021,49(1):1−11. doi: 10.3969/j.issn.1001-1986.2021.01.001 SANG Shuxun, WANG Ran, ZHOU Xiaozhi, et al. Review on carbon neutralization associated with coal geology[J]. Coal Geology & Exploration, 2021, 49(1): 1−11. doi: 10.3969/j.issn.1001-1986.2021.01.001
[10] 吴建光,叶建平,唐书恒. 注入CO2提高煤层气产能的可行性研究[J]. 高校地质学报,2004,10(3):463−367. doi: 10.3969/j.issn.1006-7493.2004.03.016 WU Jianguang, YE Jianping, TANG Shuheng. A feasibility study on CO2 injection for enhancing the coalbed methane recovery[J]. Geological Journal of China Universities, 2004, 10(3): 463−367. doi: 10.3969/j.issn.1006-7493.2004.03.016
[11] 韩思杰. 深部无烟煤储层CO2-ECBM的CO2封存机制与存储潜力评价方法[D]. 徐州,中国矿业大学,2020. [12] 代世峰,张贝贝,彭苏萍,等. 河北开滦矿区晚古生代煤对CO2和CH4气体吸附模型探讨[J]. 地质学报,2009,83(5):731−737. DAI Shifeng, ZHANG Beibei, PENG Suping, et al. Models of pure CO2 and pure CH4 adsorption on the late paleozoic coals from the Kailuan Coalfield, Hebei, China[J]. Acta Geologica Sinica, 2009, 83(5): 731−737.
[13] 李向东,冯启言,刘波,等. 注入二氧化碳驱替煤层甲烷的试验研究[J]. 洁净煤技术,2009,16(2):101−103 LI Xiangdong , FENG Qiyan , LIU Bo, et al. Study on displacement coalbed methane by carbon dioxide injection[J]. Clean Coal Technology, 2009, 16(2): 101−103.
[14] 吴迪,孙可明. 不同温度条件下型煤吸附CH4/CO2混合气的实验研究[J]. 岩石力学与工程学报,2013,32(S2):3291−3296. WU Di, SUN Keming. Experimental study for coal briquette adsorption of CH4/CO2 mixture under different temperatures[J]. Chinese Journal of Rock Mechanics and Engineering, 2013, 32(S2): 3291−3296.
[15] 王晋,王延斌,范晶晶,等. 注CO2置换煤层CH4试验研究[J]. 科技导报,2015,33(17):84−89. WANG Jin, WANG Yanbin, FAN Jingjing, et al. Experimental study on coal-bed methane displacement by carbon dioxide injection[J]. Science & Technology Review, 2015, 33(17): 84−89.
[16] 王瑞雪. 不同变质程度煤对CH4和CO2吸附特征差异的实验研究[D]. 徐州:中国矿业大学,2020. [17] 苏现波,张丽萍,林晓英. 煤阶对煤的吸附能力的影响[J]. 天然气工业,2005,25(1):19−21. doi: 10.3321/j.issn:1000-0976.2005.01.006 SU Xianbo, ZHANG Liping, LIN Xiaoying. Influence of coal rank on coal adsorption capacity[J]. Natural Gas Industry, 2005, 25(1): 19−21. doi: 10.3321/j.issn:1000-0976.2005.01.006
[18] 降文萍,崔永君. 深部煤层封存CO2的地质主控因素探讨[J]. 中国煤炭地质,2010,22(11):1−6. doi: 10.3969/j.issn.1674-1803.2010.11.01 JIANG Wenping, CUI Yongjun. A discussion on main geologic controlling factors of CO2 sequestration in deep coal seams[J]. Coal Geology of China, 2010, 22(11): 1−6. doi: 10.3969/j.issn.1674-1803.2010.11.01
[19] 于洪观. 煤对CH4、CO2、N2及其二元混合气体的吸附特性、预测和CO2驱替CH4的研究[D]. 青岛:山东科技大学,2005. [20] 朱苏阳. 煤层气的吸附—解吸机理及应用研究[D]. 成都:西南石油大学,2018. [21] 刘成龙. 煤层CO2埋藏主控地质因素及对煤储层影响研究[D]. 北京:中国矿业大学(北京),2017. [22] 郝光旭. 煤对CH4、CO2和N2混合气体等温吸附实验研究[D]. 北京:中国地质大学(北京),2021. [23] 于洪观,范维唐,孙茂远,等. 煤中甲烷等温吸附模型的研究[J]. 煤炭学报,2004,29(4):463−467. doi: 10.3321/j.issn:0253-9993.2004.04.019 YU Hongguan, FAN Weitang, SUN Maoyuan, et al. Study on fitting models for methane isotherms adsorption of coals[J]. Journal of China Coal Society, 2004, 29(4): 463−467. doi: 10.3321/j.issn:0253-9993.2004.04.019
[24] 李延河,李喜员,孙矩正,等. 煤对超临界CO2的吸附实验及不同表征模型对比[J]. 煤矿安全,2023,54(6):18−26. LI Yanhe, LI Xiyuan, SUN Juzheng, et al. Adsorption experiments of supercritical CO2 by coals and comparisons of different description models[J]. Safety in Coal Mines, 2023, 54(6): 18−26.
-
期刊类型引用(5)
1. 李鹏. 基于FLAC3D模拟的14103辅运顺槽支护参数设计. 江西煤炭科技. 2025(01): 90-92+96 . 百度学术
2. 马清水,郭瑞,韩伟,宋永明,梁燕翔,刘耀,王佳明. 基于长短期记忆网络的煤巷支护设计参数预测. 能源与环保. 2025(02): 235-240 . 百度学术
3. 李迎富,宋武,杨子仪,赵江梅,曾凡刚. 掘进巷道锚杆支护优化研究. 煤矿机械. 2024(05): 116-118 . 百度学术
4. 陈万辉,郭瑞,韩伟,宋永明,梁燕翔,刘耀,王佳明,许娜,孟波. 煤矿巷道支护方案智能设计研究. 工矿自动化. 2024(08): 76-83+90 . 百度学术
5. 彭启友,缪易辰,潘文,尚江淮. 吸能锚杆数值模拟研究现状与展望. 建筑技术. 2024(20): 2546-2551 . 百度学术
其他类型引用(1)