曹相生1，姜安玺1，依学农2
(1.哈尔滨工业大学市政环境工程学院；2.山东建筑工程学院环境工程系) 摘要：根据雨水管道的计算特点，在分析目标函数的各种影响因素之后，提出一种新的雨水管道优化设计方法，得出了较为满意的结果。
关键词：雨水管道；优化设计；枚举管径法
中图分类号：TU992.2
文献标识码：C
文章编号：1000-4602(2000)10-0040-02 雨水管道的优化设计包括三方面：①雨水径流模型的建立；②管线平面布置的优化选择；③一定管线下的管道优化设计。对第三方面，前人大多数采用线性规划法、非线性规划法、动态规划法、遗传算法等间接优化的方法并编制了相应的程序［1～2］。这些程序针对污水管道并人为忽略和简化了某些约束条件和控制变量，因此对于计算方法不同于污水管道的雨水管道不适用或得不到最优解。 1 雨水管道的优化设计模型 依据现行的《室外排水设计规范》，雨水管道采用极限强度法和均匀流公式来计算。另外，雨水管道计算应满足管径、流速、埋深等方面的约束条件。具体参见文献［2～6］。
另外，考虑到管道之间的连接，本程序还考虑了一个约束条件，即管道管内底高程应小于(接纳支管时)或大于(接至干管时)某一设定高程值。
雨水管道优化设计的目的是使整个管系在投资偿还期内运行成本和按标准偿还期折算的年度投资最低。目标函数的意义是给各方案提供一个统一标准做比较，本程序把文献［6］采用的费用函数作为目标函数。费用函数形式如下：
 (1) 式中 M ——雨水管道总费用
Di——第i管段设计管径
Ii——第i管段设计坡度
Li——第i管段管道长度
I0i——第i管段地面坡度
H0i——第i管段起点埋深
n——设计管段数
对雨水管道的优化即是求式(1)的最小值Mmin。 2 程序优化思想 当雨水管道定线以后，各管段的长度、地面坡度、第一管段的设计流量均已确定，因此 式(1)可写成如下形式：  (2)
式(2)中的自变量并不是相互独立的，他们之间的制约关系可表示如下：
H0j＝g1(D1，D2，…，Dj，I1，I2，…，Ij-1，ΔH1，ΔH2，…，ΔHj，H01) (3)
Ii＝g2(Di,Qi) (4)
Qj＝g3(D1,D2，…，…Dj-1) (5)
式中 ΔHi——第i管段跌水值
i=1,…，n j=2，…，n
综合式(2)、(3)、(4)、(5)，雨水管道总费用函数可表示如下：
M＝p(D1,D2,…,Di，…，Dn，ΔH1，ΔH2，…，ΔHi，…ΔHn，H01) (6)
从式(6)可以看出，影响雨水管道总费用的因素只有第一设计管段起点埋深，各管段跌水值和各管段设计管径。分析费用函数可知，H01和ΔHi(i=1，…，n)取最小值时，M才可能取最小值，因此式(6)又可写成如下的简单形式：
M=q(D1，D2，…，Di，…Dn) (7)
求解式(1)的最小值变成了在满足约束条件下寻求各设计管段管径的最优组合。
对式(7)直接求解其数学表达式是很困难也是没有必要的。事实上式(7)中的自变量也不是独立变量，它们之间存在着依次制约的关系。对式(7)最小值的求解最简单的方法就是遍历各设计管段可能的管径组合，进行直接寻优，这可利用计算机来完成。依据以上分析，确定了枚举管径法优化雨水管道的程序计算思路如下：
首先输入原始数据，依据管径和流速的约束条件计算第一管段可利用管径，选择一种管径后，计算管段其他设计参数，验证是否符合其他约束条件。当埋深或管内底标高不符合约束条件时，采用两种方案，一是加大管道坡度，二是设跌水。找到符合约束条件下第一管段各参数后，依据管径和流速的约束条件确定下一管段可利用管径，根据相同办法计算其管段设计参数直至计算完全部设计管段，找到第一种符合约束条件的管径组合并计算此种组合的总费用。重复上述过程找到其他的管径组合，比较各种组合的总费用，把总费用最低的组合作为最优的设计方案保留起来。 3 程序实际运行结果 为了验证程序的实用性，采用文献［7］中的例题做了比较，见表1。表中括号内的数据为程序运算结果。为节约篇幅，表中只列出了部分计算数据，同时为了增加可比性，对文献［7］采用了与本程序相同的费用函数并计算了管道造价。 表1雨水干管水力计算表管段
编号 | 设计流量(L/s) | 管径(mm) | 坡度(%) | 流速
(m/s) | 管道输水能力(L/s) | 起点埋深(m) | 终点埋深(m) | 造价(万元) | 1-2
2-3
3-5
5-9
9-10
10-11
11-12
12-16
16-17
17-18
18-19 | 94.58(94.6)
163.98(163.6)
233.78(233.1)
335.00(333.4)
503.42(499.0)
518.38(499.0)
557.21(550.6)
681.32(673.2)
681.32(673.2)
791.40(798.9)
853.41(867.7) | 400(400)
500(500)
600(600)
700(700)
800(800)
800(800)
800(800)
900(800)
900(800)
900(800)
900(800) | 0.21(0.206)
0.19(0.188)
0.15(0.144)
0.14(0.130)
0.15(0.142)
0.16(0.142)
0.18(0.173)
0.15(0.259)
0.15(0.259)
0.20(0.365)
0.23(0.431) | 0.76(0.75)
0.84(0.83)
0.84(0.82)
0.90(0.87)
1.02(0.99)
1.05(0.99)
1.12(1.10)
1.10(1.34)
1.11(1.34)
1.29(1.59)
1.37(1.73) | 96.00(94.6)
165.00(163.6)
240.00(233.1)
350.00(333.4)
520.00(499.0)
530.00(499.0)
560.00(550.6)
700.00(673.2)
700.00(673.2)
810.00(798.9)
870.00(867.7) | 1.30(1.30)
1.75(1.74)
2.04(2.03)
2.37(2.25)
2.12(2.09)
2.27(2.24)
2.43(2.38)
2.75(2.59)
2.94(2.95)
3.11(3.26)
3.41(3.80) | 1.65(1.64)
1.94(1.93)
2.27(2.15)
2.02(1.99)
2.27(2.24)
2.43(2.38)
2.65(2.59)
2.94(2.95)
3.11(3.26)
3.41(3.80)
3.74(4.43) | 0.4521(0.4512)
0.4639(0.4612)
0.6171(0.6133)
1.1019(1.0602)
0.8577(0.8496)
0.8984(0.8869)
1.1373(1.1195)
1.7447(1.5129)
1.4526(1.3158)
1.9125(1.8117)
2.0414(2.0426) | | 注总造价：文献［7］中为12.679 6万元；枚举法为12.124 9万元(共找到满足约束条件的管径组合2 563种)。 | 由表1可见，本程序运算结果比例题计算的费用减少0.5547万元，降低了4%。另外，本程序应用在山东省青州市排水规划中效果良好。实践证明，无论地形如何复杂，本程序都能满足要求。 4 结语 ① 本程序采用枚举管径法从各种符合约束条件的管径组合中寻优，因此结果应该是最优的。
② 是否跌水及跌水值完全由计算机决定，避免了人工确定时的盲目性。
③ 对管内底标高做了约束，方便了管道之间的连接。
④ 枚举管径法的思路可推广到整个管网的优化，这为雨水管网和污水管网的全局优化提供了可行的办法。 参考文献：
［1］张景国，李树平.遗传算法用于排水管道系统优化设计［J］.中国给水排水，1997，13 (3)：28-30.
［2］张绍园.排水管网计算程序中最优方案的探讨—与彭永臻同志商榷［J］.中国给水排水，1994，10(2)：21-24.
［3］彭永臻，崔福义.给水排水工程计算机程序设计［M］.北京：中国建筑工程出版社，1994.
［4］于尔捷，张杰.给水排水工程快速设计手册(2)排水工程［M］.北京：中国建筑工程出版社，1996.
［5］彭永臻，等.雨水排水系统计算机程序设计［J］.中国给水排水，1991，7(4)：38-40.
［6］依学农.计算机在给排水工程中的应用［M］.山东建工学院，1994.
［7］重庆建工学院.排水工程(上)(第三版)［M］.北京：中国建筑工程出版社，1996.
电话：(0451)6282207
收稿日期：2000-02-05 | 
