龔世飛，肖能武，丁武漢，居學(xué)海，吳平華，余永松，李虎

基于耕地承載力的十堰市畜禽養(yǎng)殖環(huán)境風險評價

龔世飛1，肖能武1，丁武漢2，居學(xué)海3，吳平華1，余永松1，李虎2

1丹江口庫區(qū)十堰生態(tài)農(nóng)業(yè)研究院，湖北十堰 442000；2中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所，北京 100081；3農(nóng)業(yè)農(nóng)村部農(nóng)業(yè)生態(tài)與資源保護總站，北京 100125

【目的】十堰市是南水北調(diào)中線工程的核心水源區(qū)，全面認識畜禽養(yǎng)殖污染現(xiàn)狀及其對耕地的潛在威脅，為加強畜禽糞便風險管控提供依據(jù)?！痉椒ā窟\用畜禽糞便環(huán)境風險評價模型，對十堰市各縣（區(qū)、市）2020年主要畜禽糞便來源結(jié)構(gòu)進行定量分析，并從種養(yǎng)平衡視角對其環(huán)境風險進行綜合評價。【結(jié)果】十堰市2020年主要畜禽養(yǎng)殖糞便豬糞當量達306.85×104t，且主要分布于鄖陽區(qū)、丹江口市和房縣，糞便中有機物的化學(xué)需氧量（COD）最大（57.00×104t），總氮（TN）（2.63×104t），總磷（TP）含量較少（0.65×104t）。十堰市畜禽養(yǎng)殖等標污染負荷總量為15.532×1010m3，且在不同地區(qū)的來源結(jié)構(gòu)差異明顯，TP是畜禽養(yǎng)殖的首要污染物，其污染等標負荷占負荷總量的41.68%，其次為TN，污染負荷率為33.85%，COD的等標污染負荷率最小，為24.47%。養(yǎng)牛是十堰地區(qū)畜禽養(yǎng)殖面源污染的首要污染源，對等標污染負荷總量的貢獻率達49.47%，其次為家禽（21.60%）。環(huán)境污染風險評價結(jié)果顯示，全市單位耕地畜禽糞便負荷強度為10.15 t·hm-2，警報值為0.338，分級指數(shù)為Ⅰ，對環(huán)境構(gòu)成污染的威脅為“無”。TN、TP的負荷強度分別為0.09和0.02 t·hm-2，環(huán)境風險指數(shù)為1.023和1.223，但區(qū)域間存在明顯的空間分異特征。全市有超過1/2地區(qū)存在潛在的畜禽養(yǎng)殖污染風險，其中以茅箭區(qū)和鄖陽區(qū)最為突出。以磷（P）為衡量標準，十堰市養(yǎng)殖總量須在現(xiàn)有基礎(chǔ)上消減98.51萬頭豬當量?！窘Y(jié)論】保持合理的畜禽養(yǎng)殖規(guī)模（≤441.00萬頭豬當量），同時做好污染物消減措施，對促進十堰市種養(yǎng)平衡和畜禽養(yǎng)殖業(yè)有序發(fā)展具有重要意義。

十堰市；耕地承載力；畜禽糞便；耕地負荷；環(huán)境風險

0 引言

【研究意義】長期以來，畜禽養(yǎng)殖業(yè)一直是我國農(nóng)業(yè)農(nóng)村經(jīng)濟發(fā)展的重要支柱產(chǎn)業(yè)。尤其是近年來，隨著社會經(jīng)濟的迅猛發(fā)展，畜禽養(yǎng)殖業(yè)的經(jīng)營方式也正由分散養(yǎng)殖向規(guī)?；?、集約化和產(chǎn)業(yè)化發(fā)展方向轉(zhuǎn)變[1]，我國的肉類、蛋類和禽類產(chǎn)量經(jīng)過連續(xù)多年穩(wěn)定增長已躍居全球首位[2-3]。然而，在保證畜禽產(chǎn)品有效供給的同時，迅速增加的養(yǎng)殖糞污對環(huán)境構(gòu)成的威脅也在不斷攀升，已經(jīng)成為農(nóng)業(yè)面源污染的重要成因之一[4-5]。根據(jù)2020年發(fā)布的第二次全國污染源普查公報結(jié)果顯示，我國畜禽養(yǎng)殖業(yè)排放的化學(xué)需氧量（COD）、總氮（TN）、總磷（TP）分別為1 000.53×104、53.63×104、11.97×104t，分別占全國排放總量的46.67%、19.61%、37.95%，占農(nóng)業(yè)源污染物排放量的93.76%、42.14%和56.46%，總體占比較高。因此，準確認識畜禽養(yǎng)殖業(yè)的污染來源、負荷強度和分布特征及其潛在環(huán)境風險，對有針對性地提出養(yǎng)殖業(yè)面源污染防控策略，進而推動區(qū)域農(nóng)業(yè)綠色發(fā)展和生態(tài)環(huán)境保護均具有十分重要的意義。【前人研究進展】圍繞畜禽養(yǎng)殖污染現(xiàn)狀和環(huán)境負荷問題，學(xué)界已開展了很多工作。有研究表明畜禽糞便可在較高程度上填補氮、磷肥的需求，具有較大替代化肥潛力[6-7]，但不同省份間畜禽糞尿磷負荷量和還田量污染風險具有較大差異，多數(shù)地區(qū)的糞尿磷污染風險較大[8]。丹江口庫區(qū)作為我國南水北調(diào)中線工程源頭，是生態(tài)功能極重要區(qū)和生態(tài)環(huán)境極敏感區(qū)，水源地生態(tài)環(huán)境狀況直接關(guān)系到京、津、冀地區(qū)的戰(zhàn)略性用水安全。研究者[9]對水源地典型流域農(nóng)業(yè)面源污染的監(jiān)測顯示，水體中TN平均質(zhì)量濃度已超過Ⅳ類水標準，TP平均質(zhì)量濃度接近Ⅱ類水標準。綜合現(xiàn)有研究結(jié)論不難發(fā)現(xiàn)，畜禽養(yǎng)殖業(yè)伴隨的污染負荷對庫區(qū)負荷總量的貢獻超過50%[10-12]，是導(dǎo)致當前庫區(qū)水體污染的首要來源?！颈狙芯壳腥朦c】十堰市地處丹江口庫區(qū)下游腹地，是國家級生態(tài)示范區(qū)和鄂西北山區(qū)國家級生態(tài)功能保護區(qū)，同時還是南水北調(diào)中線工程的純調(diào)水區(qū)和核心水源涵養(yǎng)區(qū)，地理位置尤為特殊。且隨著中線工程正式通水以來配套“限養(yǎng)禁養(yǎng)”政策的推行，庫區(qū)的畜禽養(yǎng)殖業(yè)在發(fā)展結(jié)構(gòu)、空間分布等方面也正經(jīng)歷一系列變化，目前關(guān)于庫區(qū)畜禽養(yǎng)殖污染的認識主要集中在市域尺度或大流域尺度，抑或以整個畜禽養(yǎng)殖源為研究單元進行總體探討，而在縣域?qū)用驷槍π笄蒺B(yǎng)殖污染來源、主要污染物分布特征和畜禽養(yǎng)殖環(huán)境風險評價的研究仍顯不足，以致影響了庫區(qū)畜禽污染防控措施的制定和實施效果?！緮M解決的關(guān)鍵問題】基于上述背景，本研究依據(jù)十堰市畜禽養(yǎng)殖基本數(shù)據(jù)，對畜禽糞便總量和耕地負荷進行估算，分析各區(qū)縣養(yǎng)殖污染來源及主要污染物負荷特征；同時，參照十堰市主要農(nóng)作物種植面積，綜合運用畜禽養(yǎng)殖污染環(huán)境風險評價模型，評估各地區(qū)畜禽糞便耕地污染指數(shù)及風險等級；最后，從種養(yǎng)平衡的視角提出各地區(qū)養(yǎng)殖污染控制策略。研究結(jié)果可為十堰市畜禽養(yǎng)殖結(jié)構(gòu)調(diào)整和養(yǎng)殖污染減排提供基礎(chǔ)數(shù)據(jù)，特別是對于丹江口核心水源涵養(yǎng)區(qū)水環(huán)境保護以及農(nóng)業(yè)可持續(xù)發(fā)展提供有益借鑒。

1 材料與方法

1.1 研究區(qū)概況

十堰市范圍介于31°97′—33°25′N，109°51′—111°72′E。全市版圖面積2.36×104km2，轄三區(qū)四縣一市，分別是茅箭區(qū)、張灣區(qū)、鄖陽區(qū)、鄖西縣、竹溪縣、竹山縣、房縣和丹江口市。全市2020年總?cè)丝?46.16萬，其中農(nóng)業(yè)人口244.51萬。土壤以黃棕壤居多，地形地貌多樣，以山地和丘陵為主，海拔300— 1 700 m之間。屬北亞熱帶大陸性季風濕潤氣候，年平均氣溫13—16℃，年平均日照時數(shù)1 655—1 958 h，無霜期224—255 d。年平均降水量829 mm，徑流深263.4 mm。畜禽養(yǎng)殖業(yè)是十堰市農(nóng)業(yè)經(jīng)濟的重要組成部分，由于南水北調(diào)中線工程建設(shè)的需要，全市規(guī)?；B(yǎng)殖占比較低，畜禽養(yǎng)殖業(yè)以農(nóng)戶分散式經(jīng)營為主，大多都沒有建立糞污處理設(shè)施，污染物直接或間接排放的現(xiàn)象依然存在，給當?shù)剞r(nóng)村環(huán)境治理工作帶來了現(xiàn)實挑戰(zhàn)。

1.2 分析方法

1.2.1 污染負荷估算方法 本研究采用Johns輸出系數(shù)模型[13]對十堰市畜禽養(yǎng)殖業(yè)面源污染負荷總量進行估算。計算公式如下：

式中，T為污染物的負荷總量（kg）；E為第種畜禽養(yǎng)殖伴隨的第種污染物單頭或只產(chǎn)生量（kg·a-1），A為研究區(qū)中第種畜禽類型的養(yǎng)殖規(guī)模（頭或只）。

1.2.2 污染負荷評價方法 引用等標污染負荷法，在同一尺度上對不同污染源和污染物產(chǎn)生的環(huán)境負荷強度進行綜合評價，識別主要污染來源和主導(dǎo)污染因子。本研究主要考慮TN、TP和COD 3個污染評價因子。等標污染負荷為單位時間內(nèi)某污染物的產(chǎn)生量與環(huán)境限量標準的比值，計算公式如下[14-15]：

式中，P為第個污染源的第個污染物流失量等標污染負荷（m3）；M為第個污染源的第個污染物流失量（t）；C為基于水環(huán)境功能分區(qū)的第個污染物限量標準（mg·L-1），根據(jù)《丹江口庫區(qū)及上游水污染防治和水土保持“十三五”規(guī)劃》目標，丹江口庫區(qū)水質(zhì)管理標準采用《地表水環(huán)境質(zhì)量標準（GB3838—2002）》Ⅱ類水質(zhì)指標下限值，因此，在核算等標污染負荷時，TN、TP和COD的限量標準分別取0.5、0.1和15 mg·L-1；K為第個污染源的第個污染物等標污染負荷比。

1.2.3 畜禽養(yǎng)殖源環(huán)境風險評價方法 本研究引入畜禽糞便耕地負荷量和警報值，用以評價畜禽養(yǎng)殖糞便對環(huán)境造成的風險等級。具體數(shù)據(jù)來源和計算過程如下[16-18]：

=/=(∑XT)（4）

=（5）

式中，為以豬當量計的畜禽糞便耕地負荷量（t·hm-2）；為各類畜禽糞便豬糞當量（t）；為有效耕地面積（hm2）；X為各類畜禽糞便量（t）；T為以豬糞氮含量為基準的各類畜禽糞便豬當量轉(zhuǎn)換系數(shù)；為畜禽糞便耕地負荷警報值；為單位土地有機肥最大理論適宜施用量，由于耕地對畜禽糞便的承載力受土壤質(zhì)地、自然生態(tài)條件、社會經(jīng)濟條件和農(nóng)田管理水平等因素的綜合作用影響，因此不同地區(qū)間的糞肥施用限量值可能存在一定的差異。目前權(quán)威部門和學(xué)界關(guān)于糞肥土地最大承載量的描述大多介于30—45 t·hm-2之間，研究人員通常會根據(jù)研究區(qū)域的變化而確定不同的糞肥施用限值[7,19]。張緒美等[20]指出，應(yīng)以我國南北地理劃分線800 mm等降水量線（秦嶺—淮河—喜馬拉雅山東南端一線）為界，以南以45 t·hm-2為有機肥最大承載量，以北則以30 t·hm-2為最大承載量。根據(jù)這一論述，十堰市地處秦巴山余脈，年降雨量超過800 mm，應(yīng)采用的有機肥最大理論施用量為30 t·hm-2，且這一標準已在河南[21]、湖北[22]、四川[5,19]等地得到廣泛應(yīng)用。此外，從庫區(qū)環(huán)境敏感性、山區(qū)土壤肥力背景、坡耕地肥力固持能力和研究區(qū)復(fù)種指數(shù)等角度考量，本研究認為將糞肥最大理論施用量確定為30 t·hm-2是合理的。畜禽糞便負荷警報值分級標準參照現(xiàn)有方法[2,5,19,22]，詳見表1。

1.2.4 畜禽養(yǎng)殖環(huán)境容量的測度方法 為進一步明確基于環(huán)境負荷管理角度的養(yǎng)殖業(yè)總量控制水平，以豬當量表示畜禽養(yǎng)殖量，將不同畜禽排泄物N、P產(chǎn)生量換算為豬當量，以農(nóng)田N、P承載力計算畜禽養(yǎng)殖環(huán)境容量，公式如下[10,18]：

N/P=×N/P（6）

表1 畜禽糞便負荷警報值分級

=N/P/（7）

式（6）、（7）、（8）中，N/P為耕地N（P）環(huán)境容量（t），為有效耕地面積（hm2），N/P為歐盟地區(qū)單位土地糞肥年施N（P）限量標準：N=0.17 t·hm-2，P=0.035 t·hm-2[5,16-18]。為畜禽養(yǎng)殖環(huán)境容量（萬頭豬當量），為單頭豬糞便N（P）年產(chǎn)生量（t）。為畜禽養(yǎng)殖實際數(shù)量（萬頭豬當量），N(P)i為第種畜禽糞便N（P）年產(chǎn)生量（t）。

1.3 數(shù)據(jù)來源

本研究中涉及的畜禽養(yǎng)殖信息和耕地面積信息均來源于《十堰市統(tǒng)計年鑒-2020》。由于糧食作物、茶園、蔬菜、中藥材和果園的耕作面積約占全市耕地總面積的90%，因此本文將以上5種土地利用類型作為有效負載面積。養(yǎng)殖類型主要包括豬、牛、羊和家禽四類，根據(jù)對研究區(qū)實際養(yǎng)殖習(xí)慣的調(diào)查結(jié)果，確定以上4類畜禽的養(yǎng)殖周期分別為豬300 d、牛365 d、羊365 d、家禽365 d。畜禽糞便、羊養(yǎng)殖面源污染輸出系數(shù)根據(jù)《第一次全國污染源普查——畜禽養(yǎng)殖業(yè)源產(chǎn)排污系數(shù)手冊》[23]和相關(guān)資料[16-18]予以確定，豬、牛和家禽養(yǎng)殖面源污染輸出系數(shù)查自2021年發(fā)布的《排放源統(tǒng)計調(diào)查產(chǎn)排污核算方法和系數(shù)手冊》[24]。本文以氮含量為基準對糞肥產(chǎn)生量進行轉(zhuǎn)換。各類畜禽養(yǎng)殖規(guī)模、主要污染物年輸出系數(shù)和豬糞當量換算系數(shù)取值如表2所示。各縣（區(qū)、市）的畜禽養(yǎng)殖規(guī)模和耕地面積見表3。

表2 各類畜禽養(yǎng)殖量、污染物輸出系數(shù)及豬糞當量換算系數(shù)

表3 十堰市各縣（區(qū)、市）畜禽養(yǎng)殖結(jié)構(gòu)及耕地面積

2 結(jié)果

2.1 畜禽養(yǎng)殖污染物產(chǎn)生量及分布

根據(jù)畜禽飼養(yǎng)量、飼養(yǎng)期及污染物輸出系數(shù)，計算得到2020年十堰市各類畜禽糞便資源、TN、TP和COD總量及單位耕地負荷如表4所示。全市畜禽糞便豬糞當量產(chǎn)生總量為306.85×104t，其中牛類糞便產(chǎn)生量最多，為89.86×104t，占總量的29.28%，豬的糞便量其次，為77.95×104t，占總量的25.40%，羊的糞便量為71.71×104t，占總量的23.37%，家禽的糞便資源同樣不可忽視，達到67.33×104t，約占總量的21.94%。畜禽糞便中的化學(xué)污染物以COD為主，總量達到57.00×104t，且超過90%來源于牛和家禽。全市畜禽糞便TN產(chǎn)生量為2.63×104t，其中牛產(chǎn)生的TN量占全市的41.06%，家禽和羊產(chǎn)生的TN量大致相當，分別占全市的24.33%和23.95%。TP的產(chǎn)生量相對較少，僅為0.65×104t，且其來源結(jié)構(gòu)與TN類似。十堰市單位耕地畜禽糞便豬糞當量負荷強度為10.15 t·hm-2，與單位土地有機肥最大理論適宜施用量（30 t·hm-2）存在較大差距，TN和TP的耕地負荷同樣低于歐盟的限量標準（單位面積耕地TN、TP負荷分別為0.17和0.035 t·hm-2），分別為0.09和0.02 t·hm-2，此外，研究結(jié)果還表明，十堰地區(qū)單位耕地的畜禽糞便COD負荷強度為1.88 t·hm-2。

表4 十堰市畜禽養(yǎng)殖污染物負荷總量及單位耕地負荷強度

從區(qū)域分布來看（圖1—3），由于縣域間畜禽養(yǎng)殖數(shù)量和飼養(yǎng)結(jié)構(gòu)的不同，TN、TP和COD三類主要化學(xué)污染物的空間分布存在差異。TN和TP的區(qū)域分布特征基本一致，其中鄖陽區(qū)的污染量最大且主要來源于牛，TN和TP產(chǎn)生量分別為0.673×104和0.170×104t，占全市TN和TP總量的25.59%和26.15%。其次為丹江口市，其TN和TP產(chǎn)生量分別為0.516×104和0.123×104t，且主要來源于牛和家禽。房縣位居第3，TN和TP產(chǎn)生量分別為0.445×104和0.107× 104t，且主要來源于羊，以上3個地區(qū)產(chǎn)生的TN和TP之和超出全市總量的60%。鄖西縣、竹山縣、竹溪縣畜禽糞便的TN和TP產(chǎn)生量依次減少，三者貢獻了全市約36.62%的TN和36.92%的TP。畜禽糞便產(chǎn)生COD最多的是鄖陽區(qū)和丹江口市，且均來源于牛，其COD產(chǎn)生量分別為16.039×104和13.041×104t，二者占全市COD產(chǎn)生總量的51.02%，竹山縣、鄖西縣、竹溪縣和房縣的COD產(chǎn)生量大致相當，分別為7.547×104、7.068×104、6.746×104和6.050×104t，占全市COD產(chǎn)生總量的13.24%、12.40%、11.83%和10.61%。茅箭區(qū)和張灣區(qū)受限于養(yǎng)殖規(guī)模，其化學(xué)污染物的產(chǎn)生量相對較少，兩個地區(qū)TN、TP和COD的綜合產(chǎn)生量分別為0.035×104、0.007×104和0.507×104t，僅占全市總量的1.33%、1.08%和0.89%。

圖1 不同縣（區(qū)、市）畜禽養(yǎng)殖TN產(chǎn)生量

圖2 不同縣（區(qū)、市）畜禽養(yǎng)殖TP產(chǎn)生量

圖3 不同縣（區(qū)、市）畜禽養(yǎng)殖COD產(chǎn)生量

2.2 畜禽養(yǎng)殖等標污染負荷特征

分析十堰市畜禽養(yǎng)殖主要污染物及其來源貢獻，結(jié)果如表5所示。全市畜禽養(yǎng)殖等標污染負荷總量為15.532×1010m3，其中TP是畜禽養(yǎng)殖業(yè)的首要污染物，等標污染負荷為6.474×1010m3，占總負荷的41.68%；其次為TN（5.258×1010m3），占總負荷的33.85%；COD的等標污染負荷最小為3.800×1010m3，占總負荷的24.47%。分析畜禽養(yǎng)殖污染物的來源途徑可知，牛是畜禽養(yǎng)殖業(yè)的首要污染源，其等標污染負荷量為7.684×1010m3，占養(yǎng)殖業(yè)負荷總量的49.47%。與此同時，牛也是三類污染物的首要來源，對TN、TP和COD負荷總量的貢獻率分別達到40.93%、44.27%和70.15%。家禽是第二大畜禽養(yǎng)殖污染源，其等標污染負荷為3.354×1010m3，約占總負荷的21.60%。羊的養(yǎng)殖污染負荷為2.757×1010m3，等標污染負荷比為17.75%。豬養(yǎng)殖造成的污染等標負荷最小為1.737×1010m3，僅占總負荷的11.19%。

表5 畜禽養(yǎng)殖等標污染負荷總量及分布特征

P為第種污染物的等標污染負荷總量（×1010m3）；K為第種污染物總等標污染負荷比；K為第種污染源的第種污染物占第種污染物等標污染負荷總量的比值；P為第個污染源的等標污染負荷總量（×1010m3）；K為第個污染源的總等標污染負荷比

Pis the total equivalent pollution load of the pollutant(×1010m3);Kis the contribution of the pollutant;Kis the ratio of the pollutant i in pollution sourceto the total isostandard pollution load of the pollutant;Pis the total equivalent pollution load of the source of pollution(×1010m3);Kis the contribution of the source of pollution

進一步分析不同地區(qū)畜禽養(yǎng)殖業(yè)的等標污染負荷特征，結(jié)果如表6所示。十堰市畜禽養(yǎng)殖污染具有一定的空間分異性，其中鄖陽區(qū)是等標污染負荷最大的地區(qū)，負荷量達到4.110×1010m3，占全市負荷總量的26.46%。從負荷貢獻率來看，牛養(yǎng)殖伴隨的污染是導(dǎo)致該地區(qū)等標污染負荷較大的最主要原因。丹江口市的等標污染負荷僅次于鄖陽區(qū)為3.130×1010m3，占負荷總量的20.15%，其中超過86%的污染負荷來源于牛和家禽，房縣、鄖西縣、竹山縣和竹溪縣的負荷量依次減少，分別為2.366×1010、2.139×1010、1.947×1010和1.664×1010m3，占負荷總量的比例依次為15.24%、13.77%、12.54%和10.71%，其中羊?qū)Ψ靠h等標污染負荷的貢獻率最大，達到39.38%，是該地區(qū)的首要污染源，其他幾個地區(qū)的污染源依然集中于牛。茅箭區(qū)和張灣區(qū)的等標污染負荷均不超過0.1×1010m3，對全市的污染負荷貢獻較低且主要集中于家禽和羊。

2.3 畜禽養(yǎng)殖環(huán)境風險評價

以有效耕地面積為基準，從種養(yǎng)平衡的角度對十堰市畜禽養(yǎng)殖糞便耕地負荷及環(huán)境風險進行評價，結(jié)果如表7所示。畜禽糞便耕地負荷具有明顯的空間異質(zhì)性，其中茅箭區(qū)是十堰市畜禽糞便污染最為嚴重的地區(qū)，耕地負荷量達到22.071 t·hm-2，警報值為0.736，分級指數(shù)處于Ⅲ級，會對環(huán)境構(gòu)成威脅，說明盡管茅箭區(qū)的畜禽養(yǎng)殖量不大，但是它的耕地面積相對較少，無法消納畜禽養(yǎng)殖所產(chǎn)生的糞便量，因此它的耕地負荷較大，污染較為嚴重，需格外引起重視。鄖陽區(qū)、房縣和丹江口市的畜禽養(yǎng)殖耕地承載負荷警報值在0.4—0.7之間，屬于Ⅱ級，養(yǎng)殖環(huán)境風險狀態(tài)為“稍有”，也需加以關(guān)注。此外，鄖西縣、張灣區(qū)、竹山縣和竹溪縣的畜禽糞便耕地負荷依次降低，警報值均低于0.4，處于較為理想的狀態(tài)。十堰全市的畜禽糞便耕地負荷警報值為0.338，分級指數(shù)為Ⅰ，表明當前整個十堰市的畜禽糞便不會對環(huán)境構(gòu)成污染，縣域畜禽養(yǎng)殖環(huán)境風險存在差異是不同地區(qū)間養(yǎng)殖規(guī)模和耕地面積相互作用的結(jié)果。

從保持高產(chǎn)的角度出發(fā)，由于50%比例更趨近于實際管理情況，對控制畜禽養(yǎng)殖總量及合理調(diào)整養(yǎng)殖布局更具參考價值[16]，因此將實際養(yǎng)殖量與50%環(huán)境容量進行比較，得到養(yǎng)殖業(yè)環(huán)境污染風險指數(shù)，結(jié)果見表8。參考單位土地年施糞肥N、P限量標準，得出十堰市N、P耕地承載力分別為5.141×104和1.058×104t。以N為衡量標準，十堰市實際畜禽養(yǎng)殖量為625.95萬頭豬當量，超過50%環(huán)境容量（612.00萬頭豬當量）13.95萬頭豬當量，環(huán)境污染風險指數(shù)為1.023。以P為衡量標準，十堰市實際畜禽養(yǎng)殖量為539.51萬頭豬當量，超過50%環(huán)境容量（441.00萬頭豬當量）98.51萬頭豬當量，環(huán)境污染風險指數(shù)為1.223。由于P是十堰地區(qū)畜禽養(yǎng)殖業(yè)的首要污染物，其風險指數(shù)較N更高，且當前十堰市的畜禽養(yǎng)殖規(guī)模已超出環(huán)境安全承載水平，因此，從環(huán)境風險最小化的角度出發(fā)，應(yīng)以P的環(huán)境承載力和潛在風險為依據(jù)確定合理的養(yǎng)殖規(guī)模，即十堰市的養(yǎng)殖總量應(yīng)縮減98.51萬頭豬當量。具體分析各地區(qū)的環(huán)境風險可知，茅箭區(qū)和鄖陽區(qū)是十堰市畜禽養(yǎng)殖數(shù)量超載幅度最大的地區(qū)，其P污染風險指數(shù)都大于2，已經(jīng)遠遠超過理論環(huán)境容量，急需采取養(yǎng)殖總量強制性控制措施。丹江口市、房縣和鄖西縣的畜禽養(yǎng)殖量也處于過飽和狀態(tài)，其P污染風險指數(shù)均大于1，也會對環(huán)境構(gòu)成較大壓力，表明該地區(qū)的畜禽養(yǎng)殖總量管控和養(yǎng)殖污染消減工作也不容忽視。此外，張灣區(qū)的環(huán)境風險指數(shù)已趨近安全限值，也需引起重視。

表6 不同地區(qū)畜禽糞便等標污染負荷特征

P為地區(qū)的等標污染負荷總量（×1010m3）；K為地區(qū)的總等標污染負荷比；P為地區(qū)第個污染源的等標污染負荷量（×1010m3）；K為地區(qū)第個污染源的等標污染負荷比

Pis the total equivalent pollution load in region(×1010m3);Kis the contribution of region;Pis the equivalent pollution load of the source of pollutionin region(×1010m3);Kis the contribution of the source of pollutionin region

