李中原, 王國重, 左其亭, 屈建鋼

(1.河南省水文水資源局, 河南 鄭州 450003; 2.黃河水文水資源科學(xué)研究院, 河南 鄭州 450004;3.鄭州大學(xué) 水利與環(huán)境學(xué)院, 河南 鄭州 450001; 4.河南省水土保持監(jiān)督監(jiān)測總站, 河南 鄭州 450008)

應(yīng)用分形理論估算丹江口水庫水源區(qū)總氮、總磷的流失量

李中原1, 王國重2, 左其亭3, 屈建鋼4

(1.河南省水文水資源局, 河南 鄭州 450003; 2.黃河水文水資源科學(xué)研究院, 河南 鄭州 450004;3.鄭州大學(xué) 水利與環(huán)境學(xué)院, 河南 鄭州 450001; 4.河南省水土保持監(jiān)督監(jiān)測總站, 河南 鄭州 450008)

[目的] 掌握區(qū)域農(nóng)業(yè)面源污染物的流失狀況，證明分形計(jì)算結(jié)果的可靠性。[方法] 采用分形理論對丹江口水庫水源區(qū)河南省所在的區(qū)域2013年總氮(TN)，總磷(TP)的流失量進(jìn)行估算，并與輸出系數(shù)法的值進(jìn)行比較。[結(jié)果] 研究區(qū)污染物的流失以氮為主，TN的流失量是TP的7.16倍；禽畜養(yǎng)殖引起的流失量最多，占了總量的69.93%，農(nóng)田化肥中的流失量占21.99%；灌河產(chǎn)生的污染物最多，其次是丹江、淇河、滔河，但在流失強(qiáng)度上卻是滔河最大，其次是灌河、丹江、淇河；通過配對t檢驗(yàn)分析將分形方法計(jì)算的污染物總量與輸出系數(shù)法相比，在0.1水平下有統(tǒng)計(jì)學(xué)意義，說明前者較后者更能反映實(shí)際。[結(jié)論] 用分形方法估算區(qū)域污染物的面源污染是可行的，應(yīng)當(dāng)加強(qiáng)區(qū)域禽畜糞尿和農(nóng)田化肥的管理，以確保丹江口水庫水質(zhì)。

分形理論; 水系分維; 農(nóng)業(yè)面源污染; TN,TP流失

文獻(xiàn)參數(shù)： 李中原, 王國重, 左其亭, 等.應(yīng)用分形理論估算丹江口水庫水源區(qū)總氮、總磷的流失量[J].水土保持通報(bào)，2017,37(3)：302-306.DOI:10.13961/j.cnki.stbctb.20170407.001； Li Zhongyuan, Wang Guozhong, Zuo Qiting, et al. Applying fractal theory to estimate loss of total nitrogen and total phosphorus in Danjiangkou reservoir catchment[J]. Bulletin of Soil and Water Conservation, 2017,37(3)：302-306.DOI:10.13961/j.cnki.stbctb.20170407.001

大大小小、蜿蜒曲折的河流交匯形成水系，如何反映和分析這種蜿蜒特性及其網(wǎng)狀結(jié)構(gòu)的問題一直困擾著學(xué)術(shù)界。Mandelbrot 1977年率先將分形的概念引入了地理水文學(xué)，開辟了流域地貌形態(tài)特征定量化研究的先河[1]。1982年其專著“自然界的分形幾何”的出版，標(biāo)志著分形理論的初步形成[2]。分形理論給人類研究自然界紛繁復(fù)雜的客觀規(guī)律及其內(nèi)在聯(lián)系提供了新的思路和方法，其優(yōu)勢在于能夠揭示看起來支離破碎、處處不可微、不光滑，貌似混亂、無規(guī)則、隨機(jī)現(xiàn)象的內(nèi)部規(guī)律，能夠準(zhǔn)確、完整地界定實(shí)物的主流趨勢性質(zhì)，這是傳統(tǒng)方法難以達(dá)到的，因而也成為人類解決復(fù)雜問題的有力工具。

