CN108509698B - 滑坡体土体的力学参数反演方法、装置及计算机设备 - Google Patents

Abstract

本发明适用于力学参数计算技术领域，提供了一种滑坡体土体的力学参数反演方法、装置及计算机设备，其中方法包括：获取滑坡体主滑面的实际埋深，选取与主滑线所在二维剖面等距的任意两个二维剖面，分别建立所述主滑线所在二维剖面与所述两个二维剖面的临界状态方程，计算得到主滑面埋深，若所述差值不小于所述预设阈值，重新选取两个等距剖面，以计算主滑面埋深，否则，提取当前所述容重、所述黏结力、所述摩擦系数的数值，实现了在参数反演过程中无需对各个参数影响坡体稳定性的敏感性进行分析，步骤简单，运算量少的技术效果。

Description

滑坡体土体的力学参数反演方法、装置及计算机设备

技术领域

本发明属于力学参数计算技术领域，尤其一种滑坡体土体的力学参数反演方法、装置及计算机设备。

背景技术

在滑坡体稳定性分析和工程设计中，滑带土的力学参数、尤其是黏结力和摩擦系数的合理取值直接关系着边坡工程的安全和工程造价。目前确定滑坡体土体的力学参数的方法主要有试验法(室内试验和现场试验)、工程类比法和反演方法，反演方法主要分为单参数反演方法和双参数反演方法，且在参数反演过程中需要对各个参数影响坡体稳定性的敏感性进行分析，计算量巨大。

发明内容

有鉴于此，本发明实施例提供了一种滑坡体土体的力学参数反演方法、装置及计算机设备，以解决现有技术中在参数反演过程中需要对各个参数影响坡体稳定性的敏感性进行分析，计算量巨大的问题。

本发明实施例的第一方面提供了一种滑坡体土体的力学参数反演方法，包括：获取滑坡体主滑面的实际埋深；

重复执行下述步骤，直至主滑面埋深与所述实际埋深的差值小于预设阈值：

根据滑坡体三维地形图，选取与主滑线所在二维剖面等距的任意两个二维剖面；分别建立所述主滑线所在二维剖面与所述两个二维剖面的临界状态方程，求取所述滑坡体土体的容重、黏结力与摩擦系数；根据所述容重、所述黏结力与所述摩擦系数建立所述主滑线所在二维剖面的地质力学模型，计算得到主滑面埋深；

若所述差值小于所述预设阈值，提取当前所述容重、所述黏结力、所述摩擦系数的数值。

可选地，所述获取滑坡体主滑面的实际埋深，包括：

根据深部位移监测结果获取滑坡体的主滑方向、主滑面位置与实际埋深。

可选地，所述分别建立所述主滑线所在二维剖面与所述两个二维剖面的临界状态方程，求取所述滑坡体土体的容重、黏结力与摩擦系数，包括：

分别建立所述主滑线所在二维剖面与所述两个二维剖面关于坡高、坡角、土体容重、黏结力与摩擦系数的临界状态方程；

确定所述主滑线所在二维剖面与所述两个二维剖面的坡高与坡角；

求取所述滑坡体土体的容重、黏结力与摩擦系数。

可选地，根据所述容重、所述黏结力与所述摩擦系数建立所述主滑线所在二维剖面的地质力学模型，计算得到主滑面埋深，包括：

根据所述容重、所述黏结力、所述摩擦系数建立所述主滑线所在二维剖面的地质力学模型；

根据所述地质力学模型，采用数值模拟分析法计算得到主滑面埋深。

可选地，所述临界状态方程的表达式为：

式中，γ表示土体容重；c表示黏结力；

表示摩擦系数；H表示坡高；m，n和h均为常数。

可选地，所述m，n和h通过下式确定：

式中，x表示坡角的正切。

本发明实施例的第二方面提供了一种滑坡体土体的力学参数装置，包括：

获取模块，用于获取滑坡体主滑面的实际埋深；

