高緯區阿什河面源氮和磷污染輸出特征

馬廣文，王業耀，香 寶，蘇布達，胡 鈺，于 洋，王 光

1.中國環境監測總站，國家環境保護環境監測質量控制重點實驗室，北京 1000122.中國環境科學研究院，環境基準與風險評估國家重點實驗室，北京 100012

高緯區阿什河面源氮和磷污染輸出特征

馬廣文1，王業耀1，香 寶2，蘇布達2，胡 鈺2，于 洋1，王 光1

1.中國環境監測總站，國家環境保護環境監測質量控制重點實驗室，北京 1000122.中國環境科學研究院，環境基準與風險評估國家重點實驗室，北京 100012

監測數據；高緯度地區；氮和磷污染；基流分割；輸出量

河流富營養化不僅導致河流本身的經濟、生態、環境功能下降，也是導致其受納水體污染的重要原因[1-2]，而河流作為中國重要的淡水來源之一，也是湖泊、水庫、地下水等淡水庫的補給源，對中國工農業及社會的發展起著重要作用[3]。富營養化是國內外河流、湖泊、水庫等淡水系統面臨的一個嚴重環境問題[4-6]。國內外對水體富營養化氮和磷負荷進行了大量的研究[7-11]，據估算，中國水體氮和磷污染當中面源氮和磷污染分別占總量的81%和93%[12]。近年來利用水文模型對面源產污、匯污過程進行模擬成為國內外面源研究的熱點[13-16]，然而這些方法需要大量的氣象、水文、研究區域下墊面以及農業活動等人類對自然地貌進行改造行為的數據，以及區域內點源污染排放數據等。由于面源污染監測存在困難[17]，且中國缺乏水域面源污染常規連續監測數據，如何對模型的模擬結果進行驗證成為技術難點，這些因素都限制了模型的應用范圍。因此，利用數理統計模型對面源污染進行估算，目前仍然是面源污染研究不可或缺的手段之一。在理想的情況下，可獲取連續的水質水量監測數據，可以評估流域污染的季節性變化，然而在中國以往的水環境常規監測中，連續的水量水質同步監測日序列數據較為缺乏，如何利用有限的監測數據進行分析，已成為需要解決的科學問題之一[18]，國內利用水量和水質數據計算點源與面源污染負荷已有一些研究[19-21]，這些研究的研究區處在緯度較低區域，水量水質同步監測數據比較連續，缺少對于冰封期長、監測數據存在缺失條件下的點源與面源污染負荷輸出的計算。筆者探索處于高緯度地區阿什河氮和磷污染物輸出量計算，通過對水質和水量監測分析處理，利用改進埃特金(Aitkin)、相關性擬合插值法和基流分割法等，研究阿什河氮磷污染尤其是面源污染對松花江的影響特征，為環境管理部門對高緯度地區河流面源污染的治理提供參考。

1 實驗方法

1.1 研究區概況

阿什河位于黑龍江省哈爾濱市境內，全長257 km。阿什河流域的地理坐標為45°05′～45°49′N， 126°40′～127°42′E，流域面積為3 545 km2[22]，流域氣候冬季寒冷而漫長、夏季涼爽而短促，屬中溫帶大陸性季風氣候，年均降水量為545.7 mm左右，年降水主要集中在6—9月，其中，降雨集中在每年7—8月，降雪集中在每年11月—次年1月。 每年11月中旬—次年4月上旬為冰封期。阿什河中下游干流豐水期水面寬238～360 m，水深為4.00～4.70 m；枯水期水面寬10.00～23.5 m，水深為0.20～0.25 m。多年平均入境水量為5.16億m3，多年平均徑流量為4.58億m3，最大年徑流量為9.01億m3，最小年徑流量為1.10億m3。阿什河主要有12條支流水系：頭道河、二道河、大泥黑河、黃玉河、西泉眼水庫、洼河和海溝等，見圖1。阿什河水量季節性明顯，7—8月為豐水期;4—6月、9—10月為平水期;11—12月和次年的1—3月為枯水期[23]。阿什河流域上游森林植被覆蓋良好，中下游耕地廣布，流域內從上游到下游依次包括哈爾濱市的尚志市、五常市、阿城區、香坊區和道外區。

圖1 阿什河流域位置和水系分布Fig.1 Location and water system of Ashi River Basin

1.2 監測數據及分析方法

