CN107103167B - 针对甩负荷工况的deh调速功能诊断方法及系统 - Google Patents
针对甩负荷工况的deh调速功能诊断方法及系统 Download PDFInfo
- Publication number
- CN107103167B CN107103167B CN201710479894.0A CN201710479894A CN107103167B CN 107103167 B CN107103167 B CN 107103167B CN 201710479894 A CN201710479894 A CN 201710479894A CN 107103167 B CN107103167 B CN 107103167B
- Authority
- CN
- China
- Prior art keywords
- deh
- diagnosis
- speed regulation
- regulation function
- response time
- 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.)
- Expired - Fee Related
Links
- 238000003745 diagnosis Methods 0.000 title claims abstract description 91
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000004088 simulation Methods 0.000 claims abstract description 65
- 230000004044 response Effects 0.000 claims abstract description 41
- 238000012360 testing method Methods 0.000 claims abstract description 29
- 230000001105 regulatory effect Effects 0.000 claims abstract description 19
- 230000002159 abnormal effect Effects 0.000 claims abstract description 15
- 238000004092 self-diagnosis Methods 0.000 claims abstract description 11
- 230000010354 integration Effects 0.000 claims description 6
- 238000012888 cubic function Methods 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 230000005856 abnormality Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 description 23
- 238000010586 diagram Methods 0.000 description 13
- 230000001133 acceleration Effects 0.000 description 12
- 238000005070 sampling Methods 0.000 description 11
- 230000008569 process Effects 0.000 description 8
- 230000008859 change Effects 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 3
- 238000013178 mathematical model Methods 0.000 description 3
- 238000007664 blowing Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005422 blasting Methods 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000002405 diagnostic procedure Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000007670 refining Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Economics (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Geometry (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Control Of Turbines (AREA)
Abstract
本发明公开了一种针对甩负荷工况的DEH调速功能诊断方法及系统，其中该方法包括DEH混合仿真至满负荷后，断开并网开关，开始甩负荷仿真试验，DEH转速仿真逻辑模块对输入的高压调门反馈和中压调门反馈进行运算，并输出汽轮机的仿真转速；飞升曲线特征数据自动获取逻辑模块自动采集汽轮机的仿真转速，自动捕获飞升曲线的特征数据；同时SOE系统计算出DEH电气响应时间和液压响应时间；DEH调速功能自诊断逻辑模块接收飞升曲线的特征数据以及DEH电气响应时间和液压响应时间，并按照预设的优先级从高到低分级逐级进行DEH调速功能诊断，每一级诊断都具有独立表决权，表决之后其后的诊断级失效，表决结果作为诊断结论，并输出调速功能异常的关键故障点。
Description
技术领域
本发明属于发电领域，尤其涉及一种针对甩负荷工况的DEH调速功能诊断方法及系统。
背景技术
当今电网规模日趋庞大和电力消费复杂多变的情形，导致电力系统存在着各种扰动，而且恶劣的气候变化和自然灾害的频发，使得汽轮发电机组运行过程中经常面临发生甩负荷情况。快速甩负荷是最恶劣的一种工况，它涉及机、炉、电所有系统，对机组的冲击很大并直接危害到机组的安全稳定运行。
对汽轮发电机组运行来说，最关键的是确保甩负荷后机组转速快速稳定到额定转速，避免转速飞升导致汽轮机跳闸。
DEH(即汽轮机数字电液控制系统(Digital Electric Hydraulic ControlSystem)，简称数字电调，是DCS的重要组成部分)是管理和控制大型汽轮机运行的核心。其负荷-转速控制部分的工作状态和动态性能直接影响着整个机组的正常运行，尤其是DEH在机组甩负荷状态下的工作过程和控制品质更为人们所关注。在实际运行的机组上不可能进行大规模、大扰动、多工况、多种参数匹配的试验。而且DEH逻辑回路和硬回路非常复杂繁琐，往往通过常规仿真故障点不易被发现。近年来甩负荷试验出现的问题很多，如OPC回路二极管击穿导致电磁阀不动作、OPC时转速PID不跟踪导致OPC复位后调阀开度过大超速、OPC发生后调阀指令切零存在速率、DEH控制模式不自动切换，导致中调门指令始终置全开等种种隐蔽性强的问题，很有必要对这些经验教训总结提炼形成一套较完备并持续升级的诊断系统，通过对机组DEH调速功能的预诊断，可更有效全面反映出机组甩负荷工况下汽轮机DEH转速控制情况。
目前国内外研究主要是对汽轮机的油动机建立数学模型，蒸汽容积效应建立数学模型，再根据机组并网运行时的控制回路推算出机组甩负荷的数学模型，在MATLAB等工具下模拟汽轮机发生甩负荷事故后，汽轮机阀门的关闭时间，汽轮机转速飞升速率及预估值。
国内电厂DEH转速仿真信号普遍是通过对汽机调门开度信号进行一阶惯性、速率、线性化环节简单运算得出。这种仿真模型在模拟甩负荷工况方面精确度较差，和实际工况差距很大，DEH调速性能无法得到最大考验。
如图1传统的针对甩负荷工况的DEH调速功能诊断方法的流程为：DEH混合仿真至满负荷，甩负荷仿真试验开始，DEH仿真逻辑执行，当汽轮机仿真转速为N时，手动调取历史趋势，观察仿真转速最大飞升值，最后得到诊断结论。其中，DEH仿真逻辑如图2所示，其表达式为：
其中，N表示汽轮机的转速；FH表示高调门反馈量；FI表示中调门反馈量；T1和T2分别表示一阶惯性采样时间系数；K1、K2和K3表示分别比例常系数。
但是，传统控制方案存在以下不足：
(1)传统汽轮机转速仿真逻辑中，流量的计算直接用调门开度代替，没有考虑调门流量特性的因素，导致进汽流量环节的计算线性很差。
(2)传统转速仿真逻辑为开环回路，没有考虑鼓风摩擦和轴承机械摩擦损失负反馈作用。
(3)转速通过对进汽流量一阶惯性、比例环节简单运算得出，计算出的转速与实际甩负荷转速飞升曲线差距较大，不能满足自动诊断系统对转速仿真的要求。
(4)传统DEH甩负荷仿真试验不能实现DEH甩负荷仿真试验后的特征数据自动采集和自动处理。不能实现DEH甩负荷仿真试验后的全方位(包括DEH电气软件、电气硬件和液压部分)自动诊断，这样可能使某些隐患继续隐藏，时刻威胁机组安全运行。
发明内容
为了解决现有技术的不足，本发明提供了一种针对甩负荷工况的DEH调速功能诊断系统，其实现了DEH甩负荷仿真试验后的全方位(包括DEH电气软件、电气硬件和液压部分)自动诊断，并自动给出诊断结果和处理建议，实现对DEH深层次、无死角隐患查找，保障机组安全运行。
本发明的一种针对甩负荷工况的DEH调速功能诊断方法，包括：
步骤1：DEH混合仿真至满负荷后，断开并网开关，开始甩负荷仿真试验，DEH转速仿真逻辑模块对输入的高压调门反馈和中压调门反馈进行运算，并输出汽轮机的仿真转速；
步骤2：飞升曲线特征数据自动获取逻辑模块自动采集汽轮机的仿真转速，自动捕获飞升曲线的特征数据；同时SOE系统自动录取发变组脱网信号、调门快关指令和调门关到位信号的变位时间，读取相应时间标签后，计算出DEH电气响应时间和液压响应时间；
步骤3：DEH调速功能自诊断逻辑模块接收飞升曲线的特征数据以及DEH电气响应时间和液压响应时间，并按照预设的优先级从高到低分级逐级进行DEH调速功能诊断，每一级诊断都具有独立表决权，表决之后其后的诊断级失效，表决结果作为诊断结论，并输出调速功能异常的关键故障点。
进一步的，所述步骤1中的DEH转速仿真逻辑模块的控制逻辑为：
高调门反馈FBH通过高调门流量特性反函数F1(X)和做功比例K1，计算高压缸进汽流量FH；中调门反馈FBI通过中调门流量特性反函数F2(X)和做功比例K2，计算中低压缸进汽流量FI；
高压缸进汽流量FH和中低压缸进汽流量FI经加法模块，得出汽轮机进汽流量F；
汽轮机进汽流量F经一阶惯性模块，计算得出汽轮机功率E；
汽轮机的仿真转速N经三次方函数和比例K3运算，得出鼓风摩擦损失LD1，仿真转速经比例K4运算得出轴承机械摩擦损失LD2，LD1和LD2叠加构成汽轮机负荷LD；
汽轮机功率E与汽轮机负荷LD经减法模块，得出功率-负荷不平衡量；
再经切换模块判断，如果判断条件TR为1即脱网信号且试验投入，选择功率-负荷不平衡量输入积分模块，否则，选择0输入；功率-负荷不平衡量经积分模块运算，最终得出仿真转速N；其中，FBH、FBI、K1、K2、K3和K4均为常数；F1(X)和F2(X)均为预设的函数。
在DEH转速仿真逻辑模块中，对汽轮机调门反馈信号进行流量特性反函数修正，将非线性的反馈信号变为线性的进汽流量信号，提高了汽轮机功率计算的线性度。
进一步的，所述步骤2中飞升曲线特征数据自动获取逻辑模块的控制逻辑为：
试验投入且发电机脱网后开始执行飞升曲线特征数据自动获取逻辑模块中的控制逻辑；
选取汽轮机转速的峰值P，然后把峰值输入第一微分模块，判断第一微分运算值是否等于0且保持预设时间，若是，则取汽轮机转速的谷值B，然后把谷值输入第二微分模块，判断第二运算值是否等于0且保持预设时间，若是，置周期计数器NP为1，并将波峰波谷值分别寄存到P1、B1中，直到周期计数器NP达到预设的最大周期，则运算结束；
根据获取的甩负荷飞升曲线的波峰和波谷值，计算调速衰减率和稳态最大波动值；
将发变组脱网信号、调门快关指令和调门关到位信号接入SOE系统，甩负荷仿真试验后调取SOE数据，发变组脱网信号、调门快关指令、调门关到位变位时间分别为：TO、Tf、Tc，计算DEH电气响应时间Te＝Tf-TO，液压响应时间Th＝Tc–Tf。
进一步的，在所述步骤3中的DEH调速功能自诊断逻辑模块的控制逻辑为：
按照优先级从高到低共设置五个诊断级，分别为第一诊断级～第五诊断级；
第一诊断级为判断DEH电气响应时间是否大于预设DEH电气响应时间阈值；
第二诊断级为判断液压响应时间是否大于预设液压响应时间阈值；
第三诊断级为判断汽轮机转速是否大于预设转速阈值；
第四诊断级为判断调速衰减是否大于预设调速衰减阈值；
第五诊断级为判断稳态最大波动值是否大于预设稳态最大波动值阈值。
进一步的，在第一诊断级过程中，若DEH电气响应时间大于预设DEH电气响应时间阈值，则诊断为DEH超速控制硬回路异常；否则，进入第二诊断级。
进一步的，在第二诊断级过程中，若液压响应时间大于预设液压响应时间阈值，则诊断为DEH超速控制液压回路异常；否则，进入第三诊断级。
进一步的，在第三诊断级过程中，若汽轮机转速大于预设转速阈值，则诊断为DEH超速控制软回路异常；否则，进入第四诊断级。
进一步的，在第四诊断级过程中，若调速衰减大于预设调速衰减阈值，则诊断为调速PID控制参数异常；否则，进入第五诊断级。
进一步的，在第五诊断级过程中，若稳态最大波动值大于预设稳态最大波动值阈值，则诊断为DEH转速控制稳态偏差大，PID参数需要优化；否则，给出调速功能正常的结论，表示DEH系统具备甩负荷热态试验的条件。
本发明还提供了一种针对甩负荷工况的DEH调速功能诊断系统。
本发明的针对甩负荷工况的DEH调速功能诊断系统，采用上述所述的针对甩负荷工况的DEH调速功能诊断方法来实现。
与现有技术相比，本发明的有益效果是：
(1)汽轮机转速仿真逻辑中，对汽轮机调门反馈信号进行流量特性反函数修正，将不线性的反馈信号变为线性的进汽流量信号，提高了汽轮机功率计算的线性度。
(2)本发明将各部分余汽容积综合计算，等效为一个一阶惯性环节，简化了仿真模型，这样就可根据以往同等级容量机组实测甩负荷曲线计算余汽容积时间常数T1，减小了计算误差，可仿真出更接近实际的标准甩负荷转速飞升曲线。负荷计算中量化配比鼓风摩擦和轴承机械摩擦损失。
(3)本发明采用基于汽轮机功率-负荷不平衡量计算的汽轮机转速仿真逻辑，用一个积分环节模拟汽轮机对功率和负荷偏差的响应，极大提高了转速仿真的精确度。
(4)本发明实现DEH甩负荷仿真试验后的特征数据自动采集和自动处理。实现DEH甩负荷仿真试验后的全方位(包括DEH电气软件、电气硬件和液压部分)自动诊断，并自动给出诊断结果和处理建议，实现对DEH深层次、无死角隐患查找，保障机组安全运行。
附图说明
构成本申请的一部分的说明书附图用来提供对本申请的进一步理解，本申请的示意性实施例及其说明用于解释本申请，并不构成对本申请的不当限定。
图1是传统的针对甩负荷工况的DEH调速功能诊断方法流程图。
图2是传统的针对甩负荷工况的DEH调速功能诊断系统的DEH仿真逻辑图。
图3是本发明的针对甩负荷工况的DEH调速功能诊断方法流程图。
图4是本发明的脱网状态下汽轮机转速仿真逻辑图。
图5是本发明的飞升曲线特征数据自动获取逻辑框图。
图6(a)是本发明的一个周期的汽轮机转速峰值和谷值获取逻辑图。
图6(b)是本发明的周期计数逻辑图。
图6(c)是本发明的第一个周期的汽轮机转速峰值和谷值存储逻辑图。
图6(d)是本发明的第二个周期的汽轮机转速峰值和谷值存储逻辑图。
图6(e)是本发明的第三个周期的汽轮机转速峰值和谷值存储逻辑图。
图6(f)是本发明的第四个周期的汽轮机转速峰值和谷值存储逻辑图。
图7是本发明的调速功能自诊断逻辑图。
具体实施方式
应该指出，以下详细说明都是例示性的，旨在对本申请提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本申请的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
本发明涉及一种针对甩负荷工况的DEH调速功能诊断方法，在原诊断方法的基础上进行改进，主要包含三个部分：转速仿真逻辑、转速飞升曲线特征数据自动获取逻辑、DEH调速功能自诊断逻辑。
如图3所示，本发明的一种针对甩负荷工况的DEH调速功能诊断方法，包括：
步骤1：DEH混合仿真至满负荷后，断开并网开关，开始甩负荷仿真试验，DEH转速仿真逻辑模块对输入的高压调门反馈和中压调门反馈进行运算，并输出汽轮机的仿真转速；
步骤2：飞升曲线特征数据自动获取逻辑模块自动采集汽轮机的仿真转速，自动捕获飞升曲线的特征数据；同时SOE系统自动录取发变组脱网信号、调门快关指令和调门关到位信号的变位时间，读取相应时间标签后，计算出DEH电气响应时间和液压响应时间；
步骤3：DEH调速功能自诊断逻辑模块接收飞升曲线的特征数据以及DEH电气响应时间和液压响应时间，并按照预设的优先级从高到低分级逐级进行DEH调速功能诊断，每一级诊断都具有独立表决权，表决之后其后的诊断级失效，表决结果作为诊断结论，并输出调速功能异常的关键故障点。
下面结合附图与实例进一步说明如下。
1、转速仿真逻辑
如图4所示，高调门反馈FBH通过高调门流量特性反函数F1(X)和做功比例K1计算高压缸进汽流量FH，中调门反馈FBI通过高调门流量特性反函数F2(X)和做功比例K2计算中低压缸进汽流量FI。高压缸进汽流量FH和中低压缸进汽流量FI经加法模块得出汽轮机进汽流量F。汽轮机进汽流量F经一阶惯性模块(余汽容积环节)计算得出汽轮机功率E。
仿真转速N经三次方函数和比例K3运算得出鼓风摩擦损失LD1，仿真转速经比例K4运算得出轴承机械摩擦损失LD2。LD1和LD2叠加后生成汽轮机负荷LD。
汽轮机功率E、汽轮机负荷LD经减法模块得出功率-负荷不平衡量。经切换模块判断，如果判断条件TR(脱网信号且试验投入)为1，选择功率-负荷不平衡量输入积分块；TR为0，选择0输入。功率-负荷不平衡量经积分块运算最终得出仿真转速N。
仿真逻辑相关参数计算方法：
仿真逻辑相关参数计算方法以某350MW机组甩负荷实测曲线为实例说明，下面是计算过程，进气流量F、功率E和负荷LD的量纲统一为0-100标准量，转换公式为：工程量＝(标准量/100)*额定值。
(a)汽轮机进汽流量的计算
进汽流量在用调门开度计算的基础上，采用调门流量特性反函数修正，具体公式如下：
总流量F＝高压缸功率FH+中低压缸功率FI；
高压缸和中低压缸功率比例按1:3计算。
高调门反馈FBH通过高调门流量特性反函数计算流量指令：FDEM＝F1(FBH)。
高压缸功率FH＝FDEM*0.75。
中调门反馈FBI通过中调门流量特性反函数计算流量指令：FDEM＝F2(FBI)。
中压缸功率FI＝FDEM*0.25。
(b)汽轮机负荷的计算
汽轮机转子旋转阻力包括转子叶片鼓风摩擦作用和轴承机械摩擦阻力，产生功率损失成为汽轮机负荷的一部分。鼓风摩擦：LD1＝K3*N^3，轴承机械摩擦：LD2＝K4*N。
K3、K4计算方法：转速在3000±300变化时，K3、K4变化较小。为简化计算，取3100转为计算参考点，计算出的数值适用3000±300转速变化范围。
3100转左右的功率损失可根据实测转速飞升曲线算出，转速由最大飞升转速开始降速，降速曲线接近线性。此时的加速度基本为恒定值a2，根据降速曲线算出3100左右的汽轮机负荷标准量为LD＝a2/a0。由于高转速时，叶片鼓风摩擦作用远大于轴承机械摩擦阻力，LD1≈5*LD2。这样可以计算出K3、k4。
a0和a2的算法，由飞升曲线录波数据得出：发电机脱网后第一个采样周期转速为N0，N1，采样时间T，转子最大飞升加速度a0＝(N1-N0)/T。
在转速飞升后的降速曲线录波数据中，在3100转左右选取两个取样点，N1，N2，时间间隔T，计算降速率a2＝(N1-N2)/T
算例：在转速飞升后的降速曲线录波数据中，在3100转左右选取两个取样点，N1＝3106.6，N2＝3100.9，时间间隔T＝2S，计算降速率a2＝(N1-N2)/T＝2.85，得出3100转时的汽轮机负荷为LD3100＝(a2/a0)*100＝(2.85/285)X100＝1。LD＝6*LD2＝6*K5*3100＝1，K5＝5.4E-5。
LD＝1.2LD1＝1.2*K4*31003＝1，K4＝2.8E-11。
(c)汽轮机功率计算
调门关闭后汽轮机内的余汽容积导致转速仍有一段飞升，将各部分余汽容积综合计算，等效为一个一阶惯性环节。汽轮机进汽流量经余汽容积环节计算得出汽轮机功率。对调门关闭后纯容积环节分析，调门全关时的功率变化率VE，转速加速度的变化率Va，转速加速度为a1，由实际飞升曲线可以算出转速加速度a1和转速加速度变化率Va。因为VE＝K*Va，由此可算出余汽容积时间常数T1＝a1/Va。
a1和Va的算法，由飞升曲线录波数据得出:调门全关时第一个采样周期转速为N0，依次为N1、N2、N3，采样周期为T
计算加速度a1＝(N1-N0)/T，a11＝(N2-N1)/T
计算加速度变化率Va＝(a1-a11)/T
算例：由飞升曲线录波数据得出:调门全关时第一个采样周期转速为N0＝3083.99，依次为N1＝3097.18、N2＝3107.8、N3＝3115.46，采样周期为T＝100ms
计算加速度a1＝(N1-N0)/T＝131.86，a2＝(N2-N1)/T＝106.29
计算加速度变化率Va＝(a2-a1)/T＝255.8
余汽容积时间常数T1＝a1/Va＝0.51s
(d)转速计算
功率-负荷不平衡量与转速加速度成正比，将汽轮机功率与负荷的差值经积分环节计算得出汽轮机转速。在功率和负荷差值达100％额定负荷时，转子最大飞升加速度a0，算出积分时间常数T2＝100/a0。
算例：由飞升曲线录波数据得出：发电机脱网后第一个采样周期转速为N0＝2999.8，N1＝3028.3，采样时间T＝100ms，转子最大飞升加速度a0＝(N1-N0)/T＝28.5/0.1＝285
算出积分时间T2＝100/a0＝0.35s。
2、飞升曲线特征数据自动获取逻辑
图5为飞升曲线特征数据自动获取逻辑框图。图6(a)-图6(f)为具体的逻辑图，说明如下：
如图6(a)所示，试验投入并且发电机脱网后逻辑开始执行，汽轮机转速输入取峰值P运算程序，然后把峰值输入微分运算程序，微分运算值经是否等于0判断，输入接通延时模块，接通延时模块如果为1，将转速输入取谷值B运算程序，然后把谷值输入微分运算程序，微分运算值经是否等于0判断，输入接通延时模块，接通延时模块如果为1，置周期计数器NP为1，并将波峰波谷值分别寄存到P1、B1中。周期计数器NP经是否等于4判断，之后逻辑重复执行，开始第二周期运算，直到第四周期运算结束，如图6(b)。
通过图6(c)-图6(f)分别为第一个周期的汽轮机转速峰值和谷值存储逻辑图～第四个周期的汽轮机转速峰值和谷值存储逻辑图。本发明获取甩负荷飞升曲线四个周期的波峰(P1-P4)、波谷值(B1-B4)。计算调速衰减率AR＝(P2-P3)/(P2-3000)，稳态最大波动值SW＝P4-B4。
将发变组脱网信号、调门快关指令、调门关到位信号接入SOE系统，实现方便、准确测取电气、液压延迟时间。甩负荷仿真试验后调取SOE数据，发变组脱网信号、调门快关指令、调门关到位变位时间分别为：TO、Tf、Tc。计算DEH电气响应时间Te＝Tf-TO，液压响应时间Th＝Tc–Tf。
3、调速功能自诊断方案
甩负荷转速飞升曲线特征数据获取后，用图7的调速功能自诊断逻辑对甩负荷后的DEH调速功能全方位自动诊断。自诊断逻辑采用按优先级分级诊断，共有五个诊断级，从上到下优先级逐级降低。高优先级诊断是低优先级诊断的基础，每一级诊断都具有独立表决权，表决之后其后的诊断级失效。表决结果作为诊断结论，指出调速功能异常的关键故障点。如果程序通过五级诊断，给出调速功能正常的结论，表示DEH系统具备甩负荷热态试验的条件。
例如：
如图7所示，飞升曲线特征数据自动获取逻辑执行结束，DEH调速功能自诊断逻辑开始执行第一级诊断块。通过判断电气响应时间Te是否大于100ms，如果大于，诊断为DEH超速控制硬回路异常，重点检查DEH控制机柜内硬回路是否存在故障；如果小于，进入第二级诊断块。
通过判断液压响应时间Th是否大于500ms，如果大于，诊断为DEH超速控制液压回路异常，重点检查就地液压部件是否存在故障；如果小于，进入第三级诊断块。
通过判断转速是否大于或等于3090，如果大于，诊断为DEH超速控制软回路异常，重点检查超速控制逻辑是否存在错误，如调门快关指令发出后PID指令是否快速跟踪到零；如果小于，进入第四级诊断块。
通过判断衰减率AR是否小于0.25，如果小于，诊断为调速PID控制参数异常，重点检查PID的比例是否合适；如果大于0.25(说明控制系统稳定，参数基本合适，需下一级诊断是否需要细调)，进入第五级诊断块。
通过判断稳态最大波动值SW是否大于5转，如果大于，诊断为DEH转速控制稳态偏差大，PID参数需要优化，重点检查调速PID的积分是否设置合适；如果小于，给出调速功能正常的结论，表示DEH系统具备甩负荷热态试验的条件。
本发明实现DEH甩负荷仿真试验后的特征数据自动采集和自动处理。实现DEH甩负荷仿真试验后的全方位(包括DEH电气软件、电气硬件和液压部分)自动诊断，并自动给出诊断结果和处理建议，实现对DEH深层次、无死角隐患查找，保障机组安全运行。
上述虽然结合附图对本发明的具体实施方式进行了描述，但并非对本发明保护范围的限制，所属领域技术人员应该明白，在本发明的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
Claims (9)
1.一种针对甩负荷工况的DEH调速功能诊断方法，其特征在于，包括：
步骤1：DEH混合仿真至满负荷后，断开并网开关，开始甩负荷仿真试验，DEH转速仿真逻辑模块对输入的高压调门反馈和中压调门反馈进行运算，并输出汽轮机的仿真转速；
步骤2：飞升曲线特征数据自动获取逻辑模块自动采集汽轮机的仿真转速，自动捕获飞升曲线的特征数据；同时SOE系统自动录取发变组脱网信号、调门快关指令和调门关到位信号的变位时间，读取相应时间标签后，计算出DEH电气响应时间和液压响应时间；
步骤3：DEH调速功能自诊断逻辑模块接收飞升曲线的特征数据以及DEH电气响应时间和液压响应时间，并按照预设的优先级从高到低分级逐级进行DEH调速功能诊断，每一级诊断都具有独立表决权，表决之后其后的诊断级失效，表决结果作为诊断结论，并输出调速功能异常的关键故障点；
所述步骤1中的DEH转速仿真逻辑模块的控制逻辑为：
高调门反馈FBH通过高调门流量特性反函数F1(X)和做功比例K1，计算高压缸进汽流量FH；中调门反馈FBI通过中调门流量特性反函数F2(X)和做功比例K2，计算中低压缸进汽流量FI；
高压缸进汽流量FH和中低压缸进汽流量FI经加法模块，得出汽轮机进汽流量F；
汽轮机进汽流量F经一阶惯性模块，计算得出汽轮机功率E；
汽轮机的仿真转速N经三次方函数和比例K3运算，得出鼓风摩擦损失LD1，仿真转速经比例K4运算得出轴承机械摩擦损失LD2，LD1和LD2叠加构成汽轮机负荷LD；
汽轮机功率E与汽轮机负荷LD经减法模块，得出功率-负荷不平衡量；
再经切换模块判断，如果判断条件TR为1即脱网信号且试验投入，选择功率-负荷不平衡量输入积分模块，否则，选择0输入；功率-负荷不平衡量经积分模块运算，最终得出仿真转速N；其中，FBH、FBI、K1、K2、K3和K4均为常数；F1(X)和F2(X)均为预设的函数。
2.如权利要求1所述的针对甩负荷工况的DEH调速功能诊断方法，其特征在于，所述步骤2中飞升曲线特征数据自动获取逻辑模块的控制逻辑为：
试验投入且发电机脱网后开始执行飞升曲线特征数据自动获取逻辑模块中的控制逻辑；
选取汽轮机转速的峰值P，然后把峰值输入第一微分模块，判断第一微分运算值是否等于0且保持预设时间，若是，则取汽轮机转速的谷值B，然后把谷值输入第二微分模块，判断第二运算值是否等于0且保持预设时间，若是，置周期计数器NP为1，并将波峰波谷值分别寄存到P1、B1中，直到周期计数器NP达到预设的最大周期，则运算结束；
根据获取的甩负荷飞升曲线的波峰和波谷值，计算调速衰减率和稳态最大波动值；
将发变组脱网信号、调门快关指令和调门关到位信号接入SOE系统，甩负荷仿真试验后调取SOE数据，发变组脱网信号、调门快关指令、调门关到位变位时间分别为：TO、Tf、Tc，计算DEH电气响应时间Te＝Tf-TO，液压响应时间Th＝Tc–Tf。
3.如权利要求1所述的针对甩负荷工况的DEH调速功能诊断方法，其特征在于，在所述步骤3中的DEH调速功能自诊断逻辑模块的控制逻辑为：
按照优先级从高到低共设置五个诊断级，分别为第一诊断级～第五诊断级；
第一诊断级为判断DEH电气响应时间是否大于预设DEH电气响应时间阈值；
第二诊断级为判断液压响应时间是否大于预设液压响应时间阈值；
第三诊断级为判断汽轮机转速是否大于预设转速阈值；
第四诊断级为判断调速衰减是否大于预设调速衰减阈值；
第五诊断级为判断稳态最大波动值是否大于预设稳态最大波动值阈值。
4.如权利要求3所述的针对甩负荷工况的DEH调速功能诊断方法，其特征在于，在第一诊断级过程中，若DEH电气响应时间大于预设DEH电气响应时间阈值，则诊断为DEH超速控制硬回路异常；否则，进入第二诊断级。
5.如权利要求3所述的针对甩负荷工况的DEH调速功能诊断方法，其特征在于，在第二诊断级过程中，若液压响应时间大于预设液压响应时间阈值，则诊断为DEH超速控制液压回路异常；否则，进入第三诊断级。
6.如权利要求3所述的针对甩负荷工况的DEH调速功能诊断方法，其特征在于，在第三诊断级过程中，若汽轮机转速大于预设转速阈值，则诊断为DEH超速控制软回路异常；否则，进入第四诊断级。
7.如权利要求3所述的针对甩负荷工况的DEH调速功能诊断方法，其特征在于，在第四诊断级过程中，若调速衰减大于预设调速衰减阈值，则诊断为调速PID控制参数异常；否则，进入第五诊断级。
8.如权利要求3所述的针对甩负荷工况的DEH调速功能诊断方法，其特征在于，在第五诊断级过程中，若稳态最大波动值大于预设稳态最大波动值阈值，则诊断为DEH转速控制稳态偏差大，PID参数需要优化；否则，给出调速功能正常的结论，表示DEH系统具备甩负荷热态试验的条件。
9.一种针对甩负荷工况的DEH调速功能诊断系统，其特征在于，该系统采用如权利要求1-8中任一所述的针对甩负荷工况的DEH调速功能诊断方法来实现。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710479894.0A CN107103167B (zh) | 2017-06-22 | 2017-06-22 | 针对甩负荷工况的deh调速功能诊断方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710479894.0A CN107103167B (zh) | 2017-06-22 | 2017-06-22 | 针对甩负荷工况的deh调速功能诊断方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107103167A CN107103167A (zh) | 2017-08-29 |
| CN107103167B true CN107103167B (zh) | 2020-04-03 |
Family
ID=59664064
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710479894.0A Expired - Fee Related CN107103167B (zh) | 2017-06-22 | 2017-06-22 | 针对甩负荷工况的deh调速功能诊断方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107103167B (zh) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109543222B (zh) * | 2018-10-19 | 2023-04-18 | 广西电网有限责任公司电力科学研究院 | 一种不同功率水平下压水堆核电机组甩负荷保护仿真方法 |
| CN111123770B (zh) * | 2019-12-13 | 2021-09-24 | 山东中实易通集团有限公司 | 一种fcb工况下旁路模型开度确定方法及装置 |
| CN111175536A (zh) * | 2020-03-03 | 2020-05-19 | 贵州电网有限责任公司 | 一种发电机组甩负荷测试系统及测试方法 |
| CN111220875B (zh) * | 2020-03-09 | 2025-03-25 | 中国大唐集团科学技术研究院有限公司华东电力试验研究院 | 一种静态模拟甩负荷试验集成测试系统及测试方法 |
| CN111623988B (zh) * | 2020-05-28 | 2022-11-18 | 国网河北省电力有限公司电力科学研究院 | 一种高位布置汽轮机组甩负荷试验安全性评估方法 |
| CN113361107B (zh) * | 2021-06-08 | 2023-04-07 | 中国大唐集团科学技术研究院有限公司中南电力试验研究院 | 一种汽轮机甩负荷后转速飞升的静态仿真预测方法 |
| CN113568397A (zh) * | 2021-08-09 | 2021-10-29 | 陕西能源麟北发电有限公司 | 一种汽轮机监测仪表自检测系统及自检方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH08218897A (ja) * | 1995-02-17 | 1996-08-27 | Toshiba Corp | ガスタービン制御装置 |
| CN1959067A (zh) * | 2005-11-03 | 2007-05-09 | 沈阳创思达自动化系统有限公司 | 汽轮机的调速控制系统 |
| CN103558031A (zh) * | 2013-11-11 | 2014-02-05 | 沈阳黎明航空发动机(集团)有限责任公司 | 一种燃气轮机甩负荷试验方法 |
| CN105700380A (zh) * | 2016-01-24 | 2016-06-22 | 浙江大学 | 二次再热机组汽轮机调速系统仿真模型及其建模方法 |
| CN106096155A (zh) * | 2016-06-17 | 2016-11-09 | 华电电力科学研究院 | 一种未来超大型汽轮机组的调速系统研究方法 |
-
2017
- 2017-06-22 CN CN201710479894.0A patent/CN107103167B/zh not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH08218897A (ja) * | 1995-02-17 | 1996-08-27 | Toshiba Corp | ガスタービン制御装置 |
| CN1959067A (zh) * | 2005-11-03 | 2007-05-09 | 沈阳创思达自动化系统有限公司 | 汽轮机的调速控制系统 |
| CN103558031A (zh) * | 2013-11-11 | 2014-02-05 | 沈阳黎明航空发动机(集团)有限责任公司 | 一种燃气轮机甩负荷试验方法 |
| CN105700380A (zh) * | 2016-01-24 | 2016-06-22 | 浙江大学 | 二次再热机组汽轮机调速系统仿真模型及其建模方法 |
| CN106096155A (zh) * | 2016-06-17 | 2016-11-09 | 华电电力科学研究院 | 一种未来超大型汽轮机组的调速系统研究方法 |
Non-Patent Citations (3)
| Title |
|---|
| Load Shedding Via PID Controller for an Isolated Power System;D. Tyagi 等;《2012 IEEE Fifth Power India Conference》;20130314;第1-6页 * |
| Simulation and Optimization of Load Shedding Scheme for Islanded Power System;Cuicui Wu 等;《2010 International Conference on Power System Technology》;20101213;第1-6页 * |
| 汽轮发电机组甩负荷转速飞升特性的理论与实验研究;赵志勇;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20140415(第4期);第C042-64页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107103167A (zh) | 2017-08-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107103167B (zh) | 针对甩负荷工况的deh调速功能诊断方法及系统 | |
| US6466858B1 (en) | Methods and apparatus for monitoring gas turbine engine operation | |
| CN111159844B (zh) | 一种电站燃气轮机排气温度的异常检测方法 | |
| CN112861425A (zh) | 结合机理与神经网络的双轴燃气轮机性能状态检测方法 | |
| Mehrpanahi et al. | A novel dynamic modeling of an industrial gas turbine using condition monitoring data | |
| Tsoutsanis et al. | Non-linear model calibration for off-design performance prediction of gas turbines with experimental data | |
| CN102588315A (zh) | 透平压缩机喘振的自动测试方法 | |
| Tsoutsanis et al. | Part-Load performance of gas turbines: Part II—multi-point Adaptation with compressor map generation and GA optimization | |
| Simani et al. | Adaptive fault–tolerant control design approach for a wind turbine benchmark | |
| CN105912878A (zh) | 基于热力模型与粒子群优化算法相结合的燃气轮机自适应气路部件性能诊断方法 | |
| Wang et al. | An adaptive condition monitoring method of wind turbines based on multivariate state estimation technique and continual learning | |
| CN106021757A (zh) | 基于灰色关联理论与热力模型相结合的混合型燃气轮机气路部件自适应性能诊断方法 | |
| US20160365736A1 (en) | Model-based control system and method for power production machinery | |
| Salamat | Gas path diagnostics for compressors | |
| Zhao et al. | Gas turbine exhaust system health management based on recurrent neural networks | |
| US9500136B2 (en) | Systems and methods for generating variable ramp rates for turbomachinery | |
| CN108426717A (zh) | 发动机试验参数异常状况自动判断方法及存储介质 | |
| Gayme et al. | Fault detection and diagnosis in turbine engines using fuzzy logic | |
| CN115903484A (zh) | 基于合作博弈的航空发动机多变量鲁棒控制器优化方法 | |
| CN106354999B (zh) | 一种火电机组负荷振荡与负荷突变故障的在线诊断方法 | |
| CN112051066B (zh) | 小汽轮机进汽调门流量特性测定方法、装置、设备及介质 | |
| CN116517867B (zh) | 压气机的喘振诊断和抑制方法、装置 | |
| Undrill et al. | Modeling of combined cycle plants in grid simulation studies | |
| Tsoutsanis et al. | Performance assessment of classical and fractional controllers for transient operation of gas turbine engines | |
| Gutierrez et al. | Smoothing Techniques for Real-Time Turbine Speed Sensors |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200403 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |