与设计ResearchandDesing
基于三角模糊云贝叶斯网络的盾构下穿段施工风险评估
()中铁投资集团有限公司,北京1石家庄铁道大学道路与铁道工程安全保障教育部重点实验室,河北石家庄01.00160;2.50043
陈洪波1, 王 宁2, 王建西2
摘 要:为进行不确定性数据下地铁下穿施工风险评估,构建了基于三角模糊云贝叶斯网络(的风险分析BN)模型.将三角模糊数与云理论相结合,提出了定性评判向定量概率分布转化的方法.引入层次分析法(和状态影响因子得到云数字特征,生成随机风险云与统计方法计算节点条件概率表CAHP)PT.以实际工程为依托,建立了风险评价指标体系和网络模型,分析了施工风险水平和因素敏感性.研究结果与实际相吻合,验证了模型的合理性和有效性.
关键词:盾构下穿;风险评估;三角模糊数;云模型;贝叶斯网络:/DOI10.13219.at.2020.03.004jgjgy
()中图分类号:U458.1 文献标识码:B 文章编号:1672G3953202003G0013G006
近年来,地铁施工引起的既有老旧建筑沉降、开裂、倾斜等破坏对城市公共安全构成极大威胁.前人采用传统方法针对不同情况下地铁邻近建筑物的施工安全问题开展了大量研究.然而,地铁施工风险因素往往表现出显著的不确定性和动态性特征.贝叶斯网络在描述事件不确定性、知识表达和动态推理方面具有显著优势.相关研究虽然取得了一定成果,但施工中大量风险因素往往难以监测,获取的风险数据包含较大模糊不确定性.三角模糊数和云模型是合理考虑评估数据的模糊性与不确定性的有效方法.
本文首先将三角模糊数与云理论相结合,研究定性评判向定量概率分布转化的方法.然后,基于层次分析法、状态影响因子和云模型得到贝叶斯网络条件概率.最后,构建施工风险评价贝叶斯网络模型,为盾构下穿施工安全风险预控提供理论依据.
云,x)被称为确定度,x表示该云的一个云滴.μ(
期望Ex、熵En和超熵He是云的3个数字特征.,N(EHe2)x对Z的确定度满足:n,
-
2()xE-x22E′n2
),若x为正态云的云滴,则x~N(其中EEE′′~x,nn()x)=e1μ(
期望Ex表示云滴分布的中心点,反映了对定性概念认知的统一性和稳定性.熵En是定性概念随机不确定性的度量,反映云滴分布具有模糊特征.超熵其值越大则云层越厚,He的直观反映是云滴厚度,
它由熵的随机性和模糊性共同决定.
定量的随机数据,即产生一定数量的云滴.相反,逆向云发生器可以统计得到云的数字特征,从而进行ìïEx=1xiniï=1
ïn1|xi-Ex|ïE()2ín=ni=1lnx-2i)μ(ï
ïn1i-Ex|2ïH(|x-Ee=n)ïn-1i=1lnxî-2i)μ(1.1.2 风险评价的三角模糊数表示
实践中,模糊随机事件很难用确定的数值量化,为此引入三角模糊数描述事件风险概率等级.
定义2:设论域C上的模糊集为Z,若μx)∈Z()定性评价,如式(所示:2
n),利用正态云发生器式(可将定性概念转换为1
1 三角模糊云模型与贝叶斯网络
1.1 基于三角模糊数的云模型建立
1.1.1 云模型相关概念
云模型是大数据和人工智能领域中进行不确定
[]1
:定义1设对象X的定性论域为C,Z是论域
性分析的强有力工具,可实现定性与定量的相互转换.上的一个定性概念集合,若对于C上的一次随机定],量实现x,有μ(则x在C上的分布称为x)0,1∈[
收稿日期:2019G12G18
,陈洪波(男,高级工程师,主要从事地下 第一作者简介:1977—)工程施工技术与风险管理工作.4110744@q.comq
[]是x在Z上的映射函数,表示为上限u、下限l0,1
和模态值m的线性函数,称为三角模糊隶属函数.
()国防交通工程与技术 1020年5月 18,033 2
研究与设计 基于三角模糊云贝叶斯网络的盾构下穿段施工风险评估 陈洪波等
ì,ïx-lx∈[lm]ml-ïï
()x)=íx-u3Z(μ[,]xmu∈ïm-uïï其他î0
本文依据«城市轨道交通地下工程建设风险管理规范»定义的风险发生概率P特征集将风险事件概率等级分为5级,表1给出了各标准等级的分值区间[及三角模糊数.CmCmain,x]
标准分级等级描述
表1 风险指标标准等级划分与三角模糊数
不可能
用于量化事件间的相互作用.通常,CPT记为P,在无法获取风险指标大量数据时需借助专家经验建模.为确保合理性,本文利用层次分析法和状态影响因子评判风险指标的相互影响,再利用云模型生成条件概率表CPT.
考虑到风险指标(父节点)的不同等级状态对下
一层指标(子节点)的影响程度不同,通过专家评判法分别确定每个父节点不同状态对子节点的影响值],得到状态影响因子矩阵S=S0,10ij∈[
(将n个父节点的各个状态进行遍历,依据Sin5.×j)计算公式为:
ì=WiSijïExli1=ï
í/Cma6n=(x-Cmin)ïEï
îHe=knP1
1.1.3 三角模糊数的云模型转化
利用三角模糊数能够方便地实现被评价对象的度量,使定性语言值转化为具体的数量,然而,模糊运算规则处理评价数据时容易导致信息损失.云模型通过对生成的大量云滴进行统计分析,具备强大的模糊不确定性数据推理能力.为此,将三角模糊数转化为云模型,实现风险事件不确定性数据运算.,在区间[①生成数据样本.利用公式(3)lu]内随机抽样,产生包含三角模糊数所表达的不确定).和模糊信息的n个样本点(xx②计算云数i,Z(i)μ③利用正态云发生器生成云.
[CmCm~22~44~66~88~10in,ax]0
)()()()()x)(0,1,22,3,44,5,60,1,20,1,2Z(μ罕见
P2
偶尔
P3
可能
P4
频繁
P5
节点权重和状态影响因子得到云模型的数字特征,
()5
,);式中:1,2,5j为节点的状态循环变量,j=(
l为父节点的状态组合个数;k反映了专家评估随机性,结合实际情况选取k=0.根据均值Ex、熵E2.n和超熵He,利用云发生器和概率分布计算方法即可确定节点的条件概率表.1.2.2 网络推理与敏感性分析
)字特征,利用公式(得到均值E熵E1x、n和超熵He.1.1.4 概率分布的计算
在云模型中,云滴构成的云团对论域中每个定性概念的贡献是不同的,云团越密集则对该概念的支持也越大.设一维论域C的任一小区间上的云团为△x,它对概念Z的贡献为:
/(2()△x△A≈μx)πEn)4Z(显然,上式在坐标系(上的积分等-∞,+∞)
],令A于1.对于论域C=[0,10i=
ui)、)测推理和诊断,如式(式(所示.67
贝叶斯网络可在风险分析的任意阶段,进行预
P(B=b)=
合.假设目标节点变量为Q,证据节点变量为E,利),,如式(其值越大则对相关因素越MutualInfo8)
敏感.
P(XB=b)=i=xi|P(XP(B=b|Xi=xi)i=xi)()7
P(B=b)
敏感性分析用于确定风险防控的关键因素集
,,Xn)P(X1,X2,Xn)
B=b|XP(
1
,X2,
()6
:用信息熵计算出节点变量之间关联的互信息(MI
Q)P(
Q/(2(,).可得,对x)πE)dx,i=1,2,5
∫μ(
liZn1.2 贝叶斯网络模型
/i1A,).Ai=1,2,5i=i(
5
应标准分级“的概率分布为PPPPPP1、2、3、4、5”i=
I(Q,E)=H(Q)Q|E)=-H(
1.2.1 条件概率表建立
2 地铁下穿建筑群风险评估模型
2.1 建立风险评价指标集
],根据专家咨询并参考文献[识别盾构下穿建2
11
loQ,E)lo-P(g()g(
PQPQ|E)Q,E()8
,节点表示事件,构成网络图模型G,条件概率表P)
贝叶斯网络可记为有向非循环图BN=(G,
()国防交通工程与技术 1020年5月 18,034 2
研究与设计 基于三角模糊云贝叶斯网络的盾构下穿段施工风险评估 陈洪波等 筑群的主要风险因素,按照影响关系建立4层评价指标体系,如图1所示.
M1RTM2M3表2 风险因素权重及状态影响因子评价表父节点
子节点
权重0.54990.24020.20980.350.650.550.450.700.300.600.12600.41610.45790.29700.16340.53960.350.650.40
M1M2M11M12M21M22M31M32X1X2X3X4X5X6X7X8X10X11X12X13X14X15X9M3
P1P2P3P4211120121202100220021001
321243232413233441142212
543566464644656663456535
865799778987999886879869
P5
RT1098101010899
M11M12M21M22M31M32M1M2M3M11M12
X8X13X15X3X6X9X10X7X11X12X1X2X4X5X14101010101010101010101010888
图1 盾构下穿密集老旧建筑风险评估指标
评估模型中,每个指标因素即为一个风险事件,),按照所提出的三角模糊云模型方法进行定量2011
评价.
2.2 网络模型结构与CPT确定
将风险评价指标集X1~X15个因素作为5的1节点,建立网络模型.根据表2所示的指标权重和其中,指标权重采用AH状态影响因子评价由P法,专家评判,以0~10给出不同状态下父节点对子节点的影响评分.
,重要性进行两两比较,分析得到权重W=[0.5499
].0.2402,0.2098将RT父节点各状态进行遍历,依据节点权重
和状态影响因子计算云期望EX=[,从而得到云数字特征.例如,,EX1,EX2,EX5]EX1=
以节点RT为例,将3个父节点M1、M2和M3的状态影响因子评价值,利用式(计算云数字特征.5)(可参考«地铁工程施工安全评价标准»GB50715-
M21M22M31M32
0.14290.28570.57140.400.60
计为盾构法,K28+413.349~658.416段下穿平低压缩性.地下水类型为潜水,埋深约3位于8m,隧道底板以下.
施工区周边分布多个建筑群.其中,平安小区有8栋老旧建筑物,经鉴定房屋均为D危险u级(户)建于2条形石砌基础,无圈0世纪五六十年代,管控的难点.为此,选取平安小区2#楼作为分析对象进行动态风险评估.
石家庄地铁2号线东三教站~东岗头站区间设
安小区,隧道顶距建筑物基础15.5~16.5m.该段地层为粉细砂、粉质黏土,砂层中密~密实,稍湿,中~
0.5499×2+0.2402×1+0.2098×1=1.549.利用云发生器生成10经过统计分布计00个云滴,算条件概率表.由于RT有3个父节点,每个节点有
4
因此,节点C个元5个状态,PT共包含125行5素,表3仅给出前10行分布计算结果.按以上方法计算所有节点的条件概率表.模型建立后,即可进行各评价指标实时更新和风险预测,为下穿施工安全提供决策信息.
.其中1#、级)均为3层房屋,共12#、9#楼(31梁、结构柱等抗变形结构,保证下穿施工安全是风险
3.2 证据数据收集与处理
编制施工风险指标评分表,邀请业内专家逐项评价.通过勘察、施工和监测资料,由专家给出三角模糊评价值,如表4所示.采用1.1.3节的方法将其转化为云模型,计算数字特征和各风险因素的概率分布结果.
3 案例研究
3.1 工程背景
()国防交通工程与技术 1020年5月 18,035 2
研究与设计 基于三角模糊云贝叶斯网络的盾构下穿段施工风险评估 陈洪波等
表3
节点状态
基于状态影响分析结果和云模型生成节点RT的CPT
M1M2M31234512345
1111122222
1111111111
Ex1.54982.09973.19954.84925.94901.79002.33993.43975.08946.1892
En0.330.330.330.330.330.330.330.330.330.33
He0.20.20.20.20.20.20.20.20.20.2
P10.70130.4417
000.62680.25020.0066
000
P20.29870.55830.89260.08040.37320.74980.79220.0162
00
P300.10740.81940.4764
000
P4000.10020.5236
000.10500.6389
00
P50000000000
0.20110.87880.3611
表4 风险因素专家评分结果专家1
专家2
风险因 三角模糊数( L,M,U)素代号
X1X2X3X4X5X6X7X8X10X11X12X13X14X15X9
(,,)(,,)(,,)2.03.04.02.53.54.51.53.23.8
(,,)(,,)(,,)00.20.300.51.000.20.5(,,)(,,)(,,)3.04.04.54.55.05.34.85.25.5(,,)(,,)(,,)6.06.57.05.05.56.06.57.07.5(,,)(,,)(,,)4.05.06.04.55.05.54.04.55.0(,,)(,,)(,,)3.54.04.53.64.35.02.03.55.0(,,)(,,)(,,)1.52.53.51.02.02.51.22.53.2(,,)(,,)(,,)1.02.03.01.52.53.00.52.03.0(,,)(,,)(,,)2.02.83.82.53.54.02.53.03.5(,,)(,,)(,,)0.01.52.00.00.52.00.01.01.5(,,)(,,)(,,)2.02.53.01.52.02.51.01.52.0(,,)(,,)(,,)3.04.05.03.24.05.03.34.05.2(,,)(,,)(,,)0.50.81.01.52.53.51.82.43.0(,,)(,,)(,)2.02.53.01.82.23.02.0.2.22.5(,,)(,,)(,,)1.02.23.21.62.53.52.02.32.8
专家3
3.3.1 风险发生概率分析
图3为不可观测风险因素节点各状态等级的分布概率堆叠柱状图.为便于比较,将节点状态概率计算综合风险概率()与状态等级Istateii加权求和,p等级I绘于图3中.ave,
云.根据综合云统计云团对风险概率等级的贡献度,].按照相同的方法计算其余因0.2370,0.0770,0
素的概率分布,将上述节点概率数据输入贝叶斯网络进行推理.
10.8:※(1);+(2);×(3);○图2给出了风险指标X1的三角模糊云及综合
,,得到X0.21710.46901的概率分布计算结果为p=[
图3 不可观测节点各状态等级分布概率堆叠柱状图
和M3中,观察曲线和柱状图可M1发生概率最高.
通过I在二级风险因素M1、M2ave值比较发现,
知,潜在风险隐患主要来源于M12周边建筑环境.
该下穿施工区域内建筑结构整体性较弱、建筑密集且年代久远,近距离施工增加了事故的可能性.的影响,导致M2盾构M21的风险发生概率大于M22,“况的风险概率处于较低水平,不可能”和“罕见”的
由于注浆量控制X6和土仓压力参数不合理X8
0.60.40.200施工参数风险概率增加.M3施工组织管理与应急状
246
8103.3 结果分析
图2 风险因素X1专家评判三角模糊云和综合云
,分布概率之和超过了“综合风险概率等级96.09%”
.,为1.6539M3.6241施工管理综合风险概率等级为1
高于M32监控与应急.
()国防交通工程与技术 1020年5月 18,036 2
研究与设计 基于三角模糊云贝叶斯网络的盾构下穿段施工风险评估 陈洪波等 3.3.2 盾构下穿施工风险评价
依据«规范»定义风险损失C特征集ΩC=[可
,忽略,需考虑,严重,非常严重,灾难性]分别以)表示.利用风险矩阵得到5级风险分Cl=1~5l(
级标准,如表5.为了便于施工风险决策,结合工程,风险管控实际,参考文献[建立表6所示的风险3]可接受准则.
表5 风险矩阵及分级
风险级别
风险损失等级
技术措施将其施工风险转移并控制在可接受范围.施工安全事故常与作业中的不当操作、技术措施不合理等人为因素相关,X12班组安全培训和考核的——“风险评估结果略超过Ⅲ级—可接受”水平,施工中
应强化作业班组安全培训,降低人为操作失误风险.按风险大小依次为X6、X5、X9、X8、X3、X1X14、5和
风险等级接近Ⅲ级的因素共8个,约占53.3%.
X7.X5近接施工距离和X3基础埋深与结构类型由客观设计因素决定,其风险处于可接受范围.在技术方面,X6注浆量控制是减少地层损失和预防沉降的出现X9刀盘和刀具磨损可能性也较大,这与砂卵石致.由于地表建筑附加荷X8为土仓压力参数设置,载作用,下穿时应做好优化调整,保证合适的土压力.此外,还应对X1X14危险源监控、5应急处置与风险控制、X7推进速度不合理因素做好控制防范.
盾构下穿风险节点RT的风险概率等级后验分关键,需要采取有效措施降低风险.分析结果揭示,类地层长距离掘进易发生刀具磨损的工程经验相一
风险P4
概率P3等级P2
P5
C5
R5R4R4R3R3
C4
R4R4R3R3R2
C3
R4R3R3R2R1
C2
R3R3R2R2R1
C1
R3R2R2R1R1
P1
表6 风险接受准则
风险等级Ⅴ级Ⅳ级Ⅲ级Ⅱ级Ⅰ级
接受准则及风险控制对策拒绝接受,立即停止施工,进行整改,
规避或启动应急预案不可接受,需决策,制定控制、预警措施可接受,需引起重视,做好防范、监控
可容许,加强管理和巡查可忽略,日常管理和审视
布p=[1.242%,30.345%,57.133%,11.253%,
.由于潜在风险损失等级为C非常严0.027%]4(2+0.3034×3+0.5713×3+0.1125×4+0.00027)略超Ⅲ级.根据风险可接受准则,×4=3.1004,盾构下穿建筑群段施工风险总体可接受,但需引起重视并做好防范措施.
3.3.3 敏感性分析通过Netica软件计算盾构下穿建筑群风险节点RT与所有父节点的互信息,其值越大则对RT影响越敏感,见表8.按照因素敏感性由高到低、风险施工时X5、X1和X3无法规避外均为可控风险因素,应优先采取针对性的处置措施.
16X11X214X1X13X1012XX71510X148X8X36X9X5XX122X400246810121416,重)取表5第2列计算综合风险水平R=(0.0124×
失专家评估,结果见表7.利用表5,将各因素风险损失对应列的概率分布与风险等级相乘后求和即可得到风险等级大小,如图4所示.
X1X2X3X4
C3C3C3C4
X5X6X7X8
C3C4C4C4
X9
表7 风险损失等级专家评价
针对基本风险因素和施工总体风险进行风险损
因素损失等级因素损失等级因素损失等级因素损失等级
X10X11X12
C4C4C3C4
X13X14X15RT
C4C4C4C4
等级由大到小的顺序排列,得到图5所示的分布图.显然,越靠近左下角则该风险因素越重要,除X4、
53.920442.71162.97562.6292.87523.05562.691732.9972.7562.09712.36252.3682.6815211.50410X1X2X3X4X5X6X7X8X9X10X11X12X13X14X15占总数1由于施工区域地表X13.3%.分析可知,2,
住宅建于上世纪中期,为使用超过60a的老旧建筑,需采取X4建筑使用年限的风险等级接近Ⅳ级,
由图5可见,风险等级超过Ⅲ级的因素为X4和
图4 基本风险因素风险等级评估结果
图5 敏感性与风险排序分布图
4 结论
()将三角模糊数和云理论相结合,提出了定性1
()国防交通工程与技术 1020年5月 18,037 2
研究与设计 基于三角模糊云贝叶斯网络的盾构下穿段施工风险评估 陈洪波等
表8 盾构下穿诱发事故风险敏感因素分析节点变量
/100M.I.
节点变量
/1003.983.192.751.87M.I.
()引入层次分析与状态影响因子成功考虑了2
父节点权重和不同状态对子节点的影响,并通过生成随机风险云和统计分析方法建立条件概率表.()构建了盾构下穿施工风险评价指标体系和3
三角模糊云贝叶斯网络模型.应用表明,该模型的风险推理结果与实际吻合,为现场提供了科学准确的决策信息.
参考文献
[]王景春,张 法,林佳秀.基于改进云模型的公路施工安全1
[]王 帆,覃亚伟.地铁区间施工周边环境安全风险评估方2
[,,,3]WANGN.XUC.S.DUX.L.et.al.AriskassessG
sisandexertelicitation[J].InternationalJournalofp,SstemAssuranceEnineerinndManaement2018,9yggag():02452G466.
mentmethodofdeexcavationbasedonBaesiananalGpeyy]():,法研究[土木工程与管理学报,J.20130491G96102.]():防护系统评价[中全科学学报,J.20170873G78.
盾构下穿风险1RX438T:3埋深与结构类型
M1地质与周边环境5M3561施工管理M1211X2周边建筑环境15应急处置与风险控制M1X411工程水文地质6注浆量控制
M2施工参数与设备X4建筑使用年限X1地质条件复杂M21盾构施工参数
36X.4813施工管理措施与监督0
M2360.472状态与适应性X240.4512安全培训考核20111164
M3施工管理与应急
X8土仓压力不合理X5近接施工距离M32监控与应急
X7推进速度不合理
X10刀盘主轴密封失效X9刀盘与刀具磨损X14危险源监控X11班组经验水平
0.250.210.110.03
评判转化为定量概率分布的方法.该方法充分利用了二者在模糊定性表达与定量转化方面的优势,为风险不确定性数据处理提供了有利工具.
0.02
OntheTrianularFuzzloudBaesianNetworkGBasedRiskAssessmentofgyCy
theShieldGTunnelinonstructionUnderGCrossinxistinuildinsgCgEgBg
(,;1.TheInvestmentCo.Ltd.oftheStateRailwaororationofChinaBeiin00160,China2.KeaboratorfRoadandyCpjg1yLyo,,RailwanineerinafetontroloftheEducationMinistrShiiazhuaniedaoUniversitShiiazhuan50043,China)yEggSyCyjgTyjg0
122
,CHENHonbo WANGNin WANGJianxigg,
theshieldGtunnelinetroconstructionunderGcrossinxistinuildinsintheconditionsofknowinocertaindata.Uonthegmgegbggnp,thequalitativeassessmenttotheuantitativeassessmentisputforwardinthepresentaerwiththeanaltichierarchrocessqppyypandthestatisticalmethodtoobtaintheconditionalrobabilittable(CPT)whichreresenttheeffectsofthenodes.WithpracGpypboththeriskpredictionfortheconstructionandthesensitivitfthefactorsareanalzed.Theresultsofthemodelanalsisareyoyy(AHP)andthestateGaffectinfactorsintroducedtoobtainthediitalcharacteristicsofthecloud.Usinherandomriskcloudgggt,basisofcombininhetrianularfuzzumberwiththecloudtheortheinversionmethodoftheprobabilitistributionfromgtgynyyd
:AtAbstractrianularfuzzloudBaesiannetworkGbased(BN)riskanalsismodelisestablishedfortheriskassessmentofgycyy
,,ticalroectsasthebackroundtheriskevaluationindexsstemandtheBNmodelareestablisheduonthebasisofwhichpjgyp,foundtobeinconsistencewiththoseoftheactualproectwhichverifiestherationalitndeffectivenessofthemodel.jya
:;;;;KeordsshieldtunnelinnderGcrossinxistinuildinsriskassessmenttrianularfuzzumbercloudmodelBaesiannetgugegbggynyyw
takeintoconsiderationtheeffectofeccentricloadswheneccentricloadsareexerted.ThesoftwareoftheMidasFEAisusedin
,,,,,alsesinwhichcasehoweverthebeamelementmodelwhenusinuchsoftwarehastousetheeccentricloadcoefficienttoygs
(,,,)下接第26页)analsissoftware(suchastheMidasCivilBSASandDr.Brideetc.isoftenusedforcalculationsandanGyg
thepresentpaertohelstablishasolidelementmodelforanentiregirderofthesinleGboxtwinGcellturnoutcontinuousppeg,uidethedesin.Throuhtheanalsesitisfoundthattheeccentricloadeffectofliveloadsneedstobeconsideredwhenthegggy
,deflectioncheckinndthechoiceofthebearinaacitfsuortsareperformed.Inactualdesiniftherearenosimilargagcpyoppgcurateeccentricloadcoefficient.
,bridedesinexamlestobefolloweditisrecommendedthattheplateshellorsolidelementmodelsbeusedtoobtaintheacGggp:;;;KeordsasinleboxwithtwinGcellsturnoutcontinuousbeambridestheeffectofeccentricloadssolidelementmodelsggyw
()国防交通工程与技术 1020年5月 18,038 2
beam,uonthebasisofwhichanaccurateeccentricloadcoefficientundertheactionoftheeccentricZKloadisobtainedtop
因篇幅问题不能全部显示,请点此查看更多更全内容
Copyright © 2019- 91gzw.com 版权所有 湘ICP备2023023988号-2
违法及侵权请联系:TEL:199 18 7713 E-MAIL:2724546146@qq.com
本站由北京市万商天勤律师事务所王兴未律师提供法律服务