牛學銳 韓美

摘要 [目的]研究中國肉雞規(guī)模化養(yǎng)殖在生產(chǎn)過程中對環(huán)境的影響。[方法]通過生命周期評價（LCA）的方法，研究大型肉雞養(yǎng)殖場的化肥生產(chǎn)、作物種植、飼料加工、肉雞飼養(yǎng)等階段的環(huán)境影響。[結(jié)果]整個過程中，全球變暖的主要因子是CO2和N2O，主要貢獻是化肥生產(chǎn)與種植階段；環(huán)境酸化的主要因子是NH3和NOX，其中NH3主要來源于化肥生產(chǎn)與種植和糞便處理階段；水體富營養(yǎng)化的主要污染因子是NH3；全過程主要環(huán)境污染階段是化肥生產(chǎn)與種植、肉雞飼養(yǎng)、糞便處理，主要的環(huán)境影響類型是環(huán)境酸化（78.8%），其次是全球變暖（12.8%）；結(jié)果與SimaPro軟件相比，發(fā)現(xiàn)兩者之間的共同性。[結(jié)論]肉雞飼養(yǎng)階段和環(huán)境酸化可作為減輕環(huán)境影響的目標。

關(guān)鍵詞 肉雞規(guī)模化養(yǎng)殖；生命周期評價；環(huán)境影響分析；SimaPro軟件

中圖分類號 X83文獻標識碼 A文章編號 0517-6611（2018）31-0048-05

Abstract [Objective]The research aimed to study the environmental impact of Chinese broilers largescale breeding in the production process. [Method]LCA method was applied to analyze the environmental impact of largescale broiler farms in the stages of fertilizer production，crop cultivation，feed processing and broiler raising. [Result]The main factors of global warming were CO2 and N2O， the main contributions were fertilizer production and planting stage.The main factors of environmental acidification were NH3 and NOX， the main source of NH3 was fertilizer production and cultivation and manure treatment stage.The main pollution factor of water eutrophication was NH3. In the whole process， the main environmental pollution stages were fertilizer production and planting， broiler raising and excrement disposal. The main types of environmental impact were acidification （78.8%）， followed by global warming （12.8%）.The results were compared with those of SimaPro software finding intercommunity between two.[Conclusion]The broiler feeding stage and environmental acidification can be used as targets for reducing environmental impact.

Key words Broiler largescale breeding；Life cycle assessment；Environmental impact analysis；SimaPro software

中國畜禽養(yǎng)殖業(yè)發(fā)展迅速，根據(jù)國家統(tǒng)計局1996—2015年的家禽出欄量統(tǒng)計數(shù)據(jù)分析，家禽數(shù)量迅速增長，增長幅度約100%[1]，近幾年中國肉雞產(chǎn)量的年變化正趨于平穩(wěn)，其波動取決于市場、人口穩(wěn)定緩慢增長等因素。2015年雞肉價格為2009年以來歷史最低，然而收益仍高達140億元[2]。通過對山東某養(yǎng)殖場調(diào)查，肉雞的飼養(yǎng)期糞便排泄量為5.5 kg，其中總磷含量為1.023%，總氮含量為0.413%[3]。按照2015年的專業(yè)雞肉出欄量90.44億只[2]計算，2015年肉雞糞便產(chǎn)量為4 974.2千萬t，總磷排放量為50.886萬t，總氮排放量為20.543萬t。不控制整個生產(chǎn)過程的污染，對大氣（溫室效應）、水（富營養(yǎng)化）、土壤（鹽漬化）污染嚴重。

