一种水库大坝上游边坡最优坡角确定的新方法与流程

1.本发明涉及安全风险评估技术领域，特别是涉及一种水库大坝上游边坡最优坡角确定的新方法。


背景技术：

2.随着我国经济的不断发展，水利设施的建设越来越多，大型水库的数量在快速的增长，但在过去几十年中，因水库大坝边坡失稳而发生的灾难正在逐年增加，造成了极大的人员伤亡和经济损失。
3.在水库大坝边坡实际工程中，水库大坝边坡坡角的设计通常要考量水库库容与边坡稳定性的影响。在所需水库库容、坝高一定的前提下，上游边坡设计坡角越小，库水位高度越低，反之，上游边坡设计坡角越大，库水位高度越高。在所需不同水库库容下，当上游边坡坡角变化时，由于库水位高度不同，坝体内孔隙水压力分布是不同的，当坝体内孔隙水压力较高时，此时大坝边坡将极易发生失稳。
4.目前结合水库库容与边坡稳定程度对水库大坝上游边坡坡角进行设计的方法还不够全面，水库大坝边坡上游坡角的设计仍然面临挑战，因此亟需一种确定水库大坝上游边坡最优坡角的新方法。


技术实现要素：

5.本发明提供了一种水库大坝上游边坡最优坡角确定的新方法，目的是解决现有技术中存在的问题。
6.为实现上述目的，本发明的技术方案为：
7.一种水库大坝上游边坡最优坡角确定的新方法，包括步骤1、建立水库大坝数学模型，获得不同库容下、库水位的高度随不同上游边坡设计坡角变化曲线；步骤2、利用geo
‑
studio中seep/w模块计算获取水库大坝稳态库水位浸润线位置，绘制水库大坝上游边坡设计坡角值与失效概率关系曲线；步骤3、确定水库上游边坡最优设计坡角。
8.优选的，所述的步骤1中，建立水库大坝数学模型的方法为：设计m个所需水库库容，b1,b2,b3,......b
m
，给出k个大坝上游边坡坡角设计值，记为z1,z2,z3,.......z
k
。
9.优选的，所述的步骤1中，建立水库大坝数学模型后，在某一水库库容b
i
下，改变大坝上游边坡设计坡角，获取z1,z2,z3,.......z
k
分别对应的库水位高度h
i
，利用(z
i
，h
i
)i＝1,2,
…
k，绘制库水位的高度随不同上游边坡设计坡角变化曲线，通过改变库容，重复上述计算直至获得不同库容下，库水位的高度随不同上游边坡设计坡角变化曲线。
10.优选的，所述的步骤2包括如下具体步骤：
11.具体步骤1、在给定某一水库库容b
i
和大坝上游设计坡角z
i
下，利用geo
‑
studio中seep/w模块计算获取水库大坝稳态库水位浸润线位置，将土体黏聚力c、内摩擦角视为随机变量，利用matlab软件产生n个土层参数样本a1,a2,
……
,a
n
，用响应面法计算出a1,a2,
……
,a
n
对应的安全系数fs1,fs2,
……
,fs
n
；若fs
i
<1，则称a
i
为失效样本，如此重复进行
共获取p个不稳定样本，计算出大坝失效概率p
fi
＝p/n；
12.具体步骤2、保持水库库容b
i
不变，改变大坝上游边坡设计坡角，利用具体步骤1获得水库大坝上游设计坡角值z1,z2,z3,.......z
k
分别对应的失效概率p
f1
,p
f2
,p
f3
,......p
fk
，利用(z
i
，p
fi
)i＝1,2,
…
k，绘制水库大坝上游边坡设计坡角值与失效概率关系曲线。
13.优选的，所述的步骤3中，通过改变水库库容，重复具体步骤1和具体步骤2，直至获得在不同水库库容b1,b2,b3,......b
m
下，水库大坝上游边坡设计坡角值与失效概率关系曲线，根据所需水库库容和期望失效概率，查询曲线确定水库上游边坡最优设计坡角。
14.本发明一种水库大坝上游边坡最优坡角确定的新方法的有益效果是：本发明提供了完备的水库大坝上游边坡最优坡角设计方案，经工程验证，该方法准确可靠，可有效解决大坝上游边坡失稳的问题，防止由于边坡失稳导致的灾害发生。
附图说明
15.图1、本发明流程图；
16.图2、大坝模型图(上游设计坡角为35
°
)；
17.图3、大坝模型图(上游设计坡角为30
°
)；
18.图4、大坝模型图(上、下游设计坡角为25
°
)；
19.图5、大坝库水位高度随上游边坡设计坡角变化曲线图；
20.图6、水库上游边坡设计坡角值与失效概率曲线图。
具体实施方式
21.以下所述，是以阶梯递进的方式对本发明的实施方式详细说明，该说明仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
22.本发明的描述中，需要说明的是，术语“上”“下”“左”“右”“顶”“底”“内”“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以及特定的方位构造和操作，因此不能理解为对本发明的限制。
23.请参考图1所示：
24.一种水库大坝上游边坡最优坡角确定的新方法，包括步骤1、建立水库大坝数学模型，获得不同库容下、库水位的高度随不同上游边坡设计坡角变化曲线；步骤2、利用geo
‑
studio中seep/w模块计算获取水库大坝稳态库水位浸润线位置，绘制水库大坝上游边坡设计坡角值与失效概率关系曲线；步骤3、确定水库上游边坡最优设计坡角；
25.所述的步骤1中，建立水库大坝数学模型的方法为：设计m个所需水库库容，b1,b2,b3,......b
m
，给出k个大坝上游边坡坡角设计值，记为z1,z2,z3,.......z
k
；
26.所述的步骤1中，建立水库大坝数学模型后，在某一水库库容b
i
下，改变大坝上游边坡设计坡角，获取z1,z2,z3,.......z
k
分别对应的库水位高度h
i
，利用(z
i
，h
i
)i＝1,2,
…
k，绘制库水位的高度随不同上游边坡设计坡角变化曲线，通过改变库容，重复上述计算直至获得不同库容下，库水位的高度随不同上游边坡设计坡角变化曲线；
27.所述的步骤2包括如下具体步骤：
28.具体步骤1、在给定某一水库库容b
i
和大坝上游设计坡角z
i
下，利用geo
‑
studio中seep/w模块计算获取水库大坝稳态库水位浸润线位置，将土体黏聚力c、内摩擦角视为随机变量，利用matlab软件产生n个土层参数样本a1,a2,
……
,a
n
，用响应面法计算出a1,a2,
……
,a
n
对应的安全系数fs1,fs2,
……
,fs
n
；若fs
i
<1，则称a
i
为失效样本，如此重复进行共获取p个不稳定样本，计算出大坝失效概率p
fi
＝p/n；
29.具体步骤2、保持水库库容b
i
不变，改变大坝上游边坡设计坡角，利用具体步骤1获得水库大坝上游设计坡角值z1,z2,z3,.......z
k
分别对应的失效概率p
f1
,p
f2
,p
f3
,......p
fk
，利用(z
i
，p
fi
)i＝1,2,
…
k，绘制水库大坝上游边坡设计坡角值与失效概率关系曲线；
30.所述的步骤3中，通过改变水库库容，重复具体步骤1和具体步骤2，直至获得在不同水库库容b1,b2,b3,......b
m
下，水库大坝上游边坡设计坡角值与失效概率关系曲线，根据所需水库库容和期望失效概率，查询曲线确定水库上游边坡最优设计坡角。
31.下面结合图2
‑
6进行实例说明。
32.该水库大坝高为28m,坝基厚6m，大坝底是不透水边界，上、下游边坡高度为22m，大坝是由均质土体组成，其土体抗剪强度内摩擦角平均值为28
°
，标准差为8.4，黏聚力c平均值为7kpa，标准差为2.1，容重γ为20kn/m3，图2大坝边坡上游坡角设计值为35
°
，下游坡角为28.22
°
，保持下游坡角(28.22
°
)不变，改变上游坡角，如图3、4所示，上游坡角设计值分别为30
°
、25
°
。
33.采用本发明提出的新技术:
34.利用geo
‑
studio软件中seep/w模块建立水库大坝数值模型，设计三个水库库容分别为：b1＝515m3,b2＝360.5m3,b3＝206m3,保持大坝边坡下游坡角不变，改变大坝上游边坡坡角且设计值分别为：z1＝35
°
，z2＝30
°
，z3＝25
°
，在水库库容b1＝515m3下，改变大坝上游边坡坡角，获取z1＝35
°
，z2＝30
°
，z3＝25
°
分别对应的库水位高度为：h1＝26m，h2＝19.26m，h3＝17.71m；改变水库库容b2＝360.5m3，重复上述步骤直至获得库水位高度分别为：h4＝20m，h5＝15.4m，h6＝14.24m，同理可得水库库容b3＝206m3下，重复上述步骤获得库水位高度分别为：h7＝14m，h8＝10.68m，h9＝9.95m；在所需不同水库库容下，利用(z
i
，h
s
)i＝1,2,3，s＝1,2,3...9,绘制不同水库库容下，库水位的高度随不同上游边坡设计坡角变化曲线，如图5所示；
35.在给定某一水库库容b
i
和大坝上游设计坡角z
i
下，利用geo
‑
studio中seep/w模块计算获取水库大坝稳态库水位浸润线位置，由于考虑到大坝土体参数的不确定性，将其粘聚力c和内摩擦角视为对数正态随机变量，利用matlab软件随机产生105个土体参数样本a1,a2,......,a
100000
,基于响应面法来计算105个土体参数样本对应的n＝105个安全系数，记为fs1,fs2,
……
,fs
10000
。在n＝105个安全系数中，若安全系数fs
i
<1，则称a
i
为失效样本，如此重复进行共获取p个不稳定样本，根据公式p
fi
＝p/105计算出大坝对应的失效概率，保持水库库容不变，改变大坝上游边坡设计坡角，重复上述步骤直至获得水库大坝上游边坡设计坡角值与失效概率关系曲线；改变水库库容，重复上述步骤可以获得不同水库库容下，水库大坝上游边坡设计坡角值与失效概率关系曲线。如图6所示，在水库库容b1＝515m3下，当水库大坝上游边坡设计坡角分别为z1＝35
°
，z2＝30
°
，z3＝25
°
，其对应的失效概率分别为
p
f1
＝0.07％,p
f2
＝1.79％,p
f3
＝1.01％；在水库库容b2＝360.5m3下，当水库大坝上游边坡设计坡角分别为z1＝35
°
，z2＝30
°
，z3＝25
°
，其对应的失效概率分别为p
f1
＝12.18％,p
f2
＝5.04％,p
f3
＝0.2％；在水库库容b3＝206m3下，当水库大坝上游边坡设计坡角分别为z1＝35
°
，z2＝30
°
，z3＝25
°
，其对应的失效概率分别为p
f1
＝28.95％,p
f2
＝4.83％,p
f3
＝0.11％；在不同设计水库库容下，利用(z
i
,p
fi
)i＝1,2,3，分别绘制水库大坝上游设计坡角与边坡失效概率关系曲线，根据所需水库库容和期望失效概率，查询图6关系曲线，确定水库大坝上游边坡的最优设计坡角；譬如当所需库容为360.5m3，期望失效概率为5％，水库大坝上游边坡最优坡角设计值为30.3
°
。
36.因此，对比发现：坡角对水库大坝边坡安全性的影响往往常被忽略，导致对坝体的稳定性考虑不够全面，而本发明通过设计不同的水库库容，改变水库大坝上游边坡坡角，获得大坝不同的库水位高度，考虑土体参数的变异性，来获得不同的水库库容下，大坝上游边坡坡角对应的失效概率，根据所需水库库容和期望失效概率来获得水库大坝上游边坡最优坡角。通过工程实例验证了当水库库容一定时，库水位高度随着大坝上游边坡设计坡角的增大而增大，在水库库容分别为206m3，360.5m3时，大坝失效概率随着上游边坡设计坡角增大而减小，当水库库容为515m3时，大坝失效概率随着上游边坡设计坡角的增大呈现先增大后减小的现象，这充分验证了坡角对边坡稳定性的影响，最终可以根据所需库容和期望失效概率来查询曲线确定水库大坝上游边坡最优设计坡角，通过实例对比分析验证了本发明的有效性。