CN115358650B - Flood disaster emergency risk avoiding and material real-time allocation method - Google Patents

Flood disaster emergency risk avoiding and material real-time allocation method Download PDF

Info

Publication number
CN115358650B
CN115358650B CN202211303256.0A CN202211303256A CN115358650B CN 115358650 B CN115358650 B CN 115358650B CN 202211303256 A CN202211303256 A CN 202211303256A CN 115358650 B CN115358650 B CN 115358650B
Authority
CN
China
Prior art keywords
data
flood
model
time
generating
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
CN202211303256.0A
Other languages
Chinese (zh)
Other versions
CN115358650A (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.)
Pearl River Hydraulic Research Institute of PRWRC
Original Assignee
Pearl River Hydraulic Research Institute of PRWRC
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 Pearl River Hydraulic Research Institute of PRWRC filed Critical Pearl River Hydraulic Research Institute of PRWRC
Priority to CN202211303256.0A priority Critical patent/CN115358650B/en
Publication of CN115358650A publication Critical patent/CN115358650A/en
Application granted granted Critical
Publication of CN115358650B publication Critical patent/CN115358650B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/18Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • 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/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0637Strategic management or analysis, e.g. setting a goal or target of an organisation; Planning actions based on goals; Analysis or evaluation of effectiveness of goals
    • 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
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling
    • 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
    • Y02A10/00TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Educational Administration (AREA)
  • Tourism & Hospitality (AREA)
  • Development Economics (AREA)
  • Marketing (AREA)
  • Geometry (AREA)
  • General Business, Economics & Management (AREA)
  • General Engineering & Computer Science (AREA)
  • Game Theory and Decision Science (AREA)
  • Quality & Reliability (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Operations Research (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Primary Health Care (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Alarm Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention relates to the technical field of intelligent decision, in particular to an emergency risk avoiding and material real-time allocation method for flood disasters. The method comprises the following steps: acquiring original sampling data; constructing a flood simulation model; generating dynamic simulation process data of flood according to the flood simulation model; generating a maximum flood submerging water depth grid map according to the flood dynamic simulation process data; constructing a rapid and intelligent flood inundation prediction model according to the maximum flood inundation water depth grid map, and generating simulation process data of flood inundation; generating disaster-affected data according to the flooding simulation process data; generating danger-avoiding transfer real-time data and scheduling decision data according to the disaster-suffered data; constructing a material allocation model according to risk avoidance transfer real-time analysis and scheduling decision analysis, and generating material allocation real-time data and scheduling decision data according to the material allocation model; the method takes flood simulation, flood real-time prediction, risk avoiding transfer analysis and material allocation analysis as core technical processes to support flood emergency flood control decisions.

Description

Flood disaster emergency risk avoiding and material real-time allocation method
Technical Field
The invention relates to the technical field of intelligent decision, in particular to an emergency risk avoiding and material real-time allocation method for flood disasters.
Background
Flood disasters are one of the most common natural disasters in China. The comprehensive improvement of the disaster prevention, reduction and relief capacity of flood disasters is a key construction task for water safety guarantee. Emergency danger avoiding and emergency material allocation of personnel in flood affected areas are important non-engineering measures for guaranteeing life and property safety of people. Most of the traditional risk avoiding transfer or material allocation is static, and a risk avoiding arrangement plan, a material arrangement plan and a material allocation plan are planned and compiled in advance according to the disaster influence envelope range.
In practical application, the static plan has limited guidance benefit and poor application effect due to lack of disaster dynamic development process and consideration of supply and demand factors. The method and the system combine a flood forecasting model based on hydrodynamic force, a disaster area risk avoiding transfer model and a material demand forecasting and allocating model, comprehensively realize real-time simulation and research and judgment of flood disaster situation, and realize real-time and dynamic analysis of risk avoiding transfer and material allocation under the forecast disaster situation. The technical achievement provides powerful technical support for various flood disaster emergency decision scenes and has wide practical significance.
Disclosure of Invention
The invention provides a flood disaster emergency risk avoiding and real-time material allocation method for solving the technical problems;
in one embodiment of the invention, the method comprises the following steps:
step S1: acquiring original sampling data, and constructing a basic model data set according to the original sampling data, wherein the basic model data set comprises a one-dimensional river surge model, a two-dimensional earth surface model and an underground pipe network drainage model;
step S2: constructing a flood simulation model according to the basic model data set;
and step S3: generating flood dynamic simulation process data according to the flood simulation model;
and step S4: generating a maximum flood submerging water depth grid map according to the flood dynamic simulation process data;
step S5: constructing a flood submerging rapid intelligent prediction model according to the maximum flood submerging water depth grid map, and generating flood submerging simulation process data according to the flood submerging rapid intelligent prediction model;
step S6: generating disaster-affected data according to the flooding simulation process data;
step S7: generating danger avoiding transfer real-time data and scheduling decision data according to the disaster-suffered data;
step S8: and constructing a material allocation model according to risk avoidance transfer real-time analysis and scheduling decision analysis, and generating material allocation real-time data and scheduling decision data according to the material allocation model to execute material allocation.
The method comprises the steps of acquiring original sampling data and analyzing and constructing a basic model data set according to the original sampling data; constructing flood dynamic simulation process data according to the basic model data set; generating flood dynamic simulation process data according to the flood simulation model; generating a maximum flood submerging water depth grid map according to the flood dynamic simulation process data; constructing a flood inundation rapid intelligent prediction model according to the maximum flood inundation water depth grid map, and generating flood inundation simulation process data according to the flood inundation rapid intelligent prediction model; generating disaster-stricken data according to the flooding simulation process data; generating real-time analysis and scheduling decisions of risk avoidance transfer according to the disaster-suffered data; the material allocation model is constructed according to real-time risk transfer analysis and scheduling decision analysis, and the real-time material allocation analysis and scheduling decision are generated according to the material allocation model, so that technical support is provided for flood disaster emergency risk avoiding scheduling decision and material allocation when flood disasters occur, and efficient and convenient disaster area rescue is provided.
In an embodiment of the present invention, step S1 specifically includes:
acquiring original sampling data, wherein the original sampling data comprises position information, basic geographic information data, hydraulic structure basic data and scheduling scheme data, hydrological and historical disaster data and urban rainwater pipe network data, the basic geographic information data comprises river water system data, river terrain data, surface terrain data and land utilization data, the land utilization data comprises house data, road data and railway data, the hydraulic structure basic data and scheduling scheme data comprises sluice data, pump station data, embankment data, culvert data and water conservancy junction data, the hydrological and historical disaster data comprises historically measured rainfall data, historical flood flooding water level data and historical flood flooding range data, and the urban rainwater pipe network data comprises pipe network data and water outlet data;
analyzing according to the position information, the surface topographic data and the land utilization data, and performing production convergence calculation by adopting a comprehensive unit line method so as to construct a hydrological model of the non-built area;
generating river network topological relation data, river surge section assigned data, gate dam hydraulic structure parameters and scheduling regulation data according to river water system data, river terrain data, hydraulic structure basic data and scheduling scheme data;
analyzing according to the topological relation data of the river network, the assigned value data of the river section, the parameters of the hydraulic structure of the gate dam and the scheduling regulation data, thereby constructing a one-dimensional river model;
performing irregular mesh generation on the surface terrain data to generate irregular mesh surface terrain data;
encrypting the irregular grid ground surface terrain data, road data, railway data, sluice data, pump station data, embankment data and water conservancy junction data to generate a middle two-dimensional ground surface model;
constructing a two-dimensional ground surface model according to the middle two-dimensional ground surface model, the position information, the house data, the road data, the dike data and the culvert data;
constructing an underground pipe network drainage model according to the pipe network data, the drainage port data and the pump station data;
according to the embodiment, a non-built-up area hydrological model is constructed according to basic geographic information data analysis, a one-dimensional river surge model is constructed according to basic geographic information data, hydraulic structure basic data and scheduling data analysis, a two-dimensional surface model is constructed according to surface topographic data, hydraulic structure basic data and scheduling scheme data analysis, and an underground pipe network drainage model is constructed according to urban rainwater pipe network data analysis, so that a basic model data set is constructed according to multi-source and multi-type original sampling data, and premise preparation is made for constructing a flood simulation model in the next step.
In an embodiment of the present invention, step S2 specifically includes:
constructing a middle flooding simulation model according to the hydrological model, the one-dimensional river channel model, the two-dimensional earth surface model and the underground pipe network drainage model of the non-built area;
and carrying out parameter calibration, verification and inspection on the flood simulation model according to hydrology and historical disaster data to generate the flood simulation model.
According to the embodiment, the intermediate flood simulation model is analyzed and constructed according to the basic data set, and the parameter calibration, verification and inspection are carried out on the intermediate flood simulation model, so that the flood simulation model is generated, and preconditions are provided for the next step of generating flood dynamic simulation process data.
In an embodiment of the present invention, step S3 specifically includes:
generating a flood encounter situation acquisition control;
acquiring flood encounter scenario data through the flood encounter scenario acquisition control, wherein the flood encounter scenario data comprises flood data, tide data, flood embankment data and scheduling data;
generating typical situation condition directory data of flood disasters according to the flood encounter situation data;
and generating flood dynamic simulation process data for the flood simulation model according to the typical situation condition catalog data of the flood disaster.
According to the method and the device, the flood encounter situation data is obtained through the flood encounter situation generation obtaining control, and the flood disaster typical situation condition catalogue data is generated through the flood encounter situation data analysis, so that the flood simulation model is analyzed to generate the flood dynamic simulation process data, and precondition preparation is made for generating the maximum flood submerging water depth raster map in the next step.
In an embodiment of the present invention, step S4 specifically includes:
generating space-time raster data according to the flood dynamic simulation process data;
generating a single grid point data set according to the space-time grid data;
and combining the single grid point data sets to generate a maximum flooding water depth grid map of the flooding.
According to the embodiment, the space-time raster data are generated according to the flood dynamic simulation process data, and the space-time raster data are analyzed to generate the single-raster-point data set, so that the maximum flooding water depth raster map is generated in a combined mode, and the premise preparation is made for the next generation of the flood flooding simulation process data.
In an embodiment of the present invention, step S5 specifically includes:
generating a rainstorm data acquisition control;
acquiring rainstorm data through a rainstorm data acquisition control, wherein the rainstorm data comprises typical working condition rainstorm data and atypical working condition rainstorm data;
generating a space-time rainfall grid map according to typical working condition rainstorm data;
constructing an intermediate intelligent prediction model according to the space-time rainfall grid map, the original sampling data and the maximum flooding water depth grid map of the flooding;
generating a maximum flooding water depth grid map of the atypical working condition flooding according to the atypical working condition rainstorm data;
generating a flood submerging fast intelligent prediction model for the intermediate intelligent prediction model according to the atypical working condition flood maximum submerging water depth grid map;
and generating flood inundation simulation process data according to the flood inundation rapid intelligent prediction model.
The embodiment acquires the rainstorm data by generating a rainstorm data acquisition control and acquiring the rainstorm data through the rainstorm data acquisition control; generating a space-time rainfall grid diagram according to typical working condition rainstorm data; constructing an intermediate intelligent prediction model according to the space-time rainfall grid map, the original sampling data and the maximum flooding water depth grid map of the flooding; generating a maximum flooding water depth grid map of the atypical working condition flooding according to the atypical working condition rainstorm data; generating a flood submerging fast intelligent prediction model for the intermediate intelligent prediction model according to the atypical working condition flood maximum submerging water depth grid map; and generating flood inundation simulation process data according to the flood inundation rapid intelligent prediction model, thereby preparing for generating disaster-affected data in the next step.
In an embodiment of the present invention, step S6 specifically includes:
acquiring standing population data, land utilization layer data, placement point position data and placement point capacity data, wherein the standing population data comprise village population data and community standing population, and the land utilization layer data comprise residential area data and business area data;
generating population space distribution map data according to the data of the standing population, the data of the residential area and the data of the business area;
generating disaster-stricken data according to the flooding simulation process data and the population space distribution map data, wherein the disaster-stricken data comprises disaster-stricken area data and disaster-stricken population data;
in the embodiment, population space distribution map data are generated according to the population data of the permanent residence, the land use map layer data, the position data of the placement points and the capacity data of the placement points; and disaster suffering data is generated according to the flood inundation simulation process data and the population space distribution map data, so that preparation is made for real-time analysis and scheduling decision of risk avoiding transfer generation in the next step.
In an embodiment of the present invention, step S7 specifically includes:
calculating by a risk avoidance transfer real-time analysis calculation formula according to the disaster suffered data to generate risk avoidance transfer real-time analysis and scheduling decision data, wherein the risk avoidance transfer real-time analysis and scheduling decision data comprise recommended route transfer route data, predicted transfer time consumption data and predicted emplacer data, the risk avoidance transfer real-time analysis calculation formula comprises constraint conditions, and the constraint conditions comprise a transfer process safety degree constraint condition, a emplacer point capacity constraint, a maximum distance and time consumption constraint;
the risk avoidance transfer real-time analysis and calculation formula is specifically as follows:
Figure DEST_PATH_IMAGE001
Figure DEST_PATH_IMAGE002
represents the minimum value after the summation of the transfer time of all the transfer objects,
Figure DEST_PATH_IMAGE003
represents from first to second
Figure DEST_PATH_IMAGE004
The number of individual villages/communities is,
Figure DEST_PATH_IMAGE005
represents from first to second
Figure DEST_PATH_IMAGE006
The grid cells of the individual communities are,
Figure DEST_PATH_IMAGE007
representing the time consumed by the transfer of a single transfer object;
the transfer process safety degree constraint conditions are specifically as follows:
Figure DEST_PATH_IMAGE008
Figure DEST_PATH_IMAGE009
is as follows
Figure DEST_PATH_IMAGE010
The transfer object passes through
Figure DEST_PATH_IMAGE011
The safety factor when the route is transferred;
Figure DEST_PATH_IMAGE012
the set minimum safety factor is obtained;
Figure DEST_PATH_IMAGE013
Figure DEST_PATH_IMAGE014
represents the first
Figure DEST_PATH_IMAGE015
The transfer route of the transfer object is from 1 st to
Figure DEST_PATH_IMAGE016
A set of individual sub-segments;
Figure DEST_PATH_IMAGE017
represents a flooded road segment;
Figure DEST_PATH_IMAGE018
forecasting the inundated road section in the representative inundated simulation result;
the placement point capacity constraints are specifically:
Figure DEST_PATH_IMAGE019
Figure DEST_PATH_IMAGE020
is a first
Figure DEST_PATH_IMAGE021
The number of persons who transfer the object;
Figure DEST_PATH_IMAGE022
is a first
Figure 870709DEST_PATH_IMAGE011
Each of the arrangement regions accommodates a total amount;
Figure DEST_PATH_IMAGE023
is an emergency disposal factor;
the maximum distance and time-consuming constraints are specifically:
Figure DEST_PATH_IMAGE024
Figure DEST_PATH_IMAGE025
wherein the content of the first and second substances,
Figure DEST_PATH_IMAGE026
for transferring an object
Figure DEST_PATH_IMAGE027
The total length of the bar path;
Figure DEST_PATH_IMAGE028
the length of a sub road section in a transfer route;
Figure DEST_PATH_IMAGE029
in order to allow for the maximum time consumption,
Figure DEST_PATH_IMAGE030
for transferring the object
Figure DEST_PATH_IMAGE031
The time taken for the bar route;
Figure DEST_PATH_IMAGE032
time is consumed for sub-road sections in the transfer route;
Figure DEST_PATH_IMAGE033
is a time consumption coefficient;
Figure DEST_PATH_IMAGE034
the maximum time consumption allowed.
The risk avoidance transfer real-time analysis and scheduling decision data are generated by calculating according to disaster-suffering data through a risk avoidance transfer real-time analysis calculation formula, wherein the risk avoidance transfer real-time analysis calculation formula is set to be the shortest overall transfer of all groups, and a plurality of constraint conditions are set, wherein the constraint conditions comprise a transfer process safety degree constraint condition, a placement point capacity constraint, a maximum distance constraint and a time-consuming constraint, so that the risk avoidance transfer requirements under various disaster conditions can be obtained more accurately and more stably; the embodiment fully considers the constraint problem in the real condition, thereby providing an efficient risk avoidance transfer real-time analysis and scheduling decision.
In an embodiment of the present invention, step S8 specifically includes:
calculating according to risk avoidance transfer real-time analysis and scheduling decision data through a settling point calculation formula, and constructing a material allocation model;
the arrangement point calculation formula is specifically as follows:
Figure DEST_PATH_IMAGE035
Figure DEST_PATH_IMAGE036
is a first
Figure 335319DEST_PATH_IMAGE027
The total amount of supply of materials required by each arrangement point;
Figure DEST_PATH_IMAGE037
supply required for a single person;
Figure DEST_PATH_IMAGE038
is as follows
Figure 318319DEST_PATH_IMAGE027
A predicted number of placements for each placement point;
Figure DEST_PATH_IMAGE039
is a first
Figure 933102DEST_PATH_IMAGE027
The quantity of materials at each placement point is already;
according to the above formula, the material satisfaction of each placement point can be expressed as:
Figure DEST_PATH_IMAGE040
Figure DEST_PATH_IMAGE041
is a first
Figure 950736DEST_PATH_IMAGE027
The material satisfaction degree of each arrangement point;
Figure DEST_PATH_IMAGE042
to supplement supply to
Figure 242041DEST_PATH_IMAGE027
The amount of material at each placement site;
Figure DEST_PATH_IMAGE043
is a first
Figure 243495DEST_PATH_IMAGE027
The quantity of materials is available at each arrangement point;
Figure DEST_PATH_IMAGE044
is a first
Figure 160504DEST_PATH_IMAGE027
The total amount of material supply required by each placement point;
calculating through a material allocation analysis formula according to the material allocation model to generate material allocation real-time data and scheduling decision data, wherein the material allocation analysis formula comprises material supply path constraint conditions and maximum path constraint conditions;
wherein the material allocation analysis formula is specifically as follows:
Figure DEST_PATH_IMAGE045
Figure DEST_PATH_IMAGE046
or S represents
Figure 298224DEST_PATH_IMAGE027
The supply satisfaction of the individual placement points,
Figure DEST_PATH_IMAGE047
representing the highest value of the summation of the supply satisfaction degrees of each placement point;
the material supply path constraint conditions are specifically as follows:
Figure DEST_PATH_IMAGE048
Figure DEST_PATH_IMAGE049
represents the first
Figure DEST_PATH_IMAGE050
The delivery route of the individual material supply team is from 1 st to
Figure DEST_PATH_IMAGE051
A set of individual sub-segments;
Figure DEST_PATH_IMAGE052
representing a flooded roadA segment;
Figure DEST_PATH_IMAGE053
representing forecasted flooding stretch in flooding simulation process data.
The maximum distance constraint conditions are specifically as follows:
Figure DEST_PATH_IMAGE054
Figure DEST_PATH_IMAGE055
for supply of goods and materials to the first
Figure 386528DEST_PATH_IMAGE031
The total length of the bar path;
Figure DEST_PATH_IMAGE056
is the sub-road length in the supply route;
Figure DEST_PATH_IMAGE057
to allow for maximum time consumption.
The method comprises the steps of calculating according to risk avoidance transfer real-time analysis and scheduling decision data through a settling point calculation formula, constructing a material allocation model, calculating according to the material allocation model through a material allocation analysis formula, and generating material allocation real-time data and scheduling decision data, wherein the settling point calculation formula considers the supply quantity required by an individual, the predicted settling number of settling points and the material quantity of the settling points so as to generate more accurate material supply quantity required by the settling points; the method comprises the steps that a material allocation analysis formula is calculated by setting the highest material supply satisfaction degree of arrangement points as a target, sufficient material support data are provided for people to transfer to emergency disaster avoidance in disaster, the material allocation analysis formula fully considers the supply satisfaction degree of the arrangement points and the highest value obtained by summing the supply satisfaction degree of each arrangement point under the constraints of material supply path constraint conditions and maximum distance constraint conditions, and technical support more fitting actual constraint conditions, material allocation in emergency risk avoidance and people requirements is provided.
The invention provides a flood disaster emergency refuge and material real-time allocation method by combining a multidisciplinary multi-type model, solves the problems of unreasonable refuge transfer and material allocation scheme and lack of dynamic consideration in a static pre-plan, realizes real-time analysis and scheduling of emergency refuge and material allocation in a flood disaster scene, and effectively supports flood prevention and disaster prevention emergency decision.
Drawings
Other features, objects and advantages of the present application will become more apparent upon reading of the following detailed description of non-limiting implementations with reference to the accompanying drawings in which:
FIG. 1a is a flowchart illustrating steps of an emergency risk avoiding and real-time material allocation method for flood disaster according to an embodiment of the present invention;
FIG. 1b is a flowchart illustrating steps of an emergency risk avoiding and real-time material allocation method for flood disaster according to an embodiment of the present invention;
fig. 2 is a schematic diagram of a flood inundation rapid intelligent prediction model of the flood disaster emergency risk avoiding and material real-time allocating method according to an embodiment of the present invention.
Detailed Description
The technical method of the invention is described in detail with reference to the accompanying drawings, and it is obvious that the described embodiments are a part of the embodiments of the invention, not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Furthermore, the drawings are merely schematic illustrations of the invention and are not necessarily drawn to scale. The same reference numerals in the drawings denote the same or similar parts, and thus their repetitive description will be omitted. Some of the block diagrams shown in the figures are functional entities and do not necessarily correspond to physically or logically separate entities. The functional entities may be implemented in the form of software, or in one or more hardware modules or integrated circuits, or in different networks and/or processor methods and/or microcontroller methods.
The invention provides an emergency risk avoiding and material real-time allocating method for flood disasters, which comprises the following steps of:
step S1: acquiring original sampling data, and constructing a basic model data set according to the original sampling data, wherein the basic model data set comprises a one-dimensional river surge model, a two-dimensional earth surface model and an underground pipe network drainage model;
step S2: constructing a flood simulation model according to the basic model data set;
and step S3: generating dynamic simulation process data of flood according to the flood simulation model;
and step S4: generating a maximum flood submerging water depth grid map according to the data of the flood dynamic simulation process;
step S5: constructing a flood submerging rapid intelligent prediction model according to the maximum flood submerging water depth grid map, and generating flood submerging simulation process data according to the flood submerging rapid intelligent prediction model;
step S6: generating disaster-affected data according to the flooding simulation process data;
step S7: generating danger-avoiding transfer real-time data and scheduling decision data according to the disaster-suffered data;
step S8: and constructing a material allocation model according to risk avoidance transfer real-time analysis and scheduling decision analysis, and generating material allocation real-time data and scheduling decision data according to the material allocation model to execute material allocation.
Specifically, for example, according to a specific engineering research area, multisource and multi-type data such as regional basic geographic information, hydrology, historical disasters and engineering structures are collected, a flood simulation model based on hydrology, hydrodynamic force and pipe network multi-model coupling is constructed, computational analysis of regional flood disaster-causing processes under various typical incoming water and rainstorm scenes is achieved, and a flood dynamic simulation process result under the typical scenes is output.
And rasterizing and superposing the results according to the dynamic simulation process results of the flood under the typical scene output by the flood simulation model to generate a maximum flood submerging water depth raster map. Based on data such as rainstorm, water coming, tide level and maximum flood submergence water depth grid map under typical working conditions, a ConvLSTM neural network is used for constructing a flood submergence rapid intelligent prediction model, calculating and outputting results.
And forecasting and analyzing the population suffering from the disaster according to the flooding simulation result output by the intelligent prediction model in real time. The shortest overall transfer time of the population in the disaster area is taken as a general target, multidimensional factors such as arrangement capacity, transfer routes and the like are comprehensively considered, a risk avoiding transfer model is constructed, and real-time analysis and scheduling decision of risk avoiding transfer under the flood and disaster situation are predicted.
On the basis of real-time risk avoidance transfer analysis results, a total goal of highest personnel material supply degree of a placement point is taken, multidimensional factors such as material distribution capacity and distribution speed are comprehensively considered, a material distribution model is built, and real-time material distribution analysis and scheduling decision under the situation of predicting risk avoidance transfer are achieved.
The method comprises the steps of acquiring original sampling data and analyzing and constructing a basic model data set according to the original sampling data; constructing flood dynamic simulation process data according to the basic model data set; generating dynamic simulation process data of flood according to the flood simulation model; generating a maximum flood submerging water depth grid map according to the data of the flood dynamic simulation process; constructing a flood inundation rapid intelligent prediction model according to the maximum flood inundation water depth grid map, and generating flood inundation simulation process data according to the flood inundation rapid intelligent prediction model; generating disaster-affected data according to the flooding simulation process data; generating real-time analysis and scheduling decisions of risk avoidance transfer according to the disaster-suffered data; and constructing a material allocation model according to real-time risk transfer analysis and scheduling decision analysis, and generating real-time material allocation analysis and scheduling decision according to the material allocation model, so that technical support is provided for flood disaster emergency risk avoiding scheduling decision and material allocation when flood disasters occur, and efficient and convenient disaster area rescue is provided.
In an embodiment of the present invention, step S1 specifically includes:
acquiring original sampling data, wherein the original sampling data comprises position information, basic geographic information data, hydraulic structure basic data and scheduling scheme data, hydrological and historical disaster data and urban rainwater pipe network data, the basic geographic information data comprises river water system data, river terrain data, surface terrain data and land utilization data, the land utilization data comprises house data, road data and railway data, the hydraulic structure basic data and scheduling scheme data comprises sluice data, pump station data, embankment data, culvert data and water conservancy junction data, the hydrological and historical disaster data comprises historically measured rainfall data, historical flood flooding water level data and historical flood flooding range data, and the urban rainwater pipe network data comprises pipe network data and water outlet data;
specifically, for example, the data required for analyzing flood in urban areas include basic geographic information, hydrology and flood, flood control and drainage (water) projects and structures, flood scheduling schemes and project scheduling rules, land utilization, historical flood and the like, and the data collection specifically includes the following steps: the basic geographic information data comprise vector information and elevation data of terrain and landform, river water system, river section, underwater terrain, administrative division, residential area distribution, traffic road and the like which are generated or updated recently in the calculation range. The scale bar of the basic base map is not less than 1:2000, the ratio of the river section to the river underwater topography map is not less than 1:2000, the scale of the deep sea chart is not less than 1:25000, the scale of the chart in the shallow sea is not less than 1:10000. the basic geographic information should meet the requirements of timeliness and practicability. 2) The hydrological and flood data comprise actual measurement data such as rainfall, water level, flow, tide and the like, design data such as design rainstorm, design flood, design tide level and the like, data such as water level-flow relation, water level-area-volume relation and the like of a hydrological control station reflecting storage and discharge characteristics of rivers, lakes, reservoirs and stagnant and flooded areas, related data of drainage (water) areas, spatial position information of related hydrological stations, water level stations, tide level stations and rainfall stations and the like. The hydrologic data should meet reliability, consistency and representativeness requirements. 3) The flood control and drainage engineering and structure data comprise reservoir, dike, gate dam, pump station, drainage pipe network, bridge and culvert, underground facilities and other engineering, and characteristic parameters and spatial position information of linear ground objects which are 0.5m above the ground. Engineering and construction data should meet the requirements of reality, timeliness and accuracy. 4) The flood scheduling scheme and the engineering scheduling rule comprise flood control plans of all levels, a flood scheduling scheme, a flood defense scheme, reservoir, flood storage and stagnation areas, pump stations, flood diversion channels, gate dam flood control scheduling application rules and the like. The flood scheduling scheme and the engineering scheduling rules should meet the requirements of timeliness and authority. 5) The land utilization data comprises land utilization, remote sensing images, crop types and distribution during flood, and the like. The scale of the soil utilization map is not less than 1:10000, the resolution of the remote sensing image is not lower than 2m. The land utilization data should meet the requirements of practicability and timeliness. 6) The historical flood data comprises historical flood (rainstorm, storm surge, dam break and the like) hydrological characteristics (actually measured water level or flood mark in a station flood process, river course and a submerged area, a submerged range, submerged duration, flood arrival time and the like), dam burst (overflow) condition, project and project scheduling when flood occurs and the like.
Constructing a hydrological model of the non-built area according to the position information, the surface topographic data and the land utilization data, wherein the constructing step is specifically to construct the model by adopting a comprehensive unit line method and perform production convergence calculation;
generating river network topological relation data, river surge section assigned data, gate dam hydraulic structure parameters and scheduling regulation data according to river water system data, river terrain data, hydraulic structure basic data and scheduling scheme data;
analyzing according to the topological relation data of the river network, the assigned value data of the river section, the parameters of the hydraulic structure of the gate dam and the scheduling regulation data, thereby constructing a one-dimensional river model;
performing irregular mesh generation on the surface terrain data to generate irregular mesh surface terrain data;
encrypting the irregular grid ground surface terrain data, road data, railway data, sluice data, pump station data, embankment data and water conservancy junction data to generate a middle two-dimensional ground surface model;
constructing a two-dimensional earth surface model according to the middle two-dimensional earth surface model, the position information, the house data, the road data, the embankment data and the culvert data;
and constructing an underground pipe network drainage model according to the pipe network data, the water outlet data and the pump station data.
Specifically, for example, based geographic information data such as river water systems, river terrain, surface terrain, land utilization, and the like of an area are collected according to a specific study area; collecting basic data and scheduling scheme data of hydraulic structures such as sluice, pump station, dike, hydro-junction and the like; collecting hydrological and historical disaster data such as historical measured rainfall, historical flood submerging level and range; collecting urban rainwater pipe network data; when the important river channel terrain and surface terrain data are lacked, river channel section terrain measurement and surface terrain aviation flight measurement tasks are carried out.
Constructing a hydrological model of a non-built area: hydrologic zoning is performed on non-built areas in the research area. The non-built areas are mainly mountainous areas, and a model is built by adopting a comprehensive unit line to calculate the production convergence.
Constructing a one-dimensional river model: and (4) utilizing the river and hydraulic data to complete the establishment of the topological relation of the river network, the assignment of the section of the river, the setting of the parameters and the dispatching rules of the hydraulic buildings of the gate dam and the examination of rationality.
Constructing a two-dimensional earth surface model: and (3) adopting irregular grids to subdivide the earth surface grids, and automatically encrypting the grids to adapt to the change of the terrain and the ground features according to the terrain in the calculation area, the distribution and the trend of a water system, a dike, a road, a railway and the like. Generalizing the water-blocking function of houses, roads and dikes in the area and the water passing function of culverts.
Constructing a drainage model of the underground pipe network: and checking the topological relation of the pipe network by utilizing the pipe network data, combining the parameters and the scheduling rules of hydraulic buildings such as pipe network sections, water outlet distribution, pump stations and the like, and performing reasonability examination.
Constructing a flood control and drainage scheduling model: the flood control and drainage scheduling model is based on multi-factor combined scheduling of reservoirs, gate dams, pump stations and the like.
According to the embodiment, a non-built-up area hydrological model is constructed according to basic geographic information data analysis, a one-dimensional river surge model is constructed according to basic geographic information data, hydraulic structure basic data and scheduling data analysis, a two-dimensional surface model is constructed according to surface topographic data, hydraulic structure basic data and scheduling scheme data analysis, and an underground pipe network drainage model is constructed according to urban rainwater pipe network data analysis, so that a basic model data set is constructed according to multi-source and multi-type original sampling data, and premise preparation is made for constructing a flood simulation model in the next step.
In an embodiment of the present invention, step S2 specifically includes:
constructing a middle flood simulation model according to the hydrological model, the one-dimensional river channel model, the two-dimensional surface model and the underground pipe network drainage model of the non-built area;
and carrying out parameter calibration, verification and inspection on the flood simulation model according to the hydrologic and historical disaster data to generate the flood simulation model.
Specifically, for example, multi-model coupling: the model coupling of the built-up area is completely coupled within a time step, so that the water exchange between the confluent or overflowing river and the ground surface, the drainage or reverse flow of the pipe network and the river are reflected in real time;
and water exchange occurs at the drainage port, and water exchange between the earth surface and the pipe network inspection well such as confluence or full pipe back irrigation occurs. The production convergence calculation model of the drainage basin of the non-built area can be used as a boundary condition to be coupled with the built area model, and then the flood simulation model of the whole drainage basin/area is finally built.
Model calibration, verification and inspection: parameter calibration, verification and inspection are carried out on the watershed flood simulation model according to historical flood survey data, related design reports (water surface line calculation results designed in the initial design stage) of each river channel and mutual evidence of a mathematical model, wherein the adopted actual measurement data mainly are historical water immersion point distribution data of the watershed (region).
According to the embodiment, the intermediate flood simulation model is constructed according to the basic data set analysis, and the parameter calibration, verification and inspection are carried out on the intermediate flood simulation model, so that the flood simulation model is generated, and the precondition is prepared for generating the flood dynamic simulation process data in the next step.
In an embodiment of the present invention, step S3 specifically includes:
generating a flood encounter situation acquisition control;
acquiring flood encounter scenario data through the flood encounter scenario acquisition control, wherein the flood encounter scenario data comprises flood data, tide data, crash and flood dike data and scheduling data;
generating typical situation condition catalog data of flood disasters according to the flood encounter situation data;
and generating flood dynamic simulation process data for the flood simulation model according to the typical situation condition catalog data of the flood disaster.
Specifically, for example, based on the calculation analysis of regional flood disaster-causing processes under various typical incoming water and rainstorm scenes, the dynamic simulation process results of flood under typical scenes are output. The following table shows a simulation result table at a certain time outputted by the flooding model.
Figure DEST_PATH_IMAGE058
According to the embodiment, the flood encounter situation data is obtained through the flood encounter situation obtaining control, and the flood disaster typical situation condition catalogue data is generated through the flood encounter situation data analysis, so that the flood dynamic simulation process data is analyzed and generated through the flood simulation model, and precondition preparation is made for generating the maximum flood submerging water depth raster map in the next step.
In an embodiment of the present invention, step S4 specifically includes:
generating space-time grid data according to the flood dynamic simulation process data;
generating a single grid point data set according to the space-time grid data;
and combining the single grid point data sets to generate a maximum flooding water depth grid map of the flooding.
Specifically, for example, the output dynamic simulation process results of flooding under each typical scenario are marked as vector grid data and each grid water depth data; rendering is carried out according to the grid water depth data and the submerging water depth, and the grid water depth data is converted into space-time grid data, so that the space-time grid data is sent to a ConvLSTM neural network as input to be trained to obtain space-time grid data under typical situations, wherein the rasterization operation of the achievement is realized by GIS operation software such as ArcGIS, supermap and the like, and can also be realized by researching and developing a special format conversion program by combining with space data processing open source libraries such as GDAL, dotSpatial and the like;
performing superposition analysis according to space-time raster data under typical situations, and circularly reading raster points in a raster image layer at each moment; combining the maximum pixel value (water depth value) of each grid point at each moment to generate a flooding submerging water depth grid map;
specifically, for example, the flooding model is output to the result to be subjected to rasterization processing, a flooding water depth map is generated, integration application is performed in the information system, and after the result rasterization processing is integrated, a user can visually see the dynamic flooding condition predicted at each moment in the system.
According to the embodiment, the space-time raster data are generated according to the flood dynamic simulation process data, and the space-time raster data are analyzed to generate the single-raster-point data set, so that the maximum flooding water depth raster map is generated in a combined mode, and the premise preparation is made for the next generation of the flood flooding simulation process data.
In an embodiment of the present invention, step S5 specifically includes:
generating a rainstorm data acquisition control;
acquiring rainstorm data through a rainstorm data acquisition control, wherein the rainstorm data comprises typical working condition rainstorm data and atypical working condition rainstorm data;
generating a space-time rainfall grid diagram according to typical working condition rainstorm data;
constructing an intermediate intelligent prediction model according to the space-time rainfall grid map, the original sampling data and the maximum flooding water depth grid map of the flooding;
generating an atypical working condition flood maximum submerging water depth grid map according to atypical working condition rainstorm data;
generating a flood submerging fast intelligent prediction model for the intermediate intelligent prediction model according to the atypical working condition flood maximum submerging water depth grid map;
and generating flood inundation simulation process data according to the flood inundation rapid intelligent prediction model.
Specifically, for example, a rainstorm grid map and a flood maximum flooding water depth grid map under a plurality of rainstorm scenes are used as output, all typical working condition data are used as a training set, a ConvLSTM neural network is used for training, and an intelligent prediction model is constructed. The schematic diagram of the model is shown in fig. 2. The ConvLSTM neural network is suitable for space-time sequence grid data, and mainly comprises a multi-scale feature extraction sub-network and a multi-scale flooding identification sub-network in the embodiment. Through the composition of two self-networks and black box training, the rapid prediction and generation of the maximum flooding grid map under the rainstorm situation defined by a user are realized.
The embodiment acquires the rainstorm data by generating the rainstorm data acquisition control and acquiring the rainstorm data through the rainstorm data acquisition control; generating a space-time rainfall grid diagram according to typical working condition rainstorm data; constructing an intermediate intelligent prediction model according to the space-time rainfall grid map, the original sampling data and the maximum flooding water depth grid map of the flooding; generating a maximum flooding water depth grid map of the atypical working condition flooding according to the atypical working condition rainstorm data; generating a flood submerging fast intelligent prediction model for the intermediate intelligent prediction model according to the atypical working condition flood maximum submerging water depth grid map; and generating flood inundation simulation process data according to the flood inundation rapid intelligent prediction model, thereby preparing for generating disaster-affected data in the next step.
In an embodiment of the present invention, step S6 specifically includes:
acquiring standing population data, land utilization layer data, positioning point position data and positioning point capacity data, wherein the standing population data comprises village population data and community standing population, and the land utilization layer data comprises residential area data and business area data;
generating population space distribution map data according to the data of the standing population, the data of the residential area and the data of the business area;
generating disaster-stricken data according to the flooding simulation process data and the population space distribution map data, wherein the disaster-stricken data comprises disaster-stricken area data and disaster-stricken population data;
specifically, for example, with the goal of minimizing the overall transfer time consumption, risk avoidance transfer analysis is performed on disaster-affected points and placement points (risk avoidance points) in the area. The intra-area element information is shown in the following table.
Figure DEST_PATH_IMAGE059
According to the embodiment, population space distribution map data are generated according to the population data of the regular lives, the land use map layer data, the position data of the arrangement points and the capacity data of the arrangement points; and disaster-suffering data is generated according to the flooding simulation process data and the population space distribution map data, so that preparation is made for real-time analysis and scheduling decision of risk avoidance transfer generation in the next step.
In an embodiment of the present invention, step S7 specifically includes:
calculating by a risk avoidance transfer real-time analysis calculation formula according to the disaster suffered data to generate risk avoidance transfer real-time analysis and scheduling decision data, wherein the risk avoidance transfer real-time analysis and scheduling decision data comprise recommended route transfer route data, predicted transfer time consumption data and predicted emplacer data, the risk avoidance transfer real-time analysis calculation formula comprises constraint conditions, and the constraint conditions comprise a transfer process safety degree constraint condition, a emplacer point capacity constraint, a maximum distance and time consumption constraint;
the risk avoiding transfer real-time analysis and calculation formula specifically comprises:
Figure DEST_PATH_IMAGE060
for example, it can be simplified as:
Figure DEST_PATH_IMAGE061
Figure DEST_PATH_IMAGE062
represents the minimum value after the summation of the transfer time of all the transfer objects,
Figure 531071DEST_PATH_IMAGE021
represents from first to second
Figure DEST_PATH_IMAGE063
The number of individual villages/communities is,
Figure DEST_PATH_IMAGE064
represents from first to second
Figure DEST_PATH_IMAGE065
The number of the grid cells in the community,
Figure DEST_PATH_IMAGE066
or
Figure DEST_PATH_IMAGE067
Representing the time consumed by the transfer of a single transfer object;
the transfer process safety degree constraint conditions are specifically as follows:
Figure DEST_PATH_IMAGE068
Figure DEST_PATH_IMAGE069
is as follows
Figure 159760DEST_PATH_IMAGE027
The transfer object passes through
Figure 683145DEST_PATH_IMAGE011
Safety factors when the lines are transferred;
Figure DEST_PATH_IMAGE070
the set minimum safety factor is obtained;
Figure DEST_PATH_IMAGE071
for example, it can be simplified as:
Figure DEST_PATH_IMAGE072
Figure DEST_PATH_IMAGE073
or
Figure DEST_PATH_IMAGE074
Represents the first
Figure 768782DEST_PATH_IMAGE031
The transfer route of the transfer object is from 1 st to
Figure DEST_PATH_IMAGE075
A set of sub-segments;
Figure DEST_PATH_IMAGE076
or
Figure DEST_PATH_IMAGE077
Represents a flooded road segment;
Figure DEST_PATH_IMAGE078
or
Figure DEST_PATH_IMAGE079
Forecasting an inundated road section in the representative inundated simulation result;
the placement point capacity constraints are specifically:
Figure DEST_PATH_IMAGE080
for example, it can be simplified as:
Figure DEST_PATH_IMAGE081
Figure DEST_PATH_IMAGE082
or
Figure DEST_PATH_IMAGE083
Is a first
Figure 636506DEST_PATH_IMAGE021
The number of persons who transfer the object;
Figure DEST_PATH_IMAGE084
or
Figure DEST_PATH_IMAGE085
Is a first
Figure 380471DEST_PATH_IMAGE011
Each of the arrangement regions accommodates a total amount;
Figure 492784DEST_PATH_IMAGE033
is an emergency disposal factor;
the maximum distance and time consumption constraints are specifically as follows:
Figure DEST_PATH_IMAGE086
for example, it can be simplified as:
Figure DEST_PATH_IMAGE087
Figure DEST_PATH_IMAGE088
wherein, the first and the second end of the pipe are connected with each other,
Figure DEST_PATH_IMAGE089
or
Figure DEST_PATH_IMAGE090
For transferring the object
Figure 14901DEST_PATH_IMAGE027
The total length of the bar path;
Figure DEST_PATH_IMAGE091
for sub-track length in transfer routes
Figure DEST_PATH_IMAGE092
Or
Figure DEST_PATH_IMAGE093
In order to allow for the maximum time-consuming,
Figure DEST_PATH_IMAGE094
or
Figure DEST_PATH_IMAGE095
For transferring the object
Figure 415926DEST_PATH_IMAGE027
The time consumed by the bar route;
Figure DEST_PATH_IMAGE096
or
Figure DEST_PATH_IMAGE097
Time is consumed for sub-road sections in the transfer route;
Figure DEST_PATH_IMAGE098
is a time consumption coefficient;
Figure DEST_PATH_IMAGE099
or
Figure DEST_PATH_IMAGE100
The maximum time consumption allowed.
Specifically, for example, the risk avoidance transfer real-time analysis calculation formula is solved and calculated to generate an output result, wherein the solving and calculating method can select a simulated annealing algorithm, a genetic algorithm, a tabu search algorithm, a particle swarm algorithm or other intelligent optimization algorithms, and the model output result includes the recommended risk avoidance transfer route of each risk area/point, the predicted transfer time consumption of each transfer unit and the predicted number of the persons to be placed at each placement point.
According to the method, calculation is carried out through a risk avoidance transfer real-time analysis calculation formula according to disaster-stricken data, and risk avoidance transfer real-time analysis and scheduling decision data are generated, wherein the risk avoidance transfer real-time analysis calculation formula is set to be the shortest overall transfer of all groups, and a plurality of constraint conditions are set, wherein the constraint conditions comprise a transfer process safety degree constraint condition, a placement point capacity constraint, a maximum distance constraint and a time-consuming constraint, so that more accuracy and stability are obtained, and the risk avoidance transfer requirements under various disaster conditions can be adapted; the embodiment fully considers the constraint problem in the actual condition, thereby providing an efficient risk avoidance transfer real-time analysis and scheduling decision.
In an embodiment of the present invention, step S8 specifically includes:
calculating according to the risk avoidance transfer real-time analysis and scheduling decision data through a mounting point calculation formula, and constructing a material allocation model;
the setting point calculation formula is specifically as follows:
Figure DEST_PATH_IMAGE101
for example, it can be simplified as:
Figure DEST_PATH_IMAGE102
Figure DEST_PATH_IMAGE103
is as follows
Figure 855260DEST_PATH_IMAGE031
The total amount of material supply required by each placement point;
Figure DEST_PATH_IMAGE104
the supply required for a single person;
Figure DEST_PATH_IMAGE105
or
Figure DEST_PATH_IMAGE106
Is as follows
Figure 805767DEST_PATH_IMAGE027
A predicted number of placements for individual placement points;
Figure DEST_PATH_IMAGE107
is a first
Figure 311835DEST_PATH_IMAGE031
The quantity of materials is available at each arrangement point;
according to the above formula, the material satisfaction of each placement point can be expressed as:
Figure DEST_PATH_IMAGE108
Figure DEST_PATH_IMAGE109
Figure DEST_PATH_IMAGE110
or
Figure DEST_PATH_IMAGE111
Is as follows
Figure 200157DEST_PATH_IMAGE027
The material satisfaction degree of each arrangement point;
Figure DEST_PATH_IMAGE112
or
Figure DEST_PATH_IMAGE113
To supplement supply to
Figure 48376DEST_PATH_IMAGE027
The amount of material at each placement site;
Figure DEST_PATH_IMAGE114
is as follows
Figure 135280DEST_PATH_IMAGE027
The quantity of materials is available at each arrangement point;
Figure DEST_PATH_IMAGE115
is as follows
Figure 546670DEST_PATH_IMAGE027
The total amount of material supply required by each placement point;
calculating through a material allocation analysis formula according to the material allocation model to generate material allocation real-time data and scheduling decision data, wherein the material allocation analysis formula comprises material supply path constraint conditions and maximum path constraint conditions;
wherein the material allocation analysis formula is specifically as follows:
Figure DEST_PATH_IMAGE116
Figure DEST_PATH_IMAGE117
Figure DEST_PATH_IMAGE118
or
Figure DEST_PATH_IMAGE119
Represents the first
Figure 171555DEST_PATH_IMAGE027
The supply satisfaction of the individual placement points,
Figure DEST_PATH_IMAGE120
or
Figure DEST_PATH_IMAGE121
Representing the highest value of the summation of the supply satisfaction degrees of each placement point;
the material supply path constraint conditions are specifically as follows:
Figure DEST_PATH_IMAGE122
Figure DEST_PATH_IMAGE123
Figure DEST_PATH_IMAGE124
or
Figure DEST_PATH_IMAGE125
Represents the first
Figure 483850DEST_PATH_IMAGE027
The delivery route of the individual material supply team is from 1 st to
Figure 956420DEST_PATH_IMAGE075
A set of sub-segments;
Figure DEST_PATH_IMAGE126
or
Figure DEST_PATH_IMAGE127
Represents a flooded road segment;
Figure DEST_PATH_IMAGE128
or
Figure DEST_PATH_IMAGE129
Representing forecasted flooding stretch in flooding simulation process data.
The maximum distance constraint condition is specifically as follows:
Figure DEST_PATH_IMAGE130
for example, it can be simplified as:
Figure DEST_PATH_IMAGE131
Figure DEST_PATH_IMAGE132
or
Figure DEST_PATH_IMAGE133
For material supply vehicle
Figure 194503DEST_PATH_IMAGE027
The total length of the bar path;
Figure DEST_PATH_IMAGE134
is the sub-road length in the supply route;
Figure DEST_PATH_IMAGE135
or
Figure DEST_PATH_IMAGE136
To allow for maximum time consumption.
Specifically, for example, the above calculation formula is solved and calculated, so as to generate a model output result, wherein the solution calculation method may select a simulated annealing algorithm, a genetic algorithm, a tabu search algorithm, a particle swarm algorithm or other intelligent optimization algorithms, the model output result is a result scheme set which is sorted according to the total material satisfaction degree from high to low, and the result scheme includes operation target placement points, transportation routes, material amounts of placement points and material satisfaction degrees of supply points and supply teams;
and analyzing and adjusting according to the result scheme set and preset experience conditions to generate real-time material allocation analysis and scheduling decisions under the situations of forecasting flood submergence and risk avoidance transfer.
Specifically, for example, the material allocation analysis is performed on the placement points in the region with the material supply degree as a target. The supply information of the supplies in the area is shown in the following table:
Figure DEST_PATH_IMAGE137
the method comprises the steps of calculating according to risk avoidance transfer real-time analysis and scheduling decision data through a settling point calculation formula, constructing a material allocation model, calculating according to the material allocation model through a material allocation analysis formula, and generating material allocation real-time data and scheduling decision data, wherein the settling point calculation formula considers the supply quantity required by an individual, the predicted settling number of settling points and the material quantity of the settling points so as to generate more accurate material supply quantity required by the settling points; the method comprises the steps that a material allocation analysis formula is calculated by setting the highest material supply satisfaction degree of arrangement points as a target, sufficient material support data are provided for people to transfer to emergency disaster avoidance in disaster, the material allocation analysis formula fully considers the supply satisfaction degree of the arrangement points and the highest value obtained by summing the supply satisfaction degree of each arrangement point under the constraints of material supply path constraint conditions and maximum distance constraint conditions, and technical support more fitting actual constraint conditions, material allocation in emergency risk avoidance and people requirements is provided.
The invention provides a flood disaster emergency refuge and material real-time allocation method by combining a multidisciplinary multi-type model, solves the problems of unreasonable refuge transfer and material allocation scheme and lack of dynamic consideration in a static pre-plan, realizes real-time analysis and scheduling of emergency refuge and material allocation in a flood disaster scene, and effectively supports flood prevention and disaster prevention emergency decision.
The present embodiments are therefore to be considered in all respects as illustrative and not restrictive, the scope of the invention being indicated by the appended claims rather than by the foregoing description, and all changes which come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein. Any reference signs in the claims shall not be construed as limiting the claim concerned.
The foregoing are merely exemplary embodiments of the present invention, which enable those skilled in the art to understand or practice the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the invention. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.

