西明江流域水質(zhì)時空分異特征及污染源識別

覃楠鈞 趙侶璇 劉凱 宋曉薇 黃業(yè)翔

【摘? 要】論文應(yīng)用多種多元統(tǒng)計方法對廣西南寧市西明江流域水質(zhì)進行時空分異特征分析及污染源識別。結(jié)果表明，流域水質(zhì)采樣點可分為A（農(nóng)村生活用地）、B（城市生活用地）兩個組，其中，B組區(qū)域營養(yǎng)鹽污染及有機物污染水平均高于A組區(qū)域;全年可分為a（豐水期）、b（平、枯水期）兩個時段，a時段河流水質(zhì)相對較好，b時段流量較小，水質(zhì)易受河流沿程生活污水直排口污染;A組區(qū)域主要污染來源于面源污染，B組區(qū)域主要污染來源于生活污染及底泥內(nèi)源釋放。

【Abstract】This paper uses a variety of multivariate statistical methods to analyze the spatio-temporal variation characteristics of the water quality of the Ximingjiang basin in Nanning City， Guangxi Province and identify pollution sources. The results show that the water quality sampling points of the basin can be divided into two groups： A （rural living land） and B （urban living land）. Among them， the level of nutrient pollution and organic pollution in group B area is higher than that in group A area. The whole year can be divided into two periods： a （wet period） and b （normal and dry period）. The river water quality in period a is relatively good， and the flow in period b is relatively small， and the water quality is easily polluted by the direct discharge of domestic sewage along the river. The main pollution in group A area comes from non-point source pollution， and the main pollution in group B area comes from domestic pollution and sediment endogenous release.

【關(guān)鍵詞】多元統(tǒng)計;時空分異;污染源識別;水質(zhì)管理

【Keywords】multivariate statistics; spatio-temporal variation; identification of pollution sources; water quality management

【中圖分類號】X522? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?【文獻標(biāo)志碼】A? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?【文章編號】1673-1069（2020）09-0033-05

1 引言

西明江是流入邕江的一級支流，是廣西南寧市打造“中國水城”重點綜合治理的十八條內(nèi)河之一。在西明江未開展綜合整治前，長期承擔(dān)著“城市下水道”的功能，雨水、生活污水一同直排入河道，水質(zhì)污染嚴(yán)重，惡臭不堪，一直為劣Ⅴ類，嚴(yán)重制約著城市經(jīng)濟社會可持續(xù)發(fā)展和威脅著人民生命財產(chǎn)安全。西明江流域面積相對較小，但受到土地利用、人類活動方式的影響，流域水質(zhì)在時間和空間上仍呈現(xiàn)出較大的異質(zhì)性。水環(huán)境時空分布特征研究能夠為水質(zhì)評價提供動態(tài)信息，同時，這也是河流水質(zhì)監(jiān)測和管理的手段，正確識別流域水污染時空分異特性及污染源特征是科學(xué)、有效開展流域綜合治理和水生態(tài)環(huán)境監(jiān)管的基礎(chǔ)與前提[1，2]。目前，西明江僅在入邕江口處設(shè)有一個環(huán)保水質(zhì)監(jiān)測點，監(jiān)測數(shù)據(jù)有限，且流域水文水資源、氣候氣象、土地利用等數(shù)據(jù)存在不同程度的缺失，難以建立起高精度的水質(zhì)模擬模型。在此有限條件下，本文嘗試應(yīng)用多元統(tǒng)計分析方法分析西明江流域水質(zhì)時空分布特征。

