多模式空氣質(zhì)量集成預(yù)報(bào)模型的研究

姚文強(qiáng)，王兆青，鐵治欣，丁成富，唐 敏

(1.浙江理工大學(xué)信息學(xué)院，杭州 310018；2.聚光科技(杭州)股份有限公司，杭州 310052)

?

多模式空氣質(zhì)量集成預(yù)報(bào)模型的研究

姚文強(qiáng)1，王兆青1，鐵治欣1，丁成富2，唐 敏1

(1.浙江理工大學(xué)信息學(xué)院，杭州 310018；2.聚光科技(杭州)股份有限公司，杭州 310052)

針對目前區(qū)域空氣質(zhì)量預(yù)報(bào)出現(xiàn)的預(yù)測不準(zhǔn)確、誤差較大、漏報(bào)重污染天氣等問題，采用WRF-Chem、WRF-CMAQ、MM5-CMAQ和MM5-CAMx 4種數(shù)值模式集成預(yù)報(bào)的方式,對區(qū)域空氣質(zhì)量狀況及污染物濃度進(jìn)行預(yù)測預(yù)報(bào)，提出了一個(gè)動(dòng)態(tài)權(quán)重更新模型，該模型通過評估一段時(shí)間內(nèi)與特定條件下單個(gè)模式的預(yù)測準(zhǔn)確度，動(dòng)態(tài)調(diào)整各個(gè)模式在集成預(yù)報(bào)模型中的權(quán)重因子，然后利用各個(gè)模式預(yù)測值及其權(quán)重取加權(quán)平均值作為新模式的預(yù)測值。以收集到的浙江省空氣質(zhì)量監(jiān)測網(wǎng)中的47個(gè)國家站點(diǎn)的8個(gè)月的監(jiān)測數(shù)據(jù)為測試數(shù)據(jù)，分別用已有模型和動(dòng)態(tài)權(quán)重更新模型進(jìn)行預(yù)測預(yù)報(bào)，并對預(yù)測預(yù)報(bào)結(jié)果進(jìn)行評估分析，結(jié)果表明：動(dòng)態(tài)權(quán)重分配模型預(yù)報(bào)結(jié)果具有更高的準(zhǔn)確度與更小的誤差。

空氣質(zhì)量；多模式；數(shù)值預(yù)報(bào)；動(dòng)態(tài)權(quán)重更新

0 引 言

近年來，隨著污染問題的日益突出，空氣質(zhì)量狀況已經(jīng)成為公眾越來越關(guān)心的話題，空氣質(zhì)量預(yù)測模型也成為科學(xué)研究的熱點(diǎn)。我國的大氣污染問題比較復(fù)雜，特別是最近幾年，已經(jīng)由過去單一類型的污染逐步轉(zhuǎn)變?yōu)槎喾N污染物、多種作用機(jī)制同時(shí)存在、相互影響的復(fù)合型大氣污染問題。復(fù)雜的污染問題必然會增加大氣質(zhì)量預(yù)測的難度，傳統(tǒng)的以單一類型污染為研究對象的預(yù)測方法已無法全面揭示當(dāng)前污染的特征和演變規(guī)律。所以改進(jìn)優(yōu)化空氣質(zhì)量預(yù)測方式，研究準(zhǔn)確有效的預(yù)測模型尤為重要。

