CN116720388B - 森林草原火场蔓延预测的方法、装置、设备和介质 - Google Patents
森林草原火场蔓延预测的方法、装置、设备和介质 Download PDFInfo
- Publication number
- CN116720388B CN116720388B CN202310994183.2A CN202310994183A CN116720388B CN 116720388 B CN116720388 B CN 116720388B CN 202310994183 A CN202310994183 A CN 202310994183A CN 116720388 B CN116720388 B CN 116720388B
- Authority
- CN
- China
- Prior art keywords
- fire
- spreading
- correction factor
- wind
- slope
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 230000007480 spreading Effects 0.000 claims abstract description 202
- 241000016949 Acalypha chamaedrifolia Species 0.000 claims abstract description 29
- 239000013598 vector Substances 0.000 claims abstract description 17
- 238000012545 processing Methods 0.000 claims abstract description 13
- 238000012937 correction Methods 0.000 claims description 73
- 238000004364 calculation method Methods 0.000 claims description 24
- 238000004088 simulation Methods 0.000 claims description 12
- 238000007781 pre-processing Methods 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims 1
- 238000003786 synthesis reaction Methods 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 6
- 238000004891 communication Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000012952 Resampling Methods 0.000 description 3
- 238000002485 combustion reaction Methods 0.000 description 3
- 230000001174 ascending effect Effects 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
- G06Q50/265—Personal security, identity or safety
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Tourism & Hospitality (AREA)
- Data Mining & Analysis (AREA)
- Development Economics (AREA)
- General Engineering & Computer Science (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Mathematical Physics (AREA)
- Operations Research (AREA)
- Primary Health Care (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Geometry (AREA)
- Computer Security & Cryptography (AREA)
- Educational Administration (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Entrepreneurship & Innovation (AREA)
- Algebra (AREA)
- Computational Mathematics (AREA)
- Evolutionary Computation (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computer Hardware Design (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本申请提供了一种森林草原火场蔓延预测的方法、装置、设备和介质,涉及森林草原火场预测技术领域,该方法包括:获取待预测火点信息,并对待预测火点信息与地形数据、植被数据、气象预报数据进行分辨率统一处理;根据待预测火点信息和时间信息,通过气象预报数据和地形数据获取对应时刻、对应位置的风向和坡向,并基于风向和坡向合成矢量确定火场蔓延方向;通过改进的王正非模型计算火场蔓延速度;基于火场蔓延速度计算预设时间步长内的火头、火翼和火尾对应的火场蔓延距离,并基于火场蔓延距离计算蔓延边界。本申请提高了森林草原火场蔓延预测的准确性。
Description
技术领域
本申请涉及森林草原火场预测技术领域,尤其是涉及一种森林草原火场蔓延预测的方法、装置、设备和介质。
背景技术
火灾具有突发性、灾害发生的随机性、短时间内能造成巨大损失的特点,尤其针对森林草原火灾,准确的预测火灾蔓延情况可以尽可能的减少损失。相关技术中,采用森林草原火灾蔓延模型的方式需要大量的火场信息,导致预测难度较大;并且,目前的森林草原火灾蔓延模型在进行火灾蔓延预测时,适用范围有限,无法适应复杂的蔓延场景,火灾蔓延预测的准确性不高。
发明内容
本申请的目的在于提供一种森林草原火场蔓延预测的方法、装置、设备和介质,提高了森林草原火场蔓延预测的准确性。
第一方面,本申请实施例提供了一种森林草原火场蔓延预测的方法,方法包括:获取待预测火点信息,并对待预测火点信息与地形数据、植被数据、气象预报数据进行分辨率统一处理;根据待预测火点信息和时间信息,通过气象预报数据和地形数据获取对应时刻、对应位置的风向和坡向,并基于风向和坡向合成矢量确定火场蔓延方向;通过改进的王正非模型计算火场蔓延速度;基于火场蔓延速度计算预设时间步长内的火头、火翼和火尾对应的火场蔓延距离,并基于火场蔓延距离计算蔓延边界。
可选的,基于风向和坡向合成矢量确定火场蔓延方向,包括:
其中,spread head_dir为火头蔓延方向;wind dir为风向;Aspect为下垫面坡向。
可选的,通过改进的王正非模型计算火场蔓延速度,包括:
V real=V×F burn ×F wind ×F slope ×F cover ×F tem ×F rh;
其中,V real为火头蔓延速度(m/min);V为蔓延初速度(m/min);F burn为可燃物类型校正因子,F wind为风力校正因子,F slope为地形坡度校正因子,F cover为下垫面类型校正因子,F tem为气温校正因子,F rh为湿度校正因子。
可选的,改进的王正非模型包含的气温校正因子和相对湿度校正因子如下:
气温校正因子:
其中,为气温等级划分,/>为各气温等级对应的系数因子,气温等级越高数值越大;
相对湿度校正因子:
其中,为相对湿度等级划分,/>为等级对应的系数因子,相对湿度等级越高数值越小。
可选的,方法还包括,基于如下公式分别确定火头蔓延速度、火翼蔓延速度和火尾蔓延速度:
其中,Vhead为火头蔓延速度(m/min),Vwing为火翼蔓延速度(m/min),Vrear为火尾蔓延速度(m/min),m、n为不同风力等级对应的系数因子,风力等级越高数值越小。
可选的,基于火场蔓延速度计算预设时间步长内的火头、火翼和火尾对应的火场蔓延距离,包括:
基于预设时间步长构建惠更斯椭圆,以待预测火点作为椭圆焦点,火头和火尾的蔓延方向与长轴保持一致,火翼蔓延方向与短轴保持一致,在预设时间步长内,将火头蔓延距离确定为蔓延最长距离,将火尾蔓延距离确定为蔓延最短距离,将火翼蔓延距离确定为蔓延最宽长度。
可选的,方法还包括:
基于初始待预测火点进行初次构建椭圆边界后,以火头蔓延方向为起始的360°内,以预设间隔确定的预设个数的方向确定为蔓延方向,并将每个蔓延方向再次作为顶点构建更新火点所对应的边界;
完成全部蔓延方向的边界计算后,将全部蔓延方向所对应的所有边界对应的外包络线确定为下一时间段的蔓延预测边界,以此类推进行蔓延边界的模拟预测。
第二方面,本申请实施例提供了一种森林草原火场蔓延预测的装置,装置包括:
数据获取及预处理模块,用于获取待预测火点信息,并对待预测火点信息与地形数据、植被数据、气象预报数据进行分辨率统一处理;
蔓延方向确定模块,用于根据待预测火点信息和时间信息,通过气象预报数据和地形数据获取对应时刻、对应位置的风向和坡向,并基于风向和坡向合成矢量确定火场蔓延方向;
蔓延速度计算模块,用于通过改进的王正非模型计算火场蔓延速度;
蔓延边界计算模块,用于基于火场蔓延速度计算预设时间步长内的火头、火翼和火尾对应的火场蔓延距离,并基于火场蔓延距离计算蔓延边界。
第三方面,本申请实施例提供了一种电子设备,包括处理器和存储器,存储器存储有能够被处理器执行的计算机可执行指令,处理器执行计算机可执行指令以实现前述是实施方式中任一项的森林草原火场蔓延预测的方法。
第四方面,本申请实施例提供了一种计算机可读存储介质,计算机可读存储介质存储有计算机可执行指令,计算机可执行指令在被处理器调用和执行时,计算机可执行指令促使处理器实现前述实施方式中任一项的森林草原火场蔓延预测的方法。
本申请提供的森林草原火场蔓延预测的方法、装置、设备和介质,通过基于改进的王正非模型计算蔓延速度,结合基于惠更斯模型校正蔓延边界的综合模型,实现实时模拟预测火场蔓延信息,统计-物理模型的综合使用提高了模型的准确性,进而提升了火场蔓延预测的准确性。
附图说明
为了更清楚地说明本申请具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本申请的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本申请实施例提供的一种森林草原火场蔓延预测的方法的流程图;
图2为本申请实施例提供的一种蔓延方向的示意图;
图3为本申请实施例提供的一种具体的森林草原火场蔓延预测方法的流程图;
图4为本申请实施例提供的一种蔓延速度计算的流程图;
图5为本申请实施例提供的一种蔓延距离的示意图;
图6为本申请实施例提供的一种惠更斯原理应用于火场蔓延的示意图;
图7为本申请实施例提供的一种森林草原火场蔓延预测的装置的结构图;
图8为本申请实施例提供的一种电子设备的结构图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本申请一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本申请实施例的组件可以以各种不同的配置来布置和设计。
因此,以下对在附图中提供的本申请的实施例的详细描述并非旨在限制要求保护的本申请的范围,而是仅仅表示本申请的选定实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。
森林草原火灾具有突发性、灾害发生的随机性、短时间内能造成巨大损失的特点。如何提升森林草原防火现代化水平,将森林草原火灾带来的损失减少到最小,是目前急需解决的问题。目前国内常用的比较经典的蔓延模型主要是Rothermel模型和王正非模型,Rothermel模型属于半物理模型,物理模型物理参数较多又难以确定,计算过程需要众多的火场信息支撑,较为复杂,不满足及时性的蔓延模拟及蔓延预测需求。王正非统计模型主要依据大量历史统计数据,资料充足,置信度较高,蔓延速度计算过程简化,为模型的实时使用提供了可能,但其模型中使用的气象要素较为单一。
使用单一的蔓延速度模型进行火灾蔓延预测,适用范围有限,无法适应蔓延的复杂应用场景,相关技术中,常用的蔓延边界模型是元胞自动机栅格模型,存在元胞状态表达不完全的问题,同时元胞的蔓延方向仅为相邻的八个方向,造成模拟的最终边界细节缺失,模型的准确性受到影响。
基于此,本申请实施例提供了一种森林草原火场蔓延预测的方法、装置、设备和介质,可以提升森林草原火灾蔓延预测的准确性。
本申请实施例提供了一种森林草原火场蔓延预测的方法,参见图1所示,该方法主要包括以下步骤:
步骤S110,获取待预测火点信息,并对待预测火点信息与地形数据、植被数据、气象预报数据进行分辨率统一处理。
上述待预测火点信息可以包括待预测火点的位置、名称、火点下垫面覆盖类型等属性信息。地形数据可包括坡度、坡向等数据,植被数据诸如可以包括植被指数数据、植被覆盖度数据等植被信息,气象预报数据诸如可以包括风速、风向、气温、湿度要素等信息。
在一种实施方式中,分辨率统一处理可以采用重采样的方式进行处理,通过将上述数据进行重采样,可以使得各种数据处于同一分辨率,便于后续进行处理,提升了处理结果的准确度。
步骤S120,根据待预测火点信息和时间信息,通过气象预报数据和地形数据获取对应时刻、对应位置的风向和坡向,并基于风向和坡向合成矢量确定火场蔓延方向。
在一种实施方式中,根据待预测火点地理位置和时间信息,通过气象预报数据和地形数据获取对应时刻、位置的风向及坡向,根据惠更斯原理火头蔓延的方向由风向和坡向矢量合成得到,参见图2所示。
在一种实施方式中,基于风向和坡向合成矢量确定火场蔓延方向,可以用如下公式1计算:
(1)
其中,spread head_dir为火头蔓延方向;wind dir为风向;Aspect为下垫面坡向。
步骤S130,通过改进的王正非模型计算火场蔓延速度。
上述改进的王正非模型为引入气温校正因子和相对湿度校正因子后的模型。该模型缓解了原有的王正非模型中使用的气象要素较为单一的问题,提升了模型可靠性,同时可以更好的模拟蔓延实际情况。
在一种实施方式中,通过改进的王正非模型计算火场蔓延速度,可以采用V real=V ×F burn ×F wind ×F slope ×F cover ×F tem ×F rh进行计算,其中,V real为火头蔓延速度(m/min);V为蔓延初速度(m/min);F burn为可燃物类型校正因子,F wind为风力校正因子,F slope为地形坡度校正因子,F cover为下垫面类型校正因子,F tem为气温校正因子,F rh为湿度校正因子。
上述改进的王正非模型包含的气温校正因子和相对湿度校正因子如下:
气温校正因子:
其中,为气温等级划分,/>为各气温等级对应的系数因子,气温等级越高数值越大;
相对湿度校正因子:
其中,为相对湿度等级划分,/>为等级对应的系数因子,相对湿度等级越高数值越小。
进一步,在计算得到上述V real后,可以基于如下公式分别确定火头蔓延速度、火翼蔓延速度和火尾蔓延速度:
其中,Vhead为火头蔓延速度(m/min),Vwing为火翼蔓延速度(m/min),Vrear为火尾蔓延速度(m/min),m、n为不同风力等级对应的系数因子,风力等级越高数值越小。
步骤S140,基于火场蔓延速度计算预设时间步长内的火头、火翼和火尾对应的火场蔓延距离,并基于火场蔓延距离计算蔓延边界。
在一种实施方式中,可以基于预设时间步长构建惠更斯椭圆,以待预测火点作为椭圆焦点,火头和火尾的蔓延方向与长轴保持一致,火翼蔓延方向与短轴保持一致,在预设时间步长内,将头蔓延距离确定为蔓延最长距离,将火尾蔓延距离确定为蔓延最短距离,将火翼蔓延距离确定为蔓延最宽长度。
进一步,考虑到火场(尤其是森林草原火)通常会不断蔓延,因此基于初始待预测火点进行初次构建椭圆边界后,以火头蔓延方向为起始的360°内,以预设间隔确定的预设个数的方向确定为蔓延方向,并将每个蔓延方向再次作为顶点构建更新火点所对应的边界;完成全部蔓延方向的边界计算后,将全部蔓延方向所对应的所有边界对应的外包络线确定为下一时间段的蔓延预测边界,以此类推进行蔓延边界的模拟预测。
为便于理解,本申请实施例还提供了一种具体的森林草原火场蔓延预测方法,参见图3所示,包括以下步骤1至步骤4:
步骤1: 输入实时火点信息数据;
步骤2: 确定蔓延方向;
步骤3: 根据改进后的王正非模型进行蔓延速度计算;
步骤4: 确定蔓延步长,根据惠更斯原理进行蔓延边界预测。
以下对各个步骤的具体实现进行详细说明:
步骤1,输入实时火点信息数据,包括以下步骤1A至步骤1F:
步骤1A 输入实时火点信息,解析;
步骤1B 输入火点下垫面覆盖类型数据,解析;
步骤1C 输入地形数据,解析、提取坡度、坡向;
步骤1D 输入植被指数数据(NDVI)、植被覆盖度数据(VFC),解析;
步骤1E 输入气象预报数据,解析、提取风速、风向、气温、湿度要素;
步骤1F 采用重采样方法对上述数据进行处理,保持所有输入数据处于同一分辨率。
步骤2,确定蔓延方向,包括以下步骤2A和步骤2B:
步骤2A 根据火点地理位置和时间信息,通过气象预报数据和地形数据获取对应时刻、位置的风向及坡向,依据惠更斯原理火头蔓延的方向由风向和坡向矢量合成得到;
计算公式如下:
其中,spread head_dir为火头蔓延方向;wind dir为风向,根据气象预报数据获取,取值范围为0~360°;Aspect为下垫面坡向,根据地形数据获取,取值为-1时代表平坡,取值为0~360°代表不同坡向。
步骤2B 确定火头蔓延方向后,根据以下公式获取蔓延计算方向,以火头蔓延方向为起始的360°内,以每0.1°为间隔共3600个方向作为蔓延计算方向,后续分别对每个方向进行蔓延速度和边界的计算:
步骤3,根据改进后的王正非模型进行蔓延速度计算,参见图4所示:
步骤3A 在王正非模型中蔓延速度会受到下垫面类型、地形、气象等各种条件的影响,蔓延实际速度(Vreal)由当前时刻的蔓延速度(V)进行各项影响因子的修正得到,但在蔓延所处实际环境中,气温和相对湿度也会对可燃物的燃烧状态造成影响,王正非模型中只考虑了风速、风向两个较为单一的气象因子,为增加模型的全面性与可靠性,同时更好的模拟蔓延实际情况,在王正非模型基础上增加气温因子和湿度因子,引入与燃烧过程相关的其他特征要素,其公式如下:
V real=V×F burn ×F wind ×F slope ×F cover ×F tem ×F rh;
其中,V real为火头蔓延速度(m/min);V为蔓延初速度(m/min);F burn为可燃物类型校正因子,根据土地覆盖类型数据和植被覆盖度数据获取,F wind为风力校正因子,F slope为地形坡度校正因子,F cover为下垫面类型校正因子,根据土地覆盖类型数据和植被指数数据获取,F tem为气温校正因子,F rh为湿度校正因子。
具体的,可燃物类型校正因子计算公式如下:
其中,为土地覆盖类型因子,根据土地覆盖类型数据获取,不同的地物类型对应不同的数值因子,越接近可燃物类型数值越高;/>为植被覆盖度,/>为植被覆盖度阈值。
风力校正因子计算公式如下:
其中,为风力等级,依据风速大小进行风力等级划分;/>为各风力等级对应的系数因子,风力等级越高数值越大。
地形坡度校正因子计算公式如下:
其中,为上坡,蔓延方向与坡向相反,/>为上坡时不同坡度对应的系数因子,坡度越大时对应的数值越大;/>为下坡,蔓延方向与坡向一致,/>为下坡时不同坡度对应的系数因子,坡度越大时对应的数值越小。
下垫面类型校正因子计算公式如下:
其中,为植被指数,当下垫面地物属性为不可燃烧时,火的蔓延不会发生或停止,校正因子为0;可燃时校正因子由/>数据决定。
气温校正因子计算公式如下:
其中,为气温等级划分,/>为各气温等级对应的系数因子,气温等级越高数值越大。
相对湿度校正因子计算公式如下:
其中,为相对湿度等级划分,/>为等级对应的系数因子,相对湿度等级越高数值越小。
步骤3B 火场是由火头、火翼和火尾组成的,根据火头、火翼和火尾蔓延速度之间的经验关系,可以计算出火场边界位置。关系如下:
其中,为火头蔓延速度(m/min);/>为火翼蔓延速度(m/min);/>为火尾蔓延速度(m/min);/>、/>为不同风力等级对应的系数因子,风力等级越高数值越小。
步骤4,确定蔓延步长,根据惠更斯原理进行蔓延边界预测:
步骤4A 常用的栅格像元的模拟技术将火蔓延看作像元间的相互作用,由“燃烧”像元蔓延至“未燃烧”像元,不考虑连续面的传播,受到网格像元几何形状的影响,造成蔓延边界失真从而影响蔓延预测结果。惠更斯原理则通过椭圆波动传播“燃烧”状态,通过每个边界点和时间步长构建椭圆连续边界,与基于栅格像元的模拟技术相比具有更高的精度。
确定一定的时间蔓延步长,构建惠更斯椭圆,以火点作为椭圆焦点,火头和火尾的蔓延方向与长轴保持一致,火翼蔓延方向与短轴保持一致,在一定时间步长内,火头蔓延距离即蔓延的最长距离,火尾蔓延距离即蔓延的最短距离,火翼蔓延距离即蔓延的最宽长度,如图5所示。
根据上一步骤改进后的王正非模型计算的蔓延速度计算一定的时间步长内的火头、火尾、火翼的蔓延距离,即可得到椭圆的长轴(2a)、短轴(2b)的长度,根据下列椭圆公式构建惠更斯椭圆边界。
步骤4B 火场燃烧时是以椭圆形式向四周蔓延,在基于初始火点进行初次构建椭圆边界后,根据步骤2B确定的蔓延方向可在该椭圆边界确定3600个顶点,在下一时间步长内以这些顶点作为新的火点构建新的椭圆边界,完成全部方向计算后,获取所有椭圆的外包络线即为下一时间段的蔓延预测边界,以此类推进行蔓延边界的模拟预测。
以此获取不同预测时间的蔓延结果包括:曲线闭合连续的蔓延矢量边界、蔓延速度。在后续蔓延面积估算时,使用栅格像元形状进行蔓延面积估算时会造成较大的误差累积,而使用矢量边界造成的影响较小;同时矢量边界可在web或各软件进行展示,可视化方便快捷,相较于栅格更适合三维展示方式。
图6示出了一种惠更斯原理应用于火蔓延的示意图。
上述方式通过建立基于改进后的王正非模型计算蔓延速度和基于惠更斯模型校正蔓延边界的综合模型,实现实时模拟预测火场蔓延信息,统计-物理模型的综合使用满足复杂场景需求并提高了准确性。在王正非模型基础上增加气温因子和湿度因子,更全面的考虑气象要素的影响,使模型更贴合蔓延实际情况,提高模型精度。综合模型融入惠更斯矢量模型,保留物理波动模型特性,相较于普遍使用的栅格模型蔓延边界结果更精细、流畅,获得更逼近现实的蔓延边界。
本申请实施例为增加模型可靠性,同时更好的模拟蔓延实际情况,在王正非模型基础上增加气温因子和湿度因子,采用修正后的王正非模型来计算蔓延速度。采用惠更斯波动理论来模拟蔓延趋势,惠更斯是基于矢量的模型,火线被认为是连接顶点的闭合曲线,其中每个顶点都被视为独立椭圆扩展的潜在源,根据不同的精度需求可以调整蔓延方向的数量,能更好的模拟和预测实际的蔓延边界。最终建立基于改进后的王正非模型计算蔓延速度和基于惠更斯模型校正蔓延边界的综合模型,实现实时模拟预测火场蔓延信息,统计-物理模型的综合使用提高了模型的准确性。
基于上述方法实施例,本申请实施例还提供一种森林草原火场蔓延预测的装置,参见图7所示,该装置主要包括以下部分:
数据获取及预处理模块710,用于获取待预测火点信息,并对待预测火点信息与地形数据、植被数据、气象预报数据进行分辨率统一处理;
蔓延方向确定模块720,用于根据待预测火点信息和时间信息,通过气象预报数据和地形数据获取对应时刻、对应位置的风向和坡向,并基于风向和坡向合成矢量确定火场蔓延方向;
蔓延速度计算模块730,用于通过改进的王正非模型计算火灾蔓延速度;
蔓延边界计算模块740,用于基于火灾蔓延速度计算预设时间步长内的火头、火翼和火尾对应的火灾蔓延距离,并基于火灾蔓延距离计算蔓延边界。
可选的,上述蔓延方向确定模块720,还用于:
其中,spread head_dir为火头蔓延方向;wind dir为风向;Aspect为下垫面坡向。
可选的,上述蔓延速度计算模块730,还用于:
V real=V×F burn ×F wind ×F slope ×F cover ×F tem ×F rh;
其中,V real为火头蔓延速度(m/min);V为蔓延初速度(m/min);F burn为可燃物类型校正因子,F wind为风力校正因子,F slope为地形坡度校正因子,F cover为下垫面类型校正因子,F tem为气温校正因子,F rh为湿度校正因子。
可选的,改进的王正非模型包含的气温校正因子和相对湿度校正因子如下:
气温校正因子:
其中,为气温等级划分,/>为各气温等级对应的系数因子,气温等级越高数值越大;
相对湿度校正因子:
其中,为相对湿度等级划分,/>为等级对应的系数因子,相对湿度等级越高数值越小。
可选的,上述蔓延速度计算模块730,还用于:基于如下公式分别确定火头蔓延速度、火翼蔓延速度和火尾蔓延速度:
其中,Vhead为火头蔓延速度(m/min),Vwing为火翼蔓延速度(m/min),Vrear为火尾蔓延速度(m/min),m、n为不同风力等级对应的系数因子,风力等级越高数值越小。
可选的,上述蔓延边界计算模块740,还用于:
基于预设时间步长构建惠更斯椭圆,以待预测火点作为椭圆焦点,火头和火尾的蔓延方向与长轴保持一致,火翼蔓延方向与短轴保持一致,在预设时间步长内,将头蔓延距离确定为蔓延最长距离,将火尾蔓延距离确定为蔓延最短距离,将火翼蔓延距离确定为蔓延最宽长度。
可选的,上述装置还包括,边界迭代模块,用于:
基于初始待预测火点进行初次构建椭圆边界后,以火头蔓延方向为起始的360°内,以预设间隔确定的预设个数的方向确定为蔓延方向,并将每个蔓延方向再次作为顶点构建更新火点所对应的边界;
完成全部蔓延方向的边界计算后,将全部蔓延方向所对应的所有边界对应的外包络线确定为下一时间段的蔓延预测边界,以此类推进行蔓延边界的模拟预测。
本申请实施例提供的森林草原火场蔓延预测的装置,其实现原理及产生的技术效果和前述方法实施例相同,为简要描述,森林草原火场蔓延预测的装置的实施例部分未提及之处,可参考前述森林草原火场蔓延预测的方法实施例中相应内容。
本申请实施例还提供了一种电子设备,如图8所示,为该电子设备的结构示意图,其中,该电子设备100包括处理器81和存储器80,该存储器80存储有能够被该处理器81执行的计算机可执行指令,该处理器81执行该计算机可执行指令以实现上述任一项森林草原火场蔓延预测的方法。
在图8示出的实施方式中,该电子设备还包括总线82和通信接口83,其中,处理器81、通信接口83和存储器80通过总线82连接。
其中,存储器80可能包含高速随机存取存储器(RAM,Random Access Memory),也可能还包括非不稳定的存储器(non-volatile memory),例如至少一个磁盘存储器。通过至少一个通信接口83(可以是有线或者无线)实现该系统网元与至少一个其他网元之间的通信连接,可以使用互联网,广域网,本地网,城域网等。总线82可以是ISA(IndustryStandard Architecture,工业标准体系结构)总线、PCI(Peripheral ComponentInterconnect,外设部件互连标准)总线或EISA(Extended Industry StandardArchitecture,扩展工业标准结构)总线等。所述总线82可以分为地址总线、数据总线、控制总线等。为便于表示,图8中仅用一个双向箭头表示,但并不表示仅有一根总线或一种类型的总线。
处理器81可能是一种集成电路芯片,具有信号的处理能力。在实现过程中,上述方法的各步骤可以通过处理器81中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器81可以是通用处理器,包括中央处理器(Central Processing Unit,简称CPU)、网络处理器(Network Processor,简称NP)等;还可以是数字信号处理器(Digital SignalProcessor,简称DSP)、专用集成电路(Application Specific Integrated Circuit,简称ASIC)、现场可编程门阵列(Field-Programmable Gate Array,简称FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合本申请实施例所公开的方法的步骤可以直接体现为硬件译码处理器执行完成,或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器,闪存、只读存储器,可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器,处理器81读取存储器中的信息,结合其硬件完成前述实施例的火灾蔓延预测的方法的步骤。
本申请实施例还提供了一种计算机可读存储介质,该计算机可读存储介质存储有计算机可执行指令,该计算机可执行指令在被处理器调用和执行时,该计算机可执行指令促使处理器实现上述森林草原火场蔓延预测的方法,具体实现可参见前述方法实施例,在此不再赘述。
本申请实施例所提供的森林草原火场蔓延预测的方法、装置、设备和介质的计算机程序产品,包括存储了程序代码的计算机可读存储介质,所述程序代码包括的指令可用于执行前面方法实施例中所述的方法,具体实现可参见方法实施例,在此不再赘述。
除非另外具体说明,否则在这些实施例中阐述的部件和步骤的相对步骤、数字表达式和数值并不限制本申请的范围。
所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个处理器可执行的非易失的计算机可读取存储介质中。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
在本申请的描述中,需要说明的是,术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,或者是该发明产品使用时惯常摆放的方位或位置关系,仅是为了便于描述本申请和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本申请的限制。
最后应说明的是:以上各实施例仅用以说明本申请的技术方案,而非对其限制;尽管参照前述各实施例对本申请进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本申请各实施例技术方案的范围。
Claims (8)
1.一种森林草原火场蔓延预测的方法,其特征在于,所述方法包括:
获取待预测火点信息,并对所述待预测火点信息与地形数据、植被数据、气象预报数据进行分辨率统一处理;
根据所述待预测火点信息和时间信息,通过所述气象预报数据和所述地形数据获取对应时刻、对应位置的风向和坡向,并基于所述风向和所述坡向合成矢量确定火场蔓延方向;
通过改进的王正非模型计算火场蔓延速度,包括:
V real=V×F burn ×F wind ×F slope ×F cover ×F tem ×F rh;
其中,V real为火头蔓延速度(m/min);V为蔓延初速度(m/min);F burn为可燃物类型校正因子,F wind为风力校正因子,F slope为地形坡度校正因子,F cover为下垫面类型校正因子,F tem为气温校正因子,F rh为湿度校正因子;
所述改进的王正非模型包含的气温校正因子和相对湿度校正因子如下:
气温校正因子:
其中,为气温等级划分,/>为各气温等级对应的系数因子,气温等级越高数值越大;
相对湿度校正因子:
其中,为相对湿度等级划分,/>为等级对应的系数因子,相对湿度等级越高数值越小;
基于所述火场蔓延速度计算预设时间步长内的火头、火翼和火尾对应的火场蔓延距离,并基于所述火场蔓延距离计算蔓延边界。
2.根据权利要求1所述的森林草原火场蔓延预测的方法,其特征在于,基于所述风向和所述坡向合成矢量确定火场蔓延方向,包括:
其中,spread head_dir为火头蔓延方向;wind dir为风向;Aspect为下垫面坡向。
3.根据权利要求1所述的森林草原火场蔓延预测的方法,其特征在于,所述方法还包括,基于如下公式分别确定火头蔓延速度、火翼蔓延速度和火尾蔓延速度:
其中,Vhead为火头蔓延速度(m/min),Vwing为火翼蔓延速度(m/min),Vrear为火尾蔓延速度(m/min),m、n为不同风力等级对应的系数因子,风力等级越高数值越小。
4.根据权利要求1所述的森林草原火场蔓延预测的方法,其特征在于,基于所述火场蔓延速度计算预设时间步长内的火头、火翼和火尾对应的火场蔓延距离,包括:
基于预设时间步长构建惠更斯椭圆,以待预测火点作为椭圆焦点,火头和火尾的蔓延方向与长轴保持一致,火翼蔓延方向与短轴保持一致,在预设时间步长内,将火头蔓延距离确定为蔓延最长距离,将火尾蔓延距离确定为蔓延最短距离,将火翼蔓延距离确定为蔓延最宽长度。
5.根据权利要求4所述的森林草原火场蔓延预测的方法,其特征在于,所述方法还包括:
基于初始待预测火点进行初次构建椭圆边界后,以火头蔓延方向为起始的360°内,以预设间隔确定的预设个数的方向确定为蔓延方向,并将每个蔓延方向再次作为顶点构建更新火点所对应的边界;
完成全部蔓延方向的边界计算后,将全部蔓延方向所对应的所有边界对应的外包络线确定为下一时间段的蔓延预测边界,以此类推进行蔓延边界的模拟预测。
6.一种森林草原火场蔓延预测的装置,其特征在于,所述装置包括:
数据获取及预处理模块,用于获取待预测火点信息,并对所述待预测火点信息与地形数据、植被数据、气象预报数据进行分辨率统一处理;
蔓延方向确定模块,用于根据所述待预测火点信息和时间信息,通过所述气象预报数据和所述地形数据获取对应时刻、对应位置的风向和坡向,并基于所述风向和所述坡向合成矢量确定火灾蔓延方向;
蔓延速度计算模块,用于通过改进的王正非模型计算火场蔓延速度,包括:
V real=V×F burn ×F wind ×F slope ×F cover ×F tem ×F rh;
其中,V real为火头蔓延速度(m/min);V为蔓延初速度(m/min);F burn为可燃物类型校正因子,F wind为风力校正因子,F slope为地形坡度校正因子,F cover为下垫面类型校正因子,F tem为气温校正因子,F rh为湿度校正因子;
所述改进的王正非模型包含的气温校正因子和相对湿度校正因子如下:
气温校正因子:
其中,为气温等级划分,/>为各气温等级对应的系数因子,气温等级越高数值越大;
相对湿度校正因子:
其中,为相对湿度等级划分,/>为等级对应的系数因子,相对湿度等级越高数值越小;
蔓延边界计算模块,用于基于所述火场蔓延速度计算预设时间步长内的火头、火翼和火尾对应的火场蔓延距离,并基于所述火场蔓延距离计算蔓延边界。
7.一种电子设备,其特征在于,包括处理器和存储器,所述存储器存储有能够被所述处理器执行的计算机可执行指令,所述处理器执行所述计算机可执行指令以实现权利要求1至5任一项所述的森林草原火场蔓延预测的方法。
8.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有计算机可执行指令,所述计算机可执行指令在被处理器调用和执行时,计算机可执行指令促使处理器实现权利要求1至5任一项所述的森林草原火场蔓延预测的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310994183.2A CN116720388B (zh) | 2023-08-09 | 2023-08-09 | 森林草原火场蔓延预测的方法、装置、设备和介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310994183.2A CN116720388B (zh) | 2023-08-09 | 2023-08-09 | 森林草原火场蔓延预测的方法、装置、设备和介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116720388A CN116720388A (zh) | 2023-09-08 |
CN116720388B true CN116720388B (zh) | 2023-11-03 |
Family
ID=87870025
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310994183.2A Active CN116720388B (zh) | 2023-08-09 | 2023-08-09 | 森林草原火场蔓延预测的方法、装置、设备和介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116720388B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20100087567A (ko) * | 2009-01-28 | 2010-08-05 | 서울대학교산학협력단 | 산불확산 예측시스템 및 방법 |
CN106780732A (zh) * | 2016-12-22 | 2017-05-31 | 山东神戎电子股份有限公司 | 一种基于三维地形数据的小班林火蔓延动态模拟方法 |
CN111310344A (zh) * | 2020-02-22 | 2020-06-19 | 福州大学 | 一种林火蔓延模拟中考虑风场与火场耦合作用的方法 |
CN113902981A (zh) * | 2021-12-10 | 2022-01-07 | 航天宏图信息技术股份有限公司 | 一种林火蔓延预测方法、装置、电子设备及存储介质 |
CN115577830A (zh) * | 2022-09-16 | 2023-01-06 | 浙江大华技术股份有限公司 | 火势蔓延范围预测方法、装置、计算机设备及存储介质 |
CN116305832A (zh) * | 2023-02-16 | 2023-06-23 | 浙江大学 | 一种林火蔓延实时模拟方法、灭火决策方法与预警系统 |
CN116363300A (zh) * | 2023-02-15 | 2023-06-30 | 中国建筑科学研究院有限公司 | 一种基于遥感影像下的森林火场模拟方法 |
-
2023
- 2023-08-09 CN CN202310994183.2A patent/CN116720388B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20100087567A (ko) * | 2009-01-28 | 2010-08-05 | 서울대학교산학협력단 | 산불확산 예측시스템 및 방법 |
CN106780732A (zh) * | 2016-12-22 | 2017-05-31 | 山东神戎电子股份有限公司 | 一种基于三维地形数据的小班林火蔓延动态模拟方法 |
CN111310344A (zh) * | 2020-02-22 | 2020-06-19 | 福州大学 | 一种林火蔓延模拟中考虑风场与火场耦合作用的方法 |
CN113902981A (zh) * | 2021-12-10 | 2022-01-07 | 航天宏图信息技术股份有限公司 | 一种林火蔓延预测方法、装置、电子设备及存储介质 |
CN115577830A (zh) * | 2022-09-16 | 2023-01-06 | 浙江大华技术股份有限公司 | 火势蔓延范围预测方法、装置、计算机设备及存储介质 |
CN116363300A (zh) * | 2023-02-15 | 2023-06-30 | 中国建筑科学研究院有限公司 | 一种基于遥感影像下的森林火场模拟方法 |
CN116305832A (zh) * | 2023-02-16 | 2023-06-23 | 浙江大学 | 一种林火蔓延实时模拟方法、灭火决策方法与预警系统 |
Also Published As
Publication number | Publication date |
---|---|
CN116720388A (zh) | 2023-09-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107003834B (zh) | 行人检测设备和方法 | |
CN105469446A (zh) | 点云网格简化系统及方法 | |
CN113792427A (zh) | 一种已有边界条件下森林火灾蔓延预测模拟的方法 | |
CN115329691B (zh) | 一种基于cfd与gis的超大城市风环境模拟方法 | |
CN114936502B (zh) | 一种森林火灾蔓延态势边界分析方法、系统、终端及介质 | |
CN113065974A (zh) | 一种基于动态网络表示学习的链路预测方法 | |
CN106874602A (zh) | 气象数据处理方法和装置 | |
CN115082699A (zh) | 一种轮廓形状提取方法、装置、电子设备及存储介质 | |
CN113268800B (zh) | 优化区域规划方法、装置、电子设备和可读存储介质 | |
CN116720388B (zh) | 森林草原火场蔓延预测的方法、装置、设备和介质 | |
CN116663330B (zh) | 时空同步多因素耦合的火蔓延模拟方法、装置及电子设备 | |
CN116051699B (zh) | 动捕数据的处理方法、装置、设备及存储介质 | |
CN117173593A (zh) | 基于遥感地物分布感知的火势蔓延推演方法、介质及设备 | |
KR101800789B1 (ko) | 풍속 내삽에 적용하기 위해 개선된 mk-prism 방법 | |
CN115204507B (zh) | 大气能见度的预测方法、装置、设备及存储介质 | |
CN114637305B (zh) | 一种无人机最短路径规划方法及装置 | |
CN116432093A (zh) | 风切变指数的确定方法、装置和电子设备 | |
CN112711783B (zh) | 一种基于三维软件及虚幻引擎的珠宝设计方法及系统 | |
CN111581808B (zh) | 一种污染物信息处理方法、装置、存储介质及终端 | |
Alemayehu et al. | Autonomous urban topology generation for urban flow modelling | |
CN113775482B (zh) | 风电机组的最大风速确定方法、装置及电子设备 | |
CN109102558A (zh) | 一种基于非结构三角网的水流流线快速生成方法 | |
CN113887843B (zh) | 一种风电场输出功率的群体预测方法和系统 | |
Cao et al. | Evaluation and optimization of outdoor wind environment in block based on space syntax and CFD simulation | |
CN118245770A (zh) | 多气象要素分布特征的生成方法、装置、设备及介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |