環(huán)境氣象應(yīng)用數(shù)據(jù)庫(kù)的設(shè)計(jì)及實(shí)現(xiàn)

陳晴 王俊驕 陳曄峰 楊明 孫曉燕 高婷 吳書(shū)成

摘要：目前氣象部門(mén)能夠獲取的環(huán)境監(jiān)測(cè)數(shù)據(jù)來(lái)源不同，各部門(mén)存儲(chǔ)數(shù)據(jù)的規(guī)范也存在差異，各項(xiàng)數(shù)據(jù)缺乏統(tǒng)一性和完整性。針對(duì)上述情況，討論了數(shù)據(jù)處理的方法，梳理了環(huán)境監(jiān)測(cè)數(shù)據(jù)來(lái)源，統(tǒng)一存儲(chǔ)規(guī)范將不同來(lái)源的環(huán)保數(shù)據(jù)分別進(jìn)行處理，統(tǒng)一存儲(chǔ)在環(huán)境氣象應(yīng)用數(shù)據(jù)庫(kù)中。對(duì)環(huán)境監(jiān)測(cè)信息的數(shù)據(jù)進(jìn)行初步質(zhì)量控制，利用克里金算法對(duì)數(shù)據(jù)插值形成網(wǎng)格化產(chǎn)品，為空氣質(zhì)量預(yù)報(bào)、空氣污染氣象條件預(yù)報(bào)業(yè)務(wù)和相關(guān)科研項(xiàng)目提供基礎(chǔ)數(shù)據(jù)支撐。

關(guān)鍵詞：環(huán)境氣象應(yīng)用數(shù)據(jù)庫(kù)；存儲(chǔ)規(guī)范；質(zhì)量控制；克里金算法；空氣質(zhì)量預(yù)報(bào)

中圖分類(lèi)號(hào)：TP393文獻(xiàn)標(biāo)志碼：A文章編號(hào)：1008-1739（2019）14-59-3

0引言

隨著經(jīng)濟(jì)社會(huì)的不斷發(fā)展，人們的環(huán)保意識(shí)在不斷增強(qiáng)，各級(jí)政府部門(mén)對(duì)于環(huán)境保護(hù)工作也更加重視[1]。圍繞氣象防災(zāi)減災(zāi)需求，進(jìn)一步做好環(huán)境氣象數(shù)據(jù)共享、提升氣象業(yè)務(wù)中環(huán)境氣象數(shù)據(jù)的支撐保障能力，深化在大氣環(huán)境質(zhì)量監(jiān)測(cè)預(yù)警、預(yù)報(bào)服務(wù)、應(yīng)急響應(yīng)等方面的合作，與環(huán)保廳合作實(shí)現(xiàn)數(shù)據(jù)共享，建立了環(huán)境氣象應(yīng)用數(shù)據(jù)庫(kù)。該數(shù)據(jù)庫(kù)的建設(shè)可為全省空氣質(zhì)量預(yù)報(bào)、空氣污染氣象條件預(yù)報(bào)業(yè)務(wù)和相關(guān)科研項(xiàng)目提供基礎(chǔ)數(shù)據(jù)支撐，實(shí)現(xiàn)環(huán)保與氣象大氣環(huán)境監(jiān)測(cè)資料實(shí)時(shí)共享，豐富大氣環(huán)境預(yù)報(bào)預(yù)警服務(wù)信息，有助于共同建立極端氣象條件下空氣質(zhì)量和突發(fā)大氣環(huán)境事件監(jiān)測(cè)預(yù)報(bào)、預(yù)警響應(yīng)、應(yīng)急聯(lián)動(dòng)的工作機(jī)制。

1數(shù)據(jù)庫(kù)設(shè)計(jì)及實(shí)現(xiàn)

1.1數(shù)據(jù)庫(kù)設(shè)計(jì)

環(huán)境氣象應(yīng)用數(shù)據(jù)庫(kù)包括環(huán)保站點(diǎn)信息表、環(huán)保數(shù)據(jù)日?qǐng)?bào)表和環(huán)保數(shù)據(jù)實(shí)時(shí)報(bào)表，分別如表1、表2和表3所示。環(huán)保站點(diǎn)信息表記錄環(huán)保站點(diǎn)的詳細(xì)信息，環(huán)保數(shù)據(jù)日?qǐng)?bào)表按天記錄環(huán)保數(shù)據(jù)，一般當(dāng)天的日?qǐng)?bào)于次日錄入，環(huán)保數(shù)據(jù)實(shí)時(shí)報(bào)表按小時(shí)記錄環(huán)保數(shù)據(jù)，當(dāng)前小時(shí)的數(shù)據(jù)于下個(gè)時(shí)次錄入。

1.2數(shù)據(jù)庫(kù)的實(shí)現(xiàn)

