孔 兰，陈晓宏

(1.中山大学水资源与环境研究中心，广东广州 510275;2.中山大学华南地区水循环和水安全广东普通高校重点实验室，广东广州 510275;3.中水珠江规划勘测设计有限公司，广东广州 510610)

珠江三角洲是广东省的政治、经济、文化中心，也是我国经济最发达的地区之一。水位对珠江口排污、供水安全等有直接影响，水位变化及其影响直接制约着珠江三角洲地区的经济发展。针对珠江口的潮水位变化，国内外不少专家学者均做过研究和探讨。任美锷［1］指出由于地面沉降等原因，过去30年的珠江三角洲相对海平面上升率远大于全球和全国海平面上升率(约1.5 mm/a)。黎子浩［2］认为20 世纪50 年代末至70 年代初，珠江三角洲联围筑闸初期使河道水位迅速雍高，抬高值一般为0.20 ～0.30 m。王廷华［3］认为珠江三角洲网河区各水道水位都不同程度地存在逐年上升的趋势。杨清书等［4］利用20 个验潮站的资料，得出珠江三角洲水位变化总体呈上升趋势的结论。陆永军等［5］对珠江三角洲网河区低水位变化趋势进行了研究。20 世纪80 年代以来，珠江三角洲北江、西江和东江主要河道的水文特征与河床演变趋势先后发生了新的变化，河床的大幅下切引起同流量级的水位明显下降。目前，对珠江三角洲水位长系列变化趋势的研究较多，而对珠江三角洲水位年内变化规律的研究较少。珠江口水位受径流、潮流、河道地形等因素的耦合作用，变化机理复杂［6-11］。河道径流量与上游径流来量、汊道分流比、水利枢纽建设、河口延伸有关;潮汐影响因素包括天文、气象和河口地形等;河道地形变化包括河床自然淤积、人工河床采砂、河道整治工程、桥梁码头建设、岸滩开发等。为了识别珠江口潮水位的年内变化规律，重点选取马口站、澜石站、横门站、灯笼山站和三灶站作为代表站，对比分析珠江口代表站水位年内变化不均匀系数和集中程度，并采用Spearman 秩相关系数法对水位年内变化不均匀系数和集中程度系列进行趋势检验。掌握珠江口水位的年内变化规律，对珠江口地区的防灾减灾和水资源的合理开发利用具有重要的指导意义。

1 数据来源及研究方法

1.1 数据来源

珠江经虎门等8 大口门注入南海。为了识别珠江口潮水位年内变化特征，本研究从水文年鉴、广东省水文局和佛山水务局1973—2004 年的逐月水位数据中收集了珠江三角洲马口站、澜石站、横门站、灯笼山站和三灶站5 个主要代表水文站(水位站)(各个站点具体分布信息如图1)的实测水文数据。

1.2 研究方法

1.2.1 不均匀性

图1 珠江口主要水文(水位)站分布

采用水位年内变化不均匀系数Cu来分析水位的年内变化，Cu值越大，年内各月水位差异越大，水位年内变化越不均匀，水位调控难度也越大。由于篇幅所限，具体计算方法见文献［12-13］。

1.2.2 集中程度

目前利用集中程度研究河流径流年内分配变化规律的成果较多，笔者引用集中度Cn和集中期D表达水位在年内各时段的集中程度以及最大水位出现的时段。集中度反映水位年内集中程度，集中期反映全年水位集中的重心所出现的月份。集中度和集中期的具体计算步骤见文献［14-18］。

1.2.3 Spearman 秩相关系数法

为了研究潮水位的年内变化指标的时序变化规律，本研究采用Spearman 秩相关系数法。Spearman秩相关系数是一种等级相关系数，由于应用中无须考虑有关变量的分布类型，事实上，它比大家熟知的Pearson 相关系数适用范围更广［19］。如果秩相关系数Zsp＞0，表示序列有上升趋势;Zsp＜0，表明序列有下降趋势;Zsp=0，表明序列没有变化趋势。显著性概率水平α 对应的临界值Z0.05/2=1.96。

2 计算与分析

由表1 可以看出，珠江口潮水位的年内变化在不同年代存在明显的时空差异。位于珠江口深处的马口站在20 世纪70 年代和20 世纪90 年代水位的年内变化大于80 年代，其他站在20 世纪70 年代水位的年内变化最大，20 世纪80 年代以来水位的年内变化趋于均匀;马口站和澜石站水位年内变化不均匀系数和集中度最大，近口门处的横门站、灯笼山站和三灶站水位年内变化不均匀系数和集中度较小;马口站和澜石站水位集中期最小，主要在7 月，横门站和灯笼山站较大，分别在8 月和9 月，三灶站最大，主要在10 月，说明年内水位集中期由珠江口深处向河口处有后推(或增加)趋势。初步分析马口站和澜石站水位年内变化受径流变化的影响较大，而近口门处的横门站、灯笼山站和三灶站水位年内变化受潮汐和海平面变化的影响较大［20-22］。珠江口潮水位变化机理复杂，有待于进一步深入研究。

表1 不同年代水位年内变化指标值

图2 和图3 能直观反映珠江口水位年内变化特征随时间的变化趋势:马口站年内变化不均匀系数和集中度呈上升趋势;澜石站、横门站、灯笼山站和三灶站水位年内变化不均匀系数和集中度的时序变化相似，均有下降趋势;各代表站的集中期都随时间呈波动上升趋势，马口站和澜石站集中期主要在7—8 月份波动，横门站和灯笼山站集中期主要在8—9月份波动，而三灶站集中期主要在10—11 月份波动。

表2 显示:马口站年内变化不均匀系数和集中度系列的Spearman 检验统计量为正，表明其水位年内变化不均匀程度存在较显著的上升趋势;澜石站、横门站、灯笼山站和三灶站水位年内变化不均匀系数和集中度系列的Spearman 检验统计量都为负，表明其水位年内变化不均匀程度均有下降趋势，其中横门站和三灶站年内变化不均匀系数系列的Spearman 检验统计量的绝对值均大于1.96，通过了0.05 显著性水平检验，说明近海口门处潮水位年内变化不均匀程度均有显著下降趋势;所选代表站的年内水位集中期都具有不显著的上升趋势，说明时间上是后推趋势。

表2 Spearman 检验统计量

3 结 论

图2 各站点水位年内变化指标及其趋势

图3 各站点水位集中期变化趋势

a. 珠江口潮水位的年内变化在不同年代存在明显的时空差异。马口站在20 世纪70 年代和90年代水位的年内变化大于80 年代，其他站20 世纪80 年代以来水位的年内变化趋于均匀;珠江口水位年内变化不均匀系数和集中度由珠江口深处向近河口处有减小趋势;年内水位集中期由珠江口深处向河口处有后推(或增加)趋势。

b. 珠江口水位年内变化趋势特征明显。珠江口深处的马口站年内变化不均匀系数和集中度呈上升趋势，其他站水位年内变化不均匀系数和集中度均有下降趋势，其中近海口门处横门站和三灶站的潮水位年内变化不均匀程度均有显著下降趋势;各代表站的集中期都随时间呈波动上升趋势，说明时间上是后推(或增加)趋势。

［1］任美锷.黄河长江珠江三角洲近30 年海平面上升趋势及2030 年上升量预测［J］. 地理学报，1993，48(5):385-392. (REN Meie. Relativer sea level rise in Huanghe，Changjiang and Zhujiang (Yellow，Yangtze and Pearl River)delta over the last 30 years and predication for the next 40 years (2030)［J］. Acta Geographica Sinica，1993，48(5):385-392.(in Chinese))

［2］黎子浩.珠江三角洲联围筑闸对水流及河床演变的影响［J］. 热带地理，1985(2):99-107. (LI Zihao. The effects of united embankments on the flow shift and channel changes in the Zhujiang delta［J］. Tropical Geography，1985(2):99-107.(in Chinese))

［3］王廷华.珠江三角洲水位变化趋势及其影响分析［J］.人民珠江，1998(6):38-41.(WANG Tinghua.Analysis of the trend of water level change in Pearl River delta and its effects［J］.Pearl River，1998(6):38-41.(in Chinese))

［4］杨清书，罗章仁，张修杰.珠江三角洲近几十年水位变化趋势研究［J］.热带海洋，1998，17(2):9-14.(YANG Qingshu，LUO Zhangren，ZHANG Xiujie. Secular trend of water level change in Zhujiang River delta in recent decades［J］. Tropic Oceanology，1998，17(2):9-14. (in Chinese))

［5］陆永军，贾良文，莫思平，等.珠江三角洲网河低水位变化［M］.北京:中国水利出版社，2008.

［6］汤立群，陈洁，申锦瑜，等.风暴作用下河口动力过程数值模拟研究进展［J］. 水利水电科技进展，2009，29(5):78-83. (TANG Liqun，CHEN Jie，SHENG Jinyu，et al.Advances in simulation of dynamic estuarial processes reacting to storm surges in estuaries［J］. Advances in Science and Technology of Water Resources，2009，29(5):78-83.(in Chinese))

［7］朱文谨，李瑞杰，干岳良，等.感潮河段底床阻力特性分析［J］.水利水电科技进展，2009，29(6):25-28. (ZHU Wenjin，LI Ruijie，GAN Yueliang，et al. Analysis on the riverbed resistance characteristic of tidal river［J］.Advances in Science and Technology of Water Resources，2009，29(6):25-28.(in Chinese))

［8］谭超，杨清书，刘秋海，等.1950 年代和1990 年代西北江网河区水动力特征变化［J］. 水利水电科技进展，2009，29(S1):28-32. (TAN Chao，YANG Qingshu，LIU Qiuhai，et al. Variation characteristics of hydrodynamic condition of West River and North River network in 1950s and 1990s［J］. Advances in Science and Technology of Water Resources，2009，29(S1):28-32.(in Chinese))

