Method and apparatus for determining a property or yield of a hydrocarbon product based on NIRA of the feedstock

Abstract

The properties and yield of a product from a hydrocarbon conversion or separation process may be determined from the NIR (Near Infrared) spectrum of the feedstock. This is accomplished by: (a) determining with an IR spectrometer the absorbance of the feedstock at a certain number of frequencies in the spectral range 16667 to 3840 cm.sup.-1 starting from a defined base line, and then (b) determining the property or yield by applying a correlation between the property or yield and the absorbance values. The correlation is determined experimentally by multivariate regression and is dependent upon the type of spectrometer employed, the property or yield to be determined, and the frequencies used.

Description

This invention relates to a method for the determination of the properties of the feedstock, the properties of the product and/or the yield of the product in a hydrocarbon conversion or separation process by carrying out NIR (near infra red) spectroscopic analyses of the feedstock and correlating these with the desired properties. The method is suitable for use on-line and in real time and for incorporation in process control.

Callis et al, Analytical Chemistry, Vol 59, No. 9, pp 624A-636A (May 1987) mentions the possibility of determining the octane rating of an unleaded gasoline by NIR spectroscopy and point out in this particular case, the existence of a connection between other properties of the products and the NIR spectra of the products.

However, when a product is produced by converting or separating various components, which are often themselves mixtures, major problems are encountered in predicting the properties of the final products. This arises from the fact that several properties do not follow linear laws.

For certain processes, models exist which enable the prediction of products and yields. To use these, however, it is usually necessary to know the properties of the feed in detail, the chemical composition, the molecular weights of the components etc. This information is normally obtained by laboratory analyses which involves a delay in obtaining the results.

According to the present invention there is provided a method for the determination of the properties of the product and/or the yield of the product in a hydrocarbon conversion or separation process wherein the method comprises the following steps:

(a) determining with an IR spectrometer the absorbance of the feedstock at a certain number of frequencies in the spectral range 16667 to 3840 cm.sup.-1 (0.6 to 2.6 microns), preferably 12500 to 3840 cm.sup.-1 (0.8 to 2.6 micron), most preferably 4760 to 4000 cm.sup.-1 (2.1 to 2.5 micron) starting from a defined base line, and (b) determining the property by applying a correlation between the property and the absorbance values, the correlation being determined experimentally by multivariate regression, depending on the type of spectrometer used, the property required, and the frequencies used.

The frequencies used may be selected from the following:

4670 cm.sup.-1

4640

4615

4585

4485

4405

4390

4385

4332

4305

4260

4210

4170

4135

4100

4070

4060

4040

The corresponding frequencies expressed in statutory units (Hz) are obtained by multiplying these values by 3.times.10.sup.10, the velocity of light in cm/s.

This choice is neither exhaustive nor exclusive. The choice of other frequencies will not alter the method but will require the use of other coefficients in the models enabling the calculation of the desired properties sought from the spectra.

The base line (regarded as corresponding to zero absorbance) is preferably 4780 cm.sup.-1.

BRIEF DESCRIPTION OF THE DRAWINGS

The FIGURE is a block diagram of a process and analytical instrumentation for implementing the invention.

DETAILED DESCRIPTION OF THE INVENTION

The processing equipment 1 may be computer controlled by a feed-back control system for altering process conditions in response to variations in the property and/or yield of the product from the desired value, from the determination of the NIR spectrum of the feedstock.

The spectrometer 2 may be linked to a signal processing device 3 to permit numerical treatment of the spectrum, preferably by Fourier transformation.

The spectrometer suitably has a resolution power of 4 cm.sup.-1.

The method may be used on-line and in real time.

The correlation may contain linear terms, quadratic terms and homographic (ratio) terms.

According to a specific embodiment, the hydrocarbon conversion process may be a catalytic reforming process and the properties may be selected from the following: yields of gasoline (C.sub.5.sup.+), hydrogen and gas, octane numbers (clear or leaded), vapor pressure and density.

The correlation in this case may be either linear for a feed of relatively constant composition, for example the KUOP characteristic varying between 11.8 and 12.2 (terms p of equation (1) below), or can involve quadratic terms (terms q of equation (1)) or homographic (ratio) terms (terms r of equation (1)) for a feed of varying composition, the KUOP characteristic varying between 10.8 and 12.2.

The KUOP characteristic or Watson factor is expressed by the formula: ##EQU1## Rankine scale=absolute Fahrenheit scale Density 60/60=density of the product at 60.degree. F. compared to water at 60.degree. F. ##EQU2## D.sub.i . . . D.sub.j are absorbances at the wavelengths in question.

The constant C and the various coefficients p, q and r are obtained from multivariate digital analysis techniques applied to a set of feedstocks used for preliminary calibration.

In the computer-assisted operation of a reformer, the NIR spectrum of the feed as it enters the unit is taken in real time and is treated as an information vector which continuously monitors the potential of the feed.

The abundance of the NIR spectrum and the experimental accuracy obtained, resulting from the accumulation of spectral data by fast Fourier transformation, means that this information is both accurate and relevant to the chemical operations involved in the reforming process. The NIR spectrum is therefore a numerical indicator of the suitability of the feed for reforming.

Furthermore, it has been observed that the differential spectrum, when measured in real time, constitutes a digital signal enabling computer assisted operation of the unit to be improved. Catalytic reforming results in an increase in octane number which is related to the severity of the operating conditions for a given state of activity of the catalyst.

It has been noted that the differential spectrum can be directly correlated with the gain in octane number (motor and/or research) of the clear or leaded product at different levels of tetraethyl or tetramethyl lead and is related to the severity of the process.

It is therefore possible to have a sensitive and accurate method of computer assisted control for a catalytic reformer by continuously monitoring the differential spectrum, enabling parameters to be adjusted to produce gasoline of the specified octane number, regardless of variation in the quality of the feed.

With this new method of operating the unit, variations in the differential spectrum in the NIR sensitively and rapidly optimise processing conditions using a computer programme.

Thus, according to a modification of the above embodiment, the hydrocarbon conversion process may be a catalytic reforming process wherein the method comprises the further steps of:

(c) carrying out differential spectrum measurements in the same range and at the same frequencies on the feedstock and the product, (d) applying a correlation between the differential spectrum and the octane number gain, the correlation being determined experimentally by multivariate regression analysis, (e) comparing the sum of the two values thus determined (octane number and octane number gain) with a desired objective, and (f) modifying the operating conditions of the catalytic reformer until the determined octane number gain is in agreement with the desired octane number gain.

In catalytic reforming the frequencies used may be selected from the following:

4670 cm.sup.-1

4485

4332

4305

4210

4100

4060

4040

The correlation between the differential spectrum and the octane number gain may be a linear correlation, e.g. ##EQU3## where .alpha.O and .alpha.i are constants and DDi is the input/output differential spectrum for wavelength i.

According to another embodiment, the hydrocarbon conversion process may be a selective hydrogenation process and the property may be the maleic anhydride number of the product.

The feedstock may be a C.sub.5 -200.degree. C. fraction derived from steam cracking called steam cracked gasoline.

The process is not restricted to this type of product, however, similar gasolines being derived from other types of cracking reaction.

It is possible to determine on-line and in real time the level of diolefins in steam cracked gasoline selectively hydrogenated, starting with the NIR spectrum of the steam cracked gasoline as it enters the unit. This is done under given operating conditions and for a catalyst of defined composition and activity.

The process based on the invention is not restricted to this particular property, however, it is possible to predict other properties of the product, for example:

the content of olefins the clear octane numbers and those for leaded gasolines, different levels of tetraethyl lead or tetramethyl lead

Analysis by spectrometer in the NIR range makes it possible to calculate by differing numerical methods several properties of a hydrocarbon fraction from its spectrum.

According to a further embodiment, the hydrocarbon conversion process may be a catalytic hydrocracking process and the properties may be selected from the following: total conversion of the feedstock to product and conversion to gasoline.

In this case, the frequencies used are preferably selected from the following:

4585 cm.sup.-1

4390

4332

4260

4210

4170

4100

4070

The sum of the absorbances for the eight wavelengths considered is reduced to 100 and thereafter the relative percentage of absorption for a given band will be used. Hence, if D.sub.i is the absorbance for the wavelength i, then: ##EQU4## and the percentage of relative absorption for the frequency i: p.sub.i will be equal to: ##EQU5##

The feedstock may be vacuum distillates, atmospheric residues, extracts derived from lube oil processing, hydrogenated residues, deasphalted residues, cracked distillates and cat cracker recycle oils.

According to yet another embodiment, the hydrocarbon conversion process may be a catalytic cracking process and the properties are selected from the following:

(i) the conversion of the feed to gasoline and the yield of residue in the product, (ii) the density, research octane number (RON) and motor octane number (MON) of clear and leaded gasolines from the cracking process.

For catalytic cracking, the frequencies may be selected from the following:

4670 cm.sup.-1

4615

4485

4135

4100

4060

The feed may be a distillate, for example, vacuum distillate gas oils, an atmospheric residue (deasphalted or not), a lube oil extract, a hydrogenated or deasphalted residue, a cracked distillate or a cat cracker recycle oil.

Some of the properties which may be determined include:

the research octane number (RON) and the motor octane number (MON) of the gasoline produced, whether clear or leaded, and at different levels of lead, the cetane number of gas oil, the densities of various fractions produced (gasoline, gas oil, residue), the viscosity of the residue, etc.

By combination with computer control of a fluidised bed catalytic cracker, the NIR spectrum allows complete and accurate characterisation of the feed which, together with the processing conditions of the unit (for example, inlet and outlet temperature of the heater, C/O (ratio of the catalyst flow rate to the feed flow rate), space velocity, nature and activity of the catalyst, etc) leads to an efficient, computer model.

It is thus possible to control the unit with great precision.

According to another embodiment, the hydrocarbon conversion process may be a visbreaking process and the properties may be selected from the following: the yield of the gas oil product, the viscosity, specific gravity and xylene equivalent of the residue, and the conversion limit, i.e, the maximum conversion before an unstable residue is produced.

The feedstock may be atmospheric or vacuum distillate residues. They may be used in a blend and may have various petroleum fractions added to them: e.g. lube oil extracts, hydrogenated or deasphalted residues, cracked distillate units, catalytic cracker recycled oils, etc.

In addition to the hydrocarbon conversion processes exemplified above, the method is also applicable to the control of separation processes such as distillation.

Thus according to another embodiment, the hydrocarbon separation process may be a distillation process and the properties may be selected from the following:

(i) the quantitative composition of a mixture of crude oils constituting the feedstock when the distill and is known to be a blend of crude oils, (ii) the yields of products obtained, (iii) the densities and viscosities of the various products, (iv) the cetane number, aniline point, cloud point and the KUOP characterisation factor of the gas oil, (v) the n-paraffinic, isoparaffinic, cycloparaffinic and aromatic content of the gasoline and the research octane number (RON) and motor octane number (MON) of clear and leaded gasoline with different contents of tetraethyl or tetramethyl lead.

In computer assisted operation of a processing unit, the NIR spectrum of the feedstock entering the unit is thus taken in real time and treated as an information vector which characterises the processing properties of the feedstock in the process on a continuous basis. The detail in the NIR spectrum and the experimental precision which may be obtained from the accumulation of spectral data ensures that this information is reliable and relevant to processing operations.

The analysis and data processing time is less than 1 minute.

According to another aspect of the present invention there is provided apparatus for carrying out a method for the determination of the properties of the product and/or the yield of the product in a hydrocarbon conversion or separation process as hereinbefore described, the apparatus comprising an infrared spectrometer linked to a computer programmed in such manner that the property may be determined continuously and in real time.

The apparatus may comprise a feedback control system for computer control 4 of the processing equipment in response to variations of the property and/or yield of the product from the desired value, from the determination of the NIR spectrun of the feedstock.

On-line analysis may use fibre optics as a means for transmitting light signals or may involve physical transferance of a sample to the cell of the spectrometer.

The invention is illustrated but not limited with reference to the following examples:

EXAMPLE 1

Catalytic Reforming

Operating conditions of the reformer were as follows:

reactor temperature: 510.degree. C. LHSV: 1.85 catalyst: R.sub.12 UOP pressure: 39 bar

The NIR spectrum of the feed at the entry to the reformer was determined at four frequencies, F.sub.1, F.sub.2, F.sub.3 and F.sub.4.

The product yield and RON of the product were then predicted using the following equations (5) and (6):

RON clear=2151.36 D.sub.5 +131.22 D.sub.14 +280.67 D.sub.9 -243.157 (5) Yield C.sub.5.sup.+ =-114.84 D.sub.14 +227.18 D.sub.16 +338.81 D.sub.9 -249.3 % vol (6)

D.sub.i represents the absorbance of the NIR spectrum at frequency F.sub.i, the zero base being taken as 4780 cm.sup.-1.

The predictions were then compared with the actual results obtained by conventional means.

The feedstock was a highly paraffinic naphtha having a KUOP factor of 12.2, a 90% ASTM distillation point of 135.degree. C. and a density d.sup.15 of 0.686.

NIR spectrum: ______________________________________Frequency Absorbance______________________________________F.sub.1 4485 cm.sup.-1 D.sub.1 0.0347F.sub.2 4100 D.sub.2 0.55594F.sub.3 4060 D.sub.3 0.61466F.sub.4 4305 D.sub.4 0.65502______________________________________ Yield % vol RON Clear______________________________________Predicted properties from 81.2 88.3equations (5) and (6)Actual properties by 80.4 88.0experiment______________________________________

EXAMPLE 2

Catalytic Reforming

Example 1 was repeated using a diferent feedstock.

This was a moderately paraffinic naphtha having a KUOP factor of 11.8, a 90% ASTM distillation point of 148.degree. C., and a density d.sup.15 of 0.7421.

Absorbance values at the same frequencies as Example 1 were:

______________________________________ D.sub.1 0.0364 D.sub.2 0.52127 D.sub.3 0.60595 D.sub.4 0.66828______________________________________ Yield % vol RON Clear______________________________________Predicted properties from 88.3 91.1equations (5) and (6)Actual properties by 87.6 91.4experiment______________________________________

EXAMPLE 3

Catalytic Reforming

(1) The octane number of the naphtha feedstock to the reformer was obtained from its NIR spectrum at the following frequencies:

______________________________________Frequency Absorbance______________________________________F.sub.1 4670 cm.sup.-1 D.sub.1 0.01223F.sub.2 4485 D.sub.2 0.03878F.sub.3 4332 D.sub.3 0.9339F.sub.4 4100 D.sub.4 0.6167F.sub.5 4040 D.sub.5 0.6083F.sub.6 4305 D.sub.6 0.6794F.sub.7 4210 D.sub.7 0.4800______________________________________

The research octane number values were obtained after adding 0.4 g/l tetraethyl lead.

Operating conditions of the reformer were as follows:

reactor temperature 510.degree. C. LHSV 1.85 catalyst R.sub.12 UOP pressure 39 bar

As before, the octane number ON of the product from the feed is predicted by equation (1) in the form: ##EQU6##

The constant C and the various coefficients p,q and r are obtained from multivariate digital analysis techniques applied to a set of feedstocks used for preliminary calibration.

Equation (1) becomes:

(RON-0.4)=-65.52-3.178 D.sub.4.D.sub.6 +180.3 D.sub.2 -6.32 D.sub.3 -1814 D.sub.5.D.sub.6 +1232.D.sub.5 +1165.D.sub.6 -109 D.sub.7

The RON calculated by means of Equation (1) is 83.2 compared with an experimental value of 83.5.

A similar procedure could be followed for other clear or leaded octane numbers.

(2) The octane gain is estimated by Equation (2): ##EQU7## where .alpha.O and .alpha.i are constants and D D.sub.i is the input/output differential spectrum for the frequency i.

The values of the constants are as follows:

.alpha.O=25.32

.alpha.1=136.48

.alpha.2=-268.14

.alpha.3=51.15

.alpha.4=138.756

.alpha.5=-68.40

.alpha.6=61.85

.alpha.7=-288.12

The absorbance of the NIR spectrum of the feedstock and product were as follows for the frequencies investigated:

Under these conditions Equation (2) gives an octane increase of 13.5.

(3) In view of the quality of the feedstock and the product specification an octane increase of 15 is required.

Consequently the reforming temperature must be progressively increased until the required gain is achieved. This operation is carried out by increasing the temperature by 4.degree. C.

In this operation the differential spectrum serves as a guide in increasing the temperature. This type of adjustment based on the differential spectrum remains valid even if the nature of the feedstock changes during adjustment of the temperature.

EXAMPLE 4

Selective Hydrogenation

The operating conditions of a gasoline selective hydrogenation unit were as follows:

temperature entering the reactor: 70.degree. C. temperature emerging from the reactor: 170.degree. C. through-put of hydrogen: 45 Nm.sup.3 /m.sup.3 of batch V.V.h=4 total pressure: 15 bars palladium catalyst: Procatalyse LD 265.

The unit was fed with a steam cracked gasoline having a maleic anhydride number of 40.2. The hydrogenated gasoline has a maleic anhydride number of 7.

______________________________________ Spectrum of Spectrum at Outlet Feedstock from reactorFrequencies cm.sup.-1 Absorbance Absorbance______________________________________4670 D.sub.1 0.18186 0.184764615 D.sub.2 0.14555 0.148884585 D.sub.3 0.16501 0.174284485 D.sub.4 0.11747 0.072424332 D.sub.5 0.56587 0.581864260 D.sub.6 0.46701 0.49344170 D.sub.7 0.29995 0.301844100 D.sub.8 0.35576 0.342644060 D.sub.9 0.95571 0.988124040 D.sub.10 0.50597 0.505674640 D.sub.11 0.14633 0.154064385 D.sub.12 0.41282 0.404664305 D.sub.13 0.48511 0.502434210 D.sub.14 0.29911 0.298364135 D.sub.15 0.28855 0.28183______________________________________

