201 8年第17卷第1 3期 考虑地震动参数不确定性的 边坡可靠度分析方法探索 口耿学勇 刘 哲 【内容摘要】影响边坡稳定性的因素包括岩土体强度、岩体结构、地下水特征、地震作用等,这些因素均存在不确定性。目前关 于边坡可靠度的研究主要集中在岩土体强度参数不确定性方面。本文初步探索了在考虑地震动加速度随机性特 征条件下的边坡可靠度分析方法,并以某核电厂边坡为例,开展可靠度分析计算。分析过程采用蒙特卡洛模拟法, 通过计算和分析可以发现地震动参数不确定性对边坡失效概率影响较大。 【关键词】可靠度分析;边坡稳定性;核电厂;地震动参数 【作者单位】耿学勇,刘哲;中国核电工程有限公司总图地质所 一、引言 场点收到的最大水平向地震动峰值加速度作为地震强度的 代表值。地震危险性分析概率方法是美国Cornell(1968)首 先引出的,其分析流程图见图1,_1 图2为某核电厂经过地震 危险陛评价后得到的水平向地震动峰值加速度年超越概率 曲线。需要指出的是,地震危险性分析概率方法基于以下几 个基本假定:一是潜在震源区内任何地方发生地震的可能性 是相同的;二是潜在震源区内地震的平均发生率在时间轴上 是个常数;三是地震年平均发生率符合泊松分布,即假定地 震事件是的,随机的(即地震发生的时间、震源坐标、震 级等以相互的方式出现)。四是一个地区内(潜在震源 区内)地震次数随震级增高以指数形式减少,大小震级之间 比例关系,可用古登堡一里克特震级频度关系表示。五是场 地地震动参数是震中距和震级的函数。 三、可靠度分析方法选取 随着理论和技术的不断发展,有关边坡可靠性的研究已 经成为岩土工程领域中一个重要的课题。边坡可靠度分析 就是用定量的方法分析评价影响边坡工程稳定性的不确定 性因素,其评价结果为边坡的失效概率或可靠度指标。 影响边坡稳定性的不确定性因素有很多,主要包括岩土 体强度、岩体结构、地下水特征、地震作用等。目前关于边坡 可靠度的研究主要集中在岩土体强度参数不确定性方面。 地震作用在核电厂边坡设计中具有非常大的影响权重,如果 不能考虑地震作用的不确定性,研究核电厂边坡的可靠度将 失去意义。 二、地震作用的随机特性 地震发生的时间、空间和强度都具有明显的随机性。核 电厂设计中,依据地震地质资料和历史地震资料,采用地震 危险性分析概率方法来建立地震发生的随机模型,并以工程 筋排布位置与拟绑扎梁板位置相同。植筋采用结构胶,不得 使用耐候密封胶等代替。由于梁柱节点等钢筋排布较密,部 分梁上植筋不具备可操作性,剔出框梁主筋,将新增框梁主 筋与之焊接。 第五,对隔油池进行了开挖,隔油池、排污沟重新砌筑并 可靠度分析过程中,将不确定性因素视为基本随机变量, 及附属构筑物;避免由于地基受损导致建筑不均匀沉降带来 的安全隐患甚至遭受不必要的经济损失。 第三,填充墙未做框架地梁,均坐落在素混凝土地梁之 上,承载力较差,当地基发生扰动时极易撕裂倾斜;加之厨房 用水量较多,当排污沟、隔油池等排水设施发生断裂或者渗 进行防水处理;与周边梁柱接合处进行混凝土灌浆处理防止 下沉。待地基地面处理完成后进行观察,无沉降情况继续发 生后砌筑填充墙体、抹灰以及装修处理工作。经处理后食堂 回复正常使用,后历次降雨后未发生室外积水浸泡现象,建 筑本体未再次出现沉降情况,现状良好。 六、结语 漏时,污水灌入回填土地基使得地基沉降情况加剧。对于地 基较差或者经回填处理的建筑可考虑使用框架地梁增加地 基的承载力,增加上部结构的稳定性。 【参考文献】 [1]GB5007—2002建筑地基基础技术规范[s].北京:中国建 筑工业出版社.2002 第一,#3食堂周围未设置雨排水沟,食堂周围地面未作 地面硬化;屋面雨排水直接冲刷散水渗入回填土地基是本次 事故的直接原因;因此,在设计之初应当考虑到雨排水的收 集排放,沿建筑物设计环形排水沟并进行地面硬化处理。 第二,后期建设冷库阻挡了地面雨排水的自然散流,使 得雨水在散水处积存是本次事故的次要原因。因此,在复杂 [2]JGJ79—2002建筑地基处理技术规范[s].北京:中国建 筑工业出版社。2002 [3]于兴银.地基不均匀沉降分析与加固改造措施[J].建筑 技术,2015,2:410~411 [4]曹瑞钠,张海东.建筑物不均匀沉降原因分析[J].工程勘 察,2013,1:114~119 地基建筑的施工阶段应结合场地条件合理布置临时设施以 .80. Industrial&Seience Trihune 201 8年第1 7卷第1 3期 将边坡稳定性系数视为基于这些随机变量的综合随机变量, 理论上可以根据基本随机变量的概率分布和函数形式来确 定边坡稳定性系数的概率分布。但是,在考虑地震动参数随 机特性的情况下,这种解析过程过于复杂,难以实现。因此 本文采用模特卡洛法(Monte Carlo Me ̄od)计算边坡的可 靠度。 (一)边坡物理模型简化。某核电厂场地平整后在厂区 形成多段人工边坡,人工边坡长度共计约1,200m,坡高lm~ 77m不等,边坡综合坡率约为1:1,边坡主体为岩质边坡。该 核电厂厂址进行了完整的边坡岩土工程勘察和地震危险性 评价工作,具备开展核电厂边坡可靠度分析条件。图3为边 坡代表性断面。依据边坡工程地质勘察成果,边坡岩土体可 分为四个岩组,各岩土层的物理力学参数见表1,设定岩土体 强度参数随机特征符合正态分布。_3 表1 各工程地质岩组强度、变形参数表 内聚力c(kPa) 内摩擦角‘p(。) 岩组 第一组粉 质粘土 第二组强 风化岩体 第三组中 风化岩体 第四组微 风化岩体 (均值 标准方差 均值 标准方差 kN/m ) Ixc El"C 。 仃。 18.5 26.7 9.8 25.50 2.5O 19.3 0 25.0 l8O.0 1O 72 32.O 28.O 6.4 5.6 23.图1 地震危险性评价概率方法工作流程图 1J为地震年平均发生率;b为震级频度关系式中的系数;M0为起算 震级;M 为潜在震源区震级上限;A为地震动峰值加速度;PT为T 年内超越概率 25.O 8o0 320 37.O 18.5 (二)地震参数处理。 1.地震动加速度超越概率函数。常规情况下,核电厂地 = 1.0E-2 震安全性评价中最终提供厂址地震动峰值加速度的超越概 曲线。该曲线一般仅涵盖地震动加速度17gal以上段,即大 舌 8 1 0 ̄-3 震级段,如图1。若将地震发生视为随机事件,大震级地震在 场所有地震中仅占很少部分,以图2为例,厂址水平向地震 动加速度17gal以上地震仅占约2%。因此进行概率分析必 须补充小震级段概率分布函数。对于地震动加速度大于 17gal段,经过对某核电厂安评中地震动加速度超越概率曲线 0 1 0E-4 专 T0 1UU 拟合,超越概率函数可表示为P (A>a)= Sa(gaI) e xln(a)2“176 xIn(a) · ;小于等于l7gal区段,假定分布函 图2某核电厂水平向地震动峰值加速度年超越概率曲线 数为线性,则超越概率函数可表示为PT(A>a)=1—0.05771 ×a。最终厂址地震动峰值加速度超越概率函数可表示为公 式1,厂址地震动加速度超越概率曲线见图4和图5。 fPT(A>a)=e一。·74 “( ) · ’ ln(a) o24(a>17ga1) 蒙特卡洛法,又称统计试验法或随机抽样技巧法。它适 用于随机变量的概率密度分布形式已知或符合假定的情况, 在目前可靠度计算中是一种相对精确的方法。_2 蒙特卡洛法 是从概率的角度出发来求解失效概率的,先对影响可靠度的 随机变量进行大量随机抽样,然后将这些抽样值逐个代人功 (A>a):1-0.05771 X a(0<as17 ga1)一 ( ) 能函数,累计功能函数值小于零的个数,由此确定失效频率。 四、以某核电厂边坡为对象开展可靠度分析 图4某核电厂址水平向地震动加速度超越概率曲线 2.地震动加速度概率分布函数。地震动加速度概率分 ∞ .0 ∞80 100 ,2o D…∞ “I 1∞ 180 200 220 撕 布函数可以由超越概率函数换算而来,具体见公式(2),厂址 水平向地震动加速度概率分布曲线见图6。 F(a)=1一PT(A>a) 图3边坡可靠度分析典型断面 Industrial&Science Tribune .R1. 201 8年第17卷第1 3期 第三,采用蒙特卡洛法抽样计算时,抽样次数N≥100/ 一O ~0 0 7 9 8 褂鞋悟求 0 6 O 5 0 4 P0,P0为预估的失效概率。 0 3 O 2 0 0 图5水平向地震动加速度超越概率曲线 (加速度大于17gal段) fl—e一。 “‘ “ “‘ 一 (a>17ga1) 【O.05771×a(0<a≤17ga1) (2) Fr0∞幕 v(邑 ∞ ∞ 巧 ∞ 图6 某核电厂址水平向地震动加速度概率分布曲线 3.地震动加速度概率密度函数。边坡可靠度分析计算 需要输入的是地震动参数的概率密度函数,将概率分布函数 求导即可得到概率密度函数,详见公式(3)。厂址水平向地 震动加速度概率密度曲线见图7。 F(a):F (a) r e一“ ‘ 、 … ‘ 一 ×(一0.74 x2 xln(a)+an176+a)(a>17 )… i0.05771(0<a---17ga1) { 0 04 伸 甚0 03 图7某核电厂址水平向地震动加速度概率密度曲线 (三)计算流程和结果。边坡可靠度分析流程分为:建立 计算模型;岩土、地震动参数输入;最危险滑面搜索;蒙特卡 洛抽样计算;结果统计。限于篇幅,计算具体过程不再赘述。 计算过程中部分关键因素作如下设置。 第一,最危险滑面仅搜索可能影响边坡整体稳定的圆弧 形滑面,采用摩根斯坦方法计算各滑面安全系数,取最小安 全系数滑面作为最危险滑面,之后开展可靠度分析。 第二,计算过程中同时施加水平向和竖直向地震系数, 竖直向地震系数取水平向2/3,且方向向上。 .R9. 第四,抽样计算安全系数Fs<1时判断为边坡失效。 第五,计算工况1只考虑岩土强度参数的不确定性;计 算工况2同时考虑岩土强度参数和地震动参数的不确定性。 计算结果见表2,图8为计算工况1边坡安全系数概率 密度函数图,图9为计算工况2边坡安全系数概率密度图。 表2边坡可靠度分析成果表 考虑的不确 失效 安全系数 安全系数 抽样计算 定性因素 概率 均值 标准差 次数 工况1 C、‘p 0.0020106 3.2623 0.906 100 0o0 工况2 c、‘p和地震 0动参数 .0oO027 4.6284 1.2 10o0 o00 Fr 誊 m 0 0 75 3 95 7 15 10 35 13 55 16 75 19 95 23 15 26 35 29 55 Factor of Safety 图8考虑岩土强度参数不确定性的安全系数概率密度 。。 幅oFac翟:6of s 0 34。40 46。 刖。 考虑岩土强度和地震动参数不确定性的安全系数概率密度 五、结语 本文初步探索了核电厂边坡可靠度分析方法,首次将地 震动加速度的随机性特征考虑到核电厂边坡可靠度分析中。 经过对某核电厂边坡典型断面的试算,验证了该方法是可行 的,为今后核电厂边坡可靠度分析提供了参考。 某核电厂边坡典型断面计算结果表明,在考虑地震动参 数和岩土强度参数不确定性的情况下,边坡失效概率为 0.000027;在只考虑岩土强度参数不确定性情况下,边坡失效 概率为0.002。可见核电厂边坡设计中地震动参数考虑了非 常高的保守性,对失效概率影响较大。 【参考文献】 [1]胡聿贤,地震安全性评价技术教程[M].北京:地震出版 社.1999:222~223 [2]陈祖煜,土质边坡稳定性分析[M].北京:中国水利水电 出版社,2003:292~294 [3]高谦,土木工程可靠性理论及其应用[M].北京:中国建 筑工业出版社,2007:5一一6 IndustHal&Seience Trihlille