CN115293469A - Urban flood control and drainage risk prediction method - Google Patents
Urban flood control and drainage risk prediction method Download PDFInfo
- Publication number
- CN115293469A CN115293469A CN202211231275.7A CN202211231275A CN115293469A CN 115293469 A CN115293469 A CN 115293469A CN 202211231275 A CN202211231275 A CN 202211231275A CN 115293469 A CN115293469 A CN 115293469A
- Authority
- CN
- China
- Prior art keywords
- flood
- water level
- current city
- day
- city
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 46
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 177
- 238000012937 correction Methods 0.000 claims abstract description 62
- 239000002689 soil Substances 0.000 claims abstract description 11
- 238000012545 processing Methods 0.000 claims abstract description 5
- 230000002354 daily effect Effects 0.000 claims description 40
- 238000001556 precipitation Methods 0.000 claims description 27
- 239000011159 matrix material Substances 0.000 claims description 20
- 239000010865 sewage Substances 0.000 claims description 15
- 239000000126 substance Substances 0.000 claims description 13
- 230000003203 everyday effect Effects 0.000 claims description 6
- 230000002093 peripheral effect Effects 0.000 claims description 6
- 238000010606 normalization Methods 0.000 claims description 3
- 238000012163 sequencing technique Methods 0.000 claims description 2
- 238000003672 processing method Methods 0.000 abstract 1
- 238000010521 absorption reaction Methods 0.000 description 4
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229960001948 caffeine Drugs 0.000 description 1
- 238000010924 continuous production Methods 0.000 description 1
- 230000000994 depressogenic effect Effects 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- RYYVLZVUVIJVGH-UHFFFAOYSA-N trimethylxanthine Natural products CN1C(=O)N(C)C(=O)C2=C1N=CN2C RYYVLZVUVIJVGH-UHFFFAOYSA-N 0.000 description 1
Images
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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- 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
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Data Mining & Analysis (AREA)
- Tourism & Hospitality (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Mathematical Physics (AREA)
- Development Economics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Operations Research (AREA)
- Health & Medical Sciences (AREA)
- Entrepreneurship & Innovation (AREA)
- Quality & Reliability (AREA)
- Bioinformatics & Computational Biology (AREA)
- Probability & Statistics with Applications (AREA)
- Educational Administration (AREA)
- Evolutionary Biology (AREA)
- Primary Health Care (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Algebra (AREA)
- General Health & Medical Sciences (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention relates to the technical field of electric digital data processing, in particular to a method for predicting urban flood control and drainage risks. The method is a digital processing method which is designed by computer assistance and is particularly suitable for flood disaster early warning, the method is characterized in that after the historical water storage variation and the water level variation are compared to obtain the water storage capacity corresponding to the unit water level variation of the current city, the prediction of the flood disaster occurrence probability of the current city at the future moment is completed based on the water storage capacity corresponding to the unit water level variation, then, historical factor correction parameters are additionally constructed according to the flood disaster situation and the water level situation at the same moment in the same year, land factor correction parameters are constructed according to the soil occupation ratio of the city, and the land factor correction parameters are constructed according to the land difference situation and the flood disaster occurrence situation of the current city and other cities, the flood disaster prediction is carried out after the obtained flood disaster occurrence probability is corrected by using the three correction parameters, and the flood prediction accuracy of the flood disaster is improved.
Description
Technical Field
The invention relates to the technical field of electric digital data processing, in particular to a method for predicting urban flood control and drainage risks.
Background
In recent years, along with the aggravation of climate change, extreme climate events occur frequently, wherein flood disasters not only bring huge economic losses to cities, but also seriously threaten the safety of the cities and residents. In order to resist the influence of flood disasters on cities, flood control and drainage risks need to be predicted, so that the loss caused by the flood control and drainage risks is reduced.
The conventional flood control and drainage risk prediction method generally realizes flood control and drainage risk prediction of the current city by using historical hydrological data of the current city, the prediction method is simple, and the considered data source is single, so that a space for further improving the prediction accuracy exists.
Disclosure of Invention
The invention provides a method for predicting urban flood control and drainage risks, which is used for improving the accuracy of predicting the urban flood disasters, and adopts the following technical scheme:
the invention discloses a method for predicting urban flood control and drainage risks, which comprises the following steps:
determining the water storage variable quantity of the current city for the set historical days today according to the actual precipitation, daily average sewage discharge and daily drainable quantity of the current city per day historically, and determining the water level variable quantity of the current city for the set historical days today according to the historical water level of the current city per day ending moment historically;
determining the water storage capacity corresponding to the current city unit water level change according to the water storage change and the water level change;
determining the predicted water level of the current city every day after the current day according to the water storage capacity corresponding to the unit water level change of the current city, the predicted precipitation amount of the future every day, the daily drainable amount and the daily average sewage discharge amount, thereby determining the flood occurrence probability at the end time of each day after the current day;
determining a set date section comprising the reference day by taking the number of days corresponding to the predicted flood occurrence probability as the reference day, and determining the flood disaster occurrence condition and the water level condition in the set date section in each year in the set historical number of years to obtain historical factor correction parameters for correcting the flood occurrence probability;
determining a land factor correction parameter for correcting the flood disaster occurrence probability according to the occupation ratio of the soil area in the current city in the total area of the city;
determining a terrain factor correction parameter for correcting the flood occurrence probability according to the terrain difference between the current city and other cities and the flood occurrence conditions of other cities;
and correcting the flood occurrence probability according to the determined historical factor correction parameters, the land factor correction parameters and the relief factor correction parameters, and completing flood prediction according to the corrected flood occurrence probability.
The beneficial effects of the invention are as follows:
according to the method, after the flood occurrence probability at the future moment is predicted based on the water storage capacity corresponding to the unit water level change of the current city, historical factor correction parameters are additionally constructed according to the flood disaster situation and the water level situation at the same moment in the past year, land factor correction parameters are constructed according to the soil proportion of the city, and the land difference situation of the current city and other cities and the flood occurrence situation are constructed, the flood prediction is carried out after the obtained flood occurrence probability is corrected by using the three-aspect correction parameters, and the accuracy of the flood prediction is improved.
Further, the method for determining the water storage capacity corresponding to the current city unit water level change comprises the following steps:
wherein, the first and the second end of the pipe are connected with each other,the water storage capacity corresponding to the current city unit water level change,indicating the set historical number of days since today,representing the amount of change in the impounded water for the j day before the current city today,the historical water level value representing the current city at yesterday's end time,indicating that the current city is the first to this dayThe historical water level value at the end of the day,representing the actual precipitation of the current city on day j before today,indicating the daily average sewage discharge amount of the current city,representing the daily drainable amount of the current city.
Further, the predicted water level of each day after the current city is:
wherein, the first and the second end of the pipe are connected with each other,indicating the predicted water level at the end of day I after the current city,representing the predicted impounded water change amount of the current city on the ith day after today,the historical water level value representing the current city at yesterday's end time,representing the predicted precipitation for the current city on day i after this day,represents the daily average sewage discharge amount of the current city,representing the daily drainable amount of the current city.
Further, the flood occurrence probability at the end of each day after today is:
wherein the content of the first and second substances,the flood probability of the current city at the end of the day I after the current day,、andrespectively are the height values of a normal water level line, a warning water level line and a dangerous water level line,predicting a weight for the flood probability, wherein,For the predicted water level at the end of day I after the current city today,the numerical value in parentheses is 0 when the numerical value in parentheses is not positive, and the numerical value in parentheses is self-value when the numerical value in parentheses is positive.
Further, the method for obtaining the historical factor correction parameter comprises the following steps:
wherein, the first and the second end of the pipe are connected with each other,and withRespectively representing a first historical factor correction parameter and a second historical factor correction parameter,indicating that the number of the historical year copies is set,indicating that the water level value is less than the predicted water level in the set date section in the set historical year numberAnd the number of yearly flood disasters occur,indicating that the water level value is greater than the predicted water level in the set date period in the set historical number of yearsAnd the annual number of flood disasters do not occur.
Further, the method for obtaining the land factor correction parameter comprises the following steps:
and (4) counting the proportion Q of the soil area in the total area in the current city, and taking the proportion Q as a land factor correction parameter.
Further, the method for obtaining the terrain factor correction parameter comprises the following steps:
dividing the periphery of the city into peripheral areas with the number of the set areas, sequentially making a difference between the altitude of the city area and the altitude of each peripheral area to obtain the altitude difference values of the number of the set areas, then performing normalization processing, and sequencing according to the sequence from small to large to obtain a terrain sequence;
determining the terrain sequence of the current city and the terrain sequences of other cities according to a terrain sequence acquisition method, and then calculating the terrain advantages and disadvantages parameters of the current city compared with other cities:
wherein the content of the first and second substances,representing the terrain superiority and inferiority parameters of the current city compared with other cities,andrespectively representing the current cityThe maximum and minimum values of the elements in the potential matrix,andrespectively representing the maximum and minimum values of the elements in the terrain matrix of the other cities,the number of elements whose values fall into the latter half of the value range of the elements in the terrain matrix of the current city,representing the number of elements with values falling into the latter half of the value range of the elements in the terrain matrix of other cities;
selecting other cities with the set number of Y, and counting the number of Y to obtain each cityThe number of other cities which are larger than 1 and corresponding to other cities in flood disasterAnd anThe number of the value is less than 1 and corresponds to the number of other cities without flood disastersTo in order toAndas a relief factor correction parameter for correcting the flood occurrence probability.
Further, the method for completing flood prediction according to the corrected flood occurrence probability comprises:
firstly determining the corrected flood occurrence probability:
wherein the content of the first and second substances,indicating the flood probability of the revised current city at the end time of the day I after the current day,and withRespectively representing a first historical factor correction parameter and a second historical factor correction parameter,the land factor correction parameter is represented by a land factor correction parameter,andrespectively indicating the number of other cities which are flat in terrain and are subjected to flood damage compared with the current city and the number of other cities which are not subjected to flood damage compared with the terrain of the current city in the set number Y of other cities;
Drawings
Fig. 1 is a flow chart of the urban flood control and drainage risk prediction method of the invention.
Detailed Description
The conception of the invention is as follows:
according to the method, after the historical water storage variation and the water level variation are compared to obtain the water storage capacity corresponding to the unit water level variation of the current city, the flood occurrence probability of the current city at the future time is predicted based on the water storage capacity corresponding to the unit water level variation, and then the flood occurrence probability is corrected by respectively combining three correction parameters obtained by the historical factor, the land factor and the terrain factor, so that the corrected flood occurrence probability capable of more accurately predicting the flood is obtained, the flood judgment is completed, and the accuracy of predicting whether the flood occurs or not is improved.
The following describes a method for predicting risk of flood control and drainage in cities in detail with reference to the accompanying drawings and embodiments.
The method comprises the following steps:
the embodiment of the urban flood control and drainage risk prediction method provided by the invention has the overall flow as shown in figure 1, and the specific process is as follows:
the method comprises the steps of firstly, obtaining the historical daily actual precipitation of the current city, the future daily predicted precipitation, daily drainable quantity, daily average sewage discharge quantity and historical water level of the historical end time of each day.
In order to realize the flood disaster early warning of the current city, the flood control capability condition of the current city needs to be judged by combining historical hydrological information of the current city.
Therefore, the embodiment firstly obtains the actual daily rainfall amount of the current city in history through the data issued by the official departments such as the meteorological office or the hydrological officeAnd also according to the data issued by the official department of the weather bureau, the forecast value of the precipitation of the current city today and every day thereafter is obtained. Wherein, the first and the second end of the pipe are connected with each other,representing the actual precipitation of the current city on day j before today,representing the predicted precipitation of the current city today,representing the predicted precipitation for the day i after the current city today. It is easy to understand that today's actual precipitation is not available for the entire day today, but only for prediction of precipitation today by weather prediction, since it is not completely finished today.
Then, the daily drainable quantity of the current city is obtained according to official dataThe daily drainable amount is a current embodiment of city drainage capacity and can be determined by planning data of the current city. Meanwhile, the sewage discharge amount of the current city in a certain time is obtained from historical data, and the daily average sewage discharge amount of the current city is determined according to the sewage discharge amount in the certain time。
Meanwhile, the historical water level of the current city at the end of each day in history is obtainedWherein, in the step (A),representing the historical water level value at the end of the jth day before the current city today.
And step two, determining the water storage capacity corresponding to the unit water level change of the current city according to the actual precipitation, daily dischargeable water discharge, daily average sewage discharge and historical water level of the current city in a certain historical time.
It can be understood that, the early warning of the flood disaster for different cities is required, not only because the rainfall conditions of different cities are different due to different geographical locations, but also because the rainfall bearing capacities of different cities are different, specifically, the actual water storage amount corresponding to the unit water level change of each city is different, so even if the same rainfall occurs, the situation that one city can have the flood disaster and the other city can not have the flood disaster is likely to occur for different cities.
In consideration of the above situation, the embodiment first determines the water storage capacity corresponding to the current city unit water level change according to the current city precipitation situation in a certain historical time period.
Specifically, the amount of change in the accumulated water per day in history is first determined:
wherein the content of the first and second substances,representing the amount of change in the impounded water for the j day before the current city today,representing the actual precipitation of the current city on day j before today,indicating the daily average sewage discharge amount of the current city,representing the daily drainable capacity of the current city.
Then, selecting the set historical days J before today from the historical data of the current city, thereby determining the water storage capacity corresponding to the unit water level change of the current city:
wherein, the first and the second end of the pipe are connected with each other,the water storage capacity corresponding to the current city unit water level change,indicating the set historical number of days since today,representing the amount of change in the impounded water for the j day before the current city today,the historical water level value representing the current city at the end of yesterday,indicating that the current city is the first from todayHistorical water level values at the end of day time.
In the water storage capacity calculation formula,the method is characterized in that the water storage variation from the J th day to the yesterday before today is shown, and the variation is compared with the variation of the water level value in the whole process from the J +1 th end moment before today, namely from the J th start moment to the yesterday end moment before todayAnd comparing to determine the water storage capacity corresponding to the current city unit water level change.
Therefore, the variable quantity of the water storage capacity corresponding to the whole city after the unit water level change of the current city is determined, and the water storage capacity corresponding to the unit water level change of the current city is also obtained.
And step three, determining the predicted water levels of the current day and the end time of each day in the future according to the water storage capacity corresponding to the current city unit water level change and the predicted precipitation of each day in the future.
The change of the water level is a continuous process related to the water level value at the end time of the previous day, namely the water level value at the end time of a certain day is obtained by the water storage change of the city at the day and the water level value at the end time of the previous day.
Therefore, to achieve the determination of the predicted water level for today and the end of each day in the future, the present embodiment first determines the predicted impounded water variation for the ith day after the present city:
wherein the content of the first and second substances,representing the predicted impounded water variation of the current city on day i after today,representing the predicted precipitation for the current city on day i after today,represents the daily average sewage discharge amount of the current city,representing the daily drainable amount of the current city.
Then determining the predicted water level of the current city at the end time of day I after today:
wherein the content of the first and second substances,indicating the predicted water level at the end of day I after the current city,representing the predicted impounded water change amount of the current city on the ith day after today,representing the historical water level value of the current city at the yesterday end time.
In the predicted water level calculation formula,indicating the predicted amount of change in the impounded water during the course of day to predicted day I after day, soThe predicted variation of the water level during the ith day after today is shown, and the predicted variation is summed with the historical water level value at the end time of yesterday to obtain the predicted water level at the end time of ith day after today.
It is understood that I and I are integers and both values are [0, +,infinity ], when I is 0,indicating that the precipitation is predicted today, when the value of I is 0,indicating the predicted water level at the end of today.
And step four, determining the flood probability at the end time of today and the end time of each day in the future according to the predicted water levels at the end time of today and the end time of each day in the future.
As can be known by combining with the common knowledge in the field, the water level line of three levels is usually set to mark the city water level state, specifically, the normal water level line and the alarmThe three water level lines respectively correspond to the determined height values thereof, and the heights of the normal water level line, the alarm water level line and the dangerous water level line are respectively recorded as the heights of the normal water level line, the alarm water level line and the dangerous water level line in the embodiment、Andthen according to the predicted water level of the end time of the day I after the current city todayAnd finishing the calculation of flood probability at the finishing time of the I day after the current city today according to the size relation among the heights of the three water level lines:
wherein the content of the first and second substances,the flood probability of the current city at the end of the day I after the current day,、andrespectively are the height values of a normal water level line, a warning water level line and a dangerous water level line,predicting a weight for the flood probability, wherein,For the predicted water level at the end of day I after the current city today,this indicates that the value in the parenthesis is 0 when the number in the parenthesis is not positive, and the value in the parenthesis is self value when the number in the parenthesis is positive.
Prediction weight value calculated in flood occurrence probabilityIn the determination process, as the predicted water level gradually approaches the normal water level line, the alarm water level line and the dangerous water level line, the values of the first term, the second term and the third term multiplied in the corresponding calculation process are non-linearly reached to 1 by a method that the speed increase is gradually increased from a value larger than zero to a value smaller than 1, so that the predicted weight value is gradually increased along with the gradual increase of the water levelThe water level of the flood is correspondingly increased, the increasing speed of the water level is gradually increased, and the characteristics that the increasing amount of the flood probability caused by unit water level change is different when the water level is high compared with the water level is low are reflected.
And step five, determining a set date section comprising the reference day by taking the number of days corresponding to the predicted flood occurrence probability as the reference day, determining the water level condition of each historical year on the set date section and the flood occurrence condition, comparing the water level condition with the predicted water level of the end time of the reference day, and determining historical factor correction parameters for correcting the flood occurrence probability.
In this embodiment, the flood occurrence probability at the end time of the I th day after today is obtained, then the set number of days is taken forward and backward respectively on the basis of the I th day after today in this embodiment, the set number of days is 5 days in this embodiment, in other embodiments, the set number of days may also be taken as other values according to the requirement for the prediction accuracy of the flood occurrence probability, and when the requirement for the prediction accuracy is higher, the value of the set number of days may be taken as larger. After taking the set number of days in front and back, a set date segment comprising a certain number of days is obtained, and the set date segment in this embodiment consists of 11 days. In the embodiment, the set date section is obtained by respectively taking the set days forward and backward on the basis of the day I after the day, and in other embodiments, any means for determining the set date section can be adopted, so that the set date section only needs to be ensured to comprise the day I after the day, namely the reference day.
Then, in the set historical number of years, the situation that flood disasters occur in the set date segment in each year is counted, the set historical number of years is 20 years, and other values can be taken according to the prediction accuracy requirement on the flood occurrence probability for the same set historical number of years.
In each historical year, the daily precipitation and the corresponding situation whether flood occurs in the time period corresponding to the set date segment, specifically, in the 11 days in the history which are the same as the 11 days in the set date segment in the whole year, specifically, if the predicted flood occurrence probability at the end time of the day I after today is specifically the probability of No. 9 and No. 16, the set date segment is formed by 11 days from No. 9 and No. 11 to No. 9 and No. 21, and the range of the searched daily precipitation and the corresponding situation whether flood occurs in the historical year is the daily precipitation and the corresponding situation whether flood occurs in 11 days from No. 9 and No. 11 to No. 9 and No. 21 in the historical year.
In the set historical years, namely 20 years, the water level value is counted to be less than the prediction probability in 11 days from No. 11 month 9 to No. 21 month 9 in each yearThe current city corresponding to it is todayPredicted water level at end of the last day IAnd the number of years of flood disasterAnd the water level value is greater than the predicted probabilityThe corresponding predicted water level of the ending time of the I day after the current city todayAnd the number of years of flood disasterThen determine the prediction probabilityThe first historical factor correction parameter and the second historical factor correction parameter for correction:
wherein the content of the first and second substances,andrespectively representing a first historical factor correction parameter and a second historical factor correction parameter,indicating that the number of years of history is set,indicating that the water level value is less than the predicted water level in the set date period in the set historical number of yearsAnd the number of yearly flood disasters occur,indicating that the water level value is greater than the predicted water level in the set date section in the set historical year numberAnd the number of years of flood disasters does not occur.
In setting the historical year, if the water level value is less than the predicted water levelAnd the number of years of flood disasterThe more, the more the flood probability based on the predicted water level is indicatedThe smaller the size is; if the water level value is greater than the predicted water levelAnd the number of years of flood disasterThe more, the more the flood probability based on the predicted water level is indicatedThe larger the size of the system is, the probability of flood can be correspondingly obtainedThe first and second historical factor correction parameters are corrected.
And step six, determining a land factor correction parameter for correcting the flood disaster occurrence probability according to the occupation ratio of the soil area in the current city in the total area of the city.
The daily drainable amount of the city is not only influenced by the delivery performance of the drainage network of the city for delivering precipitation outside the city, but also is substantially influenced by the absorption and storage performance of the land area in the city, and in the daily drainable amount of the city counted by the official, the daily drainable amount of the city actually implies the drainage amount corresponding to the absorption and storage performance of the land area in the city for the precipitation.
However, the water absorption of the soil is limited, and the water level in the city is generated along with the long-term occurrence of rainfall, and actually, the absorption and storage performance of the soil for the rainfall reaches the upper limit at the moment, and the rainfall is accumulated above the ground just because the soil in the city cannot absorb the rainfall any more.
Therefore, in the case of the water level, the actual daily dischargeable quantity of the city is smaller than the daily dischargeable quantity given by the authority, and the larger the area of the urban land is, the more obvious the actual daily dischargeable quantity of the city is reduced compared with the daily dischargeable quantity given by the authority in the case of the water level, and the more the predicted flood occurrence probability should be corrected to a greater extent in the flood occurrence probability prediction process.
In the embodiment, the occupation ratio Q of the urban soil area in the total area is counted, and the occupation ratio Q is used as a land factor correction parameter for correcting the flood occurrence probability.
And step seven, determining a relief factor correction parameter for correcting the flood occurrence probability according to the relief difference between the current city and other cities and the flood occurrence conditions of other cities.
Whether flood disasters occur or not is also influenced by the situation difference from the periphery of the city to the city, and the lower the city is compared with the periphery, the higher the probability of flood disasters is, and the lower the probability is.
Therefore, in the present embodiment, the city area is taken as the center, the periphery of the city is divided into the areas with the set number of areas, and the number of the areas is preferably set to 8, so that the periphery of the city is correspondingly divided into 8 areas, the altitude of each peripheral area and the altitude of the city area are obtained, the altitude of the city area is sequentially differed from the altitude of each peripheral area, the 8 altitude differences are obtained, normalization processing is performed, and the terrain sequence of the current city is obtained by sorting the differences from small to large。
Then, the water level value in the set date section and the predicted water level of the current city at the end time of the day I after the current city are selected from other citiesThe difference value of the current city is smaller than that of the preset water level difference value, and the terrain sequence of the other selected cities is obtained according to the method for obtaining the terrain sequence of the current cityAnd labeling the selected terrain sequences of other cities, wherein the label is V, when the value of V is 0, the fact that the flood disaster does not occur is shown, and when the value of V is 1, the fact that the flood disaster occurs is shown.
The terrain sequence of the current city is multiplied by itself to obtain a terrain matrix of the current cityDetermining the maximum and minimum values of the elements in the matrixAndthen will be composed ofAndthe value ranges of the elements in the determined matrix are divided into set number of parts, 8 parts are preferred in the embodiment, and the value of the set number of parts is determined according to the correction accuracy requirement. After the value ranges of the elements in the matrix are divided into the set number of parts, all the elements in the matrix are correspondingly in the set number of parts, if the set number of parts is 8, which is preferred in this embodiment, the number of the elements in the matrix, whose values fall into the second half of the value ranges, that is, the number of the elements whose values fall into the larger of the two front and rear half value ranges, can be obtained through statistics.
Similarly, the terrain sequences of other cities are multiplied to obtain the terrain matrix of other citiesAnd determining the maximum value and the minimum value of elements in the terrain matrix of other cities according to the same methodAndand correspondingly dividing all elements in the terrain matrixes of other cities into the set parts.
Then, calculating the terrain superiority and inferiority parameters of the current city compared with other cities:
wherein the content of the first and second substances,representing the terrain superiority and inferiority parameters of the current city compared with other cities,and withRespectively representing the maximum value and the minimum value of the elements in the terrain matrix of the current city,andrespectively representing the maximum and minimum values of the elements in the terrain matrix of the other cities,the number of elements whose values fall into the latter half of the value range of the elements in the terrain matrix of the current city,the number of elements whose values fall into the latter half of the value range of the elements in the terrain matrix of other cities is represented.
Terrain quality parameterThe size of (b) represents the magnitude of the terrain difference degree of the current city compared with other cities,andrespectively showing the ratio of the maximum value of the difference between the current city and the surrounding terrain to the maximum value of the difference between the other cities and the surrounding terrain, and the ratio of the minimum value of the difference between the current city and the surrounding terrain to the minimum value of the difference between the other cities and the surrounding terrain, wherein the larger the two ratios are, the more obvious the difference between the current city and the surrounding terrain is compared with the difference between the other cities and the surrounding terrain,the difference between the current city and the surrounding areas is larger than that between other cities, and finally, the terrain superiority and inferiority parametersThe larger the size, the more depressed the city is compared to other cities.
When in useIf the number of the urban areas is more than 1, the situation shows that the current city is more prone to flood disasters compared with other cities, and when the number of the urban areas is more than 1, the situation is that the current city is more prone to flood disastersAnd if the number of the urban areas is less than 1, the current city is less prone to flood disasters compared with other cities. In this embodiment, the number of other cities is set as the number Y of other cities, and the values of the number Y of other cities are set according to the accuracy requirement of flood prediction, and obviously, the accuracy requirement is correspondingly higher if the values of the number Y of other cities are set to be larger.
Counting each corresponding to the number Y of other citiesThe number of other cities which are larger than 1 and corresponding to other cities in flood disasterAnd anThe number of other cities which are less than 1 and are not flood-damaged is corresponding to other cities. To be provided withAndas a relief factor correction parameter for correcting the flood occurrence probability.
As can be readily appreciated, inIn the corresponding situation, because other cities have flood disasters, and the current city is in a low-lying area compared with other cities, the flood disasters of the current city are more likely to occur; in addition, theIn the corresponding situation, since no flood occurs in other cities, and the current city has a lower depression degree, i.e., a flatter topography compared to other cities, the probability that no flood occurs in the current city is higher.
And step eight, correcting the flood occurrence probability according to the determined historical factor correction parameters, the land factor correction parameters and the relief factor correction parameters, and finishing flood prediction according to the corrected flood occurrence probability.
According to the three correction factors, the embodiment corrects the obtained flood occurrence probability value to obtain the corrected flood occurrence probability:
wherein, the first and the second end of the pipe are connected with each other,indicating the flood probability of the revised current city at the end time of the I day after the current day,andrespectively representing a first historical factor correction parameter and a second historical factor correction parameter,the land factor correction parameter is represented by a land factor correction parameter,andin the other-city number Y, the number of other cities in which the terrain is flatter than the current city and flood occurs and the number of other cities in which the terrain is low than the current city and flood does not occur are indicated.
Probability of flood after correctionIf the threshold value is larger than the threshold value of the flood occurrence probability, the flood disaster can occur, otherwise, the flood disaster does not occur. In this embodiment, the threshold of the flood occurrence probability is preferably 0.8, and in other embodiments, the threshold of the flood occurrence probability may be set to other values according to the requirement of the flood prevention sensitivity level.
The above-mentioned embodiments are only used for illustrating the technical solutions of the present application, and not for limiting the same; although the present application has been described in detail with reference to the foregoing embodiments, it should be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; such modifications or substitutions do not cause the essential features of the corresponding technical solutions to depart from the scope of the technical solutions of the embodiments of the present application, and are intended to be included within the scope of the present application.
Claims (8)
1. A city flood control and drainage risk prediction method is characterized by comprising the following steps:
determining the amount of variation of water storage of the current city in the current set historical days today according to the actual precipitation, daily average sewage discharge and daily drainable quantity of the current city in the history of the current city, and determining the amount of variation of water level of the current city in the current set historical days today according to the historical water level of the current city in the end moment of each day in the history of the current city;
determining the water storage capacity corresponding to the current city unit water level change according to the water storage change and the water level change;
determining the predicted water level of the current city every day after the current day according to the water storage capacity corresponding to the unit water level change of the current city, the predicted precipitation amount of the future every day, the daily drainable amount and the daily average sewage discharge amount, thereby determining the flood occurrence probability at the end of each day after the current day;
determining a set date section comprising the reference day by taking the number of days corresponding to the predicted flood occurrence probability as the reference day, and determining the flood disaster occurrence condition and the water level condition in the set date section in each year in the set historical number of years to obtain historical factor correction parameters for correcting the flood occurrence probability;
determining a land factor correction parameter for correcting the flood disaster occurrence probability according to the occupation ratio of the soil area in the current city in the total area of the city;
determining a relief factor correction parameter for correcting the flood occurrence probability according to the relief difference between the current city and other cities and the flood occurrence conditions of other cities;
and correcting the flood occurrence probability according to the determined historical factor correction parameters, the land factor correction parameters and the relief factor correction parameters, and completing flood prediction according to the corrected flood occurrence probability.
2. The urban flood control and flood drainage risk prediction method according to claim 1, wherein the method for determining the water storage capacity corresponding to the current urban unit water level change comprises the following steps:
wherein the content of the first and second substances,the water storage capacity corresponding to the current city unit water level change,indicating the set historical number of days since today,representing the amount of change in the impounded water for the j day before the current city today,the historical water level value representing the current city at the end of yesterday,indicating that the current city is the first from todayThe historical water level value at the end of the day,representing the actual precipitation of the current city on day j before today,represents the daily average sewage discharge amount of the current city,representing the daily drainable capacity of the current city.
3. The method for predicting the risk of flood control and drainage in cities according to claim 1, wherein the predicted water level of the current city after today every day is:
wherein the content of the first and second substances,indicating the predicted water level at the end of day I after the current city today,representing the predicted impounded water variation of the current city on day i after today,the historical water level value representing the current city at yesterday's end time,representing the predicted precipitation for the current city on day i after this day,represents the daily average sewage discharge amount of the current city,representing the daily drainable capacity of the current city.
4. The method for predicting the urban flood control and flood drainage risk according to claim 1, wherein the flood occurrence probability at the end of each day after the current day is:
wherein the content of the first and second substances,the flood probability of the current city at the end of the day I after the current day,、andrespectively are the height values of a normal water level line, a warning water level line and a dangerous water level line,predicting a weight for the flood probability, wherein,For the predicted water level at the end of day I after the current city today,the numerical value in parentheses is 0 when the numerical value in parentheses is not positive, and the numerical value in parentheses is self-value when the numerical value in parentheses is positive.
5. The urban flood control and drainage risk prediction method according to claim 1, wherein the method for obtaining the historical factor correction parameters comprises:
wherein, the first and the second end of the pipe are connected with each other,andrespectively representing a first historical factor correction parameter and a second historical factor correction parameter,indicating that the number of the historical year copies is set,indicating that the water level value is less than the predicted water level in the set date period in the set historical number of yearsAnd the number of yearly flood disasters occur,indicating that the water level value is greater than the predicted water level in the set date period in the set historical number of yearsAnd the number of years of flood disasters does not occur.
6. The urban flood control and drainage risk prediction method according to claim 1, wherein the method for obtaining the land factor correction parameter comprises:
and (4) counting the proportion Q of the total soil area in the current city, and taking the proportion Q as a land factor correction parameter.
7. The urban flood control and drainage risk prediction method according to claim 1, wherein the method for obtaining the terrain factor correction parameters comprises:
dividing the periphery of the city into peripheral areas with the number of the set areas, sequentially making a difference between the altitude of the city area and the altitude of each peripheral area to obtain the altitude difference values of the number of the set areas, then performing normalization processing, and sequencing according to the sequence from small to large to obtain a terrain sequence;
determining the terrain sequence of the current city and the terrain sequences of other cities according to a terrain sequence acquisition method, and then calculating the terrain advantages and disadvantages parameters of the current city compared with other cities:
wherein, the first and the second end of the pipe are connected with each other,representing the terrain superiority and inferiority parameters of the current city compared with other cities,andrespectively representing the maximum value and the minimum value of elements in the terrain matrix of the current city,andrespectively representing the maximum and minimum values of the elements in the terrain matrix of the other cities,representing the number of elements whose values fall into the latter half of the value range of the elements in the terrain matrix of the current city,representing the number of elements of which the values fall into the second half of the value range of the elements in the terrain matrix of other cities;
selecting other cities with the set number of Y, and counting the number of Y to obtain each cityThe number of other cities which are larger than 1 and corresponding to other cities in flood disasterAnd anThe number of the value is less than 1 and corresponds to the number of other cities without flood disastersTo do so byAndas a relief factor correction parameter for correcting the flood occurrence probability.
8. The method for predicting the urban flood control and flood drainage risk according to claim 1, wherein the method for completing flood prediction according to the corrected flood occurrence probability comprises the following steps:
firstly determining the corrected flood occurrence probability:
wherein the content of the first and second substances,indicating the flood probability of the revised current city at the end time of the day I after the current day,and withRespectively representing a first historical factor correction parameter and a second historical factor correction parameter,a land factor-based correction parameter is represented,andrespectively indicating the number of other cities which are flat in terrain and are subjected to flood damage compared with the current city and the number of other cities which are not subjected to flood damage compared with the terrain of the current city in the set number Y of other cities;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211231275.7A CN115293469B (en) | 2022-10-10 | 2022-10-10 | Urban flood control and drainage risk prediction method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211231275.7A CN115293469B (en) | 2022-10-10 | 2022-10-10 | Urban flood control and drainage risk prediction method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115293469A true CN115293469A (en) | 2022-11-04 |
CN115293469B CN115293469B (en) | 2022-12-27 |
Family
ID=83819226
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211231275.7A Active CN115293469B (en) | 2022-10-10 | 2022-10-10 | Urban flood control and drainage risk prediction method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115293469B (en) |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101441683A (en) * | 2008-06-17 | 2009-05-27 | 浙江大学 | Prediction method of city rainstorm flood evolvement process |
CN104298841A (en) * | 2013-07-16 | 2015-01-21 | 杭州贵仁科技有限公司 | Flood forecasting method and system based on historical data |
CN106779232A (en) * | 2016-12-28 | 2017-05-31 | 合肥学院 | A kind of modeling and forecasting method of urban waterlogging |
CN109146140A (en) * | 2018-07-25 | 2019-01-04 | 佛山科学技术学院 | A kind of Urban Flood Waterlogging prediction technique |
CN110837925A (en) * | 2019-11-04 | 2020-02-25 | 浙江贵仁信息科技股份有限公司 | Urban waterlogging prediction method and device |
CN111613017A (en) * | 2020-04-26 | 2020-09-01 | 李海波 | Urban flood disaster emergency early warning processing method and system |
CN111814407A (en) * | 2020-07-28 | 2020-10-23 | 安徽沃特水务科技有限公司 | Flood forecasting method based on big data and deep learning |
CN112069378A (en) * | 2020-09-14 | 2020-12-11 | 郑州大学 | Urban disaster-causing rainfall classification determination method based on IDF curve |
CN112785053A (en) * | 2021-01-15 | 2021-05-11 | 北京市水科学技术研究院 | Method and system for forecasting urban drainage basin flood |
US20210192627A1 (en) * | 2020-02-27 | 2021-06-24 | Intercontinental Exchange Holdings, Inc. | Systems and methods for converting live weather data to weather index for offsetting weather risk |
CN113554236A (en) * | 2021-07-28 | 2021-10-26 | 武汉中科集思地理信息科技有限公司 | Smart city construction flood disaster real-time online monitoring and early warning management cloud platform |
CN113610264A (en) * | 2021-06-22 | 2021-11-05 | 国网浙江省电力有限公司电力科学研究院 | Refined power grid typhoon flood disaster prediction model |
CN113762645A (en) * | 2021-10-11 | 2021-12-07 | 昆仑(重庆)河湖生态研究院(有限合伙) | Natural disaster forecasting method and device |
CN115019477A (en) * | 2022-07-18 | 2022-09-06 | 广东广宇科技发展有限公司 | Intelligent urban waterlogging disaster real-time prediction early warning system and method |
-
2022
- 2022-10-10 CN CN202211231275.7A patent/CN115293469B/en active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101441683A (en) * | 2008-06-17 | 2009-05-27 | 浙江大学 | Prediction method of city rainstorm flood evolvement process |
CN104298841A (en) * | 2013-07-16 | 2015-01-21 | 杭州贵仁科技有限公司 | Flood forecasting method and system based on historical data |
CN106779232A (en) * | 2016-12-28 | 2017-05-31 | 合肥学院 | A kind of modeling and forecasting method of urban waterlogging |
CN109146140A (en) * | 2018-07-25 | 2019-01-04 | 佛山科学技术学院 | A kind of Urban Flood Waterlogging prediction technique |
CN110837925A (en) * | 2019-11-04 | 2020-02-25 | 浙江贵仁信息科技股份有限公司 | Urban waterlogging prediction method and device |
US20210192627A1 (en) * | 2020-02-27 | 2021-06-24 | Intercontinental Exchange Holdings, Inc. | Systems and methods for converting live weather data to weather index for offsetting weather risk |
CN111613017A (en) * | 2020-04-26 | 2020-09-01 | 李海波 | Urban flood disaster emergency early warning processing method and system |
CN111814407A (en) * | 2020-07-28 | 2020-10-23 | 安徽沃特水务科技有限公司 | Flood forecasting method based on big data and deep learning |
CN112069378A (en) * | 2020-09-14 | 2020-12-11 | 郑州大学 | Urban disaster-causing rainfall classification determination method based on IDF curve |
CN112785053A (en) * | 2021-01-15 | 2021-05-11 | 北京市水科学技术研究院 | Method and system for forecasting urban drainage basin flood |
CN113610264A (en) * | 2021-06-22 | 2021-11-05 | 国网浙江省电力有限公司电力科学研究院 | Refined power grid typhoon flood disaster prediction model |
CN113554236A (en) * | 2021-07-28 | 2021-10-26 | 武汉中科集思地理信息科技有限公司 | Smart city construction flood disaster real-time online monitoring and early warning management cloud platform |
CN113762645A (en) * | 2021-10-11 | 2021-12-07 | 昆仑(重庆)河湖生态研究院(有限合伙) | Natural disaster forecasting method and device |
CN115019477A (en) * | 2022-07-18 | 2022-09-06 | 广东广宇科技发展有限公司 | Intelligent urban waterlogging disaster real-time prediction early warning system and method |
Non-Patent Citations (1)
Title |
---|
吕苹等: "昆明城区防洪减灾问题探讨", 《人民长江》 * |
Also Published As
Publication number | Publication date |
---|---|
CN115293469B (en) | 2022-12-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Tingsanchali et al. | Flood hazard and risk analysis in the southwest region of Bangladesh | |
CN112506994B (en) | Power equipment flood hidden danger point monitoring and early warning method and related device | |
CN110991046B (en) | Drainage system waterlogging risk rapid early warning method based on response surface function | |
Afzali-Gorouh et al. | Probable maximum precipitation estimation in a humid climate | |
CN114360214B (en) | Extra-large scale glacier debris flow early warning method | |
CN113762645B (en) | Natural disaster forecasting method and device | |
Arora et al. | Spatial distribution and seasonal variability of rainfall in a mountainous basin in the Himalayan region | |
CN111815043B (en) | Flood flow prediction method and system based on storm characteristics | |
CN110619433B (en) | Rapid selection method and system for power grid heavy rain numerical mode parameterization scheme | |
CN115293469B (en) | Urban flood control and drainage risk prediction method | |
CN115619045A (en) | Estuary tidal range forecasting method and system | |
CN107944627A (en) | A kind of river basin flood forecasting procedure | |
CN115203639A (en) | Irregular grid surface rainfall calculation method and system based on matrix operation | |
CN115994602A (en) | Dynamic control method and dynamic control equipment for water level in reservoir flood period | |
CN114840989A (en) | River confluence calculation method with grid scale considering hydraulic engineering storage | |
CN114037248A (en) | River ecological risk assessment method based on ecological deficit index | |
Bacalhau et al. | Water supply reservoir operation in relation to climate variability: Pirapama river basin (Pernambuco-Brazil) | |
CN113589404A (en) | Method for predicting runoff of storm of field | |
Thoảng et al. | Assessment of Potential Rainfall Distribution Patterns and Their Relationship with Inundation in Tra Vinh Province, Vietnam | |
CN111768310A (en) | Reservoir water replenishing potential prediction method and device and electronic equipment | |
CN116805031B (en) | Prediction method and device for flood inundation risk | |
CN113642261B (en) | Ionized layer TEC disturbance determination method and system | |
CN114548559B (en) | Mountain torrent disaster dynamic early warning method based on early-stage rainfall | |
CN117436619B (en) | Cascade reservoir flood control reservoir capacity combined reservation method based on equivalent flood control effect | |
CN116679357B (en) | Method for carrying out water balance treatment on radar quantitative precipitation estimation |
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 |