ISSN 1003-8035 CN 11-2852/P
  • 中国科技核心期刊
  • CSCD收录期刊
  • Caj-cd规范获奖期刊
  • Scopus 收录期刊
  • DOAJ 收录期刊
  • GeoRef收录期刊
欢迎扫码关注“i环境微平台”

基于节理不确定性的可靠度分析以西藏某岩质边坡为例

胡康, 任光明, 常文娟, 李征征, 邹林志

胡康,任光明,常文娟,等. 基于节理不确定性的可靠度分析−以西藏某岩质边坡为例[J]. 中国地质灾害与防治学报,2022,33(2): 53-60. DOI: 10.16031/j.cnki.issn.1003-8035.2022.02-07
引用本文: 胡康,任光明,常文娟,等. 基于节理不确定性的可靠度分析−以西藏某岩质边坡为例[J]. 中国地质灾害与防治学报,2022,33(2): 53-60. DOI: 10.16031/j.cnki.issn.1003-8035.2022.02-07
HU Kang, REN Guangming, CHANG Wenjuan, et al. Reliability analysis based on joint uncertainty: A case study of a rock slope in Tibet[J]. The Chinese Journal of Geological Hazard and Control, 2022, 33(2): 53-60. DOI: 10.16031/j.cnki.issn.1003-8035.2022.02-07
Citation: HU Kang, REN Guangming, CHANG Wenjuan, et al. Reliability analysis based on joint uncertainty: A case study of a rock slope in Tibet[J]. The Chinese Journal of Geological Hazard and Control, 2022, 33(2): 53-60. DOI: 10.16031/j.cnki.issn.1003-8035.2022.02-07

基于节理不确定性的可靠度分析——以西藏某岩质边坡为例

详细信息
    作者简介:

    胡 康(1996-),男,四川内江人,硕士研究生,环境地质专业。E-mail: 898863062@qq.com

    通讯作者:

    任光明(1964-),男,四川南充人,教授,研究领域为岩土体工程特性及稳定性与环境效应、地质灾害评价等。E-mail: 635863163@qq.com

  • 中图分类号: TU457

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.
  • 我国是世界上地震活动最强烈和地震灾害最严重的国家之一。从全球地震带的划分来看,我国位于太平洋地震带和欧亚地震带之间,地质构造复杂。破坏性地震灾害对我国人民的生命财产安全和生产生活造成极大的威胁[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)构建了城市地震易损性综合评价指标体系。利用熵权法确定指标的权重,将云理论与雷达图结合起来,构建城市地震易损性评价模型,解决城市地震综合易损性评价过程中的模糊性与随机性问题,并运用雷达图法解决了城市区域同一等级难以进行优劣评定问题,为城市合理配置防灾资源提供依据,提高城市的综合抗震防灾能力。

    熵是一个热力学概念,后由申农引入信息论,称之为信息熵。按照信息论基本原理的解释,熵是系统无序程度的一个度量,而信息熵是系统有序程度的度量。一个指标的信息熵越小,则该指标的变异程度越大,在综合评价中所起的作用就越大,指标的权重越大[7]。熵权法的基本思路是根据各变量的变异程度来确定各指标的权重,通过修正得到权重,所以运用熵权法可以保证结果的客观性和准确性。假设有m个待评价项目,n个评价指标,形成原始判断矩阵R=(rij)m×n,通过式(1)与式(2)分别计算第j个指标的信息熵值与权重。

    {pij=riji=1mrijej=1lnmi=1mpijlnpij (1)
    ωj=1ejj=1n(1ej) (2)

    式中:pij——第i个项目的第j个评价指标的权重值;

    rij——第j个指标下第i个项目的评价值;

    m——待评价项目个数;

    ej——第j个指标的信息熵值;

    ωj——第j个指标的权重。

    (1)云模型的定义

    X是一个普通集合,X={x},称为论域,对于论域中的任一元素x,都存在一个稳定倾向的随机数A(x)[0,1]与之相对应,A(x)叫做xA的隶属度。如果论域中的元素是简单有序的,则X为基础变量,如若不是,根据某个法则f,将X映射到一个有序的论域X,则X为基础变量,隶属度在基础变量上的分布称为云[8]

    (2)云模型的数字特征

    云模型通过使用三个数字特征:期望Ex、熵En、超熵He来表示语言值的数学性质。其中期望是最能代表定性概念的点,熵用来反映定性概念概率和模糊度,熵越大,概念所接受的范围也越大,概念就越模糊,超熵是熵的熵,用来反映云的离散程度及厚度。

    目前,虽然有很多学者对地震易损性进行了定义,但仍未形成统一的定义。尹之潜[9]认为,地震受体的易损性是指受体在确定强度地震的作用下,设定的极限状态的发生概率,受体可以是一座城市,也可以是城市里的建筑群或单体建筑。叶清[10]指出,地震易损性是指一个确定区域内由于地震发生造成损失的程度,易损性是评定震害的一个指数,对地震预测区内未来地震造成建筑物破坏和损失程度作出的预测。综合上述观点,本文研究的城市承灾体综合易损性是指:假定该城市在地震灾害的影响下,由涉及的自然、人、生命线系统、经济状况等所组成的承灾体系统在面临地震灾害时受到灾害破坏或损害程度的性质。该程度的大小主要由承灾体面临的震害压力、承灾体状态和响应程度决定。

    PSR模型由加拿大统计学家提出,20世纪80年代到90年代,经济合作与发展组织(OECD)和联合国环境规划署将PSR发展起来用以研究生态环境问题[11]。该理论认为,人类的经济、社会活动与自然环境之间存在相互作用关系,即:人类从自然环境取得各种资源,通过生产消费又向环境排放,从而改变了资源的数量和环境的质量,进而影响了人类的经济社会活动及其福利,如此循环往复,形成了人类活动与自然环境之间的PSR关系[11]。PSR目前已被广泛应用于生态安全评价、生态系统健康状况评价。由于城市地震易损性高低涉及城市的震害压力、生命线系统等状态和响应程度的大小,所以利用PSR可以很好地表达各要素间的关系。

    根据城市地震灾害的特点及易损性的概念,查询《城市抗震防灾规划管理规定》和《城市抗震防灾规划编制标准》等相关规范,综合考虑评价指标数据的可获取性,结合现有的文献[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
    社会状态C414~65岁人口比例C4114~65岁人口比例/%
    失业率C42失业人数占人口总数比例/%+
    教育C43人均教育费用支出/元
    医疗C44人均医疗卫生费用支出/元+
    经济状态C5人均GDP C51GDP/%
    人均GDP增长率C52人均GDP增长率/%
    可支配收入C53人均可支配收入/元
    经济多样性C54第三产业构成比例/%
    响应(B3)工程抗震C6建筑物抗震能力C612002年以后的建筑物占总建筑物比例/%
    构筑物抗震能力C622013年以后的构筑物占总建筑物比例/%
    回复能力C7医疗救助能力C71病床数和医生数/10万人
    社会保障C72人均抚恤和社会福利救济费/元
    地方财政收入C73年末财政收入/万元
    应急避难所覆盖度C74学校数和广场数/个
    避震减灾知识普及率C75
    政府应急反应能力C76
      注:指标属性中,“+”代表越大越优(正向)指标,“−”代表越小越优(逆向)指标。
    下载: 导出CSV 
    | 显示表格

    城市地震综合易损性是一个错综复杂的庞大体系,传统的评价方法难以衡量易损性评价过程中的模糊性与随机性,因此,本文引入了云模型,以解决城市地震综合易损性评价中的模糊性与随机性问题,并运用熵权法赋权,确定城市地震综合易损性等级。云模型可以对某一区域的综合水平进行评价,但对于同一等级的研究对象,无法对其相对优劣进行判定,且城市各个区综合地震易损性差异较小,为了明确各区的地震综合易损性程度,本文采用雷达图法确定各个区地震综合易损性相对高低。

    根据易损性的定义,本文引入云模型来量化城市地震综合易损性的大小。应用逆向云发生器分别生成易损性评价云的数字特征,具体公式如式(3)所示。

    {Ex=1qi=1qxiEn=π2×1qi=1q|xiEx|He=|S2En2|S2=1q1i=1q(xiEx)2 (3)

    式中:Ex——期望值;

    q——研究城市区域的数量;

    xi——第i个区域的易损性指数值;

    En——熵;

    He——超熵;

    S2——城市易损性指数值方差。

    将城市地震易损性大小划分为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)
    下载: 导出CSV 
    | 显示表格
    Ex=Sjmax+Sjmin2 (4)
    En=SjmaxSjmin22ln2 (5)
    He=k (6)

    式中:ExEnHe−分别为标准云的期望值、熵、           超熵;

    Sjmin——第j个区间的下边界值;

    Sjmax——第j个区间的上边界值;

    k−常数,根据变量本身的模糊度确定[15],本文取0.01。

    根据式(3)算出二级云特征数字矩阵,然后与二级权重矩阵进行合成运算,将二级权重根据二级指标的个数进行归一化处理,进一步得到一级云特征数字矩阵,然后将得出的一级云特征数字矩阵与一级权重矩阵进行合成运算,最后得到综合云数字特征。矩阵的合成运算方法如式(7)所示。

    (ω1,ω2,ω3,ωn)(Ex1En1He1ExnEnnHen)=(Ex,En,He) (7)

    式中:ExEnHe−分别为更高等级的期望值、           熵、超熵;

    n——评价指标的个数。

    利用Matlab对正向云发生器编程生成云图,对综合云与标准云进行分析对比,初步得到评价对象整体的等级。

    考虑到评价结果存在相似性,仅仅从综合评价云图上无法准确判别综合易损性等级,提出一种综合易损性等级与标准云贴近程度的计算方法。贴近度越大,城市综合易损性等级越靠近该评价云所对应的标准云等级,计算公式如式(8)所示。

    T=1|ExEx| (8)

    式中:T——综合易损性等级的贴近度。

    (1)评价指标无量纲化处理

    对各评价指标的基础数值进行标准化处理就是使得标准化处理后的值在[0,1]范围内。指标处理如式(9)所示[16]

    rij={rijmaxrij正向指标处理1rijmaxrij逆向指标处理 (9)

    式中:rij——第j个指标下第i个项目的评价值;

    rij——第j个指标下第i个项目的评价转化值。

    maxrij——第j个指标下所有项目的最大值。

    (2)评价指标权重与角度转换

    评价指标中有m个指标,则在雷达图中绘制m条指标轴。通过式(10)将权重转化为指标轴之间的角度。

    θi=360ωi (10)

    式中:θi——夹角角度;

    ωi——权重。

    (3)雷达图绘制

    通过Excel,对评价指标基础数值进行无量纲化处理,并投影在对应的评价指标数轴上,然后将投影点连接成线,绘制雷达图。

    (4)特征量计算

    采取雷达图的平均面积和周长两个特征量,构造一个二维特征量来计算综合评价值的函数,并根据评价函数计算各评价对象的评价值。

    平均面积:Sx=mi=1n1i<jn12RiRjsinθijCn2 (11)
    周长:Lx=mi=1n1i<jnRi2+Rj22RiRjcosθijCn2 (12)

    式中:x——拟评区域;

    Sx——雷达图平均面积;

    Lx——雷达图周长;

    m——评价指标数;

    R——指标轴的边长;

    i,j——其中的第i,j个指标;

    θ——两条边的夹角。

    Sx为面积评价向量,则周长评价向量Vx为:

    Vx=Sxπ(Lx/2π)2 (13)

    最后,各区域的综合评价值Yx为:

    Yx=Sx×Vx (14)

    甘肃省兰州市是西北地区重要的工业基地和交通枢纽,西部地区重要的中心城市之一,全市土地面积13085.6 km2,市内地貌类别多样。兰州处于祁连山地震带,历史上发生多次重大地震,造成极大损失。兰州市下辖5区、3县和兰州新区。为了数据采集方便,本文以5区为研究对象建立综合评价指标体系,评定兰州市辖区的地震易损性等级,将各个区的易损性状况进行排序,从中也获得了其易损性高或低的原因。

    本文选取兰州市5个中心城区(城关区、七里河区、西固区、安宁区、红古区)作为研究对象。邀请13位专家依据城市地震灾害的特点和兰州地域特色,针对前文建立的评价指标体系各项指标的重要程度进行打分,指标越重要,相对分值越高。根据式(1)与式(2)算得各个指标的权重,结果如下。

    W{C11C76}={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}

    关于数据的采集,本文采用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
    指标C11C21C22C23C24C25C26C31C32C33
    城关区0.20099.1200.0306332.00013.1206.2304.800624.660239568.750102878.160
    七里河区0.20083.6100.3601458.00015.56015.3003.8001153.230117122.50045088.580
    西固区0.20087.2800.1101023.00017.89017.7503.6701105.180101151.25029424.000
    安宁区0.200100.0000.0003443.00020.86019.0103.190288.31042590.00023285.180
    红古区0.15078.1700.400263.00027.12026.1003.1001633.74031942.50011007.540
    指标C34C35C41C42C43C44C51C52C53C54
    城关区8950.3592041.20082.4001.9101188.0001009.00072396.0006.90036326.00086.300
    七里河区3922.685894.60082.1001.9601324.0001124.00079074.0006.00028534.00062.600
    西固区2559.875583.80081.3002.0101751.0001487.000110782.0004.30033150.00038.200
    安宁区2025.805462.00083.6001.830935.000794.00059814.0006.60032574.00061.500
    红古区957.650218.40085.6003.1501488.0001264.00092080.0008.60025720.00031.000
    指标C61C62C71C72C73C74C75C76
    城关区82.05028.9501270.0557709.000400518.000438.00085.36076.470
    七里河区55.20017.0301247.3056141.000210098.000231.00071.05066.320
    西固区48.97011.750566.3956338.000124664.000144.00068.36063.100
    安宁区36.44011.250362.8955370.000117027.00096.00066.07060.560
    红古区32.3207.260854.6754657.00052393.00054.00056.32041.180
    下载: 导出CSV 
    | 显示表格
    表  4  兰州市5个中心城区统计数据的无量纲化处理
    Table  4.  Dimensionless of statistical data of five central urban districts of Lanzhou
    指标C11C21C22C23C24C25C26C31C32C33
    城关区1.0000.9910.0751.0000.4840.2391.0000.6180.0000.000
    七里河区1.0000.8360.9000.2300.5740.5860.7920.2940.5110.562
    西固区1.0000.8730.2750.1620.6600.6800.7650.3240.5780.714
    安宁区1.0001.0000.0000.5440.7690.7280.6650.8240.8220.774
    红古区0.7500.7821.0000.0421.0001.0000.6460.0000.8670.893
    指标C34C35C41C42C43C44C51C52C53C54
    城关区0.0000.0000.0370.6060.3220.6790.3470.1980.0000.000
    七里河区0.5620.5620.0410.6220.2440.7560.2860.3020.2150.275
    西固区0.7140.7140.0500.6380.0001.0000.0000.5000.0870.557
    安宁区0.7740.7740.0230.5810.4660.5340.4600.2330.1030.287
    红古区0.8930.8930.0001.0000.1500.8500.1690.0000.2920.641
    指标C61C62C71C72C73C74C75C76
    城关区0.0000.0000.0000.0000.0000.0000.0000.000
    七里河区0.3270.4120.0180.2030.4750.4730.1680.133
    西固区0.4030.5940.5540.1780.6890.6710.1990.175
    安宁区0.5560.6110.7140.3030.7080.7810.2260.208
    红古区0.6060.7490.3270.3960.8690.8770.3400.461
    下载: 导出CSV 
    | 显示表格

    根据式(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)C10.030(0.950,0.100,0.050)C110.030(0.950,0.100,0.050)
    C20.155(0.626,0.278,0.103)C210.021(0.896,0.099,0.025)
    C220.039(0.450,0.501,0.178)
    C230.018(0.395,0.377,0.079)
    C240.038(0.697,0.188,0.067)
    C250.025(0.647,0.235,0.143)
    C260.014(0.773,0.123,0.070)
    C30.228(0.552,0.307,0.150)C310.039(0.412,0.310,0.070)
    C320.045(0.556,0.301,0.171)
    C330.050(0.588,0.308,0.165)
    C340.046(0.588,0.308,0.165)
    C350.048(0.588,0.308,0.165)
    C40.124(0.491,0.086,0.042)C410.018(0.030,0.019,0.006)
    C420.050(0.690,0.156,0.080)
    C430.038(0.236,0.162,0.069)
    C440.018(0.764,0.162,0.069)
    C50.148(0.508,0.113,0.032)C510.051(0.252,0.168,0.051)
    C520.018(0.247,0.155,0.093)
    C530.031(0.139,0.114,0.009)
    C540.048(0.352,0.248,0.060)
    C60.074(0.423,0.240,0.108)C610.039(0.378,0.215,0.105)
    C620.035(0.473,0.268,0.112)
    C70.241(0.562,0.201,0.077)C710.044(0.323,0.314,0.046)
    C720.021(0.216,0.134,0.064)
    C730.055(0.548,0.311,0.129)
    C740.028(0.560,0.325,0.123)
    C750.049(0.187,0.103,0.067)
    C760.044(0.195,0.140,0.094)
    下载: 导出CSV 
    | 显示表格
    图  1  兰州中心城区综合易损性云图
    Figure  1.  Comprehensive vulnerability cloud map of center districts in Lanzhou

    为了明确各区的易损性程度,采用雷达图法来确定各个区综合易损性相对高低。根据式(9)、式(10)、表4,进行标准化处理与角度转换,绘制雷达图(图2)。根据式(11)~式(14),分别计算出各个区域雷达图的特征量(表6)。通过表6,兰州中心城市各区域地震综合易损性高低排序为红古区>安宁区>西固区>七里河区>城关区。对比5个区域,西固、安宁与红古的老旧建筑物比重大,导致易损性压力较大;红古区的易损性最大,其生命线、社会、经济易损性皆处于较高水平,其地方状态较差;红古、安宁的工程抗震防范措施差,地方财政收入低、医疗卫生条件差、避灾减灾意识差使地方的易损性响应差。将各个区的易损性状况进行排序,从中也获得了其易损性高或低的原因,经济发展水平较高的区域其综合易损性较低,经济欠发达的区域其综合易损性较高。究其原因,地方政府的财政收入使得该区域有更为雄厚的社会财富和生产力资源,财政收入较高的区域地方政府有更多的财力与精力用于防灾减灾建设。城市地震综合易损性评价结果体现了城市系统在面临地震的压力、状态与响应状况,易损性越高,说明城市在地震灾害下受到危害的可能性和损失程度越大。通过进一步优化人口结构、加强城市建筑抗震规范实施、合理规划城市布局、完善生命线系统的建设、提高城市的医疗卫生条件,加大力度宣传防灾减灾知识等,重点提高城市的经济发展水平,从而降低城市的地震易损性。

    图  2  兰州中心城市各区域综合易损性雷达图
    Figure  2.  Radar map of comprehensive vulnerability of various regions in central districts of Lanzhou
    表  6  兰州中心城市各区域的特征量计算结果
    Table  6.  The calculation results of the characteristic quantity of various regions in central districts of Lanzhou
    城市区平均面积/km2平均周长/kmVxYx排序
    城关区0.37017.44860.08380.17615
    七里河区0.61798.52950.10670.25644
    西固区0.786210.68160.08650.26083
    安宁区0.93208.73030.15360.37812
    红古区1.223411.43040.11760.37911
    下载: 导出CSV 
    | 显示表格

    (1)以城市承灾系统为研究对象,从压力、状态、响应3个方面出发,建立了城市综合易损性评价指标体系,包含一级指标7个,二级指标28个。本文运用云模型结合雷达图分析法所建立的评价模型不仅可以解决城市易损性综合评价过程中的模糊性与随机性问题,也可以对城市区域易损性的相对高低进行对比,容易找出影响城市综合易损性高低的主要因素。

    (2)基于熵权-云模型对兰州中心城区进行易损性综合评价,结果表明:兰州中心城区的地震综合易损性等级偏向Ⅲ级,易损状况中等,勉强满足城市承受震害能力的要求,但整体水平还有待提高;运用雷达图分析法将研究区域进行易损性高低排序,结果表明经济水平较高的城关、七里河区的易损性程度较低,通过重点提高城市的经济发展水平,加强城市建筑抗震规范实施、合理规划城市布局、完善生命线系统的建设等,从而降低城市的地震易损性。

  • 图  1   Baecher节理网络模型

    Figure  1.   Baecher joint network model

    图  2   边坡出露典型节理

    Figure  2.   Exposed typical joints of slope

    图  3   节理极点等密度图

    Figure  3.   Isodensity map of joint poles

    图  4   边坡模拟剖面结构图

    Figure  4.   Simulation model

    图  5   边坡节理有限元模型

    Figure  5.   Finite element model of rock slope joints

    图  6   PEM法计算稳定性系数概率分布图

    Figure  6.   Probability distribution of stability coefficient calculated by PEM method

    图  7   最大剪应变云图(添加节理)

    Figure  7.   Nephogram of maximum shear strain(adding joints)

    图  8   最大剪应变云图(不添加节理)

    Figure  8.   Nephogram of maximum shear strain (without joints)

    图  9   边坡位移矢量图(添加节理)

    Figure  9.   Slope displacement vector diagram (add joints)

    图  10   边坡位移矢量图(不添加节理)

    Figure  10.   Slope displacement vector diagram (without joints)

    图  11   边坡沿节理剪切变形

    Figure  11.   Shear deformation of slope along joints

    图  12   屈服节理分布图(局部)

    Figure  12.   Distribution of yield joints (local)

    表  1   岩体参数取值

    Table  1   Parameter values of rock mass

    材料弹性模量/MPa泊松比φ/(°)c/MPa容重/(kN·m−3)
    均值标准差均值标准差
    强风化花岗岩100000.283820.120.0123.5
    弱风化花岗岩130000.264220.150.0225.0
    微新花岗岩200000.254930.500.0427.0
    下载: 导出CSV

    表  2   结构面几何特征值

    Table  2   Geometric characteristc values of structural plane

    边坡分层节理组产状迹长密度/(条·m−2
    平均倾向/(°)平均倾角/(°)分布类型Fisher常数均值/m分布类型
    强风化层J123567Fisher958负指数分布0.20
    J212545455.50.25
    J1245745015负指数分布0.06
    弱风化层J212558Fisher4050.03
    J3105309050.05
    微新基岩J123678Fisher306负指数分布0.02
    J21355511040.01
    下载: 导出CSV

    表  3   结构面强度特征值

    Table  3   Strength characteristic values of structural plane

    节理组号φ/(°)c/MPa法向刚度/(MPa·m−1切向刚度/(MPa·m−1抗拉强度/MPa
    均值标准差均值标准差
    J13020.040.002800025000
    J23320.080.0031000030000
    J33830.100.0051300038000
    下载: 导出CSV

    表  4   敏感性计算结果

    Table  4   Sensitivity calculation results

    节理强度参数岩体强度参数稳定性系数
    c/MPφ/(°)c/MPaφ/(°)
    0.04300.15421.35
    0.034300.15421.33
    0.046300.15421.38
    0.04240.15421.23
    0.04360.15421.45
    0.04300.144421.34
    0.04300.156421.37
    0.04300.150361.27
    0.04300.150481.44
    下载: 导出CSV

    表  5   敏感性指标计算结果

    Table  5   Sensitivity index calculation results

    节理强度参数岩体强度参数
    c/MPφ/(°)c/MPφ/(°)
    0.03700.16300.02220.1259
    下载: 导出CSV

    表  6   可靠度对比分析结果

    Table  6   Results of reliability comparison analysis

    稳定性参数JEF-PEM点估计MCS蒙托卡罗
    平均值1.3441.590
    标准差0.05730.0854
    可靠性指数6.0076.907
    破坏概率0.001<0.001
    下载: 导出CSV
  • [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)

图(12)  /  表(6)
计量
  • 文章访问数:  426
  • HTML全文浏览量:  192
  • PDF下载量:  308
  • 被引次数: 2
出版历程
  • 收稿日期:  2021-06-02
  • 修回日期:  2021-07-29
  • 网络出版日期:  2022-03-31
  • 刊出日期:  2022-04-26

目录

/

返回文章
返回