发明名称 一种基于改进遗传算法的超超临界机组主汽温PID控制方法
摘要 本发明提供了一种基于改进遗传算法的火力发电厂锅炉发电超超临界机组主汽温PID控制方法,方法中包括有改进遗传算法。本发明的改进遗传算法对简单遗传算法进行了改进,引入排序保优选择、非均匀线性交叉和高斯变异算子,采用均匀设计方法构造初始种群,通过Sigmoid函数适应性调整交叉和变异概率;采用面向DCS系统的PID控制器,该控制器在分负荷段下的参数集通过改进的遗传算法寻优得到。本发明提高了遗传算法的寻优速度和精度,可以代替工程复杂且准确性不高的工程整定方法,降低了现场调试人员的工作量,达到主汽温智能控制的目的,实现了智能优化技术在DCS系统下的应用,满足超超临界机组主汽温在深度调峰时的稳定控制需求。
申请公布号 CN102777878B 申请公布日期 2015.02.11
申请号 CN201210234722.4 申请日期 2012.07.06
申请人 广东电网公司电力科学研究院;武汉大学 发明人 陈世和;易凤飞;韩玲;方彦军
分类号 F22B35/00(2006.01)I;G05B11/42(2006.01)I 主分类号 F22B35/00(2006.01)I
代理机构 广州知友专利商标代理有限公司 44104 代理人 周克佑
主权项 一种基于改进遗传算法的超超临界机组主汽温PID控制方法,其特征是:所述的方法中包括改进遗传算法,所述的改进遗传算法包括以下步骤:Step1: 建立PID参数整定问题数学模型,包括目标函数和约束条件;Step2: 设置遗传参数,包括种群大小、进化代数、选择压力、交叉概率范围、变异概率范围;Step3: 采用实数编码方式编码PID参数<img file="2012102347224100001dest_path_image001.GIF" wi="8" he="15" />,并设置参数上、下限;Step4: 建立均匀设计表,并根据均匀表对种群初始化,产生初始种群<img file="588827dest_path_image002.GIF" wi="23" he="14" />;Step5: 判断适应度是否收敛,收敛则终止程序,不收敛则进入下一步骤;Step6: 计算每一代种群中个体的适应度<img file="941311dest_path_image003.GIF" wi="30" he="18" />;Step7: 采用线序排序算子和精英选择策略选择优秀个体;Step8: 采用Sigmoid函数公式计算交叉概率<img file="821542dest_path_image004.GIF" wi="15" he="42" />:概率调整曲线采用Sigmoid函数,其计算公式为<img file="259476dest_path_image005.GIF" wi="138" he="39" />式中,<i>P</i><sub>c1</sub>为初始交叉概率,<i>P</i><sub>c2</sub>为终止交叉概率,<img file="617777dest_path_image006.GIF" wi="9" he="18" />为形状因子,设定为20;<i>Ns</i>为分界点,设定为0.25;Step9: 采用非均匀线性交叉算子产生新的基因;Step10:采用Sigmoid函数公式计算变异概率<img file="2012102347224100001dest_path_image007.GIF" wi="20" he="42" />:概率调整曲线采用Sigmoid函数,其计算公式为<img file="28029dest_path_image008.GIF" wi="149" he="38" />式中,<i>P</i><sub>m1</sub>为初始变异概率;<i>P</i><sub>m2</sub>为终止变异概率;<img file="141479dest_path_image006.GIF" wi="9" he="18" />为形状因子,设定为20;<i>Ns</i>为分界点,设定为0.25;Step11: 采用高斯变异算子产生新的基因;Step12: 记录最个体、最佳适应度和子代种群信息;Step13: 判断是否达到最大进化代数,若否则返回Step6,重复以上运算过程直到达到最大进化代数;Step14: 获得最优解,终止遗传算法优化程序;所述的步骤Step4包括以下子步骤:Step4‑1:给定均匀设计表的水平数<img file="2012102347224100001dest_path_image009.GIF" wi="9" he="12" />,当<img file="831541dest_path_image009.GIF" wi="9" he="12" />为奇数时,生成一组由正整数向量<img file="993532dest_path_image010.GIF" wi="144" he="17" />,其中<img file="55029dest_path_image011.GIF" wi="9" he="15" />比<img file="277063dest_path_image012.GIF" wi="9" he="12" />小,且它们的最大公因素为1;当<img file="689590dest_path_image009.GIF" wi="9" he="12" />为偶数时,先构造均匀设计表<img file="655272dest_path_image013.GIF" wi="64" he="19" />和向量<img file="836854dest_path_image014.GIF" wi="33" he="17" />;Step4‑2:构造均匀设计表中第<img file="964210dest_path_image015.GIF" wi="9" he="15" />列,其个体计算公式为<img file="864033dest_path_image016.GIF" wi="129" he="14" />,其中<img file="2012102347224100001dest_path_image017.GIF" wi="30" he="15" />表示同余运算,<img file="367827dest_path_image018.GIF" wi="92" he="14" />,<img file="2012102347224100001dest_path_image019.GIF" wi="96" he="14" />,同余运算可以保证个体落在区间<img file="607178dest_path_image020.GIF" wi="30" he="18" />,从而得到了均匀设计表<img file="905436dest_path_image021.GIF" wi="36" he="17" />的<img file="292555dest_path_image022.GIF" wi="33" he="14" />矩阵,当<img file="597110dest_path_image009.GIF" wi="9" he="12" />为偶数时,产生矩阵<img file="2012102347224100001dest_path_image023.GIF" wi="56" he="17" />;Step4‑3:采用L<sub>2</sub>‑偏差函数<img file="425388dest_path_image024.GIF" wi="54" he="17" />来测试其均匀性,选出均匀性做最好的<img file="222443dest_path_image025.GIF" wi="9" he="12" />列构成均匀设计表<img file="831279dest_path_image026.GIF" wi="33" he="17" />,其函数方程为<img file="942454dest_path_image027.GIF" wi="285" he="120" />Step4‑4:通过均匀设计表产生遗传算法初始种群,其中<img file="953136dest_path_image025.GIF" wi="9" he="12" />为参数个数,需要根据参数搜索问题进行设定,种群大小为<img file="858775dest_path_image009.GIF" wi="9" he="12" />,可构造一个<img file="954907dest_path_image025.GIF" wi="9" he="12" />因素<img file="604194dest_path_image009.GIF" wi="9" he="12" />水平的均匀设计表,使初始种群中的个体均匀分布在<img file="469382dest_path_image025.GIF" wi="9" he="12" />维空间。
地址 510080 广东省广州市越秀区东风东路水均岗8号