专利摘要

一种基于数据协调的汽轮机湿蒸汽参数的估计方法，属于汽轮机性能监测领域。具体技术步骤包括：1)利用蒸汽膨胀特性构建约束方程，结合准静态测量数据，建立汽轮机系统的数据协调问题；2)对于干蒸汽，直接计算其焓值；对于湿蒸汽，初步假定湿蒸汽的干度，计算出焓值并带入方程；3)通过循环迭代求解数据协调问题，计算得到湿蒸汽参数的估计值。该方法有效利用汽轮机系统中冗余的测量信息及约束方程，并结合蒸汽膨胀过程的性质，给出更优的湿蒸汽参数的估计值，据此可以绘制真实的蒸汽膨胀过程线，并计算系统热力参数的估计值，有利于实现汽轮机性能的监测。该方法简单易于实现，不增加任何额外设备及成本，高效率、低成本。

权利要求

1.一种基于数据协调的汽轮机湿蒸汽参数的估计方法，其特征在于该方法按如下步骤进行：

1)通过准静态数据的选取方法，从电厂数据库选取准静态工况的测量数据；根据电厂系统的测点信息，选取参与热力系统平衡计算的测量变量x，依次记为x₁，x₂，…x_n，并将测量变量的测量值y依次记为y₁，y₂，…y_n；未测量变量u通过热平衡计算得到，其计算值u^*依次记为u₁^*，u₂^*，…u_p^*，其中，n为所有测量变量的个数，p为所有未测量变量的个数；

2)根据电厂热力系统中测量仪表的精度等级，按照方程(1)评定各个测量变量测量值的不确定度，分别记为σ₁，σ₂，…σ_n；

${\sigma }_i=\frac{{\mathrm{k\xi }}_i}{1.96\sqrt{N_s}}---\left(1\right)$

其中，i代表第i个测量变量，i取1，2，3...n；ξ_i代表置信区间为95％的仪表最大允许误差；N_s代表测量同一变量的仪表个数；k为修正系数；

3)根据汽轮机系统的质量平衡、能量平衡及级组的特性方程，建立关于测量变量x及未测量变量u的等式约束方程组f；并根据汽轮机内蒸汽膨胀过程的性质，建立不等式约束方程组g，如下所示：

f(x，u)＝0(2)

g(x，u)≤0(3)

4)对于汽轮机内的干蒸汽，利用实测的蒸汽压力和温度值计算蒸汽焓值；对于湿蒸汽，在蒸汽干度未知的情况下，初步假定蒸汽干度值X₀，并利用实测的蒸汽压力或温度值计算得到湿蒸汽焓值的初值h₀；

5)结合测量变量的测量值y、测量值的不确定度σ、等式约束方程组f及不等式约束方程组g，带入干蒸汽焓值或湿蒸汽焓值的初值h₀，通过循环迭代过程，利用目标函数 计算得到测量变量的协调值x^*和未测量变量的计算值u^*；x^*依次记为x₁^*，x₂^*，…x_n^*；

6)根据测量变量的协调值和未测量变量的计算值，得到湿蒸汽参数的估计值；据此在热力学性质图上绘制汽轮机内蒸汽实际膨胀的过程线。

2.根据权利要求1所述的一种基于数据协调的汽轮机湿蒸汽参数的估计方法，其特征在于：准静态数据的选取方法包括滑动窗口法和移动平均法。

3.根据权利要求1所述的一种基于数据协调的汽轮机湿蒸汽参数的估计方法，其特征在于：步骤1)中所述的测量变量包括流量、压力、温度和输出功率。

4.根据权利要求1所述的一种基于数据协调的汽轮机湿蒸汽参数的估计方法，其特征在于：级组的特性方程包括伏流格尔公式和等熵效率公式。

5.根据权利要求1所述的一种基于数据协调的汽轮机湿蒸汽参数的估计方法，其特征在于：步骤5)中所述的目标函数 包含最小二乘形式的目标函数或强鲁棒性的目标函数；最小二乘形式的目标函数如方程(5)所示，强鲁棒性的目标函数包括方程(6)所示的目标函数：

其中，x_i^*代表x^*中第i个元素，y_i代表y中第i个元素；Σ代表由测量变量的不确定度构成的协方差矩阵；上标T代表矩阵的转置；上标-1代表矩阵的逆； 代表从第1个到第n个测量变量的求和；c_F代表设定的权重参数。

6.根据权利要求1所述的一种基于数据协调的汽轮机湿蒸汽参数的估计方法，其特征在于：步骤6)中所述的湿蒸汽参数包括蒸汽焓值或蒸汽干度。

7.根据权利要求1所述的一种基于数据协调的汽轮机湿蒸汽参数的估计方法，其特征在于：步骤6)中所述的热力学性质图包括焓熵图。

说明书

技术领域

一种基于数据协调的汽轮机湿蒸汽参数的估计方法，属于汽轮机在线性能监测领域。

背景技术

汽轮机性能对整个电厂的性能、热效率及可靠性来说都非常重要。在汽轮机性能的监测过程中，通常需要计算和分析级组的等熵效率和蒸汽膨胀过程线来反映和评估级组的通流能力和做功能力。精确的在线测量数据是汽轮机性能计算及监测的重要基础，测量数据的误差可能导致计算出不合理的等熵效率值或蒸汽膨胀特性。然而，由于人为的误操作或测量仪表的老化，电厂直接测量得到的原始数据总是包含测量误差。因此，能够降低测量数据不确定度的数据处理方法显得尤为重要。此外，由于蒸汽干度在线测量技术的匮乏，汽轮机内湿蒸汽区的蒸汽干度通常无法直接得知，从而湿蒸汽的焓值无法计算。虽然有一些湿蒸汽参数的估计方法，然而可以直接应用于电厂的简单、快捷、高效的方法还很少，尤其缺乏将湿蒸汽区参数估计和提高数据精度同时结合起来的方法。

发明内容

本发明的目的是提供一种基于数据协调的汽轮机湿蒸汽参数的估计方法，能够对汽轮机湿蒸汽参数(如蒸汽干度、焓值)进行估计，并绘制出汽轮机内蒸汽实际膨胀的过程线，有利于汽轮机的性能计算及监测，同时得到更优的电厂系统热力参数的估计值。

本发明的技术方案是：

1)通过准静态数据的选取方法，从电厂数据库选取准静态工况的测量数据；根据电厂系统的测点信息，选取参与热力系统平衡计算的测量变量x，依次记为x₁,x₂,…x_n，并将测量变量的测量值y依次记为y₁,y₂,…y_n；未测量变量u通过热平衡计算得到，其计算值u^*依次记为u₁^*,u₂^*,…u_p^*，其中，n为所有测量变量的个数，p为所有未测量变量的个数；

2)根据电厂热力系统中测量仪表的精度等级，按照方程(1)评定各个测量变量测量值的不确定度，分别记为σ₁,σ₂,…σ_n；

${\mathrm{\σ}}_i=\frac{{\mathrm{k\ξ}}_i}{1.96\sqrt{N_s}}---\left(1\right)$

其中，i代表第i个测量变量，i取1,2,3…n；ξ_i代表置信区间为95％的仪表最大允许误差；N_s代表测量同一变量的仪表个数；k为修正系数；

3)根据汽轮机系统的质量平衡、能量平衡及级组的特性方程，建立关于测量变量x及未测量变量u的等式约束方程组f；并根据汽轮机内蒸汽膨胀过程的性质，建立不等式约束方程组g，如下所示：

f(x,u)＝0 (2)

g(x,u)≤0 (3)

4)对于汽轮机内的干蒸汽，利用实测的蒸汽压力和温度值计算蒸汽焓值；对于湿蒸汽在蒸汽干度未知的情况下，初步假定蒸汽干度值X₀，并利用实测的蒸汽压力或温度值计算得到湿蒸汽焓值的初值h₀；

