垃圾填埋場惡臭污染統(tǒng)計分析及預(yù)報方程建立

吳 洋，趙 玥，孫樹鵬，廖云琛

(天津市津南區(qū)氣象局 天津 300350)

0 引 言

近年來，隨著我國城市化進程的加快，城市人口不斷增加，生活垃圾產(chǎn)量劇增。與此同時，我國城市生活垃圾分類處理體系尚未健全，填埋方式粗放，垃圾散發(fā)出的惡臭氣體嚴重危害周邊居民身心健康，已成為政府和社會公眾高度關(guān)注的問題[1]。國內(nèi)外有大量學(xué)者針對垃圾填埋場惡臭氣體成分等方面開展研究討論，指出硫化氫、氨氣、甲硫醇等是產(chǎn)生惡臭的主要氣體[2-3]。惡臭污染物主要是由填埋區(qū)的暴露作業(yè)面和填埋堆體產(chǎn)生，進而導(dǎo)致垃圾填埋場周邊惡臭污染嚴重[4]。還有學(xué)者研究發(fā)現(xiàn)惡臭氣體污染與氣象條件關(guān)系密切，氣象要素不僅影響惡臭污染物的生成，往往也影響著惡臭污染物的稀釋、擴散、輸送等過程。紀華等[5]的研究表明惡臭污染主要惡臭氣體硫化氫濃度具有明顯的季節(jié)變化，在夏季及秋季硫化氫濃度明顯升高，除與垃圾成分有關(guān)外，還與氣溫具有明顯關(guān)系。除氣溫外，在惡臭氣體輸送過程中，區(qū)域風場分布狀況也是影響空氣質(zhì)量的重要因素[6]。風在邊界層內(nèi)影響污染物的稀釋與擴散，風速有助于污染物水平輸送和擴散[7]，風向則決定著大氣污染物濃度的分布[8]。地面風場的變化很大程度上取決于大氣環(huán)流形勢的變化，因此，應(yīng)將環(huán)流場分布視為另一個影響垃圾填埋場周邊惡臭污染的氣象因素。李新令[9]和鄭秀美等[10]認為如果某地長期受移動緩慢的高壓控制就會造成該地風速小和穩(wěn)定層結(jié)的天氣特點，不利于污染物的稀釋和擴散；陳朝暉等[11]對一次重污染過程的大尺度天氣型分析發(fā)現(xiàn)，持續(xù)存在的均壓場是造成重污染濃度累積的主要背景場；鄧芳等[12]研究發(fā)現(xiàn)在實際的大氣環(huán)境中除了風速、風向等對惡臭污染起重要作用外，空氣濕度也影響污染物在大氣中的擴散；Guo H等[13]研究發(fā)現(xiàn)大氣穩(wěn)定度對惡臭污染具有影響，惡臭在穩(wěn)定的大氣中及低風速條件下能跨越更長更寬的距離。

國內(nèi)外許多學(xué)者針對氣象條件對垃圾填埋場惡臭污染物的影響進行了研究，多為各氣象要素對惡臭污染物影響的定性結(jié)論，針對惡臭污染物定量研究較少。本文以天津市津南區(qū)大韓莊垃圾填埋場為例，利用2019年5月14日至8月13日OU值和相應(yīng)時段的氣象要素地面觀測數(shù)據(jù)，經(jīng)相關(guān)分析和主成分分析，詳細闡述臭氣濃度OU值與各氣象要素的關(guān)系，篩選與惡臭污染物積累、擴散有關(guān)的氣象因子，建立惡臭污染物統(tǒng)計預(yù)報方程。

1 資料與方法

1.1 資料來源

本文所用數(shù)據(jù)為2020年5月14日至9月9日大韓莊垃圾填埋場OU值；同期津南區(qū)八里臺區(qū)域氣象觀測站觀測的氣壓、風向、風速、氣溫和相對濕度等氣象要素；美國環(huán)境預(yù)報中心(NCEP)提供的全球同期02時、08時、14時和20時(北京時間，下同)高空風場、位勢高度場和海平面氣壓場，格點分辨率為2.5°×2.5°；歐洲中心提供的全球同期02時、08時、14時和20時低云量及總云量，分別率為0.25°×0.25°。

1.2 分析方法

1.2.1 帕斯圭爾-特納爾方法

根據(jù)太陽高度角、云量確定輻射等級，然后再由輻射等級和風速將大氣穩(wěn)定度分為6個等級：1為極不穩(wěn)定、2為不穩(wěn)定、3為弱不穩(wěn)定、4為中性、5為弱穩(wěn)定、6為穩(wěn)定。

