引言

全球变暖加剧了气候系统能量和水分循环相互作用的变化，水分平衡变化导致极端旱涝事件频发(张霞等，2018；倪深海等，2019)。干旱是最常见、最严重的自然灾害之一，由于其影响范围广、持续时间长、出现频率高，常常给人类生产、生活带来巨大的损失(金燕等，2018；张强等，2020)，据统计，气象灾害引起的损失在各类自然灾害造成的总损失中约占85%，而干旱灾害又占气象灾害的50%左右。中国是干旱灾害发生频率最高且影响最严重的国家之一，旱涝格局呈现为北方易受旱灾影响、南方旱涝并发的特征(马柱国等，2018)，大范围的干旱灾害连年发生，农作物每年平均受旱面积为2.09×10⁷ hm²，最高年份达4.05×10⁷ hm²，平均干旱成灾面积为8.87×10⁶ hm²，最高年份达2.68×10⁷ hm²。每年造成的粮食减产从数百万吨到3000多万吨，干旱的直接经济损失高达每年440亿元(Su et al，2018)。干旱是个错综复杂的过程，根据不同的研究目标选择适宜的干旱指标来预测、评估旱情的发生、发展和变化情况，才能更好地为水资源的合理利用和科学抗旱提供依据(Obasi，1994；谢五三等，2019；左冰洁和孙玉军，2019；杨歆雨和李栋梁，2019；张强等，2021)。考虑实际业务需求，2012年国家气候中心在业务中推行了改进的气象干旱综合指数(MCI)，是综合考虑前期不同时间段降水和蒸散对当前干旱的影响而构建的一种单站干旱指数(杨帆等，2015；廖要明和张存杰，2017；王素萍等，2020)，并发现MCI指数对干旱过程的刻画效果最好(王荣江等，2021；谢五三等，2021)。

对于干旱事件需要多个干旱特征量来表征，如干旱历时、干旱强度、强度峰值和干旱空间范围等，目前对干旱的描述多侧重于干旱历时和干旱强度两个方面，中国气象局(2021)发布的气象行业标准(以下简称“行标”)中规定了区域性干旱过程的确定方法、评估指标和等级划分，根据该标准可以识别单站和区域干旱过程的历时和强度。另外，用于评价气候事件严重程度的重要指标——重现期，也备受人们的关注。近年来，学者们使用全历时重现期方法(冯国章，1995)、多元平稳随机过程游程分析(张学成等，1996)、模糊聚类分析(朱廷举和胡和平，2001)、马尔科夫链模型(孙鹏等，2014)等方法对重现期进行研究，但是不同的干旱变量(如：干旱历时和干旱强度)往往是用不同的分布函数进行拟合的，使得利用传统的统计方法建立多变量模型变得困难。Sklar(1959)提出的Copula作为一种联合函数近年来被广泛地引入到水文干旱、气象干旱等研究中。Copula函数能较好地对干旱问题中的多特征量进行多元拟合，并对多变量系统进行条件概率和重现期分析。许多研究也表明，游程理论和Copula函数的联合应用能够为多变量干旱特征分析提供有效的研究途径(徐春晓等，2011；周玉良等，2011)，拟合得到的干旱历时和干旱强度的联合概率分布，其计算的重现期与实际受旱状况相符(李天水等，2016)。

福建地理位置特殊，天气气候复杂，属于亚热带湿润季风区，虽然地处中国东南沿海，降水充沛，但由于降水量年际间和季节间的不均匀性以及季内分布的脉冲性，使季节性连旱发生频繁，干旱季节性、区域性很强(陈家金等，2006)。21世纪以来，福建省出现多次严重或较严重的干旱事件，尤其近年来干旱呈现多发频发态势：2003年，发生自1939年以来最严重的夏秋冬连旱，全省农作物受旱面积达17.50万hm²，361个乡镇187万人出现饮水困难；2018年，发生春夏连旱，粮食作物受灾较重，特别是对早稻播种和秧苗生长不利，部分地区早稻改种旱地作物；2019年，发生秋冬连旱，导致水稻、玉米、薏米、锥栗、茶叶等减产，全省13.42万人受灾，农作物受灾面积为1.753万hm²，绝收面积为0.11万hm²，直接经济损失达1.52亿元；2020年，发生夏秋连旱和冬旱，福州、宁德等地出现明显人饮水困难和工业用水紧张情况。可见，干旱灾害已成为福建防灾减灾的重要内容之一。

