孟麗麗 王怡

摘要：分析了天津市大氣污染物的時空特征，利用數(shù)據(jù)資料重點(diǎn)研究了風(fēng)速、相對溫度、相對濕度，以及降水等對空氣污染的影響，在此基礎(chǔ)上利用逐步回歸統(tǒng)計方法建立了不同季節(jié)的PM10的濃度預(yù)報模型，并對這3種預(yù)報模型進(jìn)行了比較，結(jié)果表明：夏季的獨(dú)立樣本的檢驗(yàn)結(jié)果明顯優(yōu)于冬季，和天津當(dāng)?shù)氐臍夂虮容^融洽。同時，該模型的建立大大縮減了日常監(jiān)測的工作量和運(yùn)算過程，提高了效率，為政府采取普遍有效的污染防治措施提供了依據(jù)。

關(guān)鍵詞：空氣質(zhì)量;PM10;預(yù)報模型;逐步回歸

中圖分類號：X511

文獻(xiàn)標(biāo)識碼：A

文章編號：1674-9944（2018）6-0069-04

1引言

天津是北方重要的港口城市，處于環(huán)渤海地區(qū)和東北亞的核心區(qū)域。但是隨著現(xiàn)代化的發(fā)展，天津的白然環(huán)境卻在日益惡化。尤其近幾十年來，汽車數(shù)量成倍增長，汽車尾氣加重了空氣污染。大氣中的主要污染物是大氣氣溶膠、二氧化硫、二氧化氮、一氧化碳、氯氣和光化學(xué)煙霧?？諝馕廴締栴}日益嚴(yán)峻，人們開始對空氣質(zhì)量進(jìn)行預(yù)報，以便采取有效的污染防治措施。所以研究分析空氣污染狀況，建立空氣污染預(yù)報方案是當(dāng)前研究的重要課題。本文的研究內(nèi)容包括：天津市污染氣象條件分析，研究了風(fēng)速、溫度、相對濕度、降水與空氣的關(guān)系及對空氣污染的影響。通過上述分析，選出與PM10 濃度相關(guān)性較強(qiáng)的因子作為統(tǒng)計預(yù)報方法的輸入因子，利用逐步回歸法精選出預(yù)報模型的輸入因子，通過分季節(jié)建立PM10濃度的線性回歸模型，該模型通過了檢驗(yàn)，在此基礎(chǔ)上利用該模型進(jìn)行預(yù)報。

2天津市大氣質(zhì)量影響因素分析

2.1天津市地理位置氣候概況

天津市地處華北平原東北部，環(huán)渤海灣的中心.東臨渤海，北依燕山，位于北緯38°34'～40°15'，東經(jīng)116°43'-118°04'之間，處于國際東八時區(qū)。天津市土地面積為11919170hm2，南北長189km，東西寬117km。天津地質(zhì)結(jié)構(gòu)多樣，大部分被新陳代謝沉積物覆蓋，地勢以平原和洼地為主，海拔由南到北降低，北部最高有低山丘陵，海拔1052m，東南部最低，海拔為3.5m。

天津位于中緯度歐亞大陸東側(cè)，主要受季風(fēng)環(huán)流影響.是東亞季風(fēng)盛行的地區(qū)，屬于大陸性氣候。四季分明，春季多風(fēng).干旱少雨;夏季炎熱，降雨集中;秋季冷暖適中;冬季寒冷干燥，降雪量少。天津的年平均氣溫在11.4～12.9℃，市區(qū)平均氣溫最高值是12.9℃.1月份最冷，平均氣溫在-5-- -3℃之間;7月份最熱，平均氣溫在26 --27℃。夏季季風(fēng)盛行，冬季春季風(fēng)速最大;夏季秋季風(fēng)速最小。年平均風(fēng)速在2～4m/s，大部分為西南風(fēng)。天津年平均無霜期為196--246d，最長無霜期為267d，最短無霜期為171d。四季中，冬季最長，有156--167d，夏天有87--103d;春季56～61d，秋季最短，僅有50～56d。天津年平均降水量為520～660mm，降水日數(shù)在63--70d之間。在地區(qū)分布上山地多于平原，沿海多于內(nèi)陸。在季節(jié)分布上，6、7、8月份的降水量占全年的75%。天津的日照時間比較長，年日照時長在2500--2900h。

2.2天津市空氣中污染物主要來源

天津市空中的污染物主要為可吸入顆粒物二氧化硫，其次是二氧化氮和一氧化碳等。目前天津市常規(guī)進(jìn)項(xiàng)的空氣檢測檢測的污染物包括：二氧化硫、二氧化氮、臭氧、一氧化碳、PM10、PM2.5等。

