EP1155357A1

EP1155357A1 - Method to control a lubricant production

Info

Publication number: EP1155357A1
Application number: EP00913273A
Authority: EP
Inventors: James William Gleeson; William Francis Heaney; Eugenio Sanchez; Viswanathan Visweswaran
Original assignee: ExxonMobil Oil Corp
Current assignee: ExxonMobil Oil Corp
Priority date: 1999-01-29
Filing date: 2000-01-27
Publication date: 2001-11-21
Also published as: AU3474600A; EP1155357A4; CA2359669A1; JP2002535481A; AU756690B2; WO2000045228A1; JP5524436B2; WO2000045228A9; JP2012246490A; RU2216573C2; KR100682553B1; KR20010089767A; JP5608712B2; CA2359669C

Abstract

A compositional model (801) is run on a computer or processor (803) having attached memory (805). The processor (803) may be PC-based or main frame-based and the memory (805) may be RAM, ROM or any appropriate storage device. The compositional model (801) typically receives as input data representing base stock(s) (807) (in a blending case), or data representing crude(s) (809) (in the blending case). The compositional model (801) may also receive as input data representing the refining data (810) of the particular crude(s) (809).In the case of crude(s) (809), such compositional model (801) parameters may be predicted from a library of industry data which is also stored in the memory (805). Once presented with such data, the compositional model (801) predicts whether the inputted base stock(s) (807) or crude(s) (809) are acceptable and provides an appropriate output (811).

Description

METHOD TO CONTROL A LUBRICANT PRODUCTION

The present invention is directed to the production of petroleum-based lubricants. More specifically, the present invention relates to a method for controlling lubricant quality and the lubricant production process by refining and blending by a method which includes a predictive technique relating the quality of the final lubricant product and its performance to compositional parameters.

Lubricating oils are, for the most part, based on petroleum fractions that boil above 232°C (450°F). The molecular weight of the hydrocarbon constituents is high and these constituents display almost all conceivable structures and structure types depending in large part on the type of crude oil from which they were prepared.

The rationale in lubricant refining is that a suitable crude oil, as shown by experience or by assay, can be refined into lubricant base stock having a predetermined set of properties such as, for example, appropriate viscosity, oxidation stability, and maintenance of fluidity at low temperatures. The refining process employed to isolate the lubricant base stock currently consists of a set of subtractive unit operations which remove the unwanted components. These unit operations include distillation, solvent refining, hydroprocessing, and dewaxing, each of which is basically a separation process.

A lubricant base stock (i.e. from a refined crude oil) may be used as a lubricant component, or it may be blended with another lubricant base stock having somewhat different properties. A particular base stock, prior to use as a lubricant, is conventionally compounded with one or more additives such as antioxidants, extreme pressure additives, and viscosity index (V.I.) improvers. As used herein, the term "stock," regardless of whether the term is further qualified, refers only to a hydrocarbon oil without additives. The term "solvent-refined stock" or "raffinate" refers to an oil that has been solvent extracted, for example with furfural. The term "dewaxed stock" refers to an oil which has been treated by any method to remove or otherwise convert the wax contained therein and thereby reduce its pour point. The term "waxy" as used herein refers to an oil of sufficient wax content to result in a pour point greater than 4°C (+25°F). The term "base stock" refers to an oil refined to a point suitable for some particular end use, such as for preparing an automotive oil, marine oil, hydraulic oil, etc.

The current practice for the preparation of high grade lubricating oil base stocks is illustrated in Fig. 1. The overall process is designated as 10. The first step is to vacuum distill an atmospheric tower residuum from an appropriate crude oil (step 100). This step provides one or more raw stocks within the boiling range of 700°F to 1000°F designated as light 15 (from 700° to 810°F), medium 20 (from 810° to 890°F), and heavy 30 neutrals (from 890° to 1000°F), and a vacuum residuum 40. Each stock is characterized by a different viscosity range, e.g., light neutral stocks have the lowest viscosity range (from 3.5 to 5.5 cS KV at 100°C) and heavy neutral stocks have the highest viscosity range (from 10.0 to 15.3 cS KV at 100°C). After preparation, each raw stock is extracted with a solvent, e.g., furfural, phenol or chlorex, which is selective for aromatic hydrocarbons, and which removes undesirable components (step 140). The vacuum residuum 140 usually requires an additional step to remove asphaltic material prior to solvent extraction (step 120). The raffinate from the solvent refining is generally very waxy and typically requires a dewaxing operation (step 160). Dewaxing raffinates is generally carried out by solvent dewaxing or catalytic dewaxing under conditions which produce a predetermined or target pour point for the base stock. The art of lube base stock production may further require an additional step of hydroprocessing the dewaxed stock (step 180). Hydroprocessing is the addition of hydrogen to the stock for the purpose of removing certain impurities that may be detrimental to lubricant performance, depending upon the final lubricant product specifications. Hydroprocessing is especially useful for removing sulfur as the hydrogen combines with sulfur to form hydrogen sulfide. Either the dewaxing step (step 160) or the hydroprocessing step (step 180), or both, can yield a stock of sufficient quality to be used as a base stock for final lubricant formulation.

Suitable materials, or additives, are then added to the base stock to augment the base stock's lubricating properties and to meet the required specifications, such as viscosity index (VI) (shown in Figure 6). Historically, lubricant manufacturing has been inflexible. Base stock produced from new crude sources is considered to be a "new" base stock and must be approved for use as a lubricant. The lube base stock approval system requires that a base stock produced from a new crude undergoes a costly and time-consuming approval process before it can be used as a lubricant oil. Full approval for samples produced from a new crude or base stock requires extensive bench and engine testing which generally takes 5 to 10 months to complete. Changes in process conditions also require re-approval with additional bench and engine testing. The high cost of bench and engine testing discourages optimization of process conditions on all but the most frequently run crudes. In addition, the need to streamline this system has become increasingly important as a result of the rapid changes in crude mixes from at least some geographic locations, quality concerns over purchased base stocks, and the drive for economic benefits from increased lube flexibility. Therefore, there is a need in the art for increased flexibility in lubricant production and to streamline the approval system for lube crudes, base stocks, and blends of these stocks.

According to the present invention, a lubricant base stock for use in manufacturing a lubricant for a particular application can be selected by a particular method. A candidate lubricant base stock is characterized according to a plurality of compositional components common to all lubricant base stocks and each compositional component varies in amount with respect to each base stock. A model predicting the performance of a lubricant base stock in a final lubricant product is formed by identifying particular ones of the plurality of compositional components in a lubricant base stock that are probative of lubricant performance for a plurality of applications, in combination.

