>[0040] 球磨机磨矿过程的预测控制的结构图如图2所示。
[0041] 为了得到球磨机磨矿过程的预测输出，首先定义如下一组变量：
[0043] 公式中Ny是预测时域，NU(NU<Ny)是控制时域，Ny-般取值在15至25之间，Nu取 值一般在4至10之间。
，表示基于t时刻的值对t+i时刻的状态变量 预测值
，表示基于t时刻的值对t+j时刻的输出变量预测值。这里 假设大于控制时域以后的控制量不会发生改变，即u(t+j) =u(t+Nu-l) (j彡Nu)。这样就 可以得到在t时刻的Ny步预测的表示形式
[0045] 其中
[0050]公式（8)可以转变为模型预测输出的另外一种形式：
[0052] 定义如下控制增量△冲)和期望输出.M0 :
[0054]公式（14)中Au(i) =u(i)-u(i_l)，i=t, ? ? ?，t+Nu_l，它表示i时刻控制量的 增量；yji)，i=t+1,. . .，t+Np，表示i时刻的期望输出。
[0055] 采用如下的目标函数，该目标函数能够使实际的溢流粒度、排矿率与期望的溢流 粒度、排矿率保持一致：
[0057]公式中Q，札和R2是期望输出与模型预测输出偏差权重系数矩阵、输入u(t)权重 系数矩阵和Au(t)权重系数矩阵。
。Q的值过大会导致系统不稳定，Q的 值过小会导致稳态误差，一般在〇. 001至100之间；&的值过大会导致稳态误差，一般在 0. 0001至0. 1之间；当控制量变化过快时，可将馬调小，当控制量变化太小时，可将R2调大， 具体值根据实际情况确定。这三个系数矩阵都为对角阵。将公式（13)代入公式（15)得：
[0059]公式中
[0061] u0(t_l) = [u(t_l)Tu(t_l)T ...u(t_l)T]T，
[0063] I是2 X 2的单位矩阵；umin，umax分别表示控制量的最小值与最大值；y min，ymax分别 表示输出量的最小值与最大值；u(t-l)表示t-1时刻的控制量；Aumin,Aumax分别表示控制 变化率的最小值与最大值。
[0064] 对于公式（16)是典型的带线性不等式约束的二次规划问题，目前已有成熟的算 法进行求解，可采用有效集法、置信域法等算法进行求解出控制量在求的最优控制量 中只取第一部分u(t)作为实际的控制量即最优的给矿速度和给水流量。系统的控制 周期与采集数据的周期保持一致。
[0065] 给矿速度控制回路的原理图如图3所示，给矿速度采用PID控制，通过调整变频器 的频率来控制给料电机的转动，最终控制给矿速度。
[0066] 给水流量控制回路的原理图如图4所示，给水流量控制回路采用PID控制，通过调 整给水阀门的开度来调整最终的给水流量。
[0067] 给矿速度控制回路和给水流量控制回路是典型的PID控制，其PID参数可以通过 试凑法、临界比例度法等方法获得。
[0068] 本发明的具体实现步骤分成两大部分：离线计算与在线计算。
[0069] 离线部分主要是获得磨矿过程的模型的参数，具体如下：
[0070] 1、采集实际的给矿速度、实际的给水流量、磨矿的溢流粒度、球磨机的排矿率的历 史运行数据。
[0071] 2、选定不同的模型阶次ny，nu，通过最小二乘法求出球磨机磨矿过程ARX模型的系 数。
[0072] 3、根据公式（3)求出不同模型阶次ny，nu，下的AIC值，AIC最小对应的ny，nu为磨 矿过程的ARX模型的最佳阶次，对应的模型参数为磨矿过程模型的参数。
[0073] 在线计算部分的实现步骤如下：
[0074] 1、首先在工业控制计算机中设定控制参数Q，&，R2，预测控制长度Ny，Nu，设定控制 方法相关的上下限值umin、umax、ymin、ymax、Aumin、Aumax，设定给水流量控制回路中的PID参 数，设定给矿速率控制回路的PID参数。
[0075] 2、工业控制计算通过PLC采集当前的磨矿溢流粒度和球磨机的排矿率，当前实际 的给矿速度和实际的给水流量，根据公式（8)预测磨矿过程未来的磨矿溢流粒度与排矿 率，然后获得公式（16)中的系数，最后通过有效集法、置信域法等算法求解出最优的控制 量，即最优的给矿速度和给水流量。
[0076] 3、工业计算机通过网络将最优的控制量发送给PLC，PLC中将上述最优控制量中 的给矿速度作为给矿控制回路中的给定值，将最优控制量中的给水流量作为给水控制回路 中的给定值。
[0077] 4、皮带称称量系统控制实际的给矿速度与给矿速度的给定值一致；给水流量控制 回路通过PID控制阀门开度使实际的给水流量与给水流量的给定值一致。
[0078] 以上步骤2至步骤4按照一定的控制周期不断的循环进行。
【主权项】
1. 一种球磨机控制方法，其特征在于，求解以下目标函数J的最小值对应的最优控制 量将的第一部分U (t)作为最优的给矿速度和给水流量，从而控制球磨机的溢 流粒度和排矿率：上述公式中，A⑴="(4 "(y + V…+ 丨） ；I是2X2的单 位矩阵；umin，11_分别表示控制量u⑴的最小值与最大值，u (t) = [u i (t) U2⑴]T，U1 (t)为t时刻的球磨机的给矿速度，单位为kg/min，u2 (t)为t时刻的球磨机的 给水流量，单位为m3/h ;ymin，7_分别表示输出量y⑴的最小值与最大值，y⑴= [yi⑴y2⑴]T，Y1 (t)为t时刻的磨矿的溢流粒度，y2 (t)为t时刻的排矿率，单位 为kg/min ; Δ umin, Δ 11_分别表示控制增量Δ u (m)的最小值与最大值，Δ u (m)= u (m)-U (m-1)，m = t,. . .，t+Nu_l，Nu是控制时域；的未知元素利用ARX模型和最小二乘法辨识求得，，./ = 1，2···""，ny, nu分别为 y ⑴、 u⑴对应的阶次，e⑴是2 X 1的白噪声，且与y⑴和u⑴都线性无关；Q是期望输出与球 磨机磨矿过程ARX模型预测输出偏差权重系数矩阵，&是u(t)权重系数矩阵，R2是Au(m) 权重系数矩阵T;免(0为t时刻的期望输出。2.根据权利要求1所述的球磨机控制方法，其特征在于，Φ，A i 和Bj中的参数辨识方法如下：将需要辨识的参数集合成一向量 5 = [φ 4 4…4,.尽4…見/，并将t-i时刻的输入输出数据整合成向 量挪-1) = [1 .v(/ - l) >'(/-2)…"(/-I) w(/ - l)…Μ(?-"")] ，定义 Φ.\ =[供(Λ-1)，ρ⑷，…供(，V-.l)y , Yn= [y (k)，y (k+1)，· · · y (N)] T，k = max(nu，ny)+l，则根据 最小二乘法，向量3的估计值为5 = ((1/(.(1>、）>〈\，将6的估计值作为球磨机磨矿过程 ARX模型的最佳参数，即得到Φ，化和B j参数。。3.根据权利要求2所述的球磨机控制方法，其特征在于，n u，ny的求解方法为：先选择 多组nu，ny，根据设定函数AIC = Nlog 〇2+2d，其中〇2为噪声方差，N是采集的溢流粒度、排 矿率、给矿速度、给矿流量辨识数据的组数，N>>max(ny，nu)，d = 4*nu+4*ny+2 ;求得不同nu， ny组合下的AIC值，当AIC最小时对应的n u，ny，为最佳的nu，ny。。
【专利摘要】本发明公开了一种球磨机控制方法，采集球磨机在过去运行的数据包括实际的给矿速度，实际的给水流量，实际的磨矿的溢流粒度和排矿率，将采集的这些数据通过PLC最终通过工业网络传输给工业控制计算机，在工业控制计算机中完成控制计算，计算出最佳给矿速度和给水流量，工业控制计算机通过通网络将这些值发送给PLC，PLC控制实际的给矿速度、给水流量与给定的给矿速度、给水流量保持一致。本发明能够实现球磨机的自动控制，避免人工控制产生的不确定因素导致的磨矿效果的不确定；本发明采用的预测控制能够同时控制给矿速度和给水流量，避免了通过溢流粒度偏差分别进行PID控制控制给矿速率和给水流量存在的相互干扰的情况。
【IPC分类】G05B19/404
【公开号】CN104898563
【申请号】CN201510217335
【发明人】魏吉敏, 黄应龙, 梁军, 王宇星, 杨鸿波
【申请人】长沙有色冶金设计研究院有限公司
【公开日】2015年9月9日
【申请日】2015年4月30日