0 引言

河流是流域内工农业生产和生活的主要水源地，河流水质是影响人类生存和生态环境的重要因素，也是影响河口海湾水质的主要因素^[1-3]。南流江是广西桂东南诸河中最大的独流入海河流^[4]。随着工农业、畜禽养殖业、城市化的迅速发展，工农业污水、养殖废水和生活污水排放量增加，导致南流江水质严重恶化^[5]，治理南流江污染成为广西跨世纪绿色工程之一^[6]。鉴于此，研究近年来南流江的水质空间分布特征及其变化情况，对南流江的科学治理及其治理效果评价具有十分重要的意义。

目前已有一些研究对南流江水质进行了分析评价，使人们初步认识到南流江水质的总体分布及21世纪初的变化。代俊峰等^{[6, 7]}分析2003-2009年南流江下游常乐站和总江口的水质数据，发现南流江水质恶化的影响因素主要是禽畜养殖业、农业等面源污染。李昆明等^[8]以南流江玉林段为例，分析南流江上游玉林段的水质。林卉等^{[9, 10]}分析南流江干流亚桥、南域、江口大桥、横塘和六司桥5个断面，发现南流江下游的污染程度低于中上游，氨氮、总磷和五日生化需氧量是主要污染指标。已有的研究主要集中在南流江干流或者上游、下游部分干流断面，但对南流江包括各个支流在内的全流域水质缺乏整体性研究，因此无法全面了解南流江水质空间分布状况。

本研究利用2016年、2019年和2020年共3年的南流江干流和各支流超过60个断面的监测数据，结合单因子和综合水质标识指数法，研究南流江全流域水质空间分布，分析自2016年实施综合治理以来南流江水质的年际变化特征，探讨南流江水质时空变化原因，以期为南流江流域的精细化管理、治理和评估等工作提供支撑。

1 材料与方法 1.1 研究区域概况

南流江发源于广西北流市大容山，于北海党江镇注入廉州湾，是广西独流入海第一大河，其流域(图 1)位于广西东南部(108°51′-110°22′E，21°34′-22°52′N)^[11]。南流江流经玉州区、北流市、博白县、浦北县等10县(市、区)，干流全长287 km，流域面积为9 232 km^{2 [12, 13]}。南流江流域是广西农业、禽畜养殖业、轻工业较为发达的地区，而且近年来该流域工业发展迅速，对广西经济增长的贡献率明显提升^[13]。

1.2 数据来源

针对南流江以面源污染为主的特点，结合南流江综合治理的实践效果，本研究选取2016年5月、2019年5月和2020年5月的水质监测数据，分析南流江综合治理前后的变化。南流江干流和支流的监测断面布设如图 1所示，其中2016年5月布设42个，用编号“A1-A42”表示; 2019年5月和2020年5月均布设61个，用编号“B1-B61”表示。监测数据由广西壮族自治区海洋环境监测中心站和广西壮族自治区生态环境监测中心提供。

1.3 水质评价方法

采用单因子水质标识指数法和综合水质标识指数法对南流江水质进行评价。根据前期研究成果^[14]，选择南流江的主要超标(Ⅲ类)因子进行分析，包括高锰酸盐指数、五日生化需氧量、氨氮和总磷4项指标，结合水质标识指数，利用ArcGIS将南流江水质状况进行空间化，分析其空间变化和年际变化特征。

单因子水质标识指数既能按照《地表水环境质量标准(GB 3838-2002)》^[15]进行水质类别的定性评价，又能根据标识指数在同一类水质指标中进行定量评价^[16]。单因子水质标识指数法由一位数字和小数点后2位小数组成，表达式为^{[17, 18]}

$ P = {X_1}.{X_2}{X_3}, $

(1)

$ {X_2} = \;\frac{{{C_i} - {C_{ik下}}}}{{{C_{ik上}} - {C_{ik下}}}}\; \times 10。$

(2)

式(1)中：P为单因子水质标识指数；X₁是某项水质指标的水质类别，可与相应国家水质标准指数比较确定，当某项水质指标水质类别为Ⅰ类时，X₁=1，当某项水质指标水质类别为Ⅱ类时，X₁=2，…，以此类推；X₂为水质监测数据在X₁类水质变化区间中所处的位置；X₃为水质类别与水质功能区划类别(《地表水环境质量标准(GB 3838-2002)》)^[15]Ⅲ类标准比较的结果。式(2)中：C_i是第i项指标的实测浓度，mg/L；C_ik上、C_ik下分别为某项水质指标在第k类水质区间的上限值和下限值，mg/L。单因子水质指标判断关系如表 1所示。单因子水质标识指数法的具体计算方法见参考文献[17, 18]。