［9］杨涛，陈喜，杨红卫，等.基于线性矩法的珠江三角洲区域洪水频率分析［J］. 河海大学学报:自然科学版，2009，37(6):615-619. (YANG Tao，CHEN Xi，YANG Hongwei，et al.Regional flood frequency analysis in Pearl River delta region based on L-moments approach［J］.Journal of Hohai University:Natural Sciences，2009，37(6):615-619.(in Chinese))

［10］陈晓宏，张蕾，时钟.珠江三角洲河网区水位特征空间变异性研究［J］. 水利学报，2004(10):36-42. (CHEN Xiaohong，ZHANG Lei，SHI Zhong. Study on spatial variability of water levels in river net of Pearl River delta［J］.Journal of Hydraulic Engineering，2004(10):36-42.(in Chinese))

［11］ZHANG Q，XU C Y，CHEN Y Q，et al.Spatial assessment of hydrologic alteration across the Pearl River delta，China，and possible underlying causes ［J］. Hydrology Process，2009(23):1565-1574.

［12］李艳，陈晓宏，张鹏飞.北江流域径流序列年内分配特征及其趋势分析［J］.中山大学学报:自然科学版，2007(5):113-116. (LI Yan，CHEN Xiaohong，ZHANG Pengfei.Research on annual distribution and its trends of runoff in Beijiang River，Guangdong［J］.Acta Scientiarum Naturalium Universitatis Sunyatsen，2007(5):113-116.(in Chinese))

［13］王金星，张建云，李岩，等. 近50 年来中国六大流域径流年内分配变化趋势［J］. 水科学进展，2008(5):656-660. (WANG Jinxing，ZHANG Jianyun，LI Yan，et al.Variation trends of runoffs seasonal distribution of the six larger basins in China over the past 50 years［J］.Advances in Water Science，2008 (5):656-660. (in Chinese))

［14］汤奇成，程天文，李秀云.中国河川月径流的集中度和集中期的初步研究［J］. 地理学报，1982，37(4):383-393. (TANG Qicheng，CHENG Tianwen LI Xiuyun.Preliminary study on the degree and time of concentration of monthly runoff of Chinese streams ［J］. Acta Geographica Sinica，1982，37 (4 ):383-393. (in Chinese))

［15］杨远东.河川径流年内分配的计算方法［J］.地理学报，1984 (2):218-227. (YANG Yuandong. Calculating method of annual distribution of stream flows［J］. Acta Geographica Sinica，1984(2):218-227.(in Chinese))

［16］陆志华，夏自强，于岚岚，等.松花江干流中游段径流年内分配变化规律［J］. 河海大学学报:自然科学版，2012，40 (1):63-69. (LU Zhihua，XIA Ziqiang，YU Lanlan，et al. Changes of annual runoff distribution in middle reaches of Songhua River［J］. Journal of Hohai University:Natural Sciences，2012，40(1):63-69. (in Chinese))

［17］刘贤赵，李嘉竹，宿庆，等.基于集中度与集中期的径流年内分配研究［J］. 地理科学，2007(6):791-795. (LIU Xianzhao，LI Jiazhu，SU Qing，et al. Interannual runoff distribution based on degree and time of concentration for rivers［J］. Scientia Geographica Sinica，2007(6):791-795.(in Chinese))

［18］王兆礼，陈晓宏，杨涛.东江流域径流序列年内分配特征研究［J］. 人民黄河，2011，33(2):37-39. WANG Zhaoli，CHEN Xiaohong，YANG Tao. Study on characters of interannual runoff distribution in Dongjiang River area［J］.Yellow River，2011，33(2):37-39.(in Chinese))

［19］王钦德，冯锁江. 关于Spearman 系数计算公式的论证及应用［J］. 山西农业大学学报，1993，13(1):30-33.(WANG Qinde，FENG Suojiang. The demonstration and application about calculating formula of Spearman codfficient［J］. Journal of Shanxi Agricultural University，1993，13(1):30-33.(in Chinese))

［20］黄镇国.广东海平面变化及其影响与对策［M］. 广州:广东科技出版社.2000.

［21］孔兰，陈晓宏，陈栋为，等. 珠江三角洲水位演变分析［J］.生态环境学报，2010，19(11):2642-2646. (KONG Lan，CHEN Xiaohong，CHEN Dongwei，et al. Analysis on evolution of water level in Pearl River delta，China［J］.Ecology and Environmental Sciences，2010，19(11):2642-2646.(in Chinese))

［22］孔兰，陈晓宏，张强，等.海平面上升对珠江口水位影响的分析［J］. 生态环境学报，2010，19(2):390-393.(KONG Lan，CHEN Xiaohong，ZHANG Qiang，et al.Impacts of rising sea level on water level changes along the Pearl River Estuary， China ［J］. Ecology and Environmental Sciences，2010，19 (2):390-393. (in Chinese))