CN103995951A - Typhoon key parameter extraction method based on half-normal model - Google Patents

Typhoon key parameter extraction method based on half-normal model Download PDF

Info

Publication number
CN103995951A
CN103995951A CN201410016770.5A CN201410016770A CN103995951A CN 103995951 A CN103995951 A CN 103995951A CN 201410016770 A CN201410016770 A CN 201410016770A CN 103995951 A CN103995951 A CN 103995951A
Authority
CN
China
Prior art keywords
typhoon
normal model
sample
tracks
pmax
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201410016770.5A
Other languages
Chinese (zh)
Inventor
汪穗峰
陈碧云
何平
陈玫丽
陈绍南
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangxi University
Yangjiang Power Supply Bureau of Guangdong Power Grid Co Ltd
Original Assignee
Guangxi University
Yangjiang Power Supply Bureau of Guangdong Power Grid Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Guangxi University, Yangjiang Power Supply Bureau of Guangdong Power Grid Co Ltd filed Critical Guangxi University
Priority to CN201410016770.5A priority Critical patent/CN103995951A/en
Publication of CN103995951A publication Critical patent/CN103995951A/en
Pending legal-status Critical Current

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于半正态模型的台风关键参数提取方法,本发明将台风历史数据输入到计算机中,计算机选取与台风高发区域相对应的地图数据中的任一点为模拟点,做出模拟圆,从而得到经过模拟圆的台风关键参数;计算半正态模型特征参数,得到Al的台风关键参数的半正态模型特征参数,AR的台风关键参数的半正态模型特征参数。本发明具有可以准确提取台风的关键参数,为进一步研究台风强度、台风路径和台风对人民生活影响奠定了可靠基础;为台风风能的利用提供可靠数据基础的特点。

The invention discloses a method for extracting key parameters of a typhoon based on a half-normal model. The invention inputs typhoon historical data into a computer, and the computer selects any point in the map data corresponding to a typhoon high-incidence area as a simulation point, and draws Simulate the circle to obtain the key parameters of the typhoon passing through the simulated circle; calculate the characteristic parameters of the half-normal model to obtain the characteristic parameters of the half-normal model of the key parameters of the typhoon of Al and the characteristic parameters of the half-normal model of the key parameters of the typhoon of A R. The invention has the characteristics of being able to accurately extract key parameters of typhoons, laying a reliable foundation for further research on typhoon intensity, typhoon path and impact of typhoon on people's lives, and providing reliable data basis for utilization of typhoon wind energy.

Description

基于半正态模型的台风关键参数提取方法Extraction method of key parameters of typhoon based on half-normal model

技术领域technical field

本发明涉及台风研究技术领域,尤其是涉及一种可以准确提取台风的关健参数,为进一步研究台风强度、台风路径和台风对人民生活影响奠定了可靠基础的基于半正态模型的台风关键参数提取方法。The present invention relates to the technical field of typhoon research, in particular to a key parameter of typhoon that can be accurately extracted, which lays a reliable foundation for further research on typhoon intensity, typhoon path and impact of typhoon on people's lives. Extraction Method.

背景技术Background technique

台风的关键参数直接决定着台风强度、台风路径、台风对人民生活的影响等。台风的关键参数包括中心气压差、移动速度、移动方向等。通常根据参数的频率分布图,用某种确定的概率分布函数(如正态分布、对数正态分布、威布尔分布等)来进行模拟,然后利用χ2拟合检验法与K-S拟合检验法进行拟合检验,最终得到模拟点的台风关键参数概率模型,并根据台风关键参数概率模型对台风关键参数进行提取。但是,上述方法并不能完全模拟出关键参数的概率变化,因此,提取的关健参数准确性相对较差。The key parameters of a typhoon directly determine the intensity of the typhoon, the path of the typhoon, and the impact of the typhoon on people's lives. The key parameters of a typhoon include central air pressure difference, moving speed, moving direction, etc. Usually, according to the frequency distribution diagram of the parameters, a certain probability distribution function (such as normal distribution, lognormal distribution, Weibull distribution, etc.) is used to simulate, and then the χ2 fitting test method and KS fitting test are used Finally, the probability model of key typhoon parameters at the simulation point is obtained, and the key parameters of typhoon are extracted according to the probability model of key typhoon parameters. However, the above methods cannot completely simulate the probability changes of key parameters, so the accuracy of the extracted key parameters is relatively poor.

中国专利授权公开号:CN103177301A,授权公开日2013年6月26日,公开了一种台风灾害风险预估方法,针对指定监测区域的台风灾害所造成的损失数据进行统计分析,选择致灾因子危险性、孕灾环境敏感性、承灾体易损性和防灾减灾能力作为台风灾害风险评估指标体系,用模糊变换理论建立台风灾害风险预估模型,把台风预报结果作为预估模型的启动条件和输入条件,经过预估模型的计算和分析,得到未来一段时间被预估地区是否致灾以及致灾的灾害风险等级,从而提高气象灾害的预警能力。该发明的不足之处是,功能单一,无法提取台风关键参数。Chinese Patent Authorization Publication No.: CN103177301A, authorized publication date June 26, 2013, discloses a typhoon disaster risk estimation method, which conducts statistical analysis on the loss data caused by typhoon disasters in designated monitoring areas, and selects the risk of disaster-causing factors The typhoon disaster risk assessment index system is based on the sensitivity of disaster-pregnant environment, the vulnerability of disaster-bearing bodies and the ability of disaster prevention and mitigation. The fuzzy transformation theory is used to establish a typhoon disaster risk prediction model, and the typhoon forecast results are used as the starting conditions of the prediction model. And input conditions, through the calculation and analysis of the prediction model, whether the predicted area will cause disasters in the future and the disaster risk level of disasters will be obtained, so as to improve the early warning ability of meteorological disasters. The weak point of this invention is, single function, can't extract typhoon key parameter.

发明内容Contents of the invention

