CN102749141A - Radiation temperature measuring method and apparatus for measuring true target temperature - Google Patents

Radiation temperature measuring method and apparatus for measuring true target temperature Download PDF

Info

Publication number
CN102749141A
CN102749141A CN201210266196XA CN201210266196A CN102749141A CN 102749141 A CN102749141 A CN 102749141A CN 201210266196X A CN201210266196X A CN 201210266196XA CN 201210266196 A CN201210266196 A CN 201210266196A CN 102749141 A CN102749141 A CN 102749141A
Authority
CN
China
Prior art keywords
temperature
wavelength
target
warm
approximate
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201210266196XA
Other languages
Chinese (zh)
Inventor
肖功弼
张文生
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Automation of Chinese Academy of Science
Original Assignee
Institute of Automation of Chinese Academy of Science
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Institute of Automation of Chinese Academy of Science filed Critical Institute of Automation of Chinese Academy of Science
Priority to CN201210266196XA priority Critical patent/CN102749141A/en
Publication of CN102749141A publication Critical patent/CN102749141A/en
Priority to CN201210417758.6A priority patent/CN103048050B/en
Pending legal-status Critical Current

Links

Images

Abstract

The invention discloses a radiation temperature measuring method for measuring true target temperature, which comprises the steps of: obtaining measured brightness temperature Tsi corresponding to three wavelengths and the approximate true temperature T of a target through a three-wavelength radiation temperature measuring apparatus, wherein the three wavelengths are equally spaced; determining the emissivity of epsilon lambda 1, epsilon lambda 2 and epsilon lambda 3 of the target at the three wavelengths according to the measured brightness temperature Ts1, Ts2, Ts3 and the numerical value of the approximate true temperature T; computing brightness temperature Ts1', Ts2' and Ts3' of the target corresponding to the three wavelengths according to the computed emissivity values of the target at the three wavelengths and the approximate true temperature T; computing the sum sigma delta T2 of mean square errors according to the measured brightness temperature and the computed brightness temperature; further computing the emissivity and the brightness temperature value by using the varied approximate true temperature value; further computing the mean square errors of the measured brightness temperature and the computed brightness temperature; comparing the sizes of the two mean square errors; and finding the practical true temperature value T until the sum of the mean square errors is minimum. By using the method provided by the invention, the true temperature of the target can be measured.

Description

