System and method for providing raw mix proportioning control in a cement plant with a gradient-based predictive controller

Abstract

A system and method for providing raw mix proportioning control in a cement plant with a gradient-based predictive controller. A raw mix proportioning controller determines the correct mix and composition of raw materials to be transported to a mixer. The raw mix proportioning controller uses a gradient-based predictive controller to determine the proper mix and composition of raw materials. The gradient-based predictive controller takes targeted set points and the chemical composition of the raw material as inputs and generates the proportions of the raw material to be provided as an output for the next time step. The output is generated by using a non-linear constrained optimization.

Description

BACKGROUND OF THE INVENTION

This invention relates generally to a cement plant and more particularly to providing raw mix proportioning control in a cement plant.

A typical cement plant uses raw material such as limestone, sandstone and sweetener to make cement. Transport belts (e.g. weighfeeders) transport each of the three raw materials to a mixer which mixes the materials together. A raw mill receives the mixed material and grinds and blends it into a powder, known as a "raw mix". The raw mill feeds the raw mix to a kiln where it undergoes a calcination process. In order to produce a quality cement, it is necessary that the raw mix produced by the raw mill have physical properties with certain desirable values. Some of the physical properties which characterize the raw mix are a Lime Saturation Factor (LSF), a Alumina Modulus (ALM) and a Silica Modulus (SIM). These properties are all known functions of the fractions of four metallic oxides (i.e., calcium, iron, aluminum, and silicon) present in each of the raw materials. Typically, the LSF, ALM and SIM values for the raw mix coming out of the raw mill should be close to specified set points.

One way of regulating the LSF, ALM and SIM values for the raw mix coming out of the raw mill to the specified set points is by providing closed-loop control with a proportional controller. Typically, the proportional controller uses the deviation from the set points at the raw mill as an input and generates new targeted set points as an output for the next time step. Essentially, the closed-loop proportional controller is a conventional feedback controller that uses tracking error as an input and generates a control action to compensate for the error. One problem with using the closed-loop proportional controller to regulate the LSF, ALM and SIM values for the raw mix coming out of the raw mill is that there is too much fluctuation from the targeted set points. Too much fluctuation causes the raw mix to have an improper mix of the raw materials which results in a poorer quality cement. In order to prevent a fluctuation of LSF, ALM and SIM values for the raw mix coming out of the raw mill, there is a need for a system and a method that can ensure that there is a correct mix and composition of raw materials for making the cement.

BRIEF SUMMARY OF THE INVENTION

In a first embodiment of this invention there is a system for providing raw mix proportioning control in a cement plant. In this embodiment, there is a plurality of raw material and a plurality of transport belts for transporting the material. A raw mix proportion controller, coupled to the plurality of raw material and the plurality of transport belts, controls the proportions of the raw material transported along the transport belts. The raw mix proportion controller comprises a gradient-based predictive controller that uses a plurality of target set points and the composition of the plurality of raw material as inputs and generates a control action to each of the plurality of transport belts that is representative of the proportions of the material to be transported along the belt. A mixer, coupled to the plurality of transport belts, mixes the proportions of each of the plurality of raw material transported therefrom.

In a second embodiment of this invention there is a method for providing raw mix proportioning control in a cement plant. In this embodiment, a plurality of raw material are transported with a plurality of transport belts to a mixer. Proportions of the plurality of raw material transported along the plurality of transport belts to the mixer are controlled by obtaining a plurality of target set points and the composition of the plurality of raw material. A gradient-based predictive control is performed on the plurality of target set points and the composition of the plurality of raw material. The proportions of the plurality of raw material transported along the plurality of transport belts to the mixer are determined according to the gradient-based predictive control. The determined proportions of the plurality of raw material are sent to the mixer for mixing.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 shows a block diagram of a system for providing raw mix proportioning control in a cement plant according to this invention;

FIG. 2 shows a schematic of the gradient-based predictive control provided by the raw mix proportioning controller shown in FIG. 1 according to this invention;

FIG. 3 shows a more detailed schematic of the open-loop system shown in FIG. 2;

FIG. 4 shows a schematic depicting geometric interpretation of the gradient-based predictive control performed according to this invention; and

