CN115099073B - Real-time forest fire spreading simulation method and system - Google Patents

Real-time forest fire spreading simulation method and system Download PDF

Info

Publication number
CN115099073B
CN115099073B CN202211022678.0A CN202211022678A CN115099073B CN 115099073 B CN115099073 B CN 115099073B CN 202211022678 A CN202211022678 A CN 202211022678A CN 115099073 B CN115099073 B CN 115099073B
Authority
CN
China
Prior art keywords
fire
information
fuel
change
forest
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
Application number
CN202211022678.0A
Other languages
Chinese (zh)
Other versions
CN115099073A (en
Inventor
杨邦会
胡乔利
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhongke Haihui Beijing Technology Co ltd
Original Assignee
Zhongke Haihui Beijing Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhongke Haihui Beijing Technology Co ltd filed Critical Zhongke Haihui Beijing Technology Co ltd
Priority to CN202211022678.0A priority Critical patent/CN115099073B/en
Publication of CN115099073A publication Critical patent/CN115099073A/en
Application granted granted Critical
Publication of CN115099073B publication Critical patent/CN115099073B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/10Services
    • G06Q50/26Government or public services
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/05Geographic models
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information 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)
  • Strategic Management (AREA)
  • Human Resources & Organizations (AREA)
  • Geometry (AREA)
  • Economics (AREA)
  • Tourism & Hospitality (AREA)
  • Software Systems (AREA)
  • Marketing (AREA)
  • Development Economics (AREA)
  • General Business, Economics & Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Quality & Reliability (AREA)
  • Educational Administration (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Operations Research (AREA)
  • General Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Game Theory and Decision Science (AREA)
  • Computer Graphics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention relates to the field of virtual reality technology and fire prevention and control, in particular to a real-time forest fire spreading simulation method and system. The method comprises the following steps: determining a target forest; acquiring relevant information of the target forest through a remote sensing satellite and ground communication equipment; obtaining overall fire spreading information by using the related information; setting a prediction time length; predicting a change in fire climate information and a change in fire fuel information within the predicted time period; and obtaining the whole fire spreading situation according to the change, and simultaneously carrying out real-time simulation on the fire spreading situation. The method combines and maps meteorological factors, topographic factors and fuel factors closely related to the forest fire and the fire condition, thereby accurately and quickly simulating the spreading condition of the fire, not only providing decision support for fire departments to put out the forest fire, but also having great significance and effect on forest disaster prevention and reduction.

Description

Real-time forest fire spreading simulation method and system
Technical Field
The invention relates to the field of virtual reality technology and fire prevention and control, in particular to a real-time forest fire spreading simulation method and system.
Background
The forest is the biggest ecosystem main body on land, is the support of the life of the earth, is an important component of national resources and national wealth, is commonly and vividly compared with the lung of the earth, and plays a decisive role in maintaining and improving the global ecological environment. Due to global warming, the combustible storage in the forest is continuously increased, the risk of forest fire is aggravated, and the forest fire not only destroys the forest ecosystem and causes forest resource loss, but also threatens the life and property safety of people. The research on forest fire spreading is of great significance for extinguishing forest fires, but in the prior art, aiming at forest fires which are not generated at present during forest fire spreading, a technology for performing real-time simulation by combining meteorological factors, topographic factors and fuel factors is needed urgently.
Disclosure of Invention
Aiming at the defects in the prior art, the invention provides a real-time forest fire spreading simulation method in a first aspect, which comprises the following steps: determining a target forest; acquiring real-time fire information, fire geographical information, fire fuel information and fire climate information of the target forest through a remote sensing satellite and ground communication equipment; acquiring overall fire spreading information by using the real-time fire information, the fire geographical information, the fire fuel information and the fire climate information; setting a prediction time length; calculating the change of the fire climate information and the change of the fire fuel information within the prediction duration according to the prediction duration by combining the overall fire spread information, the fire fuel information and the fire climate information; predicting the change of the overall fire spread information after the prediction duration according to the overall fire spread information, the change of the fire climate information and the change of the fire fuel information; and simulating the fire spreading condition of the target forest according to the whole fire spreading information and the change of the whole fire spreading information. The method of the invention predicts the forest fire spreading by combining the real-time fire information, the fire geographical information, the fire fuel information and the fire climate information, namely combining and mapping the meteorological factors, the terrain factors and the fuel factors which are closely related to the forest fire and the fire condition, thereby accurately and quickly simulating the fire spreading condition, improving the real-time and efficient requirements of forest fire spreading simulation, providing decision support for fire control departments to suppress the forest fire, and simultaneously having great significance and effect in forest disaster prevention and reduction.
Optionally, the obtaining of the real-time fire information, the fire geographic information, the fire fuel information, and the fire climate information of the target forest by using a remote sensing satellite and a ground communication device includes the following steps: acquiring fire point information and a fire boundary of the target forest through the remote sensing satellite; the fire point information and the fire boundary are arranged, and the real-time fire information is obtained; acquiring the gradient, the slope direction and DEM data of the target forest through the ground communication equipment; combining the slope, the slope direction and the DEM data to obtain the fire disaster geographic information; acquiring the loading capacity, the density of dried particles, the surface area-to-volume ratio, the depth of a combustible bed, the low heat content of the combustible, the total mineral content, the effective mineral content, the water content of the combustible and the extinguishing water content of the combustible in the target forest by using the ground communication equipment; converging the dried combustible load, the dried particle density, the surface area to volume ratio, the combustible bed depth, the combustible low heat content, the total mineral content, the effective mineral content, the combustible water content, and the combustible extinction water content to obtain the fire fuel information; acquiring the wind speed, the wind direction, the dryness and the precipitation of the target forest by using the ground communication equipment; and acquiring the fire climate information through the wind speed, the wind direction, the dryness and the precipitation.
Optionally, the obtaining of overall fire spreading information by using the real-time fire information, the fire geographic information, the fire fuel information and the fire climate information includes: constructing a two-dimensional coordinate interface of the target forest by using the fire geographic information; setting fire point coordinates and fire boundary coordinates in the two-dimensional coordinate interface according to the real-time fire information; assigning corresponding combustion parameters to corresponding fire point coordinates and fire boundary coordinates according to the fire fuel information and the fire climate information; and acquiring overall fire spreading information by setting the fire point coordinates and the fire boundary coordinates and assigning a two-dimensional coordinate interface.
Optionally, the setting the predicted time period includes the following steps: and setting the predicted time length according to the time resolution of the remote sensing satellite, wherein the predicted time length is greater than the time resolution.
Optionally, the calculating, according to the predicted duration, the change of the fire climate information and the change of the fire fuel information within the predicted duration by combining the overall fire spread information, the fire fuel information, and the fire climate information includes: predicting the change of the fire climate information within the prediction duration according to the fire climate information; constructing a fuel loss function by using the overall fire spread information and the fire fuel information; and obtaining the change of the fire disaster fuel information in the prediction time length by combining the change of the fire disaster climate information through the fuel loss function.
Further optionally, the fuel loss function satisfies the following equation:
Figure 503712DEST_PATH_IMAGE001
wherein, the first and the second end of the pipe are connected with each other,
Figure 581389DEST_PATH_IMAGE002
the length of time is indicated by a time period,
Figure 738570DEST_PATH_IMAGE003
representing two-dimensional coordinates of
Figure 218093DEST_PATH_IMAGE004
The fuel loss function of the fire at the location,
Figure 920470DEST_PATH_IMAGE005
representing two-dimensional coordinates of
Figure 485443DEST_PATH_IMAGE004
The initial total reserve of combustible material at the location,
Figure 665889DEST_PATH_IMAGE006
the coefficient of fuel loss is expressed as,
Figure 514765DEST_PATH_IMAGE007
representing two-dimensional coordinates of
Figure 388043DEST_PATH_IMAGE004
The direction of the fire at the location,
Figure 705892DEST_PATH_IMAGE008
representing two-dimensional coordinates of
Figure 424449DEST_PATH_IMAGE004
The fuel distribution vector of the fire at the location,
Figure 144144DEST_PATH_IMAGE009
representing two-dimensional coordinates of
Figure 922744DEST_PATH_IMAGE004
Fuel storage function at the location.
Optionally, the predicting a change of the overall fire spread information after the predicted time period according to the overall fire spread information, the change of the fire climate information, and the change of the fire fuel information includes: constructing a fire point migration speed function by using the overall fire spread information; determining a head fire point coordinate, a tail fire point coordinate and a side fire point coordinate by combining the corresponding combustion parameters and according to the fire point coordinate and the fire boundary coordinate; selecting a plurality of coordinates on the fire boundary as auxiliary side fire point coordinates, and assigning corresponding combustion parameters to the auxiliary side fire point coordinates; taking the head fire point coordinate, the tail fire point coordinate, the side fire point coordinate and the auxiliary side fire point as a plurality of starting points, and combining the fire point migration speed function and the predicted time length to obtain the fire point coordinate after the predicted time length; according to the fire point coordinates after the predicted duration, combining the change of the fire climate information, and re-demarcating a fire boundary; and re-defining a fire boundary according to the fire point coordinates after the predicted time length, and obtaining the change of the overall fire spreading information.
Further optionally, the fire migration velocity function satisfies the following formula:
Figure 977156DEST_PATH_IMAGE010
wherein, the first and the second end of the pipe are connected with each other,
Figure 764984DEST_PATH_IMAGE011
representing two-dimensional coordinates of
Figure 73605DEST_PATH_IMAGE004
The speed of migration of the fire at the location,
Figure 288686DEST_PATH_IMAGE012
Figure 315548DEST_PATH_IMAGE013
indicating a fire response corrected intensity factor at a location in two-dimensional coordinates (x, y), wherein,
Figure 641487DEST_PATH_IMAGE014
the spread rate of the forest fire is shown,
Figure 53883DEST_PATH_IMAGE015
the wind speed correction coefficient is represented by a coefficient,
Figure 439865DEST_PATH_IMAGE016
the slope correction coefficient is represented by a slope correction coefficient,
Figure 219602DEST_PATH_IMAGE017
representing two-dimensional coordinates of
Figure 83653DEST_PATH_IMAGE004
The density of the fuel at the location(s),
Figure 366867DEST_PATH_IMAGE018
representing the amount of heat required to ignite a unit mass of fuel,
Figure 173017DEST_PATH_IMAGE019
representing the effective thermal coefficient of heat.
Optionally, the simulating the fire spreading condition of the target forest according to the overall fire spreading information and the change of the overall fire spreading information includes: constructing an initial fire condition of a target forest by using the overall fire spreading information; and simulating the fire spreading condition of the target forest by combining the initial fire condition with the change of the overall fire spreading information.
In a second aspect, the present invention further provides a real-time forest fire spreading simulation system, which is suitable for the real-time forest fire spreading simulation method in the first aspect, and includes: the information input module is used for determining a target forest; setting a prediction time length; the data acquisition module is used for acquiring real-time fire information, fire geographical information, fire fuel information and fire climate information of the target forest through a remote sensing satellite and ground communication equipment; the initial fire behavior module is used for acquiring overall fire spread information by utilizing the real-time fire information, the fire geographic information, the fire fuel information and the fire climate information; the first fire behavior simulation module is used for calculating the change of the fire climate information and the change of the fire fuel information within the prediction duration according to the prediction duration by combining the overall fire spread information, the fire fuel information and the fire climate information; the second fire behavior simulation module is used for predicting the change of the overall fire spread information after the prediction duration according to the overall fire spread information, the change of the fire climate information and the change of the fire fuel information; and the fire simulation module simulates the fire spreading condition of the target forest according to the whole fire spreading information and the change of the whole fire spreading information. The real-time forest fire spreading simulation system provided by the invention has the advantages of compact structure, accurate operation result and high efficiency, so that the real-time forest fire spreading simulation system has more practical application value and commercial value.
Drawings
FIG. 1 is a flow chart of a real-time simulation method for forest fire spreading according to the present invention;
FIG. 2 is a schematic diagram of an initial fire condition simulation of the present invention;
FIG. 3 is a schematic diagram of the interpolation of the coordinates of a lateral fire point after the spread of a fire according to the present invention;
FIG. 4 is a schematic diagram illustrating a comparison between an initial fire condition simulation and a fire spread condition according to the present invention;
FIG. 5 is a schematic diagram of a real-time forest fire spreading simulation system according to the present invention.
Detailed Description
Specific embodiments of the present invention will be described in detail below, and it should be noted that the embodiments described herein are only for illustration and are not intended to limit the present invention. In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present invention. However, it will be apparent to one of ordinary skill in the art that: it is not necessary to employ these specific details to practice the present invention. In other instances, well-known circuits, software, or methods have not been described in detail in order to avoid obscuring the present invention.
Throughout the specification, reference to "one embodiment," "an embodiment," "one example" or "an example" means: the particular features, structures, or characteristics described in connection with the embodiment or example are included in at least one embodiment of the invention. Thus, the appearances of the phrases "in one embodiment," "in an embodiment," "one example" or "an example" in various places throughout this specification are not necessarily all referring to the same embodiment or example. Furthermore, the particular features, structures, or characteristics may be combined in any suitable combination and/or sub-combination in one or more embodiments or examples. Further, those of ordinary skill in the art will appreciate that the illustrations provided herein are for illustrative purposes and are not necessarily drawn to scale.
In an alternative embodiment, referring to fig. 1, the present invention provides a real-time forest fire spreading simulation method, including the following steps:
s1, determining a target forest.
In an alternative embodiment, identifying the target forest includes determining a forest in which a fire is occurring and for which fire suppression is desired and setting the forest as the target forest, and obtaining geographic information about the forest, such as area size, latitude and longitude.
And S2, acquiring real-time fire information, fire geographical information, fire fuel information and fire climate information of the target forest through a remote sensing satellite and ground communication equipment.
Specifically, in this embodiment, the obtaining of the real-time fire information, the geographic information of the fire, the fuel information of the fire, and the climate information of the fire of the target forest by using a remote sensing satellite and a ground communication device includes the following steps: acquiring fire point information and a fire boundary of the target forest through the remote sensing satellite; the fire point information and the fire boundary are arranged, and the real-time fire information is obtained; acquiring gradient, slope direction and DEM data (Digital Elevation Model ) of the target forest by the ground communication equipment, wherein the Digital Elevation Model is used for realizing Digital simulation of ground topography, namely Digital expression of topography surface morphology through limited topography Elevation data; combining the slope, the slope direction and the DEM data to obtain the fire disaster geographic information; acquiring the loading capacity, the density of dried particles, the surface area-to-volume ratio, the depth of a combustible bed, the low heat content of the combustible, the total mineral content, the effective mineral content, the water content of the combustible and the extinguishing water content of the combustible in the target forest by using the ground communication equipment; converging the dried combustible load, the dried particle density, the surface area to volume ratio, the combustible bed depth, the combustible low heat content, the total mineral content, the effective mineral content, the combustible water content, and the combustible extinction water content to obtain the fire fuel information; acquiring the wind speed, the wind direction, the dryness and the precipitation of the target forest by using the ground communication equipment; and acquiring the fire climate information through the wind speed, the wind direction, the dryness and the precipitation.
In an optional embodiment, a high-resolution real-time fire image is obtained through known latitude and longitude information of a target forest by using remote sensing satellites, and the remote sensing satellites select polar orbit satellites which comprise Himapari-8, terra/Aqua, NPP, NOAA, landsat-8, sentinil 2 and other satellites. The polar orbit satellite can quickly read and accurately identify fire point information in forest fires. Specifically, data from a Himapari-8 satellite is selected as the source of telemetry satellite data, the Himapari-8 satellite having a temporal resolution of 10 minutes per time, while also having a spatial resolution of 500 meters. The method comprises the steps of preprocessing collected data by a Himapari-8 satellite, performing brightness temperature conversion to obtain surface brightness temperature data, generating a comprehensive abnormal heat source extraction result by combining reflectivity data of a visible light-infrared band and utilizing methods such as an absolute fire point identification algorithm, a context algorithm, a class variance, a fixed threshold value and a self-adaptive threshold value, and accordingly obtaining real-time fire information, namely real-time fire information, including parameters such as fire point information and fire boundary.
In yet another optional embodiment, the ground communication device may select the unmanned aerial vehicle to carry an infrared monitoring camera, a video camera, a real-time communication device, and the like, so as to complete the assembly, and detect the actual climate, the geographic environment, the vegetation status, and the combustion status of the target forest in real time and obtain the fire geographic information, the fire fuel information, and the fire climate information through the ground communication device according to the known longitude and latitude information of the target forest. The fire geographical information comprises the gradient, the slope direction and DEM data of the target forest. The fire fuel information comprises the capacity of the dry combustible materials in the target forest, the density of the dry particles, the surface area-volume ratio, the depth of the combustible material bed, the low heat content of the combustible materials, the total mineral content, the effective mineral content, the water content of the combustible materials and the water content for extinguishing the combustible materials; converging the oven dried combustible load, the oven dried particle density, the surface area to volume ratio, the combustible bed depth, the combustible low heat content, the total mineral content, the effective mineral content, the combustible water content, and the combustible quench water content; specifically, the dry combustible load represents the stock of dry combustible in the target forest, the unit is kilogram per square meter, and the value range is 0 to 10000; the density of the dried particles represents the storage of dry combustible particles in a target forest, the unit is kilogram per square meter, and the value range is 0 to 10000; the surface area-volume ratio represents the surface area-volume ratio of total combustible substances and combustible particles in the target forest, the unit is square meters per cubic meter, and the value range is 0 to 10000; the depth of the combustible bed indicates the depth of the combustible in the target forest below the ground, the unit is meter, and the value range is 0 to 5; the low heat content of the combustible substance represents the heat content of the combustible substance in the target forest, the unit is kilojoule per kilogram, and the value range is 0 to 10000; the total mineral content represents the ratio of the total mineral content in the target forest, and the value range is 0 to 100; the effective mineral content represents the ratio of combustible effective minerals in the total mineral content in the target forest, and the value range is 0 to 100; the combustible water content represents the water content ratio of the combustible in the target forest, and the value range is 0 to 100; the extinguishing water content of the combustible substance represents the water content of the tiny combustible substance in the target forest, the value range is 0 to 100, when the extinguishing water content of the combustible substance is more than 8 percent, the fire can not be caused, and when the extinguishing water content of the combustible substance is less than 4 percent, the fire is very easy to ignite. Meanwhile, the remote sensing satellite is used for specifically selecting the meteorological condition of the target forest obtained through the meteorological satellite, so that the fire climate information of the target forest is obtained and perfected. And the fire climate information comprises the wind speed, the wind direction, the dryness and the precipitation of the target forest.
And S3, obtaining overall fire spreading information by utilizing the real-time fire information, the fire geographic information, the fire fuel information and the fire climate information.
Specifically, in this embodiment, the obtaining of the overall fire spreading information by using the real-time fire information, the fire geographic information, the fire fuel information, and the fire climate information includes: constructing a two-dimensional coordinate interface of the target forest by using the fire geographic information; setting fire point coordinates and fire boundary coordinates in the two-dimensional coordinate interface according to the real-time fire information; assigning corresponding combustion parameters to corresponding fire point coordinates and fire boundary coordinates according to the fire fuel information and the fire climate information; and acquiring overall fire spreading information by setting the fire point coordinates and the fire boundary coordinates and assigning a two-dimensional coordinate interface.
In an optional embodiment, the fire geographical location information further includes longitude and latitude of the fire location, and the area occupied by the target forest can be subjected to two-dimensionalization, i.e., a two-dimensional coordinate plane is constructed, through the area size of the target forest and the corresponding specific longitude and latitude information, wherein each specific coordinate is mapped to an actual target forest site; and then the fire point position and the fire boundary of the fire disaster are subjected to coordinate transformation in the two-dimensional coordinate plane, so that the fire point position generates a specific fire point coordinate, and the fire boundary forms a boundary curve or a broken line, thereby facilitating subsequent analog simulation calculation, and simultaneously facilitating the coordinate to be changed back to longitude and latitude to obtain actual site information, thereby facilitating fire rescuers such as firemen to carry out actual fire suppression.
In an optional embodiment, the combustion parameters include parameters such as wind speed, wind direction, dryness, precipitation, drying combustible load, drying particle density, surface area to volume ratio, combustible bed depth, combustible low heat content, total mineral content, effective mineral content, combustible water content, combustible extinguishment water content and the like at corresponding coordinate positions, which is equivalent to constructing a small database for a coordinate point, so as to provide data required by subsequent calculation for fire spread, facilitate understanding of actual conditions in the field by workers, and help the workers to provide new ideas for implementing a fire extinguishing method.
And S4, setting the prediction time length.
Specifically, in this embodiment, the setting of the predicted time duration includes the following steps: and setting the predicted time length according to the time resolution of the remote sensing satellite, wherein the predicted time length is greater than the time resolution.
In an optional embodiment, the data of the Himapari-8 satellite is selected as a remote sensing satellite data source, so that the predicted time length is set to be 10 minutes, the preset time length represents an actual fire spreading result after the preset result is 10 minutes, and the preset time length is not the required operation time length for prediction. And (3) refreshing remote sensing satellite data after the Hiwari-8 satellite carries out 10 minutes, verifying a prediction result by using the refreshed remote sensing satellite data, and correcting the prediction result so as to carry out the next fire spread simulation.
And S5, calculating the change of the fire climate information and the change of the fire fuel information within the prediction duration according to the prediction duration and by combining the overall fire spread information, the fire fuel information and the fire climate information.
Specifically, in this embodiment, the calculating, according to the predicted duration and in combination with the overall fire spread information, the fire fuel information, and the fire climate information, a change in the fire climate information and a change in the fire fuel information within the predicted duration includes: predicting the change of the fire climate information within the prediction duration according to the fire climate information; constructing a fuel loss function by using the overall fire spread information and the fire fuel information; and obtaining the change of the fire disaster fuel information in the prediction time length by combining the change of the fire disaster climate information through the fuel loss function. And the simulation accuracy of fire spreading is improved by updating the fire fuel information and the fire climate information in real time.
In an alternative embodiment, the fuel loss function satisfies the following equation:
Figure 174471DEST_PATH_IMAGE001
wherein the content of the first and second substances,
Figure 842213DEST_PATH_IMAGE002
the length of time is indicated by the length of time,
Figure 979933DEST_PATH_IMAGE003
representing two-dimensional coordinates of
Figure 973297DEST_PATH_IMAGE004
The fuel loss function of the fire at the location,
Figure 930889DEST_PATH_IMAGE005
representing two-dimensional coordinates of
Figure 917168DEST_PATH_IMAGE004
The initial total reserve of combustible material at the location,
Figure 174974DEST_PATH_IMAGE005
and at a position with two-dimensional coordinates of (x, y)The loading capacity of the drying combustible materials, the density of drying particles, the surface area-to-volume ratio, the depth of the combustible material bed and other parameters, the total mineral content, the effective mineral content and other parameters are related,
Figure 808081DEST_PATH_IMAGE006
which represents the coefficient of fuel loss as a function of time,
Figure 784127DEST_PATH_IMAGE006
related to the water content of the combustible materials in the target forest, the water content of extinguishment of the combustible materials and other parameters,
Figure 793672DEST_PATH_IMAGE007
representing two-dimensional coordinates of
Figure 905984DEST_PATH_IMAGE004
The direction of the fire at the location,
Figure 490418DEST_PATH_IMAGE008
representing two-dimensional coordinates of
Figure 422602DEST_PATH_IMAGE004
The fuel distribution vector of the fire at the location,
Figure 501416DEST_PATH_IMAGE009
representing two-dimensional coordinates of
Figure 733815DEST_PATH_IMAGE004
Fuel storage function at the location.
And S6, predicting the change of the overall fire spread information after the prediction duration according to the overall fire spread information, the change of the fire climate information and the change of the fire fuel information.
Specifically, in this embodiment, the predicting, based on the overall fire spread information, the change in the fire climate information, and the change in the fire fuel information, a change in the overall fire spread information after the predicted time period elapses, includes: constructing a fire point migration speed function by using the overall fire spread information; determining a head fire point coordinate, a tail fire point coordinate and a side fire point coordinate by combining the corresponding combustion parameters and according to the fire point coordinate and the fire boundary coordinate; selecting a plurality of coordinates on the fire boundary as auxiliary side fire point coordinates, and assigning corresponding combustion parameters to the auxiliary side fire point coordinates; taking the head fire point coordinate, the tail fire point coordinate, the side fire point coordinate and the auxiliary side fire point as a plurality of starting points, and combining the fire point migration speed function and the predicted time length to obtain the fire point coordinate after the predicted time length; according to the fire point coordinates after the predicted duration, combining the change of the fire climate information, and re-defining the fire boundary; and obtaining the change of the whole fire spreading information according to the fire point coordinates after the predicted time length and the newly defined fire boundary.
In an optional embodiment, the plurality of coordinates are selected as the auxiliary side fire coordinates on the fire boundary, and the selection of the auxiliary side fire coordinates can be performed according to the change of the fire climate information and the change of the fire fuel information, that is, the parameters such as downwind direction, combustible storage amount and the like are positively correlated with the insertion amount of the auxiliary side fire coordinates.
In an alternative embodiment, the fire migration velocity function satisfies the following equation:
Figure 708724DEST_PATH_IMAGE010
wherein the content of the first and second substances,
Figure 659362DEST_PATH_IMAGE011
representing two-dimensional coordinates of
Figure 259977DEST_PATH_IMAGE004
The speed of migration of the fire at the location,
Figure 346882DEST_PATH_IMAGE012
Figure 492692DEST_PATH_IMAGE013
indicating a fire response corrected intensity factor at a position in two-dimensional coordinates (x, y), wherein,
Figure 665047DEST_PATH_IMAGE014
the spread rate of the forest fire is shown,
Figure 85665DEST_PATH_IMAGE015
a wind speed correction coefficient is expressed, and,
Figure 10764DEST_PATH_IMAGE016
the slope correction coefficient is represented by a slope correction coefficient,
Figure 858634DEST_PATH_IMAGE017
representing two-dimensional coordinates of
Figure 252707DEST_PATH_IMAGE004
The density of the fuel at the location of the fuel,
Figure 211435DEST_PATH_IMAGE018
representing the amount of heat required to ignite a unit mass of fuel,
Figure 272932DEST_PATH_IMAGE019
representing the effective thermal coefficient of heat.
In an optional embodiment, the obtaining the coordinates of the fire point after the predicted time period by combining the fire point migration speed function and the predicted time period includes the following steps: using said fire migration velocity function
Figure 26125DEST_PATH_IMAGE020
And the predicted duration
Figure 156761DEST_PATH_IMAGE021
Multiplying to obtain a fire spreading step S,
Figure 653601DEST_PATH_IMAGE022
(ii) a Obtaining the predicted duration by using the fire spreading step S and combining the wind direction angleFire point coordinates of
Figure 569604DEST_PATH_IMAGE023
Wherein, in the step (A),
Figure 759277DEST_PATH_IMAGE024
Figure 127942DEST_PATH_IMAGE025
wherein, in the step (A),
Figure 412161DEST_PATH_IMAGE026
and the angle of the wind direction is shown, namely the included angle between the wind direction and the positive direction of the x axis of the two-dimensional coordinate system.
And S7, simulating the fire spreading condition of the target forest according to the whole fire spreading information and the change of the whole fire spreading information.
Specifically, in this embodiment, the simulating the fire spreading condition of the target forest according to the overall fire spreading information and the change of the overall fire spreading information includes: constructing an initial fire condition of a target forest by using the overall fire spreading information; and simulating the fire spreading condition of the target forest by combining the initial fire condition with the change of the overall fire spreading information.
In an optional embodiment, the simulation is visualized to make the fire situation more intuitive. Referring to fig. 2, 3 and 4, fig. 2 is a schematic view of simulation of an initial fire condition, fig. 3 is a schematic view of interpolation of coordinates of a side fire point after propagation of a fire, fig. 4 is a schematic view of comparison between simulation of an initial fire condition and a fire spreading condition, specifically, in fig. 1, a dot a1 represents coordinates of a head fire point, a dot a2 represents coordinates of a tail fire point, a square represents coordinates of the side fire point, and a curve represents a fire boundary; in fig. 2, a dot b1 represents a head fire coordinate, a dot b2 represents a tail fire coordinate, a square represents the side fire coordinate, and a triangle represents the auxiliary side fire; in fig. 3, a dot a1 represents an initial head fire coordinate, a dot a2 represents an initial tail fire coordinate, a dot b1 represents a simulated head fire coordinate, a dot b2 represents a simulated tail fire coordinate, a square represents the side fire coordinate curve, a fire boundary dot represents the head fire coordinate and the tail fire coordinate, a square represents the side fire coordinate, a triangle represents the auxiliary side fire, a curve a represents an initial fire boundary, and a curve b represents a simulated fire boundary.
The method of the invention predicts the forest fire spreading by combining the real-time fire information, fire geographical information, fire fuel information and fire climate information, namely combining and mapping meteorological factors, topographic factors and fuel factors which are closely related to the forest fire and the fire condition, thereby accurately and rapidly simulating the fire spreading condition, improving the real-time and efficient requirements of simulation of forest fire spreading, providing decision support for fire departments to put out forest fires, and simultaneously having great significance and effect on forest fire prevention and reduction.
In a second aspect, the present invention further provides a real-time forest fire spreading simulation system, which is suitable for the real-time forest fire spreading simulation method of the present invention, and includes: the information input module is used for determining a target forest; setting a prediction time length; the data acquisition module is used for acquiring real-time fire information, fire geographical information, fire fuel information and fire climate information of the target forest through a remote sensing satellite and ground communication equipment; the initial fire behavior module is used for acquiring overall fire spreading information by utilizing the real-time fire information, the fire geographic information, the fire fuel information and the fire climate information; the first fire behavior simulation module is used for calculating the change of the fire climate information and the change of the fire fuel information within the prediction duration according to the prediction duration by combining the overall fire spread information, the fire fuel information and the fire climate information; the second fire behavior simulation module is used for predicting the change of the overall fire spread information after the prediction duration according to the overall fire spread information, the change of the fire climate information and the change of the fire fuel information; and the fire simulation module simulates the fire spreading condition of the target forest according to the whole fire spreading information and the change of the whole fire spreading information. Specifically, referring to fig. 5, fig. 5 shows a connection relationship of each functional module in the real-time forest fire spreading simulation system, and such a connection relationship enables the real-time forest fire spreading simulation system provided by the present invention to have a compact structure, an accurate operation result and high efficiency, so that the present invention has a practical application value and a commercial value.
Finally, it should be noted that: the above embodiments are only used to illustrate the technical solution of the present invention, and not to limit the same; while the invention has been described in detail and with reference to the foregoing embodiments, it will be understood by those skilled in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some or all of the technical features may be equivalently replaced; such modifications and substitutions do not depart from the spirit and scope of the present invention, and they should be construed as being included in the following claims and description.

Claims (6)

1. A real-time forest fire spreading simulation method is characterized by comprising the following steps:
determining a target forest;
acquiring real-time fire information, fire geographical information, fire fuel information and fire climate information of the target forest through a remote sensing satellite and ground communication equipment;
acquiring overall fire spreading information by using the real-time fire information, the fire geographical information, the fire fuel information and the fire climate information;
setting a prediction time length;
calculating the change of the fire climate information and the change of the fire fuel information within the prediction duration according to the prediction duration by combining the overall fire spread information, the fire fuel information and the fire climate information;
predicting the change of the overall fire spreading information after the prediction duration according to the overall fire spreading information, the change of the fire climate information and the change of the fire fuel information;
simulating the fire spreading condition of the target forest according to the overall fire spreading information and the change of the overall fire spreading information;
the method for acquiring the overall fire spreading information by utilizing the real-time fire information, the fire geographic information, the fire fuel information and the fire climate information comprises the following steps:
constructing a two-dimensional coordinate interface of the target forest by using the fire geographic information;
setting fire point coordinates and fire boundary coordinates in the two-dimensional coordinate interface according to the real-time fire information;
assigning corresponding combustion parameters to corresponding fire point coordinates and fire boundary coordinates according to the fire fuel information and the fire climate information;
acquiring overall fire spreading information by setting the fire point coordinate and the fire boundary coordinate and assigning a two-dimensional coordinate interface;
the setting of the predicted time length comprises the following steps:
setting the predicted time length according to the time resolution of the remote sensing satellite, wherein the predicted time length is greater than the time resolution;
the step of calculating the change of the fire climate information and the change of the fire fuel information within the prediction duration by combining the overall fire spread information, the fire fuel information and the fire climate information according to the prediction duration comprises the following steps:
predicting the change of the fire climate information within the prediction duration according to the fire climate information;
constructing a fuel loss function by using the overall fire spread information and the fire fuel information;
obtaining the change of the fire fuel information in the predicted time length by combining the change of the fire climate information through the fuel loss function;
the predicting the change of the overall fire spread information after the prediction duration according to the overall fire spread information, the change of the fire climate information and the change of the fire fuel information comprises the following steps:
constructing a fire point migration speed function by using the overall fire spread information;
determining a head fire point coordinate, a tail fire point coordinate and a side fire point coordinate by combining the corresponding combustion parameters and through the fire point coordinate and the fire boundary coordinate;
selecting a plurality of coordinates on the fire boundary as auxiliary side fire point coordinates according to the change of fire climate information and the change of fire fuel information, and assigning corresponding combustion parameters to the auxiliary side fire point coordinates;
taking the head fire point coordinate, the tail fire point coordinate, the side fire point coordinate and the auxiliary side fire point as a plurality of starting points, and combining the fire point migration speed function and the predicted time length to obtain the fire point coordinate after the predicted time length;
according to the fire point coordinates after the predicted duration, combining the change of the fire climate information, and re-demarcating a fire boundary;
and re-defining a fire boundary according to the fire point coordinates after the predicted time length, and obtaining the change of the overall fire spreading information.
2. The real-time forest fire spread simulation method according to claim 1, wherein the obtaining of the real-time fire information, the fire geographical information, the fire fuel information and the fire climate information of the target forest by a remote sensing satellite and a ground communication device comprises the following steps:
acquiring fire point information and a fire boundary of the target forest through the remote sensing satellite;
the fire point information and the fire boundary are arranged, and the real-time fire information is obtained;
acquiring the gradient, the slope direction and DEM data of the target forest through the ground communication equipment;
combining the slope, the slope direction and the DEM data to obtain the fire disaster geographic information;
acquiring the loading capacity, the density of dried particles, the surface area-to-volume ratio, the depth of a combustible bed, the low heat content of the combustible, the total mineral content, the effective mineral content, the water content of the combustible and the extinguishing water content of the combustible in the target forest by using the ground communication equipment;
converging said dry combustible load, said dry particle density, said surface area to volume ratio, said combustible bed depth, said combustible low heat content, said total mineral content, said effective mineral content, said combustible water content, and said combustible quench water content to obtain said fire fuel information;
acquiring the wind speed, the wind direction, the dryness and the precipitation of the target forest by using the ground communication equipment;
and acquiring the fire climate information through the wind speed, the wind direction, the dryness and the precipitation.
3. A real-time forest fire spread simulation method according to claim 1, wherein the fuel loss function satisfies the following formula:
Figure 118527DEST_PATH_IMAGE001
wherein the content of the first and second substances,
Figure 629142DEST_PATH_IMAGE002
the length of time is indicated by the length of time,
Figure 143300DEST_PATH_IMAGE003
representing two-dimensional coordinates of
Figure 7351DEST_PATH_IMAGE004
The fuel loss function of the fire at the location,
Figure 680778DEST_PATH_IMAGE005
representing two-dimensional coordinates of
Figure 972082DEST_PATH_IMAGE004
The initial total reserve of combustible material at the location,
Figure 973536DEST_PATH_IMAGE006
the coefficient of fuel loss is expressed as,
Figure 765912DEST_PATH_IMAGE007
representing two-dimensional coordinates of
Figure 169211DEST_PATH_IMAGE004
The direction of the fire at the location of the fire,
Figure 896996DEST_PATH_IMAGE008
representing two-dimensional coordinates of
Figure 244800DEST_PATH_IMAGE004
The fuel distribution vector of the fire at the location,
Figure 716233DEST_PATH_IMAGE009
representing two-dimensional coordinates of
Figure 836023DEST_PATH_IMAGE004
Fuel storage function at the location.
4. A real-time forest fire spread simulation method according to claim 1, wherein the fire migration velocity function satisfies the following formula:
Figure 288DEST_PATH_IMAGE010
wherein the content of the first and second substances,
Figure 710755DEST_PATH_IMAGE011
representing two-dimensional coordinates of
Figure 844934DEST_PATH_IMAGE004
The speed of migration of the fire at the location,
Figure 222825DEST_PATH_IMAGE012
Figure 26833DEST_PATH_IMAGE013
indicating a fire response corrected intensity factor at a location in two-dimensional coordinates (x, y), wherein,
Figure 614809DEST_PATH_IMAGE014
the spread rate of the forest fire is shown,
Figure 162465DEST_PATH_IMAGE015
the wind speed correction coefficient is represented by a coefficient,
Figure 394864DEST_PATH_IMAGE016
the slope correction coefficient is represented by a slope correction coefficient,
Figure 759986DEST_PATH_IMAGE017
representing two-dimensional coordinates of
Figure 445045DEST_PATH_IMAGE004
The density of the fuel at the location of the fuel,
Figure 61971DEST_PATH_IMAGE018
representing the amount of heat required to ignite a unit mass of fuel,
Figure 273510DEST_PATH_IMAGE019
representing the effective thermal coefficient of heat.
5. A real-time forest fire spread simulation method according to claim 1, wherein the simulating the fire spread condition of the target forest according to the overall fire spread information and the change of the overall fire spread information comprises the following steps:
constructing an initial fire condition of a target forest by using the overall fire spreading information;
and simulating the fire spreading condition of the target forest by combining the initial fire condition and the change of the overall fire spreading information.
6. A real-time forest fire spread simulation system, the real-time forest fire spread simulation system being adapted to the real-time forest fire spread simulation method of any one of claims 1 to 5, comprising:
the information input module is used for determining a target forest; setting a prediction time length;
the data acquisition module is used for acquiring real-time fire information, fire geographical information, fire fuel information and fire climate information of the target forest through a remote sensing satellite and ground communication equipment;
the initial fire behavior module is used for acquiring overall fire spread information by utilizing the real-time fire information, the fire geographic information, the fire fuel information and the fire climate information;
the first fire behavior simulation module is used for calculating the change of the fire climate information and the change of the fire fuel information within the prediction duration according to the prediction duration by combining the overall fire spread information, the fire fuel information and the fire climate information;
the second fire behavior simulation module is used for predicting the change of the overall fire spread information after the prediction duration according to the overall fire spread information, the change of the fire climate information and the change of the fire fuel information;
and the fire simulation module simulates the fire spreading condition of the target forest according to the whole fire spreading information and the change of the whole fire spreading information.
CN202211022678.0A 2022-08-25 2022-08-25 Real-time forest fire spreading simulation method and system Active CN115099073B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211022678.0A CN115099073B (en) 2022-08-25 2022-08-25 Real-time forest fire spreading simulation method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211022678.0A CN115099073B (en) 2022-08-25 2022-08-25 Real-time forest fire spreading simulation method and system

Publications (2)

Publication Number Publication Date
CN115099073A CN115099073A (en) 2022-09-23
CN115099073B true CN115099073B (en) 2022-11-08

Family

ID=83300722

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211022678.0A Active CN115099073B (en) 2022-08-25 2022-08-25 Real-time forest fire spreading simulation method and system

Country Status (1)

Country Link
CN (1) CN115099073B (en)

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101090266B1 (en) * 2009-01-28 2011-12-07 대한민국(관리부서 : 산림청 국립산림과학원장) Forecasting System and method of Forest Fire Spreading
CN103020369A (en) * 2012-12-21 2013-04-03 浙江农林大学 High-resolution forest fire forecasting method
CN103065045A (en) * 2012-12-21 2013-04-24 浙江农林大学 Forecast method of forest fire danger of coupling climate and vegetation
US11966670B2 (en) * 2018-09-06 2024-04-23 Terrafuse, Inc. Method and system for predicting wildfire hazard and spread at multiple time scales

Also Published As

Publication number Publication date
CN115099073A (en) 2022-09-23

Similar Documents

Publication Publication Date Title
US11660480B2 (en) Fire forecasting
Mandel et al. Recent advances and applications of WRF–SFIRE
De Meij et al. Impact of SRTM and Corine Land Cover data on meteorological parameters using WRF
US20220016455A1 (en) Fire monitoring
Zhao A study of AR-, TS-, and MCS-associated precipitation and extreme precipitation in present and warmer climates
JP6006902B1 (en) Fire simulation apparatus, fire simulation method, and fire simulation program
Li et al. Projecting species loss and turnover under climate change for 111 Chinese tree species
JP5922570B2 (en) Diffusion state prediction apparatus, method and program
Schunk et al. Forest fire danger rating in complex topography–results from a case study in the Bavarian Alps in autumn 2011
Rahman et al. Forest fire occurrence and modeling in Southeastern Australia
JP4880440B2 (en) Snow accretion prediction method and snow accretion prediction program
Mcmorrow et al. Moorland wildfire risk, visitors and climate change: patterns, prevention and policy
CN115099073B (en) Real-time forest fire spreading simulation method and system
Velasquez et al. Wildfire propagation simulation tool using cellular automata and GIS
Patnaik et al. Large scale urban simulations with FCT
Dobrinkova et al. WRF-fire applied in Bulgaria
Di Leo et al. Dynamic data driven ensemble for wildfire behaviour assessment: A case study
Mu et al. On the dynamics and air-quality impact of the exceptional East Asian dust outbreak in mid-March 2021
Mandel et al. New features in WRF-SFIRE and the wildfire forecasting and danger system in Israel
Butler et al. High resolution wind direction and speed information for support of fire operations
CN113139272A (en) Forest fire spreading prediction method, device, equipment and storage medium
Kumar Mapping and modeling of fires in the wildland-urban interface
KR102663384B1 (en) Method for analyzing effect of hygroscopic seeding material sprayed on ground aerosol concentration through airborne cloud seeding experiment
Majlingová et al. Evaluation of dynamic modelling applications to support the disaster risk management at local level
KR102333524B1 (en) Vegetation cooling model calculation system and method to analyze the effect of temperature reduction according to vegetation characteristics

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