海洋科学  2020, Vol. 44 Issue (4): 1-12   PDF    
http://dx.doi.org/10.11759/hykx20190826001

文章信息

郭云霞, 侯一筠, 齐鹏. 2020.
GUO Yun-xia, HOU Yi-jun, QI Peng. 2020.
中国东南沿海区域台风数值模拟与危险性分析
Typhoon wind numerical simulation and risk analysis for southeast coastal region of China
海洋科学, 44(4): 1-12
Marine Sciences, 44(4): 1-12.
http://dx.doi.org/10.11759/hykx20190826001

文章历史

收稿日期:2019-08-26
修回日期:2019-12-05
中国东南沿海区域台风数值模拟与危险性分析
郭云霞1,2,3,5, 侯一筠1,2,3,4,5, 齐鹏1,2,4,5     
1. 中国科学院海洋研究所, 山东 青岛 266071;
2. 中国科学院海洋环流与波动重点实验室, 山东 青岛266071;
3. 中国科学院大学, 北京 100049;
4. 青岛海洋科学与技术试点国家实验室海洋动力过程与气候功能实验室, 山东 青岛 266237;
5. 中国科学院海洋大科学中心, 山东 青岛 266071
摘要:通过网格定点法对我国东南沿海区域性台风危险性进行了分析。利用对各网格点有影响的历史台风数据,建立了各网格点的台风关键参数的最优概率模型。利用Monte-Carlo方法产生每个网格点1 000年间的虚拟台风事件。采用Yan Meng(YM)风场模型模拟了100个历史台风的最大风速,通过使这些最大风速与观测的最大风速误差和最小,建立了一组新的计算最大风半径Rmax和Holland气压参数B的公式,结果表明新的台风参数计算方案效果良好。利用新的参数计算方案、YM风场模型、特定地点的台风衰减模型以及合适的极值分布模型,预测了各个网格点不同重现期的极值风速,进而绘制了台风多发区域的设计风速图。最后研究了不同B模型、Rmax模型和极值分布模型对预测的极值风速的影响。可以为结构抗风设计和台风防灾减灾提供新的参考。
关键词Monte-Carlo模拟    风场模型    台风参数    极值风速    
Typhoon wind numerical simulation and risk analysis for southeast coastal region of China
GUO Yun-xia1,2,3,5, HOU Yi-jun1,2,3,4,5, QI Peng1,2,4,5     
1. Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China;
2. CAS Key Laboratory of Ocean Circulation and Waves, Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China;
3. University of Chinese Academy of Sciences, Beijing 100049, China;
4. Laboratory for Ocean and Climate Dynamics, Pilot National Laboratory for Marine Science and Technology(Qingdao), Qingdao 266237, China;
5. Center for Ocean Mega-Science, Chinese Academy of Sciences, Qingdao 266071, China
Abstract: In this paper, the grid fixed-point method was used to analyze the typhoon risk in the southeast coastal region of China. Using the historical typhoon data for each grid, we established the site-specific statistical distributions of typhoons key parameters. A Monte-Carlo method was used to generate virtual typhoons of 1000 years for each grid. Yan Meng (YM) wind field model was adopted to simulate the maximum wind speeds of 100 historical typhoons. By minimizing the errors between these maximum wind speeds and the observed maximum wind speeds, a new set of formulas was established to calculate the radius to maximum winds (Rmax) and Holland pressure profile parameter (B). The results show that the new parameter calculation scheme works well. Using the new scheme, YM wind field model, region-specific statistical models for the decay rate of typhoons after reaching land, and the extreme value distribution, we predicted the site-specific extreme wind speeds associated with various return periods and proposed a map of design wind speeds for the typhoon-prone coasts of China. Finally we investigated the effects of different B model, Rmax model, and extreme value distributions on the predicted wind speed, which can provide a new reference for wind-resistant structural design and typhoon disaster prevention and mitigation.
Key words: Monte-Carlo simulation    wind field model    typhoon parameters    extreme wind speed    

中国东南沿海区域经济发达、人口稠密, 同时也是遭受台风灾害最为严重的地区之一。统计表明台风灾害造成的损失是所有自然灾害中最多的。历史上台风给我国沿海各省造成了重大经济损失和严重伤亡。在这些地区, 分析台风灾害风险并预测台风的设计风速对于结构抗风设计和台风防灾减灾至关重要。

在中国东南沿海台风多发区域, 由于缺乏足够的台风观测数据, 采用传统的统计方法不能准确预测台风不同重现期的极值风速。为了克服这个问题, Monte-Carlo方法应运而生, 成为一种可替代而且可行的台风危险性分析方法。它利用台风风场模型和台风关键参数统计模型, 采用数值模拟方法进行台风影响区极值风速的预测(又称为台风危险性分析)。美国[1]和澳大利亚[2]使用该方法编制设计风速图。

美国的Russell[3-4]最早采用经验台风模型模拟计算台风产生的极值风速, 并用来估计墨西哥湾的台风风速。自此以后, Batts等[5]、Shapiro[6]、Georgiou[7]、Vickery等[8]、Meng等[9]、Simiu等[10]、Thompson等[11]都对该模拟方法进行了扩展和改进。正如Vickery等[12]所指出的那样, 尽管这些研究的方法是类似的, 但是在风场模型以及衰减模型方面仍然存在比较大的差异。此外在使得台风的气候态特征保持一致的研究区域的大小上以及提取对特定地点有影响的台风事件方法上也存在差异, 有的采用沿岸分段方法[4-5], 有的采用模拟圆方法[7, 12-14]

Vickery等[15]开发了一种新的台风模拟方法-经验路径模型, 该方法对台风的整个路径进行建模, 成为首个可以产生美国沿岸连续变化的台风气候的研究。自此以后, 对于该研究也有很多的改进, 主要包括改进的Holland气压参数B模型、最大风速半径Rmax模型[16-17]以及新的台风衰减模型[18]。Mudd等[19]、Lee等[20]、Rosowsky等[21]都对经验路径模型进行了一些其他的衍生研究。

在过去的几十年, 我国在台风风险分析领域也有许多的研究。欧进萍等[22]利用Batts[5]模型对我国东南沿海重点城市台风风险进行了研究。Xiao等[23]结合Thompson[11]台风模型和Monte-Carlo模拟方法对我国东南沿海地区11个重点城市进行台风风灾分析。Hong等[24]基于Monte-Carlo模拟, 绘制了中国沿海地区台风极值风速等值线图。Li等[25-26]使用全路径模型分析了我国台风灾害, 并对全路径模型进行了简化。这些研究为沿海重点城市的抗风设计提供了有益的参考。然而, 不同的风场模型、台风衰减模型、Rmax模型、B模型以及极值风速概率分布模型会导致风速预测存在差异。

在本文中, 我们采用了不同的风场模型、极值分布模型以及新拟合的Rmax模型和B模型来研究中国东南沿海区域的台风危险性。整个中国东南沿海区域被分成0.25°×0.25°的网格, 对每个网格点采用前人[8, 10, 12]构建的Monte-Carlo模拟方法预测相应的极值风速, 最后编制整个沿海区域的结构设计风速图, 这将会为结构抗风设计提供新的参考。尽管台风引起的风暴潮和洪水也会造成灾害, 但是我们在这里只考虑了最直接的风灾(最大风速)。

1 数据与方法

我们将中国东南沿海台风多发区域作为本文的研究区域。参考Hong等[24], 我们采用网格法来确定研究点, 网格的分辨率是0.25°×0.25°, 一共有579个网格点覆盖整个研究区域(如图 1所示)。

图 1 中国东南沿海区域网格研究站点分布图 Fig. 1 Distribution of grid sites in southeast coastal region of China

本研究的基本台风数据来自CMA-STI西北太平洋热带气旋最佳路径数据集(1949—2015年)。数据集包括台风每6 h间隔的位置与强度信息, 具体包括台风名称、编号、台风中心位置(经度与纬度)、中心压强、近台风中心的2 min平均最大风速。

为了解决台风风速观测数据匮乏的问题, Monte- Carlo方法作为可替代以及可行的方法被广泛用来分析台风的危险性, 采用此方法可以产生大量虚拟台风事件。

2 结果与讨论 2.1 台风关键参数与概率分布

台风模型通常用一组参数进行参数化, 称之为台风关键参数, 主要包括台风年发生率λ, 台风中心气压差Δp, 最大风速半径Rmax, 台风中心移动速度VT, 台风中心移动方向θ, 台风移动路径与研究点间的最小距离Dmin, 这些参数表征了台风的结构和强度。对于某一站点, 当台风距该站点的距离满足特定的要求时(模拟圆法), 台风被提取出来, 以这些提取的台风信息为基础建立台风关键参数的概率模型。Li等[25-26]和Vickery等[17]通过参数化实验, 建议模拟圆半径采用250 km。通过统计分析, 每个网格点有25~267个历史台风样本。Vickery等[12]、Simiu等[10]、Yasui等[27]介绍了各个台风关键参数可能符合的概率模型, 如表 1所示。本文采用卡方检验与Kolmogorov– Smirnov (KS)检验两个非参数检验方法对台风关键参数统计特征的拟合优度进行检验, 检验的显著性水平设置为5%。如果理论分布没有通过检验, 则采用经验分布。通过拟合分析, 我们确定了每个网格点各个台风关键参数的最优概率模型, 这可以为其他学者的研究提供一些有益的参考。

表 1 台风关键参数备选概率模型 Tab. 1 Alternative probability distributions of typhoon key parameters
关键参数 备选概率分布
λ 泊松分布, 二项分布, 负二项分布
Δp 对数正态分布, 伽马分布, 威布尔分布
VT 对数正态分布, 正态分布, 伽马分布
θ 正态分布, 双峰分布, Von Mises分布
Dmin 均匀分布, 梯形分布
2.1.1 年发生率

台风年发生率λ是指一年中某地受台风侵袭或影响的次数, 这些台风发生在以研究站点为中心的模拟圆里面或是与模拟圆相交。该参数用来决定Monte-Carlo方法抽取样本的容量。年发生率通常用均匀泊松分布、二项分布或者负二项分布来描述。我们对每一个研究站点的年发生率分别用3个概率模型进行拟合, 通过拟合优度检验我们得出每个网格点的最优概率模型如图 2所示。在31°N以北区域, 经验分布比较适用, 而在23°N以南区域, 泊松分布是适用的。二项分布主要适用于23°—31°N之间的区域。

