長株潭區域生物質開放燃燒的大氣污染物排放清單及不確定性分析*

劉 湛 李貝睿， 尤翔宇 黃 妍(.湖南省環境保護科學研究院，湖南 長沙 40004；.湘潭大學化工學院，湖南 湘潭 405)

早在1976年，CRUTZEN等[1]最早提出生物質燃燒會對大氣環境產生影響，近年來，由于生物質燃燒引起的區域大氣污染事件時有發生，再次引起社會和學者的廣泛關注。尹聰等[2]利用遙感火點信息，結合氣象及空氣質量監測資料，通過后向軌跡模擬，發現秸稈焚燒排放的污染物的區域尺度輸送與局地污染結合會造成城市嚴重空氣污染事件。朱佳雷等[3]通過區域大氣環境模式系統(RegAEMS)模式模擬長三角地區秸稈焚燒污染物排放情況，發現秸稈焚燒污染物排放導致大氣中PM10、CO濃度升高30%(質量分數)以上，黑炭、有機物所占消光比例有所增加。DUAN等[4]通過研究大氣顆粒物中有機碳(OC)/水溶性K+，提出生物質燃燒對北京地區灰霾污染有明顯影響。

目前我國有關生物質開放燃燒大氣污染物排放的研究仍然較少。唐喜斌等[5]利用自己設計開發的開放式燃燒源排放測試系統，測試了各類秸稈的排放因子與顆粒物的成分譜。YAN等[6]利用統計、調查、衛星火點等數據，估算了全國生物質燃燒量。對于生物質燃燒排放清單的研究，目前主要集中于宏觀層面，如洲際間相互影響[7]、亞洲區域[8]、全國區域[9-11]等，微觀的城市或城市群生物質開放燃燒排放清單研究多集中在北京[12]、長三角[13]、廣州[14-15]地區，對中部地區研究仍處于空白。清單不確定性研究多停留在定性識別、半定量評估階段，鮮有定量評估的研究。環境保護部于2014年底發布了《生物質燃燒源大氣污染物排放清單編制技術指南(試行)》(以下簡稱《指南》)等一系列規范，給出了生物質燃燒大氣污染物排放清單的編制方法。

長株潭(長沙市、株洲市、湘潭市)區域是我國典型的中部城市群，農業、工業、旅游業發達，作為國家“兩型”社會示范區域，隨著經濟持續高速發展，長株潭區域在國家“中部崛起”戰略中有著舉足輕重的地位。筆者以該區域為研究對象，廣泛收集區域內秸稈產量、森林火災消耗生物量、露天焚燒比例、燃燒效率等信息，建立了2012年長株潭區域生物質開放燃燒的大氣污染物排放清單，對清單做了時間分配和空間分布，并對結果做了定量化的不確定性分析，以期對長株潭區域空氣質量持續改善提供決策依據，為完善全國精細化排放清單提供區域數據參考，并為后續區域大氣排放清單的完善提供基礎。

1 研究方法

1.1 研究范圍

以2012年為基準年，分別計算長株潭區域23個縣區生物質開放燃燒排放的大氣污染物SO2、NOX、PM10、PM2.5、揮發性有機物(VOCs)、CO、元素碳(EC)、OC量。生物質開放燃燒源分為秸稈露天焚燒、森林火災、草原火災3類。長株潭區域地處長江流域，以丘陵和山地為主要地形，地表草原面積比例很小，因此本研究忽略該區域草原火災，分別以秸稈露天焚燒和森林火災的大氣污染物為對象建立排放清單。

1.2 計算方法

生物質開放燃燒源難以獲取固定排放位置和活動水平，統一按照面源處理。采用“自下而上”的方法，根據區域內生物質開放燃燒活動水平與排放因子數據計算，某種大氣污染物的排放量(Ei，t)的計算模型見式(1)。

(1)

式中：M為排放源活動水平，t；EF為排放因子，g/kg；i為某種大氣污染物；j為地區(市、縣)；k為生物質開放燃燒類型；m為植被帶/秸稈類型。

1.3 活動水平