Claims (5)

1. A flood disaster emergency risk avoiding and material real-time allocation method is characterized by comprising the following steps:
step S1: acquiring original sampling data, and constructing a basic model data set according to the original sampling data, wherein the basic model data set comprises a one-dimensional river surge model, a two-dimensional earth surface model and an underground pipe network drainage model;
step S2: constructing a flood simulation model according to the basic model data set;
and step S3: generating flood dynamic simulation process data according to the flood simulation model;
and step S4: generating maximum flood submerging water depth raster map data according to the flood dynamic simulation process data;
step S5: constructing a flood inundation rapid intelligent prediction model according to the maximum flood inundation water depth raster map data, and generating flood inundation simulation process data according to the flood inundation rapid intelligent prediction model;
step S6: generating disaster-affected data according to the flooding simulation process data;
step S7: generating danger avoiding transfer real-time data and scheduling decision data according to the disaster-suffered data;
step S8: constructing a material allocation model according to the risk avoidance transfer real-time data and the scheduling decision data, and generating material allocation real-time data and scheduling decision data according to the material allocation model to execute material allocation;
wherein the step S3 specifically comprises the following steps:
generating a flood encounter situation acquisition control;
acquiring flood encounter scenario data through the flood encounter scenario acquisition control, wherein the flood encounter scenario data comprises flood data, tide data, flood embankment data and scheduling data;
generating typical situation condition catalog data of flood disasters according to the flood encounter situation data;
correcting the flood simulation model according to the typical situation condition catalog data of the flood disaster to generate dynamic simulation process data of the flood;
step S4 specifically includes:
generating space-time grid data according to the flood dynamic simulation process data;
generating a single grid point data set according to the space-time grid data;
combining the single grid point data sets to generate the maximum flooding submerging water depth grid map data;
the step S5 specifically comprises the following steps:
generating a rainstorm data acquisition control;
acquiring rainstorm data through a rainstorm data acquisition control, wherein the rainstorm data comprises typical working condition rainstorm data and atypical working condition rainstorm data;
generating a space-time rainfall grid diagram according to typical working condition rainstorm data;
constructing an intermediate intelligent prediction model according to the space-time rainfall raster image, the original sampling data and the maximum flood submerging water depth raster image data;
generating a maximum flooding water depth grid map of the atypical working condition flooding according to the atypical working condition rainstorm data;
correcting the intermediate intelligent prediction model according to the atypical working condition maximum flooding water depth grid map to generate a flooding fast intelligent prediction model;
generating flood inundation simulation process data according to the flood inundation rapid intelligent prediction model;
step S8 specifically includes:
calculating according to the risk avoidance transfer real-time data and the scheduling decision data through a mounting point calculation formula, and constructing a material allocation model;
the setting point calculation formula is specifically as follows:
Figure 767760DEST_PATH_IMAGE001
Figure 435502DEST_PATH_IMAGE002
is a first
Figure 556910DEST_PATH_IMAGE003
The total amount of supply of materials required by each arrangement point;
Figure 550274DEST_PATH_IMAGE004
the supply required for a single person;
Figure 773445DEST_PATH_IMAGE005
is as follows
Figure 244878DEST_PATH_IMAGE006
A predicted number of placements for individual placement points;
Figure 237104DEST_PATH_IMAGE007
is as follows
Figure 401369DEST_PATH_IMAGE008
The quantity of materials at each placement point is already;
calculating the material satisfaction degree of each arrangement point and expressing the material satisfaction degree as follows:
Figure 594060DEST_PATH_IMAGE009
Figure 869184DEST_PATH_IMAGE010
is a first
Figure 981496DEST_PATH_IMAGE006
The material satisfaction degree of each arrangement point;
Figure 51083DEST_PATH_IMAGE011
to supplement supply to
Figure 514426DEST_PATH_IMAGE006
The amount of material at each placement site;
Figure 327661DEST_PATH_IMAGE007
is as follows
Figure 543747DEST_PATH_IMAGE006
The quantity of materials at each placement point is already;
Figure 784236DEST_PATH_IMAGE002
is as follows
Figure 469295DEST_PATH_IMAGE006
The total amount of material supply required by each placement point;
calculating through a material allocation analysis formula according to the material allocation model to generate material allocation real-time data and scheduling decision data, wherein the material allocation analysis formula comprises a material supply path constraint condition and a maximum path constraint condition;
wherein the material allocation analysis formula is specifically as follows:
Figure 86221DEST_PATH_IMAGE012
Figure 907547DEST_PATH_IMAGE013
represents the first
Figure 69669DEST_PATH_IMAGE006
The material satisfaction degree of each arrangement point is determined,
Figure 242024DEST_PATH_IMAGE014
representing the highest value after the sum of the material satisfaction degrees of each arrangement point;
the material supply path constraint conditions are specifically as follows:
Figure 397062DEST_PATH_IMAGE015
Figure 604052DEST_PATH_IMAGE016
represents the first
Figure 920764DEST_PATH_IMAGE006
The delivery route of the individual material supply team is from 1 st to
Figure 564104DEST_PATH_IMAGE017
A set of sub-segments;
Figure 522833DEST_PATH_IMAGE018
represents a flooded road segment;
Figure 584330DEST_PATH_IMAGE019
representing forecasted flooding sections in the flooding simulation process data;
the maximum distance constraint condition is specifically as follows:
Figure 337522DEST_PATH_IMAGE020
Figure 484470DEST_PATH_IMAGE021
for supply of goods and materials to the first
Figure 981310DEST_PATH_IMAGE022
The total length of the bar path;
Figure 373678DEST_PATH_IMAGE023
is the sub-road length in the supply route;
Figure 563351DEST_PATH_IMAGE024
to allow for maximum time consumption.
2. The method according to claim 1, wherein step S1 is specifically:
acquiring original sampling data, wherein the original sampling data comprises position information, basic geographic information data, hydraulic structure basic data and scheduling scheme data, hydrologic and historical disaster data and urban rainwater pipe network data, the basic geographic information data comprises river water system data, river channel topographic data, surface topographic data and land utilization data, the land utilization data comprises house data, road data and railway data, the hydraulic structure basic data and scheduling scheme data comprises sluice data, pump station data, dike dam data, culvert data and water conservancy data, the hydrologic and historical disaster data comprises historically measured rainfall data, historical flood flooding water level data and historical flood flooding range data, and the urban rainwater pipe network data comprises pipe network data and water outlet data;
analyzing according to the position information, the surface topographic data and the land utilization data, and performing production convergence calculation by adopting a comprehensive unit line method so as to construct a hydrological model of the non-built area;
generating river network topological relation data, river surge section assigned data, gate dam hydraulic structure parameters and scheduling regulation data according to river water system data, river terrain data, hydraulic structure basic data and scheduling scheme data;
analyzing according to the topological relation data of the river network, the assigned value data of the river section, the parameters of the hydraulic structure of the gate dam and the scheduling regulation data, thereby constructing a one-dimensional river model;
performing irregular mesh generation on the surface terrain data to generate irregular mesh surface terrain data;
encrypting the irregular grid ground surface topographic data, road data, railway data, sluice data, pump station data, embankment data and hydro junction data into a grid to generate a middle two-dimensional ground surface model;
constructing a two-dimensional earth surface model according to the middle two-dimensional earth surface model, the position information, the house data, the road data, the embankment data and the culvert data;
and constructing an underground pipe network drainage model according to the pipe network data, the water outlet data and the pump station data.
3. The method according to claim 2, wherein step S2 is specifically:
constructing a middle flooding simulation model according to the hydrological model, the one-dimensional river surge model, the two-dimensional earth surface model and the underground pipe network drainage model of the non-built area;
and carrying out parameter calibration, verification and inspection on the flood simulation model according to the hydrologic and historical disaster data to generate the flood simulation model.
4. The method according to claim 2, wherein step S6 is specifically:
acquiring standing population data, land utilization layer data, positioning point position data and positioning point capacity data, wherein the standing population data comprises village population data and community standing population data, and the land utilization layer data comprises residential area data and business area data;
generating population space distribution map data according to the data of the standing population, the data of the residential area and the data of the business area;
and generating disaster-stricken data according to the flood inundation simulation process data and the population space distribution map data, wherein the disaster-stricken data comprises disaster-stricken area data and disaster-stricken population data.
5. The method according to claim 1, wherein step S7 is specifically:
calculating by a risk avoidance transfer real-time analysis calculation formula according to the disaster data to generate risk avoidance transfer real-time data and scheduling decision data, wherein the risk avoidance transfer real-time data and the scheduling decision data comprise recommended route transfer route data, predicted transfer time consumption data and predicted emplacer data, the risk avoidance transfer real-time analysis calculation formula comprises constraint conditions, and the constraint conditions comprise a transfer process safety degree constraint condition, a emplacer point capacity constraint and a maximum distance and time consumption constraint;
the risk avoidance transfer real-time analysis and calculation formula is specifically as follows:
Figure 932015DEST_PATH_IMAGE025
Figure 232546DEST_PATH_IMAGE026
represents the minimum value after the summation of the transfer time of all the transfer objects,
Figure 737477DEST_PATH_IMAGE027
represents from first to second
Figure 98051DEST_PATH_IMAGE028
The number of individual villages/communities is,
Figure 468859DEST_PATH_IMAGE029
represents from first to second
Figure 41922DEST_PATH_IMAGE030
The grid cells of the individual communities are,
Figure 932518DEST_PATH_IMAGE031
representing a single transferred object from first to second
Figure 198414DEST_PATH_IMAGE032
The transfer of an individual village/community is time consuming,
Figure 276092DEST_PATH_IMAGE033
representing a single transferred object from first to second
Figure 934737DEST_PATH_IMAGE034
From first to second in individual community grid cells
Figure 679839DEST_PATH_IMAGE032
The transfer of an individual village/community is time consuming;
the transfer process safety degree constraint conditions are specifically as follows:
Figure 116637DEST_PATH_IMAGE035
Figure 947190DEST_PATH_IMAGE036
is a first
Figure 862056DEST_PATH_IMAGE037
The object to be transferred passes
Figure 727244DEST_PATH_IMAGE029
Safety factors when the lines are transferred;
Figure 584210DEST_PATH_IMAGE038
is the set lowest safety factor;
Figure 902059DEST_PATH_IMAGE039
Figure 886196DEST_PATH_IMAGE040
represents the first
Figure 605890DEST_PATH_IMAGE008
The transfer route of the transfer object is from 1 st to
Figure 384490DEST_PATH_IMAGE017
A set of individual sub-segments;
Figure 189635DEST_PATH_IMAGE041
represents a flooded road segment;
Figure 459686DEST_PATH_IMAGE042
forecasting an inundated road section in the simulation process of flood inundation;
the capacity constraint of the placement points is specifically as follows:
Figure 33887DEST_PATH_IMAGE043
Figure 983388DEST_PATH_IMAGE044
is as follows
Figure 10250DEST_PATH_IMAGE027
The number of persons who transfer the object;
Figure 336189DEST_PATH_IMAGE045
is as follows
Figure 14164DEST_PATH_IMAGE029
Each of the arrangement regions accommodates a total amount;
Figure 400146DEST_PATH_IMAGE046
is an emergency disposal factor;
the maximum distance and time-consuming constraints are specifically:
Figure 914304DEST_PATH_IMAGE047
Figure 43934DEST_PATH_IMAGE048
wherein the content of the first and second substances,
Figure 61569DEST_PATH_IMAGE049
for transferring the object
Figure 884031DEST_PATH_IMAGE050
The total length of the bar path;
Figure 370638DEST_PATH_IMAGE051
the length of a sub road section in a transfer route;
Figure 38380DEST_PATH_IMAGE052
in order to allow for the maximum time-consuming,
Figure 441680DEST_PATH_IMAGE053
the time consumed for transferring the ith route of the object;
Figure 169464DEST_PATH_IMAGE054
time is consumed for sub-road sections in the transfer route;
Figure 392635DEST_PATH_IMAGE055
is a time consumption coefficient;
Figure 113335DEST_PATH_IMAGE056
the maximum time consumption allowed.
CN202211303256.0A 2022-10-24 2022-10-24 Flood disaster emergency risk avoiding and material real-time allocation method Active CN115358650B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211303256.0A CN115358650B (en) 2022-10-24 2022-10-24 Flood disaster emergency risk avoiding and material real-time allocation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211303256.0A CN115358650B (en) 2022-10-24 2022-10-24 Flood disaster emergency risk avoiding and material real-time allocation method

