CN113806346A - 一种基于大数据分析的汽轮机劣化趋势测量方法及终端机 - Google Patents

一种基于大数据分析的汽轮机劣化趋势测量方法及终端机 Download PDF

Info

Publication number
CN113806346A
CN113806346A CN202111089301.2A CN202111089301A CN113806346A CN 113806346 A CN113806346 A CN 113806346A CN 202111089301 A CN202111089301 A CN 202111089301A CN 113806346 A CN113806346 A CN 113806346A
Authority
CN
China
Prior art keywords
data
data set
method based
value
state
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
Application number
CN202111089301.2A
Other languages
English (en)
Other versions
CN113806346B (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.)
Zhejiang Zheneng Taizhou No2 Power Generation Co ltd
Shandong Luruan Digital Technology Co Ltd
Original Assignee
Zhejiang Zheneng Taizhou No2 Power Generation Co ltd
Shandong Luruan Digital Technology 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 Zhejiang Zheneng Taizhou No2 Power Generation Co ltd, Shandong Luruan Digital Technology Co Ltd filed Critical Zhejiang Zheneng Taizhou No2 Power Generation Co ltd
Publication of CN113806346A publication Critical patent/CN113806346A/zh
Application granted granted Critical
Publication of CN113806346B publication Critical patent/CN113806346B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/21Design, administration or maintenance of databases
    • G06F16/215Improving data quality; Data cleansing, e.g. de-duplication, removing invalid entries or correcting typographical errors
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F01MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
    • F01DNON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
    • F01D21/00Shutting-down of machines or engines, e.g. in emergency; Regulating, controlling, or safety means not otherwise provided for
    • F01D21/003Arrangements for testing or measuring
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/24Querying
    • G06F16/245Query processing
    • G06F16/2458Special types of queries, e.g. statistical queries, fuzzy queries or distributed queries
    • G06F16/2462Approximate or statistical queries
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Databases & Information Systems (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Probability & Statistics with Applications (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • Fuzzy Systems (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Computational Linguistics (AREA)
  • Evolutionary Computation (AREA)
  • Mechanical Engineering (AREA)
  • Quality & Reliability (AREA)
  • Testing And Monitoring For Control Systems (AREA)

Abstract

本发明提供一种基于大数据分析的汽轮机劣化趋势测量方法及终端机,根据预设的时间段,以预设的时间间隔,从数据库中获取目标设备的相关参数构建各时间段的矩阵。再进行稳态及非稳态判定,过滤其中非稳态数据得到剩余稳态数据集;对稳态数据集通过测点逐条计算高压缸效率;将整合数据集进行分割;对数据进行非线性拟合;将同一层的数据重新组合成分层数据集,对稳定的分层数据集进行非线性拟合,得到阈值矩阵;统计合适层的均值得出劣化值。方法以数据对数据做比较,提升了计算所得劣化数值的真实性与准确性。能够建立相对准确的一次非线性拟合模型,有助于之后循环划分拟合各层非线性模型,划分各层模型间缸效分界阈值。

Description

一种基于大数据分析的汽轮机劣化趋势测量方法及终端机
技术领域
本发明涉及火电机组汽轮机发电技术领域,具体涉及一种基于大数据分析的汽轮机劣化趋势测量方法及终端机。
背景技术
汽轮机通流部分是将蒸汽热能转换为功的核心部件,其完善程度对机组能耗水平有重要影响。在机组设备投产或大修后,随着机组设备的运行,其性能势必会产生下降,高压缸效率降低,这便是高压缸性能的劣化。掌握其劣化规律,以便判断电厂设备管理水平,进行调整与优化,指导电厂节能工作。
汽轮机厂家给出的缸效率设计值,通常只涉及到几个工况,而且机组运行工况往往会偏离试验工况,随着高压调阀开度和负荷的变化,缸效率也会同步发生变化,使得缸效率设计值或试验值不能反映机组整个运行过程中通流部分的变化,以此方法难以判定每年具体劣化数值。
而现有技术中,使用ARMA模型预测方法,用历史数据预测得到的预测值与实际值做对比得到劣化数值,以此方法得出的劣化数值与预测模型的精确度有较大关系,难以保证其结果准确性。
现有技术中还有利用经验模态分解方法与特征流通面积去计算劣化数值,无法对后期微小劣化算得结果。
发明内容
本发明提供基于大数据分析的汽轮机劣化趋势测量方法,方法利用历史实际运行数值,通过多算法与统计对历史数据组合性使用,得出更为精准且真实的数值,有利于维持汽轮机设备良好运行与持续性维护。
方法包括:
步骤1、根据预设的时间段,以预设的时间间隔,从数据库中获取目标设备的相关参数构建各时间段的矩阵Ft
步骤2、对采集到的各时间段的矩阵Ft进行数据清洗,剔除不合格数据,获得正常数据集Fzt
步骤3、进行稳态及非稳态判定,过滤其中非稳态数据得到剩余稳态数据集;
步骤4、对稳态数据集通过测点逐条计算高压缸效率,并将算得的缸效率值存入稳态数据集中,合并得到整合数据集Ftb_W
步骤5、将整合数据集Ftb_W按照机组负荷和调阀开度进行分割,均匀的划分出i×j个区间;
步骤6、对各个区间下数据集F分别抽取其下高压缸效率最大最小值,计算其最大值与最小值之间的差值与均值的百分比,作为区间下数据的离散度;
步骤7、选取各个区间下离散度MOD值最大的前N个区间,将这N个区间内的数据通过高压缸效率的最大最小值等间距的进行分割出m层,并将处于同一层的数据组合成分层数据集Ftfc,共m个;
步骤8、分别对每一层的数据Ftfc进行非线性拟合;
步骤9、将各个区间内的数据分别代入每一层的非线性拟合模型中,得到其高压缸缸效效预测值,并与实际值做差比较,距离最近的点划分入于该层数据中,并再次将同一层的数据重新组合成分层数据集Ftfc
步骤10、重复步骤8、9直到分层数据集Ftfc保持不变;
步骤11、对稳定的分层数据集Ftfc分别再次进行非线性拟合,再对不同层在各个区间下的高压缸效率进行分界阈值划分,得到阈值矩阵L;
步骤12、对历年数据集Ft分别通过步骤10中的阈值矩阵L进行分层划分,分别计算不同层中历年数据的均值,最终统计合适层的均值得出劣化值。
本发明还提供一种实现基于大数据分析的汽轮机劣化趋势测量方法的终端机,包括:
存储器,用于存储计算机程序及基于大数据分析的汽轮机劣化趋势测量方法;
处理器,用于执行所述计算机程序及基于大数据分析的汽轮机劣化趋势测量方法,以实现基于大数据分析的汽轮机劣化趋势测量方法的步骤。
从以上技术方案可以看出,本发明具有以下优点:
相比较现有技术采用以历史数据做预测方案并与当前数据做差值比较,本发明提供的基于大数据分析的汽轮机劣化趋势测量方法直接历史运行数值,以数据对数据做比较,提升了计算所得劣化数值的真实性与准确性。
在机组实际运行中,缸效率会根据各类运行工况发生动态的变化,通过对同机组负荷与调阀开度下的各个高压缸效率进行分层等同于对影响缸效的其他参数进行了模糊工况划分的效果,避免了不同工况下的高压缸效率做比较的劣化分析误区。
本发明提供的基于大数据分析的汽轮机劣化趋势测量方法中,在通过选取前N个极差除以均值区间数据进行分层,提高了第一次分层的准确率,避免了样本数较少的区间对分层上的影响,能够建立相对准确的一次非线性拟合模型,有助于之后循环划分拟合各层非线性模型,划分各层模型间缸效分界阈值。
本发明提供的基于大数据分析的汽轮机劣化趋势测量方法中,在历史数据上各年份数据在经过稳态非稳态检测后通过划分出的分界阈值进行分层,再对各层、各个区间中的数据统计均值,最终将不同年份在相同层数与区间中的均值做差,剔除缺失数据,按区间样本数统计其均值的均值,可得到相对准确具体劣化数值。
附图说明
为了更清楚地说明本发明的技术方案,下面将对描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为基于大数据分析的汽轮机劣化趋势测量方法流程图;
图2为基于大数据分析的汽轮机劣化趋势测量方法实施例流程图;
图3为实例中缸效率历年变化趋势图。
具体实施方式
本发明提供的基于大数据分析的汽轮机劣化趋势测量方法中,所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、计算机软件或者二者的结合来实现,为了清楚地说明硬件和软件的可互换性,在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本发明的范围。
本发明提供的基于大数据分析的汽轮机劣化趋势测量方法中,提供许多具体细节从而给出对本发明的实施例的充分理解。然而,本领域技术人员将意识到,可以实践本发明的技术方案而没有特定细节中的一个或更多,或者可以采用其它的方法、组元、装置、步骤等。在其它情况下,不详细示出或描述公知方法、装置、实现或者操作以避免模糊本发明的各方面。
如图1和图2所示,本发明提供的基于大数据分析的汽轮机劣化趋势测量方法包括:
S101、根据预设的时间段,以预设的时间间隔,从数据库中获取目标设备的相关参数构建各时间段的矩阵Ft
其中,
构建各时间段的矩阵Ft包括:分别构建各时间段m1×n的Ft矩阵,其中m1为当前年份数据的条数,n为目标设备相关参数测点数,v为个点数值;
Figure BDA0003266679160000041
预设的时间段可以以年,季度,月,周,天等等时间段来设置。预设的时间间隔可以以秒,分,时,天等等时间间隔来设置,具体设置方式根据使用条件来进行设置,具体数据不做限定。
作为获取目标设备的相关参数包括:机组负荷、调阀开度、进汽压力、进汽温度、排汽压力以及排汽温度。
S102、对采集到的各时间段的矩阵Ft进行数据清洗,剔除不合格数据,获得正常数据集Fzt
剔除的不合格数据可以涉及停机、异常、突变等相关的质量不合格数据。
也就是说对采集到的各时间段的矩阵Ft进行数据清洗包括,利用箱线图剔除采集数据中突变、异常、停机的质量不合格数据,得到整合数据集listzh
S103、进行稳态及非稳态判定,过滤其中非稳态数据得到剩余稳态数据集;
具体来讲,参与稳态识别的参数为机组负荷和调阀开度,两者需同时满足条件判断为实时稳态,具体条件可以根据实际运行需要进行设置,这里不做限定。
进一步的讲,
步骤3.1、从矩阵Ft中抽取与机组稳态非稳态判定相关测点N的数据FtN
Figure BDA0003266679160000051
步骤3.2、对数据FtN按照时间窗口长度T1进行2阶多项式滤波,得到滤波后的测点N的数据FtN_L
Figure BDA0003266679160000052
步骤3.3、对滤波后数据FtN_L按照时间窗口长度T2进行一阶线性拟合,得到时间窗口长度T2内各段数据的斜率表Lg
步骤3.4、对测点N的原始数据FtN按照时间窗口长度T2计算其时间窗口长度内的方差表Lv
步骤3.5、通过设置合适的斜率与方差阈值分别从Lg与Lv中过滤都在阈值内的数据序号,按照序号索引Ft,得到当年处于稳态时段下的数据Ft_W
S104、对稳态数据集通过测点逐条计算高压缸效率,并将算得的缸效率值存入稳态数据集中,合并得到整合数据集Ftb_W
这里,将算得的缸效率值v存入数据集中,合并得到整合数据集Ftb_W,其中m2为稳态条件下的数据条数(m2≤m1);
Figure BDA0003266679160000061
对采集的数据计算缸效率,并将算得的缸效率放入该条数据中,其中缸效率计算公式为:
η=(H0-Hk)÷(H0-Hs)
上式中,η为缸内效率,H0表示进汽蒸汽焓,Hk表示排汽焓,Hs表示由进汽蒸汽熵和排汽压力查得等熵焓。
S105、将整合数据集Ftb_W按照机组负荷和调阀开度进行分割,均匀的划分出i×j个区间;
步骤S105中,将整合数据集Ftb_W中的各条数据根据负荷与调阀开度分别划分进各个区间内,构建分割数据集FA,其下数据集F分别代表各个区间下的数据集;
Figure BDA0003266679160000071
将整合数据集Ftb_W按照步骤1中采集的机组负荷和调阀开度进行分割,机组负荷按照a MW一段、调阀开度按照b%一段进行分隔,得到分割数据集。
S106、对各个区间下数据集F分别抽取其下高压缸效率最大最小值,计算其最大值与最小值之间的差值与均值的百分比,作为区间下数据的离散度;
区间下数据的离散度为;
Figure BDA0003266679160000072
S107、选取各个区间下离散度MOD值最大的前N个区间,将这N个区间内的数据通过高压缸效率的最大最小值等间距的进行分割出m层,并将处于同一层的数据组合成分层数据集Ftfc,共m个;
步骤7将离散度最大的N个区间的缸效率均值按照远近进行等距分层,得到分层数据集
Figure BDA0003266679160000073
组合规则如下:
Figure BDA0003266679160000074
上式中,
Figure BDA0003266679160000075
表示在区间i、j中的第m层分层数据集。
S108、分别对每一层的数据Ftfc进行非线性拟合;
在步骤8中对每一层分层数据集按照缸效率、调阀开度、机组负荷三个维度,进行非线性回归拟合,得到每层的回归方程、拟合曲面和均方误差MSE与R2决定系数。
S109、将各个区间内的数据分别代入每一层的非线性拟合模型中,得到其高压缸缸效效预测值,并与实际值做差比较,距离最近的点划分入于该层数据中,并再次将同一层的数据重新组合成分层数据集Ftfc
S110、重复步骤8、9直到分层数据集Ftfc保持不变;
S111、对稳定的分层数据集Ftfc分别再次进行非线性拟合,再对不同层在各个区间下的高压缸效率进行分界阈值划分,得到阈值矩阵L;
S112、对历年数据集Ft分别通过步骤10中的阈值矩阵L进行分层划分,分别计算不同层中历年数据的均值,最终统计合适层的均值得出劣化值。
本发明提出了一种基于大数据分析的汽轮机劣化趋势测量方法。本发明通过采用大数据分析,从多个维度划分合理且准确汽轮机工况,通过对相同工况下的汽轮机历史数据进行差值比较,最终得到具体劣化数值,对其劣化趋势的观察与测量对于汽轮机的安全性和经济性具有重要意义。
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本文通过对南方某电厂的数据进行具体实施来验证本发明的有效性与可靠性,使用从2016年1月1日00:00至2016年12月31日24:00的数据进行其模糊工况对应高压缸效率分界阈值的挖掘,2016年1月1日00:00至2020年12月31日24:00的数据进行汽轮机劣化数值的计算。
步骤1、采集的参数包括:机组负荷、高压调阀开度、主蒸汽压力、主蒸汽温度、高压缸排汽压力、高压缸排汽温度,共计21个测点。取数时间间隔为5秒,获取2016-2020共计31536000条数据,各年数据为6307200条,对每年的数据分别建立一个数据集,共建立5个,分别为F2016,F2017,F2018,F2019,F2020
Figure BDA0003266679160000081
上式中,t代表2016至2020各个年份,m1为数据行数,m1=6307200,n为测点数,n=21。
步骤2、利用箱线图剔除采集数据中突变、异常、停机等质量不合格数据5312600条,再经专业人员筛选后,获得实际运行时正常数据集26223400条数据,其中2016年数据集共存有5171904条,再次保存为Ft
Figure BDA0003266679160000091
上式中,t代表2016至2020各个年份,m1为数据行数,m1=5171904,n为测点数,n=21。
通过计算各参数数据四分位距iqr、上四分位数prctile75及下四分位数,计算其参数上下限阈值为:
threupper=prctile75+3*iqr
threlower=prctile25-3*iqr
每个参数的数据按照上限阈值threupper下限阈值threlower判断并剔除该参数存在的异常数据。
步骤3、以2016年数据的稳态非稳态判定与筛选为例,得到其中稳态时间段内数据集:
a.设定机组负荷与调阀开度作为稳态非稳态参考测点,抽取得到这两个参数的数据集FtN,;
Figure BDA0003266679160000101
其中FtN分别为2016至2020年稳态非稳态判定使用数据集,为例m2为数据行数,以2016年为例,m2=5171904。
b.分别对机组负荷与调阀开度分布进行时间窗口长度为5分钟与1分钟的二阶滤波得到新数据集FtN_L
Figure BDA0003266679160000102
c.随后对机组负荷与调阀开度测点滤波后的数据集FtN_L分别进行时间窗口长度为20分钟与5分钟一阶线性拟合,将一阶线性拟合后的斜率数据Lg作为扩维数据存储作为新列存储回步骤2的Ft中。
Figure BDA0003266679160000103
鉴于取数为5秒一条,故20分钟数据条数为240条,上式中
Figure BDA0003266679160000104
为第m2条数据至第m2-240条数据一阶线性拟合后的斜率数据。
Figure BDA0003266679160000111
d.机组负荷与调阀开度未经滤波前的原始数据进行时间窗口长度为20分钟与5分钟的方差计算,将方差扩维数据Lv作为新列存储在上个步骤中的Ft中。
Figure BDA0003266679160000112
鉴于取数为5秒一条,故20分钟数据条数为240条,上式中
Figure BDA0003266679160000113
为第m2条数据至第m2-240条数据方差数据。
Figure BDA0003266679160000114
e.通过分别设置机组负荷与调阀开度斜率阈值都为0.03,机组负荷与调阀开度方差阈值分别为7.6与1.3进行数据过滤,识别出稳态数据段并存储稳态时间段中的数据,最终从2016年数据集Ft中收集到的稳态数据集Ft_W有1084672条。
Figure BDA0003266679160000121
其中Ft_W分别为2016至2020年稳态数据集,m3为数据行数,以2016年为例,m3=1084672。
步骤4、对处理后的稳态时间段中的数据逐条计算高压缸效率,对计算得到缸效率值大于100与小于80的数据进行初步的筛选与剔除,得到975333条数据,再从原数据集Ft_W抽取高压缸效率、机组负荷与调阀开度保存为新数据集Ftp
Figure BDA0003266679160000122
其中Ft_W分别为2016至2020年稳态数据集,m4为数据行数,以2016年为例,m4=975333。
步骤5、将数据集Ftp通过每10负荷1个区间与每1%调阀开度1个区间,将400-1050负荷数据与22-95调阀开度均匀的以此划分进65*73个区间内。
步骤6、对65*73个区间下的数据集分别计算高压缸效率的最大值与最小值之间的差值除以其均值作为离散度数值,统计成表。
Figure BDA0003266679160000131
步骤7、根据选取表中MOD值大于0.02的数据,筛选得到695个区间内据,共有468623条数据,按照(最大值-最小值)/9,进行等距分为10层,划分分别属于10层数据集下的数据。
步骤8、对每层数据分别进行非线性回归拟合,这里采用设机组负荷为x,调阀开度为y,进行3阶扩维,得到新参数x2,xy,y2,x3,x2y,xy2,y3,将原始参数x,y与新参数再带入常数项后,拟合高压缸效率,最终得到各个参数的权重的方案。
返回每层的回归方程、拟合曲面和均方误差MSE与R2决定系数。拟合效果如下表:
Figure BDA0003266679160000132
Figure BDA0003266679160000141
步骤9、将步骤5所得共计975333条数据再分别代入每一层的非线性拟合模型中,将各条数据在各层拟合模型中返回的预测值与实际值做差,将偏差最小的那一层作为这条数据的归属层,再次划入该层对应的数据集中。
步骤10、重复进行步骤8与步骤9,直到分层数据集中各层的数据不再发生变化。
步骤11、最终得到的各层的回归方程、拟合曲面和均方误差MSE与R2决定系数。均方误差MSE与R2决定系数如下表:
层数 MSE R2 缸效率均值 数据条数
第一层 0.1005 0.9587 83.7626 21286
第二层 0.0402 0.9898 84.0649 39721
第三层 0.0306 0.9916 84.2658 95844
第四层 0.0213 0.9935 84.5884 130109
第五层 0.0360 0.9915 84.9289 95979
第六层 0.0747 0.9856 85.3773 54854
第七层 0.1065 0.9587 85.7087 32684
第八层 0.1218 0.9308 85.9826 16596
第九层 0.2384 0.9370 86.3715 7478
第十层 2.1106 0.8970 86.8847 1083
再对不同层在各个区间下的高压缸效率进行分界阈值划分,结果如下表:
层数 阈值
第一层至第二层 83.9226
第二层至第三层 84.1844
第三层至第四层 84.4129
第四层至第五层 84.7451
第五层至第六层 85.1836
第六层至第七层 85.5932
第七层至第八层 85.8817
第八层至第九层 86.2477
第九层至第十层 865281
步骤12、将2017年至2020年在经过步骤3与步骤4后得到的数据分别用步骤11得出的各层阈值进行分层,并计算每一层数据均值,统计2016年至2020年结果如下表。
Figure BDA0003266679160000151
由于边界工况中可能存在极端数据,故根据结果仅统计第二至九层均值,并进行差值计算得,自16年至2020年逐年劣化数值为0.3215,0.11525,0.047,0.00725,随时间变化其劣化数值逐渐趋于稳定,对其2016年至2020年第二至九层均值作图,结果如图3。
鉴于该电厂机组于2015年投产,根据劣化数值各层计算结果综合显示汽轮机的劣化是处于逐渐降低并稳定的过程,投入运行时间越短,劣化越明显,在投入运行的第4-5年,劣化逐渐可忽略不计,根据现有数据即可建立完整汽轮机劣化模型。
本发明还提供一种实现基于大数据分析的汽轮机劣化趋势测量方法的终端机,包括:存储器,用于存储计算机程序及基于大数据分析的汽轮机劣化趋势测量方法;处理器,用于执行所述计算机程序及基于大数据分析的汽轮机劣化趋势测量方法,以实现基于大数据分析的汽轮机劣化趋势测量方法的步骤。
终端机可以以各种形式来实施。例如,本发明实施例中描述的终端可以包括诸如笔记本电脑、个人数字助理(PDA,Personal Digital Assistant)、平板电脑(PAD)等等的移动终端以及诸如数字TV、台式计算机等等的固定终端。
实现基于大数据分析的汽轮机劣化趋势测量方法的终端机是结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、计算机软件或者二者的结合来实现,为了清楚地说明硬件和软件的可互换性,在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本发明的范围。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (10)

1.一种基于大数据分析的汽轮机劣化趋势测量方法,其特征在于,方法包括:
步骤1、根据预设的时间段,以预设的时间间隔,从数据库中获取目标设备的相关参数构建各时间段的矩阵Ft
步骤2、对采集到的各时间段的矩阵Ft进行数据清洗,剔除不合格数据,获得正常数据集Fzt
步骤3、进行稳态及非稳态判定,过滤其中非稳态数据得到剩余稳态数据集;
步骤4、对稳态数据集通过测点逐条计算高压缸效率,并将算得的缸效率值存入稳态数据集中,合并得到整合数据集Ftb_
步骤5、将整合数据集Ftb_按照机组负荷和调阀开度进行分割,均匀的划分出i×j个区间;
步骤6、对各个区间下数据集F分别抽取其下高压缸效率最大最小值,计算其最大值与最小值之间的差值与均值的百分比,作为区间下数据的离散度;
步骤7、选取各个区间下离散度MOD值最大的前N个区间,将这N个区间内的数据通过高压缸效率的最大最小值等间距的进行分割出m层,并将处于同一层的数据组合成分层数据集Ftfc,共m个;
步骤8、分别对每一层的数据Ftfc进行非线性拟合;
步骤9、将各个区间内的数据分别代入每一层的非线性拟合模型中,得到高压缸缸效效预测值,并与实际值做差比较,距离最近的点划分入于该层数据中,并再次将同一层的数据重新组合成分层数据集Ftfc
步骤10、重复步骤8、9直到分层数据集Ftfc保持不变;
步骤11、对稳定的分层数据集Ftfc分别再次进行非线性拟合,再对不同层在各个区间下的高压缸效率进行分界阈值划分,得到阈值矩阵L;
步骤12、对历年数据集Ft分别通过步骤10中的阈值矩阵L进行分层划分,分别计算不同层中历年数据的均值,最终统计合适层的均值得出劣化值。
2.根据权利要求1所述的基于大数据分析的汽轮机劣化趋势测量方法,其特征在于,步骤1中,构建各时间段的矩阵Ft包括:分别构建各时间段m1×n的Ft矩阵,其中m1为当前年份数据的条数,n为目标设备相关参数测点数,v为个点数值;
Figure FDA0003266679150000021
3.根据权利要求1或2所述的基于大数据分析的汽轮机劣化趋势测量方法,其特征在于,步骤3中,
步骤3.1、从矩阵Ft中抽取与机组稳态非稳态判定相关测点N的数据FtN
Figure FDA0003266679150000022
步骤3.2、对数据FtN按照时间窗口长度T1进行2阶多项式滤波,得到滤波后的测点N的数据FtN_L
Figure FDA0003266679150000023
步骤3.3、对滤波后数据FtN_L按照时间窗口长度T2进行一阶线性拟合,得到时间窗口长度T2内各段数据的斜率表Lg
步骤3.4、对测点N的原始数据FtN按照时间窗口长度T2计算其时间窗口长度内的方差表Lv
步骤3.5、通过设置合适的斜率与方差阈值分别从Lg与Lv中过滤都在阈值内的数据序号,按照序号索引Ft,得到当年处于稳态时段下的数据Ft_W
4.根据权利要求1或2所述的基于大数据分析的汽轮机劣化趋势测量方法,其特征在于,
步骤4中,
将算得的缸效率值v存入数据集中,合并得到整合数据集Ftb_W,其中m2为稳态条件下的数据条数(m2≤m1);
Figure FDA0003266679150000031
对采集的数据计算缸效率,并将算得的缸效率放入该条数据中,其中缸效率计算公式为:
η=(H0-Hk)÷(H0-HS)
上式中,η为缸内效率,H0表示进汽蒸汽焓,Hk表示排汽焓,Hs表示由进汽蒸汽熵和排汽压力查得等熵焓。
5.根据权利要求1或2所述的基于大数据分析的汽轮机劣化趋势测量方法,其特征在于,
步骤5中,
将整合数据集Ftb_W中的各条数据根据负荷与调阀开度分别划分进各个区间内,构建分割数据集FA,其下数据集F分别代表各个区间下的数据集;
Figure FDA0003266679150000041
将整合数据集Ftb_W按照步骤1中采集的机组负荷和调阀开度进行分割,机组负荷按照aMW一段、调阀开度按照b%一段进行分隔,得到分割数据集。
6.根据权利要求1或2所述的基于大数据分析的汽轮机劣化趋势测量方法,其特征在于,
步骤6中,
区间下数据的离散度为;
Figure FDA0003266679150000042
7.根据权利要求1或2所述的基于大数据分析的汽轮机劣化趋势测量方法,其特征在于,
步骤2中,对采集到的各时间段的矩阵Ft进行数据清洗包括,利用箱线图剔除采集数据中突变、异常、停机的质量不合格数据,得到整合数据集listzh
8.根据权利要求1或2所述的基于大数据分析的汽轮机劣化趋势测量方法,其特征在于,
步骤7将离散度最大的N个区间的缸效率均值按照远近进行等距分层,得到分层数据集
Figure FDA0003266679150000043
组合规则如下:
Figure FDA0003266679150000044
上式中,
Figure FDA0003266679150000045
表示在区间i、j中的第m层分层数据集。
9.根据权利要求1或2所述的基于大数据分析的汽轮机劣化趋势测量方法,其特征在于,
在步骤8中对每一层分层数据集按照缸效率、调阀开度、机组负荷三个维度,进行非线性回归拟合,得到每层的回归方程、拟合曲面和均方误差MSE与R2决定系数。
10.一种实现基于大数据分析的汽轮机劣化趋势测量方法的终端机,其特征在于,包括:
存储器,用于存储计算机程序及基于大数据分析的汽轮机劣化趋势测量方法;
处理器,用于执行所述计算机程序及基于大数据分析的汽轮机劣化趋势测量方法,以实现如权利要求1至9任意一项所述基于大数据分析的汽轮机劣化趋势测量方法的步骤。
CN202111089301.2A 2021-08-25 2021-09-16 一种基于大数据分析的汽轮机劣化趋势测量方法及终端机 Active CN113806346B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202110982153 2021-08-25
CN2021109821530 2021-08-25

Publications (2)

Publication Number Publication Date
CN113806346A true CN113806346A (zh) 2021-12-17
CN113806346B CN113806346B (zh) 2024-07-16

Family

ID=78941406

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111089301.2A Active CN113806346B (zh) 2021-08-25 2021-09-16 一种基于大数据分析的汽轮机劣化趋势测量方法及终端机

Country Status (1)

Country Link
CN (1) CN113806346B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114818127A (zh) * 2022-04-18 2022-07-29 中节能工业节能有限公司 一种汽轮机劣化程度的量化方法
CN116050674A (zh) * 2023-03-31 2023-05-02 长江空间信息技术工程有限公司(武汉) 水利工程运行趋势预测方法以及装置
CN116500898A (zh) * 2023-05-11 2023-07-28 华电国际电力股份有限公司莱城发电厂 基于特征流量辨识的火电机组agc负荷控制系统
WO2024067824A1 (zh) * 2022-09-30 2024-04-04 广州明珞装备股份有限公司 气缸劣化趋势确定方法及装置、电子设备、存储介质

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106989377A (zh) * 2015-12-28 2017-07-28 安萨尔多能源公司 用于控制热回收蒸汽产生器中热交换的方法和设备
CN107742053A (zh) * 2017-11-28 2018-02-27 国华(河北)新能源有限公司 风电机组异常识别方法及装置
CN107908865A (zh) * 2017-11-14 2018-04-13 西安热工研究院有限公司 一种基于机器算法的汽轮机低压缸排汽焓实时计算方法
KR20180126648A (ko) * 2017-05-17 2018-11-28 한국전력공사 수명 예측 방법 및 시스템
CN110007652A (zh) * 2019-03-22 2019-07-12 华中科技大学 一种水电机组劣化趋势区间预测方法与系统
CN110094251A (zh) * 2019-05-05 2019-08-06 东南大学 基于分时段多模型建模的scr催化剂性能劣化分析方法
CN111461450A (zh) * 2020-04-10 2020-07-28 河海大学 基于eemd和lstm的水电机组劣化度预测方法
CN111522858A (zh) * 2020-03-16 2020-08-11 国家电网有限公司 一种基于多维状态向量的抽蓄机组性能劣化预警方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106989377A (zh) * 2015-12-28 2017-07-28 安萨尔多能源公司 用于控制热回收蒸汽产生器中热交换的方法和设备
KR20180126648A (ko) * 2017-05-17 2018-11-28 한국전력공사 수명 예측 방법 및 시스템
CN107908865A (zh) * 2017-11-14 2018-04-13 西安热工研究院有限公司 一种基于机器算法的汽轮机低压缸排汽焓实时计算方法
CN107742053A (zh) * 2017-11-28 2018-02-27 国华(河北)新能源有限公司 风电机组异常识别方法及装置
CN110007652A (zh) * 2019-03-22 2019-07-12 华中科技大学 一种水电机组劣化趋势区间预测方法与系统
CN110094251A (zh) * 2019-05-05 2019-08-06 东南大学 基于分时段多模型建模的scr催化剂性能劣化分析方法
CN111522858A (zh) * 2020-03-16 2020-08-11 国家电网有限公司 一种基于多维状态向量的抽蓄机组性能劣化预警方法
CN111461450A (zh) * 2020-04-10 2020-07-28 河海大学 基于eemd和lstm的水电机组劣化度预测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
梁汝印: "《汽轮发电机组智能诊断大数据系统研究及应用》", 当代电力文化, no. 33, 19 April 2021 (2021-04-19), pages 1 - 3 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114818127A (zh) * 2022-04-18 2022-07-29 中节能工业节能有限公司 一种汽轮机劣化程度的量化方法
WO2024067824A1 (zh) * 2022-09-30 2024-04-04 广州明珞装备股份有限公司 气缸劣化趋势确定方法及装置、电子设备、存储介质
CN116050674A (zh) * 2023-03-31 2023-05-02 长江空间信息技术工程有限公司(武汉) 水利工程运行趋势预测方法以及装置
CN116050674B (zh) * 2023-03-31 2023-06-23 长江空间信息技术工程有限公司(武汉) 水利工程运行趋势预测方法以及装置
CN116500898A (zh) * 2023-05-11 2023-07-28 华电国际电力股份有限公司莱城发电厂 基于特征流量辨识的火电机组agc负荷控制系统
CN116500898B (zh) * 2023-05-11 2024-03-19 华电国际电力股份有限公司莱城发电厂 基于特征流量辨识的火电机组agc负荷控制系统

Also Published As

Publication number Publication date
CN113806346B (zh) 2024-07-16

Similar Documents

Publication Publication Date Title
CN113806346B (zh) 一种基于大数据分析的汽轮机劣化趋势测量方法及终端机
CN110442833A (zh) 一种基于多维度scada数据评估风电机组健康状态评估方法
CN106447210A (zh) 一种计及可信度评价的配网设备健康度动态诊断方法
EP2239642B1 (en) Analysis method
CN104794535B (zh) 一种基于主导行业的电力需求预测及预警的方法
CN114397526A (zh) 状态全息感知数据驱动的电力变压器故障预测方法及系统
Oberoi et al. SPC (statistical process control): a quality control technique for confirmation to ability of process
Evans et al. Multivariate data analytics in PV manufacturing—Four case studies using manufacturing datasets
Chen et al. Recipe-independent indicator for tool health diagnosis and predictive maintenance
CN113632025B (zh) 用于评估工业环境中的能量消耗的方法、系统和计算机程序产品
Chang Collaborative decision making algorithm for selection of optimal wire saw in photovoltaic wafer manufacture
CN116910128A (zh) 风电机组硬软故障数据分析系统及其分析方法
CN115358355B (zh) 一种主变压器油温表及顶层油温异常的判断方法和装置
Corsini et al. Assessment of a diagnostic procedure for the monitoring and control of industrial processes
CN115879680A (zh) 钢铁表面缺陷判定规则管理系统
CN111612221B (zh) 一种基于电流指纹技术的工厂工效统计方法
CN111199419B (zh) 股票异常交易的识别方法及系统
Bassetto et al. Operational methods for improving manufacturing control plans: case study in a semiconductor industry
Lin et al. Forecasting rare faults of critical components in LED epitaxy plants using a hybrid grey forecasting and harmony search approach
CN110705866A (zh) 评估综合能源系统能源生产设备配置与负荷匹配度的方法
Wang et al. Forecasting Crop Yield Variability in Response to Weather Information Using Quantile Random Forest Regression---An Application in Non-Irrigated Upland Cotton
Seifi et al. Designing different sampling plans based on process capability index
CN117235648B (zh) 基于数据处理的钢丝加工全流程集成化管理系统
CN114969017B (zh) 风功率数据清洗方法、清洗装置及预测方法
CN113819581B (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