一种基于功率值计算的汽轮机流量曲线辨识及优化方法
【技术领域】
[0001] 本发明涉及汽轮机发电技术领域，尤其涉及一种基于功率值计算的汽轮机流量曲 线辨识及优化方法。
【背景技术】
[0002] 汽轮机是现代火力发电厂中应用最广泛的原动机，其调门阀组的流量特性曲线是 根据原始设计而定，然而许多在役的老机组，尤其是经过DH1改造、汽机通流部分改造或者 经过大修的机组，经常发生调门阀组流量特性设定曲线与实际曲线不一致的情况。汽轮机 调门阀组作为DH1系统的主要执行机构，其流量特性偏差过大会导致汽轮机节流损失加 大、一次调频的响应负荷不足或者过大、AGC响应变慢、阀门切换负荷波动等问题，最终影响 机组的安全稳定运行。
[0003] 在汽轮机机组的运行过程中，DCS系统每日存储大量历史数据，这些数据来源于运 行实际，是对机组运行特性的全面、客观反映，不可改变，具备绝对真实性。如何正确运用历 史数据，提取机组运行特性，使其更好地为运行实际服务，也是摆在工程技术人员面前的一 个课题。
【发明内容】
[0004] 本发明的目的在于提供一种基于功率值计算的汽轮机流量曲线辨识及优化方法， 能够辨识并优化汽轮机调门阀组在特定主汽压力下的实际流量特性曲线，使调门阀组流量 特性设定曲线与实际曲线一致，保证汽轮机机组的安全稳定运行。
[0005] 为实现上述目的，本发明采用的技术方案是：
[0006] -种基于功率值计算的汽轮机流量曲线辨识及优化方法，包括以下步骤:A、在汽 轮机机组设置若干个采样测点，并设定采样周期和采样时长，采集汽轮机机组的实时工况 数据，分析每个采样测点的实时工况数据包含的参数类型是否完整，如果参数类型有缺失 则重新采集，如果参数类型无缺失则进入步骤B ;
[0007] B、在同一采样周期内，若一个或多个采样测点的实时工况数据中参数大小为零， 或者超出受其物理意义约束的数值范围，则清除该采样周期内所有采样测点的实时工况数 据，得到基本采样数据，然后进入步骤C ;
[0008] C、从基本采样数据中选取输出功率值、主汽压力值和总阀位值均连续稳定若干个 采样周期的实时工况数据，并对实时工况数据中各个类型的参数进行均值化处理，得到稳 定采样数据，然后进入步骤D ;
[0009] D、根据得到的稳定采样数据，采用聚类分析方法中的K-中心点算法挖掘出表征 汽轮机机组运行特性的目标数据集，得到采样工况数据，然后进入步骤E ;
[0010] E、通过对汽轮机机组的输出功率方程式求全微分，计算在总阀位值不变的情况 下，将机组的采样工况修正到特定主汽压力下时机组输出功率值的改变量，然后根据采样 工况下机组的输出功率值计算修正到特定主汽压力下时机组的输出功率值，从而辨识出特 定主汽压力下机组总阀位值与输出功率值之间的关系，即为汽轮机调门阀组流量特性曲 线，然后进入步骤F;
[0011] F、提取特定主汽压力值为额定主汽压力的汽轮机调门阀组流量特性曲线，根据机 组负荷调整的快速性和稳定性对机组流量特性曲线的要求，分析额定主汽压力下汽轮机调 门阀组流量特性曲线的线性度，然后采用分段线性方法对线性度不合理的区段进行优化， 得到优化后的机组总阀位值与输出功率值以及各阀阀位值之间的关系。
[0012] 所述的步骤E中汽轮机机组的输出功率方程式为：
[0014] 式中Pi、D分别为米样工况下机组的输出功率标么值和机组进汽量标么值，η ;为 全机相对内效率，为机组理想比焓降，通过对汽轮机机组的输出功率方程式求全微 分，得到机组的采样工况修正到特定主汽压力下时机组输出功率值的改变量A Pi
[0016] 式中ApQ= p' Q-p。，p'。为修正后特定主汽压力，p。为采样工况的主汽压力，当 汽轮机机组的主汽压力变化范围在8MPa以内时，全机相对内效率^视为不变，即偏导数 0 = 〇,汽轮机机组的总阀位值不变时，凝汽式机组或调节级为临界工况的机组进气 量与主汽压力成正比，即偏导数
将水蒸汽视为理想气体，对于非再热整机机组有
[0019] 式中R为理想气体参数，k为水蒸气绝热系数，T。、V。分别为采样工况的温度和比 体积，pz为汽轮机排汽压力，将三个偏导数代入全微分方程，得到
[0021] 故将机组的采样工况修正到特定主汽压力p'。下时，机组的输出功率值P'
[0023] 所述的步骤F中提取辨识得到额定主汽压力下的汽轮机调门阀组流量特性曲线， 进行高阶次多项式拟合和线性化分析，然后采用分段线性方法对汽轮机调门阀组流量特性 曲线线性度不合理的区段进行优化时，根据汽轮机实际运行中机组通流能力与机组输出负 荷的唯一确定关系，固定机组输出功率值与各阀阀位值的对应关系，然后对汽轮机调门阀 组流量特性曲线线性度不合理的区段进行线性优化处理。
[0024] 所述的步骤D中将稳定采样数据按照数值距离分配到相应的簇，目标函数采用平 方误差准则，采用K-中心点算法挖掘出包含150个以上采样工况的目标数据集。
[0025] 所述的步骤C中输出功率值、主汽压力值和总阀位值连续稳定的采样周期为5个 以上。
[0026] 本发明通过采集并挖掘汽轮机机组运行的历史数据，利用基于功率值计算的方法 辨识出汽轮机调门阀组在特定主汽压力下的实际流量特性曲线，并采用分段线性方法对实 际流量特性曲线线性度不合理的区段进行优化，调整机组总阀位值与输出功率值以及各阀 阀位值之间的关系，使优化后的汽轮机调门阀组流量特性曲线与设定曲线一致，改善机组 运行调节品质。
【附图说明】
[0027] 图1为本发明的流程示意图；
[0028] 图2为采用本发明方法辨识得到的汽轮机调门阀组流量特性曲线；
[0029] 图3为分段线性优化后的汽轮机调门阀组流量特性曲线；
[0030] 图4为分段线性优化后的总阀位值与各阀阀位值对应关系。
【具体实施方式】
[0031] 如图1所示，本发明所述的基于功率值计算的汽轮机流量曲线辨识及优化方法， 包括以下步骤：
[0032] A、在汽轮机机组设置若干个采样测点，并设定采样周期和采样时长，采集汽轮机 机组的实时工况数据，分析每个采样测点的实时工况数据包含的参数类型是否完整，如果 参数类型有缺失则重新采集，如果参数类型无缺失则进入步骤B。
[0033] B、针对所有采样测点，若在同一采样周期内，一个或多个采样测点的实时工况数 据中有参数大小为零，即出现空缺工况数据，或参数大小超出受其物理意义约束的数值范 围，即出现异常工况数据时，剔除该采样周期内所有采样测点的实时工况数据，得到基本采 样数据，分析基本采样数据中各个类型参数的最大值和最小值，以确定各个类型参数的值 域范围，然后进入步骤C。
[0034] C、根据汽轮机机组的控制逻辑，从基本采样数据中选取输出功率值、主汽压力值 和总阀位值均连续稳定5个以上采样周期的实时工况数据，亦即输出功率值、主汽压力值 和总阀位值在5个以上连续采样周期内的变化偏差均小于各自的规定精度，然后对实时工 况数据中各个类型的参数进行均值化处理，得到稳定采样数据，然后进入步骤D。
[0035] D、根据得到的稳定采样数据，采用聚类分析方法中的K-中心点算法，将稳定采样 数据按照数值距离分配到相应的簇，目标函数采用平方误差准则，挖掘出包含150个以上 采样工况的目标数据集，得到能够表征汽轮机机组运行特性的采样工况数据，然后进入步 骤E。
[0036] E、通过对汽轮机机组的输出功率方程式求全微分，计算在总阀位值不变的情况 下，将机组的采样工况修正到特定主汽压力下时机组输出功率值的改变量，然后根据采样 工况下机组的输出功率值计算修正到特定主汽压力下时机组的输出功率值，从而辨识出特 定主汽压力下机组总阀位值与输出功率值之间的关系，即为汽轮机调门阀组流量特性曲 线。
[0037] 汽轮机机组的输出功率方程式为：
[0039] 式中Pi、D分别为米样工况下机组的输出功率标么值和机组进汽量标么值，η ;为 全机相对内效率，为机组理想比焓降。通过对汽轮机机组的输出功率方程式求全微 分，得到机组的采样工况修正到特定主汽压力下时机组输出功率值的改变量A Pi
[0041] 式中ApQ= p' Q-p。，p'。为修正后特定主汽压力，p。为采样工况的主汽压力。
[0042] 当汽轮机机组的主汽压力变化范围在8MPa以内时，全机相对内效率n i视为不 变，即偏导数3-1./0? =〇,汽轮机机组的总阀位值不变时，凝汽式机组或调节级为临界工况 的机组进气量与主汽压力成正比，即偏导数，将水蒸汽视为理想气体，对于非再热 整机机组有
[0045] 式中R为理想气体参数，k为水蒸气绝热系数，Y。、V。分别为采样工况的温度和比 体积，pz为汽轮机排汽压力，将三个偏导数代入全微分方程，得到
[0047] 故将机组的采样工况修正到特定主汽压力p'。下时，机组的输出功率值P'
[0049] 由此得到