Publications (2)

Publication Number Publication Date
CN115358650A CN115358650A (en) 2022-11-18
CN115358650B true CN115358650B (en) 2023-01-10

Family

ID=84008707

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211303256.0A Active CN115358650B (en) 2022-10-24 2022-10-24 Flood disaster emergency risk avoiding and material real-time allocation method

Country Status (1)

Country Link
CN (1) CN115358650B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116090677B (en) * 2023-04-10 2023-06-20 湖南大学 Air-ground emergency resource planning method considering electric power-communication-traffic network coupling
CN116401473B (en) * 2023-06-07 2023-08-15 上海维智卓新信息科技有限公司 Evacuation location information recommendation method, electronic device and storage medium
CN116957303B (en) * 2023-09-20 2023-11-28 河海大学 Emergency response scheduling decision method and system for flood disaster scene
CN118153787A (en) * 2024-05-11 2024-06-07 珠江水利委员会珠江水利科学研究院 Rain and tide disaster emergency risk avoiding path and material allocation optimization method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104851360A (en) * 2014-02-14 2015-08-19 杭州贵仁科技有限公司 Method and system for generating flood risk map
CN111652777A (en) * 2020-05-15 2020-09-11 长江勘测规划设计研究有限责任公司 Flood emergency risk avoiding method
CN114240119A (en) * 2021-12-10 2022-03-25 浙江水利水电学院 Digital twin-based flood control and waterlogging prevention system and early warning method for whole elements of territorial universe

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140056645A1 (en) * 2010-07-28 2014-02-27 Richard C. Hogan Flood Control Method
CN110532952B (en) * 2019-08-30 2022-03-08 四川大学 Flood risk early warning and evacuation system based on GIS positioning technology
CN113902168B (en) * 2021-09-08 2022-11-01 长江信达软件技术(武汉)有限责任公司 Flood emergency evacuation method based on real-time crowd data fusion

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104851360A (en) * 2014-02-14 2015-08-19 杭州贵仁科技有限公司 Method and system for generating flood risk map
CN111652777A (en) * 2020-05-15 2020-09-11 长江勘测规划设计研究有限责任公司 Flood emergency risk avoiding method
CN114240119A (en) * 2021-12-10 2022-03-25 浙江水利水电学院 Digital twin-based flood control and waterlogging prevention system and early warning method for whole elements of territorial universe

