一种基于多变量约束区间预测控制的氨法脱硫优化控制方法技术

本发明专利技术公开了一种基于多变量约束区间预测控制的氨法脱硫优化控制方法，本发明专利技术通过采用多变量约束区间预测控制方法，能够更好地协调脱硫塔吸收段与氧化段的进氨水流量，在保证脱硫效率的前提下，改善系统控制品质，降低系统运行成本，促进系统经济环保地运行。

【技术实现步骤摘要】
一种基于多变量约束区间预测控制的氨法脱硫优化控制方法
本专利技术涉及热工自动控制领域，特别是涉及一种基于多变量约束区间预测控制的氨法脱硫优化控制方法。
技术介绍
控制SO2等污染气体的排放是我国当前促进火力发电行业的技术进步和可持续发展、解决环境污染问题的重大任务。2012年1月1日起实施的《火电厂大气污染物排放标准》规定：自2014年7月1日起，现有燃煤锅炉SO2排放的最高限值为200mg/m3。因此，减少火力发电行业SO2的排放不仅对减少环境污染具有重要现实意义，同时也是国家对于建立资源节约型和环境友好型企业的迫切需求。氨法脱硫工艺是采用氨作为吸收剂除去烟气中的SO2和部分NOx的工艺。相对于钙法脱硫工艺来说，反应速度快、反应完全、吸收剂利用率高，可以达到很高的脱硫效率，而且系统简单、设备体积小、能耗低。近年来国内围绕氨法烟气脱硫这一新的环保技术已经开展了大量实验和应用研究工作，然而其中绝大多数工作都集中在对脱硫工艺的改进方面，但是脱硫装置能否充分发挥其效能，达到设计的各项技术指标，还与其运行和操作的水平直接相关。尽管目前氨法脱硫装置基本都安装了分布式计算机控制系统(DCS)，但是主要用于实现设备启停和联锁保护功能，很少能投入闭环控制。造成这些情况的原因在于整个脱硫流程非常复杂，涉及到反应、吸收、传质、传热、干燥、分离及结晶等过程，循环溶液pH等被控对象具有很强的非线性和大迟延特性，操作变量之间相互关联、耦合严重，烟气流量、烟气SO2浓度、氨水浓度等干扰众多。因此以PID控制器为核心的常规控制方案很难取得满意的控制效果。此外，重要过程参数，如SO2浓度、溶液PH等不能准确测量且故障率高，也是一个重要原因。例如中国专利，申请号为：201610006777.8,公开了一种水泥窑烟气氨法脱硫系统及方法。本专利技术的水泥窑烟气氨法脱硫系统及方法，进行脱硫工艺能够实现80％以上的脱硫效率，选择合适的氨硫比不会造成额外的氨逃逸污染，而且无其他污染。运行成本主要为氨水的费用，降低了传统脱硫工艺的一次投资成本及运行成本。但是本申请中氨水的喷射量为定值，并没有根据水泥窑中SO2浓度对氨水的喷射量进行实施控制，导致氨水浪费。例如中国专利，申请号为：201210203236.6，公开了一种基于多参数的氨法烟气脱硫效预测方法，选用四种不同的人工智能计算模型.将氨法脱硫系统运行过程中所采集的多组烟气量、循环泵流量、浓缩泵流量、氨浓度、吸收液浓度、液气比、进口烟气温度、耗氨量、喷淋浆液密度、喷淋塔浆液pH值和预洗涤塔浆液pH值等参数作为四个模型的输入变量。对各模型分别进行训练，建立四个脱硫参数与脱硫效率间的非线性函数关系。再将实时监侧到的参数分别传输至己训练好的人工智能模型中，对脱硫效率做出预测。虽然本申请中采用了计算模型，将氨法脱硫系统运行过程中所采集的多组烟气量、循环泵流量、浓缩泵流量、氨浓度、吸收液浓度、液气比、进口烟气温度、耗氨量、喷淋浆液密度、喷淋塔浆液pH值和预洗涤塔浆液pH值等参数作为四个模型的输入变量，进行了实施控制，但并没有进行闭环控制，因此预测控制精度低。
技术实现思路
本专利技术的目的在于提供一种提高系统经济性的基于多变量约束区间预测控制策略的氨法脱硫优化控制方法，解决现有技术中尽管目前氨法脱硫装置基本都安装了分布式计算机控制系统(DCS)，但是主要用于实现设备启停和联锁保护功能，很少能投入闭环控制的技术问题。为实现上述专利技术目的，本专利技术采用以下技术方案：以脱硫后SO2浓度和循环浆液pH为被控量，吸收段、氧化段的进氨水流量为控制量，设计基于多变量约束区间预测控制的优化控制器，利用预测模型提前预测未来时刻被控量的变化，从而计算求解当前时刻的最优控制量。通过设置输出软约束，引入松弛变量，实现脱硫后SO2浓度的定值控制和循环浆液pH的区间控制，改善系统控制品质，提高系统经济性，所述多变量预测控制方法具体步骤如下：1)获取氨法脱硫对象的辨识模型，在稳态工况下，分别以脱硫塔吸收段和氧化段的进氨水流量为输入，对出口SO2浓度和循环浆液pH进行开环阶跃响应试验，经滤波平滑后，分别得到吸收段、氧化段进氨水流量对出口SO2浓度、循环浆液pH的阶跃响应模型系数为a11(1),…,a11(N)、a12(1),…,a12(N)、a21(1),…,a21(N)和a22(1),…,a22(N)，其中，N为建模时域，对不同输入、输出间的阶跃响应采用统一的建模时域；2)设置多变量预测控制器的相关参数，包括采样时间Ts、预测步数P、控制步数M，输出误差权矩阵Q，控制权矩阵R；Ts可以利用经验规则T95/Ts＝5～15来选取，其中，T95为过渡过程上升到95％的调节时间；P一般选为等于对象阶跃响应的上升时间；M一般选1或2；Q＝I,R＝I，I为单位矩阵；3)控制器参数确定后，采用公式(1)所述的预测模型对脱硫系统未来输出进行预测：其中，m为输入变量个数，p为输出变量个数；表示在k时刻对未来P时刻的被控量i的预测值，表示在k时刻对未来P时刻的被控量i的预测初值；Δuj,M(k)表示在k时刻对未来M个采样时刻控制量j的增量的估计值；即由单位阶跃响应系数组成的动态矩阵；为表达方便，用yt来表示即表示利用预测模型对未来P个采样时刻的所有被控量的预测值；4)控制器状态初始化，即在k时刻，检测当前时刻脱硫后SO2浓度和循环浆液pH的测量值y(k)，并令y0＝yt＝y(k),t＝1,…,P，其中，y0为利用预测模型对当前时刻被控量的预测值；5)计算输出预测偏差e＝y(k)-y0；6)进行反馈校正，yt+hte→yt,t＝1,…,P，其中，ht为校正系数，取为1；7)计算吸收段进氨水流量增量Δu1(k)和氧化段进氨水流量增量Δu2(k)，取性能指标函数为式：其中，εi为第i个输出超出上下限约束的大小，uj(k)为第j个控制量，yi(k)为第i个被控量，Δuj(k)为第j个控制量的变化量，uj,min为第j个控制量的变化下限，uj,max为第j个控制量的变化上限，Δuj,min为第j个控制量的变化率下限，Δuj,max为第j个控制量的变化率上限，yi,min为第i个被控量的输出区间下限，yi,max为第i个被控量的输出区间上限，表示被控量允许超过上下限的范围；表示第i个被控量允许超过上下限的范围；8)计算吸收段进氨水流量u1(k)＝u1(k-1)+Δu1(k)，以及氧化段进氨水流量u2(k)＝u2(k-1)+Δu2(k)；9)输出u1(k)和u2(k)，计算并更新控制后的输出预测初值：而后在每个采样周期内，重复执行第5)步到第9)步。进一步改进，所述步骤2)中，Ts可以利用经验规则T95/Ts＝5～15来选取，其中，T95为过渡过程上升到95％的调节时间；P一般选为等于对象阶跃响应的上升时间；M一般选1或2；Q＝I，R＝I，I代表单位矩阵。进一步改进，取采样时间Ts＝250秒，预测步数P＝30，控制步数M＝1，输出误差权矩阵Q＝I60，控制权矩阵R＝2×I2。进一步改进，所述吸收段、氧化段进氨水流量的约束范围均为[0，8]m3/h,变化速率均为[-0.5，0.5]m3/h；取脱硫后SO2浓度的初始设定值为35mg/Nm3，循环浆液pH的控制区间为[y2,min本文档来自技高网...

【技术保护点】
1.一种基于多变量约束区间预测控制的氨法脱硫优化控制方法，其特征在于：采用多变量约束区间预测控制方法对氨法脱硫系统进行控制，包括控制脱硫塔吸收段和氧化段的进氨水流量，并设置输出软约束，引入松弛变量，最终实现脱硫后SO2的定值控制和循环浆液pH的区间控制，改善系统的控制品质，所述多变量约束区间预测控制方法具体步骤如下：1)获取氨法脱硫对象的阶跃响应模型，在稳态工况下，分别以脱硫塔吸收段和氧化段的进氨水流量为输入，对出口SO2浓度和循环浆液pH值进行开环阶跃响应试验，经滤波平滑后，分别得到吸收段、氧化段进氨水流量对出口SO2浓度、循环浆液pH的阶跃响应模型系数为a11(1),…,a11(N)、a12(1),…,a12(N)、a21(1),…,a21(N)和a22(1),…,a22(N)，其中，N为建模时域，对不同输入输出间的阶跃响应采用统一的建模时域；2)设置多变量预测控制器的相关参数，包括采样时间Ts、预测步数P、控制步数M，输出误差权矩阵Q，控制权矩阵R；3)控制器参数确定后，采用公式(1)所述的预测模型对脱硫系统未来输出进行预测：

【技术特征摘要】
1.一种基于多变量约束区间预测控制的氨法脱硫优化控制方法，其特征在于：采用多变量约束区间预测控制方法对氨法脱硫系统进行控制，包括控制脱硫塔吸收段和氧化段的进氨水流量，并设置输出软约束，引入松弛变量，最终实现脱硫后SO2的定值控制和循环浆液pH的区间控制，改善系统的控制品质，所述多变量约束区间预测控制方法具体步骤如下：1)获取氨法脱硫对象的阶跃响应模型，在稳态工况下，分别以脱硫塔吸收段和氧化段的进氨水流量为输入，对出口SO2浓度和循环浆液pH值进行开环阶跃响应试验，经滤波平滑后，分别得到吸收段、氧化段进氨水流量对出口SO2浓度、循环浆液pH的阶跃响应模型系数为a11(1),…,a11(N)、a12(1),…,a12(N)、a21(1),…,a21(N)和a22(1),…,a22(N)，其中，N为建模时域，对不同输入输出间的阶跃响应采用统一的建模时域；2)设置多变量预测控制器的相关参数，包括采样时间Ts、预测步数P、控制步数M，输出误差权矩阵Q，控制权矩阵R；3)控制器参数确定后，采用公式(1)所述的预测模型对脱硫系统未来输出进行预测：其中，m、p分别为输入输出变量个数；表示在k时刻对未来P时刻的被控量i的预测值，表示在k时刻对未来P时刻的被控量i的预测初值；Δuj,M(k)表示在k时刻对未来M个采样时刻控制量j的增量的估计值；即由单位阶跃响应系数组成的动态矩阵；为表达方便，用yt来表示即表示利用预测模型对未来P个采样时刻的所有被控量的预测值；4)控制器状态初始化，即在k时刻，检测当前时刻脱硫后SO2浓度和循环浆液pH的测量值y(k)，并令y0＝yt＝y(k),t＝1,…,P，其中，y0为利用预测模型对当前时刻被控量的预测值；5)计算输出预测偏差e＝y(k)-y0；6)进行反馈校正，yt+hte→yt,t＝1,…,P，其中，ht为校正系数，取为1；7)计算吸收段进氨水流量增量Δu1(k)和氧化段进氨水流量增量Δu2(k)，取性能指标函数为式：其中，εi为第i个输出超出上下限约束的大小，uj(k)为第j个控制量，yi(k)为第i个被控量，Δuj(k)为第j个控制量的变化量，uj,min为第j个控制量的变化下限，uj,max...

【专利技术属性】
技术研发人员：赵东晓，宋祎，丛经纬，沙志强，宋俊，金晓，刘洪兵，孙玉芬，陈勇，
申请(专利权)人：中国石化扬子石油化工有限公司，中国石油化工股份有限公司，
类型：发明
国别省市：江苏,32