板庙子金矿深部开采留设隔离矿柱控制地压数值优化.pdf

板庙子金矿深部开采留设隔离矿柱控制地压数值优化 ∗ 赵兴东, 朱乾坤, 赵一凡 东北大学 采矿地压与控制研究中心, 辽宁 沈阳 １ １ ０ ８ １ ９ 摘 要 井下开采常引起地表产生岩体移动, 直接影响地表构 建 筑的稳定.针对板庙子 金矿深部急倾斜薄矿体开采诱发地表产生移动及其控制问题, 提出采用留设隔离矿柱与 干式充填相结合的方法, 减少或控制地表岩体移动影响范围及岩体移动速度.结合板庙 子金矿深部开采工程地质条件与生产实际情况, 分别构建无矿柱、 留１个矿柱、３个矿柱 和留两矿柱与干式充填采矿结合４种研究方案; 分别应用F L A C ３ D数值分析软件, 对４种 隔离矿柱留设方案进行数值计算, 结果表明 随着留设矿柱数量的增加, 地表沉降能够得 到有效控制, 即增加留设矿柱数量对控制地表移动能起到决定性的作用; 减少矿柱而采用 干式充填后地表沉降能得到进一步控制, 但由于充填体自身强度较低的原因仅能起到辅 助支撑作用, 通过采用隔离矿柱与干式充填的方法能更好地控制开采诱发地表产生的变 形与采场围岩的稳定. 关键词 采矿工程; 深部开采; 急倾斜矿体; 地表沉降; 数值分析 地下开采时, 矿产资源被大量采出后, 岩体原有 的平衡状态受到破坏, 上覆岩层将依次发生冒落、 断 裂、 弯曲等变形破坏, 最终波及到地表, 在采空区的 上方造成大面积的地表沉陷, 形成一个比开采面积 大得多的下沉盆地.在地表沉降影响区的外围, 会 因受拉而产生地表裂缝.矿山开采造成的地表沉降 和环境损害等一系列问题, 在我国及世界范围内越 来越引起人们的关注[ １ Ｇ ２].目前金属矿山常采用留 设矿柱与采空区充填等方法控制地表沉降, 已取得 较为理想的地表沉降控制效果, 但是目前对深部地 下金属矿床开采尤其是急倾斜矿体深部开采诱发覆 岩移动规律及地表沉降控制对策研究比较少. 关于金属矿山开采岩层移动规律方面, 夏开宗 等[ ３ Ｇ ４]以典型陡倾结构面条件下的金属矿山－程潮 铁矿西区为例, 通过对矿区的地表变形监测资料及 宏观破坏特征分析, 认为矿区的岩层移动分为采空 区顶板岩体破坏扩展至地表引起塌陷阶段和采空区 周边围岩向采空区的倾倒破坏阶段, 并得出了倾倒 滑移区的地表岩体变形规律.矿山开采地表沉降预 计通常采用概率积分法[ ５ Ｇ ６], 随着计算机技术及数值 模拟方法的引入及发展, 数值模拟方法逐渐成为地 表沉陷与岩层移动领域的主流方法之一.张国权 等[ ７]以 金 山 店 铁 矿 东 区 为 例, 利 用 离 散 元 软 件 UD E C模拟该矿区２ ＃矿体－１ ３ ０m 水平以上矿体 开挖以及－２ ７ ０m放顶工程引起的地表移动变形, 探究了两开采活动所引起的地表移动变形规律.王 彦玮等[ ８]根据望儿山金矿区浅部回采的工程地质与 采矿条件, 采用数值模拟手段对开采引起的地表沉 陷灾害进行了研究.李一帆等[ ９]利用 UD E C软件 结合某磷矿山的具体工程地质情况, 对破碎带下采 矿诱发地表沉陷的地质力学现象进行了数值模拟研 究.丁德民等[ １ ０]运用 A D I NA软件对金川镍矿不 同充填开采条件下的围岩应力场、 位移场变化特征 进行了探讨.本论文分析得出急倾斜薄矿体开采覆 岩移动破坏规律, 结合吉林省桦甸市板庙子金矿工 程地质概况, 采用F L A C ３ D软件对该矿深部开采诱 发的地表沉陷控制方法进行了研究. １ 急倾斜薄矿体开采覆岩移动诱发地表 沉降与控制方法 根据文献[１ １] 可知, 急倾斜矿体开采上覆岩层 呈现非对称破坏方式, 分为滑移变形区、 锥形破坏区 I S S N１ ６ ７ １－２ ９ ０ ０ C N４ ３－１ ３ ４ ７/T D 采矿技术 第１ ９卷 第１期 M i n i n gT e c h n o l o g y,V o l ．１ ９,N o ．１ ２ ０ １ ９年１月 J a n ． ２ ０ １ ９ ∗ 基金项目 国家自然科学基金资助项目５ １ ４ ７ ４ ０ ５ ２ ; 中央高校基本科研业务费专项资金资助项目 N １ ５ ０ １ ０ ２ ０ ０ １ ; 国家重点研发计划项 目２ ０ １ ６ Y F ０ ６ ０ ０ ８ ０ ３． 见图１ ; 由于金属矿山多为含节理裂隙岩体或完 整岩体, 开采诱发的覆岩变形破坏区域形成机制略 有不同.随着矿体采出, 采空区体积逐渐增大, 岩体 在重力作用下发生弯曲变形, 这期间顶柱及其岩体 会发生拉剪破坏而产生岩体破裂面, 然后上覆岩体 沿着破裂面或先期存在结构面产生整体滑移变形, 即为滑移变形区, 此时地表形成下沉盆地, 地表大致 可分为连续变形带和非连续变形带 见图２ . 图１ 急倾斜开采诱发岩层移动分区模型 图２ 急斜矿体开采诱发地表移动模型 在岩体上盘通常形成楔形破坏区, 楔形破坏区 尖部朝向顶板深部, 其发展演化过程为矿体上方岩 体在自重应力、 采动应力作用下发生变形, 当变形达 到一定程度后使得临空面岩体产生平行于临空面的 板状结构, 当变形持续增加, 板状结构会发生结构失 稳产生溃屈破坏乃至上盘岩体冒落, 这种破坏由上盘 围岩浅部逐渐向深部发展形成楔形破坏区 见图３ . 图３ 上盘楔形破坏区 由于岩体发生变形破坏, 会导致围岩部分区域 应力得以释放, 形成低应力区, 同时由于采动应力的 影响, 仍然存在高应力区域 见图４ . 图４ 急倾斜薄矿体开采围岩应力重新分布 针对上述覆岩移动破坏特征, 常见控制方法为 充填开采和留设隔离矿柱, 隔离矿柱设计方法有极 限跨 度 法、 经 验 公 式 法 以 及 极 限 平 衡 分 析 方 法 等[ １ ２]. ２ 板庙子金矿深部开采诱发地表沉降控制 研究 板庙子金矿位于吉林省桦甸市夹皮沟镇锦山村 板庙子屯, 位于长白山脉与张广才岭交接地带, 属丘 陵地貌, 最高峰５ ０ ３． １ ４m, 最低侵蚀基准面标高为 ３ ５ ０m, 坑口竖井标高４ ３ ０m.板庙子金矿主要开采 的矿体有４条, 即 ３ ０ ２－７号、３ ０ ２－７－１号、３ ０ ２－８ 号、 ３ ０ ３－１ ３号４条金矿体, 倾角为７ ０ ８ ９ , 矿体 为含金硫化物石英脉, 其围岩主要为花岗质片麻岩 及一系列由北西向剪切带形成的矿化破碎蚀变带. 该矿属井下开采, 采用明竖井＋两段盲竖井开拓, 明 竖井地表标高为４ ３ ５m, 井底标高为－５ ７ ５m, 井筒 开拓深度为１ ０ １ ０ m.当前矿体开 采在５ ４ ０中段 －１ ３ ０m 至９ ５ ２中段－５ ４ ２m , 共计有１ １个开 拓中段, 主要开采３ ０ ２－７号、 ３ ０ ２－８号、３ ０ ３－１ ３号 矿体; 现生产中段有５个５ ４ ０中段、５ ８ ０中段、６ ２ ０ 中段和６ ５ ３中段 ; 当前矿体开采深度为６ ７ ８m.该 矿采用上向分层干式充填采矿法, 矿房的中段高度 ４ ０m, 矿体走向长度４ ０５ ０m, 回采分层高度为３ m, 采用钢筋混凝土人工假底.深部矿房充填空区 主要分布在４号３ ３． ２号勘探线之间,５ ０ ０中段 －９ ０m标高 以上, 采空区均用废石料充填. ２．１ 板庙子金矿深部开采诱发地表沉降控制方案 针对板庙子深部开采诱发地表沉降现状, 考虑 到其现有开采方法和工程状况, 结合覆岩移动与破 ２６ 采矿技术 ２ ０ １ ９,１ ９１ 坏特征, 提出以下４种控制方案进行分析研究 方案 １为既不留设矿柱, 也不进行充填, 模拟时完全将矿 体采出; 方案２留设１隔离矿柱, 矿柱位于－２ ４ ６m 水平以上, 厚度为２ ０ m; 方案３在－１ １ ８,－２ ４ ６, －３ ８ ４m水平各留设厚度为２ ０m矿柱; 方案４在 －２ １ ０,－３ ７ ０m水平各留设厚度为２ ０m矿柱 见 图５ , 且对开采后的采场进行充填.方案１方案 ４均按从上往下开采顺序开采. 图５ 隔离矿柱留设方案４ ２．２ 数值模型建立 建模时将矿体开采模型视为平面应变模型, 假 定岩层内部为均匀连续介质.模型采用位移边界条 件, 即模型两侧和底部均为限定垂直和水平方向的 位移.在数值模拟计算过程中, 不考虑构造应力影 响, 仅考虑岩体自重引起的应力.岩体内部初始应 力状态取决于上覆岩层的重量和性质. 建模时将模型建至地表, 矿体平均厚度为２． ２ m, 倾角为７ ８ .模型尺寸为长和宽均为１ ０ ０ ０ m. 考虑到计算时间, 本次模拟采用完全弹性模型.岩 体参数来源于现场岩体质量分级和实验室实验获得 的参数 见表１ . 表１ 岩体物理力学参数 岩性 密度 / k g /m３ 弹性模量 /G P a 泊松 比 体积模 量/G P a 剪切模 量/G P a 花岗质片麻岩 ２ ８ ０ ０１ ３．８ ５０．１ ８７．２ １５．８ ７ 石英脉 矿体 ２ ８ ０ ０２ ６．６ １０．２ ２１ ５．８ ４１ ０．９ １ 采用F L A C ３ D数值模拟软件进行模拟, 并应用 T e c p l o t软件从模拟结果文件中提取位移、 应力数据 进行分析研究. ２．３ 采场围岩应力分布特征分析 从图６可以看出, 矿体采出后采出空间顶底板 垂直应力集中比较明显, 其余部位均处于或接近原 岩应力状态而未明显受开挖扰动影响, 矿柱留设数 目越多, 垂直应力受采动影响区域越小, 方案４开采 充填后垂直应力受扰动区域也明显减小.垂直应力 集中区域通常位于各采场右上角以及左下角, 矿柱 的左上角及右下角, 其它位置受开采影响垂直应力 得以释放和转移, 因此垂直应力较低.方案１中矿 体完全采出后, 在采空区左下角出现垂直应力集中 值较大, 最大可达６ ５． ０ MP a, 其它方案垂直应力集 中程度相对方案１来说降低很多, 但仍有局部小范 围垂直应力集中值接近６ ５． ０MP a. 图６ 不同方案围岩垂直应力分布云图 从图７可以看出, 矿体采出后水平应力集中区 域主要位于采空区顶底板以及矿柱附近, 且深部水 平应力普遍大于浅部围岩水平应力.方案１采空区 顶部、 方案２方案４最上部采场顶部水平应力等 值线呈现心形分布特征, 矿柱附近水平应力等值线 大致沿采场中心轴线对称分布.不同方案下水平应 力集中值无较大差别.在采场上下盘岩体中存在一 较大范围水平应力为０MP a区域, 说明在完全弹性 情况下, 矿体开挖后仅表现出对开挖空间顶底板水 平应力影响较明显, 而使上下盘岩体水平应力得以 释放或转移.方案４各应力集中区域之间也均存在 较大范围水平应力为０MP a区域, 这是由于充填体 力学参数较围岩和矿体低, 存在被压缩变形现象, 从 而使上下盘岩体水平应力能得以释放或转移, 采场 顶底板水平应力集中也比较明显. ２．４ 采场围岩位移分布特征分析 从图８可以看出, 各方案垂直位移较大处位于 采空区上下盘岩体, 垂直位移上盘岩体影响范围较 ３６ 赵兴东, 等 板庙子金矿深部开采留设隔离矿柱控制地压数值优化 大, 整个上盘岩体均出现不同程度沉降, 矿柱上方区 域垂直位移明显小于矿柱下方区域垂直位移.方案 １上盘岩体最大下沉值为１ ４ ０．０mm, 下盘岩体最大 抬升值为１ ０ ０． ０mm.方案２上盘岩体最大下沉值 为１ ０ ０． ０mm, 下盘岩体最大抬升值为８ ０．０mm.方 案３上盘岩体最大下沉值为６ ０． ０mm, 下盘岩体最 大抬升值为５ ０． ０mm.方案４上盘岩体最大下沉值 为６ ０． ０mm, 下盘岩体最大抬升值为４ ０．０mm.可 见, 方案３和方案４差别不大, 但方案４采取采场充 填开采后矿体围岩稳定性较好, 围岩变形量与影响 范围变小. 图７ 不同方案围岩水平应力分布云图 图８ 不同方案围岩垂直位移分布云图 从图９可以看, 水平位移较大处均位于上下盘 最下部采场围岩中, 其中上盘水平位移影响高度大 于下盘岩体, 采空区顶底板水平位移均比较小, 矿柱 上方区域水平位移明显小于矿柱下方区域水平的位 移.方案１水平位移最大值为１ ２ ０． ０mm, 方案２水 平位移最大值为１ ０ ０． ０mm, 方案３水平位移最大值 为７ ０． ０mm, 方案４水平位移最大值为５ ０．０mm. 从这也可以看出留设矿柱能够显著地减小上下盘岩 体的变形和移动, 与充填结合能取得更好的围岩控 制和减小地表沉陷的效果. 图９ 不同留设方案矿体围岩水平位移分布云图 ２．５ 各方案模拟结果对比分析 将４种模型模拟结果分别按水平位移、 垂直位 移、 水平应力、 垂直应力绘制曲线 见图１ ０图１ ３ . 图１ ０ ４种方案水平位移曲线图 从图１ ０可以看出, 随着矿柱留设的增加, 地表 水平位移呈现减小趋势, 尤其是上盘岩体地表位移 显著降低; 方案４采场充填后上盘水平位移减小明 显, 但下盘位移较方案３略有增大.方案１方案３ 随着矿柱数目增加地表沉陷最大值点逐渐趋于模型 垂直中心线位置, 这说明随着留设矿柱的增加, 地表 ４６ 采矿技术 ２ ０ １ ９,１ ９１ 水平位移受地下开采影响较小.方案４最大位移位 置位于模型中心线右侧, 说明方案４留设２矿柱并 充填后可使地表岩层移动范围减小, 同时可使地表 最大水平位移位置趋近于矿体上端的正上方位置. 方案１、 方案２、 方案３、 方案４的最大水平位移值分 别为２ １． ０ ７,１ ２．２ ０,５．８ ０mm和６．４ ９mm, 地表水平 移动最大值点方案１方案３分别位于其竖直中心 线左侧８． ５ １,４．２ ７,０m, 方案４则位于竖直中心线 右侧２ ３． ５ ３m处. 图１ １ ４种方案垂直位移曲线图 从图１ １可以看出, 地表下沉最大值点均处于模 型边界位置, 说明在深部开采情况下矿床开采对上 盘垂直位移影响范围较大.随着矿柱留设数目增加 地表垂直位移呈现降低趋势, 且地表沉降值降低非 常明显, 下盘地表岩体沉降影响范围逐渐减小, 同时 地表下沉曲线倾斜逐渐减小.方案４留２矿柱并充 填可使上下盘地表垂直位移减小, 下盘地表垂直方 向位移几乎无明显移动.与模型３相比, 充填相对 于增加留设矿柱数来说减沉效果不太明显, 两者结 合可有效控制地表竖直方向沉降.方案１方案４ 的垂直位移最大值分别为１ ０ ２． ５ ９,６ ４．９ ５,３ ３．６ ７, ２ ８．５ ３mm. 图１ ２ ４种方案水平应力曲线图 从图１ ２可以看出, 随着矿柱留设数目增加, 地 表岩体水平应力值逐渐降低, 且位移最大值点逐渐 趋近于模型竖直中心线, 这说明地下开采对地表水 平应力产生影响的范围逐渐减小.方案４充填后上 盘岩体地表水平应力有增加也有降低, 下盘岩体地 表水平应力既有增加也有降低, 整体来看降低位置 位于开采空间上部.方案１上盘地表岩体所受最大 拉应力为０． ９ １MP a, 位于模型中心线左侧２ ０ ７．２ ６m 处; 下盘地表岩体所受最大压应力为１． ０ ０ MP a, 位 于模型中心线右侧１ ８ ０． ３ ０m处.方案２上盘地表 岩体所受最大拉应力为０． ５ ５MP a, 位于模型中心线 左侧１ ７ ５． １ ２m处; 下盘地表岩体所受最大压应力为 ０．５ ９MP a, 位于模型中心线右侧１ ８ ３．６ ７m处.方案 ３上盘地表岩体所受最大拉应力为０．２ ６MP a, 位于 模型中心线左侧１ ５ ６． ９ １m处; 下盘地表岩体所受最 大压应力为０． ２ ７MP a, 位于模型中心线右侧１ ４ ６．７ ９ m处.方案４上盘地表岩体所受最大拉应力为０．３ ４ MP a, 位于模型中心线左侧１ ４ ６．１ ４m处; 下盘地表 岩体所受最大压应力为０． ２ ０MP a, 位于模型中心线 右侧１ ９ ３． ２ ０m处. 图１ ３ ４种模型方案应力曲线图 从图１ ３可以看出, 方案１、 方案２、 方案４地表 垂直应力比较接近, 说明在仅留一矿柱情况下和留 ２矿柱并充填对地表垂直应力影响较小.而当隔离 矿柱留设数目达到３个时, 地表所受垂直应力降低 至０． ０ ８ MP a左右, 这也反映出多矿柱开采情况下 地表垂直应力受地下开采影响较小, 充填所能起到 降低垂直应力作用有限. ３ 结 论 １针对急倾斜薄矿体开采诱发采场覆岩破坏 及移动规律进行分析, 分析得出采空区上盘岩体在 重力及采动应力作用下岩体层裂增长及失稳后破断 形成楔形破坏区, 采场上部边界以上为滑移变形区, 采空区上部存在潜在冒落区的规律.滑移变形区变 形引起地表形成下沉盆地, 地表分为连续变形区和 ５６ 赵兴东, 等 板庙子金矿深部开采留设隔离矿柱控制地压数值优化 非连续变形区. ２对于深埋急倾斜矿体, 采场岩体所受载荷 以岩体自重为主, 地表岩体移动变形以整体竖向沉 降为主, 水平方向变形移动较小, 地表沉降曲线一直 保持单沉降中心的特征, 沉降过程中沉降中心曲率 逐渐增大. ３数值模拟结果表明, 随着矿柱留设增加, 地 表沉降仍以整体竖向沉降位移为主, 水平方向变形 移动较小, 地表沉降曲线也是保持单沉降中心的尖 底形特征, 但沉降过程中尖底形曲率逐渐降低, 说明 随着矿柱的增加地表沉降得到良好的控制, 而且也 能看出增加留设矿柱数量对控制地表移动能起到决 定性的作用.在这种情况下将矿柱数减为２个并充 填采空区, 由于充填体自身强度较低的原因仅能对 地表沉降起到辅助作用, 但是能更好地控制采场上 盘围岩变形破坏.因此, 基于数值模拟结论, 将板庙 子金矿深部开采引发地表沉降控制方案定为留设两 隔离矿柱并干式充填采空区. 参考文献 [１] S e p e h r iM,A p e lDB,H a l lRA． P r e d i c t i o no fm i n i n g－ i n Ｇ d u c e ds u r f a c es u b s i d e n c ea n dg r o u n dm o v e m e n t sa taC a n a Ｇ d i a nd i a m o n dm i n eu s i n ga ne l a s t o p l a s t i c f i n i t e e l e m e n tm o d Ｇ e l[J]． I n t e r n a t i o n a lJ o u r n a lo fR o c k M e c h a n i c s M i n i n g S c i e n c e s,２ ０ １ ７１ ０ ０ ７ ３ Ｇ ８ ２． [２]付 华, 陈从新, 夏开宗, 等．地下采矿引起地表滑移变形分 析[J]．岩石力学与工程学报, ２ ０ １ ６S ２ ３ ９ ９ １ Ｇ ４ ０ ０ ０． [３]夏开宗, 陈从新, 付 华, 等．金属矿山崩落采矿法引起的岩 层移动规律分析[J]．岩土力学, ２ ０ １ ６,３ ７５ １ ４ ３ ４ Ｇ １ ４ ４ ０． [４]夏开宗, 陈从新, 付 华, 刘秀敏, 邓洋洋, 董元滨, 蒋玄苇．程 潮铁矿西区不同采矿水平下的岩体变形规律分析[J]．岩石 力学与工程学报,２ ０ １ ６,３ ５４ ７ ９ ２ Ｇ ８ ０ ５． [５]何国清, 杨 伦, 等．矿山开采沉陷学[M]．徐州 中国矿业大 学出版社,１ ９ ９ １． [６]沈 震, 徐良骥, 刘 哲, 等．基于M a t l a b的概率积分法开采 沉陷预计参数解算[J]．金属矿山, ２ ０ １ ５,４ ４９ １ ７ ０ Ｇ １ ７ ４． [７]张国权, 盛建龙, 孙明伟, 等．基于UD E C的复杂条件下放顶 引起的地表沉陷规律研究[J]．金属矿山,２ ０ １ ２,４ １６ １ ４ ５ Ｇ １ ４ ８． [８]王彦玮, 张忠辉, 邱俊刚, 等．山东望儿山金矿地表沉陷数值 模拟研究[J]．黄金科学技术, ２ ０ １ １３ ７ ２ Ｇ ７ ６． [９]李一帆, 邓 飞, 周玉斌．破碎带下磷矿开采诱发地表沉陷的 离散元数值模拟[J]．化工矿物与加工, ２ ０ ０ ６,３ ５８ ４ Ｇ ７． [１ ０] 丁德民, 马凤山, 张亚民, 等．急倾斜矿体分步充填开采对地 表沉陷的影响[J]．采矿与安全工程学报, ２ ０ １ ０,２ ７２ ２ ４ ９ Ｇ ２ ５ ４． [１ １] 王明立, 张玉卓, 张华兴．急斜煤层开采覆岩非均衡破坏机理 分析[J]．采矿与安全工程学报, ２ ０ １ ０,２ ７４ ５ ５ ８ Ｇ ５ ６ １． [１ ２] 赵兴东．谦比希矿深部开采隔离矿柱稳定性分析[J]．岩石力 学与工程学报,２ ０ １ ０,２ ９S １ ２ ６ １ ６ Ｇ ２ ６ ２ ２． 收稿日期 ２ ０ １ ８ Ｇ ０ ７ Ｇ ２ １ 作者简介 赵兴东 １ ９ ７ ５－ , 男, 辽宁辽中人, 博士生导师, 教 授, 从事采矿地压与控制研究工作,E m a i l z h a o x i n g d o n g＠ m a i l ．n e u ．e d u ．c n. 􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙􀤙 上接第５ ６页 参考文献 [１]颜宪禹, 王 迟, 周锡德．煤硫共生矿床开采特点的分析[J]． 西部探矿工程,１ ９ ９ ９２ ６ ４ Ｇ ６ ６． [２]张绍周．大红山铁矿１＃铜矿带房柱法采空区顶板矿柱稳 定性分析[D]．昆明 昆明理工大学, ２ ０ １ ４． [３]江文武, 徐国元, 马长年．基于尖点突变理论的矿房间矿柱的 稳定性分析[J]．金属矿山, ２ ０ ０ ７９ ３ ９ Ｇ ４ １＋４ ５． [４]张钦礼, 曹小刚, 王艳利, 刘洪强．基于尖点突变模型的采场 顶板Ｇ矿 柱 稳 定 性 分 析 [J]．中 国 安 全 科 学 学 报,２ ０ １ １,２ １ １ ０ ５ ２ Ｇ ５ ７． [５]陈寅聪, 乔登攀, 李克钢．基于结构可靠度矿柱稳定性分析 [J]．价值工程,２ ０ １ ７,３ ６１ １ ２ ８ Ｇ １ ３ ０． [６]郭建军, 窦源东, 杨玉泉．矿柱失稳突变学机理及对夏甸金矿 矿柱稳定性分析[J]．黄金,２ ０ ０ ８ １ １ ２ ４ Ｇ ２ ８． [７] J i a n g Y a o d o n g,W a n gH o n g w e i,Z h a oY i x i n ． T h e i n f l u e n c eo f r o a d w a yb a c k f i l lo nb u r s t i n gl i a b i l i t ya n ds t r e n g t ho fc o a l p i l l a r b yn u m e r i c a l i n v e s t i g a t i o n[C] / /T h eF i r s t I n t e r n a t i o n a l S y m p o －s i u mo nM i n eS a f e t yS c i e n c e a n dE n g i n e e r i n g 首届 矿山安全科学与工程国际学术会议论文集．北京 中国安全 生产科学研究院,２ ０ １ １． [８]周宗红, 侯克鹏, 任凤玉．分段空场崩落采矿法顶板稳定性分 析[J]．采矿与安全工程学报, ２ ０ １ ２,２ ９４ ５ ３ ８ Ｇ ５ ４ ２． [９]罗周全, 管佳林, 张旭芳, 谭浪浪．基于三维实测的复杂边界 矿柱回采崩矿方式数值优化[J]．湖南科技大学学报 自然科 学版 ,２ ０ １ ２, ２ ７１ ７ Ｇ １ １． [１ ０] 刘洪强, 张钦礼, 潘常甲, 康 虔．空场法矿柱破坏规律及稳 定性分析[J]．采矿与安全工程学报, ２ ０ １ １,２ ８１ １ ３ ８ Ｇ １ ４ ３． [１ １] 陶 明, 罗福友．宁都硫铁矿矿柱稳定性分析及回采顺序优 化[J]．采矿技术, ２ ０ １ ６,１ ６２ ２ ５ Ｇ ２ ９＋７ ０． [１ ２] 刘沐宇, 徐长佑．地下采空区矿柱稳定性分析[J]．矿冶工程, ２ ０ ０ ０１ １ ９ Ｇ ２ ２． [１ ３] 张 涛, 张 帅, 张百胜．矿柱安全留设尺寸的宽度折减法与 应用[J]．岩土力学, ２ ０ １ ４,３ ５７ ２ ０ ４ １ Ｇ ２ ０ ４ ６＋２ ０ ７ ８． 收稿日期 ２ ０ １ ８ Ｇ ０ ６ Ｇ ０ ８ 作者简介 郭军业 １ ９ ９ ２－ , 男, 甘肃武威人, 硕士研究生, 主 要从事矿业工程采矿方向的研究,E m a i l １ １ ０ ３ ８ ９ ３ ５ ４ ５＠q q ． c o m. ６６ 采矿技术 ２ ０ １ ９,１ ９１