農作物產量數據來自湖南省農村年鑒(2013)，長株潭區域秸稈以水稻、玉米為主，具體秸稈產量及露天焚燒量見表1。秸稈露天焚燒量計算見式(2)。

M秸稈=P×N×R×η

(2)

式中：M秸稈為秸稈露天焚燒量，t/a；P為農作物產量，t/a；N為草谷比；R為秸稈露天焚燒比例；η為燃燒效率。

秸稈的燃燒效率、草谷比及露天焚燒比例參考現有的相關研究成果[16-18]。其中，燃燒效率取值0.9，各類農作物平均草谷比取值見表2。不同研究給出的各省秸稈露天焚燒比例差異很大，本研究通過抽樣調查、統計分析，并結合前人在類似自然、氣候條件區域的研究成果，秸稈露天焚燒比例取值為0.2。

長沙市、株洲市、湘潭市2012年森林火災消耗生物量計算見式(3)。

M森林=A×D×η

(3)

式中：M森林為森林火災消耗的生物量，t；A為火災受害面積，hm2；D為森林干生物量，t/hm2。

2012年長株潭區域森林火災受害面積來源于國家林業局的年度森林火災統計資料，其數值見表3。由于該資料是按省級區域進行統計，根據湖南省所處的植被氣候帶，選取相應的生物量，再按照長沙市、株洲市、湘潭市的森林植被面積進行分配，因此森林火災排放清單以市一級行政單位建立。植被帶平均生物量及平均燃燒效率引用中國林業科學院的研究成果。

1.4 排放因子

排放因子的獲取方式為文獻調研法，由于缺乏可靠的本地實測資料，本研究匯總了國內外現有研究所選取的排放因子[19-23]及《指南》中排放因子數據。按照長株潭區域實際的地理、氣候情況，選取了相應的適宜本地的排放因子，見表4。

表1 2012年長株潭區域秸稈產量及露天焚燒量Table 1 Output of crop residue and biomass consumption of open-burning in Changsha-Zhuzhou-Xiangtan district in 2012 t/a

表2 各類農作物平均草谷比Table 2 Average grass valley ratio of different kinds of crops

表3 2012年長株潭區域森林火災消耗生物量Table 3 Biomass consumption of forest fire in Changsha-Zhuzhou-Xiangtan district in 2012

2 結果與討論

2.1 計算結果

根據生物質開放燃燒源排放因子的計算方法，估算得到2012年長株潭區域生物質開放燃燒大氣污染物排放清單，見表5。由表5可知，2012年生物質開放燃燒大氣污染物SO2、NOX、PM10、PM2.5、VOCs、CO、EC、OC排放量分別為783.48、4 248.00、10 325.94、10 117.29、6 882.92、76 002.99、816.09、3 478.28 t，其中秸稈露天焚燒排放量所占比例分別為96.76%、98.33%、96.00%、96.00%、94.77%、93.92%、96.40%、93.35%，為最主要排放源。

2.2 時間分配

根據實際調查及文獻調研，得到長株潭區域生物質開放燃燒的月度時間分配曲線，如圖1所示。從圖1可以看出，秸稈露天焚燒量在7、10月出現2個峰值，原因是長株潭區域水稻種植一般為一年2季，收獲時間分別為7、10月，而秋季，尤其是10月，也是大部分農作物的收獲季節，因此秸稈產量和露天焚燒量在10月達到最高值。森林火災消耗生物量呈季節性階梯變化，夏、秋季氣溫高，降水頻率小，極易發生火災，加上人類生活生產周期的影響，森林火災集中出現在夏、秋季，冬季氣溫低，春季降水多，都不易發生森林火災，因此分配占比低。由此，秸稈露天焚燒管理工作重點應放在7—10月，森林火災的防治工作應集中在夏、秋季開展。

表4 生物質開放燃燒源排放因子Table 4 Emission factors of biomass open-burning g/kg

表5 2012年長株潭區域生物質開放燃燒大氣污染物排放清單Table 5 Emission inventory of air pollutant from biomass open-burning in Changsha-Zhuzhou-Xiangtan district in 2012 t