多元統(tǒng)計分析方法主要包括主成分分析法[3]、聚類分析[4]、判別分析[5]、多維尺度分析[6]等方法，在河流、近海海域、湖泊和地下水水質(zhì)評價與分析中得到了很好的應(yīng)用。例如，陳靜生等[7]應(yīng)用多元統(tǒng)計法分析長江、黃河、松花江20世紀(jì)60年代到80年代水質(zhì)變化趨勢及其成因和機理;董慧峪等[8]應(yīng)用空間聚類分析方法得出制約南苕溪支流錦溪水質(zhì)改善的主要因素為農(nóng)業(yè)面源;田智慧等[9]利用方差分析對清潩河流域（許昌段）12個斷面進行時空間顯著差異性檢驗;韓向楠[10]通過Mann-Kendall趨勢檢驗、有序聚類分析等方法分析涇河流域多年來徑流、泥沙、降水的時空變化特征。但在以往的研究中多是分別進行水污染時間特征和空間特征分析，單獨探討時間或空間尺度的污染源特征，對流域綜合治理和水生態(tài)環(huán)境監(jiān)管指導(dǎo)意義有限。同時，過去的研究多以大河和湖泊水質(zhì)變化過程分析，對于小型河流的研究較少，在有限的數(shù)據(jù)資料中較為準(zhǔn)確掌握小流域水環(huán)境時空分異特征也是值得探討的問題。鑒于此，本研究應(yīng)用聚類分析、判別分析等方法對西明江水污染物時空變化特征進行分析，掌水污染物在時間、空間兩個尺度上的演變規(guī)律，應(yīng)用時空聯(lián)合分析、因子分析識別不同時間和區(qū)域的污染源特征，實現(xiàn)橫向和縱向?qū)Ρ?，以利于相關(guān)部門掌握流域水質(zhì)特征，有針對性地制定綜合治理方案。

2 研究區(qū)域

西明江位于邕江北岸，廣西南寧市西鄉(xiāng)塘區(qū)，發(fā)源于西鄉(xiāng)塘區(qū)大彎弓、小香爐嶺、大香爐嶺一帶，于羅文園藝場東側(cè)約500m處匯入邕江。天然河長13.9km，寬30～60m，深度7～13m，左右支流上游修建有羅文、禾屋等小水庫，控制流域面積33.6km2，河道枯水期水深0.1m左右，豐水期水深約1.5m，天然河道的常水位水域面積為27.5hm2。西明江月平均徑流量為Q=0.19m3/s（表1），河道徑流量在南寧市內(nèi)河水系中相對較小。西明江的功能定位是“邕江支流，城市排洪通道，作為城市內(nèi)的自然生態(tài)要素，結(jié)合水體形成防護林帶，與邕江交匯的河段形成公共、開敞的濱水景觀帶”，根據(jù)《南寧市水功能區(qū)劃》（南府〔2012〕107號），西明江流域執(zhí)行地表水Ⅲ類標(biāo)準(zhǔn)。

3 數(shù)據(jù)采集與研究方法

3.1 數(shù)據(jù)采集

為了對西明江流域水質(zhì)時空分布變化特征進行研究，沿西明江河段有大量廢水集中排入河段、沿岸主要居民密集區(qū)域、支流與干流匯合處的上下游等布設(shè)地表水質(zhì)采樣點位16個（圖1），其中于沿岸畜禽養(yǎng)殖區(qū)域、污水直排口和已判定為黑臭水體的河段加密布點。采樣點位避開死水區(qū)、回水區(qū)和排污口附近。各采樣點位結(jié)合水寬和水深的實際情況，按照HJ/T 91—2002《地表水和污水監(jiān)測技術(shù)規(guī)范》、HJ 494—2009《水質(zhì) 采樣技術(shù)指導(dǎo)》要求設(shè)置相應(yīng)的采樣垂線數(shù)和采樣點數(shù)。根據(jù)流域水文規(guī)律，分別在2018-2019年的豐水期（6～7月）、平水期（8～9月）和枯水期（2～3月）進行采樣。在前期調(diào)查中未發(fā)現(xiàn)有重金屬污染物，所以本次采樣監(jiān)測項目主要包括：溫度（T）、pH值、溶解氧（DO）、化學(xué)需氧量（CODGr）、五日化學(xué)需氧量（BOD5）、高錳酸鹽指數(shù)（CODMn）、氨氮（NH3-N）、總氮（TN）、總磷（TP）、硫化物（SULF）、懸浮物（SS）、電導(dǎo)率（SEC）、硫酸鹽（SO42-）、硝酸鹽（NO3-）共14項指標(biāo)，按照HJ/T 91—2002《地表水和污水監(jiān)測技術(shù)規(guī)范》規(guī)定的分析方法進行測定。西明江流域水質(zhì)采樣監(jiān)測數(shù)據(jù)統(tǒng)計描述見表2。

