麻德明，石洪華，豐愛平

(國家海洋局第一海洋研究所，青島 266061)

水環(huán)境污染問題已經(jīng)是影響人類社會(huì)可持續(xù)發(fā)展的主要制約因素，正越來越受到人們的關(guān)注。全球有30%—50%的地表水體受到非點(diǎn)源污染的影響［1］。研究表明，非點(diǎn)源污染是導(dǎo)致地表水污染的主要原因，其中又以農(nóng)業(yè)非點(diǎn)源污染貢獻(xiàn)率最大。在美國，即使點(diǎn)源污染全面控制之后，水體的質(zhì)量也并未因此而有所改善。美國的非點(diǎn)源污染約占污染總量的2/3，其中農(nóng)業(yè)非點(diǎn)源污染占非點(diǎn)源污染總量的68%—83%，農(nóng)業(yè)非點(diǎn)源污染已經(jīng)成為全美河流污染的第一污染源［2］。荷蘭農(nóng)業(yè)非點(diǎn)源提供的總氮、總磷分別占水環(huán)境污染總量的60%和40%—50%［3］。在我國，非點(diǎn)源污染問題也日益嚴(yán)重，農(nóng)業(yè)非點(diǎn)源污染已經(jīng)成為水質(zhì)惡化的主要原因之一。因此，加強(qiáng)農(nóng)業(yè)非點(diǎn)源污染研究對我國水污染防治具有重要意義。

農(nóng)業(yè)非點(diǎn)源污染是指在農(nóng)業(yè)生產(chǎn)活動(dòng)中，農(nóng)田中的土粒、氮素、磷、農(nóng)藥及其他有機(jī)或無機(jī)污染物質(zhì)，在降水或灌溉過程中，通過農(nóng)田地表徑流、農(nóng)田排水和地下滲漏，進(jìn)入水體而形成的水環(huán)境污染［4］。隨著沿海經(jīng)濟(jì)的快速發(fā)展，農(nóng)業(yè)非點(diǎn)源污染正日益成為水體富營養(yǎng)化最主要的污染源，并已成為當(dāng)前限制我國經(jīng)濟(jì)可持續(xù)發(fā)展的重要因素，控制陸源污染是防治海洋污染的關(guān)鍵，而陸源污染物定量估算則是控制陸源污染的基礎(chǔ)。

海灣區(qū)域非點(diǎn)源污染的研究在國外較早引起學(xué)者的重視［5－7］，國內(nèi)在該領(lǐng)域的研究還相對薄弱。在我國沿海城市工業(yè)、生活污水等點(diǎn)源污染治理程度逐步提高的前提下，近岸海域水環(huán)境仍然沒有明顯的改善，局部區(qū)域甚至還面臨著惡化的趨勢［8］，非點(diǎn)源污染作為海灣主要污染源之一日益受到管理部門和研究學(xué)者的重視。由于非點(diǎn)源污染的發(fā)生具有隨機(jī)性和復(fù)雜性，計(jì)算非點(diǎn)源污染負(fù)荷的難度比較大，隨著水環(huán)境實(shí)施總量控制的管理，非點(diǎn)源污染的定量化尤為重要［9］。目前，估算非點(diǎn)源污染的傳統(tǒng)方法有很多［10－23］，歸納起來主要可以分為4種:斷面實(shí)測總負(fù)荷減去點(diǎn)源負(fù)荷的方法、水文估算法、輸出系數(shù)法(單位負(fù)荷法)和模型估算法，這4種方法各有優(yōu)缺點(diǎn)。一般來說，農(nóng)業(yè)非點(diǎn)源污染包括農(nóng)田徑流、畜禽養(yǎng)殖、水產(chǎn)養(yǎng)殖和農(nóng)村生活4個(gè)主要來源。本文以流域單元?jiǎng)澐值乃枷霝榛A(chǔ)，借助GIS技術(shù)，以匯水區(qū)為單元，嘗試建立適宜的機(jī)理模型，提出一種基于流域單元的農(nóng)業(yè)非點(diǎn)源污染負(fù)荷估算方法，應(yīng)用于萊州灣入海流域，從以上4個(gè)方面來定量地估算農(nóng)業(yè)非點(diǎn)源主要污染物入?？偭俊?/p>

1 研究區(qū)概況

萊州灣位于山東半島北側(cè)，渤海南部，是山東省最大的海灣，也是山東半島城市群經(jīng)濟(jì)圈的腹地。灣口東起屺姆島的高角，西至現(xiàn)代黃河新入?？?。萊州灣總面積23370 km2，海岸線577.91 km，海域面積約11603 km2，其中淺海面積8726 km2，灘涂廣袤，海洋漁業(yè)生物資源和土地資源豐富，盛產(chǎn)蟹、蛤、蝦等，還是許多海洋生物的產(chǎn)卵場和索餌場。本區(qū)域海洋漁業(yè)、海洋油氣業(yè)、海洋鹽業(yè)、海洋化工、海洋礦業(yè)、海洋交通運(yùn)輸業(yè)等主要海洋產(chǎn)業(yè)優(yōu)勢明顯，在全省海洋經(jīng)濟(jì)中占有重要地位。本海區(qū)地處中緯度，水溫分布有明顯的季節(jié)性，冬季(2月)近海在0℃以下，有冰凍出現(xiàn);春季(5月)在14℃左右。夏季(8月)表層水溫可達(dá)26℃;秋季(11月)在11℃左右。東部海水較清，西南部較混濁，海灣東部有兩個(gè)附屬小灣，即龍口灣和太平灣。灣東南部海底有一西北東南走向的水下淺灘，即“萊州淺灘”。萊州灣的海岸由3部份組成:高角至虎頭崖為沙質(zhì)岸段，虎頭崖至淄脈溝為粉砂淤泥質(zhì)岸，淄脈溝以北為現(xiàn)代黃河三角洲的迅速淤進(jìn)型岸段。沿岸自北向南、自西向東有小島河、溢洪河、廣利河、淄脈溝、小清河、彌河、白浪河、虞河、濰河、膠萊河、沙河、王河和界河等十幾條較大河流，這些河流是造成海岸的最活躍的參與者。隨著經(jīng)濟(jì)、人口的迅速增長，萊州灣的環(huán)境問題日趨嚴(yán)重［24］。

2 匯水區(qū)和子流域劃分方法

河網(wǎng)的提取、匯水區(qū)和子流域的劃分，借助已經(jīng)成熟的美國ESRI公司的ARC/INFO軟件中的ArcHydro模塊。

2.1 ArcHydro模型介紹

ArcHydro模塊是由得克薩斯大學(xué)的水資源研究中心(Center for research in Water Resources，CRWR)基于ARC/INFO聯(lián)合ESRI開發(fā)的，其數(shù)據(jù)模型是把水文的地表水流模擬和GIS數(shù)據(jù)結(jié)構(gòu)相結(jié)合形成的。ArcHydro數(shù)據(jù)模型包括5個(gè)部分:Drainage(匯流)、Hydrography(水文地形)、Channel(河道)、Network(水文網(wǎng)絡(luò))和TimeSeries(時(shí)間序列)。

2.2 流域劃分

采用地表徑流漫流模型，模擬地表徑流在地表從高到低的流動(dòng)來生成水系，即借助ArcGIS軟件中的ArcHydro水文分析模塊，先填洼，然后根據(jù)坡降確定每個(gè)格網(wǎng)的水流方向的方法提取匯流網(wǎng)絡(luò)，最終生成匯水區(qū)。利用ArcGIS可以提取不同級別的匯流網(wǎng)絡(luò)，對同一流域可以劃分為面積由大到小不同等級的子匯流區(qū)域。

(1)數(shù)據(jù)預(yù)處理

數(shù)據(jù)源采用SRTM DEM(分辨率90 m)，首先利用ArcGIS的Mosaic To New Raster工具，把涉及研究區(qū)域的兩幅DEM進(jìn)行合并，然后利用 ArcGIS的Extract by Polygon工具，裁減出研究區(qū)域的DEM數(shù)據(jù)。如圖1所示。

圖1 研究區(qū)域DEMFig.1 the study area of DEM

(2)洼地填充

