CN114297953B - Mountain area storm disaster risk division method - Google Patents
Mountain area storm disaster risk division method Download PDFInfo
- Publication number
- CN114297953B CN114297953B CN202111599073.3A CN202111599073A CN114297953B CN 114297953 B CN114297953 B CN 114297953B CN 202111599073 A CN202111599073 A CN 202111599073A CN 114297953 B CN114297953 B CN 114297953B
- Authority
- CN
- China
- Prior art keywords
- mountain
- grid
- wind
- risk
- wind speed
- 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 18
- 239000012530 fluid Substances 0.000 claims abstract description 8
- 238000004088 simulation Methods 0.000 claims abstract description 5
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000004927 fusion Effects 0.000 claims description 2
- 230000007547 defect Effects 0.000 abstract description 3
- 230000002265 prevention Effects 0.000 abstract description 2
- 238000012876 topography Methods 0.000 abstract description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Landscapes
- Alarm Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a method for dividing the risk of a mountain storm disaster, belonging to the field of weather disaster prevention and reduction. According to the invention, the simulation of the mountain projection surface wind field is performed by using fluid dynamics CFD simulation software, and the wind speed of each grid intersecting the intersecting line of the mountain projection surface and the ground surface is obtained; reuse ofAnd calculating to obtain a global wind field of the mountain area, and finally, carrying out corresponding risk division according to the established wind disaster risk division index system. The technical method effectively overcomes the defects of the existing method for dividing the risk of the storm disaster in the mountain area, fully considers the influence of the topography of the mountain area on wind field simulation, and reasonably realizes the division of the risk of the storm disaster in the mountain area. The method of the invention is convenient to use and easy to realize.
Description
Technical Field
The invention relates to the field of weather disaster prevention and reduction, in particular to a method for dividing the risk of a mountain disaster.
Background
A digital elevation model (Digital Elevation Model), abbreviated as DEM, is a physical ground model for realizing the digital simulation of ground terrain through limited terrain elevation data, which is a set of ordinal value array form for representing ground elevation, the data can be obtained through a public channel, CFD (Computational Fluid Dynamics) is the intersection science of modern hydrodynamics, numerical mathematics and computer combination. The integral and differential terms in the fluid mechanics control equation are approximately expressed as discrete algebraic forms, so that the discrete algebraic equation sets are formed, then the discrete algebraic equation sets are solved through a computer, a numerical solution on discrete time/space points is obtained, and CFD (computational fluid dynamics) related calculation is completed through existing public software.
The current high wind risk zone is mainly divided aiming at plain areas, the research on the high wind risk zone in mountain areas is less, and the research on the high wind risk zone in mountain areas is difficult and inaccurate mainly because the wind fields in mountain areas are greatly influenced by terrain factors and weather sites in mountain areas are less and unevenly distributed.
Disclosure of Invention
In order to overcome the defects of the existing mountain area high wind disaster risk division method, the invention provides a mountain area high wind disaster risk division method.
The method comprises the following steps:
Step 1, establishing a digital elevation grid model DEM of mountain terrain;
Step 2, dividing mountain boundary lines and mountain ridge lines from mountain bodies in the digital elevation grid model;
Step 3, obtaining the wind speed U αi of each grid at the mountain boundary line; where αi represents the wind direction of each grid at the mountain boundary line;
Step 4, obtaining the wind speed of each grid in the boundary line of the mountain Wherein H aij is the elevation of the jth grid consistent with the grid wind direction at the mountain boundary, H ai is the elevation of the grid at the mountain boundary with the wind direction of αi, λ is the proportionality coefficient, and for H aij-Hai >10 meters, λ takes 3.2; for h aij-Hai to be less than or equal to 10 meters, lambda is 1.4, and wind speed calculation on the surface of the mountain is completed;
Step 5: wind speed calculation is carried out on the mountain surface in the whole digital elevation model DEM, and a mountain terrain wind speed distribution map is obtained;
Step 6: establishing a large wind disaster risk division index system:
Let v be the wind speed (unit meter/second)
A:5<V is less than or equal to 10: a low risk area;
b: v is 10< 15 or less: a risk area;
c:15< v: a high risk area;
step 7: mountain area major wind disaster risk division:
comparing the wind speed of each grid in the digital elevation model DEM with the index system of the step 6 to determine the risk level; and sequentially completing grid fusion with consistent risk levels to form high-wind disaster risk areas with different levels, and realizing regional high-wind disaster risk areas.
Preferably, the wind speed U αi at the mountain boundary line is obtained by simulating the mountain projection surface wind field by using fluid dynamics CFD simulation software.
Preferably, in the step 4, the grids on two sides of the mountain ridge line respectively calculate the wind speed by using the wind speed of each grid where the projection plane of the mountain body on one side intersects with the intersecting line of the ground surface.
Preferably, in step 4, the grids repeatedly calculated are arranged on two sides of the ridge line, and the average value of the wind speed obtained by repeated calculation is taken.
The beneficial effects are that: the method effectively overcomes the defects of the existing method for dividing the risk of the storm disaster in the mountain area, fully considers the influence of the topography of the mountain area on wind field simulation, and reasonably realizes the division of the risk of the storm disaster in the mountain area. The method of the invention is convenient to use and easy to realize.
Drawings
FIG. 1 is a partial digital elevation grid model of the present invention;
FIG. 2 is a schematic view of the intersection of the projection plane of the mountain and the ground surface according to the present invention;
FIG. 3 is a grid wind speed distribution intersecting an intersection line as simulated by the fluid dynamics CFD simulation software of the present invention;
FIG. 4 is a graph showing the wind velocity distribution of each grid within the boundary line of the present invention;
fig. 5 is a schematic diagram of a mountain wind disaster risk zone according to the present invention.
Detailed Description
The invention will be further explained with reference to the accompanying drawings, wherein a part of the digital elevation grid model of a mountain terrain is selected for detailed explanation.
Please refer to fig. 1: the area is a typical mountain terrain which is a 5 row 5 column grid pattern, where the numerical value of each grid represents the elevation of the corresponding location in meters.
Please refer to fig. 2: an intersecting grid can be determined from the intersection line (black bolded representation) of the projection surface and the ground surface of the mountain, and a ridge line (broken line representation) is also divided.
Please refer to fig. 3: wind speeds at the intersecting grids are simulated by using fluid dynamics CFD simulation software, and grid values at the moment represent the wind speeds at the mountain boundaries.
Please refer to fig. 4: obtaining wind speed of each grid inside boundary line of mountainWherein H aij is the elevation of the jth grid consistent with the grid wind direction at the mountain boundary, H ai is the elevation of the grid at the mountain boundary with the wind direction of αi, λ is the proportionality coefficient, and for H aij-Hai >10 meters, λ takes 3.2; for h aij-Hai to be less than or equal to 10 meters, lambda is 1.4, and wind speed calculation on the surface of the mountain is completed;
Please refer to fig. 5: and comparing the mountain area wind speed distribution with a mountain area major wind disaster risk division index system, thereby determining the risk level. For simplicity of explanation, high, medium, and low represent high risk areas, medium risk areas, and low risk areas of a large wind disaster in this order.
The foregoing is merely a preferred embodiment of the present invention and it should be noted that modifications and adaptations to those skilled in the art may be made without departing from the principles of the present invention, which are intended to be comprehended within the scope of the present invention.
Claims (4)
1. A method for distinguishing a mountain area from a major wind disaster risk, the method comprising the steps of:
Step 1, establishing a digital elevation grid model DEM of mountain terrain;
step 2, dividing mountain boundary lines and ridge lines for mountain bodies in the digital elevation grid model DEM;
step 3, obtaining the wind speed U αi of each grid at the mountain boundary line; wherein αi represents the wind direction of each grid at the mountain boundary line;
step 4, obtaining the wind speed of each grid in the mountain boundary line Wherein H aij is the elevation of the jth grid consistent with the grid wind direction at the mountain boundary line, H ai is the elevation of the grid at the mountain boundary line with the wind direction of alpha i, lambda is a proportionality coefficient, and for H aij-Hai >10 m, lambda takes 3.2; for h aij-Hai to be less than or equal to 10 meters, lambda is 1.4, and wind speed calculation of the mountain surface is completed;
step 5: wind speed calculation is carried out on the mountain surface in the whole digital elevation model DEM, and a mountain terrain wind speed distribution map is obtained;
Step 6: establishing a large wind disaster risk division index system:
a:5<V is less than or equal to 10: a low risk area;
b: v is 10< 15 or less: a risk area;
c:15< v: a high risk area; wherein v is wind speed;
Step 7: mountain area major wind disaster risk division: comparing the wind speed of each grid in the digital elevation model DEM with the index system to determine a risk level; and sequentially completing grid fusion with consistent risk levels, forming a disaster risk area of each level, and establishing a regional disaster risk area.
2. The method for distinguishing a mountain storm disaster risk according to claim 1, wherein the wind speed U αi of each grid at the boundary line of the mountain is obtained by performing a mountain projection surface wind field simulation by using fluid dynamics CFD simulation software.
3. The mountain area high wind disaster risk zone method of claim 1 wherein the grids on the left side of the ridge line perform wind speed calculation using wind speed of each grid at the mountain boundary line on the left side; the grid on the right side of the ridge line performs wind speed calculation using the wind speed of each grid at the boundary line of the mountain on the right side.
4. The method for distinguishing a risk of a mountain storm disaster as claimed in claim 1 wherein said grid wind speed at said ridge line is an average of wind speeds obtained by repeating calculation.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111599073.3A CN114297953B (en) | 2021-12-24 | 2021-12-24 | Mountain area storm disaster risk division method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111599073.3A CN114297953B (en) | 2021-12-24 | 2021-12-24 | Mountain area storm disaster risk division method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114297953A CN114297953A (en) | 2022-04-08 |
| CN114297953B true CN114297953B (en) | 2024-04-30 |
Family
ID=80968961
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111599073.3A Active CN114297953B (en) | 2021-12-24 | 2021-12-24 | Mountain area storm disaster risk division method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114297953B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103177301A (en) * | 2013-03-12 | 2013-06-26 | 南京信息工程大学 | Typhoon disaster risk estimate method |
| CN103455711A (en) * | 2013-08-15 | 2013-12-18 | 广州地理研究所 | Small watershed region-oriented landslide hazard risk division method based on mechanism analysis |
| CN110570107A (en) * | 2019-08-28 | 2019-12-13 | 浙江仁欣环科院有限责任公司 | mountain torrent disaster risk assessment method based on DEM |
| CN113807740A (en) * | 2021-09-30 | 2021-12-17 | 上海交通大学 | Risk assessment system for hydrological disasters of water burst runoff of mountain railway tunnel |
-
2021
- 2021-12-24 CN CN202111599073.3A patent/CN114297953B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103177301A (en) * | 2013-03-12 | 2013-06-26 | 南京信息工程大学 | Typhoon disaster risk estimate method |
| CN103455711A (en) * | 2013-08-15 | 2013-12-18 | 广州地理研究所 | Small watershed region-oriented landslide hazard risk division method based on mechanism analysis |
| CN110570107A (en) * | 2019-08-28 | 2019-12-13 | 浙江仁欣环科院有限责任公司 | mountain torrent disaster risk assessment method based on DEM |
| CN113807740A (en) * | 2021-09-30 | 2021-12-17 | 上海交通大学 | Risk assessment system for hydrological disasters of water burst runoff of mountain railway tunnel |
Non-Patent Citations (2)
| Title |
|---|
| 基于DEM的山洪灾害风险区划分研究;许小华;何雯;;中国农村水利水电;20151015(第10期);全文 * |
| 渭南市设施农业大风灾害风险区划研究;张永红;葛徽衍;韩蓓蓓;刘红;;陕西气象;20160115(第01期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114297953A (en) | 2022-04-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103363962B (en) | Remote sensing evaluation method of lake water reserves based on multispectral images | |
| Yamazaki et al. | Deriving a global river network map and its sub-grid topographic characteristics from a fine-resolution flow direction map | |
| CN108021780B (en) | Mountain torrent dynamic simulation method based on irregular unstructured grid model | |
| CN114997541B (en) | Urban waterlogging prediction method and early warning platform based on digital twin technology | |
| CN110009037B (en) | Short-term engineering wind speed prediction method and system based on physical information coupling | |
| CN110232471B (en) | Rainfall sensor network node layout optimization method and device | |
| CN111507375B (en) | Urban waterlogging risk rapid assessment method and system | |
| WO2022077767A1 (en) | Traffic flow prediction method and apparatus, computer device, and readable storage medium | |
| CN103279588B (en) | Steel bridge deck Fatigue stress calculation method under vehicular load and temperature acting in conjunction | |
| CN111767597A (en) | City model verification method, device, equipment and storage medium | |
| CN109978275B (en) | Extreme strong wind speed prediction method and system based on mixed CFD and deep learning | |
| CN114970302B (en) | Regional groundwater condition prediction method based on groundwater monitoring system | |
| CN115329691A (en) | CFD and GIS-based super-large city wind environment simulation method | |
| CN115115262A (en) | Flood risk disaster assessment method | |
| CN111177917B (en) | Slope length extraction method based on SRTM | |
| CN115081360A (en) | Wind power plant wake flow evaluation method and device based on simple actuating disc model | |
| CN115471679A (en) | Method and intelligent system for synchronously assimilating water level and flow of natural river | |
| CN115688246A (en) | Reservoir capacity simulation method and device under local coordinate system | |
| CN114297953B (en) | Mountain area storm disaster risk division method | |
| CN110457772B (en) | DEM flow direction estimation method combining plane curvature and steepest descent direction | |
| CN110263428B (en) | Riverbed evolution analysis method based on flow weighted average streamline length index | |
| CN114840896B (en) | Waterlogging and ponding simulation analysis method based on urban road BIM | |
| CN112015784B (en) | Wind condition data mining method and device, wind measuring device and data mining equipment | |
| CN111125893B (en) | Non-dispersive water flow path simulation method based on DEM and flow collection | |
| CN114492206A (en) | Calculation method for determining tunnel excavation footage based on fractured rock mass joint model |
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 |