CN107783415B - 一种dcs数据挖掘的火电机组定滑压运行曲线计算方法 - Google Patents
一种dcs数据挖掘的火电机组定滑压运行曲线计算方法 Download PDFInfo
- Publication number
- CN107783415B CN107783415B CN201710916277.2A CN201710916277A CN107783415B CN 107783415 B CN107783415 B CN 107783415B CN 201710916277 A CN201710916277 A CN 201710916277A CN 107783415 B CN107783415 B CN 107783415B
- Authority
- CN
- China
- Prior art keywords
- pressure
- main steam
- sequence
- load
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/0205—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric not using a model or a simulator of the controlled system
- G05B13/024—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric not using a model or a simulator of the controlled system in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Control Of Turbines (AREA)
Abstract
本发明一种DCS数据挖掘的汽轮机组定滑压运行曲线计算方法:a、选取与定滑压运行曲线计算相关的采样点,获得DCS数据;b、计算调节级效率、当量主蒸汽容积流量,修正当量主蒸汽流量、修正负荷并拟合前两者之间和后两者之间的关系曲线α、β;c、假设负荷以及该负荷下主蒸汽压力序列,计算修正负荷序列;d、通过β和主蒸汽压力序列计算当量主蒸汽容积流量序列;e、通过α计算调节级效率序列,拟合主蒸汽压力序列与调节级效率序列之间的关系曲线γ,关系曲线α上的拐点对应在关系曲线γ的主蒸汽压力即为最优;f、对不同负荷下最优主蒸汽压力拟合直线即为定滑压运行曲线;本发明能够利用DCS数据辨识汽轮机组定滑压运行曲线以指导机组经济运行。
Description
技术领域
本发明属于火电机组汽轮机发电技术领域,尤其涉及一种DCS数据挖掘的汽轮机组定滑压运行曲线计算方法。
背景技术
近年来,分布式控制系统(Distributed Control System,简称DCS)在绝大多数汽轮机组上得以应用,其形成的DCS数据流蕴涵着丰富而有价值的信息有待挖掘。
根据热力循环理论,汽轮机组在低负荷下滑压运行时,由于进汽节流损失小,漏汽损失也小,使得机组的相对内效率比定压运行有一定的提高,但随着主蒸汽压力低导致循环热效率也降低,而相对内效率的增加幅度补偿了循环热效率的下降幅度时,汽轮机组的滑压压力比较经济的,因此对汽轮机组在不同负荷下经济的滑压压力进行连线作图,即得到定滑压运行曲线。
发明内容
本发明的目的在于提供一种DCS数据挖掘的汽轮机组定滑压运行曲线计算方法,在汽轮机组以顺序阀运行且高压调门流量线性度较好时,该方法能够简便地获得汽轮机组定滑压运行曲线。
为实现上述目的,本发明采用的技术方案是:
一种DCS数据挖掘的汽轮机组定滑压运行曲线计算方法,包括以下步骤:
A、从汽轮机组的DCS数据库中,选取主蒸汽压力、主蒸汽温度、调节级后蒸汽压力、调节级后蒸汽温度、凝汽器真空和机组负荷这六个测点,设定采样周期小于或等于30秒,采样时长能涵盖机组运行负荷的全范围,进入步骤B;
B、由主蒸汽压力、主蒸汽温度、调节级后蒸汽压力和调节级后蒸汽温度计算调节级效率,由主蒸汽压力、主蒸汽温度和调节级后蒸汽压力计算当量主蒸汽容积流量和修正当量主蒸汽流量,由主蒸汽压力、凝汽器真空和机组负荷计算修正负荷,进入步骤C;
C、拟合当量主蒸汽容积流量与调节级效率之间的关系曲线α,拟合修正负荷与修正当量主蒸汽流量之间的关系曲线β,进入步骤D;
D、假设负荷有85%、80%、75%、70%、65%、60%、55%和50%额定负荷,在以上每一个假设负荷下假设机组运行的主蒸汽压力序列a,由机组额定工况的主蒸汽压力和序列a计算假设修正负荷序列b,通过关系曲线β计算获得修正当量主蒸汽流量序列c,与序列a一起计算获得当量主蒸汽容积流量序列d,进入步骤E;
E、对于每一个假设负荷下的当量主蒸汽容积流量序列d,在关系曲线α上插值计算获得假设调节级效率序列e,拟合序列d与序列e之间的关系曲线γ,在该关系曲线上寻找随着主蒸汽压力降低调节级效率增长开始缓慢的点,即为该假设负荷下的最优主蒸汽压力,进入步骤F;
F、利用最小二乘法对每一个假设负荷下的最优主蒸汽压力拟合直线,在该直线到达机组额定主蒸汽压力的负荷到机组额定负荷之间绘制不同负荷但同一额定主蒸汽压力的直线,该直线连同拟合直线即为汽轮机组定滑压运行曲线。
本发明进一步的改进在于,步骤B中,调节级效率η由公式(1)计算:
式中:H(p,t)、H(p,S)、S(p,t)表示以压力p和温度t来计算焓、压力p和熵S来计算焓、压力p和温度t来计算熵的函数,函数结果在水及水蒸气物理性质表中查询获得;角标0表示主蒸汽参数,角标s表示调节级后蒸汽参数。
本发明进一步的改进在于,当量主蒸汽容积流量Qv和修正当量主蒸汽流量Qc由公式(2)计算:
式中:V(p,t)表示以压力p和温度t来计算比容,角标d表示机组额定工况参数。
本发明进一步的改进在于,修正负荷Nc由公式(3)计算:
式中:N表示实际机组负荷,VC表示凝汽器真空,Pa表示当地大气压力,Pc表示汽轮机低压缸排汽绝对压力设计值,Cs表示汽轮机低压缸排汽绝对压力对功率的修正系数,在设备厂商提供的修正曲线上查询得到。
本发明进一步的改进在于,步骤D中,序列b中元素bi与对应的序列a中元素ai满足公式(4)的计算关系:
式中:Na表示假设负荷,下标i表示统计时间段内的样本个数;
序列d中元素di由公式(5)计算:
本发明进一步的改进在于,步骤E中,最优主蒸汽压力是随着主蒸汽压力降低调节级效率增长开始缓慢的点,其寻找方法是:先在关系曲线α上寻找曲线拐点对应的Y轴坐标,在关系曲线γ上查看相同Y轴坐标下对应的主蒸汽压力值,即为最优主蒸汽压力。
本发明具有如下的优点:
1、常规采用汽轮机热耗率性能试验获得机组定、滑压试验的方法,其试验费用昂贵,本发明能够以较低费用获得机组定、滑压运行曲线。
2、常规采用汽轮机热耗率性能试验获得的机组定、滑压试验方法,需要的负荷量较大,机组在不同负荷下均要求连续稳定试验时间达到7~9个小时,总试验时间一般超过5天,特别在电网负荷量偏低的时间段进行试验,会导试验困难较大、试验周期较长。本发明能够在较短时间内,获得机组定、滑压运行曲线,一般仅需要1~2天时间,且对机组运行负荷无要求。
3、常规采用汽轮机热耗率性能试验获得的机组定、滑压试验方法,因为往往使用汽轮机性能试验规程中的简化试验方法,试验不确定度较大,特别是当采用现场永久性安装的流量测量装置作为整个汽轮机性能试验基准时,试验不确定度可到达1.5%~3.0%,不适用于精度要求较高的机组定、滑压试验。本发明采用大量机组运行统计数据,并依托理论计算及分析,得到的定、滑压曲线随机性更小,可靠性更高。
4、机组在运行过程中的老化、特性改变等因素往往会导致机组预设在控制系统内的定、滑压曲线在运行一段时间后,偏离机组最佳定、滑压曲线。本发明可在机组运行控制系统中进行扩展应用,实现机组动态定、滑压曲线计算,使机组时刻都运行在最佳定、滑压曲线下。
附图说明
图1为本发明的流程示意图。
图2为关系曲线α、β,注意X轴和Y轴的数值经过标幺。
图3为关系曲线γ,注意X轴和Y轴的数值经过标幺。
图4为实例汽轮机组定滑压运行曲线,注意X轴和Y轴的数值经过标幺。
具体实施方式
下面结合附图和实例对本发明的一种DCS数据挖掘的汽轮机组定滑压运行曲线计算方法作进一步的详细说明。
如图1所示,一种DCS数据挖掘的汽轮机组定滑压运行曲线计算方法,有以下步骤。
1、从汽轮机组的DCS数据库中,选取主蒸汽压力、主蒸汽温度、调节级后蒸汽压力、调节级后蒸汽温度、综合阀位指令、凝汽器真空、机组负荷这七个测点,设定采样周期小于或等于30秒,采样时长能涵盖机组运行负荷的全范围。
2、由主蒸汽压力、主蒸汽温度、调节级后蒸汽压力、调节级后蒸汽温度计算调节级效率η:
式中:H(p,t)、H(p,S)、S(p,t)表示以压力p和温度t来计算焓、压力p和熵S来计算焓、压力p和温度t来计算熵的函数,函数结果可以在水及水蒸气物理性质表中查询获得;角标“0”表示主蒸汽参数,角标“s”表示调节级后蒸汽参数;由主蒸汽压力、主蒸汽温度、调节级后蒸汽压力计算当量主蒸汽容积流量Qv和修正当量主蒸汽流量Qc:
式中:V(p,t)表示以压力p和温度t来计算比容,角标“d”表示机组额定工况参数;由主蒸汽压力、凝汽器真空、机组负荷计算修正负荷Nc:
式中:N表示实际机组负荷,VC表示凝汽器真空,Pa表示当地大气压力,Pc表示汽轮机低压缸排汽绝对压力设计值,Cs表示汽轮机低压缸排汽绝对压力对功率的修正系数,在设备厂商提供的修正曲线上查询得到。
3、拟合当量主蒸汽容积流量与调节级效率之间的关系曲线α,拟合修正负荷与修正当量主蒸汽流量之间的关系曲线β;如图2所示,关系曲线α拐点处的Y轴坐标值为0.9790。
4、假设负荷有85%、80%、75%、70%、65%、60%、55%、50%额定负荷,在以上每一个假设负荷下假设机组运行的主蒸汽压力序列a,其中元素以0.1MPa为步长;由机组额定工况的主蒸汽压力和序列a计算假设修正负荷序列b,序列b中元素bi与对应的序列a中元素ai满足计算公式:
式中:Na表示假设负荷,下标i表示统计时间段内的样本个数;通过关系曲线β计算获得修正当量主蒸汽流量序列c,与序列a一起计算获得当量主蒸汽容积流量序列d,其中元素di满足计算公式:
5、对于每一个假设负荷下的当量主蒸汽容积流量序列d,在关系曲线α上插值计算获得假设调节级效率序列e,拟合序列d与序列e之间的关系曲线γ,其中80%额定负荷关系曲线γ如图3所示,其余假设负荷的关系曲线γ与之类似,不再赘述;关系曲线α拐点处的Y轴坐标值为0.9790,对应在80%额定负荷关系曲线γ的主蒸汽压力为0.9418,即为该假设负荷下的最优主蒸汽压力。
6、利用最小二乘法对每一个假设负荷下的最优主蒸汽压力拟合直线,在该直线到达机组额定主蒸汽压力的负荷到机组额定负荷之间绘制不同负荷但同一额定主蒸汽压力的直线,该直线连同拟合直线即为汽轮机组定滑压运行曲线,如图4所示。
Claims (6)
1.一种DCS数据挖掘的汽轮机组定滑压运行曲线计算方法,其特征在于,包括以下步骤:
A、从汽轮机组的DCS数据库中,选取主蒸汽压力、主蒸汽温度、调节级后蒸汽压力、调节级后蒸汽温度、凝汽器真空和机组负荷这六个测点,设定采样周期小于或等于30秒,采样时长能涵盖机组运行负荷的全范围,进入步骤B;
B、由主蒸汽压力、主蒸汽温度、调节级后蒸汽压力和调节级后蒸汽温度计算调节级效率,由主蒸汽压力、主蒸汽温度和调节级后蒸汽压力计算当量主蒸汽容积流量和修正当量主蒸汽流量,由主蒸汽压力、凝汽器真空和机组负荷计算修正负荷,进入步骤C;
C、拟合当量主蒸汽容积流量与调节级效率之间的关系曲线α,拟合修正负荷与修正当量主蒸汽流量之间的关系曲线β,进入步骤D;
D、假设负荷有85%、80%、75%、70%、65%、60%、55%和50%额定负荷,在以上每一个假设负荷下假设机组运行的主蒸汽压力序列a,由机组额定工况的主蒸汽压力和序列a计算假设修正负荷序列b,通过关系曲线β计算获得修正当量主蒸汽流量序列c,序列c与序列a一起计算获得当量主蒸汽容积流量序列d,进入步骤E;
E、对于每一个假设负荷下的当量主蒸汽容积流量序列d,在关系曲线α上插值计算获得假设调节级效率序列e,拟合序列d与序列e之间的关系曲线γ,在关系曲线γ上寻找随着主蒸汽压力降低调节级效率增长开始缓慢的点,即为该假设负荷下的最优主蒸汽压力,进入步骤F;
F、利用最小二乘法对每一个假设负荷下的最优主蒸汽压力拟合直线L1,在L1上计算机组额定主蒸汽压力对应的负荷到机组额定负荷之间绘制不同负荷但同一额定主蒸汽压力的直线L2,组合L1与L2即为汽轮机组定滑压运行曲线。
6.如权利要求1所述的一种DCS数据挖掘的汽轮机组定滑压运行曲线计算方法,其特征在于,步骤E中,最优主蒸汽压力是随着主蒸汽压力降低调节级效率增长开始缓慢的点,其寻找方法是:先在关系曲线α上寻找曲线拐点对应的Y轴坐标,在关系曲线γ上查看相同Y轴坐标下对应的主蒸汽压力值,即为最优主蒸汽压力。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710916277.2A CN107783415B (zh) | 2017-09-30 | 2017-09-30 | 一种dcs数据挖掘的火电机组定滑压运行曲线计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710916277.2A CN107783415B (zh) | 2017-09-30 | 2017-09-30 | 一种dcs数据挖掘的火电机组定滑压运行曲线计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107783415A CN107783415A (zh) | 2018-03-09 |
CN107783415B true CN107783415B (zh) | 2020-09-01 |
Family
ID=61434013
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710916277.2A Active CN107783415B (zh) | 2017-09-30 | 2017-09-30 | 一种dcs数据挖掘的火电机组定滑压运行曲线计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107783415B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108716424B (zh) * | 2018-03-27 | 2019-11-05 | 华电电力科学研究院有限公司 | 基于背压实时修正的汽轮机滑压曲线的计算方法 |
CN108930565A (zh) * | 2018-07-30 | 2018-12-04 | 深圳众诚联合能源科技有限公司 | 一种汽轮机组调峰优化方法及系统 |
CN109057885B (zh) * | 2018-08-09 | 2020-12-04 | 上海电气电站设备有限公司 | 一种汽轮机组启动过程中的升负荷方法 |
CN109869198B (zh) * | 2019-01-23 | 2021-06-29 | 西安热工研究院有限公司 | 一种基于阀位修正的汽轮机全寿命高压缸效率监测方法 |
CN110118107A (zh) * | 2019-06-03 | 2019-08-13 | 中国神华能源股份有限公司 | 汽轮机组的配汽方法和配汽装置 |
CN110469372B (zh) * | 2019-08-21 | 2022-04-12 | 西安热工研究院有限公司 | 一种汽轮机热力性能试验不确定度控制方法 |
CN111027748B (zh) * | 2019-11-11 | 2022-04-26 | 国网湖南省电力有限公司 | 火力发电机组滑压运行曲线的多元函数实时优化方法、系统及介质 |
CN111694275B (zh) * | 2020-06-09 | 2022-06-03 | 烟台龙源电力技术股份有限公司 | 一种设备控制方法、装置及电子设备 |
CN112302736B (zh) * | 2020-10-31 | 2023-03-24 | 河南九域恩湃电力技术有限公司 | 一种基于监视段参数的汽轮机主蒸汽压力控制方法 |
CN112765227B (zh) * | 2020-12-09 | 2023-02-24 | 华能陕西发电有限公司 | 一种基于dcs数据挖掘的燃煤发电机组背压特性曲线获取系统及方法 |
CN114943127B (zh) * | 2022-06-28 | 2023-07-11 | 西安热工研究院有限公司 | 一种汽轮机顺序阀状态修正到单阀状态下的计算方法 |
CN116500898B (zh) * | 2023-05-11 | 2024-03-19 | 华电国际电力股份有限公司莱城发电厂 | 基于特征流量辨识的火电机组agc负荷控制系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009098007A2 (de) * | 2008-02-04 | 2009-08-13 | Voith Patent Gmbh | Rotorblattgestaltung für eine wellsturbine |
CN102661176A (zh) * | 2012-06-04 | 2012-09-12 | 哈尔滨工业大学 | 一种以主蒸汽流量作为调度变量获取汽轮机滑压曲线方法 |
CN103850726A (zh) * | 2014-03-28 | 2014-06-11 | 苏州华瑞能泰发电技术有限公司 | 一种快速确定汽轮机定滑压优化曲线的方法 |
CN104915545A (zh) * | 2015-05-21 | 2015-09-16 | 华电电力科学研究院 | 一种火力发电机组定滑压运行方法 |
CN106932201A (zh) * | 2017-05-02 | 2017-07-07 | 哈尔滨工业大学 | 考虑阀门节流损失的火电机组纯凝工况滑压曲线确定方法 |
-
2017
- 2017-09-30 CN CN201710916277.2A patent/CN107783415B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009098007A2 (de) * | 2008-02-04 | 2009-08-13 | Voith Patent Gmbh | Rotorblattgestaltung für eine wellsturbine |
CN102661176A (zh) * | 2012-06-04 | 2012-09-12 | 哈尔滨工业大学 | 一种以主蒸汽流量作为调度变量获取汽轮机滑压曲线方法 |
CN103850726A (zh) * | 2014-03-28 | 2014-06-11 | 苏州华瑞能泰发电技术有限公司 | 一种快速确定汽轮机定滑压优化曲线的方法 |
CN104915545A (zh) * | 2015-05-21 | 2015-09-16 | 华电电力科学研究院 | 一种火力发电机组定滑压运行方法 |
CN106932201A (zh) * | 2017-05-02 | 2017-07-07 | 哈尔滨工业大学 | 考虑阀门节流损失的火电机组纯凝工况滑压曲线确定方法 |
Non-Patent Citations (3)
Title |
---|
Sliding pressure operation of large conventional steam power units;GEORGE DARIE 等;《Proceedings of the 5th IASME/WSEAS Int. Conference on Heat Transfer, Thermal Engineering and Environment》;20140528;全文 * |
实用机组定滑压运行曲线试验研究;邵峰 等;《动力工程学报》;20130831(第8期);全文 * |
汽轮机定滑压曲线在实际应用中的计算寻优;文乐 等;《动力工程学报》;20180131(第1期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN107783415A (zh) | 2018-03-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107783415B (zh) | 一种dcs数据挖掘的火电机组定滑压运行曲线计算方法 | |
CN103850726B (zh) | 一种快速确定汽轮机定滑压优化曲线的方法 | |
CN111352400B (zh) | 一种供热机组滑压运行优化方法及滑压运行控制系统 | |
CN111652418B (zh) | 一种火电机组动态精细化复合参数滑压曲线生成方法 | |
CN105275508A (zh) | 一种基于功率值计算的汽轮机流量曲线辨识及优化方法 | |
CN102767529B (zh) | 一种轴流压缩机节能效果的计算方法 | |
CN110608070A (zh) | 汽轮机最小安全流量控制方法 | |
CN105317476A (zh) | 一种基于特征通流面积的汽轮机流量曲线辨识及优化方法 | |
CN105512800A (zh) | 根据以热定电确定供热机组调峰调度的方法 | |
CN103235512B (zh) | 一种发电机组运行的方法 | |
CN106250617B (zh) | 一种火电机组双背压冷端系统变频泵循环水流量优化方法 | |
CN109086551B (zh) | 一种适用于空冷供热机组变工况的汽轮机最优滑压曲线获取方法 | |
CN105240058A (zh) | 一种基于喷嘴流量计算的汽轮机流量曲线辨识及优化方法 | |
CN110702175B (zh) | 一种火电厂汽轮机主蒸汽流量在线软测量装置及方法 | |
CN112100862A (zh) | 基于变比热容的新型涡轮等熵效率快速求解方法 | |
CN109869198B (zh) | 一种基于阀位修正的汽轮机全寿命高压缸效率监测方法 | |
CN104408540A (zh) | 一种汽轮机改造节能量的评估方法 | |
CN114017200B (zh) | 一种涡扇发动机进口总压重构方法及装置 | |
CN113218109B (zh) | 一种深度余热回收智能调控装置及方法 | |
CN115292946A (zh) | 一种基于变比热计算的高压涡轮效率评估方法及装置 | |
CN111507011B (zh) | 供热抽汽对汽轮机滑压运行影响量的修正方法 | |
CN106202826B (zh) | 一种火电机组双背压冷端系统双速泵循环水流量优化方法 | |
CN110298492B (zh) | 一种凝汽器抽真空系统节能改造收益的评价方法 | |
KR102234570B1 (ko) | 화력발전 플랜트의 효율예측 정확성 향상을 위한 연료의 발열량 보정방법 | |
CN108345573B (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 |