CN106124119B - 汽机抽汽压力软测量方法 - Google Patents

汽机抽汽压力软测量方法 Download PDF

Info

Publication number
CN106124119B
CN106124119B CN201610621088.8A CN201610621088A CN106124119B CN 106124119 B CN106124119 B CN 106124119B CN 201610621088 A CN201610621088 A CN 201610621088A CN 106124119 B CN106124119 B CN 106124119B
Authority
CN
China
Prior art keywords
steam turbine
pressure
data
extraction pressure
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
CN201610621088.8A
Other languages
English (en)
Other versions
CN106124119A (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.)
Guoneng Shenwan Anqing Power Generation Co ltd
China Shenhua Energy Co Ltd
Shenwan Energy Co Ltd
Original Assignee
China Shenhua Energy Co Ltd
Shenwan Energy Co Ltd
Anhui Anqing Wanjiang Power Generation 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 China Shenhua Energy Co Ltd, Shenwan Energy Co Ltd, Anhui Anqing Wanjiang Power Generation Co Ltd filed Critical China Shenhua Energy Co Ltd
Priority to CN201610621088.8A priority Critical patent/CN106124119B/zh
Publication of CN106124119A publication Critical patent/CN106124119A/zh
Application granted granted Critical
Publication of CN106124119B publication Critical patent/CN106124119B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01LMEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L15/00Devices or apparatus for measuring two or more fluid pressure values simultaneously

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Measuring Fluid Pressure (AREA)

Abstract

本发明涉及测量领域,公开了一种汽机抽汽压力软测量方法,该方法包括:按预定时间间隔采集所述汽机的主蒸汽压力P0和所述汽机的多级抽汽压力中除待测级抽汽压力Pi以外的其它各级抽汽压力Pj,其中j≠i;根据汽机负荷状态筛除所采集的所述汽机的主蒸汽压力P0和所述汽机的其它各级抽汽压力Pj中不稳定的压力值得到筛除后数据;对所述筛除后数据进行标准化处理得到标准化后的数据X(n,m);对所述标准化后的数据X(n,m)进行主元分析PCA选取主元,以确定对应的得分矩阵;以及将所述得分矩阵作为支持向量机的输入来确定所述待测级抽汽压力Pi。通过使用本发明上述的软测量方法,能够更为准确地实现对汽机抽汽压力的软测量。

Description

