• 中文核心期刊
  • 中国科技核心期刊
  • RCCSE中国核心学术期刊

基于负压衰减效应的抽采长钻孔残余瓦斯分布规律研究

蔡春林, 江成玉, 苟仁涛, 夏勋鹏, 王沉, 刘文昌

蔡春林,江成玉,苟仁涛,等. 基于负压衰减效应的抽采长钻孔残余瓦斯分布规律研究[J]. 煤矿安全,2024,55(11):61−68. DOI: 10.13347/j.cnki.mkaq.20230150
引用本文: 蔡春林,江成玉,苟仁涛,等. 基于负压衰减效应的抽采长钻孔残余瓦斯分布规律研究[J]. 煤矿安全,2024,55(11):61−68. DOI: 10.13347/j.cnki.mkaq.20230150
CAI Chunlin, JIANG Chengyu, GOU Rentao, et al. Study on residual gas distribution in long extraction borehole based on negative pressure attenuation effect[J]. Safety in Coal Mines, 2024, 55(11): 61−68. DOI: 10.13347/j.cnki.mkaq.20230150
Citation: CAI Chunlin, JIANG Chengyu, GOU Rentao, et al. Study on residual gas distribution in long extraction borehole based on negative pressure attenuation effect[J]. Safety in Coal Mines, 2024, 55(11): 61−68. DOI: 10.13347/j.cnki.mkaq.20230150

基于负压衰减效应的抽采长钻孔残余瓦斯分布规律研究

基金项目: 国家自然科学基金资助项目(52174072);贵州省科技计划资助项目(黔科合支撑[2021]一般354,黔科合基础-ZK[2021]一般264)
详细信息
    作者简介:

    蔡春林(1998—),男,贵州遵义人,硕士研究生,研究方向为矿山安全与灾害防治。E-mail:1461318562@qq.com

    通讯作者:

    江成玉(1980—),男,贵州仁怀人,副教授,硕士研究生导师,硕士,主要从事煤矿安全开采技术和瓦斯灾害防治等方面的教学与研究工作。E-mail:cyjiang1981@163.com

  • 中图分类号: TD712.6

Study on residual gas distribution in long extraction borehole based on negative pressure attenuation effect

  • 摘要:

    为研究顺层长钻孔抽采下煤层瓦斯残余含量的分布规律,基于负压衰减效应推导了长钻孔瓦斯抽采模型公式,利用COMSOL模拟软件建立了瓦斯−空气模型,分析了顺层长钻孔抽采下煤层瓦斯残余含量分布规律。结果表明:抽采负压呈现对数函数分布,孔底处负压最小,随着钻孔孔深增加,抽采负压衰减越小;由于钻孔存在抽采负压衰减现象,钻孔横向和纵向残余瓦斯含量呈现对数函数分布;在封孔段末端煤体受到封孔材料压力作用,使得煤体压实,孔隙率减小,瓦斯在附近运移更加困难,瓦斯流动阻力增大,致使封孔段末端瓦斯含量变化存在“滞后”现象;通过在贵州林华煤矿测量1#、2#钻孔不同孔深处煤层残余瓦斯含量,获得了取样深度与残余瓦斯含量的对数函数关系,拟合度分别为0.94300.9572,与理论分析结果相吻合。

    Abstract:

    In order to study the distribution law of gas residual content in coal seam under long borehole extraction, the formula of long borehole gas extraction model was deduced based on negative pressure attenuation effect. The gas-air model was established by COMSOL simulation software, and the distribution law of gas residual content in coal seam under long borehole extraction was analyzed. The results show that the extraction negative pressure presents a logarithmic function distribution, and the negative pressure at the bottom of the hole is the smallest. As the depth of the borehole increases, the attenuation of the extraction negative pressure is smaller. Due to the phenomenon of negative pressure attenuation in drilling, the horizontal and vertical residual gas content of drilling presents logarithmic function distribution. At the end of the sealing section, the coal body is subjected to the pressure of the sealing material, which makes the coal body compacted, the porosity decreases, the gas migration in the vicinity is more difficult, and the gas flow resistance increases, resulting in a “lag” phenomenon in the gas content change at the end of the sealing section. By measuring the residual gas content of coal seams at different depths of 1# and 2# boreholes in Linhua Coal Mine, Guizhou, the logarithmic function relationship between sampling depth and residual gas content was obtained. The fitting degrees are 0.943 0 and 0.957 2, respectively, which were consistent with the theoretical analysis results.

  • 煤炭是我国主要能源和重要的工业原料,煤炭安全稳定供应关系到国民经济发展全局,但煤矿瓦斯灾害长期制约着煤炭产业健康发展[1-3]。长钻孔是治理瓦斯的有效手段之一,与传统钻孔方式相比较,具有抽采体积分数大、抽放时间长等优点[4-7]。但并非钻孔长度越长,瓦斯抽采效果越好。由于钻孔长度过长,在瓦斯流动过程中存在负压衰减现象,对钻孔周围瓦斯流动造成影响,进而影响瓦斯连续高效抽采,造成残余瓦斯分布不均[8-9]。因此,研究负压衰减效应下长钻孔抽采时的煤层残余瓦斯分布规律具有重要意义。

    众多学者对这一问题进行了大量研究,取得了较为丰富的研究成果。江成等[10]分析了抽采过程邻近层瓦斯压力演化与煤层变形关系,得出邻近层瓦斯压力衰减呈近似直线;李路广等[11]为了提高矿井瓦斯抽采钻孔施工效率和煤层瓦斯的抽采效果,依据千米定向钻机轨迹可控、过程可溯的特点,通过工业试验实现了对煤与瓦斯突出风险的有效管控,提高了瓦斯抽采效率;冀超辉等[12]在乌冬煤矿进行了特厚急倾斜煤层抽采应用,建立了一套针对特厚急倾斜煤层瓦斯抽采技术;翟成等[13]在祁东煤矿采用高位钻孔抽放技术解决了瓦斯超限问题;毕慧杰等[14]采用理论计算和数值模拟的方法研究了采动裂隙演化规律,计算出合理的施工位置,达到提高了抽采纯量的效果;李宏等[15]为进一步提高采空区裂隙带瓦斯抽采效率,提出了一种大直径顶板定向长钻孔瓦斯抽采技术,取得了经济、高效的现场应用效果;王广宏[16]根据榆树田煤矿实际情况建立抽采流固模型,分析了钻孔负压、长度等钻孔参数对瓦斯抽采的影响;徐超等[17]基于COMSOL数值模拟软件和理论计算方法,建立煤层瓦斯扩散—渗流模型,分析了钻孔内各区段沿程阻力损失和抽采钻孔内负压衰减的影响因素;王兆丰等[18]、李杰[19]、宋建军[20]通过理论分析和现场试验,对顺层钻孔孔内负压分布规律也开展了大量的研究工作。

    综上所述,目前针对负压衰减效应下长钻孔抽采时煤层瓦斯残余量分布规律的研究较少。因此,采用数值模拟方法,结合瓦斯-空气模型与钻孔负压衰减方程,阐明长钻孔抽采下煤层残余瓦斯分布规律,以期为优化瓦斯抽采钻孔布置、提高瓦斯抽采效率等方面提供理论支撑。

    在理想状态煤体内部瓦斯压力下降,煤体基质将会产生收缩变形,煤体孔隙体积增大。根据相关研究可知,煤体发生瓦斯解吸收缩应变εp计算如下[21]

    $$ {\varepsilon _p} = \frac{{4ac{\rho _{\mathrm{v}}}RT}}{{9{V_{\mathrm{m}}}{K_{\mathrm{m}}}}}\ln \left( {1 + b{p_{\mathrm{m}}}} \right) $$ (1)

    式中:a为吸附常数,m3/t;b为吸附常数,MPa−1c为煤可燃部分质量,kg;ρv为煤的视密度,kg/m3R为理想气体常数;T为温度,K;Vm为标准气体摩尔体积,L/mol;Km为基质弹性模量,GPa;pm为煤基质瓦斯压力,MPa。

    煤体基质吸收瓦斯收缩变形为可逆过程,当煤体瓦斯压力由初始瓦斯压力p0降至pm,基质收缩应变$ \Delta {\varepsilon _p} $为:

    $$ \Delta {\varepsilon _p} = {\varepsilon _{{p_0}}} - {\varepsilon _{{p_{\text{m}}}}} = \frac{{4ac{\rho _{\mathrm{v}}}RT}}{{9{V_{\mathrm{m}}}{K_{\mathrm{m}}}}}\ln \left( {\frac{{1 + b{p_0}}}{{1 + b{p_{\mathrm{m}}}}}} \right) $$ (2)

    式中:$ {\varepsilon _{p_0}} $为初始应变;$ {\varepsilon _{p_{\mathrm{m}}}} $为压力为pm时的应变。

    假设瓦斯解吸不影响煤体体积,将式(2)代入孔隙率计算公式,得到瓦斯解吸后煤体孔隙率φm为:

    $$ {\varphi _{\mathrm{m}}} = {\varphi _0} + \Delta {\varepsilon _p} = {\varphi _0} + \frac{{4ac{\rho _{\mathrm{v}}}RT}}{{9{V_{\mathrm{m}}}{K_{\mathrm{m}}}}}\ln \left( {\frac{{1 + b{p_0}}}{{1 + b{p_{\mathrm{m}}}}}} \right) $$ (3)

    式中:φ0为初始孔隙率。

    Kozeny-Carman方程关于孔隙率与渗透性的方程为[22]

    $$ k = \frac{{\rho g{d^2}{\varphi ^3}}}{{180\mu {{\left( {1 - \varphi } \right)}^2}}} $$ (4)

    式中:k为动态渗透率,m²;g为重力常量,m/s²;ρ为煤的密度,kg/m³;d为固体颗粒平均直径,m;μ为气体黏度,Pa·s;$\varphi $为孔隙率。

    将式(3)代入式(4)可以得到基质渗透率km为:

    $$ {k_{\mathrm{m}}} = {k_0}{\left\{ {1 + \frac{{4ac{\rho _{\mathrm{v}}}RT}}{{9{V_{\mathrm{m}}}{K_{\mathrm{m}}}}}\ln \left( {\frac{{1 + b{p_0}}}{{1 + b{p_{\mathrm{m}}}}}} \right)} \right\}^3} $$ (5)

    式中:k0为初始渗透率,m2

    考虑有效应力对孔隙率的影响,双重介质中的孔隙率表达式可定义为:

    $$ {\varphi _{\mathrm{m}}} = \frac{{{V_{\mathrm{p}}}}}{{{V_{\mathrm{t}}}}} = \frac{{{V_{{\mathrm{p}}0}} - \Delta V}}{{{V_{{\mathrm{t}}0}} - \Delta V}} = \frac{{{\varphi _0} - {\varepsilon _{\mathrm{v}}}}}{{1 - {\varepsilon _{\mathrm{v}}}}} $$ (6)

    式中:Vp为煤变形前后的孔隙体积,m3Vt为孔隙体积,m3Vp0为煤变形前的孔隙体积,m3Vt0为初始孔隙体积,m3ɛv为体积应变;ΔV为孔隙体积变化量,m3

    煤体瓦斯压力由p0降至p1,空气压力由0升至p2,将应变代入式(5),得到多孔介质孔隙变化方程:

    $$ \begin{split} & \begin{gathered} \varphi = {\varphi _0} - \frac{{a\left( {1 - {\varphi _0}} \right)\left( {{p_0} - {p_1} + {p_2}} \right)}}{{{K_{\mathrm{f}}} - a\left( {{p_0} - {p_1} + {p_2}} \right)}} +\\ \frac{{4ac{\rho _{\mathrm{v}}}RT\left[ {\ln \left( {1 + b{p_0}} \right) - \ln \left( {1 + b{p_{\mathrm{m}}}} \right)} \right]}}{{9{V_{\mathrm{m}}}{K_{\mathrm{m}}}\left[ {{K_{\mathrm{f}}} - a\left( {{p_0} - {p_1} + {p_2}} \right)} \right]}} \end{gathered} \\[-14pt]& \end{split} $$ (7)

    考虑基质收缩与有效应力,得到基质渗透率动态演化方程:

    $$\begin{split} & \begin{gathered} k = {k_0}{\left( {\frac{\varphi }{{{\varphi _0}}}} \right)^3} = {k_0}\left\{ {1 - \frac{{a\left( {1 - {\varphi _0}} \right)\left( {{p_0} - {p_1} + {p_2}} \right)}}{{{K_{\mathrm{f}}} - a\left( {{p_0} - {p_1} + {p_2}} \right)}}}+ \right. \\ {\left. { \frac{{4ac{\rho _{\mathrm{v}}}RT\left[ {\ln \left( {1 + b{p_0}} \right) - \ln \left( {1 + b{p_{\mathrm{m}}}} \right)} \right]}}{{9{V_{\mathrm{m}}}{K_{\mathrm{m}}}\left[ {{K_{\mathrm{f}}} - a\left( {{p_0} - {p_1} + {p_2}} \right)} \right]}}} \right\}^3} \end{gathered} \\[-14pt] & \end{split} $$ (8)

    瓦斯−空气双气体系统中,只考虑气体在多孔介质扩散,控制方程与单一气体在多孔物质相似。其物质交换方程可表示为[23]

    $$ {Q_{\mathrm{m}}} = {D_t}\tau \left( {{c_{\mathrm{m}}} - {\rho _{\mathrm{f}}}} \right) = {D_t}\frac{{3{\pi ^2}}}{{{L^2}}}\left( {{c_{\mathrm{m}}} - {\rho _{\mathrm{f}}}} \right) $$ (9)

    式中:Qm为基质与孔隙裂隙的气体交换量,m³/s;Dt为流体扩散系数,m2/s;τ为形状因子;cmρf分别为基质与孔隙裂隙中的流体密度,kg/m3L为基质的孔径,m。

    流体扩散系数可表示为:

    $$ {D_t} = {D_0}{{\mathrm{e}}^{ - \xi t}} + {D_{\mathrm{r}}} $$ (10)

    式中:D0为气体初始扩散系数,m2/s;ξ为衰减系数;t为时间,s;Dr为剩余气体扩散系数,m2/s。

    将气体压强公式与式(10)代入式(9),得到质量交换方程的控制方程:

    $$ {Q_{\mathrm{m}}} = \frac{{3{\pi ^2}{M_{\mathrm{c}}}\left( {{p_{\mathrm{n}}} - {p_{\mathrm{f}}}} \right)\left[ {{D_0}{{\mathrm{e}}^{ - \xi t}} + {D_{\mathrm{r}}}} \right]}}{{{L^2}RT}} $$ (11)

    式中:Mc为甲烷摩尔质量,kg/mol;pn为基质中流体压力,MPa;pf为孔隙中流体压力,MPa。

    单位质量煤基质含气量可表示为:

    $$ {m_{\mathrm{p}}} = \frac{{abp\rho {M_{\mathrm{c}}}}}{{\left( {1 + bp} \right){V_{\mathrm{m}}}}} + {\varphi _{\mathrm{p}}}\frac{{{M_{\mathrm{c}}}p}}{{RT}} $$ (12)

    将式(10)和式(11)代入质量守恒方程得到煤体裂隙瓦斯压力分布方程,即:

    $$ \begin{split} & \begin{gathered} \frac{\partial }{{\partial t}}\left( {\frac{{{M_{\mathrm{c}}}{\rho _{\mathrm{c}}}}}{{{V_{\mathrm{m}}}}} \cdot \frac{{{V_{\mathrm{L}}}{p_{\mathrm{m}}}}}{{{p_{\mathrm{L}}} + {p_{\mathrm{m}}}}} + \phi m\frac{{{M_{\mathrm{c}}}{p_{\mathrm{m}}}}}{{RT}}} \right)= \\ - \frac{{3{\pi ^2}{M_{\mathrm{c}}}\left( {{p_{\mathrm{m}}} - {p_{\mathrm{f}}}} \right)\left[ {{D_0}{{\mathrm{e}}^{ - \xi t}} + {D_{\mathrm{r}}}} \right]}}{{{L^2}RT}} \end{gathered} \\[-14pt]& \end{split} $$ (13)

    随着埋深的增加,煤体受到地应力与瓦斯压力共同作用。煤体变形方程可用Navier-Stokes方程表示[24]

    $$ G{u_{i,jj}} + \frac{G}{{1 - 2v}}{u_{j,ji}} - {\alpha _{\mathrm{m}}}{p_{\mathrm{m}}} - {k_{\mathrm{f}}}{\varepsilon _{s,i}} + {F_i} = 0 $$ (14)

    式中:G为剪切模量,G=E/2(1+υ),GPa;E为弹性模量,GPa;υ为泊松比;kf为裂隙的Biot’s系数;αm为煤基质的Biot’s系数;εs,i为基质最大应变;Fi为体积力分量。

    假设钻孔全部处于煤体中,为方便分析,沿着钻孔长度方向将钻孔分割为若干个控制单元,顺层钻孔控制单元示意图如图1所示。

    图  1  顺层钻孔控制单元示意图
    Figure  1.  Schematic diagram of bedding drilling control unit

    钻孔中煤体压力损失由沿程阻力损失和压降损失组成,其中压降损失由加速度压降和局部压降组成,沿程损失是由流体流动与管道壁面形成的阻力;加速度压降是抽采过程流体质量不断变化,当截面相同,流体流量越大,速度越大。在速度增加过程中伴随部分能量损失;局部压降是由钻孔内部截面发生变化形成是能量损失,如钻孔变形,塌孔等。其理论控制单元压力损失公式为[23]

    $$ \begin{split} & \begin{gathered} \frac{{{\mathrm{d}}p\left( x \right)}}{{{\mathrm{d}}x}} = \Delta {p_{{\mathrm{yc}}}} + \Delta {p_{{\mathrm{yj}}}} = \frac{{2{f_i}\rho {{\left[ {2Q\left( x \right) + q\left( x \right)} \right]}^2}}}{{{\pi ^2}{D^5}}}+ \\ \frac{{16\rho \left[ {{q^2}\left( x \right) + 2Q\left( x \right)q\left( x \right)} \right]}}{{{\pi ^2}{D^4}}} \end{gathered} \\[-14pt]& \end{split} $$ (15)

    式中:pyc为钻孔沿程阻力损失;pyj为钻孔压降;p(x)为距离孔口x处钻孔内部气体压力,MPa;q(x)为距离孔口x处瓦斯涌出量,m³;Q(x)为距离孔口x处钻孔内瓦斯气体流量,m3/s;x为距孔口距离,m;fi为孔壁摩擦修正系数;D为钻孔直径,m。

    钻孔周围煤体残存瓦斯含量与煤层瓦斯压力存在对应关系,而煤体瓦斯压力与煤体受到的压力差有关。越靠近钻孔煤体瓦斯与钻孔形成的压力差越大,煤体瓦斯降低越多,使得煤体吸附-解吸平衡破坏,更多吸附瓦斯解吸为游离瓦斯,煤体瓦斯含量降低。因此,在钻孔横截面瓦斯呈对数函数趋势。由于钻孔存在抽采负压衰减现象,孔口抽采负压低于孔底,因此理论分析钻孔深度残余瓦斯含量也呈现对数函数分布。钻孔附近残余瓦斯含量示意图如图2所示。

    图  2  钻孔附近残余瓦斯含量示意图
    Figure  2.  Diagram of residual gas content near the drill hole

    为了研究负压衰减效应下抽采长钻孔残余瓦斯分布规律,根据工程地质背景,建立数值模型。三维模型尺寸为60 m×300 m ×11 m(宽×深×高),其中,中间层设置3 m厚煤层,在煤层的顶部设置4 m厚顶板,底部设置4 m厚底板。在模型内部x方向14 m处,开挖1条巷道;巷道沿y轴开挖,巷道长60 m,宽4 m。沿x轴在巷道中部设置钻孔,钻孔初始长度为150 m,直径为37.5 mm。模型网格划分示意图如图3所示。

    图  3  模型网格划分
    Figure  3.  Model schematic diagram

    由于贵州地区地表地貌不平整,不同位置上覆岩层厚度不同,受到的垂直应力也不同,为了简化计算,上覆岩层厚度取研究巷道埋深平均值。整个模型采用自由四面体进行划分,选择与达西定律和一般形式偏微分方程的2个物理场初始设置值一致,根据岩层密度简易计算出应在模型上方施加8 MPa垂直应力,模型下方采用固定约束,四周指定位移方向并施加3 MPa水平应力,巷道内部采用自由边界。初始模型中设置20 m为钻孔封孔长度,封孔材料膨胀后将对钻孔周围施加应力,因此在该位置钻孔受到1.5 MPa封孔压力,其余抽采段钻孔周围采用自由边界,并施加13 kPa抽采负压。瓦斯抽采模型参数见表1

    表  1  瓦斯抽采模型参数表
    Table  1.  Gas extraction model parameters table
    输入参数 数值 输入参数 数值
    初始孔隙率/% 3 初始渗透率/m2 1.189 9×10−17
    煤密度/(kg·m−3 1 400 瓦斯动力黏度/(Pa·s) 1.05×10−6
    煤的泊松比 0.29 初始瓦斯压力/MPa 1.4
    煤弹性模量/MPa 2 713 煤层温度/K 293.15
    初始瓦斯扩散系数/
    (m2·s−1
    5.48×10−12 煤体骨架
    弹性模量/MPa
    5 139
    吸附常数b/MPa−1 2 吸附常数a/(m3·t−1 35
    下载: 导出CSV 
    | 显示表格

    不同抽采时间煤体残余瓦斯含量如图4所示,各时间段延钻孔长度方向煤层残存瓦斯含量曲线如图5所示。

    图  4  不同抽采时间煤体残余瓦斯含量
    Figure  4.  Residual gas content of coal body at different extraction times
    图  5  各时间段延钻孔长度方向煤层残存瓦斯含量曲线
    Figure  5.  Curves of residual gas content in coal seams with drilling length at different time periods

    由数值模拟结果可知,随着抽采不断进行,钻孔周围瓦斯含量不断降低,并且范围不断扩大。观测各图钻孔处瓦斯含量,钻孔孔口处瓦斯含量要低于孔底,并且越到抽采后期越明显。这是由于钻孔内部负压受到衰减作用,负压分布不是均匀的,越靠近孔底负压衰减越多,最终造成孔底瓦斯残余含量高于孔口。

    不同测线瓦斯残余含量如图6所示。由图6可知,随着抽采的进行,钻孔附近瓦斯含量不断降低,钻孔深部瓦斯含量也存在降低现象。钻孔初期,抽采段初始位置瓦斯含量与其他位置相差不大,但在抽采后期差距不断增大。由于在封孔段末端煤体受到封孔材料压力,使得煤体压实,孔隙率减小,瓦斯在附近运移更加困难,瓦斯流动阻力增大,使得封孔段末端含量变化存在“滞后”现象。

    图  6  不同测线瓦斯残余含量
    Figure  6.  Gas residual content of different measurement lines

    通过前面理论分析可知,在长钻孔抽采过程存在抽采负压衰减现象,在短钻孔抽采过程中,负压衰减引起的瓦斯残存量差异较小,当钻孔深度达到一定长度时抽采负压衰减成为不可忽略的因素。抽采钻孔由于孔壁的粗糙性以及钻孔不断流入钻孔瓦斯等因素影响,越靠近孔底抽采负压损失越大,由于抽采动力由负压及煤层压力共同决定。根据现场实际工程情况,煤层压力在同一区域内变化不大,因此抽采负压衰减影响抽采效果。靠近封孔段瓦斯损失越小,压力梯度越大,影响范围越广。

    试验选取在贵州林华煤矿进行,该矿位于贵州省金沙县西部,9#煤层位于龙潭组的中部,上段的底部,层位稳定,是井田内厚度最大的中厚煤层,区内主要可采煤层,硫分低,为优质无烟煤,9#煤层平均厚度为2.98 m。为验证模拟抽采后煤体残余瓦斯分布趋势,在该煤矿20915瓦斯治理专用巷瓦斯抽采达标位置钻取2个钻孔,分别测试钻孔40、50、60 m处瓦斯含量,共钻取1#、2# 2个取样孔,取样孔距离抽采孔0.5 m左右,为防止2个取样孔距离太近造成的测试误差,2个取样孔间距为10 m。取样过程采用ZDY4200钻机,配合使用96 mm钻头,取样过程采用风排渣,按照井下瓦斯含量测定要求严格控制取样时间。不同取样深度残余瓦斯含量如图7所示。

    图  7  不同取样深度残余瓦斯含量
    Figure  7.  Residual gas content at different sampling depths

    图7可知,随着取样深度的增加,煤层瓦斯残余含量也随之增加。其中,1#钻孔在深度为40~60 m处,残余瓦斯含量在3.95~4.3 m3/t之间;2# 钻孔在深度为40~60 m处,残余瓦斯含量在4.3~5.7 m3/t之间;2个钻孔的残余瓦斯含量呈对数函数分布,1# 钻孔取样深度与残余量关系为y=1.0373 ln(x)+0.2708,拟合度达到0.943 0;2# 钻孔取样深度与残余瓦斯含量关系为y=2.9321×ln(x)−6.5214,拟合度达到0.9572。现场测试数据与模拟数据存在一定差异,是由于实际情况影响钻孔抽采因素过多,模拟过程忽略了部分因素,并且实际抽采过程负压不断变化,各支管负压也不相同,最终使得各位置残余瓦斯含量不同。总体而言,现场实测结果与理论分析结果相吻合。

    1)钻孔横截面残余瓦斯呈对数函数趋势。由于钻孔存在抽采负压衰减现象,孔口抽采负压低于孔底,因此理论分析钻孔深度残余瓦斯含量也呈现对数函数分布。

    2)由数值模拟可知,在封孔段末端煤体受到封孔材料压力,使得煤体压实,孔隙率减小,瓦斯在附近运移更加困难,瓦斯流动阻力增大,使得封孔段末端含量变化存在“滞后”现象。

    3)现场实测结果表明,钻孔取得的残余瓦斯含量呈对数函数分布,1#钻孔取样深度与残余量关系为y=1.0373 ln(x)+0.1708,拟合度达到0.9430;2#钻孔取样深度与残余量关系为y=2.9321 ln(x)−6.5214,拟合度达到0.9572,与理论分析结果相吻合。

  • 图  1   顺层钻孔控制单元示意图

    Figure  1.   Schematic diagram of bedding drilling control unit

    图  2   钻孔附近残余瓦斯含量示意图

    Figure  2.   Diagram of residual gas content near the drill hole

    图  3   模型网格划分

    Figure  3.   Model schematic diagram

    图  4   不同抽采时间煤体残余瓦斯含量

    Figure  4.   Residual gas content of coal body at different extraction times

    图  5   各时间段延钻孔长度方向煤层残存瓦斯含量曲线

    Figure  5.   Curves of residual gas content in coal seams with drilling length at different time periods

    图  6   不同测线瓦斯残余含量

    Figure  6.   Gas residual content of different measurement lines

    图  7   不同取样深度残余瓦斯含量

    Figure  7.   Residual gas content at different sampling depths

    表  1   瓦斯抽采模型参数表

    Table  1   Gas extraction model parameters table

    输入参数 数值 输入参数 数值
    初始孔隙率/% 3 初始渗透率/m2 1.189 9×10−17
    煤密度/(kg·m−3 1 400 瓦斯动力黏度/(Pa·s) 1.05×10−6
    煤的泊松比 0.29 初始瓦斯压力/MPa 1.4
    煤弹性模量/MPa 2 713 煤层温度/K 293.15
    初始瓦斯扩散系数/
    (m2·s−1
    5.48×10−12 煤体骨架
    弹性模量/MPa
    5 139
    吸附常数b/MPa−1 2 吸附常数a/(m3·t−1 35
    下载: 导出CSV
  • [1] 刘程,孙东玲,武文宾,等. 我国煤矿瓦斯灾害超前大区域精准防控技术体系及展望[J]. 煤田地质与勘探,2022,50(8):82−92.

    LIU Cheng, SUN Dongling, WU Wenbin, et al. Analysis and prospect of precise prevention and control technical system ahead of large areas and time for gas disasters in China[J]. Coal Geology & Exploration, 2022, 50(8): 82−92.

    [2] 赵佳佳,田世祥,杨家向,等. 基于知识图谱的煤矿热动力灾害可视化研究[J]. 煤矿安全,2023,54(3):33−39.

    ZHAO Jiajia, TIAN Shixiang, YANG Jiaxiang, et al. Study on visualization of thermodynamic disaster in coal mine based on knowledge map[J]. Safety in Coal Mines, 2023, 54(3): 33−39.

    [3] 曹垚林. 水力化技术防治煤与瓦斯突出研究现状及展望[J]. 煤矿安全,2020,51(10):60−66.

    CAO Yaolin. Research status and prospect of coal and gas outburst prevention and control by hydraulic technology[J]. Safety in Coal Mines, 2020, 51(10): 60−66.

    [4] 侯国培,郭昆明,岳茂庄,等. 高位定向长钻孔瓦斯抽采技术应用[J]. 煤炭工程,2019,51(1):64−67.

    HOU Guopei, GUO Kunming, YUE Maozhuang, et al. Application of gas drainage technology in high-level directional long drilling hole[J]. Coal Engineering, 2019, 51(1): 64−67.

    [5] 郝世俊,段会军,莫海涛,等. 大直径高位定向长钻孔瓦斯抽采技术及实践[J]. 煤田地质与勘探,2020,48(6):243−248. doi: 10.3969/j.issn.1001-1986.2020.06.032

    HAO Shijun, DUAN Huijun, MO Haitao, et al. Gas drainage technology and practice analysis of large diameter high position directional long borehole[J]. Coal Geology & Exploration, 2020, 48(6): 243−248. doi: 10.3969/j.issn.1001-1986.2020.06.032

    [6] 乔伟,王凯,程波. 顺层长钻孔瓦斯抽采工艺技术应用研究[J]. 矿业安全与环保,2021,48(6):93−98.

    QIAO Wei, WANG Kai, CHENG Bo. Study on the application of gas drainage technology in bedding long drilling[J]. Mining Safety & Environmental Protection, 2021, 48(6): 93−98.

    [7] 王兆丰,田富超,赵彬,等. 羽状千米长钻孔抽采效果考察试验[J]. 煤炭学报,2010,35(1):76−79.

    WANG Zhaofeng, TIAN Fuchao, ZHAO Bin, et al. The testing of gas drainage efficiency about feather-veined 1 000 m length-drillhole[J]. Journal of China Coal Society, 2010, 35(1): 76−79.

    [8] 张飞燕,韩颖,吕帅,等. 基于残余瓦斯含量的顺层钻孔抽采有效半径阶梯式测定法[J]. 河南理工大学学报(自然科学版),2021,40(6):10−16.

    ZHANG Feiyan, HAN Ying, LYU Shuai, et al. Staged testing method for gas drainage effective radius of hole drilled along seam based on residual gas content[J]. Journal of Henan Polytechnic University (Natural Science), 2021, 40(6): 10−16.

    [9] 张建国,李喜员,高建成,等. 长钻孔远距离定点瓦斯含量测定研究[J]. 工矿自动化,2021,47(8):33−40.

    ZHANG Jianguo, LI Xiyuan, GAO Jiancheng, et al. Study on long drilling long-distance fixed-point gas content measurement[J]. Industry and Mine Automation, 2021, 47(8): 33−40.

    [10] 江成,李奇贤,韩恩德. 瓦斯抽采过程中本煤层及邻近层储层参数演化特征[J]. 煤矿安全,2021,52(11):1−7

    JIANG Cheng, LI Qixian, HAN Ende. Evolution characteristics of reservoir parameters of drainage layer and adjacent layer in the process of gas drainage[J]. Safety in Coal Mines, 2021, 52(11): 1−7.

    [11] 李路广,李向阳,魏路浩,等. 千米定向钻机在大宁煤矿瓦斯抽采中的应用[J]. 煤炭工程,2021,53(10):84−88.

    LI Luguang, LI Xiangyang, WEI Luhao, et al. Application of directional kilometer drilling machine in gas extraction in Daning Coal Mine[J]. Coal Engineering, 2021, 53(10): 84−88.

    [12] 冀超辉,崔洪庆. 特厚急倾斜煤层瓦斯抽采关键技术与实践[J]. 矿业安全与环保,2021,48(6):99−103.

    JI Chaohui, CUI Hongqing. Key technology and practice of gas extraction in steep extra-thick coal seam[J]. Mining Safety & Environmental Protection, 2021, 48(6): 99−103.

    [13] 翟成,林柏泉,吴海进. 顶板高位钻孔抽放在瓦斯治理中的应用[J]. 煤炭工程,2005,37(9):4−6.

    ZHAI Cheng, LIN Boquan, WU Haijin. Gas drainage with high level borehole in roof applied to gas control[J]. Coal Engineering, 2005, 37(9): 4−6.

    [14] 毕慧杰,邓志刚,赵善坤,等. 高瓦斯综采工作面定向高位钻孔瓦斯抽采技术研究[J]. 煤炭科学技术,2019,47(4):134−140.

    BI Huijie, DENG Zhigang, ZHAO Shankun, et al. Research on gas drainage technology of directional high-level borehole in high gassy fully-mechanized mining face[J]. Coal Science and Technology, 2019, 47(4): 134−140.

    [15] 李宏,马金魁. 大直径顶板定向长钻孔替代高抽岩巷的瓦斯抽采效果分析[J]. 煤炭科学技术,2020,48(7):304−310.

    LI Hong, MA Jinkui. Analysis on gas drainage effect of large diameter roof directional long borehole instead of high drainage rock roadway[J]. Coal Science and Technology, 2020, 48(7): 304−310.

    [16] 王广宏. 基于定向长钻孔抽采瓦斯的数值模拟研究[J]. 矿业安全与环保,2021,48(4):38−42.

    WANG Guanghong. Numerical simulation study of gas drainage based on directional long drilling[J]. Mining Safety & Environmental Protection, 2021, 48(4): 38−42.

    [17] 徐超,王建文,杜昌昂,等. 瓦斯抽采长钻孔负压沿孔长衰减机制及影响因素模拟研究[J]. 采矿与安全工程学报,2021,38(2):419−428.

    XU Chao, WANG Jianwen, DU Chang’ang, et al. Simulation of attenuation mechanism of negative pressure along the hole length of long boreholes in gas drainage and its influencing factors[J]. Journal of Mining & Safety Engineering, 2021, 38(2): 419−428.

    [18] 王兆丰,刘慧鹏,杨宏民,等. 伏岩煤业抽采负压沿钻孔深度变化规律[J]. 煤矿安全,2013,44(10):156−157.

    WANG Zhaofeng, LIU Huipeng, YANG Hongmin, et al. Changing laws of drainage negative pressure along with the change of drilling depth in Fuyan coal industry[J]. Safety in Coal Mines, 2013, 44(10): 156−157.

    [19] 李杰. 预抽瓦斯钻孔抽采效果沿孔长变化规律研究[D]. 焦作:河南理工大学,2012.
    [20] 宋建军. 瓦斯抽采钻孔负压沿孔长动态变化特性研究[J]. 煤炭技术,2016,35(10):212−214.

    SONG Jianjun. Study on dynamic changes of negative pressure along gas drainage borehole[J]. Coal Technology, 2016, 35(10): 212−214.

    [21] 周军平,鲜学福,姜永东,等. 基于热力学方法的煤岩吸附变形模型[J]. 煤炭学报,2011,36(3):468−472.

    ZHOU Junping, XIAN Xuefu, JIANG Yongdong, et al. A model of adsorption induced coal deformation based on thermodynamics approach[J]. Journal of China Coal Society, 2011, 36(3): 468−472.

    [22] 石浩. 大直径高位定向长钻孔瓦斯抽采技术及应用[J]. 煤炭科学技术,2018,46(10):190−195.

    SHI Hao. Gas drainage technology and its application of large diameter and high directional long drilling[J]. Coal Science and Technology, 2018, 46(10): 190−195.

    [23]

    WANG K, ZANG J, WANG G, et al. Anisotropic permeability evolution of coal with effective stress variation and gas sorption: Model development and analysis[J]. International Journal of Coal Geology, 2014, 130: 53−65.

    [24] 撒占友,吴静波,杨永亮,等. 考虑Klinkenberg因子状态的CO2-ECBM模拟研究[J]. 煤田地质与勘探,2023,51(3):37−45. doi: 10.12363/issn.1001-1986.22.09.0713

    SA Zhanyou, WU Jingbo, YANG Yongliang, et al. CO2-ECBM simulation study considering Klinkenberg factor state[J]. Coal Geology & Exploration, 2023, 51(3): 37−45. doi: 10.12363/issn.1001-1986.22.09.0713

  • 期刊类型引用(1)

    1. 董志伟,杨世瑜,白帆,程晓强. 综采工作面瓦斯分源抽采技术研究与应用. 当代化工研究. 2025(03): 121-123 . 百度学术

    其他类型引用(0)

图(7)  /  表(1)
计量
  • 文章访问数:  35
  • HTML全文浏览量:  4
  • PDF下载量:  12
  • 被引次数: 1
出版历程
  • 收稿日期:  2023-02-16
  • 修回日期:  2023-06-25
  • 刊出日期:  2024-11-19

目录

/

返回文章
返回