1. 引言

自古以来，人与自然有着密不可分的联系。龙卷风作为一种自然界破坏力极强的强对流天气之一，对人民的生命财产安全具有不可忽视的影响。例如，2005年7月30日11时50分左右，安徽省灵璧县出现龙卷风，此次龙卷风共造成15人死亡、51人受伤，房屋倒塌1617间，农作物受灾面积达27万亩；2023年6月1日14时10分左右，受强对流天气影响，辽宁省开原市金沟子镇出现龙卷风，导致80余人受灾，60余处的房屋造成了不同程度的损坏，农业受损面积约4000亩。由此可见，龙卷风造成的危害极为严重，加强对龙卷风的认识并建立有效的防灾减灾措施迫在眉睫。

目前，国内外对于龙卷风的研究主要集中在龙卷风的灾情调查、形成机制、风场特性及龙卷风作用下建筑物的风荷载与风致响应四个方面。自20世纪70年代起，便有学者开始对龙卷风发生后的灾情情况进行记录分析。例如，相德勤等[1]对1978年4月14日发生在陕西干县周城公社的龙卷风的伤亡人数及财产损失进行了统计分析；张军[2]对1981年5月发生在河北沫水的一次强龙卷风的持续时间及运动轨迹进行了调查研究；霍星远等[3]分析了2005年6月10日发生在辽宁朝阳县的一次龙卷风的天气状况及农作物受损情况；毛以伟等[4]调查了武汉5·14罕见龙卷风天气发生后的人员伤亡及树木、建筑的受损情况。在龙卷风的形成机制方面，Lewellen等[5]指出近地表流入层的差异是确定现有超级单体低水平中气旋是否产生龙卷风的关键因素；王沛霖等[6]指出在对流层深厚气层中气压低、风速大、层结不稳定、低层高温高湿和中下层位势不稳定的条件下，龙卷风更易发生；刘勇等[7]认为龙卷的形成与地面双冷锋结构形成有关，副冷锋不断向冷锋补充干冷空气导致地面上出现了较强气流；金巍等[8]指出低空气流由水平涡度扭曲为垂直涡度，在上升气流中被拉伸为中气旋，地面附近的背景垂直涡度由于增加的垂直拉伸作用而迅速加强，从而产生龙卷风。在龙卷风的风场特性方面，刘式适等[9]通过理论推导方式，利用气压梯度力、惯性离心力、黏性力三力平衡的柱坐标系下的大气动力和热力学方程组，求解了龙卷风的三维速度场，从理论上绘制出龙卷风的三维漏斗型结构，并指出龙卷风由惯性离心力造成的涡旋流和水平辐合辐散而引起的急流叠加而成；Kosiba等[10]通过对发生在香港国际机场的一场龙卷风进行现场实测，得出龙卷风的三维风场以及角动量的径向高度分布，发现龙卷风达到最大强度之前低水平角动量的增加；Karen等[11]测量了龙卷风的高度及不同高度处的风力。此外，随着对龙卷风的认识不断加深，已有研究通过风洞试验和计算流体动力学(CFD)数值模拟等方法实现了龙卷风风场的模拟，例如，田玉基等[12]运用龙卷风模拟器生成了5种涡流比的龙卷风风场，总结分析了龙卷风的切向风速、径向风速、竖向风速以及气压降的空间分布规律，并验证了龙卷风模拟器的有效性；王蒙恩等[13]运用数值模拟方法构建了龙卷风物理模拟器的数值计算模型，经过风洞试验验证后，指出随着涡流比的增大，龙卷风风场最大切向风速逐渐增大，涡核中心气压降明显降低，涡核半径随之变大，涡核中心附近切向风速的标准差变小；Verma等[14]通过CFD数值模拟并结合现场实测验证，指出计算域出口直径变化对龙卷风风场的影响大于龙卷风模拟风场总高度变化；Hangan和Kim [15]建立了一个小尺度龙卷风数值模型，可以在平稳条件下模拟龙卷风的切向风速、径向风速和竖向风速风场。为研究龙卷风作用下建筑物的风荷载与风致响应特性，国内外学者通过理论分析、现场实测、风洞试验及CFD数值模拟开展了一定的研究。例如，Fortuin等[16]通过分析景观尺度格局及其与地形属性、总体气候格局和土壤特征的关系，预测了龙卷风和剧烈雷暴对美国东南部森林的危害风险；Bain等[17]通过现场实测与风洞试验方法对龙卷风荷载作用下木屋盖的设计提供了经验；王进等[18]采用龙卷风模拟器，对受静止龙卷风涡旋影响下的立方体建筑物进行了压力测量，并指出涡流比为0.72的模拟龙卷风涡旋可以代表斯宾塞(Spencer)发生的一次改良藤田级数EF3级龙卷风的涡旋；汪大海等[19]通过建立三维风场揭示了龙卷风作用下输电线路风灾破坏的机理，得出在最大风速确定的情况下，响应随涡核尺度增大而减小，各跨不平衡张力引起的纵向的位移及反力响应会随着跨数的增多而增加，随绝缘子的增长而减小；Selvam [20]采用大涡模拟湍流模型对龙卷风作用下结构的风致效应进行了模拟，为基于大涡模拟的龙卷风数值模拟提供了经验；杨庆山等[21]提出龙卷风影响下高速列车倾覆概率分析方法，并结合实际案例进行分析，得出在列车数量不变的情况下，通过调整列车发车间隔可以降低19%的倾覆风险。综上所述，目前我国关于龙卷风的研究主要针对单个龙卷风展开，对龙卷风的时空统计特性及灾害指示物的破坏情况缺少较为全面且详尽的分析。研究中国龙卷风的统计特性并分析不同指示物的破坏情况及总体灾情情况对中国龙卷风数据库的建立及灾害指示物的确定具有一定的参考意义，对减少风致灾害损失具有一定的工程价值。

