粤港澳大湾区金山湖流域水质变化规律

张万顺，张紫倩，彭 虹，李 琳，张 潇，夏 函，章 玲

(1.武汉大学资源与环境科学学院，湖北 武汉 430072；2.武汉大学水资源与水电工程科学国家重点实验室，湖北 武汉 430072；3.武汉大学中国发展战略与规划研究院，湖北 武汉 430072； 4.武汉大学水利水电学院，湖北 武汉 430072)

粤港澳大湾区是我国城镇化水平最高的城市群区域[1]，人口集聚度高、河湖水系发达、自然地貌复杂，随着人水关系逐渐增强，水资源短缺、水体黑臭、生态服务功能下降等城市河湖水环境问题日益严峻[2-5]。受社会经济过程与自然过程交叉影响，城市河湖水质变化规律愈发复杂，存在水质变化机制不确定的突出问题，成为制约粤港澳大湾区高质量发展的瓶颈[6-8]。精准揭示城市河湖独特的水文情势和强人类活动相互影响下的水质变化规律，涉及河湖水系交汇、陆面降雨径流、土地利用格局和水库调度等多过程[9-11]，亟须对气象、水文、水质等多要素全过程进行全面精细的量化研究。

近年来，国内外学者针对精细量化水环境变化规律开展了大量研究，主要集中在自然条件下大尺度流域降水、径流、水量、水质等过程的定量模拟[12-13]。城市河湖存在水文、生态、社会异质性特征复杂多变等特点，为精细刻画城市河湖水质变化过程，需要统筹陆面单元大尺度、河网中尺度、重要水域小尺度水环境过程机制差异的模型体系[14-15]。基于耦合气象数值模型、陆地模型和水动力-水质-水生态模型的“空-地-水”一体化模型逐渐成为全面揭示水质变化规律的重要工具[16-18]，可突破流域陆面空间上不同尺度单元耦合的瓶颈，实现流域自然过程与社会经济过程的耦合与响应关系的定量模拟[19-20]，能够弥补传统流域水环境模型尺度单一、唯自然过程的不足。

金山湖流域是位于粤港澳大湾区的典型城市流域，河湖水系交错，社会经济发展空间分异明显，水质不达标、水体黑臭等水环境问题较为严重。本文以金山湖流域为例，建立了耦合陆面单元、河网、湖泊等不同尺度的“空-地-水”一体化模型，模拟强人类活动、水文、气象过程等多重因素影响下流域内污染负荷以及水质变化过程，以揭示水质变化机制以及土地利用格局对水质的影响。

1 研究区概况

金山湖流域位于广东省惠州市城区南部，东江流域中下游，西枝江左岸。金山湖流域水系由金山湖与4条主要入湖河涌构成，入湖河涌为金山河、莲塘布、河桥水和冷水坑，流域集雨面积约 63.7 km2，水系如图1所示。流域上游为农业区域，中游为城乡接合部，下游为城市地区，点源排放密集，农业面源、城镇面源污染严重，人类活动影响显著。2019年，入湖4条河涌整体为劣Ⅴ类水体，水质状况均属于“轻度黑臭”；湖区水质整体为Ⅴ类，金山河入湖口至河桥水入湖口湖段水质属于“轻度黑臭”，全湖段均呈现出不同程度的富营养化现象。

图1 金山湖流域水系

2 “空-地-水”一体化模型构建

基于课题组前期研究成果[17,20-21]，综合考虑陆面降雨径流、河网及湖泊水动力演进与污染物输移，构建“空-地-水”一体化模型，包含一维河网与二维湖泊水动力水质模型以及面源模型。以流域面源模型的模拟结果作为水动力水质模型的旁侧入流条件，实现降雨径流面源模型与水动力水质模型的耦合计算；通过在一维河网和二维湖泊的连接断面处设置过渡单元，根据过渡单元上两种模型模拟的水位、流量、浓度相等的条件，实现一、二维模型的耦合。

2.1 面源模型

本文选用的流域面源模型为SWMM模型，包含水文与水质模块。选择地表汇流方式、Horton下渗模型进行地表径流模拟，选择指数函数方程进行污染物累积冲刷模拟[20-21]。采用30 m精度的DEM栅格，将金山湖流域陆面共划分为153个子汇水区，如图1所示。通过DEM数据、土地利用数据确定各子汇水区的面积、不透水面积占比、流域宽度、平均坡度等确定性参数；采用惠阳国家气象站逐日降雨、蒸发、气温、相对湿度等数据建立气象数据库。

2.2 一维河网与二维湖泊水动力水质模型

采用圣维南方程组描述水动力学过程，并充分考虑水体污染物的对流、扩散、降解作用，建立一维河网与二维湖泊水动力水质模型[17]。收集整理流域内水文地形资料，构建金山河、莲塘布、河桥水、冷水坑4条河涌的一维模型地形断面，共计75个；构建金山湖湖区的二维模型，划分正交结构网格共计46 624个，网格精度为20 m×20 m，地形断面与网格设置如图2所示。

2.3 模型参数率定与验证

根据金山湖流域水质现状，选取主要污染因子COD、TP作为模拟指标。利用金山湖2018年12月(枯水期)和2017年8月(丰水期)的GF-2号遥感影像水质反演结果(图3)做模型率定。流域内糙率取值范围为0.025～0.055，COD降解系数为0.01～0.4 d-1，TP降解系数为0.02～0.5 d-1。以2019年1月水量验证断面(图1)实测汇水水量进行流域面源模型验证，以4条河涌10个断面、湖区5个点位(图1)的水质监测结果进行一维河网与二维湖泊水动力水质模型验证。模型验证结果见图4，面源模型水量的相对误差均在5%以内，水动力水质模型COD和TP的相对误差均在15%以内，表明“空-地-水”一体化模型能够满足金山湖流域面源及水质模拟的要求。

(a) COD(2017年8月)

(a) 水量

3 水质变化规律分析

3.1 水质时空变化规律

由于降水是造成非点源污染的主要驱动因素，丰水期非点源污染尤为显著[22]，因此选择典型丰水年2017年作为研究时段。基于“空-地-水”一体化模型，模拟了2017年水文、气象条件下COD、TP质量浓度的变化过程。选取4条河涌上游、中游和下游典型断面开展水质变化分析，结果见图5和图6。金山湖2017年典型丰水月6月与典型枯水月12月的水质指标不同浓度范围面积占比如图7所示。由图5～7可见，金山湖流域水质整体较差，枯水期水质比丰水期差，金山湖在典型丰水月劣Ⅴ类面积占比均在86.14%以上，在典型枯水月劣Ⅴ类面积占比在96.66%以上。在春夏季降雨丰沛时期，河道与湖泊的流量、流速整体较大，对污染物的稀释作用较强。以地表水Ⅳ为标准，各河涌断面COD质量浓度年达标率均在16.71%以下，TP质量浓度年达标率均在2.47%以下；湖泊TP质量浓度超标严重，COD质量浓度整体较好。

(a) 金山河

(a) 金山河

(a) COD

水质波动与降水量密切相关，枯水期水质无明显波动，丰水期水质随降雨波动明显。2017年，日降水量为15～60 mm的中小雨共出现31次，在中小雨后的1～3 d，金山湖全流域的COD质量浓度易升高，TP质量浓度有所降低；降水量超过60 mm的大雨共出现8次，在大雨过后，流域的COD、TP质量浓度均有所降低。水质波动特征表明日降水量小于15 mm的降雨对陆面的冲刷作用不明显，带来的降雨径流污染影响较小，降雨期水质较为稳定；中小雨时期产生的陆面降雨径流污染程度最高，对水体的COD质量浓度影响显著；大雨冲刷作用产生更高的面源污染负荷，但由于此时降雨径流量较大，面源负荷对水质的影响程度有所降低。各河涌水质具有明显的空间分异特征，金山河、莲塘布的中上游水质较下游好，河桥水、冷水坑的中下游水质较上游好。金山湖丰水期与枯水期COD、TP质量浓度分布如图8所示，湖泊的南部出现超标污染带的污染物浓度较高、范围较广，西北部次之，东部水质相对较好。

(a) COD(2017-06-15)

