删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

基于Poisson过程的沿海地区飓风灾害评估

本站小编 Free考研考试/2020-04-15

李全旺 , 王草 , 邹阿鸣 , 庞龙 , 张龙
清华大学 土木工程系, 北京 100084

收稿日期: 2015-06-14
基金项目: 交通运输部重大专项基金项目(201332849A090)
作者简介: 李全旺(1975-),男,副教授。E-mail:li_quanwang@tsinghua.edu.cn


摘要:对沿海地区进行飓风灾害评估时,考虑到已有方法中往往忽略飓风发生的间歇性特性,导致分析结果不准确,该文提出了飓风灾害评估的新方法。将飓风发生的随机过程描述为平稳Poisson过程,并在此基础上给出了飓风灾害评估公式。选取美国Miami-Dade县1900-2010年的飓风数据进行分析,并应用该文提出的方法对该县未来50年内的飓风灾害进行评估。结果表明:Poisson过程能够准确刻画飓风发生的随机过程。如果忽略飓风的间歇性特性,会高估飓风灾害的均值和标准差,但低估其变异性。在50年的评估期内,如果忽略飓风的间歇性特性,则飓风灾害的均值和方差分别被高估41.4%和20.05%。
关键词: 飓风 灾害评估 Poisson过程 间歇性 概率模型
Hurricane damage assessments for coastal areas based on a Poisson model
LI Quanwang, WANG Cao, ZOU Aming, PANG Long, ZHANG Long
Department of Civil Engineering, Tsinghua University, Beijing 100084, China


