基于扩展的Mathews稳定图法的采场稳定性评价及参数优化_赵兴东.pdf

收稿日期2019-10-22 基金项目NSFC-山东联合基金项目 （编号 U1806208） ， 国家重点研发计划项目 （编号 2016YF0600803， 2018YFC0604401， 2018YFC0604604） 。 作者简介赵兴东 （1975） ， 男， 教授， 博士研究生导师。 基于扩展的Mathews稳定图法的采场稳定性评价 及参数优化 赵兴东牛佳安 1 （东北大学采矿地压与控制研究中心， 辽宁 沈阳 110819） 摘要为了分析大尺寸采场的稳定性并优化采场参数， 对森鑫矿业三道桥铅锌矿试验采场进行了现场地质 调查、 地质编录、 室内岩石力学试验等工作， 得到了评价矿山岩体质量的多种基础指标。对试验采场矿岩进行了岩 体质量分级， 并估算了其岩体力学参数。根据取得的结果， 利用扩展的Mathews稳定图法对试验采场在当前尺寸下 的稳定性进行了分析， 试验采场上盘的破坏概率较大， 需要进行采场结构参数的优化。基于扩展的Mathews稳定图 法对试验采场参数进行优化， 利用理论法和数值模拟对优化的结果进行验证。验证的结果与利用扩展的Mathews 稳定图法得出的结果一致， 最终确定了试验采场需保留顶柱的安全厚度为7 m。 关键词岩体质量分级Mathews稳定图采场稳定性采场参数优化数值模拟 中图分类号TD853文献标志码A文章编号1001-1250 （2020） -02-141-07 DOI10.19614/ki.jsks.202002024 Stability uation and Parameter Optimization of Stope Based on Extended Mathews Stability Graph Zhao XingdongNiu Jiaan2 （Geomechanics Research Center， Northeastern University， Shenyang 110819， China） AbstractIn order to analyze the stability of large-scale stope and optimize the parameters of the stope， field geological survey，geological catalogue and indoor rock mechanics experiment were carried out on the testing stope of Sandaoqiao Lead- Zinc Mine in Senxin Mining Company，and various basic inds for uating the rock mass quality of the mine were ob- tained. The rock mass quality classification of the rock mass in the stope is carried out，and the mechanics parameters of the rock mass were estimated. According to the results，the stability of the stope at the current size is analyzed by the extended Mathews stability graph . The probability of failure of the hanging-wall is large，so it is needed to optimize the stope structure parameters. The test stope parameters are optimized based on the extended Mathews stability graph ， and the optimization results are verified by theoretical s and numerical simulations，which are consistent with the results ob- tained by the extended Mathews stability graph . Hence， the safe thickness of the cap pillar to be retained in the stope is determined to be 7 m. KeywordsRock mass quality classification， Mathews stability graph， Stability of stope，Stope parameter optimization， Numerical simulation 在全球化经济大环境下， 矿山企业首要追求的 就是产量最大化的同时降低成本， 在合理的范围内 增大矿房尺寸能很好地适应这一发展趋势。合并采 场进行开采在增大采场尺寸的同时降低了采切工程 量和回采成本， 是金属矿山一种常见的降本增效手 段。但是， 采场尺寸的增大势必会造成采场稳定性 发生变化。 岩体是一种复杂的地质体， 通常情况下都处在 复杂的应力状态下， 承受三向应力的作用 ［1］。在采矿 过程中， 回采工作会造成明显的开采扰动， 这种开采 扰动会使采场附近一定范围内岩体的力学性质发生 变化。在一定的外部扰动作用下， 这部分岩体可能 总第 524 期 2020 年第 2 期 金属矿山 METAL MINE Series No. 524 February 2020 141 ChaoXing 发生进一步的损伤演化， 最终导致采场的失稳破坏。 因此， 合理地评价采场稳定性并优化其参数， 对于矿 山生产的安全性、 可行性、 经济性具有相当重要的意 义。 本研究在对三道桥铅锌矿大尺寸试验采场的稳 定性分析与参数优化中， 采用在国外空场法中应用 广泛的扩展的Mathews稳定图法 ［2-3］， 并借助理论法［4］ 和数值模拟对结果进行验证。基于分析与实际所得 出的结论， 为矿山提供了理论依据， 对矿山安全生 产、 提高经济效益具有重要意义并为类似工程设计 和施工提供参考。 1Mathews稳定图法及其扩展 评价采场稳定性的Mathews图表法于1980年首 次提出 ［5］， 此后， 大量的研究人员从各种的开采深度 和岩体条件中收集新的数据， 对该方法进行扩展， 并 验证了其有效性。最初的稳定图只是基于26个案例 进行研究， 经过随后几十年的扩展与改进， Mathews 稳定图的样式已经改变， 以便能适用于更多样的采 场尺寸、 岩体条件和采矿方法。 1. 1Mathews稳定图法的发展过程 自Mathews方法诞生以来， 对其进行的修改和发 展， 主要与稳定性图中各代表区域的位置和数量有 关。原始的马修斯图包含3个不同的区域， 分别为稳 定区、 不稳定区和垮落区。Potvin ［6］于1988年收集了 更多的矿山数据， 他将稳定图的区域数量减少到了2 个。该图表区域由一条过度带隔开， 分成了稳定区 和垮落区2部分。Potvin ［7］还改进了稳定图法， 使其 能够应用在受支护的采场中。Nickson ［8］和 Hadji- georgiou ［9］于1992年和1995年， 分别对Potvin改进的 稳定图法进行了完善， 他们在稳定图的数据库中添 加了更多受支护和不受支护采场的案例。Stewart和 Forsyth ［10］于1995年重新调整了Mathews稳定图， 他们 通过3条过度带将图表区域重新划分为了4部分， 分 别是可能稳定区、 可能不稳定区、 可能严重破坏区和 可能垮落区。虽然他们已经考虑到了稳定概率这个 因素， 但是没有定量的数据来表示稳定的可能性。 Trueman ［11］于2000年在稳定图的数据库中加入 了许多大尺寸采场的案例， 形成了一个包含483个案 例的综合稳定性数据库。Mawdesley ［2］于2001年基于 此数据库， 提出了扩展的 Mathews 稳定图。随后， Mawdesley ［3］于2004年验证并改进了扩展的Mathews 稳定图法。时至今日， 仍有大量的研究人员在对 Mathews稳定图法进行研究和改进。 1. 2稳定数N与水力半径R Mathews 稳定图方法的设计过程是以稳定数 N 和水力半径R这2个因子为基础进行的。其中， 稳定 数N代表岩体在给定应力条件下维持稳定的能力， 水 力半径R反映了采场的尺寸和形状。稳定数N的计 算方法如下 N Q A B C，（1） 式中，Q为假设节理水和应力折减系数均为1时计算 出的Q值， 为修正的Q系统分级法。A为岩石应力系 数， 由完整岩石单轴抗压强度与采场中线的诱导应 力的比值确定； B为节理产状调整系数， 其值由采场 面倾角与主要节理组的倾角之差来度量； C为重力调 整系数， 反映重力对采场矿岩稳定性的影响。A、 B、 C 的具体取值方法参考文献 ［2］ 。 水力半径R可以通过下式进行计算 R a l ，（2） 式中， a 为待分析采场的帮壁或者采空面的横截面 积； l为待分析采场帮壁或者采空面的周长。 1. 3扩展的Mathews稳定图法 扩展的 Mathews稳定图是基于 logistic回归分析 改进后的稳定图， Mawdesley ［3］为了消除区域划分的 主观性， 在扩展的稳定性数据库基础上， 利用Logistic 回归分析的方法确定了扩展的Mathews稳定图的稳 定破坏边界和破坏垮落边界。采用logistic回归 的优点是 可以提供基于logit模型的稳定概率预测， 利用求得的稳定概率可以在稳定图上画出稳定概率 等值线， 从而能直观反映采场的具体稳定性。 Mawdesley ［2］基于Logit模型， 采用极大似然估计 法求出Logit模型中的参数。结果如式 （3） 、 式 （4） 所 示， 该概率密度函数模型可以预测采场的稳定概 率。 z 2.960 3 - 1.442 7lnR 0.792 8lnN，（3） f z 1 1 e-z ，（4） 式中， z为预测的稳定几率值； f （z） 为预测的稳定概率 值。扩展的 Mathews 稳定图使用对数对数图 （图 1） ， 而不是传统的对数线性图。这么做的目的是 使稳定图上的边界呈线性， 图表分区更加清晰的同 时拥有显示更大采场尺寸的能力。 从图 1 可知， 由于稳定破坏边界是拟合而成 的， 所以稳定区中仍有破坏的案例。要确保采场稳 定， 采场稳定概率要尽可能高。本文根据文献 ［9］ 对 稳定破坏边界和95稳定概率等值线进行拟合， 分别求出稳定数N与水力半径R的函数表达式 lgN 1.806 5lgR - 0.651 1，（5） lgN 1.810 7lgR - 0.363 7.（6） 金属矿山2020年第2期总第524期 142 ChaoXing 赵兴东等 基于扩展的Mathews稳定图法的采场稳定性评价及参数优化 根据公式可以求出从稳定到破坏的允许水力半 径， 以及采场要达到95稳定概率所对应的水力半 径。 扩展的Mathews稳定图法在减小原方法的主观 性、 量化原方法应用中的不确定性的同时还能适用 于大尺寸采场， 是基于Mathews稳定图法改进后的一 种新方法。 2三道桥铅锌矿试验采场稳定性评价 三道桥铅锌矿位于大兴安岭山脉北段西缘， 是 集铅锌矿采、 选矿于一体的地下金属矿山。其中Ⅲ3 矿体规模最大， 为全矿区主矿体。矿体走向为286～ 345， 倾角为70～85。矿体呈脉状， 形态较规则。矿 体中部厚度较大， 深部及两侧厚度变小。 三道桥铅锌矿之前采用长50 m， 高40 m的采场 尺寸。开采至Ⅲ3 矿体 610 m 中段的 9 线穿脉附近 时， 出现地质断层使得Ⅲ3矿体发生严重错位， 导致 剩余矿体长度为80 m。该矿体长度如果布置2个采 场， 会使采切工程量和回采成本大幅增加。为达成 降低成本、 提高效益的目标， 矿山决定将原来的采场 长度由50 m提高至80 m。增大采场尺寸后， 一个采 场就能采出Ⅲ3矿体610 m中段的剩余矿石， 大大减 少了采切工程量和回采成本。但是， 随之产生的问 题是， 采场尺寸的变化会造成采场稳定性也发生变 化。对于之前尺寸的采场稳定性经验， 已不适用于 判断增大采场尺寸后的采场稳定性。因此， 需要对 大尺寸试验采场的稳定性进行合理评价并优化出合 理的采场结构参数。 2. 1试验采场 本研究的试验采场处于Ⅲ3 矿体 610 m 中段与 650 m中段之间， 如图2所示。试验采场平均跨度为 5 m， 采场高度为40 m， 采场走向长度为80 m， 采场倾 角为70。试验采场上盘围岩岩性为岩屑晶屑凝灰 岩， 下盘围岩岩性为安山岩。采矿方法选用平底结 构浅孔留矿法。 2. 2岩体质量分级 （1） RMR分类。RMR分类系统是Bieniawski ［12］于 1973年提出的确定岩体质量等级的方法。历经多次 修正， 目前广泛采用的为1989年的标准 ［13］。该方法 考虑了6个主要的因素， 即岩块单轴抗压强度， 岩石 质量指标RQD， 结构面间距， 结构面条件， 地下水条 件和结构面产状与工程走向的关系， 并以其总和值 作为岩体的RMR值。 （2） Q分类。Q分类系统由挪威岩土工程研究所 Barton等 ［14］基200个工程实例1974年提出的隧道开 挖质量分类方法。 （3） GSI分类。地质强度指标 （GSI） 分类体系由 Hoek等人 ［15］在RMR分级与Q分级的基础上提出， 目 的在于修正Hoek-Brown岩体破坏准则， 估算不同地 质条件下的岩体强度， 为工程岩体数值模拟分析提 供必要的岩体参数。 通过对三道桥铅锌矿610 m中段试验采场的工 程地质调查以及基础岩石力学试验， 得到的各项岩 体质量分级方法的结果见表1。 综合来看， 试验采场下盘岩体质量好， 上盘和矿 体的岩体质量一般。 2. 3岩体力学参数估算 针对岩石力学参数的研究来获得实际岩体力学参 数， 国内外学者在这方面进行了大量的工作， 提出了各 自的经验关系式。岩体的Hoek-Brown常数 ［16 ］、 单轴抗 压强度 ［16 ］、 单轴抗拉强度［17 ］、 弹性模量［18 ］和等效Mohr- Coulomb强度参数 ［16 ］是采场稳定性分析最重要的参 数， 在RMR、 Q和GSI岩体质量分级的基础上， 应用经 验公式估算岩体强度和变形参数。结果见表2。 2020年第2期 143 ChaoXing 2. 4采场稳定性分析 根据式 （1） 和式 （2） 分别求出试验采场的稳定数 N和水力半径R并取对数， 结果见表3。 将表3中计算得到的试验采场顶板以及上、 下盘 的稳定数的对数lgN与水力半径的对数lgR， 绘制到 扩展的稳定图中， 得出图3所示稳定图。 从图3中可以看出， 试验采场无论顶板还是上、 下盘均处于稳定破坏边界上方， 但采场上盘距离 稳定破坏边界较近， 有可能发生破坏。根据式 （3） 和式 （4） ， 即可预测出采场顶板、 上盘和下盘的稳定 几率和稳定概率。求出的具体数值见表4。 从表4可知， 试验采场的顶板、 上盘和下盘的稳 定概率约为93.3、 89.2和97.3。 从上述的稳定概率中分析可知， 采场上盘的稳 定概率小于采场下盘和顶板的稳定概率， 且不稳定 的概率相对较高， 所以按照当前的采场计划尺寸继 续回采存在安全风险， 需要优化采场结构参数。 3采场参数优化及结果验证 因 Mawdesley 概率密度函数模型得到的采场上 盘稳定概率较小， 存在破坏和垮落的风险。所以， 本 次采场参数优化的目的是调整采场结构参数， 使采 场上盘的稳定概率达到95， 从而保证矿山生产的安 全、 高效。 3. 1采场参数优化 为了进一步提高采场稳定性， 保证回采时的安 全并减少矿石损失贫化， 必须科学合理地确定采场 顶柱的安全厚度。本研究基于扩展的Mathews稳定 图法进行采场参数优化。 把采场上盘的稳定数N （见表3） 代入式 （6） 中， 可 以得到 95 稳定概率的采场上盘水力半径 R 为 12.30。将水力半径R及采上盘长度80 m代入式 （2） 中可以算出采场上盘的斜长为35.52 m。由于试验采 场倾角为 70， 则得到的优化后的采场高度为 33.38 m。即试验采场需保留厚度为7 m 的顶柱， 使整体采 场的稳定程度达到安全要求。 3. 2优化结果验证 因为这是扩展的Mathews稳定图法在三道桥铅 锌矿的首次应用， 所以对应用此法得到的结果进行 验证是十分必要的。本研究用理论法与数值模拟共 同验证优化结果， 能达到较高的准确性。 3. 2. 1理论法验证 采用理论法中应用较为广泛的极限分析法 ［4］进 行验证， 该方法具有应用方便、 过程误差小等优点。 基于极限分析法的顶柱安全厚度计算式为 ■ ■ ■ ■ ■ ■ Mp q 1 2 ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ LCnL 2 3L - CnL2 3 CL , 1 CL ＜ ■ ■ ■ ■ ■ ■ Mp q 1 2 ≤ 3 CL ll 2 3CnL2- l2 3 CCnL , ■ ■ ■ ■ ■ ■ Mp q 1 2 ＞ 3 CL ; （7） 式中，Mp为单位极限弯矩； q为顶板上的分布载荷， MPa； L 为顶板长度， 取 80 m；l为顶板跨度， 取 5 m； C、 Cn为当量系数， 不同边界条件下C和Cn的取值见 表5。 根据公式q γH可得q为4.15 MPa， 其中H为采 金属矿山2020年第2期总第524期 144 ChaoXing 场顶板距地表深度， m；γ为岩石容重， N/m3。Mp综合 反映了顶板岩性和厚度等因素， 计算式为 Mp 1 6 σ1h2，（8） 式中， h为顶板厚度， m；σ1为顶板岩体的抗拉强度， MPa； 根据表2可知σ1为0.51 MPa。 根据采场实际情况， 取C4；Cn1， 将C、Cn代入式 （7） 及式 （8） 中， 求得顶板厚度h为8.42 m。考虑到上 阶段采场保留了2 m厚的底柱， 并综合参考安全、 经 济等因素， 最终建议试验采场保留厚度为 7 m 的顶 柱。该结果与应用扩展的Mathews稳定图法所得到 的结果一致。 3. 2. 2数值模拟验证 数值模拟作为常见的验证方法有着结果直观、 适用性强等优点， 本研究利用RS2数值模拟软件验证 稳定图法的优化结果。按照实际采场尺寸、 矿体形 态和埋藏深度， 利用RS2建立数值计算模型， 模型如 图4所示。 本次模拟为了验证优化采场参数后的采场稳定 性， 因此分别对优化前的采场和优化后的采场进行 数值模拟。按现场实际矿岩条件、 采矿方法、 力学性 质等因素， 利用RS2数值模拟软件对试验采场进行开 挖计算， 参数选取见表2。 模拟结果如下 （1） 最大主应力分析。由图5可以看出， 压应力 在采场上盘与底板、 下盘与顶板交界处集中， 且优化 前的压应力较优化后更加集中。优化前， 采场围岩 受到的压应力最大值为22.48 MPa。在采场参数优化 后， 采场围岩受到的最大压应力降低至 16.55 MPa。 此时的最大压应力小于采场周围岩体的抗压强度， 采场整体呈稳定状态。模拟结果与基于扩展的 Mathews稳定图法得到的结果一致。 （2） 最小主应力分析。由图6可以看出， 拉应力 集中在采场上盘和下盘处， 且优化前的拉应力范围 较大。优化前， 采场围岩受到的拉应力最大值为 0.23 MPa。在采场参数优化后， 采场围岩受到的最大 拉应力降低至0.10 MPa。此时的最大拉应力小于采 场周围岩体的抗拉强度， 采场整体呈稳定状态。模 拟结果与基于扩展的Mathews稳定图法得到的结果 一致。 （3） 塑性区分析。由图7可以看出， 优化后的采 场塑性区体积较优化前有了明显的减小。优化前， 在采场的上盘和顶板存在较大剪切破坏的风险， 其 赵兴东等 基于扩展的Mathews稳定图法的采场稳定性评价及参数优化2020年第2期 145 ChaoXing ［1］ ［2］ ［3］ ［4］ ［5］ ［6］ ［7］ ［8］ ［9］ ［10］ ［11］ ［12］ 中采场上盘的塑性区体积最大， 塑性区深度将近 8 m， 易发生破坏。优化后， 采场的上盘和顶板发生剪 切破坏的概率大大降低， 各塑性区深度基本保持在1 m以内， 采场稳定不易破坏。模拟结果与基于扩展的 Mathews稳定图法得到的结果一致。 综上所述， 扩展的Mathews稳定图法可以用于三 道桥铅锌矿试验采场的稳定性评价及参数优化， 且 试验采场保留厚度为7 m的顶柱是安全、 合理的。 4结论 （1） 通过现场地质调查、 地质编录和室内岩石力 学试验等工作， 评价了试验采场的岩体质量。试验 采场下盘岩体质量好， 上盘和矿体的岩体质量一般。 （2） 通过扩展的Mathews稳定图法评价试验采场 的稳定性， 采场顶板和下盘较稳定， 上盘存在较大破 坏风险。利用概率密度函数模型计算出试验采场的 稳定概率， 试验采场的顶板、 上盘和下盘的稳定概率 约为93.3、 89.2和97.3。 （3） 基于扩展的Mathews稳定图法进行了采场参 数优化， 试验采场需保留厚度为7 m 的顶柱， 使整体 采场的稳定程度达到安全要求。 （4） 用理论法和数值模拟共同验证了优化结果的 准确性， 最终证明扩展的Mathews稳定图法可以用于 三道桥铅锌矿试验采场的稳定性评价及参数优化， 且 试验采场保留厚度为7 m的顶柱是安全、 合理的。 参 考 文 献 朱泽奇， 盛谦， 张勇慧， 等. 大岗山水电站地下厂房洞室群围 岩开挖损伤区研究 ［J］ . 岩石力学与工程学报， 2013， 32 （4） 734- 739. Zhu Zeqi， Sheng Qian， Zhang Yonghui，et al. Research on excava- tion damage zone of underground powerhouse of dagangshan hydro- power station［J］ . Chinese Journal of Rock Mechanics and Engi- neering， 2013， 32 （4） 734-739. Mawdesley C， Trueman R， Whiten W.Extending the Mathews stabil- ity graph for open-stope design ［J］ . Transactions of the Institution of Mining ＆Metallurgy， 2001， 10 （1） 27-39. Mawdesley C.Using logistic regression to investigate and improve an empirical design ［J］ . International Journal of Rock Me- chanics and Mining Sciences， 2004， 41 （1） 756-761. 赵兴东. 深部硬岩矿床采动地压与控制 ［M］ .北京 冶金工业出 版社， 2019 205-208. Zhao Xingdong. Deep Ming Induced Stress and Controlling at Hard Rock Deposits ［M］ . Beijing Metallurgical Industry Press， 2019 205-208. Mathews K E， Hoek E， Wyllie D C， et al.Prediction of Stable Exca- vation Spans at Depths below 1 000 Metres in Hard Rock Mines ［R］ . Ottowa Canada Centre for Mining and Energy Technology， 1981. Potvin Y. Empirical Open Stope Design in Canada ［D］ . Vancouver University of British Columbia， 1988. Potvin Y，Hudyma M R， Miller H D S.Design guidelines for open stope support ［J］ . CIM Bulletin， 1989， 82（1） 53-62. Nickson S D.Cable Support Guidelines for Underground Hard Rock Mine Operations［D］ . Vancouver University of British Columbia， 1992. Hadjigeorgiou J，Leclair J， Potvin Y. An update of the stability graph for open stope design ［C］ // CIM Rock Mechanics and Strata Control. Halifax， Nova Scotia ［s.n.］ ， 1995 14-18. Stewart S B V， Forsyth W W. The Mathews for open stope design ［J］ . CIM Bulletin， 1995， 88 （1） 45-53. Trueman R， et al. Experience in Australia with the Mathews meth- od for open stope design ［J］ . CIM Bulletin， 2000， 93 （2） 162-167. Bieniaski Z T. Engineering classification of jointed rockmasses ［J］ . 金属矿山2020年第2期总第524期 146 ChaoXing ［13］ ［14］ ［15］ ［16］ ［17］ ［18］ Transaction of the South African Institution of Civil Engineering， 1973， 15 （12） 335-344. Bieniaski Z T. Engineering Rock Mass Classification［M］ . New York Science Press， 1989 180-250. Barton N R，Lien R，Lunde J. Engineering classification of rock- masses for the design of tunnel support ［J］ . Rock Mech Rock Eng， 1974， 6 （4） 189-236. Hoek E，Marinos P，Benissi M. Applicability of the geological strength index（GSI）classification for very weak and sheared rock masses The case of the Athens Schist ation ［J］ . Bull Eng Geol Env， 1998， 57 （2） 151-160. Hoek E， Carranza-Torres C， Corkum B. Hoek-Brown failure criteri- on-2002 edition ［C］ // Hammah，Bawden WF，Curran J，et al， ed. Proceedings of the North American Rock Mechanics Society NARMS-TAC 2002. Toronto University of Toronto Press，2002 267-273. Hoek E，Wood D，Shah S. A modified Hoek-Brown criterion for jointed rock masses ［C］ // Hudson J A， ed. Proceedings of the Rock Characterization，Symposium of ISRM. London British Geotechni- cal Society， 1992 209-214. Hoek E，Diederichs M S. Empirical estimation of rock mass modu- lus ［J］ . Int J Rock Mech Min Sci， 2006， 43 203-215. 赵兴东等 基于扩展的Mathews稳定图法的采场稳定性评价及参数优化2020年第2期 147 ChaoXing