锅炉汽包水位模型辨识的方法技术

一种计算机应用技术领域的锅炉汽包水位模型辨识的方法。本发明专利技术在现有的化工锅炉汽包水位控制系统中添加带滞环的偏置继电辨识模块，通过对系统在继电特性作用下的周期振荡极限环进行分析，获得汽包水位模型中四个参数的辨识公式。本发明专利技术针对化工过程中常见的汽包水位问题，提出一种新型的在线精确继电辨识方法，能在不需要任何时滞和静态增益等先验信息的情况下，根据时域分析得到的准确解析表达式精确计算出四个待辨识参数，从而得到逆响应的时滞模型与实际系统相匹配。

Boiler drum water level model identification method

Method for identifying boiler drum water level model in computer application technology field. The bias relay identification module of the invention added with hysteresis in the available chemical boiler drum level control system, based on the periodic oscillation limit system under relay characteristic loop analysis, identification formula of four parameters in the model of drum level. The invention of drum level in chemical process common problems, this article puts forward a new online accurate relay identification method, can without any delay and static gain prior information, according to the time domain analysis of accurate analytical expressions are obtained accurately calculate the four parameters, so as to obtain the delay model and the actual system inverse in response to the match.

【技术实现步骤摘要】

本专利技术涉及的是一种计算机应用
的方法，具体是一种。
技术介绍
影响锅炉汽包水位的扰动主要有四个来源一是给水方面的扰动；二是蒸汽负荷的变化；三是燃料量的变化；四是汽包压力的变化。其中，影响汽包水位变化的因素主要是给水流量、蒸汽流量和燃料量，选定汽包水位作为输出变量y，给水流量作为控制变量u，则汽包水位对象在水流量作用下的动态特性，是由积分环节、逆响应环节、一阶惯性环节和时滞环节组成的，其传递函数模型G(s)可以表示成G(s)=Y(s)U(s)=Kp(-τ1s+1)(-τ2s+1)se-θs]]>这里Kp是静态增益，τ1和τ2是时间常数，θ是时滞。当前的突出问题是如何精确辨识出这四个模型参数。采用基于常规的辨识方法(例如阶跃响应法)时，不能同时辨识出所有模型参数，需要事先估算时滞和静态增益大小，而且辨识精度不高，特别是容易受到外界噪声干扰，这种不精确辨识偏差会导致据此设计的控制器失效，严重的会引起汽包破裂或爆炸，这对实际生产非常不利和有害。因此精确辨识上述一阶积分逆响应时滞过程是化工实践中很困难并亟待解决的问题，必须采用先进的辨识技术来解决。strm，K.J.等在1984最早提出了继电反馈辨识方法，该方法可保障被控过程的稳定闭环振荡响应，从而便捷地获取过程临界信息，为PID控制器参数整定提供基础，因此近年来已广泛应用于工业PID控制器参数自整定系统中。许多改进和扩展的继电反馈辨识方法相继被提出来。经对现有技术的文献检索发现，针对稳定时滞对象，Wang，Q.G.在文献Low-order modeling from relay feedback，(基于继电反馈的低阶建模方法，发表在Industrial & Engineering Chemistry Research，化工工业与工程研究，1997，36(2)，375-381.)中首次应用带滞环的偏置继电器获得更为精确的辨识结果。但由于锅炉汽包水位对象的动态特性比较特殊，既有积分和时滞环节，又必须考虑“虚假水位”现象所对应的逆响应环节。Wang，Q.G.的方法只适用于辨识包含三个参数的稳定时滞对象。Luyben，W.L.在文献Identification and tuning of integrating processes with deadtime andinverse response(具有时滞和反向响应的积分对象的辨识和整定，发表在Industrial & Engineering Chemistry Research，化工工业与工程研究，2003，42，3030-3035.)中提出了一种阶跃响应辨识法，但该方法不仅需要事先知道时滞，且在有干扰噪声的情况下辨识精度很差。如何设计继电反馈实验快速有效的在线辨识出四个模型参数，如何提高辨识精度，如何克服外界干扰和模型失配问题，从而缓解辨识误差带给系统的控制难题一直都没有得到彻底的解决。
技术实现思路
本专利技术的目的在于针对化工过程中常见的锅炉汽包水位控制系统中，现有辨识技术的不足，提出一种，能在不需要任何时滞和静态增益等先验信息的情况下，根据对周期振荡极限环的时域分析得到准确解析表达式，从而精确计算出待辨识参数，使逆响应时滞模型与实际系统相匹配，进而为设计汽包水位PID控制器提供基础。本专利技术是通过以下技术方案实现的，具体包括以下步骤步骤一偏置继电器输出值μ+，μ-和滞环ε已知参数送到存储单元RAM中。步骤二通过继电辨识测试，上位机记录辨识过程的输入输出函数u和y，读取多个稳定振荡极限环周期参数的平均值，即周期Pu1，Pu2，幅值Ap，Ad，时间间隔 和 将各参数送到存储单元RAM中。步骤三调用事先编制好的算法程序解析辨识出系统模型的四个参数，具体算法步骤如下A)读取内存中的μ+，μ-，ε，Pu1，Pu2， 和 代入方程(1)式，计算待辨识的时间常数τ2。-τ2ln(μ+(1-e-Pu1/τ2)-μ-(1-e-Pu2/τ2))=t~min-t~max---(1)]]> B)将τ2代入方程组(2)、(3)和(4)，读取内存中的Ap和Ad，计算待辨识的时间常数τ1。θ=t~min+τ2ln(μ+(1-e-(Pu1+Pu2)/τ2)2μ(τ1/τ2+1)(1-e-Pu2/τ2))---(2)]]>Kp=Ap-Ad(μ+Pu1+2μτ1)-τ2---(3)]]>-Kp(μ-Pu2+2μθ)-2Kpμ(τ1+τ2)=2ϵ---(4)]]>C)将τ1和τ2代入方程(2)计算出待辨识参数θ。D)将τ1和τ2代入方程(3)计算出待辨识参数Kp。在以上各式中，Kp——对象静态增益，τ1——对象时间常数一，τ2——对象时间常数二，θ——对象时滞，μ+——继电器输出高电平，μ-——继电器输出低电平，ε——偏置继电器的滞环，Pu1——继电器输出高电平的维持时间，Pu2——继电器输出低电平的维持时间，Ap——对象输出振荡曲线的正峰值，Ad——对象输出振荡曲线的负峰值， ——从继电器跳变到对象达到正峰值的时间间隔， ——从继电器跳变到对象达到负峰值的时间间隔。本专利技术在现有的化工锅炉汽包水位控制系统中添加带滞环的偏置继电辨识模块，通过对系统在继电特性作用下的周期振荡极限环进行分析，获得汽包水位模型中四个参数的辨识公式。本专利技术在现有工控计算机上可以直接运行实施，在线辨识出锅炉汽包水位模型。本专利技术应用过程可以大致分为两个阶段，第一，锅炉汽包水位调试控制器参数初期，在继电器作用下，整个闭环系统进入稳定振荡过程。在此过程中，分别记录下继电器和对象输出数组，并保存在RAM中，通过简单的数据分析，得到多个关键的量测信息。第二，在振荡周期的多个关键量测信息已经存入工控机实时数据库后，调用事先编制好的程序解析辨识出系统模型参数。此算法正是在精确时域分析的基础上创新得出的。本专利技术全套调节过程可以在工控机仿真控制界面上完成。在实际化工过程控制现场采用本专利技术方法，只需单次继电辨识实验即可根据以上解析方法精确地得到四个参数，因而操作简便迅速，并能有效克服噪声干扰，从而克服了传统辨识方法的主要缺点。本专利技术给出的辨识方法可广泛应用于能源、冶金、石化、火力发电、纺织等行业中带时滞的逆响应生产过程的辨识中。附图说明图1为本专利技术继电辨识组成原理图；图2为本专利技术所应用带滞环偏置继电器的特性；图3为系统及继电特性输出波形；图4为工控机仿真控制界面；图5为辨识算法程序流程图；图6为具体实施例子应用本专利技术继电辨识方法与现有其他辨识方法辨识结果比较图。具体实施例方式以下结合附图和实施例对本专利技术的技术方案作进一步描述。实施例对于一个化工锅炉汽包水位对象，在前向给水流量在水流量作用下本文档来自技高网...

