曾明智,韓 波,莫海萍
(北海市環(huán)境監(jiān)測(cè)中心站,廣西 北海 536000)
空氣質(zhì)量指數(shù)計(jì)算方法在Excel中VBA的實(shí)現(xiàn)
曾明智,韓 波,莫海萍
(北海市環(huán)境監(jiān)測(cè)中心站,廣西 北海 536000)
空氣質(zhì)量指數(shù)計(jì)算方法是目前全國(guó)各城市空氣質(zhì)量重要信息發(fā)布的基礎(chǔ),廣泛應(yīng)用于城市空氣自動(dòng)監(jiān)測(cè)的空氣質(zhì)量評(píng)價(jià)。依據(jù)《HJ 633-2012環(huán)境空氣質(zhì)量指數(shù)(AQI)技術(shù)規(guī)定(試行)》發(fā)布環(huán)境空氣質(zhì)量指數(shù)(AQI)計(jì)算方法,利用VBA在Excel下編寫宏程序,可以輕松實(shí)現(xiàn)空氣質(zhì)量日?qǐng)?bào)信息的發(fā)布,使用戶從手工重復(fù)性的填充和插入函數(shù)操作方式中徹底解放出來(lái)。
空氣質(zhì)量指數(shù);計(jì)算工具;Excel;VBA
2013年1月以來(lái),我國(guó)京、津、冀及東部地區(qū)出現(xiàn)大范圍、長(zhǎng)時(shí)間霧霾天氣,一些大中城市空氣質(zhì)量明顯下降,給人民群眾生產(chǎn)生活帶來(lái)嚴(yán)重影響,引起了社會(huì)各界的高度關(guān)注。自2013年1月1日起,按照國(guó)務(wù)院批準(zhǔn)的空氣質(zhì)量新標(biāo)準(zhǔn)“三步走”實(shí)施方案,全國(guó)包括國(guó)家環(huán)保重點(diǎn)城市、模范城市在內(nèi)共116個(gè)城市449個(gè)監(jiān)測(cè)點(diǎn)位,已執(zhí)行《環(huán)境空氣質(zhì)量標(biāo)準(zhǔn)》 (GB 3095-2012),并按《環(huán)境空氣質(zhì)量指數(shù)(AQI)技術(shù)規(guī)定(試行)》 (HJ 633-2012)發(fā)布環(huán)境空氣質(zhì)量指數(shù)(AQI),不再按《城市空氣質(zhì)量日?qǐng)?bào)技術(shù)規(guī)定》 (總站辦字 [2000]26號(hào))發(fā)布空氣污染指數(shù)(API)。方案的實(shí)施對(duì)于及時(shí)按空氣質(zhì)量新標(biāo)準(zhǔn)監(jiān)測(cè)并發(fā)布信息,應(yīng)對(duì)重污染天氣,保護(hù)公眾身體健康等方面發(fā)揮了重要作用。
空氣質(zhì)量指數(shù)(AQI),與空氣污染指數(shù)(API)的定義和計(jì)算方法類似,首先收集各類主要污染物數(shù)據(jù),使用線性插值法計(jì)算每一類空氣質(zhì)量分指數(shù)(IAQI),比較各類IAQI大小,并將最大值作為向公眾發(fā)布的AQI。與API相比,
AQI更新了分級(jí)方案、濃度限值,并新增了評(píng)價(jià)項(xiàng)目PM2.5。因此,現(xiàn)行的API信息發(fā)布和計(jì)算方法已不能滿足當(dāng)前AQI計(jì)算的需要[1-5],最近有人提出用Excel計(jì)算AQI[6],該法僅限于Excel內(nèi)置條件函數(shù)、公式等進(jìn)行填充和插入函數(shù)的人工、半人工操作,這種計(jì)算模式,若在數(shù)據(jù)量龐大的情況下,容易出現(xiàn)錯(cuò)誤,降低了AQI計(jì)算方法的可行性、準(zhǔn)確性。為解決上述問(wèn)題,筆者應(yīng)用VBA語(yǔ)言編寫嵌入式模塊,開發(fā)出AQI自動(dòng)化的計(jì)算方法,通過(guò)實(shí)際應(yīng)用可以看出這種方法不僅對(duì)AQI計(jì)算效率和準(zhǔn)確度高,而且可移植到空氣污染指數(shù)(API)的計(jì)算工作中,具有較強(qiáng)的通用性和擴(kuò)展性。
1.1 空氣質(zhì)量指數(shù)分級(jí)方案
AQI是定量描述空氣質(zhì)量狀況的指數(shù),分為六級(jí),相對(duì)應(yīng)空氣質(zhì)量的6個(gè)類別,其數(shù)值越大說(shuō)明空氣污染狀況越嚴(yán)重,對(duì)人體健康的危害也就越大。參與空氣質(zhì)量指數(shù)計(jì)算的主要污染物為二氧化硫(SO2)、二氧化氮(NO2)、細(xì)顆粒物(粒徑小于等于10 μm)、顆粒物(粒徑小于等于2.5 μm)、一氧化碳(CO) 以及臭氧(O3),共計(jì)6個(gè)指標(biāo)。日?qǐng)?bào)AQI的IAQI級(jí)別及對(duì)應(yīng)的污染物項(xiàng)目濃度限值見表1。
表1 空氣質(zhì)量分指數(shù)及對(duì)應(yīng)的污染物項(xiàng)目濃度限值
1.2 空氣質(zhì)量分指數(shù)的計(jì)算方法
將監(jiān)測(cè)點(diǎn)各項(xiàng)污染物濃度與相應(yīng)的污染物濃度限值相比較,確定對(duì)應(yīng)于該濃度值時(shí)空氣質(zhì)量分指數(shù)(IAQI)所在的指數(shù)區(qū)間,再按照插值法計(jì)算該污染物濃度的IAQI值。
污染物項(xiàng)目P的空氣質(zhì)量分指數(shù)按式(1)計(jì)算:
式中:IAQIp——污染物項(xiàng)目P的空氣質(zhì)量分項(xiàng)指數(shù);Cp——污染物項(xiàng)目 的質(zhì)量濃度;BPHi——表1中Cp相近的污染物限值的高位值;BPL0——表1中Cp相近的污染物限值的低位值;IAQIHi——表1中與BPHi對(duì)應(yīng)的空氣質(zhì)量分指數(shù);IAQIL0——表1中與BPL0對(duì)應(yīng)的空氣質(zhì)量分指數(shù)。
1.3 空氣質(zhì)量指數(shù)級(jí)別
空氣質(zhì)量指數(shù)級(jí)別根據(jù)表2規(guī)定進(jìn)行劃分。
表2 空氣質(zhì)量指數(shù)、級(jí)別和類別劃分
1.4 空氣質(zhì)量指數(shù)及首要污染物的確定方法
1.4.1 空氣質(zhì)量指數(shù)計(jì)算方法
空氣質(zhì)量指數(shù)按式(2)計(jì)算:AQI=max{IAQI1,IAQI2,IAQI3,…,IAQIn} (2)式中:AQI——空氣質(zhì)量分指數(shù);n——污染物項(xiàng)目。
1.4.2 首要污染物及超標(biāo)污染物的確定方法
空氣質(zhì)量分指數(shù)大于50時(shí),最大的污染物為首要污染物,若最大的污染物為兩項(xiàng)或兩項(xiàng)以上時(shí),并列為首要污染物。
2.1 Excel空氣日?qǐng)?bào)數(shù)據(jù)表設(shè)計(jì)
新建一個(gè)Excel,按《HJ 633-2012》中表3的城市空氣日?qǐng)?bào)數(shù)據(jù)格式(見圖1所示的截圖)制作一個(gè)數(shù)據(jù)表,輸入到名為“sheet1”工作表中,將表1空氣質(zhì)量分指數(shù)及對(duì)應(yīng)的污染物項(xiàng)目濃度限值數(shù)據(jù)輸入到名為“sheet2”的工作表(見圖2所示的截圖),同時(shí)將該工作薄保存一個(gè)名為“AQI日?qǐng)?bào)”的Excel文件。
圖1 空氣質(zhì)量指數(shù)日?qǐng)?bào)數(shù)據(jù)格式
圖2 空氣質(zhì)量分指數(shù)及對(duì)應(yīng)的污染物項(xiàng)目濃度限值
2.2 程序設(shè)計(jì)
空氣質(zhì)量指數(shù)需要進(jìn)行大量的數(shù)據(jù)運(yùn)算,為了簡(jiǎn)化手工處理過(guò)程,筆者使用應(yīng)用程序開發(fā)語(yǔ)言VBA(Visual Basic for Application) 在Excel下編寫宏程序,自定義Excel工具欄、菜單欄和界面。簡(jiǎn)化模板的使用,使用戶在運(yùn)用空氣質(zhì)量指數(shù)計(jì)算方法中,直接點(diǎn)擊“日?qǐng)?bào)計(jì)算”按鈕,即可得到空氣質(zhì)量指數(shù)日?qǐng)?bào)結(jié)果??諝赓|(zhì)量指數(shù)計(jì)算方法VBA實(shí)現(xiàn)程序代碼清單如圖3所示的截圖。
圖3 空氣質(zhì)量指數(shù)計(jì)算方法程序清單
在進(jìn)行空氣質(zhì)量日?qǐng)?bào)信息發(fā)布或全年空氣質(zhì)量日?qǐng)?bào)統(tǒng)計(jì)時(shí),可以直接運(yùn)行上述程序即可得到參評(píng)的7項(xiàng)質(zhì)量分指數(shù)(IAQI)、空氣質(zhì)量指數(shù)(AQI)、主要污染物、空氣質(zhì)量指數(shù)級(jí)別和空氣質(zhì)量指數(shù)類別的計(jì)算結(jié)果。例如,計(jì)算北海市2013年9月份空氣日?qǐng)?bào)監(jiān)測(cè)結(jié)果時(shí),先將各指標(biāo)濃度數(shù)據(jù)輸入或?qū)氲綀D1的數(shù)據(jù)表中(見圖4所示的截圖),然后直接點(diǎn)擊“日?qǐng)?bào)計(jì)算”按鈕,便得到用戶所需要的結(jié)果(見圖5所示的截圖)。若需要日?qǐng)?bào)信息發(fā)布,用戶可將所需發(fā)布的日數(shù)據(jù)直接復(fù)制到Word文檔或另一個(gè)Excel數(shù)據(jù)表,然后上傳即可,不需要另外制作表格。
圖4 空氣質(zhì)量日?qǐng)?bào)監(jiān)測(cè)濃度數(shù)據(jù)輸入或?qū)虢Y(jié)果
圖5 空氣質(zhì)量日?qǐng)?bào)指數(shù)計(jì)算結(jié)果
與文獻(xiàn)[6]中的計(jì)算程序相比,本程序有如下特點(diǎn):
4.1 操作簡(jiǎn)單,運(yùn)算快速
實(shí)例的10條日?qǐng)?bào)數(shù)據(jù),只需單擊“日?qǐng)?bào)計(jì)算”按鍵,立即就能得到準(zhǔn)確的計(jì)算結(jié)果。
4.2 適于各月、年度空氣質(zhì)量日?qǐng)?bào)數(shù)據(jù)統(tǒng)計(jì)
程序不僅可以處理一條日?qǐng)?bào)數(shù)據(jù),而且可以處理數(shù)以萬(wàn)計(jì)的日?qǐng)?bào)數(shù)據(jù)(Excel2003及以前版本,每張工作表最多能保存65 536條日?qǐng)?bào)數(shù)據(jù)),只需花費(fèi)幾秒鐘就能完成,不需要修改程序或做任何操作,使得用戶統(tǒng)計(jì)每月、每年城市空氣質(zhì)量?jī)?yōu)劣天數(shù)和信息發(fā)布變得十分方便。
4.3 具有通用性與普及性
程序清晰,易掌握易改寫,不僅可計(jì)算各監(jiān)測(cè)點(diǎn)每日的AQI值,若將程序稍加改動(dòng),同樣能計(jì)算出HJ 633-2012中的表4空氣質(zhì)量指數(shù)實(shí)時(shí)報(bào)數(shù)據(jù)格式的發(fā)布內(nèi)容以及空氣污染指數(shù)(API)。
4.4 具有開拓性
城市空氣監(jiān)測(cè)數(shù)據(jù)往往按HJ 633-2012中的表3格式保存在Excel,在Excel中應(yīng)用VBA宏語(yǔ)言,開發(fā)協(xié)同性較好,開發(fā)成本低,簡(jiǎn)單易行,編程工作量較小,大大提高了空氣質(zhì)量信息發(fā)布的準(zhǔn)確性和使用效率,目前在空氣質(zhì)量日?qǐng)?bào)指數(shù)計(jì)算中尚為首次。
[1]黃戰(zhàn)勝,薛平.Excel軟件在計(jì)算空氣污染指數(shù)中的應(yīng)用[J].中國(guó)產(chǎn)業(yè),2011(2):40-41.
[2]李月娥,賀曉蕾,李昌平.利用Excel軟件計(jì)算空氣質(zhì)量日?qǐng)?bào)污染指數(shù)API[J].四川環(huán)境,2005,24(2):81-83.
[3]韓波,黃雄巍.神經(jīng)網(wǎng)絡(luò)模型在北海市空氣日?qǐng)?bào)污染指數(shù)計(jì)算中的應(yīng)用[J].中國(guó)環(huán)境管理干部學(xué)院學(xué)報(bào),2012,22(3):55-58,69.
[4]吳良忠,韓波,韋華.水質(zhì)灰色預(yù)測(cè)模型在VBA中的實(shí)現(xiàn)[J].中國(guó)環(huán)境管理干部學(xué)院學(xué)報(bào),2013,23(4):45-47.
[5]牟克林,黃世芹,吳江.基于EXCEL報(bào)表格式的城市空氣質(zhì)量預(yù)報(bào)評(píng)估系統(tǒng)[J].貴州氣象,2007,31(5):22-24.
[6]林藝濱.Excel軟件在計(jì)算空氣質(zhì)量指數(shù)的應(yīng)用[J].科技資訊,2012(15):131,133.
(編輯:周利海)
The Realization of Air Quality Index Calculation Method based on VBA in Excel
Zeng Mingzhi,Han Bo,Mo Haiping
(Beihai Environmental Monitoring Center,Beihai Guangxi 536000,China)
Air quality index calculation method is the basis of the release of important urban air quality information,which is widely used in air quality assessment in urban air automatic monitoring.According to the calculation method of air quality index(AQI)from the Technical Regulation on Ambient Air Quality Index(on trial,HJ 633-2012),the macro-program based on VBA in Excel was built to achieve the easy release of important daily urban air quality information,which can relieve the users from their repeated filling and inserting function.
air quality index;calculation tool;Excel;VBA
X823
A
1008-813X(2014)01-0008-04
2013-12-13
曾明智(1970-),男,廣西北海市人,畢業(yè)于天津科技大學(xué)工業(yè)分析專業(yè),工程師,主要從事環(huán)境監(jiān)測(cè)和科研工作。
10.13358 /j.issn.1008-813x.2014.01.03