地下工程岩体质量评判是进行工程设计、灾难控制的重要依据。为了提高评判效率,利用多粒度计算的方法对工程围岩稳定性进行质量评判。首先,对多粒度形式背景的粒度树上的属性块进行组合;其次,给出了计算广义介粒度剪枝形式背景权重的方法;最后,计算影响工程岩体质量的指标权重的方法,并用实例证明了方法的可行性和科学性。 Quality evaluation of underground engineering rock mass is an important basis for engineering design and disaster control. In order to improve the evaluation efficiency, the multi-granularity calculation method is used to evaluate the stability of surrounding rock. Firstly, attribute blocks in the granularity tree of multi-granularity formal background are combined. Secondly, the weight is taken as the criterion to evaluate the performance of optimal rock mass selection in underground engineering. Finally, the effectiveness of the weight-based optimal rock mass selection metric method is analyzed.
地下工程岩体质量评判是进行工程设计、灾难控制的重要依据。为了提高评判效率,利用多粒度计算的方法对工程围岩稳定性进行质量评判。首先,对多粒度形式背景的粒度树上的属性块进行组合;其次,给出了计算广义介粒度剪枝形式背景权重的方法;最后,计算影响工程岩体质量的指标权重的方法,并用实例证明了方法的可行性和科学性。
多粒度形式背景,岩体质量评判,粒计算,权重
Qian Liu1*, Baocang Du1,2, Hua Mao1#, Mengxuan Lian1, Chang Liu1, Zhiming Zhang1, Lanzhen Yang1
1University of Hebei, Baoding Hebei
2Hebei Finance University, Baoding Hebei
Received: May 12th, 2022; accepted: Jun. 16th, 2022; published: Jun. 23rd, 2022
Quality evaluation of underground engineering rock mass is an important basis for engineering design and disaster control. In order to improve the evaluation efficiency, the multi-granularity calculation method is used to evaluate the stability of surrounding rock. Firstly, attribute blocks in the granularity tree of multi-granularity formal background are combined. Secondly, the weight is taken as the criterion to evaluate the performance of optimal rock mass selection in underground engineering. Finally, the effectiveness of the weight-based optimal rock mass selection metric method is analyzed.
Keywords:Multi-Granularity Formal Concept, Rock Mass Quality Evaluation, Granular Computing, Weight
Copyright © 2022 by author(s) and beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
在工程施工中,为了给工程设计、支护衬砌、建筑物选型和施工方法选择等提供参数和依据,需对地下工程围岩进行质量评价。
目前,针对岩体质量评判,国内外有不少学者提出了许多建设性的方法 [
粒计算是1997年有模糊数学之父的Zadeh与Lin教授提出的一个研究课题,是一种看待客观世界的世界观和方法论 [
岩体质量分析是分析工程岩体地质条件、认识岩体工程特性,为工程地质灾害原因做出评判依据的重要准备工作。多粒度的思想更适合分析影响因素较复杂的岩体质量评判。本文以文献 [
本文第一节为引言,主要介绍研究背景以及国内外研究现状;第二节为预备知识,主要介绍本文所涉及到的定理和国家标准;第三节为本文的重点研究内容,通过多粒度形式背景的方法对地下工程岩体进行评判;第四节为工程实例;第五节为结论和展望。
本节将给出后文所涉及的有关概念理论和国家标准。关于信息粒化的具体详情见文献 [
本节相关多粒度形式背景的定义详见文献 [
定义2.1 形如C = (U, A)表示信息系统,其中U表示对象的非空集合;A表示属性的非空集合,对于任意的属性 a ∈ A ,存在一个对应值集Va,满足 a : U → V a ,在给定的信息系统C = (U, A)中,若对任意的属性 a ∈ A ,值集Va不为空,则称其为完备信息系统。
定义2.2 设(U, M, I)为形式背景,其中I是对象集U与属性集M之间的二元关系。Ia表示属性a拥
有的对象,对于B Í M,若 U = ∪ a ∈ B I a ,且 I a ( a ∈ B ) 两两不相交,则称B为(U, M, I)的单粒度类属性块。
定义2.3 设(U, Mi, Ii) ( i ∈ { 1 , 2 , ⋯ , n } )为n个单粒度形式背景,Bi1, Bi2, ∙∙∙, Bis是Mi的类属性块,其中B1k, B2k, ∙∙∙, Bnk( k ∈ { 1 , 2 , ⋯ , s } )是不同粒度下的同类别属性块。若(U, M1, I1) < (U, M2, I2) < ∙∙∙ < (U, Mn, In),
则称 π = ∪ i = 1 n ( U , M i , I i ) 为多粒度形式背景。
定义2.4 设 π = ∪ i = 1 n ( U , M i , I i ) 为多粒度形式背景,Bi1, Bi2, ∙∙∙, Bis是Mi的类属性块。令
M m = B n 1 1 ∪ B n 2 1 ∪ ⋯ ∪ B n s 1 ,其中 n 1 , n 2 , ⋯ , n s ∈ { 1 , 2 , ⋯ , n } ,ImÍ U × Mm,则称(U, Mm, Im)为介粒度形式背景。
定义2.5 设(U, M, I)为形式背景,对于 a ∈ M ,若属性a及其特化属性能够形成一棵树,则称它为a的粒度树,记做Ta。
定义2.6 设 π = ∪ i = 1 n ( U , M i , I i ) 可形成s棵n + 1层的粒度树,令Bi1, Bi2, ∙∙∙, Bis是Mi的类属性块,B1k, B2k, ∙∙∙, Bnk( k ∈ { 1 , 2 , ⋯ , s } 是第k棵粒度树上2至n + 1层的属性,且每个属性块包含的子属性块个数相同(记为γk), B 1 k 1 , B 2 k 2 , ⋯ , B n k γ k 为Bik的划分, D k = ∪ t = 1 γ k B n t k t 为第k棵粒度树上的多粒度剪枝类属性块。若
M g m p = ∪ k = 1 s D k , I g m p ⊆ U × M g m p
则称S = (U, Mgmp, Igmp)为π的一个广义介粒度剪枝形式背景。
定义2.7 广义介粒度剪枝形式背景的粒度和等于各个多粒度剪枝类属性块的粒度之和,且多粒度剪枝类属性块的粒度又等于其子属性块的平均粒度。S = (U, Mgmp, Igmp)的粒度和定义为
G S ( M g m p ) = ∑ k = 1 S ∑ t = 1 γ k n t γ k
其中{n1, n2, n3, ∙∙∙, nγk} Í {1, 2, ∙∙∙, n}。
地下工程围岩划分为5个等级,即稳定、较稳定、一般、不稳定、极不稳定,分别用罗马数字I、II、III、IV和V来表示。影响地下工程围岩稳定性的因素很复杂,不但取决于岩体内部的地质结构,还取决于外部环境,是多种因素的综合作用。
等级 | RQD/% | RW/MPa | Kv | Kf | ω/[L∙(min∙10m)−1] |
---|---|---|---|---|---|
I | 90~100 | 120~200 | 0.75~1.00 | 0.8~1.0 | 0~5 |
II | 75~90 | 60~120 | 0.45~0.75 | 0.6~0.8 | 5~10 |
III | 50~75 | 30~60 | 0.30~0.45 | 0.4~0.6 | 10~25 |
IV | 25~50 | 15~30 | 0.20~0.30 | 0.2~0.4 | 25~125 |
V | 0~25 | 0~15 | 0.00~0.20 | 0.0~0.2 | 125~300 |
表1. 地下工程岩体等级分类标准
根据已有研究成果以及有关的国家规范标准 [
首先,生成多粒度形式背景的粒度树。根据属性和粒化属性分布在不同的粒度层中,将属性之间的特化和泛化关系放到粒度树上进行研究,既而可以清楚分辨出属性之间的粒度层次关系。
给定一个信息系统C = (U, A),其中U = {x1, x2, x3, ∙∙∙, xn}表示对象集;A = {y1, y2, y3, ∙∙∙, ym}表示属性集。先将信息系统转换成多个单粒度形式背景,在将多个单粒度形式背景并置形成多粒度形式背景(U, M1, I1)。
在工程中,影响地下工程岩体质量因素的结构面强度系数Kf,可以从结构面的抗剪强度、结构面的抗拉强度这两部分进行研究。通常结构面的抗拉强度非常小,所以抗剪强度和抗拉强度对Kf的影响所占权重也应不同。同样,结构面的抗剪强度又受摩擦角、粘聚力等剪应力作用的影响。
例1 给定一个信息系统C = (U, A),其中U = {1, 2, 3}表示地下工程岩体质量评判对象;A = {Kf, ω}表示影响地下工程岩体质量评判的指标,Kf为结构面强度系数;ω为地下水渗水量。见表2。
U | Kf | ω/[L∙(min∙10m)−1] |
---|---|---|
1 | 0.45 | 0 |
2 | 0.55 | 10.5 |
3 | 0.52 | 12.0 |
4 | 0.38 | 21.0 |
5 | 0.30 | 18.5 |
表2. 信息系统C = (U, A)
根据数据表的阈值,可将信息系统中的多值属性Kf分成2个子类:“[0.30, 0.45)”和“[0.45, 0.60)”,分别记为a1和a2;将多值属性ω分成2个子类:“[0, 15)”和“[15, 25)”,分别记为a3和a4。则信息系统C = (U, A)可转化为形式背景(U, M1, I1),见表3。其中,若对象满足该属性则,横纵坐标的交叉处为1;若对象不满足该属性则为0。
U | a1 | a2 | a3 | a4 |
---|---|---|---|---|
1 | 0 | 1 | 1 | 0 |
2 | 0 | 1 | 1 | 0 |
3 | 0 | 1 | 1 | 0 |
4 | 1 | 0 | 0 | 1 |
5 | 1 | 0 | 0 | 1 |
表3. 形式背景(U, M1, I1)
此外,根据数据本身,还可把a1,a2,a3和a4进一步细分,将a1细分为“[0.30, 0.37)”和“[0.37, 0.45)”,分别记为b1和b2;将a2细分为“[0.45, 0.52)”和“[0.52, 0.60)”,分别记为b3和b4;将a3细分为“[0, 5)”、“[5, 10)”和“[10, 15)”,分别记为b5、b6和b7;将a4细分为“[15, 20)”和“[20, 25)”,分别记为b8和b9。形式背景则转化为表4。
U | b1 | b2 | b3 | b4 | b5 | b6 | b7 | b8 | b9 |
---|---|---|---|---|---|---|---|---|---|
1 | 0 | 0 | 1 | 0 | 1 | 0 | 0 | 0 | 0 |
2 | 0 | 0 | 0 | 1 | 0 | 0 | 1 | 0 | 0 |
3 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 0 |
4 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 1 |
5 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 0 |
表4. 形式背景(U, M2, I2)
针对表2和表3中的数据,显然有 I a 1 ∪ I a 2 = I a 3 ∪ I a 4 = U 。因此表3有两个类属性块,分别记为M1= {a1, a2},M2= {a3, a4, a5};同样,表3中也有两个类属性块M3= {b1, b2, b3, b4},M4= {b5, b6, b7, b8, b9}。其中, I a 1 = I b 1 ∪ I b 2 , I a 2 = I b 3 ∪ I b 4 , I a 3 = I b 5 ∪ I b 6 ∪ I b 7 , I a 4 = I b 8 ∪ I b 9 ,有M1< M3和M2< M4,即
(U, M1, I1)和(U, M2, I2)之间形成了粒度粗细关系,所以表3和表4可并置形成多粒度形式背景 π = ∪ i = 1 2 ( U , M i , I i ) 。
由此可知π是由2个单粒度形式背景组成的,且每个单粒度形式背景均有2个类属性块,所以π可形成2棵3层的属性粒度树,如图1和图2。
图1. 结构面强度系数Kf的属性粒度树Ta1
图2. 地下水渗水量ω的属性粒度树Ta2
在对属性进行分块时,如果分块太细,则会增加时间复杂度;如果分块太粗,则不能很好的将数据和实际问题进行匹配。所以,分块的粗细,即信息粗细程度的刻画很重要。
定义3.1 广义介粒度剪枝形式背景(U, Mgmp, Igmp)的信息粒度定义为
D ( M g m p ) = | U | ∑ y ∈ U | U | | y ∇ ∇ | . (1)
其中,广义介粒度形式背景(U, Mgmp, Igmp)的所有属性概念为 { y ∇ , y ∇ ∇ | y ∈ U } 。由定义可知,属性粒度树剪取越细的分支组合,信息粒度的值就越大。
定义3.2 广义介粒度剪枝形式背景(U, Mgmp, Igmp)的权重定义为
W ( M g m p ) = 1 | U | ∑ y ∈ U ( | U | | y ∇ ∇ | ) . (2)
定义3.3 设(U, M g m p α , I g m p α )和(U, M g m p β , I g m p β )为多粒度形式背景 π = ∪ i = 1 n ( U , M i , I i ) 的两个广义介粒
度剪枝形式背景,若 G S ( M g m p α ) = G S ( M g m p β ) 且 W ( M g m p α ) > W ( M g m p β ) ,则称(U, M g m p α , I g m p α )的粒度组合优于(U, M g m p β , I g m p β )。
也就是说权重相同的广义介粒度剪枝形式背景,信息粒度越大的粒度组合越优。
例2 在图1和图2中,结构面强度系数Kf粒度树Ta1的剪枝方式有4中:B11= {a1, a2}, B12= {a1, b3, b4},B13= {b1, b2, a2},B14= {b1, b2, b3, b4};地下水渗水量ω粒度树Ta2的剪枝方式也有4中:B21= {a3, a4},B22= {a3, b8, b9},B23= {b5, b6, b7, a4},B24= {b5, b6, b7, b8, b9}。先将他们两两组合可得到16种广义介粒度剪枝形式背景,分别记为 S μ = ( U , M g m p μ , I g m p μ ) 。其各广义介粒度形式背景的权重和粒度和见表5。
剪枝背景 | 粒度和 | 权重 | 剪枝背景 | 粒度和 | 权重 |
---|---|---|---|---|---|
S1 | 2 | 1/32 | S9 | 2.5 | 1/100 |
S2 | 2.5 | 1/100 | S10 | 3 | 1/144 |
S3 | 2.5 | 1/144 | S11 | 3 | 1/196 |
S4 | 3 | 1/112 | S12 | 3.5 | 1/256 |
S5 | 2.5 | 1/100 | S13 | 3 | 1/100 |
S6 | 3 | 1/144 | S14 | 3.5 | 1/112 |
S7 | 3 | 1/196 | S15 | 3.5 | 1/256 |
S8 | 3.5 | 1/256 | S16 | 4 | 1/324 |
表5. 图1和图2的各广义介粒度形式背景的权重和粒度和
由定义可知当粒度和为2或4时,广义介粒度剪枝形式背景S1或S16的权重最大,即为S1和S16为该粒度和下的最优剪枝形式背景;当粒度和为2.5时,S2、S5和S9为该粒度和下的最优剪枝形式背景;当粒度和为3时,S13为最优剪枝形式背景;当粒度和为3.5时,S14为最优剪枝形式背景。
广义介粒度剪枝形式背景在岩体质量评价中的工程意义:对目标岩体的各项指标通过系统的分类计算方法,得到最优剪枝形式背景,即影响岩体质量的每个指标的最优分类区间。首先,通过生成属性粒度树的方法,建立广义介粒度剪枝形式背景;其次,给出计算广义介粒度剪枝形式背景的权重,来得到最优剪枝形式背景;最后,根据其工程意义,找到判断岩体质量的最优剪枝形式背景,既而得到岩体质量每个指标的最优分类区间。
在本文方法中,给出了计算广义介粒度剪枝形式背景权重的方法,在工程上即为,计算影响工程岩体质量的指标权重的方法;在其他文献当中,或者是通过约简概念格的方法对属性进行约简,将影响岩体质量的指标缩减;或者是通过其他新模型在岩体质量各项指标的测量上进行优化,从而达到简化数据处理的目的;或者只是针对某些特定的地区、地形或地貌使用,在指标的权重上相关研究较少。
本文通过工程岩体的广义介粒度形式背景,计算其权重,从而得到影响岩体质量每个指标的最优分类区间,将分类具体化,通过对各项数据细分的方法,得到具体各项指标细分的程度和细分的最优区间。
通过多粒度形式背景的方法,计算地下工程岩体质量指标的权重。
工程意义:希望通过计算广义介粒度形式背景的权重的方法,达到计算最优工程岩体质量指标分类的方法。
输入:地下工程岩体样本质量指标样本数据;
输出:影响地下工程岩体样本质量指标的权重。
第一步:通过对岩体样本数据细分的方法,建立专属地下工程岩体的多粒度形式背景;
第二步:根据多粒度形式背景,生成地下工程岩体的属性粒度树;
第三步:根据多粒度属性块的不同组合方式,生成基于粒度树的广义介粒度形式背景;
第四步:由广义介粒度形式背景,计算每个介粒度形式背景的权重。
采用文献 [
编号 | RQD/% | RW/MPa | Kv | Kf | ω/[L∙(min∙10m)−1] |
---|---|---|---|---|---|
样本1 | 71.8 | 90.1 | 0.57 | 0.45 | 0 |
样本2 | 51.0 | 40.2 | 0.38 | 0.55 | 10.5 |
样本3 | 52.0 | 25.0 | 0.22 | 0.52 | 12.0 |
表6. 地下工程岩体学习样本数据
下面将用本文的方法实现表6中地下工程岩体的质量评判。
首先,生成多粒度形式背景的粒度树。
根据每一列数据阈值,可将RQD细分为“[40, 65)”和“[65, 90)”,分别记为e1和e2;RW细分为“[25, 50)”和“[50, 100)”,分别记为e3和e4;Kv细分为“[0.20, 0.40)”和“[0.40, 0.60)”,分别记为e5和e6;Kf细分为“[0.40, 0.50)”和“[0.50, 0.60)”,分别记为e7和e8;ω细分为“[0, 10)”和“[10, 15)”,分别记为e9和e10。若样本i (i = 1, 2, 3)满足该项指标ej(j = 1, 2, 3, ∙∙∙, 10)记为“1”,若不满足该项指标记为“0”。将样本集记为对象集U;将指标集记为属性集M,两者之间的关系为I,见表7形式背景(U, M, I)。
U | e1 | e2 | e3 | e4 | e5 | e6 | e7 | e8 | e9 | e10 |
---|---|---|---|---|---|---|---|---|---|---|
1 | 0 | 1 | 0 | 1 | 0 | 1 | 1 | 0 | 1 | 0 |
2 | 1 | 0 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | 1 |
3 | 1 | 0 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | 1 |
表7. 形式背景(U, M, I)
根据每一列数据阈值,仍可以继续细分,将e1细分为“[40, 52)”和“[52, 55)”,分别记为f1和f2;将e2细分为“[65, 75)”,分别记为f3;将e3细分为“[25, 40)”和“[40, 50)”,分别记为f4和f5;将e4细分为“[90, 100)”,记为f6;将e5细分为“[0.20, 0.30)”和“[030, 0.40)”,分别记为f7和f8;将e6细分为“[0.50, 0.60)”,分别记为f9;将e7细分为“[0.45, 0.50)”,分别记为f10;将e8细分为“[0.50, 0.55)”和“[0.55, 0.60)”,分别记为f11和f12;将e9细分为“[0, 5)”,分别记为f13;将e10细分为“[10, 12)”和“[12, 12)”,分别记为f14和f15。见表8,记为形式背景(U, M', I')。
U | f1 | f2 | f3 | f4 | f5 | f6 | f7 | f8 | f10 | f11 | f12 | f13 | f14 | f15 |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
1 | 0 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 1 | 0 | 0 |
2 | 1 | 0 | 0 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 1 | 0 | 1 | 0 |
3 | 0 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 0 | 1 |
表8. 形式背景(U, M', I')
由表7和表8可得多粒度形式背景 π = ∪ i = 1 5 ( U , M i , I i ) 的属性粒度树,如图3。
图3. 岩石质量指标RQD的粒度树Ta3
其余4个指标属性粒度树形式和图3形式相同,这里便不再一一列举。
按预备知识中对多粒度类属性块的定义,可将图3岩石质量指标RQD的属性粒度树剪枝组合为{e1, e2},{e1, f3},{e2, f1, f2},{f1, f2, f3} 4个多粒度剪枝类属性块。
计算表6中的形式背景的信息粒度,由定义可知,表7的广义介粒度剪枝形式背景(U, M g m p 1 , I g m p 1 )是由图3岩石质量指标RQD的属性粒度树的第2层,以及其余4个指标的属性粒度树的第2层进行属性组合得到的。
图3的岩石质量指标RQD的粒度树Ta3的剪枝方式共有4种:B31= {e1, e2},B32= {e1, f3},B33= {f1, f2, e2},B34= {f1, f2, f3};同理,其余4个指标的粒度树Tai(i = 4, 5, 6, 7)的剪枝方式也分别有4种;它们两两组合便能形成表6的广义介粒度剪枝形式背景。
例:当 S 1 = ( U , M g m p 1 , I g m p 1 ) ,其中 M g m p = ∪ k = 1 5 D k ; D k = ∪ t = 1 10 B n t k t ; I g m p ⊆ U × M g m p 时, e 1 ∇ ∇ = e 3 ∇ ∇ = e 5 ∇ ∇ = e 8 ∇ ∇ = e 10 ∇ ∇ = { e 1 , e 3 , e 5 , e 8 , e 1 0 } , e 2 ∇ ∇ = e 4 ∇ ∇ = e 6 ∇ ∇ = e 7 ∇ ∇ = e 9 ∇ ∇ = { e 2 , e 4 , e 6 , e 7 , e 9 } ,则 D ( M g m p ) = 10 ( 10 5 × 10 ) = 200 ; W ( M g m p ) = 1 200 。
即广义介粒度剪枝形式背景:
S 1 = ( U , M g m p 1 , I g m p 1 ) 的 W ( M g m p ) = 1 200 。
本文方法在复杂度上虽没有大幅度提升,但在计算工程岩体各项质量指标上提出了新颖的计算思路和方法,即通过计算工程岩体广义介粒度剪枝形式背景权重的方法,对工程岩体各项指标的权重有个整体的分析,但在计算具体每项指标的具体权重上没有给出相关定义,这将在未来的工作中体现。相关工程施工方面的知识,期待更多相关专业人士的加入和指导。
河北大学2022年研究生创新资助项目(HBU2022bs028)。
刘 谦,杜宝苍,毛 华,连萌璇,刘 川,张植明,杨兰珍. 岩体质量评判的多粒度计算方法Multi-Grained Calculation Method for Rock Mass Quality Evaluation[J]. 理论数学, 2022, 12(06): 986-995. https://doi.org/10.12677/PM.2022.126108
https://doi.org/10.1016/j.enggeo.2021.106382
https://doi.org/10.1016/j.enggeo.2021.106346
https://doi.org/10.1007/s10064-018-1419-y
https://doi.org/10.1007/s12517-021-06537-1
https://doi.org/10.1155/2021/6652960
https://doi.org/10.1016/j.ins.2021.01.027
https://doi.org/10.1016/j.engappai.2020.103694
https://doi.org/10.1016/j.knosys.2021.106880