基于小波分析的廊坊市2015-2018年PM₂₅質量濃度研究

高清泉 紀詩璇 管林 劉伯華 范開宇 高清源

摘要 本文利用Morlet小波函數(shù)進行小波變換，以2015-2018年PM25日均質量濃度的時間序列為例，研究廊坊市大氣PM25質量濃度在時間序列上的周期特性。結果表明，廊坊市空氣PM2優(yōu)良等級天數(shù)保持連續(xù)增長，2015-2018年PM25年均質量濃度呈現(xiàn)逐年下降的良好趨勢。坊市PM2日均質量濃度具有顯著的周期性，以短期震蕩為主。主要周期為19d左右、35d左右以及60d左右，第一周期為19d?？蔀楹罄m(xù)廊坊市大氣污染物PM2s濃度變化周期性研究提供借鑒。

關鍵詞 小波變換;PM25;時間序列;周期性;廊坊市

中圖分類號：X513

文獻標識碼：B

文章編號：2095-3305（2021）02-116-04

大氣污染一直是擺在人類面前個不可忽視的問題，對人類健康和社會經濟發(fā)展具有嚴重的負面影響。我國北方大部分地區(qū)都遇到過具有區(qū)域性、復合性特點的比較嚴重的大氣污染過程叫，尤其是華北地區(qū)更為嚴峻。近年來，大氣污染現(xiàn)象雖然逐漸有所改善，但由于經濟、能源、交通等行業(yè)的發(fā)展，使得大氣污染的形勢依然比較嚴峻，受到政府與社會的持續(xù)關注。工業(yè)化和城市化造成大氣中可吸入顆粒物污染物増多，可吸入顆粒物極易被人體吸收，在呼吸系統(tǒng)累積而危害人體健康。PM2作為一種可吸人顆粒物，其空氣動力學當量直徑≤2.5um，相當于直徑小于人頭發(fā)絲的1/20。與其他較大的大氣顆粒物相比，PM23具有粒徑小、表面積大、強活性的特點，因其吸附能力強容易攜帶有毒物質（例如重金屬），且能夠在大氣遠距離輸送過程中滯留，可以經呼吸道到達肺泡，引起身體不適甚至多種疾病。PM2含量在地球大氣的占比微小，但是它對人民健康與大氣環(huán)境危害嚴重叫。有研究表明，PM2對人體的心血管系統(tǒng)、神經系統(tǒng)免疫系統(tǒng)有很大危害，大氣環(huán)境中PM2s濃度每天每增加10pg/m3，患腦血管疾病風險增加1.29%，因心血管疾病死亡的風險增加12%，此外相關研究表明，PM2與循環(huán)系統(tǒng)疾病死亡的相關性為正相關。因此，大氣環(huán)境中長期存在的PM2s被人體吸入，給身體健康帶來危害。

因此，對PM2s的社會關注度與相關研究也逐漸被重視。2012年國務院新修訂《環(huán)境空氣質量標準》，PM25作為新添加監(jiān)測指標，被納入污染指數(shù)項目中，并公布PM23污染等級劃分（表1）。許多專家學者也開展了針對PM2s相關研究。Zhao P S等研究了不同地區(qū)PM25質量濃度變化，宏觀分析了PM2質量濃度日變化、周變化及季節(jié)變化等污染特征及污染潛在源分析，主要采用統(tǒng)計及散點圖等直接分析的方法，此類方法在數(shù)據(jù)處理方式上相比較簡單，具有較好的可行性。但是面對大量數(shù)據(jù)，直接分析法存在顯示較為煩瑣、數(shù)據(jù)詳細變化可視化較低等不足。為更好地分析研究污染物濃度年際變化的基礎規(guī)律，小波分析因其可以精確地揭示時間序列中的頻率結構隨時間的變化規(guī)律叫，在大氣中可吸入顆粒物的時間序列特征分析應用中多有涉及?，此外，針對大中城市的PM2s的小波分析的研究取得了較好的研究進展。但是由于PM2的形成原因復雜，源頭多樣，傳輸途徑不穩(wěn)定，目前對PM23的溯源分析及影響因素未達成共識，要準確地了解其在不同時間尺度下的客觀規(guī)律性仍然很困難。對此，可以利用小波分析了解其在時間序列上的分布規(guī)律，從而進行時間序列上預測分析。研究發(fā)現(xiàn)，在時序規(guī)律變化中，PM2表現(xiàn)出較好的周期性，其質量濃度變化有顯著的周期規(guī)律和循環(huán)變化的特點。本文基于小波分析的基本原理利用小波函數(shù)，從小波變換角度出發(fā)，通過小波方差的方式，對廊坊市2015-2018年PM25質量濃度日平均數(shù)據(jù)進行時間序列上周期性分析與研究，獲取廊坊市PM2s污染的區(qū)域時空分布特點，增強認識其演變趨勢和影響因素，為廊坊市大氣污染物PM2質量濃度以后的小波分析研究提供依據(jù)，為廊坊市空氣污染監(jiān)控與防治提供參考。

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