2.3 空間分布

以PM2.5為代表污染物，長株潭區域秸稈露天焚燒污染物排放量按縣區空間分布如圖2所示，其他污染物空間分布與其類似。結合表5、圖2可以明顯看出，長株潭區域污染物排放量最小的3個地區為長沙市雨花區、芙蓉區、天心區，其他較小的區域包括長沙市開福區，湘潭市岳塘區，株洲市石峰區、荷塘區和蘆淞區，全部為城市區域，且相互接壤。另外，韶山市秸稈露天焚燒PM2.5排放量較小，其為我國典型的紅色旅游區，以第三產業為經濟支柱，農業發展相對落后，加上風景旅游區為環境敏感區，污染物排放標準高、地區控制嚴格，所以秸稈露天焚燒量小。長沙市寧鄉縣和湘潭市湘潭縣耕地面積廣、農作物產量高、秸稈產量大，因此PM2.5排放量大。總而言之，生物質開放燃燒污染物排放量最大的縣區為長沙市寧鄉縣，其次為湘潭市湘潭縣和長沙市瀏陽市，此3縣區特點為地域面積大、耕地面積廣、秸稈產量大，秸稈露天焚燒污染物排放量大。森林火災大氣污染物排放量最大的為株洲市，株洲市中部和南部森林覆蓋面積廣，人口密度小，森林火災頻率較大。而在長株潭中心區域(長沙市開福區、湘潭市岳塘區、株洲市石峰區、荷塘區和蘆淞區等城市區域)形成一片污染物排放量較小的區域。

圖1 生物質開放燃燒的月度時間分配曲線Fig.1 Time distribution of biomass open-burning

圖2 2012年長株潭區域秸稈露天焚燒PM2.5排放空間分布Fig.2 Spatial distribution of PM2.5 from field crop residue burning in Changsha-Zhuzhou-Xiangtan district in 2012

2.4 不確定性分析

2.4.1 輸入信息不確定性評價

科學評價一個數據的不確定性方法一般為：隨機抽取樣本庫中樣本，通過數學方法，獲得其概率密度函數，并包含分布形式、平均值、標準方差等信息，它的不確定性通過相對標準偏差的概念表達。我國排放清單建立過程中所選取的數據，如活動水平信息，大多來自各類統計公報，一般僅存在一個有效數值；排放因子數據多取自前人的測試或模式模擬結果，缺乏足夠數量的數據組成有效數據庫。

本研究假設輸入的農作物產量、草谷比、秸稈露天焚燒比例、火災受害面積、森林火災消耗生物量、排放因子等輸入因子均呈正態分布形式。根據美國環境保護署(USEPA)提出的基于測試次數的不確定性等級分類，對于直接源于農業年鑒、統計年鑒及國家林業局的年度森林火災統計資料等官方統計部門資料的取值，其不確定度取±30%。通過其他可靠統計數據，或者根據官方統計部門資料數據得出的相應轉化系數的取值，其不確定度取±50%。

本研究排放因子數據來自《指南》及國內外學者的實測數據，盡管測試樣本及自然條件與本區域有所差異，但是秸稈燃燒行業差別不大，秸稈類型全國各地區都較相似，相對而言，排放因子具有較高的可信度。根據USEPA提出的基于測試次數的不確定性等級分類，本研究排放因子級別大致為B～A級。

2.4.2 不確定性結果輸出

使用蒙特卡羅法，使排放因子和活動水平等各參數的不確定性通過式(1)傳遞到模型的輸出，從而定量評估各污染物排放量的不確定性范圍。通過10 000次抽樣計算，取95%置信區間，其分析結果如表6所示。由表6可知，秸稈露天焚燒SO2、NOX、PM10、PM2.5、VOCs、CO、EC、OC抽樣計算平均值分別為765.43、4 202.83、9 961.56、9 748.81、6 524.77、71 753.28、790.46、3 282.81 t，與清單計算結果接近，不確定性總體在-84%～168%。森林火災SO2、NOX、PM10、PM2.5、VOCs、CO、EC、OC抽樣計算平均值分別為25.38、70.93、413.84、403.93、359.77、4 623.36、29.39、231.34 t，與清單計算結果接近，不確定性總體在-83%～176%。

