CN110821576A - 一种基于小样本暂态数据的汽轮机高调门流量特性优化方法 - Google Patents

一种基于小样本暂态数据的汽轮机高调门流量特性优化方法 Download PDF

Info

Publication number
CN110821576A
CN110821576A CN201911132428.0A CN201911132428A CN110821576A CN 110821576 A CN110821576 A CN 110821576A CN 201911132428 A CN201911132428 A CN 201911132428A CN 110821576 A CN110821576 A CN 110821576A
Authority
CN
China
Prior art keywords
steam
steam turbine
data
small sample
valve position
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.)
Pending
Application number
CN201911132428.0A
Other languages
English (en)
Inventor
霍红岩
于海存
辛晓钢
张谦
秦成果
张国斌
殷建华
范志强
周磊
郭瑞君
杜荣华
赵炜
王立红
李荣丽
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Inner Mongolia Electric Power Research Institute of Inner Mongolia Power Group Co Ltd
Original Assignee
Inner Mongolia Electric Power Research Institute of Inner Mongolia Power Group Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Inner Mongolia Electric Power Research Institute of Inner Mongolia Power Group Co Ltd filed Critical Inner Mongolia Electric Power Research Institute of Inner Mongolia Power Group Co Ltd
Priority to CN201911132428.0A priority Critical patent/CN110821576A/zh
Publication of CN110821576A publication Critical patent/CN110821576A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
    • F01D—NON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
    • F01D17/00—Regulating or controlling by varying flow
    • F01D17/10—Final actuators
    • F01D17/12—Final actuators arranged in stator parts
    • F01D17/14—Final actuators arranged in stator parts varying effective cross-sectional area of nozzles or guide conduits
    • F01D17/141—Final actuators arranged in stator parts varying effective cross-sectional area of nozzles or guide conduits by means of shiftable members or valves obturating part of the flow path
    • F01D17/145—Final actuators arranged in stator parts varying effective cross-sectional area of nozzles or guide conduits by means of shiftable members or valves obturating part of the flow path by means of valves, e.g. for steam turbines
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
    • F01D—NON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
    • F01D17/00—Regulating or controlling by varying flow
    • F01D17/20—Devices dealing with sensing elements or final actuators or transmitting means between them, e.g. power-assisted
    • F01D17/22—Devices dealing with sensing elements or final actuators or transmitting means between them, e.g. power-assisted the operation or power assistance being predominantly non-mechanical
    • F01D17/24—Devices dealing with sensing elements or final actuators or transmitting means between them, e.g. power-assisted the operation or power assistance being predominantly non-mechanical electrical
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
    • F01D—NON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
    • F01D21/00—Shutting-down of machines or engines, e.g. in emergency; Regulating, controlling, or safety means not otherwise provided for
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
    • F01D—NON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
    • F01D21/00—Shutting-down of machines or engines, e.g. in emergency; Regulating, controlling, or safety means not otherwise provided for
    • F01D21/003—Arrangements for testing or measuring

Landscapes

  • Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Control Of Turbines (AREA)

Abstract

本发明公开了一种基于小样本暂态数据的汽轮机高调门流量特性优化方法，通过采样汽轮机变负荷过程中的暂态数据，再对所采集的数据进行清洗和聚类分析，得到各个负荷工况下的样本数据集；根据汽轮机调节级压力和进汽量之间的正比关系，对样本数据集中的汽轮机进汽量进行修正，得到相同边界条件下汽轮机进汽量；最后采用最小二乘拟合算法，对进汽流量与总阀门指令进行线性化拟合，并结合原有流量与阀门指令之间的关系，对汽轮机高调门流量特性曲线进行修正，得到线性化后的汽轮机各高调门配汽曲线。

Description

