不规则大采空区不稳定区域研究

唐敏康， 刘明荣， 张东炜 ， 赵 奎

(江西理工大学资源与环境工程学院，江西 赣州 341000)

不规则大采空区不稳定区域研究

唐敏康， 刘明荣， 张东炜 ， 赵 奎

(江西理工大学资源与环境工程学院，江西 赣州 341000)

某金银矿经过多年的地下开采，已形成暴露面积约为0.8万m2，体积约为9.5万m3不规则大采空区.为消除大规模地压隐患，需要确定不稳定区域以便进行地压处理.数值模拟研究结果表明，94 m中段顶板存在较大的拉应力及下沉位移量最大，初步判断此中段为不稳定区域.在此基础上，通过相似模拟试验发现94 m中段顶板位置首先出现裂纹，并且最先垮塌，进一步论证了此中段为不稳定区域，为地压处理方案提供了科学依据.

数值模拟；相似模拟；采空区；不稳定区域

0 引 言

空场法开采目前在我国仍广泛应用，空场法开采将遗留下采空区，采空区地压问题是矿山开采中重要安全问题[1－2].采空区处理是消除大规模地压隐患，保障开采安全的重要手段.采空区处理方案的确定需要对采空区稳定性作出评价，尤其是要确定关键不稳定性区域，以便制定合理的采空区处理方案[3－4].

某金银矿经过多年开采，目前已形成暴露面积约为0.8万m2，体积约为9.5万m3不规则大采空区.为了保护矿产资源，避免地压灾害，确保此采空区周围的安全作业.有必要确定采空区的不稳定区域，并进行地压处理.

确定不规则大采空区的不稳定区域，现场监测虽然真实、可靠，但周期长、工作量大、成本也高；理论分析与经验公式计算简单易行，但很多复杂的地质特征被忽略或简化，解决工程问题误差较大[5].数值模拟与相似模拟结合应用则可以弥补上述不足 结合某金银矿工程实例，采用数值模拟和相似模拟试验相结合的研究方法，首先对不规则大采空区进行了数值模拟研究，得到采空区围岩位移及应力分布情况，在此基础上，进行了相似模拟试验，分析了采空区围岩变形、破坏规律，最后确定此采空区的不稳定区域，为地压处理方案提供了科学依据.

1 工程背景

金银矿目前已形成高度约80 m，顶板斜长约100 m的不规则大采空区，此采空区标高从50～166 m，贯通 154 m、132 m、111 m、94 m、78 m 等中段. 顶板跨度及暴露面积比较大，影响空区周围安全作业.考虑到安全兼顾经济利益，需确定不稳定区域并进行处理.矿体围岩主要是中细粒花岗岩.围岩与矿体力学参数见表1.

表1 岩石力学参数值

2 数值模拟

根据矿山实际资料，建立了8节点四边形单元二维有限元模型.模型长 200 m，高60 m，节点19642个，单元 28571个.在模型左、右两边施加水平约束，底边施加垂直方向约束.

有限元计算结果及分析.先施加覆土层自重载荷，然后开挖.此类不规则大空区顶板破坏主要形式是拉应力破坏.X方向（水平方向）的应力云图与位移矢量图见图1、图2.

图1 X方向应力云图

图2 位移矢量图

从X方向的应力云图可知空区的顶板及底板都存在相对较大的拉应力，根据位移矢量图可以清楚显示94 m中段的顶板下沉位移最大.底板有向上的位移，推断底板可能出现微量的底鼓.根据空区顶板处的拉应力相对较大及位移量最大，是最先出现破坏的位置，可以初步确定94 m中段位置是此空区的不稳定区域.

3 相似模拟试验

有限元数值模拟适用于弹塑性介质，无法模拟顶板变形、开裂、垮塌破坏发展全过程[8－9].因此，借助相似模拟试验[10]可以准确判断开挖过程中顶板最先冒落的位置.这对不稳定区域的确定具有参考依据.

3.1 试验装置及测量手段

试验采用ZYDL-YS120/100岩体平面相似模拟试验系统，模型尺寸：长×宽×高=2000 mm×300 mm×1200 mm.该系统能控制竖向均布载荷和水平载荷，其中最大竖向载荷120 kN（最大应力0.2 MPa），最大水平载荷100 kN（最大应力为0.25 kN），变形范围5～80 mm.模型表面贴有BE120-5AA胶基应变片，通过BZ2205C静态电阻应变仪进行变形测量.

3.2 模型设计

采空区的平面应力模拟，需满足几何相似、容重相似、应力相似，具体选取如下：

（1）选取几何相似比为 1∶300，即：∂L＝LH／LM＝300；

其中 LH为实际尺寸；LM为模型尺寸.

根据现场空区大小与几何相似比，确定了模型尺寸.模型尺寸见图3.

（2）由于岩石围岩的平均容重为2.6 t/m3，而模型材料的容重大约为1.7 t/m3，故容重相似比为：∂γ=γH/γM=1.53；

其中 γH为实际容重 ；γM为模型材料容重.

图3 模型图

（3）应力相似比可由几何相似比和容重相似比计算得到，其具体表达式为：∂σ=σH／σM＝∂L·∂γ=300×1.53=459；其中σH为现场岩层应力；σM为模拟岩层应力.

根据应力相似比计算实际值与模型应力值见表2.

表2 实际值与模型应力值

相似材料由分级尾砂和水泥组成.试验前对不同浓度不同灰沙比的试样进行了单轴抗压强度试验.最终选定重量比为1∶25的水泥与尾砂作为相似模拟材料，其抗压强度为0.142 MPa.

为了解模型在开挖过程中空区的变形规律，本试验采用应变片进行变形测量.根据数值模拟位移云图，选择最佳应变测量位置.应变片布置及对应中段图，见图4.

图4 应变片位置及对应中段图

3.3 试验过程

试验过程中加压系统加载32.5 kN的压力，相当于原型的覆岩层的垂直应力，稳定压力后进行开挖.开挖过程中应变仪每分钟采集一次数据.空区完全形成后，模型无明显的破坏，继续加载至77.79 kN，空区垮塌，在此过程中了解模型破坏的规律.

3.4 试验结果及分析

模型空区形成后继续加载过程中，12#应变片位置（即94 m中段顶板）首先出现裂纹，见图5.裂纹宽度1 mm左右，长度8.7 cm左右.继续加载，25#、26#应变片位置（即132 m中段至154 m中段中间位置）出现宽度1.2 mm左右，长度约10 mm的裂纹.随后在不同位置的表面及空区内侧出现裂纹，继续加载，15#应变片位置的空区内侧（即94 m中段顶板）首次出现小块体冒落.紧接其它位置也出现块体冒落，直至完全压垮，见图6.

图5 12＃位置处裂纹

图6 模型压垮

对应变值进行统计整理，分中段把应变片在相同方向的应变作为纵坐标，对应的载荷值作为横坐标，作在同一张图上，从而相互对比确定不同中段的相对变形.部分中段载荷应变图，见图7～图12.

94 m中段在垂直方向上9#、11#测点由受压变为受拉.由此看出开挖完成后9#测点处存在微量底鼓，11#测点顶板处有微量下沉.13#测点由受拉变为受压，继续加载，又变为受拉.在水平方向上16#测点开始处于受拉状态，而后一直处于受压状态，应变基本不变.14#和18#测点均一直处于受拉状态，变化趋势相近，应变开始一直增大，到达峰值后逐渐减小.由此可见水平方向也向空区变形，测点在加载后期均有进入塑性破坏的趋势.

图7 94 m中段垂直方向载荷应变图

图8 94 m中段水平方向载荷应变图

图9 111 m中段垂直方向载荷应变图

图10 111 m中段水平方向载荷应变图

图11 132 m中段垂直方向载荷应变图

图12 132 m中段水平方向载荷应变图

111 m中段在垂直方向上，只有在15#测点即该中段底板附近应变比较大外，在其他测点应变微小，19#测点在空区形成阶段有明显受压增大趋势，但在后期增加荷载期间受压减小直至变为受拉，17#测点变形相对比较小.在水平方向，18#测点应变较大，后期进入塑性破坏，16#和20#测点在开挖阶段应变有明显的增大趋势，但是随后实验中均有减小的趋势，其中16#一直受压，20#一直受拉.

