发明名称 一种粒子群寻优的调峰炉热力站供热负荷优化调度方法
摘要 一种粒子群寻优的调峰炉热力站供热负荷优化调度方法,属于热力站供热负荷控制技术领域。本发明为了使热力站的优化调度能兼顾供热系统的能耗和经济两个方面,通过在供热高峰期和非高峰期不同调度模式,使其能达到综合最优,使优化调度的计算的速度快,结果更准确。主要步骤:构建调峰炉热力站优化调度模型<img file="DDA00002870492600011.GIF" wi="501" he="60" />确定权系数λ<sub>1</sub>,λ<sub>2</sub>;采用粒子群寻优算法求得满足约束条件下综合能耗E<sub>w</sub>最小时对应的锅炉供热量Q<sub>b</sub>和二级网锅炉供热通道流量G<sub>b</sub>,至此完成调峰炉热力站供热负荷优化调度。本方法通过建立能耗性和经济性结合调峰炉热力站优化调度模型,实现了对不同热源供热负荷的最优分配,使能耗和经济的综合最优,达到了节能的目的。
申请公布号 CN103197547A 申请公布日期 2013.07.10
申请号 CN201310063770.6 申请日期 2013.02.28
申请人 哈尔滨工业大学 发明人 齐维贵;张永明;邓盛川;于德亮
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 杨立超
主权项 1.一种粒子群寻优的调峰炉热力站供热负荷优化调度方法,其特征在于,所述方法的实现过程为:步骤一、构建调峰炉热力站优化调度模型:对能耗支出和运行费用分别加以权重,通过线性加权和法构造的评价函数为:<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>E</mi><mi>w</mi><mo>*</mo></msubsup><mo>=</mo><mi>min</mi><mrow><mo>(</mo><msub><mi>&lambda;</mi><mn>1</mn></msub><msub><mi>E</mi><mi>E</mi></msub><mo>+</mo><msub><mi>&lambda;</mi><mn>2</mn></msub><msub><mi>E</mi><mi>C</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,E<sub>w</sub>为优化调度综合能耗,λ<sub>1</sub>,λ<sub>2</sub>分别为能耗支出权重系数和运行费用权重系数,λ<sub>1</sub>、λ<sub>2</sub>由通过熵权系数的计算获得,λ<sub>1</sub>≥0,λ<sub>2</sub>≥0;<img file="FDA00002870492300012.GIF" wi="58" he="60" />为加以权重计算后的综合能耗最小值;<img file="FDA00002870492300013.GIF" wi="60" he="80" />为能耗费用,E<sub>E</sub>表示能耗费用优化结果;<img file="FDA00002870492300014.GIF" wi="54" he="67" />为运行费用消耗,E<sub>C</sub>为运行费用消耗的优化结果;<maths num="0002"><![CDATA[<math><mrow><msubsup><mi>E</mi><mi>E</mi><mo>*</mo></msubsup><mo>=</mo><mi>min</mi><mrow><mo>(</mo><msub><mi>E</mi><mn>1</mn></msub><mo>+</mo><msub><mi>E</mi><mn>2</mn></msub><mo>+</mo><msub><mi>E</mi><mn>3</mn></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,E<sub>1</sub>为热网供热负荷燃料能耗;E<sub>2</sub>为热力站锅炉供热负荷燃料能耗;E<sub>3</sub>为热力站二级网循环水泵电耗;E<sub>1</sub>=(T<sub>g1</sub>-T<sub>h1</sub>)×G<sub>1</sub>×C×k<sub>1</sub>                   (3)式中,T<sub>g1</sub>、T<sub>h1</sub>为分别为热力站一级网供水温度和回水温度;G<sub>1</sub>为力站一级网流量;C为水的比热容;k<sub>1</sub>为热电厂标准煤耗率;E<sub>2</sub>=(T<sub>b2</sub>-T<sub>h2</sub>)×G<sub>b</sub>×C×1.43/k<sub>oil</sub>η<sub>b</sub>            (4)式中,T<sub>b2</sub>、T<sub>h2</sub>为分别为热力站燃油炉供水温度和二级网回水温度;G<sub>b</sub>为热力站二级网燃油锅炉供水流量;k<sub>oil</sub>为取燃油的低位发热值;η<sub>b</sub>为燃油锅炉总的燃烧热效率;<maths num="0003"><![CDATA[<math><mrow><msub><mi>E</mi><mn>3</mn></msub><mo>=</mo><msub><mi>k</mi><mi>e</mi></msub><mo>&times;</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><mn>2.78</mn><mo>&times;</mo><msup><mn>10</mn><mrow><mo>-</mo><mn>4</mn></mrow></msup><mfrac><mrow><msub><mi>G</mi><mi>p</mi></msub><mi>H</mi></mrow><msub><mi>&rho;&eta;</mi><mi>p</mi></msub></mfrac><msub><mi>n</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,G<sub>p</sub>为循环水泵流量;H为循环水泵扬程;n<sub>0</sub>为循环水泵运行时间;η<sub>p</sub>为水泵运行效率;ρ为水密度;m为燃油锅炉运行台数;k<sub>e</sub>为电折算为标准煤的系数;<maths num="0004"><![CDATA[<math><mrow><msubsup><mi>E</mi><mi>C</mi><mo>*</mo></msubsup><mo>=</mo><mi>min</mi><mrow><mo>(</mo><msub><mi>E</mi><mn>1</mn></msub><mo>+</mo><msub><mi>E</mi><mn>2</mn></msub><msub><mi>k</mi><mi>b</mi></msub><msub><mi>k</mi><mi>c</mi></msub><mo>/</mo><mn>1.43</mn><mo>+</mo><mfrac><mrow><msub><mi>E</mi><mn>3</mn></msub><msub><mi>k</mi><mi>E</mi></msub><msub><mi>k</mi><mi>c</mi></msub></mrow><msub><mi>k</mi><mi>e</mi></msub></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,k<sub>C</sub>为标准煤的价格;k<sub>b</sub>为燃油的价格;k<sub>E</sub>为电的价格;将式(2)和式(6)代入式(1),采用n台同型号的水泵,假设所述水泵具有相同的工况,有:E<sub>w</sub>=λ<sub>1</sub>E<sub>E</sub>+λ<sub>2</sub>E<sub>C</sub>                                       (7)=J<sub>1</sub>×(T<sub>g1</sub>-T<sub>h1</sub>)×G<sub>1</sub>+J<sub>2</sub>×(T<sub>b2</sub>-T<sub>h2</sub>)×G<sub>b</sub>/η<sub>b</sub>+J<sub>3</sub>nG<sub>p</sub>式中,J<sub>1</sub>=c×k<sub>1</sub>×(λ<sub>1</sub>+λ<sub>2</sub>),J<sub>2</sub>=c×(1.43λ<sub>1</sub>+λ<sub>2</sub>k<sub>b</sub>k<sub>C</sub>)/k<sub>oil</sub>,<maths num="0005"><![CDATA[<math><mrow><msub><mi>J</mi><mn>3</mn></msub><mo>=</mo><mn>2.78</mn><mo>&times;</mo><msup><mn>10</mn><mrow><mo>-</mo><mn>4</mn></mrow></msup><mfrac><mi>H</mi><msub><mi>&rho;&eta;</mi><mi>p</mi></msub></mfrac><msub><mi>n</mi><mn>0</mn></msub><mo>&times;</mo><mrow><mo>(</mo><msub><mi>&lambda;</mi><mn>1</mn></msub><msub><mi>k</mi><mi>e</mi></msub><mo>+</mo><msub><mi>&lambda;</mi><mn>2</mn></msub><msub><mi>k</mi><mi>E</mi></msub><msub><mi>k</mi><mi>C</mi></msub><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>由此,可建立调峰炉热力站调度模型如下:minE<sub>w</sub>=minJ<sub>1</sub>×(T<sub>g1</sub>-T<sub>h1</sub>)×G<sub>1</sub>+J<sub>2</sub>×(T<sub>b2</sub>-T<sub>h2</sub>)×G<sub>b</sub>/η<sub>b</sub>+J<sub>3</sub>nG<sub>p</sub><maths num="0006"><![CDATA[<math><mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>G</mi><mn>1</mn></msub><mo>&le;</mo><msub><mi>G</mi><mrow><mn>1</mn><mi>max</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>G</mi><mi>H</mi></msub><mo>&le;</mo><msub><mi>G</mi><mrow><mi>z</mi><mn>2</mn></mrow></msub><mo>&le;</mo><msub><mi>G</mi><mrow><mi>p</mi><mi>max</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>G</mi><mrow><mi>z</mi><mn>2</mn></mrow></msub><mo>=</mo><msub><mi>nG</mi><mi>p</mi></msub><mo>=</mo><msub><mi>G</mi><mn>2</mn></msub><mo>+</mo><msub><mi>G</mi><mi>b</mi></msub></mtd></mtr><mtr><mtd><msub><mi>K</mi><mi>min</mi></msub><mo>&le;</mo><msub><mi>K</mi><mi>b</mi></msub></mtd></mtr><mtr><mtd><msub><mi>Q</mi><mi>ex</mi></msub><mo>=</mo><msub><mi>K</mi><mi>b</mi></msub><msub><mi>A&Delta;t</mi><mi>m</mi></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>T</mi><mrow><mi>g</mi><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>T</mi><mrow><mi>h</mi><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>&times;</mo><mi>c</mi><mo>&times;</mo><msub><mi>G</mi><mn>1</mn></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>T</mi><mrow><mi>g</mi><mn>2</mn></mrow></msub><mo>-</mo><msub><mi>T</mi><mrow><mi>h</mi><mn>2</mn></mrow></msub><mo>)</mo></mrow><mo>&times;</mo><mi>c</mi><mo>&times;</mo><msub><mi>G</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><msub><mi>Q</mi><mi>out</mi></msub><mo>=</mo><msub><mi>Q</mi><mi>ex</mi></msub><mo>+</mo><msub><mi>Q</mi><mi>b</mi></msub></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>式中J<sub>1</sub>,J<sub>2</sub>,J<sub>3</sub>——各项常系数;E<sub>w</sub>——优化调度综合能耗(kg标煤);T<sub>g1</sub>,T<sub>h1</sub>——分别为热力站一级网供水温度和回水温度(℃);T<sub>b2</sub>,T<sub>h2</sub>——分别为热力站燃油炉供水温度和二级网回水温度(℃);G<sub>1</sub>——热力站一级网流量(m<sup>3</sup>/h);G<sub>1max</sub>——热力站一级网流量最大值(m<sup>3</sup>/h);G<sub>b</sub>——热力站二级网燃油锅炉供水流量(m<sup>3</sup>/h);G<sub>z2</sub>——热力站二级网总输出供水流量(m<sup>3</sup>/h);G<sub>pmax</sub>——热力站二级网循环泵最大可输出流量(m<sup>3</sup>/h);G<sub>H</sub>——热力站二级网达到最小的扬程时循环泵最小输出流量(m<sup>3</sup>/h);η<sub>b</sub>——燃油锅炉总的燃烧热效率;n——循环水泵运行台数(台);G<sub>p</sub>——循环水泵流量(m<sup>3</sup>/h);Q<sub>out</sub>——供热负荷最终输出值(GJ);Q<sub>ex</sub>——换热器供出热量(GJ);Q<sub>b</sub>——锅炉供出热量(GJ);K<sub>b</sub>——板式换热器换热系数;c——水的比热容为4.187×10<sup>-3</sup>(GJ/kg·℃);步骤二、确定权系数λ<sub>1</sub>,λ<sub>2</sub>当有m个评价对象,每个对象有n个评价指标时,可构成一个评价指标矩阵R=(r<sub>ij</sub>)<sub>m×n</sub>;第i对象下的第j个评价指标的比重P<sub>ij</sub>为:<maths num="0007"><![CDATA[<math><mrow><msub><mi>P</mi><mi>ij</mi></msub><mo>=</mo><msub><mi>r</mi><mi>ij</mi></msub><mo>/</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>r</mi><mi>ij</mi></msub><mo>,</mo><mrow><mo>(</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths>则可得出第j个评价指标的熵值为:<maths num="0008"><![CDATA[<math><mrow><msub><mi>E</mi><mi>j</mi></msub><mo>=</mo><mo>-</mo><mi>ln</mi><msup><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>P</mi><mi>ij</mi></msub><mi>ln</mi><mrow><mo>(</mo><msub><mi>P</mi><mi>ij</mi></msub><mo>)</mo></mrow><mo>,</mo><mrow><mo>(</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow></math>]]></maths>第j个评价指标的权重ω<sub>j</sub>可表示为:<maths num="0009"><![CDATA[<math><mrow><msub><mi>&omega;</mi><mi>j</mi></msub><mo>=</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>E</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>/</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>E</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>,</mo><mrow><mo>(</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,0≤ω<sub>j</sub>≤1,且有<img file="FDA00002870492300034.GIF" wi="202" he="122" />在获得客观评价的权重ω<sub>j</sub>后,加入专家的主观评价权重α<sub>j</sub>,使其能全面地反映出评价指标的重要性,得到一个复合权重为:<maths num="0010"><![CDATA[<math><mrow><msub><mi>&beta;</mi><mi>j</mi></msub><mo>=</mo><msub><mi>&omega;</mi><mi>j</mi></msub><msub><mi>&alpha;</mi><mi>j</mi></msub><mo>/</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>&omega;</mi><mi>j</mi></msub><msub><mi>&alpha;</mi><mi>j</mi></msub><mo>,</mo><mrow><mo>(</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow></math>]]></maths>记R矩阵每列的最优值为<img file="FDA00002870492300036.GIF" wi="61" he="66" />即:<img file="FDA00002870492300037.GIF" wi="1268" he="181" />记r<sub>ij</sub>与<img file="FDA00002870492300038.GIF" wi="40" he="66" />接近程度为:<img file="FDA00002870492300039.GIF" wi="1234" he="313" />则获得综合评价的权重系数λ<sub>i</sub>的表示为:<maths num="0011"><![CDATA[<math><mrow><msub><mi>&lambda;</mi><mi>i</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>&beta;</mi><mi>j</mi></msub><msub><mi>d</mi><mi>ij</mi></msub><mo>,</mo><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>14</mn><mo>)</mo></mrow></mrow></math>]]></maths>综合评价的权重系数λ<sub>i</sub>可全面反映出能耗支出最小运行模式和运行费用最少运行模式的优劣和对综合优化调度模式的影响;步骤三、采用改进的粒子群寻优算法对调峰炉热力站的供热负荷进行优化调度,粒子的适应值即最小的综合能耗E<sub>w</sub>,需求解的粒子包括有锅炉供热量Q<sub>b</sub>和二级网锅炉供热通道流量G<sub>b</sub>;即采用粒子群寻优算法求得满足约束条件下综合能耗E<sub>w</sub>最小时对应的锅炉供热量Q<sub>b</sub>和二级网锅炉供热通道流量G<sub>b</sub>,至此完成调峰炉热力站供热负荷优化调度。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号