钱塘江河口七堡段氯度时空变化及与水文的关系

2019-03-26李若华夏冬梅姚凯华

水利水电科技进展 2019年2期

李若华, 夏冬梅, 姚凯华

(1.浙江省水利河口研究院,浙江杭州 310020; 2.河海大学环境学院,江苏南京 210098)

钱塘江河口潮强流急,涌潮汹涌,山潮水比值仅为0.01,是典型的强涌潮河口,在极端枯水条件下,盐水可上溯至离杭州湾口207 km处的闻家堰以上,其盐水入侵长度和强度远大于一般潮汐河口[1]。盐水入侵对河口的水动力、水质和生态具有重要影响[2],会影响河口地区的取用水[3],同时盐度作为保守的示踪剂,还经常用于研究河口的水循环和混合过程[4]。因此,盐水入侵研究具有重要学术意义和应用价值,一直是河口研究的热点[5-9]。

从20世纪60年代开始,不少学者研究了钱塘江河口的盐水入侵规律。毛汉礼等[10]分析了钱塘江河口沿程月平均氯度,认为其符合Aron和Stommel混合长度理论。韩曾萃等[11]首次采用一维水流盐度数值模型研究了钱塘江河口盐水入侵规律,预测了钱塘江水体含氯度,在此基础上建立了钱塘江河口盐水入侵中长期动床预报数值模型,并将该模型应用于富春江水库每年枯季泄水流量调度[4],同时还分析了新安江建库、河口治江缩窄对河口盐水入侵的影响[12]。史英标等[13]采用钱塘江河口段长历时实测水文、氯度资料,分析了强潮作用下盐水入侵的时空变化特征。潘存鸿等[1,14]通过实测资料及数学模型分析,研究了径流、潮汐、涌潮等对钱塘江河口盐水入侵的影响。

钱塘江河口是一个强冲淤的河口,其盐水入侵强度及盐度在河口的时空分布不仅受上游径流和外海潮汐控制,还深受河床冲淤的影响,这也是钱塘江河口的盐水入侵规律较其它河口更加复杂的原因所在。随着现场观测资料和研究成果的累积,钱塘江河口的盐水入侵规律研究也逐渐深入。钱塘江河口七堡河段是杭州市最重要的水源地,本文基于多年的水文、地形、氯度等实测资料及现有研究成果,结合盐水入侵与水文之间的关系特性,分析钱塘江河口七堡河段盐水入侵的一些新规律。

表3 2009年9月大潮定点实测氯度

注:分层系数为底表层盐度差与垂线平均盐度之比;横向差异系数为深槽平均和浅滩平均之差与深槽平均之比。

1 氯度空间分布特征

1.1 纵向分布特征

钱塘江河口长期观测氯度的水文站有闸口、七堡、仓前、盐官、澉浦5个站(图1),每日采用滴定法观测日最大、最小氯度(氯离子的质量浓度),已积累了近50年的资料。本文采用尖山河湾治江缩窄后的2000—2015年年均氯度资料,统计钱塘江河口沿程各站多年氯度特征值如表1所示,这16年中包括了丰、平、枯水年,具有较好的代表性。表2为1953年以来钱塘江河口各站潮汐特征值。

图1 钱塘江河口形势图及站点位置

表1 钱塘江河口2000—2015年沿程各站氯度特征值

注:变异系数为标准差与均值之比;采用的河床容积为闸口—盐官河段吴淞高程7 m以下的年均河床容积,下同。

表2 钱塘江河口各站潮汐特征值

由表1、表2可见,澉浦以上河段潮差快速减小,涨潮历时大幅缩短,至七堡后多年平均潮差由澉浦的5.64 m减小至0.79 m,平均涨潮历时也由澉浦的5.47 h缩短至1.42 h,说明潮汐在河口上溯过程中潮汐能量快速减弱;与此相对应的是年平均氯度也由澉浦的4.29 g/L减小至0.11 g/L,表明盐水在上溯过程中,潮量逐渐减小,受下泄淡水径流的稀释作用逐渐增强,氯度逐渐降低;同时还可以看到,越往上游,氯度的变异系数越大,其中闸口—仓前河段变异系数较大,七堡达到最大,表明闸口—仓前河段是径流、潮流势力彼消此长的区域,淡水、盐水掺混强烈,径流大则可有效压制盐水。杭州的多个取水口布置在闸口—七堡河段,该特点为杭州在秋季大潮期调度上游水库加大下泄流量抗咸提供了可行性[15]。

1.2 横向及垂向分布特征

钱塘江河口虽属强潮河口,但不同河段的混合程度仍有所差别,尤其是仓前上游的杭州段,处于河口段的上游,山潮水比值逐渐增大,涌潮逐渐消失,潮流速减小,水体紊动掺混作用较中、下游弱,加之在复式河槽中浅滩与深槽的水深相差较多,还有弯道环流的影响等,导致横向流速分布不均匀,致使氯度在江道断面横向分布上不一致。一般认为,因深槽水流流速大,盐度随潮流运动的速度较快,在涨憩前后,深槽的盐度大于浅滩,且深槽与浅滩水深之比越大越显著[16]。

表3为2009年9月大潮定点实测氯度情况,由表3可见,钱塘江河口在一个潮中的最大氯度,上、下游均呈现出深槽大于浅滩的现象,总体趋势是下游差异小,上游差异大,但横向差异最大的闸口断面的横向差异系数也小于0.3。同时可以看到,七堡—澉浦河段氯度分层差异也较小,分层系数均在0.16以下,闸口氯度已低于0.25 g/L的饮用水标准,可看作为淡水。因此,钱塘江河口氯度的差异基本在0.3以内,属强混合型河口[17]。