表7 不同地區(qū)畜禽糞便負荷警報值及環(huán)境風險

表8 不同地區(qū)畜禽養(yǎng)殖環(huán)境容量及實際總量

3 討論

3.1 畜禽養(yǎng)殖污染總量及其耕地負荷

本文依據(jù)2020年統(tǒng)計數(shù)據(jù)和國家權(quán)威部門推薦的畜禽糞便污染物輸出系數(shù)，估算十堰市畜禽糞便資源豬糞當量為306.85×104t，且4種畜禽類型的糞便總量差距相對較小，與黃美玲等[22]的測算結(jié)果總體相當。全市TN污染負荷量是TP污染負荷量的4.05倍，與現(xiàn)有研究中TN與TP的污染負荷量之比大多介于5—9的研究結(jié)論稍有出入[25]，這可能與不同研究區(qū)的養(yǎng)殖結(jié)構(gòu)存在差異有關(guān)。COD的產(chǎn)生量最大且明顯高于其他污染物，分別達到TN和TP負荷量的21.67倍和87.69倍。有研究指出湖北省畜禽養(yǎng)殖業(yè)COD總量是TN和TP總量的12.62倍和85.07倍[22]，從全國范圍來看，畜禽養(yǎng)殖業(yè)的COD總量約為TN和TP產(chǎn)生量的16.41倍和93.76倍[18]，可見，十堰市的畜禽養(yǎng)殖污染產(chǎn)生結(jié)構(gòu)與其他地區(qū)具有一致性。具體分析可知，養(yǎng)牛過程中的COD輸出系數(shù)較高和家禽養(yǎng)殖基數(shù)較大是促使十堰市COD總量顯著偏大的主要原因。全市污染物來源構(gòu)成具有明顯的空間分異性，其中鄖陽區(qū)、丹江口市和房縣是污染物產(chǎn)生量最為突出的地區(qū)，尤其值得關(guān)注。具體而言，以上3個地區(qū)貢獻了全市62.13%的TN、61.54%的TP和61.63%的COD，與之相對應(yīng)的是，以上3個地區(qū)容納了全市約58.33%的豬、60.45%的牛、61.94%的羊和66.82%的家禽養(yǎng)殖量，據(jù)此可以推斷，不同的畜禽養(yǎng)殖結(jié)構(gòu)和養(yǎng)殖規(guī)模是導(dǎo)致地域間污染總量存在顯著差異的決定性因素。

本文以糞便中TN含量為基準，對不同畜禽糞便的產(chǎn)生量進行轉(zhuǎn)換得到全市單位耕地畜禽糞便負荷強度為10.15 t·hm-2，與劉曉永等[8]以TP含量為依據(jù)提出的全國畜禽糞便單位耕地負荷（3.1 t·hm-2）存在較大差距，因此，在今后的研究工作中，很有必要對評價參數(shù)進行統(tǒng)一以便積累更多關(guān)聯(lián)性更強的參考資料。需要指出的是，在以TN含量為轉(zhuǎn)換因子的前提下，本文的糞便耕地負荷依然低于朱建春等[18]提出的湖北省平均水平（16.07 t·hm-2）和黃美玲等[22]關(guān)于湖北省的測算結(jié)果（16.2 t·hm-2），這可能與不同地區(qū)間的經(jīng)濟發(fā)展水平存在差距有關(guān)。十堰市地處鄂西北山區(qū)，畜牧養(yǎng)殖業(yè)和農(nóng)業(yè)水平較省內(nèi)其他地區(qū)相對滯后，導(dǎo)致該地區(qū)的畜禽糞便耕地負荷低于全省平均水平。本文的研究結(jié)果顯示，十堰市的單位耕地TN、TP負荷強度分別為0.09和0.02 t·hm-2，與現(xiàn)行的國際推薦標準尚有較大距離，整體負荷強度處于較低水平。通過比較其他研究者在三峽庫區(qū)[26]、海南省[7,27]、四川省[19]、湖北省[22]和全國[16,18]等不同尺度的研究結(jié)論可知，單位耕地的TN、TP負荷強度會隨著時間的推移和空間的轉(zhuǎn)換而發(fā)生較大幅度的變化，即使在同一研究區(qū)內(nèi)，不同地域間的負荷強度仍然可能存在明顯的差異。

3.2 畜禽養(yǎng)殖污染物及其來源

對十堰市畜禽養(yǎng)殖過程中主要污染物的解析結(jié)果表明，TP和TN是當?shù)匦笄蒺B(yǎng)殖環(huán)境管理中最需要關(guān)注的污染物，二者對全域等標污染負荷貢獻超過75%，COD的產(chǎn)生量最大，但其等標污染負荷率卻相對較小，這一結(jié)果與其他多數(shù)研究大致吻合。張英等[21]和郭利京等[28]對河南省和淮河生態(tài)經(jīng)濟帶的研究也揭示了畜禽養(yǎng)殖等標污染負荷貢獻最大的是TN，其次是TP，COD占比較小僅為6%左右。在丹江口庫區(qū)典型流域的研究也證實，畜禽養(yǎng)殖等標污染物以TN和TP為主，其總負荷率約為95%[10]。然而，時曉瑞等[29]對淮河流域奎濉河典型小流域的研究則指出，COD是畜禽養(yǎng)殖最重要的污染物，負荷貢獻率達到69.91%，其次是TN的23.71%，TP的等標負荷貢獻僅為6.37%，這可能與該地區(qū)養(yǎng)殖業(yè)發(fā)展過程中以產(chǎn)生COD為主的家禽養(yǎng)殖規(guī)模明顯高于其他畜禽種類密切相關(guān)。分析污染物的來源可知，牛和家禽的綜合污染量超過了十堰市畜禽養(yǎng)殖面源污染負荷總量的70%，是今后開展面源污染消減工作時需要特別關(guān)注的畜禽種類。與此同時，從整個漢江流域[30]來看，牛和豬是最主要的畜禽污染來源，其等標污染負荷率分別為35.95%和34.98%，家禽的負荷率緊隨其后為26.28%，羊的負荷率較小，僅為2.79%。而對丹江口庫區(qū)典型農(nóng)業(yè)小流域面源污染源解析的結(jié)果則證實，生豬TN、TP和COD的等標負荷占畜禽養(yǎng)殖源負荷總量的比例分別高達99.29%、92.73%和93.94%，是流域內(nèi)畜禽養(yǎng)殖污染乃至整體農(nóng)業(yè)面源污染防控工作的重點[10]，與本文的研究結(jié)果存在較大差異，這主要是因為該小流域內(nèi)建有規(guī)?；i養(yǎng)殖企業(yè)，其生豬存欄基數(shù)大幅超出其他畜禽養(yǎng)殖量，使得生豬成為流域內(nèi)的主要污染來源?？梢姡笄蒺B(yǎng)殖的主要污染物和主要污染來源與特定地域的產(chǎn)業(yè)布局息息相關(guān)，在制定污染管控策略時，只有結(jié)合不同地區(qū)的養(yǎng)殖污染特征才能精準施策，以提升污染治理的針對性和時效性。十堰市各地區(qū)間污染負荷貢獻率和污染結(jié)構(gòu)存在較大差異，一方面與區(qū)域內(nèi)畜禽養(yǎng)殖總量高度相關(guān)，顯而易見的，茅箭區(qū)和張灣區(qū)因為持有的養(yǎng)殖總量大幅低于其他地區(qū)，所以其養(yǎng)殖污染負荷率占比較小。另一方面，也受高污染畜禽養(yǎng)殖量和畜禽養(yǎng)殖結(jié)構(gòu)的共同作用。牛作為污染物輸出系數(shù)最高的畜禽種類，在鄖陽區(qū)和丹江口市的分布比例明顯高于其他地區(qū)，且均貢獻了以上兩個地區(qū)超出一半的負荷壓力，是導(dǎo)致鄖陽區(qū)和丹江口市成為全市污染負載量最大地區(qū)的重要原因。竹溪縣、竹山縣和鄖西縣的牛養(yǎng)殖規(guī)模相對接近且依次低于以上兩個地區(qū)，三者的主要污染也均來源于牛。而作為十堰市畜禽養(yǎng)殖污染第三高的地區(qū)，房縣的主要污染來源于羊，其次為牛，這是因為房縣容納了全市最大體量的羊，而牛的養(yǎng)殖規(guī)模僅列第六，固然羊的污染系數(shù)低于牛，但養(yǎng)殖規(guī)模上的巨大差距依然使得羊的污染負荷占據(jù)主動。說明地區(qū)間的污染負荷特征差異不僅取決于高污染畜禽養(yǎng)殖量，也可能受特定地區(qū)內(nèi)的養(yǎng)殖結(jié)構(gòu)所左右。此外，不同地區(qū)的養(yǎng)殖面源污染特征還受生態(tài)環(huán)境、社會經(jīng)濟水平、種養(yǎng)配套能力以及養(yǎng)殖區(qū)和水源地之間的距離差異等因素影響，需要在今后的工作中進一步深入研究。

3.3 畜禽養(yǎng)殖環(huán)境風險評價

