钱塘江河口洪水特性及动床数值预报模型论文.docVIP

　　钱塘江河口洪水特性及动床数值预报模型论文 摘要：在分析钱塘江河口水文、地形及洪水位实测资料的基础上，得到了钱塘江河口段潮汐、径流、河床三者的非线性关系对河口段的洪水位有显著影响的规律，以此为基础建立了洪水预报的一维动床数值预报模型。模型复演了钱塘江河口一次典型的洪水位过程，并利用动床模型和定床模型的两种计算模式对各站洪水位进行了比较，两种模式计算得到的沿程洪水位可差0.3～2.0m，但动床模拟与实测过程十分吻合,表明在象钱塘江河口这样大冲大淤的洪水位预报中考虑动床冲淤过程是必要的。 关键词：钱塘江河口 洪水特性 动床模型 0前 言 钱塘江河口是一个冲积性河口(图1)，洪冲潮淤，年内基本平衡。每年汛期前4月份.freelm，泥沙起动流速0.5～0.7m/s，极易起动，致使河床在洪水期变化较快，一般历时3～5天的洪水可将前期7个月淤积的泥沙搬运至河口段的下游。这表明在钱塘江河口洪水预报中必须考虑河床的冲淤过程。 图1钱塘江河口示意图 Fig.1 Sketch map of Qiantang Estuary 洪水预报的主要任务是根据河段上游的已知流量过程，预测下游沿程重要城镇的流量、流速及水位过程。洪水系长波，其运动规律可用圣维南方程组来描述。钱塘江河口段河床由细粉沙组成，易冲易淤，河床惯性较小，特别是在洪水期较短的时间内，河床冲刷剧烈，同一洪峰流量的洪水位可相差1m之多，这主要是由于在钱塘江河口段潮汐、径流、河床地形三者的非线性关系对河口段的洪水位有显著的影响。几十年来，科学工作者在分析钱塘江河口大量水文、潮汐及地形的实测资料的基础上，得到了影响河口洪水位的几种因素，并根据河口段潮流、泥沙运动与河床变形理论进行了数值预报模型的研究，取得了不少的研究成果12。但限于当时计算技术条件，对洪水前期的江道地形采用经验修正的方法得到，再用经验合轴相关曲线作闸口站洪水位预报，它的优点是方便快速，缺点是不能动态预报洪水位沿时空的变化过程。本文在前人工作的基础上，建立了一维动床数值预报模型，进行了钱塘江河口的一次典型洪水过程的数值预报，其速度和精度可满足预报要求，为实时预报奠定了基础。 1钱塘江河口洪水位变化特性 1.1钱塘江河口洪水特性钱塘江河口七堡以上河段最高水位受梅汛期洪水控制（仅个别年份由台风暴潮引起）。由于钱塘江河口具有潮大流急，河床冲淤幅度大、速度快的特点，使得该河段的洪水位变化规律既不同于无潮河流，也不同于冲淤变化较小的潮汐河流。具体地说，它除了与汛期洪水流量大小直接有关外，还受江道冲淤面貌、下游潮汐、河口下段尖山河湾主槽走向、长度及人类活动（包括新安江水库的兴建及治江围垦）等诸多因素控制。一般地，闸口洪水位（Z）可由闸口洪峰流量（Q）、闸口至盐官江道容积（V）、下游澉浦高潮位（ZZ）三者的非线性回归关系确定,由此建立经验预报合轴相关图2。 1.2 泥沙输移对洪水位影响钱塘江河口泥沙输移对沿程洪水位影响比较敏感。实测资料分析表明，在枯水大潮期间，盐官以下河段落潮含沙量大于涨潮含沙量，而盐官以上河段涨潮含沙量大于落潮，净泥沙输移是指向上游，从而造成钱塘江河口自上而下的淤积，导致该河段年内潮淤洪冲，年际间表现为丰水年冲刷，枯水年淤积的格局。泥沙输移的最终结果使闸口至盐官河段河床发生冲淤变化，进而影响行洪面积。该河段的冲淤面貌可用钱塘江河口沙坎高程及闸口至盐官河段7m下容积等指标表征，它们均可用闸口站低水位集中反映。从表1的几次洪水资料可以看出江道冲淤面貌对洪水位的影响起决定作用。比较1955年6月与1968年7月两次闸口洪水位，由于汛前江道面貌不同，1955年汛前江道冲刷，面积、容积均很大，汛前闸口断面低水位仅3.87m，沙坎高程2.1m，而1968年汛前遭遇枯水而淤积，江道面积、容积较小，仅2.68亿m3，闸口低水位高达6.14m，沙坎高程4.0m，尽管闸口洪水流量相差2倍多，但其最高水位非常接近，分别为9.07m和9.01m，出现小流量高水位的现象主要在于江道淤积所致。而1968年7月与1955年4月另一次洪水相比闸口洪峰流量相近，分别为12000m3/s和11000m3/s，但由于江道面貌不同，两者的洪水位差达2m之多。 表1 实测典型洪水比较 Tab.1 parison of measured typical floods in the Qiantang river 年.月.日 洪峰 流量/m3.s-1 闸口 洪水位/m 汛前闸口低水位/m 河口沙坎 高程/m 澉浦 高潮位/m 闸口至盐官 容积/亿m3 1955.6.22 31000 9.07 3.87 2.10 6.61 9.31 1968.7.10 12000 9.01 6.14 4.0 6.60 2.68 1955.4.17 11000 7.00 4.47 3