北江干流及主要支流水质及其时空变化特征.pdf

第4 3卷第5期2 0 2 3年1 0月水土保持通报B u l l e t i no fS o i l a n dW a t e rC o n s e r v a t i o nV o l.4 3,N o.5O c t.,2 0 2 3 收稿日期:2 0 2 3-0 1-1 3 修回日期:2 0 2 3-0 3-1 4 资助项目:广东省省级科技计划项目“农业面源污染绿色防控生态农场应用与示范”(2 0 2 2 B 0 2 0 2 1 6 0 0 0 7);广东省自然科学基金项目“沼渣生物炭限 域 单 原 子 铁 活 化 过 硫 酸 盐 降 解 抗 生 素 的 效 能 及 机 制”(2 0 2 2 A 1 5 1 5 0 1 0 8 0 7);广 东 省 重 点 领 域 研 发 计 划 项 目(2 0 2 0 B 1 1 1 1 3 5 0 0 0 1)第一作者:周铭浩(1 9 9 0),男(汉族),湖南省常德市人,博士,主要从事水生态环境质量监测与预警预报等研究。E m a i l:9 4 2 4 1 4 8 8 2 q q.c o m。通信作者:林学明(1 9 8 9),男(汉族),广东省汕尾市人,博士,副教授,主要从事水生态环境修复等研究。E m a i l:x m.l i n s c a u.e d u.c n。北江干流及主要支流水质及其时空变化特征周铭浩1,李 彤1,林学明2,叶四化1,张雪容1,邓 滢1(1.广东省生态环境监测中心,广东 广州5 1 0 3 0 8;2.华南农业大学 资源环境学院,广东 广州5 1 0 6 4 2)摘 要:目的分析2 0 2 0年北江干流及主要支流监测断面的水质情况,为北江水环境监控及保护提供参考依据。方法采用单因子评价法和主成分分析法,研究北江水质时空变化趋势及驱动因素。结果2 0 2 0年北江监测断面整体水质良好,仅水车头、漫水河河口、黄浪水桥、龙归断面水质为类,其他断面均为类。北江水温和溶解氧受季节影响较大,水体全年基本保持弱碱性,下游营养盐和有机物污染程度高于上、中游。结论由空间分析结果可知,下游水质劣于上、中游,其主要的影响因素是营养盐和溶解氧,p H值单独作用不会对水质产生明显影响。季节污染程度排序为:夏季秋季春季冬季,其中全河段的水质主导因素为溶解氧,而下游河段的主导因素为营养盐和有机物。关键词:北江;水质评价;主成分分析;时空变化文献标识码:A 文章编号:1 0 0 0-2 8 8 X(2 0 2 3)0 5-0 1 1 1-1 0 中图分类号:X 5 2 2文献参数:周铭浩,李彤,林学明,等.北江干流及主要支流水质及其时空变化特征J.水土保持通报,2 0 2 3,4 3(5):1 1 1-1 2 0.D O I:1 0.1 3 9 6 1/j.c n k i.s t b c t b.2 0 2 3.0 5.0 1 4;Z h o uM i n g h a o,L iT o n g,L i nX u e m i n g,e ta l.W a t e rq u a l i t ya n ds p a t i a l-t e m p o r a lv a r i a t i o nc h a r a c t e r i s t i c sf r o m m a i ns t r e a ma n d m a i nt r i b u t a r i e so fN o r t hR i v e rJ.B u l l e t i no fS o i l a n dW a t e rC o n s e r v a t i o n,2 0 2 3,4 3(5):1 1 1-1 2 0.W a t e rQ u a l i t ya n dS p a t i a l-t e m p o r a lV a r i a t i o nC h a r a c t e r i s t i c sf r o m M a i nS t r e a ma n dM a i nT r i b u t a r i e so fN o r t hR i v e rZ h o uM i n g h a o1,L iT o n g1,L i nX u e m i n g2,Y eS i h u a1,Z h a n gX u e r o n g1,D e n gY i n g1(1.G u a n g d o n gE c o l o g i c a la n dE n v i r o n m e n t a lM o n i t o r i n gC e n t e r,G u a n g z h o u,G u a n g d o n g5 1 0 3 0 8,C h i n a;2.C o l l e g eo fN a t u r a lR e s o u r c e sa n dE n v i r o n m e n t,S o u t hC h i n aA g r i c u l t u r a lU n i v e r s i t y,G u a n g z h o u,G u a n g d o n g5 1 0 6 4 2,C h i n a)A b s t r a c t:O b j e c t i v eT h ew a t e rq u a l i t yo ft h em a i ns t r e a ma n dm a i nt r i b u t a r ym o n i t o r i n gs e c t i o n so ft h eN o r t hR i v e r i n2 0 2 0w a sa n a l y z e d i no r d e r t op r o v i d ear e f e r e n c ea n db a s i s f o rm o n i t o r i n ga n dp r o t e c t i n gt h ew a t e re n v i r o n m e n t o f t h eN o r t hR i v e r.M e t h o d sT h es i n g l e f a c t o re v a l u a t i o nm e t h o da n dp r i n c i p a lc o m p o n e n ta n a l y s i sw e r eu s e dt od e t e r m i n et h es p a t i a l a n dt e m p o r a l t r e n d sa n dd r i v i n gf a c t o r so fw a t e rq u a l i t yf o rt h eN o r t hR i v e r.R e s u l t sI n2 0 2 0,t h eo v e r a l lw a t e rq u a l i t yo f t h eN o r t hR i v e rm o n i t o r i n gs e c t i o n sw a sg o o d.O n l yt h es e c t i o n so fS h u i c h e t o u,M a n s h u iR i v e re s t u a r y,H u a n g l o n g s h u iB r i d g e,a n dL o n g g u im e tt h es t a n d a r d sf o rw a t e rq u a l i t yg r a d e s,a n dt h eo t h e rs e c t i o n sw e r ed e t e r m i n e dt ob eg r a d e.W a t e rt e m p e r a t u r ea n dd i s s o l v e do x y g e n i nt h eN o r t hR i v e rw e r e i n f l u e n c e db yt h es e a s o n s.T h ew a t e rb o d yb a s i c a l l yr e m a i n e dw e a k l ya l k a l i n e t h r o u g h o u t t h ey e a r,a n dt h ep o l l u t i o nl e v e l so fn u t r i e n t sa n do r g a n i cm a t t e rw e r eh i g h e rd o w n s t r e a mt h a na tu p s t r e a ma n d m i d s t r e a ml o c a t i o n s.C o n c l u s i o nS p a t i a la n a l y s i ss h o w e dt h a td o w n s t r e a m w a t e rq u a l i t yw a sw o r s et h a nu p s t r e a m a n d m i d s t r e a m w a t e rq u a l i t y.T h e m a i ni n f l u e n c i n gf a c t o r sw e r en u t r i e n t sa n dd i s s o l v e do x y g e n.p Hv a l u e a l o n ed i dn o th a v eas i g n i f i c a n t i m p a c to nw a t e rq u a l i t y.T h es e a s o n a l p o l l u t i o nd e g r e e f o l l o w e dt h eo r d e ro fs u mm e r a u t u m n s p r i n g w i n t e r.T h ed o m i n a n tf a c t o ro n w a t e rq u a l i t yf o rt h ee n t i r er i v e r w a sd i s s o l v e do x y g e n,w h i l et h ed o m i n a n tf a c t o r sf o rt h ed o w n s t r e a mr i v e rs e c t i o nw e r en u t r i e n t sa n do r g a n i cm a t t e r.K e y w o r d s:N o r t hR i v e r;w a t e rq u a l i t ye v a l u a t i o n;p r i n c i p a l c o m p o n e n t a n a l y s i s;s p a t i a l-t e m p o r a l v a r i a t i o n 作为广东省5大河流之一,北江在该省的社会经济发展和生态安全方面具有重要地位。北江水资源量丰富,年平均径流量达到4.8 21 01 0m3,是省内农业、工业、市政用水的重要水源1。北江流域是广东省“一核一带一区”区域发展格局中北部生态发展区的组成部分,其功能定位是省内的重要生态屏障2。近年来,北江流域(尤其下游接近珠江三角洲地区)的水域、林地、草地等面积逐渐下降,耕地面积逐渐上升,带来水环境污染风险3。北江干流繁多的水利工程也对河流的水生态环境造成不利影响4。了解北江水质的时空变化特征、分析水体污染源有助于提升水资源保护和管理的水平。目前,水体水质分析评价方法主要包括单因子评价法、综合污染指数法、主成分分析法、灰色关联分析法、模糊数学评价法、人工神经网络等。国内外学者采用多种方法分析水体水质变化特征。李晓玉等5综合运用单因子评价法、综合污染指数法和主成分分析法对温榆河北运河段进行水质分析,提取主导水质变化的影响因子;刘宇等6运用灰色关联分析法和主成分分析法分析了云南省抚仙湖、星云湖和杞麓湖的水质时空分布特征;武俐等7采用主成分分析法和模糊数学评价法综合分析了小浪底水库水质状况和影响因素;G a r c i a-A l b a等8运用人工神经网络分析了西班牙西北海岸河口水质时空分布特征。结合多种评价方法能够更全面地评价北江水质状况。本研究采用单因子评价法和灰色关联分析法对北江干流及主要支流水质进行分析评价。基于北江沿程3 2个监测断面2 0 2 0年11 2月的逐月水质数据,使用主成分分析法探究了北江水质的空间和时间变化趋势,并结合北江沿岸环境特点讨论这些变化的驱动因素,以期为加强北江水环境监控、治理与保护提供科学依据。1 材料与方法1.1 研究区概况北江是珠江流域的第二大水系,发源于江西省信丰县石碣大茅坑,经大余县进入广东,在广东境内流经南雄市、始兴县、仁化县、浈江区、武江区、曲江区、英德市、清新区、清城区、三水区、四会市,在思贤滘注入珠江三角洲网河区。本次北江研究区域位于广东省境内(1 1 2 4 8 5 3 1 1 4 4 1 0 9 N,2 3 0 9 1 3 2 5 1 9 1 5 E),主要支流包括墨江、锦江、武江、南水、滃水、连江、潖江、滨河、大燕河、漫水河、绥江等。其中,韶关市沙洲尾以上河段为北江上游、沙洲尾至清远市飞来峡河段为北江中游、飞来峡以下河段为北江下游。北江流域多年平均降水量约为18 0 0mm,强降水大部分发生在汛期(58月),占年降水量6 0%7 0%9。目前,北江流域已建成孟洲坝水利枢纽、濛里水利枢纽、白石窑水利枢纽、飞来峡水利枢纽、清远水利枢纽、锦江水库、小坑水库、乐昌峡水利枢纽、南水水库、长湖水库、潭岭水库、锦潭水库共1 2宗大型水库,总库容为5.5 61 09m3。其中,孟州坝、濛里、白石窑、飞来峡、清远水利枢纽位于北江干流,锦江水库、小坑水库、乐昌峡水利枢纽、南水水库、长湖水库、潭岭水库、锦潭水库分别位于北江支流锦江、枫湾河、武江、南水河、滃江、连江、黄洞河。此外,北江流域已建成中型水库5 0宗,总库容为1.6 11 09m3;已建成小(一)型水库1 8 4宗,小(二)型 水 库7 2 1宗,总 库 容 为2.6 11 09m3。北江流域水库群优化调度可带来防洪、供水、航运、旅游、环保等综合效益。1.2 数据采集本研究选取了北江干流及主要支流沿程3 2个断面作为研究对象,各监测断面位置和基本信息见图1和表1。图1 北江监测断面位置示意图F i g.1 S c h e m a t i cd i a g r a mo fm o n i t o r i n gs e c t i o nl o c a t i o ni nN o r t hR i v e r2 0 1 22 0 2 1年北江流域降雨量范围在15 3 923 8 2mm,平均降雨量为19 5 9.7mm;地表水资源量范围在3.6 51 01 06.4 51 01 0m3,平均地表水资源量为5.2 11 01 0m31 0。2 0 2 0年北江流域降雨量为18 6 0.4 mm,地表水资源量为5.0 41 01 0m3,接近211 水土保持通报 第4 3卷2 0 1 22 0 2 1年降雨量和地表水资源量的平均值,由此选择2 0 2 0年水质数据研究北江水质状况的变化特征。监测断面2 0 2 0年水质数据来源于广东省生态环境监测中心,水样采集方法依据 地表水环境质量监测技术规范(H J9 1.2-2 0 2 2)1 1,每月采样1次,水质监测指 标 包 括 地 表 水 环 境 质 量 标 准(G B 3 8 3 8-2 0 0 2)中地表水环境质量标准基本项目2 4项。其中水温(WT)、溶 解 氧(D O)、p H值、高 锰 酸 盐 指 数(C O DM n)、生化需氧量(B O D5)、氨氮(NH+4-N)、总氮(T N)、总磷(T P)共8项指标均有完整检出值,作为主要分析指标。根据 水和废水监测分析方法 要求,WT测定采用温度计法,D O测定采用电化学探头法,p H值测定采用电极法,C O DM n测定采用酸性高锰酸盐法,B O D5测定采用稀释与接种法,NH+4-N测定采用纳氏试剂光度法,T N测定采用碱性过硫酸钾消解紫外分光光度法,T P测定采用钼酸铵分光光度法1 2。由于北江流域所在的韶关市和清远市属于广东省“北部生态发展区”,控制建设重金属污染相关企业,且2 0 2 0年研究断面的重金属铜、锌、汞、镉、六价铬、铅检出率仅3 8.8%,所以本文主要研究居民生活及农业生产排放的相关水体污染物指标。上述8项水质监测指标中,WT,D O,p H值属于基础理化指标,NH+4-N,T P,T N属于营养盐指标,C O DM n,B O D5属于有机物指标。表1 北江监测断面基本信息T a b l e1 B a s i c i n f o r m a t i o no f i n v e s t i g a t e ds e c t i o n s i nN o r t hR i v e r编号河流断面名称 编号河流断面名称 1北江干流河坪(上游)1 7北江干流沙口(中游)2北江干流南雄(上游)1 8北江干流英德(中游)3北江干流古市(上游)1 9滃 江滃江大坝(中游)4北江干流小古菉(上游)2 0连 江高道(中游)5墨 江墨江出口(上游)2 1北江干流黎溪(中游)6锦 江瑶山电站(上游)2 2北江干流飞来峡水库(中游)7北江干流长坝(上游)2 3潖 江潖江口(下游)8北江干流黄浪水桥(上游)2 4北江干流七星岗(下游)9北江干流新韶水文站(上游)2 5滨 河飞水桥(下游)1 0武 江武江桥(上游)2 6大燕河水车头(下游)1 1北江干流九公里(中游)2 7北江干流石角(下游)1 2北江干流白土(中游)2 8漫水河漫水河河口(下游)1 3南水龙归(中游)2 9北江干流北江水厂(下游)1 4北江干流白沙(中游)3 0北江干流芦苞(下游)1 5北江干流大坑口(中游)3 1绥 江马房(下游)1 6北江干流高桥(中游)3 2北江干流梁村(下游)1.3 数据分析及评价采用单因子评价法对研究断面水质状况进行分析评价。单因子评价法根据 地表水环境质量标准(G B 3 8 3 8-2 0 0 2)指标限制对水质进行分类,选取指标中最高的类别作为水质评价结果1 3-1 4。采用灰色关联分析法评价研究断面的水质状况。步骤为:选择能体现水质变化情况的因子(D O,C O DM n,B O D5,NH+4-N,T P)作为参考序列,选择水质标准类别(类)作为比较序列;归一化及无量纲化处理所选水质数据;根据沈珍瑶等1 5提出的方法计算绝对差;使用平均值法求得关联度,并将关联度最大值所对应的标准作为研究断面的水质类别。采用主成分分析法进一步分析北江研 究区域水质的主导影响因素。主成分分析法是通过重组和综合变量的方式,将相关性高的变量分在同一类因子中,这些新的公共因子保留了原来变量中存在的大部分变化,达到全面反映水质状况和污染因素的目的1 6。本文利用软件E x c e l 2 0 2 1处理所得试验数据,采用软件A r c G I SP r oV 2.5和O r i g i n2 0 2 1进行图形分析与绘制,利用软件I BM S P S SS t a t i s t i c s2 7对数据进行主成分分析。2 结果与分析2.1 北江干流及主要支流水质状况分析由2 0 2 0年北江干流及主要支流逐月监测数据可知水质指标波动较大,其中WT范围为1 03 4,311第5期 周铭浩等:北江干流及主要支流水质及其时空变化特征D O范围为4.6 51 1.6 3m g/L,p H值范围为6.0 08.7 6,C O DM n范围为0.65.2 m g/L,B O D5范围为0.25.9m g/L,NH+4-N范围为0.0 21.5 4 m g/L,T P范围为0.0 0 80.3 9 5m g/L,T N范围为0.2 5 39.6 5 0m g/L。基于地表水质量标准(G B 3 8 3 8-2 0 0 2),采用单因子评价法分析北江监测断面的水质类别。根据目前的评价规范,总氮不作为日常水质类别评价指标,仅作为参考指标单独评价。从全年平均值来看,水车头和漫水河河口为类,关键影响指标为总磷;黄浪水桥、龙归和潖江口水质类别为类,其中黄浪水桥和潖江口的关键影响指标为总磷,龙归的关键影响指标为氨氮和总磷;其余监测断面均属于类。为了进一步分析2 0 2 0年研究断面的水质状况,采用灰色关联分析法对北江流域各断面水质等级进行评价(图2)。如图2所示,在全年中,类水断面比例达到8 4%,其他超类断面主要集中在北江支流,所得结果与单因子评价法一致。图2 北江流域全年各断面水质等级评价结果F i g.2 A n n u a lw a t e rq u a l i t yg r a d ee v a l u a t i o nr e s u l t so fe a c hs e c t i o n i nN o r t hR i v e rb a s i n根据2 0 2 0年北江监测断面不同季节、不同河段的水质数据制作箱线图(图3),以获取直观的水质变化特征。由图3可知,WT均值基本呈现:冬季春季秋季春季秋季夏季,这是因为水温升高不仅降低了水体饱和溶解氧浓度,而且促进了水体中微生物的耗氧过程1 9-2 0。北江监测断面的D O普遍处于类标准,说明研究区域水体含氧量比较充足。不同季节、河段监测断面的p H值集中于78,总体表现为弱碱性。上游p H均值小于中、下游,可能的原因是中、下游水体含有更丰富的营养盐和浮游藻类,光合 作用过程中 吸收C O2并导致p H值 升高2 1。各监测断面的C O DM n和B O D5大部分时间段都处于类标准限值以内,且这两个指标的中游均值小于上、下游,表明北江上、下游的有机物输入量多于中游。此外,下游监测断面B O D5数据波动大于上、中游,反映下游不同断面B O D5差异性较大。对于NH+4-N和T P而言,除了少量偏高的异常数据达到类、类标准,大部分数据处于类标准以内,甚至达到类 标 准。营 养 盐 指 标NH+4-N和T P质 量总体良好。监测断面T N均值呈现:上游中游秋季春季冬季。由季节变化考察水质污染主导成分,可知夏、秋两季的主导成分是D O,由夏转秋的过程中D O的影响力逐步下降,营养盐和有机物的作用逐渐上升。如前所述,夏季WT在四季中偏高(图3 a),导致水体中D O的含量偏低,进而影响水质2 9。WT与D O之间的负相关关系也在图6 a中得到证实。北江流域的降雨主要集中于夏季,不仅会形成地表径流,使陆地上的营养物质和有机物进入水体,也会增加水体水量,稀释其中的污染物浓度3 0。由此可知,北江全河段夏季降雨的稀释作用降低了陆源营养物和有机物的污染性,秋季降雨量减弱则提高了水体营养物和有机物的浓度。相比于夏、秋季,冬、春季的水质明显提高。如表5所示,冬季的主导污染成分是营养盐,而春季的主导成分则变为有机污染物,但是低温带来的高D O条件提升了水体的质量。冬季水温下降,水生植物的衰亡不仅降低营养盐的吸收能力,而且会释放氮、磷物质,导致水体中营养盐浓度上升3 1;春季水温有所回升,水生植物逐渐恢复,并通过光合作用把营养物质转化为有机物3 2。综上所述,对北江全河段而言,营养物和有机物对水体水质的影响较弱,季节、温度变化引起的D O浓度差异是水质变化的主导因素。图6 b6 e置信椭圆表示不同季节各河段水质样本的协方差结构。结合表5可知各季节水体污染排序大致为:下游中游上游,仅夏季中游水质稍好于上游。研究发现北江中游两岸植被林地和草地较多,耕地和居民生活区相对较少3。夏季较大的雨量所产生的稀释作用进一步优化了该河段的水体质量。各季节上、中游河段水质的主导成分是D O,与全河段基本一致;下游河段夏、秋季体现营养盐、有机物与D O共同污染作用,冬、春季则体现营养盐和有机物的主导污染作用。3 讨 论北江流域地形多为山地与丘陵,地势北高南低,经过飞来峡后河道逐渐平缓。2 0世纪9 0年代以来,北江干流由上至下开发建设孟州坝、濛里、白石窑、飞来峡、清远等梯级水坝,坝址多年平均径流量依次为1.3 01 01 0,1.6 11 01 0,1.6 81 01 0,3.4 71 01 0,4.0 41 01 0m3。梯级水库群建设降低了水流流速,有助于提高北江干流堤防稳定性,但也可能影响陆源污染物在水体中的稀释扩散作用,降低水体纳污降解能力。北江干流流经韶关、清远、肇庆、佛山中的1 1个区县,总面积2.0 81 04k m2,人口约6.1 41 06人,上游、中游和 下 游 沿 岸 人 口 密 度 分 别 为2 0 0,2 2 0,7 3 8人/k m2。北江流域生产活动基本以种植业、畜禽养殖业、水产养殖业等农业生产为主。流域降雨强度较大,一方面会导致种植业农药化肥、养殖业粪污以及811 水土保持通报 第4 3卷农村生活污染物等面源污染进入河流,降低水体质量,另一方面又会增加河流水量,降低污染水平。由于北江流域水量较大,且人口密度较低,生产生活带来的污染强度较低,总体来看水质良好。图6 北江不同季节水质分布特征F i g.6 W a t e rq u a l i t yd i s t r i b u t i o no fN o r t hR i v e r i nv a r y i n gs e a s o n s 本文分析发现,北江下游水质总体劣于上中游,主要的污染指标为营养盐和D O,表明下游较高的人口密度提升了生产生活污染强度,污染物带有明显的农业面源特征。此外,研究发现下游断面水车头及漫水河河口指示的大燕河及漫水河有机物污染较为显著,表明这两条北江支流沿岸污染物排放量较大、收集能力较弱,亟需建设污水管网等收集设施,降低污染物入河率。广东省的主要粮食作物是水稻。早稻种植时间段为47月,晚稻种植时间段为81 1月,种植过程中包含基肥、蘖肥、穗肥、粒肥4个施肥阶段。尿素和磷酸二氢钾是重要的氮磷肥料。由此可知,夏季、秋季及春季部分时间会产生较强的农业面源氮磷污染。北江流域主要位于韶关和清远,可根据两市的月降雨量计算北江流域不同季节降雨情况(图7)。北江流域春季、夏 季、秋 季、冬 季 降 雨 量 分 别 为13 6 4.9,13 1 5.3,3 9 9.2,4 2 1mm。由表5可知,北江全河段污染程度季节排序为:夏季秋季春季冬季。夏季降雨量高于秋季,叠加种植业施肥过程影响,导致北江夏季污染程度高于秋季;尽管春季降雨量居于四季首位,但只有部分时段与施肥时间段重合,所以春季污染程度低于夏秋两季;冬季降雨量较少,且处于农闲期,污染程度最低。图7 2 0 2 0年韶关和清远月降雨量F i g.7 M o n t h l yr a i n f a l l o fS h a o g u a na n dQ i n g y u a ni n2 0 2 04 结 论(1)2 0 2 0年北江干流及主要支流整体水质良好,3 2个监测断面中,按单因子评价法分析,水车头及漫水河河口断面达到类水质标准,黄浪水桥及龙归断面达到类水质标准,其他干支流断面水质均为类。根据箱线图可知,季节变化对WT和D O的影响911第5期 周铭浩等:北江干流及主要支流水质及其时空变化特征较为显著,水体总体保持弱碱性,上、中游营养盐和有机物指标普遍优于下游。(2)使用主成分分析法研究8项主要监测指标,北江研究区域主要受到4个主成分的影响,其方差贡献率分别为2 8.9 9 9%,2 2.2 8 2%,1 9.0 7 1%,1 3.3 5 1%,其中第一主成分反映水体营养盐信息,第二主成分反映水体有机物信息,第三、四主成分反映水体理化性质。(3)通过北江监测断面的主成分/综合得分,得出水质空间变化状况,总体呈现下游比上、中游污染大的趋势,营养盐和D O是下游水质的重要影响因素。由4个主成分得分的系统聚类分析可得出不同监测断面之间污染因素的相似性,p H值单独作用不会对水质产生明显影响。(4)根据综合得分模型分析季节变化对北江河段的影响,可知全河段夏季和冬季的水质差异最大,污染程度季节排序为:夏季秋季春季冬季,水质主导因素是D O;下游河段夏、秋季水质由营养盐、有机物和D O共同作用,冬、春季水质则主要受到营养盐和有机物的影响。参考文献1 赵清贺,冀晓玉,丁圣彦,等.北江干流河岸带植物物种多样性的纵向梯度效应J.生态 学杂 志,2 0 1 8,3 7(1 2):3 6 5 4-3 6 6 0.2 吴锦泽,王明旭,叶晓颖.广东省北部生态发展区保护与发展协同共进研究J.环境与可持续发展,2 0 2 1,4 6(6):6 9-7 4.3 黎明.基于SWA T的北江流域土地覆盖及气象条件变化的水文响应模拟研究D.广东 广州:中国科学院研究生院(广州地球化学研究所),2 0 1 5.4 高雪山,刘树锋,贾建辉,等.基于生态水力半径法的武江流域生态流量研究J.水资源与水工程学报,2 0 2 0,3 1(2):8-1 2.5 李晓玉,王鹏,邵光艺,等.基于水面无人船快速监测的温榆河北运河水质分析与评价J.环境化学,2 0 2 2,4 1(5):1 5 6 8-1 5 7 8.6 刘宇,吉正元,刘淑娟,等.三种方法在高原湖泊水质评价中的应用与比较J.海洋湖沼通报,2 0 2 0(2):1 6 6-1 7 4.7 武俐,王祖恒,王亮,等.基于主成分分析和模糊数学的黄河小浪底水质监测与评价J.水土保持通报,2 0 2 0,4 0(5):1 1 8-1 2 4.8 G a r c a-A l b aJ,B r c e n aJF,U g a r t e b u r uC,e t a l.A r t i f i c i a ln e u r a l n e t w o r k sa se m u l a t o r so fp r o c e s s-b a s e d m o d e l st oa n a l y s eb a t h i n g w a t e rq u a l i t yi ne s t u a r i e sJ.W a t e rR e s e a r c h,2 0 1 9,1 5 0:2 8 3-2 9 5.9 王亚雄,黄淑娴,刘祖发,等.变化环境下北江下游年径流量的加权马尔可夫链预测J.生态环境学报,2 0 1 1,2 0(4):7 5 4-7 6 0.1 0 广东省 水利 厅.水资 源公 报 E B/O L2 0 2 3-7-3 0.h t t p:s l t.g d.g o v.c n/s z y g b/i n d e x.h t m l.1 1 中华人民共和国生态环境部.地表水环境质量监测技术规范:H J 9 1.22 0 2 2S.北京:中国环境科学出版社,2 0 2 2.1 2 国家环境保护总局水和废水监测分析方法编委会.水和废水监测分析方法M.北京:中国环境科学出版社,2 0 0 2.1 3 国家环境保护总局,国家质量监督检验检疫总局.地表水环境质量标准:G B3 8 3 82 0 0 2S.北京:中国环境科学出版社,2 0 0 2.1 4 罗平平,武阳,王双涛,等.沣河流域水质的时空对比分析J.水资源与水工程学报,2 0 2 1,3 2(5):3 5-4 1.1 5 沈珍瑶,谢彤芳.一种改进的灰关联分析方法及其在水环境质量评价中的应用J.水文,1 9 9 7,1 7(3):1 3-1 5.1 6 丁启振,周金龙,季彦桢,等.石河子昌吉地区地下水水质时空变化及 污 染源 解析 J.环 境 科学,2 0 2 3,4 4(3):1 4 4 0-1 4 5 1.1 7 M o o r eR,S p i t t l e h o u s eD,S t o r yA.R i p a r i a nm i c r o c l i-m a t ea n ds t r e a mt e m p e r a t u r e r e s p o n s e t o f o r e s t h a r v e s-t i n g:ar e v i e wJ.J o u r n a lo ft h e Am e r i c a n W a t e rR e s o u r c e sA s s o c i a t i o n,2 0 0 5,4 1(4):8 1 3-8 3 4.1 8 豆高飞.长江水热情势变化及其鱼类影响研究D.河南郑州:华北水利水电大学,2 0 2 1.1 9 黄炜惠,马春子,李文攀,等.我国地表水溶解氧时空变化及其对全球变暖的响应J.环境科学学报,2 0 2 1,4 1(5):1 9 7 0-1 9 8 0.2 0 杨艾琳,杨芳,李少斌,等.福建省流域近海溶解氧时空格局与低氧调控机制J.环境科学,2 0 2 2,4 3(1 1):4 9 5 0-4 9 6 0.2 1 岑言霸,苏斌,冯泽波,等.滇池入湖河流p H时空分异特征及其变化机理J.科学技术与工程,2 0 2 1,2 1(2 6):1 1 4 3 2-1 1 4 4 2.2 2 R a v i s hS,S e t i aB,D e s w a lS.G r o u n d w a t e rq u a l i t ya n a l y s i so f N o r t h e a s t e r n H a r y a n a u s i n g m u l t i v a r i a t es t a t i s t i c a lt e c h n i q u e sJ.J o u r n a lo ft h e G e o l o g i c a lS o c i e t yo f I n d i a,2 0 2 0,9 5(4):4 0 7-4 1 6.2 3 房丽晶,高瑞忠,贾德彬,等.草原流域地下水化学时空特征及环境驱动因素:以内蒙古巴拉格尔河流域为例J.中国环境科学,2 0 2 1,4 1(5):2 1 6 1-2 1 6 9.2 4 Z h a n gS h u a n g,Y a n gD o n g,L i F a n g l i n,e t a l.D e t e r m i-n a t i o no fr e g i o n a ls o i lg e o c h e m i c a lb a s e l i n e sf o rt r a c em e t a l s w i t h p r i n c i p a lc o m p o n e n tr e g r e s s i o n:ac a s es t u d y i nt h eJ i a n g h a np l a i n,C h i n aJ.A p p l i e dG e o-c h e m i s t r y,2 0 1 4,4 8:1 9 3-2 0 6.2 5 X uP a n p a n,Z h a n gQ i y i n g,Q i a nH u i,e ta l.C h a r a c t e r i-z a t i o no fg e o t h e r m a lw a t e ri nt h ep i e d m o n tr e g i o no fQ i n l i n gM o u n t a i n s a n dL a n t i a n-B a h eG r o u p i nG u a n z h o n g B a s i n,C h i n aJ.E n v i r o n m e n t a l E a