1.2.2 Lamb-Jenkinson分型方法

利用Lamb-Jenkinson 客觀分型方法[14-16](下文簡稱為L-J分型法)對津南區(qū)每日02時、08時、14時和20時海平面氣壓場進行客觀定量的環(huán)流分型。本文在每隔5個緯度、10個經(jīng)度的網(wǎng)格上取一個差分格點，共計16個格點。本文以(117°E，39°N)為中心點，選取100°—130°E、30°—50°N為研究范圍。

在該研究范圍內(nèi)，利用差分公式由16個網(wǎng)格點的日平均海平面氣壓值計算中心點的地轉(zhuǎn)渦度和地轉(zhuǎn)風，并根據(jù)地轉(zhuǎn)風和地轉(zhuǎn)渦度關(guān)系將大氣環(huán)流分為26類。

2 結(jié)果分析

2.1 OU值日變化

惡臭污染除受垃圾自身產(chǎn)生的惡臭氣體有關(guān)外，還與氣象要素密切聯(lián)系，氣象要素影響惡臭氣體的生成、輸送和擴散。氣溫、相對濕度等氣象要素具有明顯日變化，本文為探究OU值是否具有明顯日變化進行了研究。圖1為OU值日變化情況，可見津南區(qū)大韓莊垃圾填埋場OU值具有明顯日變化，呈明顯單峰型，峰值出現(xiàn)在每日00時至06時，07時OU值迅速下降，14時達到最低值，15時OU值迅速升高。

圖1 OU值日變化 Fig.1 Diurnal change of OU value

2.2 L-J環(huán)流分型分析

本文采用L-J分型方法對津南區(qū)2020年5月 14日至8月13日每日02時、08時、14時和20時大氣環(huán)流進行客觀分型，如圖2所示，W型、SW型、A 型和NW型等10類大氣環(huán)流型出現(xiàn)頻率超過3%，占所有環(huán)流型的83.2%。因此，通過對前10類大氣環(huán)流型進行研究可以較好地代表津南區(qū)該時段內(nèi)大氣環(huán)流。

圖2 2020年5月14日至8月13日津南區(qū)大氣環(huán)流型發(fā)生頻率 Fig.2 Frequency of atmospheric circulation types in Jinnan from May 14 to August 13，2020

為進一步研究前10類大氣環(huán)流對OU值的影響，分別統(tǒng)計02時、08時、14時和20時前10類大氣環(huán)流型出現(xiàn)時OU值的變化情況。由圖3可見，各大氣環(huán)流型均在02時OU值最高，14時OU值最低，08時和20時OU值介于14時和02時之間，與OU值日變化一致。02時SW型和A型OU值較高，14時ASW旋型和SE型OU值較高，08時和20時分別為AW型和ASW型OU值較高。

圖3 02時、08時、14時和20時各大氣環(huán)流型出現(xiàn)時OU值 Fig.3 OU value in different circulation types at 02:00,08:00, 14:00 and 20:00

不同的大尺度大氣環(huán)流對污染物擴散具有不同影響，為進一步研究大氣環(huán)流對惡臭氣體擴散的影響，對前10類大氣環(huán)流對惡臭氣體擴散的影響進行研究，篩選出不宜于臭氣擴散的大氣環(huán)流型。圖4給出了前10類大氣環(huán)流型下OU平均值，可見當出現(xiàn)W型、SW型、A型、C型、ASW和AW型6類大氣環(huán)流型時，OU值高于平均值，表示該6類大氣環(huán)流型不利于惡臭氣體擴散，易造成本地累積形成惡臭 污染。

圖4 各大氣環(huán)流型下OU平均值及OU值 Fig.4 Average OU values and OU values under different circulation types

2.3 相關(guān)性分析

國內(nèi)外學(xué)者針對氣象要素對垃圾填埋場惡臭污染物開展了很多研究，為進一步揭示氣象要素對津南區(qū)大韓莊垃圾填埋場惡臭污染物擴散的影響，統(tǒng)計了OU值與各氣象要素的相關(guān)性。

表1列出了大韓莊垃圾填埋場OU值與氣象要素相關(guān)性。氣溫、風速與OU值具有明顯負相關(guān)關(guān)系，OU值隨著氣溫升高、風速加大而減?。幌鄬穸燃按髿夥€(wěn)定度與OU值具有明顯正相關(guān)，表明OU值隨著相對濕度和大氣穩(wěn)定度的增加而增大。氣溫、相對濕度、風速和大氣穩(wěn)定度均通過0.001信度檢驗。風向和氣壓與OU值相關(guān)性不顯著；相對濕度與OU值相關(guān)系數(shù)最高(0.62)，大氣穩(wěn)定度和氣溫與OU值相關(guān)系數(shù)次之(分別為0.52和-0.48)，風向和氣壓與OU值相關(guān)系數(shù)相對較低(分別為0.08和0.09)。

