img

QQ群聊

img

官方微信

  • CN 62-1112/TF 
  • ISSN 1005-2518 
  • 创刊于1988年
高级检索

黄金科学技术, 2022, 30(3): 404-413 doi: 10.11872/j.issn.1005-2518.2022.03.037

采选技术与矿山管理

基于正态隶属度—属性区间识别模型的岩爆倾向等级预测

贡力,1,2, 陆丽丽,1, 靳春玲1, 梁栋3, 周汉国3, 谢平3

1.兰州交通大学土木工程学院,甘肃 兰州 730070

2.兰州交通大学调水工程及输水安全研究所,甘肃 兰州 730070

3.中国中铁科学研究院有限公司,四川 成都 610000

Propensity Grade Prediction of Rockburst Based on Normal Membership-Attribute Interval Recognition Model

GONG Li,1,2, LU Lili,1, JIN Chunling1, LIANG Dong3, ZHOU Hanguo3, XIE Ping3

1.College of Civil Engineering, Lanzhou Jiaotong University, Lanzhou 730070, Gansu, China

2.Institute of Water Diversion Engineering and Water Security, Lanzhou Jiaotong University, Lanzhou 730070, Gansu, China

3.China Railway Academy Co. , Ltd. , Chengdu 610000, Sichuan, China

通讯作者: 陆丽丽(1994-),女,甘肃白银人,硕士研究生,从事土木工程建造与管理方面的研究工作。lullin@qq.com

收稿日期: 2021-03-05   修回日期: 2021-08-29  

基金资助: 国家自然科学基金项目“西北盐渍干寒地区多因素耦合作用下长距离明渠的劣化机理研究”.  51969011
甘肃省科技计划资助项目“ 黄河流域突发水污染应急调控关键技术与应用”.  20JR10RA274

Received: 2021-03-05   Revised: 2021-08-29  

作者简介 About authors

贡力(1977-),男,江苏丹阳人,教授,博士生导师,从事深埋隧道及输水工程方面的研究工作gongl@mail.lzjtu.cn , E-mail:gongl@mail.lzjtu.cn

摘要

岩爆是地下工程开挖面临的关键问题之一,为了准确预测深埋隧洞中岩爆烈度倾向等级,提出了正态隶属度—属性区间识别模型的岩爆预测方法。针对岩爆倾向等级属于典型的多属性有序分割类问题,构建了属性区间识别模型,并将岩爆倾向等级划分为4个等级进行预测。根据岩爆发生的成因和机理,选取应力系数、脆性系数、弹性应变指数和岩石完整性系数作为预测指标,考虑各指标之间、指标与标准等级之间的交互关系,采用正态隶属度函数和Jousselme距离计算评价指标权重。结合13个深埋隧洞工程对该预测模型进行准确性测试,并以双江口水电站SPD9厂房为例进行工程实例验证,该模型预测结果与实际相吻合,证明该模型用于具体工程实践中是可行且有效的,研究结果可为类似深埋隧洞岩爆倾向等级预测提供新的思路。

关键词: 深埋隧洞 ; 岩爆倾向等级 ; 岩爆预测 ; 正态隶属度函数 ; 属性区间识别模型 ; 均化系数

Abstract

Rockburst is one of the key problems in underground engineering excavation,in order to accurately predict the grade of rockburst propensity in deep-buried tunnels,a rockburst prediction method based on the normal membership degree-attribute interval recognition model was proposed.Due to the rockburst propensity is a typical multi-attribute orderly segmentation problem,the attribute interval recognition model was constructed to divide the rockburst propensity into four grades for prediction.The occurrence of rockburst is affected by engineering geology,and the geological parameter is usually an interval value.The attribute interval recognition model can be better applied to the problem that each grade index value is an interval value.According to the mechanism of rockburst,the stress coefficient,brittleness coefficient,elastic strain index and rock integrity coefficient were selected as predictive indicators,considering the interaction between the indicators and the standard grade,the normal membership function and Jousselme distance were used to calculate the weight of the index.The method had different weights for the measured values of different indexes under the same index system,which can avoid the drawback that the traditional method didn’t consider the interaction relationship between the indexes,and improve the robustness of the model.The accuracy of the prediction model was tested with 13 deep-buried tunnel projects.Since the value of the averaging coefficient had a greater impact on the prediction accuracy of the model,in order to obtain the optimal value of the averaging coefficient,the step size was 0.1 and selected within the interval [0.05,0.95].The analysis shows that when the averaging coefficient is 0.65,the prediction accuracy of the model is the highest,which is 92.31%.The left bank of SPD9 in the Huangjiangkou hydropower project was used to verify the model.The prediction results are consistent with the actual rockburst propensity grade.The prediction results prove that the model is feasible and effective in specific engineering practice,and it can provide new ideas for predicting the rockburst propensity grade in similar deep-buried tunnels.

