周齊 彭亞綿 黃少騫 閆苗苗

【摘要】 隨著社會的發(fā)展，環(huán)境問題也日益嚴(yán)重，而大氣污染問題尤為突出. 本文首先通過歷史數(shù)據(jù)，建立空氣主要污染物關(guān)于時間的灰色預(yù)測模型，通過誤差分析，得到良好的預(yù)測結(jié)果，對未來空氣狀況做出預(yù)判；然后分析氣象對污染物濃度的影響，根據(jù)污染物濃度與污染物本身擴(kuò)散程度大致成反比的關(guān)系，得到氣象對空氣污染物擴(kuò)散的影響.

【關(guān)鍵詞】 灰色預(yù)測；誤差分析；氣象；擴(kuò)散程度；層次分析

【資助項目】 華北理工大學(xué)以升大學(xué)生創(chuàng)新實驗中心項目（Z2015017）

前 言

大氣問題目前日益突出，工程院院士鐘南山曾經(jīng)說過：中國大氣污染比沙斯更可怕，而我國作為工業(yè)大國，在經(jīng)濟(jì)高速發(fā)展的同時，大氣問題尤其突出，而建立污染物隨時間變化的函數(shù)模型，對于預(yù)測污染物濃度、預(yù)判大氣質(zhì)量以及采取相應(yīng)改良措施都具有極強(qiáng)必要性；而研究氣象對于污染物擴(kuò)散的影響對于了解自然、尋找與自然和諧相處之道具有重大價值.

一、濃度隨時間變化的灰色預(yù)測模型

查詢數(shù)據(jù)得到唐山市1—9月SO2濃度，并將其作為原始序列X（0），經(jīng)過一次累加弱化隨機(jī)性后得到一次累加序列X（1），如表1示：

按照回歸分析的思想，對總體X（1）進(jìn)行預(yù)測往往采用X（1）的條件數(shù)學(xué)期望值，又由：

+ af（t） = b

解得：f（t） = ce-at+

令t = ti時，f（t） = x（0）（ti），整理得：

X（1） = x（0）（ti）- e-a（t-1） + + ε

根據(jù)最小二乘法：

min （x（1）（ti） - f（ti））2

可得：

（1）x（0）（ti） - e +

代入數(shù)據(jù)可得：

x（1）（k） = x（0）（1） - e-0.110346（k-1） + （k=1，2，3，…）

進(jìn)行誤差分析，預(yù)測方程可用.

分別找到唐山市1—9月NO2濃度、PM10濃度、PM2.5濃度以及CO濃度通過同樣的灰色預(yù)測建模過程，分別得到灰色預(yù)測函數(shù)：

x（1）（k） = x（0）（1） - e-0.047453（k-1） +

x（1）（k） = x（0）（1） - e-0.113621（k-1） +

x（1）（k） = x（0）（1） - e-0.105408（k-1） +

二、氣象對擴(kuò)散影響的層次分析模型

首先通過查閱文獻(xiàn)得出風(fēng)、氣壓、氣溫和降水主要氣象指標(biāo)對污染物擴(kuò)散的影響過程，建立層次結(jié)構(gòu). 由于風(fēng)對大氣中污染物擴(kuò)散起決定性的作用，即風(fēng)在空氣中污染物濃度擴(kuò)散中的作用較大；云和氣壓有一部分原因是影響氣溫和降水，影響污染物擴(kuò)散. 由以上可以得到判斷矩陣，計算層次單排序的權(quán)向量和一致性檢驗成對比較矩陣A的最大特征根λ = 5.127. 進(jìn)行一致性檢驗，求解權(quán)重向量為：

ω = {0.541，0.223，0.121，0.070，0.045}

進(jìn)行一致性檢驗：

CI = = 0.032；RI = 1.12

表明通過了一致性檢驗.

【參考文獻(xiàn)】

[1]崔立志.灰色預(yù)測技術(shù)及其應(yīng)用研究[D].南京：南京航空航天大學(xué)，2010.

[2]楊華龍，劉金霞，鄭斌.灰色預(yù)測GM（1，1）模型的改進(jìn)及應(yīng)用[J].數(shù)學(xué)的實踐與認(rèn)識，2011（23）：39-46.

[3]張永波.基于灰色系統(tǒng)理論的預(yù)測模型的研究[D].哈爾濱：哈爾濱工程大學(xué)，2005.

[4]郭金玉，張忠彬，孫慶云.層次分析法的研究與應(yīng)用[J].中國安全科學(xué)學(xué)報，2008（05）：148-153.

[5]鄧雪，李家銘，曾浩健，陳俊羊，趙俊峰.層次分析法權(quán)重計算方法分析及其應(yīng)用研究[J].數(shù)學(xué)的實踐與認(rèn)識，2012（07）：93-100.

[6]蘭繼斌，徐揚，霍良安，劉家忠.模糊層次分析法權(quán)重研究[J].系統(tǒng)工程理論與實踐，2006（09）：107-112.