貴志成 王艷瑩 趙小平 甘思舊

（1 63796部隊(duì)，西昌 615000；2 民航海南空管分局，?？?570100）

0 引言

淺層風(fēng)是近地面層中的風(fēng)，近地面層是大氣邊界層最接近地表面的部分，約為距離地表10～100 m的大氣。航天發(fā)射任務(wù)中，火箭在垂直轉(zhuǎn)運(yùn)和加注發(fā)射階段，遇到強(qiáng)淺層風(fēng)天氣，可能會(huì)影響航天發(fā)射任務(wù)成敗，因此淺層風(fēng)預(yù)報(bào)一直是航天氣象保障的重點(diǎn)。熱帶氣旋是最具破壞力的自然災(zāi)害之一[1]，文昌航天發(fā)射場(chǎng)（以下簡(jiǎn)稱(chēng)場(chǎng)區(qū)）位于海南島東北海岸，年平均受到3.58個(gè)熱帶氣旋影響，大多數(shù)臺(tái)風(fēng)都是從東南方向移近海南島，海南島東部海岸受臺(tái)風(fēng)影響最嚴(yán)重，受到臺(tái)風(fēng)影響時(shí)間也較長(zhǎng)[2]。而出現(xiàn)強(qiáng)淺層風(fēng)天氣個(gè)例中，熱帶氣旋造成的淺層風(fēng)風(fēng)速大、大風(fēng)影響時(shí)間長(zhǎng)、風(fēng)速變化規(guī)律較復(fù)雜。夏秋季場(chǎng)區(qū)頻繁受到臺(tái)風(fēng)影響，受臺(tái)風(fēng)影響時(shí)淺層風(fēng)變化特征一直是淺層風(fēng)預(yù)報(bào)研究的重點(diǎn)和難點(diǎn)。日常淺層風(fēng)預(yù)報(bào)保障中，臺(tái)風(fēng)影響文昌場(chǎng)區(qū)時(shí)淺層風(fēng)風(fēng)速的定量化研究較不足，臺(tái)風(fēng)開(kāi)始影響的具體時(shí)間也較難把握，如2016年臺(tái)風(fēng)莎莉嘉登陸萬(wàn)寧時(shí)，24 h風(fēng)速預(yù)報(bào)最大偏差為14.5 m/s。場(chǎng)區(qū)受臺(tái)風(fēng)影響時(shí)，淺層風(fēng)強(qiáng)風(fēng)是造成設(shè)備設(shè)施的損壞、危害人員安全的主要因素，提高淺層風(fēng)預(yù)報(bào)準(zhǔn)確率可以為航天發(fā)射任務(wù)氣象保障提供支撐。

目前對(duì)臺(tái)風(fēng)機(jī)理研究和臺(tái)風(fēng)引發(fā)區(qū)域大風(fēng)研究的主要方法是基于動(dòng)力學(xué)的數(shù)值模擬分析[3-5]。孫學(xué)金等[6]利用WRF-LES模式對(duì)大氣邊界層內(nèi)大尺度湍渦進(jìn)行解析，提高了干旱湖區(qū)近地面風(fēng)場(chǎng)數(shù)值模擬的真實(shí)性和準(zhǔn)確性。馮慧敏等[7]利用數(shù)值預(yù)報(bào)資料進(jìn)行鄭州市地面風(fēng)的多模式預(yù)報(bào)試驗(yàn)，發(fā)現(xiàn)在數(shù)值預(yù)報(bào)的基礎(chǔ)上，通過(guò)統(tǒng)計(jì)降尺度預(yù)報(bào)能進(jìn)一步提高地面風(fēng)預(yù)報(bào)準(zhǔn)確率。