金山湖流域污染物浓度较高的区域与点源排放密集区域一致，各区域点源负荷如表1所示，金山河、莲塘布的上游区域，河桥水、冷水坑的中下游区域，以及湖泊的南部与西北部区域点源负荷较高，表明点源是影响水质空间分布特征的主要原因。此外，湖区南部是河桥水、冷水坑、莲塘布等河涌集中入汇区域，河涌入汇携带大量污染负荷至湖泊，对湖泊水质造成不利影响。

表1 金山湖流域各区域点源负荷

3.2 土地利用格局对水质的影响

金山湖流域陆面以城镇用地为主，占全流域陆面面积的37.38%，其次为林地、耕地和农村居民点，占比分别为25.29%、24.27%和9.12%，荒地与草地占比仅分别为3.21%和0.75%。城镇用地主要分布在各支流的中下游，林地主要分布在各支流的中上游，耕地主要分布在河桥水中游、冷水坑上游以及莲塘布中上游，农村居民点主要分布在河桥水中游、冷水坑中下游以及莲塘布中上游。

从土地利用格局尺度分析，不同水质断面受区间入汇直接影响[23]，因此本文将每个断面的子流域范围定义为断面区间范围内符合产汇流机制的所有子流域。分别对4条河涌选取典型干流断面开展子流域土地利用格局分析，4条河涌所选断面以及不同断面子流域土地利用类型占比如图9所示。各支流从上游至下游，金山河5个断面子流域城镇用地、农村居民点面积占比逐渐增大，林地、耕地面积占比逐渐减小，草地、荒地面积占比低于2%；河桥水4个断面子流域城镇用地面积占比逐渐增大，农村居民点、耕地面积占比先增大后减小，林地占比逐渐减小，草地、荒地占比低于1%；冷水坑4个断面子流域城镇用地面积占比先增大后减小，林地、耕地面积占比逐渐减小，农村居民点面积占比逐渐增大，草地、荒地面积占比低于4%；莲塘布7个断面子流域城镇用地面积占比逐渐增大，耕地面积占比先增大后减小，林地面积占比逐渐减小，农村居民点面积先增大后减小，草地、荒地面积占比低于5%。

图9 不同断面子流域土地利用类型

将不同水质断面COD、TP质量浓度的年均值、典型丰水月6月均值和典型枯水月12月均值与对应子流域土地利用类型结构进行Pearson相关性分析，结果如表2所示。COD质量浓度与耕地、草地面积占比呈显著正相关关系，与林地面积占比呈显著负相关关系；TP质量浓度与耕地、农村居民点面积占比均呈显著正相关关系；其余土地利用类型面积占比与水质指标质量浓度的相关性不显著。耕地面积占比与水质指标质量浓度均呈正相关关系，耕地区域农业耕作活动频繁，化肥流失、农业面源污染严重，导致水体中的污染加重；丰水月COD质量浓度和TP质量浓度与耕地面积占比相关性高于枯水月和年均值，表明降水会加重农业面源污染；COD质量浓度的相关性较TP质量浓度更高，表明农业面源污染对COD质量浓度的影响程度更大。林地面积占比与水质指标质量浓度均呈负相关关系，林地区域人类活动较弱，植被覆盖度高，水土保持能力强，植被对污染物的截流以及净化作用可以有效减少污染物。农村居民点面积占比与水质指标质量浓度均呈正相关关系，农村居民点承载了较多农村居民生活、畜禽养殖活动，给水环境带来大量农村生活污染。城镇用地面积占比与水质指标质量浓度相关性不显著，城镇带来的污染以工业废水和市政污水点源为主，而流域内点源在支流各区域分布不均。流域内草地与荒地由于面积占比极小，与水质指标质量浓度相关性的不确定性较大。

表2 水质指标与土地利用类型面积占比的相关系数

4 结 语

本研究建立的耦合陆面单元、河网、湖泊等不同尺度的“空-地-水”一体化模型能够实现对粤港澳大湾区强人类活动影响下的金山湖流域河湖水质变化过程的精细化模拟。模型综合考虑了陆面降雨径流、河网及湖泊水动力演进与污染物输移的影响，能有效反映流域水质变化机制。模拟结果表明，金山湖流域水质与流域降雨和土地利用类型密切相关，为推进流域水环境精细化管理提供有力的技术支撑，助力形成“人水和谐”的新时代生态文明格局。