区间二型模糊免疫PID在环己烷无催化氧化温度控制系统中的应用

区间二型模糊免疫PID在环己烷无催化氧化温度控制系统中的应用化工生产在国民经济生产活动中占有较大比重，但由于其规模大、工艺复杂、生产原料及产品具有易燃易爆性等特点，其生产过程具有高度危险性。一旦发生化工生产事故，将对国民经济、人员安全和周边环境造成巨大损失。环己烷无催化氧化过程是一个容易发生危险事故的复杂反应过程。环己烷无催化氧化工艺流程如图1所示，整个反应器由五个气升式环流反应器串联组成，首先苯和氢气在催化剂的助推下发生气相加氢反应并生产出高质量的液态环己烷。然后液相环己烷经过换热器加热到160℃左右由第一个反应器的底部自下而上进入反应器，按顺序依次经过剩下的四个反应器。氧气、空气和氮气组成含有一定量氧气的空气先进入空气分布器，通过空气分布器的调节作用使其自下而上流入各个反应器，同时推动反应器底部的液体向上流动，在反应器上部发生化学反应，反应结束后生成的反应物在反应器顶部进行气液分离。这时未参加反应的氧气、环己烷蒸气，以及产生的一氧化碳和二氧化碳气体进入尾气管道，将气态环己烷冷凝分离，而液相则流入下一个反应器底部继续反应。反应器内有夹套和蛇形管。反应过程中，蒸汽流经反应器内的夹套进行加热，蛇形管内有循环冷却水对反应器进行冷却[1-4]。

图1

图1环己烷无催化氧化工艺流程图

Fig.1Flowchartofuncatalysedoxidationofcyclohexane

环己烷无催化氧化过程是一个放热过程，反应温度一般在430~460℃之间。反应温度过高，中间产物环己基过氧化物的热不稳定性增强，副产物产量增加；如果反应温度过低，反应器尾气中的氧气浓度会升高，给生产过程带来很大的危险。此外，环己烷无催化氧化是一个复杂化学过程，具有非线性、多变量耦合和大时滞的特点，这使得传统的控制方案难以达到预期的效果，一旦反应温度失控会带来严重的安全隐患。因此，采用智能控制方案来进行复杂化工过程的控制是一个趋势。

近年来，区间二型模糊逻辑系统（intervaltype-2fuzzylogicsystem，IT2FLS）在过程控制等领域得到了广泛应用[5-8]。IT2FLS在模糊化及模糊推理过程中使用了区间二型模糊集合（intervaltype-2fuzzyset，IT2FS），它通过不确定性足迹（footprintofuncertainty，FOU）扩展了维数运算，使IT2FLS扩展了系统设计的自由度和模糊推理的能力，其处理不确定性的能力也有所提高[9-11]。Zhao等[12]设计了一种区间二型模糊逻辑控制器（intervaltype-2fuzzylogiccontroller，IT2FLC），并将其成功地应用到了双容水箱液位控制中，验证了IT2FLC的抗干扰性能。Galluzzo等[13]设计了一种IT2FLC，并将其成功地应用到了连续生物反应器混合废物的生物降解问题中，将仿真结果与传统的模糊控制器进行了比较。王永富等[14]提出了一种自适应IT2FLC，应用到了质子交换膜燃料电池的过氧比控制问题中，仿真结果显示提出的控制器具有良好的鲁棒性、跟踪性能及自适应特性。Kumbasar等[15]设计了一种基于IT2FLS的逆控制器，利用IT2FLS来逼近过程模型，最后将提出的控制器应用到了pH中和过程问题中，取得了比基于传统模糊系统的逆控制器更好的控制效果。Miccio等[16]提出了一种区间二型模糊PID控制器并成功地用于精馏塔的控制问题中，将仿真结果与模糊PID控制器进行了比较，取得了良好的控制效果。Han等[17]提出一种二型模糊广义学习控制器，该控制器可以通过自适应学习律在线调整参数，提升了控制器的鲁棒性，最后应用到了污水处理过程溶解氧浓度控制问题中，取得了理想的控制效果。何青等[18]提出了一种基于二型模糊变积分的PID控制器并用于无刷直流电机的转速控制，仿真结果表明该方法具有良好的鲁棒性能。李天等[19]提出了一种基于区间二型模糊逆模型的控制方案，并在乙烯裂解炉出口温度系统中进行了验证，实现了对出口温度平稳高效的控制。Liao等[20]提出了一种基于二型T-S模糊模型的分散控制方案，可以把多输入多输出的过程转化为多个独立的单回路进行控制器设计，并且无须知道系统的精确数学模型，利用二型T-S模糊模型逼近被控过程，构建了模糊鲁棒控制器，最后在多级蒸发器制冷系统中验证了该方法的有效性。

