中钢集团
搜索 用户 菜单

(5积分下载)

如不能正常显示 ,请点击 此处 访问电脑版

2020-07-16
Series No. 528  
June 2020  
METAL MINE  
总第 528 期  
2020 年第 6 期  
厚大破碎矿体采场矿岩冒落规律及影响因素研究  
李晓斌1 勒治华22 李新成11  
1. 中冶沈勘秦皇岛工程设计研究总院有限公司,河北 秦皇岛066004;2. 东北大学资源与土木工程学院,  
辽宁 沈阳 110819)  
自然崩落法矿岩冒落是拉底卸荷、矿岩冒落、直至形成平衡拱的过程,其本质则是不断重复该过程。矿  
岩冒落规律受其地应力、采场结构参数、矿岩物理力学性质等地质因素和工程因素的影响,目前的研究多是假定为二  
维模型,侧重于采用理论方法或数值模拟分析冒落过程和冒落拱形态,较少涉及三维空间的冒落过程研究。根据矿岩  
3
D
冒落过程的主要特征参数,建立了描述该过程力学特征的本构模型,并采用 FLAC 实现数值模拟,模拟研究了矿体赋  
存深度、侧压力系数、拉底形态及采场尺寸参数对矿岩冒落规律的影响,研究结论有助于指导现场实践。  
关键词 地应力 采矿结构参数 冒落高度 冒落体积  
中图分类号 TD853.36  
文献标志码 A  
文章编号 1001-1250(2020)-06-172-05  
DOI 10.19614/j.cnki.jsks.202006027  
The Caving Process and Influential Factors in the Stope of Thick and Large Fractured Ore Body  
1
2
2
12  
Li Xiaobin Jin Zhihua Fan Shuo Li Xincheng  
1. Shenkan Qinhuangdao General Engineering Design&Research Institute CorporationMCCQinhuangdao 066004Chi?  
na2. School of Resources and Civil EngineeringNortheastern UniversityShenyang 110819China)  
Abstract Ore rock caving in the natural caving method is a process of undercutting and unloadingore-rock caving  
and even forming a balance arch. Its essence is to repeat the process continuously. The regularity of ore rock caving is affected  
by geological and engineering factors such as ground stressstope structural parametersphysical and mechanical properties  
of ore rock. At presentmost of the research is assumed to be two-dimensional modelfocusing on the theoretical method or  
numerical simulation analysis of the caving process and the shape of the caving archless involving the research of the caving  
process in three-dimensional space. According to the main characteristic parameters of the ore rock caving processa consti?  
3
D
tutive model is established to describe the mechanical characteristics of the processand FLAC is used to realize the numer?  
ical simulation on the effects of the ore body occurrence depthside pressure coefficientbottom form and stope size parame?  
ters on the regularity of ore rock caving. The research conclusions are helpful to guide the field practice.  
Keywords GeostressStructural parameters of stopeFalling heightFalling volume  
我国金属矿产资源的特点是规模较小、开采成  
本较高、市场核心竞争力弱。研究结果表明,大型金  
属矿山采矿成本的 30%~60%与岩体破碎有关。通过  
研究矿岩冒落的规律,调控采动应力诱导矿岩自然  
冒落,实现利用地压破碎岩体,可大幅度节约采准工  
程量与落矿费用,由此大幅度降低生产成本并提高  
生产效率,从而为国内金属矿山高效开采提供一条  
拟软件、理论分析等方法对不同构造应力场中采空区  
上覆岩层的破坏规律进行研究,结果表明随侧向应力  
的增加,上覆岩层塑性区高度和体积先减小后增大;  
[3]  
范广勤等 研究了在地应力较高的情况下矿岩冒落  
区域的轴比(或纵横比)与原始岩体的应力场的侧压  
力系数成正相关关系,研究发现,拉底的大小、形状和  
布置,以及切帮和预裂工程的存在,都直接影响岩体  
[1]  
[4]  
新的思路 。  
的崩落。徐腊明 利用三维有限元数值模拟程序  
近年来,人们逐渐认识到构造应力对采场上覆  
ADINA 模拟了程潮铁矿拉底和割帮工艺作用下矿体  
的应力分布,在此基础上,优选出合理的拉底方案。  
[2]  
3D  
岩层移动起重要作用。王军等 利用 FLAC 数值模  
收稿日期  
·
172 ·  
李晓斌等:厚大破碎矿体采场矿岩冒落规律及影响因素研究  
2020年第6期  
影响矿岩可冒性的因素有很多,主要包括矿岩  
的强度、岩体结构面特征、地应力的相对大小、采场  
结构参数等,本项目重点研究地应力与采场结构参  
数对矿岩冒落的影响。为实现借助拉底等工程因素  
诱导矿岩应力发生改变,利用地压破碎岩体的目的,  
就需要深入研究“破碎厚大矿体采场矿岩冒落规律”  
这一科学问题,这对于合理确定诱导冒落工程,充分  
利用地压,实现低成本、高效率采矿,具有重要的理  
论价值和实际意义。  
岩。实际上,利用自然崩落法开采矿岩体质量均稍  
差,在冒落过程中采场边帮围岩也会发生剪切破坏,  
但由于堆积在采场内的冒落矿岩散体支撑作用,往  
3
D
往导致剪切片帮破坏并不明显。但采用 FLAC 模拟  
时并不能考虑该作用,因此选用基于摩尔—库伦准  
则的理想弹塑性模型描述采场围岩,围岩可以发生  
[7]  
屈服,但并不降低承载力 。在两个准则中,拉伸准  
则优先。  
实际上,到目前为止,还没有一种数值模拟方  
法,可以完整地模拟矿岩冒落、堆积、并对采场边帮  
1
矿岩冒落数值模型  
3D  
3
D
起到支撑作用的全过程。本研究利用 FLAC 模拟矿  
岩冒落也是在一些假定条件下实现的。  
1
. 1 矿岩冒落的FLAC 系统  
3D  
FLAC 是基于连续介质力学理论的数值分析方  
(1)拉伸破坏准则  
法,虽然岩体中会产生塑性区,但岩体不会发生冒  
σ 3 σt,  
(1)  
落,不能模拟冒落的整个过程。为了模拟矿岩的冒  
3
D
式中,σ3σt 分别为岩体受到的拉应力与岩体的单轴  
落过程,引入了人为的假定。在 FLAC 数值模型中,  
如果数值网格的单元位于拉底空间的上部塑性区  
中,且该单元中存在拉应力,则认为该单元在重力的  
抗拉强度。  
2)摩尔—库伦破坏准则  
3
D
fs = σ1 -σ3 11 + sin? -2c 1 + sin?  
作用下可以冒落,利用 FLAC 的建模命令,该单元被  
杀死(“NULL”掉),被 NULL 掉单元的模型又处于非  
平衡状态,再次进行迭代计算,计算平衡后如果塑性  
区内依然有拉应力单元,则重复该过程,直至模型中  
不再出现拉应力区、平衡为止,计算停止。以此来近  
似模拟矿岩冒落过程。在计算过程中很难做到塑性  
区内完全不存在拉应力,所以在塑性区内基本无拉  
应力即可将计算停止。具体计算步骤如图 1。  
,
(2)  
-sin?  
1 -sin?  
式中,σ1σ3 分别是最大和最小主应力;c?分别为黏  
聚力和内摩擦角。当fs > 0 时,材料将发生剪切破坏。  
1
. 3 数值计算模型与参数  
由冒落拱理论可知,矿体埋藏深度、侧应力以及  
冒落拱跨度等因素是影响矿岩冒落高度的重要因素。  
3
D
基于上述 FLAC 的数值模拟方法和矿岩本构关系,  
对不同埋藏深度、不同侧压力系数、拉底形态及采场  
尺寸参数条件下的矿岩冒落规律进行了数值模拟研  
究。各种工况条件的数值计算模型略有不同,现以矿  
体埋藏深度对矿岩冒落的影响为例,简要介绍其数值  
计算模型以及矿、岩参数。建立三维数值模型如图 2  
所示,数值模型长 300 m,宽 150 m,高 350 m。模型共  
分 3 个部分,分组 1 为围岩部分、分组 2 为矿体部分,  
分组 3 为拉底部分。矿体部分长 100 m,宽 50 m,高  
1
00 m。矿体部分的底部距离模型底部 50 m,保证其  
不受边界效应的影响。岩体力学参数如表1。  
1
. 2 矿岩冒落本构关系  
研究表明,利用自然崩落法开采矿石时,矿岩冒  
落以拉伸破坏为主,因此采用拉伸准则描述冒落矿  
·
173 ·  
总第528期  
2020年第6期  
2
数值模拟方案及边界条件设计  
数值模拟分析了矿体不同埋藏深度、不同侧压  
力系数、拉底形态及采场尺寸参数 4种工况条件下矿  
岩冒落规律。各工况条件下的模拟方案、边界条件  
和拉底尺寸等如表 2。  
数值计算过程中矿岩冒落达到稳定状态判定:  
以矿体埋深 700 m计算分析为例,如图3所示,过拉底  
区域中心点,沿矿体长度方向的截面,可以看出,在  
将塑性区开挖 3次之后形成的冒落拱,拱线的上方基  
本无拉应力,可以判定冒落拱达到稳定状态。其他  
工况条件下矿岩冒落是否达到稳定状态判定完全相  
同。  
此,在实际生产中,对于厚大破碎矿体,小埋深反而  
有利于矿岩的冒落;矿体埋深较大时,可通过割帮等  
诱导工程使矿体水平应力减小,促进冒落拱的发育。  
同时矿体埋深较大时应注意冒落过程的片帮现象,  
制定合理的拉底范围,以免矿体开采过程中上下盘  
围岩的片帮,造成矿石贫化率增加。  
3
. 2 侧压力系数对采场矿岩冒落的影响  
不同侧压力系数下的矿体冒落部分如图 6所示,  
3
矿岩冒落规律分析  
图中下部为拉底区域,上部为冒落部分。发现侧压  
力系数为 0.4 时冒落拱发育最好,矿体冒落高度与冒  
落体积在5个方案中最大。  
3
. 1 矿体埋藏深度对采场矿岩冒落的影响  
各埋藏深度矿体拉底后矿体冒落部分如图 4 所  
示,其中下部区域为拉底部分,上部区域为塑性区部  
分。  
矿体的冒落高度和冒落体积随侧压力系数变化  
如图 7所示,矿岩的冒落高度和冒落体积的变化趋势  
基本一致。当侧应力系数从 0.4 增加到 1 时,矿体的  
冒落高度仅减小了 15.7%,但是矿体的冒落体积却减  
少了 43%,表现出在这一过程中冒落拱的发育变差。  
图 5 为矿岩冒落高度和冒落体积随矿体埋藏深  
度变化曲线。可知,矿岩的冒落高度和冒落体积随  
着埋藏深度的增加而减小,但减小幅度均较小。因  
·
174 ·  
李晓斌等:厚大破碎矿体采场矿岩冒落规律及影响因素研究  
2020年第6期  
变化曲线,矿体的冒落高度和冒落体积随拉底长宽  
比增大均呈现逐渐减小趋势。拉底形状为正方形时  
长宽比为 1),冒落高度和冒落体积达到最大,然而,  
在拉底长宽比为 2时,冒落体积有小幅度的增加。因  
此,在拉底设计时应考虑拉底形状为正方形(长宽比  
为 1),在地质条件或开采条件的约束下,拉底的长宽  
比最好不大于2。  
侧压力系数超过 1 后,冒落高度保持稳定,冒落体积  
出现缓慢增长。  
在实际生产中,如果矿体的侧应力较大,拉底完  
成后片帮问题会较为严重,应合理规划拉底的范围。  
因此,若使上赋矿岩冒落得更加充分,可以通过切帮  
等辅助工序来尽可能降低矿块的侧应力,从而达到  
使冒落高度和冒落体积增加的目的。  
3
. 4 采场尺寸参数对矿岩冒落的影响  
如图 10所示,随着拉底跨度的增加,采场冒落高  
度在前期增长较快,拉底的跨度接近拉底宽度时,冒  
落高度的增长趋势开始减缓,随着拉底跨度的持续  
增加,矿体的冒落高度保持稳定。拉底跨度较小时,  
不同宽度的矿体冒落高度较为接近,随着拉底跨度  
的增加,30 m 宽度矿体的冒落高度增长速度首先降  
低下来,50 m 和 70 m 宽度的矿体的冒落高度增长速  
率在拉底跨度不断增加的过程中也相继降低,最后  
都趋于稳定。图 11为不同矿体宽度的矿岩冒落高度  
与采场拉底跨度比值(高跨比)随拉底跨度变化曲  
线,从图中可以看出,拉底跨度相同时,高跨比随矿  
体宽度的增加而增加。相同宽度的采场,高跨比随  
拉底跨度的增加而减小。  
3
. 3 拉底长宽比对采场矿岩冒落的影响  
图 8为不同拉底长宽比条件下矿体冒落部分,下  
部为拉底区域,上部为上部塑性区。由塑性区分布  
图可知,拉底形状不同,各冒落拱的形状差异较大,  
冒落高度随着拉底长宽比的增大而逐渐减小。  
图 9 为矿体冒落高度和冒落体积随拉底长宽比  
4
结 论  
3D  
以 FLAC 数值模拟软件为工具,研究矿体在不  
同地质条件和结构参数条件下的冒落规律,利用所  
得到的冒落规律指导生产。主要结论如下:  
·
175 ·  
总第528期  
2020年第6期  
最后冒落高度趋于平衡。如果是厚大矿体,矿块尺  
寸设计时适当增大矿块的宽度以增加冒落高度和冒  
落体积。  
参 考 文 献  
1] 何荣兴,任凤玉,谭宝会,等 .论诱导冒落与自然崩落[J.金属矿  
山,20173:9-14.  
He RongxingRen FengyuTan Baohuiet al. Discussion on in?  
duced caving and block cavingJ.Metal Mine20173:9-14.  
2]  
军,宁建国,刘学生,等 . 构造应力下采空区上覆岩层破坏  
规律数值模拟[J. 煤矿安全,2017481:183-186.  
Wang JunNing JianguoLiu Xueshenget al. Numerical simula?  
tion of failure law of overlying strata in goaf under tectonic stress  
J. Safety in Coal Mines2017481:183-186.  
3] 范广勤,张千里 . 高地应力地下工程仿围岩自然崩落外廓特征  
的洞形结构模拟试验研究[C//15 届全国结构工程学术会议  
论文集(第册).焦作[: s.n.],2006 :377-384.  
Fan GuangqinZhang Qianli. Study on simulation test of cave-  
shaped structure of underground engineering with high in-situ  
stress and imitation of surrounding rock′s natural collapse outline  
C//Proceedings of the 15th National Conference on Structural En?  
gineeringThe 2nd book. Jiaozuo:s.n.],2006 :377-384.  
4] 徐腊明 . 自然崩落法拉底上部矿体应力分布有限元分析[J. 金  
属矿山,19962:11-15.  
Xu Laming. Analysls of orebody stress distribution over undercut  
dreft for natural block caving method by finite element methodJ.  
Metal Mine19962:11-15.  
(1)矿体在不同地质赋存条件下呈现不同的冒  
落规律。拉底条件相同,埋藏深度增加矿岩冒落高  
度,冒落体积略有减小,但总体变化不大;冒落高度  
随侧压力系数的增加呈现先减小后趋于平衡,冒落  
体积先减小后略有增加,侧应力系数为 1 时,冒落高  
度和冒落体积最小。  
5]  
利,陈亚娟,高 谦 . 拉底诱导的岩体崩落过程区数值模拟  
J. 有色金属:矿山部分,2009613:44-48.  
Wang LiChen YajuanGao Qian. Numerical on caving process  
zone induced by undeercutJ.Nonferrous Metals:Mining Section,  
2
009613:44-48.  
(2)拉底面积相同,拉底长宽比越大,矿体冒落  
6] 金川集团有限公司,长沙矿山研究院 . 自然崩落法矿岩崩落规  
高度和冒落体积越小。因此在设计拉底形状时,矿  
体拉底的长宽比应最好小于 2,保证拉底后有足够的  
冒落高度和冒落体积,以减小冒落辅助工程的工作  
量。  
律研究(研究报告之四)[R. 长沙:长沙矿山研究院,2006.  
Jinchuan Co. Ltd.Changsha Institute of Mining Research. Study on  
the Law of Natural Caving Ore and Rock CavingResearch Report  
4)[R. Changsha: Changsha Institute of Mining Research2006.  
7Hoek EBrown E T. Practical estimates of rock mass strengthJ.  
(3)采场尺寸确定后,冒落高度随拉底长度的增  
International Journal of Rock Mechanics Mining Sciences1997,  
加而增大,增长速度随拉底长度的增加而逐渐减小,  
3
48:1165-1186.  
·
176 ·  
骚小妹 电影

|

相关文章

导航 | 关于我们 | 联系客服