WO1991010177A1

WO1991010177A1 - Method and device for the consideration of varying volume and flow in the control of a continuous flow process

Info

Publication number: WO1991010177A1
Application number: PCT/FI1990/000315
Authority: WO
Inventors: Antti Johannes Niemi
Original assignee: Antti Johannes Niemi
Priority date: 1989-12-28
Filing date: 1990-12-28
Publication date: 1991-07-11
Also published as: FI896298A0; EP0507794A1; FI88343B; US5428527A; JP3133061B2; FI88343C; EP0507794B1; FI896298A; DE69030106T2; JPH05504857A; DE69030106D1

Abstract

The control method is directed to a quality characteristic of the process material, typically to its outlet concentration. In control of the actuator which is a part of feedforward or feedback control of outlet concentration, one takes account of changes of the material volume and flow, in addition to the inlet concentration or control deviation. With this aim, at least two of the quantities inlet flow, outlet flow and material volume are measured or determined and by means of them a variable is further determined with regard to which the control method is invariant. Two alternative variables are presented, one being typically applicable to control of a process with finite distribution of the flow velocity and the other to control of a process with vigorous mixing. By means of a time-variable control method which has been derived from the method presented above, the controls can also be implemented with the time as the independent variable, and two PID control methods are presented corresponding to this. The invention relates also to the devices which implement the control methods presented.

Description

Method and device for the consideration of varying volume and flow in the control of a continuous flow process

The present invention relates to a control method which takes account of changes occurring in the material volume and mass flow of a process to be controlled, at the manipulation of an actuator, like a control valve, for control of an other process quantity. The invention also relates to apparatus for implementing the method.

In control of the concentration or other quantity which describes the instantaneous state of the process material, changes of the amount of material contained in the process and changes of the material flow make the control task difficult. For observation of variation of the continuous flow, a method and the corresponding apparatus have been disclosed earlier, and they have received a patent in Finland (59494), United Kingdom (2051424) and United States (4358821). In the referred patent documents, common disturbances and changes of flow and the benefit gained through their observation are illustrated by examples. In the method disclosed, the material volume of the process has been considered constant and thereby one has completely ignored its possible variation, the disturbances caused by this and incorporation of its variation in control of other process quantities. On the basis of the above method, one cannot either make any obvious conclusions on how changes of material volume could be taken into account in the pursuance of better control results. However, the amount of material, like the volume of a liquid, varies often in the industrial equipment. Typical examples of this are buffer vessels charged with the task of decreasing the changes which are present in both the quality and the flow of the incoming material. In such a case, both the material volume and flow, in addition to the input concentration, affect the concentrations in the output channel.

The controllers of a production process of variable volume have to be tuned for the worst case occurring in the practice. If e.g. the process exhibits a time delay which is essentially dependent on the volume, the feedback controller of concentration has to be tuned for the greatest volume present in practice. This is so, because the control would get worse and the control loop could run into a state of oscillation, if the volume could become greater than the volume under tuning conditions.

The applicant is not informed of examples of cases in which the variable volume would have been incorporated in control of material properties, the less so of cases in which both the volume and flow would have been incorporated. It may be possible to construct such a control system by means of conventional control components whereby the system would consist of e.g. cascade control loops. An inaccurate operation would be characteristic to such a controller including the fact that the control parameters would deviate almost always from their optimal values. The aim of the present invention is therefore to provide a method and apparatus for process control which takes the changes of process volume and flow into account. The invention expands the conventional, known control methods and apparatuses with this feature. Its principal characteristics appear from the patent claims enclosed. Especially an apparatus used for feedback control may typically be a proportional, integral and derivative controller whereby it incorporates an improvement with regard to the conventional PID controller in the respect that it performs its duties equally well at variation of the volume and flow as while these keep to constant values, while the operation of a loop provided with a conventional PID controller becomes worse, if they deviate from the constant values corresponding to the conditions of tuning of the controller.

Description of the method

A basis of the invention is the observation of the fact that, at a change of the liquid volume, the input and output flows or at least one of these must change. The concentration dynamics of a vessel do therefore not depend unambiguously on the volume. Corresponding to this, the control of the output concentration has to take account of at least two disturbance quantities, while the input concentration is considered the controlling quantity. This is not possible with previously known methods.

