Model-free online recursive optimization method for batch process based on variable period decomposition

Abstract

The present invention discloses a model-free online recursive optimization method for a batch process based on variable period decomposition. Variable operation data closely related to product quality is acquired, optimization action on each subset is integrated on the basis of time domain variable division on the process by utilizing a data driving method and a global optimization strategy is formed, based on which an online recursive error correction optimization strategy is implemented. According to the method, the online optimization strategy is formed completely based on the operation data of the batch process without needing prior knowledge or a model of a process mechanism. Meanwhile, the optimized operation locus line has better adaptability by using the online recursive correction strategy, and thus the anti-interference requirement of the actual industrial production is better met.

Claims

1. A method of performing a batch reaction in a reactor, the batch reaction comprising N time periods, N being an integer larger than one, the method comprising: for i-th time period of the N time periods, i=2, . . . , N, determining, using a computer, a target value J.sub.0(i) of a temperature in the reactor as a sum of: a difference E(i1) between a reference value J(i1) of the temperature for (i1)-th time period and an actual value RV(i1) of the temperature in the (i1)-th time period, and a reference value J(i) of the temperature for the i-th time period; controlling the reactor so that the temperature in the i-th time period is at the target value J.sub.0(i).

2. The method of claim 1, further comprising, for k-th time period of the N time periods, k=1, . . . , N, determining the reference value J(k), without using a model characterizing the batch reaction.

3. The method of claim 1, wherein the N time periods have the same length.

4. The method of claim 1, further comprising, for k-th time period of the N time periods, k=1, . . . , N, determining the reference value J(k) by: measuring actual values AV(m, k) of the temperature in the k-th time period and an actual value of yield Q(m) of the batch reaction, in m-th run among M actual runs of the batch reaction, m=1, . . . , M; determining the reference value J(k) based on one or more statistical parameters of AV(m, k), m=1, . . . , M.

5. The method of claim 4, further comprising removing actual values AV(m, k) measured in one of the M actual runs by performing a principal component analysis on the actual values AV(m, k) with respect to m.

6. The method of claim 4, wherein determining the reference value J(k) based on the one or more statistical parameters of AV(m, k), m=1, . . . , M, comprises computing the reference value J(k) as M(k)+sign(k)3(k), wherein M(k) is the mean value of AV(m, k), m=1, . . . , M, wherein (k) is the standard deviation of AV(m, k), m=1, . . . , M, and wherein sign(k) has values 1, 1 and 0 respectively when AV(m, k) and Q(m), m=1, . . . , M, have a positive correlation, a negative correlation and no correlation.

7. The method of claim 4, further comprising smoothing J(k) with respect to k.

Description

BRIEF DESCRIPTION OF FIGURES

(1) FIG. 1 is a temperature curve example of a batch process.

(2) FIG. 2 is a principal component mode diagram indicating that the temperature of a batch process is an optimization variable.

(3) FIG. 3 is a composition diagram of period variables.

(4) FIG. 4 is a principal component load diagram of period variables and indicator variables.

(5) FIG. 5 is an action classification diagram of the period variables on the indicator variables.

(6) FIG. 6 is a comparison diagram of an optimized temperature curve and an original temperature curve of a batch process.

(7) FIG. 7 is a generation diagram of an online recursive error correction strategy.

(8) FIG. 8 is a block diagram of implementation steps of the present invention.

(9) FIG. 9 shows an optimized curve after moving average filtering and an original optimized curve.

(10) FIG. 10 is an optimization result (partial) diagram of a batch crystallization process.

DETAILED DESCRIPTION

(11) A batch crystallization process is taken as the example, and the method does not limit the scope of the present invention.

(12) This implementation method is divided into four parts. The first part is data acquisition and preprocessing. The second part is construction of a combined data matrix. The third part is calculation of a basic optimization strategy. The fourth part is establishment of a recursive error correction online optimization strategy.

(13) The block diagram of the implementation steps of the present method is shown as FIG. 8, and the specific implementation steps and algorithms are as follows.

(14) Step 1: For operating a complete batch crystallization process, operation temperature closely related to product yield is selected as a variable to be optimized, and 50 groups of temperature variables and final yield indicator data are acquired in batches. The acquisition time interval of the data is 1 minute. FIG. 1 is a temperature curve data acquisition example of a batch crystallization process, and for the sake of clarity, only the temperature curves of 2 batches are drawn in the figure.

(15) Step 2: For all the acquired 50 batches of temperature data, principal component analysis is performed on the temperature variables in batches, and singular points are removed from a principal component mode diagram, so that all data points are within one degree of credibility. FIG. 2 is a principal component mode diagram indicating that the temperature of a batch process is an optimization variable, and it can be seen from the figure that a batch of temperature data on the right is greatly different from the overall data mode and thus should be removed.

(16) Step 3: The remaining 49 batches of temperature data are divided into 300 periods at equal intervals on a time axis to constitute 300 period variables C1, C2, . . . , C300. For the sake of clarity, FIG. 3 shows C40 to C70 period variables.

(17) Step 4: Each corresponding batch of yield indicator data in step 3 forms an indicator variable Q.

(18) Step 5: The 300 period variables C1, C2, . . . , C300 and one indicator variable Q formed in step 3 and step 4 are combined to generate a 49301-dimensional combined data matrix L.

(19) Step 6: Principal component analysis is performed on the combined matrix L to form a principal component load diagram. For the sake of clarity, FIG. 4 shows a principal component load diagram example generated by combining 25 period variables C36 to C60 and the indicator variable Q.

(20) Step 7: The action directions and magnitudes of the period variables on the indicator variable are classified for the principal component load diagram in step 6. FIG. 5 is a classification example, it can be seen from FIG. 5 that the actions of C154, C155, C156 and C273 on the indicator variable Q are maximum, wherein C154, C155 and C156 are reverse actions, and C273 is a positive action. C66, C111 and the like having an included angle of about 90 degrees with the indicator variable Q in the directions nearly do not act on the indicator variable Q.

(21) Step 8: Mean value and standard deviation of each period variable are calculated respectively. For example, the mean value of C154 having a reverse action on the indicator variable Q is 134.58 DEG C., and the standard deviation is 6.08 DEG C.

(22) Step 9: The optimization target value of the ith period variable is acquired according to the following perturbation calculation formula:
J(i)=M(i)+sign(i)3(i)

(23) wherein J(i), M(i) and (i) herein are respectively optimization target value, mean value and standard deviation of the ith period variable; and sign(i) is a cosine symbol of an included angle formed by the ith period variable and the indicator variable. On the classification diagram of FIG. 5, the sign(i) is +1 when the included angle is smaller than 90 degrees, 1 when the included angle is greater than 90 degrees, and 0 when the included angle is equal to 90 degrees.

(24) Step 10: The optimization target values of all periods obtained in step 9 constitute a basic optimization variable curve according to a period sequence i=1, 2, . . . , 300.

(25) Step 11: Moving average filtering is performed on the basic optimization curve, so that the filtered optimization curve is relatively smooth and facilitates later tracking control design. FIG. 6 is a comparison of the optimized temperature curve and the original temperature curve, and FIG. 9 shows an optimization curve after moving average filtering and an original optimization curve. It can be seen from FIG. 9 that the filtered optimization curve is smoother and facilitates implementation of a tracking controller.

(26) Step 12: When the basic optimization control locus obtained by the above series of steps is used on line, recursive error correction is performed in each time period:

(27) (1) for the (i1)th time period, the error of the offline basic optimization target value J(i1) and the actual measured value RV(i1) is calculated:
E(i1)=J(i1)RV(i1).

(28) (2) on the offline basic optimization strategy, a new optimization target value of next period is constituted:
J.sub.o(i)=J(i)+E(i1).

(29) Step 12 is sequentially calculated according to the period sequence i=1, 2, . . . , 300, till the operation of the whole batch process is over. FIG. 7 is a generation calculation schematic diagram of an online recursive error correction strategy.

(30) FIG. 10 is an optimization result example of a batch crystallization process. It can be seen from the result in the figure that the optimization-free yield is 90.25%, whereas under the recursive correction optimization strategy, the actual operation yield is 94.88% and is substantially close to the theoretical optimal yield 95.51%. This result shows the effectiveness and practicability of the method of the present invention.

(31) While the present invention has been described in some detail for purposes of clarity and understanding, one skilled in the art will appreciate that various changes in form and detail can be made without departing from the true scope of the invention. All figures, tables, appendices, patents, patent applications and publications, referred to above, are hereby incorporated by reference.