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

本发明专利技术公开了一种基于二进制多目标粒子群算法的供水管网关阀方法。本发明专利技术首先按异常位置划定搜索区域。其次找出搜索区域内阀门，包括所有的阀门位置和个数。然后确立二进制多目标粒子群算法优化模型，令每个阀门对应一个粒子，阀门的位置对应粒子的位置，进行初始化种群及编码、确定粒子初始值和粒子速度和位置更新。最后求解上述优化模型，生成关阀方案。本发明专利技术能在异常情况下快速、有效的确定供水管网关阀方案以隔离异常点，相比于调度员查询图纸、依靠经验确定关阀方案，提高了关阀速度和准确性，且避免GIS软件专门针对关阀而建立网络拓扑机制所耗费大量的人工校核。

【技术实现步骤摘要】
基于二进制多目标粒子群算法的供水管网关阀方法
本专利技术属于城市供水领域，具体是一种基于二进制多目标粒子群算法(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_b本文档来自技高网...

【技术保护点】
基于二进制多目标粒子群算法的供水管网关阀方法，其特征在于该方法具体是：步骤(一)按异常位置划定搜索区域：根据实际出现爆管、漏损的异常情况，在供水管网水力模型中根据二维坐标确定异常位置，然后以异常位置为中心，按得到搜索区域；步骤(二)找出搜索区域内阀门：根据划定的搜索区域内各个树状结构、环状结构，找出所有的阀门位置和个数；步骤(三)确立二进制多目标粒子群算法优化模型：根据步骤(二)找出搜索区域内所有阀门，令每个阀门对应一个粒子，阀门的位置对应粒子的位置，执行以下过程：初始化种群及编码；确定粒子初始值；粒子速度和位置更新；步骤(四)求解上述优化模型，生成关阀方案：根据步骤(三)中对粒子位置和速度的更新进行迭代，将粒子位置对应的阀门状态在水力模型中的阀门开关状态进行设置，并进行单步离线模拟；判断是否达到给定的迭代次数，若未达到，继续更新粒子速度和位置，给出新的非劣解再进行单步离线模拟；若达到迭代次数得到的非劣解后，根据目标函数计算非劣解中关阀位置流量以及关阀数量和，设置流量最小情况下关阀数量和最小为优选条件，得到最终关阀位置。

【技术特征摘要】
1.基于二进制多目标粒子群算法的供水管网关阀方法，其特征在于该方法具体是：步骤(一)按异常位置划定搜索区域：根据实际出现爆管、漏损的异常情况，在供水管网水力模型中根据二维坐标确定异常位置，然后以异常位置为中心，按得到搜索区域；步骤(二)找出搜索区域内阀门：根据划定的搜索区域内各个树状结构、环状结构，找出所有的阀门位置和个数；步骤(三)确立二进制多目标粒子群算法优化模型：根据步骤(二)找出搜索区域内所有阀门，令每个阀门对应一个粒子，阀门的位置对应粒子的位置，执行以下过程：初始化种群及编码；确定粒子初始值；粒子速度和位置更新；步骤(四)求解上述优化模型，生成关阀方案：根据步骤(三)中对粒子位置和速度的更新进行迭代，将粒子位置对应的阀门状态在水力模型中的阀门开关状态进行设置，并进行单步离线模拟；判断是否达到给定的迭代次数，若未达到，继续更新粒子速度和位置，给出新的非劣解再进行单步离线模拟；若达到迭代次数得到的非劣解后，根据目标函数计算非劣解中关阀位置流量以及关阀数量和，设置流量最小情况下关阀数量和最小为优选条件，得到最终关阀位置。2.根据权利要1所述的基于二进制多目标粒子群算法的供水管网关阀方法，其特征在于：所述的搜索区域为矩形，其长度、宽度至少大于400米。3.根据权利要1所述的基于二进制多目标粒子群算法的供水管网关阀方法，其特征在于：树状结构搜索阀门的方式为：根据关阀位置所在管段向管段两端查找阀门，然后再通过树状的分支查找下一管段上的阀门，以此类推，直到找出搜索区域内所有分支上的阀门为止。4.根据权利要1所述的基于二进制多目标粒子群算法的供水管网关阀方法，其特征在于：环状结构搜索阀门的方式为：根据关阀位置所在管段向环状管段两侧查找阀门，依次...

【专利技术属性】
技术研发人员：徐哲，张慧，熊晓锋，何必仕，
申请(专利权)人：杭州电子科技大学，
类型：发明
国别省市：浙江,33