Using the spectrum of the feedstock, the value of the maleic anhydride number of the hydrogenated gasoline emerging from the reactor was determined by means of the following relation:

Maleic anhydride number=-325.05-1744.6.D.sub.1 +649.D.sub.2 -565.5.D.sub.3 +813.1.D.sub.4 +349.3.D.sub.5 -356.9.D.sub.6 +845.8.D.sub.7 +528.9.D.sub.8 +285.7.D.sub.9 +116.1.D.sub.10 +124.3.6.D.sub.11 +445.6.D.sub.12 -333.9.D.sub.13 +182.2.D.sub.14 -1773.4.D.sub.15

The calculated value is 5.6 as against 7 by measurement.

EXAMPLE 5

Selective Hydrogenation

A similar exercise was carried out on the feedstock to the selective hydrogenation unit of Example 4.

The maleic anhydride number values were measured on the feedstock, using the feedstock spectrum and on the hydrogenated gasoline product using the feedstock spectrum according to the following relationship:

maleic anhydride number=-283.05-704.2.D.sub.1 +633.9.D.sub.2 -552.4.D.sub.3 +794.2.D.sub.4 +341.2.D.sub.5 -348.6.D.sub.6 +826.2.D.sub.7 +516.7.D.sub.8 +279.1.D.sub.9 +113.4.D.sub.10 +1214.8.D.sub.11 +435.25.D.sub.12 -326.1.D.sub.13 +178.1.D.sub.14 -1732.3.D.sub.15

The following results are obtained:

______________________________________ Calculated Measured______________________________________Maleic anhydride number on entry 40 40.2Maleic anhydride number at exit 8 7.8______________________________________

This example shows that it is also possible to determine the value of the maleic anhydride number directly, using the product spectrum.

EXAMPLE 6

Hydrocracking

The feedstock was a distillate derived from crude petroleum from the Cameroons, KOLE.

The relative absorbance spectrum as a percentage of the total absorbance is as follows:

D.sub.4 18.43%

D.sub.5 13.06%

D.sub.6 13.98%

D.sub.7 12.12%

D.sub.8 11.87%

The KUOP factor of the feedstock was obtained from the following equation:

KUOP=0.618 D.sub.4 +0.36 D.sub.5 +1.07 D.sub.6 -0.248 D.sub.7 +1.032 D.sub.8 -28.5

The KUOP factor as measured was 12.31 compared with a calculated value of 12.2.

EXAMPLE 7

Hydrocracking

A heavy aromatic distillate of KUOP factor 10.9 was used.

The relative absorbance spectrum of this feedstock was as follows:

D.sub.4 16.93%

D.sub.5 12.76%

D.sub.6 12.77%

D.sub.7 10.94%

D.sub.8 12.50%

The value for Conradson carbon residue may be determined by the following equation:

C % weight=36.84-3.77.D.sub.4 -2.42.D.sub.5 +8.444 D.sub.6 -6.473.D.sub.7 +2.023.D.sub.8

The value determined in this way is 4.5% by weight as against 5.5% by weight by experiment. This agreement is satisfactory, taking into account experimental accuracy.

Other correlations can be obtained for the nitrogen content and density of the feed.

EXAMPLE 8

Hydrocracking

The operating conditions used in this example were as follows:

Reactor - entry temperature 400.degree. C. Reactor - exit temperature 430.degree. C. Total pressure 155 bars Hydrogen flow rate 1300 N m.sup.3 /m.sup.3 Average temperature 423.degree. C. Catalyst: molybdenum cobalt on aluminium silicate.

The following correlations ccan be established for these operating conditions: ##EQU8##

A, B, C and D are constants linked to the operating conditions of the reactor.

For the conditions described above, the values of these constants are:

A=123.12 B=277.14 C=35.893 D=8.55

The unit is fed with a distillate, the characteristics of which are as follows:

Density: 0.952 KUOP factor: 11.5 nitrogen content 1735 ppm Conradson carbon residue: 0.42% by weight Saturated hydrocarbons: 37.20% by weight Aromatic hydrocarbons: 61.30% by weight

The absorption spectrum of the feedstock gives the following relative absorbance factors:

D.sub.8 =12.2 D.sub.4 =17.75

The application of equations (3) and (4) gives the following results by comparison with experimental figures:

______________________________________ Calculated Measured______________________________________Conversion into gasoline % by weight 16.1 15Total conversion 78 79.5______________________________________

The reactor operating in the same conditions is fed with a different batch whose characteristics are as follows:

Density 0.920 KUOP factor 11.7 Nitrogen content 1300 ppm Carbon residue 0.26% by weight

The relative absorbance factors are as follows:

D.sub.8 =12.19 D.sub.4 =18.05

The equations (7) and (8) give the following results:

______________________________________ Calculated Measured______________________________________Conversion into gasoline % by weight 15.7 16.5Total conversion 75.9 75______________________________________

EXAMPLE 9

Catalytic Cracking

In this and the following Examples 10, 11 and 12 operating conditions were as follows:

Temperature at entry to riser: 250.degree. C. Temperature at exit from riser: 525.degree. C. Mass flow rate: 78 kg/h per kg C/O: 6.6 Catalyst MAT activity ASTM-D 3907: 65

The feedstock had a KUOP factor of 11.95 and had the following absorption spectrum:

F.sub.1 =4615 cm.sup.-1 D.sub.1 =0.02395 F.sub.2 =4135 D.sub.2 =0.55726 F.sub.3 =4100 D.sub.3 =0.5193 F.sub.4 =4485 D.sub.4 =0.04015

The gasoline yield had the following characteristics: ASTM distillation: initial boiling point 38.degree. C., 90% distilled at 190.degree. C., being 47.7% by weight.

The yield of residue following the ASTM distillation provided 10% distilled at 385.degree. C., being 9% by weight.

The yields can be obtained from the feedstock spectrum by the following relationships, for the control conditions for the unit as given below: ##EQU9##

