发明名称 一种基于改进神经网络与遗传算法的燃煤锅炉优化方法
摘要 本发明涉及人工智能应用领域,具体是一种基于改进神经网络与遗传算法的燃煤锅炉运行优化预测的方法。主要原理是:通过BP神经网络的优化、预测模型的建立及训练,从而实现燃煤锅炉运行优化。具体步骤是:应用基于遗传算法优化的神经网络对锅炉参数(燃料量、送风量、氧量、温差量、煤种特性)进行实验数据建模,进而利用遗传算法对某一工况进行锅炉可调参数(燃料量、送风量、氧量、温差量)的优化。该方法可以利用优化后的神经网络和遗传算法对电站锅炉进行优化。
申请公布号 CN103324862B 申请公布日期 2016.02.24
申请号 CN201310290049.0 申请日期 2013.07.11
申请人 中国石油大学(华东) 发明人 李克文;高俊苓
分类号 G06F17/50(2006.01)I;G06N3/02(2006.01)I;G06N3/12(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 代理人
主权项 一种基于改进BP神经网络与遗传算法的锅炉生产运行优化方法,其特征在于,包括以下步骤:A.BP神经网络的优化采用遗传算法对神经网络权值进行优化:首先将神经网络的权值进行编码,然后利用遗传算法的选择、交叉、变异操作,得到新的权值,比较训练目标值是否达到0.05,若未达到,则继续利用遗传算法进行神经网络权值寻优;其中:1)适应度函数为F=1/Σ(t<sub>i</sub>‑x<sub>i</sub>)<sup>2</sup>,其中t为实际输出值,x为期望值;2)第i个个体选择概率ps(i)为ps(i)=r×(1‑q)<sup>M‑1</sup>,其中<img file="FDA0000888795500000011.GIF" wi="349" he="127" />M为种群数,q为选择系数;3)交叉操作:利用实数型的算术交叉,根据选择概率选择一对父代个体y(j,i<sub>1</sub>)和y(j,i<sub>2</sub>),进行如下交叉运算产生两个子代个体y(j+1,i<sub>1</sub>)和y(j+1,i<sub>2</sub>):<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mi>y</mi><mo>(</mo><mi>j</mi><mo>+</mo><mn>1</mn><mo>,</mo><msub><mi>i</mi><mn>1</mn></msub><mo>)</mo><mo>=</mo><mi>&alpha;</mi><mi>y</mi><mo>(</mo><mi>j</mi><mo>,</mo><msub><mi>i</mi><mn>2</mn></msub><mo>)</mo><mo>+</mo><mo>(</mo><mn>1</mn><mo>-</mo><mi>&alpha;</mi><mo>)</mo><mi>y</mi><mo>(</mo><mi>j</mi><mo>,</mo><msub><mi>i</mi><mn>1</mn></msub><mo>)</mo></mtd></mtr><mtr><mtd><mi>y</mi><mo>(</mo><mi>j</mi><mo>+</mo><mn>1</mn><mo>,</mo><msub><mi>i</mi><mn>2</mn></msub><mo>)</mo><mo>=</mo><mi>&alpha;</mi><mi>y</mi><mo>(</mo><mi>j</mi><mo>,</mo><msub><mi>i</mi><mn>1</mn></msub><mo>)</mo><mo>+</mo><mo>(</mo><mn>1</mn><mo>-</mo><mi>&alpha;</mi><mo>)</mo><mi>y</mi><mo>(</mo><mi>j</mi><mo>,</mo><msub><mi>i</mi><mn>2</mn></msub><mo>)</mo></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000888795500000012.GIF" wi="868" he="172" /></maths>其中,α是在区间(0,1)之间的随机数;4)变异操作:产生的下一代变异个体为<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><msubsup><mi>x</mi><mi>k</mi><mi>j</mi></msubsup><mo>+</mo><mo>(</mo><msubsup><mi>U</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mi>k</mi></msubsup><mo>-</mo><msubsup><mi>x</mi><mi>k</mi><mi>j</mi></msubsup><mo>)</mo><mo>&times;</mo><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>c</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mfrac><mi>j</mi><mi>T</mi></mfrac><mo>)</mo><mo>&times;</mo><mi>b</mi></mrow></msup><mo>)</mo><mo>,</mo><mi>R</mi><mi>a</mi><mi>n</mi><mi>d</mi><mi>o</mi><mi>m</mi><mo>&lsqb;</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>&rsqb;</mo><mo>=</mo><mn>1</mn></mtd></mtr><mtr><mtd><msubsup><mi>x</mi><mi>k</mi><mi>j</mi></msubsup><mo>+</mo><mo>(</mo><msubsup><mi>x</mi><mi>k</mi><mi>j</mi></msubsup><mo>-</mo><msubsup><mi>U</mi><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mi>k</mi></msubsup><mo>)</mo><mo>&times;</mo><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>c</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mfrac><mi>j</mi><mi>T</mi></mfrac><mo>)</mo><mo>&times;</mo><mi>b</mi></mrow></msup><mo>)</mo><mo>,</mo><mi>R</mi><mi>a</mi><mi>n</mi><mi>d</mi><mi>o</mi><mi>m</mi><mo>&lsqb;</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>&rsqb;</mo><mo>=</mo><mn>0</mn></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000888795500000013.GIF" wi="974" he="222" /></maths>其中,第j代父代个体<img file="FDA0000888795500000014.GIF" wi="487" he="102" />中的<img file="FDA0000888795500000015.GIF" wi="57" he="69" />为变异点,取值范围<img file="FDA0000888795500000016.GIF" wi="299" he="101" />c是区间[0,1]上符合均匀概率分布的随机数,T为终止代数,b为调整步长参数;B.预测模型的建立及训练采用BP神经网络建立锅炉高效率低污染模型和生产运行模型,并将其应用到锅炉生产运行预测与优化系统中:用燃料量、送风量、氧量、温差量、煤种特性五个参数作为输入变量;用锅炉效率、NO<sub>x</sub>排放两个参数作为输出变量,来训练神经网络模型;C.燃煤锅炉运行优化采用遗传算法作为优化算法,根据燃烧优化目标和燃烧过程神经网络模型将基于神经网络的燃煤锅炉预测模型的输出构造为目标函数进行工况参数寻优,应用遗传算法快速、准确地计算得到当前工况下燃料量、送风量、氧量、温差量四个操作量的最优设定值;其中:1)目标函数:目标函数定义为minf=a(η<sub>c</sub>/η<sub>fc</sub>)+(1‑a)([NOx]<sub>fc</sub>/[NOx]<sub>c</sub>),式中,η<sub>c</sub>,η<sub>fc</sub>为当前炉效及优化后炉效;[NOx]<sub>c</sub>,[NOx]<sub>fc</sub>为NO<sub>x</sub>排放物当前值及优化后的预测值;a、1‑a为锅炉效率和NO<sub>x</sub>浓度的加权;2)适应度函数:第i个个体的适应度函数值Fit(i)=exp(‑f(i)),其中f(i)目标函数;3)选择操作:第i个个体选择概率ps(i)为ps(i)=r×(1‑q)<sup>M‑1</sup>,其中<img file="FDA0000888795500000021.GIF" wi="349" he="126" />M为种群数,q为选择系数;4)交叉操作:利用实数型的算术交叉,根据选择概率选择一对父代个体y(j,i<sub>1</sub>)和y(j,i<sub>2</sub>),进行如下交叉运算产生两个子代个体y(j+1,i<sub>1</sub>)和y(j+1,i<sub>2</sub>):<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mi>y</mi><mo>(</mo><mi>j</mi><mo>+</mo><mn>1</mn><mo>,</mo><msub><mi>i</mi><mn>1</mn></msub><mo>)</mo><mo>=</mo><mi>&alpha;</mi><mi>y</mi><mo>(</mo><mi>j</mi><mo>,</mo><msub><mi>i</mi><mn>2</mn></msub><mo>)</mo><mo>+</mo><mo>(</mo><mn>1</mn><mo>-</mo><mi>&alpha;</mi><mo>)</mo><mi>y</mi><mo>(</mo><mi>j</mi><mo>,</mo><msub><mi>i</mi><mn>1</mn></msub><mo>)</mo></mtd></mtr><mtr><mtd><mi>y</mi><mo>(</mo><mi>j</mi><mo>+</mo><mn>1</mn><mo>,</mo><msub><mi>i</mi><mn>2</mn></msub><mo>)</mo><mo>=</mo><mi>&alpha;</mi><mi>y</mi><mo>(</mo><mi>j</mi><mo>,</mo><msub><mi>i</mi><mn>1</mn></msub><mo>)</mo><mo>+</mo><mo>(</mo><mn>1</mn><mo>-</mo><mi>&alpha;</mi><mo>)</mo><mi>y</mi><mo>(</mo><mi>j</mi><mo>,</mo><msub><mi>i</mi><mn>2</mn></msub><mo>)</mo></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000888795500000022.GIF" wi="889" he="174" /></maths>式中,α是在区间(0,1)之间的随机数;5)变异操作,产生的下一代变异个体为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><msubsup><mi>x</mi><mi>k</mi><mi>j</mi></msubsup><mo>+</mo><mo>(</mo><msubsup><mi>U</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mi>k</mi></msubsup><mo>-</mo><msubsup><mi>x</mi><mi>k</mi><mi>j</mi></msubsup><mo>)</mo><mo>&times;</mo><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>c</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mfrac><mi>j</mi><mi>T</mi></mfrac><mo>)</mo><mo>&times;</mo><mi>b</mi></mrow></msup><mo>)</mo><mo>,</mo><mi>R</mi><mi>a</mi><mi>n</mi><mi>d</mi><mi>o</mi><mi>m</mi><mo>&lsqb;</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>&rsqb;</mo><mo>=</mo><mn>1</mn></mtd></mtr><mtr><mtd><msubsup><mi>x</mi><mi>k</mi><mi>j</mi></msubsup><mo>+</mo><mo>(</mo><msubsup><mi>x</mi><mi>k</mi><mi>j</mi></msubsup><mo>-</mo><msubsup><mi>U</mi><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mi>k</mi></msubsup><mo>)</mo><mo>&times;</mo><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>c</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mfrac><mi>j</mi><mi>T</mi></mfrac><mo>)</mo><mo>&times;</mo><mi>b</mi></mrow></msup><mo>)</mo><mo>,</mo><mi>R</mi><mi>a</mi><mi>n</mi><mi>d</mi><mi>o</mi><mi>m</mi><mo>&lsqb;</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>&rsqb;</mo><mo>=</mo><mn>0</mn></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000888795500000023.GIF" wi="1004" he="231" /></maths>其中,第j代父代个体<img file="FDA0000888795500000024.GIF" wi="485" he="102" />中的<img file="FDA0000888795500000025.GIF" wi="55" he="76" />为变异点,取值范围<img file="FDA0000888795500000026.GIF" wi="302" he="111" />c是区间[0,1]上符合均匀概率分布的随机数,T为终止代数,b为调整步长参数。
地址 266580 山东省青岛市青岛经济技术开发区长江西路66号