Exploring early warning and forecasting of meteorological risk of landslide and rockfall induced by meteorological factors by the approach of machine learning
-
摘要: 在划分气象风险等级时,传统地质灾害气象风险预警方法忽略了承灾体脆弱性因素,且气象风险预报等级整体偏高,导致高等级风险区空报率较高。基于此,提出基于机器学习的滑坡、崩塌灾害气象风险预警方法。利用信息量法,分析气象因素影响程度。选取坐标点、降雨量、易发生等级,将其作为机器学习人工神经网络的输入节点,判断是否发生崩塌、滑坡灾害;针对地质灾害区域,根据影响程度计算气象引发因子指数,结合滑坡、崩塌灾害潜势度G和承灾体脆弱性M,确定气象风险预警指数R,划分预警级别,完成滑坡、崩塌灾害气象风险预警。实验结果表明,设计方法有效降低了三级预报和四级预警空报率,提升了预警精细化程度。
-
关键词:
- 机器学习 /
- 地质灾害;滑坡灾害;崩塌灾害 /
- 气象风险 /
- 预警
Abstract: In the traditional methods of meteorological risk early warning and forecasting, the vulnerability factors of disaster bearing bodies are ignored when classifying the meteorological risk level, and the meteorological risk prediction level is relatively high, which leads to the high air report rate in high-level risk areas. Based on this, a method of early warning and forecasting of meteorological risk of landslide and collapse geological disasters based on machine learning is proposed. By using the information quantity method, the influence degree of meteorological factors is analyzed, and the coordinate point, rainfall and prone level are selected as input nodes of machine learning artificial neural network to judge whether geological disaster occurs; for the area of ground damage, the meteorological cause sub index is calculated according to the influence degree. Combined with the potential degree of geological disaster and vulnerability of disaster bearing body, the meteorological risk warning index is determined, divide the warning and forecast level, and complete the forecast of geological disaster meteorological risk. The experimental results show that the proposed method can effectively reduce the three-level forecast air report rate and the fourth level air alarm rate, and improve the precision of the early warning forecast. -
0. 引言
滑坡、崩塌灾害是地质灾害中的重要灾种,对社会经济、生命财产、地理环境都产生了重要影响[1]。贵州是地质灾害高发易发地区,常见诱发因素主要为降雨[2]。因此,研究地质灾害气象风险预警方法,对地质灾害防治具有重要作用。研究者对该领域进行了很多研究。文献[3]应用GIS数据处理功能,分析地质灾害区域的工程地质,针对复杂地质的地理环境,绘制出各气象风险等级下区域地质灾害分布图,但该方法对气象因素提取不全面,造成气象风险预警等级整体偏高。文献[4]分析气象环境与降水量之间的规律性,包括降雨阈值和临界降雨量等,通过地质灾害调查统计,综合评价气象风险,但该方法对降雨量的相关性分析较差,划分的气象风险预警等级同样较高。
针对以上问题,提出基于机器学习的滑坡、崩塌灾害气象风险预警方法。借助机器学习算法中人工神经网络实现贵州省滑坡、崩塌灾害气象风险预警方法的研究。实验结果表明:采用设计方法有效降低了滑坡、崩塌三级、四级预警空报率,提升了预警精细化程度。
1. 贵州省滑坡、崩塌灾害气象风险预警方法设计
1.1 气象因素对滑坡、崩塌灾害影响程度计算
为实现贵州省滑坡、崩塌灾害气象风险预警,需要提取贵州省滑坡、崩塌灾害相关气象因素,计算气象因素对滑坡、崩塌灾害的影响程度。汇总贵州省各区域的地质灾害详细调查报告、气象局实测降雨量数据,采用资料查阅、实地访问调查等方式,以月报形式统计气象引发灾害事件,同时综合考虑GPS、天气雷达、闪电定位、自动雨量站等数据来源,采集非常规观测资料,获取灾害隐患点数据和灾害易发分区数据。筛选与滑坡、崩塌灾害相关的气象因素。将滑坡、崩塌灾害的发生看作气象因素和下垫面相互作用结果。其中,气象因素主要为降水,包括降雨诱发作用、降雨滞后作用等因素。利用信息量法,客观反映预警单元的地质条件,分析气象因素对滑坡、崩塌灾害提供的信息量[5]。计算气象因素和地质灾害的相关函数
,公式为:(1) 式中:
——贵州省地质灾害; ——第 种地质灾害相关气象因素; ——因素数量。 是与 有关的变量, 为各变量 条件下的条件概率[6]。计算单个气象因素 的信息量值 ,公式为:(2) 将整个贵州省区域划分为单元网格,利用频率统计,估算条件概率
,确定省内地质灾害敏感性[7]。则地质灾害相关气象因素的总信息量值 ,计算公式为:(3) 式中:
——第 种气象因素引发的地质灾害点数量; ——地质灾害点总数; ——第 种气象因素引发的地质灾害点面积; ——贵州省内总面积。通过总信息量值,定量分析气象因素对地质灾害的影响程度,
值为正时,判定气象因素利于地质灾害发生, 值为负时,则判定不利于灾害发生,且气象因素影响程度大小与 值大小呈正相关[8]。至此完成气象因素对地质灾害影响程度的计算。1.2 地质灾害发生的判断
在确定气象因素对地质灾害影响程度基础上,利用机器学习中的人工神经网络,判断贵州省各区域是否发生地质灾害。首先,参照采集地质环境数据,结合岩性特征、地形地貌、以及气候条件因素,对贵州省各区域进行条件比较,利用GIS空间分析功能,定量划分贵州省地质灾害易发生区,并明确其易发生等级[9]。划分结果见图1。
根据贵州省近20年历史降雨量数据,以及记录的滑坡、崩塌灾害数据,明确各区域的降雨量数据,包含当日临界雨量和5日临界雨量。同时,根据滑坡、崩塌灾害野外调查结果,对临界雨量值进行适当调整,以此减小历史统计数据误差[10]。最后,将预报区域中易发生区等级、实际降雨量、坐标点X坐标和Y坐标,作为机器学习的人工神经网络4个输入节点,利用线性函数,激活神经网络的输入层和输出层,再利用Sigmoid函数,激活隐含层,输入前对数据进行归一化处理,使各数据处于同等水平,消除量纲影响,避开Sigmoid函数训练数据的饱和区[11]。神经网络结构见图2。
图2中,设置网络误差收敛到最小时,其相对应的隐含层节点数为4,将滑坡、崩塌灾害性发生可能性,作为神经网络的1个输出节点。根据机器学习输出结果,判定坐标点区域是否发生滑坡、崩塌灾害,完成滑坡、崩塌灾害发生的判断。
1.3 滑坡、崩塌灾害气象风险等级划分
针对贵州省滑坡、崩塌灾害发生区域,根据气象因素影响程度,计算气象引发因子指数,结合该区域的地质灾害潜势度、承灾体脆弱性,划分气象风险的预警等级。
对预警区域进行单元编号,计算单元区域
内的有效降雨量 ,公式为:(4) 式中:
——有效降雨日数; ——预警当日向前计算的天数; ——区域 的当日和5日预报雨量值; ——当日和5日临界雨量值;计算单元区域
的气象引发因子指数 ,公式为:(5) 式中:
——有效降雨系数; ——代表单元区域滑坡、崩塌灾害气象因素分量。计算滑坡、崩塌灾害潜势度
,公式为:(6) 式中:
——地质环境因子总个数; ——地质环境因子 的权重; ——地质环境因子 的量化值。采用评价指标方式,在承灾体范围内,提取一级指标和二级指标,计算承灾体脆弱性
,公式为:(7) 式中:
——评价因子总个数; ——评价因子 的权重; ——脆弱性评价因子 的量化值。气象风险可概化公式为:
(8) 其中,
为区域 的气象风险预警指数,取值介于0~1之间,预警指数越大,判定其地质灾害越易发生[14-17]。表 1 滑坡、崩塌灾害高易发区气象风险预警级别Table 1. Early warning level of meteorological risk in high areas prone to geological disasters累积降水
/mm预报小雨
0.01~10预报中雨
10~25预报大雨
25~50预报暴雨
50~100预报大暴雨
≥100≤30 蓝色 黄色 橙色 红色 30~50 蓝色 黄色 橙色 红色 红色 50~100 黄色 橙色 红色 红色 红色 ≥100 橙色 红色 红色 红色 红色 表 2 滑坡、崩塌灾害中易发区气象风险预警级别Table 2. Warning level of meteorological risk in areas prone to geological disasters累积降水
/mm预报小雨
0.01~10预报中雨
10~25预报大雨
25~50预报暴雨
50~100预报大暴雨
≥100≤30 蓝色 黄色 橙色 30~50 蓝色 黄色 橙色 红色 50~100 蓝色 黄色 橙色 红色 红色 ≥100 黄色 橙色 红色 红色 红色 表 3 滑坡、崩塌灾害低易发区气象风险预警级别Table 3. Early warning level of meteorological risk in low areas prone to geological disasters累积降水
/mm预报小雨
0.01~10预报中雨
10~25预报大雨
25~50预报暴雨
50~100预报大暴雨
≥100≤30 蓝色 黄色 30~50 蓝色 黄色 橙色 50~100 蓝色 黄色 橙色 红色 ≥100 蓝色 黄色 橙色 红色 红色 注:其中预报降水为24 h预报降雨量,累积降水为最近五天累计降雨量。 在滑坡、崩塌灾害气象风险预警级别中,滑坡、崩塌灾害气象风险预警的等级为:
(1)蓝色预警(一级):有一定风险,关注降雨;
(2)黄色预警(二级):风险较高,关注降雨,做好监控;
(3)橙色预警(三级):风险高,注意降雨,做好监控及应急准备;
(4)红色预警(四级):风险很高,注意降雨,做好监控与应急撤离准备。
将气象风险预警指数
,与预警临界值相比较,确定该区域是否发布预警,以及相应的预警级别,完成基于机器学习的贵州省滑坡、崩塌灾害气象风险预警方法设计。2. 实验分析
选取两种常规滑坡、崩塌灾害气象风险预警方法,与此次设计方法进行对比实验,比较各预警等级的空报率大小。
2.1 实验准备
将滑坡、崩塌灾害,作为贵州省地质灾害研究范围,采集降雨量数据和地质灾害数据,作为实验数据源,在样本中剔除不符合降雨诱发地质灾害个例、以及不匹配区域站降雨资料的降雨量数据。其降雨量历史信息见图3。
各区域的当日临界雨量和5日临界雨量,其各级预警的具体数值见表4。
表 4 贵州省当日临界雨量和5日临界雨量Table 4. Critical rainfall and mm rainfall of 5 th Day of Guizhou Province灾害易发区域 一级 二级 三级 四级 当日临界雨量
/m不易发区 92 55 37 28 低易发区 110 67 45 34 中易发区 132 79 53 40 高易发区 254 151 101 76 5日临界雨量
/m不易发区 223 133 89 67 低易发区 243 157 103 77 中易发区 262 157 105 79 高易发区 304 181 121 91 统计可得2014—2020年之间,贵州省地质灾害共发生1204处,发生地质灾害具体数据见表5。
表 5 贵州省典型地质灾害统计数据Table 5. Statistical data of typical geological disasters in Guizhou Province灾害点类型 灾害点数量/个 分布市镇数量/个 占灾害点总数比例/% 滑坡 1032 58 85.7% 崩塌 111 19 9.2% 泥石流 29 12 2.4% 地面塌陷 25 8 2.1% 地裂缝 7 3 0.5% 可见滑坡、崩塌占灾害总数的94.9%,三组预警方法分别根据以上历史数据中的滑坡、崩塌灾害,对贵州省地质灾害气象风险进行预警,并以2020年地质灾害作为参照,对比检验三组预警结果。
2.2 实验结果
2.2.1 崩塌灾害空报率测试结果
2020年崩塌灾害隐患点数量共29处,三组方法均可准确预测出该类地质灾害,其预警等级见图4。
由图4可知,两组常规方法三级预报数量和四级预警数量要明显多于设计方法,隐患点崩塌预警的严重程度整体偏高。进一步统计所有年份中,各预警级别的空报率,实验对比结果见表6。
表 6 崩塌预警空报率Table 6. Empty reporting rate of collapse early warning and forecast设计方法 常规方法1 常规方法2 一级预报/% 0 0 0 二级预报/% 0 0 0 三级预报/% 8.27 14.92 17.92 四级预警/% 7.26 13.29 19.26 由表6可知,相比常规方法1和常规方法2,设计方法对崩塌的三级预报空报率分别降低了6.65%和9.65%,四级预警空报率分别降低了6.03%和12.0%。
2.2.2 滑坡灾害空报率测试结果
2020年滑坡灾害隐患点数量共33处,三组方法都准确预测出该类地质灾害,其预警等级见图5。
由图5可知,针对滑坡这一地质灾害,两组常规方法的三级预警数量和四级预警数量,同样多于设计方法,隐患点预警的严重程度仍整体偏高。进一步统计所有年份中,各预警级别的空报率,实验对比结果见表7。
表 7 滑坡预警空报率Table 7. Empty reporting rate of landslide early warning and forecast设计方法 常规方法1 常规方法2 一级预报/% 0 0 0 二级预报/% 0 0 1.21 三级预报/% 9.92 14.96 16.92 四级预警/% 6.12 14.63 17.29 由表7可知,相比常规方法1和常规方法2,设计方法对滑坡的三级预警空报率分别降低了5.04%和7%,四级预警空报率分别降低了8.51%和11.17%,且常规方法2的二级预警仍存在空报率。
3. 结论
针对现有滑坡、崩塌地质灾害预警方法中存在的不足,本文提出采用机器学习算法对地质灾害气象风险进行预警的方法。
(1)该方法通过采用机器学习神经网络中节点的输入,有效预测地质灾害发生。
(2)采用设计方法与常规方法对比中,本文方法对崩塌预警的三级预报空报率分别降低了5.04%和7%,四级预警空报率分别降低了8.51%和11.17%;对滑坡预警的空报率中,三级预报空报率分别降低了5.04%和7%,四级预警空报率分别降低了8.51%和11.17%,验证了本文方法适用于贵州省滑坡、崩塌灾害气象风险预警。
但此次研究仍存在一定不足,在今后研究中,会持续验证每年实际发生的地质灾害,修正该模型,进一步提高地质灾害预警精度。
-
表 1 滑坡、崩塌灾害高易发区气象风险预警级别
Table 1 Early warning level of meteorological risk in high areas prone to geological disasters
累积降水
/mm预报小雨
0.01~10预报中雨
10~25预报大雨
25~50预报暴雨
50~100预报大暴雨
≥100≤30 蓝色 黄色 橙色 红色 30~50 蓝色 黄色 橙色 红色 红色 50~100 黄色 橙色 红色 红色 红色 ≥100 橙色 红色 红色 红色 红色 表 2 滑坡、崩塌灾害中易发区气象风险预警级别
Table 2 Warning level of meteorological risk in areas prone to geological disasters
累积降水
/mm预报小雨
0.01~10预报中雨
10~25预报大雨
25~50预报暴雨
50~100预报大暴雨
≥100≤30 蓝色 黄色 橙色 30~50 蓝色 黄色 橙色 红色 50~100 蓝色 黄色 橙色 红色 红色 ≥100 黄色 橙色 红色 红色 红色 表 3 滑坡、崩塌灾害低易发区气象风险预警级别
Table 3 Early warning level of meteorological risk in low areas prone to geological disasters
累积降水
/mm预报小雨
0.01~10预报中雨
10~25预报大雨
25~50预报暴雨
50~100预报大暴雨
≥100≤30 蓝色 黄色 30~50 蓝色 黄色 橙色 50~100 蓝色 黄色 橙色 红色 ≥100 蓝色 黄色 橙色 红色 红色 注:其中预报降水为24 h预报降雨量,累积降水为最近五天累计降雨量。 表 4 贵州省当日临界雨量和5日临界雨量
Table 4 Critical rainfall and mm rainfall of 5 th Day of Guizhou Province
灾害易发区域 一级 二级 三级 四级 当日临界雨量
/m不易发区 92 55 37 28 低易发区 110 67 45 34 中易发区 132 79 53 40 高易发区 254 151 101 76 5日临界雨量
/m不易发区 223 133 89 67 低易发区 243 157 103 77 中易发区 262 157 105 79 高易发区 304 181 121 91 表 5 贵州省典型地质灾害统计数据
Table 5 Statistical data of typical geological disasters in Guizhou Province
灾害点类型 灾害点数量/个 分布市镇数量/个 占灾害点总数比例/% 滑坡 1032 58 85.7% 崩塌 111 19 9.2% 泥石流 29 12 2.4% 地面塌陷 25 8 2.1% 地裂缝 7 3 0.5% 表 6 崩塌预警空报率
Table 6 Empty reporting rate of collapse early warning and forecast
设计方法 常规方法1 常规方法2 一级预报/% 0 0 0 二级预报/% 0 0 0 三级预报/% 8.27 14.92 17.92 四级预警/% 7.26 13.29 19.26 表 7 滑坡预警空报率
Table 7 Empty reporting rate of landslide early warning and forecast
设计方法 常规方法1 常规方法2 一级预报/% 0 0 0 二级预报/% 0 0 1.21 三级预报/% 9.92 14.96 16.92 四级预警/% 6.12 14.63 17.29 -
[1] 申欣凯, 吕义清, 张静, 等. 基于ArcGIS的山西省暴雨灾害风险评估[J]. 国土资源科技管理,2020,37(1):61 − 73. [SHEN Xinkai, LYU Yiqing, ZHANG Jing, et al. Risk assessment of rainstorm disaster in Shanxi Province based on ArcGIS[J]. Scientific and Technological Management of Land and Resources,2020,37(1):61 − 73. (in Chinese with English abstract) [2] 罗贵东, 于竹娟, 张驹. 基于概率分布的达州市地质灾害发生频次与降水量级的关系[J]. 成都信息工程大学学报,2020,35(6):697 − 703. [LUO Guidong, YU Zujuan, ZHANG Ju. The relationship between the frequency of geological disasters and the precipitation in Dazhou based on probability distribution[J]. Journal of Chengdu University of Information Technology,2020,35(6):697 − 703. (in Chinese with English abstract) [3] 李宇梅, 杨寅, 狄靖月, 等. 全国地质灾害气象风险精细化网格预报方法及其应用[J]. 气象,2020,46(10):1310 − 1319. [LI Yumei, YANG Yin, DI Jingyue, et al. Meteorological risk assessment method of geological disaster in China and its mesh refinement application[J]. Meteorological Monthly,2020,46(10):1310 − 1319. (in Chinese with English abstract) DOI: 10.7519/j.issn.1000-0526.2020.10.005 [4] 兰腾达. 宁德市地质灾害气象风险预警系统应用与展望[J]. 福建地质,2020,39(2):155 − 160. [LAN Tengda. Application and prospect of early warning of meteorological risk system for geological hazards, Ningde City[J]. Geology of Fujian,2020,39(2):155 − 160. (in Chinese with English abstract) DOI: 10.3969/j.issn.1001-3970.2020.02.008 [5] 薛俊卓, 詹辉, 张国良, 等. 长输油气管道河沟段水毁灾害气象风险评价研究[J]. 安全与环境工程,2020,27(2):168 − 174. [XUE Junzhuo, ZHAN Hui, ZHANG Guoliang, et al. Meteorological risk assessment of river channel washout disaster of long-distance oil and gas pipelines[J]. Safety and Environmental Engineering,2020,27(2):168 − 174. (in Chinese with English abstract) [6] 黎力, 朱永亮, 陈亮, 等. 重庆市地质灾害气象风险预警问题研究[J]. 安徽农业科学,2019,47(24):234 − 236. [LI Li, ZHU Yongliang, CHEN Liang, et al. Study on geological hazard meteorological risk early warning in Chongqing[J]. Journal of Anhui Agricultural Sciences,2019,47(24):234 − 236. (in Chinese with English abstract) DOI: 10.3969/j.issn.0517-6611.2019.24.069 [7] 罗鸿东, 李瑞冬, 张勃, 等. 基于信息量法的地质灾害气象风险预警模型:以甘肃省陇南地区为例[J]. 地学前缘,2019,26(6):289 − 297. [LUO Hongdong, LI Ruidong, ZHANG Bo, et al. An early warning model system for predicting meteorological risk associated with geological disasters in the Longnan area, Gansu Province based on the information value method[J]. Earth Science Frontiers,2019,26(6):289 − 297. (in Chinese with English abstract) [8] 杨竹云, 李华宏, 胡娟, 等. 昭通市地质灾害特征及气象风险预警模型研究[J]. 云南大学学报(自然科学版),2019,41(4):753 − 764. [YANG Zuyun, LI Huahong, HU Juan, et al. Research on features of geological hazard and meteorological risk early-warning model in Zhaotong[J]. Journal of Yunnan University(Natural Sciences Edition),2019,41(4):753 − 764. (in Chinese with English abstract) [9] 方然可, 刘艳辉, 苏永超, 等. 基于逻辑回归的四川青川县区域滑坡灾害预警模型[J]. 水文地质工程地质,2021,48(1):181 − 187. [FANG Ranke, LIU Yanhui, SU Yongchao, et al. A early warning model of regional landslide in Qingchuan County, Sichuan Province based on logistic regression[J]. Hydrogeology & Engineering Geology,2021,48(1):181 − 187. (in Chinese with English abstract) [10] 杨寅, 包红军, 彭涛. 台风“鲇鱼”强降水引发的地质灾害气象风险预警检验与分析[J]. 暴雨灾害,2019,38(3):221 − 228. [YANG Yin, BAO Hongjun, PENG Tao. Verification and analysis of meteorological early warning of geological hazards during precipitation of Typhoon "MEGI"[J]. Torrential Rain and Disasters,2019,38(3):221 − 228. (in Chinese with English abstract) [11] 包红军, 曹勇, 林建, 等. 山洪灾害气象预警业务技术进展[J]. 中国防汛抗旱,2020,30(1):40 − 47. [BAO Hongjun, CAO Yong, LIN Jian, et al. Operational technology advances in meteorological early warning for flash flood disasters[J]. China Flood & Drought Management,2020,30(1):40 − 47. (in Chinese with English abstract) [12] 方楠, 黄清瀚, 丁雨鑫, 等. 浙江省近3 a气象类灾害预警信息特征[J]. 浙江农业科学,2020,61(6):1246 − 1250. [FANG Nan, HUANG Qinghan, DING Yuxin, et al. The characteristics of meteorological disaster early warning information in Zhejiang Province in recent 3 years[J]. Journal of Zhejiang Agricultural Sciences,2020,61(6):1246 − 1250. (in Chinese with English abstract) [13] 高树静, 董廷坤, 王程龙. 基于机器学习的车道线检测系统仿真与优化[J]. 计算机仿真,2020,37(2):140 − 143. [GAO Shujing, DONG Tingkun, WANG Chenglong. Simulation and optimization of lane detection system based on machine learning[J]. Computer Simulation,2020,37(2):140 − 143. (in Chinese with English abstract) DOI: 10.3969/j.issn.1006-9348.2020.02.029 [14] 李爽爽. 精细化暴雨监测预报及风险预警系统的开发与应用[J]. 浙江气象,2020,41(3):29 − 35. [LI Shuangshuang. Development and application of refined rainstorm monitoring and forecasting and risk early warning system[J]. Journal of Zhejiang Meteorology,2020,41(3):29 − 35. (in Chinese with English abstract) [15] 刘云, 康卉君. 江西崩塌滑坡泥石流灾害空间时间分布特征分析[J]. 中国地质灾害与防治学报,2020,31(4):107 − 112. [LIU Yun, KANG Huijun. Spatial-temporal distribution of landslide, rockfall and debris flow hazards in Jiangxi Province[J]. The Chinese Journal of Geological Hazard and Control,2020,31(4):107 − 112. (in Chinese with English abstract) [16] 程素珍, 路璐, 翟淑花, 等. 2004—2018年北京市突发地质灾害时空分布特点和监测预警状况[J]. 中国地质灾害与防治学报,2020,31(6):38 − 46. [CHENG Suzhen, LU Lu, ZHAI Shuhua, et al. Temporal-spatial distribution and monitoring and early warning of sudden geological disasters in Beijing during the period of 2004 to 2018[J]. The Chinese Journal of Geological Hazard and Control,2020,31(6):38 − 46. (in Chinese with English abstract) [17] 李滨, 殷跃平, 高杨, 等. 西南岩溶山区大型崩滑灾害研究的关键问题[J]. 水文地质工程地质,2020,27(4):5 − 13. [LI Bin, YIN Yueping, GAO Yang, et al. Critical issues in rock avalanches in the karst mountain areas of southwest China[J]. Hydrogeology & Engineering Geology,2020,27(4):5 − 13. (in Chinese with English abstract) -
期刊类型引用(1)
1. 苏娜,徐林荣,李永威,李佳珅,汤玉兰. 汶川地震震后七盘沟泥石流动态物源危险性评价. 中国地质灾害与防治学报. 2025(01): 16-27 . 本站查看
其他类型引用(0)