民航黑龍江空中交通管理分局 黃紅兵

?

哈爾濱機(jī)場(chǎng)相似預(yù)報(bào)系統(tǒng)的設(shè)計(jì)與研制

民航黑龍江空中交通管理分局 黃紅兵

【摘要】采用NECP資料建立歷史樣本數(shù)據(jù)庫(kù)，采用T639數(shù)值預(yù)報(bào)產(chǎn)品為當(dāng)前數(shù)據(jù)，對(duì)各數(shù)據(jù)場(chǎng)都進(jìn)行了標(biāo)準(zhǔn)化處理；采用相似系數(shù)、歐氏距離、相似離度等3個(gè)相似判據(jù)，經(jīng)過值域轉(zhuǎn)換后定義了3個(gè)相似指數(shù)；采用7個(gè)氣象要素、間隔12h的2個(gè)相鄰時(shí)刻的數(shù)據(jù)場(chǎng)，設(shè)計(jì)了1個(gè)建立在3個(gè)相似指數(shù)基礎(chǔ)上、用以描述天氣形勢(shì)演變過程的相似程度的綜合相似指數(shù)，并設(shè)計(jì)了哈爾濱機(jī)場(chǎng)24h短期天氣相似預(yù)報(bào)模式；利用哈爾濱機(jī)場(chǎng)現(xiàn)有的氣象網(wǎng)絡(luò)設(shè)備，研制了哈爾濱機(jī)場(chǎng)相似預(yù)報(bào)系統(tǒng)，實(shí)現(xiàn)自動(dòng)讀入數(shù)據(jù)、自動(dòng)分析、自動(dòng)發(fā)布預(yù)報(bào)等功能；本文中演示了該系統(tǒng)在哈爾濱機(jī)場(chǎng)的應(yīng)用效果。

【關(guān)鍵詞】相似分析；多元判據(jù)；短期天氣；預(yù)報(bào)系統(tǒng)

1 引言

數(shù)值預(yù)報(bào)是天氣預(yù)報(bào)的發(fā)展主流方向，是現(xiàn)代天氣預(yù)報(bào)業(yè)務(wù)的重要基礎(chǔ)[1]。雖然數(shù)值預(yù)報(bào)產(chǎn)品具有準(zhǔn)確性好、時(shí)空分辨率高、連續(xù)性優(yōu)良、應(yīng)用廣泛等特點(diǎn)，但是目前仍不能充分滿足日常業(yè)務(wù)需求，因此開發(fā)數(shù)值預(yù)報(bào)產(chǎn)品的釋用技術(shù)成為當(dāng)前的重要研究課題，地方氣象部門已對(duì)此進(jìn)行了多方面研究，其中，相似預(yù)報(bào)技術(shù)[2-8]成為其采用的主要釋用方法之一。由于種種原因，相似預(yù)報(bào)技術(shù)在民航氣象界還未得到有效開發(fā)，為將該技術(shù)應(yīng)用到航空天氣預(yù)報(bào)中，民航黑龍江空管分局組成一個(gè)科研組，于2013年申請(qǐng)獲得了民航局空管局研發(fā)哈爾濱機(jī)場(chǎng)相似預(yù)報(bào)系統(tǒng)的科研項(xiàng)目?？蒲薪M基于前期技術(shù)試驗(yàn)成果[9]，重新設(shè)計(jì)了多元判據(jù)相似預(yù)報(bào)模式，研發(fā)了哈爾濱機(jī)場(chǎng)相似預(yù)報(bào)系統(tǒng)，并于2015年7月在哈爾濱機(jī)場(chǎng)投入24h短期預(yù)報(bào)試用。

2 預(yù)報(bào)模式設(shè)計(jì)

2.1 基本設(shè)計(jì)思路

