CN108824593B - 一种针对复杂雨水管网的快速多目标工程优化设计方法 - Google Patents

一种针对复杂雨水管网的快速多目标工程优化设计方法 Download PDF

Info

Publication number
CN108824593B
CN108824593B CN201810281014.3A CN201810281014A CN108824593B CN 108824593 B CN108824593 B CN 108824593B CN 201810281014 A CN201810281014 A CN 201810281014A CN 108824593 B CN108824593 B CN 108824593B
Authority
CN
China
Prior art keywords
pipe network
rainwater pipe
rainwater
optimization
formula
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810281014.3A
Other languages
English (en)
Other versions
CN108824593A (zh
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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201810281014.3A priority Critical patent/CN108824593B/zh
Publication of CN108824593A publication Critical patent/CN108824593A/zh
Application granted granted Critical
Publication of CN108824593B publication Critical patent/CN108824593B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • E—FIXED CONSTRUCTIONS
    • E03—WATER SUPPLY; SEWERAGE
    • E03F—SEWERS; CESSPOOLS
    • E03F1/00—Methods, systems, or installations for draining-off sewage or storm water
    • E—FIXED CONSTRUCTIONS
    • E03—WATER SUPPLY; SEWERAGE
    • E03F—SEWERS; CESSPOOLS
    • E03F3/00—Sewer pipe-line systems
    • E03F3/02—Arrangement of sewer pipe-lines or pipe-line systems
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Hydrology & Water Resources (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • Theoretical Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Sewage (AREA)

Abstract

本发明提供一种针对复杂雨水管网的快速多目标工程优化设计方法，包括：选定暴雨强度公式及其参数，给定暴雨重现期，选取典型降雨过程曲线；建立雨水管网水力模型，并将降雨过程曲线导入模型中的降雨模块；确定充满度解集；计算雨水管网的管径，调整计算管径，确定雨水管网初始优化设计方案集，建立雨水管网多目标优化模型，启动优化算法，以快速获得雨水管网帕累托优化解集。本发明具有原创性，是对雨水管网优化设计研究领域的一个重要补充，可为城市内涝防控和海绵城市建设提供了重要的技术支撑，因此具有很好的推广和实际工程应用价值。

Description

一种针对复杂雨水管网的快速多目标工程优化设计方法
技术领域
本发明涉及市政工程和城市雨水管网领域。
背景技术
近年来，全国多处城市出现内涝现象，导致城市交通瘫痪、电力供应中断等问题，给人民日常生产和社会经济造成巨大影响。造成城市内涝的因素很多，如气候变化、城市化下垫面硬化、下游水位顶托等。雨水管网的不合理设计也是城市内涝的一个重要因素，如雨水管道设计标准偏低，往往不能在有效时间内将雨水排出，造成城市内涝灾害。因此，雨水管网的合理设计对于城市内涝的有效防控具有重要意义。
针对雨水管网设计问题，国内工程界普遍采用推理公式法，并基于管道满管流设计，虽然该方法具有简单方便的优点，但未考虑真实降雨过程，导致所设计的管网往往不能满足预定的重现期排水目标。这个问题在气候变化背景下变的更加突出，这是因为暴雨峰值在气候变化驱动下呈现显著上升的趋势。为解决这个问题，有专家学者提出了雨水管网的多目标优化，即在雨水管网设计中不仅考虑经济性目标，也要考虑基于降雨过程曲线的设计充满度优化目标。管网设计的充满度越低，其应对气候变化导致的暴雨不确定性能力越强，城市内涝风险越小，但其管网造价越高，这就需要设计者在多目标优化中根据实际情况选出最合理的设计方案。
尽管多目标优化已开始在城市雨水管网设计中尝试应用，但其普遍存在优化效率低下的问题，尤其是在解决大型复杂雨水管网的设计问题时，优化效率已经是一个重要的瓶颈问题。这是因为考虑降雨过程线需要反复的调用雨水管网水力模型进行目标函数计算，导致优化效率低下。例如，使用智能算法对某城市局部地区雨水管网进行优化设计，该雨水管网包含102个节点、102根管道和1个排水口。如用计算机(2.9GHz Dell PC(InterR))基于SWMM排水管网水力模型迭代100000次目标函数(即调用100000次雨水管网水力模型)，每次运行需要7秒钟，则总共需要194.4小时(约8天)，如此低效的优化计算使算法很难进行实际应用。由此可见，如何提高复杂雨水管网多目标优化设计效率，已成为排水研究领域的关键科学瓶颈问题。
发明内容
本发明所要解决的技术问题：提供一种针对复杂雨水管网的快速多目标工程优化设计方法，该方法结合推理公式法和过程优化法，可大幅度提高优化求解效率。本方法的核心技术方案如下：
(1)根据雨水管网所在地区选定其对应的暴雨强度公式及其参数，依据雨水管网设计要求给定暴雨重现期P，并从降雨历史记录数据选取与之相对应的典型降雨过程曲线。
(2)建立雨水管网水力模型，并将步骤(1)中的降雨过程曲线导入该模型中的降雨模块。
(3)确定充满度解集E，可表达如下:
E＝[e1,e2,...,en]T 式1
式中，ei为管道设计充满度，取值范围为0.60～1.00，i＝1,2,…,n，表示n个不同的设计充满度。
(4)根据步骤(1)的推理公式参数及确定的管网布局，可计算管道设计流量Q，进而通过管道非满流水力计算公式，为每一个设定的充满度ei计算雨水管网中每一根管道的管径dk(ei),k＝1,2,…,K(K为雨水管道总数)，结果可用公式表达如下：
ψ(ei)＝[d1(ei),d2(ei),...,dk(ei)]T 式2
式中，ψ(ei)为设计充满度ei时优化管径设计方案。
(5)调整ψ(ei)的计算管径，可表达如下：
S(ei)＝adjust{ψ(ei)} 式3
由于采用推理公式法计算管径，ψ(ei)中可存在连续的管径值，而实际工程中的管径是非连续的，因此采用adjust函数对ψ(ei)中dk(ei)调整为其最接近的工程管径。
(6)确定雨水管网初始优化设计方案集。循环步骤(4)和(5)直到充满度解集E的所有值完成计算，即可获得在不同设计充满度下的雨水管网优化设计方案集，可表达为：
S＝{S(e1),S(e2),...,S(en)} 式4
(7)建立雨水管网多目标优化模型。本发明方法以雨水管道管径d为决策变量，以管道投资成本最低与管道平均充满度最低为优化目标函数，并选定多目标优化算法求解。
(8)将步骤(6)中初始优化管径解集作为步骤(7)中优化算法的初始种群，启动优化算法，以快速获得雨水管网帕累托优化解集。
为解决大型复杂雨水管网优化设计耗时长的瓶颈问题，本发明首次提出了一种针对复杂雨水管网系统的快速多目标优化设计方法。本发明与现有雨水管网优化设计方法相比具有的优点是：现有雨水管网优化方法是基于算法的随机搜索，计算速度慢，很难应用于大型复杂雨水管网的设计，而在本发明中，通过推理公式法为优化算法提供初始解，作为优化算法的多个目标，通过对多目标的优化，节省了算法的计算过程，有机耦合了推理公式法和优化算法的优点，大幅度提高了算法的优化效率。本发明具有原创性，是对雨水管网优化设计研究领域的一个重要补充，可为城市内涝防控和海绵城市建设提供了重要的技术支撑，因此具有很好的推广和实际工程应用价值。
附图说明
图1是典型降雨过程曲线图。
图2是本发明快速多目标工程优化设计方法总流程图。
图3是某城市雨水管网示意图。
图4是本发明设计方法迭代5000次与传统优化方法迭代150000次的对比结果。
具体实施方式
参见图2，本发明的具体实施步骤如下：
(1)根据雨水管网所在地区选定其对应的暴雨强度公式
Figure BDA0001614635900000031
其中的参数A＝57.694，C＝0.93，b＝31.546，n＝1.008，t＝t1+t2，t1为地表集水时间，取t1＝6min，t2为管内雨水流行时间。依据设计要求给定暴雨重现期P，在本案例中P＝5，并从降雨历史记录数据选取与之相对应的典型降雨过程曲线，如图1所示。
(2)收集排水地区相关信息并建立雨水排水模型。收集排水地区的土地利用信息、管网参数信息、排放口的边界条件及历史降雨数据与检查井液位数据。在排水建模软件中，根据管网参数信息确定管网布局，依据地表高程、建筑物的分布划分节点的汇水面积及输入面积上相关参数，如面积大小F，用历史降雨数据及检查井液位数据校正模型。将步骤(1)中的降雨过程曲线导入校正后的排水模型作为降雨模块，并根据土地利用信息确定该地区平均降雨径流系数
Figure BDA0001614635900000041
(3)根据式1确定充满度解集E，如E＝[0.70,0.75,0.80,0.85,0.90,0.95,1.00]T；
(4)根据步骤(1)所确定的暴雨强度公式、步骤(2)中节点汇水面积大小F，及平均降雨径流系数
Figure BDA0001614635900000042
来计算管道设计流量：
Figure BDA0001614635900000043
选定E中不同管道设计充满度ei，由管道非满流水力计算公式：θ＝2cos-1(1-2ei)、
Figure BDA0001614635900000044
(n＝0.013；I为管道坡度，由管网布局得到)，以计算在选定充满度下雨水管网的管径，结果可用式2表示。
(5)调整ψ(ei)的计算管径，由于，ψ(ei)中可存在连续的管径值，而实际工程中的管径是非连续的，因此采用式2将ψ(ei)调整为其最接近的工程管径，例如连续管径值可能为248，而实际中只能取300。
(6)将充满度解集E的所有值完成计算，即可获得在不同设计充满度下的雨水管网初始优化设计方案集，结果表达如式3。
(7)建立雨水管网多目标优化算法模型。本发明方法以雨水管道管径d为决策变量，以管道投资成本最低与管道平均充满度最低为优化目标函数，并选定Borg为多目标优化算法进行求解。
(8)将步骤(6)中初始优化管径解集作为步骤(7)中Borg算法的初始种群，启动优化算法，并可快速获得雨水管网的帕累托优化解集。
将本发明的方法应用到某城市部分雨水管网中，该雨水管网系统服务面积约为3.2平方公里，服务人口约3万人左右。由于设计不合理以及气候变化的影响，使得该区域在雨季经常出现管道节点溢流和部分地区内涝情况，现需对该地区的雨水管网进行重新设计。如图3所示，该区域的水力模型包含102个节点，102个子汇水面积，102根管道和1个排水口。图4给出了本发明方法迭代5000次(约9.7h)和传统优化方法迭代150000次(约291.7h)结果，由此图可见，在找到较为接近的帕累托优化解集时，本发明方法比传统优化方法要快约30倍。因此，本发明方法对雨水管网进行优化设计时，优化效率得到了极大的提高，该优化效率优势在大型复杂雨水管网优化设计中更加明显。

Claims (1)

1.一种针对复杂雨水管网的快速多目标工程优化设计方法，其特征在于：包括以下步骤：
(1)根据雨水管网所在地区选定其对应的暴雨强度公式及其参数，并根据给定的暴雨重现期P，从降雨历史记录数据选取与之相对应的典型降雨过程曲线；
(2)建立雨水管网水力模型，并将步骤(1)中的降雨过程曲线导入该模型中的降雨模块；收集排水地区的土地利用信息、管网参数信息，在排水建模软件中，根据管网参数信息确定管网布局；
(3)确定充满度解集E，可表达如下：
E＝[e1,e2,...,en]T 式1
式中，ei为管道设计充满度，取值范围为0.6～1，i＝1,2,…,n，表示n个不同的设计充满度；
(4)根据步骤(1)的暴雨强度公式及其参数及步骤(2)确定的管网布局，计算管道设计流量Q，进而通过管道非满流水力计算公式，为每一个设定的管道设计充满度ei计算雨水管网中每一根管道的管径dk(ei),k＝1,2,…,K；K为雨水管道总数,结果可用公式表达如下：
ψ(ei)＝[d1(ei),d2(ei),...,dk(ei)]T 式2
式中ψ(ei)为设计充满度ei时优化管径设计方案；
(5)采用adjust函数对ψ(ei)中dk(ei)调整为其最接近的工程管径，表达如下：
S(ei)＝adjust{ψ(ei)} 式3
(6)确定雨水管网初始优化设计方案集；循环步骤(4)和(5)直到充满度解集E的所有值完成计算，即可获得在不同管道设计充满度下的雨水管网优化设计方案集，可表达为：
S＝{S(e1),S(e2),...,S(en)} 式4
(7)建立雨水管网多目标优化模型；
(8)将步骤(6)中雨水管网初始优化设计方案集作为步骤(7)中雨水管网多目标优化模型的初始种群，启动雨水管网多目标优化模型，以快速获得雨水管网帕累托优化解集。
CN201810281014.3A 2018-04-02 2018-04-02 一种针对复杂雨水管网的快速多目标工程优化设计方法 Active CN108824593B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810281014.3A CN108824593B (zh) 2018-04-02 2018-04-02 一种针对复杂雨水管网的快速多目标工程优化设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810281014.3A CN108824593B (zh) 2018-04-02 2018-04-02 一种针对复杂雨水管网的快速多目标工程优化设计方法

Publications (2)

Publication Number Publication Date
CN108824593A CN108824593A (zh) 2018-11-16
CN108824593B true CN108824593B (zh) 2020-07-07

Family

ID=64155244

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810281014.3A Active CN108824593B (zh) 2018-04-02 2018-04-02 一种针对复杂雨水管网的快速多目标工程优化设计方法

Country Status (1)

Country Link
CN (1) CN108824593B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110033139A (zh) * 2019-04-22 2019-07-19 扬州市职业大学(扬州市广播电视大学) 一种基于海绵城市理念的城市雨水优化系统和方法
CN111898222B (zh) * 2020-06-01 2024-02-09 浙江大学 排水系统设计方法、存储介质及系统
RU2746766C1 (ru) * 2020-08-29 2021-04-20 Общество с ограниченной ответственностью "Ассоциация инженеров и ученых по водоснабжению и водоотведению" Саморегулируемая система водоотведения поверхностного стока

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09259186A (ja) * 1996-03-27 1997-10-03 Hitachi Ltd 管網状態把握システム
CN103902783A (zh) * 2014-04-11 2014-07-02 北京工业大学 一种基于广义逆向学习差分算法的排水管网优化方法
CN107067119A (zh) * 2017-05-18 2017-08-18 上海宏波工程咨询管理有限公司 一种基于多目标分层选择的多工况给水管网优化方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150187025A1 (en) * 2014-01-01 2015-07-02 International Business Machines Corporation Optimization platform for water distribution networks

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09259186A (ja) * 1996-03-27 1997-10-03 Hitachi Ltd 管網状態把握システム
CN103902783A (zh) * 2014-04-11 2014-07-02 北京工业大学 一种基于广义逆向学习差分算法的排水管网优化方法
CN107067119A (zh) * 2017-05-18 2017-08-18 上海宏波工程咨询管理有限公司 一种基于多目标分层选择的多工况给水管网优化方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
基于关联矩阵的污水管网优化研究;史义雄;《中国农村水利水电》;20090915(第09期);全文 *
给水管网耦合优化调度模型的建立及应用;郑飞飞等;《中国给水排水》;20080901(第17期);全文 *
雨水管网优化设计探讨;高嘉谦;《山西建筑》;20131010(第29期);全文 *

Also Published As

Publication number Publication date
CN108824593A (zh) 2018-11-16

Similar Documents

Publication Publication Date Title
CN108824593B (zh) 一种针对复杂雨水管网的快速多目标工程优化设计方法
Hou et al. Spatial simulation of the ecological processes of stormwater for sponge cities
CN112001010B (zh) 一种控制分流制系统径流污染的雨水调蓄设施的设计方法
CN109033517B (zh) 一种提高雨水管网设计方案适用性的工程优化方法
Gao et al. A novel multi-objective optimization framework for urban green-gray infrastructure implementation under impacts of climate change
CN111476676B (zh) 一种城镇小流域分类分级智能调度方法
Timbadiya et al. HEC-RAS based hydrodynamic model in prediction of stages of lower Tapi River
CN113502893B (zh) 一种基于调蓄设施的下游雨水管道过流量的调控方法及系统
CN103967105A (zh) 居住组团雨水生态排放系统
CN102409747B (zh) 带浅层雨水调蓄结构的雨水排水系统
Chao-Hsien et al. Designing a rainwater harvesting system for urban green roof irrigation
CN116796480A (zh) 基于nsga-ⅲ算法的海绵城市防涝体系缓解城市内涝方法
Chen et al. Dimensionless parameter method for evaluating decentralized water reuse systems in buildings
CN103366309B (zh) 基于低影响开发的城市雨水利用规模确定方法
Arvand et al. Evaluating the impact of LID-BMPs on urban runoff reduction in an urban sub-catchment
CN116681564A (zh) 一种指导海绵城市建设的次暴雨总量控制率评估方法
CN116050705A (zh) 一种城市排水系统溢流与内涝控制最优方案分析方法
CN111898222B (zh) 排水系统设计方法、存储介质及系统
CN113722945A (zh) 一种用于辅助cso调蓄池规模设计的动态模拟和成本效益分析方法
Gao et al. Multi-objective optimization of sponge facility layout in built-up urban areas
CN202202394U (zh) 带浅层雨水调蓄结构的雨水排水系统
CN105549396A (zh) 一种基于在线动态规划的城市明渠排水系统控制方法
BENTALHA Evaluation of the hydraulic and hydrology performance of the green roof by using SWMM
CN114444184B (zh) 一种基于高精度水动力模型的lid设施优化设计方法
Liu et al. Determination of the minimum soil infiltration rate of sunken green space considering the annual runoff collection ratio, sunken depth and sunken green space area of Hefei city, China

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