本发明的发明目的是为了克服现有技术中的方法无法提取的关键参数准确性差的不足,提供了一种可以准确提取台风的关键参数,为进一步研究台风强度、台风路径和台风对人民生活影响奠定了可靠基础的基于半正态模型的台风关键参数提取方法。The purpose of the invention is to overcome the deficiency of poor accuracy of key parameters that cannot be extracted by the methods in the prior art, and provide a key parameter that can accurately extract typhoons, so as to further study the typhoon intensity, typhoon path and typhoon’s impact on people’s lives The key parameter extraction method of typhoon based on half normal model has laid a reliable foundation.

一种基于半正态模型的台风关键参数提取方法,包括如下步骤:A method for extracting key parameters of a typhoon based on a half-normal model, comprising the following steps:

(1-1)将台风历史数据输入到计算机中,计算机选取与台风高发区域相对应的地图数据中的任一点为模拟点,以模拟点为圆心,R为半径做模拟圆,将经过模拟圆的台风的中心气压差Δp、台风移动方向、台风移动速度VT设为所述模拟点的台风关健参数;(1-1) Input the typhoon historical data into the computer, the computer selects any point in the map data corresponding to the typhoon high-incidence area as a simulation point, takes the simulation point as the center, and R as the radius to make a simulation circle, and will pass through the simulation circle The central air pressure difference Δp of the typhoon, the typhoon moving direction, and the typhoon moving speed V T are set as the key parameters of the typhoon at the simulation point;

台风中心气压差Δp是指台风中心气压与台风外围气压(一般取1010hPa)之差,通常用对数正态分布及Weibull分布进行来描述;The typhoon center pressure difference Δp refers to the difference between the typhoon center pressure and the typhoon peripheral pressure (generally 1010hPa), which is usually described by lognormal distribution and Weibull distribution;

台风移动方向由台风中心经纬度的位置计算得到,通常用正态分布或双正态分布描述;The moving direction of the typhoon is calculated from the latitude and longitude of the center of the typhoon, which is usually described by a normal distribution or a double normal distribution;

台风移动速度VT由台风中心前后记录点的经纬度求得,通常用正态分布或对数正态分布描述。The typhoon moving velocity V T is obtained from the latitude and longitude of the record points before and after the typhoon center, which is usually described by normal distribution or lognormal distribution.

(1-2)样本划分:(1-2) Sample division:

计算机用公式计算台风中心气压差Δp、台风移动方向和台风移动速度VT的频率;computer formula Calculate the frequency of typhoon center pressure difference Δ p , typhoon moving direction and typhoon moving speed V T ;

其中,x为台风中心气压差、台风移动方向或台风移动速度VT的任一个样本值,N为台风中心气压差、台风移动方向或台风移动速度VT的样本总数量,f为台风中心气压差、台风移动方向或台风移动速度VT的频率;Among them, x is any sample value of typhoon center pressure difference, typhoon moving direction or typhoon moving speed V T , N is the total number of samples of typhoon center pressure difference, typhoon moving direction or typhoon moving speed V T , f is typhoon center pressure Difference, frequency of typhoon moving direction or typhoon moving speed V T ;

计算机画出台风中心气压差Δp、台风移动方向和台风移动速度VT的概率分布直方图;The computer draws the probability distribution histogram of typhoon center pressure difference Δp, typhoon moving direction and typhoon moving speed V T ;

(1-2-1)选择各个概率分布直方图的分界点:(1-2-1) Select the cut-off point of each probability distribution histogram:

台风移动速度的直方图选择双正态分布的第二个分布的均值θ作为半正态模型的第一个分界点;以分界点θ为界,将台风移动方向划分为左半区与右半区;将左边区直方图的峰值Vpmax作为左边概率分布的分界点,将右边区直方图的峰值Vpmax作为右边概率分布的分界点;For the histogram of the typhoon’s moving speed, the mean θ of the second distribution of the binormal distribution is selected as the first dividing point of the half-normal model; with the dividing point θ as the boundary, the moving direction of the typhoon is divided into the left half and the right half area; the peak value V pmax of the histogram of the left area is used as the boundary point of the left probability distribution, and the peak value V pmax of the histogram of the right area is used as the boundary point of the right probability distribution;

计算机选择台风中心气压差Δp、台风移动方向直方图的峰值Vpmax作为概率分布的分界点;The computer selects the pressure difference Δp at the center of the typhoon and the peak value V pmax of the typhoon movement direction histogram as the dividing point of the probability distribution;

(1-2-2)对各个概率分布直方图均进行下述处理:(1-2-2) Perform the following processing on each probability distribution histogram:

以分界点为界,将台风中心气压差Δp、台风移动方向的直方图划分为左、右两个半区,得到左半区样本集合Al={xl|x≤Vpmax,x∈A}和右半区样本集合Ar={xr|x>Vpmax,x∈A};将台风移动速度直方图划分为左一区、左二区、右一区、右二区四个半区;Taking the boundary point as the boundary, divide the histogram of typhoon center pressure difference Δp and typhoon moving direction into left and right half areas, and obtain the left half area sample set A l ={x l |x≤V pmax , x∈A } and the sample set A r in the right half area = {x r |x>V pmax , x∈A}; divide the typhoon movement speed histogram into four halves: the first left area, the second left area, the first right area, and the second right area district;

其中,A为台风中心气压差样本集、台风移动方向样本集或台风移动速度VT样本集,xl为Al的任一个样本值,xr为Ar的任一个样本值;Among them, A is the sample set of pressure difference at the typhoon center, the sample set of typhoon moving direction or the sample set of typhoon moving velocity V T , x l is any sample value of A l , x r is any sample value of A r ;

(1-3)分别构造每个半区的对称映像样本集,构造得到半正态模型:(1-3) Construct the symmetrical image sample set of each half area separately, and construct the half normal model:

设定Al′={xl′|xl′=2Vpmax-xl,xl∈Al}为Al={xl|x≤Vpmax,x∈A}的对称映像样本集,Al与Al′组成第一半正态模型;Set A l ′={x l ′|x l ′=2V pmax -x l , x l ∈ A l } as the symmetrical image sample set of A l ={x l |x≤V pmax , x∈A}, A l and A l ′ form the first half-normal model;

设定Ar′={xr′|xr′=2Vpmax-xr,xr∈Ar}为Ar={xr|x>Vpmax,x∈A}的对称映像样本集,Ar′和Ar组成第二半正态模型;Set A r ′={x r ′|x r ′=2V pmax -x r , x r ∈ A r } as the symmetrical image sample set of A r ={x r |x>V pmax , x∈A}, A r ′ and A r form the second half normal model;

设定AL为左半区对称化样本概率密度集,AR右半区对称化样本概率密度集;其中,AR=Ar∪Ar′,AL=Al∪Al′;∪为并集运算符;Set AL as the left-half symmetric sample probability density set, and A R the right-half symmetric sample probability density set; where, A R =A r ∪A r ′, AL =A l ∪A l ′;∪ is the union operator;

(1-4)计算半正态模型特征参数:(1-4) Calculate the characteristic parameters of the half-normal model:

利用参数估计方法计算获得AL的风速数据的期望值μL与标准差δL,AR的风速数据的期望值μR与标准差δRUsing the parameter estimation method to calculate the expected value μ L and standard deviation δ L of the wind speed data of AL, and the expected value μ R and standard deviation δ R of the wind speed data of AR;

得到AL的台风关键参数的半正态模型特征参数AR的台风关键参数的半正态模型特征参数 Obtain the characteristic parameters of the half-normal model of the typhoon key parameters of AL The Characteristic Parameters of the Half-Normal Model of Typhoon Key Parameters in AR

其中,p为比例系数,p=ml/(ml+mr),ml为左区域台风关键参数样本数量,mr为右区域台风关键参数样本数量,N表示正态分布,(p,1]和[0,p]为比例系数的区间。Among them, p is the proportionality coefficient, p=m l /(m l +m r ), m l is the number of samples of key parameters of typhoon in the left region, m r is the number of samples of key parameters of typhoon in the right region, N represents the normal distribution, (p , 1] and [0, p] are the intervals of proportional coefficients.

本发明首先将台风历史数据输入到计算机中,计算机选取与台风高发区域相对应的地图数据中的任一点为模拟点,做出模拟圆,从而得到经过模拟圆的台风关键参数;The present invention firstly inputs typhoon historical data into the computer, and the computer selects any point in the map data corresponding to the typhoon high-incidence area as a simulation point, and makes a simulation circle, thereby obtaining key parameters of the typhoon passing through the simulation circle;

计算机画出台风中心气压差Δp、台风移动方向和台风移动速度VT的概率分布直方图;以分界点为界,将台风中心气压差Δp、台风移动方向的直方图划分为左、右两个半区,得到左半区样本集合Al={xl|x≤Vpmax,x∈A}和右半区样本集合Ar={xr|x>Vpmax,x∈A};将台风移动速度直方图划分为左一区、左二区、右一区、右二区四个半区;分别构造每个半区的对称映像样本集,构造得到半正态模型;计算半正态模型特征参数,得到AL的台风关键参数的半正态模型特征参数AR的台风关键参数的半正态模型特征参数 The computer draws a histogram of the probability distribution of typhoon center pressure difference Δp, typhoon moving direction and typhoon moving speed V T ; with the dividing point as the boundary, the histogram of typhoon center pressure difference Δp and typhoon moving direction is divided into left and right half area, the left half area sample set A l = {x l |x≤V pmax , x∈A} and the right half area sample set A r ={x r |x>V pmax , x∈A}; the typhoon The moving speed histogram is divided into four half-areas: the first area on the left, the second area on the left, the first area on the right, and the second area on the right; respectively construct a symmetrical image sample set for each half area, and construct a half-normal model; calculate the half-normal model Characteristic parameter, get the half-normal model characteristic parameter of the typhoon key parameter of AL The Characteristic Parameters of the Half-Normal Model of Typhoon Key Parameters in AR

本发明的半正态模型能够准确反映出关键参数之间的模糊性、关联性,从而可以更准确地实现对关键参数规则的提取。The half-normal model of the invention can accurately reflect the ambiguity and correlation between key parameters, so that the key parameter rules can be extracted more accurately.

本发明的台风历史数据来自《CMA-STI西北太平洋热带气旋最佳路径数据集》。The typhoon historical data of the present invention comes from "CMA-STI Northwest Pacific Tropical Cyclone Best Track Dataset".

作为优选,R为200至500KM。Preferably, R is 200 to 500KM.

作为优选,台风移动方向由台风中心经纬度的位置计算得到,在计算机中设定4个标准方向,4个标准方向分别为北行0°、东行90°、南行180°、西行-90°。As a preference, the moving direction of the typhoon is calculated from the latitude and longitude of the center of the typhoon, and four standard directions are set in the computer. The four standard directions are 0° to the north, 90° to the east, 180° to the south, and -90° to the west .

作为优选,参数估计方法包括最小二乘法、极大似然法、极大验后法、最小风险法和极小化极大熵法。Preferably, the parameter estimation methods include least square method, maximum likelihood method, maximum posterior method, minimum risk method and minimization maximum entropy method.

作为优选,Δp为0至135hPa;台风移动速度VT为2km/h至65km/h。Preferably, Δp is 0 to 135 hPa; typhoon moving speed V T is 2 km/h to 65 km/h.

为了实现上述目的,本发明采用以下技术方案:In order to achieve the above object, the present invention adopts the following technical solutions:

因此,本发明具有如下有益效果:(1)可以准确提取台风的关键参数,为进一步研究台风强度、台风路径和台风对人民生活影响奠定了可靠基础;(2)为台风风能的利用提供可靠数据基础。Therefore, the present invention has the following beneficial effects: (1) can accurately extract the key parameters of the typhoon, laying a reliable foundation for further research on typhoon intensity, typhoon path and typhoon's impact on people's lives; (2) provide reliable data for the utilization of typhoon wind energy Base.

附图说明Description of drawings

图1是本发明的模拟圆的一种结构示意图;Fig. 1 is a kind of structural representation of simulation circle of the present invention;

图2是本发明的实施例的一种流程图;Fig. 2 is a kind of flowchart of the embodiment of the present invention;

图3是常规方法的台风移动方向拟合图;Figure 3 is a fitting diagram of typhoon movement direction by conventional method;

图4是本发明的A区台风移动方向拟合图;Fig. 4 is the fitting figure of typhoon moving direction in A district of the present invention;

图5是本发明的B区台风移动方向拟合图;Fig. 5 is the fitting diagram of typhoon movement direction in B area of the present invention;

图6是本发明的C区台风移动方向拟合图;Fig. 6 is a fitting diagram of typhoon moving direction in C zone of the present invention;

图7是本发明的D区台风移动方向拟合图;Fig. 7 is the fitting figure of typhoon movement direction in D area of the present invention;

图8是常规方法的台风移动速度拟合图;Figure 8 is a fitting diagram of typhoon moving speed by conventional method;

图9是本发明的左区域台风移动速度拟合图;Fig. 9 is a fitting diagram of typhoon moving speed in the left region of the present invention;

图10是本发明的右区域台风移动速度拟合图;Fig. 10 is a fitting diagram of typhoon moving speed in the right region of the present invention;

图11是常规方法的台风中心气压差拟合图;Figure 11 is a fitting diagram of the typhoon center air pressure difference by the conventional method;

图12是本发明的左区域台风中心气压差拟合图;Fig. 12 is a fitting diagram of the left regional typhoon center air pressure difference of the present invention;

图13是本发明的右区域台风中心气压差拟合图。Fig. 13 is a fitting diagram of the air pressure difference at the center of the typhoon in the right region of the present invention.

图中:模拟点1、模拟圆2、海岸线3、台风路径4。In the figure: simulation point 1, simulation circle 2, coastline 3, typhoon track 4.

具体实施方式Detailed ways

下面结合附图和具体实施方式对本发明做进一步的描述。The present invention will be further described below in conjunction with the accompanying drawings and specific embodiments.

如图2所示的实施例是一种基于半正态模型的台风关键参数提取方法,包括如下步骤:The embodiment shown in Figure 2 is a method for extracting key parameters of a typhoon based on a half-normal model, comprising the following steps:

步骤100,将台风历史数据输入到计算机中,计算机选取与台风高发区域相对应的地图数据中的任一点为模拟点1,以模拟点为圆心,R为半径做模拟圆2,将经过模拟圆的台风的中心气压差Δp、台风移动方向、台风移动速度VT设为所述模拟点的台风关键参数;Step 100, input the typhoon history data into the computer, the computer selects any point in the map data corresponding to the typhoon high-incidence area as the simulation point 1, takes the simulation point as the center, and R as the radius to make the simulation circle 2, and will pass through the simulation circle The typhoon's central air pressure difference Δp, typhoon moving direction, and typhoon moving speed V T are set as the typhoon key parameters of the simulation point;

如图1所示,本实施例中,计算机选取与海岸线3相对应的地图数据中的任一点为模拟点,以模拟点为圆心,R=250KM为半径做模拟圆,将台风路径4经过模拟圆的各个台风的中心气压差Δp、台风移动方向、台风移动速度VT设为模拟点的台风关键参数;As shown in Figure 1, in the present embodiment, the computer selects any point in the map data corresponding to the coastline 3 as the simulation point, takes the simulation point as the center, and R=250KM as the radius to do the simulation circle, and the typhoon path 4 is simulated The central air pressure difference Δp of each typhoon in the circle, the typhoon moving direction, and the typhoon moving speed V T are set as the key parameters of the typhoon at the simulation point;

如图2所示,步骤200,样本划分:As shown in Figure 2, step 200, sample division:

计算机用公式计算台风中心气压差Δp、台风移动方向和台风移动速度VT的频率;computer formula Calculate the frequency of typhoon center pressure difference Δp, typhoon moving direction and typhoon moving speed V T ;

其中,x为台风中心气压差、台风移动方向或台风移动速度VT的任一个样本值,N为台风中心气压差、台风移动方向或台风移动速度VT的样本总数量,f为台风中心气压差、台风移动方向或台风移动速度VT的频率;Among them, x is any sample value of typhoon center pressure difference, typhoon moving direction or typhoon moving speed V T , N is the total number of samples of typhoon center pressure difference, typhoon moving direction or typhoon moving speed V T , f is typhoon center pressure Difference, frequency of typhoon moving direction or typhoon moving speed V T ;

计算机画出台风中心气压差Δp、台风移动方向和台风移动速度VT的概率分布直方图;The computer draws the probability distribution histogram of typhoon center pressure difference Δp, typhoon moving direction and typhoon moving speed V T ;

计算机选择台风移动速度的直方图选择双正态分布的第二个分布的均值θ作为半正态模型的第一个分界点;以分界点θ为界,将台风移动方向划分为左半区与右半区;将左边区直方图的峰值Vpmax作为左边概率分布的分界点,将右边区直方图的峰值Vpmax作为右边概率分布的分界点;The computer selects the histogram of the typhoon’s moving speed and selects the mean θ of the second distribution of the double normal distribution as the first cut-off point of the half-normal model; with the cut-off point θ as the boundary, the direction of the typhoon’s movement is divided into the left half and Right half area; the peak value V pmax of the histogram of the left area is used as the boundary point of the left probability distribution, and the peak value V pmax of the histogram of the right area is used as the boundary point of the right probability distribution;

计算机选择台风中心气压差Δp、台风移动方向直方图的峰值Vpmax作为概率分布的分界点;The computer selects the pressure difference Δp at the center of the typhoon and the peak value V pmax of the typhoon movement direction histogram as the dividing point of the probability distribution;

计算机根据不同的分界点,将每个直万图划分不同的区域,例如,将台风中心气压差Δp、台风移动方向划分为左、右两个半区,得到左半区样本集合Al={xl|x≤Vpmax,x∈A}和右半区样本集合Ar={xr|x>Vpmax,x∈A};将台风移动速度直方图划分为四个半区;The computer divides each graph into different regions according to different demarcation points. For example, divide the typhoon center pressure difference Δp and typhoon moving direction into left and right half regions, and obtain the left half region sample set A l ={ x l |x≤V pmax , x∈A} and the right half-area sample set A r ={x r |x>V pmax , x∈A}; divide the typhoon velocity histogram into four half-areas;

其中,A为台风中心气压差样本集、台风移动方向样本集或台风移动速度VT样本集,xl为Al的任一个样本值,xr为Ar的任一个样本值;Among them, A is the sample set of pressure difference at the typhoon center, the sample set of typhoon moving direction or the sample set of typhoon moving velocity V T , x l is any sample value of A l , x r is any sample value of A r ;

步骤300,分别构造每个半区的对称映像样本集,构造得到半正态模型:Step 300, respectively constructing a symmetrical image sample set for each half region, and constructing a half normal model:

设定Al′={xl′|xl′=2Vpmax-xl,xl∈Al}为Al={xl|x≤Vpmax,x∈A}的对称映像样本集,Al与Al′组成第一半正态模型;Set A l ′={x l ′|x l ′=2V pmax -x l , x l ∈ A l } as the symmetrical image sample set of A l ={x l |x≤V pmax , x∈A}, A l and A l ′ form the first half-normal model;

设定Ar′={xr′|xr′=2Vpmax-xr,xr∈Ar}为Ar={xr|x>Vpmax,x∈A}的对称映像样本集,Ar′和Ar组成第二半正态模型;Set A r ′={x r ′|x r ′=2V pmax -x r , x r ∈ A r } as the symmetrical image sample set of A r ={x r |x>V pmax , x∈A}, A r ′ and A r form the second half normal model;

设定AL为左半区对称化样本概率密度集,AR右半区对称化样本概率密度集;其中,AR=Ar∪Ar′,AL=Al∪Al′;Set AL as the probability density set of symmetric samples in the left half area, and the symmetric sample probability density set in the right half area of A R ; among them, A R =A r ∪A r ′, AL =A l ∪A l ′;

步骤400,计算半正态模型特征参数:Step 400, calculating the characteristic parameters of the half-normal model:

利用参数估计方法计算获得AL的风速数据的期望值μL与标准差δL,AR的风速数据的期望值μR与标准差δRUsing the parameter estimation method to calculate the expected value μ L and standard deviation δ L of the wind speed data of AL, and the expected value μ R and standard deviation δ R of the wind speed data of AR;

得到AL的台风关键参数的半正态模型特征参数AR的台风关键参数的半正态模型特征参数 Obtain the characteristic parameters of the half-normal model of the typhoon key parameters of AL The Characteristic Parameters of the Half-Normal Model of Typhoon Key Parameters in AR

其中,p为比例系数,p=ml/(ml+mr),ml为左区域台风关健参数样本数量,mr为右区域台风关健参数样本数量,N表示正态分布,(p,1]和[0,p]为比例系数的区间。Among them, p is the proportional coefficient, p=m l /(m l + m r ), m l is the number of samples of key typhoon parameters in the left region, m r is the number of samples of key parameters of typhoon in the right region, and N represents the normal distribution, (p, 1] and [0, p] are intervals of proportional coefficients.

仿真试验:Simulation test:

一、生成台风的中心气压差Δp样本集、台风移动方向样本集和台风移动速度VT样本集中的混合半云风速样本集:1. Generate the typhoon’s central air pressure difference Δp sample set, typhoon moving direction sample set, and typhoon moving speed V T sample set and the mixed semi-cloud wind speed sample set:

(1)利用MATLAB自带randn(1)函数随机产生随机数a;(1) Use the randn(1) function that comes with MATLAB to randomly generate a random number a;

(2)判断a值的大小,如果0≤a≤p,利用MATLAB自带normrnd指令生成符合分布且位于左区域数值范围的随机数;如果p<a≤1,利用MATLAB自带normrnd指令生成符合分布且位于右区域数值范围的的随机数;(2) Determine the value of a, if 0≤a≤p, use the normrnd command that comes with MATLAB to generate Random numbers that are distributed and located in the value range of the left area; if p<a≤1, use the normrnd command that comes with MATLAB to generate A random number distributed and located in the value range of the right area;

(3)重复步骤(1)(2),获得N个数据,将N个数据按照由小至大排序并组成混合半云风速样本集;(3) Repeat steps (1) and (2) to obtain N data, sort the N data from small to large and form a mixed semi-cloud wind speed sample set;

二、对生成台风的中心气压差Δp样本集、台风移动方向样本集和台风移动速度VT样本集的半正态模型进行校验:2. Verify the half-normal model of the typhoon’s center pressure difference Δp sample set, typhoon moving direction sample set, and typhoon moving speed V T sample set:

(1)计算残差值(1) Calculate the residual value

将台风的中心气压差Δp样本集、台风移动方向样本集或台风移动速度VI样本集中的各个样本值xi由小至大排序,利用公式对生成的混合半云风速样本集进行残差值求解;Sort the sample values x i of the typhoon central pressure difference Δp sample set, typhoon moving direction sample set or typhoon moving speed V I sample set from small to large, using the formula Solve the residual value of the generated mixed semi-cloud wind speed sample set;

