顶板支护条件下的切顶工程优化研究及应用.pdf

第 卷第期 年月有色金属(矿山部分)N o n f e r r o u sM e t a l s(M i n i n gS e c t i o n)V o l ,N o M a y d o i:/j i s s n 顶板支护条件下的切顶工程优化研究及应用高博,李广涛,罗正良,王俊,龙赣,李聪(昆明理工大学 国土资源工程学院,昆明 ;云南黄金矿业集团股份有限公司 保山金厂河矿业有限公司,云南 保山 )摘要:以金厂河Z n V 矿体阶段空场嗣后充填采场顶部切顶条柱稳定性为研究对象,通过分析采场在不同阶段应力状态的分布规律,结合顶板破碎围岩锚固支护的现状,为提升条柱的稳定性,采用理论分析和数值模拟相结合的方法进行研究分析.研究结果表明:切顶联络道两侧交替布置m宽凿岩硐室、m宽条柱的优化方案,能有效分担应力集中作用对条柱的影响,极大提高了条柱的稳定性;优化方案改变了顶板支护区域,进一步提高了空区后顶板的稳定性.根据研究成果开展了工业试验,验证了优化方案的可行性,对安全生产具有重要意义,可在矿山进行应用.切顶工程的优化,在提高条柱稳定性的情况下,进一步提高了顶板的稳定性,可为类似矿山提供借鉴.关键词:阶段空场嗣后充填法;切顶工程优化;破碎顶板支护;数值模拟;稳定性分析中图分类号:T D 文献标志码:A文章编号:()O p t i m i z a t i o nr e s e a r c ha n da p p l i c a t i o no f r o o f c u t t i n ge n g i n e e r i n gu n d e r r o o f s u p p o r t c o n d i t i o nGAOB o,L IG u a n g t a o,L UOZ h e n g l i a n g,WANGJ u n,L ONGG a n,L IC o n g(F a c u l t yo fL a n dR e s o u r c e sE n g i n e e r i n g,K u n m i n gU n i v e r s i t yo fS c i e n c ea n dT e c h n o l o g y,K u n m i n g ,C h i n a;B a o s h a nJ i n c h a n g h eM i n i n gC o,L t d,Y u n n a nG o l dM i n i n gG r o u pC o,L t d,B a o s h a nY u n n a n ,C h i n a)A b s t r a c t:I nt h i sp a p e r,t h es t a b i l i t yo f t h er o o f c u t t i n gs t r i p s h a p e dm i n ec o l u m na t t h et o po f t h eS t a g eo p e n s t o p eb a c k f i l l i n g m e t h o d o ft h eJ i n c h a n g h e Z n V o r e b o d yi st a k e na st h er e s e a r c h o b j e c t B ya n a l y z i n gt h ed i s t r i b u t i o n l a wo ft h es t r e s ss t a t eo ft h es t o p ea td i f f e r e n ts t a g e s,c o m b i n e d w i t ht h ec u r r e n ts i t u a t i o no ft h ea n c h o r i n gs u p p o r to f t h eb r o k e ns u r r o u n d i n gr o c ko f t h er o o f,i no r d e rt oi m p r o v et h es t a b i l i t yo ft h es t r i p s h a p e dm i n ec o l u m n,t h e m e t h o do fc o m b i n i n gt h e o r e t i c a la n a l y s i sa n dn u m e r i c a ls i m u l a t i o ni su s e df o rr e s e a r c ha n da n a l y s i s T h er e s e a r c hr e s u l t ss h o wt h a t t h eo p t i m i z a t i o ns c h e m eo f a l t e r n a t e l ya r r a n g i n gm w i d ed r i l l i n gc h a m b e ra n dmw i d es t r i p s h a p e dm i n ec o l u m no nb o t hs i d e so ft h er o o f c u t t i n gc o n t a c tc h a n n e lc a ne f f e c t i v e l ys h a r et h ei n f l u e n c eo fs t r e s sc o n c e n t r a t i o no nt h es t r i p s h a p e d m i n ec o l u m na n dg r e a t l yi m p r o v et h es t a b i l i t yo ft h es t r i p s h a p e dm i n e c o l u m n T h eo p t i m i z a t i o ns c h e m e c h a n g e s t h e r o o f s u p p o r t a r e aa n d f u r t h e r i m p r o v e s t h e s t a b i l i t yo f t h er o o f a f t e r t h eg o a f A c c o r d i n gt ot h e r e s e a r c hr e s u l t s,i n d u s t r i a l e x p e r i m e n t sw e r ec a r r i e do u t t ov e r i f yt h e f e a s i b i l i t yo f t h eo p t i m i z a t i o ns c h e m e,w h i c hi so fg r e a ts i g n i f i c a n c et os a f e t yp r o d u c t i o na n dc a nb ea p p l i e di n m i n e s T h eo p t i m i z a t i o no f t h e r o o f c u t t i n gp r o j e c t f u r t h e r i m p r o v e s t h e s t a b i l i t yo f t h e r o o f u n d e r t h e c o n d i t i o no f i m p r o v i n g t h es t a b i l i t yo f t h es t r i p s h a p e dm i n ec o l u m n,w h i c hc a np r o v i d ear e f e r e n c e f o r s i m i l a rm i n e s K e yw o r d s:s t a g eo p e n s t o p eb a c k f i l l i n g m e t h o d;o p t i m i z a t i o no fr o o fc u t t i n gp r o j e c t;b r o k e nr o o fs u p p o r t;n u m e r i c a l s i m u l a t i o n;s t a b i l i t ya n a l y s i s收稿日期:基金项目:云南省基础研究专项青年项目(AU );昆明理工大学人培基金(KK Z ,KK S Y )作者简介:高博(),男,硕士研究生,采矿工程专业,主要研究方向为采矿工艺理论及技术,E m a i l:q q c o m.通信作者:李广涛(),男,高级工程师,博士,硕士生导师,采矿工程专业,主要从事采矿工艺理论研究工作,E m a i l:q q c o m.矿产资源是国家经济发展的重要基础,安全高效地进行开采尤为重要 .随着我国经济进入高质量发展时期,对环境保护的要求也更加严格.矿产资源的开发也向着安全、高效、绿色的方向发展,充填法作为一种安全、高效、环保的采矿方法,得到了极大的推广应用 .有色金属(矿山部分)第 卷大直径阶段空场嗣后充填法作为一种高效的采矿方法,为保证切顶层凿岩落矿作业的安全性,在切顶层需布置一定空间结构和数量的条柱作为对采场顶板的支撑.目前,已有众多研究人员对矿柱稳定性开展了研究.史秀志等运用突变理论推导了条柱失稳的必要条件,并结合数值模拟研究确定了大跨度采场上部硐室条柱宽度.姚高辉等基于岩体力学参数的确定,采用理论分析和数值模拟方法,对破碎围岩条件下采场内留存矿柱的稳定性进行了分析.黄敏等对传统矿柱安全系数计算过于保守,采用理论计算和数值模拟进行了分析,通过拟合得到了新的矿柱安全系数计算公式.赵国彦等通过分析房柱法开采过程中矿柱的承载机理,采用普式拱理论对矿柱尺寸进行优化,通过数值模拟验证了矿柱尺寸的合理性.叶海旺等 采用极限强度理论和数值模拟方法,分析了矿柱埋深、结构尺寸、间距等因素对矿柱稳定性的影响规律,提出了指导生产的结构参数.史秀志等 对深部大跨度采场盘区的矿柱采用能量突变分析法,导出盘区矿柱能量突跳量的稳定性判据,对盘区矿柱的稳定性进行了分析.黄慎等 基于采场顶板与矿柱力学模型的基础上,采用突变理论推导出矿柱突变失稳的充分必要条件,为优化顶板和矿柱的参数提供了借鉴.上述前人的研究,多是对于矿柱尺寸的稳定性研究,鲜有对矿柱空间结构布置下的稳定性,以及矿柱回采后空区在无矿柱支撑下顶板的稳定性进行研究.本文对切顶层结构下的应力状态分布规律进行分析,在矿山已有的矿岩体物理力学参数的成果上,基于前人对矿柱稳定性的研究成果,提出切顶层优化方案,结合采场顶板锚固支护现状,采用数值模拟方法进行计算分析,确定了切顶层的结构参数,并开展了工业应用.采矿方法及存在的问题 采矿方法金厂河Z n V 矿体赋存标高 m,距地表埋深 m.A、B采区采用阶段空场嗣后充填法两步骤连续开采,沿南北走向将采场一分为三,中间为 m胶结矿柱,南北两侧为矿房,顶板围岩破碎,采用预应力“锚杆锚索”联合支护,按照锚杆mm、锚索mm的网度布置.采场切顶层高 m,沿采场走向中线布置不连续条柱,条柱宽m、长 m不等,条柱沿走向m间隔布置,两侧凿岩硐室宽m,南北两侧采场长 m,宽均为 m.采场上部采用大直径深孔落矿,下部采用中深孔拉底,m铲运机出矿.采矿方法如图所示.图采矿方法图F i g M i n i n gm e t h o dd i a g r a m 问题分析)切顶层沿采场走向中线间断布置条柱,受切顶工程开挖扰动影响,应力再平衡过程中,条柱顶板围岩呈明显的应力集中现象并作用于条柱,导致条柱发生破坏,对条柱上部围岩稳定性也造成一定影响,如图所示.第期高博等:顶板支护条件下的切顶工程优化研究及应用)切顶工程的空间结构布置直接决定采场顶板支护区域.现有切顶工程的结构布置,造成采场顶板沿走向中线区域不能实现连续的锚固支护,加之问题中条柱位置顶板围岩出现的应力集中对顶板的稳定性造成一定影响,不利于采场空区后顶板应力平衡过程中的稳定.图采场各阶段应力状态F i g T h e s t r e s s s t a t eo f e a c hs t a g eo f t h e s t o p e 矿柱稳定性分析 矿柱荷载矿柱面积承载理论认为矿柱所承受的荷载是在开挖扰动影响下,顶板上部围岩所形成的塑性区内岩体自重荷载,单一矿柱承载的范围包括矿柱自身所占面积和相邻矿柱之间所分摊的空间顶板面积.根据Z n V 矿体阶段空场法切顶工程结构布置(图),由此得到切顶层条柱承载的平均应力为:pv(aa)(bb)ab()式中:v H,为顶板上部塑性区围岩的自重应力,MP a;塑性区围岩容重,N/m,g,g取 m/s,围岩密度,g/c m;H塑性区高度,m.a为条柱长度,m;b为条柱宽度,m;b为采场凿岩硐室的宽度,m;a为两相邻条柱之间的距离,m.图切顶层工程结构图F i g T o pc u t t i n gp r o j e c t s t r u c t u r ed i a g r a m根据普式拱理论(图),受切顶工程开挖扰动影响,采场范围内上部围岩会形成塑性区,通过解卡斯特纳方程 得到塑性区半径:RpR(pcc o t)(s i n)cc o t s i n s i n()图条柱承载机理F i g L o a d b e a r i n gm e c h a n i s mo f s t r i p s h a p e dm i n ec o l u m n有色金属(矿山部分)第 卷式中:R为开挖半径,m;P开挖处的自重应力,MP a,P H,H为顶板距地表埋深,m;c为塑性区围岩的内聚力,MP a;为塑性区围岩的内摩擦角,().切顶工程暴露后等效开挖半径为:RLh()式中:h为条柱高度,m;L为采场长度,m.为便于分析计算,将塑性区半径Rp视为顶板上部塑性区围岩的均布荷载高度,由此得到v(Rph/),从而得到条柱所承载的平均应力为:pLh(Hcc o t)(s i n)cc o t s i n s i nh(aa)(bb)ab()矿柱强度各国学者通过长期的试验和测试对矿柱强度开展了研究,取得了丰硕的成果.其中,L u n d e r 公式用于硬 岩矿柱研究,且形状效 应的精度高,与Z n V 矿体的矽卡岩相适用,公式如下:sp U Ka()()式中:U为完整矿柱单轴抗压强度,MP a;Ka为矿柱摩擦系数.Kat a nc o s CpCp()式中:Cp矿柱平均强度系数.Cp l g(b/h )b/h()由此得到矿柱安全系数F,计算公式如下:Fsppab U Ka()Lh(Hcc o t)(s i n)cc o t s i n s i nh(aa)(bb)()切顶硐室及矿柱尺寸设计矿山对破碎的顶板围岩采用预应力“锚杆锚索”联合支护,通过锚杆、锚索的协同支护作用,保证了顶板围岩的稳定性.采场空区后,顶板承受最大拉应力位置为平面几何中心处,结合前述切顶层应力状态的分析,为保证矿柱的稳定性,并使切顶工程优化后能更好地与现有的支护方式相适用,提出切顶工程沿采场中线布置切顶联道,两侧凿岩硐室和条柱垂直采场走向交替布置的优化方案,如图所示.岩体物理力学参数如表所示.凿岩硐室;条柱;切顶联道;锚固承载结构;浅层围岩破碎区;螺纹锚杆;锚索图优化后的切顶工程与锚固承载结构示意图F i g T h e s c h e m a t i cd i a g r a mo f t h eo p t i m i z e dr o o f c u t t i n gp r o j e c t a n da n c h o r a g eb e a r i n gs t r u c t u r e第期高博等:顶板支护条件下的切顶工程优化研究及应用表岩体物理力学参数T a b l eP h y s i c a l a n dm e c h a n i c a lp a r a m e t e r so f r o c km a s s名称弹性模量/G P a泊松比密度/(gc m)内摩擦角/()内聚力/MP a抗拉强度/MP a抗压强度/MP a顶板围岩 矿体 底板围岩 Z n V 矿体属微倾斜至缓倾斜矿体,矿体埋深变化幅度不大,采场顶板至地表距离按照 m考虑.矿山采用T 钻 机施工下向 大直径 深 孔,C Y T J (A)掘进台车进行顶板支护的作业,深孔间排距为m.为满足施工设备的通行,确定切顶联道宽 m,条柱长 m.凿岩硐室内采用T 钻机施工深孔时距离条柱帮壁需保持 m的施工距离,同时考虑顶板支护网度,选用m的炮孔排距,由此确定凿岩硐室宽为m.考虑条柱沿采场走向两侧炮孔排距,提出条柱分别宽 m和m两种优化方案,条柱两侧间的炮孔的排距分别为 m和m.通过理论计算,方案条柱的结构尺寸可满足稳定性,具体结果如表所示.表各方案条柱安全系数T a b l eC o e f f i c i e n t o f s a f e t yo f s t r i p s h a p e dm i n ec o l u m nf o re a c hp r o g r a m方案凿岩硐室宽度/m条柱宽度/m安全系数原方案 优化方案 优化方案 采用理论计算只能得到条柱的结构尺寸,下面通过数值模拟方法对条柱和顶板各作业时段的稳定性进行分析.数值模拟计算分析 模拟方案数值模拟按照优化方案和原方案的切顶工程结构参数布置,分别依次进行切顶工程开挖、顶板支护、矿房开挖三个模拟步骤.其中,矿房分三次进行开挖,由于三种方案中条柱沿走向的长度不同,三种方案前两次开挖步距在 m不等,使前两次开挖步距相近,便于分析比较,第三次完成矿房的剩余开挖.模型建立采用F L A C D软件内置建模,矿房模型长 m,宽 m,高 m.原方案切顶工程高 m,条柱宽m,长m,条柱沿采场走向中线间隔m布置,凿岩硐室宽m,长度为矿房长度.优化方案矿房尺寸与原方案基本一致,切顶工程按照各方案尺寸建模.考虑开挖对周边矿岩体扰动影响,建立三倍围岩,模型尺寸长 m,宽 m,高 m.锚杆、锚索支护采用C a b l e单元,按照矿山支护实际情况,锚杆 长 m,直 径 mm,锚 索 长 m,直 径 mm,分别对锚杆、锚索预紧段施加 k N和 k N预紧力,如图所示.支护材料参数见表.图矿房数值模型F i g N u m e r i c a lm o d e l o f t h em i n er o o m 模型赋值 初始地应力Z n V 矿体属浅部埋深,模型顶部边界距离地表 m,采用上部岩体自重应力,根据侧压力系数确定水平应力,计算式如下所示:v H()hv()有色金属(矿山部分)第 卷式中:为岩体容重,N/m,;v 为模型顶部距地表范围内岩体的自重应力,MP a;h为水平应力,MP a;为围岩 泊松比;H为模型顶部 距 地 表 埋深,m.表支护材料参数T a b l eS u p p o r tm a t e r i a lp a r a m e t e r s材料弹性模量/G P a直径/mm抗拉强度/(Nmm)锚固段长度/m锚杆 锚索 本构模型及边界条件数值模拟计算将矿岩体视为各向同性连续均质体,选取M o h r C o u l o m b准则.将模型四周及底部边界固定速度位移,顶部为自由边界.计算结果对比分析为全面分析数值模拟结果,更好地参照对比分析,对优化前后模拟结果按图位置截取剖面.图模拟结果方位截面图F i g A z i m u t hc r o s s s e c t i o no f s i m u l a t i o nr e s u l t s 切顶层开挖模拟结果切顶工程开挖后引起应力变化,条柱及上部顶板出现了不同程度的塑性区破坏(图).原方案受开挖扰动影响,条柱出现大范围破坏,塑性区几乎贯通,有效宽度仅有 m.其对应位置的顶板围岩受应力集中影响发生一定范围的剪切破坏.优化后,方案中条柱帮壁出现剪切破坏,条柱内部的塑性区呈现一定的贯通趋势,有效宽度为 m不等.方案中条柱帮壁虽出现剪切破坏区,但范围小,条柱有效宽度均匀控制在m,未出现塑性区贯通趋势.其次,优化方案均只在条柱端部顶板出现少量剪切破坏区,条柱顶板围岩均未出现塑性破坏区.采用F L A C D软件中“v_s h e a r_n o w”等塑性区统计指令,对切顶工程开控后的条柱塑性区体积占比(图(g)统计得到.切顶工程开控后,原方案条柱的塑性区体积占比,方案占比,方案占比.结合条柱的塑性区图和塑性区体积占比统计图来看,在满足大直径深孔孔网参数下,方案与方案相比,条柱更稳定.顶板支护下矿房开挖模拟结果顶板锚固支护后依次对矿房进行三次开挖.限于篇幅原因,同时又为保证模拟结果具有足够的代表性,这里仅展示原方案(图)与关键性方案(图)的第二次开挖结果,该步骤矿房顶板具有足够大的暴露面积,又能呈现条柱受到的应力状态.)原方案模拟结果顶板支护后,对矿房进行开挖.从A线可知,矿房空区后,顶板最小主应力为拉应力,锚固支护位置顶板拉应力值为 MP a,条柱顶板未支护位置表层围岩出现 m厚度拉应力区,应力值 MP a,上部围岩拉应力值在 MP a.矿房空区顶板最大主应力为 MP a的压应力区,靠近空区自由面条柱出现明显的压应力集中,应力值达到 MP a.由图(a)塑性区图来看,条柱顶板未支护位置在空区后,与图(a)对比,通过预应力扩散和锚索减跨作用,表层围岩由剪切破坏转变为拉破坏,塑性破坏区没有向深部围岩延伸,临近空区自由面条柱发生剪切破坏,条柱塑性区占比达(图).由A线得到,矿房边部顶板在均匀锚固支护下,顶板围岩出现均匀分布的应力区,最大主应力和最小主应力分别为 MP a的压应力和 MP a的拉应力、压应力混合区,顶板表层围岩出现 m厚度的拉剪塑性区.矿房空区位移在 mm.总体来看,除临近空区自由面的条柱发生破坏,空区采场顶板保持稳定.第期高博等:顶板支护条件下的切顶工程优化研究及应用图切顶工程开挖后的塑性区图F i g P l a s t i cz o n ed i a g r a ma f t e r e x c a v a t i o no f c u t t i n g t o pp r o j e c t有色金属(矿山部分)第 卷图原方案第二步开挖模拟结果F i g T h e s e c o n ds t e pe x c a v a t i o ns i m u l a t i o nr e s u l t so f t h eo r i g i n a l s c h e m e第期高博等:顶板支护条件下的切顶工程优化研究及应用图 方案第二步开挖模拟结果F i g T h e s e c o n ds t e pe x c a v a t i o ns i m u l a t i o nr e s u l t so fS c h e m e)优化方案模拟结果从矿房第二步开挖模拟结果来看,沿矿房走向中线B可知,在均匀锚固支护下,顶板最大主应力和最小主应力为 MP a的压应力区和 MP a的拉应力区,顶板表层围岩出现 m厚度的剪切塑性区.有色金属(矿山部分)第 卷从矿房边部B线模拟结果可得到,硐室顶板锚固区域顶板呈现 MP a拉应力,条柱对应的顶 板 未 锚 固 区 域 表 层 围 岩 拉 应 力 值 达 到 MP a,上部围岩拉应力值在 MP a.空区顶板最大压应力值为MP a,同样在临近空区自由面的条柱受到明显的压应力集中作用,应力值为 MP a.条柱顶板未支护位置受预应力扩散和锚索的减跨作用,表层围岩出现拉应力和剪应力塑性区,锚固支护位置顶板表层围岩只发生剪切破坏,上部围岩保持稳定,临近空区自由面的条柱受应力集中影响发生剪切破坏,空区顶板位移在 mm,此 时 条 柱 塑 性 区 占 比 为(图).除临近空区的条柱出现剪切破坏,空区后矿房保持稳定.图 矿房第二步开挖后的条柱塑性区占比F i g T h ep r o p o r t i o no f t h ep l a s t i cz o n eo f t h e s t r i p s h a p e dm i n ec o l u m na f t e r t h e s e c o n ds t e pe x c a v a t i o no f t h em i n e r o o m 优化前后对比分析切顶工程开挖后,优化后与优化前相比,条柱塑性区占比大幅度减少,有效控制了应力集中作用对条柱及顶板围岩的影响.条柱塑性区占比分别为、和.从矿房第二步开挖结果来看,方案与原方案相比,通过切顶层结构参数优化,改变了顶板锚固支护区域,有效控制了顶板塑性区破坏范围,采场顶板的稳定性进一步得到提高.其次,分担了矿房空区后应力集中作用对临空面条柱的影响,条柱承受的应力集中值仅为原方案的一半左右.此时,三个方案中条柱的塑性区占比分别为、和.综合来看,通过优化切顶工程结构参数,方案在满足条柱稳定的情况下,进一步提高了矿房空区顶板的稳定性,由此确定方案为最佳方案.工业应用根据方案的优化结果,对B东采区 采场(优化后)开展工业应用,并与A东采区采场(未优化)进行对比,如图 所示.随着采场回采落矿的进行,原方案切顶工程的条柱受应力集中作用的影响,条柱出现明显的“X”型剪切破坏,优化后的条柱表面没有呈现明显的结构性破坏,能更好地承受应力变化带来的影响.图 现场优化前后对比F i g C o m p a r i s o nb e f o r ea n da f t e r f i e l do p t i m i z a t i o n 结论本文以金厂河Z n V 矿体阶段空场嗣后充填法切顶工程条柱的稳定性为研究对象,根据应力状态分析和顶板锚固支护现状,通过理论分析与数值模拟相结合的方法进行了分析研究,并开展了工业应用,得出以下结论:)切顶工程凿岩硐室m宽,条柱m宽的优化方案,在满足大直径深孔孔网参数和顶板锚固支护网度的前提下,有效控制了应力集中作用对条柱的影响,极大提高了条柱的稳定性.)通过切顶工程结构参数的优化,避免了原方案中应力集中作用对沿中线位置上的条柱及顶板围岩的不利影响,实现了该区域顶板连续的锚固支护,使空区后顶板的稳定性进一步提高.)根据优化成果开展了工业试验,现场应用验证了优化结果的可行性,对进一步提升矿山的生产安全具有重要意义.第期高博等:顶板支护条件下的切顶工程优化研究及应用参考文献吴爱祥,王勇,张敏哲,等金属矿山地下开采关键技术新进展与展望J金属矿山,():WU A i x i a n g,WAN G Y o n g,Z HANG M i n z h e,e ta l N e wd e v e l o p m e n ta n dp r o s p e c to fk e yt e c h n o l o g yi nu n d e r g r o u n dm i n i n go fm e t a lm i n e sJ M e t a lM i n e,():李夕兵,周健,王少锋,等深部固体资源开采评述与探索J中国有色金属学报,():L IX i b i n g,Z HOU J i a n,WAN G S h a o f e n g,e ta l R e v i e wa n dp r a c t i c eo fd e e p m i n i n gf o rs o l i d m i n e r a lr e s p i r c e sJ T h eC h i n e s e J o u r n a l o f N o n f e r r o u s M e t a l s,():吴爱祥,王勇,王洪江膏体充填技术现状及趋势J金属矿山,():WU A i x i a n g,WAN G Y o n g,WAN G H o n g j i a n g S t a t u sa n dp r o s p e c t so ft h ep a s t eb a c k f i l lt e c h n o l o g yJ M e t a l M i n e,():于润沧我国胶结充填工艺发展的技术创新J中国矿山工程,():YUR u n c a n g D e v e l o p m e n ta n di n n o v a t i o no fc e m e n t e df i l l i n gt e c h n o l o g yi n C h i n aJ C h i n a M i n i n g E n g i n e e r i n g,():程海勇,吴爱祥,吴顺川,等金属矿山固废充填研究现状与发展趋势J工程科学学报,():C HE N G H a i y o n g,WU A i x i a n g,WU S h u n c h u a n,e t a l R e s e a r c hs t a t u sa n dd e v e l o p m e n tt r e n do fs o l i dw a s t eb a c k f i l li n m e t a l m i n e sJ C h i n e s eJ o u r n a lo fE n g i n e e r i n g,():史秀志,柯武阳,苟永刚,等大跨度采场上部硐室条形矿柱尺寸优化 研 究 J采 矿 与 安 全 工 程 学 报,():S H IX i u z h i,K E W u y a n g,G OU Y o n g g a n g,e ta l O p t i m i z a t i o no f s t r i pp i l l a rs i z ei nu p p e rc h a m b e ro fl o n g s p a ns t o p eJJ o u r n a lo f M i n i n g a n d S a f e t y E n g i n e e r i n g,():姚高辉,吴爱祥,王贻明,等破碎围岩条件下采场留存矿柱稳定性分析J北京科技大学学报,():YA O G a o h u i,WU A i x i a n g,WAN G Y i m i n g,e ta l S t a b i l i t ya n a l y s i s o f s t o p e r e t e n t i o np i l l a r s i nb r o k e n r o c kc o n d i t i o n sJJ o u r n a l o fU n i v e r s i t yo fS c i e n c e a n dT e c h n o l o g yB e i j i n g,():黄敏,KU L A T I L AK EPHS W,罗嗣海,等传统矿柱安全系数计算公式优化研究J矿冶工程,():HUAN G M i n,KU L A T I L AK EP H S W,L UOS i h a i,e ta l O p t i m i z a t i o no f t r a d i t i o n a l c a l c u l a t i o n f o r m u l a f o r s a f e t y f a c t o ro fp i l l a r sJ M i n i n ga n d M e t a l l u r g i c a lE n g i n e e r i n g,():赵国彦,周礼,李金跃,等房柱法矿柱合理尺寸设计及矿块结构参数优选J中南大学学报(自然科学版),():Z HA OG u o y a n,Z HOU L i,L IJ i n y u e,e ta l R e a s o n a b l ep i l l a rs i z ed e s i g na n dn u g g e ts t r u c t u r a lp a r a m e t e r so p t i m i z a t i o ni nr o o m a n d p i l l a rm i n i n gJ J o u r n a l o fC e n t r a l S o u t hU n i v e r s i t y(S c i e n c ea n dT e c h n o l o g y),():叶海旺,龙贵川,王其洲,等层状倾斜矿体矿柱稳定性特征及矿房参数优化J金属矿山,():Y EH a i w a n g,L ON GG u i c h u a n,WAN GQ i z h o u,e ta l S t a b i l i t yc h a r a c t e r i s t i c so f t h eo r e p i l l a r s o f l a y e r e d i n c l i n e do r e b o d ya n do p t i m i z a t i o no fm i n er o o mp a r a m e t e r sJ M e t a lM i n e,():史秀志,吝学飞,郭霆深部大 跨度采场盘区 矿柱的失稳研究J有色金属工程,():S H IX i u z h i,L I N X u e f e i,GUO T i n g S t u d yo ni n s t a b i l i t yo fp a n e l p i l l a r si nd e e ps u p e rl o n gs p a ns t o p eJ N o n f e r r o u sM e t a l sE n g i n e e r i n g,():黄慎,黄茂钧,周平矿柱的突变失稳研究J有色金属(矿山部分),():HUAN GS h e n,HUAN G M a o j u n,Z HOU P i n g C a t a s t r o p h ea n dd e s t a b i l i z a t i o no fo r ep i l l a rJ N o n f e r r o u sM e t a l s(M i n i n gS e c t i o n),():蔡美峰,何满潮,刘冬燕岩石力学与工程M北京:科学出版社,:C A IM e i f e n g,HE M a n c h a o,L I U D o n g y a n R o c k M e c h a n i c sa n dE n g i n e e r i n gM B e i j i n g:S c i e n c eP r e s s,:L UN D E RPJ,P AKA L N I SRC D e t e r m i n a t i o no f t h e s t r e n g t ho fh a r d r o c km i n ep i l l a r sJ W o r l dM i n i n gE x p r e s s,:罗基伟,张顶立,房倩,等超大跨度隧道预应力锚杆 锚索协同支护机理J中国铁道科学,():L UOJ i w e i,Z HAN G D i n g l i,F AN G Q i a n,e ta l C o m b i n e ds u p p o r tm e c h a n i s mo fp r e t e n s i o n e dr o c kb o l ta n da n c h o rc a b l ef o r s u p e r l a r g e s p a nt u n n e lJ C h i n aR a i