专利摘要

本发明公开了一种预测函数控制优化的废塑料炼油裂解炉炉膛压力控制方法。本发明方法首先基于裂解炉炉膛压力对象的阶跃响应数据建立炉膛压力对象的模型，挖掘出基本的对象特性；然后依据预测函数控制的特性去整定相应PI控制器的参数；最后对裂解炉炉膛压力对象实施PI控制。本发明将预测函数控制的性能赋给了PI控制，有效地提高了传统控制方法的性能，同时也促进了先进控制方法的应用。

权利要求

1.预测函数控制优化的废塑料炼油裂解炉炉膛压力控制方法，其特征在于该方法的具体步骤是：

步骤(1).通过裂解炉炉膛压力对象的实时阶跃响应数据建立被控对象的模型，具体方法是：

Ⅰ.将过程的比例积分控制器停留在手动操作状态，操作拨盘使其输出有个阶跃变化，由记录仪表记录实际过程的输出值，将实际过程输出值y_p(k)的响应曲线转换成无量纲形式y_p^*(k)，具体是：

y_p^*(k)＝y_p(k)/y_p(∞)

其中，y_p(∞)是比例积分微分控制器的输出有阶跃变化时的实际过程输出y_p(k)的稳态值；

Ⅱ.选取满足y_p^*(k₁)＝0.39以及y_p^*(k₂)＝0.63的两个计算点k₁和k₂，依据下式计算过程对象的模型参数K_m、T和τ：

K_m＝y_p(∞)/q

T＝2(k₁-k₂)

τ＝2k₁-k₂

最后得到的过程对象的传递函数为：

$G\left(s\right)=\frac{K_m}{\mathrm{Ts}+1}{e}^{-\mathrm{\tau s}}$

其中，q为过程的比例积分微分控制器输出的阶跃变化幅度，G(s)为过程对象的传递函数，s为拉普拉斯变换算子，K_m为模型的增益系数，T为模型的时间常数，τ为模型的滞后时间参数；

步骤(2).设计过程对象的PI控制器，具体方法是：

a.对得到的传递函数在采样时间T_s下加一个零阶保持器离散化，得到离散模型为

y_m(k)＝a_my_m(k-1)+K_m(1-a_m)u(k-1-L)

y_m(k)为k时刻的过程对象模型预测输出， u(k-1-L)为k-1-L时刻的过程对象的控制输入，L为离散传递函数模型的时滞，L＝τ/T_s；

b.计算过程对象去掉纯滞后以后在预测函数控制下的第P步预测输出，形式如下：

y_mav(k)＝a_my_mav(k-1)+K_m(1-a_m)u(k-1)

y_mav(k+P)＝a_m^Py_mav(k)+K_m(1-a_m^P)u(k)

其中，P为预测步长，y_mav(k+P)为k时刻去掉纯滞后的过程对象在预测函数控制下的第P步预测输出，y_mav(k)为k时刻去掉纯滞后的过程模型输出；

c.修正当前时刻的实际输出得到包含未来预测信息的新的过程际输出值，形式如下：

y_Pav(k)＝y_P(k)+y_mav(k)-y_mav(k-L)

其中，y_Pav(k)为校正得到的k时刻包含未来预测信息的新的过程输出值，y_P(k)为k时刻的实际输出值；

d.选取预测函数控制方法的参考轨迹y_r(k+P)以及目标函数J，形式如下：

y_r(k+P)＝β^Py_p(k)+(1-β^P)c(k)

J＝min(y_r(k+P)-y_mav(k+P)-e(k))²

e(k)＝y_pav(k)-y_mav(k)

其中，β为参考轨迹柔化系数，c(k)为k时刻的设定值，e(k)为k时刻校正的误差值；

e.依据步骤d中的目标函数求解PI控制器中的参数，这里将控制量u(k)进行变换：

u(k)＝u(k-1)+K_p(e₁(k)-e₁(k-1))+K_ie₁(k)

e₁(k)＝βy_p(k-1)+(1-β)c(k-1)-y_p(k)

进一步化简为：

u(k)＝u(k-1)+w(k)^ΤE(k)

w(k)＝[w₁(k),w₂(k)]^Τ

w₁(k)＝K_p+K_i,w₂(k)＝-K_p

E(k)＝[e₁(k),e₁(k-1)]^Τ