图 2 每个网格点台风年发生率的最优概率分布 Fig. 2 The optimal probability distribution of the typhoon annual occurrence rate for each grid
2.1.2 中心压差

台风中心压差Δp定义为台风系统最外围一条闭合等压线的气压值(Holland[28]建议西北太平洋区域是1 010.0 hPa)与台风中心最低气压值之差, 是描述台风强度的主要参数。通常用对数正态分布、伽马分布和三参数威布尔分布来进行描述。KS检验方法用来检验每个研究点概率模型的拟合优度。图 3给出了每个网格点台风中心压差的最优概率模型。由图 3可以看出对数正态分布适合28°N以北的区域, 而28°N以南的区域主要采用经验分布。

图 3 每个网格点台风中心压差的最优概率分布 Fig. 3 The optimal probability distribution of central pressure difference for each grid
2.1.3 移动速度

台风中心移动速度VT可根据CMA数据集中台风6 h间隔的前后中心位置计算得到。备选概率分布有正态分布、对数正态分布和伽马分布。通过拟合优度检验, 我们发现伽马分布主要是适合24°N以北的区域, 经验分布适合剩余的大部分区域, 如图 4所示。

图 4 每个网格点台风移动速度的最优概率分布 Fig. 4 The optimal probability distribution of translation speed for each grid
2.1.4 移动方向

台风移动方向θ由台风数据中台风相邻的两个中心位置(经纬度)连线计算得到, 规定取正北方向为0°, 顺时针方向旋转为正。我们检验了正态分布、双峰分布和Von Mises分布对每个网格点台风移动方向的拟合优度。图 5展示了每个网格点台风移动方向最优的概率分布模型。对大多数的研究点的台风移动方向都有双峰分布的特征, 这也许是因为台风来自不同的热带气旋种群[12]。剩下的极少数的研究点适合用经验分布来描述。

图 5 每个网格点台风移动方向的最优概率分布 Fig. 5 The optimal probability distribution of storm heading for each grid
2.1.5 最小距离

台风移动路径与研究点间的最小距离Dmin是台风直线路径与研究点之间的垂直距离, 是基于台风中心每6 h间隔的台风记录计算得到。当台风中心的移动方向在研究点的左侧时取正值, 在右侧时取负值。Dmin取值的限值是模拟圆半径。Dmin可以用均匀分布或梯形分布来建模。基于拟合优度检验, 每个网格点的最小距离的最优概率分布模型如图 6所示。24°N以北区域适合采用均匀分布, 24°N以南区域适合采用梯形分布。

图 6 每个网格点台风最小距离的最优概率分布 Fig. 6 The optimal probability distribution of the minimum approach distance for each grid
2.2 Monte-Carlo模拟 2.2.1 虚拟台风的生成

本文对每个网格点模拟了1 000年的台风事件, 在此期间台风发生数目由各网格点的年发生率概率分布(图 2)来确定。首先采用Monte-Carlo方法对每个研究点的ΔpVTθDmin的概率分布(图 3图 6)进行随机抽样, 以此来产生虚拟台风起始点的关键参数。从最优概率分布随机抽取的关键参数需要在合理的范围内: Δp为[0, 135] hpa, VT为[2, 65] km/h, θ为[−180°, 180°], |Dmin|≤250 km。然后确定台风的起始点, 此点是台风路径与模拟圆的交点, 用Dminθ、模拟圆半径以及网格点的位置来确定。接着用移动速度VT与移动方向θ来确定台风每小时间隔的中心位置, 在此过程中台风路径被假定为直线, 由抽取的Dminθ确定。在台风登陆前, 中心压强保持不变, 登陆后采用衰减模型来描述中心压差Δp的变化。最后当台风中心离开模拟圆范围则认为台风消失。该过程示意图如图 7所示。

图 7 虚拟台风过程示意图 Fig. 7 A schematic graph for simulation of the virtual typhoon
2.2.2 衰减模型

一旦台风登陆, 强度将会减弱, 台风的中心压差和风速都会减小。台风衰减的速率与台风本身、登陆的位置以及登陆时的中心压差有关。采用合适的衰减模型描述台风的衰减规律, 对于预测台风登陆后的风速有重要影响。Batts等[5]、Georgiou[7]以及Vickery等[8]采用中心气压差的衰减来描述台风的衰减。他们利用台风中心气压差和登陆时间或登陆以后向内陆移动的距离的相关性表述衰减规律。在这些衰减模型中应用最广泛的是Vickery等[8]发展的衰减函数。本文采用了类似的衰减模型, 基于CMA数据来描述每个网格点的台风衰减规律。衰减率模型如公式(1)所示:

$ \Delta p\left( t \right) = \Delta {p_0}\exp \left( { - at} \right), a = {a_0} + {a_1}\Delta {p_0} + \varepsilon $ , (1)

式中Δp(t)是台风登陆后t时刻的中心压差(hPa), Δp0是登陆时刻的中心压差(hPa), a是特定点衰减系数, ε是正态分布的误差项, 均值为0, 标准偏差为σε

通过拟合登陆的台风数据, 我们得到了每个网格点的台风衰减模型的系数。对于每个网格点, 用于拟合衰减模型的台风数量从16到118不等。对于个别的在回归分析中没有通过显著性检验(显著性水平是5%)或是登陆的台风数目非常少的网格点采用离它最近的网格点的衰减模型来代替。每个网格单元衰减模型的系数分布如图 8所示, 作为个例, 文中给出了福州、汕头、广州、湛江的衰减模型拟合示意图, 如图 9所示。

