引言

设施农业生产是农业增产、农民增收的主要手段，对于保障农产品有效供给、增强农业综合生产力以及实施乡村振兴战略具有重要意义。20世纪80年代以来，中国以日光温室、塑料大棚为代表的设施农业发展迅速(陈思宁等，2014)，据农业农村部统计，2008—2016年期间，全国(不含港澳台)设施农业总面积由81.312万hm²增加至208.288万hm²，以平均每年14.108万hm²的幅度增长，其中日光温室和塑料大棚总面积占比达95%以上。尽管中国设施农业取得了长足进步，但农业设施抵御自然灾害能力依然薄弱，尤其是在气候变化背景下，极端灾害性天气频发，给设施农业生产带来的损失显而易见(李秀芬等，2007；苗克红等，2016)，严重制约设施农业可持续发展。国外发达国家设施农业标准化程度高，设施环境综合调控及农业机械化技术发达，设施农业受气象灾害影响较小，其研究主要侧重于经济风险、社会风险和环境风险(Piers，1998；Carter，1995)。我国对于设施农业的研究主要集中在气候适应性和气候区划(李德等，2013；张亚红，2003；张明洁和赵艳霞，2013)、气象灾害风险分布(蔡冰等，2011；杨再强等，2012；陈思宁等，2017；黄川容等，2012)及灾害影响和防御措施方面(柳军等，2018；刘佳等，2018)，而对于灾害指标的研究甚少。随着天气预报业务的发展，数值天气预报为气象灾害快速、精准预警服务提供了更广阔的发展空间。大风天气时有发生，准确的预报预警信息可使生产者提前做好防范，有效减少风灾损失，因此，制定科学合理的风灾指标尤为重要。设施农业风灾指标常规的获取方法有：现场观测法、计算流体力学法(CFD法)和风洞试验法(Lukiantchuki et al，2020；任兰红等，2019；张波，2013；杨振超，2006)。现场观测法需要耗费大量的时间和人力，风洞试验法受试验条件所限，试验模型、灾害发生环境与实际情况存在明显差异，且因耗资较大，研究也较少，而CFD法主要用于温室内部环流场的模拟。目前针对设施农业风灾的预报预警服务，大多采用蒲福风力等级作为风灾指标(孙治贵等，2018；孙智辉等，2011)，但此标准未考虑各地气候特点差异，也未与农业设施抗风设计标准及风灾影响程度相结合，不利于设施农业生产者有针对性地采取防御措施，导致风灾预警和防御的有效性大打折扣。可见，通过以上方法确定的风灾指标在实际应用中均存在一定局限性。

本文基于风对农业设施致灾原理及其影响程度划分风灾等级，以我国各地气候特点及主要农业生产设施结构设计使用年限为依据，探讨了日光温室和塑料大棚风灾指标确定方法，并利用实际灾情资料对风灾指标进行验证，以期为各地农业设施风灾预报预警和风险评估、农业防灾减灾以及设施农业科学管理、减少生产成本提供参考依据。

1 资料及方法 1.1 资料及来源

本研究所用气象资料为1989—2018年全国2467个气象观测站10 m高度处日最大风速，即从给定时间段(1 d)内选取任意10 min平均风速的最大值，用于分析计算不同重现期最大风速；2009—2018年全国1714个气象观测站10 m高度处日极大风速，即给定时间段(1 d)内3 s最大瞬时风速，用于1.5节中构建日最大风速与极大风速间关系模型，以上风速资料均来自国家气象信息中心。2009—2018年和2020年灾情资料用于风灾指标验证，来自中国气象局气象灾害管理系统。

1.2 风灾等级划分

在设施农业生产中，当风速达到一定量级，对日光温室和塑料大棚的直接影响是掀翻棚膜、棚体倾斜，随着风速的增加，超过设施结构设计的风荷载，设施主体结构将受到损坏，甚至垮塌(龚伟等，2017；陈松来，2009)。因此，本文基于以上风灾对农业设施的影响，将其划分为两个等级。由于日光温室建造时，墙体的承载力一般都大于其可能承受的荷载量，因此对于日光温室，主体结构主要考虑钢骨架结构。

1.3 设施结构设计使用年限确定

