台风台风森拉克生成

第19卷第3期2012年6月水土保持研究ResearchofSoilandWaterConservationVol.
19,No.
3Jun.
,2012收稿日期:2011-10-17修回日期:2011-12-04资助项目:温州市科技局重大项目"温州台风暴雨预测预报及灾情模型研究"(S20080017);浙江省教育厅项目"地貌形态多重分形特征与区域滑坡灾害空间分布关系研究"(20070527)作者简介:封毅(1978—),女,陕西宝鸡人,硕士,助理研究员,研究方向为台风预报及灾情评估.
E-mail:wanfe@126.
com通信作者:崔灵周(1971—),男,陕西西安人,博士,副教授,研究方向为环境灾害与GIS应用.
E-mail:clingzhou@126.
com基于BP神经网络的台风降雨量预测研究封毅1,武博强3,崔灵周2(1.
温州大学学报编辑部,浙江温州325035;2.
温州大学生命与环境科学学院,浙江温州325035;3.
长安大学地质工程学院,西安710054)摘要:台风登陆所带来的强降雨可引发严重的灾害.
依据BP神经网络基本原理,结合台风及降雨特征,建立了台风降雨量BP神经网络预测模型,利用200509号台风"麦莎"期间温州地区88个雨量观测站的连续6h雨量和台风特征参数,对该模型进行了训练和检验.
模型检验表明,6h降雨量预测相对误差小于30%的数据组所占比例达到71.
1%.
为台风期间降雨量时空分布快速预测提供了有益的尝试,对于沿海地区台风灾害预警、损失评估和应急决策具有重要现实意义.
关键词:BP神经网络;台风;降水量;预测中图分类号:P457.
8文献标识码:A文章编号:1005-3409(2012)03-0289-05StudyonForecastingTyphoonRainfallBasedonBPNeuralNetworkFENGYi1,WUBo-qiang3,CUILing-zhou2(1.
EditorialDepartmentofJournal,WenzhouUniversity,Wenzhou,Zhejiang325035,China;2.
CollegeofLifeandEnvironmentalScience,WenzhouUniversity,Wenzhou,Zhejiang325035,China;3.
CollogeofGeologicalEngineering,Chang′anUniversity,Xi′an710054,China)Abstract:Theheveayrainfallresultingfromtyphoonoftencausedseveredisasters.
BasedontheprincipleofBPneuralnetworkandrainfallcharacteroftyphoon,themodelofforecastingtyphoonrainfallwassetup.
Byusingtherainfallof6hoursfrom88rainfallobservationstationsduringtyphoonMatsa(No:200509)andpa-rameteroftyphooncharacter,themodelwastrainedandvalidated.
Theresultsshowedthattheproportionofrelativeerrorlowerthan30%ofper6hoursrainfallwas71.
1%.
Thestudytriedrapidpredictionofthetemporalvariationandspatialdistributionoftyphoonrainfall,whichwasmoreimportantfortyphoondisas-terforecast,thedamageassessmentanddecision-making.
Keywords:BPneuralnetwork;typhoon;rainfall;forecasting台风作为影响我国东南沿海的强烈天气系统,往往会引发持续性的特大暴雨降水过程,造成洪涝突发、农田受淹、城市内涝和路毁车阻等灾害,带来严重的经济和人员损失[1-2].
如何准确预报台风暴雨降水的空间分布,成为近年来倍受关注的重点领域之一.
黄永玉等[3]、钮学新[4]利用中尺度MM5模式分别对0418号"艾利"台风和0216号"森拉克"台风降水进行了数值模拟.
曾欣欣等[5]对0716号"罗莎"台风造成的浙江大暴雨到特大暴雨过程中,大环流天气形势演变、动力条件、水汽输送及物理量特征的影响进行了诊断研究.
有些学者[6-9]建立了适用于登陆台风的定量降水估计方法,实现了登陆台风未来0~3h的短时定量降水预报;利用多普勒雷达体扫和自动雨量站资料,结合概率配对法对"海棠"和"麦莎"两个台风的小时降水量进行了定量估测;采用多种物理量综合诊断分析方法,对登陆台风"云娜"在西进途中水汽来源和不稳定层结的维持等方面进行研究;用综合多级相似预报技术,对"浙东南沿海"台风过程所引发的暴雨天气进行研究.
邵月红[10]、邵利民等[11]分别利用BP神经网络方法,结合多普勒雷达体扫的回波强度及雨量计观测资料,估测了临沂地区暴雨过程的小时降雨量,对我国沿海的热带气旋的移动路径进行了预报.
可以看出,台风暴雨降水的研究主要集中在中尺度MM5模式、基于多普勒雷达体扫资料的概率配对和多级相似预报等方面,但由于模式分辨率较低和雷达探测范围有限等原因,台风暴雨降水预报的准确性和实际应用受到明显限制.
BP神经网络由于结构简单、自适应和自学习能力强等特点,在气象水文等领域得到了广泛应用[1,10-11].
本文以200509号台风"麦莎"在温州地区登陆所形成的降雨量为研究对象,利用雨量站实测降水数据,建立台风期间各雨量站6h降水量BP神经网络预测模型,实现台风暴雨期间重点区域降水量空间分布快速、准确估算,为防台减灾预案制定及综合决策提供科学依据.
1200509号台风"麦莎"概况1.
1台风路径及强度200509号"麦莎"台风于7月31日20时在菲律宾以东洋面上生成后向西北方向移动,8月2日8时加强为强热带风暴,并继续向西北方向移动,于8月3日2时加强为台风.
8月6日3时40分在温州地区东侧的玉环县干江镇登陆,登陆后继续沿西北方向移动进入温州东北的乐清和永嘉等地区,8月6日17时进入诸暨境内并减为强热带风暴.
1.
2台风中心气压、风力及风速从表1可以看出,200509号台风"麦莎"在登陆以前为超强台风,中心气压均保持在950hPa,最大风力和最大风速分别达到14级和45m/s.
台风登陆后,其中心气压逐步升高,即从登陆前的950hPa升高到960hPa,直至8月6号17时的980hPa;最大风力及风速均相应减小,分别从登陆前的14级和45m/s降低到13级和40m/s,至8月6号17点分别降低到11级和30m/s,此时的台风减弱为强热带风暴.
表1200509号台风"麦莎"中心气压、最大风力及风速变化表时间/月-日时间/时中心气压/hPa最大风力/级最大风速/(m·s-1)08-0514:00950144508-0520:00950144508-0604:00950144508-0605:00960134008-0608:00970123508-0611:00975123308-0614:00975123308-0617:009801130注:数据来源于温州气象局1.
3降水量时空变化图1—4是200509号台风"麦莎"登陆前8月5日8时至登陆后8月6日8时温州地区88个雨量站6h观测雨量空间分布图.
可见,温州地区的6h雨量随着台风逐步逼近和过境,其降水量呈现先增加后减小的变化趋势,其中6h降水量峰值出现在台风登陆后的8月6日3时至8时,北部地区的李家山站达到单站6h观测雨量最大值,为213mm;随着台风登陆各雨量站点的6h观测雨量快速下降.
从降雨量的空间分布来看,台风"麦莎"登陆前后温州地区的降水量主要分布于东北部的雁荡山区,中南部降水量显著偏少.
由于台风"麦莎"移动路径一直在温州地区东侧沿向西北方向行进,在玉环县登陆后继续沿西北方向穿过温州地区东北部的永嘉和乐清两地,台风路径变化成为导致降雨量在温州地区分布严重不均衡的主要原因.
图18月5日14:00-20:00温州地区降水量分布图28月5日21:00-8月6日2:00温州地区降水量分布092水土保持研究第19卷图38月6日3:00-8:00温州地区降水量分布图48月6日9:00-14:00温州地区降水量分布2台风降雨量预测模型建立BP神经网络也称误差反传(ErrorBackPropa-gation)前向网络,一般由输入层、输出层和若干隐层组成.
当信息输入网络时,先从输入层传至隐层节点,经激活函数作用和联接权重加权后,再传至下一隐层,然后经处理后由输出层产生计算结果.
将输出层所得到的计算结果与期望输出进行比较,若二者相差较大,则将误差进行反向传播,通过调整连接权值进行反复训练,直到误差达到预设的允许范围之内或训练次数;否则,停止训练.
完成训练后,该BP神经网络就可以进行预报等实际应用.
2.
1BP神经网络结构设计本文采用典型的三层BP神经网络结构(图5),由输入层、隐层和输出层构成,每个数据层包括若干数据节点.
输入数据层主要为台风降量影响因子层,具体包括雨量观测站点的经度(°)、纬度(°)、高程(m)、某时刻站点距台风中心距离(km)和台风中心气压5个输入数据节点,输出层仅包括相应时刻前6小时该站点的降雨量观测值(mm)1个输出数据节点,隐层的数据节点主要根据经验公式和网络训练效果进行调整,最后节点数确定为17个.
图5三层神经网络结构2.
2隐层激活函数和网络训练参数选取由于S型函数具有非线性放大系数功能,可把(-∞,+∞)变化范围的数据变换到(-1,+1)间输出,常被用作隐层的激活函数,本文BP神经网络隐层的激活函数选取Sigmoid函数,即:f(x)=ex-e-xex+e-x(1)该函数是连续可微的,便于误差反向传播过程节点权值的调节.
期望误差(训练目标)、学习速率、最大训练次数、初次权值和阈值等参数值设置对于BP网络训练和模型精度控制具有重要影响.
通过反复调整和比较,确定所建立的BP神经网络训练目标为0.
01、学习速率为0.
1、最大训练次数为1000、初始输入层至隐层及隐层至输出层的权值均设为0.
5、隐层节点及输出节点阈值分别设为0.
75和0.
2.
3输入层样本数据预处理BP神经网络输入层节点的数据由于量纲、数量级和单位等具有较大差别,在网络训练前须对原始数据进行归一化处理,以提高BP神经网络训练速度和灵敏性,同时可有效避开Sigmoid函数的饱和区.
本文采用的归一化公式如下,即:X′=2(X-Xmin)Xmax-Xmin-1(2)式中:X′———归一化后的数值;X———处理前的数值;Xmax,Xmin———该样本数据的最大值和最小值,归一化后的样本数据输出范围变为[-1,1].
2.
4BP神经网络的MATLAB实现MATLAB是美国MathWorks公司开发的商业数学软件,主要用于算法开发、数据可视化、数据分析以及数值计算.
该软件提供了神经网络工具箱,利用192第3期封毅等:基于BP神经网络的台风降雨量预测研究该工具箱,通过编写程序和相关函数调用,实现了该BP神经网络建立、网络训练参数设定与调整、输入层样本数据预处理、网络训练和模型应用.
3模型应用利用200509号台风"麦莎"期间温州地区的黄山、张溪、碧莲、上塘、潘山、曹村、埭头、泰顺等74个雨量站(见图6)连续6h观测总雨量和相应台风信息为训练数据,对本文所建立的台风降雨量BP神经网络预测模型进行训练,然后利用该场台风期间温州地区的雁荡、叶山、应坑、永嘉石柱、中保和朱涂等14个雨量站点(见图6)连续6h降雨总量进行了检验.
3.
1模型训练表2和表3分别是用于前文所建立的台风降雨量BP神经网络预测模型训练的部分雨量站(黄山站和碧莲站)及相关数据和模型训练结果误差分析.
从表2可见,模型训练的输入数据包括雨量站的经纬度、雨量站高程、初始时刻雨量站距台风中心距离及台风中心气压,输出数据为初始时刻后6h观测雨量.
受台风"麦莎"影响,黄山站和碧莲站均从8月5日3:00出现17.
2mm和19.
9mm的降雨,降雨分别持续至8月7日2:00和8月6日20:00.
以6h为时间单元,黄山站和碧莲站共有15组数据参与模型训练,其中黄山站8组数据、碧莲站7组数据.
全部参与模型训练的数据共74个雨量站的152组.
模型训练结果显示(见表3),不同雨量站各时段6h雨量预测相对误差均值为28.
4%、最大相对误差值为90%、相对误差小于30%的数据组所占比例达到75%,各雨量站总雨量、最大雨量预测相对误差和平均雨量预测相对误差分别为10.
2%,11.
7%,0.
9%.
本次模型训练精度较好.
图6BP模型训练和检验雨量站点分布表2台风降雨量BP神经网络预测模型训练部分数据表雨量站名称及经纬度/(°)日期(月-日)时间训练输入值训练输出值雨量站高程/m距台风中心距离/km台风中心气压/hPa6h降雨量/mm黄山(120.
7510E,28.
5654N)08-0503:00—08:0054049095017.
209:00—14:0054043495068.
615:00—20:0054034995053.
708-05—08-0621:00—02:00540216950131.
808-0603:00—08:00540111950168.
809:00—14:005403497095.
415:00—20:005405497521.
808-06—08-0721:00—02:005401919854.
0碧莲(120.
5830E,28.
3182N)08-0503:00—08:008048295019.
909:00—14:008042895036.
815:00—20:00803489508.
108-05—08-0621:00—02:008021795012.
608-0603:00—08:008011295023.
409:00—14:00803597042.
615:00—20:0080689755.
8表3台风降雨量BP神经网络预测模型训练结果及误差分析雨量站数/个数据量/组6h雨量相对误差均值/%6h雨量最大相对误差/%6h雨量相对误差≤30%比例/%单站总雨量相对误差均值/%单站最大雨量相对误差均值/%站点平均雨量相对误差/%7415228.
4907510.
211.
70.
93.
2模型检验由表4可见,模型检验输入输出数据同模型训练,即雨量站经纬度、雨量站高程、雨量站距台风中心距离、台风中心气压和6h降雨量.
雁荡山站从8月5日3:00至8月7日2:00,以6h为时间单元,共有8组数据参与模型检验;同样,叶山站从8月4日21:00至8292水土保持研究第19卷月6日14:00,以6h为时间单元,有9组数据参与模型检验;全部参与模型检验共计14个雨量站的56组数据.
模型检验结果表明(表5),各雨量站不同时段6h降雨量预测相对误差均值、最大相对误差和相对误差小于30%比例分别为38.
2%,95%,71.
1%,分别比模型训练提高9.
8%,5%和降低4.
9%.
各雨量站总雨量预测相对误差均值、最大雨量预测相对误差均值和平均雨量相对误差分别为29.
5%,33.
3%和29.
9%,分别比模型训练提高19.
3%,21.
6%和29%.
总体来看,模型检验结果可以接受.
表4台风降雨量BP神经网络预测模型检验部分数据表雨量站名称及经纬度/(°)日期(月-日)时间检验输入值检验输出值雨量站高程/m距台风中心距离/km台风中心气压/hPa6h降雨量/mm雁荡(121.
0570E,28.
3593N)08-0503:00—08:007045195014.
308-0509:00—14:007039495035.
708-0515:00—20:007030895051.
708-05—08-0621:00—02:007017595076.
408-0603:00—08:00707195055.
108-0609:00—14:00701997044.
508-0615:00—20:00709597560.
208-06—08-0721:00—02:00702309853.
5叶山(120.
9390E,28.
3608N)08-04—08-0521:00—02:0049056195012.
408-0503:00—08:0049045995024.
009:00—14:0049040395046.
615:00—20:0049031995061.
708-05—08-0621:00—02:00490185950113.
008-0603:00—08:0049081950115.
79:00—14:00490897069.
215:00—20:004908597515.
808-06—08-0721:00—02:004902229853.
8表5台风降雨量BP神经网络预测模型检验结果及误差分析雨量站数/个数据量/组6h雨量相对误差均值/%6h雨量最大相对误差/%6h雨量相对误差≤30%比例/%单站总雨量相对误差均值/%单站最大雨量相对误差均值/%站点平均雨量相对误差/%145638.
29571.
129.
533.
329.
94结论(1)基于BP神经网络的台风降雨量预测模型可通过台风中心距离、气压、空间点经纬度及高程信息快速实现台风影响区域6h降雨量空间分布预测,该模型经过200509号台风"麦莎"期间温州地区14个雨量站的56组数据的检验,6h降雨量预测相对误差小于30%的比例达到71.
1%.
(2)模型检验仅采用了同场次台风未参与模型训练的14个雨量站56组数据,使模型的适用性受到一定影响;同时,由于数据限制,本模型没有考虑局地环流特征对台风降雨空间分布的影响,导致模型训练及检验精度不高.
参考文献:[1]程正泉,陈联寿,徐祥德,等.
近年中国台风暴雨研究进展[J].
气象,2005,31(12):3-9.
[2]李江南,龚志鹏,王安宇,等.
近十年来台风暴雨研究的若干进展与讨论[J].
热带地理,2004,42(2):113-117.
[3]黄永玉,沈桐立,沈新勇,等.
0418号"艾利"台风暴雨过程的数值模拟[J].
台湾海峡,2006,25(1):102-109.
[4]钮学新,杜惠良,刘建勇.
0216号台风降水及其影响降水机制的数值模拟试验[J].
气象学报,2005,63(1):58-69.
[5]曾欣欣,黄新晴,滕代高.
罗莎台风造成浙江特大暴雨的过程分析[J].
海洋学研究,2010,28(1):63-72.
[6]岳彩军,陈佩燕,雷小途,等.
登陆台风短时定量降水预报方法初探[J].
气象科技,2006,34(1):7-11.
[7]冀春晓,陈联寿,徐祥德,等.
多普勒雷达资料动态定量估测台风小时降水量的研究[J].
热带气象学报,2008,42(2):147-155.
[8]施望芝,毛以伟,王建生,等.
台风降水云区中单站强降水诊断分析和预报[J].
气象科学,2006,26(6):668-75.
[9]李博,赵思雄.
用SMAT建立台风暴雨预报模型的试验研究[J].
气象,2009,35(6):3-12.
[10]邵月红,张万昌,刘永和,等.
BP神经网络在多普勒雷达降水量的估测中的应用[J].
高原气象,2009,28(4):846-853.
[11]邵利民,傅刚,曹祥村,等.
BP神经网络在台风路径预报中的应用[J].
自然灾害学报,2009,18(6):104-111.
392第3期封毅等:基于BP神经网络的台风降雨量预测研究