图 8 衰减模型中衰减系数a0a1的分布图 Fig. 8 Illustration of obtained coefficients a0 and a1 for decay model

图 9 衰减常数a与登陆时中心压差Δp0的回归关系 Fig. 9 Decay constant a versus the central pressure difference Δp0 at landfall 注: Mean表示均值; SD表示标准偏差
2.3 风场模型与新的台风参数计算方案

台风风场模型是进行台风数值求解的关键。一般来说, 台风数值模拟精度高, 但是计算量也相应变大, 目前还不适合与随机路径结合使用进行风险评估。而参数化风场模型计算简便, 精度也可以达到要求, 在台风灾害风险评估领域有广泛的应用[29]。本文采用了Yan Meng(YM)风场模型, 是由Meng等[9]提出的解析模型, 该模型是半经验半数值的风场模型, 其数值形式如下:

$ \frac{{\partial {\rm{ }}V}}{{\partial t}} + {\rm{ }}V\nabla {\rm{ }}V = - \frac{1}{p}\nabla p - fkV + F $ , (2)

式中, V代表台风风速; ρ代表空气密度; k代表垂向单位向量; f代表科氏参数; F代表边界层摩擦力。该模型以扰动平衡方程和经过边界层摩擦力修正的压力梯度平衡方程为基础, 引入了可以考虑地貌等复杂因素的“等效粗糙长度”, 能够很好地表示出台风风场中各种力之间的相互关系, 并且能够比较直观地求解梯度风方程, 得到完整的解析解, 这种方法相对简单易解, 而且对于分析台风危险性具有足够的精度, 该模型得到了广泛应用, 如Matsui等[30]、Okazaki等[31]、谢汝强等[32]、赵林等[33]。关于风场模型更详细的内容请参考文献Meng等[9]

2.3.1 新的台风参数计算方案

气压是台风风场模型的重要输入参数, YM风场模型采用Holland[28]气压模型, 其形式如下:

$ p(r) = {p_0} + \Delta p\exp [ - {(\frac{{{R_{\max }}}}{r})^B}], $ (3)

式中, p(r)是距台风中心径向距离为r处的海平面压强; p0是台风中心的气压值; Δp是台风中心的气压差。式中有两个未知的参数, 分别是最大风速半径Rmax和气压剖面参数B。它们对台风风场的模拟起重要作用,但是由于缺乏相应的观测数据, 很难得到它们的精确值, 在我国目前还没有统一的计算它们的公式。本文中我们提出了一个新的计算台风参数RmaxB的方案。基于之前的研究[12, 16-17, 34], Rmax与Δp和纬度ψ相关, BRmax和Δp相关, 或是与Rmaxψ相关。我们总结这些关系列出了4类计算台风风场参数的方案,每个方案都带有需要进一步确定的系数, 如公式(4)—(7)所示:

$ \ln {R_{\max }} = {b_1} + {c_1}\Delta p + {d_1}\psi , B = {e_1} + {f_1}\Delta p + {g_1}{R_{\max }}, $ (4)
$ \ln {R_{\max }} = {b_2} + {c_2}\Delta {p^2} + {d_2}{\psi ^2}, B = {e_2} + {f_2}\Delta p + {g_2}{R_{\max }}, $ (5)
$ \ln {R_{\max }} = {b_3} + {c_3}\Delta {p^2} + {d_3}{\psi ^2}, B = {e_3} + {f_3}{R_{\max }} + {g_3}\psi , $ (6)
$ \ln {R_{\max }} = {b_4} + {c_4}\Delta {p^2} + {d_4}{\psi ^2}, B = {e_4} + {f_4}\Delta p + {g_4}{R_{\max }}, $ (7)