我国温室、大棚由于其使用材料、功能以及用户需求的差异，在重要性、安全度、实际使用年限等方面各不相同(孙德发等，2002)，因此，针对不同设施结构设计使用年限，一般规定基本风压的重现期为5~30 a，即风速的保证概率为80.0%~96.7%。本研究中日光温室指三面为保温围护墙体，单坡面以钢管为骨架支撑、以塑料薄膜覆盖的农业生产设施，使用年限一般不超过20 a(李中明等，2011)；塑料大棚指以钢或钢、竹木混合为支撑，以塑料薄膜覆盖的拱形农业生产设施，使用年限一般不超过10 a。综合以上日光温室和塑料大棚结构设计使用标准，以及实际使用寿命和风速保证率，本文对于日光温室和塑料大棚主体结构的设计使用年限，分别采用20 a和10 a的重现期标准(王笃利等，2005；闫俊月等，2014)，相应风速保证率分别为95%和90%。

1.4 T年重现期最大风速估算方法

基于1989—2018年逐日最大风速，计算各站点年最大风速x(单位：m·s^-1)，采用广义极值分布函数(Coles，2001)对逐年最大风速x的分布进行拟合：

$ \begin{aligned} &g(x)= \begin{cases}\frac{1}{\delta} \mathrm{e}^{-(1+k z)-\frac{1}{k}(1+k z)-1-\frac{1}{k}} & k \neq 0 \\ \frac{1}{\delta} \mathrm{e}^{-z-\mathrm{e}^{-z}} & k=0\end{cases}\\ &z=\frac{x-\mu}{\delta} \end{aligned} $

(1)

式中：g(x)为拟合函数, k为形状参数, δ为尺度参数, μ为位置参数。参数k、δ、μ的估计采用最大似然法，函数拟合优度检验利用Kolmogorov-Smirnov(K-S)检验法。

基于式(1)对参数k、δ、μ的估计结果，T年重现期最大风速U通过式(2)做估算：

$ U= \begin{cases}\mu-\frac{\delta}{k}\left\{1-\left[-\ln \left(1-\frac{1}{T}\right)\right]^{-k}\right\} & k \neq 0 \\ \mu-\delta \ln \left[-\ln \left(1-\frac{1}{T}\right)\right] & k=0\end{cases} $

(2)

式中：U为T年一遇最大风速(单位：m·s^-1)；T为重现期，根据1.3节的规定取值。

1.5 风灾指标阈值计算方法

农业设施抵御风灾的能力由其结构本身的风荷载决定，风荷载的大小取决于基本风压，而各地基本风压是一个关于重现期、风速时距和离地高度三个变量的复合单调函数(王笃利等，2005)，表示如下：

$ W=\frac{1}{2} r U^2 \beta^2 \omega $

(3)

式中: W为农业设施结构设计基本风压(单位：kN·m^-2)；r为空气密度(单位：t·m^-3)；U为不同重现期最大风速(单位：m·s^-1)；ω为标准地貌下的风压高度变化系数；β为时距换算系数，对于日光温室和塑料大棚，此类轻型结构设计宜采用3 s时距的最大瞬时风速(即极大风速)作为基本风压取值的基础要素，因此本研究β即为极大风速(3 s时距)与最大风速(10 min时距)的换算系数，通过构建最大风速(x)与极大风速(y)间关系模型y=f(x)得到，因此式(3)可改写为：

$ W=\frac{1}{2} r f U^2 \omega $

(4)

根据伯努利方程，实际来流风压可表示为：

$ W_0=\frac{1}{2} r V^2 $

(5)

式中：W₀为实际来流风压(单位：kN·m^-2)，V为实际风速(单位：m·s^-1)。

基于式(4)和式(5)，当W₀=W，即实际来流风压达到农业设施结构设计的基本风压时，设施主体结构可能受到损坏，此时的风速V即为造成设施主体结构受损的风灾指标阈值，即：

$ V_j=\sqrt{f\left(U_j\right)^2 \omega_j} $

(6)

式中：f(U_j)为不同重现期最大风速对应的极大风速(单位：m·s^-1)；j为设施类型，j=1代表日光温室，j=2代表塑料大棚；ω_j为风压高度变化系数，按照《农业温室结构荷载规范》(中华人民共和国住房和城乡建设部，2016b)中4 m(日光温室)和3 m(塑料大棚)高度标准取值，ω1=0.76，ω₂=0.7。

根据1.2节风灾等级划分标准，各等级风灾指标阈值按照式(7)计算：

$ V_{i, j}=V_j \alpha_i $

(7)