農(nóng)業(yè)面源污染具有隨機(jī)性、分散性、隱蔽性、廣泛性和不確定性的特點(diǎn)，難以監(jiān)測，一般需借助模型來定量分析，如SWAT，AGNPS，SWMM等[3]。這些模型含有大量反映土壤、水文、氣象、土地利用等方面的參數(shù)，且參數(shù)的率定過程較為繁瑣，其模擬的精度就會大打折扣[4-5]。輸出系數(shù)法所需參數(shù)少、操作便捷、有一定的精度，而且不計(jì)污染產(chǎn)生與發(fā)展的過程，因此頗受研究者的青睞[6]。

農(nóng)業(yè)面源污染帶有非線性的特征[7]，就可以采用分形理論來研究區(qū)域的農(nóng)業(yè)面源污染，但是這方面的報(bào)道目前還比較少見。河流水系分形是目前分形研究的熱點(diǎn)，大量研究表明：水系分維數(shù)與地貌的發(fā)育、地貌的類型、土地的利用等因素密切相關(guān)，水系分維數(shù)不僅反映了水系的發(fā)育程度，還代表它所處流域的地貌侵蝕發(fā)育的階段[8-9]。但許多研究只是通過GIS軟件計(jì)算水系的分維數(shù)，而沒有作進(jìn)一步的研究[10]，如將水系分維數(shù)與環(huán)境污染、流域生態(tài)相聯(lián)系。

農(nóng)業(yè)面源污染是影響丹江口水庫水質(zhì)的首要污染源[11]，總氮(TN)，總磷(TP)又是導(dǎo)致農(nóng)業(yè)面源污染的主要因素[12]。本研究擬以丹江口水庫水源區(qū)河南省所在的區(qū)域?yàn)檠芯繉ο?，根?jù)ArcGIS軟件計(jì)算的水系分維值，結(jié)合對該區(qū)域農(nóng)業(yè)面源污染情況的調(diào)查，運(yùn)用分形理論來估算該區(qū)域TN，TP的污染情況，并與輸出系數(shù)法的結(jié)果相比較，以期說明分形計(jì)算結(jié)果的可靠性。

1 研究區(qū)概況

丹江口水庫由位于湖北省境內(nèi)的漢庫和位于河南省境內(nèi)的丹庫組成，丹庫正是所要研究的區(qū)域，地處河南省西南部，包含南陽市的西峽縣、淅川縣，以及內(nèi)鄉(xiāng)縣、鄧州市、三門峽盧氏縣、洛陽欒川縣的部分鄉(xiāng)鎮(zhèn)，流域總面積為8 047 km2,位于東經(jīng)110°52′—112°，北緯32°54′—34°，海拔121～2 212.5 m。

區(qū)域水系主要為丹江及其支流灌河、淇河和滔河。丹江全長378.6 km，源于陜西鳳凰山，經(jīng)商南縣月亮灣入淅川，在淅川境內(nèi)長116.6 km，自滔河鎮(zhèn)申明鋪以下為丹江口庫區(qū)；老灌河發(fā)源于欒川縣小廟嶺，經(jīng)盧氏和西峽縣，在淅川縣注入丹江，上集鎮(zhèn)奕子營以下為庫區(qū)；淇河發(fā)源于盧氏縣花園寺西，至淅川縣寺灣鄉(xiāng)老君洞附近匯入丹江，全長150 km；滔河發(fā)源于陜西商南縣白魯?shù)A鄉(xiāng)白龍洞，流經(jīng)商南縣、鄖縣，在淅川縣滔河鄉(xiāng)滔河街注入丹江，河流全長155.2 km。根據(jù)河南省2014年的統(tǒng)計(jì)年鑒，結(jié)合走訪調(diào)查，確定了各水系的人口、耕地面積、畜禽養(yǎng)殖情況(表1)。

表1 丹江口水庫2013年基本情況

2 計(jì)盒維數(shù)的計(jì)算與輸出系數(shù)法簡介

2.1 計(jì)盒維數(shù)的計(jì)算