基于前人的研究可以给定每一个未知系数(b, c, d, e, f, g )合理的变化范围。我们提取了不同气象站观测得到的100个台风风速的最大值, 通过使其模拟值与观测值的误差和最小来确定各个待定系数在其变化范围内的最优值。关注最大风速的误差是因为最大风速模拟的准确性直接影响我们预测的极值风速的准确性。100个观测台风的最大风速主要来自“中国热带气旋气候图集”[35]和香港天文台(http://www.hko.gov.hk/contentc.htm)。每个观测台风的路径以及相应观测站点的信息如图 10所示。

图 10 100个观测台风的路径和观测气象站点的位置 Fig. 10 The tracks of 100 observed typhoons and the locations of meteorological stations

我们采用YM风场模型模拟每一个观测台风的风速, 并且计算模拟风速与观测风速最大值之间的误差。以公式(4)为例, 我们首先对每一个系数b1, c1, d1, e1, f1, g1设定了一个合理的变化范围。然后利用六层嵌套循环来依次描述各系数值的变化, 最里面的循环是g1的变化, 最外面的循环是b1的变化。这样我们可以得到20多万组b1, c1, d1, e1, f1, g1的不同值的组合。接下来, 我们将各组合值依次代入公式(4)结合YM风场模型计算各个台风的最大风速, 并与观测的最大风速对比得到相对误差, 进而计算100个观测台风的相对误差之和, 由此我们可以得到二十多万个台风的最大风速误差和, 从中选取误差和最小的那一组系数(b1, c1, d1, e1, f1, g1)作为公式(4)中各个系数的最优值。以此类推, 我们得到了公式(4)—(7)的系数的最优值, 结果如下:

$ \begin{array}{l} \ln {R_{\max }} = 2.8 - 0.002\Delta p + 0.034\psi ,\\ B = 1.2 + 0.01\Delta p - 0.02{R_{\max }}, \end{array} $ (8)
$ \begin{array}{l} \ln {R_{\max }} = 3.4 - 4 \times {10^{ - 5}}\Delta {p^2} + 4 \times {10^{ - 4}}{\psi ^2},\\ B = 1.2 + 0.01\Delta p - 0.002{R_{\max }}, \end{array} $ (9)
$ \begin{array}{l} \ln {R_{\max }} = 3.4 - 4 \times {10^{ - 5}}\Delta {p^2} + 4 \times {10^{ - 4}}{\psi ^2},\\ B = 1.4 - 0.03{R_{\max }} - 0.07\psi , \end{array} $ (10)
$ \begin{array}{l} \ln {R_{\max }} = 3.5 - 3 \times {10^{ - 5}}\Delta {P^2} + 0.02\psi ,\\ B = 1.4 + 0.001\Delta p - 0.002{R_{\max }}. \end{array} $ (11)

确定公式(4)—(7)中系数最优值时100个台风的模拟最大风速误差和的最小值分别是29.33、28.16、28.18、26.06。由此可知, 在公式(8)—(11)中, 采用公式(11)计算的100个台风的最大风速误差和最小。因此我们选取公式(11)作为计算台风参数的最优方案, 这为计算气压剖面参数以及最大风速半径提供了新的参考。为了进一步验证这个新的参数计算方案, 我们还采用了RmaxB的另外两个常用的模型, 一个来自Vickery等[15], 模型形式如下:

$ \begin{array}{l} \ln {R_{\max }} = 2.636 - 0.0000508\Delta {P^2} + 0.0394\psi ,\\ B = 1.38 + 0.0184\Delta p - 0.00309{R_{\max }}, \end{array} $ (12)

另一个模型来自Vickery等[16], 形式如下:

$ \begin{array}{l} \ln {R_{\max }} = 3.015 - 6.291 \times {10^{ - 5}}\Delta {P^2} + 0.0337\psi + {\varepsilon _{\ln {R_{\max }}}},\\ B = 1.833 - 0.326\sqrt {1000{f_c}{R_{\max }}} + {\varepsilon _B}, \end{array} $ (13)

式中Δp单位是hPa, ${\varepsilon _{\ln {R_{\max }}}}$是标准偏差, 当$\Delta p \le 87{\rm{hPa, }}{\sigma _{\ln {R_{\max }}}} = 0.448$; 当$87{\rm{hPa < }}\Delta p \le 120{\rm{hPa}}, {\sigma _{\ln {R_{\max }}}} = 1.137 - 0.00792\Delta p$, 当$\Delta p \le 120{\rm{hPa, }}{\sigma _{\ln {R_{\max }}}} = 0.186;\psi $是纬度(单位是度); fc是科氏参数; ${\sigma _B} = 0.221$。同样我们采用公式(12)与(13)结合YM风场模型计算了100个台风的最大风速误差和, 分别为29.06、27.88。可以看出采用公式(12)或(13)得到的误差和也大于采用公式(11)的结果。在图 11中, 我们绘制了采用公式(12)计算得到的每个台风最大风速的相对误差与采用公式(11)所得结果的差值以及公式(13)和公式(11)所得结果的差值。可以看出, 对于大多数台风, 采用公式(12)或公式(13)计算的最大风速的相对误差大于公式(11)的结果, 说明我们的新方案可以有效的降低台风最大风速的模拟误差。

图 11 公式(12)与公式(11), 公式(13)与公式(11)计算的各台风最大风速相对误差的差值 Fig. 11 The difference of relative error of maximum wind speed for each typhoon between the Eq. (12) and Eq. (11) and that between the Eq. (13) and Eq. (11)

为了对公式(11)进行进一步验证, 我们将其带入YM模型, 模拟了上述100个台风中York(1999)、Usagi(2013)、Hagupit(2008)、Hato(2017)的风速, 模拟和实测风速的对比结果如图 12所示。由此可见, 新的台风参数计算方案效果良好。

图 12 将式(11)代入YM模型时, 台风York、Usagi、Hagupit以及Hato的模拟风速和实测风速的对比 Fig. 12 The simulated and measured wind speeds for typhoons York, Usagi, Hagupit, and Hato when Eq. (11) was brought into the Yan Meng (YM) model
2.3.2 粗糙度的确定

在Meng等[9]的研究中, 地形效应用等效地面粗糙度来代替。粗糙度对于风场模型也是非常重要的参数。基于之前的研究[2, 36-37]不同地形地貌的粗糙度的取值如表 2所示。我们根据实际地形标记了每个网格点的粗糙度, 如图 13所示。

表 2 不同地形地貌的粗糙度的值 Tab. 2 Roughness length of different geomorphology
地形等级 下垫面特征 粗糙度/m
I 海面、泥滩、冰雪覆盖的平原, 无障碍海岸地区等 0.000 5~0.003
II 平坦开阔的田野、乡村及丛林等(气象学标准) 0.003~0.2
III 丘陵和房屋比较稀疏的乡镇及城市郊区 0.2~1.0
IV 密集建筑群的城市 1.0~2.0
V 密集建筑群且房屋较高的城市 2.0~4.0

图 13 不同研究点的粗糙度 Fig. 13 The roughness length for each site of interest
2.4 台风危险性分析

基于本文得到的台风关键参数概率分布, 采用Monte-Carlo模拟方法得到的虚拟台风事件, 结合采用最优台风参数计算方案的YM风场模型, 可以计算得到每个网格点的台风极值风速序列。YM风场模型计算的风速为10 m高度处每小时平均风速, 采用因子1.06将其转换成10 min平均的风速。极值风速序列通常会用理论分布来进行描述, 本文中采用了应用广泛的Weibull分布和Gumbel分布, 并进行了KS拟合优度检验, 检验显著性水平为5%。如果某个研究点的理论分布都没有通过显著性检验则采用经验分布。如果两个理论分布都通过了显著性检验则根据p值选择最优的理论分布。

首先, 我们利用经验分布来预测每个网格单元不同重现期的极值风速。图 14显示了预测的每一个网格点50年及100年重现期的极值风速。从风速等值线可以看出, 预测的极值风速从沿海地区到内陆地区逐渐减小。在浙江南部、广东中部和南部沿海地区、珠江三角洲和长江三角洲地区预测的极值风速普遍比较大, 我们需要对这些地区的台风灾害给予更多的关注, 并投入更多的资源研究这些地区高层建筑的抗风性能。

图 14 采用经验分布预测的50年与100年重现期的极值风速 Fig. 14 The predicted extreme wind speed of the 50a and 100a return period using empirical distribution

在我国《建筑结构荷载规范》[38]中, 将许多地方的参考风压列制成表。Hong等[24]对设计规范中给出的等值线图进行数字化处理, 得到了50年重现期的极值风速图, 并通过Monte-Carlo模拟预测了中国大陆沿海地区年最大风速的50年重现期值。然而Hong等[24]使用的是改进的Thompson风场模型[11]以及Vickery等[16]发展的Rmax模型和B 模型。将Hong等[24]提出的等值线与图 14所示的等值线进行对比, 可以看出设计规范、Hong等[24]和本研究得到的设计风速的空间趋势是相似的。本研究预测风速比Hong等[24]的结果大1~2 m/s。与设计规范的结果相比, 本研究预测风速在沿海地区比规范值大1~2 m/s, 在内陆地区比规范值大4~5 m/s。从海岸到内陆, 本研究预测的极值风速的下降速度要低于设计规范中风速的下降速度, 这与Hong等[24]的结果相似。由于设计规范不仅考虑了台风, 还考虑了其他较弱的风场, 所以很容易低估极值风速。

其次, 为了研究不同Rmax模型和B模型对预测风速的影响, 我们也采用了Vickery等[15]和Vickery等[16]中的Rmax模型和B模型来计算台风风场, 如公式(12)和(13)。基于公式(12)和(13)计算的风速, 我们预测了中国东南沿海区域50年重现期的极值风速。图 15a是基于式(12)预测的极值风速减去基于式(11)预测的风速的结果。图 15b是基于式(13)预测的风速减去基于式(11)预测的风速的结果。由图 15a可以看出, 采用式(12)预测的沿海地区风速比采用式(11)预测的风速大0~2 m/s。在内陆地区, 采用式(12)预测风速略小于采用式(11)预测的风速。从图 15b, 可以看出30°N以北区域, 采用式(13)预测的风速通常比用式(11)预测的风速小0~3 m/s, 而26°N以南区域, 采用式(13)预测风速通常比用式(11)预测的风速大0~2 m/s。

图 15 基于公式(12)和(13)预测的50年重现期极值风速减去基于式(11)预测的风速的差值分布图(a, b) Fig. 15 The results that the predicted-wind-speeds of 50-year return period based on the Eq. (12) (a) and Eq. (13) (b) minus those based on the Eq. (11)

最后, 我们研究了不同极值分布对预测风速的影响。通过拟合最优极值分布, 我们发现每个网格研究点的最优极值分布不是Weibull分布就是Gumbel分布。因此, 我们比较了理论分布与经验分布预测风速的差异, 如图 16所示。图 16中正值表示理论分布预测风速大于经验分布预测风速。我们可以看到, 在中间区域(24°—30°N, 主要是Gumbel分布)用理论分布预测的风速通常比使用经验分布预测的结果大0~2 m/s。在两边区域(主要是Weibull分布), 使用理论分布预测的风速通常比使用经验分布预测的结果小, 尤其是在23°N以南的区域。

图 16 采用理论分布预测的50年重现期极值风速减去经验分布预测的风速的差值分布图 Fig. 16 The results that the predicted-wind-speeds of 50-year return period using theoretical distribution minus those using empirical distribution
3 结论

采用基于新的台风参数计算方案的YM风场模型以及Monte-Carlo模拟我们分析了中国东南沿海区域台风的危险性。整个东南沿海台风多发区被分成0.25°×0.25°的网格, 每个网格点作为我们的研究点。基于历史数据, 我们对每个研究点的每个台风关键参数建立了最优概率模型。基于每个研究点台风关键参数的概率模型, 我们用Monte-Carlo方法产生了每个研究点1 000年的虚拟台风事件。通过使模拟与观测的台风最大风速的误差和达到最小, 我们确定了台风参数计算的新方案。采用YM风场模型结合台风参数计算的新方案、特定点的台风衰减模型以及最优的极值分布, 我们预测了每个研究点不同重现期的极值风速, 并绘制了新的中国东南沿海区域的设计风速图。最后我们研究了不同B模型、Rmax模型和极值分布对极值风速预测的影响, 为抗风结构设计提供了新的参考。

参考文献
[1]
ASCE/SEI 7-05, Minimum Design Loads for Buildings and other Structures[S]. http://faculty.arch.tamu.edu/media/cms_page_media/4198/NS13-6asce7_v0JLEKW.pdf
[2]
AS/NZS 1170.2: 2002, Structural Design Actions part2: Wind Actions[S]. http://www.doc88.com/p-33973362933901.html
[3]
Russell L B. Probability distributions for Texas gulf coast hurricane effects of engineering interest[D]. Stanford: Stanford University, 1969. http://www.oalib.com/references/17239495
[4]
Russell L R. Probability distributions for hurricane effects[J]. Journal of the Waterways, Harbors and Coastal Engineering Division, 1971, 97(1): 139-154.
[5]
Batts M E, Simiu E, Russell L R. Hurricane wind speeds in the United States[J]. Journal of the Structural Division, 1980, 106(10): 2001-2016.
[6]
Shapiro L J. The Asymmetric boundary layer flow under a translating hurricane[J]. Journal of the Atmospheric Sciences, 1983, 40(8): 1984-1998. DOI:10.1175/1520-0469(1983)040<1984:TABLFU>2.0.CO;2
[7]
Georgiou P N. Design wind speeds in tropical cyclone-prone regions[D]. London, Canada: University of Western Ontario, 1986. https://www.sciencedirect.com/science/article/pii/0167610583901368
[8]
Vickery P J, Twisdale L A. Wind-field and filling models for hurricane wind-speed predictions[J]. Journal of Structural Engineering, 1995b, 121(11): 1700-1709. DOI:10.1061/(ASCE)0733-9445(1995)121:11(1700)
[9]
Meng Y, Matsui M, Hibi K. An analytical model for simulation of the wind field in a typhoon boundary layer[J]. Journal of Wind Engineering and Industrial Aerodynamics, 1995, 56(2-3): 291-310.
[10]
Simiu E, Scanlan R H. Wind Effects on Structures:Fundamental and Applications to Design[M]. New York: John Wiley & SONS, INC, 1996.
[11]
Thompson E F, Cardone V J. Practical modeling of hurricane surface wind fields[J]. Journal of Waterway, Port, Coastal, and Ocean Engineering, 1996, 122(4): 195-205. DOI:10.1061/(ASCE)0733-950X(1996)122:4(195)
[12]
Vickery P J, Twisdale L A. Prediction of hurricane wind speeds in the United States[J]. Journal of Structural Engineering, 1995, 121(11): 1691-1699. DOI:10.1061/(ASCE)0733-9445(1995)121:11(1691)
[13]
Georgiou P N, Davenport A G, Vickery B J. Design wind speeds in regions dominated by tropical cyclones[J]. Journal of Wind Engineering and Industrial Aerodynamics, 1983, 13(1-3): 139-152. DOI:10.1016/0167-6105(83)90136-8
[14]
Neumann C J. The National Hurricane Center Risk Analysis Program (HURISK)[R]. Miami: U.S. Department of Commerce, National Oceanic and Atmospheric Administration, National Hurricane Center, 1987. https://www.nhc.noaa.gov/data/publications.php?text
[15]
Vickery P J, Skerlj P F, Twisdale L A. Simulation of hurricane risk in the U.S. using empirical track mode[J]. Journal of Structural Engineering, 2000, 126(10): 1222-1237. DOI:10.1061/(ASCE)0733-9445(2000)126:10(1222)
[16]
Vickery P J, Wadhera D. Statistical models of Holland pressure profile parameter and radius to maximum winds of hurricanes from flight-level pressure and H*Wind data[J]. Journal of Applied Meteorology and Climatology, 2008, 47(10): 2497-2517. DOI:10.1175/2008JAMC1837.1
[17]
Vickery P J, Wadhera D, Twisdale L A, et al. U.S. hurricane wind speed risk and uncertainty[J]. Journal of Structural Engineering, 2009, 135(3): 301-320. DOI:10.1061/(ASCE)0733-9445(2009)135:3(301)
[18]
Vickery P J. Simple empirical models for estimating the increase in the central pressure of tropical cyclones after landfall along the coastline of the United States[J]. Journal of Applied Meteorology, 2005, 44(12): 1807-1826. DOI:10.1175/JAM2310.1
[19]
Mudd L, Wang Y, Letchford C, et al. Assessing climate change impact on the U.S. east coast hurricane hazard:Temperature, frequency, and track[J]. Natural Hazards Review, 2014, 15(3): 04014001. DOI:10.1061/(ASCE)NH.1527-6996.0000128
[20]
Lee K H, Rosowsky D V. Synthetic hurricane wind speed records:development of a database for hazard analyses and risk studies[J]. Natural Hazards Review, 2007, 8(2): 23-34. DOI:10.1061/(ASCE)1527-6988(2007)8:2(23)
[21]
Rosowsky D V, Mudd L, Letchford C. Assessing climate change impact on the joint wind-rain hurricane hazard for the northeastern U.S. coastline[C]//Gardoni P, Murphy C, Rowell A. Risk Analysis of Natural Hazards: Interdisciplinary Challenges and Integrated Solutions. Cham: Springer, 2016: 79-92. https://link.springer.com/chapter/10.1007%2F978-3-319-22126-7_8
[22]
欧进萍, 段忠东, 常亮. 中国东南沿海重点城市台风危险性分析[J]. 自然灾害学报, 2002, 11(4): 9-17.
Ou Jinping, Duan Zhongdong, Chang Liang. Typhoon risk analysis for key coastal cities in southeast China[J]. Journal of Natural Disasters, 2002, 11(4): 9-17. DOI:10.3969/j.issn.1004-4574.2002.04.002
[23]
Xiao Y F, Duan Z D, Xiao Y Q, et al. Typhoon wind hazard analysis for southeast China coastal regions[J]. Structural Safety, 2011, 33(4-5): 286-295. DOI:10.1016/j.strusafe.2011.04.003
[24]
Hong H P, Li S H, Duan Z D. Typhoon wind hazard estimation and mapping for coastal region in mainland China[J]. Natural Hazards Review, 2016, 17(2): 4016001. DOI:10.1061/(ASCE)NH.1527-6996.0000210
[25]
Li S H, Hong H P. Observations on a hurricane wind hazard model used to map extreme hurricane wind speed[J]. Journal of Structural Engineering, 2014, 141(10): 04014238.
[26]
Li S H, Hong H P. Typhoon wind hazard estimation for China using an empirical track model[J]. Natural Hazards, 2016, 82(2): 1009-1029. DOI:10.1007/s11069-016-2231-2
[27]
Yasui H, Ohkuma T, Marukawa H, et al. Study on evaluation time in typhoon simulation based on Monte Carlo method[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2002, 90(12-15): 1529-1540. DOI:10.1016/S0167-6105(02)00268-4
[28]
Holland G J. An Analytic model of the wind and pressure profiles in hurricanes[J]. Monthly Weather Review, 1980, 108(8): 1212-1218. DOI:10.1175/1520-0493(1980)108<1212:AAMOTW>2.0.CO;2
[29]
方伟华, 林伟. 面向灾害风险评估的台风风场模型研究综述[J]. 地理科学进展, 2013, 32(6): 852-867.
Fang Weihua, Lin Wei. A review on typhoon wind field modeling for disaster risk assessment[J]. Progress in Geography, 2013, 32(6): 852-867.
[30]
Matsui M, Ishihara T, Hibi K. Directional characteristics of probability distribution of extreme wind speeds by typhoon simulation[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2002, 90(12-15): 1541-1553.
[31]
Okazaki T, Watabe H, Ishihara T. Development of typhoon simulation model for insurance risk estimation[C]//International Association for Wind Engineering (IAWE). The Sixth Asia-Pacific Conference on Wind Engineering. Seoul, Korea: IAWE, 2005: 1790-1802. https://www.researchgate.net/publication/252545338_Development_of_Typhoon_Simulation_Model_for_Insurance_Risk_Estimation
[32]
谢汝强, 李利孝, 王艳华, 等. 广东省台风数值模拟与危险性分析[J]. 安徽建筑大学学报, 2015, 23(4): 51-55.
Xie Ruqiang, Li Lixiao, Wang Yanhua, et al. Typhoon wind numerical simulation and hazard analysis for Guangdong Province[J]. Journal of Anhui Jianzhu University, 2015, 23(4): 51-55.
[33]
赵林, 朱乐东, 葛耀君, 等. 上海地区台风风特性Monte-Carlo随机模拟研究[J]. 空气动力学学报, 2009, 27(1): 25-31.
Zhao Lin, Zhu Ledong, Ge Yaojun, et al. Monte-Carlo simulation about typhoon extreme value wind characteristics in Shanghai region[J]. Acta Aerodynamica Sinica, 2009, 27(1): 25-31. DOI:10.3969/j.issn.0258-1825.2009.01.005
[34]
Powell M, Soukup G, Cocke S, et al. State of Florida hurricane loss projection model:Atmospheric science component[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2005, 93(8): 651-674. DOI:10.1016/j.jweia.2005.05.008
[35]
中国气象局上海台风研究所.中国热带气旋气候图集(1951-2000)[M].北京: 科学出版社, 2007.
Shanghai Typhoon Institute, China Meteorological Administration. Climatological Atlas for Tropical Cyclones Affecting China (1951-2000)[M]. Beijing: Science Press, 2007.
[36]
Counihan J. Adiabatic atmospheric boundary layers:A review and analysis of data from the period 1880-1972[J]. Atmospheric Environment, 1976, 9(10): 871-905.
[37]
Choi H, Kanda J. Characteristics of the vertical wind profile for wind load estimation[J]. Journal of Wind Engineering, 2010, 1990(45): 23-43.
[38]
GB50009-2012, 建筑结构荷载规范[S].
GB 50009-2012, Load code for the design of building structures[S]. http://www.bookask.com/book/130062.html