极分散深部矿体开采对地表及建构筑物的影响分析.pdf

第 卷第期 年月矿冶M I N I N GAN D ME TA L L UR G YV o l ,N o A u g u s t d o i:/j i s s n 极分散深部矿体开采对地表及建构筑物的影响分析宋洪宇马宏伟邵琳赵培显王金洋(中矿金业股份有限公司,山东 招远 )摘要:地下矿山大规模开采会对周围岩体造成扰动,会出现围岩的变形、开裂或者移动,逐渐对地表及其建构筑物造成威胁.以张家岭整合矿区深部大规模开采对地表及其建构筑物的稳定性影响为研究背景,通过辨识得到张家岭整合矿区影响范围内主要建构筑物的保护等级,进而分析影响地表变形的主要因素,结合实际工程地质条件,采用数值模拟方法,利用 D m i n e F L A C D联合建模分析方法,对张家岭整合矿区深部矿体一期及二期开采后的地表变形及地表构建筑物稳定性进行分析.结果表明:一期开采地表移动变形较小,随着二期开采的展开,地表最大竖直位移及最大水平位移逐步增大,其中地表最大沉降为 mm;分析在深部开采结束后,G 国道、尾矿库、竖井、村落、诸河流的地表倾斜、地表曲率及地表水平变形值均满足相关规范要求,因此得出张家岭矿区地下开采对该区域的地表变形影响较弱的结论.关键词:深部开采;稳定性分析;数值模拟;地表变形中图分类号:TU 文献标志码:A文章编号:()A n a l y s i so f i n f l u e n c eo fm i n i n go f e x t r e m e l yd i s p e r s e dd e e po r e b o d i e so ns u r f a c ea n ds t r u c t u r e sS ONG H o n g y u MA H o n g w e i S HAOL i n Z HAOP e i x i a n WANGJ i n y a n g(Z h o n g k u a n gG o l dI n d u s t r yC o L t d,Z h a o y u a n ,S h a n d o n g,C h i n a)A b s t r a c t:L a r g e s c a l eu n d e r g r o u n d m i n i n g w i l lc a u s ed i s t u r b a n c et ot h es u r r o u n d i n gr o c k m a s s,a n dd e f o r m a t i o n,c r a c k i n go rm o v e m e n to fs u r r o u n d i n gr o c k w i l lo c c u r,w h i c h w i l lg r a d u a l l yt h r e a t e nt h es u r f a c ea n di t ss t r u c t u r e s T a k i n gt h e i n f l u e n c eo fd e e pl a r g e s c a l em i n i n go nt h es t a b i l i t yo f t h es u r f a c ea n d i t ss t r u c t u r e si nt h eZ h a n g j i a l i n gi n t e g r a t e dm i n i n ga r e aa st h er e s e a r c hb a c k g r o u n d,t h ep r o t e c t i o nl e v e l so ft h em a i ns t r u c t u r e sw i t h i nt h ei n f l u e n c er a n g eo ft h eZ h a n g j i a l i n gi n t e g r a t e d m i n i n ga r e aa r ei d e n t i f i e d,a n dt h e nt h em a i nf a c t o r sa f f e c t i n gt h es u r f a c ed e f o r m a t i o na r ea n a l y z e d C o m b i n e dw i t ht h ea c t u a le n g i n e e r i n gg e o l o g i c a lc o n d i t i o n s,n u m e r i c a ls i m u l a t i o ni sa d o p t e d U s i n g Dm i n e F L A C Dj o i n tm o d e l i n ga n da n a l y s i sm e t h o d,t h es u r f a c ed e f o r m a t i o na n ds u r f a c es t r u c t u r es t a b i l i t ya f t e r t h em i n i n go ft h e f i r s t a n ds e c o n dd e e po r eb o d i e s i nZ h a n g j i a l i n g i n t e g r a t e dm i n i n ga r e aw e r ea n a l y z e d T h ec a l c u l a t i o nr e s u l t ss h o wt h a tt h es u r f a c e m o v e m e n td e f o r m a t i o ni ss m a l li nt h ef i r s tp h a s eo fm i n i n g W i t ht h ed e v e l o p m e n to ft h es e c o n dp h a s eo f m i n i n g,t h e m a x i m u m v e r t i c a ld i s p l a c e m e n ta n dt h e m a x i m u mh o r i z o n t a l d i s p l a c e m e n to f t h es u r f a c ea r e i n c r e a s e dg r a d u a l l y,a n dt h em a x i m u ms u r f a c es e t t l e m e n t i s mm A f t e rd e e p m i n i n g,t h es u r f a c ei n c l i n a t i o n,s u r f a c ec u r v a t u r ea n ds u r f a c eh o r i z o n t a ld e f o r m a t i o nv a l u e so fG N a t i o n a lR o a d,t a i l i n g sr e s e r v o i r,s h a f t,v i l l a g ea n dr i v e ra l lm e e t t h er e l e v a n ts t a n d a r dr e q u i r e m e n t s,s o i t i s c o n c l u d e d t h a tu n d e r g r o u n dm i n i n g i nZ h a n g j i a l i n gm i n i n ga r e ah a s aw e a k i n f l u e n c eo ns u r f a c ed e f o r m a t i o n i nt h i s r e g i o n K e yw o r d s:d e e pm i n i n g;s t a b i l i t ya n a l y s i s;n u m e r i c a l s i m u l a t i o n;s u r f a c ed e f o r m a t i o n收稿日期:第一作者:宋洪宇,高级工程师,主要从事矿山企业安全生产管理、机电技术管理等工作.E m a i l:q q c o m目前,我国发展亟待解决的问题之一就是开发深部地下的矿产资源及处理相关风险,但地下开采会伴随着岩体变形、错位、地面垮塌等问题,严重影响地表建筑物及生态环境等的稳定 .矿冶国内外对开采深部地下矿产资源时周围地表建筑物的稳定性进行了一系列深入研究.刘明宇等利用F L A C D构建原始地应力场,通过监测采场开挖顺序 对地表道 路及周围河 流 的 力 学 响应,分析开采对地表带来的影响.余正方通过理论加模拟圈定深部开采下岩石的移动规律,验证开采碎石充填塌落区域可降低表层边坡变形.陈嘉生运用离散元软件 D E C构建真实开采场景,研究其对地表建筑物的沉降影响.徐钊利用F L A C D软件,构建多种开采模型,研究河流下开采所产生的影响.张斌等针对缓倾斜中厚矿体开采,运用数值模拟监测开采下的地表稳定性,评选三种采矿方法,验证房柱嗣后充填法对地表影响最小.宋霁洪等利用F L A C D建立水文地质模型,发现开采产生的扰动受地表河流水位渗流影响较小.闫照存 利用地表实测数据,结合理论分析与数值模拟分析地下多工作面开采对地表的变形和影响.为进一步探究深部开采对地表建筑物的影响,本文以张家岭整合矿区深部大规模开采为研究背景,利用 Dm i n e F L A C D建立三维工程地质模型,通过基本力学试验及岩体质量评价对模型进行赋值,分析深部开采所产生的围岩变形规律,研究其对地表建筑物的稳定性影响.矿区基本情况张家岭整合矿区面积 k m,开采标高为 m至 m.矿区分两期开采,一期充分利用现有工程,总体采用自上而下开采,开采 m水平以上 矿 体.二 期 采 用 自 下 而 上 开 采,开 采 至 m标高位置.整合矿区的采矿方法依据矿体水平厚度来划分,矿体水平厚度在 m的采用浅孔留矿嗣后充填采矿法,矿体水平厚度大于m的采用上向水平分层充填采矿法.矿区开采区域周围分布有村落、学校、工厂、高速公路、尾矿库及竖井.矿区地层岩性较简单.矿体主要位于断裂中部或两侧节理发育带中,以黄铁绢英岩化花岗质碎裂岩为主,其饱和单轴抗压强度可达 MP a,岩石质量可达到中等.围岩以二长花岗岩为主,属未风化基岩坚硬岩组,岩性简单,为块状结构,其抗压强度为 MP a,抗拉 强 度 为 MP a,黏 聚 力 MP a,岩体质量等级属类(多数为类),工程地质条件良好.地表变形形式及稳定性影响因素当进行大规模开采活动时,上覆岩体应力会进行多次改变,扩挖到一定范围后,地表及其建筑物可能会受到扰动,引起地表变形及建筑物倾斜等灾害.引起地表变形的主要形式包括下沉、水平变形、倾斜和曲率等.地表水平、倾斜变形代表了地表移动变形的相对量,是引起建筑物破坏的重要因素.本文通过计算地表倾斜、地表曲率及地表水平变形更直观地表征地表变形.地表倾斜、曲率、水平变形计算示意图如图.图地表倾斜、曲率、水平变形计算示意图F i g S c h e m a t i cd i a g r a mf o r c a l c u l a t i o no f s u r f a c e i n c l i n a t i o n,c u r v a t u r e,a n dh o r i z o n t a l d e f o r m a t i o no fb u i l d i n g su n d e r t h e i n f l u e n c eo f s u r f a c ec u r v a t u r e)地表倾斜计算:地表下沉且沿某一方向有一定坡度值,其值为下沉量与水平间距之比:iA BWBWAA B Wl()地表曲率计算:地表下沉呈盆地后剖面的弯曲度,其值为相邻两线段倾斜差 除以两线段地表水平长度的平均值(l/m):kBiBCiA B(B CA B)i(ll)()地表水平变形计算:一线段两端的水平差与长度之比(mm/m):uEuFE F ul()矿区内岩体变形破坏和地表移动受多种因素影响和控制,这些因素共同决定着岩层移动扩展的速度、形成时间和地表变形的大小、形状及其扩展过程.结合现场工况,分析归纳张家岭整合矿区深部开采引起岩层移动和地表塌陷产生的影响因素主要有:)采矿方法与结构参数;)充填时机与回采顺宋洪宇等:极分散深部矿体开采对地表及建构筑物的影响分析序;)采空区治理;)关键地质结构;)地下水.为此,通过数值模拟进对开采影响范围内的岩层移动和地表变形进行预测计算.地表稳定性模拟研究 数值模拟建立根据矿山提供的矿体勘探线剖面图,DM i n e软件构建精确的地质模型,利用其自带接口软件导入F L A C D完成建模.三维地质模型包括第四系风化层、蚀变带、未风化基岩以及矿体,矿体赋存标高为 m,共圈定矿体 个,赋存F 蚀变带、F 蚀变带、F 蚀变带、F 蚀变带.矿体倾角为 ,厚度为 m.根据地质报告,第四系厚度为 m,风化层厚度为 m,数值模型中以保守厚度 m建立第四系风化层模型;蚀变带依据矿山地质平面图绘制建立;最终建立三维地质模型尺寸长宽高为 m m m,包含 个节点,个网格,如图所示.模型中地应力根据实测结果进行施加,采用M o h r C o u l o m b强度屈服准则,物理力学参数见表.图计算模型及矿体模型图F i g C a l c u l a t i o nm o d e l a n do r eb o d ym o d e l d r a w i n g整个计算模型通过模拟对矿体的分段开挖和充填来研究地下采矿行为对岩体移动及地表变形的影响.采用数值模拟对矿体开采过程进行简化,将矿区内 处主要矿体按分段高度进行划分,模拟实际采矿的分段下行式开采,首先将同一个水平的所有矿体进行开挖,然后将开挖的中段进行充填,计算平衡后,再进行下一个中段的矿体开挖和充填,如此向下循环.表物理力学参数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 s岩性密度/(gc m)内聚力/MP a内摩擦角/()抗拉强度/MP a弹性模量/G P a泊松比第四系风化层 蚀变带(绢英岩化花岗岩)基岩(二长花岗岩)矿体 充填体 开采对地表影响分析 矿岩位移发展规律分析为了研究矿体回采过程中上覆岩层的移动变形规律,分析随着开采活动的进行矿岩的位移变化趋势.如图所示,A A 剖面线经过点(,),倾 向 ,B B 剖 面 位 置 经 过 点(,),倾向 .通过提取典型剖面云图可大致反映围岩内部的变形移动情况,由图可知,矿体围岩变形集中在开采扰动周围的上、下盘处,一期开采上盘最大水平位移及最大竖直沉降分别为 、c m,下盘最大水平位移及最大隆起分别为 、c m.矿体下盘处由于开挖卸荷作用而出现不同程度的底鼓现象,这体现在水平位移及竖直位移均为正值,矿体上盘的水平位移及竖直位移均为负值,表明矿体上盘的变形特征为下沉.距离矿体越远,围岩变形越小,可以清楚地看出,围岩变形移动范围成拱形,并且随着与矿体距离的增加而逐渐减小.此外,一期开采阶段,靠近地表处水平变形及竖直变形均较小,这表明一期开采对地表的影响较弱.矿冶图剖面位置F i g S e c t i o np o s i t i o n图一期开采矿岩位移F i g D i s p l a c e m e n to fo r ea n dr o c kd u r i n gp h a s e Im i n i n g二期开采矿体附近围岩体的变形特征与一期开采期间围岩体的变形特征类似,同样表现出矿体下盘底鼓、矿体上盘下沉的一般规律.上盘最大水平位移及最大竖直沉降分别为 、c m,下盘最大水平位移及最大上升分别为 、c m.十分明显,由于开挖回采强度及面积增大,矿体附近围岩的变形移动强度增大.从图中还可看出,由矿体开挖导致的岩层及地表变形集中在上盘岩体及模拟区域东侧,而模拟区域西侧的岩层移动及地表变形较小,并且距离矿体越远,围岩变形越小,可以清楚地看出,围岩变形移动范围成拱形,并且随着与矿体距离的增加而逐渐减小.因此可以得出,模拟区域西侧的地表重要建(构)筑物受地下开采的影响较弱.宋洪宇等:极分散深部矿体开采对地表及建构筑物的影响分析图二期开采矿岩位移F i g D i s p l a c e m e n to fo r ea n dr o c ki np h a s e I Im i n i n g 深部开采引起的地表影响分析将一期及 二 期 开 采 结 束 后 的 地 表 沉 降 值 进行模拟与 数 据 分 析,如 表所 示,根 据 有 色金属采矿设计规范(G B ),矿区范围内地表塌陷所允许的地表倾斜i、地表曲率k和地表水平变形均满足保护等级与等级的要求.图及图为一 期及二 期 的 地 表 沉 降等值线图.表不同开采阶段的地表最大竖直及最大水平位移T a b l eT h em a x i m u mv e r t i c a l a n dh o r i z o n t a ld i s p l a c e m e n to f t h es u r f a c ea td i f f e r e n tm i n i n gs t a g e s开采阶段一期开采二期开采最大沉降/mm X方向倾斜i/(mmm)Y方向倾斜i/(mmm)X方向曲率k/(m)Y方向曲率k/(m)X方向水平变形/(mmm)Y方向水平变形/(mmm)图一期开采地表沉降等值线图(单位:mm)F i g C o n t o u rm a po f s u r f a c es u b s i d e n c ed u r i n gp h a s e m i n i n g(U n i t:mm)图二期开采地表沉降等值线图(单位:mm)F i g C o n t o u rm a po f s u r f a c es u b s i d e n c ed u r i n gp h a s e m i n i n g(U n i t:mm)矿冶 建筑物影响性分析一期开采地表沉降位移量以原家岭、张家岭主竖井为中心,向附近国道村落逐渐衰减,其中最大沉降量为mm;二期开采沉降范围同样集中在原家岭主竖井附近,开采结束后沉降量进一步增大,最大沉降量达到 mm,均满足保护等级与等级的要求.为了探究地表重要建筑物受深部开采的影响,通过提取数据对国道、竖井及地表村落的倾斜变形、曲率和水平变形值计算,如表、表所示.表一期开采地表建筑物的倾斜、曲率及水平变形T a b l eT h e i n c l i n a t i o n,c u r v a t u r e,a n dh o r i z o n t a l d e f o r m a t i o no f s u r f a c eb u i l d i n g sd u r i n gp h a s e m i n i n g建筑物X向倾斜i/(mmm)Y向倾斜i/(mmm)X向曲率k/(m)Y向曲率k/(m)X向水平变形/(mmm)Y向水平变形/(mmm)G 国道 青叶耩尾矿库 魏家沟主竖井 张家岭主竖井 魏家沟村 老翅张家 诸流河 表二期开采地表建筑物的倾斜、曲率及水平变形T a b l eT h e i n c l i n a t i o n,c u r v a t u r e,a n dh o r i z o n t a l d e f o r m a t i o no f s u r f a c eb u i l d i n g sd u r i n gp h a s e m i n i n g建筑物X向倾斜i/(mmm)Y向倾斜i/(mmm)X向曲率k/(m)Y向曲率k/(m)X向水平变形/(mmm)Y向水平变形/(mmm)G 国道 青叶耩尾矿库 魏家沟主竖井 张家岭主竖井 魏家沟村 老翅张家 诸流河 G 国道及尾矿库影响分析G 国道及青叶耩尾矿库的保护等级为级,其所允许的倾斜变形值为mm/m,水平变形为mm/m,曲率为 mm/m.通过数据分析可知,G 国道移动带内 地表的最大倾斜变形为 mm/m,最大曲率为 /m,最 大 水 平 变 形 为 mm/m;青叶耩尾矿库 移动带内地 表的最大倾 斜 变 形 为 mm/m,最大曲率为 /m,最大水平变形为 mm/m;分析结果能够满足 有色金属采矿设计规范(G B )相关要求.竖井的影响分析选取魏家沟及张家岭两个具有代表性的主竖井分析,如图、图所示,其中魏家沟主竖井围岩图魏家沟主竖井围岩倾斜、曲率及水平变形曲线图F i g C u r v eo f i n c l i n a t i o n,c u r v a t u r e,a n dh o r i z o n t a l d e f o r m a t i o no f t h es u r r o u n d i n gr o c ko f t h em a i ns h a f t i nW e i j i a g o u宋洪宇等:极分散深部矿体开采对地表及建构筑物的影响分析图张家岭主竖井围岩倾斜、曲率及水平变形曲线图F i g C u r v eo f i n c l i n a t i o n,c u r v a t u r e,a n dh o r i z o n t a l d e f o r m a t i o no f t h es u r r o u n d i n gr o c ko fZ h a n g j i a l i n gm a i ns h a f tX向倾斜最 大 值 为 mm/m,Y向 倾 斜最大值为 mm/m,整体倾斜值较小.井筒围岩X向曲率最大值约为 /m,Y向 曲 率最大值约为 /m,井筒围岩曲率变形值很小.张家岭主竖井围岩X向倾斜最大值为 mm/m,Y向倾斜最大值为 mm/m,整体倾斜值较小.井筒围岩X向曲率最大值约为 /m,Y向曲率最大值约为 /m,井筒围岩曲率变形值很小.两者分析结果均满足要求.村落的影响分析选取魏家 沟 村 和 老 翅 张 家 村 两 个 具 有 代 表性的村落分析.这类构筑物的保护等级为级,其所允许的倾斜变形值为mm/m,水平变形为mm/m,曲率为 mm/m,通过数据计算 分 析 可 知,魏 家 沟 村X向 倾 斜 最 大 值 为 mm/m,Y向倾斜最大值为 mm/m,整体倾斜值较小.X向曲率最大值约为 /m,Y向曲率最大值约为 /m.老翅张家村X向倾斜最大值为 mm/m,Y向倾斜最大值为 mm/m,整体倾斜值较小.这些区域在不同开采阶段诱发倾斜、水平变形及曲率值均较小.两者分析结果均满足规范要求.诸流河根据规范,诸流河的保护等级为级,模拟结果可知诸流河 区域地表 的沉降 数 据 为 一 期 开 采 mm、二期开采mm.通过计算该区域的倾斜变形、水平变形和曲率等数值,该区域在不同开采阶段诱发倾斜变形、水平变形及曲率值均较小.图 为一期、二期诸流河剖面的塑性区分布图,从图 中可以看出,开采扰动塑性圈未影响构造裂隙含水层岩体,因此矿体回采不会对诸流河造成影响.图 塑性区分布图F i g d i s t r i b u t i o nd i a g r a mo fp l a s t i cz o n e经模拟与计算,一期及二期开采后所有村庄建筑物的结果均能够满足 有色金属采矿设计规范(G B )相关要求.张家岭矿区地下开采对该区域的地表变形影响较弱.结论)通过辨识,张家岭矿区地表影响范围内主要的建(构)筑物及对象保护等级为G 国道、青叶耩尾矿库、魏家沟主竖井、张家岭主竖井为级;魏家沟村、老翅张家村、诸流河为级.)深部开采都会伴随着上覆岩体的应力扰动,进而发生开裂、变形、位移甚至塌落,逐步影响至地表及其建筑物.地表变形影响因素主要包括采矿方法与结构参数、采空区治理、充填时机与回采顺序、关键地质结构、地下水.)数值模拟计算表明矿体围岩变形集中在开采扰动周围的上、下盘处,距离矿体越远,围岩变形越小,可以清楚地看出,围岩变形移动范围成拱形,并且随着与矿体距离的增加而逐渐减小.同时矿冶张家岭矿区一期,地表移动变形较小;随着二期开采的展开,地表最大竖直位移及最大水平位移逐步增大,地表最大沉降为 mm.)通过计算地表各位置的倾斜变形、水平变形和曲率等数值,地表及其构筑物在不同开采阶段诱发倾斜、水 平 变 形 及 曲 率 值 均 较 小,满 足 规 范要求.参考文献李云涛,付煜,任海龙,等眼前山铁矿露天转地下开采边坡破坏过程及破坏机理研究J采矿技术,():L IY T,F U Y,R E N H L,e ta l S t u d yo nt h ef a i l u r e p r o c e s s a n d m e c h a n i s m o f s l o p e i n t h ec o n v e r s i o n o f o p e n p i tt o u n d e r g r o u n d m i n i n g o fQ i a n y a n s h a ni r o n m i n eJ M i n i n g T e c h n o l o g y,():牛淑慧迤纳厂铁铜矿采空区及地表稳定性现状评价分析 J昆 明 冶 金 高 等 专 科 学 校 学 报,():N I USHE v a l u a t i o na n da n a l y s i so ft h es t a t u sq u oo ft h e m i n e d o u t a r e a a n d s u r f a c e s t a b i l i t y o fY i n a c h a n gi r o n a n d c o p p e r m i n eJJ o u r n a lo fK u n m i n gM e t a l l u r g yC o l l e g e,():张康顺瓮福磷矿穿岩洞矿段地下开采对边坡稳定性的影响研究D贵阳:贵州大学,Z HAN GKS S t u d yo nt h e i n f l u e n c eo fu n d e r g r o u n dm i n i n go ns l o p es t a b i l i t y i nC h u a n y a n d o n go fW e n g f up h o s p h a t em i n eD G u i y a n g:G u i z h o uU n i v e r s i t y,刘明宇,王荣林,刘发平某铁矿地下开采对地表设施的 稳 定 性 影 响 分 析 J现 代 矿 业,():L I U M Y,WAN GRL,L I UFPA n a l y s i so f t h es t a b i l i t yi m p a c to fu n d e r g r o u n d m i n i n g o ns u r f a c ef a c i l i t i e so f ac e r t a i ni r o nm i n eJM o d e r n M i n i n g,():余正方大红山铁矿深部开采F 断层对岩层移动控制研究J采矿技术,():YUZFR e s e a r c ho nt h ec o n t r o l o fF f a u l t o nr o c km o v e m e n td u r i n gd e e p m i n i n go fD a h o n g s h a ni r o nm i n eJ M i n i n g T e c h n o l o g y,():陈嘉生水域动载荷条件下复杂矿体开采安全技术 D长沙:中南大学,CHE NJS As t u d yo fm i n i n gt e c h n o l o g yf o ro r ed e p o s i tb e l o w d y n a m i cl o a da n d w a t e rb o d yDC h a n g s h a:C e n t r a lS o u t hU n i v e r s i t y,徐钊某矿山地下开采对地表河流安全性影响的研究J黄金,():XU ZR e s e a r c h o n t h ei m p a c t o f u n d e r g r o u n dm i n i n go nt h es a f e t yo fs u r f a c er i v e r si nac e r t a i nm i n eJG o l d,():张斌,王孟来,李树建,等缓倾斜中厚磷矿体采矿方法设计及采动效应研究J采矿技术,():,Z HANGB,WAN G M L,L ISJ,e ta l D e s i g no fm i n i n gm e t h o d s a n ds t u d yo fm i n i n ge f f e c t s f o r g e n t l yi n c l i n e d m e d i u m t h i c kp h o s p h o r u so r eb o d i e sJM i n i n gT e c h n o l o g y,():,宋霁洪,于灯凯,孙洋,等开采扰动下三圣铁矿渗流特性 及 涌 水 风 险 分 析J黄 金,():S ONGJH,YU D K,S UN Y,e ta l A n a l y s i so fs e e p a g e c h a r a c t e r i s t i c s a n dw a t e r i n r u s h r i s ko f S a n s h e n gi r o nm i n eu n d e rm i n i n gd i s t u r b a n c eJG o l d,():闫照存黄土沟壑区多工作面开采地表形变破坏分析与评价D西安:西安科技大学,YAN Z C A n a l y s i s a n d e v a l u a t i o n o f s u r f a c ed e f o r m a t i o nd a m a g ea n df a i l u r ei n m u l t i f a c em i n i n gi n l o e s sg u l l yr e g i o nDX i a n:X i a nU 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,(编辑:周叶)