创刊于1987年, 双月刊
主管:

江西理工大学

主办:

江西理工大学
江西省有色金属学会

ISSN:1674-9669
CN:36-1311/TF
CODEN YJKYA9

基于断层影响下的充填体强度需求优化研究

李广波, 盛宇航, 吴再海, 宋泽普

李广波, 盛宇航, 吴再海, 宋泽普. 基于断层影响下的充填体强度需求优化研究[J]. 有色金属科学与工程, 2023, 14(5): 684-691. DOI: 10.13264/j.cnki.ysjskx.2023.05.011
引用本文: 李广波, 盛宇航, 吴再海, 宋泽普. 基于断层影响下的充填体强度需求优化研究[J]. 有色金属科学与工程, 2023, 14(5): 684-691. DOI: 10.13264/j.cnki.ysjskx.2023.05.011
LI Guangbo, SHENG Yuhang, WU Zaihai, SONG Zepu. Research on strength requirement optimization of infill bodies under the influence of faults[J]. Nonferrous Metals Science and Engineering, 2023, 14(5): 684-691. DOI: 10.13264/j.cnki.ysjskx.2023.05.011
Citation: LI Guangbo, SHENG Yuhang, WU Zaihai, SONG Zepu. Research on strength requirement optimization of infill bodies under the influence of faults[J]. Nonferrous Metals Science and Engineering, 2023, 14(5): 684-691. DOI: 10.13264/j.cnki.ysjskx.2023.05.011

基于断层影响下的充填体强度需求优化研究

基金项目: 

山东省重大科技创新工程项目 2019SDZY05

国家“十三五”重点研发计划项目 2018YFC0604600

详细信息
    通讯作者:

    李广波(1989— ),工程师,主要从事矿山充填工艺技术研究。E-mail:985315860@qq.com

  • 中图分类号: TD853.34