現(xiàn)有的空氣質(zhì)量預(yù)測預(yù)報(bào)方法主要分為統(tǒng)計(jì)預(yù)報(bào)與數(shù)值預(yù)報(bào)。統(tǒng)計(jì)預(yù)報(bào)是根據(jù)統(tǒng)計(jì)學(xué)原理，用數(shù)理統(tǒng)計(jì)的方法對空氣質(zhì)量的變化趨勢進(jìn)行預(yù)測。數(shù)值預(yù)報(bào)是根據(jù)大氣實(shí)際情況，在一定的初值條件下，運(yùn)用計(jì)算機(jī)作數(shù)值計(jì)算，預(yù)測未來一定時(shí)間大氣污染物的濃度與變化趨勢。數(shù)值預(yù)報(bào)起源于國外，迄今已經(jīng)發(fā)展到第三代[1]。我國有關(guān)大氣數(shù)值預(yù)報(bào)的研究起步較晚，但最新一代的空氣質(zhì)量數(shù)值預(yù)報(bào)模式已經(jīng)在國內(nèi)空氣質(zhì)量預(yù)報(bào)及污染物預(yù)報(bào)領(lǐng)域得到了廣泛的應(yīng)用[2-7]。多模式集成預(yù)報(bào)日益成為數(shù)值預(yù)報(bào)的發(fā)展趨勢，國內(nèi)學(xué)者運(yùn)用基于多模式的空氣質(zhì)量集成預(yù)報(bào)技術(shù)，實(shí)現(xiàn)對當(dāng)前大氣復(fù)合污染的預(yù)測預(yù)報(bào)。王自發(fā)等[8]以NAQPMS模式、CMAQ模式[9]及CAMx模式為核心，采用統(tǒng)一的模式網(wǎng)格、統(tǒng)一的排放源、統(tǒng)一的氣象場，構(gòu)建了EMS-Beijing空氣質(zhì)量多模式集成預(yù)報(bào)系統(tǒng)，對于預(yù)測北京及周邊污染物起到了重要作用。王茜等[10]為上海市建立了包含NAQPMS、CMAQ、CAMx、WRF-Chem等多個(gè)國內(nèi)外著名數(shù)值模型的集合業(yè)務(wù)預(yù)報(bào)模式系統(tǒng)，為上海市提供了豐富的氣象和污染預(yù)報(bào)產(chǎn)品，經(jīng)過評估后發(fā)現(xiàn)集合預(yù)報(bào)系統(tǒng)對于PM10、SO2和NO2日均濃度的預(yù)測值與實(shí)測值相關(guān)系數(shù)達(dá)0.5～0.6，并且在2010年上海世博會開幕前夕及2010年5月19日兩次污染事件中，該集成預(yù)報(bào)系統(tǒng)為上海市環(huán)保局啟動(dòng)世博預(yù)警聯(lián)動(dòng)方案提供了重要的決策依據(jù)。謝濤等[11]將地理信息系統(tǒng)(GIS)技術(shù)與多模式集合預(yù)報(bào)系統(tǒng)結(jié)合的方式，在保障預(yù)測準(zhǔn)確度的前提下，充分發(fā)揮GIS在數(shù)據(jù)管理、信息可視化、輔助決策支持等方面的優(yōu)點(diǎn)，為廣州市在環(huán)境質(zhì)量預(yù)測方面提供了保障。黃思等[12]運(yùn)用多模式集合預(yù)報(bào)與多元線性回歸的方法，將數(shù)值預(yù)報(bào)與統(tǒng)計(jì)預(yù)報(bào)的技術(shù)相結(jié)合使用，對北京地區(qū)的PM10預(yù)報(bào)結(jié)果進(jìn)行修正，并且大幅提高了對污染過程的預(yù)報(bào)能力。

本文通過研究國內(nèi)外現(xiàn)有較好的數(shù)值預(yù)報(bào)模式，并且結(jié)合浙江省的氣候、地形、污染源排放、大氣污染等方面的特點(diǎn)，建立起包括WRF-Chem、WRF-CMAQ、MM5-CMAQ[13]和MM5-CAMx 4套模式組合在內(nèi)的多模式空氣質(zhì)量集成預(yù)報(bào)系統(tǒng)，并設(shè)計(jì)了一個(gè)基于預(yù)測準(zhǔn)確度的動(dòng)態(tài)權(quán)重更新模型，通過評估單個(gè)模式在一段時(shí)間與特定條件下的預(yù)測準(zhǔn)確度，來動(dòng)態(tài)調(diào)整各個(gè)模式在集成預(yù)報(bào)過程中的權(quán)重因子，以提高集成預(yù)報(bào)的準(zhǔn)確度。

1 多模式預(yù)報(bào)系統(tǒng)的總體框架設(shè)計(jì)

1.1 系統(tǒng)框架