Differing from the controller according to the earlier patent which is aimed to control of a process of variable flow but constant volume, no such single variable can be pointed out at a variation of the volume which would apply to modelling of all flow characteristics and would have such a difference of its mutually independent values which could be used as the argument of the weighting function of the concentration process. At least the general nature of the flow characteristics has to be known and correspondingly in the following, the control of the vessel of finite distribution of the flow velocity will be discussed first and subsequently the control of the perfectly mixing vessel.

In many vessels, the velocity of motion of liquid elements is distributed to both sides of a finite average value. It has been shown (A. J. Niemi, Mathl Comput. Modelling 11 (1988) pp. 32-37) that, if the flow pattern keeps unchanged, the average axial, relative displacement x_av/L of the liquid elements, i.e. the average displacement x_av proportional to the axial length L of the vessel, in the interval ν > t, depends on the average velocity of motion and therefore on the gross flow Q and volume V. The following result is then valid, irrespective if the variation of the volume V(t) proceeds in the axial direction, perpendicularly to it or simultaneously in all directions. μ

( 1 )

It has been shown in the above reference that, if the velocity of motion is distributed, the corresponding relative distribution x/x_av of the location x depends at each time t on the quantity x_av/L. It has further been shown that the residence of liquid in a vessel of length L between the time of entrance (V) and time of exit (t) depends on the same quantity (1) which will be symbolized by z in the following.

It has been shown that many industrial continuous flow processes can be described by means of residence distributions p(z) which are based on this variable whereby the explicit time variable is eliminated from said descriptions. Examples of processes of this kind are the plug flow and laminar flow and, further on, the Gaussian velocity profile and, on certain conditions, the dispersion flow vessel. The description can be generalized to such vessels of arbitrary flow pattern in which the relative velocity profile may be different in different cross sections, especially in the neighbourhood of the vessel ends, but in which it anyway is continuously invariable in each individual cross section. Such a process needs not to have a closed mathematical model, while its flow pattern, i.e. typically its residence distribution, may be known e.g. numerically, through experiments and measurements.

The residence time distribution of the process, and its residence distribution presented in terms of z are, assuming that the other general conditions are satisfied, i.e. the flow pattern is unchanged and the variables and parameters which have not been singled out here are constant, also weighting functions of the concentration process. With their aid, the outlet concentration of the process can therefore in principle be computed assuming that the said distributions are known on the basis of theoretical or logical inference or of accomplished experiments or measurements. However, one cannot directly perceive, how such a computation is performed. Two of the three variables, i.e. of the volume, inlet flow and outlet flow, may change freely, while the third one is determined by them. In the invented method, one requires that two of the said quantities are measured continuously or repeatedly. If the measured variables are the volume and inlet flow, the value of the variable z can be determined for each pair (t, V), and further on, the value of the outlet concentration C at the time t(z) on the basis of the inlet concentration C_in (i=1):

z ( ) (z )

In the equation (2), z has been presented as the difference of its two factors:

1₁ )

η= fixed origin of time

Both (2) and (3) are applicable to the computation of the outlet concentration. It is essential, that the storing of a single function p only which is a function of one variable, in the memory is sufficient, in addition to that of the measured values of C_in and Q_in. Especially if the model (2) is used, it is practical to perform the integration with regard to ζ₁ with constant increments, in the decreasing direction. The corresponding increment of time is thereby obtained from the equation (5) and further on, v by means of these increments, while the corresponding C_in ( ν) [=C_in (ζ₁) which have not been distinguished above from each other] can be read or interpolated from the memory. The integration is finished as such low values of the time are reached which do not any more essentially influence the result. ( 5 )

The above discussion of the process model was necessary, because the feedforward control method according to the invention exploits it. It is known from previous art that the effect of one disturbance variable can be eliminated completely by feedforward means, if the process model is accurately known. A method has been presented in the referenced patent for simultaneous elimination of another disturbance quantity. The previous art does not know a method which, in addition to the main disturbance quantity, would eliminate the effects of the two other disturbance quantities in the process, at least if this is to be made accurately.

Fig. 1 presents the feedforward control loop of a continuous flow process in which the present value and earlier values ζ₁ of the variable z_i are determined. In the case of a process with a finite distribution of the flow velo&ity i=1, and in the case of a well mixing process i=2. The controller determines the controlling quantity u as a function of this variable on the basis of the measured outlet concentration and directs this control to the process. Fig. 2 presents the feedback control loop of a continuous flow process. The variable z_i (ζ₁) is determined equally to the case of Fig. 1. The controller determines the controlling quantity u on the basis of the difference of the prescribed value and measured value of the outlet concentration.

The feedforward control method according to the invention presented in Fig. 1 thus requires that the main disturbance quantity, i.e. typically the variation of the inlet concentration, is measured by means of a suitable measuring device or analyzer and stored in the control device, i.e. in the control computer. Similarly it requires that the inlet flow and liquid volume of the vessel are measured and stored. On the basis of the these ones, the values of z₁ and ζ₁ are then computed, and further on that of C(t), as presented above. The correction corresponding to the value -C(t) is directed to the input or output point of the process or to a point between them, by means of a proportioning pump or other suitable actuator. With regard to the last alternative above, a reference is made to an earlier publication (A. J. Niemi, Proc. ISA Conf. 1978, Part II, pp. 63-68 & Proc. JACC 1978 Part I, pp. 37-42. ISA, Pittsburgh) in which a method has been presented for control of a process of constant parameters. This method can be logically augmented with the features according to the invention which take the variable inlet flow and volume into account.

The two quantities to be measured can be chosen also differently from the set of the volume and inlet and outlet flows, because these are connected by the following relationship:

( 6 )

= Q_in (t) - Q (t) V (η) = V_η

If the flow Q_in or volume V which can be obtained also indirectly from the level height of liquid are not measured, it can always be calculated by means of the two other, measurable quantities from the equation (6) or from its integrated form. - If the inlet and outlet flows are measured continuously or repeatedly, the volume must anyway be also measurable, at least at longer intervals. Otherwise the systematic errors possibly present in the flow measurements may accumulate and gradually result in an essentially incorrect value of the volume. The updating of the initial value of the volume in the equation (6) corresponds to the above.

The invention relates also to the feedback control method. In this method, one controls, on the basis of the measured value of the output quantity which is typically the outlet concentration, the input quantity which is typically the input concentration, in such a manner that the previous one follows as accurately and stably as possible the prescribed value which can be chosen and adjusted by the user. This controller does not usually contain a process model according to the function p. Despite of this, it can be expediently constructed for execution of the control actions as functions of z , without an explicit dependence on the time t.

In this method according to Fig. 2, one measures at least two of the quantities inlet flow, outlet flow and volume, simultaneously as they appear or are directed to the process. The measured signals are brought into the controller and the dependence (6) is used for calculation of Q_in or V, as may be needed. The controller first calculates the needed values of z₁ (and ζ₁). By means of them and the deviation e which is the difference of the prescribed and outlet concentrations, it further on determines the needed control u and directs it to the process through a suitable actuator. The controller may have the structure of any known feedback controller assuming this is able to operate with z as the argument instead of the time t. It may then be useful to determine the control with regard to z₁ with equal increments Δζ₁, and to accomplish the needed integrations in the direction of the decreasing value of the variable. Thus e.g. the algorithm according to the general principle of the proportional, integral and derivative control or PID control determines the controlling variable as follows, with i=1: z z K_I _D

_i

As now a process with constant parameters is controlled with a controller of constant coefficients, the closed control loop is invariant with regard to z. At the same time, it is linear due to the characters of the process and controller. For the design and tuning of the controller and i.a. for selection of values of its parameters K_p, K_I and K_D, one may therefore use the design methods of linear constant coefficient systems which have been described in textbooks on control engineering, which are in general use and which have been described i.a. in the referred earlier patent. The controller can be tuned also by experimentation, e.g. whenever the process model is not sufficiently accurately known. It is then useful although not necessary that the volume and flows are unchanged at the time of tuning of the controller. Whenever they afterwards change, the controller takes the changes into account, and the result of control is equally good as at a time they keep to constant values.

The method can be implemented also with the time as the independent variable. Especially in the case of PID control, one obtains, by means of the control function (7) and equation (4) (i=1), the control function:

