三峡工程运用初期石首河弯河势演变三维数值模拟 3d numerical simulation of fluvial processes in the shishou bend during the early filling of the three gorges reservoir.pdfVIP

VoL 第2l卷第1期 水科学进展 21，No．1 2010年1月 ADVANCESINWATERSCIENCE Jan．，2010 三峡工程运用初期石首河弯河势演变三维数值模拟 假冬冬，邵学军，王 虹，肖 毅，周 刚 (清华大学水沙科学与水利水电工程国家重点实验室，北京100084) 摘要：采用考虑河岸坍塌的三维水沙数值模型，对三峡工程运用初期石首河弯的河势演变进行了模拟分析。数值模 型中采用基于非正交网格的局部网格可动技术，将传统水沙模型与二元结构河岸坍塌力学模型相结合。利用三峡 水库蓄水前后2002～2(104年石首河段冲淤演变资料对模型进行了验证，模拟结果与实测结果吻合较好；模型不仅 能模拟出河道的垂向冲淤过程，而且还可模拟出由河岸坍塌所引起的河道横向摆动过程。计算分析了三峡工程运 行至2016年末河段的河势演变情况，结果可为河段的河道整治提供参考。 关键词：三峡工程；石首河弯；河势演变；河岸坍塌；三维；水沙数值模拟 147．5 中图分类号：TV 文献标志码：A 文章编号：1001-6791(2010)01-0043-07 三峡工程运行后，出库泥沙量将大为减小，下游河道将发生长时期、长距离的冲刷，而荆江河段更是 “首当其冲”。对于下荆江，河岸抗冲性较弱、坍塌严重，由此引起剧烈的河势变化，并对防洪、航运及人 们的生产生活产生一定影响，对其河势变化的研究具有重要实际意义。陆永军等…采用二维水沙数值模型， 对上荆江瓦口子至马家咀河段进行了模拟研究；谢作涛等[21则对宜昌至杨家脑河段的冲淤演变进行了二维 水沙数值模拟。随着计算技术的进步，三维水沙数学模型得到了一定的发展和应用[3剖。但对三峡水库运用 后，下游干流河道河势演变的三维水沙数值模拟研究则较少。尤其对下荆江等蜿蜒型河道来说，其水沙运动 具有显著三维性。为此本文采用考虑河岸坍塌的三维水沙数值模型【8】，对三峡工程运用初期石首河弯的河 势演变进行了模拟分析。模型中采用基于非正交网格的局部网格可动技术，将三维水沙数学模型与二元结构 河岸坍塌力学模型有效结合。模型不仅能模拟出河道的垂向冲淤过程，而且还可模拟出由河岸坍塌引起的横 向摆动过程。利用三峡水库蓄水前后2002～2004年石首河段的冲淤演变资料对模型进行了验证，并计算分 析了三峡工程运行至2016年石首河段的河势演变情况，可为河段的河道整治提供参考。 1 考虑河岸坍塌的三维水沙数值模型 1．1三维水沙模型 三维水沙模型可参见文献[8]，其中本文采用非均匀沙计算模式。 1．2崩岸模型 河床沿垂向冲淤变形以及边岸的淤长过程通过传统的泥沙冲淤计算即可模拟出，但边岸的后退过程，是 无法通过沿垂向冲淤变形体现的，需对崩岸过程进行模拟才能实现旧o。 下荆江河岸由二元结构组成，在垂向存在明显的分层结构，下层是较厚的非粘性泥沙，而上层则为较薄 的粘性土层。通常下部的非粘性岸坡坡度接近泥沙颗粒的摩擦角，而上部粘性土体则几乎为直立的。崩岸为 收稿日期：2009-02-23 2006BAB05802) 作者简介：假冬冬(1982一)，男，广西桂林人，博士研究生，主要从事水力学及河流动力学研究。 E—lIlail cn jdd05@mails．tsinghua．edIL 万方数据 水科学进展 第2l卷 坍塌形式，其冲蚀机理可分作两个相互作用的过程阳]，一个是水流的直接淘刷下层非粘性土层的过程，另 一个就是上层粘性土层在重力作用下的失稳过程。当河岸下层的非粘性土层冲蚀宽度达到极限状态后，上部 挂空的粘性土层将会受拉崩裂，其受力状态如图I所示。 +1)以及断裂面的拉力；根据力矩平衡，处于坍塌I临界状态时有如下受力关系： A形