一种基于小样本暂态数据的汽轮机高调门流量特性优化方法
技术领域
本发明专利属于汽轮机运行优化技术领域，涉及一种基于小样本暂态 数据的汽轮机高调门流量特性优化技术。以机组小样本暂态数据为依据，通过 数据清洗、聚类分析及蒸汽流量修正等技术，对汽轮机高调门线性度进行优化， 可以有效改善汽轮机的动态调节性能，提升机组负荷控制精度和一次调频控制 品质。
背景技术
随着风能、太阳等新能源电力的大量接入，为了保证电网负荷和频率的稳 定，平抑随机波动性电源对电网安全造成的影响，常规燃煤发电机组需频繁的 参与到电网调峰和调频中，电网对火电机组的控制要求不断提高，机组需能够 快速跨出调节死区，并实现快速高精度的负荷调节。
机组负荷控制通常由汽轮机数字电液控制系统(DEH系统)来完成，DEH 系统的阀门流量特性曲线线性度偏差过大，其不仅会增加机组对负荷的响应时 间，而且在快速调节过程中容易造成高调门的频繁震荡，导致机组负荷的调节 速率和调节精度性能变差、主蒸汽参数波动，阀门执行机构的疲劳损伤增加， 影响电网两个细则的考核指标，不利于机组安全、经济运行。
在机组实际运行过程中，其每天频繁的参与负荷调节，在分散控制系统 (DCS系统)中存在大量机组暂态数据，这些数据能够客观、全面的反映出机 组当前的运行特性和性能。如何从小样本的暂态数据中分析汽轮机高调门的流 量特性，使其更好地为运行实际服务，也是摆在工程技术人员面前的一个课题。
发明内容
为解决上述问题，本发明专利提供了一种基于小样本暂态数据的汽轮机高 调门流量特性优化技术。本发明专利的目的在于提供一种基于小样本暂态数据 的汽轮机高调门流量特性优化技术，通过小样本的机组暂态数据，辨识并优化 汽轮机总阀位指令与主蒸汽流量之间的线性关系，使得汽轮机高调门配汽函数 更加合理，提升机组的变负荷动态性能。
为实现上述目的，本发明专利采用的技术方案是：
第一步：小样本暂态数据采集。利用DCS中存储的历史数据，选取机组负 荷在50％-100％区间内的连续变负荷数据段；
第二步：数据清洗。采用孤立深林算法对小样本数据中的异常值进行检查 并剔除，并采用线性差值方法进行补充；
第三步：数据聚类。以机组负荷作为质心划分的依据，每变化1MW为一个 质心，采用PAM聚类算法对清洗后的数据进行聚类计算，得到表征汽轮机运 行特性的目标数据集；
第四步：主蒸汽流量修正。选取目标数据集中最高负荷工况下对应的主蒸 汽压力和调节级压力作为参考值，根据机组调节级压力和蒸汽流量之间的正比 关系，将目标数据集中汽轮机进汽量修正至参考值下对用的汽轮机进汽量；
第五步：阀门流量特性线性化拟合。根据修正后的汽轮机进汽量与总阀位 指令，采用最小二次算法进行线性化拟合，得到期望的汽轮机进汽量与总阀位 指令间的线性关系；
第六步：高调门配汽曲线修正。依据拟合得到的汽轮机进汽量与总阀位之 前的线性曲线，结合原始的汽轮机高阀门配汽函数，对总阀位指令分配至各单 阀指令函数进行修正，得到新的汽轮机阀门配汽函数。
附图说明
图1为本发明实施例实际流量特性曲线与修正后流量特性曲线。
图2为本发明实施例优化后总阀位指令与GV1阀位指令对应关系曲线示例 图。
具体实施方式
本发明专利所述的一种基于小样本暂态数据的汽轮机高调门流量特性优化 技术，包括以下步骤：
第一步：小样本暂态数据采集。利用DCS中存储的历史数据，选取机组负 荷在50％-100％区间内的连续变负荷数据段；如表1所示，本发明专利采集参数 包括：机组功率、主蒸汽流量、主蒸汽压力、主蒸汽温度、调节级压力、高压 缸排汽压力、高压缸排汽温度、总阀位开度、GV1阀位开度、GV2阀位开度、 GV3阀位开度、GV4阀位开度。
采用周期：1秒。
第二步：采用孤立深林算法对小样本数据中的异常值进行检查并剔除，并 采用线性差值方法进行补充。
孤立森林由iTree组成，采用多次迭代方式构建二叉搜索树，并将这些二叉 树组成森林。其构建过程如下：
(1)从n个训练数据中随机抽取m个子样本，放入一颗孤立树的根节点；
(2)随机指定一个维度，在当前节点数据范围内，随机产生一个切割点p—— 切割点产生于当前节点数据中指定维度的最大值与最小值之间；
(3)此切割点的选取生成了一个超平面，将当前节点数据空间切分为2个 子空间：把当前所选维度下小于p的点放在当前节点的左分支，把大于等于p 的点放在当前节点的右分支；
(4)把节点的左分支和右分支节点递归步骤2、3，不断构造新的叶子节点， 直到叶子节点上只有一个数据(无法再继续切割)或树已经生长到了所设定的 高度。
(5)记录每个数据从根节点到叶子节点经过边的数量，即路径长度h(x)， 以此来计算该数据的异常指数S(x,n):
c(n)＝2H(n-1)-(2(n-1)/n)
其中，H(n)＝ln(n)+ξ,ξ＝0.5772156649为欧拉常数。S(x,n)越接近1表示该点是异 常点的可能性越高，越接近0表示该点为正常点的可能性越高。
第三步：数据聚类。以机组负荷作为质心划分的依据，每变化1MW为一个 质心，随机抽取初始聚类中心，并采用PAM聚类算法对清洗后的数据进行聚 类计算，得到表征汽轮机运行特性的目标数据集。以机组负荷变化范围为 [150MW,300MW]为例，以1MW间隔作为质心划分依据，确定聚类质心为150 个，采用PAM聚类算法获取目标数据集；
第四步：汽轮机进汽流量修正。选取目标数据集中最高机组负荷工况下对 应的主蒸汽压力和调节级压力作为参考值，将目标数据集中的蒸汽流量修正至 同一蒸汽压力水平下，计算方法如下：
其中，G、p分别表示汽轮机进汽量t/h、压力MPa；下标A表示参考工况； 小标B表示实际工况；下标0表示主蒸汽参数；下标1表示调节级蒸汽参数。
第五步：高调门流量线性化拟合。根据修正后的主蒸汽流量特征值RealValue 与总阀位指令，采用最小二次算法进行线性化拟合，得到期望的汽轮机进汽量 与总阀位指令间的线性关系，如图2所示。
第六步：高调门配汽曲线修正。依据拟合得到的汽轮机进汽流量与总阀位 之前的线性曲线，结合原始的汽轮机高调门配汽函数，对总阀位指令分配至各 单阀指令函数进行修正，得到新的汽轮机高调门配汽函数。如图2所示为优化 后总阀位指令与GV1阀位指令对应关系曲线。