Keywords: deep-buried tunnels ; rockburst propensity grade ; rockburst prediction ; normal membership function ; attribute interval recognition model ; averaging coefficient

PDF (1004KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

贡力, 陆丽丽, 靳春玲, 梁栋, 周汉国, 谢平. 基于正态隶属度—属性区间识别模型的岩爆倾向等级预测[J]. 黄金科学技术, 2022, 30(3): 404-413 doi:10.11872/j.issn.1005-2518.2022.03.037

GONG Li, LU Lili, JIN Chunling, LIANG Dong, ZHOU Hanguo, XIE Ping. Propensity Grade Prediction of Rockburst Based on Normal Membership-Attribute Interval Recognition Model[J]. Gold Science and Technology, 2022, 30(3): 404-413 doi:10.11872/j.issn.1005-2518.2022.03.037

岩爆是一种发生在深埋隧洞中的常见地质灾害,是隧洞开挖工程中由于岩体本身积聚的应变能量超过岩体破坏时消耗的能量,岩石受到应力重分布作用将多余能量释放出来致使岩石发生爆裂、剥离的现象(李任豪等,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为全体待评价隧洞,DiDi=1,2,…,m)为若干典型的待评价隧洞,每个隧洞有n个评价指标I1I2,…,Inxij 为第i个标段的第j个评价指标的测量值。假定GD上的某类属性空间,评价集(E1E2Ek )为属性测度区间空间G的有序分割类,根据已知的分级标准构造标准的分级矩阵为

E1E2EkI1a11,b11a12,b12a1k,b1kI2a21,b21a22,b22a2k,b2kInan1,bn1an2,bn2ank,bnk

式中:ajp,bjp1jn,1pk)为第j个评价指标在评价集Ep 上的端点值,且满足aj1<aj2<<ajk,bj1<bj2<<bjk 或者aj1>aj2>>ajk,bj1>bj2>>bjk

1.1 单指标属性测度区间的计算

xij 具有属性Ep (1≤p≤k)的属性测度区间记为μijp=μ̲ijp,μ¯ijp,则有:

xijaj1(xijaj1时,有μ̲ij1=1,μ̲ij2==μ̲ijk=0;
xijbj1(xijbj1)时,有μ¯ij1=1,μ¯ij2==μ¯ijk=0;
xijajk(xijajk时,有μ̲ijk=1,μ̲ij1==μ̲ij(k-1)=0;
xijbjk(xijbjk时,有μ¯ijk=1,μ¯ij1==μ¯ij(k-1)=0;

ajyxijaj(y+1)aj(y+1)xijajy)(1yk-1)时,则有:

μ̲ijy=xij-aj(y+1)ajy-aj(y+1)              μ̲ijp=0(p<yp>y+1)μ̲ij(y+1)=xij-ajyajy-aj(y+1)         

bjyxijbjy+1bj(y+1)xijbjy)(1yk-1)时,则有:

μ¯ijy=xij-bj(y+1)bjy-bj(y+1)               μ¯ijp=0(p<yp>y+1)μ¯ij(y+1)=xij-bjybjy-bj(y+1)          

1.2 综合属性测度的计算

将均化系数(邹强等,2011)引入属性测度区间,计算各指标的综合属性测度μijp,表示为

μijp=εμ̲ijp+1-εμ¯ijp

式中:ε为均化系数,ε∈(0,1)。

1.3 待评价隧洞等级识别

通过综合属性测度计算第i个评价标段属于等级Ep 的综合属性测度μip,表示为

μip=j=1nwjμijp(1im,1pk)

式中:wj 为第j个评价指标的权重。

按照置信度准则,取置信度λ(0.5≤λ≤1),若p满足式(10),则可确定评价对象xi 属于第p个评价等级Ep

pi=minp:y=1pμiyλ,1pk

2 权重的确定

2.1 正态隶属度函数

评价对象的性质决定评价指标,模糊理论的关键是建立指标隶属度函数(林大涌等,2018)。正态隶属度函数具有模糊性,且能够减少主观因素的影响,因此本文采用正态隶属度函数来确定岩爆评价指标属于某评价等级的隶属度,具体公式如下:

u(x)=e-(x-x¯)/c2
c=ak+-ak-2-ln (0.5)

式中:x为指标的实测值;x¯为某等级区间的均值;ak+ak-为评价等级的端点值。

各评价等级指标边界值与相邻2个评价等级隶属度相同,因此计算某评价等级端点值时,式(11)中的c应为2c。假设各评价等级kak-,a¯k,ak+表示,其中,a¯k表示该等级区间的均值,ak-ak+表示该等级区间的上限值和下限值,x为岩爆某指标的实测值。用任意区间k对应岩爆评价等级中的非边界评价等级,则以k=1对应第一个评级等级,k+1对应最后一个评价等级,将正态隶属度函数推广到一般形式,可得到各区间的一般表达式(以效益型指标为例):

uk+1x=1  (x<a¯k+1) e--ln 0.5(x-a¯k+1ak+1+-a¯k+1)2  (a¯k+1xak+1+)1-e--ln 0.5(x-a¯ka¯k-ak+1+)2  (a k+1+xa¯k  )
ukx=0 (xa¯k+1,xa¯k-1)1-e--ln 0.5(x-a¯k+1a¯k-a¯k+1)2 (a¯k+1<xak-)e-2-ln 0.5(x-a¯kak+-ak-)2 (ak-x<ak+)1-e--ln 0.5(x-a¯k-1a¯k-1-ak+)2 (ak-x<ak+)
u1x=0 (xa¯k)1-e--ln 0.5(x-a¯ka1--a¯k)2 (a¯1<xa1-)e--ln 0.5(x-a¯1a1--a¯1)2 (a1-<xa¯1)1 (x>a¯1)

根据式(13)~式(15)确定各评价指标对应各评价等级的隶属度函数。

2.2 指标权重的确定

评价指标的权重值是度量指标重要性的方法。考虑各指标间以及各等级间的相互影响,采用指标间的相互距离计算指标权重。正态隶属度函数确定的指标隶属度又称为指标概率密度,由概率密度计算指标距离进而确定权重,计算步骤如下:

IiIj 为指标集合中的2个评价指标,则其对应的各等级概率密度为mimj,计算指标IiIj 之间的Jousselme距离dijShen et al.,2019),表示为

dij=0.5mi2+mj2-2mi,mj

式中:mimj;mi,mj为2个向量内积。

由指标间的距离可以计算指标IiIj 的信任度,表示为

αij=1-dij(i,j=1,2,,n)

指标Ii 相对于指标Ij 的可信度为

sup Mi=j=1,jinαij(i=1,2,,n)

将各指标的可信度进行归一化处理,得到指标权重,表示为

wi=    sup Mii=1nsup Mi(i=1,2,,n)

3 模型验证

3.1 构建岩爆倾向等级评价指标体系

内因和外因共同作用造成岩爆灾害。其中,内因是指岩石的物理力学性质,外因是指不同的开挖方式造成岩石应力重分布。目前,关于岩爆等级评价指标的选择国内外没有统一的理论标准,因此本文在参考大量岩爆实例和已有岩爆判据(Wang et al.,2020Russenes,1974)的基础上,选择4个评价指标:(1)应力系数,即岩石最大切向应力与单轴抗压强度之比σθ/σc,比值大小代表岩体的强度条件,而岩爆灾害多发生于硬岩条件下;(2)脆性系数,即岩石单轴抗压强度与抗拉强度之比σc/σt,比值大小代表岩体的稳定性条件,稳定性越差则岩爆程度越严重;(3)弹性应变指数Wet,反映了岩石储存弹性变形势能的能力,数值越大则发生岩爆灾害的倾向性越大;(4)岩石完整性系数kv,岩石完整性越好,越容易储存应变能量,发生岩爆灾害的倾向性也越大。