DEM是比較光滑的地形表面的模擬，但由于內(nèi)插的原因以及一些真實(shí)地形(如喀斯特地貌)的存在，使得DEM表面存在著一些凹陷的區(qū)域。在進(jìn)行地表水流模擬時(shí)，會(huì)得到不合理的或錯(cuò)誤的水流方向。因此，在進(jìn)行水流方向的計(jì)算之前，利用ArcHydro模塊中Fill功能，把洼地點(diǎn)的高程升高，與洼地邊緣的最低點(diǎn)的高程相同，生成無洼地DEM。

(3)流向計(jì)算

流向計(jì)算基于無洼地DEM，通過計(jì)算中心柵格與鄰域柵格的最大距離權(quán)落差來確定(D8算法)，應(yīng)用ArcHydro模塊下的Flow Direction功能，生成格網(wǎng)流向圖，每個(gè)格網(wǎng)值表示它流向相鄰像元的方向。

(4)匯流累積量計(jì)算

基于無洼地DEM的水流方向的計(jì)算，執(zhí)行ArcHydro模塊中Flow Accumulation命令計(jì)算出每個(gè)格網(wǎng)上游累積匯流數(shù)，得到匯流累積量分布圖。

(5)河網(wǎng)提取

基于匯流累積矩陣，首先執(zhí)行map algebra工具集中的Raster Calculator工具中的con命令;然后設(shè)置閾值，閾值的設(shè)定在河網(wǎng)的提出過程中非常關(guān)鍵，直接影響到河網(wǎng)的提取結(jié)果，本研究經(jīng)多次試驗(yàn)，設(shè)定閾值為10000，其單位是柵格數(shù)，可獲得基于DEM的水系柵格圖;最后進(jìn)行柵格河網(wǎng)失量化，執(zhí)行ArcHydro模塊中的stream to feature命令，得到矢量形式的河網(wǎng)分布圖(圖2)。

圖2 萊州灣入海河網(wǎng)分布圖Fig.2 The distribution map of river network into Laizhou Bay

(6)流域劃分

流域又稱集水區(qū)域，是指流經(jīng)其中的水流和其它物質(zhì)從一個(gè)公共的出水口排出從而形成一個(gè)集中的排水區(qū)域，是給定點(diǎn)的匯水區(qū)域。匯水區(qū)劃分基于水流方向數(shù)據(jù)，執(zhí)行basin命令，獲得匯水區(qū)，如圖3所示。

圖3 萊州灣入海河流匯水區(qū)Fig.3 Rivers into Laizhou Bay catchments

(7)子流域劃分

子流域的劃分是水文模擬的重要一步，也是流域農(nóng)業(yè)污染的重要研究步驟。本研究基于DEM，根據(jù)河流入?？诘奈恢?，結(jié)合實(shí)際情況，從北向南，從東到西，劃分為十一大子流域，包括小島河流域;溢洪河、廣利河和滋脈河(溝)/支脈河流域;小清河流域;彌河流域;白浪河流域;虞河流域;濰河流域;膠萊河流域;沙河流域;王河流域和界河流域。各子流域如圖4所示。

圖4 萊州灣子流域分布圖Fig.4 The distribution map of Laizhou Bay sub-watershed

2.3 子流域匯水區(qū)面積計(jì)算

在劃分出匯水區(qū)和子流域的基礎(chǔ)上，經(jīng)等面積投影，計(jì)算出各入海河流子流域的匯水區(qū)面積，如表1所示。

表1 萊州灣各子流域面積Table 1 The area of Laizhou Bay sub-watershed

3 污染排放系數(shù)與入河系數(shù)

根據(jù)調(diào)查，認(rèn)為海灣農(nóng)業(yè)非點(diǎn)源污染主要來源于農(nóng)業(yè)種植、畜禽養(yǎng)殖污染物、農(nóng)村生活污水排放和海水養(yǎng)殖。因此，在劃分匯水區(qū)和子流域的基礎(chǔ)上，重點(diǎn)研究該區(qū)內(nèi)由農(nóng)田耕作、畜禽飼養(yǎng)、海水養(yǎng)殖、農(nóng)村居民生活所排放的農(nóng)業(yè)非點(diǎn)源污染中COD、NH3－N、TN和TP污染負(fù)荷的流域分配情況。

匯水區(qū)內(nèi)各污染物總量為:

式中，Wi為第 i種(1、2、3、4 分別代表 COD、NH3－N、TN、TP)農(nóng)業(yè)污染負(fù)荷總量，Wij是土地利用類型為j種(1、2、3、4分別代表農(nóng)田、畜禽飼養(yǎng)用地、海水養(yǎng)殖用地、農(nóng)村居住地)污染物的污染負(fù)荷;Si為第j種土地利用類型的面積，其中j=4為匯水中農(nóng)村人口總數(shù);Kij為第j種土地利用類型對應(yīng)第i種污染物的排污系數(shù)［25－29］，其中，農(nóng)田耕作和農(nóng)村居民部分排污系數(shù)參照《全國水環(huán)境容量核定技術(shù)指南》提供的系數(shù)［30］;畜禽飼養(yǎng)用地內(nèi)各畜禽種類排污系統(tǒng)采用國家環(huán)?？偩汁h(huán)發(fā)［2004］43號文件《關(guān)于減免家禽業(yè)排污費(fèi)等有關(guān)問題的通知》中提供的排污系數(shù)［26］;海水養(yǎng)殖排污系數(shù)參見黃渤海海水養(yǎng)殖污染排放成果［31－32］(表 2)。

農(nóng)業(yè)非點(diǎn)源污染入河量是指一定時(shí)期內(nèi)，由地表徑流攜帶進(jìn)入河流等地表水體的污染負(fù)荷。要估算入河量，需要確定入河系數(shù)。入河系數(shù)需要較長期對水質(zhì)和水量的同步監(jiān)測，本文利用已有的研究成果［33－36］(表 3)。

表2 排污系數(shù)Table 2 Pollutant emission coefficients

表3 入河系數(shù)Table 3 Into the river coefficients

4 污染排放總量估算與精度評定

利用2005年的山東省土地利用圖和2010年的海水養(yǎng)殖分布圖(圖5)，土地利用現(xiàn)狀圖包括居民用地、草地、林地、耕地等類型，由于本研究是針對農(nóng)業(yè)非點(diǎn)源污染進(jìn)行估算，因此只考慮農(nóng)用地，而農(nóng)村居住地則通過農(nóng)村居民數(shù)進(jìn)行估算，畜禽飼養(yǎng)用地通過畜禽種類的數(shù)量進(jìn)行估算，匯水區(qū)內(nèi)2005年到2010年間耕地的面積變化不大，對估算結(jié)果的影響在允許誤差范圍內(nèi)，故可作為2010年的污染物估算結(jié)果。把土地利用現(xiàn)狀圖和海水養(yǎng)殖分布圖與匯水區(qū)進(jìn)行疊加，分別計(jì)算不同土地利用類型的匯水區(qū)面積，利用上面的污染排放系數(shù)，估算出各子流域污染排放量。

圖5 萊州灣匯水區(qū)土地利用分布圖Fig.5 the land-use map of Laizhou Bay catchments

經(jīng)過疊加分析，計(jì)算得到匯水區(qū)內(nèi)耕地面積約為2392810 hm2;海水養(yǎng)殖面積為1.7 hm2;根據(jù)山東省2011年統(tǒng)計(jì)年鑒［37］，截至到2010年底，匯水區(qū)內(nèi)農(nóng)村人口數(shù)約為360.85萬人;?？倲?shù)約為529.66萬頭，豬為1106.27萬頭;雞、鴨等家禽約為44825.05萬只。

由上述的公式計(jì)算可得，COD的每年污染總量約為3095461 t;NH3－N的每年污染總量約為297917 t;TN的每年污染總量約為698870 t;TP的每年污染總量約為217118 t。再根據(jù)表3中的各污染物的入河系數(shù)，計(jì)算得到 COD、NH3－N、TN 和 TP每年的入河總量分別為:236933、23956、53684 t和 15922 t。利用子流域所占的面積，分別計(jì)算出各子流域的COD、NH3－N、TN 和 TP 每年的入河總量(表4)。

