代俊峰　楊藝　王璐瑜　王媛

(1.桂林理工大學(xué)廣西環(huán)境污染控制理論與技術(shù)重點實驗室　廣西桂林 541004；2.桂林理工大學(xué)廣西巖溶地區(qū)水污染控制與用水安全保障協(xié)同創(chuàng)新中心　廣西桂林 541004)

?

河流監(jiān)測斷面幾種污染物的點源、非點源污染負(fù)荷分割*

代俊峰1,2楊藝1王璐瑜1王媛1

(1.桂林理工大學(xué)廣西環(huán)境污染控制理論與技術(shù)重點實驗室廣西桂林 541004；2.桂林理工大學(xué)廣西巖溶地區(qū)水污染控制與用水安全保障協(xié)同創(chuàng)新中心廣西桂林 541004)

以廣西入海河流南流江2003—2011年水質(zhì)數(shù)據(jù)為例，通過水文估算法和數(shù)字濾波法對點源和非點源污染負(fù)荷分別進(jìn)行分割估算。結(jié)果表明，2003—2011南流江流域監(jiān)測斷面的非點源污染較為嚴(yán)重，采用水文估算法計算的多年平均硫酸鹽、氯化物、硝酸鹽氮、亞硝酸鹽氮的非點源污染負(fù)荷占總負(fù)荷的比例分別為0.82，0.81，0.82，0.72；基于數(shù)字濾波法計算的多年平均硫酸鹽、氯化物、硝酸鹽氮、亞硝酸鹽氮的非點源污染負(fù)荷占總負(fù)荷的比例分別為0.72，0.54，0.68，0.63。結(jié)果顯示，這兩種方法的污染負(fù)荷分割結(jié)果變化趨勢大致相同。

非點源污染點源污染數(shù)字濾波水文估算法

0　引言

近年來對非點源污染的控制成為污染治理的關(guān)鍵和難點之一，河口、湖泊、水庫富營養(yǎng)化的原因主要是由于非點源污染[1]。非點源污染的發(fā)生在時空上具有隨機性、發(fā)生方式具有間斷性、遷移過程具有繁雜性、排放途徑及排放量具有不確定性、監(jiān)測模擬與控制具有困難性等特點。針對非點源污染的上述特征，許多專家對于非點源量化問題也進(jìn)行了相關(guān)的探討與研究，且根據(jù)實際資料與現(xiàn)有條件提出了相對適合的非點源污染負(fù)荷估算方法，例如：水量水質(zhì)類方法、平均濃度法、降雨量差值法、徑流分割法、模型法等[2]。

廣西北部灣經(jīng)濟區(qū)水資源豐富，但水資源時空分布十分不均，非點源污染嚴(yán)重，河流都受到了不同程度污染。本文以南流江為例，分別運用水文估算法和數(shù)字濾波法對流域內(nèi)2003—2011年以硫酸鹽、氯化物、硝酸鹽氮、亞硝酸鹽氮為例的污染負(fù)荷進(jìn)行分割，定量描述點源污染和非點源污染對水生態(tài)環(huán)境的影響。

1　研究區(qū)域概況

南流江干流長285 km，平均坡降為0.035%，多年平均入海徑流量74.96億m3，流域內(nèi)地形起伏平緩，但其河道蜿蜒，河水流動速度較慢，自我凈化水平低，生態(tài)水環(huán)境極其脆弱。

筆者收集了南流江常樂水文站2003—2011年的系列水量資料和部分水質(zhì)監(jiān)測資料，水量資料和水質(zhì)資料來自水文站和水環(huán)境監(jiān)測站。南流江和常樂水文站位置見圖1。

圖1南流江和常樂水文站位置示意圖

2　點源、非點源污染物負(fù)荷的估算方法

2.1水文估算法

水文估算法是依據(jù)水文學(xué)原理，根據(jù)點源和非點源污染的產(chǎn)生、運動與轉(zhuǎn)移搬運規(guī)律，在已有的監(jiān)測資料及數(shù)據(jù)前提下來劃定、區(qū)分點源與非點源污染負(fù)荷量，并對非點源污染負(fù)荷量實行量化[3-5]，為進(jìn)一步了解、剖析水環(huán)境中生命活動及其污染現(xiàn)象與源頭提供參考和依據(jù)。

點源污染排入水環(huán)境的污染物數(shù)量相對比較平穩(wěn)固定，河川基流量Qu可通過解析徑流的形成以及運用直線斜割法直接從年徑流過程中獲得，且依據(jù)已測得的降水?dāng)?shù)據(jù)，采用枯水期的污染物濃度，對其進(jìn)行加權(quán)平均，以此作為點源污染物的平均濃度Cu，根據(jù)公式推求。非點源污染負(fù)荷可分為2個部分，一部分是因降雨徑流溶解在水中隨地表徑流匯入水體的可溶性污染物，另一部分是侵蝕土壤及其吸附的固態(tài)污染物，因侵蝕被運送入水體。