综合水质标识指数法是基于单因子水质标识指数法，结合河流和监测点总体的水质信息、水环境功能区等相关内容进行合理评价的水质评价方法^{[19, 20]}。综合水质标识指数法由一位整数和小数点后3位小数组成，其表达式为^{[19, 20]}

$ I = {X_1}.{X_2}{X_3}{X_4}, $

(3)

$ {X_1}.{X_2} = \frac{1}{m}\sum\limits_{i = 1}^m {\left( {{P_1} + {P_2} + \cdots + {P_m}} \right)} 。$

(4)

式(3)中：X₁.X₂为综合水质类别；X₃是参评的单项水质指标中比水体功能区目标差的水质指标个数；X₄是综合水质类别与水体功能区类别的比较结果，如果X₄=0，说明综合水质类别达到或者优于功能区，如果X₄=1，说明综合水质类别比功能区差1个类别，以此类推。按照南流江水质功能区《地表水环境质量标准(GB 3838-2002)》^[15]Ⅲ类标准进行比较评价。式(4)中：m表示参与综合水质评价的单项指标个数；P₁、P₂、…、P_m分别为第1，2，…，m个水质因子的单因子水质标识指数，对应的是单因子水质标识指数中的X₁.X₂。X₁.X₂可以判断综合水质类别，具体判断标准如表 2所示。综合水质标识指数法的具体计算方法见参考文献[19, 20]。

2 结果与分析 2.1 单因子水质标识指数空间变化

根据单因子水质标识指数(以下称“单因子指数”)结果显示(图 2)，2016年4个指标的单因子指数空间变化显著，空间变化规律各异(2019年和2020年部分断面的高锰酸盐指数和五日生化需氧量未监测)。高锰酸盐指数在中上游出现劣Ⅴ类、Ⅳ类水质较多，下游以Ⅲ类、Ⅱ类为主。五日生化需氧量在整个流域内以Ⅳ类、Ⅲ类为主。氨氮和总磷指数在中游、上游局部出现劣Ⅴ类、Ⅳ类水质较多，下游以Ⅰ类和Ⅱ类为主，但在水鸣河、马江、张黄江、武利江的部分监测断面中出现劣Ⅴ类、Ⅳ类。

2019年的高锰酸盐指数、五日生化需氧量指数空间差异不明显，氨氮、总磷指数的单因子指数空间差异较大。高锰酸盐指数和五日生化需氧量指数均达到或优于南流江Ⅲ类水质功能区类别。氨氮指数以Ⅲ类、Ⅱ类为主，上游局部出现Ⅳ类，中游博白县部分监测断面出现劣Ⅴ类，下游北海石康镇断面出现劣Ⅴ类。总磷指数由上游的Ⅲ类逐渐向中游Ⅳ类、Ⅴ类变化，再逐渐向下游变为Ⅲ类。

2020年的高锰酸盐指数、五日生化需氧量指数差异不显著，但均达到南流江Ⅲ类水质功能区类别。氨氮指数在中上游主要以Ⅲ类、Ⅱ类为主，在兴业县、玉州区、博白县部分断面出现Ⅳ类、Ⅴ类，下游主要以Ⅰ类、Ⅱ类为主，在石康镇断面出现了Ⅴ类。总磷指数则由上游的Ⅲ类逐渐向中游增加为Ⅳ类、Ⅴ类，再逐渐向下游减少至Ⅲ类、Ⅱ类。

2.2 综合水质标识指数空间变化

综合水质标识指数空间分布结果(图 3)显示，2016年南流江总体水质类别为Ⅳ类，主要分布在北流市、玉州区、福绵区、陆川县。Ⅴ类水质分布在兴业县、玉州区、福绵区、博白县。劣Ⅴ类水质则出现在玉州区，Ⅰ类和Ⅱ类水质则主要分布在浦北县、灵山县和合浦县。2019年和2020年南流江总体水质类别为Ⅲ类，2019年和2020年Ⅳ类水质主要分布在兴业县、玉州区、博白县，Ⅴ类水质主要分布在福绵区、博白县。Ⅰ类和Ⅱ类水质则主要集中在中下游的博白县县城往下的流域以及灵山县、合浦县。与2016年相比，2019年和2020年玉林市境内的综合水质类别为Ⅰ-Ⅲ类的水质有所增加。

2.3 单因子标识指数年际变化

由于2019年5月和2020年5月未对南流江支流部分监测断面的高锰酸盐指数和五日生化需氧量进行监测，因此下文只分析氨氮和总磷的单因子水质标识指数的年际变化。

