CRACKING PROPERTY DETERMINATION

Abstract: A method of determining or predicting a value Px of a property of a feed X to a steam cracking process or yield of said process, which method comprises measuring the absorption Dix of said material at more than one wavelength in the region 600-2600 nm, comparing signals indicative of said absorptions or a mathematical function thereof with signals indicative of absorptions Dim or mathematical function thereof at the same wavelength for a number of standards S in a bank for which the said property or yield P is known, and choosing from the bank at least one standard Sm with property Pm said standard having the smallest average value of the absolute difference at each wavelength i between the signal for the material and the signal for the standard Sm to obtain Px, with averaging of said properties or yields Pm when more than one standard Sm is chosen.

Description: Case 8499( 2 ) CRACKING PROPERTY DETERMINATION This invention relates to a method of determining or predicting by near infra red (NIR) spectroscopy properties of feeds or products and/or yields in steam cracking processes in particular involving hydrocarbons, especially in hydrocarbon refineries NIR spectroscopy has many advantages over other methods of analysis in refineries and can cover a large number of repetitive applications accurately, quickly and on line. The NIR region between 800 and 2500nm contains the totality of molecular information in the form of combinations and overtones from polyatomic vibrations, but Mathematical techniques are needed to exploit this information and to calculate the desired parameters. EP-A-285251, 304232, 305090, and 345182 the disclosure of which is hereby incorporated by reference, describe respectively the use of NIR for determining octane number, yields and/or properties of a product of a chemical process or separation process from analysis on the feeds to that process, and yields and/or properties of a product of a blending operation again from analysis on the feed thereto and for determining at least one condition of operation of the steam cracking reaction. At present, numerical methods described for modelling physicochemical properties based on NIR spectra all are of a correlative nature and involve relations of a regressional character between the property(ies) studied. Among these multivariable analyses are multilinear regression (MLR), Principle Component Regression (PCR), Canonic regression, and regression by Partial Least Squares (PLS). In all cases there is sought between the property and the NIR spectrum a relation which may be linear but is usually quadratic or of higher algebraic form involving regression coefficients applied to each absorption. The establishment of any regression requires a progressive calibration, as the approach is empirical and not supported by a theory. These techniques have disadvantages, the chief of which is the need for establishing a strong correlation between the spectrum and the property, and their difficulty in dealing with positive or negative synergy between components contributing to that property. Furthermore there are practical difficulties mainly in the need to identify sample families having the same kind of relation between the spectra and the properties to be modelled. Thus the model may be limited especially with a non linear relation between spectrum and property. Especially when at the edges of the available data the accuracy of the model diminishes. The stability of the model is also a problem, as is the need when adding new standards to do laborious revisions to give the new model, especially when adjusting to a new feedstock for a process; thus testing the yields on 4 products leaving a steam cracker requires 4 models, each of which has to be changed for each change of the feed not included in the calibration. We have discovered a new approach avoiding the above problems with correlations, and regression calculations, and being capable of being expanded automatically with use of a new feed of different quality. The present invention provides a method of determining or predicting a value Px, of a property of a feed to a steam cracking process or yield of said process which method comprises measuring the absorption DJX of said feed at more than one wavelength in the region 600-2600nm, comparing signals indicative of said absorptions or mathematical functions thereof with signals indicative of absorptions Djm at the same wavelengths or mathematical functions thereof for a number of standards S in a bank for which the said property or yield P is known, and choosing from the bank at least one, and preferably at least 2 standard(s) S m with properties or yields P m said standard S m having the smallest average values of the absolute difference at each wavelength i between the signal for the feed and the signal for the standard S m to obtain value P x and with averaging of said properties or yields Pm when more than 1 standard S m is chosen The above method can be performed without regression or correlation techniques, e.g. between the absorptions at any wavelength of the material and the property/yield. Thus for the performance of the method of the invention, a bank is prepared in which the NIR spectra are recorded at many wavelengths for a large number of standard materials, together with their properties determined by alternative techniques e.g. gaschromatography for chemical compositions or yields determined by known methods. The standards are chosen to cover the area in which the method is to be used, so for determination of the yield of ethylene, a range of feeds to the steam cracker is chosen and under constant cracking conditions eg temperature pressure and catalyst the yield of ethylene is measured The number of wavelengths chosen may be 2-1000 e g 5-200 or 20-90 such as 60-80 while the number of standards can be at least 100 or 1000, or 100,000 up to 5 million depending on property(ιes) chosen The wavelengths chosen may be at regular intervals such as each 1-50 or 15- 35nm, (or each 1-5 nm or each nanometer) or may be at n regular intervals e g with intervals of 1- l OOnm, e g 1-50 especially l-35nm which may be random or chosen because of a change in the shape of the spectral curve at that wavelength e g a peak, trough or shoulder or chosen by chemical or statistical criteria such as factor analysis The wavelengths may be in the region 600-2600nm, eg 1000-2500nm but preferably 1500-2600 such as 2000-2550nm while the wavenumbers may be 16,600-3840cm- 1 , e g 10000-4000CIT.- 1 , 6660-3840, or 5000-3900 cm" 1 Frequencies in Hertz are obtained by multiplying the wavenumbers bv 3 x 10 '0 cm/sec The signals eg absorptions (oi derivatives) tor the unknown sample are compared with the signals eg absoi ptions (or derivatives) at the same wavelength of the standards, and those standaids chosen having the smallest differences The properties of those chosen standaids aie then averaged to determine the property or yield of the unknown sample The absorptions at more than one wavelength may be chosen, e g 2- 1000 such as 5- 100 or 10-20 In the method of the invention the standards chosen are those with the smallest average values of the absolute difference at each wavelength i between the signal exemplified by the absorption/optical density (or a derivative theieof) D 1X for the unknown material and the corresponding signal/function eg absorption/optical density (or derivative thereof) D ιm for the standard The averages may be in respect of the mean value of D ιx -D ιm (whatevei its sign le the absolute differences) or (D ιx -D ιm )2 and may be the simple mean value oi the differences may be weighted to take account of the different sensitivity of the absorption to the property at that wavelength or the different sensitivity of the spectrometer at that wavelength For each standard in the bank of standards for the type of mateual in question, the average difference is found as described and the standards with the smallest average differences chosen, e g at least 1 but preferably at least 2 such as up to 1000 smallest such as 1 (or 2) to 100 or 1 (or 2)-20 but in particular 1 (or 2) -10 and especially 2-6 smallest Advantageously the average differences chosen and hence the standard (or standards S m ) chosen for the property or yield wanted are such that in relation to the unknown material X and each chosen standard S m the following function is met ι xm /ΣD lx <expeπmental erroi wherein ι xm is the proximity index and is defined by ι^( x ) = ∑(D ] χ-D ιrn )2 and the experimental error is in determining said property or yield in the standard The value P x of the property or yield is the same as P m or the average of Pm when more than one standard S m is chosen In order to aid the choice of the appropriate standards, especially in relation to a large number of wavelengths for a complex unknown mixture, it is preferred to limit the choice to those defined by means of a minimal index For the chosen standard the minimal index is at least the same as the differences between the absorptions of the unknown and the standards Mathematically, this may be expressed as where iM is the minimal index for the property, and lab is a measure of the deviation (called the proximity index) at all the chosen wavelengths between absorption of the unknown and a chosen standard b That measure is defined by ι(ab)2 = ∑, (D ιa -D lb )2 (I) where D ιa is the optical density (or absorbence) of unknown a at wavelength i (or a derivative theieof e g a first, second oi third derivative of that density), and D,^ is the optical density (or absorbence) of standard b at that wavelength i (or a derivative thereof e g a first, second or thud derivative of that density) The value of Dj is the optical density or the optical density difference with tespect to the baseline of the spectrum at that wavelength, oi the baseline interpolated between 2 wavelengths on either side thereof If desired signals conesponding to mathematical function of the absorption eg after Fourier transformation oi derivatives or spectral subtraction or division may be used to provide corresponding proximity and Minimal Indices If desired instead of the optical density D, a normalized density Wj may be used where Wj =D,/∑D, This normalization avoids errors due to small electronic fluctuations in the apparatus and compensates for small differences in the optical path between the optical cells In this case the proximity index is defined by I(ab)2 = ∑, (W ιa -W lb )2 (2) The indices can be weighted as desned for increased resolution One approach is to define the indices as follows I(ab)m = ∑ Abs val ( X, a -X,b) σ ι n (3) where Xj is D, or W, or a mathematical combination thereof, σ, is the standard deviation of X for the set of samples considei ed (at that wavelength) and each of m and n, which are the same different is the weighting factor which is positive but can be a whole numbei or a fraction Other vai lants can be used with other weighting factors such as those involving the spectral experimental error e,, where e, is the reproducibihty of the spectral measurement at wavelength i The choice between the different options for the weighted indices may be dictated by numerical efficiency The reproducibihty of the expenmental measurements in the standards may be at least 90% oi 94% or 95% The minimal index may be obtained fi om a reference standard samples set according to the following piocedure, hereafter called the Minimal Index Procedure The NIR spectra for the standard samples A and B and their property P e g density are determined By means of equation ( 1 ), (2) or (3), the value of the proximity index ι ab is determined via the absorptions at a series of wavelengths, this index is applicable to the difference in properties P 'Pb called EP ab This process is repeated with other pairs of standards c and d, e and f ete to obtain a series of Proximity Indices ι cc j etc with corresponding property differences EP cc j etc For different values of a parameter L which is greater than the indices ι ab etc, the corresponding values of EP ab etc are averaged to give an average EP,, for that value of L, the different values of EP,, + tσ/Vκ are then plotted on a graph against L σ is the accuracy of the property deteimination and K is the number of pairs of samples for which l is inferior to a given L t is the Student factoi at a given level of confidence The intercept is then measured between the curve obtained and a line usually horizontal which is the reproducibihty of the propertv level at an appropriate confidence interval e g 90% or more usually 95% The abcissa portion of the intei cept gives the minimal index ι mm which is the minimum value ot ι ab for which P a =P b within the frame of experimental en oi From this minimal index by Pi ceduie I, the standards can be chosen which have values of |2 a ^ < |2 mιn wheie in this case a is the unknown and b is a standard, as in this case the difference between Pi operty a and Property b is less than or equal to σV 2, where σ is the experimental error in measuring the property Then from the property P values of the chosen standards, the property of the unknown is obtained by averaging those values, usually the arithmetic mean, but optionally with weighting The method of the invention may be used to determine more than one Property P at once, e g at least 2, such as 1-30 e g 2-10 properties at once Each property of the standards has a particular minimal index, which may lie in the region 0-10-^ eg 10- to 10" 8 , 10" 5 to 10' 8 in particulai l O' 8 to 5xl 0" 7 or 5x l 0" 7 to 5xl 0" 4 especially for Minimal Indices derived fiom absoi bances coπesponding Minimal Indices may be obtained for other signals/functions If the Minimal Index chosen is the smallest for all the properties desired then the same one may be used foi all the properties and the standards chosen will be suitable for all the properties The Minimal Index for each property may be used separately, with different numbers of standards chosen for each property (assuming different Minimal Indices) If desired the same Minimal Index may be used, which is not the smallest, resulting in some of the chosen standards (with less than their Minimal Index) giving some properties of high accuracy and some (with a higher Minimal Index) giving properties of less high accuracy. The property to be determined may be of the sample being analyzed or the yield of a product obtained from that sample i e a product of cracking the sample, as the property value obtained is derived from the standards, and they will have been determined as needed for the eventual use EP3-34232, 305090 and 345 182 referred to above describe such techniques when applied to use of NIR with correlation to blending or cracking operation, the same principles apply in the present method If the density of the standards in the bank is sufficient to have -- a ^ ≤ i m j n as is usually the case, the above procedure is very satisfactory But there are occasions when the bank is incomplete, because of shortage of data of properties in a particular area i.e. a low density of standards or the sensitivity of the property to changes in absorption is so small, that a very small Minimal Index is required and there may be few standards with proximity indices meeting it It is possible simply to use a larger Minimal Index with e.g 1-5 times such as 1 5-2 times the Minimal Index; the results may be less accurate than those from a smaller Minimal Index However, a more accurate approach with a low density of standards involves a special densification process of Procedure 2, in which random or semi random densification of the neighbourhood of the unknown is achieved by generation of synthetic standards, based on standards already in the bank Each new synthetic standards may be obtained from combinations of standards taken at random from the bank but preferably it is obtained from the other standaids by the constraint of choosing only a mixture of N standards for which (Min)Cj - iij <Cjj <(Max)Cj+uj (4) and ∑Cjj = I (5) where C;; is the fraction of component j in the sample j Min C; is the minimum amount of ; in the initial calibration, i.e standards in the bank, which covers the correct area of ; in the industrial mixture, or in the samples for which the method is to be used Max Cj is the maximum amount of ; in the initial calibration, i.e standards in the bank, which covers the correct area of ; in the industrial mixture or in the samples for which the method is to be used uj is usually between 1 and 0 01 (or 1 and 0 05) preferably between 0 5 and 0 1 and can be fixed for each property The constraints over the choice of such mixtures of N standards can also be equally fixed in the spectral area from which the samples will be drawn in order to remain in the areas of similar chemical nature The number of samples effectively drawn into the bank in this densification can be of several thousand generally 1000-2000 The calculation time is extended without significant deterioration in the results If no further neighbours are found, the trawl of new samples drawn in is enlarged The spectrum of each mixture is calculated by the combination of the spectra of the standards used according to the formula S Mi = Σ CijXS, (6) where S; is the spectrum in the mixture of component ; in the calibration matrix The properties of each mixture PMi can be calculated by a generally linear combination of the properties of the standards according to the formula where P; is the property of component | In the case of non linear additive properties, appropriate mixing factors can be applied e g by blending factors or similar for density and viscosity Having obtained the spectrum and the properties of the synthetic mixtures, these can be used as "standards" to help determine the properties of an unknown sample in the same way as a conventional standard Instead of using either of the two above approaches, 1-7, a third type Procedure 3 may be used as follows The Q nearest samples to unknown X can be found from a selection from the bank samples for which the proximity index to the unknown sample is (V) X i mjn ) where v is 0 I < v < 10, (8) preferably 0 5< v <2 or 1< v <5 Then by the method of least squares is found a generally linear combination of the standard products which are (the Q nearest samples to reproduce the spectrum of X according to the equation S x = ∑ CR X S,- (9) where C r is the coefficient for sample R in the total Q and SR is the spectrum of sample R The coefficient CR which can be normalized to CR = I or not and/or optimized by the least squares route, allows an estimation of the property P x according to the equation P x = ∑ CR X PR ( 10) where PR is the property of sample R The eventual size of the estimation error can be derived by application of Gaussian theory, also called the propagation error (see Eq 10) The above third approach can only be applied if the product X is situated inside the maximum extension of the standard products defined by equation (8) I e within the range of bank samples defined in equation (8) If this is not the case, X is outside the field of the actual bank of products and escapes from the area of knowledge of the method into the area of learning The densification process described in relation to equations 4-7, or 9 or 10 is usually applied to the method of the invention involving no correlation or regression techniques Howevei, if desired the densification process may be applied to increase the number of "standards" for consideration in an NIR analytical technique involving the correlation or regression techniques as described above e g MLR The present invention also provides a method for adding an extra synthetic standard to a bank of known standards, each of which relates at least one absorption in the 600-2600nm region (or a signal indicative thereof oi of a mathematical function of said absorption eg a derivative thereof) of a feed to a steam ciacking process to a known property related to that feed or yield of said pi ocess, which method compnses choosing from the bank at least 2 standards foi which equations 4 and 5 above are met, considering mixing the chosen standards in at least one proportion to produce at least one mixture for use as a synthetic standard, and estimating the spectrum and property/yield of said mixture according to equation 6 and 7 l espectively The spectrum and property/yield of each "mixture" can then be added to the bank and may be used to develop models through the known correlation/regression approach eg as described in EP345182 The method of this invention is applicable to determination of the properties of feeds to a steam cracking, as well as the determination from the feed to the steam cracker of the potential yields of C2- ^ alkenes e g ethylene, propylene, butenes and a C5 olefin fraction, or benzene in the steam cracking under standard conditions Other "properties" of the feed which relate to its suitability for use in steam cracking under standard conditions may be measuied such as its potential yield of methane and its cokability, le the tendency to form coke in that ciacking Othei properties of the feed that can be determined in addition include its chemical composition eg percentages of linear paraffinic hydrocai bons. isoparaffins, naphthenics compounds, aromatics and benzene, and its density and mean molecular weight These parameters are important for the smooth operation of a steam ciackei as variations affect that operation, especially the olefin production Such variations may result from changes intentional or otherwise in the feeds, which may be naphtha, condensates, liquefied gas and/or gas oil The steam cracking may involve amounts of steam of 30-50% eg 40% (relative to hydrocarbon feed) riser inlet tempei atuies of I00-400°C eg 250°C and riser outlet temperature of 800-900°C eg 830°C and l esidence time of 0 1-5 sec eg 0 5 sec In each case the method may be applied to determine a property of the feed tested or yield from that feed prefeiably on line and especially with contmous feed back from the results to control the production process In each of the above piocesses the property of a teed or yield determined or predicted by the method of the invention can be compaied to the desired figure and notice taken ot any deviations by ad|ustιng the parameters of the cracking process e g proportion or nature of feed le steam oi hydrocaibon (e g via operation of control valves) and/or temperature/pressure etc to bring the property back to the desired figure This control of the cracking pi ocess, is usually performed with a micro computer which is linked to the spectrometei and also pei forms the search for the standards S m The inline control ot the process is very efficient and very fast The present invention also piovides an apparatus suitable for carrying out the method of the invention compi ising an infia red spectrometer and a computer wherein the infra red spectrometei is linked to the computer program in such manner that the property may be determined continuously and in real time The spectrometer is suitable for measuring spectia at least partly in the 600-2600nm wavelength lange and can be linked to a signal processing device to allow numerical treatment of the spectrum, preferably by Fouπei transtoi niation The spectiometer receives at least one signal from a vessel containing feed oi ti om a feed line The information obtained can be used as an information vectoi foi the computei which is programmed to determine the property oi yield eg via calculations on the proximity indices in relation to standards Conveniently in relation to a process, the computer may be used in a closed loop feed-back or -forward control system for controlling processing equipment eg changing the process parameters in iesponse to variations in the property of feed and/or yield of product from the desned value, from measurement of more than one absorption in the NIR spectrum of the teed The present invention also pi ovides a computer programmed to perform the method of the invention of detei mining oi predicting the value P x of the property or yield, or programmed to perform the method ot the invention of adding an extra synthetic standard to the bank of known standards The apparatus for use with the former method of the invention compπses an NIR spectrometer receiving at least one signal from a feed line in said process and being coupled to a computer to effect continuous measurement of the spectra of the feed and provide feed-back or -forward control of the process The present invention also provides a computer implemented method for a system including a spectrometer linked to a process line containing a material X, which is a feed to a steam cracking process, a computer linked to the spectrometer, and a controller linked to the computer and the process line, the computer including databanks having stoi ed therein signals indicative of absorptions of standard materials or mathematical functions thereof and corresponding properties of said materials, or products of said piocess of which X is a feed or yield of said process, the method comprises steps ot measui ing absoi ption at moie than one wavelength in the region 600-2600nm at the pi ocess line and producing absorption signals (or mathematical functions eg derivatives thereof) by the spectrometer in accordance therewith, accessing the databanks of the computei in accordance with the absorption signals (oi functions thereof), compaπng, by the computer, the signals (or functions thereof) to the signals (or functions theieot) of the standaid materials stored in the databanks, choosing at least one standard based on the comparing, said standard having the smallest average value of the absolute difference at each wavelength I between the signal (or functions thei eot) foi the material and the signal (or function thereof) for the standards, with avei aging of said pioperties or yields when more than one standard is chosen, and controlling said process in accordance with the outputted property/yield, preferably each signal is an absorption value or derivative thereof The benefits of invention allow improvements in modelling with the following areas, identification and classification of novel feeds, simultaneous estimation of all of P properties on a sample without the need foi generating P different models, and with the option of automatic upgrading of the model, the method being self learning or adjusting The method of the invention overcomes the difficulties with the classical regressional approach, in particulai avoiding all difficulties with numerical stability of the models, allowing easy and rapid identification and classification of a sample of a feed analyzed by spectral recognition and then instant conclusions as to whether the sample is known oi unknown, allowing simultaneous deteimination of many properties The method also allows an extension of the field of application of the method without the need to rewrite the model, apart fiom the need to integrate the new samples which are inside or outside the pi evious field of validity of the method This possibility of automatic learning, which is not possessed by traditional regression techniques, is a decisive advantage in the fi amewoi k of continuous inline industrial control processes, because it allows the leturn of the industrial plant operations to the model in a certain and rapid manner in a minimum time and with all the properties considered in the model In contrast classical regression methods would necessitate the redevelopment of all the models, which is long and laborious without being able to guarantee the result of the new model obtained, because a new validation period is necessary, in addition during the redevelopment of the model any commercial refinery use of the model is very limited Furtheimore, the method of invention allows equally the easy extension to a number of properties, which are simply incorporated into the known bank This remarkable possibility is ti ne not only for conventional properties such as physical chemical and/or rheological pi operties, but also for complex ones The methods of the invention equally allow application of the models from one apparatus to another and from one spectral region to anothei , where conventional regressive method cannot give satisfactoiy solutions This apparatus portability is made possible by the fact that the distinctions between different spectra ai e the same in one apparatus as another, for the same type of specti onietei being considered (e g network scatter, Fourier transform, accousto optical system AOTS, diode etc) This portability between spectral region depends on the tact that as the spectial regions are intercorrelated, the relations between the spectra are maintained between one another The invention is illustrated in the accompanying Figures in which Figure 1 represents a schematic diagi am showing apparatus for use in the invention. Figure 2 represents a schematic block flow diagram for the method of the invention In Figure 1, an optical fibie 3 links a spectrometer 2 and a probe 6 in or at process line 1 The spectrophotometei 2 pioduces absorbance signals at more than 1 wavelength, which signals as such (oi after mathematical treatment e g to form derivative signals) are passed via line 4 to computer 5, where the signals as such or after conversion e g to one oi more deuvative signals, are used to enable the computer to access the databank 7 of standard absorptions and properties/yields therein The signals are compared to those of one oi more standaid absorptιon(s) as described above and its/their corresponding propeιty(ιes) or yιeld(s) The output of the computer 5 is in the form of a property of the material in line 1 or yield of the process from that material and may be printed in hard copy Preferably however, the output property/yield is used to control the piocess involved with the material in line 1, le for which line 1 is a feed, in this case the computei 5 is linked to and instructs the controller 9 which, via 10, controls that piocess eg via valves/temperature and/or pressure contiols in line 1 or in i elation to line 1 By this means the property of the material in line 1 or yield or property ot product of the process therefrom can be optimized can be optimised In Figure 2, the initial opeiation l l is to measure the absorption of the unknown, aftei which in the second step 12 the absorptions are compared to absorptions in spectra of standards, and in the third step 13, the spectra of the standards Sm are chosen according to cuteria described above, and then in step 14, the property(ιes) of the standard(s) Sm chosen is used to obtain the desired property or yield If the spectrum of only l standaid Sm is chosen, then the value P x of the unknown is the same as that of that standard Pm If more than 1 spectrum is chosen, the value P x of the unknown is the average of the values Pm of the standards If desired in an optional step 16, the value P x is compared to the desired value for the unknown and in step 16 the process involving the unknown is adjusted to make the value P x the same as the desn ed value The invention is illusti ated in the following Examples in which the Minimal Index is calculated accoi ding to the Minimal Index Procedure descnbed above Mathematically the steps concerned ai e as follows For each couple of standai d samples i, j, the Proximity Index i,, is determined from the NIR spectra by use of equation 1 , 2, or 3 and the pi operties aie measured For each Proximity Index is calculated the absolute difference EP„ between the properties of the samples The Minimal Index for property P is obtained from the average EM p (L) of EPJ; foi different values of L when L ≥ij Thus the EM p (L) = (1/K) ∑, ∑, EP,; for each of K samples for which ι,, <L EMp(L) + tσ iK is plotted against the pioximity index and in addition there is plotted the reproducibihty of the standaid method at a given confidence level as defined in the Minimal Index Procedure above The intercept of the curve from EMpL and the reproducibihty give the uppei limit I e the Minimal Index For the Examples the data is expressed in Tables in a form as shown below in which the data is as follows The wavelengths chosen ai e shown in columns 1 and 2 Column 3 gives the weight loading associated with each wavelength for the proximity index for the standai ds 1 denotes no loading Column 4 shows for the unknown sample the absoi ption at the various wavelengths and at the bottom the pioperties ot that sample detei mined by standard methods Column 5 shows for the unknown sample the estimated values of the properties and the absorptions using the method of the invention based on the properties and absorptions of the chosen standai ds Columns 6, 7 etc show the values of the absoi ptions and properties for the standards chosen from the bank Line 2 give the value of the pioximity index between the unknown sample and each of the chosen standards Example 1 Determination of the pioperties of a mixtuie of hydiocaibons for feeding to a steam crackei and the yields of products theieti om The pi operties were sought ot a mixtui e of hydrocaibons, a naphtha, a petroleum hydrocaibon cut which was a feed to a steam thermal cracker for producing olefins, such as ethylene, propylene and butenes, as well as ci acked gasoline The key properties of the feed are its density and its "potential content of ethylene", called CPC2, (le yield of ethylene product) which is its capacity to product ethylene under standard steam cracking conditions In addition, other important properties are its cokability Index ( 1 C), which is the potential of the feed to generate coke [under the standard conditions of steam cracking] as well as the chemical composition eg the percentages of linear paraffinic hydrocarbon (LIN%), isoparaffins ISO%), naphthenics (NAPHT%) aromatics (ARO%) and benzene (BENZ%) The NIR spectaim of a naphtha F on line to a steam cracker furnace was measured with a Fourier Transfoim specti onietei of 4800-4000cm ~ ' (with baseline at 4720cm-'), the absorbances were normalized, but not weighted The conditions of the steam cracking were % steam 40%, outlet temperature 830°C, residence time 0 5 sec A bank of data on known naphthas with their NIR spectra and properties had been previously prepared, the properties being detei mined with respect to a steam cracking under the same standard conditions By the Minimal Index Procedure and non weighting of the absorbances and Equation 2 the Minimal Index was determined as 3 7 x l O" 7 By the method of this invention the pioximitv indices between the unknown naphtha IF and the bank of known naphthas were determined, and those known naphthas with proximity indices less than 3 7 \ 10 ~ ' chosen Five such naphthas were found, designated l A- I E and the pi opeities ot naphtha l F obtained by averaging those of 1 A-1E The Table 1 1 shows the absoi ba ces foi the napthas 1 A- I F, as well as the estimated absoi bances foi I F and the pi pei ties of I A- I F, and the estimated properties of I F The table shows the accui acv of the method for the calculated properties, which are in agieement with the standard methods, any differences between the measured and estimated pi opeities (columns 3 and 4) being within the limits of reproducibihty of the reference methods Other properties of the naphtha such as average molecular weight, may be obtained in a similar way Example 2 The process of Example I was repeated with a different naphtha feed, whose analysis changed with time The method was applied as in Ex 1 (with Minimal Index The feed was analyzed on line with time and the potential yield of ethylene and the density monitored, the results ai e shown in the Table below Other properties were measured once with the following l esults, potential yield of methane 14 55% and propylene 16 08% and 4 caibon olefins 9 86%, the coking index 87 21 , mean molecular weight 88 68, and chemical composition 42 89% LIN, 34 60% of ISO, 17 62% NAP, 4 88% ARD and 1 79% of BEN O 96/11400 TABLE l. l Determination of (lie properties of a mixture of hydrocarbons lor feeding to a steam cracker TABLE 1.1 (cont) TABLE 1.1 (cont) TABLE 1.1 (cont) TABLE 1.1 (cont) TABLE 1.1 (cont) In the Table 3.8 E-08 means 3 8 x 10 "s