畜禽養(yǎng)殖環(huán)境污染風險的評價結(jié)果顯示，當前十堰市畜禽養(yǎng)殖糞便帶來的環(huán)境風險尚處于較為理想的分級區(qū)間之內(nèi)，但縣域之間存在一定的空間異質(zhì)性。全市有1/2以上的地區(qū)存在畜禽養(yǎng)殖超載現(xiàn)象，其中尤以茅箭區(qū)和鄖陽區(qū)的超載幅度最為突出，其P污染風險指數(shù)分別達到2.429和2.133，說明其養(yǎng)殖規(guī)模已達環(huán)境安全上限的兩倍以上，處于嚴重超載區(qū)間，對養(yǎng)殖規(guī)模采取總量減半的措施迫在眉睫。進一步分析，由于茅箭區(qū)的耕地面積較小，環(huán)境承載量相對有限，較低水平的過飽和養(yǎng)殖量就會導(dǎo)致風險指數(shù)偏高，因此其養(yǎng)殖規(guī)模絕對消減量并不突出，而鄖陽區(qū)在耕地面積相對較大的情況下，其環(huán)境風險指數(shù)依然超過2，說明該地區(qū)的超載現(xiàn)象已十分嚴重，是十堰地區(qū)養(yǎng)殖規(guī)模管控的重中之重。丹江口市、房縣和鄖西縣的環(huán)境風險指數(shù)也已超過安全臨界值，說明這些地區(qū)的畜禽養(yǎng)殖也會對環(huán)境造成不同程度的污染，考慮到這些地區(qū)的畜禽養(yǎng)殖基數(shù)較大，且相對充足的有效耕地面積理論上能承載更大的養(yǎng)殖量，說明這些地區(qū)超載的養(yǎng)殖規(guī)模絕對值也已處于很高水平，從全市層面看，十堰市的畜禽糞便環(huán)境負荷警報值為0.338，N、P環(huán)境風險指數(shù)分別為1.023和1.223，基于P平衡視角的十堰市畜禽養(yǎng)殖規(guī)模須在當前水平上削減98.51萬頭豬當量。本文的結(jié)論與前人的研究不盡一致，有研究得到2010年湖北省畜禽養(yǎng)殖的N、P環(huán)境風險指數(shù)分別為1.51和1.12[16]，2011年分別為1.7和1.4[22]，其中2011年十堰市的糞便負荷警報值為0.52，N、P環(huán)境風險指數(shù)分別為1.5和1.3[22]。湖北省三峽庫區(qū)的N、P風險指數(shù)分別達到3.3和2.5[26]?？梢?，十堰市的畜禽養(yǎng)殖環(huán)境風險較全國和省域內(nèi)其他地區(qū)整體較低，但不同學(xué)者對不同時段的研究結(jié)論也存在差異。本文關(guān)于十堰市畜禽養(yǎng)殖環(huán)境負荷和環(huán)境風險的描述較其他研究者提出的結(jié)果更低，這主要與近年來十堰市養(yǎng)殖規(guī)模的大幅縮減有關(guān)，一方面由于庫區(qū)蓄水需要，2009—2011年，十堰市部分農(nóng)戶移民外地，養(yǎng)殖主體減少促使養(yǎng)殖量下降。另一方面，更重要的是自2014年12月南水北調(diào)中線工程正式通水以來，丹江口庫區(qū)沿線關(guān)于畜禽養(yǎng)殖業(yè)的發(fā)展規(guī)劃和限養(yǎng)政策進一步收緊，致使庫區(qū)畜禽養(yǎng)殖規(guī)模大幅縮減。查閱資料可知，2012—2020年，十堰市的生豬養(yǎng)殖量降幅達66.77%，牛的養(yǎng)殖量減少了41.76%，羊的養(yǎng)殖規(guī)模降低了17.75%，家禽養(yǎng)殖量縮小了34.57%。相對于大幅縮減的畜禽養(yǎng)殖規(guī)模，同期的全市有效耕地面積降幅相對較小，約為8.64%。更小的養(yǎng)殖體量意味著更少的畜禽糞便排放，而相對較小的耕地面積降幅則意味著具有相對穩(wěn)定的環(huán)境承載能力，正是以上兩方面原因促使十堰市的養(yǎng)殖環(huán)境風險較以往出現(xiàn)下降。

需要指出的是，單位耕地畜禽糞便年養(yǎng)分施用限量標準往往需要根據(jù)特定地域的土壤質(zhì)地、坡度、降雨、氣候等因素來確定，基于此，歐盟將糞肥年施用養(yǎng)分量分別確定為0.17 t·hm-2和0.035 t·hm-2，超過這個限值即認為會造成硝酸鹽和磷污染[31-32]。由于國內(nèi)外在種植制度、作物類型、生產(chǎn)方式、土壤質(zhì)地、農(nóng)業(yè)氣候等方面存在明顯差異，因此直接引用國外標準可能會引起一定的誤差。對此，部分國內(nèi)學(xué)者針對北京平谷區(qū)[33]、上海[34]、河北[35]、海南[7]等地區(qū)的耕地或果園等做了相關(guān)研究，并根據(jù)研究區(qū)域的變化提出了不同的糞肥施用限量標準?？梢?，我國學(xué)者對不同地區(qū)和不同土地利用方式下的養(yǎng)分合理需求量也尚未形成共識。本研究在環(huán)境容量計算時采用歐盟地區(qū)的單位土地糞肥年施氮磷限量標準，主要出于以下幾點考慮。其一，這一歐盟地區(qū)限值在我國被廣泛應(yīng)用[5,10,16,18,22]，具有比較成熟的參考價值；其二，十堰市糧食作物面積占比較低，有效耕地中以茶葉、林果、中藥材等多年生植被覆蓋為主，整體復(fù)種指數(shù)小，對養(yǎng)分的消耗和需求較低；其三，十堰山區(qū)土地破碎，農(nóng)田以坡耕地為主，水土流失和徑流強度相對較大，過多的養(yǎng)分投入可能會造成大量氮磷流失，從環(huán)境風險管理的角度出發(fā)，需采取更加安全的養(yǎng)分投入水平。因此，本文認為采取歐盟的氮磷限量標準是比較合理的，縱然如此，如何根據(jù)庫區(qū)農(nóng)業(yè)特點確定更加準確的氮磷投入安全閾值依然是今后的一項重要工作。

此外，本文所述的耕地負荷警報值和環(huán)境風險指數(shù)是在假設(shè)畜禽糞便全量還田的情況下得出的，然而在實際情況下，一方面由于糞尿的就地揮發(fā)和滲漏、糞便制成商品有機肥后跨境流通、零散糞污的隨意拋棄以及運輸損失等因素的影響，畜禽糞污不可能全部被應(yīng)用于當?shù)馗?。另一方面，是否所有的耕地均具有滿足畜禽糞便還田的現(xiàn)實條件依然需要加以考證。因此，由公式推導(dǎo)出來的理論值可能無法足夠準確地反映出畜禽糞便的實際環(huán)境風險，這也是今后相關(guān)研究工作中無法回避的科學(xué)問題。但從丹江口庫區(qū)這一特殊地域水環(huán)境保護的角度考量，庫區(qū)的環(huán)境污染風險評價工作應(yīng)趨于保守，且本文計算的十堰市污染物產(chǎn)生總量是客觀存在的，文中列出的耕地面積也是在剔除撂荒地的基礎(chǔ)上考慮了復(fù)種指數(shù)之后的實際耕種面積，具備較高的還田利用行為可及性。因此，本文關(guān)于畜禽糞便潛在環(huán)境風險的評價結(jié)論依然能為有關(guān)部門做出謹慎的宏觀決策提供重要參考。

3.4 畜禽養(yǎng)殖環(huán)境風險防控策略