干流中氨氮、总磷年际变化明显，与2016年相比，2019年、2020年氨氮分别下降27.78%、30.55%，总磷分别下降27.78%、25%。中上游支流中氨氮、总磷年际变幅较大，与2016年相比，2019年、2020年氨氮都下降31.75%，总磷分别下降25.40%、44.44%；下游各支流中，氨氮、总磷的年际变幅不明显(图 2，4)。

从图 2和图 4来看，2016年南流江干流中，北流市、博白县、玉州区、福绵区的大部分监测断面氨氮浓度超标，超标断面单因子指数为4.0-8.0，下游合浦县单因子指数为1.5-2.5。2019年和2020年玉林市氨氮单因子指数有所下降，特别是北流市、玉州区、福绵区，但是博白县部分断面氨氮单因子指数仍较高。2016年干流中，北流市、福绵区和玉州区监测断面总磷浓度超过南流江Ⅲ类水质功能区类别，博白县超过一半的监测断面超过南流江Ⅲ类水质功能区类别，5个区(县)超标断面单因子指数为4.5-8.5。2019年和2020年，福绵区、北流市、玉州区、博白县、合浦县总磷指数均有所下降，但仍有些监测断面总磷浓度超过南流江Ⅲ类水质功能区类别，单因子指数为4.5-8.5。

氨氮浓度在上游支流中定川江流域下降明显，2016年单因子指数为3.5-7.5，2019年和2020年为1.0-4.0。清湾江流域氨氮浓度下降明显，2016年单因子指数为4.5-6.0，2019年为1.0-5.0，到2020年转为2.0-3.5。仁东河流域氨氮浓度有所下降，2016年单因子指数为6.5-8.5，到2019年和2020年为4.5-6.5，但氨氮含量依然超过南流江Ⅲ类水质功能区类别。在所调查的3年里，沙生江流域、丽江流域、沙田河流域的氨氮浓度均优于水质功能区类别，单因子指数均小于4.0。在中游支流中，绿珠江流域在这3年中的氨氮浓度均优于水质功能区类别。水鸣河流域这3年氨氮浓度达到Ⅲ类水质功能区类别。鸦山江流域氨氮单因子指数变动较大，由2016年的5.0-6.0变为2019年的3.5-8.0，再变为2020年的5.5-6.0。合江流域氨氮单因子指数由2016年的4.0变为2019年和2020年的2.5-4.0。在下游支流中，武利江流域、张黄江流域、洪潮江流域3年的氨氮浓度优于水质功能区类别，单因子指数均为1.5-3.0。

总磷浓度在上游支流中定川江流域显著下降，2016年单因子指数为4.5-7.5，2019年和2020年为3.0-4.5。清湾江流域总磷浓度下降明显，2016年单因子指数为5.0-6.0，2019年为1.0-5.0，到2020年转为2.0-3.5。仁东河流域总磷浓度超标，2016年单因子指数为6.0-9.0，到2019年和2020年为4.3-6.0。沙生江流域总磷单因子指数较小，水质类别均达到水质功能区类别。沙田河流域总磷浓度下降明显，由2016年的4.5-5.5转为2019年和2020年的3.0-4.0。在中游支流中，绿珠江流域总磷浓度变化反复，由2016的3.6变为2019年的4.3，到2020年转为2.7。水鸣河流域、鸦山江流域、合江流域总磷浓度超标，3年单因子指数均大于4.5。在下游支流中，利江流域、张黄江流域、洪潮江流域在3年中的氨氮浓度均优于Ⅲ类水质功能区类别，单因子指数均小于4.0。

2.4 综合水质标识指数年际变化

根据综合水质标识指数结果(图 3，5)分析：与2016年相比，2019年和2020年南流江水质总体趋于改善，Ⅰ-Ⅲ类水质比例分别上升21.98%、25.25%。

2016年南流江干流区域的北流市、福绵区、玉州区和博白县部分监测断面水质较差(I＞4)，下游的合浦县、浦北县总体水质较好(I≤4)。到2019年和2020年，北流市、玉州区、福绵区的水质均有所好转。虽然博白县部分监测断面水质有所好转，但部分断面水质仍较差，主要集中在博白镇、水鸣镇、东平镇、菱角镇、亚山镇。合浦县除了石康镇水质较差(I≥6)外，其他水质监测断面保持良好。

在上游支流中，定川江流域在2016年总体水质较差(I＞4)，而在2019年和2020年总体水质呈现好转趋势(I≤3)。仁东河流域在2016年的水质差(I＞5)，2019年和2020年总体水质得以改善(I＞4)。清湾江流域和丽江流域总体水质呈现好转趋势，水质由Ⅳ类逐渐转为Ⅱ类和Ⅲ类。沙生江流域和沙田河流域水质在2016年、2019年和2020年保持良好，均在Ⅲ类及以上(I≤3)。

在中游支流中，马江流域水质较好，绿珠江流域水质在这3年里均达到Ⅲ类及以上。鸦山江流域水质污染加重，由2016年的Ⅳ类逐渐转为2019年和2020年的劣Ⅴ类(I≥6)。水鸣河流域水质呈现好转趋势，由Ⅴ类逐渐转为Ⅳ类。合江流域水质状况反复，2016年水质类别为Ⅳ类，2019年达到水质功能区要求，2020年水质又恶化为Ⅳ类。

在下游的支流中，张黄江流域、洪潮江流域水质保持良好，均达到Ⅲ类及以上(I≤4)。武利江流域水质呈现好转趋势，2016年有部分监测断面水质为Ⅳ类，到2019年和2020年水质均达到Ⅱ类。

3 讨论

南流江水质较差的区域集中在玉林市玉州区、福绵区、博白县，主要污染物为氨氮、总磷。根据已有的研究可知，禽畜养殖污染、城镇及农村生活污染是南流江水质污染的主要来源^[21]。南流江流域存在大量禽畜养殖，是广西生猪养殖业的集中地之一，2016年南流江流域生猪出栏838.09万头，主要分布在玉州区、福绵区、博白县、浦北县、合浦县等地^{[14, 22]}。但由于养殖模式粗放，流域环境污染防治管理水平较差等原因，南流江流域养殖业污染问题越来越突出。根据监测数据显示，2015年南流江生猪养殖排放的氨氮、总磷分别为673.29，256.32 t^[23]，这些养殖废水直排南流江，导致南流江氨氮、总磷超标。此外，南流江流域人口众多，根据研究显示，2015年玉林市、钦州市、北海市城镇生活污水排放量分别为6 428.68，635.11，3 025.94 t，农村生活污水排放量分别为4 055.60，1 899.24，813.57 t^[22]。然而，南流江流域城镇和农村生活污水处理设施建设滞后，无法有效处理日益增多的生活污水，甚至部分生活污水未经处理直排入江。因此，生活污水及禽畜养殖业污水的排放是南流江氨氮、总磷的主要来源^[9]。

2016年、2019年、2020年3年里南流江水质总体趋于改善，玉林市范围南流江干支流中氨氮、总磷年际变化明显。这可能与政府对南流江流域综合整治相关。近年来，南流江水质污染较为严重，甚至因整治不力被生态环境部约谈。为综合治理南流江流域水环境，广西各级政府针对南流江实际情况出台污染防治方案，制定计划和措施，2016-2018年全区共投入约28.6亿元。例如，为治理禽畜养殖污染，玉林市完成畜禽养殖场截污池、生态养殖场等建设；博白县截污建池，收运还田^[24]。截至2018年，玉林市完成41个镇级污水处理厂的建设和设备安装，用以治理生活污水，北海市建成6个乡镇污水处理厂并投入试运行。为治理工业废水，玉林市7个工业园区已建成并投入运行^[24]。2020年5月南流江流域国控断面和区控断面水质均达Ⅲ类，说明综合治理已见成效。然而，在南流江水质总体好转的情况下，依然还有部分地区的水质严重超标，例如2019年在玉林市博白县博白镇、亚山镇，北海市合浦县石康镇出现劣Ⅴ类水质；玉林市福绵区福绵镇和玉州区南江街道出现Ⅴ类水质。2020年在博白县博白镇、东平镇、菱角镇，合浦县石康镇出现劣Ⅴ类水质，意味着还要继续推进南流江“两年显著见效，三年大见成效”的工作目标。

值得注意的是，南流江水质变好可能受生猪数量减少等因素影响。与2016年相比，2019年玉林市出栏生猪460.03万头，同比下降21.9%^[25]。南流江生猪养殖数量的下降，在一定程度上从源头减少了禽畜养殖污染。

4 结论

本研究利用水质标识指数法，分析南流江2016年、2019年和2020年水质空间分布及年际变化特征，并探讨其原因。主要结论如下：

① 南流江水质较差的区域集中在中上游的玉州区、福绵区、博白县，主要污染因子为氨氮、总磷，这与南流江流域的禽畜养殖业及生活污水排放有关。

② 经过综合整治，与2016年相比，2019年和2020年南流江流域的水质趋于改善。

③ 本文对南流江包括各个支流在内的全流域水质进行空间分布特征及年际变化情况分析，可为下一步南流江精细化治理提供科学依据。