3.2 研究方法

3.2.1 聚類分析

聚類分析是一種研究“物以類聚”問題的多元統(tǒng)計方法，根據(jù)樣本（或研究對象）的數(shù)值特征，對其進行分類，判定各

樣本的親疏關(guān)系。本文所采用的層次聚類分析又叫作系統(tǒng)聚類分析，是目前使用最多的一種聚類方法。本次分析采用平方歐式距離度量樣本間的距離，運用離差平方和法（Ward法）生成具有層次結(jié)構(gòu)的聚類樹，分別在空間尺度和時間尺度上對西明江水質(zhì)進行分類分析[11-13]。

3.2.2 判別分析

判別分析是在已知樣本（或研究對象）分成若干類，并已知各類樣品觀測數(shù)據(jù)的基礎(chǔ)上，根據(jù)一定的準(zhǔn)則建立判別函數(shù)，然后根據(jù)判別函數(shù)對位置所屬類別的樣本進行分類的一種分析方法[14]。本文采用后退式逐步判別法對聚類分析結(jié)果進行驗證分析，相應(yīng)的判別函數(shù)表達式為：

式（1）中，i為組別G的序號，ki為各組的固有常數(shù)，n為各組用于分類的變量個數(shù)，Wj為特定變量（Pj）的權(quán)重系數(shù)，Pj為某變量的數(shù)值，f為判別函數(shù)值。

3.2.3 因子分析

因子分析是在以丟失最少信息為前提，將具有一定相關(guān)關(guān)系的多個變量綜合為少數(shù)幾個指標(biāo)（即因子）的分析方法[15]。本文在聚類分析和判別分析基礎(chǔ)上，采用因子分析逐步識別各個時間與空間組合的污染源特征，數(shù)據(jù)分析采用KMO和Bartlett球形檢驗進行適宜性檢驗，通過最大方差法進行旋轉(zhuǎn)，以特征值大于1為原則提取因子變量[16，17]。

3.2.4 數(shù)據(jù)處理方法

為了消除量綱的影響，在進行聚類分析和因子分析前，需對數(shù)據(jù)進行標(biāo)準(zhǔn)化預(yù)處理，主要采用Z得分法，即將各變量值減去均值后除以標(biāo)準(zhǔn)差，標(biāo)準(zhǔn)化后的Z分?jǐn)?shù)平均值為0，標(biāo)準(zhǔn)差為1。數(shù)據(jù)預(yù)處理和多元統(tǒng)計過程采用Microsoft Excel 2013和IBM公司的SPSS 24。

4 結(jié)果與分析

4.1 空間尺度聚類分析

由空間尺度聚類分析結(jié)果可知（圖2），當(dāng)20<聚類距離<25時，采樣點可分為2組;當(dāng)18<聚類距離<20時，采樣點可分為3組。應(yīng)用后退式判別分析分別對空間兩種分組方式的原始數(shù)據(jù)進行判別，其中2組分類的判別正確率為89.6%，3組分類的判別正確率為81.3%（表3）。因此，采樣點應(yīng)當(dāng)分為2組：A組{2、4、1、3、10、11、13、12};B組{5、8、14、15、6、9、16、7}。基本上是以大學(xué)西路為分界，以北區(qū)域為羅文村所在地，人類活動影響較小，主要受農(nóng)業(yè)面源、養(yǎng)殖源影響，A組解釋為農(nóng)村生活用地;以南為城市建設(shè)所在地，人類活動較為頻繁，主要受城鎮(zhèn)生活源、點源影響，B組解釋為城市生活用地。通過水質(zhì)指標(biāo)標(biāo)準(zhǔn)化箱式圖分析空間污染特征差異（圖3），可以看出A組DO、TN濃度明顯高于B組，而B組CODGr、BOD5、CODMn、NH3-N、TP、SS濃度明顯高于A組，即大學(xué)西路以南區(qū)域有機物污染水平均高于以北區(qū)域，建議重點加強監(jiān)測與管控。