氣象部門(mén)能獲取的環(huán)境數(shù)據(jù)分別來(lái)源于4種不同途徑：①中國(guó)氣象信息中心內(nèi)網(wǎng)下載的全國(guó)環(huán)保站的數(shù)據(jù)；②省環(huán)保監(jiān)測(cè)中心獲取的全省環(huán)保站的數(shù)據(jù)；③省內(nèi)自建的大氣成分站上傳錄入數(shù)據(jù)庫(kù)的本省自建氣象站數(shù)據(jù)；④CMACAST廣播通信系統(tǒng)下發(fā)的報(bào)文中獲取錄入數(shù)據(jù)庫(kù)的周邊省份氣象站數(shù)據(jù)。

通過(guò)分析不同來(lái)源的數(shù)據(jù)屬性，將數(shù)據(jù)按不同方式、相同規(guī)則統(tǒng)一錄入到戶數(shù)據(jù)庫(kù)，方便環(huán)保與氣象大氣環(huán)境監(jiān)測(cè)資料實(shí)時(shí)共享。具體數(shù)據(jù)流程如圖1所示。

（1）全國(guó)環(huán)保站數(shù)據(jù)

通過(guò)網(wǎng)絡(luò)抓取技術(shù)獲取從中國(guó)氣象局內(nèi)網(wǎng)上的全國(guó)環(huán)保站數(shù)據(jù)。網(wǎng)頁(yè)抓取技術(shù)是一個(gè)自動(dòng)抓取網(wǎng)頁(yè)信息的程序，它是搜索引擎的重要組成部分，為搜索引擎從Web上采集信息資源[2]。傳統(tǒng)網(wǎng)頁(yè)抓取技術(shù)從包含一個(gè)或若干個(gè)原始的種子Web頁(yè)面的URL開(kāi)始，最先抓取種子Web頁(yè)面并提取該頁(yè)面的子鏈接，然后將提取的鏈接加入到待抓取的URL隊(duì)列中，在Web頁(yè)面的抓取過(guò)程中，網(wǎng)頁(yè)抓取程序不斷從當(dāng)前抓取的Web頁(yè)面上提取新的子鏈接URL加入到待抓取的URL隊(duì)列中去[2]。重復(fù)上述過(guò)程，直到網(wǎng)頁(yè)抓取系統(tǒng)的停止條件得到滿足時(shí)才停止抓取網(wǎng)頁(yè)。分析網(wǎng)頁(yè)屬性將有用的數(shù)據(jù)信息錄入數(shù)據(jù)庫(kù)。

（2）全省環(huán)保站數(shù)據(jù)

通過(guò)WebService接口調(diào)用獲取數(shù)據(jù)，采用SOAP協(xié)議登錄服務(wù)平臺(tái)的Web Service接口后，實(shí)時(shí)向服務(wù)平臺(tái)發(fā)出數(shù)據(jù)調(diào)用請(qǐng)求，經(jīng)驗(yàn)證通過(guò)后返回所需要的數(shù)據(jù)內(nèi)容。上述基于Web Service方法的數(shù)據(jù)交換均按照XML規(guī)范進(jìn)行。

（3）本省自建氣象站數(shù)據(jù)

從省內(nèi)自建的大氣成分站上傳入庫(kù)至數(shù)據(jù)庫(kù)，編寫(xiě)程序?qū)?shù)據(jù)實(shí)時(shí)同步到環(huán)境氣象應(yīng)用數(shù)據(jù)庫(kù)中。

（4）周邊省份氣象站數(shù)據(jù)

從國(guó)家局的CMACAST廣播系統(tǒng)下發(fā)的報(bào)文中獲取報(bào)文，編寫(xiě)程序?qū)?shí)時(shí)數(shù)據(jù)錄入環(huán)境氣象應(yīng)用數(shù)據(jù)庫(kù)中。

2數(shù)據(jù)質(zhì)量控制

為提高資料的準(zhǔn)確率，對(duì)環(huán)境監(jiān)測(cè)信息的數(shù)據(jù)進(jìn)行初步質(zhì)量控制。大氣中的觀測(cè)數(shù)據(jù)與時(shí)間顯著相關(guān)，具有良好的時(shí)間一致性[3]。根據(jù)這一特性，將數(shù)據(jù)與其時(shí)間上前、后的測(cè)值相比較，來(lái)判斷其數(shù)據(jù)是否發(fā)生異常[4]。時(shí)變檢查主要是根據(jù)要素在某一時(shí)段內(nèi)可能變化范圍判斷要素值質(zhì)量[5]。

在實(shí)際過(guò)程中，設(shè)某一要素0時(shí)次的觀測(cè)值為0，0-1時(shí)次的觀測(cè)值為-1，0+1時(shí)次的觀測(cè)值為+1。要檢驗(yàn)0的值是否在正常范圍內(nèi)，設(shè)定：avg=（-1++1）/2，0>=2avg，初步判定0為異常；0<2avg，初步判定0為正常。對(duì)于初步判定是異常的數(shù)據(jù)，將該數(shù)據(jù)設(shè)為缺測(cè)。

