随着浅部资源的逐渐枯竭, 矿产资源开采向深部推进已成为未来发展趋势[1].全尾砂充填采矿法不仅可以有效控制深部开采特有的安全难题 (如控制岩爆、降低工作面温度、保障作业安全), 还可以解决地面尾矿库建设占用大量土地, 恶化人类赖以生存的自然环境的难题, 因此, 绝大部分深部开采矿山均首选全尾砂充填采矿法[2-3].
目前深部开采相关充填理论体系尚未形成, 深部开采充填体的强度设计均类比浅部充填理论.但是, 深部开采技术条件与浅部开采差异大, 按照浅部充填理论选择充填体强度指标, 往往造成不必要的水泥浪费或者强度达不到工程实际要求.针对上述问题, 国内外学者进行了大量研究工作:Berndt, Fall和Belem等从力学角度对充填体的破坏特征进行了大量研究[4-6]; Liu等基于损伤力学理论研究了深部开采高阶段尾砂充填体破坏的能量耗损规律及强度的非线性优化设计[7]; 马斐等基于变形量控制对充填体力学参数进行了研究[8]; 周士霖根据深部开采充填体的不同破坏形式建立了剪切控制型和强度控制型破坏模型[9].但上述研究均将充填体视为线弹性介质, 较难真实反映充填体破坏前的实际力学特性, 因为充填体是一种多相复合非线性介质;而且相对于浅部开采, 深部开采围岩应力高, 聚集的应变能大, 充填体的主要作用是控制和减轻应变能瞬间发生释放 (岩爆), 这也从另一方面对充填体的质量和强度设计提出了更高的要求.
本文根据某矿全尾砂充填体单轴压缩实验的应力-应变曲线在峰值应力前呈三次多项式函数变化趋势的特征, 建立了5种不同配比全尾砂充填体峰值应力前的非线性本构模型.在此基础上, 根据深部开采过程中全尾砂充填体与围岩耦合作用下能量损耗相近的原则, 通过线性回归分析, 得出了深部开采全尾砂充填体强度指标设计公式, 为确定深部开采全尾砂充填体强度指标提供了理论依据.
1 全尾砂充填体力学实验与分析某矿全尾砂, 粒径组成 (粒径/mm, 质量分数/%) 为:>0.5, 1.7;0.5~0.25, 22.6;0.25~0.075, 39.1;0.075~0.05, 12.3;0.05~0.005, 21.3; < 0.005, 3.采用全尾砂与32.5普通硅酸盐水泥制备质量分数为68%, 灰砂比分别为1:4, 1:6, 1:8, 1:10和1:12的5组试块.20 ℃标准恒温养护箱养护28 d后, 用MTS Landmark实验机对试块进行单轴压缩实验.5种不同配比全尾砂充填体的单轴抗压应力-应变曲线如图 1中实线所示.
由图 1可知, 峰值应力前, 全尾砂充填体变形主要分为三个阶段:
1) 加载初期, 应力-应变曲线呈明显的下凹形 (如图 1中曲线1的OA段), 曲线斜率随着应力的增加而快速增大, 这主要是因为充填体内的微孔隙、微裂纹被逐渐压密, 充填体抵抗变形的能力快速增强, 充填体弹性模量快速增大.充填配比越高, 压密过程越明显.
2) 随着加载的增大, 应力-应变曲线呈微微下凹形 (如图 1中曲线1的AB段), 曲线斜率随着应变的增加而稍有增加, 表明充填体进入弹性阶段, 充填体弹性模量增加不大.
3) 加载继续增大, 应力-应变曲线呈明显的上凸形 (如图 1中曲线1的BC段), 曲线斜率随着应力的增加而逐渐减小至零, 表明充填体进入屈服阶段, 充填体内原有微裂纹开始扩展, 并产生新裂纹, 充填体抵抗变形的能力逐渐减小, 充填体弹性模量逐渐减小.充填配比越高, 屈服过程越明显.
由上述分析可知, 峰值应力前, 随着应力的增大, 全尾砂充填体应力-应变曲线由明显的下凹形逐渐变为上凸形, 应力-应变关系并不呈线性正比关系, 而是呈明显的三次多项式函数变化趋势.
2 全尾砂充填体非线性本构模型由上述实验与分析可知, 全尾砂充填体在峰值应力σM(即抗压强度) 前应力σ和应变ε之间存在良好的三次多项式函数关系, 设其表达式为
(1) |
式中a, b和c为常数.
根据实验得到的全尾砂充填体应力-应变曲线及其边界条件有
(2) |
式中:εM为相应峰值应力的应变; εS为相应屈服应力的应变 (曲线拐点处).
将式 (1) 代入式 (2) 可得出有关a, b和c的3个方程, 解方程组得到:
(3) |
根据实验数据 (图 1中实线), 可得到不同配比全尾砂充填体的峰值强度、峰值应变和屈服应变参数, 将其分别代入式 (3) 和式 (1), 可得到不同配比全尾砂充填体应力-应变的非线性本构方程 (见表 1).
根据尾砂充填体应力-应变的非线性本构方程得出的不同配比全尾砂充填体峰值应力前的应力-应变的计算曲线 (图 1中虚线) 与实验测得曲线吻合.
3 深井全尾砂充填体强度指标深部矿石开采前, 矿岩处于原始平衡状态, 矿石开采形成采空区, 围岩卸载并释放能量; 全尾砂充填体充入采空区后受到围岩变形的压缩而产生抗力, 抑制围岩变形, 达到新的平衡状态, 同时蓄积变形能.也就是说,深部开采过程中, 全尾砂充填体与围岩在耦合作用下存在能量交换, 可以从能量角度分析全尾砂充填体的强度指标.
3.1 全尾砂充填体蓄能分析考虑到深部开采过程中, 全尾砂充填体至少暴露一侧, 分析最危险情况并对模型加以简化, 即考虑一维情况.单轴压缩下, 单位体积全尾砂充填体达到峰值应力时蓄积的能量 (即峰值比能VF) 可按下式计算[10]:
(4) |
根据实验结果和式 (4), 计算出配比分别为1:4, 1:6, 1:8, 1:10和1:12时, 全尾砂充填体峰值比能 (MJ/m3) 分别为0.034 2, 0.022 5, 0.015 3, 0.008 3和0.004 9.可见, 全尾砂充填体配比越高, 峰值比能越大.
3.2 全尾砂充填体强度指标对不同配比全尾砂充填体抗压强度VF和峰值比能σM进行线性回归分析, 如图 2所示.
VF和σM回归方程如下:
(5) |
经回归检验, 式 (5) 的相关系数为95.03%, 满足工程精度要求[7, 9].
由于深部围岩的刚度远大于全尾砂充填体的刚度, 围岩可采用线弹性模型.文献[10]中详细给出了围岩比能VR(矿石开采后单位体积围岩释放的能量) 的计算过程 (不考虑围岩的构造应力), 即
(6) |
式中:σ0为原岩应力; E0为围岩弹性模量; H为上覆岩层厚度; γ为上覆岩层容重.
根据深部开采过程中全尾砂充填体与岩体耦合作用下能量损耗相近的原则, 并结合式 (4) 和式 (6), 全尾砂充填体的强度指标设计公式为
(7) |
式中K为安全系数.
式 (7) 表明深部开采全尾砂充填体强度指标除与充填体本身的力学性能有关外, 还与围岩原岩应力和弹性模量有关.围岩原岩应力越大、弹性模量越小, 所需的全尾砂充填体强度越高.
不考虑围岩的构造应力, 则式 (7) 可写成:
(8) |
设围岩容重γ=27 kN/m3, 围岩弹性模量E0分别为15, 20, 25和30 GPa, 安全系数K取1.05, 由式 (8) 计算的不同开采深度所要求的全尾砂充填体强度指标如图 3所示.
从图 3可看出, 围岩弹性模量越低或埋藏越深, 深部开采所要求的全尾砂充填体强度越高.
4 实例运用某深部开采矿山最大开采深度为1 050 m, 设计采用全尾砂阶段空场嗣后充填采矿法, 围岩容重27 kN/m3、弹性模量28 GPa, 安全系数取1.02.按照线性本构模型确定的全尾砂充填体强度指标为1.56 MPa[7], 根据式 (8) 计算的强度指标为1.45 MPa.
矿山生产实践表明, 按照全尾砂充填体非线性本构模型确定的全尾砂充填体强度指标1.45 MPa是合理的, 实现了矿山的安全低成本开采[11].
5 结论1) 全尾砂充填体峰值应力前的应力-应变关系并不呈线性正比关系, 而是呈非线性的三次多项式函数变化趋势.
2) 全尾砂充填体配比越高, 峰值比能越大, 充填体承载力越强.
3) 深部开采全尾砂充填体强度指标除与充填体本身的力学性能有关外, 还与围岩原岩应力和弹性模量有关.围岩原岩应力越大、弹性模量越小, 所需的全尾砂充填体强度越高.
4) 利用全尾砂充填体非线性本构模型确定的某深部开采矿山所需的全尾砂充填体强度指标为1.45 MPa, 实现了矿山的安全低成本开采, 为矿山全尾砂充填体强度设计提供了一定理论依据.
[1] |
董陇军, 李夕兵, 李萍萍.
深井开采灾害应对决策技术综述[J]. 中国有色金属学报, 2009, 61(1): 116–119.
( Dong Long-jun, Li Xi-bing, Li Ping-ping. Crucial understanding and decision-making technology for deep mining[J]. The Chinese Journal of Nonferrous Metals, 2009, 61(1): 116–119. ) |
[2] | Sheshpari M. A review of underground mine backfilling methods with emphasis on cemented paste backfill[J]. Electronic Journal of Geotechnical Engineering, 2015, 20(13): 5183–5208. |
[3] | Sheshpari M. A review on influential parameters on the strength of cemented paste backfill:part Ⅰ[J]. Electronic Journal of Geotechnical Engineering, 2015, 20(22): 12443–12465. |
[4] | Berndt C C, Rankine K J, Sivakugan N N. Materials properties of barricade bricks for mining applications[J]. Geotechnical and Geological Engineering, 2007, 25(4): 449–471. DOI:10.1007/s10706-007-9122-y |
[5] | Fall M, Benzaazoua M, Saa E G. Mix proportioning of underground cemented tailings backfill[J]. Tunnelling and Underground Space Technology, 2008, 23(1): 80–90. DOI:10.1016/j.tust.2006.08.005 |
[6] | Belem T, Benzaazoua M. Design and application of underground mine paste backfill technology[J]. Geotechnical and Geological Engineering, 2008, 26(2): 147–174. DOI:10.1007/s10706-007-9154-3 |
[7] | Liu Z X, Lan M, Xiao S Y, et al. Damage failure of cemented backfill and its reasonable match with rock mass[J]. Transactions of Nonferrous Metals Society of China, 2015, 25(3): 954–959. DOI:10.1016/S1003-6326(15)63684-6 |
[8] |
马斐, 张东升, 张晓春.
基于变形量控制的充填体力学参数研究[J]. 岩土力学, 2007, 28(sup1): 545–548.
( Ma Fei, Zhang Dong-sheng, Zhang Xiao-chun. Research on properties of the filled body based on the deformation control[J]. Rock and Soil Mechanics, 2007, 28(sup1): 545–548. ) |
[9] |
周士霖. 深部开采充填体稳定性及与岩体智能匹配研究[D]. 长沙: 中南大学, 2012.
( Zhou Shi-lin.Research on stability of filling and intelligent matching between filling and rock mass in deep mining [D].Changsha:Central South University, 2012. ) |
[10] | Du K, Li X B, Li D Y, et al. Failure properties of rocks in true triaxial unloading compressive test[J]. Transactions of Nonferrous Metals Society of China, 2015, 25(2): 571–581. DOI:10.1016/S1003-6326(15)63639-1 |
[11] |
鄯善联合矿业有限公司. 新疆鄯善县细沙沟铁矿充填采矿关键技术综合研究[R]. 长沙: 中南大学, 2014.
( Shanshan Mining Corporation Limited.Comprehensive research on key technology of filling mining in Xinjiang Shanshan Xishagou Iron Mine [R].Changsha:Central South University, 2014. ) |