Claims: Case 8499( 21 Claims 1 A method of determining or pi edictmg a value P x of a property of a feed X to a steam cracking process or yield of said piocess, which method comprises measuring the absorption D j X of said feed at moie than one wavelength in the region 600- 2600 nm, comparing signals indicative of said absorptions or a mathematical function thereof with signals indicative of absoi ptions D,m or a mathematical function thereof at the same wavelengths for a numbei of standards S in a bank for which the said property or yield P is known, and choosing fiom the bank at least one standard S m with property oi yield P m having the smallest average value of the absolute difference at each wavelength i between the signal for the material and the signal for the standard S m to obtain value P x with averaging of said pioperties or yields P m when more than one standard S m is chosen 2 A method according to claim l comprising comparing absorption DJX (or a derivative thereof) with absorptions D,m oi a derivative thereof 3 A method according to claim I oi 2 wheiein the standard S m chosen for the property or yield wanted is such that m ι elation to the unknown feed X and each chosen standard S m the following function is met i xm /ΣD lx <expeι ιmental ei roi wherein ι xm is the proximity index and is defined by ι^( xm ) = ∑(Dιχ-D ιm )2 and the experimental error is in determining said property or yield in the standard 4 A method according to claim 3 wheiein the proximity index is less than the minimal index ι m which has been detei mined from the standards Sa, Sb, Sc by (a) calculating for each pair of standards Sa/Sb, Sa/Sc the value of ι^(a,b) etc, (b) relating the values of i^(a,b) etc to the corresponding differences EP (ab) in properties Pa, Pb etc (c) calculating for each value L for which \ L (ab) is <L, the average of the corresponding differences EPab, (cl) calculating Minimal index from the value of minimal index ι-(ab) wheie eι age EPab is the same as l epioducibi ty standard for the property 5 A method according to any one ot the pi eceding claims, wherein the properties of synthetic standards, which aie mixtures, and then spectra for consideration for possible choice foi S m are estimated fiom existing standards in the bank for which, in respect of each existing standaid for use m said mixture equation (4) and (5) are met, MinCj - ui < C,, < Max Cj + uj (4) wherein Cj, is fraction of component j in the sample i Min Cj is the minimum ofj in the standaids m the bank or in the samples for which the method is to be used Max Cj is the maximum ofj in the standaids in the bank or in the samples for which the method is to be used and uj is between l 0 and 0 05 6 A method accoi ing to any one ot claims 2-^ whei ein pi operties of Standards and spectra for consideiation toi possible choice aie estimated by interpolation from measured properties of Standaids and specti toi which the pioximity index with respect to the unknown X is not mote than I times the Minimal Index 7 A method according to claim 5 oi 6 wherein at least one ot (i) the estimated Standards and the corresponding spectia and (n) the property P x of the unknown material and its spectrum, ai e added to the bank 8 A method accoiding to any one ot the pieceding claims wheiein the property P or yield is compaied to the desued value and any deviations used in a closed loop control system to control the pi cessing equipment in i elation to a process for which the material is a teed 9 A method accoi ding to any one of the preceding claims wherein the property is a physicochemical property of feed X 10 A method according to any one of the pieceding claims which comprises determining at least one property of said feed in particular the percentage of at least one component thereof, its density and/or its mean molecular weight 1 1 A method for adding an exti synthetic standaid to a bank of known standards, each of which relates at least one absoi ption in the 600-2600nm region (or signal indicative thereof oi of a mathematical function ot said absot ption) ot a known material to a known piopeity related to that material, wheiein said property is of said material, which is a hydrocai bon teed to a steam cracking process, or is a yield of said process which method comprises choosing ti om the bank at least 2 standards for which equations 4 and 5 (as defined in claim 4) ai e met, considering mixing the chosen standards in at least one proportion to pioduce at least cue mixture for use in a synthetic standard, and estimating the spectrum of said mixture according to equation 6 and estimating a property of said mixtuie according to equation 7 and then adding the spectrum and property of each "mixtuie" to the bank , and using them in at least one model involving a correlation/regiession approach to relate NIR spectra to at least one said property 12 A method according to any one of claims 1-7 and 1 1, which comprises determining from the absorptions of the feed to said process, the yield of at least one product from the process 13 A method according to claim 12 which comprises determining the yield of at least one of alkene of 2-4 carbons, methane and coke (as expressed by a cokability index) 14 A method accoiding to any one ot the preceding claims which is computer implemented 15 A computer implemented method according to claim 14 involving a spectrometer linked to a process line containing a material X which is a feed to a steam cracking process, a computer linked to the spectrometer, and a controller linked to the computer and the process line the computei including databanks having stored therein signals indicative of absorptions of standai matei ials (or mathematical functions thereof) and corresponding pioperties ot said mateiials, or products of said process of which X is a feed oi yield of said process, the method compiises steps of measuring absorption at more than one wavelength in the region 600-2600nm at the process line and producing absoi ption signals or mathematical functions thereof by the spectrometer in accordance therewith, accessing the databanks of the computer in accordance with the absorption signals (or functions thereof), comparing, by the computer, the absorption signals (or functions thereof) to the absorptions (or functions thereof) of the standard materials stored in the databanks, choosing at least one standard based on the comparing, said standard having the smallest average value of the absolute diffeience at each wavelength I between the signal (or function thereof) foi the material and the signal (or function thereof) for the standards, with averaging of said pioperties oi yields when more than one standard is chosen, and controlling said process line in accoi dance with the outputted property/yield 16. A computer programmed to perform the method of any of claims 1-15. 17. Apparatus suitable for use in the method of any one of claims 1-15 which comprises an NIR spectrometer receiving at least one signal from a feed or product line in said process and being coupled to a computer to effect continuous measurement of the spectrum of the feed and provide control of the process.

Inventors: DESCALES BERNARD, GRANZOTTO CLAUDE, LAMBERT DIDIER, LLINAS JEAN-RICHARD, MARTENS ANDRE, DESCALES, BERNARD, GRANZOTTO, CLAUDE, LAMBERT, DIDIER, LLINAS, JEAN-RICHARD, MARTENS, ANDRE
Applicants: BP CHEM INT LTD, DESCALES BERNARD, GRANZOTTO CLAUDE, LAMBERT DIDIER, LLINAS JEAN RICHARD, MARTENS ANDRE, NAPHTACHIMIE SA, BP CHEMICALS LIMITED, DESCALES, BERNARD, GRANZOTTO, CLAUDE, LAMBERT, DIDIER, LLINAS, JEAN-RICHARD, MARTENS, ANDRE, NAPHTACHIMIE S.A.
Application year: 1995
Main classification: G01N33
IPC codes: C10G9, G01N21, G01N33