快速检索
  气象   2019, Vol. 45 Issue (2): 251-262.  DOI: 10.7519/j.issn.1000-0526.2019.02.010

论文

引用本文 [复制中英文]

陈雯超, 刘爱君, 宋丽莉, 等, 2019. 不同强风天气系统风特性的个例分析[J]. 气象, 45(2): 251-262. DOI: 10.7519/j.issn.1000-0526.2019.02.010.
[复制中文]
CHEN Wenchao, LIU Aijun, SONG Lili, et al, 2019. Case Study of Wind Characteristics of Different Strong Wind Systems[J]. Meteorological Monthly, 45(2): 251-262. DOI: 10.7519/j.issn.1000-0526.2019.02.010.
[复制英文]

资助项目

广东省气象局科学技术研究项目(2015B32)、科技部科技基础性工作专项(2013FY112400)和国家自然科学基金青年基金项目(41305007)共同资助

第一作者

陈雯超,从事气候可行性论证和风特性研究.Email:wenchao678@163.com

文章历史

2017年8月24日收稿
2018年9月30日收修定稿
不同强风天气系统风特性的个例分析
陈雯超 1, 刘爱君 1, 宋丽莉 2, 王丙兰 2    
1. 广东省气象防灾技术服务中心,广州 510080
2. 中国气象局公共气象服务中心,北京 100081
摘要:利用广东省东莞市的80 m气象塔在台风莫拉菲、强对流和强冷空气期间的观测数据,分析不同强风天气系统近地面边界层的平均风和脉动风特性参数。研究发现:(1)1号强对流的风速变化最为剧烈,风速的10 min升幅可分别达台风和强冷空气大风的1.8和3.7倍。强对流的10 min平均风向变幅也最剧烈,可达到120°·(10 min)-1。(2)台风强风廓线的幂指数值为0.177,要大于建筑结构荷载规范(GB 5009—2012)给出的B类下垫面0.15的参考值,而强对流的幂指数值接近规范给出的A类下垫面的参考值0.12,强冷空气的幂指数值则仅为0.10。(3)1号强对流强风平均湍流强度分别是台风和强冷空气大风平均湍流强度的1.4~2.1和1.7~2.0倍。(4)台风和强对流的风攻角在-3°~3°范围内, 强冷空气大风的风攻角在0°附近变动。(5)强风影响期间湍流空间积分尺度均有增大的现象,台风强风经过时,水平(纵向和横向)湍流空间积分尺度明显增大了一个量级。(6)台风强风的湍流功率谱不满足在惯性子区湍流能谱密度的-5/3次方与频率成正比的规律和各向同性假设。台风强风的湍流能谱值要显著高于强对流和强冷空气,其中,在桥梁敏感频域台风的平均湍流能谱值约为强对流和强冷空气的3倍。
关键词台风    强对流    强冷空气    风特性    
Case Study of Wind Characteristics of Different Strong Wind Systems
CHEN Wenchao1, LIU Aijun1, SONG Lili2, WANG Binglan2    
1. Guangdong Meteorology Disaster Prevention Technology Service Center, Guangzhou 510080;
2. CMA Public Meteorological Service Centre, Beijing 100081
Abstract: Using the wind data observed from the 80 m meteorological tower in Dongguan of Guangdong Province during the passes of Typhoon Molave, severe convection and severe cold air, the mean wind and fluctuating wind characteristics of different strong wind weather systems in the near-surface boundary layer are analyzed. The results show that: (1) the wind speed variation of the severe convection No.1 is most dramatic and the 10 min mean wind speed amplitude of severe convection can go up to 1.8 times and 3.7 times as much as that of typhoon and severe cold air. The 10 min mean wind direction amplitude can be 120°·(10 min)-1. (2) The power law exponent of typhoon strong wind profile is 0.177 which is larger than 0.15 that is recommended as the power law exponent of exposure class B in Chinese National Code (GB 5009-2012). The power law exponent of severe convection is near 0.12 which is recommended as the exposure class A in the code and the power law exponent of severe cold air is only 0.10. (3) The 10 min average turbulence intensity of No.1 severe convection can go up to 1.4-2.1 times and 1.7-2.0 times as much as that of typhoon and severe cold air. (4) The attack angles of typhoon and severe convection change between -3°-3° while the attack angle of severe cold air changes near 0°. (5) The turbulence integral length scales increase due to the impact of strong wind processes. The horizontal (longitudinal and transverse) turbulence integral length scales of typhoon can obviously increase by an order of magnitude. (6) The typhoon strong wind spectra does not satisfy the -5/3 law and the isotropy assumption. The power spectral density (PSD) va-lue of typhoon is significantly larger than that of severe convection and severe cold air. The PSD value of typhoon is about 3 times as large as that of severe convection and severe cold air.
Key words: typhoon    severe convection    severe cold air    wind characteristics    
引言