1.1數(shù)據(jù)來源與處理

數(shù)據(jù)來源于廊坊市氣象信息共享平臺（htp：/10.48.36.193/），以廊坊市國家氣象觀測站內大氣成分站PM25質量濃度濃度日平均數(shù)據(jù)為數(shù)據(jù)源。數(shù)據(jù)起止時間：2015年1月1日至2018年12月31日。其中因大氣成分觀測設備故障停機造成的10d數(shù)據(jù)缺失，利用前后24h的逐時數(shù)據(jù)滑動平均的方法填補。對數(shù)據(jù)利用周期法進行邊界效應的消除。

1.2小波分析

1.2.1小波函數(shù)用小波函數(shù)表示或逼近某一信號或函數(shù)是小波分析的基本原理。小波分析過程的關鍵主要通過小波函數(shù)實現(xiàn)。小波函數(shù)是指類函數(shù)，具有震蕩性且能夠迅速衰減到0的特征。數(shù)學上，定義小波函數(shù)v（）∈Z（R），滿足：

滿足式（1）式的v（t）為基小波函數(shù)（或基本小波），基小波函數(shù)通過尺度的伸展、縮小和時間軸上的平移可以組成簇函數(shù)系列，見公式（2）：

其中，a，b∈R，a≠0.公式（2）中a（）為子小波函數(shù);a為尺度因子，可以反映小波的周期長度;b為平移因子反映在時間上的平移。

1.2.2小波変換小波變換是一種時域一頻域局部化方法，其時間窗和頻率窗都可以改變。大氣污染物日數(shù)據(jù)雖然不是原始的觀測數(shù)據(jù)，但是日數(shù)據(jù)也是種污染指標數(shù)據(jù)，可以看作是一種一維信號，適應于小波分析的基本原理，通過小波變換方程完成時間序列的特征研究。對于子小波函數(shù)yab（r），如果定義一個大氣污染物質量濃度有限時間序列可積函數(shù)f）∈Z（R），則の的連續(xù)小波變換為：

式（3）中，F(xiàn)（ab）為小波變換系數(shù);a為污染物質量濃度時間序列尺度因子;b為污染物質量濃度時間序列平移因子。

在大氣污染物實際的數(shù)據(jù)處理過程中，常見的時間序列是離散的。若定義時間序列函數(shù)f（k△1）（k=1，2，3，…，N;△為時間間隔），則公式（3）的離散函數(shù)為：

公式（4）中，小波變換系數(shù)Ⅳ（a，b）可以反映時域參數(shù)b與頻域參數(shù)a的特性。

Fa，b）受參數(shù)a和b的影響二變化。若以b為橫坐標，a為縱坐標，繪制Fa，b）的二維等值線圖，可得到小波變化系數(shù)圖。小波變換系數(shù)圖可以表現(xiàn)出在時間序列上的小波變化特征。若a固定，小波變換隨時間的變化過程可以發(fā)現(xiàn)系統(tǒng)在a尺度下的變換特征：小波變換系數(shù)絕對值越大，表明該時間尺度變化越顯著。因此，頻域-時域特性表現(xiàn)為：頻域參數(shù)變小，代表對頻域的分辮率低，引起對時域的分辨率變高;頻域參數(shù)變大，代表對頻域的分辨率高，從而能引起對時域的分辨率變低，小波變換能夠完成固定窗口大小、形狀可變的時-頻局部化處理。

小波函數(shù)的選擇是進行時間序列分析的關鍵。在實際處理數(shù)據(jù)時，小波類型的選取依靠自相似原則、判別函數(shù)以及支集長度等。本文選取Morlet小波變換的模和實部是信號序列分析的關鍵。模的大小表示特征時間尺度信號的強弱，實部表示不同特征時間尺度信號在不同時間上的分布和位相兩方面的信息，因此，Morlet小波比實數(shù)小波在數(shù)據(jù)處理時效果更佳，可以用作連續(xù)大氣污染物時間序列的周期分析1。

1.2.3小波方差為了判斷信號序列中的不同尺度變化的相對強度以及時間尺度，即主周期，可將小波系數(shù)的平方值在b域上的積分，得到小波方差，見式（5）：