回顾已有的工作发现，前人通过日雨量≤2 mm连续日数(连旱日数)统计，对福建干旱气候态已有足够的认识，但由于连旱日数不属于随机变量，无法开展重现期的研究。另外，张容焱等(2019)分析改进的MCI指数在福建的适应性，认为该指数适用于福建干旱事件的研究，构建了适合福建的气象干旱过程监测预警和强度评估指标(以下简称“本地算法”)，并开展了业务应用。本文根据“行标”及MCI指数的干旱等级(中华人民共和国国家质量监督检验检疫总局和中国国家标准化管理委员会，2017)，识别1961年1月至2020年12月福建省各站点的干旱过程和全省干旱过程。但是在研究中发现，根据“行标”规定方法，用66个站的MCI均值来表征全省区域的MCI值，识别出来的干旱过程与《中国气象灾害大典福建卷》(温克刚等，2007)及《福建水旱灾害》(福建省人民政府防汛抗旱指挥部办公室，2018)中记载的历史过程存在一定差距，而“本地算法”建立的全省MCI综合指数序列，识别出来的干旱过程更加接近历史过程，因此本文在研究全省区域干旱时使用“本地算法”，并进行重现期特征研究，以期对福建干旱风险的进一步认识提供参考。

1 资料和方法 1.1 资料

根据GB/T 20481—2017气象干旱等级(中华人民共和国国家质量监督检验检疫总局和中国国家标准化管理委员会，2017)干旱综合指数计算方法，利用福建省66个国家气象站，1961年1月至2020年12月逐日气温、降水等数据，计算逐日MCI指数，并划分逐日干旱等级。

1.2 方法介绍 1.2.1 干旱过程的确定

(1) 单站干旱过程

根据单站逐日MCI指数，采用“行标”中单站干旱过程监测方法，计算单站干旱历时、干旱累计强度和干旱过程强度：当单站(区域)日干旱强度达轻旱及以上等级，且连续15天及以上，至少有1天干旱强度达中旱及以上等级，则发生一次单站干旱过程。单站干旱时段内第一次出现轻旱的日期，为干旱开始日；干旱过程发生后，但出现连续5天干旱等级为无旱时，则单站干旱过程结束，干旱过程结束前最后一天气象干旱等级达轻旱及以上的日期为结束日。开始到结束日(含结束日)的总天数为干旱过程持续日数，即干旱历时。

累积干旱过程强度是干旱指数强度与持续时间的综合，而：

$ D(n)=n^{a-1} \sum\limits_{i=1}^n I_{\mathrm{d}}(i) $

(1)

式中：D(n)为累积干旱强度，n为干旱过程日数，I_d(i)为干旱过程内第i天的单站日干旱强度绝对值，a为权重系数，一般取0.5~1.0，推荐取0.5。

通过滑动单站干旱过程内持续干旱天数计算单站累积干旱强度，取单站干旱过程内最强累积干旱强度作为单站干旱过程强度，而：

$ Z=\operatorname{MAX}_{k=1, m}\left\{\operatorname{MAX}_{n=1, k}[D(n)]\right\} $

(2)

式中: m为干旱过程总天数，n为干旱过程内干旱持续天数1≤n≤m；D(n)为单站累积干旱强度；MAX_{n=1, k}()为通过不断滑动单站干旱过程内持续天数n，比较寻找n=1, 2, …, k天的最大单站累积干旱强度；MAX_{k=1, m}( )为通过不断滑动单站干旱过程内干旱日k(1≤k≤m)，比较寻找包含第k天的最大累积干旱强度。

(2) 区域干旱过程

采用“本地算法”，建立全省区域MCI综合指数序列，分析区域干旱等级和过程判断指标，具体定义和算法等详见张容焱等(2019)。

(3) 干旱过程命名规则

干旱过程严格按照福建气候季节命名，以3—6月为春旱，7—9月为夏旱，10—11月为秋旱，12月至翌年2月为冬旱，冬旱年份命名以12月的年份为冬旱年。

跨季节连旱年份命名规则：以起始一季干旱的年份命名。如：秋冬连旱时间虽然跨入新年，年份命名为上一年的秋冬旱，冬春连旱也命名为上一年的冬春旱。

