一种基于均匀设计的边坡影响因素确定方法、装置与流程

1.本技术涉及数据处理技术领域，具体而言，涉及一种基于均匀设计的边坡影响因素确定方法、装置。


背景技术：

2.公路、铁路、水利等基础设施的兴建，形成了大量的边坡，日益增多的边坡事故给人民造成了巨大的损失。边坡的破坏是在多因素复合作用下形成的，各因素对边坡稳定性的影响程度决定滑动机制与破坏模式，对边坡稳定性影响因素进行敏感性分析对于边坡治理、监测和预报有重要意义。
3.边坡稳定性敏感性分析是目前的热点问题，很多专家学者对其进行了较深入的研究，归纳起来分为单因素分析法和多因素分析法。
4.单因素分析方法本质上都是选定一个指标值，变化其中一个因素，同时假定其它因素保持不变，比较基准值随因素变化的大小，然后在基准值
‑
因素变化曲线中直观的反映各因素的敏感性大小。该方法能够比较直观的反映各因素对基准值的影响，但是需一定的假设前提，与实际情况不符。
5.多因素分析法中有正交设计，基于正交设计的人工神经网络等方法，它们的共同特点是采用概率统计原理与电子计算机技术使试验次数和计算工作量有了一定程度的减小，但当因素和水平都较多时，试验次数和计算量依然很大，给试验和数值分析带来困难。


技术实现要素：