2.4.3 單個污染物不確定性模擬

選取秸稈露天焚燒PM2.5排放量為代表，對生物質露天焚燒源單個污染物做深入分析。通過10 000次重復計算，得到的秸稈露天焚燒PM2.5排放量的概率密度分布，突破了傳統的排放清單計算僅獲得一個計算數值的局限性，其詳細結果見圖3。由于假設各輸入數據為正態分布，PM2.5排放量概率密度函數呈對數正態分布。計算得到PM2.5排放量平均值為9 748.81 t，中位值比平均值略低，為平均值的84.85%。輸出結果顯示，PM2.590%置信區間的不確定性為-76.08%～125.76%，95%置信區間的不確定性為-83.07%～167.18%。

2.4.4 不確定性貢獻分析

在分析排放清單的不確定性時，研究對不確定性貢獻最大的一個來源具有較大的實際意義。因此，比較各不確定性來源對排放清單不確定性的貢獻率非常重要，常用敏感性分析方法實現。研究利用蒙特卡羅法計算秸稈露天焚燒PM2.5排放量過程中所產生的10 000組模擬數值，進一步分析輸入變量對輸出結果的敏感程度，結果見圖4。從圖4可知，露天焚燒比例是結果不確定性的最大來源，貢獻率為24.8%，其次，燃燒效率、水稻產量、水稻草谷比等也是不確定性的重要來源，貢獻率為23.7%、20.1%和20.0%，這4個因子均決定了活動水平的輸入，因此活動水平數據的精準程度會大幅影響排放清單的可靠性。PM2.5排放因子對不確定性的貢獻率為11.0%，對結果也會產生較大影響。

表6 生物質開放燃燒大氣污染物排放不確定性分析結果Table 6 The result of uncertainty analysis of pollutant emissions from biomass open-burning

圖3 2012年秸稈露天焚燒PM2.5排放量不確定性模擬結果Fig.3 Simulation results of uncertainty of PM2.5 emissions from field burning of crop residue for year 2012

圖4 各因子對PM2.5排放量不確定性貢獻率Fig.4 Contributions of various parameters to variability of estimated PM2.5 emission from field crop residue burning

3 結 論

(1) 2012年長株潭區域生物質開放燃燒大氣污染物SO2、NOX、PM10、PM2.5、VOCs、CO、EC、OC排放量分別為783.48、4 248.00、10 325.94、10 117.29、6 882.92、76 002.99、816.09、3 478.28 t。

(2) 秸稈露天焚燒大氣污染物排放量在7、10月形成2個峰值，森林火災集中出現在夏、秋季。

(3) 生物質開放燃燒污染物排放量最大的縣區為長沙市寧鄉縣，其次為湘潭市湘潭縣和長沙市瀏陽縣。而在長株潭中心區域(長沙市開福區、湘潭市岳塘區、株洲市石峰區、荷塘區和蘆淞區等城市區域)形成一片污染物排放量較小的區域。

(4) 采用蒙特卡羅法，計算得到2012年長株潭區域秸稈露天焚燒源和森林火災大氣污染物排放量95%置信區間的不確定性分別為-84%～168%、-83%～176%。通過定量模擬得到秸稈露天焚燒PM2.5排放量概率密度函數呈對數正態分布。

(5) 以PM2.5為代表污染物，對其排放量的不確定性貢獻率最大的是露天焚燒比例，其次，燃燒效率、水稻產量、水稻草谷比等也是不確定性的重要來源，活動水平數據的精準程度會大幅影響排放清單的可靠性。

[1] CRUTZEN P J，HEIDT L E，KRASNEC J P，et al.Biomass burning as a source of atmospheric gases CO，H2，N2O，NO，CH3Cl and COS[J].Nature，1979，282(5736)：253-256.

[2] 尹聰，朱彬，曹云昌，等.秸稈焚燒影響南京空氣質量的成因探討[J].中國環境科學，2011，31(2)：207-213.

