摘要：為掌握區(qū)域大氣污染物排放特征，系統(tǒng)分析2023年臨沂市人為源大氣污染物排放數(shù)據(jù)。采用自下而上的方法編制排放清單，從時空維度進(jìn)行分布特征研究。結(jié)果表明，工業(yè)過程源和化石燃料燃燒是污染物的重要來源。其中，工業(yè)過程源貢獻(xiàn)可吸入顆粒物（PM10）、細(xì)顆粒物（PM2.5）和揮發(fā)性有機(jī)物（Volatile Organic Compounds，VOCs）排放量的最大份額，分別為273.15萬、193.86萬、116.83萬t；化石燃料燃燒在SO2和氮氧化物（NOx）排放中占主導(dǎo)地位，分別為78.32萬、129.87萬t。污染物排放呈現(xiàn)顯著的季節(jié)性變化，SO2和NOx在冬季達(dá)到年均水平的1.37倍，夏季降至年均的0.85倍；空間分布上具有明顯的區(qū)域集聚特征，高排放區(qū)主要集中在重工業(yè)城市和沿海地區(qū)。研究結(jié)果為臨沂市制定差異化的大氣污染防治策略提供科學(xué)依據(jù)，有助于改善區(qū)域空氣質(zhì)量。

關(guān)鍵詞：大氣污染物；人為源；排放清單；時空分布特征；臨沂市

中圖分類號：X781.5 文獻(xiàn)標(biāo)識碼：A 文章編號：1008-9500（2025）01-0-03

Emission Inventory and spatiotemporal distribution characteristics of anthropogenic air pollutants in Linyi City

Abstract： In order to understand the characteristics of regional air pollutant emissions， a systematic analysis of anthropogenic air pollutant emissions data in Linyi city in 2023 is conducted. The bottom-up approach is adopted to compile emission inventories， and study distribution characteristics from a spatiotemporal perspective. The results indicate that industrial process sources and fossil fuel combustion are important sources of pollutants. Among them， industrial process sources contribute the largest share of inhalable particulate matter （PM10）， fine particulate matter （PM2.5）， and Volatile Organic Compounds （VOCs） emissions， with emissions of 2.731 5 million， 1.938 6 million， and 1.168 3 million tons， respectively; fossil fuel combustion accounts for the dominant emissions of SO2 and nitrogen oxides （NOx）， with emissions of 0.783 2 million and 1.298 7 million tons， respectively. The emissions of pollutants show significant seasonal changes， with SO2 and NOx reaching 1.37 times the annual average level in winter and decreasing to 0.85 times the annual average level in summer; there are obvious regional agglomeration characteristics in spatial distribution， with high emission areas mainly concentrated in heavy industrial cities and coastal areas. The research results provide scientific basis for the development of differentiated air pollution prevention and control strategies in Linyi city， which will help improve regional air quality.

Keywords： air pollutants; anthropogenic source; emission inventory; spatiotemporal distribution characteristics; Linyi city

大氣污染已成為制約臨沂市經(jīng)濟(jì)社會可持續(xù)發(fā)展的重要因素之一，準(zhǔn)確掌握大氣污染物排放特征是制定有效防控措施的基礎(chǔ)。臨沂市位于我國東部地區(qū)，東臨黃海，重化工業(yè)發(fā)達(dá)，鋼鐵、石化、電力等高耗能產(chǎn)業(yè)占比較大，為大氣污染物排放的主要來源。下面以臨沂市2023年人為源大氣污染物排放數(shù)據(jù)為依據(jù)編制清單，并根據(jù)當(dāng)?shù)禺a(chǎn)業(yè)結(jié)構(gòu)、經(jīng)濟(jì)發(fā)展?fàn)顩r等因素分析其時空分布特征，探討影響排放特征的關(guān)鍵因素。研究結(jié)果將為臨沂市大氣污染防治提供科學(xué)依據(jù)，對改善區(qū)域空氣質(zhì)量具有重要意義。

1 數(shù)據(jù)與方法

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

排放源活動水平數(shù)據(jù)來自臨沂市歷年發(fā)布的統(tǒng)計年鑒、能源統(tǒng)計年報和各地統(tǒng)計公報，數(shù)據(jù)包括工業(yè)、交通、農(nóng)業(yè)等部門的能源消費量、工業(yè)產(chǎn)品產(chǎn)量、機(jī)動車保有量等關(guān)鍵信息[1]。采用本地化排放因子，并結(jié)合臨沂市實際情況進(jìn)行適當(dāng)調(diào)整。對于缺乏本地化排放因子的污染源，使用美國發(fā)布的排放因子數(shù)據(jù)[2]。

1.2 排放清單編制方法

采用自下而上的方法編制臨沂市人為源大氣污染物排放清單，按照污染源類別將排放源劃分為點源、面源和線源。排放清單涵蓋SO2、氮氧化物（NOx）、CO、可吸入顆粒物（PM10）、細(xì)顆粒物（PM2.5）、揮發(fā)性有機(jī)物（Volatile Organic Compounds，VOCs）等主要污染物，采用式（1）計算排放量。不同類型污染源采用不同的計算方法。一是工業(yè)點源。工業(yè)點源主要包括電力、鋼鐵、水泥等重點行業(yè)的固定污染源，采用式（2）計算排放量。二是面源排放。面源包括分散式供熱、餐飲油煙等，采用式（3）計算排放量。三是機(jī)動車線源。采用式（4）計算機(jī)動車尾氣排放量。

E=A×F×（1-η）（1）

Eij=Ci×Qj×t×（1-ηj）（2）

Es=∑[Ak×Fk×（1-ηk）]（3）

Ev=∑（Pi×Ti×Fi）（4）

式中：E為主要污染物排放量；A為活動水平；F為未控制排放因子；η為污染控制效率；Eij為第i個污染源第j種污染物的排放量，t/a；Ci為污染物濃度，mg/m3；Qj為廢氣排放量，m3/h；t為年運行時間，h；ηj為污染物去除率，%；Es為面源排放量，t/a；Ak為第k類面源的活動水平；Fk為相應(yīng)的排放因子；ηk為污染控制效率，%；Ev為機(jī)動車排放量，t/a；Pi為第i類車輛保有量，輛；Ti為年均行駛里程，km；Fi為單位行駛里程的排放因子，g/km。

1.3 時空分布特征分析方法

時空分布特征分析采用多種方法相結(jié)合的策略。時間分布分析基于月度和季度排放量數(shù)據(jù)，運用時間序列分解技術(shù)，將排放量序列分解為趨勢、季節(jié)和隨機(jī)成分，揭示排放的長期變化趨勢和周期性特征。同時，利用自相關(guān)函數(shù)和偏自相關(guān)函數(shù)分析排放量的時間依賴性，識別關(guān)鍵的時間尺度?？臻g分布分析先采用地理信息系統(tǒng)，將排放清單數(shù)據(jù)與行政區(qū)劃圖層疊加，生成污染物排放的空間分布圖。隨后，應(yīng)用空間自相關(guān)分析方法，計算全局莫蘭指數(shù)和局部統(tǒng)計量，識別排放的空間聚集模式和熱點區(qū)域[3]。為深入理解空間分布的影響因素，采用地理加權(quán)回歸模型，探究經(jīng)濟(jì)發(fā)展水平、產(chǎn)業(yè)結(jié)構(gòu)、能源消費等因素對排放空間格局的影響。結(jié)合多尺度空氣質(zhì)量模型，模擬污染物的輸送擴(kuò)散過程，評估排放源對周邊區(qū)域空氣質(zhì)量的影響。

2 臨沂市人為源大氣污染物排放清單

2.1 主要污染物種類及排放量

基于排放清單編制結(jié)果，臨沂市人為源大氣污染物主要包括SO2、NOx、CO、PM10、PM2.5和VOCs。采集原始數(shù)據(jù)后，進(jìn)行嚴(yán)格的數(shù)據(jù)清洗。采用箱線圖法識別和處理異常值。缺失值使用多重插補(bǔ)法進(jìn)行填補(bǔ)，該方法考慮數(shù)據(jù)的不確定性，能更好地保持?jǐn)?shù)據(jù)的統(tǒng)計特性。在插值過程中，保留95%的置信區(qū)間內(nèi)的有效數(shù)據(jù)點，主要污染物種類及排放量如表1所示。

2.2 各行業(yè)部門排放貢獻(xiàn)

工業(yè)過程源是多種污染物的主要來源，尤其在PM10、PM2.5和VOCs的排放中占據(jù)主導(dǎo)地位，分別貢獻(xiàn)273.15萬、193.86萬、116.83萬t，顯示工業(yè)污染控制仍是重中之重?；剂先紵赟O2和NOx排放中貢獻(xiàn)最大，分別為78.32萬、129.87萬t，反映出能源結(jié)構(gòu)調(diào)整的重要性[4]。道路移動源對NOx和CO的排放貢獻(xiàn)顯著，分別為51.86萬、191.23萬t，突出交通污染治理的必要性。農(nóng)業(yè)源是NH3排放的絕對主力，達(dá)98.26萬t，說明農(nóng)業(yè)面源污染不容忽視。從污染物總量來看，CO排放量最高，達(dá)754.13萬t，其次是PM10（405.25萬t）和PM2.5（246.74萬t）。VOCs和NOx排放量相當(dāng)，分別為194.82萬、204.43萬t。