比较模块，用于重复执行下述步骤，直至主滑面埋深与所述实际埋深的差值小于预设阈值：根据滑坡体三维地形图，选取与主滑线所在二维剖面等距的任意两个二维剖面；分别建立所述主滑线所在二维剖面与所述两个二维剖面的临界状态方程，求取所述滑坡体土体的容重、黏结力与摩擦系数；根据所述容重、所述黏结力与所述摩擦系数建立所述主滑线所在二维剖面的地质力学模型，计算得到主滑面埋深；若所述差值小于所述预设阈值，提取当前所述容重、所述黏结力、所述摩擦系数的数值。

可选地，所述获取模块具体用于：

根据深部位移监测结果获取滑坡体的主滑方向、主滑面位置与实际埋深。

本发明实施例的第三方面提供了一种计算机设备，包括：存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现第一方面任一项所述方法的步骤。

本发明实施例的第四方面提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现第一方面任一项所述方法的步骤。

本发明实施例与现有技术相比存在的有益效果是：通过获取滑坡体主滑面的实际埋深，重复执行下述步骤，直至主滑面埋深与所述实际埋深的差值小于预设阈值，根据滑坡体三维地形图，选取与主滑线所在二维剖面等距的任意两个二维剖面，分别建立所述主滑线所在二维剖面与所述两个二维剖面的临界状态方程，求取所述滑坡体土体的容重、黏结力与摩擦系数；根据所述容重、所述黏结力与所述摩擦系数建立所述主滑线所在二维剖面的地质力学模型，计算得到主滑面埋深，若所述差值小于所述预设阈值，提取当前所述容重、所述黏结力、所述摩擦系数的数值，实现了在参数反演过程中无需对各个参数影响坡体稳定性的敏感性进行分析，运算量少的技术效果。

附图说明

为了更清楚地说明本发明实施例中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

图1是本发明实施例一提供的滑坡体土体的力学参数反演方法的流程图；

图2是本发明实施例二提供的滑坡体土体的力学参数反演方法的流程图；

图3是本发明实施例三提供的滑坡体土体的力学参数反演方法的流程图；

图4是本发明实施例四提供的滑坡体土体的力学参数反演装置的示意图；

图5是本发明实施例五提供的计算机设备的示意图。

具体实施方式

以下描述中，为了说明而不是为了限定，提出了诸如特定系统结构、技术之类的具体细节，以便透彻理解本发明实施例。然而，本领域的技术人员应当清楚，在没有这些具体细节的其它实施例中也可以实现本发明。在其它情况中，省略对众所周知的系统、装置、电路以及方法的详细说明，以免不必要的细节妨碍本发明的描述。

为了说明本发明所述的技术方案，下面通过具体实施例来进行说明。

实施例一

本发明实施例一提供一种滑坡体土体的力学参数反演方法。图1为本发明实施例一提供的滑坡体土体的力学参数反演方法的流程图。如图1所示，本实施例中的方法，可以包括：

步骤101、获取滑坡体主滑面的实际埋深；

滑坡是指斜坡上的土体或者岩体，受河流冲刷、地下水活动、雨水浸泡、地震及人工切坡等因素影响，在重力作用下，沿着一定的软弱面或软弱带，整体地或者分散地顺坡向下滑动的自然现象。滑坡体指滑坡的整个滑动部分，简称滑体。所述软弱面就是滑坡体的滑面，所述软弱带就是滑坡体的滑带，滑面通常是上陡下缓，近似圆弧形。滑面有时有一个，有时可有多个，根据滑动速度的快慢，将滑动最快的滑面作为主滑面，其余滑面作为次滑面。主滑面埋深是指滑坡体表面与主滑面之间的铅直距离。主滑线位于主滑面上，是指滑坡体滑动速度最快的纵向线，也叫滑坡轴。主滑线代表着一个滑坡整体滑动的方向，它位于滑坡体上推力最大、滑坡床凹槽最深的纵断面上。

在本实施例中，滑面可以替换为滑带。

可选地，获取滑坡体主滑面的实际埋深，可选定滑坡体主滑线的中部某一点或多点作为测试点，获取所述测试点的实际埋深，若滑坡体表面面积大、坡长较长，可以在所述主滑线的两端增加测试点，以增加计算结果的准确率。

步骤102、根据滑坡体三维地形图，选取与主滑线所在二维剖面等距的任意两个二维剖面；

本实施例中的二维剖面均为与水平面垂直的剖面。

可选地，所述滑坡体三维地形图采用测绘局提供的数据。

可选地，所述选取与主滑线所在二维坡面等距的任意两个二维剖面为随机选取，也可以设定初始值，按照递增或递减的顺序进行选取。

例如，第一次选取与主滑线所在二维剖面均距离n米的两个二维剖面，第二次选取与主滑线所在二维剖面均距离n+a米的两个二维剖面，第三次选取与主滑线所在二维剖面均距离n+2a米的两个二维剖面，第i次选取与主滑线所在二维剖面均距离n+(i-1)a米的两个二维剖面，n为正数，a为正数或负数。

本实施例在参数反演过程中，通过应用滑坡体三维地形图，考虑到了滑坡体形态的三维效应，使计算结果更合理、可靠。

步骤103、分别建立所述主滑线所在二维剖面与所述两个二维剖面的临界状态方程，求取所述滑坡体土体的容重、黏结力与摩擦系数；

现有技术中单参数反演假定一个参数已知的前提下，反算另外一个参数，需要选择对滑坡稳定性影响较敏感的参数作为未知参数。在本实施例中，针对每个剖面建立关于所述滑坡体土体的容重、黏结力与摩擦系数等三个未知参数的方程，针对所述主滑线所在二维剖面与所述两个二维剖面等三个二维剖面可以建立三个方程。未知参数的数量与方程数量相等，可以得到所述滑坡体土体的容重、黏结力与摩擦系数等三个未知数的解，无需进行参数敏感性分析。

步骤104、根据所述容重、所述黏结力与所述摩擦系数建立所述主滑线所在二维剖面的地质力学模型，计算得到主滑面埋深；

可选地，采用CAE(Computer Aided Engineering，计算机辅助工程)软件中的ABAQUS有限元软件根据所述容重、所述黏结力与所述摩擦系数建立所述主滑线所在二维剖面的地质力学模型，采用数值模拟分析方法计算得到主滑面埋深。

步骤105、比较计算得到的主滑面埋深与所述实际埋深的差值是否小于预设阈值，若所述差值小于预设阈值，执行步骤106，若所述差值不小于预设阈值，重新执行步骤102至步骤105；

可选地，所述预设阈值等于所述实际埋深乘以7％，也就是若计算得到的主滑面埋深与所述实际埋深的差值不小于所述实际埋深的7％时，重新执行步骤102至步骤105直至计算得到的主滑面埋深与所述实际埋深的差值小于所述实际埋深的7％时，执行步骤106。

步骤106、提取当前所述容重、所述黏结力、所述摩擦系数的数值。

当前所述容重、所述黏结力、所述摩擦系数的数值就是本实施例通过反演方法同时计算得到的三个力学参数的结果。

本实施例提供的滑坡体土体的力学参数反演方法，通过获取滑坡体主滑面的实际埋深，根据滑坡体三维地形图，选取与主滑线所在二维剖面等距的任意两个二维剖面，分别建立所述主滑线所在二维剖面与所述两个二维剖面的临界状态方程，求取所述滑坡体土体的容重、黏结力与摩擦系数，进而计算得到主滑面埋深，若所述差值不小于所述预设阈值，重新选取与主滑线所在二维剖面等距的两个二维剖面，否则，提取当前所述容重、所述黏结力、所述摩擦系数的数值，实现了在参数反演过程中无需对各个参数影响坡体稳定性的敏感性进行分析，运算量少的技术效果。

实施例二

本发明实施例二提供一种滑坡体土体的力学参数反演方法。本实施例是在实施例一提供的技术方案的基础上，通过对深部位移监测结果的分析，获取滑坡体的主滑方向、主滑面位置与实际埋深。

图2为本发明实施例二提供的滑坡体土体的力学参数反演方法的流程图。

如图2所示，本实施例中的方法，可以包括：

步骤201、根据深部位移监测结果获取滑坡体的主滑方向、主滑面位置与实际埋深。

具体地，所述深部位移监测结果为通过采用钻孔测斜仪得到的变形曲线与位移曲线，通过对两条曲线的对比分析确定滑坡体的主滑方向、主滑面位置与实际埋深，提高了数据的准确率。

