1. 引言
在村镇供水以及工矿企业给水系统的建设初期,一般都采用树状管网,因为树状给水管网一般投资比较小,以后会随着用水规模的扩大,由树状管网逐渐发展为环状管网。树状管网形式较多,不同的形式有不同的特点,也会满足不同的情形。本文主要对水源位置较高时,采用的重力树状管网系统进行优化。单水源重力树状管网系统要解决的问题是:在保证各节点所需流量和服务水压的情况下,确定管网中各个管段尺寸的最优值,使得整个管网投资最少。目前,对管网系统的优化已经引起了专家和学者的重视。王文芬 [1] 提出将管网优化问题转化成管网投资线性子模型问题,并采用大系统分解-动态聚合方法对模型进行求解,为泵站加压式树状供水管网系统优化设计提供了一种可行而有效的新方法;许晓明 [2] 以油气集输管网系统投资极小化为目标,综合考虑多级星式管网拓扑结构特点、生产运行参数等多种影响因素,建立了多级星式油气集输系统模糊优化数学模型;郑凤 [3] 建立了天然气管网系统输配气运行方案优化的数学模型,采用线性逼近法对模型进行了求解,应用于实际管网系统运行调度分析。本文主要对单水源重力树状管网系统进行了优化设计,建立了线性规划数学模型,应用MATLAB编制优化设计计算程序。
2. 优化设计数学模型
首先,我们先对涉及的一些术语进行解释说明。在给水管网系统中,节点水压标高是指管道上节点压力水头与该点地面高程之和。给水管网系统中具有一定的水压,用来保证用户的用水,我们将这个水压称为最小服务水压。服务水压标高就是指管道上节点最小服务水压与该点地面高程之和 [4] 。地面坡度、管段流量和管段承压力都相同的管道称为管段。
一般来说,工程上采用标准管径,也就是说每一管段可以由一种或多种标准管径组合成。由于标准管径档次太多,计算量太大,为方便计算,可以通过流速来限制管径。为防止水击破坏,管道内的流速应该限制不高于3 m/s;为防止淤积,管道内的流速也不应该低于0.6 m/s [5] 。管道流速(m/s)用Vij表示,i表示管段序号,j表示备选标准管径序号,
(1)
(2)
其中
表示管段设计流量(m3/s),dj表示管道计算内径(m)。
在第i管段中,可以根据(1)式的条件来选择标准管径。若这样的标准管径档次有
个,则这
个标准管径组成第i管段的备选管径组。
以单水源重力树状管网投资最小为目标函数
(3)
其中W表示重力树状管网投资(元),n表示管段数,
表示第i管段的备选管径数,
表示第i管段中的第j种标准管径管道单价(元/m),
表示第i管段中选用的第j种标准管径的管长(m)。
每一管段各标准管径管长之和等于该管段长,也就是
(4)
其中
表示管段长(m)。要求各个节点水压标高均不得低于该节点服务水压标高,即
(5)
其中r表示节点序号,R表示节点总数,
表示水源至第r点的管段数,
表示水源水面高程(m),
表示节点服务水压标高(m),pij表示管道水力坡度,对于混凝土管及钢筋混凝土管,可用巴甫洛夫斯基公式 [6]
(6)
其中
表示水力半径(m),
,流速系数。
(7)
其中N表示粗糙系数,混凝土管及钢筋混凝土管一般取0.013~0.014。另外,
。
由以上分析,可建立单水源重力树状管网系统的数学模型
3. 数值实验
这是一个带约束的线性规划模型,优化变量为lij,采用单纯形法求解,下面根据例子对数学模型进行进一步的解释。单水源重力给水树状管网布置如图1所示。采用钢筋混凝土管,管道单价见表1,管道的粗糙系数
,当管道起点与终点高差为5 m (
),确定单水源重力树状管网优化设计方案 [7] 。
Step 1:划分管段
Step 2:确定各管段备选管径组
根据式(2)确定各个管段的备选管径组见表2。
Step 3:计算水力坡度
根据巴氏公式(6),计算各管段水力坡度见表2。
Figure 1. Single-source gravity tree pipe network calculation diagram
图1. 单水源重力树状管网计算图
Table 2. Pipeline alternative diameter group and its hydraulic gradient
表2. 管段备选管径组及其水力坡度
Step 4:确定目标函数
根据式(3)和表1,确定目标函数如下:
Step 5:确定约束条件
根据式(4)和图1,确定管长约束条件如下:
根据式(5)和表2,确定压力约束条件如下:
非负约束条件如下:
。
通过整理以上式子,得到该问题的数学模型如下:
Step 6:优化计算
调用优化程序LP_simplex_general,
编写程序代码:
其中,C是目标函数的系数向量,Aeq是等式约束的系数矩阵,beq是等式约束对应的右端向量,A是不等式约束的系数矩阵,b是不等式约束对应的右端向量。
在MATLAB环境中运行上面的程序代码,得如下结果:
以上结果表明,重力输水管道系统优化设计线性规划模型的最优解为:
根据计算结果显示,管段1~2应选用管径0.4 m的管长500 m;管段2~3应选用管径为0.25 m的管长336 m;选用管径0.2 m的管长314 m。管道投资(目标函数)为68,796元。
基金项目
大学生创新创业训练计划项目(2018)。