基于蚁群算法的环状燃气管网布置优化研究

摘 要

摘 要:将环状燃气管网布置优化问题转换为旅行商问题,采用蚁群算法进行求解,列举了算例。关键词:环状燃气管网 布置优化 旅行商问题 蚁群算法Layout Optimization of Loop Ga

摘 要:将环状燃气管网布置优化问题转换为旅行商问题,采用蚁群算法进行求解,列举了算例。

关键词:环状燃气管网  布置优化  旅行商问题  蚁群算法

Layout Optimization of Loop Gas Pipe Network Based on Ant Colony Algorithm

AbstractThe layout optimization problem of loop gas pipe network is converted into a traveling salesman problemand solution is performed by ant colony algorithmThe calculation examples are enumerated

Keywordsloop gas pipe networklayout optimizationtraveling salesman problemant colony algorithm

 

在设计燃气管网时,设计人员通常根据经验确定布置形式,这使得管网的合理性有待探讨。此外,国内外学者对枝状管网的布置优化问题研究较多,对环状管网的布置优化研究较少[1-5]。本文将环状燃气管网布置优化问题转换为旅行商问题(Traveling Salesman Problem,简称TSP),采用蚁群算法进行求解。

1 TSP

TSP是图论领域中著名问题之一。形象地说,TSP是指一名旅行售货员从一座城市出发,访问每座城市恰好一次,再回到出发城市,最终使得旅费最低。对于燃气管网,将调压站、调压装置视为TSP中的“城市(节点)”,节点之间敷设的燃气管道造价视为“旅费”,要求从起始节点出发,访问剩余的每个节点恰好一次,再回到起始节点,要求燃气管道造价最低。由此可知,环状燃气管网的布置优化问题与TSP十分相似,因此可以将环状燃气管网布置优

化问题转换为TSP

2 蚁群算法

蚁群算法(Ant Colony AlgorithmACA)1991MDorigo提出的一种智能算法。该算法来源于蚂蚁种群的觅食行为,其本质是一个庞大的智能体系,具有很强的容错性及并行计算能力,可与其他算法无缝衔接。

2.1 蚁群算法的数学模型

为了方便讨论,设bi(t)表示t时刻位于节点i的蚂蚁数量,则蚂蚁总数量m的表达式为:

 

式中m——蚂蚁总数量

n——节点总数

n个节点的集合为:

C{c1c2c3cn}

式中C——n个节点的集合

c1…cn——集合C中的元素

节点ij路径长度的集合为:

L{lijúcicjÌC}

 

式中L——集合C中节点ij路径长度的集合

lij——节点ij路径长度,m

xiyi——节点i在二维坐标系(xOy)中的坐标

xiyi——节点歹在二维坐标系(xOy)中的坐标

t时刻节点ij路径上残留信息量集合G的表达式为:

G{tij(t)ú cicjÌC}

式中G——t时刻节点ij路径上残留信息量集合

tij(t)——t时刻节点ij路径上的残留信息量,在初始时刻各条路径上残留信息量相等

蚂蚁k(123m)的运动方向是根据各条路径上的残留信息量来确定的,t时刻蚂蚁k由节点i转移到节点j的状态转移概率为:

 

式中Pij(t)——t时刻蚂蚁k由节点i转移到节点j的状态转移概率

a——信息启发式因子,表示轨迹的相对重要性

hij(t)——启发函数

b——期望启发式因子,表示能见度的相对重要性

A——蚂蚁k下一步允许选择的节点集合,A{C-B}B表示蚂蚁k已经走过的节点集合

需要注意的是,在每只蚂蚁走完一步或者完成对所有n个节点的遍历后,要对残留信息进行更新处理,否则过多的残留信息会淹没启发信息。因此,在蚂蚁经过一个节点或经过所有节点完成一个循环所需的时间T后,即t+T时刻在节点ij路径上的残留信息量tij(t+T)可按以下规则进行调整:

 

式中tij(t+T)——t+T时刻在节点ij路径上的残留信息量

r——信息素挥发系数,取值范围为[01)

