基于有限元的坡沟系统稳定性分析.docxVIP

基于有限元的坡沟系统稳定性分析 青藏高原有27万多条沟。沟很深，倾斜着破碎。沟密度为3.5km/km2。沟渠占土壤总面积的20%40%。破碎的地形主要由沟蚀引起。沟蚀的发展使沟面变得越来越大，耕地也越来越小。它是世界上水土流失最严重的地区之一。黄土高原的各类沟壑中以沟头前进、沟底下切、沟岸扩张3种形式的沟蚀危害最为严重。沟蚀加剧了面蚀的发展,造成了更多的陡峭临空面,加剧了重力侵蚀。重力侵蚀作为黄土区主要侵蚀类型, 其研究一直处于滞后状态。理论研究方面, 重力侵蚀的微观及宏观力学机理尚不明晰, 计算模型的构建刚起步;试验及观测方面, 至今仍未形成重力侵蚀系统观测和试验研究体系。根据黄河水利委员会西峰、天水、绥德3个水土保持站在典型小流域的调查,黄土高原沟壑区的西峰南小河沟,重力侵蚀面积占流失面积的9.1%,重力侵蚀量占总流失量的57.5%;黄土丘陵沟壑区第Ⅲ副区的天水吕二沟,重力侵蚀面积占流失面积的30.7%,重力侵蚀量占总流失量的68.0%;黄土丘陵沟壑区第Ⅰ副区的绥德韭园沟,重力侵蚀面积占总流失面积的12.9%,重力侵蚀量占总流失量的20.2%。在黄土塬区和丘陵区,沟头前进多以土体崩塌形式进行,沟岸扩张是崩塌与滑坡共同作用的结果。 淤地坝系作为治理水土流失的主要工程措施,在蓄水拦沙、防洪保收等方面起到了重要作用,有机统一了当地致富和治河的关系,深受黄土高原地区群众的喜爱,同时又为治河部门所关注。打坝淤地、蓄水拦沙是流域水土保持综合治理的一项重要措施,目前国内外有关淤地坝的研究主要集中在坝系相对稳定、减水减沙和生态环境效应等方面,而对淤地坝减轻坡沟系统重力侵蚀特别是滑坡侵蚀方面的研究则少有人提及。淤地坝通过抬高沟底侵蚀基准面,提高了沟坡的稳定性,减少了沟坡重力滑坡侵蚀发生的可能性。本研究采用数值模拟方法对随着坝地逐渐淤高沟道坡沟系统的稳定性、滑塌概率以及坡沟系统滑坡侵蚀破坏的部位等方面进行研究,以期为坡沟系统水土保持工程措施及生物措施的配置提供有益的参考,并为评价坡沟系统的稳定性提供一定可靠度的依据。 1 地形地貌、地貌条件 本研究选取的关地沟是王茂沟流域左岸的一条典型支沟。王茂沟是陕北绥德县韭园沟中游左岸的一条试验性治理支沟,是无定河的一条2级支沟,地理坐标为东经110°20′26″—110°22′46″,北纬37°34′13″—37°36′03″,海拔高度940～1 188 m,流域面积5.97 km2,主沟长3.75 km,沟道平均比降为2.7%,沟谷地面积2.97 km2,占流域总面积的46.7%,沟壑密度4.3 km/km2。流域内共有沟道46条,其中0.1 km2的沟道24条,0.1～0.5 km2的沟道18条,0.5～1.0 km2的沟道2条,1.0 km2的沟道2条。地面坡度一般在20°以上,坡度10°的面积仅占2.2%,坡度45°以上的面积占34.7%。地形地貌主要由梁、峁以及分割梁峁的沟谷组成,地形破碎,沟壑纵横,坡陡沟深。梁顶和峁顶面积不大,均略成穹形,坡度在8°～10°,梁、峁坡的坡度一般在20°～35°。梁、峁坡以下为沟谷,沟谷横断面在上游及支沟均呈“V”字型,在下游略呈“U”字型。沟谷坡极陡,一般都在35°以上。按土壤侵蚀区划,属黄土高原丘陵沟壑区第Ⅰ副区。 流域内地质构造比较单纯,表层多被质地匀细、组织疏松的黄绵土覆盖,厚度20～30 m,在梁、峁顶,梁、峁坡上均有分布。其下为红色黄土,厚度50～100 m,多出露于沟谷。再下为基岩,基岩主要为三迭纪的砂页岩,岩层倾角甚小,接近水平,其露头只在干沟中下游沟床及其两侧有所见。 流域属北温带干旱大陆性气候,年平均气温10.2 ℃,无霜期160 d左右,夏季多东南风,春秋多西北风。流域多年平均降水量513.1 mm,降雨量的年际变化率大,年最大降雨量是年最小降水量的3.5倍;降雨的年内分配极不均匀,年内降雨量主要集中在汛期7—9这3个月,汛期降水量占年降水量的73.1%,且多以暴雨形式出现,一次暴雨产沙量为全年产沙量的60%以上,土壤侵蚀以水蚀及重力侵蚀为主,治理前流域年平均侵蚀模数为18 000 t/(km2·a),属剧烈侵蚀区,在黄土丘陵沟壑区具有一定的代表性。 2 分析和计算倾斜系统侵蚀的数值 2.1 破坏概率的估计 对于土质边坡,根据其土体结构、破坏机理和受力状况,可以建立坡体地质条件和环境因素的状态函数:Z=g(X1,X2,…,Xm)其中X1,X2,…,Xm为m个具有一定分布、独立统计的随机变量,假定它们的统计量已知。如果把状态函数定义为安全系数,且随机地从诸随机变量Xi的全体中抽取同分布的变量X′1,X′2,…,X′m,则由上述状态函数可求得安全系数的一个随机样本。如此重复,直到达到预期精度的充分次数N,则可得到N个相对独立的安全系数Z1,