ESTIMATION OF THE TEMPERATURE OF A STEEL PRODUCT
20240118139 ยท 2024-04-11
Inventors
Cpc classification
G01J5/0003
PHYSICS
International classification
G01J5/52
PHYSICS
Abstract
A method for estimating the temperature of a steel product including a calibration step wherein the intensities at 5 wavelengths ranging from 0.9 to 2.1 ?m are recorded for several measurement condition and spectral attenuation coefficients are computed, a measurement step wherein the intensities at said 5 wavelengths are recorded and spectral attenuation coefficients are computed for several temperatures and a comparison step wherein a probability test is performed to estimate the steel product temperature.
Claims
1-12. (canceled)
13. A method for estimating a temperature T.sub.REAL of a steel product, having a temperature from 300? C. to 1600? C., the method comprising: A. a calibration step including the steps of i. measuring intensities (I), at 5 wavelengths (?) ranging from 0.9 to 2.1 ?m, wherein one is from 0.91 ?m to 1.35 ?m, one is from 1.35 ?m to 1.55 ?m, one is from 1.55 ?m to 1.85 ?m, one is from 1.85 ?m to 2.05 ?m and one is from 2.05 ?m to 2.1 ?m, by a sensor, of the radiation emitted by a reference having a known temperature (T.sub.REF) in measurement conditions characterized by an emissivity of the reference (?.sub.REF) and a transmittance of a medium between the reference and the sensor (?.sub.REF), wherein the reference is a steel product, ii. computing a spectral attenuation coefficient C.sub.CALIB using the measured intensities (I) at the 5 wavelengths,
14. The method as recited in claim 13 wherein the cooling treatment is performed during or following a hot rolling and the steel product has a temperature from 300? C. to 1100? C. and wherein in step B, T.sub.J ranges from 300? C. to 1100? C.
15. The method as recited in claim 13 wherein the cooling treatment is performed during or following a continuous casting and said steel product has a temperature from 800? C. to 1600? C. and wherein in step B, T.sub.J ranges from 800? C. to 1600? C.
16. The method as recited in claim 13 wherein in steps A)i. and B)i., the radiation intensities of 8 wavelengths (?) ranging from 0.9 to 2.1 ?m, wherein one is from 0.9 ?m to 1.11 ?m, one is from 1.11 ?m to 1.15 ?m, one is from 1.15 ?m to 1.35 ?m, one is from 1.35 ?m to 1.55 ?m, one is from 1.55 ?m to 1.85 ?m, one is from 1.85 ?m to 2.05 ?m, one is from 2.05 ?m to 2.07 ?m and one is from 2.07 ?m to 2.1 ?m are measured and in steps A)ii. and B)ii., spectral attenuation coefficients for the 8 wavelengths are computed.
17. The method as recited in claim 16 wherein in steps A)i. and B)i., the radiation intensities at 5 additional wavelengths ranging from 0.9 to 2.1 ?m are measured and in steps A)ii. and B)ii., spectral attenuation coefficients for the 8 wavelengths and the 5 additional wavelengths are computed.
18. The method as recited in claim 16 wherein in steps A)i. and B)i., the radiation intensities at 42 additional wavelengths ranging from 0.9 to 2.1 ?m are measured and in steps A)ii. and B)ii., spectral attenuation coefficients for the 8 wavelengths and the 42 additional wavelengths are computed.
19. The method as recited in claim 16 wherein in steps A)i. and B)i., the radiation intensities at 92 additional wavelengths ranging from 0.9 to 2.1 ?m are measured and in steps A)ii. and B)ii., spectral attenuation coefficients for the 8 wavelengths and the 92 additional wavelengths are computed.
20. The method as recited in claim 13 wherein N.sub.CALIB is an integer from 2 to 1000.
21. The method as recited in claim 20 wherein N.sub.CALIB is an integer from 20 to 1000.
22. The method as recited in claim 13 wherein in step C)i., the probability test includes a dimensionality reduction on said C.sub.CALIB defining main components.
23. The method as recited in claim 22 wherein in step C)i., the dimensionality reduction is performed with a principal component analysis.
24. The method as recited in claim 13 wherein in step C)i., the probability test includes a projection of the C.sub.CALIB in a probabilistic model.
25. The method as recited in claim 24 wherein in step C)i., said probabilistic model is a Gaussian mixture model.
26. A method for cooling a steel strip comprising employing the method as recited in claim 13.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
[0024] Other characteristics and advantages will become apparent from the following detailed description of the invention.
[0025]
[0026]
[0027]
[0028]
[0029]
DETAILED DESCRIPTION
[0030] The invention relates to a method for estimating the temperature T.sub.REAL of a steel product, having a temperature from 300? C. to 1600? C. comprising:
[0031] A. A calibration step comprising the steps of [0032] i. Measuring intensities (I), [0033] at 5 wavelengths (?) ranging from 0.9 to 2.1 ?m, wherein one is from 0.9 ?m to 1.35 ?m, one is from 1.35 ?m to 1.55 ?m, one is from 1.55 ?m to 1.85 ?m, one is from 1.85 ?m to 2.05 ?m and one is from 2.05 ?m to 2.1 ?m, by means of a sensor, [0034] of the radiation emitted by a reference having a known temperature (T.sub.REF) in measurement conditions characterised by an emissivity of said reference (?.sub.REF) and a transmittance of a medium between said reference and said sensor (?.sub.REF), wherein said reference is a steel product, [0035] ii. Computing a spectral attenuation coefficient C.sub.CALIB using said measured intensities (I) at said 5 wavelengths,
[0038] B. A measurement step comprising the steps of [0039] i. Measuring intensities of the radiation emitted by said steel product, I, at said 5 wavelengths (?) ranging from 0.9 to 2.1 ?m, [0040] ii. Computing N.sub.T spectral attenuation coefficients C.sub.COMPUTETj, for N.sub.T temperatures (Tj) ranging from 300 to 1600? C. and for said 5 wavelengths, N.sub.T being an integer from 2 to 1300,
[0043] C. A comparison step comprising the steps of [0044] i. Performing a probability test for finding the most likely C.sub.COMPUTETj among the C.sub.CALIB [0045] ii. Estimating the temperature of said steel product, T.sub.REAL as being equal to the temperature T.sub.J of said most likely C.sub.COMPUTET.sub.J.
[0046] The steps of the process are illustrated in
[0047] The steel product can be any type of product such as a strip, a band or a slab. The temperature of the steel product is not known however, depending on the process steps wherein the measurement is performed, the skilled in the art knows the range in which the temperature should be. For example, in a cooling treatment of a steel strip following a hot rolling, the temperature of the steel strip is generally comprised between 300? C. and 1100? C.
[0048] In the calibration step A.i., the intensities of the radiation emitted by the reference steel product can be measured by any suitable means. They are preferably measured by a hyperspectral camera.
[0049] In the calibration step A.i., the reference steel product has preferably a composition similar to the steel product having its temperature estimated. Even more preferably, the reference steel product has the same grade as the steel product being analysed.
[0050] In the calibration step A.i., the temperature of the reference steel product can be measured by any means. Preferably, this temperature is measured using thermocouples.
[0051] In the calibration step A.ii., a spectral attenuation coefficient C.sub.CALIB can be computed by dividing each of the recorded intensities I by P(?, Tj) at the measured temperature.
[0052] The step A.i. and A.ii. are repeated for N.sub.CALIB times, each time corresponding to new measurement conditions. Measurement conditions are defined by a combination of the reference emissivity (?.sub.REF) and the transmittance of the media between the reference and the sensor (?.sub.REF) It permits to get various C.sub.CALIB for different measurement conditions. The greater N.sub.CALIB, greater is the accuracy of the estimation.
[0053] In other words, the steps i. and ii. of the calibration step A. are repeated for various combinations of steel product emissivity (?.sub.REF) and measurement condition transmittance (?.sub.REF). ?.sub.REF varies as a function of several factors such as the temperature of the reference steel product, the surface properties (such as presence of oil). ?.sub.REF depends on the medium between the sensor and the reference steel product such as the thickness of a water layer on the reference steel product.
[0054] However, the value of ?.sub.REF and ?.sub.REF do not need to be known to compute the C.sub.CALIB.
[0055] In the measurement step B.i., the intensities emitted by said steel product are preferably measured using a hyperspectral camera. The wavelength at which the intensities are measured in step B.i. are the same as the one in step A.i.
[0056] In the measurement step B.ii., greater is the number of temperature (N.sub.T) for which a spectral attenuation coefficient is computed, greater is the accuracy of the estimation.
[0057] Moreover, in both the calibration and the measurement steps, the measured intensities can be adjusted using a transfer function.
[0058] Measuring the radiation intensities at said 5 wavelengths, one from 0.9 ?m to 1.35 ?m, one from 1.35 ?m to 1.55 ?m, one from 1.55 ?m to 1.85 ?m, one from 1.85 ?m to 2.05 ?m and one from 2.05 ?m to 2.1 ?m permits to describe the spectrum shape of a combination of the product emissivity and of the medium transmittance using a reduced number of intensities. Indeed, in a range from 0.9 to 2.1, the absorption spectrum of water, as represented in
[0059] In the comparison step C., a probability test is performed to find the most likely C.sub.COMPUTEDT.sub.J among the C.sub.CALIB. Any method permitting to find the most likely C.sub.COMPUTEDT.sub.J among the C.sub.CALIB can be performed.
[0060] This method permits to improve the accuracy of the temperature estimation of a steel product because the estimation considers the influence of the measurement condition, such as the presence of water on the steel product.
[0061] Preferably, said cooling treatment is done during or following a hot rolling and said steel product has a temperature from 300? C. to 1100? C. and wherein in step B, T.sub.J ranges from 300? C. to 1100? C. This cooling is usually performed in a run-out-table.
[0062] Preferably, said cooling treatment is done during or continuous casting and said steel product has a temperature from 800? C. to 1600? C. and wherein in step B, T.sub.J ranges from 800? C. to 1600? C.
[0063] Preferably, in steps A)i. and B)i., the radiation intensities of 8 wavelengths (?) ranging from 0.9 to 2.1 ?m, wherein one is from 0.9 ?m to 1.11 ?m, one is from 1.11 ?m to 1.15 ?m, one is from 1.15 ?m to 1.35 ?m, one is from 1.35 ?m to 1.55 ?m, one is from 1.55 ?m to 1.85 ?m, one is from 1.85 ?m to 2.05 ?m, one is from 2.05 ?m to 2.07 ?m and one is from 2.07 ?m to 2.1 ?m are measured and in steps A)ii. and B)ii., spectral attenuation coefficients for said 8 wavelengths are computed. Measuring the radiation intensities of those said 8 wavelengths improves the accuracy of the measurement when a medium, such as a vapor, is present between the steel product and the sensor. The ranges, from 1.11 ?m to 1.15 ?m and from 2.05 to 2.07 ?m corresponds to peaks of the absorption spectrum of vapor, as illustrated in
[0064] Preferably, in steps A)i. and B)i., the radiation intensities at 5 additional wavelengths ranging from 0.9 to 2.1 ?m are measured and in step steps A)ii. and B)ii., spectral attenuation coefficients for said 8 wavelengths and said 5 additional wavelengths are computed. The 13 wavelengths are preferably evenly distributed in said ranges which means that the 13 wavelengths are spaced by an interval of 0.1 ?m (0.9, 1.0, . . . , 2.0, 2.1).
[0065] Preferably, in steps A)i. and B)i., the radiation intensities at 42 additional wavelengths ranging from 0.9 to 2.1 ?m are measured and in step steps A)ii. and B)ii., spectral attenuation coefficients for said 8 wavelengths and said 42 additional wavelengths are computed.
[0066] Preferably, in steps A)i. and B)i., the radiation intensities at 92 additional wavelengths ranging from 0.9 to 2.1 ?m are measured and in step steps A)ii. and B)ii., spectral attenuation coefficients for said 8 wavelengths and 92 additional wavelengths are computed.
[0067] Preferably, N.sub.CALIB is an integer from 2 to 1000 and preferably from 20 to 1000.
[0068] Preferably, in step C)i., the probability test comprises a dimensionality reduction on said C.sub.CALIB defining main components. Even more preferably, in step C)i., the dimensionality reduction is performed with a principal component analysis.
[0069] Preferably, in step C)i., the probability test comprises the projection of said C.sub.CALIB in a probabilistic model. Even more preferably, in step C)i., said probabilistic model is a Gaussian mixture model.
Embodiment of the Invention
[0070] In order to assess the accuracy of the claimed method, the temperature of a steel strip being cooled in a run-out table has been measured by a thermocouple and estimated by the present method.
A. Calibration Step
[0071] During the calibration step, the radiation intensity emitted by a steel strip, being cooled in a run-out table, has been recorded by a hyperspectral camera. More precisely, the intensities at 256 wavelengths ranging from 1.1 to 2.1 ?m have been recorded. The 256 wavelengths were evenly distributed in said range. The temperature of the steel strip has been measured by thermocouples.
[0072] The calibration steps comprised 6 trials lasting between 5 and 15 minutes wherein the radiation intensities were measured between 10 and 100 times per second. So, few thousands measurement conditions were recorded. They took place during the cooling (presence of droplet, fog and water) and after/prior the cooling (no presence of water/fog) and for steel strip having various oxidation level.
[0073] Then, several spectral attenuation coefficient C.sub.CALIB have been computed using the spectral density of electromagnetic radiation emitted by a black body in thermal equilibrium, the measured temperature by the thermocouples and the transfer function of the hyperspectral camera. Each of the C.sub.CALIB had 256 values, one for each of said 256 wavelengths. The C.sub.CALIB are plotted in
B. Measurement
[0074] During the measurement step, the radiation intensity emitted by a steel strip, being cooled in a run-out table, has been recorded by a hyperspectral camera. More precisely, the intensities, at the same 256 wavelengths as in the calibration step, have been recorded.
[0075] Moreover, 261 temperatures ranging from 300 to 1600? C. have been defined as possible temperatures for the measured steel strip, so N.sub.T=261. Those temperatures were spaced by an interval of 5? C.: (300, 305, 310 . . . 1595, 1600? C.).
[0076] Then, the measured radiation intensities at said 256 wavelengths have been divided by the value of the Planck's law and multiplied by the transfer function for each of the 261 temperatures defined previously. So, 261 spectral attenuation coefficients C.sub.COMPUTET.sub.J were computed. Each of said 261 spectral attenuation coefficients had 256 values, one for each measured wavelength, and a temperature associated.
C. Comparison
[0077] The probability test of this embodiment comprises the following steps [0078] 1) A principal component analysis (PCA) is performed on all the spectral attenuation coefficients C.sub.CALIB, obtained during the calibration step, to get 3 main PCA components (PC.sub.1, PC.sub.2, PC.sub.3). Then, each C.sub.CALIB of the database is approximated using said three main components. [0079] 2) The spectral attenuation coefficients are projected into a PCA space defined by the 3 main components previously defined. [0080] 3) The density of points in the PCA space of step 2) is projected in a Gaussian mixture model. [0081] 4) Each of said 261 C.sub.COMPUTET.sub.J is projected into the PCA space defined by the 3 main components previously defined. [0082] 5) For each of the 261 C.sub.COMPUTET.sub.J, an approximated spectral attenuation coefficient is reconstructed from its coordinates in the PCA space defined in 4) and a first likelihood factor (Like1) is defined by comparing the original spectral attenuation and its reconstruction from the PCA coefficients. [0083] 6) Each of said 261 C.sub.COMPUTET.sub.J is projected in the Gaussian mixture model defined in 3) and a second likelihood factor (Like2) is determined. Like2 represents the local density of the learning data projected in the PCA space. [0084] So, two likelihoods (Like1 and Like2) are associated to each of said 261 C.sub.COMPUTET.sub.J. [0085] 7) Then, the product of Like1 and Like2 is calculated for each of said 261 C.sub.COMPUTET.sub.J.
[0086] Finally, the temperature associated to the C.sub.COMPUTET.sub.J maximizing the product of 7) is defined as being the temperature of the steel strip.
[0087] Comparative results of the temperature obtained through the thermocouples measures and via the claimed estimation are plotted in