CN105741045B - A kind of flood risk dynamic analysis and display systems and method - Google Patents
A kind of flood risk dynamic analysis and display systems and method Download PDFInfo
- Publication number
- CN105741045B CN105741045B CN201610081204.1A CN201610081204A CN105741045B CN 105741045 B CN105741045 B CN 105741045B CN 201610081204 A CN201610081204 A CN 201610081204A CN 105741045 B CN105741045 B CN 105741045B
- Authority
- CN
- China
- Prior art keywords
- flood risk
- flood
- real
- user
- model
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Economics (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Marketing (AREA)
- Educational Administration (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
A kind of flood risk dynamic analysis and display systems and method, for building user oriented flood risk in-situ analysis and displaying, it is related to computer realm, the present invention is when structure is towards flood risk in-situ analysis and displaying, create flood risk dynamic analysis subsystem and flood risk real-time exhibition subsystem based on the system, initially set up data model, by the data interaction between web service Middleware implementation different layers, then pass through the flood risk real-time exhibition unit users Dynamic Display based on FLEX.The computational efficiency of existing flood risk mapping simulation model software is improved using the present invention, quick processing of the implementation model by the space expression of hour grid, and model data numerical procedure management data system is established, realize user's real-time online dynamic analysis flood risk, and to user's real-time exhibition.
Description
Technical Field
The invention relates to a flood risk dynamic analysis and display system and a flood risk dynamic analysis and display method, and belongs to the technical field of computers.
Background
The flood risk graph is compiled through scientific development and observation, a new idea of water control is practiced, flood prevention work is converted from flood control to flood management, a risk management system is established, important basic support of flood risk management is developed, and the flood risk graph is forward-looking, strategic and basic work of the country. The method has the advantages of actively promoting the establishment of flood risk maps in China, reasonably developing and utilizing the lands in flood control areas, enhancing flood consciousness of the whole people, promoting the scientization of flood control decisions, taking disaster reduction countermeasures according to local conditions, coordinating the water relations of people, and ensuring the stable development of society and national economy, and has important significance.
In order to promote the establishment work of flood risk maps in China, the water conservancy department starts 'national flood risk map establishment project trial (first stage)' work in 2008 and 'national flood risk map establishment project trial (second stage)' work in 2011, and the second stage trial project develops a large amount of work on the aspects of technical standards of flood risk map establishment, flood risk analysis software systems, flood risk map making and management platforms and the like of China, explores various flood risk map establishment methods, and initially forms the technical standards of flood risk map establishment and an organization management system. On the basis, a flood risk map compilation project of national key areas is started in 2013, and flood risk map compilation work is started in the range of each river basin of each province in the country. Despite the substantial research and discussion in pilot projects, and the substantial work foundation, there are still some problems in the large scale nationwide risk mapping.
Aiming at the problems that the foundation of flood risk map compilation and application is still imperfect, practical application is insufficient, compilation work is still irregular and the like, according to the requirements of comprehensively developing flood risk map compilation and promoting flood risk map application, management regulations, a standard system and a technical platform are further improved and perfected, the data management and application infrastructure of the flood risk map is enhanced, and a dynamic flood risk analysis and display system and a dynamic flood risk analysis and display method are provided.
Disclosure of Invention
The flood risk graph generating method and the flood risk graph generating system aim to dynamically generate the flood risk graph according to user requirements, carry out influence analysis and drawing in real time and quickly display results to users. The method breaks through the static property of the risk map compiled by the conventional method, improves the compiling efficiency of the flood risk map, simplifies the conventional compiling process, and achieves the effect of integrating flood risk analysis software, influence analysis software and drawing software which are necessary in the conventional compilation. The method has the advantages that the calculation efficiency of simulation model software in the conventional flood risk map compilation process is improved on the whole, an efficient grid geospatial processing program is developed, the rapid processing of the spatial expression of the grid on a small-hour basis by the model is realized by the parallel operation of the grid geospatial processing program and the flood simulation model software, a model data calculation scheme management data system is established, and parameter management, model calculation and real-time calculation result management are realized.
The technical scheme adopted by the invention for solving the technical problems is as follows: a flood risk dynamic analysis and display system and method.
In a first embodiment of the present invention, the flood risk dynamic analysis and display system includes a flood risk dynamic analysis subsystem and a flood risk real-time display subsystem.
The flood risk dynamic analysis subsystem adopts a distributed architecture, decomposes the calculation task of the flood risk model in a distributed mode, cooperates integrally, and drives the system to realize real-time online analysis and calculation of the model according to flood scheme conditions set online by a user.
The flood risk real-time display subsystem is used for rapidly displaying a two-dimensional result of dynamic flood risk analysis to a user in a chart mode in real time, and providing online information browsing query and resource sharing services for the user.
In another embodiment of the present invention, the flood risk dynamic analysis subsystem includes:
a service layer unit: the system is constructed on the basis of hardware, and a group of basic services are formed by a file or database server and an authentication server;
platform layer unit: the method mainly aims at a flood risk simulation model calculation processing platform and a data access engine based on an Oracle database, wherein the flood risk simulation model calculation processing platform covers a flood risk simulation model module and an hourly data geographic processing module, and the support of calculation and data access is provided for a business logic layer unit through the deployment of the flood risk simulation model calculation processing platform and the hourly data geographic processing module;
service logic layer unit: basic data management, model calculation management, model parameter management, statistical retrieval and thematic charting management are realized through a real-time scheme management component, and a flood risk simulation and evaluation display system of the business logic layer unit comprises a UI (user interface) interacting with a user;
a client layer unit: the system mainly provides a security access mechanism for a user to realize the security access of the system through a unified login mode and an authentication mode.
Preferably, the flood risk simulation model module realizes simulation of flood risk simulation by adopting a two-dimensional coupling mode, and the two-dimensional coupling mode is a mode combining a river network mode and a two-dimensional mode so as to simulate the whole dam-break flood risk development process.
In any of the above schemes, preferably, the Service logic layer unit includes a Web Service intermediate layer module, where the Web Service intermediate layer module includes six components, namely a real-time scheme management component, a basic data management component, a model calculation management component, a model parameter management component, a statistical retrieval component, and a thematic drawing management component, and message transmission, task scheduling, and data interaction between the flood risk dynamic analysis subsystem and the flood risk real-time display subsystem are realized through interaction of the six components.
In any of the above schemes, preferably, the flood risk simulation and evaluation presentation system includes a FLEX-based flood risk real-time presentation unit, and the real-time presentation unit is configured to provide, for a user, retrieval of a flood plan customized by the user, automatic submission of a real-time calculation result, and dynamic presentation of a real-time analysis calculation result.
In any of the above schemes, preferably, the flood risk dynamic analysis subsystem implements secure access of the system by using an automatic dumping and recovering mechanism of data, a multi-level security control mechanism, a user access database authorization mechanism, and a record (document) access security control.
In another embodiment of the present invention, a flood risk dynamic analysis and display method includes the following steps:
(1) A user sets parameters of a breach and a flood risk model through a flood risk dynamic analysis subsystem;
(2) Calling Web Service middleware through the B/S terminal, and storing the parameters into a database;
(3) The Web Service middleware calls and starts a grid data geographic processing module; the grid data processing module reads the parameters from the database, processes the parameters hour by hour, and stores the processed data into the database;
(4) And the B/S terminal polls the data in the database and dynamically displays the processed data in the flood risk real-time display subsystem.
The flood risk model is implemented using a two-dimensional coupling mode, wherein,
(1.1) when an upstream-downstream type coupling interface is configured, controlling iterative computation at the coupling interface through two indexes, namely maximum iterative computation times and a flow error control coefficient, wherein the flow error control coefficient is provided with a maximum value and a minimum value;
(1.2) when a lateral coupling interface is configured, configuring between every two adjacent sections;
(1.3) when the dike and the breach thereof are configured, if the height column number is 1, the development time of the breach is set to be 0; if the number of elevation rows is N (N > 2), the number of time for which the breach has progressed is N-2.
The Web Service middleware is used for realizing data interaction between a front end and a background of a system, and specifically comprises the following steps when the Web Service middleware is called:
(2.1) the B/S front end calls the Web Service middleware, and the computing Service of the Web Service middleware is automatically started;
(2.2) after starting, notifying a front-end user of the state of the calculation task according to the load condition of the current server, wherein the state comprises waiting, running and ending;
(2.3) executing corresponding command operation for transmitting the calculation task according to the parameters input by the user;
and (2.4) after the operation is finished, automatically closing the computing service and releasing the resources.
And the hourly processing adopts a harmonic inverse distance weight (AIDW) interpolation method to realize image interpolation.
The dynamic display comprises the following steps:
(4.1) displaying risk elements output by the model in real time according to model parameters input by a user, wherein the risk elements comprise a flow process and a water level of a specified cross section, a water surface line along the way, a basic risk thematic map and a submerging process;
(4.2) the user can perform real-time retrieval and browsing according to the calculation result, and the calculated inundation process is played back in a map interaction mode;
and (4.3) the user views, browses and plays back the inundation process by means of the customized flood scheme.
The technical scheme of the invention can achieve the following beneficial effects: in the technical means, the online real-time flood analysis evaluation capability is realized mainly through a distributed flood risk analysis model and the development of the high-efficiency exchange capability of flood risk analysis result data. The flood risk display subsystem is a platform for displaying flood risk real-time analysis and flood influence analysis results for users. The user can inquire and apply for various flood schemes in the platform and corresponding influence analysis and drawing results in real time on line by accessing the subsystem.
Drawings
FIG. 1 is a hierarchical diagram of a preferred embodiment of a flood risk dynamic analysis and presentation system according to the present invention;
FIG. 2 is a framework diagram of a preferred embodiment of a flood risk dynamic analysis subsystem in accordance with the present invention;
FIG. 3 is a block diagram of a preferred embodiment of a flood risk real-time presentation subsystem in accordance with the present invention;
FIG. 4 (a) is an overall analysis diagram of a preferred embodiment of the flood risk dynamic analysis subsystem according to the present invention;
FIG. 4 (b) is a parameter setting diagram of a preferred embodiment of the flood risk dynamic analysis subsystem according to the present invention;
FIG. 4 (c) is a real-time effect diagram of a preferred embodiment of the flood risk dynamic analysis subsystem according to the present invention;
FIG. 5 is an interface diagram of a preferred embodiment of the flood risk real-time presentation subsystem in accordance with the present invention;
Detailed Description
The invention is further described with reference to the drawings and the detailed description.
Referring to fig. 1, in a first embodiment of the present invention, a flood risk dynamic analysis and display system is first provided. The system comprises a flood risk dynamic analysis subsystem and a flood risk real-time display subsystem.
The flood risk dynamic analysis subsystem adopts a distributed architecture, decomposes the calculation task of the flood risk model in a distributed mode, cooperates integrally, and drives the system to realize real-time online analysis and calculation of the model according to flood scheme conditions set online by a user.
The flood risk real-time display subsystem is used for rapidly displaying a two-dimensional result of dynamic flood risk analysis to a user in a chart mode in real time, and providing online information browsing query and resource sharing services for the user.
Referring to fig. 2, the flood risk dynamic analysis subsystem is sequentially divided into a service layer unit, a platform layer unit, a business logic layer unit and a client layer unit from bottom to top.
A service layer unit: the system is constructed on the basis of hardware, and a group of basic services are formed by a file or database server and an authentication server;
platform layer unit: the method mainly aims at a flood risk simulation model calculation processing platform and a data access engine based on an Oracle database, wherein the flood risk simulation model calculation processing platform covers a flood risk simulation model module and an hourly data geographic processing module, and the support of calculation and data access is provided for a business logic layer unit through the deployment of the flood risk simulation model calculation processing platform and the hourly data geographic processing module;
service logic layer unit: basic data management, model calculation management, model parameter management, statistical retrieval and thematic charting management are realized through a real-time scheme management component, and a flood risk simulation and evaluation display system of the business logic layer unit comprises a UI (user interface) which interacts with a user;
a client layer unit: the system mainly provides a security access mechanism for realizing the security access of the system by a user through a unified login mode and an authentication mode.
In a more preferred embodiment, the flood risk dynamic analysis and presentation system first selects a model: in the scheme of the application, a two-dimensional coupling mode is adopted to realize simulation of flood risk simulation. The evolution process of flood occurrence is simulated and reproduced, and key hydrodynamic information required by flood disaster analysis, such as flood peak position, maximum submergence water depth, flow velocity, flow and other numerical information, is provided.
The two-dimensional coupling mode is used for combining a river network mode with a two-dimensional mode so as to simulate the whole dam-break flood risk development process. The one-dimensional river network mode provides river water level and flow velocity information, when the water level reaches a critical value, dam break occurs, flood enters a flood discharge area, flood discharge flow information is calculated according to a river network simulation result, and the calculation mode is divided into two modes: lateral type coupling, upstream and downstream type coupling. And controlling the water level parameter when the dam break occurs, wherein the shape of the dam can be set by a user. The two-dimensional flood risk analysis model mainly comprises a one-dimensional river network hydrodynamic model, a two-dimensional hydrodynamic model simulation and a two-dimensional connection.
The river network flow simulation selects a one-dimensional holy-Venn equation set as a control equation, wherein the control equation comprises two equations, namely a continuity equation and a momentum balance equation of water flow, and the expression is as follows:
in the above two formulas, Q is the flow rate, m 3 S; z is water level, m; q is the side inflow per unit water length, m 2 S; a is the area of the water passing cross section, m 2 (ii) a R is hydraulic radius, m; n roughness is obtained; a the momentum correction coefficient can be 1.0 in general; g is gravity acceleration, and is 9.8m/s 2 。
In the above two formulas, Q is the flow rate, m 3 S; z is water level, m; q is the side inflow per unit water length, m 2 S; a is aWater cross-sectional area, m 2 (ii) a R is hydraulic radius, m; n is roughness; a is a momentum correction coefficient, and can be 1.0 in general; g is gravity acceleration, and is 9.8m/s 2 。
Equations (2.1) and (2.2) are two first-order nonlinear equations, and the analytic solutions cannot be considered, but only numerical methods can be adopted for solving. Numerical solution first discretizes the equation, and common methods have explicit format and implicit format discretization. The explicit format has a large limitation on the time advance step length, and the change of the boundary information cannot affect all the grid points in one step, so that the explicit format is not suitable for solving the river network problem. The relation of the next time of the adjacent lattice point information is given by the hidden format dispersion, and the change of the boundary information can affect all lattice points in one step, so that the method is suitable for river simulation.
Adopting a Preissmann hidden format to disperse the two formulas and solving a domain [ x ] j ,x j+1 ]×[t n ,t n+1 ]The Taylor expansion is carried out on the function at a certain point in the space, and the first-order truncation is taken, so that the relation between the function value and the partial derivative value of the point and the function values of four surrounding points can be obtained. In short, the partial derivative value of a certain point in a small solution domain is determined by adopting a weighted average method, and the weight can be freely set. Here, the weight in the spatial direction is taken to be 1/2, and the time direction is taken to be 1/2 to ensure the stability of the calculation<θ<1。
Two-dimensional hydrodynamic models remain a hotspot problem in the research community, where the main problem is how to construct a suitable reconstruction method that enables the dry-wet boundary problem to be handled self-consistently. The model adopts a simple non-negative water depth reconstruction method, so that the reconstructed water depth and flow has stability at a dry-wet boundary and quality conservation is ensured. For the source term, the slope is changed into the flux term of the boundary, so that the processing is simpler and more convenient. The flux at the boundary is constructed using the Roe method, which allows for efficient capture of discontinuities. The time advance may be of selectable first and second order precision. And processing the resistance term by adopting an operator splitting method. Preliminary results show that the model has good still water simulation performance.
The conservation form of the two-dimensional shallow water equation is as follows,
in the above formula, t is time, and x and y are space coordinates; q is a conservation variable, f and g represent the flux of the conservation in the x, y directions, S is a source term, without taking into account the coriolis force term caused by the earth' S rotation, the viscous term caused by turbulence, and the shear stress term caused by wind. The source items comprise a land slope source item and a friction item caused by uneven terrain. These vectors are of a particular form such that,
wherein h is water depth and u is x Average velocity in the direction, v is the average velocity in the y direction, z b Is the ground elevation; q. q.s x And q is y Is the flux gradient in the x and y directions, q x Is equal to uh, q y Is equal to vh, g represents the acceleration of gravity, 9.8m/s 2 ；C f Is the coefficient of surface roughness, controlled by the Manning coefficient n and the depth of water h, C f ＝gn 2 /h 1/3 。
Two-dimensional hydrodynamic simulation is a common finite volume method, so the above formula is divided into volumes
Applying the Gaussian theorem, the above formula can be written as
In the above formula, Ω and Γ represent controlThe control volume unit is a polygon under the two-dimensional condition, and generally takes a triangle or a quadrangle for calculation, wherein the control volume unit is taken as the triangle. Each cell is labeled with i, n is the outward normal from the cell boundary, and has a component of (n) x ,n y ) The flux stream at the boundary can be written as
After river simulation and two-dimensional flood simulation can be carried out respectively, the key is how to connect the two models, namely the exchange of data. The river channel embankment is connected with the boundary of the flood discharge area, once the breach happens, the river channel provides water flow values for the boundary of the flood discharge area, and the key of model connection is to accurately describe information interaction of water flow inside and outside the breach.
The specific flux exchange formula is as follows:
wherein h is 1 ＝max(Z u ,Z d )-Z b ；h 2 ＝min(Z u ,Z d )-Z b And q is the flow rate.
A two-dimensional coupling flood risk simulation model software is mainly realized based on a file basis. Developed by Fortran. The Fortran development is more convenient and faster for developing a two-dimensional coupling model for scientific calculation.
And a two-dimensional joint solution calculation mode comprises upstream and downstream type interface configuration, lateral type coupling interface configuration, dike and break configuration thereof, output result configuration files and finally output result description, wherein the upstream and downstream type coupling interfaces are configured through a linker1 file, the lateral type coupling interfaces are configured through a linker2 file, the dike and break information thereof are configured through a bike file, and the coupling interface calculation result output configuration is carried out through a linker.
The one-two-dimensional joint solution calculation result comprises a one-dimensional model calculation result, a two-dimensional model calculation result and a coupling interface calculation result. Wherein, the calculation result of the one-dimensional model is stored in an OUT1D folder; the two-dimensional model calculation result is stored in an OUT2D folder; the results of the coupled interface calculations are stored in the OUT12D folder.
The coupling interface calculation results include idname. Linker1.Txt, idname. Linker2.Txt, and warning12d. Txt. Wherein Warning12D.txt records the relevant real-time calculation information which does not meet the flow convergence condition when the upstream and downstream coupling interfaces are subjected to iterative calculation; linker1.Txt records the water level and flow process of the id upstream and downstream coupling interface (named as name), wherein the flow is more than 0 and flows from two dimensions to one dimension, and the flow is less than 0 and flows from one dimension to two dimensions; linker2.Txt records the upstream section water level, the downstream section water level, the upstream section flow, the downstream section flow and the flood diversion flow process of the id-th lateral coupling interface combination (named as name), wherein the flood diversion flow is larger than 0 and indicates that the flow flows from one dimension to two dimensions, and the flood diversion flow is smaller than 0 and indicates that the flow flows from two dimensions to one dimension.
Preferably, a common two-dimensional coupling hydrodynamic model is selected to simulate the process of flood generation and development, and the selected model method is also a general model for flood risk simulation and analysis and has good practicability. The measurement of software efficiency mainly adopts two indexes of the capacity of processing grids and the data processing efficiency based on 1 hour flood interval to develop, improve and perfect the model.
The flood risk simulation model software is fully utilized to output the calculated interval based on 1-hour flood interval data, the data are cleaned, processed and compressed according to the real-time display requirement of the front end, the large data volume and the front end processing capacity are balanced through the software, so that the complex and large-amount calculation is finished in the background, and the data are subjected to front-end and back-end data communication through cleaning and compression processing.
Hourly data geo-processing software is primarily the geo-processing of current model calculation data. Data access (DataAccess), task management (TaskManager), a two-dimensional coupled flood risk simulation model (HydroMPM) and a GeoProcessor are mainly involved in the hourly data geographic processing process. In the geographic processing process, images simulating flood development at the front end are mainly generated in the background. Since the grids and nodes are discontinuous in information expression, a main link in the geographic processing is how to perform image interpolation. The image interpolation mainly comprises two aspects, namely: and performing boundary calculation and interpolating the feature graphs of the grids and the nodes in the boundary. On the interpolation algorithm, a harmonic inverse distance weight interpolation method is mainly selected and realized. The subsequent adoption of a new method depends on the performance of the interpolation method in the real-time flood risk evaluation.
The harmonic inverse distance weight algorithm has the following specific formula:
in the formula, Z i ，d i N, p are as defined for IDW; i is a sequence number (the sequence is from near to far) of the sample points in the interpolation search neighborhood from the interpolation points; k is a radical of formula i Is the azimuth harmonic weight coefficient of the sample point of serial number i, which represents the combined effect of the sample point of serial number i being masked by other sample points. Since the sample with sequence number 1 is closest to the interpolation point, no sample will have a masking effect on it, so k is set 1 ＝1。sin p θ ij Is a calculation formula of the shielding degree of the sample point of sequence number i by the sample point of sequence number j
Wherein j is more than or equal to 1 and less than or equal to i-1, theta ij Is the included angle (acute angle or right angle) between the connecting line of two sample points with serial numbers i and j and the central line (the over-interpolation point) thereof; a is a ij Is the connecting line angle between the two sample points of the serial number i, j and the interpolation point. According to the basic assumption of AIDW, when a ij At 360 DEG/n or more, the sequence number j samples have no masking effect on the sequence number samples, so that sin is specified at this time p θ ij The value is 1.
The Service layer based on Web Service mainly comprises a real-time scheme management component, a basic data management component, a model calculation management component, a model parameter management component, a statistical retrieval component and a thematic drawing management component. The information transmission, task scheduling and data interaction of the flood risk dynamic analysis subsystem and the flood risk real-time display subsystem are realized through the interaction of the six components
The real-time scheme management is mainly used for connecting flood risk graph UI interface requests in series and forwarding the requests to the flood risk graph model for calculation. The real-time scheme management is used as a circulation tandem connection person of the whole system and plays a core command role in the system.
The management of the breach mainly realizes the management of the information of the position of the breach. The method mainly comprises the following steps:
(1) Automatically acquiring the coordinate information of the breach node along the river channel according to the information of the breach position and the breach length, and forming a breach node table which can be identified by a model;
(2) Automatically identifying the upstream and downstream sections according to the position and the length of the burst opening;
(3) According to the different shapes of the burst generated by the burst development process (t 1, t2, t3, t4 …) and the burst node, the formation and development of different burst shapes such as rectangular, trapezoidal and triangular burst are supported.
The monitoring point management mainly comprises the setting of one-dimensional and two-dimensional grid monitoring points. The monitoring point setting transmits the coordinates of the monitoring points to the flood risk graph through scheme management so as to be output by the model. The monitoring point can be arranged in a plurality of numbers. The specific process is as follows:
(1) The scheme management stores the monitoring point information into a monitoring point (MonitorPoint) database through a flood risk graph.
(2) And starting a model scheduling program by the scheme, reading monitoring point information by the model pair through the scheme ID, and then generating the model hour by hour.
The management of monitoring sections is similar to the management of monitoring points. The monitoring section has a scheme manager to store the serial number in the MonitorSection. The specific method is consistent with the monitoring point management.
Grid management is mainly to process grid data of dam break areas and flood discharge areas. The method mainly comprises the following steps:
(1) And vectorizing grid data. The mesh data is the fundamental data for model analysis. After grid data is subjected to grid processing through a special grid processing tool, the grid data needs to be subjected to vectorization processing, and meanwhile, textual grid files which can be identified by a model are processed to enable the model to perform simulation calculation.
(2) And (4) classified management of grid data. The density of the grids is directly related to the simulation efficiency, and the project realizes data management on the grids with different resolutions so as to save and reference different grids according to the selection of a user.
(3) Data management and interconversion of linear and polygonal meshes. Both the linear grid and the polygonal grid are basic background data of the current system and participate in system data calculation and information extraction.
The grid boundary data is extracted by scanning the grid data boundary conditions. Is the most accurate boundary for delineating the gridded flood discharge area. The mesh boundary data is polygon data and may be composed of a plurality of polygon data.
The river bank line is basic background data of a project and is directly put in storage as initialization data. River section is the key data for a project. To facilitate real-time flood simulation of the system. The node data includes a set of node data along the bank in addition to a set of mesh node data. Both data sets are used as basic background and participate in real-time flood simulation calculation.
Furthermore, the flood risk dynamic analysis and display system comprises a flood risk real-time display unit based on FLEX, wherein the real-time display unit is used for providing retrieval of a user customized flood scheme, automatic submission of real-time calculation results and dynamic display of real-time analysis calculation results for users.
Furthermore, the flood risk dynamic analysis subsystem adopts an automatic dumping and recovering mechanism of data, a multi-level security control mechanism, a user access database authorization mechanism and record (literature) access security control to realize the security access of the system.
An important task of dynamic flood risk analysis and real-time display system construction is to integrate various related data resources by taking a real-time flood risk scheme as a basis and combining flood risk graph compilation specification so as to form a uniform, complete and standard flood risk graph database. The database comprises various spatial information and non-spatial information such as grids, grid base maps and the like, and corresponding metadata description information, and has the function of providing standard and reasonable data for a flood risk map query analysis system.
Referring to fig. 3, a frame diagram of a flood risk real-time display subsystem of the flood risk dynamic analysis and display system provided by the present invention,
the flood risk real-time display subsystem comprises a business knowledge module, an information announcement module, a map service module, an application case module, a risk map achievement module, a data downloading module, a technical service module and a user management module.
The system access users are divided into non-registered users, registered users and administrator users. The user management mainly aims at that map service and partial download resources in the system are internal data, is not suitable for being fully opened, and can be checked and downloaded by users with corresponding authorities. The registered user inputs a user name and a password through the login entrance, and can check the online service resources and download the required technical data after authority authentication. And the system administrator user is responsible for the flood risk map management and user audit of the website. Non-registered users may access functions other than map services and download public material.
And registering a user account of the login system. And inputting account login information and user identity data such as a user name and a password for registering a real name, a user name, a password, a mailbox, a telephone and an online map service resource. The newly registered account can be used only after the audit of the flood risk map administrator is passed. The auditing result is sent to the mailbox filled in during registration through the mail. The registered user is required to comply with the "terms of service" and "user must know" of the resource servicing system. Wherein, each module of the system has the following functions
The user management module is used for providing a login platform for a user, the platform has corresponding use and management authority according to the type of the user, the user can edit and modify user information and reset a password, the user is divided into a non-registered user, a registered user and an administrator user, and the user can check information such as online service resources, downloaded data and the like only after obtaining authority authentication.
The business knowledge module is used for providing related basic knowledge and browsing and downloading operation of legal and legal information for users, and an administrator can manage, release, edit and delete information bulletins;
the information bulletin module is used for providing browsing and downloading operations of headlines, bulletins and dynamic information of bulletin release for a user, providing a news information list related to the water conservancy industry for the user, displaying all water conservancy important news and supporting fuzzy search of titles;
the map service module is used for providing service use, service preview, service registration, service application and service retrieval for a user, and displaying the distribution condition of each risk point by adopting an arcgis map, wherein the related information of the risk comprises the following steps: the names of the risk points, the belonged watersheds, the belonged provinces and the number of the schemes.
The application case module is used for providing the application cases for the user to carry out operations such as retrieval, viewing and the like;
the risk map achievement module is used for providing a basin, risk points and a scheme for a user, wherein the basin is divided into a Yangtze river basin, a Taihu lake basin, a yellow river basin, a Zhujiang river basin, a sea river basin, a Songliao basin and a Huaihe river basin, and the scheme is divided into: land utilization planning, flood control and disaster reduction, and flood insurance; the flood is divided into: river and lake flood, flood in stagnant flood storage area, reservoir flood and urban flood; the risk map is divided into: a flood submergence range map, a flood submergence depth map, a flood submergence duration map, and a flood flow velocity map.
The data downloading module is used for providing downloading of software resources and technical standard resources for a user;
the technical service module provides relevant knowledge, training materials, data summarization, an expert database and a service guide for a user.
In another embodiment of the present invention, the present invention provides a flow chart of a flood risk dynamic analysis and display method, wherein the method comprises the following steps:
(1) A user sets parameters of a breach model and a flood risk model through a flood risk dynamic analysis subsystem;
(2) Calling Web Service middleware through the B/S terminal, and storing the parameters into a database;
(3) The Web Service middleware calls and starts a grid data geographic processing module; the grid data processing module reads the parameters from the database, processes the parameters hour by hour, and stores the processed data into the database;
(4) And the B/S terminal polls the data in the database and dynamically displays the processed data in the flood risk real-time display subsystem.
Further, the flood risk model is implemented using a two-dimensional coupling mode, wherein,
(1.1) when an upstream and downstream type coupling interface is configured, controlling iterative computation at the coupling interface through two indexes, namely maximum iterative computation times and a flow error control coefficient, wherein the flow error control coefficient is provided with a maximum value and a minimum value;
(1.2) when a lateral coupling interface is configured, configuring between every two adjacent sections;
(1.3) when the embankment and the breach thereof are configured, if the height column number is 1, the time for developing the breach is set to be 0; if the number of elevation columns is N (N > 2), the number of the time for which the breach develops is N-2.
Furthermore, the Web Service middleware is configured to implement data interaction between a front end of the system and a background, and when the Web Service middleware is called, the method specifically includes:
(2.1) the B/S front end calls the Web Service middleware, and the computing Service of the Web Service middleware is automatically started;
(2.2) after starting, notifying a front-end user of the state of the calculation task according to the load condition of the current server, wherein the state comprises waiting, running and ending;
(2.3) executing corresponding command operation of transmitting the calculation task according to the parameters input by the user;
and (2.4) after the operation is finished, automatically closing the computing service and releasing the resources.
Preferably, the dynamic presentation includes:
(4.1) displaying risk elements output by the model in real time according to model parameters input by a user, wherein the risk elements comprise a flow process and a water level of a specified cross section, a water surface line along the way, a basic risk thematic map and a submerging process;
(4.2) the user can perform real-time retrieval and browsing according to the calculation result, and the calculated inundation process is played back in a map interaction mode;
and (4.3) the user views, browses and plays back the inundation process by means of the customized flood scheme.
(4.4) outputting one-dimensional information and two-dimensional information of the flood risk model in real time, wherein,
the one-dimensional information includes: a breach flow process, a section flow process, a single section water level, and a section on-way water level.
The two-dimensional information includes: clicking a section on a map, and displaying a section flow process line and a water level change process line; and (5) selecting and displaying submerged water depth maps and flow velocity maps at different moments according to time axes on a two-dimensional surface result.
The technical scheme of the invention can achieve the following beneficial effects: in the technical means, the online real-time flood analysis evaluation capability is realized mainly through a distributed flood risk analysis model and the development of the high-efficiency exchange capability of flood risk analysis result data. The flood risk display subsystem is a platform for displaying flood risk real-time analysis and flood influence analysis results for users. The user can inquire and apply for various flood schemes in the platform and corresponding influence analysis and drawing results in real time on line by accessing the subsystem.
Claims (5)
1. A flood risk dynamic analysis and display method comprises the following steps:
(1) Setting flood risk model parameters by a user through a flood risk dynamic analysis subsystem; the flood risk dynamic analysis subsystem adopts a distributed architecture, decomposes the calculation task of a flood risk model in a distributed mode, cooperates integrally, and realizes real-time online analysis and calculation of the model according to the flood risk dynamic analysis subsystem arranged online by a user;
the flood risk dynamic analysis subsystem comprises:
a service layer unit: the system is constructed on the basis of hardware, and comprises a group of basic service databases consisting of a database server and an authentication server;
platform layer unit: the objects contained in the platform layer unit comprise a flood risk simulation model calculation processing platform and a data access engine based on an Oracle database; the flood risk simulation model calculation processing platform covers a flood risk simulation model module and an hourly data geographic processing module, and provides calculation and data access support for a business logic layer unit through the deployment of the flood risk simulation model module and the hourly data geographic processing module;
service logic layer unit: the system comprises a WebService intermediate layer module and a flood risk simulation and evaluation display system; the flood risk simulation and evaluation display system comprises a UI (user interface) for interacting with a user and a flood risk real-time display unit based on FLEX (flexible exchange); the flood risk real-time display unit is used for providing retrieval of a user customized flood scheme, automatic submission of real-time calculation results and dynamic display of real-time analysis calculation results for users; the WebService intermediate layer module comprises a real-time scheme management component, a basic data management component, a model calculation management component, a model parameter management component, a statistical retrieval component and a thematic charting management component, and the real-time scheme management component is used for realizing the interaction of the six components of basic data management, model calculation management, model parameter management, statistical retrieval and thematic charting management and realizing the message transmission, task scheduling and data interaction of the flood risk dynamic analysis subsystem and the flood risk real-time display subsystem; the flood risk real-time display subsystem is used for rapidly displaying a two-dimensional result of dynamic flood risk analysis to a user in a chart mode in real time, and providing online information browsing query and resource sharing services for the user;
a client layer unit: a security access mechanism of a unified login mode and a verification mode is provided for a user, and the security access of the system is realized;
(2) Calling a WebService intermediate layer module through a B/S terminal, and storing the flood risk model parameters into an Oracle database;
(3) The WebService intermediate layer module calls and starts an hourly data geographic processing module; the hourly data geographic processing module reads the flood risk model parameters from the Oracle database, performs hourly processing on the flood risk model parameters, and stores the processed data into the Oracle database;
(4) The B/S terminal polls the data in the Oracle database and dynamically displays the processed data in the flood risk real-time display subsystem;
the flood risk model is implemented using a two-dimensional coupling mode, wherein,
(1.1) when an upstream-downstream type coupling interface is configured, controlling iterative computation at the coupling interface through two indexes, namely maximum iterative computation times and a flow error control coefficient, wherein the flow error control coefficient is provided with a maximum value and a minimum value;
(1.2) when a lateral coupling interface is configured, configuring between every two adjacent sections;
(1.3) when the dike and the breach thereof are configured, if the height column number is 1, the development time of the breach is set to be 0; if the number of elevation columns is N (N > 2), the time for the development of the breach is N-2.
2. The flood risk dynamic analysis and presentation method according to claim 1, wherein the WebService intermediate layer module is configured to implement data interaction between a front end and a background of a system, and when the WebService intermediate layer module is called, the method specifically includes:
(2.1) the B/S front end calls the WebService middle layer module, and the calculation service of the WebService middle layer module is automatically started;
(2.2) after starting, notifying a front-end user of the state of the calculation task according to the load conditions of the current database server and the current authentication server, wherein the state comprises waiting, running and ending;
(2.3) executing corresponding command operation for transmitting the calculation task according to the flood risk model parameters input by the user;
and (2.4) after the operation is finished, automatically closing the computing service and releasing the resources.
3. The flood risk dynamic analysis and presentation method according to claim 1, wherein the hourly processing employs a harmonic inverse distance weight interpolation method to realize image interpolation, wherein the harmonic inverse distance weight interpolation algorithm formula is:
wherein Z is an interpolation point estimate; zi is the ith sample point observation; di is the euclidean distance between the interpolation point and the ith sample point; n is the number of samples used to estimate the interpolated point value; p is a power exponent; i is an interpolation searchSequencing serial numbers of the sample points in the neighborhood far from the interpolation point are arranged from near to far; k is a radical of i Is the azimuth harmonic weight coefficient of the sample point of the serial number i, which represents the comprehensive effect of the sample point of the serial number i shielded by other sample points; since the sample with sequence number 1 is closest to the interpolation point, no sample will have a masking effect on it, so k is set 1 ＝1；sin p θ ij Is a calculation formula of the shielding degree of the sample point of sequence number i by the sample point of sequence number j
Wherein j is more than or equal to 1 and less than or equal to i-1, theta ij The included angle between the connecting line of the two sample points of the serial numbers i and j and the central line of the interpolated point is acute angle or right angle; a is ij Is the connecting line angle between two sample points of sequence number i, j and the interpolation point, according to the basic assumption of harmonic and inverse distance weight, when a ij When the value is more than or equal to 360 DEG/n, the sample point of sequence number j has no shielding effect on the sample point of sequence number, so that sin at this time is specified p θ ij The value is 1.
4. The flood risk dynamic analysis and presentation method of claim 2, wherein the dynamic presentation comprises:
(4.1) displaying risk elements output by the flood risk model in real time according to flood risk model parameters input by a user, wherein the risk elements comprise a flow process and a water level of a specified cross section, a water surface line along the way, a basic risk thematic map and a submerging process;
(4.2) the user can perform real-time retrieval and browsing according to the calculation result, and the calculated inundation process is played back in a map interaction mode;
and (4.3) the user views, browses and plays back the flooding process on the customized flooding scheme through the flood risk real-time display subsystem.
5. The flood risk dynamic analysis and presentation method of claim 4, wherein the dynamic presentation further comprises:
(4.4) outputting one-dimensional information and two-dimensional information of the flood risk model in real time, wherein,
the one-dimensional information includes: a breach flow process, a section flow process, a single section water level and a section on-way water level;
the two-dimensional information includes: clicking a section on a map, and displaying a section flow process line and a water level change process line; and (5) selecting and displaying submerged water depth maps and flow velocity maps at different moments according to time axes on a two-dimensional surface result.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610081204.1A CN105741045B (en) | 2016-02-05 | 2016-02-05 | A kind of flood risk dynamic analysis and display systems and method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610081204.1A CN105741045B (en) | 2016-02-05 | 2016-02-05 | A kind of flood risk dynamic analysis and display systems and method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105741045A CN105741045A (en) | 2016-07-06 |
| CN105741045B true CN105741045B (en) | 2018-03-20 |
Family
ID=56244975
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610081204.1A Expired - Fee Related CN105741045B (en) | 2016-02-05 | 2016-02-05 | A kind of flood risk dynamic analysis and display systems and method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105741045B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106599140A (en) * | 2016-12-06 | 2017-04-26 | 河海大学 | GIS-based fast processing method of flood risk factors |
Families Citing this family (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107229742A (en) * | 2017-06-20 | 2017-10-03 | 华中科技大学 | A kind of method that city easily flood point is determined based on remote sensing big data |
| CN107895076B (en) * | 2017-11-10 | 2021-04-27 | 长江大学 | Non-penetrating river channel modeling method and system based on target |
| CN108256193B (en) * | 2018-01-10 | 2018-11-16 | 中国水利水电科学研究院 | A kind of river flood due to dam-break one, two-dimensional numerical model coupling analytical method |
| CN108062453B (en) * | 2018-01-12 | 2021-06-29 | 河南省水利勘测设计研究有限公司 | Flood high-efficiency simulation and high-fidelity visual dynamic display method for water conservancy information system |
| CN109242300B (en) * | 2018-08-31 | 2022-08-26 | 创新先进技术有限公司 | Method and device for displaying service index value and computer equipment |
| CN109933637B (en) * | 2019-03-19 | 2021-06-18 | 北京航天宏图信息技术股份有限公司 | Flood risk dynamic display and analysis system |
| CN110070619B (en) * | 2019-04-09 | 2021-07-06 | 长江空间信息技术工程有限公司(武汉) | Three-dimensional visualization method for river water surface considering water surface gradient |
| CN110111538B (en) * | 2019-04-30 | 2021-06-25 | 中国电建集团北京勘测设计研究院有限公司 | Dynamic monitoring, early warning and analyzing system for mountain torrent disasters |
| CN111027175B (en) * | 2019-11-06 | 2021-09-24 | 中国地质大学（武汉） | Method for evaluating social and economic influences of flood based on coupling model integrated simulation |
| CN111080060A (en) * | 2019-11-11 | 2020-04-28 | 天津大学 | Method for constructing embankment flood risk evaluation index system |
| CN112084280B (en) * | 2020-09-04 | 2023-07-21 | 广州南方智能技术有限公司 | Multi-scale terrain cutting and splicing method |
| CN117114428B (en) * | 2023-10-25 | 2024-01-30 | 国网山西省电力公司电力科学研究院 | Meteorological disaster analysis and early warning method for power equipment |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104408900A (en) * | 2014-11-10 | 2015-03-11 | 柳州师范高等专科学校 | Dynamic optimization based neural network flood warning device and method |
| CN104851360A (en) * | 2014-02-14 | 2015-08-19 | 杭州贵仁科技有限公司 | Method and system for generating flood risk map |
-
2016
- 2016-02-05 CN CN201610081204.1A patent/CN105741045B/en not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104851360A (en) * | 2014-02-14 | 2015-08-19 | 杭州贵仁科技有限公司 | Method and system for generating flood risk map |
| CN104408900A (en) * | 2014-11-10 | 2015-03-11 | 柳州师范高等专科学校 | Dynamic optimization based neural network flood warning device and method |
Non-Patent Citations (4)
| Title |
|---|
| 基于Web Services的洪水风险图信息管理系统设计与实现;张强;《中国学位论文全文数据库库》;20131231;第14-39页 * |
| 实时动态耦合模型及其在洪水风险图中的应用;付成威等;《水利水运工程学报》;20131031(第5期);第32-38页 * |
| 洪水风险图制图方法及系统;徐美等;《中国防汛抗旱》;20150831;第25卷(第4期);第8-13页 * |
| 顾及方向遮蔽性的反距离权重插值法;李正泉等;《测绘学报》;20150131;第44卷(第1期);第91-96页 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106599140A (en) * | 2016-12-06 | 2017-04-26 | 河海大学 | GIS-based fast processing method of flood risk factors |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105741045A (en) | 2016-07-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105741045B (en) | A kind of flood risk dynamic analysis and display systems and method | |
| Pollock | User guide for MODPATH Version 7—A particle-tracking model for MODFLOW | |
| Sukoriansky et al. | An analytical theory of the buoyancy–Kolmogorov subrange transition in turbulent flows with stable stratification | |
| Oliveira et al. | OPENCoastS: An open-access service for the automatic generation of coastal forecast systems | |
| Varduhn et al. | Multi-resolution models: Recent progress in coupling 3D geometry to environmental numerical simulation | |
| Yang et al. | A digital information system for cultural landscapes: The case of Slender West Lake scenic area in Yangzhou, China | |
| Anselin et al. | Open source software for spatial data science | |
| Whiteaker et al. | From a NEXRAD rainfall map to a flood inundation map | |
| Du et al. | Ensemble data assimilation applied to an adaptive mesh ocean model | |
| Nessel | The place for information models in landscape architecture, or a place for landscape architects in information models | |
| Ribault et al. | Simulation processes in the cloud for emergency planning | |
| Minocha et al. | Reservoir Assessment Tool version 3.0: a scalable and user-friendly software platform to mobilize the global water management community | |
| Röber et al. | Visualization of climate science simulation data | |
| Rustom et al. | Estimating productivity using simulation: A case study of Gaza beach embankment protection project | |
| Lahayon et al. | GENEALOGICAL AGGREGATION AND ANCESTRAL DOMAIN MANAGEMENT: A WEB-BASED GIS IMPLEMENTATION IN A DEVELOPING COUNTRY | |
| Stainforth et al. | Climateprediction. net: a global community for research in climate physics | |
| Shaw et al. | A nonhydrostatic primitive-equation model for studying small-scale processes: an object-oriented approach | |
| Kurzbach | Flood modeling in spatial data and grid infrastructures | |
| da Silva Marques | RiverCure Portal: Exploring Geographic Features on Context Definition and Integration with the HiSTAV Hydrometric Tool | |
| Saghafian | Implementation of a distributed hydrologic model within Geographic Resources Analysis Support System (GRASS) | |
| Chen | [Retracted] Construction of a Large Visual Landscape Design Platform considering Low‐Carbon Environmental Protection Based on Nanomaterials | |
| Hluchy et al. | Grid based flood prediction virtual organization | |
| Miura et al. | Development and Validation of a Three‐Dimensional Variably Saturated Flow Model for Future Water Resource Assessment at a Global Scale—Targeting Saturated Groundwater Flow in Plains | |
| Vieira | Towards a decision support system for flood management in a river basin | |
| Bournaski et al. | Incorporation of computational models to support water resources management |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB03 | Change of inventor or designer information |
Inventor after: Xiang Liyun Inventor after: Zeng He Inventor after: Li Kaifeng Inventor after: Zhang Xichuan Inventor after: Liu Shu Inventor after: Xu Mei Inventor after: Zhang Hongping Inventor after: Hu Changwei Inventor after: Zheng Jingwei Inventor after: Sha Ruihua Inventor after: Song Lixiang Inventor after: Lan Hualin Inventor before: Xiang Liyun Inventor before: Zeng He Inventor before: Liu Shu Inventor before: Xu Mei Inventor before: Song Lixiang Inventor before: Zhang Xichuan Inventor before: Sha Ruihua Inventor before: Lan Hualin Inventor before: Zhang Hongping Inventor before: Hu Changwei |
|
| COR | Change of bibliographic data | ||
| CB03 | Change of inventor or designer information |
Inventor after: Xiang Liyun Inventor after: Cang Wenbin Inventor after: Zhang Xichuan Inventor after: Zeng He Inventor after: Li Kaifeng Inventor after: Liu Shu Inventor after: Xu Mei Inventor after: Zhang Hongping Inventor after: Hu Changwei Inventor after: Yao Yonghui Inventor after: Zheng Jingwei Inventor after: Song Lixiang Inventor after: Sha Ruihua Inventor before: Xiang Liyun Inventor before: Zeng He Inventor before: Li Kaifeng Inventor before: Zhang Xichuan Inventor before: Liu Shu Inventor before: Xu Mei Inventor before: Zhang Hongping Inventor before: Hu Changwei Inventor before: Zheng Jingwei Inventor before: Sha Ruihua Inventor before: Song Lixiang Inventor before: Lan Hualin |
|
| CB03 | Change of inventor or designer information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180320 Termination date: 20210205 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |