王建平 李伟华 孙红兵 郭华
张家口市京新高速公路管理处
摘 要: 为了解决公路边坡稳定性分析的不确定性、随机性和模糊性问题,采用模糊数学的理论建立二级模糊综合评价模型,采用降半梯型分布函数确定连续型变量的隶属度,并利用层次分析法确定各因素指标权重,采用专家现场评定的方法确定层次分析法的判断矩阵,可以减少公路边坡稳定性分析过程中主观因素造成的误差;结合实例分析,验证了方法的可行性。
关键词: 道路工程;模糊综合评价;层次分析法;边坡稳定性;
作者简介: 王建平(1969—),男,河北阳原人,本科,正高级工程师,主要从事道路工程建设管理与研究工作;
基金: 河北省交通运输厅科技计划(编号:T-2012129);
近年来,我国公路的建设逐渐由平原区到地形复杂山区建设过渡,公路等级也更多的变为高等级公路建设,与此同时为保证道路施工以及运营阶段的安全,边坡稳定性问题显得尤为重要[1],公路路基边坡的稳定性评价与计算过程中,涉及到许多定量和定性因素边坡稳定性分析至关重要。我国已有较多的研究者将模糊数学方法用于边坡稳定性的研究,洪海春[2]在考虑影响边坡稳定性的10个主控因素基础上建立二级模糊综合评价模型并用Visual FoxPro编程语言将其模型可视化。张春宇[3]运用模糊综合评价法并结合湖南高速工程实例总结了4项20类边坡稳定指标,计算得到公路边坡稳定性的预测结果与实际结果吻合。胡田飞[4]结合熵权法确定边坡因素指标的权重进行模糊综合评价。但是,边坡稳定性问题是一个多因素耦合的一种不确定性的复杂系统,传统的定性分析定量计算的方法难以精确求解,本文采用层次分析法与模糊综合评价相结合的方法对公路路基边坡进行系统分析评价,可以有效解决公路边坡稳定性分析的不确定性、随机性和模糊性问题,得到相对科学合理的结果。
1模糊综合评价模型
1.1 评价指标
根据模糊数学中的隶属度理论,综合评价法将定性评价转变成为定量的评价,适用于解决各种不确定性问题。
首先确定评语集 V :本文设定 V =( v 1 v 2 v 3 v 4)表征路基边坡的稳定性状态,即:(稳定,较稳定,较不稳定,不稳定)。
评判的因素集 U :采用二级模糊综合评价,分工程地质特征、气象水文特征、地形地貌特征和其他因素4大项21个指标建立评价因素集,如图1所示。

图1 边坡二级模糊评价指标 *载下**原图
1.2 隶属度的确定
正确地构造隶属度函数是用好模糊控制的关键,本文采用专家评定法来确定离散型变量的隶属度,离散型变量的隶属度指标如表1所示。
表1 定性指标特征及隶属度取值 导出到EXCEL
|
指标 |
特征 |
隶属度μ |
|||
|
Ⅰ |
Ⅱ |
Ⅲ |
Ⅳ |
||
|
岩体类型a1 |
坚硬岩石 |
0.75 |
0.25 |
0 |
0 |
|
胶结好的半坚硬岩石 |
0.25 |
0.65 |
0.1 |
0 |
|
|
胶结差的半坚硬岩石 |
0 |
0.2 |
0.6 |
0.2 |
|
|
软弱岩类 |
0 |
0 |
0.25 |
0.75 |
|
|
土体类型a2 |
砂性土 |
0.7 |
0.3 |
0 |
0 |
|
黏性土 |
0.4 |
0.5 |
0.1 |
0 |
|
|
粉性土 |
0 |
0.2 |
0.4 |
0.4 |
|
|
特殊土 |
0 |
0 |
0.25 |
0.75 |
|
|
岩体结构a5 |
均质结构 |
0.75 |
0.25 |
0 |
0 |
|
块状结构 |
0.2 |
0.5 |
0.2 |
0.1 |
|
|
层状结构 |
0.1 |
0.2 |
0.5 |
0.2 |
|
|
松散结构 |
0 |
0.05 |
0.25 |
0.7 |
|
|
节理面状况a6 |
不发育 |
0.8 |
0.2 |
0 |
0 |
|
平行路线发育 |
0.1 |
0.4 |
0.4 |
0.1 |
|
|
斜交路线发育 |
0.05 |
0.1 |
0.6 |
0.25 |
|
|
垂直路线发育 |
0 |
0.1 |
0.2 |
0.7 |
|
|
降雨强度b3 |
小雨 |
0.9 |
0.1 |
0 |
0 |
|
中雨 |
0.2 |
0.7 |
0.1 |
0 |
|
|
大雨 |
0 |
0.1 |
0.6 |
0.3 |
|
|
暴雨及特大暴雨 |
0 |
0 |
0.3 |
0.7 |
|
|
斜坡类型c3(软弱结构面与斜坡坡面关系) |
逆向坡 |
0.65 |
0.2 |
0.1 |
0.05 |
|
横向坡 |
0.25 |
0.5 |
0.2 |
0.05 |
|
|
斜交坡 |
0.05 |
0.2 |
0.5 |
0.25 |
|
|
顺向坡 |
0.05 |
0.1 |
0.2 |
0.65 |
|
|
边坡类型c4 |
凸形坡 |
0.75 |
0.2 |
0.1 |
0 |
|
直形坡 |
0.25 |
0.65 |
0.1 |
0 |
|
|
凹形坡 |
0 |
0.1 |
0.65 |
0.25 |
|
|
S形坡 |
0 |
0.1 |
0.2 |
0.7 |
|
|
已有地质现象c5 |
极稳定 |
0.8 |
0.2 |
0 |
0 |
|
稳定 |
0.1 |
0.8 |
0.1 |
0 |
|
|
较不稳定 |
0.05 |
0.1 |
0.6 |
0.25 |
|
|
不稳定 |
0.05 |
0.1 |
0.2 |
0.65 |
|
|
边坡防护d2 |
强 |
0.8 |
0.2 |
0 |
0 |
|
中 |
0.1 |
0.7 |
0.2 |
0 |
|
|
弱 |
0 |
0.2 |
0.7 |
0.1 |
|
|
无 |
0 |
0 |
0.2 |
0.8 |
|
对于连续性变量的隶属度,采用“降半梯型”分布函数计算[5], j 分别为1、2、3、4,隶属度函数公式如公式(1)~公式(4)所示。
rij=⎧⎩⎨⎪⎪1 x≤x1x2−xx2−x1 x1<x<x20 x>x2 (1)rij=⎧⎩⎨⎪⎪⎪⎪0,x<x1或x>3−x1−xx2−x1 x1<x≤x2x3−xx3−x2 x2<x≤x3 (2)rij=⎧⎩⎨⎪⎪⎪⎪0,x<x2或x>4−x2−xx3−x2 x2<x≤x3x4−xx4−x3 x3<x≤x4 (3)rij=⎧⎩⎨⎪⎪1 x≤x1x2−xx2−x1 x1<x<x20 x>x2 (4)rij={1 x≤x1x2-xx2-x1 x1<x<x20 x>x2 (1)rij={0,x<x1或x>3-x1-xx2-x1 x1<x≤x2x3-xx3-x2 x2<x≤x3 (2)rij={0,x<x2或x>4-x2-xx3-x2 x2<x≤x3x4-xx4-x3 x3<x≤x4 (3)rij={1 x≤x1x2-xx2-x1 x1<x<x20 x>x2 (4)
式中: x 为指标变量值; x 1、 x 2、 x 3、 x 4为分别是各定量指标所对应的4级分级界限值,见表2。
表2 定量指标及分级界限 导出到EXCEL
|
指标 |
分级 |
|||
|
Ⅰ |
Ⅱ |
Ⅲ |
Ⅳ |
|
|
岩体软化系数a3 |
>0.9 |
0.7~0.9 |
0.5~0.7 |
<0.5 |
|
土体粘聚力a4/kPa |
>45 |
45~20 |
20~5 |
<5 |
|
初始地应力a7/MPa |
<2 |
2~10 |
10~20 |
>20 |
|
年平均降雨量b1 /mm |
<500 |
500~1 000 |
1 000~1 500 |
>1 500 |
|
雨水冲刷深度b2/m |
<0.1 |
0.1~0.3 |
0.3~0.5 |
>0.5 |
|
地下水位b4/m |
<20 |
20~50 |
50~80 |
>80 |
|
岩土体饱水率b5/% |
<0.25 |
0.25~0.5 |
0.5~0.75 |
>0.75 |
|
坡度c1/(°) |
<15 |
15~30 |
30~50 |
>50 |
|
坡高c2/m |
<50 |
50~100 |
100~150 |
>150 |
|
植被覆盖率d1/% |
>30 |
30~15 |
15~5 |
<5 |
|
岩土风化程度d3/% |
<5 |
5~10 |
10~30 |
>30 |
|
地震烈度d4/度 |
<3 |
3~5 |
5~8 |
>8 |
1.3 综合评价
综合评估模型由因素判断矩阵的模糊运算和指标确定的权重矩阵来构建 M = w · R ,评判矩阵 R 为
R=[rij]m×n=⎡⎣⎢⎢⎢⎢r11r21⋮rm1r12r22⋮rm2⋯⋯⋮⋯r1nr2n⋮rmn⎤⎦⎥⎥⎥⎥ (5)R=[rij]m×n=[r11r12⋯r1nr21r22⋯r2n⋮⋮⋮⋮rm1rm2⋯rmn] (5)
式中: w 为各因素的权重分配向量; R 为由各个单因素评价组合成的评判矩阵。
做模糊变化后得到综合隶属度 B ,便可根据最大隶属度原则进行边坡稳定性的评判,得到具体边坡稳定等级。
2权重的确定
本文采用层次分析法确定各指标的权重,根据各项指标对于路基边坡稳定性的影响的重要程度,通过因素之间的两两比较重要性构造判断矩阵,进行9标度相对重要性判断。计算得出各指标权重,如表3所示。
表3 各因素指标权重值 导出到EXCEL
|
指标 |
一级权重 |
二级权重 |
指标 |
一级权重 |
二级权重 |
||
|
a1 |
0.14 |
||||||
|
a2 |
0.17 |
c1 |
0.12 |
||||
|
a3 |
0.15 |
c2 |
0.28 |
||||
|
工程地质特征 |
a4 |
0.14 |
0.41 |
地形地貌特征 |
c3 |
0.2 |
0.13 |
|
a5 |
0.12 |
c4 |
0.15 |
||||
|
a6 |
0.22 |
c5 |
0.25 |
||||
|
a7 |
0.06 |
||||||
|
b1 |
0.14 |
d1 |
0.21 |
||||
|
b2 |
0.19 |
d2 |
0.32 |
||||
|
气象水文特征 |
b3 |
0.40 |
0.38 |
其他因素特征 |
d3 |
0.37 |
0.08 |
|
b4 |
0.15 |
d4 |
0.10 |
||||
|
b5 |
0.12 |
||||||
3工程实例分析
京新高速公路张家口段地处寒温带大陆性季风气候,年平均降雨量350~500 mm, 地面标高800~1 030 m, 地面坡度3°~22°,地震影响烈度对高速公路的最大影响烈度为Ⅶ,地层上部2.1~6.0 m为中等湿陷性黄土,局部覆盖有2~25 m左右厚冲洪积碎石、粘性土等,坡段有少量低矮灌木,岩石风化作用强烈,裂隙率为1.1%,坡高54~60 m。以京新高速公路中张家口路段K23+500~K23+780边坡实际情况为例,构造出模糊关系矩阵如下。
A 工程地质特征模糊关系矩阵、对应权重矩阵
wa =[0.14 0.17 0.15 0.14 0.12 0.22 0.06] (6)
Ra=⎡⎣⎢⎢⎢⎢0.250.650.100.000.000.000.250.750.650.350.000.000.200.800.000.000.100.200.500.200.050.100.600.250.3750.6250.000.00⎤⎦⎥⎥⎥⎥T (7)Ra=[0.250.000.650.200.100.050.3750.650.000.350.800.200.100.6250.100.250.000.000.500.600.000.000.750.000.000.200.250.00]Τ (7)
B 气象水文特征模糊关系矩阵、对应权重矩阵
wb =[0.14 0.19 0.40 0.15 0.12] (8)
Rb=⎡⎣⎢⎢⎢⎢1.00.00.00.01.00.00.00.00.00.00.30.71.00.00.00.00.00.20.80.0⎤⎦⎥⎥⎥⎥T (9)Rb=[1.01.00.01.00.00.00.00.00.00.20.00.00.30.00.80.00.00.70.00.0]Τ (9)
C 地形地貌特征模糊关系矩阵、对应权重矩阵:
wc =[0.12 0.28 0.20 0.15 0.25] (10)
Rc=⎡⎣⎢⎢⎢⎢0.000.450.550.000.860.140.000.000.050.100.200.650.250.650.100.000.050.100.650.05⎤⎦⎥⎥⎥⎥T(11)Rc=[0.000.860.050.250.050.450.140.100.650.100.550.000.200.100.650.000.000.650.000.05]Τ(11)
D 其他因素特征模糊关系矩阵、对应权重矩阵:
wd =[0.21 0.32 0.37 0.10] (12)
Rd=⎡⎣⎢⎢⎢⎢0.000.000.500.500.000.000.200.800.000.001.000.000.000.670.330.00⎤⎦⎥⎥⎥⎥T (13)Rd=[0.000.000.000.000.000.000.000.670.500.201.000.330.500.800.000.00]Τ (13)
将上述模糊关系矩阵与相应一级权重进行模糊评价,为体现权数作用以及各因素综合作用程度采用算子 M (·,⊕)得到二级模糊关系矩阵:
S =[ SaSbScSd ] T (14)
Sa =[0.206 0 0.339 0 0.248 5 0.206 5] (15)
Sb =[0.480 0.024 0.216 0.280] (16)
Sc =[0.300 8 0.235 7 0.283 5 0.192 5] (17)
Sd =[0.000 0.067 0.572 0.361] (18)
然后模糊运算二级权重,归一化处理其结果,得到最后的综合的模糊评价矩阵 DM =[0.270 5 0.183 4 0.300 5 0.245 6]。
由最大隶属度原则可知该段边坡段处于“较不稳定”状态,需要增加边坡防护措施。这与现场施工时的实际结果相吻合。
4结 论
(1)采用二级模糊综合评价,工程地质特征、气象水文特征、地形地貌特征和其他因素4大项21个指标,建立了公路边坡稳定性评价指标体系。
(2)采用专家评定法和“降半梯型”分布函数计算分别确定了离散型变量和连续型变量的隶属度,并应用层次分析法原理确定了不同评价指标的权重。
(3)结合工程实例,基于层次分析-模糊评价法对公路边坡稳定性进行了计算分析,可以减少公路边坡稳定性分析过程中主观因素造成的误差,提高边坡分析的准确性与科学性。
参考文献
[1] 周创兵,李典庆.暴雨诱发滑坡致灾机理与减灾方法研究进展[J].地球科学进展,2009,(5):477-487.
[2] 洪海春,徐卫亚,叶明亮.基于模糊综合评判的边坡稳定性分析[J].河海大学学报,2005,33(5):557-562.
[3] 张春宇,高燕希,蒋明锋,等.二级模糊综合评判在公路边坡稳定性分析中的应用[J].公路工程,2008,33(5):110-113.
[4] 胡田飞,朱本珍.基于熵权法和层次分析法的复杂边坡稳定性模糊综合评价方法[J].铁道建筑,2013,(12):69-73.
[5] 夏卜敬.基于模糊综合评价和神经网络对边坡稳定性的分析研究[D].武汉:武汉科技大学,2008.

声明: 我们尊重原创,也注重分享。有部分内容来自互联网,版权归原作者所有,仅供学习参考之用,禁止用于商业用途,如无意中侵犯了哪个媒体、公司、企业或个人等的知识产权,请联系删除(邮箱:glyhzx@126.com),另本头条号推送内容仅代表作者观点,与头条号运营方无关,内容真伪请读者自行鉴别,本头条号不承担任何责任。