其中,xi′为与xi相对应的混合半云风速样本集的样本值。Among them, xi is the sample value of the mixed semi-cloud wind speed sample set corresponding to xi .

(2)拟合优度检验(2) Goodness of fit test

拟合优度检验是用于检验某种理论分布与原数据分布是否相一致的统计方法。The goodness of fit test is a statistical method used to test whether a certain theoretical distribution is consistent with the original data distribution.

RR 22 == 11 -- ee ww &Sigma;&Sigma; ii == 11 Mm (( xx ii -- xx &OverBar;&OverBar; )) 22 &times;&times; 100100 %%

式中,ew为残差值;为样本集中各个样本值xi的均值;R2值越大,模型拟合效果越佳。In the formula, e w is the residual value; is the mean value of each sample value x i in the sample set; the larger the value of R 2 , the better the model fitting effect.

三、仿真结果对比3. Comparison of simulation results

本实施例中,分别利用常规方法、半正态提取出的某模拟点(E111.83°,N21.58°)台风关键参数对比。In this embodiment, the key parameters of a typhoon at a simulated point (E111.83°, N21.58°) extracted by conventional methods and half-normal were compared.

1.台风移动方向1. Typhoon movement direction

使用正态分布、双正态分布对台风移动方向概率分布分别进行拟合,如图3所示。The normal distribution and double normal distribution are used to fit the probability distribution of typhoon movement direction respectively, as shown in Figure 3.

由如图3拟合效果,可直观看出,双正态分布模型拟合效果优于正态分布。由台风移动方向的概率直方图知,该分布具有双峰特性,可于两波峰与两波峰间的波谷处进行区域划分。将频率直方图划分为A、B、C及D区域,并分别构建四区域的对称数据,使用正态模型对其进行拟合。拟合效果如图4、图5、图6、图7所示。From the fitting effect shown in Figure 3, it can be seen intuitively that the fitting effect of the binormal distribution model is better than that of the normal distribution. According to the probability histogram of typhoon moving direction, the distribution has bimodal characteristics, and regions can be divided at the two peaks and the trough between the two peaks. The frequency histogram is divided into A, B, C and D areas, and the symmetrical data of the four areas are respectively constructed, and the normal model is used to fit them. The fitting effect is shown in Figure 4, Figure 5, Figure 6, and Figure 7.

由四个正态分布的数据参数得到相应的半正态模型数字特征,使用正态模型、双正态模型及四个半正态模型随机生成风速数据,拟合度计算结果如表1所示。The digital characteristics of the corresponding half-normal model were obtained from the four normally distributed data parameters, and the wind speed data were randomly generated using the normal model, the double normal model and the four half-normal models. The results of the fitting degree calculation are shown in Table 1 .

表1不同模型台风移动方向拟合度对比Table 1 Comparison of fitting degree of typhoon movement direction of different models

模型Model 正态分布normal distribution 双正态分布binormal distribution 半正态half normal 拟合度%Fit % 92.480 792.480 7 99.341 099.341 0 99.508 399.508 3

由表1结果知,半正态模型拟合度最高,较正态分布、双正态分布分别高出7.027 6%和0.167 3%,表明半正态模型对台风移动方向具有很好的拟合效果。From the results in Table 1, the half-normal model has the highest fitting degree, which is 7.027 6% and 0.167 3% higher than the normal distribution and double normal distribution, respectively, indicating that the half-normal model has a good fit for the direction of typhoon movement Effect.

2.台风速度2. Typhoon speed

根据统计结果,台风速度概率分布具有单峰特性,使用正态分布、对数正态分布模型分别拟合,拟合效果如图8所示。According to the statistical results, the probability distribution of typhoon speed has unimodal characteristics, and the normal distribution and lognormal distribution models are used to fit them respectively. The fitting effect is shown in Figure 8.

由于台风速度概率分布具有单峰特性,将概率分布划分为左、右两区域,并构建左右区域的映像样本,使用正态模型分别拟合。拟合效果如图9、图10所示。Since the typhoon speed probability distribution has a unimodal characteristic, the probability distribution is divided into left and right regions, and image samples of the left and right regions are constructed, and the normal model is used to fit them respectively. The fitting effect is shown in Figure 9 and Figure 10.

由两个正态分布的数据参数得到相应的半正态模型数字特征,使用正态模型、对数正态模型及两个半正态模型随机生成风速数据,拟合度计算结果如表2所示。The digital characteristics of the corresponding half-normal model were obtained from the data parameters of two normal distributions, and the wind speed data were randomly generated using the normal model, the lognormal model and two half-normal models, and the fitting degree calculation results are shown in Table 2 Show.

表2不同模型台风移动速度拟合度对比Table 2 Comparison of typhoon movement speed fitting degree of different models

模型Model 正态分布normal distribution 对数正态分布lognormal distribution 半正态half normal 拟合度%Fit % 92.117592.1175 96.159396.1593 97.250197.2501

由表2结果知,半正态模型拟合度最高,较正态分布、对数正态分布分别高出5.132 6%和1.090 8%,表明半正态模型对台风移动速度具有很好的拟合效果。From the results in Table 2, the half-normal model has the highest fitting degree, which is 5.132 6% and 1.090 8% higher than the normal distribution and the lognormal distribution, respectively, indicating that the half-normal model has a good fit for typhoon moving speed. combined effect.

3.台风中心气压差3. Air pressure difference at the typhoon center

根据统计结果,台风中心气压差概率分布具有单峰特性,使用对数正态分布、Weibull分布模型分别拟合,拟合效果如图11所示。According to the statistical results, the probability distribution of pressure difference in the center of the typhoon has a unimodal characteristic, and the lognormal distribution and Weibull distribution models are used to fit them respectively. The fitting results are shown in Figure 11.