An acceptable range is next determined for a combination of identified compositional components from a plurality of performance tests required of a plurality of lubricant products where the lubricant products each have a different viscosity. After the acceptable region is determined, adjustments are preferably made for viscosity variations of each lubricant base stock by predicting the amounts of the compositional components from at least one other acceptable range for a base stock of a different viscosity wherein the acceptable range is determined for a combination of the identified compositional components from the plurality of performance tests. A candidate base stock is then analyzed to determine the amounts of the identified components. It is then determined whether the amount of each identified component, in combination, in the candidate base stock is within the acceptable performance range.

The lubricant base stocks whose composition falls outside of the acceptable region predicted by the compositional model can also be blended to form a lubricant of satisfactory performance. This method begins with determining the acceptable region and adjusting for viscosity variations as described above. Next, two or more lubricant base stocks with complementary compositions are selected. The amount of each required to be combined to define the final lubricant composition within the acceptable region as predicted by the compositional model is calculated. The lubricant base stocks are then blended.

The predictive model can also be applied to the lubricant manufacturing process, making use of the effect of known refining parameters such as hydroprocessing conditions, solvent treat rate and temperature, on a given range of refinery stocks to be used in the manufacture of the lubricants. In this variant of the application of the predictive model, a lubricant base stock which is to be used in manufacturing a lubricant for a particular application from a refinery stream is characterized according to a plurality of compositional components common to the base stock and to the refinery stream. This method begins with determining the acceptable region and adjusting for viscosity variations as described above. Next, the amounts of the identified compositional components are identified for a lubricant base stock that, in combination, are required for acceptable lubricant performance and the amounts of identified components, in combination, in the refinery stream are determined to be outside of the acceptable range for such identified components. The refinery stream is then refined to produce a lubricant base stock having a combination of identified components within the acceptable range as determined by the model. By extension of this technique, combination of petroleum refinery streams which are to be refined into lubricant base stock can be blended into a refinery stream which can then be processed into the desired final stock or stocks. This method begins with determining the acceptable region and adjusting for viscosity variations as described above. Next, the amounts of the identified compositional components are identified for a lubricant base stock which is to be refined from the refinery streams and the composition of the blend is determined in accordance with the combination of identified components within the acceptable range as determined by the model. The streams are then blended in accordance with the determined composition and refined. In one variation, the refinery streams are crudes. In another variation of this embodiment, the refinery streams are selected from any one of the unit processes in the refinery. In the drawings:

Fig. 1 is a flow diagram showing an overview of a current lube base stock manufacturing process;

Fig. 2 is a flow diagram showing the method of the present invention; Fig. 3 is a graph showing how a compositional specification model is made for a plurality of lubricant performance tests according to one embodiment of the present invention;

Fig. 4 is a graph showing the overlap of the compositional specification models developed for light and heavy neutral base stocks according to one embodiment of the present invention; Fig. 5 is a graph showing a compositional specification model for light neutral base stocks, where such graph is used in connection with a method to blend a plurality of base stocks to yield a resulting blended base stock whose composition is within the acceptable region predicted by the compositional model;

Fig. 6 is a graph similar to the graph of Fig. 5, and shows the compositional specification model of Fig. 5 used with extraction severity data to obtain a base stock whose composition is within the acceptable region predicted by the compositional model; and

Fig. 7 is a block diagram illustrating the operation of an embodiment of the present invention. In the present invention, a computer model (801 in Fig. 7) is created and used to predict the quality of a lubricant base stock for use in a plurality of products having a plurality of viscosities. The model is based upon a quantitative analysis of key compositional parameters and performance criteria for a plurality of products representing a plurality of viscosities. Such model is herein referred to as comprehensive compositional model and, as will be shown, is normally and preferably comprised of compositional models developed for specific viscosities. The compositional models predict, among other things, whether the composition of a new lubricant base stock is similar or dissimilar to known lubricant base stock compositions that have previously passed or failed performance tests. New lubricant base stock compositions are, thus, evaluated in light of lubricant base stocks that have demonstrated acceptable performance in a plurality of products. Statistical analysis of base stock properties has shown that lubricant performance may be predicted from the chemical composition of the base stock. In particular, the various chemical compositional parameters are combined to obtain a comprehensive prediction of lubricant performance. Fig. 2 is a process flow diagram illustrating a method, to predict a base stock's lubricant performance in accordance with one embodiment of the present invention. The compositional parameters incorporated into the compositional model are chemical components of the base stock that have been identified as significantly contributing to lubrication performance (step 200). Examples of such components are, for example, total sulfur, aliphatic sulfur, basic nitrogen, aromatic distribution, nitrogen, aliphatic compound distribution, degree of aliphatic chain branching, molecular weight distribution or any combination thereof. The preferred embodiment includes total sulfur, aliphatic sulfur, basic nitrogen, total aromatics, aromatic ring distribution, or a combination of these parameters.

At a minimum, only two of the factors are used to create the compositional model to accurately predict a base stock's lubricant performance in a final lubricant formulation. All of the components are quantitatively determined according to standard analytical chemistry techniques. Total aromatics, for example, are determined from a chemometric correlation of the UV spectrum. Specifically, a Perkin-Elmer Lambda 19 spectrometer is used. The aromatic ring distribution is determined by conventional techniques known to those skilled in the art for classifying aromatic rings. Basic nitrogen is determined according to ASTM Method D2896. Total sulfur is determined by ASTM method D2622. The aliphatic sulfur content may be determined by UV spectrophotometry based on the quantitation of an iodine complex of the aliphatic sulfides according to the technique described by Drushel and Miller in Anal. Chem. 27, 495 (1955); and Anal Chem. 39, 1819 (1967). This method, which suitably measures absorbance of the complex at 310 nm, provides a reliable determination of aliphatic sulfur, being free of interferences by nitrogen compounds and the aromatic hydrocarbons indigenous to lubestocks.

Determining the compositional model's performance parameter (step 210) comprises evaluating the performance of a plurality of products having different viscosities. Variation of the compositional model's parameters as described below in accordance with variations in the viscosity of the lubricant products provides superior prediction of lubricant performance as compared to scaling compositional parameters in accordance with viscosity. The viscosities of the products are, thus, a function of the products' end use. For example, the compositional model's performance parameters include engine oil performance of automotive engine lubricants, marine engine performance of marine engine lubricants, industrial equipment performance of industrial equipment lubricants, and hydraulic performance of hydraulic lubricants. As each application requires a product of a different viscosity, the compositional model thereby incorporates a plurality of kinetic viscosities ranging from 3.5 to at least 20 cS (100°C).

Adjustments for viscosity variations may be made by simple linear interpolation derived from similar compositional models, having different determined parameters from different viscosities. It is also possible although less preferred to adjust for viscosity changes by scaling composition rather than model parameters for viscosity variations between samples processed to similar viscosity index from the same crude. The composition factors scale with viscosity according to a power law where the composition at a selected temperature and viscosity is related to the composition at a reference viscosity scaled by a power of the ratio between the viscosities at the selected temperature and the reference temperature; in this relationship, the exponent value varies typically from -2 to 2 (depending on measurement unit and possibly other empirically determined factors) for the compositional factors referred to above, with the value in many cases being in the range 0.1 to 1.0. The greatest sensitivity in the scaling relationship has been found to exist with respect to the basic nitrogen content and the least with total aromatics, polycyclic aromatics (two or more rings) and sulfur (total and aliphatic) occupying comparable intermediate rankings.

Minor variations in viscosity outside the range determined by the compositional model may also be extrapolated. It has been found satisfactory to generate compositional models for each lubricant type such as, for example, light and heavy neutrals, using a limited range of viscosity values rather than to use a single value of viscosity for each type. As will be shown, the comprehensive compositional model will incorporate the separate compositional models.

The limits of the performance parameters are determined by industry standard performance tests specific to the particular products and corresponding viscosities. For example, for products made from light neutrals, the following performance tests (accompanied by the particular application and acceptable performance level) may be used as examples of performance tests by which lubricant performance may be assessed. Other performance tests may be used according to standard requirements of the lubricant manufacturer.

Test Performance Level

Sequence IIIE (Automotive oils) ACEA A3-96

Sequence VE (Automotive Oils) ACEA A3-96

Mack T-8 (Automotive oils) Mack EOL TOST life (Hydraulic oils) >1500 hr

For products made from heavy neutrals, the following performance tests

(accompanied by the particular application and acceptable performance level) may be used as examples of performance tests by which lubricant performance may be assessed. Other performance tests may be used according to standard requirements of the lubricant manufacturer:

Test Performance Level

Cat 1-G2/1M-PC (Marine oils) CF

Cat 1M-PC (Automotive oils) CF

TOST life (Hydraulic oils) >1000 hr Such tests are known to those skilled in the art of performance testing.

For the purpose of making the compositional model, a lubricant product, when tested according to any one of the above tests, is designated as either "pass" or "fail." Next, an acceptable region in a multi-dimensional space, indicating good lubricant performance, is determined by defining the performance boundaries (i.e., "pass" or "fail") with respect to key compositional parameters. This is accomplished over a broad range of product applications and viscosities. Such acceptable region is then represented by the compositional model. The end result is that the compositional model allows lubricant performance to be predicted over all possible base stock compositions and viscosities.

To develop a comprehensive compositional model (step 220), separate compositional models are first developed for specific viscosities. The model parameters for a test viscosity are found by interpolating between the parameters of the model at the nearest higher viscosity and the parameters of the model at the nearest viscosity lower than the test viscosity. The following describes how to develop the comprehensive compositional model by first developing separate compositional models for a light neutral viscosity (kinematic viscosity at 100°C of 4,5 cSt) and a heavy neutral viscosity (kinematic viscosity at 100°C of 12.5 cSt). Small variations in the viscosity of the samples tested from these viscosities, such as within the viscosity grade, were accounted for by scaling the composition according to typical trends of composition with viscosity. These can be readily determined by those skilled in the art by fitting the trends in composition with viscosities that are produced from the same crude:

Composition Scaled = (composition at y KV 100°C) x (y/reference viscosity)^Λz wherein ^Λ is "taken to the power of;" the reference kinematic viscosity at 100°C is 4.5 or 12.5 cSt; and z is typically between 0 and 1.

In view of the above, developing the acceptable region (step 220) for the individual compositional models are the same for base stocks of all viscosities. A series of base stocks are performance tested for a plurality of product formulations that employ base stocks of the same viscosity range. Preferably, each of the lubricant product formulations are lubricant products for different applications, i.e., automotive, hydraulic, industrial, etc. As herein described, regions are drawn around the points that are a "pass" for each formulation, where such regions exclude the "failed" base stocks. A method employed to determine the acceptable regions, i.e., the regions drawn around the points that are a "pass" for each formulation, is preferably a technique involving cluster analysis of the performance test data. One such technique is, for example, Soft Independent Classification Analysis (SIMCA) as described by Muhammad A. Sharaf, Deborah L. Illman, and Bruce R. Kowalski's Chemometrics, John Wiley & Sons, Inc., New York 1986. SIMCA models were constructed to model base stock compositions that have either passed or failed engine tests. The SIMCA models first depend, however, on various feature or variable selection techniques used to screen the base stock properties. The screening techniques determine the optimal inputs for a classification model. Two univariate strategies, variance and Fischer weights, are employed to select the variables that discriminate between passing and failing engine tests.

The variance weight is the ratio of inter category variance and the sum of the intra-category variances for a given property. It approximates the distance between two distributions. For example the total variance for the percent aromatic fraction of the base stock is calculated for passing Caterpillar 1-G engine tests and is divided by the sum of the variance for the passing tests. Overall variance weights for n>2 categories are the geometric mean of the individual category variance weights. The variance weights are given by the formula:

where:

w = Variance weight for two categories; x = Base stock compositional parameter, i.e. sulfur, aliphatic sulfur, basic nitrogen, or aromatic distribution;

N = Number of samples in a given class

P or F = Denotes the Pass or Fail data The following is an example of variance weights for two base stock clusters:

Total 2+ Basic Sulfur Aliphatic Thiophen 2+-

Aromati Aromati Nitrog Sulfur es Thio

CS CS en

weight 7.72 10.90 2.08 11.44 4.99 15.05 3.71

