河网一维水流模型在赣江二级航道整治中的应用.pdfVIP

河网一维水流模型在赣江二级航道整治中的应用 马中伟，陈界仁，汪迎春 河海大学水利水电工程学院，江苏南京 (210098) E-mail：zhwma@163.com 摘 要：针对河网地区复杂的地形条件及防洪评价中存在的问题，建立河网一维水流数学 模型，对赣江南昌～湖口段水流运动以及二级航道整治工程对防洪的影响进行了数值模 拟。结果表明：工程建成后，由于丁坝的影响，水位壅高，在整治线范围内，流速明显增 加，工程的整治效果明显。 关键词：数学模型；河网；赣江航道 在珠江三角洲、荆江－洞庭湖区、太湖以及鄱阳湖流域等地区，众多河流与湖泊交织在 一起，形成复杂河网。这些河网中各河段纵横交错，相互干扰，在实测资料不足或难以取得 的情况下，对各单一河段的计算相当困难[1] 。本文在已有研究成果的基础上，对赣江东河、 西河河网一维水流数学模型进行初步研究探讨，对赣江二级航道整治工程建设前后的河道水 流流速、水位、流态的变化情况进行准确模拟和计算，并由此分析工程建设对防洪的影响， 为该工程提供技术支持。 1. 工程河段概况 赣江干流在南昌市裘家洲头首次分为东西两河后又分四汊入湖：东河于礁矶头分汊为中 支、东支，中支经南新、大口湖在朱港入湖，河长43km ；东支是东河的主流，经叶楼、滁 槎，在三江口与抚河、信江两支交汇后入湖，河长56km；西河在樵舍分汊为北支和西支， 北支干流经蒋埠在朱港同中支汇合，长28km，仅在中洪水分流，枯水过流少；西支经樵舍、 昌邑于吴城入湖，是赣江入湖进长江的主要通道。吴城至湖口属鄱阳湖区，长约81km，枯 水平均比降 0.047‰，吴城以下20km 的褚溪河口是鄱阳湖五河来水的总汇口。河道形势见 图1。 图1 河道形势图 -1- 2. 模型建立及相关问题处理 2.1 基本方程 河网水流运动基本方程一般采用 Saint-Venant 方程组[2]，包括水流连续方程和动力方 程。该方程组广泛应用于河道、河网水流运动的模拟。 Q A q ∂ ∂ L x t + 水流连续方程 ∂ ∂ （1） ∂Q ∂ Q 2 ∂Z Q Q + (α ) + gA + gA 2 q L ν x 动力方程 ∂t ∂x A ∂x K （2 ） 式中:x，t—距离（米）和时间（秒），为自变量； A—过水面积（米2 ）； Q—断面流量(米3/秒) ； Z—水位（米）； α—动量修正系数； K—流量模数； qL—旁侧入流（米2/秒）；入流为正，出流为负； ν ν 0 x—入流沿水流方向的速度（米/秒），若旁侧入流垂直于主流，则： x 。 2.2 方程的离散 基本方程的离散格式是Preissmann 加权四点线性隐式差分格式[3]，见图2 。函数及其时 间导数和空间导数的离散公式是 图2 差分格式示意图