農(nóng)業(yè)生產(chǎn)的環(huán)境污染一直備受人們關(guān)注，有關(guān)農(nóng)業(yè)的生命周期評價（LCA）于20世紀90年代開始在國外開展[4-6]，歐洲一些國家在同時期開展LCA在農(nóng)業(yè)方面的研究，如Brentrup等[6]利用LCA方法分析小麥生產(chǎn)系統(tǒng)中氮肥施用對環(huán)境的影響；Hospido等[7]運用LCA研究西班牙奶牛生產(chǎn)的環(huán)境影響；Thomassen等[8]也運用LCA研究荷蘭奶牛場的環(huán)境影響；Van等[9]關(guān)于法國奶牛生產(chǎn)過程的LCA。國外研究為中國LCA在農(nóng)業(yè)的應用提供方法經(jīng)驗，如中國開展農(nóng)業(yè)LCA研究的意義和應采取的對策[10]；牛奶生產(chǎn)的LCA[11]；玉米的LCA[12]；奶牛的LCA[13]；肉牛養(yǎng)殖的LCA[14]。有關(guān)肉雞規(guī)?；B(yǎng)殖的LCA相關(guān)研究僅有蛋雞養(yǎng)殖和肉雞屠宰場，肉雞與蛋雞在飼養(yǎng)周期及飼料配比上的區(qū)別，研究結(jié)果不適用于肉雞，而N.Pelletier[15]已經(jīng)對美國肉雞供應鏈進行過生命周期評價，為美國家禽業(yè)可持續(xù)發(fā)展提出改善建議，美國的飼養(yǎng)方式、養(yǎng)殖規(guī)模等方面與中國存在差異，應按照中國肉雞規(guī)?；B(yǎng)殖的例子進行具體研究。

肉雞養(yǎng)殖中糞便等污染物不合理利用和處理，造成周邊的水、大氣、土壤的污染[16]。筆者通過對肉雞養(yǎng)殖過程的生命周期評價，發(fā)現(xiàn)肉雞養(yǎng)殖過程中溫室效應、富營養(yǎng)化、酸化的環(huán)境影響以及環(huán)境負荷，并闡明肉雞養(yǎng)殖階段的環(huán)境效應，將分析結(jié)果與Simapro軟件計算結(jié)果相比較，印證計算準確性，為我國畜禽養(yǎng)殖的可持續(xù)發(fā)展提供科學理論與技術(shù)支持。

1 資料與方法

1.1 系統(tǒng)邊界

以山東某肉雞養(yǎng)殖場為目標，對肉雞規(guī)?；B(yǎng)殖進行生命周期評價，首先確定系統(tǒng)邊界[17]（圖1）。有關(guān)肉雞生產(chǎn)的系統(tǒng)邊界包括作物種植（化肥生產(chǎn)、玉米和大豆的種植）、豆粕生產(chǎn)、飼料加工、飼料運輸、肉雞飼養(yǎng)和糞便處理在內(nèi)的整個生產(chǎn)過程。肉雞生產(chǎn)的主要輸出產(chǎn)品為雞肉，功能單位為1 kg雞肉[18]。

1.2 生命周期評價數(shù)據(jù)清單

1.2.1 作物種植?；噬a(chǎn)過程的污染排放系數(shù)參考胡志遠等[19]的研究，僅計算生產(chǎn)過程中CO、NOX、SOX、CO2的排放。

玉米與大豆種植為并列階段，作物種植階段肉雞的飼料原料由玉米、豆粕、磷酸氫鈣、石粉等組成。根據(jù)實地考察，飼料原料僅考慮玉米和豆粕。作物種植階段主要包括玉米和大豆的種植階段。大豆、玉米的種植期間，主要的污染是化肥的釋放，其他污染農(nóng)藥和地膜等的生產(chǎn)數(shù)據(jù)和污染物排放數(shù)據(jù)無法收集，僅考慮化肥使用的環(huán)境污染[19]，根據(jù)研究目標，此過程考慮CO、NOX、SOX、CO2的排放量。

研究中飼料產(chǎn)地為山東，玉米和大豆種植期間的化肥投入相關(guān)數(shù)據(jù)取自高強等[20]、王偉[21]的研究。玉米、大豆的年平均單產(chǎn)量來自《中國農(nóng)業(yè)年鑒》[1]。玉米和大豆施肥過程中的污染排放量來自國外公開發(fā)表的文獻。經(jīng)查閱，目前國內(nèi)并沒有種植過程中磷肥和鉀肥的排放數(shù)據(jù)，只考慮施用氮肥對環(huán)境的影響。玉米種植期間，N2O、NOX、NO3排放量占氮肥總量的4.9%、5.3%、6.5%[22]。種植1 hm2的玉米NH3排放量是3.0kg[22]。大豆種植期間，N2O排放量占氮肥總量的1.25%，種植1 kg的大豆NH3排放量是2.9 g，種植1 hm2大豆NOX排放量是3.5 kg[23]。

