主权项 |
一种多变电站设备运行维护优化方法,包括以下步骤:步骤1:输入变量信息素的计算,依据变电站的缺陷分布、停电时间分布、运维效果量化、运维工作量量化四个输入量进行信息素的量化;<img file="dest_path_FDA0000666451130000011.GIF" wi="1091" he="381" />式一式一中:D(i,l,k,m)表示每个变电站在特定月份l序号为i的信息素,Data(i,l,q,p,t)代表对应i序号的输入信息度量,即缺陷分布、停电时间分布、运维效果量化、运维工作量的度量;i为量化的输入信息内容序号,缺陷分布、停电时间分布、运维效果量化、运维工作量分别为1,2,3,4;r为黄金分割数0.618;l为月份的序号,取值范围1~12;m为变电站序号;q为变电站数量;n为设备当前运行的年份;式一左半部分为该变电站本身的输入变量信息素特征,而右边为所有变电站该输入变量的信息素特征;步骤2:根据步骤1的结果,构建根据月份的变电站信息素<img file="dest_path_FDA0000666451130000012.GIF" wi="474" he="145" />式二式二中D<sub>S</sub>(l,m)为变电站的维护信息素;步骤3:路径信息素<img file="dest_path_FDA0000666451130000013.GIF" wi="792" he="141" />式三式三中,D<sub>R</sub>(l,Idx)为路径的信息素;l的含义参见式一,Idx为线路的序号;Time(l,Idx)为路线Idx在工作时间内的平均通过时间;TimeMin(l)以及TimeMax(l)分别为所有路径中最小以及最大的平均通过时间;步骤4:采用蚁群算法进行信息素的更新;将所有变电站与路径进行编号,则(c,d)代表从c点到d点的路径:<img file="dest_path_FDA0000666451130000021.GIF" wi="688" he="88" />式四其中,路径(c,d)在时刻h的信息素轨迹强度为τ<sub>cd</sub>(h),而蚂蚁e在该路径上留下的单位信息素数量<img file="dest_path_FDA0000666451130000022.GIF" wi="172" he="77" />的持久性为ρ(0≤ρ<1);h为迭代次数;<img file="dest_path_FDA0000666451130000023.GIF" wi="174" he="75" />的初始值为D<sub>R</sub>(l,Idx);τ<sub>cd</sub>(h)的初始值为(D<sub>S</sub>(l,c)+D<sub>S</sub>(l,d))/2;蚁群蚂蚁的个数为E,与运维队伍个数相同;对应的蚂蚁转移概率为:<img file="dest_path_FDA0000666451130000024.GIF" wi="682" he="271" />式五其中,U为可行项顶点集,α及β分别表示蚂蚁在运动过程中所积累的信息素及启发式因子的权重因子,均设置为1;而η<sub>cd</sub>表示选择元素d的期望程度,取值为待选路径信息素的倒数,以正态分布函数为基础设定ρ值,如式六,<img file="dest_path_FDA0000666451130000025.GIF" wi="626" he="173" />式六其中x为正整数,σ为变电站节点的个数,而<img file="dest_path_FDA0000666451130000026.GIF" wi="154" he="71" />是ρ值取0.382时x的界限;步骤5:遗传算法对于交叉算子Pc,其典型取值为0.6~0.95,在本应用中采用稳态选择方式,见式七,P<sub>c</sub>=0.6+ρ×(0.95‑0.6) 式七=0.6+0.35×ρ而对于变异操作,变异率Pm通常取值为0.001~0.01,式八给出变异算子,P<sub>m</sub>=0.01‑0.009×ρ 式八遗传算法的目标函数定义为:<img file="dest_path_FDA0000666451130000031.GIF" wi="679" he="107" />式九约束函数为:<img file="dest_path_FDA0000666451130000032.GIF" wi="530" he="106" />式十步骤6:算法边界从蚁群算法到遗传算法存在一个优化边界,这里采用ρ值作为判据,将整个优化过程分为3个阶段,见式十一,<img file="dest_path_FDA0000666451130000033.GIF" wi="555" he="191" />式十一第一阶段和第二阶段的分界点划分蚁群算法和遗传算法的连接,第二阶段和第三阶段的分界点划分遗传算法的精度判断,一旦达到该分界点,则本次优化结束,为避免陷入死循环,最大优化次数为变电站的节点数,根据正态分布函数的特征,σ在第1阶段时ρ值变化非常缓慢,在算法执行初期尽量保证信息的完整,体现了慢启动的特点,而在σ这个拐点后急剧下降,在后面的优化中速度加快,符合算法的设计思想;步骤7:将当前正在优化的变电站从节点集中去掉,返回步骤4进行下一个优化节点的判断,直至所有变电站的优化顺序都排列。 |