基于前期試驗(yàn)成果，采用3個(gè)相似判據(jù)，經(jīng)值域轉(zhuǎn)換后定義3個(gè)新相似指數(shù)；采用7個(gè)氣象要素?cái)?shù)據(jù)場(chǎng)，以其3個(gè)相似指數(shù)為基礎(chǔ)，定義1個(gè)綜合相似指數(shù)；計(jì)算代表天氣演變過程的2個(gè)相鄰時(shí)刻的綜合相似指數(shù)，并取其平均值為天氣演變的過程綜合相似指數(shù)；建立以起始時(shí)刻、結(jié)束時(shí)刻和演變過程的3個(gè)綜合相似指數(shù)為依據(jù)的相似樣本確定原則，選取當(dāng)前最為相似的歷史樣本；基于相似樣本，計(jì)算并發(fā)出機(jī)場(chǎng)未來(lái)24h天氣預(yù)報(bào)信息。

2.2 樣本數(shù)據(jù)資料

2.2.1 數(shù)據(jù)資料來(lái)源

科研組采用了美國(guó)國(guó)家環(huán)境預(yù)報(bào)中心（NCEP）的數(shù)據(jù)來(lái)建立歷史樣本數(shù)據(jù)庫(kù)，這主要是因?yàn)镹CEP歷史數(shù)據(jù)資料齊全；同時(shí)，考慮到發(fā)布航空天氣預(yù)報(bào)時(shí)效性的要求，科研組采用了中國(guó)氣象局T639數(shù)值預(yù)報(bào)產(chǎn)品作為當(dāng)前預(yù)報(bào)樣本數(shù)據(jù)，這主要是因?yàn)楫?dāng)前最近的T639數(shù)據(jù)能夠在較短時(shí)間內(nèi)獲取。由于NCEP與T639數(shù)據(jù)格式、數(shù)據(jù)范圍、網(wǎng)格分辨率不同，因此在建立歷史樣本數(shù)據(jù)庫(kù)過程中，科研組參照T639數(shù)值預(yù)報(bào)產(chǎn)品目前的數(shù)據(jù)格式，對(duì)NCEP歷史樣本數(shù)據(jù)進(jìn)行了統(tǒng)一格式的轉(zhuǎn)換。

2.2.2 采用氣象要素

目前采用的氣象要素場(chǎng)有：500hPa、700hPa和850hPa的高度場(chǎng)和溫度場(chǎng)，海平面氣壓場(chǎng)等共7類氣象要素的數(shù)值預(yù)報(bào)場(chǎng)。計(jì)劃在下一階段的升級(jí)改造中嘗試引進(jìn)渦度場(chǎng)、散度場(chǎng)、水汽通量場(chǎng)、垂直速度場(chǎng)、水汽通量散度場(chǎng)等的數(shù)值預(yù)報(bào)產(chǎn)品。

2.2.3 樣本年限時(shí)刻

當(dāng)前預(yù)報(bào)樣本：采用預(yù)報(bào)有效時(shí)段的起始時(shí)刻前12 ~ 24小時(shí)發(fā)布的12 ~ 36h數(shù)值預(yù)報(bào)的格點(diǎn)場(chǎng)數(shù)據(jù)，其中以12小時(shí)前發(fā)布的12、24h預(yù)報(bào)為主用數(shù)據(jù)，以18小時(shí)前發(fā)布的18、30h預(yù)報(bào)為備用數(shù)據(jù)，以24小時(shí)前發(fā)布的24、36h預(yù)報(bào)為應(yīng)急數(shù)據(jù)。每天使用的數(shù)據(jù)時(shí)間包括00、06、12、18UTC共4個(gè)時(shí)刻。