浙江省大氣質(zhì)量多模式預(yù)報(bào)系統(tǒng)運(yùn)行在統(tǒng)一的計(jì)算平臺和運(yùn)行環(huán)境中，以WRF/CHEM、Models-3/CMAQ、CAMx為核心模式構(gòu)成WRF-Chem、WRF-CMAQ、MM5-CMAQ、MM5-CAMx預(yù)報(bào)組合，并基于4套模式組合的預(yù)報(bào)結(jié)果建立集成預(yù)報(bào)技術(shù)，整合不同模式的預(yù)報(bào)性能優(yōu)勢來實(shí)現(xiàn)對污染過程影響的分析和污染來源的量化。

在區(qū)域空氣質(zhì)量模式本地化的基礎(chǔ)上，采用模塊化設(shè)計(jì)、編程技術(shù)和網(wǎng)絡(luò)技術(shù)實(shí)現(xiàn)全球預(yù)報(bào)場的自動(dòng)下載、初始數(shù)據(jù)的自動(dòng)預(yù)處理、空氣質(zhì)量模式的自動(dòng)運(yùn)行、預(yù)報(bào)的結(jié)果的后處理及產(chǎn)品的自動(dòng)生成等功能。多模式預(yù)報(bào)系統(tǒng)的構(gòu)成框架如圖1所示。

圖1 浙江省多模式預(yù)報(bào)系統(tǒng)框架

1.2 模式模擬網(wǎng)格設(shè)置

水平模擬網(wǎng)格采用LambertConformal投影，投影參數(shù)為北緯25°、北緯47°以及東經(jīng)108°55’；坐標(biāo)系原點(diǎn)取中國大地坐標(biāo)原點(diǎn)東經(jīng)108°55’、北緯34°32’。預(yù)報(bào)系統(tǒng)采用36、12、4km三重嵌套網(wǎng)格：36km網(wǎng)格覆蓋中國大部分地區(qū)，為更高分辨率的嵌套網(wǎng)格提供氣象和空氣質(zhì)量的邊界條件；12km網(wǎng)格覆蓋浙江省和上海市全境、江蘇省和安徽省絕大部分地區(qū)以及江西省東北部和福建省北部，實(shí)現(xiàn)浙江省及周邊地區(qū)的空氣質(zhì)量預(yù)報(bào)，能夠以較高的空間分辨率反映浙江省及周邊地區(qū)的污染態(tài)勢，并可用于相互影響分析；4km網(wǎng)格覆蓋浙江省全境，提供浙江省各城市更高分辨率的空氣質(zhì)量預(yù)報(bào)結(jié)果。設(shè)置的網(wǎng)格區(qū)域如圖2所示。垂直模擬網(wǎng)格分層為16層，最底層高度距離地面約40m。

圖2 浙江省多模式預(yù)報(bào)系統(tǒng)嵌套網(wǎng)格設(shè)置

1.3 輸入數(shù)據(jù)

氣象模式最外層的氣象初始和邊界條件采用全球預(yù)報(bào)系統(tǒng)GFS每天的預(yù)報(bào)結(jié)果，系統(tǒng)設(shè)定每天北京時(shí)間12:30開始自動(dòng)下載全球預(yù)報(bào)模式在世界時(shí)當(dāng)天0時(shí)的預(yù)報(bào)結(jié)果。最外層空氣質(zhì)量模式的邊界條件采用全球模式輸出的污染物濃度分布廓線，初始場采用預(yù)報(bào)系統(tǒng)在前一天的預(yù)報(bào)結(jié)果。

對于源排放輸入數(shù)據(jù)，浙江省范圍內(nèi)采用本項(xiàng)目研制的高分辨率污染源排放清單的處理結(jié)果，浙江省以外的地區(qū)采用包括中國在內(nèi)的東亞地區(qū)近年的排放清單資料作為輸入。

1.4 運(yùn)行環(huán)境

預(yù)報(bào)系統(tǒng)的運(yùn)行環(huán)境的操作系統(tǒng)為RedHatEnterpriseLinux5，安裝有GNUC/C++編譯器、GNUFortran77編譯器、IntelC/C++編譯器、IntelFortran77/90編譯器以及MPICH2并行計(jì)算環(huán)境。

2 預(yù)報(bào)性能評估

2.1 評估方法

根據(jù)統(tǒng)計(jì)學(xué)中關(guān)于評估實(shí)測與預(yù)測相關(guān)性原理，本文通過特定的預(yù)報(bào)評估模塊對預(yù)報(bào)系統(tǒng)預(yù)報(bào)結(jié)果的準(zhǔn)確性進(jìn)行評價(jià)，評估模塊采用標(biāo)準(zhǔn)化平均偏差(NMB)、標(biāo)準(zhǔn)化平均誤差(MME)、相關(guān)系數(shù)(COR)來評估污染物預(yù)測結(jié)果與實(shí)測值的吻合程度。計(jì)算公式如下：

a)標(biāo)準(zhǔn)化平均偏差：

(1)

b)標(biāo)準(zhǔn)化平均誤差：

(2)

c)相關(guān)系數(shù)：

(3)

其中：N代表評估的總月份數(shù)，Mi代表第i個(gè)月模式的預(yù)測值，Ma代表總月份的模式預(yù)測平均值，Oi為第i個(gè)月的實(shí)際測量值，Oa代表總月份的實(shí)際測量值的平均值。標(biāo)準(zhǔn)化平均偏差NMB可以反應(yīng)模式預(yù)測的正負(fù)相關(guān)性大小，正數(shù)代表總體預(yù)測較高于實(shí)際測量值，負(fù)數(shù)則相反，標(biāo)準(zhǔn)化平均偏差NME可以反應(yīng)模式預(yù)測的偏差大小，而相關(guān)系數(shù)COR可以直接反應(yīng)模式預(yù)測的準(zhǔn)確度，相關(guān)系數(shù)大，則表示模式預(yù)測值的吻合程度比較高。

2.2 污染物濃度的預(yù)報(bào)評估

收集了浙江省空氣質(zhì)量監(jiān)測網(wǎng)中47個(gè)國家站點(diǎn)2014年1月1日至8月31日期間每天各類污染物監(jiān)測數(shù)據(jù)及各個(gè)預(yù)測模型的預(yù)測值，通過式(1)-(3)計(jì)算各預(yù)測模型對各類污染物的預(yù)報(bào)效果統(tǒng)計(jì)評估參數(shù)，如表1所示，其中表1中“平均集成”模型的各類污染物數(shù)值是由4種模型相對應(yīng)的污染物數(shù)值求平均值后，通過式(1)-式(3)計(jì)算平均集成模型的評估參數(shù)。

表1 各個(gè)模型對不同污染物的預(yù)測評估參數(shù)指標(biāo)

由表1可知，從單個(gè)模式預(yù)測各個(gè)污染物角度出發(fā)，參考各污染物的預(yù)報(bào)評估參數(shù)可以看出，不同模式系統(tǒng)在預(yù)報(bào)過程中對不同污染物的預(yù)報(bào)效果存在一定差異。對SO2的預(yù)測表現(xiàn)最佳的模型系統(tǒng)為MM5-CMAQ和WRF-CMAQ，MM5-CMAQ模型系統(tǒng)對SO2預(yù)測的偏差及誤差為最低，雖然WRF-CMAQ預(yù)測的偏差及誤差大于MM5-CMAQ，但是一致性系數(shù)優(yōu)于MM5-CMAQ；WFR-CMAQ預(yù)測NO2的偏差及誤差均為最低，同時(shí)一致性也最優(yōu)，因此其對NO2的預(yù)測效果最佳，MM5-CMAQ與MM5-CAMx對NO2的預(yù)測偏差及誤差略高于WRF-CMAQ，對NO2預(yù)測的效果較優(yōu)；綜合各項(xiàng)評估參數(shù)，對O3預(yù)測最佳的模型系統(tǒng)為MM5-CMAQ與WRF-CMAQ，這兩套模型對O3的預(yù)測偏差和誤差均較低并且一致性系數(shù)較高，其中MM5-CMAQ的預(yù)測結(jié)果偏低而WRF-Chem的預(yù)測結(jié)果偏高；對CO的預(yù)測表現(xiàn)較好的模型系統(tǒng)為WRF-Chem和WRF-CMAQ，MM5-CMAQ與MM5-CAMx對CO的預(yù)測偏低顯著，可能由于排放清單的低估導(dǎo)致。各套模型系統(tǒng)對顆粒物的預(yù)測均存在偏低現(xiàn)象，其中WRF-Chem的偏低水平最不顯著，但是其預(yù)測誤差高于WRF-CMAQ，綜合來看WRF-CMAQ對顆粒物的預(yù)測效果最佳。