對(duì)基層氣象預(yù)報(bào)員來(lái)說(shuō)，統(tǒng)計(jì)學(xué)方法在計(jì)算效率和應(yīng)用層面上較數(shù)值模式更具優(yōu)勢(shì)，在數(shù)據(jù)趨勢(shì)描述方面，線(xiàn)性回歸具有直觀性和簡(jiǎn)潔性。黃奕武等[8]利用T639數(shù)據(jù)線(xiàn)性回歸方法，建立了臺(tái)風(fēng)中心預(yù)報(bào)誤差與環(huán)境風(fēng)整層垂直切變的24～120 h線(xiàn)性預(yù)估模型，且預(yù)估模型有較好的預(yù)估效果。農(nóng)吉夫[9]也通過(guò)條件數(shù)分析方法，建立了西北太平洋臺(tái)風(fēng)強(qiáng)度多元回歸模型?？蒲袉挝恢饕芯恐攸c(diǎn)為臺(tái)風(fēng)發(fā)生發(fā)展機(jī)理和臺(tái)風(fēng)結(jié)構(gòu)分析，地方業(yè)務(wù)單位主要進(jìn)行地面風(fēng)預(yù)報(bào)，對(duì)淺層風(fēng)的預(yù)報(bào)技術(shù)研究較少。風(fēng)力發(fā)電站有淺層風(fēng)精確預(yù)報(bào)的需求，但是大多是向氣象公司尋求商業(yè)服務(wù)，利用統(tǒng)計(jì)學(xué)方法進(jìn)行臺(tái)風(fēng)影響時(shí)淺層風(fēng)特征及其預(yù)報(bào)技術(shù)的研究相對(duì)較少，本文針對(duì)業(yè)務(wù)工作中的實(shí)際需求，以提高預(yù)報(bào)準(zhǔn)確率為目的，通過(guò)統(tǒng)計(jì)分析的方法研究臺(tái)風(fēng)影響場(chǎng)區(qū)時(shí)淺層風(fēng)預(yù)報(bào)偏差的影響因子，進(jìn)行因子設(shè)計(jì)試驗(yàn)，建立受臺(tái)風(fēng)影響時(shí)場(chǎng)區(qū)淺層風(fēng)24 h預(yù)報(bào)回歸方程。

1 數(shù)據(jù)來(lái)源及統(tǒng)計(jì)范圍

淺層風(fēng)實(shí)況數(shù)據(jù)是架設(shè)在文昌場(chǎng)區(qū)內(nèi)的三個(gè)淺層風(fēng)測(cè)風(fēng)塔探測(cè)2015—2017年50 m高度10 min平均風(fēng)向風(fēng)速。資料統(tǒng)計(jì)范圍為2015年6月—2017年10月影響場(chǎng)區(qū)的9個(gè)臺(tái)風(fēng)個(gè)例。

場(chǎng)區(qū)淺層風(fēng)受臺(tái)風(fēng)影響日的定義：

1）臺(tái)風(fēng)開(kāi)始影響日：淺層風(fēng)風(fēng)速比上一日平均風(fēng)速增大4 m/s或以上。

2）臺(tái)風(fēng)影響結(jié)束日：當(dāng)日淺層風(fēng)平均風(fēng)速比前一日淺層風(fēng)平均風(fēng)速減小4 m/s且最大風(fēng)速小于10 m/s，則前一日為臺(tái)風(fēng)影響結(jié)束日。

2 淺層風(fēng)預(yù)報(bào)偏差分析

2.1 預(yù)報(bào)概況

場(chǎng)區(qū)24 h淺層風(fēng)預(yù)報(bào)共分為四段，分別為18：00—24：00、00：00—06：00、06：00—12：00、12：00—18：00。預(yù)報(bào)每個(gè)時(shí)段50 m高度最大風(fēng)速和最大風(fēng)對(duì)應(yīng)風(fēng)向。根據(jù)臺(tái)風(fēng)影響場(chǎng)區(qū)的定義，共篩選出受臺(tái)風(fēng)影響日數(shù)26 d，每天數(shù)據(jù)分為4段，共整理了臺(tái)風(fēng)影響場(chǎng)區(qū)共26 d（104段）實(shí)況和預(yù)報(bào)數(shù)據(jù)。26 d淺層風(fēng)24 h預(yù)報(bào)準(zhǔn)確率為67 %，離淺層風(fēng)精準(zhǔn)精細(xì)預(yù)報(bào)的要求還有一定差距。