The results obtained are 47.27% for the yield of gasoline and 9.39% for the yield of residue, using these proportions.

EXAMPLE 10

Catalytic Cracking

A feedstock having a KUOP factor of 11.6 was used with the following spectrum measured under the same conditions as in Example 9:

______________________________________F.sub.1 = 4615 cm.sup.-1 D.sub.1 = 0.3357F.sub.2 = 4135 D.sub.2 = 0.52092F.sub.3 = 4100 D.sub.3 = 0.48834F.sub.4 = 4485 D.sub.4 = 0.04989______________________________________

The observed and calculated yields using the proportions as in Example 9 are:

______________________________________ Observed Calculated______________________________________Yield of gasoline % by weight 40.6 41.25Yield of residue % by weight 15 14.47______________________________________

EXAMPLE 11

Catalytic Cracking

A feed of KUOP factor=11.95 was used. This had the following absorption spectrum:

______________________________________F.sub.1 = 4670 cm.sup.-1 D.sub.1 = 0.00977F.sub.2 = 4485 D.sub.2 = 0.04015F.sub.3 = 4100 D.sub.3 = 0.5193F.sub.4 = 4060 D.sub.4 = 0.52665______________________________________

The following feed properties are obtained by the following relationships:

Density of the feed=0.96497+4.1D.sub.1 -0.1935D.sub.3 Conradson carbon residue of the feed as % by weight=-7.744+16.07D.sub.3 +48.175D.sub.1 Aniline point of feed in .degree.C.=70-500 (D.sub.4 /D.sub.3)-1.066

The experimentally determined values and those calculated with the aid of the above correlations are as follows:

______________________________________ Calculated Found Experimentally______________________________________Density at 15.degree. C. 0.905 0.908% weight carbon residue 1.07 1.10Aniline point .degree.C. 96 94.4______________________________________

EXAMPLE 12

Catalytic Cracking

A feed of KUOP factor=11.6 was used. The gasoline obtained by cracking had the following characteristics:

Distillation: initial boiling point 38.degree. C., 90% distilled at 190.degree. C. Density at 15.degree. C.: 0.749 Research octane number after addition of lead tetraethyl at 0.4 g/l: 92.4

The absorption spectrum for the feed is measured under the following conditions:

______________________________________Frequencies Absorbance______________________________________F.sub.1 = 4670 cm.sup.-1 D.sub.1 = 0.0155F.sub.2 = 4485 D.sub.2 = 0.05079F.sub.3 = 4100 D.sub.3 = 0.48692F.sub.4 = 4060 D.sub.4 = 0.51066______________________________________

The characteristics of the gasoline produced are obtained from the spectrum of the feed from the following relationships:

Density of the gasoline=0.7553-0.0355 log.sub.10 (D.sub.3 /6.2D.sub.2)

Research octane number after addition of lead tetraethyl at 0.4 g/l (RON 0.4)=93.6-4.82 log.sub.10 (D.sub.3 /6.2D.sub.2)

The calculated values are d=0.7486 and octane number 92.7 for experimentally determined values of d=0.749 and octane number 92.4.

The process used therefore accurately predicts the experimental results obtained from the unit operated under the conditions established above.

EXAMPLE 13

Visbreaking

A feed was used consisting of a vacuum distillation residue obtained from a crude of Kuwait origin distilling above 565.degree. C.

The product temperature on leaving the visbreaker was 481.4.degree. C.

The following were determined from the NIR spectrum of the unit feed:

the visbreaking conversion, the xylene equivalent of the visbreaking residue, the gasoline yield.

The absorbances were determined at the following frequencies from the feed spectrum:

______________________________________Frequency (cm.sup.-1) Absorbance______________________________________4670 D.sub.1 0.009644640 D.sub.2 0.023014615 D.sub.3 0.029614585 D.sub.4 0.035664485 D.sub.5 0.053144385 D.sub.6 0.34854332 D.sub.7 0.69714305 D.sub.8 0.47104260 D.sub.9 0.56924210 D.sub.10 0.40194170 D.sub.11 0.39744135 D.sub.12 0.35464100 D.sub.13 0.31854060 D.sub.14 0.33814040 D.sub.15 0.26954405 D.sub.16 0.255______________________________________

The xylene equivalent of the visbroken residue was determined by the following equation which links the absorbances and the reaction temperature T:

Xylene equivalent=-1553.9+1.849 T-162 D.sub.1 +6824.8 (D.sub.3)+1027.9 (D.sub.5)+2480.9 (D.sub.6)+598.3 (D.sub.7)-1717.2 (D.sub.10)+7368.3 (D.sub.11)-8313 (D.sub.14)-573.4 (D.sub.15)-262.6 (D.sub.16)

The experimentally determined xylene equivalent had a value of 79 as against a calculated value of 78.9.

The gasoline yield (C.sub.5 cut - 185.degree. C.) was obtained by the following equation:

Gasoline yield=2.73+1166.2 (D.sub.5)+46.59 (D.sub.8)-253.03 (D.sub.13)

The residue yield (cut point above 350.degree. C.) was obtained by the following equation:

Residue yield=20.46+2934.4 (D.sub.1)-3905.6 (D.sub.2)+173.39 (D.sub.7)

The gasoline and residue yields determined by the equations above are respectively 6.07 and 79.75% weight.

A comparison with the values noted is shown below:

______________________________________ Calculated Actual______________________________________Xylene equivalent 79 78.9Gasoline yield % weight 6.07 7.7Residue yield % weight 79.75 78______________________________________

EXAMPLE 14

Visbreaking

The reaction temperature which will give an xylene equivalent of 70 and obtain the corresponding gas oil yield was determined.

The temperature may be obtained by the following equation which links the xylene equivalent and the reaction temperature RT in .degree.C.

RT=840.4+0.54 xylene equivalent+87.6 (D.sub.1)-3691.1 (D.sub.3)-555.9 (D.sub.5)-1341.75 (D.sub.6)-323.6 (D.sub.7)+928.7 (D.sub.10)-3985 (D.sub.11)+4495.9 (D.sub.14)+310.1 (D.sub.15)+142 (D.sub.16)