目前還沒有見到針對萊州灣污染總量的估算，相應(yīng)的污染監(jiān)測也沒有形成每年入灣的農(nóng)業(yè)非點(diǎn)源污染總量數(shù)據(jù)，不能與實(shí)際的污染總量進(jìn)行直接的對比來評定本文的估算精度，但從已有對小清河流域氮素污染總量的研究成果來看，朱虎［38］計(jì)算2006年小清河流域年均氮失潛力為10.44×103—36.86×103t，平均23.65×103t，本文估算的小清河流域全氮和氨氮的總?cè)牒Ａ繛?5477 t，在其范圍值內(nèi)，基本與平均量相同，比較接近實(shí)際情況，其精度具有一定的可靠性。

表4 萊州灣各子流域入海主要污染物排放量/tTable 4 Each sub-watershed emissions of major pollutants into Laizhou Bay

5 結(jié)論

本文借助GIS技術(shù)，提出一種基于流域單元的農(nóng)業(yè)非點(diǎn)源污染負(fù)荷的估算方法，并以萊州灣為例，分別估算了由農(nóng)業(yè)污染引起的非點(diǎn)源污染物COD、NH3－N、TN和TP入海污染總量。上述4種農(nóng)業(yè)非點(diǎn)源污染物的定量化估算表明，萊州灣農(nóng)業(yè)非點(diǎn)源污染主要是由農(nóng)業(yè)生產(chǎn)過程所產(chǎn)生的;從主要污染物角度來看，4種主要污染物中COD入海量最大，TN次之，然后是 NH3－N 和 TP，COD、NH3－N、TN 和 TP 的排放量占總排放量的比例分別為71.69%、7.25%、16.24%和4.82%，這為萊州灣環(huán)境污染原因診斷和調(diào)控措施的制定更有針對性。該研究通過探索適合的流域農(nóng)業(yè)非點(diǎn)源污染負(fù)荷計(jì)算方法，來識(shí)別流域非點(diǎn)源污染的負(fù)荷量及其空間分布，從而為相應(yīng)流域管理措施的制定和海洋環(huán)境的調(diào)控提供依據(jù)。

雖然對萊州灣入海非點(diǎn)源負(fù)荷進(jìn)行了定量估算，但由于缺乏詳細(xì)的、長序列的、以及農(nóng)業(yè)非點(diǎn)源污染專項(xiàng)監(jiān)測資料，并且海灣污染是一個(gè)非常復(fù)雜的機(jī)理過程，其農(nóng)業(yè)非點(diǎn)源污染影響因素眾多，各污染系數(shù)主要采用文獻(xiàn)中的方法確定，并不一定完全適應(yīng)于萊州灣區(qū)域，今后應(yīng)更多地通過典型區(qū)域的污染發(fā)生學(xué)試驗(yàn)研究來確定各類系數(shù)的值;其次由于DEM的分辨率為90 m，缺乏部分特殊地形的細(xì)部特征，在匯水區(qū)劃分上不夠精細(xì)?？梢灶A(yù)計(jì)，隨著人們對非點(diǎn)源污染機(jī)制、過程的深入理解，海灣農(nóng)業(yè)非點(diǎn)源污染必將獲得進(jìn)一步發(fā)展，為流域非點(diǎn)源的管理、控制提供更好的技術(shù)支持。

我國海灣區(qū)域農(nóng)業(yè)非點(diǎn)源污染控制還處于探索階段，國家尚未進(jìn)行全面調(diào)控，缺乏對造成海灣污染的農(nóng)業(yè)非點(diǎn)源進(jìn)行全面、系統(tǒng)的認(rèn)識(shí)和研究，但非點(diǎn)源污染對于近岸海灣主體的影響日益突出，因此，正確的引導(dǎo)政府管理部門和公眾提高對非點(diǎn)源污染的重視并采取有效的調(diào)控措施，對保護(hù)海灣水環(huán)境意義重大。期望進(jìn)一步開展海灣非點(diǎn)源污染總量控制、管理計(jì)劃等方面的研究。

［1］Dennis L C，Peter J V，Keith L.Modeling non－point source pollution in vadose zone with GIS.Environmental Science and Technology，1997，8:2157－2175.

［2］Miller G T.Living in the Environment:An Introduction to Environmental Science.17th ed.Belmont:Wadsworth Publishing Company，1992:602－611.

［3］Boers P C M. Nutrient emissions from agriculture in the netherlands: causes and remedies. Water Science and Technology，1996，33(4/5):22－25.

