Influence characteristics of regional micrometeorology on macroscopic scale of external dump
-
摘要:
为了研究露天煤矿外排土场对周围微气象因子的影响关系，进而为矿区生态修复提供基础，以矿区周围气流运动及水分分布为影响因子，进一步探索外排土场尺度及形态的生态效应。以红沙泉露天煤矿为例，运用Fluent流体模拟软件，采用大涡模拟和组分传输模型，分别建立不同高度与角度的外排土场模型进行模拟。通过红沙泉露天煤矿所处地区气象站点，得到该地区2020年逐月温度与露点数据，从而得出该地区逐月凝结高度。分析不同条件下得到的气流速度云图，水分质量分数云图，获得不同外排土场尺寸对周围气流运动及水分分布的影响规律，并与该地区凝结高度进行对比。数值模拟结果表明：随着外排土场高度增加，气流爬坡速度缓慢增加，爬升最大高度也随之增长，迎风坡以及外排土场顶部的水分质量分数均增大，背风坡形成越来越明显的卡门涡街现象，涡流数量增加且影响范围变大，涡流区域内气流速度、水分质量分数均增加并达到最大值，外排土场角度增大对气流最大爬升高度影响较大，对周围微气象因子影响不大。在极限高度360 m条件下，当外排土场角度达到26°临界值时，会使气流最大爬升高度达到该地区夏季凝结高度；或在外排土场极限角度22°条件下，当外排土场高度达到380 m临界值时，会使气流最大爬升高度达到该地区夏季凝结高度，从而对矿区周围降水产生促进作用。
Abstract:To study the relationship of external dump of opencast coal mine impacting on surrounding micrometeorological factors, and then provide the foundation acknowledge for the ecological restoration in mining area, the airflow movement and water distribution around the mine area were regarded as the influencing factors to furtherly explore the ecological effects of the scale and form of the external dump. Taking Hongshaquan open-pit coal mine as an example, using Fluent fluid simulation software, adopting large eddy simulation and component transfer model, the external dump models with different heights and angles were established for simulation. The monthly temperature and dew point data of the area in 2020 were obtained through the regional meteorological station where Hongshaquan open-pit coal mine located, thereby reaching the monthly condensation height of the area. By analyzing the cloud diagrams of air velocity and moisture mass fraction under different conditions, the influence law of different sizes of external dump on the surrounding air flow movement and water distribution was obtained, and compared with the condensation height in the area. The numerical simulation results showed that: with the increasing dump height, the air climbing speed increased slowly, while the maximum climbing height increased, as well as the moisture mass fraction on the windward slope and the dump top ascended. More and more obvious Karman vortex street phenomenon was formed along the leeward slope. The number of vortices increased and the influence range became larger. Both of the airflow velocity and moisture mass fraction in the vortex area increased and reached their maximum values. The increase in the angle of the external dump presented a significant impact on the maximum climb height of the airflow, but showed little impact on the surrounding micrometeorological factors. Under the condition of 360 m limit height, when the external dump Angle reached the critical value of 26°, the maximum airflow climbing height would reach the summer condensation height in this area. Alternatively, under the condition of 22° limit Angle of the outer dump, when the height of the outer dump reached the critical value of 380 m, the maximum climbing height of airflow would reach the summer condensation height of the region, thereby promoting the precipitation around the mining area
-
0. 引 言
我国中东部煤炭资源逐渐濒临枯竭，煤炭开采战略重心逐渐向西部地区转移[1]，随着采煤技术与装备水平的发展，矿井开采强度逐渐增大，开采扰动强度过大时产生的动力灾害更具冲击性和突然性[2]，致使冲击地压已成为目前威胁煤矿安全最为突出的动力灾害之一[3]。而影响冲击地压的因素众多，矿井开采强度作为众多开采技术影响因素之一，近些年来随着对冲击地压研究工作的深入与细化，逐渐受到重视。开采强度的不同主要体现在煤层采高、工作面宽度、推进速度等参数的差异。在采高、工作面宽度等确定的情况下，主要受推进速度的影响，特别是近距离煤层开采存在实体煤下单次采动与采空区下重复采动的复杂空间影响，煤岩赋存条件特殊，不同采动强度下动静载叠加诱发强矿压显现的机制复杂，因此亟需开展冲击地压矿井近距离特厚煤层综放工作面推进速度研究。近距离煤层指井田开采范围内相邻两层煤的层间距离较小，开采时相互影响的煤层[4]，一般采用下行开采，在上层煤开采结束后，其上部覆岩垮落，下层煤二次开采造成重复采动，对覆岩结构的整体性造成严重影响[5]。在工作面推进速度研究方面，现场实践发现，工作面高速推进易诱发冲击地压等动力灾害，诸多学者根据煤层资源赋存条件等取得不同的研究成果，其中张宏伟等[6]认为随推进速度提高，煤体应力升高速率加快，煤体脆性增强，进而导致煤岩冲击倾向性增强。王家臣等[7]认为工作面推进速度提高，基本顶断前积聚的弹性应变能及初始动能增加，发生破断失稳的可能性越大。谭云亮等[8]分析不同开采因素与所产生动能之间关系，得到开采进尺与所产生的动能之间呈正相关关系。冯龙飞等[9]研究发现推进速度越快围岩单位时间内积聚的能量越多，结合微震监测分析得到大采高工作面安全临界推进速度为4 m/d。笔者[10-11]通过对采动影响的微震监测分析，得到宽沟煤矿B4-1煤层推进速度为4 m/d时较为合宜。谢广祥等[12]认为推进速度越快，超前支承压力峰值越大，且峰值位置距离煤壁越近。闫宪磊[13]指出推进速度越快，煤岩发生冲击地压危险性越高。刘金海等[14]揭示了工作面高速推采、非匀速开采时冲击危险性较大，易诱发冲击地压。赵同彬等[15]通过微震活动表现出的低能高频现象认为推进速度为4.0 m/d左右时，更有利于覆岩能量释放。马海峰等[16]发现随着推进速度的增加，煤岩的峰值强度呈不断增大趋势。崔峰等[17]基于加卸载响应比研究认为冲击地压矿井急倾斜巨厚煤层工作面推进速度为4.8 m/d时能较好实现安全高效开采。
上述学者通过理论分析、模拟试验、现场监测等不同方法，对工作面推进速度影响下的采场围岩矿压及冲击危险性进行了有益的探索，为后续相关研究奠定了坚实的基础，但针对冲击地压矿井近距离煤层综放工作面推进速度的研究相对较少。因此，采用数值模拟、理论分析等方法，对冲击地压矿井近距离煤层综放工作面合理推进速度展开综合研究，为类似工程条件下冲击地压灾害防治提供理论基础。
1. 工程概况及数值模型建立
1.1 工程概况
宽沟煤矿隶属国能新疆宽沟矿业有限责任公司，为典型的冲击地压矿井，行政区隶属昌吉州呼图壁县管辖，井田东西长约9.7 km，南北宽约3.15 km，井田面积约20.132 5 km2，矿井主要赋存B4-1和B2两层煤，其中B4-1煤层具有强冲击倾向性，顶板具有弱冲击倾向性，煤层平均厚度3.0 m，平均倾角14°，其主采的W1145工作面倾向宽度约175 m，走向长度约746 m，目前已回采结束。在B4-1煤层下方约43.8 m的B2煤层具有弱冲击倾向性，顶板具有强冲击倾向性，煤层平均埋深约392 m，煤层平均厚度9.5 m，煤层倾角为12°～14°，顶板多为粉砂质泥岩、粉砂岩；底板多为粉砂岩、中砂岩，结构简单，不含夹矸，属稳定特厚煤层。其主采的W1123工作面倾向宽度约192 m，走向长度约1468 m，采用综放开采，采高3.2 m，放煤高度6.3 m，采放比为1∶1.97，B4-1煤层与B2煤层工作面分布情况如图1所示。
由图1可知，W1123工作面回采面临着在上覆B4-1煤层实体煤下单次采动与W1145工作面采空区下重复采动的复杂环境，上覆B4-1煤层实体煤与采空区交界作为边界支承点，承受着采空区覆岩破断后传递的应力，在实体煤一定区域内形成支承压力升高区，在其下方岩层一定范围内形成应力集中区，随着B2煤层W1123工作面的开采，工作面的采动应力和W1145采空区交界实体煤引起的应力集中区叠加，形成一个新的应力集中区，该区域的形成在一定程度上会加剧下部W1123工作面矿压显现强度，煤层交界覆岩结构特征如图2所示。
![]()
图 2 采动覆岩结构及应力分布Figure 2. Overlying rock structure and stress distribution during mining1.2 数值分析模型的构建及方案设计
鉴于3DEC数值模拟分析能够很好反映工作面覆岩运移垮落、能量积聚特征等，根据宽沟煤矿地质条件，借助3DEC数值模拟软件，建立并基于现场实测数据校核了数值模型，采用理想弹塑性本构模型莫尔−库仑准则判断岩体的破坏。
通过数值模拟试验，对W1123工作面不同推进速度下的采场围岩矿压规律及冲击危险性进行研究，其W1123工作面区域钻孔局部综合柱状（ZK201）如图3所示。根据现场地质调查与岩石力学试验结果综合确定煤岩物理力学参数，数值模拟实验中煤层顶底板附近区域主要煤岩物理力学参数见表1。
![]()
图 3 W1123工作面局部钻孔综合柱状（ZK201）Figure 3. Comprehensive histogram of local drilling in W1123 working face (ZK201)表 1 部分岩层煤岩力学参数Table 1. Mechanical parameters of rock strata and coal rocks岩性 密度/(kg·m−3) 抗拉强度/MPa 弹性模量/MPa 抗压强度/MPa 泊松比 体积模量/GPa 剪切模量/GPa 黏聚力/MPa 内摩擦角/(°) 砂砾岩 2 532 2.33 25 400 34.32 0.27 14.17 8.50 8.92 28 B4-1煤 1 303 1.97 7 485 24.80 0.35 6.25 2.88 3.81 21 粗砂岩 2 541 4.48 25 820 48.23 0.19 13.45 10.94 9.73 29 泥岩 2 533 2.28 15 420 37.65 0.20 7.42 3.83 4.39 30 细砂岩 2 592 3.17 33 420 57.60 0.18 19.57 14.07 6.58 27 B2煤 1 640 2.15 4 923 26.76 0.36 4.92 2.74 4.30 24 建立宽沟煤矿走向数值模拟模型的尺寸为：长×宽×高为1 600 m×60 m×378 m，采用摩尔−库伦准则判断煤岩体的破坏，将模型顶部未模拟岩层换算成5.4 MPa的等效载荷施加于模型顶部，四周和底边限制变形，得到W1123工作面走向数值计算模型如图4所示，为去除边界效应，B1煤层及B2煤层边界煤柱均留设65 m。
设计模拟方案采用一次采全高的开采方法，对于冲击地压矿井的推进速度研究，业内学者[7,9-10,15]大多采用4 m/d作为安全或中等推进速度。随着矿井机械化水平的提高及智能监测设备的应用，一定程度上实现了对冲击地压等事故的精准防控，在保证安全生产的前提下，部分矿井将提高工作面推进速度作为提升矿井生产能力的首选措施。笔者对此加以充分考虑，在综合考虑上述因素的基础上，此次将推进速度暂定为4.0 m/d以下、4.0～8.0 m/d，8.0 m/d以上3个区域范围，将推进速度划分为低速、中速、高速3种推进速度。同时为对比工作面推进至同一位置时推进速度对采动应力、能量积聚、覆岩破断特征及煤岩冲击危险性的影响规律，结合宽沟煤矿采煤机滚筒0.8 m的单次截深，设计3.2、6.4、9.6 m/d作为低速、中速、高速推进速度，根据设计的3种推进速度，将模型W1123工作面煤层走向开挖区域网格划分采用3.2、6.4、9.6 m/d推进速度下的最大公因子3.2 m，以此来保证同一数值模型进行不同推进速度下试验效果分析，进一步提高了数值模拟的准确性，在靠近煤层附近区域的顶板岩层中，为了使数值模拟效果更接近于实际，更好地模拟出覆岩运移破断等效果，块体尺寸相对较小，随着逐渐远离煤层，块体尺寸略微增大。并据此展开近距离煤层综放工作面不同推进速度对采动应力、能量积聚、覆岩破断特征及冲击危险性分析。
2. 不同推进速度对采动应力变化与覆岩能量积聚特征的影响
地下煤体开采对煤体来说是卸载过程，相反对周围岩体是加载过程，推进速度作为开采强度的主要影响因素之一，其对煤岩体内部应力、能量演化影响显著，为此本次采用数值模拟试验分析不同推进速度下的超前支承压力分布、覆岩的弹性能积聚等规律。
2.1 推进速度对采动应力的影响
由于W1123工作面煤层赋存空间复杂，地下开采活动的进行造成采场围岩应力的转移及重新分布，破坏了原有的应力平衡状态，推进速度越快，周围岩体的加载速率越大[18]。推进速度提高导致周围岩体加卸载活动剧烈，加剧了围岩局部应力集中现象，易诱发覆岩发生破断失稳。
为研究工作面在实体煤下回采和采空区下回采的超前支承压力演化规律，分析了工作面以不同推进速度推进至相同位置时工作面前方的顶板应力演化规律，并选取工作面推进至280、600、800、1 000、1 300 m五个位置，绘制出超前支承压力演化趋势如图5所示。
![]()
图 5 不同推进速度超前支承压力演化规律Figure 5. Evolution law of advance abutment pressure at different advancing speeds由图5可知，工作面分别在实体煤与采空区下回采时，超前支承压力峰值均随工作面向前推进不断增大，其影响范围随工作面的推进不断向前推移，峰值应力随之向前推移，且随推进速度的提高，超前支承压力峰值也随之不断增大。实体煤下超前支承压力明显大于采空区下，在接近上覆B1煤层W1145工作面开切眼后方100 m左右范围内出现了较为明显的应力集中现象，超前支承分布近似呈现“双峰值”现象。主要是由于上覆W1145工作面回采后，其开切眼左侧实体煤承载了部分未垮落覆岩结构以及边界未压实的铰接结构的应力，W1123回采形成的采动应力与W1145工作面采空区的残余支承压力叠加，致使应力水平升高，当工作面推进至上覆W1145开切眼附近时，超前支承压力峰值快速上升，峰值应力达到最大。
工作面在实体煤下回采时，随着工作面推进速度的不断提高，超前支承压力峰值不断增大。从图5a可以看出，当工作面以3.2 m/d回采时，超前支承压力逐渐由40.85 MPa增加至55.56 MPa，增加了14.71 MPa，增幅约36%；从图5b推进速度提高至6.4 m/d时可以看出，超前支承压力逐渐由46.94 MPa增加至74.48 MPa，增加了27.54 MPa，增幅约58.67%；从图5c推进速度提高至9.6 m/d时可以看出，超前支承压力逐渐由50.35 MPa增加至111.36 MPa，增加了61.01 MPa，增幅约121.17%。结合图5a、图5b、图5c可以看出工作面推进至同一位置时的超前支承压力随推进速度的提高明显升高，当工作面推进至280、600、800 m时，超前支承压力峰值随推进速度的提高分别从40.85 MPa增加至50.35，44.39 MPa加至55.50，55.56 MPa增加至111.36 MPa，3个不同位置的超前支承压力峰值随推进速度的提升其幅值分别增加了23.26%、25.03%、100.43%。
而在采空区下回采时，由于W1145采空区的存在，对采空区与W1123工作面中间的岩层起到一定程度的卸荷作用，使得W1123工作面上覆岩层承担的应力相对比较小，当工作面以3.2 m/d回采时，超前支承压力逐渐由24.86 MPa增加至32.90 MPa，增加了8.04 MPa，增幅约32.34%；当推进速度提高至6.4 m/d时，超前支承压力逐渐由37.49 MPa增加至41.18 MPa，增加了3.69 MPa，增幅约9.84%；当推进速度提高至9.6 m/d时，超前支承压力逐渐由42.33 MPa增加至43.26 MPa，增加了0.93 MPa，增幅约2.2%。
工作面推进至不同位置时，超前支承压力峰值与不同推进速度的关系如图6所示。由图6可知，工作面以3.2、6.4、9.6 m/d的推进速度推进至同一位置时的平均峰值应力分别为39.75、50.35、60.56 MPa，中速（6.4 m/d）、高速（9.6 m/d）推进速度下的峰值应力比低速分别增加约26.67%、52.35%。峰值应力随着推进速度的增加呈明显的非线性递增趋势，且其增幅逐渐增加，其非线性表达式为：
$ y=335.101\;57- 306.337\;71 {\rm{e}}^{(-0.011\;42 v)}$ 。![]()
图 6 峰值应力与不同推进速度的关系Figure 6. Relationship between peak stress and different advancing speed从上述分析可知，推进速度越快，超前支承压力峰值愈大，峰值应力随推进速度增加呈明显的非线性递增趋势，且其增幅逐渐增加。在实体煤下回采时，超前支承压力峰值随推进速度的提高不断增大，其增幅随之不断增加，推进速度较低时，采场围岩内部应力水平相对较为缓和。而在采空区下回采时，超前支承压力较实体煤下明显降低，随推进速度的提高，峰值应力水平呈现略微升高的趋势，但与在实体煤下回采时的应力水平相比还相差甚远，其增幅随推进速度的提高呈现逐渐降低的趋势。
2.2 推进速度对覆岩能量积聚的影响
推进速度的不同造成上覆顶板岩层的应力有所不同，进而导致覆岩内能量变化不一，煤岩体内的能量聚集和突然释放容易诱发动力灾害事故。相关研究表明覆岩破断一般超前于工作面发生[19]。因此，绘制如图7所示工作面回采过程中的覆岩结构演化特征。
![]()
图 7 工作面覆岩结构特征Figure 7. Structural characteristics of overlying strata in working face同时建立如图8所示的顶板岩层简支梁模型。
其弯矩M(x)为
$$ M(x)=\frac{q b x}{2}(L-x) $$ (1) 式中：q为上覆岩层均布荷载，MPa；x为截面到支点的距离，m；L为极限长度，m；b为岩梁宽度，取1 m。
通过积分计算得到该简支梁积聚的弹性能W为
$$ W=\int_0^L \frac{M^2(X)}{2 E I} {\rm{d}} x=\frac{q^2 L^5 b^2}{240 E I} $$ (2) 式中：E为岩层的弹性模量，GPa；I为转动惯量kg·m2。
由式（2）可知，覆岩积聚的能量与覆岩悬空长度呈正相关性，而覆岩的极限悬空长度与推进速度密切相关，推进速度越快，覆岩悬空长度越长，随着覆岩悬空长度的增加，覆岩积聚的弹性能显著增加。结合表1所示岩层物理力学参数，根据3DEC数值模拟试验研究得出不同推进速度下顶板的初次垮落破断距分别为54.4、64、76.8 m，取其上方均布载荷q=8 MPa，转动惯量I=h3/12，并计算覆岩初次破断时的弹性能见表2。
表 2 不同推进速度覆岩初次破断距及弹性能特征Table 2. Initial breaking distance and elastic properties of overlying strata at different advancing speeds推进速度/（m•d−1） 岩层 初次破断距/m 弹性能/kJ 3.2 厚硬直接顶 54.4 11.79 6.4 厚硬直接顶 64.0 26.57 9.6 厚硬直接顶 76.8 66.11 从表2中得出，随着推进速度的不断提高，覆岩顶板初次垮落步距不断增加，覆岩内积聚的弹性能随之不断增加，工作面以3.2、6.4、9.6 m/d推进时覆岩初次破断时覆岩内积聚的弹性能分别为11.79、26.57、66.11 kJ。中速（6.4 m/d）、高速（9.6 m/d）推进速度下的覆岩弹性能分别是低速（3.2 m/d）时的2.25倍、5.61倍。
同时结合三向应力状态下岩体应变能密度公式[20]，计算煤岩体单位积内积聚的弹性能：
$$ U=\frac{1}{2 E}\left[\sigma_1{ }^2+\sigma_2{ }^2+\sigma_3{ }^2-2 \mu\left(\sigma_1 \sigma_2+\sigma_2 \sigma_3+\sigma_1 \sigma_3\right)\right] $$ (3) 式中：U为岩体弹性势能，J/m3；
$ \sigma_1 $ ，$ \sigma_{2} $ ，$ \sigma_{3} $ 为岩体三向应力，Pa；μ为岩体的泊松比。将W1123工作面不同推进速度推进至覆岩初次破断时覆岩应力及煤岩物理力学参数代入式（3），计算并绘制出如图9所示不同推进速度覆岩初次破断时覆岩应变能密度演化趋势。
![]()
图 9 不同推进速度覆岩初次破断时覆岩应变能密度演化趋势Figure 9. Evolution trend of strain energy density of overlying strata during initial fracture of overlying strata at different advancing speeds从图9中可以看出，不同推进速度推进至覆岩初次破断时，应变能密度演化趋势呈相似性，且实体煤下应变能密度远高于上覆W1145采空区下。采空区内由于应力释放导致应变能密度极低，当W1123工作面以3.2 m/d推进时在工作面煤壁前方应变能密度出现激增现象，应变能密度为43.82 kJ/m3，而后出现缓慢下降趋势，在上覆W1145工作面开切眼后方约100 m位置时应变能不断增加，在距离开切眼约50 m时应变能再次出现激增现象，应变能密度为30.27 kJ/m3，与实体煤下超前支承压力分布规律的“双峰值”现象相对应，在实体煤下应变能密度呈现出先减后增的演化趋势；当推进速度提高至6.4 m/d时，煤壁前方应变能密度激增至53.10 kJ/m3，在距离W1145开切眼约50 m时应变能密度激增至36.10 kJ/m3；而当推进速度提高至9.6 m/d时，煤壁前方应变能密度激增至64.71 kJ/m3，在距离W1145开切眼约50 m时应变能密度激增至37.17 kJ/m3。在上覆W1145开切眼附近，应变能密度均出现骤降现象，而后缓慢增长，在停采线附近出现缓慢下降趋势，在W1145采空区下应变能密度均呈现出先增→平稳过渡→后减的演化趋势。
综上可知，随着工作面推进速度的提高，覆岩初次破断距越大，顶板悬空距离越长，覆岩内积聚的能量不断增加，初次破断时应变能密度随之不断增加，同时工作面不同推进速度推进至覆岩初次破断时的应变能密度演化趋势可划分为采空区应变能密度降低阶段、工作面前方应变能密度激增、实体煤下先降低后缓慢增长阶段、W1145采空区下先增后减的演化趋势。
3. 不同推进速度的覆岩运移破断特征及冲击危险性分析
不同推进速度对采动应力与覆岩能量积聚影响较为显著，随着煤岩体采出后，采场应力演化与能量积聚造成覆岩集中破断，引起工作面的矿压显现，其显现的强度与覆岩运移特征密切相关，覆岩运移破断是引起动力失稳、造成冲击的主要原因，因此分析不同推进速度下覆岩破断特征及煤岩冲击危险性对于矿井安全、高效生产尤为重要。
3.1 推进速度对覆岩运移破断特征的影响
随着煤层采出后，其覆岩空间演化特征普遍发生变化，主要表现在围岩开始活动，引起覆岩结构上的变形、破坏和相对位置上的改变，随着开采的不断扰动与破坏，进一步造成覆岩运移破断，而工作面的矿压显现主要是由采动覆岩集中破断运动引起的，为研究推进速度对覆岩破断的影响，选取3种速度推进至相同位置（以推进至141.8 m时）覆岩演化情况进行分析。笔者以位移云图分析覆岩运移与贯通情况，进而得到覆岩变形破坏情况。
由图10工作面不同推进速度的覆岩破坏特征可以看出，当工作面以3.2 m/d的速度推进至141.8 m时，采空区覆岩中部向上产生近似三角形状的发育趋势，最大位移约2.75 m，距离顶板愈远，位移变化量越小；当推进速度提高至6.4 m/d时，采空区左侧覆岩出现明显下沉，覆岩下沉进一步向上扩展，区域进一步加大，采空区中部下沉较为明显的区域不断向上扩展，最大位移为4.68 m；而当速度提高至9.6 m/d时，采空区左侧覆岩下沉愈加强烈，除了采空区覆岩顶板位移变化强烈外，上覆岩层位移变化不断增大，并与周围岩体相互贯通，最大变化量为6.86 m；从位移云图可以看出，采空区右侧边界与上覆W1145工作面开切眼形成近似半椭球形包络线，当推进速度为3.2、6.4 m/d时，采空区覆岩右侧端部与W1145工作面开切眼形成相连通的2个层次分明的半椭球形包络线，当推进速度提高至9.6 m/d时，采空区覆岩上部仅剩一个与W1145工作面开切眼相连通的包络线，且覆岩变化量大、运移破断剧烈，呈现出明显的非连续变形发育特征，采空区下部覆岩运移演化快速向上传递并逐渐贯通。
![]()
图 10 工作面不同推进速度的覆岩破坏特征Figure 10. Failure characteristics of overlying strata at different advancing speeds in working face综上分析，随着推进速度的增加，覆岩位移变化量不断增大，即覆岩破坏高度及范围进一步扩大，呈现出明显的非连续变形发育特征。由于W1123工作面从开切眼开始，需经历实体煤下−上覆W1145工作面采空区下2个阶段回采过程，因此，选取不同推进速度实体煤下的覆岩破断形态进行对比分析，不同推进速度下W1123工作面在实体煤下回采时覆岩破断形态如图11所示。
![]()
图 11 工作面在实体煤下回采覆岩破坏特征Figure 11. Failure characteristics of overlying strata in mining face under solid coal从图11可知，工作面在实体煤下回采结束时，覆岩破断下沉均呈现出梯形结构，在采空区覆岩两端中下部位移变化量交替区域（红、绿色变化交替区）裂隙较为发育，随着推进速度的增加，覆岩破断范围进一步加大，其主要体现在采空区右上方覆岩颜色交替变化区域内，推进速度由3.2 m/d逐渐增加到9.6 m/d时，采空区右侧区域内（方框及椭圆圈内）的位移进一步加大，范围进一步扩展，说明随着推进速度提高，覆岩破断范围进一步加大，导致采空区覆岩变化量进一步加大，覆岩运移愈加强烈，直至W1123工作面回采结束，覆岩破断形态如图12所示。
![]()
图 12 工作面回采结束时覆岩结构破断形态Figure 12. Fracture morphology of overlying rock structure at the end of mining in working face从图12中可以看出，覆岩破断下沉呈现出梯形结构，同时由于上覆B1煤层W1145工作面的存在，使得模型右侧区域覆岩的变形明显高于实体煤下，在采空区两侧中下部位移变化量交替区域，裂隙离层较为发育。随着推进速度的提高，覆岩右侧端部及W1145开切眼附近区域上覆岩层位移波动较大的范围及变化量进一步加大，当推进速度由3.2 m/d增加至9.6 m/d时，右侧端部覆岩位移变化量逐渐加大，W1145开切眼附近上覆岩层位移变化量也随推进速度的增加不断增大，且由于开切眼的存在，导致变化量较大的区域不断向左边实体煤上方覆岩扩展，进一步说明随着推进速度的增加，覆岩位移变化量大、破断范围广，覆岩运移剧烈的特点。
3.2 不同推进速度冲击危险性分析
煤岩冲击倾向性是评估及预测冲击地压危险程度的主要依据之一[21]。由于宽沟煤矿B2煤层顶板具有强冲击倾向性，通常煤岩体冲击倾向性越强，煤层开采引起采场应力重新分布，易在局部形成高度应力集中现象，越易诱发冲击地压。
业界学者公认的冲击危险性工程判据可概括为“引起冲击的应力至少应大于煤岩单轴抗压强度且煤岩具有冲击倾向性”，通常将围岩外应力与围岩抗压强度的比值记为Im，引起冲击地压发生的必要条件是围岩外应力超过煤岩单轴抗压强度，即Im>1。将应力比Im在1～1.5，1.5～3、大于3时相对应的冲击危险性划分为弱、中等、强冲击3种危险性[22]。为方便表述，将应力比值Im称为冲击危险性指数。计算公式如下：
$$ I_{\mathrm{m}}=\frac{\sigma_{\mathrm{p}}}{\sigma_{\rm{c}}}$$ (4) 式中：
$ \sigma_{\mathrm{p}} $ 为前方煤体支承压力，MPa；$\sigma_{{\rm{c}}}$ 为煤样抗压强度，MPa。通过数值模型监测工作面推进至不同区域的超前支承压力峰值，结合B2煤体的单轴抗压强度（26.76 MPa），根据式（4）计算并绘制如图13所示冲击危险性指数分布。
由图13可知，在工作面不同回采区域，由于工作面推进过程中扰动程度不一，造成不同区域冲击危险性指数呈现波浪式分布，在初始采动及采空区下回采时，冲击危险性较为缓和。当工作面以3.2 m/d的速度推进时，冲击危险性逐渐从无冲击−弱冲击−中等冲击区域过渡，在冲击高危区域，最大冲击危险性指数约为2.0；当推进速度提高至6.4 m/d时，冲击危险性从弱冲击−中等冲击区域过渡，回采全过程冲击危险性主要以中等冲击为主，在高危险冲击区域，最大冲击危险性指数接近3，即将触碰到强冲击的临界值；当推进速度提高至9.6 m时，冲击危险性明显升高，逐步从弱冲击−中等冲击−强冲击区域过渡，在接近上覆W1145开切眼附近应力集中区域时，最大冲击危险性指数约为4.2，冲击危险性异常升高。工作面在实体煤下回采时，随着推进速度的提高，冲击危险性不断增加，在距离W1145工作面开切眼约100 m时冲击危险性开始稳步上升，变化幅度较大，上覆W1145采空区与实体煤的交界区域为明显的高危险区域，工作面距离上覆W1145开切眼愈近冲击危险性越高。工作面过渡到W1145采空区下回采时，冲击危险性呈现缓慢下降的趋势，相较于实体煤下回采，采空区下回采时冲击危险性明显较低，冲击危险性随推进速度的变化幅度也明显小于实体煤下。
由上分析可知，随推进速度的提高，工作面冲击危险性明显增加，尤其表现在距W1145开切眼约100 m区域，冲击危险性稳步上升，在开切眼附近冲击危险性快速升高，随推进速度由3.2 m/d逐渐递增至9.6 m/d的过程中，冲击危险性指数随之由2.0递增至4.2，由中等冲击危险性跃升至强冲击危险性。
4. W1123工作面推进速度合理确定
通过上述推进速度对超前支承压力、覆岩能量积聚、覆岩运移破断特征及冲击危险性分析的研究发现：推进速度越快，超前支承压力峰值愈大，采动影响范围越广，峰值应力随推进速度增加呈明显的非线性递增趋势，且其增幅逐渐增加。在实体煤下回采超前支承压力峰值随推进速度的提高不断增大，其增幅随之不断增加，而在采空区下回采超前支承压力随推进速度提高呈现略微升高的趋势，其增幅随推进速度提高呈现逐渐降低的趋势；同时推进速度提高覆岩初次破断距越大，破断时释放的能量随之增加，将覆岩初次破断时应变能密度演化趋势划分为采空区应变能降低阶段、应变能激增阶段、实体煤下先降后增过渡阶段、W1145采空区下先增后降的演化趋势。同时推进速度越快，覆岩位移变化量愈大、破断范围广，运移演化剧烈，呈现出非连续变形发育特征，采空区下部覆岩的运移演化很快向上传递并相互贯通。
工作面冲击危险性随推进速度提高显著增加，实体煤下冲击危险性指数显著高于采空区下，冲击危险性随推进速度的变化也明显大于采空区下，尤其表现在距W1145开切眼约100 m区域，冲击危险性稳步上升，在开切眼附近冲击危险性快速升高，当推进速度大于6.4 m/d时，冲击危险性指数大于3，处于强冲击危险性区域，相较于中等冲击危险性，当工作面煤岩体处于强冲击危险性范畴时，煤岩体内部应力高度集中，通常造成采场围岩难以支护等情况，此时需要加大围岩卸压力度，加强支护强度，以防出现冲击地压等动力灾害。通过上文分析，汇总出近距离煤层综放工作面推进速度综合影响统计表见表3。
表 3 近距离煤层推进速度综合影响统计Table 3. Statistical of comprehensive influence of advancing speed of contugous coal seam推进速度/(m•d−1) 3.2 6.4 9.6 采动应力/MPa 39.75 50.35 60.56 初次破断能量积聚/kJ 11.79 26.57 66.11 初次破断应变能密度/(kJ·m−3) 43.82 53.10 64.71 位移变形/m
（推进至141.8 m）2.75 4.68 6.86 危险性指数 2.0 2.8 4.2 针对宽沟煤矿近距离特厚煤层综放工作面推进速度效应展开综合分析，通过研究B2煤层W1123工作面不同推进速度下超前支承压力、能量积聚、覆岩破断、煤岩冲击危险性等演化特征，综合研究认为W1123工作面在保证安全高效生产前提下的相对适宜的推进速度应不大于6.4 m/d。
5. 结 论
1）随着推进速度的提高，峰值应力不断增大，呈明显的非线性递增趋势，其增幅也随之逐渐增加。随推进速度的提高，覆岩初次破断距不断增大，覆岩能量积聚快速增加，应变能密度也随之显著增加，在煤壁前方呈现明显激增现象。
2）工作面推进至同一位置时覆岩位移变化量随推进速度增加急剧增加，破坏范围广，呈现出明显的非连续变形特征，推进速度愈快，覆岩破断范围越广，运移演化愈剧烈。
3）随推进速度的提高工作面冲击危险性随之显著增加，在实体煤下回采时的冲击危险性指数远大于采空区下，变化幅度随推进速度提高明显增大，上覆W1145采空区与实体煤的交界区域为明显的高危险区域，工作面距离上覆W1145开切眼愈近冲击危险性越高。数值模拟研究结果综合认为W1123工作面推进速度应以不大于6.4 m/d。
-
![]()
图 1 越山气流类型与风速垂直廓线关系
Figure 1. Relationship between the type of cross-mountain airflow and the vertical profile of wind speed
![]()
图 5 外排土场稳定系数与高度关系
Figure 5. Relationship between stability factor of external dump with dump height
![]()
图 7 不同外排土场高度对气流以及水分质量分数影响模拟云图
Figure 7. Simulated cloud map of the influence of different external dump heights on airflow and moisture content
![]()
图 8 不同外排土场高度下气流速度与水分质量分数关系
Figure 8. Relationship between air velocity and moisture mass fraction at different external dump heights
![]()
图 9 外排土场高度与气流最大爬升高度关系
Figure 9. Relationship between height of external dump and maximum climbing height of air flow
![]()
图 10 不同外排土场角度对气流以及水分质量分数影响模拟云图
Figure 10. Simulated cloud map of the influence of different external dump angles on airflow and moisture content
![]()
图 11 不同外排土场角度下气流速度与水分质量分数关系
Figure 11. Relationship between air velocity and water mass fraction at different external dump angles
![]()
图 12 外排土场角度与气流最大爬升高度关系
Figure 12. Relationship between the angle of external dump site and the maximum climbing height of air flow
表 1 外排土场高度变化的边坡稳定系数
Table 1 Slope stability coefficient of external dumping Yard Height Variation
高度/m 角度/(°) 各计算方法下外排土场稳定系数 稳定系数均值 Ordinary Bishop Janbu 80 24 1.39 1.44 1.39 1.41 120 23 1.45 1.53 1.45 1.48 160 22 1.24 1.26 1.14 1.21 200 22 1.22 1.24 1.22 1.23 240 22 1.21 1.24 1.21 1.22 280 22 1.21 1.22 1.20 1.21 320 22 1.20 1.21 1.20 1.20 360 22 1.19 1.21 1.19 1.20 400 22 1.19 1.20 1.19 1.19 
下载: 导出CSV
-
[1] 钟水新. 地形对降水的影响机理及预报方法研究进展[J]. 高原气象,2020,39(5):1122−1132. doi: 10.7522/j.issn.1000-0534.2019.00083 ZHONG Shuixin. Research progress on the mechanism and forecast methods of the influence of topography on precipitation[J]. Plateau meteorological,2020,39(5):1122−1132. doi: 10.7522/j.issn.1000-0534.2019.00083
[2] COLLE Brian. Sensitivity of orographic precipitation to changing ambient conditions and terrain geometries: An idealized modeling perspective[J]. Journal of the Atmospheric Sciences,2004,61(5):588−606. doi: 10.1175/1520-0469(2004)061<0588:SOOPTC>2.0.CO;2
[3] 孟英杰,李丽平,王珊珊,等. 中尺度暴雨过程中地形抬升作用分析[J]. 安徽农业学,2010,38(12):6333−6336,6402. MENG Yingjie,LI Liping,WANG Shanshan,et al. Analysis of topographic uplift during mesoscale rainstorm[J]. Anhui Agriculture,2010,38(12):6333−6336,6402.
[4] 把 黎. 祁连山地区大气背景场及地形云降水过程的特征分析[D]. 兰州: 兰州大学, 2020: 39−44 BA Li. Characteristic analysis of atmospheric background field and orographic cloud precipitation process in Qilian Mountains[D]. Lanzhou: Lanzhou university, 2020: 39−44
[5] KATZFEY Jack. Simulation of extreme New Zealand precipitation events. part Ⅰ: sensitivity to orography and resolution[J]. Monthly Weather Review,1995,123(3):737−754. doi: 10.1175/1520-0493(1995)123<0737:SOENZP>2.0.CO;2
[6] 安小艳. 天山山区地形对降水空间分布的影响研究[D]. 石河子: 石河子大学, 2015: 31−38 AN Xiaoyan. Effect of topography on spatial distribution of precipitation in Tianshan Mountains [D]. Shihezi: Shihezi University, 2015: 31−38
[7] 廖 菲, 洪延超, 郑国光. 地形对降水的影响研究概述[J]. 气象科技, 2007, 35(3): 309−316. LIAO Fei, HONG Yanchao, ZHENG Guoguang. Research overview of the effect of topography on precipitation[J]. Meteorological Science and Technology, 2007, 35(3): 309−316.
[8] COLLE Brain,ZENG Yanguang. Bulk microphysical sensitivities within the MM5 for orographic precipitation. part II: Impact of Barrier Width and Freezing Level[J]. Monthly Weather Review,2005,132(12):2802−2815.
[9] 马英武. 准东矿区露天煤矿排土场水土保持防治措施分析[C]// 首届干旱半干旱区域水土保持生态保护论坛论文选编. 2013: 271−273. MA Yingwu. Analysis of soil and water conservation prevention and control measures of dump in Zhundong Mine area[C]// Selected papers of the first Forum on Soil and Water Conservation and Ecological Protection in arid and Semi-arid Regions. 2013: 271−273.
[10] 翟朝阳,陈高安,杨新峰,等. 微地形对大西沟野杏幼苗生境的气候和土壤温湿度的影响[J]. 中国农学通报,2019,35(27):114−120. ZHAI Chaoyang,CHEN Gaoan,YANG Xinfeng,et al. Effects of microtopography on climate and soil temperature and humidity of daxigou Apricot seedling habitat[J]. Chinese Agricultural Science Bulletin,2019,35(27):114−120.
[11] 孟 清,白红英,赵 婷,等. 秦岭山地气候变化的地形效应[J]. 山地学报,2020,38(2):180−189. MENG Qing,BAI Hongying,ZHAO Ting,et al. Terrain effects of climate change in the Qinling Mountains[J]. Journal of Mountain Science,2020,38(2):180−189.
[12] 黄菊梅, 邹用昌, 赵光平, 等. 中小尺度地形对降水影响研究[J]. 人民长江, 2013, 44(S1): 28-31. HUANG Jumei, ZOU Yongchang, ZHAO Guangping, et al. Impact of mesoscale topography on precipitation. [J] Yangtze River, 2013, 44(S1): 28-31.
[13] 文 迁,谭国良,罗嗣林. 降水分布受地形影响的分析[J]. 水文,1997(S1):64−66. WEN Qian,TAN Guoliang,LUO Silin. Analysis of precipitation distribution influenced by topography[J]. Hydrology,1997(S1):64−66.
[14] 孙继松. 气流的垂直分布对地形雨落区的影响[J]. 高原气象,2005,24(1):62−69. doi: 10.3321/j.issn:1000-0534.2005.01.010 SUN Jisong. Influence of Vertical Distribution of Air Flow on Topographic Rainfall Area[J]. Plateau Meteorology,2005,24(1):62−69. doi: 10.3321/j.issn:1000-0534.2005.01.010
[15] SCHNEIDER Linda,BARTHLOTT Christian,BARRETT Andrew,et al.Corinna Hoose The precipitation response to variable terrain forcing over low mountain ranges in different weather regimes[J]. Quarterly Journal of the Royal Meteorological Society,2018,144(713):970−989.
[16] 范广洲,吕世华. 地形对华北地区夏季降水影响的数值模拟研究[J]. 高原气象,1999(4):659−667. doi: 10.3321/j.issn:1000-0534.1999.04.023 FAN Guangzhou,Lv Shihua. Numerical simulation of the effect of Topography on summer precipitation in North China[J]. Plateau Meteorology,1999(4):659−667. doi: 10.3321/j.issn:1000-0534.1999.04.023
[17] 郭 璇,田 苗,洪程之,于瑞波,陈权亮. 川西高原地形特征对降水的影响[J]. 新疆农垦科技,2017,40(12):40−42. doi: 10.3969/j.issn.1001-361X.2017.12.017 GUO Xuan,TIAN Miao,HONG Chengzhi,et al. Effects of topographic characteristics on precipitation in western sichuan plateau[J]. Xinjiang agricultural reclamation technology,2017,40(12):40−42. doi: 10.3969/j.issn.1001-361X.2017.12.017
[18] 陈 印,杨 溢,刘 磊,等. 排土场边坡稳定性分析及台阶高度的确定[J]. 中国锰业,2016,34(2):21−24. CHEN Yin,YANG Yi,LIU Lei,et al. Stability analysis of dump slope and determination of step Height[J]. China Manganese Industry,2016,34(2):21−24.
[19] 高登义. 中国山地环境气象学[M]. 河南科技出版社, 2005: 90−92. [20] 于晓晶,赵 勇. 地形对天山夏季降水影响的模拟[J]. 中国沙漠,2016,36(4):1133−1143. YU Xiaojing,Zhao Yong. Simulation of the effect of topography on summer precipitation in tianshan Mountains[J]. The Chinese desert,2016,36(4):1133−1143.
[21] 王沛霖,许丽章. 抬升凝结高度的精确计算[J]. 中山大学学报论丛,1993(1):146−149. WANG Peilin,Xu Lizhang. Precise calculation of uplift condensation height[J]. Journal of Sun Yat-sen University,1993(1):146−149.
[22] 陈创买, 郭英琼. 气象常用参数和物理量查算表 [M]. 北京: 气象出版社, 1980. 3−4. CHEN Chuangmai, Guo Yingqiong. Checking table of Common Meteorological Parameters and Physical Quantities [J]. Beijing: China Meteorological Press, 1980. 3-4.
[23] 张兆顺, 崔桂香, 许春晓. 湍流大涡数值模拟的理论和应用[M]. 北京: 清华大学出版社, 2008: 57-60. [24] 张兆顺, 崔桂香, 许春晓. 湍流理论与模拟[M]. 北京: 清华大学出版社, 2017: 236-248. [25] 唐家鹏. ANSYS FLUENT 16.0超级学习手册[M]. 北京: 人民邮电出版社, 2016: 395-397. [26] 凌桂龙. Fluent 19.0流体计算从入门到精通[M]. 北京: 电子工业出版社, 2015: 221-222. -
期刊类型引用(5)
1. 赵兵朝，赵阳，冯杰，王京滨，冯欣怡. 采动地表与输电线塔大板基础协同变形规律研究. 采矿与安全工程学报. 2024(03): 472-480 . 
百度学术
2. 张建，亓佳利，李长青. 厚煤层开采上覆岩层移动与应力分布特征研究. 山东煤炭科技. 2024(06): 139-143+154 . 百度学术
3. 李文福，马小辉，刘咸富，王向阳，何伟，王岩，白晓敏，刘万成. 冲击性煤层综放面采空区漏风区域识别及控制方法研究. 煤炭技术. 2024(09): 133-137 . 百度学术
4. 彭宝山，王永乐，杨学孟. 特厚弱胶结顶板煤层综放开采覆岩破坏特征与强矿压机理. 煤炭技术. 2024(10): 75-80 . 百度学术
5. 李志，杨栋，贾毅超. 深厚煤层开采过程中覆岩运移演变机理及矿压管控研究. 矿业研究与开发. 2024(12): 68-75 . 百度学术
其他类型引用(1)
计量
- 文章访问数: 118
- HTML全文浏览量: 23
- PDF下载量: 30
- 被引次数: 6