式中: V_i，j为不同等级日光温室和塑料大棚风灾指标阈值(单位：m·s^-1)；i为风灾等级(分别为Ⅰ和Ⅱ)；α_i为不同等级风灾影响系数，根据张波(2013)，日光温室和塑料大棚顶部和侧面局部往往承受的风压较大，即当风速达到一定量级，以上部位棚膜首先被损坏，而风在设施表面产生的风压可通过风荷载体型系数体现，根据中华人民共和国住房和城乡建设部(2016b)相关规范，0°风向时设施顶部和垂直风向的设施表面风荷载体型系数为0.8，同时参考蒲福风力等级影响以及风力为6~8级时相邻等级间风速阈值比值，由此以Ⅰ级风灾的影响系数α_Ⅰ=1.0为标准，Ⅱ级风灾影响系数α_Ⅱ=α_Ⅰ×0.8=0.8。

2 结果与分析 2.1 不同重现期最大风速

利用1.4节中介绍的方法计算不同重现期最大风速，并检验函数模拟拟合优度，结果显示：全国2 467个气象观测站点中，共1 714个站点数据拟合结果通过检验，基于上述1 714个站点计算20 a、10 a重现期最大风速，结果见图 1和图 2。可见，最大风速高值区主要分布在西藏中东部、新疆中东部、青海、内蒙古、黑龙江、吉林和辽宁大部、山西中北部、河北南部、山东半岛、华南和东南沿海地区等地，其中西藏北部、内蒙古北部、青海西部、新疆北部等地受地形等因素影响，华南沿海地区受台风影响，20 a、10 a一遇最大风速分别在21 m·s^-1和17 m·s^-1以上；最大风速低值区主要分布在长江以南的非沿海地区，其中湖北南部、湖南、江西、华南西部和北部、四川盆地等地20 a、10 a一遇最大风速均在14 m·s^-1以下；其他地区属于风速较高值区，大部地区20 a、10 a一遇最大风速在14~17 m·s^-1。可见，在不同地区，不考虑设施体型、高度等风荷载参数不同的情况下，即使是相同设计使用年限的温室或大棚，其基本风压、抗风能力也会有所差别，在温室和大棚建造时，应加以考虑。

2.2 最大风速与极大风速换算系数

考虑到各地气候特点不同，以及近年来气候较20世纪80年代之前的变化(孟丹等，2019；敖雪等，2020)，如最大风速与极大风速换算系数按照经验常数(朱瑞兆，1976)取值，在时间和空间上将会存在一定差异，为使指标更加准确，本文利用1.5节中介绍的β计算方法，对各站点最大风速与极大风速换算系数进行估算。分析结果显示，两者间存在f(x)=kx+b的线性关系，且所有站点线性方程均可通过P＜0.05的显著性水平检验，全国85%的站点参数b的绝对值在[0，1]，99%站点参数b的绝对值在[0，2]。因此，拟将参数k作为最大风速U与极大风速换算系数，即极大风速为Uk。经方差分析检验，f(x)与Uk无明显差异，因此，k可作为最大风速与极大风速换算系数，尤其是在极大风速资料获取困难、无法构建极大风速与最大风速转换方程的情况下使用，由此图 3给出了k的空间分布。可见，k值在长江以南地区和新疆较大，大部在1.6~1.9；在东北地区、内蒙古、青海、河南、安徽、广西等地最小，普遍在1.1~1.6，与重现期风速大小的空间分布大致相反。从全国来看，77%站点的k在1.4~1.8，全国平均值为1.6。

2.3 风灾指标计算结果

由于日光温室以冬季保温为目的，主要分布在长江以北和高海拔地区，塑料大棚则在全国范围内均有分布。根据1.5节方法，计算日光温室和塑料大棚风灾指标阈值，Ⅰ级指标阈值分布见图 4和图 5。可见，风灾指标阈值高值区主要分布在西藏、新疆、青海、甘肃中北部、宁夏、内蒙古、黑龙江、吉林、辽宁、山西中北部、河北中南部、山东半岛等地，以上地区日光温室和塑料大棚Ⅰ级风灾指标分别为极大风速≥24.0 m·s^-1和≥21.9 m·s^-1；指标阈值低值区主要分布在河南西部以及长江以南的湖北南部、贵州东北部、重庆东南部、湖南西北部、江西、广西中北部等地，日光温室和塑料大棚(以塑料大棚为主)Ⅰ级风灾指标分别为极大风速在[15.2，19.0] m·s^-1和[13.5，16.9] m·s^-1；其他地区属于指标较高值区，主要为黄淮海地区、西北地区东南部和长江以南非沿海的大部分地区，Ⅰ级风灾指标位于高值区和低值区的风灾指标之间。当极大风速达到或超过以上标准，日光温室和塑料大棚主体结构可能损坏或倒塌。