1.2.2 豆粕生產(chǎn)。豆粕生產(chǎn)階段的豆粕作為大豆煉油過程中的輸出產(chǎn)品，生產(chǎn)豆油的過程中用煤燃燒加熱，大豆煉油過程中產(chǎn)生豆粕率80%[15]。根據(jù)董進寧等[24]的計算可知生產(chǎn)1 kg豆油消耗標煤0.0591 g，排放SO2 1.412 kg、CO2 176550 g、NOX 0.636 g。大豆在煉油過程中豆粕的環(huán)境負擔占過程總負擔的64%，豆油的環(huán)境負擔占過程總負擔的36%[15]，可推算1 kg豆粕生產(chǎn)過程中消耗標煤0.014 78 kg，污染物的排放量是0.627 6 g SO2、78.74 g CO2、0.283 g NOX（1 kg標準煤的熱值為29.27 MJ）[24]。

1.2.3 飼料加工。飼料加工階段主要消耗電能，實地調(diào)查平均每噸飼料消耗電能28.5（kW·h）/t。飼料加工過程中電能的污染物排放系數(shù)來自狄向華等[25]的研究，根據(jù)研究目標，此過程僅考慮CO、NOX、SOX、COX的排放量。

1.2.4 飼料運輸。整個生命周期的車輛運輸只包括飼料廠到養(yǎng)殖場的飼料運輸，實際調(diào)查中具體公里數(shù)無法統(tǒng)計。因1 kg雞肉中運輸階段所產(chǎn)生的環(huán)境影響很小，所以可采用調(diào)查養(yǎng)雞場經(jīng)驗數(shù)據(jù)，運輸距離160 km，主要考慮CO2、CO、SO2、NOX的排放[26]。污染物排放計算公式：

1.2.5 肉雞飼養(yǎng)。肉雞飼養(yǎng)階段向環(huán)境釋放的污染物主要分為肉雞自身和其他能源消耗。根據(jù)實地調(diào)研，整合飼養(yǎng)階段數(shù)據(jù)得出：飼養(yǎng)周期為42 d；飼料為省內(nèi)固定飼料廠購買，運輸距離為160 km；每只雞消耗飼料4.5 kg；體重2 kg，出肉率70%；玉米占總飼料量的58.6%，豆粕27.5%。計算得一只肉雞整個飼養(yǎng)周期，消耗玉米2.637 kg、豆粕1.238 kg，大豆煉油過程中產(chǎn)生豆粕率80%[15]，消耗大豆1.548 kg，生產(chǎn)1 kg雞肉消耗飼料3.21 kg（其中消耗玉米1.8 kg、大豆1.1 kg）。

飼養(yǎng)過程中的照明、取暖、添水等產(chǎn)生電費與煤費，每只雞的電力與燃煤費用分別為0.13、0.22元[27]，電價以0.64元/（kW·h），煤價以600元/t計算，推算出1 kg雞肉對玉米、大豆、電力、燃煤的消耗量分別為1.88 kg、1.10 kg、0.15 kW·h、0.26 kg。養(yǎng)雞生產(chǎn)過程中的煤電生產(chǎn)和燃燒的污染物排放系數(shù)來自胡志遠等[19]的研究。

根據(jù)Coufal等[28]的研究可知，1 000 kg肉雞在飼養(yǎng)過程中分別向空氣中排放1 kg N2O、11 kg NH3和1 kg CH4，計算得1 kg雞肉在飼養(yǎng)過程分別向空氣排放1.43 g N2O、15.70 g NH3、1.43 g CH4。