點源污染負(fù)荷推算采用式(1)。

Wu=Qu×Cu×Δtu

(1)

式中，Wu為點源污染負(fù)荷量；Qu為河川徑流量中地下水滲透補給河水的水量，m3/s；Cu為污染物的流量加權(quán)平均濃度，mg/L；Δtu為第u時段的時間，s。

由于南流江的懸移質(zhì)輸沙量相對較小，且缺乏連續(xù)性的河流懸移質(zhì)泥沙資料。所以，本文非點源污染負(fù)荷量的估算不考慮暴雨徑流中侵蝕土壤所吸附的污染物濃度負(fù)荷量，僅考慮地表徑流過程中溶解態(tài)污染物的非點源負(fù)荷。

實際估算中，非點源污染物負(fù)荷計算采用式(2)。

Wg=Cg×Qg×Δtg

(2)式中，Wg為非點源污染負(fù)荷量；Qg為第g時段的平均流量，m3/s；Δtg為相對應(yīng)時間段的時間，s；Cg為污染物在地表徑流過程中的平均濃度，推算采用式(3)。

(3)

式中，Ci為第i次監(jiān)控測得的污染物質(zhì)量濃度，mg/L；Qi為第i次監(jiān)控檢測獲得的流量，m3/s；n為每年監(jiān)測污染物濃度的次數(shù)。

2.2數(shù)字濾波法

2.2.1濾波原理

數(shù)字濾波法以傅立葉分析為基本理論依據(jù)，而傅立葉分析的關(guān)鍵又在于傅立葉變換。換而言之，傅立葉變換創(chuàng)設(shè)了以頻率為自變量的頻譜函數(shù)和以時間為自變量的時間函數(shù)之間的轉(zhuǎn)換關(guān)系。一般點源污染負(fù)荷在某流域出口斷面的出流時間序列是相對平穩(wěn)的，對其出流時間序列進(jìn)行傅立葉變換后，點源污染負(fù)荷的出流時間序列頻域下的信號較弱，主要是低頻信號。而非點源污染負(fù)荷的出流時間序列易受降水、強降雨等外界條件影響，因此變化幅度較大，對其出流時間序列進(jìn)行傅立葉變換后，非點源污染負(fù)荷的出流時間序列頻域下的信號較強，主要是高頻信號，顯然強于點源污染負(fù)荷出流的信號。以上表明，對于點源、非點源污染負(fù)荷的分割，能夠通過分離高低頻信號來實現(xiàn)[6]。

實際在傅里葉變換的眾多形式中，常采用離散傅立葉變換(DFT)，是由水文時間序列出現(xiàn)形式?jīng)Q定的(常以離散形式給出)。它針對的主要是周期序列或有限長序列存在的，離散傅立葉變換的一般形式見式(4)。

(4)

式中，N為周期序列一個周期的抽樣點數(shù)，或有限長序列(時域及頻域)的抽樣點數(shù)，X(k)為{x(n)}的離散傅立葉變換。

2.2.2水文學(xué)的實際應(yīng)用

NATHAN R J和MCMAHON T A首次采用Lyne-Hollick算法對流量過程進(jìn)行解析和運算[7]，從流量過程中分割出基流。濾波法較為客觀、易掌握、實施簡便、速率快，且參數(shù)較少。濾波方程為式(5)。

(5)

式中，qt表示t時刻過濾出的直接徑流信號(以d為時間步長) ；Qt，Qt-1表示相應(yīng)時段的實測河川的總徑流量；β為濾波參數(shù)。從總徑流中過濾出快速響應(yīng)，即可得出基流bt。

bt=Qt-qt

(6)

由式(5)和式(6)可計算出點源、非點源污染負(fù)荷的出流過程，進(jìn)而根據(jù)污染負(fù)荷的出流過程線對出流污染負(fù)荷進(jìn)行分割。

2.2.3數(shù)字濾波應(yīng)用分析

理論上說，假設(shè)β1、β2為兩個不同的濾波參數(shù)，進(jìn)行濾波多次，它們得出的結(jié)果近似。而在實際應(yīng)用中卻發(fā)現(xiàn)，如果選取相對較大濾波參數(shù)時，濾波進(jìn)行2～3次就可得出合理的點源、非點源污染負(fù)荷分割量。但如果選取的是相對較小濾波參數(shù)時，則需要進(jìn)行濾波 5～6次才能夠得到合理的污染負(fù)荷分割量。黎坤等以東江流域為例[8]，運用數(shù)字濾波法對點源、非點源污染負(fù)荷進(jìn)行分割，最后和平均濃度法計算結(jié)果進(jìn)行對比并得出結(jié)論，結(jié)論為采用比較大的濾波參數(shù)(β=0.925～0.950)，進(jìn)行3次濾波，得到的分割結(jié)果與平均濃度法計算的平均結(jié)果近似。由以上分析，本文取濾波參數(shù)β為0.925。

