A SYSTEM FOR ON LINE INFERENCE OF PHYSICAL AND CHEMICAL PROPERTIES AND SYSTEM FOR ON LINE CONTROL

Abstract: A system for on line inference and control of physical and chemical properties of polypropylene and its copolymers is described. The system comprises models for the inference of physical and chemical properties that are not continuously measured and relevant models to control these properties as well as production rate, density of the reaction medium and other process variables of interest. The described control system allows to maximize production rate as well as catalyst yield in the producing process.

Description: A SYSTEM FOR ON LINE INFERENCE OF PHYSICAL AND CHEMICAL PROPERTIES AND SYSTEM FOR ON LINE CONTROL FIELD OF THE INVENTION The present invention relates to a system for on line inference and control, of the physical and chemical properties of polypropylene and its copolymers with other olefins produced in loop reactors, and, optionally, other gas-phase reactors, either isolated or combined in series, with the aid of mathematical models. More specifically, said system consists of models to infer physical and chemical properties that are not continuously measured and of relevant models to control these properties as well as other process variables of interest in the process being studied. The control system further enables the maximization of production rate and catalyst yield of producing process. BACKGROUND INFORMATION Usually, control of process variables in petrochemical plants is manually carried out by periodically taking product samples to be tested, while maintenance and correction actions concerning operating conditions are taken in order to obtain a product having desired characteristics. This procedure involves delays related to required corrections, since sampling and analysis are usually delayed as compared to the on line process, not to mention any human mistakes. On the other hand, rigorous and empirical modeling techniques may be used to build process models. These mathematical models are able to infer the value of certain process variables which are measured periodically based on the values of other process variables which are measured continuouslv. Besides, the mathematical models can also be used to predict the future behavior of process variables resulting from modifications ofthe plant operating conditions. Rigorous modeling techniques are based in the natural laws which dictate the fundamental relationships among process variables. Rigorous models are difficult to obtain and demand a high level of knowledge about the process. Besides, the complexity ofthe equations that constitute the rigorous model make it inadequate for on line implementation due to the long time necessary for its execution, even using computers. In compensation, the techniques of empirical modeling do not require such a deep knowledge about the process and result in simpler mathematical models, that can be executed quickly being more adequate for real time execution. The disadvantage of empirical models is the fact they cannot be used in operating conditions different from those used in its identification. Linear and nonlinear regression models and neural nets are examples of empirical modeling techniques commonly cited in the open literature. In the systems for inference and control described in present invention, the preferred modeling technique comprises the use of neural networks. Neural networks are networks of elements interconnecting themselves in a unique way. Typically, networks consist of input neurons that receive signals or information from the outside of the network, output neurons that transmit signals or information to the outside of the network, and, at least, one intermediate layer of neurons that receive and pass information on to other neurons. A system for controlling polymeric intrinsic viscosity during its production was described in the US Patent 3878379. However this technology was specific for the production of polyethylene terephtalate polyester and it comprised just a single output, excluding the control of others polymer properties. The open literature contains a series of publications that deal with inference and control of processes in a general way using neural networks. US Patent 51 1 1531 teaches how to implement an on line predictive control system applied to a continuous process where a neural network predicts the value of a controlled variable, the values of manipulated variables being adjusted so that the value predicted by the network closely approaches the desired value. US Patent 5175678 teaches how to implement an on line feedforward control system using neural networks. US Patent 5282261 teaches how to implement an inference and control system in a continuous process using neural networks. US Patent 5479573 teaches how to develop a predictive network to operate in real time, or in training mode. Preprocessed data is used as inputs of the model for the system, which provides a predicted output that is used in the distributed control system. US Patent 5675253 teaches how to use a neural network to develop a model, which correlates process variables together with the results of an on line nuclear magnetic resonance system to predict a desired property ofthe polymer. US Patent 5729661 teaches a method and a device to preprocess input data for a neural network. Data is preprocessed to replace wrong data, or to fill in missing data and then used to train a model. Examples of inference systems that use neural networks applied to specific processes are also found. US Patent 5548528 teaches how to build a system to continuously monitor pollutant emissions. Such system uses the neural networks to replace instruments that perform pollutant analysis in a plant. Then, the value predicted by the neural network is used as a controlled variable in a control system. WO Patent 98/26336 teaches a method to infer the blade temperature of a vapor turbine. This invention uses temperature and pressure from other places of the turbine as inputs of a neural network that predicts blade temperature. However, no published Patents, neither isolated nor in combination with other published literature, explains the on line inference of physical and chemical properties of polypropylene and its copolymers with other olefins produced in loop reactors and, optionally, in other gas- phase reactors, either isolated, or combined in series from measured process variables using mathematical models. Further, no Patent teaches or suggests the on line control of process variables such as production rate, density of reaction medium, and physical and chemical properties of polypropylene and its copolymers with other olefins produced in loop reactors, and, optionally, in other gas-phase reactors, either isolated or combined in series using mathematical models, such on line control and inference being enabled by the system described and claimed in present invention. SUMMARY OF THE INVENTION The present invention, as disclosed and claimed herein, is directed to a system for on line inference of physical and chemical properties of polypropylene and its copolymers produced in a plant using loop reactors, and, optionally, gas-phase reactor(s) with the aid of mathematical models, such system comprising the steps of: - collecting data of all variables measured in the plant, so as to build a large historical database; - out of said large database, selecting a subset including historical data on the physical and chemical properties to be inferred, as well as potential input variables for mathematical models so as to build a small database; - treating data of said small database, removing non- representative data of usual operating conditions, applying noise removal filters on signals and complementing said database with new calculated variables based on theoretical concepts such as energy and mass balances; - from the data of said small database, selecting those variables that will be the input variables ofthe mathematical model; - using said input variables to identify these models off line; - applying the so identified mathematical models to the realtime measurements ofthe input variables to the model used to calculate the on line inference of desired physical and chemical properties. The present invention also provides an on line control system for a plant of polypropylene, and its copolymers produced in loop reactors, and, optionally, gas-phase reactor(s) using mathematical models, the system being based on a control matrix which comprises controlled variables, constrained controlled variables, manipulated variables, and disturbance variables. Thus, the present invention provides a system for the on line inference and control of physical and chemical properties of polypropylene and its copolymers. The present invention also provides a system for the on line control of the production rate of polypropylene and its copolymers, of the density ofthe reaction medium as well as other process variables. Further, the present invention provides an on line control system for a plant of polypropylene and its copolymers with the aid of mathematical models based on a control matrix that enables the maximization ofthe production rate, as well as of the catalyst yield. BRIEF DESCRIPTION OF THE DRAWINGS FIGURE 1 represents the categories of process variables pertaining to a chemical process. FIGURE 2 represents a control matrix. FIGURE 3 represents a simplified plant for the production of polypropylene, to which the present invention is applied. FIGURE 4 represents schematically a system for the on line inference according to the invention. FIGURE 5 represents schematically a control system according to the invention. FIGURE 6 represents a control matrix that characterizes a control system according to the invention. FIGURE 7 shows the increase in production rate of the first loop reactor due to the decrease in the variability. FIGURE 8 shows the decrease in the variability of the melt flow index of the polymer produced in the first loop reactor due to the distribution of the difference between the desired value and the process variable. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS In the present specification, the expression on line means that the physical and chemical properties, which are desired to infer, are obtained in real time. The expression polymer quality refers to the physical and chemical properties of the polymer such as melt flow index, viscosity, chemical composition, etc. that affect end-use properties ofthe resin. The expression SISO (Single Input Single Output) control strategy means that a single manipulated variable is adjusted to control a given controlled variable. The expression MIMO (Multiple Input Multiple Output) control strategy means that multiple manipulated variables are adjusted at the same time so as to control a set of controlled variables. The expression gas-phase reactor means that, at least, one gas phase and one solid phase are present in the reactor. The preferred gas- phase reactor is a fluidized-bed reactor, but other reactors can be used. The expression loop reactor refers to a tubular reactor with a high recycle rate, where the produced polymer is entrained by the flow of a liquid phase. Broadly, the present process applies to reactor arrangements as presented in FIGURE 3. However, other arrangements are also possible as follows: - one or more loop rectors connected in series; and - one or more loop reactors connected in series, followed by one or more gas-phase reactor(s) connected in series. The present process is specially directed to the polypropylene homopolymer, copolymer, or terpolymer. The comonomers, in the case of copolymers or terpolymers, may be selected, besides propylene, among the group of reactive olefins, such as: ethylene; 1-butene; 2-methyl propylene; 1-pentene; 1-hexene; 1-heptene; 1-octene; 1-decene; 3-methyl 1-butene; 4- methyl 1-pentene; cyclopentene, the comonomers being polymerized alone or in admixture. The preferred catalyst is a Ziegler-Natta type catalyst, but others may be used as well. In order to establish a control strategy for polymerization reactors, it is convenient that process variables be classified according to FIGURE 1 as: controlled variables (8), manipulated variables (1), constrained controlled variables, or simply, constraints (6), and disturbance variables, or simply, disturbances (7). Controlled Variables (CV's) are those variables the value of which must be maintained as close as possible to a desired value, or set point. In a polymerization plant, controlled variables can be divided into two categories: physical and chemical properties of the polymer (4), and process variables that define operating conditions ofthe plant (5). The physical and chemical properties of the polymer (4), such as melt flow index, viscosity, chemical composition, etc., define polymer quality. Process variables, which define operating conditions of the plant (5), or simply, operating conditions, involve such properties as temperature, pressure, flow rate, etc. Constrained Controlled Variables (CCV's) are variables that despite being controlled do not require to be maintained close to a set point. However, constrained controlled variables must be controlled so as not to exceed certain limits. Manipulated Variables (MV's) are variables that must be adjusted so that controlled variables are maintained close to a set point or within certain limits. To be used as a manipulated variable in a control strategy, a process variable must affect at least one controlled variable, or a constrained controlled variable in a certain way. Monomer flow rates, catalysts and cooling water are examples of common MVs. Generally, operators who follow up the operation of a chemical process know, at least qualitatively, the way a certain MV affects a given CV. Disturbance Variables (DV's) are process variables that also affect controlled variables, but operators' actions can not affect them. Examples of disturbance variables are: content of impurities in the monomer stream fed to the reactor, and catalyst quality. Disturbances are usually caused by known sources, but it is not easy to eliminate them as they are generated by factors that are external to the process. Some disturbances can be measured, but others can not. A control system may take into consideration the measured disturbances (2), promoting thus an adjustment of the manipulated variables before the effects of disturbances are evidenced in the product quality. It is important to emphasize that from the point of view of control, if a process input is not used as a MV, then it can be considered as being a DV. The key for the successful use of the neural network in inference systems for specific applications is the adequate selection of a set of process variables used as inputs of the neural network. In the case of a process control system, the adequate selection of variables that comprise the control matrix may be more important than the selection of the methodology of process modeling. In case the neural networks are used to represent process models, network inputs and outputs to be used by the control system are defined by the control matrix of the system itself, since the manipulated variables and disturbance variables must be the network inputs, while the controlled variables and constrained variables must be the outputs. As an input for controllers, the present invention uses such values as predicted by the neural networks. The classification of a process variable in one of the above- cited categories is by no means trivial. For example, the hydrogen concentration in the reactor feed may be manipulated to fix the melt flow index of the polymer, and, at the same time, may be the controlled variable in a loop for the control ofthe concentration which uses the hydrogen flow rate as a manipulated variable. In this case, the manipulated variable used to control the melt flow index is the desired value for hydrogen concentration. Such approach is known as cascade control. The function of those models used in process control is really the storage of cause/effect relationships among CV's, MV's, and, in some cases, DV's in a format of mathematical equations . The physical and chemical properties of a polymer as well as the economic features involved in the process, such as production rate and catalyst yield depend on the raw materials used in polymer manufacturing besides the operating conditions of the polymerization reaction. Therefore, once the raw materials to be used in the polymerization are determined, it is very important that the process inputs are manipulated so as to control the physical and chemical properties of the polymer, while keeping high production rates and high catalyst yield. The majority of the physical and chemical properties that define product quality are not measured on line. In many cases, there are no analyzers that can be directly connected to the process, and even when the technology required to measure these variables is available, the cost of sensors may hinder their use. Therefore, analyses are carried out in a laboratory in one-hour intervals. Thus, during the period between analyses, there is no indication of the product quality, causing a delay between possible modifications in such variables and the implementation of any correction actions as required. A further relevant point regarding process control is the fact that, generally, economic aspects require the process to be operated close to plant capacity limits. Thus, catalyst flow rate could be manipulated to control the production rate of the reactor, however, when polymerization heat exceeds the limit of the thermal exchange capacity of the reactor, the objective of increasing production even more should be restricted in favor of thermal stability. In view of In view of this example, the importance of considering CCV's in a multivariable control strategy is clear, as these variables will determine the limits under which plant optimization is achieved. In order to overcome the scarcity of measurements related to polymer quality, process modeling techniques can be used to allow continuous inference of properties that are not continuously available. In the inference and control systems, described in present invention, the preferred modeling technique comprises the use of neural networks as the process models. Optionally, empirical or rigorous modeling techniques or a combination of rigorous and empirical models, including neural networks could be used instead of neural networks to create said process models. Neural networks, as mentioned herein, are networks which consist of interconnected elements that perform elementary mathematical operations. Each element receives either a single, or multiple inputs so as to yield an output used as input for different elements in a process which mimics nervous transmission by neurons in the human brain. Among the several structures suggested for the neural networks, perhaps the one which finds more practical applications is the structure that consists of a layer of neurons that receive external signals or inputs, an intermediate layer that processes a combination of external signals, and an output layer that presents the result of calculations. Each interconnection among neurons is weighted; that is, each interconnection is associated to a weight the value of which is an indication of the degree of importance ofthe information being transmitted. The training process of the neural network comprises the use of the data stored in a historical database to adjust network weights. At first, weight values are randomly chosen, and then adjustments are made so as to minimize the difference between the neural network output (also called prediction), and the actual value of the process variable being examined. As described herein before, polymerization processes present a huge number of inputs and outputs. Cause/effect relationships among different variables are too complex to be assimilated, even by a trained and experienced operator. It could be observed that the most important relationships are assimilated and used in the control of systems composed by only one input and one output (SISO). One example is the control ofthe melt flow index (MFI) ofthe polymer. The variable which mostly affects the melt flow index is the hydrogen concentration in the reactor feed, but there are other variables, such as the reactor temperature, the effect of which on the MFI is not so evident. Typically, upon receiving the results of MFI laboratory analyses, the operator adjusts the hydrogen concentration, without however taking into consideration the temperature of the reactor, since the temperature effect is less evident. Despite the fact that the temperature effect is somehow less important, this type of control is not as effective as it would be if every pertinent cause/effect relationship could be taken into consideration. Therefore, in order to obtain significant benefits, the control system of those processes should be multivariable (MIMO), comprising various controlled variables as well as various manipulated variables. Besides, the control system must take into consideration the physical limits of the process under study as well as any disturbance which could possibly be measured. A control technique that presents these characteristics is the Model Predictive Control which is based on process models to determine the best set of actions to be taken so that the controlled variables reach desired values. US Patent 5282261 teaches the use of predicted values by a neural network in real time so as to replace a sensor or a laboratory analysis as a controller input, implementing a system for inference and control in a continuous process using neural networks. On the other hand, US Patent 5282261 neither teaches how to apply these values to polymerization processes, in special processes for the production of polypropylene and its copolymers, nor teaches which process variables would influence physical and chemical properties of polypropylene and its copolymers. Once the neural network "has learned" the relationships among variables and historical data, it can be used for both inference and control. In the case of inference, the role of the neural network is to provide real-time estimates for variables that are not continuously available. In case the neural network is used for predictive control, the role of the network is to predict the future behavior of the CVs and CCVs if certain adjustments on MVs are done and if DVs vary. A multivariable control system may be precisely characterized by a set of manipulated variables, controlled variables, constraints, and disturbances considered in the modeling of the system. The set of variables determines what may be encompassed by the control system as well as its capacity to be used for optimizing a process. A control matrix may represent a multivariable control system as the one presented in FIGURE 2. In this matrix, columns represent controlled variables (CV's) and constraints (CCV's), while the lines represent manipulated variables (MV's) and measured disturbances (DV's). Marked positions within the control matrix indicate that a given MV, or DV is used in the modeling of a given CV, or CCV. FIGURE 3 illustrates the main unit operations carried out in a polypropylene plant that is the object of the present invention. Stream 9 represents feed streams of catalyst, cocatalyst, electron donor, propylene, and hydrogen to the pre-polymerization reactor, PR. Stream 10 represents the transfer of pre-polymerized catalyst to the first loop reactor, LR1, to which propylene, hydrogen, and optionally, comonomer(s) are fed, represented by stream 1 1. Stream 12 represents the stream which exits LR1, to which is added stream 13, said stream 13 representing the feed of propylene, hydrogen, and, optionally, comonomer(s); both streams, 12 and 13, represent the feed of the second loop reactor, LR2. Stream 14 represents the exit of LR2 to the monomer recovery cyclone, C 1 , which separates the polymer formed in the reactors from excess monomer. Excess monomers are recovered. Stream 15 feeds the gas-phase reactor, GPR. To the GPR is directed stream 16, made up of the recycle stream of the GPR, by stream 19 and by stream 20. Stream 20 may be made up of propylene, comonomer(s), and hydrogen. Stream 17 effects the discharge from the GPR to cyclone C2, from where stream 18 and a polymer stream are separated. The polymer stream is sent to degassing, deactivation, and extrusion. Stream 18 is fed to the separation tower, ST. Two streams leave the ST, a propylene stream, and stream 19 which is fed to the GPR. FIGURE 3 illustrates one process configuration, however, the present invention is not restricted to this configuration only, as it applies to other arrangements of loop reactor(s), and, optionally, to gas-phase reactor(s). Among the main objectives of a plant which produces polypropylene and its copolymers is the manufacturing of polymers with physical and chemical properties within certain specifications, and the optimization of economic objectives such as specific consumption of catalyst, and plant production rate, that is, the purpose is to minimize catalyst consumption per ton of polymer produced, and to maximize the ton/h ratio ofthe polymer produced. Traditionally, in such a plant, reactors are operated according to pre-established procedures, and laboratory analyses indicate if it is necessary to make adjustments in operating conditions so that product quality is maintained within specifications. Laboratory analyses involve measurements such as: melt flow index, percentage(s) of comonomer(s) incorporated in the copolymer, percentage of xylene-soluble matter, intrinsic viscosity, and percentage(s) of comonomer(s) incorporated in the fraction of polymer produced in the gas-phase reactor. Loop reactors contain a reaction medium made up of liquid propylene, and solid polymer. The density of this reaction medium is measured on line by a density gage, and it is controlled by a SISO PID control algorithm. The density of the reaction medium is inversely proportional to the catalyst consumption per ton of polymer produced. The reaction medium is made to circulate within the loop reactors with the aid of pumps. The solid phase contains a small fraction of catalyst that remains intimately connected to the polymer. Small amounts of hydrogen are added to the liquid phase of the reaction medium so as to modify the molecular weight of the polymer. The use of different concentrations of hydrogen in each loop reactor results in a bimodal polymer. The production rate and the polymer quality are controlled in an open loop by means of SISO systems. The control of the physical and chemical properties of the polymer when the gas-phase reactor is operating is based on the open-loop adjustment of the concentrations of comonomer(s), propylene, and hydrogen within the reactor, or on the open-loop adjustment of operating conditions of the loop reactors. The gas composition within the gas-phase reactor is measured at intervals of a few minutes by gas chromatography, said composition being controlled by SISO PID control algorithms. A heat exchanger removes reaction heat from the recycle stream. As in the case of loop reactors, the heat which is exchanged depends on the flow rate of the cooling water which flows through the heat exchanger. Therefore, the state-of-the-art control strategy in a plant that produces homopolymer polypropylene and its copolymers with varied olefins is based on an open-loop SISO control of the production rate of reactors, as well as of the variables that define/affect the quality of the polymer. In other words, to each controlled variable corresponds a manipulated variable that is adjusted according to the judgment and experience of an individual. Thus, this configuration is not very effective, as it does not take into consideration the interaction among the several objectives of control. The changes in hydrogen flow rate, with the purpose of controlling the melt flow index of the polymer, cause changes in the production rate and in the density of the reaction medium. A further example of interaction is how the reaction medium density is affected when changes in catalyst flow rate are done in order to control the production rate. Besides the interaction problem, analyses involved in the measurement of the polymer properties result in delays which render more difficult any required corrective actions. In view of the considerations above, according to the present invention, and aiming to an effective operation of the plant at issue, it is required: • an on line system that infers the properties that are not continuously available; • a control system that can evaluate the interaction among different control objectives as well as the different process variables. Thus, it is possible to minimize problems related to the low frequency of laboratory analyses carried out on the variables that affect polymer quality, and to the interaction among control loops. TABLE 1 below lists the physical and chemical properties of the polymer that are inferred in each reactor. TABLE 1 TABLE 2 lists the physical and chemical properties of the polymer which are controlled in each reactor. This table indicates whenever a given variable is a controlled variable (CV), or a constraint (CCV). Besides those models, two further models were also built. One of such models was built to predict the production rate of the arrangement of loop reactors, while the other one was built to predict the ratio between the production rates of two consecutive loop reactors (split). TABLE 2 The input variables of each model composing the systems of the present invention are presented below: Melt flow index CMF This property is a controlled variable that is usually correlated to the molecular weight and polymer composition. In the state-of-the-art technique, the melt flow index is evaluated through periodical sampling, which causes time gaps between successive measurements, so as to mask what actually occurs in the plant. In case the MFI model is designed to infer or to control said property in loop reactors, the relevant parameters are: the hydrogen concentration in the feed stream(s) of the reactor(s); • the flow rate of catalyst fed to the reactor(s); • the flow rate of propylene fed to the reactor(s); • the flow rate of corήonomer(s) fed to the reactor(s); • the ratio of cocatalyst flow rate to electron donor flow rate; • the reactor(s) temperature; In case the MFI model is designed to infer said variable in gas-phase reactors, the relevant parameters are: • the hydrogen concentration in the feed stream(s) of the loop reactor(s) preceding said gas phase reactor; • flow of comonomer(s) fed to said gas-phase reactor; • the production rate of said gas phase reactor; • the production rates of the reactor(s) preceding said gas phase reactor; In case the MFI model is designed to control this variable in gas-phase reactors, the relevant parameters are: • the hydrogen concentration in the feed stream(s) of the loop reactor(s) preceding the gas phase reactor; • flow of comonomer(s) fed to gas-phase reactor; • the flow rate of catalyst fed to the reactor(s); • the reactor(s) temperature; Intrinsic viscosity dV) This property is only measured during operation of the gas- phase reactor(s). In case the IV model is designed to infer or to control said variable in gas-phase reactors, the relevant parameter is: • the ratio of hydrogen concentration to comonomer(s) concentration within the gas-phase reactor. Percentage of xylene-soluble matter In case the model ofthe percentage of xylene-soluble matter is designed to infer said property in loop reactors, the relevant parameters are: • the flow rate of comonomer(s) fed to the reactor(s); - the ratio of cocatalyst flow rate to the electron donor flow rate; • the reactor(s) temperature; • the density ofthe reaction medium within the reactor(s). In case the model ofthe percentage of xylene-soluble matter is designed to infer said variable in a gas-phase reactor, the relevant parameters are: • the percentage of xylene-soluble matter present in the polymer fed to said gas-phase reactor; ■ the percentage of comonomer(s) in the xylene-soluble matter present in the polymer; • the percentage of comonomer(s) present in the polymer. Percentagefs of comonomerfs) incorporated in the copolymer This variable is measured both in loop reactors and gas-phase reactors. In case the model(s) of percentage(s) of comonomer(s) incorporated in copolymer is(are) designed to infer and to control said variable in loop reactors, the relevant parameters are: • flow of comonomer(s) fed to the reactor; • production rate ofthe reactor; • flow of comonomer(s) fed to the preceding reactor(s); • production rate ofthe preceding reactor(s); In case the model(s) of percentage of comonomer(s) incorporated in the copolymer is(are) designed to infer said variable in gas- phase reactors, the relevant parameters are: • flow of comonomer(s) fed to gas-phase reactor; • production rate of gas phase reactor; • production rate of the reactor(s) preceding the gas phase reactor; In case the model(s) of percentage(s) of comonomer(s) incorporated in the copolymer is(are) designed to control this variable in gas-phase reactors, the relevant parameters are: • the hydrogen concentration in the feed stream(s) of the loop reactor(s) preceding the gas phase reactor; • flow of comonomer(s) fed to gas-phase reactor; • the flow rate of catalyst fed to the reactor(s); • the reactor(s) temperature; Percentage of comonomerfs incorporated in bipolymer It should be understood that bipolymer is the fraction of the polymer produced within gas-phase reactors. This property is only controlled during operation of the gas-phase reactor(s). In case the model(s) of the percentage of comonomer(s) incorporated in bipolymer is designed to infer or to control said variable, the parameter relevant to modeling is: • the ratio of the concentration of a given comonomer to the sum of the concentration of comonomer(s) and the concentration of propylene within the gas-phase reactor. Production rate and opening of the temperature control valve of loop reactors These variables are of utmost importance for the process economics. Lower variance resulting from the use of a more advanced control method such as that described and claimed in the present invention makes that average figures for these variables be increased relative to state-of-the-art control systems, provided the upper limits of the reactors heat exchange capacity, which are reflected in the temperature control valve opening, are not violated. The models for control of production rate and temperature control valve opening use the following variables as inputs to the mathematical models to predict the future behavior of said variables in loop reactors: • the ratio of flow rate of electron donor to flow rate of cocatalyst that comprise catalytic system; • the hydrogen concentration in feed stream(s) of loop reactor(s); • the flow rate of catalyst fed to loop reactor(s); • the flow rate of propylene fed to loop reactor(s); • the flow rate comonomer(s) fed to loop reactor(s); • the temperature ofthe loop reactor(s); Density of the reaction medium and power of the recirculation pump of loop reactors The models for density of reaction medium and recirculation pump power use the following variables as inputs to the mathematical models to predict the future behavior of said variables: • the flow rate of monomer fed to loop reactor; • the production rate ofthe loop reactor; • the density of reaction medium ofthe preceding loop reactor; Difference between the temperature and the bubble point of the liquid within loop reactors The model for the difference between the temperature and the bubble point of the liquid within loop reactors uses the following variables as inputs to the mathematical models to predict the future behavior of said variable: • the hydrogen concentration in feed stream(s) of loop reactor(s) and gas-phase reactor(s); • the flow rate of propylene fed to loop reactor(s); ■ the flow rate of comonomer(s) fed to loop reactor(s); • the temperature ofthe loop reactor; ■ the production rate ofthe loop reactor(s). Pressure of gas phase reactors The model for pressure of gas phase reactors uses the following variables as inputs to the mathematical models to predict the future behavior of said variable: • the flow rate of catalyst fed to the reactor(s); • the density of reaction medium ofthe loop reactor(s); • the flow rate comonomer(s) fed to the gas phase reactor; • the gas-phase reactor bed level; • the flow rate of gas returning from separation tower; Opening ofthe temperature control valve of gas phase reactors The model for the opening ofthe temperature control valve of gas phase reactors uses the following variables as inputs to the mathematical models to predict the future behavior of said variable: • the flow rate of catalyst fed to the reactor(s); • the density of reaction medium ofthe loop reactor(s); • the flow rate of comonomer(s) fed to the gas phase reactor; • the gas-phase reactor bed level; Ratio between the production rates of two consecutive loop reactors (split) The model for the ratio between the production rates of two consecutive loop reactors (split) uses the following variables as inputs to the mathematical models to predict the future behavior of said variable: • the ratio of flow rate of electron donor to flow rate of cocatalyst that comprise catalytic system; • the hydrogen concentration in feed stream(s) of loop reactor(s) and gas-phase reactor(s); • the flow rate of catalyst fed to loop reactor(s); • the flow rate of propylene fed to loop reactor(s); • the flow rate comonomer(s) fed to loop reactor(s); • the temperature of loop reactor(s); The several steps which make up the system of the present invention will be further discussed below. Collection of historical data of the plant of polypropylene and its copolymers with other olefins comprises storing the values of all variables measured during the daily operation of the plant, so as to build a large database. Selection of variables that will constitute the small database comprises the analysis of a great number of variables that are potentially relevant to the modeling of a given variable to be inferred or controlled. Not every variable is useful when building the mathematical models. The criteria to determine the usefulness of a variable is based on correlation analyses, as well as on the evaluation of models obtained from the use of said variable. Thus, the determination of the adequacy is a cyclic process, obtained conclusions being a feedback to the collection of historical data. Once the mathematical model is identified, many variables will not be useful to the modeling process, and will be discarded. Treatment of data ofthe small database comprises the removal of non-representative data due to failures of measurement instruments, or to abnormal operating conditions, as well as application of filters to remove noises which could be found in signals. The database is complemented with new variables calculated from theoretical concepts and mathematical correlations such as energy and mass balances, thermodynamic correlations, or any empirical correlations. These new variables, in turn, can be used as input for the mathematical models. Mathematical model identification comprises calculation of the model output based on the variables pertinent to the determination of the desired physical and chemical property using the values contained in the smaller database so as to determine the value of the model parameters which minimize the difference between the model prediction and the property actual values. In the case of neural nets, the identification step consists of the determination of the set of weights to be attributed to its connections, so as to obtain a neural mathematical model able to represent the system of interest. Checking the adequacy of the mathematical model obtained through the evaluation of a qualitative cause/effect relationship between each input variable and the estimated physical and chemical property. The purpose of the inference system of the variables of the process according to the present invention is to estimate the above cited polymer physical and chemical properties, and comprises a Distributed Digital Control System (DCS), a process computer connected to the DCS, an information system to store the historical data of the plant, and a system fed by data from laboratory analyses. FIGURE 4 is a schematic representation of the system for on line inference. When it is desired to infer a property according to the principles of this invention, the inputs of the mathematical models (which in this case are neural networks) are measured in real time by measurement instruments connected to the DCS. The values of said measurements are sent by the DCS to the computer program installed in the process computer responsible for the calculation of the mathematical model. The program that performs the calculations of the model also performs any preprocessing ofthe input signals if required. The model output is directed to the program that maintains the historical database of the process so that it can be used afterwards. The inference of the desired property is obtained from the multiplication of the output of the mathematical model by a factor, here called ratio, calculated from laboratory analyses. The purpose of the ratio is to correct possible errors on the estimate done by the model when compared with the results of laboratory analyses. This factor is calculated as the ratio of the value provided by the laboratory to the value estimated by the model at the very moment of the collection of the analyzed sample, and should be filtered before multiplication to avoid sudden changes in the inferred value, such changes being harmful to the control system that receives those signals. The ratio is calculated by the same program responsible for the maintenance of the historical database of the process, since, as mentioned before, this calculation takes into consideration the value estimated by the mathematical model at the moment ofthe collection of the analyzed sample, a piece of information only present in the historical database. A modification related to the preferred embodiment consists in the use of another parameter, called bias here, to perform the correction of the output of the mathematical model. The difference is that the bias must be added to the value of the output of the model, and is calculated as the difference between the value provided by the laboratory and the value estimated by the neural mathematical model at the moment of the collection of the analyzed sample. The bias must also be filtered before being added to the output of model. FIGURE 5 represents, according to the invention, how the use of an inference system together with a system of predictive control can be implemented. The values estimated for the controlled variables are used as inputs for the predictive controller. Then, the predictive controller solves a problem of non linear optimization to determine the values of the manipulated variables (set points of MV's), so that the controlled variables are maintained close to the desired values, without violating the limits of the plant. PID controllers receive the set points of the MV's calculated by the predictive controller, and are in charge of maintaining the variables within these values. The solution for the problems related to control comprises the use of the control matrix represented in FIGURE 6. The kind of the process variables shown in FIGURE 6 are described in TABLE 3. The maximization ofthe production rate and catalyst yield can be carried out through the selection of set points, which are arbitrarily increased, for the controlled variables production rate and density of the reaction medium. This procedure drives the process towards its physical constraints. Thus, the production rate is limited by the capacity of thermal exchange of the system, the control of said parameter being carried out through the opening of the valves for temperature adjustment, while the density of the reaction medium is limited by the power of the circulation pumps of the loop reactors. Thus, a balanced increase in the flow rates of catalyst and monomer is required so that the increase in production rate and density do not cause a violation of such limits as established for the opening of the valves for temperature adjustment, and for the power of the circulation pumps ofthe loop reactors. Benefits from the present invention could be noticed from the reduction in standard deviation of controller error calculated as the difference between the value of a process variable and its desired value. Oscillations in production rate and density were reduced thus allowing an increase in the average values of these variables. As a consequence of the increase in these average values, an increase in the production of polymer per kg of catalyst was observed. TABLE 3 EXAMPLE One of the main benefits obtained from the advanced control relates to the decrease in the variability of the product quality and in process variables such as the density of the reaction medium and production rate. If the maximum values usually observed for variables such as density ofthe reaction medium and production rate are maintained, the decrease in the variability is reflected in an increase of the average values for these variables, this being desirable under the economic point of view. FIGURE 7, where SP means set point and PV represents the process variables, shows an increase in the average production rate ofthe first loop reactor that resulted from the application of the control strategy, according to the present invention, in the polypropylene plant at issue. It is observed that the decrease in the variability of the production rate allows that the average of this variable is shifted 2.4% towards the maximum value usually obtained. FIGURE 8, on its turn, presents a decrease in the variability that is observed in the melt flow index of the polymer produced in the first loop reactor. Standard deviation of the error between the desired value for the property and obtained value was decreased from 7.2% to 3.2%, a 55% reduction of the standard deviation, equivalent to an 80% decrease in error variance.