FIG. 5 shows a flow chart setting forth the steps of using gradient-based predictive control to provide raw mix proportioning according to this invention.

DETAILED DESCRIPTION OF THE INVENTION

FIG. 1 shows a block diagram of a system 10 for providing raw mix proportioning control in a cement plant according to this invention. The raw mix proportioning control system 10 comprises a plurality of raw material 12 such as limestone, sandstone and sweetener to make cement. In addition, moisture can be added to the raw materials. While these materials are representative of a suitable mixture to produce a cement raw mix, it should be clearly understood that the principles of this invention may also be applied to other types of raw material used for manufacturing cement raw mix. Containers 14 of each type of raw material move along a transport belt 16 such as a weighfeeder. A raw mix proportioning controller 18 controls the proportions of each raw material 12 transported along the transport belts 16. A mixer 20 mixes the proportions of each raw material 12 transported along the transport belts 16. A raw mill 22 receives mixed material 24 from the mixer 20 and grinds and blends it into a raw mix. The raw mill 22 feeds the raw mix to a kiln 26 where it undergoes a calcination process.

As mentioned above, it is necessary that the raw mix produced by the raw mill 22 have physical properties with certain desirable values. In this invention, the physical properties are the LSF, ALM and SIM. These properties are all known functions of the fractions of four metallic oxides (i.e., calcium, iron, aluminum, and silicon) present in each of the raw materials. A sensor 28, such as an IMA QUARCON.TM. sensor, located at one of the transport belts 16 for conveying the limestone, measures the calcium, iron, aluminum and silicon present in the limestone. Those skilled in the art will recognize that more than one sensor can be used with the other raw materials if desired. Typically, the LSF, ALM and SIM values for the raw mix coming out of the raw mill should be close to specified target set points. Another sensor 30 such as an IMA IMACON.TM. sensor located before the raw mill 22 measures the calcium, iron, aluminum and silicon present in the mix 24. Although this invention is described with reference to LSF, ALM and SIM physical properties, those skilled in the art will recognize that other physical properties that characterize the raw mix are within the scope of this invention.

The raw mix proportioning controller 18 continually changes the proportions of the raw material 12 in which the material are mixed prior to entering the raw mill 22 so that the values of LSF, ALM and SIM are close to the desired set points and fluctuate as little as possible. The raw mix proportioning controller 18 uses gradient-based predictive control to continually change the proportions of the raw material. In particular, the gradient-based predictive control uses targeted set points and the chemical composition of the raw material as inputs and generates control actions to continually change the proportions of the raw material. The mixer 20 mixes the proportions of the raw material as determined by the gradient-based predictive control and the raw mill 22 grinds the mix 24 into a raw mix.

FIG. 2 shows a schematic of the gradient-based predictive control provided by the raw mix proportioning controller 18. There are two main components to the gradient-based predictive control provided by the raw mix proportioning controller; a gradient-based predictive controller 32 and an open-loop system 34. The gradient-based predictive control takes S* and P as inputs and generates S as an output, where S* is the targeted set points, P is the process composition matrix, and S is the actual set points. A more detailed discussion of these variables is set forth below. At each time step, the gradient-based predictive control attempts to eliminate the tracking error, which is defined as;

.DELTA.S(t)=S*-S(t) (1) by generating .DELTA.U(t), the change in control action, which results in proper control action for the next time step which is defined as: U(t+1)=.DELTA.U(t)+U(t) (2)

More specifically, the gradient-based predictive controller 32 uses gradient information to produce change in control to compensate the tracking error. In FIG. 2, a subtractor 31 performs the operation of equation 1 and a summer 33 performs the operation of equation 2.

FIG. 3 shows a more detailed diagram of the open-loop system 34 shown in FIG. 2. The open-loop system 34 receives P and U as inputs and generates S as an output, where P is a process composition matrix of size 4 by 3, U is a control variable matrix of size 3 by 1, S is the actual set point matrix of size 3 by 1, and R is a weight matrix of size 4 by 1.

