CN105320167B - 火电机组主汽压力在线优化方法 - Google Patents
火电机组主汽压力在线优化方法 Download PDFInfo
- Publication number
- CN105320167B CN105320167B CN201410854217.9A CN201410854217A CN105320167B CN 105320167 B CN105320167 B CN 105320167B CN 201410854217 A CN201410854217 A CN 201410854217A CN 105320167 B CN105320167 B CN 105320167B
- Authority
- CN
- China
- Prior art keywords
- value
- point
- vapour pressure
- data
- cut
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 71
- 238000005457 optimization Methods 0.000 title claims abstract description 40
- 238000010977 unit operation Methods 0.000 claims abstract description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 36
- 238000005070 sampling Methods 0.000 claims description 17
- 230000004907 flux Effects 0.000 claims description 13
- 238000000605 extraction Methods 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 10
- 239000007921 spray Substances 0.000 claims description 10
- 230000000153 supplemental effect Effects 0.000 claims description 9
- 238000012935 Averaging Methods 0.000 claims description 4
- 239000000203 mixture Substances 0.000 claims description 4
- 238000009825 accumulation Methods 0.000 claims description 3
- YNBADRVTZLEFNH-UHFFFAOYSA-N methyl nicotinate Chemical compound COC(=O)C1=CC=CN=C1 YNBADRVTZLEFNH-UHFFFAOYSA-N 0.000 claims description 3
- 230000011218 segmentation Effects 0.000 claims description 3
- 241001269238 Data Species 0.000 claims description 2
- 238000003303 reheating Methods 0.000 claims 1
- 238000004364 calculation method Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 2
- 238000010998 test method Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 101100129500 Caenorhabditis elegans max-2 gene Proteins 0.000 description 1
- 241000254158 Lampyridae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 102220049405 rs147669920 Human genes 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
Landscapes
- Feedback Control In General (AREA)
Abstract
本发明提供一种火电机组主汽压力在线优化方法,包括:获取火电机组运行参数实时数据;预处理运行参数实时数据;对预处理后的数据进行稳定性判断,得到准稳态工况数据;计算准稳态工况数据对应的热耗值;用最小二乘法对准稳态工况数据中的机组负荷特征值下的主汽压力值及其对应的热耗值进行曲线拟合,确定机组负荷特征值下的最优主汽压力值;对获得的最优主汽压力值及其对应的机组负荷值进行曲线拟合,得到机组主汽压力优化目标曲线。本发明能够提高主汽压力优化结果的准确度,降低寻找主汽压力优化值的复杂度。
Description
技术领域
本发明涉及一种工程优化方法,特别是一种火电机组主汽压力在线优化方法。
背景技术
随着电力工业和社会经济的不断发展,电网的峰谷差不断变大,火电机组越来越多地参与到电网调峰中。与定压运行相比,滑压运行机组在部分负荷下,具有良好的经济性和安全性。随着主汽压力降低,机组相对内效率不断提高,而循环热效率不断降低,因此,在给定运行条件下机组存在最优主汽压力值。
目前,获得火电机组最优运行初压的方法主要有两种:基于机组变工况理论的计算方法和机组主蒸汽压力寻优试验方法。由于变工况计算过程过多依赖理论模型,且机组设备及系统状态偏离设计工况,因此,基于机组变工况理论计算的方法不确定度较高,难以实际推广。而机组主蒸汽压力寻优试验方法,需要安排大量的现场试验,同时,随着机组的性能变化,需要重新试验,费时费力。
发明内容
本发明提出了一种火电机组主汽压力在线优化方法,解决现有技术中主汽压力优化复杂,优化结果不准确的缺陷,本发明能够实现对最优主汽压力的在线跟踪,为现场运行人员提供帮助。
本发明提供一种火电机组主汽压力在线优化方法,其中,在线优化方法包括:
获取火电机组运行参数实时数据;
预处理运行参数实时数据;
对预处理后的数据进行稳定性判断,得到准稳态工况数据;
计算准稳态工况数据对应的热耗值;
用最小二乘法对准稳态工况数据中的机组负荷特征值下的主汽压力值及其对应的热耗值进行曲线拟合,确定机组负荷特征值下的最优主汽压力值;
对所述最优主汽压力值及其对应的机组负荷值进行曲线拟合,得到机组主汽压力优化目标曲线;
其中,确定机组负荷特征值下的最优主汽压力值,包括:将热耗最小值对应的主汽压力值确定为最优主汽压力值;
其中,若求得的热耗最小值等于离散点中的热耗值,则该最小热耗值对应的主汽压力值为该机组负荷特征值下的最优主汽压力值;若求得的热耗最小值不等于离散点中的热耗值,则在用于拟合曲线用的离散点中选取与热耗最小值差值最小的离散点的热耗值对应的主汽压力值为该机组负荷特征值下的最优主汽压力值。
一实施例中,预处理运行参数实时数据的过程为:
对机组运行参数实时数据进行上下限判断,将越限的参数实时数据用其默认值代替,从而保证系统的稳定运行;
若参数有多个测点,则对参数的多个测点数据进行平均处理或累加处理。
一实施例中,稳定性判断过程为:
设定一滑动数据窗口,所述滑动数据窗口的时长为N个采样周期;
每隔一个采样周期读取一组预处理后的数据至滑动数据窗口中;
当滑动数据窗口内预处理后的数据都满足稳定性判据时,则将滑动数据窗口内最后一个采样周期读取的预处理后的数据作为当前滑动数据窗口下的准稳态工况数据输出;
当滑动数据窗口内有一组参数实时数据不满足稳定性判据时,则当前滑动数据窗口下没有准稳态工况数据输出;
更新滑动数据窗口内的数据,移出滑动数据窗口内最早读入的预处理后的数据,读取下一采样周期的预处理后的数据至滑动数据窗体内继续进行稳定性判断。
一实施例中,当滑动数据窗体内预处理后的数据都满足稳定性判据时,当前滑动数据窗口下的准稳态工况数据输出还可为滑动数据窗体内的预处理后的每种参数数据的加权处理值。
一实施例中,稳定性判据为:
判断滑动数据窗体内预处理后的特征参数数据的最大值与该特征参数数据的最小值的差值是否小于设定的阈值,若差值小于所述的阈值,则当前滑动数据窗口内的预处理后的数据满足稳定性,若差值大于或等于设定的阈值,则当前预处理后的数据不满足稳定性。
一实施例中,所述参数包括:机组负荷、主汽压力、给水泵运行方式、主蒸汽流量、主给水流量、过热器减温水量、冷再热蒸汽流量、热再热减温水量、汽轮机输出功率、供热抽汽量。
一实施例中,火电机组热耗计算公式为:
其中,Fms为主蒸汽流量;Hms为主蒸汽焓;Ffw为主给水流量;Hfw为主给水焓;Fshsp为过热器减温水量;Hshop为过热器减温水焓;Fcrh为冷再热蒸汽流量;Hcrh为冷再热蒸汽焓;Hhrh为热再热蒸汽焓;Frhsp为热再热减温水流量;Hrhsp为热再热减温水焓;P为汽轮机输出功率;Fcq供热抽汽量;Hcq为供热抽汽焓;Hhs为供热抽汽回水焓。
一实施例中,计算完准稳态工况数据对应的热耗值之后,还包括按给水泵运行方式不同对准稳态工况数据进行分类;
在同一机组负荷特征值下,用最小二乘法对每种给水泵运行方式的主汽压力值及其对应的热耗值进行曲线拟合,确定该机组负荷特征值下不同给水泵运行方式的一个最优主汽压力值;
对同一给水泵运行方式下的最优主汽压力值及各其对应的机组负荷值进行曲线拟合,得到不同给水泵运行方式的机组主汽压力优化目标曲线。
一实施例中,定义主汽压力优化值为纵坐标,定义机组负荷为横坐标,机组主汽压力优化目标曲线拟合过程为:
从横坐标中选取六个分割点:第一分割点、第二分割点、第三分割点、第四分割点、第五分割点及第六分割点,
所述六个分割点将机组负荷分为五个区间:大于第一分割点小于或等于第二分割点为第一区间,大于第二分割点小于或等于第三分割点为第二区间,大于第三分割点小于或等于第四分割点为第三区间,大于第四分割点小于或等于第五分割点为第四区间,大于第五分割点小于或等于第六分割点为第五区间;
定义主汽压力目标优化曲线组成:第一区间内为第一直线、第二区间内为第一圆弧、第三区间内为第二直线、第四区间内为第二圆弧、第五区间内为第三直线;用最小二乘法拟合第一直线、第二直线及第三直线,得到第一、第二及第三直线的表达式;
根据第一及第二直线的表达式得到第一及第二直线间的第一角分线表达式,根据所述第二及第三直线的表达式得到第二及第三直线间的第二角分线表达式;
第一圆弧的圆心位于第一角分线上,设置第一圆心,以第一圆心至第一直线或第二直线的垂直距离为半径画弧,与第一直线及第二直线的交点分别记为第七分割点,第八分割点,交点第七、第八分割点取代第二分割点,第三分割点作为新的分割点,交点第七、第八分割点之间的圆弧为拟合的第一圆弧;
第二圆弧的圆心位于第二角分线上,设置第二圆心,以第二圆心至第二直线或第三直线的垂直距离为半径画弧,与第二直线及第三直线的交点分别记为第九分割点,第十分割点,交点第九、第十分割点取代第四分割点,第五分割点作为新的分割点,交点第九、第十分割点之间的圆弧为拟合的第二圆弧;
根据新的分割点,调整第一圆心及第二圆心位置,重新拟合第一直线、第一圆弧、第二直线、第二圆弧及第三直线;
新的分割点变动小于预定值或迭代次数达到设定值时,停止拟合,最后一次拟合的第一直线、第一圆弧、第二直线、第二圆弧及第三直线为主汽压力优化目标曲线。
综上所述技术方案,采用本发明提供的火电机组主汽压力在线优化方法能够准确获得机组主汽压力优化目标曲线,本发明提供的主汽压力目标曲线拟合方法能够得到光滑的主汽压力优化目标曲线,使得优化压力曲线拐角处的主汽压力值与实际最优的主汽压力值更加接近。能为后续控制过程提供准确的主汽压力目标值。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例的火电机组主汽压力在线优化方法的流程图;
图2为本发明一实施例的拟合主汽压力目标曲线的原理图。
附图符号说明:
X0~X5、X1'~X4':分割点;
L1:第一直线;
L2:第二直线;
L3:第三直线;
L4:第一角分线;
L5:第二角分线;
C1:第一圆弧;
C2:第二圆弧。
具体实施方式
为了使本发明的技术特点及效果更加明显,下面结合附图对本发明的技术方案做进一步说明,本发明也可有其他不同的具体实例来加以说明或实施,任何本领域技术人员在权利要求范围内做的等同变换均属于本发明的保护范畴。
如图1所示,图1为本发明实施例的火电机组主汽压力在线优化方法的流程图,其中,该在线优化方法包括:步骤S11,获取火电机组运行参数实时数据;步骤S12,预处理所述运行参数实时数据;步骤S13,对预处理后的数据进行稳定性判断,得到准稳态工况数据;步骤S14,计算准稳态工况数据对应的热耗值;步骤S15,用最小二乘法对准稳态工况数据中的机组负荷特征值下的主汽压力值及其对应的热耗值进行曲线拟合,确定机组负荷特征值下的最优主汽压力值;步骤S16,对最优主汽压力值及其对应的机组负荷值进行曲线拟合,得到机组主汽压力优化目标曲线。
本发明一实施例中,从DCS中获取火电机组的运行数据,本发明另一实施例中,从SIS中获取机组的运行数据,获取运行参数实时数据的方式视各电厂的具体情况而定,本发明对此并不做限制。
优选的,获取的运行参数实时数据包括机组负荷、主汽压力、主汽温度、给水泵运行方式、主蒸汽流量、主给水流量、过热器减温水量、冷再热蒸汽流量、热再热减温水量、汽轮机输出功率、供热抽汽量。本发明其它实施例中,获取的数据还包括给水温度、给水压力、再冷温度等,需要说明的是,本发明对获取的实时参数数据并不做限制,凡是影响主汽压力值的参数都可包括于获取的运行参数实时数据中,可根据实际工程需要选择采样的运行参数。获取的运行参数实时数据可先存储至一实时数据平台上,以便后续调用。
步骤S12中的预处理运行参数实时数据的过程为:
对机组运行参数实时数据进行上下限判断,将越限的参数实时数据用其默认值代替,从而保证系统的稳定运行。
若某一或某些参数有多个测点,则需对该种参数实时数据进行平均处理或累加处理。
在本实施例中,根据系统配置文件,先对获取的实时数据进行上下限判断,保留满足上下限要求的数据,用于后续计算。若同一参数的数据有多个测点,则需对同种参数数据进行平均处理或累加处理,例如,主汽压力包括主汽压力1、主汽压力2及主汽压力3三个测点,减温水量包括减温水量1、减温水量2、减温水量3及减温水量4四个测点,预处理后的主汽压力等于主汽压力1、主汽压力2及主汽压力3的平均值,减温水量等于减温水量1、减温水量2、减温水量3及减温水量4的累加和。通过对获取的参数实时数据进行上下限判断,可以有效的去除不符合条件的数据,以免影响后续热耗计算的准确度,对同种参数的多个测点进行四则运算处理,使得各参数的数值更加接近真实值,从而提高后续过程的准确度。
优选的,步骤S13中的稳定性判断过程为:
设定一个滑动数据窗口,用于筛选预处理后的参数实时数据,滑动数据窗口的时长为N个采样周期,需要说明的是,本发明对采样周期及滑动数据窗口的时长并不做限制,可根据获取的参数实时数据及优化需要适当调节采样周期及滑动数据窗口的时长。
每隔一个采样周期读取一组预处理后的数据至滑动数据窗口中,当滑动数据窗口内预处理后的数据都满足稳定性判据时,则将滑动数据窗口内最后一个采样周期读取的预处理后的数据作为当前滑动数据窗口下的准稳态工况数据输出,当滑动数据窗口内有一组参数实时数据不满足稳定性时,则当前滑动数据窗口下没有准稳态工况数据输出,更新滑动数据窗口内的数据,移出滑动数据窗口内最早读入的预处理后的数据,读取下一采样周期的预处理后的数据至滑动数据窗体内继续进行稳定性判断。
需要说明的是,当滑动数据窗体内预处理后的数据都满足稳定性判据时,当前滑动数据窗口下的稳态参数数据输出还可以为滑动数据窗体内的预处理后的每种参数数据的加权处理值。以机组负荷为例,滑动数据窗口的时长为10个采样周期,当前滑动数据窗口下的机组负荷为X1、X2、X3、X4、X5、X6、X7、X8、X9、X10,各机组负荷对应的权值为A1、A2、A3、A4、A5、A6、A7、A8、A9、A10,对机组负荷值进行加权处理后得X=A1*X1+A2*X2+A3*X3+A4*X4+A5*X5+A6*X6+A7*X7+A8*X8+A9*X9+A10*X10,其中,A1+A2+A3+A4+A5+A6+A7+A8+A9+A10=1,X为当前滑动数据窗口下的机组负荷的准稳态输出,同理,滑动数据窗口内的其他数据的处理方式与上述机组负荷的处理方式相同,此处不再赘述,同时,本发明对加权处理时的权值并不做限定,可根据实际需要选择合适的权值。
优选的,稳定性判据为:判断滑动数据窗体内预处理后的特征参数数据的最大值与该特征参数数据的最小值的差值是否小于该特征参数设定的阈值,若差值小于设定的阈值,则当前滑动数据窗口内的预处理后的数据满足稳定性,若差值大于或等于设定的阈值,则当前滑动数据窗口内的预处理后的数据不满足稳定性。
以主汽压力及主汽温度为例,定义滑动数据窗体内主汽压力的最大值为Max1,主汽压力最小值为Min1,主汽压力满足稳定性的门限值为0.2MPa,主汽温度最大值为Max2,主汽温度最小值为Min2,主汽温度满足稳定性的门限值为5℃,则滑动数据窗体内数据满足稳定性判断条件为Max1-Min1<0.2MPa,Max2-Min2<5℃,若滑动数据窗体内的主汽压力及主汽温度都满足,则当前滑动数据窗口内的参数数据满足稳定性,若滑动数据窗体内有一组参数不满足条件,则当前滑动数据窗口内的参数数据不满足稳定性,将滑动数据窗口中最早读取的预处理后的数据移出,读取下一采样周期的预处理后的数据至滑动数据窗体内,再重复进行稳定性判断。
当然,滑动数据窗体内数据的稳定性判断还可以选取其他稳定性判断方法,对不同的参数也可选取不同的稳定性判断方法,本发明并不以上述实施例为限,同时,本发明对设定的阈值具体取值也不做限制。预处理后的参数实时数据通过上述稳定性判断,能够保证其稳定性,去除那些波动幅度较大的点,进一步提高热耗计算的准确度,便于准确找出最优的主汽压力值点。
本发明一实施例中,火电机组热耗计算公式为:
其中,Fms为主蒸汽流量;Hms为主蒸汽焓;Ffw为主给水流量;Hfw为主给水焓;Fshsp为过热器减温水量;Hshop为过热器减温水焓;Fcrh为冷再热蒸汽流量;Hcrh为冷再热蒸汽焓;Hhrh为热再热蒸汽焓;Frhsp为热再热减温水流量;Hrhsp为热再热减温水焓;P为汽轮机输出功率;Fcq供热抽汽量;Hcq为供热抽汽焓;Hhs为供热抽汽回水焓。
优选的,确定同一机组负荷特征值下的一个最优主汽压力值的方法可以为最小值法或萤火虫寻优法,热耗最小值对应的主汽压力值即为最优的主汽压力值。
本发明还可以对同一机组负荷特征值下拟合曲线的方式求热耗最小值,若求得的热耗最小值等于离散点(如图2所示离散点※)中的热耗值,则该最小热耗值对应的主汽压力值为该机组负荷特征值下的最优主汽压力值;若求得的热耗最小值不等于离散点中的热耗值,则在用于拟合曲线用的离散点中选取与热耗最小值差值最小的离散点的热耗值对应的主汽压力值为该机组负荷特征值下的最优主汽压力值。需要说明的是,机组负荷特征值并非一个机组负荷值,而是一定范围内的机组负荷值,本发明对机组负荷特征值的范围大小并不做限定,可视获得的准稳态工况数据具体情况而定。
如图2所示,图2为本发明一实施例的拟合主汽压力目标曲线的原理图,本实施例中,定义主汽压力优化值为纵坐标,定义机组负荷为横坐标,图2中的离散点※为各区域内的最优的主汽压力值,机组主汽压力优化目标曲线拟合过程为:
从横坐标中选取六个分割点:X0~X5,分割点X0~X5将机组负荷分为五个区间:X0<X<=X1为第一区间,X1<X<=X2为第二区间,X2<X<=X3为第三区间,X3<X<=X4为第四区间,X04<X<=X5为第五区间。
定义主汽压力优化目标曲线组成:第一区间内为第一直线L1、第二区间内为第一圆弧C1、第三区间内为第二直线L2、第四区间内为第二圆弧C2、第五区间内为第三直线L3。
定义主汽压力优化曲线组成:第一区间内为第一直线、第二区间内为第一圆弧、第三区间内为第二直线、第四区间内为第二圆弧、第五区间内为第三直线。
用最小二乘法拟合第一直线L1、第二直线L2及第三直线L3,得到第一、第二及第三直线的表达式。
根据第一直线L1及第二直线L2的表达式得到第一直线L1及第二直线L2之间的第一角分线L4表达式,根据第二直线L2及第三直线L3的表达式得到第二直线L2及第三直线L3之间的第二角分线L5表达式。
第一圆弧C1的圆心位于第一角分线L4上,设置第一圆心,以第一圆心至第一直线L1或第二直线L2的垂直距离为半径画弧,与第一直线L1及第二直线L2的交点分别记为X1'、X2',交点X1'、X2'取代分割点X1、X2作为新的分割点,交点X1'、X2'之间的圆弧为拟合的第一圆弧C1。
第二圆弧C2的圆心位于第二角分线L5上,设置第二圆心,以第二圆心至第二直线L2或第三直线L3的垂直距离为半径画弧,与第二直线L2及第三直线L3的交点分别记为X3'、X4',交点X3'、X4'取代分割点X3、X4作为新的分割点,交点X3'、X4'之间的圆弧为拟合的第二圆弧C2。
根据新的分割点,调整第一圆心及第二圆心位置,重新拟合第一直线L1、第一圆弧C1、第二直线L2、第二圆弧C2及第三直线L3。
当新的分割点变动小于预定值或迭代次数达到设定值时,停止拟合,最后一次拟合的第一直线、第一圆弧、第二直线、第二圆弧及第三直线为主汽压力优化目标曲线。
需要说明的是,本发明对分割点的选取并不做限制,可根据离散点的分散程度参考选取。通过本实施例提供的拟合主汽压力优化曲线的方法,能够得到光滑的主汽压力优化曲线,使得优化压力曲线拐角处的主汽压力值与实际最优的主汽压力值更加接近。本发明得到的主汽压力优化目标曲线可通过标准硬件接线方式(4-20mA)提供DCS使用,能实时提供主汽压力目标值。
本发明一实施例中,计算完准稳态工况数据对应的热耗值之后,还包括按给水泵运行方式不同对准稳态工况数据进行分类;
在同一机组负荷特征值下,用最小二乘法对每种给水泵运行方式的主汽压力值及其对应的热耗值进行曲线拟合,确定该机组负荷特征值下不同给水泵运行方式的一个最优主汽压力值。
对同一给水泵运行方式下的最优主汽压力值及其对应的机组负荷值进行曲线拟合,得到不同给水泵运行方式的机组主汽压力优化目标曲线。由于不同工况下,给水泵运行方式对机组整体热耗具有较大影响,为了保证本优化方法的优化结果对机组运行具有指导意义,因此,根据当前给水泵的运行方式,给出不同的优化主汽压力曲线。
以上所述仅用于说明本发明的技术方案,任何本领域普通技术人员均可在不违背本发明的精神及范畴下,对上述实施例进行修饰与改变。因此,本发明的保护范围应视权利要求范围为准。
Claims (5)
1.一种火电机组主汽压力在线优化方法,其特征在于,所述在线优化方法为:
获取火电机组运行参数实时数据;
预处理所述运行参数实时数据;
对预处理后的数据进行稳定性判断,得到准稳态工况数据;
计算准稳态工况数据对应的热耗值;
用最小二乘法对准稳态工况数据中的机组负荷特征值下的主汽压力值及其对应的热耗值进行曲线拟合,确定机组负荷特征值下的最优主汽压力值;
对所述最优主汽压力值及其对应的机组负荷值进行曲线拟合,得到机组主汽压力优化目标曲线;
计算完准稳态工况数据对应的热耗值之后,还包括按给水泵运行方式不同对准稳态工况数据进行分类;
在同一机组负荷特征值下,用最小二乘法对每种给水泵运行方式的主汽压力值及其对应的热耗值进行曲线拟合,确定该机组负荷特征值下不同给水泵运行方式的最优主汽压力值;
对不同给水泵运行方式下的所述最优主汽压力值及其对应的机组负荷值进行曲线拟合,得到不同给水泵运行方式的机组主汽压力优化目标曲线;
预处理所述运行参数实时数据的过程为:
对机组运行参数实时数据进行上下限判断,将越限的参数实时数据用其默认值代替,从而保证系统的稳定运行;
若参数有多个测点,则对参数的多个测点数据进行平均处理或累加处理;
其中,稳定性判断过程为:
设定一滑动数据窗口,所述滑动数据窗口的时长为N个采样周期;
每隔一个采样周期读取一组预处理后的数据至滑动数据窗口中;
当所述滑动数据窗口内预处理后的数据都满足稳定性判据时,则将所述滑动数据窗口内最后一个采样周期读取的预处理后的数据作为当前滑动数据窗口下的准稳态工况数据输出;
当所述滑动数据窗口内有一组参数实时数据不满足稳定性时,则当前滑动数据窗口下没有准稳态工况数据输出;
更新滑动数据窗口内的数据,移出所述滑动数据窗口内最早读入的预处理后的数据,读取下一采样周期的预处理后的数据至滑动数据窗体内继续进行稳定性判断;
其中,当所述滑动数据窗体内预处理后的数据都满足稳定性判据时,当前滑动数据窗口下的准稳态工况数据输出还可为所述滑动数据窗体内的预处理后的每种参数数据的加权处理值;
其中,确定机组负荷特征值下的最优主汽压力值,包括:将热耗最小值对应的主汽压力值确定为最优主汽压力值;
其中,若求得的热耗最小值等于离散点中的热耗值,则该最小热耗值对应的主汽压力值为该机组负荷特征值下的最优主汽压力值;若求得的热耗最小值不等于离散点中的热耗值,则在用于拟合曲线用的离散点中选取与热耗最小值差值最小的离散点的热耗值对应的主汽压力值为该机组负荷特征值下的最优主汽压力值。
2.如权利要求1所述的火电机组主汽压力在线优化方法,其特征在于,所述稳定性判据为:
判断滑动数据窗体内预处理后的特征参数数据的最大值与该特征参数数据的最小值的差值是否小于设定的阈值,若差值小于所述设定的阈值,则当前滑动数据窗口内的预处理后的数据满足稳定性,若差值大于或等于所述设定的阈值,则当前滑动数据窗口内的预处理后的数据不满足稳定性。
3.如权利要求1所述的火电机组主汽压力在线优化方法,其特征在于,所述参数包括:机组负荷、主汽压力、给水泵运行方式、主蒸汽流量、主给水流量、过热器减温水量、冷再热蒸汽流量、热再热减温水量、汽轮机输出功率、供热抽汽量。
4.如权利要求3所述的火电机组主汽压力在线优化方法,其特征在于,火电机组热耗计算公式为:
其中,Fms为主蒸汽流量;Hms为主蒸汽焓;Ffw为主给水流量;Hfw为主给水焓;Fshsp为过热器减温水量;Hshop为过热器减温水焓;Fcrh为冷再热蒸汽流量;Hcrh为冷再热蒸汽焓;Hhrh为热再热蒸汽焓;Frhsp为热再热减温水流量;Hrhsp为热再热减温水焓;P为汽轮机输出功率;Fcq供热抽汽量;Hcq为供热抽汽焓;Hhs为供热抽汽回水焓。
5.如权利要求1所述的火电机组主汽压力在线优化方法,其特征在于,定义主汽压力优化值为纵坐标,定义机组负荷为横坐标,机组主汽压力优化目标曲线拟合过程为:
从横坐标中选取六个分割点:第一分割点、第二分割点、第三分割点、第四分割点、第五分割点及第六分割点,
所述六个分割点将机组负荷分为五个区间:大于第一分割点小于或等于第二分割点为第一区间,大于第二分割点小于或等于第三分割点为第二区间,大于第三分割点小于或等于第四分割点为第三区间,大于第四分割点小于或等于第五分割点为第四区间,大于第五分割点小于或等于第六分割点为第五区间;
定义主汽压力优化目标曲线组成:第一区间内为第一直线、第二区间内为第一圆弧、第三区间内为第二直线、第四区间内为第二圆弧、第五区间内为第三直线;
用最小二乘法拟合第一直线、第二直线及第三直线,得到第一、第二及第三直线的表达式;
根据第一及第二直线的表达式得到第一及第二直线间的第一角分线表达式,根据所述第二及第三直线的表达式得到第二及第三直线间的第二角分线表达式;
第一圆弧的圆心位于第一角分线上,设置第一圆心,以第一圆心至第一直线或第二直线的垂直距离为半径画弧,与第一直线及第二直线的交点分别记为第七分割点,第八分割点,交点第七、第八分割点分别取代第二分割点,第三分割点作为新的分割点,交点第七、第八分割点之间的圆弧为拟合的第一圆弧;
第二圆弧的圆心位于第二角分线上,设置第二圆心,以第二圆心至第二直线或第三直线的垂直距离为半径画弧,与第二直线及第三直线的交点分别记为第九分割点,第十分割点,交点第九、第十分割点分别取代第四分割点,第五分割点作为新的分割点,交点第九、第十分割点之间的圆弧为拟合的第二圆弧;
根据新的分割点,调整第一圆心及第二圆心位置,重新拟合第一直线、第一圆弧、第二直线、第二圆弧及第三直线;
新的分割点变动小于预定值或迭代次数达到设定值时,停止拟合,最后一次拟合的第一直线、第一圆弧、第二直线、第二圆弧及第三直线为主汽压力优化目标曲线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410854217.9A CN105320167B (zh) | 2014-12-31 | 2014-12-31 | 火电机组主汽压力在线优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410854217.9A CN105320167B (zh) | 2014-12-31 | 2014-12-31 | 火电机组主汽压力在线优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105320167A CN105320167A (zh) | 2016-02-10 |
CN105320167B true CN105320167B (zh) | 2019-08-02 |
Family
ID=55247739
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410854217.9A Active CN105320167B (zh) | 2014-12-31 | 2014-12-31 | 火电机组主汽压力在线优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105320167B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111352400B (zh) * | 2020-03-13 | 2021-05-07 | 中国大唐集团科学技术研究院有限公司华东电力试验研究院 | 一种供热机组滑压运行优化方法及滑压运行控制系统 |
CN112182779B (zh) * | 2020-10-29 | 2023-08-29 | 同济大学 | 一种基于模拟退火算法的汽油机在线标定方法 |
CN112576326B (zh) * | 2020-12-07 | 2023-10-13 | 润电能源科学技术有限公司 | 一种火电机组滑压优化运行控制方法、装置及设备 |
CN118070083A (zh) * | 2024-01-11 | 2024-05-24 | 太极计算机股份有限公司 | 基于人工智能算法的汽轮机诊断方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101719230A (zh) * | 2009-12-18 | 2010-06-02 | 东南大学 | 大型汽轮机组主汽压力优化方法 |
CN102661176A (zh) * | 2012-06-04 | 2012-09-12 | 哈尔滨工业大学 | 一种以主蒸汽流量作为调度变量获取汽轮机滑压曲线方法 |
CN103646162A (zh) * | 2013-11-12 | 2014-03-19 | 国电科学技术研究院 | 一种汽轮机主蒸汽压力运行曲线及其优化方法 |
-
2014
- 2014-12-31 CN CN201410854217.9A patent/CN105320167B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101719230A (zh) * | 2009-12-18 | 2010-06-02 | 东南大学 | 大型汽轮机组主汽压力优化方法 |
CN102661176A (zh) * | 2012-06-04 | 2012-09-12 | 哈尔滨工业大学 | 一种以主蒸汽流量作为调度变量获取汽轮机滑压曲线方法 |
CN103646162A (zh) * | 2013-11-12 | 2014-03-19 | 国电科学技术研究院 | 一种汽轮机主蒸汽压力运行曲线及其优化方法 |
Non-Patent Citations (2)
Title |
---|
Forecasting of turbine heat rate with online least squares support vector machine based on gravitational search algorithm;Zhang W,等;《Knowledge-Based Systems》;20131231;第39卷;第34-44页 |
基于偏最小二乘算法的热耗率回归分析;张文琴,等;《现代电力》;20091031;第26卷(第5期);第56-59页 |
Also Published As
Publication number | Publication date |
---|---|
CN105320167A (zh) | 2016-02-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105224735B (zh) | 发电机组能效分析方法 | |
CN105320167B (zh) | 火电机组主汽压力在线优化方法 | |
CN102063566B (zh) | 一种水力热力耦合仿真模型的多气源蒸汽管网计算系统 | |
CN105303032B (zh) | 影响发电机组能效的客观因素分析方法 | |
CN103336891B (zh) | 一种用于配电网状态估计的伪量测生成方法 | |
CN108509752A (zh) | 一种结合现场数据的过热器机理模型动态参数整定方法 | |
CN105317476B (zh) | 一种基于特征通流面积的汽轮机流量曲线辨识及优化方法 | |
CN105275508A (zh) | 一种基于功率值计算的汽轮机流量曲线辨识及优化方法 | |
CN110059388B (zh) | 核电厂冷凝器沿工质流动方向的参数分布建模仿真方法 | |
CN110941186A (zh) | 基于神经网络与万有引力搜索算法的汽温控制寻优方法 | |
CN112417764A (zh) | 一种面向锅炉特种设备蒸汽量预测的k近邻回归预测方法 | |
CN105240058B (zh) | 一种基于喷嘴流量计算的汽轮机流量曲线辨识及优化方法 | |
CN106761967A (zh) | 机侧蒸汽参数测量偏差对机组耗煤成本的评估方法及系统 | |
CN104122291A (zh) | 超超临界火电机组水冷壁向工质传热速率的实时辨识方法 | |
CN108595723A (zh) | 一种锅炉暖风器回热量计算方法及装置 | |
CN109002741A (zh) | 一种压水堆核电机组一、二回路系统传递功率模拟方法及系统 | |
CN109709911B (zh) | 一种火电机组循环工质外漏在线测量方法及测量系统 | |
CN113885607B (zh) | 一种汽温控制方法、装置、电子设备和计算机存储介质 | |
CN108615121B (zh) | 一种基于多因素影响的热电负荷分配方法及系统 | |
CN107451304B (zh) | 再热器传递函数模型的机理建模计算方法 | |
CN109973159A (zh) | 确定汽轮机低压缸流量的方法和装置及机器可读存储介质 | |
CN108919642A (zh) | 一种炉跟机协调控制系统控制器参数优化整定方法 | |
CN105276562B (zh) | 一种母管制锅炉协调的方法和系统 | |
CN112395817B (zh) | 一种基于实时数据在线计算电厂管道效率的方法 | |
CN105158007B (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 |