埃特金插值法主要原理是將一個高次多項式逼近缺值點的過程轉換成多個線性插值的線性組合：一個復雜的高次插值計算由多個線性計算逐步插值來實現。如果實測數據序列(xj,yj)，j=0, 1, 2, …,n, 選取其中m個(m

(k=2,3,…,m;i=k,k+1,…,m)

(1)

在指定的x=R處，運用線性組合求得近似函數值，若精度不夠，可再增加選點數(m值)。在河流的封凍期和解凍期，缺失若干個數據，因此，需要對原有埃特金插值法適當改進，原指定R點作為流動點(Rj)，取出m個(可取m=n/2)實測數據，成為補差數據列，該列數在實測數據中的序號，由數據的可調指針LB(0≤LB≤m)來確定，需要使Rj位于XLB-P(P=1, 2, …,m) 之間，然后用線性組合，求得一個等距序列(其中包括實測點)。因此，式(2)為

(j=1,2,…,n;k=2,3,…,m;

i=k,k+1,…,m;p=1,2,…,m)

(2)

式中：j為補差點序列號(含缺值點與實測點)；變量i為實測數據列中補差時線性組合的次數；變量k為取出實測數據列中的數據序列號，p為補差函數值的序號[24]。阿什河口2008—2010年TN、NH+4-N和TP濃度插值后結果見圖2。

圖2 2008—2010年阿什河口TN、NH+4-N和TP月均濃度結果Fig.2 TN, NH+4-N and TP average monthly concentration of Ashi River estuary from 2008 to 2010

改進后埃特金插值法計算結果精度較高，是由于每次插值組合都是線性的，所以計算簡單、數據可靠，高次插值不會振蕩，增加插值組合的點數可以進一步提高插值函數的精度。

流量數據為阿什河水文站流量監測數據，由于流量的監測數據上下游之間的相關性比較好，對于空間點位缺失的流量數據，利用相關性進行數據插值，依據回歸方程對缺失數據擬合插值。一般擬合曲線需要6個以上監測數據，呈冪函數相關。 利用Microsoft Excel軟件圖表向導和平滑線功能繪制過程線擬合。依照擬合公式再計算缺失的數據，得出缺失點的數據。阿什河口2008—2010年月均流量擬合結果見圖3。

圖3 2008—2010年阿什河口月均流量擬合結果Fig.3 Average monthly flow fitting results of Ashi River estuary from 2008 to 2010

1.3 面源污染輸出量估算方法

阿什河面源污染輸出量主要是在水質水量的監測并補差缺值點數據，運用基流分割法進行計算。河流不同時間的水質和水量監測數據，經過水文分割處理，用于估算進入水體的面源污染負荷。

1.3.1 基流分割法

根據水文學原理可知總徑流分為基流和地表徑流2部分，基流一般指來源于地下水或其他延遲部分的徑流，或者定義為下滲水到達地下水面并注入河道的部分[25]。 一般認為基流是河道內常年出現的那部分水流，主要受流域的氣候、植被、地質和土壤等的影響。國外的學者一般將總徑流分為地面徑流和基流2個部分[26]，該基流是指除地面徑流以外的徑流量。在中國一般認為總徑流由深層地下徑流、淺層地下徑流和地面徑流3個部分構成，其中深層地下徑流成分比較穩定，也稱基流[27]。目前國際上普遍采用的是這種直接將總徑流分為基流和地面徑流2部分的方法[28]，本文也采用這一概念的方法。該方法考慮了點源、面源污染形成規律，基流輸送的污染物代表了點源負荷和自然背景值，而地表徑流輸送的污染物為面源污染負荷。

基流分割方法中直線分割法具有直觀、易于操作等特點，應用較多。直線分割法分為平割法和斜割法2種。平割法(枯水期直線分割法)是指用直線連接流量過程線中不同拐點進行基流分割，在對年流量過程進行基流分割時，一般以枯水期月均流量的最小值作為基準，進行基流的分割。平割法一般適用于山區閉合流域的基流分割，對于洪水歷時短、地下徑流小的流域比較適用，適合阿什河流域的水文特征。

1.3.2 輸出量核算方法

根據水文分割法原理，基流分割是在河流基流和地表徑流劃分基礎上，將地表徑流狀態下水體中的污染物視為面源豐(平)水期天然背景值，將基流狀態下水體中的污染負荷視為點源及枯水期天然背景值，河流年輸出污染總負荷可表示為

(3)

式中：t為時間；Cp(t)為t時刻點源污染物質量濃度；Cn(t)為t時刻面源污染物質量濃度；Qp(t)為河流基流量；Qn(t)為地表徑流量；Wt為輸出污染總負荷。Wt可由監測斷面的水質、水量數據直接求出。點源污染負荷通過枯水期實測污染物質量濃度來計算，河流基流由河流的徑流分割來劃分。通過對式(3)進行離散化并移項處理，面源污染負荷由總污染負荷與點源污染負荷之差來計算[29-30]：

(4)

式中：Wn為面源污染負荷；Cni為第i次監測的面源污染物質量濃度；Qni為第i次監測的流量；Δt為第i次監測的時間段；Ci為第i次監測的污染物質量濃度；Qi為第i次監測的流量；Cpi為第i次監測的點源污染物質量濃度；Qpi為第i次監測點源時的流量，即河流基流量。

2 結果與討論

2.1 污染物年際間輸出量

計算阿什河2008—2010年氮磷污染輸出量，見表1。

表1 2008—2010年阿什河輸出污染負荷計算結果Table 1 Output pollution loading of Ashi River from 2008 to 2010

2.2 污染物年內輸出量分布

2008年各月TN、NH+4-N和TP負荷隨時間變化過程見圖4。

圖4 2008年TN、NH+4-N、TP污染負荷隨時間變化過程Fig.4 Change of TN, NH4+-N and TP pollutional loading of Ashi River in 2008

從圖4可以看出，5月TP輸出量最大，達129.2 t，面源TP輸出量為91.5 t，TP輸出量主要由面源污染造成，面源TP輸出量占總量的70.8%，其次，6月TP輸出量較大，再次為8月；5月TN輸出量最大，達886.9 t，面源TN輸出量為568.3 t，TN負荷主要由面源污染造成，面源TN輸出量占總量的64.1%，其次，6月TN輸出量較大，再次為4月；5月NH+4-N輸出量最大，達517.7 t，面源NH+4-N輸出量為420.5 t，NH+4-N輸出量主要由面源污染造成，面源NH+4-N輸出量占總量的81.2%，其次，6月NH+4-N輸出量較大，再次為8月。

2009年各月TP、NH+4-N、TN輸出量隨時間變化過程見圖5。

圖5 2009年TN、NH+4-N和TP隨時間變化過程Fig.5 Change of TN, NH+4-N and TP loading of Ashi River in 2009

從圖5可以看出，10月TP輸出量最大，達到38.8 t，面源TP輸出量為22.2 t，TP輸出量主要由面源污染造成，面源TP輸出量占總量的57.1%，其次，3月TP輸出量較大；9月TN輸出量最大，達450.6 t，面源TN輸出量為177.7 t，TN輸出量主要由點源污染造成，面源TN輸出量占總量的39.4%，其次，6月TN輸出量較大；3月NH+4-N輸出量最大，達到335.7 t，面源NH+4-N負荷為149.2 t，面源NH+4-N負荷占總量的44.5%，其次，2月NH+4-N輸出量較大。

2010年各月TP、NH+4-N、TN負荷隨時間變化過程見圖6。

圖6 2010年TN、NH+4-N和TP負荷隨時間變化過程Fig.6 Change of TN, NH+4-N and TP loading of Ashi River in 2010

3 結論

[1] THIEU V, BILLEN G, GARNIER J. Nutrient transfer in three contrasting NW European watersheds: the Seine,Somme,and Scheldt Rivers. A comparative application of the Seneque/ Riverstrahler modle[J].Water Research,2009,43(6):1 740-1 754.

[2] WITHERS P J A,SHARPLEY A N.Characterization and apportionment of nutrient and sediment sources in catchments[J]. Journal of Hydrology,2008,350(3-4):127-130.

[3] 李躍飛,夏永秋,李曉波,等. 秦淮河典型河段總氮總磷時空變異特征[J].環境科學,2013,34(1):91-97.

LI Y F, XIA Y Q, LI X B, et al. Temporal and spatial variations of total nitrogen and total phosphorus in the typical reaches of Qinhuai River[J]. Environmental Science,2013,34(1):91-97.

[4] 陳軍，權文婷，孫記紅. 太湖氮磷濃度與水質因子的關系[J].中國環境監測,2011,27(3):79-83.

CHEN J，QUAN W T，SUN J H．Relationship between Nitrogen，Phosphorus concentration and water quality factors[J]. Environmental Monitoring in China, 2011,27(3):79-83.

[5] SMITH V H,JOYE S B,HOWARTH R W. Eutrophication of freshwater and marine ecosystems[J].Limnol Oceanogr,2006(51):351-355.

[6] 蔡履冰. 太湖流域水體富營養化成因及防治對策的初步研究[J].中國環境監測,2003,19(3):52-54.

CAI L B. Preliminary study on the cause and the prevent ing and cur ing of the water body eutrophication of the Taihu Basin[J]. Environmental Monitoring in China,2003,19(3):52-54.

[7] STEINMAN A, CHU X F, OGDAHL M. Spatial and temporal variability of internal and external phosphorus loads in Mona Lake,Michigan[J]. Aqual Ecol,2009(43):1-18.

[8] 朱冠霖，劉德富，謝濤，等. 香溪河庫灣氮的時空分布及影響因子分析[J].中國環境監測,2013,29(6):73-78.

ZHU G L，LIU D F，XIE T, et al. The temporal and spatial distribution of Nitrogen in Xiangxi Bay and the analyze of the impact factor [J]. Environmental Monitoring in China,2013,29(6):73-78.

[9] CHEN Y X, LIU R M, SUN C C, et al. Spatial and temporal variations in nitrogen and phosphorous nutrients in the Yangtze River Estuary[J]. Marine Pollution Bulletin,2012(64):2 083-2 089.

[10] HOU D K, HE J, LYU C W,et al. Spatial variations and distributions of phosphorus and nitrogen in bottom sediments from a typical north-temperate lake, China[J]. Environmental Earth Sciences, 2014,71(7):3 063-3 079.

[11] 曾凡海，張晟，熊強，等. 三峽水庫干流氮和磷含量的季節變化[J].中國環境監測,2012,28(5):29-32.

ZENG F H，ZHANG S，XIONG Q，et al．Seasonal variation of Nitrogen and Phosphorus in Three-Gorge Reservoir [J]. Environmental Monitoring in China, 2012,28(5):29-32.

[12] EDWIN D O, ZHANG X L,YU T. Current status of agricultural and rural non-point source pollution assessment in China [J]. Environmental Pollution,2010(158):1 159-1 168.

[13] STEVEN W,THOMAS S, TREVOR F, et al. A modeling approach to water quality management of an agriculturally dominated watershed, kansas, USA [J]. Water Air Siol Pollut,2009(203):193-206.

[14] SERDAR G, ERDEM A. Modeling climate change effects on streams and reservoirs with HSPF [J]. Water Resour Manage,2010(24):707-726.

[15] QIU Z, CHEN Y X, GHULAM J L, et al. Model AVSWAT apropos of simulating non-point source pollution in Taihu lake basin [J]. Journal of Hazardous Materials,2010(174):824-830.

[16] SHEN Z Y, GONG Y W, LI Y H, et al. A comparison of WEPP and SWAT for modeling soil erosion of the Zhangjiachong Watershed in the Three Gorges Reservoir Area [J]. Agricultural Water Management,2009(96):1 435-1 442.

[17] MA G W, WANG S R. Temporal and spatial distribution changing characteristics of exogenous pollution load into Dianchi Lake, Southwest of China[J]. Environmental Earth Sciences,2015,74(5):3 781-3 793.

[18] 劉潔,陳曉宏,周純,等. 非點源污染在東江河流水環境中的貢獻比例估算[J].中國人口·資源與環境,2014,24(11):79-82.

LIU J, CHEN X H, ZHOU C, et al. Study on estimation method of the pollutant loads ratio from non-point source in the Dongjiang River[J]. China Population, Resources and Environment,2014,24(11):79-82.

[19] 張革,劉德富,宋林旭,等. 不同基流分割方法在香溪河流域的應用對比研究[J].長江流域資源與環境,2013,22(2):164-171.

ZHANG G, LIU D F, SONG L X, et al. Different baseflow separation methods in the Xiangxi River watershed [J]. China Population, Resources and Environment, 2013,22(2):164-171.

[20] 喬飛,孟偉,鄭丙輝,等．長江流域污染負荷核算及來源分析[J].環境科學研究,2013,26(1):80-87.

QIAO F, MENG W, ZHENG B H, et al. Pollutant accounting and source analyses in the Yangtze River Basin[J]. Research of Environmental Sciences,2013, 26(1):80-87.

[21] 荊紅衛，華蕾，郭婧，等. 北京市水環境非點源污染監測與負荷估算研究[J].中國環境監測,2012,28(6):106-111.

JING H W，HUA L，GUO J，et al. Non-point source pollution monitoring and its loads estimation study on surface water environment in Beijing[J]. Environmental Monitoring in China,2012,28(6):106-111.

[22] MA G W, WANG Y Y, BAO X, et al. Nitrogen pollution characteristics and source analysis using the stable isotope tracing method in Ashi River, northeast China[J]. Environmental Earth Sciences,2015,73(8):4 831-4 839.

[23] 楊育紅. 吉林省地表水非點源氨氮污染負荷研究[D].長春:吉林大學,2007.

YANG Y H. Study on surface water nonpoint sources NH3-N pollution loading in Jilin province [D].Changchun: Jilin University,2007.

[24] 任德記, 趙大友, 朱忠榮.改進埃特金插值法在監測數據補差中的應用[J].人民長江, 2001, 32(1): 44-46.

REN D J,ZHAO D Y,ZHU Z R.Application of improved Aitkin method in the makeup of the monitoring data[J].Yangtze River,2001,32(1): 44-46.

[25] 徐磊磊,劉敬林,金昌杰,等. 水文過程的基流分割方法研究進展[J].應用生態學報, 2011,22(11):3 073-3 080.

XU L L, LIU J L, JIN C J, et al. Baseflow separation methods in hydrological process research: A review[J]. Chinese Journal of Applied Ecology, 2011,22(11):3 073-3 080.

[26] SMAKHTIN V U. Low flow hydrology: A review[J]. Journal of Hydrology,2001,240(3):147-186.

[27] 董曉華,鄧霞,薄會娟,等．平滑最小值法與數字濾波法在流域徑流分割中的應用比較[J].三峽大學學報:自然科學版,2010,32(2):1-4.

DONG X H, DENG X, BO H J,et al.A comparison between smoothed minimma and digital filtering methods applied to catchment baseflow separation[J]. China Three Gorges Univ. (Natural Sciences),2010,32(2):1-4.

[28] 黨素珍，王中根，劉昌明. 黑河上游地區基流分割及其變化特征分析[J].資源科學,2011,33(12):1-8.

DANG S Z, WANG Z G, LIU C M. Baseflow separation and its characteristics analysis in the Upstream of Heihe River Basin[J]. Resources Science,2011,33(12):1-8.

[29] 鄭丙輝, 王麗婧, 龔斌. 三峽水庫上游河流入庫面源污染負荷研究[J].環境科學研究, 2009, 22(2): 125-131.

ZHENG B H, WANG L J, GONG B. Load of non-point source pollutants from Upstream Rivers into Three Gorges Reservoir[J]. Research of Environmental Sciences,2009,22(2):125-131.

[30] 涂安國,李英,莫明浩,等. 基于水文分割法的鄱陽湖入湖面源污染研究[J].人民長江, 2012, 43(1):63-66.

TU A G, LI Y, MO M H, et al. Study on non-point source pollution from rivers into Poyang Lake based on hydrological base flow segmentation method[J]. Yangtze River,2012,43(1):63-66.

Non-point Source Pollution Output Characteristics of Nitrogen and Phosphorus in the Ashi River at High Latitudes

MA Guangwen1,WANG Yeyao1,XIANG Bao2,SUBUDA2,HU Yu2,YU Yang1,WANG Guang1

1.State Environment Protection key Laboratory of Environmental Monitoring Quality Control, China National Environmental Monitoring Centre,Beijing 100012,China2.State Key Laboratory of Environmental Criteria and Risk Assessment, Chinese Research Academy of Environmental Sciences,Beijing 100012,China

The improved Aitkin and correlation fitting interpolation method were used to complemented missing data by limited monitoring conditions, the base flow segmentation method were then used to estimated non-point source of nitrogen and phosphorus pollution output load of the Ashi River into the Songhua River, its pollution characteristics were analyzed, based on monitoring data of water quality and flow in the Ashi River from 2008 to 2010. The results showed that total nitrogen (TN) about 5 436.6 t, ammonia (NH+4-N) about 3 057.8 t, total phosphorus (TP) about 554.8 t of the Ashi River were discharged into the Songhua River per year, which point source pollution were 2 775.8, 1440.7 and 250.8 t, non-point source pollution were 2 660.7, 1 630.7, 304.1 t respectively, non-point source pollution accounts for about 50% of the total output. Ashi River pollution output amount emerged TP reduce, TN and NH+4-N added feature from 2008 to 2010. The TP, TN and NH+4-N amount of the Ashi River into the Songhua River were comprehensively analyzed each month from 2008 to 2010, the pollution output peaked in March-May or July-September, mainly caused by the spring flood or summer flood, when the peak pollution source was non-point source pollution generally.

monitoring data;high latitudes;nitrogen and phosphorus pollution;base flow separation;output

2016-04-20；

2016-05-25

國家自然科學基金(31400449)；國家水體污染控制與治理科技重大專項(2014ZX07201-009-01)

馬廣文(1980-),男,內蒙古呼和浩特人，博士,高級工程師。

王業耀

X824

A

1002-6002(2017)02- 0047- 08

10.19316/j.issn.1002-6002.2017.02.08