SYSTEM AND METHOD FOR IMPROVED FLOW DATA RECONCILIATION _________

Abstract: The present invention is a system (100) and method for improved processing of flow data signals to provide output flow data signals which better approximate the true flows being measured. The general approach on which the invention is based is commonly known as flow data reconciliation. A dynamic simulation (110) is run in parallel with a flow data reconciliation, and model predicted flows are used to complete balances wherever sensor measurements are unavailable. Model predicted compositions are also used to allow the computation of stream enthalpies so that enthalpy balances can be used in the reconciliation. The use of model predicted values of changes in inventory allow dynamic material balances to be used, making the reconciliation much more effective for process with slow dynamic behavior. Weighting factor on sensors (106), computed based on the current value of the error of the sensor (106), reduce the impact of sensors (106) with high errors. A multiplicative correction factor distributes corrections more evenly among streams with large and small flows. Statistical behavior of historical performance of sensors (106) is used to reduce the adverse impact of rapid shifts in sensor accuracy and to detect unconverged reconciliation solutions.

Description: SYSTEM AND METHOD FOR IMPROVED FLOW DATA RECONCILIATION CROSS REFERENCE TO COPENDING PATENT APPLICATION The present application is related to the following concurrently filed patent application, assigned to the Assignee of the present application: U.S. Patent Application (Attorney Docket 1100.1220000) by David V- Cardner, entitled "SYSTEM AND METHOD FOR IMPROVING MODEL PRODUCT PROPERTY ESTIMATES", which is hereby incorporated by reference in its entirety. BACKGROUND OF THE INVENTION 1. Field of the Invention The invention relates to flow data reconciliation, the estimation of true flow values from measured flow values, and to the application of flow data reconciliation in process control. 2. Discussion of Related Art Anytime sensors are used to measure flow in a chemical process, the resulting sensor measurements are susceptible to errors. Flow meters vary in their accuracy depending on the material in the stream, the conditions of the stream, and the sensing technology used in the sensor; and they vary in reproducibility as their performance deteriorates due to fouling, wear, or failure. Plants can improve or reduce sensor performance through frequency of maintenance. Thus, flow sensors have both fixed and variable performance factors. Anytime that flow measurements are used for critical control or monitoring tasks in a chemical process, it is desirable to have the true flow value rather than the measured value, with its attendant uncertainty. This is especially true for two key uses of data: as input to process models and as input to closed loop control schemes. In these uses, inaccurate data can cause the generation of erroneous predictions or control moves, both of which could adversely affect performance of the process. Although a process operator cannot have perfect knowledge of flow rates, there are techniques for deriving more accurate flow values. In many cases, all flows around a unit operation will be measured. Applying a fundamental engineering principle of mass conservation, a steady state material balance can be applied to a unit operation, which states that the flows in minus the flows out must equal zero. When all the flows around a unit are measured, the measured flow values can be substituted into this material balance equation. However, due to inaccuracies in the measured flow values, the flows will not equate to zero. The difference between the actual sum and correct value of zero is the error (or residual) around the unit Each measured flow value can then be adjusted to reduce the error in the balance equation to zero. The corrected flow values are estimates of the true flow values. Conventional techniques teach that for systems with random measurement errors it is preferred to adjust the flows by minimizing the sum of the squared differences between each measured and estimated flow value. When a flow is unmeasured around a unit, there is not enough information to calculate the residual error, and thus no basis on which to estimate true flows. However, in some cases, the missing value is an output stream from one unit which is also an input stream to another unit If a balance equation is written around both units, the flow value disappears from the balance. When the common flow value is the only missing flow for both units, a residual error can be calculated around the two units together, and thus new flow estimates can be made. The missing common flow value can then be estimated by completing the balances with the new flow estimates. Conventional techniques teach that this process can be expanded, making balances around all combinations of units. By this process it may be possible to make improved estimates for all flows in the process, even if insufficient flows are measured to make single-unit balances to correct all the flows. In addition to the use of conservation of mass to develop balances around each unit, the art suggests that the principle of conservation of energy can also be applied. However, to compute energy flows around a unit, the enthalpy of each stream must be computed, plus direct energy flows (such as steam heating) must be estimated. In multicomponent streams, which are typical in chemical processes, stream enthalpy can only be computed if both the temperature and the composition of the stream are known. Since composition information is usually known for only some components in some streams and temperatures may be missing, this method is not generally useful in commercial practice. Some conventional techniques teach that steady state balances be constructed for computing the errors. This means, however, that the technique cannot be applied to a process in which flows or inventories are changing. Since many process have dynamic flows or inventories, these techniques fail to address many important process situations. Morover, some conventional technqiques teach that each error derives from the magnitude of the corrected versus measured value. This approach has the limitation that absolute errors in large flows tend to cause relatively large corrections to small flows. Although in principle a complete set of reconciled flows could be obtained by conventional methods, a plant does not typically measure enough flows to build the balances around all the units in the plant The absence of important flow measurement can severely limit the ability to derive corrected estimates. To compensate for this, additional balances can be constructed by combining two or more units into a single material balance. This allows additional balances to be constructed, providing more opportunity to find balances in which all flows are measured. However, these techniques are very difficult to implement automatically. It requires a complicated search procedure to identify which multi-unit balances will improve the resulting estimates. It is thus very difficult to apply this technique where full sensor data is lacking. This approach is particularly difficult to implement in real world situations where analyzers both fail and are returned to service in a random manner across the process. Furthermore, some conventional techniques teach summing the weighted squared errors of the corrected versus the raw flows, and minimizing the sum, where the weighs are selected based on prior knowledge of sensor performance and flow importance. However, this approach uses fixed weights, even though some flow sensors can and do undergo significant changes in accuracy. Thus a sensor whose behavior changes due to failure or deterioration could drive many related sensors to have incorrect adjusted estimates. Other techniques teach a simple minimizing of the global error. However, there may be cases where the result of the process used to find the minimum error leads to corrections that are many times different than the expected corrections for one or more sensors. These techniques are incapable of detecting such problems. Finally, some conventional techniques teach that stream compositions must be known to compute stream enthalpies. Since compositions are not generally known except for some components in some streams, conventional techniques are not practical for computing stream enthalpies. SUMMARY OF THE INVENTION The present invention is a system and method for improved processing of flow data signals to provide output flow data signals which better approximate the true flows being measured. The general approach on which the invention is based is commonly known as flow data reconciliation. Sensors are used to generate raw flow measurement signals from a chemical process. The signals are indicative of the actual flow rates in the process. The signals are converted to computer form, and input to the flow data reconciliation system, along with measured properties or compositions. A first principles simulation of the process is carried out on the computer, using the pre-defined estimates of flows as starting values, to define as many process flows as possible. Unmeasured process flows are predicted by the simulation, resulting in predicted flow measurement signals. The predicted flow measurement signals are combined with the raw flow measurement signals to form a complete set of flow measurement signals. The simulation may also produce predicted product property signals. The complete set of flow measurement signals is processed by a flow data reconciliation module running on the computer, using steady state flow data reconciliation techniques. In this process, steady state total material balances are constructed around each process unit. Whenever the residual error around a unit is non-zero, each of the flows around the unit is adjusted by a correction factor to generate corrected flow signals. An error is computed between each raw/predicted flow signal and the corrected flow signal. The error for each raw flow measurement signal can be assigned an a priori weight, which can increase or decrease the impact of error from that signal. The correction factors are chosen so that the balance equation is satisfied, and so that the sum of the squared errors is minimized. This basic reconciliation process of the invention is improved in several ways. First, a dynamic simulation can be used, and the accumulation term computed by the simulation for each process unit can be used to perform dynamic flow data reconciliation. In this process, the residual error is the difference between the sum of all flows around a unit and the accumulation within the unit In some cases the rate of accumulation can be measured directly. This allows accurate corrected flow signals to be generated even when the process is not in a steady state. Second, the dynamic simulation can be used to predict the concentration (such as a mole fraction or weight fraction) of individual components in each stream. By combining this information with signals indicative of temperatures and direct heat inputs/outputs from a unit, an enthalpy (or energy content) can be computed for each stream. This allows an enthalpy reconciliation to be performed. In this process, errors in enthalpy balances rather than mass balances are used to improve the corrected flow signal estimates. Third, the distribution of residual error among the various flows around a unit is improved by generating corrected flow, signals by multiplying each raw or predicted flow measurement by a correction factor rather than adding a correction term, then minimizing the weighted squared differences of these from 1.0. This avoids large adjustments in flows which are small relative to other flows around a unit Fourth, sensor failure and degradation can be better handled by two improvements. The weighting factor on the error of a signal can be reduced when that error is large. This tends to reduce the impact of sensors with large errors. Also, a statistical history of the error of each signal can be used to improve the determination of when a raw flow measurement signal has a high error. This reduces the impact of dramatic sensor failures, and also identifies situations where the reconciled flow signals are erroneous due to computational errors in the reconciliation process. The utility of the invention is further improved by combining the signal processing of the invention with other known or novel functions to produce an overall process control system. Regulatory controllers (PTLF controllers), product property controllers, statistical filters, simulation models, alignment and calibration functions, and comparator functions can be combined to form a robust and effective embodiment of the invention for controlling properties such as composition in a chemical process. The various aspects of the invention can be used individually or in combination. BRIEF DESCRIPTION OF THE DRAWINGS The above and further advantages of this invention may be better understood by referring to the following description taken in conjunction with the accompanying drawings, in which: Figure 1 shows a schematic diagram of the present invention; Figure 2 shows a schematic diagram of the present invention implemented as a computer-based flow data reconciliation system and product property control system; Figure 3 Shows a schematic diagram of the invention implemented as model alignment and product property control system; Figure 4 Shows a schematic diagram of the invention implemented as computer based model calibration and product property correction and product property control system; and Figure 5 shows a schematic diagram of the preferred embodiment of the invention implemented as a computer based process stream composition and product property control system. DETAILED DESCRIPTION OF THE INVENTION The preferred embodiment of the present invention is a computer- based system referred to here as "the process autopilot". The process autopilot is a system for performing closed-loop composition and product property control in any manufacturing process handling flowable materials such as liquids, vapors, gases and flowable solids. Examples are chemical manufacturing process plants, refineries, or polymer processing plants. For closed loop stream composition or product property control to be successful, each stream component or product property to be controlled must be estimated frequently and accurately, and with a high level of reliability. Various aspects of the present invention are used to process raw flow measurement signals to produce corrected flow measurement signals, which can then be used as inputs to simulations which provide the needed frequent, accurate, reliable estimates of stream compositions or product properties. The preferred embodiment is named process autopilot because it performs closed-loop stream composition or product property control in a plant with one or more unit operations. This fully automates the process so that it no longer requires operator intervention during normal operation. The preferred embodiment is best employed in control of multi- unit continuous process operations. However, the system can be applied to batch operations, which are continuously varying operations that are usually performed in a single process unit As employed herein, "product properties" are measured or computed physical or chemical characteristics of any process material. Properties might be inherent in the material or indicative of their end-use performance. The term "product properties" generally means stream composition when applied to liquid or vapor streams in the chemical industry. In petroleum refining operations where the number and types of components may become too numerous to define, product properties may be characterized by property measurements such as: cut-point temperature, flash point, and Reid vapor pressure. Product properties for solid streams usually refer to properties such as: melt index, molecular weight distribution, and particle size distribution. Product properties can also be end-use characteristics of a product, such as the sheen produced when a pigment is used to make a paint, or the knocking behavior of gasoline. The equipment regulatory variables - pressure, temperature, level, and flow (designated PTLF variables) - are controlled to maintain process equipment operating at a designated stationary state, regardless of disturbances to the inputs. Certain PTLF variables may also be adjusted by the operator to move the process equipment to another stationary state (e.g., to increase production) while other PTLF variables are held constant in order to maintain the desired stream composition and/or product properties. Figure 1 illustrates a flow data reconciliation system 100. The system 100 operates on a process 104. The present invention is applicable to any process 104, including chemical processes, refineries, and polymer processes. Often the material in the process can be treated as a fluid, including viscous polymers and finely divided solids, which are moved within the process in streams. A process is normally made up of more than one unit, where each unit carries out some specific processing function, such as reaction, distillation, or heat exchange. Units interconnect via streams. In the presently preferred embodiment of the present invention, the process simulation models 110 and the flow data reconciliation module 108 are implemented on a digital computer. In one interpretation, block 102 may be viewed as the computer platform on which these two modules are implemented. In addition, as shown in Figure 2, the data flow reconciliation system 102 is used in combination with a product property control system 270. In the presently preferred embodiment, the product property control system includes the product property controllers 240 and PTLF controllers 2S0, both of which are implemented in a distributed control system (DCS). The product property controllers 240 and PTLF controller 250 are both well known inethe art Block 270 can be interpreted as showing the DCS as a platform for these modules. Figure 5 shows a schematic diagram of the preferred embodiment of the present invention (process autopilot) as a computer-based stream composition and product property control system. In the. process autopilot, all of the computer-based modules 402, 4? 410, 414, 406, 316, 110, 130, 108, 308 and 312 are implemented on a common computer, a PDP-11. It should be understood, however, that any number and nature of computer platforms, and any distribution of modules among them is contemplated by the present invention, as long as the signals in the invention (odd numbered in the drawings) can be communicated via analog or digital means between the modules. I. Flow Data Reconciliation In the following discussion, the terms model, simulation, and simulation model may be used. These terms have the same meaning, and refer to computer-based predictors 110 which utilize engineering and scientific principles to compute desired values from known values. These models are also referred to as first principle models. Referring again to Figure 1, the process simulation models 110 compute predicted product property signals 117 that can be used to provide closed-loop product property or composition control of the process, as well as for other procedures within the process autopilot The process simulation models 110 also provide estimates of PTLF variables 113 used by other process autopilot functions. The predicted flow values 113 are used in completing balances in the flow data reconciliation procedure. Input to the simulation models 110 consists of reconciled (corrected) flows, temperatures, pressures, and levels 115 and input property values, which can be measured 119 or estimated (not shown), and may include reconciled composition values 115. In addition, to initialize the process autopilot, a priori estimated flow measurement values are used to generate initial simulation results. This is the first step used to initialize the process autopilot The preferred method of simulation is with a dynamic, first- principles model for each unit operation in the process. First principles simulation of chemical processes is well known in the art Any method of modeling, either analytic or numerical, can be used. Numerically solved models are typically used because they can model a much wider range of process units. The simulations should use a computation cycle much shorter than the process response time. When the process dynamics are faster than the fastest possible computation cycle, the process appears to be at steady- state and a steady-state first-principles model is adequate. In this case all inventoiy rate-of-change factors are zero. The more general case of dynamic simulation will be discussed next The models 110 can use three types of input First, they must be provided with the most recent measurement signals of pressure, temperature, level, and flow 115. The flow values are preferably reconciled (corrected) signals 115 (except for initialization, when raw signals 109 or estimates are used), and preferably the pressure and temperature values should be limit checked. The second type of input is feed composition or feed properties 119. These can be provided from lab or analyzer 112, or computed sources. Third, the simulation models 110 may optionally accept adjustments to model calibration and to alignment factors. Model calibration adjustments change model parameters, such as tray efficiency for distillation columns or catalyst activity for reactors, or biases model flows. This is done to make the model more accurately reproduce measured product property or stream composition values. Continual recalibration of the models compensates for changes in the process behavior relative to the model behavior. Alignment keeps the models in synchronization with the process by matching equipment PTLF regulatory variables. The models are run on a fixed frequency, in which each run simulates a given time period of operation of the plant For example, the models may be run every five minutes to simulate five minutes of plant operation. A process may have several process steps and be represented by several models, which can be run independently of each other with respect to flows. Each model is independently initialized with respect to flows, temperatures and pressures for each computation cycle using reconciled values from ihe process. Since flow values are not transferred between models, they can be manipulated by the calibration and alignment function without affecting downstream models. Usually, only composition values are transferred between the models. This method of decoupling the models greatly enhances the ability to accurately and reliably track product property values across a multi-unit process. The validity check module 130 first inputs the raw sensor measurement signals 109 (which include flow measurement signals and signals from other types of sensors such as pressure). These signals aie screened against high and low limits to detect gross errors, and default values are substituted should they occur. This produces validated raw sensor measurement signals 111 (which include validated raw flow measurement signals), which are used by the flow data reconciliation module 108. The flow data reconciliation module 108 provides improved flow measurement signals for use by the process simulation models 110. This is important since an inaccurate flow measurement signal value could cause substantial inaccuracy in the computed product properties, resulting in improper control, production of off-specification product, and resulting ingredients and energy waste, and lost production opportunity. A number of aspects of the flow data reconciliation are configured in accordance with the known process. The balances used in the reconciliation are configured to match the process connections for each process unit Also, an error weight may be entered for each flow measurement signal. A low weight indicates an expectation that a signal will be inaccurate. A high weight indicates an expectation that a signal will be accurate. The flow data reconciliation module 108 uses validated raw flow measurement signals 111, simulated (predicted) pressure, temperature, level, and flow measurement signals 113 and rate of change of inventory signals 113 from the dynamic simulation 110. The reconciliation module 108 uses dynamic mass, energy, and/or component balances around each unit (of process 104) that describe the overall energy and material balances for that unit For those process steps where the rate of change of equipment inventory is significant, validated measurements 111 or estimates 113 of the rate-of-change of inventory are incorporated into the energy and material balance relations. The flow data reconciliation module 108 first inputs predicted pressure, temperature, level, and flow measurement signals 113 and optionally predicted product property signals 117 from the process simulation models 110 sufficient to satisfy all of the data requirements of the reconciliation method. Then, all available validated raw flow measurement signals 111 are used to replace the corresponding predicted signals. Unmeasured flows and other unmeasured variables are thus estimated directly by the process simulation model 110. A calculation to determine the enthalpy of each stream is then performed. This produces a complete flow data reconciliation model data set, which is used by the flow data reconciliation module 108. Moreover, the flow data reconciliation module 108 corrects flow values using the mass and energy balances around each unit, using the following steps. First, a starting correction factor (close to 1.0) is chosen for each flow signal and each inventory accumulation term, such that the sum of all flows and inventory accumulation terms times their correction factors equals zero (i.e., the balance is satisfied around the unit). The correction factor for each flow signal is subtracted from 1.0 to produce an error signal, squared, and then multiplied by two weights; a weight reflecting the degree of confidence in the measurement, and a weight that reflects the degree of correction required for the flow. The weighted squared errors are summed to give a sum of weighted error values. The mass and energy balances constitute the equality constraints, and the sum of weighted error values is the objective function to be minimized by adjusting the correction factors for each flow signal. The solution may be obtained using quadratic or non-linear programming algorithms or other techniques described in the art If component balances with corrections to component concentrations are included as additional balance equations, and/or correction factors are applied to stream enthalpies, then the equality constraints become non-linear and non-linear programming methods or other spedal methods in the art must be utilized to minimize the sum of weighted error values, while continuing to satisfy the balance relationships (equality constraints) around each unit A non-linear programming algorithm is utilized in the preferred embodiment to minimize the sum of weighted error values. When a converged solution is obtained, each correction factor is compared with its previous moving average correction factor using the below described procedure for statistical filter 404. If this comparison detects a major change in one or more correction factors, then new variable weights for those terms in the sum of weighted error values are generated and the reconciliation process repeated until all correction factors fall within a prescribed range of the moving average correction factor (normally 5 deviations), or the maximum number of iterations is reached, or the time limit is exceeded. When the final converged solution is obtained, then new estimates of the moving average correction factors are generated by combining the old moving average correction factor with the new value in an exponentially weighted filter. If a flow signal is found to have an extreme correction factor, meaning the flow is being adjusted significantly, the variable weight on the weighted squared error term is reduced. This effectively discounts the sensor signal, and causes that correction factor to absorb more correction while other correction factors are changed little. On completion of the reconciliation, new moving average values of the correction factor on each flow signal and new moving average variances are computed and saved, along with the updated values for the variable weight applied to the terms in the sum of weighted error values. This procedure is performed for each computational cycle of the process autopilot, which may be several times per hour. The frequency depends on the time constant of property dynamics of the process, and must be frequent relative to the time constant For example, qne estimate of the time constant is how long it takes the product property, in response to a step change in a regulatory setpoint (which affects that property), to move halfway to its new steady state value. The flow data reconciliation module 108 produces a corrected set of pressure, temperature, and flow signals 115 as input to other process autopilot procedures. If composition reconciliation is also performed, reconciled input compositions 115 are also available. Referring to Figure 2, PTLF equipment regulatory controllers 250 can operate as single-loop, cascade or other types of controllers as known in the art Fluid handling equipment can operate over relatively wide ranges of pressure, temperature, level, and flow. These variables, as measured by sensor signals 109, must be carefully controlled to provide stable, predictable operation of the process equipment at the correct operating rate. Normally, automated control of these variables is easily achieved, primarily because accurate sensors 106 are available that provide continuous signals to automatic controllers 250. The PTLF controllers 250 regulate pressure, temperature, level, and flow to provide stable, predictable operation of the process equipment at their setpoint values 245, 247. Some PTLF controllers receive controller setpoint adjustments 247 only from the operator, while others can also receive their setpoint adjustments 245 from product property controllers. The PTLF controllers 250 send output signals to actuators 260, such as valve loadings to control valves, to maintain the PTLF variables as close to setpoint as possible. Examples of PTLF variables that do not significantly influence product properties are condensate pot level, storage tank pressure, and condenser exit water temperature. Examples of PTLF variables that often affect product properties are fixed- bed reactor inlet temperature, calandria steam flow on a distillation column, and reflux ratio on a distillation column. In the preferred embodiment, PTLF controllers 250 are implemented using a distributed control system. However, any implementation, such as analog pneumatic controllers, could be used. The setpoints 245 for some PTLF controllers 250 may be adjusted by an operator after examination of analytical results 119 from a process analyzer or laboratory analysis 112. Alternatively, automatic product property control can be implemented by providing measured product property signals 119 or predicted product property signals 117 or 401 to product property controllers 240 that automatically adjust the setpoints 245 of the PTLF controllers 250. Effective automation requires product property measurements 119 or estimates 117 or 401, which are accurate, reliable, and produced frequently relative to the speed of process changes. The operator adjusts, as shown in reference number 235, the product property setpoints to match the product property specification of the product that is being produced. In the preferred embodiment, product property controllers 240 are implemented using a distributed control system. However, any implementation, such as analog pneumatic controllers, or a computer- based supervisory control system could be used. π. Calibration and Alignment System and Method One feature of the present invention is a system and method for improving the accuracy and precision of stream composition and/or product property estimates produced by process simulation models run in parallel with the process 104. The general approach of .the present invention is to use an assortment of procedures singly, or in combination, on each model and each product stream where appropriate. One procedure keeps the model in alignment with the process by minimizing the errors between computed and measured values of key regulatory variables. A second procedure calibrates the model by adjusting the model parameters, or input flow or composition signals, to minimize the difference between measured and calculated output quality signals. A third procedure adjusts the model output property signals 117 directly. A fourth procedure adjusts alignment controllers 316. The objective of all four procedures is to minimize the error between measured and computed stream compositions and/or product properties. A. Calibration Figure 4 shows a block diagram of a computer-based model calibration and product property correction and product property control system. The calibration -system 408 is based on comparing computed and measured stream compositions and/or product properties. Adjustments based on the difference are then applied to the models 110, the model estimates 117, and/or the alignment system 318. Product property signals 119 are made by on-line or in-line analyzers 112, or by extracting a sample from the process and analyzing it in a laboratory (also shown as 112). Both methods produce a product property measurement signal 119. In-line analyses produce continuous measurements. Analyses from on-line sample-based analyzers typically vary in frequency from every few minutes to several hours. Laboratoiy analyses typically vary in frequency from several hours to once each day. The product property signals 119 are indicative of the actual stream compositions and product properties of process 104. Any analog product property signals 119 are converted to digital signals and input with other digital product property signals to the calibration and product property corrector system 408. A first principles simulation 110 of the process is carried out on the computer using the reconciled flows and validated pressure, temperature, and level measurements 115, and feed composition values in signal 119. The simulation 110 produces stream composition and/or product property estimates 117 that are introduced into the product property comparison 410. In processes with several process operations, not all streams exiting each process operation are analyzed. Usually final product streams (main product and waste streams) are analyzed. It is often possible, using simple material balance calculations, to estimate the concentration of a single component for one of the interior streams based on the final product stream analysis. The product property transform module 414 performs this function. The product property transform module 414 utilizes product property measurement signals 119, with reconciled flows 115 and predicted product property signals 117 to produce artificial analyses 407 for streams that normally are not analyzed using component and flow ratios and simple material balance relationships. These are processed along with the normal analyses (product property signals 119) by product property comparison module 410 to produce error estimates 405, which can then be used to automatically recalibrate the appropriate model. Use of artificial analysis greatly expands the capability to maintain the simulation models in calibration. The product property comparison module 410 computes instantaneous product property errors 405 between measured analytical signals 119 and/or artificial analytical signals 407 and the predicted product property signals 117 computed by the simulation models 110 at the time of the sample. The instantaneous product property errors 405 can be expressed as arithmetic differences, as arithmetic differences divided by the range of the measured variables, as arithmetic differences divided by the computed or measured value, as the natural logarithm of the computed and measured values, or any combination of the above. The natural log of the ratio is preferred for low concentrations of stream components so that a symmetry in the error is established around the computed value. Either relative error or the log of the ratio can be utilized for higher concentrations. Proper selection of the method for expressing the error can many times linearize the relationships between the error and the response to the corrective action; thus making corrective actions more accurate over wider ranges of error. The best method often eliminates the possibility of producing unrealistic product property estimates (e.g., negative mole fractions). The statistical filter module 404 calculates a corresponding product property moving average error 403 each time an instantaneous product property error 405 is calculated. Each instantaneous product property error 405 is treated independently of the others. The statistical filter module 404 compares the instantaneous product property error signal 405 with the previously determined value for the moving average error signal 403. The difference between the error 405 and the moving average error 403, when squared, is the instantaneous variance. The deviation of the error from the average error, measured in multiples of the average variance, is computed and applied to the normal frequency distribution to determine weighting factors for an exponentially weighted filter. The filter is used to update the moving average error. The deviation can be applied to the normal frequency distribution to produce an estimate of the probability of occurrence of a larger deviation. If the deviation is large, then the probability of occurrence of an even larger value is low. If the deviation is very small, then the probability of occurrence of a larger deviation may approach 1.0, and the signal is considered to be very accurate. The probability of occurrence of a larger value is used as a multiplier of the instantaneous product property error 405, and the probability of a smaller value is used as the multiplier of the previous moving average error 403. The two products are added to produce the new moving average error 403. The statistical filter module 404 calculates a moving average variance each time a new moving average error is calculated. The new average variance is obtained by combining the current value of variance with the previous moving average value of variance in an exponentially weighted filter where the fixed weights are selected by the user. On occasion, the property signals 119 or 407 may shift quickly to new and substantially different values due to instrument or process changes. This would cause the error signal 405 to also change quickly. Under normal circumstances, the sudden large difference between the error signal 405 and the previous value of the moving average error 403 would cause the new error signal to be treated as being erroneous and would be given little weight in generating a new moving average error 403. The instantaneous error signals 405 would continue to be heavily discounted until the moving average variance had increased as a results of the repeated large differences. This would normally require many cycles before the moving average error 403 reflected the true error. During this period, the corrections made by the product property corrector 402 would be inaccurate. This problem is overcome by using a modification of the cusum (cumulative summation) statistical procedure that involves counting the number of consecutive times the deviation exceeds 1.0 (or other threshold value) in the same direction from the moving average error. This count is used as a divisor of the deviation to produce an artifiάal deviation that is applied to the normal frequency distribution to produce probability factors that weight the current value of error 405 higher and the previous moving average error 403 lower in the filter. This method moves the moving average error 403 to the true value of the error in much fewer steps, while the moving average variance remains relatively undistorted. Occasionally, the normal variability of process analyzer or laboratory analyses is large and produces a large moving average variance. When this occurs, filtering is inadequate to prevent significant changes in the moving average error 403 with each new analysis, which when used to correct product property estimates causes erratic control of the process. To reduce this variability, an amplification factor set by the user is applied to the deviation to produce a larger artificial deviation which produces filter weighting factors that weight the previous moving average error 403 higher and the instantaneous error signal 405 lower. This reduces the variability of the new moving average error 403. The preferred embodiment utilizes an artifidal deviation developed by application of both the cusum divisor and the amplification factor, when required. Occasionally, the calculations will fail for one reason or another. This may produce extremely large errors that do not represent normal operation. It is important to prevent these errors from causing abnormal increases in the moving average variance. One extremely large error can result in distortion of the moving average variance that may last for weeks. To prevent distortion of the moving average variance due to model failures, the fixed weight multiplying the current value of variance in the filter is divided by the deviation if the deviation is greater than a threshold value. The preferred embodiment specifies a limit of 5. One minus the reduced weight multiplies the previous moving average variance in the filter calculation. In the case of flow data reconciliation, the variable weighting factors applied to each term in a sum of weighted squared errors represent the accuracy of the measurement relative to its true value. The further from the true value, the smaller the weighting factor. The function of the variable weighting factor is to allow signals from failed or highly inaccurate sensors to be discounted relative to the accurate measurements. This allows the correction factors for those erroneous signals to absorb essentially the entire correction and thus prevents distortion of the other correction factors. The preferred embodiment utilizes a variable weighting factor that is the redprocal of one plus the square of the product of an adjustable constant and the error. Proper choice of the adjustable constant produces a weighting factor that remains relatively high within a band around the average, and then declines rapidly outside that band. Within the preferred embodiment, the adjustable factor is 4. The variable weighting factors are set prior to the computation that minimizes the sum of weighted error values and are not varied during the computation. Moving average variable weighting factors may be calculated from the instantaneous weighting factors using an exponentially weighted moving average filter with the weighting factors determined as in the case of the moving average error, but without applying cusum counts or amplification factors. However, in this case, the application of the weighting factors is reversed so that a large deviation will result in heavily weighting the current value, rather than the moving average value, so as to respond immediately to a sudden failure of a sensor. The filtered difference 403 is introduced into the product property corrector 402 along with the stream composition/product property estimates 117 from the simulation 110 to produce a corrected stream composition product property estimate 401. The product property corrector module 402 corrects the predicted product property signals 117 by application of the product property average error 403 according to the manner by which it was derived, to produce the corrected product property signals 401, which in turn are used as input to the product property controllers 240 and by the flow data reconciliation module 108. In the event there is no filtered difference 403 corresponding to a stream composition product property estimate 117, then the property estimate 117 is passed unchanged through product property corrector 402 as property estimate 401. The calibration module 406 adjusts the model calibration factors 411 to drive the product error 405 and the product property average error 403 toward zero. The calibration module 406 may also generate setpoint signals 509 for adjusting setpoints of alignment controllers 316. The preferred method for adjusting calibration factors is through use of the standard PID control procedure, which can be tuned to avoid changing the calculated product properties to the extent that they alter the operation of the plant significantly, while still obtaining a reasonable adjustment rate. Other types of control procedures, such as lead-lag and dead time compensation, may be employed if their use is appropriate for the situation. The control procedure can be triggered two ways. Usually, the model calibration parameters are updated only when a new analytical value is measured. In some cases it may be desirable to trigger the control procedure each calculation cycle to move the model toward the last measured value much more quickly, but without upsetting the process. B. Alignment Figure 3 shows a block diagram of a model alignment system 318 and property control system 270. The alignment system 318 utilizes comparisons between measured and computed values of the process variables, normally temperature and pressure to produce corrective action that drives the differences toward zero. The corrective action usually involves adjusting a selected flow, or model parameter. Sensors 106 are used to generate raw pressure, temperature, level and flow measurement signals 109 from process 104. Temperature and pressure are normally the variables sensitive to stream composition and/or product property variations. The alignment system 318 provides a means for keeping a model in alignment with the process whenever pressure and/or temperature are composition sensitive. Either temperature is maintained by the model and the alignment function maintains pressure by biasing flows or feed compositions, or vice versa, whenever analyses are infrequent or non existent Analog PTLF signals are converted to digital signals and along with digital PTLF signals are input to the alignment system 318. The signals 109 are first limit checked and bounded at validity check 130. A first principles simulation 110 of the process is carried out on the computer using the reconciled flows and validated temperature, pressure, and level measurements 115, along with analytical determinations of the feed composition 119. The process variable comparison module 312 determines the difference between the reconciled pressure and temperature signals 115 and the computed regulatory pressure and temperature signals 113, usually as an arithmetic difference, relative difference, or relative difference based on the logarithm of the ratio of the raw process variable signals 115 to the predicted process variable signals 113. This comparison produces a differential signal 315. If the differential signal 315 is highly variable, then a moving average value of the difference is maintained by the filter 308 and updated each compute cycle. If the difference is fairly stable, then the difference 315 may pass unchanged through filter 308. A second embodiment of the present invention uses the validated signal 111 as a direct input into process variable comparison 312, completely bypassing the flow data reconciliation 108. The statistical filter 308 updates the value of the average error 309 and average variance for each co; .utation cycle for each variable processed, using the same methods as the product property statistical filter module 404 (i.e., the filter first subtracts the previous filtered signal from the current signal to produce a measure of the variance, which is combined with the average variance to determine the deviation). The average error value 309 is used by the alignment module 318. The statistical filter 308 receives a signal each compute cycle. The filtered difference 309 is then processed by an alignment module 316. The alignment module 316 provides a regulatory control adjustment 313 to the simulation models 110 so that they remain synchronized with the process operations by driving the process variable difference 315 and the filtered difference 309 toward zero. The preferred method for adjusting the model is through the use of a PID control procedure, although other types of control procedures may be used, if appropriate. A simple example can be used to illustrate how the alignment function operates. Consider a binary distillation separation where there is a nearly equal molar mixture of the two components in the base of the column. The two components have significantly different boiling points, so that changes in composition at constant pressure are readily detected as changes in temperature. On the actual column the base pressure is held constant, while a controller adjusts the tails rate to maintain a constant base temperature, and thus, a constant composition. The model that simulates this column sets the base temperature of the model equal to the base temperature of the process. The simulated base composition and pressure are then determined from the simulation. If the simulated composition is in error, then the calculated base pressure will not match the actual base pressure. To maintain the base pressure, and hence the base composition, aligned with the actual process an alignment controller is used to adjust the model tails flow to drive the difference between the base temperature of the model and that of the actual process to zero. Since the base pressure measurement may be quite noisy and have a significant fixed error, it is processed through the statistical filter to eliminate the noise and obtain a reliable estimate of the true base pressure. Thus, with the model's base temperature and pressure both closely ma ching process conditions, and the estimated tails concentrations should closely track the actual concentrations. At the top of the column the concentration of the high boiling component is so low that concentration changes cannot be detected by temperature changes. Overhead stream composition values from the model are provided directly to the overhead product property (stream composition) controller. In this case no alignment functionality is present and without a periodic analysis of the overhead stream, there is no way to calibrate the model or directly correct model composition estimates. If, however, there are periodic measurements made of overhead and bottoms streams, then the calibration function can be used to both correct simulation model estimates and to adjust elements of the model so as to drive the difference between computed and measured concentrations toward zero. In the case of the overhead stream, the estimated concentration of the high boiling component may be directly adjusted with the moving average composition error term before being introduced into the overhead product property controllers. The same error term can be utilized as input to a calibration controller that adjusts the model value of tray efficiency above the feed point on the column as a way to cause measured and computed composition values in the overhead stream to match. Likewise, the measure of tails composition can be compared with the computed value and a moving average error determined that can be used directly to correct the estimated composition before sending it to the bottoms product property (composition) control. The same moving average error could also be used as input to a calibration controller that adjusts the tray effidency below the feed tray to drive the moving average error towards zero. Finally, the bottoms moving average composition error might be used to adjust the setpoint of the alignment controller. If for example, the pressure reading has a significant fixed error, a calibration controller can be used to place a bias in the setpoint of the alignment controller to correct for the error. Thus, a continuous, but noisy pressure measurement with a fixed error can be used to maintain the model in alignment with the process. This approach is particularly useful if the base of the column is sampled infrequently, so that the samples are inadequate to provide the degree of calibration necessary to successfully track the tails composition between analyses. This example illustrates how the models in conjunction with the calibration and alignment functions can be used in various combinations to improve the estimates of stream composition and/or product properties provided the product property controllers. The combination of the alignment 318 and the calibration and product property corrector 408 systems is to improve the composition/product property estimates in several ways. First, the alignment function 318, where process measurements are linked to product properties and/or stream compositions, maintains models in alignment with the process where composition analyses are unavailable or infrequent Second, the calibration and product property corrector system produces artificial analyses to augment the normal analysis set This increases the number of calibration points throughout the model and enhances the ability of the models to adapt to changing conditions automatically. Third, the product property comparison can utilize alternative methods to determine errors. The use of errors determined in different - 31 - forms can be utilized to linearize the relationship between the error and the response of the model to the variable or parameter manipulated by the calibration controllers, and thus provide more accurate and effective control of the models. Fourth, the statistical filter provides a means for effectively utilizing laboratory and analyzer data to detect and eliminate the impact of bad analytical results. By effectively smoothing laboratory results relative to model results, the latest filtered difference can be applied directly to correcting model output with much less risk of upsetting process operations due to a bad analysis. Fifth, by providing the filter the ability to adapt to process shifts and protect itself against gross errors, greatly enhances the accuracy and reliability of product property estimates. Sixth, The calibration controllers can be tuned so that they recalibrate the model slowly enough to minimize the impact of calibration effects on plant operations. The various elements of the invention can be used individually or in combination. Not shown in the figures is an optional historical database. The historical database maintains a history of raw flow sensor signals 109, analytical values 119, and may optionally store a history of corrected flow values 115, predicted product property signals 117, or any of the other signal values of the present invention. The historical database can serve as a data buffer providing data to the invention. It can also provide historical data which may be used by the present invention for various enhanced signal processing tasks. An historical database is used in the preferred embodiment, but functions largely as a data buffer, and thus is not essential to the invention or the preferred embodiment III. Sequence of Modules The process autopilot carries out all its functions in a cycle, executed at a fixed time interval, with idle time passing until the next cycle. During each cycle, the processes of the various modules are carried out in sequence. On the first cycle, when no previous computed values are available, the autopilot initializes by executing the process simulation models 110 using pre-defined estimates for all sensor values. This provides a starting point from which the simulation will converge cycle-by-cycle to eventually align with the process. Then, whether the first or later cycle, the validity check 130 executes. Then the flow data reconciliation 108 is executed using validated raw sensor measurements 111 with previously determined corrected product property signals 401, computed regulatory variables 113, and predicted product property signals 117, to produce corrected flow signals 115. On each cycle product property transform 414 executes if new lab or analyzer data is available or if stream composition and/or product property values are to be synthesized from reconciled flow measurements. Then product property comparison 410, statistical filter for product properties 404 and calibration 406 are carried out This corrects for drift of predicted product property signals 117 from the product property signals 119. Also, at each cycle, process variable comparison 312, statistical filter for process variables 308 and alignment 316 are carried out This also corrects for drift of predicted product property signals 117 from the product property signals 119. The product property corrector 402 is then executed to produce the desired corrected product property signals 401. Then, product property controllers 240 are provided with the predicted product property signals 401. The simulation models 110 are then executed to produce estimates for the next cycle. Product property controllers 240 operate on their own cycle, and may be synchronized with the provision of updated predicted product property signals 401, or may function synchronously using the latest predicted product property signals 401, repeatedly if needed. PTLF controllers 250 operate independently of the autopilot or the product property controllers 240. Whenever a product property controllers 240 provides a new setpoint 245, the PTLF controller 250 continue their independent operation using the new setpoint While the invention has been particularly shown and described with reference to preferred embodiments thereof, it will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the spirit and scope of the invention.

Claims: What is claimed is: 1. A computer-based method of processing flow signals, said method employing data reconciliation, comprising the steps of: (1) measuring, in a process comprising one or more units and two or more streams, flow rates and temperatures of two or more of said streams to produce raw flow signals; (2) predicting, using dynamic simulation, predicted flow signals for one of more of said streams not measured in step (1); (3) constructing a complete mass balance around at least one of said units, by using said raw flow signals and said predicted flow signals; (4) correcting one or more of said raw flow signals and/or predicting flow signals by using said complete mass balance to produce mass corrected flow signals. 2. The computer-based method of processing flow signals of claim 1, further comprising the steps of: predicting, using said dynamic simulation, predicted values of rate of change of inventory in one or more of said units; and constructing said complete mass balance(s) of step (3) using a dynamic balance(s). 3. The computer-based method of processing flow signals of claim 1, further comprising the steps of: predicting, using said dynamic simulation, predicted composition values in one or more of said streams; and constructing a complete enthalpy balance around at least one of said one or more units, by computing stream enthalpies using said mass corrected flow signals, said predicted composition values, and said raw temperature signals; and further correcting one of more of said mass corrected flow signals to produce enthalpy corrected flow signals. 4. The computer based method of claim 1, further comprising the steps: operating said process to produce a product; and using said mass corrected flow signals to directly or indirectly change the state of a final control element of said process in accordance with an operating objective. 5. The computer-based method of processing flow signals of claim 1, further comprising the steps: operating the process to create said flow rates; changing the state of a final control element of said process in accordance with a predicted product property signal; and wherein said predicting step (2) further predicts said predicted product property signal. 6. A computer-based data recondliation system, comprising: (a) two or more sensors for producing raw flow signals indicative of flow rates of two or more streams in a process comprising one or more units and two or more streams; (b) a computer-based dynamic simulation, connected to use mass corrected flow signals and said raw flow signals as inputs, for producing predicted flow signals for one or more said streams; and (c) a computer-based flow data reconciliation module, connected to use said predicted flow signals and said raw flow signals as inputs, for correcting said raw flow signals to minimize a mass error sum of a complete mass balance around at least one of said one or more units, to produce said mass corrected flow signals. 7. The computer-based data reconciliation system claim 6, wherein said computer-based dynamic simulation further produces predicted rate of change of inventory signals in one or more of said units; and wherein said computer-based flow data reconciliation module is further connected to use said predicted rate of change of inventory signals, for correcting said raw flow signals in accordance with a dynamic complete mass balance. 8. The computer-based data reconciliation system of claim 6, wherein: said computer-based dynamic simulation further produces predicted composition signals for one or more of said streams; and further comprising: a computer-based flow data reconciliation module, connected to use said mass corrected flow signals said raw temperature signals and said predicted composition signals as inputs, for further correcting said mass corrected flow signals to minimize an enthalpy error sum of a complete enthalpy balance around at least one of said one or more units. 9. The computer-based data reconciliation system of claim 6, further comprising: a controller, for changing a state of a final control element of said process, connected to directly or indirectly use said corrected flow signals as inputs, and connected to change said state in accordance with an operating objective. 10. The computer-based data reconciliation system of claim 6, wherein said dynamic simulation further produces a predicted product property signal; and further comprising: a product property controller,connected to directly or indirectly use said predicted product property signal as input, for producing a controller output signal in accordance with said input and a control objective; and an actuator, for changing a state of the process, connected to change said state in accordance with said controller output signal. 11. A computer-based method of processing flow signals, said method employing data reconciliation, comprising the steps of: (1) measuring, in a process comprising one or more units and two or more streams, flow rates of two or more of said streams to produce raw flow signals; (2) constructing a complete mass balance around at least one of said units, using said raw flow signals; (3) adjusting each of said raw flow signals by multiplying each signal by a respective constant; and (4) adjusting said respective constants to minimize a sum of values, each value indicative of the deviation of each said respective constant from a value of one. 12. The computer-based method of processing flow signals of claim 11, wherein: said values of step (4) are computed by multiplying each said deviation times a respective weight, and further comprising a step of computing each said value of step (4), and reducing said respective weight of each said value whenever said value exceeds a threshold. 13. The computer-based method of processing flow signals of claim 11, wherein said values of step (4) are computed by multiplying each said deviation times a respective weight, and further comprising the steps of: computing a statistical characteristic of said value; and computing each said value of step (4), and reducing said respective weight of each said value whenever said value is high relative to said statistical characteristic. 14. A computer-based data reconciliation system, comprising: (a) two or more sensors for producing raw flow signals indicative of flow rates of two or more streams in a process comprising one or more units; (b) a computer-based flow data reconciliation module, connected to use said raw flow signak as inputs, for correcting said raw flow signals to produce mass corrected flow signals by multiplying each signal by a respective constant and adjusting said respective constants to minimize a sum of values, each value indicative of the deviation of each said respective constant from a value of one. 15. The computer-based data reconciliation system of claim 14, wherein said computer-based flow data reconciliation module corrects said raw flow signals by multiplying each said deviation times a respective weight, and by reducing said respective weight of each said value whenever said value exceeds a threshold. 16. The computer-based data reconciliation system of claim 14, wherein said computer-based flow data reconciliation module corrects said raw flow signals by computing a statistical characteristic of said value and reducing said respective weight of each said value whenever said value is high relative to said statistical characteristic.

Inventors: CARDNER DAVID V, CARDNER, DAVID, V.
Applicants: DU PONT, E.I. DU PONT DE NEMOURS & CO. (INC.)
Application year: 1993
Main classification: G06D5
IPC codes: G01F1, G01F15, G01F25, G05B13, G05B17, G06F19