CN115270525B - Method and system for constructing dual EXP function runoff coefficient design model - Google Patents
Method and system for constructing dual EXP function runoff coefficient design model Download PDFInfo
- Publication number
- CN115270525B CN115270525B CN202211179274.2A CN202211179274A CN115270525B CN 115270525 B CN115270525 B CN 115270525B CN 202211179274 A CN202211179274 A CN 202211179274A CN 115270525 B CN115270525 B CN 115270525B
- Authority
- CN
- China
- Prior art keywords
- runoff
- model
- node
- rainfall
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
-
- 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)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Sewage (AREA)
Abstract
The method and the system for constructing the dual EXP function runoff coefficient design model provided by the embodiment of the application comprise the steps of obtaining an accumulated rainfall index and a regional average population density index; grading the accumulated rainfall index to obtain graded accumulated rainfall; grading the area average population density index to obtain the grading area average population density; constructing a dual EXP function runoff coefficient model, combining the graded accumulated rainfall and the graded region average population density, and determining a model independent variable based on the obtained multiple node runoff calculation indexes; determining instantaneous comprehensive coefficients of runoff calculation indexes of each node based on a double EXP function runoff coefficient model, taking weighted mean values of the instantaneous comprehensive coefficients as node mean value samples, and screening out key node mean value samples from the node mean value samples according to the feasibility of national standard control; and when the weighted mean value of the target instantaneous comprehensive runoff coefficient obtained based on the key node mean value sample accords with the national standard control value, constructing a target model based on the currently adjusted model parameters.
Description
Technical Field
The application relates to the technical field of water supply and drainage engineering, in particular to a method and a system for constructing a dual EXP function runoff coefficient design model.
Background
With the rapid development of urbanization in China, the service range of urban rainwater engineering is continuously expanded, and the risk of inland inundation and disaster prevention is obviously increased. The runoff coefficient value for determining the flow of the rainwater runoff has important influence on the disaster bearing capacity of an urban waterlogging prevention system and a rainwater project. At present, the runoff coefficient method proposed by national standards is mostly adopted in domestic rainwater engineering to determine the designed rainwater runoff of a pipe duct, wherein a comprehensive runoff coefficient value taking method is mostly adopted according to the urban building density grade. However, the elasticity range of the comprehensive runoff coefficient value interval corresponding to the building density of different towns reaches 15% -80%, and the national standard does not further clearly define the detailed value. Therefore, when the urban construction density and the value of the corresponding comprehensive runoff coefficient value interval are judged, the existing artificial influence factors are large, and the quantitative design level needs to be improved.
Disclosure of Invention
The embodiment of the application aims to provide a method and a system for constructing a dual EXP function runoff coefficient design model, so that the reasonability of rainfall runoff calculation in the whole rainfall process can be improved.
The embodiment of the application further provides a method for constructing a dual EXP function runoff coefficient design model, which comprises the following steps:
s1, obtaining a plurality of node runoff calculation indexes, wherein the node runoff calculation indexes comprise accumulated rainfall indexes and regional average population density indexes;
s2, grading the accumulated rainfall index according to the value of the rainfall recurrence period to obtain graded accumulated rainfall;
s3, according to the maximum range of the average population density of the region, carrying out grade division on the average population density index of the region to obtain the average population density of the classified region;
s4, constructing a double-EXP function runoff coefficient model based on a double-EXP function form, combining graded accumulated rainfall and graded region average population density, and determining a model independent variable based on the obtained multiple node runoff calculation indexes;
s5, determining instantaneous comprehensive runoff coefficients of runoff calculation indexes of all nodes based on the double EXP function runoff coefficient model, taking weighted average values of the instantaneous comprehensive runoff coefficients as node average value samples, and screening key node average value samples from the node average value samples according to feasibility of national standard control;
s6, when the weighted mean value of the target instantaneous comprehensive runoff coefficient correspondingly obtained based on the key node mean value sample is determined to accord with the national standard control value, constructing the target model based on the model parameter in the current adjusting stage, otherwise, entering the next parameter adjusting stage, and continuously approaching the weighted mean value of the target instantaneous comprehensive runoff coefficient to the national standard control value by adjusting the model parameter.
In a second aspect, an embodiment of the present application further provides a system for constructing a dual EXP function runoff coefficient design model, where the system includes an index obtaining module, a ranking module, and a target model constructing module, where:
the index acquisition module is used for acquiring a plurality of node runoff calculation indexes, and the node runoff calculation indexes comprise accumulated rainfall indexes and regional average population density indexes;
the grading module is used for grading the accumulated rainfall index according to the value of the rainfall recurrence period to obtain graded accumulated rainfall;
the grading module is further configured to grade the area average population density index according to the area average population density maximum range to obtain a graded area average population density;
the target model building module is used for building a double-EXP function runoff coefficient model based on a double-EXP function form, combining the graded accumulated rainfall and the graded region average population density, and determining a model independent variable based on the obtained multiple node runoff calculation indexes;
the target model building module is further used for determining an instantaneous comprehensive coefficient of each node runoff calculation index based on the double EXP function runoff coefficient model, taking a weighted average of the instantaneous comprehensive coefficients as a node average sample, and screening key node average samples from the node average sample according to feasibility of national standard control;
the target model building module is further used for building a target model based on the model parameters in the current adjusting stage when the weighted mean value of the target instantaneous comprehensive runoff coefficient correspondingly obtained based on the key node mean value sample is determined to accord with the national standard control value, otherwise, entering the next parameter adjusting stage, and enabling the weighted mean value of the target instantaneous comprehensive runoff coefficient to continuously approach the national standard control value by adjusting the model parameters.
In a third aspect, an embodiment of the present application further provides a readable storage medium, where the readable storage medium includes a program of a method for constructing a dual EXP function runoff coefficient design model, and when the program of the method for constructing a dual EXP function runoff coefficient design model is executed by a processor, the steps of the method for constructing a dual EXP function runoff coefficient design model as described in any one of the above are implemented.
From the above, the method, the system and the readable storage medium for constructing the dual EXP function runoff coefficient design model provided by the embodiment of the application use the national standard control value (i.e. the national standard comprehensive runoff coefficient control index) as a basis, adopt the regression analysis method, and establish the runoff coefficient design model based on the dual EXP function form and the runoff calculation index with large influence and strong operability, so that under the condition of better connection with the national standard, the traditional fixed value of the comprehensive runoff coefficient can be converted into the instantaneous dynamic value, the rationality of rainfall whole-course rainfall runoff calculation is improved, the design system for determining the rainwater flow by using the mathematical model method is facilitated and perfected, and the design scientificity of the urban drainage waterlogging prevention engineering is improved.
Additional features and advantages of the present application will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by the practice of the embodiments of the present application. The objectives and other advantages of the application may be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings that are required to be used in the embodiments of the present application will be briefly described below, it should be understood that the following drawings only illustrate some embodiments of the present application and therefore should not be considered as limiting the scope, and that those skilled in the art can also obtain other related drawings based on the drawings without inventive efforts.
Fig. 1 is a flowchart of a method for constructing a dual EXP function runoff coefficient design model according to an embodiment of the present disclosure;
FIG. 2 is a graph comparing two functions (i.e., single and double EXP functions) of the integrated runoff coefficient;
FIG. 3 is a schematic diagram of model mean and node mean sample fitting (time step 5 min);
fig. 4 is a schematic structural diagram of a system for constructing a dual EXP function runoff coefficient design model according to an embodiment of the present application.
Detailed Description
The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are only a part of the embodiments of the present application, and not all of the embodiments. The components of the embodiments of the present application, as generally described and illustrated in the figures herein, could be arranged and designed in a wide variety of different configurations. Thus, the following detailed description of the embodiments of the present application, presented in the accompanying drawings, is not intended to limit the scope of the claimed application, but is merely representative of selected embodiments of the application. All other embodiments, which can be derived by a person skilled in the art from the embodiments of the present application without making any creative effort, shall fall within the protection scope of the present application.
It should be noted that: like reference numbers and letters refer to like items in the following figures, and thus, once an item is defined in one figure, it need not be further defined and explained in subsequent figures. Meanwhile, in the description of the present application, the terms "first", "second", and the like are used only for distinguishing the description, and are not to be construed as indicating or implying relative importance.
Referring to fig. 1, fig. 1 is a flowchart of a method for constructing a dual EXP function runoff coefficient design model according to some embodiments of the present disclosure. The method is exemplified by being applied to a computer device (the computer device may specifically be a terminal or a server, and the terminal may specifically be but is not limited to various personal computers, notebook computers, smart phones, tablet computers and portable wearable devices, the server may be an independent server or a server cluster composed of a plurality of servers), and the method includes the following steps:
step S1, obtaining a plurality of node runoff calculation indexes, wherein the node runoff calculation indexes comprise accumulated rainfall indexes and regional average population density indexes.
And S2, grading the accumulated rainfall index according to the value of the rainfall recurrence period to obtain graded accumulated rainfall.
And S3, according to the maximum range of the average population density of the region, carrying out grade division on the average population density index of the region to obtain the average population density of the classified region.
And S4, constructing a double EXP function runoff coefficient model based on a double EXP function form, and combining the graded accumulated rainfall and the graded region average population density to determine a model independent variable based on the obtained multiple node runoff calculation indexes.
And S5, determining instantaneous comprehensive runoff coefficients of all node runoff calculation indexes based on the double EXP function runoff coefficient model, taking the weighted average of the instantaneous comprehensive runoff coefficients as a node average sample, and screening key node average samples from the node average sample according to the feasibility of national standard control.
And S6, when the weighted average of the target instantaneous comprehensive runoff coefficient correspondingly obtained based on the key node average sample is determined to accord with the national standard control value, constructing a target model based on the model parameters in the current adjusting stage, otherwise, entering the next parameter adjusting stage, and continuously approaching the weighted average of the target instantaneous comprehensive runoff coefficient to the national standard control value by adjusting the model parameters.
Therefore, the method for constructing the dual EXP function runoff coefficient design model disclosed by the application is based on a national standard control value (namely a national standard comprehensive runoff coefficient control index, wherein the national standard adopts GB 50318-2017 urban drainage engineering planning standard and GB 50014-2021 outdoor drainage design standard), adopts a regression analysis method, and establishes the runoff coefficient design model based on the dual EXP function form and the runoff calculation index with large influence and strong operability, so that under the condition of better connection with the national standard, the traditional fixed value of the comprehensive runoff coefficient can be converted into the instantaneous dynamic value, the reasonability of rainfall whole-journey rainwater runoff calculation is improved, the design system for determining the rainwater flow by using a mathematical model method is facilitated and perfected, and the scientificity drainage waterlogging prevention engineering design is improved.
In one embodiment, in step S2, the step of performing level division on the accumulated rainfall index according to a value of a rainfall recurrence period to obtain a graded accumulated rainfall includes:
and S21, based on the value specification of the drainage waterlogging prevention design recurrence period of the urban area based on the outdoor drainage design standard, grading the accumulated rainfall index according to the value of the rainfall recurrence period to obtain a plurality of graded rainfall indexes, wherein the rainfall is designed according to the rainstorm of the preset time duration of the area, and the accumulated rainfall of each rainfall recurrence period can be determined.
In one embodiment, in step S3, the step of performing rank division on the index of average population density according to the maximum range of average population density of a region to obtain the graded average population density of the region includes:
and S31, carrying out grade division on the area average population density index based on a maximum interval range obtained by average population density in a research area to obtain a plurality of graded population density indexes, wherein the minimum value in the maximum interval range is used as a starting unit, and each grade is gradually increased according to a preset length unit until the maximum value in the maximum interval range is reached.
It should be noted that:
(1) The definition of runoff coefficient includes: according to the hydrology principle, the runoff coefficient is the ratio of the runoff depth (or the total runoff) in any time period to the precipitation depth (or the total precipitation) in the same time period, and the proportion of the precipitation quantity converted into the runoff quantity is expressed.
The runoff coefficient is an important technical index commonly used by rain flood control and utilization theoretical analysis and engineering design, and has great influence on the design scale of river channels, pipe ducts and pump stations.
(2) Factors that influence the change in runoff coefficient include: the runoff is greatly influenced by natural environment and construction environment. Therefore, the runoff coefficient is mainly influenced by factors such as rainfall, ground type composition, soil classification, wherein:
(a) The runoff coefficient is influenced by rainfall, mainly because in the early stage of rainfall, the water content in soil is relatively low, the soil infiltration capacity is strong, the runoff production is less, and simultaneously, in the early stage of rainfall, part of rainwater is consumed in plant interception and depression, so that the actual runoff coefficient is less, and as the rainfall continues, the soil water content is increased, the rainwater infiltration rate is reduced, the actual runoff coefficient is gradually increased, and finally, the actual runoff coefficient tends to be stable. Under the influence of rainfall variables, the runoff coefficient in the rainfall process is at a variable value.
(b) The runoff coefficient is influenced by the types of the ground, mainly because the area proportion of urban houses, roads and squares is large, the permeability of the ground is poor, the runoff coefficient is high, and if the area proportion of green lands and non-paved road surfaces of parks is large, the permeability of the ground is good, the runoff coefficient is small.
(c) The runoff coefficient is affected by the soil class, mainly because different soil classes will affect the soil permeability. Wherein, the clay and silt particles with smaller water permeability have larger runoff coefficient, and the gravel and sand particles with larger water permeability have smaller runoff coefficient.
In the present embodiment, the computer device takes the rainfall amount, the ground water permeability as independent variables of the model. The rainfall is designed to be the accumulated rainfall from the rainfall starting point to the calculation node, and the rainfall also represents the rainfall intensity and the rainfall time of the total rainfall duration. The ground water permeability can be combined with the urban construction environment to adopt the partition calculation indexes of average greening rate, average hardening rate and the like.
It should be noted that the urban public greening rate, urban construction hardening rate and other diameter flow calculation indexes can better reflect rainwater infiltration conditions determined by urban ground type composition, but the existing urban information platform is not perfect urban construction statistical indexes and cannot be used. Therefore, in the current embodiment, based on that the average population density of each area in the city can reflect the building density and greening level condition of the city area, the area average population density is adopted to replace the greening rate of the city as another runoff calculation index.
Specifically, in order to meet the comparability of the model construction sample and the national standard under corresponding conditions, the computer equipment correspondingly marks different grades for two runoff calculation indexes, namely the determined comprehensive rainfall index and the regional average population density index, according to the urban building intensity and rainfall grading conditions defined by the comprehensive runoff coefficient values specified by the national standard, so that the calculation result is more accurate.
In one embodiment, the obtained graded rainfall index is divided into 8 graded accumulated rainfall amounts according to the rainfall recurrence period of 2 years, 3 years, 5 years, 10 years, 20 years, 30 years, 50 years and 100 years, and the maximum range of the average population density of the graded population density index in each subarea area of a city is 50 persons/hm 2 496 people/hm 2 In order to ensure that the value sample can adapt to the interval range, the index is started and graded by 50 units, each grade is increased by 50 units to 500 units, and when 10 grading area average population densities are obtained through division, after the 8 grading accumulated rainfall and the 10 grading area average population densities are combined with each other, 80 node runoff calculation indexes are formed.
In one embodiment, in step S4, the node runoff calculation index is expanded as a function in a form of a three-parameter power function, and in step S3, the calculation formula of the dual EXP function runoff coefficient model includes:
wherein the content of the first and second substances,
in order to obtain the instantaneous comprehensive runoff coefficient,Pis the average population density in people in the classification areahm 2 ；HThe rainfall is accumulated in grades, and the unit is mm,
、
、
for the metering parameters covered in the resulting three-parameter power function developed through the rank region average population density,
、
、
the three-parameter power function has strong adaptability to regression line type change for the metering parameters covered in the three-parameter power function obtained by the hierarchical accumulated rainfall expansion.
Specifically, the runoff coefficient is the ratio of the runoff depth (or the runoff total amount) in any time period to the precipitation depth (or the precipitation total amount) in the same time period, so that the value is in a 0-1 change interval, and the condition that the value of the dependent variable of the model always meets the control condition of being more than or equal to 0 and less than or equal to 1 no matter what the value of the independent variable in the model is.
In the current embodiment, the computer device analyzes the numerical characteristic of the comprehensive runoff coefficient, and further reflects the bidirectional change characteristic that the curve change is increased from zero to zero along with the increase of the accumulated rainfall and obeys the increase of the curve particle derivative. Wherein:
(1) The expression for the single EXP function includes:
(2) The expression of the dual EXP function includes:
wherein the content of the first and second substances,
the instantaneous comprehensive runoff coefficient;Pthe average population density of the centralized construction area of the town is expressed in units of peoplehm 2 ；HThe unit of the accumulated rainfall is mm for designing the rainfall duration of the rain type.
It should be noted that according to the increase attribute of the runoff coefficient value from 0 to 1, the function curves (as shown in fig. 2) of the above two formulas express the change characteristics that the runoff coefficient may appear, and the dual EXP function model can reflect the duality of the comprehensive runoff coefficient. The comprehensive runoff coefficient obeys the characteristic of slow increase, namely double EXP function curves, in the early stage of rainfall instead of the characteristic of quick increase of a single EXP function curve, and the double EXP functions adopt accumulated rainfall and area average population density as runoff calculation indexes, and besides, the slow increase form of the function curves can express the influence characteristics of initial loss of rainfall and rainfall infiltration rate on the runoff coefficient. Therefore, the double EXP functions are selected as a runoff coefficient design model, so that the design scale of the rainwater pipe canal and the pump station can be further reasonably quantized and improved, and the safety of the urban drainage waterlogging prevention system is enhanced.
In one embodiment, in step S5, the type of the key node mean sample includes at least one of a corresponding minimum comprehensive runoff coefficient type, a corresponding maximum comprehensive runoff coefficient type of a sparse area of an urban building, and a corresponding comprehensive runoff coefficient type applicable to urban rainwater engineering design; and the value of the key node mean value sample is determined according to corresponding specified critical conditions in the national standard and the intermediate value of the urban building dense area.
Specifically, if the Wuhan city is taken as a target to research a city, the minimum comprehensive runoff coefficient such as an average sample obtained at 2-50 key nodes is controlled according to the corresponding critical condition of 0.20 in the national standard table according to the standard specified in the national standard table. The maximum comprehensive runoff coefficient of the sparse area of the urban building is controlled according to the corresponding critical condition of 0.60 in the national standard table, such as an average sample obtained at 100-50 key nodes. General rainwater engineering design in Wuhan City is suitable for comprehensive runoff coefficients such as an average sample obtained at 2-400 key nodes, the comprehensive runoff coefficients generally adopt 0.60-0.70, and the corresponding values are controlled according to the intermediate value of 0.65 of the urban building compact region in the national standard table. The maximum comprehensive runoff coefficient can reach 1.0 based on the instantaneous comprehensive runoff coefficient as an average sample obtained at a key node of 100-500, and the value is controlled according to the intermediate value of the urban building dense area in the national standard table, wherein the value of the intermediate value is about 0.90.
In one embodiment, the accumulated rainfall index related to each computing node is based on the design rainfall of the city for a long time of 24h, and is distinguished according to a plurality of accumulated rainfall generated in preset time steps.
In one embodiment, in step S6, the constructing step of the target model includes:
and S61, constructing a dual EXP function runoff coefficient model according to the formula (1).
And S62, in the initial parameter adjustment stage, the weighted mean value of the instantaneous comprehensive runoff coefficients generated at each calculation node is obtained through preliminary calculation by randomly setting the value of each metering parameter.
Specifically, the computer device calculates the value of each parameter in the arbitrary setting formula (1) to obtain the instantaneous comprehensive runoff coefficient of each node runoff calculation index, and further determines the weighted average value of each node instantaneous comprehensive runoff coefficient, wherein the weight of the weighted average value is determined by the rainfall of different time step lengths corresponding to each instantaneous comprehensive runoff coefficient.
And S63, when the weighted mean value of the target instantaneous comprehensive runoff coefficient obtained at the key calculation node is determined not to meet the national standard control value, readjusting the value of each metering parameter by a minimum error iterative approximation method until the weighted mean value of the target instantaneous comprehensive runoff coefficient obtained in the corresponding parameter adjustment stage meets the national standard control value, and constructing a target model according to the model metering parameters obtained by adjustment in the current adjustment stage.
Specifically, the computer equipment adjusts the values of all parameters by a minimum error iterative approximation method, so that the weighted average values of instantaneous comprehensive runoff coefficients of 4 key nodes, such as 2-50, 100-50, 2-400, 100-500 and the like, calculated by the double-EXP function runoff coefficient model all accord with the national standard control value (namely, the error is zero), and the set parameter values are the optimal sample values for constructing the target runoff coefficient design model.
In one embodiment, the calculation formula of the dual EXP function runoff coefficient model obtained thereby includes:
(1) The sample calculation for a time step of 5min includes:
(2) The sample calculation for a time step of 60min includes:
determining a node mean sample (see table 1 below) by an arithmetic mean according to the 24h uniform rainfall instantaneous comprehensive runoff coefficients calculated by the formulas (2) and (3):
TABLE 1
In one embodiment, the computer device calculates instantaneous comprehensive runoff coefficients of different runoff calculation indexes with a time step of 5min and a time step of 60min respectively by giving initial values of parameters of the model, then determines weighted mean comprehensive runoff coefficients of nodes to form fitting values corresponding to node mean value samples, further adjusts the parameters of the model by a minimum error iterative approximation method to minimize the average value of absolute values of corresponding data of the two relative error rates, and the calculation formula of the determined double EXP function runoff coefficient design model (namely, the target model) is as follows:
(1) Designing an application calculation formula of a rain type time step length of 5min, comprising the following steps:
(2) Designing an application calculation formula of the rain type time step length of 60min, comprising the following steps:
wherein, in the formula (4) and the formula (5),Ψin order to obtain the instantaneous comprehensive runoff coefficient,Pthe average population density of the urban concentrated construction area,Hand designing the rainfall duration of the rain type for 24h to accumulate the rainfall.
The mean values of the designed rainfall integrated runoff coefficient models calculated by the formulas (4) and (5) are compared with the mean value samples of the uniform rainfall integrated runoff coefficient calculated by the formulas (2) and (3), and the obtained graph distribution comparison result can be exemplarily referred to fig. 3. The model fits the node mean and the error results with the node mean sample, as shown in tables 2 and 3 below. The overall performance is as follows: the graph distribution goodness of fit is high, the relative error rate of the model fitting node mean value total average absolute value is 2.73 percent and is less than 5 percent of the maximum relative error rate standard of engineering design, and the model design meets the requirement of regression precision.
TABLE 2
TABLE 3
In the embodiment, the calculated mean value of the model has a larger elastic range than the value specified by the national standard, the comprehensive runoff coefficient mean value of different rainfall recurrence periods and construction density partitions is provided, the correction deviation range reaches-10% -26.7% compared with the national standard, most of the design conditions belong to the state of value increase, and the correction effect is obvious.
The instantaneous comprehensive runoff coefficient adopted by the rain peak design runoff calculated by the model is remarkably improved in amplitude compared with the value specified in the national standard. In a typical representative area in a city concentrated construction area, the drainage system of the rainwater pipe canal is improved by 30.4 percent on average, and the pumping station pumping drainage waterlogging prevention system is improved by 18.9 percent on average.
Referring to fig. 4, a system 400 for constructing a dual EXP function runoff coefficient design model disclosed in the present application includes an index obtaining module 401, a ranking module 402, and a target model constructing module 403, where:
the index obtaining module 401 is configured to obtain a plurality of node runoff calculation indexes, where the node runoff calculation indexes include a cumulative rainfall index and a regional average population density index.
The grading module 402 is configured to grade the accumulated rainfall index according to a value of a rainfall recurrence period to obtain a graded accumulated rainfall.
The ranking module 402 is further configured to rank the area average population density index according to the area average population density maximum range to obtain a ranked area average population density.
The target model building module 403 is configured to build a dual EXP function runoff coefficient model based on a dual EXP function form, and combine the graded accumulated rainfall and the graded region average population density to determine a model independent variable based on the obtained multiple node runoff calculation indexes.
The target model building module 403 is further configured to determine an instantaneous comprehensive runoff coefficient of each node runoff calculation index based on the dual EXP function runoff coefficient model, use a weighted average of the instantaneous comprehensive runoff coefficients as a node average sample, and screen a key node average sample from the node average sample according to feasibility of national standard control.
The target model building module 403 is further configured to build a target model based on the model parameters in the current adjustment stage when it is determined that the weighted average of the target instantaneous comprehensive runoff coefficient obtained based on the key node average sample corresponds to the national standard control value, otherwise, enter the next parameter adjustment stage, and continuously approach the weighted average of the target instantaneous comprehensive runoff coefficient to the national standard control value by adjusting the model parameters.
In one embodiment, each module in the system is further configured to implement the method in any optional implementation manner of the foregoing embodiment.
Therefore, the system for constructing the dual EXP function runoff coefficient design model provided by the embodiment of the application adopts a regression analysis method based on the national standard control value (namely the national standard comprehensive runoff coefficient control index), and establishes the runoff coefficient design model based on the dual EXP function form and the runoff calculation index with large influence and strong operability, so that under the condition of better connection with the national standard, the traditional fixed value of the comprehensive runoff coefficient can be converted into the instantaneous dynamic value, the reasonability of rainfall whole-course rainfall runoff calculation is improved, the design system for determining the rainwater flow by using a mathematical model method is facilitated and perfected, and the design scientificity of the urban drainage and waterlogging prevention engineering is improved.
The computer program provided by the embodiment of the present application, when executed by a processor, performs the method in any optional implementation manner of the above embodiment. The storage medium may be implemented by any type of volatile or nonvolatile storage device or combination thereof, such as a Static Random Access Memory (SRAM), an Electrically Erasable Programmable Read-Only Memory (EEPROM), an Erasable Programmable Read-Only Memory (EPROM), a Programmable Read-Only Memory (PROM), a Read-Only Memory (ROM), a magnetic Memory, a flash Memory, a magnetic disk, or an optical disk.
The readable storage medium is based on a national standard control value (namely a national standard comprehensive runoff coefficient control index), adopts a regression analysis method, and establishes a runoff coefficient design model based on a double EXP function form and a runoff calculation index with large influence and strong operability, so that under the condition of better connection with the national standard, the traditional fixed value of the comprehensive runoff coefficient can be converted into an instantaneous dynamic value, the reasonability of rainfall whole-course rainwater runoff calculation is improved, a design system for determining rainwater flow by using a mathematical model method is facilitated and perfected, and the scientificity of urban drainage and waterlogging prevention engineering design is improved.
In the embodiments provided in the present application, it should be understood that the disclosed apparatus and method may be implemented in other ways. The above-described embodiments of the apparatus are merely illustrative, and for example, the division of the units is only one logical division, and there may be other divisions when actually implemented, and for example, a plurality of units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection of devices or units through some communication interfaces, and may be in an electrical, mechanical or other form.
In addition, units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of the embodiment.
In addition, functional modules in the embodiments of the present application may be integrated together to form an independent part, or each module may exist alone, or two or more modules may be integrated to form an independent part.
In this document, relational terms such as first and second, and the like may be used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions.
The above description is only an example of the present application and is not intended to limit the scope of the present application, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, improvement and the like made within the spirit and principle of the present application shall be included in the protection scope of the present application.
Claims (7)
1. A method for constructing a dual EXP function runoff coefficient design model is characterized by comprising the following steps:
s1, obtaining a plurality of node runoff calculation indexes, wherein the node runoff calculation indexes comprise accumulated rainfall indexes and regional average population density indexes;
s2, grading the accumulated rainfall index according to the value of the rainfall recurrence period to obtain graded accumulated rainfall;
s3, according to the maximum range of the average population density of the region, carrying out grade division on the average population density index of the region to obtain the average population density of the classified region;
s4, constructing a double-EXP function runoff coefficient model based on a double-EXP function form, combining graded accumulated rainfall and graded region average population density, and determining a model independent variable based on the obtained multiple node runoff calculation indexes;
s5, determining instantaneous comprehensive runoff coefficients of runoff calculation indexes of all nodes based on the double EXP function runoff coefficient model, taking weighted average values of the instantaneous comprehensive runoff coefficients as node average value samples, and screening key node average value samples from the node average value samples according to feasibility of national standard control;
s6, when the weighted mean value of the target instantaneous comprehensive runoff coefficient correspondingly obtained based on the key node mean value sample is determined to accord with the national standard control value, constructing a target model based on the model parameter in the current adjusting stage as the parameter of the model sample calculation formula, otherwise, entering the next parameter adjusting stage, and continuously approaching the weighted mean value of the target instantaneous comprehensive runoff coefficient to the national standard control value by adjusting the model parameter;
in the step S4, the node runoff calculation index takes a three-parameter power function form as a function expansion;
the calculation formula of the dual EXP function runoff coefficient model comprises:
wherein the content of the first and second substances,
in order to obtain the instantaneous comprehensive runoff coefficient,Pis the average population density of the classified area with the unit of people/hm 2 ，HThe rainfall is accumulated in grades and the unit is mm,
、
、
for the metering parameters covered in the resulting three-parameter power function developed through the rank region average population density,
、
、
(ii) a metric parameter covered in a three-parameter power function developed by the graded accumulated rainfall;
in step S6, the constructing step of the target model includes:
s61, constructing a dual EXP function runoff coefficient model according to a formula (1);
s62, in the initial parameter adjusting stage, the weighted mean value of the instantaneous comprehensive runoff coefficients generated at each calculation node is obtained through preliminary calculation by randomly setting the value of each metering parameter;
and S63, when the weighted average of the target instantaneous comprehensive runoff coefficient obtained at the key calculation node is determined not to meet the national standard control value, readjusting the value of each metering parameter by a minimum error iterative approximation method until the weighted average of the target instantaneous comprehensive runoff coefficient obtained in the corresponding parameter adjustment stage meets the national standard control value, and constructing a target model according to the model metering parameters obtained by adjustment in the current adjustment stage.
2. The method according to claim 1, wherein in step S2, the step of performing a grading on the index of accumulated rainfall according to the value of the rainfall recurrence period to obtain a graded accumulated rainfall includes:
s21, setting values of drainage waterlogging prevention design recurrence periods of urban areas based on outdoor drainage design standards, and grading the accumulated rainfall indexes according to values of rainfall recurrence periods to obtain a plurality of graded rainfall indexes, wherein the rainfall is designed according to rainstorm of a preset time length of an area, so that the accumulated rainfall of each rainfall recurrence period can be determined.
3. The method according to claim 1, wherein the step S3 of ranking the area average population density index according to the area average population density maximum range to obtain the graded area average population density comprises:
and S31, carrying out grade division on the area average population density index based on a maximum interval range obtained by average population density in a research area to obtain a plurality of graded population density indexes, wherein the minimum value in the maximum interval range is used as a starting unit, and each grade is gradually increased according to a preset length unit until the maximum value in the maximum interval range is reached.
4. The method according to claim 1, wherein in step S5, the type of the key node mean value sample includes at least one of a corresponding minimum comprehensive runoff coefficient type, a corresponding maximum comprehensive runoff coefficient type, a corresponding urban building sparse area maximum comprehensive runoff coefficient type, and a corresponding urban rainwater engineering design applicable comprehensive runoff coefficient type;
and the value of the key node mean value sample is determined according to the corresponding specified critical condition in the national standard and the intermediate value of the urban building dense area.
5. The method of claim 4, wherein the indication of cumulative rainfall involved at each computing node is based on a 24h design rain profile for a city, differentiated by a plurality of cumulative rainfall events occurring in predetermined time steps.
6. The system for constructing the dual EXP function runoff coefficient design model is characterized by comprising an index acquisition module, a grading module and a target model construction module, wherein:
the index acquisition module is used for acquiring a plurality of node runoff calculation indexes, and the node runoff calculation indexes comprise accumulated rainfall indexes and regional average population density indexes;
the grading module is used for grading the accumulated rainfall index according to the value of the rainfall recurrence period to obtain graded accumulated rainfall;
the grading module is further used for grading the area average population density index according to the area average population density maximum range to obtain the grading area average population density;
the target model building module is used for building a double EXP function runoff coefficient model based on a double EXP function form, combining graded accumulated rainfall and graded region average population density, and determining a model independent variable based on the obtained multiple node runoff calculation indexes;
the target model building module is further used for determining an instantaneous comprehensive runoff coefficient of each node runoff calculation index based on the double EXP function runoff coefficient model, taking a weighted average value of the instantaneous comprehensive runoff coefficients as a node average value sample, and screening key node average value samples from the node average value sample according to feasibility of national standard control;
the target model building module is also used for building a target model based on the model parameters in the current adjusting stage when the weighted mean value of the target instantaneous comprehensive runoff coefficient correspondingly obtained based on the key node mean value sample is determined to accord with the national standard control value, otherwise, entering the next parameter adjusting stage, and enabling the weighted mean value of the target instantaneous comprehensive runoff coefficient to continuously approach the national standard control value by adjusting the model parameters;
the node runoff calculation index takes a three-parameter power function form as a function expansion;
the calculation formula of the dual EXP function runoff coefficient model comprises:
wherein the content of the first and second substances,
in order to obtain the instantaneous comprehensive runoff coefficient,Pis the average population density of the classified area with the unit of people/hm 2 ，HThe rainfall is accumulated in grades and the unit is mm,
、
、
for the metering parameters covered in the resulting three-parameter power function developed through the rank region average population density,
、
、
(ii) a metric parameter covered in a three-parameter power function developed by the graded accumulated rainfall;
the specific implementation of the construction of the target model is as follows:
constructing a double EXP function runoff coefficient model according to the formula (1);
in the initial parameter adjustment stage, the weighted mean value of the instantaneous comprehensive runoff coefficients generated at each calculation node is obtained through preliminary calculation by randomly setting the value of each metering parameter;
and when the weighted mean value of the target instantaneous comprehensive runoff coefficient obtained at the key calculation node is determined not to meet the national standard control value, readjusting the value of each metering parameter by a minimum error iterative approximation method until the weighted mean value of the target instantaneous comprehensive runoff coefficient obtained in the corresponding parameter adjustment stage meets the national standard control value, and constructing a target model according to the model metering parameter obtained by adjustment in the current adjustment stage.
7. A readable storage medium, characterized in that the readable storage medium comprises a method program for constructing a dual EXP function runoff coefficient design model, which when executed by a processor implements the steps of the method according to any one of claims 1 to 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211179274.2A CN115270525B (en) | 2022-09-27 | 2022-09-27 | Method and system for constructing dual EXP function runoff coefficient design model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211179274.2A CN115270525B (en) | 2022-09-27 | 2022-09-27 | Method and system for constructing dual EXP function runoff coefficient design model |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115270525A CN115270525A (en) | 2022-11-01 |
| CN115270525B true CN115270525B (en) | 2023-01-24 |
Family
ID=83756051
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211179274.2A Active CN115270525B (en) | 2022-09-27 | 2022-09-27 | Method and system for constructing dual EXP function runoff coefficient design model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115270525B (en) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114971987A (en) * | 2022-05-31 | 2022-08-30 | 武汉市政工程设计研究院有限责任公司 | Typical year comparison and selection method and system for urban rainfall and readable storage medium |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105808948B (en) * | 2016-03-08 | 2017-02-15 | 中国水利水电科学研究院 | Automatic correctional multi-mode value rainfall ensemble forecast method |
| CN106651211B (en) * | 2016-12-30 | 2021-02-02 | 吉林师范大学 | Flood disaster risk assessment method for different scale areas |
| CN108108838B (en) * | 2017-12-18 | 2021-08-27 | 华电福新能源股份有限公司福建分公司 | Season regulation reservoir optimal scheduling method with high water utilization rate |
| CN109657948A (en) * | 2018-12-04 | 2019-04-19 | 昆山市建设工程质量检测中心 | Plot sponge urban construction performance comprehensive estimation method |
| CN110188323B (en) * | 2019-05-15 | 2020-09-22 | 华南理工大学 | Water surface rate value range measuring and calculating method suitable for urban water ecological protection |
| EP3839870A1 (en) * | 2019-12-20 | 2021-06-23 | Universidade de Trás-os-Montes e Alto Douro | Methodology to determine locations of systems for rainwater harvesting in hydrographic basins |
| CN112801342A (en) * | 2020-12-31 | 2021-05-14 | 国电大渡河流域水电开发有限公司 | Adaptive runoff forecasting method based on rainfall runoff similarity |
| CN114372685B (en) * | 2021-12-28 | 2024-06-11 | 长江生态环保集团有限公司 | Urban storm waterlogging risk assessment method based on SWMM model |
-
2022
- 2022-09-27 CN CN202211179274.2A patent/CN115270525B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114971987A (en) * | 2022-05-31 | 2022-08-30 | 武汉市政工程设计研究院有限责任公司 | Typical year comparison and selection method and system for urban rainfall and readable storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115270525A (en) | 2022-11-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111582755B (en) | Mountain torrent disaster comprehensive risk dynamic assessment method based on multi-dimensional set information | |
| Lettenmaier et al. | Climatic sensitivity of California water resources | |
| Rajsekhar et al. | Regionalization of drought characteristics using an entropy approach | |
| Kim et al. | Hydrological modeling and evaluation of rainwater harvesting facilities: case study on several rainwater harvesting facilities in Korea | |
| Hussain et al. | Modeling the effects of land-use and climate change on the performance of stormwater sewer system using SWMM simulation: case study | |
| CN111581767B (en) | Calibrating method for checking characteristic parameters of pipe network-river coupling model | |
| Ruiz-Villanueva et al. | Climate change impacts on discharges of the Rhone River in Lyon by the end of the twenty-first century: model results and implications | |
| Fewkes et al. | Method of modelling the performance of rainwater collection systems in the United Kingdom | |
| CN115130396A (en) | Distributed hydrological model modeling method for riverway type reservoir area | |
| Ledingham et al. | Contrasting seasonality of storm rainfall and flood runoff in the UK and some implications for rainfall-runoff methods of flood estimation | |
| CN112904458B (en) | Hydrological forecasting method and system for super-seepage-full-storage mixed runoff yield mode | |
| CN115759459A (en) | Intelligent drainage basin flood flow set probability forecasting method based on rainwater engineering condition coupling | |
| Mattos et al. | Improving urban flood resilience under climate change scenarios in a tropical watershed using low-impact development practices | |
| Ding et al. | Assessment of the impact of climate change on urban flooding: A case study of Beijing, China | |
| Saedi et al. | Optimal water allocation of the Zayandeh-Roud Reservoir in Iran based on inflow projection under climate change scenarios | |
| CN115270525B (en) | Method and system for constructing dual EXP function runoff coefficient design model | |
| Lai et al. | Development of a polder module in the SWAT model: SWATpld for simulating polder areas in south‐eastern China | |
| Zhang et al. | Runoff prediction under different precipitation scenarios based on SWAT model and stochastic simulation of precipitation | |
| Liu et al. | Urban waterlogging control: A novel method to urban drainage pipes reconstruction, systematic and automated | |
| Wang et al. | Application of SWAT model with CMADS data for hydrological simulation in western China | |
| Gan et al. | Hydrological process simulation of sluice-controlled rivers in the plains area of China based on an improved SWAT model | |
| Singson et al. | Modeling climate change impact on the inflow of the Magat reservoir using the Soil and Water Assessment Tool (SWAT) model for dam management | |
| Moravej et al. | Site-scale Urban Water Mass Balance Assessment (SUWMBA) Tool (Version 2) | |
| Yigzaw et al. | Land use and land cover impact on probable maximum flood and sedimentation for artificial reservoirs: case study in the Western United States | |
| CN111782711A (en) | Water quantity data statistical method |
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 |