3環(huán)境氣象數(shù)據(jù)庫(kù)的應(yīng)用

3.1基于克里金算法的網(wǎng)格化產(chǎn)品生成

由于站點(diǎn)的數(shù)據(jù)量有限，空間分布是離散的，利用插值通過(guò)已知樣點(diǎn)的數(shù)據(jù)估算未知點(diǎn)的數(shù)據(jù)，通過(guò)一次或多次插值將離散點(diǎn)數(shù)據(jù)插值成連續(xù)的表面數(shù)據(jù)，進(jìn)一步形成等值線或者色斑圖后顯示[6]。對(duì)已有的環(huán)保監(jiān)測(cè)資料進(jìn)行克里金算法的空間插值，將離散的環(huán)保站點(diǎn)數(shù)據(jù)合成連續(xù)的表面數(shù)據(jù)，進(jìn)一步形成色斑圖后顯示。目前產(chǎn)品已在省環(huán)境氣象預(yù)報(bào)業(yè)務(wù)平臺(tái)等業(yè)務(wù)系統(tǒng)中得到應(yīng)用，豐富大氣環(huán)境預(yù)報(bào)預(yù)警服務(wù)信息，為政府決策和社會(huì)公眾生活提供環(huán)境預(yù)報(bào)預(yù)警服務(wù)，為生態(tài)省建設(shè)工作提供數(shù)據(jù)保障。

3.2克里金算法的插值過(guò)程

克里金插值的前提是根據(jù)空間場(chǎng)的結(jié)構(gòu)，選擇適當(dāng)?shù)淖儺惡瘮?shù)模型，并求出變異函數(shù)。其計(jì)算步驟如下：

①網(wǎng)格化，選擇區(qū)域的范圍及網(wǎng)格的大小。

②計(jì)算被估點(diǎn)（網(wǎng)格節(jié)點(diǎn)的）的坐標(biāo)。

③根據(jù)搜索策略（近點(diǎn)距離搜索和方位搜索）選擇合適的參估點(diǎn)。

④根據(jù)已經(jīng)求出的變異函數(shù)，求出方程組的系數(shù)。

⑦重復(fù)步驟②～⑥，直到網(wǎng)格節(jié)點(diǎn)的值全部求出。

4結(jié)束語(yǔ)

建立環(huán)境氣象數(shù)據(jù)的入庫(kù)機(jī)制，梳理環(huán)境監(jiān)測(cè)數(shù)據(jù)來(lái)源，統(tǒng)一存儲(chǔ)規(guī)范，將不同來(lái)源的環(huán)保數(shù)據(jù)按照統(tǒng)一的存儲(chǔ)規(guī)范錄入到環(huán)境氣象應(yīng)用數(shù)據(jù)庫(kù)中，實(shí)現(xiàn)數(shù)據(jù)的集中管理。利用時(shí)變檢查環(huán)境監(jiān)測(cè)數(shù)據(jù)進(jìn)行初步質(zhì)量控制，生成環(huán)境監(jiān)測(cè)要素的網(wǎng)格化產(chǎn)品，為全省空氣質(zhì)量預(yù)報(bào)提供技術(shù)支撐，豐富大氣環(huán)境預(yù)報(bào)預(yù)警服務(wù)信息，為政府決策和社會(huì)公眾生活提供環(huán)境預(yù)報(bào)預(yù)警服務(wù)。

參考文獻(xiàn)

[1]張霞.基層環(huán)境監(jiān)測(cè)站存在的問(wèn)題與改革措施[C]//2008中國(guó)環(huán)境科學(xué)學(xué)會(huì)學(xué)術(shù)年會(huì)優(yōu)秀論文集（下卷），2008： 319 -320.

[2]謝治軍.垂直搜索引擎的主題網(wǎng)頁(yè)抓取策略研究[D].重慶：重慶理工大學(xué)，2012：8-31.

[3]范文波.地面氣象觀測(cè)數(shù)據(jù)綜合質(zhì)量控制方法研究與實(shí)現(xiàn)[D].江蘇：南京信息工程大學(xué)，2016：88-93.

[4]周青，張樂(lè)堅(jiān)，李峰，等.自動(dòng)站實(shí)時(shí)數(shù)據(jù)質(zhì)量分析及質(zhì)控算法改進(jìn)[J].氣象科技，2015，43（5）：814-822.

[5]王海軍，楊志彪，楊代才，等.自動(dòng)氣象站實(shí)時(shí)資料自動(dòng)質(zhì)量控制方法及其應(yīng)用[J].氣象，2007，33（10）：102-109.

[6]吳煥萍，羅兵，王維國(guó)，等.GIS技術(shù)在決策氣象服務(wù)系統(tǒng)建設(shè)中的應(yīng)用[J].應(yīng)用氣象學(xué)報(bào)，2008，19（3）：380-384.