成都市人為源揮發(fā)性有機(jī)物排放清單及特征

周子航，鄧 也，陸成偉，田 紅，宋丹林，馮 淼，黃鳳霞

成都市環(huán)境保護(hù)科學(xué)研究院，四川 成都 610072

成都市人為源揮發(fā)性有機(jī)物排放清單及特征

周子航，鄧 也，陸成偉，田 紅，宋丹林，馮 淼，黃鳳霞

成都市環(huán)境保護(hù)科學(xué)研究院，四川 成都 610072

基于成都市實(shí)地調(diào)查和環(huán)境統(tǒng)計(jì)等活動(dòng)水平數(shù)據(jù)，采用排放因子法和計(jì)算模型等，編制了2014年成都市人為源VOCs排放清單，并完成了空間分配和不確定性分析。成都市人為源VOCs排放量為15.8×104t，其中化石燃料固定燃燒源、工藝過(guò)程源、溶劑使用源、移動(dòng)源、儲(chǔ)存運(yùn)輸源、其他源排放量分別為0.5×104、3.8×104、6.0×104、4.9×104、0.4×104、2.2×104t，溶劑使用源為最大人為排放源，其次是移動(dòng)源和工藝過(guò)程源。木材加工業(yè)為最大工業(yè)貢獻(xiàn)源，然后依次是醫(yī)藥制造業(yè)、非金屬礦物制品業(yè)、化學(xué)原料、化學(xué)制品制造業(yè)、汽車制造業(yè)等。成都市人為源82%的VOCs排放量分布于二、三圈層的工業(yè)園區(qū)，而中心城區(qū)主要為移動(dòng)源和建筑施工所貢獻(xiàn)，其排放分布已隨建成區(qū)聯(lián)片發(fā)展而形成整體。排放清單活動(dòng)水平數(shù)據(jù)可靠性較高，而排放因子存在一定不確定性。

成都市；揮發(fā)性有機(jī)物；排放清單；空間分配

Abstract:Based on survey data, environmental statistics and other relevant activity data, VOCs anthropogenic sources emission inventory of Chengdu in 2014 was compiled by adopting the emission factor method and computing model. Besides, spatial distribution and uncertainties of the anthropogenic sources emission inventory were discussed. The total VOCs emission from these sources in Chengdu was 15.8×104tons, in which the stationary combustion source of fossil fuels, industrial processes, solvent utilization, mobile source, storage and distribution of fuel and other sources contributed 0.5×104tons, 3.8×104tons, 6.0×104tons, 4.9×104tons, 0.4×104tons, 2.2×104tons, respectively. Solvent utilization was the largest contributor compared to the other sources in Chengdu, followed by mobile source and industrial processes. Meanwhile, wood processing, pharmaceutical manufacturing industry, nonmetal mineral products, chemical industry and automobile industry were the major emission sectors of the industrial source. Industrial park in the suburbs contributed 82% of the total VOCs. The spatial distribution of VOCs emission in central city was integrated due to the urban development, and it was mainly contributed by the mobile source and architecture wall painting. The activity data tended to be reliable, whereas emission factors had some uncertainties.

Keywords:Chengdu;volatile organic compounds(VOCs);emission inventory;spatial distribution

揮發(fā)性有機(jī)物(VOCs)是影響大氣區(qū)域性復(fù)合型污染的重要前體物和參與物，在對(duì)流層大氣物理化學(xué)過(guò)程中體現(xiàn)出重要的作用[1]。VOCs和氮氧化物在光照條件下發(fā)生光化學(xué)反應(yīng)，生成臭氧及其他光化學(xué)污染物[2]。同時(shí)，大氣中的二次有機(jī)氣溶膠(SOA)是人為源或天然源排放的VOCs在大氣中氧化而生成[3]。隨著成都市經(jīng)濟(jì)高速發(fā)展，近年成都市臭氧濃度呈上升趨勢(shì)[4]，且明顯受人為活動(dòng)等因素影響[5]。成都市交通居民混合區(qū)和工業(yè)區(qū)大氣中VOCs濃度較高[6]，由機(jī)動(dòng)車、工業(yè)等污染源排放大量VOCs[7]，呈現(xiàn)為大氣復(fù)合污染特征。國(guó)內(nèi)VOCs排放清單研究日益增多[8-10]，成都市VOCs的排放清單多為以區(qū)域或省為研究對(duì)象的研究成果[1,7，11-13]，且主要以全國(guó)或省為尺度“自上而下”[14]的方法獲取基礎(chǔ)數(shù)據(jù)，甚少針對(duì)成都市VOCs排放源的特征進(jìn)行研究。因此，本研究基于成都市環(huán)境統(tǒng)計(jì)數(shù)據(jù)、實(shí)地調(diào)查數(shù)據(jù)和其他部門資料，著重通過(guò)實(shí)地調(diào)研取得成都市VOCs排放源基礎(chǔ)數(shù)據(jù)，以自下而上為主，自上而下為輔，建立成都市2014年人為源VOCs排放清單，并完成時(shí)空分配和不確定性分析，作為重要基礎(chǔ)數(shù)據(jù)[15]，最終應(yīng)用于空氣質(zhì)量模型的預(yù)警預(yù)報(bào)與控制對(duì)策等研究中[16]。

1 材料與方法

1.1研究范圍與估算方法

圖1為研究范圍，以成都市全域，包括中心城區(qū)(錦江區(qū)、青羊區(qū)、金牛區(qū)、武侯區(qū)、成華區(qū))、二圈層(龍泉驛區(qū)、青白江區(qū)、溫江區(qū)、新都區(qū)、雙流縣、郫縣)、三圈層(蒲江縣、大邑縣、金堂縣、新津縣、都江堰市、彭州市、邛崍市、崇州市)。

圖1 研究范圍示意圖Fig.1 Research scope

根據(jù)環(huán)境保護(hù)部相關(guān)指南，研究將VOCs人為源分為化石燃料固定燃燒源、工藝過(guò)程源、移動(dòng)源、溶劑使用源、儲(chǔ)存運(yùn)輸源、其他源(包括生物質(zhì)燃燒源、廢棄物處理源、餐飲油煙排放源)。研究基準(zhǔn)年為2014年。估算方法以排放因子法為主，VOCs排放量估算公式：

Ei=EFi,k×Ai,k

(1)

式中：i為排放源；k為技術(shù)類型；Ei為i排放源的排放量；EF為排放系數(shù)；A為活動(dòng)水平。

道路移動(dòng)源VOCs排放量使用IVE機(jī)動(dòng)車排放模型進(jìn)行計(jì)算，IVE模型由國(guó)際可持續(xù)發(fā)展研究中心(ISSRC)和加州大學(xué)河邊分校(UCR)開發(fā)，適用于發(fā)展中國(guó)家[17]。IVE模型利用其內(nèi)嵌的基礎(chǔ)排放因子和相應(yīng)修正系數(shù)的計(jì)算，最終得出不同技術(shù)水平機(jī)動(dòng)車的修正后排放因子[18]。該公式為

(2)

式中：Qj為j車型的修正后排放因子，g/km；Bj為j車型的基礎(chǔ)排放因子，g/km；i為待修正的各類修正系數(shù)；j為車輛類型；Ki,j為j車型的第i種修正系數(shù)，無(wú)量綱。

部分企業(yè)涉及到化學(xué)溶劑儲(chǔ)存的，使用美國(guó)環(huán)保署(USEPA)開發(fā)的Tanks 4.0.9d模型進(jìn)行計(jì)算[19-20]。

1.2活動(dòng)水平與排放因子

活動(dòng)水平數(shù)據(jù)通過(guò)多種渠道獲取，其中通過(guò)實(shí)地調(diào)查和填表方式，共獲取1 600余家企業(yè)數(shù)據(jù)、800余家加油站和油庫(kù)數(shù)據(jù)，并從交管部門獲取包含324萬(wàn)余輛機(jī)動(dòng)車信息的數(shù)據(jù)庫(kù)。排放因子則主要基于環(huán)保部指南[21]和國(guó)內(nèi)外研究成果。

1.2.1 化石燃料固定燃燒源

化石燃料固定燃燒源是指利用化石燃料燃燒時(shí)產(chǎn)生熱量，為發(fā)電、工業(yè)生產(chǎn)和生活提供熱能和動(dòng)力的燃燒設(shè)備。本研究將化石燃料固定燃燒源分為電力及熱力生產(chǎn)和供應(yīng)業(yè)、制造業(yè)、民用部門3部分。前兩部分活動(dòng)水平數(shù)據(jù)包括燃料信息、鍋爐和窯爐信息、排污設(shè)施經(jīng)緯度、排放筒信息等，主要來(lái)源于環(huán)境統(tǒng)計(jì)數(shù)據(jù)、排污申報(bào)數(shù)據(jù)、安監(jiān)部門鍋爐使用企業(yè)數(shù)據(jù)、企業(yè)調(diào)查數(shù)據(jù)、統(tǒng)計(jì)數(shù)據(jù)。民用部門活動(dòng)水平數(shù)據(jù)包括燃料使用信息和鍋爐信息，主要來(lái)源于環(huán)境統(tǒng)計(jì)數(shù)據(jù)、實(shí)地調(diào)查企業(yè)數(shù)據(jù)、統(tǒng)計(jì)數(shù)據(jù)。排放因子基于環(huán)保部指南[21]，詳見表1。

表1 化石燃料固定燃燒源VOCs排放因子Table 1 Emission factors of VOCs for stationary combustion source of fossil fuels

注：“*”表示排放因子單位為g/m3，其他排放因子單位為g/kg。

1.2.2 工藝過(guò)程源

工藝過(guò)程源是指工業(yè)生產(chǎn)和加工過(guò)程中，以對(duì)工業(yè)原料進(jìn)行物理和化學(xué)轉(zhuǎn)化為目的的工業(yè)活動(dòng)。研究主要考慮鋼鐵制造、水泥制造、玻璃制造、磚瓦制造、農(nóng)副食品加工、食品制造、制酒、化學(xué)原料和化學(xué)制品制造、醫(yī)藥制造、橡膠和塑料制品制造、造紙、石油加工等行業(yè)?；顒?dòng)水平數(shù)據(jù)包括原輔料使用信息、排污設(shè)施經(jīng)緯度、排放筒信息、VOCs處理效率和儲(chǔ)罐信息等，主要來(lái)源于環(huán)境統(tǒng)計(jì)數(shù)據(jù)、排污申報(bào)數(shù)據(jù)、經(jīng)信部門提供的四川石化基地2014年運(yùn)行數(shù)據(jù)、四川石化基地環(huán)評(píng)報(bào)告、實(shí)地調(diào)查企業(yè)數(shù)據(jù)、統(tǒng)計(jì)數(shù)據(jù)。排放因子主要借鑒了環(huán)保部指南[21]和國(guó)內(nèi)外研究[11,22-24]，詳見表2。

