A calculation method of reasonable size of coal pillar in large mining height section based on elastic theory
-
摘要:
不同工况下区段煤柱两侧支承压力分布及岩体变形存在显著差异，考虑煤柱两侧不同支承压力对煤柱整体稳定性的影响，基于大采高区段煤柱的弹性力学计算模型，分析了支承压力下煤柱任一单元岩体的应力应变分量。通过建立大采高煤柱弹塑性界面上岩体的柱条模型，确定在0.65倍煤柱高度处单元岩体将首先发生水平拉伸破坏，利用虎克定律提出了该单元岩体极限拉应变与煤柱极限平衡区宽度的关系式。依据煤柱破裂区岩体的受力特征，运用摩尔库伦准则推导了煤柱破裂区宽度的计算公式。结果表明：①煤柱极限平衡区宽度与岩体极限拉应变和弹性模量反相关，与煤柱埋深和煤柱高度正相关；②煤柱高度及其与顶底板的界面摩擦角是影响破裂区宽度的关键性因素；③煤柱两侧不同工况下，煤柱岩体极限拉应变与其所受侧压呈正变关系，区段煤柱采空区侧所受侧压较巷道侧偏大，采空区侧岩体的极限拉应变也相应较大，表现为采空侧极限平衡区宽度较巷道侧偏小。最后，将上述理论公式应用于陕北某矿30109工作面大采高区段煤柱极限平衡区和破裂区宽度的分析计算，给出了该工作面两侧区段煤柱的合理宽度及其支护方案。工程应用表明，30109工作面区段巷道围岩变形控制效果良好，满足现场生产需求。
Abstract:There are significant differences in abutment pressure distribution and rock mass deformation on both sides of section coal pillar under different working conditions, the influence of different abutment pressures on the overall stability of coal pillar is considered, based on the elastic mechanics calculation model of coal pillar in large mining height section, the stress-strain components of any unit rock mass of coal pillar under abutment pressure are analyzed. Through the establishment of large mining height pillar elastic-plastic interface on the rock column model, it is determined that unit rock mass will first undergo horizontal tensile failure at 0.65 times the height of coal pillar, based on Hooke's law, the relationship between the ultimate stretching strain of rock mass and the width of limit equilibrium zone of coal pillar is proposed. According to the mechanical characteristics of the rock mass in the coal pillar fracture zone, the calculation formula of the width of the coal pillar fracture zone is derived by using the Mohr-Coulomb criterion. The results show that: ①The width of limit equilibrium zone of coal pillar is inversely related to ultimate tensile strain and elastic modulus of rock mass and positively related to buried depth and height of coal pillar;②The height of coal pillar and the interface friction angle between coal pillar and roof and floor are the key factors affecting the width of fracture zone; ③Under different working conditions on both sides of coal pillar, due to the positive relationship between the ultimate tensile strain of coal pillar rock mass and its lateral pressure, the lateral pressure on the gob side of the section coal pillar is larger than that on the roadway side, and the ultimate tensile strain of the rock mass on the gob side is correspondingly larger, which shows that the width of the limit equilibrium area on the gob side is smaller than that on the side of the gateway. Finally, the above theoretical formula is applied to the analysis and calculation of the width of limit equilibrium zone and fracture zone of coal pillar in large mining height section of 30109 working face in a mine in Northern Shaanxi, and the reasonable width of coal pillar and its supporting scheme under different working conditions on both sides of the working face are given. The engineering application shows that the deformation control effect of surrounding rock along the working face is good, which can meet the demand of field production.
-
0. 引 言
近年来，煤柱型冲击地压灾害事故频率不断增加，严重威胁矿井的生产安全[1-2]。例如2016年8月15日0时33分，山东能源肥矿集团梁宝寺能源有限责任公司发生了煤柱型冲击地压事故，造成2人死亡。2019年8月2日12时24分，开滦(集团)有限责任公司唐山矿业分公司也发生了煤柱型冲击地压事故，造成了12人伤亡，直接经济损失高达614.024万元。实践表明，煤柱型冲击地压一般具有隐蔽性强、无明显前兆特征等特点，造成设备破坏、巷道变形和人员伤亡等事故，给矿井的安全生产带来严重威胁。因此，有必要开展煤柱动力失稳机理研究，以期对煤柱型冲击地压预警与防治提供理论依据。
国内学者针对煤柱稳定性问题进行了大量的研究。部分学者从静力学角度进行了大量的试验。殷志强等[3]通过分离式霍普金森压杆试验，分析了不同静载下煤样的动态力学特性，得出峰值应力随静载增大而增大，峰值应变随静载增大而减小的结论。何永琛[4]通过对煤岩组合体试件进行单轴压缩试验，获得了不同煤岩组合体试件的破坏形态及峰值应变等变化特征。艾迪昊等[5]通过对不同粒度的型煤进行单轴压缩试验，得出型煤的破坏可分为压实，初始破裂，加速破裂和卸压4个阶段。李成武等[6]通过室内试验和现场试验研究了煤体静载破坏中低频磁场时、频谱特征，并结合微震信号提出了低频磁信号的产生机制。
也有学者从动力学的角度，探讨了煤柱稳定性问题。徐青云等[7]采用理论分析、数值模拟和现场实测等方法，指出窄煤柱在受到剧烈采动影响时，主要在其中上部区域发生变形破坏。姜学伟[8]通过分析房采区煤柱在采动影响下围岩受力、煤柱破坏变形状态等得出了煤柱保持稳定性所必须的弹性区宽度计算公式，并作为判断煤柱是否破坏失稳的依据。李峰等[9]基于LS-DYNA有限元分析软件的HJC本构模型，模拟不同动载强度下煤柱的动态损伤过程。结果表明，应变率相同时，煤柱受动载荷越大，破坏越严重。王方田等[10]将动载荷以应力波的形式作用在煤体上，来分析不同动载荷烈度下煤柱坝体的动力响应机制，发现巷道围岩塑性区发育面积及深度显著增加，竖向应力积聚强度增大，并与动载荷烈度成正比。丛利等[11]通过数值模拟对相邻工作面回采中和回采后3次高强度开采扰动对煤柱内应力变化情况的研究，得出煤柱内布置的联络巷会加剧巷道围岩的高应力集中，在外部动载扰动下，极易达到冲击载荷并诱发冲击的结论。因煤柱在外部动载扰动下的破坏失稳机理极其复杂，当前的理论研究和实践仍需进一步完善。因此针对不同尺寸深井护巷煤柱在外部动载扰动下的破坏特征和突变失稳机制开展了详细研究。
首先，采用数值模拟分析了动载冲击下煤柱失稳特征；其次，建立了动载作用下高静载煤柱突变失稳力学模型，推导了煤柱稳定性判别准则，揭示了动载作用下深井煤柱突变失稳机理，以期对煤柱稳定性控制提供参考。
1. 深井护巷煤柱数值模型的建立
1.1 模型建立与参数选择
利用FLAC3D软件建立“顶板-煤柱-底板”数值计算模型，模型宽高比分别设置为4、6、8和10，煤层高度保持为3 m。图1为宽高比为10的数值计算模型，模型共计9层，7种岩性，模型尺寸为34 m×34 m×39 m，共划分62820个单元网格。参考相关数值计算经验[12-13]，顶底板岩层采用各向同性弹性模型，煤层采用应变软化模型。煤层节理、裂隙发育，室内试验获得的煤样物理力学参数很难直接用于数值计算。因此，采用试错法拟合煤柱强度公式以校核煤柱力学参数。相关学者提出了多种煤柱强度公式[14]，其中Salamon公式能够准确描述宽高比为2～20时的煤柱强度，计算结果精确度能够满足现场要求，可用于计算现场保护煤柱的强度[15]。因此，采用Salamon公式校正数值模型计算结果。
$$ {S_{\rm{p}}} = 7.716\frac{{{W^{0.46}}}}{{{H^{0.66}}}} $$ (1) 式中：Sp为保护煤柱强度， MPa；W为保护煤柱宽度， m；H为保护煤柱高度， m。
采用试错法时，首先建立“顶板-煤柱-底板”数值计算模型；其次，在模型顶部施加速度边界条件以模拟单轴压缩，获得煤柱单轴压缩强度；再者，通过调整煤层各项力学参数，使数值计算获得煤柱单轴压缩强度与Salamon公式匹配。图2为数值计算获得的煤柱强度与Salamon公式对比图，由图2可知，数值计算获得的煤柱强度与Salamon公式基本保持一致，因此，采用表1所示的岩性参数模拟煤柱的应变软化行为。
![]()
图 2 数值模拟和Salamon公式计算峰值强度对比Figure 2. Comparison of peak strength between numerical simulation and Salamon formula calculation表 1 煤柱模拟参数Table 1. Simulation parameters of protective coal pillar岩性 粗砂岩 细砂岩 细粉砂岩互层 泥岩 3煤 粉砂岩 体积模量/GPa 13.56 8.5 14 3 2.59 16 剪切模量/GPa 8.2 5.6 8.7 1.3 1.22 10 抗拉强度/MPa 5.7 4.7 3.4 1.5 5.1 2.3 黏聚力/MPa 8.9 5.52 4.3 1.7 4.2 2.4 内摩擦角/(°) 31 35 34 29 34 32 1.2 边界条件及模型方案
以新河煤矿6302工作面南部轨道集中巷，6301运输巷和6302轨道巷之间的保护煤柱为研究对象。煤层埋深980 m，煤的容重为24 kN/m³，故上覆岩层的自重应力为24 MPa。
数值计算分为2步。①在模型顶板施加恒定载荷，开展静载计算；②待静载计算平衡后，在模型顶板施加速度边界条件以模拟动载。
静载计算时，模型顶部设置为应力边界，上覆岩层自重采用均布荷载代替，底部为固定边界，顶部垂直应力取上覆岩层的自重应力(24 MPa)。在静载条件下，模型达到平衡状态，再在动载条件下开展计算。动力学计算时，数值模型顶部边界自由，底部设置为固定边界。
为了模拟动载扰动，假设上覆顶板岩层断裂产生含有多种子波的压缩波，而任何复杂的波形都可以由正弦波、余弦波或三角波合成。为简化动载扰动，选择正弦应力波作为输入波形。DUBINSKI[16]等现场监测发现应力波PPV为0.03 m/s。而HADJIGEORGIOU等[17]通过使用炸药诱发岩体振动所记录的应力波 PPV不能大于3.3 m/s，因此，为了模拟不同震动强度的应力波，振幅设定为0.6、0.8、1.0、1.2、1.4 m/s。
2. 动载扰动下煤柱失稳破坏变化规律
选取宽高比分别为4、6、8、10四种模型，静载应力下平衡后，在煤柱模型的顶部分别施加振幅为0.6、0.8、1.0、1.2、1.4 m/s的应力波（以应力波振幅表示动载强度），当模型再次达到平衡后，研究应力波输入保护煤柱破坏失稳规律。
2.1 煤柱尺寸对煤柱动力失稳的影响
由图3可知，动载冲击后煤柱塑性区由表及里向深部扩展最终导致煤柱整体失稳。相比而言，随着煤柱尺寸增加，煤柱发生整体失稳所需的动载冲击强度具有一定的差异。以应力波振幅为0.6 m/s的模拟结果为例，以阐述此种差异。
![]()
图 3 不同应力波振幅及不同尺寸保护煤柱塑性区分布规律Figure 3. Distribution law of plastic zone of protective coal pillars with different sizes and different stress wave amplitudes煤柱宽高比为6时，静载作用后，塑性区宽度为7～8 m，内部存在弹性核心。受到振幅为0.6 m/s应力波冲击后，煤柱弹性核心完全破坏，塑性区贯穿整个煤柱，意味着煤柱整体失稳。相比而言，煤柱宽高比增加至8与10时，煤柱塑性区仅向深度扩展，内部依旧存在稳定的弹性核心，且弹性核心所占整个煤柱的比例由21%增加至47%。此外，对比发现，在相同静载作用下，不同尺寸的煤柱失稳启动所需的动载强度亦不同。模拟发现，随着煤柱宽高比由6增加至10时，煤柱破坏失稳所需的动载冲击强度（以应力波振幅表征）由0.6 m/s增加至1.2 m/s，提高了1倍。由此说明，增加煤柱尺寸有助于提高煤柱的稳定性，减少动力失稳的几率。
2.2 动载冲击强度对煤柱动力失稳的影响
以煤柱宽高比为10为例，分析动载冲击强度为0.6、0.8、1.0、1.2、1.4 m/s时，煤柱失稳破坏特征。如图3c所示，当应力波振幅从0.6 m/s增加至1.0 m/s，煤柱塑性区由浅部逐步向着中部扩展，范围由8 m增加至13.5 m，煤柱内部弹性核区逐步减少。当动载冲击强度增加至1.2 m/s时，煤柱整体产生塑性破坏，意味着煤柱整体失稳。由此说明，动载冲击强度越大，煤柱发生失稳破坏的概率也相对提高。
通过对比不同振幅应力波0.6、0.8、1.0、1.2、1.4 m/s作用下，煤柱的失稳破坏特征，发现煤柱的稳定性与煤柱尺寸和动载冲击强度有关。相同条件下，尺寸越大，煤柱发生失稳破坏的几率越小；动载强度越大，煤柱破坏范围及程度越大，即保护煤柱失稳概率越大。
3. 动载作用下煤柱突变失稳机制
3.1 煤柱承载力学模型
一般而言，煤柱载荷主要来源于上覆岩层重量和动载扰动。上覆岩层重量作为静载，使煤柱产生塑性区和弹性核区[18]，而后续动载冲击使塑性区向煤柱深处扩展，导致煤柱发生整体性破坏，内部弹性区减小而发生失稳，煤柱承载力学模型如图4所示。
煤柱埋深为H，上覆岩层容重为γ，煤柱高度为Hm，煤柱宽度为B。假设外部动载q0均匀施加至煤柱，动荷载在某一时刻的运动状态可简化为静荷载[19]，动载简化系数可以表示为kd，则煤柱承受的载荷为：
$$ p = \left( {\gamma H + {k_{\rm{d}}}{q_0}} \right)B $$ (2) 煤柱塑性区与弹性区本构关系有所不同，弹性区煤体呈线弹性关系，塑性区内呈非线性关系，具有应变软化特征，煤柱塑性区内本构关系可以表示为
$$ \sigma {\text{ = }}E\varepsilon \exp \left[ { - \frac{1}{2}{{\left( {\frac{\varepsilon }{{{\varepsilon _0}}}} \right)}^2}} \right] $$ (3) 式中：E为煤柱的初始弹性模量；ε为应变；ε0为极限应变。
假设煤柱塑性区宽度为B0，则煤柱塑性区内载荷ps与煤柱压缩量u的关系可以表示为
$$ {p_{\rm{s}}} = \frac{{{B_0}Eu}}{{{H_{\rm{m}}}}}{{\rm{e}}^{\left[ { - \tfrac{1}{2}{{\left( {\tfrac{u}{{{u_0}}}} \right)}^2}} \right]}} $$ (4) 式中：u0为煤柱最大压缩量， m。
弹性区内，煤体载荷pt与煤柱压缩量u的关系为
$$ {p_{\rm{t}}} = \frac{{\left( {B - {B_0}} \right)Eu}}{{{H_{\rm{m}}}}} $$ (5) 3.2 煤柱失稳的尖点突变模型
整个煤柱系统的总势能由煤柱自身的应变能和外部载荷所做的功组成。由式（4）可得煤柱塑性区应变能
${V_1}\left( u \right) $ 为$$ {V_1}\left( u \right) = \frac{{{B_0}E}}{{{H_{\rm{m}}}}}\int_0^u {u{{\rm{e}}^{\left[ { - \tfrac{1}{2}{{\left( {\tfrac{u}{{{u_0}}}} \right)}^2}} \right]}}} {\rm{d}}u $$ (6) 由式（5）可得煤柱弹性区应变能
${V_2}\left( u \right) $ 为$$ {V_2}\left( u \right) = \frac{{\left( {B - {B_0}} \right)E}}{{{H_{\rm{m}}}}}\int_0^u {u{\rm{d}}u} $$ (7) 上覆岩层静载重力势能和外部动载荷做功的势能函数
${V_3}\left( u \right) $ 为$$ {V_3}\left( u \right) = pu = \left( {\gamma H + {k_{\rm{d}}}{q_0}} \right)Bu $$ (8) 根据上述力学模型的分析，动载作用下煤柱系统的总势能函数
${V}\left( u \right) $ 可以表示为$$ \begin{gathered} V\left( u \right){\text{ = }} = \frac{{{B_0}E}}{{{H_{\rm{m}}}}}\int_0^u {u{{\rm{e}}^{\left[ { - \tfrac{1}{2}{{\left( {\tfrac{u}{{{u_0}}}} \right)}^2}} \right]}}} {\rm{d}}u{\text{ + }} \\ \frac{{\left( {B - {B_0}} \right)E}}{{{H_{\rm{m}}}}}\int_0^u {u{\rm{d}}u} - \left( {\gamma H + {k_{\rm{d}}}{q_0}} \right)Bu \\ \end{gathered} $$ (9) 以u为状态变量进行突变理论分析，对V(u)求导得：
$$ u = {u_1} = \sqrt 3 {u_0} $$ (10) 引入无量纲常量x，令
$ x=u-u_{1} $ ，求解得：$$\begin{split} & {x^3} + \left[ {\frac{{{{\rm{e}}^{\tfrac{3}{2}}}\left( {B - {B_0}} \right){u_0}^2}}{{4{B_0}}} - \frac{{{u_0}^2}}{2}} \right]x + \frac{{\sqrt 3 u_0^3}}{4}\left[ {1 + \frac{{{{\rm{e}}^{\tfrac{3}{2}}}\left( {B - {B_0}} \right)}}{{{B_0}}}} \right] - \\&\qquad\qquad \frac{{B{H_{\rm{m}}}{u_0}^2\left( {\gamma H + {k_{\rm{d}}}{q_0}} \right){{\rm{e}}^{\tfrac{3}{2}}}}}{{4E{B_0}}} = 0 \end{split} $$ (11) 进一步简化，令：
$$ \left\{ \begin{gathered} m = \left[ {\frac{{{{\rm{e}}^{\tfrac{3}{2}}}\left( {B - {B_0}} \right){u_0}^2}}{{4{B_0}}} - \frac{{{u_0}^2}}{2}} \right] \\ n = \frac{{\sqrt 3 u_0^3}}{4}\left[ {1 + \frac{{{{\rm{e}}^{\tfrac{3}{2}}}\left( {B - {B_0}} \right)}}{{{B_0}}}} \right] - \\\qquad \frac{{B{H_{\rm{m}}}{u_0}^2\left( {\gamma H + {k_{\rm{d}}}{q_0}} \right){{\rm{e}}^{\tfrac{3}{2}}}}}{{4E{B_0}}} \\ \end{gathered} \right. $$ (12) 将式（12）代入式（11）则有：
$$ {x^3} + mx + n = 0 $$ (13) 以x为状态参量，m、n为控制参量。由尖点突变理论可知，当m、n满足分叉集方程时，顶板-煤柱系统处于临界平衡状态。联立上式与其拐点方程
$ 3{x^2} + m = 0 $ ，消去状态参量x，得到分叉集方程：$$ \varDelta = 4{m^3} + 27{n^2} = 0 $$ (14) 由突变理论可知，当Δ＞0时，系统处于稳定状态，煤柱不会失稳；当Δ=0时，系统处于临界状态，受到微小扰动会发生失稳；当Δ＜0时，是系统发生失稳的充要条件。将式（12）代入式（14）得：
$$ \begin{gathered} \qquad\quad \varDelta = 4{\left[ {\frac{{{{\rm{e}}^{\tfrac{3}{2}}}\left( {B - {B_0}} \right){u_0}^2}}{{4{B_0}}} - \frac{{{u_0}^2}}{2}} \right]^3} + \\ 27{\left\{ \begin{gathered} \frac{{\sqrt 3 u_0^3}}{4}\left[ {1 + \frac{{{{\rm{e}}^{\tfrac{3}{2}}}\left( {B - {B_0}} \right)}}{{{B_0}}}} \right] - \frac{{B{H_{\rm{m}}}{u_0}^2\left( {\gamma H + {k_{\rm{d}}}{q_0}} \right){{\rm{e}}^{\tfrac{3}{2}}}}}{{4E{B_0}}} \\ \end{gathered} \right\}^2} \\ \end{gathered} $$ (15) 煤柱失稳则方程（13）有解，必须满足Δ≤0，由此得到煤柱失稳的必要条件为：
$$ \begin{gathered} \qquad\quad \varDelta = 4{\left[ {\frac{{{{\rm{e}}^{\tfrac{3}{2}}}\left( {B - {B_0}} \right){u_0}^2}}{{4{B_0}}} - \frac{{{u_0}^2}}{2}} \right]^3} + \\ 27{\left[ \begin{gathered} \frac{{\sqrt 3 u_0^3}}{4}\left( {1 + \frac{{{{\rm{e}}^{\tfrac{3}{2}}}\left( {B - {B_0}} \right)}}{{{B_0}}}} \right) - \frac{{B{H_{\rm{m}}}{u_0}^2\left( {\gamma H + {k_{\rm{d}}}{q_0}} \right){{\rm{e}}^{\tfrac{3}{2}}}}}{{4E{B_0}}} \\ \end{gathered} \right]^2} \leqslant 0 \\ \end{gathered} $$ (16) 3.3 煤柱最大压缩量计算
煤柱失稳判别式Δ中还含有一个与煤柱相关的未知量u0（煤柱最大压缩量），下面对煤柱最大压缩量u0进行求解。
煤柱的最大压缩量是由顶板的最大下沉量决定，如图5a和5b所示。考虑到煤柱的尺寸较大，拟采用弹性薄板理论研究顶板的最大下沉量[20-21]。图5c所示为直接顶-煤柱弹性薄板模型，其中，2a为煤柱长度，2b为煤柱宽度，q为直接顶上覆动静载等效均布载荷。直接顶可以视为四边固支的矩形板，其力学模型如图5d所示。
该模型满足以下边界条件：
$$\left \{\begin{array}{l}W|{}_{x = -a} = W|{}_{x = a} = 0\text{，}\dfrac{\partial W}{\partial x}|{}_{x = -a} = \dfrac{\partial W}{\partial x}|{}_{x = a} = 0\\ W|{}_{y = -a} = W|{}_{y = b} = 0\text{，}\dfrac{\partial W}{\partial y}|{}_{y = -b} = \dfrac{\partial W}{\partial y}|{}_{y = b} = 0\end{array}\right. $$ (17) 根据弹性力学中的伽辽金解法，将板的挠度展开为以下表达式：
$$ W = \displaystyle\sum\limits_{i = 1}^\infty {{A_{{i}}} {f_{{i}}}} \left( {x,y} \right) $$ (18) 式中：W为直接顶挠度；Ai为挠度系数；fi(x,y)为挠度函数，要求能满足薄板挠曲面的形状和几何边界条件，而且是连续的。
为便于计算，取一级近似函数：
$$ W = {A_1}f\left( {x,y} \right){\text{ = }}{A_1}{\left( {{x^2} - {a^2}} \right)^2}{\left( {{y^2} - {b^2}} \right)^2} $$ (19) 显然该式能满足薄板挠曲面的形状和所有几何边界条件。根据伽辽金法对上式进行求解，得到挠度的表达式：
$$\begin{gathered} W = \frac{{7q}}{{128D\left( {{a^4} + {b^4} + \frac{4}{7}{a^2}{b^2}} \right)}}{\left( {{x^2} - {a^2}} \right)^2}{\left( {{y^2} - {b^2}} \right)^2} \\ \qquad\qquad\qquad D=\dfrac{E h^{3}}{12\left(1-v^{2}\right)} \end{gathered} $$ (20) 式中，D为板的抗弯刚度，
$D=\dfrac{E h^{3}}{12\left(1-v^{2}\right)} $ ，$v$ 为泊松比；h为直接顶高度。对挠度方程求极值，在板的中心处取到挠度的最大值：
$$ {W}_{\mathrm{max}}=W\left(0\text{，}0\right)\text=\frac{7q{a}^{4}{b}^{4}}{128D\left({a}^{4}+{b}^{4}+\dfrac{4}{7}{a}^{2}{b}^{2}\right)} $$ (21) 煤柱最大压缩量等于直接顶最大挠度，如图5e所示，即
$ u_{0}=W_{\text {max }} $ ，至此得到了煤柱失稳的最终判别式：$$ \begin{gathered} \qquad \varDelta = 4{\left[ {\frac{{{{\rm{e}}^{\tfrac{3}{2}}}\left( {B - {B_0}} \right){W_{\max }}^2}}{{4{B_0}}} - \frac{{{W_{\max }}^2}}{2}} \right]^3} +27 \times \\ {\left[ \frac{{\sqrt 3 {W_{\max }}^3}}{4}\left( {1 + \frac{{{{\rm{e}}^{\tfrac{3}{2}}}\left( {B - {B_0}} \right)}}{{{B_0}}}} \right) - \frac{{B{H_{\rm{m}}}{W_{\max }}^2\left( {\gamma H + {k_{\rm{d}}}{q_0}} \right){{\rm{e}}^{\tfrac{3}{2}}}}}{{4E{B_0}}} \right]^2} \\ \end{gathered} $$ (22) 4. 深井煤柱稳定性判别
4.1 深井煤柱稳定性判别方法
通过现场调研获得煤柱的物理力学参数，计算煤柱塑性区范围。
$$ {B_0}{\text{ = }}\frac{{{H_{\rm{m}}}A}}{{\tan\; {\varphi _0}}}\ln \left( {1 + \frac{{kH\tan \;{\varphi _0}}}{{{c_0}}}} \right) $$ (23) 式中：k为应力集中系数；c0为黏聚力；
$ {\varphi _0} $ 为内摩擦角；A为侧压系数。计算煤柱最大压缩量：
$$ {u_0} = {W_{\max }}{\text{ = }}\frac{{7q{a^4}{b^4}}}{{128D\left( {{a^4} + {b^4} + \dfrac{4}{7}{a^2}{b^2}} \right)}} $$ (24) 最后，利用式（21）计算Δ值，进行稳定性判别。
4.2 工程背景
新河煤矿位于山东省济宁市唐口勘探区西南部边缘，目前开采水平为−980 m延伸水平，延伸水平位于矿井的东南部，东西宽2 km，南北长2.4 km。新河煤矿正在开采6302工作面，集中大巷分布在6302工作面以南，大巷之间的煤柱对大巷的稳定起到支撑保护作用。利用前文建立的煤柱稳定性判别方法，对新河煤矿6302工作面南部轨道集中巷，6301运输巷和6302轨道巷之间的煤柱进行稳定性进行判别，如图6所示。
1) 计算煤柱塑性区范围：煤柱高2.9 m，黏聚力1.7 MPa，内摩擦角29°，埋深980 m，容重24 kN/m³，侧压系数1.12，应力集中系数1.97。将煤柱参数带入煤柱塑性区范围计算公式：
$$ {B_0}{\text{ = }}\frac{{{H_{\rm{m}}}A}}{{\tan\; {\varphi _0}}}\ln \left( {1 + \frac{{k\gamma H\tan \;{\varphi _0}}}{{{c_0}}}} \right) =16.4\;{\rm{m}} $$ 2) 计算煤柱最大压缩量：大巷保护煤柱长度46 m，宽度32 m，等效荷载3.75×107 N，泥岩弹性模量6.5 GPa，泊松比0.4，顶板厚度6.63 m。
$$ {u_0} = {W_{\max }}{\text{ = }}\frac{{7q{a^4}{b^4}}}{{128D\left( {{a^4} + {b^4} + \dfrac{4}{7}{a^2}{b^2}} \right)}} =0.473\;2 $$ 3) 利用尖点突变模型计算Δ值，并进行煤柱稳定性判别：
$$ \begin{gathered} \qquad\quad \varDelta = 4{\left[ {\frac{{{{\rm{e}}^{\tfrac{3}{2}}}\left( {B - {B_0}} \right){W_{\max }}^2}}{{4{B_0}}} - \frac{{{W_{\max }}^2}}{2}} \right]^3} +27 \times \\ {\left[ \begin{gathered} \frac{{\sqrt 3 {W_{\max }}^3}}{4}\left( {1 + \frac{{{{\rm{e}}^{\tfrac{3}{2}}}\left( {B - {B_0}} \right)}}{{{B_0}}}} \right) - \frac{{B{H_{\rm{m}}}{W_{\max }}^2\left( {\gamma H + {k_{\rm{d}}}{q_0}} \right){{\rm{e}}^{\tfrac{3}{2}}}}}{{4E{B_0}}} \\ \end{gathered} \right]^2} = \\ \qquad\quad\qquad\quad\qquad\quad\qquad\quad 1.58 \\ \end{gathered} $$ 根据煤柱稳定性判别公式可知，Δ＝1.58＞0，该煤柱能保持稳定。
5. 结 论
1) 数值模拟研究表明，煤柱受到外部动载扰动后，塑性区向中部扩展，弹性核心占比减小，导致煤柱破坏失稳。煤柱动力失稳与煤柱尺寸和外界动载冲击强度有关。增大煤柱的宽度能够提高煤柱的稳定性，减小煤柱破坏失稳的概率。相同条件下，煤柱破坏失稳的几率随着动载强度的增大而增大。
2）基于弹性薄板理论，建立了直接顶-煤柱弹性薄板模型，推导了直接顶的挠度方程，获得了煤柱的最大压缩量表达式。基于尖点突变理论，建立了煤柱尖点突变模型，推导了主要受塑性区宽度和煤柱压缩量影响的煤柱突变失稳判别式，为煤柱突变失稳提供了力学判据。
3) 基于煤柱突变失稳判别式，对新河煤矿6302工作面南部轨道集中巷，6301运输巷和6302轨道巷之间的保护煤柱进行了验证。结果表明，新河煤矿6302工作面南部集中大巷保护煤柱能够保持稳定，与工程现场结果一致。
-
![]()
图 1 不同工况下煤柱支承压力分布模型
Figure 1. Distribution model of coal pillar supporting pressure under different working conditions
![]()
图 3 两侧不同工况下煤柱受力变形简化模型
Figure 3. Simplified model of mechanical deformation of coal pillar under different conditions on both sides
![]()
图 4 左侧极限平衡区线性荷载计算模型
Figure 4. Linear load calculation model for limit equilibrium zone
![]()
图 6 弹性区均布荷载计算模型
Figure 6. Calculation model of uniformly distributed load in elastic zone
![]()
图 9 区段煤柱岩体压缩柱条分析模型
Figure 9. Compression column bar analysis model of section coal pillar rock mass
![]()
图 10 煤柱极限平衡区柱条分析模型
Figure 10. Column bar analysis model for limit equilibrium zone of coal pillars
![]()
图 11 煤柱弹塑性界面岩体柱条模型
Figure 11. Fracture-flexural mechanical model of rock mass with elastic-plastic interface of coal pillar
![]()
图 12 煤柱极限平衡区宽度计算模型
Figure 12. Calculation model of the width of coal pillar limit equilibrium zone
![]()
图 13 煤柱破裂区简化计算模型
Figure 13. Simplified calculation model for coal pillar fracture zone
![]()
图 14 煤柱极限平衡区与煤体弹性模量的关系
Figure 14. Relationship between coal pillar limit equilibrium zone and coal elastic modulus
![]()
图 15 煤柱极限平衡区宽度与煤体极限拉应变的关系
Figure 15. Relationship between the ultimate equilibrium zone of coal pillar and ultimate tensile strain of coal
![]()
图 16 煤柱极限平衡区宽度与煤柱埋深的关系
Figure 16. Relationship between coal pillar limit equilibrium zone and buried depth of roadway
![]()
图 17 煤柱极限平衡区宽度与煤柱高度的关系
Figure 17. Relationship between coal pillar limit equilibrium zone and coal pillar height
![]()
图 20 工作面巷道围岩变形监测
Figure 20. Deformation monitoring results of surrounding rock along the channel of working face
表 1 煤帮各计算参数取值
Table 1 Value of calculating parameters of coal wall
$ \gamma $/（kN·m−3） K/GPa $ \lambda $ $ \varphi $/（º） $ c $/MPa $ {\varphi _{\rm{u}}} $/（º） $ {\varphi _{\rm{d}}} $/（º） $ E $/GPa $ \mu $ ${\varepsilon _{ {\rm{t} ,}\max } }$/10−3 D 25 0.15 0.8 30 1.8 25 25 3.6 0.3 0.68 1.6 注：D为直接顶与煤层弹量比。 
下载: 导出CSV
表 2 煤帮各计算参数取值
Table 2 Value of calculating parameters of coal wall
$ \gamma $/（kN·m−3） K/GPa $ \lambda $ $ \varphi $/(º) $ c $/MPa $ {\varphi _{\rm{u}}} $/（º） $ {\varphi _{\rm{d}}} $/（º） $ E $/GPa $ \mu $ $ {\varepsilon _{{\rm{t}}\;\max }}_{,{\rm{BC}}} $/×10−3 ${\varepsilon _{ {\rm{t} },\max } }_{,{\rm{EH} } }$/10−3 D 25 0.18 0.8 32 2.0 27 27 3.4 0.3 0.63 0.68 1.6
下载: 导出CSV
-
[1] 侯朝炯,马念杰. 煤层巷道两帮煤体应力和极限平衡区的探讨[J]. 煤炭学报,1989(4):21−29. doi: 10.3321/j.issn:0253-9993.1989.04.001 HOU Chaojiong,MA Nianjie. Stress in in-seam roadway sides and limit equilibrium zone[J]. Journal of China Coal Society,1989(4):21−29. doi: 10.3321/j.issn:0253-9993.1989.04.001
[2] 高 玮. 倾斜煤柱稳定性的弹塑性分析[J]. 力学与实践,2001,23(2):23−26. doi: 10.3969/j.issn.1000-0879.2001.02.005 GAO Wei. Elastoplastic analysis of the stability of inclined coal pillars[J]. Mechanics and Practice,2001,23(2):23−26. doi: 10.3969/j.issn.1000-0879.2001.02.005
[3] 于远祥,洪 兴,陈方方. 回采巷道煤体荷载传递机理及其极限平衡区的研究[J]. 煤炭学报,2012,37(10):1630−1636. doi: 10.13225/j.cnki.jccs.2012.10.017 YU Yuanxiang,HONG Xing,CHEN Fangfang. Study on load transmission mechanism and limit equilibrium zone of coal-wall in roadway[J]. Journal of China Coal Society,2012,37(10):1630−1636. doi: 10.13225/j.cnki.jccs.2012.10.017
[4] 翟 锦. 倾斜煤层区段煤柱宽度留设研究[D]. 西安: 西安科技大学, 2012. ZHAI Jin. Research on the width of coal pillar in inclined coal seam section[D]. Xi'an: Xi'an University of Science and Technology, 2012.
[5] 王德超,李术才,王 琦,等. 深部厚煤层综放沿空掘巷煤柱合理宽度试验研究[J]. 岩石力学与工程学报,2014,33(3):539−548. doi: 10.13722/j.cnki.jrme.2014.03.012 WANG Dechao,LI Shucai,WANG Qi,et al. Experimental study of reasonable coal pillar width of fully mechanized top coal caving in deep thick seam[J]. Chinese Journal of Rock Mechanics and Engineering,2014,33(3):539−548. doi: 10.13722/j.cnki.jrme.2014.03.012
[6] 孔德中,王兆会,李小萌,等. 大采高综放面区段煤柱合理留设研究[J]. 岩土力学,2014,35(S2):460−466. doi: 10.16285/j.rsm.2014.s2.040 KONG Dezhong,WANG Zhaohui,LI Xiaomeng,et al. Study of reasonable width of full-mechanized top-coal caving with large mining height[J]. Rock and Soil Mechanics,2014,35(S2):460−466. doi: 10.16285/j.rsm.2014.s2.040
[7] 张念超,孙元田,蔡胜海,等. 基于统一强度理论的护巷煤柱尺寸与支护技术研究[J]. 煤矿安全,2016,47(6):209−213. doi: 10.13347/j.cnki.mkaq.2016.06.057 ZHANG Nianchao,SUN Yuantian,CAI Shenghai,et al. Research on coal pillar size and support technology of roadway protection based on unified strength theory[J]. Coal Mine Safety,2016,47(6):209−213. doi: 10.13347/j.cnki.mkaq.2016.06.057
[8] 魏 臻, 李晋平, 何富连, 等. 综放面单侧采空煤柱稳定性研究及实测[J]. 矿业科学学报, 2017, 2(4): 371-378. WEI Zhen, LI Jinping, HE Fulian, et al. Stability study and field observation of coal pillar between goaf and unmined top-coal caving face[J]. Journal of Mining Science and Technology, 2017, 2(4). 371-378.
[9] 赵 宾,王方田,梁宁宁,等. 高应力综放面区段煤柱合理宽度与控制技术[J]. 采矿与安全工程学报,2018,35(1):19−26. doi: 10.13545/j.cnki.jmse.2018.01.003 ZHAO Bin,WANG Fangtian,LIANG Ningning,et al. Reasonable width and control technology of coal pillars in high-stress fully mechanized caving face[J]. Journal of Mining & Safety Engineering,2018,35(1):19−26. doi: 10.13545/j.cnki.jmse.2018.01.003
[10] 邸 帅. 8.5 m超大采高综采面回采巷道合理煤柱宽度分析[J]. 地下空间与工程学报,2019,15(2):614−621. DI Shuai. Analysis of the reasonable coal pillar width of the mining roadway in a fully mechanized mining face with an ultra-large cutting height of 8.5m[J]. Chinese Journal of Underground Space and Engineering,2019,15(2):614−621.
[11] 宁 静. 深部大采高综采工作面区段煤柱宽度优化研究[J]. 煤炭工程,2019,51(3):13−17. NING Jing. Tranquility Study on optimization of coal pillar width in deep and large mining height fully mechanized mining face[J]. Coal Engineering,2019,51(3):13−17.
[12] 顾士坦,魏宝贞,蒋邦友,等. 双侧采空不规则煤柱稳定性分析[J]. 煤矿安全,2020,51(3):216−220. doi: 10.13347/j.cnki.mkaq.2020.03.046 GU Shitan,WEI Baozhen,JIANG Bangyou,et al. Stability analysis of irregular coal pillars in bilateral goaf[J]. Safety in Coal Mines,2020,51(3):216−220. doi: 10.13347/j.cnki.mkaq.2020.03.046
[13] 王东星. 浅埋煤层大采高综采面区段煤柱宽度留设理论及试验研究[D]. 西安: 西安科技大学, 2017. WANG Dongxing. Theoretical and experimental research on coal pillar width in the section of a fully mechanized mining face with high mining height in shallow coal seams[D]. Xi’an: Xi’an University of Science and Technology, 2017.
[14] 王宝石. 区段煤柱宽度合理留设研究[D]. 邯郸: 河北工程大学, 2013. WANG Baoshi. Study on the reasonable retainning of the width of the section coal pillar[D]. Handan: Hebei University of Engineering, 2013.
[15] 薛雄飞. 杭来湾煤矿工作面区段煤柱合理尺寸研究[D]. 西安: 西安科技大学, 2020. XUE Xiongfei. Research on reasonable size of coal pillars in working face section of Hanglaiwan Coal Mine[D]. Xi’an: Xi’an University of Science and Technology, 2020.
[16] 白进龙. 浅埋煤层大采高综采面区段煤柱宽度优化研究[D]. 西安: 西安科技大学, 2018. BAI Jinlong. Optimization of Section Coal Pillar Width on Full-mechanized Mining Face in Shallow Coal Seam[D]. Xi’an: Xi’an University of Science and Technology, 2018.
[17] 张艳军. 浅埋煤层大采高综采面区段煤柱合理宽度留设研究[D]. 西安: 西安科技大学, 2016. ZHANG Yanjun. Research on the reasonable width of coal pillars in the section of fully-mechanized mining face with large mining height in shallow seam[D]. Xi’an: Xi’an University of Science and Technology, 2016.
[18] 钱鸣高, 石平五. 矿山压力与岩层控制[M]. 徐州: 中国矿业大学出版社, 2003. [19] 王卫军,黄成光,侯朝炯,等. 综放沿空掘巷底鼓的受力变形分析[J]. 煤炭学报,2002,27(1):26−30. doi: 10.3321/j.issn:0253-9993.2002.01.006 WANG Weijun,HUANG Chengguang,HOU Chaojiong,et al. Mechanical deformation analysis of floor strata of roadway driving along next goaf in fully mechanized sub-level caving face[J]. Journal of China Coal Society,2002,27(1):26−30. doi: 10.3321/j.issn:0253-9993.2002.01.006
[20] 薛 强. 弹性力学[M]. 北京: 北京大学出版社, 2006. [21] 秦忠诚,王同旭. 深井孤岛综放面支承压力分布及其在底板中的传递规律[J]. 岩石力学与工程学报,2004,23(7):1127−1131. QIN Zhongcheng,WANG Tongxu. Abutment pressure distribution and its transfer law in floor of deep isolated fully-mechanized mining faces using sub-level caving[J]. Chinese Journal of Rock Mechanics and Engineering,2004,23(7):1127−1131.
[22] 黄庆享,刘建浩. 浅埋大采高工作面煤壁片帮的柱条模型分析[J]. 采矿与安全工程学报,2015,32(2):187−191. doi: 10.13545/j.cnki.jmse.2015.02.003 HUANG Qingxiang,LIU Jianhao. Analysis of the column model of the coal wall slab in the shallow burying and large mining height face[J]. Journal of Mining & Safety Engineering,2015,32(2):187−191. doi: 10.13545/j.cnki.jmse.2015.02.003
[23] 孙训方, 方孝淑. 材料力学(I)[M]. 北京: 高等教育出版社, 2009. [24] 尹希文,闫少宏,安 宇. 大采高综采面煤壁片帮特征分析与应用[J]. 采矿与安全工程学报,2008,25(2):50−52. doi: 10.3969/j.issn.1673-3363.2008.02.021 YIN Xiwen,YAN Shaohong,AN Yu. Characters of the rib spalling in fully mechanized caving face with great mining height[J]. Journal of Mining & Safety Engineering,2008,25(2):50−52. doi: 10.3969/j.issn.1673-3363.2008.02.021
[25] 宁 宇. 大采高综采煤壁片帮冒顶机理与控制技术[J]. 煤炭学报,2009,34(1):50−52. doi: 10.3321/j.issn:0253-9993.2009.01.010 NING Yu. Mechanism and control technique of the rib spalling in fully mechanized mining face with great mining height[J]. Journal of China Coal Society,2009,34(1):50−52. doi: 10.3321/j.issn:0253-9993.2009.01.010
[26] 沈明荣. 岩体力学[M]. 上海: 同济大学出版社, 2002 [27] 陈建君. 厚硬顶板特厚煤层孤岛煤柱应力集中程度及错层防冲研究[D]. 徐州: 中国矿业大学, 2016. CHEN Jianjun. Study on the stress concentration and the anti-scouring of the isolated island coal pillar in the thick and hard roof extra-thick coal seam[D]. Xuzhou: China University of Mining and Technology, 2016.
[28] 董方庭. 巷道围岩松动圈支护理论及应用技术[M]. 北京: 煤炭工业出版社, 2001. [29] 于远祥. 矩形巷道围岩变形破坏机理及在王村矿的应用研究[D]. 西安: 西安科技大学, 2013. YU Yuanxiang. Study on deformation mechanism of surrounding rock in rectangular roadway and its application in Wangcun coal mine [D]. Xi'an: Xi'an University of science and technology, 2013.
-
期刊类型引用(11)
1. 罗波远，张涛，安彦成，姚新宇，刘志刚，李树文. 厚煤层大巷保护煤柱合理尺寸优化. 煤炭技术. 2025(01): 48-51 . 
百度学术
2. 郭东明，肖博丰，赵志峰，叶贵川，刘嘉华，唐耿福. 新型叠片式吸能锚杆与普通高强度锚杆抗冲击性能对比研究. 煤矿安全. 2025(02): 109-117 . 百度学术
3. 杨敬轩，张瑞，卢硕，孙兴平，陈建本. 工作面顶板超前预裂应力阻隔护巷机理. 中国矿业大学学报. 2024(02): 264-276 . 百度学术
4. 左常清，刘亚楠，田晓伟. 义桥煤矿边界煤柱留设合理性研究. 山东煤炭科技. 2024(05): 31-36 . 百度学术
5. 马斌文，谢和平，张修峰，周宏伟，陈洋，郑福润，孙文斌，朱建波. 动载扰动下巷道围岩冲击破坏与能量释放规律研究. 采矿与岩层控制工程学报. 2024(04): 5-22 . 百度学术
6. 李庆祥，张晓广，刘道园，赵旭晔. 不同岩性煤岩组合体力学特性模拟研究. 煤. 2024(11): 53-56 . 百度学术
7. 马小利. 双陷落柱孤岛工作面托顶煤巷关键区域防控研究. 煤炭与化工. 2024(11): 16-19+23 . 百度学术
8. 杨敬轩，朱乐章，周萌，卢硕. 顶板应力作用下的覆岩切顶孔间距确定及切顶线位置选择. 采矿与安全工程学报. 2024(06): 1170-1178 . 百度学术
9. 谷拴成，张炜，文嘉豪，康恽博，呼嘉龙. 防隔水煤柱合理留设宽度分析. 西安科技大学学报. 2024(06): 1060-1070 . 百度学术
10. 李怀展，唐超，郭广礼，陈福，李伟，周华安，黄建勇. 热力耦合作用下煤炭地下气化地表沉陷预测方法. 煤炭科学技术. 2023(10): 242-251 . 本站查看
11. 韦四江，王生柱，李鑫鹏. 煤柱留设宽度对密闭墙稳定性的影响研究. 煤炭技术. 2023(12): 1-5 . 百度学术
其他类型引用(3)
计量
- 文章访问数: 96
- HTML全文浏览量: 8
- PDF下载量: 78
- 被引次数: 14