Abstract:Previous hurricane damage assessments for coastal areas have ignored the intermittent nature of hurrucanes. This paper presents an hurricane damage assessment method that accounts for the intermittency. A Poisson model is used in this paper to describe the hurricane stochastic process with explicit formulas for the hurricane damage assessment. Hurricane data for Miami-Dade County, USA, from 1900 to 2010 were used to illustrate the method. The Poisson model provides a reasonable description of the hurricane stochastic characteristics. The mean and the variance of the hurricane damage are overestimated while the coefficient of variance is underestimated if the intermittency is ignored. For example, for a 50-year service period, the mean cumulative hurricane damage is overestimated by 41.4% while the variance is over estimated by 20.05% if the intermittency is ignored.
Key words: hurricanesdamage assessmentPoisson processintermittencyprobabilistic model
对于沿海地区(如中国的东南沿海、美国的东南海岸等),飓风灾害往往会导致巨大的经济损失和人员伤亡。2005年美国Katrina飓风带来的直接经济损失超过1亿美元[1]。据福建省民政厅的统计结果,2013年7月的台风“苏力”造成福建省9个地级市及平潭综合实验区80个县(市、区)103万人受灾,直接经济损失达17.7亿元[2]。对飓风灾害进行风险评估,近年来得到了众多****的关注[3-7]。 Bjarnadottir等[5]对美国佛罗里达州Miami-Dade县的飓风灾害进行了评估,该方法首先根据每年最大风速的概率模型,计算得到每年的飓风灾害,然后对其累加得到评估期内的累积灾害。 Li等[6]对澳大利亚North Queensland地区由飓风引起的经济损失进行了评估,该方法也是先求得每年的灾害损失,然后对其进行累加。
对某一特定地区,飓风的发生存在着间歇性(例如,每隔3—5年发生一次)[8-9]。这表明,根据每年最大风速的概率分布来计算每年的飓风灾害,得到的结果是不准确的。为了体现飓风过程的间歇性特性,本文提出采用Poisson过程模型来描述飓风的发生。在此基础上,提出了飓风灾害评估新的计算方法。
1 飓风灾害评估方法1.1 已有方法考虑时间(0,T]内的累积飓风灾害DT,已有研究[5-6]中采用下式计算:
${{D}_{T}}={{d}_{1}}+{{d}_{2}}+\ldots +{{d}_{T}}.$ (1)
其中,di表示第i年内的飓风灾害,且
$E{{d}_{i}}=\int {{F}_{D}}\left( v \right)\cdot {{f}_{v}}\left( v,\text{ }i \right)dv,$ (2)
$\begin{align} & Var{{d}_{i}}= \\ & \int {{\left[ {{F}_{D}}\left( v \right)-E{{d}_{i}} \right]}^{2}}\cdot {{f}_{v}}\left( v,i \right)dv. \\ \end{align}$ (3)
其中: E()和Var() 分别表示括号中变量的均值和方差; FD(v)为飓风灾害函数; fv(v,i)为第i年的最大风速的概率密度函数PDF。假定各di互相独立,则由式(1)—(3)得:
$\begin{align} & E\left( {{D}_{T}} \right)=\sum\limits_{i=1}^{T}{{}}E{{d}_{i}}= \\ & \sum\limits_{i=1}^{T}{{}}\int {{F}_{D}}\left( v \right)\cdot {{f}_{v}}\left( v,\text{ }i \right)dv. \\ \end{align}$ (4)
$\begin{align} & Var\left( {{D}_{T}} \right)=\sum\limits_{i=1}^{T}{{}}Var{{d}_{i}}= \\ & \sum\limits_{i=1}^{T}{{}}\int {{\left[ {{F}_{D}}\left( v \right)-E{{d}_{i}} \right]}^{2}}\cdot {{f}_{v}}\left( v,\text{ }i \right)dv. \\ \end{align}$ (5)
式(4)和(5)即为已有的评估某一地区未来一段时间内累积飓风灾害的方法。然而,对于沿海地区,飓风的发生是罕见的,可能存在着一年中不发生飓风的情形[8-9]。因此,式(4)和(5)中的fv(v,i)是“虚拟”的年最大风速。这表明,利用已有方法计算累积飓风灾害DT时,得到的结果是不准确的。
1.2 本文方法为了更加准确地评估时间(0,T]内的累积飓风灾害DT,本文提出将飓风发生的随机过程描述为平稳的Poisson过程,其Poisson强度记为λ(即单位时间内平均有λ次飓风发生)。将时间段(0,T]划分为N个等长的小区间(N足够大),则根据Poisson过程的性质,在每个区间上,发生一次飓风的概率为$\frac{\lambda T}{N}$, 不发生飓风的概率为1-$\frac{\lambda T}{N}$。 若每次飓风风速v的PDF为${\tilde{f}}$(v),将第j个区间上的飓风灾害记作${\tilde{d}}$jj=1,2,…,N,则
$E\left( {{{\tilde{d}}}_{j}} \right)=\frac{\lambda T}{N}\cdot \int {{F}_{D}}\left( v \right)\cdot \tilde{f}\left( v \right)dv,$ (6)
$\begin{align} & Var\left( {{{\tilde{d}}}_{j}} \right)=\frac{\lambda T}{N}\cdot \int {{\left[ {{F}_{D}}\left( v \right)-E\left( {{{\tilde{d}}}_{j}} \right) \right]}^{2}}\cdot \tilde{f}\left( v \right)dv,+ \\ & \left( 1-\frac{\lambda T}{N} \right)\cdot {{E}^{2}}\left( {{{\tilde{d}}}_{i}} \right).~ \\ \end{align}$ (7)
根据式(6)和(7),假设各d~j互相独立,有
$\begin{align} & E\left( {{D}_{T}} \right)=\sum\limits_{j=1}^{N}{{}}E\left( {{{\tilde{d}}}_{j}} \right)= \\ & \lambda T\cdot \int {{F}_{D}}\left( v \right)\cdot \tilde{f}\left( v \right)dv, \\ \end{align}$ (8)
$\begin{align} & Var{{D}_{T}}=\sum\limits_{j=1}^{N}{{}}Var\left( {{{\tilde{d}}}_{j}} \right)= \\ & \lambda T\cdot \int {{\left[ {{F}_{D}}\left( v \right)-E\left( {{{\tilde{d}}}_{j}} \right) \right]}^{2}}\cdot \tilde{f}\left( v \right)dv+ \\ & E\left( {{{\tilde{d}}}_{j}} \right)E\left( {{D}_{T}} \right)-\lambda T\cdot {{E}^{2}}\left( {{{\tilde{d}}}_{j}} \right). \\ \end{align}$ (9)
N足够大时,式(9)可改写为
$\begin{align} & Var\left( {{D}_{T}} \right)=\underset{N\to \infty }{\mathop{\lim }}\,\sum\limits_{j=1}^{N}{{}}Var\left( {{{\tilde{d}}}_{j}} \right)= \\ & \lambda T\cdot \int F_{D}^{2}\left( v \right)\cdot \tilde{f}\left( v \right)dv.~ \\ \end{align}$ (10)
式(8)和(10)即为本文提出的飓风灾害评估的新方法。与已有方法对比,本文方法中考虑到了飓风发生的间歇性特性,因此更能准确地反映真实情况。
将飓风发生随机过程描述为Poisson过程,则相邻2次飓风的时间间隔Δ也是随机变量,且服从指数分布,其累积密度函数CDF为
$~F\delta =1-exp-\lambda \cdot \delta $ (11)
2 算例分析2.1 数据本文选取美国佛罗里达州Miami-Dade县1900—2010年的登陆飓风数据进行分析[10],如表 1所示(数据来源: http://coast.noaa.gov/hurricanes)。在110年里共登陆27次飓风,其轨迹如图 1所示。
表 1 美国Miami-Dade县1900—2010年登陆飓风记录
时期2010-072006-082005-082004-091999-101992-081976-081972-091971-08
MWS18.00420.57636.00812.86033.43674.58818.00415.43212.860
日期1968-091952-021948-101947-101945-091941-101936-071936-061935-11
MWS15.43223.14846.29638.58054.01241.15225.72020.57641.152
日期1932-081926-091924-101916-081916-111907-091906-061906-091904-10
MWS28.29264.30033.43620.57625.72015.43238.58054.01230.864
?注: MWS为1 min持续最大风速(maximum wind speed),单位为m/s。


表选项






图 1 美国Miami-Dade县1900—2010年登陆飓风轨迹
图选项





将飓风发生随机过程看作平稳Poisson过程,其Poisson强度λ=27/110=0.245。 根据式(11),相邻2次飓风的时间间隔Δ的CDF为
$~F\delta =1-exp\left( -0.245\delta \right).$ (12)
式(12)中的F(δ)如图 2所示。为了验证本文提出的Poisson过程模型的准确性,计算表 1中每相邻2次飓风的发生时间间隔δi,并按升序排序δ1 <δ<2 <…<δ26; 然后将点对(δii /27)绘制于图 2中(i=1,2,…,26)。 借助Chi-Square的一致性检验,得到p值为0.969[11]。 p值的数学意义是指“出现更不符合假定概率模型的样本的几率”,通常当p>0.200时,便认为假定概率分布模型是合理的。由此可知,飓风的发生时间间隔Δ可用指数分布描述其离散性,且拟合效果良好。
图 2 飓风时间间隔的累积密度函数
图选项





计算得到表 1中27次风速的均值为31.8 m/s,标准差为16.19 m/s。 用Weibull分布来拟合${\tilde{v}}$的分布,则${\tilde{v}}$的CDF,${\tilde{F}}$(v)为
$\tilde{F}\left( v \right)=1-exp{{\left[ -\left( \frac{v}{35.9} \right) \right]}^{2.06}}.$ (13)
表 1中的风速按升序排列v1<v2<…<v27; 然后将点对(vii /28)绘制于图 3中(i=1,2,…,27)。 为方便比较,在图 3中也绘制了式(13)中的 ${\tilde{F}}$(v)。 对Weibull分布拟合进行Chi-Square的一致性检验[11],得到p值为0.663>0.200。 这表明,实际记录的飓风风速(1 min持续风速)可以用Weibull分布描述,且拟合效果良好。
图 3 飓风风速的累积密度函数
图选项





根据表 1中的数据可推算出式(4)和(5)中的“虚拟”年最大风速的PDF,fv(v,i)。 用Poisson过程描述飓风发生的随机过程,在时间段(0,T0]内,记飓风发生次数为Q,则
$\begin{align} & Pr\left( Q=k \right)=\frac{\left( \lambda T_{0}^{k} \right)}{k!}exp\left( -\lambda {{T}_{0}} \right), \\ & k=0,1,2,\ldots \\ \end{align}$ (14)
根据式(14)得:
$\begin{align} & F\left( v,{{T}_{0}} \right)=\sum\limits_{k=0}^{\infty }{{}}Pr\left( Q=k\cdot \right){{{\tilde{F}}}^{k}}\left( v \right)= \\ & exp\left[ -\lambda {{T}_{0}}\left( 1-\tilde{F}\left( v \right) \right) \right].~ \\ \end{align}$ (15)
其中,F(vT0)为时间(0,T0]内的最大风速vT0的CDF。 由式(15),vT0的均值为
$~E{{v}_{{{T}_{0}}}}=\int_{0}^{\infty }{{}}v\cdot \frac{\partial F(v,{{T}_{0}})}{\partial v}\cdot dv.$ (16)
用Weibull分布刻画年最大风速,其CDF记为Fv(v,i),则
${{F}_{v}}\left( v,i \right)=1-exp\left[ -{{\left( \frac{v}{{{u}_{1}}} \right)}^{{{\alpha }_{1}}}} \right].$ (17)
为拟合式(17)中的参数u1α1,令式(16)中的T0=30,40,50,…,100,可分别计算E(v30),E(v40),E(v50),…,E(v100); 然后线性拟合点对ln[E(vT0)],ln(lnT0),则斜率为α1,截距为-α1ln u1。 如图 4所示,拟合得到u1=20.55 m/s,α1=1.27。 进而,由式(17)得:
图 4 拟合式(17)中的参数u1和α1
图选项





$0.0618{{\left( \frac{v}{20.55} \right)}^{0.27}}\cdot exp{{\left[ -\left( \frac{v}{20.55} \right) \right]}^{1.27}}.$ (18)
2.2 灾害评估考虑从当前时刻开始50年内的累积飓风灾害D1D2,…,D50。根据本文提出的式(8)和(10)计算得到DT的均值和标准差,如图 5所示。其中,飓风灾害函数FD(v)采用文[12]中的模型,如式(19)所示。FD(v)的含义是风速为v的飓风引起的损伤资产占该地区所有资产的比重,
$\eqalign{ & {F_D} = \cr & \left\{ \matrix{ 0.01exp\left( {0.252v - 5.823} \right),v \le 35; \hfill \cr 0.01\left[ {20 + 11.43\left( {v - 35} \right)} \right],35 < v \le 42; \hfill \cr 1,v > 42. \hfill \cr} \right. \cr} $ (19)
其中,v表示10 min持续风速(m/s)。 表 1中的风速为1 min持续风速,因此,需要将1 min风速乘以转换系数0.9来转化为10 min风速。
为方便对比,分别根据式(4)和(5)计算得到不考虑飓风间歇性特性的累积灾害DT,其均值和标准差也绘制于图 5中。从图 5(a)可以看到,50年内的飓风累积灾害随着时间的增加而线性增加。如果忽略飓风发生过程的间歇性,则会显著高估飓风灾害的均值。以E(D50)为例,由式(8)计算得E(D50)=3.332,而由式(4)(即不考虑飓风间歇性)计算得E(D50)=4.713,高估41.4%。 从图 5(b)可以看到,如果忽略飓风的间歇性特性,则会高估飓风灾害的方差。比如,对于Var(D50),由本文方法计算得Var(D50)=2.823,而由式(5)计算得Var(D50)=3.389,高估20.05%。 2种方法计算得到的DT的变异系数COV,即标准差与均值的比值如图 5(c)所示。可以看到,随着T增加,DT的变异性变小。这是因为2种方法中,DT的均值随T线性增加,但其标准差则随T平方根线性增加。此外,如果忽略飓风的间歇性特性,则会低估飓风灾害的变异性。
图 5 50年内的累积飓风灾害
图选项





3 结 论本文考虑到飓风的间歇性特征,采用平稳Poisson过程来描述飓风发生的随机过程,提出了飓风灾害评估的新方法。对美国Miami-Dade县1900—2010年的飓风数据进行了分析,结果表明: 平稳Poisson过程可以准确描述飓风发生的随机过程; 此外,每次飓风的风速可以用Weibull分布来刻画。应用本文方法对Miami-Dade县未来50年内的飓风灾害进行了评估,结果表明: 利用已有方法(忽略飓风的间歇性特性),则会高估飓风灾害的均值和方差,但会低估其变异性。在50年内的评估期内,已有方法会对飓风灾害的均值和方差分别高估41.4%和20.05%。

参考文献
[1] Journal of Central South University(Science and Technology), 41(2):649-654.-->Lin N, Emanuel K, Oppenheimer M, et al. Physically based assessment of hurricane surge threat under climate change[J]. Nature Climate Change, 2012, 2(6) : 462–467.DOI:10.1038/nclimate1389
[2] Journal of Central South University(Science and Technology), 41(2):649-654.-->杨林, 曹春荣, 林秋, 等. 台风"苏力"灾害风险与防御行为效益评估[J]. 气象与环境学报, 2015, 31(1) : 106–111.YANG Lin, CAO Chunrong, LIN Qiu, et al. Study on disaster risk of typhoon Soulik and its prevention benefit assessment[J]. Journal of Meteorology and Environment, 2015, 31(1) : 106–111.(in Chinese)
[3] Journal of Central South University(Science and Technology), 41(2):649-654.-->Huang Z, Rosowsky D V, Sparks P R. Long-term hurricane risk assessment and expected damage to residential structures[J]. Reliability Engineering and System Sa[KG-0.2mm]fety, 2001, 74 : 239–249.DOI:10.1016/S0951-8320(01)00086-2
[4] Journal of Central South University(Science and Technology), 41(2):649-654.-->Li Y, Ellingwood B R. Hurricane damage to residential construction in the US:Importance of uncertainty modeling in risk assessment[J]. Engineering Structures, 2006, 28 : 1009–1018.DOI:10.1016/j.engstruct.2005.11.005
[5] Journal of Central South University(Science and Technology), 41(2):649-654.-->Bjarnadottir S, Li Y, Stewart M G. A probabilistic-based framework for impact and adaptation assessment of climate change on hurricane damage risks and costs[J]. Structural Safety, 2011, 33 : 173–185.DOI:10.1016/j.strusafe.2011.02.003
[6] Journal of Central South University(Science and Technology), 41(2):649-654.-->Li Y, Stewart M G. Cyclone damage risks caused by enhanced greenhouse conditions and economic viability of strengthened residential construction[J]. Natural Hazard Review, 2011, 12 : 9–18.DOI:10.1061/(ASCE)NH.1527-6996.0000024
[7] Journal of Central South University(Science and Technology), 41(2):649-654.-->Mo H M, Hong H P, Fan F. Estimating the extreme wind speed for regions in China using surface wind observations and reanalysis data[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2015, 143 : 19–33.DOI:10.1016/j.jweia.2015.04.005
[8] Journal of Central South University(Science and Technology), 41(2):649-654.-->Klima K, Lin N, Emanuel K, et al. Hurricane modification and adaptation in Miami-Dade county,Florida[J]. Environmental Science and Technology, 2012, 46 : 636–642.DOI:10.1021/es202640p
[9] Journal of Central South University(Science and Technology), 41(2):649-654.-->Emanuel K. Divine Wind:The History and Science of Hurricanes[M].New York: Oxford University Press, 2005.
[10] Journal of Central South University(Science and Technology), 41(2):649-654.-->Vickery P J, Twisdale L A. Prediction of hurricane wind speed in the United States[J]. Journal of Structural Engineering, 1995, 121(11) : 1691–1709.DOI:10.1061/(ASCE)0733-9445(1995)121:11(1691)
[11] Journal of Central South University(Science and Technology), 41(2):649-654.-->葛余博. 概率论与数理统计[M].北京: 清华大学出版社, 2005.GE Yubo. Probability and Statistics.Beijing:Tsinghua University Press[M].2005.(in Chinese)
[12] Journal of Central South University(Science and Technology), 41(2):649-654.-->Stewart M G, Rosowsky D V, Huang Z. Hurricane risks and economic viability of strengthened construction[J]. Natural Hazards Review, 2003, 4 : 12–19.DOI:10.1061/(ASCE)1527-6988(2003)4:1(12)

相关话题/过程 计算

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 快速的飞行全过程航迹预测
    邹翔,程朋,程农清华大学自动化系,北京100084收稿日期:2015-11-25基金项目:国家自然科学基金资助项目(61179052);“九七三”基础科学研究项目(2010CB731805)作者简介:邹翔(1988-),男,博士研究生摘要:该文提出一种面向空中交通管理的飞行全过程快速4-D航迹预测方 ...
    本站小编 Free考研考试 2020-04-15
  • 平面柔性铰链机构的柔度计算方法
    杜云松,李铁民,姜峣,张京雷清华大学机械工程系,北京100084收稿日期:2015-07-02基金项目:国家科技重大专项(2015ZX04001002)作者简介:杜云松(1983-),男,博士研究生。通讯作者:李铁民,副研究员,E-mail:litm@mail.tsinghua.edu.cn摘要:作 ...
    本站小编 Free考研考试 2020-04-15
  • Fuzzing过程中的若干优化方法
    马金鑫1,张涛1,李舟军2,张江霄31.中国信息安全测评中心,北京100085;2.北京航空航天大学计算机学院,北京100191;3.邢台学院数学与信息技术学院,邢台054001收稿日期:2016-01-22基金项目:国家自然科学基金资助项目(61502536);河北省高等学校青年拔尖人才计划项目( ...
    本站小编 Free考研考试 2020-04-15
  • 现代教育技术专业是否有计算机方面要求
    提问问题:现代教育技术专业是否有计算机方面要求学院:教育学部提问人:15***09时间:2019-09-2309:07提问内容:请问老师,现代教育技术专业,是否有计算机方面的要求,需要计算机二级证书吗?回复内容:看初试考试科目吧 ...
    本站小编 天津师范大学 2019-11-27
  • mpacc最终成绩计算方式以及复试科目
    提问问题:mpacc最终成绩计算方式以及复试科目学院:管理学院提问人:18***83时间:2018-09-2111:10提问内容:老师您好,我想问下mpacc最终成绩计算方式以及复试科目回复内容:请咨询88364386管理学院 ...
    本站小编 山东大学 2019-11-26
  • mpacc最终成绩计算方式以及复试科目
    提问问题:mpacc最终成绩计算方式以及复试科目学院:提问人:18***83时间:2018-09-2111:11提问内容:老师您好,我想问下mpacc最终成绩计算方式以及复试科目回复内容:复试方案可咨询管理学院 ...
    本站小编 山东大学 2019-11-26
  • 专业目录没找到非全的计算机技术专业
    提问问题:专业目录没找到非全的计算机技术专业学院:计算机科学与技术学院提问人:18***48时间:2018-09-1915:44提问内容:我想查非全日制的计算机技术专业,但是在研招网的系统中没有找到该专业,何解?回复内容:你好,请查看计算机学院专业目录:http://www.yz.sdu.edu.c ...
    本站小编 山东大学 2019-11-26
  • 请问上海工作的个人能否报读计算机专业的单独考试
    提问问题:请问上海工作的个人能否报读计算机专业的单独考试学院:提问人:54***om时间:2017-09-1916:28提问内容:请问上海工作的个人能否报读计算机专业的单独考试回复内容:招生简章中有报考条件 ...
    本站小编 山东大学 2019-11-26
  • 计算机技术非全日制双证工程硕士
    提问问题:计算机技术非全日制双证工程硕士学院:提问人:ah***om时间:2016-09-2312:49提问内容:老师,您好!贵校招收计算机技术专业“非全日制双证工程硕士”么?如果招收,上课地点是在哪?山大本部么?上课方式是周末上课还是短期集中上课?我在江苏南通。谢谢回复内容:详情请参考招生专业目录 ...
    本站小编 山东大学 2019-11-26
  • 计算机类非全日制招生
    提问问题:计算机类非全日制招生学院:提问人:18***66时间:2016-09-1916:16提问内容:您好,请问计算机类非全日制今年有招生计划吗,有的话收费是怎么样的?回复内容:有,学费还不确定 ...
    本站小编 山东大学 2019-11-26