US20030065462A1 - Multivariate statistical process analysis systems and methods for the production of melt polycarbonate - Google Patents
Multivariate statistical process analysis systems and methods for the production of melt polycarbonate Download PDFInfo
- Publication number
- US20030065462A1 US20030065462A1 US09/682,273 US68227301A US2003065462A1 US 20030065462 A1 US20030065462 A1 US 20030065462A1 US 68227301 A US68227301 A US 68227301A US 2003065462 A1 US2003065462 A1 US 2003065462A1
- Authority
- US
- United States
- Prior art keywords
- variables
- predetermined
- process variables
- predetermined process
- melt
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B23/00—Testing or monitoring of control systems or parts thereof
- G05B23/02—Electric testing or monitoring
- G05B23/0205—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults
- G05B23/0218—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults
- G05B23/0243—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults model based detection method, e.g. first-principles knowledge model
- G05B23/0245—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults model based detection method, e.g. first-principles knowledge model based on a qualitative model, e.g. rule based; if-then decisions
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2135—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on approximation criteria, e.g. principal component analysis
Definitions
- the present invention is directed to systems and methods for the control of chemical manufacturing processes and, more specifically, to multivariate statistical process analysis systems and methods for the production of melt polycarbonate.
- Manufacturing process variables (X i ) such as flows, pressures, concentrations, temperatures, and others have traditionally been subject to statistical process control (SPC) strategies.
- SPC statistical process control
- These SPC strategies include a variety of statistical methods designed to maintain process quality and productivity.
- the statistical methods focus on a single variable X i at a time, using univariate controls such as Shewhart charts, cumulative sum charts, and exponentially-weighted moving average charts.
- Such charts are used to monitor the performance of a process, such as a chemical manufacturing process, over time to verify that the process consistently operates within the production specifications of a given product.
- MSPA multivariate statistical process analysis
- process chemometrics The application of multivariate statistical analysis methods to industrial process data characterized by a large number of correlated chemical process measurements is the area of process chemometrics.
- the objectives of process chemometrics include the determination of key process variables, the generation of inference models used to forecast and optimize product quality, the detection and diagnosis of faults and potential process abnormalities, and the overall monitoring of chemical processes to ensure production control. Achieving these goals is often difficult with regard to the production of melt polycarbonate, however, as the determination of key process variables may be an inexact and time consuming process, and accurate and reliable inference models may be difficult to generate.
- the present invention is directed to automated multivariate statistical process analysis systems and methods for the production of melt polycarbonate.
- a computerized system for the production of melt polycarbonate includes a plurality of sensors for obtaining a plurality of measurements relating to a plurality of predetermined process variables, a preprocessor for preprocessing each of the plurality of measurements for multivariate statistical analysis, an identifier for identifying which of the plurality of predetermined process variables affect each of a plurality of predetermined product variables, a correlator for correlating the plurality of predetermined process variables and the plurality of predetermined product variables, a model generator for modeling the relationship between the plurality of predetermined process variables and the plurality of predetermined product variables, and an analyzer for analyzing the plurality of predetermined process variables to predict polymer performance and/or to ensure that the value of each of the plurality of predetermined process variables is within a predetermined range.
- a computerized method for the production of melt polycarbonate includes the steps of obtaining a plurality of measurements relating to a plurality of predetermined process variables, preprocessing each of the plurality of measurements for multivariate statistical analysis, identifying which of the plurality of predetermined process variables affect each of a plurality of predetermined product variables, correlating the plurality of predetermined process variables and the plurality of predetermined product variables, modeling the relationship between the plurality of predetermined process variables and the plurality of predetermined product variables, and analyzing the plurality of predetermined process variables to predict polymer performance and/or to ensure that the value of each of the plurality of predetermined process variables is within a predetermined range.
- FIG. 1 is a flow chart of a computerized multivariate statistical process analysis method for the production of melt polycarbonate
- FIG. 2 is a plot of the eigenvalues of process variables, measured using the method of FIG. 1, after autoscaling;
- FIG. 3 is a graph of the percent of captured variance for each process variable for the first two principal components (PCs) of the principal components analysis (PCA) model of the present invention
- FIG. 4 is a loadings plot from the PCA model of FIG. 3;
- FIG. 5 is a loadings plot for the first PC from the PCA model of FIG. 3;
- FIG. 6 is a loadings plot for the fourth PC from the PCA model of FIG. 3;
- FIG. 7 is a plot of the first two scores for the PCA model of FIG. 3;
- FIG. 8 is a plot of the result of the prediction of pellet intrinsic viscosity (IV) using only process variables
- FIG. 9 is a plot of the result of the prediction of Fries concentration using only process variables
- FIG. 10 is a Q control chart for the multivariate statistical process analysis of the production of melt polycarbonate resin
- FIG. 11 is a T 2 control chart for the multivariate statistical process analysis of the production of melt polycarbonate resin.
- FIG. 12 is a functional block diagram of a computerized multivariate statistical process analysis system for the production of melt polycarbonate.
- Polycarbonates are typically prepared from dihydric phenol compounds and carbonic acid derivatives.
- melt polycarbonate may be prepared via the melt polymerization of diphenyl carbonate and Bisphenol A (BPA). The reaction is conducted at high temperatures, allowing the starting monomers and product to remain molten while the reactor pressure is staged in order to more effectively remove phenol, the by-product of the polycondensation reaction.
- BPA Bisphenol A
- melt polycarbonate manufacturing process data may be collected via sensors in order to monitor process performance. Using this collected information, the relative importance of various process variables (X i ) may be determined and an inference model may be developed to predict such outcomes as Fries concentration, pellet intrinsic viscosity (IV), melt polycarbonate grade, and other product parameters. Further, process variables X i causing abnormal performance may be detected and identified. As a result, a manufacturing plant staff may better understand process performance and make sound business decisions.
- a multivariate statistical process analysis (MSPA) method 10 for the consistent production of melt polycarbonate begins with the collection of sensor data over a predetermined period of time 12 .
- data from one melt polycarbonate manufacturing plant included information about process variables X i and product variables Y i from about 343 polymerization runs. For each polymerization run, process variables X i and product variables Y i were recorded once.
- a list of process variables X i and product variables Y i used for multivariate analysis is presented in Table 1. TABLE 1 Process and Product Variables Process Variables X i Product Variables Y i 1. Molar Ratio 1. Fries Concentration 2. Flow Rate 2. Melt Flow Ratio (MFR) 3.
- R 3 and R 4 refer to specific polymerization stages within the physical melt polycarbonate manufacturing system. An exemplary description of such polymerization stages is presented in International Patent Application WO 00/37531.
- Viscosity (IV) of a sample in methylene chlorine End Cap (EC) % Calculated from measured concentration of terminal OH groups Yellowness Index — Yellowness index of pellets (YI) Melt Polycarbonate — Based upon Fries concentration, MFR, Grade and EC
- gathered data may be preprocessed 14 .
- the process data may be arranged as a matrix of samples and product variables Y i .
- the product variables Y i may be used to label data points after applying pattern recognition tools.
- the process variables X i and product variables Y i in the data set may, however, have different physical measurement units. This difference in units between the variables may be eliminated by autoscaling the data 16 .
- Autoscaling is typically the application of both variance scaling and mean centering.
- the data may be autoscaled by subtracting the mean and dividing by the standard deviation for each variable:
- each column of the autoscaled matrix has zero mean and unit standard deviation.
- PCA principal components analysis
- the preprocessed matrix X a is decomposed into a set of scores T (I ⁇ K) and loadings P (J ⁇ K), where K is the number of principal components.
- K is the number of principal components.
- the number of PCs is chosen to be as few as possible to explain important variation in the data set.
- the matrix X a is then expressed as a linear combination of orthogonal vectors along the directions of the PCs:
- the percent variance captured by the PCA model may be analyzed (see Table 3 below) in combination with a plot of eigenvalues as a function of PCs 21 (see FIG. 2).
- four PCs captured about 80.8% of the variance in the data.
- Higher-order PCs had eigenvalues less than unity and did not describe any systematic variance in the system.
- four PCs were determined to be adequate for the PCA model.
- PC 1 22 describes about 43% of the total variance in the data, but no single process variable X i provides an exclusive contribution.
- the five process variables X i that individually contributed more than about 50% of the captured variance in PC 1 22 were discharge pressure in R 3 , temperature in R 4 , stirring speed in R 4 , discharge pressure in R 4 , and melt viscosity in R 4 .
- the two process variables X that individually contributed more than about 50% of the captured variance in PC 2 24 were flow rate and throughput.
- a loadings plot 26 may also be used as a diagnostic tool for the PCA model and process variables X i .
- the loadings plot 26 determines which variables are important for describing variation in the data set. When a given process variable X i contributes a significant variation to the PC, the absolute value of the loading of the variable will be close to unity. When the process variable X i does not contribute a significant variation to the PC, the absolute value of the loading of the variable will be close to zero.
- the different signs of the loadings indicate that the process variables inversely contribute to the variance described by the PC. In the above example, discharge pressure in R 4 contributes the most to PC 1 , while vacuum in R 3 contributes the least. Most of the contribution to PC 2 is provided by flow rate and throughput, while vacuum in R 4 and adjusted molar ratio contribute the least to PC 2 .
- FIG. 5 presents the values of loadings of process variables X i to PC 1 28 and FIG. 6 presents the values of loadings of process variables X i to PC 4 30 .
- vacuum in R 4 has a large loading of about ⁇ 0.5 to PC 4 30 and a relatively small loading of about ⁇ 0.2 to PC 1 28 .
- PC 1 28 describes about seven times more variation in the data set.
- Another step in the multivariate statistical process analysis method 10 of the present invention is establishing a correlation between process variables X i 32 .
- Multivariate methods are capable of detecting changes in the correlation structure of a group of variables that may not be detected by univariate methods.
- the correlation structure in the data set may be visualized using, for example, a correlation analysis.
- R is the correlation coefficient and N is the number of data points.
- Results of the correlation analysis of process variables X i may be presented as pseudocolor maps, and may, optionally, be reordered using, for example, a k-nearest neighbor (KNN) cluster analysis. This hierarchical cluster analysis determines the similarity of process variables X i based on their measured properties.
- KNN k-nearest neighbor
- a further step in the multivariate statistical process analysis method 10 of the present invention is establishing a correlation between process variables X i and product variables Y i 34 .
- An immediate benefit of this step 34 may include the ability to forecast product quality based upon measurements involving only process variables X i .
- material quality may be monitored as a function of several product variables X i .
- Initial analysis of the correlation structure in the combined data set of process variables X i and product variables Y i may also be performed using correlation analysis.
- a pseudocolor correlation map of the process variables X i and product variables Y i may be reordered using KNN cluster analysis.
- a strong inverse correlation (R ⁇ 1) was found between the stirring speed in R 4 and pellet IV.
- FIG. 7 shows the plot of the first two scores for PCA 36 .
- the plot 36 illustrates a general relationship between samples of melt polycarbonate resin manufacture and different process variable X i conditions.
- Several possible clusters of data may be assigned on the scores plot 36 . Intuitively, polycarbonate produced under similar process conditions should have similar properties.
- a more in-depth understanding of the relationships between process variables X i and product variables Y i may be obtained given additional knowledge of the data.
- the data points on the scores plots may be labeled according to information about product variables Y i .
- Histogram plots may be used to determine the distribution of product variables Y i and to assign labels to respective portions of the distribution. For example, the correlation between process variables X i captured by the first two PCs and product variables Y i such as melt flow ratio, pellet IV, Fries concentration, and melt polycarbonate grade may be plotted.
- Analysis of histograms of these product variables Y i may, for example, indicate that melt flow ratio and pellet IV are more rightly arranged in clusters on the histogram plots than, for example, Fries concentration. This variation in distribution may be explained by more pronounced effects of the process variables X i on Fries formation.
- the multivariate statistical process analysis method 10 of the present invention further includes the generation of an inference model for predicting and analyzing the melt polycarbonate production process and its performance 38 .
- pellet IV and Fries concentration of manufactured polycarbonate resin may be predicted using “virtual analyzers,” analyzing only the information from process variables X i .
- a multivariate calibration method such as partial least-squares (PLS) regression may be used.
- PLS partial least-squares
- the quality of the developed PLS models for quantitation of pellet IV and Fries concentration may be evaluated using, for example, the root mean squared error of calibration (RMSEC).
- RMSEC root mean squared error of calibration
- Performance of the PLS models developed for quantitation of pellet IV and Fries concentration may also be validated using a leave-one-out cross-validation algorithm.
- the root mean squared error of cross-validation (RMSECV) may be used to estimate the ability of the models to predict pellet IV and Fries concentration.
- RMSECV is essentially the standard deviation of the predicted values minus laboratory estimated values (i.e. the standard deviation of the test set residuals). A large RMSECV indicates poor correlation with the reference method and/or poor precision.
- the quality of collected process variables X i may be evaluated using statistical tools such as multivariate control charts and multivariate contribution plots, among others. This allows for the detection of faults and the diagnosis of problems in the process variables X i 40 .
- Multivariate control charts use two statistical indicators of the PCA model, such as Q and T 2 values plotted as a function of manufactured sample. The significant principal components of the PCA model are used to develop the T 2 -chart and the remaining PCs contribute to the Q-chart. The Q residual is the squared prediction error and describes how well the PCA model fits each sample. It is a measure of the amount of variation in each sample not captured by K principal components retained in the model:
- t i is the ith row of T k
- ⁇ ⁇ 1 is the diagonal matrix containing the inverse of the eigenvalues associated with the k eigenvectors (PCs) retained in the model.
- the T 2 -chart monitors the multivariate distance of a new sample from a target value in the reduced PCA space.
- the Q and T 2 control charts may be statistical indicators in the multivariate statistical process control of the production of melt polycarbonate resin 42 and 44 . These charts illustrate samples that exceed some predetermined confidence limit, for example the 95% confidence limit, described by the PCA model.
- the contribution plots of Q residuals and Hotelling′′s T statistic may provide an indication of which process variables X i cause problems in a given sample. A large Q residual may occur, for example, due to data collection errors or process disturbances.
- a multivariate statistical process analysis (MSPA) system 42 for the consistent production of melt polycarbonate includes a plurality of sensors 44 for collecting manufacturing process data and a computer 46 for determining the relative importance of various process variables X i and developing an inference model to predict such outcomes as Fries concentration, pellet intrinsic viscosity (IV), and melt polycarbonate grade.
- the computer 46 may also be used to detect and identify process variables X i causing abnormal performance on-line.
- the system 42 of the present invention thus allows a manufacturing plant staff to monitor process performance, better understand it, and make sound business decisions.
- the computer 46 typically includes inputs/outputs, a memory, and a processor for receiving, sending, storing, and processing signals and data to operate, monitor, record, and otherwise functionally control the operation of the system 42 .
- the computer 46 may include software, hardware, firmware, and other similar components for functionally controlling the operation of the system 42 .
- the computer 46 may be a single device, or it may be a plurality of devices working in concert.
- the computer 46 is preferably in communication with all of the other components of the system 42 .
- the input/output devices may include, for example, a keyboard and a mouse for entering data and instructions into the computer 46 .
- a video display allows the user or process operator to view what the computer 46 has accomplished.
- the memory typically includes a random-access memory (RAM) and a read-only memory (ROM).
- the memory may also include other types of memory, such as programmable read-only memory (PROM), erasable programmable read-only memory (EPROM), and electrically erasable programmable read-only memory (EEPROM).
- PROM programmable read-only memory
- EPROM erasable programmable read-only memory
- EEPROM electrically erasable programmable read-only memory
- the memory also preferably includes an operating system that executes on the processor. The operating system performs basic tasks which include recognizing input from input devices, sending output to output devices, keeping track of files and directories, and controlling various peripheral devices.
- the memory may also contain one or more databases.
- the processor accepts data and instructions from the memory and performs various calculations.
- the processor may include an arithmetic logic unit (ALU), which performs arithmetic and logical operations, and a control unit, which extracts instructions from the memory.
- ALU arithmetic logic unit
- the computer 46 may also include a modem or other network connection, a mass storage device, and any other suitable peripheral.
- the above-described computer 46 may take the form of a hand-held digital computer, a personal computer, a workstation, a mainframe computer, and a supercomputer.
- the computer's memory preferably contains a number of programs or algorithms for functionally controlling the operation of the system 42 , including a preprocessor 48 , for preprocessing collected process variable X i and product variable Y i data, and an identifier 50 , for identifying the process variables X i of importance in, for example, the melt polycarbonate manufacturing process.
- the preprocessor 48 may, for example, include an algorithm for scaling each of the measurements related to the process variables X i and product variables Y i
- the identifier 50 may, for example, include an algorithm for performing principal components analysis (PCA).
- PCA principal components analysis
- the computer 46 may also contain a correlator 52 , for performing a correlation analysis between process variables X i and process variables X i and product variables Y i , and a model generator 54 , for generating a “virtual analyzer,” used to understand and predict the performance of the manufacturing process.
- the model generator 54 may, for example, include an algorithm for performing partial least-squares (PLS) regression.
- the computer 46 may contain an analyzer 56 for detecting faults and diagnosing problems with monitored process variables X i that may lead to inferior polycarbonate product. This may be accomplished through the use of multivariate control charts.
- Process and product data may be analyzed using multivariate techniques included in, for example, a chemometrics software package, such as PLS_Toolbox (Version 2.0, Eigenvector Research, Inc., Manson Wash.). This software package may operate, for example, with MATLAB (Version 5.3, The Mathworks, Inc., Natick, Mass.).
Abstract
Description
- The present invention is directed to systems and methods for the control of chemical manufacturing processes and, more specifically, to multivariate statistical process analysis systems and methods for the production of melt polycarbonate.
- Manufacturing process variables (Xi) such as flows, pressures, concentrations, temperatures, and others have traditionally been subject to statistical process control (SPC) strategies. These SPC strategies, introduced in the 1930's, include a variety of statistical methods designed to maintain process quality and productivity. The statistical methods focus on a single variable Xi at a time, using univariate controls such as Shewhart charts, cumulative sum charts, and exponentially-weighted moving average charts. Such charts are used to monitor the performance of a process, such as a chemical manufacturing process, over time to verify that the process consistently operates within the production specifications of a given product.
- As the number of monitored variables Xi affecting the behavior of a manufacturing process increases, however, univariate SPC methods become inadequate. The application of these techniques in such situations may result in misleading information being presented to a process operator, leading him or her to take unnecessary or erroneous control actions.
- An alternative approach is to employ multivariate statistical process analysis (MSPA) methods to extract more relevant information from measured data. MSPA methods provide the staff of a manufacturing plant, for example, with a greater understanding of process performance, allowing them to make sound business decisions. Thus, the application of multivariate methodologies to industrial manufacturing processes has experienced increasing popularity in recent years. For example, MSPA methods have been utilized in emulsion polymerization, low-density continuous polyethylene polymerization, batch polymerization, and pilot-scale penicillin fermentation processes. Similarly, MSPA methods have been utilized to improve the productivity of a titanium dioxide plant, monitor the processing conditions of a nuclear waste storage tank, and control the performance of chromatographic instrumentation.
- The application of multivariate statistical analysis methods to industrial process data characterized by a large number of correlated chemical process measurements is the area of process chemometrics. The objectives of process chemometrics include the determination of key process variables, the generation of inference models used to forecast and optimize product quality, the detection and diagnosis of faults and potential process abnormalities, and the overall monitoring of chemical processes to ensure production control. Achieving these goals is often difficult with regard to the production of melt polycarbonate, however, as the determination of key process variables may be an inexact and time consuming process, and accurate and reliable inference models may be difficult to generate.
- Thus, the present invention is directed to automated multivariate statistical process analysis systems and methods for the production of melt polycarbonate.
- These systems and methods allow process variables causing abnormal performance to be detected and identified. As a result, a manufacturing plant staff may better understand process performance and make sound business decisions.
- In one embodiment, a computerized system for the production of melt polycarbonate includes a plurality of sensors for obtaining a plurality of measurements relating to a plurality of predetermined process variables, a preprocessor for preprocessing each of the plurality of measurements for multivariate statistical analysis, an identifier for identifying which of the plurality of predetermined process variables affect each of a plurality of predetermined product variables, a correlator for correlating the plurality of predetermined process variables and the plurality of predetermined product variables, a model generator for modeling the relationship between the plurality of predetermined process variables and the plurality of predetermined product variables, and an analyzer for analyzing the plurality of predetermined process variables to predict polymer performance and/or to ensure that the value of each of the plurality of predetermined process variables is within a predetermined range.
- In another embodiment, a computerized method for the production of melt polycarbonate includes the steps of obtaining a plurality of measurements relating to a plurality of predetermined process variables, preprocessing each of the plurality of measurements for multivariate statistical analysis, identifying which of the plurality of predetermined process variables affect each of a plurality of predetermined product variables, correlating the plurality of predetermined process variables and the plurality of predetermined product variables, modeling the relationship between the plurality of predetermined process variables and the plurality of predetermined product variables, and analyzing the plurality of predetermined process variables to predict polymer performance and/or to ensure that the value of each of the plurality of predetermined process variables is within a predetermined range.
- FIG. 1 is a flow chart of a computerized multivariate statistical process analysis method for the production of melt polycarbonate;
- FIG. 2 is a plot of the eigenvalues of process variables, measured using the method of FIG. 1, after autoscaling;
- FIG. 3 is a graph of the percent of captured variance for each process variable for the first two principal components (PCs) of the principal components analysis (PCA) model of the present invention;
- FIG. 4 is a loadings plot from the PCA model of FIG. 3;
- FIG. 5 is a loadings plot for the first PC from the PCA model of FIG. 3;
- FIG. 6 is a loadings plot for the fourth PC from the PCA model of FIG. 3;
- FIG. 7 is a plot of the first two scores for the PCA model of FIG. 3;
- FIG. 8 is a plot of the result of the prediction of pellet intrinsic viscosity (IV) using only process variables;
- FIG. 9 is a plot of the result of the prediction of Fries concentration using only process variables;
- FIG. 10 is a Q control chart for the multivariate statistical process analysis of the production of melt polycarbonate resin;
- FIG. 11 is a T2 control chart for the multivariate statistical process analysis of the production of melt polycarbonate resin; and
- FIG. 12 is a functional block diagram of a computerized multivariate statistical process analysis system for the production of melt polycarbonate.
- Polycarbonates are typically prepared from dihydric phenol compounds and carbonic acid derivatives. For example, one important polycarbonate, melt polycarbonate, may be prepared via the melt polymerization of diphenyl carbonate and Bisphenol A (BPA). The reaction is conducted at high temperatures, allowing the starting monomers and product to remain molten while the reactor pressure is staged in order to more effectively remove phenol, the by-product of the polycondensation reaction.
- During the melt polycarbonate manufacturing process, data may be collected via sensors in order to monitor process performance. Using this collected information, the relative importance of various process variables (Xi) may be determined and an inference model may be developed to predict such outcomes as Fries concentration, pellet intrinsic viscosity (IV), melt polycarbonate grade, and other product parameters. Further, process variables Xi causing abnormal performance may be detected and identified. As a result, a manufacturing plant staff may better understand process performance and make sound business decisions.
- Referring to FIG. 1, a multivariate statistical process analysis (MSPA)
method 10 for the consistent production of melt polycarbonate begins with the collection of sensor data over a predetermined period oftime 12. For example, data from one melt polycarbonate manufacturing plant included information about process variables Xi and product variables Yi from about 343 polymerization runs. For each polymerization run, process variables Xi and product variables Yi were recorded once. A list of process variables Xi and product variables Yi used for multivariate analysis is presented in Table 1.TABLE 1 Process and Product Variables Process Variables Xi Product Variables Y i1. Molar Ratio 1. Fries Concentration 2. Flow Rate 2. Melt Flow Ratio (MFR) 3. Adjusted Molar Ratio 3. Pellet Intrinsic Viscosity (IV) 4. Vacuum in R3 4. End Cap (EC) 5. Temperature in R3 5. Yellowness Index (YI) 6. Torque in R3 6. Melt Polycarbonate Grade 7. Discharge Pressure in R3 8. Melt Viscosity in R3 9. Vacuum in R4 10. Temperature in R4 11. Stirring Speed in R4 12. Torque in R4 13. Discharge Pressure in R4 14. Melt Viscosity in R4 15. Throughput 16. Pellet Speed - In the above table, “R3” and “R4” refer to specific polymerization stages within the physical melt polycarbonate manufacturing system. An exemplary description of such polymerization stages is presented in International Patent Application WO 00/37531.
- The variables used for multivariate analysis are further described in Table 2.
TABLE 2 Description of Variables Variable Units Description Molar Ratio — Ratio of moles of DPC to moles of BPA Flow Rate kg/hour Flow rate of the monomer mix into R1 oligomerization reactor Adjusted Molar — Ratio of moles of DPC (and additional Ratio moles of DPC) to moles of BPA Vacuum in R3 torr Measured in headspace of R3 stage using a pressure gauge Temperature in R3 degrees C. Temperature of reaction components in R3 Torque in R3 N*m Torque on a stirrer in R3 reactor Discharge Pressure torr Pressure after the gear pump in the 3rd in R3 reactor stage Melt Viscosity in poise Melt viscosity of reaction components R3 in R3 Vacuum in R4 torr Measured in headspace of R4 stage using a pressure gauge Temperature in R4 degrees C. Temperature of reaction components in R4 Stirring Speed in R4 RPM Speed of a shaft in R4 reactor Torque in R4 — Torque on a shaft in R4 reactor Discharge Pressure torr Pressure in the 4th reactor in R4 Melt Viscosity in poise Melt viscosity of reaction components R4 in R4 Throughput kg/h Amount of polycarbonate material per unit of time coming from the last reactor stage Pellet Speed m/s Speed of a polycarbonate strand entering the pelletizer Fries Concentration ppm Concentration of Fries product measured by LC Melt Flow Ratio g/10 min Measure of weight in grams extruded (MFR) through a capillary for a 10 min test Pellet Intrinsic dL/g Measured at 20 degrees C. on a solution Viscosity (IV) of a sample in methylene chlorine End Cap (EC) % Calculated from measured concentration of terminal OH groups Yellowness Index — Yellowness index of pellets (YI) Melt Polycarbonate — Based upon Fries concentration, MFR, Grade and EC - Prior to multivariate analysis, gathered data may be preprocessed14. For example, the process data may be arranged as a matrix of samples and product variables Yi. The product variables Yi may be used to label data points after applying pattern recognition tools. The process variables Xi and product variables Yi in the data set may, however, have different physical measurement units. This difference in units between the variables may be eliminated by autoscaling the
data 16. Autoscaling is typically the application of both variance scaling and mean centering. The data may be autoscaled by subtracting the mean and dividing by the standard deviation for each variable: - x ij a=(x y −{overscore (x)} j )/σ
j ′ (1) - where Xa :: is the autoscaled matrix X of size I×J and
- and σj are the mean and standard deviation of the ith column of the original matrix X. Thus, each column of the autoscaled matrix has zero mean and unit standard deviation.
- Following the
data preprocessing step 14, as part of the determination of key process variables step 18, information about the relations between the process variables Xi and their effects on the product variables Yi may be obtained using principal components analysis (PCA)tools 20. PCA is a multivariate data analysis tool that projects a data set onto a subspace of lower dimensionality. In this reduced space, data is represented with reduced colinearity. PCA achieves this objective by describing the variance of the preprocessed data matrix Xa (Eq. 1) in terms of the weighted sums of the original variables with no significant loss of information. These weighted sums of the original variables are called principal components (PCs). The preprocessed matrix Xa is decomposed into a set of scores T (I×K) and loadings P (J×K), where K is the number of principal components. Preferably, the number of PCs is chosen to be as few as possible to explain important variation in the data set. The matrix Xa is then expressed as a linear combination of orthogonal vectors along the directions of the PCs: - x a =t 1 p 1 r +t 2 p 2 4 +. . . +t A p K T +E, (2)
- where ti and pi are the score and loading vectors, respectively, E is a residual matrix that represents random error, and T is the transpose of the matrix.
- To determine the number of principal components to retain in the PCA model, the percent variance captured by the PCA model may be analyzed (see Table 3 below) in combination with a plot of eigenvalues as a function of PCs21 (see FIG. 2). In one example, four PCs captured about 80.8% of the variance in the data. Higher-order PCs had eigenvalues less than unity and did not describe any systematic variance in the system. Thus, four PCs were determined to be adequate for the PCA model.
TABLE 3 Percent Variance Captured by PCA Model Principal Eigenvalue % Variance % Variance Component of CoV (X) this PC Cumulative 1 6.85e+000 42.83 42.83 2 3.58e+000 22.37 65.21 3 1.50e+000 9.37 74.58 4 9.96e−001 6.22 80.80 5 8.57e−001 5.36 86.16 6 5.95e−001 3.72 89.88 7 4.71e−001 2.94 92.82 8 3.49e−001 2.18 95.00 9 2.93e−001 1.83 96.84 10 2.46e−001 1.54 98.37 11 8.81e−002 0.55 98.92 12 5.96e−002 0.37 99.29 13 5.12e−002 0.32 99.61 14 3.42e−002 0.21 99.83 15 1.46e−002 0.09 99.92 16 1.29e−002 0.08 100.00 - Information regarding the amount of variance for each process variable Xi captured by individual PCs in the PCA model of this example is presented in FIG. 3.
PC1 22 describes about 43% of the total variance in the data, but no single process variable Xi provides an exclusive contribution. The five process variables Xi that individually contributed more than about 50% of the captured variance inPC1 22 were discharge pressure in R3, temperature in R4, stirring speed in R4, discharge pressure in R4, and melt viscosity in R4. The two process variables X that individually contributed more than about 50% of the captured variance inPC2 24 were flow rate and throughput. - Referring to FIG. 4, a
loadings plot 26 may also be used as a diagnostic tool for the PCA model and process variables Xi. Theloadings plot 26 determines which variables are important for describing variation in the data set. When a given process variable Xi contributes a significant variation to the PC, the absolute value of the loading of the variable will be close to unity. When the process variable Xi does not contribute a significant variation to the PC, the absolute value of the loading of the variable will be close to zero. The different signs of the loadings indicate that the process variables inversely contribute to the variance described by the PC. In the above example, discharge pressure in R4 contributes the most to PC1, while vacuum in R3 contributes the least. Most of the contribution to PC2 is provided by flow rate and throughput, while vacuum in R4 and adjusted molar ratio contribute the least to PC2. - It is also important to note the amount of variation described by a PC when interpreting loadings. A variable with a large loading value contributes significantly to a particular PC. However, the variable may not be truly important if the PC does not describe a large amount of the variation in the data set. FIG. 5 presents the values of loadings of process variables Xi to
PC1 28 and FIG. 6 presents the values of loadings of process variables Xi toPC4 30. For example, vacuum in R4 has a large loading of about −0.5 toPC4 30 and a relatively small loading of about −0.2 toPC1 28. However,PC1 28 describes about seven times more variation in the data set. - Another step in the multivariate statistical
process analysis method 10 of the present invention is establishing a correlation betweenprocess variables X i 32. Multivariate methods are capable of detecting changes in the correlation structure of a group of variables that may not be detected by univariate methods. The correlation structure in the data set may be visualized using, for example, a correlation analysis. The correlation between a pair of variables x and y is defined as: - where R is the correlation coefficient and N is the number of data points. The correlation coefficient R is between −1 and 1 and is independent of the scale of x and y values. For an exact linear relation between x and y, R=−1 if increasing x values correspond to increasing y values and R=−1 if increasing x values correspond to decreasing y values. R=0 if the variables are independent.
- Results of the correlation analysis of process variables Xi may be presented as pseudocolor maps, and may, optionally, be reordered using, for example, a k-nearest neighbor (KNN) cluster analysis. This hierarchical cluster analysis determines the similarity of process variables Xi based on their measured properties.
- In the example discussed above, two pairs of process variables Xi, temperature in R3 and temperature in R4 and flow rate and throughput, had correlation coefficients close to 1. Stirring speed in R4 and both temperature in R4 and discharge pressure in R4 had correlation coefficients close to −1. A further step in the multivariate statistical
process analysis method 10 of the present invention is establishing a correlation between process variables Xi andproduct variables Y i 34. An immediate benefit of thisstep 34 may include the ability to forecast product quality based upon measurements involving only process variables Xi. As discussed above, during the manufacturing of melt polycarbonate, material quality may be monitored as a function of several product variables Xi. - Initial analysis of the correlation structure in the combined data set of process variables Xi and product variables Yi may also be performed using correlation analysis. For example, a pseudocolor correlation map of the process variables Xi and product variables Yi may be reordered using KNN cluster analysis. In the example discussed above, none of the individual product variables Yi correlated with process variables Xi with a correlation coefficient of about 1. A strong inverse correlation (R˜−1) was found between the stirring speed in R4 and pellet IV.
- A detailed analysis of the correlation between process variables Xi and product variables Yi may be performed using PCA of the process variables Xi and analyzing a PC scores plots. FIG. 7 shows the plot of the first two scores for
PCA 36. Theplot 36 illustrates a general relationship between samples of melt polycarbonate resin manufacture and different process variable Xi conditions. Several possible clusters of data may be assigned on thescores plot 36. Intuitively, polycarbonate produced under similar process conditions should have similar properties. - A more in-depth understanding of the relationships between process variables Xi and product variables Yi may be obtained given additional knowledge of the data. In particular, the data points on the scores plots may be labeled according to information about product variables Yi . Histogram plots may be used to determine the distribution of product variables Yi and to assign labels to respective portions of the distribution. For example, the correlation between process variables Xi captured by the first two PCs and product variables Yi such as melt flow ratio, pellet IV, Fries concentration, and melt polycarbonate grade may be plotted. Analysis of histograms of these product variables Yi may, for example, indicate that melt flow ratio and pellet IV are more rightly arranged in clusters on the histogram plots than, for example, Fries concentration. This variation in distribution may be explained by more pronounced effects of the process variables Xi on Fries formation.
- Referring again to FIG. 1, the multivariate statistical
process analysis method 10 of the present invention further includes the generation of an inference model for predicting and analyzing the melt polycarbonate production process and itsperformance 38. For example, pellet IV and Fries concentration of manufactured polycarbonate resin may be predicted using “virtual analyzers,” analyzing only the information from process variables Xi. For this application, a multivariate calibration method such as partial least-squares (PLS) regression may be used. The quality of the developed PLS models for quantitation of pellet IV and Fries concentration may be evaluated using, for example, the root mean squared error of calibration (RMSEC). Performance of the PLS models developed for quantitation of pellet IV and Fries concentration may also be validated using a leave-one-out cross-validation algorithm. The root mean squared error of cross-validation (RMSECV) may be used to estimate the ability of the models to predict pellet IV and Fries concentration. RMSECV is essentially the standard deviation of the predicted values minus laboratory estimated values (i.e. the standard deviation of the test set residuals). A large RMSECV indicates poor correlation with the reference method and/or poor precision. - In the example discussed above, the results of the prediction of pellet IV and Fries concentration using only process variables Xi are presented in FIGS. 8 and 9, respectively. For quantitation of pellet IV, a six-factor PLS regression model was found to be adequate. This model accounted for about 85% of the variance in process variables Xi and about 98% of the pellet IV variance. For quantitation of Fries concentration, a four-factor PLS regression model was found to be optimal. It accounted for about 75% of the variance in process variables Xi and about 90% of Fries concentration variance. A summary of the performance of the PLS models for the prediction of pellet IV and Fries concentration from the measured process variables Xi is presented in Table 4.
TABLE 4 Summary of Calibration Model Performances Root Mean Root Mean Squared Error Squared Error Product Variable of Calibration of Cross-Validation R2 Pellet IV 0.00703 dL/g 0.00758 dL/g 0.99 Fries Concentration 0.0146% Fries 0.0152% Fries 0.90 - To ensure normal manufacturing plant operation, the quality of collected process variables Xi may be evaluated using statistical tools such as multivariate control charts and multivariate contribution plots, among others. This allows for the detection of faults and the diagnosis of problems in the
process variables X i 40. Multivariate control charts use two statistical indicators of the PCA model, such as Q and T2 values plotted as a function of manufactured sample. The significant principal components of the PCA model are used to develop the T2-chart and the remaining PCs contribute to the Q-chart. The Q residual is the squared prediction error and describes how well the PCA model fits each sample. It is a measure of the amount of variation in each sample not captured by K principal components retained in the model: - Q i =e i e i T =x i(I−P k P k T)X i T, (4)
- where ei is the ith row of E, xi is the ith sample in X, Pk is the matrix of the k loading vectors retained in the PCA model (where each vector is a column of Pk) and I is the identity matrix of appropriate size (n×n). The Q residual chart monitors the deviation from the PCA model for each sample. The sum of normalized squared scores, known as Hotelling″s T2 statistic, is a measure of the variation within the PCA model and allows for the identification of statistically anomalous samples. T2 is defined as:
- T i 2 =t iλ−1 t i T= x i Pλ − P T x i T, (5)
- where ti is the ith row of Tk, the matrix of k scores vectors from the PCA model, and λ−1 is the diagonal matrix containing the inverse of the eigenvalues associated with the k eigenvectors (PCs) retained in the model. The T2-chart monitors the multivariate distance of a new sample from a target value in the reduced PCA space.
- Referring to FIGS. 10 and 11, the Q and T2 control charts may be statistical indicators in the multivariate statistical process control of the production of
melt polycarbonate resin - Referring to FIG. 12, a multivariate statistical process analysis (MSPA)
system 42 for the consistent production of melt polycarbonate includes a plurality ofsensors 44 for collecting manufacturing process data and acomputer 46 for determining the relative importance of various process variables Xi and developing an inference model to predict such outcomes as Fries concentration, pellet intrinsic viscosity (IV), and melt polycarbonate grade. Thecomputer 46 may also be used to detect and identify process variables Xi causing abnormal performance on-line. Thesystem 42 of the present invention thus allows a manufacturing plant staff to monitor process performance, better understand it, and make sound business decisions. - Structurally, the
computer 46 typically includes inputs/outputs, a memory, and a processor for receiving, sending, storing, and processing signals and data to operate, monitor, record, and otherwise functionally control the operation of thesystem 42. Thecomputer 46 may include software, hardware, firmware, and other similar components for functionally controlling the operation of thesystem 42. Thecomputer 46 may be a single device, or it may be a plurality of devices working in concert. Thecomputer 46 is preferably in communication with all of the other components of thesystem 42. The input/output devices may include, for example, a keyboard and a mouse for entering data and instructions into thecomputer 46. A video display allows the user or process operator to view what thecomputer 46 has accomplished. Other output devices may include, for example, a printer, a plotter, a synthesizer, and speakers. The memory typically includes a random-access memory (RAM) and a read-only memory (ROM). The memory may also include other types of memory, such as programmable read-only memory (PROM), erasable programmable read-only memory (EPROM), and electrically erasable programmable read-only memory (EEPROM). The memory also preferably includes an operating system that executes on the processor. The operating system performs basic tasks which include recognizing input from input devices, sending output to output devices, keeping track of files and directories, and controlling various peripheral devices. The memory may also contain one or more databases. The processor accepts data and instructions from the memory and performs various calculations. The processor may include an arithmetic logic unit (ALU), which performs arithmetic and logical operations, and a control unit, which extracts instructions from the memory. Optionally, thecomputer 46 may also include a modem or other network connection, a mass storage device, and any other suitable peripheral. The above-describedcomputer 46 may take the form of a hand-held digital computer, a personal computer, a workstation, a mainframe computer, and a supercomputer. - The computer's memory preferably contains a number of programs or algorithms for functionally controlling the operation of the
system 42, including apreprocessor 48, for preprocessing collected process variable Xi and product variable Yi data, and anidentifier 50, for identifying the process variables Xi of importance in, for example, the melt polycarbonate manufacturing process. Thepreprocessor 48 may, for example, include an algorithm for scaling each of the measurements related to the process variables Xi and product variables Yi Theidentifier 50 may, for example, include an algorithm for performing principal components analysis (PCA). Thecomputer 46 may also contain acorrelator 52, for performing a correlation analysis between process variables Xi and process variables Xi and product variables Yi, and amodel generator 54, for generating a “virtual analyzer,” used to understand and predict the performance of the manufacturing process. Themodel generator 54 may, for example, include an algorithm for performing partial least-squares (PLS) regression. Further, thecomputer 46 may contain ananalyzer 56 for detecting faults and diagnosing problems with monitored process variables Xi that may lead to inferior polycarbonate product. This may be accomplished through the use of multivariate control charts. Process and product data may be analyzed using multivariate techniques included in, for example, a chemometrics software package, such as PLS_Toolbox (Version 2.0, Eigenvector Research, Inc., Manson Wash.). This software package may operate, for example, with MATLAB (Version 5.3, The Mathworks, Inc., Natick, Mass.). - The present invention has been described with reference to examples and preferred embodiments. Other examples and embodiments may achieve the same results. Variations in and modifications to the present invention will be apparent to those skilled in the art and the following claims are intended to cover all such equivalents.
Claims (33)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US09/682,273 US6549864B1 (en) | 2001-08-13 | 2001-08-13 | Multivariate statistical process analysis systems and methods for the production of melt polycarbonate |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US09/682,273 US6549864B1 (en) | 2001-08-13 | 2001-08-13 | Multivariate statistical process analysis systems and methods for the production of melt polycarbonate |
Publications (2)
Publication Number | Publication Date |
---|---|
US20030065462A1 true US20030065462A1 (en) | 2003-04-03 |
US6549864B1 US6549864B1 (en) | 2003-04-15 |
Family
ID=24738970
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US09/682,273 Expired - Fee Related US6549864B1 (en) | 2001-08-13 | 2001-08-13 | Multivariate statistical process analysis systems and methods for the production of melt polycarbonate |
Country Status (1)
Country | Link |
---|---|
US (1) | US6549864B1 (en) |
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040058081A1 (en) * | 2001-10-23 | 2004-03-25 | Potyrailo Radislav Alexandrovich | Systems and methods for the deposition and curing of coating compositions |
US20050128703A1 (en) * | 2003-12-15 | 2005-06-16 | Chang Shye-Bin S. | Controlling cooling air intake for air cooled equipment |
WO2006034179A2 (en) * | 2004-09-17 | 2006-03-30 | Mks Instruments, Inc. | Method and apparatus for multivariate control of semiconductor manufacturing processes |
US20060095232A1 (en) * | 2004-11-02 | 2006-05-04 | Purdy Matthew A | Fault detection through feedback |
EP1654600A1 (en) * | 2003-08-13 | 2006-05-10 | Cargill, Incorporated | Computer-aided modeling and manufacture of products |
WO2006051096A2 (en) * | 2004-11-12 | 2006-05-18 | Schott Ag | Device and method for optimizing technical production processes by compacting, consolidating and analyzing process/quality data |
EP1665076A2 (en) * | 2003-09-12 | 2006-06-07 | Tokyo Electron Limited | Method and system of diagnosing a processing system using adaptive multivariate analysis |
US20070021859A1 (en) * | 2005-07-07 | 2007-01-25 | Mks Instruments, Inc. | Self-correcting multivariate analysis for use in monitoring dynamic parameters in process environments |
EP1751628A1 (en) * | 2004-05-20 | 2007-02-14 | McMaster University | Method for controlling the appearance of products and process performance by image analysis |
WO2007064860A1 (en) * | 2005-12-01 | 2007-06-07 | Siemens Corporate Research, Inc. | Robust sensor correlation analysis for machine condition monitoring |
US20080082195A1 (en) * | 2006-09-29 | 2008-04-03 | Fisher-Rosemount Systems, Inc. | Univariate method for monitoring and analysis of multivariate data |
WO2008042757A2 (en) * | 2006-09-29 | 2008-04-10 | Fisher-Rosemount Systems, Inc. | Univariate method for monitoring and analysis of multivariate data |
WO2008042758A2 (en) * | 2006-09-29 | 2008-04-10 | Fisher-Rosemount Systems, Inc. | Multivariate monitoring and diagnostics of process variable data |
US20080109090A1 (en) * | 2006-11-03 | 2008-05-08 | Air Products And Chemicals, Inc. | System And Method For Process Monitoring |
WO2008106071A1 (en) | 2007-02-27 | 2008-09-04 | Exxonmobil Research And Engineering Company | Method and system of using inferential measurements for abnormal event detection in continuous industrial processes |
US20100057237A1 (en) * | 2008-09-02 | 2010-03-04 | Mks Instruments, Inc. | Automated model building and batch model building for a manufacturing process, process monitoring, and fault detection |
US20100191361A1 (en) * | 2009-01-23 | 2010-07-29 | Mks Instruments, Inc. | Controlling a Manufacturing Process with a Multivariate Model |
US20100209307A1 (en) * | 2007-08-07 | 2010-08-19 | Drabish Jeffrey B | method for improving the prediction of polymer properties and a system having improved polymer property prediction capabilities |
US8271103B2 (en) | 2007-05-02 | 2012-09-18 | Mks Instruments, Inc. | Automated model building and model updating |
US8855804B2 (en) | 2010-11-16 | 2014-10-07 | Mks Instruments, Inc. | Controlling a discrete-type manufacturing process with a multivariate model |
US20150365251A1 (en) * | 2014-06-16 | 2015-12-17 | Azbil Corporation | Gateway device and packet capturing method |
US9429939B2 (en) | 2012-04-06 | 2016-08-30 | Mks Instruments, Inc. | Multivariate monitoring of a batch manufacturing process |
US20160371137A1 (en) * | 2015-06-22 | 2016-12-22 | Hitachi Power Solutions Co., Ltd. | Preprocessor of Abnormality Sign Diagnosing Device and Processing Method of the Same |
US9541471B2 (en) | 2012-04-06 | 2017-01-10 | Mks Instruments, Inc. | Multivariate prediction of a batch manufacturing process |
CN107541716A (en) * | 2016-06-27 | 2018-01-05 | 台湾积体电路制造股份有限公司 | The parameter regulation means of deposition manufacture process |
WO2018063805A1 (en) * | 2016-09-27 | 2018-04-05 | General Electric Company | Method for detection and isolation of faulty sensors |
KR20180134696A (en) * | 2017-06-09 | 2018-12-19 | 주식회사 엘지화학 | Method for predicting content of residual monomer and preparation method of polycarbonate resin |
Families Citing this family (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1197130C (en) * | 2000-07-04 | 2005-04-13 | 东京毅力科创株式会社 | Operation monitoring method for treatment apparatus |
TWI266249B (en) * | 2002-05-16 | 2006-11-11 | Mosel Vitelic Inc | Statistical process control method and system thereof |
EP1543392A1 (en) * | 2002-06-28 | 2005-06-22 | Umetrics AB | Method and device for monitoring and fault detection in industrial processes |
JP2004047885A (en) * | 2002-07-15 | 2004-02-12 | Matsushita Electric Ind Co Ltd | Monitoring system and monitoring method of semiconductor manufacturing apparatus |
US20050144274A1 (en) * | 2003-12-12 | 2005-06-30 | General Electric Company | Apparatus for monitoring the performance of a distributed system |
US7079984B2 (en) | 2004-03-03 | 2006-07-18 | Fisher-Rosemount Systems, Inc. | Abnormal situation prevention in a process plant |
WO2006001416A1 (en) * | 2004-06-29 | 2006-01-05 | Nikon Corporation | Management method, management system, and program |
CN101156119B (en) | 2005-04-04 | 2011-04-13 | 费希尔-罗斯蒙德系统公司 | Diagnostics in industrial process control system |
KR100640663B1 (en) | 2005-08-12 | 2006-11-01 | 삼성전자주식회사 | Method and apparatus for modeling multivariate parameters having constants and same patterns and method for fabricating semiconductor using the same |
JP2007250748A (en) * | 2006-03-15 | 2007-09-27 | Omron Corp | Apparatus, method and program of analyzing process abnormality |
US7954072B2 (en) * | 2006-05-15 | 2011-05-31 | Taiwan Semiconductor Manufacturing Company, Ltd. | Model import for electronic design automation |
US7580768B2 (en) * | 2006-09-22 | 2009-08-25 | Texas Instruments Deutschland Gmbh | Method of adjusting process variables in a processing flow |
US8145337B2 (en) * | 2007-05-04 | 2012-03-27 | Taiwan Semiconductor Manufacturing Company, Ltd. | Methodology to enable wafer result prediction of semiconductor wafer batch processing equipment |
US7974728B2 (en) * | 2007-05-04 | 2011-07-05 | Taiwan Semiconductor Manufacturing Company, Ltd. | System for extraction of key process parameters from fault detection classification to enable wafer prediction |
US7783999B2 (en) * | 2008-01-18 | 2010-08-24 | Taiwan Semiconductor Manufacturing Company, Ltd. | Electrical parameter extraction for integrated circuit design |
CN101504543B (en) * | 2008-02-05 | 2011-04-20 | 台湾积体电路制造股份有限公司 | Extraction of key process parameter |
US8037575B2 (en) * | 2008-02-28 | 2011-10-18 | Taiwan Semiconductor Manufacturing Company, Ltd. | Method for shape and timing equivalent dimension extraction |
US8001494B2 (en) | 2008-10-13 | 2011-08-16 | Taiwan Semiconductor Manufacturing Company, Ltd. | Table-based DFM for accurate post-layout analysis |
US8806386B2 (en) * | 2009-11-25 | 2014-08-12 | Taiwan Semiconductor Manufacturing Company, Ltd. | Customized patterning modulation and optimization |
US8745554B2 (en) * | 2009-12-28 | 2014-06-03 | Taiwan Semiconductor Manufacturing Company, Ltd. | Practical approach to layout migration |
US20110288808A1 (en) | 2010-05-20 | 2011-11-24 | International Business Machines Corporation | Optimal test flow scheduling within automated test equipment for minimized mean time to detect failure |
US8855959B2 (en) | 2010-08-30 | 2014-10-07 | International Business Machines Corporation | Integrated cross-tester analysis and real-time adaptive test |
US10896388B2 (en) | 2011-07-25 | 2021-01-19 | Prevedere, Inc. | Systems and methods for business analytics management and modeling |
US10176533B2 (en) | 2011-07-25 | 2019-01-08 | Prevedere Inc. | Interactive chart utilizing shifting control to render shifting of time domains of data series |
US10740772B2 (en) | 2011-07-25 | 2020-08-11 | Prevedere, Inc. | Systems and methods for forecasting based upon time series data |
US8762301B1 (en) | 2011-10-12 | 2014-06-24 | Metso Automation Usa Inc. | Automated determination of root cause |
JP5648157B2 (en) * | 2011-12-28 | 2015-01-07 | 株式会社日立ハイテクノロジーズ | Semiconductor manufacturing equipment |
US9378197B1 (en) * | 2012-08-16 | 2016-06-28 | Gmg Holdings, Llc | Statistical analysis method for automatically selecting a statistical analysis algorithm based on data values and input parameters |
Family Cites Families (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2628905B2 (en) | 1988-12-06 | 1997-07-09 | 日本ジーイープラスチックス 株式会社 | Method for producing polycarbonate |
US5212765A (en) | 1990-08-03 | 1993-05-18 | E. I. Du Pont De Nemours & Co., Inc. | On-line training neural network system for process control |
US5257206A (en) | 1991-04-08 | 1993-10-26 | Praxair Technology, Inc. | Statistical process control for air separation process |
US5339257A (en) | 1991-05-15 | 1994-08-16 | Automated Technology Associates Inc. | Real-time statistical process monitoring system |
US5351202A (en) | 1992-02-27 | 1994-09-27 | International Business Machines Corporation | Evaluation and ranking of manufacturing line non-numeric information |
JP3187957B2 (en) | 1992-07-16 | 2001-07-16 | 日本ジーイープラスチックス株式会社 | Manufacturing method of aromatic polycarbonate |
US5477444A (en) | 1992-09-14 | 1995-12-19 | Bhat; Naveen V. | Control system using an adaptive neural network for target and path optimization for a multivariable, nonlinear process |
DE9212406U1 (en) * | 1992-09-15 | 1992-12-17 | Roehm Gmbh, 6100 Darmstadt, De | |
US5311759A (en) | 1992-12-14 | 1994-05-17 | Ford Motor Company | Method and system for real-time statistical process monitoring |
US5423199A (en) | 1992-12-14 | 1995-06-13 | Ford Motor Company | Method and apparatus for monitoring stamping press process |
CA2118885C (en) | 1993-04-29 | 2005-05-24 | Conrad K. Teran | Process control system |
GB9318875D0 (en) | 1993-09-11 | 1993-10-27 | Jankovic Ljubomir | Self-learning method and apparatus for dynamic system characterisation and simulation |
US5408405A (en) | 1993-09-20 | 1995-04-18 | Texas Instruments Incorporated | Multi-variable statistical process controller for discrete manufacturing |
US5442562A (en) * | 1993-12-10 | 1995-08-15 | Eastman Kodak Company | Method of controlling a manufacturing process using multivariate analysis |
JPH07175535A (en) | 1993-12-16 | 1995-07-14 | Nec Corp | Power circuit for fet amplifier |
GB9326354D0 (en) | 1993-12-23 | 1994-02-23 | British Aerospace | Methods and apparatus for the testing,monitoring and improvement of manufacturing process effectiveness |
US5440478A (en) | 1994-02-22 | 1995-08-08 | Mercer Forge Company | Process control method for improving manufacturing operations |
JP3161921B2 (en) | 1994-10-27 | 2001-04-25 | ファナック株式会社 | Product quality influence factor analysis method and equipment, molding condition adjustment method, product quality judgment item selection method |
DE4446966A1 (en) | 1994-12-28 | 1996-07-04 | Itt Ind Gmbh Deutsche | Production control information system |
US5864773A (en) * | 1995-11-03 | 1999-01-26 | Texas Instruments Incorporated | Virtual sensor based monitoring and fault detection/classification system and method for semiconductor processing equipment |
US5658423A (en) | 1995-11-27 | 1997-08-19 | International Business Machines Corporation | Monitoring and controlling plasma processes via optical emission using principal component analysis |
US5864483A (en) | 1996-08-01 | 1999-01-26 | Electronic Data Systems Corporation | Monitoring of service delivery or product manufacturing |
JP3576335B2 (en) | 1996-10-31 | 2004-10-13 | 松下電器産業株式会社 | Method and apparatus for extracting abnormalities in process processing steps |
US5850339A (en) | 1996-10-31 | 1998-12-15 | Giles; Philip M. | Analysis of data in cause and effect relationships |
US5841676A (en) | 1996-11-04 | 1998-11-24 | General Electric Company | System and method for estimating a change point time in a manufacturing process |
US5905989A (en) | 1996-11-27 | 1999-05-18 | Bently Nevada Corporation | Knowledge manager relying on a hierarchical default expert system: apparatus and method |
US5987398A (en) | 1998-04-30 | 1999-11-16 | Sony Corporation | Method and apparatus for statistical process control of machines and processes having non-constant mean of a response variable |
ATE395373T1 (en) | 1998-12-18 | 2008-05-15 | Gen Electric | OPTICAL DATA CARRIER CONTAINING POLYCARBONATE WITH REDUCED FRIESE CONTENT |
-
2001
- 2001-08-13 US US09/682,273 patent/US6549864B1/en not_active Expired - Fee Related
Cited By (66)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040058081A1 (en) * | 2001-10-23 | 2004-03-25 | Potyrailo Radislav Alexandrovich | Systems and methods for the deposition and curing of coating compositions |
EP1654600A1 (en) * | 2003-08-13 | 2006-05-10 | Cargill, Incorporated | Computer-aided modeling and manufacture of products |
EP1654600A4 (en) * | 2003-08-13 | 2008-06-18 | Cargill Inc | Computer-aided modeling and manufacture of products |
US7328126B2 (en) | 2003-09-12 | 2008-02-05 | Tokyo Electron Limited | Method and system of diagnosing a processing system using adaptive multivariate analysis |
EP1665076A2 (en) * | 2003-09-12 | 2006-06-07 | Tokyo Electron Limited | Method and system of diagnosing a processing system using adaptive multivariate analysis |
EP1665076A4 (en) * | 2003-09-12 | 2007-06-20 | Tokyo Electron Ltd | Method and system of diagnosing a processing system using adaptive multivariate analysis |
US20050128703A1 (en) * | 2003-12-15 | 2005-06-16 | Chang Shye-Bin S. | Controlling cooling air intake for air cooled equipment |
US7003377B2 (en) * | 2003-12-15 | 2006-02-21 | Lucent Technologies Inc. | Controlling cooling air intake for air cooled equipment |
EP1751628A1 (en) * | 2004-05-20 | 2007-02-14 | McMaster University | Method for controlling the appearance of products and process performance by image analysis |
US8131020B2 (en) * | 2004-05-20 | 2012-03-06 | Mcmaster University | Method for controlling the appearance of products and process performance by image analysis |
US20080013821A1 (en) * | 2004-05-20 | 2008-01-17 | Macgregor John F | Method For Controlling The Appearance Of Products And Process Performance By Image Analysis |
EP1751628A4 (en) * | 2004-05-20 | 2007-12-05 | Univ Mcmaster | Method for controlling the appearance of products and process performance by image analysis |
WO2006034179A3 (en) * | 2004-09-17 | 2006-05-18 | Mks Instr Inc | Method and apparatus for multivariate control of semiconductor manufacturing processes |
US7151976B2 (en) | 2004-09-17 | 2006-12-19 | Mks Instruments, Inc. | Multivariate control of semiconductor processes |
US20060111804A1 (en) * | 2004-09-17 | 2006-05-25 | Mks, Instruments, Inc. | Multivariate control of semiconductor processes |
WO2006034179A2 (en) * | 2004-09-17 | 2006-03-30 | Mks Instruments, Inc. | Method and apparatus for multivariate control of semiconductor manufacturing processes |
US20060095232A1 (en) * | 2004-11-02 | 2006-05-04 | Purdy Matthew A | Fault detection through feedback |
US8676538B2 (en) * | 2004-11-02 | 2014-03-18 | Advanced Micro Devices, Inc. | Adjusting weighting of a parameter relating to fault detection based on a detected fault |
WO2006051096A3 (en) * | 2004-11-12 | 2007-02-08 | Schott Ag | Device and method for optimizing technical production processes by compacting, consolidating and analyzing process/quality data |
WO2006051096A2 (en) * | 2004-11-12 | 2006-05-18 | Schott Ag | Device and method for optimizing technical production processes by compacting, consolidating and analyzing process/quality data |
US20070021859A1 (en) * | 2005-07-07 | 2007-01-25 | Mks Instruments, Inc. | Self-correcting multivariate analysis for use in monitoring dynamic parameters in process environments |
US7809450B2 (en) | 2005-07-07 | 2010-10-05 | Mks Instruments, Inc. | Self-correcting multivariate analysis for use in monitoring dynamic parameters in process environments |
WO2007064860A1 (en) * | 2005-12-01 | 2007-06-07 | Siemens Corporate Research, Inc. | Robust sensor correlation analysis for machine condition monitoring |
US20070162241A1 (en) * | 2005-12-01 | 2007-07-12 | Siemens Corporate Research, Inc. | Robust Sensor Correlation Analysis For Machine Condition Monitoring |
US7769561B2 (en) | 2005-12-01 | 2010-08-03 | Siemens Corporation | Robust sensor correlation analysis for machine condition monitoring |
US20080082195A1 (en) * | 2006-09-29 | 2008-04-03 | Fisher-Rosemount Systems, Inc. | Univariate method for monitoring and analysis of multivariate data |
US8489360B2 (en) | 2006-09-29 | 2013-07-16 | Fisher-Rosemount Systems, Inc. | Multivariate monitoring and diagnostics of process variable data |
WO2008042758A2 (en) * | 2006-09-29 | 2008-04-10 | Fisher-Rosemount Systems, Inc. | Multivariate monitoring and diagnostics of process variable data |
US7917240B2 (en) | 2006-09-29 | 2011-03-29 | Fisher-Rosemount Systems, Inc. | Univariate method for monitoring and analysis of multivariate data |
US20080082194A1 (en) * | 2006-09-29 | 2008-04-03 | Fisher-Rosemount Systems, Inc. | On-line multivariate analysis in a distributed process control system |
WO2008042758A3 (en) * | 2006-09-29 | 2008-07-24 | Fisher Rosemount Systems Inc | Multivariate monitoring and diagnostics of process variable data |
US20080082302A1 (en) * | 2006-09-29 | 2008-04-03 | Fisher-Rosemount Systems, Inc. | Multivariate detection of abnormal conditions in a process plant |
US20080082181A1 (en) * | 2006-09-29 | 2008-04-03 | Fisher-Rosemount Systems, Inc. | Statistical signatures used with multivariate analysis for steady-state detection in a process |
WO2008042757A3 (en) * | 2006-09-29 | 2009-01-29 | Fisher Rosemount Systems Inc | Univariate method for monitoring and analysis of multivariate data |
WO2008042757A2 (en) * | 2006-09-29 | 2008-04-10 | Fisher-Rosemount Systems, Inc. | Univariate method for monitoring and analysis of multivariate data |
US20080082299A1 (en) * | 2006-09-29 | 2008-04-03 | Fisher-Rosemount Systems, Inc. | On-line monitoring and diagnostics of a process using multivariate statistical analysis |
US8014880B2 (en) | 2006-09-29 | 2011-09-06 | Fisher-Rosemount Systems, Inc. | On-line multivariate analysis in a distributed process control system |
US20080082297A1 (en) * | 2006-09-29 | 2008-04-03 | Fisher-Rosemount Systems, Inc. | Multivariate monitoring and diagnostics of process variable data |
US7966149B2 (en) | 2006-09-29 | 2011-06-21 | Fisher-Rosemount Systems, Inc. | Multivariate detection of transient regions in a process control system |
US7937164B2 (en) | 2006-09-29 | 2011-05-03 | Fisher-Rosemount Systems, Inc. | Multivariate detection of abnormal conditions in a process plant |
US7853339B2 (en) | 2006-09-29 | 2010-12-14 | Fisher-Rosemount Systems, Inc. | Statistical signatures used with multivariate analysis for steady-state detection in a process |
US7853431B2 (en) | 2006-09-29 | 2010-12-14 | Fisher-Rosemount Systems, Inc. | On-line monitoring and diagnostics of a process using multivariate statistical analysis |
US20080109090A1 (en) * | 2006-11-03 | 2008-05-08 | Air Products And Chemicals, Inc. | System And Method For Process Monitoring |
WO2008106071A1 (en) | 2007-02-27 | 2008-09-04 | Exxonmobil Research And Engineering Company | Method and system of using inferential measurements for abnormal event detection in continuous industrial processes |
US8285513B2 (en) * | 2007-02-27 | 2012-10-09 | Exxonmobil Research And Engineering Company | Method and system of using inferential measurements for abnormal event detection in continuous industrial processes |
EP2132605A4 (en) * | 2007-02-27 | 2011-05-18 | Exxonmobil Res & Eng Co | Method and system of using inferential measurements for abnormal event detection in continuous industrial processes |
US20080281557A1 (en) * | 2007-02-27 | 2008-11-13 | Emigholz Kenneth F | Method and system of using inferential measurements for abnormal event detection in continuous industrial processes |
EP2132605A1 (en) * | 2007-02-27 | 2009-12-16 | ExxonMobil Research and Engineering Company | Method and system of using inferential measurements for abnormal event detection in continuous industrial processes |
US8271103B2 (en) | 2007-05-02 | 2012-09-18 | Mks Instruments, Inc. | Automated model building and model updating |
US8862438B2 (en) | 2007-08-07 | 2014-10-14 | Dow Global Technologies Llc | Method for improving the prediction of polymer properties |
US20100209307A1 (en) * | 2007-08-07 | 2010-08-19 | Drabish Jeffrey B | method for improving the prediction of polymer properties and a system having improved polymer property prediction capabilities |
US8494798B2 (en) | 2008-09-02 | 2013-07-23 | Mks Instruments, Inc. | Automated model building and batch model building for a manufacturing process, process monitoring, and fault detection |
US20100057237A1 (en) * | 2008-09-02 | 2010-03-04 | Mks Instruments, Inc. | Automated model building and batch model building for a manufacturing process, process monitoring, and fault detection |
US20100191361A1 (en) * | 2009-01-23 | 2010-07-29 | Mks Instruments, Inc. | Controlling a Manufacturing Process with a Multivariate Model |
US9069345B2 (en) | 2009-01-23 | 2015-06-30 | Mks Instruments, Inc. | Controlling a manufacturing process with a multivariate model |
US8855804B2 (en) | 2010-11-16 | 2014-10-07 | Mks Instruments, Inc. | Controlling a discrete-type manufacturing process with a multivariate model |
US9541471B2 (en) | 2012-04-06 | 2017-01-10 | Mks Instruments, Inc. | Multivariate prediction of a batch manufacturing process |
US9429939B2 (en) | 2012-04-06 | 2016-08-30 | Mks Instruments, Inc. | Multivariate monitoring of a batch manufacturing process |
US20150365251A1 (en) * | 2014-06-16 | 2015-12-17 | Azbil Corporation | Gateway device and packet capturing method |
US20160371137A1 (en) * | 2015-06-22 | 2016-12-22 | Hitachi Power Solutions Co., Ltd. | Preprocessor of Abnormality Sign Diagnosing Device and Processing Method of the Same |
US10496466B2 (en) * | 2015-06-22 | 2019-12-03 | Hitachi Power Solutions Co., Ltd. | Preprocessor of abnormality sign diagnosing device and processing method of the same |
CN107541716A (en) * | 2016-06-27 | 2018-01-05 | 台湾积体电路制造股份有限公司 | The parameter regulation means of deposition manufacture process |
WO2018063805A1 (en) * | 2016-09-27 | 2018-04-05 | General Electric Company | Method for detection and isolation of faulty sensors |
US11940318B2 (en) | 2016-09-27 | 2024-03-26 | Baker Hughes Energy Technology UK Limited | Method for detection and isolation of faulty sensors |
KR20180134696A (en) * | 2017-06-09 | 2018-12-19 | 주식회사 엘지화학 | Method for predicting content of residual monomer and preparation method of polycarbonate resin |
KR102047981B1 (en) | 2017-06-09 | 2019-11-22 | 주식회사 엘지화학 | Method for predicting content of residual monomer and preparation method of polycarbonate resin |
Also Published As
Publication number | Publication date |
---|---|
US6549864B1 (en) | 2003-04-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6549864B1 (en) | Multivariate statistical process analysis systems and methods for the production of melt polycarbonate | |
CN108062565B (en) | Double-principal element-dynamic core principal element analysis fault diagnosis method based on chemical engineering TE process | |
Huang et al. | Quality relevant and independent two block monitoring based on mutual information and KPCA | |
US20070088550A1 (en) | Method for predictive maintenance of a machine | |
Said et al. | Machine learning technique for data-driven fault detection of nonlinear processes | |
CN101169623B (en) | Non-linear procedure fault identification method based on kernel principal component analysis contribution plot | |
US6804563B1 (en) | Multidimensional method and system for statistical process control | |
Jiang et al. | Weighted kernel principal component analysis based on probability density estimation and moving window and its application in nonlinear chemical process monitoring | |
Chen et al. | Probabilistic contribution analysis for statistical process monitoring: A missing variable approach | |
EP2287620B1 (en) | Automatic analyzer comprising means for assisting judgment of abnormality of reaction process data | |
Shang et al. | Recursive dynamic transformed component statistical analysis for fault detection in dynamic processes | |
EP3422222A1 (en) | Method and state machine system for detecting an operation status for a sensor | |
WO2022020642A1 (en) | Predicting equipment fail mode from process trace | |
Yu et al. | An improved just-in-time learning scheme for online fault detection of nonlinear systems | |
CN109683594B (en) | Method for accurately identifying and positioning abnormal variable | |
US20220236144A2 (en) | Apparatus for equipment monitoring | |
Oliveira-Santos et al. | Submersible motor pump fault diagnosis system: A comparative study of classification methods | |
US20100161132A1 (en) | System and method for monitoring an industrial production process | |
EP3926427A1 (en) | Work efficiency evaluating method, work efficiency evaluating device, and program | |
Fei et al. | Online process monitoring for complex systems with dynamic weighted principal component analysis | |
Sena et al. | Multivariate statistical analysis and chemometrics | |
CN116661410A (en) | Large-scale industrial process fault detection and diagnosis method based on weighted directed graph | |
US20230384780A1 (en) | Construction method of abnormality diagnosis model, abnormality diagnosis method, construction device of abnormality diagnosis model, and abnormality diagnosis device | |
Rusinov et al. | Real time diagnostics of technological processes and field equipment | |
CN113052232A (en) | Motor fault analysis method and system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: GENERAL ELECTRIC COMPANY, NEW YORK Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:POTYRAILO, RADISLAV A.;REEL/FRAME:011858/0514 Effective date: 20010504 |
|
REMI | Maintenance fee reminder mailed | ||
FPAY | Fee payment |
Year of fee payment: 4 |
|
SULP | Surcharge for late payment | ||
AS | Assignment |
Owner name: SABIC INNOVATIVE PLASTICS IP B.V., NETHERLANDS Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:GENERAL ELECTRIC COMPANY;REEL/FRAME:020820/0578 Effective date: 20070831 |
|
AS | Assignment |
Owner name: CITIBANK, N.A., AS COLLATERAL AGENT, NEW YORK Free format text: SECURITY AGREEMENT;ASSIGNOR:SABIC INNOVATIVE PLASTICS IP B.V.;REEL/FRAME:021423/0001 Effective date: 20080307 Owner name: CITIBANK, N.A., AS COLLATERAL AGENT,NEW YORK Free format text: SECURITY AGREEMENT;ASSIGNOR:SABIC INNOVATIVE PLASTICS IP B.V.;REEL/FRAME:021423/0001 Effective date: 20080307 |
|
REMI | Maintenance fee reminder mailed | ||
LAPS | Lapse for failure to pay maintenance fees | ||
STCH | Information on status: patent discontinuation |
Free format text: PATENT EXPIRED DUE TO NONPAYMENT OF MAINTENANCE FEES UNDER 37 CFR 1.362 |
|
FP | Lapsed due to failure to pay maintenance fee |
Effective date: 20110415 |