近年珠江三角洲洪季余水位和潮差的演变特征研究

杨易; 贺涛; Yi YANG; Tao HE

中国农村水利水电 ›› 2021 ›› (7) : 19-23.

作者信息 +

Research on the Evolution of Residual Water Level and Tidal Range of Flood Season in the Pearl River Delta in Recent Years

Author information +

稿件信息 +

摘要

珠江三角洲径潮动力受高强度人类活动的干扰，发生了与自然演变截然不同的变化。为探明近年来珠江三角洲河网内径流与潮汐动力特征的演变，基于珠江三角洲“997”和“177”两场洪水的逐时水位流量资料，采用傅里叶变换方法，获取水位时间序列中的余水位（低频）和潮汐（高频）成分进行对比分析。结果表明：马口、三水断面余水位-流量关系曲线下移，同流量下余水位明显降低，而三角洲下游口门附近余水位略有抬升，使得河道沿程余水位比降变缓，河网内水深增大导致潮汐向上游传播的阻力变小，潮差衰减率也相应减小。近年来三角洲上游水库建设和河网内大范围的航道整治工程导致的河网河道地形下切是径潮动力发生上述变化的主要原因。

Abstract

The river and tidal dynamics in the Pearl River Delta （PRD） are disturbed by strong human activities and changes different from the natural evolution. To explore the evolution of river and tidal dynamic characteristics in the river network of the PRD in recent years， based on the hourly water level and discharge data of two floods in the PRD， “99.7” and “17.7”， this paper adopts the Fourier transform method to obtain the residual water level （low frequency） and tidal （high frequency） components in the water level time series for comparative analysis. The results show that the relation curves of residual water level and discharge at Makou and Sanshui stations move down， with lower residual water level at the same discharge； whereas the residual water levels near the lower mouth of the PRD rise slightly， making the decrease in the residual water level slope； the increase in water depth results in the decrease in the resistance of tidal propagation， which leads to the attenuation rate of tidal range also decreases accordingly. The riverbed down-cutting caused by reservoir construction and waterway regulation is the main reason for the change of river and tidal dynamics of the PRD in recent years.

关键词

珠江三角洲 / 余水位 / 潮差 / 人类活动

Key words

Pearl River Delta / residual water level / tidal range / human activities

基金

国家自然科学基金项目(51879172)

四川省水利厅专项(SKY2019-05)

引用本文

EndNote

Ris (Procite)

Bibtex

导出引用

杨易 , 贺涛. 近年珠江三角洲洪季余水位和潮差的演变特征研究[J].中国农村水利水电, 2021(7): 19-23

Yi YANG , Tao HE. Research on the Evolution of Residual Water Level and Tidal Range of Flood Season in the Pearl River Delta in Recent Years[J].China Rural Water and Hydropower, 2021(7): 19-23

河口三角洲作为连接海洋和陆地的平台，受到流域径流和海洋潮汐两大驱动力的共同作用^［1］。余水位（即潮周期平均水位）和潮差的变化便是河口径潮动力非线性作用的典型结果，而余水位可以用来量化径流动力强度，潮差则是衡量潮汐动力强度的重要指标^［2］。目前国内外学者针对河口区径潮动力非线性耦合作用的研究成果表明，动量守恒方程中余水位梯度主要与非线性摩擦相平衡^［2-4］。因此，研究余水位及其梯度变化成为定量分析径潮动力相互作用的有效切入点，为探讨河口潮波衰减提供了新的思路。

随着社会经济的飞速发展，珠江三角洲在径潮相互作用以及强人类活动的驱动下发生了与自然演变截然不同的变化，研究其河网径潮动力的演变对当地防洪、航运和供水具有重要意义^［5，6］。朱金格等^［7］通过对珠江河网不同年代的水位流量进行数值模拟发现河网区水位流量关系发生明显变化，同流量下水位降低，水面比降变缓；黄伟民^［8］指出，河网顶点马口、三水站分流比在20世纪90年代发生突变，三水分流比增大，北江河网分流增加；张蔚等^［9，10］则通过对河网1955-2005年长序列实测潮位资料进行Mann-Kendall检验发现河网区潮差在20世纪80年代后开始显著增大，涨潮历时也有所增加。虽然近年珠江河网径潮动力演变得到了许多学者的关注，但聚焦到余水位和潮差变化的研究较少。本文以1999年和2017年珠江河网洪季实测水位资料为基础，对比分析河网主要河道余水位和潮差的变化，并探讨了两者的关系，以期为珠江河网径潮动力非线性作用的相关研究提供新的认识。

1 研究区域概况

珠江三角洲位于中国南部沿海，周边从城市经济发达，人类活动频繁。珠江三角洲河网由西江、北江、东江和诸多支流汇聚而成，可大致分为西江河网、北江河网和东江河网三大部分，水流经崖门、虎跳门、鸡蹄门、磨刀门、横门、洪奇门、蕉门、虎门八大口门汇入中国南海，形成“三江汇流、八口出海”的水系格局。

如图1所示，本文关注的研究河段分别为：河网顶端马口站起，经西江干流至磨刀门水道灯笼山站处，全长约136 km，是西江河网的主要排水通道；河网顶端三水站起，经东平水道和顺德水道至洪奇沥水道冯马庙站处，全长约104 km，是北江河网的主要排水通道。所选取的两条河段均为感潮河段，其水位变化既受上游径流的影响，又受外海潮波的作用。河段上游来水季节变化较大，马口、三水站的洪季多年平均流量分别为 10 512和2 989 m³/s，枯季多年平均流量分别为3 099和669 m³/s^［11］。河段下游口门潮汐属不规则半日潮，多年平均潮差0.85~1.62 m，其中磨刀门和洪奇门多年最大平均潮差分别可达2.04和2.94 m^［2，12］。

图1 研究区域简图

Full size|PPT slide

2 数据及方法

2.1 水文数据

本文收集了珠江三角洲两场洪季洪水的实测水位流量资料，分别为1999年7月16日至7月25日的“997”洪水，时长共计192 h；以及2017年7月8日至7月18日的“177”洪水，时长共计240 h。其中“997”期间河网顶端马口和三水站得平均流量之和为31 500 m³/s，“177”期间马口和三水站的平均流量之和为 27 700 m³/s，两者相差不大，“177”期间马口和三水站的平均流量之和相比“997”期间减小约9%。

2.2 傅里叶变换

为了消除潮汐振荡的影响以便获取研究河段上各站点的余水位，采用傅里叶变换的方法对各站点实测逐时潮位数据进行滤波处理。傅里叶变换的公式如下：

Y (k) = \sum_{n = 0}^{N - 1} X (n) e^{- j 2 π n f} k = 0,1, 2, . . ., N - 1

（1）

式中：X（n）为离散的物理量时间序列；Y（k）为离散的傅里叶变换结果，即X（n）的不同频率f的成分。Y（k）通常为一个复杂的数学多项式，它的模代表相应频率f（或周期T）下的振幅。将周期小于（大于）T₁ 的振幅设置为0，则可消除周期小于（大于）T ₁的高频（低频）成分，获得物理量时间序列的低通（高通）滤波X ₁（n）^［13］。

X_{1} (n) = \frac{2}{N} \sum_{n = 0}^{N - 1} Y (k) e^{- j 2 π n f} k = 0,1, 2, . . ., N - 1

（2）

珠江河口一个完整的潮周期大约为一个太阳日的时长25 h，本文以25 h为截断周期（T ₁）对潮位时间序列进行傅里叶变换，其低频成分则为余水位

\bar{Z}

（上划线表示潮周期平均值，下同），高频成分即为潮波。而本文关注的另一个特征参数最大潮差则是一个潮周期内潮波的最大值最小值的差值。

2.3 余水位坡度和余水位方程

余水位

\bar{Z}

沿河道方向的一阶导数

\frac{\partial \bar{Z}}{\partial x}

即为余水位坡度。在不考虑密度梯度影响条件下，一维动量守恒方程表达如下^［14］：

\frac{\partial U}{\partial t} + U \frac{\partial U}{\partial x} + g \frac{\partial Z}{\partial x} + g \frac{U |U|}{K^{2} h^{4 / 3}} = 0

（3）

式中：x表示从口门向上游延伸的距离；t表示时间；U表示断面平均流速；Z表示自由水面高程；h表示水深；g表示重力加速度；K表示曼宁摩擦系数的倒数，上式第四项代表了非线性摩擦项。

可将变化前和变化后的动量守恒方程分别表示为：

\frac{\partial U_{1}}{\partial t} + U_{1} \frac{\partial U_{1}}{\partial x} + g \frac{\partial Z_{1}}{\partial x} + g F_{1} = 0

（4）

\frac{\partial U_{2}}{\partial t} + U_{2} \frac{\partial U_{2}}{\partial x} + g \frac{\partial Z_{2}}{\partial x} + g F_{2} = 0

（5）

式中：用下标1表示变化前的量，下标2表示变化前的量；F则代表动量方程中的非线性摩擦项，即

F = \frac{U |U|}{K^{2} h^{4 / 3}}

。

式（4）与（5）相减可得：

(\frac{\partial U_{1}}{\partial t} - \frac{\partial U_{2}}{\partial t}) + (U_{1} \frac{\partial U_{1}}{\partial x} - U_{2} \frac{\partial U_{2}}{\partial x}) + g (\frac{\partial Z_{1}}{\partial x} + F_{1} - \frac{\partial Z_{2}}{\partial x} + F_{2}) = 0

（6）

假设断面平均流速U是由一个稳定的速度U_r （径流作用）和一个随时间变化的速度U_t （潮流作用）组成：

U = U_{t} - U_{r} = v s i n (ω t) - U_{r}

（7）

式中：v表示潮波的流速振幅；ω为潮波频率。

基于以下假定和近似处理：①断面平均流速U呈周期性变化；②上游下泄流量变化不大，径流作用的变化可忽略；③潮波振幅相对于平均水深为极小量，其沿程变化可忽略^［15］。在一个潮周期内对一维动量守恒方程式（6）进行积分，可得余水位坡度的变化为：

Δ \frac{\partial \bar{Z}}{\partial x} = - Δ \frac{\bar{U |U|}}{K^{2} h^{4 / 3}}

（8）

由式（8）可知，动量守恒方程中余水位梯度的变化与河口河网内非线性摩擦项的变化趋于一致，即余水位坡度的变化很大程度上反映了潮波向上游传播的阻力变化，这为探讨河口潮波衰减提供了新的思路。

3 结果与讨论

3.1 余水位-流量关系

如图2所示，本文统计了“997”和“177”洪水期间河网顶端的马口站和三水站的流量及其对应的水位过程。可以发现，相比于“997”洪水，“177”洪水期间马口站和三水站的余水位-流量关系绳套曲线均发生明显下移，表明马口和三水站在同流量下的余水位明显下降。另外，两场洪水期间马口站同流量下的余水位降幅可达1.5 m左右，大于三水站同流量下的余水位降幅（约0.8 m）。

图2 河网顶部马口、三水站余水位-流量关系

Full size|PPT slide

河网顶端马口、三水余水位-流量关系的变化与其断面地形近年来的显著下切有关。如图3所示，马口断面在1998-2015年间，河床底部高程平均下切约0.83 m，断面面积增大约 3 000 m²，增幅达17.9%；而三水断面在1997-2015年间，河床底部高程平均下切约0.2 m，断面面积增大约1 600 m²，增幅达17.4%。马口、三水断面地形发生剧烈下切，过水面积增大，导致了其同流量下的余水位大幅度下降，而由于马口断面下切幅度比三水断面更大，马口断面余水位-流量关系曲线发生的下移现象也比三水断面更为明显。

图3 河网顶部马口、三水断面地形变化

Full size|PPT slide

3.2 余水位坡度变化

图4是珠江三角洲河网主要河道余水位的沿程变化情况。可以发现，相比与“997”洪水，“177”洪水期间河网顶部余水位大幅度下降，发生这种现象的一部分原因是“177”洪水来流量要略小于“997”洪水（平均流量相差约9%），主要原因还是河网顶部河道发生剧烈下切导致的河床高程下降、过水断面增大。而在河网下游近口门处，余水位因为受到外海潮汐的顶托作用，并未发生明显下降，甚至还由于近年来海平面的升高而略有上升。因此，总体来看，河网内部河道余水位沿程坡度变缓，西江河网主要河道余水位坡度由“997”洪水期间的3.78×10^-5减小为“177”洪水期间的2.37×10^-5，北江河网主要河道余水位坡度则从“997”洪水期间的4.29×10^-5减小为“177”洪水期间的2.69×10^-5。

图4 河网内部主要河道余水位变化

Full size|PPT slide

3.3 河网内潮汐传播阻力及潮差衰减率的变化

如前文所述，余水位坡度的变化与潮波向上游传播所遭遇非线性摩擦作用的变化密切相关。近年珠江河网地形下切导致的余水位坡度变化，也反映出珠江河网地区潮汐的陆向运动时遭遇的摩擦阻力发生变化。一维动量方程中阻力项［gU|U|/（K ² h ^4/3）］与水深h的4/3次方成反比，即水深越大，潮波向上游传播的阻力越小。地形的下切使珠江河网河道断面趋于窄深（图3），河网平均水深整体上呈增大的趋势，导致地形变化后珠江河网内潮汐传播的阻力有所减小。

