随着工业化和城市化进程的逐步加快,我国在输水管道的建造和管理上的投入越来越大。本文主要对单水源重力树状管网系统进行了优化设计,建立了线性规划数学模型,并应用MATLAB编制优化设计计算程序,给出了投资少、能耗低的输水管网的优化设计方案。最后,给出了数值实验,进一步验证了模型的合理性、实用性,具有非常大的经济和社会效益。 With the gradual acceleration of industrialization and urbanization, China’s investment in the construction and management of the pipeline is getting bigger and bigger. In this paper, on the single water gravity tree pipe network system optimization design, the linear programming ma-thematics model was established, and using the MATLAB procedures to optimize design calcula-tion, an optimal design scheme for water pipe network was given. Finally, a numerical experiment was given to further verify the rationality and practicability of the model, which has great economic and social benefits.
陈昊东,杜家康,孙洪春,任庆军
临沂大学数学与统计学院,山东 临沂
收稿日期:2018年5月11日;录用日期:2018年5月22日;发布日期:2018年5月30日
随着工业化和城市化进程的逐步加快,我国在输水管道的建造和管理上的投入越来越大。本文主要对单水源重力树状管网系统进行了优化设计,建立了线性规划数学模型,并应用MATLAB编制优化设计计算程序,给出了投资少、能耗低的输水管网的优化设计方案。最后,给出了数值实验,进一步验证了模型的合理性、实用性,具有非常大的经济和社会效益。
关键词 :输水管道,线性规划,单纯形法,优化设计
Copyright © 2018 by authors and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
在村镇供水以及工矿企业给水系统的建设初期,一般都采用树状管网,因为树状给水管网一般投资比较小,以后会随着用水规模的扩大,由树状管网逐渐发展为环状管网。树状管网形式较多,不同的形式有不同的特点,也会满足不同的情形。本文主要对水源位置较高时,采用的重力树状管网系统进行优化。单水源重力树状管网系统要解决的问题是:在保证各节点所需流量和服务水压的情况下,确定管网中各个管段尺寸的最优值,使得整个管网投资最少。目前,对管网系统的优化已经引起了专家和学者的重视。王文芬 [
首先,我们先对涉及的一些术语进行解释说明。在给水管网系统中,节点水压标高是指管道上节点压力水头与该点地面高程之和。给水管网系统中具有一定的水压,用来保证用户的用水,我们将这个水压称为最小服务水压。服务水压标高就是指管道上节点最小服务水压与该点地面高程之和 [
一般来说,工程上采用标准管径,也就是说每一管段可以由一种或多种标准管径组合成。由于标准管径档次太多,计算量太大,为方便计算,可以通过流速来限制管径。为防止水击破坏,管道内的流速应该限制不高于3 m/s;为防止淤积,管道内的流速也不应该低于0.6 m/s [
0.6 ≤ V i j ≤ 3 (1)
V i j = 4 Q i π d j 2 (2)
其中 Q i 表示管段设计流量(m3/s),dj表示管道计算内径(m)。
在第i管段中,可以根据(1)式的条件来选择标准管径。若这样的标准管径档次有 m ( i ) 个,则这 m ( i ) 个标准管径组成第i管段的备选管径组。
以单水源重力树状管网投资最小为目标函数
min W = ∑ i = 1 n ∑ j = 1 m ( i ) a i j l i j (3)
其中W表示重力树状管网投资(元),n表示管段数, m ( i ) 表示第i管段的备选管径数, a i j 表示第i管段中的第j种标准管径管道单价(元/m), l i j 表示第i管段中选用的第j种标准管径的管长(m)。
每一管段各标准管径管长之和等于该管段长,也就是
∑ j = 1 m ( i ) l i j = L i ( i = 1 , 2 , ⋅ ⋅ ⋅ , n ) (4)
其中 L i 表示管段长(m)。要求各个节点水压标高均不得低于该节点服务水压标高,即
∑ i = 1 H ( r ) ∑ j = 1 m ( i ) p i j l i j ≤ E 0 − E r ( r = 1 , 2 , ⋅ ⋅ ⋅ , R ) (5)
其中r表示节点序号,R表示节点总数, H ( r ) 表示水源至第r点的管段数, 表示水源水面高程(m), E r 表示节点服务水压标高(m),pij表示管道水力坡度,对于混凝土管及钢筋混凝土管,可用巴甫洛夫斯基公式 [
p i j = V i j 2 e j 2 R j (6)
其中 R j 表示水力半径(m), R j = d j / 4 ,
e j = 1 N R j 1 / 6 (7)
其中N表示粗糙系数,混凝土管及钢筋混凝土管一般取0.013~0.014。另外, l i j ≥ 0 。
由以上分析,可建立单水源重力树状管网系统的数学模型
min W = ∑ i = 1 n ∑ j = 1 m ( i ) a i j l i j s .t . ∑ j = 1 m ( i ) l i j = L i ( i = 1 , 2 , ⋅ ⋅ ⋅ , n ) ∑ i = 1 H ( r ) ∑ j = 1 m ( i ) p i j l i j ≤ E 0 − E r ( r = 1 , 2 , ⋅ ⋅ ⋅ , R ) l i j ≥ 0
这是一个带约束的线性规划模型,优化变量为lij,采用单纯形法求解,下面根据例子对数学模型进行进一步的解释。单水源重力给水树状管网布置如图1所示。采用钢筋混凝土管,管道单价见表1,管道的粗糙系数 N = 0.013 ,当管道起点与终点高差为5 m ( E 0 − E 3 = 5 ),确定单水源重力树状管网优化设计方案 [
Step 1:划分管段
Step 2:确定各管段备选管径组
根据式(2)确定各个管段的备选管径组见表2。
Step 3:计算水力坡度
根据巴氏公式(6),计算各管段水力坡度见表2。
图1. 单水源重力树状管网计算图
管径(m) | 单价(元/m) |
---|---|
0.1 | 20.3 |
0.15 | 29.1 |
0.2 | 36.0 |
0.25 | 46.7 |
0.3 | 56.7 |
0.4 | 83.6 |
0.5 | 103.7 |
0.6 | 144.3 |
表1. 管道单价
管段序号i | 管段名称 | 管段流量(m3/s) | 备选管径序号j | 备选管径(m) | 流速(m/s) | 水力坡度 |
---|---|---|---|---|---|---|
1 | 1~2 | 0.150 | 1 | 0.5 | 0.76 | 0.001578 |
2 | 0.4 | 1.19 | 0.005182 | |||
3 | 0.3 | 2.12 | 0.02401 | |||
2 | 2~3 | 0.025 | 1 | 0.25 | 0.51 | 0.00176 |
2 | 0.2 | 0.80 | 0.00579 | |||
3 | 0.15 | 1.42 | 0.02683 |
表2. 管段备选管径组及其水力坡度
Step 4:确定目标函数
根据式(3)和表1,确定目标函数如下:
min W = 103.7 l 11 + 83.6 l 12 + 56.7 l 13 + 46.7 l 21 + 36 l 22 + 29.1 l 23
Step 5:确定约束条件
根据式(4)和图1,确定管长约束条件如下:
l 11 + l 12 + l 13 = 500 l 21 + l 22 + l 23 = 650
根据式(5)和表2,确定压力约束条件如下:
0.001578 l 11 + 0.005182 l 12 + 0.02401 l 13 + 0.00176 l 21 + 0.00579 l 22 + 0.02683 l 23 ≤ 5
非负约束条件如下: l i j ≥ 0 ( i = 1 , 2 , 3 ; j = 1 , 2 , 3 ) 。
通过整理以上式子,得到该问题的数学模型如下:
min W = 103.7 l 11 + 83.6 l 12 + 56.7 l 13 + 46.7 l 21 + 36 l 22 + 29.1 l 23 s .t . l 11 + l 12 + l 13 = 500 l 21 + l 22 + l 23 = 650 0.001578 l 11 + 0.005182 l 12 + 0.02401 l 13 + 0.00176 l 21 + 0.00579 l 22 + 0.02683 l 23 ≤ 5 l i j ≥ 0
Step 6:优化计算
调用优化程序LP_simplex_general,
编写程序代码:
其中,C是目标函数的系数向量,Aeq是等式约束的系数矩阵,beq是等式约束对应的右端向量,A是不等式约束的系数矩阵,b是不等式约束对应的右端向量。
在MATLAB环境中运行上面的程序代码,得如下结果:
以上结果表明,重力输水管道系统优化设计线性规划模型的最优解为:
l 12 = 500 m , l 21 = 336 m , l 22 = 314 m
根据计算结果显示,管段1~2应选用管径0.4 m的管长500 m;管段2~3应选用管径为0.25 m的管长336 m;选用管径0.2 m的管长314 m。管道投资(目标函数)为68,796元。
大学生创新创业训练计划项目(2018)。
陈昊东,杜家康,孙洪春,任庆军. 单水源重力树状管网系统优化设计 Optimal Design of Single Source Gravity Tree Pipe Network System[J]. 运筹与模糊学, 2018, 08(02): 72-76. https://doi.org/10.12677/ORF.2018.82009