专利名称：混烧低热值煤气机组煤耗计算方法
技术领域：
本发明涉及一种混烧低热值煤气机组煤耗计算方法，特别是一种大型电站机组混 烧低热值高炉煤气机组的煤耗计算方法，属于电站节能降耗技术领域。
背景技术：
近年来，煤炭价格走高以及节能减排的要求，促使电厂进行改革以减少能耗，提高 经济性，因此对于电厂性能分析提出了很高的要求。发电厂常规的机组性能分析是离线状 态进行的，利用数据采集系统(DCS)采集数据，基于机组热量平衡和质量平衡的计算方法， 主要步骤是焓熵图表程序化、各级抽汽量的计算、热耗计算和能损分析，这样的分析需要的 计算参数众多，尤其是对于混烧低热值煤气的气煤混烧的机组，计算煤耗就更加复杂。当这 些测点的传感器出现故障时，其计算结果将出现较大偏差，一般来说会无法进行。因此采用 热平衡法等常规计算方法，存在原始数据数量大、计算复杂且维护不当时性能计算误差大 或无法进行等问题，不太适合在线的节能减排管理。用现有的技术方法计算一次机组煤耗需要54个数据的测点(见表1)，如果其中有 一个测点的传感器损坏，将无法进行计算，如果对所有测点的传感器设置硬件冗余，造成硬 件费用增加，成本提高。如果硬件冗余不足，一旦有测点的传感器损坏，性能分析系统就不 能正常运行，影响机组节能减排的评估。表1常规方法煤耗计算所需数据的测点
发明内容
本发明针对上述现有技术的不足，提出了一种混烧低热值煤气机组煤耗计算方 法，即在由现场数据传感器组采集的数据中选取少量测点的数据，然后利用本方法计算即可得到机组煤耗。本发明不需要复杂的焓熵图计算程序，也不需要对所有数据的测点传感 器设置硬件冗余，只是利用人工神经网络建立起输入和输出的非线性映射关系，就可以在 线计算出精度较高的发电煤耗。本发明是通过如下技术方案实现的，本发明方法包括步骤一、原始数据准备。对电厂混烧低热值煤气机组进行热力试验，获得各种发电功率下的54个测点数 据和用传统的热平衡计算方法计算得到的发电煤耗数据；步骤二、人工控制下的学习建模(1)原始数据选取选用步骤一中计算得到的发电煤耗数据作为目标输出煤耗数据，然后根据经验和 敏感性分析得到在54个测点中选择15个测点(见表2)的数据作为原始数据，选取相同时 刻的一组目标输出煤耗数据和原始数据作为一个样本，选取能够覆盖所有运行功率的N组 (N ^ 250)样本构成一个原始数据集。表2原始输入数据测点清单 针对具有N组样本的原始数据集中的每一个测点的数据列，找出其最大值和最小 值。对于15个测点的原始数据和一个目标输出煤耗数据，就有15+1个最大值和15+1个最小值。(2)原始数据集归一化由于网络计算数据范围为
，因此，在进行训练之前必须将原始数据集进行归 一化，以减小计算量和避免误差的扩大。归一化采用公式(1)计算 Xjp = 其中
min X，ip表示原始数据集中第P组样本第i个测点的原始数据，xiD为X’ iD归一化后
的值，称为原始输入数据；i = 1，2，. . .，15，表示第i个测点的原始数据；i = 0时，Vtlp表示原始数据集中第ρ组样本目标输出煤耗数据，目标输出数据Xtlp 为Vtlp归一化后的值；ρ= 1,2,..., N，表示第ρ组原始数据集的样本，N为原始数据集的样本总数；χ'大值。
分别表示原始数据集中第i个测点数据列的N个数据中的最小和最
χ，0fflin, χ，0fflax表示原始数据集中的目标输出煤耗数据的最小值和最大值。(3)初始化权矩阵Vij和Wj在训练学习开始，需要对权矩阵Vu和Wi进行赋初值，这里Vu是mX 15阶矩阵,Wj 是IXm阶矩阵，计算之前先对Vu和^矩阵内每个元素赋初始值为0 1之间的随机数。 其中m为该模型的中间隐层节点数；(4)针对ρ = 1，2，...，N，学习训练，获得权矩阵Vij和Wj ①计算中间隐层节点yw.；对于第ρ组学习样本，有 其中 ②计算中间目标数据dp: ( 其中 ③计算误差对每一个中间目标数据和目标输出数据进行误差计算，为 ④调整修正计算获得新的连接权矩阵Vij和Wj，调整方法计算中间变量δ。S0 =^(dp-x0p)(l-xlp)2dy =\{·1 — y2pj)Sowi j = 1，2，· · ·，m令p= p+l，计算Affjp= η 5oy(p_1)J+a Δ Wj^)Δ Vjip = η δ yXi (ρ_1} +α Δ Vji (^1) j = 1，2，· · ·，m(5)Wj' = Wj+ Δ Wjpi = 1· 2·…，15Vji' =Yji+AYjip为了进行下一轮的迭代计算，需要进行数据准备令Vij = Vi/Wj = W/;式中AWjp、八、为计算第ρ组到第p+1组学习样本时调整权矩阵的增量， AWj(p_D、AVjip^为前一组调整的权矩阵的增量，式(5)中的Vij^W/表示第P+1组学习样 本的权矩阵，y(p_lh.、Xifp-D为前一组样本的中间隐层节点和第i个输入数据的原始输入参 数；α、η是为了加快收敛速度而设置的两个系数，一般取0 1之间的任意数，开始运算
7时取0. 5，之后根据收敛状态可以进行增加5%或降低5%的调整。⑤判断N组样本全计算结束了吗？如果没有计算完，则重复进行① ④的计算， 直到N组样本全计算结束，得到中间目标数据dp，然后计算总输出误差E 如果E不满足精度要求，则重新针对ρ = 1,2,...，N，再进行(2) (6)式的计算， 这样重复进行，直到误差达到要求。通过调整各Vu、Wj,可以求得输出的误差平方和E达到要求精度时的V『Wj，将该Vu、 Wpm以及χ'imax作为混烧低热值煤气机组煤耗计算方法模型保存在数据库模块内。步骤三、在线混烧低热值煤气机组煤耗计算方法实现通过上述人工控制下的在线学习建模过程，得到保存在数据库模块中混烧低热值 煤气机组煤耗计算方法模型。实际使用时，将现场实测的如表2所示的15个测点的数据作 为混烧低热值煤气机组煤耗计算方法模型的原始数据，即X’ χ’ 2、…、χ’ 15，采用数据库 模块中保存的χ'x' imax经过⑴式归一化后得到原始输入数据Χι、χ2>…、X15，将 数据库模块中保存的Vu、Wpm根据公式(7)和⑶计算，即可得到所需要的经过归一化的 目标输出数据d，采用数据库模块中保存的χ'-,经过反归一化式(9)后就可得 到目标输出煤耗数据χ’ ο，用来作为机组经济性指标输出。
χ' ο = d · (x' 0max-x' 0min)+x' 0min (9)有益效果发电厂常规的机组性能分析是离线状态进行的，基于机组热量平衡和质量平衡的 计算方法，当这些测点的传感器出现故障时，其计算结果将出现较大偏差，一般来说会无法 进行。因此采用热平衡法等常规计算方法，存在原始数据数量大、计算复杂且维护不当时性 能计算误差大或无法进行等问题，不太适合。本方法基于现场实测数据和人工神经网络的 有机结合，利用数据采集系统(DCS)采集数据，可以在不需要焓熵图表程序和只需要15个 测点参数下计算煤耗，这对于混烧低热值煤气的气煤混烧的机组实现在线系统煤耗计算和 在线的节能减排管理成为可能，提高了在线系统的可靠性。由于采用了大量现场实际信息， 故所得的煤耗精度一般在0. 1 %以内，符合工程应用的要求。
图1是本发明混烧低热值煤气机组煤耗计算方法的系统原理图
8
图2是本发明混烧低热值煤气机组煤耗计算方法的在线学习建模示意3是本发明混烧低热值煤气机组煤耗计算方法示意图
具体实施例方式下面结合附图对本发明的具体实施作进一步的描述。如图1、图2和图3所示，本发明在现场电站机组及测量系统和电厂经济性管理系 统之间嵌入一个混烧低热值煤气机组煤耗计算系统。该系统由数据库模块、数据选取模块、 人工控制下的在线学习建模模块和在线混烧低热值煤气机组煤耗计算模块组成。本发明实施方法包括步骤一、原始数据准备对电厂混烧低热值煤气机组进行热力试验，获得各种发电功率下的54个测点数 据和用传统的热平衡计算方法计算得到的发电煤耗数据；步骤二、人工控制下的在线学习建模(1)原始数据选取选用步骤一中计算得到的发电煤耗数据作为目标输出煤耗数据，然后根据表2中 15个测点的数据作为原始数据，选取相同时刻的一组目标输出煤耗数据和原始数据作为一 个样本，共取250组样本数据构成一个原始数据集，如表3所示。表3 将表3中的每一列数据下面的一列共有250个数据，选择其中的最大值记为 X’ imax，将其中的最小值记为 X’ imin，得到表 4 的 X’ imax、X’ imin，(i = 0，1，. . .，15)。表 4 (2)原始数据集归一化
将表3数据运用式(1)和表4的最大值和最小值进行归一化计算，其结果如表5 所示。表 5 (3)建立初始权矩阵Vij和Wj在训练学习开始，根据经验选m = 12，建立初始权矩阵Vij和Wj,这里Vij是12 X 15 阶矩阵，W」是1 X 12阶矩阵，对初始值Vu和^矩阵内每个元素赋初始值为0 1之间的随 机数。(4)针对ρ = 1，2，...，250，学习训练，获得模型Vij和Wj ①利用原始数据集中的原始输入数据，对ρ = 1,2,...，250组学习样本，根据公式 (2)计算中间隐层节点；②根据公式(3)计算未经误差分析的中间目标数据dp ；③根据公式(4)计算中间目标数据与目标输出数据的误差；④根据公式(5)调整修正计算获得新的权矩阵Vij和Wj ；⑤判断250组样本全计算结束了吗？如果没有计算完，则重复进行① ④计算， 直到250组样本全计算结束，得到中间目标数据dp，并根据公式(6)计算得到总输出误差如果总输出误差不满足精度要求(le-3)，则重新对250个样本进行(2) (6)式 的计算，这样重复进行，直到误差达到要求。将通过以上计算误差达到要求时的V^Wpm以及χ'imax(见表4)作为 混烧低热值煤气机组煤耗计算模型保存在数据库模块内供在线混烧低热值煤气机组煤耗 计算时使用。步骤三、在线混烧低热值煤气机组煤耗计算过程实际使用时，将现场实测的如表2所示的15个测点的数据作为混烧低热值煤气机 组煤耗计算方法模型的原始数据，即X’ X’ 2、…、X’ 15，采用数据库模块中保存的X' imin 和X' imax经过(1)式归一化后得到原始输入数据Xl、X2>…、X15，将数据库模块中保存的 VipWj、!!!根据公式(7)和⑶计算，即可得到所需要的经过归一化的目标输出数据d，采用 数据库模块中保存的X' Clmin和X' Clmax经过反归一化式(9)后就可得到目标输出煤耗数据 X’ 0，用来作为机组经济性指标输出。
权利要求
一种混烧低热值煤气机组煤耗计算方法，其特征在于该方法包括下述步骤；步骤一、原始数据准备对电厂混烧低热值煤气机组进行热力试验，获得各种发电功率下的有54个测点数据和用传统的热平衡计算方法计算得到的发电煤耗数据；步骤二、人工控制下的在线学习建模(1)原始数据选取选用步骤一中计算得到的发电煤耗数据作为目标输出煤耗数据，选择15个测点的数据作为原始数据，选取相同时刻的一组目标输出煤耗数据和原始数据作为一个样本，选取能够覆盖所有运行功率的N组(N≥250)样本构成一个原始数据集。(2)原始数据集归一化采用公式(1)对原始数据集进行归一化计算 <mrow><msub> <mi>x</mi> <mi>ip</mi></msub><mo>=</mo><mfrac> <mrow><msubsup> <mi>x</mi> <mi>ip</mi> <mo>′</mo></msubsup><mo>-</mo><msubsup> <mi>x</mi> <mrow><mi>i</mi><mi>min</mi> </mrow> <mo>′</mo></msubsup> </mrow> <mrow><msubsup> <mi>x</mi> <mrow><mi>i</mi><mi>max</mi> </mrow> <mo>′</mo></msubsup><mo>-</mo><msubsup> <mi>x</mi> <mrow><mi>i</mi><mi>min</mi> </mrow> <mo>′</mo></msubsup> </mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo></mrow> </mrow>(1)式中x’ip表示原始数据集中第p组样本第i个测点的原始数据，xip为x’ip归一化后的值，称为原始输入数据；i＝1，2，...，15，表示第i个测点的原始数据；i＝0时，x’0p表示原始数据集中第p组样本目标输出煤耗数据，目标输出数据x0p为x’0p归一化后的值；p＝1，2，...，N，表示第p组原始数据集的样本，N为原始数据集的样本总数；x’imin、x’imax分别表示原始数据集中第i个测点所测得的N个数据列中的最小和最大值；(3)初始化权矩阵Vij和WjVij是m×15阶矩阵，Wj是1×m阶矩阵，开始学习时对Vij和Wj矩阵内每个元素赋初始值为0～1之间的随机数；其中m为该模型的中间隐层节点数；(4)针对p＝1，2，...，N，学习训练，获得权矩阵Vij和Wj①计算中间隐层节点ypj；对于第p组学习样本，有 <mrow><msub> <mi>y</mi> <mi>pj</mi></msub><mo>=</mo><mi>f</mi><mrow> <mo>(</mo> <msub><mi>λ</mi><mi>pj</mi> </msub> <mo>)</mo></mrow><mo>=</mo><mfrac> <mn>2</mn> <mrow><mn>1</mn><mo>+</mo><msup> <mi>e</mi> <msub><mrow> <mo>-</mo> <mi>λ</mi></mrow><mi>pj</mi> </msub></msup> </mrow></mfrac><mo>-</mo><mn>1</mn> </mrow> j＝1，2，，m，i＝1，...，l；(2)其中 <mrow><msub> <mi>λ</mi> <mi>pj</mi></msub><mo>=</mo><munderover> <mi>Σ</mi> <mrow><mi>i</mi><mo>=</mo><mn>1</mn> </mrow> <mi>l</mi></munderover><msub> <mi>V</mi> <mi>ij</mi></msub><msub> <mi>x</mi> <mi>ip</mi></msub> </mrow>②计算中间目标数据dp <mrow><msub> <mi>d</mi> <mi>p</mi></msub><mo>=</mo><mi>f</mi><mrow> <mo>(</mo> <msub><mi>λ</mi><mi>p</mi> </msub> <mo>)</mo></mrow><mo>=</mo><mfrac> <mn>2</mn> <mrow><mn>1</mn><mo>+</mo><msup> <mi>e</mi> <mrow><mrow> <mo>-</mo> <mi>λ</mi></mrow><mi>p</mi> </mrow></msup> </mrow></mfrac><mo>-</mo><mn>1</mn> </mrow>(3)其中 <mrow><msub> <mi>λ</mi> <mi>p</mi></msub><mo>=</mo><munderover> <mi>Σ</mi> <mrow><mi>j</mi><mo>=</mo><mn>1</mn> </mrow> <mi>m</mi></munderover><msub> <mi>W</mi> <mi>j</mi></msub><msub> <mi>y</mi> <mi>pj</mi></msub> </mrow>③计算误差对每一个中间目标数据和目标输出数据进行误差计算，为 <mrow><msub> <mi>E</mi> <mi>p</mi></msub><mo>=</mo><mfrac> <mn>1</mn> <mn>2</mn></mfrac><msup> <mrow><mo>(</mo><msub> <mi>d</mi> <mi>p</mi></msub><mo>-</mo><msub> <mi>x</mi> <mrow><mn>0</mn><mi>p</mi> </mrow></msub><mo>)</mo> </mrow> <mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo></mrow> </mrow>④调整修正计算获得新的连接权矩阵Vij和Wj，调整方法计算中间变量δ0、δyj <mrow><msub> <mi>δ</mi> <mi>o</mi></msub><mo>=</mo><mfrac> <mn>1</mn> <mn>2</mn></mfrac><mrow> <mo>(</mo> <msub><mi>d</mi><mi>p</mi> </msub> <mo>-</mo> <msub><mi>x</mi><mrow> <mn>0</mn> <mi>p</mi></mrow> </msub> <mo>)</mo></mrow><msup> <mrow><mo>(</mo><mn>1</mn><mo>-</mo><msubsup> <mi>x</mi> <mrow><mn>0</mn><mi>p</mi> </mrow> <mn>2</mn></msubsup><mo>)</mo> </mrow> <mn>2</mn></msup> </mrow> <mrow><msub> <mi>δ</mi> <msub><mi>y</mi><mi>j</mi> </msub></msub><mo>=</mo><mfrac> <mn>1</mn> <mn>2</mn></mfrac><mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msubsup><mi>y</mi><mi>pj</mi><mn>2</mn> </msubsup> <mo>)</mo></mrow><msub> <mi>δ</mi> <mi>o</mi></msub><msub> <mi>w</mi> <mi>j</mi></msub> </mrow>j＝1，2，...，m令p＝p+1，计算ΔWjp＝ηδoy(p 1)j+αΔWj(p 1)ΔVjip＝ηδyxi(p 1)+αΔVji(p 1)j＝1，2，...，m(5)Wj′＝Wj+ΔWjp i＝1.2....，lVji′＝Vji+ΔVjipVij＝Vij’Wj＝Wj’；式中ΔWjp、ΔVjp为计算第p组到第p+1组学习样本时调整权矩阵的增量，ΔWj(p 1)、ΔVj(p 1)为前一组调整的权矩阵的增量，式(5)中的Vij’、Wj’表示第P+1组学习样本的权矩阵，y(p 1)j、xi(p 1)为前一组样本的中间隐层节点和第i个输入传感器的原始输入参数；α、η是为了加快收敛速度而设置的两个系数；⑤判断N组样本全计算结束了吗？如果没有计算完，则重复进行①～④的计算，直到N组样本全计算结束，得到中间目标数据dp，然后计算总输出误差E <mrow><mi>E</mi><mo>=</mo><mfrac> <mn>1</mn> <mn>2</mn></mfrac><munderover> <mi>Σ</mi> <mrow><mi>P</mi><mo>=</mo><mn>1</mn> </mrow> <mi>N</mi></munderover><msup> <mrow><mo>(</mo><msub> <mi>d</mi> <mi>p</mi></msub><mo>-</mo><msub> <mi>x</mi> <mrow><mn>0</mn><mi>p</mi> </mrow></msub><mo>)</mo> </mrow> <mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo></mrow> </mrow>如果E不满足精度要求，则重新针对p＝1，2，...，N，再进行(2)～(6)式的计算，这样重复进行，直到误差达到要求；通过调整各Vij、Wj，可以求得输出的误差平方和E达到要求精度时的Vij、Wj，将该Vij、Wj、m以及x′imin和x′max作为混烧低热值煤气机组煤耗计算模型保存在数据库模块内；步骤三、在线混烧低热值煤气机组煤耗计算实际使用时，将现场实测的15个测点的数据作为混烧低热值煤气机组煤耗计算方法模型的原始数据，即x’1、x’2、…、x’15，采用数据库模块中保存的x′imin和x′imax经过(1)式归一化后得到原始输入数据x1、x2、…、x15，将数据库模块中保存的Vij、Wj、m根据公式(7)和(8)计算，即可得到所需要的经过归一化的目标输出数据d，采用数据库模块中保存的x′0min和x′0max经过反归一化式(9)后就可得到目标输出煤耗数据x’0，用来作为机组经济性指标输出； <mrow><msub> <mi>λ</mi> <mi>j</mi></msub><mo>=</mo><munderover> <mi>Σ</mi> <mrow><mi>i</mi><mo>=</mo><mn>1</mn> </mrow> <mi>l</mi></munderover><msub> <mi>V</mi> <mi>ij</mi></msub><msub> <mi>x</mi> <mi>i</mi></msub> </mrow>j＝1，2，…，m；(7) <mrow><msub> <mi>y</mi> <mi>j</mi></msub><mo>=</mo><mi>f</mi><mrow> <mo>(</mo> <msub><mi>λ</mi><mi>j</mi> </msub> <mo>)</mo></mrow><mo>=</mo><mfrac> <mn>2</mn> <mrow><mn>1</mn><mo>+</mo><msup> <mi>e</mi> <msub><mrow> <mo>-</mo> <mi>λ</mi></mrow><mi>j</mi> </msub></msup> </mrow></mfrac><mo>-</mo><mn>1</mn> </mrow> <mrow><mi>λ</mi><mo>=</mo><munderover> <mi>Σ</mi> <mrow><mi>j</mi><mo>=</mo><mn>1</mn> </mrow> <mi>m</mi></munderover><msub> <mi>W</mi> <mi>j</mi></msub><msub> <mi>y</mi> <mi>j</mi></msub> </mrow>j＝1，2，…，m；(8) <mrow><mi>d</mi><mo>=</mo><mi>f</mi><mrow> <mo>(</mo> <mi>λ</mi> <mo>)</mo></mrow><mo>=</mo><mfrac> <mn>2</mn> <mrow><mn>1</mn><mo>+</mo><msup> <mi>e</mi> <mrow><mo>-</mo><mi>λ</mi> </mrow></msup> </mrow></mfrac><mo>-</mo><mn>1</mn> </mrow>x′0＝d·(x′0max x′0min)+x′0min(9)
2.根据权利要求1所述的在线混烧低热值煤气机组煤耗计算方法，其特征是在步骤 ⑶初始化权矩阵、和^中所述的中间隐层节点数m = 12。
全文摘要
在线混烧低热值煤气机组煤耗计算方法，属于电站节能降耗技术领域。本发明在由现场电站机组及测量系统采集的数据和电厂经济性管理系统之间增加一个混烧低热值煤气机组煤耗计算系统。该系统由数据库模块、人工控制下的在线学习建模模块和在线混烧低热值煤气机组煤耗计算模块组成。本发明不需要复杂的焓熵图计算程序，也不需要对所有数据的测点传感器设置硬件冗余，只是利用人工神经网络建立起输入和输出的非线性映射关系，就可以在线计算出精度较高的发电煤耗。
文档编号F23N1/00GK101929684SQ20101024078
公开日2010年12月29日 申请日期2010年7月30日 优先权日2010年7月30日
发明者廖立, 张东欧, 忻建华, 杜光程, 杨昆, 苏明, 蔡青 申请人:上海交通大学