随着高耸和大跨等风敏感建筑物、构筑物的不断出现,边界层强风特性也越来越成为风工程界关注的热点(郭冬艳等,2011黄林宏等,2016廖菲等,2017唐明晖等,2017周福等,2017),风廓线、湍流强度、湍流积分尺度和湍流功率谱等反映强风风特性的参数也在工程设计中起着关键的作用(Cao et al, 2009; Wang et al, 2011; 钱燕珍等,2013)。台风、强对流大风和冷空气大风由于生成机理不一样,因此其风特性也会有差别。现场测风是研究近地层风特性的最直接可靠的手段,近年来,我国在近地层风特性研究方面取得了一些重要研究成果。宋丽莉等(20052006)和Song et al(201020122016)采用强风型超声测风仪对登陆广东的多个台风个例进行观测,对比分析了这几个台风的中心区域和其外围风况的平均风廓线和脉动风参数等特征,发现台风过程的中心强风区和外围风况的风廓线和湍流特征差异显著,并建立了粗糙长度与风廓线幂指数的数学模型;庞加斌等(2002)根据在上海观测的“派比安”台风实测风速资料得到了台风的水平风速功率谱和垂直风速功率谱,并对实测水平脉动风谱和垂直脉动风谱进行了曲线拟合。赵小平等(2016)基于强台风海鸥登陆期间的近地层测风资料,对台风过程不同位置的风廓线和湍流特征参数进行了细致的对比分析,但没有进一步考虑台风经过前后下垫面差异的影响。班欣等(2012)基于连云港沿海气象站的观测资料,分析了强冷空气、强对流和台风的湍流强度分布规律和特点,发现不同的大气层结稳定度和下垫面均会对湍流强度产生影响。

本文基于广东省珠江口的气象观测塔的测风数据,分析台风、强对流天气和强冷空气条件下,工程抗风最为关注的近地层的强风特性,以期为工程抗风设计提供参考。

1 观测塔与强风过程描述 1.1 沙田塔描述与观测设置

沙田塔是位于中国广东省东莞市的一个气象测风塔(22.853°N、113.574°E,塔基海拔高度8 m)。该塔位于江边,珠江入海口处,西面邻江,下垫面平坦,以耕地为主,周围地面长有稀疏杂草和低矮树木(图 1)。气象塔在10、20、40、60、80 m高度层安装有美国Renewable NRG Systems公司的NRG-#40c型杯式风速计和#200p型风向传感器,风速数据采集的时间步长为1 s,测量范围为1~96 m·s-1,仪器记录每10 min的平均风速、最大阵风(3 s)风速和1 s风速相对10 min平均风速的标准差,风向记录每10 min平均风向,仪器安装之前均在标准风洞中进行统一标定。在75 m高度层安装WindMaster Pro型三维超声风速仪,该仪器最大风速量程为65 m·s-1,采用10 Hz的采样频率,每0.1 s记录uvw三维方向的瞬时风速值、虚温和判别码。在观测期间,沙田塔测得了多次强风过程,包括台风莫拉菲、强对流强风和冷空气大风。

图 1 沙田塔地理位置图 Fig. 1 Geographic location of Shatian Tower
1.2 台风莫拉菲

台风莫拉菲于2009年7月19日00:50在深圳市大鹏半岛(南澳镇)沿海地区登陆,登陆时中心最低气压为965 hPa,中心附近最大风力为13级,达到38 m·s-1。“莫拉菲”在深圳市登陆后持续向西北偏西方向移动,强度逐渐减弱,先后穿过东莞、佛山、广州和肇庆等市。19日05:30左右,“莫拉菲”减弱为热带风暴,从沙田塔西南面约12 km处经过。图 2为台风路径图及观测塔距离台风中心的相对位置。在19日的00—04时,沙田塔测得的台风风向在321°~344°变化,为离岸风。随后风向随顺时针变化,直至06时开始,风向开始稳定在157°~175°,为向岸风。

图 2 台风莫拉菲路径图及观测塔距离台风中心的位置 Fig. 2 Path of Typhoon Molave and observation tower location relative to the typhoon center
1.3 强对流天气

本文选取了沙田塔测得的两次强对流天气的观测资料进行研究。

1.3.1 1号强对流过程

2009年8月23日,沙田塔观测到短时风速、风向剧烈变化的强对流飑线过程,该过程持续约1 h。图 3给出沙田塔东南面约3.3 km处自动气象站在23日13时至24日01时观测到的10 min平均风速、10 min平均风向和气温、气压时程曲线。由图 3可见,在19时前后有强对流过程影响,自动气象站19时测得10 min平均风速为9.5 m·s-1,比18时的10 min平均风速增加8.6 m·s-1,对应的10 min平均风向也发生显著变化;同时,气温急速降低,18时的气温为33.1℃,19时气温降至26.7℃;气压则从18时的1001.2 hPa增加到19时的1002.9 hPa;该站在19—20时还测得9.7 mm的累积雨量。此次强对流过程简称为1号强对流。