3　實例應(yīng)用的結(jié)果與分析

3.1水文估算法的實例應(yīng)用

本文以南流江常樂水文站2003—2011年的月、年平均流量統(tǒng)計表為基礎(chǔ)，采用枯水季最小月平均流量法，選擇年內(nèi)最小月的平均流量乘上全年總時間得到全年總基流量，而后用全年總徑流量減去基流量得出地表徑流量。南流江年徑流的分割見表1。

表1　南流江年徑流的分割　億m3

得到基流分割結(jié)果之后，根據(jù)已經(jīng)獲得的數(shù)據(jù)，用河川基流量乘上監(jiān)測的流量最小月污染物濃度求得點源污染負(fù)荷；計算求得非點源污染在地表徑流的流量加權(quán)平均濃度，得出的平均濃度再乘上地表徑流量得到非點源污染負(fù)荷。

3.2數(shù)字濾波法的實例應(yīng)用

分別將硫酸鹽、氯化物、硝酸鹽氮和亞硝酸鹽氮的數(shù)據(jù)存儲為txt文件，根據(jù)Baseflow Filter Program 的要求處理數(shù)據(jù)，然后運行bflow.exe[9-10]，得到的數(shù)據(jù)即為點源污染負(fù)荷量和指數(shù)Fr(點源污染負(fù)荷占總污染負(fù)荷比值的平均值)，并依據(jù)實際情況選取與水文估算法中最相近的Fr值。

3.3兩種方法分割結(jié)果對比

水文估算法計算出的南流江多年平均硫酸鹽、氯化物、硝酸鹽氮、亞硝酸鹽氮的非點源污染負(fù)荷占總負(fù)荷的比例分別為0.82，0.81，0.82，0.72；數(shù)字濾波法的計算結(jié)果分別為0.72，0.54，0.68，0.63。

水文估算法是一種簡單、易用、易理解的方法，利用常規(guī)的河道斷面監(jiān)測資料對污染物濃度、流量等進(jìn)行估算，從而求出點源、非點源污染物負(fù)荷。相比較而言，數(shù)字濾波法是根據(jù)計算機編程直接得出點源污染負(fù)荷量以及指數(shù)Fr，最后再根據(jù)已知量求得非點源污染負(fù)荷量，其操作更加便捷快速且易修改，在一定程度上簡化了計算的復(fù)雜性。兩種方法的具體對比結(jié)果見圖2、圖3、圖4和圖5。

圖2　兩種方法硫酸鹽污染負(fù)荷估算比較

圖3　兩種方法氯化物污染負(fù)荷估算比較

圖4　兩種方法硝酸鹽氮污染負(fù)荷估算比較

圖5　兩種方法亞硝酸鹽氮污染負(fù)荷估算比較

對比分析結(jié)果如下。

(1)數(shù)字濾波法和水文估算法的污染負(fù)荷分割結(jié)果變化趨勢大致相同，且近年總體污染得到一定控制。

(2)點源污染負(fù)荷的出流時間序列相對穩(wěn)定，波動較?。环屈c源污染負(fù)荷由于易受降水、強降雨等外界因素影響，從而變化幅度比較明顯，波動較大。

(3)南流江流域監(jiān)測斷面的非點源污染較為嚴(yán)重。

4　結(jié)語

(1)水文估算法不僅經(jīng)濟適用，而且可以估算不同年、某些特殊年份(如枯水年) 、不同頻率代表年、次洪水的非點源污染負(fù)荷量。與此同時，數(shù)字濾波法因具有簡便和可適用性強等優(yōu)點，能夠再現(xiàn)分割結(jié)果，讓使用者能在分割污染負(fù)荷時自己根據(jù)以往的經(jīng)驗靈活地選取濾波參數(shù)。

在可行性上，水文估算法已經(jīng)應(yīng)用于多地區(qū)研究污染負(fù)荷的實際問題中，使用較為廣泛。數(shù)字濾波法由于是近年結(jié)合計算機編程的新型研究方法，仍處于嘗試狀態(tài)，需要與其他方法得出的結(jié)論進(jìn)行對比，判斷其合理性，最后確定濾波參數(shù)。