4.2 時間尺度聚類分析

由時間尺度聚類分析結(jié)果可知（圖4），可分為兩個時段：a時段{6月、7月};b時段{10月、11月、2月、3月}，對比西明江徑流表（表1），分組情況基本符合流域水文規(guī)律，a時段可以解釋為豐水期，b時段可以解釋為平、枯水期。通過水質(zhì)指標(biāo)標(biāo)準(zhǔn)化箱式圖分析時間污染特征差異（圖5），b時段pH、CODGr、BOD5、CODMn、NH3-N、TN、TP濃度明顯高于a時段，結(jié)合表1可知在P=90%標(biāo)準(zhǔn)條件下b時段最枯季流量多在0.02～0.07m3/s，水體環(huán)境容量小，稀釋能力弱。根據(jù)前期污染源調(diào)查，西明江沿程污水直排口平均流量約為0.01m3/s，污水流量接近b時段最枯季流量，易造成水體污染。

4.3 因子分析

4.3.1 A組各時段因子分析

將A組原始數(shù)據(jù)分為A組-a時段、A組-b時段進行因子分析，識別不同時間段的污染源特征。由因子分析結(jié)果可知（表4），在A組-a時段抽取了3個因子，A組-b時段抽取了4個因子，累計方差分別為85.3%和81.2%。

在A組-a時段，因子1方差貢獻42.8%，特征因子為TN、TP;因子2方差貢獻30.6%，特征因子為CODMn;因子3方差貢獻11.9%，特征因子為BOD5。此類特征因子表明，污染主要來源于過量的營養(yǎng)鹽和有機物，A組區(qū)域沿河有大片蔬菜地和養(yǎng)鴨場，a時段雨水沖刷土壤中氮、磷進入水體以及畜禽養(yǎng)殖污染物排入河道引起耗氧有機物和營養(yǎng)鹽污染。

在A組-b時段，因子1方差貢獻28.3%，特征因子為DO、SEC;因子2方差貢獻25.5%，特征因子為CODMn、T;因子3方差貢獻15.3%，特征因子為SS;因子4方差貢獻12.1%，特征因子為TN。此類特征因子表明，b時段河流水量小，水中含氧量不足，溫度升高極易加速有機物和營養(yǎng)鹽消耗氧氣。

4.3.2 B組各時段因子分析

將B組原始數(shù)據(jù)分為B組-a時段、B組-b時段進行因子分析，識別不同時間段的污染源特征。由因子分析結(jié)果可知（表5），在B組-a時段抽取了4個因子，B組-b時段抽取了5個因子，累計方差分別為87.9%和85.7%。

在B組-a時段，因子1方差貢獻35.9%，特征因子為NH3-N、TN;因子2方差貢獻30.8%，特征因子為CODMn;因子3方差貢獻11.6%，特征因子為BOD5、SULF;因子4方差貢獻9.5%，特征因子為DO。此類特征因子表明，B組區(qū)域因左右兩側(cè)支流交匯且a時段河流流量大，因此，DO增加，但由于該區(qū)域有生活污水直排口匯入，導(dǎo)致水體有機污染物濃度升高。

在B組-b時段，因子1方差貢獻29.2%，特征因子為CODGr、CODMn、NH3-N、T;因子2方差貢獻19.5%，特征因子為SS、SEC;因子3方差貢獻16.4%，特征因子為pH;因子4方差貢獻11.4%，特征因子為BOD5、SULF;因子5方差貢獻9.3%，特征因子為TN。此類特征因子表明，B組-b時段流量小，來自左右兩側(cè)支流的底質(zhì)不斷堆積，此時不僅有生活污水帶來大量的有機污染物，還有底泥釋放的有機污染和氨氮。

4 結(jié)論

①在空間尺度上，西明江流域水質(zhì)采樣點分成A（農(nóng)村生活用地）、B（城市生活用地）兩個組。其中，大學(xué)西路以南的B組區(qū)域營養(yǎng)鹽污染及有機物污染水平均高于以北的A組區(qū)域，建議重點加強監(jiān)測與管控。