指标高值区的日光温室和塑料大棚Ⅱ级风灾指标分别为极大风速在[19.2，24.0) m·s^-1和[17.5，21.9) m·s^-1，低值区指标分别为极大风速在[12.2，15.2] m·s^-1和[10.8，13.5] m·s^-1，其他地区属于风灾指标较高值区，当极大风速达到或超过以上标准，日光温室和塑料大棚棚膜可能被损坏或掀翻。由于风灾Ⅱ级指标阈值与Ⅰ级指标阈值分布相同，只存在量级上的差异，因此不再给出空间分布图。可见，风灾指标存在空间差异，如采用同一抗风标准建造日光温室或塑料大棚，一定程度上可能会造成抗风能力不足而影响使用，或者抗风标准过高浪费生产成本。

3 风灾指标检验 3.1 合理性检验

从2009—2018年全国风灾灾情数据中共提取出754个农业设施受灾，根据风灾发生时间、地点，从灾情记录和观测数据中提取每个风灾过程中的日极大风速，将其与风灾指标进行相关性分析。结果显示，各风灾样本的极大风速观测值与风灾指标显著相关(n=754，r=0.372，P＜0.05)。说明各地风灾发生时的极大风速与该地风灾指标具有较好的一致性。此外，754个检验样本中极大风速观测值≥10.8 m·s^-1(塑料大棚风灾Ⅱ级指标阈值下限)的占比为97.6%，一定程度也反映了指标的合理性。

3.2 准确性检验

因历史灾情资料中大多未明确记录承灾体类型是日光温室或是塑料大棚，受灾程度为棚膜受损还是主体结构受损，因此很难对灾害指标的准确性进行验证，这也是目前很多指标缺少验证的主要原因。因此，本文选择南、北方设施农业主要分布和适宜发展地区(张亚红，2003)中的代表省份山东省(黄淮海地区)、甘肃省(西北地区)，以及江西省、湖北省和江苏省(长江中下游地区)，从2020年风灾灾情中提取有详细记录的24个农业设施受灾样本，进一步对风灾指标和等级进行验证，并规定风灾样本中观测极大风速与风灾等级及其对应的指标范围完全吻合即为“符合”，相差一个等级为“基本符合”，否则为“不符合”，验证结果见表 1。表 1显示，24个风灾样本中，17个样本的验证结果为“符合”；4个样本极大风速观测值与日光温室、塑料大棚的Ⅱ级风灾指标吻合，但未达到1级指标即出现“部分结构损坏”，验证结果为“基本符合”；3个样本在极大风速观测值未达到风灾Ⅱ级指标即出现“棚膜受损”，验证结果为“不符合”，但指标阈值与极大风速观测值相差较小，均在2.0 m·s^-1以内。可见，本文确定的日光温室和塑料大棚风灾指标较为准确。

4 结论与讨论

根据我国主要农业生产设施(日光温室和塑料大棚)结构抗风设计标准，基于气候特点及风对农业设施致灾原理，确定了日光温室和塑料大棚风灾等级及相应等级风灾指标。日光温室和塑料大棚风灾指标阈值呈北方、沿海地区及高海拔地区大、南方及低海拔地区小的特点。指标阈值高值区主要集中在西北地区中西部、内蒙古、东北地区和华北部分地区以及沿海地区，日光温室和塑料大棚风灾Ⅰ级指标分别为极大风速≥24.0 m·s^-1和≥21.9 m·s^-1，Ⅱ级指标分别为[19.2，24.0) m·s^-1和[17.5，21.9) m·s^-1，上述地区农业设施一般宜为抗风和保温性能好、厚度较厚的钢、镀锌管等材料；湖北南部、江西、湖南、广西和四川盆地东部等地属于指标阈值低值区，塑料大棚Ⅰ级风灾指标为极大风速[13.5，16.9] m·s^-1，Ⅱ级指标为[10.8，13.5] m·s^-1，上述地区农业设施主体结构一般为厚度较薄的经济型材料，农业设施以塑料大棚为主，如设施内栽培经济价值较高的作物，则多采用玻璃或大型连栋温室等抗风性能更好的设施；全国其余大部地区为指标阈值较高值区，也是我国设施农业主要分布区域，设施主体结构一般为可综合考虑抗风性能和经济投入的标准日光温室和塑料大棚。张波(2013)通过风洞试验得到塑料大棚和日光温室的致灾临界极大风速分别为14.5 m·s^-1和18.8 m·s^-1，在本研究塑料大棚和日光温室风灾Ⅱ级指标较高区域范围之内；张淑杰等(2015)对东北地区2013年大风灾害统计结果显示，80%以上风灾样本导致掀棚的极大风速在20 m·s^-1以上；郭小芹等(2016)给出了甘肃省武威市设施农业轻度、中度和重度风灾指标的极大风速分别为13.8~17.1、17.1~20.7和20.7~24.4 m·s^-1，与本文的研究结果较为一致。此外，现有地方标准(山西省质量技术监督局，2016；辽宁省质量技术监督局，2015；浙江省市场监督管理局，2021)中，山西省塑料大棚和日光温室棚膜被吹破的指标分别为风力7~8级和8~9级，辽宁省风力达8~9级时启动风灾橙色预警，浙江省塑料中棚和单体钢架大棚中度和重度受损的指标分别为风力6~7级和8级以上，均与本文研究结果较为吻合。

