CN114896762A - Forest fire spreading speed calculation method based on trial algorithm - Google Patents
Forest fire spreading speed calculation method based on trial algorithm Download PDFInfo
- Publication number
- CN114896762A CN114896762A CN202210399850.8A CN202210399850A CN114896762A CN 114896762 A CN114896762 A CN 114896762A CN 202210399850 A CN202210399850 A CN 202210399850A CN 114896762 A CN114896762 A CN 114896762A
- Authority
- CN
- China
- Prior art keywords
- fire
- fire point
- acquiring
- spreading
- 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.)
- Pending
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 25
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 12
- 238000000034 method Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 abstract description 3
- 238000004088 simulation Methods 0.000 abstract description 2
- 241000282414 Homo sapiens Species 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005469 granulation Methods 0.000 description 1
- 230000003179 granulation Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
Images
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
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
-
- 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
-
- 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)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Tourism & Hospitality (AREA)
- General Engineering & Computer Science (AREA)
- Marketing (AREA)
- Databases & Information Systems (AREA)
- Development Economics (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Remote Sensing (AREA)
- Data Mining & Analysis (AREA)
- Game Theory and Decision Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Educational Administration (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a forest fire spreading speed calculation method based on a trial algorithm, which comprises the following steps of: s1, acquiring the fire point (fx) position in the forest area; s2, reading a GIS (geographic information System) space database according to the fire point (fx) position, and meshing the influence factors of the forest area; s3, acquiring grid data (gr) of the fire point (fx); s4, acquiring the spreading rate (v) of the fire point (fx) in the direction (d) according to the grid data (gr); s5, setting the calculation interval as unit time t, and calculating the time (tv) when the fire point (fx) reaches the boundary of the grid (gr) in the direction (d) at the speed (v); s6, repeating the steps S2-S4 until the time (tv) is less than or equal to the unit time (t); s7, acquiring the position of the fire point (fx) in the unit time (t) and the direction d. The invention alternately calculates the spreading speed in different grids by a trial algorithm, and the calculation result is more accurate, thereby achieving better simulation effect.
Description
Technical Field
The invention relates to the technical field of forest fire spreading, in particular to a forest fire spreading speed calculation method based on a trial algorithm.
Background
Forest plays more and more important role in the development of human life and economic society, and forest fire plays a significant factor in the damage to the forest. In order to protect the forest, people need to study the spreading rule of the forest fire, so as to reduce the loss of the forest fire to human beings when the forest fire occurs.
The prior art has the defects that factor information is too general, influence factors are induced by methods such as block averaging and uniformization, geographic factor mining is not detailed enough, and regularization results are obtained when the fire speed and the fire are calculated.
Disclosure of Invention
The present invention is directed to solving, at least to some extent, one of the technical problems in the related art. Therefore, the invention provides a forest fire spreading speed calculation method based on a trial algorithm, which is characterized in that the spreading speed and the spreading time of an individual grid are accumulated to obtain the running track of a fire point in a grid GIS space database in unit time and predict the fire behavior.
The technical scheme adopted by the invention is as follows:
a forest fire spreading speed calculation method based on a trial algorithm comprises the following steps: s1, finding out the fire, and acquiring the position of a fire point (fx) in the forest area; s2, reading a GIS (geographic information System) space database according to the fire point (fx) position, and meshing the influence factors of the forest area; s3, acquiring grid data (gr) of the fire point (fx), wherein the grid data (gr) comprise terrain gradient, combustible type, combustible humidity and wind speed; s4, acquiring the spreading rate (v) of the fire point (fx) according to the grid data (gr) in the directions (d), wherein the directions (d) comprise east, south, west, north, southeast, northeast, southwest and northwest directions; s5, setting a calculation interval as a unit time t, and calculating the time (tv) when the fire point (fx) reaches the boundary of the grid (gr) in the direction (d) at the speed (v); s6, repeating the steps S2-S4 until the time (tv) is less than or equal to the unit time (t); s7, acquiring the position of the fire point (fx) in the unit time (t) and the direction d.
Further, the spreading speed is calculated by adopting a Wangzheng non-model, and the calculation formula is as follows: v is V 0 *K s *K w *K φ Wherein V is the propagation rate, V 0 To initial propagation rate, K s Correction factor for combustible type, K w As a factor of wind correction, K φ Correcting the coefficient for the terrain gradient; v 0 The calculation formula is as follows: v 0 ＝1.0372*e -0.057*m Wherein m is combustible humidity.
The invention has the beneficial effects that:
the method fully utilizes the GIS spatial database to simulate factors which seriously affect the spreading speed, such as complex and changeable combustible information, terrain gradient information and the like in the forest so as to achieve data granulation, alternately calculates the spreading speed in different grids through trial algorithm, and realizes more accurate calculation of the spreading speed of the area so as to achieve better simulation effect.
Drawings
Fig. 1 is a schematic flow chart of an embodiment of a forest fire spreading speed calculation method based on a trial algorithm.
Detailed Description
To facilitate an understanding of the invention, the invention will now be described more fully hereinafter with reference to the accompanying drawings. Preferred embodiments of the present invention are shown in the drawings. This invention may, however, be embodied in many different forms and should not be construed as limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete.
Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. The terminology used in the description of the invention herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. As used herein, the term "and/or" includes any and all combinations of one or more of the associated listed items.
Referring to fig. 1, fig. 1 is a schematic flow chart of a forest fire spread rate calculation method based on a trial algorithm according to an embodiment of the present invention. As shown in fig. 1, the method comprises the steps of:
s1, finding out a fire, and acquiring the position of a fire point (fx) in the forest area, wherein the fire point is the position of the fire point corresponding to the electronic map;
s2, reading a GIS (geographic information system) spatial database according to the position of the fire point (fx), and meshing the influence factors of the forest area according to a specified scale for the geographical area taking the fire point (fx) as the center;
s3, respectively generating a spreading point in the east, south, west, north, southeast, northeast, southwest and northwest directions by taking the position of the fire point (fx) as the center of a circle, and respectively obtaining eight spreading pointsGrid data (gr) of the location of the propagation point, wherein the grid data (gr) comprises terrain slope (Pi), combustible type, combustible humidity m, wind speed V Wind power ；
Preferably, the terrain slope (Pi) ranges between-60 and 60 degrees.
S4, calculating the spreading rate of all spreading points, and calculating the spreading rate by adopting a Wangzheng non-model, wherein the calculation formula is as follows:
V＝V 0 *K s *K w *K φ
wherein V is the propagation rate, V 0 To initial propagation rate, K s Correction factor for combustible type, K w As a factor of wind correction, K φ Correcting the coefficient for the terrain gradient;
V 0 the calculation formula is as follows:
V 0 ＝1.0372*e -0.057*m
wherein m is the combustible humidity, and the percentage value is taken;
when wind speed V Wind power When changed, the wind power correction coefficient K w When the terrain slope Pi changes along with the change, the terrain slope correction coefficient K φ Changes also occur; specifically, the wind force correction coefficient K w Correction coefficient K of terrain gradient φ And the specific calculation formula is within the understanding of the person skilled in the art, and it is not expanded here that the correction factor K for combustible types s The determination is made based on the vegetation type in the database, with river set to 0 and building set to 0.05.
S5, calculating the time (tv) for the fire point (fx) to reach the boundary of the grid (gr) in the direction (d) at the speed (v), wherein the calculation formula is as follows:
tv-grid assigned dimension/V
S6, setting the calculation interval as unit time t, repeating the steps S2-S4 until the time (tv) is less than or equal to the unit time (t), and updating the spreading rate V' of each spreading point;
and S7, acquiring the positions of the fire points (fx) in the unit time (t) and the direction d, wherein the spread distance of each spread point in the unit time (t) is Delta S-V-t.
Preferably, a section of circular arc line is drawn at the current spreading position of each spreading point, and adjacent circular arc lines are connected to form a wavy irregular circle integrally so as to simulate the forest fire spreading effect. Wherein, the arc line of each segment uses the ignition point as the center of a circle, and the spreading distance of each trigger point is the radius to draw an arc.
In this embodiment, eight spread points are generated from the initial fire point position (f1) and spread in eight directions, after a certain calculation period, the end point spread position of each spread point is set as a new fire point position (f2), and new eight spread points are generated and spread, and after a certain calculation period, the end point spread position of each new spread point is set as a new fire point position (f3), and a new spread point is generated and spread … ….
While the preferred embodiments of the present invention have been illustrated and described, it will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the spirit and scope of the invention as defined by the appended claims.
Claims (2)
1. A forest fire spreading speed calculation method based on a trial algorithm is characterized by comprising the following steps:
s1, finding out the fire, and acquiring the position of a fire point (fx) in the forest area;
s2, reading a GIS (geographic information System) space database according to the fire point (fx) position, and meshing the influence factors of the forest area;
s3, acquiring grid data (gr) of the fire point (fx), wherein the grid data (gr) comprise terrain gradient, combustible type, combustible humidity and wind speed;
s4, acquiring the spreading rate (v) of the fire point (fx) according to the grid data (gr) in the directions (d), wherein the directions (d) comprise east, south, west, north, southeast, northeast, southwest and northwest directions;
s5, setting a calculation interval as a unit time t, and calculating the time (tv) when the fire point (fx) reaches the boundary of the grid (gr) in the direction (d) at the speed (v);
s6, repeating the steps S2-S4 until the time (tv) is less than or equal to the unit time (t);
s7, acquiring the position of the fire point (fx) in the unit time (t) and the direction d.
2. The forest fire spreading speed calculating method based on the trial algorithm as claimed in claim 1, wherein the spreading speed is calculated by adopting a Wangzhong non-model, and the calculation formula is as follows: v is V 0 *K s *K w *K φ Wherein V is the propagation rate, V 0 To initial propagation rate, K s Correction factor for combustible type, K w As a factor of wind correction, K φ Correcting the coefficient for the terrain gradient; v 0 The calculation formula is as follows: v 0 ＝1.0372*e -0.057*m Wherein m is combustible humidity.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210399850.8A CN114896762A (en) | 2022-04-15 | 2022-04-15 | Forest fire spreading speed calculation method based on trial algorithm |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210399850.8A CN114896762A (en) | 2022-04-15 | 2022-04-15 | Forest fire spreading speed calculation method based on trial algorithm |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114896762A true CN114896762A (en) | 2022-08-12 |
Family
ID=82716814
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210399850.8A Pending CN114896762A (en) | 2022-04-15 | 2022-04-15 | Forest fire spreading speed calculation method based on trial algorithm |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114896762A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115512510A (en) * | 2022-10-11 | 2022-12-23 | 深圳市狮子王科技有限公司 | Charging pile fire intelligent processing system and method |
-
2022
- 2022-04-15 CN CN202210399850.8A patent/CN114896762A/en active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115512510A (en) * | 2022-10-11 | 2022-12-23 | 深圳市狮子王科技有限公司 | Charging pile fire intelligent processing system and method |
| CN115512510B (en) * | 2022-10-11 | 2023-09-29 | 深圳市狮子王科技有限公司 | Intelligent fire disaster treatment system and method for charging pile |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Kitoh et al. | First super-high-resolution model projection that the ancient “Fertile Crescent” will disappear in this century | |
| CN109916788B (en) | PM pair for distinguishing emission change and meteorological condition change of different areas2.5Method for influencing concentration | |
| Mohr | Interannual, monthly, and regional variability in the wet season diurnal cycle of precipitation in sub-Saharan Africa | |
| Phan et al. | Extreme climatic events over Vietnam from observational data and RegCM3 projections | |
| Earley-Spadoni | Landscapes of warfare: Intervisibility analysis of Early Iron and Urartian fire beacon stations (Armenia) | |
| CN108280849B (en) | Prediction correction and leakage rate estimation method for gas leakage concentration field of comprehensive pipe gallery | |
| Kadioğlu et al. | Variations and trends in Turkish seasonal heating and cooling degree-days | |
| Xu et al. | Projected changes in climate extremes over China in the 21st century from a high resolution regional climate model (RegCM3) | |
| CN107704641A (en) | Fine simulation of wind method based on outdoor scene vegetation spatial distribution roughness | |
| CN114781169B (en) | Forest fire spreading simulation method based on cellular automaton | |
| CN103292784B (en) | Underwater surface topography visuality analyzing method based on sound ray track | |
| CN114896762A (en) | Forest fire spreading speed calculation method based on trial algorithm | |
| Pierce et al. | Creating climate projections to support the 4th California climate assessment | |
| Jeong et al. | Influences of atmospheric blocking on North American summer heatwaves in a changing climate: a comparison of two Canadian Earth system model large ensembles | |
| Guo et al. | Irrigated cropland expansion exacerbates the urban moist heat stress in northern India | |
| Photiadou et al. | Modeling European hot spells using extreme value analysis | |
| Sono et al. | Spatiotemporal evolution of West Africa’s urban landscape characteristics applying harmonized DMSP-OLS and NPP-VIIRS nighttime light (NTL) data | |
| Zhang et al. | Reconstruction of cropland spatial patterns and its spatiotemporal changes over the 20th century on the Songnen Plain, Northeast China | |
| Kim et al. | Effects of the density of meteorological observations on the diagnostic wind fields and the performance of photochemical modeling in the greater Seoul area | |
| Barros et al. | Urban-biased trends in Buenos Aires’ mean temperature | |
| CN113642215A (en) | Simulation method for gridding forest fire spreading trend in any direction | |
| CN109061719B (en) | Method for predicting earthquake | |
| Tošić et al. | Landslide susceptibility zonation: a case study of the Municipality of Banja Luka (Bosnia and Herzegovina) | |
| Sobouti | The relationship between unique geomagnetic and auroral events | |
| Van der Kamp et al. | Evaluation of the Monthly Drought Code as a metric for fire weather in a region of complex terrain and un-certainties in future projections |
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 |