分形維數(shù)是分形理論的主要參數(shù)，有多種定義和計(jì)算方法，如：相似維數(shù)(Ds)、豪斯道夫維數(shù)(DH)、計(jì)盒維數(shù)(Db)、信息維數(shù)(Di)、關(guān)聯(lián)維數(shù)(Dg)、容量維數(shù)(DC)等。其中，計(jì)盒維數(shù)(Db)因計(jì)算簡便而最為常用，反映的是分形體對空間的占據(jù)程度，用邊長為r的小方盒子去覆蓋分形體，非空小方盒的數(shù)量記為N(r)，不斷縮小r的取值，相應(yīng)地得到一系列的N(r)值，當(dāng)r→0時，則計(jì)盒維數(shù)定義為：

(1)

實(shí)際的工作中，一般是借助ArcGIS軟件利用網(wǎng)格法來計(jì)算計(jì)盒維數(shù)[13]。GIS技術(shù)可以獲得研究區(qū)的總體數(shù)據(jù)，比傳統(tǒng)的人工計(jì)算或計(jì)算機(jī)編程計(jì)算更為客觀、準(zhǔn)確、經(jīng)濟(jì)、高效。具體步驟如下：

(1) 先由ArcGIS 工具箱中的conversion tools工具設(shè)置需要輸出的柵格尺寸，將矢量河網(wǎng)數(shù)據(jù)轉(zhuǎn)換成柵格文件；

(2) 打開相應(yīng)柵格文件的屬性表，統(tǒng)計(jì)柵格數(shù)量；

(3) 由Excel軟件將相應(yīng)數(shù)據(jù)點(diǎn)繪在雙對數(shù)坐標(biāo)圖上，相應(yīng)直線的斜率即為水系分維數(shù)。

由中國科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心數(shù)據(jù)平臺獲得該區(qū)域的DEM數(shù)據(jù)，采用ArcGIS軟件和網(wǎng)格法提取并計(jì)算各水系的計(jì)盒維數(shù)。水系分維數(shù)乘以不同污染物的平均流失率，就得到了考慮下墊面作用的面源污染物的流失率[14]。

2.2 輸出系數(shù)法簡介

輸出系數(shù)法就是利用污染物輸出系數(shù)估算流域輸出的面源污染負(fù)荷，20世紀(jì)70年代初北美國家在評價土地利用和湖泊富營養(yǎng)之間的關(guān)系時提出的，經(jīng)不斷完善與發(fā)展，該法在大的流域尺度上也有一定的精度，其表達(dá)式一般為：

(2)式中：Lj——流域內(nèi)污染物j的總負(fù)荷〔t/(hm2·a)〕；j——污染物類型；i——流域內(nèi)土地、禽畜或人口分類，共有n類；Cij——污染物j在第i類土地、牲畜、人口的輸出系數(shù)；Ai——第i類土地的面積(hm2)或牲畜(頭)、人口(人)的數(shù)目；P——降雨輸入的污染物數(shù)量〔t/(hm2·a)〕，這里不考慮此項(xiàng)的影響。

3 結(jié)果與分析

3.1 各水系面源污染物TN，TP的流失率

中國農(nóng)用化肥流失率普遍較高,其中氮肥流失率10%～25%,磷肥流失率4%～5%，這是因?yàn)榱自啬鼙煌寥缽?qiáng)烈吸附。根據(jù)對研究區(qū)西峽縣、淅川縣施肥及其流失情況的調(diào)查，結(jié)合統(tǒng)計(jì)年鑒以及前人的研究[15]，確定研究區(qū)2013年氮、磷的流失率分別取該年施肥量的12.5%和3%。全國范圍內(nèi)畜禽排泄物中糞便的流失率為2%～8%，液體排泄物的流失率則可達(dá)到50%[16]。

研究區(qū)養(yǎng)殖較多的是大牲畜、豬、羊和家禽，根據(jù)對養(yǎng)殖業(yè)污染調(diào)查結(jié)合相關(guān)研究結(jié)果[17]，確定了畜禽排泄物相應(yīng)的流失率(表2)。根據(jù)研究區(qū)農(nóng)民的生活習(xí)慣和生活方式，按5%計(jì)算村民糞尿流入水體的量，生活污水的流失率為50%[18]。

表2 丹江口水庫禽畜糞污的流失率

注:大牲畜糞、尿值為牛、馬、騾、驢的均值；家禽糞為雞、鴨、鵝糞的均值。下同。