Claims (5)

1.一种基于小样本暂态数据的汽轮机高调门流量特性优化方法，其特征在于，包括以下步骤：
步骤一：小样本暂态数据采集，基于DCS系统中存储的历史数据，以秒级采样周期，采集机组在50％-100％负荷工况下连续变负荷的小样本暂态数据；
步骤二：历史数据清洗，采用孤立森林算法对步骤一所采集的小样本暂态数据中的异常值进行检查并剔除，并采用线性差值方法进行补充；
步骤三：历史数据聚类，以机组负荷作为聚类质心的划分依据，从步骤二清洗后的历史数据中，根据负荷的变化幅度，每隔1MW选择初始聚类中心点，采用PAM聚类算法得到不同负荷工况下的机组蒸汽参数及阀位指令的目标数据集；
步骤四：汽轮机主蒸汽进汽量修正，选取步骤三所得到的目标数据集中最高负荷工况下对应的主蒸汽压力和调节级作为参考值，根据机组调节级压力和蒸汽流量之间的正比关系，将目标数据集中汽轮机进汽量修正至参考值下对应的汽轮机进汽量；
步骤五：高调门配汽曲线优化，将步骤四修正后的汽轮机进汽量与总阀位指令采用最小二乘算法进行线性化拟合，得到期望的汽轮机进汽量与总阀位指令之间的线性关系曲线；
步骤六：高调门配汽曲线修正，利用步骤五拟合得到的汽轮机进汽量与对应的总阀位指令之间的线性关系，结合原汽轮机高调门流量特性曲线，对总阀位指令分配至各单阀的指令函数进行修正优化，得到新的汽轮机阀门配汽函数。
2.根据权利要求1所述的一种基于小样本暂态数据的汽轮机高调门流量特性优化方法，其中，步骤一中的采样周期为1秒，所采集参数包括：机组功率、主蒸汽流量、主蒸汽压力、主蒸汽温度、调节级压力、高压缸排汽压力、高压缸排汽温度、总阀位开度、GV1阀位开度、GV2阀位开度、GV3阀位开度、GV4阀位开度。
3.根据权利要求1所述的一种基于小样本暂态数据的汽轮机高调门流量特性优化方法，其中，步骤二所述孤立森林算法，由iTree组成，采用多次迭代方式构建二叉搜索树，并将这些二叉树组成森林，其构建过程包括：
(1)从n个训练数据中随机抽取m个子样本，放入一颗孤立树的根节点，其中；
(2)随机指定一个维度，在当前节点数据范围内，随机产生一个切割点p——切割点产生于当前节点数据中指定维度的最大值与最小值之间；
(3)此切割点的选取生成了一个超平面，将当前节点数据空间切分为2个子空间：把当前所选维度下小于p的点放在当前节点的左分支，把大于等于p的点放在当前节点的右分支；
(4)把节点的左分支和右分支节点递归步骤2、3，不断构造新的叶子节点，直到叶子节点上只有一个数据(无法再继续切割)或树已经生长到了所设定的高度；
(5)记录每个数据从根节点到叶子节点经过边的数量，即路径长度h(x)，以此来计算该数据的异常指数S(x,n):
Figure FDA0002278691710000031
c(n)＝2H(n-1)-(2(n-1)/n)
其中，H(n)＝ln(n)+ξ,ξ＝0.5772156649为欧拉常数，S(x,n)越接近1表示该点是异常点的可能性越高，越接近0表示该点为正常点的可能性越高。
4.根据权利要求1所述的一种基于小样本暂态数据的汽轮机高调门流量特性优化方法，其中，步骤三中选取50％～100％连续变负荷数据段，以1MW间隔作为质心划分依据，确定聚类质心个数。
5.根据权利要求1所述的一种基于小样本暂态数据的汽轮机高调门流量特性优化方法，其中，步骤四中所述的将目标数据集中汽轮机进汽量修正至参考值下对应的汽轮机进汽量所采用的计算方法为：
Figure FDA0002278691710000032
其中，ReaValue代表将目标数据集中汽轮机进汽量修正至参考值下所对应的汽轮机进汽量；
G、p分别表示汽轮机进汽量t/h、压力MPa；GA表示参考工况进汽量；GB表示实际工况进汽量；下标G0表示主蒸汽参数；P1B表示实际工况下调节级蒸汽参数，P0A表示表示参考工况下主蒸汽压力，P0B表示表示实际工况下主蒸汽压力，P1A表示表示参考工况下调节级蒸汽压力，P1B表示表示实际工况下调节级蒸汽压力。
CN201911132428.0A 2019-11-19 2019-11-19 一种基于小样本暂态数据的汽轮机高调门流量特性优化方法 Pending CN110821576A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911132428.0A CN110821576A (zh) 2019-11-19 2019-11-19 一种基于小样本暂态数据的汽轮机高调门流量特性优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911132428.0A CN110821576A (zh) 2019-11-19 2019-11-19 一种基于小样本暂态数据的汽轮机高调门流量特性优化方法

Publications (1)

Publication Number Publication Date
CN110821576A true CN110821576A (zh) 2020-02-21

Family

ID=69556539

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911132428.0A Pending CN110821576A (zh) 2019-11-19 2019-11-19 一种基于小样本暂态数据的汽轮机高调门流量特性优化方法

Country Status (1)

Country Link
CN (1) CN110821576A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113359450A (zh) * 2021-06-08 2021-09-07 国网湖南省电力有限公司 一种阀门流量特性曲线拟合方法及系统
CN113446073A (zh) * 2021-08-11 2021-09-28 内蒙古京隆发电有限责任公司 一种汽轮机高压调节汽门的故障诊断方法
CN114593117A (zh) * 2021-11-02 2022-06-07 杭州远视智能科技有限公司 基于举升液压缸压力测量的叉车载荷重量测量系统及方法
CN115031956A (zh) * 2022-06-15 2022-09-09 河北工业大学 一种节流调节型汽轮机调速汽门特性图谱绘制方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6550253B2 (en) * 2001-09-12 2003-04-22 General Electric Company Apparatus and methods for controlling flow in turbomachinery
JP2004005250A (ja) * 2002-05-31 2004-01-08 Toshiba Corp プロセス最適化制御システム
CN105275508A (zh) * 2015-11-06 2016-01-27 国网河南省电力公司电力科学研究院 一种基于功率值计算的汽轮机流量曲线辨识及优化方法
CN105317476A (zh) * 2015-11-06 2016-02-10 国网河南省电力公司电力科学研究院 一种基于特征通流面积的汽轮机流量曲线辨识及优化方法
CN107035430A (zh) * 2017-05-26 2017-08-11 国家电网公司 一种基于大数据分析的汽轮机阀门流量特性参数优化方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6550253B2 (en) * 2001-09-12 2003-04-22 General Electric Company Apparatus and methods for controlling flow in turbomachinery
JP2004005250A (ja) * 2002-05-31 2004-01-08 Toshiba Corp プロセス最適化制御システム
CN105275508A (zh) * 2015-11-06 2016-01-27 国网河南省电力公司电力科学研究院 一种基于功率值计算的汽轮机流量曲线辨识及优化方法
CN105317476A (zh) * 2015-11-06 2016-02-10 国网河南省电力公司电力科学研究院 一种基于特征通流面积的汽轮机流量曲线辨识及优化方法
CN107035430A (zh) * 2017-05-26 2017-08-11 国家电网公司 一种基于大数据分析的汽轮机阀门流量特性参数优化方法

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113359450A (zh) * 2021-06-08 2021-09-07 国网湖南省电力有限公司 一种阀门流量特性曲线拟合方法及系统
CN113359450B (zh) * 2021-06-08 2022-06-28 国网湖南省电力有限公司 一种阀门流量特性曲线拟合方法及系统
CN113446073A (zh) * 2021-08-11 2021-09-28 内蒙古京隆发电有限责任公司 一种汽轮机高压调节汽门的故障诊断方法
CN114593117A (zh) * 2021-11-02 2022-06-07 杭州远视智能科技有限公司 基于举升液压缸压力测量的叉车载荷重量测量系统及方法
CN114593117B (zh) * 2021-11-02 2023-09-08 杭州远视智能科技有限公司 基于举升液压缸压力测量的叉车载荷重量测量系统及方法
CN115031956A (zh) * 2022-06-15 2022-09-09 河北工业大学 一种节流调节型汽轮机调速汽门特性图谱绘制方法
CN115031956B (zh) * 2022-06-15 2023-02-17 河北工业大学 一种节流调节型汽轮机调速汽门特性图谱绘制方法