2 氯度时间变化特征

2.1 年际变化

以仓前站为例,由表1可知,其年均氯度最大可达0.85 g/L,最小仅0.16 g/L,同时也可看到影响盐水入侵的径流量及河床容积变化也很大。绘制仓前站年均氯度与年均径流量、潮差、河床容积的关系如图2～4所示。

图2 仓前站年均氯度与年均径流量

图3 仓前站年均氯度与年均潮差

图4 仓前站年均氯度与年均河床容积

由图2～4可见,丰水年(2002年、2010年、2012年和2015年)时河道冲刷,河床容积和潮差大,但因径流量大压制了盐水入侵,仓前站氯度较小;连续枯水年(2005—2009年)后,河床淤高容积变小,潮差小,潮流上溯能力弱,虽然此时径流量小,但仓前站氯度也较小;径流从丰转枯的第一年(2003年),虽径流量大幅减小,因河床容积和潮差还较大,氯度也较高。同时也可看到,潮差与河床容积的变化趋势基本是一致的,河床容积大,则潮流上溯能力强,导致潮差大,由此可见,潮差大小可在一定程度上反映河床的冲淤程度。

2.2 季节变化

因径流和潮汐的季节性变化,导致钱塘江河口氯度在年内变化十分明显。统计2000—2015年各月平均日最大氯度,绘制各月氯度与月均径流、潮差、河床容积(钱塘江河口每年4月、7月、11月进行3次地形测验)的关系见图5～7。从图5～7可见,3—7月为丰水期,径流较大,且潮汐相对较弱,氯度较小;8—11月径流大幅减小,平均不足丰水期的一半,对盐度的稀释能力大幅减弱,且河床经过汛期的强烈冲刷,河床容积较4月大幅增加,另外秋季是大潮期,外海潮汐较强,故河口段潮流上溯能力大幅增强,导致氯度较大;12月至次年2月虽为枯水期,径流量较小,但河口段的河床经过秋季大潮期4个月的回淤,容积大幅减小,潮流上溯能力弱,河口进潮量将大幅减少,氯度也较小。

图5 仓前站多年月均日最大氯度与月均径流量

图6 仓前站多年月均日最大氯度与月均潮差

图7 仓前站多年月均日最大氯度与河床容积

2.3 月相变化

2007年10—11月,在七堡断面进行了连续半个月逐时的氯度、潮位、潮流观测,其结果如图8所示。由图8可见,小潮期潮位低,潮差小,随着小潮期向大潮期过渡,潮差逐渐增大,高、低潮位均逐渐抬升,10月23—29日低潮位抬升1.3 m,高潮位抬升2.1 m,表明潮差同步增大了0.8 m。小潮时氯度很低,水体为淡水,进入大潮期后,氯度快速升高,10月26—30日氯度波动剧烈,最高氯度达到1.3 g/L。大潮期过后,氯度又恢复至低值。由此可见,七堡河段的盐水入侵主要集中在大潮期的这几天,其他时间基本不受盐水入侵的影响。

图8 2007年七堡断面氯度和潮位月相变化

3 大潮期氯度变化与潮位、流速的关系

图9 大潮期七堡断面氯度与潮位的关系

图10 大潮期七堡断面氯度与流速的关系

图11 大潮期七堡断面水深与潮差的关系

由上述分析可知,钱塘江河口盐水入侵主要发生在大潮期。图9～11展示了2007年10月25—29日大潮期间七堡断面的潮位、垂线平均流速(设定涨潮流速为正,落潮流速为负)、垂线平均氯度的变化。由图9～10可见,由于涌潮的存在[18],涨潮初期潮位急剧抬高,同时水流迅速从落潮流转为涨潮流,并在半小时内涨潮流速达到极值,大潮期间涨潮最大流速可以达到1.8 m/s,而落潮最大流速仅0.8 m/s,由此可见,该河段涨潮历时虽然小于落潮历时,但涨潮期的水动力远大于落潮期的水动力。同时可看到,因盐分需由下游随水流带至上游,故氯度抬升滞后于潮位和流速,高潮位和涨潮最大流速过后氯度才出现最大值,而泥沙含量基本与潮流同步达到最大值[19],这是因为泥沙不仅可被涨潮流从下游带至上游,本地河床的泥沙还会被涌潮掀起悬浮,这也是钱塘江河口盐度和泥沙输运特征不同的主要机制之一。由图11可见,在前5个潮(约60 h内)期间,该河段处于蓄潮阶段,潮平均水深逐渐增加,潮位、潮差、涨落潮流速逐渐增大,氯度也逐渐增大。在第5个潮,高潮位、潮差、涨潮最大流速、氯度均达到最大值。之后低潮位基本保持不变,高潮位呈现一高一低交替的现象,相应的潮平均水深、潮差和涨潮最大流速也呈现一大一小交替的现象。在前60 h内,每个潮周期的氯度峰值逐渐升高,60 h后氯度出现最大值,为1.28 g/L,之后与高潮位(或潮差)对应,每个潮周期的氯度峰值呈现一高一低交替的现象,但氯度峰值逐渐减小。整个大潮期间每个潮周期的最小氯度基本相同,约0.05 g/L,说明该河段在落潮后期已为淡水。