The process composition matrix P represents the chemical composition (in percentage) of the input raw material (i.e., limestone, sandstone and sweetener) and is defined as: ##EQU1## Column 1 in matrix P represents the chemical composition of limestone, while columns 2 and 3 in P represent sandstone and sweetener, respectively. This invention assumes that only column 1 in P varies over time, while columns 2 and 3 are considered constant at any given day. Row 1 in matrix P represents the percentage of the chemical element CaO present in the raw material, while rows 2, 3, and 4 represent the percentage of the chemical elements S.sub.i O.sub.2, Al.sub.2 O.sub.3 and Fe.sub.2 O.sub.3, respectively, present in the raw materials.

The control variable vector U represents the proportions of the raw material (i.e., limestone, sandstone and sweetener) used for raw mix proportioning. The matrix U is defined as: ##EQU2##

The set point vector S contains the set points LSF, SIM and ALM and is defined as: ##EQU3##

The weight matrix R is defined as: ##EQU4## wherein C, S, A and F are the weight of CaO, S.sub.i O.sub.2, Al.sub.2 O.sub.3 and Fe.sub.2 O.sub.3, respectively, and R is derived by multiplying U by P. A function .function. takes R as input and generates S as output. The function .function. comprises three simultaneous non-linear equations defined as follows: ##EQU5## where:

C=c.sub.1 .multidot.u.sub.1 +c.sub.2 .multidot.u.sub.2 +c.sub.3 .multidot.(1-u.sub.1 -u.sub.2) S=s.sub.1 .multidot.u.sub.1 +s.sub.2 .multidot.u.sub.2 +s.sub.3 .multidot.(1-u.sub.1 -u.sub.2) A=a.sub.1 .multidot.u.sub.1 +a.sub.2 .multidot.u.sub.2 +a.sub.3 .multidot.*1-u.sub.1 -u.sub.2) F=.function..sub.1 .multidot.u.sub.1 +.function..sub.2 .multidot.u.sub.2 +.function..sub.3 .multidot.(1-u.sub.1 -u.sub.2) (10)(11)(12)(13) and u.sub.1, u.sub.2 and u.sub.3 =1-u.sub.1 -u.sub.2 are the dry basis ratio of limestone, sandstone and sweetener, respectively. Furthermore, c.sub.i, s.sub.i, a.sub.i and f.sub.i are the chemical elements of process matrix P defined in equation 3.

Referring back to FIG. 2, the gradient-based predictive controller 32 maps .DELTA.S to .DELTA.U. The mapping of .DELTA.S to .DELTA.U is defined as:

.DELTA.U=.function..smallcircle..DELTA.S (14) wherein .smallcircle. is a function mapping and .function. performs the function mapping of the gradient-based predictive controller 32. The open-loop system gradient can then be defined as: ##EQU6##

The open-loop system gradient defined in equation 15 can then be reorganized and represented as follows:

R=PxU S=.function.f(R) (16) (17) wherein x represents matrix multiplication, R is the weight matrix, P is the process composition matrix, U is the control variable matrix, and S is the set point matrix.

Assuming that P changes insignificantly between two control iterations, then two gradients, G.sub.1 and G.sub.2, which calculate the partial derivatives of .DELTA.S with respect to .DELTA.R and .DELTA.R with respect to .DELTA.U, respectively, can be derived and are defined as:

.DELTA.S=G.sub.1 x.DELTA.R .DELTA.R=G.sub.2 x.DELTA.U (18)(19) wherein .DELTA.S is a 3 by 1 matrix defined as: ##EQU7## .DELTA.R is a 4 by 1 matrix defined as: ##EQU8## .DELTA.U is a 2 by 1 matrix defined as: ##EQU9## G.sub.1 is a 3 by 4 matrix defined as: ##EQU10##

G.sub.2 is a 4 by 2 matrix defined as: ##EQU11## Note that C.sub.i, s.sub.i, a.sub.i and f.sub.i are the chemical elements as defined in equation 3. Further, the partial derivatives of .DELTA.S with respect to .DELTA.U can be obtained from: ##EQU12## where G=G.sub.1 x G.sub.2 is a 3 by 2 matrix.

As mentioned above, the functionality of the gradient-based controller 32 is to invert the G matrix. More specifically, the functionality is to find G.sup.-1 such that:

.DELTA.U=G.sup.-1 x.DELTA.S (25) However, G is not a square matrix and is generally not directly invertible. This results in an over-constrained problem and usually renders its solution to pseudo-inversion or optimization techniques such a least means squares.

