Reliability analysis based on joint uncertainty: A case study of a rock slope in Tibet
-
摘要: 边坡稳定性一直是边坡安全的重点研究对象,针对边坡评价中常见的不确定性因素,可靠度分析是值得利用的方法。为评价某节理发育的岩质岸坡稳定性,通过有限元计算软件,结合现场勘探测绘数据,建立以边坡节理强度参数c、φ为输入变量,安全系数为输出变量的点估计(PEM)计算概率模型,计算结果表明:节理发育对该边坡变形具有明显控制作用;边坡整体可靠性较好,破坏概率极低。最后,通过蒙托卡罗法对可靠度结果进行验证,结果表明两种方法的计算结果不存在显著性差异。研究结果表明节理对岩质边坡稳定具有良好的敏感性,基于节理不确定性的点估计法分析边坡可靠度是一种有效的方法。Abstract: Slope stability has always been the key research object of slope safety. In view of the common uncertain factors in slope evaluation, reliability analysis is a valuable method. To evaluate the stability of a jointed rock slope, a point estimation (PEM) calculation probability model is established by using the finite element analysis software and the field survey data, The PEM method takes the joint strength parameters c, φ and the safety factor as input and output variables. The results suggest that, the joints have a significant control effect on the deformation of the slope; The overall reliability of slope is great, and the failure probability is extremely low. Finally, the reliability results are verified by Monte Carlo method, and the results show that there is no significant difference between the two methods. The results indicate that the joints have good sensitivity to the stability of rock slope, and the point estimation method based on joint uncertainty is an effective method to analyze the slope reliability.
-
Keywords:
- Jointed rock slope /
- reliability /
- baecher model /
- point estimate /
- Monte carlo
-
0. 引言
我国是世界上地震活动最强烈和地震灾害最严重的国家之一。从全球地震带的划分来看,我国位于太平洋地震带和欧亚地震带之间,地质构造复杂。破坏性地震灾害对我国人民的生命财产安全和生产生活造成极大的威胁[1]。近年来我国地震活动频繁,给受灾地区的人们带来了沉痛的灾难,严重地阻碍和制约了我国的经济发展。尤其是2008年5月的汶川8.0级大地震,致使87406人遇难或失踪,直接经济损失达8451亿元;2014年4月的玉树7.1级地震,致使2968人遇难或失踪,造成实际经济损失超过125亿元。由此可见,研究抗震防灾具有明显的必要性和紧迫性。伴随着城市化进程的加快,在人口、建筑与财富聚集的同时也聚集了风险,一旦遭受地震,后果将不堪设想。对城市地震易损性研究是城市地震问题的重要研究内容,具有重要的科学实用价值与现实意义。
近年来,国内外学者越来越关注地震灾害的易损性研究。DUZGUN等[2]提出了集成的城市地震脆弱性评估框架,采用脆弱性指数全面考虑城市区域评估的结构易损性、社会经济易损性等。PASTICIER等[3]利用SAP2000对古建筑进行IDA分析,并利用分析结果对其地震易损性进行了评估。苏亮等[4]利用层次分析法,提出了一种基于模糊相似理论的地震易损性分析方法,通过类比分析实现目标结构的易损性评估,极大地提高了结构地震易损性分析的计算效率。臧石磊等[5]利用模糊综合评判法计算地震易损性矩阵,并对辽宁省的老旧砌体房屋进行震害预测;毕可为[6]将模糊数学中的熵权法确定权重和相似理论引入到群体易损性中,可快速准确评估出群体建筑物的震害情况。虽然国内外在城市地震易损性方面的研究取得了诸多成果,但大部分是针对某一类承灾体的地震易损性分析,而对城市承灾系统的研究却十分罕见,对于城市易损性评价指标体系的建立和评价方法的研究尚处于不成熟的阶段,评价体系还难以达到统一标准。
在既有研究的基础上,本文综合考虑了城市地震易损性评价中指标选取存在的问题,引入压力-状态-响应模型(Pressure State Response, PSR)构建了城市地震易损性综合评价指标体系。利用熵权法确定指标的权重,将云理论与雷达图结合起来,构建城市地震易损性评价模型,解决城市地震综合易损性评价过程中的模糊性与随机性问题,并运用雷达图法解决了城市区域同一等级难以进行优劣评定问题,为城市合理配置防灾资源提供依据,提高城市的综合抗震防灾能力。
1. 熵权、云模型的基本理论
1.1 熵权
熵是一个热力学概念,后由申农引入信息论,称之为信息熵。按照信息论基本原理的解释,熵是系统无序程度的一个度量,而信息熵是系统有序程度的度量。一个指标的信息熵越小,则该指标的变异程度越大,在综合评价中所起的作用就越大,指标的权重越大[7]。熵权法的基本思路是根据各变量的变异程度来确定各指标的权重,通过修正得到权重,所以运用熵权法可以保证结果的客观性和准确性。假设有
个待评价项目, 个评价指标,形成原始判断矩阵 ,通过式(1)与式(2)分别计算第 个指标的信息熵值与权重。(1) (2) 式中:
——第 个项目的第 个评价指标的权重值; ——第 个指标下第 个项目的评价值;m——待评价项目个数;
——第 个指标的信息熵值; ——第 个指标的权重。1.2 云模型的基本理论
(1)云模型的定义
设
是一个普通集合, ,称为论域,对于论域中的任一元素 ,都存在一个稳定倾向的随机数 与之相对应, 叫做 对 的隶属度。如果论域中的元素是简单有序的,则 为基础变量,如若不是,根据某个法则 ,将 映射到一个有序的论域 ,则 为基础变量,隶属度在基础变量上的分布称为云[8]。(2)云模型的数字特征
云模型通过使用三个数字特征:期望
、熵 、超熵 来表示语言值的数学性质。其中期望是最能代表定性概念的点,熵用来反映定性概念概率和模糊度,熵越大,概念所接受的范围也越大,概念就越模糊,超熵是熵的熵,用来反映云的离散程度及厚度。2. 城市地震综合易损性评价指标体系的建立
2.1 城市地震综合易损性的概念
目前,虽然有很多学者对地震易损性进行了定义,但仍未形成统一的定义。尹之潜[9]认为,地震受体的易损性是指受体在确定强度地震的作用下,设定的极限状态的发生概率,受体可以是一座城市,也可以是城市里的建筑群或单体建筑。叶清[10]指出,地震易损性是指一个确定区域内由于地震发生造成损失的程度,易损性是评定震害的一个指数,对地震预测区内未来地震造成建筑物破坏和损失程度作出的预测。综合上述观点,本文研究的城市承灾体综合易损性是指:假定该城市在地震灾害的影响下,由涉及的自然、人、生命线系统、经济状况等所组成的承灾体系统在面临地震灾害时受到灾害破坏或损害程度的性质。该程度的大小主要由承灾体面临的震害压力、承灾体状态和响应程度决定。
2.2 压力-状态-响应PSR模型构建
PSR模型由加拿大统计学家提出,20世纪80年代到90年代,经济合作与发展组织(OECD)和联合国环境规划署将PSR发展起来用以研究生态环境问题[11]。该理论认为,人类的经济、社会活动与自然环境之间存在相互作用关系,即:人类从自然环境取得各种资源,通过生产消费又向环境排放,从而改变了资源的数量和环境的质量,进而影响了人类的经济社会活动及其福利,如此循环往复,形成了人类活动与自然环境之间的PSR关系[11]。PSR目前已被广泛应用于生态安全评价、生态系统健康状况评价。由于城市地震易损性高低涉及城市的震害压力、生命线系统等状态和响应程度的大小,所以利用PSR可以很好地表达各要素间的关系。
2.3 城市地震综合易损性评价指标体系
根据城市地震灾害的特点及易损性的概念,查询《城市抗震防灾规划管理规定》和《城市抗震防灾规划编制标准》等相关规范,综合考虑评价指标数据的可获取性,结合现有的文献[12],遵循指标体系构建的基本原则:全面性、代表性、层次性与可操作性,基于PSR模型,遴选、总结和整合出28个变量构成评价指标体系(表1)。
表 1 城市地震综合易损性评价指标体系Table 1. Evaluation index system for comprehensive vulnerability of urban earthquake属性 一级指标 二级指标 指标计算说明 指标属性 压力(B1) 自然因素C1 城市近源地震等效等级C11 反映地震灾害危险度指数g + 人为因素C2 城市化率C21 城镇人口占人口总数比例/% + 年平均自然增长率C22 年平均自然增长率/% + 人口密度C23 人口密度/(人·km−2) + 建筑物密度C24 建筑物的基底面积占规划建设面积比例/% + 老旧建筑物比例C25 老旧建筑物数量占建筑物数量比例/% + 容积率C26 地上总建筑面积占可建设用地面积比例/% + 状态(B2) 生命线系统C3 交通系统C31 公路网综合能力/(104 m2) − 通讯系统C32 电信业务收入/万元 − 供电系统C33 供电量/万千瓦时 − 供水系统C34 供水总量、用水总量/(104 m3) − 供气系统C35 供气总量/m3 − 社会状态C4 14~65岁人口比例C41 14~65岁人口比例/% − 失业率C42 失业人数占人口总数比例/% + 教育C43 人均教育费用支出/元 − 医疗C44 人均医疗卫生费用支出/元 + 经济状态C5 人均GDP C51 GDP/% − 人均GDP增长率C52 人均GDP增长率/% − 可支配收入C53 人均可支配收入/元 − 经济多样性C54 第三产业构成比例/% − 响应(B3) 工程抗震C6 建筑物抗震能力C61 2002年以后的建筑物占总建筑物比例/% − 构筑物抗震能力C62 2013年以后的构筑物占总建筑物比例/% − 回复能力C7 医疗救助能力C71 病床数和医生数/10万人 − 社会保障C72 人均抚恤和社会福利救济费/元 − 地方财政收入C73 年末财政收入/万元 − 应急避难所覆盖度C74 学校数和广场数/个 − 避震减灾知识普及率C75 − − 政府应急反应能力C76 − − 注:指标属性中,“+”代表越大越优(正向)指标,“−”代表越小越优(逆向)指标。 3. 评价模型的构建
城市地震综合易损性是一个错综复杂的庞大体系,传统的评价方法难以衡量易损性评价过程中的模糊性与随机性,因此,本文引入了云模型,以解决城市地震综合易损性评价中的模糊性与随机性问题,并运用熵权法赋权,确定城市地震综合易损性等级。云模型可以对某一区域的综合水平进行评价,但对于同一等级的研究对象,无法对其相对优劣进行判定,且城市各个区综合地震易损性差异较小,为了明确各区的地震综合易损性程度,本文采用雷达图法确定各个区地震综合易损性相对高低。
3.1 确定易损性评价云
根据易损性的定义,本文引入云模型来量化城市地震综合易损性的大小。应用逆向云发生器分别生成易损性评价云的数字特征,具体公式如式(3)所示。
(3) 式中:
——期望值; ——研究城市区域的数量; ——第 个区域的易损性指数值; ——熵; ——超熵; ——城市易损性指数值方差。3.2 构建评判等级的标准云
将城市地震易损性大小划分为5个等级[13]。将区间[0,1]等分为5个子区间,其中第j个子区间表示为[Sjmin,Sjmax]。按照式(4)~式(6)确定标准云的数字特征(表2)[14]。
表 2 评价等级数字特征表Table 2. Digital characteristic of evaluation level评价等级 取值区间 综合易损性 数字特征 Ⅰ [0.0,0.2) 低 (0.1,0.085,0.01) Ⅱ [0.2,0.4) 较低 (0.3,0.085,0.01) Ⅲ [0.4,0.6) 中等 (0.5,0.085,0.01) Ⅳ [0.6,0.8) 较高 (0.7,0.085,0.01) Ⅴ [0.8,1.0] 高 (0.9,0.085,0.01) (4) (5) (6) 式中:
、 、 −分别为标准云的期望值、熵、 超熵; ——第 个区间的下边界值; ——第 个区间的上边界值; −常数,根据变量本身的模糊度确定[15],本文取0.01。3.3 综合云
根据式(3)算出二级云特征数字矩阵,然后与二级权重矩阵进行合成运算,将二级权重根据二级指标的个数进行归一化处理,进一步得到一级云特征数字矩阵,然后将得出的一级云特征数字矩阵与一级权重矩阵进行合成运算,最后得到综合云数字特征。矩阵的合成运算方法如式(7)所示。
(7) 式中:
、 、 −分别为更高等级的期望值、 熵、超熵;n——评价指标的个数。
3.4 综合评价云图
利用Matlab对正向云发生器编程生成云图,对综合云与标准云进行分析对比,初步得到评价对象整体的等级。
3.5 贴近度
考虑到评价结果存在相似性,仅仅从综合评价云图上无法准确判别综合易损性等级,提出一种综合易损性等级与标准云贴近程度的计算方法。贴近度越大,城市综合易损性等级越靠近该评价云所对应的标准云等级,计算公式如式(8)所示。
(8) 式中:
——综合易损性等级的贴近度。3.6 雷达图分析
(1)评价指标无量纲化处理
对各评价指标的基础数值进行标准化处理就是使得标准化处理后的值在[0,1]范围内。指标处理如式(9)所示[16]。
(9) 式中:
——第j个指标下第i个项目的评价值; ——第j个指标下第i个项目的评价转化值。 ——第j个指标下所有项目的最大值。(2)评价指标权重与角度转换
评价指标中有
个指标,则在雷达图中绘制 条指标轴。通过式(10)将权重转化为指标轴之间的角度。(10) 式中:
——夹角角度; ——权重。(3)雷达图绘制
通过Excel,对评价指标基础数值进行无量纲化处理,并投影在对应的评价指标数轴上,然后将投影点连接成线,绘制雷达图。
(4)特征量计算
采取雷达图的平均面积和周长两个特征量,构造一个二维特征量来计算综合评价值的函数,并根据评价函数计算各评价对象的评价值。
(11) (12) 式中:
——拟评区域; ——雷达图平均面积; ——雷达图周长; ——评价指标数; ——指标轴的边长; ——其中的第 个指标; ——两条边的夹角。令
为面积评价向量,则周长评价向量 为:(13) 最后,各区域的综合评价值
为:(14) 4. 案例分析
甘肃省兰州市是西北地区重要的工业基地和交通枢纽,西部地区重要的中心城市之一,全市土地面积13085.6 km2,市内地貌类别多样。兰州处于祁连山地震带,历史上发生多次重大地震,造成极大损失。兰州市下辖5区、3县和兰州新区。为了数据采集方便,本文以5区为研究对象建立综合评价指标体系,评定兰州市辖区的地震易损性等级,将各个区的易损性状况进行排序,从中也获得了其易损性高或低的原因。
4.1 研究对象及权重确定
本文选取兰州市5个中心城区(城关区、七里河区、西固区、安宁区、红古区)作为研究对象。邀请13位专家依据城市地震灾害的特点和兰州地域特色,针对前文建立的评价指标体系各项指标的重要程度进行打分,指标越重要,相对分值越高。根据式(1)与式(2)算得各个指标的权重,结果如下。
W{C11~C76}={0.030,0.021,0.039,0.018,0.038,0.025,0.014,0.039,0.045,0.050,0.046,0.048,0.018,0.050,0.038,0.018,0.051,0.018,0.031,0.048,0.039,0.035,0.044,0.021,0.055,0.028,0.049,0.044}
4.2 数据采集及处理
关于数据的采集,本文采用2018年《兰州市统计年鉴》和地方政府各部门公开报告中的统计数据。由于城市综合易损性的复杂性,针对表1中统计困难的指标(C75,C76),本文采用问卷调查法确定;由于C25, C61,C62难以量化[17],本文以建筑业占第二产业的投资额比例衡量,其余指标按照统计实际数据计算。表3是5个中心城区的基础统计数据,将表3按照式(9)进行无量纲化处理(表4)。
表 3 兰州市5个中心城区的基础数据统计Table 3. Basic data statistics of five central urban districtsin Lanzhou City指标 C11 C21 C22 C23 C24 C25 C26 C31 C32 C33 城关区 0.200 99.120 0.030 6332.000 13.120 6.230 4.800 624.660 239568.750 102878.160 七里河区 0.200 83.610 0.360 1458.000 15.560 15.300 3.800 1153.230 117122.500 45088.580 西固区 0.200 87.280 0.110 1023.000 17.890 17.750 3.670 1105.180 101151.250 29424.000 安宁区 0.200 100.000 0.000 3443.000 20.860 19.010 3.190 288.310 42590.000 23285.180 红古区 0.150 78.170 0.400 263.000 27.120 26.100 3.100 1633.740 31942.500 11007.540 指标 C34 C35 C41 C42 C43 C44 C51 C52 C53 C54 城关区 8950.359 2041.200 82.400 1.910 1188.000 1009.000 72396.000 6.900 36326.000 86.300 七里河区 3922.685 894.600 82.100 1.960 1324.000 1124.000 79074.000 6.000 28534.000 62.600 西固区 2559.875 583.800 81.300 2.010 1751.000 1487.000 110782.000 4.300 33150.000 38.200 安宁区 2025.805 462.000 83.600 1.830 935.000 794.000 59814.000 6.600 32574.000 61.500 红古区 957.650 218.400 85.600 3.150 1488.000 1264.000 92080.000 8.600 25720.000 31.000 指标 C61 C62 C71 C72 C73 C74 C75 C76 城关区 82.050 28.950 1270.055 7709.000 400518.000 438.000 85.360 76.470 七里河区 55.200 17.030 1247.305 6141.000 210098.000 231.000 71.050 66.320 西固区 48.970 11.750 566.395 6338.000 124664.000 144.000 68.360 63.100 安宁区 36.440 11.250 362.895 5370.000 117027.000 96.000 66.070 60.560 红古区 32.320 7.260 854.675 4657.000 52393.000 54.000 56.320 41.180 表 4 兰州市5个中心城区统计数据的无量纲化处理Table 4. Dimensionless of statistical data of five central urban districts of Lanzhou指标 C11 C21 C22 C23 C24 C25 C26 C31 C32 C33 城关区 1.000 0.991 0.075 1.000 0.484 0.239 1.000 0.618 0.000 0.000 七里河区 1.000 0.836 0.900 0.230 0.574 0.586 0.792 0.294 0.511 0.562 西固区 1.000 0.873 0.275 0.162 0.660 0.680 0.765 0.324 0.578 0.714 安宁区 1.000 1.000 0.000 0.544 0.769 0.728 0.665 0.824 0.822 0.774 红古区 0.750 0.782 1.000 0.042 1.000 1.000 0.646 0.000 0.867 0.893 指标 C34 C35 C41 C42 C43 C44 C51 C52 C53 C54 城关区 0.000 0.000 0.037 0.606 0.322 0.679 0.347 0.198 0.000 0.000 七里河区 0.562 0.562 0.041 0.622 0.244 0.756 0.286 0.302 0.215 0.275 西固区 0.714 0.714 0.050 0.638 0.000 1.000 0.000 0.500 0.087 0.557 安宁区 0.774 0.774 0.023 0.581 0.466 0.534 0.460 0.233 0.103 0.287 红古区 0.893 0.893 0.000 1.000 0.150 0.850 0.169 0.000 0.292 0.641 指标 C61 C62 C71 C72 C73 C74 C75 C76 城关区 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 七里河区 0.327 0.412 0.018 0.203 0.475 0.473 0.168 0.133 西固区 0.403 0.594 0.554 0.178 0.689 0.671 0.199 0.175 安宁区 0.556 0.611 0.714 0.303 0.708 0.781 0.226 0.208 红古区 0.606 0.749 0.327 0.396 0.869 0.877 0.340 0.461 4.3 结果分析
根据式(3)和表4,分别求得兰州中心城区二级指标的数字特征,运用式(7)将二级指标的数字特征与二级熵权矩阵进行合成运算,然后将一级指标的数字特征与一级熵权矩阵进行合成运算,得到兰州市中心城市区的易损性云综合数字特征(表5)。应用正向云发生器生成综合易损性与标准云对比云图(图1)。从图1知,兰州中心城区综合易损性在0.480处隶属度最高,云滴也最为集中,然而,综合云的跨度大于标准云,说明评价结果具有一定的模糊性,综合云的厚度大于标准云,说明评价结果存在一定的随机性。根据式(8)计算贴近度。结果表明:兰州市中心城区的综合易损性等级偏向Ⅲ级,易损状况中等,该评价结果说明兰州中心城区的综合易损性勉强满足城市承受震害能力的要求,但整体水平还有待提高,需要加强城市的管理与建设。
表 5 兰州中心城区地震综合易损性评价云数字特征Table 5. Cloud characteristic for comprehensive evaluation of earthquake vulnerability in center districts of Lanzhou总指标 数字特征 一级指标 权重 数字特征 二级指标 权重 数字特征 城市综合易损性评价 (0.554,0.210,0.088) C1 0.030 (0.950,0.100,0.050) C11 0.030 (0.950,0.100,0.050) C2 0.155 (0.626,0.278,0.103) C21 0.021 (0.896,0.099,0.025) C22 0.039 (0.450,0.501,0.178) C23 0.018 (0.395,0.377,0.079) C24 0.038 (0.697,0.188,0.067) C25 0.025 (0.647,0.235,0.143) C26 0.014 (0.773,0.123,0.070) C3 0.228 (0.552,0.307,0.150) C31 0.039 (0.412,0.310,0.070) C32 0.045 (0.556,0.301,0.171) C33 0.050 (0.588,0.308,0.165) C34 0.046 (0.588,0.308,0.165) C35 0.048 (0.588,0.308,0.165) C4 0.124 (0.491,0.086,0.042) C41 0.018 (0.030,0.019,0.006) C42 0.050 (0.690,0.156,0.080) C43 0.038 (0.236,0.162,0.069) C44 0.018 (0.764,0.162,0.069) C5 0.148 (0.508,0.113,0.032) C51 0.051 (0.252,0.168,0.051) C52 0.018 (0.247,0.155,0.093) C53 0.031 (0.139,0.114,0.009) C54 0.048 (0.352,0.248,0.060) C6 0.074 (0.423,0.240,0.108) C61 0.039 (0.378,0.215,0.105) C62 0.035 (0.473,0.268,0.112) C7 0.241 (0.562,0.201,0.077) C71 0.044 (0.323,0.314,0.046) C72 0.021 (0.216,0.134,0.064) C73 0.055 (0.548,0.311,0.129) C74 0.028 (0.560,0.325,0.123) C75 0.049 (0.187,0.103,0.067) C76 0.044 (0.195,0.140,0.094) 为了明确各区的易损性程度,采用雷达图法来确定各个区综合易损性相对高低。根据式(9)、式(10)、表4,进行标准化处理与角度转换,绘制雷达图(图2)。根据式(11)~式(14),分别计算出各个区域雷达图的特征量(表6)。通过表6,兰州中心城市各区域地震综合易损性高低排序为红古区>安宁区>西固区>七里河区>城关区。对比5个区域,西固、安宁与红古的老旧建筑物比重大,导致易损性压力较大;红古区的易损性最大,其生命线、社会、经济易损性皆处于较高水平,其地方状态较差;红古、安宁的工程抗震防范措施差,地方财政收入低、医疗卫生条件差、避灾减灾意识差使地方的易损性响应差。将各个区的易损性状况进行排序,从中也获得了其易损性高或低的原因,经济发展水平较高的区域其综合易损性较低,经济欠发达的区域其综合易损性较高。究其原因,地方政府的财政收入使得该区域有更为雄厚的社会财富和生产力资源,财政收入较高的区域地方政府有更多的财力与精力用于防灾减灾建设。城市地震综合易损性评价结果体现了城市系统在面临地震的压力、状态与响应状况,易损性越高,说明城市在地震灾害下受到危害的可能性和损失程度越大。通过进一步优化人口结构、加强城市建筑抗震规范实施、合理规划城市布局、完善生命线系统的建设、提高城市的医疗卫生条件,加大力度宣传防灾减灾知识等,重点提高城市的经济发展水平,从而降低城市的地震易损性。
表 6 兰州中心城市各区域的特征量计算结果Table 6. The calculation results of the characteristic quantity of various regions in central districts of Lanzhou城市区 平均面积/km2 平均周长/km 排序 城关区 0.3701 7.4486 0.0838 0.1761 5 七里河区 0.6179 8.5295 0.1067 0.2564 4 西固区 0.7862 10.6816 0.0865 0.2608 3 安宁区 0.9320 8.7303 0.1536 0.3781 2 红古区 1.2234 11.4304 0.1176 0.3791 1 5. 结论
(1)以城市承灾系统为研究对象,从压力、状态、响应3个方面出发,建立了城市综合易损性评价指标体系,包含一级指标7个,二级指标28个。本文运用云模型结合雷达图分析法所建立的评价模型不仅可以解决城市易损性综合评价过程中的模糊性与随机性问题,也可以对城市区域易损性的相对高低进行对比,容易找出影响城市综合易损性高低的主要因素。
(2)基于熵权-云模型对兰州中心城区进行易损性综合评价,结果表明:兰州中心城区的地震综合易损性等级偏向Ⅲ级,易损状况中等,勉强满足城市承受震害能力的要求,但整体水平还有待提高;运用雷达图分析法将研究区域进行易损性高低排序,结果表明经济水平较高的城关、七里河区的易损性程度较低,通过重点提高城市的经济发展水平,加强城市建筑抗震规范实施、合理规划城市布局、完善生命线系统的建设等,从而降低城市的地震易损性。
-
表 1 岩体参数取值
Table 1 Parameter values of rock mass
材料 弹性模量/MPa 泊松比 φ/(°) c/MPa 容重/(kN·m−3) 均值 标准差 均值 标准差 强风化花岗岩 10000 0.28 38 2 0.12 0.01 23.5 弱风化花岗岩 13000 0.26 42 2 0.15 0.02 25.0 微新花岗岩 20000 0.25 49 3 0.50 0.04 27.0 表 2 结构面几何特征值
Table 2 Geometric characteristc values of structural plane
边坡分层 节理组 产状 迹长 密度/(条·m−2) 平均倾向/(°) 平均倾角/(°) 分布类型 Fisher常数 均值/m 分布类型 强风化层 J1 235 67 Fisher 95 8 负指数分布 0.20 J2 125 45 45 5.5 0.25 J1 245 74 50 15 负指数分布 0.06 弱风化层 J2 125 58 Fisher 40 5 0.03 J3 105 30 90 5 0.05 微新基岩 J1 236 78 Fisher 30 6 负指数分布 0.02 J2 135 55 110 4 0.01 表 3 结构面强度特征值
Table 3 Strength characteristic values of structural plane
节理组号 φ/(°) c/MPa 法向刚度/(MPa·m−1) 切向刚度/(MPa·m−1) 抗拉强度/MPa 均值 标准差 均值 标准差 J1 30 2 0.04 0.002 8000 2500 0 J2 33 2 0.08 0.003 10000 3000 0 J3 38 3 0.10 0.005 13000 3800 0 表 4 敏感性计算结果
Table 4 Sensitivity calculation results
节理强度参数 岩体强度参数 稳定性系数 c/MP φ/(°) c/MPa φ/(°) 0.04 30 0.15 42 1.35 0.034 30 0.15 42 1.33 0.046 30 0.15 42 1.38 0.04 24 0.15 42 1.23 0.04 36 0.15 42 1.45 0.04 30 0.144 42 1.34 0.04 30 0.156 42 1.37 0.04 30 0.150 36 1.27 0.04 30 0.150 48 1.44 表 5 敏感性指标计算结果
Table 5 Sensitivity index calculation results
节理强度参数 岩体强度参数 c/MP φ/(°) c/MP φ/(°) 0.0370 0.1630 0.0222 0.1259 表 6 可靠度对比分析结果
Table 6 Results of reliability comparison analysis
稳定性参数 JEF-PEM点估计 MCS蒙托卡罗 平均值 1.344 1.590 标准差 0.0573 0.0854 可靠性指数 6.007 6.907 破坏概率 0.001 <0.001 -
[1] 陈祖煜, 陈立宏, 王玉杰, 等. 滑坡和建筑物抗滑稳定分析中的可靠度分析和分项系数设计方法[C]//陈祖煜. 水利水电工程风险分析及可靠度设计技术进展. 北京: 中国水利水电出版社, 2010: 27–39 CHEN Zuyu, CHEN Lihong, WANG Yujie, et al. Reliability analysis and partial factor design methods for slope stability[C]// CHEN Zuyu. Proceedings of the symposium on risk analysis and reliability based design for water resources and hydropower projects. Beijing: China Water Power Press, 2010: 27–39. (in Chinese)
[2] 刘超, 袁颖, 左朝晖, 等. 考虑土参数自相关距离影响的单桩可靠性设计方法[J]. 水文地质工程地质,2020,47(3):122 − 127. [LIU Chao, YUAN Ying, ZUO Zhaohui, et al. Pile reliability design considering the influence of soil auto-correlation distance[J]. Hydrogeology & Engineering Geology,2020,47(3):122 − 127. (in Chinese with English abstract) [3] 李涛, 刘国栋, 王聪. 基于可靠度理论的尾矿坝失稳概率及敏感性分析[J]. 中国地质灾害与防治学报,2019,30(3):81 − 86. [LI Tao, LIU Guodong, WANG Cong. Instability probability and sensitivity analysis of tailings dam based on reliability theory[J]. The Chinese Journal of Geological Hazard and Control,2019,30(3):81 − 86. (in Chinese with English abstract) [4] JIANG S H, LI D Q, ZHANG L M, et al. Slope reliability analysis considering spatially variable shear strength parameters using a non-intrusive stochastic finite element method[J]. Engineering Geology,2014,168:120 − 128. DOI: 10.1016/j.enggeo.2013.11.006
[5] LI Y J, HICKS M A, NUTTALL J D. Comparative analyses of slope reliability in 3D[J]. Engineering Geology,2015,196:12 − 23. DOI: 10.1016/j.enggeo.2015.06.012
[6] ZHAO L H, YU C H, CHENG X, et al. A method for seismic stability analysis of jointed rock slopes using Barton-Bandis failure criterion[J]. International Journal of Rock Mechanics and Mining Sciences,2020,136:104487. DOI: 10.1016/j.ijrmms.2020.104487
[7] LU R L, WEI W, SHANG K W, et al. Stability analysis of jointed rock slope by strength reduction technique considering ubiquitous joint model[J]. Advances in Civil Engineering,2020,2020:8862243.
[8] 冯开帅, 姜谙男, 吴洪涛, 等. 顺层非贯通节理边坡稳定性及破坏模式研究[J]. 公路工程,2020,45(6):52 − 58. [FENG Kaishuai, JIANG Annan, WU Hongtao, et al. Study on failure mode and stability of slope based on intermittent joints[J]. Highway Engineering,2020,45(6):52 − 58. (in Chinese with English abstract) [9] CHANG Y, CHANG L S, REN F Q. Strength anisotropy of jointed rock slope considering mining damage: A numerical analysis method[J]. Geomatics, Natural Hazards and Risk,2020,11(1):2587 − 2614. DOI: 10.1080/19475705.2020.1856200
[10] 陈国良, 廖国华. 岩体节理网络的计算机模拟[J]. 有色金属(矿山部分),1989,41(6):23 − 28. [CHEN Guoliang, LIAO Guohua. Computer simulation of rock mass joint network[J]. Nonferrous Metals (Mining),1989,41(6):23 − 28. (in Chinese) [11] 王宏, 陶振宇. 边坡稳定分析的节理网络模拟原理及工程应用[J]. 水利学报,1993,24(10):20 − 27. [WANG Hong, TAO Zhenyu. The principle of joint network simulation in the stability analysis of rock slope and its engineering application[J]. Journal of Hydraulic Engineering,1993,24(10):20 − 27. (in Chinese with English abstract) DOI: 10.3321/j.issn:0559-9350.1993.10.003 [12] 严豪, 宋彦辉, 陈子玉. 基于Voronoi节理模型的碎裂岩边坡稳定性分析[J]. 中国地质灾害与防治学报,2018,29(1):34 − 39. [YAN Hao, SONG Yanhui, CHEN Ziyu. Stability analysis of broken rock slope based on Voronoi joint net model[J]. The Chinese Journal of Geological Hazard and Control,2018,29(1):34 − 39. (in Chinese with English abstract) [13] 李源亮, 任光明, 刘艳领, 等. 基于节理网络有限元的岩质边坡稳定性分析[J]. 长江科学院院报,2019,36(1):78 − 83. [LI Yuanliang, REN Guangming, LIU Yanling, et al. Stability analysis of rock slope based on finite elements with joint network[J]. Journal of Yangtze River Scientific Research Institute,2019,36(1):78 − 83. (in Chinese with English abstract) DOI: 10.11988/ckyyb.20170823 [14] 肖术, 吴顺川, 高永涛, 等. 基于PEM-JFEM方法的节理岩质边坡稳定性评价[J]. 工程科学学报,2015,37(7):844 − 850. [XIAO Shu, WU Shunchuan, GAO Yongtao, et al. Jointed rock slope stability evaluation based on PEM- JFEM method[J]. Chinese Journal of Engineering,2015,37(7):844 − 850. (in Chinese with English abstract) [15] 张宜杰, 任光明, 常文娟, 等. 节理岩质边坡稳定性概率分析[J]. 成都理工大学学报(自然科学版),2021,48(2):235 − 241. [ZHANG Yijie, REN Guangming, CHANG Wenjuan, et al. Probabilistic analysis of stability of jointed rock slope[J]. Journal of Chengdu University of Technology (Science & Technology Edition),2021,48(2):235 − 241. (in Chinese with English abstract) [16] PRZEWLOCKI J, ZABUSKI L, WINKELMANN K. Reliability analysis of sea cliff slope stability by point estimate method[J]. IOP Conference Series:Materials Science and Engineering,2019,471:042003. DOI: 10.1088/1757-899X/471/4/042003
[17] 许湘华, 曲广琇, 方理刚. 基于节理几何参数不确定性的边坡可靠度分析[J]. 中南大学学报(自然科学版),2010,41(3):1139 − 1145. [XU Xianghua, QU Guangxiu, FANG Ligang. Reliability analysis of rock slope based on uncertainty of joint geometric parameters[J]. Journal of Central South University (Science and Technology),2010,41(3):1139 − 1145. (in Chinese with English abstract) [18] 王双, 陈征宙, 吴强, 等. 基于节理产状不确定性的边坡稳定性及敏感度分析[J]. 岩土工程学报,2013,35(2):348 − 354. [WANG Shuang, CHEN Zhengzhou, WU Qiang, et al. Stability and sensitivity analysis of slopes based on uncertainty of joint orientations[J]. Chinese Journal of Geotechnical Engineering,2013,35(2):348 − 354. (in Chinese with English abstract) [19] 王双. 节理产状概率模型研究及其在产状分组和岩坡不确定分析中的应用[D]. 南京: 南京大学, 2013 WANG Shuang. Research on distribution model of joint orientations and its application to joint set clustering and rock slope uncertainty analysis[D]. Nanjing: Nanjing University, 2013. (in Chinese with English abstract)
[20] 申矫健. 基于参数不确定性的岩质边坡稳定性分析方法研究[D]. 赣州: 江西理工大学, 2014 SHEN Jiaojian. Study on stability analysis method of rock slope based on parameter uncertainty[D]. Ganzhou: Jiangxi University of Science and Technology, 2014. (in Chinese with English abstract)
[21] BAECHER G B, LANNEY N A, EINSTEIN H H. Statistical description of rock properties and sampling[J]. 18th U S Symposium on Rock Mechanics, USRMS 1977,1977:5 − 8.
[22] 王宇, 李晓, 刘帅, 等. 岩体离散裂隙网络稳定性计算的节理有限元法[J]. 岩石力学与工程学报,2013,32(增刊 2):3337 − 3345. [WANG Yu, LI Xiao, LIU Shuai, et al. Stability calculation for discrete fracture network of rock masses based on joint finite element method[J]. Chinese Journal of Rock Mechanics and Engineering,2013,32(Sup 2):3337 − 3345. (in Chinese with English abstract) [23] CACAS M C, LEDOUX E, DE MARSILY G, et al. Modeling fracture flow with a stochastic discrete fracture network: calibration and validation: 1. The flow model[J]. Water Resources Research,1990,26(3):479 − 489.
[24] DERSHOWITZ W S, EINSTEIN H H. Characterizing rock joint geometry with joint system models[J]. Rock Mechanics and Rock Engineering,1988,21(1):21 − 51. DOI: 10.1007/BF01019674
[25] ROSENBLUETH E. Point estimates for probability moments[J]. Proceedings of the National Academy of Sciences of the United States of America,1975,72(10):3812 − 3814. DOI: 10.1073/pnas.72.10.3812
[26] 李侃, 巨能攀. 基于蒙特卡洛方法的边坡可靠性评价[J]. 中国地质灾害与防治学报,2014,25(1):23 − 27. [LI Kan, JU Nengpan. Integrated application of Monte-Carlo simulation for landslide reliability analysis[J]. The Chinese Journal of Geological Hazard and Control,2014,25(1):23 − 27. (in Chinese with English abstract) -
期刊类型引用(2)
1. 陈忠源,戴自航,简文彬. 基于因子权重反分析的新近失稳土质边坡稳定性评价云模型. 中国地质灾害与防治学报. 2023(04): 125-133 . 本站查看
2. 郄子君,郭文林. 常态化防控阶段公民新冠疫情应对行为及引导策略研究——基于PSR模型的实证分析. 中国卫生事业管理. 2022(12): 881-888+898 . 百度学术
其他类型引用(0)