對表1中通過平均集成后獲得的評估參數(shù)與其他單個(gè)模式進(jìn)行比較，發(fā)現(xiàn)就浙江全省的平均預(yù)測效果而言，“平均集成”模型的預(yù)測誤差優(yōu)于任何一個(gè)單模式的預(yù)測結(jié)果，對SO2、NO2和O3的預(yù)測效果提升較為明顯，優(yōu)于或至少相當(dāng)于任一單一模式的預(yù)測結(jié)果。

2.3 空氣質(zhì)量指數(shù)(AQI)與空氣質(zhì)量等級的預(yù)測評估

空氣質(zhì)量指數(shù)(air quality index，AQI)是定量描述空氣質(zhì)量狀況的無量綱指數(shù)。AQI值可以直觀反映出空氣質(zhì)量的等級狀況，所以我們運(yùn)用特定的評估參數(shù)對多模式集成預(yù)報(bào)系統(tǒng)的AQI預(yù)報(bào)進(jìn)行評估。

表2給出了各套模型系統(tǒng)及平均集成模型的AQI預(yù)報(bào)效果的統(tǒng)計(jì)評估結(jié)果。觀察評估參數(shù)后發(fā)現(xiàn)各“平均集成”模型的結(jié)果相對單個(gè)模式預(yù)報(bào)誤差有所減少，與觀測值的相關(guān)性有了明顯提高。

表2 日均AQI預(yù)測效果評估結(jié)果

對每日預(yù)報(bào)的污染等級與實(shí)際觀測等級進(jìn)行比較和逐月統(tǒng)計(jì)，表3是2014年1月至8月各個(gè)模式與“平均集成”模型預(yù)報(bào)的污染等級預(yù)報(bào)準(zhǔn)確率的逐月統(tǒng)計(jì)，比較后發(fā)現(xiàn)就浙江省整體預(yù)報(bào)的情況而言，“平均集成”模型的污染等級預(yù)報(bào)結(jié)果準(zhǔn)確率方面有了顯著提高。

表3 2014年1—8月污染等級預(yù)測準(zhǔn)確率逐月統(tǒng)計(jì) %

3 多模式動(dòng)態(tài)權(quán)重分配模型

近幾十年來，空氣質(zhì)量模型有了快速發(fā)展，雖然在模式結(jié)構(gòu)、數(shù)值算法和物理化學(xué)過程等方面都有很大改進(jìn)，大氣污染模式對主要的污染過程處理還存在很多不確定性。這主要在于不同模式因開發(fā)時(shí)所關(guān)注的主要科學(xué)問題不同，重點(diǎn)考慮的大氣物理、化學(xué)過程有所側(cè)重，而簡化其他一些過程，這使得即使是評估中最優(yōu)的模式也存在模擬、預(yù)報(bào)上的某些缺陷，單一模式并不能完全模擬所有大氣物理、化學(xué)過程，且針對不同地域及各天氣過程各模式預(yù)報(bào)效果差異顯著。這種不確定性的存在，使得目前空氣質(zhì)量模式在空氣質(zhì)量預(yù)報(bào)和空氣質(zhì)量管理政策制定上存在一定的風(fēng)險(xiǎn)。而多模式集合預(yù)報(bào)正成為改進(jìn)預(yù)報(bào)和模擬的一種重要手段，通過對不同模式預(yù)報(bào)結(jié)果，采用合理的集成方法能有效改進(jìn)模式的整體模擬效果。

基于以上評估分析，發(fā)現(xiàn)多模式空氣質(zhì)量預(yù)報(bào)取算術(shù)平均值的引入，相對單個(gè)模式而言提高了浙江省空氣質(zhì)量預(yù)測的準(zhǔn)確度，但是數(shù)值預(yù)報(bào)模式有不同的輸入條件與計(jì)算規(guī)則，模式在預(yù)報(bào)過程中的不同會得到不同的預(yù)報(bào)結(jié)果，通常情況下單個(gè)模式不可能完全準(zhǔn)確地預(yù)報(bào)出空氣質(zhì)量狀況，往往會有偏差，但是單個(gè)模式在不同條件下會呈現(xiàn)出不同的預(yù)測準(zhǔn)確度，如上文分析得到WRF-CMAQ在浙江省空氣質(zhì)量預(yù)測中具有較高的準(zhǔn)確度，只是簡單取各個(gè)模式的算術(shù)平均值，往往不能發(fā)揮單個(gè)模式在特定條件下具有較高準(zhǔn)確度的優(yōu)勢，于是本文提出并且設(shè)計(jì)了一個(gè)動(dòng)態(tài)權(quán)重分配模型，動(dòng)態(tài)權(quán)重分配的優(yōu)勢在于可以通過分析各個(gè)模式在一段時(shí)間與特定條件下的預(yù)測準(zhǔn)確度，來分配集成預(yù)報(bào)的權(quán)重因子。基于以上思路，系統(tǒng)通過分析集成預(yù)報(bào)與單個(gè)模式的準(zhǔn)確度來動(dòng)態(tài)調(diào)整權(quán)重因子。

3.1 模型構(gòu)建

動(dòng)態(tài)權(quán)重分配模型的數(shù)學(xué)模型如下：

定義單個(gè)模式的權(quán)重Wi，各個(gè)模式的權(quán)重因子相加等于1，那么Wi應(yīng)該滿足：

(4)

其中n表示參與集成預(yù)報(bào)的模式數(shù)量。

偏差率是預(yù)測值與實(shí)際測量值做減法運(yùn)算后求絕對值與實(shí)際測量值的除法運(yùn)算結(jié)果。偏差率可以反應(yīng)出單個(gè)模式預(yù)測的誤差大小。設(shè)第i個(gè)預(yù)報(bào)模式第j天預(yù)測值Mij，實(shí)際觀測值Qj，則該i模式第j天的偏差率Rij可用下式計(jì)算：

(5)

根據(jù)模式每天的預(yù)測偏差率，做算數(shù)平均值求單個(gè)模式一段時(shí)間內(nèi)的預(yù)測偏差率。定義第i個(gè)模式一段時(shí)間內(nèi)，總共m天內(nèi)的偏差率Rim為：

(6)

過渡因子是該模型定義的一個(gè)變量，定義過渡因子與權(quán)重因子呈線性關(guān)系并且與預(yù)測偏差率成反比。設(shè)第i個(gè)模式的權(quán)重過度因子Vi，可用下式計(jì)算：

(7)

因?yàn)檫^渡因子與權(quán)重因子呈線性關(guān)系，并且滿足各個(gè)模式權(quán)重因子和為1。第i個(gè)模式預(yù)報(bào)的權(quán)重因子Wi可用下式計(jì)算：

(8)

定義Mi為第i個(gè)模式的預(yù)測結(jié)果，F(xiàn)in為最后的集成結(jié)果，求解公式為：

(9)

其中n表示參與集成預(yù)報(bào)的模式數(shù)，集成結(jié)果通過每個(gè)模型預(yù)測結(jié)果Mi乘以對應(yīng)的權(quán)重因子Wi之后相加得到。

3.2 模型評估

將動(dòng)態(tài)權(quán)重分配模型運(yùn)用到多模式空氣質(zhì)量集成預(yù)報(bào)系統(tǒng)中，并對浙江省空氣質(zhì)量監(jiān)測網(wǎng)中47個(gè)站點(diǎn)2014年1月1日至8月31日的AQI數(shù)據(jù)進(jìn)行評估分析，用動(dòng)態(tài)權(quán)重分配模型與平均集成計(jì)算出來的AQI預(yù)測結(jié)果的評估參數(shù)指標(biāo)對比結(jié)果列于表4。經(jīng)過動(dòng)態(tài)權(quán)重分配模型處理后的AQI的逐日變化情況與觀測結(jié)果的時(shí)間序列對比結(jié)果見圖3。

表4 動(dòng)態(tài)權(quán)重分配模型處理后的預(yù)測結(jié)果的評估參數(shù)