6.本技术实施例的目的在于提供一种基于均匀设计的边坡影响因素确定方法、装置，以改善目前试验次数和计算量大，给试验的和数值分析带来困难的问题。
7.本发明是这样实现的：
8.第一方面，本技术实施例提供一种基于均匀设计的边坡影响因素确定方法，包括：确定试验因素；其中，所述试验因素来源于待试验边坡结构模型；基于所述试验因素以及所述待试验边坡结构模型确定每个所述试验因素的取值范围；基于所述试验因素以及每个所述试验因素的取值范围，构建均匀设计表；其中，所述均匀设计表中包括试验号以及与每个所述试验号中的试验因素的数值；通过所述均匀设计表对所述试验因素进行试验，得到每个所述试验号对应的稳定性系数；基于每个所述试验号对应的稳定性系数，确定出所述试验因素中的边坡影响敏感因素。
9.在本技术实施例中，在确定试验因素以及每个试验因素的取值范围之后，构建均匀设计表；均匀设计能够很好的对影响边坡稳定性的各因素进行敏感性分析，特别是对于多因素多水平情况非常适用，可以大大减少试验次数和计算量，也使试验点在试验范围内分布更均匀，达到了较理想的效果。同时，每个试验因素的取值范围由试验因素以及待试验边坡结构模型所确定，可以提高后期分析结论的正确性和合理性。
10.结合上述第一方面提供的技术方案，在一些可能的实现方式中，通过如下步骤构
建所述均匀设计表；确定出试验次数n；获取向量组h＝(h1,h2,...h
m
)；其中，h1,h2,...h
m
均为小于n，且与n的最大公约数为1的数；生成所述均匀设计表的第j列，直至得到所述均匀设计表；其中，第j列的表达式为：u
ij
＝(i
×
h
j
)[mod n]；其中，[modn]表示同余运算；若i
×
h
j
超过n，则减去n的一个预设倍数，使落差在[1，n]之中；u
ij
通过递推算法生成；所述均匀设计表记为u
n
(n
m
)；对于确定的n和h所对应的m的个数由欧拉函数e(n)求得；当n为素数时，e(n)＝n
‑
1；当n为素数幂时，即n＝p
l
时，p为素数，l为正整数；当n既不是素数，也不是素数幂时，n表示为不同素数的方幂积；即此时，p1,p2,l p
m
为不同的素数，l1,l2,l l
m
为整数。
[0011]
在本技术实施例中，首先确定出试验次数n，然后获取小于n，且与n的最大公约数为1的数的集合，然后通过公式u
ij
＝(i
×
h
j
)[mod n]生成均匀设计表中的每一列，最后得到均匀设计表。通过该方式得到的均匀设计表使得每个因素的每个水平做一次且仅做一次试验；任何两个因素的试验点可标在平面的格子上，每行每列有且仅有一个试验点均匀设计表任意两列组成的试验方案均匀性一般不等价；当因素的水平增加时，试验次数按水平数的增加而增加。
[0012]
结合上述第一方面提供的技术方案，在一些可能的实现方式中，所述试验因素包括：容重、内摩擦角、黏聚力、坡角以及地震系数。
[0013]
结合上述第一方面提供的技术方案，在一些可能的实现方式中，所述试验次数n为所述试验因素的个数的两倍。
[0014]
在本技术实施例中，试验次数n为试验因素的个数的两倍，使得在保证减少试验次数和计算量的同时，也使试验点在试验范围内分布更均匀，达到了较理想的效果。
[0015]
结合上述第一方面提供的技术方案，在一些可能的实现方式中，所述通过所述均匀设计表对所述试验因素进行试验，得到每个所述试验号对应的稳定性系数，包括：通过所述均匀设计表对所述试验因素进行试验；通过简化毕肖普法得到每个所述试验号对应的稳定性系数；其中，所述稳定性系数的表达式为：
[0016]
q
hi
＝k
h
c
z
α
i
w
i
；其中，k为
[0017]
所述稳定性系数；c
i
为第i次试验对应的黏聚力；φ
i
为第i次试验对应的内摩擦角；w
i
为第i次试验分块自重；u
i
表示孔隙水压力；b
i
为第i次试验分块宽度；β
i
为第i次试验对应的坡角；q
hi
为水平地震力；k
h
为水平地震系数；c
z
为综合影响系数；α
i
为地震加速度分布系数。
[0018]
在本技术实施例中，通过简化毕肖普法得到每个试验号对应的稳定性系数，进而使得稳定性系数更加准确，能够更加直观可靠的表现出在各项试验因素下的边坡稳定性。
[0019]
结合上述第一方面提供的技术方案，在一些可能的实现方式中，所述基于每个所述试验号对应的稳定性系数，确定出所述试验因素中的边坡影响敏感因素，包括：采用回退
法对所述试验因素进行回归分析，建立回归方程；对所述回归方程进行显著性检验；基于所述显著性检验的结果，确定出所述试验因素中的边坡影响敏感因素。
[0020]
在本技术实施例中，采用回退法对试验因素进行回归分析，建立回归方程；对回归方程进行显著性检验；基于显著性检验的结果，确定出试验因素中的边坡影响敏感因素，通过该方式，能够对各项试验因素下的边坡稳定性进行有效的分析，进而提高了后续确定出准确的影响边坡稳定性敏感性的敏感因素。
[0021]
结合上述第一方面提供的技术方案，在一些可能的实现方式中，述对所述回归方程进行显著性检验，包括：对所述回归方程进行显著性检验，确定出所述试验因素中的不显著因素；将所述不显著因素剔除，重新构建新的回归方程；对所述新的回归方程进行显著性检验；相应的，所述基于所述显著性检验的结果，确定出所述试验因素中的边坡影响敏感因素，包括：基于对所述新的回归方程进行显著性检验的结果，确定出所述试验因素中的边坡影响敏感因素。
[0022]
在本技术实施例中，通过对回归方程进行显著性检验，确定出试验因素中的不显著因素；将不显著因素剔除，重新构建新的回归方程；对新的回归方程进行显著性检验，进一步的提高了显著性检验的可靠性。
[0023]
结合上述第一方面提供的技术方案，在一些可能的实现方式中，所述确定试验因素，包括：建立边坡数值以构建试验边坡结构模型；其中，所述模型为均质土坡；基于所述待试验边坡结构模型的各项数值，确定出所述试验因素。
[0024]
在本技术实施例中，先建立边坡数值构建试验边坡结构模型，然后通过，然后基于待试验边坡结构模型的各项数值，确定出试验因素；其中该模型为均质土坡，通过构建上述模型，以及基于上述模型确定试验因素来进行试验，对边坡成灾机理及防治具有重要意义，特别是对于从事三峡库区库岸边坡防治研究与应用具有重要的理论价值和实践价值。
[0025]
第二方面，本技术实施例提供一种基于均匀设计的边坡影响因素确定装置，包括：第一确定模块，用于确定试验因素；其中，所述试验因素来源于待试验边坡结构模型；第二确定模块，用于基于所述试验因素以及所述待试验边坡结构确定每个所述试验因素的取值范围；构建模块，用于基于所述试验因素以及每个所述试验因素的取值范围，构建均匀设计表；其中，所述均匀设计表中包括试验号以及与每个所述试验号中的试验因素的数值；处理模块，用于通过所述均匀设计表对所述试验因素进行试验，得到每个所述试验号对应的稳定性系数；第三确定模块，用于基于每个所述试验号对应的稳定性系数，确定出所述试验因素中的边坡影响敏感因素。
[0026]
第三方面，本技术实施例提供一种电子设备，包括：处理器和存储器，所述处理器和所述存储器连接；所述存储器用于存储程序；所述处理器用于调用存储在所述存储器中的程序，执行如上述第一方面实施例和/或结合上述第一方面实施例的一些可能的实现方式提供的方法。
[0027]
第四方面，本技术实施例提供一种存储介质，其上存储有计算机程序，所述计算机程序在被处理器运行时执行如上述第一方面实施例和/或结合上述第一方面实施例的一些可能的实现方式提供的方法。
附图说明
[0028]
为了更清楚地说明本技术实施例的技术方案，下面将对本技术实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本技术的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
[0029]
图1为本技术实施例提供的一种电子设备的结构框图。
[0030]
图2为本技术实施例提供的一种基于均匀设计的边坡影响因素确定方法的步骤流程图。
[0031]
图3为本技术实施例提供的一种待试验边坡结构模型的示意图。
[0032]
图4为本技术实施例提供的一种基于均匀设计的边坡影响因素确定装置的模块框图。
[0033]
图标：100
‑
电子设备；10
‑
处理器；11
‑
存储器；200
‑
基于均匀设计的边坡影响因素确定装置；201
‑
第一确定模块；202
‑
第二确定模块；203
‑
构建模块；204
‑
处理模块；205
‑
第三确定模块。
具体实施方式
[0034]
下面将结合本技术实施例中的附图，对本技术实施例中的技术方案进行描述。
[0035]
鉴于目前的边坡稳定性敏感性分析，当因素和水平都较多时，试验次数和计算量依然很大，给试验和数值分析带来困难等问题，本技术发明人经过研究探索，提出以下实施例以解决上述问题。
[0036]
请参阅图1，本技术实施例提供的一种基于均匀设计的边坡影响因素确定方法方法及装置的电子设备100的示意性结构框图。本技术实施例中，电子设备100可以是，但不限于个人计算机(personal computer，pc)、智能手机、平板电脑、个人数字助理(personal digital assistant，pda)、移动上网设备(mobile internet device，mid)等。在结构上，电子设备100可以包括处理器10和存储器11。
[0037]
处理器10与存储器11直接或间接地电性连接，以实现数据的传输或交互，例如，这些元件相互之间可通过一条或多条通讯总线或信号线实现电性连接。基于均匀设计的边坡影响因素确定装置包括至少一个可以软件或固件(firmware)的形式存储在存储器11中或固化在电子设备100的操作系统(operating system，os)中的软件模块。处理器10用于执行存储器11中存储的可执行模块，例如，基于均匀设计的边坡影响因素确定装置所包括的软件功能模块及计算机程序等，以实现基于均匀设计的边坡影响因素确定方法。处理器10可以在接收到执行指令后，执行计算机程序。
[0038]
其中，处理器10可以是一种集成电路芯片，具有信号处理能力。处理器10也可以是通用处理器，例如，可以是中央处理器(central processing unit，cpu)、数字信号处理器(digital signal processor，dsp)、专用集成电路(application specific integrated circuit，asic)、分立门或晶体管逻辑器件、分立硬件组件，可以实现或者执行本技术实施例中的公开的各方法、步骤及逻辑框图。此外，通用处理器可以是微处理器或者任何常规处理器等。
[0039]
存储器11可以是，但不限于，随机存取存储器(random access memory，ram)、只读
存储器(read only memory，rom)、可编程只读存储器(programmable read
‑
only memory，prom)、可擦可编程序只读存储器(erasable programmable read
‑
only memory，eprom)，以及电可擦编程只读存储器(electric erasable programmable read
‑
only memory，eeprom)。存储器11用于存储程序，处理器10在接收到执行指令后，执行该程序。
[0040]
需要说明的是，图1所示的结构仅为示意，本技术实施例提供的电子设备100还可以具有比图1更少或更多的组件，或是具有与图1所示不同的配置。此外，图1所示的各组件可以通过软件、硬件或其组合实现。
[0041]
请参阅图2，图2为本技术实施例提供的基于均匀设计的边坡影响因素确定方法的步骤流程图，该方法应用于图1所示的电子设备100。需要说明的是，本技术实施例提供的基于均匀设计的边坡影响因素确定方法不以图2及以下所示的顺序为限制，该方法包括：步骤s101
‑
步骤s105。
[0042]
步骤s101：确定试验因素；其中，试验因素来源于待试验边坡结构模型。
[0043]
步骤s102：基于试验因素以及待试验边坡结构模型确定每个试验因素的取值范围。
[0044]
步骤s103：基于试验因素以及每个试验因素的取值范围，构建均匀设计表；其中，均匀设计表中包括试验号以及与每个试验号中的试验因素的数值。
[0045]
步骤s104：通过均匀设计表对试验因素进行试验，得到每个试验号对应的稳定性系数。
[0046]
步骤s105：基于每个试验号对应的稳定性系数，确定出试验因素中的边坡影响敏感因素。
[0047]
在本技术实施例中，在确定试验因素以及每个试验因素的取值范围之后，构建均匀设计表；均匀设计能够很好的对影响边坡稳定性的各因素进行敏感性分析，特别是对于多因素多水平情况非常适用，可以大大减少试验次数和计算量，也使试验点在试验范围内分布更均匀，达到了较理想的效果。同时，每个试验因素的取值范围由试验因素以及待试验边坡结构模型所确定，可以提高后期分析结论的正确性和合理性。
[0048]
为了便于理解，下面首先对本技术实施例所提供的均匀设计表进行说明。具体的，通过如下步骤构建均匀设计表，包括：确定出试验次数n；获取向量组h＝(h1,h2,...h
m
)；其中，h1,h2,...h
m
均为小于n，且与n的最大公约数为1的数；生成均匀设计表的第j列，直至得到均匀设计表。
[0049]
其中，第j列的表达式为：
[0050]
u
ij
＝(i
×
h
j
)[mod n]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
[0051]
公式(1)中，[modn]表示同余运算；若i
×
h
j
超过n，则减去n的一个预设倍数，使落差在[1，n]之中。
[0052]
u
ij
通过递推算法生成；均匀设计表记为u
n
(n
m
)；对于确定的n和h所对应的m的个数由欧拉函数e(n)求得。具体的，欧拉函数e(n)的求取条件为：
[0053]
当n为素数时，e(n)＝n
‑
1。
[0054]
当n为素数幂时，即n＝p
l
时，p为素数，l为正整数。
[0055]
当n既不是素数，也不是素数幂时，n表示为不同素数的方幂积；即此时，p1,p2,l p
m
为不同的素数，l1,l2,l l
m
为整数。
[0056]
需要说明的是，上述三种情形，n为素数的情况最好，此时最多能获得n
‑
1列，即可安排n
‑
1个因素做试验，而当n为非素数时，在均匀设计表的结构中不可能有n
‑
1列，甚至可能的列数会更少。因此，可以将u
n
(n
m
)的最后一行去掉来构造的方法。由于比u
n
(n
m
)具有更好的均匀性，故一般会着重考虑
[0057]
此外，每个均匀表都对应有一个使用表。使用表示保证从u
n
(n
m
)或中所选的s列具有很好的均匀性，而均匀性是用偏差来衡量的。设x1,x2,l x
n
是c
m
中的n个点，任一向量x＝(x1,x2,l x
m
)∈c
m
，记为v(x)＝x1l x
m
为矩形[0,x]的体积，n
x
为x1,x2,l x
n
中落入[0,x]的点数，则称为点集{x1,x2,l x
n
}在c
m
中的偏差。
[0058]
下面对上述步骤进行举例说明。
[0059]
在步骤s101中确定试验因素，具体包括：建立边坡数值以构建试验边坡结构模型；其中，模型为均质土坡；基于待试验边坡结构模型的各项数值，确定出试验因素。
[0060]
于本技术实施例中，试验因素可以包括容重、内摩擦角、黏聚力、坡角以及地震系数中的任意一种或多种，当然，在其他实施例中，试验因素还可以是其他的任意因素，比如含水率等等，本技术不作限定。
[0061]
可以理解的是，该步骤中，先根据待试验的边坡构建对应的试验边坡结构模型，然后赋予对应的数值。之后，根据构建的待试验边坡结构模型的各项数值，确定出试验因素。也即，可以根据构建的待试验边坡结构模型的各项数值筛选出所需试验的因素，也可以是包含全部数值对应的因素，本技术不作限定。
[0062]
综上，在本技术实施例中，先建立边坡数值构建试验边坡结构模型，然后通过，然后基于待试验边坡结构模型的各项数值，确定出试验因素；其中该模型为均质土坡，通过构建上述模型，以及基于上述模型确定试验因素来进行试验，对边坡成灾机理及防治具有重要意义，特别是对于从事三峡库区库岸边坡防治研究与应用具有重要的理论价值和实践价值。
[0063]
可选地，本技术实施例中，试验次数n为试验因素的个数的两倍，使得在保证减少试验次数和计算量的同时，也使试验点在试验范围内分布更均匀，达到了较理想的效果。
[0064]
可以理解的是，在其他实施例中，试验次数n为试验因素的个数的三倍，本技术不作限定。
[0065]
在步骤s102中，基于试验因素以及待试验边坡结构模型确定每个试验因素的取值范围，可以以对试验边坡结构模型赋予的数值为基准值，分别在其上下浮动预设范围来确定出的。
[0066]
在步骤s104，通过均匀设计表对所述试验因素进行试验，得到每个试验号对应的稳定性系数，包括：通过均匀设计表对试验因素进行试验；通过简化毕肖普法(bisshop)得
到每个试验号对应的稳定性系数；
[0067]
其中，稳定性系数的表达式为：
[0068][0069]
公式(2)中，k为稳定性系数；c
i
为第i次试验对应的黏聚力；φ
i
为第i次试验对应的内摩擦角；w
i
为第i次试验分块自重；u
i
表示孔隙水压力；b
i
为第i次试验分块宽度；β
i
为第i次试验对应的坡角；q
hi
为水平地震力；k
h
为水平地震系数；c
z
为综合影响系数；α
i
为地震加速度分布系数。
[0070]
需要说明的是，毕肖普法是土坡稳定分析考虑土条间相互作用力的圆弧滑动分析法。仍然是基于极限平衡原理，把滑裂土体当作刚体绕圆心旋转，并分条计算其滑动力与抗滑力，最后求出稳定安全系数，计算时考虑了土条之间的相互作用力，是一种改进的圆弧滑动法。
[0071]
在本技术实施例中，通过简化毕肖普法得到每个试验号对应的稳定性系数，进而使得稳定性系数更加准确，能够更加直观可靠的表现出在各项试验因素下的边坡稳定性。
[0072]
在步骤s105中，基于每个试验号对应的稳定性系数，确定出试验因素中的边坡影响敏感因素，具体包括：采用回退法对试验因素进行回归分析，建立回归方程；对回归方程进行显著性检验；基于显著性检验的结果，确定出试验因素中的边坡影响敏感因素。
[0073]
在本技术实施例中，采用回退法对试验因素进行回归分析，建立回归方程；对回归方程进行显著性检验；基于显著性检验的结果，确定出试验因素中的边坡影响敏感因素，通过该方式，能够对各项试验因素下的边坡稳定性进行有效的分析，进而提高了后续确定出准确的影响边坡稳定性敏感性的敏感因素。
[0074]
可选地，上述对回归方程进行显著性检验，具体包括：对回归方程进行显著性检验，确定出试验因素中的不显著因素；将不显著因素剔除，重新构建新的回归方程；对新的回归方程进行显著性检验；相应的，上述的基于显著性检验的结果，确定出试验因素中的边坡影响敏感因素，包括：基于对新的回归方程进行显著性检验的结果，确定出试验因素中的边坡影响敏感因素。
[0075]
上述的回归方程的表达式为：
[0076]
y＝b0 b1×
γ b2×
φ b3×
c b4×
β b5×
α
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3)
[0077]
公式(3)中，b0＝1.75、b1＝
‑
2.5
×
10
‑2、b2＝2.14
×
10
‑2、b3＝4.62
×
10
‑2、b4＝
‑
4.34
×
10
‑2、b5＝
‑
3.18
×
10
‑4；α为显著性水平，取值为0.05。
[0078]
在本技术实施例中，通过对回归方程进行显著性检验，确定出试验因素中的不显著因素；将不显著因素剔除，重新构建新的回归方程；对新的回归方程进行显著性检验，进一步的提高了显著性检验的可靠性。
[0079]
为了便于理解，下面以一个具体的示例对本技术实施例所提供的一种基于均匀设计的边坡影响因素确定方法进行说明。
[0080]
请参阅图3，图3为建立边坡数值以构建出的试验边坡结构模型。其中，坡高h＝20m，土容重γ＝20kn/m3；内摩擦角φ＝17
°
，黏聚力c＝42kpa，地震系数α＝0.05g。
[0081]
然后，基于试验因素以及待试验边坡结构模型确定每个试验因素的取值范围。假设以上述的参数为基准值，范围为在其上下浮动预设数值。
[0082]
示例性的，每个试验因素的取值范围为：
[0083]
土容重γ的范围为18～27(kn/m3)；
[0084]
内摩擦角φ的范围为15～33(
°
)；
[0085]
黏聚力c的范围为27～54(kpa)；
[0086]
坡角β的范围为31～49(
°
)；可根据图示计算得到坡角的基准值；
[0087]
地震系数α的范围为10～100(cm/s2)。
[0088]
在确定出试验因素以及每个试验因素的取值范围之后，构建均匀设计表其中，试验次数n为试验因素的个数的两倍，即试验次数为10，然后通过均匀设计表对所述试验因素进行试验，得到每个试验号对应的稳定性系数。(具体的请参阅表1～3)
[0089]
表1均匀设计表
[0090][0091]
表2均匀设计表的使用表
[0092][0093]
表3均匀设计计算结果
[0094][0095]
最后基于每个试验号对应的稳定性系数，确定出试验因素中的边坡影响敏感因素。于本技术实施例中，可通过简化毕肖普法得到每个试验号对应的稳定性系数，然后采用回退法对试验因素进行回归分析，建立回归方程；对回归方程进行显著性检验，进而确定出试验因素中的边坡影响敏感因素。
[0096]
具体的，采用后退法进行回归分析，取显著性水平α为0.05。
[0097]
回归方程的表达式为：
[0098]
y＝b0 b1×
γ b2×
φ b3×
c b4×
β b5×
α
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(4)
[0099]
公式(4)中，b0＝1.75、b1＝
‑
2.5
×
10
‑2、b2＝2.14
×
10
‑2、b3＝4.62
×
10
‑2、b4＝
‑
4.34
×
10
‑2、b5＝
‑
3.18
×
10
‑4；该公式(4)与公式(3)相同。
[0100]
回归方程显著性检验可参考表4。
[0101]
表4变量分析表
[0102][0103][0104]
样本容量n＝10，检验值ft＝52.20，临界值f(0.05,5,4)＝6.265，ft＞f(0.05,5,4)，回归方程显著。
[0105]
f检验值：f(1)＝5.992，f(2)＝17.50，f(3)＝184.3，f(4)＝72.27，f(5)＝9.707
×
10
‑
2，各方程项对回归的贡献大小依次为c、β、φ、γ、α，方程项α对回归的贡献最小，对其进行显著性检验，临界值f(0.05,1,4)＝7.709，f(5)≤f(0.05,1,4)，此方程项不显著，需要剔
除。
[0106]
剔除不显著方程项，新建回归方程：
[0107]
y＝b0 b1×
γ b2×
φ b3×
c b4×
β
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)
[0108]
公式(5)中，b0＝1.73、b1＝
‑
2.42
×
10
‑2、b2＝2.10
×
10
‑2、b3＝4.65
×
10
‑2、b4＝
‑
4.38
×
10
‑2；
[0109]
回归方程显著性检验可参考表5。
[0110]
表5变量分析表
[0111][0112]
样本容量n＝10，检验值ft＝79.61,临界值f(0.05,4,5)＝5.192,ft＞f(0.05,4,5),回归方程显著。
[0113]
f检验值：f(1)＝7.312，f(2)＝21.95，f(3)＝242.6，f(4)＝95.81，各方程项对回归的贡献大小依次为：c、β、φ、γ，方程项γ对回归的贡献最小，对其进行显著性检验，临界值f(0.05,1,5)＝6.608，f(1)＞f(0.05,1,5)，此方程项显著。
[0114]
残差分析可参考表6。
[0115]
表6残差分析表
[0116][0117]
由回归分析结果可知，抗震设防烈度7度以下地区地震作用α为不显著因素，土容重γ、内摩擦角φ、黏聚力c、坡角β均为显著因素，其显著性大小依次为：黏聚力c，坡角β，内摩擦角φ，土容重γ。
[0118]
请参阅图4，基于同一发明构思，本技术实施例还提供一种基于均匀设计的边坡影响因素确定装置200，该装置包括：第一确定模块201、第二确定模块202、构建模块203、处理模块204以及第三确定模块。
[0119]
第一确定模块201，用于确定试验因素；其中，所述试验因素来源于待试验边坡结构模型。
[0120]
第二确定模块202，用于基于所述试验因素以及所述待试验边坡结构确定每个所述试验因素的取值范围。
[0121]
构建模块203，用于基于所述试验因素以及每个所述试验因素的取值范围，构建均匀设计表；其中，所述均匀设计表中包括试验号以及与每个所述试验号中的试验因素的数值。
[0122]
处理模块204，用于通过所述均匀设计表对所述试验因素进行试验，得到每个所述试验号对应的稳定性系数。
[0123]
第三确定模块205，用于基于每个所述试验号对应的稳定性系数，确定出所述试验因素中的边坡影响敏感因素。
[0124]
可选地，构建模块203具体用于确定出试验次数n；获取向量组h＝(h1,h2,...h
m
)；其中，h1,h2,...h
m
均为小于n，且与n的最大公约数为1的数；生成所述均匀设计表的第j列，直至得到所述均匀设计表；其中，第j列的表达式为：u
ij
＝(i
×
h
j
)[modn]；其中，[modn]表示同余运算；若i
×
h
j
超过n，则减去n的一个预设倍数，使落差在[1，n]之中；u
ij
通过递推算法生成；所述均匀设计表记为u
n
(n
m
)；对于确定的n和h所对应的m的个数由欧拉函数e(n)求得；当n为素数时，e(n)＝n
‑
1；当n为素数幂时，即n＝p
l
时，p为素数，l为正整数；当n既不是素数，也不是素数幂时，n表示为不同素数的方幂积；即此时，p1,p2,l p
m
为不同的素数，l1,l2,l l
m
为整数。
[0125]
可选地，处理模块204具体用于通过所述均匀设计表对所述试验因素进行试验；通过简化毕肖普法得到每个所述试验号对应的稳定性系数。
[0126]
可选地，第三确定模块205具体用于采用回退法对所述试验因素进行回归分析，建立回归方程；对所述回归方程进行显著性检验；基于所述显著性检验的结果，确定出所述试验因素中的边坡影响敏感因素。
[0127]
可选地，第三确定模块205还具体用于对所述回归方程进行显著性检验，确定出所述试验因素中的不显著因素；将所述不显著因素剔除，重新构建新的回归方程；对所述新的回归方程进行显著性检验；基于对所述新的回归方程进行显著性检验的结果，确定出所述试验因素中的边坡影响敏感因素。
[0128]
可选地，第一确定模块201具体用于建立边坡数值以构建试验边坡结构模型；其中，所述模型为均质土坡；基于所述待试验边坡结构模型的各项数值，确定出所述试验因素。
[0129]
需要说明的是，由于所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统、装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
[0130]
基于同一发明构思，本技术实施例还提供一种计算机可读存储介质，其上存储有计算机程序，计算机程序在被运行时执行上述实施例中提供的方法。
[0131]
该存储介质可以是计算机能够存取的任何可用介质或者是包含一个或多个可用介质集成的服务器、数据中心等数据存储设备。所述可用介质可以是磁性介质，(例如，软盘、硬盘、磁带)、光介质(例如，dvd)、或者半导体介质(例如固态硬盘solid state disk(ssd))等。
[0132]
在本技术所提供的实施例中，应该理解到，所揭露装置和方法，可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，又例如，多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些通信接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
[0133]
另外，作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
[0134]
再者，在本技术各个实施例中的各功能模块可以集成在一起形成一个独立的部分，也可以是各个模块单独存在，也可以两个或两个以上模块集成形成一个独立的部分。
[0135]
在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。
[0136]
以上所述仅为本技术的实施例而已，并不用于限制本技术的保护范围，对于本领域的技术人员来说，本技术可以有各种更改和变化。凡在本技术的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本技术的保护范围之内。