［4］Huang H，Zou CW，Chen X G.Review of the researches on nonpoint source pollution in China.Ecology and Environment，2004，13(2):255－257.

［5］Donigian A S，Bichnell B R，Patwardhan A S，Linker L C，Chang C H，Reynolds R，Carsel R.Chesapeake Bay Program Watershed Model Application to Calculate Bay Nutrient Loadings－Final Facts and Recommendations.Report No.EPA 903－R－94－042.Annapolls，Maryland:USEPA Chesapeake Bay Program Office，1994.

［6］Thomton J A.Assessment and Control of Non－Point Source Pollution of Aquatic Ecosystems:A practical Approach.Pearl River，New York:The Parthenon Parthenon Publishing Group，1999.

［7］Prink C R.Estimating nutrient exports to estuaries.Journal of Environ Quality，1999，20(4):717－724.

［8］Liu J H，Lai G Y.Research progress of the agricultural non－point source pollution.Journal of Water Resources＆ Water Engineering，2007，18(1):29－32.

［9］Liang B，Wang X Y，Cao L P.Water environment non－point source pollution loading estimation methods in China.Journal of Jilin Normal University:Natural Science Edition，2004，25(3):58－61.

［10］Gilliland M W，Baxter－Potter W.A geographic information system to predict non－point source pollution potential.Water Resources Bulletin，1987，23(2):281－291.

［11］ Yong R A，Onstad CA，Bosch D D，Anderson WP.AGNPS:A nonpoint－source pollution model for evaluating agricultural watersheds.Journal of Soil Water Conservation，1989，44(2):168－173.

［12］Bouraoui F，Dllaha T A.Answers－2000:runoff and sediment transport model.Journal of Environmental Engineering，1996，122(6):493－502.

［13］Sivertun ?，Prange L.Non－point source critical area analysis in the Gissel?watershed using GIS.Environmental Modelling＆Software，2003，18(10):887－898.

［14］Cai M，Li H N，Zhuang Y T.Rainfall deduction method for estimating non－point source pollution load for watershed.Journal of Northwest Sci－Tech University of Agriculture and Forestry:Natural Science Edition，2005，33(4):102－106.

［15］Wang L J.Application of partial least square regress method to estimate the annual load of non－point source pollution.Resources and Environment in the Yangtze Basin，2007，16(6):810－813.

［16］Chen Y Y，Hui E Q，Jin CJ，Qiu H X，Wu D X.A hydrological method for estimation of non－point source pollution loads and its application.Research of Environmental Sciences，2003，16(1):10－13.

［17］Li H E.Mean concentration method for estimation of nonpoint source load and its application.Acta Scientiae Circumstantiae，2000，20(4):397－400.

［18］Hong X K，Li H E.Correlation method of water quality and quantity and its application to load estimation of nonpoint source pollution.Journal of Xi'an University of Technology，2000，16(4):384－386.

［19］Yang Y H，Shen W B.Preliminary study on estimating surface water non－point source pollution loads.Journal of Jilin University:Earth Science Edition，2006，36(S1):105－108.

［20］Li QK，Li H E，Hu Y W，Sun J.Unit－analysis－based assessment of agricultural non－point source pollution in Qingtongxia Irrigation Area.Journal of Ecology and Rural Environment，2007，23(4):33－36.

［21］Li H E，Zhuang Y T.The export coefficient modeling approach for load prediction of nutrient from nonpoint source and its application.Journal of Xi'an University of Technology，2003，19(4):307－310.

［22］Yu T，Meng W，Ongley E，Zheng B H，Deng Y X.Problems and recommendations for non－point source pollution identification in China.Acta Scientiae Circumstantiae，2008，28(3):401－407.

［23］Cai M，Li H E，Zhuang Y T，Wang QH.Application of modified export coefficient method in polluting load estimation of non－point source pollution.Journal of Hydraulic Engineering，2004，(7):40－45.

［24］Shandong Provincial Oceanic and Fishery Information Network.The Laizhou Bay ecological demonstration project planning.2011.http://www.mingbo－aquatic.com/showarticle.asp?newsid=160