【技术保护点】
一种锅炉汽包水位模型辨识的方法，其特征在于，包括如下步骤：步骤一：偏置继电器输出值μ↓［＋］，μ↓［－］和滞环ε已知参数送到存储单元ＲＡＭ中；步骤二：通过继电辨识测试，上位机记录辨识过程的输入输出函数ｕ和ｙ，读取多个稳定振荡 极限环周期参数的平均值，即周期Ｐ↓［ｕ１］，Ｐ↓［ｕ２］，幅值Ａ↓［ｐ］，Ａ↓［ｄ］，时间间隔＊↓［ｍｉｎ］和＊↓［ｍａｘ］，将各参数送到存储单元ＲＡＭ中；步骤三：调用事先编制好的算法程序解析辨识出系统模型的四个参数，具体算法步骤如 下：Ａ）读取内存中的μ↓［＋］，μ↓［－］，ε，Ｐ↓［ｕ１］，Ｐ↓［ｕ２］，＊↓［ｍｉｎ］和＊↓［ｍａｘ］代入方程（１）式，计算待辨识的时间常数τ↓［２］：－τ↓［２］ｌｎ（μ↓［＋］（１－ｅ↑［－Ｐ↓［ｕ１］／τ↓［２］］ ）／－μ↓［－］（１－ｅ↑［－Ｐ↓［ｕ２］／τ↓［２］］））＝＊↓［ｍｉｎ］－＊↓［ｍａｘ］（１）Ｂ）将τ↓［２］代入方程组（２）、（３）和（４），读取内存中的Ａ↓［ｐ］和Ａ↓［ｄ］，计算待辨识的时间常数τ↓［１］：θ＝ ＊↓［ｍｉｎ］＋τ↓［２］ｌｎ（μ↓［＋］（１－ｅ↑［－（Ｐ↓［ｕ１］＋Ｐ↓［ｕ２］）／τ↓［２］］）／２μ（τ↓［１］／τ↓［２］＋１）（１－ｅ↑［－Ｐ↓［ｕ２］／τ↓［２］］））（２）＊＊＊（３）＊＊＊（４） Ｃ）将τ↓［１］和τ↓［２］代入方程（２）计算出待辨识参数θ；Ｄ）将τ↓［１］和τ↓［２］代入方程（３）计算出待辨识参数Ｋ↓［ｐ］；在以上各式中，Ｋ↓［ｐ］－对象静态增益，τ↓［１］－对象时间常数一，τ↓［２ ］－对象时间常数二，θ－对象时滞，μ↓［＋］－继电器输出高电平，μ↓［－］－继电器输出低电平，ε－偏置继电器的滞环，Ｐ↓［ｕ１］－继电器输出高电平的维持时间，Ｐ↓［ｕ２］－继电器输出低电平的维持时间，Ａ↓［ｐ］－对象 输出振荡曲线的正峰值，Ａ↓［ｄ］－对象输出振荡曲线的负峰值，＊↓［ｍａｘ］－从继电器跳变到对象达到正峰值的时间间隔，＊↓［ｍｉｎ］－从继电器跳变到对象达到负峰值的时间间隔。...

【技术特征摘要】

【专利技术属性】
技术研发人员：顾诞英，张卫东，王萍，蔡云泽，欧林林，
申请(专利权)人：上海交通大学，
类型：发明
国别省市：31[中国|上海]