一种边坡滑面搜索方法、系统、设备及可读存储介质

摘要

本发明提供了一种边坡滑面搜索方法、系统、设备及可读存储介质，涉及边坡安全性研究技术领域，包括划分边坡得到若干个截面，依次对每个截面进行建模和网格划分，得到若干个网格节点；计算每个截面中的目标滑面得到目标控制点组合，所述目标控制点均位于网格节点上；划分边坡得到若干个标准区间集，一个标准区间集由若干个截面的标准区间构成；将所有目标控制点坐标的平均值所在位置作为当前标准区间的集合点；连接所有标准区间集的曲线得到边坡的目标滑面。本发明用于解决现有技术在搜索滑面时需要生成大量随机参数，计算时间较长，且局限于简单的二维滑面搜索，难以进行三维滑面搜索和非圆弧任意滑面的搜索的技术问题。

权利要求

1.一种边坡滑面搜索方法，其特征在于，包括：划分边坡得到若干个截面，依次对每个截面进行建模和网格划分，得到若干个网格节点；

计算每个截面中的目标滑面得到目标控制点组合，所述目标控制点均位于网格节点上；

具体包括：

设定初始控制点的数量范围i∈{i1，i2，…iq}，依次选取所述范围中的值i作为初始控制点的数量；

选取i个网格节点作为初始控制点，由i个初始控制点(Cix1,Ciy1)，(Cix2,Ciy2),…，(Cixi,Ciyi)构成新的控制点组合，得到由所述控制点组合构成的滑面；

随机移动一个初始控制点(Cix,Ciy)至相邻的网格节点，得到新的控制点组合和由新的控制点组合构成的滑面；

根据滑面计算得到截面的安全系数CiKs，将安全系数CiKs最小的滑面作为目标滑面，并将构成目标滑面的控制点组合作为目标控制点组合；

划分边坡得到若干个标准区间集，一个标准区间集由若干个截面的标准区间构成；

具体包括：

确定标准区间集的数量；

根据标准区间集的数量计算边坡划分的间距： 其中，L表示边坡的长度，e表示间距，m表示标准区间集的数量；

划分边坡得到若干个标准区间集 其中，每个标准区间集均包含若干个截面的标准区间；

根据每个标准区间集的集合点拟合曲线，将所有目标控制点坐标的平均值所在位置作为当前标准区间的集合点；

根据每个标准区间集的集合点拟合曲线，连接所有标准区间集的曲线得到边坡的目标滑面。

2.根据权利要求1所述的边坡滑面搜索方法，其特征在于,所述根据滑面计算得到截面的安全系数CiKs，将安全系数CiKs最小的滑面作为目标滑面，并将构成目标滑面的控制点组合作为目标控制点组合，具体包括：根据初始滑面计算得到截面的初始安全系数，令初始安全系数为目标安全系数CiKs；

根据当前控制点组合构成的滑面计算得到截面的当前安全系数判断当前安全系数 与目标安全系数CiKs的大小关系：若 则令 并继续随机移动其他初始控制点至相邻的网格节点得到当前控制点组合；

若 则再次移动所述初始控制点(Cix,Ciy)至其他相邻的网格节点得到当前控制点组合；

若 则判断初始安全系数CiKs不变的次数是否超过预设稳定迭代次数：若超过，则将当前控制点组合作为目标控制点组合；否则，继续随机移动其他初始控制点至相邻的网格节点得到当前控制点组合；

根据当前控制点组合构成的滑面计算得到截面的当前安全系数 并重复判断当前安全系数 与目标安全系数CiKs的大小关系。

3.根据权利要求1所述的边坡滑面搜索方法，其特征在于,所述根据每个标准区间集的集合点拟合曲线，将所有目标控制点坐标的平均值所在位置作为当前标准区间的集合点，具体包括：获取位于标准区间内的目标控制点和目标控制点的坐标；

计算当前标准区间内所有目标控制点的横坐标平均值 和纵坐标平均值得到当前标准区间的集合点坐标为

4.根据权利要求1所述的边坡滑面搜索方法，其特征在于,所述根据每个标准区间集的集合点拟合曲线，连接所有标准区间集的曲线得到边坡的目标滑面,具体包括：获取任意一个标准区间集中的所有集合点；

将一个标准区间集中任意相邻的三个集合点拟合为曲线段，连接所有曲线段后得到当前标准区间集的曲线；

连接所有标准区间集的曲线得到边坡的目标滑面。

5.根据权利要求4所述的边坡滑面搜索方法，其特征在于,所述任意相邻的三个集合点拟合为曲线段，具体包括：2

设拟合曲线为f＝Ax+Bx+C，式中，A、B、C均为拟合参数；

获取任意相邻的三个集合点的坐标，利用三个集合点的坐标计算得到A、B、C的值，得到由所述相邻的三个集合点构成的拟合曲线。

6.一种边坡滑面搜索系统，其特征在于，包括：网格划分模块：用于划分边坡得到若干个截面，依次对每个截面进行建模和网格划分，得到若干个网格节点；

计算模块：用于计算每个截面中的目标滑面得到目标控制点组合，所述目标控制点均位于网格节点上；

具体包括：

设定初始控制点的数量范围i∈{i1,i2,…iq}，依次选取所述范围中的值i作为初始控制点的数量；

选取i个网格节点作为初始控制点，由i个初始控制点(Cix1,Ciy1)，(Cix2,Ciy2),…，(Cixi,Ciyi)构成新的控制点组合，得到由所述控制点组合构成的滑面；

随机移动一个初始控制点(Cix，Ciy)至相邻的网格节点，得到新的控制点组合和由新的控制点组合构成的滑面；

根据滑面计算得到截面的安全系数CiKs，将安全系数CiKs最小的滑面作为目标滑面，并将构成目标滑面的控制点组合作为目标控制点组合；

标准区间划分模块：用于划分边坡得到若干个标准区间集，一个标准区间集由若干个截面的标准区间构成；

具体包括：

确定标准区间集的数量；

根据标准区间集的数量计算边坡划分的间距： 其中，L表示边坡的长度，e表示间距，m表示标准区间集的数量；

划分边坡得到若干个标准区间集 其中，每个标准区间集均包含若干个截面的标准区间；

集合点确定模块：用于根据每个标准区间集的集合点拟合曲线，将所有目标控制点坐标的平均值所在位置作为当前标准区间的集合点；

目标滑面确定模块：用于根据每个标准区间集的集合点拟合曲线，连接所有标准区间集的曲线得到边坡的目标滑面。

7.一种边坡滑面搜索设备，其特征在于，包括：存储器，用于存储计算机程序；

处理器，用于执行所述计算机程序时实现如权利要求1至5任一项所述边坡滑面搜索方法的步骤。

8.一种可读存储介质，其特征在于：所述可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现如权利要求1至5任一项所述边坡滑面搜索方法的步骤。

说明书

一种边坡滑面搜索方法、系统、设备及可读存储介质

技术领域

[0001] 本发明涉及边坡安全性研究技术领域，具体而言，涉及一种边坡滑面搜索方法、系统、设备及可读存储介质。

背景技术

[0002] 覆盖型边坡是常见的一种边坡类型，广泛分布在我国中西部地区，严重威胁铁路建设、运营的安全性。当前，对于边坡的滑面搜索方法大多是圆弧形搜索方法，该方法的思路为：假设地基土的滑动面呈圆弧形，沿该圆弧滑动面发生滑动破坏的安全系数为总剪切力与总抗剪力之比。取不同的圆弧面，可得到不同的安全系数值，通过试算可以找到最危险的圆弧滑动面，并可确定最小的安全系数值。但该方法在搜索滑面时需要生成大量随机参数，计算时间较长，且局限于简单的二维滑面搜索，难以进行三维滑面搜索和非圆弧任意滑面的搜索。

发明内容

[0003] 本发明的目的在于提供一种边坡滑面搜索方法、系统、设备及可读存储介质，以改善上述问题。为了实现上述目的，本发明采取的技术方案如下：

[0004] 第一方面，本申请提供了一种边坡滑面搜索方法，包括：

[0005] 划分边坡得到若干个截面，依次对每个截面进行建模和网格划分，得到若干个网格节点；

[0006] 计算每个截面中的目标滑面得到目标控制点组合，所述目标控制点均位于网格节点上；

[0007] 划分边坡得到若干个标准区间集，一个标准区间集由若干个截面的标准区间构成；

[0008] 根据每个标准区间集的集合点拟合曲线，将所有目标控制点坐标的平均值所在位置作为当前标准区间的集合点；

[0009] 根据每个标准区间集的集合点拟合曲线，连接所有标准区间集的曲线得到边坡的目标滑面。

[0010] 进一步的，所述计算每个截面中的目标滑面得到目标控制点组合，所述控制点均位于网格节点上，具体包括：

[0011] 设定初始控制点的数量范围 ，依次选取所述范围中的值 作为初始控制点的数量；选取 个网格节点作为初始控制点，由 个初始控制点 ，，… ， 构成新的控制点组合，得到由所述控制点组合构成的滑面；

[0012] 随机移动一个初始控制点至 相邻的网格节点，得到新的控制点组合和由新的控制点组合构成的滑面；

[0013] 根据滑面计算得到截面的安全系数 ，将安全系数 最小的滑面作为目标滑面，并将构成目标滑面的控制点组合作为目标控制点组合。

[0014] 进一步的，所述根据滑面计算得到截面的安全系数 ，将安全系数 最小的滑面作为目标滑面，并将构成目标滑面的控制点组合作为目标控制点组合，具体包括：

[0015] 根据初始滑面计算得到截面的初始安全系数，令初始安全系数为目标安全系数；

[0016] 根据当前控制点组合构成的滑面计算得到截面的当前安全系数 ；

[0017] 判断当前安全系数 与目标安全系数 的大小关系：

[0018] 若 ，则令 ，并继续随机移动其他初始控制点至相邻的网格节点得到当前控制点组合；

[0019] 若 ，则再次移动所述初始控制点 至其他相邻的网格节点得到当前控制点组合；

[0020] 若 ，则判断初始安全系数 不变的次数是否超过预设稳定迭代次数：若超过，则将当前控制点组合作为目标控制点组合；否则，继续随机移动其他初始控制点至相邻的网格节点得到当前控制点组合；

[0021] 根据当前控制点组合构成的滑面计算得到截面的当前安全系数 ，并重复判断当前安全系数 与目标安全系数 的大小关系。

[0022] 进一步的，所述划分边坡得到若干个标准区间集，一个标准区间集由若干个截面的标准区间构成，具体包括：

[0023] 确定标准区间集的数量；

[0024] 根据标准区间集的数量计算边坡划分的间距： ，其中， 表示边坡的长度，表示间距， 表示标准区间集的数量；

[0025] 划分边坡得到若干个标准区间集 ， ，…， ，其中，每个标准区间集均包含若干个截面的标准区间。

[0026] 进一步的，所述根据每个标准区间集的集合点拟合曲线，将所有目标控制点坐标的平均值所在位置作为当前标准区间的集合点，具体包括：

[0027] 获取位于标准区间内的目标控制点和目标控制点的坐标；

[0028] 计算当前标准区间内所有目标控制点的横坐标平均值 和纵坐标平均值 ；

[0029] 得到当前标准区间的集合点坐标为 ；

[0030] 重复上述步骤，得到所有标准区间中的集合点。

[0031] 进一步的，所述根据每个标准区间集的集合点拟合曲线，连接所有标准区间集的曲线得到边坡的目标滑面,具体包括：

[0032] 获取任意一个标准区间集中的所有集合点；

[0033] 将一个标准区间集中任意相邻的三个集合点拟合为曲线段，连接所有曲线段后得到当前标准区间集的曲线；

[0034] 重复上述步骤，得到所有标准区间集的曲线；

[0035] 连接所有标准区间集的曲线得到边坡的目标滑面。

[0036] 进一步的，所述将任意相邻的三个集合点拟合为曲线段，具体包括：

[0037] 设拟合曲线为 ，式中， 均为拟合参数；

[0038] 获取任意相邻的三个集合点的坐标，利用三个集合点的坐标计算得到的值，得到由所述相邻的三个集合点构成的拟合曲线。

[0039] 第二方面，本申请还提供了一种边坡滑面搜索系统，包括：

[0040] 网格划分模块：用于划分边坡得到若干个截面，依次对每个截面进行建模和网格划分，得到若干个网格节点；

[0041] 计算模块：用于计算每个截面中的目标滑面得到目标控制点组合，所述目标控制点均位于网格节点上；

[0042] 标准区间划分模块：用于划分边坡得到若干个标准区间集，一个标准区间集由若干个截面的标准区间构成；

[0043] 集合点确定模块：用于根据每个标准区间集的集合点拟合曲线，将所有目标控制点坐标的平均值所在位置作为当前标准区间的集合点；

[0044] 目标滑面确定模块：用于根据每个标准区间集的集合点拟合曲线，连接所有标准区间集的曲线得到边坡的目标滑面。

[0045] 第三方面，本申请还提供了一种边坡滑面搜索设备，包括：

[0046] 存储器，用于存储计算机程序；

[0047] 处理器，用于执行所述计算机程序时实现所述边坡滑面搜索方法的步骤。

[0048] 第四方面，本申请还提供了一种可读存储介质，所述可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现上述基于边坡滑面搜索方法的步骤。

[0049] 本发明的有益效果为：

[0050] 本发明通过Advancing Front算法划分覆盖型边坡各层的网格，并提取网格节点，作为覆盖型边坡的滑面控制点的备选坐标，利用控制点对单个截面的进行滑面搜索与判定，并综合多个截面的滑面形成三维滑面的搜索，得到边坡中最危险的滑面。同时，通过搜索出的滑面可计算出边坡的安全系数，根据安全系数给出相应的边坡工程处置方案。本发明适用于任意水位线、任意地形、任意土层参数的二维截面或三维截面，计算方法合理、快捷，无需设定随机裂缝，打破了传统的滑面搜索法的限制。

[0051] 本发明的其他特征和优点将在随后的说明书阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明实施例了解。本发明的目的和其他优点可通过在所写的说明书、权利要求书、以及附图中所特别指出的结构来实现和获得。

附图说明

[0052] 为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。

[0053] 图1为本发明实施例中所述的边坡滑面搜索方法流程示意图;

[0054] 图2为本发明实施例中所述的控制点分布示意图;

[0055] 图3为本发明实施例中所述的拟合曲线图;

[0056] 图4为本发明实施例中所述的边坡滑面搜索系统结构示意图；

[0057] 图5为本发明实施例中所述的边坡滑面搜索设备结构示意图。

[0058] 图中标记：800、边坡滑面搜索设备；801、处理器；802、存储器；803、多媒体组件；804、I/O接口；805、通信组件。

具体实施方式

[0059] 为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。因此，以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

[0060] 应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。同时，在本发明的描述中，术语“第一”、“第二”等仅用于区分描述，而不能理解为指示或暗示相对重要性。

[0061] 实施例1

[0062] 本实施例提供了一种边坡滑面搜索方法。

[0063] 参见图1，图中示出了本方法：

[0064] S1.划分边坡得到若干个截面，依次对每个截面进行建模和网格划分，得到若干个网格节点；

[0065] 所述步骤S1具体包括：