歷史樣本數(shù)據(jù)：歷史樣本數(shù)據(jù)年限越遠(yuǎn)久，數(shù)據(jù)量越豐富，對(duì)找到最為相似的歷史樣本越有利。但鑒于哈爾濱太平國(guó)際機(jī)場(chǎng)24h氣象觀測(cè)始于2000年，之前為13h觀測(cè)、機(jī)場(chǎng)觀測(cè)數(shù)據(jù)缺測(cè)較多，因此只能使用2000年以來(lái)、截止2015年共16a的歷史數(shù)據(jù)來(lái)建立歷史樣本數(shù)據(jù)庫(kù)，并逐年累加。歷史樣本數(shù)據(jù)時(shí)刻與當(dāng)前預(yù)報(bào)樣本數(shù)據(jù)時(shí)刻同步。

2.2.4 數(shù)據(jù)場(chǎng)的范圍

數(shù)據(jù)場(chǎng)數(shù)據(jù)范圍的劃分既不宜太大，又不宜太小，如果劃分太大可能摻雜進(jìn)一些遠(yuǎn)距離無(wú)關(guān)信息，如果劃分太小可能造成一些重要信息丟失，這兩種情況都將影響相似分析效果。

鑒于預(yù)報(bào)對(duì)象為24h短期天氣，科研組經(jīng)過前期技術(shù)試驗(yàn)[9]，將用于大尺度形勢(shì)分析的高空數(shù)據(jù)場(chǎng)經(jīng)緯度范圍劃定為100~147.5°E、32.5~60°N（分辨率2.5°×2.5°），將海平面氣壓場(chǎng)關(guān)鍵區(qū)經(jīng)緯度范圍劃定為110~137.5°E、37.5~55°N（分辨率2.5°×2.5°）。

2.2.5 數(shù)據(jù)場(chǎng)標(biāo)準(zhǔn)化

考慮到不同氣象要素場(chǎng)量級(jí)、值域的差異和各相似判據(jù)的可比性要求，在進(jìn)行多元判別相似分析前，科研組首先對(duì)預(yù)報(bào)樣本和歷史樣本的原始數(shù)據(jù)采取了標(biāo)準(zhǔn)化處理，計(jì)算方法略。經(jīng)過標(biāo)準(zhǔn)化處理后，數(shù)據(jù)場(chǎng)各點(diǎn)數(shù)據(jù)Xk（k = 1, 2, 3, … , K）的值域調(diào)整為[0, 1]，基于標(biāo)準(zhǔn)化處理后的兩個(gè)樣本數(shù)據(jù)場(chǎng)的相似判據(jù)，其值域也都是可確定的，這對(duì)下一步建立多元判據(jù)相似預(yù)報(bào)模式有利。

2.3 相似分析判據(jù)

常用相似預(yù)報(bào)判據(jù)主要有相似系數(shù)、相關(guān)系數(shù)、歐氏距離、海明距離、相似離度、羅陽(yáng)相似量等?？蒲薪M在經(jīng)過大量試驗(yàn)后，確定使用相似系數(shù)、歐氏距離、相似離度等3個(gè)相似判據(jù)。

2.3.1 相似系數(shù)

相似系數(shù)Ri,j的計(jì)算表達(dá)式略。Ri,j主要描述i、j兩個(gè)樣本“形”的相似程度，其值域?yàn)閇-1, 1]，數(shù)值越大表示這兩個(gè)樣本越相似，數(shù)值越小則越不相似。

2.3.2 歐氏距離

歐氏距離Oi,j的計(jì)算表達(dá)式略。Oi,j主要描述i、j兩個(gè)樣本“值”的相似程度，其數(shù)值越小表示兩個(gè)樣本相似程度越高，數(shù)值越大則相似程度越低。如果樣本數(shù)據(jù)沒有經(jīng)過標(biāo)準(zhǔn)化處理，那么Oi,j值域應(yīng)該為[0, d），其中d的數(shù)值不確定。在樣本數(shù)據(jù)經(jīng)過前期標(biāo)準(zhǔn)化處理后其值域調(diào)整為[0, 1]，Oi,j值域隨之調(diào)整為[0, K1/2]，K為數(shù)據(jù)場(chǎng)數(shù)據(jù)點(diǎn)數(shù)。

2.3.3 相似離度

相似離度Ci,j的計(jì)算表達(dá)式略。Ci,j考慮了i、j兩個(gè)樣本“形”的相似程度，又考慮了“值”的相似程度，其數(shù)值越小表示兩個(gè)樣本相似程度越高，數(shù)值越大則相似程度越低。如果樣本數(shù)據(jù)沒有經(jīng)過標(biāo)準(zhǔn)化處理，那么Ci,j的值域應(yīng)該為[0, d），其中d數(shù)值不確定。在樣本數(shù)據(jù)經(jīng)過前期標(biāo)準(zhǔn)化處理后其值域調(diào)整為[0, 1]，Ci,j值域隨之調(diào)整為[0, 1]。

2.4 多元判據(jù)分析

2.4.1 定義相似指數(shù)

在應(yīng)用多元相似判據(jù)過程中，由于各判據(jù)值域、量級(jí)不一致，不具備直接可比性，因此不宜直接接合起來(lái)使用[2]。為此，科研組基于相似系數(shù)Ri,j、歐氏距離Oi,j和相似離度Ci,j，經(jīng)過再處理后設(shè)計(jì)定義了3個(gè)新的相似指數(shù)。

1）基于相似系數(shù)的相似指數(shù)SRi,j，其計(jì)算表達(dá)式為：

式（1）中：Rmax、Rmin分別為相似系數(shù)最大值和最小值，可取理論上的最大值和最小值，分別為Rmax= 1、Rmin= -1。

2）基于歐氏距離的相似指數(shù)SOi,j，其計(jì)算表達(dá)式為：

式（2）中：K為數(shù)據(jù)場(chǎng)數(shù)據(jù)點(diǎn)數(shù)。

3）基于相似離度的相似指數(shù)SCi,j，其計(jì)算表達(dá)式為：

經(jīng)過式（1）~（3）的再處理后得到的相似指數(shù)SRi,j、SOi,j、SCi,j，其值域、量級(jí)是一致的，其中值域都調(diào)整為[0, 1]，而且都轉(zhuǎn)換為數(shù)值越大相似程度越高、數(shù)值越小相似程度越低，這樣有利于下一步建立多元判據(jù)綜合相似預(yù)報(bào)模式。

2.4.2 綜合相似分析

為了綜合評(píng)估各相似判據(jù)及基于各相似判據(jù)的相似指數(shù)的綜合作用，科研組在上述3個(gè)相似判據(jù)和基于這些判據(jù)設(shè)計(jì)的3個(gè)相似指數(shù)基礎(chǔ)上，定義了1個(gè)多元判據(jù)綜合相似指數(shù)Si,j，其計(jì)算表達(dá)式為：

式（4）中：N為氣象要素場(chǎng)數(shù)量，N = 7；n為序號(hào)，n = 1 ~ 7分別代表500hPa高度場(chǎng)、溫度場(chǎng)，700hPa高度場(chǎng)、溫度場(chǎng)，850hPa高度場(chǎng)、溫度場(chǎng)，海平面氣壓場(chǎng)；WRn、WOn、WCn分別為各氣象要素場(chǎng)的相似指數(shù)SRi,j、SOi,j、SCi,j在綜合相似指數(shù)Si,j中所占有的權(quán)重系數(shù)。

借鑒多元決策加權(quán)法，科研組經(jīng)過大量試驗(yàn)，在分別驗(yàn)證各氣象要素場(chǎng)的各相似指數(shù)的相似分析效果后，統(tǒng)計(jì)出各權(quán)重系數(shù)（數(shù)據(jù)表略），其中權(quán)重系數(shù)之和為1。

綜合相似指數(shù)Si,j是在綜合分析了i、j兩個(gè)樣本的各氣象要素場(chǎng)的各相似判據(jù)、相似指數(shù)后，綜合給出的一個(gè)描述整體相似程度的相似指標(biāo)，其值域?yàn)閇0, 1]，且數(shù)值越大表示兩個(gè)樣本的相似程度越高，數(shù)值越小表示相似程度越低。

2.5 相似預(yù)報(bào)方案

2.5.1 相似時(shí)域確定

由于哈爾濱機(jī)場(chǎng)部分重要天氣年際差異較大，比如在2000~2013年記錄中，初雷最早曾發(fā)生于3月底，最晚曾發(fā)生于5月下旬，因此相似時(shí)域不宜太短。經(jīng)科研組反復(fù)試驗(yàn)證明，將相似時(shí)域確定在當(dāng)前日期±15d的范圍內(nèi)較為適宜。

2.5.2 相似樣本確定

首先計(jì)算綜合相似指數(shù)Sn,m（n = 0, 1, 2；m = 1, 2, 3, … M，M為歷史樣本數(shù)），包括起始時(shí)刻的綜合相似指數(shù)S0,m、結(jié)束時(shí)刻的綜合相似指數(shù)S1,m、從起始到結(jié)束時(shí)刻的天氣演變過程的過程綜合相似指數(shù)S2,m（暫定S2,m= (S0,m+ S1,m)/2）。然后將Sn,m按照從大到小重新排序，同時(shí)調(diào)整Sn,m對(duì)應(yīng)的日期序列Dn,m(n、m與前述相同)。最后依據(jù)天氣形勢(shì)演變過程的最大指數(shù)相似原則、近似最大相似原則、高值指數(shù)相似原則、樣本數(shù)量限制原則，依次篩選、確定當(dāng)前的最佳相似歷史樣本（詳細(xì)過程略）。

2.5.3 相似預(yù)報(bào)方法

（1）天氣形勢(shì)分析：直接把第1個(gè)歷史樣本、也就是過程綜合相似指數(shù)最大的歷史樣本的天氣形勢(shì)前期分析結(jié)果，作為當(dāng)前最為相似的天氣形勢(shì)分析結(jié)果予以輸出。

（2）重要天氣現(xiàn)象：考慮到民航氣象服務(wù)工作的特點(diǎn)，對(duì)重要天氣現(xiàn)象的預(yù)報(bào)，科研組采用計(jì)算并集的方法來(lái)確定，即對(duì)于前m個(gè)最為相似的歷史樣本，它們?cè)谥付ǖ哪硞€(gè)時(shí)段內(nèi)出現(xiàn)的各種重要天氣現(xiàn)象，是當(dāng)前氣象條件下同時(shí)段內(nèi)最可能出現(xiàn)的重要天氣現(xiàn)象。

（3）氣象要素預(yù)報(bào)：

1）對(duì)日降水量、溫度、風(fēng)速、相對(duì)濕度等氣象要素的預(yù)報(bào)，科研組用設(shè)計(jì)的以下預(yù)報(bào)公式來(lái)計(jì)算：

2）對(duì)低能見度、低云高度等氣象要素的預(yù)報(bào)，可以采用公式（5）來(lái)計(jì)算?？紤]到民航氣象服務(wù)工作的特點(diǎn)，科研組采用前m個(gè)最為相似的歷史樣本的最小值或最低值，作為當(dāng)前氣象條件下可能出現(xiàn)的最小或最低預(yù)報(bào)值。

2.5.4 相似預(yù)報(bào)信度

借鑒鐘元的方法[2]，根據(jù)天氣演變過程第1綜合相似指數(shù)，將相似預(yù)報(bào)信度劃為3類：如果S2,1≥0.90，則相似預(yù)報(bào)信度較高；如果0.80 ≤S2,1< 0.90，則相似預(yù)報(bào)信度一般；如果S2,1< 0.80，則相似預(yù)報(bào)信度較差。

3 預(yù)報(bào)系統(tǒng)研制

3.1 硬件結(jié)構(gòu)布局

科研組利用哈爾濱機(jī)場(chǎng)現(xiàn)有的氣象設(shè)備、設(shè)施，將預(yù)報(bào)系統(tǒng)的相關(guān)硬件作如圖1所示的結(jié)構(gòu)布局，其中：（1）當(dāng)前預(yù)報(bào)樣本數(shù)據(jù)，由預(yù)報(bào)系統(tǒng)訪問氣象信息網(wǎng)的地方信息終端，讀取省氣象局轉(zhuǎn)發(fā)的T639數(shù)值預(yù)報(bào)產(chǎn)品格點(diǎn)場(chǎng)數(shù)據(jù)；（2）歷史樣本數(shù)據(jù)，已事先保存在當(dāng)前電腦指定的數(shù)據(jù)庫(kù)位置，并逐年累積；（3）機(jī)場(chǎng)觀測(cè)數(shù)據(jù)，由預(yù)報(bào)系統(tǒng)訪問氣象信息網(wǎng)的觀測(cè)終端，讀取機(jī)場(chǎng)觀測(cè)數(shù)據(jù)；（4）配置的音響和打印機(jī)，以備告警、提醒和預(yù)報(bào)產(chǎn)品輸出之用。

圖1 預(yù)報(bào)系統(tǒng)硬件結(jié)構(gòu)示意圖

3.2 軟件研發(fā)情況

科研組在研發(fā)相似預(yù)報(bào)系統(tǒng)軟件過程中先后推出了HAFS1.01 ~ 1.06共6個(gè)版本，其中最初的5個(gè)版本均屬于測(cè)試版本，HAFS1.06為當(dāng)前最新的應(yīng)用版本。預(yù)報(bào)系統(tǒng)軟件始研于2013年，初步完成于2015 年7月并聯(lián)網(wǎng)試用，經(jīng)2015年7 ~ 12月試用、再測(cè)試和不斷完善，各項(xiàng)工作最后完成于2015年12月。

預(yù)報(bào)系統(tǒng)軟件設(shè)計(jì)在Windows XP及以上環(huán)境下運(yùn)行，開發(fā)語(yǔ)言使用VB6.0，相關(guān)支持平臺(tái)主要有Oracle、Office Excel，其主窗體界面見圖2。該系統(tǒng)軟件基本結(jié)構(gòu)按其功能和研發(fā)的技術(shù)需要，設(shè)計(jì)劃分為網(wǎng)絡(luò)和數(shù)據(jù)監(jiān)控單元、數(shù)據(jù)采集和分析單元、相似計(jì)算和預(yù)報(bào)單元、查詢打印設(shè)置等綜合單元。限于論文篇幅，預(yù)報(bào)系統(tǒng)軟件及其各單元的研發(fā)技術(shù)原理、技術(shù)措施和操作方法等不再詳細(xì)介紹。

3.3 主要技術(shù)指標(biāo)

1）具有標(biāo)準(zhǔn)、美觀、友好的Windows人機(jī)界面：相似預(yù)報(bào)系統(tǒng)實(shí)現(xiàn)簡(jiǎn)體中文顯示，具有標(biāo)準(zhǔn)、美觀、友好和便捷的人機(jī)互動(dòng)功能；操作方法簡(jiǎn)便易懂，只需要短時(shí)間培訓(xùn)，就能夠很快掌握操作和維護(hù)方法；

2）具有完備的網(wǎng)絡(luò)和數(shù)據(jù)監(jiān)控、告警功能：監(jiān)控單元連續(xù)監(jiān)控需要訪問的氣象信息網(wǎng)和有關(guān)終端，定時(shí)監(jiān)控T639數(shù)值預(yù)報(bào)和機(jī)場(chǎng)觀測(cè)的數(shù)據(jù)，如果網(wǎng)絡(luò)正常、數(shù)據(jù)齊全、數(shù)據(jù)格式合格則發(fā)出可以讀取數(shù)據(jù)的信息，否則就以屏顯和音頻告警方式發(fā)出告警信息；

3）具有自動(dòng)讀取、分析綜合數(shù)據(jù)的功能：當(dāng)監(jiān)控單元發(fā)出網(wǎng)絡(luò)和數(shù)據(jù)正常、可以讀取數(shù)據(jù)的信息時(shí)，數(shù)據(jù)單元將按照預(yù)置時(shí)間段，依據(jù)“一主、二備、三應(yīng)急”的方式，讀取、分析所需各類數(shù)據(jù)，并發(fā)出可以制作預(yù)報(bào)的信息；

4）具有自動(dòng)制作、發(fā)布并輸出全方位相似預(yù)報(bào)產(chǎn)品的功能：當(dāng)數(shù)據(jù)單元發(fā)可以制作預(yù)報(bào)的信息時(shí)，預(yù)報(bào)單元隨即啟動(dòng)自動(dòng)發(fā)布預(yù)報(bào)的程序，完成之后自動(dòng)音頻提醒、打印輸出；

5）具有的其它綜合功能：綜合單元可隨時(shí)進(jìn)行參數(shù)設(shè)置、打印輸出、查詢相似分析數(shù)據(jù)和預(yù)報(bào)產(chǎn)品、查看對(duì)比相似圖形、提供幫助信息等。

圖2 預(yù)報(bào)系統(tǒng)主窗體示意圖

4 試用預(yù)報(bào)效果

2015年7~12月，相似預(yù)報(bào)系統(tǒng)在哈爾濱太平國(guó)際機(jī)場(chǎng)投入試運(yùn)行，試運(yùn)行期間，根據(jù)民航局空管局關(guān)于發(fā)布航空天氣預(yù)報(bào)時(shí)間和內(nèi)容的相關(guān)要求，預(yù)報(bào)系統(tǒng)每天于04:00 ~ 06:00、10:00 ~ 12:00、16:00 ~ 18:00、20:00 ~ 24:00時(shí)段（北京時(shí)，下同），分別發(fā)布機(jī)場(chǎng)08:00 ~ 08:00、14:00 ~ 14:00、20:00 ~ 20:00、02:00 ~ 02:00航空天氣預(yù)報(bào)，預(yù)報(bào)項(xiàng)目包括：天氣形勢(shì)分析、重要天氣預(yù)報(bào)、重要?dú)庀笠仡A(yù)報(bào)、主要要素逐時(shí)預(yù)報(bào)，取得了預(yù)期的效果。下面以在20:00 ~ 24:00時(shí)段內(nèi)利用12小時(shí)前發(fā)布的12、24h T639數(shù)值預(yù)報(bào)產(chǎn)品，來(lái)發(fā)布哈爾濱機(jī)場(chǎng)次日02:00 ~ 02:00預(yù)報(bào)為例，相似預(yù)報(bào)系統(tǒng)對(duì)出現(xiàn)的重要天氣現(xiàn)象的預(yù)報(bào)效果見表1，對(duì)達(dá)到中雨量級(jí)及以上的日降水預(yù)報(bào)效果見表2。

表1 2015年6~12月重要天氣預(yù)報(bào)統(tǒng)計(jì)表

表2 2015年6~12月≥10mm日降水量預(yù)報(bào)統(tǒng)計(jì)表

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

1）相似預(yù)報(bào)技術(shù)至今已發(fā)展成為一種成熟的預(yù)報(bào)技術(shù)?？蒲薪M基于相似系數(shù)、歐氏距離和相似離度3個(gè)相似判據(jù)，定義了3個(gè)相似指數(shù)、1個(gè)綜合相似指數(shù)和1個(gè)過程綜合指數(shù)，建立了綜合相似分析預(yù)報(bào)模式，研發(fā)了相似預(yù)報(bào)系統(tǒng)，經(jīng)試用證明，對(duì)短期重要天氣和主要?dú)庀笠鼐哂幸欢ǖ念A(yù)報(bào)效果，能夠?yàn)轭A(yù)報(bào)員提供一個(gè)可供參考的預(yù)報(bào)結(jié)果。

2）目前雖取得一些效果，但也發(fā)現(xiàn)一些問題：一是歷史樣本數(shù)據(jù)年限、數(shù)據(jù)量偏少；二是利用500hPa、700hPa要素場(chǎng)的相似分析效果比850hPa和地面要素場(chǎng)的效果差，這可能與數(shù)據(jù)場(chǎng)范圍大小有關(guān)；三是部分重要天氣例如雷暴大風(fēng)的預(yù)報(bào)效果較差，這是因?yàn)槟壳斑€沒設(shè)計(jì)專門針對(duì)小概率事件的預(yù)報(bào)方案。歷史樣本可以逐年積累，其它問題有待試驗(yàn)、改進(jìn)。

3）科研組計(jì)劃在運(yùn)行滿一年后，擇機(jī)升級(jí)相似預(yù)報(bào)系統(tǒng)，屆時(shí)除了繼續(xù)研究、完善前期科研成果外，還將嘗試引進(jìn)新的氣象要素場(chǎng)，嘗試采用多級(jí)、逐步相似過濾法等來(lái)改進(jìn)預(yù)報(bào)模式，使我們的相似預(yù)報(bào)系統(tǒng)逐步趨于完善、準(zhǔn)確率不斷提高。

參考文獻(xiàn)

[1]陳德輝,薛繼善.數(shù)值天氣預(yù)報(bào)業(yè)務(wù)模式現(xiàn)狀與展望[J].氣象學(xué)報(bào),2004,62(5):623-633.

[2]鐘元.天氣與短期氣候的多元判據(jù)綜合相似預(yù)報(bào)方法[M].北京:氣象出版社,2003.

[3]閻惠芳,李社宗,黃躍青等.常用相似性判據(jù)的檢驗(yàn)和綜合相似系數(shù)的使用[J].氣象科技,2003,31(4): 211-215.

[4]萬(wàn)日金,何溪澄,林鋼.用動(dòng)力相似方法預(yù)報(bào)廣東省區(qū)域暴雨預(yù)報(bào)試驗(yàn)[J].熱帶氣象學(xué)報(bào),2006,22(2): 198-202.

[5]李博,趙思雅,陸漢城等.綜合多級(jí)相似預(yù)報(bào)技術(shù)在暴雨短期預(yù)報(bào)中的檢驗(yàn)[J].應(yīng)用氣象學(xué)報(bào), 2008,19(3):307-314.

[6]陳錄元,周海,王式功等.逐步相似過濾法在青海省4~10天降水預(yù)報(bào)中的應(yīng)用[J].蘭州大學(xué)學(xué)報(bào), 2010, 46(2): 33-38.

[7]羅陽(yáng),聶新旺,王廣山.幾種統(tǒng)計(jì)相似方法的適用性比較[J].氣象,2011,37(11):1443-1447.

[8]劉艷華,張少鋒,張治.利用歐洲數(shù)值預(yù)報(bào)產(chǎn)品制作哈爾濱市七日天氣預(yù)報(bào)[J].自然災(zāi)害學(xué)報(bào), 2013,20(6):185-190.

[9]黃紅兵,鐘偉,趙凱等.基于多元判據(jù)的短期天氣相似預(yù)報(bào)研究[J].空中交通,2014,9(S): 29-33.

黃紅兵（1965—），男，重慶忠縣人，教授級(jí)高級(jí)工程師，從事航空氣象研究。

作者簡(jiǎn)介：

基金項(xiàng)目：中國(guó)民用航空局空管局科研項(xiàng)目(民航空發(fā)2013[323]號(hào))。