表2 工藝過(guò)程源VOCs排放因子Table 2 Emission factors of VOCs for industrial processes g/kg

1.2.3 移動(dòng)源

移動(dòng)源是指由發(fā)動(dòng)機(jī)牽引、能夠移動(dòng)的各種客運(yùn)、貨運(yùn)交通設(shè)施和機(jī)械設(shè)備，包括道路移動(dòng)源和非道路移動(dòng)源。

道路移動(dòng)源劃分為輕型客車、大型客車、貨車、公交車、出租車、摩托車?；顒?dòng)水平數(shù)據(jù)包括機(jī)動(dòng)車保有量、排放標(biāo)準(zhǔn)、最大總重量、行駛工況、行駛特征、行駛里程、燃料信息和氣象地形參數(shù)等，主要來(lái)源于車輛管理部門數(shù)據(jù)庫(kù)、各車型道路行駛工況測(cè)試、主要道路車流量調(diào)查、氣象部門以及相關(guān)技術(shù)文檔[25-27]。將上述數(shù)據(jù)，作為輸入數(shù)據(jù)經(jīng)IVE模型計(jì)算后，得到各類道路移動(dòng)源的排放因子，詳見表3。

由表3可知，非道路移動(dòng)源包括工程機(jī)械、農(nóng)業(yè)機(jī)械(包括農(nóng)用運(yùn)輸車)、小型通用機(jī)械、柴油發(fā)電機(jī)組、船舶、鐵路內(nèi)燃機(jī)車。

表3 移動(dòng)源VOCs排放因子Table 3 Emission factors of VOCs for mobile source g/kg

注：“*”表示由本研究建立的IVE模型計(jì)算得到；飛機(jī)排放因子單位為kg/起降周期(起降周期指飛機(jī)1個(gè)完整降落與起飛周期，即1個(gè)LTO循環(huán))，為A320、A321、B747、B737等26種國(guó)內(nèi)民航主要機(jī)型VOCs排放因子均值。

本研究的活動(dòng)水平數(shù)據(jù)主要參考國(guó)內(nèi)研究[28]的獲取方法：由成都市農(nóng)業(yè)機(jī)械總動(dòng)力占四川省農(nóng)用機(jī)械總動(dòng)力的比重，得到成都市農(nóng)用柴油使用量，進(jìn)而估算農(nóng)業(yè)機(jī)械和農(nóng)用運(yùn)輸車柴油使用量；由成都市建筑施工面積占四川省的比例，得到成都市建筑業(yè)柴油消耗量。成都雙流國(guó)際機(jī)場(chǎng)官方網(wǎng)站每日公布詳細(xì)的航班進(jìn)出港信息，本研究使用Python語(yǔ)言，結(jié)合urllib2庫(kù)[29]，從其網(wǎng)站上獲取航班號(hào)、始發(fā)地、目的地等詳細(xì)航班信息，并建立起飛著陸循環(huán)(LTO)動(dòng)態(tài)排放模型，通過(guò)doParallel包[30-31]實(shí)現(xiàn)了并行化計(jì)算，以提高效率。LTO動(dòng)態(tài)排放模型開源自由傳播，可以訪問(wèn)網(wǎng)址為https://github.com/Chengwei-Lu/dynaLTO獲取。

1.2.4 溶劑使用源

溶劑使用源是指生產(chǎn)、使用有機(jī)溶劑的工業(yè)生產(chǎn)和生活部門?；顒?dòng)水平數(shù)據(jù)主要來(lái)源于實(shí)地調(diào)查企業(yè)數(shù)據(jù)和統(tǒng)計(jì)數(shù)據(jù)，其中工業(yè)企業(yè)的溶劑使用作為點(diǎn)源數(shù)據(jù)，建筑噴涂、日用品使用、干洗、汽車修理、醫(yī)院等作為面源數(shù)據(jù)。建筑涂料使用量由成都市統(tǒng)計(jì)年鑒中建筑竣工面積通過(guò)經(jīng)驗(yàn)換算系數(shù)計(jì)算后獲取[1]。干洗、汽車修理和醫(yī)院則通過(guò)抽樣調(diào)查獲取成都市不同區(qū)域的干洗溶劑用量、修補(bǔ)漆使用量、醫(yī)用溶劑使用量后，由經(jīng)驗(yàn)系數(shù)計(jì)算后獲取。排放因子主要借鑒了環(huán)保部指南[21]和國(guó)內(nèi)外研究[32-40]，其中汽車制造排放因子由本地?cái)?shù)家大型汽車制造企業(yè)的涂料使用數(shù)據(jù)經(jīng)過(guò)物料衡算得到[32-33]，詳見表4。