Also Published As

Publication number Publication date
CN115358650A (en) 2022-11-18

Similar Documents

Publication Publication Date Title
CN115358650B (en) Flood disaster emergency risk avoiding and material real-time allocation method
Coles et al. Beyond ‘flood hotspots’: Modelling emergency service accessibility during flooding in York, UK
Gül et al. A combined hydrologic and hydraulic modeling approach for testing efficiency of structural flood control measures
CN104851360B (en) The generation method and system of a kind of flood risk mapping
Paz et al. Large‐scale modelling of channel flow and floodplain inundation dynamics and its application to the Pantanal (Brazil)
Mani et al. Flood hazard assessment with multiparameter approach derived from coupled 1D and 2D hydrodynamic flow model
CN114240119A (en) Digital twin-based flood control and waterlogging prevention system and early warning method for whole elements of territorial universe
CN108896117A (en) A kind of method of Remote Sensing To Hydrology station monitoring discharge of river
CN116070918B (en) Urban flood safety assessment and flood disaster prevention and control method
Ahmad et al. Hydrological modelling and flood hazard mapping of Nullah Lai
Seenu et al. Visualisation of urban flood inundation using SWMM and 4D GIS
Filipova Urban flooding in Gothenburg-A MIKE 21 study
Mangukiya et al. Integrating 1D and 2D hydrodynamic models for semi-arid river basin flood simulation
CN114528672A (en) Urban hydrological station network layout method and system based on 3S technology
Kobayashi et al. Development of a distributed rainfall‐run‐off/flood‐inundation simulation and economic risk assessment model
Athira et al. Flood modelling and inundation mapping of Meenachil River using HEC-RAS and HEC-HMS software
CN116523189A (en) Soil moisture content site planning method, device and storage medium considering hydrologic characteristics
Hashimoto et al. Multi-scale flooding hazards evaluation using a nested flood simulation model: case study of Jamuna River, Bangladesh
Kaoje Application of Geographical Information System Techniques in Urban Flood Risk Assessment and Vulnerability Mapping. A Case Study of Cardiff, Wales
Zaifoglu et al. Developing flood mitigation measures for the northern part of Nicosia
Konan Urban Flood modelling and Floodplain Mapping using ArcGIS, HEC-HMS and HEC-RAS in Abidjan city, Côte D’Ivoire–West Africa: Case study of the watershed of Bonoumin-Rivièra Palmeraie
Мартинюк et al. Identification of areas with potential significant flood risk using specialized software in the Vistula river basin within Ukraine
Ali Floodplain Determination and Management Case study: Wadi Watir-Sinai
Stone et al. Geospatial database and preliminary flood hydrology model for the Lower Colorado Basin
Zhou et al. Flood submergence simulation and risk analysis based on GIS in plain area

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