②在時間尺度上，全年分a（豐水期）、b（平、枯水期）兩個時段。其中，a時段河流水質(zhì)相對較好，b時段流量較小，河流沿程生活污水直排口流量接近該時段最枯季流量，易造成水體污染。

③綜合時空聯(lián)合分析與因子分析的結(jié)果，A組區(qū)域主要污染來源于面源污染，a時段雨水沖刷土壤以及畜禽養(yǎng)殖污染物排入河道引起耗氧有機物和營養(yǎng)鹽污染;b時段水中含氧量不足，溫度升高極易加速有機物和營養(yǎng)鹽消耗氧氣。B組區(qū)域主要污染來源于生活污染及底泥內(nèi)源釋放，a時段河流流量大，但生活污水直排口匯入導(dǎo)致水體有機污染物濃度升高;b時段不僅有生活污水帶來大量的有機污染物，還有底泥釋放的有機污染物和氨氮。

【參考文獻】

【1】伍冠星，楊永宇，劉暢.水環(huán)境因子時空分布特征研究進展[J].民營科技，2016，199（10）：21-23.

【2】王翠榆，楊永輝，周豐，等.沁河流域水體污染物時空分異特征及潛在污染源識別[J].環(huán)境科學(xué)學(xué)報，2012，32（9）：2267-2278.

【3】任澤，楊順益，汪興中，等.洱海流域水質(zhì)時空變化特征[J].生態(tài)與農(nóng)村環(huán)境學(xué)報，2011，27（4）：14-20.

【4】李一平，嚴(yán)瑩，韓廣毅.太湖水質(zhì)時空相關(guān)性分析[J].河海大學(xué)學(xué)報（自然科學(xué)版），2005（5）：505-508.

【5】郭蓉，吳長年，張寧紅，等.京杭運河（蘇州段）水質(zhì)時空分布特征分析[J].四川環(huán)境，2010，29（5）：10-13，29.

【6】李然然.查干湖濕地水環(huán)境演變及生態(tài)風(fēng)險評估[D].長春：中國科學(xué)院研究生院（東北地理與農(nóng)業(yè)生態(tài)研究所），2014.

【7】陳靜生，夏星輝，張利田，等.長江、黃河、松花江60—80年代水質(zhì)變化趨勢與社會經(jīng)濟發(fā)展的關(guān)系[J].環(huán)境科學(xué)學(xué)報，1999（5）：500-505.

【8】董慧峪，強志民，李庭剛，等.南苕溪支流錦溪水質(zhì)時空變化特征分析[J].環(huán)境工程學(xué)報，2012，6（3）：772-778.

【9】田智慧，張丹丹，赫曉慧，等.清潩河流域水質(zhì)時空分異與評價研究[J].人民長江，2017，48（7）：26-30.

【10】韓向楠.涇河流域水沙時空分異特征及其影響因素分析[D].重慶：西南大學(xué)，2019.

【11】任巖，張飛，王娟，等.新疆艾比湖流域地表水豐水期和枯水期水質(zhì)分異特征及污染源解析[J].湖泊科學(xué)，2017，29（5）：1143-1157.

【12】邸琰茗，白文榮，郭逍宇.密云水庫上游控制點水質(zhì)時空交叉統(tǒng)計分析[J].水資源開發(fā)與管理，2019（12）：21，27-34.

【13】薛薇.統(tǒng)計分析與SPSS的應(yīng)用[M].北京：中國人民大學(xué)出版社，2017.

【14】王昱，盧世國，馮起，等.黑河上中游水質(zhì)時空分異特征及污染源解析[J].中國環(huán)境科學(xué)，2019，39（10）：4194-4204.

【15】孫娜.表觀污染特征物質(zhì)及污染水體分類研究[D].蘇州：蘇州科技學(xué)院，2015.

【16】黃金良，黃亞玲，李青生，等.流域水質(zhì)時空分布特征及其影響因素初析[J].環(huán)境科學(xué)，2012，33（4）：1098-1107.

【17】魏彬.地表水表觀污染特性及其污染物分布規(guī)律研究[D].蘇州：蘇州科技大學(xué)，2019.