1.2.6 糞便處理?；谡g氣候變化專門委員會（IPCC）專家組的裁決，在雞糞便處理過程中的N2O直接排放量占總氮排放量的0.1%，從雞糞中揮發(fā)出的NH3和排放的NOX總量占總氮的40%，從雞糞中揮發(fā)出的NH3量占總氮的304%，則排放的NOX總量占總氮的9.6%[29]。飼養(yǎng)階段中糞便的日均產(chǎn)生量為0.1 kg/a，總氮含量1.023%[3]。糞便處理階段中，CH4的排放因子（平均氣溫為15～25 ℃）為每只肉雞生命周期內(nèi)排放的糞便產(chǎn)生0.02 kg CH4。結(jié)合飼養(yǎng)階段數(shù)據(jù)可計算1 kg雞肉在糞便處理階段的污染排放。

1.3 環(huán)境影響計算

1.3.1 特征化。選擇特征化中3種類別的影響因子[30-31]，由各影響類型及其污染物的影響當量系數(shù)，結(jié)合清單分析結(jié)果，特征化數(shù)值計算公式：

1.3.2 標準化。標準化有2個目的：一是通過規(guī)定一組環(huán)境影響潛值標準人當量基準值，進而比較污染物對多種影響方式的貢獻大?。欢菍τ谙乱徊降脑u價估算數(shù)據(jù)支撐。標準化計算公式：

1.3.3 加權(quán)。加權(quán)評估是將各影響類別經(jīng)專家討論賦予不同的權(quán)重值與標準化數(shù)值計算。加權(quán)計算的公式：

2 結(jié)果與分析

根據(jù)清單分析結(jié)果，在整個生命周期中，溫室氣體、富營養(yǎng)化、酸化作為該研究的環(huán)境影響分析類型。

2.1 特征化

2.1.1 全球變暖影響潛值分析。按公式（2）計算各階段全球變暖影響潛值，結(jié)果發(fā)現(xiàn)（圖2），整個生命周期中，全球變暖影響潛值方面貢獻最大的是作物種植階段，數(shù)值為1 852.04 g（CO2）eq/kg，占總量的65.9%；第二的是肉雞飼養(yǎng)階段，數(shù)值為674.37 g（CO2）eq/kg，占總量的24%。影響全球變暖的主要污染因子是N2O和CO 分別占總溫室氣體的56%、41%。CO2來源于化肥生產(chǎn)與種植階段。N2O主要來源于作物種植，原因是作物種植中化肥生產(chǎn)會排放大量的溫室氣體；而肉雞飼養(yǎng)階段考慮到飼養(yǎng)過程中溫室氣體的排放，所以總量占比第二。

2.1.2 酸化影響潛值分析?；谇鍐畏治鼋Y(jié)果，按上述計算各階段的酸化影響潛值，結(jié)果發(fā)現(xiàn)（圖3），整個生命周期在酸化影響方面貢獻最大的是肉雞飼養(yǎng)，數(shù)值為29.72 g（SO2）eq/kg，占總量的40.0%；其次是作物種植，數(shù)值為23.54 g（SO2）eq/kg，占總量的31.7%；糞便處理的貢獻第三，數(shù)值為19.77 g（SO2）eq/kg，占總量的26.6%。環(huán)境酸化的主要影響因子是NH3和NOX，從清單分析NH3的主要來源是作物種植和糞便處理階段，NOX主要來源是種植過程中氮肥的施加以及糞便中含氮物質(zhì)的揮發(fā)。

2.1.3 富營養(yǎng)化環(huán)境影響潛值分析。根據(jù)清單分析數(shù)據(jù)，計算各階段的富營養(yǎng)化環(huán)境影響潛值，結(jié)果發(fā)現(xiàn)（圖4），在富營養(yǎng)化方面，作物種植階段、肉雞飼養(yǎng)階段、糞便處理階段的影響較大，分別占總量的38.5%、36.5%、24.3%。富營養(yǎng)化的主要污染因子是NH3，清單分析中富營養(yǎng)化的主要來源是肉雞飼養(yǎng)與糞便處理。

2.2 標準化

經(jīng)查閱文獻，采用的環(huán)境影響標準是人當量基準值[31]。根據(jù)2015年中國肉雞產(chǎn)業(yè)的檢測，中國專業(yè)性肉雞產(chǎn)量為1 255萬t[2]，與已計算的3種環(huán)境影響潛值的特征化數(shù)值和當量基準值帶入公式（3），可計算有關(guān)肉雞規(guī)?；a(chǎn)的全生命周期影響潛值的標準化數(shù)值，結(jié)果見表1。

從表1可以看出，酸化對環(huán)境影響最大，為252.49×108；其次是全球變暖，為40.00×108。受肉雞養(yǎng)殖污染影響的地區(qū)，可從環(huán)境酸化方向入手改善當?shù)丨h(huán)境。

有關(guān)肉雞規(guī)?；B(yǎng)殖各階段標準化數(shù)值比較結(jié)果基于清單分析計算得出，結(jié)果見表2。肉雞飼養(yǎng)、作物種植、糞便處理是對環(huán)境影響較大的前3個階段；在環(huán)境影響最嚴重的環(huán)境酸化中，也是肉雞飼養(yǎng)、作物種植、糞便處理占前三。

2.3 加權(quán)

經(jīng)查閱文獻得到環(huán)境影響的權(quán)重值[28]。由肉雞各生產(chǎn)階段標準值（表2），計算各階段加權(quán)結(jié)果見表3。肉雞生命周期的環(huán)境影響總負荷（EI）為74.34×108，作物種植占38.8%，肉雞飼養(yǎng)占36.7%，糞便處理占21.7%。結(jié)合標準值分析肉雞飼養(yǎng)階段全球變暖和環(huán)境酸化的貢獻為913%；作物種植占各階段總影響的36.6%，結(jié)合標準化分析，該階段污染物的環(huán)境酸化和全球變暖貢獻為90.5%；飼料加工和運輸僅占各階段的1.6%。

2.4 SimaPro軟件計算結(jié)果比較

SimaPro軟件作為全球范圍內(nèi)權(quán)威的生命周期評價軟件，具有計算的科學性[32]，以SimaPro軟件的計算結(jié)果為輔助，分析該研究的合理性。通過SimaPro軟件計算結(jié)果為：①特征化。全球變暖（1699 7 kg GWP）、環(huán)境酸化（0.034 2 kg AP）、水體富營養(yǎng)化（0.013 3 kg NP）。②標準化。全球變暖（0.000 185）、環(huán)境酸化（0.000 975）、水體富營養(yǎng)化（0.000 375）。

經(jīng)計算，各類型影響潛值在3個方面均高于SimaPro軟件的計算結(jié)果。原因是此研究考慮到飼養(yǎng)階段氣體排放和糞便的污染物揮發(fā)，所以結(jié)果高于SimaPro軟件，但是2種計算結(jié)果存在共同性：全球變暖的影響潛值最高，其次是環(huán)境酸化，最后是水體富營養(yǎng)化。

SimaPro軟件對于數(shù)據(jù)的標準化分析的結(jié)果（圖5），環(huán)境酸化、水富營養(yǎng)化、全球變暖（溫室氣體）分別占總影響的64%、24%、12%，與表2的標準化結(jié)果相比較后發(fā)現(xiàn)，SimaPro軟件中環(huán)境酸化與全球變暖的影響占比低于該研究的計算，因為USLCI數(shù)據(jù)庫[27]在限定條件下沒有考慮肉雞飼養(yǎng)階段的氣體排放以及糞便的污染物揮發(fā)，且這2個方面的主要環(huán)境貢獻為環(huán)境酸化，所以在結(jié)果上全球變暖的占比升高，環(huán)境酸化的占比降低?？偟膩碚f，2個數(shù)據(jù)結(jié)果雖有不同，但科學可靠，該研究可作為科學理論數(shù)據(jù)用來指導政策和生產(chǎn)。

從各階段來看，作物種植、肉雞飼養(yǎng)、糞便處理不管是在特征化、標準化、加權(quán)計算都占比較高；主要的環(huán)境影響類型是環(huán)境酸化，其次是全球變暖。

3 小結(jié)與討論