步骤202、根据滑坡体三维地形图，选取与主滑线所在二维剖面等距的任意两个二维剖面；

步骤203、分别建立所述主滑线所在二维剖面与所述两个二维剖面的临界状态方程，求取所述滑坡体土体的容重、黏结力与摩擦系数；

步骤204、根据所述容重、所述黏结力与所述摩擦系数建立所述主滑线所在二维剖面的地质力学模型，计算得到主滑面埋深；

步骤205、比较计算得到的主滑面埋深与所述实际埋深的差值是否小于预设阈值，若所述差值小于预设阈值，执行步骤206，若所述差值不小于预设阈值，重新执行步骤202至步骤205；

步骤206、提取当前所述容重、所述黏结力、所述摩擦系数的数值。

本实施例中，步骤202至步骤206与实施例一中的步骤102至步骤106类似，此处不再赘述。

本实施例提供的滑坡体土体的力学参数反演方法，通过根据深部位移监测结果，获取滑坡体的主滑方向、主滑面位置与实际埋深，重复执行下述步骤，直至主滑面埋深与所述实际埋深的差值小于预设阈值，根据滑坡体三维地形图，选取与主滑线所在二维剖面等距的任意两个二维剖面，分别建立所述主滑线所在二维剖面与所述两个二维剖面的临界状态方程，求取所述滑坡体土体的容重、黏结力与摩擦系数，根据所述容重、所述黏结力与所述摩擦系数建立所述主滑线所在二维剖面的地质力学模型，计算得到主滑面埋深，若所述差值小于所述预设阈值，提取当前所述容重、所述黏结力、所述摩擦系数的数值，提高了数据的准确率，实现了在参数反演过程中无需对各个参数影响坡体稳定性的敏感性进行分析，步骤简单，运算量少的技术效果。

实施例三

本发明实施例三提供一种滑坡体土体的力学参数反演方法。本实施例是在实施例一提供的技术方案的基础上，分别建立所述主滑线所在二维剖面与所述两个二维剖面关于坡高、坡角、土体容重、黏结力与摩擦系数的临界状态方程，确定所述主滑线所在二维剖面与所述两个二维剖面的坡高与坡角，进而求取所述滑坡体土体的容重、黏结力与摩擦系数。

图3为本发明实施例三提供的滑坡体土体的力学参数反演方法的流程图。

如图3所示，本实施例中的方法，可以包括：

步骤301、获取滑坡体主滑面的实际埋深；

步骤302、根据滑坡体三维地形图，选取与主滑线所在二维剖面等距的任意两个二维剖面；

本实施例中，步骤301至步骤302与实施例一中的步骤101至步骤102类似，此处不再赘述。

步骤303、分别建立所述主滑线所在二维剖面与所述两个二维剖面关于坡高、坡角、土体容重、黏结力与摩擦系数的临界状态方程；

具体地，所述临界状态方程的表达式为：

式中，γ表示土体容重；c表示黏结力；

表示摩擦系数；H表示坡高；m，n和h均为常数。

所述m，n和h通过下式确定：

式中，x表示坡角的正切。

上述公式(1)与(2)是通过多个算例得到的经验公式，确定了容重、黏结力、摩擦系数与坡高、坡角的关系式，根据公式(1)与(2)为所述主滑线所在二维剖面与所述两个二维剖面建立方程组，计算过程简单，可快速得出结果。

步骤304、确定所述主滑线所在二维剖面与所述两个二维剖面的坡高与坡角；

具体地，根据滑坡体三维地形图确定所述主滑线所在二维剖面的坡高与坡角的数值，所述主滑线所在二维剖面的坡高为主滑线坡顶到坡底的距离，所述主滑线所在二维剖面的坡角为主滑线与水平面的平均夹角；根据滑坡体三维地形图确定所述两个二维剖面的坡高与坡角的数值，所述两个二维剖面分别与主滑面相交于两条滑线上，所述两个二维剖面的坡高分别为所述两条滑线的坡顶到坡底的距离，所述两个二维剖面的坡角分别为所述两条滑线与水平面的平均夹角。

步骤305、求取所述滑坡体土体的容重、黏结力与摩擦系数。