汽机抽汽压力软测量方法
技术领域
本发明涉及测量领域,具体地,涉及一种汽机抽汽压力软测量方法。
背景技术
在过程工业中,存在一些无法直接测量或测量有很大滞后的变量,需要通过软测量技术建立模型对其进行估计。建模的成功与否依赖于正确的建模数据,但从现场采集的数据由于各种原因不可避免地存在误差。误差分为随机误差和过失误差,随机误差是噪声,过失误差是指失效数据,在实际过程中,过失数据出现的几率很小,但它的存在会使数据品质严重恶化,因此对建模数据进行过失误差侦破,剔除过失数据样本,是建立软测量模型的前提。
汽轮机抽汽压力时常会出现因为前几级抽汽量过大,造成后面抽汽压力降低或者前几级叶片结垢严重,甚至于没有末级抽汽测点而造成压力信息的不可靠性的情况。因此,需要对现有的软测量方法进行改进,以解决上述问题。
发明内容
本发明的目的是提供一种汽机抽汽压力软测量方法,以解决上述现有技术中的问题。
为了实现上述目的,本发明提供一种汽机抽汽压力软测量方法,其中,该方法包括:按预定时间间隔采集所述汽机的主蒸汽压力P0和所述汽机的多级抽汽压力中除待测级抽汽压力Pi以外的其它各级抽汽压力Pj,其中j≠i;根据汽机负荷状态筛除所采集的所述汽机的主蒸汽压力P0和所述汽机的其它各级抽汽压力Pj中不稳定的压力值得到筛除后数据;对所述筛除后数据进行标准化处理得到标准化后的数据X(n,m);对所述标准化后的数据X(n,m)进行主元分析PCA选取主元,以确定对应的得分矩阵;以及将所述得分矩阵作为支持向量机的输入来确定所述待测级抽汽压力Pi
通过上述技术方案,在需要对汽机(汽机回热系统)的某级抽汽压力进行软测量时,可以利用执行了筛除操作后的不同时刻的汽机的主蒸汽压力和汽机其余各级抽汽压力作为主元分析的输入数据,该输入数据在用于主元分析之前可以被标准化,然后对该标准化后的数据进行主元分析选取预定数量的主元,以确定对应的得分矩阵,并将所述得分矩阵作为支持向量机的输入进而确定所述待测级抽汽压力。由此,能够更为准确地实现对汽机抽汽压力的软测量。
本发明的其它特征和优点将在随后的具体实施方式部分予以详细说明。
附图说明
附图是用来提供对本发明的进一步理解,并且构成说明书的一部分,与下面的具体实施方式一起用于解释本发明,但并不构成对本发明的限制。在附图中:
图1是根据本发明的一种实施方式的汽机抽汽压力软测量方法的流程图。
具体实施方式
以下结合附图对本发明的具体实施方式进行详细说明。应当理解的是,此处所描述的具体实施方式仅用于说明和解释本发明,并不用于限制本发明。
在本发明实施方式中,对现有的主元分析方法和支持向量机进行了利用。
其中,统计学习理论是由Vapnik建立的一种专门研究小样本情况下机器学习规律的理论,支持向量机是在这一理论基础上发展而来的一种新的通用学习方法。支持向量机通过结构风险最小化原理来提高泛化能力,它较好地解决了小样本、非线性、高维数、局部极小点等实际问题,已在模式识别、信号处理、函数逼近等领域得到了应用。
主元分析方法是多元统计分析的主要方法之一,通过对过程数据的多元统计分析,可以从大量过程变量的变化中找到影响过程主要监测参数变化或指标降低的主要原因。把主元分析方法用于热工对象或过程,就是利用主元分析的寻找影响过程变化的主要因素的特点,找到在机组工况变化过程中,引起监测参数或运行指标变化的主导因素。
图1是根据本发明的一种实施方式的汽机抽汽压力软测量方法的流程图。
如图1所示,本发明一种实施方式提供的汽机抽汽压力软测量方法包括:
S100,按预定时间间隔采集所述汽机的主蒸汽压力P0和所述汽机的多级抽汽压力中除待测级抽汽压力Pi以外的其它各级抽汽压力Pj,j≠i;
S102,根据汽机负荷状态筛除所采集的所述汽机的主蒸汽压力P0和所述汽机的其它各级抽汽压力Pj中不稳定的压力值得到筛除后数据;
S104,对所述筛除后数据进行标准化处理得到标准化后的数据X(n,m);
S106,对所述标准化后的数据X(n,m)进行主元分析PCA选取主元,以确定对应的得分矩阵;以及
S108,将所述得分矩阵作为支持向量机的输入来确定所述待测级抽汽压力Pi
在本发明实施方式中,以汽机回热系统具有8级抽汽为例,此时i=1,2,…8,j=1,2,…8,但不会出现i与j相等的情况。也就是,假如末级抽汽压力为待测级抽汽压力,则此时i的取值为8,而j的取值为1-7(即,其余各级抽汽压力分别为P1、P2、P3、P4、P5、P6、P7。其他级抽汽压力为待测级抽汽压力的情况与上述情况类似,本发明不再赘述。
其中,本领域技术人员可以根据实际情况对预定时间间隔进行设定,例如,该预定时间间隔可以为1分钟,但本发明不限于此。
通过上述技术方案,在需要对汽机(汽机回热系统)的某级抽汽压力进行软测量时,可以利用执行了筛除操作后的不同时刻的汽机的主蒸汽压力和汽机其余各级抽汽压力作为主元分析的输入数据,该输入数据在用于主元分析之前可以被标准化,然后对该标准化后的数据进行主元分析选取预定数量的主元,以确定对应的得分矩阵,并将所述得分矩阵作为支持向量机的输入进而确定所述待测级抽汽压力。由此,能够更为准确地实现对汽机抽汽压力的软测量。
根据本发明一种实施方式,在该方法中,步骤S102可以包括:
将所述汽机负荷状态为不稳定时对应的所述汽机的主蒸汽压力P0和所述汽机的其它各级抽汽压力Pj筛除,得到所述筛除后数据。
也就是,由于在不同的时刻采集得到了多组主蒸汽压力P0和所述汽机的其它各级抽汽压力Pj,相应地不同的时刻也就可能对应不同的汽机负荷状态,因此可以以不稳定的汽机负荷状态作为不稳定数据的筛除依据利用过滤器删除其中不稳定的数据(即,不稳定的汽机负荷状态所对应的时刻采集的那一组或多组数据被删除,剩余的为稳定数据),以确保后续软测量过程中数据的可靠性。
其中,本领域技术人员可以采用现有技术中已有的汽机负荷状态判断方法对汽机负荷状态进行判断,本发明不对此进行限定。例如,汽机负荷上下波动不超过2MW的情况可以被认为是稳定负荷状态(稳定工况),反之则被认为是不稳定负荷状态(不稳定工况)。
根据本发明一种实施方式,在该方法中,通过下述等式对所述筛除后数据进行标准化处理得到标准化后的数据X(n,m):
在该等式中,Xold(n×m)表示预定时间间隔所采集的筛除后的所述汽机的主蒸汽压力P0和所述汽机的其它各级抽汽压力Pj的矩阵集合,n表示采样次数,而m表示所采集的属性变量的数量;
在该等式中,i=1,2,...,n,j=1,2...,m,average(xold(:,j))代表第j个变量下采样点的均值,std(xold(:,j))代表第j个变量下采样点的标准差。
根据本发明一种实施方式,在该方法中,步骤S106可以包括:
S1060,计算X(n,m)的协方差矩阵、特征值λi和特征向量pi
S1062,对协方差矩阵进行奇异值SVD分解;以及
S1064,选取满足预定累积贡献率的特征向量pi作为主元,以确定对应的得分矩阵。
本领域技术人员可以根据实际情况对预定累积贡献率进行设定,例如所述预定累积贡献率可以为80%,但本发明不限于此。
根据本发明一种实施方式,在该方法中,可以通过下述等式计算X(n,m)的协方差矩阵、特征值λi和特征向量pi
COV(X)pi=λipi
根据本发明一种实施方式,在步骤S108中,将所述得分矩阵中的多组数据划分成两部分:一部分作为训练数据,另一部分作为测试数据。划分比例例如可以为等比例划分,但本发明不限于此。
在本发明中,支持向量机测试部分的计算值即为待测级抽汽压力的软测量所得值。
并且,在步骤S108中可以通过下述函数实现支持向量机的运算。
具体地,在高维特征空间中构造最优线性决策函数y(x)=sgn[wψ(x)+b],采取下述等式目标函数:
在该等式中,约束条件i=1,2,…,n,w为权重因子,C为罚参数,b为偏差值,为映射函数。含约束条件的目标函数通过拉格朗日法转化为无约束的目标函数,如下式所示。
其中,y=[y1,…,yn]T;Iv=[1,…,1]T;a=[a1,…,an]T;Ω={Ωij|i,j=1…n}。K(·)为核函数。优选地,此处可选用径向基核函数。
下面结合实例对本发明一种实施方式提供的汽机抽汽压力软测量方法进行表述。
在该实例中,以某电站600MW机组为例,末级抽汽压力P8为待测级抽汽压力,采集SIS系统中2013年11月12日10点至2013年11月12日21点的P0,...,P7数据值,采集的预定时间间隔为1分钟。
首先,以预定时间间隔1分钟采集P0,...,P7数据值,依据负荷状态采用过滤器删除其中不稳定数据,获得稳定工况下的样本数据,剩余492组样本数据。
接着,对剩余492组样本数据执行标准化处理,并将标准化处理后所得到的样本数据作为PCA的输入数据。
其次,依据PCA算法计算输入数据的协方差矩阵、特征值及特征向量,选取满足累积贡献率80%的特征向量作为主元,以确定对应的得分矩阵M(492,k),其中k表示主元数量。
最后,将得到的矩阵M(492,k)作为支持向量机的输入,支持向量机的输出为末级抽汽压力P8的软测量所得值。其中,前246组样本用作训练数据,后246组样本作为测试数据。
优选地,在该实例中,参数C可以取100,径向基核函数中和宽度可以取0.01,时误差在工程误差的接受范围内。
本领域技术人员应当理解,上述实例仅仅是示例性的,并非用于限定本发明。
以上结合附图详细描述了本发明的优选实施方式,但是,本发明并不限于上述实施方式中的具体细节,在本发明的技术构思范围内,可以对本发明的技术方案进行多种简单变型,这些简单变型均属于本发明的保护范围。
另外需要说明的是,在上述具体实施方式中所描述的各个具体技术特征,在不矛盾的情况下,可以通过任何合适的方式进行组合。为了避免不必要的重复,本发明对各种可能的组合方式不再另行说明。
此外,本发明的各种不同的实施方式之间也可以进行任意组合,只要其不违背本发明的思想,其同样应当视为本发明所公开的内容。

Claims (6)

1.一种汽机抽汽压力软测量方法,其中,该方法包括:
按预定时间间隔采集所述汽机的主蒸汽压力P0和所述汽机的多级抽汽压力中除待测级抽汽压力Pi以外的其它各级抽汽压力Pj,j≠i;
根据汽机负荷状态筛除所采集的所述汽机的主蒸汽压力P0和所述汽机的其它各级抽汽压力Pj中不稳定的压力值得到筛除后数据;
对所述筛除后数据进行标准化处理得到标准化后的数据X(n,m),n表示采样次数,m表示所采集的属性变量的数量;
对所述标准化后的数据X(n,m)进行主元分析PCA选取主元,以确定对应的得分矩阵;以及
将所述得分矩阵作为支持向量机的输入来确定所述待测级抽汽压力Pi
2.根据权利要求1所述的方法,其中,根据汽机负荷状态筛除所采集的所述汽机的主蒸汽压力P0和所述汽机的其它各级抽汽压力Pj中不稳定的压力值得到筛除后数据包括:
将所述汽机负荷状态为不稳定时对应的所述汽机的主蒸汽压力P0和所述汽机的其它各级抽汽压力Pj筛除,得到所述筛除后数据。
3.根据权利要求1或2所述的方法,其中,通过下述等式对所述筛除后数据进行标准化处理得到标准化后的数据X(n,m):
在该等式中,Xold(n×m)表示预定时间间隔所采集的筛除后的所述汽机的主蒸汽压力P0和所述汽机的其它各级抽汽压力Pj的矩阵集合,n表示采样次数,而m表示所采集的属性变量的数量;
在该等式中,i=1,2,...,n,j=1,2...,m,:表示i=1,2,...,n中所有元素,X(i,j)代表第i次采样所采集的第j个属性变量下标准化处理得到标准化后的数据,Xold(i,j)代表第i次采样所采集的第j个属性变量下标准化处理前的数据,Xold(:,j)代表所有次数采样所采集的第j个属性变量下标准化处理前的数据,average(xold(:,j))代表第j个变量下采样点的均值,std(xold(:,j))代表第j个变量下采样点的标准差。
4.根据权利要求3所述的方法,其中,对所述标准化后的数据X(n,m)进行主元分析PCA选取预定数量的主元、以确定对应的得分矩阵包括:
计算X(n,m)的协方差矩阵、特征值λi和特征向量pi
对协方差矩阵进行奇异值分解;以及
选取满足预定累积贡献率的特征向量pi作为主元,以确定对应的得分矩阵。
5.根据权利要求4所述的方法,其中,通过下述等式计算X(n,m)的协方差矩阵、特征值λi和特征向量pi
COV(X)pi=λipi
6.根据权利要求4所述的方法,其中,所述预定累积贡献率为80%。
CN201610621088.8A 2016-08-01 2016-08-01 汽机抽汽压力软测量方法 Active CN106124119B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610621088.8A CN106124119B (zh) 2016-08-01 2016-08-01 汽机抽汽压力软测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610621088.8A CN106124119B (zh) 2016-08-01 2016-08-01 汽机抽汽压力软测量方法

Publications (2)

Publication Number Publication Date
CN106124119A CN106124119A (zh) 2016-11-16
CN106124119B true CN106124119B (zh) 2019-02-12

Family

ID=57254613

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610621088.8A Active CN106124119B (zh) 2016-08-01 2016-08-01 汽机抽汽压力软测量方法

Country Status (1)

Country Link
CN (1) CN106124119B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102004460A (zh) * 2010-11-24 2011-04-06 东北电力大学 一种汽轮机通流部分结垢程度的在线监测方法
CN102692875A (zh) * 2012-06-07 2012-09-26 王卫良 汽轮机主蒸汽压力跟踪优化方法
CN103670536A (zh) * 2013-05-30 2014-03-26 甘肃大唐国际连城发电有限责任公司 一种火力发电厂汽轮机调门流量的调节方法
CN103759769A (zh) * 2014-01-24 2014-04-30 烟台龙源电力技术股份有限公司 燃煤电站再热蒸汽流量软测量方法和装置
CN104048842A (zh) * 2014-05-29 2014-09-17 华中科技大学 一种基于软测量技术的汽轮机热耗率在线监测方法
CN104634489A (zh) * 2014-12-05 2015-05-20 广东电网有限责任公司电力科学研究院 汽轮机的功率检测方法和系统
CN104748807A (zh) * 2014-12-12 2015-07-01 东南大学 一种基于流量校正的电站主蒸汽流量在线计算方法
CN105201564A (zh) * 2015-08-26 2015-12-30 国网河南省电力公司电力科学研究院 一种基于主蒸汽流量的汽轮机滑压优化的控制方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100325071A1 (en) * 2007-08-17 2010-12-23 Institutt For Energiteknikk System and method for empirical ensemble-based virtual sensing of gas emission
KR101360790B1 (ko) * 2009-04-22 2014-02-11 한국전력공사 발전소 계측기 성능 감시 예측 방법

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102004460A (zh) * 2010-11-24 2011-04-06 东北电力大学 一种汽轮机通流部分结垢程度的在线监测方法
CN102692875A (zh) * 2012-06-07 2012-09-26 王卫良 汽轮机主蒸汽压力跟踪优化方法
CN103670536A (zh) * 2013-05-30 2014-03-26 甘肃大唐国际连城发电有限责任公司 一种火力发电厂汽轮机调门流量的调节方法
CN103759769A (zh) * 2014-01-24 2014-04-30 烟台龙源电力技术股份有限公司 燃煤电站再热蒸汽流量软测量方法和装置
CN104048842A (zh) * 2014-05-29 2014-09-17 华中科技大学 一种基于软测量技术的汽轮机热耗率在线监测方法
CN104634489A (zh) * 2014-12-05 2015-05-20 广东电网有限责任公司电力科学研究院 汽轮机的功率检测方法和系统
CN104748807A (zh) * 2014-12-12 2015-07-01 东南大学 一种基于流量校正的电站主蒸汽流量在线计算方法
CN105201564A (zh) * 2015-08-26 2015-12-30 国网河南省电力公司电力科学研究院 一种基于主蒸汽流量的汽轮机滑压优化的控制方法

Also Published As

Publication number Publication date
CN106124119A (zh) 2016-11-16

Similar Documents

Publication Publication Date Title
CN110596492B (zh) 一种基于粒子群算法优化随机森林模型的变压器故障诊断方法
US8868985B2 (en) Supervised fault learning using rule-generated samples for machine condition monitoring
CN102829974B (zh) 一种基于lmd和pca的滚动轴承状态辨识方法
CN109193650B (zh) 一种基于高维随机矩阵理论的电网薄弱点评估方法
CN102778355B (zh) 一种基于emd和pca的滚动轴承状态辨识方法
CN109344517A (zh) 一种新能源汽车的高压绝缘故障诊断方法
CN107704992A (zh) 输电线路雷击风险评估的方法及装置
CN105784556A (zh) 一种基于自组织模糊神经网络的空气细颗粒物pm2.5软测量方法
CN108090628A (zh) 一种基于pso-lssvm算法的粮情安全检测分析方法
CN107065842A (zh) 一种基于粒子群优化核独立元分析模型的故障检测方法
CN110442911B (zh) 一种基于统计机器学习的高维复杂系统不确定性分析方法
CN106570633A (zh) 一种基于权重设计的指标体系评估方法
CN105068928A (zh) 一种基于复杂网络理论的软件测试用例生成方法
CN110175682A (zh) 一种基于混沌粒子群的优化核主元分析故障监测方法
CN112926625A (zh) 一种卫星辐射数据的偏差影响因子分析方法
CN105956318A (zh) 基于改进分裂 h-k 聚类方法的风电场机群划分方法
CN108171271B (zh) 一种设备劣化早期预警方法和系统
CN101446828A (zh) 一种非线性过程质量预报方法
Bayarri et al. A Bayesian sequential look at u-control charts
CN106124119B (zh) 汽机抽汽压力软测量方法
CN114580828A (zh) 数控机床运输过程安全性评价方法
CN105590167A (zh) 电场多元运行数据分析方法及装置
CN103389360A (zh) 基于概率主元回归模型的脱丁烷塔丁烷含量软测量方法
Gao et al. Fused machine learning based on Boosting algorithm for noise recognition of moving parts in sealed cavity
CN107179758A (zh) 一种动态信号参数辨识方法及系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CP03 Change of name, title or address

Address after: 100011 Beijing Dongcheng District, West Binhe Road, No. 22

Patentee after: China Shenhua Energy Co.,Ltd.

Patentee after: SHENWAN ENERGY Co.,Ltd.

Patentee after: Guoneng Shenwan Anqing Power Generation Co.,Ltd.

Address before: 100011 Shenhua building, 22 West Binhe Road, Dongcheng District, Beijing

Patentee before: China Shenhua Energy Co.,Ltd.

Patentee before: SHENWAN ENERGY Co.,Ltd.

Patentee before: ANHUI ANQING WANJIANG POWER GENERATION Co.,Ltd.

CP03 Change of name, title or address