以上评价指标相互独立又相互关联,参照王元汉等(1998)殷欣等(2020)刘晓悦等(2021)的研究,将岩爆等级划分为4个区间,即:无岩爆(Ⅰ级)、弱岩爆(Ⅱ级)、中等岩爆(Ⅲ级)和剧烈岩爆(Ⅳ级),并确定了岩爆等级与评价指标之间的关系,如式(20)所示。由式(20)可以看出,每个评价指标集的端点值均满足aj1<aj2<<ajk,bj1<bj2<<bjk,因此可以判定岩爆倾向等级预测属于典型的有序分割类问题。

σθ/σθ0.1,0.30.3,0.50.5,0.70.7,0.9σo/σt40.0,46.026.7,40.014.5,26.70.0,14.5Wet0.0,2.02.0,3.53.5,5.05.0,6.5kv0.0,0.550.55,0.650.65,0.750.75,1.0

在工程实践中常采用微震监测法和声发射监测法进行岩爆监测,岩爆等级的特征描述如表1所示。

表1   岩爆倾向等级的特征描述

Table 1  Characteristic description of rockburst tendency grade

岩爆倾向等级分级特征描述
无岩爆(Ⅰ级)表层围岩没有出现爆炸脱落、剥离等现象,对施工基本无影响
弱岩爆(Ⅱ级)表层围岩出现爆裂脱落、剥离等现象,并且内部不时出现轻微噼啪、撕裂声,对施工影响较小
中等岩爆(Ⅲ级)围岩爆裂脱落、剥离现象较严重,出现少量弹射、偶尔抛射现象,时常发出尖锐爆裂声音,破坏范围明显,对施工有一定影响
强岩爆(Ⅳ级)围岩大片爆裂脱落,出现强烈弹射现象,持续发出尖锐爆裂声音,围岩急剧变形,严重影响施工

新窗口打开| 下载CSV


3.2 深埋隧洞工程数据

为了验证本文所用模型的准确性,选取国内外13个典型深埋隧洞原始数据进行数据处理,指标实测值见表2

表2   典型深埋隧洞原始数据

Table 2  Original index of typical deep-buried tunnels

深埋隧洞编号工程名称评价指标实测值文献来源
σθcσctWetkv
1天生桥二级水电站引水隧洞0.3024.06.600.73陈祥等(2009)、张彪等(2017)
2瀑布沟地下水电站地下洞室0.3324.65.000.80徐琛等(2017)
3二滩水电站2#支洞0.4129.77.300.64殷欣等(2020)、张彪等(2017)
4鲁布革水电站引水隧洞0.2327.87.800.69陈祥等(2009)、张彪等(2017)
5李家峡水电站地下洞室0.1023.05.700.34张研等(2011)、殷欣等(2020)
6鱼子溪水电站引水隧洞0.5314.89.000.71陈祥等(2009)、张彪等(2017)
7Vietas水电站引水隧洞0.4426.75.500.73殷欣等(2020)
8太平驿水电站引水隧洞0.3817.69.000.75陈祥等(2009)、张研等(2011)
9锦屏一级水电站引水隧洞0.6712.03.500.53徐琛等(2017)
10猴子岩水电站引水隧洞0.5428.54.710.55贾哲强等(2016)
11锦屏二级水电站引水隧洞0.8218.53.800.80陈祥等(2009)、张彪等(2017)
12龙羊峡水电站地下洞室0.1131.27.400.42刘冉等(2019)、殷欣等(2020)
13挪威Sima水电站地下厂房0.2721.75.000.81李鹏程等(2019)

新窗口打开| 下载CSV


3.3 指标权重的确定

采用正态隶属度函数计算各指标的基本概率分配值,利用Jousselme距离计算指标间距离,然后通过指标间距离值确定指标的信任度,最后通过指标信任度计算每个指标的权重,具体计算过程如下:

(1)根据表2式(20),利用式(13)~式(15)的指标隶属度函数计算各指标基本概率分配值。以天生桥二级水电站引水隧洞为例,计算结果列于表3中。

表3   指标概率分配表

Table 3  Probability assignment of each index

评价指标岩爆倾向等级基本概率分配
Ⅰ级Ⅱ级Ⅲ级Ⅳ级
σθ /σc0.50000.50000.00000.0000
σc /σt0.00000.19300.80700.0000
Wet0.00000.00000.00001.0000
kv0.00000.00000.78000.2200

新窗口打开| 下载CSV


(2)根据式(16)计算指标间距离,计算结果见表4

表4   各指标间距离

Table 4  Distance between each index

评价指标σθ /σcσc /σtWetkv
σθ /σc0.00000.70550.86600.7605
σc /σt0.70550.00000.91880.2078
Wet0.86600.91880.00000.7800
kv0.76050.20780.78000.0000

新窗口打开| 下载CSV


(3)根据式(17)和式(18)计算指标间信任度。指标间距离越小,指标相互支持程度越大,可信度越高,13个深埋隧洞指标权重见表5

表5   深埋隧洞指标权重值

Table 5  Index weight of deep buried tunnels

深埋隧洞编号工程名称指标权重值
σθ /σcσc /σtWetkv
1天生桥二级水电站引水隧洞0.18960.33150.12350.3554
2瀑布沟地下水电站地下洞室0.16930.25340.32060.2567
3二滩水电站2#支洞0.29400.34510.04090.3200
4鲁布革水电站引水隧洞0.15810.44690.10480.2902
5李家峡水电站地下洞室0.44010.07150.03100.4574
6鱼子溪水电站引水隧洞0.28780.32650.12060.2651
7Vietas水电站引水隧洞0.22770.35220.14940.2707
8太平驿水电站引水隧洞0.05240.34390.26200.3417
9锦屏一级水电站引水隧洞0.28110.23350.28370.2017
10猴子岩水电站引水隧洞0.28770.26220.25010.2000
11锦屏二级水电站引水隧洞0.21070.27910.21790.2923
12龙羊峡水电站地下洞室0.36640.11530.04560.4726
13挪威Sima水电站地下厂房0.14870.21040.37210.2688

新窗口打开| 下载CSV


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


ε值为0.65时,各深埋隧洞的综合属性测度结果见表6。从表6中可以看出,除隧洞编号为8的工程之外,其他深埋隧洞岩爆等级的预测结果与实际情况一致,准确率达到92.31%。对于8号太平驿水电站引水隧洞,根据陈祥等(2009)研究记录,该隧洞共发生400多例岩爆,岩爆等级不定,本文将其判定为强岩爆(Ⅳ),可以采取应力释放钻孔、超前导洞等措施进行超前应力解除,降低岩爆发生的概率。因此,从工程安全角度考虑,将其判定为高等级岩爆是允许的。

表6   岩爆预测评价等级

Table 6  Rockburst prediction evaluation level

深埋隧洞编号综合属性测度本文预测结果实际级别
μ1μ2μ3μ4
10.06640.31620.28360.3338
20.04750..22490.31620.4113
30.10910.46470.38530.0409
40.11450.40750.29770.1803Ⅰ或Ⅱ
50.57420.23660.17450.0147
60.00230.23900.50670.2520
70.25280.24730.21480.2551
80.01650.08840.33250.5626Ⅳ(误判)400多例,等级不定
90.09010.33670.46540.1078
100.09350.38700.35090.1686Ⅱ或Ⅲ
110.00000.12090.39370.4854Ⅲ或Ⅳ
120.54780.27450.13720.0404
130.02860.15730.34470.4693Ⅲ或Ⅳ

注:μ1~μ4表示各深埋隧洞属于Ⅰ~Ⅳ等级的综合属性测度

新窗口打开| 下载CSV


通过模型验算,说明本文提出的正态隶属度—属性区间识别模型对于岩爆倾向等级的预测有一定的实践意义。该模型有以下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平硐处岩体坚硬完整,裂隙发育稀疏,极易发生岩爆等围岩破坏,因此需要提前进行岩爆倾向等级预判。

通过本文提出的正态隶属度—属性区间识别模型,结合李航(2020)提供的SPD9厂房指标数据实测值对指标权重进行分析,结果见表7

表7   SPD9洞室指标原始数据及指标权重值

Table 7  Original index and index weight of SPD9 cavern

桩号评价指标实测值指标权重值
σθ /σcσc /σtWetkvσθ /σcσc /σtWetkv
0+1230.5010.25.00.700.25100.19280.32930.2269
0+2000.629.45.50.750.12610.28780.29310.2931
0+3010.5710.245.00.700.28890.15420.28890.2681
0+4061.0910.565.50.750.25820.27480.27480.1923
0+5700.5610.614.50.700.30460.08170.31340.3003

新窗口打开| 下载CSV


ε取值为0.65时,根据指标权重值和式(1)~式(9)计算指标的综合属性测度,同时,根据置信度识别准则,取置信度λ=0.5,得到每个待评价洞段的岩爆预测等级,结果见表8。将本文预测结果与李航(2020)研究结果进行比较,发现本文预测结果与现场实际相吻合,说明本文提出的正态隶属度—属性区间识别模型运用于具体工程实践中是可行且有效的。

表8   SPD9洞室岩爆预测评价等级

Table 8  Rockburst prediction evaluation level of SPD9 cavern

桩号综合属性测度本文预测结果实际级别
μ1μ2μ3μ4
0+1230.00000.12760.48000.3924
0+2000.00000.01770.54190.4404
0+3010.00000.11260.46330.4241
0+4060.00000.00000.26150.7385
0+5700.00000.14800.45240.3996

新窗口打开| 下载CSV


5 结论

(1)针对深埋隧洞开挖过程中遇到的岩爆问题,选取应力系数、脆性系数、弹性应变指数和岩石完整性系数作为岩爆倾向等级预测指标,建立了岩爆等级与评价指标之间的对应关系。

(2)采用正态隶属度函数确定指标概率密度,并利用Jousselme距离计算指标间相互距离,计算得到4个不同指标的权重值。与传统赋权方法相比,该方法充分考虑了指标间、指标与标准等级间的联系,对同一指标体系不同工程实例赋予不同权重值,使权重的计算更具科学性。

(3)建立了正态隶属度—属性区间识别模型,将该模型应用到13个典型深埋隧洞工程进行模型准确性测试,并以双江口水电站SPD9厂房为例进行工程实例验证,验证结果精确度较高,表明该模型在实际工程中有一定的应用价值。

中国矿业报

http://www.goldsci.ac.cn/article/2022/1005-2518/1005-2518-2022-30-3-404.shtml

参考文献

Adoko A CGokceoglu CWu Let al2013.

Knowledge-based and data-driven fuzzy modeling for rockburst prediction

[J].International Journal of Rock Mechanics and Mining Sciences,6186-95.

[本文引用: 1]

Chen XiangSun JinzhongZhang Jiekunet al2009.

Judgment indexes and classification criteria of rock-burst with the extension judgment method

[J].China Civil Engineering Journal,429):82-88.

Dong L JLi X BPeng K2013.

Prediction of rockburst classification using random forest

[J].Transactions of Nonferrous Metals Society of China,232):472-477.

[本文引用: 1]

Jia ZheqiangZhang RuZhang Yanfeiet al2016.

Integrated prediction for rockburst of underground powerhouse of Houziyan hydropower station

[J].Chinese Journal of Geotechnical Engineering,38Supp.2):110-116.

Kidybinski A1981.

Bursting liability indices of coal

[J].International Journal of Rock Mechanics and Mining Sciences,184):295-304.

[本文引用: 1]

Li Hang2020.

Research on Tunnel Rockburst Prediction Method Based on Combination Weight Ideal Point Method-Database

[D].BeijingChina University of Geosciences (Beijing).

Li PengchengYe YichengWang Qihuet al2019.