表1 氣象要素與OU值相關(guān)系數(shù) Tab.1 Correlation coefficient between meteorological elements and OU valu e

2.4 建立統(tǒng)計預(yù)報方程

主成分分析方法是一種將多個指標化為少數(shù)幾個不相關(guān)的綜合指標的統(tǒng)計分析方法。它對于分析多指標的大量數(shù)據(jù)以了解數(shù)據(jù)間的關(guān)系及趨勢是一種很有效的方法。在進行主成分分析前需進行KMO檢驗和Bartlett球性檢驗。KMO檢驗是用于比較變量間簡單相關(guān)系數(shù)和偏相關(guān)系數(shù)的指標。通過計算KMO檢驗值為0.848，適合進行主成分分析；Bartlett球性檢驗結(jié)果為P＜0.001，說明數(shù)據(jù)呈球形分布，各個變量在一定程度上相關(guān)，適合進行主成分分析。

主成分分析按氣象要素和OU值進行分析，其特征值見表2。第一主成分的特征值為2.615，前4個主成分的累積貢獻率達79%以上，故選取前4個主成分。用表3中列出的特征向量值即可對選取的主成分進行解釋。

表2 相關(guān)系數(shù)矩陣特征值及貢獻率 Tab.2 Eigenvalues and contribution rates

表3 特征值對應(yīng)的特征向量 Tab.3 Eigenvectors corresponding to eigenvalues

第一主成分反映的是溫度、濕度的綜合，第二主成分代表風速指標，第三主成分代表大氣穩(wěn)定度指標，第四主成分代表氣壓指標。主成分分析表明溫度、濕度、風速、大氣穩(wěn)定度和氣壓對大韓莊垃圾填埋場OU值的影響比較大。通過上一節(jié)分析可知，氣溫、相對濕度、風速和大氣穩(wěn)定度與OU值具有明顯正負相關(guān)性。因此，綜合考慮選取了氣溫、相對濕度、風速和大氣穩(wěn)定度作為回歸方程自變量，OU值作為因變量，回歸方程如表4所示。

表4 統(tǒng)計預(yù)報方程 Tab.4 Statistical prediction equation

表中，1X～2X分別為氣溫、相對濕度、風速和大氣穩(wěn)定度。對各系數(shù)進行t檢驗，顯著性概率P值均小于0.001，說明各主成分對因變量OU值的影響作用均很顯著。

通過2020年5月14日至8月13日氣象要素驗證，回歸方程準確率為92.15%，從圖5中可以看出回代值可以較好地反映OU值變化趨勢，但是針對OU值極值模擬較差。為進一步驗證回歸方程預(yù)報準確性，本文對2021年8月14日至9月9日OU值進行預(yù)測，回歸方程預(yù)報準確率為85.42%，但是OU值的 預(yù)報值明顯偏高，如圖6所示。

圖5 2020年5月14日至8月13日OU觀測值和回代值 Fig.5 OU observation and regression from May 14 to August 13，2020

圖6 2020年8月14日至9月9日OU觀測值和預(yù)報值 Fig.6 OU observation and predicted value from August 14 to September 9，2020

3 結(jié) 論

①津南區(qū)大韓莊垃圾填埋場OU值具有明顯日變化規(guī)律，OU值呈明顯單峰型，峰值出現(xiàn)在每日 00時至06時。

②西風型、西南風型、反氣旋型和西北風型等 10類大氣環(huán)流型出現(xiàn)頻率占所有環(huán)流型的83.2%。

③當出現(xiàn)西風型、西南風型、反氣旋型、氣旋型、西南風反氣旋型和西風反氣旋型6類大氣環(huán)流型時不利于惡臭氣體擴散，易造成本地累積。

④氣溫、風速與OU值具有明顯負相關(guān)關(guān)系，相對濕度及大氣穩(wěn)定度與OU值具有明顯正相關(guān)關(guān)系。 ⑤通過主成分分析和相關(guān)性檢驗選取氣溫、相對濕度、風速和大氣穩(wěn)定度等氣象要素建立OU值的回歸方程，該方程預(yù)報準確率為85.42%，OU值的預(yù)報值明顯偏高。