CN105697074A - 一种汽轮机组高压调门流量优化装置及其优化方法 - Google Patents
一种汽轮机组高压调门流量优化装置及其优化方法 Download PDFInfo
- Publication number
- CN105697074A CN105697074A CN201610073044.6A CN201610073044A CN105697074A CN 105697074 A CN105697074 A CN 105697074A CN 201610073044 A CN201610073044 A CN 201610073044A CN 105697074 A CN105697074 A CN 105697074A
- Authority
- CN
- China
- Prior art keywords
- valve
- flow
- fdem
- curve
- flow characteristic
- 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
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 64
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000004088 simulation Methods 0.000 claims abstract description 98
- 230000001105 regulatory effect Effects 0.000 claims description 89
- 239000013598 vector Substances 0.000 claims description 72
- 238000012937 correction Methods 0.000 claims description 43
- 238000005259 measurement Methods 0.000 claims description 28
- 238000002474 experimental method Methods 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000012360 testing method Methods 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 7
- 238000012986 modification Methods 0.000 claims description 7
- 230000004048 modification Effects 0.000 claims description 7
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 230000001276 controlling effect Effects 0.000 claims description 3
- 230000003828 downregulation Effects 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims 1
- 239000002994 raw material Substances 0.000 claims 1
- 230000008569 process Effects 0.000 abstract description 8
- 230000007547 defect Effects 0.000 abstract description 2
- 230000033228 biological regulation Effects 0.000 description 3
- 238000012546 transfer Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000012887 quadratic function Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000002791 soaking Methods 0.000 description 1
- 238000003466 welding Methods 0.000 description 1
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
- F01D—NON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
- F01D17/00—Regulating or controlling by varying flow
- F01D17/10—Final actuators
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Details Of Valves (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种汽轮机组高压调门流量优化装置,包括指令控制模块,用于发出各个调门的开度信号;执行模块,用于接收指令并控制各个调门的开度变化;数据采集模块,用于采集汽轮机的调节级压力、主汽压等实时参数;数据仿真优化模块,用于根据调门的开度以及数据采集模块采集的汽轮机参数进行数据仿真和以及对流量特性曲线进行优化。本发明还公开了一种汽轮机组高压调门流量的优化方法。本发明能够改进现有技术的不足,消除切换过程中负荷扰动现象,进而保障机组阀切换状态下安全、稳定、经济运行。
Description
技术领域
本发明涉及发电汽轮机组技术领域,尤其是一种汽轮机组高压调门流量优化装置及其优化方法。
背景技术
目前,随着电力建设的快速进行,电网结构越发复杂,为了保障电网的安全稳定运行,必须使发电机组具有快速响应电网自动发电控制指令(AGC,AutomaticGenerationControl)变化和一次调频需求的调节能力。而实现方法主要是控制汽轮机高压调门动作,快速调节汽轮机的进汽量,进而实现机组负荷的快速响应。由此可见,数字电液控制系统(DEH,DigitalElectricHydraulicControlSystem)及其高压调门整体流量特性的优劣将直接关系到发电机的控制品质、调节性能、经济效益等重要问题。
汽轮机高压阀门运行方式分为单阀运行方式和顺序阀(或多阀)运行方式。单阀运行时,汽轮机总的流量信号平均分配到各个高压调节门;顺序阀控制时,流入汽轮机的蒸汽流量是各阀门流量的总和,它将按顺序依次分配到各个高压调节门,而阀门流量特性曲线就是流量与阀门开度的转换函数,如果单阀运行方式下整体流量特性曲线与顺序阀运行方式下整体流量特性曲线啮合度较差的情况下,则阀门切换过程中负荷扰动就较大,严重时会影响机组响应电网AGC指令变化需求。单-顺序阀切换目的是为了提高机组的经济性和稳定性,其实质是实现节流调节与喷嘴调节的无扰切换,解决变负荷过程中均热要求与部分负荷经济性的矛盾。
因此,需要发明一种针对单-顺序阀整体流量特性曲线啮合度较差的情况,单独优化顺序阀方式下阀门管理相关参数无法消除切换过程中负荷扰动情况的优化整定策略。使机组单-顺序阀切换时无限接近无扰切换,消除切换过程中负荷扰动现象,进而保障机组阀切换状态下安全、稳定、经济运行。
发明内容
本发明要解决的技术问题是提供一种汽轮机组高压调门流量优化装置及其优化方法,能够解决现有技术的不足,消除切换过程中负荷扰动现象,进而保障机组阀切换状态下安全、稳定、经济运行。
为解决上述技术问题,本发明所采取的技术方案如下。
一种汽轮机组高压调门流量优化装置,包括,
指令控制模块,用于发出各个调门的开度信号;
执行模块,用于接收指令并控制各个调门的开度变化;
数据采集模块,用于采集汽轮机的调节级压力、主汽压等实时参数;
数据仿真优化模块,用于根据调门的开度以及数据采集模块采集的汽轮机参数进行数据仿真和以及对流量特性曲线进行优化。
一种使用上述的汽轮机组高压调门流量优化装置的优化方法,包括以下步骤:
A、建立顺序阀方式下整体流量特性仿真模型;
指令控制模块和执行模块对调门开度进行调整,数据采集模块实测单个调门流量特性曲线和调门整体流量特性曲线,数据仿真优化模块结合阀门控制程序建立顺序阀方式下整体流量特性仿真模型;
B、建立单阀方式下整体流量特性仿真模型;
指令控制模块和执行模块对调门开度进行调整,数据采集模块实测单阀整体流量特性曲线,数据仿真优化模块结合阀门控制程序建立单阀方式下整体流量特性仿真模型;
C、单-顺序阀切换参数整定及优化;
数据仿真优化模块根据单阀整体流量特性仿真模型,调整最佳阀位系数,优化调整单阀整体流量特性曲线的斜率,使之向多阀整体流量特性曲线;
数据仿真优化模块根据多阀整体流量特性仿真模型,调整流量比例因子,优化调整多阀整体流量特性曲线的斜率,使之向多阀整体流量特性曲线;
数据仿真优化模块通过优化重叠度函数,提高多阀整体流量特性曲线的重叠区线性度;
数据仿真优化模块通过优化流量开度修正函数,提高多阀整体流量特性曲线的非重叠区线性度。
作为优选,步骤A中,单个调门流量特性实测曲线通过如下步骤进行计算,
通过调门流量特性实验得到调门开度反馈、调节级压力、主汽压力、向量的维度,其中表示调门全关,表示调门全开,则调节级压力和主汽压力比值
与调门流量成正比,即分别对应开度为0和100时的调门流量百分比0%和100%,将归一化0-100之间的数,得到调门流量百分比
其中以向量为横坐标、向量为纵坐标构成的曲线即为第i个调门的流量特性实测曲线,该曲线用线性分段函数表示,记作
其中,为调门开度反馈,为调节级压力,为主汽压力,为向量的维度。
作为优选,步骤A中,多阀方式下整体流量特性实测曲线通过如下步骤进行计算,
进行多阀方式下整体流量特性曲线测试实验,实验测量获得数据向量:调门总流量指令、调节级压力、主汽压力、向量的维度,且FDEMseq(1)>0,FDEMseq(nseq)=100,表示顺序阀方式下调门全部开启;调节级压力和主汽压力比值Hseq(k)与调门流量FLseq(k)成正比
则Hseq(nseq)对应调门全开时的调门整体流量百分比100%。由于调门总流量指令为0时,调门全关,整体流量百分比必为0%,对Hseq(k)进行归一化,得到调门流量百分比FLseq(k)
且FLseq(nseq)=100,以向量FDEMseq为横坐标、向量FLseq为纵坐标构成的曲线即为顺序阀方式调门整体流量特性实测曲线;
其中,FDEMseq为调门总流量指令、PIseq为调节级压力、TPseq为主汽压力、nseq×1为向量的维度。
作为优选,步骤A中,顺序阀方式下整体流量特性仿真的步骤如下,
流量总指令FDEM经背压修正函数f1(x)、流量比例偏置因子、重叠度函数修正和流量开度修正函数修正分配到各个调门上,获得各调门的开度指令,根据实测的各调门流量特性曲线计算各调门自身的流量百分比。为线性分段函数;
设当前流量总指令为FDEM,FDEM为向量,取值范围为0-100,每隔0.1取一个值,经背压修正函数f1(x)修正得到流量指令FL1
FL1=f1(FDEM)(5)
表示对FDEM中的每个元素做函数f1(x)映射,得到的元素构成数据向量FL1,将流量指令FL1经过流量比例偏置因子得到分配到各个调门上的流量指令
其中为分配到调门i的流量指令,ki、bi为调门i的比例偏置因子,调门流量指令经各自的重叠度函数修正得到新的流量指令向量
流量指令经过各调门流量开度修正函数得到各调门开度指令向量OPENCMDi,为保证单阀方式下各调门开度一致,各调门流量开度修正函数相同
在获得各调门开度OPENCMDi的基础上,根据调门i的实测流量特性曲线计算各调门流量百分比向量
调门实测流量特性曲线是在其他调门全开时测试获得,顺序阀方式下调门依次开启,因此,整体流量特性中各调门流量百分比和实测时流量百分比之间有一定线性关系,计算顺序阀方式下仿真整体流量特性FLseq
其中,a、b、c、d为仿真参数,表示调门仿真流量在整体流量中的占比,为式(9)中各调门仿真流量,FL′seq为整体流量,以向量FDEM为横坐标、向量FLseq为纵坐标构成的曲线即为顺序阀方式调门整体流量特性仿真曲线;
将顺序阀方式整体流量特性实测曲线和仿真曲线绘制在同一坐标系中,调整仿真参数a、b、c、d使两条曲线拟合,则公式(10)即为顺序阀方式下整体流量特性仿真模型。
作为优选,所述步骤B中,
对单阀整体流量特性实测曲线计算步骤如下,
进行单阀方式下整体流量特性曲线测试实验,实验测量获得数据向量:调门总流量指令FDEMsig,调门开度指令OPENsig,调节级压力PIsig,主汽压力TPsig,向量的维度为nsig×1,且有FDEMsig(1)>0,FDEMsig(nsig)=100,表示单阀方式下调门全部开启,使用步骤A中多阀方式下整体流量特性实测曲线的计算步骤,获得单阀方式调门整体流量百分比向量FLsig,以向量FDEMsig为横坐标、向量FLsig为纵坐标构成的曲线即为单阀方式整体流量特性实测曲线;
此外,单阀方式下各调门开度指令相同,以单阀方式下调门开度OPENsig为横坐标,整体流量FLsig为纵坐标,得到整体流量关于调门开度的曲线,用线性分段函数h(x)表示;
单阀整体流量特性仿真的步骤如下,
单阀方式下流量总指令FDEM经过背压修正函数f1(x)修正,乘以单阀方式下最佳阀位系数kVLVPT,再经过流量开度修正函数修正,得到阀门开度指令;
设当前流量总指令为FDEM,FDEM为向量,取值范围为0-100,每隔0.1取一个值,FDEM经背压修正函数f1(x)修正得到流量指令FL1,流量指令FL1乘以最佳阀位系数kVLVPT,得到流量指令FL′1,再经流量开度函数修正得到各调门开度指令OPEN′sig
FL′1=FL1*kVLVPT(11)
根据单阀方式下调门开度与整体流量的函数关系h(x),计算单阀方式下整体仿真流量FLsig
FL′sig=h(OPEN′sig)(13)
以向量FDEM为横坐标、向量FLsig为纵坐标构成的曲线即为单阀方式调门整体流量特性仿真曲线,公式(13)为单阀方式调门整体流量仿真模型。
作为优选,步骤C中,
单阀整体流量特性曲线斜率优化的步骤如下,
在同一坐标轴下绘制单阀和多阀方式仿真流量特性曲线,调整kVLVPT使单阀仿真流量特性曲线向多阀方式仿真流量特性曲线靠近,从而减小单多阀切换时的扰动。
作为优选,步骤C中,
多阀整体流量特性曲线优化的步骤如下,
1)背压修正函数调整
单阀方式下当总流量指令FDEM=100时,经背压函数和最佳阀位系数后的流量指令FL′1满足FL′1=100,单阀整体流量特性曲线斜率优化的步骤中调整最佳阀位系数kVLVPT后,背压修正函数根据式(5)和式(11)调整;
2)流量比例偏置因子优化
流量比例偏置因子ki、bi(i=1,2,…,6)控制相应调节阀在总流量指令FDEM分别为某两个值时开始开启和全开;在不考虑重叠度函数的情况下,流量指令分别取值为0和100时对应的FDEM值即为调节阀i的开启位置FDEM0和全开位置FDEM100,根据式(5)和式(6)反算得出
f1(FDEM0)*ki+bi=0
f1(FDEM100)*ki+bi=100(14)
保持偏置因子bi的值不变,增大比例偏置因子ki,调节阀的开启位置FDEM0和全开位置FDEM100均减小,且调节阀所占的FDEM区间缩小,整理流量特性曲线的斜率增大;相反减小ki,整理流量特性曲线的斜率减小;
在同一坐标轴下绘制单阀和多阀方式仿真流量特性曲线,调整流量比例因子ki使多阀仿真流量特性曲线向单阀方式仿真流量特性曲线靠近,从而减小单多阀切换时的扰动。
作为优选,步骤C中,
重叠度函数的优化步骤为,
重叠度函数控制相邻两个调门重叠区域的大小;根据顺序阀调门整体流量仿真曲线中重叠区域流量特性的线性度,调整重叠度函数使重叠区呈现线性关系。
作为优选,步骤C中,
调门流量开度修正函数的优化步骤为,
调门流量开度修正函数是将各调门的流量指令折算成调门开度,为保证单阀方式下各调门的开度保持一致,设置各调门的流量开度修正函数一样;根据多阀方式下单个调门流量特性实测曲线,做一条流量开度曲线使之与各调门的实测曲线的距离之和最小,将此曲线作为优化后的调门流量开度函数。
采用上述技术方案所带来的有益效果在于:本发明可以使机组单-顺序阀切换时无限接近无扰切换,消除切换过程中负荷扰动现象,进而保障机组阀切换状态下安全、稳定、经济运行。
附图说明
图1是本发明一个具体实施方式的硬件示意图。
图2是本发明一个具体实施方式中DEH阀门控制过程的示意图。
图3是高压调门流量特性优化方法框图
图4是优化前单多阀整体流量特性曲线。
图5是优化后单多阀整体流量特性曲线。
图中:1、指令控制模块;2、执行模块;3、数据采集模块;4、数据仿真优化模块。
具体实施方式
本发明中使用到的标准零件均可以从市场上购买,异形件根据说明书的和附图的记载均可以进行订制,各个零件的具体连接方式均采用现有技术中成熟的螺栓、铆钉、焊接、粘贴等常规手段,在此不再详述。
参照图1-5,本发明一个具体实施方式包括指令控制模块1,用于发出各个调门的开度信号;
执行模块2,用于接收指令并控制各个调门的开度变化;
数据采集模块3,用于采集汽轮机的调节级压力、主汽压等实时参数;
数据仿真优化模块4,用于根据调门的开度以及数据采集模块3采集的汽轮机参数进行数据仿真和以及对流量特性曲线进行优化。
一种使用上述汽轮机组高压调门流量优化装置的优化方法,包括以下步骤:
A、建立顺序阀方式下整体流量特性仿真模型;
指令控制模块1和执行模块2对调门开度进行调整,数据采集模块3实测单个调门流量特性曲线和调门整体流量特性曲线,数据仿真优化模块4结合阀门控制程序建立顺序阀方式下整体流量特性仿真模型;
B、建立单阀方式下整体流量特性仿真模型;
指令控制模块1和执行模块2对调门开度进行调整,数据采集模块3实测单阀整体流量特性曲线,数据仿真优化模块4结合阀门控制程序建立单阀方式下整体流量特性仿真模型;
C、单-顺序阀切换参数整定及优化;
数据仿真优化模块4根据单阀整体流量特性仿真模型,调整最佳阀位系数,优化调整单阀整体流量特性曲线的斜率,使之向多阀整体流量特性曲线;
数据仿真优化模块4根据多阀整体流量特性仿真模型,调整流量比例因子,优化调整多阀整体流量特性曲线的斜率,使之向多阀整体流量特性曲线;
数据仿真优化模块4通过优化重叠度函数,提高多阀整体流量特性曲线的重叠区线性度;
数据仿真优化模块4通过优化流量开度修正函数,提高多阀整体流量特性曲线的非重叠区线性度。
步骤A中,单个调门流量特性实测曲线通过如下步骤进行计算,
通过调门流量特性实验得到调门开度反馈、调节级压力、主汽压力、向量的维度,其中表示调门全关,表示调门全开,则调节级压力和主汽压力比值
与调门流量成正比,即分别对应开度为0和100时的调门流量百分比0%和100%,将归一化0-100之间的数,得到调门流量百分比
其中以向量为横坐标、向量为纵坐标构成的曲线即为第i个调门的流量特性实测曲线,该曲线用线性分段函数表示,记作
其中,为调门开度反馈,为调节级压力,为主汽压力,为向量的维度。
步骤A中,多阀方式下整体流量特性实测曲线通过如下步骤进行计算,
进行多阀方式下整体流量特性曲线测试实验,实验测量获得数据向量:调门总流量指令、调节级压力、主汽压力、向量的维度,且FDEMseq(1)>0,FDEMseq(nseq)=100,表示顺序阀方式下调门全部开启;调节级压力和主汽压力比值Hseq(k)与调门流量FLseq(k)成正比
则Hseq(nseq)对应调门全开时的调门整体流量百分比100%。由于调门总流量指令为0时,调门全关,整体流量百分比必为0%,对Hseq(k)进行归一化,得到调门流量百分比FLseq(k)
且FLseq(nseq)=100,以向量FDEMseq为横坐标、向量FLseq为纵坐标构成的曲线即为顺序阀方式调门整体流量特性实测曲线;
其中,FDEMseq为调门总流量指令、PIseq为调节级压力、TPseq为主汽压力、nseq×1为向量的维度。
步骤A中,顺序阀方式下整体流量特性仿真的步骤如下,
流量总指令FDEM经背压修正函数f1(x)、流量比例偏置因子、重叠度函数修正和流量开度修正函数修正分配到各个调门上,获得各调门的开度指令,根据实测的各调门流量特性曲线计算得到各调门自身的流量百分比。f1(x)、为线性分段函数;
设当前流量总指令为FDEM,FDEM为向量,取值范围为0-100,每隔0.1取一个值,经背压修正函数f1(x)修正得到流量指令FL1
FL1=f1(FDEM)(5)
表示对FDEM中的每个元素做函数f1(x)映射,得到的元素构成数据向量FL1,将流量指令FL1经过流量比例偏置因子得到分配到各个调门上的流量指令
其中为分配到调门i的流量指令,ki、bi为调门i的比例偏置因子,调门流量指令经各自的重叠度函数修正得到新的流量指令向量
流量指令经过各调门流量开度修正函数得到各调门开度指令向量OPENCMDi,为保证单阀方式下各调门开度一致,各调门流量开度修正函数相同
在获得各调门开度OPENCMDi的基础上,根据调门i的实测流量特性曲线计算各调门流量百分比向量
调门实测流量特性曲线是在其他调门全开时测试获得,顺序阀方式下调门依次开启,因此,整体流量特性中各调门流量百分比和实验时流量百分比之间有一定线性关系,计算顺序阀方式下仿真整体流量特性FL′seq
其中,a、b、c、d为仿真参数,表示调门仿真流量在整体流量中的占比,为式(9)中各调门仿真流量,FL′seq为整体流量,以向量FDEM为横坐标、向量FL′seq为纵坐标构成的曲线即为顺序阀方式调门整体流量特性仿真曲线;
将顺序阀方式整体流量特性实测曲线和仿真曲线绘制在同一坐标系中,调整仿真参数a、b、c、d使两条曲线拟合,则公式(10)即为顺序阀方式下整体流量特性仿真模型。
步骤B中,
对单阀整体流量特性实测曲线计算步骤如下,
进行单阀方式下整体流量特性曲线测试实验,实验测量获得数据向量:调门总流量指令FDEMsig,调门开度指令OPENsig,调节级压力PIsig,主汽压力TPsig,向量的维度为nsig×1,且有FDEMsig(1)>0,FDEMsig(nsig)=100,表示单阀方式下调门全部开启,使用步骤A中多阀方式下整体流量特性实测曲线的计算步骤,获得单阀方式调门整体流量百分比向量FLsig,以向量FDEMsig为横坐标、向量FLsig为纵坐标构成的曲线即为单阀方式整体流量特性实测曲线;
此外,单阀方式下各调门开度指令相同,以单阀方式下调门开度OPENsig为横坐标,整体流量FLsig为纵坐标,得到整体流量关于调门开度的曲线,用线性分段函数h(x)表示;
单阀整体流量特性仿真的步骤如下,
单阀方式下流量总指令FDEM经过背压修正函数f1(x)修正,乘以单阀方式下最佳阀位系数kVLVPT,再经过流量开度修正函数修正,得到阀门开度指令;
设当前流量总指令为FDEM,FDEM为向量,取值范围为0-100,每隔0.1取一个值,FDEM经背压修正函数f1(x)修正得到流量指令FL1,流量指令FL1乘以最佳阀位系数kVLVPT,得到流量指令FL′1,再经流量开度函数修正得到各调门开度指令OPEN′sig
FL′1=FL1*kVLVPT(11)
根据单阀方式下调门开度与整体流量的函数关系h(x),计算单阀方式下整体仿真流量FL′sig
FL′sig=h(OPEN′sig)(13)
以向量FDEM为横坐标、向量FLsig为纵坐标构成的曲线即为单阀方式调门整体流量特性仿真曲线,公式(13)为单阀方式调门整体流量仿真模型。
步骤C中,
单阀整体流量特性曲线斜率优化的步骤如下,
在同一坐标轴下绘制单阀和多阀方式仿真流量特性曲线,调整kVLVPT使单阀仿真流量特性曲线向多阀方式仿真流量特性曲线靠近,从而减小单多阀切换时的扰动。
步骤C中,
多阀整体流量特性曲线优化的步骤如下,
1)背压修正函数调整
单阀方式下当总流量指令FDEM=100时,经背压函数和最佳阀位系数后的流量指令FL′1满足FL′1=100,单阀整体流量特性曲线斜率优化的步骤中调整最佳阀位系数kVLVPT后,背压修正函数根据式(5)和式(11)调整;
2)流量比例偏置因子优化
流量比例偏置因子ki、bi(i=1,2,…,6)控制相应调节阀在总流量指令FDEM分别为某两个值时开始开启和全开;在不考虑重叠度函数的情况下,流量指令分别取值为0和100时对应的FDEM值即为调节阀i的开启位置FDEM0和全开位置FDEM100,根据式(5)和式(6)反算得出
f1(FDEM0)*ki+bi=0
f1(FDEM100)*ki+bi=100(14)
保持偏置因子bi的值不变,增大比例偏置因子ki,调节阀的开启位置FDEM0和全开位置FDEM100均减小,且调节阀所占的FDEM区间缩小,整理流量特性曲线的斜率增大;相反减小ki,整理流量特性曲线的斜率减小;
在同一坐标轴下绘制单阀和多阀方式仿真流量特性曲线,调整流量比例因子ki使多阀仿真流量特性曲线向单阀方式仿真流量特性曲线靠近,从而减小单多阀切换时的扰动。
步骤C中,
重叠度函数的优化步骤为,
重叠度函数控制相邻两个调门重叠区域的大小;根据顺序阀调门整体流量仿真曲线中重叠区域流量特性的线性度,调整重叠度函数使重叠区呈现线性关系。
另外,在优化重叠度函数时,还可以利用调门流量百分比FLseq(k)与调门总流量指令FDEMseq的比例变化率,对重叠度函数设置与上述比例变化率呈线性关系的斜率参数,然后通过调门总流量指令DDEMseq在全开和关闭两个极点上的重叠度函数数值作为重叠度函数的端点值,形成重叠度函数的预设值。最后使用重叠区流量特性采样点对重叠度函数进行比例常数校准,得到优化后的重叠度函数
步骤C中,
调门流量开度修正函数的优化步骤为,
调门流量开度修正函数是将各调门的流量指令折算成调门开度,为保证单阀方式下各调门的开度保持一致,设置各调门的流量开度修正函数一样;根据多阀方式下单个调门流量特性实测曲线,做一条流量开度曲线使之与各调门的实测曲线的距离之和最小,将此曲线作为优化后的调门流量开度修正函数。
另外,在优化调门流量开度修正函数时,每段线性函数的斜率通过传递函数确定。传递函数为二次函数,传递函数的二次项系数、一次项系数和常数系数与用调节级压力PIseq和主汽压力TPseq的拟合函数最高次系数线性相关。其中二次项系数与四次拟合函数的最高次系数线性相关,一次项系数与二次拟合函数的最高次系数线性相关,常数系数与一次拟合函数的常数项线性相关。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (10)
1.一种汽轮机组高压调门流量优化装置,其特征在于:包括,
指令控制模块(1),用于发出各个调门的开度信号;
执行模块(2),用于接收指令并控制各个调门的开度变化;
数据采集模块(3),用于采集汽轮机的调节级压力、主汽压等实时参数;
数据仿真优化模块(4),用于根据调门的开度以及数据采集模块(3)采集的汽轮机参数进行数据仿真和以及对流量特性曲线进行优化。
2.一种使用权利要求1所述的汽轮机组高压调门流量优化装置的优化方法,其特征在于包括以下步骤:
A、建立顺序阀方式下整体流量特性仿真模型;
指令控制模块(1)和执行模块(2)对调门开度进行调整,数据采集模块(3)实测单个调门流量特性曲线和调门整体流量特性曲线,数据仿真优化模块(4)结合阀门控制程序建立顺序阀方式下整体流量特性仿真模型;
B、建立单阀方式下整体流量特性仿真模型;
指令控制模块(1)和执行模块(2)对调门开度进行调整,数据采集模块(3)实测单阀整体流量特性曲线,数据仿真优化模块(4)结合阀门控制程序建立单阀方式下整体流量特性仿真模型;
C、单-顺序阀切换参数整定及优化;
数据仿真优化模块(4)根据单阀整体流量特性仿真模型,调整最佳阀位系数,优化调整单阀整体流量特性曲线的斜率,使之向多阀整体流量特性曲线;
数据仿真优化模块(4)根据多阀整体流量特性仿真模型,调整流量比例因子,优化调整多阀整体流量特性曲线的斜率,使之向多阀整体流量特性曲线;
数据仿真优化模块(4)通过优化重叠度函数,提高多阀整体流量特性曲线的重叠区线性度;
数据仿真优化模块(4)通过优化流量开度修正函数,提高多阀整体流量特性曲线的非重叠区线性度。
3.根据权利要求2所述的优化方法,其特征在于:步骤A中,单个调门流量特性实测曲线通过如下步骤进行计算,
通过调门流量特性实验得到调门开度反馈、调节级压力、主汽压力、向量的维度,其中表示调门全关,表示调门全开,则调节级压力和主汽压力比值
与调门流量成正比,即分别对应开度为0和100时的调门流量百分比0%和100%,将归一化0-100之间的数,得到调门流量百分比
其中 以向量为横坐标、向量为纵坐标构成的曲线即为第i个调门的流量特性实测曲线,该曲线用线性分段函数表示,记作
其中,为调门开度反馈,为调节级压力,为主汽压力,为向量的维度。
4.根据权利要求3所述的优化方法,其特征在于:步骤A中,多阀方式下整体流量特性实测曲线通过如下步骤进行计算,
进行多阀方式下整体流量特性曲线测试实验,实验测量获得数据向量:调门总流量指令、调节级压力、主汽压力、向量的维度,且FDEMseq(1)>0,FDEMseq(nseq)=100,表示顺序阀方式下调门全部开启;调节级压力和主汽压力比值Hseq(k)与调门流量FLseq(k)成正比
Hseq(k)∝FLscq(k),k=1,2,…,nseq
则Hseq(nseq)对应调门全开时的调门整体流量百分比100%。由于调门总流量指令为0时,调门全关,整体流量百分比必为0%,对Hseq(k)进行归一化,得到调门流量百分比FLseq(k)
且FLseq(nseq)=100,以向量FDEMseq为横坐标、向量FLseq为纵坐标构成的曲线即为顺序阀方式调门整体流量特性实测曲线;
其中,FDEMseq为调门总流量指令、PIseq为调节级压力、TPseq为主汽压力、nseq×1为向量的维度。
5.根据权利要求4所述的优化方法,其特征在于:步骤A中,顺序阀方式下整体流量特性仿真的步骤如下,
流量总指令FDEM经背压修正函数f1(x)、流量比例偏置因子、重叠度函数修正和流量开度修正函数修正分配到各个调门上,获得各调门的开度指令,根据实测的各调门流量特性曲线计算各调门自身的流量百分比。f1(x)、i=1,2,…,6为线性分段函数;
设当前流量总指令为FDEM,FDEM为向量,取值范围为0-100,每隔0.1取一个值,经背压修正函数f1(x)修正得到流量指令FL1
FL1=f1(FDEM)(5)
表示对FDEM中的每个元素做函数f1(x)映射,得到的元素构成数据向量FL1,将流量指令FL1经过流量比例偏置因子得到分配到各个调门上的流量指令
其中为分配到调门i的流量指令,ki、bi为调门i的比例偏置因子,调门流量指令经各自的重叠度函数修正得到新的流量指令向量
流量指令经过各调门流量开度修正函数得到各调门开度指令向量OPENCMDi,为保证单阀方式下各调门开度一致,各调门流量开度修正函数相同
在获得各调门开度OPENCMDi的基础上,根据调门i的实测流量特性曲线计算各调门流量百分比向量
调门实测流量特性曲线是在其他调门全开时测试获得,顺序阀方式下调门依次开启,因此,整体流量特性中各调门流量百分比和实测时流量百分比之间有一定线性关系,计算顺序阀方式下仿真整体流量特性
其中,a、b、c、d为仿真参数,表示调门仿真流量在整体流量中的占比,为式(9)中各调门仿真流量,为整体流量,以向量FDEM为横坐标、向量为纵坐标构成的曲线即为顺序阀方式调门整体流量特性仿真曲线;
将顺序阀方式整体流量特性实测曲线和仿真曲线绘制在同一坐标系中,调整仿真参数a、b、c、d使两条曲线拟合,则公式(10)即为顺序阀方式下整体流量特性仿真模型。
6.根据权利要求5所述的优化方法,其特征在于:步骤B中,
对单阀整体流量特性实测曲线计算步骤如下,
进行单阀方式下整体流量特性曲线测试实验,实验测量获得数据向量:调门总流量指令FDEMsig,调门开度指令OPENsig,调节级压力PIsig,主汽压力TPsig,向量的维度为nsig×1,且有FDEMsig(1)>0,FDEMsig(nsig)=100,表示单阀方式下调门全部开启,使用步骤A中多阀方式下整体流量特性实测曲线的计算步骤,获得单阀方式调门整体流量百分比向量FLsig,以向量FDEMsig为横坐标、向量FLsig为纵坐标构成的曲线即为单阀方式整体流量特性实测曲线;
此外,单阀方式下各调门开度指令相同,以单阀方式下调门开度OPENsig为横坐标,整体流量FLsig为纵坐标,得到整体流量关于调门开度的曲线,用线性分段函数h(x)表示;
单阀整体流量特性仿真的步骤如下,
单阀方式下流量总指令FDEM经过背压修正函数f1(x)修正,乘以单阀方式下最佳阀位系数kVLVPT,再经过流量开度修正函数修正,得到阀门开度指令;
设当前流量总指令为FDEM,FDEM为向量,取值范围为0-100,每隔0.1取一个值,FDEM经背压修正函数f1(x)修正得到流量指令FL1,流量指令FL1乘以最佳阀位系数kVLVPT,得到流量指令再经流量开度函数修正得到各调门开度指令
根据单阀方式下调门开度与整体流量的函数关系h(x)计算单阀方式下整体仿真流量
以向量FDEM为横坐标、向量为纵坐标构成的曲线即为单阀方式调门整体流量特性仿真曲线,公式(13)为单阀方式调门整体流量仿真模型。
7.根据权利要求6所述的优化方法,其特征在于:步骤C中,
单阀整体流量特性曲线斜率优化的步骤如下,
在同一坐标轴下绘制单阀和多阀方式仿真流量特性曲线,调整kVLVPT使单阀仿真流量特性曲线向多阀方式仿真流量特性曲线靠近,从而减小单多阀切换时的扰动。
8.根据权利要求7所述的优化方法,其特征在于:步骤C中,
多阀整体流量特性曲线优化的步骤如下,
1)背压修正函数调整
单阀方式下当总流量指令FDEM=100时,经背压函数和最佳阀位系数后的流量指令满足单阀整体流量特性曲线斜率优化的步骤中调整最佳阀位系数kVLVPT后,背压修正函数根据式(5)和式(11)调整;
2)流量比例偏置因子优化
流量比例偏置因子ki、bi(i=1,2,…,6)控制相应调节阀在总流量指令FDEM分别为某两个值时开始开启和全开;在不考虑重叠度函数的情况下,流量指令分别取值为0和100时对应的FDEM值即为调节阀i的开启位置FDEM0和全开位置FDEM100,根据式(5)和式(6)反算得出
f1(FDEM0)*ki+bi=0
f1(FDEM100)*ki+bi=100(14)
保持偏置因子bi的值不变,增大比例偏置因子ki,调节阀的开启位置FDEM0和全开位置FDEM100均减小,且调节阀所占的FDEM区间缩小,整理流量特性曲线的斜率增大;相反减小ki,整理流量特性曲线的斜率减小;
在同一坐标轴下绘制单阀和多阀方式仿真流量特性曲线,调整流量比例因子ki使多阀仿真流量特性曲线向单阀方式仿真流量特性曲线靠近,从而减小单多阀切换时的扰动。
9.根据权利要求8所述的优化方法,其特征在于:步骤C中,
重叠度函数的优化步骤为,
重叠度函数控制相邻两个调门重叠区域的大小;根据顺序阀调门整体流量仿真曲线中重叠区域流量特性的线性度,调整重叠度函数使重叠区呈现线性关系。
10.根据权利要求9所述的优化方法,其特征在于:步骤C中,
调门流量开度修正函数的优化步骤为,
调门流量开度修正函数是将各调门的流量指令折算成调门开度,为保证单阀方式下各调门的开度保持一致,设置各调门的流量开度修正函数一样;根据多阀方式下单个调门流量特性实测曲线,做一条流量开度曲线使之与各调门的实测曲线的距离之和最小,将此曲线作为优化后的调门流量开度函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610073044.6A CN105697074B (zh) | 2016-02-03 | 2016-02-03 | 一种汽轮机组高压调门流量优化装置及其优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610073044.6A CN105697074B (zh) | 2016-02-03 | 2016-02-03 | 一种汽轮机组高压调门流量优化装置及其优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105697074A true CN105697074A (zh) | 2016-06-22 |
CN105697074B CN105697074B (zh) | 2017-05-10 |
Family
ID=56228901
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610073044.6A Active CN105697074B (zh) | 2016-02-03 | 2016-02-03 | 一种汽轮机组高压调门流量优化装置及其优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105697074B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107133433A (zh) * | 2017-07-10 | 2017-09-05 | 中国大唐集团科学技术研究院有限公司华东分公司 | 一种基于模型自适应汽轮机调门流量特性优化方法 |
CN108278134A (zh) * | 2017-12-28 | 2018-07-13 | 清华大学 | 一种母管制直流锅炉发电机组协调控制方法 |
CN110081410A (zh) * | 2018-04-26 | 2019-08-02 | 广东电网有限责任公司 | 一种汽动给水泵小汽轮机的控制方法 |
CN110118107A (zh) * | 2019-06-03 | 2019-08-13 | 中国神华能源股份有限公司 | 汽轮机组的配汽方法和配汽装置 |
CN110761852A (zh) * | 2019-09-27 | 2020-02-07 | 国网江西省电力有限公司电力科学研究院 | 一种汽轮机调节阀行程仿真建模方法 |
CN111505943A (zh) * | 2020-06-03 | 2020-08-07 | 国电科学技术研究院有限公司 | 基于全行程建模的汽轮机流量特性优化方法 |
CN111931298A (zh) * | 2020-10-19 | 2020-11-13 | 国网江西省电力有限公司电力科学研究院 | 一种汽轮机配汽端仿真计算方法 |
CN111946405A (zh) * | 2020-08-20 | 2020-11-17 | 中国大唐集团科学技术研究院有限公司华东电力试验研究院 | 汽轮机阀序切换优化控制方法、装置、终端及存储介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5257405A (en) * | 1975-11-04 | 1977-05-11 | Westinghouse Electric Corp | Optimum sequence valve position indicating system of turbine generating station |
US4418285A (en) * | 1972-11-15 | 1983-11-29 | Westinghouse Electric Corp. | System and method for controlling a turbine power plant in the single and sequential valve modes with valve dynamic function generation |
JPS6296702A (ja) * | 1985-10-24 | 1987-05-06 | Mitsubishi Heavy Ind Ltd | 親子弁群式流量加減装置 |
CN101871840A (zh) * | 2010-07-12 | 2010-10-27 | 山东电力研究院 | 汽轮机高压调门流量特性在线测试方法 |
CN103046972A (zh) * | 2012-12-13 | 2013-04-17 | 哈尔滨工业大学 | 汽轮机单阀或多阀的一种非线性自动无扰切换方法 |
CN103513640A (zh) * | 2013-10-11 | 2014-01-15 | 国家电网公司 | 一种燃煤机组自动发电控制系统整体优化方法及系统 |
CN103758583A (zh) * | 2014-01-03 | 2014-04-30 | 广东电网公司电力科学研究院 | 基于deh的汽轮机的调门配汽曲线转换装置 |
-
2016
- 2016-02-03 CN CN201610073044.6A patent/CN105697074B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4418285A (en) * | 1972-11-15 | 1983-11-29 | Westinghouse Electric Corp. | System and method for controlling a turbine power plant in the single and sequential valve modes with valve dynamic function generation |
JPS5257405A (en) * | 1975-11-04 | 1977-05-11 | Westinghouse Electric Corp | Optimum sequence valve position indicating system of turbine generating station |
JPS6296702A (ja) * | 1985-10-24 | 1987-05-06 | Mitsubishi Heavy Ind Ltd | 親子弁群式流量加減装置 |
CN101871840A (zh) * | 2010-07-12 | 2010-10-27 | 山东电力研究院 | 汽轮机高压调门流量特性在线测试方法 |
CN103046972A (zh) * | 2012-12-13 | 2013-04-17 | 哈尔滨工业大学 | 汽轮机单阀或多阀的一种非线性自动无扰切换方法 |
CN103513640A (zh) * | 2013-10-11 | 2014-01-15 | 国家电网公司 | 一种燃煤机组自动发电控制系统整体优化方法及系统 |
CN103758583A (zh) * | 2014-01-03 | 2014-04-30 | 广东电网公司电力科学研究院 | 基于deh的汽轮机的调门配汽曲线转换装置 |
Non-Patent Citations (2)
Title |
---|
武海澄: "DEH单顺阀切换控制参数的试验研究", 《吉林电力》 * |
武海澄: "汽轮机DEH阀门流量特性仿真方法研究", 《安徽电力》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107133433A (zh) * | 2017-07-10 | 2017-09-05 | 中国大唐集团科学技术研究院有限公司华东分公司 | 一种基于模型自适应汽轮机调门流量特性优化方法 |
CN107133433B (zh) * | 2017-07-10 | 2020-09-01 | 中国大唐集团科学技术研究院有限公司华东分公司 | 一种基于模型自适应汽轮机调门流量特性优化方法 |
CN108278134B (zh) * | 2017-12-28 | 2020-07-14 | 清华大学 | 一种母管制直流锅炉发电机组协调控制方法 |
CN108278134A (zh) * | 2017-12-28 | 2018-07-13 | 清华大学 | 一种母管制直流锅炉发电机组协调控制方法 |
CN110081410A (zh) * | 2018-04-26 | 2019-08-02 | 广东电网有限责任公司 | 一种汽动给水泵小汽轮机的控制方法 |
CN110081410B (zh) * | 2018-04-26 | 2020-06-09 | 广东电网有限责任公司 | 一种汽动给水泵小汽轮机的控制方法 |
CN110118107A (zh) * | 2019-06-03 | 2019-08-13 | 中国神华能源股份有限公司 | 汽轮机组的配汽方法和配汽装置 |
CN110761852A (zh) * | 2019-09-27 | 2020-02-07 | 国网江西省电力有限公司电力科学研究院 | 一种汽轮机调节阀行程仿真建模方法 |
CN111505943A (zh) * | 2020-06-03 | 2020-08-07 | 国电科学技术研究院有限公司 | 基于全行程建模的汽轮机流量特性优化方法 |
CN111505943B (zh) * | 2020-06-03 | 2022-08-16 | 国电科学技术研究院有限公司 | 基于全行程建模的汽轮机流量特性优化方法 |
CN111946405A (zh) * | 2020-08-20 | 2020-11-17 | 中国大唐集团科学技术研究院有限公司华东电力试验研究院 | 汽轮机阀序切换优化控制方法、装置、终端及存储介质 |
CN111931298A (zh) * | 2020-10-19 | 2020-11-13 | 国网江西省电力有限公司电力科学研究院 | 一种汽轮机配汽端仿真计算方法 |
CN111931298B (zh) * | 2020-10-19 | 2021-03-02 | 国网江西省电力有限公司电力科学研究院 | 一种汽轮机配汽端仿真计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105697074B (zh) | 2017-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105697074B (zh) | 一种汽轮机组高压调门流量优化装置及其优化方法 | |
CN107133433B (zh) | 一种基于模型自适应汽轮机调门流量特性优化方法 | |
CN110161841A (zh) | 一种适用于暂冲式跨声速风洞的前馈-模糊pid控制方法 | |
CN102841539A (zh) | 基于多模型预测控制的亚临界协调控制方法 | |
CN104849052A (zh) | 一种喷嘴配汽汽轮机流量特性试验方法 | |
DE102010037368A1 (de) | Verfahren zur Steuerung einer Turboladeranordnung eines Verbrennungsmotors sowie Turboladeranordnung | |
CN104734588B (zh) | 一种生物质气内燃发电机组转速控制方法 | |
CN103322647B (zh) | 一种中央空调的冷却水供水温度预测控制方法 | |
CN110879620A (zh) | 一种核电站立式蒸汽发生器液位控制方法以及系统 | |
CN105425612A (zh) | 一种水轮机调节系统控制参数的优选方法 | |
CN106682735A (zh) | 基于pid调节的bp神经网络算法 | |
CN103324828A (zh) | 基于运行方式灵敏度的电力系统低频振荡辅助决策方法 | |
CN103516274A (zh) | 一种包含有多类型限制单元的发电机励磁系统 | |
CN104456513B (zh) | 一种火力发电厂主汽温预估优化控制方法 | |
CN111322885B (zh) | 一种间冷系统百叶窗控制装置及其控制方法 | |
CN109541935A (zh) | 一种参数自适应分数阶自抗扰自动发电控制方法 | |
CN108131238B (zh) | 一种抑制水锤压力波动的pid控制方法 | |
CN105449698A (zh) | 一种新型水电机组负荷与频率控制器 | |
CN115236977A (zh) | 一种风洞连续变速压流场控制方法 | |
CN107906499B (zh) | 火电机组智能控制系统 | |
CN111946405B (zh) | 汽轮机阀序切换优化控制方法、装置、终端及存储介质 | |
Burlibaşa et al. | Frequency domain design of gain scheduling control for large wind systems in full-load region | |
CN105156159B (zh) | 一种基于调门最佳重叠度的汽轮机流量特性优化方法 | |
CN107687668A (zh) | 一种供热系统自适应动态热负荷调控方法 | |
CN108631368B (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 |