其中，K_p、K_i分别为PI控制器的比例、积分参数，e₁(k)为k时刻参考轨迹值与实际输出值之间的误差，Τ为矩阵的转置符号；

结合上述式子，可以求得：

$w\left(k\right)=\frac{(y_r(k+P)-{a_m}^P{y}_{\mathrm{mav}}\left(k\right)-K_m(1-{a_m}^P)u(k-1))E}{K_m(1-{a_m}^P)E^{T}E}$

进一步可以得到：

K_p＝-w₂(k)

K_i＝w₁(k)-K_P

f.得到PI控制器的参数K_p、K_i以后构成控制量u(k)作用于被控对象，u(k)＝u(k-1)+K_p(e₁(k)-e₁(k-1))+K_ie₁(k)；

g.在下一时刻，依照b到f中的步骤继续求解PI控制器新的参数K_p、K_i，依次循环。

说明书

技术领域

本发明属于自动化技术领域，涉及一种基于预测函数控制优化的废塑料炼油裂解炉炉膛温度比例积分(PI)控制方法。

背景技术

在实际工业对象的控制中，由于硬件、成本、实施难度等方面的限制，一些先进控制方法虽然得到了一定程度的应用，但是目前占主流的仍是PID控制。炼油裂解炉是石油化工生产过程中的重要装置，其中炉膛压力对裂解过程有着很重要的影响，炉膛压力过大可能会带来危险，炉膛压力过低会导致裂解效率变低。由于先进控制方法应用受限，目前裂解炉炉膛压力的控制通常采用比例积分(PI)控制。预测函数控制作为先进控制方法的一种，在裂解炉炉膛压力的控制中相比PI控制拥有更好的控制性能，如果能将预测函数控制的性能赋给PI控制，那将进一步推进先进控制方法的应用，同时也能获得更好的实际控制效果。

发明内容

本发明的目的是针对现有先进控制方法的应用不足之处，提供一种基于预测函数控制优化的裂解炉炉膛压力PI控制方法，以获得更好的实际控制性能。该方法通过结合预测函数控制和PI控制，得到了一种带有预测函数控制性能的PI控制方法。该方法在继承预测函数控制优良性能的同时也保证形式简单并满足实际工业过程的需要。

本发明方法首先基于裂解炉炉膛压力对象的阶跃响应数据建立炉膛压力对象的模型，挖掘出基本的对象特性；然后依据预测函数控制的特性去整定相应PI控制器的参数；最后对裂解炉炉膛压力对象实施PI控制。

本发明的技术方案是通过数据采集、模型建立、预测机理、优化等手段，确立了一种基于预测函数控制优化的PI控制方法，利用该方法可有效提高控制的精度与稳定性。

本发明方法的步骤包括：

步骤(1).通过裂解炉炉膛压力对象的实时阶跃响应数据建立被控对象的模型，具体方法是：

a.将过程的比例积分控制器停留在手动操作状态，操作拨盘使其输出有个阶跃变化，由记录仪表记录实际过程的输出值，将实际过程输出值y_p(k)的响应曲线转换成无量纲形式y_p^*(k)，具体是：

y_p^*(k)＝y_p(k)/y_p(∞)

其中，y_p(∞)是比例积分微分控制器的输出有阶跃变化时的实际过程输出y_p(k)的稳态值。

b.选取满足y_p^*(k₁)＝0.39以及y_p^*(k₂)＝0.63的两个计算点k₁和k₂，依据下式计算过程对象的模型参数K_m、T和τ：

K_m＝y_p(∞)/q

T＝2(k₁-k₂)

τ＝2k₁-k₂

最后得到的过程对象的传递函数为：

G(s)=KmTs+1e-τs]]>

其中，q为过程的比例积分微分控制器输出的阶跃变化幅度，G(s)为过程对象的传递函数，s为拉普拉斯变换算子，K_m为模型的增益系数，T为模型的时间常数，τ为模型的滞后时间参数。

步骤(2).设计过程对象的PI控制器，具体方法是：

a.对得到的传递函数在采样时间T_s下加一个零阶保持器离散化，得到离散模型为

y_m(k)＝a_my_m(k-1)+K_m(1-a_m)u(k-1-L)

y_m(k)为k时刻的过程对象模型预测输出， u(k-1-L)为k-1-L时刻的过程对象的控制输入，L为离散传递函数模型的时滞，L＝τ/T_s。

b.计算过程对象去掉纯滞后以后在预测函数控制下的第P步预测输出，形式如下：