［25］Dai X Y，Guo Z Y，Shi C，Wu J P.Study on watershed distribution of agricultural non－point source pollution in typical areas of Shanghai.Environmental Pollution ＆ Control，2004，26(4):278－280.

［26］The State Environmental Protection Administration of China.For the Reduction of the Poultry Industry Sewage Related Issues Notice. Beijing: The State Environmental Protection Administration of China，2004.

［27］Lin Z L，Yu X G，Liu X，Chen B.Estimation of COD discharge from non－point pollution into Luoyuan Bay. Journal of Oceanography in Taiwan Strait，2006，25(2):209－215.

［28］Yang SJ，Zhang A P，Yang Z L，Yang S Q.Agricultural nonpoint source pollution in Ningxia irrigation district and preliminary study of load estimation methods.Scientia Agricultura Sinica，2009，42(11):3947－3955.

［29］Yuan X T，Zhang M J，Liang Y B，Liu D，Guan D M.Selfpollutant loading from a suspension aquaculture system of Japanese scallop(Patinopecten yessoensis) in the Changhai sea area，Northern Yellow Sea of China.Aquaculture，2010，304(1/4):79－87.

［30］Chinese Academy for Environmental Planning.National Water Environmental Capacity Verification Manual.Beijing:Chinese Academy for Environmental Planning，2003.

［31］Cui Y，Chen B J，Chen JF.Evaluation of self－pollution of marine culture in the Yellow Sea and Bohai Sea.Chinese Journal of Applied Ecology，2005，16(1):180－185.

［32］Huang X Q，Yao Y M，Wang JH，Jiang X S，Shi SH，Chen L.Leqing Bay Marine Environmental Capacity and Total Pollutants Control Research.Beijing:China Ocean Press，2011.

［33］Funge－Smith S J，Briggs M R P.Nutrient budgets in intensive shrimp ponds:implications for sustainability.Aquaculture，1998，164(1/4):117－133.

［34］Xin Z W，Sun R，Lu X Q.Comprehensive Analysis and Management Strategy for Regional Water Environment GEF in Tianjin.Beijing:China Environmental Science Press，2009.

［35］Zhu M.Study on Agricultural NPS Loads of Haihe Basin and Assessment on Its Environmental Impact［D］.Beijing:Chinese Academy of Agricultural Sciences，2011.

［36］Ren X S.The Haihe River Basin Water Resources Evaluation.Beijing:China Water Power Press，2007.

［37］Shandong Provincial Bureau of Statistics，National Bureau of Statistics in Shandong Investigated Teams.Shandong Statistical 2011 Year Book.Beijing:China Statistics Press，2011.

［38］Zhu H.Quantifying Analysis on Non－Point Source Nitrogen Pollution in Xiaoqing River Basin of Shandong Province［D］.Beijing:Chinese Academy of Agricultural Sciences，2009.

參考文獻(xiàn):

［4］黃虹，鄒長偉，陳新庚.中國非點(diǎn)源污染研究評述.生態(tài)環(huán)境，2004，13(2):255－257.

［8］劉紀(jì)輝，賴格英.農(nóng)業(yè)非點(diǎn)源污染研究進(jìn)展.水資源與水工程學(xué)報(bào)，2007，18(1):29－32.

［9］梁博，王曉燕，曹利平.我國水環(huán)境非點(diǎn)源污染負(fù)荷估算方法研究.吉林師范大學(xué)學(xué)報(bào):自然科學(xué)版，2004，25(3):58－61.

［14］蔡明，李懷恩，莊詠濤.估算流域非點(diǎn)源污染負(fù)荷的降雨量差值法.西北農(nóng)林科技大學(xué)學(xué)報(bào):自然科學(xué)版，2005，33(4):102－106.

［15］王玲杰.偏最小二乘回歸法在非點(diǎn)源污染年負(fù)荷量估算中的應(yīng)用.長江流域資源與環(huán)境.2007，16(6):810－813.

［16］陳友媛，惠二青，金春姬，邱漢學(xué)，吳德星.非點(diǎn)源污染負(fù)荷的水文估算方法.環(huán)境科學(xué)研究，2003，16(1):10－13.

［17］李懷恩.估算非點(diǎn)源污染負(fù)荷的平均濃度法及其應(yīng)用.環(huán)境科學(xué)學(xué)報(bào)，2000，20(4):397－400.