图 3 2009年8月23日1号强对流过程的(a)10 min平均风速和风向,(b)气温和气压 (虚线表示强对流经过时刻) Fig. 3 Variations of meteorological elements in convection process No.1 on 23 August 2009 (a) 10 min wind speed and wind direction, (b) temperature and air pressure (vertical dotted line: the time with severe convection)
1.3.2 2号强对流过程

2009年7月11日,沙田塔观测到持续约半小时的风速、风向剧烈变化的强对流飑线过程。由于该次强对流影响范围小,周边自动站未测得相关气象现象。下文将给出此次过程沙田塔的10 min平均风速、10 min平均风向并进行分析。此次强对流过程简称为2号强对流过程。

1.4 强冷空气

2008年12月21—23日一次强冷空气过程影响沙田现场测风区域,出现大风降温天气,此次强冷空气具有南下速度快、风力大的特点。图 4给出沙田塔东南面约3.3 km处的自动气象站在21日18时至23日00时观测到的10 min平均风速、10 min平均风向和气温、气压时程曲线。由图 4中可见,22日00时开始,10 min平均风速显著增加,风向转为偏北风,22日11时测得10 min平均风速最大值7 m·s-1;气温呈显著下降趋势,在22日00—11时,气温从20.1℃降为14.7℃,随后继续下降;气压则显著增加。

图 4图 3,但为2008年12月21—23日强冷空气过程 Fig. 4 Same as Fig. 4, but for severe cold air process from 21 to 23 December 2008
1.5 强风数据分析说明

台风、2号强对流和强冷空气过程均有杯式测风仪和超声测风仪的观测记录,1号强对流过程仅有杯式测风仪的观测记录,超声风速仪数据缺测。因此,文中风速、风向、风廓线、湍流强度的对比是基于四个强风过程的杯式测风仪数据进行的分析;而基于超声数据的风攻角、湍流空间积分尺度和湍流功率谱分析则针对台风、2号强对流和强冷空气过程进行。

2 平均风特性 2.1 风速、风向

图 5给出沙田塔80 m测得的台风、强对流和冷空气过程的强风10 min平均风速、10 min平均风向随时间的变化。从图 5a可见,台风莫拉菲过程10 min平均风速呈显著的“M”型双峰分布,双峰之间最低10 min平均风速为9.1 m·s-1,8级强风的10 min平均风向持续转变208°,根据Song et al(2016)的台风强风代表性判别标准可知,沙田塔测得的强风过程采集到了包括台风眼区风、眼壁和外围强风的完整台风风况数据。台风过程中80 m高度层测得最大10 min平均风速为25.5 m·s-1,3 s阵风风速为31 m·s-1图 5b5c显示受强对流影响, 10 min平均风速有显著的增加,1号强对流过程的80 m的10 min平均风速在18:30—19:00从5.5 m·s-1增加到25.3 m·s-1,对应3 s阵风风速为32.9 m·s-1,2号强对流过程的80 m的10 min平均风速在01:00—01:20从5 m·s-1增加到13.6 m·s-1,对应3 s阵风风速为18.8 m·s-1, 同时伴随着风向的显著突变;图 5d显示强冷空气过程中,80 m的10 min平均风速从01:20开始缓慢增大,在09:20测得最大的10 min平均风速(12.1 m·s-1),对应3 s阵风风速为14.6 m·s-1,风向以东北偏北风为主。台风、1号强对流、2号强对流和强冷空气过程风速的10 min升幅最大值分别是5.4、9.6、5.9和2.6 m·s-1。1号强对流过程的风速变化最为剧烈,风速的10 min升幅可分别达台风和强冷空气大风的1.8和3.7倍。

图 5 强风过程10 min平均风速和风向随时间的变化(a)台风过程,(b)1号强对流过程,(c)2号强对流过程,(d)冷空气过程 Fig. 5 Variations of 10 min wind speed and wind direction with time in strong wind processes (a) typhoon process, (b) severe convection process No.1, (c) severe convection process No.2, (d) severe cold air

计算各强风过程的10 min平均风向时间变率(指后一个10 min平均风向和前一个的差值),当差值为正时表示风向顺时针偏转,为负时表示风向逆时针偏转。结果显示,台风过程的最大10 min平均风向时间变率为28°·(10 min)-1,伴随着10 min平均风速从25.5 m·s-1降至21.0 m·s-1;强对流强风出现时风向显著突变,1号和2号强对流过程的最大10 min平均风向时间变率分别为109和120°·(10 min)-1,均在风速急剧回落后出现;强冷空气过程的最大10 min平均风向时间变率为14°·(10 min)-1

2.2 强风样本的选取