[0066] S11.通过现场钻孔获取边坡的地层构造和水位线分布情况，并对边坡的土样进行力学实验分析得到土样的力学参数：含水率、粘聚力、内摩擦角、密度等。

[0067] S12.将边坡划分为k个截面。其中，边坡包括平滑区域和地形急剧变化区域。

[0068] 平滑区域的划分标准为：以5m为间距对边坡进行划分；

[0069] 地形急剧变化区域的划分标准为：以边坡坡度改变10°、地层线角度改变15°、水位线角度改变10°所处的位置为划分线对边坡进行划分。

[0070] S13.依次选取每一个边坡截面，使用ABAQUS软件对截面进行建模，并使用ABAQUS软件中的Advancing Front算法对边坡进行网格划分。需要注意的是，截面中包含有多种土层，每种土层至少划分为三层，为网格划分提供划分标准，并且网格尺寸需尽可能保持一致。

[0071] S2.计算每个截面中的目标滑面得到目标控制点组合，所述目标控制点均位于网格节点上；

[0072] 所述步骤S2具体包括：

[0073] S21.设定初始控制点的数量范围 ，依次选取所述范围中的值 作为初始控制点的数量；选取 个网格节点作为初始控制点，由 个初始控制点 ，，… ， 构成新的控制点组合，得到由所述控制点组合构成的滑面；

[0074] S23. 随机移动一个初始控制点至 相邻的网格节点，得到新的控制点组合和由新的控制点组合构成的滑面；

[0075] S24根据滑面计算得到截面的安全系数 ，将安全系数 最小的滑面作为目标滑面，并将构成目标滑面的控制点组合作为目标控制点组合。

[0076] 所述步骤S24具体包括：

[0077] S241.根据初始滑面计算得到截面的初始安全系数，令初始安全系数为目标安全系数 ；

[0078] S242.根据当前控制点组合构成的滑面计算得到截面的当前安全系数 ；

[0079] S243.判断当前安全系数 与目标安全系数 的大小关系：

[0080] 若 ，则令 ，并继续随机移动其他初始控制点至相邻的网格节点得到当前控制点组合；

[0081] 若 ，则再次移动所述初始控制点 至其他相邻的网格节点得到当前控制点组合；

[0082] 若 ，则判断初始安全系数 不变的次数是否超过预设稳定迭代次数：若超过，则将当前控制点组合作为目标控制点组合；否则，继续随机移动其他初始控制点至相邻的网格节点得到当前控制点组合；

[0083] S244.根据当前控制点组合构成的滑面计算得到截面的当前安全系数 ，并重复判断当前安全系数 与目标安全系数 的大小关系。

[0084] 以第N个截面为例，计算所述截面的目标控制点组合：

[0085] 1）设定第N个截面的控制点的数量范围 ，其中 ，式中，表示边坡的高度。

[0086] 2）在第一轮选取 个网格节点作为初始控制点，即第一轮设置有5个初始控制点，将所述5个初始控制点构成的面视作初始滑面；

[0087] 需要注意的是，请参阅图2，其中一个控制点应该位于边坡顶部，一个控制点应该位于坡脚，其余三个控制点应该按水平方向均匀分布在边坡上，边坡从上往下依次的控制点的坐标为： ， , ， ，并始终控制， 。

[0088] 3）根据初始滑面计算得到当前截面的初始安全系数 ，令所述初始安全系数为目标安全系数：

[0089] ；（1）

[0090] 式中， 表示第 块滑体、第 块滑体，所述滑体的底面为滑面，根据滑面的位置可得知滑体所在的位置， 表示第n个滑块的抗滑力， 表示第n个滑块的下滑力，表示第 块滑体的上一块滑体的抗滑力， 表示第 块滑体的上一块滑体的下滑力,  表示第 块滑体的传递系数。

[0091] 其中，下滑力的计算公式为：

[0092] ；（2）

[0093] 式中， 表示下滑力， 为滑块重量， 为滑块倾角， 为滑块上侧静水压力， 为滑块下侧静水压力， 为滑块上侧超孔隙水压力， 为滑块下侧超孔隙水压力， 为动水压力， 为水平地震力， 为竖向地震力。

[0094] 抗滑力的计算公式为：

[0095] ；（3）

[0096] 式中， 表示抗滑力， 为第 块滑体的上一块滑块的剩余下滑力， 为第 块滑体的上一块滑块的滑面倾角， 为滑面的静水压力， 为滑面的超孔隙水压力， 为粘聚力， 为滑面的长度。

[0097] 4）随机移动控制点 的至其相邻的网格节点 ，并利用公式（1）‑（3）计算当前控制点组合下截面的当前安全系数 ，

[0098] 5）比较 与 的大小：

[0099] 若 ，则 令 ，将 作 为，并继续随机移动初始控制点 至相邻的网格节点
得到新的控制点组合；

[0100] 若 ，则再次移动所述初始控制点 至其他相邻的网格节点得到当前控制点组合；

[0101] 若 ，则判断初始安全系数 不变的次数是否超过预设稳定迭代次数：

[0102] 若超过，则将当前控制点组合作为目标控制点组合；

[0103] 否则，继续随机移动其他初始控制点至相邻的网格节点得到当前控制点组合；

[0104] 优选的，令所述预设稳定迭代次数为20，若连续移动控制点20次所述目标安全系数 仍然不发生改变，则表示控制点数量 的情况下，当前的目标安全系数为当前截面的最小安全系数，其所代表的滑面为最危险滑面，即目标滑面，其对应的5个控制点组合为目标控制点组合。