［18］洪小康，李懷恩.水質(zhì)水量相關(guān)法在非點(diǎn)源污染負(fù)荷估算中的應(yīng)用.西安理工大學(xué)學(xué)報(bào)，2000，16(4):384－386.

［19］楊育紅，沈萬斌.地表水非點(diǎn)源污染負(fù)荷計(jì)算方法探討.吉林大學(xué)學(xué)報(bào):地球科學(xué)版，2006，36(S1):105－108.

［20］李強(qiáng)坤，李懷恩，胡亞偉，孫娟.基于單元分析的青銅峽灌區(qū)農(nóng)業(yè)非點(diǎn)源污染估算.生態(tài)與農(nóng)村環(huán)境學(xué)報(bào)，2007，23(4):33－36.

［21］李懷恩，莊詠濤.預(yù)測非點(diǎn)源營養(yǎng)負(fù)荷的輸出系數(shù)法研究進(jìn)展與應(yīng)用.西安理工大學(xué)學(xué)報(bào)，2003，19(4):307－310.

［22］于濤，孟偉，Ongley E，鄭丙輝，鄧義祥.我國非點(diǎn)源負(fù)荷研究中的問題探討.環(huán)境科學(xué)學(xué)報(bào)，2006，28(3):401－407.

［23］蔡明，李懷恩，莊詠濤，王清華.改進(jìn)的輸出系數(shù)法在流域非點(diǎn)源污染負(fù)荷估算中的應(yīng)用.水利學(xué)報(bào)，2004，(7):40－45.

［24］山東省海洋與漁業(yè)廳.萊州灣生態(tài)整治示范工程規(guī)劃.2011.http://www.mingbo－aquatic.com/showarticle.asp?newsid=160

［25］戴曉燕，過仲陽，石純，吳建平.上海市青浦區(qū)農(nóng)業(yè)非點(diǎn)源污染的流域分配研究.環(huán)境污染與防治，2004，26(4):278－280.

［26］國家環(huán)境保護(hù)總局.關(guān)于減免家禽業(yè)排污費(fèi)等有關(guān)問題的通知.北京:國家環(huán)境保護(hù)總局，2004.

［27］林志蘭，余興光，劉希，陳彬.羅源灣區(qū)域入海非點(diǎn)源污染COD排放量估算.臺(tái)灣海峽，2004，25(2):209－215.

［28］楊淑靜，張愛平，楊正禮，楊世琦.寧夏灌區(qū)農(nóng)業(yè)非點(diǎn)源污染負(fù)荷估算方法初探.中國農(nóng)業(yè)科學(xué)，2009，42(11):3947－3955.

［30］中國環(huán)境規(guī)劃院.全國水環(huán)境容量核定技術(shù)指南.北京:中國環(huán)境規(guī)劃院，2003.

［31］崔毅，陳碧鵑，陳聚法.黃渤海海水養(yǎng)殖自身污染的評估.應(yīng)用生態(tài)學(xué)報(bào)，2005，16(1):180－185.

［32］黃秀清，姚炎明，王金輝，蔣曉山，石少華，陳雷.樂清灣海洋環(huán)境容量及污染物總量控制研究.北京:海洋出版社，2011.

［34］辛志偉，孫韌，盧學(xué)強(qiáng).區(qū)域水環(huán)境綜合解析與管理策略－GEF在天津.北京:中國環(huán)境科學(xué)出版社，2009.

［35］朱梅.海河流域農(nóng)業(yè)非點(diǎn)源污染負(fù)荷估算與評價(jià)研究［D］.北京:中國農(nóng)業(yè)科學(xué)院，2011.

［36］任憲韶.海河流域水資源評價(jià).北京:中國水利水電出版社，2007.

［37］山東省統(tǒng)計(jì)局，國家統(tǒng)計(jì)局山東調(diào)查總隊(duì).山東統(tǒng)計(jì)年鑒2011.北京:中國統(tǒng)計(jì)出版社，2011.

［38］朱虎.小清河流域農(nóng)田非點(diǎn)源氮污染定量評價(jià)研究［D］.北京:中國農(nóng)業(yè)科學(xué)院，2009.