基于二进制多目标粒子群算法的供水管网关阀方法与流程

本发明属于城市供水领域，具体是一种基于二进制多目标粒子群算法(BinaryMulti-ObjectiveParticleSwarmOptimization,BMOPSO)的供水管网关阀方法。
背景技术：
：爆管、漏损等异常情况在供水管网中普遍存在，爆管的不确定性、突发性、规模大以及漏损的持续时间长、不易察觉、规模小等特点会对日常工业生产、居民用水产生巨大影响。因此，如何找到一种快速、有效的关阀方法对供水管网异常情况下的处理以及平稳运行有着重要的意义和作用。然而实际的供水管网拓扑结构复杂，管道盘根交错，阀门的种类众多，因此确定快速、有效的关阀方案难度较大。目前，供水管网的关阀方法主要有两种：一种是调度员查询图纸、依靠经验确定关阀方案；另一种是依靠地理信息系统(GIS)在管网资料基础上确定关阀方案。前者速度慢、工作量大、准确性差，容易因为人为因素对管网关阀进行误判,导致关阀方案对实际情况不适用；后者由于许多常用的GIS软件没有专门建立网络拓扑机制，因此，关阀分析效率低，且欠准确。如果GIS软件专门建立网络拓扑机制，则需要耗费大量的人工进行校核。有鉴于此，本发明针对供水管网关阀方法的需要，提出了一种基于二进制多目标粒子群算法的供水管网关阀方法，减少了传统关阀的工作量，保证了关阀的快速性、有效性。技术实现要素：本发明的目的在于克服现有方法的不足，提出了一种基于二进制多目标粒子群算法的供水管网关阀方法，能够实现在拓扑结构供水管网的众多阀门中快速、有效的确定关阀方案，从而降低对管网的影响。本发明采取了以下步骤：(一)按异常位置划定搜索区域根据实际出现爆管、漏损等异常情况，在供水管网水力模型中根据二维坐标(横、纵坐标)确定异常位置，然后以异常位置为中心，按长度L宽度W得到搜索区域。L、W初始值须大于400m。(二)找出搜索区域内阀门根据划定的搜索区域内各个树状结构、环状结构，找出所有的阀门位置和个数。(1)树状结构根据关阀位置所在管段向管段两端查找阀门，然后再通过树状的分支查找下一管段上的阀门，以此类推，直到找出搜索区域内所有分支上的阀门为止。(2)环状结构根据关阀位置所在管段向环状管段两侧查找阀门，依次遍数阀门的个数，直到在搜索区域内查找的路径完全闭合。(三)确立关阀BMOPSO优化模型根据步骤(二)找出搜索区域内所有阀门，令每个阀门对应一个粒子，阀门的位置对应粒子的位置。(1)初始化种群及编码种群是粒子的集合，初始化种群Pt＝{p1,p2,…,pn}，对种群中粒子p位置x＝{x1,x2,…,xi,…,xm}进行二进制编码，其中1表示阀门状态为关闭，反之0表示开启，位置xi编码方式如下：其中：n为种群数，t为种群代数，m为位置数，rand为0～1之间的随机数。(2)确定粒子初始值分别确定粒子初始速度、粒子初始最好位置p_best，粒子初始非劣解g_best。粒子初始速度一般设为0的向量，粒子初始最好位置p_best、粒子初始非劣解g_best由适应度函数f(x)得到，其中Ci＝f(xi)，Ci表示粒子p位置xi对应的适应度。适应度函数包括两部分f＝[f1f2]。其中，f1表示关阀位置流量flow，f2表示粒子位置xi的和。即：minf1＝flow(3)粒子速度和位置更新更新每个粒子的速度、位置以及最好位置p_best，得到更新后的Pt+1。具体公式如下：vi(t+1)＝wvi(t)+c1rand×(p_besti-xi)+c2rand×(p_besti-xi)w＝w×wd其中：w和wd分别为惯性常数和惯性常数衰减率，c1和c2为正的常数。为更好的保持粒子的多样性，BMOPSO采用自适应变异，变异公式如下：其中：l表示当前迭代数，S表示当前迭代次数，ρ表示当前迭代次数的变异概率，μ表示变异概率为常数。(四)求解BMOPSO优化模型，生成关阀方案根据(三)中对粒子位置和速度的更新进行迭代，将粒子位置对应的阀门状态在水力模型中的阀门开关状态进行设置(1表示阀门状态为关闭，反之0表示开启)，并进行单步离线模拟，然后判断是否达到给定的迭代次数，若未达到，继续更新粒子速度和位置，给出新的非劣解再进行单步离线模拟；若达到迭代次数得到的非劣解后，根据目标函数(适应度函数)计算非劣解中关阀位置流量以及关阀数量和，设置流量最小情况下关阀数量和最小为优选条件，得到最终关阀位置。若根据目标函数(适应度函数)值无法确定关阀位置，则扩大搜索长度L和宽度W，重新划定搜索区域，重新执行步骤(二)，输出关阀方案。本发明有益效果：本发明供水管网关阀方法能在异常情况下快速、有效的确定供水管网关阀方案以隔离异常点，相比于调度员查询图纸、依靠经验确定关阀方案，提高了关阀速度和准确性，且避免GIS软件专门针对关阀而建立网络拓扑机制所耗费大量的人工校核。附图说明图1：本发明方法流程图；图2：本发明实施例待关阀位置分布图；图3：本发明实施例非劣解分布图。具体实施方式为使本发明实现的技术手段与创作特征易于明白，下面结合附图和实施例，对本发明的实现方式进一步详述，并不限制本发明的权利范围。本实施例以S市核心城区供水管网模型为例，为使实例更具普遍性，通过选取供水管网模型中四根管径不同的管段作为异常位置，异常位置在管网中位置见图二中的空白方块区，具体位置信息见表一所示。表一异常位置信息实验编号管道ID模拟位置横坐标(m)模拟位置纵坐标(m)管径(mm)175432758817.0819848.951000278783655448.7423191.67800375130957636.7619395.64600478826656177.3020216.28400(一)按异常位置划定搜索区域根据表一异常位置的信息，结合异常位置在供水管网模型中的坐标分布，为具有普遍性，四个实验根据异常率变化明显的测点分别划定出大小不同的搜索范围，实验一划定搜索范围为1000×1000的矩形(正方形)，实验二划定搜索范围为1200×400的矩形，实验三划定搜索范围为600×800的矩形，实验四划定搜索范围为400×1200的矩形。具体坐标范围见表二。表二搜索范围坐标(二)找出搜索区域内阀门结合供水管网模型以及EPANET软件，在相应的搜索范围内，根据树状结构或环状结构在搜索区域内的分布情况，获取符合范围坐标的阀门的位置以及阀门的数量，其中在实验1搜索区域内阀门数量为204个，实验2搜索区域内阀门数量为71个，实验3搜索区域内阀门数量为623个，实验4搜索区域内阀门数量为652个。(三)确立关阀BMOPSO优化模型(1)初始化种群及编码根据(二)在各个搜索区域内的阀门数量，实验1中的204个阀门作为一个种群，实验2中71个阀门作为一个种群，实验3中623个阀门作为一个种群，实验4中652个阀门作为一个种群。对每个种群内阀门的位置xi进行编码，编码方式如下：说明：n为种群数，t为种群代数，m为位置数，rand为0～1之间的随机数。其中：种群数n＝300,种群代数t＝100，位置数m为每个实验搜索出的阀门位置数。(2)确定粒子初始值粒子初始速度为向量，为是后续的速度易于表示，初始速度定为0的向量。粒子初始最好位置p_best、粒子初始非劣解g_best由适应度函数f(x)得到，其中Ci＝f(xi)，Ci表示粒子p位置xi对应的适应度。适应度函数包括两部分f＝[f1f2]。其中，f1表示关阀位置流量flow，f2表示粒子位置xi的和。即：minf1＝flow其中：初始最好位置p_best，粒子初始非劣解g_best见图三黑色圆圈分布。(3)粒子速度和位置更新根据(2)对粒子初始值设定，为达到种群数，对粒子位置和速度以及最好位置p_best进行更新。公式如下：vi(t+1)＝wvi(t)+c1rand×(p_besti-xi)+c2rand×(p_besti-xi)w＝w×wd说明：w和wd分别为惯性常数和惯性常数衰减率，c1和c2为正的常数。为更好的保持粒子的多样性，BMOPSO采用自适应变异，变异公式如下：其中，l表示当前迭代数，S表示当前迭代次数，ρ表示当前迭代次数的变异概率，μ表示变异概率为常数。这里，迭代次数S＝200，惯性常数w＝0.5和惯性常数衰减率wd＝0.9，常数c1＝1和c2＝2，自适应变异概率μ＝0.1，非劣解最大个数为Np＝100。(四)求解BMOPSO优化模型，生成关阀方案根据(三)中对粒子位置和速度的更新进行迭代，将粒子位置对应的阀门状态在水力模型中进行设置(1表示阀门状态为关闭，反之0表示开启)，并通过进行单步离线模拟得到关阀位置的流量，达到迭代次数后，得到四组实验的非劣解(如图3)，图中横轴表示关阀位置流量值f1，纵轴表示关闭阀门个数f2，圆圈黑色表示所有样本的适应度值，星型标记表示非劣解的适应度值。由于采用的粒子群算法为多目标算法，得到的最终解并不是唯一的最优解而是得到非劣解，所以根据目标函数(综合适应度)设置流量最小前提下粒子位置数量和最小为优选条件，得到最终的关阀方案见表三。表三关阀方案由此，本发明专利解决了以往关阀方法的工作量大、速度慢、关阀准确性差等问题，并且本发明方法不仅适用于本实施例，在供水管网关阀方法上具有普遍性。当前第1页1 2 3