基于动态规划的雨水泵站最优服务半径研究-李翠梅.pdf

《基于动态规划的雨水泵站最优服务半径研究-李翠梅.pdf》由会员分享，可在线阅读，更多相关《基于动态规划的雨水泵站最优服务半径研究-李翠梅.pdf（6页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第33卷第4期2016年7月深圳大学学报理工版JOURNAL OF SHENZHEN UNIVERSITY SCIENCE AND ENGINEERINGV0133 No4July 2016【土木建筑工程Architecture and Civil Engineering】基于动态规划的雨水泵站最优服务半径研究李翠梅，龚希博苏州科技大学环境科学与工程学院，江苏苏州215009摘要：针对雨水提升泵站最优服务半径的研究缺失，通过建立雨水管网系统满流状态下的数学模型。以雨水管网系统运行参数为约束条件，以雨水管段埋深构建系统状态转移方程，研究了应用动态规划方法求解雨水提升泵站最优服务半径过程，并对雨水

2、系统初始管段展开了实例研究结果表明，应用动态规划方法求解雨水提升泵站最优服务半径可以实现雨水系统的全局优化以中国不同地区最大允许埋深为约束条件，应用所建立数学模型计算了相应的最优服务半径，可为雨水系统优化设计与规划提供参考关键词：市政工程；雨水系统；提升泵站；最优服务半径；动态规划；埋深中图分类号：TU 992 文献标志码：A doi：103724SPJ1249201604388The optimal service radius of stormwater pumpstation based on dynamic programming methodLi Cuimei+and Gong Xi

3、boSchool of Environmental Science and Engineering，Suzhou University of Science and Technology，Suzhou 215009，Jiangsu Province，PRChinaAbstract：By taking the stormwater pipe operation parameters as constraint conditions and using the stormwater pipesection depth to build system state transition equatio

4、n，we establish the mathematical model for stormwater pipe system and study the optimal service radius of stormwater pump stationWe use the dynamic programming method tosolve the optimal model and carry out a case studyThe results show that the application of dynamic programmingmethod for solving the

5、 optimal service radius could achieve global optimization in stormwater systemCombined withthe maximum allowable depths in different regions of China，we calculate the corresponding optimal serviee radius ofstorm water pump station which can provide a referenee for improving the planning of stormwate

6、r systemKey words：municipal engineering；stormwater system；pump station；optimal service radius；dynamic programmingmethod；buried depth随着我国城镇化进程的快速发展，雨水排放系统暴露出越来越多的问题，每年58月，大规模的城市内涝似乎正在成为常态2012年7月的北京，2014年5月的深圳，不论是古老的历史文化名城，还是年轻的现代大都市，似乎都无法经受雨水的考验对比国内外雨水系统设计规范与设计过程发现，在雨水系统规划设计阶段，主观人为的参数设置往往导致雨水系统建成后在关键

7、时刻不能发挥Received：20151212：Accepted：20160420Foundation：National Natural Science Foundation of China(51109153)：Science and Teehnolo静r Foundation of Housing and Construction Department of Jiangsu Province(20142I)96)十Corresponding author：Professor Li CuimeiEmail：cuimeili163cornCitation：Li CuimeiGong XiboT

8、he optimal service radius of stormwater pump station based On dynamic programming methodJJournal of Shenzhen University Science and Engineering，2016，33(4)：388-393(in Chinese)万方数据第4期 李翠梅，等：基于动态规划的雨水泵站最优服务半径研究 389排涝作用，其中雨水提升泵站最优服务半径的研究基本处于空白，不论是规范上，还是实际规划设计过程中，都是人为主观确定，缺少科学数据发达国家在经历排水系统的创建、发展和治理之后，自20

9、世纪70年代起就开始致力于暴雨防涝系统的管理，重点集中在限制雨水和合流制系统溢流次数和流量，以及研究溢流污水的处理途径由于城市排水管理和排水过程的复杂性，要想对雨水径流过程有清晰的认识，必须借助现代信息方法和技术手段进行准确模拟和预测因此，世界上涌现了上百种城市暴雨排放模型，用于雨水管网的规划、设计以及运行管理，以达到减少溢流量、减少管道超载和地面淹水以及减少溢流的污染负荷等多个目标，降低暴雨造成的城市危害1。71城市降雨与径流过程模型包括径流过程模拟(水文和水力学模拟)和径流水质模拟两部分陈礼国等。2。自20世纪80年代起对雨水管网汇流模型进行了系列研究，先后对瞬时单位线法求雨水管网系统人流

10、流量过程线的数值计算方法、城市排水系统非恒定流模拟模型、雨水道变径流系数、瞬时单位线法推求雨水管网人流流量过程线、城市雨水管网系统地面径流损失规律、牛顿迭代逼近法求解雨水管道非线性运动波方程、城市雨水管网非线性运动波法模拟技术、特殊情况下雨水管道设计流量的确定方法、瞬时单位线法求雨水管网系统人流流量过程线的数值计算方法、遗传算法推求瞬时单位线参数并计算流量过程线、城市排水系统非恒定流模拟模型、雨水管渠设计参数折减系数的理论、节点递归算法优化雨水管网和城市暴雨强度公式推求系统的开发等进行了深入的研究近20年来，随着城市径流问题的日益突出，各种精度较高的城市水文和水力计算模型已被广泛用于雨水管道系

11、统的规划、设计和管理目前最具影响的有：英国环境部及全国水资源委员会的沃林福特模型(Wallingford model)、美国陆军工程师兵团水文学中心的“暴雨”模型(storage，treatment，overkqow，runoff mode storm)和美国环保局的雨水管理模型(storm water management mode，SWMM)等这些模型可对整个城市降雨、径流过程进行较为准确的模拟虽然国内外学者对城市雨水径流过程模拟已展开了大量研究并取得了显著成果，但关于雨水泵站最优服务半径的研究尚未见报道本研究采用动态规划的数学方法，以全局最优为目标，建立并求解了不同地质条件下雨水提升泵站

12、的最优服务半径，为雨水系统规划提供了一个新的参考1 雨水提升泵站最优服务半径数学模型雨水提升泵站最优服务半径是指在既定最大允许埋深情况下，雨水提升泵站能够服务的最大半径范围，在这个范围内，降雨经地表径流、汇流后进入地下雨水管网系统有组织收集和输配，随着管道长度的增加，坡降越来越大，埋深越来越大，当埋深达到最大允许值时，必须设置提升泵站，在泵站内将机械能转变为雨水的位能和动能，提升到埋深较浅的位置，继续向下游管网输送或者就近排放到水体中在这个过程中，可以清楚地看到泵站的排涝作用，上游汇集的雨水在泵站开启的情况下，被不断地排放到水体或者下游管网，使上游汇水区的地表径流和汇流过程得以持续如果提升泵站

13、规划的位置比较好，其服务的范围比较大，整个城市所需要设置的提升泵站的数量就会减少，雨水提升系统的效率就会提高，这是雨水泵站最优服务半径的理念根据我国地质结构，地下水水位较高、气候较为温暖的地区，例如江南地区，考虑到施工因素，地下管线最大允许埋深一般为45 m；在地下水水位较低、气候偏冷的地区，地下管线最大允许埋深一般为810 m，基本不会超过10 m所需提升水泵的最大扬程也就是十几米，对水泵来说可以满足那么为什么还会经常发生雨水经提升泵站排不出去的现象?要想回答这个问题，还需要从雨水提升泵站规划设计的位置、最优服务半径值的选取等方面出发，进行数学建模，展开系统研究11 数学模型及约束条件表达雨

14、水泵站最优服务半径数学模型为F=max+1) (1)=J其中，F为提升泵站的服务半径；凡为管道管段数；乙，为从矗阶段到k+1阶段对应的设计管段长度结合我国室外排水设计规范(GB 50014一万方数据390 深圳大学学报理工版 第33卷2006)m 3相关规定，模型必须满足以下约束条件：stAH=(iizi)H。一H。i。i=l，。i。Ii，。王，mi。i11。axH。H。H。 (2)矿。矽l+l(告)(告)i(告)一D。D，Di亿其中，。i。为最小允许设计坡度；，。为最大允许设计坡度；，。为管段i(江1，2，n)的设计坡度；矽。油为最小允许设计流速；。，为最大允许设计流速；i为管段i的设计流速

15、；Di为管段i的管径；亿为标准规格管径集；p和D分别为与管段i相邻的上游管段的流速和管径；H为最小允许埋深；H。为最大允许埋深(视当地情况而定)；Hi为节点i埋深；(hD)。j。为最小允许设计充满度；(hD)。为最大允许设计充满度；(hD)i为管段i的设计充满度12初始条件雨水管道支管一般从小区内部或者街区开始铺设，起端管段因汇流面积小，径流量较小为了降低起端管段的坡降，大部分定义为不计算管段，管径直接采用最小管径200 mm，坡度采用最小设计坡度，之后随着汇水面积增大而逐级增大管径，与之对应管段流量也相应增大图1为雨水管道系统剖面示意图，其中，H，为起端埋深；也为最大允许埋深；1-7为雨水节

16、点编号管道系统由多个管段组成，每一管段有不同的长度，每一管段最大允许坡降与最小允许坡降可以结合规范先推算出来，如表1所示图1 雨水管道系统剖面示意图Fig1 Sketch of the profde ofstormwater pipeline system表1不同管径设计坡度范围Table 1 Grade ranges for different diametersDmm ， ，。 Omm mi|。200 0400 2463 500 0113 0 0625250 0350 1829 600 0088 5 0490300 0300 1434 700 0072 l 0399350 0182 104

17、4 800 0060 3 0334400 0152 0873 900 0051 6 0285450 0130 0746 l 000 0050 0 02422动态规划方法21动态规划方法原理动态规划(dynamic programming，DP)是运筹学的一个分支，是求解决策过程最优化的数学方法20世纪50年代初Bellman等在研究多阶段决策过程(multistep decision process，MDP)的优化问题时，提出了著名的最优化原理(principle of optimality)，把多阶段过程转化为一系列单阶段问题，利用各阶段之间的关系，逐个求解，创立了解决这类过程优化问题的新方

18、法动态规划动态规划的特点是在局部不是最优的情况下可以获得全局最优，善于分析解决多阶段最优化问题所谓多阶段最优化问题，是指这样一类活动的过程，即可以将全局活动过程划分为若干个相互联系的阶段，在它的每一个阶段都需要做出决策，并且一个阶段的决策确定以后，会影响下一阶段的决策，从而影响整个活动的过程各个阶段所确定的决策就构成一个决策系列，通常称为一个策略由于每一个阶段可供选择的决策往往不止一个，因而就形成许多策略可供选取一般来说，对应每一个活动过程，策略会有其相应的活动效果，该效果经常可以用一个数量指标来衡量不同的策略，效果一般也不相同，多阶段决策问题就是要在允许选择的那些策略中间，选择一个使在预定的

19、标准下达到最佳效果的最优决策序列，即最优策略动态规划基本方程的逆序形式为(xI)= opt n(戈，“一1)+“一IDk(xk)丘一(x一。)，k=凡，凡一l，2，1 (3)其中，边界条件有L+(戈川)=0瓤(戈。)=K(z。，万方数据第4期 李翠梅。等：基于动态规划的雨水泵站最优服务半径研究 391u。)第k阶段的状态为X。，其决策变量u。表示状态处于戈的决策，状态转移方程为工=瓦(x。，2。)，k阶段的允许决策集合记为D。(z。)，K(戈。，2。)为指标函数求解过程中，由边界条件从k=n开始，由后向前逆推，逐阶段求出最优决策和过程的最优值，直到最后求出工(x，)，即得到问题最优解22雨水管

20、网动态规划模型雨水管网由多个管段组成，每一管段可以选择不同的管径、坡度、流速和充满度等设计参数，因而各个管段有多个策略可以选择，不同的策略选择，最终会影响从起始管段到雨水提升泵站之间的埋深和总长度，即雨水提升泵站的服务范围如何在既定的最大允许埋深下，在各个管段选择合理(局部不一定最优)的设计参数，最终实现雨水提升泵站服务半径最大，是应用动态规划来求解的核心问题221 划分阶段、确定状态变量与决策变量雨水管网系统基本为枝状管网，从起始段到雨水提升泵站之间由多个管段组成，每一个设计管段可以划为一个阶段，多阶段划分如图2选取管径作为雨水管网的状态变量，选取管段的水力坡度为决策变量，决策变量与管段长度

21、的积即为坡降，当各个阶段坡降的代数和接近最大允许埋深时，需要图2 雨水管道系统多阶段划分示意图Fig2 Schematic diagram of rainmulti-phase pipeline systemY6(S6)=maxL67+f7(S7)A(S5)=maxL56+八S6)五(S。)=maxL4，+六(S，)厶(S3)=maxL34+f4(S4)五(S：)=maxL：，+六(S，)Z(SI)=maxL2+正(S2)设置提升泵站，此时就确定了从起始管段到提升泵站之间的设计管段数量并可以得出相应的各个管段长度代数和，即为泵站的最大服务范围222 系统状态转移方程系统中，如果第k个阶段状态变

