基于综合决策云模型的围岩稳定性分级方法研究
Research on Classification Method of Surrounding Rock Stability Based on Comprehensive Decision Cloud Model
通讯作者:
收稿日期: 2019-12-05 修回日期: 2020-03-12 网络出版日期: 2020-07-01
基金资助: |
|
Received: 2019-12-05 Revised: 2020-03-12 Online: 2020-07-01
作者简介 About authors
周科平(1964-),男,湖南衡阳人,教授,博士生导师,从事矿业工程方面的研究工作
关键词:
Keywords:
本文引用格式
周科平, 侯霄峰, 林允.
ZHOU Keping, HOU Xiaofeng, LIN Yun.
目前国内外关于工程围岩稳定性分级的方法有很多,从单指标的普氏系数法[3]、Terzaghi分级[4]、RQD法[5]到多指标的RMR法[6]、Q系统分类法[7]等,考虑的因素逐渐丰富,从定性评价转向半定量和定量评价,评价结果更加准确。近年来,一些新的理论和数学方法被运用到围岩稳定性评价中。贾超等[8]在综合分析物元理论、可拓集合论和关联函数运算理论的基础上,成功建立了洞室岩体质量评价的物元模型;盛继亮[9]建立了围岩稳定性分级的模糊综合评价模型,并将该模型应用于工程实际中,取得了较好的结果;Zhang等[10]基于15组工程案例建立了巷道围岩稳定性评价的距离判别模型。此外,还有学者将粗糙集理论、人工神经网络原理和支持向量机等方法引入稳定性评价中[11,12,13,14,15,16]。上述多数方法选取岩石质量指标、单轴饱和抗压强度及岩体完整性系数等因素作为评价指标,应用性能参差不齐,稳定性和普适性不足。以往方法多采用确定性理论分析,忽略了围岩稳定性分级的随机性和模糊性的影响。因此,需要引入新的评价方法,以提高围岩稳定性分级的准确性。
云理论是在不确定性人工智能理论的基础上发展起来的一种智能评价方法,能够适配围岩稳定性分级中存在的随机性和模糊性问题,解决上述方法的不足。因此,本文将云理论运用到围岩稳定性评价中,引入DEMATEL决策模型和熵权法组合确定各指标权重,进而建立围岩稳定性评价的云模型,并通过工程实例来验证模型的可行性和有效性,以期为地下工程围岩稳定性评价提供新思路。
1 围岩稳定性分级的综合决策云模型
1.1 云理论
云模型是李德毅院士基于概率测度空间提出的认知模型[17],该模型在不确定环境下能够实现定性概念与定量数据之间的相互转换,提供了数据与知识之间的双向认知通道。该模型已成功运用于行业竞争力评价、能源模式提取、生态安全评估、长短期变量预测、图像处理和岩土工程等众多领域。
(1)云模型概念。设U ={(x1,x2,…,xm)}是一个用于表示所研究范围的定量集合,称为论域。C是U上的定性概念。论域U中的任意元素x是具有稳定趋势的随机数,
云模型引入了期望
式中:
如果云分布满足:
则称在论域U上呈正态云模型分布。李德毅院士已经通过详细分析云分布和钟形隶属函数的普遍性论证了云模型的普适性[17],因此本文在围岩稳定性分级中采用云模型。
(2)X条件云发生器。云发生器[21]是云模型应用到实际的关键,主要包括正向云发生器和逆向云发生器。当给定论域同时包含3个数字特征(Ex,En,He)和特定值x时,云发生器称为X条件云发生器,生成云滴drop[x,μ(x)]。X条件云发生器的特点是输入相同的x,每次计算获得的确定度均略有变化,反映人们对于相同数据属于某个概念程度的不同观点,使得不确定性得到了较好的展现。由于围岩稳定性分级具有随机性和模糊性,且属于定性到定量的转换,因此本文采用X条件云发生器进行围岩稳定性分级。具体算法为
步骤1:确定3个数字特征(Ex,En,He)和特定数据x。
步骤2:产生1个以En为期望、He为方差的正态随机数
步骤3:迭代至生成N个云滴。
1.2 围岩稳定性评价指标的选取
围岩稳定性的影响因素较多,根据国内现行规范(GB/T 50218-2014工程岩体分级标准),并结合陈昌彦等[22]在岩体质量评价分级方法方面的研究,综合考虑地质、力学和工程3个方面的因素,选取了各种岩体质量分级方法中应用较为广泛的岩石单轴饱和抗压强度Rc、完整性系数Kv、岩石基本质量指标RQD、地下水影响修正系数K1和软弱结构面产状影响修正系数K2作为围岩稳定性分级指标。参照前人的研究成果[20],结合云模型的特点,将地下工程围岩划分为5个等级,即:Ⅰ级(稳定)、Ⅱ级(较稳定)、Ⅲ级(基本稳定)、Ⅳ级(不稳定)和Ⅴ级(极不稳定)。围岩稳定性级别与各评价指标的对应关系[23,24,25]如表1所示。
表1 围岩稳定性指标分级标准
Table 1
类别 | Rc/MPa | Kv | RQD/% | K1 | K2 |
---|---|---|---|---|---|
Ⅰ | >200 | 1.00~0.75 | 90~100 | 0~0.2 | 0~0.1 |
Ⅱ | 100~200 | 0.50~0.75 | 75~90 | 0.2~0.4 | 0.1~0.2 |
Ⅲ | 50~100 | 0.30~0.50 | 50~75 | 0.4~0.6 | 0.2~0.4 |
Ⅳ | 25~50 | 0.15~0.30 | 25~50 | 0.6~0.8 | 0.4~0.5 |
Ⅴ | 0~25 | 0~0.15 | 0~25 | 0.8~1.0 | >0.5 |
1.3 DEMATEL决策模型
DEMATEL(Decision Making Trial and Evaluation Laboratory,决策试验和评价实验法)决策模型是20世纪70年代由欧洲科学家提出的,能够在现实生活环境中研究解决复杂影响因素之间的关系。该方法可用于处理图论和矩阵论中的因素,分析各种影响因素互相影响关系矩阵,得到每个因素的影响度和被影响度,由此确定因素间的关联性,辨析出系统的主要因素,进而可以计算因素的权重[26]。目前,DEMATEL决策模型被广泛应用于现实情况下复杂因素之间的关系研究和权重确定,建模步骤如下:
Step 1:建立初始影响矩阵。通过专家学者打分确定各因素之间互相影响程度值,
Step 2:将初始影响矩阵归一化,得到直接关联矩阵X。
Step 3:计算综合影响矩阵T。
式中:tij表示指标i对指标j的影响度;I为同维度的单位矩阵,
Step 4:计算影响度fi和被影响度ej。
式中:fi称为影响度,表示因素i对其他因素的总影响值。ej称为被影响度,表示元素j受其他因素的总影响值。
Step 5:计算中心度Mj和原因度Rj。
中心度Mj表示指标j在系统中的重要程度,中心度越大表明该指标的重要程度越高。原因度Rj表示指标j与其他指标间的因果逻辑关系,Rj>0表示该指标对其他指标的影响大,j为原因因素,Rj<0表示该指标受其他指标的影响大,j为结果因素。
Step 6:确定权重。根据DEMATEL模型确定的中心度
1.4 基于熵权法的指标权重确定
熵权法是一种客观赋权方法[27]。客观评价方法能够很好地消除人为因素的影响,因此本文选用熵权法作为计算评价指标客观权重的方法,计算步骤如下:①根据评价对象和评价指标构造原始评价指标数据矩阵X=(xij)m*n;②对矩阵X进行归一化处理;③计算各指标的熵值
采用离差标准化法对原始数据进行归一化处理,公式如下:
各评价指标熵值和熵权的计算公式如下:
式中:m为评价对象的个数;n为每个评价对象中评价指标的个数;
1.5 组合权重综合决策模型
虽然DEMATEL决策模型在处理复杂情况下有优势,但难免受到主观因素的影响。因此,根据更为客观的熵权法得到的评价指标权重
根据各指标组合权重,计算确定度进行最终决策[20]。
2 综合决策云模型的建立与检验
2.1 模型建立
表2 标准化云模型数字特征
Table 2
类别 | Ex | En | He |
---|---|---|---|
Ⅰ | 0.1 | 0.0849 | 0.001 |
Ⅱ | 0.3 | 0.0849 | 0.001 |
Ⅲ | 0.5 | 0.0849 | 0.001 |
Ⅳ | 0.7 | 0.0849 | 0.001 |
Ⅴ | 0.9 | 0.0849 | 0.001 |
图1
2.2 模型检验
表3 隧道围岩分析数据
Table 3
样本 | Rc/MPa | Kv | RQD/% | K1 | K2 | 实际级别 |
---|---|---|---|---|---|---|
1 | 38.1 | 0.16 | 18.7 | 0 | 0.3 | Ⅴ |
2 | 37.1 | 0.18 | 27.4 | 0.2 | 0.1 | Ⅴ |
3 | 38.0 | 0.45 | 90.0 | 0.4 | 0.1 | Ⅳ |
4 | 38.0 | 0.35 | 91.0 | 0.2 | 0.2 | Ⅲ |
5 | 9.2 | 0.30 | 75.0 | 0.4 | 0.2 | Ⅴ |
6 | 17.4 | 0.34 | 38.4 | 0 | 0.2 | Ⅳ |
7 | 27.5 | 0.37 | 39.7 | 0 | 0.2 | Ⅳ |
8 | 56.8 | 0.50 | 58.9 | 0 | 0 | Ⅲ |
9 | 72.3 | 0.53 | 62.8 | 0 | 0 | Ⅲ |
10 | 67.5 | 0.45 | 58.9 | 0 | 0.2 | Ⅳ |
11 | 82.3 | 0.66 | 79.8 | 0 | 0 | Ⅱ |
12 | 56.1 | 0.29 | 41.5 | 0 | 0.2 | Ⅳ |
13 | 15.8 | 0.27 | 39.3 | 0 | 0.1 | Ⅴ |
14 | 49.2 | 0.19 | 29.8 | 0 | 0.6 | Ⅴ |
15 | 38.0 | 0.80 | 93.0 | 0.3 | 0.2 | Ⅲ |
16 | 21.4 | 0.17 | 25.9 | 0.3 | 0.1 | Ⅴ |
17 | 38.0 | 0.60 | 92.0 | 0.2 | 0.2 | Ⅳ |
18 | 65.9 | 0.58 | 69.3 | 0 | 0.2 | Ⅲ |
19 | 54.3 | 0.30 | 33.6 | 0.1 | 0 | Ⅳ |
20 | 47.1 | 0.25 | 27.5 | 0.2 | 0.3 | Ⅴ |
采用组合权重综合决策模型计算各评价指标的权重,得到各评价指标的主观权重
表4 各评价指标组合权重
Table 4
指标 | Rc/MPa | Kv | RQD/% | K1 | K2 |
---|---|---|---|---|---|
0.2066 | 0.2290 | 0.2255 | 0.1769 | 0.1620 | |
0.0858 | 0.0666 | 0.1572 | 0.4938 | 0.1966 | |
0.0945 | 0.0812 | 0.1889 | 0.4656 | 0.1698 |
计算得到20组隧道围岩工程数据的稳定性分级结果,并与现场围岩稳定性级别进行对照,如表5所示。
表5 围岩稳定性分级结果
Table 5
样本 | 综合确定度 | 本文结果 | 现场级别 | ||||
---|---|---|---|---|---|---|---|
1 | 0.2290 | 0.0100 | 0.2131 | 0.1041 | 0.2931 | Ⅴ | Ⅴ |
2 | 0.1233 | 0.0467 | 0.1882 | 0.1436 | 0.3777 | Ⅴ | Ⅴ |
3 | 0.2693 | 0.0452 | 0.2870 | 0.1433 | 0.0598 | Ⅲ | Ⅳ |
4 | 0.1337 | 0.1309 | 0.3354 | 0.2441 | 0.0181 | Ⅲ | Ⅲ |
5 | 0.0503 | 0.2055 | 0.0238 | 0.1507 | 0.2372 | Ⅴ | Ⅴ |
6 | 0.2346 | 0.1205 | 0.0358 | 0.4427 | 0.2371 | Ⅳ | Ⅳ |
7 | 0.2332 | 0.1205 | 0.0631 | 0.5833 | 0.0676 | Ⅳ | Ⅳ |
8 | 0.3142 | 0.2242 | 0.4467 | 0.0118 | 0 | Ⅲ | Ⅲ |
9 | 0.2979 | 0.2090 | 0.3678 | 0.0028 | 0 | Ⅲ | Ⅲ |
10 | 0.2802 | 0.2538 | 0.4109 | 0.0524 | 0.0001 | Ⅲ | Ⅳ |
11 | 0.3141 | 0.4237 | 0.0019 | 0 | 0 | Ⅱ | Ⅱ |
12 | 0.2348 | 0.2616 | 0.0828 | 0.3325 | 0.1251 | Ⅳ | Ⅳ |
13 | 0.3510 | 0.0372 | 0.0088 | 0.2874 | 0.3697 | Ⅴ | Ⅴ |
14 | 0.2295 | 0.0365 | 0.1534 | 0.0508 | 0.4360 | Ⅴ | Ⅴ |
15 | 0.1373 | 0.1220 | 0.2035 | 0.1963 | 0.0262 | Ⅲ | Ⅲ |
16 | 0.1237 | 0.0380 | 0.0023 | 0.1680 | 0.5104 | Ⅴ | Ⅴ |
17 | 0.1187 | 0.1668 | 0.2416 | 0.3035 | 0.0008 | Ⅳ | Ⅳ |
18 | 0.2767 | 0.6532 | 0.0874 | 0.0002 | 0 | Ⅱ | Ⅲ |
19 | 0.1820 | 0.2101 | 0.0720 | 0.2550 | 0.1920 | Ⅳ | Ⅳ |
20 | 0.0001 | 0.0362 | 0.4181 | 0.0891 | 0.4258 | Ⅴ | Ⅴ |
由表5可知,本文所建模型分级结果与现场级别吻合度达到85%以上,说明将综合决策云模型应用于围岩稳定性分级具有一定的可行性。同时,应用该模型可通过定量指标确定度反映围岩稳定性的随机性和不确定性,具有一定的优越性,在实际工程项目中具有一定的应用价值意义。
3 工程应用
表6 小红石砬子玉石矿工程实测数据和围岩分级结果
Table 6
序号 | Rc/MPa | Kv | RQD/% | K1 | K2 | 综合确定度 | 本文 结果 | KNN | RF | 现场 级别 | ||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
1 | 37.10 | 0.37 | 90.0 | 0.1 | 0.2 | 0.1308 | 0.3937 | 0.1778 | 0.0750 | 0.0795 | Ⅱ | Ⅳ | Ⅳ | Ⅱ |
2 | 56.10 | 0.45 | 60.0 | 0 | 0.1 | 0.3021 | 0.0453 | 0.3713 | 0.0177 | 0 | Ⅲ | Ⅲ | Ⅳ | Ⅲ |
3 | 23.28 | 0.19 | 25.9 | 0.3 | 0 | 0.0800 | 0.0006 | 0 | 0.0030 | 0.4129 | Ⅴ | Ⅴ | Ⅴ | Ⅴ |
4 | 54.30 | 0.34 | 36.8 | 0.1 | 0.1 | 0.1838 | 0.2580 | 0.3378 | 0.0821 | 0.1584 | Ⅲ | Ⅳ | Ⅳ | Ⅲ |
5 | 27.50 | 0.35 | 39.3 | 0.1 | 0.2 | 0.0063 | 0.3008 | 0.1175 | 0.1545 | 0.3730 | Ⅴ | Ⅳ | Ⅳ | Ⅳ |
由表6可知,综合决策云模型在小红石砬子玉石矿工程围岩分级中表现较好,准确率达到80%,优于K-最邻近算法和随机森林法的分级结果。研究表明,综合决策云模型在工程实际中具有一定的应用价值,对工程项目具有参考意义。
4 结论
(1)选取了各种岩体质量分级方法中应用较为广泛的岩石单轴饱和抗压强度Rc、完整性系数Kv、岩石基本质量指标RQD、地下水影响修正系数K1和软弱结构面产状影响修正系数K2作为围岩稳定性分级指标,并划定分级标准。
(2)考虑到围岩的复杂性、不确定性以及影响因素之间的互相影响关系,引入DEMATEL决策模型评估因素之间的影响关系并进行权重确定,并结合更为客观的熵权法对DEMATEL决策模型得到的指标权重进行修正,得到组合权重。
(3)基于云模型理论,运用潼荣高速公路曾家山隧道和西攀、攀田高速公路隧道工程的20组实测数据进行回检,结果表明本文所建模型分级结果与现场级别吻合度达到85%以上,与实际情况基本吻合,证明了评判模型的准确性和可靠性。
(4)运用围岩稳定性分级的综合决策云模型,对吉林省磐石市小红石砬子玉石矿巷道工程的围岩稳定性进行了分级,正判率达到80%,结果与现场结果吻合度较高,并优于K-最邻近算法和随机森林法,证明基于综合决策云模型的围岩稳定性分级方法为围岩稳定性分级提供了定量方法,具有一定的实际应用价值。
需要注意的是,本文所建立的综合决策云模型尚有不足之处,例如需要引入更多的评价指标,搜集更多的地下工程围岩稳定性分级的案例,从而有助于进一步完善模型,提高其适用性和泛化能力。
http://www.goldsci.ac.cn/article/2020/1005-2518/1005-2518-2020-28-3-372.shtml
参考文献
地下工程建设安全面临的挑战与对策
[J].
Challenges faced by underground projects construction safety and countermeasures
[J].
The index system of stability of surrounding rock based on PF-tree and cloud transform
[J].
铁路隧道围岩分级方法研究及发展
[J].
Research and development on the method of rockmasses classification in railway tunnel
[J].
粘弹性圆形巷道支护结构上的荷载及其围岩应力的解析解
[J].
Analytical solution of loads on supporting structur for circular tunnel and stress in viscoelastic surrounding rock
[J].
RQD在隧道围岩分类中的应用
[J].
Application of RQD in tunnel surrounding rock classification
[J].
RMR围岩分级法与中国公路隧道围岩分级方法对比
[J].
Comparison between RMR surrounding rock classification method and China highway tunnel surrounding rock classification method
[J].
特大采空区围岩稳定性评价
[C]//
Evaluation of surrounding rock stability in extra large goaf
[C]//
可拓学理论在洞室岩体质量评价中的应用
[J].
Application of extenics theory to evaluation of tunnel rock quality
[J].
地下工程围岩稳定性模糊综合评价模型研究
[J].
Research on fuzzy synthetic evaluation model for the stability of surrounding rock of underground project
[J].
Stability classification model of mine-lane surrounding rock based on distance discriminant analysis method
[J].
Deformation prediction model of surrounding rock based on GA-LSSVM-markov
[J].
Intelligent prediction of surrounding rock deformation of shallow buried highway tunnel and its engineering application
[J].
地下工程围岩稳定性分类的突变级数法研究
[J].
Catastrophe progression method for stability classification of underground engineering surrounding rock
[J].
基于RS-TOPSIS法的地下工程岩体质量评价
[J].
Rock mass quality evaluation of underground engineering based on RS-TO-PSIS method
[J].
围岩稳定性分类的BP网络模型研究
[D].
Study of the BP Neural Network on the Stability Classification of Surrounding Rocks
[D].
支持向量机在围岩稳定性分类中的应用
[J].
Application of support vector machine in classification of surrounding rock stability
[J].
不确定性人工智能
[J].
Artificial intelligence with uncertainty
[J].
基于RS-云模型的硫化矿石自燃倾向性综合评价
[J].
Comprehensive evaluation on spontaneous combustion tendency of sulfide ore based on RS-cloud model
[J].
隶属云和隶属云发生器
[J].
Membership clouds and membership cloud generators
[J].
正态云模型在煤巷围岩稳定性分类中的应用
[J].
Application of normal cloud model in surrounding rock classification of coal roadway
[J].
基于数据场和云模型的洪水灾害风险等级评估
[J].
Flood risk evaluation methods based on date field and cloud model
[J].
各类岩体质量评价方法的相关性探讨
[J].
Discussion on the interrelation of various rock mass quality classification systems at home and abroad
[J].
基于云模型及粗糙集理论的围岩稳定性分级方法研究
[J].
Stability evaluation for the surrounding rock structure based on the normal cloud model and the rough set theory
[J].
潼荣高速公路曾家山隧道围岩分级及稳定性评价
[D].
The Classification and Stability Evaluation of the Surrounding Rock of Chongqing Tongrong High-speed Zengjiashan Tunnel
[D].
基于熵权—云模型的隧道围岩分级方法研究
[J].
Classification method for tunnel surrounding rock based on the entropy-cloud model
[J].
基于因子分析-DEMATEL定权法的期刊综合评价研究
[J].
Research on journal comprehensive evaluation based on factor analysis and DEMATEL weight method
[J].
熵权法在重大危险源应急救援评估中的应用
[J].
Entropy method for major hazards emergency rescue
[J].
主客观权重的组合方式及其合理性研究
[J].
Research on a combined method of subjective-objective weighing and the its rationality
[J].
山区高速公路隧道围岩分级方法及应用研究
[D].
The Methods of Surrounding Rock Classification and Their Application in the Expressway Tunnel of Mountainous Area
[D].
A comparative study of K-nearest neighbour,support vector machine and multi-layer perceptron for thalassemia screening
[J].
Mining data with random forests:A survey and results of new tests
[J].
/
〈 |
|
〉 |
