S. Howes, I. Mohler, N. Bolf
PiControl Solutions, Houston, TX, USA, e-mail: firstname.lastname@example.org
Thispaperillustrates a novel, modern and practical method for tuninglevelPID controllers (LCs). The example is based on a distillation column sump.The new PID tuning algorithm works entirely in the time domain(no Laplace or Z-discrete domain). It neitherrequiresdata preprocessing nordata preconditioning. The process modelparametersarefirst identified followed by simulation and optimization of the PID tuning parameters. The process identification andtuning techniques are performedusing Pitops and Simcetsoftware from PiControl Solutions LLC.
Many level controllers (LCs) in chemical plants do not work wellbecause their PID tuning parametersare notoptimally tuned.Poor LC tuning can cause cycling or sluggish control response.Cyclingof LCoutputcan causemanydownstream flows to oscillate. This in turn can cripple the performance of higher level Advanced Process Control (APC) systems resulting in reduced profit marginsand fail to deliver the expected payback.
Many authors use different commonly known tuning techniques for controller parametersoptimization.These include: Ziegler–Nichols tuning formulaslike:Paor et al. 1989,Hang et al. 1991, Astrom et al. 2004,Hagglundet al. 2008, Lambdamethod like: Ingimundarson et al. 2005,Lennartson et al. 2009, Cohen Coon tuning method like:Chen 1989, Astrom et al. 1993, Kilian 2000and Shen 2002or Internal Model Control (IMC) tuning method like Seborg et atl. 2004.Many other tuning methods and their comparisons are available in various open literaturesbut are not the subjectof this paper.
This paper shows a new, quick, simple and powerful new PID tuningmethodwhich is simpler and more practicalthan currently known methodologies. Past efforts in this area include contributions from: Mehra et al. 1972, Gabriele et al. 1977and Perdew et al. 1996.With this new method described in this paper,it is possibleto tune LC parameters effectively and improve the overall plant control qualityquickly and easily.
2.Chemical Process Dynamics
Most chemical processdynamics can be well represented by the transfer functions shown in Figure 1.
The first and second transfer functions are typical for flow, pressure and temperature controllers. The zero order transfer function characterizes level controllers (LCs).Detailedexplanation on transfer functions isavailable in Seborg et al. 2004.Figure 2illustrateszero order processes. If flowisgoing into a tank and exact flowisleaving the tank, then the level in the tank will beheld constant. Now,if the flow out from the tank is reduced by ΔF, then the level starts to rise at a steady rate (called Ramp Rate) and will continue to rise until the tank overflows. Let’s say, after time Δt, the change in level is ΔL. The Ramp Rate is calculated as:
Ramp Rate= ΔL/ Δt /ΔF
Note that the engineering units of Ramp Rate in this example will be%/minutes/m3/h.This is the zero order transfer function behavior, where there is no final new steady state like a first, and secondor higher order transfer functions. The zero order transfer function is also called an integrating type transfer function. At a first glance, the zero-order transfer function appears the simplest since it has only two parameters:Process Dead Time and Ramp Rate. There are no Time Constants involved as in the case of first and secondorder transfer functions. However, the zero-order transfer function interestingly, poses some unique challenges when it comes to determining optimal tuning parameters forthe LC PID.
Figure 3, shows a distillation column bottoms sump LCtrends. The LC isacascade controlPID. LC is the master loop and itsslave is the bottoms product flow controller (FC) connected directly to the valve.
The top window shows the LC setpoint (SP), the horizontal line at a fixed SP of 60. The cyclical trend in the top window is the actual process value(PV). Notice that the level is cycling about ±5% from the LC’s SP. The bottom trend shows the LC’s output, OP. The OP directly manipulates the SP of the FC. If the LC’s PV increases, the LC’s OP increases the FC’s SP to take more liquid out of the sump and vice versa.
Level is the controlled variable (CV) because it is controlled by the flow. The bottoms flow is the manipulated variable (MV) because of it the level is kept on the desired value. The level (CV) range is 0 –100% and the FC (MV) range is 0 –200 t/hr.
The LC’s PID equation is:
OP(t)= OP(t-1)+ P [ΔE + E(t)Δt/ I + D Δ(ΔPV)/Δt ]
OP(t)= PID controller Outputin the time t;
P = Proportional Gain;I = Integral Constant;
D = Derivative Constant;
E(t)= Controller Error (PV(t)–SP(t))inthe timet;
ΔE = Controller Error (E(t)–E(t-1));
ΔPV = Process Value (PV(t)–PV(t-1));
Δt = Scan Time of PID.
The initial LC’s tuning parameters that are causing the oscillations in Figure 3are:
P = 0.5,
I = 3 min,
D = 0min.
To scientifically calculate optimum PID parameters for the LC, first put the LC in Manual mode and set the LC’s OP at approximately the average value to try and keep the level reasonably stable, like inFigure 4.
Now bump the flow for at least 10 minutes to see anoticeable change in the level like in Figure 5which shows a bump in the FC’s SP by 5 ton/h.
The level is allowed to drop all the way to about 30%. This is a rather large change and is shown here only for clarity of the illustration, but the level could have been allowed to drop to only 45% in as horter time period without much loss of information.
The Ramp Rate is calculated as(57 –33) /(160 –60) / 5= 0.05 %/(ton/h)/min. The process dead time can be seen by visual inspection of the data and this is approximately 2 minutes in this case. Now we can build a LC simulation in the Pitops software with the following additional configuration information:
MV range = 0 –200 ton/h(PV range of the slave FC)
CV range =0 –100%(PV range of the LC)
New tuning PID parameters are optimized by Pitops software from PiControl Solutions LLC. A simulation illustration is shown in Figure 6.The new optimized tuning PID parameterscorrespondto the mathematical optimizationsolution that minimizes the error for the following custom controller challenges:
2.Typical noise in the LC PV
3.Injection of a Ramp disturbance simulating a production rate change.
4.Injection of a Pulse disturbance simulating minor upsets.
New tuning PID parameters are:
P = 6.2,
I = 22min,
D = 0.25min.
The optimized tuning PID parameters might be a little aggressive for the real process. The optimized tuning sets upper limits on the level of aggression and tuning tightness, beyond which incipient oscillations will start andmight grow, making the controller possibly unstable.To strike a balance between tight control and smoothness of downstream flow, the final tuning PID parameters aresomewhat detuned to be:
P = 2,
D = 0min,
Gapgain = 0.5,
Gaphigh/low = 2%.
If the level is within 2% of the SP, the gap action detunes the controller action by 50% and to further smooth out the downstream flow changes.Thiscontrol action is shown in Figure 7.Now, as typical disturbance come, the control action is just about right. The levelis controlled nicely without excessive oscillations or jerking around of the downstream flow. Note that the ripples seen are because of external oscillatory disturbances causing the level to cycle a little, and this could be because ofupstream valve problems on interacting controllers or process issues which are often unavoidable.
All illustrations above were generated using Pitopsand Simcetprocess control software from PiControl Solutions Company. The use of thissoftware and the RGR techniques allow fast and precise tuning of critical controllers and can reduce oscillations and improve control performance in any industrial process. The software benefits tuning of all PID controllers encountered in the industry:FC, PC, TC, LC, AC, compressor surge control, motor control turbine control, power control, robotics and all related control problems.