(a)一月份AQI變化趨勢對比

(b)二月份AQI變化趨勢對比

(c)三月份AQI變化趨勢對比

(d)四月份AQI變化趨勢對比

(e)五月份AQI變化趨勢對比

(f)六月份AQI變化趨勢對比

(g)七月份AQI變化趨勢對比

(h)八月份AQI變化趨勢對比

由表4可知，將動(dòng)態(tài)權(quán)重分配模型與平均集成模型預(yù)報(bào)結(jié)果的評估參數(shù)指標(biāo)進(jìn)行比較，發(fā)現(xiàn)動(dòng)態(tài)權(quán)重分配模型的標(biāo)準(zhǔn)化平均偏差(NMB)變小，雖然在整體預(yù)報(bào)方面還存在偏小的趨勢，但相對算術(shù)平均值有了一定的提高，標(biāo)準(zhǔn)化平均誤差(NME)變小，說明預(yù)測誤差有所減小，而相關(guān)系數(shù)(COR)優(yōu)化后的結(jié)果相比平均集成的結(jié)果大，說明在預(yù)測準(zhǔn)確度上有明顯的提高。

圖3的觀測值變化趨勢與預(yù)測變化趨勢，發(fā)現(xiàn)動(dòng)態(tài)權(quán)重分配模型的預(yù)測值與實(shí)際觀測值具有較好的擬合性與較相同的變化趨勢。綜上可以看出，說明動(dòng)態(tài)權(quán)重分配模型預(yù)測準(zhǔn)確度有所提高，比平均集成模型要好。

4 結(jié) 語

本文根據(jù)空氣質(zhì)量數(shù)值預(yù)報(bào)技術(shù)，建立了多模式集成預(yù)報(bào)模型，以WRF-Chem、WRF-CMAQ、MM5-CMAQ和MM5-CAMx 4種數(shù)值預(yù)報(bào)模式為基礎(chǔ)，提出了一個(gè)動(dòng)態(tài)權(quán)重更新模型，對區(qū)域空氣質(zhì)量狀況及污染物濃度進(jìn)行預(yù)測預(yù)報(bào)。以浙江省47個(gè)國家站點(diǎn)8個(gè)月的監(jiān)測數(shù)據(jù)為測試數(shù)據(jù)，分別用WRF-Chem、WRF-CMAQ、MM5-CMAQ、MM5-CAMx、平均集成模型及動(dòng)態(tài)權(quán)重更新模型對浙江省空氣質(zhì)量狀況及污染物濃度進(jìn)行預(yù)測預(yù)報(bào)，并對它們的預(yù)報(bào)結(jié)果進(jìn)行分析。實(shí)驗(yàn)結(jié)果表明，動(dòng)態(tài)權(quán)重更新模型具有較高的準(zhǔn)確度與較小的誤差，可以準(zhǔn)確預(yù)測空氣質(zhì)量狀況與污染物濃度分布情況，為浙江省的區(qū)域空氣質(zhì)量預(yù)測預(yù)報(bào)業(yè)務(wù)提供了有力的技術(shù)支持與保障。

[1] 譚成好,陳昕,趙天良,等.空氣質(zhì)量數(shù)值模型的構(gòu)建及應(yīng)用研究進(jìn)展[J].環(huán)境監(jiān)控與預(yù)警,2014(6):1-7.

[2] 陳彬彬,林長城.基于CMAQ模式產(chǎn)品的福州市空氣質(zhì)量預(yù)報(bào)系統(tǒng)[J].中國環(huán)境科學(xué),2012,32(10):1744-1752.

[3] 程興宏,徐祥德,丁國安.CMAQ模式衛(wèi)星產(chǎn)品源同化模型及其在空氣質(zhì)量預(yù)報(bào)中的應(yīng)用研究[J].中國科學(xué):地球科學(xué),2010,40(4):511-522.

[4] 陳煥盛,王自發(fā),吳其重,等.亞運(yùn)時(shí)段廣州大氣污染物來源數(shù)值模擬研究[J].環(huán)境科學(xué)學(xué)報(bào)，2010,30(11):2145-2153.

