一种基于广义逆向学习差分算法的排水管网优化方法技术

本发明专利技术涉及一种基于广义逆向学习差分算法的排水管网优化方法。本发明专利技术针对排水管网优化过程中存在的计算复杂、且不能彻底排除不合理或者不可行方案，计算精度不高的问题，利用改进的自适应差分算法，得到管网造价函数最优值，包括满足水力约束条件的各管段的流速、坡度、埋深等水力参数。本发明专利技术采用逆向学习方法进行自适应算法的初始化过程，对算法的初始化种群在进行变异之前进行处理，利用概率论的知识，在不增加种群个数的前提下精炼了初始解，增加了种群的多样性，提高了找到全局最优解的可能性。且该算法涉及的变量较少，自适应选择变异策略，算法操作简单，收敛速度快，寻优精度高。

【技术实现步骤摘要】
一种基于广义逆向学习差分算法的排水管网优化方法
本专利技术属于排水管网优化领域，涉及一种基于广义逆向学习差分算法的排水管网优化方法，采用广义逆向学习方法来改进自适应差分算法，用于排水管网的造价优化设计。
技术介绍
排水系统通常由排水管网和污水处理厂组成，排水管网是收集和输送生活污水、工业废水和雨雪水的设施，并负责将这些废弃物分别输送到污水厂或其他出水口。作为城市基础设施建设中不可或缺的部分，系统主要由排水管道、排水泵站、检查井等工程设施组成。为了满足人民的物质生活需要，进一步改善生活环境，还将建设或改造大量的排水管道。于是，在满足规定的水力约束条件下，降低排水管道的基建费用成为一个重要研究方向。从上世纪60年代起，国际上的学者结合数理分析和经验总结，建立了许多排水管网系统的数学模型，将实际中的复杂管网在满足设计规范和管径、流速、坡度、充满度之间的水力关系的前提下，去掉次要管道，抽象得到由节点和管线组成的管网模型，以此作为研究基准，进行排水管网优化设计。由《室外排水规范》GB50014-2006规定，排水管网是以造价最低为目标，在均匀流公式下计算，满足管径、流速、埋深、衔接方式等条件下，以充满度和流速为决策性变量，按照污水在管道中流动的方向，采用先计算支管后计算干管、主干管，从上游到下游依次对管段进行计算，进而完成整个排水管网的计算。在管线布置形式给定条件下，其优化设计一般包括不同管径、坡度两个相互关联的方面组合优化。对于某个具体设计管段，设计流量确定以后，满足要求的管径和埋深通常有很多组，如果选择较大管径，管材费用就较高，但是坡度可以很小，管道埋深较浅，施工费用则较低；相反，如果管径选择较小，施工费用高而管材费用低。因此，对于每一段管道，总存在一组较优的管径和埋深组合，使管网投资最小。目前在管网优化中主要采用的方法有线性规划法、非线性规划法、动态规划法、电子表格法、两相优化法等，但是这些方法需要将所有的目标函数和约束条件严格线性化，因此在管段设计的优化过程中，无法完全排除不合理或不可行方案，并且计算时间较长，对计算机内存要求高，参数设置复杂，易陷入局部最优解，计算精度不能保证。
技术实现思路
针对在排水管网优化中存在的预处理工作量大、求解过程复杂等问题，本专利技术提出一种基于广义逆向学习差分算法的排水管网优化方法。利用广义的逆向学习方法进行多策略自适应差分算法的初始化策略调整，求出各个候选解的相应逆向点，并在候选解和其逆向点中选择所需要的最优初始种群，然后再进行自适应变异、杂交、选择操作。具有收敛速度快、寻优精度高、参数设置简单的优点。根据概率理论，一个点有50％的可能性能比它相应的逆向点获得更好的适应值，将广义逆向学习方法应用到自适应差分算法中，可以有效利用群体和逆向群体的信息，提高对种群的原搜索空间的利用能力，更加方便高效的解决排水管网优化问题。本专利技术的目的是选择合适的敷设方式降低排水管网的造价。排水管网作为一种复杂系统，设计的约束变量较多，主要有设计充满度、流速、最小管径、最小设计坡度、埋设深度和管道衔接方式等。且这些约束条件是相互联系与制约的：在流量确定后，管径若增加，坡度、充满度则减少；相同管径时，坡度若减少，充满度则增加；相同坡度时，管径若增加，充满度则减少。解决此类问题的决策变量通常有三种选择方式：坡度，管径，还有将坡度和管径同时作为决策变量。如果采用坡度作为决策变量，虽然可以直接选择标准管径，优化结果与初始解无关，便于控制优化精度，但要求设置状态点的埋深间隔很小，造成存储量过大，使计算时间过大。所以，在这里选用管径为决策变量，优化设计决定各个管段的最佳管径后，通过水力计算公式求出其他所需的变量。一般只需确定出管径和埋深，由排水管道造价指标公式即可求出所需的费用。一种基于广义逆向学习差分算法的排水管网优化方法，包括以下步骤：步骤1，编码与初始化。步骤1.1，采用实数编码，种群P中有NP个个体，每个个体代表管网的一种敷设方式；自变量有D维，即有D条管段需要优化，则群体中的第G代第i个个体Xi,G可表示为：其中，i＝1,2,…,NP，j＝1,2,…D，为一实数，代表管道的管径的大小，且在自变量范围[a,b]内均匀初始化；步骤1.2，计算点Xi,G的逆向点：GOPi,G＝k[a+b]-Xi,G，组成NP个个体的种群GOP，其中，k＝rand(0,1)，a,b分别表示自变量的取值范围的上下界。判断个体是否在边界[a,b]内。如果不在边界内，采用截断的方法使GOPi,G在边界内，即GOPi,G＝rand(a,b)，以管网敷设造价为目标函数，计算(P,GOP)中每一个个体的适应值，即目标函数的值，并挑选适应值最优的NP个个体；步骤1.3，判断是否满足rand(0,1)≤p0，p0为概率值，通常取0.05。如果满足，执行1.2更新当前种群；否则，执行下一步。步骤2，变异操作。差分算法中最重要的部分是的差分变异算子的选用。用“DE/c/d”的形式区别不同的算子，其中“DE”表示差分算法；“c”代表基向量的选择方式，一般有rand和best两种；“d”表示算子中差分向量的个数。比较常用的变异策略有：DE/rand/1、DE/best/1、DE/rand-to-best/1、DE/best/2、DE/rand/2，产生变异向量通过控制变异概率F完成对差分量的放大和缩小，由此来控制算法的搜索步长的大小。设置的越大，变异向量对后代个体的影响越大，算法收敛越慢。但是过小，又会使种群多样性低，容易出现“早熟”收敛现象。在解决不同的最优化问题时，实验向量采用不同的DE策略进行迭代将产生不同的效果。将不同特性的策略放在一起形成一个候选策略库，策略库中包含四种变异策略：DE/rand/1/bin，DE/rand-to-best/2/bin,DE/rand/2/bin和DE/current-to-rand/1，其中，DE/rand/1/bin,DE/rand-to-best/2/bin在很多DE算法的改进工作中频繁使用，是比较经典的两种变异策略；DE/rand/2/bin由于加上了高斯扰动使算法具有很好的探测能力；DE/current-to-rand/1可以更有效地解决多目标函数优化问题。四种变异策略的公式为：DE/rand/1/bin：Vi,G＝Xr1,G+F(Xr2,G-Xr3,G)DE/rand-to-best/2/bin：Vi,G＝Xr1,G+F(Xbest,G-Xi,G)+F(Xr1,G-Xr2,G)+F(Xr3,G-Xr4,G)DE/rand/2/bin：Vi,G＝Xr1,G+F(Xr2,G-Xr3,G)+F(Xr4,G-Xr5,G)DE/current-to-rand/1：Vi,G＝Xi,G+K(Xr1,G-Xi,G)+F(Xr2,G-Xr3,G)其中，每个个体i的变异概率F＝rndn(0.5,0.3)为均值为0.5、方差为0.3的正态分布随机数，可以使变异步长和搜索方向根据管网造价作自适应调整，提高了算法的精度。K＝rand(0,1)为一随机数，Xbest为当前群体的最优个体，Xi为父个体，下标r1，r2，r3，r4，r5为群体中随机选择的5个不同于Xi的个体。在迭代的过程中，根据一定概率pk从策略库中选本文档来自技高网...

【技术保护点】
一种基于广义逆向学习差分算法的排水管网优化方法，其特征在于包括以下步骤：步骤1，编码与初始化；步骤1.1，采用实数编码，种群P中有NP个个体，每个个体代表管网的一种敷设方式；自变量有D维，即有D条管段需要优化，则群体中的第G代第i个个体Xi,G可表示为其中，i=1,2,…,NP，j=1,2,…D，为一实数，代表管道的管径的大小，且在自变量范围[a,b]内均匀初始化；步骤1.2，计算点Xi,G的逆向点GOPi,G＝k[a+b]‑Xi,G，组成NP个个体的种群GOP，其中，k=rand（0,1），a、b分别表示自变量的取值范围的上下界；判断个体是否在边界[a,b]内；如果不在边界内，采用截断的方法使GOPi,G在边界内，即GOPi,G＝rand(a,b)；以管网敷设造价为目标函数，计算(P,GOP)中每一个个体的适应值，即目标函数的值，并挑选适应值最优的NP个个体；步骤1.3，判断是否满足rand(0,1)≤p0，p0为概率值，通常取0.05；如果满足，执行1.2更新当前种群；否则，执行下一步；步骤2，变异操作；用“DE/a/b”的形式区别不同的算子，其中“DE”表示差分算法；“a”代表基向量的选择方式，一般有rand和best两种；“b”表示算子中差分向量的个数；比较常用的变异策略有：DE/rand/1、DE/best/1、DE/rand‑to‑best/1、DE/best/2、DE/rand/2，产生变异向量通过控制变异概率F完成对差分量的放大和缩小，由此来控制算法的搜索步长的大小；设置的越大，变异向量对后代个体的影响越大，算法收敛越慢；但是过小，又会使种群多样性低，容易出现“早熟”收敛现象；在解决不同的最优化问题时，实验向量采用不同的DE策略进行迭代将产生不同的效果；将不同特性的策略放在一起形成一个候选策略库，策略库中包含四种变异策略：DE/rand/1/bin，DE/rand‑to‑best/2/bin，DE/rand/2/bin和DE/current‑to‑rand/1；四种变异策略的公式为：DE/rand/1/bin：Vi,G＝Xr1,G+F(Xr2,G‑Xr3,G)DE/rand‑to‑best/2/bin：Vi,G＝Xr1,G+F(Xbest,G‑Xi,G)+F(Xr1,G‑Xr2,G)+F(Xr3,G‑Xr4,G)DE/rand/2/bin：Vi,G＝Xr1,G+F(Xr2,G‑Xr3,G)+F(Xr4,G‑Xr5,G)DE/current‑to‑rand/1：Vi,G＝Xi,G+K(Xr1,G‑Xi,G)+F(Xr2,G‑Xr3,G)其中，每个个体i的变异概率F＝rndn(0.5,0.3)为均值为0.5、方差为0.3的正态分布随机数，可以使变异步长和搜索方向根据管网造价作自适应调整，提高了算法的精度；K=rand(0,1)为一随机数，Xbest为当前群体的最优个体，Xi为父个体，下标r1，r2，r3，r4，r5为群体中随机选择的5个不同于Xi的个体；在迭代的过程中，根据一定概率pk从策略库中选择一种变异策略，在之前的迭代中越容易达到最优解的DE策略，越容易在此迭代操作中被选择；算法参数的设置如下：选择概率pk初始化为1/K，采用随机广义选择的方式对每一个目标向量选择策略，经过LP次迭代演化以后，pk按下式进行更新：pk=Sk,GΣk=1KSk,G]]>Sk,G=Σg=G-LPG-1nsk,gΣg=G-LPG-1nsk,g+Σg=G-LPG-1nfk,g+ϵ]]>其中：G＞LP是当前的演化代数；nsk,g和nfk,g分别是过去LP次迭代中第k个策略所产生子个体成功或者失败进入下一代群体的数量；Sk,g为第k个策略所产生子个体成功进入下一代群体的成功率；ε是一个很小的常数，可避免出现零值；步骤3：杂交操作；此步操作是为了增强新种群的离散程度，设置合适的CR值可以加快算法的收敛，尽快找出管网的最优造价，提高了算法的收敛速度；通过变异算子产生的变异向量Vi,G与父向量Xi,G进行离散杂交操作得到试验向量通过二项式杂交算子得到个体的第j维个体可标记为：ui,Gj=vi,Gj,if(rand[0,1)≤CR)or(j=jrand)xi,Gj,otherwise]]>其中，j=1,…,D，CR为交叉概率因子，jrand是介于维数[1,D]之间的一个均匀分布随机整数，用来保证试验向量Ui,G中至少有一维来自变异向量Vi,G，避免与父个体向量Xi,G相同；种群中每个个体选择第k个策略时的杂交概率CRi,k＝rndni(CRmk,0.1)，如其值不在[0,1]范围内，则采用截断方法限制在这一范围内，CRmk初始为0.5，CRMemoryk为前LP代中第k个策略产生的...