基于生命周期評價對肉禽規(guī)模化養(yǎng)殖進行全球變暖、環(huán)境酸化、水體富營養(yǎng)化3方面環(huán)境影響分析，經(jīng)過特征化、標準化、加權(quán)分析后看出，整個生命周期對環(huán)境貢獻從大到小依次為環(huán)境酸化、全球變暖、水體富營養(yǎng)化。從階段來看，作物種植在加權(quán)評估階段占比最高，肉雞飼養(yǎng)次之；肉雞飼養(yǎng)在標準化分析占比最高，作物種植階段次之，原因是該階段全球變暖貢獻占比高于肉雞飼養(yǎng)階段中的，所以加權(quán)評估后化肥生產(chǎn)與種植階段會高于肉雞飼養(yǎng)階段，但很明顯這2個階段最主要的環(huán)境貢獻都是環(huán)境酸化。

減輕肉雞規(guī)?；B(yǎng)殖的環(huán)境影響，在生產(chǎn)過程中最可行的方法是加強飼養(yǎng)階段的污染物排放控制；治理過程中，對受肉雞養(yǎng)殖污染的環(huán)境，治理角度從環(huán)境酸化入手效果明顯。

該研究不足之處是清單收集不夠完善，磷肥、鉀肥的污染排放缺少可靠數(shù)據(jù)無法計算；飼料中非主要原料的環(huán)境影響數(shù)據(jù)缺失，也因研究條件沒有計算在內(nèi)；運輸距離采用調(diào)查過程中的經(jīng)驗值，運輸環(huán)境影響相對很小，對研究結(jié)果的影響不大。我國農(nóng)業(yè)生產(chǎn)方式與其他國家存在差異，現(xiàn)很多排放系數(shù)與統(tǒng)計數(shù)據(jù)使用國外的研究成果，研究成果受到影響，期望在原材料和能源為基礎，逐步向上建立數(shù)據(jù)清單。很多工藝與原材料投加量屬于商業(yè)機密，在研究過程中清單數(shù)據(jù)的獲取十分困難，數(shù)據(jù)清單的建立需要行業(yè)內(nèi)專業(yè)人員或企業(yè)自身來完成。

參考文獻

[1] 中國農(nóng)業(yè)年鑒編輯委員會.中國農(nóng)業(yè)年鑒[M].北京：中國農(nóng)業(yè)出版社，2010.

[2] 鄭麥青，宮桂芬，高海軍，等.2015年我國肉雞產(chǎn)業(yè)發(fā)展監(jiān)測報告[J].中國家禽，2016，38（7）：67-70.

[3] 侯世忠.山東畜禽糞便產(chǎn)生量及其環(huán)境效應分析[J].山東畜牧獸醫(yī)，2013，34（2）：1-3.

[4] VAN DER WERF H M G，PETIT J.Evaluation of the environmental impact of agriculture at the farm level：A comparison and analysis of 12 indicator-based methods[J].Agriculture，ecosystems & environment，200 93（1/2/3）：131-145.

[5] KRAMER K J，MOLL H C，NONHEBEL S.Total greenhouse gas emissions related to the Dutch crop production system[J].Agriculture，ecosystems & environment，1999，72（1）：9-16.

[6] BRENTRUP F，KUSTERS J，KUHLMANN H，et al.Environmental impact assessment of agricultural production systems using the life cycle assessment methodology：Ⅰ.Theoretical concept of a LCA method tailored to crop production[J].European journal of agronomy，2004，20（3）：247-264.

[7] HOSPIDO A，MOREIRA M T，F(xiàn)EIJOO G.Simplified life cycle assessment of galician milk production[J].International dairy journa 2003，13（10）：783-796.

[8] THOMASSEN M A，DOLMAN M A，VAN CALKER K J，et al.Relating life cycle assessment indicators to gross value added for Dutch dairy farms[J].Ecological economics，2009，68（6）：2278-2284.

[9] VAN DER WERF H M G，KANYARUSHOKI C，CORSON M S.An operational method for the evaluation of resource use and environmental impacts of dairy farms by life cycle assessment[J].Journal of environmental management，2009，90（11）：3643-3652.

[10] 楊印生，盛國輝，呂廣宏.我國開展農(nóng)業(yè)LCA研究的對策建議[J].中國軟科學，2003（5）：7-11.

[11] 秦鳳賢.生命周期評價在牛奶生產(chǎn)中的應用研究[J].乳業(yè)科學與技術(shù)，2006，29（5）：224-226.