(p/

Larger numbers indicate a greater likelihood that the category discriminates between the passing and failing tests. In the current example, the largest number for thiophenes indicates that it is the most significant. In this way, one can choose to use the most significant categories first in the analysis.

A Fischer weight uses the distance between mean values in the two distributions divided by the sum of variances to determine the discriminating power of the variables. Overall, the Fischer weights for n>2 categories are the mean of the individual category fisher weights. The fisher weight is given by:

w_k(P,F) = \ X_P - Xp \ (A -2)

(1/ N_P) Σ(X_P- x + (1/ N_F) Σ (χ_F - χ_Ff

The SIMCA method uses principal component analysis to construct a model for each class, i.e. pass/fail. Factor analysis is used to calculate the significant chemical patterns for each class. The significant patterns are known as the principal components. The principal component analysis provides a convenient method for data compression. It also provides a rotation of the data (base stock composition) to an orthonormal basis, removing any co-linearities in the data. The principal components form a new set of axes for the data (base stock properties).

Two key statistical tests, the F-test and the Mahalanobis Distance, are used to evaluate the compositional models. The F-test is a measure of the residual error which describes the distance between a sample and the space defined by the principal components. The Mahalanobis Distance describes the distance from the compositional centroid of the compositional model. The F-test determines if the error associated with a sample is statistical error or whether it is due to the inability of the principal components to model the data. The Mahalanobis distance indicates whether the prediction is an extrapolation of the compositional model.

The compositional models are used by projecting unknown samples into each model. The chemical patterns for each model are applied to the new base stock. The residual error and the distance from the centroid of each class form a basis for determining the similarity of unknown samples. The unknown samples may belong to one or several classes.

Soft models, i.e., the local clusters of data independently modeled and independent models, are used to accurately predict or classify samples. For example, base stock compositions that have demonstrated passing performance in the Caterpillar 1-G2, for example, can be considered as one cluster; failing base stock compositions are therefore another cluster. The passing samples will form the basis for one model and the failing samples form the basis for a second model. SIMCA combines cluster analysis with principal components regression. This technique involves supervised learning; the samples used to build the compositional models must be classified. SIMCA uses principal component analysis (PCA) to calculate the chemical patterns for each class. The PCA provides a statistical basis for selecting the significant chemical patterns. The PCA provides a technique for re-expressing the original data. The original data matrix, the base stock properties, Xι, is decomposed to a linear combination of scores, Li (weights) and eigenvector patterns Vj^T: Xi = LiVi^T (A-3)

Several techniques are available for decomposing a matrix. One technique for decomposing a matrix is the Singular Value Decomposition, SVD as described by Golub, Gene H., and Charles Van Loan, in Matrix Computations, 2^nd Edition, p.70, The Johns Hopkins University Press, Baltimore (1989):

Xι = UiSiVi^T (A-4)

The SVD provides two orthogonal matrices U and V that can be used to build the SIMCA models. The diagonal of the S matrix contains the singular values, which are the positive square roots of the eigenvalues of the matrix X^TX. The singular values provide a direct indication of the rank of the matrix X. The number of significant components used is determined by the cumulative variance calculated from the singular values. The matrix V contains the eigenvectors of X with the eigenvectors used as the chemical patterns for samples in the data matrix X. The matrix multiplication of UjSj gives the scores, Lj.

The terms "fitted" and "predicted" will now be used to identify various kinds of base stocks. These terms are not interchangeable. The term "fitted" denotes that a base stock sample was used to build the compositional model. "Predicted" indicates that a base stock sample is a true prediction and was not used in building the compositional model. All of the equations that follow are based on matrix calculations to facilitate the translation to Matlab. Matlab is commercial computation software available from Math Works, Inc. (Natick, MA) and provides a implementation of the singular value decomposition, SVD for matrices.

The SIMCA model requires that the residual error of each class be calculated. The principal components analysis provides some error filtering of the original data. Therefore, the residual of each base stock property is also examined. Residual error gives an indication of the "tightness" of each compositional box. The SIMCA model uses the measured base stock properties and their chemical patterns, or eigenvectors, to calculate the residual error of each sample. The residual error can be estimated using linear regression techniques. The fitted residual error is defined as the difference between the observed values Xj and the fitted, i.e., projected, values e = X, -χ, (A -6)

XX HX, (A-7)

The fitted values of X can be calculated from:

The matrix H is commonly referred to as the hat matrix and transforms observed responses to fitted responses. The H matrix is estimated from

= X(χ^τ Xf χ^τ (A -8)

as described by Weisberg, Sanford, Applied Linear Regression, 2"^tf Edition, p.47, John Wiley & Sons, Inc., New York, 1985. The fitted residual error can be re- expressed as

e = (l-H)χ, (A -9)

The orthogonal matrix U, truncated down to its most significant principle components, provides the basis for an important projection associated with the SVD:

UiU,^T = projection on to range (X) (A-10)

The projection onto the range of X provides a method for transforming observed responses to fitted responses. The projection onto the range of X from the SVD of data matrix X provides a principal component based method for estimating the hat matrix:

The residual error of the fitted data is then estimated from the equation:

= (l-U,U^r,)X, (A- 12)

The residual variance of each class is calculated from the equation:

NS NV -₂ s, Σ (A- 13)

J'lΣ(NS-NC-1)(NV-NC)

where NS is the number of samples, NV is the number of variables (base stock properties), and NC is the number of principal components. The residual error of an unknown sample, Xu projected into class "i" is estimated from the significant eigenvectors and singular values:

and e = χ_u-u_uU,^τX, (A- 15)

The residual variance of the unknown sample is estimated by

e = X.-X,V,S?U^T,X, (A- 16)

B, = V,S-;U X, ( '¹⁷)

e-X -B,) (A- 18)

An F-test is used to determine if the sample is outside the compositional box (volume):

F - si (A -20)

Sf

As was discussed above, a second important statistic used in the SIMCA models is the Mahalanobis distance as described by Gordon, A.D., Classification,

Chapman and Hall, 1981, herein incorporated by reference. The Mahalanobis distance (M-distance) indicates whether a projection into the compositional model space is an extrapolation. The M-distance is estimated by

The covariance matrix l~ is an unbiased estimate of the covariance matrix and can be calculated from:

_{Γ =} X_X± (A - 22)

NS -1

The covariance matrix can be re-expressed in terms of the SVD of Xi:

_r - V>S1V' (A - 23)

NS -J The Mahalanobis distance for fitted samples becomes

V = (NS-1)(V,S²,V,^T)^'' (A -24)

where:

(V,S,²Fi^T)-¹ = (V.V (S,²)-¹(V,)-¹ (A-25)

T,' = (NS-l)(V,S V,SUPT) (A- 27)

Δ„ = J(NS - l)xώagonal(χ, V, Sf V,^τ X, ) (A -28)

Δ„ = sj(NS-l)xdιagonal(τj,U^τ,) (A - 29)

and similarly the Mahalanobis (M) distance for predicted samples becomes Δ„ = J(NS - l)xdiagonal(H) (A - 30)

Thus, a plot of the M distance versus the calculated F value provides a graphical representation of the base stock composition. The maximum M distance is a direct measure of the size of the principal component space. Base stock samples should not exceed the M distance calculated from the data used to construct the compositional model. The F statistic indicates the type of residuals that remain after the sample is projected into the principal component space. Low values of F, for a model developed with a statistically significant number of base oil samples, F<4,>95% confidence limits indicate the base stock property residuals are null. The predictive errors are due to normal errors associated with collecting analytical data. If the F value is greater than 4, the principal components do not accurately model the data, and the sample does not fit into the compositional model.

To arrive at the comprehensive compositional model, separate compositional models for both light and heavy neutrals are constructed as described above. The compositional models for light and heavy neutrals are developed by repeating the above process for each of the respective base stock's performance tests such as those listed above.

Fig. 3 is a simplified illustration of how the acceptable region for each viscosity grade is generally developed. In Fig. 3, the performance of a series of lubricant base stocks of the same viscosity grade were measured in two performance tests. The performance tests could be performed on either the same lubricant product formulation or on product formulations for different applications, i.e., automotive oil lubricants, industrial oil lubricants, marine oil lubricants, etc.

Still referring to Fig. 3, the "pass" regions, x and z, for each test are plotted on one graph. For simplicity, x and y represent the "pass" regions for only two performance tests of products made by either light or heavy neutral base stock. The overlapping "pass" regions become the acceptable region defining the respective compositional model for light or heavy neutral base stocks and is, thus, predictive of lubricant performance across a wide range of lubricant products.

Adjustments for viscosity variations (step 230) may be made at this point by, for example, simple linear interpolation derived from similar compositional models having different determined parameters for different viscosities. Minor variations in viscosity outside the range determined by the compositional model may also be extrapolated. It has been found satisfactory to generate compositional models for each lubricant type such as, for example, light and heavy neutrals, using a limited range of viscosity values rather than to use a single value of viscosity for each type with adjustments for minor variations in the viscosities of the samples within each viscosity grade made by scaling composition according to typical trends of compositional variation with viscosity as described above. Thus, once the compositional model for the respective light and heavy neutral base stocks are developed, the comprehensive compositional model allows the accurate interpolation of the acceptable performance regions of base stocks of intermediate viscosity from the light and heavy neutral models. The interpolation is performed by another manipulation.

In order to interpolate the compositional model, the following compositional modeling principle is employed: a base stock of intermediate viscosity is approved if it is possible to mimic its composition as a blend of two acceptable light and heavy neutrals.

In one embodiment of the present invention, the compositional model employs ellipses or ellipsoids to define the acceptable region, however, other geometrical figures can also be used. A simple linear interpolation of the principal component ends of the "ellipsoid" designated as the acceptable performance region, for example, could be used such that at every viscosity a different ellipse is used. Linear interpolations of orthogonal axes, however, are not orthogonal, and the desirable properties of the principal components would be lost.

Thus, the interpolations are accomplished by "morphing" one ellipsoid into another, i.e., linearly translating the center of the ellipse as a function of viscosity, while gradually rotating the ellipsoid, and expanding/contracting the ellipsoid axes according to their singular-value scales. This can be achieved if the compositional models are in the same "space", and they have the same number of principal components.

In matrix notation, there is a square non-symmetric rotation transformation matrix T, that applied to a matrix V₀ of mean-centered ortho-normal principal components, will rotate it exactly into the orthonormal principal components of another model, Vi:

TVc^ V_! (B1)

Since V₀ is ortho-normal, T is simply expressed as

TvoVo = T = Wo . (B2) Intermediate axes between model 0 and model 1 can be expressed as fractional rotations of T. Fortunately, well-developed matrix algebra exists for fractional transformations. If we define z=0 for model 0 (light neutral), and z=1 for model 1 (heavy neutral), any intermediate set of rotated axes V₂ can be expressed as:

V_z = T^z Vo = E D^z E Vo (B3) where T^z is the z-matrix power of T, computed from Ts eigen-analysis T E - E D, where E is ortho-normal and D is diagonal, and D^z is also diagonal, with the diagonal elements of D raised to the z power. Since T is not symmetric, both E and D will contain conjugate pairs of complex eigenvectors-eigenvalues, requiring complex algebra. Additional complications need to be addressed before computing V_z for a given Vi and V₀. One is chirality, and the other is alignment. Chiral objects are mirror images of each other. Therefore, it is impossible to rotate a chiral object into its mirror image. Likewise, a set of orthogonal vectors in multidimensional space has a certain characteristic or "sense" that is analogous to chirality. The only valid rotations that do not drastically deform the vectors into imaginary planes in intermediate steps are those that transform between sets of vectors with the same chirality, or sense. Thus, to overcome this problem, the Vi columns are multiplied by -1 until the intermediary transformation vectors V₀.₅ are in the real plane.

The last complication of alignment is solved by computing all possible rotations in the real plane, and comparing their rotation "distance" for the 1^st, 2^nd, . . ., N ^h principal components, until the rotation with the minimal angular change as a function of z change is found. Once the proper rotation is found, the rest of the compositional model parameters are interpolated accordingly. The average, axis scaling, and Mahalanobis distance are interpolated linearly, while the f-ratio and the residual variance are interpolated quadratically. This technique can be extended to a viscosity grade higher or lower than heavy or light neutrals, respectively. For example, by developing a compositional model for bright stocks (i.e. - VI of approximately 30 centastokes), one skilled in the art can readily interpolate for viscosity ranges between heavy neutrals and bright stocks. Prediction of product performance for products of all viscosities are made and compared to the actual performance at those viscosities. If the accuracy of the predictions at viscosities between those for which the compositional models were developed is less than that at the viscosities for which the compositional models were developed, additional compositional models can be developed at more viscosities until the accuracy is satisfactory. Once the comprehensive compositional model is determined across a plurality of products, for example, PVL, CVL, marine, and industrial lubricants, the comprehensive compositional model necessarily incorporates the performance levels of each finished product formulation. Thus, any lubricant base stock, selected according to any variety of ways, having any viscosity, whose lubricant performance is unknown in a finished lubricant product formulation, merely has to be analyzed for the chemical concentrations of the key compositional parameters (step 240). Once evaluated against the compositional model, the lubricant base stock either falls into the acceptable region or it does not. Thus, the compositional model predicts whether the lubricant base stock is acceptable for use in all lubricant products or not acceptable for use as a lubricant, at least in its present form.

Fig. 4 shows the compositional model comprised of the acceptable compositional regions for the light and heavy neutral models. While there is considerable overlap between the acceptable regions for the light and heavy neutral compositional models, this is not the most efficient way to use the compositional model. Rather, the most efficient way to use the compositional model is in its comprehensive form as such comprehensive compositional model employs all reasonable interpolations for viscosity in between these extremes. Since both light and heavy neutrals are produced from a given crude, the compositional model, as shown in Fig. 4, indicates that many potential crudes are suitable to be used for lubricant products over a plurality of viscosities, as long as there is overlap in the acceptable areas as determined by the compositional model for the plurality of viscosities. The lower sulfur and aromatics region, however, is acceptable for light neutrals but not acceptable for heavy neutrals. As shown, such a region generally provides good Sequence VE, Mack T-8, and CEO Caterpillar 1M-PC performance.

As described above, the compositional model is based on lubricant base stock compositions known to yield good lubricants. In other embodiments of the invention, options are available for those base stock compositions that have been designated by the compositional model as "failed." For example, in one embodiment of the invention, the "failed" sample may be blended with a base stock that has

"passed" to yield an acceptable composition. This embodiment is illustrated in Fig. 5.

Fig. 5 shows an approximate two-dimensional representation of the acceptable region of the compositional model for a light neutral with the concentration of total aromatics versus total sulfur as the compositional parameters. The light neutral base stocks of a plurality of crudes are plotted on the graph according to the relative concentrations of total aromatics and total sulfur. Some of the crudes produce light neutrals whose compositional components fall within the acceptable region of the compositional model, and some do not. Those that do not fall within the acceptable region, however, may be blended with a base stock composition that does fall within the acceptable region. Base stock compositions can be blended only if the compositions between the two or more base stocks to be blended (step 260) are complementary. For example, in Fig. 5, the light neutral base stock produced from crude J contains a high amount of aromatics and a high concentration of total sulfur. Such base stock produced from crude J, could be blended with the light neutral base stock produced from crude D, for example, which contains a lower level of both total aromatics and sulfur such that the combined product is well within the acceptable region predicted by the compositional model. Knowing the concentrations of the compositional parameters allows a calculation of the amounts of each base stock required to achieve a resultant base stock composition whose compositional parameters are within the acceptable region (steps 270 and 280). Two "failing" base stocks can also be combined in this manner as long as their compositions are complementary in such a way that their combination will result in a base stock whose compositional parameters fall within the compositional model's predicted acceptable region. Referring again to Fig. 5, the light neutral base stock produced from crude A, for example, can be combined with the light neutral base stock from crude X to form the "passing" light neutral base stock AX. In contrast, the light neutral base stock produced from crude A cannot be combined with the light neutral base stock from crude Y to yield a "passing" combination as the bias in composition of each base stock is not favorable enough to achieve a "passing" resultant base stock.

The compositional model allows the integration of known refining data with the prediction of lubricant performance. For example, the compositions of the base stocks produced from the world's crudes can be readily determined by conventional analytical techniques by those skilled in the art. Thus, knowing either a selected crude or selected base stock's refining data allows one to determine whether the selected crude or base stock can be refined to within the acceptable region predicted by the compositional model (step 285). In the case of a crude, knowing the selected crude's refining data allows one to determine whether the selected crude can be refined to produce a base stock whose composition is within the acceptable region predicted by the compositional model. For example, it is known in the art of lube manufacturing that as the solvent, e.g. furfural, extraction severity is increased (step 140 in Fig. 1), certain key compositional components, such as, for example, total aromatics and aliphatic sulfur, decrease when a base stock is solvent extracted.

For this extraction step, any commercial refinery extraction process can be employed but normally, however, the extraction process is a continuous flow process. In most continuous flow commercial processes, the oil is introduced into the lower part of the vessel with a multiplicity of mixers inside. The multiplicity of mixers allows the refinery to control the severity of the extraction. A solvent is introduced near the top of the vessel. The solvent is typically liquid sulfur dioxide mixed with benzene, furfural, propane, phenol, N-methylpyrrolidone, or the like. A furfural extraction is preferred. The solvent then works its way towards the bottom of the vessel dissolving the extract as it goes along. The remaining solvent-refined stock, or raffinate, rises to the top of the vessel and is separated from the solvent-containing extract.

In addition, hydroprocessing (step 180 in Fig. 1) may also be used to remove sulfur and other heteroatoms from the base stock. Any commercial refinery hydroprocessing method can be employed. Most commercial hydroprocessing operations entail mixing the stream of base stock with hydrogen at temperatures from 260° to 425°C. The oil combined with the hydrogen is then charged to a vessel filled with a catalyst so that the required reactions can occur. The oil is then sent to a flash tank to remove inorganic heteroatoms such as hydrogen sulfide and any light hydrocarbons that may have resulted during the catalytic reaction.

The effect of the commercial extraction processes on, for example, light neutral base stocks is shown relative to the light neutral compositional model with the concentration of total aromatics versus total sulfur as the compositional parameters in Fig. 6. To determine whether a base stock can be refined to yield a composition within the acceptable region, a series of base stocks produced from different crudes is plotted to illustrate the effects of extraction severity on the base stocks. The effects of the extraction processes on the viscosity index (VI) is also shown. Here, it is evident that if a base stock produced from a certain crude has a high concentration of, for example, total aromatics and total sulfur such that it is outside of the acceptable region predicted by the compositional model, such as, for example, the base stock produced by crude K, the base stock may be further refined (step 290) to a point where its composition does fall within the acceptable region predicted by the compositional model. Likewise, the extracted base stocks can be combined in the manner described above to yield a composition within the acceptable region predicted by the compositional model.

The effect of refining of a number of different of crudes can enable a prediction of blend parameters to be made, based on the composition of the base stocks produced from such crudes. Likewise, knowing the composition and viscosity of any stream in the refinery, such as, for example, those shown in Fig. 1 , a composition can be predicted and blended that will further refine into a base stock having a composition within the acceptable region determined by the model. The compositional changes that depend on changes in refining, either solvent or hydroprocessing, are well known and can easily be modeled for the purposes of selecting suitable refinery streams from known crude sources and for selecting suitable crudes themselves. In particular, the effects of unit processes such as, for example, distillation, solvent extraction, dewaxing, and hydroprocessing can be modeled and used according to the present understanding. This modeling includes such elements as aromatics level and type, sulfur level and types, nitrogen level and types, and may include others as required. The effects of unit processes should include such parameters as distillation cut points, solvent dosage and temperatures, contact times, dewaxing dosages and temperature profiles, hydrogen pressure, temperature, catalyst type, conversion and hydrogen treat rate.

Referring now to Fig. 7, it is seen that in one embodiment of the present invention, the compositional model 801 developed as discussed above is run on a computer or processor 803 having attached memory 805. The processor 803 and memory 805 may be any appropriate processor and memory, respectively, without departing from the spirit and scope of the present invention. For example, the processor may be PC-based or main-frame-based and the memory may be RAM, ROM, or any appropriate storage device such as a hard drive and the like.

Still referring to Fig. 7, the compositional model typically receives as input data representing a base stock 807 or base stocks 807 (in the blending case), or data representing a crude 809 or crudes 809 (in a blending case). Regarding the data representing the crude(s) 809, the model may also receive as input data representing the refining data 810 of the particular crude(s) 809. Typically, the compositional parameters are, in fact, input in the base stock(s) case, or are predicted in the case of crude(s). Typically, appropriate information on the compositional parameters is available as assay data which is stored in the memory 805 in the case of the base stock(s). In the case of crudes, such compositional parameters may be predicted from a library of industry data which is also stored in the memory 805. Once presented with such data, the compositional model predicts whether the inputted base stock(s) 807 or crude(s) 809 are acceptable in accordance with the method as described above and provides an appropriate output 811.

Claims

CLAIMS;

1. A method for predicting the lubricant performance of a lubricant base stock for use in manufacturing a lubricant for a particular application, the base stock being characterizable according to a plurality of compositional components common to a plurality of lubricant base stocks, each compositional component varying in amount with respect to the base stocks, the method comprising the steps of: identifying particular ones of the plurality of compositional components in a lubricant base stock that, in combination, are probative of lubricant performance for a plurality of applications; determining an acceptable range for a combination of the identified compositional components at two or more viscosities from a plurality of performance tests on a plurality of lubricant products having different viscosities, to define a compositional model at each of the viscosities; determining the amounts of the identified compositional components in a candidate lubricant base stock to derive an acceptable range of compositional components, in combination, at the viscosity of the candidate lubricant base stock; and determining whether the amount of each identified component in the candidate base stock, in combination, is within the acceptable range for the identified components.

2. The method of claim 1 wherein the determined acceptable range for each identified component is selected from the group consisting of an acceptable range that is independent of any other identified component and an acceptable range that is inter-related to an acceptable range of another identified component.

3. The method of claim 1 wherein the identifying step comprises selecting compositional components from the group consisting of basic nitrogen, sulfur, aliphatic sulfur, aromatics, and a combination thereof.

4. The method of claim 1 wherein an acceptable range for a combination of the identified compositional components is determined at two or more viscosities from a plurality of performance tests on a plurality of lubricant products having different viscosities, to define a compositional model at each of the viscosities; and an adjustment is made for viscosity variations for lubricant base stocks whose viscosities differ from that of the selected viscosities by adjustment of the parameters of the compositional model from at least two compositional models at values of viscosity differing from those of a candidate lubricant base stock, to derive an acceptable range of compositional components, in combination, at the viscosity of the candidate lubricant base stock

5. A method according to claim 4 in which the adjustment for viscosity variations is made by interpolation of compositional model parameters between compositional model parameters for a compositional model at a viscosity nearest above the viscosity of the candidate lubricant base stock and the compositional model parameters at the viscosity nearest below the viscosity of the candidate lubricant base stock.

6. A method according to claim 5 in which the compositional models for at least one of the compositional models at the next higher viscosity value and the next lower viscosity value are generated with variations about a viscosity grade with adjustments for viscosity variations within the viscosity grade made by scaling compositional variations with viscosity.

7. A method according to claim 1 in which a plurality of lubricant base stocks are selected from among a plurality of lubricant base stocks for use in manufacturing a lubricant for a particular by determining a required amount of each of the plurality of selected base stocks to combine to define the combination to be used in a final lubricant composition that will perform within the acceptable range for each identified component; and blending the required amount of each selected base stock required to form the final lubricant composition.

8. A method according to claim 1 in which the lubricant base stock is selected from a number of lubricant base stocks in accordance with the determination made by the method.

9. A method according to claim 1 in which the lubricant base stock is selected from a binary blend of two base stocks, neither of which has a composition which is acceptable for the selected base stock for the particular application.

10. A method for blending a combination of petroleum refinery streams from among a plurality of such refinery streams for manufacturing a lubricant base stock for a particular lubricant application, each of the streams being capable of being useful for lubricant base stocks characterizable according to a plurality of compositional components common to a plurality of the base stocks, each compositional component varying in amount with respect to each base stock, the method comprising the steps of: identifying particular ones of the plurality of compositional components in a base stock that, in combination, are probative of lubricant performance for the particular application; determining an acceptable range for a combination of the identified compositional components at at least two viscosities from a plurality of performance tests on a plurality of lubricant products of different viscosities, thereby defining a compositional model at each of the at least two viscosities; deriving an acceptable range of compositional components, in combination, at the viscosity of the candidate lubricant base stock; identifying the amounts of the identified compositional components for a lubricant base stock which is to be refined from the plurality of refinery stream; determining the composition of the blend of the streams in accordance with the combination of identified components within the acceptable range for such identified components for the particular application, and blending the streams in accordance with the determined composition.

11. A method according to claim 10 an adjustment is made for viscosity variations for lubricant base stocks whose viscosities differ from that of the two viscosities by adjustment of the parameters of the compositional model from at least two compositional models at values of viscosity differing from those of the candidate lubricant base stock.

12. A method according to claim 10 in which a plurality of lubricant base stocks are blended to make a lubricant base stock for a particular application, where none of the plurality of base stocks is acceptable as the selected lubricant base stock, the method including determining a composition for a blend of selected base stocks such that the composition of the blended combination of base stocks with respect to the identified components is within the acceptable range for such identified components for the base stock for the particular application, and blending the base stocks in accordance with the composition determined for the blend.

13. A method according to claim 12 in which the blend is a binary blend of two base stocks, neither of which has a composition which is acceptable for the selected base stock for the particular application.

14. A method according to claim 13 in which the compositions of the base stocks which are not acceptable for the selected base stock for the particular application are complementary for at least one of the identified components.

15. A method according to claim 14 in which at least one of the base stocks which is not acceptable for the selected base stock has a composition in which the amount of at least one of the identified components is below the acceptable range for that component in the selected base stock and the amount of the same at least one component in the base stock having the complementary composition is above the acceptable range for that component.

16. A method according to claim 11 wherein the adjustment for viscosity variations is made by interpolation of compositional model parameters between the compositional model parameters for the compositional model at the viscosity next higher than the viscosity of the candidate lubricant base stock and the compositional model parameters at the viscosity next below the viscosity of the candidate lubricant base stock.

17. A method according to claim 16 wherein the adjustment for viscosity variations is made by interpolation of compositional model parameters between the compositional model parameters for the compositional model at the viscosity next higher than the viscosity of the candidate lubricant base stock and the compositional model parameters at the viscosity next below the viscosity of the candidate lubricant base stock.

18. A method according to claim 16 wherein the adjustment for viscosity variations is made by extrapolation of compositional model parameters from compositional model parameters at viscosities differing from the viscosity of the candidate lubricant base stock.

19. The method of claim 10 wherein the determined acceptable range for each identified component is selected from the group consisting of an acceptable range that is independent of any other identified component and an acceptable range that is inter-related to an acceptable range of another identified component.

20. A method according to claim 13 in which the determination of the blend of composition of the blend of streams is made in accordance with the composition of the streams and a refining processing model which correlates stream composition with the composition of the lubricant base stock according to the identified components.

21. A method for making a lubricant base stock to be used in manufacturing a lubricant for a particular application from a petroleum refinery stream by refining, the lubricant base stock and the refinery stream being characterizable according to a plurality of compositional components common to the base stock and to the refinery stream, each compositional component varying in amount with respect to different lubricant base stocks and to different refinery streams, the method comprising the steps of: identifying particular ones of the plurality of compositional components in a base stock that, in combination, are probative of lubricant performance for the particular application; determining an acceptable range for a combination of the identified compositional components at two or more viscosities from a plurality of performance tests on a plurality of lubricant products, each lubricant product having a viscosity, the plurality of lubricant products having a plurality of different viscosities, thereby defining a compositional model at each of the two viscosities; determining the amounts of the identified components which are required in the base stock at the viscosity selected for the lubricant base stock for acceptable performance in the particular application; determining the amounts of the identified components in the refinery stream, in combination, are not within the acceptable range for such identified components; and refining the refinery stream to produce the lubricant base stock having the identified components in combination in the compositional ranges for acceptable performance of the lubricant base stock in the particular application.

22. A method according to claim 21 which includes adjusting for viscosity variations for lubricant base stocks whose viscosities differ from that of the two viscosities by adjustment of the compositional model's parameters from two compositional models at values of viscosity differing from those of a candidate lubricant base stock to derive an acceptable range of compositional components, in combination, at the viscosity of the candidate lubricant base stock.

23. The method of claim 21 wherein the determined acceptable range for each identified component is selected from the group consisting of an acceptable range that is independent of any other identified component and an acceptable range that is inter-related to an acceptable range of another identified component.

24. The method of claim 21 wherein the identifying step comprises selecting compositional components from the group consisting of basic nitrogen, sulfur, aliphatic sulfur, aromatics, and combinations of these components.

25. A method according to claim 21 in which the refinery stream is refined in accordance with the composition of the stream of the composition required of the lubricant base stock and a refining processing model which correlates stream composition with the composition of the lubricant base stock according to the identified components of the refinery stream and the base stock.

26. A method according to claim 21 in which the refinery streams are crude oils or gas oils.

27. A method according to claim 25 in which the crude oil is subjected to a refining process model comprised of unit subtractive refining which accounts for compositional changes that depend on changes in the unit subtractive refining operations.

28. A method according to claim 21 wherein the adjustment for viscosity variations is made by interpolation of compositional model parameters between the compositional model parameters for the compositional model at the viscosity next higher than the viscosity of the candidate lubricant base stock and the compositional model parameters at the viscosity next below the viscosity of the candidate lubricant base stock.

29. A method according to claim 28 in which the compositional models for at least one of the compositional models at the next higher viscosity value and the next lower viscosity value are generated with variations about a viscosity grade with adjustments for viscosity variations within the viscosity grade made by scaling compositional variations with viscosity.

30. A method according to claim 21 wherein the adjustment for viscosity variations is made by extrapolation of compositional model parameters from compositional model parameters at viscosities differing from the viscosity of the candidate lubricant base stock.

31. A method according to claim 21 in which none of the refining streams is capable of being refined under acceptable refining conditions into the selected lubricant base stock, and a composition for the blended refinery streams is determined such that the blended combination will refine into the lubricant base stock whose amount of each identified component, in combination, is within the acceptable range for such identified components for the particular application, and the streams are blended in accordance with the determined composition.

32. A method according to claim 30 in which the determination of the composition of the blend of streams is made in accordance with the composition of the streams and a refining processing model which correlates stream composition with the composition of the selected lubricant base stock according to the identified components.

33. The use in connection with a method for selecting a base stock from among a plurality of base stocks for use in manufacturing a lubricant for a particular application, each base stock being characterizeable according to a plurality of compositional components common to all of the base stocks, each of which varies in amount with respect to each base stock, of a computer-readable medium having computer-implementable instructions for performing the steps of: identifying particular ones of the plurality of compositional components in a base stock that are probative of lubricant performance for a plurality of applications; determining an acceptable range for a combination of the identified compositional components at two or more viscosities from a plurality of performance tests on a plurality of lubricant products, each lubricant product having a viscosity, the plurality of lubricant products having a plurality of different viscosities, thereby defining a compositional model at each of the two viscosities; adjusting for viscosity variations for lubricant base stocks whose viscosities differ from that of the two viscosities to derive an acceptable range of compositional components, in combination, at the viscosity of the candidate lubricant base stock; determining the amounts of the identified compositional components in a candidate base stock; and determining whether the amount of each identified component in the candidate base stock is within the acceptable range for such identified component.

34. The use according to claim 33 of a computer-readable medium having computer-implementable instructions for performing the steps of adjusting for viscosity variations of the lubricant base stocks whose viscosities differ from that of the at least two viscosities by adjustment of the parameters of the compositional model from at least two compositional models at values of viscosity differing from those of the candidate lubricant base stock to derive an acceptable range of compositional components, in combination, at the viscosity of the candidate lubricant base stock.

35. The use according to claim 34 of a computer-readable medium having computer-implementable instructions for performing the further steps of: selecting one of the base stocks as a candidate base stock; determining the amounts of the identified composition components in a candidate base stock; and determining whether the amount of each identified component in candidate stock is within the acceptable range for the identified component.