The unit used operates under the same conditions as Example 13 and was fed with a vacuum distillation residue of Brent origin, the NIR spectrum of which has the following characteristics:

______________________________________Frequency (cm.sup.-1) Absorbance______________________________________4670 D.sub.1 0.018984640 D.sub.2 0.033184615 D.sub.3 0.038894585 D.sub.4 0.045434485 D.sub.5 0.064894385 D.sub.6 0.37894332 D.sub.7 0.76124305 D.sub.8 0.52164260 D.sub.9 0.65014210 D.sub.10 0.46134170 D.sub.11 0.46464135 D.sub.12 0.42654100 D.sub.13 0.3894060 D.sub.14 0.40414040 D.sub.15 0.3324405 D.sub.16 0.2688______________________________________

The temperature determined by means of the above equation and the spectrum above was 480.4.degree. C.

The gas oil yield was obtained by the following equation, linking the reaction temperature in .degree.C. and the absorbances:

Gas oil yield, % weight=-290.45+0.34 (RT)+626.5 (D.sub.6)-234.34 (D.sub.14)

The gas oil yield obtained for the calculated temperature of 480.4.degree. C. was 15.6% weight.

The actual and predicted results are given below:

______________________________________ Predicted Actual______________________________________Xylene equivalent 70 60RT .degree.C. 480.4 480Gas oil yield, % weight 15.6 14.8______________________________________

EXAMPLE 15

Visbreaking

The conversion limits of the Kuwait and Brent feedstocks of Examples 13 and 14 were determined from the spectra given in Examples 13 and 14 under the operating conditions of Example 13.

The conversion limit value was determined by the following equation:

CL=-53.05+917.83 (D.sub.5)+1903.94 (D.sub.12)-2090.1 (D.sub.13)+55.86 (D.sub.16)

The actual and calculated values were as follows:

______________________________________ Calculated Actual______________________________________Brent 20.5 21.1Kuwait 19.4 19.9______________________________________

EXAMPLE 16

Distillation

In this Example and the following Examples 17-23, the experimental distillation procedure was carried out on a distillation column having an efficiency of 14 theoretical plates.

The NIR spectra of three crude oils was determined. The values of the absorbances for these three crude oils are as follows at the frequencies used, the base line being taken as 4870 cm.sup.-1 :

______________________________________ Absorbances D.sub.1 MixtureCrude Oil Iranian Heavy Zakum Djeno 50-35-15______________________________________FrequenciesF.sub.i cm.sup.-1F.sub.3 4615 0.02736 0.0248 0.01758 0.0251F.sub.4 4585 0.02671 0.02562 0.01895 0.0255F.sub.5 4485 0.03811 0.03981 0.02172 0.0379F.sub.7 4332 0.93621 0.96899 0.94434 0.9494F.sub.9 4260 0.77561 0.79598 0.80701 0.7876F.sub.11 4170 0.55478 0.56599 0.57862 0.5622F.sub.13 4100 0.49855 0.50971 0.51222 0.5050F.sub.14 4060 0.53475 0.54304 0.52624 0.5361F.sub.15 4040 0.41833 0.42156 0.41788 0.4196______________________________________

A homogeneous mixture of the three crude oils was produced in the following proportions: Iranian Heavy: 50% by weight - Zakum: 35% by weight - Djeno: 15% by weight.

The last column in the above table indicates the absorbances obtained by analysis of the mixture.

By a least square method the nine equations are solved in relation to the three unknowns (% by weight of each of the crude oils in the mixture).

The values thus calculated are as follows, in % by weight: Iranian Heavy: 51.4% - Zakum: 34.2% - Djeno: 14.4%

The accuracy, in relation to the control mixture, is therefore satisfactory, the differences observed being acceptable and due to the repeatability of the absorbance measurements.

This method enables the composition of a mixture of known crude oils to be obtained in real time.

It is therefore possible, if the distillation unit is supplied from tanks containing different crude oils, to know, in real time, the composition of the mixture of crude oils which is supplied to the distillation unit.

This method may be applied to measuring the supplies of mixed crude oils.

EXAMPLE 17

Distillation

A Djeno crude oil was used and the absorbances D.sub.i were measured for 12 frequencies F.sub.i. The results are summarised below:

______________________________________F.sub.2 4640 cm.sup.-1 D.sub.2 = 0.01364F.sub.3 4615 D.sub.3 = 0.01758F.sub.4 4584 D.sub.4 = 0.01895F.sub.5 4485 D.sub.5 = 0.03172F.sub.6 4385 D.sub.6 = 0.45629F.sub.7 4332 D.sub.7 = 0.94434F.sub.8 4305 D.sub.8 = 0.63842F.sub.9 4260 D.sub.9 = 0.80701F.sub.11 4170 D.sub.11 = 0.57862F.sub.13 4100 D.sub.13 = 0.51222F.sub.14 4060 D.sub.14 = 0.52624F.sub.15 4040 D.sub.15 = 0.41788______________________________________

The yields of the defined fractions previously obtained by multivariate regression are expressed as a function of the absorbances D.sub.i of the different frequencies by the following expressions:

Residue output=359.88-1385.51 D.sub.2 -441.15 D.sub.6 -291.84 D.sub.7 +225.1 D.sub.8 +128.33 D.sub.15 Light distillate output=-113.6+541.26 D.sub.2 -108.75 D.sub.6 +122.71 D.sub.7 -55.49 D.sub.8 -31.64 D.sub.15 Gas oil output=-54.03+373.77 D.sub.2 +199.13 D.sub.6 +67.77 D.sub.7 -101.61 D.sub.8 -57.92 D.sub.15

The calculation carried out on the basis of the above formulae and of the absorbances gives the following results, compared with the experimental values observed for the Djeno crude oil examined here:

______________________________________ Calculated Observed______________________________________Yield of residue, % by weight 61.4 63Light distillate 10.6 9.1Gas oil 16.8 16.2______________________________________

If this type of regression method is extended to the entire distillation range, it is then possible to calculate the distillation curve and the mean distillation temperature of the feed in real time.

This method is also applicable to mixtures and feedstocks of variable compositions.

In a computer-assisted process in a distillation unit, the NIR spectrum of the feedstock entering the unit is thus taken in real time and treated as an information vector which characterises the distillation properties of the feedstock in the distillation process on a continuous basis. The amount of detail in the near infra red spectrum and the experimental accuracy which may be obtained as required from the accumulation of spectral data, using fast Fourier transformation, ensures that this information is reliable and relevant to distillation.

In the following examples it is demonstrated that for a fixed set of conditions in the unit, the variations of the quantity of the products formed may be correlated, by means of a numerical treatment, with variations in the NIR spectrum of the feedstock.

EXAMPLE 18

Distillation

The Conradson carbon residue according to Standard T60-116) of a crude petroleum of Mexican ISTHMUS origin was 4.1% by weight.

The absorbances of a test sample of this crude petroleum were measured for 3 frequencies, and had the following values:

______________________________________F.sub.2 4640 cm.sup.-1 D.sub.2 = 0.0186F.sub.7 4332 D.sub.7 = 0.9435F.sub.8 4305 D.sub.8 = 0.6374______________________________________

The value for the carbon residue was calculated from the following equation:

C=84.306-64.1794 D.sub.7 -252.437 D.sub.2 -22.432 D.sub.8 /D.sub.7

In the case of the crude petroleum, the figure of 3.9% by weight was obtained, which was satisfactory, taking into account the accuracy of the test.

EXAMPLE 19

Distillation

The asphaltene content of a crude petroleum of Suez Blend origin, measured by the AFNOR T60 - 115 method was 1.5% by weight.

The NIR spectrum of this crude was determined at the following frequencies and gave the following absorbances:

______________________________________F.sub.2 4640 cm.sup.-1 D.sub.2 = 0.01868F.sub.7 4332 D.sub.7 = 0.93972F.sub.8 4305 D.sub.8 = 0.6358F.sub.2 4640 cm.sup.-1 D.sub.2 = 0.0210F.sub.7 4332 D.sub.7 = 0.9362F.sub.8 4305 D.sub.8 = 0.6427______________________________________

The asphaltene content was obtained by the following equation correlating the absorbances:

Asphaltenes (% by weight)=22.31 (3.758-2.861 D.sub.7 -11.253 D.sub.2 -D.sub.8 /D.sub.7).sup.1.5 -0.465

The calculated value was 1.3% by weight which was satisfactory.

EXAMPLE 20

Distillation

A gas oil fraction was used which was taken from a Iranian Heavy crude petroleum. The cetane number of this fraction was 49.6 (Standard MD7 - 035).

The NIR spectrum was obtained for the crude petroleum. Absorbances D.sub.i were measured at the following frequencies:

The absorbance values were used in the following equation in order to calculate the cetane number of the gas oil from the NIR spectrum of the crude:

Cetane number=155.529-12.053 D.sub.7 -47.408 D.sub.2 -135.217 D.sub.8 /D.sub.7

The value obtained was 50.4.

Thus, the above equation enables the cetane number of a gas oil to be determined as a function of the NIR spectrum of the feedstock.

EXAMPLE 21

Distillation

The aniline point of a gas oil fraction from a light crude petroleum of Lokele, Cameroons origin, was determined (Standard ASTM-D611). This was 52.degree. C.

The values for the absorbances measured at the following 4 frequencies of the NIR spectrum of this crude petroleum were:

______________________________________F.sub.6 4385 cm.sup.-1 D.sub.6 = 0.5181F.sub.7 4332 D.sub.7 = 0.8268F.sub.8 4305 D.sub.8 = 0.6786F.sub.15 4040 D.sub.15 = 0.4293______________________________________

The following equation enabled the aniline point to be calculated from the absorbances of the crude:

Aniline point of the gas oil produced=27.1+53.57 D.sub.7 +113.8 D.sub.15 +199.6 D.sub.8 -391.2 D.sub.6

The value thus calculated was 53.degree. C.

EXAMPLE 22

Distillation

The cloud point of a gas oil fraction taken from an Egyptian Belayim Crude was -10.degree. C. (Standard T60 - 105).

The values for the absorbances measured at the following four frequencies of the NIR spectrum of a sample of this crude petroleum were:

______________________________________F.sub.5 4485 cm.sup.-1 D.sub.5 = 0.0335F.sub.6 4385 D.sub.6 = 0.45821F.sub.13 4100 D.sub.13 = 0.49486F.sub.15 4040 D.sub.15 = 0.41039______________________________________

EXAMPLE 23

The viscosity at 60.degree. C. of the 369.degree.-509.degree. C. fraction of a Congolese Djeno crude petroleum was 27.2 cSt.

The values for the absorbances measured at four frequencies of the NIR spectrum of a sample of this crude petroleum were:

______________________________________F.sub.2 4640 cm.sup.-1 D.sub.2 = 0.01364F.sub.7 4332 D.sub.7 = 0.94434F.sub.8 4305 D.sub.8 = 0.63842F.sub.14 4060 D.sub.14 = 0.52624______________________________________

The value calculated by introducing the measured absorbances in the following equation was:

log.sub.10 viscosity at 60.degree. C. =1.4468-5.662 D.sub.7 +0.874 D.sub.14 =29.980 D.sub.2 +8.310 D.sub.8 =28.5 cSt

The equation enabled the viscosity of the residue to be calculated from the NIR spectrum of the feedstock.

Claims

1. The method for determining at least one of a property of a product and a yield of a product from a process selected from the group consisting of a hydrocarbon conversion process and a hydrocarbon separation process, comprising the following steps:

(a) determining with an infrared spectrometer an absorbance of a feedstock to the process at a certain number of frequencies in the spectral range 16667 to 3840 cm.sup.-1, starting from a defined base line, and

(b) determining at least one of (i) the property by applying a correlation between the property of the product and the absorbance values of the feedstock and (ii) the yield by applying a correlation between the yield of the product and the absorbance values of the feedstock wherein the correlation is determined experimentally by multivariate regression and is dependent upon the spectrometer used, the property or yield to be determined, and the frequencies used.