免疫是生物体体内的一系列特性生理反应组成的防御系统，通过协调主反馈机理与抑制机理的互动作用，完成对外来侵犯的保护与免疫。生物免疫系统是一个复杂系统，具有显著抵抗抗原的适应性能力，在不确定及强干扰的环境中具有较强的适应性和鲁棒性。生物免疫系统所具备的这些智能行为，从理论和技术上为解决一些科学和工程领域问题提供了一种新的思路。王玉勤等[21]提出了一种免疫PID控制算法，并成功地将其应用到了吊车-双摆控制系统，实验结果与线性二次型最优控制器相比，稳定性与鲁棒性都更优。李琦等[22]提出了一种改进免疫PID算法，利用遗传算法对免疫PID中的参数进行离线优化，并成功地应用到了双容水箱液位控制中，取得了预期的控制效果。任重昕等[23]将传统PID控制器与免疫控制器串联起来，设计了一种免疫非线性PID控制器，并成功地将其应用在了热磨机料位控制系统。Peng等[24]提出了一种串级免疫PID控制系统并应用到了主蒸汽温度控制系统中。Wang等[25]基于免疫反馈规则设计了免疫PID控制方案并将其成功地应用于干燥器温度控制系统。

本文结合免疫PID算法和IT2FLS的优点，提出一种区间二型模糊免疫PID（intervaltype-2fuzzyimmunePID，IT2FIPID）控制器。该控制器本质上是一种基于免疫PID的非线性控制器，利用IT2FLS的万能逼近特性来逼近免疫反馈控制律中的非线性函数，这样可以有效地处理具有时变和大延迟的复杂非线性系统。最后，将所提出的IT2FIPID控制器应用于环己烷无催化氧化温度控制系统。

1区间二型模糊免疫PID系统的结构

1.1免疫PID控制器

免疫是生物体特有的一种生理现象。生物免疫系统可以产生相应的抗体来抵抗外部入侵的抗原。当抗体与抗原相遇结合后，会激发一系列的反应机制，这时抗体会产生一种特殊酶来破坏抗原或者直接通过吞噬作用消灭抗原。淋巴细胞和抗体分子构成了生物免疫系统，淋巴细胞包括胸腺产生的T细胞（主要是辅助细胞TH和抑制细胞TS）以及骨髓产生的B细胞。当抗原进入生物体后，辅助细胞TH和抑制细胞TS收到信息，B细胞激活产生抗体以抵御抗原的入侵。基于上述生物免疫原理，有学者提出了免疫PID控制器结构。如果以ε(k)表示第k代抗原浓度，TH(k)和TS(k)分别表示辅助细胞TH和抑制细胞TS的浓度，则对B细胞的刺激S(k)表示如式(1)所示：

S(k)=TH(k)-TS(k)(1)

其中，TH(k)=K1ε(k)，TS(k)=K2f[S(k-d),ΔS(k-d)]ε(k)，d表示纯滞后时间，f([])为一个非线性激活函数，用于表示第k-d代抗体及抗体变化和抗原浓度对抑制细胞TS的影响。

这里偏差e(k)由抗原浓度ε(k)表示，控制输入u(k)由传入B细胞的刺激S(k)表示，免疫反馈控制律则可用式(2)表示：

u(k)=K1e(k)-K2f[u(k-d),Δu(k-d)]e(k)=K{1-ηf[u(k-d),Δu(k-d)]}e(k)(2)

其中，反应速度由K=K1来控制，系统的稳定性由η=K2/K1来决定[26]。

由式(2)可以看出，基于免疫反馈原理的控制实际上是一个非线性P控制器，通过与常规PID控制相结合进一步改进控制器的性能。给出常规PID控制算法如式(3)所示：

uPID(k)=Kp1+Kiz-1+Kdz-1ze(k)(3)

结合免疫控制律式(2)与常规PID控制算法式(3)，可以得到免疫PID控制器的输出如式(4)所示：

u(k)=u(k-1)+K{1-ηf[u(k-d),Δu(k-d)]}×1+kiz-1+kdz-1ze(k)(4)

其中，比例系数Kp=K1-ηf[u(k-d),Δu(k-d)]随着控制器输出的变化而变化，K表示增益。当满足0≤ηf[u(k-d),Δu(k-d)]≤1时，系统为负反馈控制。当满足η=0时，相当于传统的增量式PID控制算法。很明显，参数K、ki、kd、η以及非线性函数f()的选择对控制器的性能有重大的影响。目前，非线性函数f()的选择也是一个难点问题，还没有统一的标准或者规则用来指导f()的选择，大部分学者都选择神经网络或者模糊逻辑系统（FLS）来构造f()，两者都可以实现逼近任意的非线性函数。

1.2区间二型模糊免疫PID（IT2FIPID）系统