1.2.2 Copula函数

当选用干旱历时和干旱强度两个干旱特征变量共同来描述干旱过程时，则需要解析两者之间的联系，即计算它们的联合概率分布函数。目前，应用较广泛的是Archimedean Copula函数族(Shiau and Modarres, 2009)。利用样本经验概率和Copula函数联合概率的离差平方和最小可选取最优Copula函数。本文选取Gumbel Copula函数来模拟干旱强度和干旱历时的依存关系。

1.2.3 边缘分布函数

干旱历时和干旱强度分别服从指数分布和Gamma分布，其模型参数可用极大似然估计法得到。干旱历时、干旱强度的联合分布函数和联合超越概率的计算方法详见杨歆雨和李栋梁(2019)。

1.2.4 重现期确定

重现期是用于评价气候事件严重程度的一项重要指标，近年来也受到学者们的广泛关注。干旱重现期定义为在一定年代的资料统计期间内，大于或等于某干旱持续时间(或强度)的干旱出现一次的平均间隔时间，通俗的来讲就是这么严重的干旱，多少年出现一次。Shiau and Shen(2001)推导出了干旱历时和干旱强度单变量重现期公式，但在干旱特征的分析中往往需考虑多个干旱特征量来分析干旱的重现期，因此本文引入干旱历时和干旱强度的联合属性。Shiau(2006)提出基于Copula函数的干旱历时和干旱强度联合分布的重现期可以分为两种情况：干旱历时或干旱强度其中一个达到某一阈值，干旱历时和干旱强度同时达到某一阈值。本文主要分析干旱的联合重现期，即第一种情况，其计算公式如下：

$ T_{\mathrm{DS}}=E(L) /\left[1-F_{\mathrm{DS}}(d, s)\right] $

(3)

式中: T_DS为干旱历时和干旱强度的联合重现期，E(L)为干旱间隔的数学期望值，F_DS(d, s)为干旱历时和干旱强度的联合分布函数。另外，对于建立好的Copula函数，其边际函数概率分布值的不同组合可能产生相同的累计概率，在干旱风险评估中，往往对随机事件的超临界状态更感兴趣，因此引进第二重现期的概念(刘晓云等，2015)。联合重现期与第二重现期不同在于，联合重现期是多元变量的联合属性，第二重现期则是将多元变量投影到一元变量的属性。本文研究发现，由联合重现期计算得到的干旱典型年份多年一遇小于历史记载的干旱重现期，而第二重现期获得的多年一遇接近历史真实情况，因此最终选取第二重现期对福建干旱过程进行描述。

2 干旱气候特征 2.1 干旱过程年季变化特征 2.1.1 全省干旱过程年频数

60年间，共发生69次干旱，年度干旱次数多为1次(共计36年，占52.17%)，一年中发生干旱次数最多为3次的有5年(1964、1971、1977、2007、2011年)，未出现干旱的有9年(1961、1965、1972、1975、1976、1997、2001、2010、2016年)。

2.1.2 全省干旱季节分布

福建气象干旱以单季旱为多，占比达58%，其中又以夏旱为最多，冬旱为次；两季连旱占比达36%，以秋冬旱为多，夏秋旱次之，冬春旱最少；三季连旱占比仅6%，皆为夏秋冬连旱，其他季节组合的三季连旱从未出现过。60年间共发生4起夏秋冬连旱，分别出现在1967、1978、1995和2003年。

2.1.3 历年全省干旱过程强度和持续时间

由于指数累加的干旱强度包含有时间分量，不能充分体现干旱时间长强度弱，或干旱时间短强度强的过程差异，故采用平均强度分析历次全省干旱过程强度演变特征。

干旱平均强度最大的是2011年的春旱(图 1)，其次是1962年的冬春旱，2003年的夏秋冬旱仅位列第三。持续时间最长的是2003年夏秋冬旱，其次是1967年的夏秋冬旱。按照95%、80%、50%百分位数划分强度，特强干旱过程有3次(2011、1962和2003年)，强干旱过程有10次，较强干旱过程有21次，弱干旱过程有35次。同样采用百分位数划分持续时间，特长过程有2003、1967和2020年共3次。多年来，干旱持续时间变化不大(图 1中绿虚线)，但强度变化有所增加(图 1中黑虚线)。

