способ пирометрических измерений

Формула изобретения

1. Способ пирометрических измерений, включающий измерение истинной температуры поверхности по величине яркостных температур не менее чем на трех длинах волн, отличающийся тем, что, измеряя не менее трех компонент спектра, вычисляют логарифм произведения коэффициента спектральной излучательной способности на коэффициент пропускания промежуточной среды для каждой длины волны, при этом каждые три значения длин волн выбираются настолько близкими, что зависимость логарифма названного произведения от длины волны будет линейна, вычитают значения логарифма названного произведения из измеряемых обратных значений яркостных температур и полученные выражения суммируют для получения обратного значения искомой истинной температуры.

2. Способ по п.1, отличающийся тем, что увеличивают количество длин волн и спектральных интервалов при сохранении их ширины.

3. Способ по п.1, отличающийся тем, что увеличивают количество длин волн и спектральных интервалов с сужением величины последних.

Описание изобретения к патенту

Изобретение относится к измерительной технике, а более конкретно к оптическим бесконтактным способам измерения излучательной способности и истинных температур различных объектов.

В методах оптической пирометрии результат измерений зависит не только от температуры, но и от свойств поверхности излучающего тела, пропускания промежуточной среды. Априорная информация об этих параметрах при измерении истинной температуры обычно используется в виде коэффициентов излучательной способности и пропускания промежуточной среды, табулированных в справочной литературе. Неполнота данных по этим коэффициентам, несоответствие их значений в реальных условиях табулированным, изменение свойств излучающей поверхности и промежуточной среды в процессе измерения приводят к возникновению методических погрешностей, превышающих инструментальные в десятки раз. С другой стороны, без определенной априорной информации об излучательных характеристиках объекта контроля сама задача об измерении его температуры методами оптической пирометрии оказывается беспочвенной.

Из предыдущего уровня техники известны различные варианты способов измерения излучательной способности и истинной температуры, основанные на многоволновом принципе.

Смотри, например, патент Д.Я.Света RU 2162210 C1, 7 G 01 J 5/50, недостаток описанного в нем способа заключается в невозможности его применения для полупрозрачных тел и необходимости в дополнительном двухволновом источнике излучения.

Известны также многоволновые способы спектрального отношения измерения истинных температур, например, описанные Hombeck G.A. Temperature: Its. Measurement and Control in Science and Industry, N.Y.Reinolds, 1962; Bourianes R. Morea M., Review Phys. Appl. 1977 v.42, p.893 и др., в частности, более поздний патент США 577323 Int. cl. G 01 J 5/00, G 01 J 5/08.

Все варианты способа спектрального отношения включают исключение “влияния” излучательной способности на результат измерения температуры при корректности, принятой аппроксимации зависимости спектральной излучательной способности от длины волны.

Недостатком всех способов отношения является весьма большая погрешность, за счет увеличения значения эквивалентной длины волны.

Это объясняется тем, что при расширении спектрального интервала принятая аппроксимация, например линейная, становится некорректной - теряет смысл. При сужении же спектрального интервала в методах/отношения возрастает эквивалентная длина волны, это приводит к уменьшению полезного сигнала и увеличивается погрешность измерения истинной температуры от влияния излучательной способности и пропускания промежуточной среды.

Известно, что величина логарифма сигнала InU(,T) в оптической пирометрии согласно закона Вина обратно пропорциональна эквивалентной длине волны, а погрешность измеряемой температуры пропорциональна величине последней.

Так, основная формула, определяющая погрешность в оптической пирометрии, гласит:

Т/Т= / · _эТ/С₂

Здесь: Т/Т - погрешность измерения температуры,

- коэффициент излучательной способности, ее произведение на коэффициент пропускания промежуточной среды,

/ - погрешность в значении излучательной способности (пропускания промежуточной среды),

Т - величина измеряемой температуры,

С₂ - 14388 мкм град. - 2-я пирометрическая константа,

_экв - эквивалентная длина волны пирометрических систем.

Так, например, в способе двойного спектрального отношения при ₁=0,72 мкм, ₂=0,74 мкм и ₃=0,76 мкм

^* _экв= _{1 2 3}/(2 _{1 3}- _{1 2}- _{2 3})=506,6 мкм

В способе одинарного спектрального отношения (двухволновой пирометрии) величина эквивалентной длины волны также велика:

^**= _{1 2}/( ₁- ₂)=0,72· 0,74/0,02=26,64

