爆破震动对采空区稳定性影响的 FLAC3D 分析.pdf

返回 相似 举报
爆破震动对采空区稳定性影响的 FLAC3D 分析.pdf_第1页
第1页 / 共6页
爆破震动对采空区稳定性影响的 FLAC3D 分析.pdf_第2页
第2页 / 共6页
爆破震动对采空区稳定性影响的 FLAC3D 分析.pdf_第3页
第3页 / 共6页
亲,该文档总共6页,到这儿已超出免费预览范围,如果喜欢就下载吧!
资源描述:
第24卷 第16期 岩石力学与工程学报 Vol.24 No.16 2005 年 8 月 Chinese Journal of Rock Mechanics and Engineering Aug.,2005 收稿日期收稿日期2005–04–21;修回日期修回日期2005–06–14 基金项目基金项目国家自然科学基金重大项目50490272,50490274;国家“十五”科技攻关项目2003BA612A–10–2;中南大学博士学位论文创新工程 项目040109 作者简介作者简介 闫长斌1980–, 男, 2001 年毕业于中南工业大学建筑工程专业, 现为博士研究生, 主要从事岩土工程与地下工程方面的研究工作。 E-mail yanchangbin_2001。 爆破震动对采空区稳定性影响的爆破震动对采空区稳定性影响的 FLAC3D分析分析 闫长斌,徐国元,李夕兵 中南大学 资源与安全工程学院,湖南 长沙 410083 摘要摘要针对厂坝铅锌矿受到乱采乱挖破坏的情况,运用 FLAC3D对爆破震动作用下采空区的稳定性进行了数值分 析。根据现场实测数据,选取了 2 条不同的实测速度时程曲线作为爆破震动输入,分别进行了静力和动力计算。 分析结果表明1 采空区开挖后,最大位移量为 4.75 mm;2 施加爆破动载荷后,最大位移量增至 5.47 mm, 增幅与爆破震动速度成正比; 3 围岩位移分布的变化受震动波形的影响。 爆破震动作用下围岩应力场发生二次重 分布,塑性区分布面积明显扩大,采空区右上方露天边坡脚亦受到影响。分析结果与实际观测结果相吻合,并给 出了避免爆破震动诱发采空区失稳塌陷的几条建议。 关键词关键词爆破工程;采矿工程;爆破震动;采空区;稳定性;FLAC3D;数值模拟 中图分类号中图分类号TD 235;TD 85 文献标识码文献标识码A 文章编号文章编号1000–6915200516–2894–06 STABILITY ANALYSIS OF MINED-OUT AREAS INFLUENCED BY BLASTING VIBRATION WITH FLAC3D YAN Chang-bin,XU Guo-yuan,LI Xi-bing School of Resources and Safety Engineering,Central South University,Changsha 410083,China AbstractFor the reason that Changba lead-zinc mine is severely destroyed by disordered mining,the numerical analysis of mined-out areas under blasting vibration is carried out with FLAC3D . Two different measured velocity time-dependent curves are selected as blasting vibration according to the in-situ test data;and static calculation and dynamic calculation are carried out,respectively. The analytical results show that the maximum displacement of surrounding rock is 4.75 mm after mined-out areas are excavated,while the maximum displacement of surrounding rock increases to 5.47 mm after the blasting dynamic loading is applied,and amplitude of increasing displacement is proportional to the velocity of blasting vibration. The shape of vibration wave influences the variation of displacement distribution for surrounding rock. When the blasting vibration dynamic loading is applied, the second redistribution of stress field of surrounding rock occurs; and the distribution area of plastic zone is expanded obviously. The foot of the right open slope above the mined-out area is also influenced as well. The analytical results are accordant with the observed results. Several suggestions to avoid instability and subside of mined-out areas induced by blasting vibration are brought forward. Key wordsblasting engineering;mining engineering;blasting vibration;mined-out areas;stability;FLAC3D; numerical simulation 1 引引 言言 爆破震动是爆破公害之一,不但对地面建构 筑物的安全产生影响,而且对采空区等地下洞室的 围岩产生损伤和扰动,甚至引起失稳破坏[1 ~7]。对 爆破震动危害的研究,主要是进行现场测试,根据 实测数据,结合爆破安全规程,确定最大震速和安 第 24 卷 第 16 期 闫长斌等. 爆破震动对采空区稳定性影响的 FLAC3D分析 2895 全距离。该方法能直观地反映爆破震动危害,但无 法对爆破震动引起的围岩动应力场、位移和塑性区 的变化以及动力失稳机制进行研究。而数值模拟方 法可以系统地研究爆破作用下围岩的动应力场分布 规律及动静应力场的叠加作用机理,因此成为近年 来研究爆破震动作用的有效方法[5 ~8]。鉴于 FLAC3D 在岩土力学分析中的优势和求解动力问题的特点, 使得其可很好地解决非线性动力分析问题。FLAC3D 已成功应用于岩土开挖、边坡稳定分析及地震动力 响应分析等许多领域[9 ~12]。 厂坝铅锌矿区位于甘肃省成县的黄渚镇,由厂 坝、小厂坝和李家沟组成,是我国特大型铅锌矿床 之一。20 世纪 80 年代后期以来,当地乱采乱挖留 下的数百万平方米未处理群采空区,给矿山安全生 产构成了严重威胁。群采空区诱发的地质灾害时有 发生,致使矿山曾连续发生重大事故[13 ~15]。特别是 近年来,爆破地震波作用下大范围采空区发生动力 失稳的趋势日益明显。为保障露天转地下过渡层的 安全生产,进行了爆破地震效应现场测试研究,获 得了较准确的基础数据。为进一步研究复杂采空区 在爆破震动作用下产生动力失稳的过程和机理,本 文在现场测试的基础上,同时考虑开挖过程和爆破 震动对采空区稳定性的影响,利用 FLAC3D进行了 三维数值模拟研究。 2 FLAC3D动力计算步骤动力计算步骤 FLAC3D是二维的有限差分程序 FLAC2D的扩 展,能够进行岩石、土质和其他材料在达到屈服极 限后经历塑性变形的三维空间行为分析,为岩土工程 领域求解三维问题提供了一种理想的分析工具[9]。 运用 FLAC3D进行动力计算,必须首先进行静 力分析。在完成静力分析的基础上,才能施加动力 载荷进行动力分析。FLAC3D动力计算大致可以分为 以下几个步骤 1 确定计算区域,并进行网格划分; 2 选择动力计算模式,定义本构模型和材料 的物理力学参数; 3 定义计算所需的边界条件和初始条件; 4 进行计算,获得初始平衡状态,即开挖前 的原岩应力状态; 5 进行工程采空区开挖计算分析, 得到开挖 后的静力计算结果; 6 检查静力计算结果,认为满意后,设置动 力计算边界条件和所需的阻尼; 7 施加动力载荷,进行动力计算分析,得到 动力计算结果。一般的计算流程如图 1 所示。 图 1 FLAC3D 动力分析流程图 Fig.1 Flow chart of FLAC3D dynamic analysis 3 计算分析计算分析 3.1 计算模型计算模型 3.1.1 模型范围 根据前期采空区探测成果,2矿体存在范围较 大、互相贯通的采空区。因此,本次计算区域选择 2矿体及周边围岩。计算区域尺寸为长 120 m,宽 180 m,高 125 m。采空区位于该区域的中部,即 1 202~1 250 分段。贯通的采空区尺寸为高 24 m, 横断面为下宽上窄,底宽 36 m,顶宽 11 m,沿矿体 延伸 32 m。上边界取至地表,下边界距离采空区底 板 70 m, 左右边界踞离采空区分别为 90 m 和 60 m, 沿矿体方向Y 方向取 88 m。根据圣维南原理,模 型大小满足计算精度要求。 3.1.2 模型介绍 采空区开挖后,模型的 FLAC3D剖分网格如图 2 所示。该模型共有 70 200 个单元和 75 733 个节点; 2896 岩石力学与工程学报 2005 年 图 2 开挖后的 FLAC3D模型 Fig.2 FLAC3D model after excavation 按岩性参数的不同,可将模型分为 3 个区域,即采 空区、矿体和围岩。模型中的采空区的尺寸、形态 和位置以前期探测结果为准。根据现场调查,模型 范围内岩性比较简单,可分为围岩和矿石 2 大类, 矿体上下盘围岩均为黑云母片岩,矿石为铅锌矿石。 为便于计算,模型简化为准三维情况,即不考虑断 面形状和尺寸在长度方向Y 方向上的变化。 3.1.3 初始应力场和边界条件 初始地应力场由于采空区埋深较浅,且矿区 构造应力不明显,因此未进行现场地应力测试,初 始地应力场仅按自重应力场考虑。根据弹性力学原 理可知,竖向应力和水平应力分别为 ⎭ ⎬ ⎫ vh v σσ γσ k H 1 式中H为埋置深度;k为侧压力系数,且k 1/−,为泊松比。 边界条件采用位移边界条件,即模型的左右 X方向边界、前后Y方向边界和底边界均施加位 移约束条件,上边界为自由边界。 3.2 岩体物理力学参数岩体物理力学参数 本次计算模型中岩性较简单,仅分为矿石和围 岩 2 种。计算所需的岩体物理力学参数是根据室内 试验确定的岩石物理力学参数[16]折减得来的,具体 参数结果见表 1。 表表 1 厂坝铅锌矿岩体物理力学参数厂坝铅锌矿岩体物理力学参数 Table 1 Physico-mechanical parameters of the rocks in Changba lead-zinc mine 岩石种类 抗压 强度 /MPa 抗拉 强度 /MPa 变形 模量 /GPa 粘聚力 /MPa 内摩 擦角 / 泊 松 比 密度 /kgm -3 铅锌矿石 7.716 1 0.205 2 13.220.629 2 71.47 0.262 720 黑云母片岩 8.569 3 0.245 9 12.980.725 8 70.77 0.222 720 在 FLAC3D计算中,岩体变形参数采用的是体 积模量K和剪切模量G。 因此, 必须将弹性模量或 变形模量E和泊松比转化成体积模量K和剪 切模量G,转化公式如下 21 3− E K 2 1 2 E G 3 3.3 计算方案计算方案 为研究采空区开挖形成的围岩应力重分布和产 生的位移变化,找出爆破震动作用下采空区围岩应 力发生二次重分布的规律和位移变化,确定爆破震 动对采空区稳定的危害程度,运用 FLAC3D对开挖 和爆破震动作用下采空区的稳定性进行了较为系统 的模拟分析。分析结果可为预防采空区失稳塌陷和 减震、加固、治理提供依据。根据震动载荷形式, 共有以下 3 种计算方案1 自重应力场中,采空 区开挖过程计算;2 采空区在中深孔爆破载荷作 用下的动力计算;3 采空区在斜坡道掘进爆破载 荷作用下的动力计算。本次计算采用 Mohr-Coulomb 准则,即弹塑性模型,并选取大应变计算模式。 4 爆破震动输入爆破震动输入 4.1 爆破震动速度时程爆破震动速度时程 为获得第一手资料,对厂坝铅锌矿厂坝矿区进 行了爆破震动现场测试研究。测试分两个部分,一 是中深孔爆破强制崩矿地震效应测试,二是斜坡道 掘进爆破地震效应测试。有关研究和实践证明,地 震波的震动强度可用质点振动位移、速度、加速度 3 个物理量来表达。而震动速度既能反映震动能量 大小,又与爆破产生的动应力成正比,且测试方便 可靠,因此观测物理量以质点振动速度为主。 测试仪器为成都中科动态仪器有限公司研制的 IDTS3850 二通道和三通道爆破振动记录仪,传感器 采用 SCD–Z 型垂直和 SCD–P 型水平速度传感器。 本次计算所选用的速度时程曲线共两条,一条 来自中深孔爆破强制崩矿地震效应测试,作为动载 荷一,如图 3 所示;另一条来自斜坡道掘进爆破斜 坡道掘进爆破地震效应测试,作为动载荷二,如图 4 所示。爆破震动速度时程曲线的主要参数为1 动 载荷一的主频为 5.64 Hz,震动时间为 1.00 s,最大 幅值为 0.181 2 m/s; 2 动载荷二的主频为 12.28 Hz, 第 24 卷 第 16 期 闫长斌等. 爆破震动对采空区稳定性影响的 FLAC3D分析 2897 图 3 中深孔爆破震动速度时程曲线 Fig.3 Velocity time-dependent curves of medium-depth bore blasting concussion 图 4 斜坡道掘进爆破震动速度时程曲线 Fig.4 Velocity time-dependent curves of decline driving blasting vibration 振动时间为 5.15 s,最大幅值为 0.460 9 m/s。 4.2 爆破震动载荷施加爆破震动载荷施加 在 FLAC3D动力计算中,动载荷输入可以采用 加速度时程、速度时程、位移时程和应力时程 4 种 方式。若采用粘滞边界条件,则必须输入速度时程 进行分析[2]。由于测试数据为质点振动速度,因此 动载荷输入选用爆破震动速度时程。由于测试数据 均来自采场或巷道的底板,因此将爆破震动载荷施 加于采空区的底板,沿Z轴负方向传播。 5 模拟结果分析模拟结果分析 整个计算过程分为静力计算和动力计算。动力 分析以静力计算为前提,即在完成静力计算的基础 上,再施加动力时程进行计算。静力分析时,整个 采空区分 8 个步骤开挖形成,即沿Y方向进行开挖, 每一步开挖 4 m。每一步开挖后都计算至平衡状态, 直到整个空区开挖完成为止。将两条速度时程通过 表的形式施加于采空区底板,分别进行动力计算, 得到动力分析结果。 5.1 应力分析应力分析 图 5,6 分别为采空区开挖后和施加爆破动载后 最大主应力、最小主应力分布图。可以看出,采空 区开挖后应力发生重分布,边界和拐角处出现较大 应力集中,采空区右上方围岩和底板出现较大的拉 应力,且延伸至露天坑的右边坡脚附近;围岩中最 大拉应力为 1.096 MPa,大于岩体的抗拉强度 0.245 9 MPa,右上方围岩将发生拉裂破坏;施加动载后, 采空区底板处的拉应力集中已经不明显,但是右上 方围岩中仍然存在较大拉应力集中,且最大拉应力 仍大于岩体的抗拉强度, 故应采取有效的支护措施, 预防采空区右上方围岩失稳塌陷。动载荷一和动载 荷二对最大主应力、最小主应力分布图的影响基本 类似,所以没有在图中具体进行区分。 a 开挖后 b 动载作用后 图 5 最大主应力分布单位MPa Fig.5 Distribution of maximum principal stressesunitMPa a 开挖后 b 动载作用后 图 6 最小主应力分布单位MPa Fig.6 Distribution of minimum principal stressesunitMPa 5.2 位移分析位移分析 图 7 为采空区开挖和施加动载后围岩位移分布 图。可见,最大位移出现在采空区上边界,该边界 形状极不规则,且已基本形成“危岩”形势。采空 区开挖后围岩最大位移量大约为 4.75 mm,最大位 移发生在右上方边界围岩处,同时底板和包括右上 方露天边坡在内的区域也出现了较大的位移量;采 空区受动载荷作用后,最大位移仍然发生在右上方 边界围岩处,最大位移量较静力计算开挖计算结 果,有所增大。采空区受动载荷一作用后,围岩最 大位移量为 5.40 mm,而受动载荷二作用后,围岩 最大位移量达到 5.47 mm,增幅与最大震动速度成 2898 岩石力学与工程学报 2005 年 a 开挖后 b 动载一作用后 c 动载二作用后 图 7 位移分布单位m Fig.7 Distribution of displacementsunitm 正比。动载荷一作用后,位移量较大的区域范围主 要集中在采空区右上方,包括整个露天边坡;动载 荷二作用后,位移量较大的区域范围主要集中在采 空区右上方和底板,但底板处的范围较静力计算小。 这说明,施加动载荷后,不仅增大了围岩的最大位 移量,而且改变了位移的分布形式和范围。 5.3 塑性区分析塑性区分析 图 8 为采空区开挖后和施加爆破动载后塑性区 分布图。从塑性区分布情况来看,开挖后塑性区分 布较广,以剪切塑性区为主,集中在采空区右上方 和底板附近,且露天坑的右边坡脚处也出现塑性区。 很明显,采空区受爆破动载荷作用后,塑性区分布 范围比静力计算结果大,且动载荷二作用后比动载 荷一作用后,塑性区范围大。说明塑性区分布面积 也与动载荷幅值成正比。 a 开挖后 b 动载一作用后 c 动载二作用后 图 8 塑性区分布 Fig.8 Distribution of plastic zones 6 结论与建议结论与建议 1 采空区开挖后,引起围岩应力重分布和应 力集中;施加爆破动载荷后,围岩应力场发生二次 重分布;围岩最大位移量增大,且增幅与最大震动 速度成正比;爆破震动作用不仅增大了围岩最大位 移量,而且改变了位移的分布形式和范围。 2 已连通、范围大、且形状不规则的采空区 开挖后塑性区范围较大,且延伸至右上方露天坑边 坡脚;施加爆破动载荷后,塑性区会进一步扩大, 扩大幅度与动载荷幅值成正比。 3 复杂采空区在爆破震动荷载作用下发生失 稳塌陷的危险性增加。 4 避免爆破动载荷引发采空区顶板塌陷的几 条建议 ① 尽量减少最大段药量,并采用微差爆破等减 震措施。 ② 对出现的塑性区, 应采取喷射混凝土、 锚杆、 锚索等有效支护措施。 ③ 应对露天坑右边坡进行监测,必要时采取加 固措施。 参考文献参考文献References [1] 言志信,王永和,江 平,等. 爆破地震测试及建筑结构安全标 准研究[J]. 岩石力学与工程学报, 2003, 2211 1 907–1 911.Yan Zhixin,Wang Yonghe,Jiang Ping,et al. Study on measurement of blast-induced seism and building safety criteria[J]. Chinese Journal of Rock Mechanics and Engineering,2003,22111 907–1 911.in Chinese [2] Singh P K. Blast vibration damage to underground coal mines from adjacent open-pit blasting[J]. International Journal of Rock Mechanics and Mining Sciences,2002,39959–973. [3] Yang R L,Rocque P, Katsabanis P, et al. Measurement and analysis 第 24 卷 第 16 期 闫长斌等. 爆破震动对采空区稳定性影响的 FLAC3D分析 2899 of near-field blast vibration and damage[J]. Geotechnical and Environmental Engineering,1994,123169–182. [4] 言志信, 吴德伦, 王 漪, 等. 地震效应及安全研究[J]. 岩土力学, 2002, 232 201–203.Yan Zhixin, Wu Delun, Wang Yi, et al. The study of blast vibration effect and safety[J]. Rock and Soil Mechanics,2002,232202–204.in Chinese [5] 谭忠盛,杨小林,王梦恕. 复线隧道施工爆破对既有隧道的影响 分析[J]. 岩石力学与工程学报,2003,222281–285.Tan Zhongsheng, Yang Xiaolin, Wang Mengshu. Effect of blast in double line tunnel on existing tunnel[J]. Chinese Journal of Rock Mechanics and Engineering,2003,222281–285.in Chinese [6] 刘亚群, 李海波, 赵 坚. 确定岩质边坡安全阈值的新方法[J]. 爆 炸与冲击, 2004, 245 448–452.Liu Yaqun, Li Haibo, Zhao Jian. A new to determine safety threshold velocity of slope[J]. Explosion and Shock Waves,2004,245448–452.in Chinese [7] 李宏男, 王炳乾, 林 皋. 爆破地震效应苦干问题的探讨[J]. 爆炸 与冲击,1996,16166–72.Li Hongnan,Wang Bingqian,Lin Gao. Studies on some problems on explosion earthquake responses[J]. Explosion and Shock Waves,1996,16161–67.in Chinese [8] Kuhlemeyer R L,Lysmer J. Finite element accuracy for wave propagation problems[J]. J. Soil Mech. and Foundations Div.,ASCE, 1973,99SM5421–427. [9] Itasca Consulting Group,Inc. FLAC3DFast Lagrangian Analysis of Continua in 3D Version2.10,Users Manual[M]. USAItasca Consulting Group,Inc.,2002. [10] 李廷春,李术才,邱祥波. 三维快速拉格朗日法在安全顶板厚度 研究中的应用[J]. 岩土力学, 2004, 256 935–939.Li Tingchun, Li Shucai,Qiu Xiangbo. Application of fast Lagrangian analysis of continua to researching on safe rock covers of Xiamen subsea tunnel[J]. Rock and Soil Mechanics,2004,256935–939.in Chinese [11] 朱维申,李晓静,郭彦双,等. 地下大型洞室群稳定性的系统性 研究[J]. 岩石力学与工程学报,2004,23101 689–1 693.Zhu Weishen,Li Xiaojing,Guo Yanshuang,et al. Systematical study on stability of large underground caverns[J]. Chinese Journal of Rock Mechanics and Engineering,2004,23101 689–1 693.in Chinese [12] 刘春玲,祁生文,童立强,等. 利用 FLAC3D分析某边坡的稳定 性[J]. 岩石力学与工程学报,2004,23162 730–2 733.Liu Chunling,Qi Shengwen,Tong Liqiang,et al. Stability analysis of slope under earthquake with FLAC3D[J]. Chinese Journal of Rock Mechanics and Engineering,2004,23162 730–2 733.in Chinese [13] 董玉书,徐友宁. 厂坝铅锌矿区民采开发引发的问题思考[J]. 华 南地质与矿产, 2002, 4 23–27.Dong Yushu, Xu Youning. Some problems resulted from disordered mining of lead-zinc ores in Changba orefield[J]. Geology and Mineral Resource of South China, 2002,423–27.in Chinese [14] 刘松伟,尹贤刚,陈黎明,等. 厂坝铅锌矿 II矿体过渡层的工程 地质调查与评价[J]. 采矿技术,2003,341–4.Liu Songwei, Yin Xiangang,Chen Liming,et al. Engineering geology research and uation of transitional layer of Changba lead-zinc orefield No.II mine[J].Mining Technology,2003,341–4.in Chinese [15] 袁积余,陈双世. 群采破坏后矿山的技术改造与可持续发展[J]. 矿业快报, 2002, 3 10–12.Yuan Jiyu, Chen Shuangshi. Technical alteration and persistent development of mine after disordered mining destroy[J]. Express Ination of Mining Industry,2002,3 10–12.in Chinese [16] 尹贤刚,李庶林,唐海燕,等. 厂坝铅锌矿岩石物理力学性质测 试研究[J]. 矿业研究与开发,2003,23512–14.Yin Xiangang, Li Shulin,Tang Haiyan,et al. The tests and studies of physical mechanical properties of the rocks from Changba lead-zinc mine[J]. Mining Research and Development,2003,23512–14.in Chinese
展开阅读全文

资源标签

最新标签

长按识别或保存二维码,关注学链未来公众号

copyright@ 2019-2020“矿业文库”网

矿业文库合伙人QQ群 30735420