5)结合测量变量的测量值y、测量值的不确定度σ、等式约束方程组f及不等式约束方程组g，带入干蒸汽焓值或湿蒸汽焓值的初值h₀，通过循环迭代过程，利用目标函数 计算得到测量变量的协调值x^*和未测量变量的计算值u^*；x^*依次记为x₁^*,x₂^*,…x_n^*；

6)根据测量变量的协调值和未测量变量的计算值，得到湿蒸汽参数的估计值；据此在热力学性质图上绘制汽轮机内蒸汽实际膨胀的过程线。

本发明所述的方法的技术特征在于：所述准静态数据的选取方法包括滑动窗口法和移动平均法。步骤(1)中所述的测量变量包括流量、压力、温度和输出功率。所述级组的特性方程包括伏流格尔公式和等熵效率公式。步骤(6)中所述的湿蒸汽参数包括蒸汽焓值或蒸汽干度。步骤(6)中所述的热力学性质图包括焓熵图。

本发明的技术方案中，其特征还在于：

步骤(5)中所述的目标函数 包含最小二乘形式的目标函数或强鲁棒性的目标函数；最小二乘形式的目标函数如方程(5)所示，强鲁棒性的目标函数包括方程(6)所示的目标函数：

其中，x_i^*代表x^*中第i个元素，y_i代表y中第i个元素；Σ代表由测量变量的不确定度构成的协方差矩阵；上标T代表矩阵的转置；上标-1代表矩阵的逆； 代表从第1个到第n个测量变量的求和；c_F代表设定的权重参数。

本发明具有以下优点及突出性的技术效果：

①在湿蒸汽干度在线测量技术匮乏的现状下，本发明可以在不改变测量仪表和设备的前提下，仅采用数据协调循环迭代结合系统建模的方法实现湿蒸汽参数的估计，低成本，简单可行，高效便捷。

②本发明克服了传统方法不能同时实现湿蒸汽区参数估计和提高数据精度的缺点，利用数据协调方法结合汽轮机膨胀过程特性循环迭代的方法，对湿蒸汽参数进行估计，同时给出更优的热力参数的估计值；根据湿蒸汽参数可以绘制出汽轮机内蒸汽实际膨胀的过程线，有利于汽轮机的性能计算及监测；更优的热力参数的估计值有利于整个电厂系统的性能监测。

附图说明

图1为本发明方法的步骤流程示意图。

图2为实施例中某1000MW火电机组汽轮机系统的示意图。

图3为实施例中应用本发明方法前后分别绘制的蒸汽膨胀过程线图。

具体实施方式

下面结合附图对本发明的原理和具体实施作进一步的说明。

1)通过准静态数据的选取方法，比如滑动窗口法等，从电厂数据库选取汽轮机系统在准静态工况下的测量数据，包括流量、压力、温度和输出功率。图2展示了某1000MW火电机组汽轮机系统的示意图。除了图2中所示的汽轮机系统外，参与实施例中数据协调计算的还有电厂其他的热力系统，比如锅炉系统、给水回热系统和凝结水系统。根据电厂系统的测点信息，选取参与热力系统平衡计算的测量变量x，依次记为x₁,x₂,…x₉₇，其测量值y依次记为y₁,y₂,…y₉₇，未测量变量u通过热平衡计算等方法得到，其计算值u^*依次记为u₁^*,u₂^*,…u₈₄^*，其中，测量变量的总个数为97，未测量变量的总个数为84。

2)根据电厂热力系统中测量仪表的精度等级，按照方程(1)评定各个测量变量测量值的不确定度，分别记为σ₁,σ₂,…σ₉₇；

${\mathrm{\σ}}_i=\frac{{\mathrm{k\ξ}}_i}{1.96\sqrt{N_s}}---\left(1\right)$

其中，i代表第i个测量变量，i取1,2,3…97；ξ_i代表置信区间为95％的仪表最大允许误差；N_s代表测量同一变量的仪表个数，通常为1；k为修正系数，通常为1；

3)根据汽轮机系统的质量平衡、能量平衡及级组的特性方程，建立关于测量变量x及未测量变量u的等式约束方程组f，如下所示：

f(x,u)＝0 (2)

级组的特性方程包括伏流格尔公式和等熵效率公式，分别如方程(3)和方程(4)所示：

伏流格尔公式：

等熵效率公式：

其中， 代表设计工况下蒸汽的折合流量；p_f代表设计工况下级组入口的蒸汽压力；p_d代表设计工况下级组出口的蒸汽压力；T_f代表设计工况下级组入口的蒸汽温度，单位为℃，需要转换成绝对温度K；下标0代表非设计工况下的蒸汽参数；η_s代表级组的等熵效率；h_in、h_out和h_s分别代表级组入口、出口和等熵膨胀的蒸汽焓值；

根据汽轮机内蒸汽膨胀过程的性质，建立热力参数的不等式约束方程组g，如下所示：

g(x,u)≤0 (5)

根据美国机械工程师学会给出的汽轮机内蒸汽膨胀过程的性质，蒸汽膨胀过程线应该平滑，等熵效率应该介于0-1之间且逐级递减，因此得到如下的不等式约束关系：

$0\<\<{\left(\frac{h_{in}-h_{out}}{h_{in}-h_s}\right)}_j\≤1---\left(6\right)$

${\left(\frac{\∂h}{\∂s}\right)}_j\≥{\left(\frac{\∂h}{\∂s}\right)}_{j+1}---\left(7\right)$

其中，j代表汽轮机的第j个级组，以图2为例，j取1,2,3…9；

4)对于汽轮机内的干蒸汽，利用实测的蒸汽压力和温度值，通过水蒸汽热力性质表可以计算蒸汽焓值；对于汽轮机内的湿蒸汽，虽然可以测量得到蒸汽温度及压力，但是无法直接测量蒸汽干度，所以蒸汽焓值无法确定；此时，按照正常运行工况下湿蒸汽干度应遵守的范围，比如X＝0.85～0.95，初步假定蒸汽干度值X₀＝0.9，并结合实测的蒸汽压力或温度值，计算得到湿蒸汽焓值的初值h₀；

5)结合测量变量的测量值y、测量值的不确定度σ、等式约束方程组f及不等式约束方程组g，构建汽轮机系统的数据协调问题；带入湿蒸汽焓值的初值h₀，通过循环迭代过程，利用目标函数 计算得到测量变量的协调值x^*和未测量变量的计算值u^*；

其中，x^*依次记为x₁^*,x₂^*,…x₉₇^*；

目标函数 包含最小二乘形式的目标函数或强鲁棒性的目标函数；最小二乘形式的目标函数如方程(9)所示，强鲁棒性的目标函数包括方程(10)所示的目标函数：

其中，x_i^*代表x^*中第i个元素，y_i代表y中第i个元素；Σ代表由测量变量的不确定度构成的协方差矩阵；上标T代表矩阵的转置；上标-1代表矩阵的逆； 代表从第1个到第97个的求和；c_F代表设定的权重参数。

6)根据测量变量的协调值x₁^*,x₂^*,…x₉₇^*和未测量变量的计算值u₁^*,u₂^*,…u₈₄^*，得到湿蒸汽参数的估计值；并根据汽轮机各级级组的入口和出口蒸汽参数，在热力学性质图上绘制汽轮机内蒸汽实际膨胀的过程线。

在实施例中，数据协调计算后直接得到了湿蒸汽焓值的估计值，再通过下式计算得到真实的蒸汽干度X：

$X=\frac{h_{wet}-h_l}{h_g-h_l}---\left(11\right)$

其中，h_wet代表湿蒸汽的焓值；h_g代表饱和水蒸汽的焓值；h_l代表饱和水的焓值；

实施例中，选取的热力学性质图为焓熵图。

实施例1：