表4 溶劑使用源VOCs排放因子Table 4 Emission factors of VOCs for solvent utilization g/kg

注：“*”表示通過(guò)企業(yè)物料衡算；“**”表示乙醇以100%揮發(fā)計(jì)。

1.2.5 儲(chǔ)存運(yùn)輸源

儲(chǔ)存運(yùn)輸源是指揮發(fā)性油氣產(chǎn)品被收集、儲(chǔ)存、運(yùn)輸和銷售的過(guò)程。活動(dòng)水平數(shù)據(jù)包括加油站、油庫(kù)的燃料吞吐量、加油槍和裝卸裝置等油氣回收改造情況、儲(chǔ)罐結(jié)構(gòu)構(gòu)造等數(shù)據(jù)，數(shù)據(jù)來(lái)源于成都市境內(nèi)加油站和油庫(kù)的調(diào)查表格填報(bào)，以及經(jīng)信部門統(tǒng)計(jì)數(shù)據(jù)。排放因子主要參考國(guó)內(nèi)標(biāo)準(zhǔn)[41-42]、AP-42排放因子庫(kù)[43]和專家咨詢，詳見表5。

表5 儲(chǔ)存運(yùn)輸源VOCs排放因子

注：“*”表示為專家咨詢，“/”表示管道運(yùn)輸排放極少，可忽略；“—”表示不考慮其有控制技術(shù)下的排放因子。

1.2.6 其他源

表6為其他源VOCs排放因子。其他源包括生物質(zhì)燃燒源、廢棄物處理源、餐飲油煙源。生物質(zhì)鍋爐、秸稈露天燃燒等活動(dòng)水平數(shù)據(jù)來(lái)源于主管部門數(shù)據(jù)和實(shí)地調(diào)查；廢水、固體廢棄物處理數(shù)據(jù)從環(huán)境統(tǒng)計(jì)數(shù)據(jù)中獲??；社會(huì)餐飲數(shù)據(jù)則源于對(duì)700余家餐館的現(xiàn)場(chǎng)抽樣調(diào)查，并由經(jīng)驗(yàn)系數(shù)計(jì)算后獲取。排放因子主要參考環(huán)保部指南[21]。

表6 其他源VOCs排放因子Table 6 Emission factor of VOCs for other source g/kg

1.3時(shí)空分配方法

時(shí)間分配：工業(yè)企業(yè)月變化系數(shù)主要以企業(yè)調(diào)查數(shù)據(jù)為主，按企業(yè)產(chǎn)品產(chǎn)量的月變化和涂料、有機(jī)溶劑等使用量進(jìn)行分配，其他污染源則按月平均分配；周變化系數(shù)和小時(shí)變化系數(shù)使用城市大氣污染物排放清單編制技術(shù)手冊(cè)(賀克斌主編)中提供的周變化系數(shù)。

空間分配：工業(yè)企業(yè)、儲(chǔ)存運(yùn)輸源、廢水處理和固廢處理均作為點(diǎn)源，利用經(jīng)緯度信息定位；移動(dòng)源排放主要與車流量和路網(wǎng)信息有關(guān)，進(jìn)行空間分配時(shí)則基于實(shí)際道路車流量和路網(wǎng)；餐飲油煙等居民生活排放分布基于人口因素，使用中國(guó)科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心提供的2010年1 km分辨率人口柵格數(shù)據(jù)(RESDC)，并結(jié)合2014年成都市各區(qū)(市)縣人口進(jìn)行修正。最后利用ARCGIS建立1 km×1 km的網(wǎng)格，將各污染源VOCs排放量分配到網(wǎng)格中。

1.4不確定性分析方法

研究采用不確定性定量分析軟件工具AuvToolPro(analysis of variability and uncertainty tool progress)進(jìn)行定量評(píng)估分析，該軟件工具由華南理工鄭君瑜教授提供，該軟件工具已經(jīng)獲取國(guó)家計(jì)算機(jī)軟件著作版權(quán)登記(登記號(hào)：2009SR09716)，其以自展模擬和蒙特卡羅模擬為基礎(chǔ)，完整地在計(jì)算機(jī)平臺(tái)上實(shí)現(xiàn)基于數(shù)值模擬的排放源清單不確定性定量分析過(guò)程[16]。

2 結(jié)果與討論

2.1成都市人為源VOCs排放量分析

基于上述方法，成都市2014年人為源VOCs排放量為158 442 t。其中，化石燃料固定燃燒源排放量為5 114 t，占總排放量的3%；工藝過(guò)程源排放量為38 127 t，占總排放量的24%；溶劑使用源排放量為59 636 t，占總排放量的38%；移動(dòng)源排放量為49 105 t，占總排放量的31%；儲(chǔ)存運(yùn)輸源、其他源，分別占總排放量的3%、1%。詳見表7。

