[发明专利]一种燃煤锅炉系统混合控制方法

权利要求书

1、一种燃煤锅炉系统混合控制方法，其特征在于该方法包括以下步骤：

(1)利用实时数据驱动的方法建立过程模型，具体方法是：

首先，建立燃煤锅炉实时运行数据库，通过数据采集装置采集实时过程运行数据，将采集的实时过程运行数据作为数据驱动的样本集合，表示为{X_i，y(i)}_i＝1^N，其中X_i表示第i组工艺参数的输入数据，y(i)表示第i组工艺参数的输出值；

其次，以该实时过程运行数据集合为基础建立基于最小二乘法的离散差分方程形式的局部受控自回归滑动平均模型：

y_L(k)＝Φ^TX，Φ＝[a′₁，a′₂，…，a′_n，b′₀，b′₁，…，b′_m-1]^T

X＝[y(k-1)，…，y(k-n)，u(k-d-1)，…，u(k-d-m)]^T

其中，y_L(k)表示当前时刻过程模型的工艺参数的输出值，X表示过程模型的工艺参数的过去时刻的输入和输出数据的集合，u(k)表示过程模型工艺参数对应的控制变量，Φ表示通过辨识得到的模型参数的集合，n，m，d+1分别为对应实际过程的输出变量阶次、输入变量阶次、时滞；

采用的辨识手段为：

Φk=Φk-1+K‾(k)[y(k)-ΦkTXk]]]>

K‾(k)=P(k-1)Xk[XkTP(k-1)Xk+γ]-1]]>

其中，γ为遗忘因子，为单位阵；

(2)采用典型的响应曲线法设计过程模型的比例积分微分控制器，具体方法是：

a.将过程模型的比例积分微分控制器停留在手动操作状态，操作拨盘使其输出有阶跃变化，由记录仪表记录过程模型的输出值，将过程模型输出值y_L(k)的响应曲线转换成无量纲形式y_L^*(k)，具体是：yL*(k)=yL(k)/yL(∞)]]>其中，y_L(∞)是过程模型的比例积分微分控制器的输出有阶跃变化时的过程模型输出y_L(k)的稳态值；

b.选取满足yL*(k1)=0.39,]]>yL*(k2)=0.63]]>的两个计算点k₁和k₂，，依据下式计算比例积分微分控制器所需要的参数K、T和τ：

K＝y_L(∞)/q

T＝2(k₁-k₂)

τ＝2k₁-k₂

其中，q为过程模型的比例积分微分控制器输出的阶跃变化幅度；

c.计算过程模型的比例积分微分控制器的参数，具体是：

K_c＝1.2T/Kτ

T_i＝2τ

T_d＝0.5τ

其中K_c，T_i，T_d分别为比例积分微分控制器的比例参数，积分参数，微分参数；

(3)设计预测比例积分微分控制器，具体步骤是：

d.将过程模型的比例积分微分控制器停留在自动操作状态，操作拨盘使其输入有阶跃变化，，由记录仪表记录实时过程的输出，将过程输出值y(k)的响应曲线转换成无量纲形式y^*(k)，具体是：y^*(k)＝y(k)/y(∞)

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

e.选取满足y(k₃)＝0.39，y(k₄)＝0.63的两个计算点，依据下式计算预测比例积分微分控制器所需要的参数K₁，T₁和τ₁：

K₁＝y(∞)/q₁

T₁＝2(k₃-k₄)

τ₁＝2k₃-k₄

其中，q₁为过程模型的比例积分微分控制器输入的阶跃变化幅度；

f.将步骤e得到的参数转化为离散差分方程形式的局部受控自回归滑动平均模型：

y(k)＝Φ^TX，Φ＝[a₁，a₂，…，a_n，b₀，b₁，…，b_m-1]^T

X＝[y(k-1)，…，y(k-n)，u(k-d-1)，…，u(k-d-m)]^T

其中，y(k)表示当前时刻过程模型的输出值，X表示过去时刻的输入和输出数据的集合，u(k)表示过程模型对应的控制变量，Φ表示通过转换得到的模型参数的集合，n，m，d+1分别为对应实际过程的输出变量阶次、输入变量阶次、时滞；

g.依据步骤f计算出的模型参数整定预测控制器的参数，具体方法是：

①建立多步最优预测输出Y^=Ypast+GU+E]]>

其中，

Y^=(y^(k+d+1/k),y^(k+d+2/k),...,y^(k+d+p/k))T]]>

Y_past＝(y_past(k+d+1)，y_past(k+d+2)，…，y_past(k+d+p))^T

U＝(Δu(k)，Δu(k+1)，…，Δu(k+p-1))^T

E=(y(k)-y^(k),y(k)-y^(k),...,y(k)-y^(k))T]]>

G=B1,0B2,0B1,00......Bp,0Bp-1,0...B1,0]]>

Y_past是依据过程过去时刻过程模型的输入输出数据计算得到的自由响应输出，G是预测模型的参数矩阵，u是过程控制增量的集合，E是预测模型的校正项；

②建立预测控制器的参考轨迹Y_r和目标函数J；

Y_r＝(y_r(k+d+1)，y_r(k+d+2)，…，y_r(k+d+p))^T

J=min{(Yr-Y^)T(Yr-Y^)+λ2UTU}]]>

③依据步骤②的目标函数得到当前的控制参数值u(k)

u(k)＝u(k-1)+q^T(Y_r-Y_past-E)

其中q^T是预测控制器参数矩阵的第一行向量。