具体地，将所述主滑线所在二维剖面与所述两个二维剖面对应的三组坡高与坡角代入公式(1)与(2)，求取所述滑坡体土体的容重、黏结力与摩擦系数。

步骤306、根据所述容重、所述黏结力与所述摩擦系数建立所述主滑线所在二维剖面的地质力学模型，计算得到主滑面埋深；

步骤307、比较计算得到的主滑面埋深与所述实际埋深的差值是否小于预设阈值，若所述差值小于预设阈值，执行步骤308，若所述差值不小于预设阈值，重新执行步骤302至步骤307；

步骤308、提取当前所述容重、所述黏结力、所述摩擦系数的数值。

本实施例中，步骤306至步骤308与实施例一中的步骤104至步骤106类似，此处不再赘述。

本实施例通过获取滑坡体主滑面的实际埋深，重复执行下述步骤，直至主滑面埋深与所述实际埋深的差值小于预设阈值，根据滑坡体三维地形图，选取与主滑线所在二维剖面等距的任意两个二维剖面；分别建立所述主滑线所在二维剖面与所述两个二维剖面关于坡高、坡角、土体容重、黏结力与摩擦系数的临界状态方程，确定所述主滑线所在二维剖面与所述两个二维剖面的坡高与坡角，进而求取所述滑坡体土体的容重、黏结力与摩擦系数，根据所述容重、所述黏结力与所述摩擦系数建立所述主滑线所在二维剖面的地质力学模型，计算得到主滑面埋深，若所述差值小于所述预设阈值，提取当前所述容重、所述黏结力、所述摩擦系数的数值，结合了滑坡体的实际滑面动态与滑坡体的三维效应，使反演得到的力学参数更准确，实现了简单直接的求取所述滑坡体土体的容重、黏结力与摩擦系数等力学参数，省去了参数敏感性分析的步骤。

实施例四

本发明实施例四提供一种滑坡体土体的力学参数反演装置。图4为本发明实施例四提供的滑坡体土体的力学参数反演装置的示意图。如图4所示，本实施例中的滑坡体土体的力学参数反演装置，可以包括：

获取模块，用于获取滑坡体主滑面的实际埋深；

比较模块，用于重复执行下述步骤，直至主滑面埋深与所述实际埋深的差值小于预设阈值：根据滑坡体三维地形图，选取与主滑线所在二维剖面等距的任意两个二维剖面；分别建立所述主滑线所在二维剖面与所述两个二维剖面的临界状态方程，求取所述滑坡体土体的容重、黏结力与摩擦系数；根据所述容重、所述黏结力与所述摩擦系数建立所述主滑线所在二维剖面的地质力学模型，计算得到主滑面埋深；若所述差值小于所述预设阈值，提取当前所述容重、所述黏结力、所述摩擦系数的数值。

可选地，所述获取模块具体用于：

根据深部位移监测结果获取滑坡体的主滑方向、主滑面位置与实际埋深。

本实施例中的滑坡体土体的力学参数反演装置，可以用于执行上述任一实施例所述的滑坡体土体的力学参数反演方法，其具体实现原理可以参见上述实施例，此处不再赘述。

本实施例提供的滑坡体土体的力学参数反演装置，通过获取滑坡体主滑面的实际埋深，重复执行下述步骤，直至主滑面埋深与所述实际埋深的差值小于预设阈值，根据滑坡体三维地形图，选取与主滑线所在二维剖面等距的任意两个二维剖面，分别建立所述主滑线所在二维剖面与所述两个二维剖面的临界状态方程，求取所述滑坡体土体的容重、黏结力与摩擦系数；根据所述容重、所述黏结力与所述摩擦系数建立所述主滑线所在二维剖面的地质力学模型，计算得到主滑面埋深，若所述差值小于所述预设阈值，提取当前所述容重、所述黏结力、所述摩擦系数的数值，实现了在参数反演过程中无需对各个参数影响坡体稳定性的敏感性进行分析，运算量少的技术效果。

实施例五