2.2 預(yù)報(bào)偏差的過(guò)程能力分析

臺(tái)風(fēng)影響期間，由于淺層風(fēng)風(fēng)速較大，風(fēng)向變化相對(duì)連續(xù)[10]，淺層風(fēng)風(fēng)向預(yù)報(bào)準(zhǔn)確率高，按照風(fēng)向的評(píng)分標(biāo)準(zhǔn)基本不會(huì)出現(xiàn)風(fēng)向預(yù)報(bào)錯(cuò)誤的情況，因此本文主要對(duì)風(fēng)速預(yù)報(bào)偏差進(jìn)行分析。

通過(guò)正態(tài)檢驗(yàn)發(fā)現(xiàn)，淺層風(fēng)預(yù)報(bào)偏差服從正態(tài)分布，且各數(shù)據(jù)間相互獨(dú)立，可以進(jìn)行過(guò)程能力分析。過(guò)程能力指數(shù)反映過(guò)程在受控狀態(tài)下時(shí)，客戶(hù)要求與過(guò)程表現(xiàn)（產(chǎn)品質(zhì)量或服務(wù)質(zhì)量的變動(dòng)程度）的比值，如果過(guò)程表現(xiàn)越能滿(mǎn)足顧客要求，則過(guò)程能力指數(shù)值越大[11]。

根據(jù)淺層風(fēng)預(yù)報(bào)評(píng)定標(biāo)準(zhǔn)，設(shè)定目標(biāo)為淺層風(fēng)預(yù)報(bào)準(zhǔn)確率大于80 %，假設(shè)24 h 4段淺層風(fēng)預(yù)報(bào)和實(shí)況風(fēng)速偏差均小于等于3 m/s時(shí)，預(yù)報(bào)準(zhǔn)確率為84 %（在風(fēng)向全部正確的情況下），滿(mǎn)足目標(biāo)要求。因此將上下限規(guī)定為±3 m/s，即USL=3 m/s，LSL=?3 m/s。

通過(guò)計(jì)算得到各段淺層風(fēng)預(yù)報(bào)偏差的能力直方圖（圖1），其中，Cp值反映過(guò)程的潛在能力，Cp值越大，則過(guò)程受控程度越高，Cpk反映過(guò)程的實(shí)際表現(xiàn)，能反映出工序偏離規(guī)格中心的問(wèn)題，Cpk值越大，則過(guò)程偏移量越小。

淺層風(fēng)預(yù)報(bào)平均絕對(duì)偏差為3.46 m/s，預(yù)報(bào)偏差組內(nèi)標(biāo)準(zhǔn)差為4.72，平均偏差Cp值為0.21，Cpk值為0.13，4段淺層風(fēng)預(yù)報(bào)偏差都存在正偏倚情況，即實(shí)況風(fēng)速大于預(yù)報(bào)風(fēng)速出現(xiàn)次數(shù)較多。

3 臺(tái)風(fēng)關(guān)鍵因子的識(shí)別

從臺(tái)風(fēng)特點(diǎn)出發(fā)，考慮臺(tái)風(fēng)強(qiáng)度、臺(tái)風(fēng)相對(duì)場(chǎng)區(qū)位置和臺(tái)風(fēng)與場(chǎng)區(qū)距離三個(gè)影響因素與預(yù)報(bào)偏差的關(guān)系[12]，通過(guò)統(tǒng)計(jì)發(fā)現(xiàn)，臺(tái)風(fēng)與場(chǎng)區(qū)距離對(duì)預(yù)報(bào)偏差有正相關(guān)關(guān)系，即臺(tái)風(fēng)距離場(chǎng)區(qū)越近，越容易出現(xiàn)預(yù)報(bào)偏差。接下來(lái)將對(duì)臺(tái)風(fēng)強(qiáng)度、臺(tái)風(fēng)相對(duì)場(chǎng)區(qū)位置兩個(gè)要素與預(yù)報(bào)偏差的關(guān)系進(jìn)行分析。