圖2列出了各污染源的主要排放貢獻(xiàn)子源。化石燃料固定燃燒源VOCs排放主要由燃煤使用造成，雖然電力生產(chǎn)和供應(yīng)的燃煤消耗量與制造業(yè)相差不大，但由于電力生產(chǎn)和供應(yīng)部門的鍋爐更先進(jìn)、燃燒效率更高，所以制造業(yè)占化石燃料固定燃燒源排放量的70%，而電力生產(chǎn)和供應(yīng)僅占8%。工藝過(guò)程源VOCs主要來(lái)源于醫(yī)藥制造業(yè)、非金屬礦物制品業(yè)和化學(xué)制品制造業(yè)，分別占工藝過(guò)程源排放量的36%、31%、15%，其中非金屬礦物制品以水泥和玻璃制造為主。溶劑使用源VOCs主要來(lái)源于木材加工業(yè)、建筑裝飾、汽車制造業(yè)、家具制造業(yè)，分別占溶劑使用源排放量的36%、29%、7%、6%，其中木材加工業(yè)主要為配制和使用脲醛樹脂等膠粘劑時(shí)造成的排放。移動(dòng)源以道路移動(dòng)源的排放為主，由于成都市輕型客車、摩托車保有量分別達(dá)到280、60余萬(wàn)輛，且主要以揮發(fā)性較強(qiáng)的汽油為燃料，因此兩者貢獻(xiàn)了移動(dòng)源79%的排放量。儲(chǔ)存運(yùn)輸源主要來(lái)源于加油站的汽油油氣排放和儲(chǔ)油庫(kù)發(fā)油過(guò)程中的汽油油氣排放，而柴油油氣排放量較小。餐飲油煙貢獻(xiàn)了其他源67%的VOCs排放量。

表7 成都市人為源VOCs排放量Table 7 VOCs emissions of different sectors in Chengdu t

圖2 成都市人為源VOCs排放分擔(dān)率Fig.2 Sector contributions of VOCs for different sources

2.2成都市人為源VOCs排放清單空間分布

成都市19個(gè)區(qū)(市)縣中，VOCs排放量較大的為雙流縣(25 746 t)、龍泉驛區(qū)(14 440 t)、青白江(11 362 t)。成都市中心城區(qū)VOCs排放量為27 832 t，占總排放量的18%；二圈層VOCs排放量為79 171 t，占總排放量的50%；三圈層VOCs排放量為51 439 t，占總排放量的32%。

圖3、圖4為人為源和移動(dòng)源VOCs的排放分布。

圖3 成都市人為源VOCs排放空間分布Fig.3 Spatial distribution of VOCs from anthropogenic sources in Chengdu

圖4 成都市移動(dòng)源VOCs排放空間分布Fig.4 Spatial distribution of VOCs from mobile source in Chengdu

成都市中心城區(qū)工業(yè)企業(yè)數(shù)量較少，但機(jī)動(dòng)車保有量和建筑施工面積較大，其VOCs排放分布與移動(dòng)源和建筑工地分布有關(guān)，主要位于建成區(qū)和主要道路，同時(shí)VOCs排放量分布相對(duì)均勻，排放分布已隨各城區(qū)連成一片。二、三圈層VOCs排放以工業(yè)源為主，同時(shí)移動(dòng)源和建筑裝飾的貢獻(xiàn)亦較大，其VOCs排放主要分布于建成區(qū)、工業(yè)園區(qū)以及主要道路。

圖5、圖6為成都市工業(yè)源VOCs排放空間分布和行業(yè)貢獻(xiàn)。成都市中心城區(qū)工業(yè)VOCs排放主要分布：高新區(qū)科技工業(yè)園，以醫(yī)藥制造企業(yè)為主；高新西區(qū)，以塑料泡沫制造企業(yè)為主；武侯區(qū)金花鞋城，以制鞋企業(yè)為主；錦江區(qū)，以印刷企業(yè)為主。

二圈層工業(yè)VOCs排放主要分布：郫縣靠近高新西區(qū)一側(cè)，以印刷企業(yè)、電子設(shè)備制造企業(yè)和金屬制品企業(yè)為主；溫江海峽兩岸科技園，以塑料制造業(yè)和人造板制造企業(yè)為主；雙流縣西南航空港經(jīng)濟(jì)開發(fā)區(qū)，以人造板制造企業(yè)為主；龍泉驛區(qū)成都市經(jīng)濟(jì)開發(fā)區(qū)，以汽車制造企業(yè)為主；新都區(qū)內(nèi)新都工業(yè)園區(qū)，以塑料制品、汽車制造企業(yè)為主；青白江區(qū)工業(yè)園區(qū)，以鋼鐵制造和玻璃制造企業(yè)為主。

三圈層工業(yè)VOCs排放主要分布：蒲江縣較為分散，以印刷業(yè)為主；大邑縣經(jīng)濟(jì)開發(fā)區(qū)、沙渠工業(yè)園，以人造板制造為主；金堂工業(yè)園區(qū)、成阿工業(yè)園區(qū)，以玻璃制造為主；新津縣工業(yè)集中發(fā)展區(qū)，以化學(xué)制品為主；都江堰市經(jīng)濟(jì)開發(fā)區(qū)，以水泥制造為主；彭州市工業(yè)開發(fā)區(qū)及西南部，以水泥制造為主；邛崍市工業(yè)集中發(fā)展區(qū)，以及臥龍鎮(zhèn)等，以人造板制造為主；崇州市較為分散，以家具制造為主。

2.3排放清單不確定性定量分析結(jié)果