FIG. 4 shows a schematic depicting geometric interpretation of the gradient-based predictive control performed according to this invention. In particular, FIG. 4 shows three lines lying on a plane in a two-dimensional space. The three lines represent .DELTA.LSF, .DELTA.SIM and .DELTA.ALM and are in the two-dimensional space spanned by .DELTA.u.sub.1 and .DELTA.u.sub.2. In FIG. 4 the lines for .DELTA.LSF, .DELTA.SIM and .DELTA.ALM are labeled as .DELTA.LSF, .DELTA.SIM and .DELTA.ALM, respectively. The points on .DELTA.LSF represent the change in control action which is able to bring the system to the change in set point, .DELTA.LSF. Similarly, the points on .DELTA.SIM and .DELTA.ALM represent the change in control actions which are able to bring the system to the change in set points, .DELTA.SIM and .DELTA.ALM, respectively.

Reaching the change in three set points simultaneously means that there exists a point on the plane which is on .DELTA.LSF, .DELTA.SIM and .DELTA.ALM. This can be interpreted as where the sum of distance from the point to the three lines is minimized. Similarly, there will be only two lines on the plane if there are two change in set points. To find a change in control action to reach the two change in set points at the same time is equivalent to finding the point on the plane at which the two lines meet. This again could be interpreted as where the sum of distance from the point to the two lines is minimized. In general, the distance from a point (a change in control action) to a line (a change in set point) can be interpreted as the degree of unreachability for the change in control action to reach the change in set point. The shorter the distance, the greater the degree of reachability. The longer the distance, the less the degree of reachability. In this context, to what degree a change in control action (a point on the plane) drives the system to a specific change in set point (a line on the plane) depends on how far the point is from the line.

In order to provide the gradient-based predictive control according to this invention, the gradient-based predictive controller 32 uses an optimization algorithm to determine a sequence of future controller outputs over a control horizon, such that a specified objective function is minimized. In this invention, the objective function is a modification of the quadratic function which is defined as: ##EQU13## wherein J is the objective function to be minimized, S* and S are the target and the actual set points, respectively, H and H.sub.c are the prediction and control horizons, respectively, .alpha. and .beta. define the weighting of the tracking error and the control effort with respect to each other and with respect to time, k is the current time step, i is the time step index, .DELTA.S is defined in equation 25, .function.(.multidot.) represents the functionality of the open-loop system, P is the process composition matrix, and U.sup.1 and U.sup.u are the lower and upper bounds of .DELTA.U, respectively.

In essence, the first term of J minimizes the tracking error and the second term of J minimizes control jockeying in order to provide smooth changes in control as opposed to abrupt changes. Thus, J seeks a balance between minimizing tracking error and maintaining smooth control, while the tradeoff is controlled by .alpha. and .beta.. Note that it is not possible to satisfy all three equations in equation 25 simultaneously. At most, two out of the three can be satisfied at the same time. It is therefore, up to the choice of the user, which depends on the priority of the set points. Furthermore, only the first control U(k) is applied to the system and the optimization is repeated at the next time step k+1, which known as the receding horizon principle.

In this invention, MATLAB, a well-known scientific computing software, is used for fast prototyping and simulation of the constrained optimization. MATLAB's non-linear constrained optimization routines use a Sequential Quadratic Programming (SQP) method which is a form of gradient descent, which finds a local optima to the problem. To find the local optima it is assumed that the objective function and constraints are non-linear. The explicit constraints are assumed to be inequality constraints since the parameters are bounded from below and above. The objective function is approximated by a quadratic function. This is done by approximating its Hessian at the current point. The non-linear constraints are linearly approximated locally. The approximation produces a quadratic programming problem, which can be solved by any of several standard methods. The solution is used to form a new iterate for the next step. The step length to the next point is determined by a line search, such that a sufficient decrease in the objective function is obtained. The Hessian and constraint planes are then updated appropriately and this method is iterated until there is no appropriate non-zero step length to be found.