A kind of radiation temperature measurement method and instrument of measurement target true temperature
Technical field
The present invention relates to utilize radiation to carry out the technical field of thermometric, particularly a kind of under the situation of contact target not the method and the instrument of measurement target true temperature.
Background technology
The thermal radiation temperature meter is through measuring target emanation intensity, measuring the thermometer of target temperature.Its theoretical foundation is about kirchhoff (kirchhoff) law of object radiation intensity and absorptive power relation (absorptive power of target equals its emissivity when the local thermodynamic equilibrium), and about Planck (Planck) law of blackbody radiation.Radiation temperature measurement has many good qualities, and (for example, it does not measure ceiling restriction, the temperature field of not disturbance measured target and thermal equilibrium; Its measuring accuracy is high, and dynamic response is good or the like); Become the main method of non-contact temperature measuring.Countries in the world have the radiation thermometer product of a lot of models at present, and are used widely.Various countries' high temperature mete-wand and standard thermometer all are radiation thermometer (and all being the brightness temperature meters).But the disadvantage of radiation thermometer is, the temperature of the realistic objective that it records is not the true temperature (being called for short very temperature) of target, but the apparent temperature of target.Apparent temperature mainly comprises three kinds of brightness temperature, radiation temperature and colour temperatures (or claiming colour temperature to be called for short colour temperature) at present.Corresponding these three kinds of apparent temperatures three kinds of basic radiation temperature measurement method and instrument of can extending out; Be brightness method temperature instrumentation, total radiation temperature instrumentation (abbreviation radiation thermometer) and colourimetry temperature instrumentation.Wherein brightness method temperature instrumentation is an instrument most important, the most universal in the radiation temperature measurement.It is in theory very strict after introducing the effective wavelength notion, has become various countries' high temperature mete-wand and standard meter, and is playing a leading role for a long time aspect metering transmission of quantity value and the commercial Application.
The thermal radiation temperature meter all is to be foundation with the blackbody furnace in calibrating (calibration), black matrix is carried out calibration handle.But when measuring realistic objective (non-black-body) temperature, what obtain is the apparent temperature (brightness temperature, radiation temperature or colour temperature) of target.Know the true temperature of target, then will know the emissivity (λ .T) of target.The emissivity of target is the function of the temperature of target and wavelength still not, but also relevant with factors such as the material character of target, surface state.So just big difficulty and uncertain is come in the true temperate zone of confirming target.In the last hundred years, this significant drawback of radiation temperature measurement is perplexing thermometric circle always.After the eighties in 20th century, multinational scientist has carried out big quantity research to the multi-wavelength thermometer, and attempt is through using the very temperature that methods such as multi-wavelength thermometer record realistic objective.Begin further to develop six wavelength thermometers from the three-wavelength thermometer, in addition the thermometer of 35 wavelength (specifically can be referring to Li Jilin. Xiao Gongbi. Yu Lunpeng writes; " radiation temperature measurement and calibrating/collimation technique "; The .2009 of China Measuring Press publishes August; P105-P106.) etc. all develop.But up to the present, these multi-wavelength thermometers all can not be measured the true temperature of target.
Summary of the invention
The present invention is in order to overcome the above-mentioned defective that prior art exists; Promptly to the problem that can not measure the target true temperature; Propose a kind of specific three-wavelength radiation temperature measurement instrument and utilize this instrument to carry out the method for thermometric; The foundation warm mathematical model of looking for the truth, and through high speed interative computation method, the true temperature of The real time measure target.
The radiation temperature measurement method of the measurement target true temperature that the present invention proposes comprises: step 1, obtain the corresponding bright temperature T of actual measurement of three wavelength through the three-wavelength Radiation Temperature Measurement Instrument SiAnd the approximate very warm T of target, wherein the value of i is 1,2 or 3, three wavelength for uniformly-spaced; Step 2 is according to the bright temperature Ts of actual measurement 1, Ts 2, Ts 3And the numerical value of approximate very warm T, confirm the emissivity of target under three wavelength λ 1, ε λ 2, ε λ 3Step 3 according to the emissivity value of target under three wavelength and the approximate very warm T value that step 2 calculates, is calculated the bright temperature Ts of corresponding three wavelength of target 1', Ts 2', Ts 3'; Step 4 according to the bright temperature of calculating in bright temperature of the actual measurement in the step 1 and the step 3, is tried to achieve mean square deviation sum ∑ Δ T 2As follows: ∑ Δ T 2=(Ts 1-Ts 1') 2+ (Ts 2-Ts 2') 2+ (Ts 3-Ts 3') 2 Step 5, the approximate very warm T ' of order=T+1, and according to the bright temperature Ts of actual measurement 1, Ts 2, Ts 3Calculate the emissivity of three wavelength of target according to the mode of step 2 λ 1', ε λ 2', ε λ 3'; Step 6, the emissivity that calculates according to step 5 λ 1', ε λ 2', ε λ 3' reach the value of T '=T+1, calculate three brightness temperature value Ts of target according to the mode of step 3 1", Ts 2", Ts 3"; Step 7 is according to calculating bright temperature Ts 1", Ts 2", Ts 3" value and the bright temperature Ts of actual measurement 1, Ts 2, Ts 3, try to achieve the mean square deviation sum and be: (∑ Δ T 2) '=(Ts 1-Ts 1") 2+ (Ts 2-Ts 2") 2+ (Ts 3-Ts 3") 2If step 8 is (∑ Δ T 2) '<(∑ Δ T 2), then approximate very warm T '=T+1 than T more near the very warm T of reality Very, the approximate very warm T of order then "=T '+1=T+2, continue repeating step 5-7, until the mean square deviation sum hour, promptly think and find out the very warm T of reality VeryIf value is (∑ Δ T 2) '>(∑ Δ T 2), then approximate very warm T '=T+1 is away from the very warm T of reality Very, realistic border is very warm in the other direction is worth, even T '=T-1, continuing set by step, 5-7 calculates ∑ Δ T 2, until the mean square deviation sum hour, promptly find very temperature value T of reality Very
The present invention also provides a kind of radiation temperature measurement instrument of measurement target true temperature, and this temperature measuring device is the three-wavelength temperature measuring device, and it comprises: storer, the bright temperature scale of three wavelength of storage is decided curve Tsi=f (I λ i) corresponding numerical tabular.And the numerical tabular of
Figure BDA00001949229900031
calibration curve.The supercomputing unit, it further comprises: first module is used for the bright temperature Ts of actual measurement that obtains according to from storer 1, Ts 2, Ts 3And the numerical value of approximate very warm T, confirm the emissivity of target under three wavelength λ 1, ε λ 2, ε λ 3Unit second according to the emissivity value of target under three wavelength and the approximate very warm T value that first module calculates, calculates the bright temperature Ts of corresponding three wavelength of target 1', Ts 2', Ts 3'; Mean square deviation sum ∑ Δ T according to surveying the bright temperature of calculating that try to achieve bright temperature and Unit second, is tried to achieve in Unit the 3rd 2As follows: ∑ Δ T 2=(Ts 1-Ts 1') 2+ (Ts 2-Ts 2') 2+ (Ts 3-Ts 3') 2Unit the 4th, the approximate very warm T ' of its order=T+1, and according to the bright temperature Ts of actual measurement 1, Ts 2, Ts 3Calculate the emissivity of three wavelength of target according to the mode of first module λ 1', ε λ 2', ε λ 3'; Unit the 5th, the emissivity that calculates according to Unit the 4th λ 1', ε λ 2', ε λ 3' reach the value of T '=T+1, calculate three brightness temperature value Ts of target according to the mode of Unit second 1", Ts 2", Ts 3"; Unit the 6th is according to calculating bright temperature Ts 1", Ts 2", Ts 3" value and the bright temperature Ts of actual measurement 1, Ts 2, Ts 3, try to achieve the mean square deviation sum and be: (∑ Δ T 2) '=(Ts 1-Ts 1") 2+ (Ts 2-Ts 2") 2+ (Ts 3-Ts 3") 2If Unit the 7th is (∑ Δ T 2) '<(∑ Δ T 2), then approximate very warm T '=T+1 than T more near the very warm T of reality Very, the approximate very warm T of order then "=T '+1=T+2, continue the processing procedure of repetition Unit the four-six, until the mean square deviation sum hour, promptly think and find out the very warm T of reality VeryIf value is (∑ Δ T 2) '>(∑ Δ T 2), then approximate very warm T '=T+1 is away from the very warm T of reality Very, realistic border is very warm in the other direction is worth, even T '=T-1 continues the processing procedure of repetition Unit the four-six and calculates ∑ Δ T 2, until the mean square deviation sum hour, promptly find very temperature value T of reality Very
Description of drawings
Fig. 1 is the radiation temperature measurement instrument structural drawing of measurement target true temperature of the present invention;
Fig. 2 is the software flow pattern of the radiation temperature measurement method of measurement target true temperature of the present invention.
Embodiment
For making the object of the invention, technical scheme and advantage clearer, below in conjunction with specific embodiment, and with reference to accompanying drawing, to further explain of the present invention.
Fig. 1 is the radiation temperature measurement instrument structural drawing of measurement target true temperature of the present invention.
Temperature measuring device of the present invention is three-wavelength (be three the narrow wave bands in fact) thermometer of specific wavelength, and three specific wavelengths are chosen according to following principle: (1) three wavelength be spaced apart equidistance, and at interval should not be too big.The correctness (seeing after) of the unique supposition of the too big the present invention of influence in interval.As in the near-infrared region, to be advisable about 0.1~0.3 μ m (micron).For example get λ 1=0.8 μ m, λ 2=0.9 μ m, λ 3Three wavelength of=1.0 μ m come thermometric.Certainly at interval also should not be too little, because of too hour wavelength differences is little at interval, influence measuring accuracy.Wherein this principle is the most important, is the key distinction place of the present invention and prior art, and other three principles are identical with general thermometer.(2) absorb measuring influence in order to reduce intermediate gas, above-mentioned three-wavelength should be avoided water vapor common in the air, CO 2Deng the gas molecule absorption band.(3) in order to improve measuring accuracy, above-mentioned three-wavelength should be chosen in the higher zone of photovalve spectrum sensitivity curve medium sensitivity commonly used.For example above-mentioned choose 0.8~1.0 between three wavelength, promptly choose the higher zone of silicon photocell sensitivity commonly used.(4) for convenient data processing, target light change in radiation intensity scope (corresponding to target temperature range) should be chosen in the photovalve linear zone working range.Therefore the radiation intensity of the target when the photocurrent signal of the photovalve of thermometer output and this wavelength is proportional.
With reference to Fig. 1, the chopper wheel 5 that the radiation temperature measurement instrument of measurement target true temperature of the present invention (following abbreviation temperature measurer) comprises the catoptron 2, aiming catoptron 3, target imaging lenslet 4 of object lens (adjustable focus) 1, band central small hole, driven by motor M, measuring light fulgurite 6, LED 7, control with synchronizable optical fulgurite 8, observation with eyepiece 9.
Object lens 1 adjustable focus of temperature measurer, with target imaging in the band central small hole catoptron 2 on.After catoptron 3 reflect again, supply eye-observations aiming usefulness by the light of catoptron 2 reflection through eyepiece 9 imagings.Behind the light transmission catoptron 2 of the central small hole of catoptron 2, by target imaging lenslet 4 focal imagings on measuring light fulgurite 6.Therefore the diameter of the central small hole of catoptron 2 determines the diameter of measurement point; The focal distance f of object lens 1 is divided by the central small hole diameter d, and promptly the numerical value of f/d is the important indicator of temperature measurer, is called the distance coefficient of temperature measurer.Behind the light process target imaging lenslet 4 that thermometric is used, need the chopper wheel 5 that drives through by modulating motor M.Three narrow band pass filters are housed on the chopper wheel 5, let λ respectively 1, λ 2, λ 3The light timesharing of three wavelength is radiated on the measuring light fulgurite 6 through chopper wheel 5.Measuring light fulgurite 6 produces the photocurrent signal that is directly proportional with target emanation intensity through photoelectric effect.An aperture is arranged on the chopper wheel 5, and when the light transmission aperture of the LED 7 that is installed in chopper wheel 5 left sides, short irradiation is on synchronously with photoelectric tube 8.The synchronizable optical fulgurite produces synchronizable optical current impulse signal through photoelectric effect to the LED light that sees through chopper wheel, and this pulse signal is used to distinguish the order of three wavelength.When chopper wheel 5 revolutions move a week, produce a synchronization pulse, distinguish three signals of measuring wavelength when synchronization pulse is used for signal Processing.The signal of measuring light fulgurite 6 is input to storer with synchronization pulse jointly after prime amplifier and A/D conversion.
Because of prior (with blackbody furnace calibration temperature measurer the time) in the storer stored the demarcation pitch curve and R=f (T) curve (explanation of seeing after) of three wavelength; So can find out the bright temperature T of corresponding three wavelength rapidly through storer S1, T S2, T S3And the numerical value of approximate very warm T.With T S1, T S2, T S3The numerical value that reaches T is input to the supercomputing unit, by interative computation flow process shown in Figure 2, can calculate the true temperature T of target fast VeryNumerical value.Adopt modern computing technique fast, can in 100 milliseconds, accomplish interative computation, can think to measure in real time true temperature.
In the structure of temperature measurer shown in Figure 1, three wavelength of gating rotate through chopper wheel 5 and realize.In fact three wavelength of gating can also be through prismatic decomposition, and method such as beam divider realizes.The method that realizes three wavelength of gating is a lot, and it does not influence realization of the present invention and effect.
When demarcating temperature measurer with blackbody furnace, the photoelectricity flow valuve I of the photovalve output of corresponding certain wavelength λFor:
I λ = k λ M ( λ . T ) = k λ C 1 λ - 5 [ e C 2 / λT - 1 ] - 1 - - - ( 1 )
In the formula: I λPhotoelectricity flow valuve for the photovalve of corresponding wavelength λ temperature measurer; k λPhotovalve photoelectric conversion factors (wavelength confirms that the back is a constant) for corresponding wavelength λ temperature measurer; M (λ .T) is a planck formula, C in the formula 1, C 2Be respectively first radiation constant and second radiation constant; T is the blackbody furnace temperature; λ is a wavelength.At most of practical occasion C 2>>λ T can enough accurately substitute planck formula with Wien formula, and then above-mentioned (1) formula is reduced to:
I λ = k λ C 1 λ - 5 e - C 2 / λT - - - ( 1 ) '
When with the temperature of brightness temperature instrumentation amount realistic objective (non-black-body), that record is the brightness temperature Ts (being called for short bright temperature Ts) of target.And the true temperature of target (being called for short very temperature) T and its bright temperature Ts have following relation (being drawn by bright temperature definition):
M(λ.Ts)=ε λM(λ.T) (2)
With planck formula substitution following formula, another expression-form of (2) formula is arranged:
Ts = C 2 λ ln [ e c 2 / λ i T - 1 ϵ λ + 1 ] - - - ( 2 )
In the formula: T is the true temperature of target; Ts is the brightness temperature of target; ε λFor target in the fashionable emissivity of wavelength (between 0~1).By Wien formula substitution (2) formula, can draw the true kelvin relation simplified style of bright gentleness and do
Ts = C 2 T C 2 - λ T ln ϵ λ - - - ( 2 ) '
Behind selected three wavelength, then can make a three-wavelength temperature measurer by mentioned above principle.This temperature measurer can be measured three brightness temperature Ts of target simultaneously l, Ts 2, Ts 3
When above-mentioned three-wavelength temperature measurer is examined and determine with blackbody furnace, can obtain the calibration curve of three brightness temperatures simultaneously, and calibration curve is stored in the storer.If blackbody furnace calibrating point temperature is T 0The time, corresponding photovalve photocurrent is I 0(three corresponding I of wavelength difference 01, I 02, I 03Three values).When using the temperature of this temperature measurer measurement target (non-black-body), can obtain three photoelectricity flow valuve I of three wavelength simultaneously λ 1, I λ 2, I λ 3, can measure three brightness temperature Ts of three wavelength of target respectively according to these three photoelectricity flow valuves 1, Ts 2, Ts 3Value.According to planck formula can derive the formula of embodying as follows (can be referring to Li Jilin. Xiao Gongbi. Yu Lunpeng writes; " radiation temperature measurement and calibrating/collimation technique "; The .2009 of China Measuring Press publishes August; P64-P65):
Tsi = C 2 λ i ln { I 0 i I λi [ exp ( C 2 λ i T 0 ) - 1 ] + 1 } - - - ( 3 )
In the formula: three brightness temperature values that corresponding three wavelength produced when Tsi was measurement target (have three value Ts 1, Ts 2, Ts 3); λ iBe three selected wavelength X of thermometer 1, λ 2, λ 3I λ iThe photoelectricity flow valuve that produces during for measurement target, three wavelength are to there being I λ 1, I λ 2, I λ 3Three values; I 0 λWhen examining and determine, be T in blackbody temperature with blackbody furnace 0The time corresponding three wavelength the photoelectricity flow valuve (have three value I 01, I 02, I 03).
When replacing planck formula with Wien formula is approximate, then above-mentioned (3) formula is reduced to
Tsi = C 2 λ i ln I 0 i I λi + C 2 T 0 - - - ( 3 ) ,
For most solids or liquid (containing metal or nonmetal and metal oxide etc.), their emissivity with wavelength change all mild (can be referring to Li Jilin. Xiao Gongbi. Yu Lunpeng writes; " radiation temperature measurement and calibrating/collimation technique "; The .2009 of China Measuring Press publishes August; P36-P37), that is to say the emissivity of most solids or liquid, is that approximately linear changes among a small circle with wavelength change.Radiation thermometer is used for measure solid or fluid temperature mostly, seldom is used for the temperature of measurement gas.
We suppose that the emissivity of measurement target (non-black-body) is with wavelength change, at λ 1To λ 3It among a small circle linear change.This is unique supposition that the present invention derives, and this supposition all is approximate correct to most solids and liquid.According to this approximate correct supposition, can draw and work as λ 1, λ 2, λ 3Be spaced apart equidistance the time have:
ε λ1×ε λ3≈ε 2 λ2 (4)
In the formula: ε λ 1Be wavelength X 1The emissivity of corresponding target; ε λ 2Be wavelength X 2The emissivity of corresponding target; ε λ 3Be wavelength X 3The emissivity of corresponding target.
Because be the three-wavelength temperature measurer, it can measure the photoelectricity flow valuve of corresponding three wavelength of target simultaneously, corresponds to I λ 1, I λ 2, I λ 3, they can have following relation (can draw with reference to (1) formula and (2) formula:
I λ 1 × I λ 3 I 2 λ 2 = ϵ λ 1 · k λ 1 · M ( λ 1 . T ) × ϵ λ 3 · k 3 · M ( λ 3 . T ) [ ϵ λ 2 · k λ 2 · M ( λ 2 . T ) ] 2 - - - ( 5 )
In the formula: ε λ 1, ε λ 2, ε λ 3Be respectively the emissivity of corresponding three wavelength of corresponding target; k λ 1, k λ 2, k λ 3Be respectively the photoelectric conversion factors (corresponding each wavelength is constant) of corresponding three wavelength of photovalve of thermometer; M (λ 1.T), M (λ 2.T), M (λ 3.T) be respectively the radiation intensity that planck formula is temperature corresponding three wavelength when being T (T is the very temperature of target).
With approximate expression (4) substitution (5) formula, and make
Figure BDA00001949229900092
to get:
R = K λ 1 · M ( λ 1 . T ) × K λ 3 · M ( λ 3 . T ) K λ 2 2 [ M ( λ 2 . T ) ] 2 ≈ f ( T ) - - - ( 6 )
Through (6) formula can find out the R value with the spectral emittance ε of target λIrrelevant, it is the function of target true temperature T.And the R value can be examined and determine blackbody furnace through the three-wavelength thermometer fully and actually measured and be stored in the storer.And when examining and determine with blackbody furnace, (4) formula and (6) formula accurately equate (because of the emissivity of black matrix is 1).When (6) formula of derivation, adopted unique approximate correct supposition of the present invention, promptly the emissivity of target is at λ 1To λ 3Be linear change among a small circle, and derive approximate expression (4) thus.In fact the target emissivity is at λ 1To λ 3Not necessarily be linear change among a small circle, and (4) formula is approximate correct.The target of measuring by (6) so very temperature has certain error.Therefore we claim the target of measuring by (6) formula very warm approximate very warm T for target.
The reality that makes target very temperature is used T VeryExpression.Approximate very warm T and the very warm T of reality VeryHave certain error.But approximate very warm T is certainly at T VeryNeighbouring (promptly difference is not too large between the two, common used in industry with temperature range in, generally 100 the degree about).(mathematical model) so as follows is starting point with approximate very warm T, carries out interative computation and seeks very temperature of reality.
Fig. 2 is the radiation temperature measurement method flow diagram of measurement target true temperature of the present invention, and main embodiment is according to the approximate very warm T value of target among the figure, seeks the very warm T of target reality VeryThe interative computation process, going down comprises:
Step 1 is according to three bright temperature Ts of the actual target of measuring 1, Ts 2, Ts 3Numerical value (seeing (3) formula expression formula for details) and the approximate very warm T value of obtaining by (6) formula of target, it is following to calculate the emissivity of target under three wavelength according to (2) formula:
ϵ λi = M ( λ i . T si ) M ( λ i . T ) = I λi I ( λ i . T ) - - - ( 7 )
In the formula: ε λ iBe respectively wavelength X 1, λ 2, λ 3Corresponding target emissivity λ 1, ε λ 2, ε λ 3I λ iBe respectively wavelength X 1, λ 2, λ 3Corresponding target (non-black-body) the photocurrent I that spectral radiance produced λ 1, I λ 2, I λ 3I (λ i.T) be respectively wavelength X 1, λ 2, λ 3Article three, calibration curve is similar to the corresponding photocurrent of very warm T.(can from the calibration curve memory table, find).Be the temperature photocurrent I (λ that blackbody radiation intensity is produced when being T 1.T), I (λ 2.T), I (λ 3.T).
So tentatively confirmed the emissivity numerical value (also having certain error certainly) under corresponding three wavelength of target because approximate very temperature has error.
Step 2 according to the emissivity value of target under three wavelength and the approximate very warm T value that (7) formula is calculated, can be calculated the brightness temperature value Ts of corresponding three wavelength of target by (2) formula 1', Ts 2', Ts 3' as follows (above Ts, add ' represent that this brightness temperature calculates for the first time):
T si ′ = C 2 λ i ln [ e C 2 / λ i T - 1 ϵ λ i + 1 ] - - - ( 8 )
In the formula: T measures C by (6) formula for approximate very temperature 2Be the Planck second radiation constant.C 2=1.4388 * 10 4μ m.K, λ iBe wavelength (total λ 1, λ 2, λ 3Three wavelength), ε λ iBe wavelength X iThe time the target emissivity.Tentatively calculate approximate value by (7) formula.
Step 3, the difference of three brightness temperature values of target of calculating according to (8) formula and three brightness temperature values of the actual target of measuring is zero, explains that then unique supposition of the present invention and correlation computations are entirely true.Promptly approximate very warm T equals the very warm T of reality VerySome errors (containing supposition correctness error, (4) formula approximation, the substantial measurement errors and the error of calculation etc.) are in fact always arranged.Can obtain three wavelength brightness temperature mean square deviation sum ∑ Δ T by following formula 2As follows: (being least square method)
∑ΔT 2=(Ts 1-Ts 1’) 2+(Ts 2-Ts 2’) 2+(Ts 3-Ts 3’) 2 (9)
In the formula: Ts 1, Ts 2, Ts 3Be respectively three brightness temperature values of actual measurement target, Ts 1', Ts 2', Ts 3' be respectively three brightness temperature values of the target of calculating by (8) formula.
Obvious ∑ Δ T 2Approximate very warm T of more little explanation and the very warm T of reality VeryBetween error more little, as long as find ∑ Δ T 2The approximate very temperature value T that the minimum value place is corresponding has promptly found very temperature value T of reality Very
Step 4 near approximate very warm T value, is sought the very warm T of reality VeryProcess is following:
Step 401, the approximate very warm T ' of order=T+1, and three bright temperature Ts of the actual target of measuring 1, Ts 2, Ts 3The emissivity of calculating three wavelength of target by (7) formula is following,
ϵ ′ λi = M ( λ i . T si ) M ( λ ′ 1 . T ′ ) = I λi I ( λ i . T ′ ) - - - ( 7 ) ,
The symbol implication is identical with (7) formula in the formula, adds on the emissivity of calculating ' represent it is the further approximate emissivity of extrapolating, press the emissivity that (7) formula is calculated first with difference.
Step 402 is according to (7) ' ε that calculates of formula λ 1', ε λ 2', ε λ 3' value and T '=T+1 value, it is following to calculate three brightness temperature values of target by (2):
T si ′ ′ = C 2 λ i ln [ e C 2 / λ i T - 1 ϵ λ i + 1 ] - - - ( 8 ) ,
On Tsi, add " the bright temperature calculated for the second time of expression.
Step 403 is according to (8) ' Ts that calculates of formula 1", Ts 2", Ts 3" value and the actual Ts that measures 1, Ts 2, Ts 3, (9) formula of pressing is obtained the mean square deviation sum and is done
(∑ΔT 2)’=(Ts 1-Ts 1”) 2+(Ts 2-Ts 2”) 2+(Ts 3-Ts 3”) 2 (9)’
If (∑ Δ T 2) '<(∑ Δ T 2), then the approximate very warm T ' of explanation=T+1 than T more near the very warm T of reality VeryThen can make approximate very warm T fully "=T '+1=T+2, continue by (7) ', (8) ', (9) ' formula goes down, until the mean square deviation sum hour, promptly think to seek out the very warm T of reality VeryValue.
If (∑ Δ T 2) '>(∑ Δ T 2), then the approximate very warm T ' of explanation=T+1 is away from the very warm T of reality Very, should seek very temperature value of reality in the other direction.', (8) ', (9) even T '=T-1 continues by (7) ' formula calculates ∑ Δ T 2, until the mean square deviation sum hour, promptly find very temperature value T of reality VeryAbove interative computation software flow calcspar is as shown in Figure 2.
In addition, the radiation temperature measurement instrument of the measurement target true temperature that the present invention proposes is used to carry out above-mentioned temp measuring method, and this temperature measuring device is the three-wavelength temperature measuring device, and it comprises: the temperature acquiring unit is used to obtain the corresponding bright temperature T of actual measurement of three wavelength SiAnd the approximate very warm T of target, wherein the value of i is 1,2 or 3, three wavelength for uniformly-spaced; The supercomputing unit, it further comprises: first module is used for according to the bright temperature Ts of actual measurement that obtains 1, Ts 2, Ts 3And the numerical value of approximate very warm T, confirm the emissivity of target under three wavelength λ 1, ε λ 2, ε λ 3Unit second according to the emissivity value of target under three wavelength and the approximate very warm T value that first module calculates, calculates the bright temperature Ts of corresponding three wavelength of target 1', Ts 2', Ts 3'; Mean square deviation sum ∑ Δ T according to surveying the bright temperature of calculating that try to achieve bright temperature and Unit second, is tried to achieve in Unit the 3rd 2As follows: ∑ Δ T 2=(Ts 1-Ts 1') 2+ (Ts 2-Ts 2') 2+ (Ts 3-Ts 3') 2Unit the 4th, the approximate very warm T ' of its order=T+1, and according to the bright temperature Ts of actual measurement 1, Ts 2, Ts 3Calculate the emissivity of three wavelength of target according to the mode of first module λ 1', ε λ 2', ε λ 3';
Unit the 5th, the emissivity that calculates according to Unit the 4th λ 1', ε λ 2', ε λ 3' reach the value of T '=T+1, calculate three brightness temperature value Ts of target according to the mode of Unit second 1", Ts 2", Ts 3"; Unit the 6th is according to calculating bright temperature Ts 1", Ts 2", Ts 3" value and the bright temperature Ts of actual measurement 1, Ts 2, Ts 3, try to achieve the mean square deviation sum and be: (∑ Δ T 2) '=(Ts 1-Ts 1") 2+ (Ts 2-Ts 2") 2+ (Ts 3-Ts 3") 2If Unit the 7th is (∑ Δ T 2) '<(∑ Δ T 2), then approximate very warm T '=T+1 than T more near the very warm T of reality Very, the approximate very warm T of order then "=T '+1=T+2, continue the processing procedure of repetition Unit the four-six, until the mean square deviation sum hour, promptly think and find out the very warm T of reality VeryIf value is (∑ Δ T 2) '>(∑ Δ T 2), then approximate very warm T '=T+1 is away from the very warm T of reality Very, realistic border is very warm in the other direction is worth, even T '=T-1 continues the processing procedure of repetition Unit the four-six and calculates ∑ Δ T 2, until the mean square deviation sum hour, promptly find very temperature value T of reality Very
Above interative computation adopts modern computing technique fast, can accomplish with interior at 100 milliseconds, can think to measure in real time true temperature.
Above-described specific embodiment; The object of the invention, technical scheme and beneficial effect have been carried out further explain, and institute it should be understood that the above is merely specific embodiment of the present invention; Be not limited to the present invention; All within spirit of the present invention and principle, any modification of being made, be equal to replacement, improvement etc., all should be included within protection scope of the present invention.