Similar Documents

Publication Publication Date Title
CN110821576A (zh) 一种基于小样本暂态数据的汽轮机高调门流量特性优化方法
CN106703904B (zh) 一种基于数据挖掘技术的汽轮机配汽曲线优化方法
CN103850726B (zh) 一种快速确定汽轮机定滑压优化曲线的方法
CN108470249B (zh) 一种耦合聚类分析和决策树的梯级水电站群短期发电调度方法
CN110837226B (zh) 基于智能优化算法的火电机组操作寻优方法及相关装置
CN105275508A (zh) 一种基于功率值计算的汽轮机流量曲线辨识及优化方法
CN105317476A (zh) 一种基于特征通流面积的汽轮机流量曲线辨识及优化方法
WO2004111446A1 (en) Method and installation for extracting energy from a flowing fluid
Saarinen et al. Field measurements and system identification of three frequency controlling hydropower plants
CN110989360A (zh) 一种基于全量数据的火电机组稳态历史寻优方法
CN108363295B (zh) 基于系统辨识的火电机组agc性能考核指标计算与预测方法
CN111159844A (zh) 一种电站燃气轮机排气温度的异常检测方法
CN105240058A (zh) 一种基于喷嘴流量计算的汽轮机流量曲线辨识及优化方法
CN111754030A (zh) 一种基于hac和rf-ga的火电机组供电煤耗优化方法
CN110821575A (zh) 一种基于动态数据的汽轮机配汽曲线优化方法
CN113489024A (zh) 热电联供机组多模式抽汽辅助调峰调频控制系统及方法
CN111931967B (zh) 一种风电场短期功率预测方法
CN111488712B (zh) 一种基于迁移学习的风力发电机功率曲线建模方法
CN113359435A (zh) 用于火电机组动态工况数据的修正方法
CN110457800B (zh) 考虑机械惯性的水平轴风机风速出力折算方法
CN109869198B (zh) 一种基于阀位修正的汽轮机全寿命高压缸效率监测方法
CN107591847B (zh) 一种利用变参数的方式调节水电机组agc的方法
CN113111314B (zh) 一种基于多机组协同的并网风电机组健康状态评估方法
CN114856901A (zh) 一种水轮机调速器导叶开度调节方法及系统
CN113505909B (zh) 一种面向短期风电功率趋势预测的误差补偿方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20200221

RJ01 Rejection of invention patent application after publication