北京市農業面源污染負荷特征分析及控制分區

謝培，宮健,陳誠

1.環境基準與風險評估國家重點實驗室，中國環境科學研究院 2.首都師范大學資源環境與旅游學院

在城鎮化大背景下，水環境中來自城市和工業點源的污染問題逐步得到解決，而農業面源污染由于具有來源廣泛性、隨機性、排放不確定性、時空分布不均性，已經成為我國當前水環境污染的重要原因之一[1-2]。進行農業面源污染源解析，準確掌握農業面源污染負荷，并開展面源污染控制區劃具有重要意義。已有學者對太湖流域[3]、三峽庫區[4]、密云水庫流域[5-6]等的農業面源污染進行了定性和定量研究，圍繞污染負荷估算和空間分布特征[7]、評價方法及源解析[8-9]、污染類型劃分及防治[10-11]等開展了研究。污染負荷主要的估算方法包括監測法[12]、輸出系數法[13]、清單分析法[14]及模型模擬法[15]等，通常根據研究區數據資料的完整程度加以選擇。其中，清單分析法具有簡單、方便、易于操作和實用等優點，適用于不同尺度區域污染負荷核算，已得到了廣泛的應用[16]。

大量文獻調研發現，北京市各區自然環境和經濟社會發展各具特色，農業面源污染特征差異較大。由于調查難度大，已有對北京市農業面源污染的研究多集中于小流域，對整個北京市農業面源污染時空特征的研究相對缺乏。筆者運用清單分析、等標負荷和聚類分析等研究方法，在區級尺度上對北京市農業面源污染負荷及其空間特征、污染分類控制等進行綜合研究，確定各區主要污染源、污染物與污染空間分布特征，并進行面源污染控制二級分區，以期為北京市農業面源污染防治提供參考。

1 數據來源與處理方法

1.1 研究區概況

北京市地處海河流域，位于華北平原北部，毗鄰渤海灣，總面積為16 410.54 km2，其中山區面積為10 418 km2，約占全市總面積的62%。屬暖溫帶半濕潤季風型大陸性氣候區，多年平均降水量為595 mm，人均水資源較為缺乏；地勢西北高東南低，過渡較為急劇；土壤類型主要為褐土和潮土，地帶性植被類型為落葉闊葉林。境內有北運河、潮白河、大清河、永定河、薊運河五大水系及密云水庫，長期以來高強度的工農業生產造成地表水水質惡化，其中北運河、潮白河水系部分河段污染較為嚴重。北京市下轄16個區，其中朝陽、豐臺、海淀、房山、通州、順義、昌平、大興、門頭溝、懷柔、平谷、密云、延慶這13個區涉及農業，2016年農村人口為291.7萬人，農業總產值為145.2億元；北京市各郊縣種植面積較大，農業生產發達，2016年糧食播種面積為87 328.7 hm2，糧食產量達54.17萬t。

1.2 數據來源

北京市涉及農業的13個區的基礎數據包括國土面積、耕地面積、地表水資源量、農作物產量、農作物播種面積、鄉村人口數、畜禽養殖產量等來自《北京統計年鑒2016》[17]、《北京區域統計年鑒2016》[18]。

1.3 數據處理方法

1.3.1污染負荷估算方法

采用清單分析法建立產排污清單，以此估算各類污染源的污染物排放量，具體計算公式如下：

E=∑SUi×ρi×LCi

(1)

式中：E為進入水體的農業面源污染物總氮(TN)、總磷(TP)、化學需氧量(COD)的排放總量，t/a；SUi為i污染單元污染物產生基數，t/a；ρi為i污染單元產污強度系數；SUi×ρi表示污染單元農業面源污染物產生量，即不考慮資源綜合利用和管理因素時農業生產和農村生活造成的最大潛在農業面源污染物產生量，t/a；LCi為考慮資源綜合利用和管理因素時i污染單元污染物排放系數。上述參數通過相關文獻[19-23]分析并結合北京市典型區域內的調研結果確定。

1.3.2污染負荷特征分析

采用等標污染負荷衡量各類污染物(TN、TP、COD)或污染源對環境潛在污染能力的大小，以確定主要污染源或污染物。計算公式如下：

Pj=Cj/C0

(2)

式中：Pj為污染物j的等標排放量，m3；Cj為污染物j的絕對排放量，t/a；C0為污染物j基于水環境功能分區的水質目標控制類別(GB 3838—2002《地表水環境質量標準》)標準值，mg/L。

等標污染指數指以各地區污染源的等標排放量除以該地區地表水資源量，用以反映該地區水體存在的潛在污染程度。等標污染指數參照文獻[24]進行計算，公式如下：

Mj=PjV

(3)

式中：Mj為污染源或污染物j的等標污染指數；V為地表水資源量，m3。

等標相對排放系數的大小可以反映農業面源污染的排放強度，基于耕地面積、農村人口以及農業產值計算北京市5類農業面源污染源(化肥施用、作物秸稈、畜禽養殖、水產養殖、農村生活)的等標相對排放系數。不同污染源等標相對排放系數的計算公式如下：

QI=QabS

(4)

式中：QI為第a區第b類污染源的等標相對排放系數，t/(km2·a)；Qab為第a區第b類污染源的等標絕對排放量，t/a；S為第a區的耕地面積，km2。

2 結果與討論

2.1 北京市農業面源污染負荷及空間分布

2.1.1農業面源污染源負荷估算

北京市5類農業面源污染源TN、TP、COD排放量及占比如表1所示。由表1可知，北京市農業面源污染物TN、TP、COD的總排放量分別為18 169.00、2 668.45和70 211.10 t/a。農村生活為COD的主要來源，其次為畜禽養殖、作物秸桿；畜禽養殖、農村生活和化肥施用為TN的主要來源；農村生活為TP的主要來源，其次是畜禽養殖，這與TN的來源排放特征較為一致。與上海市農業面源污染源排放特征[11]相比，北京市農業面源TP排放量低于上海市(4 400 t/a)，而TN排放量高于上海市(11 300 t/a)。2016年北京市城市生活污水COD排放量為79 396 t/a，工業廢水COD排放量為4 738 t/a[17]，而農業面源污染COD排放量約占北京市COD總排放量的45.49%，可見農業面源污染負荷排放對水體的貢獻很大。

表1 北京市5類農業面源污染源TN、TP、COD排放量及占比

2.1.2農業面源污染負荷空間分布

北京市13個區農業面源污染物TN、TP與COD排放量及占比見表2，由地理信息系統(GIS)軟件做出污染物排放量空間分布如圖1所示。

由表2、圖1可知，13個區農業面源污染負荷呈現出較強的區域分異特征。從污染物排放總量來看，北運河水系周邊，潮白河下游段和大清河水系附近的房山、通州、順義、大興、平谷等區較高，合計占北京市農業面源污染排放總量的50%以上。從污染物排放空間分布來看，北運河和永定河水系周邊的通州、順義、大興區的TN排放量較高；COD排放量高值區主要集中分布在永定河水系、北運河周邊、潮白河下游段、薊運河水系附近，包括房山、通州、順義、大興、平谷區；TP排放量低值區出現在豐臺、朝陽、海淀、門頭溝區，高值區出現在順義、通州、大興、平谷區。

表2 北京市13個區農業面源污染物TN、TP、COD排放量及占比

(續表2)

圖1 北京市13個區農業面源污染物TN、TP及COD排放量空間分布Fig.1 Spatial distribution of TN, TP and COD emissions from agricultural NPS in 13 districts of Beijing City

2.1.3不同污染源絕對排污貢獻率

北京市13個區農業面源污染物TN、TP、COD不同污染源排放占比如圖2所示。由圖2結合表1可知，畜禽養殖、農村生活對農業面源中TN和COD的排放貢獻率均較高。從TN來看，畜禽養殖、農村生活和化肥施用是主要來源，朝陽、豐臺、海淀、昌平、門頭溝等區TN排放主要來自農村生活；從TP來看，農村生活和畜禽養殖是主要來源，朝陽、豐臺、海淀、昌平、門頭溝等區TP排放主要來自于農村生活；從COD來看，農村生活、畜禽養殖和作物秸稈是主要來源，豐臺、海淀、朝陽、門頭溝等區的COD排放幾乎全部來自于農村生活，占比均超過80%。除平谷和朝陽區水產養殖對TN和TP排放量貢獻較高外，其他各區水產養殖對COD、TN、TP的排放貢獻率均較小。順義、通州、大興、懷柔、延慶等區TN排放主要來自化肥施用和畜禽養殖，這與王磊等[25]指出的順義、大興、通州等區的化肥污染風險較高的結果相一致。

圖2 北京市13個區農業面源污染物TN、TP及COD不同污染源排放貢獻率Fig.2 Discharge contribution rate of TN, TP and COD of agricultural NPS from different pollution sources in 13 districts of Beijing City

2.2 北京市農業面源污染負荷分析

2.2.1等標負荷特征

北京市農業面源污染物TN、TP、COD的等標排放量與等標排放比如表3所示。由表3可知，北京市農業面源的主要污染物為TN和TP，其等標排放量分別為1 5176.15和9 393.78 t/a，分別占農業面源污染物等標排放量的55.15%和34.14%，遠高于COD的等標排放量。等標排放量、等標排放比大小均為畜禽養殖>農村生活>化肥施用>水產養殖>作物桔桿。可見，北京市農業面源主要污染源為畜禽養殖、農村生活、化肥施用，其等標排放比分別為33.30%、27.33%和19.55%，三者累計污染負荷比高達80%以上。

表3 北京市農業面源污染物的等標排放量與等標排放比

北京市13個區農業面源不同污染源的等標排放量、等標排放比分別如表4、表5所示。由表4可知，北京市13個區農業面源等標排放量平均為2 243.40 t/a，其中朝陽、豐臺、海淀、門頭溝、懷柔5個區的等標排放量低于北京市的平均水平，豐臺區的等標排放量最小，為77.18 t/a；順義區的等標排放量最大，為4 283.43 t/a。由表5可知，北京市13個區化肥施用、農村生活、畜禽養殖、水產養殖、作物秸稈的等標排放比平均值分別為17.30%、9.24%、28.06%、10.98%、35.05%，其中畜禽養殖和農村生活的等標排放比最大。從等標排放比來看，各區的農業面源污染來源差異較大，呈現不同的來源特征：海淀、門頭溝這2個區的農業面源污染主要來自農村生活，其等標排放比分別為53.68%、73.65%；順義、懷柔、密云、平谷、延慶等區的農業面源污染主要來自于畜禽養殖、農村生活，二者的等標排放比之和達50%以上；通州、大興等區的農業面源污染主要來自于農村生活、化肥施用，二者的等標排放比超過50%。

表5 北京市13個區農業面源不同污染源的等標排放比

表4 北京市13個區農業面源不同污染源的等標排放量

2.2.2污染源排放強度

盡管等標排放量可以將不同污染物在同一尺度上進行比較，但它考慮的仍然是污染負荷總量。農業面源污染受區域大小、農村人口、農業產值等因素的影響，采用基于耕地面積的等標相對排放系數更能體現農業面源污染排放強度的可比性[26-28]。計算出北京市基于單位耕地面積的等標相對排放系數為3.14 t/(km2·a)，由此計算北京市不同污染源的等標相對排放系數，結果如表6所示。由表6可知，北京市5類污染源的等標相對排放系數表現為農村生活>畜禽養殖>化肥施用>水產養殖>作物秸稈，農村生活的等標相對排放系數最大，為主要污染源，其次為畜禽養殖。可見，在排除了區域耕地面積對農業面源污染物絕對等標排放量的影響后，北京市農業面源主要污染源為農村生活和畜禽養殖。

表6 北京市13個區基于耕地面積的不同污染源的等標相對排放系數

2.2.3等標污染指數

北京市13個區農業面源污染源和污染物的等標污染指數分別如圖3、圖4所示。由圖3可知，5類污染源等標污染指數平均值以農村生活為最高，作物秸稈為最低。各區的等標污染指數差異較大，農村生活的等標污染指數在房山區最高，達6.95；化肥施用等標污染指數在通州、房山、順義、延慶、大興等區較高，分別為3.18、3.14、2.79、2.35、2.33；畜禽養殖等標污染指數在房山、順義、平谷等區較高，分別為6.80、5.97、4.86；水產養殖等標污染指數在平谷區最高，達2.13，水產養殖和作物秸稈的等標污染指數在各區均較低。由圖4可知，TN、TP、COD的等標污染指數的平均值分別為4.22、3.19、0.87。各區TN、TP、COD的等標污染指數均以房山區為最大，分別為9.82、7.51、2.10。通州、房山、順義、大興、平谷、延慶等區的TN、TP、COD的等標污染指數均超過全市平均水平，屬于污染程度較為嚴重的地區。

圖3 北京市13個區農業面源污染源的等標污染指數Fig.3 Equivalent pollution index of agricultural NPS pollution sources in 13 districts of Beijing City

圖4 北京市13個區農業面源污染物的等標污染指數Fig.4 Equivalent pollution index of agricultural NPS pollutants in 13 districts of Beijing City

2.3 北京市農業面源污染空間特征及控制分區

對北京市13個區農業面源污染等標污染指數進行空間聚類分析，結果如圖5所示，由GIS軟件做出北京市農業面源污染風險程度分區如圖6所示。由圖5、圖6可知，將北京市農業面源污染風險程度劃分為3個等級，分別為高度、中度、輕度風險區。高度風險區集中分布在北京市通州、順義、房山和平谷區，其平均等標污染指數為15.33；中度風險區主要分布在大興、延慶、昌平、懷柔、密云區，其平均等標污染指數為7.71；低風險區分布在朝陽、豐臺、門頭溝、海淀區，其平均等標污染指數為1.98。

圖5 北京市13個區農業面源污染風險評價聚類分析Fig.5 Cluster analysis of agricultural NPS pollution risk assessment in 13 districts of Beijing City

圖6 北京市農業面源污染風險程度分區Fig.6 Risk zoning of agricultural NPS pollution in Beijing City

基于一級風險分區結果，為實現污染風險分級控制和污染源的分類管理，結合各區農業面源主要污染源類型(表5)，進行農業面源污染的二級分區，結果如表7所示。由表7可知，北京市農業面源污染控制可分為3個一級區和8個二級區(表7)。

表7 北京市農業面源污染控制分區

2.4 北京市農業面源污染控制對策

2.4.1農村生活污染型控制對策

由農村生活引起的農業面源污染覆蓋北京市大部分郊區，這些區域地勢起伏，鄉村人口眾多，聚居地分布較為分散，沒有完善的排污管網設施；并且農村居民的水環境保護意識淡薄，污水和生活垃圾隨意排放或丟棄現象嚴重。因此，應加快農村地區污水處理廠建設，因地制宜地規劃建設污水處理廠，建廠規模應與鄉村規模和人口總量相匹配，同時建設配套污水管網，確保污水主管網全覆蓋；加快推進農村改廁工作，結合末端集中治污，推廣和實施綠色鄉村清潔工程建設，加快農村生活污水及廢物的資源轉化與綜合利用。

2.4.2畜禽養殖污染型控制對策

北京市大部分區的畜禽養殖業較為發達，為北京市農業面源污染貢獻率最高的污染源類型，尤其是北部和中部(順義、房山、密云、平谷、昌平、懷柔等區)地區的畜禽養殖污染物排放量大。因此，應建立改善畜禽污染、促進畜牧業持續發展的良性養殖機制預防和治理畜禽養殖污染。相關職能部門應通過各種科普形式，宣傳引導養殖戶科學養殖，重點建設循環養殖鏈，實行清潔養殖、生態養殖，促進農村畜禽養殖場污染的減量化和資源化，實現集中統一處置。

2.4.3農田化肥污染型控制對策

北京市農田主要分布在通州、延慶、順義、大興等區，約占全市面積的34.30%。由于種植面積較大，農田化肥施用強度高及施用比例不合理，農產品種植結構與布局不規范，而農田化肥的有效利用率較低，造成農田土壤中肥料養分流失率較高。可實行測土配方施肥，提前測試分析土壤肥力狀態，按照農田作物產量和作物生長的營養物質需要，決定施用肥料的種類和數量。另外，應根據土壤性質和肥料的特點，合理選擇施肥方式，大范圍普及農田測土配方施肥技術和土壤有機質提升技術，制定合理的施肥方式，精確施肥，保證化肥施用效力最大化和化肥養分流失最小化。

3 結論

(1)2016年北京市農業面源污染物TN、TP、COD排放量分別為18 169.00、2 668.45、70 211.10 t/a，TN、TP是主要污染物；北京市農業面源首要污染源為農村生活，其次為畜禽養殖，建議北京市未來農業面源管控中加強對農村生活和畜禽養殖污染的控制。

(2)北京市13個區農業面源污染呈現較強的區域分異特征。水產養殖等標污染指數以平谷區為最高，農村生活等標污染指數以房山區為最高。通州、房山、順義、大興和平谷等區的TN、TP和COD等標污染指數超過全市平均水平，房山區的等標污染指數最大。

(3)可將北京市劃分為高度、中度和輕度三類農業面源污染風險一級分區，根據各區農業面源污染類型進一步劃分為8個二級分區，由此可實現對不同區不同污染源類型的分類管控。