图 6给出了各个强风过程来风方位角的范围示意图。由图 6可见,台风过程的强风主要来自326°~24°和170°~174°方位角,1号强对流和2号强对流分别来自63°~73°和35°~53°方位角,强冷空气则来自19°~36°方位角。台风北面来风和强对流、冷空气的来风下垫面状况较一致,基本为平坦的城郊下垫面,以大片平坦耕地为主,而台风南面来风下垫面陆地面积较少,以江河水道下垫面为主。由于下垫面对风特性有显著的影响,为了减小不同强风来向下垫面的差异,在下文给出的测风参数比对中,台风强风样本只选取与强对流和冷空气来风下垫面状况较为一致的台风北面来风,且60 m高度的10 min平均风速达到9级及以上(≥20.8 m·s-1)的强风作对比分析,选取的时段为02:50—04:20。强对流的强风样本选自风速显著增大的时段,1号强对流的强风样本时段为18:50—19:10,2号强对流的强风样本时段为01:10—01:20;强冷空气过程的强风样本选取60 m高度的10 min平均风速增加到10 m·s-1并在附近波动的最大风时段,冷空气强风样本时段为07:40—14:00。GB 50009—2012(中华人民共和国住房和城乡建设部,2012)将下垫面类型分为4类,各类下垫面描述如下:A类指近海海面和海岛、海岸、湖岸及沙漠地区;B类指田野、乡村、丛林、丘陵以及房屋比较稀疏的乡镇;C类指有密集建筑群的城市市区;D类指有密集建筑群且房屋较高的城市市区。上述所选的沙田塔测得的强风样本的来风下垫面基本属于B类下垫面。

图 6 各强风过程来风方位角范围示意图 Fig. 6 Wind azimuth range of strong wind processes
2.3 风廓线

风廓线用以描述风随高度的变化情况,在此指风速随高度的变化形态。近地边界层中风速的垂直分布特征与大气层结、下垫面粗糙度和天气系统特性等有关。工程应用中常以幂指数或对数率来量化拟合描述风廓线,目前,国内外大部分相关规范(中华人民共和国住房和城乡建设部,2012; Hong Kong Building Department, 2004; Architectural Institute of Japan, 2004; American Society of Civil Engineers, 2006)推荐幂指数形式,表达式如下:

$ u = {u_1}{(\frac{\mathit{z}}{{{z_1}}})^\alpha } $ (1)

式中,uz高度处的风速,u1z1高度处的风速(单位:m·s-1),α为该高度层内的风速廓线指数。

图 7给出各强风过程的每10 min平均风廓线的变化图(其中由于冷空气时间较长,每20 min间隔给出10 min平均风廓线),横坐标表示该风廓线测得的时间。结合图 5可见,台风过程中03—04时左右的台风眼壁强风以及1号强对流过程的19时左右、2号强对流过程的01:20左右的强风经过沙田塔时,风廓线的形态均有显著变化,10 min平均风速随高度的增加非常显著,而冷空气强风经过时,风廓线的变化不明显。

图 7 各强风过程10 min风廓线随时间的变化(a)台风过程,(b)1号强对流过程,(c)2号强对流过程,(d)强冷空气过程 Fig. 7 Variations of 10 min mean wind speed profiles during strong wind processes (a) typhoon process, (b) severe convection process No.1, (c) severe convection process No.2, (d) severe cold air

根据式(1),计算各强风过程的平均风廓线(过程中的10 min风廓线的平均值)和最大的10 min风廓线的幂指数值。图 8给出各强风过程的最大的10 min风廓线及其幂指数拟合线。结果显示,强冷空气过程的风廓线幂指数值(α值)最小,平均风廓线和最大的10 min风廓线的α值分别为0.102和0.112;1号强对流过程的平均风廓线和最大的10 min风廓线的α值分别为0.121和0.122;2号强对流过程的平均风廓线和最大的10 min风廓线的α值分别为0.125和0.135;台风强风的平均风廓线和最大的10 min风廓线的α值分别为0.177和0.183。可见,在强冷空气和强对流过程中,风速随高度增加而增加的速度没有在台风条件下风速的增长速度快。台风强风平均风廓线的α值分别是强对流强风平均风廓线和冷空气大风平均风廓线的α值的1.4和1.7倍(表略)。GB 50009—2012(中华人民共和国住房和城乡建设部,2012)中给出A、B、C、D类的风廓线幂的α取值为0.12、0.15、0.22和0.30。可见,沙田塔B类下垫面下的台风强风廓线的α值要大于规范给出的0.15的参考值,而强对流的α值接近规范给出的A类下垫面的参考值0.12,而强冷空气的α值则更小。

图 8 强风过程的最大10 min风廓线及其幂指数拟合线 Fig. 8 Maximum 10 min wind profiles and fitting curves of the strong wind processes
3 脉动风特性 3.1 湍流强度

湍流强度反映了风的脉动强度,是确定结构脉动风荷载的关键参数,定义湍流强度为10 min时段的脉动风速标准差与水平平均风速的比值(盛裴轩等,2003)。本文的湍流强度(I)是基于杯式测风仪的1 s风速记录的标准差(σ,单位:m·s-1)相对10 min平均风速(U,单位:m·s-1)的比值计算所得结果。

$ I = \frac{\sigma }{U} $ (2)

图 9给出强风过程的80 m湍流强度与10 min平均风速随时间的变化。由图 9可见,在台风中心前后的强风区都有显著湍流强度增大的现象,尤其是04:00—04:20,台风强风叠加陆地的影响,湍流强度显著增大,台风中心经过后的台风强风区由于下垫面是江河水面为主,湍流强度增大幅度明显较小;强对流期间湍流强度在最大的10 min平均风速出现(19:00和01:20)前后都有显著的增大,主要是风速的骤然升降造成;而强冷空气的湍流强度变化不明显。台风、1号强对流、2号强对流和强冷空气期间测得的80 m高度层的最大的10 min平均风速分别为25.5、25.3、13.6和12.1 m·s-1,最大风时刻对应的80 m高度测得的湍流强度分别为0.094、0.146、0.199和0.091。表 1给出了各天气系统强风时段的湍流强度平均值,从表中可见,1号强对流天气下的强风湍流强度平均值分别是台风和强冷空气强风湍流强度平均值的1.4~2.1和1.7~2.0倍。在短时强对流过程中,风的组织结构可能较大尺度的强台风的涡旋结构更为紊乱,并且由于同样具有较大的风速,因而强对流对结构起着不容忽视的破坏作用。

图 9图 5,但为强风过程湍流强度与10 min平均风速 Fig. 9 Same as Fig. 5, but for turbulence intensity and 10 min wind speed

表 1 强风过程各高度层的平均湍流强度 Table 1 Average turbulence intensities of different heights in strong wind processes

基于三维超声观测数据,计算不同天气系统的强风在纵向、横向、竖向方面的湍流强度比值(Iu:Iv:Iw)。该比值有较大差异,75 m高度层台风强风的Iu:Iv:Iw=1:0.84:0.51,2号强对流过程大风的Iu:Iv:Iw=1:0.78:0.45,冷空气过程大风的Iu:Iv:Iw=1:0.93:0.61,与《公路桥梁抗风设计指南》(项海帆等,1996)中给出的1:0.88:0.5有差异。

3.2 风攻角

风攻角指风的来流方向与水平面的夹角。风攻角对结构物特别是柔性结构物的影响比较突出,其长期作用会加速或加重结构的疲劳损坏。强风时,风攻角的改变可能对结构造成突然损伤甚至破坏。风攻角主要由不均匀地形致使气流强迫抬升或下沉而产生,另外,不同的天气系统,尤其是涡旋结构的(如热带气旋、龙卷风等)强烈天气系统也可以导致风攻角的变化。在桥梁动力抗风设计中,应考虑高风速时风的平均攻角,根据《公路桥梁抗风设计指南》(项海帆等,1996),一般推荐采用±3°,平静海面和开阔平坦地形大约为0°。风攻角的计算公式如下:

$ {\rm{tan}}\theta = \frac{W}{U} $ (3)

式中,θ为风攻角(单位:°),W为垂直风风速平均值(单位:m·s-1),U为主风向风速(单位:m·s-1)。

图 10给出强风过程的10 min平均风速和风攻角随时间的变化。由图 10可见,由于台风中心经过前后来风从离岸风转为向岸风,离岸强风攻角在-1.4°~0.1°变化,而受地形抬升影响,台风中心经过后风攻角显著增加。但仍在-3°~3°范围内。冷空气大风的风攻角在0°附近变动。而2号强对流影响前,风攻角有显著增加的现象,测得最大风攻角为4°,但此时的10 min平均风速仅为4.3 m·s-1,风速急速增大时,风攻角仍在-3°~3°以内。

图 10 强风过程75 m高度层风攻角与10 min平均风速随时间的变化(a)台风过程,(b)2号强对流过程,(c)强冷空气过程 Fig. 10 Variations of attack angle and 10 min mean wind speed at the height of 75 m with time during strong wind processes (a) typhoon process, (b) severe convection process No.2, (c) severe cold air
3.3 湍流空间积分尺度

湍流空间积分尺度是表征湍涡大小的物理量。由于结构风荷载对湍流尺度特性的敏感性,湍流积分长度通常是一项重要的风特性指标。湍流空间积分尺度的大小决定了脉动风对结构的影响范围。与纵向脉动速度有关的涡旋的纵向、横向和垂直方向的平均尺寸用LuxLuyLuz表示,在数学上可定义Lux为:

$ L_u^x{\rm{ = }}\frac{1}{{\sigma _u^2}}\int_0^\infty {{R_{{u_1}{u_2}}}\left(x \right){\rm{d}}x} $ (4)