各水系所在的地域和環(huán)境條件不同，污染物的流失率必然有差異。分形作用下農(nóng)田、禽畜養(yǎng)殖、農(nóng)村生活污染物的流失率如表3—4所示。

表3 分形對丹江口水庫農(nóng)田化肥、村民生活中TN，TP流失率的影響

表4 分形對丹江口水庫農(nóng)村禽畜養(yǎng)殖中TN，TP流失率的影響

水系分維反映了水系發(fā)育的復(fù)雜程度，分維數(shù)越大說明水系發(fā)育的越好。由表4可知，由于引入了分形維數(shù)，各水系的TN，TP流失率隨著分維值的不同而不同，淇河水系TN，TP的流失率稍大，因其分維值稍大，相應(yīng)的滔河的流失率較低。

3.2 面源污染物TN，TP的流失量

污染物的流失率乘以其產(chǎn)生量既為污染物的流失量。據(jù)實(shí)地調(diào)查，河南省2014年統(tǒng)計(jì)年鑒和相關(guān)文獻(xiàn)[19-20]，結(jié)合河南省農(nóng)村生活用水定額和生活污水排放標(biāo)準(zhǔn)，確定了研究區(qū)2013年各水系農(nóng)田化肥中TN，TP的施用量、禽畜養(yǎng)殖中TN，TP的排泄量、生活污染物中TN，TP的產(chǎn)生量見表5所示。相應(yīng)污染物的流失量見表6。

從表6可以看出，研究區(qū)2013年共流失污染物6 864.458 t，其中TN總量為6 022.855 t，TP總量為841.603 t，前者是后者的7.16倍，說明研究區(qū)產(chǎn)生的污染物以氮的流失為主；禽畜流失的污染物量為4 800.116 t，占研究區(qū)污染物總量的69.93%，這表明應(yīng)當(dāng)加強(qiáng)對禽畜糞尿的管理；農(nóng)田中化肥產(chǎn)生的污染物量為1 509.608 t，占研究區(qū)污染物總量的21.99%，說明農(nóng)田中的養(yǎng)分流失也是農(nóng)業(yè)面源污染的重要來源。

表5 丹江口水庫污染物全氮(TN)，全磷(TP)的產(chǎn)生量 t/a

表6 丹江口水庫污染物全氮(TN)，全磷(TP)的流失量

從水系角度，灌河產(chǎn)生的污染物最多，對研究區(qū)污染的貢獻(xiàn)也最大，達(dá)43.97%；其次是丹江；污染物流失最少的是滔河，對環(huán)境污染的貢獻(xiàn)率僅為3.53%。從流失強(qiáng)度上來看，滔河>灌河>丹江(0.004 05)>淇河(0.004 02)，而水系分維數(shù)上卻是：淇河>丹江>灌河>滔河，水系分維數(shù)越大其污染物的流失強(qiáng)度越小，流失強(qiáng)度小并不代表其產(chǎn)生的流失量就少，這取決于流域地形地貌、土壤、植被、人類活動等的綜合作用。

3.3 與輸出系數(shù)法計(jì)算結(jié)果的比較

本文在用分形方法計(jì)算的同時，也用輸出系數(shù)法估算了研究區(qū)的農(nóng)業(yè)面源污染物TN，TP的流失量，結(jié)果如表7所示。2種方法計(jì)算的污染物流失量適合用配對t檢驗(yàn)進(jìn)行分析，其總量差值的t檢驗(yàn)結(jié)果詳見表8。

表7 兩種方法計(jì)算的污染物流失量結(jié)果

表8 兩種方法計(jì)算的污染物流失量的配對t檢驗(yàn)

由表7—8不難看出，2種方法計(jì)算的結(jié)果不管是TN還是TP，這幾條河流的差值不大，雖然數(shù)據(jù)樣本量不大，兩者總量差值配對t檢驗(yàn)的p=0.096，在顯著性水平0.1下具有統(tǒng)計(jì)學(xué)意義，也能夠說明分形方法的計(jì)算結(jié)果相對于輸出系數(shù)法是可靠的。此外，各個水系在地形地貌、土壤、植被、巖性等方面存在差別，分形方法能夠反映這些因素的影響：灌河、淇河、丹江的分維數(shù)均大于1，相應(yīng)的TN，TP的流失量大于輸出系數(shù)的結(jié)果；滔河的分維數(shù)小于1，相應(yīng)的流失量也小于輸出系數(shù)的值。

4 結(jié) 論

(1) 研究區(qū)污染物的流失主要表現(xiàn)為氮的流失，TN的流失量是TP的7.16倍，其中禽畜養(yǎng)殖產(chǎn)生的污染物占了污染物總量的69.93%，農(nóng)田中化肥的流失占了總污染量的21.99%，這表明：應(yīng)當(dāng)加強(qiáng)對禽畜糞尿的管理，采取措施減少農(nóng)田中養(yǎng)分的流失。

(2) 灌河產(chǎn)生的污染物最多，其次是丹江、淇河、滔河，但在流失強(qiáng)度上卻是滔河最大，其次是灌河、丹江、淇河，這表明分形方法能夠反映水系對農(nóng)業(yè)面源污染的影響，因?yàn)榱饔虻匦蔚孛病⑼寥?、植被、人類活動等因素與水系分維息息相關(guān)。

(3) 通過配對t檢驗(yàn)分析，分形方法計(jì)算的污染物總量與輸出系數(shù)法的差值，在0.1水平下具有統(tǒng)計(jì)學(xué)意義，說明前者的結(jié)果較后者可靠；此外，分形法考慮了下墊面因素的影響，并以水系分維數(shù)來表征，而輸出系數(shù)法沒有考慮這些差異，將之平均化，認(rèn)為其流失率相同，使得分形法的結(jié)果要優(yōu)于后者，更能反映實(shí)際。

[1] 江滔,武偉,劉洪斌.基于DEM的重慶長江流域水系分維估算[J].中國農(nóng)學(xué)通報(bào),2010,26(15):365-368.

[2] Mandelbrot B B, Wheeler J A. The Fractal Geometry of Nature[M]. New York: Freeman W H, 1982.

[3] 孫本發(fā),馬友華,胡善寶,等.農(nóng)業(yè)面源污染模型及其應(yīng)用研究[J].農(nóng)業(yè)環(huán)境與發(fā)展,2013,30(3):1-5.

[4] 耿潤哲,王曉燕,焦帥,等.密云水庫流域非點(diǎn)源污染負(fù)荷估算及特征分析[J].環(huán)境科學(xué)學(xué)報(bào),2013,33(5):1484-1492.

[5] 楊立夢,付永勝,高紅濤.四川省茫溪河流域非點(diǎn)源污染負(fù)荷研究[J].重慶理工大學(xué)學(xué)報(bào):自然科學(xué)版,2014,28(11):57-63.

[6] 周躍龍,汪懷建,余輝,等.應(yīng)用輸出系數(shù)模型對太湖流域面源污染負(fù)荷測算研究[J].江西農(nóng)業(yè)大學(xué)學(xué)報(bào),2014,36(3):678-683.

[7] 肖新成,倪九派,何炳輝,等.三峽庫區(qū)重慶段農(nóng)業(yè)面源污染負(fù)荷的區(qū)域分異與預(yù)測[J].應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào),2014,22(4):634-646.

[8] 王平,劉少峰,鄭洪波,等.四川盆地東部弧形構(gòu)造控制的地形和水系發(fā)育[J].第四紀(jì)研究,2013,33(3):461-470.

[9] Turcotte D L. Self-organized complexity in geomorphology: Observations and models[J]. Geomorphology, 2007,91(3):302-310.

[10] 王碹,陳雯靜,徐璐.基于GIS的小流域地貌分形維數(shù)測定方法研究[J].沈陽農(nóng)業(yè)大學(xué)學(xué)報(bào),2011,42(4):500-503.