FLS的优点是将技术人员或专家的经验转化为规则，系统可以根据规则进行推理和决策，实现自动控制、分类、辨识等特定任务。IT2FLS不仅具有FLS的优点，而且在处理不确定性和复杂非线性问题方面更有优势。免疫PID具有结构简单、易于实现、鲁棒性和自适应能力强等特点。本文利用IT2FLS具备的万能逼近特性来逼近免疫PID系统的非线性函数f()，结合IT2FLS和免疫PID的优点，提出了一种IT2FIPID控制器，它的控制结构如图2所示。

图2

图2IT2FIPID控制系统框图

Fig.2BlockdiagramofIT2FIPIDcontrolsystem

IT2FLS由五部分组成：模糊化（fuzzification）、规则库（rulebase）、推理机（inferenceengine）、降型（reduction）和解模糊（defuzzifier），如图3所示。与传统的FLS相比，由于在模糊化和模糊推理过程中使用了IT2FS，因此，在得到系统的最终输出之前，必须经历将IT2FS转化为普通集合的过程，即IT2FLS的降型过程[27-30]。下面具体给出IT2FLS各部分的功能与作用。

图3

图3区间二型模糊系统的结构框图

Fig.3TheblockdiagramofIT2FLS

模糊化：本文中的IT2FLS是一个两输入单输出的系统。输入变量是控制变量u(k-d)和它的变化量Δu(k-d)，输出变量是免疫反馈律中的调节函数f()。每个输入变量分别划分为正（P）和负（N）两个模糊子集，而输出变量划分为正（P）、零（Z）和负（N）三个模糊子集。输入变量u(k-d)和Δu(k-d)以及输出变量f()均采用IT2FS表示，它们的隶属函数分别如图4~图6所示，图中的阴影区域为IT2FS的FOU，它由上、下隶属度函数包围组成。

图4

图4控制变量的隶属函数

Fig.4Membershipfunctionofcontrolvariable

图5

图5控制变量变化的隶属函数

Fig.5Membershipfunctionofcontrolvariablechange

图6

图6非线性函数的隶属函数

Fig.6Membershipfunctionofnonlinearfunction

规则库：本文中IT2FLS规则的确定与细胞接受的刺激有关，刺激大则抑制能力弱，刺激小则抑制能力强。最终确定IT2FLS的规则库由4条模糊规则组成：

(1)Ifu(k-d)isPandΔu(k-d)isPthenf[u(k-d),Δu(k-d)]isN

(2)Ifu(k-d)isPandΔu(k-d)isNthenf[u(k-d),Δu(k-d)]isZ

(3)Ifu(k-d)isNandΔu(k-d)isPthenf[u(k-d),Δu(k-d)]isZ

(4)Ifu(k-d)isNandΔu(k-d)isNthenf[u(k-d),Δu(k-d)]isP

推理机与降型：推理机的作用是将输入变量的模糊集合通过与规则表示的模糊关系进行合成运算。这里IT2FLS的输入变量首先经过模糊化的作用转化为IT2FS，然后与模糊规则前件作用产生激励强度，最后与模糊规则后件部分作用产生输出，得到的也是一个IT2FS，再经过降型过程处理得到普通模糊集合。目前，最为常见的中心集降型法的计算如式(5)所示：

Ycos(x)=∪fn∈Fn(x)yn∈Yn∑n=1Nfnyn∑n=1Nfn=[yl,yr](5)

其中，Fn(x)=[f̲n,f¯n]表示模糊规则前件激励强度；Yn=[y̲n,y¯n]表示模糊规则后件模糊集合的质心；Ycos(x)表示降型过程的输出；yl和yr表示降型集的两个端点值，具体的计算如式(6)、式(7)所示：

yl=mink∈[1,N-1]∑n=1kf¯ny̲n+∑n=k+1Nf̲ny̲n∑n=1kf¯n+∑n=k+1Nf̲n=∑n=1Lf¯ny̲n+∑n=L+1Nf̲ny̲n∑n=1Lf¯n+∑n=L+1Nf̲n(6)yr=maxk∈[1,N-1]∑n=1kf̲ny¯n+∑n=k+1Nf¯ny¯n∑n=1kf̲n+∑n=k+1Nf¯n=∑n=1Rf̲ny¯n+∑n=R+1Nf¯ny¯n∑n=1Rf̲n+∑n=R+1Nf¯n(7)

式中的转折点L和R必须满足下面的约束条件：

y̲L≤yl≤y̲L+1(8)y¯R≤yr≤y¯R+1(9)