式中,Ru1u2(x)是两个不同空间位置上的纵向脉动速度u1=u(x1, y1, z1, t)和u2=u(x1+x, y1, z1, t)的互协方差函数,t为时间,σu2为脉动速度u的方差,同样的定义也适用于LuyLuz。假设气流以U的速度迁移,则脉动速度u(x1, t+τ)可以定义为u(x1, -x, τ),x=Ut,即Taylor假设,那么式(4)可以改写为:

$ L_u^x = \frac{U}{{\sigma _u^2}}\int_0^\infty {{R_u}\left(\tau \right){\rm{d}}\tau } $ (5)

式中,Ru(τ)是脉动量u(x1, t)的自相关函数,Ru(τ)=E[u(t)u(t+τ)],一般积分上限取到Ru(τ)=0.05σu2。基于测风塔上的采样频率为10 Hz的三维超声测风数据,计算纵向、横向和垂直方向脉动风速的平均空间尺度,下文分别用LuLvLw表示(单位:m)。

图 11给出强风过程的10 min平均风速和湍流空间积分尺度随时间的变化。当台风强风经过时,水平(纵向和横向)湍流空间积分尺度明显增大了一个量级,由几十米变为几百米,但垂直方向没有明显变化,台风强风的平均三维湍流空间积分尺度分别为318、111、37 m。2号强对流过程在风速显著增加时,三维方向的湍流空间积分尺度均有增大的现象,平均三维湍流空间积分尺度分别为238、108、68 m。强冷空气大风经过观测塔时,三维方向的湍流空间积分尺度也有增大的现象,平均三维湍流空间积分尺度分别为115、81、40 m。

图 11图 10,但为湍流空间积分尺度与10 min平均风速 Fig. 11 Same as Fig. 10, but for turbulence integral length scale and 10 min mean wind speed
3.4 湍流功率谱

大气中某个空间点的脉动风速是由许多不同尺度的涡在该点形成的各种频率的脉动叠加组合而成的。所以,湍流的脉动动能可以被认为是各种频率的涡的贡献。脉动功率谱是湍流脉动能在各个频率上的分布,用来描述湍流中不同尺度的涡的湍流脉动动能的贡献。湍流功率谱函数Si(i=u, v, w)能够准确地描述脉动风的特性,它们在频域上的全积分等于脉动风对应方向上的湍流动能,即

$ \int_0^\infty {{S_i}\left(n \right){\rm{d}}n} = \sigma _i^2 $ (6)

式中,i=u, v, w(单位:m·s-1);n为频率(单位:Hz),Si在频域上的分布可以描述湍流动能在不同尺度水平上的比例(单位:m2·s-2)。功率谱函数采用傅里叶变换计算得到。以随机信号x(t)为例,取有限时段子样本xk(t)(0≤tT),定义

$ {S_x}\left({f, T, k} \right) = \frac{1}{T}X_k^*\left({f, T} \right){X_k}\left({f, T} \right) $ (7)

式中,${X_k}\left({f, T} \right) = \int_0^T {{x_k}\left(t \right){{\rm{e}}^{ - i2{\rm{ \mathsf{ π} }}ft}}{\rm{d}}t}, {X_k}(f, T)$表示xk(t)的有限傅里叶变换,Xk*(f, T)为Xk(f, T)的共轭,T为有限长度时间。

图 12给出各强风时刻10 min样本的湍流功率谱图。拟合1~4 Hz频率的区间,分析惯性子区(图中拟合区域)的湍流谱特征,可以发现,台风强风的湍流功率谱样本不满足在惯性子区湍流能谱密度的-5/3次方与频率成正比的规律和各向同性假设,而2号强对流和强冷空气大风样本则基本符合-5/3律和各向同性假设。表 2给出强风样本在3个频段内三维方向的平均湍流能谱密度。由表 2可见,台风强风的湍流能谱值要显著高于强对流和强冷空气过程,其中,在桥梁敏感频域台风的平均湍流能谱值约为强对流和强冷空气的3倍。

图 12 强风75 m高度层过程湍流功率谱图 (a)台风,(b)2号强对流,(c)强冷空气[b(u)、b(v)、b(w)分别表示纵向、横向和竖向湍流功率谱在惯性子区的线性拟合斜率值] Fig. 12 Spectra of strong wind processes at the height of 75 m (a) typhoon process, (b) severe convection process No.2, (c) severe cold air [b(u), b(v), b(w) separately indicate the slopes of linear fittings of longitudinal, lateral and vertical spectra in inertial sub-area]

表 2 强风过程75 m高度层三维方向平均湍流能谱值(单位:m2·s-1) Table 2 Average power spectral density at the height of 75 m during strong wind processes (unit: m2·s-1)
4 结论