За прототип нами взят способ яркостной взаимно корреляционной пирометрии, описанной нами в авторском свидетельстве №368500 (опубликовано в Б.И. №9 за 1973 год).

Во взаимно корреляционной пирометрии температура определяется не по величине отношения, а по величине произведения спектральных интенсивностей излучения.

Так для произведения двух спектральных интенсивностей излучения: U( ₁,T) и U( ₂,T) величина логарифма их произведения будет:

lnU( ₁,T)U( ₂,T)=In _{1 2}-C₂/ _эТ

и величина эквивалентной длины волны:

_экв= _{1 2}/( ₁+ ₂)

При весьма близких значениях ₁ и ₂, _э /2, т.е. при

где _э /n

Таким образом, в методе взаимно корреляционной пирометрии эквивалентная длина волны не увеличивается с увеличением числа спектральных интенсивностей, как в методах отношения, а уменьшается, т.е. величина полезного сигнала увеличивается.

Недостаток способа взаимно корреляционной пирометрии является измерение не истинной, а взаимно корреляционной температуры.

Поставленная задача измерения истинной температуры поверхности по величине яркостных температур не менее чем на трех длинах волн решена тем, что, измеряя не менее трех компонент спектра, вычисляют логарифм произведения коэффициента спектральной излучательной способности на коэффициент пропускания промежуточной среды для каждой длины волны, при этом каждые три значения длины волны выбираются настолько близкими, что зависимость логарифма названного произведения от длины волны будет линейна; вычитают значения названного произведения из измеренных обратных значений яркостных температур, и полученные выражения суммируют для получения обратного значения искомой истинной температуры.

Таким образом, эквивалентная длина волны измеренной истинной взаимно корреляционной температуры уменьшается минимум в 3 раза. За счет чего величина полезного пирометрического сигнала увеличивается также в 3 раза.

В случае если требуемая величина погрешности не достигается из-за недостаточной величины сигнал/шум, следует увеличивать количество длин волн, сохраняя расстояние между ними, т.е. сохраняя линейность зависимости In ()=а ₀+а₁ на интервалах:

In ( )=a_0I+a_1,I на интервале ₁- ₃

In ( )=а_0II+a_1,II на интервале ₂- ₄

... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ...

... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ...

In ( )=a_0N+a_1,N на интервале _N-1- _N

Если же причиной погрешности является немонотонность спектральной зависимости произведения спектральной излучательной способности на пропускание промежуточной среды, следует сближать значения используемых длин волн, т.е. уменьшать расстояние между ними.

Рассмотрим применение предлагаемого способа при достаточно монотонном характере зависимости от длины волны произведения спектральной излучательной способности и пропускания промежуточной среды, имеющем, в частности, место в металлургии и других областях.

В этом случае каждые три значения длин волн выбираются настолько близкими, что зависимость логарифма ( ) от длины волны будет линейна, т.е.:

In ( ) =a_0, +a_1,

где: a_0, , а_1, - коэффициенты полинома первой степени, аппроксимирующего указанную зависимость и

( ) - произведение спектральной излучательной способности на пропускание промежуточной среды, на интервале ₁- ₃.

Обратные значения измеряемых трех яркостных температур запишем в обычном виде:

Для простоты дальнейших вычислений положим, что длины волн выбраны эквидистантно, т.е.

₁= ; ₂= + ; ₃= +2 , и введем следующие обозначения:

Тогда:

Следовательно, величина коэффициента спектральной излучательной способности (его произведения на коэффициент пропускания промежуточной среды) на длине волны будет определена как:

Вычитая значение In ( ) из соответствующего обратного значения измеренной яркостной температуры (1), (2) и (3), определим обратную величину истинной температуры на спектральном интервале ₁- ₂:

+[I_1,2/ -I_2,3/2 ( ₁+ ₂+ )]C₂=Т_n ^-1

Определяя, таким образом, Т_n из ₁, ₂ и ₃ и суммируя результаты, мы уменьшаем в три раза погрешность измерения истинной температуры Т_n за счет укорочения эквивалентной длины волны.

Таким образом, в предлагаемом способе на спектральном интервале ( ₁- ₃):

т.е. _э ^-1=1/ ₁+1/ ₂+1/ ₃

и при << , _э 1/3 ₁.

Очевидно, что тоже самое мы получим, используя вместо яркостных температур спектральные интенсивности: U( ,T), т.к.

lnU( ,T)=In +In ( ,T)-C₂/ T_n