Research on strength requirement optimization of infill bodies under the influence of faults

  • 摘要: 为分析断层影响下的充填体强度需求问题,以某矿的断层赋存特点为工程背景,采用力学推导的方法,重点对临近断层处两步采进行了力学分析,提出了充填体强度需求解析公式。同时,采用FLAC3D对断层处的充填体与矿体稳固性进行数值模拟研究,得到了充填体强度需求模拟解。对比分析了2种计算结果,确定了合理的充填体强度需求,保证了断层处矿体安全高效开采。本研究可为同类矿山充填体强度需求优化计算提供参考。
    Abstract: To analyze the strength requirements of infill bodies under the influence of faults, with the fault occurrence characteristics of a certain mine as the engineering background, the mechanics derivation method was used to focus on the mechanical analysis of the two-step mining near the fault, and the calculation formula for the strength requirements of infill bodies was proposed. Meantime, the FLAC3D software was used to conduct numerical simulation research on the stability of the infill bodies and ore body at the fault, and the simulation solution of the strength requirements of the infill bodies was obtained. By comparing and analyzing the two calculation results, the reasonable strength requirements of the infill bodies were determined, ensuring the safety and high-efficiency mining of the ore body at the fault. This study can provide reference for the calculation of infill body strength requirement optimization in similar mines.
  • 随着金属矿山逐渐进入深井开采,充填采矿法得到广泛应用,充填体与矿体的变形破坏过程变得越来越复杂,研究其相互作用具有重要意义,特别是开采至断层处存在应力释放的风险[1-4]。当工作面向断层推进时,断层的构造应力与工作面超前支承压力叠加,形成的高应力集中作用于两侧的充填体及矿体上[5-6]。目前,研究者对充填体与矿体相互作用、断层影响下的围岩失稳机理等开展了大量研究。吴爱祥等[7]利用Belem模型分析计算了充填体3种应力状态下的目标强度,对比拱架效应引入前后充填体垂直应力降低比率。刘光生等[8-9]基于实际采充时序过程,考虑充填体与矿体接触作用,提出了充填体强度要求修正解析方法,并采用FLAC3D进行验证。还有研究者采用3种力学理论模型进行了分层胶结充填体强度需求解析计算,得到了胶结充填体强度与充填次数之间的数学模型及折减系数[10-11];利用数值模拟、物理模型、解析公式和现场监测等方法对充填体成拱作用和应力分布进行研究[12-15]。目前的研究大多只考虑单一因素,并没有考虑断层对于充填体强度需求的影响,使得采场回采至上盘断层附近,因断层不稳固极易导致塌落[16-20]。现场施工时,经常出现2种情况,一是为保证安全偏离上盘断层较远,导致矿体损失;二是沿上盘断层施工时,上盘暴露面积过大,矿柱位置选择不合理,导致上盘塌落,矿体难以采出。

    本文针对以上不足,以山东某矿的断层赋存特点为研究对象,首先根据该矿采场尺寸和充填体与矿体的物理力学参数,开展断层处充填体与矿体的两步采力学分析,确定出合理的充填体强度需求解析解;再利用FLAC3D建立采场三维数值模型,对两步采过程中断层处的充填体与矿体稳固性进行模拟研究,得到采充过程中两步采采场充填体强度需求数值解,对比研究充填体强度需求解析解、数值解和实际生产值,校验修正解析模型的可靠性,保证断层处矿体安全高效开采。

    该矿原设计的充填体强度及矿柱尺寸,主要是利用同类矿山经验类比选取的强度指标,导致矿山存在充填强度设计不合理、充填成本居高不下等问题。为了有效解决矿山面临的生产难题,本文选择该矿100线以北、116线以南矿体进行断层影响下的充填体强度需求研究。

    该矿床属于第四系松散沉积物广泛分布,主断层倾角55°,落差2 m,矿体上盘主要为绢英岩,下盘主要为闪长岩,在断裂面上发育有一层厚约0.2 m的黑色断层泥。矿体上盘与围岩为断层接触关系,界线明显;而矿体下盘与围岩呈渐变过渡关系,无明显的边界线。断层两侧矿体受其影响,强度明显降低,破碎程度增大,并有一定的似层状特点,施工至断层处时极易发生冒顶、片帮等安全事故,所以亟需解决充填体强度需求问题。

    由于矿区周边临近农田和村庄,同时结合矿体赋存条件,经综合考虑采用上向水平分层充填采矿法,按隔一采一顺序沿矿体走向先采矿房,出矿后胶结充填,再按顺序依次进行相邻两侧两步采采场回采,如图1所示。采场埋深600 m,走向长度120 m,倾向长度为80~100 m。

    图  1  盘区式上向水平分层充填法
    1.斜坡道;2.脉外溜井;3.分段巷道;4.分段联络巷道;5.脉内溜矿井;6.人行井;7.充填通风井;8.阶段运输平巷。
    Figure  1.  Panel upward cut and fill method

    目前,对充填体力学特性的研究逐渐增多,但国内矿山仍然以经验类比进行充填强度计算研究,导致国内同类型矿山的充填水泥单耗远大于国外。针对该矿地质赋存条件,矿体整体较为破碎,特别是临近断层时两侧矿体受断层影响,强度明显降低,破碎程度增大,极易发生冒顶、片帮等安全事故。

    本文进行充填体强度要求解析计算时,重点研究两步采回采过程中断层处力学模型。在两步采采场回采过程中,断层与相邻矿体、充填体并不是孤立存在的,受上盘岩体和断层两侧胶结充填体的夹制,以此进行两步采采场揭露前后胶结充填体的应力计算,如图2所示。

    图  2  两步采回采时断层力学分析
    Figure  2.  Mechanical analysis of fault in two step mining

    考虑充填体与矿体接触作用以及断层受到一步采采场充填体的水平挤压力,断层两侧壁接触面产生相应摩擦力,对断层进行应力平衡分析,得到:

    W上盘+W=2f1sinϕ (1)

    式(1)中:W上盘为上盘岩体作用于两步采断层处的挤压力,单位kN;W为断层自身重力,单位kN;f1为断层受到两侧壁的一步采采场胶结充填体水平挤压力产生的摩擦力,单位kN;ϕ为断层的倾角,单位(°)。

    上盘矿体的重力W上盘为:

    W上盘=γ1BhH/tanϕ (2)

    式(2)中:γ1是上盘岩体的容重;h为矿体埋深高度;BH为矿房的宽度和高度。由于上盘岩层较厚大,只考虑其投影等效面积。

    断层的局部自重为:

    W=γ2BhH/sinϕ (3)

    式(3)中:γ2是断层的容重;h为断层的厚度。

    断层受到两侧壁的一步采采场胶结充填体摩擦力f1计算如下:

    f1=0HcsH-htanϕdh (4)

    式(4)中:cs为接触面内聚力;H为矿房的高度。

    假设采场岩壁十分粗糙,接触面内聚力cs与充填体内聚力c相等[21-23],则胶结充填体内聚力c为:

    c=cs=1H-hγ1Bh2sinφ+γ2Bhsin2φ (5)

    假设充填体服从直线型摩尔库仑准则,则充填体内聚力c与单轴抗压强度σc的关系满足:

    σc=2ccosϕc1-sinϕc (6)

    式(6)中:ϕc为一步采胶结充填体的内摩擦角,单位(°)。

    将式(5)代入式(6)后,计算得出充填体单轴抗压强度σc

    σc=2cosϕc1-sinϕcH-hγ1Bh2sinϕ+γ2Bhsin2ϕ (7)

    将矿山采场尺寸及充填体参数γ1=26.87 kN/m3γ2=15.47 kN/m3B=4 m、H=5 m、ϕ=55°、h=0.2 m、h= 600 m、ϕc=30°代入式(7),可得极限平衡状态两步采采场回采时充填体强度需求解析解σc=2.85 MPa。

    为进一步对比采场两步采过程中断层处充填体强度需求,需要利用FLAC3D建立采场三维数值模型,对两步采过程中断层处的充填体与矿体稳固性进行数值模拟研究,得到采场两步采过程中充填体强度需求数值解,将充填体强度需求解析解和数值解进行对比研究。

    该矿经过岩石分级评价、钻探取芯和室内力学实验等后,得到各岩体的物理力学参数,如表1所列;同时,利用该矿的分级尾砂和胶凝材料开展充填体配比实验(质量浓度为68%),得到不同灰砂比充填体物理力学参数,如表2所列。

    表  1  各岩体物理力学参数表
    Table  1.  Table of physical and mechanical parameters of rock masses
    岩性容重(γ)/( kN/m3弹性模量(E)/ GPa泊松比(μ内聚力(c)/ MPa内摩擦角φ/(°)抗压强度/ MPa抗拉强度/ MPa
    上盘绢英岩26.8721.90.335.1633.817.90.86
    下盘闪长岩26.0726.30.2412.0345.227.90.50
    矿体24.740.30.1221.6342.766.02.07
    断层泥15.473.20.310.1615.00.850.12
    下载: 导出CSV 
    | 显示表格
    表  2  不同灰砂比充填体物理力学参数表
    Table  2.  Physical and mechanical parameters of filling body with different lime sand ratio
    灰砂比容重(γ)/(kN/m3弹性模量(E)/MPa泊松比(μ内聚力(c)/MPa内摩擦角(φ/(°)
    1∶417.326440.281.07731
    1∶1017.092930.300.35930
    1∶2016.951270.310.10029
    下载: 导出CSV 
    | 显示表格

    在前人研究的基础上,结合该矿的地质赋存条件和物理力学参数,建立断层与采场的数值计算模型,模型尺寸为300 m×200 m×100 m,如图3所示。计算采用摩尔库仑准则。

    图  3  断层计算三维模型
    Figure  3.  3D model of fault calculation

    边界条件为:模型Z方向上部为自由面,施加15 MPa的均布荷载,模拟上覆岩层600 m的自重荷载;模型X方向施加最大水平主应力17.72 MPa;Y方向施加最小水平主应力13.45 MPa。模型Z方向下部限制垂直方向位移,模型X、Y方向限制水平移动,工作面布置在断层下盘中。采场规格为长×宽×高=100 m×4 m×5 m,断层倾角55°。

    开挖顺序为:①一步采:按照隔一采一顺序沿X轴正方向推进,分步开挖20、40、60、80、90 m,逐步向断层推进,记录断层处应力和位移变化情况,达到计算平衡后,将模型中的弹性位移清零(实际采场开挖后围岩的弹性位移在充填之前己基本释放完成)。②充填:开挖空间充填模拟时,采用多层充填加载方式,模拟矿山实际充填情况对开挖空间重新赋值(底部1 m采用质量浓度68%、灰砂质量比1∶4进行充填;接顶0.5 m采用质量浓度68%、灰砂质量比1∶10进行充填;中间层采用质量浓度68%、灰砂质量比1∶20进行充填)。待模型重新计算平衡后,进行两步采回采。③两步采:按顺序依次进行相邻两侧两步采采场的回采,分步开挖20、40、60、80、90 m,逐步向断层推进,记录断层处应力和位移变化情况,如图4所示。

    图  4  开挖示意图
    Figure  4.  Schematic Diagram of Excavation

    一步采回采过程中,工作面推进至距断层不同距离应力演化规律如图5所示。

    图  5  一步采过程中工作面推进至距断层不同距离处应力云图:(a) 80 m;(b) 60 m;(c) 40 m;(d) 20 m;(e) 10 m
    Figure  5.  Stress nephogram of working face advancing to different distances from the fault in one-step mining process:(a) 80 m;(b) 60 m;(c) 40 m;(d) 20 m;(e) 10 m

    图5可知,一步采回采矿体时,原岩应力的平衡状态被打破,应力重新分布,造成相邻的矿体应力集中。随着一步采工作面的推进,工作面顶底板垂直应力集中程度进一步增加,在工作面端部及两侧帮1~5 m处形成规则的应力集中区,当工作面距断层距离分别为80、60、40、20、10 m时,工作面端部的垂直应力值依次分别为2.72、3.65、4.57、5.29、4.21 MPa,两侧帮的垂直应力值依次分别为1.98、2.57、2.92、3.53、4.15 MPa。

    由此可以看出,随着工作面的推进,工作面端部及两侧帮应力逐渐增加,当工作面推进至距断层10 m处,应力集中程度有减弱的趋势,工作面端部应力集中区与断层贯通,说明此处极易发生断层垮落等事故。

    两步采回采过程中,工作面推进至距断层不同距离应力演化规律如图6所示。

    图  6  两步采过程中工作面推进至距断层不同距离处应力云图:(a) 80 m;(b) 60 m;(c) 40 m;(d) 20 m;(e) 10 m
    Figure  6.  Stress nephogram of working face advancing to different distances from the fault in two-step mining process:(a) 80 m;(b) 60 m;(c) 40 m;(d) 20 m;(e) 10 m

    图6可知:一步采采空区胶结充填后,两侧的充填体与矿体形成应力的二次平衡状态,此时胶结充填体既能够限制矿体的变形,又能够与矿体共同承载。随着两步采回采矿体,造成应力的二次分布。随着两步采工作面的推进,工作面顶底板及两侧帮垂直应力较开挖前明显减小,应力集中区向外移动。当工作面距断层距离分别为80、60、40、20、10 m时,工作面端部的垂直应力值依次分别为1.65、2.37、3.25、4.43、2.91 MPa,两侧帮的垂直应力值依次分别为1.47、1.95、2.61、2.88、2.97 MPa。与一步采应力值相比可知,一步采胶结充填体可减轻采场断层处因应力集中而导致的破坏,进而为采场安全施工及矿柱回采提供便利。

    对比图6(c)和图6(d)可知:随着工作面的推进,工作面端部和两侧帮应力同样逐渐增加,当工作面推进至距断层20 m处,工作面端部应力集中区与断层应力集中区贯通,说明此时充填体可能已处于失稳状态,据此得出了该矿两步采的充填体强度需求解析解2.97 MPa。同时,由图6(d)和图6(e)可以明显看出,断层对于应力传递起到阻隔作用,在工作面上方30 m处形成应力释放区。

    通过上述研究分别得到该矿断层影响下两步采采场采充过程中充填体强度需求的解析解和数值解,对比可知解析解和数值解吻合度较好,说明基于断层影响下的充填体强度需求模型和数值模拟可以很好地计算两步采采场充填体的强度需求。经矿山调研,该矿充填质量管理规定根据经验公式求得,要求采场充填体强度大于两者,如表3所列。

    表  3  3种方式充填体强度需求
    Table  3.  Strength requirements of three filling methods
    不同计算方式解析解数值解生产实际
    单轴抗压强度/MPa2.852.974.00
    下载: 导出CSV 
    | 显示表格

    为了进一步验证优化满足矿山充填开采技术现状的采场充填体实际强度需求,下一步将与矿方沟通进行试验采场验证。

    1)通过对断层处的充填体和矿体进行力学分析,建立基于断层影响下的充填体强度需求力学模型,推导获得了极限平衡条件下充填体强度需求计算公式,得到了该矿二步采的充填体强度需求解析解为2.85 MPa。

    2)利用FLAC3D构建了两步采回采充填体与矿体的数值计算模型,可以看出:随着工作面的推进,工作面端部及两侧帮应力逐渐增加,当工作面推进至距断层10~20 m处,应力集中程度有减弱的趋势,工作面前端应力集中区与断层贯通,得到了该矿两步采的充填体强度需求解析解为2.97 MPa。

    3)通过进行充填体强度需求解析解和数值解对比,进一步验证了基于断层影响下的充填体强度需求力学模型的有效性;但矿山实际充填强度要求大于两者,需进一步验证优化。

    4)虽然本文所建立的力学模型可以反映充填体内应力演化规律,但以上研究只涉及了力学推导和数值模拟等理论验证,未开展现场充填体应力监测验证,下一步通过开展现场应力监测来进一步验证力学模型的有效性。

  • 图  1   盘区式上向水平分层充填法

    1.斜坡道;2.脉外溜井;3.分段巷道;4.分段联络巷道;5.脉内溜矿井;6.人行井;7.充填通风井;8.阶段运输平巷。

    Fig  1.   Panel upward cut and fill method

    图  2   两步采回采时断层力学分析

    Fig  2.   Mechanical analysis of fault in two step mining

    图  3   断层计算三维模型

    Fig  3.   3D model of fault calculation

    图  4   开挖示意图

    Fig  4.   Schematic Diagram of Excavation

    图  5   一步采过程中工作面推进至距断层不同距离处应力云图:(a) 80 m;(b) 60 m;(c) 40 m;(d) 20 m;(e) 10 m

    Fig  5.   Stress nephogram of working face advancing to different distances from the fault in one-step mining process:(a) 80 m;(b) 60 m;(c) 40 m;(d) 20 m;(e) 10 m

    图  6   两步采过程中工作面推进至距断层不同距离处应力云图:(a) 80 m;(b) 60 m;(c) 40 m;(d) 20 m;(e) 10 m

    Fig  6.   Stress nephogram of working face advancing to different distances from the fault in two-step mining process:(a) 80 m;(b) 60 m;(c) 40 m;(d) 20 m;(e) 10 m

    表  1   各岩体物理力学参数表

    Table  1   Table of physical and mechanical parameters of rock masses

    岩性容重(γ)/( kN/m3弹性模量(E)/ GPa泊松比(μ内聚力(c)/ MPa内摩擦角φ/(°)抗压强度/ MPa抗拉强度/ MPa
    上盘绢英岩26.8721.90.335.1633.817.90.86
    下盘闪长岩26.0726.30.2412.0345.227.90.50
    矿体24.740.30.1221.6342.766.02.07
    断层泥15.473.20.310.1615.00.850.12
    下载: 导出CSV

    表  2   不同灰砂比充填体物理力学参数表

    Table  2   Physical and mechanical parameters of filling body with different lime sand ratio

    灰砂比容重(γ)/(kN/m3弹性模量(E)/MPa泊松比(μ内聚力(c)/MPa内摩擦角(φ/(°)
    1∶417.326440.281.07731
    1∶1017.092930.300.35930
    1∶2016.951270.310.10029
    下载: 导出CSV

    表  3   3种方式充填体强度需求

    Table  3   Strength requirements of three filling methods

    不同计算方式解析解数值解生产实际
    单轴抗压强度/MPa2.852.974.00
    下载: 导出CSV
  • [1] 刘杰. 充填采矿法的应用现状及发展[J]. 当代化工研究,2020(7): 8-9.
    [2] 王锐. 探析充填采矿法的应用现状及发展方向[J]. 科技展望, 2015, 25(28): 160.
    [3] 张国柱, 杨春立. 金属矿山充填采矿方法[J]. 科技创新导报, 2015, 12(2): 251.
    [4] 吕建波. 断层影响下的工作面矿压显现规律[J]. 中国矿山工程, 2022, 51(4): 62-68.
    [5] 王同旭, 曹明辉, 江东海. 采动影响下断层活化失稳及能量释放规律研究[J]. 煤炭科学技术, 2022, 50(7): 75-83.
    [6] 蒋金泉, 武泉林, 曲华. 硬厚岩层下逆断层采动应力演化与断层活化特征[J]. 煤炭学报, 2015(2): 267-277.
    [7] 吴爱祥, 沈慧明, 姜立春, 等. 窄长型充填体的拱架效应及其对目标强度的影响[J]. 中国有色金属学报, 2016, 26(3): 648-654.
    [8] 刘光生, 杨小聪, 郭利杰. 基于采充时序的空场嗣后充填体强度要求优化研究[J]. 有色金属工程, 2019, 9(10): 85-94.
    [9] 刘光生. 充填体与围岩接触成拱作用机理及强度模型研究[D]. 北京: 北京科技大学, 2017.
    [10] 熊朝辉, 谭玉叶, 楚立申, 等. 高阶段嗣后分层胶结充填体强度及结构优化研究[J]. 金属矿山, 2022(3): 35-42.
    [11] 魏晓明, 郭利杰, 李长洪, 等. 高阶段胶结充填体强度空间变化规律研究[J]. 岩土力学, 2018(增刊2):45-52.
    [12]

    LI L. A: New concept of backfill design application of wick drains in backfilled stopes[J]. International Journal of Mining Science and Technology, 2013, 23(5): 763-770.

    [13]

    SOBHI M A, LI L, AUBERTIN M. Numerical investigation of the earth pressure coefficient along the central line of backfilled stopes[J]. Canadian Geotechnical Journal, 2017, 54(1): 138-145.

    [14]

    SOBHI M A, LI L, AUBERTIN M. Numerical investment of the lateral earth pressure coefficient along the VCL of vertical backfilled stopes[J]. Canadian Geotechnical Journal, DOI: 10.1139/cgj-2016-0165.

    [15] 赵研, 秦洪岩, 商佳新. 基于数值计算的断层群小断面巷道稳定性分析[J]. 煤炭与化工, 2022, 45(8): 26-28.
    [16] 付豪, 吴爱祥, 阮竹恩, 等. 高硫尾矿充填体强度演化规律及其机理分析[J]. 金属矿山, 2022(4):46-53.
    [17] 吉坤, 韩斌, 胡亚飞, 等. 基于强度折减与智能算法的井下充填体强度预测[J]. 中国有色金属学报, 2021(3): 796-805.
    [18] 常宝孟, 杜翠凤, 魏丁一, 等. 基于库仑摩擦原理的充填体强度力学模型[J]. 中南大学学报(自然科学版),2020(3): 777-782.
    [19] 张盛友, 孙伟, 宋登学. 上向进路全尾砂胶结充填体强度现场测试及优化研究[J]. 中国金属通报, 2018(8):271-272.
    [20] 饶运章, 徐灵彬. 某铜矿逆断层对围岩稳定性影响数值分析[J]. 有色金属科学与工程, 2012, 3(6): 55-60.
    [21]

    MITCHELL R, OLSEN R, SMITH J. Model studies on cemented tailings used in mine backfill[J]. Canadian Geotechnical Journal, 1982, 19(1): 14-28.

    [22]

    LI L.Generalized solution for minging backfill design[J]. International Journal of Geomechanics, 2014,14(3):1-11.

    [23]

    LIU G S, LI L, YANG X C, et al. Numerical analysis of stress distribution in backfilled stopes considering interfaces between the backfill and rock walls[J]. International Journal of Geomechanics, 2017, 17(1): 1-9.

  • 期刊类型引用(1)

    1. 刘杰,白腾飞,李广波,郭加仁,盛宇航. 不同胶凝材料充填体强度与水化产物含量的量化关系研究. 黄金科学技术. 2024(06): 1046-1055 . 百度学术

    其他类型引用(0)

图(6)  /  表(3)
计量
  • 文章访问数:  65
  • HTML全文浏览量:  12
  • PDF下载量:  6
  • 被引次数: 1
出版历程
  • 收稿日期:  2022-09-21
  • 修回日期:  2022-12-18
  • 网络出版日期:  2023-11-07
  • 刊出日期:  2023-10-30

目录

/

返回文章
返回