由于台风中心气压差概率分布具有单峰特性,将概率分布划分为左、右两区域,并构建左右区域的映像样本,使用正态模型分别拟合。拟合效果如图12、图13所示。Since the probability distribution of pressure difference in the typhoon center has a unimodal characteristic, the probability distribution is divided into left and right regions, and the image samples of the left and right regions are constructed, and the normal model is used to fit them respectively. The fitting effect is shown in Figure 12 and Figure 13.

由两个正态分布的数据参数得到相应的半正态模型数字特征,使用对数正态模型、Weibull模型及两个半正态模型随机生成风速数据,拟合度计算结果如表3所示。The digital characteristics of the corresponding half-normal model were obtained from two normally distributed data parameters, and the lognormal model, Weibull model, and two half-normal models were used to randomly generate wind speed data. The fitting results are shown in Table 3 .

表3不同模型台风中心气压差拟合度对比Table 3 Comparison of fitting degree of typhoon center pressure difference of different models

由表3结果知,半正态模型拟合度最高,较对数正态分布、Weibull分布分别高出12.455 1%和0.029%,表明半正态模型对台风中心气压差具有很好的拟合效果。From the results in Table 3, the half-normal model has the highest fitting degree, which is 12.455 1% and 0.029% higher than the lognormal distribution and Weibull distribution, respectively, indicating that the half-normal model has a good fit for the typhoon center pressure difference Effect.

由以上模型的计算结果对比可知,将数据概率分布进行区域划分后,使用正态分布分别拟合,得到半正态分布的参数后,对原数据进行数据拟合具有高拟合度、通用性强的特点,因此提取的台风关键参数准确性更高,为进一步研究台风强度、台风路径和台风对人民生活影响奠定了可靠基础;为台风风能的利用提供可靠数据基础。From the comparison of the calculation results of the above models, it can be seen that after the data probability distribution is divided into regions, the normal distribution is used to fit them separately, and after the parameters of the half-normal distribution are obtained, the data fitting of the original data has a high degree of fitting and versatility Therefore, the accuracy of the key parameters of the typhoon extracted is higher, laying a reliable foundation for further research on typhoon intensity, typhoon track and typhoon’s impact on people’s lives; providing a reliable data basis for the utilization of typhoon wind energy.

应理解,本实施例仅用于说明本发明而不用于限制本发明的范围。此外应理解,在阅读了本发明讲授的内容之后,本领域技术人员可以对本发明作各种改动或修改,这些等价形式同样落于本申请所附权利要求书所限定的范围。It should be understood that this embodiment is only used to illustrate the present invention but not to limit the scope of the present invention. In addition, it should be understood that after reading the teachings of the present invention, those skilled in the art can make various changes or modifications to the present invention, and these equivalent forms also fall within the scope defined by the appended claims of the present application.

Claims (5)

1. the typhoon key parameter extracting method based on half normal model, is characterized in that, comprises the steps:
(1-1) typhoon historical data is input in computing machine, any point in the computer selecting map datum corresponding with typhoon region occurred frequently is simulation points (1), taking simulation points as the center of circle, R is that radius does simulation circle (2), by central gas pressure reduction Δ p, Typhoon Tracks direction, the Typhoon Tracks speed V of the typhoon through simulation circle tbe made as the typhoon crux parameter of described simulation points;
(1-2) sample is divided:
Computing machine formula calculate the poor Δ p of typhoon central pressure, Typhoon Tracks direction and Typhoon Tracks speed V tfrequency;
Wherein, x is center of typhoon draught head, Typhoon Tracks direction or Typhoon Tracks speed V tany sample value, N is center of typhoon draught head, Typhoon Tracks direction or Typhoon Tracks speed V ttotal sample number amount, f is center of typhoon draught head, Typhoon Tracks direction or Typhoon Tracks speed V tfrequency;
Computing machine draws center of typhoon draught head Δ p, Typhoon Tracks direction and Typhoon Tracks speed V tprobability distribution histogram;
(1-2-1) select the histogrammic separation of each probability distribution:
The average θ of second distribution of the two normal distributions of histogram selection of Typhoon Tracks speed is as first separation of half normal model; Taking separation θ as boundary, Typhoon Tracks direction is divided into left half-court and right half-court; By histogrammic left side district peak value v pmaxas the separation of left side probability distribution, by histogrammic the right district peak value v pmaxas the separation of the right probability distribution;
Computing machine is selected the peak value v of center of typhoon draught head Δ p, Typhoon Tracks direction histogram pmaxas the separation of probability distribution;
(1-2-2) each probability distribution histogram is all carried out to following processing:
Taking separation as boundary, the histogram of center of typhoon draught head Δ p, Typhoon Tracks direction is divided into left and right Liang Geban district, obtain left half-court sample set A l={ x l| x≤V pmax, x ∈ A} and right half-court sample set A r={ x r| x > V pmax, x ∈ A}; Typhoon Tracks velocity histogram is divided into the first from left district, the second from left district, You Yiqu, Si Geban district of You Er district;
Wherein, A is center of typhoon draught head sample set, Typhoon Tracks direction sample set or Typhoon Tracks speed V tsample set, X lfor A lany sample value, x rfor A rany sample value;
(1-3) construct respectively the symmetry reflection sample set in each halfth district, structure obtains half normal model:
Set A l'={ x l' | x l'=2V pmax-x l, x l∈ A lbe A l={ x l| x≤V pmax, the symmetry reflection sample set of x ∈ A}, A lwith A l' composition the first half normal models;
Set A r'={ x r' | x r'=2V pmax-x r, x r∈ A rbe A r={ x r| x > V pmax, the symmetry reflection sample set of x ∈ A}, A r' and A rform the second half normal models;
Set A lfor left half-court symmetrization sample probability density collection, A rright half-court symmetrization sample probability density collection; Wherein, A r=A r∪ A r', A l=A l∪ A l';
(1-4) calculate half normal model characteristic parameter:
Utilize method for parameter estimation to calculate and obtain A lthe expectation value μ of air speed data lwith standard deviation δ l, A rthe expectation value μ of air speed data rwith standard deviation δ r;
Obtain A lhalf normal model characteristic parameter of typhoon key parameter a rhalf normal model characteristic parameter of typhoon key parameter
Wherein, p is scale-up factor, p=m l/ (m l+ m r), m lfor left region typhoon key parameter sample size, m rfor right region typhoon key parameter sample size, N represents normal distribution, (p, 1] and [0, the p] interval that is scale-up factor.
2. the typhoon key parameter extracting method based on half normal model according to claim 1, is characterized in that, R is 200 to 500KM.
3. the typhoon key parameter extracting method based on half normal model according to claim 1, Typhoon Tracks direction is obtained by the position calculation of center of typhoon longitude and latitude, in computing machine, set 4 reference directions, 4 reference directions are respectively 0 ° of northern row, eastbound 90 °, 180 ° of southern row, head west-90 °.
4. the typhoon key parameter extracting method based on half normal model according to claim 1, is characterized in that, method for parameter estimation comprises least square method, maximum-likelihood method, Maximum Verified Method, minimum risk method and minimization Maximum entropy method.
5. according to the typhoon key parameter extracting method based on half normal model described in claim 1 or 2 or 3 or 4, it is characterized in that, Δ p is 0 to 135hPa; Typhoon Tracks speed V tfor 2km/h to 65km/h.
CN201410016770.5A 2014-01-14 2014-01-14 Typhoon key parameter extraction method based on half-normal model Pending CN103995951A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410016770.5A CN103995951A (en) 2014-01-14 2014-01-14 Typhoon key parameter extraction method based on half-normal model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410016770.5A CN103995951A (en) 2014-01-14 2014-01-14 Typhoon key parameter extraction method based on half-normal model