Claims: CLAIMS WE CLAIM: 1. A system for on line inference of physical and chemical properties of polypropylene, its copolymers, and teφolymers produced in a plant with loop reactors, and, optionally, gas-phase reactor(s) with the aid of mathematical models, which comprises the steps of: • collecting data of all variables measured in the plant, so as to build a large historical database; ■ out of said large database, selecting a subset including historical data on the physical and chemical properties to be inferred, as well as potential input variables for the mathematical models so as to create a small database; • treating data of said small database; ■ from the treated data of said small database, selecting those variables that will be the input variables ofthe mathematical models; • using said input variables of the mathematical models, identifying said mathematical models off line; • applying the so-identified mathematical models to real-time measurements of said input variables to the model used to calculate the on line inference of desired physical and chemical properties. 2. A system according to claim 1, wherein, in the plant, the variables are measured by sensors and laboratory analyses. 3. A system according to claim 1, wherein the mathematical models may be empirical or rigorous or a combination of rigorous and empirical models. 4. A system according to claim 3, wherein the empirical mathematical models are neural networks. 5. A system according to claim 1, wherein data treatment of said small database comprises deleting non-representative data, as well as filtering noisy signals. 6. A system according to claim 1 , wherein data treatment of said small database comprises the generation of new variables from theoretical concepts and mathematical correlations that can be used as inputs for mathematical models. 7. A system according to claim 6, wherein the theoretical concepts and mathematical correlations comprise mass and energy balance data, thermodynamic correlations, or any empirical correlations. 8. A system according to claim 1, wherein the selection of input variables for the mathematical models comprises the correlation analysis among the potential input variables, and the variables that are sought to be inferred, as well as acquired knowledge on the polymerization process (know-how). 9. A system according to claim 1, wherein propylene is polymerized together with one or more olefin comonomers. 10. A system according to claim 9, wherein the weight percentage of propylene in the copolymer is at least 50. 1 1. A system according to claim 9, wherein the olefin comonomers are ethylene; 1-butene; 2-methyl propylene; 1-pentene; 1- hexene; 1-heptene; 1-octene; 1-decene; 3-methyl 1-butene; 4-methyl 1- pentene, and cyclopentene. 12. A system according to claim 1, wherein the physical and chemical property of polypropylene and its copolymers that is inferred is the melt flow index. 13. A system according to claim 12, wherein the mathematical model for inference of the melt flow index of polypropylene and its copolymers refers to the polymer produced in loop reactors and uses the following input variables: • the hydrogen concentration in the feed stream(s) of the reactor(s); • the flow rate of catalyst fed to the reactor(s); • the flow rate of propylene fed to the reactor(s); • the flow rate of comonomer(s) fed to the reactor(s); • the ratio of flow rate of cocatalyst to electron donor flow rate; • the reactor(s) temperature; 14. A system according to claim 12, wherein the mathematical model for inference of the melt flow index of polypropylene and its copolymers refers to the polymer produced in gas-phase reactors and uses the following input variables: ■ the hydrogen concentration in the feed stream(s) of the loop reactor(s) preceding the gas phase reactor; • flow of comonomer(s) fed to gas-phase reactor; • production rate of gas phase reactor; • production rate of the reactor(s) preceding the gas phase reactor; 15. A system according to claim 1, wherein the physical and chemical property of polypropylene and its copolymers that is inferred is the intrinsic viscosity 16. A system according to claim 15, wherein the mathematical model for inference of the intrinsic viscosity of polypropylene and its copolymers uses the ratio of hydrogen concentration to comonomer(s) concentration within the gas-phase reactor as input variable. 17. A system according to claim 1, wherein the physical and chemical property of polypropylene and its copolymers that is inferred is the percentage of xylene-soluble matter. 18. A system according to claim 17, wherein the mathematical model for inference of the percentage of xylene-soluble matter of polypropylene and its copolymers refers to the polymer produced in loop reactors and uses the following input variables: • the flow rate ofthe comonomer(s) fed to the reactor(s); • the ratio of flow rate of cocatalyst to flow rate of electron donor; • the reactor(s) temperature; • the density of reaction medium within the reactor(s). 19. A system according to claim 17, wherein the mathematical model for inference of the percentage of xylene-soluble matter for polypropylene and its copolymers refers to the polymer produced in gas- phase reactors and uses the following input variables: • the percentage of xylene-soluble matter of the polymer fed to the gas-phase reactor(s); the percentage of comonomer(s) in xylene-soluble matter of polymer; • the percentage of comonomer(s) in polymer. 20. A system according to claim 1, wherein the physical and chemical property of polypropylene and its copolymers that is inferred is(are) the percentage(s) ofthe comonomer(s) incoφorated in polymer. 21. A system according to claim 20, wherein the mathematical model(s) for inference of percentage(s) of comonomer(s) incoφorated in polymer refers(refer) to the polymer produced in loop reactors and uses(use) the following input variables: • the flow of comonomer(s) fed to the reactor; • the production rate ofthe reactor; • the flow of comonomer(s) fed to the preceding reactor(s); • the production rate ofthe preceding reactor(s); 22. A system according to claim 20, wherein the mathematical model(s) for inference of percentage(s) of comonomer(s) incoφorated in polymer refers(refer) to the polymer produced gas-phase reactors and uses(use) the following input variables: ■ flow of comonomer(s) fed to gas-phase reactor; • production rate of said gas phase reactor; • production rate of the reactor(s) preceding the gas phase reactor; 23. A system according to claim 1, wherein the physical and chemical property of polypropylene and its copolymers that is inferred is the percentage(s) of comonomer(s) incoφorated in bipolymer. 24. A system according to claim 23, wherein the mathematical model(s) for inference of percentage(s) of comonomer(s) incoφorated in bipolymer uses(use) the ratio of the concentration of a given comonomer to the sum of comonomer(s) concentration and propylene concentration within the gas-phase reactor(s) as input variable. 25. A system for the on line control of a plant of polypropylene and its copolymers produced in loop reactors, and, optionally, gas-phase reactor(s) with the aid of mathematical models, wherein said system is based on a control matrix, said matrix comprising controlled variables, constrained controlled variables, manipulated variables, and disturbance variables. 26. A system according to claim 25, wherein the mathematical models may be empirical or rigorous or a combination of rigorous and empirical models. 27. A system according to claim 26, wherein the empirical mathematical models are neural networks. 28. A system according to claim 25, wherein the controlled variables are: • the production rate ofthe arrangement of loop reactors; • the production rate of each loop reactor; • the ratio between the production rate of consecutive loop reactors; • the melt flow index of the polymer produced in the arrangement of loop reactors; • the melt flow index of the polymer produced in each loop reactor; • the melt flow index of the polymer produced in the arrangement of loop and gas-phase reactors; • the density ofthe reaction medium within loop reactors; ■ the percentage(s) of comonomer(s) incoφorated in polymer in the arrangement of loop reactors; • the percentage(s) of comonomer(s) incoφorated in polymer in each loop reactor ; • the percentage(s) of comonomer(s) incoφorated in the polymer produced in the arrangement of loop and gas-phase reactors; • the percentage(s) of comonomer(s) incoφorated in the fraction ofthe polymer produced in gas-phase reactors; • the intrinsic viscosity ofthe polymer. 29. A system according to claim 25, wherein the constrained controlled variables are: • the power of the pump that promotes the circulation of the reaction medium within each loop reactor; • the opening ofthe valve that controls the temperature of each loop reactor; • the opening ofthe valve that controls the temperature of each gas-phase reactor; • the pressure of each gas-phase reactor; ■ the difference between the temperature and the bubble point ofthe liquid within each loop reactor. 30. A system according to claim 25, wherein the manipulated variables are: • the ratio of the electron donor flow rate to cocatalyst flow rate; the hydrogen concentration in the feed stream of each loop reactor; • the flow rate ofthe catalyst fed to the reactor arrangement; • the flow rate of propylene fed to the loop reactor(s); • the flow rate of comonomer(s) fed to the loop reactor(s); • the ratio of hydrogen concentration to comonomer(s) concentration within gas-phase reactors; • the ratio of each comonomer concentration to the sum of propylene concentration and comonomer(s) concentration within each gas- phase reactor; • the flow rate of comonomer(s) fed to each gas-phase reactor. 31. A system according to claim 25, wherein the disturbance variables are: • the bed level in the gas-phase reactor(s); • the density ofthe reaction medium within the loop reactor(s); • the temperature of loop reactor(s); ■ the production rates of loop reactor(s); ■ the flow rate of the stream that returns from the comonomer/propylene separation tower; 32. A system according to claim 25, wherein the physical and chemical property of polypropylene and its copolymers that is controlled is the flow index. 33. A system according to claim 32, wherein the mathematical model for control of the flow index of polypropylene and its copolymers refers to the polymer produced in loop reactors and uses the following input variables: • the hydrogen concentration in the feed stream(s) of the reactor(s); • the flow rate of catalyst fed to the reactor(s); • the flow rate of propylene fed to the reactor(s); • the flow rate of comonomer(s) fed to the reactor(s); • the ratio of flow rate of cocatalyst to electron donor flow rate; • the temperature ofthe reactor(s); 34. A system according to claim 32, wherein the mathematical model for control of the flow index of polypropylene and its copolymers refers to the polymer produced in gas-phase reactors and uses the following input variables: • the hydrogen concentration in the feed stream(s) of the loop reactor(s) preceding the gas phase reactor; ■ flow of comonomer(s) fed to gas-phase reactor; • the flow rate of catalyst fed to the reactor(s); • the temperature ofthe reactor(s); 35. A system according to claim 25, wherein the physical and chemical property of polypropylene and its copolymers that is inferred is the intrinsic viscosity 36. A system according to claim 35, wherein the mathematical model for control of the intrinsic viscosity of polypropylene and its copolymers uses the ratio of hydrogen concentration to comonomer(s) concentration inside the gas-phase reactor as input variable. 37. A system according to claim 25, wherein the physical and chemical property of polypropylene and its copolymers that is controlled is(are) the percentage(s) ofthe comonomer(s) incoφorated in polymer. 38. A system according to claim 37, wherein the mathematical model(s) for control of percentage(s) of comonomer(s) incoφorated in polymer refers(refer) to the polymer produced in loop reactors and uses(use) the following input variables: • flow of comonomer(s) fed to the reactor; • production rate ofthe reactor; • flow of comonomer(s) fed to the preceding reactor(s); • production rate ofthe preceding reactor(s); 39. A system according to claim 37, wherein the mathematical model(s) for control of percentage(s) of comonomer(s) incoφorated in polymer refers(refer) to the polymer produced in gas-phase reactors and uses(use) the following input variables: • the hydrogen concentration in the feed stream(s) of the loop reactor(s) preceding the gas phase reactor; • flow of comonomer(s) fed to gas-phase reactor; • the flow rate of catalyst fed to the reactor(s); • the temperature ofthe reactor(s); 40. A system according to claim 25, wherein the physical and chemical property of polypropylene and its copolymers that is controlled is(are) the percentage(s) of comonomer(s) incoφorated in bipolymer. 41. A system according to claim 40, wherein the mathematical model(s) for control of percentage(s) of comonomer(s) incoφorated in bipolymer uses(use) the ratio of the concentration of a given comonomer to the sum of comonomer(s) concentration and propylene concentration inside a gas-phase reactor as input variable. 42. A system according to claim 25, wherein the variable to be controlled is the production rate of loop reactors. 43. A system according to claim 42, wherein the mathematical model for control ofthe production rate of loop reactors uses the following as input variables: • the ratio of electron donor flow rate to cocatalyst flow rate; • the hydrogen concentration in feed stream(s) of loop reactor(s); • the flow rate of catalyst fed to loop reactor(s); ■ the flow rate of propylene fed to loop reactor(s); • the flow rate comonomer(s) fed to loop reactor(s); • the temperature ofthe loop reactor(s); 44. A system according to claim 25, wherein the variable to be controlled is the opening ofthe temperature control valve of loop reactors. 45. A system according to claim 44, wherein the mathematical model for control of the opening of the temperature control valve of loop reactors uses the following as input variables: ■ the ratio of flow rate of electron donor to flow rate of cocatalyst that comprise catalytic system; ■ the hydrogen concentration in feed stream(s) of loop reactor(s) and gas-phase reactor(s); • the flow rate of catalyst fed to loop reactor(s); • the flow rate of propylene fed to loop reactor(s); • the flow rate comonomer(s) fed to loop reactor(s); • the temperature of loop reactor(s); 46. A system according to claim 25, wherein the variable to be controlled is the density of reaction medium within loop reactors. 47. A system according to claim 46, wherein the mathematical model for control of the density of reaction medium within loop reactors uses the following as input variables: • the flow rate of monomer fed to loop reactor; • the production rate ofthe loop reactor; • the density of reaction medium ofthe preceding loop reactor; 48. A system according to claim 25, wherein the variable to be controlled is the power ofthe recirculation pump of loop reactors. 49. A system according to claim 48, wherein the mathematical model for control of the power of the recirculation pump of loop reactors uses the following as input variables: ■ the flow rate of monomer fed to loop reactor; • the production rate ofthe loop reactor; • the density of reaction medium ofthe preceding loop reactor; 50. A system according to claim 25, wherein the variable to be controlled is the difference between the temperature and the bubble point ofthe liquid within loop reactors. 51. A system according to claim 50, wherein the mathematical model for control ofthe difference between the temperature and the bubble point of the liquid within loop reactors uses the following as input variables: • the hydrogen concentration in feed stream(s) of loop reactor(s) and gas-phase reactor(s); • the flow rate of propylene fed to loop reactor(s); • the flow rate comonomer(s) fed to loop reactor(s); • the temperature ofthe loop reactor; • the production rate ofthe loop reactor(s). 52. A system according to claim 25, wherein the variable to be controlled is the pressure of gas phase reactors. 53. A system according to claim 52, wherein the mathematical model for control of the pressure of gas phase reactors uses the following as input variables: • the flow rate of catalyst fed to the reactor(s); • the density of reaction medium ofthe loop reactor(s); • the flow rate comonomer(s) fed to the gas phase reactor; • the gas-phase reactor bed level; • the flow rate of gas returning from separation tower; 54. A system according to claim 25, wherein the variable to be controlled is the opening of the temperature control valve of gas phase reactors. 55. A system according to claim 54, wherein the mathematical model for control of the opening of the temperature control valve of gas phase reactors uses the following as input variables: • the flow rate of catalyst fed to the reactor(s); • the density of reaction medium ofthe loop reactor(s); • the flow rate comonomer(s) fed to the gas phase reactor; • the gas-phase reactor bed level; 56. A system according to claim 25, wherein the variable to be controlled is the ratio between the production rates of two consecutive loop reactors (split). 57. A system according to claim 56, wherein the mathematical model for control of the ratio between the production rates of two consecutive loop reactors uses the following as input variables: • the ratio of flow rate of electron donor to flow rate of cocatalyst that comprise catalytic system; • the hydrogen concentration in feed stream(s) of loop reactor(s) and gas-phase reactor(s); • the flow rate of catalyst fed to loop reactor(s); ■ the flow rate of propylene fed to loop reactor(s); ■ the flow rate comonomer(s) fed to loop reactor(s); • the temperature of loop reactor(s);

Inventors: AXELRUD CELSO, DEMORO ESDRAS PIRAGUACY, OLIVEIRA ARTUR TOLEDO MARTINS DE, AXELRUD, CELSO, DE OLIVEIRA, ARTUR, TOLEDO, MARTINS, DEMORO, ESDRAS, PIRAGUACY
Applicants: AXELRUD CELSO, DEMORO ESDRAS PIRAGUACY, OLIVEIRA ARTUR TOLEDO MARTINS DE, OPP PETROQUIMICA S A, AXELRUD, CELSO, DE OLIVEIRA, ARTUR, TOLEDO, MARTINS, DEMORO, ESDRAS, PIRAGUACY, OPP PETROQUIMICA S.A.
Application year: 1999
Main classification: G05B17
IPC codes: B01J19, C08F10, C08F110, G01N33, G05B17, G05B23, G06F19