
在石灰窑炉的建模与控制中的应用
神经网络应用广泛,尤其在系统建模与控制方面,都有很好应用。下面简要介绍神经网络在石灰窑炉的建模与控制中的应用,以便更具体地了解神经网络在实际应用中的具体问题和应用效果。
1 石灰窑炉的生产过程和数学模型
石灰窑炉是造纸厂中一个回收设备,它可以使生产过程中所用的化工原料循环使用,从而降低生产成本并减少环境污染。其工作原理和过程如图1所示,它是一个长长的金属圆柱体,其轴线和水平面稍稍倾斜,并能绕轴线旋转,所以又称为迴转窑。含有大约30%水分的(碳酸钙)泥桨由左端输入迴转窑,由于窑的坡度和旋转作用,泥桨在炉内从左向右慢慢下滑。而燃料油和空气由右端喷入燃烧,形成气流由右向左流动,以使泥桨干燥、加热并发生分解反应。迴转窑从左到右可分为干燥段、加热段、煅烧段和泠却段。最终生成的石灰由右端输出,而废气由左端排出。
图1石灰窑炉示意图
这是一个连续的生产过程,原料和燃料不断输入,而产品和废气不断输出。在生产过程中首先要保证产品质量,包括的含量、粒度和多孔性等指标,因此必须使炉内有合适的温度分布,温度太低碳酸钙不能完全分解,会残留在产品中,温度过高又会造成生灰的多孔性能不好,费燃料又易损坏窑壁。但是在生产过程中原料成分、含水量、进料速度、燃油成分和炉窑转速等生产条件经常会发生变化,而且有些量和变化是无法实时量测的。在这种条件下,要做到稳定生产、高质量、低消耗和低污染,对自动控制提出了很高的要求。
以前曾有人分析窑炉内发生的物理-化学变化,并根据传热和传质过程来建立窑炉的数学模型,认为窑炉是一个分布参数的非线性动态系统,可以用二组偏微方程来描述,给出的数学模型是
(1)
(2)
在这些方程中,是固体的第个状态变量,是固体沿轴线的运动速度,是与空间、固体的状态变量和气体的状态变量有关的一个非线性函数。z代表沿轴方向的位置,t代表时间。而分别是气体的状态变量、速度和非线性函数。该系统具有分离的边界条件,也就是说,固体状态变量的值在入料处(冷端)是已知的,而气体状态变量的值在出料处(热端)是已知的,即已知
其中,L是窑长度。初始条件是系统在扰动前正常工作时状态变量的值,即
因为固体的状态变化慢而气体的状态变化快,可以忽略气体状态的变化,将式(2)简化为
这些方程中有很多参数,必须通过机理分析、假设或大量实验来确定,应用该数学模型需要测所有的状态变量,而这在实际情况下很难做到。
2石灰窑炉的神经网络建模
上面讲的是用机理分析建立数学模型的方法,它要搞清系统内部的物理-化学变化规律,并用严格的数学方程加以描述。从过程控制的角度来看,这种建模方法不仅是很难实现的,而且也不是十分必要的。因此大多采用系统辨识方法,将对象看作一个“黑箱”,不去分析其内部的反应机理,而只研究对象主要控制变量和输出变量之间的相互关系,神经网络系统辨识方法就是其中一种。
对石灰窑炉来说,主要的控制量有两个,一个是燃料流速,另一个是风量流速,这是生产中的主要调节手段。被控量选为炉窑热端的温度和炉窑冷端的温度,因为这两点的温度决定了炉内的温度分布曲线,而温度分布曲线又是影响产品质量和能耗的最关键的因素。从实现角度看,这四个变量也是容易实时测量的。因此,石灰窑可近似为一个二输入二输出的非线性动态系统,其中与是相关的。
应用前面介绍的神经网络辨识方法,石灰窑NARMA方程是
其中
采样周期定为1.125min,根据经验和试验结果,选择已能足够精确地描述系统动态特性。神经网络模型的结构与系统的结构相同,可以用下式表示:
定义系统输出与模型输出之间的误差矢量为
则系统辨识的指标函数是
神经网络选为有两个隐层的四层前馈网络,即隐层的节点数是根据试验确定的,隐层各神经元的激活函数均是。图2给出了用神经网络辨识石灰窑模型的系统结构。
图2石灰窑的神经网络系统辨识
训练和检验神经网络模型需要大量能充分反映系统非线性的特性的输入输出样本,为了不影响正常生产,这里输入输出样本是一个已被验证的机理模型上做仿真实验得到的。为了对系统充分激励,使训练能覆盖其全部工作范围,对系统分别输入正弦信号,阶跃信号和伪随机二进制信号,得到4000组输入输出数据。将这些数据归一化,使它们都处于[-1,1]范围内。再用BP学习算法训练神经网络,直至均方根误差定义为
训练以后神经网络输出结果如图3所示。
3石灰窑炉的神经网络控制
在石灰窑炉正向模型辨识的基础上,下面着重讨论石灰窑炉神经网络内模控制,包括逆模型的辨识、内模控制的设计和仿真等问题。
1)逆模型的辨识
由前分析,石灰窑炉可以看作二输入二输出的非线性动态系统(见石灰窑的NARMA模型)。则石灰窑的逆模型可用下面的方程描述
这里用来代替,因为的作用结果,在时刻尚不知道,表示神经网络已辨识的石灰窑炉的正模型输出结果。
图3相同输入下系统输出(实线)与神经网络输出(点线)
用间接法训练石灰窑炉逆模型的系统结构如图4所示。为了克服间接法可能会使系统不稳定的缺点,在此之前先用直接法训练逆模型,获得较好的初值,逆模型仍采用二个隐层的前馈网络,其结构为,激励函数仍为,学习算法也与训练正模型相同。
图4间接法训练石灰窑炉的逆模型
2)内模控制的设计
石灰窑炉的神经网络内模控制系统原理如图5所示。图中是神经网络逆模型,作为控制器,是神经网络正模型。滤波器是一个积分器,目的在于消除静态误差。
图5石灰窑炉的神经网络内模控制系统
为了测试石灰窑炉的内模控制效果,进行以下仿真实验。仿真所采用模型及神经网络训练样本的产生是基于一个被验证机理模型。首先让参考输入在工作点附近可以随机变化,观测输出的变化,实验结果如图6所示。图(a)中实线代表参考输入的增量,点线代表输出的增量(点线基本与实践重合)。图(b)是当参考输入变化时相应的的变化。由仿真结果可以看出,神经网络内模控制具有很好的跟踪性能。
图6神经网络内模控制的跟踪性能
为了进一步考查神经网络内模控制的鲁棒性,在正常工作状态下,保持参考输入不变,依次加入三种扰动:
(1)时,使进料速度增大10%。
(2)时,使进料中固体的含量由74%降为55%,即进料中水的含量增加。
(3)时,使窑炉转速降低5%。
内模控制系统抗扰性的实验结果示于图7,图7(a)是的变化曲线,图7(b)是的变化曲线,图7(c)是(粗线)和(细线)的变化曲线。在三种扰动下,系统的输出最终能恢复到正常值,这说明内模控制具有良好的抗扰性。
在计算机实现时,所有神经网络模型,包括正向模型和逆模型以及相应的学习算法都在计算机内部实现,并通过A/D转换从现场采集信息,通过D/A将神经网络产生的控制信号送到被控对象。
图7内模控制的抗扰性能
4离散论域模糊控制器设计
考虑九管还原炉的自动控温系统,其结构和控制任务描述如下:
钨的熔点为3000℃以上,目前尚不能采用通常的冶炼法,而只能用粉末冶金法进行加工。九管还原炉就是用来对氧化钨粉末还原去氧的装置。它先将三氧化钨还原成一氧化钨,此过程称为一次还原;然后再将一氧化钨还原成钨,此过程为二次还原。九管还炉的结构示意图如图8所示。其中(a)为正视图;图(b)为侧视图。在焙烧管道(共9根,上层5根,下层4根)的上、下部,分别装有电热丝。共6个温区,各温区的温度根据工艺要求必须稳定在给定值附近,偏差不能超过℃。每个温区的几何中心装有一个热电偶,以检测本温区的温度,每两个温区由一相电源供电加热。装有三氧化钨粉末的料舟,从管道的一端进入,经过各温区后从另一端出来,一次还原时,每30min出两舟,二次还原时,每15min出两舟。
图8 九管还原炉结构示意图
本系统原靠人工控温,后来改为双位开关式控制,其偏差波动很大,往往大于℃,影响钨粉产品质量。
自动控温的任务是,将各温区的温度稳定在给定值范围内(第一温区550℃,第二温区650℃……第六温区850℃),偏差不准大于±5℃。
由于要求九管还原炉采用恒值控温,各温区的温度偏移范围很小,而本温区受其他温区的影响虽然貌视很大,但这种影响只是一个恒定值,对系统的动态特性没有影响,对静态特性却有一定的偏移作用。如果将这一恒定的偏移视为本温区给定量的一部分(恒定分量),则六个温区就可分别视为六个的电炉。这样,被控对象就成了六个的单输入、单输出环节了。
由于确定九管还原炉精确数学模型非常困难,再加上各种非线性和分布参数的影响,使得传统控制方法难于满足要求。因此下面将采用离散论域的模糊逻辑进行控制,原理如图9所示。
图9 还原炉温度的模糊控制系统
图中,r为温度设定值,e=y-r为温度误差,是差分算子,是误差差分,是论域变换因子,u是控制量,它对应电加热器装置的晶闸管导通角的变化量。离散论域的模糊控制器主要设计如下:
1)输入/输出变量论域的离散化
e,和u都是实数上的连续变量,在该还原炉问题中它们变化范围分别为:。将 e变换到离散论域,上,将变换到,上,将u变换到上,得到离散论域上的输入/输出变量。
2) 定义模糊集合及其隶属函数表
对定义8个模糊集合,分别代表PL(正大)、PM(正中)PS(正小)、PZ(正零)、NZ(负零)、NS(负小)NM(负中)、NL(负大)。对定义7个模糊集合,分别代表PL,PM,PS,ZE,NS,NM,NL。对定义7个模糊子集,分别代表PL,PM,PS,ZE,NS,NM,NL。
对的模糊集合分别定义隶属函数,如表1,表2,表3所示。
表1 的隶属函数表
| 隶属度 | 的论域 | |
| -6 –5 –4 –3 –2 –1 –0 +0 1 2 3 4 5 6 | ||
| 模糊集合 | 0 0 0 0 0 0 0 0 0 0 0 0.2 0.7 1.0 0 0 0 0 0 0 0 0 0 0.2 0.7 1.0 0.7 0.2 0 0 0 0 0 0 0 0.1 0.7 1.0 0.7 0.1 0 0 0 0 0 0 0 0 0 1.0 0.7 0.1 0 0 0 0 0 0 0 0 0.1 0.7 1.0 0 0 0 0 0 0 0 0 0 0.1 0.7 1.0 0.7 0.1 0 0 0 0 0 0 0 0.2 0.7 1.0 0.7 0.2 0 0 0 0 0 0 0 0 0 1.0 0.7 0.2 0 0 0 0 0 0 0 0 0 0 0 | |
| 隶属度 | 的论域 | |
| -6 –5 –4 –3 –2 –1 0 1 2 3 4 5 6 | ||
| 模糊集合 | 0 0 0 0 0 0 0 0 0 0 0.2 0.7 1.0 0 0 0 0 0 0 0 0 0.2 0.8 1.0 0.8 0.2 0 0 0 0 0 0 0 0.8 1.0 0.8 0.2 0 0 0 0 0 0 0 0.5 1.0 0.5 0 0 0 0 0 0 0 0.2 0.8 1.0 0.8 0 0 0 0 0 0 0 0.2 0.8 1.0 0.8 0.2 0 0 0 0 0 0 0 0 1.0 0.7 0.2 0 0 0 0 0 0 0 0 0 0 | |
| 隶属度 | 的论域 | |
| -6 –5 –4 –3 –2 –1 0 1 2 3 4 5 6 | ||
| 模糊集合 | 0 0 0 0 0 0 0 0 0 0 0.2 0.7 1.0 0 0 0 0 0 0 0 0 0.2 0.8 1.0 0.8 0.2 0 0 0 0 0 0 0.1 0.8 1.0 0.8 0.1 0 0 0 0 0 0 0 0.5 1.0 0.5 0 0 0 0 0 0 0 0.1 0.8 1.0 0.8 0.1 0 0 0 0 0 0 0.2 0.8 1.0 0.8 0.2 0 0 0 0 0 0 0 0 1.0 0.7 0.2 0 0 0 0 0 0 0 0 0 0 | |
总结操作人员的经验,给出52条控制规则,归纳于表4。控制规则的一般形式为
表4 控制规则表
| NL | NM | NS | NZ | PZ | PS | PM | PL | ||
| PL PM PS ZE NS NM NL | PL | PM | NM | NM | NM | NL | NL | × | |
| PL | PM | NM | NM | NM | NS | NS | × | ||
| PL | PM | NS | NS | NS | NS | NM | NL | ||
| PL | PM | PS | Z | Z | NS | NM | NL | ||
| PL | PM | PS | PS | PS | PS | NM | NL | ||
| × | PL | PS | PS | PM | PM | NM | NL | ||
| × | PL | PL | PM | PM | PM | NM | NL | ||
4)求模糊关系矩阵和模糊控制表
由于论域是离散的,模糊控制规则集可以表示为一个模糊关系矩阵R
其中R是一个(14×13)×13的大矩阵,占2366个内存单元。
算出R以后,便可以对和设不同的值,例如,将它们看成是模糊单点,这时的输入是模糊矢量,可求出输出模糊矢量
其中是一个14×13的矩阵,将其按行排成矢量,再与R合成。由于此时只有第一个元素为1,其他元素皆为零,即是R的第一行。中的隶属度最大对应的即为此时输出。用同样方法,对每对输入,都可以求出相应的输出,将它们整理成模糊控制表,如表5所示。
表5 模糊控制查询表
| -6 | -5 | -4 | -3 | -2 | -1 | -0 | +01 | 1 | 2 | 3 | 4 | 5 | 6 | ||
| -6 | 6 | 5 | 6 | 5 | 3 | 3 | 3 | 3 | 2 | 1 | 0 | 0 | 0 | 0 | |
| -5 | 5 | 5 | 5 | 5 | 3 | 3 | 3 | 2 | 1 | 1 | 0 | 0 | 0 | 0 | |
| -4 | 6 | 5 | 6 | 5 | 3 | 3 | 3 | 3 | 2 | 1 | 0 | 0 | 0 | 0 | |
| -3 | 5 | 5 | 5 | 5 | 4 | 4 | 4 | 4 | 2 | -1 | -1 | -1 | -1 | -1 | |
| -2 | 6 | 5 | 6 | 5 | 3 | 3 | 1 | 1 | 0 | -2 | -3 | -3 | -3 | -3 | |
| -1 | 6 | 5 | 6 | 5 | 3 | 3 | 1 | 0 | 0 | -2 | -2 | -3 | -3 | -3 | |
| 0 | 6 | 5 | 6 | 5 | 3 | 1 | 0 | 0 | -1 | -3 | -5 | -6 | -5 | -6 | |
| 1 | 3 | 3 | 3 | 2 | 0 | 0 | 0 | -1 | -3 | -3 | -5 | -6 | -5 | -6 | |
| 2 | 3 | 3 | 3 | 1 | 0 | 0 | -1 | -1 | -3 | -3 | -5 | -6 | -5 | -6 | |
| 3 | 1 | 1 | 1 | 0 | 0 | 0 | -1 | -1 | -2 | -2 | -5 | -5 | -5 | -5 | |
| 4 | 0 | 0 | 0 | -1 | -1 | -2 | -3 | -3 | -3 | -3 | -5 | -6 | -5 | -6 | |
| 5 | 0 | 0 | 0 | -1 | -1 | -2 | -3 | -3 | -3 | -3 | -5 | -5 | -5 | -6 | |
| 6 | 0 | 0 | 0 | -1 | -1 | -1 | -3 | -3 | -3 | -3 | -5 | -5 | -5 | -6 | |
控制表可离线算出,在实时控制时得到和后,先变换成离散量和,查表得到,再论域反变换求出真正的连续控制量。
以上设计的模糊控制器已成功应用于还原炉的温度控制中,并取得了很好的控制效果。温度控制的上升时间短,超调量小、控制性能稳定、对环境变化不敏感,而且算法非常简单、执行快。图10和图11给出了模糊逻辑控制的仿真结果。
图10 双位开关式与Fuzzy控制特性比较 图11 时延对Fuzzy控制器特性的影响
1—Fuzzy 控制,2—双位开关式控制
图11 时延对Fuzzy控制器特性的影响
1—=0.7,2—=1.6
该系统实现的硬件结构如图12所示,其中单片采有8031,并扩充了8KEPROM。由此可见,该系统实现容易,且对系统配置在求不高。
图12 系统硬件结构图