阻力的减小势必会引起河网内潮汐动力的响应，其中之一就是使得潮差沿程衰减的趋势变缓。图5是珠江河网两条主要河段潮汐由口门向上游传播时最大潮差的衰减情况。整体来看，北江河网潮汐向上游传播时的潮差衰减速度要大于西江河网。“997”洪水期间，西江河网灯笼山-马口的潮差衰减率为0.96×10^-5，北江河网冯马庙-三水段的潮差衰减率为1.63×10^-5；而“177”洪水期间灯笼山-马口段的潮差衰减率减小为0.85×10^-5，北江河网冯马庙-三水的潮差衰减率也有所减小，约为1.59×10^-5。上述变化与“177”洪水期间余水位坡度相比于“997”洪水期间有所减小反映出的阻力减小趋势是一致的。

图5 河网内部主要河道最大潮差变化

Full size|PPT slide

3.4 强人类活动的影响

3.4.1 上游建库

自1949年新中国成立以来，珠江流域水库大坝建设总量超过8 900 座，建成及正在建设的库容超过1 亿m³的大型水库超过33座^［16］。这些水库总库容超过753 亿m³，占珠江平均年径流量的22.4%，可很大程度地调节珠江流域的径流过程^［17］。而这些水库拥有巨大的死库容，可以拦蓄大量泥沙，大幅度削减流域下游珠江三角洲的泥沙来源。图6是珠江三角洲上游高要和石角两站的来水来沙变化情况，可以看到，上游两站年平均流量在1960-2016年间没有发生显著变化，但输沙率在1998年前后发生突变，由1998年前的2 500 kg/s骤降为1998年后的940 kg/s。输沙率大幅度下降，主要原因就是上游建设的大量水库对泥沙的拦蓄作用。而下游三角洲来沙量的减少，引起了河道断面地形冲刷下切，河道过水断面增大，而这种下切现象在三角洲河网上游河段更加显著，进而导致河网余水位-流量关系的改变，余水位坡度的下降。

图6 珠江三角洲上游高要站和石角站的水沙年际变化

Full size|PPT slide

3.4.2 河道采砂和航道整治

与上游建库从而引起下游来水来沙条件改变而导致的河床冲刷不同，采砂和航道整治工程直接作用于河床表面改变河床地形，对局部河段的河床演变产生巨大影响。20世纪80年代初至90年代末，珠江三角洲采砂活动达到高峰，仅1984-1999年间，采砂量就超过8.7 亿m³，相当于珠江三角洲100年左右的河道泥沙自然淤积量^［18］。但1999年后，珠江三角洲地区开始严格控制采砂，采砂地点逐渐转移到下游口门区和外海伶仃洋附近，采砂活动对珠江河网河道的直接改造作用已经很小。

近年来，珠江三角洲开始了大范围的航道整治工程，拟建成“三纵、三横、三线”的高等级航道网^［19］。其中西江下游出海航道、陈村水道、广州港出海航道、东平水道、潭江-莲沙蓉水道、小榄横门出海航道、崖门出海航道、虎跳门水道等河段已经通过航道整治满足规划要求。这些航道整治工程，使得整治河段河床地形进一步下切，引起余水位-流量关系的改变和余水位的降低。

3.4.3 口门围垦

珠江三角洲口门的滩涂围垦改变了沿海地区海岸线形态，直接或间接地影响到珠江河口的潮汐动力。珠江三角洲河口区的滩涂开发利用已有2 000 年的历史，海岸线逐步向海推进。而新中国成立以来，特别是改革开放后，由于城市化进程的快速推进，滩涂围垦的面积也急速增长，海岸线向海推进的速度大幅度提升，仅1966-1996年间围垦面积就达334 km²，海岸线平均向海推进3 km以上^［20］。但随着磨刀门和内伶仃洋的围垦逐渐趋于饱和，近年来围垦的速度又有所放缓。1999后珠江三角洲河口区围垦主要集中在黄茅海东、西岸两侧，围垦面积达30 km²。大部分围垦工程都缩窄了潮汐向上游传播的通道，增大了潮汐传播的阻力，一定程度上削弱了进入河网内部的潮汐动力。

4 结论

（1）近年珠江三角洲顶端马口、三水断面地形显著下切，引起余水位-流量关系曲线下移，同流量下的余水位明显降低，且马口断面降低的程度比三水断面要大。

（2）珠江三角洲上游由于断面地形下切而导致余水位大幅下降，下游余水位受潮汐顶托略有抬升，总体来看，河网马口-灯笼山和三水-冯马庙段的沿程余水位坡度变缓。

（3）地形下切导致珠江三角洲河网水深增大，河网内潮汐传播的阻力减小，潮差衰减率也相应减小，阻力减小的趋势也反映在马口-灯笼山和三水-冯马庙段沿程余水位坡度减缓的变化上。