图5为本发明实施例五提供的计算机设备的示意图。如图5所示，该实施例的计算机设备5包括：处理器50、存储器51以及存储在所述存储器51中并可在所述处理器50上运行的计算机程序52。所述处理器50执行所述计算机程序52时实现上述以计算机设备为执行主体的实施例中的步骤，例如图1所示的步骤101至106。或者，所述处理器50执行所述计算机程序52时实现实施例四中各模块/单元的功能，例如图4所示模块401至402的功能。

示例性的，所述计算机程序52可以被分割成一个或多个模块/单元，所述一个或者多个模块/单元被存储在所述存储器51中，并由所述处理器50执行，以完成本发明。所述一个或多个模块/单元可以是能够完成特定功能的一系列计算机程序指令段，该指令段用于描述所述计算机程序52在计算机设备5中的执行过程。例如，所述计算机程序52可以被分割成获取模块和比较模块，各模块具体功能如下：

获取模块，用于获取滑坡体主滑面的实际埋深；比较模块，用于重复执行下述步骤，直至主滑面埋深与所述实际埋深的差值小于预设阈值：根据滑坡体三维地形图，选取与主滑线所在二维剖面等距的任意两个二维剖面；分别建立所述主滑线所在二维剖面与所述两个二维剖面的临界状态方程，求取所述滑坡体土体的容重、黏结力与摩擦系数；根据所述容重、所述黏结力与所述摩擦系数建立所述主滑线所在二维剖面的地质力学模型，计算得到主滑面埋深；若所述差值小于所述预设阈值，提取当前所述容重、所述黏结力、所述摩擦系数的数值。

所述计算机设备5可以是桌上型计算机、笔记本、掌上电脑及云端服务器等计算设备。所述服务器可包括，但不仅限于，处理器50、存储器51。本领域技术人员可以理解，图5仅仅是计算机设备5的示例，并不构成对计算机设备5的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件，例如所述服务器还可以包括输入输出设备、网络接入设备、总线等。

所称处理器50可以是中央处理单元(Central Processing Unit，CPU)，还可以是其他通用处理器、数字信号处理器(Digital Signal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现成可编程门阵列(Field-Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。

所述存储器51可以是所述计算机设备5的内部存储单元，例如计算机设备5的硬盘或内存。所述存储器51也可以是所述计算机设备5的外部存储设备，例如所述计算机设备5上配备的插接式硬盘，智能存储卡(Smart Media Card,SMC)，安全数字(Secure Digital,SD)卡，闪存卡(Flash Card)等。进一步地，所述存储器51还可以既包括所述计算机设备5的内部存储单元也包括外部存储设备。所述存储器51用于存储所述计算机程序以及所述计算机设备5所需的其他程序和数据。所述存储器51还可以用于暂时地存储已经输出或者将要输出的数据。

本发明实施例还提供一种分布式光伏接入量计算机设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述任一项所述方法的步骤。

本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现上述任一项所述方法的步骤。

应理解，上述实施例中各步骤的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本发明实施例的实施过程构成任何限定。

所属领域的技术人员可以清楚地了解到，为了描述的方便和简洁，仅以上述各功能单元、模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能单元、模块完成，即将所述装置的内部结构划分成不同的功能单元或模块，以完成以上描述的全部或者部分功能。实施例中的各功能单元、模块可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中，上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。另外，各功能单元、模块的具体名称也只是为了便于相互区分，并不用于限制本申请的保护范围。上述系统中单元、模块的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。

在上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述或记载的部分，可以参见其它实施例的相关描述。

本领域普通技术人员可以意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。

在本发明所提供的实施例中，应该理解到，所揭露的装置/终端设备和方法，可以通过其它的方式实现。例如，以上所描述的装置/终端设备实施例仅仅是示意性的，例如，所述模块或单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通讯连接可以是通过一些接口，装置或单元的间接耦合或通讯连接，可以是电性，机械或其它的形式。

所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。

另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。

所述集成的模块/单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明实现上述实施例方法中的全部或部分流程，也可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一计算机可读存储介质中，该计算机程序在被处理器执行时，可实现上述各个方法实施例的步骤。其中，所述计算机程序包括计算机程序代码，所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。所述计算机可读介质可以包括：能够携带所述计算机程序代码的任何实体或装置、记录介质、U盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(Read-Only Memory，ROM)、随机存取存储器(RandomAccess Memory，RAM)、电载波信号、电信信号以及软件分发介质等。需要说明的是，所述计算机可读介质包含的内容可以根据司法管辖区内立法和专利实践的要求进行适当的增减，例如在某些司法管辖区，根据立法和专利实践，计算机可读介质不包括电载波信号和电信信号。

以上所述实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围，均应包含在本发明的保护范围之内。

Claims (9)

1.一种滑坡体土体的力学参数反演方法，其特征在于，包括：

获取滑坡体主滑面的实际埋深；

重复执行下述步骤，直至主滑面埋深与所述实际埋深的差值小于预设阈值：

根据滑坡体三维地形图，选取与主滑线所在二维剖面等距的任意两个二维剖面；分别建立所述主滑线所在二维剖面与所述两个二维剖面的临界状态方程，求取所述滑坡体土体的容重、黏结力与摩擦系数；根据所述容重、所述黏结力与所述摩擦系数建立所述主滑线所在二维剖面的地质力学模型，计算得到主滑面埋深；

若所述差值小于所述预设阈值，提取当前所述容重、所述黏结力、所述摩擦系数的数值；

其中，所述分别建立所述主滑线所在二维剖面与所述两个二维剖面的临界状态方程，包括，

分别建立所述主滑线所在二维剖面与所述两个二维剖面关于坡高、坡角、土体容重、黏结力与摩擦系数的临界状态方程；所述临界状态方程的表达式为：

式中，γ表示土体容重；c表示黏结力；

表示摩擦系数；H表示坡高；m，n和h均为常数。

2.如权利要求1所述的滑坡体土体的力学参数反演方法，其特征在于，所述获取滑坡体主滑面的实际埋深，包括：

根据深部位移监测结果获取滑坡体的主滑方向、主滑面位置与实际埋深。

3.如权利要求1所述的滑坡体土体的力学参数反演方法，其特征在于，所述求取所述滑坡体土体的容重、黏结力与摩擦系数，包括：

确定所述主滑线所在二维剖面与所述两个二维剖面的坡高与坡角；

求取所述滑坡体土体的容重、黏结力与摩擦系数。

4.如权利要求1所述的滑坡体土体的力学参数反演方法，其特征在于，根据所述容重、所述黏结力与所述摩擦系数建立所述主滑线所在二维剖面的地质力学模型，计算得到主滑面埋深，包括：

根据所述容重、所述黏结力、所述摩擦系数建立所述主滑线所在二维剖面的地质力学模型；

根据所述地质力学模型，采用数值模拟分析法计算得到主滑面埋深。

5.如权利要求1所述的滑坡体土体的力学参数反演方法，其特征在于，所述m，n和h通过下式确定：

式中，x表示坡角的正切。

6.一种滑坡体土体的力学参数反演装置，其特征在于，包括：

获取模块，用于获取滑坡体主滑面的实际埋深；

比较模块，用于重复执行下述步骤，直至主滑面埋深与所述实际埋深的差值小于预设阈值：根据滑坡体三维地形图，选取与主滑线所在二维剖面等距的任意两个二维剖面；分别建立所述主滑线所在二维剖面与所述两个二维剖面的临界状态方程，求取所述滑坡体土体的容重、黏结力与摩擦系数；根据所述容重、所述黏结力与所述摩擦系数建立所述主滑线所在二维剖面的地质力学模型，计算得到主滑面埋深；若所述差值小于所述预设阈值，提取当前所述容重、所述黏结力、所述摩擦系数的数值；

其中，所述比较模块，还用于：

分别建立所述主滑线所在二维剖面与所述两个二维剖面关于坡高、坡角、土体容重、黏结力与摩擦系数的临界状态方程；所述临界状态方程的表达式为：

式中，γ表示土体容重；c表示黏结力；

表示摩擦系数；H表示坡高；m，n和h均为常数。

7.如权利要求6所述的滑坡体土体的力学参数反演装置，其特征在于，所述获取模块具体用于：

根据深部位移监测结果获取滑坡体的主滑方向、主滑面位置与实际埋深。

8.一种计算机设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现如权利要求1至5任一项所述方法的步骤。

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

Images