y_mav(k)＝a_my_mav(k-1)+K_m(1-a_m)u(k-1)

y_mav(k+P)＝a_m^Py_mav(k)+K_m(1-a_m^P)u(k)

其中，P为预测步长，y_mav(k+P)为k时刻去掉纯滞后的过程对象在预测函数控制下的第P步预测输出，y_mav(k)为k时刻去掉纯滞后的过程模型输出。

c.修正当前时刻的实际输出得到包含未来预测信息的新的过程际输出值，形式如下：

y_Pav(k)＝y_P(k)+y_mav(k)-y_mav(k-L)

其中，y_Pav(k)为校正得到的k时刻包含未来预测信息的新的过程输出值，y_P(k)为k时刻的实际输出值。

d.选取预测函数控制方法的参考轨迹y_r(k+P)以及目标函数J，形式如下：

y_r(k+P)＝β^Py_p(k)+(1-β^P)c(k)

J＝min(y_r(k+P)-y_mav(k+P)-e(k))²

e(k)＝y_pav(k)-y_mav(k)

其中，β为参考轨迹柔化系数，c(k)为k时刻的设定值，e(k)为k时刻校正的误差值。

e.依据步骤d中的目标函数求解PI控制器中的参数，这里将控制量u(k)进行变换：

u(k)＝u(k-1)+K_p(e₁(k)-e₁(k-1))+K_ie₁(k)

e₁(k)＝βy_p(k-1)+(1-β)c(k-1)-y_p(k)

进一步化简为：

u(k)＝u(k-1)+w(k)^ΤE(k)

w(k)＝[w₁(k),w₂(k)]^Τ

w₁(k)＝K_p+K_i,w₂(k)＝-K_p

E(k)＝[e₁(k),e₁(k-1)]^Τ

其中，K_p、K_i分别为PI控制器的比例、积分参数，e₁(k)为k时刻参考轨迹值与实际输出值之间的误差，Τ为矩阵的转置符号。

结合上述式子，可以求得：

w(k)=(yr(k+P)-amPymav(k)-Km(1-amP)u(k-1))EKm(1-amP)ETE]]>

进一步可以得到：

K_p＝-w₂(k)

K_i＝w₁(k)-K_P

f.得到PI控制器的参数K_p、K_i以后构成控制量u(k)作用于被控对象，u(k)＝u(k-1)+K_p(e₁(k)-e₁(k-1))+K_ie₁(k)。

g.在下一时刻，依照b到f中的步骤继续求解PI控制器新的参数K_p、K_i，依次循环。

本发明提出的一种基于预测函数控制优化的废塑料炼油裂解炉炉膛压力PI控制方法将预测函数控制的性能赋给了PI控制，有效地提高了传统控制方法的性能，同时也促进了先进控制方法的应用。

具体实施方式

以废塑料炼油裂解炉炉膛压力过程控制为例：

裂解炉炉膛压力是裂解炉裂解过程中的重要参数，调节手段采用烟道挡板的开度。

步骤(1).通过裂解炉炉膛压力对象的实时阶跃响应数据建立被控对象的模型，具体方法是：

a.将裂解炉炉膛压力过程的比例积分控制器停留在手动操作状态，操作拨盘使其输出有个阶跃变化，由记录仪表记录炉膛压力过程的实际输出值，将实际输出值y_p(k)的响应曲线转换成无量纲形式y_p^*(k)，具体是：

y_p^*(k)＝y_p(k)/y_p(∞)

其中，y_p(∞)是比例积分控制器的输出有阶跃变化时的炉膛压力过程输出y_p(k)的稳态值。

b.选取满足y_p^*(k₁)＝0.39以及y_p^*(k₂)＝0.63的两个计算点k₁和k₂，依据下式计算炉膛压力过程的模型参数K_m、T和τ：

K_m＝y_p(∞)/q

T＝2(k₁-k₂)

τ＝2k₁-k₂

最后得到的过程模型的传递函数为：

G(s)=KmTs+1e-τs]]>

其中，q为炉膛压力过程的比例积分微分控制器输出的阶跃变化幅度，G(s)为炉膛压力过程的传递函数，s为拉普拉斯变换算子，K_m为炉膛压力过程模型的增益系数，T为炉膛压力过程模型的时间常数，τ为炉膛压力过程模型的滞后时间参数。

步骤(2).设计炉膛压力过程的PI控制器，具体方法是：

a.对得到的炉膛压力过程模型的传递函数在采样时间T_s下加一个零阶保持器离散化，得到离散模型为

y_m(k)＝a_my_m(k-1)+K_m(1-a_m)u(k-1-L)

y_m(k)为k时刻的炉膛压力过程模型预测输出， u(k-1-L)为k-1-L时刻的炉膛压力过程模型的控制输入，L为离散传递函数模型的时滞，L＝τ/T_s。

b.计算炉膛压力过程模型去掉纯滞后以后在预测函数控制下的第P步预测输出，形式如下：

y_mav(k)＝a_my_mav(k-1)+K_m(1-a_m)u(k-1)

y_mav(k+P)＝a_m^Py_mav(k)+K_m(1-a_m^P)u(k)

其中，P为预测步长，y_mav(k+P)为k时刻去掉纯滞后的炉膛压力过程模型在预测函数控制下的第P步预测输出，y_mav(k)为k时刻去掉纯滞后的炉膛压力过程模型输出。

c.修正当前时刻的炉膛压力过程实际输出值得到包含未来预测信息的新的过程际输出值，形式如下：

y_Pav(k)＝y_P(k)+y_mav(k)-y_mav(k-L)

其中，y_Pav(k)为校正得到的k时刻包含未来预测信息的炉膛压力过程的新的过程输出值，y_P(k)为k时刻的炉膛压力过程实际输出值。

d.选取预测函数控制方法的参考轨迹y_r(k+P)以及目标函数J，形式如下：

y_r(k+P)＝β^Py_p(k)+(1-β^P)c(k)

J＝min(y_r(k+P)-y_mav(k+P)-e(k))²

e(k)＝y_pav(k)-y_mav(k)

其中，β为参考轨迹柔化系数，c(k)为k时刻的炉膛压力过程的设定值，e(k)为k时刻校正的误差值。

e.依据步骤d中的目标函数求解PI控制器中的参数，这里将炉膛压力过程的烟道挡板开度控制量u(k)进行变换：

u(k)＝u(k-1)+K_p(e₁(k)-e₁(k-1))+K_ie₁(k)

e₁(k)＝βy_p(k-1)+(1-β)c(k-1)-y_p(k)

进一步化简为：

u(k)＝u(k-1)+w(k)^ΤE(k)

w(k)＝[w₁(k),w₂(k)]^Τ

w₁(k)＝K_p+K_i,w₂(k)＝-K_p

E(k)＝[e₁(k),e₁(k-1)]^Τ

其中，K_p、K_i分别为炉膛压力过程PI控制器的比例、积分参数，e₁(k)为k时刻炉膛压力过程的参考轨迹值与实际输出值之间的误差，Τ为矩阵的转置符号。

结合上述式子，可以求得：

w(k)=(yr(k+P)-amPymav(k)-Km(1-amP)u(k-1))EKm(1-amP)ETE]]>

进一步可以得到：

K_p＝-w₂(k)

K_i＝w₁(k)-K_P

f.得到PI控制器的参数K_p、K_i以后构成控制量u(k)作用于裂解炉炉膛的烟道挡板开度阀门，u(k)＝u(k-1)+K_p(e₁(k)-e₁(k-1))+K_ie₁(k)。

g.在下一时刻，依照b到f中的步骤继续求解裂解炉炉膛压力过程PI控制器新的参数K_p、K_i，依次循环。

预测函数控制优化的废塑料炼油裂解炉炉膛压力控制方法专利购买费用说明

Q：办理专利转让的流程及所需资料

A：专利权人变更需要办理著录项目变更手续，有代理机构的，变更手续应当由代理机构办理。

1：专利变更应当使用专利局统一制作的“著录项目变更申报书”提出。

2：按规定缴纳著录项目变更手续费。

3：同时提交相关证明文件原件。

4：专利权转移的，变更后的专利权人委托新专利代理机构的，应当提交变更后的全体专利申请人签字或者盖章的委托书。

Q:专利著录项目变更费用如何缴交

A：（1）直接到国家知识产权局受理大厅收费窗口缴纳,（2）通过代办处缴纳，（3）通过邮局或者银行汇款,更多缴纳方式

Q：专利转让变更，多久能出结果

A：著录项目变更请求书递交后，一般1-2个月左右就会收到通知，国家知识产权局会下达《转让手续合格通知书》。