Var（a）為小波方差，反映能量隨尺度a的分布，小波方差隨尺度a的變化過程，成為小波方差圖。因此，可以通過小波方差圖來査找并驗證一個時間序列中起主要作用的尺度。

2結果

2.1PM2s質量濃度時間序列變化特征

首先利用Morlet小波對2015-2018年PM2s日均質量濃度做時間序列變化處理。由于日均質量濃度數(shù)據(jù)為有限時間序列，為消除（或減?。?shù)據(jù)起點與終點附件可能會產生的“邊界效應”，先用信號延伸功能對兩端數(shù)據(jù)進行延伸處理，最后對“邊界效應”處理后的時間序列（圖1）進行小波變換。

PM2日均質量濃度在2015-2018年存在顯著的季節(jié)周期變化特征（圖1）。

其中，每年夏季PM5日均質量濃度最低，每年冬季均是一年期間污染最嚴重的季節(jié)，且存在峰值。PM2日均質量濃度污染天數(shù)呈逐年下降趨勢，根據(jù)PM污染等級劃分統(tǒng)計得知，2018年污染天數(shù)最少，2015年污染天數(shù)最為嚴重（表2）2015年PM2年平均濃度為85ug/mi2016年PM2年平均濃度為66pgm22017年PM2年平均濃度為60μgm2018年PM2年平均濃度為52pgm。

2015年嚴重污染天氣中PM2s質量濃度峰值，為490μg/m3，出現(xiàn)時間為2015年11月30日：2016年嚴重污染天氣中PM2質量濃度峰值，為393μg/m，出現(xiàn)時間為12月19日：2017年PM2s質量濃度峰值在12月13日，達到峰值358ugm;2018年的PM2s質量濃度峰值在11月26日大幅減小，為289ug/m。

2.2PM25質量濃度時間序列小波變換

根據(jù)公式（4），利用Morlet小波計算PM25質量濃度時間序列，得到小波系數(shù)實部等值線圖，即以b為橫坐標，a為縱坐標的時域一頻域分布圖（圖2）。

小波系數(shù)實部等值線圖既可以得到PM25質量濃度在時間域上的分布規(guī)律，也可以獲得PM2質量濃度時間序列上的周期性的變化特征，從而達到估測PM25質量濃度在不同時間尺度上的趨勢變化研究。等值線不同顏色表示不同大小的小波系數(shù)實部值，小波系數(shù)實部大于零，表示PM25質量濃度較大;波系數(shù)實部小于零，表示PM2質量濃度較小。小波系數(shù)實部近似于零時表示PM25質量濃度變化的突變點在時頻分布圖上為位相變換點，對分析PM2質量濃度變化規(guī)律具有很重要的意義。

在2015-2018年期間，PM2質量濃度值整體呈現(xiàn)大小交替的變化規(guī)律（圖2）。在不同時間尺度上的濃度變化結構也有所不同。在a=10～20、30～40、5060時，PM2s質量濃度的周期變化最為明顯。其中在a=10～20時，20152018年均存在顏色最深的等值線，說明有較強的周期變化，周期性明顯，且具有全域性，說明波動較大;a-3040時，2015、2016、2018年周期性比較強，2015與2016年尤為明顯，依然存在明顯的周期性。2017與2018年周期性較弱;當=5060時，PM2s質量濃度變化相對較小，等值線密度及顏色變化不明顯，周期性不強。

2.3PM2s質量濃度時間序列小波方差

為了可以更深入地探究PM2質量濃度的周期變化規(guī)律，利用繪制小波方差圖（圖3）的方法査找并驗證PM23質量濃度時間序列波動能量最大、起主要作用的時間尺度。

小波方差圖存在3個明顯的波峰（圖3），通過查看曲線波峰可以獲得時間序列的周期特性。即a=19、35和60，小波方差均達到極大值。小波方差圖總體呈現(xiàn)下降趨勢，在=60之后依然呈現(xiàn)下降的趨勢（圖3）。綜合小波系數(shù)實部等值線圖分析，=19時，PM2s質量濃度時間序列的波動能量最大，周期性最強，為第一周期;a=35時為第二周期，a=60時為第三周期。由此可以得知，2015年2018年PM2質量濃度時間序列變化規(guī)律可以通過第一周期、第二周期、第三周期完成周期性變化規(guī)律研究。

3結論

（1）郎坊市空氣PM2、優(yōu)良等級天

數(shù)保持連續(xù)增長，2015-2018年PM25年均質量濃度呈現(xiàn)逐年下降的良好趨勢。2015年PM2年平均濃度為85pgm，2016年PM2s年平均濃度為66μgm3，2017年PM2，年平均濃度為60pgm3，2018年PM2年平均濃度為52pgm。PM2日均質量濃度變化呈現(xiàn)U型變化?？諝釶M2污染主要分布在秋冬季（12月前后），2015年11月30日PM2s質量濃度高490μg/m，2018年的PM25質量濃度峰值在11月26日降至289gm。冬季污染較為嚴重，夏季污染程度相對減弱很多。

（2）利用Morlet小波對廊坊市20152018年PM2日均質量濃度時間序列分析，通過小波系數(shù)實部等值線圖與小波方差圖，研究在不同時間尺度下PM2s質量濃度變化的周期特性，結果發(fā)現(xiàn)，2015-2018年PM25質量濃度序列存在3個不同的周期，即19d左右、35d左右與60d左右，主要周期為19d左右。從而能說明廊坊市以PM2為代表的大氣污染物成一個月左右的周期變化。

（3）Morlet/小波能夠完成在不同時間尺度上對污染物周期性及局部極值變化的準確刻畫，相比較于其他分析方法具有科學而直觀優(yōu)勢。這對于更好地識別大氣污染物周期變化特征，建立基于小波分析時間序列預測預報模型，為環(huán)境監(jiān)測與預報具有十分重要的科學依據(jù)。

參考文獻

[1]任陣海，萬本太，蘇福慶，等.當前我國大氣環(huán)境質量的幾個特征[J]環(huán)境科學研究，2004，17（1）：1-6.

[2]曹明，林星春，劉剛.小波分解在空氣污染指數(shù)分析[J]中的應用中國環(huán)保產業(yè)，2006（10）：13-17.

[3]呂勇，付玉，辛士，等PM2s對人體的危害[J]遼寧化工.2017，46（60：618-620

[4]霍雨佳.“PM25”對人體健康的危害生命與災害，2011（12）：15

[5]郭新彪，魏紅英.大氣PM2對健康影響的研究進展[J]科學通報.2013，58（13）1171-1177

[6]董風鳴，莫運政，李國星，等，大氣顆粒物（PM1o/PM2）與人群循環(huán)系統(tǒng)亡系的病例交叉研究[J]北京大學學報：醫(yī)學版，2013，45（3）：398404.

[7] Pui D Y H， Chen S C， Zuo Z PM， in China： Measurements. sour and healtheffects， and mitigation[J]Particuology，2014，13（1）：1-26.

[8]顧吳元，肖翔，袁陳晨，等，基于小波神經網絡的松江區(qū)PM25濃度預測[J]上海工程技術大學學報，2015，29（2）：175-178

[9]毅，白靜文.承德市PM2成因，現(xiàn)狀及治理[J]對策科技資訊，2013（29）：14

[10] Zhao PS， Zhang X L， Xu X L， et al. Long-term visibility trends andcharacteristics inthe region of Beijing， Tianjin and HebeiChina[J]Atmospheric Research， 2011， 101（3）：711-718

[11]周偉基于MATLAB的小波分析應用[M].西安：西安電子科技大學出版社，2010

[12]陳柳，馬廣大，小波分析在PM濃度時間序列分析中的應用.監(jiān)測與評價，2006，24（1）：61-63

[13]楊書申，陳兵，邵龍義，北京市PM1o濃度變化規(guī)律的小波分析[J]中原工學院學報，2014，25（1）：39-43

[14]徐敬，丁國安，顏鵬，等，北京地區(qū)PM2s的成分特征及來源分析[J]應用氣象學報，2007，18（5）：645-654.

[15]楊鵬，劉杰，大氣污染物時空變化規(guī)律及其智能優(yōu)化算法研究[M].北京：科學出版社，2016.

[16]劉艷菊，張美根，張仁健，等，200年北京市春季大氣顆粒PM2和黑碳濃度變化特征[J]中國粉體技術，2010，16（1）18-22

[17]王文圣.Hydrology wavedet analysism北京：化學工業(yè)出版社，2005

[18]徐鳴，趙柳生，王斌PM1o濃度時間序列多時間尺度分析的小波方法[J]環(huán)境科學與技術，2008，31（4）：57-59

[19]吳小玲，張斌，艾南山，等，基于小波變換的上海市近10年SO2污染指數(shù)的變化[J]環(huán)境科學，2009，30（8）：21932198.

[20]張佩，聶義波，小波方差分析伊河流域降雨徑流周期特征[J]水資源與水工程學報，2009，20（2）：156

責任編輯：黃艷飛