中国水文区划及水文分区方法综述
2024-03-11孙周亮刘艳丽刘翠善王国庆
孙周亮,王 婕,刘艳丽,3,4,陈 鑫,刘翠善,3,4,王国庆,3,4
(1.南京水利科学研究院 水灾害防御全国重点实验室,江苏 南京 210029; 2.水资源工程与调度全国重点实验室(武汉大学),湖北 武汉 430072; 3.水利部应对气候变化研究中心,江苏 南京 210029; 4.长江保护与绿色发展研究院,江苏 南京 210098)
分区是地理学研究中的常用方法之一,是认识地理要素空间分布差异的重要手段[1]。地理学中的分区通常按照不同的标准划分,如 中国按地势可划分为三级阶梯;按自然地理可划分为东部季风区、西北干旱半干旱区和青藏高寒区;按气候可划分为热带季风气候区、亚热带季风气候区、温带季风气候区、温带大陆性气候区、高原山地气候区;中国陆地生态系统可划分为生态调节、产品提供与人居保障三类生态功能一级区,这些自然地理分区为社会经济发展提供了理论支撑[2-3]。水文区划是自然地理区划的重要组成部分[4],指根据水文指标对特定的研究区域开展空间上的划分与分类,常用于水文站网建设、水资源规划、水文资料及参数移植等,也是解决水文缺资料问题的重要方法[5-6]。在实际工作中,对“水文区划”和“水文分区”(或“流域分类”)常未作区分,为避免混淆,文中将“水文区划”界定为一项综合性研究工作,包含资料收集整理、水文分区、成果合理性检查等系列内容,将“水文分区”界定为水文区划工作的核心内容。本文首先梳理了中国水文区划的发展历程,总结水文分区方法的发展和应用现状,并展望未来水文分区的研究重点方向,为水文站网建设及水文资料移用提供方法参考。
1 中国水文区划发展历程
中国水文区划工作开始于1949 年后,主要用于水资源规划和水文站网布设规划[6-7]。水文区划研究是认识水文规律的重要手段,更是水文研究全面发展的必要环节[8]。由于当时的资料积累和技术限制,所得水文区划成果主要用于规划与研究工作,公开出版和发表的资料较少。根据水文区划成果的水平及完整性总体上可分为以下三个阶段:
(1)第一阶段:全国性水文区划初探时期。全国性的水文区划成果最早见于《中国自然区划草案》中的“中国水文区划草案”[9]和《中国水文区划(初稿)》[10]。其中《中国自然区划草案》在苏联欧亚大陆水文区划的基础上进行改进,根据流域、水流形态、冰情、含沙量等特征将全国分为3 级9 区。第一级分区将全国划分为内流区和外流区。外流区第二级按冬季河流结冰情况分为冰冻区和不冻区,冰冻区第三级按含沙量划分为高含沙区和低含沙区,不冻区第三级按相对流量划分为高流量区和低流量区。内流区第二级按水流形态划分为西藏和蒙新,在西藏区按潴水形态划分第三级区,在蒙新区按有无河流划分第三级区。作为第一次全国性的水文区划,《中国自然区划草案》的划分具有十分明显的地理特色,即划分原则主要依据河流水系的自然属性,对水文特征的考虑较少。
《中国水文区划(初稿)》主要根据水文特征和水利条件将中国划分为13 个水文区、46 个水文地带和89 个水文省。第一级划分指标为径流深,第二级划分指标为河流的季节变化特征,第三级划分指标为水利条件。在这三级区域划分体系之外,还将全国按照径流补给情况划分为雨水补给区、雨水融水补给区和融水雨水补给区,称为“0 级”区域。第二次的区划相对第一次更完善,各分区之间的水文特征异同点更加突出,尤其是径流丰枯特征及年内变化特征的空间差异性得到了较好展现。
(2)第二阶段:区域性水文区划细化时期。《中国水文区划(初稿)》的发布吸引国内多个单位开展了区域性水文区划研究工作。随着时间推移和水利工作推进,已有水文区划成果的不足逐渐凸显,如分区特征的区域代表性不足、各项分区研究成果之间水平参差不齐难以统一等[6],已无法满足中国社会经济发展需求。中国科学院地理研究所专门开展了水文区划的更新研究工作[11],将全国按径流量分为丰水、多水、平水、少水和干涸5 个径流带,再根据径流年内分配和动态特征分为56 个二级区,并对径流极值特征进行专门分区,总结了径流年内分配、季节变化、丰枯变化等水文特征分区规律。这次的水文区划成果除了吸纳前两次的经验外,相比以前的成果应用了更多的观测资料(7 万余站),在反映水文规律的区域性特征方面更加完善。其中第一级分区揭示了全国水量的空间差异,第二级分区反映了河流水量的年内分配。这主要源于本次区划采用了综合性、相似性与差异性、成因分析等基本原则。可见,中国水文区划工作一直在持续发展中,初期以开展全国性水文区划为主,但分区结果较为粗糙,后期主要以省级行政区为尺度开展水文分区研究,分区指标及方法也在不断完善中,目前还缺乏针对全国的较精细的研究成果。
(3)第三阶段:区域性水文区划改进时期。水文区划的综合性特征决定了其难以同时满足多方面应用需求。比如水文站网规划中希望分区内的水文过程及暴雨洪水特性具有一致性,水利规划中则希望分区内的水资源开发条件具有较好的空间均匀性,而在水文参数移植方面则希望分区内具有水文响应的空间相似性。在一次水文区划中难以同时兼顾到多方面的需求,如基于气候、地形划分的水文区划在暴雨特性方面差异较大[12]。因而没有一种分区是普适性的,应要根据研究需求进行更新和修正。水文工作者在前期水文区划成果的基础上又开展了大量研究,冷荣梅[13]为了合理规划四川省水文站网,以水文气象特征为主、以自然地理特征为辅,将四川省划分为2 个水文大区、7 个水文亚区、15 个水文分区;赵宝君[14]先按干旱指数、再按地形、地貌、植被和水文特征将内蒙古分为7 个大区、37 个分区,改进了单指标水文分区;徐磊等[15]为了研究天山北坡暴雨洪水空间规律,以年最大流量一阶样本矩为聚类因子,将天山北坡分为两个分区;苏乃友[16]为了研究无资料地区的设计年径流,根据年径流变差系数CVR与流域气候特点和下垫面条件将山西省划分为4 个分区,并拟合了各分区的CVR地理综合数学模型;李珂等[17]以产、汇流物理过程为依据选取降雨径流和下垫面因子将秦岭北麓划分为3 个水文相似子区,为水文相似分区内的资料移用问题提供了科学依据。
总体而言,水文区划的标准或原则是多样化的,即使在同一区域也会有不同的分区方案,这取决于研究者的分区目的,具体来说包括水资源开发利用、水文站网布设与优化、工程设计参数移用、水文资料及参数移用、生态环境分级保护等。针对不同分区目的采用的分区方法一般无明显区别,但采用的分区指标则有所侧重,如生态保护研究以生态健康类指标为主要指标,水文无资料问题则以水文过程参数的关键表征和影响要素为主要指标。综上所述,在应用已有水文区划成果时,应首先根据成果的研究方法论证其适用性,重点了解分区原则和指标,以及各分区间的差异性、分区内的相似性。鉴于水文区划的多样性属性和水文学的多元化发展趋势特征,未来将出现更多的分区方案以满足不同的研究需求。
2 水文分区方法研究进展
2.1 水文分区指标
水文分区方法体系主要包括两部分:确定分区指标、确定分类方法(表1)。分区指标(因子)是决定分区方案的关键因素。如前所述,不同的研究者会根据不同的需求选择分区指标。一般而言,所用的指标不外乎气候、水文、下垫面三大类,这些是决定水文特征的主要因素。早期的水文分区主要用于水文站网规划和水资源规划等[18],重点关注分区的流域特征相似性,因而分区指标主要有降雨量、径流量、蒸发量、土地利用类型等,在青藏高原地区还需要考虑气候和环境的垂直差异[19]。Hughes 等[20]采用年径流、月径流、低流量、高流量等16 个指标,发现低流量无法用于水文分区。随着缺资料地区水文研究的广泛开展,水文模型参数逐渐成为主要的分区指标。胡凤彬等[21]较早地将新安江水文模型参数用于水文分区,并构建了各分区的参数地理关系;颜梅春等[22-23]结合下垫面信息与土壤地形指数(STI)进行水文分区,并通过在各分区内分别率定SCS 模型的最大土壤蓄水能力参数CN 模拟降雨径流,从而验证了分区的合理性。
表1 常用水文分区指标及分区分类方法Tab.1 Frequently employed indicators and classification techniques in hydrological zoning
除了前述常用的流域特征及水文变量外,还发展了一些新的水文分区指标。Pagliero 等[36]研究发现基于偏最小二乘回归相似度指标的聚类分区方法较好;薄岩等[37]构建了基于水量平衡的水文响应模型,发现基于水文响应的分区方法优于基于水文特征相似和基于流域结构相似的分区方法。在区域洪水研究中,常将洪水频率分布曲线参数相同的区域称为“均匀水文分区”,以表明该分区的洪水特性相同[38]。徐磊等[15]将洪水频率一阶样本矩作为聚类因子,基于模糊聚类法的分区结果通过了一致性检验和均匀性检验。径流特征是水文分区中常用的指标,Bouma 等[39]研究认为对于流域特性分类而言,土壤特征比流量特征更合适,原因在于土壤特征属于永久性特征,而流量特征具有时间动态性。Trancoso 等[40]研究发现长期主导径流特性的驱动因素按重要性排序如下:干旱指数、光合作用有效辐射吸收比、饱和导水率、土壤深度、最大坡度、木本植被覆盖度。可见,不同研究中采用的水文分区指标类别基本一致(即气候、水文、下垫面),区别在于具体的因子,产生差异的原因可能与其研究目的相关。
2.2 水文分区方法
水文分区的分类在本质上属于模式识别的问题[6],在早期的水文区划研究中主要依靠研究人员的经验,即根据等值线图[41]或水文地理因子空间分布[42]进行人工判别。后期主要应用人工智能技术进行分区识别,包括聚类方法和神经网络方法。聚类方法通常包括主成分聚类、模糊聚类和主成分模糊聚类等。自组织特征映射人工神经网络(SOFM 网络)具有模拟人脑学习功能、自动识别分区个数的优势,能够提高水文分区的合理性[43]。
在水文分区方法比较方面,房晶等[24]对黑河的水文分区研究表明:模糊聚类法和集对分析法结果相近、但前者更合理,主要体现在模糊聚类法对流域差异更敏感;姬海娟等[19]首先应用主成分分析剔除了因子之间的相关性信息,然后用K均值聚类分析方法对雅江542 个子流域进行归类,分区结果与实际水文特征分布与前期研究成果吻合,基于斯米尔诺夫检验方法的结果也表明了分区的合理性;Agarwal 等[33]针对K均值聚类对初始分布敏感的问题,提出小波-自组织映射耦合的聚类方法,结果表明该方法比K均值聚类法更优,可以作为小波-K均值聚类法的替代方法。
分区方法的适用性与可靠性是分区结果的重要指标。Tasker[44]较早关注到了水文分区方法的稳定性比较,提出了数据分割方法,该方法将大量子流域及其数据分为“估算集”和“预测集”(每个流域只会属于1 个集,两个数据集的子流域所覆盖的空间区域基本重叠且统计特征基本相同),通过应用不同数据集分别划分水文分区以检验聚类和回归两类分区方法的效果。在水文频率分区研究中,常用S值法和H值法检验水文分区的均匀性,陈永勤等[45]采用H值法检验了枯水频率分区,从假定的相似区中剔除了异常子区域;魏兆珍等[46]研究了土地利用和流域尺度对水文分区的影响,发现土地利用转移变化越大则水文分区转移也越大,流域尺度越小时水文分区越离散;Di 等[47]在应用自组织映射神经网络(SOM)进行流域分类时发现,预先应用主成分分析和典型相关分析(Canonical Correlation Analysis,CCA)可以提高SOM 应用效果,降低不确定性。
分区方法总体上经历了人工判别、简单人工智能方法、组合人工智能方法等阶段,后两种方法目前应用较多。相比于人工判别方法,人工智能方法能够快速完成大批量分区任务,基于给定分区指标实现快速分类。然而其缺点也很明显,即无法考虑流域的实际水文规律。在分区方法的适用性方面,当前研究主要是案例分析,尚无一般性结论,方法的适用性主要通过检验分区结果的合理性体现。
从方法上看,与水文分区相似的为流域相似性评价,后者指基于水文相似理论评价流域之间的水文相似性,常用于相似流域选择,涉及相似性指标、相似度评价及流域分类方法等内容[48]。从理论上说,同一个水文分区的流域应具有水资源时空分布、洪峰特征、产汇流过程等方面的相似性,而水文分区之间的流域则存在较大差异。因此,可以应用流域相似性评价方法对水文分区结果进行合理性检验,目前还少有研究。
2.3 水文分区的应用进展
水文分区主要有两大方面的应用:水文水资源规划、无资料流域水文研究。前者属早期的主要应用领域,后者是当前的主要应用领域,因此本文重点关注后者。在PUB 计划(Prediction in Ungauged Basins,2003—2012)收官之际,Hydrology and Earth System Sciences期刊出版了“Catchment classification and PUB”专辑(https://hess.copernicus.org/articles/special_issue136.html),其中绝大多数研究与水文分区(即流域分类)直接相关,可见水文分区对PUB 研究的重要意义。
水文分区在无资料流域水文研究中的应用可分为理论和实践两个层面,理论层面是指发展和认识新的水文理论与规律,如水文响应特征的区域性分布、水文参数的气候-地貌规律等。Laaha 等[35]开展了Q95低流量的流域分类,发现基于季节性区域的分组效果最好,表明低流量过程特征存在季节性差异和区域性差异;陈旭等[26-27]将海河流域划为5 个水文分区,分析了各分区的洪水演变特性与暴雨洪水响应关系,发现人类活动对海河流域洪水具有显著削减作用。实践层面是指获取无资料流域的水文参数方面,如水文站网优化、分区内水文参数区域化等应用。姬海娟等[19]通过水文分区研究发现雅江流域的水文站网布设尚无法满足资料与参数移植;李萍等[49]将山东省划分为2 个1 级水文区、6 个2 级水文亚区,为公路建设安全性提供了参考信息;李捷等[50]通过改进青海省水文分区,提高了缺资料地区公路桥涵工程设计流量合理性;武光明[51]将山西省分为5 个产沙分区,建立了可用于无资料地区的基于水文下垫面的地区雨沙模型;伍远康等[52]将浙江省分为2 个洪水计算区,成果可用于省内缺资料地区洪水计算;余江游等[53]将淮河流域分为4 个水文相似区,从而实现了区域内洪水频率曲线拟合和设计洪水计算,克服了资料缺乏问题。以上研究解决了水文观测稀缺地区水文模拟、工程水文设计等方面的资料难题。
分区原则是水文分区应用中的指导思想体系,陈异植等[54]按照比较水文学的观点提出综合性、相对一致性、区域共轭性、发生学、主导因素、河系完整性等分区原则;张静怡等[6]则将其归纳为综合分析、相似性与差异性、成因分析、主导因素和河系完整性等原则;赵宝君[14]提出以“集体效应”为标准、遵循“从主、从众、从源”的原则,即根据实际与需求进行差异化分区。综上,没有一概而论的普适性分区原则,只有符合研究需求的适用性分区原则。
由于水文指标与生态指标存在大量交叉,且要素之间具有较强的相关性,因此生态水文分区常作为生态分区的替代方案,用于生态环境评价与保护。由于生态水文分区在本质上与水文分区同属于模式识别问题,因此其研究与应用方法与水文分区也类似,主要差异在于应用领域。如尹民等[55]构建了全国河流生态水文分区,综合揭示了不同区域水资源条件、生态系统功能及不同河流生态系统的需水特征;杨爱民等[56]首先提出了生态水文分区及其构建方法,将全国划分为3 个生态水文大区、36 个生态水文小区,揭示了生态水文复合系统的空间分异规律。还有其他诸多区域性生态水文分区研究,可认为是水文分区与生态分区的拓展应用,为中国不同区域定向生态环境保护与水资源管理提供了科学依据,此处不再赘述。
3 结 语
本文首先回顾了中国水文区划的发展历程,并按照水文分区研究水平及涉及的区域划分为全国性水文区划初探时期、区域性水文区划细化时期、区域性水文区划改进时期3 个阶段。已有的水文分区研究呈现多样性趋势。从水文分区指标、分区方法及应用3 个方面综述了水文分区的研究现状,当前研究中采用的水文分区指标主要包括气候、水文、下垫面三类,也有研究者将基于水文变量的统计参数等衍生指标作为水文分区指标,具体采用的指标一般取决于研究目的。水文分区的分类方法以聚类相关方法为主,也可见其他智能算法用于水文分区研究并获得了较好的结果,但针对分区方法的适用性研究较少。在应用上,当前主要用于缺资料流域水文研究,如水文水资源空间规律、水文参数区域化、工程设计参数移用等;同时还结合生态分区发展了生态水文分区。基于水文分区的研究与应用现状,其未来发展趋势和研究重点如下:
(1)加强水文分区的相似理论基础。水文分区与相似流域是一脉相承的,水文分区的研究中通常较少考虑水文相似性,有时过于注重将区域划分为“几类”。研究表明基于水文相似性的流域分类通常能获得较好的结果,因此在水文分区中应充分纳入水文相似理论。
(2)加强水文区域性规律研究与应用。大量研究表明,基于水文分区的水文水资源、生态与环境科学、自然地理学、工程规划与设计等研究及应用均取得了更好的结果。根据研究需求挖掘出相应的水文区域性规律,既能提高对目标问题的水文背景认识,也能促进上述研究及应用中涉及到的环境要素对水文水资源时空演变的影响研究。
(3)加强水文分区结果的检验与科学解释。水文分区研究得到的是多个流域群组,每个群组的成员流域在理论上应具有某些方面的相似特征,应根据最初的研究目标对各个群组进行检验。除了群组内应具有一致性以外,还要检验群组间的差异是否能够反映分区目标。对比群组间差异的成因也是检验分区结果可靠性的途径。