基于广东省东莞市的沙田气象塔在台风莫拉菲、强对流和强冷空气期间的观测数据,分析不同强风天气系统的近地边界层的平均风和脉动风特性参数的差异,主要表现如下:

(1) 台风、1号强对流、2号强对流和强冷空气过程中风速的10 min升幅最大值分别是5.4、9.6、5.9和2.6 m·s-1。1号强对流的风速变化最为剧烈,风速的10 min升幅分别可达台风和冷空气大风的1.8和3.7倍。台风过程的最大10 min平均风向时间变率为28°·(10 min)-1;1号和2号强对流过程的最大10 min平均风向时间变率分别为109和120°·(10 min)-1;强冷空气过程的最大10 min平均风向时间变率为14°·(10 min)-1

(2) 台风强风平均风廓线幂指数值分别是强对流强风平均风廓线和强冷空气大风平均风廓线的幂指数值的1.4和1.7倍。沙田塔测得的台风强风廓线的幂指数值要大于规范给出的B类下垫面风廓线0.15的幂指数参考值,而强对流的幂指数值接近规范给出的A类下垫面的参考值0.12,强冷空气的幂指数值则更小。

(3) 台风、1号强对流、2号强对流和强冷空气期间测得的80 m高度层的最大的10 min平均风速对应的湍流强度分别为0.094、0.146、0.199和0.091。综合10~80 m高度的观测值,1号强对流天气下的强风湍流强度平均值分别是台风和强冷空气强风湍流强度平均值的1.4~2.1和1.7~2.0倍。

(4) 台风和强对流的风攻角基本在公路桥梁抗风设计指南建议的-3°~3°范围内。强冷空气大风的风攻角变幅较小,在0°附近变动。

(5) 台风、强对流和强冷空气大风影响期间湍流空间积分尺度均有增大的现象,尤其是当台风强风经过时,水平(纵向和横向)湍流空间积分尺度明显增大了一个量级。

(6) 台风强风的湍流功率谱样本不满足在惯性子区湍流能谱密度的-5/3次方与频率成正比的规律和各向同性假设。台风强风的湍流能谱值要显著高于强对流和强冷空气,其中,在桥梁敏感频域(0.1~0.5 Hz)台风的平均湍流能谱值约为强对流和强冷空气的3倍。

参考文献
班欣, 冯还岭, 祁欣, 等, 2012. 连云港沿海近地层湍流强度特征[J]. 气象科技, 40(2): 285-292.
Ban X, Feng H L, Qi X, et al, 2012. Analysis of surface layer turbulence intensity along Lianyun-gang coastland[J]. Meteor Sci Technol, 40(2): 285-292 (in Chinese). DOI:10.3969/j.issn.1671-6345.2012.02.025
郭冬艳, 姜涛, 陈有龙, 等, 2011. 琼州海峡沿岸大风分布规律及影响系统分析[J]. 气象, 37(11): 1372-1379.
Guo D Y, Jiang T, Chen Y L, et al, 2011. Analysis on distribution pattern of high wind events in Qiongzhou Strait coastal zones and its westher systems[J]. Meteor Mon, 37(11): 1372-1379 (in Chinese). DOI:10.7519/j.issn.1000-0526.2011.11.007
黄林宏, 宋丽莉, 李刚, 等, 2016. IEC的风机选型参数估算方法在中国的适用性研究[J]. 气象, 42(12): 1522-1530.
Huang L H, Song L L, Li G, et al, 2016. Applicability of estimation methods of IEC WTGS selecting parameters in China[J]. Meteor Mon, 42(12): 1522-1530 (in Chinese). DOI:10.7519/j.issn.1000-0526.2016.12.010
廖菲, 邓华, 李旭, 2017. 基于风廓线雷达的广东登陆台风边界层高度特征研究[J]. 大气科学, 41(5): 949-959.
Liao F, Deng H, Li X, 2017. A study on boundary layer height characteristics of lan-ding typhoons by wind profilers in Guangdong Province[J]. Chinese J Atmos Sci, 41(5): 949-959 (in Chinese).
庞加斌, 林志兴, 葛耀君, 2002. 浦东地区近地强风特性观测研究[J]. 流体力学实验与测量, 16(3): 32-39.
Pang J B, Lin Z X, Ge Y J, 2002. Field measurements of strong wind characteristics near ground in Pudong District[J]. Experim Measurem Fluid Mech, 16(3): 32-39 (in Chinese). DOI:10.3969/j.issn.1672-9897.2002.03.006
钱燕珍, 高拴柱, 黄思源, 等, 2013. 强台风海葵登陆前后强度变化的观测分析[J]. 气象, 39(10): 1265-1274.
Qian Y Z, Gao S Z, Huang S Y, et al, 2013. Analysis on intensity variation of Haikui before and after its landing[J]. Meteor Mon, 39(10): 1265-1274 (in Chinese). DOI:10.7519/j.issn.1000-0526.2013.10.004
盛裴轩, 毛节泰, 李建国, 等, 2003. 大气物理学[M]. 北京: 北京大学出版社: 225.
Sheng P X, Mao J T, Li J G, et al, 2003. Atmospheric Pysics[M]. Beijing: Peking University Press: 225 (in Chinese).
宋丽莉, 毛慧琴, 黄浩辉, 等, 2005. 登陆台风近地层湍流特征观测分析[J]. 气象学报, 63(6): 915-921.
Song L L, Mao H Q, Huang H H, et al, 2005. Analysis on boundary layer turbulent features of landfalling typhoon[J]. Acta Meteor Sinica, 63(6): 915-921 (in Chinese). DOI:10.3321/j.issn:0577-6619.2005.06.008
宋丽莉, 毛慧琴, 钱光明, 等, 2006. 热带气旋对风力发电的影响分析[J]. 太阳能学报, 27(9): 961-965.
Song L L, Mao H Q, Qian G M, et al, 2006. Analysis on the wind power by tropical cyclone[J]. Acta Energiae Solaris Sin, 27(9): 961-965 (in Chinese). DOI:10.3321/j.issn:0254-0096.2006.09.020
唐明晖, 姚秀萍, 王强, 等, 2017. 不同强度台风相伴随的内陆台前飑线对比分析[J]. 气象, 43(8): 912-923.
Tang M H, Yao X P, Wang Q, et al, 2017. Contrastive analysis of inland Pre-TC squall line accompanied by typhoons with different intensities[J]. Meteor Mon, 43(8): 912-923 (in Chinese).
项海帆, 林志兴, 鲍卫刚, 等, 1996. 公路桥梁抗风设计指南[M]. 北京: 人民交通出版社: 15.
Xiang H F, Lin Z X, Bao W G, et al, 1996. Guidelines for Wind Resistance Design of Highway and Bridges[M]. Beijing: China Communications Press: 15 (in Chinese).
赵小平, 朱晶晶, 樊晶, 等, 2016. 强台风海鸥登陆期间近地层风特性分析[J]. 气象, 42(4): 415-423.
Zhao X P, Zhu J J, Fan J, et al, 2016. Analysis on wind characteristics in surface layer during landfall of Typhoon Kalmaegi[J]. Meteor Mon, 42(4): 415-423 (in Chinese). DOI:10.3969/j.issn.1000-6362.2016.04.005
中华人民共和国住房和城乡建设部, 2012.建筑结构荷载规范: GB 50009—2012[S].北京: 中国建筑工业出版社: 220-221.
Ministry of Housing and Urban-Rural Development of the People's Republic of China, 2012.GB 50009-2012 Load Code for the Design of Building Structures[S].Beijing: China Architecture & Building Press: 220-221(in Chinese).
周福, 蒋璐璐, 涂小萍, 等, 2017. 浙江省几种灾害性大风近地面阵风系数特征[J]. 应用气象学报, 28(1): 119-128.
Zhou F, Jiang L L, Tu X P, et al, 2017. Near-surface gust factor characteristics in several disastrous winds over Zhejiang Porvince[J]. J Appl Meteor Sci, 28(1): 119-128 (in Chinese).
Architectural Institute of Japan, 2004.Recommendations for Loads on Buildings: AIJ-RLB-2004[S].Tokyo: Architectural Institute of Japan: 6-10-6-11.
American Society of Civil Engineers, 2006.Minimum Design Loads for Buildings and Other Structures: ASCE/SEI 7-05[S].Reston, VA: ASCE: 251.
Hong Kong Building Department, 2004.Code of Practice on Wind Effects in Hong Kong[S].Hong Kong: Building Department: 19.
Cao S Y, Tamura Y, Kikuchi N, et al, 2009. Wind characteristics of a strong typhoon[J]. J Wind Eng Ind Aerod, 97(1): 11-21. DOI:10.1016/j.jweia.2008.10.002
Song L L, Chen W C, Wang B L, et al, 2016. Characteristics of wind profiles in the landfalling typhoon boundary layer[J]. J Wind Eng Ind Aerod, 149: 77-88. DOI:10.1016/j.jweia.2015.11.008
Song L L, Li Q S, Chen W C, et al, 2012. Wind characteristics of a strong typhoon in marine surface boundary layer[J]. Wind Struct, 15(1): 1-15.
Song L L, Pang J B, Jiang C L, et al, 2010. Field measurement and analysis of turbulence coherence for Typhoon Nuri at Macao Friendship Bridge[J]. Sci China Technol Sci, 53(10): 2647-2657. DOI:10.1007/s11431-010-4084-3
Wang B L, Hu F, Cheng X L, 2011. Wind gust and turbulence statistics of typhoons in South China[J]. Acta Meteor Sinica, 25(1): 113-127. DOI:10.1007/s13351-011-0009-8