Dtij(t)——本次循环中节点ij路径上的信息素增量,初始时刻为0

Dtij,k(t)——本次循环中第k只蚂蚁在节点ij路径上的信息素增量

MDorigo提出了3种不同的信息素更新模型,分别为:Ant-Cycle模型、Ant-Quantity模型、Ant-Density模型,差别在于对Dtij,k(t)的求解方法不同。在求解Dtij,k(t)时,Ant-Quantity模型、Ant-Density模型采用局部信息(即蚂蚁完成一步便更新当前路径上的信息素),而Ant-Cycle模型采用整体信息(即蚂蚁完成一个循环后才更新所有路径上的信息素)Ant-Cycle模型更适用于求解TSP,因此本文采用Ant-Cycle模型来更新信息素。

对于Ant-Cycle模型:

 

式中Q——信息素强度

Lk——第五只蚂蚁在本次循环中所走路径的总长度

2.2 参数确定

蚁群算法中的abrmQ等参数的设定尚无严格的理论依据,至今还没有确定最优参数的一般方法,已公布的参数设置成果都是针对不同蚁群算法模型的。当采用Ant-Cycle模型时,理想的参数设定范围为:0≤a50≤b≤50.1≤r≤0.9910≤Q≤10000

在面对具体问题时,可以按照段海滨[6]总结出的方法确定最优参数,具体步骤如下:第一步,确定m,按照

 

确定所需的蚂蚁数量。第二步,调整设定范围大的参数:abQ。第三步,调整设定范围小的参数:r。反复按照以上步骤进行调整,直至选择出一组比较满意的参数为止。

3 优化目标

引入当量费用长度概念,以环状燃气管网造价最低为目标,建立优化目标函数Lmin为:

 

式中Lmin——优化目标函数

lw,ij——当量费用长度,m

uij——判定系数

文献[7—8]认为燃气管网造价近似为管长的线性函数,基于这种理论,当量费用长度lw,ij的计算式为:

 

式中Wij——节点ij路径间管道的权系数

¦ij——燃气管道在实际敷设方式下节点ij路径间管道单位管长造价,元/m

¦——燃气管道在标准敷设方式下的单位管长造价,元/m

4 算例

选取某小型环状燃气管网(中压)作为优化对象,气源拟采用单气源供气。为了方便讨论,进行以下设定:所有的节点间管道的权系数相等,且为1;只考虑形成单环状管网。中压管网的初步布置形式见图l,节点4与气源连接,各节点在二维坐标系中的坐标见表1

 

 

利用C语言编制蚁群算法程序,经过多次试算确定蚁群算法的参数为:m15a1.0b5.0r0.1Q=100。迭代最大次数为200,防止出现死循环。优化后的环状管网布置形式见图2,燃气管网总长度为15.143km

 

参考文献:

[1]CHANG S KThe generation of minimal tree with a steiner topology[J]Journal of Association of Computing Machinery1972(4)699-711

[2]FLANIGN OConstrained derivatives in natural gas pipeline system optimization[J]Journal of Petroleum Technology1972(5)549-556

[3]吕木英.基于遗传算法的城市燃气管网最优化布局研究(硕士学位论文)[D].武汉:武汉理工大学,200956-59

[4]王垣,段常贵.改进遗传算法在燃气管网布局优化中的应用[J].哈尔滨工业大学学报,200638(1)46-48

[5]SALZBORN BOptimal design of gas pipeline networks[J]The Journal of the Operational Research Society1979(12)1047-1060

[6]段海滨.蚁群算法原理及其应用[M].北京:科学出版社,200534-36

[7]刘洪,胡昌权,胡攀峰,等.用遗传算法进行输气管网布置的优化设计研究[J].天然气工业,200626(10)127-129

[8]聂廷哲,段常贵.基于Hopfield神经网络的输气管网布线优化[J].天然气工业,200525(2)155-157

 

本文作者:李自力  赵峰  李佳  杨立业

作者单位:中国石油大学(华东)储运与建筑工程学院

  武汉炼化工程设计有限责任公司

  胜利油田胜利勘察设计研究院有限公司