為進一步管控和降低水源區(qū)畜禽養(yǎng)殖的環(huán)境風險，筆者認為應(yīng)從以下幾個方面考慮采取操作性較強的防控措施。一是要根據(jù)環(huán)境容量適當調(diào)整各地區(qū)畜禽養(yǎng)殖規(guī)模，建議考慮生產(chǎn)布局并加強區(qū)域間協(xié)同發(fā)展[36]，具體來看，可以將鄖陽區(qū)、丹江口市、房縣、鄖西縣和茅箭區(qū)的養(yǎng)殖規(guī)?？s減75.065、41.941、35.385、10.509和2.136萬頭豬當量，同時要鼓勵環(huán)境污染低風險區(qū)適當擴大畜禽養(yǎng)殖規(guī)模，將竹溪縣和竹山縣的養(yǎng)殖量擴增空間分別設(shè)定為35.558和30.718萬頭豬當量，以達到在環(huán)境安全閾值內(nèi)保持全域畜禽養(yǎng)殖總量的相對穩(wěn)定。二是要提升養(yǎng)殖廢棄物循環(huán)利用率。本研究是在50%畜禽糞便還田率的基礎(chǔ)上計算得到十堰市耕地糞便承載力和養(yǎng)殖規(guī)模上限的，為達到這一還田比例，應(yīng)重點從以下三個方面考慮：一方面要綜合作物—畜禽系統(tǒng)的發(fā)展[37]，提高種養(yǎng)協(xié)同性和可持續(xù)性。例如，可以在竹溪縣和竹山縣等茶葉和林果主產(chǎn)區(qū)，大力發(fā)展“茶園養(yǎng)雞”“林下養(yǎng)雞”模式，將其他地區(qū)的家禽養(yǎng)殖份額轉(zhuǎn)移至這些地區(qū)；在張灣區(qū)和茅箭區(qū)等蔬菜主產(chǎn)區(qū)，推廣應(yīng)用豬場廢水[38]回田技術(shù)等。另一方面，要大力推行養(yǎng)分替代利用技術(shù)。例如可以在鄖陽區(qū)、鄖西縣等玉米-小麥主產(chǎn)區(qū)推廣有機肥替代氮肥利用技術(shù)[39]，在規(guī)模化生豬養(yǎng)殖區(qū)周邊構(gòu)建“豬-沼-茶/果/菌”生態(tài)循環(huán)模式等。此外，要更加注重政策引導(dǎo)。可以通過強化農(nóng)業(yè)綠色發(fā)展理念，制定有機肥回田利用以獎代補政策等方式，鼓勵種養(yǎng)主體深度融合，實現(xiàn)養(yǎng)殖廢棄物就地就近消納。三是要注重加強對糞便的處理和資源化利用。尤其要加強對規(guī)?；B(yǎng)殖的雨污分流和固液分離處理基礎(chǔ)設(shè)施建設(shè)，減少糞污無序排放。在資源化利用方面，厭氧發(fā)酵處理是一種可同時實現(xiàn)畜禽廢棄物減量化、資源化和能源化的高效資源化技術(shù)[6]，應(yīng)引起重視。堆肥手段具有對糞污無害化處理比較徹底、糞便附加值高、經(jīng)濟效益好等優(yōu)點，是目前較為主流的畜禽糞便處理模式[40]。此外，畜禽糞便的能源化和飼料化[41]利用潛力巨大，也是可行的生態(tài)循環(huán)處理方式。四是要進一步完善政府部門對畜禽養(yǎng)殖風險的跟蹤管理[42-43]。加強政府監(jiān)管是倒逼畜禽養(yǎng)殖廢棄物合規(guī)處理的主要手段，也是推動畜禽養(yǎng)殖污染防治與種養(yǎng)結(jié)合規(guī)劃編制工作的現(xiàn)實需要?？傊笄蒺B(yǎng)殖環(huán)境風險防控是一項系統(tǒng)工程，受特定地域的資源條件、產(chǎn)業(yè)結(jié)構(gòu)、利用現(xiàn)狀、糞便處理方式等多種因素的綜合影響。在制定防控策略時，應(yīng)根據(jù)不同地區(qū)養(yǎng)殖污染的來源結(jié)構(gòu)及風險特征，因地制宜地從技術(shù)、觀念、政策等角度入手進行分類指導(dǎo)，針對性地提出強制管控策略和柔性引導(dǎo)方式，以實現(xiàn)種養(yǎng)平衡、綠色可持續(xù)發(fā)展的目標。

4 結(jié)論

4.1 十堰市畜禽糞便豬糞當量總計為306.85×104t，糞便中主要污染物產(chǎn)生量表現(xiàn)為COD＞TN＞TP，且主要來源于鄖陽區(qū)、丹江口市和房縣3個地區(qū)。

4.2 十堰市畜禽養(yǎng)殖等標污染負荷總量為15.532×1010m3，且不同地區(qū)間的污染負荷結(jié)構(gòu)存在明顯差異。養(yǎng)殖污染等標負荷率表現(xiàn)為TP＞TN＞COD，養(yǎng)牛是庫區(qū)首要污染源。

4.3 十堰市畜禽糞便耕地負荷整體較低，不會對環(huán)境構(gòu)成壓力，但存在明顯的空間分異性。總體來看，茅箭區(qū)和鄖陽區(qū)是養(yǎng)殖污染風險較高的地區(qū)，其次為丹江口市和房縣，鄖西縣的養(yǎng)殖污染也需引起重視。

4.4 以P為衡量標準，十堰市的大部分地區(qū)養(yǎng)殖規(guī)模已經(jīng)超載，需嚴格控制畜禽養(yǎng)殖總量，同時做好廢棄物的循環(huán)利用。

[1] 劉曉永, 李書田. 中國畜禽糞尿養(yǎng)分資源及其還田的時空分布特征. 農(nóng)業(yè)工程學(xué)報, 2018, 34(4): 1-14, 316.

LIU X Y, LI S T. Temporal and spatial distribution of nutrient resource from livestock and poultry feces and its returning to cropland. Transactions of the Chinese Society of Agricultural Engineering, 2018, 34(4): 1-14, 316. (in Chinese)

[2] 國家環(huán)境保護總局自然生態(tài)保護司. 全國規(guī)?；笄蒺B(yǎng)殖業(yè)污染情況調(diào)查及防治對策. 北京: 中國環(huán)境科學(xué)出版社, 2002.

Department of Natural Ecology Protection, State Environmental Protection Administration. Investigation and Countermeasures on Pollution of Large-Scale Livestock and Poultry Industry in China. Beijing: China Environment Science Press, 2002. (in Chinese)

[3] 高定, 陳同斌, 劉斌, 鄭袁明, 鄭國砥, 李艷霞. 我國畜禽養(yǎng)殖業(yè)糞便污染風險與控制策略. 地理研究, 2006, 25(2): 311-319.

GAO D, CHEN T B, LIU B, ZHENG Y M, ZHENG G D, LI Y X. Releases of pollutants from poultry manure in China and recommended strategies for the pollution prevention. Geographical Research, 2006, 25(2): 311-319. (in Chinese)

[4] 胡靜鋒. 重慶市農(nóng)業(yè)面源污染測算與空間特征解析. 中國農(nóng)業(yè)資源與區(qū)劃, 2017, 38(1): 135-144.

HU J F. Measurement and spatial feature analysis of agricultural non-point source pollution in Chongqing. Chinese Journal of Agricultural Resources and Regional Planning, 2017, 38(1): 135-144. (in Chinese)

[5] 吳永勝, 孫越鴻, 楊雪, 許禎瑩, 齊桂蘭, 邱時秀, 李娟. 基于種養(yǎng)平衡的成都市畜禽養(yǎng)殖環(huán)境效應(yīng)分析. 中國農(nóng)業(yè)資源與區(qū)劃, 2018, 39(1): 195-203.

WU Y S, SUN Y H, YANG X, XU Z Y, QI G L, QIU S X, LI J. Analysis of environmental effect of livestock and poultry breeding based on planting-breeding balance in Chengdu. Chinese Journal of Agricultural Resources and Regional Planning, 2018, 39(1): 195-203. (in Chinese)

[6] 吳浩瑋, 孫小淇, 梁博文, 陳家斌, 周雪飛. 我國畜禽糞便污染現(xiàn)狀及處理與資源化利用分析. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 2020, 39(6): 1168-1176.

WU H W, SUN X Q, LIANG B W, CHEN J B, ZHOU X F. Analysis of livestock and poultry manure pollution in China and its treatment and resource utilization. Journal of Agro-Environment Science, 2020, 39(6): 1168-1176. (in Chinese)

[7] 楊旭, 黃艷艷, 劉海林, 茶正早, 羅微, 林清火. 海南省畜禽養(yǎng)殖環(huán)境承載力及有機肥替代化肥潛力分析. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 2019, 38(11): 2609-2618.

YANG X, HUANG Y Y, LIU H L, CHA Z Z, LUO W, LIN Q H. Potential for the environmental carrying capacity of livestock and poultry breeding to substitute organic fertil? izers in Hainan Province, China. Journal of Agro-Environment Science, 2019, 38(11): 2609- 2618. (in Chinese)

[8] 劉曉永, 王秀斌, 李書田. 中國農(nóng)田畜禽糞尿磷負荷量及環(huán)境風險分析. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 2019, 38(11): 2594-2608.

LIU X Y, WANG X B, LI S T. Phosphorus loading rates from livestock and poultry faeces, and environmental evaluation in China. Journal of Agro-Environment Science, 2019, 38(11): 2594-2608. (in Chinese)

[9] 龔世飛, 丁武漢, 肖能武, 郭元平, 葉青松, 王巍, 李虎. 丹江口水庫核心水源區(qū)典型流域農(nóng)業(yè)面源污染特征. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 2019, 38(12): 2816-2825.

GONG S F, DING W H, XIAO N W, GUO Y P, YE Q S, WANG W, LI H. Characteristics of surface runoff and agricultural non-point source pollution in the core water source area of the Danjiangkou Reservoir. Journal of Agro-Environment Science, 2019, 38(12): 2816-2825. (in Chinese)

[10] 龔世飛, 丁武漢, 居學(xué)海, 肖能武, 葉青松, 黃進, 李虎. 典型農(nóng)業(yè)小流域面源污染源解析與控制策略: 以丹江口水源涵養(yǎng)區(qū)為例. 中國農(nóng)業(yè)科學(xué), 2021, 54(18): 3919-3931.

GONG S F, DING W H, JU X H, XIAO N W, YE Q S, HUANG J, LI H. Source analysis and control strategies of non-point source pollution in typical agricultural small watershed: a case study of Danjiangkou water conservation area. Scientia Agricultura Sinica, 2021, 54(18): 3919-3931. (in Chinese)

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

JIANG S Y, HAN P, JIA Z B, MAO X L, SI J, YE F F. Evaluation with PSR model and GIS analysis of agricultural non-point source pollution in Danjiangkou reservoir of the mid-route of the south-to-north water transfer project. Journal of Agro-Environment Science, 2010, 29(11): 2153-2162. (in Chinese)

