一种基于下游构筑物受损程度的边坡安全度评价方法与流程

本发明属于边坡安全评价与风险量化领域，尤其涉及一种基于下游构筑物受损程度的边坡安全度评价方法。
背景技术：
：滑坡是一种非常常见的地质灾害之一，尤其在夏季，我国每年都会遭受降雨引发的滑坡灾害。发生滑坡之后，滑动的土体和岩体对边坡下游的建筑物或构筑物产生冲击作用，造成不同程度的损害。发生在我国比较大型的滑坡有，千将坪古滑坡，重庆武隆滑坡和深圳光明新区余泥渣土滑坡，给我国造成了不同程度的损失。因此，考虑边坡下游构筑物受损程度的边坡安全评价就显得尤为重要。在边坡安全评价领域，一般采用极限平衡方法和有限元强度折减法进行最小安全系数的搜索，利用最小安全系数指标来评价边坡安全程度。然而，对于下游存在构筑物或建筑物的边坡而言，在评价其安全程度时，因为没有结合滑动体对建筑物或构筑物的损害程度以及边坡构筑材料的不确定性，而仅依靠安全系数，故不能直观地进行评价，不利于滑坡风险的防治。所以目前在评价下游存在构筑物或者建筑物的边坡安全程度时，在评价结果的合理和有效性上亟待提高。技术实现要素：根据以上现有技术的不足，本发明提供了一种基于下游构筑物受损程度的边坡安全度评价方法，其能有效解决现有技术中在评价下游存在构筑物或建筑物的边坡的安全度时不够有效、合理、直观，不利于滑坡风险防治的问题。本发明解决的技术问题采用的技术方案为：一种基于下游构筑物受损程度的边坡安全度评价方法，一种基于下游构筑物受损程度的边坡安全度评价方法，包括如下步骤：步骤1：针对所研究的边坡，因边坡材料力学参数的不确定性，基于材料力学参数的平均值和标准差σ，并根据材料力学参数3σ变化法则，确定力学参数的上下限值，并在此限值内按照均匀分布假定生成一系列参数组合值步骤2：基于第i组参数组合ci，i＝1,2...,n，计算其边坡的安全系数fsi，并采用光滑粒子流体动力学方法计算边坡滑动体的运动轨迹，并记录每个粒子所代表土体单元的最终位置；步骤3：根据下游边坡构筑物的位置和步骤2中得到的边坡滑动体的运动轨迹中最远点粒子位置的相对关系，计算两者之间的距离di，i＝1,2...,n；步骤4：重复步骤2和步骤3，得到一系列距离值和一系列的安全系数步骤5：计算一系列距离值的平均值，记为dm，计算一系列安全系数的平均值，记为fsm；步骤6：利用距离值的平均值dm和安全系数的平均值fsm的乘积对下游存在构筑物边坡的安全度进行评价。通过光滑粒子流体动力学方法确定滑动体轨迹，能够确定构筑物受损程度，并综合考虑边坡构筑材料的不确定性以及安全系数，最终评价边坡安全程度，相对于现有技术而言，能更加直观的评价下游存在构筑物或建筑物的边坡的安全程度，有利于滑坡风险的防治。其中，优选方式为：在步骤2中，利用极限平衡方法中的简化毕晓普法计算边坡的安全系数fsi，i＝1,2...,n。步骤3的具体实现为，下游边坡构筑物的位置坐标记为(xg,yg)，选取运动轨迹中最远点粒子，并确定其滑动后的位置坐标(xm,ym)；计算下游边坡构筑物和运动轨迹中最远点粒子两者之间的距离值di，其中，若运动轨迹中最远点粒子超过构筑物，则距离值di取负值，即否则，距离值di保持不变。在步骤5中，距离值的平均值dm具体为安全系数的平均值fsm具体为本发明具有以下有益效果：本发明首先通过光滑粒子流体动力学方法计算滑坡滑动体的运动轨迹；其次，判断滑动体最远点至构筑物的距离，并根据超过和未超过构筑物来判断距离的正负，并用该正负距离值初步评价边坡的安全度；再次，考虑边坡材料参数的不确定性，变换不同的材料参数重复以上两步计算，并保存每次计算时滑动体运动最远点至构筑物的正负距离值；最后利用该正负距离值的均值和安全系数来综合评价边坡安全度，相比于现有技术，评价下游存在构筑物或者建筑物的边坡安全度时更加有效、合理、直观，有利于滑坡风险的防治。附图说明图1是本发明所述基于下游构筑物受损程度的边坡安全度评价方法的流程图；图2是本发明所提供实施例的粘性土边坡与构筑物位置示意图；图3是本发明所提供实施例的光滑粒子流体动力学方法计算所得运动轨迹图。具体实施方式下面结合附图对本发明做进一步描述。实施例一：如图1～图3所示，本发明所述的一种基于下游构筑物受损程度的边坡安全度评价方法，包括一种基于下游构筑物受损程度的边坡安全度评价方法，其特征在于：包括如下步骤，步骤s1：针对所研究的边坡，因边坡材料力学参数的不确定性，基于材料力学参数的平均值和标准差σ，并根据材料力学参数3σ变化法则，确定力学参数的上下限值，并在此限值内按照均匀分布假定生成一系列参数组合值步骤s2：基于第i组参数组合ci，i＝1,2...,n，计算其边坡的安全系数fsi，并采用光滑粒子流体动力学方法计算边坡滑动体的运动轨迹，并记录每个粒子所代表土体单元的最终位置；步骤s3：根据下游边坡构筑物的位置和步骤2中得到的边坡滑动体的运动轨迹中最远点粒子位置的相对关系，计算两者之间的距离di，i＝1,2...,n；步骤s4：重复步骤2和步骤3，得到一系列距离值和一系列的安全系数步骤s5：计算一系列距离值的平均值，记为dm，计算一系列安全系数的平均值，记为fsm；步骤s6：利用距离值的平均值dm和安全系数的平均值fsm的乘积对下游存在构筑物边坡的安全度进行评价。通过光滑粒子流体动力学方法确定滑动体轨迹，能够确定构筑物受损程度，并综合考虑边坡构筑材料的不确定性以及安全系数，最终评价边坡安全程度，相对于现有技术而言，能更加直观的评价下游存在构筑物或建筑物的边坡的安全程度，有利于滑坡风险的防治。在步骤s2中，利用极限平衡方法中的简化毕晓普法计算边坡的安全系数fsi，i＝1,2...,n。步骤s3的具体实现为，下游边坡构筑物的位置坐标记为(xg,yg)，选取边坡滑动体的运动轨迹中最远点粒子，并确定其滑动后的位置坐标(xm,ym)；计算下游边坡构筑物和运动轨迹中最远点粒子两者之间的距离值di，其中，若运动轨迹中最远点粒子超过构筑物，则距离值di取负值，即否则，距离值di保持不变。在步骤s5中，距离值的平均值dm具体为安全系数的平均值fsm具体为以下是本发明所述方法在具体应用场景中的应用。步骤s1，如图2所示，某粘性土边坡，坡高4m，土的重度γ＝20kn/m3，内摩擦角土的粘聚力c均值为15.0kpa，考虑粘聚力c的不确定性，其标准差σ为1.67kpa，根据3σ法则，在10kpa到20kpa区间内，均匀生成n＝11个c值，即ci＝i+9.0,i＝1,2,....,11；步骤s2，在第一组参数c1下，即γ＝20kn/m3，c1＝10kpa，按照极限平衡方法中的简化毕晓普法计算其安全系数fs1为0.73；将图2所示边坡离散为9020个直径为0.1m的圆形粒子，将以上三个参数输入光滑粒子流体动力学方法程序中，得到滑动体运动轨迹图，如图3所示；步骤s3，从滑动体运动轨迹图中可以读出，滑动体运动轨迹上最远点粒子的位置大致位于xm＝5，ym＝2.0；水塔坐标为xg＝10.0，yg＝2.0，因此滑动体运动轨迹上的最远点与水塔之间的距离值d1＝5.0，又因为该滑动体运动轨迹上的最远点粒子超过了构筑物(即水塔)，因此距离值取负值，即d1＝-5.0。本发明中，利用距离值d1的值来量化下游构筑物(即水塔)的受损程度，该值越小表明构筑物(即水塔)的受损程度越大，反之亦然。步骤s4，针对其余10组参数组合，采取相同的计算方法和计算步骤，可以得到相应的安全系数值和相应的距离值(即受损程度)表1汇总了安全系数和距离值的详细信息；i1234567891011di-5.0-4.6-3.0-1.7-1.0-0.54-0.152.02.02.02.0fsi0.730.80.870.941.021.091.161.221.281.351.43表1步骤s5，基于上述所得安全系数值和距离值，分别计算距离平均值dm＝-0.73和安全系数的平均值fsm＝1.08。步骤s6，综合评定后的边坡安全度为-0.79。为了证明本发明方法的有效性，现进行对比分析，传统的边坡安全评价方法不考虑下游建筑物或构筑物的受损程度，仅仅用传统的安全系数来评价边坡的安全程度，对于本例，简化毕晓普法给出的安全系数平均值为1.08；本发明方法利用光滑粒子流体动力学方法确定滑动轨迹，通过与构筑物位置的比对判断其受损程度，受损程度平均值为-0.78，说明构筑物存在一定被撞击的风险，滑动体运动轨迹上的最远粒子距离其中心位置有0.78m，在边坡防护过程中需要引起重视。综上，通过对比发现，传统方法仅依靠安全系数来评价边坡安全程度，没有直观地反映对下游构筑物或建筑物的危害程度，不利于滑坡风险的防治，通过实例对比分析验证了本发明的有效性。以上所述为本发明的实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书以及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其他相关的
技术领域：
，均同理包括在本发明的专利保护范围内。当前第1页12