（4）受近年来水库建设、航道整治等人类活动的影响，珠江三角洲河道地形剧烈下切，是河网内部余水位坡度变缓和潮汐动力增强的主要原因，而口门围垦虽然一定程度上削弱了进入河网内部的潮汐动力，但其削弱强度要小于地形下切所引起的潮汐动力增大的幅度。

参考文献

原文顺序 | 文献年度倒序 | 文中引用次数倒序

1	BIANCHI T S， ALLISON M A. Large-river delta-front estuaries as natural “recorders” of global environmental change ［J］. Proc. Natl. Acad. Sci，2009，106（20）：8 085-8 092. 本文引用 [1]

2	蔡华阳，杨昊，郭晓娟，等. 珠江磨刀门河口径潮动力耦合条件下余水位的多时空尺度分析［J］. 海洋学报，2018，40（7）：55-65. 本文引用 [3]

3	SASSI M， HOITINK A J F. River flow controls on tides and tide-mean water level profiles in a tidal freshwater river： Subtidal Flow in a Tidal River ［J］. Journal of Geophysical Research Atmospheres，2013，118（9）：4 139-4 151.

4	CAI H， SAVENIJE H H G， JIANG C，et al. Analytical approach for determining the mean water level profile in an estuary with substantial fresh water discharge ［J］. Hydrology and Earth System Sciences，2016，12（8）：8 381-8 417. 本文引用 [1]

5	郑国栋，顾立忠，李虎成，等. 珠江三角洲河道地貌变化对网河水情影响研究［J］. 中国农村水利水电， 2010（7）：33-36. 本文引用 [1]

6	易灵，许景锋，尹开霞，等. 西江流域干支流水沙时空分布变化分析［J］. 中国农村水利水电， 2020，447（1）：68-71，76. 本文引用 [1]

7	朱金格，包芸，胡维平，等. 近50年来珠江河网区水动力对地形的响应［J］. 中山大学学报（自然科学版），2010（4）：129-133. 本文引用 [1]

8	黄伟民. 西、北江及思贤滘水文特性变化的分析研究［J］. 人民珠江，2011，32（1）：23-24. 本文引用 [1]

9	张蔚，严以新，郑金海，等. 珠江三角洲年际潮差长期变化趋势［J］. 水科学进展，2010，21（1）：77-83. 本文引用 [1]

10	ZHANG W， RUAN X， ZHANG J，et al. Long-term change in tidal dynamics and its cause in the Pearl River Delta， China ［J］. Geomorphology，2010，120（3）：209-223. 本文引用 [1]

11	王鑫，黄剑威，徐辉荣. 思贤滘近期水文情势分析［J］. 广东水利水电，2015（11）：26-30. 本文引用 [1]

12	LIU F， YUAN L， YANG Q，et al. Hydrological responses to the combined influence of diverse human interventions in the Pearl River Delta， China ［J］. Catena，2014，113：41-55. 本文引用 [1]

13	杨清书，吴超羽. 傅氏变换在短序列确定海平面变化趋势中的应用［J］. 热带地理，1996，16（2）：107-113. 本文引用 [1]

14	DRONKERS J J. Tidal Computations in Rivers and Coastal Waters ［M］. New York：Elsevier，1964：1-518. 本文引用 [1]

15	张萍，谢梅芳，杨昊，等. 潮优型河口动力对水深变化的响应机制研究：以葡萄牙Guadiana河口为例［J］. 热带海洋学报， 2020（1）：1-11. 本文引用 [1]

16	罗宪林. 珠江三角洲网河河床演变［M］. 广州：中山大学出版社，2002. 本文引用 [1]

17	马志鹏，王森，李善综，等. 龙滩水库分类洪水防洪调度规则研究［J］. 中国农村水利水电， 2018（9）：116-120. 本文引用 [1]

18	LUO X， ZENG E， JI R，et al. Effects of in-channel sand excavation on the hydrology of the Pearl River Delta， China［J］. Journal of Hydrology，2007，343（3-4）：230-239. 本文引用 [1]

19	赵学问，车进胜，林超明. 珠江三角洲高等级航道网规划与建设概述［J］. 中国水运（下半月），2011，11（9）：15-16. 本文引用 [1]

20	刘岳峰，韩慕康，邬伦. 珠江三角洲口门区近期演变与围垦远景分析［J］. 地理学报，1998（6）：492-500. 本文引用 [1]

PDF(1705 KB)

854

访问

引用

详细情况

段落导航

收稿日期	出版日期
2020-10-26	2021-07-15
发布日期
2021-08-22

选择文件类型/文献管理软件名称

选择包含的内容