[5] 黃蕊珠,陳煥盛,葛寶珠,等.京津冀重霾期間PM2.5來源數(shù)值模擬研究[J].環(huán)境科學(xué)學(xué)報(bào),2015,35(9):2670-2680.

[6] CHENG N，HONGXIA L I，F(xiàn)AN M，et al. Overview of urban PM 2.5 numerical forecast models in China[J]. Asian Agricultural Research,2015,7(10):47-53

[7] 王自發(fā),李杰,王哲,等.2013年1月我國中東部強(qiáng)霾污染的數(shù)值模擬和防控對策[J].中國科學(xué):地球科學(xué),2014,44(11):3-14.

[8] 王自發(fā),吳其重,晏平仲,等.北京空氣質(zhì)量多模式集成預(yù)報(bào)系統(tǒng)的建立及初步應(yīng)用[J].南京信息工程大學(xué)學(xué)報(bào)(自然科學(xué)版),2009(3):19-26.

[9] 王占山,李曉倩,王宗爽,等.空氣質(zhì)量模型CMAQ的國內(nèi)外研究現(xiàn)狀[J].環(huán)境科學(xué)與技術(shù),2013,36(s):386-391.

[10] 王茜,伏晴艷,王自發(fā),等.集合數(shù)值預(yù)報(bào)系統(tǒng)在上海市空氣質(zhì)量預(yù)測預(yù)報(bào)中的應(yīng)用研究[J].環(huán)境監(jiān)控與預(yù)警,2010，2(4):1-6.

[11] 謝濤,姚新,孫世友.基于GIS的空氣質(zhì)量多模式預(yù)報(bào)系統(tǒng)[J].地理信息世界，2011，9(1):11-15.

[12] 黃思,唐曉,徐文帥,等.利用多模式集合和多元線性回歸改進(jìn)北京PM10預(yù)報(bào)[J].環(huán)境科學(xué)學(xué)報(bào),2015,35(1):56-64.

[13] SHRESTHA K L，KONDO A，AKIKAZU K，et al. High-resolution modeling and evaluation of ozone air quality of Osaka using MM5-CMAQ system[J]. Journal of Environmental Sciences，2009，21(6):782-789.

(責(zé)任編輯： 陳和榜)

Research on Multi-Model Air Quality Integrated Forecasting Model

YAOWenqiang1，WANGZhaoqing1，TIEZhixin1，DINGChengfu2，TANGMin1

(1.School of Information Science and Technology，Zhejiang Sci-Tech University，Hangzhou 310018，China; 2.Focused Photonics (Hangzhou) Inc.，Hangzhou 310052，China)

Aiming at the problems existing in the field of air quality forecast,such as inaccurate prediction,larger error and missed to report the heavy pollution weather. In this paper，we used the technology named numerical forecast in the field of air quality forecast,and proposed to use multi-model ensemble forecasting method to forecast the regional air quality and pollutant concentration.The experiments used WRF-Chem、WRF-CMAQ、MM5-CMAQ and MM5-CAMx to do integration. And the dynamic weight updating model is proposed，according to the prediction accuracy of these four kinds of numerical model in a period of time and under certain conditions，weights are dynamically assigned to these models’ prediction value in the new model. The predictive value of the new model is obtained by the weighted average of the predictive value of these four existing prediction model. And the model was used to forecast the air quality of Zhejiang province,eight months for monitoring data of 47 national stations which are of air quality monitoring network in this province was collected.And the forecast results obtained in the multi-model forecasting system were evaluated and analyzed by parameters specified in paper. The experiment show that the Dynamic Weight updating Model has higher prediction accuracy and less error.

air quality;multi-model;numerical forecast;dynamic weight updating

10.3969/j.issn.1673-3851.2017.05.022

2016-08-15 網(wǎng)絡(luò)出版日期： 2017-01-19

浙江省公益技術(shù)應(yīng)用研究項(xiàng)目(2014C31G2060072)

姚文強(qiáng)(1991-)，男，山西省運(yùn)城市人，碩士研究生，主要從事環(huán)保軟件研發(fā)方面的研究。

王兆青，E-mail：zilggjb@163.com

TP311.5

A

1673- 3851 (2017) 03- 0444- 07