(2)本文數(shù)字濾波法研究分析的只是污染物濃度大小的變化，且在頻譜分析中，點源污染負(fù)荷信號以低頻信號為主，非點源污染負(fù)荷的信號以高頻信號為主。而在對濾波實際操作中，雖不可防止相位失真，但對其進(jìn)行多次濾波后可在一定程度上進(jìn)行修正。而運用水文估算法進(jìn)行數(shù)據(jù)處理時涉及污染物濃度、徑流流量等多重因素影響，所以相對數(shù)字濾波法來說可能誤差較大。通過對污染負(fù)荷分割結(jié)果對比圖的分析可以看出，兩種方法的污染物負(fù)荷估算結(jié)果整體變化趨勢相近。其中，有個別年份出現(xiàn)例外，有相對較大的差異，如2006年硫酸鹽和氯化物的分割結(jié)果，兩種計算方法估算其非點源污染負(fù)荷量所占比例分別相差了0.44、0.36，這種情況的產(chǎn)生，多是由于當(dāng)年強降雨次數(shù)或者降雨強度等外界條件造成的。

(3)上述污染物負(fù)荷估算結(jié)果對比分析說明，廣西入海河流南流江流域監(jiān)測斷面的非點源污染較為嚴(yán)重，造成了不可忽視的環(huán)境污染與環(huán)境安全問題。在今后的污染防治和管理方面，應(yīng)立法明確控制措施，制定有效的管理規(guī)劃，加強對非點源污染的認(rèn)識，重視流域的整體管理。

[1]ARHONDITSIS G, TSIRTSIS G, ANGELIDIS M O, et al. Quantification of the effects of non-point nutrient sources to coastal marine eutrophication: application to a sem-inclosed gulf in the Mediterranean Sea[J]. Ecological Modelling, 2000,129(2):209-227.

[2]樊在義，宋兵魁，楊勇，等. 非點源污染負(fù)荷估算方法探討[J]. 環(huán)境科學(xué)導(dǎo)刊，2011,30(3)：1-6.

[3]李懷恩. 估算非點源負(fù)荷的平均濃度法及其應(yīng)用[J]. 環(huán)境科學(xué)學(xué)報,2000，20(4)：397-400.

[4]施為光, 凌文州. 用實測資料計算流域非點源負(fù)荷[J]. 長江流域資源與環(huán)境，1996，5(3)：273-277.

[5]鮑全盛, 曹利軍. 密云水庫非點源污染評價研究[J].水資源保護, 1997(1)：8-11.

[6]喬繼平,代俊峰. 河流污染的點源和非點源負(fù)荷分割研究[J]．中國農(nóng)村水利水電，2015(6):17-20.

[7]NATHAN R J，MCMAHON T A． Evaluation of automated techniques for baseflow and recession analyses [J]. Water Resources Research，1990(26):1465-1473．

[8]黎坤，林凱榮，江濤，等．?dāng)?shù)字濾波法在點源和非點源污染負(fù)荷分割中的應(yīng)用[J]．環(huán)境科學(xué)研究，2010(3)：298-303．

[9]ARNOLD J G, ALLEN P M．Automated methods for estimating baseflow and ground water recharge from stream flow [J]．Journal of the American Water Resources Association, 1999,35(2):411-424．

[10]ARNOLD J G, ALLEN P M, MUTTIAH R S, et al. Automated baseflow separation and recession analysis techniques [J]. Ground Water, 1995,33(6):1010-1018．

Segmentation of the Point Source and Non-point Source Pollution Load of Several Pollutants in the River

DAI Junfeng1,2YANG Yi1WANG Luyu1WANG Yuan1

(1.GuangxiKeyLaboratoryofEnvironmentalPollutionControlTheoryandTechnology,GuilinUniversityofTechnologyGuilin,Guangxi541004)

Hydrological estimation method and the digital filtering method are used to segment the point source pollution and non-point source pollution of sulfate, chloride, nitrate nitrogen and nitrite nitrogen in Nanliu River from 2003-2011.The results show that the non-point source pollution has become a relatively serious problem in Nanliu River. The calculations of hydrological estimation method show that the ratio of non-point source pollution to total pollution load of these four indexes is 0.82, 0.81, 0.82 and 0.72 respectively. And calculations of the digital filtering method indicate that ratio of four indexes is 0.72, 0.54, 0.68 and 0.63 respectively. The results show that the segmentation change trends of the pollution with two methods are close.

non-point source pollutionpoint source pollutiondigital filtering methodhydrological estimation method

國家自然科學(xué)基金(51569007)，廣西自然科學(xué)基金(2015GXNSFCA139004)，廣西高等學(xué)校高水平創(chuàng)新團隊及卓越學(xué)者計劃(002401013001)，廣西礦冶與環(huán)境科學(xué)實驗中心項目(KH2012ZD004)。

代俊峰，男，1980年生，博士，教授，主要從事水資源高效利用與水環(huán)境研究。

2015-09-16)