132 m中段在垂直方向上应力变化趋势一致，由受压变为受拉，在开挖阶段均有受拉增大趋势，但后来又随着荷载逐渐增大，19#和23#应变稳定，25#测点变形一直增大.在水平方向变化趋势一致，除在开挖阶段有明显应变增大趋势，后阶段都有减小趋势，但20#一直处于受拉状态，24#和26#由受拉变为受压状态，其中24#测点变形较大.

4 结 论

采用数值模拟与相似模拟2种方法来确定某金银矿采空区不稳定区域.技术上可优势互补，结果上可相互印证，不失为一种理想的分析手段.通过2种方法分析得出了以下主要结论：

（1）数值模拟结果显示：94 m中段顶板位置位移量最大并存在较大的拉应力.相似模拟试验结果显示：该区域最先出现裂纹及块体冒落，并最先垮塌.由此可以确定该区域为不稳定区域.地压处理时，应优先考虑此区域.

（2）开挖过程中围岩不仅在竖直方向往空区变形，水平方向也向空区变形，底板存在微量底鼓，这都与现场实际相吻合.不规则采空区同一位置开挖过程中受力状态不同，不同位置的受力状态也不同.

[1]周宗红，任凤玉，袁国强.桃冲铁矿采空区处理方法研究[J].中国矿业，2005，14（12）:15-16.

[2]李俊平.缓倾斜采空区处理新方法及采场地压控制研究[D].北京：北京理工大学，2003.

[3]解世俊.金属矿床地下开采[M].北京：冶金工业出版社，2008.

[4]张军胜，李俊平，宋士生.采空区处理新方法研究与应用[J].金属矿山，2009（8）：19-22.

[5]李向阳，李俊平，周创兵，等.采空场覆岩变形数值模拟与相似模拟比较研究[J].岩土力学，2005，26（12）：1907-1912.

[6]蒋树屏，刘洪洲，鲜学福.大跨度扁坦隧道动态施工的相似模拟与数据分析研究[J].岩石力学与工程学报，2000，19（5）：567-572.

[7]邵小平，石平五，张嘉凡.亚急倾斜煤层采空区覆岩变形破坏数值与相似模拟[J].西安科技大学学报，2008，28（2）：249-253.

[8]李俊平，卢忠瑜.ANSYS在采空区稳定性安全评价中的应用[J].中国钼业，2006，30（5）：13-17.

[9]韩贵雷，韩立军，王延宁，等.厚层复合石灰岩顶板垮塌特性试验研究[J].采矿与安全工程学报，2007，24（1）：60-64.

[10]张明建，郜进海，魏世义，等.倾斜岩层平巷围岩破坏特征的相似模拟试验研究[J].岩石力学与工程学报，2010，29（增 1）：3259-3264.

On the unstable region of the large-sized irregular mine-out areas

TANG Min-kang,LIU Ming-rong,ZHANG Dong-wei,ZHAO Kui

(School of Resource and Environmental Engineering,Jiangxi University of Science and Technology,Ganzhou 341000,China)

The big irregular mined-out area was formed after years of underground mining in a gold-silver mine.The exposed area is approximately 0.8 million square meters with the volume of 95,000 cubic meters.To remove the hidden danger of large-scaled ground pressure,the instability district must be determined for processing ground pressure.The numerical simulation results show that the 94 m phase roof exist big tensile stress and the submersion displacement quantity is biggest,this phase roof is preliminary judged as the instability district.On these bases,the similarity simulation experiment discovered that the 94 m phase roof position first presents the crack,and collapses first.This phase roof is further proved as the instability district.It provides scientific basis for processing ground pressure.

numerical simulation;similarity simulation;mined-out areas;instability district

TD326

A

1674-9669（2011）06-0043-04

2011-09-12

唐敏康（1956－ ），男，教授，主要从事静电除尘器理论及应用研究，E－mail:tmkjxust@126.com.