2.3天津市氣象條件對空氣污染的影響

2.3.1風(fēng)速與空氣污染

風(fēng)對于空氣中的污染物有傳遞和運(yùn)輸?shù)淖饔?，風(fēng)速在一定程度上影響著空氣中污染物的濃度。天津處于東亞季風(fēng)盛行的地區(qū)，受到季風(fēng)環(huán)流支配，全年各季主導(dǎo)風(fēng)向都受到了季風(fēng)的很大影響。天津夏季多南風(fēng)、偏南風(fēng)、東南風(fēng)，受東南偏南風(fēng)的主導(dǎo)，平均風(fēng)速在1.7m/s，10%的大風(fēng)風(fēng)速可達(dá)到3m/s;過渡季多北風(fēng)、偏北風(fēng)、南風(fēng)和偏南風(fēng)，受西南風(fēng)的主導(dǎo)，平均風(fēng)速在2.2m/s，10%的大風(fēng)風(fēng)速可達(dá)到5m/s;冬季多北風(fēng)、東風(fēng)、南風(fēng)、西南風(fēng)、受東風(fēng)主導(dǎo)，平均風(fēng)速在2.9m/s，10%的大風(fēng)風(fēng)速可達(dá)到5m/s。

2.3.2相對溫度與空氣污染

近地表層的空氣溫度對空氣污染物有顯著影響，空氣中的污染物與近地表層的空氣成呈負(fù)相關(guān)關(guān)系。空氣溫度越高相對的大氣對流活動越強(qiáng)，空氣層的結(jié)構(gòu)也就越不穩(wěn)定，這就有利于空氣中的污染物傳播、擴(kuò)散，空氣中的污染物濃度得到稀釋，空氣質(zhì)量就會變好。反之，空氣溫度越低，空氣中的污染物不能被稀釋，污染物在空中沉積，就容易造成空氣污染物。夏季的太陽輻射強(qiáng)，相應(yīng)的近地層空氣溫度高，大氣對流活動強(qiáng)，空氣質(zhì)量較好。冬季太陽輻射弱，近地層空氣溫度低，大氣對流活動弱，空氣中的污染物得不到稀釋，空氣質(zhì)量差。

2.3.3相對濕度與空氣污染

空氣濕度也會影響空氣質(zhì)量的好壞，天津市隨著城市的快速發(fā)展，城區(qū)面積擴(kuò)大，水域面積相對減少，導(dǎo)致氣溫上升，同時地面的公共設(shè)施，城市馬路，人行道，建筑，多為不透水層，并具有良好的排水設(shè)施.降雨后雨水迅速流失，地面干燥;天津市內(nèi)地表植被覆蓋少，自然蒸發(fā)量和蒸騰量少，是天津天氣比較干燥。天津的濕度類型屬于干燥型，年平均相對濕度在59%左右，冬季和春季多風(fēng)，降水量少往往是天津最干燥的季節(jié)，春季的相對濕度在44%-53%，冬季的相對濕度在55%-69%;夏季季是天津雨水比較集中的季節(jié)，空氣中的濕度比較大，所以夏季的相對濕度在67%-74%，秋季的相對濕度在62%-68%。

2.3.4降水與空氣污染

降水能夠有效地清除、沖刷空氣中的各種污染物，降低空氣中各種污染物的濃度，使空氣中的污染物得到自然凈化，對減少污染物和抑制揚(yáng)塵有很好的作用，按季節(jié)來看，夏季更能體現(xiàn)降水對空氣質(zhì)量的影響，夏季降水可使空氣中的可吸入顆粒物明顯降低，使空氣質(zhì)量好轉(zhuǎn);可是冬季這種影響卻很少體現(xiàn)。主要是因?yàn)榻邓谙募据^集中且量多，冬季降水次數(shù)少而且量少。

3天津市大氣PM10預(yù)報模型及其應(yīng)用

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

本文選用的主要資料包括環(huán)境監(jiān)測資料和氣象資料。環(huán)境監(jiān)測資料：2013年12月～2014年11月天津市監(jiān)測站的主要空氣污染PM10、SO₂、NO₂的每日濃度資料，每日AQI指數(shù);月度平均資料。2012年12月～2014年12月，201 5年4月，2015年8月天津市每日PM10濃度。

氣象資料：主要是2012.12～2015.11天津基準(zhǔn)站地面氣象逐日觀測數(shù)據(jù)（數(shù)據(jù)由web站點(diǎn)“Reliable Prognosis”提供）。

