CN115510600A - Urban drainage network optimization design method - Google Patents

Urban drainage network optimization design method Download PDF

Info

Publication number
CN115510600A
CN115510600A CN202211408100.9A CN202211408100A CN115510600A CN 115510600 A CN115510600 A CN 115510600A CN 202211408100 A CN202211408100 A CN 202211408100A CN 115510600 A CN115510600 A CN 115510600A
Authority
CN
China
Prior art keywords
drainage
pipeline
urban
drainage network
model
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202211408100.9A
Other languages
Chinese (zh)
Other versions
CN115510600B (en
Inventor
王颖
兰云龙
李进
刘峰
阳利锋
任祥
何斌
李�浩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chuang Hui Da Design Co ltd
Original Assignee
Chuang Hui Da Design Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Chuang Hui Da Design Co ltd filed Critical Chuang Hui Da Design Co ltd
Priority to CN202211408100.9A priority Critical patent/CN115510600B/en
Publication of CN115510600A publication Critical patent/CN115510600A/en
Application granted granted Critical
Publication of CN115510600B publication Critical patent/CN115510600B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/04—Constraint-based CAD
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • 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
    • Y02A30/00—Adapting or protecting infrastructure or their operation
    • Y02A30/60—Planning or developing urban green infrastructure

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Computational Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Sewage (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention relates to the technical field of drainage network optimization, and discloses an urban drainage network optimization design method and device, wherein the method comprises the following steps: collecting urban drainage network data to form an urban drainage network structure, and correcting the urban drainage network structure; constructing an urban drainage network model based on the corrected urban drainage network structure; constructing an urban drainage network optimization model, and determining an urban drainage network design constraint specification expression; and optimizing and solving the constructed urban drainage network optimization model by using a gradient optimization algorithm, and placing the optimized drainage network in the urban drainage network model for operation. The method comprises the steps of determining the pressure difference between two ends of a drainage pipeline based on the flow to be discharged of the drainage pipeline, further constructing a drainage network optimization design model, and solving the model to obtain an urban drainage network optimization design scheme, namely reducing the operating pressure of the urban drainage network and reducing the frequency of urban inland inundation disasters by newly adding a plurality of drainage pipelines at specific positions.

Description

Urban drainage network optimization design method
Technical Field
The invention relates to the technical field of drainage network optimization, in particular to an urban drainage network optimization design method.
Background
Along with the acceleration of the urbanization process, the proportion of impermeable substances such as asphalt, concrete and the like in the urban underlying surface is continuously increased, the area of permeable soil is reduced, more soil is replaced by roads, and if heavy rainfall occurs, the rainfall and rainfall can only be discharged through drainage pipelines, so that the urban hydrological condition is influenced, the operating pressure of an urban drainage network is increased, and the urban inland inundation disaster occurrence frequency is obviously improved.
Disclosure of Invention
In view of this, the present invention provides an optimal design method and device for an urban drainage network, which aims to improve the utilization rate and safety of the urban drainage network and improve the urban drainage capacity.
The invention provides an urban drainage network optimization design method and device, aiming at: 1) Forming a drainage pipeline structure in the urban drainage network based on a search query mode, deleting an isolated drainage pipeline, determining an urban drainage network model based on urban precipitation, soil precipitation underwater seepage, precipitation evaporation capacity and maximum water storage capacity, wherein most of precipitation is infiltrated into soil and is stored and evaporated in a depression, and the rest of precipitation needs to be discharged by utilizing the drainage pipeline in the urban drainage network, and calculating to-be-discharged flow of drainage pipelines in different areas according to the determined urban drainage network model; 2) The method comprises the steps of determining the pressure difference between two ends of different drainage pipelines based on the flow to be discharged of drainage pipelines in different areas, if the pressure difference between two ends of the drainage pipelines is too large, equivalently increasing the operating pressure of an urban drainage network, the drainage pipelines are likely to be damaged, the frequency of urban waterlogging disasters is obviously improved, constructing a drainage network optimization design model based on the pressure difference between two ends of the drainage pipelines, solving the model by using a gradient optimization algorithm to obtain an urban drainage network optimization design scheme, namely, reducing the operating pressure of the urban drainage network and reducing the frequency of the urban waterlogging disasters by newly adding a plurality of drainage pipelines in specific positions.
The invention provides an optimal design method of an urban drainage network, which comprises the following steps:
s1: collecting urban drainage network data to form an urban drainage network structure, correcting the urban drainage network structure, and delaying the connection of a network topological relation;
s2: constructing an urban drainage network model based on a modified urban drainage network structure, wherein the urban drainage network model comprises the flow of different urban drainage pipelines and the pressure difference between two ends of the pipelines, and screening parameters in the model by using a modified Morris screening method to form an example of the urban drainage network model;
s3: constructing an urban drainage network optimization model, and determining an urban drainage network design constraint specification expression, wherein the solution result of the urban drainage network optimization model is the position of a newly added drainage pipeline in the urban drainage network;
s4: and optimizing and solving the constructed urban drainage network optimization model by using a gradient optimization algorithm, optimizing the urban drainage network structure according to the solving result, placing the optimized drainage network in the urban drainage network model for operation, implementing if the effect reaches a specified threshold value, and otherwise, solving again, wherein a rank 2 quasi Newton method based on a non-precise Armijo criterion is the gradient optimization method.
As a further improvement of the method of the invention:
optionally, the collecting urban drainage network data in the step S1 forms an urban drainage network structure, including:
collecting urban drainage network data, wherein the urban drainage network data comprise position parameters of urban drainage pipelines, pipeline lengths and pipeline wall roughness coefficients, constructing an urban drainage network structure based on the urban drainage network data, the urban drainage network structure comprises position distribution of the urban drainage pipelines and pipeline parameters, the pipeline parameters are the pipeline lengths and the pipeline wall roughness coefficients, and connecting pipe orifice positions and pipe bottom positions of the same pipeline by straight lines to form position distribution of the pipeline;
the drainage pipeline position parameters comprise the pipe orifice position and the pipe bottom position of the urban drainage pipeline, the urban drainage pipeline is straight, and the pipe orifice position is higher than the pipe bottom position;
the urban drainage network is provided with n drainage pipelines in total, and the data set of the n drainage pipelines is as follows:
Figure 451122DEST_PATH_IMAGE001
wherein:
Figure 650154DEST_PATH_IMAGE002
indicating any of the second place in a municipal drainage network
Figure 961049DEST_PATH_IMAGE003
A water draining pipeline is arranged on the upper portion of the water tank,
Figure 887417DEST_PATH_IMAGE004
is shown as
Figure 498527DEST_PATH_IMAGE003
The position coordinates of the pipe orifice of the strip drainage pipeline,
Figure 602880DEST_PATH_IMAGE005
denotes the first
Figure 401072DEST_PATH_IMAGE003
The pipe bottom position coordinates of the strip drainage pipeline are the central position coordinates of the pipe opening or the pipe bottom area;
in the embodiment of the invention, the city is constructed into a three-dimensional coordinate system, the center of the surface of the city is the origin of coordinates of the three-dimensional coordinate system, the east-west direction is the X-axis direction, the south-north direction is the Y-axis direction, and the direction from the ground bottom to the surface of the city is the Z-axis direction;
Figure 131131DEST_PATH_IMAGE006
is shown as
Figure 65589DEST_PATH_IMAGE003
The length of the pipeline of the strip drainage pipeline,
Figure 606423DEST_PATH_IMAGE007
denotes the first
Figure 626331DEST_PATH_IMAGE003
The pipe wall roughness coefficient of the strip drainage pipe.
Optionally, the step S1 of correcting the urban drainage network structure includes:
the correction process of the urban drainage network structure comprises the following steps:
s11: deleting the position coordinates of the isolated pipe orifice or the position coordinates of the pipe bottom in the constructed urban drainage network structure;
s12: selecting a drainage pipeline with a pipe orifice position coordinate on the ground surface for inquiry;
s13: if the pipe orifice position coordinate of only one drainage pipeline is connected with the pipe bottom position coordinate of the selected drainage pipeline after being inquired, connecting the inquired drainage pipeline and the selected drainage pipeline, taking the inquired drainage pipeline as the selected pipeline, and repeating the step S13 until the drainage pipeline cannot be inquired;
if the distance between the pipe orifice position coordinates of the inquired drainage pipeline and the selected drainage pipeline and the pipe bottom position coordinates is smaller than a distance threshold value, judging that the two position coordinates are connected;
s14: if the pipe orifice position coordinates of a plurality of drainage pipelines are inquired to be connected with the pipe bottom position coordinates of the selected drainage pipeline, respectively connecting the inquired drainage pipeline and the selected drainage pipeline, respectively taking the inquired drainage pipeline as the selected pipeline, and returning to the step S13 until the drainage pipeline cannot be inquired;
s15: deleting the drainage pipeline of which the position coordinate of the pipe orifice is not positioned on the ground surface and the position coordinate of the pipe orifice is not connected with the position coordinate of the pipe bottom of other drainage pipelines in the urban drainage network structure;
based on the correction process of the urban drainage network structure
Figure 160081DEST_PATH_IMAGE008
The modified urban drainage network structure of the drainage pipeline,
Figure 949045DEST_PATH_IMAGE008
indicating the number of modified urban drainage pipelines.
Optionally, in the step S2, a city drainage network model is constructed based on the modified city drainage network structure, and the method includes:
dividing the urban surface into a plurality of different sub-areas based on the modified urban drainage network structure, wherein each sub-area is provided with a drainage pipeline with a pipe orifice position coordinate positioned on the ground surface, the number of the divided sub-areas is m, and the number of the first sub-areas is
Figure 660780DEST_PATH_IMAGE009
Sub-regions are
Figure 433564DEST_PATH_IMAGE010
，
Figure 239846DEST_PATH_IMAGE010
The position coordinate of the middle pipe orifice is positioned on the drainage pipeline on the ground surface
Figure 414476DEST_PATH_IMAGE011
Said sub-region
Figure 297112DEST_PATH_IMAGE012
The drainage network structure of (1) is:
Figure 291613DEST_PATH_IMAGE013
wherein:
Figure 901586DEST_PATH_IMAGE014
indicating water flow from the drain line
Figure 930722DEST_PATH_IMAGE011
Starting from, can pass through
Figure 249839DEST_PATH_IMAGE015
A strip drainage pipeline;
the drainage pipeline
Figure 466056DEST_PATH_IMAGE011
The deepest drain pipe connected to the city drain network structure is integrated into
Figure 145299DEST_PATH_IMAGE016
The deepest drainage pipeline is a pipeline with the position coordinate of the bottom of the pipeline not connected with other drainage pipelines, and water on the ground surface can pass through the drainage pipelines
Figure 763363DEST_PATH_IMAGE011
Arrive at
Figure 265100DEST_PATH_IMAGE016
The middle and any drainage pipelines and the deepest drainage pipeline are finally communicated with a sewage treatment plant;
the urban drainage network model comprises drainage network models of different sub-areas, and the sub-areas
Figure 234193DEST_PATH_IMAGE012
Is drainedThe network model is as follows:
Figure 185968DEST_PATH_IMAGE017
Figure 189696DEST_PATH_IMAGE018
Figure 585037DEST_PATH_IMAGE019
Figure 41426DEST_PATH_IMAGE020
wherein:
Figure 796892DEST_PATH_IMAGE021
representing an area per unit time
Figure 389548DEST_PATH_IMAGE012
The amount of rainfall of;
Figure 955789DEST_PATH_IMAGE022
representing an area per unit time
Figure 899474DEST_PATH_IMAGE012
The average infiltration capacity of the soil in the middle permeable area;
Figure 458632DEST_PATH_IMAGE023
representing an area per unit time
Figure 171373DEST_PATH_IMAGE012
Maximum water storage capacity of the middle impervious area;
Figure 642937DEST_PATH_IMAGE024
representing an area per unit time
Figure 73918DEST_PATH_IMAGE012
The rainfall evaporation capacity of (d);
Figure 171187DEST_PATH_IMAGE025
representing an area per unit time
Figure 4014DEST_PATH_IMAGE012
The precipitation production rate, i.e. the drainage pipeline
Figure 646479DEST_PATH_IMAGE026
The flow rate of the discharged water;
Figure 564756DEST_PATH_IMAGE027
indicating drainage pipe
Figure 465716DEST_PATH_IMAGE011
The pressure difference between the two ends of the pipeline,
Figure 887470DEST_PATH_IMAGE028
which is indicative of the density of the water,
Figure 966416DEST_PATH_IMAGE029
which represents the acceleration of the force of gravity,
Figure 106410DEST_PATH_IMAGE030
indicating precipitation flowing into drainage pipe from surface
Figure 811061DEST_PATH_IMAGE011
The length of the distance the orifice is to flow through,
Figure 352901DEST_PATH_IMAGE031
indicating drainage pipe
Figure 602748DEST_PATH_IMAGE011
The pipe wall roughness coefficient of (a);
Figure 230038DEST_PATH_IMAGE032
connecting drainage pipeline
Figure 738380DEST_PATH_IMAGE033
All upstream drainage pipes of (2) are to the drainage pipe in unit time
Figure 885458DEST_PATH_IMAGE034
Total flow of discharged water, said upstream discharge pipe being positioned higher than said discharge pipe
Figure 555474DEST_PATH_IMAGE033
Can supply water to the drainage pipeline
Figure 670061DEST_PATH_IMAGE033
A drain pipeline for discharging sewage;
Figure 450935DEST_PATH_IMAGE035
indicating area
Figure 701788DEST_PATH_IMAGE012
The ratio of the road area of (a),
Figure 559016DEST_PATH_IMAGE036
indicating area
Figure 160899DEST_PATH_IMAGE012
The area of the forest land is used up,
Figure 745464DEST_PATH_IMAGE037
indicating area
Figure 850823DEST_PATH_IMAGE012
The area ratio of the grass;
in the embodiment of the invention, if no branch drainage pipeline exists between the connected drainage pipelines, the drainage flow of the connected drainage pipelines is the same, and if the connected drainage pipelines do not exist, the drainage flow of the connected drainage pipelines is the sameBetween are
Figure 878953DEST_PATH_IMAGE038
With branched drainage pipes, the flow of deeper positioned drainage pipes being the total flow of the upstream drainage
Figure 702553DEST_PATH_IMAGE039
。
Optionally, the step S2 of screening parameters in the model by using a modified Morris screening method to determine the urban drainage network model includes:
the model parameters comprise the average infiltration capacity and the maximum water storage capacity of the soil in different subregions, the rainfall evaporation capacity and the rainfall in unit time;
screening parameters in the model by using a modified Morris screening method, wherein the model parameter screening process of the modified Morris screening method comprises the following steps:
s21: randomly selecting one parameter from model parameters to modify for many times, wherein the modification mode of the model parameters is to collect rainfall conditions of the same city at different times as the model parameters;
s22: inputting the parameter result after each modification into the urban drainage network model to obtain a plurality of groups of surface drainage pipeline flow results in the urban drainage network model;
s23: calculating sensitivity of modified parameters
Figure 90809DEST_PATH_IMAGE040
：
Figure 66986DEST_PATH_IMAGE041
Wherein:
the number of times of modifying the parameters is K-1;
Figure 249706DEST_PATH_IMAGE042
representing a drainage pipeline flow result obtained based on the kth modified parameter;
Figure 826181DEST_PATH_IMAGE043
representing the change rate of the k-th modified parameter compared with the initial parameter value;
s24: and repeating the steps to obtain the sensitivity of different model parameters, and deleting the parameters with the sensitivity absolute value lower than 0.02 in the model to obtain the urban drainage network model suitable for the selected city.
Optionally, the constructing the urban drainage network optimization model in the step S3 includes:
the constructed urban drainage network optimization model B is as follows:
Figure 752548DEST_PATH_IMAGE044
wherein:
Figure 583232DEST_PATH_IMAGE045
indicates the added first
Figure 202433DEST_PATH_IMAGE046
The coordinate position of the bottom of the strip drainage pipeline,
Figure 624DEST_PATH_IMAGE047
indicating the coordinate location of any sewage treatment plant,
Figure 730683DEST_PATH_IMAGE048
the distance between the bottom of the h-th drainage pipeline and the nearest sewage treatment plant is shown;
Figure 415873DEST_PATH_IMAGE049
indicates the increased second
Figure 205975DEST_PATH_IMAGE050
The method comprises the following steps that the original drainage flow of the opening position of a drainage pipeline is shown, and the original drainage flow represents the drainage flow when the drainage pipeline is not increased;
Figure 491463DEST_PATH_IMAGE051
showing the current drainage flow of the pipe orifice after the h-th drainage pipeline is added;
h represents the number of newly added drainage pipelines in the urban drainage network optimization design scheme;
the solving result of the urban drainage network optimization model is the position of a newly added drainage pipeline in the urban drainage network;
determining a standard expression of urban drainage network design constraints, wherein the drainage network design constraints are as follows:
Figure 759633DEST_PATH_IMAGE052
Figure 564909DEST_PATH_IMAGE053
wherein:
Figure 791491DEST_PATH_IMAGE008
representing the number of drainage pipes in the urban drainage network to be optimized.
In the specific embodiment of the invention, the pipe orifice position of the newly added drainage pipeline is the pipe bottom position of other drainage pipelines.
Optionally, the step S4 of performing optimization solution on the constructed urban drainage network optimization model by using a gradient optimization algorithm includes:
the method comprises the following steps of utilizing a gradient optimization algorithm to carry out optimization solution on a constructed urban drainage network optimization model, wherein a rank 2 quasi Newton method based on a non-precise Armijo criterion is a main method of the gradient optimization algorithm, and the solution process of the urban drainage network optimization model comprises the following steps:
s41: constructing a penalty function form of an urban drainage network optimization model and a constraint specification expression:
Figure 564275DEST_PATH_IMAGE054
Figure 386868DEST_PATH_IMAGE055
wherein:
Figure 827077DEST_PATH_IMAGE056
the position distribution of newly added H drainage pipelines in the urban drainage network optimization design scheme is shown, and in the embodiment of the invention, H is set to be
Figure 958981DEST_PATH_IMAGE008
；
Figure 969794DEST_PATH_IMAGE057
To the penalty function coefficient, set it to 10;
s42: position distribution of initially generated H newly-added drainage pipelines
Figure 845346DEST_PATH_IMAGE058
And setting an allowable error of
Figure 140061DEST_PATH_IMAGE059
The initial positive definite matrix is
Figure 181880DEST_PATH_IMAGE060
，
Figure 929256DEST_PATH_IMAGE061
The method is characterized in that the method is an identity matrix, the current iteration number of the algorithm is w, and the initial value of w is 0;
s43: calculating a penalty function in
Figure 608499DEST_PATH_IMAGE062
Gradient of (2)
Figure 508453DEST_PATH_IMAGE063
If, if
Figure 982160DEST_PATH_IMAGE064
Then the position distribution obtained by the w-th iteration
Figure 216832DEST_PATH_IMAGE062
That is, the optimization solution result is obtained, otherwise, the step S44 is turned to;
s44: computing
Figure 919340DEST_PATH_IMAGE065
；
S45: updated using the Armijo criteria as follows
Figure 188647DEST_PATH_IMAGE066
：
Figure 567676DEST_PATH_IMAGE067
Figure 24065DEST_PATH_IMAGE068
Wherein:
Figure 530264DEST_PATH_IMAGE069
is composed of
Figure 388498DEST_PATH_IMAGE070
The update step length of (2);
Figure 204008DEST_PATH_IMAGE071
to update the parameters, it is set to 0.5;
Figure 147693DEST_PATH_IMAGE072
to satisfy
Figure 192003DEST_PATH_IMAGE073
The minimum trace distribution matrix of (a), T represents transposition;
s46: is updated to
Figure 904745DEST_PATH_IMAGE074
：
Figure 891155DEST_PATH_IMAGE075
S47: order to
Figure 322136DEST_PATH_IMAGE076
The process returns to step S43.
Optionally, in the step S4, the optimized drainage network is placed in an urban drainage network model to operate, and if the effect reaches a specified threshold, the implementation is performed, otherwise, the re-solution is performed, including:
the optimized urban drainage network topological structure is placed in an urban drainage network model, precipitation data of different seasons of a city are collected and input into the model, the model outputs the flow of each drainage pipeline in the urban drainage network and the pressure difference between two ends of each pipeline, if the pressure difference between two ends of each drainage pipeline is smaller than a specified threshold value, an optimized design scheme of the urban drainage network can be implemented, namely the load pressure of the drainage network is reduced by adding a plurality of drainage pipelines, and otherwise, a gradient optimization algorithm is reused for solving.
In order to solve the above problems, the present invention provides an optimal design device for an urban drainage network, comprising:
the drainage network construction device is used for collecting urban drainage network data to form an urban drainage network structure, correcting the urban drainage network structure, constructing an urban drainage network model based on the corrected urban drainage network structure, and screening parameters in the model by using a corrected Morris screening method to form the urban drainage network model;
the optimization model construction device is used for constructing an urban drainage network optimization model and determining an urban drainage network design constraint specification expression;
and the drainage network optimization module is used for optimizing and solving the constructed urban drainage network optimization model by using a gradient optimization algorithm, optimizing the urban drainage network structure according to a solving result, placing the optimized drainage network in the urban drainage network model for operation, implementing a corresponding urban drainage network optimization design scheme if the effect reaches a specified threshold value, and otherwise, re-solving.
Compared with the prior art, the invention provides an optimal design method of an urban drainage network, which has the following advantages:
firstly, the scheme provides an urban drainage network model, wherein the urban surface is divided into a plurality of different sub-regions, each sub-region is provided with a drainage pipeline with a pipe orifice position coordinate positioned on the ground surface, the urban drainage network model comprises drainage network models of different sub-regions, and the sub-regions are provided with drainage network models of different sub-regions
Figure 170138DEST_PATH_IMAGE077
The drainage network model is as follows:
Figure 737385DEST_PATH_IMAGE017
Figure 894697DEST_PATH_IMAGE018
Figure 812975DEST_PATH_IMAGE019
Figure 199088DEST_PATH_IMAGE020
wherein:
Figure 886421DEST_PATH_IMAGE021
representing an area per unit time
Figure 214634DEST_PATH_IMAGE012
The amount of rainfall of;
Figure 370940DEST_PATH_IMAGE022
representing an area per unit time
Figure 810012DEST_PATH_IMAGE012
The average infiltration capacity of the soil in the middle permeable area;
Figure 351852DEST_PATH_IMAGE023
represents the area in a unit time
Figure 850966DEST_PATH_IMAGE012
Maximum water storage capacity of the middle impervious area;
Figure 478257DEST_PATH_IMAGE024
representing an area per unit time
Figure 206172DEST_PATH_IMAGE012
The rainfall evaporation capacity of (d);
Figure 868098DEST_PATH_IMAGE078
representing an area per unit time
Figure 538113DEST_PATH_IMAGE012
The precipitation yield of, i.e. said drainage pipe
Figure 652700DEST_PATH_IMAGE026
The flow rate of the discharged water;
Figure 184307DEST_PATH_IMAGE027
indicating drainage pipe
Figure 700739DEST_PATH_IMAGE011
The pressure difference between the two ends of the pipeline,
Figure 541656DEST_PATH_IMAGE028
which is indicative of the density of the water,
Figure 159850DEST_PATH_IMAGE029
which represents the acceleration of the force of gravity,
Figure 478836DEST_PATH_IMAGE030
indicating precipitation flowing into drainage pipe from surface
Figure 115354DEST_PATH_IMAGE011
The length of the distance the orifice is to flow through,
Figure 127172DEST_PATH_IMAGE031
indicating drainage pipe
Figure 701504DEST_PATH_IMAGE011
The pipe wall roughness coefficient of (a);
Figure 824181DEST_PATH_IMAGE079
connecting drainage pipeline
Figure 49626DEST_PATH_IMAGE034
All upstream drainage pipes of (2) are to the drainage pipe in unit time
Figure 232345DEST_PATH_IMAGE033
The total flow of discharged water is that the upstream water discharge pipeline is positioned higher than the water discharge pipeline
Figure 808820DEST_PATH_IMAGE034
Can supply water to the drainage pipeline
Figure 485920DEST_PATH_IMAGE034
A drainage pipeline for draining sewage;
Figure 831451DEST_PATH_IMAGE035
indicating area
Figure 185072DEST_PATH_IMAGE012
The ratio of the road area of (a),
Figure 999575DEST_PATH_IMAGE036
indicating area
Figure 464055DEST_PATH_IMAGE012
The area of the forest land is in proportion,
Figure 664092DEST_PATH_IMAGE037
indicating area
Figure 454193DEST_PATH_IMAGE012
The ratio of the grass land area; compared with the traditional scheme, the method has the advantages that the structure of the drainage pipelines in the urban drainage network is formed based on the search query mode, the isolated drainage pipelines are deleted, the urban drainage network model is determined based on the precipitation amount, the soil precipitation underwater seepage amount, the precipitation evaporation amount and the maximum water storage amount of the city, most of the precipitation is infiltrated into the soil, is stored in the hollow land and is evaporated, the rest of the precipitation needs to be discharged by utilizing the drainage pipelines in the urban drainage network, and the flow to be discharged of the drainage pipelines in different areas is obtained through calculation according to the determined urban drainage network model.
Meanwhile, the scheme provides an urban drainage network optimization model, and the constructed urban drainage network optimization model B is as follows:
Figure 224834DEST_PATH_IMAGE044
wherein:
Figure 24163DEST_PATH_IMAGE045
indicates the increased second
Figure 78707DEST_PATH_IMAGE046
The coordinate position of the bottom of the strip drainage pipeline,
Figure 790442DEST_PATH_IMAGE047
indicating the coordinate location of any sewage treatment plant,
Figure 828805DEST_PATH_IMAGE080
the distance between the bottom of the h-th drainage pipeline and the nearest sewage treatment plant is shown;
Figure 900666DEST_PATH_IMAGE081
indicates the increased second
Figure 103326DEST_PATH_IMAGE050
The method comprises the following steps that the original drainage flow of the opening position of a drainage pipeline is shown, and the original drainage flow represents the drainage flow when the drainage pipeline is not increased;
Figure 500809DEST_PATH_IMAGE051
showing the drainage flow of the pipe orifice position after the h-th drainage pipeline is added: h represents the number of newly added drainage pipelines in the urban drainage network optimization design scheme; the solving result of the urban drainage network optimization model is the position of a newly added drainage pipeline in the urban drainage network; determining a standard expression of urban drainage network design constraints, wherein the drainage network design constraints are as follows:
Figure 495310DEST_PATH_IMAGE052
Figure 121595DEST_PATH_IMAGE053
wherein:
Figure 416310DEST_PATH_IMAGE082
representing the number of drainage pipes in the urban drainage network to be optimized. And carrying out optimization solution on the constructed urban drainage network optimization model by using a gradient optimization algorithm, wherein a rank 2 quasi Newton method based on a non-precise Armijo criterion is used as the gradient optimization method. According to the scheme, the pressure difference between two ends of different drainage pipelines is determined through the flow to be discharged based on drainage pipelines in different areas, if the pressure difference between two ends of the drainage pipelines is too large, the running pressure of an urban drainage network is increased equivalently, the drainage pipelines are likely to be damaged, the frequency of urban inland inundation disasters is obviously improved, a drainage network optimization design model is built based on the pressure difference between two ends of the drainage pipelines, the model is solved by utilizing a gradient optimization algorithm, the urban drainage network optimization design scheme is obtained, namely, the urban drainage network optimization design scheme is obtained through the steps ofAnd a plurality of drainage pipelines at specific positions are additionally arranged, so that the running pressure of an urban drainage network is reduced, and the frequency of urban waterlogging disasters is reduced.
Drawings
Fig. 1 is a schematic flow chart of a method for optimally designing an urban drainage network according to an embodiment of the present invention;
fig. 2 is a functional block diagram of an apparatus for optimally designing a municipal drainage network according to an embodiment of the present invention.
Detailed Description
It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
The embodiment of the application provides an optimal design method for an urban drainage network. The execution subject of the urban drainage network optimization design method includes, but is not limited to, at least one of electronic devices such as a server and a terminal, which can be configured to execute the method provided by the embodiment of the present application. In other words, the urban drainage network optimization design method may be executed by software or hardware installed in a terminal device or a server device, and the software may be a block chain platform.
Example 1:
s1: and collecting urban drainage network data to form an urban drainage network structure, correcting the urban drainage network structure, and delaying the connection of a network topological relation.
The step S1 of collecting urban drainage network data to form an urban drainage network structure comprises the following steps:
collecting urban drainage network data, wherein the urban drainage network data comprise position parameters of urban drainage pipelines, pipeline lengths and pipeline wall roughness coefficients, constructing an urban drainage network structure based on the urban drainage network data, the urban drainage network structure comprises position distribution of the urban drainage pipelines and pipeline parameters, the pipeline parameters are the pipeline lengths and the pipeline wall roughness coefficients, and connecting pipe orifice positions and pipe bottom positions of the same pipeline by straight lines to form position distribution of the pipeline; the drainage pipeline position parameters comprise pipe orifice positions and pipe bottom positions of urban drainage pipelines, the urban drainage pipelines are all straight lines, and the pipe orifice positions are higher than the pipe bottom positions;
the urban drainage network is provided with n drainage pipelines in total, and the data set of the n drainage pipelines is as follows:
Figure 719115DEST_PATH_IMAGE001
wherein:
Figure 217224DEST_PATH_IMAGE083
indicating any of the second in a municipal drainage network
Figure 630888DEST_PATH_IMAGE084
A water draining pipeline is arranged on the top of the water tank,
Figure 514530DEST_PATH_IMAGE085
is shown as
Figure 988237DEST_PATH_IMAGE084
The position coordinates of the pipe orifice of the strip drainage pipeline,
Figure 708062DEST_PATH_IMAGE086
is shown as
Figure 659837DEST_PATH_IMAGE084
The pipe bottom position coordinates of the strip drainage pipeline are the central position coordinates of the pipe opening or the pipe bottom area;
in the embodiment of the invention, the city is constructed into a three-dimensional coordinate system, the center of the surface of the city is the coordinate origin of the three-dimensional coordinate system, the east-west direction is the X-axis direction, the south-north direction is the Y-axis direction, and the direction from the ground bottom to the surface of the ground is the Z-axis direction;
Figure 397986DEST_PATH_IMAGE006
is shown as
Figure 42594DEST_PATH_IMAGE003
The length of the pipeline of the strip drainage pipeline,
Figure 249716DEST_PATH_IMAGE007
is shown as
Figure 5182DEST_PATH_IMAGE003
The pipe wall roughness coefficient of the strip drainage pipe.
The step S1 of correcting the urban drainage network structure includes:
the correction process of the urban drainage network structure comprises the following steps:
s11: deleting the position coordinates of the isolated pipe orifice or the position coordinates of the pipe bottom in the constructed urban drainage network structure;
s12: selecting a drainage pipeline with a pipe orifice position coordinate on the ground surface for inquiry;
s13: if the pipe orifice position coordinate of only one drainage pipeline is connected with the pipe bottom position coordinate of the selected drainage pipeline, connecting the searched drainage pipeline and the selected drainage pipeline, taking the searched drainage pipeline as the selected pipeline, and repeating the step S13 until the drainage pipeline cannot be searched;
if the distance between the pipe orifice position coordinates of the inquired drainage pipeline and the selected drainage pipeline and the pipe bottom position coordinates is smaller than a distance threshold value, judging that the two position coordinates are connected;
s14: if the pipe orifice position coordinates of a plurality of drainage pipelines are inquired to be connected with the pipe bottom position coordinates of the selected drainage pipeline, respectively connecting the inquired drainage pipeline and the selected drainage pipeline, respectively taking the inquired drainage pipeline as the selected pipeline, and returning to the step S13 until the drainage pipeline cannot be inquired;
s15: deleting the drainage pipeline of which the position coordinate of the pipe orifice is not positioned on the ground surface and the position coordinate of the pipe orifice is not connected with the position coordinate of the pipe bottom of other drainage pipelines in the urban drainage network structure;
based on the correction process of the urban drainage network structure
Figure 863417DEST_PATH_IMAGE087
Modified urban drainage network of drainage pipelinesIn the structure of the utility model, the utility model has the advantages of simple structure,
Figure 678926DEST_PATH_IMAGE087
indicating the number of modified urban sewerage pipelines.
S2: and constructing an urban drainage network model based on the modified urban drainage network structure, wherein the urban drainage network model comprises the flow of different urban drainage pipelines and the pressure difference between two ends of the pipelines, and screening parameters in the model by using a modified Morris screening method to form an example of the urban drainage network model.
In the step S2, based on the modified urban drainage network structure, an urban drainage network model is constructed, including:
dividing the urban surface into a plurality of different sub-areas based on the modified urban drainage network structure, wherein each sub-area is provided with a drainage pipeline with a pipe orifice position coordinate positioned on the ground surface, the number of the divided sub-areas is m, and the jth sub-area is
Figure 107765DEST_PATH_IMAGE088
，
Figure 666922DEST_PATH_IMAGE088
The position coordinate of the middle pipe orifice is positioned on the drainage pipeline on the ground surface
Figure 645242DEST_PATH_IMAGE089
Said sub-region
Figure 366073DEST_PATH_IMAGE088
The drainage network structure of (1) is:
Figure 813366DEST_PATH_IMAGE013
wherein:
Figure 645056DEST_PATH_IMAGE090
indicating water flow from the drain line
Figure 477883DEST_PATH_IMAGE089
Starting from, can pass through
Figure 635195DEST_PATH_IMAGE091
A strip drainage pipeline;
the drainage pipeline
Figure 38626DEST_PATH_IMAGE089
The deepest drainage pipe set capable of being connected in the urban drainage network structure is
Figure 939585DEST_PATH_IMAGE016
The deepest drainage pipeline is a pipeline with the position coordinate of the bottom of the pipeline not connected with other drainage pipelines, and water on the ground surface can pass through the drainage pipelines
Figure 626919DEST_PATH_IMAGE089
Arrive at
Figure 955132DEST_PATH_IMAGE016
The medium and any drainage pipelines and the deepest drainage pipeline are finally communicated with a sewage treatment plant;
the urban drainage network model comprises drainage network models of different sub-areas, and the sub-areas
Figure 845859DEST_PATH_IMAGE088
The drainage network model is as follows:
Figure 284930DEST_PATH_IMAGE017
Figure 826770DEST_PATH_IMAGE018
Figure 325884DEST_PATH_IMAGE019
Figure 703907DEST_PATH_IMAGE020
wherein:
Figure 946670DEST_PATH_IMAGE092
representing an area per unit time
Figure 608595DEST_PATH_IMAGE088
The amount of rainfall of;
Figure 278611DEST_PATH_IMAGE093
representing an area per unit time
Figure 878351DEST_PATH_IMAGE088
The average infiltration capacity of the soil in the middle permeable area;
Figure 924804DEST_PATH_IMAGE094
representing an area per unit time
Figure 441236DEST_PATH_IMAGE088
Maximum water storage capacity of the middle impervious area;
Figure 282153DEST_PATH_IMAGE095
representing an area per unit time
Figure 369189DEST_PATH_IMAGE088
The rainfall evaporation capacity of (d);
Figure 953754DEST_PATH_IMAGE096
representing an area per unit time
Figure 324693DEST_PATH_IMAGE088
The precipitation yield of, i.e. said drainage pipe
Figure 336511DEST_PATH_IMAGE089
The flow rate of the discharged water;
Figure 176422DEST_PATH_IMAGE097
indicating drainage pipe
Figure 299099DEST_PATH_IMAGE089
The pressure difference between the two ends of the pipeline,
Figure 55702DEST_PATH_IMAGE098
which is indicative of the density of the water,
Figure 972843DEST_PATH_IMAGE099
which represents the acceleration of the force of gravity,
Figure 34471DEST_PATH_IMAGE100
indicating precipitation flowing into drainage pipe from surface
Figure 960839DEST_PATH_IMAGE089
The length of the distance the orifice is to flow through,
Figure 571948DEST_PATH_IMAGE101
indicating drainage pipe
Figure 659990DEST_PATH_IMAGE089
The pipe wall roughness coefficient of (a);
Figure 208914DEST_PATH_IMAGE102
connecting drainage pipeline
Figure 938973DEST_PATH_IMAGE103
All upstream drainage pipes of (2) are to the drainage pipe in unit time
Figure 404589DEST_PATH_IMAGE103
The total flow of discharged water is that the upstream water discharge pipeline is positioned higher than the water discharge pipeline
Figure 679844DEST_PATH_IMAGE103
Can supply and drain water
Figure 965332DEST_PATH_IMAGE103
A drainage pipeline for draining sewage;
Figure 764661DEST_PATH_IMAGE104
indicating area
Figure 819204DEST_PATH_IMAGE088
The ratio of the road area of (a),
Figure 519221DEST_PATH_IMAGE105
indicating area
Figure 557584DEST_PATH_IMAGE088
The area of the forest land is in proportion,
Figure 645757DEST_PATH_IMAGE106
indicating area
Figure 820386DEST_PATH_IMAGE088
The grass area ratio of (2).
In the step S2, parameters in the model are screened by using a modified Morris screening method to determine an urban drainage network model, including:
the model parameters comprise the average infiltration amount and the maximum water storage amount of the soil in different subregions, the rainfall evaporation amount in unit time and the rainfall amount;
screening parameters in the model by using a modified Morris screening method, wherein the model parameter screening process of the modified Morris screening method comprises the following steps:
s21: randomly selecting one parameter from model parameters to modify for many times, wherein the modification mode of the model parameters is to collect rainfall conditions of the same city at different times as the model parameters;
s22: inputting the parameter result after each modification into the urban drainage network model to obtain a plurality of groups of surface drainage pipeline flow results in the urban drainage network model;
s23: calculating sensitivity of modified parameters
Figure 217869DEST_PATH_IMAGE107
：
Figure 228682DEST_PATH_IMAGE041
Wherein:
the number of times of parameter modification is K-1;
Figure 104234DEST_PATH_IMAGE108
representing a drainage pipeline flow result obtained based on the kth modified parameter;
Figure 133370DEST_PATH_IMAGE109
representing the change rate of the k-th modified parameter compared with the initial parameter value;
s24: and repeating the steps to obtain the sensitivity of different model parameters, and deleting the parameters with the sensitivity absolute value lower than 0.02 in the model to obtain the urban drainage network model suitable for the selected city.
S3: and constructing an urban drainage network optimization model, and determining an urban drainage network design constraint standard expression, wherein the solution result of the urban drainage network optimization model is the position of a newly added drainage pipeline in the urban drainage network.
And the step S3 of constructing an urban drainage network optimization model comprises the following steps:
the constructed urban drainage network optimization model B is as follows:
Figure 186908DEST_PATH_IMAGE044
wherein:
Figure 668705DEST_PATH_IMAGE110
showing the tube bottom coordinate position of the added h-th drainage pipeline,
Figure 347948DEST_PATH_IMAGE111
indicating the coordinate location of any sewage treatment plant,
Figure 231590DEST_PATH_IMAGE112
the distance between the bottom of the h-th drainage pipeline and the nearest sewage treatment plant is shown;
Figure 456029DEST_PATH_IMAGE113
showing the original drainage flow of the added h-th drainage pipeline at the pipe orifice position, wherein the original drainage flow shows the drainage flow when the drainage pipeline is not added;
Figure 425122DEST_PATH_IMAGE114
showing the current drainage flow of the pipe orifice position after the h-th drainage pipeline is added;
h represents the number of newly added drainage pipelines in the urban drainage network optimization design scheme;
the solving result of the urban drainage network optimization model is the position of a newly added drainage pipeline in the urban drainage network;
determining a standard expression of urban drainage network design constraints, wherein the drainage network design constraints are as follows:
Figure 642477DEST_PATH_IMAGE052
Figure 380626DEST_PATH_IMAGE053
wherein:
Figure 775966DEST_PATH_IMAGE087
representing the number of drainage pipes in the urban drainage network to be optimized.
S4: and optimizing and solving the constructed urban drainage network optimization model by using a gradient optimization algorithm, optimizing the urban drainage network structure according to the solving result, placing the optimized drainage network in the urban drainage network model for operation, implementing if the effect reaches a specified threshold value, and otherwise, solving again, wherein a rank 2 quasi Newton method based on a non-precise Armijo criterion is the gradient optimization method.
And in the step S4, the constructed urban drainage network optimization model is optimized and solved by using a gradient optimization algorithm, and the method comprises the following steps:
the solving process of the urban drainage network optimization model comprises the following steps:
s41: constructing a penalty function form of an urban drainage network optimization model and a constraint specification expression:
Figure 232355DEST_PATH_IMAGE054
Figure 253401DEST_PATH_IMAGE055
wherein:
Figure 846056DEST_PATH_IMAGE115
the position distribution of newly added H drainage pipelines in the urban drainage network optimization design scheme is shown, and in the embodiment of the invention, H is set to be
Figure 412298DEST_PATH_IMAGE087
；
Figure 90404DEST_PATH_IMAGE116
Set it to 10 for the penalty function coefficient;
s42: position distribution of initially generated H newly-added drainage pipelines
Figure 915140DEST_PATH_IMAGE058
And setting an allowable error of
Figure 378614DEST_PATH_IMAGE059
The initial positive definite matrix is
Figure 365024DEST_PATH_IMAGE060
，
Figure 530427DEST_PATH_IMAGE061
The method is characterized in that the method is an identity matrix, the current iteration number of the algorithm is w, and the initial value of w is 0;
s43: calculating a penalty function in
Figure 627696DEST_PATH_IMAGE062
Gradient of (2)
Figure 211255DEST_PATH_IMAGE117
If, if
Figure 368567DEST_PATH_IMAGE064
Then the position distribution obtained by the w-th iteration
Figure 21265DEST_PATH_IMAGE062
That is, the optimization solution result is obtained, otherwise, the step S44 is turned to;
s44: computing
Figure 922225DEST_PATH_IMAGE118
；
S45: updated by the Armijo criterion
Figure 360290DEST_PATH_IMAGE119
：
Figure 688504DEST_PATH_IMAGE067
Figure 828498DEST_PATH_IMAGE068
Wherein:
Figure 267570DEST_PATH_IMAGE120
is composed of
Figure 560142DEST_PATH_IMAGE121
The update step size of (c);
Figure 59256DEST_PATH_IMAGE122
to update the parameter, it is set to 0.5;
Figure 686547DEST_PATH_IMAGE123
to satisfy
Figure 929309DEST_PATH_IMAGE124
The minimum trace distribution matrix of (a), T represents transposition;
s46: is updated to
Figure 76388DEST_PATH_IMAGE125
：
Figure 746404DEST_PATH_IMAGE075
S47: order to
Figure 860990DEST_PATH_IMAGE126
The process returns to step S43.
In the step S4, the optimized drainage network is placed in an urban drainage network model to operate, if the effect reaches a specified threshold value, the operation is implemented, and if not, the solution is carried out again, and the method comprises the following steps: the optimized urban drainage network topological structure is placed in an urban drainage network model, precipitation data of different seasons of a city are collected and input into the model, the model outputs the flow of each drainage pipeline in the urban drainage network and the pressure difference between two ends of each pipeline, if the pressure difference between two ends of each drainage pipeline is smaller than a specified threshold value, an optimized design scheme of the urban drainage network can be implemented, namely the load pressure of the drainage network is reduced by adding a plurality of drainage pipelines, and otherwise, a gradient optimization algorithm is reused for solving.
Example 2:
fig. 2 is a functional block diagram of an apparatus for optimally designing a municipal drainage network according to an embodiment of the present invention, which can implement the method for optimally designing a municipal drainage network according to embodiment 1.
The urban drainage network optimization design device 100 can be installed in electronic equipment. According to the realized functions, the urban drainage network optimization design device can comprise a drainage network construction device 101, an optimization model construction device 102 and a drainage network optimization module 103. The module of the present invention, which may also be referred to as a unit, refers to a series of computer program segments that can be executed by a processor of an electronic device and that can perform a fixed function, and that are stored in a memory of the electronic device.
The drainage network construction device 101 is used for collecting urban drainage network data to form an urban drainage network structure, correcting the urban drainage network structure, constructing an urban drainage network model based on the corrected urban drainage network structure, and screening parameters in the model by using a corrected Morris screening method to form the urban drainage network model;
the optimization model construction device 102 is used for constructing an urban drainage network optimization model and determining an urban drainage network design constraint specification expression;
and the drainage network optimization module 103 is used for optimizing and solving the constructed urban drainage network optimization model by using a gradient optimization algorithm, optimizing the urban drainage network structure according to the solution result, placing the optimized drainage network in the urban drainage network model for operation, implementing the corresponding urban drainage network optimization design scheme if the effect reaches a specified threshold value, and otherwise, solving again.
In detail, when the modules in the device 100 for optimally designing a municipal drainage network according to the embodiment of the present invention are used, the same technical means as the method for optimally designing a municipal drainage network described in fig. 1 are adopted, and the same technical effects can be produced, which is not described herein again.
It should be noted that the above-mentioned numbers of the embodiments of the present invention are merely for description, and do not represent the merits of the embodiments. And the terms "comprises," "comprising," or any other variation thereof, herein are intended to cover a non-exclusive inclusion, such that a process, apparatus, article, or method that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, apparatus, article, or method. Without further limitation, an element defined by the phrase "comprising a … …" does not exclude the presence of another identical element in a process, apparatus, article, or method that comprises the element.
Through the above description of the embodiments, those skilled in the art will clearly understand that the method of the above embodiments can be implemented by software plus a necessary general hardware platform, and certainly can also be implemented by hardware, but in many cases, the former is a better implementation manner. Based on such understanding, the technical solutions of the present invention may be embodied in the form of a software product, which is stored in a storage medium (e.g., ROM/RAM, magnetic disk, optical disk) and includes instructions for enabling a terminal device (e.g., a mobile phone, a computer, a server, or a network device) to execute the method according to the embodiments of the present invention.
The above description is only a preferred embodiment of the present invention, and not intended to limit the scope of the present invention, and all modifications of equivalent structures and equivalent processes, which are made by using the contents of the present specification and the accompanying drawings, or directly or indirectly applied to other related technical fields, are included in the scope of the present invention.

Claims (7)

1. A method for optimally designing an urban drainage network is characterized by comprising the following steps:
s1: collecting urban drainage network data to form an urban drainage network structure, and correcting the urban drainage network structure;
s2: constructing an urban drainage network model based on the modified urban drainage network structure, wherein the urban drainage network model comprises the flow of different urban drainage pipelines and the pressure difference between two ends of the pipelines, and screening parameters in the model by using a modified Morris screening method to form the urban drainage network model;
s3: constructing an urban drainage network optimization model, and determining an urban drainage network design constraint specification expression, wherein the solution result of the urban drainage network optimization model is the position of a newly added drainage pipeline in the urban drainage network;
s4: optimizing and solving the constructed urban drainage network optimization model by using a gradient optimization algorithm, optimizing the urban drainage network structure according to the solving result, placing the optimized drainage network in the urban drainage network model for operation, implementing if the effect reaches a specified threshold value, otherwise, solving again, wherein a rank 2 quasi-Newton method based on a non-precise Armijo criterion is used as the gradient optimization method and comprises the following steps:
the solving process of the urban drainage network optimization model comprises the following steps:
s41: constructing a penalty function form of an urban drainage network optimization model and a constraint specification expression:
Figure 525363DEST_PATH_IMAGE001
Figure 526949DEST_PATH_IMAGE002
wherein:
Figure 367866DEST_PATH_IMAGE003
indicating new additions in the optimum design of urban drainage network
Figure 969748DEST_PATH_IMAGE004
The positions of the strip drainage pipelines are distributed;
Figure 819893DEST_PATH_IMAGE005
set it to 10 for the penalty function coefficient;
s42: initialization generation
Figure 941564DEST_PATH_IMAGE004
Position distribution of newly-added drainage pipeline
Figure 953382DEST_PATH_IMAGE006
And setting an allowable error of
Figure 42561DEST_PATH_IMAGE007
The initial positive definite matrix is
Figure 193268DEST_PATH_IMAGE008
，
Figure 684292DEST_PATH_IMAGE008
The method is characterized in that the method is an identity matrix, the current iteration number of the algorithm is w, and the initial value of w is 0;
s43: calculating a penalty function in
Figure 132591DEST_PATH_IMAGE009
Gradient of (2)
Figure 443487DEST_PATH_IMAGE010
If at all
Figure 120587DEST_PATH_IMAGE011
Then the position distribution obtained by the w-th iteration
Figure 466117DEST_PATH_IMAGE009
That is, the optimization solution result is obtained, otherwise, the step S44 is turned to;
s44: computing
Figure 819738DEST_PATH_IMAGE012
；
S45: updated based on Armijo criteria
Figure 634242DEST_PATH_IMAGE013
：
Figure 98721DEST_PATH_IMAGE014
Figure 564337DEST_PATH_IMAGE015
Wherein:
Figure 88860DEST_PATH_IMAGE016
is composed of
Figure 125080DEST_PATH_IMAGE017
The update step length of (2);
Figure 393250DEST_PATH_IMAGE018
to update the parameter, it is set to 0.5;
Figure 447794DEST_PATH_IMAGE019
to satisfy
Figure 408797DEST_PATH_IMAGE020
The distribution matrix of the smallest trace of (a),
Figure 932313DEST_PATH_IMAGE021
representing a transpose;
s46: is updated to
Figure 738595DEST_PATH_IMAGE022
：
Figure 913224DEST_PATH_IMAGE023
S47: order to
Figure 779549DEST_PATH_IMAGE024
The process returns to step S43.
2. The method according to claim 1, wherein the step S1 of collecting the municipal drainage network data to form a municipal drainage network structure comprises:
collecting urban drainage network data, wherein the urban drainage network data comprise position parameters of urban drainage pipelines, pipeline lengths and pipeline wall roughness coefficients, constructing an urban drainage network structure based on the urban drainage network data, the urban drainage network structure comprises position distribution of the urban drainage pipelines and pipeline parameters, the pipeline parameters are the pipeline lengths and the pipeline wall roughness coefficients, and connecting pipe orifice positions and pipe bottom positions of the same pipeline by straight lines to form position distribution of the pipeline;
the urban drainage network is provided with n drainage pipelines in total, and the data set of the n drainage pipelines is as follows:
Figure 790362DEST_PATH_IMAGE025
wherein:
Figure 400335DEST_PATH_IMAGE026
indicating any of the second in a municipal drainage network
Figure 429470DEST_PATH_IMAGE027
A water draining pipeline is arranged on the top of the water tank,
Figure 732276DEST_PATH_IMAGE028
is shown as
Figure 964805DEST_PATH_IMAGE027
The position coordinates of the pipe orifice of the strip drainage pipeline,
Figure 378469DEST_PATH_IMAGE029
denotes the first
Figure 262111DEST_PATH_IMAGE030
The pipe bottom position coordinates of the strip drainage pipeline are the central position coordinates of the pipe opening or the pipe bottom area;
Figure 735818DEST_PATH_IMAGE031
is shown as
Figure 721223DEST_PATH_IMAGE032
The length of the pipeline of the strip drainage pipeline,
Figure 672998DEST_PATH_IMAGE033
is shown as
Figure 676726DEST_PATH_IMAGE030
The pipe wall roughness coefficient of the strip drainage pipe.
3. The method according to claim 2, wherein the step S1 of modifying the urban drainage network structure comprises:
the correction process of the urban drainage network structure comprises the following steps:
s11: deleting the position coordinates of the isolated pipe orifice or the pipe bottom in the constructed urban drainage network structure;
s12: selecting a drainage pipeline with a pipe orifice position coordinate on the ground surface for inquiry;
s13: if the pipe orifice position coordinate of only one drainage pipeline is connected with the pipe bottom position coordinate of the selected drainage pipeline, connecting the searched drainage pipeline and the selected drainage pipeline, taking the searched drainage pipeline as the selected pipeline, and repeating the step S13 until the drainage pipeline cannot be searched;
if the distance between the pipe orifice position coordinates of the inquired drainage pipeline and the selected drainage pipeline and the pipe bottom position coordinates is smaller than a distance threshold value, judging that the two position coordinates are connected;
s14: if the pipe orifice position coordinates of a plurality of drainage pipelines are inquired to be connected with the pipe bottom position coordinates of the selected drainage pipeline, respectively connecting the inquired drainage pipeline and the selected drainage pipeline, respectively taking the inquired drainage pipeline as the selected pipeline, and returning to the step S13 until the drainage pipeline cannot be inquired;
s15: deleting the drainage pipeline of which the pipe orifice position coordinate is not positioned on the ground surface and is not connected with the pipe bottom position coordinate of other drainage pipelines in the urban drainage network structure;
based on the correction process of the urban drainage network structure
Figure 321334DEST_PATH_IMAGE034
The modified urban drainage network structure of the drainage pipeline,
Figure 794035DEST_PATH_IMAGE034
indicating the number of modified urban sewerage pipelines.
4. The method according to claim 3, wherein the step S2 of constructing the urban drainage network model based on the modified urban drainage network structure comprises:
dividing the urban surface into a plurality of different sub-areas based on the modified urban drainage network structure, wherein each sub-area is provided with a drainage pipeline with a pipe orifice position coordinate positioned on the ground surface, the number of the divided sub-areas is m, and the number of the first sub-areas is
Figure 549501DEST_PATH_IMAGE035
Sub-regions are
Figure 142157DEST_PATH_IMAGE036
，
Figure 708398DEST_PATH_IMAGE036
The position coordinate of the middle pipe orifice is positioned on the drainage pipeline on the ground surface
Figure 652084DEST_PATH_IMAGE037
Said sub-region
Figure 476820DEST_PATH_IMAGE036
The drainage network structure of (1) is as follows:
Figure 923982DEST_PATH_IMAGE038
wherein:
Figure 661125DEST_PATH_IMAGE039
indicating water flow from the drain line
Figure 92106DEST_PATH_IMAGE037
Starting from, can pass through
Figure 189375DEST_PATH_IMAGE040
A strip drainage pipeline;
the drainage pipeline
Figure 22202DEST_PATH_IMAGE037
The deepest drainage pipe set capable of being connected in the urban drainage network structure is
Figure 664667DEST_PATH_IMAGE041
The deepest drainage pipeline is a pipeline with the position coordinate of the bottom of the pipeline not connected with other drainage pipelines, and water on the ground surface can pass through the drainage pipelines
Figure 582945DEST_PATH_IMAGE037
Arrive at
Figure 483905DEST_PATH_IMAGE041
The middle and any drainage pipelines and the deepest drainage pipeline are finally communicated with a sewage treatment plant;
the urban drainage network model comprises drainage network models of different sub-areas, and the sub-areas
Figure 905659DEST_PATH_IMAGE036
The drainage network model is as follows:
Figure 984604DEST_PATH_IMAGE042
Figure 124599DEST_PATH_IMAGE043
Figure 829249DEST_PATH_IMAGE044
Figure 371089DEST_PATH_IMAGE045
wherein:
Figure 620936DEST_PATH_IMAGE046
representing an area per unit time
Figure 248226DEST_PATH_IMAGE036
The amount of rainfall of;
Figure 490989DEST_PATH_IMAGE047
representing an area per unit time
Figure 887335DEST_PATH_IMAGE036
The average infiltration capacity of the soil in the middle permeable area;
Figure 308083DEST_PATH_IMAGE048
representing an area per unit time
Figure 422670DEST_PATH_IMAGE036
Maximum water storage capacity of the middle impervious area;
Figure 203544DEST_PATH_IMAGE049
representing an area per unit time
Figure 719976DEST_PATH_IMAGE036
The rainfall evaporation amount of (2);
Figure 311626DEST_PATH_IMAGE050
representing an area per unit time
Figure 913508DEST_PATH_IMAGE036
The precipitation yield of, i.e. said drainage pipe
Figure 498073DEST_PATH_IMAGE037
The flow rate of the discharged water;
Figure 869012DEST_PATH_IMAGE051
indicating drainage pipe
Figure 885423DEST_PATH_IMAGE037
The pressure difference between the two ends of the pipeline,
Figure 709023DEST_PATH_IMAGE052
which is indicative of the density of the water,
Figure 97279DEST_PATH_IMAGE053
which represents the acceleration of the force of gravity,
Figure 57144DEST_PATH_IMAGE054
indicating precipitation flowing into drainage pipe from surface
Figure 256176DEST_PATH_IMAGE037
The length of the distance the orifice is to flow through,
Figure 832650DEST_PATH_IMAGE055
indicating drainage pipe
Figure 759018DEST_PATH_IMAGE037
The pipe wall roughness coefficient of (a);
Figure 589702DEST_PATH_IMAGE056
connecting drainage pipeline
Figure 208902DEST_PATH_IMAGE057
All upstream drainage pipes of (2) are to the drainage pipe in unit time
Figure 7094DEST_PATH_IMAGE058
The total flow of discharged water is that the upstream water discharge pipeline is positioned higher than the water discharge pipeline
Figure 737153DEST_PATH_IMAGE059
Can supply water to the drainage pipeline
Figure 422343DEST_PATH_IMAGE060
A drain pipeline for discharging sewage;
Figure 212444DEST_PATH_IMAGE061
indicating area
Figure 497932DEST_PATH_IMAGE036
The ratio of the road area of (a),
Figure 766103DEST_PATH_IMAGE062
indicating area
Figure 571379DEST_PATH_IMAGE036
The area of the forest land is in proportion,
Figure 532381DEST_PATH_IMAGE063
indicating area
Figure 305165DEST_PATH_IMAGE036
The grass area ratio of (2).
5. The method of claim 4, wherein the step S2 of determining the model of the municipal drainage network by screening parameters in the model using a modified Morris screening method comprises:
the model parameters comprise the average infiltration amount and the maximum water storage amount of the soil in different subregions, the rainfall evaporation amount in unit time and the rainfall amount;
screening parameters in the model by using a modified Morris screening method, wherein the model parameter screening process of the modified Morris screening method comprises the following steps:
s21: randomly selecting one parameter from model parameters to modify for many times, wherein the modification mode of the model parameters is to collect rainfall conditions of the same city at different times as the model parameters;
s22: inputting the parameter result after each modification into the urban drainage network model to obtain a plurality of groups of surface drainage pipeline flow results in the urban drainage network model;
s23: calculating sensitivity of modified parameters
Figure 377027DEST_PATH_IMAGE064
：
Figure 302388DEST_PATH_IMAGE065
Wherein:
the number of times of modifying the parameters is K-1;
Figure 434292DEST_PATH_IMAGE066
representing a drainage pipeline flow result obtained based on the kth modified parameter;
Figure 694372DEST_PATH_IMAGE067
representing the change rate of the k-th modified parameter compared with the initial parameter value;
s24: and repeating the steps to obtain the sensitivity of different model parameters, and deleting the parameters with the sensitivity absolute value lower than 0.02 in the model to obtain the urban drainage network model suitable for the selected city.
6. The method according to claim 1, wherein the step S3 of constructing an urban drainage network optimization model comprises:
the constructed urban drainage network optimization model B is as follows:
Figure 55078DEST_PATH_IMAGE068
wherein:
Figure 84214DEST_PATH_IMAGE069
indicates the added first
Figure 121440DEST_PATH_IMAGE070
The coordinate position of the bottom of the strip drainage pipeline,
Figure 868816DEST_PATH_IMAGE071
indicating the coordinate location of any sewage treatment plant,
Figure 33212DEST_PATH_IMAGE072
the distance between the bottom of the h-th drainage pipeline and the nearest sewage treatment plant is shown;
Figure 916855DEST_PATH_IMAGE073
indicates the increased second
Figure 390561DEST_PATH_IMAGE074
The original drainage flow of the opening position of the strip drainage pipeline represents the drainage flow when the drainage pipeline is not increased;
Figure 359654DEST_PATH_IMAGE075
showing the current drainage flow of the pipe orifice after the h-th drainage pipeline is added;
h represents the number of newly added drainage pipelines in the urban drainage network optimization design scheme;
the solving result of the urban drainage network optimization model is the position of a newly added drainage pipeline in the urban drainage network;
determining a standard expression of urban drainage network design constraints, wherein the drainage network design constraints are as follows:
Figure 62162DEST_PATH_IMAGE076
Figure 331469DEST_PATH_IMAGE077
wherein:
Figure 710498DEST_PATH_IMAGE078
representing the number of drainage pipes in the urban drainage network to be optimized.
7. The method of claim 1, wherein the step S4 of putting the optimized drainage network into the urban drainage network model for operation is performed if the effect reaches a specified threshold, and otherwise, the step of solving again includes:
the optimized urban drainage network topological structure is placed in an urban drainage network model, precipitation data of different seasons of a city are collected and input into the model, the model outputs the flow of each drainage pipeline in the urban drainage network and the pressure difference between two ends of each pipeline, if the pressure difference between two ends of each drainage pipeline is smaller than a specified threshold value, an optimized design scheme of the urban drainage network can be implemented, namely the load pressure of the drainage network is reduced by adding a plurality of drainage pipelines, and otherwise, a gradient optimization algorithm is reused for solving.
CN202211408100.9A 2022-11-10 2022-11-10 Urban drainage network optimization design method Active CN115510600B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211408100.9A CN115510600B (en) 2022-11-10 2022-11-10 Urban drainage network optimization design method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211408100.9A CN115510600B (en) 2022-11-10 2022-11-10 Urban drainage network optimization design method