[11] 姜世英,韓鵬,賈振邦,等.南水北調(diào)中線丹江口庫區(qū)農(nóng)業(yè)面源污染PSR評價與基于GIS的空間特征分析[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2010,29(11):2153-2162.

[12] 陳英旭.總氮總磷應(yīng)納入重點(diǎn)流域水體污染物總量減排指標(biāo)[J].中國產(chǎn)業(yè),2011(4):6-6.

[13] 李精忠,劉劍煒,楊澤龍.DEM數(shù)據(jù)谷地分維值的估算[J].武漢大學(xué)學(xué)報(bào):信息科學(xué)版,2014,39(11):1277-1281.

[14] 姚素珍.基于分形理論的清河水庫流域非點(diǎn)源污染研究[D].遼寧 阜新：遼寧工程技術(shù)大學(xué),2012.

[15] 夏穎,汪榮勇,高立,等.湖北省油菜原棉花輪作系統(tǒng)地表徑流氮磷流失特征[J].湖北農(nóng)業(yè)科學(xué),2014,53(23):5751-5754.

[16] 李曉光,周其文,胡梅,等.中國畜禽糞便污染現(xiàn)狀及防治對策[J].中國農(nóng)學(xué)通報(bào),2008,24(S):77-80.

[17] 全國規(guī)?；笄蒺B(yǎng)殖污染調(diào)查辦公室.全國規(guī)?；笄蒺B(yǎng)殖業(yè)污染情況調(diào)查技術(shù)報(bào)告[R],2001.

[18] 孫瑞敏.湖北省農(nóng)村生活污水水量水質(zhì)調(diào)查與分析[D].武漢：武漢理工大學(xué),2010.

[19] 劉欽普.河南省化肥施用強(qiáng)度地理分布及其環(huán)境風(fēng)險(xiǎn)評價[J].河南農(nóng)業(yè)科學(xué),2014,43(6):66-70.

[20] 國務(wù)院第一次全國污染源普查辦公室.第一次全國污染源普查畜禽養(yǎng)殖業(yè)源產(chǎn)排污系數(shù)手冊[R],2009.

Applying Fractal Theory to Estimate Loss of Total Nitrogen and Total Phosphorus in Danjiangkou Reservoir Catchment

LI Zhongyuan1, WANG Guozhong2, ZUO Qiting3, QU Jiangang4

(1.HydrologyandWaterResourcesBureauinHenanProvince,Zhengzhou,He’nan450003,China;2.HydrologyandWaterResourcesofYellowRiverScientificResearchInstitute,Zhengzhou,He’nan450004,China,China; 3.CollegeofWaterConservancyandEnvironment,ZhengzhouUniversity,Zhengzhou,He’nan450001,China;4.SoilandWaterConservationSupervisionandInspectionStationinHenanProvince,Zhengzhou,He’nan450008,China)

[Objective] To grasp the situation of agricultural non-point source pollution in Danjiangkou reservoir catchment areas in He’nan Province, and to prove the reliability of the fractal theory calculation results. [Methods] The fractal theory was applied to estimate the loss of total nitrogen(TN) and total phosphorus(TP) in the Danjiangkou reservoir region of He’nan Province in 2013. The estimated values were compared with the ones from output coefficient method. [Results] Pollutants loss was given priority to with nitrogen in the area， the TN loss was 7.16 times of TP loss. The livestock farming contributed the largest amount of loss, accounting for 69.93% of the total loss. The second was the loss from farmland fertilizer, taking up 21.99% of the total one. The pollutant yield of Guanhe River was the most, followed by Danjiang River, Qihe River and Taohe River. The pollutant erosion intensity had a rank as Taohe River>Guanhe River>Danhe River>Qihe River. The significant paired T test(<0.1) showed that the pollutant total amounts calculated by fractal method was coincided with the ones by output coefficient method, which indicated the former could reflect the actual situation more than that of the latter. [Conclusions] Fractal method is feasible to estimate logical agricultural non-point source pollution. It is time to strengthen management of livestock manure and chemical fertilizer to ensure the water quality of Danjiangkou Reservoir.

fractal theory; fractal dimension of water system; agricultural non-point source pollution; total nitrogen and total phosphorus loss

2016-11-06

2016-11-25

河南省科技攻關(guān)計(jì)劃“基于分形理論的丹江口水庫水源區(qū)農(nóng)業(yè)面源污染研究及防治措施”(GG201412)

李中原(1964—),男(漢族),河南省鄭州市人,碩士,教授級高級工程師，從事水土保持方面的研究工作。E-mail:lizy0723@163.com。

B

1000-288X(2017)03-0302-05

X524