CN101949716A - 火电机组能耗数据有效性识别及重构方法 - Google Patents
火电机组能耗数据有效性识别及重构方法 Download PDFInfo
- Publication number
- CN101949716A CN101949716A CN 201010275728 CN201010275728A CN101949716A CN 101949716 A CN101949716 A CN 101949716A CN 201010275728 CN201010275728 CN 201010275728 CN 201010275728 A CN201010275728 A CN 201010275728A CN 101949716 A CN101949716 A CN 101949716A
- Authority
- CN
- China
- Prior art keywords
- measuring point
- time data
- operational factor
- parameter
- data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000005265 energy consumption Methods 0.000 title claims abstract description 14
- 238000004364 calculation method Methods 0.000 claims abstract description 13
- 238000000605 extraction Methods 0.000 claims description 80
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 15
- 241001269238 Data Species 0.000 claims description 10
- 238000006467 substitution reaction Methods 0.000 claims description 5
- 230000001172 regenerating effect Effects 0.000 claims description 4
- 238000005259 measurement Methods 0.000 abstract description 3
- 230000002159 abnormal effect Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000010977 unit operation Methods 0.000 description 2
- 230000005856 abnormality Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及火电机组能耗数据有效性识别及重构方法，首先判断周期内采集的现场测点数据与运行参数的对应关系，并进行相应的有效性识别及标志位赋值；然后对运行参数进行热力系统变工况计算，根据各运行参数的热力学相关性，得到变工况运行参数基准值；最后，将各运行参数与变工况运行参数基准值进行比较，对超出合理性范围的运行参数进行标志位处理，完成数据的重构；本发明采用参数状态来表示机组运行参数的品质，可依据参数状态来决定内部的计算流程，进行容错。同时依据机组参数状态对与其相关的计算结果进行结果状态的标注，在保证整个系统运算正常的基础上，能迅速对异常数据进行跟踪、分类，使用户在第一时间发现异常原因。
Description
技术领域
本发明涉及一种火电机组能耗数据有效性识别及重构方法，属于发电数据监测技术领域。
背景技术
燃煤火电机组运行数据的测量和采集是个复杂的过程，需要将传感器所获得的被测状态参量从物理量转化为电量，再经过放大、采样、模数转换等过程最后经DCS（分散控制系统）传送到数据库服务器。整个过程中许多环节都可能引起测量误差：例如传感器的动态误差，机电转换过程中的零位误差，各种电磁干扰，模数转换器噪声等。而我们实际观测到的数据，往往是真实信号和各种各样的干扰或误差噪音等成分叠加在一起的结果。若不进行数据的效性识别及重构，将严重影响建模的效果，甚至出现模型与实际工业过程相矛盾的现象，使系统功能和结果有较大的偏差。
发明内容
本发明的目的是提供一种火电机组能耗数据有效性识别及重构方法，以解决由于测量或者干扰等引起的数据误差影响建模效果，出现模型与实际工业过程相矛盾的问题。
本发明火电机组能耗数据有效性识别及重构方法步骤如下：
（1）判断周期内采集的现场测点数据与运行参数的对应关系，若是运行参数与现场测点数据一一对应，则直接对该运行参数对应的现场测点数据进行有效性识别及标志位赋值；若是运行参数对应多个现场测点数据，则根据预先设定对运行测点对应的各现场测点数据的平均值或和进行有效性识别及标志位赋值；
（2）对运行参数进行热力系统变工况计算，根据各运行参数的热力学相关性，得到变工况运行参数基准值；
（3）将各运行参数与变工况运行参数基准值进行比较，对超出合理性范围的运行参数进行标志位处理，完成数据的重构。
进一步的，所述步骤（1）中对采集的现场测点数据进行有效性识别是对数据的状态正常、超限、中断、样本点少、数据波动大的判断。
进一步的，所述步骤（1）中若是数据状态正常则运行参数的标志位设为0；若是数据超限则运行参数的标志位设为1；若是数据中断则运行参数的标志位设为2；若是数据样本点少则运行参数的标志位设为3；若是数据波动大则运行参数的标志位设为4。
进一步的，所述步骤（1）若是运行参数与现场测点数据一一对应，按规定的计算周期对实时数据库进行检索，未检索到数据，运行参数=Double.MinValue，标志位=2；
按规定的计算周期对实时数据库进行检索，检索到的数据全部超出上下限，运行参数=实时数据算术平均值，参数状态=1；
按规定的计算周期对实时数据库进行检索，有≥1 and ≤20%周期内总个数的实时数据在上下限范围内，运行参数=在上下限范围内实时数据算术平均值，参数状态=3；
按规定的计算周期对实时数据库进行检索，有≥20%周期内总个数的实时数据在上下限范围内，方差≥规定范围，运行参数=在上下限范围内实时数据算术平均值，参数状态=4；
按规定的计算周期对实时数据库进行检索，有≥20%周期内总个数的实时数据在上下限范围内，方差＜规定范围，运行参数=在上下限范围内实时数据算术平均值，参数状态=0。
进一步的，所述步骤（1）若是运行参数对应多个现场测点数据，则根据预先设定对运行测点对应的各现场测点数据的平均值进行有效性识别及标志位赋值是按规定的计算周期对实时数据库进行检索，该参数对应的所有测点均未检索到数据，运行参数=Double.MinValue，参数状态=2；
按规定的计算周期对实时数据库进行检索，检索到的该参数对应的所有测点数据全部超出上下限，运行参数=实时数据算术平均值，参数状态=1；
按规定的计算周期对实时数据库进行检索，有≥1 and ≤20%周期内该参数对应的所有测点总个数的实时数据在上下限范围内，运行参数=在上下限范围内实时数据算术平均值，参数状态=3；
按规定的计算周期对实时数据库进行检索，有≥20%周期内该参数对应的所有测点总个数的实时数据在上下限范围内，方差≥规定范围，运行参数=在上下限范围内实时数据算术平均值，参数状态=4；
按规定的计算周期对实时数据库进行检索，有≥20%周期内该参数对应的所有测点总个数的实时数据在上下限范围内，方差＜规定范围，运行参数=在上下限范围内实时数据算术平均值，参数状态=0。
进一步的，所述步骤（1）的运行参数对应多个现场测点数据，则根据预先设定对运行测点对应的各现场测点数据的和进行有效性识别及标志位赋值是分两步处理：
首先进行单测点平均值计算，在规定的计算周期对实时数据库进行检索，
未检索到数据，单测点平均值=Double.MinValue，状态=2；
检索到的数据全部超出上下限，单测点平均值=实时数据算术平均值，状态=1；
有≥1 and ≤20%周期内总个数的实时数据在上下限范围内，单测点平均值=在上下限范围内实时数据算术平均值，状态=3；
有≥20%周期内总个数的实时数据在上下限范围内，方差≥规定范围，单测点平均值=在上下限范围内实时数据算术平均值，状态=4；
有≥20%周期内总个数的实时数据在上下限范围内，方差＜规定范围，单测点平均值=在上下限范围内实时数据算术平均值，状态=0；
第二步，依据单测点平均值的情况，进行求和计算，
该参数对应的所有单测点平均值状态=2，运行参数=Double.MinValue，参数状态=2；
该参数对应的所有单测点平均值状态=1，运行参数=Σ（单测点实时数据算术平均值），参数状态=1；
该参数对应的各测点中，有部分单测点平均值状态=3，运行参数=Σ（状态≠1、2的单测点实时数据算术平均值），参数状态=3；
该参数对应的各测点中，有部分单测点平均值状态=4 and ≠3，运行参数=Σ（状态≠1、2的单测点实时数据算术平均值），参数状态=4；
该参数对应的各测点中，有部分单测点平均值状态=0 and≠3、4，运行参数=Σ（状态≠1、2的单测点实时数据算术平均值），参数状态=0。
进一步的，所述步骤（2）中通过功率插值得到某一负荷下：
高压缸效率： 
，
中压缸效率：
，
低压缸效率：
，
流压比：
，
加热器端差：
，
其中：P——实时功率值，
P1——插值功率下沿，
P2——插值功率上沿；
取实时各段抽汽压力值，利用功率插值得到的加热器端差Dt，同时考虑各段抽汽压损、给水泵焓升、轴加焓升值。调用饱和压力函数可得到：
n号加热器出口温度：Tlphno = Tns – TTDn，
其中：Tns——该段抽汽压力下对应的饱和水温度，
TTDn——该加热器的上端差；
取实时主汽流量值，进行回热系统热平衡计算，可得到：
第n级抽汽级后流量：
，
其中：Ffw——主流量，
Hhphno——该加热器出口水焓值，
Hhphni——该加热器入口水焓值，
Fd——疏水流量，
Hd——疏水焓值，
Hi——加热器进汽焓值；
根据流压比基准值FD，得到各段抽汽压力计算值：
Pei=Fin/ FD，
其中：Fin——抽汽流量，
FD——流压比；
该压力与最初设定的抽汽压力作比较，如计算得到的各段抽汽压力与最初设定的抽汽压力相差大于0.1%，则利用目前计算得到的抽汽压力重新代入进行迭代计算，直到相差小于0.1%，迭代结束，此时的各段抽汽压力、抽汽温度、加热器进出口水温、疏水温度对应的值作为该工况的基准值。
进一步的，所述步骤（3）中对超出合理性范围的运行参数进行标志位处理，是将标志位设为5。
本发明的火电机组能够管理数据有效性识别及重构方法，采用参数状态来表示机组运行参数的品质，可依据参数状态来决定内部的计算流程，进行容错。同时依据机组参数状态对与其相关的计算结果进行结果状态的标注，在保证整个系统运算正常的基础上，能迅速对异常数据进行跟踪、分类，使用户在第一时间发现异常原因。
附图说明
图1是本发明的流程图；
图2是本发明实施例的流程图。
具体实施方式
本发明的火电机组能耗管理数据有效性识别及重构方法流程图如图1所示，首先判断周期内采集的现场测点数据与运行参数的对应关系，若是运行参数与现场测点数据一一对应，则直接对该运行参数对应的现场测点数据进行有效性识别及标志位赋值；若是运行参数对应多个现场测点数据，则根据预先设定对运行测点对应的各现场测点数据的平均值或和进行有效性识别及标志位赋值；然后对运行参数进行热力系统变工况计算，根据各运行参数的热力学相关性，得到变工况运行参数基准值；最后将各运行参数与变工况运行参数基准值进行比较，对超出合理性范围的运行参数进行标志位处理，完成数据的重构。
下面进行分别介绍，主要包括如下步骤：
1、判断运行参数与现场测点对应关系：
（1）运行参数与现场测点一一对应的数据有效性识别及标志位处理；
按规定的计算周期对实时数据库进行检索，未检索到数据，运行参数=Double.MinValue，参数状态=2；
按规定的计算周期对实时数据库进行检索，检索到的数据全部超出上下限，运行参数=实时数据算术平均值，参数状态=1；
按规定的计算周期对实时数据库进行检索，有≥1 and ≤20%周期内总个数的实时数据在上下限范围内，运行参数=在上下限范围内实时数据算术平均值，参数状态=3；
按规定的计算周期对实时数据库进行检索，有≥20%周期内总个数的实时数据在上下限范围内，方差≥规定范围，运行参数=在上下限范围内实时数据算术平均值，参数状态=4；
按规定的计算周期对实时数据库进行检索，有≥20%周期内总个数的实时数据在上下限范围内，方差＜规定范围，运行参数=在上下限范围内实时数据算术平均值，参数状态=0；。
（2）运行参数对应多个现场测点，取各测点平均值的数据有效性识别及标志位处理；
按规定的计算周期对实时数据库进行检索，该参数对应的所有测点均未检索到数据，运行参数=Double.MinValue，参数状态=2；
按规定的计算周期对实时数据库进行检索，检索到的该参数对应的所有测点数据全部超出上下限，运行参数=实时数据算术平均值，参数状态=1；
按规定的计算周期对实时数据库进行检索，有≥1 and ≤20%周期内该参数对应的所有测点总个数的实时数据在上下限范围内，运行参数=在上下限范围内实时数据算术平均值，参数状态=3；
按规定的计算周期对实时数据库进行检索，有≥20%周期内该参数对应的所有测点总个数的实时数据在上下限范围内，方差≥规定范围，运行参数=在上下限范围内实时数据算术平均值，参数状态=4；
按规定的计算周期对实时数据库进行检索，有≥20%周期内该参数对应的所有测点总个数的实时数据在上下限范围内，方差＜规定范围，运行参数=在上下限范围内实时数据算术平均值，参数状态=0；
（3）运行参数对应多个现场测点，取各测点和的数据有效性识别及标志位处理；
对这种情况需要分两步处理：
首先进行单测点平均值计算，在规定的计算周期对实时数据库进行检索
未检索到数据，单测点平均值=Double.MinValue，状态=2；
检索到的数据全部超出上下限，单测点平均值=实时数据算术平均值，状态=1；
有≥1 and ≤20%周期内总个数的实时数据在上下限范围内，单测点平均值=在上下限范围内实时数据算术平均值，状态=3；
有≥20%周期内总个数的实时数据在上下限范围内，方差≥规定范围，单测点平均值=在上下限范围内实时数据算术平均值，状态=4；
有≥20%周期内总个数的实时数据在上下限范围内，方差＜规定范围，单测点平均值=在上下限范围内实时数据算术平均值，状态=0；
第二步，依据单测点平均值的情况，进行求和计算
该参数对应的所有单测点平均值状态=2，运行参数=Double.MinValue，参数状态=2；
该参数对应的所有单测点平均值状态=1，运行参数=Σ（单测点实时数据算术平均值），参数状态=1；
该参数对应的各测点中，有部分单测点平均值状态=3，运行参数=Σ（状态≠1、2的单测点实时数据算术平均值），参数状态=3；
该参数对应的各测点中，有部分单测点平均值状态=4 and ≠3，运行参数=Σ（状态≠1、2的单测点实时数据算术平均值），参数状态=4；
该参数对应的各测点中，有部分单测点平均值状态=0 and ≠3、4，运行参数=Σ（状态≠1、2的单测点实时数据算术平均值），参数状态=0；
2、首批计算数据获得完毕，对其进行热力系统变工况计算。根据各参数的热力学相关性，得到参数基准值。
读取数据库并通过功率插值，可以得到两部分数据：一部分为该负荷下高中低压缸效率、流压比、加热器端差、过热器减温水流量、再热器减温水流量等变工况计算所需要的数据；另一部分为给水压力、给水温度、主汽压力、主汽温度等不参与变工况计算的数据。
通过功率插值得到某一负荷下：
高压缸效率：
，
中压缸效率：
，
低压缸效率：
，
流压比：
，
加热器端差： ，
其中：P——实时功率值；P1——插值功率下沿；P2——插值功率上沿；
Effh1——高压缸效率下沿；Effh2——高压缸效率上沿；
Effi1——中压缸效率下沿；Effi2——中压缸效率上沿；
Effl1——低压缸效率下沿；Effl2——低压缸效率上沿；
FP1——流压比下沿；FP2——流压比上沿；
Dt1——加热器端差下沿；Dt2——加热器端差上沿。
取实时各段抽汽压力值，利用功率插值得到的加热器端差Dt，同时考虑各段抽汽压损、给水泵焓升、轴加焓升值。调用饱和压力函数可得到：
n号加热器出口温度：Tlphno = Tns – TTDn，
其中：Tns——该段抽汽压力下对应的饱和水温度，
TTDn——该加热器的上端差，
取实时主汽流量值，进行回热系统热平衡计算，可得到：
第n级抽汽级后流量：
，
其中：Ffw——主流量，
Hhphno——该加热器出口水焓值，
Hhphni——该加热器入口水焓值，
Fd——疏水流量，
Hd——疏水焓值，
Hi——加热器进汽焓值，
根据流压比基准值FD，得到各段抽汽压力计算值：
Pei=Fin/ FD，
其中：Fin——抽汽流量，
FD——流压比。
该压力与最初设定的抽汽压力（即：取实时各段抽汽压力值）作比较，如计算得到的各段抽汽压力与最初设定的抽汽压力相差大于0.1%，则利用目前计算得到的抽汽压力重新代入进行迭代计算。直到相差小于0.1%，迭代结束，此时的各段抽汽压力、抽汽温度、加热器进出口水温、疏水温度对应的值作为该工况的基准值。
3、各实际参数与基准值进行比较，超出合理性范围（如5%）的，对该参数进行标志位处理，并完成数据的重构。
读取数据库中合理性判断上限、下限信息，对实时数据进行判断。如超出合理范围的数据，代入第二步中计算出的基准值，同时进行标志位处理，完成数据的重构。
如图2给出了数据有效性识别及重构流程方法具体的流程图，首先获取实时数据，当未检索到数据，即检测到数据个数的和
0时，标志位=2；检索到的数据全部超出上下限Y min 、Y max ，标志位=1；有≥1 and ≤20%周期内总个数的实时数据在上下限范围内，标志位=3；有≥20%周期内总个数的实时数据在上下限范围内，方差
（Y 1 ……Y N ）≥规定范围
，标志位=4；同时，有≥20%周期内总个数的实时数据在上下限范围内，方差＜规定范围，标志位=0，表示数据正常，首批计算数据获得完毕。通过功率插值得到某一负荷下高压缸效率、中压缸效率、低压缸效率、流压比、加热器端差。取实时各段抽汽压力值，利用功率插值得到的加热器端差，同时考虑各段抽汽压损、给水泵焓升、轴加焓升值。调用饱和压力函数可得到n号加热器出口温度，取实时主汽流量值，进行回热系统热平衡计算，可得到第n级抽汽级后流量。根据流压比基准值FD，得到各段抽汽压力计算值。压力与最初设定的抽汽压力作比较，如计算得到的各段抽汽压力与最初设定的抽汽压力相差大于0.1%，则利用目前计算得到的抽汽压力重新代入进行迭代计算。直到相差小于0.1%，迭代结束，此时的各段抽汽压力、抽汽温度、加热器进出口水温、疏水温度对应的值作为该工况的基准值。各实际参数与基准值进行比较，超出合理性范围（如5%）的，对该参数进行标志位处理，并完成数据的重构。
图2中
为检测到数据个数，Y min 为参数下限，Y max 为参数上限，NUM min为20%周期内总个数，
（Y 1 ……Y N ）为参数某一时段内的方差，
为方差的确定范围，Speacal是程序内位识别方便而赋给该参数的一个特殊值，如9999999；Y status 为标志位，N为功率值，XX为合理性判断的边界条件下限，SX为合理性判断的边界条件上限。
表1 默认值计算方法表结构
说明：默认值计算方法表1中信息用于标示某机组有默认值的参数包括哪些、是否需要进行合理性判断、合理性判断的范围等。
表2 默认值表结构
说明：默认值表2中信息为某机组计算所需的参数的与工况相关的一组数据，当实时数据出现异常时，用于容错。
以C250/N300-16.67/537/537/0.4型汽轮机为例，来具体说明火电机组能耗管理数据有效性识别及重构方法。
实时数据：
2010-7-11 08:30:00 主汽压力，16.722；主汽温度，536.58；……#1高加疏水温度，235.4；……电功率，151.21
2010-7-11 08:30:05主汽压力，16.728；主汽温度，536.51；……#1高加疏水温度，235.6；……电功率，151.21
2010-7-11 08:30:10主汽压力，16.728；主汽温度，536.51；……#1高加疏水温度，235.6；……电功率，151.22
开始处理：
1、电功率检查，不超限（0——175）。
2、按规定的计算周期对实时数据库进行检索，检索到数据。（5秒钟一条，1分钟共检索到12条数据———这里的“条”指的是“主汽压力……电功率”的一批数据）。
3、按规定的计算周期对实时数据库进行检索，检索到的数据未全部超出上下限。（主汽压力，16.722，大于下限8，小于上限20；
主汽温度，536.51，大于下限300，小于上限600；
……；
电功率，151.22，大于下限0，小于上限175）。
4、方差小于规定范围。
主汽压力，16.62；主汽温度，535.92；……#1高加疏水温度，234.1；……电功率，150.89。
5、全部测点接收完毕，标志位=0。
6、进行变工况计算：
高压缸效率：
=85.02；
中压缸效率：
=92.3；
低压缸效率：
=87.09；
流压比：
；
高缸排汽点后流压比 243.82；
中压自动主汽门前流压比 241.44；
1段抽汽点后150.13；
2段抽汽点后243.87；
3段抽汽点后464.35；
4段抽汽点后772.1；
5段抽汽点后1325.38；
6段抽汽点后2327.15；
7段抽汽点后5780.65；
8段抽汽点后6024.12；
加热器端差：
；
1号高加上端差-1.07；
1号高加下端差5.6；
2号高加上端差0；
2号高加下端差5.6；
3号高加上端差0；
3号高加下端差5.6；
5号低加上端差2.8；
5号低加下端差5.6；
6号低加上端差2.8；
6号低加下端差5.6；
7号低加上端差2.8；
7号低加下端差5.6；
8号低加上端差2.8；
8号低加下端差5.6；
n号加热器出口温度：Tlphno = Tns – TTDn；
1号高加出口温度283.34；
2号高加出口温度246.89；
3号高加出口温度205.84；
5号低加出口温度149.61；
6号低加出口温度127.59；
7号低加出口温度110.12；
8号低加出口温度75.62；
第n级抽汽级后流量：；
第1级抽汽级后流量953.32；
第2级抽汽级后流量895.53；
第3级抽汽级后流量856.66；
第4级抽汽级后流量824.41；
第5级抽汽级后流量792.54；
第6级抽汽级后流量763.56；
第7级抽汽级后流量747.13；
第8级抽汽级后流量710.23；
Pei=Fin/ FD，
1段抽汽压力6.35，
比较得到#1高加疏水温度实时值6.2与基准值相差2.3%；
2段抽汽压力3.67，
比较得到#1高加疏水温度实时值4.2与基准值相差14%；
3段抽汽压力1.84，
比较得到#1高加疏水温度实时值1.7与基准值相差7.8%；
4段抽汽压力1.07，
比较得到#1高加疏水温度实时值0.93与基准值相差13%；
5段抽汽压力0.59，
比较得到#1高加疏水温度实时值0.51与基准值相差14%；
6段抽汽压力0.32，
比较得到#1高加疏水温度实时值0.28与基准值相差15%；
7段抽汽压力0.13，
比较得到#1高加疏水温度实时值0.1与基准值相差22%；
8段抽汽压力0.11，
比较得到#1高加疏水温度实时值0.09与基准值相差23%；
除1段抽汽压力实时值6.2与基准值相差2.3%外，其余各点均超出4%（收敛条件）的范围。因此需要第二次迭代。将新的各段抽汽压力作为初值，进行第二次的计算，直至各段抽汽压力值均在5%的范围内，迭代结束。
此时，各参数的基准值为：
高压缸效率=84.8；
中压缸效率=92.12；
低压缸效率=87.03；
1段抽汽压力6.21；
2段抽汽压力3.82；
3段抽汽压力1.84；
4段抽汽压力1.12；
5段抽汽压力0.59；
6段抽汽压力0.32；
7段抽汽压力0.12；
8段抽汽压力0.11；
……
变工况计算结束。
将每一段抽汽压力的基准值与实时值相比，对超出合理性范围（如5%）的参数进行标志位处理，赋值为5，完成数据的重构。
1段抽汽压力基准值6.21，实时值6.18，相差0.48%，标志位=0；
2段抽汽压力基准值3.82，实时值3.8，相差0.52%，标志位=0；
3段抽汽压力基准值1.84，实时值1.76，相差4.34%，标志位=0；
4段抽汽压力基准值1.12，实时值1.16，相差3.57%，标志位=0；
5段抽汽压力基准值0.59，实时值0.61，相差3.38%，标志位=0；
6段抽汽压力基准值0.32，实时值0.31，相差3.12%，标志位=0；
7段抽汽压力基准值0.12，实时值0.12，相差0%，标志位=0；
8段抽汽压力基准值0.11，实时值0.09，相差18.18%，标志位=5。
Claims (8)
1.一种火电机组能耗数据有效性识别及重构方法，其特征在于，该方法的步骤如下：
（1）判断周期内采集的现场测点数据与运行参数的对应关系，若是运行参数与现场测点数据一一对应，则直接对该运行参数对应的现场测点数据进行有效性识别及标志位赋值；若是运行参数对应多个现场测点数据，则根据预先设定对运行测点对应的各现场测点数据的平均值或和进行有效性识别及标志位赋值；
（2）对运行参数进行热力系统变工况计算，根据各运行参数的热力学相关性，得到变工况运行参数基准值；
（3）将各运行参数与变工况运行参数基准值进行比较，对超出合理性范围的运行参数进行标志位处理，完成数据的重构。
2.根据权利要求1所述的火电机组能耗数据有效性识别及重构方法，其特征在于：所述步骤（1）中对采集的现场测点数据进行有效性识别是对数据的状态正常、超限、中断、样本点少、数据波动大的判断。
3.根据权利要求2所述的火电机组能耗数据有效性识别及重构方法，其特征在于：所述步骤（1）中若是数据状态正常则运行参数的标志位设为0；若是数据超限则运行参数的标志位设为1；若是数据中断则运行参数的标志位设为2；若是数据样本点少则运行参数的标志位设为3；若是数据波动大则运行参数的标志位设为4。
4.根据权利要求3所述的火电机组能耗数据有效性识别及重构方法，其特征在于：所述步骤（1）若是运行参数与现场测点数据一一对应，按规定的计算周期对实时数据库进行检索，未检索到数据，运行参数=Double.MinValue，标志位=2；
按规定的计算周期对实时数据库进行检索，检索到的数据全部超出上下限，运行参数=实时数据算术平均值，参数状态=1；
按规定的计算周期对实时数据库进行检索，有≥1 and ≤20%周期内总个数的实时数据在上下限范围内，运行参数=在上下限范围内实时数据算术平均值，参数状态=3；
按规定的计算周期对实时数据库进行检索，有≥20%周期内总个数的实时数据在上下限范围内，方差≥规定范围，运行参数=在上下限范围内实时数据算术平均值，参数状态=4；
按规定的计算周期对实时数据库进行检索，有≥20%周期内总个数的实时数据在上下限范围内，方差＜规定范围，运行参数=在上下限范围内实时数据算术平均值，参数状态=0。
5.根据权利要求1所述的火电机组能耗数据有效性识别及重构方法，其特征在于：所述步骤（1）若是运行参数对应多个现场测点数据，则根据预先设定对运行测点对应的各现场测点数据的平均值进行有效性识别及标志位赋值是按规定的计算周期对实时数据库进行检索，该参数对应的所有测点均未检索到数据，运行参数=Double.MinValue，参数状态=2；
按规定的计算周期对实时数据库进行检索，检索到的该参数对应的所有测点数据全部超出上下限，运行参数=实时数据算术平均值，参数状态=1；
按规定的计算周期对实时数据库进行检索，有≥1 and ≤20%周期内该参数对应的所有测点总个数的实时数据在上下限范围内，运行参数=在上下限范围内实时数据算术平均值，参数状态=3；
按规定的计算周期对实时数据库进行检索，有≥20%周期内该参数对应的所有测点总个数的实时数据在上下限范围内，方差≥规定范围，运行参数=在上下限范围内实时数据算术平均值，参数状态=4；
按规定的计算周期对实时数据库进行检索，有≥20%周期内该参数对应的所有测点总个数的实时数据在上下限范围内，方差＜规定范围，运行参数=在上下限范围内实时数据算术平均值，参数状态=0。
6.根据权利要求1所述的火电机组能耗数据有效性识别及重构方法，其特征在于：所述步骤（1）的运行参数对应多个现场测点数据，则根据预先设定对运行测点对应的各现场测点数据的和进行有效性识别及标志位赋值是分两步处理：
首先进行单测点平均值计算，在规定的计算周期对实时数据库进行检索，
未检索到数据，单测点平均值=Double.MinValue，状态=2；
检索到的数据全部超出上下限，单测点平均值=实时数据算术平均值，状态=1；
有≥1 and ≤20%周期内总个数的实时数据在上下限范围内，单测点平均值=在上下限范围内实时数据算术平均值，状态=3；
有≥20%周期内总个数的实时数据在上下限范围内，方差≥规定范围，单测点平均值=在上下限范围内实时数据算术平均值，状态=4；
有≥20%周期内总个数的实时数据在上下限范围内，方差＜规定范围，单测点平均值=在上下限范围内实时数据算术平均值，状态=0；
第二步，依据单测点平均值的情况，进行求和计算，
该参数对应的所有单测点平均值状态=2，运行参数=Double.MinValue，参数状态=2；
该参数对应的所有单测点平均值状态=1，运行参数=Σ（单测点实时数据算术平均值），参数状态=1；
该参数对应的各测点中，有部分单测点平均值状态=3，运行参数=Σ（状态≠1、2的单测点实时数据算术平均值），参数状态=3；
该参数对应的各测点中，有部分单测点平均值状态=4 and ≠3，运行参数=Σ（状态≠1、2的单测点实时数据算术平均值），参数状态=4；
该参数对应的各测点中，有部分单测点平均值状态=0 and≠3、4，运行参数=Σ（状态≠1、2的单测点实时数据算术平均值），参数状态=0。
7.根据权利要求1所述的火电机组能耗数据有效性识别及重构方法，其特征在于：所述步骤（2）中通过功率插值得到某一负荷下：
高压缸效率： 
，
中压缸效率： ，
低压缸效率：
，
流压比：
，
加热器端差：
，
其中：P——实时功率值，
P1——插值功率下沿，
P2——插值功率上沿；
取实时各段抽汽压力值，利用功率插值得到的加热器端差Dt，同时考虑各段抽汽压损、给水泵焓升、轴加焓升值；调用饱和压力函数可得到：
n号加热器出口温度：Tlphno = Tns – TTDn，
其中：Tns——该段抽汽压力下对应的饱和水温度，
TTDn——该加热器的上端差；
取实时主汽流量值，进行回热系统热平衡计算，可得到：
第n级抽汽级后流量：
，
其中：Ffw——主流量，
Hhphno——该加热器出口水焓值，
Hhphni——该加热器入口水焓值，
Fd——疏水流量，
Hd——疏水焓值，
Hi——加热器进汽焓值；
根据流压比基准值FD，得到各段抽汽压力计算值：
Pei=Fin/ FD，
其中：Fin——抽汽流量，
FD——流压比；
该压力与最初设定的抽汽压力作比较，如计算得到的各段抽汽压力与最初设定的抽汽压力相差大于0.1%，则利用目前计算得到的抽汽压力重新代入进行迭代计算，直到相差小于0.1%，迭代结束，此时的各段抽汽压力、抽汽温度、加热器进出口水温、疏水温度对应的值作为该工况的基准值。
8.根据权利要求1所述的火电机组能耗数据有效性识别及重构方法，其特征在于：所述步骤（3）中对超出合理性范围的运行参数进行标志位处理，是将标志位设为5。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2010102757287A CN101949716B (zh) | 2010-09-08 | 2010-09-08 | 火电机组能耗数据有效性识别及重构方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2010102757287A CN101949716B (zh) | 2010-09-08 | 2010-09-08 | 火电机组能耗数据有效性识别及重构方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101949716A true CN101949716A (zh) | 2011-01-19 |
| CN101949716B CN101949716B (zh) | 2012-07-25 |
Family
ID=43453294
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2010102757287A Expired - Fee Related CN101949716B (zh) | 2010-09-08 | 2010-09-08 | 火电机组能耗数据有效性识别及重构方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101949716B (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102542109A (zh) * | 2011-12-29 | 2012-07-04 | 西安交通大学 | 一种考虑负荷特性的火电机组热力系统优化设计方法 |
| CN103034209A (zh) * | 2012-12-18 | 2013-04-10 | 湖南大唐先一科技有限公司 | 一种在线测量数据准确性甄别方法 |
| CN104574729A (zh) * | 2015-01-28 | 2015-04-29 | 无锡物联网产业研究院 | 一种报警方法、装置及系统 |
| CN112799447A (zh) * | 2021-01-06 | 2021-05-14 | 国家电投集团河南电力有限公司 | 一种高压加热器的控制方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006065646A (ja) * | 2004-08-27 | 2006-03-09 | Toshiba Corp | 警報評価方法および装置 |
| CN101187804A (zh) * | 2006-11-15 | 2008-05-28 | 北京华电天仁电力控制技术有限公司 | 基于数据挖掘的火电机组运行优化规则提取方法 |
| CN101725998A (zh) * | 2009-12-14 | 2010-06-09 | 贵州电力试验研究院 | 煤耗在线监测系统中异常数据的判定及替换系统 |
| CN101738972A (zh) * | 2009-12-14 | 2010-06-16 | 贵州电力试验研究院 | 检测煤耗在线监测系统监测准确性的测试方法 |
-
2010
- 2010-09-08 CN CN2010102757287A patent/CN101949716B/zh not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006065646A (ja) * | 2004-08-27 | 2006-03-09 | Toshiba Corp | 警報評価方法および装置 |
| CN101187804A (zh) * | 2006-11-15 | 2008-05-28 | 北京华电天仁电力控制技术有限公司 | 基于数据挖掘的火电机组运行优化规则提取方法 |
| CN101725998A (zh) * | 2009-12-14 | 2010-06-09 | 贵州电力试验研究院 | 煤耗在线监测系统中异常数据的判定及替换系统 |
| CN101738972A (zh) * | 2009-12-14 | 2010-06-16 | 贵州电力试验研究院 | 检测煤耗在线监测系统监测准确性的测试方法 |
Non-Patent Citations (1)
| Title |
|---|
| 《中国电机工程学报》 20040217 于达仁等 基于分布信息融合的直流锅炉燃料量信号重构 第191-195页 1-8 第24卷, 第2期 2 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102542109A (zh) * | 2011-12-29 | 2012-07-04 | 西安交通大学 | 一种考虑负荷特性的火电机组热力系统优化设计方法 |
| CN103034209A (zh) * | 2012-12-18 | 2013-04-10 | 湖南大唐先一科技有限公司 | 一种在线测量数据准确性甄别方法 |
| CN103034209B (zh) * | 2012-12-18 | 2016-02-17 | 湖南大唐先一科技有限公司 | 一种在线测量数据准确性甄别方法 |
| CN104574729A (zh) * | 2015-01-28 | 2015-04-29 | 无锡物联网产业研究院 | 一种报警方法、装置及系统 |
| CN104574729B (zh) * | 2015-01-28 | 2020-07-24 | 无锡感知金服物联网科技有限公司 | 一种报警方法、装置及系统 |
| CN112799447A (zh) * | 2021-01-06 | 2021-05-14 | 国家电投集团河南电力有限公司 | 一种高压加热器的控制方法 |
| CN112799447B (zh) * | 2021-01-06 | 2022-03-04 | 国家电投集团河南电力有限公司 | 一种高压加热器的控制方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101949716B (zh) | 2012-07-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| WO2017050207A1 (zh) | 发电机组能效分析方法 | |
| CN111159844B (zh) | 一种电站燃气轮机排气温度的异常检测方法 | |
| CN105317476B (zh) | 一种基于特征通流面积的汽轮机流量曲线辨识及优化方法 | |
| CN109211439B (zh) | 一种汽轮机低压缸排汽焓值在线监测系统及方法 | |
| CN101949716B (zh) | 火电机组能耗数据有效性识别及重构方法 | |
| CN106936627A (zh) | 一种基于大数据分析挖掘的火电设备性能监测方法 | |
| CN108463622A (zh) | 设施分析装置、设施分析方法以及程序 | |
| CN112417764A (zh) | 一种面向锅炉特种设备蒸汽量预测的k近邻回归预测方法 | |
| CN108153255B (zh) | 一种基于dcs的火电机组性能监测方法及装置 | |
| JP3614640B2 (ja) | 火力発電プラントの熱効率診断方法および装置 | |
| CN115014454A (zh) | 一种火电机组主蒸汽流量软测量方法、系统、设备及介质 | |
| JP6803788B2 (ja) | 情報処理装置、情報処理方法およびプログラム | |
| RU2621422C2 (ru) | Система и способ тестирования показателя работы паровой турбины | |
| CN111027186B (zh) | 汽水分离再热器运行状态下性能诊断方法及系统 | |
| CN111723331B (zh) | 一种联合循环二拖一机组汽轮机负荷的权益分配计算方法 | |
| Li et al. | Study on the effect of regenerative system on power type relative internal efficiency of nuclear steam turbine | |
| CN104459542A (zh) | 再热回热复合循环机组的热耗率测定方法 | |
| JP2015230576A (ja) | プラント診断用データ作成システム | |
| CN116702373A (zh) | 一种汽轮机在线寻优及分析系统 | |
| TWI352360B (zh) | ||
| KR100752765B1 (ko) | 복합화력 발전소의 운전제어가능 파라미터의 실시간영향산출 시스템 및 방법 | |
| CN110032791B (zh) | 基于广义回归神经网络的汽轮机低压缸效率实时计算方法 | |
| CN111581787B (zh) | 一种实时筛选汽轮机热耗率分析数据的方法及系统 | |
| CN109189029B (zh) | 火力发电厂低温省煤器节能量在线监测系统及方法 | |
| CN208578597U (zh) | 汽轮发电机组回热系统加热器性能指标监测系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| ASS | Succession or assignment of patent right |
Owner name: ELECTRIC POWER RESEARCH INSTITUTE OF HEBEI ELECTRI Free format text: FORMER OWNER: HENAN INSTITUTE OF POWER EXPERIMENT Effective date: 20121206 Owner name: STATE GRID CORPORATION OF CHINA Effective date: 20121206 |
|
| C41 | Transfer of patent application or patent right or utility model | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20121206 Address after: 450052 Songshan South Road, Henan, No. 85, No. Patentee after: Henan Electric Power Corporation Electric Power Science Research Institute Patentee after: State Grid Corporation of China Address before: 450052 Songshan South Road, Henan, No. 85, No. Patentee before: Henan Institute of Power Experiment |
|
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120725 Termination date: 20150908 |
|
| EXPY | Termination of patent right or utility model |