发明名称 燃煤锅炉炉膛压力系统混合控制方法
摘要 本发明涉及一种燃煤锅炉炉膛压力系统混合控制方法。本发明方法首先基于燃煤锅炉炉膛压力实时过程数据建立过程模型,挖掘出基本的过程特性;然后基于该过程模型建立比例积分微分控制回路;最后通过计算预测PI控制器的参数,将比例积分微分控制与燃煤锅炉炉膛压力对象整体实施预测PI控制。本发明方法弥补了传统控制的不足,并有效地方便了控制器的设计,保证控制性能的提升,同时满足给定的生产性能指标。本发明提出的控制技术可以有效减少理想炉膛压力工艺参数与实际炉膛压力工艺参数之间的误差,进一步弥补了传统控制器的不足,同时保证控制装置操作在最佳状态,使生产过程的炉膛压力工艺参数达到严格控制。
申请公布号 CN101709863B 申请公布日期 2012.01.11
申请号 CN200910155792.9 申请日期 2009.12.18
申请人 杭州电子科技大学 发明人 张日东;薛安克;葛铭;王俊宏;李春富
分类号 F22B35/00(2006.01)I 主分类号 F22B35/00(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 杜军
主权项 1.燃煤锅炉炉膛压力系统混合控制方法,其特征在于该方法包括以下步骤:(1)利用燃煤锅炉炉膛压力实时过程数据建立过程模型,具体方法是:首先建立燃煤锅炉炉膛压力实时运行数据库,通过数据采集装置采集N组实时过程运行数据,将采集的实时过程运行数据作为数据驱动的样本集合,表示为{x<sub>i</sub>,y(i)}<sub>i=1</sub><sup>N</sup>,i=1,2,…,N,其中x<sub>i</sub>表示第i组工艺参数的输入数据,y(i)表示第i组工艺参数的输出值;然后以该炉膛压力实时过程运行数据集合为基础建立基于最小二乘法的离散差分方程形式的局部受控自回归滑动平均模型:<img file="F2009101557929C00011.GIF" wi="1016" he="61" />X=[y(k-1),…,y(k-n),u(k-d-1),…,u(k-d-m)]<sup>T</sup>其中y<sub>L</sub>(k)表示当前时刻过程模型的工艺参数的输出值,x表示过程模型的工艺参数的过去时刻的输入和输出数据的集合,u(k)表示当前过程模型工艺参数对应的控制变量,k为当前的递推步数,Φ表示通过辨识得到的模型参数的集合,T表示矩阵的转置,n,m,d+1分别为对应实际过程的输出变量阶次、输入变量阶次、实际过程的时滞;采用的辨识手段为:<img file="F2009101557929C00013.GIF" wi="643" he="59" /><img file="F2009101557929C00014.GIF" wi="783" he="59" /><img file="F2009101557929C00015.GIF" wi="625" he="61" />其中k和P为辨识中的两个矩阵,<img file="F2009101557929C00016.GIF" wi="364" he="54" />γ为遗忘因子,<img file="F2009101557929C00017.GIF" wi="19" he="41" />为单位矩阵;(2)采用典型的响应曲线法设计炉膛压力过程模型的比例积分微分控制器,具体方法是:a.将过程模型的比例积分微分控制器停留在手动操作状态,操作拨盘使其输出有阶跃变化,由记录仪表记录过程模型的输出值,将过程模型输出值y<sub>L</sub>(k)的响应曲线转换成无量纲形式y<sub>L</sub><sup>*</sup>(k),具体是:<img file="F2009101557929C00018.GIF" wi="453" he="58" />其中,y<sub>L</sub>(∞)是过程模型的比例积分微分控制器的输出有阶跃变化时的过程模型输出y<sub>L</sub>(k)的稳态值;b.选取满足<img file="F2009101557929C00019.GIF" wi="578" he="60" />的两个计算点k<sub>1</sub>和k<sub>2</sub>,依据下式计算比例积分微分控制器所需要的参数K、T和τ:K=y<sub>L</sub>(∞)/qT=2(k<sub>1</sub>-k<sub>2</sub>)τ=2k<sub>1</sub>-k<sub>2</sub>其中q为过程模型的比例积分微分控制器输出的阶跃变化幅度; c.计算过程模型的比例积分微分控制器的参数,具体是:K<sub>c</sub>=1.2T/KτT<sub>i</sub>=2τT<sub>d</sub>=0.5τ其中K<sub>c</sub>为比例积分微分控制器的比例参数,T<sub>i</sub>为比例积分微分控制器的积分参数,T<sub>d</sub>分别为比例积分微分控制器的微分参数;(3)设计预测比例积分比例积分微分控制器,具体步骤是:d.将过程模型的比例积分微分控制器停留在自动操作状态,操作拨盘使其输入有阶跃变化,由记录仪表记录实时过程的输出,将过程输出值y(k)的响应曲线转换成无量纲形式y<sup>*</sup>(k),具体是:y<sup>*</sup>(k)=y(k)/y(∞)其中,y(∞)是过程模型的比例积分微分控制器的输入有阶跃变化时的过程模型输出y(k)的稳态值;e.选取满足y(k<sub>3</sub>)=0.39,y(k<sub>4</sub>)=0.63的另两个计算点k<sub>3</sub>和k<sub>4</sub>,依据下式计算预测比例积分比例积分微分控制器所需要的参数K<sub>1</sub>,T<sub>1</sub>和τ<sub>1</sub>:K<sub>1</sub>=y(∞)/q<sub>1</sub>T<sub>1</sub>=2(k<sub>3</sub>-k<sub>4</sub>)τ<sub>1</sub>=2k<sub>3</sub>-k<sub>4</sub>其中q<sub>1</sub>为过程模型的比例积分微分控制器输入的阶跃变化幅度;f.将步骤e得到的参数转化为拉普拉斯形式的局部受控传递函数模型:<img file="F2009101557929C00021.GIF" wi="381" he="125" />其中s为拉普拉斯变换算子,λ<sub>1</sub>为局部受控传递函数模型的时间常数,L<sub>1</sub>为局部受控传递函数模型的时滞,y(s)表示当前时刻过程模型的输出值的拉普拉斯变换,q<sub>1</sub>(s)表示过程模型的比例积分微分控制器输入的拉普拉斯变换;λ<sub>1</sub>=T<sub>1</sub>L<sub>1</sub>=τ<sub>1</sub>g.依据步骤f计算出的模型参数整定预测比例积分比例积分微分控制器的参数,具体方法是:①对该对象设计预测比例积分控制器;选期望的闭环传递函数模型为G<sub>q2</sub>(s)<img file="F2009101557929C00022.GIF" wi="444" he="131" />λ<sub>2</sub>为期望的闭环传递函数模型的时间常数,L<sub>2</sub>为期望的闭环传递函数模型的时滞,L<sub>2</sub>=L<sub>1</sub>;②预测比例积分比例积分微分控制器的传递函数G<sub>c1</sub>(s)可由下式来表示 <img file="F2009101557929C00031.GIF" wi="536" he="136" />③依据步骤②得到当前的预测比例积分比例积分微分控制器参数值u(s)<img file="F2009101557929C00032.GIF" wi="621" he="150" />。
地址 310018 浙江省杭州市江干区下沙高教园区2号大街