[0105] 6）至此，第一轮计算完毕，进入第二轮、第三轮…第( 轮计算，重复上述步骤，计算出当控制点数量 的最小安全系数 、控制点数量 的最小安全系数，…，控制点数量 的最小安全系数 ，以及分别对应的目标滑面和目标控制点组合，最终得到( 个目标滑面和( 个目标控制点组合。

[0106] 同理，所有截面的目标滑面的计算方法如上步骤1）‑6），得到不同控制点组合下的目标滑面。

[0107] S3.划分边坡得到若干个标准区间集，一个标准区间集由若干个截面的标准区间构成；

[0108] 所述步骤S3具体包括：

[0109] S31.确定标准区间集的数量 ；

[0110] 式中， 表示标准区间集的数量。

[0111] S32.根据标准区间集的数量计算边坡划分的间距： ，其中， 表示边坡的长度，表示间距， 表示标准区间集的数量；

[0112] S33.划分边坡得到若干个标准区间集 ， ，…， ，其中，每个标准区间集均包含k个截面的标准区间，例如标准区间集 包含第1个截面的标准区间、第2个截面的标准区间 、…、第k个截面的标准区间 ，以此类推。

[0113] S4.根据每个标准区间集的集合点拟合曲线，将所有目标控制点坐标的平均值所在位置作为当前标准区间的集合点；

[0114] 所述步骤S4具体包括：

[0115] S41.获取位于标准区间内的目标控制点和目标控制点的坐标；

[0116] 具体的，由于每个截面包含( 个目标控制点组合，包含 个控制点，所述控制点分布于上述划分的标准区间上。

[0117] S42.计算当前标准区间内所有目标控制点的横坐标平均值 和纵坐标平均值 ；

[0118] S43.得到当前标准区间的集合点坐标为 ；

[0119] 需要说明的是，若当前区间内没有控制点时：

[0120] 令当前区间为第N个区间；

[0121] 此时，若第N‑1个区间的集合点为（ ）、第N+1个区间的集合点为（）；

[0122] 那么令第N个区间的集合点为 ；

[0123] 若第N‑1个区间的集合点为（ ）、第N+1个区间也没有集合点，那么将N个区间的集合点和N+1个区间的集合点设定在标准区间的中间位置，并计算出集合点的坐标。

[0124] S44.重复上述步骤，得到所有标准区间中的集合点。

[0125] S5.根据每个标准区间集的集合点拟合曲线，连接所有标准区间集的曲线得到边坡的目标滑面。

[0126] 所述步骤S5具体包括：

[0127] S51.获取任意一个标准区间集中的所有集合点，即第1个截面的集合点 、第2个截面的集合点 …第k个截面的集合点。

[0128] S52.将一个标准区间集中任意相邻的三个集合点拟合为曲线段，连接所有曲线段后得到当前标准区间集的曲线；

[0129] 所述步骤S52具体包括：

[0130] S521.设拟合曲线为 ，式中， 均为拟合参数；

[0131] S522.获取任意相邻的三个集合点的坐标 、 、 ，利用三个集合点的坐标计算得到 的值，得到由所述相邻的三个集合点构成的拟合曲线。

[0132] 同理，往后顺延获取三个集合点的坐标 、 、 ，利用三个集合点的坐标重新计算得到 的值，得到由所述相邻的三个集合点构成的拟合曲线。

[0133] 需要说明的是，由于b、c两点之间被拟合了两次，得到两个不同的曲线段1和曲线段2，请参阅图3，此时，需要将两条曲线段的平均值构成的最终曲线段作为b、c两点之间的拟合曲线，后续重复拟合的两点均采用上述方法进行处理。

[0134] S53.重复上述步骤，得到所有标准区间集的曲线；

[0135] S54.连接所有标准区间集的曲线得到边坡的目标滑面,所述滑面为本发明所搜索到的边坡的最危险滑面。

[0136] 实施例2：

[0137] 获取所有截面计算出的安全系数，其中一个截面包含( 个安全系数；

[0138] 选取最小安全系数作为边坡安全系数 ；

[0139] 根据边坡安全系数的划分标准，对边坡采用对应的防护措施：

[0140] 当边坡安全系数 ≥2时，不用采取特殊工程措施；

[0141] 当边坡安全系数2＞ ≥1.5时，监测边坡的变形速率，定期安排人员现场查看；

[0142] 当边坡安全系数1.5＞ ≥1.3时，加强监测边坡的变形速率，定期安排人员现场查看；

[0143] 当边坡安全系数1.3＞ ≥1.15时，采取边坡生态护坡措施，加强监测边坡的变形速率，定期安排人员现场查看；

[0144] 当边坡安全系数 ＜1.15时，采取边坡工程防护措施，如灌浆及勾缝、抹面、喷浆及喷射混凝土、喷锚网、干砌片石、浆砌片石、挡土墙以及土工合成材料防护，加强监测边坡的变形速率，定期安排人员现场查看。

[0145] 实施例3

[0146] 如图4所示，本实施例提供了一种边坡滑面搜索系统，所述系统包括：

[0147] 网格划分模块：用于划分边坡得到若干个截面，依次对每个截面进行建模和网格划分，得到若干个网格节点；

[0148] 计算模块：用于计算每个截面中的目标滑面得到目标控制点组合，所述目标控制点均位于网格节点上；

[0149] 标准区间划分模块：用于划分边坡得到若干个标准区间集，一个标准区间集由若干个截面的标准区间构成；

[0150] 集合点确定模块：用于根据每个标准区间集的集合点拟合曲线，将所有目标控制点坐标的平均值所在位置作为当前标准区间的集合点；

[0151] 目标滑面确定模块：用于根据每个标准区间集的集合点拟合曲线，连接所有标准区间集的曲线得到边坡的目标滑面。

[0152] 需要说明的是，关于上述实施例中的系统，其中各个模块执行操作的具体方式已经在有关该方法的实施例中进行了详细描述，此处将不做详细阐述说明。

[0153] 实施例4

[0154] 相应于上面的方法实施例，本实施例中还提供了一种边坡滑面搜索设备，下文描述的一种边坡滑面搜索设备与上文描述的一种边坡滑面搜索方法可相互对应参照。

[0155] 图5是根据示例性实施例示出的一种边坡滑面搜索设备800的框图。如图5所示，该边坡滑面搜索设备800可以包括：处理器801，存储器802。该边坡滑面搜索设备800还可以包括多媒体组件803， I/O接口804，以及通信组件805中的一者或多者。

[0156] 其中，处理器801用于控制该边坡滑面搜索设备800的整体操作，以完成上述的边坡滑面搜索方法中的全部或部分步骤。存储器802用于存储各种类型的数据以支持在该边坡滑面搜索设备800的操作，这些数据例如可以包括用于在该边坡滑面搜索设备800上操作的任何应用程序或方法的指令，以及应用程序相关的数据，例如联系人数据、收发的消息、图片、音频、视频等等。该存储器802可以由任何类型的易失性或非易失性存储设备或者它们的组合实现，例如静态随机存取存储器(Static Random Access Memory，简称SRAM)，电可擦除可编程只读存储器(Electrically Erasable Programmable Read‑Only Memory，简称EEPROM)，可擦除可编程只读存储器(Erasable Programmable Read‑Only Memory，简称EPROM)，可编程只读存储器(Programmable Read‑Only Memory，简称PROM)，只读存储器(Read‑Only Memory，简称ROM)，磁存储器，快闪存储器，磁盘或光盘。多媒体组件803可以包括屏幕和音频组件。其中屏幕例如可以是触摸屏，音频组件用于输出和/或输入音频信号。例如，音频组件可以包括一个麦克风，麦克风用于接收外部音频信号。所接收的音频信号可以被进一步存储在存储器802或通过通信组件805发送。音频组件还包括至少一个扬声器，用于输出音频信号。I/O接口804为处理器801和其他接口模块之间提供接口，上述其他接口模块可以是键盘，鼠标，按钮等。这些按钮可以是虚拟按钮或者实体按钮。通信组件805用于该边坡滑面搜索设备800与其他设备之间进行有线或无线通信。无线通信，例如Wi‑Fi，蓝牙，近场通信(Near FieldCommunication，简称NFC)，2G、3G或4G，或它们中的一种或几种的组合，因此相应的该通信组件805可以包括：Wi‑Fi模块，蓝牙模块，NFC模块。

[0157] 在一示例性实施例中，边坡滑面搜索设备800可以被一个或多个应用专用集成电路(Application  Specific  Integrated Circuit，简称ASIC)、数字信号处理器(DigitalSignal Processor，简称DSP)、数字信号处理设备(Digital Signal Processing Device，简称DSPD)、可编程逻辑器件(Programmable Logic Device，简称PLD)、现场可编程门阵列(Field Programmable Gate Array，简称FPGA)、控制器、微控制器、微处理器或其他电子元件实现，用于执行上述的边坡滑面搜索方法。

[0158] 在另一示例性实施例中，还提供了一种包括程序指令的计算机可读存储介质，该程序指令被处理器执行时实现上述的边坡滑面搜索方法的步骤。例如，该计算机可读存储介质可以为上述包括程序指令的存储器802，上述程序指令可由边坡滑面搜索设备800的处理器801执行以完成上述的边坡滑面搜索方法。

[0159] 实施例5

[0160] 相应于上面的方法实施例，本实施例中还提供了一种可读存储介质，下文描述的一种可读存储介质与上文描述的一种边坡滑面搜索方法可相互对应参照。

[0161] 一种可读存储介质，可读存储介质上存储有计算机程序，计算机程序被处理器执行时实现上述方法实施例的边坡滑面搜索方法的步骤。

[0162] 该可读存储介质具体可以为U盘、移动硬盘、只读存储器(Read‑Only Memory，ROM)、随机存取存储器(Random Access Memory，RAM)、磁碟或者光盘等各种可存储程序代码的可读存储介质。

[0163] 以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

[0164] 以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以权利要求的保护范围为准。