Background of the Invention
1. Field of the Invention
This invention relates to a method for preparing lubricating oils by predicting performance, perceptual, chemical or physical properties of streams entering or exiting units in the lubes manufacturing process using a combination of gas chromatography and mass spectrometry.
2. Description of the Related Art
Refineries and chemical plants typically control processing of various component streams and certain additives through the use of both on-line analyzers and off-line laboratory analyses to provide quality information. These quality parameters (chemical composition, physical or perceptual or performance properties) are then fed back into a process control computer which contains control software or algorithms which control refinery or chemical plant hardware (distillation towers or proportional flow control valves). The control programs are typically executed more than three times per hour and their output is used to control either processes or proportional flow valves to vary the quality of the finished product which can either go to tankage or directly to pipelines and terminals, trucks, or to ship-loading facilities. Multiple on-line analyzers are typically required for this process control.
Gas chromatography has been used to predict physical and performance properties of hydrocarbon mixtures boiling in the gasoline range. Crawford and Hellmuth, Fuel, 1990, 69, 443-447, describe the use of gas chromatography and principal components regression analysis to predict the octane values for gasolines blended from different refinery streams. Japanese laid-open patent application JP 03-100463 relates to a method of estimating the cetane number for fuel oils by separating an oil sample into its components using gas chromatography, measuring the signal strength of ion intensities at characteristic masses in the mass spectrum, and correlating these ion intensities to cetane number using multiple regression analysis.
Combined gas chromatography/mass spectrometry (GS/MS) analysis has been done on crude oils. U.S. Pat. No. 5,119,315 discloses a method for aligning sample data such as a mass chromatogram with reference data from a known substance. Williams et al, 12th European Assoc. Organic Geochem., Organic Geochem. Int. Mtg. (Germany 09/16-20/85); Organic Geochemistry 1986, Vol. 10 (1-3) 451-461, discusses the biodegradation of crude oils as measured by GC/MS analysis.
It would be desirable to have a single analyzer means for rapidly measuring chemical or physical properties of intermediate streams in the manufacturing process and/or lube oil products and using these properties to control units within the manufacturing process and/or product quality.
Summary of the Invention
This invention relates to a process for controlling the manufacture of lubricating oils from a feed having a boiling point of about 350.degree. C. or greater by the steps of introducing the feed into a distillation unit, separating the feed into light, medium and heavy distillates, introducing the distillates into an extraction unit, extracting the distillates to produce raffinates, introducing the raffinates into a dewaxing unit and dewaxing the raffinates to produce lubricating oils, said process comprising:
(a) obtaining samples of at least one of lubricating oil, feed, distillate and raffinate;
(b) selecting at least one physical, chemical, perceptual or performance property of at least one of lubricating oil, feed, distillate and raffinate;
(c) selecting reference samples, said reference samples containing characteristic compound types present in the at least one lubricating oil, feed, distillate and raffinate and which have known values for the property or properties selected in step (b);
(d) producing a training set by the steps of:
(1) injecting each reference sample into a gas chromatograph which is interfaced to a mass spectrometer thereby causing at least a partial separation of the hydrocarbon mixture into constituent chemical components;
(2) introducing the constituent chemical components of each reference sample into the mass spectrometer, under dynamic flow conditions;
(3) obtaining for each reference sample a series of time resolved mass chromatograms;
(4) calibrating the mass chromatograms to correct retention times;
(5) selecting a series of corrected retention time windows;
(6) selecting within each retention time window a series of molecular and/or fragment ions, said ions being representative of characteristic compounds or compound classes expected within the retention time window;
(7) recording the total amount of each characteristic compound or compound group selected in step d(6);
(8) forming the data from steps d(6) and d(7) into a X-block matrix;
(9) forming the property data selected in (a) for reference samples selected in (b) into a Y-block matrix;
(10) analyzing the data from steps d(8) and d(9) by multivariate correlation techniques including Partial Least Squares, Principal Component Regression, or Ridge Regression to produce a series of coefficients;
(e) subjecting at least one unknown sample of at least one of lubricating oil, feed, distillate and raffinate to steps d(1) to d(3) in the same manner as the reference sample to produce a series of time resolved mass chromatograms;
(f) repeating steps d(4) and d(8) for each mass chromatogram from step (e);
(g) multiplying the matrix from step (f) by the coefficients from step d(10) to produce a predicted value for the property or properties of the at least one lubricating oil, and/or at least one feed, distillate and raffinate sample or samples; and
(h) using the predicted values of the property or properties of the lubricating oil and/or at least one of the feed, distillate and raffinate sample or samples to control operation of at least one of the distillation unit, extraction unit and dewaxing unit.
Another embodiment of the invention relates to a method for controlling or monitoring chemical or refinery processes which utilize feedstocks and/or produce products having boiling points greater than about 350.degree. C. which comprises
(a) obtaining at least one sample of a refinery or chemical feedstock or product;
(b) selecting at least one physical, chemical, perceptual or performance property of at least one refinery or chemical process feedstock or product;
(c) selecting reference samples, said reference samples containing characteristic compound types present in the at least one refinery or chemical process feedstock or product and which have known values for the property or properties selected in step (b);
(d) producing a training set by the steps of:
(1) injecting each reference sample into a gas chromatograph which is interfaced to a mass spectrometer thereby causing at least a partial separation of the hydrocarbon mixture into constituent chemical components;
(2) introducing the constituent chemical components of each reference sample into the mass spectrometer, under dynamic flow conditions;
(3) obtaining for each reference sample a series of time resolved mass chromatograms;
(4) calibrating the mass chromatograms to correct retention times;
(5) selecting a series of corrected retention time windows;
(6) selecting within each retention time window a series of molecular and/or fragment ions, said ions being representative of characteristic compounds or compound classes expected within the retention time window;
(7) recording the total amount of each characteristic compound or compound group selected in step d(6);
(8) forming the data from steps d(6) and d(7) into a X-block matrix;
(9) forming the property data selected in (a) for reference samples selected in (b) into a Y-block matrix;
(10) analyzing the data from steps d(8) and d(9) by multivariate correlation techniques including Partial Least Squares, Principal Component Regression, or Ridge Regression to produce a series of coefficients;
(e) subjecting at least one unknown refinery or chemical process feedstock or process sample to steps d(1) to d(3) in the same manner as the reference samples to produce a series of time resolved mass chromatograms;
(f) repeating steps d(4) and d(8) for each mass chromatogram from step (e); (g) multiplying the matrix from step (f) by the coefficients from step d(10) to produce a predicted value of the property or properties for the refinery or chemical process feedstocks or products sample or samples; and
(h) using the predicted values of the property or properties of the refinery or chemical feedstocks, intermediate or products sample or samples to control the refinery or chemical process.
The gas chromatography/mass spectrometry (GC/MS) method described above can be used to rapidly predict a wide range of perceptual, performance, chemical and/or physical properties of complex mixtures such as intermediate streams in the lubes manufacturing process and lubricating oil products. Such properties include distillation characteristics, pour point, density, aniline point, feed stock quality, cloud point, haze, viscosity and the like as well as specific chemical constituents. The multivariate correlation methods can also treat the collinear data generated by the GC/MS analyses.
Brief Description of the Drawing
FIG. 1 is a schematic flow diagram of a manufacturing process for a lube oil.
Detailed Description of the Invention
The typical lubricating oil manufacturing process involves a distillation unit, an extraction unit, a dewaxing unit and optionally, a hydrofining unit. The lubricating oil products are commonly used as base stocks for producing finished oils or may be specialty products such as transformer oils and agricultural oils.
FIG. 1 is a schematic flow diagram showing an on-line manufacturing process. Feed stocks 10 are complex hydrocarbon mixtures such as crudes or heavy distillate streams. The feedstock is fed to distillation unit 20 through line 12. Feedstock may be sampled through sampling port 14 which is connected to sampling line 16 interrupted by valve 18.
Distillation unit 20 is typically a pipestill which distills the feed stock into a plurality of distillation cuts. For purposes of simplicity, these distillation cuts are shown exiting unit 20 as line 22 which represents a plurality of lines where n is the number of distinct distillation cuts. Each of the lines 22n may be sampled through sampling ports 24n which are connected to sampling lines 26n interrupted by valves 28n.
The distillate cuts are fed to extraction unit 30 where they are extracted using solvents such as phenol, N-methylpyrrolidone and the like. The extraction process results in raffinates and bottoms. Again for purposes of simplicity, these raffinates are shown exiting the extraction unit 30 as a single line 32n which represents a plurality of lines where n is the number of different raffinate extracts. Each of the lines 32n may be sampled through sampling ports 34n which are connected to sampling lines 36n interrupted by valves 38n. The raffinate extracts are usually stripped of solvent in a stripping unit (not shown) prior to entering the dewaxing unit.
In the dewaxing unit 40, wax is removed by dewaxing aids such as methyl ethyl ketone, propane and the like. As in the previous units, for purposes of simplicity, the dewaxed raffinates are shown exiting the dewaxing unit 40 as a single line 42n where n is the number of dewaxed raffinates. Each of lines 42n may be sampled through sampling ports 44n which are connected to sampling lines 46n interrupted by valves 48n.
The dewaxed raffinates exiting unit 40 through lines 42n may be used as lubricating oils (following removal of dewaxing aid). Typical examples of such lube oils are known as light neutrals, heavy neutrals, bright stocks and the like. These lube oils may be sampled for product quality through sampling ports 44n connected to sampling lines 46n interrupted by valves 48n. These lube oils may optionally be subjected to a hydrofinishing step in a hydrofiner unit 50. Hydrofiners use catalysts and conditions well known in the art. The hydrofinished products exit hydrofiner 50 through lines 52m where m represents the number of lube oils subjected to hydrofinishing. The hydrofinished oils may be sampled through sampling ports 54m connected to sampling lines 56m interrupted by valves 58m.
Each of lines 16, 26n, 36n, 46n and 56m are connected to sampling manifold 60. Sampling manifold 60 also includes an optional provision for at-line operation wherein a remote location 70 is connected to the sampling manifold 60 through line 62 interrupted by valve 66 and sampling port 64. Sampling manifold 60 is a central location where all samples to be analyzed are gathered prior to GC/MS analysis. Manifold 60 can be purged by sample to be analyzed to remove contaminants and purged material recycled through line 72 interrupted by valve 74.
Sample to be analyzed is fed through line 76 to sample injection valve 80. Injection valve 80 can likewise be purged by sample to be analyzed through line 82 interrupted by valve 84. Once sample manifold 60 and sample injection valve 80 are purged, sample to be analyzed is injected into GC/MS analyzer. Samples from off-line operation can be directly injected through valve 80. Off-line samples are those collected at remote locations and brought to a central facility for analysis. The raw data generated by GC/MS analyzer 90 is fed to computer 100. Spent sample and purging materials are sent to recycle through line 92 interrupted by valve 94. Computer 100 mathematically treats raw GC/MS data. Output from computer 100 is then fed as an electronic signal through line 102 to control computer 110. Control computer 110 is connected through line 112 to units 20, 30, 40 and 50 and controls operating conditions within said units 20, 30, 40 and 50. The operating conditions in turn controls the quality of lube oil product exiting units 40 and/or 50.
The operation of the GC/MS analyzer 90 and computer 100 is described in further detail as follows. The feedstock, distillate cuts and lubricating oil products are all complex hydrocarbon mixtures involving thousands of individual chemical species. The use of GC/MS to analyze such complex mixtures for chemical and physical properties involves the generation of large amounts of collinear data. Multiple regression analysis may be employed for treating normal linear data. However, this type of analysis cannot be used for collinear data.
The process according to the invention involves a means for predicting chemical, performance, perceptual and physical properties of feedstock, distillates, raffinates and lubricating oils by quantitative identification of components using a combination of retention times from a GC analysis coupled with target fragment and/or molecular ions produced by the MS. The MS information is compared against a set of known properties from reference samples which form a training set. By mathematically comparing the experimental data with that of the training set, one may predict the desired properties of other feedstocks, distillates, raffinates and lubricating oils.
GC/MS utilizes a gas chromatograph interfaced with a mass spectrometer. While a chromatographic method such as supercritical fluid chromatography, liquid chromatography or size exclusion chromatography may be used to separate the mixture into components or mixtures of components, gas chromatography, especially capillary gas chromatography is the preferred means for interfacing with a mass spectrometer. Both GC and MS use computer software for instrument control, data acquisition and data reduction. The computer platform should be capable of acquiring at least 2000 mass spectra in about 7 minutes.
The sample mixture to be analyzed is first injected into a GC where the mixture components are separated as a function of retention time, often on the basis of boiling points. Only partial chromatographic resolution of mixture components is necessary. The GC oven temperature is usually programmed for samples with a wide boiling range. Separated components may also be identified by a detector such as a flame ionization detector, atomic emission detector, thermal conductivity detector or electron capture detector.
The separated or partially separated components are then transferred to the mass spectrometer under dynamic flow conditions. Since a GC operates under atmospheric pressure and a MS under vacuum conditions (about 10.sup.-3 kPa), the instrument interface requires a coupling device such as a molecular separator (e.g., jet, membrane, etc.), open split coupler or capillary direct interface to efficiently transfer sample while minimizing carrier gas effects.
Depending on the nature of the sample, the mixture may be introduced directly into a MS using a direct insertion probe without a prior GC separation step. Other thermal separation techniques not involving a GC may be used to introduce the sample into the mass spectrometer.
In the MS, sample molecules are bombarded with high energy electrons thereby creating molecular ions which fragment in a pattern characteristic of the molecular species involved. A continuous series of mass spectra are obtained over a scan range of at least 10 to at least 800 daltons. The mass spectral data may also be acquired under selected ion monitoring (SIM) mode. In the selected ion mode, care must be taken to select ions representative of the components of interest and to operate under repeatable conditions. A variety of MS instruments may be used including low resolution, high resolution, MS/MS (hybrid, triple quadrupole, etc.), ion cyclotron resonance and time of flight. Any ionization technique may be used, such as electron ionization, chemical ionization, multiphoton ionization, field desorption, field ionization, etc., provided that the technique provides either molecular or fragment ions which are suitable for use in the analysis procedure.
The results of sample analysis are a series of e mass spectra. The mass spectra are divided into n time intervals where n is an integer from 1 to l. At least one diagnostic ion is chosen from each of m time intervals where m is an integer from 1 to n. The term "diagnostic ion" refers to an ion which is representative of a compound, a chemical class or a performance, perceptual or physical property correlated thereto. Regardless of whether mass spectra are obtained in the scan or selected ion monitoring mode, it is important that the mass spectra be obtained under repeatable conditions.
If the mass spectral data are acquired in the scan mode, the mass range covered during acquisition should be sufficient to provide acquisition of all of the ions which could be used as diagnostic ions during mathematical treatment of each mass spectral scan. If the mass spectral data are acquired in the selected ion monitoring mode, then care must be taken that the ions selected for monitoring are suitable for use in measuring the components of interest.
The sample mass spectral data are then compared to mass spectral data from a series of reference samples with known performance, perceptual, physical and/or chemical properties. In order to compare reference mass spectral data with sample mass spectral data, it may be desirable to time align the sample data to help ensure the integrity of the comparison. There are commercially available computer programs for such data alignment, for example, Hewlett-Packard GC/MS software G 1034C version C.01.05.
The reference mass spectral data, and associated properties data, are arranged in matrix form for mathematical treatment as described below. In the case of chemical composition information, one matrix is formed of reference sample ion intensities at given masses and the other matrix contains known ion intensities for molecular fragment ions of known components. The training set for chemical composition data is thus made up of mass spectral data for different components characteristic of compounds or molecular series expected to be found in the sample mixtures. Similar training sets can be formed for other chemical or perceptual or performance or physical properties of interest. These training sets form one block or matrix of data (Y-block or properties matrix). The actual sample mass spectral data (which may have been time aligned) form the other block (X-block) or matrix of data. These two matrices are subjected to mathematical treatments known as Partial Least Squares (PLS), or Principal Component Regression (PCR), or Ridge Regression (RR) to obtain a mathematically describable relationship between the properties data and mass spectral data, known as a model. Coefficients provided by this model are mathematically combined with the suitably treated mass spectral data from samples with unknown desired properties to:
a) predict desired properties,
b) assess the suitability of the model for such predictions, and
c) diagnose the stability and general correctness of the process that yielded the mass spectral data
PLS/PCR/RR are described in the literature, e.g., Wold S., A. Ruhe, H. Wold, and W. J. Dunn, "The Collinearity Problem in Linear Regression. The Partial Least Squares (PLS) Approach to Generalized Inverses", SIAM J. Sci. Stat. Comput., 1984 5(3), 735-743, or Geladi P., and B. R. Kowalki, "Partial Least Squares Regression: A Tutorial", Anal. Chim. Acta, 1986, 185, 1-17, or Hokuldsson A., "PLS Regression Methods", J. Chemometrics, 1988, 2, 211-228, or in many other articles in journals such as the Journal of Chemometrics or Intelligent Laboratory Systems; Frank, I. and J. Friedman, "A Statistical View of Some Chemometrics Regression Tools", Technometrics, 1993, Vol. 35, No. 2; Jackson, J. E., "A User's Guide To Principal Components", Wiley-Interscience, New York, 1991; Montogomery, D. C. and E. A. Peck, "Introduction To Linear Regression Analysis", Wiley-Interscience, New York, 1990; and Martens, H., and T. Naes, "Multivariable Calibration", Wiley-Interscience, New York,1989.
When dealing with a complex mixture, it is necessary to select appropriate masses or groups of masses at specific retention times for a particular compound or classes of compounds. The selection of such masses is the basis for a set a rules which then forms the data for the training set. There are no set procedures for such a selection process. The researcher must select appropriate masses for compounds of interest. For example, feed stocks and lube oils contain a wide range of compound types such as paraffins, cycloparaffins, aromatics, and olefins. It is known that paraffinic hydrocarbons yield fragment ions at masses 43, 57, 71, 85, . . . daltons, and these masses may be used as diagnostic of this class of compounds. Moreover, when coupled with retention time data, it is possible to identify concentrations of specific compounds within this class of compounds. In a similar manner, training sets for other chemical, performance, perceptual or physical properties may be developed by correlating compositional data with other properties of interest, e.g., distillation characteristics, viscosities and viscosity index, volatility, refractive index and the like. The result of a mathematical treatment such as PLS/PCR/RR of the training set is a set of coefficients for the properties of interest.
These coefficients are then multiplied by the data matrix for the sample. The result is a prediction of the desired property or properties which information can be used to control operating conditions in the distillation unit, extraction unit, dewaxing unit and/or hydrofiner as well as control lube oil quality. Moreover, variances in properties serve to promptly alert operators of upsets/changes in operating conditions which could influence product quality. The method of the invention is further illustrated by the following examples.
Example 1
The method for predicting the physical or chemical properties of a complex hydrocarbon mixture is demonstrated in this example using density at 15.degree. C. as the specific property for purposes of illustration. The method is generally applicable to a range of other physical properties as well as performance or perceptual or chemical properties of such mixtures.
The initial consideration is to establish a set of standard GC/MS operating parameters so that the GC/MS analytical data used for predicting properties is obtained under consistent operating conditions. The GC/MS instrument used in this example is a Hewlett-Packard 5972 Mass Selective Detector interfaced to a Hewlett-Packard 5890 Series II Gas Chromatograph fitted for used with microbore columns.
The GC/MS operating conditions are summarized in Table 1.
The analysis summarized in Table 4 is done for each reference sample. The results from these respective analyses form a training set which is subjected to mathematical treatment. The goal is develop a model that can be used to predict the unknown properties of future samples using their mass spectral data only. The mathematical treatments are described by multivariate correlation techniques such as Projection to Latent Structures (PLS) or otherwise known as Partial Least Squares (PLS), Principal Component Regression (PCR), and Ridge Regression (RR). These techniques are superior to ordinary multiple linear regression in their ability to treat collinearity amongst variables in the X-block or GC/MS data matrix (and Y-block or properties matrix for PLS), and in their ability to handle the quantity of data generated by the Chemist's Rules. Ordinary Multiple Linear Regression cannot be used to treat collinear variables.
PLS/PCR/RR are numerical analysis techniques for detecting and formulating a mathematical structure (model) within a data set comprising observations associated with multiple objects. Each object has associated with it observations for multiple variables, the latter being common to all objects. These multiple variables are assigned into two categories, known as X-block and Y-block (GC/MS data in this case). Observations associated with all variables in the X-block are realized from a common process. Observations associated with variables in the Y-block (known properties in this case) are realized from processes that may be different for each variable. The data set used to construct this mathematical model is referred to as the model calibration data set.
The common use of PLS/PCR/RR is to apply the model developed from the calibration data set to X-block observations realized for new objects (not in the calibration data set) to predict values for the corresponding variables in the Y-block for these new objects, without having to execute the Y-block processes used in the calibration data set. Using diagnostics that are simultaneously generated by the PLS/PCR/RR model, assessment of whether the new objects can be adequately described by the model, and whether the model is used in an extrapolation mode versus interpolation mode can be performed.
PLS/PCR addresses the collinearity features in the X-block and Y-block, by suitably reducing the dimensionality in both X- and Y-blocks (for PLS), and X-block only (for PCR) to form the model. Collinearity is a term referring to the existence of relationships between variables within the block itself. In the PLS modeling algorithm a number of independent dimensions in the X- and Y-blocks are identified by forming pseudo-variables known as principal components or latent vectors through different sets of linear combinations of original variables in each block. Each set of such combinations constitutes an independent dimension. It comprises a set of coefficients that each value associated with each variable in the block is to be weighted by to arrive at the new value for this dimension. The values for the new, reduced dimensions in the Y-block are regressed onto their counterparts in the new, reduced dimensions of the X-block to arrive at the most parsimonious dimension size (number of latent vectors) and their associated weights, with the final goal of one linear equation generated to permit prediction of Y-block variables using X-block variables. The number of dimensions used to construct the model is determined through optimization of a criterion known as PRESS (Prediction Error Sum of Squares), cumulated by a Cross Validation (CV) technique using the training data set, and, following the general model parsimony principle.
For PCR, the number of independent dimensions in the X-block are first selected and identified in a similar fashion as in PLS by forming principal components. Then, for each variable in the Y-block, a model is obtained by performing ordinary multiple linear regression using the Principal Components as "Prediction Variables".
For Ridge Regression, the collinearity problem is dealt with in a different manner than PLS/PCR. Here a diagonal matrix known as the Lambda Matrix is added to the Covariance Matrix of the X-block with the net effect of stabilizing the numerical computation needed to obtain the model coefficients. The selection of Lambda values is through optimization of PRESS criterion using cross validation of the training set.
Thus, the Chemist's Rule data for the various reference samples derived from GC/MS analysis form the X-block variables. PLS/PCR/RR treatment may require preliminary reorganization of the X-block data, such as transposition and removal of redundant data and constants or mathematical transformations. The Y-block variables are the property (or properties) to be predicted, and may also require mathematical transformations such as logarithmic or geometric, as well as reorganization. The data blocks may be represented by: ##EQU1##
The Y-block may be a single observation per set of Chemist's Rules as shown above, or it may be a n.times.m matrix of observations where there are m different properties to be predicted.
The results of PLS treatment of the training set data are a series of coefficients. Raw GC/MS data from an unknown sample (or samples) are then treated by the Chemist's Rules first to correct the retention times and then to form the ion summations. Each value for the Chemist's Rule ion summation result is then multiplied by the training set coefficients and summed to generate the prediction of the desired property. Table 5 illustrates the quality of the predicted density values for both the training set and the unknown test set for density at 15.degree. C.
Example 2
The procedure of Example 1 was repeated for predicting refractive index @75.degree. C. for the suite of 150N grade waxy raffinates. The same 22 samples were used for the training set, while the remaining five samples were used as a test set. The results for the refractive index at 75.degree. C. are summarized in Table 6.
Example 3
The procedure of Example 1 was repeated for predicting viscosity at 100.degree. C. for 150N grade waxy raffinates. The same 22 samples were used for the training set, while the remaining five samples were used as a test set. The results for density prediction are summarized in Table 7.
Example 4
The procedure of Example 1 was repeated for predicting weight percent wax content for 150N grade waxy raffinates. The same 22 samples were used for the training set, while the remaining five samples were used as a test set. The results for density prediction are summarized in Table 8.
Other properties of various grades of lubricating intermediate streams or finished products, such as, cloud point, flash point, viscosity index, pour point, basic nitrogen content, distillation and volatility characteristics, viscosity at 40.degree. C., aniline point, sulfur content, etc., could be predicted using the method according to the invention.
The subject method can also be applied for predicting properties of other types of complex hydrocarbon mixtures in similar boiling ranges, such as catalytic cracker feedstock quality, feed to extraction or dewaxing units, feedstocks to distillation or stripping towers, and to blending considerations.