где: - аппаратная постоянная, определяемая по градуировке.

При высоком уровне шумов, вызывающем недопустимую величину погрешности измерения истинной температуры в предлагаемом способе, следует увеличивать количество длин волн при сохранении выбранной ширины спектрального интервала ( ).

При существенной немонотонности названой зависимости количество используемых длин волн также следует увеличить, но при этом сужать величину интервала .

Рассмотрим примеры применения предлагаемого способа.

Пример 1

Определение истинной температуры вольфрама.

Вольфрам является одним из наиболее селективных металлов. Выбираем 3 длины волны достаточно близкими: ₁=0,9; ₂=0,92 и ₃=0,94.

Очевидно, что на таком узком спектральном интервале - 0,04 мкм линейная зависимость In () _w от длины волны будет вполне корректна.

С помощью таблицы из справочника “Излучательные свойства металлов” под ред. А.Е.Шейндлина, Энергия, 1974, для вольфрама разность обратных значений яркостных температур - формулы (2) и (3) при значении истинной температуры 1600К соответственно будут:

I _1,2=2,337· 10^-6 I_3,4=2,3207· 10^-6

Из формулы (4) и (5) определяем коэффициенты полинома а₀ и a₁. Затем по формуле (6) находим значение логарифма: In (0,9,Т_и)=0,951. Потенцируя, находим (0,9)=0,386. Табличное значение (0,9)_табл=0,399, т.е. погрешность в определении излучательной способности < / <3,5%, и погрешность в истинной температуре на одной длине волны будет:

Т/Т= / · ₁Т_и/С₂=0,035· 0,9· 1600/14388 0,38%

За счет определения ( ₂,T) и ( ₃,Т) и обратных значений: ( ₂Т_и)^-1 и ( ₃Т_и)^-1 на длинах волн ₂=0,92 мкм и ₃=0,94 мкм можно еще уменьшить значение эквивалентной длины волны, т.е. уменьшить погрешность определения истинной температуры еще в 3(!) раза, т.е. составит Т/Т 0,12%.

Пример 2

Элиминирование влияния дыма при измерении температуры в черной металлургии.

Требуется измерить температуру в условиях достаточно сильного задымления.

Положим, что пропускание дыма составляет ()=50%.

Расчет параметров дыма ведется, как обычно, по формуле Ангстрема (см., например, Ван де Хюлст “Рассеяние света малыми частицами”, М.: Иностранная литература, 1961 г.)

( )=ехр-(SI/ )

Здесь:

SI - произведение “дисперсности” дыма на “толщу” дыма;

=1,3 - коэффициент для “бурого” (Принятый в черной металлургии официальный термин) дыма.

(см. Velon R.G.Joul. Appi, Phys. 1962, #33, р.2140)

Так как зависимость пропускания дыма от длины волны явно нелинейная (экспоненциальная, по формуле Ангстрема), то для корректности линейной аппроксимации зависимости логарифма ( ₁) от длины волны, то спектральный интервал выберем более узким (чем в 1-ом примере).

Например: ₁=0,9; ₂=0,905 и ₃=0,91, т.е. =0,005 мкм.

Определив для ()=50% на волне ₁=0,9 мкм величину SI

SI= In ()=0,9 ^1,3ln0,5=0,793

Элементарно найдем логарифмы пропускания дыма для трех длин волн:

In ^*(0,9)=-0,6914

In ^*(0,905)=-0,698

In ^*(0,91)=-0,7019

Определяя предлагаемым способом, коэффициент пропускания дыма на длине волны ₁=0,9 получим:

In (0,9)=-0,686

Потенцируя, имеем (0,9)=0,5036

И погрешность в определении коэффициента пропускания дыма будет:

[ ^*(0,9)- (0,9)]/ ^*(0,9)=0,0072=0,72%

При Т=2500K это даст погрешность в температуре:

Т/Т=0,0072· 2500/14388· 0,9=0,12%

Используя еще две длины волны ₂=0,905 и ₃=0,91, можно получить погрешность ~0,005%.

Таким образом, в данном случае спектральный интервал =0,05 можно увеличить.

Принципиальная схема практической реализации предлагаемого способа изображена на чертеже.

Здесь:

1 - концентрирующая и визирная оптика;

2 - устройство для выделения необходимых спектральных составляющих;

3 - приемники излучения в виде отдельных элементов или линейки;

4 - усилители;

5 - аналого-цифровые преобразователи;

6 - микропроцессор или компьютер.