[12] 王國重, 李中原, 左其亭, 屈建鋼, 李曉宇. 丹江口水庫水源區(qū)農(nóng)業(yè)面源污染物流失量估算. 環(huán)境科學(xué)研究, 2017, 30(3): 415-422.

WANG G Z, LI Z Y, ZUO Q T, QU J G, LI X Y. Estimation of agricultural non-point source pollutant loss in catchment areas of Danjiangkou reservoir. Research of Environmental Sciences, 2017, 30(3): 415-422. (in Chinese)

[13] 劉亞瓊, 楊玉林, 李法虎. 基于輸出系數(shù)模型的北京地區(qū)農(nóng)業(yè)面源污染負荷估算. 農(nóng)業(yè)工程學(xué)報, 2011, 27(7): 7-12.

LIU Y Q, YANG Y L, LI F H. Estimation of pollution loads from agricultural nonpoint sources in Beijing region based on export coefficient modeling approach. Transactions of the Chinese Society of Agricultural Engineering, 2011, 27(7): 7-12. (in Chinese)

[14] 許策, 李曄, 束繼年, 李松炳, 李明高. 漢江流域荊門段面源污染負荷時空分布與污染現(xiàn)狀評價. 水土保持通報, 2017, 37(4): 63-68.

XU C, LI Y, SHU J N, LI S B, LI M G. Spatial distribution and evaluation of non-point pollution source in Jingmen region of Hanjiang watershed. Bulletin of Soil and Water Conservation, 2017, 37(4): 63-68. (in Chinese)

[15] 謝經(jīng)朝, 趙秀蘭, 何丙輝, 李章安. 漢豐湖流域農(nóng)業(yè)面源污染氮磷排放特征分析. 環(huán)境科學(xué), 2019, 40(4): 1760-1769.

XIE J C, ZHAO X L, HE B H, LI Z G. Analysis of the characteristics of nitrogen and phosphorus emissions from agricultural non-point sources on Hanfeng Lake Basin. Environmental Science, 2019, 40(4): 1760-1769. (in Chinese)

[16] 耿維, 胡林, 崔建宇, 卜美東, 張蓓蓓. 中國區(qū)域畜禽糞便能源潛力及總量控制研究. 農(nóng)業(yè)工程學(xué)報, 2013, 29(1): 171-179, 295.

GENG W, HU L, CUI J Y, BU M D, ZHANG B B. Biogas energy potential for livestock manure and gross control of animal feeding in region level of China. Transactions of the Chinese Society of Agricultural Engineering, 2013, 29(1): 171-179, 295. (in Chinese)

[17] 王方浩, 馬文奇, 竇爭霞, 馬林, 劉小利, 許俊香, 張福鎖. 中國畜禽糞便產(chǎn)生量估算及環(huán)境效應(yīng). 中國環(huán)境科學(xué), 2006, 26(5): 614-617.

WANG F H, MA W Q, DOU Z X, MA L, LIU X L, XU J X, ZHANG F S. The estimation of the production amount of animal manure and its environmental effect in China. China Environmental Science, 2006, 26(5): 614-617. (in Chinese)

[18] 朱建春, 張增強, 樊志民, 李榮華. 中國畜禽糞便的能源潛力與氮磷耕地負荷及總量控制. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 2014, 33(3): 435-445.

ZHU J C, ZHANG Z Q, FAN Z M, LI R H. Biogas potential, cropland load and total amount control of animal manure in China. Journal of Agro-Environment Science, 2014, 33(3): 435-445. (in Chinese)

[19] 郭珊珊, 張涵, 楊汝馨. 基于耕地承載力的畜禽養(yǎng)殖污染負荷及環(huán)境風險研究: 以四川省為例. 水土保持通報, 2019, 39(1): 226-232, 325.

GUO S S, ZHANG H, YANG R X. Pollutant loads and environmental risks of livestock production based on arable land carrying capacity—A case study in Sichuan Province. Bulletin of Soil and Water Conservation, 2019, 39(1): 226-232, 325. (in Chinese)

[20] 張緒美, 董元華, 王輝, 沈旦. 中國畜禽養(yǎng)殖結(jié)構(gòu)及其糞便N污染負荷特征分析. 環(huán)境科學(xué), 2007, 28(6): 1311-1318.

ZHANG X M, DONG Y H, WANG H, SHEN D. Structure of livestock and variation of fecal nitrogen pollution load in China. Environmental Science, 2007, 28(6): 1311-1318. (in Chinese)

[21] 張英, 武淑霞, 劉宏斌, 雷秋良, 王洪媛, 翟麗梅. 基于種養(yǎng)平衡的河南省畜禽養(yǎng)殖分析及其環(huán)境污染風險研究. 中國土壤與肥料, 2019(4): 24-30, 52.

ZHANG Y, WU S X, LIU H B, LEI Q L, WANG H Y, ZHAI L M. Analysis on rationality and pollution of livestock industry in Henan Province based on balance between crop-planting and livestock- raising. Soils and Fertilizers Sciences in China, 2019(4): 24-30, 52. (in Chinese)

[22] 黃美玲, 夏穎, 范先鵬, 黃敏, 吳茂前, 劉冬碧, 張富林. 湖北省畜禽養(yǎng)殖污染現(xiàn)狀及總量控制. 長江流域資源與環(huán)境, 2017, 26(2): 209-219.

HUANG M L, XIA Y, FAN X P, HUANG M, WU M Q, LIU D B, ZHANG F L. Pollution status and total amount control of livestock and poultry breeding in Hubei Province. Resources and Environment in the Yangtze Basin, 2017, 26(2): 209-219. (in Chinese)

[23] 中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)環(huán)境與可持續(xù)發(fā)展研究所, 環(huán)境保護部南京環(huán)境科學(xué)研究所. 第一次全國污染源普查畜禽養(yǎng)殖業(yè)源產(chǎn)排污系數(shù), 2009.

Institute of Environment and Sustainable Development in Agriculture CAAS, Nanjing Institute of Environmental Sciences, MEP. The first national survey of pollution sources on pollution discharge coefficients of livestock and poultry industry, 2009. (in Chinese)

[24] 中華人民共和國生態(tài)環(huán)境部. 排放源統(tǒng)計調(diào)查產(chǎn)排污核算方法和系數(shù)手冊. [2021-06-11]. https://www.mee.gov.cn/xxgk2018/xxgk/ xxgk01/202106/t20210618_839512.html

Ministry of Ecology and Environment of the People’s Republic of China. Emission source statistics survey production and emission accounting method and coefficient manual. [2021-06-11]. https:// www.mee.gov.cn/xxgk2018/xxgk/xxgk01/202106/t20210618_839512. Html. (in Chinese)

[25] 盧少勇, 張萍, 潘成榮, 彭書傳, 劉曉暉. 洞庭湖農(nóng)業(yè)面源污染排放特征及控制對策研究. 中國環(huán)境科學(xué), 2017, 37(6): 2278-2286.

LU S Y, ZHANG P, PAN C R, PENG S C, LIU X H. Agricultural non-point source pollution discharge characteristic and its control measures of Dongtinghu Lake. China Environmental Science, 2017, 37(6): 2278-2286. (in Chinese)

[26] 黃美玲, 夏穎, 范先鵬, 黃敏, 吳茂前, 劉冬碧, 張富林. 湖北省三峽庫區(qū)畜禽養(yǎng)殖污染現(xiàn)狀評價. 生態(tài)學(xué)雜志, 2017, 36(3): 725-733.

HUANG M L, XIA Y, FAN X P, HUANG M, WU M Q, LIU D B, ZHANG F L. Evaluation of livestock pollution status in the Three Gorges Reservoir Area of Hubei Province. Chinese Journal of Ecology, 2017, 36(3): 725-733. (in Chinese)

[27] 丁尚, 付陽, 郭浩浩, 宋晨陽, 李博玲, 趙洪偉. 海南島1988—2018年畜禽糞尿氮磷負荷量及環(huán)境效應(yīng). 中國農(nóng)業(yè)科學(xué), 2020, 53(18): 3752-3763.

DING S, FU Y, GUO H H, SONG C Y, LI B L, ZHAO H W. Nitrogen and phosphorus loads and their environmental effects of animal manure in Hainan Island from 1988 to 2018. Scientia Agricultura Sinica, 2020, 53(18): 3752-3763. (in Chinese)

[28] 郭利京, 黃振英. 淮河生態(tài)經(jīng)濟帶農(nóng)業(yè)面源污染空間分布及治理研究. 長江流域資源與環(huán)境, 2021, 30(7): 1746-1756.

GUO L J, HUANG Z Y. Study on spatial distribution and control of agricultural non-point source pollution in Huaihe ecological economic belt. Resources and Environment in the Yangtze Basin, 2021, 30(7): 1746-1756. (in Chinese)

[29] 時曉瑞, 武升, 黃瑜, 歐三青, 王樹仁, 馬友華, 姜家生. 淮河流域奎濉河典型小流域農(nóng)業(yè)面源污染調(diào)查與評價. 生態(tài)科學(xué), 2021, 40(5): 109-115.

SHI X R, WU S, HUANG Y, OU S Q, WANG S R, MA Y H, JIANG J S. Investigation and evaluation of agricultural non-point source pollution in typical small watershed of Kuisui River in Huaihe Basin in Suzhou City, Anhui Province. Ecological Science, 2021, 40(5): 109-115. (in Chinese)