SO2排放量較低，為107.97萬t，可能反映出近年來脫硫措施的成效。

3 排放時空分布特征

3.1 時間分布

為分析臨沂市人為源大氣污染物排放的時間分布特征，采用基于活動水平的時間分配方法。收集2023年臨沂市各行業(yè)月度生產(chǎn)數(shù)據(jù)、能源消費數(shù)據(jù)、交通流量數(shù)據(jù)以及氣象數(shù)據(jù)等[5]，采用相應(yīng)方法計算時間分配系數(shù)。一是工業(yè)源。利用月度工業(yè)增加值和能源消費數(shù)據(jù)，計算月度分配系數(shù)。二是電力行業(yè)?；谠露劝l(fā)電量數(shù)據(jù)進(jìn)行時間分配。三是交通源。根據(jù)交通部門提供的月度和小時交通流量數(shù)據(jù)，計算道路移動源的時間分配系數(shù)。四是農(nóng)業(yè)源。利用農(nóng)業(yè)活動季節(jié)性特征和月度氣象數(shù)據(jù)，構(gòu)建NH3排放的月度分配模型。五是采暖排放。基于日均溫度數(shù)據(jù)，計算采暖天數(shù)，并結(jié)合能源消費數(shù)據(jù)，確定采暖期排放的時間分布。日變化和周變化分析采用典型行業(yè)和污染源的連續(xù)監(jiān)測數(shù)據(jù)，結(jié)合活動水平數(shù)據(jù)進(jìn)行模擬計算。

臨沂市大氣污染物排放呈現(xiàn)顯著的時間分布特征。季節(jié)性變化方面，SO2和NOx的排放量在冬季（12月至次年2月）達(dá)到全年最高水平，其中1月的月排放量約為月均排放量的1.37倍，而夏季（6—8月）則為全年最低，約為月均值的0.85倍。這主要歸因于冬季采暖需求增加及大氣擴(kuò)散條件較差。PM2.5的季節(jié)波動較小，但仍呈現(xiàn)冬季較高、夏季較低的趨勢。NH3的排放量在春季（3—5月）達(dá)到峰值，約為月均值的1.28倍，反映春耕施肥等農(nóng)業(yè)活動的影響。

3.2 空間分布

臨沂市人為源大氣污染物排放的空間分布特征分析采用規(guī)格1 km×1 km的網(wǎng)格系統(tǒng)，結(jié)合地理信息系統(tǒng)和空間分配模型實現(xiàn)。人為源大氣污染物排放呈現(xiàn)出區(qū)域差異與集聚特征。SO2和NOx排放高值區(qū)集中在濟(jì)南市、淄博市和濟(jì)寧市等重工業(yè)城市及沿海地區(qū)，原因是煤炭消費量大和高耗能產(chǎn)業(yè)密集。PM2.5和PM10排放分布與SO2和NOx相似。VOCs排放高值區(qū)位于東營市、濱州市等石化產(chǎn)業(yè)集中的地區(qū)，集中于濟(jì)南市、青島市等地的工業(yè)園區(qū)。NH3排放高值區(qū)集中在農(nóng)業(yè)發(fā)達(dá)的德州市、聊城市等地的農(nóng)村地區(qū)。由于產(chǎn)業(yè)結(jié)構(gòu)合理和氣候條件適宜，青島市、煙臺市等沿海城市排放較低。

4 結(jié)論

編制臨沂市人為源大氣污染物排放清單，可以揭示主要污染物的排放特征及其時空分布規(guī)律。研究發(fā)現(xiàn)，工業(yè)過程源和化石燃料燃燒是污染物的重要來源，排放呈現(xiàn)顯著的季節(jié)性變化和空間集聚特征。地理探測分析表明，工業(yè)產(chǎn)值密度、機(jī)動車保有量等社會經(jīng)濟(jì)因素對污染物空間分布有重要影響。

參考文獻(xiàn)

1 崔金夢，叢春華，鄭 怡，等.2013—2021年臨沂大氣污染特征及氣象條件分析[J].氣象與環(huán)境學(xué)報，2024（3）：46-54.

2 王治非，王在峰，呂 晨，等.濟(jì)南市冬季大氣污染時空分布特征及潛在污染源區(qū)分析[J].臨沂科學(xué)，2023（4）：114-121.

3 李 波，王 冰，李 婷.大氣污染的成因及治理措施分析[J].清洗世界，2022（7）：143-145.

4 孫 剛，張 超，鐘麗華.探析大氣污染環(huán)境監(jiān)測技術(shù)及治理方案[J].低碳世界，2021（11）：23-24.

5 李 樂，劉旻霞，肖仕銳，等.臨沂半島近地面O3濃度時空變化及潛在源區(qū)解析[J].環(huán)境科學(xué)，2022（3）：1256-1267.