Publications (2)

Publication Number Publication Date
CN115510600A true CN115510600A (en) 2022-12-23
CN115510600B CN115510600B (en) 2023-08-01

Family

ID=84513954

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211408100.9A Active CN115510600B (en) 2022-11-10 2022-11-10 Urban drainage network optimization design method

Country Status (1)

Country Link
CN (1) CN115510600B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117647963A (en) * 2024-01-29 2024-03-05 四川速聚智联科技有限公司 Intelligent liquid discharge control method and system for natural gas square well pool
CN117892560A (en) * 2024-03-14 2024-04-16 长沙市海图科技有限公司 High-geographic-precision intelligent city safety drainage pipe network simulation and prediction method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090106173A1 (en) * 2007-10-17 2009-04-23 Microsoft Corporation Limited-memory quasi-newton optimization algorithm for l1-regularized objectives
CN109657841A (en) * 2018-11-27 2019-04-19 杭州师范大学 A kind of urban rainstorm waterlogging depth of accumulated water extracting method
CN114741782A (en) * 2022-03-28 2022-07-12 大连理工大学 Conical shell reinforcement structure optimization method and device, computer and storage medium
CN115186515A (en) * 2022-08-26 2022-10-14 中国长江三峡集团有限公司 Method and device for predicting residual life of concrete drainage pipeline

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090106173A1 (en) * 2007-10-17 2009-04-23 Microsoft Corporation Limited-memory quasi-newton optimization algorithm for l1-regularized objectives
CN109657841A (en) * 2018-11-27 2019-04-19 杭州师范大学 A kind of urban rainstorm waterlogging depth of accumulated water extracting method
CN114741782A (en) * 2022-03-28 2022-07-12 大连理工大学 Conical shell reinforcement structure optimization method and device, computer and storage medium
CN115186515A (en) * 2022-08-26 2022-10-14 中国长江三峡集团有限公司 Method and device for predicting residual life of concrete drainage pipeline

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘建美等: ""基于 BFGS 公式的改进截断拟牛顿法在 随机用户均衡问题上的应用"", 《山东大学学报（工学版）》, vol. 48, no. 1 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117647963A (en) * 2024-01-29 2024-03-05 四川速聚智联科技有限公司 Intelligent liquid discharge control method and system for natural gas square well pool
CN117647963B (en) * 2024-01-29 2024-04-09 四川速聚智联科技有限公司 Intelligent liquid discharge control method and system for natural gas square well pool
CN117892560A (en) * 2024-03-14 2024-04-16 长沙市海图科技有限公司 High-geographic-precision intelligent city safety drainage pipe network simulation and prediction method
CN117892560B (en) * 2024-03-14 2024-06-11 长沙市海图科技有限公司 High-geographic-precision intelligent city safety drainage pipe network simulation and prediction method