Grey evaluation rockburst prediction model based on normal whitening weight function

[J].Industrial Minerals & Processing,485):16-2226.

Li RenhaoGu HelongLi Xibinget al2020.

A PSO-RBF neural network model for rockburst tendency prediction

[J].Gold Science and Technology,281):134-141.

Lin DayongLei MingfengCao Haoronget al2018.

Study on construction risk of shield tunnel under-passing existing railway by fuzzy comprehensive evaluation

[J].Journal of Railway Science and Engineering,155):1347-1355.

Liu RanYe YichengZhang Guangquanet al2019.

Grading prediction model of rockburst based on rough set-multidiensional normal cloud

[J].Metal Mine,(3):48-55.

Liu XiaoyueYang WeiZhang Xuemei2021.

Rockburst prediction of multi-dimensional cloud model based on improved hierarchical analytic method and critic method

[J].Journal of Hunan University(Natural Sciences),482):118-124.

Pei QitaoLi HaiboLiu Yaqunet al2013.

Rockburst prediction based on a modified grey

[J].Chinese Journal of Rock Mechanics and Engineering,3210):2088-2093.

Russenes B F1974.

Analysis of Rock Spalling for Tunnels in Steep Valley Sides

[D].TrondheimNorwegian Institute of Technology.

[本文引用: 1]

Shen Y NZhang J B2019.

Equipment management performan-ce assessment based on cloud model and evidence theory

[J].Systems Engineering and Electronics,415):1049-1055.

[本文引用: 1]

Sun Chensheng2019.

A prediction model of rock burst in tunnel based on the improved MATLAB-BP neural network

[J].Journalof Chongqing Jiaotong University(Natural Science),3810):41-49.

Tian RuiMeng HaidongChen Shijianget al2020.

Prediction of intensity classification of rockburst based on deep neural network

[J].Journal of China Coal Society,45Supp.1):191-201.

Wang M WLiu Q YWang Xet al2020.

Prediction of rockburst based on multidimensional connection cloud model and set pair analysis

[J].International Journal of Geomechanics,201).DOI:10.1061/(ASCE)GM.1943-5622.0001546 .

[本文引用: 1]

Wang YingchaoJing HongwenJi Xianweiet al2014.

Model for classification and prediction of rock burst intensity in a deep underground engineering with rough set and efficacy coefficient method

[J].Journal of Central South University(Science and Technology),456):1992-1997.

Wang YingchaoJing HongwenZhang Qianget al2015.

A normal cloud model-based study of grading prediction of rockburst intensity in deep underground engineering

[J].Rock and Soil Mechanics,364):1189-1194.

Wang YuanhanLi WodongLi Qiguanget al1998.

Method of fuzzy comprehensive evaluations for rockburst prediction

[J].Chinese Journal of Rock Mechanics and Engineering,(5):15-23.

Xu ChenLiu XiaoliWang Enzhiet al2017.

Prediction and classification of strain mode rockburst based on five-factor crierion and combined weight-ideal point method

[J].Chinese Journal of Geotechnical Engineering,3912):2245-2252.

Yin XinLiu QuanshengWang Xinyuet al2020.

Prediction model of rockburst intensity classification based on combined weighting and attribute interval recognition theory

[J].Journal of China Coal Society,4511):3772-3780.

Zhang BiaoDai Xingguo2017.

rockburst intensity grade based on index distance and uncertainty measure

[J].Rock and Soil Mechanics,38Supp.2):257-265.

Zhang GuangcunGao QianDu Juqianget al2013.

Rockburst criterion based on artificial neural networks and nonlinear regression

[J].Journal of Central South University(Science and Technology),447):2977-2981.

Zhang JingjianFu Bingjun2008.

Rock burst and its criteria and control

[J].Chinese Journal of Rock Mechanics and Engineering,2710):2034-2042.

Zhang LewenZhang DeyongQiu Daohong2010.

Application of extension evaluation based on rough set in rockburst prediction

[J] .Coal Coal Journal,359):1461-1465.

ZhangYanSu GuoshaoYan Liubin2011.

Method of identifying rockburst grades based on Gaussian process machine learning

[J].Chinese Journal of Underground Space and Engineering,72):392-397.

Zou QiangZhou JianzhongYang Xiaolinget al2011.

A comprehensive assessment model for severity degree of dam failure impact based on attribute interval recognition theory

[J].Journal of Sichuan University (Engineering Science Edition),432):45-50.

陈祥孙进忠张杰坤2009.

岩爆的判别指标和分级标准及可拓综合判别方法

[J].土木工程学报,429):82-88.

[本文引用: 1]

贾哲强张茹张艳飞2016.

猴子岩水电站地下厂房岩爆综合预测研究

[J].岩土工程学报,38增2):110-116.

李航2020.

基于组合权重理想点法—数据库的隧洞岩爆预测方法研究

[D].北京中国地质大学(北京).

[本文引用: 3]

李鹏程叶义成王其虎2019.

基于正态白化权函数的灰评估岩爆预测模型

[J].化工矿物与加工,485):16-2226.

李任豪顾合龙李夕兵2020.

基于PSO-RBF神经网络模型的岩爆倾向性预测

[J].黄金科学技术,281):134-141.

[本文引用: 1]

林大涌雷明锋曹豪荣2018.

盾构下穿运营铁路施工风险模糊综合评价方法研究

[J].铁道科学与工程学报,155):1347-1355.

[本文引用: 1]

刘冉叶义成张光权2019.

岩爆分级预测的粗糙集—多维正态云模型

[J].金属矿山,(3):48-55.

刘晓悦杨伟张雪梅2021.

基于改进层次法与CRITIC法的多维云模型岩爆预测

[J].湖南大学学报(自然科学版),482):118-124.

[本文引用: 1]

裴启涛李海波刘亚群2013.

基于改进的灰评估模型在岩爆中的预测研究

[J].岩石力学与工程学报,3210):2088-2093.

[本文引用: 1]

孙臣生2019.

基于改进MATLAB-BP神经网络算法的隧道岩爆预测模型

[J].重庆交通大学学报(自然科学版),3810):41-49.

[本文引用: 1]

田睿孟海东陈世江2020.

基于深度神经网络的岩爆烈度分级预测研究

[J].煤炭学报,45增1):191-201.

[本文引用: 1]

王迎超靖洪文吉咸伟2014.

深埋地下工程岩爆烈度分级预测的RS-功效系数模型

[J].中南大学学报(自然科学版),456):1992-1997.

[本文引用: 1]

王迎超靖洪文张强2015.

基于正态云模型的深埋地下工程岩爆烈度分级预测研究

[J].岩土力学,364):1189-1194.

[本文引用: 1]

王元汉李卧东李启光1998.

岩爆预测的模糊数学综合评判方法

[J].岩石力学与工程学报,(5):15-23.

[本文引用: 2]

徐琛刘晓丽王恩志2017.

基于组合权重—理想点法的应变型岩爆五因素预测分级

[J].岩土工程学报,3912):2245-2252.

殷欣刘泉声王心语2020.

基于组合赋权和属性区间识别理论的岩爆烈度分级预测模型

[J].煤炭学报,4511):3772-3780.

[本文引用: 2]

张彪戴兴国2017.

基于指标距离与不确定度量的岩爆云模型预测研究

[J].岩土力学,38增2):257-265.

张光存高谦杜聚强2013.

基于人工神经网络及非线性回归的岩爆判据

[J].中南大学学报(自然科学版),447):2977-2981.

[本文引用: 1]

张镜剑傅冰骏2008.

岩爆及其判据和防治

[J].岩石力学与工程学报,2710):2034-2042.

[本文引用: 1]

张乐文张德永邱道宏2010.

基于粗糙集的可拓评判在岩爆预测中的应用

[J].煤炭学报,359):1461-1465.

[本文引用: 1]

张研苏国韶燕柳斌2011.

基于高斯过程机器学习的岩爆等级识别方法

[J].地下空间与工程学报,72):392-397.

邹强周建中杨小玲2011.

属性区间识别模型在溃坝后果综合评价中的应用

[J].四川大学学报(工程科学版),432):45-50.

[本文引用: 1]

/