2024年5月17日发(作者:电脑硬盘修复多少钱)
2023年第3期/第44卷
矿业工程
其次为方案3力最大,,方案1的最大主应力最小。
2)数值模拟计算结果表明,主要岩层移动位置
出现在采空区的上方,越靠近采空区,围岩的变形值
越大;反之,距离采空区距离越远,围岩的变形值越小;
且上盘水平位移值大于下盘的水平位移值,位移云图
呈不对称的形状;方案1、方案2和方案3的地表最大
沉降值分别为124.16mm、260.47mm和184.94mm,
地表最大水平位移分别为29.74mm、46.14mm和
65.08mm。
3)当采用分段空场嗣后充填采矿法只开采
-100~-400m矿体,或采用上向水平分层充填采
40~-100m矿体进行开采时,导水裂隙带矿法对-
距离地表高度大于100m,地表河流处于相对稳定状
态;当采用分段空场嗣后充填采矿法对-40~
塑性区范围显著增加,导水-100m矿体进行开采,
裂隙带距离地表距离相对较近,因此,建议-40~
-100m矿体不开采或采用上向水平分层充填采矿
法进行开采。
4)工程应用结果表明,目前矿山开采过程中,地
表仍处于相对稳定的状态,但后期应加强监测,实时
对地表的安全性进行分析。
[参考文献]
[1] 李云涛,付煜,任海龙,等.眼前山铁矿露天转地下开采边坡破
J].采矿技术,2022,22(2):61-70.坏过程及破坏机理研究[
[2] 牛淑慧.迤纳厂铁铜矿采空区及地表稳定性现状评价分析[J].
昆明冶金高等专科学校学报,2021,37(3):50-53.
[3] 张康顺.瓮福磷矿穿岩洞矿段地下开采对边坡稳定性的影响研
D].贵阳:贵州大学,2021.究[
[4] 逄铭璋,张爱民,任鹏召,等.某矿深部开采对地表及竖井稳定
性影响研究[J].中国矿山工程,2021,50(6):16-20.
[5] 张超,宋卫东,付建新,等.深部复杂金矿体充填开采对地表建
筑物稳定性影响[J].东北大学学报(自然科学版),2021,42
(8):1143-1151,1158.
[6] 刘明宇,王荣林,刘发平.某铁矿地下开采对地表设施的稳定性
影响分析[J].现代矿业,2021,37(7):254-257.
[7] 丁飞,付俊,周罕,等.河流下开采岩层移动及地表变形预测研
究[J].矿冶工程,2019,39(5):14-16,21.
[8] 刘仁冬,解联库,于世波.某硫铁矿地下开采地表变形预测及影
响评价[J].有色金属(矿山部分),2013,65(1):40-43.
[9] 吕文玉.河流下条带充填开采可行性分析[J].煤矿安全,2013,44
(11):201-203.
[10] 陈嘉生.水域动载荷条件下复杂矿体开采安全技术[D].长沙:
中南大学,2010.
[11] 闫照存.黄土沟壑区多工作面开采地表形变破坏分析与评价[D].
西安:西安科技大学,2021.
[12] 陈玉平,邓喀中.皖北矿区河流下煤炭安全开采技术[J].煤矿
安全,2012,43(7):182-184.
[13] 崔景文.地表河流水体煤柱回采可行性研究及开采措施[J].
煤炭技术,2009,28(6):70-72.
[14] 宋霁洪,于灯凯,孙洋,等.开采扰动下三圣铁矿渗流特性及涌
水风险分析[J].黄金,2022,43(11):43-46.
Researchontheinfluenceofundergroundmining
ofamineonthesafetyoftheriversonthesurface
XuZhao
牗ShenzhenShiheSafeTechnologyCounselingCo.牞Ltd.牘
Abstract牶Inacasestudyofoneundergroundmine牞basedonthecomprehensiveanalysisofthegeologicalcondi
tionsandminingconditions牞thetypicalexplorationlinesectionistakenastheresearchsubject.Anumericalmodelis
3D
establishedaccordingtothefinitedifferencesoftwareFlac牞andtheinfluenceofundergroundminingwithdifferent
schemesonthesafetyoftheriversonthesurfaceisstudied.Theresultsareappliedonsiteandverified.Theresearch
resultsshowthattheorebodiesburiedbetween-100mand-400mcanbeextractedbysublevelopenstopeand
subsequentminingmethodsforsafemining牞andtheorebodiesbetween-40mand-100mcanbeextractedbythe
cutandfillmethodorsecuritypillarreservationminingmethod.Atpresent牞thestabilityofriversonthesurfaceis
good.Lateron牞themonitoringofthesafetyoftheriversonthesurfaceshouldbeenhanced.Theresearchresultscan
provideatheoreticalbasisforthesafeandefficientminingofsimilarmines牞alsoitisofsignificanceinprotectingthe
safetyoftheriversonthesurface.
Keywords牶fillingminingmethod牷riversonthesurface牷numericalsimulation牷plasticzone牷safety
Copyright©博看网. All Rights Reserved.
矿业工程
黄 金
GOLD
2023年第3期/第44卷
基于模糊数学的岩体力学参数折减方法
何斌全,李仲泽,陈国辉,陈 功
(湖南柿竹园有色金属有限责任公司)
摘要:为评价柿竹园多金属矿地下露天协同开采后的边坡稳定性,需要全面了解边坡的地质情
况与岩石力学参数。通过开展单轴抗压强度试验、劈裂拉伸试验、剪切试验,测量出标准试件力学
特性参数,包括抗压强度、抗拉强度、剪切强度及内聚力、内摩擦角、弹性模量等。针对传统岩体参
数折算的主观性和模糊性问题,采用基于室内岩石样本试验和工程地质资料相结合的模糊数学法
进行折减,得到了可用于实际工程分析和建模计算的较为合理、准确的岩石力学参数。
关键词: 露天开采;协同开采;边坡稳定性;模糊数学;岩体力学参数;折减
中图分类号:TD679
文献标志码:A
文章编号:1001-1277(2023)03-0026-05
doi:10.11792/hj20230306
引 言
湖南柿竹园有色金属有限责任公司(下称“柿竹
园多金属矿”)是一座以钨、钼、铋、锡等战略资源为
主的特大型矿山,拥有非均衡品位的多金属矿石资源
储量3.90亿t,是中国重要的有色金属资源基地。面
对历史遗留的采空区群高品位矿柱资源和采空区群
周边低品位资源难以回采的难题,该矿山提出了高危
1]
。考虑到采空区群下地下露天协同开采技术方案
[
露天开采爆破频率高、一次爆破量大、爆破振动作用
强、扰动范围广,日积月累的爆破振动可能会引发露
天开采境界外的高陡边坡失稳垮塌,进而诱发滚石或
滑坡灾害阻塞甘港河河流、阻断县道,严重危害车辆
与行人安全。因此,亟须对其边坡稳定性进行研究。
在对边坡稳定性进行研究时,不管是采用数值模
拟,还是理论分析计算方法,都取决于岩体力学参数
2-3]
选取的可靠性与合理性
[
。由此可见,矿岩物理力
学性质的研究是边坡稳定性研究的基础。为了实现
柿竹园多金属矿后续大规模露天开采的安全,需要对
露天采坑境界和高陡边坡进行稳定性分析与安全性
评价,并提出针对性的高陡边坡安全治理措施。然
而,该矿山近期没有专题的岩石力学报告,难以开展
露天采坑境界外高陡边坡稳定性分析与安全性评价。
力学参数的确定是岩体力学工程界的一大难
4-5]
。岩体力学参数的选取决定了数值模拟、工程题
[
6]
。当前,岩体力学参数的设计的准确性与可靠性
[
7]8]
确定方法主要有工程类比法
[
、原位试验法
[
和室
9-10]
内试验法
[
等。其中,工程类比法基于相似工程
不足;原位试验法最直接、最可靠,但周期长、成本高。
因此,本次研究基于室内岩石样本试验和工程地质资
料相结合的模糊数学法进行折减,对露天开采最终境
界工程条件进行深入和全面的了解,为安全高效开采
提供理论支持。
1 矿区地质概况
矿区出露地层主要为泥盆系中—上统棋梓桥组
及佘田桥组碳酸盐岩。此外,矿区东部小面积出露震
旦系浅变质碎屑岩及泥盆系中统跳马涧组石英砂岩。
各时代地层自老至新分别为:震旦系下组、泥盆系中
统跳马涧组、泥盆系中—上统棋梓桥组及佘田桥组。
根据矿区岩体性质、矿化类型、矿化强度、矿物组合、
结构构造及有用元素含量,矿石的工业类型划分为
Ⅰ
、
Ⅱ
、
Ⅲ
、
Ⅳ
4种类型:
Ⅰ
类矿石为网脉状大理岩矿
石,
Ⅱ
类矿石为矽卡岩钨铋矿石,
Ⅲ
类矿石为花岗岩
网脉矿岩,露
Ⅳ
类矿石为云英岩钨钼铋矿石。其中,
天边坡岩体主要为大理岩、矽卡岩、花岗岩。
2 边坡岩体力学参数确定
为了研究该矿山边坡岩体的力学参数,本次试验
前期先在矿山取矽卡岩、大理岩和花岗岩岩样进行加
工打磨,制作标准试件(见图1)。本次现场采取的岩
样来自西南工作面及矿石堆中的大块。在岩样选取
过程中,遵守以下4项原则:
1)岩石类型具有代表性(矿体、上下盘围岩、地
质构造带等)。
2)岩样符合空间分布的规律性(垂直方向、水平
方向、岩层分布、强度变化规律)。对参数进行取值,精度差且受限于工程经验,可靠性
收稿日期:2022-11-08;修回日期:2023-01-10
基金项目:国家自然科学基金项目(41972283)
作者简介:何斌全(1965—),男,正高级工程师,硕士,从事矿山生产管理和智能矿山建设等采矿相关的工作;Email:2859560277@qq.com
Copyright©博看网. All Rights Reserved.
2023年第3期/第44卷
剪切试验。
2.1 抗压强度试验
矿业工程
1)试验原理。当试样在轴向压力作用下出现压
缩破坏时,单位面积上所承受的载荷称为岩石的单轴
抗压强度,即试样破坏时的最大载荷与垂直于加载方
向的截面积之比。
图1 岩体力学试验试样
3)尽量避免爆破开挖等对岩样的损伤与破坏。
4)所取钻孔岩芯应尽量靠近矿体,大小应保证
符合加工尺寸要求。
采用实验室用英国INSTRON公司的电液伺服材
料控制机及大型压力试验机对试件进行抗压、抗拉、
)计算公式如下:抗压强度(
σ
c
p
σ
c
=
A
2
(m)。
(1)
式中:p为试件的破坏载荷(kN);A为试件截面积
2)试验结果。3种岩石单轴压缩应力-应变曲
,单轴抗压强度试验结果见表1。线见图2
图2 3种岩石单轴压缩应力-应变曲线
表1 3种岩石单轴抗压强度试验结果
编号
Mc1
Mc2
Mc3
平均
Gc1
Gc2
Gc3
平均
Sc1
Sc2
Sc3
Sc4
Sc5
平均
峰值载荷/kN抗压强度/MPa弹性模量/GPa
126.96
134.51
110.80
124.09
333.03
329.14
334.86
332.34
339.52
355.77
379.54
270.28
401.43
349.31
65.62
68.81
56.96
63.80
171.12
169.79
172.06
170.99
175.29
183.46
195.79
139.09
206.76
180.08
44.67
44.01
36.01
41.56
52.91
55.04
55.85
54.60
68.08
62.02
66.29
70.12
62.65
65.83
泊松比
0.207
0.291
0.306
0.268
0.1621
0.2669
0.2853
0.2381
0.2667
0.2075
0.2140
0.2218
0.2092
0.2238
发生拉伸破坏时的载荷计算岩石的抗拉强度。假定
岩石是均质、各向同性的线弹性体,根据岩石试样受
力分析,进行拉应力计算。
1)试验原理。抗拉强度()计算公式为:
σ
t
2p
σ
t
=
Dh
π
(2)
式中:D为试样的直径(mm);h为试样的高度(mm)。
2)试验结果。3种岩石劈裂力-位移曲线见
,抗拉强度试验结果见表2。图3
2.3 剪切强度试验
1)试验原理。室内岩石剪切强度测定最常用的
是测定岩石的抗剪断强度。本次试验剪切角度为
50°、60°、70°。按式(3)可求得作用于剪切面上总法
向载荷(N)和总剪切载荷(Q)。
N=p(cosfsin
α
)
α
+
Q=p(sin
α
+fcos)
α
的摩擦系数。
由式(4)可以求得作用于剪切面上的法向应力
()和剪应力()。
στ
Np
cosfsin
α
)
σ
==(
α
+
AA
Qp
==(sin
α
+fcos)
τα
AA
2.2 抗拉强度试验
测定岩石抗拉强度的方法较多,有直接拉伸法、
劈裂拉伸试验法、弯曲试验法、离心机法、圆柱体或球
体的径向压裂法等。其中,以劈裂拉伸试验法最为简
易,其试样制作简单。劈裂拉伸试验法是将岩石加工
成圆盘形试样,沿试样轴面平行粘贴2根合金钢丝,
然后将试样置于试验机上平行该轴面加压,借助合金
钢丝将集中载荷转变为线载荷,从而产生垂直于该轴
面的拉应力,最后导致试样发生拉伸破坏,根据试样
}
(3)
式中:°);f为圆柱形滚子与上下压板
α
为剪切角度(
}
4)(
Copyright©博看网. All Rights Reserved.
矿业工程
黄 金
图3 3种岩石劈裂力-位移曲线
表2 岩石抗拉强度试验结果
编号
Mt1
Mt2
Mt3
Mt4
Mt5
Mt6
平均
Gt1
Gt2
Gt3
Gt4
Gt5
Gt6
平均
St1
St2
St3
St4
St5
St6
平均
直径/mm
49.59
49.56
49.57
49.58
49.53
49.56
49.57
49.66
49.68
49.62
49.68
49.59
49.72
49.66
49.63
49.69
49.66
49.70
49.55
49.73
49.66
高度/mm
25.66
25.66
25.51
25.62
25.52
25.61
25.60
25.65
25.92
25.67
25.73
25.67
25.11
25.63
25.67
25.65
25.57
25.68
25.64
25.65
25.64
峰值载荷/kN
12.91
8.311
9.83
9.97
5.34
5.61
8.662
18.67
16.69
19.56
20.80
12.20
14.87
17.13
19.94
24.88
17.41
24.71
29.46
25.42
23.64
抗拉强度/MPa
6.45
4.16
4.95
5.00
2.69
2.81
4.34
9.33
8.25
9.78
10.36
6.22
7.58
8.59
9.96
12.43
8.72
12.16
14.76
12.69
11.79
试验中摩擦系数f可忽略不计,根据以上公式计
算可得到试件在不同剪切角度作用下的剪应力值和
法向应力值。然后根据摩尔-库仑定律(
τ
=
C+tan
φ
)作图,利用EXCEL办公软件自动线性回
σ
归求出岩块的内聚力(
C)和内摩擦角(),根据试验
φ
-得到
τσ
关系曲线。
2)试验结果。3种岩石剪切回归曲线见图4,剪
切强度试验结果见表3。
2.4 岩石试样物理力学性质
柿竹园多金属矿岩石物理力学性质见表4。
3 岩体力学参数折减
矿区实际岩石力学研究中包括4种岩石材料:矽
卡岩、花岗岩、花岗斑岩和大理岩。岩体主要由结构
面和结构体组成,岩体在结构面切割弱化的作用下,
其力学参数与完整岩体的力学参数有较大差别,考虑
到岩石力学试验中所用的试件与工程中岩体的差别,
故室内岩石力学试验测得的岩石物理力学参数不能
直接用于工程岩体计算中,需要对试验数据进行折
减。
图4 3种岩石剪切回归曲线
岩体力学参数的折减方法很多,利用经验公式确
定岩体力学参数是最常用的方法。但是,经验公式是
以一定数量的室内试验和现场试验资料为依据,通过
回归等手段求出的,未能把较多的地质描述引入其
中,各个经验公式计算同一岩体时,离散性很大。而
影响岩体力学的因素很多,各个因素的影响程度又不
一样,所以岩体强度和岩石强度的关系不能确定,即
很难确定其折减程度。
所以,对于岩体力学参数的工程处理可以引进模
糊数学对岩石强度进行模糊折减,也就是要经过模糊
评价获得一个折减系数(),使得:SS,其中S
R
=
rR
、
S
r
分别为岩石和岩体强度指标。岩石力学参数折减
系数的模糊计算步骤如下:
把所有因素的全体记为因素论域(F)。
F={f,f,…,f}
12n
建立在该论域上的因素模糊集(A)记为:
A=afaf…+af 0 (i=1,2,…,n)
≤
a
≤
1
1
/
1
+
2
/
2
+
n
/
ni
(6)
(5)
Copyright©博看网. All Rights Reserved.
2023年第3期/第44卷
矿业工程
Ms1
Ms2
Ms3
Gs1
Gs2
Ss1
Ss2
Ss3
Ss4
Ss5
Ss6
Ss7
Ss8
表3 3种岩石剪切强度试验结果
50
60
70
50
70
50
50
60
60
60
70
70
70
31.8
11.63
6.01
100.68
13.38
132.29
108.74
45.81
43.35
62.2
12.35
17.76
12.58
37.89
20.16
16.51
119.99
36.76
157.66
129.59
79.34
75.00
107.61
33.93
48.80
34.56
式中:f的隶属度,它是因素f
i
为对A
i
在评价因素中
所起作用的大小量度,也就是因素评定折减系数的能
力,具有权重的意义。将供选择的折减系数的全体记
为备择论域(V):
V={,,…,}
12n
为:
B=bb…+b
1
/
1
+
2
/
2
+
n
/
n
属度。
则权向量为:
A=(a,a,…,a)
12n
理论中提出的1~9标度法,见表5。
(9)
其中权向量的确定是采用SaatyT.L在层次分析
(8)
b式中:
i
为折减系数
i
对综合评定模糊子集B的隶
(7)
B,评判别集)建立在该论域上的备择模糊子集(
编号峰值载荷/kN剪切角度/(°)法向应力/MPa剪应力/MPa
123.41
58.203
43.45
386.56
95.10
514.85
420.20
228.39
217.52
310.53
90.11
128.98
.1491
表4 岩石物理力学性质
类别
矽卡岩
大理岩
花岗岩
抗压强度/MPa
180.08
63.80
170.99
抗拉强度/MPa
11.79
4.34
8.59
弹性模量/GPa
65.83
41.56
54.60
泊松比
0.2238
0.268
0.2381
内摩擦角/(°)
44.77
40.15
43.63
内聚力/MPa
29.45
10.97
24.04
-1
波速/(m·s)
4966.2
3796.75
3890.67
表5 比率标度及其定义
标度
1
3
5
7
9
含义
表示两个因素相比,具有同样重要性
表示两个因素相比,一个因素比另一个因素稍微重要
表示两个因素相比,一个因素比另一个因素明显重要
表示两个因素相比,一个因素比另一个因素强烈重要
表示两个因素相比,一个因素比另一个因素极其重要
因素i与因素j比较得判断a,则因素j与因素i比较
ij
倒数
1
得判断a
ji
=
a
ij
比率标度
f
1
f
2
f
3
f
4
表6 因素权向量判断矩阵
f
1
1
0.333
0.2
0.143
f
2
3
1
0.333
0.2
f
3
5
3
1
0.333
f
4
7
5
3
1
2,4,6,8 上述两相邻判断的中值
通过有关文献估计确定岩体抗压强度的折减系
0.5,0.4,0.3,0.2,0.1}(见表7)。数的备择集为{
表7 f
1
评价的判断矩阵
比率
标度
0.5
1
2
4
5
7
0.4
0.5
1
2
4
5
折减系数
0.3
0.25
0.5
1
3
4
0.2
0.2
0.25
0.333
1
2
0.1
0.143
0.2
0.25
0.5
1
根据矩阵定律求出权矩阵的特征值和特征向量
A=(a,a,…,a),即为所需的权向量;用类似的方
12n
法可以求出折减系数
i
对综合评定模糊子集B的隶
属度b;最后利用隶属度作为权的加权平均值来确定
i
折减系数的综合评价值,即:
n
0.5
0.4
0.3
0.2
.10
i
=
b
∑
ii
1i=
n
若要计算出以上所有矩阵的特征向量,手工计算
(10)
量太大,既费时又费力,所以本次计算借用Matlab软
件里的[E,]=eig(a)进行计算。其中,a为矩阵的
λ
名称,计算时只要输入以上数据,即可很快得到计算
结果;并取特征值为
λ
。其中,n为矩阵的维数。
≈
n
T
经计算,因素集的权向量为
Α
=(0.8880,0.4120,
0.18470,0.0869)。
折减系数
i
对综合评定模糊子集B的隶属度,
即:
b
∑
i
1i=
9标度法,对影响岩石力学首先根据表5的1~
性质的4个因素:岩体结构面(f)、岩体结构完整程
1
度(f)、地下水(f)、风化程度(f),求出对岩体的抗
234
剪强度的权向量。因素权向量的判断矩阵见表6。
根据以上影响因素的分析,岩体结构面在岩体力学性
质评价中起着较为重要的作用。
Copyright©博看网. All Rights Reserved.
发布者:admin,转转请注明出处:http://www.yc00.com/num/1715931673a2693842.html
评论列表(0条)