1、基于多模型切换控制的煤气化工业过程先进控制张天辰，赵众(北京化工大学信息科学与技术学院，北京100029)摘要：煤气化是煤化工生产过程的重要组成工艺。煤气化过程以煤和氧气作为主要反应原料，在高温高压条件下发生化学反应生成粗合成气。相较石油化工行业，煤气化工业过程的先进控制(APC)实施有较大难度，原因在于煤气化过程中煤质不断变化，导致传统的基于单一模型的先进控制方法会引起控制器模型失配问题，引发生产过程中关键参数(氧煤比、气化炉温度)的较大波动。针对上述问题，基于水煤浆气化炉的生产实际，提出了一种基于多模型切换的动态矩阵控制方法。在离线状态下，利用不同煤质工况下的过程数据，构造多工况模型集。在
2、线控制时，采取累积平方误差总平方波动(ISE-TSV)作为控制器性能指标表征控制器性能与模型失配程度，同时使用多模型预测值作为模型切换准则。通过多模型切换动态矩阵控制，实现了煤气化单元的先进控制。根据所提方法开发了多模型切换控制软件 Wisdom-Controller，在 UniSim 平台上进行了控制模拟，并进行了工业应用测试，模拟及工业应用测试结果证实了所提方法能够在煤质工况波动的情况下，实现氧煤比及气化炉温度的精准控制。对比传统的人工手动操作，基于多模型切换控制的先进控制方法降低了氧煤比及气化炉温度均方控制偏差。同时减少了比煤耗，增加了合成气体产量，显著提高了装置的经济效益。从而为实现煤
3、气化工业过程的先进控制提供了一条新的有效途径。关键词：煤气化；工况变化；多模型切换；动态矩阵控制；先进控制中图分类号：TQ546文献标志码：A文章编号：02539993(2023)04174712Advanced process control of coal gasification industrial process based onmultiple models switching controlZHANGTianchen,ZHAOZhong(College of Information Science and Technology,Beijing University of Chem
4、ical Technology,Beijing100029,China)Abstract:Coalgasificationisacrucialunitofcoalchemicalproductionprocess.Coalandoxygenareusedasthemainre-actionrawmaterialsinthecoalgasificationprocessandthecrudesyngasisgeneratedbychemicalreactionunderhightemperatureandhighpressure.Comparedwiththepetrochemicalindus
5、try,itismoredifficulttoimplementadvancedpro-cesscontrol(APC)inthecoalgasificationindustrialprocessbecauseofthetimevaryingdisturbanceofcoalquality.ThetimevaryingdisturbancemayleadtothemismatchoftheAPCmodelthatcouldcausesomelargefluctuationsofkeypro-cessindices(oxygencoalratioandgasifiertemperature).I
6、nviewofaboveproblems,basedontheactualproductionofcoal-waterslurrygasifier,adynamicmatrixcontrolmethodbasedonmultiplemodelsswitchingisproposedinthispaper.Theoff-lineprocessdataunderdifferentcoalqualityconditionsareusedtoconstructamultipleworkingmodelsetforon-linedynamicmatrixcontroller.Theintergralsq
7、uarederror-totalsquaredvariation(ISE-TSV)isusedasthecontrollerper-收稿日期：20220420修回日期：20220713责任编辑：钱小静DOI：10.13225/ki.jccs.2022.0549基金项目：北京市自然科学基金资助项目(4172044)作者简介：张天辰(1997)，男，福建福州人，硕士研究生。E-mail：通讯作者：赵众(1970)，男，河南新乡人，教授，博士生导师。E-mail：引用格式：张天辰，赵众.基于多模型切换控制的煤气化工业过程先进控制J.煤炭学报，2023，48(4)：17471758.ZHANGTian
8、chen，ZHAOZhong.AdvancedprocesscontrolofcoalgasificationindustrialprocessbasedonmultiplemodelsswitchingcontrolJ.JournalofChinaCoalSociety，2023，48(4)：17471758.第48卷第4期煤炭学报Vol.48No.42023年4月JOURNALOFCHINACOALSOCIETYApr.2023formanceindextomonitorthecontrollerperformanceandmodelmismatchandthemultiplemodelp
9、redictionvalueisusedasthemodelswitchingcriterion.Theadvancedcontrolofthecoalgasificationunitisrealizedthroughthemultiplemodelsswitchingdynamicmatrixcontrol.Accordingtotheproposedmethod,amultiplemodelswitchingcontrolsoft-wareWisdom-Controllerhasbeendeveloped.TheWisdom-ControllerhasbeentestedontheUniS
10、imsimulationplatformandappliedtotherealindustrialgasifier.Thesimulationandindustrialapplicationresultshaveverifiedthattheproposedmethodcanaccuratelycontrolthechangeofoxygen-coalratioandgasifiertemperatureundertheconditionoffluctuatingcoalqualityconditions.Comparedwiththetraditionalmanualoperation,th
11、emeansquarecontroldeviationofoxygen-coalratioandthegasifiertemperaturehavebeenreducedobviouslywiththeproposedadvancedprocesscontrolmethodbasedonmultiplemodelsswitching.Also,thespecificcoalconsumptionhasbeenreduced,thesyntheticgasoutputhasbeenin-creasedandtheeconomicbenefitoftheunithasbeensignificant
12、lyimproved.Theindustrialapplicationresultshaveveri-fiedthattheproposedmethodprovidesanewandeffectivewaytorealizetheadvancedprocesscontrolofcoalgasifica-tionindustrialprocess.Key words:coalgasification；workingconditionchange；multiplemodelsswitching；dynamicmatrixcontrol；ad-vancedprocesscontrol煤气化是一个热化
13、学加工过程，以煤和氧气作为主要反应原料，在高温高压条件下发生化学反应生成粗合成气1。气化过程伴随着煤质的多变，引起了生产的波动性，主要表现在关键变量氧煤比及气化炉温度的波动2-3。依托先进控制(APC)技术对气化系统进行先进控制和实时优化，结合稳态和动态模型预测生产过程等智能化方式4-5可以帮助操作人员快速响应煤质及负荷等工况切换，使煤气化操作始终保持接近最优操作点的“卡边”优化运行，能为企业创造更大的效益已成为煤化工企业的共识6。相较于石油化工行业，煤化工行业在工业智能化方面发展相对落后。尤其在煤气化工业领域，APC 的应用相对较少。2015 年，神华包头煤化工公司在聚乙烯和聚丙烯装置上引进
14、 APC 先进控制系统7，使工艺系统的自动化程度大大提高，但也指出了 APC 系统必须在相对稳定条件下使用，仍存在改进空间。2018年，江凤月等8将 APC 技术运用在煤制氢装置，通过单元优化煤耗降低了 1.83%。2019 年，晁澄9将 APC技术运用在 GE 单喷嘴水煤浆气化中，这次 APC 首次运用在水煤浆气化炉上，单位煤的有效气产量提高了0.24%。同年，张成学等10在四喷嘴水煤浆气化装置上运用 APC，有效减少了关键变量的波动幅度。2020 年，MIGHANIM 等11将 APC 应用于气化过程，减少了生产每单位有效合成气所需的氧气 1%以上、煤炭 3%以上。2021 年，方薪晖12
15、介绍了 APC 技术在日投料 3000t 的超大型水煤浆气化炉上的应用。工业上常用的先进控制算法主要是模型预测控制(MPC)13-14，其有较好的闭环稳定性和鲁棒性(抗干扰能力)。然而传统的基于线性模型的预测控制难以适应工况的大幅变化，即应用于煤气化过程，难以适应煤质的变化，存在一定的局限性。笔者提出了一种基于多模型切换的动态矩阵控制方法。在离线状态下，利用不同煤质工况下的过程数据，构造多工况动态模型集。在线控制时，采取累积平方误差总平方波动(ISE-TSV)作为控制性能指标15，用于判别煤质变化，同时使用多模型预测值作为模型切换准则。通过多模型切换动态矩阵控制，实现了煤气化单元的先进控制。所
16、提方法在 UniSim 平台上进行了模拟，并进行了工业应用测试，模拟及工业应用测试结果证实了所提方法能够在煤质工况波动的情况下，实现对煤气化炉关键变量的精准控制。1煤气化工艺及控制需求1.1煤气化工艺在气化炉内，煤中的碳元素主要发生如下反应16，其中目标产品是 CO 和 H2。C+O2 CO2(1)C+CO2 2CO(2)2C+O2 2CO(3)C+H2O CO+H2(4)C+2H2O CO2+2H2(5)C+2H2 CH4(6)2C+2H2O CH4+CO2(7)1.2水煤浆气化技术常见的水煤浆气化技术是德士古水煤浆加压气化技术17，该工艺采用水煤浆及氧气作为原料，对其1748煤炭学报202
17、3年第48卷进行加压处理，生成合成气。其工艺流程如图 1 所示。煤浆槽中的水煤浆通过高压煤浆泵加压，结合高压氧气。2 者按照一定比例混合进入气化炉上方的燃烧室，水煤浆被雾化成细小的颗粒。在高温高压条件下进行一系列反应生成粗合成气及少量残渣。反应生成的合成气进入在气化炉下部的激冷室内被激冷水冷却后从气化炉侧面的气体出口管处排出；熔渣沉降于激冷室底部，通过气化炉底部锁斗进行排渣，气化炉中下部排出的黑水排至灰水处理系统。合成气排出后进入碳洗塔内与洗涤水混合，进行洗涤除尘工艺，最终得到水煤浆气化工艺的目标产品。水煤浆氧气德士古合成气激冷水合成气碳洗塔气化炉洗涤水洗涤水洗涤水灰水黑水黑水排渣图1德士古煤
18、气化工艺流程Fig.1TechnologicalprocessofTexacocoalgasification1.3煤质对气化的影响煤质是影响煤气化安全性和操作平稳性的关键指标18-20。根据水煤浆气化的工艺特点，本文研究煤质的变化指标包括水分、固定碳成分、灰分、挥发分。(1)水分。水分主要考虑煤中的内在水分。在水煤浆制备过程中，水分影响煤浆浓度，煤的内在水分越低，煤浆浓度越高。煤浆浓度越高，进入煤气化炉内的水分越少，从而减少了蒸发水分所需的热量，气化炉温度也会有一定程度的变化。(2)灰分。煤的灰分是指煤完全燃烧后的固体残渣，一般来源于矿物质。煤的灰分越高，其发热量越低，不仅影响了气化炉的经济
19、效应，还增大了排灰的环境压力。低灰熔点的煤也会影响煤气化炉的渣口压差及气化炉的温度，对煤气化炉的安全造成较大影响。(3)挥发分。煤的挥发分是指煤样在隔绝空气条件下加热至一定温度后除掉外部水分，分解出的液体和气体产物。煤的挥发分影响煤气的热值，在一定程度上影响气化炉的温度。(4)固定碳。煤中的固定碳成分是去除水分、灰分、挥发分后剩余的残留物。煤的发热量主要由固定碳燃烧产生，固定碳成分变化对气化炉的温度将产生较大影响。1.4煤气化操作难点及先进控制器结构煤气化过程中，氧煤比作为关键操作变量显著影响了工艺的安全性与产品的经济性21-22。当前煤气化工艺下，氧煤比大都通过手动控制氧气流量以及煤浆流量来
20、实现其稳定23，但由于阀门的滞后性与精确性，氧气流量和煤浆流量往往波动较大，造成氧煤比随之波动。氧煤比不仅受装置影响，还与煤质特性有关，主要包括固定碳成分、灰分、挥发分。为了解决煤质对气化炉的影响，笔者提出了多模型切换动态矩阵控制的方法。煤气化炉系统是一个煤质多变的非线性对象，单一的模型难以概述这一非线性系统24。因此需要在不同煤质的工况下建立模型，利用动态矩阵控制方法实现多工况煤气化过程的多变量协调先进控制。根据水煤浆气化的生产实际与工业控制需求，设计的先进控制器结构见表 1。被控变量中，氧煤比和气化炉炉温是安全性约束指标，其中氧煤比的权重更大，被控变量权重越大，代表该变量对于工艺安全来说越
21、重要。干扰变量则选取煤质特性中不可测的固定碳成分、挥发分、灰分。表 1中的干扰变量煤质难以实时测量25，用离线化验分析值来表征工况的变化。表 1 煤气化炉先进控制器结构Table 1 Advanced controller structure of coal gasifier变量类型变量名称变量描述单位操纵变量(MV)MV1氧气流量Nm3/hMV2煤浆流量t/h被控变量(CV)CV1氧煤比Nm3/m3CV2气化炉温度干扰变量(DV)DV1煤质特性固定碳成分%DV2煤质特性挥发分%DV3煤质特性灰分%煤气化炉系统的控制目标为通过调整氧气回路和煤浆回路设定值，以保证氧煤比及气化炉温度稳定在所需的控
22、制区间内。2煤气化炉多模型切换控制常用的预测控制算法是动态矩阵控制算法，它是工业上应用最广的先进控制算法之一26，以稳态操作点附近的阶跃响应模型作为预测模型。煤气化炉系统作为一个非线性动态系统，具有明显的操作区域变化特征。预测控制方法具有前馈补偿性能，可以有效针对炉温热电偶的滞后特性。2.1多煤质模型的非线性动态表述考虑如图 2 所示的非线性动态系统，其中，Go为第4期张天辰等：基于多模型切换控制的煤气化工业过程先进控制1749rduye非线性动态对象；Gc为控制器模型；为系统设定值；为系统扰动；为控制器输出；为系统实际输出；为系统的误差。reuydGcGo+图2非线性动态模型Fig.2Non
23、lineardynamicmodel考虑线性时不变系统，即系统的参数不随时间变化，非线性对象及控制器的动态模型分别为Go:x(k+1)=f(x(k),u(k),d(k)y(k)=g(x(k),u(k)(8)Gc:w(k+1)=p(w(k),e(k)u(k)=q(w(k),e(k)(9)kx(k)Rn1w(k)Rm1f()g()p()q()式中，为系统的采样时刻；为非线性对象的状态变量；为控制器的状态变量；、为对应的非线性稳态函数。忽略闭环非线性系统的静态误差，即e=ry=0(10)(x0,u0,w0,d0,0)对给定的设定值 r0和扰动 d0，根据式(8)、(10)可以得到一个静态工作点满足f
24、(x0,u0,d0)=0(11)在静态工作点泰勒展开，可以得到在该点处线性化后的模型x(k+1)=fxx(k)+fuu(k)y(k)=gxx(k)+guu(k)(12)x(k)其中，为状态变量的变化量。同理得到静态工作点线性化后的控制器模型w(k+1)=pww(k)+pee(k)u(k)=qww(k)+qee(k)(13)由被控对象线性化模型(式(12)，可以利用线性控制理论设计线性控制器(式(13)。rS但仅通过一个静态工作点设计控制器会丢失大量 Go的非线性静态和动态信息，使得控制器在设定值和扰动变化情况下的性能下降27。因此，需要将设定值和干扰信息考虑到控制器模型中。定义设定值的变化范围
25、，在设定值范围内选择 个离散点，可以得到一组闭环工作点：Zr ri,i=1,2,S(14)ridd dmin,dmaxNriSi=1NriSi=1Nri对于每个工作点，定义扰动 的变化范围，并选择个在变化范围内的扰动点，用线性子模型来描述系统在操作点附近的系统行为。因此，可以用线性系统来表示该非线性系统 Go。线性系统即可使用线性控制理论设计线性控制器。上述离散化工作点的方法可以看作一个工艺操作区域的识别问题运用于煤气化炉对象的多煤质模型建立。2.2动态矩阵控制(DMC)动态矩阵控制算法是一种基于被控对象阶跃响应模型的预测算法，具有易于建模、鲁棒性强的特点28。动态矩阵控制算法包括预测模型、滚
26、动优化、反馈校正 3 个原则。2.2.1基于对象阶跃响应的预测模型通过对象的传递函数或状态空间关系可以得到系统的阶跃响应系数矩阵，也称动态矩阵。Su=Su100Su2Su100.SuMSuMSu1.SuPSuP1SuPM+1(15)Su RPMSu1SuPM+1PMP MPM其中，,为对象的阶跃响应系数；和分别为预测步长和控制步长，需定义为。预测步长 和控制器的稳定性和快速性有关,控制步长影响系统的灵敏度。k+1P基于阶跃响应系数矩阵得到被控变量在时刻的 步预测值：Y(k+1)=Y(k)+SuU(k)(16)Y(k)RP1U(k)RM1式中，为被控变量的预测值；为操纵变量的变化量。2.2.2滚
27、动优化策略控制器的理想控制目标是通过尽可能小的调整，实现被控变量实际值跟踪设定值的效果。由此定义优化求解的目标函数为minU(k)Y(k)R(k)2y+U(k)2u(17)R(k)RP1y RPPu RMM其 中，为 被 控 变 量 设 定 值；和分别为跟踪误差作用的加权对角矩阵和控制作用的加权对角矩阵，用于表征目标函数的权值。2.2.3反馈校正方式阶跃响应模型不能完全准确地实现控制系统的1750煤炭学报2023年第48卷Y(k+1)Y0(k+1)E(k+1)RP1动态描述，因此需要通过反馈校正，实时对输出反馈进行处理，改善预测值，提高控制器的鲁棒性。将系统的预测值与实际输出比较，得到系统的反
28、馈校正偏差矩阵为E(k+1)=Y0(k+1)Y(k+1)(18)将反馈校正部分代入预测值，校正预测方程得到带有反馈校正的模型预测值为Y(k+1)=Y(k)+SuU(k)+HE(k+1)(19)H RPP式中，为反馈校正对角矩阵。2.2.4在线求解U(k)U(k)Y(k)控制作用的优化求解过程中，考虑以下 3 种变量约束：操纵变量；操纵变量的变化量；被控变量。将约束引入控制作用求解的线性二次规划问题中，实时在线求解 DMC 的控制作用。引理引理若被控对象为渐近稳定，且多个静态动作点能够表征整个模型，则存在一个控制器是渐近收敛的29。即可通过多模型切换策略筛选控制器，使当前控制器保持收敛状态。2.
29、3多煤质控制模型切换准则对于煤气化工作过程，可以由一组局部的线性模型来表述。煤质干扰对控制器的影响可以看作一个工况变化的过程。控制器模型切换的目标是从气化炉模型集中挑选出尽可能满足当前工况的控制器模型。煤质模型数据集的构建是在不影响正常工艺操作的基础上，离线测量扰动变化，采集离线数据构造的多工况控制器模型集。在线控制时，通过控制器性能指标监测及模型切换准则选取合适的控制器模型。本文提出的基于多煤质控制模型的动态矩阵控制器结构如图 3 所示。控制器模型煤气化炉对象煤质模型集ISE-TSV指标监测ruy筛选最优控制模型图3多煤质控制模型的动态矩阵控制器结构Fig.3Dynamicmatrixcon
30、trollerstructureofmultiplecoalqualitycontrolmodels采用累积平方误差(ISE)和总平方波动(TSV)作为综合性能指标评价先进控制器性能30。累积平方误差用于表征控制器的抗干扰性能，它的定义为被控变量在一段时间的控制偏差，即ISE:=J1=nCVi=11iJ1i(20)J1i=wT00ei(t)2dt=wT00yi(t)ri(t)2dt TT0k=0yi(k)ri(k)2i=1,2,nCV(21)nCVT0Tei(t)iei(t)=yi(t)ri(t)i1i式中，为被控变量个数；为控制器评价时长；为采样周期；为第 个被控变量的控制偏差，定义为，即第
31、 个被控变量实际值与设定值的偏差；为对应被控变量的 ISE 权值。总平方波动用于表征控制器的控制消耗，它的定义为操纵变量在一段时间的输出波动，即TSV:=J2=nMVi=12iJ2i(22)J2i=wT00(dui(t)dt)2dt 1TT0k=0ui(k+1)ui(k)2i=1,2,nMV(23)nMVui(t)i2i式中，为操纵变量个数；为第 个操纵变量输出；为对应操纵变量的 TSV 权值。将上述 2 个性能指标加权后得到模型切换的指标判别为J=J1+J2(24)1i、2iJ根据控制需求及 2 个性能指标之间的重要程度进行加权给定，越小代表控制器性能越优。T0JmaxJmax在煤气化炉性能
32、监测过程中，需定义 ISE-TSV 的评价时长，以及允许接受的性能指标上限即可接受的最大模型失配程度。当 APC 投用过程中，监测到的性能指标大于时，启动模型切换算法。ncj考虑到煤气化过程的模型失配主要由煤质改变造成，对于模型预测控制来说，模型失配将产生较大的预测偏差。因此本文采取最小预测偏差作为模型切换准则。预测偏差定义为一定时间内被控变量当前值与预测值之间的累积偏差，也就是预测模型与实际模型的匹配程度。需要对 个控制器模型分别求取预测偏差，其中对第 个控制器模型求解偏差过程如下：E(j)=nCVi=11iepre,ii=1,2,nCVj=1,2,nc(25)epre,i=wT00?yi(
33、t)ypre,i(t)?dt TT0k=0?yi(k)ypre,i(k)?i=1,2,nCV(26)第4期张天辰等：基于多模型切换控制的煤气化工业过程先进控制1751epre,iiT0ypre,i(t)iyi(t)其中，为第 个被控变量的预测偏差，定义为时间内的实际值与预测值的累积偏差；为第 个被控变量的预测值；为实际值。式(25)中的权值同式(22)中 ISE 权值。Jmaxnc当监测到性能指标高于性能指标上限时，启动控制器切换算法。对 个控制器模型分别求取预测偏差，在煤质模型数据集中选取最小预测偏差的控制器，即Emin=minE(j),j=1,2,nc(27)从而得到对应的控制增量U(t)
34、=UEmin(t)(28)Ypm RnCV1控制器切换过程中，为了防止被控变量超限需要提前定义约束条件31。定义切换过程中引起的被控变量的最大允许预测控制误差，当该误差过大时，说明当前工况模型切换到另一工况模型引起的装置偏差过大，该切换将被停止，此时需要等工况稳定后再次进行切换。限定了上述约束，可实现控制器模型切换过程的平顺，使得控制器投用及切换过程的波动在预设范围内。同时根据引理与性能指标权值的限定，控制器在正常工况运行时不会因 ISE-TSV性能指标的异常而导致频繁切换。基于多模型切换的动态矩阵控制算法运行框如图 4 所示。根据动态矩阵控制及模型切换理论，北京化工大学过程监测与先进控制实验
35、室开发了多模型动态矩阵控制软件 Wisdom-Controller(登记号：2017SR610595)。3应用研究3.1UniSim 模型根据煤气化原理及德士古气化工艺，采用 Honey-well 的 UniSimDesign 软件，搭建了水煤浆气化模型，如图 5 所示。流程模拟了水煤浆制作、气化以及灰水处理等关键过程。气化炉内部基于物料平衡原理和能量平衡原理，进行相关的动力学和热力学反应，得到粗合成气等产品。UniSim 平台提供了调整气化过程水煤浆煤质成分的模块，如图 6 所示，可以得到不同质量指标的煤质。由于煤质特性需要离线分析得到，因此建模过程为离线进行。使用 OPC 通信技术，让 U
36、niSim 模型点位与 Wis-dom-Controller 程序点位相连接。OPC 通信技术是工选用当前工况控制器模型完成控制器切换?返回NYY计算 ISE-TSV 性能指标 J以当前控制器继续运行N以当前控制器继续运行定义控制器评价时长 T0,定义 ISE与 TSV 指标权重,定义允许接受的性能指标上限 Jmax,定义模型切换的最大允许预测误差 YpmJ 是否大于 Jmax?启动控制器性能切换算法,计算数据库中各模型预测偏差选择最小预测偏差对应的控制器增量 UEmin判断控制器切换至 UEmin后对应的预测误差是否大于Ypm图4多模型切换的动态矩阵控制算法框图Fig.4Blockdiagr
37、amofdynamicmatrixcontrolalgorithmformultiplemodelsswitching图5UniSim 水煤浆气化模型Fig.5UniSimcoalwaterslurrygasificationmodel1752煤炭学报2023年第48卷业领域常用的数据传输技术，相当于一个工业标准接口，能够使不同通信协议的下层现场设备中的数据得以与上位机的应用程序进行交互，主要通过 OPC-Cli-ent 与 OPC-Server 进行数据交互。采用 Matrikon 软件作为仿真的 OPC 服务器，UniSim 与 Wisdom-Con-troller 之间的数据传输过程如图
38、 7 所示。MatrikonOPC-ServerWisdom-ControllerOPC-ClientUnisimOPC-Client读写读写图7OPC 数据传输过程Fig.7OPCDatatransmissionprocess3.2模型切换控制器的设计与实施3.2.1煤气化炉模型的建立选取气化炉部分工况下的运行数据，通过聚类分析的方式32，建立 5 种煤质指标的煤气化炉运行模型，见表2。选取目的是以5 种指标表征数据范围内的工况。以煤质 1 为例，根据 1.4 节中的控制需求，通过阶跃测试建立模型。即在被控变量处于稳态的前提下，保持其他操作变量不变，人为使一个 MV 发生阶跃改变(阶跃比例需
39、要根据工艺来决定，一般为 3%5%)，从而引起 CV 发生动态变化，直至其达到新的稳态的过程。表 2 气化炉模型煤质类型Table 2 Gasifier model coal quality type煤质编号煤质指标/%固定碳挥发分灰分15635.88.225833.58.536031.58.546330.07.055535.010.0即在煤气化炉系统稳定时，在氧气回路 FIC-100的操作面板中，提高氧气流量的 SP(设定值)大约 3%，如图 8 所示，待被控变量稳定后截取该时间段的数据进行辨识处理。煤浆回路 FIC-101 进行相近比例的阶跃操作。图8氧气回路操作面板Fig.8Oxygen
40、circuitoperationpanel煤质 1 工况下氧气流量(MV1)的阶跃测试结果如图 9 所示，煤浆流量(MV2)的阶跃测试结果如图 10所示。工业上先进控制实施周期一般为 1min；CV1表示氧煤比；CV2 表示煤气化炉炉膛温度。图6UniSim 煤质质量指标调整Fig.6UniSimcoalqualityindexadjustment氧气流量(MV1)/(103 Nm3h1)5453525150520510500490480491 5001 4501 4001 3501 3001 25048MV1CV1CV2采样周期/min0100200300氧煤比(CV1)/(Nm3m3)炉膛
41、温度(CV2)/C图9MV1 氧气流量阶跃测试Fig.9MV1oxygenflowsteptest第4期张天辰等：基于多模型切换控制的煤气化工业过程先进控制1753通过 Wisdom-Controller 的离线辨识模块，将阶跃跳变数据辨识为时域下的传递函数模型，得到在当前煤质工况下，该静态工作点线性化后的煤气化炉传递函数模型，离线重复上述建模过程，最终得到基于不G12(s)同煤质模型的煤气化炉模型集，见表 3，并嵌入 Wis-dom-Controller 软件的煤质模型数据集。其中，为 MV1 到 CV2 的传递函数；s 为传递函数频域表达式的自变量。表 3 气化炉煤质 15 工况下的传递函
42、数模型Table 3 Transfer function model of gasifier under coal quality 15 condition煤质MV1MV2CV1CV2CV1CV21G11(s)=0.009 20.137s+1G12(s)=0.048 50.005s2+0.14s+1G21(s)=0.6210.348s+1e0.97sG22(s)=1.407s19.90.005s2+0.234s+12G11(s)=0.009 40.086s+1G12(s)=0.0510.14s+1G21(s)=4.24s3.750.037s2+1.18s+1G22(s)=19.50.007s2
43、+0.185s+13G11(s)=0.002s+0.0090.003s2+0.816s+1G12(s)=0.0490.735s+1G21(s)=3.570.444s+1G22(s)=11.97s19.30.329s2+1.25s+14G11(s)=0.002s+0.0090.002s2+0.839s+1G12(s)=0.011s+0.0520.135s2+0.991s+1G21(s)=3.240.345s+1G22(s)=3.24s19.50.075s2+0.779s+15G11(s)=0.002s+0.0090.003s2+0.845s+1G12(s)=0.008s+0.050.101s2+
44、0.95s+1G21(s)=0.75s3.420.003s2+0.669s+1G22(s)=2.97s18.70.067s2+0.764s+13.2.2煤气化炉模型切换控制的实施气化炉运行过程中煤质未变化，即模型无失配情况的控制效果如图 11 所示。其中 PV 代表实际值，SP 代表设定值。由图 11 可以看出，在对应控制器控制效果下，调整操纵变量氧气和煤浆流量，使被控变量氧煤比和气化炉温度达到追踪设定值的效果，并且控制器可以达到稳定控制的目的。在气化炉运行过程中，在线调整 UniSim 煤质属性，模拟煤质变化过程。煤质变化产生模型失配，导致控制器性能下降，甚至发散，未采用模型切换算法即单一模
45、型的控制效果如图 12 所示。由图 12 可以看出，模型失配导致以下几个现象：阀门长周期波动、MV2CV1CV2采样周期/min0100200300煤浆流量(MV2)/(th1)140130120110520540500480460440氧煤比(CV1)/(Nm3m3)1 6001 5001 4001 3001 2001 100炉膛温度(CV2)/C图10MV2 煤浆流量阶跃测试Fig.10MV2coalslurryflowsteptestMV1MV2CV1-PVCV2-PVCV1-SPCV2-SP50.0550.0049.9549.9049.8549.80煤浆流量(MV2)/(th1)130
46、128126124氧气流量(MV1)/(103 Nm3h1)0100200300采样周期/min4904954854804754701 3601 3801 3401 3201 3001 2801 260氧煤比(CV1)/(Nm3m3)炉膛温度(CV2)/C图11无煤质变化控制效果Fig.11Controleffectofnocoalqualitychange1754煤炭学报2023年第48卷控制器连续朝一个方向变化无法拉回、被控变量波动过大。这些情况在实际工艺过程中都有可能导致煤气化全工段波动，甚至会触发联锁停车问题。J11=0.009、12=0.003、21=0.001、22=10Jmax=
47、50Ypm=530 T采取本文提出的多模型切换方法应对煤质变化过程，定义控制器评价时长 T0为 15 个采样周期 T；ISE-TSV 性能指标 中权值定义如下：；性能指标上限；模型切换最大允许误差。上述权值及指标根据仿真测试，结合控制要求，试凑得到。采取多模型动态矩阵切换策略后的气化炉控制效果如图 13 所示。图 14 展示了单一模型控制与模型切换控制的性能指标对比，可以看出，采取控制器切换算法应对煤质变化的工况时，性能指标首先会监测上升，并且在一段时间内进行监控。当达到性能指标监控周期后，启动模型切换算法，通过筛选最小预测偏差的控制器模型，得到最适用的模型输出，减少系统的模型失配程度，使控制
48、器性能指标下降，代表控制器性能得到改善。煤浆流量(MV2)/(th1)1601401001206080706050403020100氧气流量(MV1)/(103 Nm3h1)100200300400500采样周期/minMV1MV2CV1-PVCV2-PV540520500480460440420氧煤比(CV1)/(Nm3m3)1 6001 7001 5001 4001 3001 200炉膛温度(CV2)/C图12单一模型算法应对煤质变化控制效果Fig.12Controleffectofsinglemodelalgorithmoncoalqualitychange煤浆流量(MV2)/(th1)
49、1401301001101209050.5050.0049.5049.00氧气流量(MV1)/(103 Nm3h1)MV1MV2CV1-PVCV2-PV520500480460440氧煤比(CV1)/(Nm3m3)0100200300400500采样周期/min1 6001 5001 4001 3001 200炉膛温度(CV2)/C图13模型切换算法应对煤质变化控制效果Fig.13Controleffectofmodelswitchingalgorithmoncoalqualitychange模型切换算法控制器性能指标单一模型算法5004003002001000100200300400采样周期
50、/min图14单一模型算法与模型切换算法控制器性能指标对比Fig.14Comparisonofcontrollerperformanceindexbetweensinglemodelalgorithmandmodelswitchingalgorithm3.3工业测试在某化工公司的煤气化炉装置上进行了工业应用标定测试，系统架构如图 15 所示，操作画面如图 16 所示。测试过程人为改变了煤样的进料情况，控制效果如图 17 所示。APC 投用前后氧煤比及气化炉温度的方差对比、针对装置经济性分析，APC 投用前后比煤耗对比及相同负荷下的合成气流量及有效气含量(干合成气中 CO 和 H2所占的体积比)