FIG. 5 shows a flow chart describing the raw mix proportioning control of this invention. Initially, the raw mix proportioning controller obtains a plurality of target set points S* at 36. Next, the raw mix proportioning controller obtains the process composition matrix P at 38. The raw mix proportioning controller then performs the gradient-based predictive control by using the above described optimization at 40. The raw mix proportioning controller then outputs the control matrix U at 42 which is the proportion of raw materials. The raw mix proportioning controller then sets the speed of each of the transport belts to provide the proper proportion of raw material at 44 which is in accordance with the control matrix U. These steps continue until the end of the production shift. If there is still more time left in the production shift as determined at 46, then steps 36-44 are repeated, otherwise, the process ends.

It is therefore apparent that there has been provided in accordance with the present invention, a system and method for providing raw mix proportioning control in a cement plant with a gradient-based predictive controller that fully satisfy the aims and advantages and objectives hereinbefore set forth. The invention has been described with reference to several embodiments, however, it will be appreciated that variations and modifications can be effected by a person of ordinary skill in the art without departing from the scope of the invention.

Claims

1. A system for providing raw mix proportioning control in a cement plant, comprising:

a plurality of raw material;

a plurality of transport belts for transporting the plurality of raw material;

a measuring device that measures the composition of the plurality of raw material transported by the plurality of transport belts;

a raw mix proportioning controller, coupled to the plurality of transport belts and the measuring device, for controlling the proportions of the plurality of raw material transported along the plurality of transport belts, wherein the raw mix proportioning controller comprises a gradient-based predictive controller that uses a plurality of target set points and the composition of the plurality of raw material as inputs and generates a control action to each of the plurality of transport belts that is representative of the proportions of the material to be transported along the belt, wherein the gradient-based predictive controller determines a sequence of future control outputs over a control horizon to generate the control action; and

a mixer, coupled to the plurality of transport belts, for mixing the proportions of each of the plurality of raw material transported therefrom.

2. The system according to claim 1, wherein the plurality of raw material comprise limestone, sandstone and sweetener.

3. The system according to claim 1, wherein the plurality of target set points are physical properties comprising lime saturation factor, alumina modulus and silica modulus.

4. The system according to claim 1, wherein the gradient-based predictive controller performs a non-linear constrained optimization.

5. The system according to claim 4, wherein the gradient-based predictive controller minimizes a specified objective function to minimize tracking error and control jockeying.

6. The system according to claim 1, wherein the system further comprises a raw mill, coupled to the mixer for grinding and blending the mix of the plurality of raw material into a raw mix.

7. The system according to claim 6, wherein the system further comprises a kiln, coupled to the raw mill for burning the raw mix.

8. The system according to claim 1, wherein the gradient-based predictive controller performs a geometric interpretation between the plurality of target set points and the composition of the plurality of raw material.

9. A method for providing raw mix proportioning control in a cement plant, comprising:

providing a plurality of raw material;

transporting the plurality of raw material with a plurality of transport belts to a mixer;

controlling the proportions of the plurality of raw material transported along the plurality of transport belts to the mixer, comprising:

obtaining a plurality of target set points;

obtaining the composition of the plurality of raw material;

performing gradient-based predictive control on the plurality of target set points and the composition of the plurality of raw material, wherein the gradient-based predictive control comprises determining a sequence of future controller outputs over a control horizon; and

determining the proportions of the plurality of raw material transported along the plurality of transport belts to the mixer according to the gradient-based predictive control; and

mixing the determined proportions of the plurality of raw material with the mixer.

10. The method according to claim 9, further comprising providing the mix of the plurality of raw material from the mixer to a raw mill and generating a raw mix therefrom.

11. The method according to claim 10, further comprising providing the raw mix from the raw mill to a kiln.

12. The method according to claim 9, wherein performing the gradient-based predictive control comprises performing a non-linear constrained optimization.

13. The method according to claim 12, wherein the performing of the gradient-based predictive control further comprises minimizing a specified objective function to minimize tracking error and control jockeying.

14. The method according to claim 9, wherein the plurality of raw material comprise limestone, sandstone and sweetener.

15. The method according to claim 9, wherein the plurality of target set points are physical properties comprising lime saturation factor, alumina modulus and silica modulus.

16. The method according to claim 9, wherein the performing of the gradient-based predictive control comprises performing a geometric interpretation between the plurality of target set points and the composition of the plurality of raw material.