鉴于此，本文基于中国气象局编著的《中国气象灾害年鉴》记录的有关龙卷风的描述，对2003~2019年发生在我国境内的648次龙卷风的空间和时间统计特性进行分析，探讨建筑物、树木、农作物的受损情况及人员伤亡情况，利用熵权法综合评定的思想综合评价龙卷风灾情，为中国龙卷风数据库的建立及风致灾害防治提供参考。

2. 数据来源

本文采用的数据来源于《中国气象灾害年鉴》(2004~2020)，该年鉴在第2章气象灾害分述章节中，详细记录了上一年中我国龙卷风的发生情况。值得说明的是，虽然一些强度极弱、持续时间极短的龙卷风有尚未被记录的可能，但该年鉴为我国目前关于龙卷风较为完整的记录，对龙卷风的统计分析具有重要的参考价值。为了区分不同强度等级的龙卷风，本文采用日本改良藤田级数(JEF等级)进行划分。JEF等级将龙卷风强度划分为6个等级，即JEF0、JEF1、JEF2、JEF3、JEF4和JEF5，其中，JEF0等级最低，JEF5等级最高，随着等级的提高，龙卷风的破坏性增强。

表1给出了2003~2019年我国不同强度的龙卷风的数量统计情况。值得一提的是，在年鉴中，有250起龙卷风无具体描述，本文通过抽样调查分析了上述龙卷风的相关新闻报道及影像资料，发现其破坏程度均不明显，故本文将上述250起龙卷风全部归为JEF0级。由表1可知，在2003~2019年发生在我国境内的648次龙卷风中，JEF0级的龙卷风有307起，占比最大；JEF1、JEF2和JEF3级的龙卷风分别有126、195和19起；JEF4级的龙卷风有1起，为2016年6月23日发生在江苏省盐城市的龙卷风，该龙卷风为我国有记录以来强度等级最大的龙卷风；目前国内暂无JEF5级的龙卷风。

Table 1. Statistical table of the number of tornadoes with different intensity scales in China from 2003 to 2019

