深圳位于亚热带季风气候区，因为地处广东省南部海滨地区，海洋性气候显著。季风、强对流以及台风等恶劣天气带来的大风常常引起港口集装箱倾翻，影响船舶靠岸，引发港口高架作业安全事故，严重干扰了港湾码头的正常业务以及威胁到作业人员人身安全，给区域安全生产带来重大影响。强风是港口码头安全作业、轮船安全行驶的最敏感的致灾因素(顾国贤，2007；郑永光等，2016a；Li et al, 2016；王延青等，2017)，加强不同季节条件下的风力预测和风险评估，为强天气影响下大风的防范提供指引，是防风减灾的必要措施，更是港口地区安全运行的重要保障。

在港口码头的防风管理规范中的强风多指阵风(顾国贤，2007；Bos，2017)。本文中平均风定义是10 min内的平均，阵风定义是3 s内的平均风(Bureau of Meteorology, 2016；Li et al, 2016)。林良勋等(2006)指出，热带气旋中的大风具有明显的阵性特点，不少热带气旋会以较大的阵风造成破坏；澳大利亚气象局研究表明：空旷陆地上，一般来讲，阵风会比平均风大40%，在海洋上阵风会比平均风大20%~30%；对建筑物等财产造成最大破坏的往往是因为更强的阵风所引起(Bureau of Meteorology, 2016)。对于港口码头来讲，大型港机设备结构高大，迎风面积大，很容易遭受台风或突发性阵风的袭击，造成巨大损失(顾国贤，2007)，因此需要对强风严格防范；新的欧洲码头起重机标准(EN3001)明确指明港口码头的起重机设计需按照阵风值来确定开展(Bos, 2017)。

台风、冷空气及强对流天气都能给受灾地区带来大风(林良勋等，2006；郑永光等，2016a; 2016b; 2017；Li et al, 2016；李佰平等，2016；唐文苑等，2017；王延青等，2017；杨新林等，2017)。深圳地区一年四季都会受强风影响，冬春之际冷空气频繁；3—8月多强对流天气，出现雷雨大风；夏秋季节则常受登陆台风的影响。以往港湾码头在强天气影响下(特别是台风影响)的风场研究以定性结论居多，定量的分析较少。在港口风场预测和大风风险评估方面，中国尚无有效的客观方法可供业务应用参考，因此及时准确地做好强天气条件下的港口风力预报预警是一项极需要加强研究的工作。本研究将分析深圳地区8个沿海(本文指台站距离海岸线1000 m以内)及内陆(指距离海岸线1000 m以上)气象自动站的历史小时数据，分析各站在不同季节，即春、夏、秋及冬季条件下，阵风大小及方位特征；计算不同季节条件下各气象站的阵风系数，并分析不同下垫面对阵风系数的影响及特征。最后，本研究将对典型的台风、冷空气及强对流天气条件下，不同下垫面台站的阵风及平均风记录进行个例比较分析，以厘清深圳地区不同下垫面台站在不同季节条件受强风影响下的风力特征。本研究工作对于沿海地区的防风减灾工作可以提供有价值的参考。

本论文研究的区域主要为深圳的沿海地区，包括6个沿海自动站：背仔角、盐田港、东角头、蛇口码头、左炮台和妈湾港；为与深圳其他地区进行对比分析，本研究还选取分析深圳国家基本站竹子林站和龙岗基地自动气象站。8个自动站的详细信息及地理位置见表 1及图 1。本文选取的数据为深圳地区自动气象站的小时观测数据，由深圳市气象局提供。自动气象站数据的更新频率为每分钟，即分钟数据。该分钟数据中包含即时3 s平均风，1、2和10 min平均风，以及过去1 h内最大10 min平均风和过去1 h内最大3 s平均风。历史小时数据是在自动站每分钟输出的数据基础上经过计算处理而得。本研究使用的数据为小时内最大10 min平均风速(小时最大风速)和小时内最大3 s平均风速(小时极大风速)。

表 1显示了各气象台站开始业务运行的时间，因各站数据长度不同，本研究选取2007年6月1日至2017年6月1日共11年数据进行分析计算。

根据蒲福风级标准，3级风以下，无论是对海洋的船舶运行，还是对陆地的人群活动，影响都很小(Huler, 2004)，因此本文主要是针对小时内最大10 min平均风速3级以上的情形进行分析，即小时内最大10 min平均风速≥5 m·s-1的情况，也就是说如果小时最大平均风速 < 5 m·s-1，则该小时的数据不在季节阵风特征及季节阵风系数分析数据库中；但在极端天气大风个例比较分析部分，则采用了完整的逐小时观测数据，没有进行任何筛选。

由于大陆和海洋在一年之中增温和冷却程度不同，风向会随季节的变化而有规律的改变(Ramage, 1971; 汪品先，2009)，本文将按不同季节对各站进行风场的特征研究：春季为每年的3—5月，夏季为6—8月，秋季为9—11月，冬季为12月至次年2月。

本文将利用风玫瑰图对各站的阵风大小及方位进行分析描述。风玫瑰图也叫风向频率玫瑰图，它根据某一地区多年统计的各个风向和风速的值进行绘制，玫瑰图上所表示的风向(即风的来向)，是指从外面吹向地区中心的方向。按照笛卡尔坐标系，0°为正北风，90°为正东风，等。玫瑰叶瓣长度为阵风的大小。风玫瑰图长期以来在气象服务信息系统中得到广泛应用(高加嵘, 1977; 郁云, 2010; Oktal and Yildirim, 2016)，对于涉及城市规划、环保、风力发电、机场跑道建设等领域有着重要的意义。

阵风系数被定义为瞬时风与持续风力的比值(Tyner et al, 2015)：

式中，GF为阵风系数；G为阵风风速(单位：m·s-1)，即3 s平均风速值；U为10 min平均风速值。本文计算小时内阵风系数，G取小时内最大阵风风速，即该小时内出现的最大3 s平均风速值；U为小时内出现的最大10 min平均风速值(Shu et al, 2015; 魏晓琳等，2016)。阵风系数对于风电场的结构设计、空气污染控制、空气流通评估以及前景预计至关重要(Shu et al, 2015)。不同气象站因为下垫面条件的不同，阵风系数会有不同。探讨沿海及内陆不同气象台站的阵风系数特征，可以帮助了解不同粗糙度的下垫面条件对阵风系数的影响(Král, 2007;陈雯超等，2011; 王巍竹，2014)。研究沿海地区阵风的变化特征对于保障港口、码头行业的作业安全可以发挥重要指导作用(魏晓琳等, 2016；周福等，2017)。

按照第1.2节中介绍的方法，各站历史数据中先剔除10 min平均风速 < 5 m·s-1的样本，然后各站按四季进行阵风的特征分析。图 2列出背仔角(BZJ)、妈湾港(MWG)、蛇口码头(SKMT)、盐田港(YTG)、左炮台(ZPT)以及竹子林(ZZL)、龙岗基地(LGJD)的四季阵风风向以及风速大小的风玫瑰图，由于图篇幅原因，且东角头站(DJT)四季的阵风较其他海滨气象站的阵风要小，风力危害不大，所以图 2中没有列出。从该图中可以看出BZJ与YTG作为东部沿海自动站，相邻距离不远，但是四季的阵风具有不同的特点：YTG在春、夏、秋季的三个季节里的大风，东、南、西、北四个方向都会有，但是在冬季，西南方向的大风明显少于其他象限的大风频率及速度。BZJ的阵风则是所有的季节都发生在第一、二和四象限，第三象限鲜有大风发生；这有可能是BZJ的西南方向有地形或建筑物阻挡，对西南向来的风有缓冲，所以受西南方大风的影响很小；BZJ的春、冬季的大风多来自于北风，夏季大风多来自于东南风，秋季大风则多来自于北风及东南风。

与东部海滨气象站相比，位于西部海滨的MWG、SKMT、ZPT气象站的大风无论是频率还是风速都小于东部的海滨气象站，这是因为西部海滨各站的东南方向有香港岛存在(Li et al, 2016; 王延青等，2017)，因而缓减了西向移动台风对西部各站的风力影响。而YTG与BZJ的南面及东南面为海洋，极易遭受西向移动及在广东省珠江口附近登陆台风的大风影响。西部的几个滨海台站中，除了夏季，其他三个季节中大于28.5 m·s-1的大风频率都非常小。7个气象站中，相比沿海海滨站，位置偏内陆的ZZL与LGJD站，大风的频次普遍少于沿海的其他站。

为了更进一步地把多年间各气象站在不同季节因为不同的天气过程所引发的最大三个阵风风速筛查出来，设定这三个最大阵风值间隔必须为24 h以上，即为受不同的天气系统影响。以YTG夏季为例，经计算，多年来最大三个阵风皆为2012年7月24日因台风韦森特引发，此天气过程只保留最大的阵风观测，即2012年7月24日02时的阵风值35.7 m·s-1，排名第二和第三的24日04和03时的阵风值因为与排第一的阵风值属于相同的天气过程，则进行剔除。YTG由不同天气过程影响的排名第二的夏季阵风为2008年8月22日21时的32.7 m·s-1；排名第三的阵风值为2009年7月19日02时的32.3 m·s-1。以这种计算方法，各站不同季节的最大三个阵风值描绘在图 3中。把不同季节各气象站排名第二的阵风用直线连接，绿色三角线、红点线、黄色矩形线、蓝色菱形线分别代表春、夏、秋、冬季的最大阵风观测。各台站最大阵风值点上绘制的短细实线代表该阵风观测的方位，正北方为0°方向，正东方为90°方向。

从图 3可以看出，各台站在夏、秋季受台风影响的最大阵风普遍大于冬春之际受冷空气影响的最大阵风。8个气象站中，BZJ站四季最大阵风较其他气象站风力普遍要大，YTG夏季与秋季受台风影响时的阵风大小与BZJ相当，历史最大阵风都有12级以上；但是在冬春之际，YTG的风力影响明显小于BZJ。BZJ冬春季节最大的三个阵风都为北风，是冬春之际冷空气影响的典型方向特征；BZJ、YTG夏秋之际的最大阵风，有北风也有南风，这是与台风不同的位置以及气象站本身的地理环境有关。台风登陆前，台风低压区在海洋上；台风登陆后，低压区会位于各气象站的北部，因此台风引发的气象站的风力方向会发生变化。比较沿海气象站，内陆站LGJD、ZZL四季的风影响稍小(与图 2结论一致)。ZZL虽然没有邻海，但距离海岸线不远(图 1)，因此在夏秋之际，也可能受台风影响，出现大风。

按照第1.2节中介绍的方法计算各气象站不同季节的阵风系数，各站四季的阵风系数见图 4。进一步，计算各站不同季节阵风系数的90%，50%，10%分位数(图 5)。