[3] 朱佳雷，王體健，鄧君俊，等.長三角地區秸稈焚燒污染物排放清單及其在重霾污染天氣模擬中的應用[J].環境科學學報，2012，32(12)：3045-3055.

[4] DUAN Fengkui，LIU Xiande，YU Tong，et al.Identification and estimate of biomass burning contribution to the urban aerosol organic carbon concentrations in Beijing[J].Atmospheric Environment，2004，38(9)：1275-1282.

[5] 唐喜斌，黃成，樓晟榮，等.長三角地區秸稈燃燒排放因子與顆粒物成分譜研究[J].環境科學，2014，35(5)：1623-1632.

[6] YAN Xiaoyuan，OHARA T，AKIMOTO H.Bottom-up estimate of biomass burning in mainland China[J].Atmospheric Environment，2006，40(27)：5262-5273.

[7] 秦世廣，丁愛軍，王韜.歐亞大陸生物質燃燒氣團的輸送特征及對中國的影響[J].中國環境科學，2006，26(6)：641-645.

[8] STREETS D，YARBER K，WOO J H，et al.Biomass burning in Asia：annual and seasonal estimates and atmospheric emissions[J].Global Biogeochemical Cycles，2003，17(4)：1759-1768.

[9] 曹國良，鄭方成，王亞強.中國大陸生物質燃燒排放的TSP，PM10，PM2.5清單[J]. 過程工程學報,2004(增刊1):700-704.

[10] 陸炳，孔少飛，韓斌，等.2007年中國大陸地區生物質燃燒排放污染物清單[J].中國環境科學，2011，31(2)：186-194.

[11] 李興華，王書肖，郝吉明.民用生物質燃燒揮發性有機化合物排放特征[J].環境科學，2011，32(12)：3515-3521.

[12] YUAN Cheng，ENGLING G，HE Kebin，et al.The characteristics of Beijing aerosol during two distinct episodes：impacts of biomass burning and fireworks[J].Environmental Pollution，2014，185(4)：149-157.

[13] 朱彬，蘇繼鋒，韓志偉，等.秸稈焚燒導致南京及周邊地區一次嚴重空氣污染過程的分析[J].中國環境科學，2010，30(5)：585-592.

[14] 徐艷.廣東典型生物質燃燒及煙氣排放特性研究[D].廣州：華南理工大學，2012.

[15] HE Min，ZHENG Junyu，YIN Shasha，et al.Trends，temporal and spatial characteristics，and uncertainties in biomass burning emissions in the Pearl River Delta，China[J].Atmospheric Environment，2011，45(24)：4051-4059.

[16] 曹國良，張小曳，王丹，等.中國大陸生物質燃燒排放的污染物清單[J].中國環境科學，2005，25(4)：389-393.

[17] 王書肖,張楚瑩.中國秸稈露天焚燒大氣污染物排放時空分布[J].中國科技論文,2008,3(5):329-333.

[18] HAO Weimin，LIU M H.Spatial and temporal distribution of tropical biomass burning[J].Global Biogeochemical Cycles，1994，8(4)：495-503.

[19] 田賀忠，郝吉明，陸永琪，等.中國生物質燃燒排放SO2、NOX量的估算[J].環境科學學報，2002，22(2)：204-208.

[20] ANDREAE M O，MERLET P.Emission of trace gases and aerosols from biomass burning[J].Global Biogeochemical Cycles，2001，15(4)：955-966.

[21] CAO Guoliang，ZHANG Xiaoye，GONG Sunling，et al.Investigation on emission factors of particulate matter and gaseous pollutants from crop residue burning[J].Journal of Environmental Sciences，2008，20(1)：50-55.

[22] ZHANG Hefang，YE Xingnan，CHENG Tiantao，et al.A laboratory study of agricultural crop residue combustion in China：emission factors and emission inventory[J].Atmospheric Environment，2008，42(36)：8432-8441.

[23] WANG Shuxiao，WEI Wei，DU Li，et al.Characteristics of gaseous pollutants from biofuel-stoves in rural China[J].Atmospheric Environment，2009，43(27)：4148-4154.