高级检索

低渗煤层CO2-ECBM过程的CO2和CH4的扩散特征

刘峻麟, 刘会虎, 张琨, 薛生

刘峻麟,刘会虎,张 琨,等. 低渗煤层CO2-ECBM过程的CO2和CH4的扩散特征[J]. 煤炭科学技术,2023,51(S1):112−121

. DOI: 10.13199/j.cnki.cst.2022-1296
引用本文:

刘峻麟,刘会虎,张 琨,等. 低渗煤层CO2-ECBM过程的CO2和CH4的扩散特征[J]. 煤炭科学技术,2023,51(S1):112−121

. DOI: 10.13199/j.cnki.cst.2022-1296

LIU Junlin,LIU Huihu,ZHANG Kun,et al. Diffusion characteristics of CO2 and CH4 in CO2-ECBM process of low permeability coal seam[J]. Coal Science and Technology,2023,51(S1):112−121

. DOI: 10.13199/j.cnki.cst.2022-1296
Citation:

LIU Junlin,LIU Huihu,ZHANG Kun,et al. Diffusion characteristics of CO2 and CH4 in CO2-ECBM process of low permeability coal seam[J]. Coal Science and Technology,2023,51(S1):112−121

. DOI: 10.13199/j.cnki.cst.2022-1296

低渗煤层CO2-ECBM过程的CO2和CH4的扩散特征

基金项目: 

国家自然科学基金资助项目(42277483);能源研究院合肥综合性国家科学中心资助项目(21KZS218);安徽省自然科学基金资助项目(2008085MD121)

详细信息
    作者简介:

    刘峻麟: (1998—),男,重庆璧山人,硕士研究生。E-mail:603093502@qq.com

    通讯作者:

    刘会虎: (1978—),男,江西九江人,博士生导师,教授。E-mail:xixiinformation@163.com

  • 中图分类号: TD712

Diffusion characteristics of CO2 and CH4 in CO2-ECBM process of low permeability coal seam

Funds: 

National Natural Science Foundation of China (42277483); Energy Research Institute Hefei Comprehensive National Science Center Project (21KZS218); Natural Science Foundation of Anhui Province (2008085MD121)

  • 摘要:

    气体扩散能力的强弱是决定CO2-ECBM成败的关键因素之一。以中国两淮低渗含煤区刘庄矿13煤和祁东矿7煤的煤储层孔隙特征为研究对象,开展CO2-ECBM过程中CO2和CH4气体的扩散行为研究,以期获取CO2和CH4气体在储层条件下的扩散规律,为研究低渗煤层CO2-ECBM过程的扩散规律提供理论依据。研究结果表明:刘庄矿13煤和祁东矿7煤2种煤样中气体扩散类型以Fick型扩散为主,占比50%以上,其次为过渡型扩散,占比接近30%,CH4的Knudsen型扩散约占10%,CO2的Knudsen型扩散接近20%;CO2和CH4气体进行Fick型扩散和过渡型扩散时,其扩散系数与温度变化关系不明显,在压力2 MPa,刘庄矿13煤中的气体扩散系数约为10 m2/s,祁东矿7煤中的气体扩散系数达到13 m2/s以上,且随压力增加2种类型的扩散系数均开始降低,在0~5 MPa的压力区间内下降速率最快,当压力值达到20 MPa后,气体2种扩散类型的扩散系数基本趋于稳定,同2 MPa时的扩散性系数相比,缩小了近10倍;当气体扩散为Knudsen型扩散时,温度每增加6 ℃,其扩散系数增加0.2 m2/s,呈较缓的线性增长;对不同扩散类型的气体扩散系数进行比较,可知CH4比CO2在Fick型扩散中更具优势,CO2则在另外2种扩散类型中更具优势;通过比较CH4在不同扩散类型中的扩散速率,发现CH4气体在低渗煤中以Fick型或Knudsen型扩散时扩散较快,以过渡型扩散为扩散类型为主时扩散速度较慢;CO2-ECBM过程中,CO2在微孔及小孔中相比CH4扩散更快,能够更加有效地置换出微孔中的CH4,提高CH4的产量。

    Abstract:

    The ability of gas diffusion is one of the key factors that determine the success or failure of CO2-ECBM. Taking the pore characteristics of coal reservoirs of No. 13 Coal Seam of Liuzhuang Coal Mine and No. 7 Coal Seam of Qidong Coal Mine in the low permeability coal bearing area of Huaihe River and Huaihe River in China as the research object, the diffusion behavior of CO2 and CH4 gas in the process of CO2-ECBM was studied in order to obtain the diffusion law of CO2 and CH4 gas under the reservoir conditions and provide a theoretical basis for the study of the diffusion law of the low permeability coal seam in the process of CO2-ECBM. The research results show that the Fick type diffusion is the main type of gas diffusion, accounting for more than 50%, followed by the transitional type diffusion, accounting for nearly 30%, the Knudsen type diffusion of CH4 accounts for about 10%, and the Knudsen type diffusion of CO2 accounts for nearly 20%; When CO2 and CH4 gas undergo Fick type diffusion and transitional type diffusion, their diffusion coefficients have no obvious relationship with temperature change. At a pressure of 2 MPa, the gas diffusion coefficient in Coal 13 of Liuzhuang Mine is about 10 m2/s, and the gas diffusion coefficient in Coal 7 of Qidong Mine is above 13 m2/s. With the increase of pressure, the diffusion coefficients of both types begin to decrease, and the drop rate is the fastest in the pressure range of 0~5 Mpa. When the pressure reaches 20 MPa, The diffusion coefficients of the two types of gas diffusion basically tend to be stable, which is nearly 10 times smaller than that at 2 MPa; When the gas diffusion is Knudsen type, the diffusion coefficient increases by 0.2 m2/s for every 6 ℃ increase in temperature, showing a slow linear growth; Comparing the gas diffusion coefficients of different diffusion types, CH4 has more advantages than CO2 in Fick type diffusion, and CO2 has more advantages in the other two diffusion types; By comparing the diffusion rate of CH4 in different diffusion types, the diffusion rate of CH4 gas in low permeability coal is faster when Fick type or Knudsen type is used as diffusion type, and slower when transition type is used as diffusion type; In the process of CO2-ECBM, CO2 diffuses faster in micropores and small pores than CH4, which can more effectively replace CH4 in micropores and improve CH4 production.

  • 面对能源供需格局新变化、国际能源发展新趋势,CO2-ECBM作为集温室气体减排和新能源开发于一体的煤层气强化开采手段,符合国家的能源安全新战略要求。在CO2-ECBM技术中,CO2注入煤层后会在气体浓度梯度的驱动下扩散流动,气体扩散能力的强弱也是决定CO2-ECBM成败的关键因素,了解气体在煤层中的扩散行为有助于理解CO2-ECBM过程中CO2是如何进行贮存和CH4如何进行采收[1]

    煤的扩散系数是指当浓度梯度为一个单位时,单位时间内通过单位面积煤中气体量;其大小主要取决于气体成分、煤岩煤质及其温度–压力–应力等条件[2-9]。张时音[3]和张登峰[4]等应用吸附扩散理论模型计算了不同煤阶煤的扩散系数,认为扩散系数与总孔容正相关,扩散系数与吸附量、比表面积成正比。CHARRERE等[10]在283~333 K范围内研究了温度对扩散的影响,结果表明随着温度的升高,CO2和CH4吸附速率均增大。同样地,BUSCH等[11]在2011年的综述文章中也得到了相同的结论。STAIB等[12]采用湿煤研究了不同气体在煤中的扩散,结果发现当温度由300 K升温至313 K,CH4和CO2吸附速率增大。李志强等[13]研究了升温和恒温情况下瓦斯在湿煤粒中的扩散,发现同压不同温(303~373 K、4 MPa)气体吸附后的恒温扩散系数随温度的升高呈现先增大后减小的趋势,在333 K时扩散系数达到极大值,而同温同压(293~333 K、0.5 MPa)气体吸附后升温扩散时的扩散系数随温度的升高呈指数关系增大。此外,相同条件下CO2在煤中的扩散系数明显大于CH4

    在压力对于气体在煤层中的扩散影响方面,CIEMBRONIEWICZ等[14]使用简化的unipore模型研究了CO2在硬煤中的动力学特性,结果发现随着压力的增加扩散系数增大。CHARRIERE等[10]使用简化的unipore模型研究了压力对扩散系数的影响,结果发现与0.1 MPa相比,相同温度下压力为5 MPa时CO2和CH4在煤中的扩散系数更大。

    国内外学者对煤的渗透性及其控制因素方面已开展了较多的研究[15-16],对于煤中气体扩散行为的研究相对较少[17],在以往的煤层气体扩散研究中[18-21],压力条件较局限,压力范围较小,压力值较低,未能更加实际地反映储层条件下气体的扩散特征,且煤层中不同孔隙条件下的气体扩散特征没能直观的体现。

    煤层气开采效益和产量在近年来的生产中没有达到预期效果,刘庄煤矿与祁东煤矿作为两淮地区典型的低渗煤层气田,以该地区煤样为研究对象,能够为研究两淮地区低渗煤层气田的高效开发利用提供理论依据。笔者研究了储层条件(24~48 ℃、2~20 MPa)下刘庄矿13煤和祁东矿7煤2个地区煤样在不同孔隙条件下的扩散模式,模拟不同位置储层压力条件下气体在煤样中的扩散方式,以揭露CO2和CH4气体在低渗煤样中的扩散规律为目的,进一步探究2种气体在煤样中扩散的差异性,以期为研究低渗煤层CO2-ECBM过程的扩散规律提供理论依据。

    刘庄煤矿和祁东煤矿为两淮地区典型的低渗煤层气田,其中刘庄煤矿处于淮南复向斜中次一级褶皱(陈桥背斜)的南翼,煤层主要分布在谢桥向斜北翼及西转折部位,本区分别以板集断裂、陈桥断裂为东、西部边界,按构造形态及断层分布情况(图1[22],煤矿中二叠系的山西组与上下石盒子组为主要含煤层段,刘庄13煤层为主要可采煤层之一,煤种类以气煤为主,平均总厚度18.51 m。祁东煤矿位于淮北煤田宿县矿区宿南向斜内,宿南向斜的大地构造位置属徐淮隆起的徐宿坳陷区的南部(图2[23],祁东矿主要含煤地层为二叠系山西组、下石盒子组、上石盒子组,其中祁东7煤层位于中~下部富煤段,煤种类为1/3焦煤和肥煤。

    图  1  刘庄煤矿构造纲要图[22]
    Figure  1.  Structural outline map of Liuzhuang Coal Mine[22]
    图  2  祁东煤矿构造纲要图[23]
    Figure  2.  Structural outline map of Qidong Coal Mine[23]

    试验对象为刘庄13煤和祁东7煤的煤岩样品,煤样工业分析结果见表1

    表  1  煤样工业分析表
    Table  1.  Analysis table of coal sample industry
    样品名刘庄13煤祁东7煤
    煤类气煤肥煤
    采样深度/m7631 044
    水分Mad/%1.250.88
    挥发分Vdaf/%39.3438.59
    固定碳FCad/%39.8149.51
    灰分Aad/%19.6011.02
    平均真相密度/(g·cm-31.361.45
    镜质组反射率RO/%0.72520.7558
    下载: 导出CSV 
    | 显示表格

    气体在煤层中的扩散,其本质是气体分子的不规则热运动。煤作为一种典型的多孔介质,可将气体在煤层中的扩散总体分为:气相扩散、吸附相扩散、溶解相扩散、固溶体扩散4种类型。其中气相扩散分为:Fick型扩散、Knudsen型扩散、过渡型扩散3种扩散类型。气体扩散类型可以用表示孔隙直径和分子运动平均自由程相对大小的诺森(Knudsen)数表示[20]

    $$ {K_{\text{n}}} = \frac{d}{\lambda } $$ (1)

    式中,d为孔隙平均直径,nm;λ为气体分子平均自由程,nm。

    理想气体下的分子平均自由程可以通过式(2)[19]计算:

    $$ \lambda = \frac{{{k_B}{T}}}{{\sqrt 2 \pi {d_0}^2p}} $$ (2)

    式中,kB为玻尔兹曼常数,1.38×10−23,J/K;T为绝对温度,K;p为气体压力,MPa;d0为分子有效直径,nm。

    根据诺森数的取值,可以划分气体的扩散类型(表2),分子自由程的值可以判断气体分子在煤层不同孔隙直径中的扩散类型,该类型也决定了气体在煤层中的扩散速率,根据CO2和CH4的分子直径,模拟不同埋深的温度压力条件,可计算得到气体分子的平均自由程(表3),根据表3所示进一步计算能够得到关于不同埋深条件下的不同扩散类型所对应的孔隙直径(表4)。

    表  2  扩散类型及其对应诺森数
    Table  2.  Diffusion type and its counterpart
    扩散类型Kn
    Fick型扩散Kn≥10
    过渡型扩散0.1≤Kn≤10
    Knudsen型扩散Kn≤0.1
    下载: 导出CSV 
    | 显示表格
    表  3  模拟不同埋深条件下的气体分子自由程
    Table  3.  Simate the molecular free path of gas under different burial depths
    模拟埋深

    条件/m
    温度T/K压力p/MPa分子有效直径d0/nm平均分子自由程λ/nm
    CH4CO2CH4CO2
    −1 500321.15150.3800.3300.4610.611
    −1 200315.15120.5650.749
    −900309.1590.7390.980
    −600303.1561.0871.441
    −300297.1561.0651.413
    下载: 导出CSV 
    | 显示表格
    表  4  不同扩散类型的煤体孔分布情况
    Table  4.  Distribution of coal body holes in different diffusion types
    埋深条件
    /m
    扩散孔径分布/nm
    Knudsen型扩散过渡型扩散Fick型扩散
    CH4CO2CH4CO2CH4CO2
    −1 5000.0460.0610.046~4.610.061~6.114.616.11
    −1 2000.0570.0750.057~5.70.075~7.455.657.49
    −9000.0740.0980.074~7.390.098~9.807.399.80
    −6000.1080.1440.108~10.870.144~14.4210.8714.42
    −3000.1070.1410.107~10.650.141~14.1310.6514.13
    下载: 导出CSV 
    | 显示表格

    表3表4中数据显示,埋藏深度为300 m时,CO2和CH4的Fick型扩散类型分别发生在直径14.13 nm以及10.65 nm以上的孔隙中;埋藏深度为1 500 m时,CO2和CH4的Fick型扩散类型分别发生在直径4.61 nm以及6.11 nm以上的孔隙中。该结果表明随温度及压力的增加即埋藏深度的增加,Fick型扩散类型会发生在更小的孔隙直径中,代表气体分子随着压力及温度的增加,同煤层孔隙壁所产生的碰撞会减小,这将更有利于气体的在煤层中的扩散。同时,根据表中数据可看出,随埋深增加,CH4的Knudsen型扩散所处孔隙直径分布由原来的小于0.107 nm缩小到了小于0.046 nm;过渡型扩散所处孔隙直径分布分别由原来0.107~10.6 nm缩小到了0.046~4.61 nm,产生该种扩散类型的孔隙直径区间进一步缩小。

    以刘庄矿13煤和祁东矿7煤两处的煤样作为试验样品,取温度条件24、30、36、42 、48 ℃,取压力条件2、4、6、8、10、12、14、16、18、20 MPa,该温压条件能模拟2个矿区各深度煤层的储层情况,通过计算后获取两处煤矿的煤样在不同孔隙结构中气体的扩散类型分布情况以及不同扩散类型的扩散系数。

    根据式(2)计算可得物理模拟试验条件下(温度24~48 ℃,压力2~20 MPa))对应的分子自由程变化范围(表5)。

    表  5  气体分子平均自由程结果表
    Table  5.  Table of the average free course results of the gas molecules
    气体种类λ/nm
    CO20.424~4.580
    CH40.295~3.454
    下载: 导出CSV 
    | 显示表格

    根据储层条件(24~48 ℃、2~20 MPa),计算分子自由程变化范围为CO2(0.424~4.580 nm)、CH4(0.295 ~ 3.454 nm),根据诺森数Kn的值将分子的扩散类型分为:Fick型扩散、Knudsen型扩散、过渡型扩散,由分类标准(表2)可知:当Kn≤0.1,即5 nm以下的孔隙中,CO2在煤岩中以Knudsen型扩散为主;当0.1≤Kn≤10,即5 nm以上40 nm以下的孔隙中,CO2在煤岩中以过渡型扩散为主;当Kn≥10,即在40 nm以上的孔隙中,CO2在煤岩中以Fick型扩散为主。同理,当Kn≤0.1,孔隙大小在3nm以下时,CH4在煤岩中以Knudsen型扩散为主;当0.1≤Kn≤10,孔隙条件在3~30 nm时,CH4在煤岩中以过渡型扩散为主;当Kn≥10,孔隙条件在30 nm以上时,CH4在煤岩中以Fick型扩散为主。

    根据本次研究煤岩样品的压汞试验得到的孔隙结构分布,通过式(2)计算所得诺森数Kn,对不同孔径孔隙进行占比计算,具体结果见表6

    表  6  不同孔径中CO2和CH4气体的不同扩散类型占比
    Table  6.  Results of different diffusion types of CO2 and CH4 gas in different pore sizes
    气体类型煤岩样品总孔容/(cm3·g−1孔隙体积分布/%以诺森数分类计/%
    微孔
    (<10 nm)
    小孔
    (10~102 nm)
    中孔
    (102~103 nm)
    大孔
    (103~104 nm)
    Fick型扩散过渡型
    扩散
    Knudsen型
    扩散
    CO2刘庄0.01534.78038.35310.28016.58852.43630.17417.390
    祁东0.01936.45737.96610.93414.64450.88830.88418.228
    CH4刘庄0.01534.78038.35310.28016.58856.69732.86910.434
    祁东0.01936.45737.96610.93414.64455.10633.95710.937
    下载: 导出CSV 
    | 显示表格

    煤岩样品的孔隙分布及扩散类型占比,可以分析得出在该煤岩样品中气体的主要扩散类型,考虑气体在煤层中的运移时,便可采用该类型的扩散渗流理论来分析气体的流动规律。在本次的煤岩样品刘庄矿13煤和祁东矿7煤中,微孔及小孔孔容占比最高,分别约占总孔容的35%和38%。同时,储层条件下(24~48 ℃、2~20 MPa),样品中的CO2和CH4气体主要以Fick型扩散类型为主,约占总扩散类型的50%以上,过渡型扩散约占30%,CO2的Knudsen型扩散约占18%,CH4的Knudsen型扩散最少,仅占总体的约10%。

    不同的扩散类型有不同的扩散系数计算方式,Fick型扩散系数(Df)主要取决于温度及压力,同煤岩孔径大小无关,Knudsen型扩散系数(Dk)主要取决于温度及压力,同煤岩孔径大小无关;过渡型扩散系数(Dp)处于Fick型扩散以及Knudsen型扩散两者之间,同时受到温度、压力以及孔隙大小的共同影响。

    1)Fick型扩散系数。Fick型扩散系数Df可表示为[24]

    $$ {D_{\rm{f}}} = \frac{1}{3}\sqrt {\frac{{k_{\rm{B}}^3}}{{M{\pi ^3}}}} \frac{{{T^{\frac{3}{2}}}}}{{pd_0^2}} $$ (3)

    有效Fick扩散系数由式(4)计算[20]

    $$ {D_{{\text{fe}}}} = \frac{{{D_{\rm{f}}}\theta }}{\tau } $$ (4)

    式中,$ M $为气体分子摩尔质量,g/mol;$ \theta $为多孔介质孔隙率,%;$ \tau $为曲折因子,为修缮扩散路径变化而引入。

    2)Knudsen型扩散系数。Knudsen型扩散系数由式(5)计算[24]

    $$ {D_{\rm{k}}} = \frac{1}{3}d\sqrt {\frac{{8{k_{\rm{B}}}T}}{{\pi M}}} $$ (5)

    有效Knudsen型扩散系数由式(6)计算[20]

    $$ {D_{{\text{ke}}}} = \frac{{{D_{\rm{k}}}\theta }}{\tau } $$ (6)

    3)过渡型扩散系数。过渡型扩散受Fick型扩散和Knudsen型扩散2种扩散类型制约,恒压下的有效过渡型扩散系数同Fick型扩散和Knudsen扩散2种类型有关,其有效扩散系数由式(7)计算[20]

    $$ {D_{{\text{pe}}}} = {\left(\frac{1}{{{D_{{\rm{fe}}}}}} - \frac{1}{{{D_{{\rm{ke}}}}}}\right)^{ - 1}} $$ (7)

    根据上述公式,分别计算刘庄矿13煤和祁东矿7煤样品在不同温度压力(24~48 ℃、2~20 MPa)下的Fick型扩散系数、Knudsen型扩散系数以及过渡型扩散系数,具体计算结果如图37所示,显示了在不同压力温度条件下Fick扩散系数、Knudsen型扩散系数以及过渡型扩散系数的变化趋势。

    图  3  研究区样品中CO2的Fick型扩散系数随压力变化
    Figure  3.  Change diagram of the Fick-type diffusion coefficient of CO2 in the study area sample along with the pressure
    图  4  研究区样品中CH4的Fick型扩散系数随压力变化
    Figure  4.  Change diagram of the Fick-type diffusion coefficient of CH4 in the study area sample along with the pressure
    图  5  研究区样品中气体的Knudsen型扩散系数随温度变化
    Figure  5.  Change diagram of Knudsen-type diffusion coefficient of gas in the study area
    图  6  研究区样品中CO2的过渡型扩散系数随压力变化
    Figure  6.  Change diagram of the transition-type diffusion coefficient of CO2 along with the pressure in the study area samples

    图3图4可以看出,当温度条件48 ℃,压力条件20 MPa时,刘庄矿13煤样品中CO2的Fick型扩散系数为0.82 m2/s,CH4的Fick型扩散系数为0.97 m2/s;祁东矿7煤样品中CO2的Fick型扩散系数为1.36 m2/s,CH4的Fick型扩散系数为1.62 m2/s。温度条件48 ℃,压力条件2 MPa时,刘庄矿13煤样品中CO2的Fick型扩散系数为8.21 m2/s,CH4的Fick型扩散系数为9.74 m2/s;祁东矿7煤样品中CO2的Fick型扩散系数为13.61 m2/s,CH4的Fick型扩散系数为16.16 m2/s。

    同时可以看出,同温度条件下,随着压力的增加,气体的Fick型扩散系数逐渐减小,且初始时扩散系数随压力变化更明显,当压力条件增加至20 MPa后,达到一个平衡值,基本不再随压力增加而增加;同压力条件下,扩散系数的大小随温度的变化并不明显,并且当压力增加时这种影响会变得更小,但温度的增加仍与扩散系数呈正相关,如图3所示。由此分析可知,储层中的气体的扩散系数主要取决与储层中压力,温度的影响性较小。同时,在同一样品中温压条件相同的条件下,CO2的扩散系数明显低于CH4的扩散系数,可见气体分子量也是影响Fick型扩散系数的重要条件,气体分子量小的CH4同煤岩孔壁的接触更少,Fick扩散更加占优势。

    图5展示了研究区样品中气体的Knudsen型扩散系数与温度的关系。

    图5中可以看出,CO2和CH4气体的Knudsen型扩散系数随温度增加每6 ℃,扩散系数增加0.2 m2/s,受温度影响变化不明显,但其受煤层孔隙直径及气体类型影响明显。

    取刘庄矿13煤煤样孔隙直径为5 nm,祁东矿7煤煤样孔隙直径为3 nm,由于Knudsen型扩散系数同压力无关,因此取不考虑压力条件,当温度条件为48 ℃时,刘庄矿13煤样品中CO2的Knudsen型扩散系数为2.35 m2/s,CH4的Knudsen型扩散系数为1.41 m2/s,祁东矿7煤样品中CO2的Knudsen型扩散系数为2.79 m2/s,CH4的Knudsen型扩散系数为1.68 m2/s;温度条件24 ℃时,刘庄矿13煤样品中CO2的Knudsen型扩散系数为2.26 m2/s,CH4的Knudsen型扩散系数为1.36 m2/s,祁东矿7煤样品中CO2的Knudsen型扩散系数为2.69 m2/s,CH4的Knudsen型扩散系数为1.61 m2/s。

    根据数据比较分析,同一煤岩样品中气体的Knudsen型扩散系数变化会随温度增加有较小的增长。同温度条件下,气体的类型对于Knudsen型扩散系数有着更大的影响,CO2的扩散系数明显高于CH4的扩散系数。同时,样品的孔隙直径越大反而不利于气体在煤层中进行Knudsen型扩散。

    图6图7显示,研究区样品中气体的过渡型扩散系数与压力变化的关系。

    图  7  研究区样品中CH4的过渡型扩散系数随压力变化
    Figure  7.  Change diagram of the transition-type diffusion coefficient of CH4 along with the pressure in the study area samples

    图6图7可以看出,当温度条件48 ℃,压力条件20 MPa时,刘庄矿13煤样品中CO2的过渡型扩散系数为0.61 m2/s,CH4的扩散系数为0.69 m2/s;祁东矿7煤样品中CO2的扩散系数为0.72 m2/s,CH4的扩散系数为0.82 m2/s。当温度条件48 ℃,压力条件2 MPa时,刘庄矿13煤样品中CO2的扩散系数为1.83 m2/s,CH4的扩散系数为1.28 m2/s;祁东矿7煤样品中CO2的扩散系数为2.17 m2/s,CH4的扩散系数为1.52 m2/s。

    对比图6图7,能够明显看出气体类型对于过渡型扩散系数的影响,分子直径更小的气体,即CO2在压力较低的条件下能够有更大的扩散系数。但当压力开始下降时,分子直径更小的气体其敏感性更高,下降速率在压力2~10 MPa区间内尤其明显,而当压力达到10 MPa以后下降速率降低,压力达到20 Mpa后基本趋于平缓。当压力达到一定值时,气体分子直径小的过渡型扩散系数反而会低于气体分子直径更大气体。

    本文扩散速率变化特征以CH4气体为研究对象,温度选取42 ℃以及48 ℃,压力条件为2~20 MPa。

    气体由煤基质系统向裂隙系统扩散过程符合Fick第一定律,其中扩散通量即扩散速率(cm3/(g·s))表达式[25]可表示为:

    $$ {{{q}}_{\text{m}}} = \frac{1}{{D\sigma }}[C({{p}}) - C] $$ (8)

    平衡条件下基质边缘气体质量浓度(cm3/g),可表示为[25]

    $$ C({{p}}) = ({C_P} + {C_r}) $$ (9)

    煤基质孔隙系统中吸附气体质量浓度(cm3/g),可表示为[25]

    $$ {C_{\text{p}}} = \frac{{{V_L} p}}{{{{{p}}_L} + p}} $$ (10)

    游离气质量浓度(cm3/g)可表示为[25]

    $$ {C_{\text{r}}} = \frac{{{V_{{\rm{STP}}}}}}{{\text{m}}} $$ (11)

    所计算气体体积受温度及压力影响,需要将气体校正到标准状态下,才能进行同一尺度的比较和计算,标准状态下的气体体积(cm3),可表示为[25]

    $$ {V_{{\rm{STP}}}} = \frac{{{T_{}} {P_0} {V_{\text{m}}}}}{{{P_0} T}} $$ (12)

    将扩散系数D与基质形状因子σ结合起来,定义一个时间常数ε,如下[17]

    $$ \varepsilon = \frac{1}{{D\sigma }} $$ (13)

    完全推导后的扩散通量可表示为[17]

    $$ {{{q}}_{\rm{m}}} = \frac{1}{\varepsilon }[C({{p}}) - C] ^{ } $$ (14)

    式中,$ D $为扩散系数,m2/s;$ \sigma $为基质形状因子,m-2;$ C $为煤中气体平均质量浓度,cm3/g;$ {C}_{\left(p\right)} $为煤颗粒在压力p达到吸附平衡时,煤基质边界气体质量浓度,cm3/g;

    $ {C}_{p} $为煤基质孔隙系统中吸附气体质量浓度,cm3/g;$ {C}_{r} $为游离气体质量浓度,cm3/g;$ {V}_{STP} $为标准状态下气体体积,cm3;$ {P}_{0} $为标准大气压,101.33 KPa;$ {P}_{{\rm{m}}} $为介质压力,MPa;$ {P}_{{\rm{L}}} $为朗格缪尔压力,MPa;$ {V}_{{\rm{L}}} $为朗格缪尔体积,cm3/g;$ {V}_{{\rm{m}}} $为孔隙体积,cm3;$ m $为样品质量,g。

    形状因子选取根据基质形状,通过表7选取。

    表  7  煤基质单元形态因子[26]
    Table  7.  Coal matrix unit morphfactor
    基质单元几何形状特征参数形态因子
    板状板宽L(π/L)2
    柱状底面半径R5.78/R2
    球状球体半径R(π/R)2
    下载: 导出CSV 
    | 显示表格

    通过式(8)—式(13),取温度条件42 、48 ℃,压力条件2 ~20 MPa,计算CH4气体在刘庄矿13煤和祁东矿7煤样品中不同孔隙类型里的扩散速率,即气体不同扩散类型的扩散速率,结果如图810所示。

    图  8  研究区样品中CH4的Fick型扩散速率随压力变化
    Figure  8.  Change diagram of the Fick-type diffusion rate of CH4 with pressure in the study area samples
    图  9  研究区样品中CH4的Knudsen型扩散速率随压力变化
    Figure  9.  Change diagram of the Knudsen-type diffusion coefficient of CH4 in the study area sample along with the pressure
    图  10  研究区样品中CH4的过渡型扩散系数随压力变化
    Figure  10.  Change diagram of the transition-type diffusion coefficient of CH4 along with the pressure in the study area samples

    图810可以看出,温度条件48 ℃,压力条件20 MPa下,刘庄矿13煤样品中Fick型扩散类型的扩散速率为0.370 cm3/(g·s),祁东矿7煤样品中Fick型扩散类型的扩散速率为0.370 cm3/(g·s);刘庄矿13煤样品中Knudsen型扩散类型的扩散速率为0.310 cm3/(g·s),祁东矿7煤样品中Knudsen型扩散类型的扩散速率为0.277 cm3/(g·s);刘庄矿13煤样品中过渡型扩散类型的扩散速率为0.161 cm3/(g·s),祁东矿7煤样品中过渡型扩散类型的扩散速率为0.152 cm3/(g·s)。CH4气体在煤基质中Fick型扩散的扩散速率明显高于Knudsen型扩散的扩散速率,过渡型扩散的扩散速率最低,表明在刘庄矿13煤和祁东矿7煤样品中,CH4气体扩散类型以Fick型扩散为主时能够够较快的扩散,以过渡型扩散为扩散类型为主时会降低扩散速度。

    同时,还可以看出,42和48 ℃扩散速率随压力变化曲线基本重合,表明温度对于扩散速率的影响并不大。

    根据表6的数据显示,储层条件(24~48 ℃、2~20 MPa)下,刘庄矿与祁东矿的煤样中的CH4与CO2气体均以Fick型扩散为主,而根据图37的气体扩散系数的数据显示,CH4相比CO2在Fick型扩散中更具优势,但在Knudsen型扩散和过渡型扩散中CO2气体扩散系数要高于CH4

    在CO2-ECBM过程中,CO2的扩散路径为压裂裂缝-宏观孔隙-内生裂隙-显微孔隙-大孔-中孔-微孔,CH4的扩散路径刚好相反,CO2在微孔中的扩散能够快速扩散,不仅能够提高CO2的注气速度,还能够更快地同CH4进行置换,提高CH4的产量。

    1)刘庄矿13煤和祁东矿7煤样品,孔隙分布中以微孔及小孔孔容占比最高,储层条件下2个煤样中气体均主要以Fick型扩散为主,过渡型扩散次之,Knudsen型扩散占比最低。

    2)温度对于Fick型扩散系数以及过渡扩散类型系数的影响远远小于压力的影响,2种扩散类型的扩散系数都随压力的增加而降低,且气体扩散系数在0~5 Mpa内随压力下降变化速率最大。当压力值达到20 MPa后,2种扩散类型的扩散系数基本都趋于平衡,相较于2 MPa时的扩散系数,缩小了近10倍。

    3)温度对于Knudsen型扩散系数的影响比较明显地高于其他2种扩散类型,温度每增加6 ℃,其扩散系数增加0.2 m2/s,Knudsen型扩散系数基本呈较缓的线性增长。

    4)CH4相比CO2在Fick型扩散中更具优势,但在Knudsen型扩散和过渡型扩散中CO2气体扩散系数要高于CH4

    5)CH4气体以Fick型扩散或Knudsen型扩散为扩散类型时能够较快的在煤孔隙中扩散,以过渡型扩散为扩散类型为主时扩散速度会较慢;

    6)CO2-ECBM过程中,CO2在微孔及小孔中相比CH4更具有扩散优势,增加煤层的微孔及小孔占比能够提高CO2的注气速度以及CH4的产量,增强CO2-ECBM的产气效果,对实现CO2-ECBM技术突破提供理论依据。

  • 图  1   刘庄煤矿构造纲要图[22]

    Figure  1.   Structural outline map of Liuzhuang Coal Mine[22]

    图  2   祁东煤矿构造纲要图[23]

    Figure  2.   Structural outline map of Qidong Coal Mine[23]

    图  3   研究区样品中CO2的Fick型扩散系数随压力变化

    Figure  3.   Change diagram of the Fick-type diffusion coefficient of CO2 in the study area sample along with the pressure

    图  4   研究区样品中CH4的Fick型扩散系数随压力变化

    Figure  4.   Change diagram of the Fick-type diffusion coefficient of CH4 in the study area sample along with the pressure

    图  5   研究区样品中气体的Knudsen型扩散系数随温度变化

    Figure  5.   Change diagram of Knudsen-type diffusion coefficient of gas in the study area

    图  6   研究区样品中CO2的过渡型扩散系数随压力变化

    Figure  6.   Change diagram of the transition-type diffusion coefficient of CO2 along with the pressure in the study area samples

    图  7   研究区样品中CH4的过渡型扩散系数随压力变化

    Figure  7.   Change diagram of the transition-type diffusion coefficient of CH4 along with the pressure in the study area samples

    图  8   研究区样品中CH4的Fick型扩散速率随压力变化

    Figure  8.   Change diagram of the Fick-type diffusion rate of CH4 with pressure in the study area samples

    图  9   研究区样品中CH4的Knudsen型扩散速率随压力变化

    Figure  9.   Change diagram of the Knudsen-type diffusion coefficient of CH4 in the study area sample along with the pressure

    图  10   研究区样品中CH4的过渡型扩散系数随压力变化

    Figure  10.   Change diagram of the transition-type diffusion coefficient of CH4 along with the pressure in the study area samples

    表  1   煤样工业分析表

    Table  1   Analysis table of coal sample industry

    样品名刘庄13煤祁东7煤
    煤类气煤肥煤
    采样深度/m7631 044
    水分Mad/%1.250.88
    挥发分Vdaf/%39.3438.59
    固定碳FCad/%39.8149.51
    灰分Aad/%19.6011.02
    平均真相密度/(g·cm-31.361.45
    镜质组反射率RO/%0.72520.7558
    下载: 导出CSV

    表  2   扩散类型及其对应诺森数

    Table  2   Diffusion type and its counterpart

    扩散类型Kn
    Fick型扩散Kn≥10
    过渡型扩散0.1≤Kn≤10
    Knudsen型扩散Kn≤0.1
    下载: 导出CSV

    表  3   模拟不同埋深条件下的气体分子自由程

    Table  3   Simate the molecular free path of gas under different burial depths

    模拟埋深

    条件/m
    温度T/K压力p/MPa分子有效直径d0/nm平均分子自由程λ/nm
    CH4CO2CH4CO2
    −1 500321.15150.3800.3300.4610.611
    −1 200315.15120.5650.749
    −900309.1590.7390.980
    −600303.1561.0871.441
    −300297.1561.0651.413
    下载: 导出CSV

    表  4   不同扩散类型的煤体孔分布情况

    Table  4   Distribution of coal body holes in different diffusion types

    埋深条件
    /m
    扩散孔径分布/nm
    Knudsen型扩散过渡型扩散Fick型扩散
    CH4CO2CH4CO2CH4CO2
    −1 5000.0460.0610.046~4.610.061~6.114.616.11
    −1 2000.0570.0750.057~5.70.075~7.455.657.49
    −9000.0740.0980.074~7.390.098~9.807.399.80
    −6000.1080.1440.108~10.870.144~14.4210.8714.42
    −3000.1070.1410.107~10.650.141~14.1310.6514.13
    下载: 导出CSV

    表  5   气体分子平均自由程结果表

    Table  5   Table of the average free course results of the gas molecules

    气体种类λ/nm
    CO20.424~4.580
    CH40.295~3.454
    下载: 导出CSV

    表  6   不同孔径中CO2和CH4气体的不同扩散类型占比

    Table  6   Results of different diffusion types of CO2 and CH4 gas in different pore sizes

    气体类型煤岩样品总孔容/(cm3·g−1孔隙体积分布/%以诺森数分类计/%
    微孔
    (<10 nm)
    小孔
    (10~102 nm)
    中孔
    (102~103 nm)
    大孔
    (103~104 nm)
    Fick型扩散过渡型
    扩散
    Knudsen型
    扩散
    CO2刘庄0.01534.78038.35310.28016.58852.43630.17417.390
    祁东0.01936.45737.96610.93414.64450.88830.88418.228
    CH4刘庄0.01534.78038.35310.28016.58856.69732.86910.434
    祁东0.01936.45737.96610.93414.64455.10633.95710.937
    下载: 导出CSV

    表  7   煤基质单元形态因子[26]

    Table  7   Coal matrix unit morphfactor

    基质单元几何形状特征参数形态因子
    板状板宽L(π/L)2
    柱状底面半径R5.78/R2
    球状球体半径R(π/R)2
    下载: 导出CSV
  • [1] 邢万丽. 煤中CO2、CH4、N2及多元气体吸附/解吸、扩散特性研究[D]. 大连: 大连理工大学, 2016.

    XING Wanli. Study on adsorption / desorption and diffusion characteristics of CO2, CH4, N2 and multicomponent gases in coal[D]. Dalian: Dalian University of Technology, 2016.

    [2]

    RUCKENSTEINE,VAIDYANATHAN S A,YOUNGQUIST R G. Sorption by solids with bidisperse pore structures[J]. Chemical Engineering Science,1971,26:1305−1318. doi: 10.1016/0009-2509(71)80051-9

    [3] 张时音,桑树勋. 不同煤级煤层气吸附扩散系数分析[J]. 中国煤炭地质,2009,21(3):24−27. doi: 10.3969/j.issn.1674-1803.2009.03.007

    ZHANG Shiyin,SANG Shuxun. Adsorption-diffusion coefficient analysis of Coal-bed methane in different rank coals[J]. Coal Geology of China,2009,21(3):24−27. doi: 10.3969/j.issn.1674-1803.2009.03.007

    [4] 张登峰,崔永君,李松庚,等. 甲烷及二氧化碳在不同煤阶煤内部的吸附扩散行为[J]. 煤炭学报,2011,36(10):1693−1698. doi: 10.13225/j.cnki.jccs.2011.10.001

    ZHANG Dengfeng,CUI Yongjun,LI Songgeng,et al. Adsorption and diffusion behaviors of methane and carbon dioxide on various rank coals[J]. Journal of China Coal Society,2011,36(10):1693−1698. doi: 10.13225/j.cnki.jccs.2011.10.001

    [5]

    LI Xiangchun,NIE Baisheng,ZHANG Ruming,et al. Experiment of gas diffusion and its diffusion mechanism in coal[J]. International Journal of Mining Science and Technology,2012,22(6):885−889. doi: 10.1016/j.ijmst.2012.12.007

    [6]

    MENG Ya,LI Zhiping. Experimental study on diffusion property of methane gas in coal and its influencing factors[J]. Fuel,2016,185:219−228. doi: 10.1016/j.fuel.2016.07.119

    [7]

    LI Guoqing,MENG Zhaoping. A preliminary investigation of CH4 diffusion through gas shale in the Paleozoic Longmaxi Formation, southern Sichuan basin, China[J]. Journal of Natural Gas Science and Engineering,2016,36:1220−1227. doi: 10.1016/j.jngse.2016.05.069

    [8]

    WANG Yucang,XUE Sheng,XIE Jun. A general solution and approximation for the diffusion of gas in a spherical coal sample[J]. International Journal of Mining Science Technology,2014,24(3):345−348. doi: 10.1016/j.ijmst.2014.03.010

    [9]

    MENG Ya,LI Zhiping. Experimental comparisons of gas adsorption, sorption induced strain, diffusivity and permeability for low and high rank coals[J]. Fuel,2018,234:914−923. doi: 10.1016/j.fuel.2018.07.141

    [10]

    CHARRIERE D,POKRYSZKA Z,BEHRA P. Effect of pressure and temperature on diffusion of CO2 and CH4 into coal from the Lorraine basin (France)[J]. International Journal of Coal Geology,2010,81(4):373−380. doi: 10.1016/j.coal.2009.03.007

    [11]

    BUSCH A,GENSTERBLUM Y. CBM and CO2-ECBM related sorption processes in coal: A review[J]. International Journal of Coal Geology,2011,87(2):49−71. doi: 10.1016/j.coal.2011.04.011

    [12]

    STAIB G,SAKUROVS R,GRAY E M A. Dispersive diffusion of gases in coals. Part II: An assessment of previously proposed physical mechanisms of diffusion in coal[J]. Fuel,2015,143:620−629. doi: 10.1016/j.fuel.2014.11.087

    [13] 李志强,段振伟,景国勋. 不同温度下煤粒瓦斯扩散特性试验研究与数值模拟[J]. 中国安全科学学报,2012,22(4):38−42. doi: 10.3969/j.issn.1003-3033.2012.04.007

    LI Zhiqiang,DUAN Zhenwei,JING Guoxun. Experimental study on gas diffusion characteristics from coal at different temperatures and their numerical simulation[J]. China Safety Science Journal,2012,22(4):38−42. doi: 10.3969/j.issn.1003-3033.2012.04.007

    [14]

    CIEMBRONIEWICZ A,MARECKA A. Kinetics of CO2 sorption for two Polish hard coals[J]. Fuel,1993,72(3):405−408. doi: 10.1016/0016-2361(93)90062-7

    [15]

    MENG Ya,LI Zhiping,LAI Fengpeng. Experimental study on porosity and permeability of anthracite coal under different stresses[J]. Journal of Petroleum Science and Engineering,2015,133:810−817. doi: 10.1016/j.petrol.2015.04.012

    [16]

    PAN Z J,CONNELL L D. Modelling permeability for coal reservoirs: A review of analytical models and testing data[J]. International Journal of Coal Geology,2012,92:1−44. doi: 10.1016/j.coal.2011.12.009

    [17] 孟召平,张贵元,刘金融,等. 低煤阶煤中甲烷扩散性能分析[J]. 煤田地质与勘探,2019,47(2):84−89. doi: 10.3969/j.issn.1001-1986.2019.02.014

    MENG Zhaoping,ZHANG Guiyuan,LIU Jinrong,et al. Analysis of diffusion properties of methane in low rank coal[J]. Coal Geology & Exploration,2019,47(2):84−89. doi: 10.3969/j.issn.1001-1986.2019.02.014

    [18] 韦重韬,周荣福. 煤层气多煤层扩散逸失地质历史模型及数值模拟[J]. 高校地质学报,2003,9(3):390−395. doi: 10.3969/j.issn.1006-7493.2003.03.008

    WEI Chongtao,ZHOU Rongfu. Geological history model and its numerical simulation of coalbed methane emission for multi-Coal seams[J]. Ceological Journal of China Universities,2003,9(3):390−395. doi: 10.3969/j.issn.1006-7493.2003.03.008

    [19] 何学秋,聂百胜. 孔隙气体在煤层中扩散的机理[J]. 中国矿业大学学报,2001,30(1):3−6. doi: 10.3321/j.issn:1000-1964.2001.01.001

    HE Xueqiu,NIE Baisheng. Diffusion mechanism of porous gases in coal seams[J]. Journal of China University of Mining & Technology,2001,30(1):3−6. doi: 10.3321/j.issn:1000-1964.2001.01.001

    [20] 聂百胜,何学秋,王恩元. 瓦斯气体在煤层中的扩散机理及模式[J]. 中国安全科学学报,2000,10(6):27−31. doi: 10.3969/j.issn.1003-3033.2000.06.006

    NIE Baisheng,HE Xueqiu,WANG Enyuan. Mechanism and modes of gas diffusion in coal seams[J]. China Safety Science Journal,2000,10(6):27−31. doi: 10.3969/j.issn.1003-3033.2000.06.006

    [21] 曹成润,牛 伟,张遂安等. 煤层气在煤储层中的扩散及其影响因素[J]. 世界地质,2004,9(3):266−269. doi: 10.3969/j.issn.1004-5589.2004.03.011

    CAO Chengrun,NIU Wei,ZHANG Suian,et al. Diffusion and its influences of coalbed methane in coal reservoirs[J]. Global Geology,2004,9(3):266−269. doi: 10.3969/j.issn.1004-5589.2004.03.011

    [22] 张磊. 刘庄煤矿突水水害危险性评价与涌水量预测[D]. 淮南: 安徽理工大学, 2021. 000653.

    ZHANG Lei. Water inrush risk assessment and water inflow prediction in Liuzhuang Coal Mine[D]. Huainan: Anhui University of Science and Technology , 2021, 000653.

    [23] 田强国. 祁东煤矿地质构造对煤与瓦斯突出的影响研究[D]. 淮南: 安徽理工大学, 2009.

    TIAN Guoqiang. The Study on The Influence of Geologic Structure on Gas Outburst in Qidong Coal Mine[D]. Huainan : Anhui University of Science and Technology , 2009.

    [24] 闫宝珍,王延斌,倪小明. 地层条件下基于纳米级孔隙的煤层气扩散特征[J]. 煤炭学报,2008,33(6):657−660. doi: 10.3321/j.issn:0253-9993.2008.06.013

    YAN Baozhen,WANG YanBin,NI Xiaoming. Coalbed methane diffusion characters based on nano-scaled pores under fromation conditions[J]. Journal of China Coal Society,2008,33(6):657−660. doi: 10.3321/j.issn:0253-9993.2008.06.013

    [25] 杨兆彪,秦 勇,王兆丰,等. 钻井液条件下煤芯煤层气解吸-扩散模型及逸散量求取[J]. 中国科学:地球科学,2010,40(2):171−177.

    YANG Zhaobiao,QIN Yong,WANG Zhaofeng,et al. Desorption-diffusion model and lost gas quantity estimation of coalbed methane from coal core under drilling fluid medium[J]. Science in China:Earth Sciences,2010,40(2):171−177.

    [26]

    PRICE H S,MCCULOCH R C,EDWARDS J C,et al. A computer model study of methane migration in coal beds[J]. Canad Min Metallurg Bull,1973,66:103−112.

  • 期刊类型引用(1)

    1. 韩文彬,张迪,郑义,刘东洋. 煤层裂隙粗糙度与瓦斯开采效率耦合数值模型. 煤矿安全. 2024(08): 51-61 . 百度学术

    其他类型引用(0)

图(10)  /  表(7)
计量
  • 文章访问数:  55
  • HTML全文浏览量:  2
  • PDF下载量:  28
  • 被引次数: 1
出版历程
  • 收稿日期:  2022-10-19
  • 网络出版日期:  2023-08-13
  • 刊出日期:  2023-05-31

目录

/

返回文章
返回