【技术特征摘要】
1.一种基于广义逆向学习差分算法的排水管网优化方法，其特征在于包括以下步骤：步骤1，编码与初始化；步骤1.1，采用实数编码，种群P中有NP个个体，每个个体代表管网的一种敷设方式；自变量有D维，即有D条管段需要优化，则群体中的第G代第i个个体Xi,G可表示为其中，i＝1,2,…,NP，j＝1,2,…D，为一实数，代表管道的管径的大小，且在自变量范围[a,b]内均匀初始化；步骤1.2，计算点Xi,G的逆向点GOPi,G＝k[a+b]-Xi,G，组成NP个个体的种群GOP，其中，k＝rand(0,1)，a、b分别表示自变量的取值范围的上下界；判断个体是否在边界[a,b]内；如果不在边界内，采用截断的方法使GOPi,G在边界内，即GOPi,G＝rand(a,b)；以管网敷设造价为目标函数，计算(P,GOP)中每一个个体的适应值，即目标函数的值，并挑选适应值最优的NP个个体；步骤1.3，判断是否满足rand(0,1)≤p0，p0为概率值，通常取0.05；如果满足，执行1.2更新当前种群；否则，执行下一步；步骤2，变异操作；用“DE/c/d”的形式区别不同的算子，其中“DE”表示差分算法；“c”代表基向量的选择方式，有rand和best两种；“d”表示算子中差分向量的个数；常用的变异策略有：DE/rand/1、DE/best/1、DE/rand-to-best/1、DE/best/2、DE/rand/2，产生变异向量通过控制变异概率F完成对差分量的放大和缩小，由此来控制算法的搜索步长的大小；设置的越大，变异向量对后代个体的影响越大，算法收敛越慢；但是过小，又会使种群多样性低，容易出现“早熟”收敛现象；在解决不同的最优化问题时，实验向量采用不同的DE策略进行迭代将产生不同的效果；将不同特性的策略放在一起形成一个候选策略库，策略库中包含四种变异策略：DE/rand/1/bin，DE/rand-to-best/2/bin，DE/rand/2/bin和DE/current-to-rand/1；四种变异策略的公式为：DE/rand/1/bin：Vi,G＝Xr1,G+F(Xr2,G-Xr3,G)DE/rand-to-best/2/bin：Vi,G＝Xr1,G+F(Xbest,G-Xi,G)+F(Xr1,G-Xr2,G)+F(Xr3,G-Xr4,G)DE/rand/2/bin：Vi,G＝Xr1,G+F(Xr2,G-Xr3,G)+F(Xr4,G-Xr5,G)DE/current-to-rand/1：Vi,G＝Xi,G+K(Xr1,G-Xi,G)+F(Xr2,G-Xr3,G)其中，每个个体i的变异概率F＝rndn(0.5,0.3)为均值为0.5、方差为0.3的正态分布随机数，可以使变异步长和搜索方向根据管网造价作自适应调整，提高了算法的精度；K＝...

【专利技术属性】
技术研发人员：乔俊飞，刘昌芬，韩红桂，武利，王超，李瑞祥，
申请(专利权)人：北京工业大学，
类型：发明
国别省市：北京;11