22、量为戈。，作出的决策为“。，那么第k+l阶段的状态变量也可被确定，那么可以用状态转移方程表示这种演变规律选择上下游管段间埋深建立管段间的状态转移方程，表示上下游管段之间相互影响：S“l=S一L(+1)M (4)其中，S川为k阶段的下游埋深，即k+1阶段的上游埋深；S。为k阶段的上游埋深；2。为第k设计管段的坡度223指标函数指标函数是系统执行某一策略所产生结果相应指标的数量表示，是衡量策略优劣的指标，本研究中指标函数指雨水泵站的最优服务半径：K(S，u)=L+1) (5)最优指标函数为五(戈。)= optK(x。，u。)+以+。(x川)(6)“IEOk(。)根据以上最优指标函数和状态转移方程，

23、结合最大埋深要求，可求得中途提升泵站的最优服务半径、相应的管段长度、埋深及管径等参数组合23模型求解以图2为例，假设雨水系统最大允许埋深S，为5 m，将管段最小设计坡度、起始不计算管段最小坡度和最小管径代人式(4)、式(5)和式(6)，可以得出系统指标函数组为=3 846769S6=maxL56+3 846769S6=3 846769S40169L5604L45(7)=3 846769S30169L5604L45一i307L34=3 846769S20169L5604L451307LM一169L23=3846769SI一0169L5604L451307L34169L232076L12万方数据深

24、圳大学学报理工版 第33卷将雨水管道最小允许埋深|s。代入式(7)，可得出各个管段的指标函数值，最终得到雨水提升泵站最优服务半径以苏州地区为例，按照上述方法可得出雨水泵站最优服务半径为2 km左右24数据校验选取苏州地区多个雨水管网设计实例，结合室外排水设计规范(GB 5014-2006)3j相关规定，在主要设计参数满足规范的条件下，测算埋深为5 m时，雨水泵站的最大服务范围为193252km，与本研究结果较为接近3 讨 论我国室外排水设计规范 (GB 50142006)113i中对各个地区雨水管网最大允许埋深有所规定和建议采用值，但是雨水提升泵站规划设置的位置及其所能保障的最优服务半径一直未

25、给出建议值或者规定这一情况易造成雨水系统规划时，雨水提升泵站规划位置主观性强，科学性低，且易造成局部区域雨水得不到提升，形成排涝死角，故而发生有地下雨水管网系统的同时，还造成大面积积水或内涝的情况41 5I我国给水排水专业相关的多本规范中，对涉及到人身财产安全的系统一般都设置了强制条例，例如，防火规范中消防系统有许多强制条例，设计时必须遵守，目的是提高系统的安全性人们以前认为雨水系统不会直接造成人身财产安全问题，但近年来，因雨水系统排放不畅而造成的城市内涝事件中，人员伤亡常有发生，且造成的整体社会损失有明显上升趋势因此，建议在制定雨水系统相关规范时参考消防系统的做法，设置强制性规范鉴于我国雨水

26、管道最大允许埋深的情况，应用本动态规划方法，可以得出我国雨水提升泵站最优服务半径，如表2所示表2 雨水提升泵站最优服务半径Table 2 The optimal service radius ofstormwater pumping station最大允许埋深m 雨水提升泵站最优服务半径km20253540结 语针对我国近年来城市内涝频发及室外排水设计规范(GB 5014-2006)。1 3。中雨水提升泵站设置存在的问题，提出了雨水提升泵站最优服务半径的概念，并建立了数学模型，应用动态规划原理与方法对模型进行了求解研究表明，应用动态规划方法求解雨水提升泵站最优服务半径，可简化求解过程，缩短计算

27、次数和时间，提高计算效率以6管段组成的雨水管网为例，常规工程计算次数为66次，而动态规划方法仅为36次应用研究方法，推算了我国雨水管道最大允许埋深从58 m的雨水提升泵站最优服务半径，供设计和管理决策参考鉴于城市雨水系统的重要性，建议进一步展开城市雨水提升泵站最优服务半径的研究，并组织专家论证，逐步形成具有强制性的规范条文，以提高规划设计的科学性和合理性基金项目：国家自然科学基金资助项目(5l109153)；江苏省住建厅科技资助项目(2014ZD96)作者简介：李翠梅(1974一)，女，苏州科技大学教授研究方向：城市水资源与给水排水工程设计运行的最优化Email：cuimeili163com引

28、 文：李翠梅，龚希博基于动态规划的雨水泵站最优服务半径研究J深圳大学学报理工版，2016，33(4)：388393参考文献References：1Ki S J，Kang J，Lee S W，et a1Advancing assessmentand design of stormwater monitoring programs using a self-organizing map：characterization of trace metal concentration profiles in stormwater runoffJWater Research，2011，45(14)：4183

29、41972Hayes D FDynamic programming methodology for prioritizing sewerage projectsJJournal of Water ResourcesPlanning&Management，2014，137(2)：193-2043Barbosa A E，Fernandes J N，David L MKey issues forsustainable urban stormwater managementJWaterResearch，2012，46(20)：678767984Zgheib S，Moilleron R，Saad M，e

30、t a1Partition of pollu-tion between dissolved and particulate phases：what aboutemerging substances in urban stormwater catchmentsJWater Research，2011，45(2)：913-9255Xiao F，Simcik M F，Gulliver J SPeffluoroalkyl acids in万方数据第4期 李翠梅，等：基于动态规划的雨水泵站最优服务半径研究 393urban stormwater runoff：influence of land use

31、J IWater Research，2012，46(20)：660166086Vezzaro L，Ledin A，Mikkelsen P SIntegrated modellingof priority pollutants in stormwater systemsJPhysicsand Chemistry of the Earth(Parts ABC)，2012，4243州(2)：42-517Trowsdale S A，Simcock RUrban stormwater treatmentusing bioretentionJJournal of Hydrology，201 1，397(3

32、4)：1671748陈礼国，刘遂庆污水管网优化计算中动态规划变异形式和比较J同济大学学报自然科学报，1999，27(3)：355360Chen Liguo，Liu SuiqingThe variation and comparison ofdynamic programming in optimization calculation of sewagepipe networkJJournal of Ton舀i University NaturalScience Edition，1999，27(3)：355360(in Chinese)9刘遂庆排水管网优化坡度设计计算J给水排水，2008，34(

33、8)：115116Liu SuiqingSlope design and calculation about optimization of drainage pipe networkJWater Supply andDrainage，2008，34(8)：115一116(in Chinese)10周玉文，翁窈瑶，汪明明流量折减系数对雨水系统设计标准的影响J同济大学学报自然科学版，201 1，39(10)：15061527Zhou Yuwen，Weng Yaoyao，Wang MingmingFlow foldreduction coefficient have effect of rainw

34、ater system designcriteriaJJournal of Tongji University Natural ScienceEdition，2011，39(10)：1506-1527(in Chinese)11周玉文，翁窈瑶，汪明明雨水管渠设计参数折减系数m的理论研究J北京工业大学学报自然科学报，2011，37(9)：13301337Zhou Yuwen，Weng Yaoyao，Wang MingmingThe theo-retical research of parameters reduction coef!ficient m ofdesign rainwater pip

35、e drainageJJournal of BeijingUniversity of Technology Natural Science Edition，201 1，37(9)：13301337(in Chinese)12车伍，刘燕，欧岚，等城市雨水径流面污染负荷的计算模型J中国给水排水，2004，20(7)：5658Che Wu，Liu Yan，Ou Lan，et a1Calculation model forarea source pollution load of urban rain water runoffJChina Water and Drainage，2004，20(7)：5

36、6-58(inChinese)13GB 50014-2006室外排水设计规范(2014修订版)SGB 50014-2006 Code for design of outdoor drainage(2014 revision)S(in Chinese)14徐祖信，汪玲玲，尹海龙基于水质特征因子和MonteCarlo理论的雨水管网混接诊断方法J同济大学学报自然科学版，2015，43(11)：17151721Xu Zuxin，Wang Lingling，Yin HailongA hybrid methodof rainwater pipe network based on water quality

37、 charaete-ristie factor and Monte Carlo theoryJJournal of TongjiUniversity Natural Science Edition，2015，43(11)：1715-1721(in Chinese)15李爽，邢国平，刘洪海，等基于低影响开发的天津市绿色雨水基础设施应用研究J水土保持通报，2013，33(2)：147150Li Shuang，Xing Guoping，Liu Honghai，et a1Study onthe application of green rainwater infrastructure in Tianjincity based on low impact developmentJBulletin ofSoil and Water Conservation，2013，33(2)：147150(inChinese)【中文责编：坪梓；英文责编：之聿】万方数据