在医学、环境及气象领域，常常计算数据样本的10%，50%，90%百分位数来代表样本的最小值、中间值及最大值的统计特征(钱幼琼等，1988；Curry et al, 2016; Roshan et al, 2017)。本研究也计算各气象站四季阵风系数的10%，50%，90%百分位数(图 4，图 5)可以看出，一般情况下，各站中YTG的阵风系数最小，其阵风系数的中位数值在春、夏、秋季都小于1.5，夏季最小。其原因是，深圳地区的夏季盛行南风，大风均来自于海上，作为港口码头站，YTG的南方，尤其东南方直接是广阔的海域，下垫面摩擦系数小，因此该站在台风影响下，其大风会较其他站更具持续性，即其平均风也会很大。夏季大风过程，按各台站阵风系数50%百分位数大小顺序：YTG < BZJ < MWG < SKMT < ZZL < DJT=LGJD=ZPT，各台站阵风系数10%百分位数、90%百分位数的顺序也表现出与50%百分位数相似的特征，说明夏季各台站在台风影响下的大风情况稳定；秋季各台站阵风系数10%，50%，90%百分位数，仍然是YTG最小，但其他台站的阵风系数顺序则表现不稳定。春季，YTG的阵风系数各百分位数仍然是各台站中最小，说明该站在春天冷空气影响下，风力影响稳定。LGJD在各季节中都表现出大的阵风系数，说明该站在任何季节大风影响都非常具有瞬时性，即阵风与平均风的比值较其他近海台站的要高。从图 5可以看出最邻海的三个港口码头站：YTG、MWG、SKMT的阵风系数较其他站的阵风系数要小；BZJ站坐落于深圳东部，其南方为广阔的海洋，但是其北面邻近梧桐山，因此地面平均粗糙度会加大，这应该是BZJ比YTG阵风系数要大的主要原因。总之，邻海的台站(特别是港口码头站)，因离海近，平均粗糙度较距海稍远的台站要小，因此阵风系数小，距离海岸线远的台站，因陆地粗糙度增加，因此阵风系数变大，这种空间分布特征与前人研究结论基本一致(Král, 2007;陈雯超等，2011; 王巍竹，2014；Bureau of Meteorology, 2016)。

从前面计算得到各站阵风中，BZJ与YTG的阵风最大，通过阵风系数计算，得到邻海最近的YTG的阵风系数最小。进一步比较BZJ与YTG历史各季节不同天气过程的最大三个平均风值(图 6)。可以看出，历史上两个台站的最大平均风均发生在夏季，其次为秋季，即因台风引发的最大平均风大于冬、春季因季风或强对流天气引发的最大平均风。YTG因台风原因造成的最大平均风一般来讲大于BZJ的最大平均风；而冬、春季，BZJ的最大平均风则大于YTG的最大平均风。BZJ与YTG夏秋之际因台风影响的最大平均风分别达到9级和10级。

从前文分析中可见，各站在夏、秋季台风条件下，YTG具有最大的阵风及平均风，而在冬、春季，BZJ具有最大的阵风和平均风。本节分析在台风、冷空气及强对流天气极端个例条件下的YTG和BZJ的阵风、平均风时间序列，并选择ZZL作为内陆代表站进行比较，得到不同下垫面环境条件下气象站的阵风、平均风及阵风系数特征。

YTG最大阵风发生在2012年7月24日02时，为35.7 m·s-1，因台风韦森特引发。图 7a显示在台风韦森特登陆前后(23日00时至26日00时)三站各小时的最大阵风及最大平均风时间序列，以及相应的阵风系数值。从图中可以清晰看出受韦森特影响，沿海两站的平均风及阵风都大于内陆的ZZL；YTG最大阵风略大于BZJ(35.7 m·s-1比34.4 m·s-1)，但是YTG的最大平均风(24.9 m·s-1)则明显大于BZJ的最大平均风值(19.7 m·s-1)。从图 7a看出，23—24日，YTG的阵风系数普遍小于2，其风力达到最大时的阵风系数为1.51；相比较，BZJ 23—24日，阵风系数在1.5~4.0，阵风最大(34.4 m·s-1)时刻的阵风系数为1.75。说明YTG在台风韦森特影响下，风更具持续性，而BZJ的风力阵性特点较大。稍内陆的ZZL则具有比YTG、BZJ要小的阵风及平均风，但因其近内陆，下垫面较沿海粗糙，其阵风系数普遍较YTG及BZJ大。另外，从图中可以看出，各站最大阵风系数发生时间对应的都是平均风观测很小的时刻，比如，YTG阵风系数最大时刻发生在台风登陆后并远离深圳的25日07时，此时对应的YTG最大小时平均风为1.8 m·s-1，最大小时阵风为7.9 m·s-1。