不確定性定量分析結(jié)果列于表8中。主要定量分析了化石燃料固定燃燒源、部分工業(yè)過(guò)程源、部分溶劑使用源和移動(dòng)源的排放量不確定性?？傮w來(lái)看，化石燃料固定燃燒源中的制造業(yè)和民用源，以及大部分道路移動(dòng)源的不確定性較小，而工業(yè)過(guò)程源和溶劑使用源相對(duì)較大。出現(xiàn)這種結(jié)果是由于工業(yè)、企業(yè)生產(chǎn)工藝的多樣化、復(fù)雜化，成都市本地研究的較少，且國(guó)內(nèi)研究成果的排放差異較大，排放因子選擇不確定性較高所造成。

圖5 成都市工業(yè)源VOCs排放空間分布Fig.5 Spatial distribution of VOCs from industrial sources in Chengdu

圖6 成都市各區(qū)(市)縣工業(yè)VOCs排放分擔(dān)率Fig.6 Sector contributions of VOCs for different districts in Chengdu

排放源計(jì)算值/t模擬運(yùn)算平均值/t95%置信區(qū)間化石燃料-電力生產(chǎn)419694(-82,149)化石燃料-制造業(yè)35854875(-51,51)化石燃料-民用源1110718(-81,80)工業(yè)過(guò)程-造紙544695(-97,251)工業(yè)過(guò)程-化學(xué)原料藥1390211492(-86,137)溶劑使用-制鞋339742(-74,125)溶劑使用-汽車制造41197292(-85,161)機(jī)動(dòng)車-輕型客車2369223807(-36,51)機(jī)動(dòng)車-大型客車20491722(-32,37)機(jī)動(dòng)車-輕型貨車18102108(-52,83)機(jī)動(dòng)車-中型貨車689527(-41,57)機(jī)動(dòng)車-重型貨車25472107(-65,112)機(jī)動(dòng)車-摩托車1474211711(-47,72)機(jī)動(dòng)車-公交車14661821(-29,30)機(jī)動(dòng)車-出租車14411830(-60,116)

3 結(jié)論

1)2014年，成都市人為源VOCs排放量為15.8×104t，其中溶劑使用源貢獻(xiàn)了38%的排放量，為最大人為排放源，其次是移動(dòng)源的31%和工藝過(guò)程源的24%。

2)從工業(yè)源的排放分擔(dān)率來(lái)看，木材加工業(yè)為最大工業(yè)貢獻(xiàn)源，然后依次是醫(yī)藥制造業(yè)、非金屬礦物制品業(yè)、化學(xué)原料和化學(xué)制品制造業(yè)、汽車制造業(yè)。

3)成都市人為源82%的VOCs排放量分布于二、三圈層的工業(yè)園區(qū)。而中心城區(qū)主要為移動(dòng)源和建筑施工所貢獻(xiàn)，其排放分布已隨建成區(qū)聯(lián)片發(fā)展而形成整體。

4)活動(dòng)水平數(shù)據(jù)可靠性較高，而排放因子存在一定不確定性。編制指南和其他研究提供的排放因子與成都市實(shí)際情況存在差異，排放因子的本地化工作有待進(jìn)一步開展。

[ 1] 魏巍.中國(guó)人為源揮發(fā)性有機(jī)化合物的排放現(xiàn)狀及未來(lái)趨勢(shì)[D].北京:清華大學(xué),2009.

[ 2] 唐孝炎,張遠(yuǎn)航,邵敏,等.大氣環(huán)境化學(xué)[M].北京:高等教育出版社,2006.

[ 3] 郝吉明,呂子峰,楚碧武,等.大氣二次有機(jī)氣溶膠污染特征及模擬研究[M].北京:科學(xué)出版社,2015.

[ 4] 錢駿,廖瑞雪,鄧?yán)?等.成都市大氣中O3污染水平及來(lái)源分析[J].四川環(huán)境,2011,30(3):20-23.

QIAN Jun,LIAO Ruixue,DENG Liqun,et al.Analysis of pollution degree and source of O3in the atmosphere of urban Chengdu[J].Sichuan Environment,2011,30(3):20-23.

[ 5] 崔蕾,倪長(zhǎng)健,王超,等.成都市O3濃度的時(shí)間變化特征及相關(guān)因子分析[J].中國(guó)環(huán)境監(jiān)測(cè),2015,31(6):41-46.

CUI Lei,NI Changjian,WANG Chao,et al.Temporal variation of ozone concentration and impact factors in Chengdu[J].Environmental Monitoring of China,2015,31(6):41-46.

[ 6] 印紅玲,袁樺蔚,葉芝祥,等.成都市大氣中揮發(fā)性有機(jī)物的時(shí)空分布特征及臭氧生成潛勢(shì)研究[J].環(huán)境科學(xué)學(xué)報(bào),2014,35(2):386-393.

YIN Hongling,YUAN Huawei,YE Zhixiang,et al.The temporal and spatial distribution of VOCs and their OFP in the atmosphere of Chengdu[J].Acta Scientiae Circumstantiae,2014,35(2):386-393.

[ 7] 韓麗,王幸銳,何敏,等.四川省典型人為污染源VOCs排放清單及其對(duì)大氣環(huán)境的影響[J]. 環(huán)境科學(xué),2013,34(12):4 535-4 542.

