基于平衡运动方程的台风风场模型评述.docxVIP

基于平衡运动方程的台风风场模型评述 台风灾害是自然界最严重的灾难之一。从世界范围来看,台风致灾因子1(风压引起的极值风速和风暴潮以及降雨引起的洪水和山体滑坡)对承灾体(城市或区域,土木工程基础设施系统或土木工程结构)造成的灾害及带来的损失惨重。就中国而言,沿海城市因热带气旋造成的损失高达25亿元人民币。2 008年最强台风黑格比(Hagupit)以15m/s登陆广东省,多个城市潮位为百年一遇,造成广东、广西和海南三省(区)11 00余万人受灾,因灾直接经济损失134亿元。台风风险分析是城市区域进行抗风设计和规划、防台减灾的必要途径。 台风风险分析,包括台风危险性分析及承灾体的易损性分析、台风灾害的损失评估等。其中危险性分析是指,根据风的气象特征以及风的历史记录预测某一特定地点或区域在一定时期内可能遭受的强风的概率。风危险性分析具有很强的统计特性,分析结果受样本的数量和质量的影响。对于良态风,某地点大量的可靠的观测资料是危险性分析的基础。然而对于台风,这点无法实现。其一,低频度、小尺度台风特征,导致某地点遭遇台风为小概率事件;其二,台风条件下观测仪器的破坏导致某地点观测数据的质量无法保证和数据丢失;而且气象台站分布间距大,风速记录间距长,风速观测资料序列较短。 鉴于以上两个原因,为了寻求台风气候下风危险性分析的基础,目前国际上普遍流行的做法是,首先利用物理模型来构建台风风场,并通过有限的可靠的观测资料验证构建的台风风场;然后根据台风的关键参数(表征台风主要特征的参数)概率分布,利用MonteCarlo模拟方法产生大量的虚拟台风,最后将以关键参数表达的虚拟台风输入台风风场模型,通过数值模拟来扩充某地点台风的数量。此做法已经在多国风荷载规范中得到了应用,即美国(ANSIA5 8.1-1 982;ASCE7-8 8;ASCE7-9 3;ASCE7-9 5;ASCE7-9 8;ASCE7-02)和澳大利亚(SAA,1989)。尽管此法已经得到应用,各国学者为了改进分析方法一直在努力。自Russell在1968年首次发表文章以来,基于数值模拟的台风危险性分析方法的总体思路一直没有改变,只是随着台风观测数据的质量和数量的提升,台风风场模型方面有所改进,主要特点是采用越来越复杂的物理模型即大气运动方程。 台风风场模型是基于数值模拟的台风危险性分析方法的重要组成部分。气象台风模式中考虑了比较复杂的物理过程,其模拟预测的风场更为真实,并能够模拟得到细致的台风风场分布结构。但是,这种模型用于MonteCarlo数值模拟来进行台风危险性分析则过于复杂,计算量极大;而相对简化但仍然可以描述台风风场主要特征的模型则在工程台风数值模型(相对于气象台风而言,指工程抗风减灾及保险等领域应用的台风模型)中受到欢迎。台风风场包括参数化的台风风场模型和边界层模型,本文将对它们进行评述,探讨模型关键参数的标定,期望为我国沿海地区台风危险性分析的数值台风风场提供技术基础。 1 成熟台风风场的模型构建 考虑到台风危险性分析的计算及工程应用的方便,有必要用有限的参数来描述台风的结构和台风的发生过程,同时保证简化模型能够描述成熟台风风场的主要特征。各国学者相继提出了一系列台风风场模型,包括基于台风梯度平衡方程的Batts模型,基于Shapiro数值风场的Georgiou模型,基于Chow数值模型的CE模型以及Vickery风场模型。 1.1 风场中的空气微团运动平衡方程 1.1.1 柱系温度风速控制方程 对于无加速度的气流,在足够高的高度处(通常被定义为梯度高度),地面摩擦对其影响可以忽略,其水平运动由压力梯度、科氏力和离心力的平衡关系决定,在柱坐标系下的方程可以表示为: 式中:f是科氏参数;ρ是空气密度;Vgr是梯度高度风速;压梯度;r是空气微团的曲率半径。 1.1.2 大气平均运动方程 在大气边界层内,以地球为参照系的大气平均运动方程的矢量表达式为 式中:是质点导数;是局部导数;▽是二维微分算子,也是梯度算子;v是边界层内垂直平均的水平风速矢量;是对流项;f为科氏参数(s-1);k是垂直方向的单位矢量;ρ是平均空气密度;p是大气压;KH是水平涡旋粘性系数;CD是阻力系数;h为边界层厚度。其中大气压p(r)通常采用Holland气压场,即 式中:p0是中心气压;B是Holland气压剖面参数。 大气平均运动方程式是时间和空间坐标的非线性偏微分方程组,得到解析解很困难,通常采用数值求解,其中最常用的是有限差分法和谱方法。有限差分法相对简单,容易理解,应用灵活。谱方法具有无穷阶精度和指数阶收敛性,其作为一种高精度数值方法备受关注。然而,谱方法对偏微分方程真解的“光滑性”要求太强,不能灵活地适应复杂求解区域。Shapiro将运动方程转换为原点固连于台风中心的柱坐标下的运动方程,