+

n

From the control functions (7) and (8) one observes that if the flow Q_in and volume V do not change but keep continuously to their nominal values, the controller in fact perfoms the integration and differentiation with regard to the normal time variable and thus operates as controller of other process quantities in the fully same manner as the conventional PID controller. One observes further on, that the K-coefficients of the control functions can be calculated in a simple manner from the corresponding coefficients of the conventional control or from the coefficients which possibly have been used for control already earlier in which case they need not any more be retuned. In a perfectly mixing vessel, the entering material is immediately mixed within the total material volume of the vessel. It has no finite velocity distribution and no such average value of the speed of propagation which could be defined. Instead of that, it can be described by an ordinary differential equation, also under variation of its parameters.

= Q_in(t)C_in(t)-Q(t)C(t) ( 9 )

The solution of the equation can be presented in the form (2), while i=2 and p has the meaning according to (10) and z ' that according to (11) which differs from that characterizing the process type (4) discussed earlier.

p(z')dz' = e^-z' dz'

( 10 )

The presented residence time distribution is simultaneously the weighting function of the concentration process. With its aid, one may determine the outlet concentration of the perfect mixer. In the method according to the invention, one requires that at least two of the quantities inlet concentration, outlet concentration and volume are measured continuously or repeatedly; the missing one can be determined by means of equation (6) as needed. The output quantity can be determined with either the time or the quantity z ' (11) as the independent variable. Especially if the model (2) is used, it is useful to perform the integration with regard to z₂ with equal intervals.

In the feedforward control (Fig.1) of also the ideal mixer or of a vessel with closely similar characteristics, one requires that the main disturbance quantity, typically the variation of the inlet concentration, is measured and stored in the control computer. Similarly one requires that the inlet flow and liquid volume are known through either mesurements or calculations. The output quantity C(t) is computed as presented above and the correction corresponding to its value is directed through an actuator to the inlet or outlet point of the process or to a point between them. - The same additional comments which were above presented for the case of processes which had an invariable relative velocity distribution, apply also to the feedforward control of the perfectly mixing vessel.

The invention relates also to the feedback control method of the outlet concentration of the ideal mixer. Also this kind of controller which usually does not contain a process model p can be constructed to perform its control actions as functions of z ' (11), without an explicit dependence on the time t. Also in this method according to Fig.2 one measures at least two of the variables inlet flow, outlet flow and volume and brings the measured singnals into the controller. After this has determined z₂, it determines the needed control u by means of z₂ and difference quantity e, and directs it to the process. The controller may have a structure of any such known feedback controller which is able to operate with z ' (11) as the argument. instead of t. E.g. the PID control or PID control algorithm determines the controlling quantity according to equation (7) while i=2.

The feedback control of the ideal mixer is generally subject to the same comments which were above earlier presented in connection with feedback control of a vessel with invariable relative velocity distribution. Also this controller can be constructed to accomplish its control actions as functions of the time whereby the PID algorithm obtains the following form on the basis of (11): I

The presented control methods eliminate the effects of changes of volume and flow completely, in control of the content of material or its other quality characteristic, if the process model is accurately known. With variation of volume, one can match the variation of the pressure in a gas tank in which the concentration of its outlet flow is being controlled. Hereby it may be beneficial to use mass flows instead of volumetric flows. With it one can match also the variation of a fraction of solid material and thus of the mass of solid material in the solid/liquid suspension in a vessel of constant volume whereby the quantity to be controlled is a quality characteristic of solid material in the outlet flow.

If also other physical processes take place in the control object, in addition to variation of the concentration due to effects of volume, flow and mixing, the presented controllers eliminate only the effects of these stated variables. Such other process of the subject type may be e.g. a thermal process or chemical reaction which proceeds depending on kinetic factors and is typically bound with the absolute time variable. Also in this kind of a case, the presented controllers bring an essential improvement with regard to the conventional control, since they eliminate the effects of the volume and flow. If the control is performed with regard to z, one has in any case often to redetermine values of the time ν in which case the determination of the kinetic member for different increments of integration does not any more reguire an essential additional work.

The presented control methods can be implemented in a straightforward manner, if a computer is used as the controller. The presented computations are easily programmed, including the integration of primary measured variables and integration and differentiation of the control deviation or of other input quantities of the controller. The needed measuring devices and likewise the converters which bring the measured signals to digital form are known from previous art. Also a high number of actuators applicable to control of feed of chemicals is previously known, as well as digital-to-analog converters which transmit their control signals. The control computer which operates in real time may naturally serve also such additional control and other tasks which are independent of the controls presented here.

Claims

1. A method for incorporation of varying volume and flow in the control of other process quantities characterized in that, for determination of inlet flow (Q_in) and volume (V), at least two of the process variables inlet flow (Q_in), outlet flow (Q) and volume

(V) are measured, the inlet flow (Q_in) either is integrated with regard to time and then divided by the volume (V) in order to determine a new variable (z) or is divided by the volume (V) and the quotient integrated with regard to time in order to determine a new variable (z'), and the output quantity (u) of the controller is determined with a method which is invariant with regard to this variable (z or z ' ) .

2. An apparatus according to Claim 1 for implementation of feedback control in which the difference (e) between the prescribed value (C_ref) set to the output quantity and the measured value (C_o) of the output quantity of the process is determined, characterized in that, in order to determine the output quantity (u) of the control unit,

a) said difference is multiplied with a first coefficient,

