南水北调安阳河段倒虹吸水沙二维数值模拟.pdfVIP

第32卷第9期 人 民 黄 河 Vo1．32．No．9 2010年9月 YELLOW RIVER Sep．．2010 【水利水 电工程 】 南水北调安阳河段倒虹吸水沙二维数值模拟 王勤香 ，潘继顺 ，王鹏涛 (1．黄河水利职业技术学院，河南 开封475001；2．华北水利水电学院，河南 郑州450011) 摘 要 ：南水北调中线总干渠与安阳河交叉的倒虹吸工程，将改变天然河道水流条件，为了使倒虹吸交叉工程布设既安 全又经济，需进行平面二维水沙数值模型设计，对交叉河段进行分析。以河流泥沙动力学为基础。采用非均 匀矩形网格 及交替方向隐式差分逐行求解的ADI有限差分格式进行数学模型设计，依据历史洪水调查对模型进行调试，从洪水水 位、流速、流态及冲淤变化等方面对数学模型计算结果与物理模型试验结果进行比较及验证；采用水边线上切向、法向流 速为零的方法确定水边线，采用水边界全区自动跟踪的方法进行动边界处理，改善了有限差分法不能较好地模拟不规则 岸线的不足 。 关 键 词：河流泥沙动力学；二雏水沙数学模型；ADI有限差分；非均匀网格；动边界处理 中图分类号：TV672 文献标识码 ：A doi：10．3969／j．issn．1000—1379．2010．09．05l 南水北调中线总干渠与安阳河交叉：【程没计采用倒虹吸 形式，倒虹吸工程修建后原河道形成一口门，将改变天然河道 g喏 H(筹+雾) 水沙特点及上、下游河势，因此需对该河段河床演变进行分析 式中： 为水位；Ⅳ为总水深；“、分别为 、Y方向的水深平均流 计算。河床演变研究常用的方法有资料分析法、物理模型法和 速；C为谢才系数；；为水深平均涡黏系数，e=Au。H，A=I．O； 数学模型法，数学模型法对河床变形的理论计算更为灵活、精 “。为摩阻流速，“ = “／c；为水深平均对流项修正系数， 细、快速。因此，笔者采用数学模型法分析安阳河与总于渠交 口 = 1．0。 叉河段洪水演进特点及河床演变，研究倒虹吸_丁程附近不同洪 (2)河床变形模型控制方程 。河床变形方程为 水条件下的洪水位及工程附近局部冲淤变形等，为合理确定口 yy，■ ：OttO(JJ)s—一S。·)， 门位置、倒虹吸水平段尺寸及其埋置深度提供依据。 水流挟沙力方程为 1 模型的基本理论 ，，3 s ‘ ) 南水北调中线总干渠与安阳河交叉河段为宽浅河道 ，河床 床面泥沙起动方程为 以中粗沙为主，局部有卵石滩。该河段流速、水深在垂直方向 U一：Ad H 的变化远小于水平方向的变化 ，流速沿水深分布 比较均匀，水 式中： 为泥沙干密度；to为悬沙沉速；O／为不饱和系数；Js、S。 流流态、水力要素可用沿水深的平均值表示。根据交叉河段水 分别为含沙量和水流挟沙力；z 为河床高程；A为与流速分布、 沙、河床特点，采用平面二维全沙非饱和非耦合的水沙数学模