3.1 臺(tái)風(fēng)強(qiáng)度與預(yù)報(bào)偏差的分析

從影響場(chǎng)區(qū)淺層風(fēng)臺(tái)風(fēng)強(qiáng)度對(duì)預(yù)報(bào)偏差的卡方檢驗(yàn)來(lái)看（圖2），熱帶風(fēng)暴級(jí)出現(xiàn)的頻數(shù)最高，且對(duì)預(yù)報(bào)偏差的貢獻(xiàn)率最大，其次為強(qiáng)熱帶風(fēng)暴級(jí)影響場(chǎng)區(qū)時(shí)容易出現(xiàn)淺層風(fēng)預(yù)報(bào)偏差。熱帶低壓級(jí)雖然出現(xiàn)頻數(shù)第二，但是對(duì)預(yù)報(bào)偏差的貢獻(xiàn)率較小，即熱帶低壓影響場(chǎng)區(qū)時(shí)淺層風(fēng)預(yù)報(bào)偏差較小。總的來(lái)說(shuō)不同強(qiáng)度臺(tái)風(fēng)對(duì)預(yù)報(bào)偏差有顯著差異，但二者并不是呈正/負(fù)相關(guān)關(guān)系。

圖1 淺層風(fēng)預(yù)報(bào)偏差的過(guò)程能力分析Fig. 1 Analysis of process forecast capability based on shallow layer wind forecast deviations

3.2 臺(tái)風(fēng)方位與預(yù)報(bào)偏差的分析

影響海南島臺(tái)風(fēng)生成源地為南海和西太平洋，臺(tái)風(fēng)主要從海上向西或西北方向移動(dòng)影響場(chǎng)區(qū)[13]，不同方位的臺(tái)風(fēng)對(duì)場(chǎng)區(qū)淺層風(fēng)的影響有一定差別，為了將不同方向移近并影響場(chǎng)區(qū)的臺(tái)風(fēng)進(jìn)行區(qū)分，對(duì)受臺(tái)風(fēng)影響時(shí)淺層風(fēng)風(fēng)向進(jìn)行劃分，將水平方向風(fēng)向劃分為八個(gè)象限，角度分別為（0°≤第一象限＜45°、45°≤第二象限＜90°、90°≤第三象限＜135°、135°≤第四象限＜180°、180°≤第五象限＜225°、225°≤第六象限＜270°、270°≤第七象限＜315°、315°≤第八象限＜360°），按照八個(gè)象限對(duì)淺層風(fēng)風(fēng)向進(jìn)行劃分，不同的象限反映出臺(tái)風(fēng)相對(duì)于場(chǎng)區(qū)的方位。通過(guò)統(tǒng)計(jì)風(fēng)向和預(yù)報(bào)偏差的卡方擬合度檢驗(yàn)（圖3），發(fā)現(xiàn)場(chǎng)區(qū)受到臺(tái)風(fēng)影響的風(fēng)向主要在第一至第六象限，未出現(xiàn)第七、第八象限風(fēng)向。場(chǎng)區(qū)淺層風(fēng)風(fēng)向?yàn)槠摺讼笙迺r(shí)，臺(tái)風(fēng)位于場(chǎng)區(qū)西北方向，很有可能已經(jīng)登陸大陸地區(qū)，臺(tái)風(fēng)登陸后強(qiáng)度減弱明顯，大風(fēng)區(qū)影響范圍也顯著減小，因此對(duì)場(chǎng)區(qū)影響也不明顯。從圖中分析可以發(fā)現(xiàn)，風(fēng)向?yàn)榈谝?、第三象限時(shí)對(duì)淺層風(fēng)偏差貢獻(xiàn)率最大，且出現(xiàn)頻率也較高，第六象限卡方值的貢獻(xiàn)較大，但主要為負(fù)貢獻(xiàn)。臺(tái)風(fēng)影響場(chǎng)區(qū)時(shí)，南風(fēng)（包括西南風(fēng)、南風(fēng)、東南風(fēng)）比北風(fēng)（包括西北風(fēng)、北風(fēng)、東北風(fēng)）出現(xiàn)偏差的可能性更高，淺層風(fēng)為北風(fēng)時(shí)臺(tái)風(fēng)位于場(chǎng)區(qū)東側(cè)海上，淺層風(fēng)為南風(fēng)時(shí)，臺(tái)風(fēng)位于場(chǎng)區(qū)西側(cè)，很有可能已經(jīng)登陸海南島，因此臺(tái)風(fēng)中心距離場(chǎng)區(qū)位置相對(duì)于北風(fēng)時(shí)更近，風(fēng)速變化幅度更大，所以也更容易出現(xiàn)淺層風(fēng)預(yù)報(bào)偏差。

