The disclosure relates to a method and an apparatus for control of a fluid conveying system, comprising at least one pump, at least one consumer as well as at least one armature as an actuator, wherein pressure and volume flow rate of the consumer are controllable.
The closed-loop-control of the process variables of the volume flow rate and thus also of the pressure is the standard task of final controlling devices in technical process installations. Valves or armatures with for example electric or pneumatic drives, also known as control valve or armature with actuator, are preferably used as final controlling devices. Their adjustable flow resistances affect the volume flow rate and the pressure within the installation.
Besides the valves, pumps are the most important components of an installation, as they are causing the movement of fluid. Among the wide range of possible pump designs, the centrifugal pump with drive, in most cases an electric motor with a frequency converter, is the standard solution in many areas of application. Closed-loop-control of process variables by means of a pump can be achieved via the rotation speed of the pump. Just like the stroke or the valve/armature position in case of a valve, in case of a pump the volume flow rate and thus the pressure is affected by changing the rotation speed. Even though in new technical process installations today the portion of speed controlled drives amounts to about 20% to 25%, these are rarely integrated actively into the process control but are rather employed for stationary correction of the pump operating point.
A large number of applications include closed-loop-control tasks requiring for example a large adjustment range owing to the process variables. This task cannot be realized through closed-loop-control by means of the pump only on the one hand and through closed-loop-control with an armature as actuator only on the other hand. The combined use of pump and valve with associated controller opens up new possibilities in process design. However, by combining the devices the controller design becomes more complex since a multi-variable system with 2 inputs must be dealt with. Besides the coupling of the process variables, dead-times often occur in technical process installations, which additionally complicates the closed-loop-control task.
It is an object to develop a closed-loop-control concept for an installation that allows to independently act at a consumer upon the two process variables present in the installation, that is pressure and volume flow rate. The installation provides an arrangement of at least one pump, at least one consumer and at least one armature as actuator.
In a method or apparatus for closed-loop-control of a fluid conveying system, at least one pump, at least one consumer, at least one control valve, and at least one armature as an actuator of said control valve are provided. Pressure and volume flow rate of the consumer are controlled independently of each other by means of a decoupling controller.
For the purposes of promoting an understanding of the principles of the invention, reference will now be made to the preferred exemplary embodiment/best mode illustrated in the drawings and specific language will be used to describe the same. It will nevertheless be understood that no limitation of the scope of the invention is thereby intended, and such alterations and further modifications in the illustrated embodiment and such further applications of the principles of the invention as illustrated as would normally occur to one skilled in the art to which the invention relates are included.
The control target achieved in the controller design enables a high control dynamic, notwithstanding the occurring dead-time, an unprecedentedly low oscillation tendency of the controlled variables, as well as a high stationary control accuracy. Additionally, the controlled variables pressure and volume flow rate at a consumer are decoupled.
When designing the model it was assumed that a flow is balanced at an entry e and an exit a of a fluid conveying system, i.e. no significant changes of the flow variables occur across the cross-sectional area of a connecting duct system. This applies when in the areas concerned a duct cross section is constant and a frictional effect of duct walls is low. In the area between inlet and outlet a three-dimensional flow profile is allowed for. Under these assumptions the flow can be described by the flow variables along a mean flow line according to flow line theory. In this context the Bernoulli equation
applies, with:
g=gravitational acceleration z=height coordinate
c=flow speed ρ=fluid density
s=trajectory coordinate p=pressure
The flow speed is described by the quotient of volume flow rate Q over the cross-sectional area A through which the flow passes. If the fluid is incompressible and the flow cross section is constant, then c=f(t) applies, and the solution to the integral can be expressed in general terms:
L denominates the conduit length along the flow line between the points e and a.
Since the Bernoulli equation only describes flow without friction, frictional effects are described phenomenologically within flow line theory by means of pressure sinks Δpc at the valve and Δpc at the consumer. Analogously a pressure increase within the pump Δpp is described as pressure source. Thus, the following relationship results for the Bernoulli-equation between the points e and a:
The model of a test bench used for the experiments assumes for the following simplifications:
With these assumptions it follows for the fluid dynamic
Besides the fluid dynamic the dynamic of the final controlling devices, i.e. of pump and valve play a decisive roll for the total dynamic of the installation. The used final controlling devices dispose of subordinate controls that adjust in stationary fashion the control variables stroke H of the valve and rotational speed n of the pump exactly to the predetermined setpoint values. The dynamic behavior of the final controlling devices is modelled as a first order delay element, respectively. It follows for the dynamic of the final controlling devices
The control variables H and n influence the relationship between the process variables.
The relationship between differential pressure Δpv and volume flow rate Q of the valve is dependent on the medium used, the flow coefficient Kv,v and on installation and flow conditions. Water is used as a fluid, and a turbulent flow as well as standardized installation conditions without fittings are assumed. Then the relationship
applies, wherein Δpref denominates the reference pressure of one Bar. The dependency of the flow coefficient of the valve Kv,v on the valve stroke H is described via a non-linear valve characteristic. The relationship
known from equation (7), also applies in relation to the differential pressure of a consumer Δpc and the flow Q, wherein, however, the Kv,c-value of a consumer is usually constant. The relationship between the process variables rotational speed n, flow Q and differential pressure Δpp of a pump is described by the so-called throttling curve
Δpp(t)=hnnn2(t)+hnvn(t)Q(t)+hvvQ2(t). (9)
The derivation of equation (9) was given in the standard work of C. Pfleiderer, Die Kreiselpumpen, Springer-Verlag, 4th edition, 1955, analogously. The denomination of the variables is modelled on a publication by R. Isermann, Mechatronische Systeme, Springer-Verlag, 4th edition, 2008.
From the equations (4), (5) and (6) the state space model of the installation can be derived with the states Q(t), H(t) and n(t) and the inputs Hsoll(t) and nsoll(t). For this purpose the pressures in (4) are eliminated by substituting (7), (8) and (9). According to the definition of the control targets,
y1(t)=Q(t) (10)
y2(t)=p1(t)=pe(t)+Δpp(t) (11)
applies for the output values y.
In summary it follows the state space model (12). For reasons of clarity, in the following the time dependence of the variables is omitted.
{dot over (x)}=f(x)+G·n
y=h(x) (12)
with
x=(Q H n)T
In practice dead-times occur owing to the signal processing in the final controlling devices and the measurement instruments. These dead-times cannot be neglected, compared to the process dynamic. Likewise, if an inductive flow sensor is used, dead-times must be accounted for in the model. It was realized that an improvement of the control can be achieved by means of a decoupling of the controlled variables.
A control path is assumed to be coupled when an output yi is controlled by several manipulated variables uj. Ideally, the decoupling controller is comprised of a coupling between the setpoint variables wi and the manipulated variables uj, the coupling being inverse to the control path. In the controlled entire system each output yi is thus dependant on only one setpoint variable wi, owing to the effect of the controller.
The couplings in the system cannot be derived directly from the state space model (12) since the manipulated variable u generally only acts upon one of the higher derivatives of the output y. Thus, the first step in the controller design is the calculation of the derivatives, which explicitly depend on the system inputs. For this purpose, the Lie-derivative is introduced for simplification. The Lie-derivative Lf hi(x) describes the derivative of the function hi(x) along the vector field f(x)
With this the derivatives of a system output yi can be given as follows.
Thus, the inputs u act directly onto the ri-th derivative of the output yi. ri referred to as relative degree of the output yi. If the ri-th derivative of yi is defined as a new output y*i, then the system description between u and y* given below
follows. The general non-linear rule
u=r(x)+V(x) W (18)
then follows
y*=c*(x)+D*(x)r(x)+D*(x)V(x)·w (19)
for the controlled system. Since y* describes the derivative over time of the output y, the original system (12) is decoupled when the system (17) is decoupled. The setpoint variable w acts via the matrix D*(x)V(x) upon the outputs y*. By choosing
D*(x)V(x)=diag(ki)V(x)=D*(x)−1·diag(ki), (20)
each setpoint variable wi acts only upon the output variable y*i assigned to it, and equation (19) can be given line by line.
Therein, i is the index of the decoupled subsystem, and m is the number of inputs and outputs. In order to calculate the controller parameters rj, the desired transfer function
is defined with the associated time domain representation
yw,i(r
wherein the desired transfer function is stationary accurate, linear and has minimum-phase characteristics. If ki=ai, 0, is selected and
yw,i(r
yw,i(r
is observed, then equating (21) and (23) results in the following relationship between the controller parameters r and the coefficients a of the denominator polynomial of (22)
From (26) the rule
can be derived, which, together with the preliminary filter (20), guarantees a linear, decoupled input and output behavior with a defined pole position for each output y. The decouplability of a control path depends on two system characteristics.
Number of Inputs and Outputs
Since each controlled variable shall be influenced independently from the others, the number of inputs and outputs must be equal.
Invertibility of the Control Path
The decoupling of the control path requires the invertibility of the decouplability matrix D* (17).
Further it is to be noted that due to the decoupling a part of the system can become unobservable. The unobservable system portion is referred to as internal dynamic and is a system characteristic the stability of which is necessary for the realization of the controller. It was realized that an unstable internal dynamic can lead to an unlimited increase of the internal states and thus to a violation of control variable limits or to the destruction of the installation.
The application of the method to the model of the installation (12) results in a first order delay element (relative degree r1=1) for the pressure p1 and a second order delay element (relative degree r2=2) for the volume flow rate Q. The relative degree of the system is determined from the sum of the relative degrees of the partial systems and corresponds here to the system order n=3. Thus, no internal dynamic occurs, and the system is decouplable in a stable manner.
The dead-times in the system are accounted for by using a Smith-Predictor. The basic structure of the Smith-Predictor is comprised of a model that is connected in parallel to the path. This enables feedback of the calculated controlled variable before it can be measured. For the control deviation e
applies. Assuming an ideal model (˜G(s)=G(s) and ˜Td=Td) the dead-time is thus neglected in the controller design. However, since in reality no error free model exists, the control loop is subsequently examined for its robustness with respect to model errors. In particular errors in modelling the dead-times are often described as critical. Stability examinations and criteria for linear systems can be found in the publication of Z. Palmor, Stability properties of Smith dead-time compensator controllers, Int. J Control, 32-6:937-949, 1980. The extension of the Smith-Predictor to linear systems in state space description as well as to a large class of non-linear systems was presented by C. Kravaris und R. A. Wright, Deadtime compensation for nonlinear processes, in the journal AlChE, 35-9:1535-1541, 1989. This work stated that the stability of the uncontrolled system as well as a stable zero-dynamic of the dead-time free system components are required as limitations.
Given the above conditions a state-space controller is designed for a dead-time free SISO-system (single input, single output) for the linearization of the I/O-behavior. The state-space variables necessary for the control are determined by means of a dead-time free model of the open path, referred to by C. Kravaris und R. A. Wright also as “Overall State Predictor”. Some extensions are performed owing to the MIMO-system (multiple input, multiple output) on hand as well as to the fact that the dead-times of the open control path are not equal for all outputs.
Owing to the use of the path model the manipulated variables can be calculated without interference of dead-times. The decoupling of the outputs generally requires a synchronous modification of the system inputs. Since the components of the system under consideration have different dead-times, see
Furthermore, a modification of the feedback is performed in the outer loop. Similar to the classical Smith-Predictor, not the output variable but the difference between measured and calculated output is used for feedback. Here too the dead-time of the two real outputs must be “synchronized” with the model outputs. The advantage of this structure lies in the fact that for a suitable choice of the inner control loop a compensation of disturbances is possible. The inner control loop consists of a decoupling controller and a model. It is designed for stationary accuracy with respect to the setpoint variable ŵ by means of the preliminary filter V(x)
ym(t→∞)=ŵ(t→∞) (29)
Owing to model uncertainties and disturbances the output variable of the path ys will differ from the one of the model ym. The difference e of the two outputs
e=ys−ym (30)
is fed back to the controller input. Then the setpoint variable ŵ can be specified in general terms by
ŵ=w−e=w−(ys−ym) (31).
Through transformation of (31) and substitution of (29) we find
for the stationary terminal value of the output variable. It becomes apparent that an integrating behavior with respect to disturbances is present in the outer loop owing to the stationary accuracy of the inner loop.
The design of the controller was chosen such that the two decoupled transfer functions have the same sum time constants and the dampening of the PT2—system corresponds to the aperiodic limiting case. In this way a similar dynamic behavior of both outputs is achieved despite a different relative degree. In consideration of the manipulated variable limits the following transfer functions of the closed control loop were realized:
At first the controlled system is tested for its response behavior. Through the combination of the final control devices valve and pump, new degrees of freedom have emerged. The controlled variables pressure and volume flow rate of a consumer can be controlled independently from each other. The controller concept was successfully implemented and verified on a test bench with dead-times. The decoupling within the entire working range was successful, and disturbances due to model errors could be compensated by means of the modified Smith-Predictor.
In
In
With the preferred embodiment a tested concept for the independent control of pressure and volume flow rate of a processing installation is presented. The dynamic of such installations is strongly non-linear, wherein the controlled variables are dynamically coupled to each other. Dead-times occur due to the cycle times of the used instruments. For the control of such a system an extended Smith-Predictor is used in combination with a non-linear decoupling controller. The controlled variables and thus the setpoint variables can be chosen independently of each other. This opens new ways of process control.
Although a preferred exemplary embodiment is shown and described in detail in the drawings and in the preceding specification, it should be viewed as purely exemplary and not as limiting the invention. It is noted that only the preferred exemplary embodiment is shown and described, and all variations and modifications that presently or in the future lie within the protective scope of the invention should be protected.
Number | Date | Country | Kind |
---|---|---|---|
10 2010 013 568 | Mar 2010 | DE | national |
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/EP2011/001604 | 3/30/2011 | WO | 00 | 11/2/2012 |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2011/120689 | 10/6/2011 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
6076951 | Wang | Jun 2000 | A |
6155790 | Pyotsia et al. | Dec 2000 | A |
6185468 | Wells | Feb 2001 | B1 |
6474951 | Stephan et al. | Nov 2002 | B2 |
20010054395 | Stephan et al. | Dec 2001 | A1 |
20070212230 | Stavale et al. | Sep 2007 | A1 |
Number | Date | Country |
---|---|---|
10007088 | Aug 2001 | DE |
69911306 | Jul 2004 | DE |
Entry |
---|
Costas Kravaris et al., Deadtime Compensation for Nonlinear Processes, AICHE Journal, Sep. 1989, vol. 35, No. 9, pp. 1535-1542. |
Gerhard H. Nygaard et al., Automatic Coordinated Control of Pump Rates and Choke Valve for Compensating Pressure Fluctuations during Surge and Swab Operations, IADC/SPE 108344, Mar. 28-29, 2007, pp. 1-6. |
Z. Palmor, Stability properties of Smith dead-time compensator controllers, Int. J. Control, 1980, vol. 32, No. 6, pp. 937-949. |
Number | Date | Country | |
---|---|---|---|
20130037121 A1 | Feb 2013 | US |