[12] 張麗麗，劉鳴達.遼西地區(qū)兩種玉米生產(chǎn)方式的生命周期評價：以遼寧省建平縣為例[J].沈陽農(nóng)業(yè)大學學報，201 42（3）：300-305.

[13] 王效琴，梁東麗，王旭東，等.生命周期評價方法評估奶牛養(yǎng)殖系統(tǒng)溫室氣體排放量[J].農(nóng)業(yè)工程學報，201 28（13）：179-184.

[14] 程瓊儀，張言，陳昭輝，等.肉牛場生命周期估計及環(huán)境影響評價[J].家畜生態(tài)學報，2015，36（1）：63-69.

[15] PELLETIER N.Environmental performance in the US broiler poultry sector：Life cycle energy use and greenhouse gas，ozone depleting，acidifying and eutrophying emissions[J].Agricultural systems，2013，98（2）：67-73.

[16] STEINFELD H，GERBER P，WASSENAAR T D，et al.Livestock's long shadow：Environmental issues and options[M].Rome：Food and Agriculture Organization of the United Nations，2006.

[17] ISO.Environmental managementlife cycle assessmentprinciples and framework：ISO 14040[S].Geneva，Switzerland：International Standard Organization，1997.

[18] Layer chicken Sciences Group.Kwantitatieve informatie veehouderij 2006-2007[M].The Netherlands：Lelystad，2006.

[19] 胡志遠，譚丕強，樓狄明，等.不同原料制備生物柴油生命周期能耗和排放評價[J].農(nóng)業(yè)工程學報，2006，22（11）：141-146.

[20] 高強，馮國忠，王志剛.東北地區(qū)春玉米施肥現(xiàn)狀調(diào)查[J].中國農(nóng)學通報，2010，26（14）：229-231.

[21] 王偉.黑龍江省大豆合理施肥參數(shù)的研究[D].哈爾濱：東北農(nóng)業(yè)大學，2009.

[22] KNUDSEN M T，QIAO Y H，LUO Y，et al.Environmental assessment of organic soybean（Glycine max.）imported from China to Denmark：A case study [J].Journal of cleaner production，2010，18（14）：1431-1439.

[23] RENOUF M A，WEGENER M K，NIELSEN L K.An environmental life cycle assessment comparing Australian sugarcane with US com and UK sugar beet as producers of sugars for fermentation[J].Biomass and bioenergy，2008，32（12）：1144-1155.

[24] 董進寧，馬曉茜.生物柴油項目的生命周期評價[J].現(xiàn)代化工，2007，27（9）：59-63.

[25] 狄向華，聶祚仁，左鐵鏞.中國火力發(fā)電燃料消耗的生命周期排放清單[J].中國環(huán)境科學，2005，25（5）：632-635.

[26] 楊建新，徐成，王如松.產(chǎn)品生命周期評價方法及應用[M].北京：氣象出版社，2002.

[27] 國家發(fā)展和改革委員會價格司.全國農(nóng)產(chǎn)品成本收益資料匯編[M].北京：中國統(tǒng)計出版社，2014.

[28] COUFAL C D，CHAVEZ C，NIEMEYER P R，et al.Nitrogen emissions from broiler measured by mass balance over eighteen consecutive flocks[J].Poult Sci，2006，85（3）：384-391.

[29] DONG H M，MANGINO J，MCALLISTER T A，et al.Emissions from livestock and manure management[M]//EGGLESTON S，BUENDIA L，MIWA K，et al.2006 IPCC guidelines for national greenhouse gas inventories.Kanagawa，Japan：Institute for Global Environmental Strategies，2006.

[30] 高雪.黑龍江省集約化蛋雞場的生命周期評價[D].哈爾濱：東北農(nóng)業(yè)大學，2012.

[31] 馬宗虎，王美芝，丁露雨，等.規(guī)?；馀Ｓ蕡鰷厥覛怏w排放的生命周期評估[J].農(nóng)業(yè)環(huán)境科學學報，2010，29（11）：2244-2252.

[32] YANG J X，NIELSEN P H.Chinese life cycle impact assessments factors[J].Journal of environmental sciences，200 13（2）：205-209.