BZJ从有气象记录以来冬春之季的最大阵风发生在2010年12月16日06时，为36 m·s-1，因北方强冷空气南下引发。冷空气影响前后(15日00时至18日00时)BZJ、YTG及ZZL的小时最大阵风及最大平均风时间序列见图 7b。可以看出，在冷空气影响下，沿海站比内陆具有更大的风影响，BZJ在冷空气影响下具有比YTG更大的阵风及平均风。相比较图 7a与7b，在冷空气影响下，各站的阵风系数要大于在台风影响下的阵风系数，说明受冷空气影响，各站的风力比台风影响下更具阵性特点。

2013年5月21日前后的强雷雨是深圳近年来最强的过程之一。20—22日，受到切变线影响，深圳全市普遍出现暴雨降水过程，并伴有雷暴，雷雨时普遍记录到8~9级短时阵风。图 7c显示了强天气发生期间(20日00时至23日00时)BZJ、YTG和ZZL的平均风、阵风和阵风系数的变化情况。从图中可以看到，本次强天气过程平均风序列和阵风序列存在两个峰值，分别发生在20日17—20时和21日12时。强天气不仅在沿海站点引发强风(YTG最大阵风为22.1 m·s-1；BZJ最大阵风为21.8 m·s-1)，在城市稍内陆地区的ZZL，也会引发同等强度的大风(22.8 m·s-1)，对城市居民区具有很大的潜在威胁。另一方面，对比三个过程的阵风系数，强对流天气条件下的阵风系数远大于台风和冷空气条件下的阵风系数，说明强对流天气影响下的大风最具阵性特点。

更进一步绘制了三个不同天气过程中全深圳市各站的最大阵风空间分布(图 8)。从图中可以看出，三种天气条件下，深圳地区最大风影响：台风>冷空气>强对流。虽然风力分布不是完全相同，但是在台风和冷空气条件下，最大阵风的发生地点都是在东部沿海和高地地区；在强对流天气条件下，最大阵风的发生地点则在中部地区零散分布，非常具有随机性，不一定是发生在东部沿海和高地区域，因此，强对流天气下的风力定时定点精细化预报难度大。

本研究基于深圳地区不同地域分布的BZJ、DJT、MWG、SKMT、YTG、ZPT、LGJD、ZZL 8个自动气象站2007年6月1日至2017年6月1日11年的小时观测数据，分析讨论了不同季节条件下各站的阵风大小及方位特征；并计算不同季节条件下各气象站的阵风系数，分析不同下垫面对阵风系数的影响。研究结果发现，8个气象台站在夏、秋季受台风影响的最大阵风普遍大于冬春之际受冷空气影响的最大阵风；各站中，位于深圳东部沿海的BZJ与YTG是各站中大风频次最多的台站。BZJ各季节中的阵风比其他站普遍要大，其次为YTG；YTG夏、秋季受台风影响时的阵风大小与BZJ相当，历史最大阵风都达到12级以上；但是在冬春之际，YTG的风力影响明显小于BZJ。相比沿海气象站，内陆站LGJD、ZZL四季的风影响稍小。

另一方面，通过阵风系数的计算表明，各台站中，距离海岸最近的港口码头气象站，YTG、MWG、SKMT具有最小的阵风系数；随着离海岸距离的增加，单位面积下垫面的摩擦系数也会随之增加，阵风系数变大。比较阵风大风频率高的YTG与BZJ，BZJ的阵风系数要大于YTG，即BZJ的大风比YTG更具阵性的特点；夏秋之际，受台风影响下，YTG的平均风会稍大于BZJ的平均风影响。距离海岸线稍远的ZZL与LGJD气象站，则具有比沿海站更大的阵风系数，即其大风比沿海站更具阵性的特点。

最后的极端天气个例分析结果表明，相比较台风、冷空气条件下的极大风，强对流天气引发大风最具有阵性的特征，大风发生地点更随机，虽然强对流天气引发的大风量级不及台风和冷空气，但灾害发生定点定时精细化预报难度大，需要额外关注。