从图 1可见，1962年和2003年属于干旱强度强、持续时间长的类型，2011年春旱属于干旱强度强、持续时间短的类型，2007年春旱属于强度弱、持续时间短的类型，等等。

2.2 干旱空间分布 2.2.1 单站最强干旱过程强度及持续时间空间分布

各站历史最强干旱过程累计MCI指数最大是28.7(平潭，2003年)，最小是14.5(周宁，1967年)，从强度的空间分布(图 2a)看，宁德南部沿海至漳州沿海，以及南平地区闽江流域河谷盆地和九龙江流域河谷盆地强度强，鹫峰山区和闽西强度弱。

最强干旱过程持续时间也是以沿海和武夷山为长(图 2b)，其中宁德沿海大部、福州内陆县市、泉州—厦门—漳州接壤县市、武平县超过220 d，即干旱时长超过7个月，最长的是长泰343 d(1962年7月8日至1963年6月15日)，最短的有119 d(上杭，2019年10月8日至2020年2月3日)。值得一提的是，全省干旱过程最长持续时间224 d(2003年6月19日至2004年1月28日)，远短于单站干旱最长持续时间，对全省而言没有四季连旱，符合气候区概念，但对于局地而言，由于降水年际间、地域间分配不均，局地先后出现极端的四季连旱属正常，但没有出现一年以上的连旱，总有一段时间有足够的降水，通常间断于5—6月的雨季或台风季。

最强干旱发生年份有一定的空间规律(图 2c)，华安、大田、漳平、仙游，以及它们包围的泉州内陆县市(晋江和九龙江流域)、漳州南部发生年份最早(1962年)，其中九龙江流域最强过程具有强度强、持续时间最长的特点；宁德、福州沿海，及内陆大部县市(闽江流域和汀江下游)发生在2002年以后，特别是2003年的典型特旱年在这些区域成片出现。

2.2.2 单站不同干旱等级频数空间分布

统计各站历年MCI指数表征的各类干旱强度(轻旱、中旱、重旱、特旱)年平均次数(图 3)，轻旱和中旱的空间分布相似，频数由中南部沿海向北部山区递减，其中轻旱30.8 d(周宁)~65.7 d(崇武)，中旱18.0 d(周宁)~40.9 d(崇武)；重旱以中南部沿海、三明南部和龙岩东部为多，年平均超过12 d，最多是崇武的15.4 d，最少是寿宁的5.1 d；特旱年平均出现很少，分布散，长汀最少仅为0.1 d，最多不超过5.1 d(平潭和延平)，超过3.5 d的县(市)还有诏安、南靖、永定和晋江。

综上分析，鹫峰山高海拔区域干旱少且轻，中南部沿海县市干旱多且重。

3 多年一遇干旱过程

多年一遇是指小概率事件，在决策服务中更关注某个干旱过程历时和强度是否偏离了常态，出现对自然界有显著不利影响的极端现象。为了进一步认识福建干旱小概率事件的气候特征，统计了全省和66个国家气象站2、5、10、20、50和100年重现期下干旱强度、干旱历时，以及干旱历时和干旱强度联合重现期，各种重现期的阈值可用来衡量干旱过程的极端程度。

3.1 多年一遇计算模型参数确定

由于干旱历时和干旱强度分别服从指数分布和Gamma分布，其中Gamma分布的密度函数是指在序列的有序性、事件发生率的齐次性、计数特征具有独立增量和平稳增量情况下，可以得到事件发生i次时间的概率密度。而指数分布实际上是Gamma分布的特殊情况，Gamma分布可以看成是N个指数分布的数学期望，由此可以认为Gamma分布的模型参数是一种平均态(气候态)下估计出的参数。因此，本文在确定模型参数时，先识别“行标”法中全省干旱序列(即66个站逐日MCI指数逐日均值)的干旱过程，然后用极大似然法估计出气候态下全省干旱过程的模型参数，再利用“本地算法”的全省MCI序列识别用于计算多年一遇干旱过程，采用气候态下的模型参数和最优Copula模型，得到全省干旱过程的重现期。

为了验证确定模型参数方法使用的合理性，对“行标”法和“本地算法”形成的2套历史干旱过程，分别采用Copula函数确定多年一遇模型参数，结果“行标”法历史干旱过程多年一遇与史料记载的严重程度很不一致，多年一遇的年数皆偏小。究其原因，表示平均态的“行标”过程弱化了对福建有影响的干旱过程，而以“本地算法”过程估计参数是将有影响的干旱过程看做气候常态了，所以最终选择的方法明显提高了多年一遇的量值，特别是2003年福建百年一遇的特旱年，多年一遇值更接近历史。

3.2 区域干旱过程的重现期特征 3.2.1 全省干旱过程重现期时序特征

福建省严重干旱的年份较少，超50年一遇的干旱仅有1967年夏秋冬旱和2003年夏秋冬旱；超30年一遇的干旱为1962年冬春旱和1995年夏秋冬旱；超20年一遇的干旱为2020年夏秋旱；超10年一遇的干旱为1991年春夏旱、2002年春旱、2018年春夏旱及2019年秋冬旱；其余年份干旱重现期基本在1~3年(图 4)。

3.2.2 不同重现期下的干旱强度空间分布

根据单站干旱过程强度指数，分别统计了单站干旱过程强度的多年一遇阈值，如图 5所示，无论多少年一遇，鹫峰山区都是最弱区；对于2 a、5 a、10 a一遇的干旱，其强度的空间分布较为相似，中南部沿海一带为干旱强区；对于20 a、50 a、100 a一遇的干旱，除了中南部沿海地区的干旱强区外，南平西北部(光泽、邵武、武夷山)也是一个强区。50年一遇强区干旱强度指数超过18，最强的东山达到20.1；百年一遇强区干旱强度指数超过20，最强的仍然是东山，为21.9。

3.2.3 不同重现期下的干旱历时空间分布

同样的，根据单站干旱过程历时，统计了单站干旱历时多年一遇阈值，干旱历时的阈值随着年数的增加，空间分布与强度阈值类似(图 6)，鹫峰山区属于干旱历时短的区域，中南部沿海区域历时长。分析图 6a，平均2年就会遇到历时40~80 d的持续气象干旱；50年一遇除了鹫峰山区，基本上都在200 d以上，持续气象干旱接近4个月(图 6e)；而百年一遇全省都达到200 d以上，其中沿海区域接近或超过300 d，历时近10个月(图 6f)。

结合图 5可以看出，对于2 a和5 a这类重现期较短(频发)的干旱，中南部沿海一带有干旱历时长、干旱强度强的特点，越往北干旱历时短、干旱强度弱；但对于重现期长，特别是50年和100年一遇的干旱，中南部沿海一带仍有干旱历时长、干旱强度强的特点，另一个强度较强、历时较长的区域在南平西北部(光泽、邵武、武夷山)；鹫峰山区属于干旱历时最短、干旱强度最弱的区域。

3.3 典型干旱年的重现期特征

为了分析1961年以来各年干旱的空间分布，将单站年内间隔小于1个月的干旱过程合并，重新统计干旱强度、历时和重现期，各站干旱起止时间差异很大，允许跨年统计。

首先分析1962、1967、1971、2003、2004年和2020年共6年，它们的共同特点是干旱强度之强(特强等级)、范围之广(大半个省域)、持续时间之长(特长等级)为历史罕见(图 7)，史料记载的灾害严重的区域都发生了超50年一遇的干旱，局部超百年，而其余年份皆未出现超百年一遇的干旱。据史料记载，1962年冬春连旱，加上1963年秋旱，中南部沿海成为全省最严重的旱区，九龙江、晋江流域发生了50年至超100年一遇的干旱，可以徒步过河。1967年夏旱连秋旱，沿海和龙岩、三明干旱严重，其中闽江口周围和内陆局部发生了50年至超100年一遇的干旱。1971年春旱、夏旱，接着又出现秋冬旱，严重干旱区主要分布在内陆，特别是南平地区，闽江中游流域发生了50年至超100年一遇的干旱，古田水库无水发电，福州照明用电受到影响。2003年夏秋冬连旱是突破历史纪录的特旱年，宁德地区受灾最为严重，中北部沿海和三明、南平局部发生了50年至超100年一遇的干旱。续2003年冬季结束(跨年至2004年2月)的特旱年后，接着2004年发生的春夏旱，局部区域水资源短缺导致严重的饮用水紧张，龙岩西部出现了超百年一遇的干旱。2020年夏秋旱主要发生在沿海和龙岩南部，霞浦、崇武、南安、厦门、同安和永定发生超百年一遇的气象干旱。