Also Published As

Publication number Publication date
CN115510600B (en) 2023-08-01

Similar Documents

Publication Publication Date Title
CN115510600A (en) Urban drainage network optimization design method
CN109657841B (en) Deep extraction method for urban rainstorm waterlogging
CN108182543A (en) One kind becomes more meticulous grid waterlogging water logging forecasting procedure
Ren et al. Evaluating the stormwater management model to improve urban water allocation system in drought conditions
CN108388714A (en) The plain river network city flood simulation method of basin water system and urban pipe network coupling
Yang et al. Linking a storm water management model to a novel two-dimensional model for urban pluvial flood modeling
CN111199347A (en) Drainage basin pollution control unit partitioning method
Hassan et al. Development of sediment load estimation models by using artificial neural networking techniques
Yu et al. A new LID spatial allocation optimization system at neighborhood scale: Integrated SWMM with PICEA-g using MATLAB as the platform
Ciou et al. Optimization model for BMP placement in a reservoir watershed
Pan et al. Hydrologically enhanced distributed urban drainage model and its application in Beijing city
CN112699610A (en) Sponge city optimization design method based on high-dimensional multi-objective evolutionary algorithm
CN112242003B (en) City sub-catchment area division method considering land type and flow direction
CN115758886A (en) Optimal distribution method of storage tanks based on rain flood pipe network model and decision tree algorithm
CN106156949A (en) A kind of rainwater based on waterlogging risk assessment is low affects development project method
CN112907047A (en) Sponge city performance evaluation system
Yu et al. New framework for assessing urban stormwater management measures in the context of climate change
CN112699606A (en) Sponge city layout optimization multi-objective decision-making method
He et al. Constrained optimization model of the volume of initial rainwater storage tank based on ANN and PSO
CN117012004B (en) Urban waterlogging risk early warning method and system
CN111368397B (en) Method and device for predicting waterlogging risk
Zhou et al. Establishment and validation of an amended phosphorus index: refined phosphorus loss assessment of an agriculture watershed in Northern China
CN108920571B (en) Indoor hydrological experiment model scale determination method based on distributed hydrological model
CN117150600A (en) Roof water storage and drainage control method for urban waterlogging prevention
Zhang et al. Simulation of urban rainstorm waterlogging and pipeline network drainage process based on SWMM

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