IT2FLS降型的目的就是计算得到yl和yr。目前，学者们提出了很多种降型算法，其中KM算法是最常见、使用最普遍的降型算法，但是其存在计算效率低、不利于实时应用的缺陷[30]。本文采用EKM降型算法，是对KM算法的改进，通过对初始化方法、计算迭代及迭代停止条件的优化，不但可以保持和KM算法一致的计算精度，还可以至少节省39%的计算成本，有利于IT2FLS的实时应用[29]。EKM降型算法[31]计算yl的步骤如下：

（1）将y̲n按照升序排列，得到y̲1≤y̲2≤⋯≤y̲N，令k=[N/2.4]（接近N/2.4的整数）；

（2）计算a=∑n=1ky̲nf¯n+∑n=k+1Ny̲nf̲n和b=∑n=1kf¯n+∑n=k+1Nf̲n；

（3）找到k∈[1,N-1]，使其满足y̲k≤a/b≤y̲k+1；

（4）当k≠k时，计算s=sign(k-k)；

（5）更新a=a+s∑n=min(k,k)+1max(k,k)y̲n(f¯n-f̲n)和b=b+s∑n=min(k,k)+1max(k,k)(f¯n-f̲n)；

（6）更新k=k，找到k，使其满足y̲k≤a/b≤y̲k+1；

（7）令a=a，b=b，然后返回步骤（3）。

yr的计算步骤与yl相似，这里不详细介绍。

解模糊：降型过程的输出Ycos(x)是一个普通模糊集合，IT2FLS的最终输出一般取降型集两个端点的均值即可，如式(10)所示：

f[u(k-d),Δu(k-d)]=(yl+yr)/2(10)

2仿真研究

2.1环己烷无催化氧化温度控制系统的数学模型

环己烷无催化氧化过程是一个具有非线性、时变、大滞后的复杂化学变化。许多复杂化工过程都可以用一个一阶惯性环节加上纯滞后来近似表示。这里，采用阶跃测试的方法，通过对某石化公司尼龙盐厂现场数据的测试与分析得到环己烷无催化氧化温度控制系统的数学模型。首先，根据环己烷无催化氧化过程工艺特点，反应温度主要由蒸汽流量和冷却水流量控制，冷却水的流量相对于蒸汽流量受温度、压力等因素影响小，故本文中将冷却水流量近似认为恒定，所以整个控制系统的操纵变量为蒸汽流量，被控变量为反应温度。其次，在做阶跃测试时，具体施加的阶跃信号幅度，将充分考虑工艺操作要求及装置的实际生产情况，并与相关工艺操作人员协商后确定。在不影响正常生产且保障安全的条件下，施加的阶跃信号幅度应尽可能地激励出过程的动态和稳态特性，一般流量的阶跃信号幅度应保持在量程的5%~10%，并且在测试过程中每次施加的测试信号包括一个正向阶跃和一个负向阶跃。最后，根据系统产生的响应以及输入输出数据，对其进行预平滑和滤波处理，利用系统辨识的方法进行模型拟合。为了保证模型的精度和可靠性，在测试过程中，除非发生工艺参数超出安全范围的情况，否则测试人员不准干扰测试过程，一旦遇到各种未知干扰的影响，要对实验数据做详细记录，为删除不合理的数据提供依据，然后对不满足要求的测试结果，进行多次的阶跃测试，直至取得满意的测试结果。最终得到环己烷无催化氧化温度控制系统的数学模型如式（11）所示：

G(s)=0.9e-30s150s+1(11)

2.2仿真及性能比较

环己烷无催化氧化过程的正常反应温度在430~460℃之间。当反应温度低于380℃时，尾氧浓度会超过5%，达到高危水平，如果反应温度过高，环己烷选择性降低，副产物增多，达到反应器内环己烷的自燃点会造成严重事故。因此，本文中温度控制系统的设定值设置在440℃。仿真硬件环境为联想ThinkpadT480计算机，操作系统为Windows10，处理器为IntelCorei5-8250U（1.6GHz），运行内存8GB；仿真软件使用MATLAB，其版本为R2022b。为了公平地与其他方法进行比较，引入误差绝对值积分（integralofabsoluteerror，IAE）和误差平方积分（integralofsquareerror，ISE）2个性能指标，具体计算公式如式（12）和式（13）所示：

IAE=∫0∞e(t)dt(12)ISE=∫0∞e2(t)dt(13)

式中，e(t)为控制器输出与设定值之间的偏差。

提出的IT2FIPID控制器中，参数Kp=K1-ηf[u(k-d),Δu(k-d)]中的K和η的初值分别设置为0.05和0.8，ki和kd的初值分别设置为0.4和0.18。环己烷无催化氧化温度控制系统的仿真结果如图7所示，控制误差如图8所示，参数Kp的自整定结果如图9所示。

图7

图7温度控制系统仿真结果

Fig.7Simulationresultsoftemperaturecontrolsystem

图8

图8控制误差曲线

Fig.8Curvesofc