Claims (10)

1. the radiation temperature measurement method of a measurement target true temperature, this method comprises:
Step 1 is obtained the corresponding bright temperature T of actual measurement of three wavelength through the three-wavelength Radiation Temperature Measurement Instrument SiAnd the approximate very warm T of target, wherein the value of i is 1,2 or 3, three wavelength for uniformly-spaced;
Step 2 is according to the bright temperature Ts of actual measurement 1, Tss 2, Ts 3And the numerical value of approximate very warm T, confirm the emissivity of target under three wavelength λ 1, ε λ 2, ε λ 3
Step 3 according to the emissivity value of target under three wavelength and the approximate very warm T value that step 2 calculates, is calculated the bright temperature Ts of corresponding three wavelength of target 1', Ts 2', Ts 3';
Step 4 according to the bright temperature of calculating in bright temperature of the actual measurement in the step 1 and the step 3, is tried to achieve mean square deviation sum ∑ Δ T 2As follows: ∑ Δ T 2=(Ts 1-Ts 1') 2+ (Ts 2-Ts 2') 2+ (Ts 3-Ts 3') 2
Step 5, the approximate very warm T ' of order=T+1, and according to the bright temperature Ts of actual measurement 1, Ts 2, Ts 3Calculate the emissivity of three wavelength of target according to the mode of step 2 λ 1', ε λ 2', ε λ 3';
Step 6, the emissivity that calculates according to step 5 λ 1', ε λ 2', ε λ 3' reach the value of T '=T+1, calculate three brightness temperature value Ts of target according to the mode of step 3 1", Ts 2", Ts 3";
Step 7 is according to calculating bright temperature Ts 1", Ts 2", Ts 3" value and the bright temperature Ts of actual measurement 1, Ts 2, Ts 3, try to achieve the mean square deviation sum and be: (∑ Δ T 2) '=(Ts 1-Ts 1") 2+ (Ts 2-Ts 2") 2+ (Ts 3-Ts 3") 2
If step 8 is (∑ Δ T 2) '<(∑ Δ T 2), then approximate very warm T '=T+1 than T more near the very warm T of reality Very, the approximate very warm T of order then "=T '+1=T+2, continue repeating step 5-7, until the mean square deviation sum hour, promptly think and find out the very warm T of reality VeryValue,
If (∑ Δ T 2) '>(∑ Δ T 2), then approximate very warm T '=T+1 is away from the very warm T of reality Very, realistic border is very warm in the other direction is worth, even T '=T-1, continuing set by step, 5-7 calculates ∑ Δ T 2, until the mean square deviation sum hour, promptly find very temperature value T of reality Very
2. method according to claim 1 is characterized in that, in the near-infrared region, the interval between three wavelength is at 0.1~0.3 μ m.
3. method according to claim 1 is characterized in that, in step 2, and emissivity λ 1, ε λ 2, ε λ 3Try to achieve through following formula:
Figure FDA00001949229800021
Wherein, ε λ iBe respectively wavelength X 1, λ 2, λ 3Corresponding target emissivity λ 1, ε λ 2, ε λ 3I λ iBe respectively wavelength X 1, λ 2, λ 3The photocurrent I that corresponding target optical spectrum radiation intensity is produced λ 1, I λ 2, I λ 3I (λ i.T) be respectively wavelength X 1, λ 2, λ 3Article three, calibration curve is similar to the corresponding photocurrent of very warm T.
4. method according to claim 3 is characterized in that, in step 3, calculates the bright temperature Ts of corresponding three wavelength of target through following formula 1', Ts 2', Ts 3':
Figure FDA00001949229800022
Wherein, T is approximate very temperature, C 2Be Planck second radiation constant, λ iBe wavelength, ε λ iBe wavelength X iThe time the target emissivity.
5. the radiation temperature measurement instrument of a measurement target true temperature, this temperature measuring device is the three-wavelength temperature measuring device, it comprises:
The temperature acquiring unit is used to obtain the corresponding bright temperature T of actual measurement of three wavelength SiAnd the approximate very warm T of target, wherein the value of i is 1,2 or 3, three wavelength for uniformly-spaced;
The supercomputing unit, it further comprises:
First module is used for according to the bright temperature Ts of actual measurement that obtains 1, Ts 2, Ts 3And the numerical value of approximate very warm T, confirm the emissivity of target under three wavelength λ 1, ε λ 2, ε λ 3
Unit second according to the emissivity value of target under three wavelength and the approximate very warm T value that first module calculates, calculates the bright temperature Ts of corresponding three wavelength of target 1', Ts 2', Ts 3';
Mean square deviation sum ∑ Δ T according to surveying the bright temperature of calculating that try to achieve bright temperature and Unit second, is tried to achieve in Unit the 3rd 2As follows: ∑ Δ T 2=(Ts 1-Ts 1') 2+ (Ts 2-Ts 2') 2+ (Ts 3-Ts 3') 2
Unit the 4th, the approximate very warm T ' of its order=T+1, and according to the bright temperature Ts of actual measurement 1, Ts 2, Ts 3Calculate the emissivity of three wavelength of target according to the mode of first module λ 1', ε λ 2', ε λ 3';
Unit the 5th, the emissivity that calculates according to Unit the 4th λ 1', ε λ 2', ε λ 3' reach the value of T '=T+1, calculate three brightness temperature value Ts of target according to the mode of Unit second 1", Ts 2", Ts 3";
Unit the 6th is according to calculating bright temperature Ts 1", Ts 2", Ts 3" value and the bright temperature Ts of actual measurement 1, Ts 2, Ts 3, try to achieve the mean square deviation sum and be: (∑ Δ T 2) '=(T 1-Ts 1") 2+ (Ts 2-Ts 2") 2+ (Ts 3-Ts 3") 2
If Unit the 7th is (∑ Δ T 2) '<(∑ Δ T 2), then approximate very warm T '=T+1 than T more near the very warm T of reality Very, the approximate very warm T of order then "=T '+1=T+2, continue the processing procedure of repetition Unit the four-six, until the mean square deviation sum hour, promptly think and find out the very warm T of reality VeryIf value is (∑ Δ T 2) '>(∑ Δ T 2), then approximate very warm T '=T+1 is away from the very warm T of reality Very, realistic border is very warm in the other direction is worth, even T '=T-1 continues the processing procedure of repetition Unit the four-six and calculates ∑ Δ T 2, until the mean square deviation sum hour, promptly find very temperature value T of reality Very
6. temperature measuring device according to claim 5 is characterized in that, in the near-infrared region, the interval between three wavelength is at 0.1~0.3 μ m.
7. temperature measuring device according to claim 6 is characterized in that,
In first module, emissivity λ 1, ε λ 2, ε λ 3Try to achieve through following formula:
Figure FDA00001949229800031
Wherein, ε λ iBe respectively wavelength X 1, λ 2, λ 3Corresponding target emissivity λ 1, ε λ 2, ε λ 3I λ iBe respectively wavelength X 1, λ 2, λ 3The photocurrent I that corresponding target optical spectrum radiation intensity is produced λ 1, I λ 2, I λ 3I (λ i.T) be respectively wavelength X 1, λ 2, λ 3Article three, calibration curve is similar to the corresponding photocurrent of very warm T.
8. temperature measuring device according to claim 7 is characterized in that, in step 3, calculates the bright temperature Ts of corresponding three wavelength of target through following formula 1', Ts 2', Ts 3':
Figure FDA00001949229800041
Wherein, T is approximate very temperature, C 2Be Planck second radiation constant, λ iBe wavelength, ε λ iBe wavelength X iThe time the target emissivity.
9. according to each described temperature measuring device of claim 5-8; It is characterized in that chopper wheel, measuring light fulgurite, LED, control that this instrument further comprises the catoptron, aiming catoptron, target imaging lenslet of object lens, band central small hole, driven by motor M are used eyepiece with synchronizable optical fulgurite, observation.
10. temperature measuring device according to claim 9 is characterized in that, wherein the object lens adjustable focus of temperature measurer; With target imaging in the band central small hole catoptron on; After catoptron reflects again, supply eye-observation aiming usefulness through the eyepiece imaging, behind this catoptron of light transmission of the central small hole of the catoptron of band central small hole through aiming by the light of the mirror reflects of band central small hole; Behind the target imaging lenslet; Chopper wheel through being driven by modulating motor M is equipped with three narrow band pass filters on the chopper wheel, let λ respectively 1, λ 2, λ 3Chopper wheel is passed through in the light timesharing of three wavelength; Be radiated on the measuring light fulgurite; An aperture is arranged on the chopper wheel, and when the light transmission aperture of the LED that is installed in chopper wheel left side, short irradiation is on synchronously with photoelectric tube;, the chopper wheel revolution produces a synchronization pulse when moving a week; Distinguish three signals of measuring wavelength when synchronization pulse is used for signal Processing, the signal of measuring light fulgurite is input to storer with synchronization pulse jointly after prime amplifier and A/D conversion.
CN201210266196XA 2012-07-30 2012-07-30 Radiation temperature measuring method and apparatus for measuring true target temperature Pending CN102749141A (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201210266196XA CN102749141A (en) 2012-07-30 2012-07-30 Radiation temperature measuring method and apparatus for measuring true target temperature
CN201210417758.6A CN103048050B (en) 2012-07-30 2012-10-26 Radiation temperature measurement method and instrument for measuring real temperature of target

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210266196XA CN102749141A (en) 2012-07-30 2012-07-30 Radiation temperature measuring method and apparatus for measuring true target temperature

Publications (1)

Publication Number Publication Date
CN102749141A true CN102749141A (en) 2012-10-24

Family

ID=47029514

Family Applications (2)

Application Number Title Priority Date Filing Date
CN201210266196XA Pending CN102749141A (en) 2012-07-30 2012-07-30 Radiation temperature measuring method and apparatus for measuring true target temperature
CN201210417758.6A Active CN103048050B (en) 2012-07-30 2012-10-26 Radiation temperature measurement method and instrument for measuring real temperature of target

Family Applications After (1)

Application Number Title Priority Date Filing Date
CN201210417758.6A Active CN103048050B (en) 2012-07-30 2012-10-26 Radiation temperature measurement method and instrument for measuring real temperature of target

Country Status (1)

Country Link
CN (2) CN102749141A (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105333962A (en) * 2014-06-18 2016-02-17 中微半导体设备(上海)有限公司 Temperature measurement method and system for correcting dual-wave-band temperature measurement error
CN106066340A (en) * 2016-05-24 2016-11-02 南京航空航天大学 A kind of analogue experiment method measuring solid target infrared signature in normal temperature environment
CN106768366A (en) * 2017-01-12 2017-05-31 杭州电子科技大学 A kind of utilization micro spectrometer measures the method and device of high temperature
CN111044565A (en) * 2019-12-30 2020-04-21 西安近代化学研究所 True temperature multispectral measurement method in process of converting combustion into detonation
CN111562019A (en) * 2020-04-08 2020-08-21 太原理工大学 Multispectral radiation temperature measurement method and system
CN112220460A (en) * 2020-09-25 2021-01-15 电子科技大学 In-vivo temperature measurement patch based on microwave radiation technology
CN113175998A (en) * 2021-03-26 2021-07-27 合肥工业大学 Metal material surface temperature measurement method based on colorimetric temperature measurement
CN114112061A (en) * 2021-10-11 2022-03-01 黑龙江业博科技有限公司 Method and system for rapidly measuring temperature of welding head of transistor type welding power supply on line

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103411940B (en) * 2013-08-08 2014-06-18 哈尔滨工业大学 Detection method and test device for catalytic properties of heatproof material based on emission spectrum
CN105181169B (en) * 2015-05-12 2019-01-15 中国科学院苏州生物医学工程技术研究所 Thermometry, temperature measurement system and temperature acquisition device
CN105509901B (en) * 2015-12-29 2018-11-02 温州大学 High-precision measuring method of the illuminator with respect to colour temperature

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0294747B1 (en) * 1987-06-06 1994-02-16 Minolta Camera Kabushiki Kaisha Pyrometer
GB2344883B (en) * 1998-12-16 2003-10-29 Graviner Ltd Kidde Flame monitoring methods and apparatus
US7034302B2 (en) * 2002-09-19 2006-04-25 Battelle Energy Alliance, Llc Optical steam quality measurement system and method
RU2251669C2 (en) * 2003-07-08 2005-05-10 Открытое акционерное общество "КАМАЗ" Method of measuring color temperature of metal
EP1635163B1 (en) * 2004-09-09 2017-05-31 Volkswagen Aktiengesellschaft Motor vehicle comprising a device for determining a surface condition of a roadway
US20120002035A1 (en) * 2010-06-30 2012-01-05 General Electric Company Multi-spectral system and method for generating multi-dimensional temperature data
CN201867924U (en) * 2010-11-15 2011-06-15 上海翼捷工业安防技术有限公司 Three-wavelength infrared flame detector
CN102426753B (en) * 2011-11-21 2013-07-17 上海大学 Infrared multiparameter fire and gas detector

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105333962A (en) * 2014-06-18 2016-02-17 中微半导体设备(上海)有限公司 Temperature measurement method and system for correcting dual-wave-band temperature measurement error
CN105333962B (en) * 2014-06-18 2018-06-22 中微半导体设备(上海)有限公司 A kind of thermometry and system for correcting two waveband temperature measurement error
CN106066340A (en) * 2016-05-24 2016-11-02 南京航空航天大学 A kind of analogue experiment method measuring solid target infrared signature in normal temperature environment
CN106768366A (en) * 2017-01-12 2017-05-31 杭州电子科技大学 A kind of utilization micro spectrometer measures the method and device of high temperature
CN111044565A (en) * 2019-12-30 2020-04-21 西安近代化学研究所 True temperature multispectral measurement method in process of converting combustion into detonation
CN111562019A (en) * 2020-04-08 2020-08-21 太原理工大学 Multispectral radiation temperature measurement method and system
CN112220460A (en) * 2020-09-25 2021-01-15 电子科技大学 In-vivo temperature measurement patch based on microwave radiation technology
CN113175998A (en) * 2021-03-26 2021-07-27 合肥工业大学 Metal material surface temperature measurement method based on colorimetric temperature measurement
CN113175998B (en) * 2021-03-26 2022-09-09 合肥工业大学 Metal material surface temperature measurement method based on colorimetric temperature measurement
CN114112061A (en) * 2021-10-11 2022-03-01 黑龙江业博科技有限公司 Method and system for rapidly measuring temperature of welding head of transistor type welding power supply on line
CN114112061B (en) * 2021-10-11 2024-02-02 黑龙江业博科技有限公司 Quick online measuring method and system for temperature of welding head of transistor welding power supply

Also Published As

Publication number Publication date
CN103048050A (en) 2013-04-17
CN103048050B (en) 2014-09-24

Similar Documents

Publication Publication Date Title
CN103048050B (en) Radiation temperature measurement method and instrument for measuring real temperature of target
CN103076107B (en) Terahertz pulse measurement-based burning temperature sensing device and method
CN103063312B (en) Measuring system and method for measuring object emissivity
CN112378616B (en) High-speed flow field multi-parameter measurement system and method based on wavelength modulation absorption spectrum
CN203216634U (en) High temperature measurement and temperature field reconstruction device based on turnable diode laser absorption spectroscopy
CN102889934A (en) Method for measuring temperature in real time
CN107490397B (en) High-accuracy self-adaptation filters the quick Peak Search Method of FBG spectrum
CN108981953B (en) Laser absorption spectrum temperature measurement method and system based on interferometric modulation principle
CN108562237B (en) Device and method for performing spectrum calibration in optical frequency domain reflection sensing system by adopting HCN (hydrogen cyanide) air chamber
CN102062636A (en) Portable site laser energy measuring device
CN104864979A (en) Correction method of errors measured by distributed raman optical fiber temperature measuring system
CN101000264B (en) Method for spectral temperture measuring by Wien''s displacement law
CN106979832A (en) A kind of optical fibre light splitting temp measuring system and its temp measuring method
CN103557960A (en) Fabry-Perot fiber-optic temperature sensing system and method
CN103675019A (en) Method for measuring material surface emissivity by virtue of thermal infrared imager rapidly
CN103743489A (en) Infrared radiometer calibration method on basis of standard plane source black body
CN111947785A (en) Temperature measurement thermal infrared imager noise equivalent temperature difference calibration method
CN1601244B (en) Focus detecting unit, and refractive index measuring apparatus and non-contact thermometer using the same
CN106092320A (en) A kind of spectrum calibration method of LONG WAVE INFRARED spectrogrph
CN109580033A (en) A kind of concrete dam distributed optical fiber temperature measurement data error compensation method
CN102589714A (en) Temperature measuring device based on high-pressure gas Rayleigh-Brillouin scattering spectrum
CN108801496B (en) Path temperature histogram measurement system and method based on overlapping absorption spectrum
CN101915618A (en) Device and method for calibrating emissivity of high-temperature fuel gas
CN105806491A (en) Three-wavelength two-dimensional temperature field measuring device and method
CN201892573U (en) Near-infrared radiation thermometer

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20121024