b) said difference is multiplied with the value of the inlet flow, the product is integrated with regard to time, the result of integration is divided with the value of the volume and the result is multiplied with a second coefficient, c) said difference is differentiated with regard to time, the derivative is multiplied with the value of the volume, the product is divided with the value of the inlet flow and the quotient is multiplied with a third of coefficient, and

d) the results obtained either in a) and b) or in a) and c) or in a), b) and c) are summed up, the coefficients stated in a), b) and c) being constants or slowly changing quantities.

3. An apparatus according to Claim 1 for implementation of feedback control in which the difference (e) between the prescribed value (C_ref) set to the output quantity and the measured value (C_o) of the output quantity of the process is determined, characterized in that, in order to determine the output quantity (u) of the control unit,

a) said difference is multiplied with a first coefficient,

b) said difference is multiplied with the value of the inlet flow, the product is divided with the value of the volume, the quotient is integrated with regard to time and the result of integration is multiplied with a second coefficient,

c) said difference is differentiated with regard to time, the derivative is multiplied with the value of the volume, the product is divided with the value of the inlet flow and the quotient is multiplied with a third coefficient,

4. A method for incorporation of the varying mass and flow in the control of other process quantities, characterized in that, in order to determine the inlet flow and mass of gas or a quantity proportional to the latter, at least two of the process variables inlet flow, outlet flow and mass are measured, the inlet flow either is integrated with regard to time and then divided with a quantity which is proportional to the mass in order to determine a new variable or divided with a quantity which is proportional to the mass and the quotient integrated with regard to time in order to determine a new variable, and the output quantity of the controller is determined with a method wich is invariant with regard to this variable.

5. An apparatus according to Claim 4 for implementation of feedback control in which the difference between the prescribed value set to the output quantity and the measured value of the output quantity of the process is determined, characterized in that, in order to determine the output quantity of the control unit,

a) said difference is multiplied with a first coefficient,

b) said difference is multiplied with the value of the inlet flow, the product is integrated with regard to time, the result of integration is divided with the value of the pressure and the result is multiplied with a second coefficient,

c) said difference is differentiated with regard to time, the derivative is multiplied with the value of the pressure, the product is divided with the value of the inlet flow and the quotient is multiplied with a third coefficient, and

6. An apparatus according to Claim 4 for implementation of feedback control in which the difference (e) between the prescribed value set to the output quantity and the measured value of the output quantity of the process is determined, characterized in that, in order to determine the output quantity (u) of the control unit,

a) said difference is multiplied with a first coefficient,

b) said difference is multiplied with the value of the inlet flow, the product is divided with the value of the pressure, the quotient is integrated with regard to time and the result of integration is multiplied with a second coefficient,

c) said difference is differentiated with regard to time, the derivative is multiplied with the value of the pressure, the product is divided with the value of the inlet flow and the quotient is multiplied with a third coefficient,