Publications (1)

Publication Number Publication Date
CN103995951A true CN103995951A (en) 2014-08-20

Family

ID=51310114

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410016770.5A Pending CN103995951A (en) 2014-01-14 2014-01-14 Typhoon key parameter extraction method based on half-normal model

Country Status (1)

Country Link
CN (1) CN103995951A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107330583A (en) * 2017-06-09 2017-11-07 哈尔滨工业大学深圳研究生院 A kind of complete trails typhoon risk analysis method based on statistical dynamics
CN111921192A (en) * 2020-08-31 2020-11-13 网易(杭州)网络有限公司 Control method and device of virtual object

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
周旋等: "基于星载SAR数据的台风参数估计及风场构建", 《中国科学:地球科学》 *
唐建等: "基于台风风场模型的台风浪数值模拟", 《海洋湖沼通报》 *
李瑞龙: "基于改进的台风关键参数的台风极值风速预测", 《万方学位论文》 *
李茜等: "Shapiro台风风场模型及其数值模拟", 《自然灾害学报》 *
肖玉凤等: "基于数值模拟的台风危险性分析综述(Ⅰ)-台风风场模型", 《自然灾害学报》 *
陈朝晖等: "台风极值风速的数值模拟及分布模型", 《重庆大学学报》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107330583A (en) * 2017-06-09 2017-11-07 哈尔滨工业大学深圳研究生院 A kind of complete trails typhoon risk analysis method based on statistical dynamics
CN107330583B (en) * 2017-06-09 2020-06-19 哈尔滨工业大学深圳研究生院 An all-track typhoon risk analysis method based on statistical dynamics
CN111921192A (en) * 2020-08-31 2020-11-13 网易(杭州)网络有限公司 Control method and device of virtual object
CN111921192B (en) * 2020-08-31 2024-02-23 网易(杭州)网络有限公司 Virtual object control method and device

Similar Documents

Publication Publication Date Title
CN104573198B (en) Based on the theoretical digital cores of random fractal and pore network model reconstructing method
CN104064051B (en) A kind of passenger&#39;s portable mobile terminal and a locating information dynamic matching method that rides in a bus
CN102332210B (en) Method for extracting real-time urban road traffic flow data based on mobile phone positioning data
CN111079859B (en) A Passive Multi-Site Multi-Target DF Cross-Location and False Point Removal Method
CN102253426B (en) Method for predicting lightning motion by space density clustering
CN108830331A (en) A kind of Ground Penetrating Radar object detection method based on full convolutional network
CN106355166A (en) Monitoring video and remote sensing image-based dust-haze spreading path drawing and source determination method
CN113821986B (en) Vortex sea surface signal and underwater key parameter inversion method based on random forest
CN110503721B (en) Fractured Terrain Preservation Method Based on Weighted Radial Basis Function Interpolation
CN111985552B (en) Method for detecting diseases of thin strip-shaped structure of airport pavement under complex background
CN110852243B (en) Road intersection detection method and device based on improved YOLOv3
CN103530500B (en) Debris fans plane net is formatted method and dangerous measure
CN104504709A (en) Feature ball based classifying method of three-dimensional point-cloud data of outdoor scene
CN106412841A (en) Indoor location method based on DBSCAN algorithm
CN105678047A (en) Wind field characterization method with empirical mode decomposition noise reduction and complex network analysis combined
CN114966845A (en) Unsupervised clustering seismic waveform real-time classification detection method
CN101789048A (en) Method for quickly extracting critical area of layout
CN105938655A (en) Real-time traffic state evaluation method based on Gaussian mixture model
CN103995951A (en) Typhoon key parameter extraction method based on half-normal model
CN103093044A (en) Electric transmission line icing galloping distribution diagram surveying and mapping method
CN106872579B (en) Method for Classifying Rock Mass Quality by Fitting Rock Mass Wave Velocity with Normal Distribution
CN108009134A (en) A kind of Information of Ancient Human Activity method for digging
CN105005649A (en) Potential underground water distribution mapping method
CN109543236A (en) Method is determined based on the rock structural plane roughness statistical sample number of variation lines several levels score analysis
CN102928871A (en) Rotation diamond based attribute extraction and fault description method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20140820

RJ01 Rejection of invention patent application after publication