4 臺(tái)風(fēng)影響時(shí)淺層風(fēng)預(yù)報(bào)方程建立與檢驗(yàn)

圖3 影響場(chǎng)區(qū)淺層風(fēng)預(yù)報(bào)偏差和風(fēng)向的卡方擬合度檢驗(yàn)Fig. 3 The goodness of Chi-Square fitting test between shallow wind forecast deviation and typhoon wind direction

臺(tái)風(fēng)的預(yù)報(bào)主要依靠數(shù)值預(yù)報(bào)模式，數(shù)值預(yù)報(bào)模式對(duì)未來(lái)24 h內(nèi)臺(tái)風(fēng)強(qiáng)度的預(yù)報(bào)準(zhǔn)確率可達(dá)80 %左右[14]，24 h內(nèi)臺(tái)風(fēng)路徑預(yù)報(bào)的偏差不超過(guò)100 km。淺層風(fēng)預(yù)報(bào)中主要參考場(chǎng)區(qū)附近數(shù)值預(yù)報(bào)低層風(fēng)速，在研究中發(fā)現(xiàn)，建立回歸方程時(shí)選取925 hPa高度風(fēng)場(chǎng)作為回歸方程自變量有較好的顯著性。因此在建立回歸方程時(shí)考慮四個(gè)自變量，分別為925 hPa數(shù)值預(yù)報(bào)風(fēng)速、臺(tái)風(fēng)強(qiáng)度、臺(tái)風(fēng)象限、臺(tái)風(fēng)相對(duì)場(chǎng)區(qū)距離。

由于數(shù)值預(yù)報(bào)產(chǎn)品中925 hPa高度風(fēng)場(chǎng)分辨率為0.25°×0.25°，場(chǎng)區(qū)位置沒(méi)有處于數(shù)值預(yù)報(bào)網(wǎng)格點(diǎn)中，因此選取場(chǎng)區(qū)周?chē)?個(gè)網(wǎng)格點(diǎn)進(jìn)行數(shù)據(jù)分析（圖4），4個(gè)點(diǎn)的經(jīng)緯度分別為場(chǎng)區(qū)左上（19.75°N，110.75°E）、左下（19.50°N，110.75°E）、右上（19.75°N，111.00°E）、右下（19.50°N，111.00°E），下文分別以a、b、c、d代表4個(gè)網(wǎng)格點(diǎn)，對(duì)4點(diǎn)數(shù)據(jù)進(jìn)行單因子試驗(yàn)設(shè)計(jì)。

4.1 數(shù)值預(yù)報(bào)風(fēng)速的單因子檢驗(yàn)

分別對(duì)選取的4點(diǎn)數(shù)據(jù)進(jìn)行一天4段的相關(guān)性分析（表1），從表中可以看到，a、b、c、d點(diǎn)和四點(diǎn)平均與實(shí)況的相關(guān)系數(shù)均在0.6以上，且均通過(guò)水平為95 %的可信度檢驗(yàn)，縱向來(lái)看，數(shù)值預(yù)報(bào)與第一段實(shí)況的相關(guān)系數(shù)最大為0.81，與第四段實(shí)況的相關(guān)系數(shù)最小為0.68，相關(guān)系數(shù)由第一段實(shí)況向第四段實(shí)況逐漸遞減。橫向比較a、b、c、d點(diǎn)和四點(diǎn)平均與實(shí)況的相關(guān)系數(shù)，發(fā)現(xiàn)各點(diǎn)之間相關(guān)系數(shù)差異較小，a點(diǎn)和四點(diǎn)平均與實(shí)況的相關(guān)系數(shù)略大。

圖4 場(chǎng)區(qū)位置及網(wǎng)格點(diǎn)的選取Fig. 4 Location of the site and grid selection

表1 數(shù)值預(yù)報(bào)與實(shí)況序列的相關(guān)系數(shù)Table 1 Correlation coefficient between numerical prediction and observed data

分析各點(diǎn)數(shù)據(jù)與實(shí)況風(fēng)速偏差的基本統(tǒng)計(jì)量（表2），可以看到，b點(diǎn)的平均偏差值最小，四點(diǎn)平均偏差數(shù)據(jù)的標(biāo)準(zhǔn)差、標(biāo)準(zhǔn)誤差最小。綜合考慮相關(guān)系數(shù)與基本統(tǒng)計(jì)量結(jié)果，選取相關(guān)系數(shù)最大，且偏差的標(biāo)準(zhǔn)差、標(biāo)準(zhǔn)誤差最小的四點(diǎn)平均數(shù)據(jù)作為擬合方程中的基本數(shù)據(jù)。

表2 各點(diǎn)數(shù)值預(yù)報(bào)值與實(shí)況偏差的基本統(tǒng)計(jì)量 Table 2 Basic statistics of numerical prediction and observed data deviation in different grids

通過(guò)相關(guān)性檢驗(yàn)可以發(fā)現(xiàn)，數(shù)值預(yù)報(bào)風(fēng)速與實(shí)況風(fēng)速有很高的相關(guān)性，而臺(tái)風(fēng)強(qiáng)度、臺(tái)風(fēng)方位、臺(tái)風(fēng)距離與實(shí)況風(fēng)速的相關(guān)性均不超過(guò)±0.3，因此在因子設(shè)計(jì)中主要分析這三個(gè)要因與實(shí)況風(fēng)速的檢驗(yàn)。

4.2 因子設(shè)計(jì)分析

為了確定因子對(duì)響應(yīng)變量的影響，確定關(guān)鍵因素的適當(dāng)設(shè)置以取得響應(yīng)變量的變異最小化。對(duì)925 hPa數(shù)值預(yù)報(bào)風(fēng)速、臺(tái)風(fēng)強(qiáng)度、臺(tái)風(fēng)方位（象限）、臺(tái)風(fēng)與場(chǎng)區(qū)距離四個(gè)因子進(jìn)行因子設(shè)計(jì)?？紤]到進(jìn)行因子設(shè)計(jì)、選定回歸模型后要進(jìn)行模型效果檢驗(yàn)，因此在進(jìn)行因子分析時(shí)選取2015—2016年場(chǎng)區(qū)受臺(tái)風(fēng)影響的17 d×4段共68段數(shù)據(jù)，在模型效果檢驗(yàn)階段選擇2017年5次臺(tái)風(fēng)的9 d×4段共36段數(shù)據(jù)進(jìn)行效果檢驗(yàn)。

因子設(shè)計(jì)階段共設(shè)計(jì)三組對(duì)照試驗(yàn)（表3），通過(guò)因子逐步剔除的方法，將回歸方程中不顯著項(xiàng)剔除，方程中自變量能通過(guò)顯著性檢驗(yàn)，則方程的穩(wěn)定性較高，因此試驗(yàn)1至試驗(yàn)3的因子數(shù)量逐漸減少。在建立回歸方程過(guò)程中，首先考慮回歸方程的穩(wěn)定性，其次考慮回歸方程擬合度。

對(duì)比三組試驗(yàn)回歸效果（表4），其中S為殘差的標(biāo)準(zhǔn)差，R-Sq為回歸模型誤差占總誤差的百分比，取值范圍為0～100%，數(shù)值越大，表明回歸模型數(shù)據(jù)溫和越好。R-Sq（調(diào)整）為調(diào)整后的R-Sq，與R-Sq數(shù)值越接近，表明回歸模型越可靠。各試驗(yàn)間S、R-Sq和R-Sq（調(diào)整）差距不大，但從因子的顯著性來(lái)看，試驗(yàn)3中各因子都能通過(guò)顯著性檢驗(yàn)，而且在回歸效果差距不大的情況下，選擇更少的因子可以節(jié)約計(jì)算量。

表3 試驗(yàn)設(shè)計(jì)及各試驗(yàn)回歸方程Table 3 Tests designs and the regression equations

表4 各試驗(yàn)效果對(duì)比Table 4 Comparison of test results

4.3 試驗(yàn)3可靠性分析

通過(guò)逐步剔除法找到試驗(yàn)3為最理想回歸方程，對(duì)試驗(yàn)3回歸方程的可靠性進(jìn)行分析，回歸方程的顯著性檢驗(yàn)P值為0.0，可以通過(guò)可信度為99 %的顯著性檢驗(yàn)。從回歸方程的殘差圖中可以看到（圖5），圖5a為殘差正態(tài)性概率圖，累計(jì)的概率點(diǎn)列沿對(duì)角線(xiàn)分布，則說(shuō)明殘差服從正態(tài)分布，且概率點(diǎn)越接近對(duì)角線(xiàn)分布，說(shuō)明回歸效果越好。圖5b為響應(yīng)變量對(duì)應(yīng)的殘差圖，殘差點(diǎn)在橫軸上下隨機(jī)波動(dòng)，即不存在某種連續(xù)性的波動(dòng)趨勢(shì)。圖5c為殘差的正態(tài)分布圖，可以看到殘差分布集中的區(qū)域?yàn)?4～4。圖5d為按觀測(cè)值順序的殘差分布，殘差分布為隨機(jī)波動(dòng)。通過(guò)對(duì)殘差圖的分析，認(rèn)為試驗(yàn)3回歸方程（Y=12.2+0.523×風(fēng)速?0.0133×距離）穩(wěn)定性較好，且方程結(jié)果不發(fā)散，是逐步剔除試驗(yàn)中的最優(yōu)選擇。

4.4 回歸方程效果檢驗(yàn)

將2017年9d 36段風(fēng)速數(shù)據(jù)代入試驗(yàn)3回歸方程Y=12.2+0.523×風(fēng)速?0.0133×距離，算出預(yù)報(bào)風(fēng)速后與實(shí)況對(duì)比，對(duì)偏差進(jìn)行過(guò)程能力分析（圖6），偏差的組內(nèi)標(biāo)準(zhǔn)差為3.74，整體標(biāo)準(zhǔn)差為3.68，Cp值為0.27，Cpk值為0.19。

將試驗(yàn)3的回歸驗(yàn)證結(jié)果與淺層風(fēng)預(yù)報(bào)偏差的過(guò)程能力分析進(jìn)行對(duì)比（表5），回歸方程得到預(yù)報(bào)風(fēng)速有8次偏差大于等于4 m/s。改進(jìn)前預(yù)報(bào)偏差有12次風(fēng)速大于等于4 m/s。建立回歸方程后得到的回歸風(fēng)速標(biāo)準(zhǔn)差比改進(jìn)前減小21 %，Cp值提高28 %，Cpk值提高46 %，說(shuō)明運(yùn)用回歸方程的預(yù)報(bào)準(zhǔn)確率有提升、穩(wěn)定性更好。

5 結(jié)論

本文通過(guò)分析2015—2017年臺(tái)風(fēng)影響場(chǎng)區(qū)的26 d淺層風(fēng)實(shí)況、24 h預(yù)報(bào)數(shù)據(jù)，分析預(yù)報(bào)偏差與臺(tái)風(fēng)強(qiáng)度、臺(tái)風(fēng)與場(chǎng)區(qū)方位的關(guān)系，找出建立回歸方程的4個(gè)因子（數(shù)值預(yù)報(bào)925 hPa風(fēng)速、臺(tái)風(fēng)強(qiáng)度、臺(tái)風(fēng)方位、臺(tái)風(fēng)與場(chǎng)區(qū)距離），通過(guò)因子試驗(yàn)設(shè)計(jì)，找到顯著性和穩(wěn)定性最高的淺層風(fēng)預(yù)報(bào)回歸方程，并對(duì)最優(yōu)回歸方程進(jìn)行數(shù)據(jù)檢驗(yàn)，通過(guò)研究發(fā)現(xiàn)：

1）不同臺(tái)風(fēng)的強(qiáng)度、臺(tái)風(fēng)的方位對(duì)淺層風(fēng)預(yù)報(bào)偏差影響不同，這種影響為非線(xiàn)性關(guān)系。

2）場(chǎng)區(qū)周?chē)?點(diǎn)算數(shù)平均的數(shù)值預(yù)報(bào)925 hPa風(fēng)場(chǎng)數(shù)據(jù)與實(shí)況偏差最小，最適合作為建立回歸方程的自變量之一。

3）利用2015—2016年72段數(shù)據(jù)，通過(guò)因子分析試驗(yàn)發(fā)現(xiàn)，臺(tái)風(fēng)強(qiáng)度、臺(tái)風(fēng)方位兩個(gè)因子在回歸方程中顯著性不明顯，選取顯著性、穩(wěn)定性最高，且標(biāo)準(zhǔn)差最小的試驗(yàn)3方程作為淺層風(fēng)預(yù)報(bào)回歸方程。

4）利用2017年36段數(shù)據(jù)對(duì)試驗(yàn)3回歸方程效果進(jìn)行檢驗(yàn)，發(fā)現(xiàn)回歸方程能有效減小預(yù)報(bào)偏差，且減少較大誤差出現(xiàn)次數(shù)，回歸效果較好。

圖5 試驗(yàn)3風(fēng)速的殘差四合一圖Fig. 5 Four-in-one diagram of wind speed residuals in test 3

圖6 試驗(yàn)3回歸方程的效果檢驗(yàn)Fig. 6 The verification of regression equation in test 3

表5 試驗(yàn)效果對(duì)比Table 5 Comparison of verified results

通過(guò)建立風(fēng)速預(yù)報(bào)回歸方程，可以有效提高24 h預(yù)報(bào)準(zhǔn)確率，但由于場(chǎng)區(qū)受臺(tái)風(fēng)影響時(shí)淺層風(fēng)數(shù)據(jù)積累不夠多，還需進(jìn)行資料收集，通過(guò)增大樣本數(shù)據(jù)個(gè)數(shù)改善回歸方程。后續(xù)研究中還應(yīng)進(jìn)行回歸方程效果評(píng)定，針對(duì)臺(tái)風(fēng)影響后期風(fēng)速明顯減小時(shí)擬合效果不好的情況進(jìn)行改進(jìn)提高。