Dual-rate Interval Control of Pump Pool Level and Feeding Pressure During Regrinding
-
摘要: 赤铁矿再磨过程是以矿浆泵频率为输入、以给矿压力为内环输出、以泵池液位为外环输出的强非线性串级工业过程.当赤铁矿粒度分布大范围变化,导致一段磨矿与磁选矿浆流量和再磨排矿流量频繁波动,使泵池液位频繁波动,造成内外环频繁波动.本文将上述动态特性变化用未建模动态来描述,通过设计消除前一时刻未建模动态补偿信号叠加于基于线性模型设计的反馈控制器,采用一步最优前馈控制律和提升技术,提出了泵池液位和给矿压力双速率区间控制算法.给出了所提控制算法的稳定性和收敛性分析,通过半实物仿真实验表明所提出的控制算法可以将处于频繁随机波动的泵池液位和给矿压力变化率控制在目标值范围内.Abstract: The regrinding process of hematite is a strong nonlinear cascade process with frequency of slurry pump as the input, feeding pressure as the inner loop output and sump level as the outer loop output. During its operation, frequent fluctuations of first stage grinding, magnetic separation slurry flowrate and discharge regrinded flowrate generated from the distribution of hematite particles which changes in a large range will cause large variations in the sump level and frequent fluctuation between inner and outer loop. In this paper, a novel dual-rate interval controller of sump level and feeding pressure is developed by representing the variations of dynamics characters as the unmodeled dynamics. In the proposed controller design, a compensation signal is constructed and added onto the control signal obtained from the linear deterministic model based on feedback control design. Such a compensation signal aims at eliminating the previous moment unmodeled dynamics. Meantime, the stability and convergence analysis of the proposed algorithm is given. A simulated experiment on the hardware-in-the-loop simulation system is conducteded, and it is shown that the sump level and the changing rate of feeding pressure in frequent random fluctuations can be well controlled inside their targeted ranges under the proposed control algorithm.
-
Key words:
- Regrinding process /
- unmodeled dynamic compensation /
- PI control /
- dual-rate control
-
选矿再磨过程是将一段磨矿浆和不合格的再磨矿浆进行再次磨矿分级的过程.在再磨的生产过程中, 稳定泵池液位是保证再磨过程的安全运行的重要条件[1], 稳定旋流器给矿压力是保证再磨分级效率的先决条件[2], 因此再磨过程高效稳定运行必须保证泵池液位和旋流器的给矿压力都要达到控制目标范围内.
国外磨矿基本采用棒磨开路的定量给矿和配比给水的磨矿工艺[3].由于其矿石性质较国内稳定, 使一段磨矿浆流量稳定, 泵池液位波动小, 采用定值闭环控制方式对泵池液位控制也不会使旋流器的给矿压力发生较大波动.文献[4]针对铜矿的再磨过程是以泵池补加水流量和矿浆泵的频率为输入, 以泵池液位和旋流器给矿浓度为输出, 采用多变量控制方法对泵池液位和给矿浓度进行定值控制.文献[5]针对镍矿再磨过程以矿浆泵流量为输入, 以泵池液位为输出, 采用预测控制方法对泵池液位进行定值控制.上述控制方法[4-5]通过工业应用实验说明文献[6-7]针对金矿的再磨过程方法的有效性.以矿浆泵频率为输入, 以泵池液位为输出, 分别采用PID控制方法和模糊控制方法对泵池液位进行定值控制.文献[8]针对金矿磨矿过程, 以泵池出口流量为输入, 以泵池液位为输出, 采用分数阶PI控制方法对泵池液位进行定值控制.文献[9-10]针对金矿磨矿过程, 以矿浆泵频率和泵池补加水流量为输入, 以泵池液位和旋流器给矿浓度为输出, 分别采用非线性模型预测控制方法和鲁棒非线性模型预测方法对泵池液位进行定值控制.文献[11]针对铂矿磨矿过程, 以泵池出口流量和泵池补加水流量为输入, 以泵池液位和旋流器给矿浓度为输出, 采用多变量控制方法实现对泵池液位定值控制.上述控制方法[6-11]通过仿真实验说明方法的有效性.
我国赤铁矿不仅品位低而且具有嵌布粒度细且粒度分布不均等特点, 需要两段闭路磨矿才能获得粒度合格的矿浆.由于一段磨矿矿浆流量波动大和再磨不合格矿浆-排矿流量等随机干扰影响, 导致泵池液位发生较大波动, 如果采用液位定值闭环控制方法[4-11], 会造成旋流器给矿压力较大波动, 使给矿压力内环的闭环系统处于动态之中.文献[12]提出了泵池液位区间与给矿压力模糊切换控制, 外环采用静态模型设定、模糊补偿器和切换机制组成的控制方法.由于再磨过程产生的不合格矿浆-排矿流量与给矿浓度相关, 给矿浓度波动引起排矿流量波动, 为了提高再磨过程控制系统的动态性能, 文献[13]在文献[12]基础上设计了消除给矿浓度波动的区间智能控制器.当赤铁矿的粒度分布大范围变化时, 使一段磨矿与磁选矿浆流量和再磨排矿流量频繁波动, 导致给矿压力内环控制器设定值频繁波动, 使给矿压力内环闭环控制系统处于动态之中.由于上述控制方法[12-13]的外环控制器设计没有考虑内环闭环控制系统的动态特性, 基于规则的切换机制难以给出合适的给矿压力控制器的设定值, 难以建立控制系统的稳定性和收敛性分析, 因此, 难以将泵池液位和给矿压力变化率在运行时间内控制在目标值范围内.
本文将上述随机干扰造成的动态特性未知变化用未建模动态来描述, 设计了基于消除前一时刻未建模动态补偿信号的PI控制器, 采用一步最优前馈控制律和提升技术, 提出了泵池液位和给矿压力双速率区间控制算法, 给出了所提算法的稳定性和收敛性分析, 采用工业现场实际数据的半实物仿真实验, 结果表明了所提算法的有效性.
1. 控制问题描述
1.1 赤铁矿再磨过程简介
赤铁矿再磨过程的结构图如图 1所示.由磁选机、球磨机、水力旋流器、泵池、泵池补加水阀门、矿浆泵和污水池、污水泵、冲洗水阀门等组成.再磨生产过程的来料为一段磨矿和磁选产生的矿浆, 最终产品为旋流器溢流矿浆.
一段磨矿矿浆和磁选矿浆以流量 $q_1(t)$ 进入到泵池, 通过控制泵池补加水阀门开度 $u_1(t)$ 使得进入旋流器的矿浆浓度 $D_H(t)$ 达到工艺规定的目标值, 稀释后的矿浆在矿浆泵作用下以一定的给矿压力 $y_1(t)$ 输送给旋流器, 在离心力的作用下分级出粒度合格的溢流矿浆, 粒度不合格的矿浆则形成沉砂返回球磨机进行再磨, 球磨机排矿以流量 $q_2(t)$ 进入到泵池, 并和一段磨矿与磁选矿浆流量 $q_1(t)$ 、不定期工作的污水与冲洗水流量 $q_3(t)$ 以及泵池补加水流量 $q_4(t)$ 充分混合后在泵池矿浆泵的输送下, 以一定的频率 $u(t)$ 从泵池抽出, 从而维持泵池液位 $y_2(t)$ 在工艺要求的安全范围内.
1.2 赤铁矿再磨过程动态模型分析
根据文献[14], 以矿浆泵频率 $u(t)$ 为输入, 以泵池液位 $y_2(t)$ 为输出的动态模型为
$ \begin{align} \begin{aligned} \frac{{{\rm d}{y_2}}}{{{\rm d}t}} &= \frac{1}{A} \left[\vphantom{\frac{{{K_d}}}{{\sqrt {{K_p}} }}\sqrt {{y_1}(t)}} {q_1}(t) + {q_2}({y_1}(t), {D_H}, {f_d}, {\phi _H})+ \right. \\ & \left. {q_3}(t) + {q_4}(t)-\frac{{{K_d}}}{{\sqrt {{K_p}} }}\sqrt {{y_1}(t)} \right] \end{aligned} \end{align} $
(1) $ \begin{align} \frac{{{\rm d}{y_1}}}{{{\rm d}t}} = - \frac{{{y_1}(t)}}{\tau } + \frac{{{K_p}}}{\tau }{u^2}(t) \end{align} $
(2) 式(1) 和式(2) 模型参数描述, 如表 1所示.
表 1 模型参数表Table 1 Model parameters变量 描述 $q_1(t)$ 一段磨矿矿浆和磁选矿浆流量 ${q_2}({y_1}(t),{D_H},{f_d},{\phi _H}) $ 再磨排矿流量 $q_3(t)$ 污水与冲洗水流量 $q_4(t)$ 泵池补加水流量 $D_H(t)$ 旋流器给矿浓度 $f_d(t)$ 旋流器给矿粒度分布 $\phi_H$ 旋流器结构参数 $A$ 泵池横截面积 $K_p$ 矿浆泵比例系数 $K_d$ 矿浆泵流量比例系数 $\tau$ 惯性时间常数 由式(1) 和式(2) 可知, 再磨过程的泵池液位和给矿压力是以矿浆泵频率 $u(t)$ 为输入, 以给矿压力 $y_1(t)$ 为内环输出, 以泵池液位 $y_2(t)$ 为外环输出的强非线性串级过程.其中, 再磨过程不合格矿浆-排矿流量 ${q_2}({y_1}(t), {D_H}, {f_d}, {\phi _H}) $ 是与旋流器给矿浓度 $D_H(t)$ 、给矿压力 $y_1(t)$ 、旋流器给矿粒度分布 $f_d(t)$ 和旋流器结构参数 $\phi_H$ 相关的未知非线性函数.当赤铁矿粒度分布大范围变化时, 使旋流器给矿浓度 $D_H(t)$ 和粒度分布 $f_d(t)$ 发生波动, 造成 ${q_2}({y_1}(t), {D_H}, {f_d}, {\phi _H})$ 波动, 同时一段磨矿与磁选矿浆流量 $q_1(t)$ 和污水与冲洗水流量 $q_4(t)$ 的干扰, 使矿浆液位频繁波动, 导致给矿压力控制器设定值波动, 使给矿压力内环闭环控制系统处于动态之中, 泵池液位外环和给矿压力内环相互影响, 为了消除未建模动态变化的影响, 利用再磨过程运行在工作点附近的特点, 将式(1) 和式(2) 在工作点处线性化可得:
$ \begin{align} {A_2}({z^{ - 1}}){y_2}(T + 1) = {B_2}({z^{ - 1}}){y_1}(T) + {v_2}(T) \end{align} $
(3) 式中, ${A_2}({z^{ - 1}}) = 1 + {a_{21}}{z^{ - 1}}$ , ${B_2}({z^{ - 1}}) = {b_{20}}$ .
$ \begin{align} {A_1}({z^{ - 1}}){y_1}(k + 1) = {B_1}({z^{ - 1}})u(k) + {v_1}(k) \end{align} $
(4) 式中, ${A_1}({z^{ - 1}}) = 1 + {a_{11}}{z^{ - 1}}$ , ${B_1}({z^{ - 1}}) = {b_{10}}$ .
泵池液位 $y_2(T)$ 和给矿压力 $y_1(k)$ 具有不同的时间尺度, 其中, 泵池液位采样周期为 $T$ , 给矿压力采样周期为 $k$ .利用输入输出数据通过实验确定 ${A_2}({z^{ - 1}})$ 、 ${B_2}({z^{ - 1}})$ 、 ${A_1}({z^{ - 1}})$ 和 ${B_1}({z^{ - 1}})$ 的阶次和参数, 未建模动态 ${v_2}(T)$ 和 ${v_1}(k)$ 为高阶非线性项, 且 ${v_2}(T)$ 和 ${v_1}(k)$ 有界, 即: $|{v_2}(T)|\leq N$ ; $|{v_1}(k)|\leq M$ .
再磨过程的控制目标是:在所有运行时间内, 满足矿浆泵的约束下, 即: ${u_{\min}}\leq u(k)\leq {u_{\max}}$ , 将泵池液位 $y_2(T)$ 和给矿压力变化率 $\Delta {y_1}(k)$ 控制在目标值范围内, 即:
$ \begin{align} \left| {{y_2}(T) - {y_{2sp}}(T)} \right| \le {\delta _1}, \;\;\;\;{\kern 1pt} 1 \le T \le \infty \end{align} $
(5) $ \begin{align} \left| {{y_1}(k) - {y_1}(k - 1)} \right| < {\delta _2}, \;\;\;\;1 \le k \le \infty \end{align} $
(6) 其中, $y_{2sp}(T)$ 为泵池液位的目标值, $\delta _1$ 为泵池液位波动幅度的上限值; $\delta _2$ 为给矿压力波动的上限值.
2. 泵池液位和给矿压力双速率区间控制方法
由于泵池液位和给矿压力是串级过程, 因此首先设计给矿压力内环控制器.为了克服给矿压力内环闭环系统动态波动对泵池液位外环的影响, 采用提升技术[17], 将给矿压力内环闭环控制系统的动态特性引入到泵池液位模型式(3), 得到以给矿压力设定值 $y_1(T)$ 为输入, 以泵池液位 $y_2(T)$ 为输出的采样周期为外环采样周期的动态模型, 以此模型设计泵池液位外环控制器, 从而实时产生给矿压力内环控制的设定值 $y_{1sp}(T)$ , 通过给矿压力内环控制器, 使给矿压力 $y_1(T)$ 跟踪设定值 $y_{1sp}(T)$ , 实现对泵池液位的控制.
$q_1(t)$ 、 $q_2(t)$ 、 $q_3(t)$ 和 $q_4(t)$ 的波动引起的泵池液位和给矿压力的特性变化, 由未建模动态 $v_2(T)$ 和 $v_1(k)$ 表示.由于 $v_2(T)$ 和 $v_1(k)$ 始终处于动态波动之中, 导致控制器积分作用失效, 因此必须设计消除 $v_2(T)$ 和 $v_1(k)$ 影响的控制器.虽然未建模动态 $v_2(T)$ 和 $v_1(k)$ 未知, 但是 $v_2(T)$ 和 $v_1(k)$ 可以表示为 $v_2(T)=v_2(T-1)+\Delta v_2(T)$ 和 $v_1(k)=v_1(k-1)+\Delta v_1(k)$ .因此利用前一时刻的未建模动态, 设计消除前一时刻未建模动态 $v_2(T-1)$ 和 $v_1(k-1)$ 的影响的补偿信号, 叠加到基于确定线性模型设计的反馈控制器, 使闭环系统跟踪误差变小, 将泵池液位和给矿压力变化率在运行时间内控制在目标值范围内.
被控对象模型式(3) 和式(4) 可以表示为
$ \begin{align} {A_2}({z^{ - 1}})&{y_2}(T + 1) = {B_2}({z^{ - 1}}){y_1}(T)+\nonumber\\ & {v_2}(T - 1) + \Delta {v_2}(T) \end{align} $
(7) $ \begin{align} {A_1}({z^{ - 1}})& {y_1}(k + 1) ={B_1}({z^{ - 1}})u(k)+\nonumber\\ & {v_1}(k - 1) + \Delta {v_1}(k) \end{align} $
(8) 由式(7) 和式(8) 可得 $v_2(T-1)$ 和 $v_1(k-1)$ ,
$ \begin{align} {v_2}(T - 1) &= {y_2}(T) + A_2^*({z^{ - 1}}){y_2}(T)- \nonumber\\ & {B_2}({z^{ - 1}}){y_1}(T - 1)=\nonumber\\ & {y_2}(T) - y_2^*(T) \end{align} $
(9) $ \begin{align} {v_1}(k - 1) &= {y_1}(k) + A_1^*({z^{ - 1}}){y_1}(k) -\nonumber\\ & {B_1}({z^{ - 1}})u(k - 1) =\nonumber\\ & {y_1}(k) - y_1^*(k) \end{align} $
(10) 式中, $A_2^*({z^{ - 1}})={A_2}({z^{ - 1}})-1$ , $A_1^*({z^{ - 1}})={A_1}({z^{ - 1}})-1$ .
$ \begin{align} y_2^*(T) = - A_2^*({z^{ - 1}}){y_2}(T) + {B_2}({z^{ - 1}}){y_1}(T - 1) \end{align} $
(11) $ \begin{align} y_1^*(k) = - A_1^*({z^{ - 1}}){y_1}(k) + {B_1}({z^{ - 1}})u(k - 1) \end{align} $
(12) 式(11) 和式(12) 为泵池液位和给矿压力的控制器驱动模型[18]的输出.采用式(9)~(12) 可以求得 $v_2(T-1)$ 和 $v_1(k-1)$ , 同时该方法在考虑当一段磨矿与磁选矿浆流量和再磨排矿流量频繁波动的基础上, 采用一步最优前馈补偿律设计PI控制器, 以此设计未建模动态补偿的PI控制器.
本文提出的双速率区间控制结构图如图 2所示.给矿压力内环控制器和泵池液位外环控制器设计如下.
2.1 未建模动态补偿的给矿压力PI控制器
基于未建模动态补偿的给矿压力PI控制器的结构如图 3所示, 其中, $u_1(k)$ 为基于给矿压力模型式(4) 中的确定线性部分模型设计PI控制器输出, $u_2(k)$ 为前一时刻未建模动态 $v_1(k-1)$ 补偿器的输出, 即
$ \begin{align} u(k) = {u_1}(k) + {u_2}(k) \end{align} $
(13) 以式(8) 的确定线性部分模型设计的PI控制器为
$ \begin{align} {H_1}({z^{ - 1}}){u_1}(k) = {G_1}({z^{ - 1}}){e_1}(k) \end{align} $
(14) 式中, ${H_1}({z^{ - 1}})=1-z^{ - 1}$ 、 ${G_1}({z^{ - 1}})=g_{10}-g_{11}z^{ - 1}$ , $g_{10}$ 和 $g_{11}$ 为PI控制参数, $e_1(k)$ 为跟踪误差, 即: $e_1(k)=y_{1sp}(k)-y_1(k)$ .
未建模动态 ${v_1}(k - 1)$ 补偿器为
$ \begin{align} {u_2}(k) = - {K_1}({z^{ - 1}}){v_1}(k - 1) \end{align} $
(15) 式中, ${K_1}({z^{ - 1}})$ 为补偿器的参数.
采用一步最优前馈补偿律来设计 ${G_1}({z^{ - 1}})$ 和 ${K_1}({z^{ - 1}})$ 的参数, 将式(14) 中的 ${u_1}(k)$ 和式(15) 中的 ${u_2}(k)$ 代入式(13) 中得到 $u(k)$ 为
$ \begin{align} {H_1}({z^{ - 1}})u(k) &= {G_1}({z^{ - 1}})\left[{{y_{1sp}}(k)-{y_1}(k)} \right]- \nonumber\\ &{H_1}({z^{ - 1}}){K_1}({z^{ - 1}}){v_1}(k - 1) \end{align} $
(16) 引入下列性能指标[19]:
$ \begin{align} J_1=\, &\min \left[P_1(z^{-1})y_1(k+1)-R_1(z^{-1})y_{1sp}(k)+\right.\nonumber\\ &\left.Q_1(z^{-1})u(k)+K(z^{-1})v_1(k)\right]^2 \end{align} $
(17) 式中, $P_1(z^{-1})$ 、 $R_1(z^{-1})$ 、 $Q_1(z^{-1})$ 和 $K(z^{-1})$ 均为关于 $z^{-1}$ 的加权多项式.
引入广义输出 ${\phi _1}(k + 1)$ 为
$ \begin{align} {\phi _1}(k + 1) = {P_1}({z^{ - 1}}){y_1}(k + 1) \end{align} $
(18) 定义广义理想输出 $\phi _1^*(k + 1)$ 为
$ \begin{align} \phi _1^*(k + 1) =\, & {R_1}({z^{ - 1}}){y_{1sp}}(k) - {Q_1}({z^{ - 1}})u(k)- \nonumber\\ & K({z^{ - 1}}){v_1}(k) \end{align} $
(19) 定义式(17) 中的 ${P_1}({z^{ - 1}})$ 为
$ \begin{align} {P_1}({z^{ - 1}}) = {A_1}({z^{ - 1}}) + {z^{ - 1}}{G_1}({z^{ - 1}}) \end{align} $
(20) 由式(8) 和式(20) 可得
$ \begin{align} {P_1}({z^{ - 1}}){y_1}(k + 1) =\, & {G_1}({z^{ - 1}}){y_1}(k)+ {B_1}({z^{ - 1}})u(k)+\nonumber\\ & {v_1}(k - 1) + \Delta {v_1}(k) \end{align} $
(21) 将式(21) 代入式(17) 中, 使 $J_{\min}=\Delta v_1(k)$ 可得带有未建模动态补偿的一步最优控制律为
$ \begin{align} &\left[{{B_1}({z^{-1}}) + {Q_1}({z^{-1}})} \right]u(k) = {R_1}({z^{ - 1}}){y_{1sp}}(k)-\nonumber\\ & \qquad {G_1}({z^{ - 1}}){y_1}(k) -\left[{1 + K({z^{-1}})} \right]{v_1}(k - 1) \end{align} $
(22) 由式(16) 和式(22) 可得 ${Q_1}({z^{ - 1}})$ 、 ${R_1}({z^{ - 1}})$ 、 $K({z^{ - 1}})$ 为
$ \begin{align} \left\{ \begin{array}{l} {Q_1}({z^{ - 1}}) = {H_1}({z^{ - 1}}) - {B_1}({z^{ - 1}})\\ {R_1}({z^{ - 1}}) = {G_1}({z^{ - 1}})\\ K({z^{ - 1}}) = {H_1}({z^{ - 1}}){K_1}({z^{ - 1}}) - 1 \end{array} \right. \end{align} $
(23) 将式(22) 和式(23) 代入到给矿压力被控对象式(8) 中得到给矿压力闭环系统方程为
$ \begin{align} \big[(1 + &{a_{11}}{z^{-1}})(1-{z^{-1}})+\nonumber\\ & {z^{ - 1}}{b_{10}}({g_{10}} + {g_{11}}{z^{ - 1}}) \big]{y_1}(k + 1) = \nonumber\\& {b_{10}}({g_{10}} + {g_{11}}{z^{ - 1}}){y_{1sp}}(k) + (1 - {z^{ - 1}})\Delta {v_1}(k)+\nonumber\\ & (1 - {z^{ - 1}})\left[{1-{b_{10}}{K_1}({z^{-1}})} \right]{v_1}(k - 1) \end{align} $
(24) 选择 $g_{10}$ 和 $g_{11}$ 满足:
$ \begin{align} &(1 + {a_{11}}{z^{ - 1}})(1 - {z^{ - 1}})+ \nonumber\\ & \qquad {z^{ - 1}}{b_{10}}({g_{10}} + {g_{11}}{z^{ - 1}}) \ne 0, \;\left| z \right| > 1 \end{align} $
(25) 由式(24) 可知, 为实现对 ${v_1}(k - 1)$ 的动态和静态补偿, 选择 ${K_1}({z^{ - 1}})$ 使 $1-b_{10}{K_1}({z^{ - 1}}) =0$ , 即:
$ \begin{align} {K_1}({z^{ - 1}}) = \frac{1}{{{b_{10}}}} \end{align} $
(26) 未建模动态补偿的给矿压力PI控制器输出 $u(k)$ 可以由式(13) 求得, 其中, $u_1(k)$ 由式(14) 和式(25) 求得, $u_2(k)$ 由式(15) 和式(26) 求得.
2.2 未建模动态补偿的泵池液位PI控制算法
将式(16) 代入式(8) 中, 得到给矿压力闭环控制系统为
$ \begin{align} {T_1}({z^{ - 1}}){y_1}(k + 1) =\, & {B_1}({z^{ - 1}}){G_1}({z^{ - 1}}){y_{1sp}}(k)+ \nonumber\\ & {H_1}({z^{ - 1}})\Delta {v_1}(k) \end{align} $
(27) 式中, ${T_1}({z^{ - 1}}) = (1 + {a_{11}}{z^{ - 1}})(1 - {z^{ - 1}}) + {z^{ - 1}}{b_{10}}({g_{10}} + {g_{11}}{z^{ - 1}})\; = 1 + {t_1}{z^{ - 1}} + {t_2}{z^{ - 2}}$
一段磨矿与磁选矿浆流量和再磨排矿流量频繁波动, 导致泵池液位 $y_2(T)$ 发生频繁波动, 造成给矿压力内环控制器设定值 $y_{1sp}(T)$ 波动, 给矿压力闭环控制系统式(27) 处于动态之中, 泵池液位外环控制器设计需要建立基于给矿压力内环闭环控制系统特性的泵池液位外环动态模型, 采用提升技术[19]将式(27) 转化为状态方程形式, 选择状态变量: $x_1(k)=-t_2y_1(k-1)+b_{10}g_{11}y_{1sp}(k-1)-\Delta v_1(k-1)$ 、 $x_2(k)=y_1(k)$ .将式(27) 表示为
$ \begin{align} \left\{ \begin{array}{l} {{{x}}}(k + 1) = {{{Ax}}}(k) + {{{B}}}{y_{1sp}}(k) + {{{D}}}\Delta {v_1}(k)\\ {y_1}(k) = {\textit{Cx}}(k) \end{array} \right. \end{align} $
(28) 式中, .
由式(28) 可得:
$ \begin{align} {{x}}(5k + i + 1) = \, &{{Ax}}(5k + i) + {{B}}{y_{1sp}}(5k + i)+ \nonumber\\ & {{D}}\Delta {v_1}(5k + i) \end{align} $
(29) 当 $i=0, 1, 2, \cdots, 4$ 时, 将式(29) 进行迭代, 可得:
$ \begin{align} {{x}}(5k + 5) = \, &{{{A}}^5}{{x}}(5k) + \sum\limits_{i = 0}^4 {{{{A}}^{4 - i}}} ({{B}}{y_{1sp}}(5k + i) +\nonumber\\ & {{D}}\Delta {v_1}(5k + i)) \end{align} $
(30) 利用零阶保持器, 使
$ \begin{align} {y_{1sp}}(5k) = {y_{1sp}}(5k + i) \end{align} $
(31) 其中 $i=0, 1, 2\cdots, 4$ .
由 $T=5k$ , 状态方程式(28) 和式(29) 中的输出方程分别转化为
$ \begin{align} {{x}}(T + 1) = \, &{{{A}}^5}{{x}}(T) + \nonumber\\ &\left(\mathop \sum \limits_{i = 0}^4 {{{A}}^i}{{B}}\right){y_{1sp}}(T) + \Delta {{\tilde v}}(T) \end{align} $
(32) 式中, $\Delta {{\tilde v}}(T) = \sum_{i = 0}^4 {{{{A}}^{4 - i}}} {{D}}\Delta {v_1}(5k + i)$ .
$ \begin{align} {y_1}(T) = {{Cx}}(T) \end{align} $
(33) 将式(32) 和式(33) 转换成输入输出的形式, 即:
$ \begin{align} {T_0}({z^{ - 1}}){y_1}(T) = {D_0}({z^{ - 1}}){y_1}_{sp}(T) + \Delta {\tilde v_1}(T) \end{align} $
(34) 式中
$ {T_0}({z^{ - 1}}) = \det ({{I}} - {z^{ - 1}}{{{A}}^5}) $
$ {D_0}({z^{ - 1}}) = {{C}}{\rm{adj}}({{I}} - {z^{ - 1}}{{{A}}^5})\left(\mathop \sum \limits_{i = 0}^4 {{{A}}^i}{{B}}\right) $
$ \Delta {\tilde v_1}(T) = {{C}}{\rm{adj}}({{I}} - {z^{ - 1}}{{{A}}^5})\Delta {{\tilde v}}(T) $
又因
$ \begin{align} \begin{aligned} \left| {\Delta {v_1}(k)} \right| = \, &\left| {{v_1}(k) - {v_1}(k - 1)} \right| \le \left| {{v_1}(k)} \right|+ \\ & \left| {{v_1}(k - 1)} \right| \end{aligned} \end{align} $
(35) 由 $\left| {{v_1}(k)} \right|\le M$ 和式(35) 可得:
$ \begin{align} \left| {\Delta {v_1}(k)} \right| \le 2M \end{align} $
(36) 由式(32)、式(34) 和式(36) 得, ${\tilde v_1}(T)$ 有界, 即:
$ \begin{align} \left| {{{\tilde v}_1}(T)} \right| \le {M_1} \end{align} $
(37) 式中, $M_1$ 为 ${\tilde v_1}(T)$ 的上界.
将式(34) 代入式(3) 中, 得到以给矿压力设定值 $y_{1sp}(T)$ 为输入, 以泵池液位 $y_2(T+1)$ 为输出的动态模型, 即:
$ \begin{align} A({z^{ - 1}}){y_2}(T + 1) = B({z^{ - 1}}){y_{1sp}}(T) + v(T) \end{align} $
(38) 式中
$ \begin{align} & A({{z}^{-1}})={{T}_{0}}({{z}^{-1}}){{A}_{2}}({{z}^{-1}})= \\ & \ \ \ \ \ \ \ \ \ \ 1+{{a}_{1}}{{z}^{-1}}+{{a}_{2}}{{z}^{-2}}+{{a}_{3}}{{z}^{-3}} \\ \end{align} $
$ B({z^{ - 1}}) = {D_0}({z^{ - 1}}){B_2}({z^{ - 1}}) = {b_0} + {b_1}{z^{ - 1}} $
$ v(T) = {B_2}({z^{ - 1}})\Delta {\tilde v_1}(T) + {z^{ - 1}}{T_0}({z^{ - 1}}){v_2}(T) $
$v(T)$ 由给矿压力内环未建模动态变化率 $\Delta v_1(k)$ 和泵池液位外环未建模动态 $v_2(T)$ 组成.由 $v_2(T)$ 和 $v_1(k)$ 有界, 可知 $v(T)$ 有界, 即: $\left| {v(T)} \right|\le W$ .因此通过设计对 $v(T-1)$ 的补偿信号, 消除 $v(T-1)$ 对泵池液位 $y_2(T)$ 的影响.
采用未建模动态补偿的给矿压力PI控制器的设计方法, 设计泵池液位控制器, 其输出 $y_{1sp}(T)$ 为
$ \begin{align} {y_{1sp}}(T) = {y_{1sp1}}(T) + {y_{1sp2}}(T) \end{align} $
(39) $y_{1sp1}(T)$ 为确定线性部分模型设计的PI控制器输出, 即:
$ \begin{align} {H_2}({z^{ - 1}}){y_{1sp1}}(T) = {G_2}({z^{ - 1}}){e_2}(T) \end{align} $
(40) 式中, ${H_2}({z^{ - 1}})=1-z^{ - 1}$ 、 ${G_2}({z^{ - 1}})=g_{20}+g_{21}({z^{ - 1}})$ , $g_{20}$ 和 $g_{21}$ 为PI控制器参数、 $e_2(T)$ 为跟踪误差, $e_2(T)=y_{2sp}(T)-y_2(T)$ .
$y_{1sp2}(T)$ 为 $v(T-1)$ 的补偿器输出, 即:
$ \begin{align} {y_{1sp2}}(T) = - {K_2}({z^{ - 1}})v(T - 1) \end{align} $
(41) 式中, $v(T - 1) = A({z^{ - 1}}){y_2}(T) - B({z^{ - 1}}){y_{1sp}}(T - 1)$ .
采用一步最优前馈补偿律来设计 $G_2({z^{ - 1}})$ 和 $K_2({z^{ - 1}})$ 的参数, 将式(40) 中的 $y_{1sp1}(T)$ 和式(41) 中的 $y_{1sp2}(T)$ 代入式(39) 中得到 $y(T)$ 为
$ \begin{align} {y_{1sp}}(T) =\, & {y_{1sp}}(T - 1)+\nonumber\\ & ({g_{20}} + {g_{21}}{z^{ - 1}})\left[{{y_{2sp}}(T)-{y_2}(T)} \right]-\nonumber\\ & (1 - {z^{ - 1}}){K_2}({z^{ - 1}})v(T - 1) \end{align} $
(42) 将式(42) 代入式(38) 中得到泵池液位闭环系统方程为
$ \begin{align} \big[(1 +\, & {a_1}{z^{-1}} + {a_2}{z^{-2}} + {a_3}{z^{-3}})(1 - {z^{ - 1}})+\nonumber\\ & {z^{ - 1}}({g_{20}} + {g_{21}}{z^{ - 1}})({b_0} + {b_1}{z^{ - 1}}) \big]{y_2}(T + 1)=\nonumber\\ & ({g_{20}} + {g_{21}}{z^{ - 1}})({b_0} + {b_1}{z^{ - 1}}){y_{2sp}}(T) +\nonumber\\ & (1 - {z^{ - 1}})\Delta v(T) + (1 - {z^{ - 1}})\big[1-\nonumber\\ &({b_0} + {b_1}{z^{-1}}){K_2}({z^{-1}}) \big]v(T - 1) \end{align} $
(43) 选择 $g_{20}$ 和 $g_{21}$ 使下式成立
$ \begin{align} (1 +& {a_1}{z^{ - 1}} + {a_2}{z^{ - 2}} + {a_3}{z^{ - 3}})(1 - {z^{ - 1}})+\nonumber\\ & {z^{ - 1}}({g_{20}} + {g_{21}}{z^{ - 1}})({b_0} + {b_1}{z^{ - 1}}) \ne 0, \;\left| z \right| > 1 \end{align} $
(44) 由 ${B}({z^{ - 1}})$ 稳定, 根据式(43) 可得未建模动态补偿参数 ${K_2}({z^{ - 1}})$ 为 $1/(b_0+b_1{z^{ - 1}})$ , 将 ${K_2}({z^{ - 1}})$ 代入式(42) 可得:
$ \begin{align} {y_{1sp}}(T) =\, & \frac{{{b_0} - {b_1}}}{{{b_0}}}{y_{1sp}}(T - 1) + \frac{{{b_1}}}{{{b_0}}}{y_{1sp}}(T - 2)- \nonumber\\ & \frac{{1 - {z^{ - 1}}}}{{{b_0}}}v(T - 1)+ \nonumber\\ & \frac{{({b_0} + {b_1}{z^{ - 1}})({g_{20}} + {g_{21}}{z^{ - 1}})}}{{{b_0}}}\times\nonumber\\ & \left[{{y_{2sp}}(T)-{y_2}(T)} \right] \end{align} $
(45) 其中, $g_{20}$ 和 $g_{21}$ 通过凑式满足式(44).
3. 稳定性和收敛性分析
为了证明本文所提出泵池液位和给矿压力双速率串级控制方法可以使闭环系统具有稳定性和收敛性, 引入引理.
引理1.当控制器式(42) 作用于被控对象式(38) 时, 系统的输入输出方程可以由式(46) 和式(47) 表示.
$ \begin{align} &\left[{A({z^{-1}}){H_2}({z^{-1}}) + {z^{-1}}{G_2}({z^{ - 1}})B({z^{ - 1}})} \right]\times\nonumber\\ &\qquad {y_2}(T + 1)= {G_2}({z^{ - 1}})B({z^{ - 1}}){y_{2sp}}(T) +\nonumber\\ & \qquad{H_2}({z^{ - 1}})\Delta v(T) \end{align} $
(46) $ \begin{align} &\left[A({z^{-1}}){H_2}({z^{-1}}) + {z^{-1}}{G_2} ({z^{ - 1}})B({z^{ - 1}}) \right]{y_{1sp}}(T) =\nonumber\\ &\qquad A({z^{ - 1}}){G_2}({z^{ - 1}}) {y_{2sp}}(T)- \nonumber\\ &\qquad {z^{ - 1}}\left[{A({z^{-1}}){K_2}({z^{-1}}) + {G_2}({z^{-1}})} \right]v(T) \end{align} $
(47) 证明.将式(42) 代入式(38), 可以证明泵池液位输入输出方程式(46) 成立.将式(46) 代入式(38) 可以证明泵池液位输入输出方程式(47) 成立.
定理1.被控对象由式(7)~(10) 表示, 当 ${T\to\infty}$ 时, 未建模动态变化率 $\Delta v_1({\infty})$ 和 $\Delta v_2({\infty})$ 为常数.若采用控制律式(16) 和式(42), 则被控对象的闭环系统的输入输出信号一致有界(BIBO (Bounded input bounded)稳定), 即
$ \begin{align} \left| {{y_2}(T)} \right| < \infty, \;\;\;\left| {u(k)} \right| < \infty \end{align} $
(48) 并且, 被控对象的输出 ${y_2}(T)$ 和 ${y_1}(k)$ 与参考输入 ${y_{2sp}}(T)$ 和 ${y_{1sp}}(k)$ 之间的稳态误差小于预先设定值 ${\varepsilon _2}$ 和 ${\varepsilon _1}$ , 即
$ \begin{align} \mathop {\lim }\limits_{T \to \infty } \left| {{e_2}(T)} \right| = \mathop {\lim }\limits_{T \to \infty } \left| {{y_{2sp}} - {y_2}(T)} \right| \le {\varepsilon _2} \end{align} $
(49) $ \begin{align} \mathop {\lim }\limits_{k \to \infty } \left| {{e_1}(k)} \right| = \mathop {\lim }\limits_{k \to \infty } \left| {{y_{1sp}}(k) - {y_1}(k)} \right| \le {\varepsilon _1} \end{align} $
(50) 证明.因为
$ \begin{align} \left| {\Delta v(T)} \right| = \left| {v(T) - v(T - 1)} \right| \le \left| {v(T)} \right| + \left| {v(T - 1)} \right| \end{align} $
(51) 由式(51) 和 $v(T)$ 有界可得 $\Delta v(T)$ 有界, 即
$ \begin{align} \left| {\Delta v(T)} \right| \le 2W \end{align} $
(52) 由式(46) 和式(47) 可知, $y_{2sp}(T)$ 、 $v(T)$ 和 $\Delta v(T)$ 有界, 采用文献[20]类似方法可证:
$ \left| {{y}_{2}}(T) \right|\le {{c}_{1}}+{{c}_{2}}\underset{0\le \tau \le T}{\mathop{\max }}\, \left| \Delta v(\tau ) \right|\le {{c}_{1}}+{{c}_{2}}W<\infty $
(53) $ \left| {{y}_{1sp}}(T) \right|\le {{c}_{3}}+{{c}_{4}}\underset{0\le \tau \le T}{\mathop{\max }}\, \left| v(\tau ) \right|\le {{c}_{3}}+{{c}_{4}}W<\infty $
(54) 其中, $c_1$ 、 $c_2$ 、 $c_3$ 和 $c_4$ 为正常数.
由式(16) 作用于被控对象式(8) 可得给矿压力的输入输出方程, 由式(27) 和式(55) 表示.
$ \begin{align} &\left[{{A_1}({z^{-1}}){H_1}({z^{-1}}) + {z^ {-1}}{B_1}({z^{ - 1}}){G_1}({z^{ - 1}})} \right]u(k) =\nonumber\\ &\qquad {A_1}({z^{ - 1}}){G_1}({z^{ - 1}}){y_{1sp}}(k)- \nonumber\\ &\qquad {z^{ - 1}}\left[{{A_1}({z^{-1}}){K_1}({z^{-1}}) + {G_1}({z^{-1}})} \right]{v_1}(k) \end{align} $
(55) 由式(27) 和式(55) 可知, $y_{1sp}(k)$ 、 $v_1(k)$ 和 $\Delta v_1(k)$ 有界, 采用文献[20]类似方法可证:
$ \begin{align} & \left| {{y}_{1}}(k) \right|\le {{d}_{1}}+{{d}_{2}}\underset{0\le {k}'\le k}{\mathop{\max }}\, \left| \Delta {{v}_{1}}({k}') \right|\le {{d}_{1}}+ \\ & \ \ \ \ \ {{d}_{2}}M<\infty \\ \end{align} $
(56) $ \left| u(k) \right|\le {{d}_{3}}+{{d}_{4}}\underset{0\le {k}'\le k}{\mathop{\max }}\, \left| {{v}_{1}}({k}') \right|\le {{d}_{3}}+{{d}_{4}}M<\infty $
(57) 其中, $d_1$ 、 $d_2$ 、 $d_3$ 和 $d_4$ 为正常数.
由式(57) 和式(53) 可知闭环系统的输入 $u(k)$ 和输出 $y_2(T)$ 均有界.
由式(46), 当 ${T\to\infty}$ 时, 泵池液位的稳态跟踪误差 $e_2(T)$ 为
$ \begin{align} \begin{aligned} &\mathop {\lim }\limits_{T \to \infty } \left| {{e_2}(T)} \right| = \mathop {\lim }\limits_{T \to \infty } \left| {{y_{2sp}}(T) - {y_2}(T)} \right|= \\&\qquad \mathop {\lim }\limits_{T \to \infty } \left| {\frac{{{H_2}({z^{ - 1}})\Delta v(T)}}{{{G_2}({z^{ - 1}})B({z^{ - 1}})}}} \right|=\\ &\qquad \left| {\frac{1}{{{G_2}(1)B(1)}}} \right|\mathop {\lim }\limits_{T \to \infty } \left| {(1 - {z^{ - 1}})\Delta v(T)} \right|=\\ &\qquad \left| {\frac{1}{{({g_{20}} + {g_{21}})({b_0} + {b_1})}}} \right|\mathop {\lim }\limits_{T \to \infty } \left| {(1 - {z^{ - 1}})\Delta v(T)} \right| \end{aligned} \end{align} $
(58) 当 ${T\to\infty}$ 时, 由于 $\Delta v_1({\infty})$ 和 $\Delta v_2({\infty})$ 为常数, 根据式(37) 和式(38), $\Delta v({\infty})$ 为常数, 得:
$ \begin{align} \mathop {\lim }\limits_{T \to \infty } \left| {{e_2}(T)} \right| = 0 \le {\varepsilon _2} \end{align} $
(59) 由式(27), 当 ${k\to\infty}$ 时, 给矿压力的稳态跟踪误差 $e_1(k)$ 为
$ \begin{align} \begin{aligned} &\mathop {\lim }\limits_{k \to \infty } \left| {{e_1}(k)} \right| = \mathop {\lim }\limits_{k \to \infty } \left| {{y_{1sp}}(k) - {y_1}(k)} \right|=\\ &\qquad \mathop {\lim }\limits_{k \to \infty } \left| {\frac{{{H_1}({z^{ - 1}})\Delta {v_1}(k)}}{{{G_1}({z^{ - 1}}){B_1}({z^{ - 1}})}}} \right| =\\ &\qquad \left| {\frac{1}{{{G_1}(1){B_1}(1)}}} \right|\mathop {\lim }\limits_{k \to \infty } \left| {(1 - {z^{ - 1}})\Delta {v_1}(k)} \right|=\\ &\qquad \left| {\frac{1}{{{b_{10}}({g_{10}} + {g_{11}})}}} \right|\mathop {\lim }\limits_{k \to \infty } \left| {(1 - {z^{ - 1}})\Delta {v_1}(k)} \right| \end{aligned} \end{align} $
(60) 当 ${k\to\infty}$ 时, 由于 $\Delta v_1({\infty})$ 为常数, 得:
$ \begin{align} \mathop {\lim }\limits_{k \to \infty } \left| {{e_1}(k)} \right| = 0 \le {\varepsilon _1} \end{align} $
(61) 4. 再磨过程半实物仿真系统实验
将本文提出的再磨过程的泵池液位和给矿压力双速率区间控制方法在再磨过程的半实物平台进行实验研究, 以验证其有效性和实用性.
4.1 半实物仿真的软硬件平台
为了说明本文提出的双速率区间控制方法的有效性, 研发了赤铁矿再磨过程的半实物仿真实验系统.硬件平台由对象计算机、监控计算机、控制器设计计算机、控制系统和虚拟仪表及执行机构组成, 如图 4所示.
软件平台主要由Matlab、控制系统软件、OPC-Scout软件和相关系统软件组成, 主要实现系统的整体设计、可视化人机界面设计、被控对象模型设计、控制器设计和系统之间的实时通讯等功能, 如图 5~7所示.
4.2 被控对象仿真模型
采用欧拉法将式(1) 和式(2) 进行离散化为
$ \begin{align} {y_{\rm{2}}}(T + 1){\rm{ = }}{y_2}(T) + \frac{{{Q_{in}}(T)}}{A} - \frac{{{K_d}}}{{A\sqrt {{K_p}} }}\sqrt {{y_1}(T)} \end{align} $
(62) $ \begin{align} {y_1}(k + 1) = {y_1}(k) - \frac{1}{\tau }{y_1}(k) + \frac{{{K_p}}}{\tau }{u^2}(k) \end{align} $
(63) 式中, $Q_{in}(T)=q_1(T)+q_2(T)+q_3(T)+q_4(T)$ 为未知非线性函数, ${K_d}{\rm{/}}\left( {A\sqrt {{K_p}} } \right)$ 、 $1/\tau $ 、 ${K_p}/\tau $ 为常数.
由于式(62) 中的 ${Q_{in}}(T)/A$ 和 ${K_d}{\rm{/}}\left( {A\sqrt {{K_p}} } \right)$ 未知, 采用未知常数 $F$ 和 $K_0$ 代替, 得到式(62) 的近似模型:
$ \begin{align} {\bar y_{\rm{2}}}(T + 1) = {y_2}(T) + F - {K_0}\sqrt {{y_1}(T)} \end{align} $
(64) 未建模动态 $\Delta {y_2}(T)$ 可以表示为
$ \begin{align} \Delta {y_2}(T) = {y_{\rm{2}}}(T + 1) - {\bar y_{\rm{2}}}(T + 1) \end{align} $
(65) 由式(64) 和式(65) 得到泵池液位辨识方程为
$ \begin{align} {\tilde y_2}(T) = F - {K_0}{\tilde y_1}(T - 1) + \Delta {y_2}(T - 1) \end{align} $
(66) 式中, ${\tilde y_2}(T) = {y_2}(T{\rm{ + 1}}) - {y_2}(T)$ , ${\tilde y_1}(T) = \sqrt {{y_1}(T)} $ .
采用工业过程实际的给矿压力和泵池液位数据, 采用文献[21]的交替辨识算法获得式(66) 中 $F$ 和 $K_0$ 的估计值为 $\hat F$ 和 ${\hat K_0}$ 、未建模动态 $\Delta {y_2}(T)$ 的估计值为 $\Delta {\hat y_2}(T)$ , 泵池液位仿真模型为
$ \begin{align} {y_2}(T + 1) = {y_2}(T) + \hat F - {\hat K_0}\sqrt {{y_1}(T)} + \Delta {\hat y_2}(T) \end{align} $
(67) 采用工业过程实际的矿浆泵频率和给矿压力数据, 采用最小二乘辨识方法估计式(63) 中模型参数, 可得给矿压力仿真模型为
$ \begin{align} {y_1}(k + 1) = {y_1}(k) - {\hat K_1}{y_1}(k) + {\hat K_2}{u^2}(k) \end{align} $
(68) 其中, ${\hat K_1}$ 和 ${\hat K_2}$ 分别为 $1/\tau $ 和 ${K_p}/\tau $ 的估计值.
将式(67) 和式(68) 作为半实物仿真系统的被控对象模型.
4.3 控制目标及控制器参数选择
工艺要求在矿浆泵频率约束 $30\le u(k)\le 45$ Hz前提下, 泵池液位的控制目标范围为
$ \begin{align} 0.3 \le {y_2}(T) \le 1.3 \end{align} $
(69) 给矿压力工艺要求的最大波动范围为 $ \pm 6$ kPa, 可以表示为
$ \begin{align} \left| {\Delta {y_1}(k)} \right| < 6 \end{align} $
(70) 根据泵池区间液位控制目标范围, 确定泵池液位的参考值为 $y_{2ref}=(1.3+0.3)/2=0.8$ m.
泵池液位和给矿压力控制器设计模型参数为
$ \begin{align} \left\{ \begin{array}{l} {A_2}({z^{ - 1}}) = 1 - {z^{ - 1}}\\ {B_2}({z^{ - 1}}) = - 0.0263\\ {A_1}({z^{ - 1}}) = 1 - 0.3697{z^{ - 1}}\\ {B_1}({z^{ - 1}}) = 6.4285 \end{array} \right. \end{align} $
(71) 由式(25) 和式(26) 确定基于未建模动态补偿的给矿压力回路控制的参数为
$ \begin{align} \left\{ \begin{array}{l} {G_1}({z^{ - 1}}) = {\rm{0}}{\rm{.1661}} - {\rm{ 0}}{\rm{.0106}}{z^{ - 1}}\\ {K_1}({z^{ - 1}}) = 0.157 \end{array} \right. \end{align} $
(72) 由式(38) 确定基于提升技术给矿压力内环闭环控制系统的动态特性引入到泵池液位控制器设计模型参数为
$ \begin{align} \left\{\begin{array}{l} A({z^{ - 1}}) = 1 - 0.49{z^{ - 1}} - 0.42{z^{ - 2}} - 0.09{z^{ - 3}}\\ B({z^{ - 1}}) = - 0.0347 - 0.0074{z^{ - 1}} \end{array} \right. \end{align} $
(73) 由式(44) 和式(45) 确定基于未建模动态补偿的泵池液位回路控制的参数为
$ \begin{align} \left\{ \begin{array}{l} {G_2}({z^{ - 1}}) = - 14.20 + 10.20{z^{ - 1}}\\ {K_2}({z^{ - 1}}) = \dfrac{1}{- 0.0347 - 0.0074{z^{ - 1}}} \end{array} \right. \end{align} $
(74) 泵池液位区间控制和给矿压力回路控制组成的模糊切换控制方法的控制参数与文献[12]中的参数相同.
串级PI控制参数与本文的PI控制参数相同.
4.4 实验研究
为了验证本文提出的控制算法的动态性能, 在相同的频繁干扰下, 将本文所提控制方法、文献[12]所提控制方法以及串级PI控制方法进行半实物仿真对比实验.在21:31:50~21:36:50时间内, $Q_{in}(t)$ 在[227m $^3$ /h, 253m $^3$ /h]内频繁波动.在21:36:50~21:41:50时间内, $Q_{in}(t)$ 在[117m $^3$ /h, 263m $^3$ /h]内频繁波动, 波动曲线如图 8所示.
采用被控对象式(67) 和式(68) 作为半实物仿真对象, 对比仿真实验结果如图 9~12所示.
从图 9~12中分析可知, 在21:31:50~21:41:50时间内, 采用本文提出的控制方法可以将泵池液位和给矿压力的波动控制在式(69) 和式(70) 的目标范围内.但是明显可以看出采用文献[12]提出的控制方法和串级PI控制方法在21:31:50~21:36:50时间内, 可以将泵池液位控制在式(69) 目标范围内; 但是在21:36:50~21:41:50时间内, 使泵池液位超出式(69) 的目标范围, 同时也造成旋流器的给矿压力波动超出式(70) 的目标范围.以下是将上述三种控制器进行性能比较, 如表 2和表 3所示.
由上表可以分析出, 在21:31:50~21:41:50时间内, 采用本文控制方法对式(67) 和式(68) 的被控对象进行控制时, 使泵池液位 $y_2$ 和给矿压力变化率 $\Delta y_1$ 控制在工艺要求目标区间内.在21:31:50~21:41:50时间内, 采用文献[12]控制方法对式(67) 和式(68) 的被控对象进行控制时, 泵池液位 $y_2$ 超过被控对象的控制目标式(69) 的区间最大值为0.42, 超过该目标区间的绝对累积和为15.62;旋流器的给矿压力变化率 $\Delta y_1$ 超过被控对象的控制目标式(70) 的区间最大值为5.73, 超过该目标区间的绝对累积和为165.58.在21:31:50~21:41:50时间内, 采用串级PI控制方法对式(67) 和式(68) 的被控对象进行控制时, 泵池液位 $y_2$ 超过被控对象的控制目标式(69) 的区间最大值为0.50, 超过该目标区间的绝对累积和为22.92;旋流器的给矿压力变化率 $\Delta y_1$ 超过被控对象的控制目标式(70) 的区间最大值为6.15, 超过该目标区间的绝对累积和为84.88.半实物仿真结果表明, 本文方法在受到大范围随机干扰时, 不需要文献[12]基于规则的切换机制, 可将泵池液位 $y_2$ 和给矿压力变化率 $\Delta y_1$ 控制在工艺规定的范围内.
5. 结论
本文提出的再磨过程双速率区间控制由内环未建模动态补偿给矿压力PI控制器和基于给矿压力内环闭环控制系统动态特性的泵池液位外环动态模型设计的未建模动态补偿PI控制器组成, 采用一步最优前馈补偿律设计补偿器和PI控制器的参数.稳定性与收敛性分析和半实物仿真实验结果表明, 再磨过程受到频繁随机干扰时, 所提出的控制方法可将泵池液位和给矿压力变化率控制在目标值范围内.本文所提的双速率区间控制器设计方法对受到大范围随机干扰的非线性串级工业过程的控制器设计具有参考价值.
-
表 1 模型参数表
Table 1 Model parameters
变量 描述 $q_1(t)$ 一段磨矿矿浆和磁选矿浆流量 ${q_2}({y_1}(t),{D_H},{f_d},{\phi _H}) $ 再磨排矿流量 $q_3(t)$ 污水与冲洗水流量 $q_4(t)$ 泵池补加水流量 $D_H(t)$ 旋流器给矿浓度 $f_d(t)$ 旋流器给矿粒度分布 $\phi_H$ 旋流器结构参数 $A$ 泵池横截面积 $K_p$ 矿浆泵比例系数 $K_d$ 矿浆泵流量比例系数 $\tau$ 惯性时间常数 -
[1] Li X, McKee D J, Horberry T, Powell M S. The control room operator: the forgotten element in mineral process control. Minerals Engineering, 2011, 24(8): 894-902 doi: 10.1016/j.mineng.2011.04.001 [2] van Vuuren M J J, Aldrich C, Auret L. Detecting changes in the operational states of hydrocyclones. Minerals Engineering, 2011, 24(14): 1532-1544 doi: 10.1016/j.mineng.2011.08.002 [3] Wei D H, Craig I K. Grinding mill circuits——a survey of control and economic concerns. International Journal of Mineral Processing, 2009, 90(1-4): 55-56 [4] Duarte M, Sepúlveda F, Castillo A, Contreras A, Lazcano V, Giménez P, Castelli L. A comparative experimental study of five multivariable control strategies applied to a grinding plant. Powder Technology, 1999, 104(1): 1-28 doi: 10.1016/S0032-5910(98)00210-1 [5] 楚云飞, 徐文立, 王峻, 万维汉.基于切换控制的均匀液位控制.清华大学学报(自然科学版), 2005, 45(1): 107-110 http://www.cnki.com.cn/Article/CJFDTOTAL-QHXB200501028.htmChu Yun-Fei, Xu Wen-Li, Wang Jun, Wan Wei-Han. Averaging level control based on switching control. Journal of Tsinghua University (Science and Technology), 2005, 45(1): 107-110 http://www.cnki.com.cn/Article/CJFDTOTAL-QHXB200501028.htm [6] Pomerleau A, Hodouin D, Desbiens A, Gagnon É. A survey of grinding circuit control methods: from decentralized PID controllers to multivariable predictive controllers. Powder Technology, 2000, 108(2-3): 103-115 doi: 10.1016/S0032-5910(99)00207-7 [7] 梁蕾, 李振国.选矿过程矿浆液位的模糊控制算法研究.金属矿山, 2009, 39(7): 103-105, 135 http://www.cnki.com.cn/Article/CJFDTOTAL-JSKS200907034.htmLiang Lei, Li Zhen-Guo. Fuzzy-intelligent control arithmetic for the ore pulp level in mineral separation process. Metal Mine, 2009, 39(7): 103-105, 135 http://www.cnki.com.cn/Article/CJFDTOTAL-JSKS200907034.htm [8] Aguila-Camacho N, Le Roux J D, Duarte-Mermoud M A, OrchardM E. Control of a grinding mill circuit using fractional order controllers. Journal of Process Control, 2017, 53: 80-94 doi: 10.1016/j.jprocont.2017.02.012 [9] le Roux J D, Padhi R, Craig I K. Optimal control of grinding mill circuit using model predictive static programming: a new nonlinear MPC paradigm. Journal of Process Control, 2014, 24(12): 29-40 doi: 10.1016/j.jprocont.2014.10.007 [10] Coetzee L C, Craig I K, Kerrigan E C. Robust nonlinear model predictive control of a run-of-mine ore milling circuit. IEEE Transactions on Control Systems Technology, 2010, 18(1): 222-229 doi: 10.1109/TCST.2009.2014641 [11] Matthews B, Craig I K. Demand side management of a run-of-mine ore milling circuit. Control Engineering Practice, 2013, 21(6): 759-768 doi: 10.1016/j.conengprac.2013.02.005 [12] 赵大勇, 柴天佑.再磨过程泵池液位区间与给矿压力模糊切换控制.自动化学报, 2013, 39(5): 556-564 http://www.aas.net.cn/CN/abstract/abstract17816.shtmlZhao Da-Yong, Chai Tian-You. Fuzzy switching control for sump level interval and hydrocyclone pressure in regrinding process. Acta Automatica Sinica, 2013, 39(5): 556-564 http://www.aas.net.cn/CN/abstract/abstract17816.shtml [13] Zhao D Y, Chai T Y, Wang H, Fu J. Hybrid intelligent control for regrinding process in hematite beneficiation. Control Engineering Practice, 2014, 22: 217-230 doi: 10.1016/j.conengprac.2013.02.015 [14] Sanchis R, Romero J A, Martín J M. A new approach to averaging level control. Control Engineering Practice, 2011, 19(9): 1037-1043 doi: 10.1016/j.conengprac.2011.04.010 [15] 王泽红, 陈晓龙, 袁致涛, 于福家, 李丽匣.选矿数学模型.北京:冶金工业出版社, 2015.Wang Ze-Hong, Chen Xiao-Long, Yuan Zhi-Tao, Yu Fu-Jia, Li Li-Xia. Mathematical Model of Mineral Processing. Beijing: Metallurgical Industry Press, 2015. [16] 吴学娟, 郎朗.模糊自适应控制在变频恒压供水系统中的应用.工业控制计算机, 2010, 23(11): 53-54 doi: 10.3969/j.issn.1001-182X.2010.11.025Wu Xue-Juan, Lang Lang. Fuzzy Adaptive control of constant pressure water supplying system with frequency conversion. Industrial Control Computer, 2010, 23(11): 53-54 doi: 10.3969/j.issn.1001-182X.2010.11.025 [17] Liu F Z, Gao H J, Qiu J B, Yin S, Fan J L, Chai T Y. Networked multirate output feedback control for setpoints compensation and its application to rougher flotation process. IEEE Transactions on Industrial Electronics, 2014, 61(1): 460-468 doi: 10.1109/TIE.2013.2240640 [18] Chai T Y, Zhang Y J, Wang H, Su C Y, Sun J. Data-based virtual unmodeled dynamics driven multivariable nonlinear adaptive switching control. IEEE Transactions on Neural Networks, 2011, 22(12): 2154-2172 doi: 10.1109/TNN.2011.2167685 [19] Jia Y, Chai T Y. A data-driven dual-rate control method for a heat exchanging process. IEEE Transactions on Industrial Electronics, 2017, 64(5): 4158-4168 doi: 10.1109/TIE.2016.2608878 [20] Chai T Y, Zhai L F, Yue H. Multiple models and neural networks based decoupling control of ball mill coal-pulverizing systems. Journal of Process Control, 2011, 21(3): 351-366 doi: 10.1016/j.jprocont.2010.11.007 [21] Zhang Y J, Chai T Y, Wang D H. An alternating identification algorithm for a class of nonlinear dynamical systems. IEEE Transactions on Neural Networks and Learning Systems, 2016, PP(99): 1-12 https://www.researchgate.net/publication/301278679_An_Alternating_Identification_Algorithm_for_a_Class_of_Nonlinear_Dynamical_Systems?_sg=0GrGvPEsUQO8xrHfwsP-YyaZoz3IXJzhTLh9wx6sBEYsfD1rteepDimRIMqZ-HKCo1BilrLdiewJkJagj94BSA 期刊类型引用(7)
1. 代伟,张政煊,杨春雨,马小平. 基于SCN数据模型的SISO非线性自适应控制. 自动化学报. 2024(10): 2002-2012 . 本站查看
2. 张威,张强,王玲晓,刘江,张俊飞. 塔磨分级过程的智能控制方法及应用. 现代矿业. 2023(04): 135-139 . 百度学术
3. 李伟涛,苟晓东. 未建模动态补偿驱动的重介质选煤自适应控制. 洁净煤技术. 2023(12): 128-135 . 百度学术
4. 曾喆昭,刘文珏. 自耦PID控制器. 自动化学报. 2021(02): 404-422 . 本站查看
5. 任双兴,侯磊,李凌峰,靳宇宏,昌泽元,姚迟森. 一种用于转子裂纹故障检测的区间控制策略. 动力学与控制学报. 2021(06): 76-82 . 百度学术
6. 姜艺,范家璐,贾瑶,柴天佑. 数据驱动的浮选过程运行反馈解耦控制方法. 自动化学报. 2019(04): 759-770 . 本站查看
7. 郑杰辉. MCU存储器双倍速率控制技术研究. 计算机测量与控制. 2018(11): 133-136 . 百度学术
其他类型引用(15)
-