除了上述大范围的典型干旱过程外，1963、1964、1977、1978、1979、1980、1991、1993、1996、2002、2007年和2011年共12年局部出现了超50年一遇干旱，其中1991年和2011年干旱范围虽然大，但是强度和历时明显较典型特旱年逊色(图 8，仅列出1991年和2011年)，这些年份受旱面积基本上都在200万亩(1亩=666.67 m²)以上。有趣的是，1977—1980年连续4年在不同的局地发生小范围严重的干旱：严重旱区1977年在福州和漳州；1978年在三明和南平，以及中部沿海；1979年在福州和宁德，以及龙岩局部；1980年在福州及漳州局部，持续4年少雨导致1980年水资源特别紧张，全省最大受旱面积高达548.14万亩，粮食减产6亿kg。除以上列出的18年，加上9年无旱年外，余下的33年干旱重现期均未超过50年。

4 结论与讨论 4.1 结论

本文基于MCI干旱指数，利用“行标”法与福建“本地算法”，识别了1961—2020年福建省66个国家气象站及全省区域的干旱过程，分析了干旱的变化特征及其重现规律，主要结论如下：

(1) 福建省基本每年都会发生一次干旱，最多的可达三次；福建主要受单季干旱和两季连旱影响，两季连旱主要为夏秋旱和秋冬旱，全省区域无四季连旱；干旱历时和干旱强度的变化趋势基本一致；

(2) 福建大部分干旱过程存在约为1~3 a的周期，严重干旱较少，超50年一遇的干旱仅有1967年夏秋冬旱和2003年夏秋冬旱；超30年一遇的干旱为1962年冬春旱和1995年夏秋冬旱；超20年一遇的干旱为2020年夏秋旱；超10年一遇的干旱为1991年春夏旱、2002年春旱、2018年春夏旱及2019年秋冬旱；

(3) 对于重现期较短(频发)的干旱，中南部沿海一带有干旱历时长、干旱强度强的特点，越往北干旱历时短、干旱强度弱；但对于重现期长，特别是50年和100年一遇的干旱，中南部沿海一带仍有干旱历时长、干旱强度强的特点，另一个强度较强、历时较长的区域在南平西北部(光泽、邵武、武夷山)；鹫峰山区属于干旱历时最短、干旱强度最弱的区域；

(4) 1962、1967、1971、2003、2004年和2020年干旱强度之强、范围之广、持续时间之长为历史罕见，灾害严重的区域都发生了超50年一遇的干旱，局部超百年；1963、1964、1977、1978、1979、1980、1991、1993、1996、2002、2007年和2011年共12年局部出现了超50年一遇干旱；除上述18年及9年无旱年之外，余下的33年干旱重现期均未超过50年。

4.2 讨论

福建是一个干旱比较频发的省份，每年各地或多或少都会发生气象干旱，这种干旱气候格局已经造就了当地人类和生态景观的布局，对于一般的干旱灾害已具有防御能力。由于每年干旱空间差异非常大，导致针对性干旱防御策略存在很大的不确定性。本文仅对福建的干旱分布特征等做初步的分析与讨论，为干旱灾害风险评估奠定基础。实际上干旱灾害程度与承灾体、孕灾环境及天气系统息息相关，“沿海干旱强(多)于内陆”这一气候现象，究其原因，福建年雨量多呈双峰型，雨季主要为锋面系统性降水，集中在内陆地区；夏季主要是台风雨，集中在沿海地区。而沿海地区下垫面多为岩土砂石，植被相对内陆稀少，地下储水能力差，加上台风年际变化大，极易形成气象干旱。内陆植被覆盖率高，地下水蕴含量丰富，雨季降水相对台风雨充沛，因此气象干旱相对沿海弱。另外，干旱灾害的防御能力，还要考虑影响区域人口、经济状况，比如：人口、工业密集地区，对于用水的需求量大，干旱对其的影响相对较大；经济较好的区域对于干旱的防御水平更高，水库较多的地区对于干旱时用水的调节能力也更强等等，仍需要开展深入的探索。目前已分析了福建干旱的致灾因子危险性，今后将结合福建地形、环流背景等开展进一步的灾害风险评估研究。