HAN Li,WANG Xingrui,HE Min,et al.Inventory and environmental impact of VOCs emission from the typical anthropogenic sources in Sichuan province[J].Environmental Science,2013,34(12): 4 535-4 542.

[ 8] 葉賢滿,徐昶,洪盛茂,等.杭州市大氣污染物排放清單及特征[J].中國(guó)環(huán)境監(jiān)測(cè),2015,31(2):5-11.

YE Xianman,XU Chang,HONG Shengmao,et al.Emission inventory and characteristics of air pollutants in Hangzhou[J].Environmental Monitoring of China,2015,31(2):5-11.

[ 9] 吳洪杰,劉玲英,蔡慧華,等.珠江三角洲制鞋行業(yè)揮發(fā)性有機(jī)化合物排放系數(shù)研究[J].中國(guó)環(huán)境監(jiān)測(cè),2013,29(4):74-78.

WU Hongjie,LIU Lingying,CAI Huihua,et al.Study on the emission coefficient of VOCs from shoe-making industry in the Pearl River Delta region[J].Environmental Monitoring of China,2013,29(4):74-78.

[10] 劉靜,王靜,宋傳真,等.青島市港口船舶大氣污染排放清單的建立及應(yīng)用[J].中國(guó)環(huán)境監(jiān)測(cè),2011,27(3):50-53.

LIU Jing,WANG Jing,SONG Chuanzhen,et al.The Establishment and application of ship emissions inventory in Qingdao Port[J].Environmental Monitoring of China,2011,27(3):50-53.

[11] WEI W,WANG S X,CHATANI S,et al.Emission and speciation of non-methane volatile organic compounds from anthropogenic sources in china[J].Atmospheric Environment,2008,42(20):4 976-4 988.

[12] BO Y,CAI H,XIE S D.Spatial and temporal variation of historical anthropogenic NMVOCs emission inventories in China[J].Atmospheric Chemistry and Physics,2008,8(23):7 297-7 316.

[13] CAI H,XIE S D.Estimation of vehicular emission inventories in China from 1980 to 2005[J].Atmospheric Environment,2007,41(39):8 963-8 979.

[14] WIKIPEDIA.Top-down and bottom-up design[EB/OL].[2008-10-20].http://en.wikipedia.org/wiki/Top-down_and_bottom-up_design,2008-10-20.

[15] 劉娟.長(zhǎng)三角區(qū)域環(huán)境空氣質(zhì)量預(yù)測(cè)預(yù)警體系建設(shè)的思考[J].中國(guó)環(huán)境監(jiān)測(cè),2012,28(4):135-140.

LIU Juan.Discussions on the system of regional air quality forecasting and early warning in Yangtze Delta area[J].Environmental Monitoring of China,2012,28(4):135-140.

[16] 鄭君瑜,王水勝,黃志炯,等.區(qū)域高分辨率大氣排放源清單建立的技術(shù)方法與應(yīng)用[M].北京:科學(xué)出版社,2014.

[17] 姚志良,賀克斌,王歧東,等.IVE機(jī)動(dòng)車排放模型應(yīng)用研究[J].環(huán)境科學(xué),2006,27(10):1 928-1 933.

YAO Zhiliang,HE Kebin,WANG Qidong,et al.Application study of IVE vehicle emission model[J].Enviromental Science,2006,27(10):1 928-1 933.

[18] International Sustainable Systems Research Center.IVE Model Users Manual[EB/OL].[2008-05-01].http://www.issrc.org/ive/downloads/manuals/UsersManual.pdf.

[19] 李靖,王敏燕,張健,等.基于Tanks4.0.9d模型的石化儲(chǔ)罐VOCs排放定量方法研究[J].環(huán)境科學(xué),2013,34(12):4 718-4 723.

LI Jing,WANG Minyan,ZHANG Jian,et al.Study on the quantitative estimation method for VOCs emission from petrochemical storage tanks based on Tanks4.0.9d model[J]. Enviromental Science, 2013,34(12):4 718-4 723.

[20] Environmental Protection Agency.User’s Manual for TANKS4.0.9d[EB/OL].[1999-09-30].http://www3.epa.gov/ttn/chief/software/tanks/tank4man.pdf.

[21] 環(huán)境保護(hù)部.大氣揮發(fā)性有機(jī)物源排放清單編制技術(shù)指南[EB/OL].[2014-08-20].http://www.mep.gov.cn/gkml/hbb/bgg/201408/W020140828351293705457.pdf

[22] 陳穎.我國(guó)工業(yè)源VOCs行業(yè)排放特征及未來(lái)趨勢(shì)研究[D].廣州:華南理工大學(xué),2013.

[23] European Environment Agency.EMEP/EEA air pollutant emission inventory guidebook-2013[EB/OL].[2013-08-29].http://www.eea.europa.eu//publications/emep-eea-guidebook-2013.

[24] 王宇楠,葉代啟,林俊敏,等.漆包線行業(yè)揮發(fā)性有機(jī)物(VOCs)排放特征研究[J].中國(guó)環(huán)境科學(xué),2012,32(6):980-987.

WANG Yunnan,YE Daiqi,LIN Junmin,et al.A study on emission characteristics of volatile organic compounds(VOCs) from enameled wire industry[J].China Environmental Science,2012,32(6):980-987.

[25] International Sustainable Systems Research Center.Beijing,China[EB/OL].[2008-09-29].http://www.issrc.org/ive/downloads/data/Beijing.zip.

[26] International Sustainable Systems Research Center.Shanghai,China[EB/OL].[2005-06-20].http://www.issrc.org/ive/downloads/data/Shanghai.zip.

[27] International Sustainable Systems Research Center.Tianjing,China[EB/OL].[2008-09-30].http://www.issrc.org/ive/downloads/data/Tianjin.zip.

[28] 張禮俊,鄭君瑜,尹沙沙,等.珠江三角洲非道路移動(dòng)源排放清單開發(fā)[J].環(huán)境科學(xué),2010,31(4):886-891.

ZHANG Lijun,ZHENG Junyu,YIN Shasha,et al.Development of non-road mobile source emission inventory for the Pearl River Delta region[J].Environmental Science,2010,31(4):886-891.

[29] KONG Q,CAI Y,ZHU Q.The case study for the basic information service of job post resource based on web mining[C]//Computer Science & Service System (CSSS),2012 International Conference on. IEEE, 2012:498-501.

[30] Revolution Analytics and Steve Weston (2015).Doparallel:foreach parallel adaptor for the ‘parallel’ package.R package version 1.0.10[EB/OL].[2015-10-14].http://CRAN.R-project.org/package=doParallel.

[31] Steve Weston,Rich Calaway(2014).Getting started with doparallel and foreach[EB/OL].[2015-10-13].https://cran.r-project.org/web/packages/doParallel/vignettes/gettingstartedParallel.pdf.

[32] 高廣亮.汽車涂裝揮發(fā)性有機(jī)化合物排放量的計(jì)算與分析[J].汽車工藝與材料,2013(1):46-54.

GAO Guangliang.Calculation and analysis of volatiole organic compounds emissons for car coating[J].Automobile Technology & Material,2013(1):46-54.

[33] 王錫春,楊必暖,林鳴玉,等.最新汽車涂裝技術(shù)[M].北京:機(jī)械工業(yè)出版社,1999.

[34] Environmental Protection Agency.Emissions factors & AP-42,pressure sensitive tapes and labels[EB/OL].[2009-10-14].http://www3.epa.gov/ttnchie1/ap42/ch04/final/c4s02_2j.pdf.

[35] FU X,WANG S X,ZHAO B,et al.Emission inventory of primary pollutants and chemical speciation in 2010 for the Yangtze River Delta region, China[J].Atmospheric Environment,2013,70:39-50.

[36] 楊利嫻,黃萍,趙建國(guó),等.我國(guó)印刷業(yè)VOCs污染狀況與控制對(duì)策[J].包裝工程,2012,33(3):125-131.

YANG Lixian,HUANG Pin,ZHAO Jianguo,et al.Current status and control strategies of VOC emission from printing industry[J].Packaging Engineering,2012,33(3):125-131.

[37] 葉代啟.全國(guó)VOCs污染控制指導(dǎo)培訓(xùn)手冊(cè)[EB/OL].[2010-10-21].http://www.docin.com/p-115802014.html.

[38] 臺(tái)灣環(huán)境保護(hù)署.2007.空氣污染物排放量推估手冊(cè)[EB/OL].[2010-02-06].http://teds.epa.gov.tw/new_main2.htm.

[39] 魏巍,王書肖,郝吉明.中國(guó)涂料應(yīng)用過(guò)程揮發(fā)性有機(jī)物的排放計(jì)算及未來(lái)發(fā)展趨勢(shì)預(yù)測(cè)[J].環(huán)境科學(xué),2009,30(10):2 809-2 815.

WEI Wei,WANG Shuxiao,HAO Jiming.Estimation and forecast of volatile organic compounds emitted from paint uses in China[J].Environmental Science,2009,30(10):2 809-2 815.

[40] 陳穎,葉代啟,劉秀珍,等.我國(guó)工業(yè)源VOCs排放的源頭追蹤和行業(yè)特征研究[J].中國(guó)環(huán)境科學(xué),2012,32(1):48-55.

CHEN Ying,YE Daiqi,LIU Xiuzhen.Source tracing and characteristics of industrial VOCs emissions in China[J].China Environmental Science,2012,32(1):48-55.

[41] 國(guó)家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局.散裝液態(tài)石油產(chǎn)品損耗：GB 11085—89[S].北京:中國(guó)標(biāo)準(zhǔn)出版社,1989.

[42] 國(guó)家石油和化學(xué)工業(yè)局.石油庫(kù)節(jié)能設(shè)計(jì)導(dǎo)則：SH/T 3002—2000[S].北京:中國(guó)石化出版社，2000.

EmissionInventoryandCharacteristicsofVolatileOrganicCompoundsfromAnthropogenicSourcesinChengdu

ZHOU Zihang, DENG Ye, LU Chengwei, TIAN Hong, SONG Danlin, FENG Miao, HUANG Fengxia

Chengdu Academy of Environmental Sciences, Chengdu 610072, China

X823

A

1002-6002(2017)03- 0039- 10

10.19316/j.issn.1002-6002.2017.03.07

2016-04-28;

2016-05-06

成環(huán)科研ky2014第010號(hào)

周子航(1985-)，男，重慶人，學(xué)士，工程師。