本研究风灾指标可为不同地区日光温室和塑料大棚风灾预报预警和评估提供依据和参考，同时生产者可根据预报风速对应的不同风灾等级，采取相应防御措施，提高防灾避险的有效性，减少不必要的生产成本浪费。当预报的极大风速达到Ⅱ级风灾指标时，生产者可通过适时关闭通风口、提前修补破损棚膜，增加设施密封性，并通过采取压严和抻平棚膜、增加压膜线等措施，提高设施防风性，北方日光温室在没有降雨、降雪的大风天气，可适当增加草帘厚度，同时用绳索、纱袋等固定、压牢，减少棚膜受灾风险；当预报的极大风速达到Ⅰ级风灾指标时，除采取Ⅱ级风灾相应防御措施外，还应采取加固或增加棚室立柱和骨架等措施，提高设施主体结构抗风能力，必要时可适时揭膜或破膜保棚，减少因垮棚而造成的经济损失，对于北方日光温室，还可采取在后墙设置防风膜等风灾防御措施。此外，基于本文风灾指标，可计算全国不同区域日光温室和塑料大棚基本风压，其结果也可为设施农业建设投资和合理布局提供参考，以防止因设施抗风能力不足导致受灾风险增大，或因风压标准过大导致建造成本浪费。

在进行指标合理性检验时，2.4%的灾害样本在未达到风灾指标的情况下出现了不同程度灾情，主要原因之一可能是在风灾发生期间，同时伴随其他灾害性天气，之二可能是承灾体为闲置的年久失修大棚或者简易拱棚等，基于以上情况，今后风灾指标在构建时还应进一步考虑利用其他致灾因子和周期因子对指标加以修正；还有部分样本的观测风速明显偏大，可能是受客观条件所限，该地区样本较少，有记录的灾情多为受灾极为严重的情况。利用2020年灾情资料对风灾指标准确性的检验结果总体较好，但甘肃省风灾指标较样本观测值总体偏大，导致风速未达到风灾指标，即出现设施受灾的情况，其原因可能是样本中设施使用时间较久或疏于管理，也可能是当地设施结构设计标准与本研究存在一定差异，具体原因还需进一步调研确定。

本研究是以我国各地气候特点及主要农业生产设施抗风结构设计标准为依据，在正常生产和管理的条件下，给出了北方日光温室和全国塑料大棚不同等级风灾指标构建方法，虽然一定程度上解决了目前农业设施风灾预报预警、评估指标针对性不足、农业设施结构复杂风灾指标难以确定的问题，但由于农业设施种类和用途繁多，即使是同种设施类型，其在不同地区间的规模、覆盖材料、管理水平等方面仍存在差异，因此，其他类型和抗风设计标准的农业设施风灾指标可参考本文方法进行计算，在应用时也应根据当地设施农业生产管理水平、设施结构材料本身特点等实际情况进行适当调整。此外，根据工程技术规范(中华人民共和国住房和城乡建设部，2016a)，农业设施用地一般选择在背风、非冰雹频发、平整的地块，地形坡度不宜大于10°，而实际生产中，受客观条件所限，部分农业设施并不能够按照规范建造，且在风灾发生时还可能伴随冰雹、暴雨等灾害性天气，更易造成棚膜破损、主体结构不稳，一定程度上会降低设施的抗风能力，针对以上情况，建议在实际应用时对风灾指标进行适当折减。