[30] 唐肖陽, 唐德善, 魯佳慧, 常文倩, 唐新玥. 漢江流域農(nóng)業(yè)面源污染的源解析. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 2018, 37(10): 2242-2251.

TANG X Y, TANG D S, LU J H, CHANG W Q, TANG X Y. Source apportionment of agricultural nonpoint source pollution in the Hanjiang River Basin. Journal of Agro-Environment Science, 2018, 37(10): 2242-2251. (in Chinese)

[31] HENKENS P L C M, KEULEN H V. Mineral policy in the Netherlands and nitrate policy within the European Community. NJAS - Wageningen Journal of Life Sciences, 2001, 49(2/3): 117-134.

[32] OENEMA O, VAN LIERE L, PLETTE S, PRINS T, VAN ZEIJTS H, SCHOUMANS O. Environmental effects of manure policy options in The Netherlands. Water Science and Technology, 2004, 49(3): 101-108.

[33] 潘瑜春, 孫超, 劉玉, 唐秀美, 任艷敏. 基于土地消納糞便能力的畜禽養(yǎng)殖承載力. 農(nóng)業(yè)工程學(xué)報, 2015, 31(4): 232-239.

PAN Y C, SUN C, LIU Y, TANG X M, REN Y M. Carrying capacity of livestock and poultry breeding based on feces disposal volume of land. Transactions of the Chinese Society of Agricultural Engineering, 2015, 31(4): 232-239. (in Chinese)

[34] 沈根祥, 錢曉雍, 梁丹濤, 顧海蓉, 姚政, 錢永清. 基于氮磷養(yǎng)分管理的畜禽場糞便匹配農(nóng)田面積. 農(nóng)業(yè)工程學(xué)報, 2006, 22(S2): 268-271.

SHEN G X, QIAN X Y, LIANG D T, GU H R, YAO Z, QIAN Y Q. Research on suitable area for cropland application of animal manure based on nutrients management of nitrogen and phosphorus. Transactions of the Chinese Society of Agricultural Engineering, 2006, 22(S2): 268-271. (in Chinese)

[35] 盧樹昌, 陳清, 張福鎖, 賈文竹. 河北果園主分布區(qū)土壤磷素投入特點及磷負荷風險分析. 中國農(nóng)業(yè)科學(xué), 2008, 41(10): 3149-3157.

LU S C, CHEN Q, ZHANG F S, JIA W Z. Characteristics of soil phosphorus input and phosphorus load risk in major orchards region of Hebei. Scientia Agricultura Sinica, 2008, 41(10): 3149-3157. (in Chinese)

[36] MA L, GUO J H, VELTHOF G L, LI Y M, CHEN Q, MA W Q, OENEMA O, ZHANG F S. Impacts of urban expansion on nitrogen and phosphorus flows in the food system of Beijing from 1978 to 2008. Global Environmental Change, 2014, 28: 192-204.

[37] GARRETT R D, NILES M T, GIL J D B, GAUDIN A, CHAPLIN- KRAMER R, ASSMANN A, ASSMANN T S, BREWER K, DE FACCIO CARVALHO P C, CORTNER O, DYNES R, GARBACH K, KEBREAB E, MUELLER N, PETERSON C, REIS J C, SNOW V, VALENTIM J. Social and ecological analysis of commercial integrated crop livestock systems: current knowledge and remaining uncertainty. Agricultural Systems, 2017, 155: 136-146.

[38] 付莉, 王貴云, 杜連柱, 楊柳, 張克強, 杜會英. 豬場廢水施用對設(shè)施白菜地土壤氮盈余的影響. 農(nóng)業(yè)現(xiàn)代化研究, 2021, 42(2): 215-222.

FU L, WANG G Y, DU L Z, YANG L, ZHANG K Q, DU H Y. Effects of swine wastewater manure application on soil nitrogen surplus characteristics of facility cabbage field. Research of Agricultural Modernization, 2021, 42(2): 215-222. (in Chinese)

[39] 楚天舒, 韓魯佳, 楊增玲. 考慮種養(yǎng)平衡的黃淮海小麥-玉米模式下畜禽承載量估算. 農(nóng)業(yè)工程學(xué)報, 2019, 35(11): 214-222.

CHU T S, HAN L J, YANG Z L. Livestock carrying capacity estimation in wheat-corn production model of Huang-Huai-Hai Region considering planting-raising balance. Transactions of the Chinese Society of Agricultural Engineering, 2019, 35(11): 214-222. (in Chinese)

[40] SáEZ J A, CLEMENTE R, BUSTAMANTE M á, YA?EZ D, BERNAL M P. Evaluation of the slurry management strategy and the integration of the composting technology in a pig farm - Agronomical and environmental implications. Journal of Environmental Management, 2017, 192: 57-67.

[41] 李文哲, 徐名漢, 李晶宇. 畜禽養(yǎng)殖廢棄物資源化利用技術(shù)發(fā)展分析. 農(nóng)業(yè)機械學(xué)報, 2013, 44(5): 135-142.

LI W Z, XU M H, LI J Y. Prospect of resource utilization of animal faeces wastes. Transactions of the Chinese Society for Agricultural Machinery, 2013, 44(5): 135-142. (in Chinese)

[42] 鄧遠遠, 郭焱, 朱俊峰. 政府規(guī)制下畜禽養(yǎng)殖廢棄物處理合規(guī)行為選擇. 中國農(nóng)業(yè)資源與區(qū)劃, 2021, 42(10): 95-102.

DENG Y Y, GUO Y, ZHU J F. Disposal of waste in livestock and poultry breeding compliance behaviors on the government regulations. Chinese Journal of Agricultural Resources and Regional Planning, 2021, 42(10): 95-102. (in Chinese)

[43] 佘磊, 姜珊, 姜彩紅, 李想, 柳王榮, 蘇文幸, 吳根義. 我國畜禽養(yǎng)殖環(huán)境管理進程及展望. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 2021, 40(11): 2277- 2282.

SHE L, JIANG S, JIANG C H, LI X, LIU W R, SU W X, WU G Y. Progress and prospect for environmental management of livestock and poultry breeding in China. Journal of Agro-Environment Science, 2021, 40(11): 2277-2282. (in Chinese)

Environmental Risks Assessment of Livestock and Poultry Non- Point Source Pollution in Shiyan City Based on Arable Land Carrying Capacity

GONG ShiFei1, XIAO NengWu1, DING WuHan2, JU XueHai3, WU PingHua1, YU YongSong1, LI Hu2

1Danjiangkou Reservoir of Shiyan Ecological Agriculture Institute, Shiyan 442000, Hubei;2Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences, Beijing 100081;3Rural Energy and Environment Agency, Ministry of Agriculture and Rural Affairs, Beijing 100125

【Objective】 Shiyan city was the core water source area of the middle route of the South-to-North Water Transfer Project, a comprehensive understanding of the current situation of livestock and poultry pollution and its potential threat to farmland would provide a basis for strengthening environmental risk control. 【Method】This paper used the livestock manure environmental risk assessment model to quantitatively analyze the source structure of livestock and poultry manure in different areas of Shiyan city in 2020, and the environmental risk from the perspective of the balance between planting and raising was comprehensively evaluated. 【Result】The total amount of livestock and poultry manure reached 306.85×104t in 2020, and it was mainly distributed in Yunyang district, Danjiangkou city and Fangxian county. COD production was the highest (57.00×104t), followed by TN (2.63×104t), and TP production was less (0.65×104t). The total equivalent standard pollution load was 15.532×1010m3, and the source structure was obviously inconsistent in different regions. TP was the primary pollutant in livestock and poultry breeding, and its equivalent standard pollution load accounts for 41.68% of the total load, followed by TN with the pollution load rate of 33.85%, and the lowest equivalent standard pollution load rate of COD was 24.47%. Cattle was the primary source of non-point source pollution in livestock and poultry production in Shiyan city, contributing 49.47% to the total pollution load of the equivalent standard, followed by poultry (21.60%). According to the result of environmental pollution risk assessment, the load of animal excretion per hectare of arable land in Shiyan city was 10.15 t·hm-2, the alarm value was 0.338 and the risk level was Ⅰ, indicating the environmental pollution was “no risk”. The load intensity of TN and TP was 0.09 and 0.02 t·hm-2, respectively, and the environmental risk index was 1.023 and 1.223, respectively, but there were obvious spatial differences among regions. More than half of the whole city had the potential pollution risk of livestock and poultry breeding, among which Maojian district and Yunyang District were the most serious. With phosphorus as the measurement standard, the total amount of breeding in Shiyan city was close to saturation, the total amount of livestock and poultry must be reduced by 98.51×104pigs equivalent on the existing basis.【Conclusion】It was of great significance to maintain a reasonable scale of poultry breeding (≤441.00×104pigs equivalent) and take measures of pollutant abatement, in order to promote the balance of planting and breeding and orderly development of livestock and poultry breeding in Shiyan city.

Shiyan city; arable land carrying capacity; livestock and poultry manure; cropland load; environmental risk

10.3864/j.issn.0578-1752.2023.05.009

2022-02-07；

2022-06-17

中國農(nóng)業(yè)科學(xué)院科技創(chuàng)新工程協(xié)同創(chuàng)新任務(wù)（CAAS-XTCX2016015）

龔世飛，E-mail：gongsfsf@163.com。通信作者肖能武，E-mail：413740790@qq.com。通信作者李虎，E-mail：lihu0728@sina.com

（責任編輯 李云霞）