下面以一个某1000MW火电机组的汽轮机系统为例，对本发明进行说明。如图2所示，该汽轮机系统包括一个汽轮机高压缸(HPT1,HPT2)、一个中压缸(IPT1,IPT2)、两个低压缸(LPT1,LPT2,LPT3,LPT4,LPT5，图中只画出一个低压缸)、8个抽蒸汽管道(EP1,EP2,EP3,EP4,EP5,EP6,EP7,EP8)、高压蒸汽管道(Hpipe)和中压蒸汽管道(Ipipe)。此外，参与实施例中数据协调计算的还有电厂其他热力系统，如锅炉系统、给水回热系统和凝结水系统。所有测量变量的个数为97个，未测量变量的个数为84个。

根据汽轮机系统的质量平衡、能量平衡及级组的特性方程，建立关于测量变量及未测量变量的等式约束方程组f，等式方程的总个数为110个。并根据汽轮机内蒸汽膨胀过程的性质，建立热力参数的不等式约束方程组g，不等式方程的个数为10个。不等式方程组见表1所示。

表1汽轮机系统的不等式约束方程组

其中，eta代表级组的等熵效率。

按照正常运行工况下干度X应遵守的范围，比如X＝0.85～0.95，初步假定蒸汽干度值X₀＝0.9，并结合实测的压力或温度值，计算得到湿蒸汽焓值的初值h₀，并将湿蒸汽焓值的初值带入数据协调计算中。

结合测量变量的测量值及其不确定度，等式及不等式约束方程组构建汽轮机系统数据协调问题并进行计算，通过一系列循环迭代过程，最终使得协调值满足等式及不等式约束方程，从而得到了测量变量的协调值，由此可以得到湿蒸汽参数的估计值，并可绘制实际过程的蒸汽膨胀过程线。

图3为实施例中应用本发明方法前后分别绘制的蒸汽膨胀过程线图(焓熵图)。可以看到，在应用本方法前，采用实测数据绘制得到的蒸汽膨胀过程线，出现了等熵效率远远大于1的情况，违背了真实的物理特性。而应用本方法后绘制得到的蒸汽膨胀过程线，满足了等熵效率的约束关系，更符合真实的物理特性，修正了蒸汽膨胀过程线。

同时，还可以计算湿蒸汽参数的估计值的不确定度。以汽轮机低压缸的排汽焓值及排汽湿度为例，分别计算基于测量数据和协调数据的不确定度，结果显示，采用本方法后排汽焓值的不确定度降低22.63％，蒸汽湿度的不确定度降低15.61％。具体数据见表2所示。

表2数据协调前后排汽焓值及蒸汽湿度的不确定度

以上结果说明由于数据协调方法及冗余测量信息的引入，本方法可以给出更符合物理规律的湿蒸汽参数估计值，并同时降低汽轮机系统未测量的热力参数的不确定度。

基于以上结果，还可以计算整个电厂热力系统的热耗率和发电效率，通过和热力实验值或历史经验值对比，进一步说明本发明方法的效果。

一种基于数据协调的汽轮机湿蒸汽参数的估计方法专利购买费用说明

Q：办理专利转让的流程及所需资料

A：专利权人变更需要办理著录项目变更手续，有代理机构的，变更手续应当由代理机构办理。

1：专利变更应当使用专利局统一制作的“著录项目变更申报书”提出。

2：按规定缴纳著录项目变更手续费。

3：同时提交相关证明文件原件。

4：专利权转移的，变更后的专利权人委托新专利代理机构的，应当提交变更后的全体专利申请人签字或者盖章的委托书。

Q:专利著录项目变更费用如何缴交

A：（1）直接到国家知识产权局受理大厅收费窗口缴纳,（2）通过代办处缴纳，（3）通过邮局或者银行汇款,更多缴纳方式

Q：专利转让变更，多久能出结果

A：著录项目变更请求书递交后，一般1-2个月左右就会收到通知，国家知识产权局会下达《转让手续合格通知书》。