表1. 2003~2019年我国不同强度等级龙卷风的数量统计表

3. 龙卷风时空统计特性分析

3.1. 空间分布特性

Figure 1. Spatial distribution of tornadoes with different intensity scales in China from 2003 to 2019

图1. 2003~2019年我国不同强度等级龙卷风的空间分布情况

图1为2003~2019年我国不同强度等级龙卷风发生位置的空间分布图。由图可知，我国龙卷风大部分发生在东南沿海地区，部分发生在西北内陆地区。其中，JEF0、JEF1级龙卷风发生位置较广泛，JEF2级龙卷风在江苏省和湖北省的分布较多，JEF3级主要分布在河南省和江苏省，唯一的一起JEF4级龙卷风发生在江苏省。

图2所示为我国各省级行政区不同强度龙卷风发生频数的统计图。截至2019年，除北京、贵州未见龙卷风灾情外，其他30个省、自治区、直辖市在2003~2019年间均有龙卷风案例的记录。从数量上看，广东省龙卷风发生频数较多，共发生了49起JEF0级龙卷风、33起JEF1级龙卷风和16起JEF2级龙卷风；从强度等级上看，大多数省份JEF0级龙卷风占比最大，如山西省发生的龙卷风均为JEF0级。

Figure 2. A statistical map of the frequency of tornadoes with different intensities in various provinces of China from 2003 to 2019

图2. 2003~2019年我国各省级行政区不同强度龙卷风发生频数统计图

3.2. 时间分布特性

为了分析2003~2019年我国全国范围及各个省份(本文以广东和江苏为例进行说明)龙卷风的频数变化情况，本文将每年中不同强度等级龙卷风的发生频数汇总，并将其按发生时间进行划分，结果如图3所示。从总体上看，我国龙卷风发生的频率呈下降趋势，龙卷风发生频数最高的一年在2005年，累计有74起，而频数最低的一年在2019年，对应有14起，两年相差60起；其中，广东和江苏两地的龙卷风频数变化趋势与总体类似。值得一提的是，虽龙卷风发生次数呈现一定的下降趋势，但这并不意味着龙卷风产生的破坏程度减小，例如，2015年7月1日发生在黑龙江哈尔滨市呼兰区的龙卷风最大风速达到了70 m/s，导致2.6万人受灾；2016年6月23日发生在江苏盐城的JEF4级龙卷风造成了99人死亡，因此，龙卷风及其产生的危害仍应受到格外重视。

Figure 3. The frequency change diagram of tornadoes in the whole China and some provinces from 2003 to 2019

图3. 2003~2019年我国全国范围及部分省份龙卷风发生频数变化图

为了分析一年中不同月份龙卷风的发生情况，本文统计了2003~2019年各个月份我国龙卷风发生的频数，结果如图4所示。由图可知，龙卷风频发的月份主要集中在6~8月份，其中7月份的发生频数最多，这是由于7月份高温及强降水的气候特点为龙卷风的形成提供了条件[22]；而在12~2月期间，龙卷风的发生频数相对较少。

Figure 4. Variation diagram of occurrence frequency of tornadoes in each month from 2003 to 2019

图4. 2003~2019年各月份龙卷风发生频数变化图

在各月份的基础上，《中国气象灾害年鉴》(2004~2020)对222起龙卷风记录了其具体发生时间。基于此，本文对该222起龙卷风的发生时间进行了统计分析，如图5所示。由图可知，大多数龙卷风的发生时间集中在14~20时，而在其他时间段，龙卷风发生的次数相对较少。

Figure 5. Statistical graph of occurrence frequency of tornadoes in different time periods from 2003 to 2019

图5. 2003~2019年不同时段龙卷风发生频数统计图

4. 龙卷风灾情分析

龙卷风造成的灾害损失程度往往与其袭击的区域的特征紧密相关，包括地理位置、经济发展水平、人口密集程度等因素。即便是强度相当的龙卷风，其对不同地区带来的破坏和影响也存在显著差异。为了深入理解我国不同省份遭受龙卷风灾害时的具体损失情况，本文采用宏观视角，对各地龙卷风对建筑、树木、农作物以及人口所造成的损害进行了量化统计和全面评价。通过这一分析，我们旨在揭示龙卷风灾害的影响规模，并为灾害风险管理和应急响应策略的制定提供科学依据。

4.1. 建筑物受损情况分析

龙卷风过境后，建筑物的受损情况可以分为三类，即一般损坏、严重损坏和倒塌[23]，其破坏实例如图6所示。经统计分析，我国2003~2019年房屋遭受龙卷风后出现一般损坏和严重损坏的数量累计282,349间，倒塌的房屋数量累计86,263间。图7列举了不同结构形式的建筑物遭受龙卷风后的破坏情况。对比砖混结构、钢筋混凝土、钢结构和木结构建筑的破坏特点发现，钢筋混凝土遭遇龙卷风后破坏程度较低，主要以围护结构破坏为主，其抗风性能优于轻钢结构和砖木结构。因此，对龙卷风多发地区的不同类型建筑物加固也应有所侧重[24]。

(a) (b) (c)

Figure 6. Examples of damaged buildings after tornadoes: (a) Generally damaged (2020.5.18, Weifang, Shandong); (b) Seriously damaged (2017.8.21, Nenjiang, Heilongjiang); (c) Collapsed (2021.5.14, Caidian, Wuhan)

图6. 龙卷风发生后建筑物受损情况举例：(a) 一般损坏(2020.5.18，山东潍坊)；(b) 严重损坏(2017.8.21，黑龙江嫩江)；(c) 倒塌(2021.5.14，武汉蔡甸)

(a) (b)

(c) (d)

Figure 7. Examples of damaged buildings of different structure types after tornadoes: (a) Damage of brick-concrete house (2018.8.13, Jinghai, Tianjin); (b) Damage of reinforced concrete building (2021.8.25, Huludao, Liaoning); (c) Damage of steel structure facilities (2022.7.4, Huangpu, Guangdong); (d) Damage of wooden house (2023.2.26, Oklahoma, USA)

图7. 龙卷风发生后不同结构类型的建筑物破坏举例：(a) 砖混结构民房破坏(2018.8.13，天津静海)；(b) 钢筋混凝土楼房破坏(2021.8.25，辽宁葫芦岛)；(c) 钢结构设施破坏(2022.7.4，广东黄埔)；(d) 木结构房屋破坏(2023.2.26，美国俄克拉荷马)

为分析我国龙卷风导致的建筑物受损数量的空间分布情况，本文将各个省份龙卷风造成的一般损坏、严重损坏和倒塌的建筑物数量累加，得到各个省份建筑物受损数量总数，如图8所示。从整体上看，我国东部地区的建筑物受损数量大于西部，其中，江苏、安徽、黑龙江等地建筑物的破坏数量极多。

Figure 8. The distribution map of the number of damaged buildings after tornadoes (2003~2019)

图8. 龙卷风发生后建筑物受损数量分布图(2003~2019年)

4.2. 树木受损情况分析

龙卷风过境后，森林树木及城市行道树也会发生不同程度的损坏。统计《中国气象灾害年鉴》龙卷风导致的树木受损相关描述，发现2003~2019年间树木发生不同受损情况的总数为2,390,457棵，其中，树木被连根拔起(如图9(a))和拦腰折断(如图9(b))的数量分别为144,130和113,841棵，分别占6.03%和4.67%。除了年鉴中提到的上述两种受损模式外，结合龙卷风灾后相关新闻报道，树木还可能发生树干倾斜和枝条折断破坏，如图9(c)和图9(d)所示。

(a) (b)

(c) (d)

Figure 9. Examples of damaged types of trees after tornadoes: (a) Uprooted trees (2020.6.12, Yangzhou, Jiangsu); (b) Broken trees (2023.8.4, Tongliao, Inner Mongolia); (c) Tilted trunks (2021.5.14, Caidian, Wuhan); (d) Trees with broken branches (2020.7.6, Chongming, Shanghai)

图9. 龙卷风发生后树木受损类型举例：(a) 连根拔起(2020.6.12，江苏扬州)；(b) 拦腰折断(2023.8.4，内蒙古通辽)；(c) 树干倾斜(2021.5.14，武汉蔡甸)；(d) 枝条折断(2020.7.6，上海崇明)

图10中展示了2003~2019年间我国各个省份树木受损数量总数的空间分布，由图可知，我国大部分省份树木均遭受到龙卷风灾破坏，其中，四川、河南、江苏的树木破坏情况尤为严重，应引起足够重视。

Figure 10. The distribution map of the number of damaged trees after tornadoes (2003~2019)

图10. 龙卷风发生后树木受损数量分布图(2003~2019年)

4.3. 农作物受损情况分析

多起龙卷风的相关记录表明，龙卷风对农作物的生长会产生致命的影响。根据龙卷风灾后调查结果，笔者总结出农作物的失效模式主要有根部被吹离土壤和茎部倒伏两类，如图11(a)和图11(b)所示。农作物倒伏与其茎秆的刚度、稳定性密切相关[25]。根据龙卷风过境后农作物产量的降低程度，《中国气象灾害年鉴》统计了龙卷风发生后农作物的受灾、成灾及绝收面积[23]，分别为1296920.8、337,529和233,010公顷。根据2003~2019年间年鉴的相关记录，农作物受灾、成灾和绝收面积最大可分别达到18万、4.2万和15万公顷，农作物受灾、成灾和绝收面积位居前三的相关信息如表2所示。

(a) (b)

Figure 11. Examples of damaged crops after tornadoes: (a) Root uprooted crops; (b) Stem lodging crops

图11. 龙卷风发生后农作物受损情况举例：(a) 根部被吹离土壤；(b) 作物茎部倒伏

Table 2. Statistical table of relevant information on the top three areas of crop disaster, disaster and crop failure after tornadoes from2003 to 2019

表2. 2003~2019年龙卷风发生后农作物受灾、成灾和绝收面积位居前三的相关信息统计表

图12展示了我国各个省份2003~2019年间农作物受灾、成灾和绝收面积总数的空间分布。由图可知，受损农作物主要分布在胡焕庸线以东，其中安徽、黑龙江、吉林的农作物受损情况较严重。

4.4. 人员伤亡情况分析

龙卷风造成的建筑物倒塌、输电线路破坏及车辆、船只翻起会对人的生命造成一定威胁。据统计，2003~2019年龙卷风共造成388人死亡、7485人受伤、受灾人数¹达20,033,080人次。将2003~2019年每起龙卷风造成的死亡及受伤人数累加，并统计伤亡人数在各个省份的空间分布情况，结果如图13所示。由图可知，湖北、江苏和安徽的伤亡人数较多。

Figure 12. Distribution map of crop damage after tornadoes (2003~2019)

图12. 龙卷风发生后农作物受损情况分布图(2003~2019年)

Figure 13. Distribution map of casualties after tornadoes (2003~2019)

图13. 自龙卷风人员伤亡情况分布图(2003~2019年)

5. 龙卷风灾情综合评价

5.1. 龙卷风灾情综合评价模型

为评估龙卷风造成的灾情严重程度，本文选择熵权法结合综合评定的思想[27]，通过确定灾害评价指标及其权重，对龙卷风灾情进行分析及评价。具体过程如下：

设$X_{ij}$ 为龙卷风灾情指标记录的真实值，$X_{\min }$ 、$X_{\max }$ 分别为某起龙卷风案例中此指标的最小值和最大值，对$X_{ij}$ 进行无量纲化得到$Y_{ij}$ ，即

$Y_{ij}=\frac{X_{ij}-X_{\min }}{X_{\max }-X_{\min }}$ (1)

为确定该指标的变异情况，计算该指标在相同案例、所有指标数据中的占比$p_{ij}$ ，即

$p_{ij}=\frac{Y_{ij}}{{\displaystyle {\sum }_{i=1}^n{Y}_{ij}}}\left(i=1,\cdots ,n;j=1,\cdots ,m\right)$ (2)

式中，n为指标数，m为案例总数。

根据各指标在不同案例中的比值，分别计算各个指标的信息熵$H_j$ 及权重${\omega }_j$ ，即

$H_j=-\left({\displaystyle {\sum }_{i=1}^n{p}_{ij}\ln p_{ij}}\right)/\ln n$ (3)

${\omega }_j=\frac{1-H_j}{n-{\displaystyle {\sum }_{j=1}^m{H}_j}}$ (4)

根据各指标的权重值可进一步求出综合受灾指数$G_i$ ，即

$G_i={\displaystyle {\sum }_{j=1}^m{\omega }_j}\times Y_{ij}\left(j=1,2,\cdots ,m\right)$ (5)

考虑龙卷风指标记录不全的情况，计算龙卷风灾情指数$S_i$ ：

$S_i=\frac{G_i}{{\displaystyle {\sum }_{j=1}^m{\omega }_j}}$ (6)

5.2. 龙卷风灾情综合评价结果

Figure 14. Spatial distribution map of tornado disaster index

图14. 龙卷风灾情指数空间分布图

采用熵权法对2003~2019年发生的648起龙卷风灾情案例进行计算，得到我国龙卷风灾害评价指标及权重，如表3所示。

Table 3. Disaster evaluation index and weight of tornadoes

表3. 龙卷风灾害评价指标及权重

将2003~2019年发生的648起龙卷风灾情信息代入到龙卷风综合评价模型中，得到了不同龙卷风造成的灾情指数，其空间分布如图14所示。由图可知，龙卷风灾情主要集中在胡焕庸线以东，其中，华东、华南和东北地区灾情最为严重，华北，华中地区次之，西北、西南和港澳台地区灾情较轻。

6. 结论

本文根据对2003~2019年《中国气象灾害年鉴》有记录的648起龙卷风灾情进行统计分析，得出以下结论：

从空间上看，我国龙卷风主要分布在东南沿海及东北的平原和高原地区，江苏、广东的分布数量较多；从时间上看，我国龙卷风的发生次数呈逐年递减的趋势；在一年内，龙卷风主要发生在6、7、8月份的夏季，其中7月份发生次数最多；在一天内，龙卷风主要发生在14~16时，此时段内龙卷风发生次数明显高于其它时段。

在2003~2019年间，我国房屋遭受龙卷风破坏的累计损坏数为282,349间，累计倒塌数为86,263间，我国东部房屋破坏数量明显大于西部，江苏、安徽和黑龙江等地房屋的破坏数量极多；我国树木遭受龙卷风破坏数量为2,390,457棵，主要失效模式有连根拔起、树干折断、树干倾斜及枝条折断，四川、河南和江苏的树木受龙卷风破坏较多；我国农作物遭受龙卷风破坏受灾、成灾和绝收面积分别为1296920.8公顷、337,529公顷和233,010公顷，主要失效模式有根部被吹离土壤和茎部倒伏，安徽、黑龙江和吉林的农作物受损情况较严重；龙卷风共造成388人死亡，7485人受伤，受灾人数达20,033,080人次，相同强度的龙卷风对不同地区造成的影响不同。通过熵权法综合评价的思想对龙卷风灾情进行评价得出，龙卷风灾情主要集中在胡焕庸线以东，其中，华东、华南和东北地区灾情最为严重。

基金项目

国家电网有限公司科技项目(5200-202318147A-1-1-ZN)。

NOTES

^*通讯作者。

¹受灾人数：行政区域内因自然灾害遭受损失的人口数量(含非常住人口) [26]。

参考文献