基于正态隶属度—属性区间识别模型的岩爆倾向等级预测
1.
2.
3.
Propensity Grade Prediction of Rockburst Based on Normal Membership-Attribute Interval Recognition Model
1.
2.
3.
通讯作者:
收稿日期: 2021-03-05 修回日期: 2021-08-29
基金资助: |
|
Received: 2021-03-05 Revised: 2021-08-29
作者简介 About authors
贡力(1977-),男,江苏丹阳人,教授,博士生导师,从事深埋隧道及输水工程方面的研究工作
关键词:
Keywords:
本文引用格式
贡力, 陆丽丽, 靳春玲, 梁栋, 周汉国, 谢平.
GONG Li, LU Lili, JIN Chunling, LIANG Dong, ZHOU Hanguo, XIE Ping.
岩爆是一种发生在深埋隧洞中的常见地质灾害,是隧洞开挖工程中由于岩体本身积聚的应变能量超过岩体破坏时消耗的能量,岩石受到应力重分布作用将多余能量释放出来致使岩石发生爆裂、剥离的现象(李任豪等,2020)。岩爆是高地应力地区常见的地质现象,最早发生在矿井中,后来随着大量深埋隧洞修建于高地应力区域,岩爆问题受到了矿山和隧道工程领域研究人员的重视。为了降低岩爆灾害对作业人员及机械的损害程度,探寻科学、合理的岩爆倾向等级预测方法对岩爆防控具有极其重要的意义。
多年来,众多学者对岩爆预测展开了诸多研究。早期岩爆预测研究只考虑单一因素的影响,如:岩石刚度、强度和能量等,且国内外学者提出了不同的判据,如:Russense判据(张镜剑等,2008)和Kidybinski判据(Kidybinski,1981)等。然而,岩爆发生受多种因素共同制约,单一判据具有片面性。随着岩爆灾害研究工作的深入,岩爆发生的影响因素与岩爆倾向等级的非线性关系逐渐引起学者们的关注,新的数学方法及智能算法模型开始应用到岩爆预测过程中。数学方法主要有模糊综合评判法(Adoko et al.,2013;王元汉等,1998)、物元可拓理论(张乐文等,2010)、灰色系统理论(裴启涛等,2013)、功效系数法(王迎超等,2014)和云模型(王迎超等,2015)等。智能算法主要有BP神经网络(孙臣生,2019)、人工神经网络(张光存等,2013)、深度神经网络(田睿等,2020)和随机森林方法(Dong et al.,2013)等。
上述方法在岩爆倾向等级预测方面取得了一定的成果,但是由于岩爆预测的复杂性,这些方法在实际应用中仍存在2个问题:(1)各种方法自身具有一定的局限性,例如:模糊综合评价法在确定指标权重时,人为因素的影响较多,评判结果不够客观;BP神经网络学习速率的选择缺乏有效的方法且不能保证收敛到全局最小点;随机森林法对于小样本数据可信度不高。(2)运用非线性关系评价方法的核心是指标权重的计算,其中应用较广的是采用主观权重与客观权重相结合来消除单一权重的影响,但其核心是针对评价指标体系的构建,忽略了指标实测值的改变以及其余指标对权重的影响,导致同一指标体系下的不同工程实例和不同指标实测值所赋权重均相同。
岩爆倾向等级预测属于多属性有序分割类问题,通过分析现有数学方法和智能算法,发现智能算法应用简单,计算过程全部由软件驱动,但是部分智能算法由于初始化随机权重使得计算结果不稳定。而数学方法可以避免该缺点,其中属性区间识别模型在解决多属性有序分割类问题上具有明显的优越性(殷欣等,2020)。因此,本文引入属性区间识别模型进行岩爆倾向等级预测。为了避免权重计算出现前述问题,本文采用正态隶属度函数计算指标概率密度,并利用Jousselme公式计算指标间距离,进而确定指标权重。该方法从各指标所属各等级的距离出发,综合考虑各指标各等级间的相互影响程度,对同一指标体系下不同工程实例赋予不同权重值,使得权重计算相对更加合理。最后,应用13个典型深埋引水隧洞案例进行模型准确性测试,并以双江口水电站SPD9厂房为例进行了工程实例验证,结果表明该模型预测结果能够客观反映岩爆实际结果,可为预测类似深埋引水隧洞岩爆倾向等级提供新思路。
1 属性区间识别模型
设D为全体待评价隧洞,Di ∈D(i=1,2,…,m)为若干典型的待评价隧洞,每个隧洞有n个评价指标I1,I2,…,In,xij 为第i个标段的第j个评价指标的测量值。假定G为D上的某类属性空间,评价集(E1,E2,…,Ek )为属性测度区间空间G的有序分割类,根据已知的分级标准构造标准的分级矩阵为
式中:
1.1 单指标属性测度区间的计算
将xij 具有属性Ep (1≤p≤k)的属性测度区间记为
当
当
1.2 综合属性测度的计算
将均化系数(邹强等,2011)引入属性测度区间,计算各指标的综合属性测度
式中:ε为均化系数,ε∈(0,1)。
1.3 待评价隧洞等级识别
通过综合属性测度计算第i个评价标段属于等级Ep 的综合属性测度μip,表示为
式中:wj 为第j个评价指标的权重。
按照置信度准则,取置信度λ(0.5≤λ≤1),若p满足
2 权重的确定
2.1 正态隶属度函数
评价对象的性质决定评价指标,模糊理论的关键是建立指标隶属度函数(林大涌等,2018)。正态隶属度函数具有模糊性,且能够减少主观因素的影响,因此本文采用正态隶属度函数来确定岩爆评价指标属于某评价等级的隶属度,具体公式如下:
式中:
各评价等级指标边界值与相邻2个评价等级隶属度相同,因此计算某评价等级端点值时,
根据式(13)~
2.2 指标权重的确定
评价指标的权重值是度量指标重要性的方法。考虑各指标间以及各等级间的相互影响,采用指标间的相互距离计算指标权重。正态隶属度函数确定的指标隶属度又称为指标概率密度,由概率密度计算指标距离进而确定权重,计算步骤如下:
设Ii 和Ij 为指标集合中的2个评价指标,则其对应的各等级概率密度为mi 和mj,计算指标Ii 和Ij 之间的Jousselme距离dij (Shen et al.,2019),表示为
式中:mi 和mj
由指标间的距离可以计算指标Ii 和Ij 的信任度,表示为
指标Ii 相对于指标Ij 的可信度为
将各指标的可信度进行归一化处理,得到指标权重,表示为
3 模型验证
3.1 构建岩爆倾向等级评价指标体系
内因和外因共同作用造成岩爆灾害。其中,内因是指岩石的物理力学性质,外因是指不同的开挖方式造成岩石应力重分布。目前,关于岩爆等级评价指标的选择国内外没有统一的理论标准,因此本文在参考大量岩爆实例和已有岩爆判据(Wang et al.,2020;Russenes,1974)的基础上,选择4个评价指标:(1)应力系数,即岩石最大切向应力与单轴抗压强度之比σθ/σc,比值大小代表岩体的强度条件,而岩爆灾害多发生于硬岩条件下;(2)脆性系数,即岩石单轴抗压强度与抗拉强度之比σc/σt,比值大小代表岩体的稳定性条件,稳定性越差则岩爆程度越严重;(3)弹性应变指数Wet,反映了岩石储存弹性变形势能的能力,数值越大则发生岩爆灾害的倾向性越大;(4)岩石完整性系数kv,岩石完整性越好,越容易储存应变能量,发生岩爆灾害的倾向性也越大。
以上评价指标相互独立又相互关联,参照王元汉等(1998)、殷欣等(2020)和刘晓悦等(2021)的研究,将岩爆等级划分为4个区间,即:无岩爆(Ⅰ级)、弱岩爆(Ⅱ级)、中等岩爆(Ⅲ级)和剧烈岩爆(Ⅳ级),并确定了岩爆等级与评价指标之间的关系,如
在工程实践中常采用微震监测法和声发射监测法进行岩爆监测,岩爆等级的特征描述如表1所示。
表1 岩爆倾向等级的特征描述
Table 1
岩爆倾向等级 | 分级特征描述 |
---|---|
无岩爆(Ⅰ级) | 表层围岩没有出现爆炸脱落、剥离等现象,对施工基本无影响 |
弱岩爆(Ⅱ级) | 表层围岩出现爆裂脱落、剥离等现象,并且内部不时出现轻微噼啪、撕裂声,对施工影响较小 |
中等岩爆(Ⅲ级) | 围岩爆裂脱落、剥离现象较严重,出现少量弹射、偶尔抛射现象,时常发出尖锐爆裂声音,破坏范围明显,对施工有一定影响 |
强岩爆(Ⅳ级) | 围岩大片爆裂脱落,出现强烈弹射现象,持续发出尖锐爆裂声音,围岩急剧变形,严重影响施工 |
3.2 深埋隧洞工程数据
为了验证本文所用模型的准确性,选取国内外13个典型深埋隧洞原始数据进行数据处理,指标实测值见表2。
表2 典型深埋隧洞原始数据
Table 2
深埋隧洞编号 | 工程名称 | 评价指标实测值 | 文献来源 | |||
---|---|---|---|---|---|---|
σθ/σc | σc/σt | Wet | kv | |||
1 | 天生桥二级水电站引水隧洞 | 0.30 | 24.0 | 6.60 | 0.73 | 陈祥等(2009)、张彪等(2017) |
2 | 瀑布沟地下水电站地下洞室 | 0.33 | 24.6 | 5.00 | 0.80 | 徐琛等(2017) |
3 | 二滩水电站2#支洞 | 0.41 | 29.7 | 7.30 | 0.64 | 殷欣等(2020)、张彪等(2017) |
4 | 鲁布革水电站引水隧洞 | 0.23 | 27.8 | 7.80 | 0.69 | 陈祥等(2009)、张彪等(2017) |
5 | 李家峡水电站地下洞室 | 0.10 | 23.0 | 5.70 | 0.34 | 张研等(2011)、殷欣等(2020) |
6 | 鱼子溪水电站引水隧洞 | 0.53 | 14.8 | 9.00 | 0.71 | 陈祥等(2009)、张彪等(2017) |
7 | Vietas水电站引水隧洞 | 0.44 | 26.7 | 5.50 | 0.73 | 殷欣等(2020) |
8 | 太平驿水电站引水隧洞 | 0.38 | 17.6 | 9.00 | 0.75 | 陈祥等(2009)、张研等(2011) |
9 | 锦屏一级水电站引水隧洞 | 0.67 | 12.0 | 3.50 | 0.53 | 徐琛等(2017) |
10 | 猴子岩水电站引水隧洞 | 0.54 | 28.5 | 4.71 | 0.55 | 贾哲强等(2016) |
11 | 锦屏二级水电站引水隧洞 | 0.82 | 18.5 | 3.80 | 0.80 | 陈祥等(2009)、张彪等(2017) |
12 | 龙羊峡水电站地下洞室 | 0.11 | 31.2 | 7.40 | 0.42 | 刘冉等(2019)、殷欣等(2020) |
13 | 挪威Sima水电站地下厂房 | 0.27 | 21.7 | 5.00 | 0.81 | 李鹏程等(2019) |
3.3 指标权重的确定
采用正态隶属度函数计算各指标的基本概率分配值,利用Jousselme距离计算指标间距离,然后通过指标间距离值确定指标的信任度,最后通过指标信任度计算每个指标的权重,具体计算过程如下:
表3 指标概率分配表
Table 3
评价指标 | 岩爆倾向等级基本概率分配 | |||
---|---|---|---|---|
Ⅰ级 | Ⅱ级 | Ⅲ级 | Ⅳ级 | |
σθ /σc | 0.5000 | 0.5000 | 0.0000 | 0.0000 |
σc /σt | 0.0000 | 0.1930 | 0.8070 | 0.0000 |
Wet | 0.0000 | 0.0000 | 0.0000 | 1.0000 |
kv | 0.0000 | 0.0000 | 0.7800 | 0.2200 |
(2)根据
表4 各指标间距离
Table 4
评价指标 | σθ /σc | σc /σt | Wet | kv |
---|---|---|---|---|
σθ /σc | 0.0000 | 0.7055 | 0.8660 | 0.7605 |
σc /σt | 0.7055 | 0.0000 | 0.9188 | 0.2078 |
Wet | 0.8660 | 0.9188 | 0.0000 | 0.7800 |
kv | 0.7605 | 0.2078 | 0.7800 | 0.0000 |
(3)根据
表5 深埋隧洞指标权重值
Table 5
深埋隧洞编号 | 工程名称 | 指标权重值 | |||
---|---|---|---|---|---|
σθ /σc | σc /σt | Wet | kv | ||
1 | 天生桥二级水电站引水隧洞 | 0.1896 | 0.3315 | 0.1235 | 0.3554 |
2 | 瀑布沟地下水电站地下洞室 | 0.1693 | 0.2534 | 0.3206 | 0.2567 |
3 | 二滩水电站2#支洞 | 0.2940 | 0.3451 | 0.0409 | 0.3200 |
4 | 鲁布革水电站引水隧洞 | 0.1581 | 0.4469 | 0.1048 | 0.2902 |
5 | 李家峡水电站地下洞室 | 0.4401 | 0.0715 | 0.0310 | 0.4574 |
6 | 鱼子溪水电站引水隧洞 | 0.2878 | 0.3265 | 0.1206 | 0.2651 |
7 | Vietas水电站引水隧洞 | 0.2277 | 0.3522 | 0.1494 | 0.2707 |
8 | 太平驿水电站引水隧洞 | 0.0524 | 0.3439 | 0.2620 | 0.3417 |
9 | 锦屏一级水电站引水隧洞 | 0.2811 | 0.2335 | 0.2837 | 0.2017 |
10 | 猴子岩水电站引水隧洞 | 0.2877 | 0.2622 | 0.2501 | 0.2000 |
11 | 锦屏二级水电站引水隧洞 | 0.2107 | 0.2791 | 0.2179 | 0.2923 |
12 | 龙羊峡水电站地下洞室 | 0.3664 | 0.1153 | 0.0456 | 0.4726 |
13 | 挪威Sima水电站地下厂房 | 0.1487 | 0.2104 | 0.3721 | 0.2688 |
3.4 预测结果及分析
采用置信度识别准则对13个典型深埋隧洞岩爆工程进行等级识别,取置信度λ=0.5。均化系数ε的目的是各指标的综合属性测度在单指标属性测度区间的上限值和下限值取一个均衡的结果,因此ε取值不同,各评价指标属性测度也发生不同程度的变化,对预测准确度会产生一定的影响。为了获得最优ε值,通过收集实际深埋隧洞岩爆工程指标数据进行计算,使ε步长为0.1并在区间[0.05,0.95]选值,最终选择岩爆预测最正确的ε值。由图1可知,样本数据为13个,当ε取值为0.05、0.45时,正确预测岩爆数量为9个,此时预测准确度最低;当ε取值为0.65时,正确预测岩爆数量为12个,此时预测准确度最好,因此本文选择最优ε值为0.65。
图1
图1
正确预测岩爆个数与ε值的关系
Fig.1
Relationship between the number of rockburst prediction correctly and ε value
表6 岩爆预测评价等级
Table 6
深埋隧洞编号 | 综合属性测度 | 本文预测结果 | 实际级别 | |||
---|---|---|---|---|---|---|
1 | 0.0664 | 0.3162 | 0.2836 | 0.3338 | Ⅲ | Ⅲ |
2 | 0.0475 | 0..2249 | 0.3162 | 0.4113 | Ⅲ | Ⅲ |
3 | 0.1091 | 0.4647 | 0.3853 | 0.0409 | Ⅱ | Ⅱ |
4 | 0.1145 | 0.4075 | 0.2977 | 0.1803 | Ⅱ | Ⅰ或Ⅱ |
5 | 0.5742 | 0.2366 | 0.1745 | 0.0147 | Ⅰ | Ⅰ |
6 | 0.0023 | 0.2390 | 0.5067 | 0.2520 | Ⅲ | Ⅲ |
7 | 0.2528 | 0.2473 | 0.2148 | 0.2551 | Ⅱ | Ⅱ |
8 | 0.0165 | 0.0884 | 0.3325 | 0.5626 | Ⅳ(误判) | 400多例,等级不定 |
9 | 0.0901 | 0.3367 | 0.4654 | 0.1078 | Ⅲ | Ⅲ |
10 | 0.0935 | 0.3870 | 0.3509 | 0.1686 | Ⅲ | Ⅱ或Ⅲ |
11 | 0.0000 | 0.1209 | 0.3937 | 0.4854 | Ⅲ | Ⅲ或Ⅳ |
12 | 0.5478 | 0.2745 | 0.1372 | 0.0404 | Ⅰ | Ⅰ |
13 | 0.0286 | 0.1573 | 0.3447 | 0.4693 | Ⅳ | Ⅲ或Ⅳ |
通过模型验算,说明本文提出的正态隶属度—属性区间识别模型对于岩爆倾向等级的预测有一定的实践意义。该模型有以下2个优势:
(1)岩爆的发生与工程地质有很大的关系,在实践中,地质参数通常为一个区间值,与传统模型相比,属性区间识别模型能够更好地适用于各等级指标值为区间的问题,且该模型的准确率有所提升。
(2)指标权重是判定分类结果的重要因素,基于Jousselme距离和指标信任度确定的权重方法考虑了各指标之间、指标与标准等级之间的交互关系的计算方法,因此对于同一指标体系下不同的指标实测值有不同的权重,提高了模型的稳健性。
4 实例验证
双江口水电站位于四川省阿坝藏族羌族自治州马尔康县与金川县交界处,坝高312 m(李航,2020)。地层岩性主要为可尔因花岗岩杂岩体,由钾长石、斜长石、石英、黑云母和白云母等矿物组成,具细—中粒结构,片麻状构造,岩体坚硬完整,含水量少,局部存在断层和挤压破碎带。其中,地下厂房(SPD9)处硐口高程为2 267.6 m,主厂房水平埋深400~640 m,垂直埋深320~500 m,在硐深400 m处测出最大地应力为37.82 MPa,属于高地应力区。SPD9平硐处岩体坚硬完整,裂隙发育稀疏,极易发生岩爆等围岩破坏,因此需要提前进行岩爆倾向等级预判。
表7 SPD9洞室指标原始数据及指标权重值
Table 7
桩号 | 评价指标实测值 | 指标权重值 | ||||||
---|---|---|---|---|---|---|---|---|
σθ /σc | σc /σt | Wet | kv | σθ /σc | σc /σt | Wet | kv | |
0+123 | 0.50 | 10.2 | 5.0 | 0.70 | 0.2510 | 0.1928 | 0.3293 | 0.2269 |
0+200 | 0.62 | 9.4 | 5.5 | 0.75 | 0.1261 | 0.2878 | 0.2931 | 0.2931 |
0+301 | 0.57 | 10.24 | 5.0 | 0.70 | 0.2889 | 0.1542 | 0.2889 | 0.2681 |
0+406 | 1.09 | 10.56 | 5.5 | 0.75 | 0.2582 | 0.2748 | 0.2748 | 0.1923 |
0+570 | 0.56 | 10.61 | 4.5 | 0.70 | 0.3046 | 0.0817 | 0.3134 | 0.3003 |
表8 SPD9洞室岩爆预测评价等级
Table 8
桩号 | 综合属性测度 | 本文预测结果 | 实际级别 | |||
---|---|---|---|---|---|---|
μ1 | μ2 | μ3 | μ4 | |||
0+123 | 0.0000 | 0.1276 | 0.4800 | 0.3924 | Ⅲ | Ⅲ |
0+200 | 0.0000 | 0.0177 | 0.5419 | 0.4404 | Ⅲ | Ⅲ |
0+301 | 0.0000 | 0.1126 | 0.4633 | 0.4241 | Ⅲ | Ⅲ |
0+406 | 0.0000 | 0.0000 | 0.2615 | 0.7385 | Ⅳ | Ⅳ |
0+570 | 0.0000 | 0.1480 | 0.4524 | 0.3996 | Ⅲ | Ⅲ |
5 结论
(1)针对深埋隧洞开挖过程中遇到的岩爆问题,选取应力系数、脆性系数、弹性应变指数和岩石完整性系数作为岩爆倾向等级预测指标,建立了岩爆等级与评价指标之间的对应关系。
(2)采用正态隶属度函数确定指标概率密度,并利用Jousselme距离计算指标间相互距离,计算得到4个不同指标的权重值。与传统赋权方法相比,该方法充分考虑了指标间、指标与标准等级间的联系,对同一指标体系不同工程实例赋予不同权重值,使权重的计算更具科学性。
(3)建立了正态隶属度—属性区间识别模型,将该模型应用到13个典型深埋隧洞工程进行模型准确性测试,并以双江口水电站SPD9厂房为例进行工程实例验证,验证结果精确度较高,表明该模型在实际工程中有一定的应用价值。
http://www.goldsci.ac.cn/article/2022/1005-2518/1005-2518-2022-30-3-404.shtml
参考文献
Knowledge-based and data-driven fuzzy modeling for rockburst prediction
[J].
Judgment indexes and classification criteria of rock-burst with the extension judgment method
[J].
Prediction of rockburst classification using random forest
[J].
Integrated prediction for rockburst of underground powerhouse of Houziyan hydropower station
[J].
Bursting liability indices of coal
[J].
Research on Tunnel Rockburst Prediction Method Based on Combination Weight Ideal Point Method-Database
[D].
Grey evaluation rockburst prediction model based on normal whitening weight function
[J].
A PSO-RBF neural network model for rockburst tendency prediction
[J].
Study on construction risk of shield tunnel under-passing existing railway by fuzzy comprehensive evaluation
[J].
Grading prediction model of rockburst based on rough set-multidiensional normal cloud
[J].
Rockburst prediction of multi-dimensional cloud model based on improved hierarchical analytic method and critic method
[J].
Rockburst prediction based on a modified grey
[J].
Analysis of Rock Spalling for Tunnels in Steep Valley Sides
[D].
Equipment management performan-ce assessment based on cloud model and evidence theory
[J].
A prediction model of rock burst in tunnel based on the improved MATLAB-BP neural network
[J].
Prediction of intensity classification of rockburst based on deep neural network
[J].
Prediction of rockburst based on multidimensional connection cloud model and set pair analysis
[J].
Model for classification and prediction of rock burst intensity in a deep underground engineering with rough set and efficacy coefficient method
[J].
A normal cloud model-based study of grading prediction of rockburst intensity in deep underground engineering
[J].
Method of fuzzy comprehensive evaluations for rockburst prediction
[J].
Prediction and classification of strain mode rockburst based on five-factor crierion and combined weight-ideal point method
[J].
Prediction model of rockburst intensity classification based on combined weighting and attribute interval recognition theory
[J].
rockburst intensity grade based on index distance and uncertainty measure
[J].
Rockburst criterion based on artificial neural networks and nonlinear regression
[J].
Rock burst and its criteria and control
[J].
Application of extension evaluation based on rough set in rockburst prediction
[J] .
Method of identifying rockburst grades based on Gaussian process machine learning
[J].
A comprehensive assessment model for severity degree of dam failure impact based on attribute interval recognition theory
[J].
岩爆的判别指标和分级标准及可拓综合判别方法
[J].
猴子岩水电站地下厂房岩爆综合预测研究
[J].
基于组合权重理想点法—数据库的隧洞岩爆预测方法研究
[D].
基于正态白化权函数的灰评估岩爆预测模型
[J].
基于PSO-RBF神经网络模型的岩爆倾向性预测
[J].
盾构下穿运营铁路施工风险模糊综合评价方法研究
[J].
岩爆分级预测的粗糙集—多维正态云模型
[J].
基于改进层次法与CRITIC法的多维云模型岩爆预测
[J].
基于改进的灰评估模型在岩爆中的预测研究
[J].
基于改进MATLAB-BP神经网络算法的隧道岩爆预测模型
[J].
基于深度神经网络的岩爆烈度分级预测研究
[J].
深埋地下工程岩爆烈度分级预测的RS-功效系数模型
[J].
基于正态云模型的深埋地下工程岩爆烈度分级预测研究
[J].
岩爆预测的模糊数学综合评判方法
[J].
基于组合权重—理想点法的应变型岩爆五因素预测分级
[J].
基于组合赋权和属性区间识别理论的岩爆烈度分级预测模型
[J].
基于指标距离与不确定度量的岩爆云模型预测研究
[J].
基于人工神经网络及非线性回归的岩爆判据
[J].
岩爆及其判据和防治
[J].
基于粗糙集的可拓评判在岩爆预测中的应用
[J].
基于高斯过程机器学习的岩爆等级识别方法
[J].
属性区间识别模型在溃坝后果综合评价中的应用
[J].
/
〈 |
|
〉 |