3.2逐步回歸法

在多元回歸方程中，選取合適的變量是十分重要的。逐步回歸法是一種選擇最優(yōu)變量子集的方法，可以從提供的可供挑選的變量中，根據(jù)一定的顯著性檢驗(yàn)，每一步只選人一個變量進(jìn)入方程。

逐步回歸法是一種考慮全部解釋變量對被解釋變量作用大小，顯著程度大小，然后從大到小的引入回歸方程，對被解釋變量作用小的變量可能被排除在回歸方程外。對已經(jīng)被引入回歸方程的解釋變量在引入新的解釋變量后也有可能失去原來的重要性而被從回歸方程中提出。引入一個變量和從方程中提出一個變量都是逐步回歸的一步.每進(jìn)行一步都要進(jìn)行顯著性檢驗(yàn)，以確?；貧w方程中只含有對被解釋變量顯著影響的解釋變量。

3.3預(yù)報因子的選取

使用逐步回歸法對冬季預(yù)報因子進(jìn)行精選，得到6個預(yù)報因子，對春季預(yù)報因子精選得到4個預(yù)報因子，對夏秋預(yù)報因子精選得到8個預(yù)報因子。

3.4 PM10統(tǒng)計預(yù)報模型的建立

3.4.1多元線性回歸方程建立

利用所選取的預(yù)報因子，將天津市的PM10日均濃度作為預(yù)報對象建立回歸分析的資料長度分別是：春季5個月139d，夏秋季17個月518d，冬季9個月329d，使用SPSS統(tǒng)計軟件分別建立冬季，春季，夏秋季的PM10濃度的線性模型。其模型如下：

冬季：

y=415.215-4.667X1-4.089X2-45.556X3-23.838X4+0.216X5-l.239X6

春季：

y=5491.69-3.124X1-6.899X2+5.83X3-2.708X4

夏秋：

y=1679.62+0.559X1-1.737X2-l.079X3-1.307X4-2.063X5-4.667X6+0.802X7+3.105X8

利用回歸方程的顯著性檢驗(yàn)的方法對上述方程進(jìn)行顯著性檢驗(yàn)，由F檢驗(yàn)臨界值表可以查出，F(xiàn)（0.01）（8120） =2.16，F(xiàn)（0.01）（4120）-3.47，冬季、春季、夏秋季的F檢驗(yàn)值分別是F（冬）=2.907、F（春）=3.48、F（夏秋）=2.52。冬季和夏秋季F檢驗(yàn)值均大于F（0.01）（ 8120）=2.16，春季的F檢驗(yàn)值大于F（0.01）（4120）-3.47。所以分季建立的PM10預(yù)報方程是顯著的。方程式可行的。

3.4.2模型對獨(dú)立樣本檢驗(yàn)

利用該預(yù)報模型，對于沒有參與模型建立的2014年12月，2015年4月，8月的這3個月作為獨(dú)立樣本做逐日預(yù)報。從表2中可以看出2014年12月31d的相對誤差平均值72.20μg/m?，相對誤差在61.95%，對PM10 的濃度有一定的預(yù)報能力。

表3通過對2015年4月30d的數(shù)據(jù)進(jìn)行預(yù)測，得到結(jié)果其相對誤差在，絕對誤差控制在56.95%。對PM10濃度處于144～188之間的天數(shù)，有較好的預(yù)報效果，對PM10濃度較小時反應(yīng)不靈敏，一定范圍內(nèi)可以有較好的預(yù)報PM10濃度。

表4對2015年8月31 d進(jìn)行預(yù)測，得到結(jié)果其相對誤差在，絕對誤差是18.09%。在8月3日出現(xiàn)一次較大誤差，其相對誤差為83.21%。其他天數(shù)誤差均小于50%，維持在很低的誤差水平上，說明此模型可以很好的預(yù)測夏秋PM10濃度。

通過對3個獨(dú)立月份的預(yù)測，得到的結(jié)果。分季節(jié)建立的PM10濃度預(yù)報模型，模型表現(xiàn)出季節(jié)的差異，從預(yù)測效果上來看：夏秋季最好、春季次之、冬季最次。其原因是冬季一方面由于大風(fēng)天氣帶來的揚(yáng)塵，另一方面冬季氣溫偏低不利于污染物擴(kuò)散。夏季影響污染物濃度的因素變化活躍，各種因素的影響是回歸模型包含的因子豐富，更易得到準(zhǔn)確的預(yù)測。在實(shí)際的防治中應(yīng)該在冬季加以重視，保障人體健康安全不受到威脅.