CN115031956B - 一种节流调节型汽轮机调速汽门特性图谱绘制方法 - Google Patents

一种节流调节型汽轮机调速汽门特性图谱绘制方法 Download PDF

Info

Publication number
CN115031956B
CN115031956B CN202210679167.XA CN202210679167A CN115031956B CN 115031956 B CN115031956 B CN 115031956B CN 202210679167 A CN202210679167 A CN 202210679167A CN 115031956 B CN115031956 B CN 115031956B
Authority
CN
China
Prior art keywords
valve
steam flow
array
valve opening
steam
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.)
Active
Application number
CN202210679167.XA
Other languages
English (en)
Other versions
CN115031956A (zh
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.)
Hebei University of Technology
Original Assignee
Hebei University of Technology
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 Hebei University of Technology filed Critical Hebei University of Technology
Priority to CN202210679167.XA priority Critical patent/CN115031956B/zh
Publication of CN115031956A publication Critical patent/CN115031956A/zh
Application granted granted Critical
Publication of CN115031956B publication Critical patent/CN115031956B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00—Testing of machine parts
    • G01M13/003—Machine valves
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N25/00—Investigating or analyzing materials by the use of thermal means
    • G01N25/20—Investigating or analyzing materials by the use of thermal means by investigating the development of heat, i.e. calorimetry, e.g. by measuring specific heat, by measuring thermal conductivity
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/15—Correlation function computation including computation of convolution operations

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Software Systems (AREA)
  • Chemical & Material Sciences (AREA)
  • Algebra (AREA)
  • General Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Databases & Information Systems (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Control Of Turbines (AREA)

Abstract

本发明涉及一种节流调节型汽轮机调速汽门特性图谱绘制方法,包括如下步骤：步骤S1：根据原始数据，按照阀门开度递增顺序，绘制若干条不同蒸汽压力下的阀门开度—蒸汽流量曲线；步骤S2：确定阀门线性段的初始数组D0；步骤S3：确定各阀门开度—蒸汽流量曲线线性段终点数组，形成线性段的终点数组集D1；步骤S4：根据终点数组集D1，确定线性段终点的拟合函数F1(px)；步骤S5：根据初始数组D0和拟合函数F1(px)确定阀门线性段蒸汽流量拟合函数G1(p)。本发明在已有数据基础上，通过阀门线性段和阀门非线性段特性函数的拟合计算，将不同蒸汽压力下的阀门特性(开度‑流量)图谱化，以利于汽轮机精准的进行一次调频或配合其他技术研究。

Description

一种节流调节型汽轮机调速汽门特性图谱绘制方法
技术领域
本发明涉及一种节流调节型汽轮机调速汽门特性图谱绘制方法。
背景技术
随着当前清洁能源发电的高速发展，火力发电在调峰、一次调频等方面的作用更加突显。蒸汽轮机在参与一次调频过程中，调速汽门在低于额定负荷时的动作准确度直接影响机组调频性能。所以，如何准确试验(计算)调速汽门特性图谱是一个重点研究方向。
目前，节流调节型汽轮机的调速汽门试验以热力试验方法为主，此方法受制于试验参数要求高、试验效率低、试验结果精度受多方参数影响等不利因素，试验过程繁琐，需要各主要参数均能够在额定附近运行，现场操作的难度大。
理论上，节流调节型汽轮机可以使用弗留格尔公式拟合调速汽门特性曲线，但由于此类型机组不设置首级压力等重要参数，同时由于补气阀的影响，导致整个试验过程不可控，这也决定了当前仍使用热力学方式测试阀门特性的实际情况。
为保证经济性，汽轮机机组变工况均以滑压运行方式为主，机组从冲转到带满负荷过程中的主蒸汽压力不断变化，加之过程中回热系统与真空系统变化的影响，使得相同蒸汽参数下，汽轮机蒸汽流量与发电负荷呈现带状离散分布关系。因此，热力试验需要的数据采集、数据计算量巨大，过程繁琐。
发明内容
本发明所要解决的技术问题是提供一种精度较高并且计算较简便的节流调节型汽轮机调速汽门特性图谱绘制方法。
本发明采用如下技术方案：
本发明包括如下步骤：
步骤S1：根据原始数据，按照阀门开度递增顺序，绘制若干条不同蒸汽压力下的阀门开度—蒸汽流量曲线；
步骤S2：确定阀门线性段的初始数组D0；
步骤S3：确定各阀门开度—蒸汽流量曲线线性段终点数组，形成线性段的终点数组集D1；
步骤S4：根据终点数组集D1，确定线性段终点的拟合函数F1(px)；其中，px为自变量，代表蒸汽压力；
步骤S5：根据初始数组D0和拟合函数F1(px)确定阀门线性段蒸汽流量拟合函数G1(x)；
步骤S6：在其中一条阀门开度—蒸汽流量曲线上选取阀门开度不小于其线性段终点阀门开度的五个数据组，并且确定各条阀门开度—蒸汽流量曲线的阀门非线性段蒸汽流量拟合函数G2(x)；
步骤S7：根据各条阀门开度—蒸汽流量曲线的阀门非线性段蒸汽流量拟合函数G2
(x)确定各条阀门开度—蒸汽流量曲线的非线性段终点数组，形成非线性段的终点数组集D2；
步骤S8：根据终点数组集D2，确定非线性段终点的拟合函数F2(px)；
步骤S9：按照终点数组集D1到终点数组集D2，将各条阀门开度—蒸汽流量曲线的阀门非线性段蒸汽流量拟合函数G2(x)分为等分的五段，并且根据每段的节点形成Gn、Gn+1、Gn+2以及Gn+3四个数组集；
步骤S10：根据步骤S7和步骤S8，并且根据Gn、Gn+1、Gn+2以及Gn+3四个数组集确定四条拟合函数；
步骤S11：输入需求的蒸汽压力p1，根据拟合函数F1(px)，得到蒸汽压力p1时的线性段终点数组，根据非线性段终点的拟合函数F2(px)，得到得到蒸汽压力p1时的非线性段终点数组，根据对Gn、Gn+1、Gn+2以及Gn+3四个数组集确定的四条拟合函数，确定阀门非线性段蒸汽流量拟合函数G2(x)部分的五等分点，进而利用计算机程序绘制蒸汽压力p1下的阀门开度—蒸汽流量曲线。
本发明所述步骤S2中，初始数组D0为阀门开度逐渐增大，对应的蒸汽流量由0跃迁至某一值，此时的阀门开度和蒸汽流量组成的数组。
本发明所述步骤S3中，对其中一条阀门开度—蒸汽流量曲线中原始数据所在的点进行求导，当相邻两点中存在后一点导数值比前一点导数值大3％时，则将前一点对应的数据组设置为阀门线性段的终点数组。
本发明所述拟合函数F1(px)、拟合函数F2(px)以及根据Gn、Gn+1、Gn+2以及Gn+3四个数组集确定的四条拟合函数均为参数方程，参数为蒸汽压力px，因变量分别为阀门开度和蒸汽流量；
所述阀门线性段蒸汽流量拟合函数G1(x)为线性方程，自变量为阀门开度，因变量为蒸汽流量。
本发明所述阀门非线性段蒸汽流量拟合函数G2(x)为四阶多项式函数，即G2(x)＝a0+a1x+a2x2+a3x3+a4x4，其中x为阀门开度，根据步骤S6的五个数据组进行拟合，当阀门非线性段蒸汽流量拟合函数G2(x)与原始数据之间的拟合优度R2≥0.995时，判定拟合成功。
本发明各条阀门开度—蒸汽流量曲线最大值的97.3％设置为非线性段终点的蒸汽流量值，进而确定各阀门开度—蒸汽流量曲线非线性段终点的阀门开度，进而确定非线性段终点数组和终点数组集D2。
本发明所述步骤S11中，根据零点、初始数组D0、线性段终点数组、阀门非线性段蒸汽流量拟合函数G2(x)部分的五等分点以及非线性段终点数组，绘制蒸汽压力p1时的阀门开度—蒸汽流量曲线。
本发明积极效果如下：
本发明采用全新的方法，以电厂长期运行数据为原料，利用各项参数的相互耦合特性，进行机组稳定运行数据的清洗与筛选，形成排除系统误差的阀门开度、蒸汽流量、机组负荷等特征数据。在已有数据基础上，通过阀门线性段和阀门非线性段特性函数的拟合计算，将不同蒸汽压力下的阀门特性(开度-流量)图谱化，以利于汽轮机精准的进行一次调频或配合其他技术研究。
附图说明
图1为本发明阀门结构示意图；
图2为本发明原理图；
图3为本发明阀门非线性段蒸汽流量拟合函数G2(x)确定流程图；
图4为本发明阀门开度—蒸汽流量曲线示意图。
具体实施方式
下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。以下对至少一个示例性实施例的描述实际上仅仅是说明性的，决不作为对本申请及其应用或使用的任何限制。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本申请的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
除非另外具体说明，否则在这些实施例中阐述的部件和步骤的相对布置、数字表达式和数值不限制本申请的范围。同时，应当明白，为了便于描述，附图中所示出的各个部分的尺寸并不是按照实际的比例关系绘制的。对于相关领域普通技术人员已知的技术、方法和设备可能不作详细讨论，但在适当情况下，所述技术、方法和设备应当被视为授权说明书的一部分。在这里示出和讨论的所有示例中，任何具体值应被解释为仅仅是示例性的，而不是作为限制。因此，示例性实施例的其它示例可以具有不同的值。应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步讨论。
在本申请的描述中，需要理解的是，方位词如“前、后、上、下、左、右”、“横向、竖向、垂直、水平”和“顶、底”等所指示的方位或位置关系通常是基于附图所示的方位或位置关系，仅是为了便于描述本申请和简化描述，在未作相反说明的情况下，这些方位词并不指示和暗示所指的装置或元件必须具有特定的方位或者以特定的方位构造和操作，因此不能理解为对本申请保护范围的限制；方位词“内、外”是指相对于各部件本身的轮廓的内外。
为了便于描述，在这里可以使用空间相对术语，如“在……之上”、“在……上方”、“在……上表面”、“上面的”等，用来描述如在图中所示的一个器件或特征与其他器件或特征的空间位置关系。应当理解的是，空间相对术语旨在包含除了器件在图中所描述的方位之外的在使用或操作中的不同方位。例如，如果附图中的器件被倒置，则描述为“在其他器件或构造上方”或“在其他器件或构造之上”的器件之后将被定位为“在其他器件或构造下方”或“在其他器件或构造之下”。因而，示例性术语“在……上方”可以包括“在……上方”和“在……下方”两种方位。该器件也可以其他不同方式定位(旋转90度或处于其他方位)，并且对这里所使用的空间相对描述作出相应解释。
此外，需要说明的是，使用“第一”、“第二”等词语来限定零部件，仅仅是为了便于对相应零部件进行区别，如没有另行声明，上述词语并没有特殊含义，因此不能理解为对本申请保护范围的限制。
通过研究多台不同时期、不同厂商的节流调节方式汽轮机的阀门物理结构，可以发现明显的特征：曲线由两段或三段形成，即调节死区、线性调节区以及非线性调节区。调节死区为一个较小阀门开度，低于此开度是阀门截止状态；阀门开度超过死区后，进入线性调节区域，蒸汽流量变化与阀门开度为一次线性方程关系式；超越线性调节区域(拐点)后，阀门调节呈非线性关系，蒸汽流量变化与阀门开度为4次多项式关系。所述的汽轮机阀门结构如附图1，蒸汽从左侧进气口1进入，再从下端排气口2排出，阀杆3带动阀塞4上下运动，阀塞4的底面为向下弯曲的弧形面。阀塞4处于初始状态时与排气口2上端的弧形边沿5顶紧，阀门开度为0。当阀门开度增大但处于预启状态时，阀塞4与排气口2上端的弧形边沿5保持接触，此时处于调节死区；当阀门开度增大时，阀门有效通流面积为阀塞4底面的弧形面与排气口2上端弧形边沿5之间的铅垂线组成的面，此时处于线性调节区；当阀门开度继续增大时，阀门有效通流面积为阀塞4底面的弧形面与弧形边沿5外围处之间的铅垂线组成的面，此时为非线性调节区。
本发明数组均包括阀门开度和蒸汽流量两个量。
实施例1
如附图1—3所示，本发明包括如下步骤：
步骤S1：根据原始数据，按照阀门开度递增顺序，绘制若干条不同蒸汽压力下的阀门开度—蒸汽流量曲线；
步骤S2：确定阀门线性段的初始数组D0；
步骤S3：确定各阀门开度—蒸汽流量曲线线性段终点数组，形成线性段的终点数组集D1；
步骤S4：根据终点数组集D1，确定线性段终点的拟合函数F1(px)；
步骤S5：根据初始数组D0和拟合函数F1(px)确定阀门线性段蒸汽流量拟合函数G1(x)；
步骤S6：在其中一条阀门开度—蒸汽流量曲线上选取阀门开度不小于其线性段终点阀门开度的五个数据组，并且确定各条阀门开度—蒸汽流量曲线的阀门非线性段蒸汽流量拟合函数G2(x)；
步骤S7：根据各条阀门开度—蒸汽流量曲线的阀门非线性段蒸汽流量拟合函数G2(x)确定各条阀门开度—蒸汽流量曲线的非线性段终点数组，形成非线性段的终点数组集D2；
步骤S8：根据终点数组集D2，确定非线性段终点的拟合函数F2(px)；
步骤S9：按照终点数组集D1到终点数组集D2，将各条阀门开度—蒸汽流量曲线的阀门非线性段蒸汽流量拟合函数G2(x)分为等分的五段，并且根据每段的节点形成Gn、Gn+1、Gn+2以及Gn+3四个数组集；
步骤S10：根据步骤S7和步骤S8，并且根据Gn、Gn+1、Gn+2以及Gn+3四个数组集确定四条拟合函数；
步骤S11：输入需求的蒸汽压力p1，根据拟合函数F1(px)，得到蒸汽压力p1时的线性段终点数组，根据非线性段终点的拟合函数F2(px)，得到得到蒸汽压力p1时的非线性段终点数组，根据对Gn、Gn+1、Gn+2以及Gn+3四个数组集确定的四条拟合函数，确定阀门非线性段蒸汽流量拟合函数G2(x)部分的五等分点，进而利用计算机程序绘制蒸汽压力p1下的阀门开度—蒸汽流量曲线。
本发明所述步骤S2中，初始数组D0为阀门开度逐渐增大，对应的蒸汽流量由0跃迁至某一值，此时的阀门开度和蒸汽流量组成的数组。
本发明所述步骤S3中，对其中一条阀门开度—蒸汽流量曲线中原始数据所在的点进行求导，当相邻两点中存在后一点导数值比前一点导数值大3％时，则将前一点对应的数据组设置为阀门线性段的终点数组。
本发明所述拟合函数F1(px)、拟合函数F2(px)以及根据Gn、Gn+1、Gn+2以及Gn+3四个数组集确定的四条拟合函数均为参数方程，参数为蒸汽压力px，因变量分别为阀门开度和蒸汽流量；
所述阀门线性段蒸汽流量拟合函数G1(x)为线性方程，自变量为阀门开度，因变量为蒸汽流量。
本发明所述阀门非线性段蒸汽流量拟合函数G2(x)为四阶多项式函数，即G2(x)＝a0+a1x+a2x2+a3x3+a4x4，其中x为阀门开度，根据步骤S6的五个数据组进行拟合，当阀门非线性段蒸汽流量拟合函数G2(x)与原始数据之间的拟合优度R2≥0.995时，判定拟合成功。
实施例2
如附图1—3所示，基于实施例1，本发明各条阀门开度—蒸汽流量曲线最大值的97.3％设置为非线性段终点的蒸汽流量值，进而确定各阀门开度—蒸汽流量曲线非线性段终点的阀门开度，进而确定非线性段终点数组和终点数组集D2。
实施例3
如附图1—3所示，基于实施例1和实施例2，本发明所述步骤S11中，根据零点、初始数组D0、线性段终点数组、阀门非线性段蒸汽流量拟合函数G2(x)部分的五等分点以及非线性段终点数组，绘制蒸汽压力p1时的阀门开度—蒸汽流量曲线。
图4具体为蒸汽压力为额定压力和冲转压力时的阀门开度—蒸汽流量曲线，以及蒸汽压力p1时的阀门开度—蒸汽流量曲线。
本发明提供的计算思路具有独创性，对原始数据格式或数据结构未进行明确要求，应用本发明计算调速汽门特性图谱的各种常规数据处理方法均在本发明保护范围内。
本实例以X轴为阀门开度、Y轴为蒸汽流量、曲线代表某个压力值形式进行计算；由于工程实际情况，更多的阀门特性曲线是X轴代表蒸汽流量，Y轴代表阀门开度；相应的拐点判据也应同步调整。
如何合理选取公式拟合方法均在本发明技术路线的保护范围。
改变拟合函数的方次数仅是实际情况计算与数据量需要，均在本发明保护范围内。
利用矩阵反算多项式系数以及通过快速傅里叶转换并结合坐标变换求取最终多项式函数等其他方法。所以，现有成熟的多项式系数求解方法均属于本发明计算过程的保护范围。
本发明提及的数组与函数，均以(阀门开度x--蒸汽流量G)形式成对出现，在对应求解目标压力p时的参数方程时，采用了x、G分别计算的方法；类似的采用复变量(含三角函数、极坐标方法)函数结合积分变换(图形映射方法)同样可以达到相应的目的，均在本发明保护范围。
目前，本申请的技术方案已经进行了中试，即产品在大规模量产前的小规模实验；中试完成后，在小范围内开展了用户使用调研，调研结果表明用户满意度较高；现在已经着手准备产品正式投产进行产业化(包括知识产权风险预警调研)。
最后应说明的是：以上所述仅为本发明的优选实施例而已，并不用于限制本发明，尽管参照前述实施例对本发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (6)

1.一种节流调节型汽轮机调速汽门特性图谱绘制方法，其特征在于：包括如下步骤：
步骤S1：根据原始数据，按照阀门开度递增顺序，绘制若干条不同蒸汽压力下的阀门开度—蒸汽流量曲线；
步骤S2：确定阀门线性段的初始数组D0；
步骤S3：确定各阀门开度—蒸汽流量曲线线性段终点数组，形成线性段的终点数组集D1；
步骤S4：根据终点数组集D1，确定线性段终点的拟合函数F1（px）；
步骤S5：根据初始数组D0和拟合函数F1（px）确定阀门线性段蒸汽流量拟合函数G1（x）；
步骤S6：在其中一条阀门开度—蒸汽流量曲线上选取阀门开度不小于其线性段终点阀门开度的五个数据组，并且确定各条阀门开度—蒸汽流量曲线的阀门非线性段蒸汽流量拟合函数G2（x）；
步骤S7：根据各条阀门开度—蒸汽流量曲线的阀门非线性段蒸汽流量拟合函数G2（x）确定各条阀门开度—蒸汽流量曲线的非线性段终点数组，形成非线性段的终点数组集D2；
步骤S8：根据终点数组集D2，确定非线性段终点的拟合函数F2（px）；
步骤S9：按照终点数组集D1到终点数组集D2，将各条阀门开度—蒸汽流量曲线的阀门非线性段蒸汽流量拟合函数G2（x）分为等分的五段，并且根据每段的节点形成Gn 、Gn+1、Gn+2以及Gn+3四个数组集；
步骤S10：根据步骤S7和步骤S8，并且根据Gn 、Gn+1、Gn+2 以及Gn+3四个数组集确定四条拟合函数；
步骤S11：输入需求的蒸汽压力p1 ，根据拟合函数F1（px），得到蒸汽压力p1时的线性段终点数组，根据非线性段终点的拟合函数F2（px），得到蒸汽压力p1时的非线性段终点数组，根据对Gn 、Gn+1、Gn+2 以及Gn+3四个数组集确定的四条拟合函数，确定阀门非线性段蒸汽流量拟合函数G2（x）部分的五等分点，进而利用计算机程序绘制蒸汽压力p1下的阀门开度—蒸汽流量曲线；
所述拟合函数F1（px）、拟合函数F2（px）以及根据Gn 、Gn+1、Gn+2 以及Gn+3四个数组集确定的四条拟合函数均为参数方程，参数为蒸汽压力px，因变量为阀门开度和蒸汽流量；
所述阀门线性段蒸汽流量拟合函数G1（x）为线性方程，自变量为阀门开度，因变量为蒸汽流量。
2.根据权利要求1所述的一种节流调节型汽轮机调速汽门特性图谱绘制方法，其特征在于：所述步骤S2中，初始数组D0为阀门开度逐渐增大，对应的蒸汽流量由0跃迁至某一值，此时的阀门开度和蒸汽流量组成的数组。
3.根据权利要求2所述的一种节流调节型汽轮机调速汽门特性图谱绘制方法，其特征在于：所述步骤S3中，对其中一条阀门开度—蒸汽流量曲线中原始数据所在的点进行求导，当相邻两点中存在后一点导数值比前一点导数值大3%时，则将前一点对应的数据组设置为阀门线性段的终点数组。
4.根据权利要求3所述的一种节流调节型汽轮机调速汽门特性图谱绘制方法，其特征在于：所述阀门非线性段蒸汽流量拟合函数G2（x）为四阶多项式函数，即G2（x）=
Figure DEST_PATH_IMAGE001
，其中x为阀门开度，根据步骤S6的五个数据组进行拟合，当阀门非线性段蒸汽流量拟合函数G2（x）与原始数据之间的拟合优度R2≥0.995时，判定拟合成功。
5.根据权利要求4所述的一种节流调节型汽轮机调速汽门特性图谱绘制方法，其特征在于：各条阀门开度—蒸汽流量曲线最大值的97.3%设置为非线性段终点的蒸汽流量值，进而确定各阀门开度—蒸汽流量曲线非线性段终点的阀门开度，进而确定非线性段终点数组和终点数组集D2。
6.根据权利要求5所述的一种节流调节型汽轮机调速汽门特性图谱绘制方法，其特征在于：所述步骤S11中，根据零点、初始数组D0、线性段终点数组、阀门非线性段蒸汽流量拟合函数G2（x）部分的五等分点以及非线性段终点数组，绘制蒸汽压力p1时的阀门开度—蒸汽流量曲线。
CN202210679167.XA 2022-06-15 2022-06-15 一种节流调节型汽轮机调速汽门特性图谱绘制方法 Active CN115031956B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210679167.XA CN115031956B (zh) 2022-06-15 2022-06-15 一种节流调节型汽轮机调速汽门特性图谱绘制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210679167.XA CN115031956B (zh) 2022-06-15 2022-06-15 一种节流调节型汽轮机调速汽门特性图谱绘制方法

Publications (2)

Publication Number Publication Date
CN115031956A CN115031956A (zh) 2022-09-09
CN115031956B true CN115031956B (zh) 2023-02-17

Family

ID=83124715

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210679167.XA Active CN115031956B (zh) 2022-06-15 2022-06-15 一种节流调节型汽轮机调速汽门特性图谱绘制方法

Country Status (1)

Country Link
CN (1) CN115031956B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0004415A1 (en) * 1978-03-24 1979-10-03 Westinghouse Electric Corporation System for minimizing valve throttling losses in a steam turbine power plant
CN103743560A (zh) * 2014-01-14 2014-04-23 国家电网公司 汽轮机deh系统高压调节阀门流量特性测试及整定方法
CN104200121A (zh) * 2014-09-19 2014-12-10 国家电网公司 基于机组运行数据的阀门流量函数优化方法
CN106932201A (zh) * 2017-05-02 2017-07-07 哈尔滨工业大学 考虑阀门节流损失的火电机组纯凝工况滑压曲线确定方法
CN110821576A (zh) * 2019-11-19 2020-02-21 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 一种基于小样本暂态数据的汽轮机高调门流量特性优化方法
CN111505943A (zh) * 2020-06-03 2020-08-07 国电科学技术研究院有限公司 基于全行程建模的汽轮机流量特性优化方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040073400A1 (en) * 2002-10-11 2004-04-15 Mitsubishi Heavy Industries, Ltd. Method for managing lifespans of high temperature gas turbine components and its implementation through computer program
US20090178468A1 (en) * 2008-01-10 2009-07-16 General Electric Company Systems and methods for determining steam turbine operating efficiency
US10387775B2 (en) * 2015-09-09 2019-08-20 Emerson Process Management Power & Water Solutions, Inc. Model-based characterization of pressure/load relationship for power plant load control
CN113962140A (zh) * 2021-09-14 2022-01-21 国网河北省电力有限公司电力科学研究院 基于ga-lstm的汽轮机阀门流量特性函数的优化方法
CN114491378A (zh) * 2021-12-22 2022-05-13 华电电力科学研究院有限公司 一种基于迭代方法的汽机背压计算方法及电子装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0004415A1 (en) * 1978-03-24 1979-10-03 Westinghouse Electric Corporation System for minimizing valve throttling losses in a steam turbine power plant
CN103743560A (zh) * 2014-01-14 2014-04-23 国家电网公司 汽轮机deh系统高压调节阀门流量特性测试及整定方法
CN104200121A (zh) * 2014-09-19 2014-12-10 国家电网公司 基于机组运行数据的阀门流量函数优化方法
CN106932201A (zh) * 2017-05-02 2017-07-07 哈尔滨工业大学 考虑阀门节流损失的火电机组纯凝工况滑压曲线确定方法
CN110821576A (zh) * 2019-11-19 2020-02-21 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 一种基于小样本暂态数据的汽轮机高调门流量特性优化方法
CN111505943A (zh) * 2020-06-03 2020-08-07 国电科学技术研究院有限公司 基于全行程建模的汽轮机流量特性优化方法

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
P. K. Bhowmik ; et al.Flow characteristics analysis of ultra supercritical overload turbine control valve for power plant.《The 2nd IEEE Conference on Power Engineering and Renewable Energy (ICPERE) 2014》.2015, *
Yanna Xi ; et al.Improvement of an Elman Neural Network Model for Once-through Boiler Unit and Its Verification.《2016 12th World Congress on Intelligent Control and Automation (WCICA)》.2014, *
基于数据拟合的汽轮机调节阀流量特性优化;尚星宇等;《热力发电》;20170331;第46卷(第3期);第121-131页 *
基于数据挖掘的汽轮机主蒸汽调节阀特性建模研究;曲广浩等;《计算技术与自动化》;20200328;第39卷(第01期);第43-47页 *
汽轮发电机组流量特性试验分析与研究;张曦等;《控制工程》;20130520;第20卷(第03期);第461-465页 *
汽轮机阀门流量特性分析与优化;张曦等;《南方电网技术》;20101130;第4卷;第72-75页 *
汽轮机阀门流量特性分析与优化研究;张曦等;《全国电站辅机及汽轮机热力系统节能降耗技术论坛论文集》;20110401;第233-237页 *
汽轮机阀门流量特性曲线优化;华兴鲁等;《河北电力技术》;20200425(第02期);第51-53页 *

Also Published As

Publication number Publication date
CN115031956A (zh) 2022-09-09

Similar Documents

Publication Publication Date Title
CN105631140B (zh) 一种变循环发动机稳态性能分析及优化方法
CN106096223B (zh) 一种五孔探针数据处理方法
CN111859746B (zh) 一种基于流场重构的透平机械变工况性能预测方法
CN104533821A (zh) 一种压气机控制方法
CN108448585B (zh) 一种基于数据驱动的电网潮流方程线性化求解方法
CN108107723A (zh) 非线性间歇过程的2d最优模糊控制器设计方法
CN108932392A (zh) 基于改进三重互易边界元法的瞬态温度计算方法
CN110442911B (zh) 一种基于统计机器学习的高维复杂系统不确定性分析方法
CN110219736A (zh) 基于非线性模型预测控制的航空发动机直接推力控制方法
CN108649562A (zh) 一种基于苍狼算法的电力系统负荷建模参数辨识方法
WO2021237910A1 (zh) 一种基于查表法的主动磁轴承控制器的构造方法
CN111783343A (zh) 非结构网格方向模板选择方法与流场求解方法
Merrikh-Bayat et al. Extension of the root-locus method to a certain class of fractional-order systems
CN110489877A (zh) 一种适用于航空发动机实时模型的插值方法
CN107146241A (zh) 一种基于差分进化算法和TrimmedICP算法的点云配准方法
CN115031956B (zh) 一种节流调节型汽轮机调速汽门特性图谱绘制方法
CN105952499B (zh) 一种基于蚁群算法获取汽轮机高压调门组流量的方法
CN108181617A (zh) 一种基于张量积模型变换的非线性调频系统的滤波方法
CN107301268A (zh) 一种船舶燃气轮机可调静叶压气机偏转角度优化方法
CN107203364A (zh) 一种用于压气机全工况特性的预测和辨识方法
CN103742276B (zh) 割草机汽油发动机空燃比跟踪控制方法及其装置
CN108123434B (zh) 一种计算pv曲线斜率以求取pv曲线运行点的方法
CN110137967B (zh) 一种针对关键节点的大规模电力系统潮流收敛性调整方法
CN108919642B (zh) 一种炉跟机协调控制系统控制器参数优化整定方法
CN111291968A (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
GR01 Patent grant
GR01 Patent grant