2. The method according to claim 1 wherein the frequencies are in the spectral range 12500 to 3840 cm.sup.-1.

3. The method according to claim 2 wherein the frequencies are in the spectral range 4760 to 4000 cm.sup.-1.

4. The method according to claim 1 wherein the frequencies used are selected from the following:

4670 cm.sup.-1

4640 cm.sup.-1

4615 cm.sup.-1

4585 cm.sup.-1

4485 cm.sup.-1

4405 cm.sup.-1

4390 cm.sup.-1

4385 cm.sup.-1

4332 cm.sup.-1

4305 cm.sup.-1

4260 cm.sup.-1

4210 cm.sup.-1

4170 cm.sup.-1

4135 cm.sup.-1

4100 cm.sup.-1

4070 cm.sup.-1

4060 cm.sup.-1, and

4040 cm.sup.-1.

5. The method according to claim 1 wherein the base line is taken as 4780 cm.sup.-1.

6. The method according to claim 1 wherein the spectrometer is linked to a signal processing device to permit numerical treatment of the spectral range.

7. The method according to claim 6 wherein the treatment is by Fourier transformation.

8. The method according to claim 1 wherein the method is on-line and in real time.

9. The method according to claim 1 wherein the correlation contains linear terms, quadratic terms and homographic terms.

10. The method according to claim 1 wherein the process is a hydrocarbon conversion process which is a catalytic reforming process and the properties are selected from the following: yields of gasoline (C.sub.5.sup.+), hydrogen and gas, octane numbers (clear or leaded), vapor pressure and density.

11. The method according to claim 1 wherein the process is a hydrocarbon conversion process which is a catalytic reforming process and wherein the method comprises the further steps of:

(c) carrying out differential spectrum measurements on the feedstock and the product in a range and at frequencies which are the same for the feedstock and the product,

(d) applying a correlation between the differential spectrum and an octane number gain determined for the product, the correlation being determined experiment by multivariate regression analysis,

(e) determining an octane number and octane number gain for the product and comparing a sum of the octane number and octane number gain with a desired objective, and

(f) modifying at least one operating condition of the catalytic reforming process until the determined octane number gain is in agreement with a desired octane number gain.

12. The method according to claim 10 wherein the frequencies used are selected from the following:

4670 cm.sup.-1

4485 cm.sup.-1

4332 cm.sup.-1

4305 cm.sup.-1

4210 cm.sup.-1

4100 cm.sup.-1

4060 cm.sup.-1, and

4040 cm.sup.-1.

13. The method according to claim 1 wherein the process is a hydrocarbon conversion process which is a selective hydrogenation process and the property is a maleic anhydride number of the product.

14. The method according to claim 1 wherein the process is a hydrocarbon conversion process which is a catalytic hydrocracking process and the properties are selected from the following: total conversion of the feedstock to product and conversion of the feedstock to gasoline.

15. The method according to claim 14 wherein the frequencies used are selected from the following:

4585 cm.sup.-1

4390 cm.sup.-1

4332 cm.sup.-1

4260 cm.sup.-1

4210cm.sup.-1

4170 cm.sup.-1

4100 cm.sup.-1, and

4070 cm.sup.-1.

16. The method according to claim 1 wherein the process is a hydrocarbon conversion process which is a catalytic cracking process and the properties are selected from the following:

(i) conversion of the feedstock to gasoline and yield of residue in the product, and

(ii) density, research octane number (RON) and motor octane number (MON) of clear and leaded gasoline from the catalytic cracking process.

17. The method according to claim 16 wherein the frequencies used are selected from the following:

4670 cm.sup.-1

4615 cm.sup.-1

4485 cm.sup.-1

4135 cm.sup.-1

4100 cm.sup.-1, and

4060 cm.sup.-1.

18. The method according to claim 1 wherein the process is a hydrocarbon conversion process which is a visbreaking process and the properties are selected from the following: yield of a gas oil product, viscosity, specific gravity and xylene equivalent of residue, and conversion limit.

19. The method according to claim 1 wherein the process is a hydrocarbon separation process which is a distillation process and the properties are selected from the following:

(i) a quantitative composition of a mixture of crude oils constituting the feedstock when the product is known to be a blend of crude oils,

(ii) yields of products obtained,

(iii) densities and viscosities of the products,

(iv) cetane number, aniline point, cloud point and KUOP characterisation factor of a gas oil product, and

(v) n-paraffinic, isoparaffinic, cycloparaffinic and aromatic content of a gasoline product and a research octane number (RON) and motor octane number (MON) of clear and leaded gasoline.

20. The method according to claim 1 wherein a value of at least one of a property or yield of a product is determined by an NIR spectrum of the feedstock, and wherein the process is computer controlled by a feedback control system which will alter process conditions in response to variations in that value from a desired value.

21. An apparatus for carrying out a method for determining at least one of a property of a product and a yield of a product from a process selected from the group consisting of a hydrocarbon conversion process and a hydrocarbon separation process, wherein the method comprises the following steps:

(a) determining with an infrared spectrometer an absorbance of a feedstock to the process at a certain number of frequencies in the spectral range 16667 to 3840 cm.sup.-1, starting from a defined base line, and

(b) determining at least one of (i) the property by applying a correlation between the property of the product and the absorbance values of the feedstock, and (ii) the yield by applying a correlation between the yield of the product and the absorbance values of the feedstock, wherein the correlation is determined experimentally by multivariate regression and dependent upon the spectrometer used, the property or yield to be determined, and the frequencies used;

said apparatus comprises an infrared spectrometer and a computer wherein the computer is programmed to receive and utilize the absorbance from the infrared spectrometer to determine the at least one of the property and yield continuously and in real time.

22. The apparatus according to claim 21 wherein a value of at least one of a property or yield of a product is determined by an NIR spectrum of the feedstock, and wherein the process is computer controlled by a feedback control system which will alter process conditions in response to variations in that value from a desired value.