CC BY
1Ц.Ц. Дамбиев, д-р техн. наук, проф.
Т.В. Еремина, д-р техн. наук А.Ф. Мункуев, аспирант Восточно-Сибирский государственный университет технологий и управлений
УДК 536:574(571.54)
ОРГАНИЗАЦИЯ МОНИТОРИНГА ВОЗДУШНОГО БАССЕЙНА ТЕПЛОВОЙ ЭЛЕКТРИЧЕСКОЙ СТАНЦИИ С ПОМОЩЬЮ ЛАЗЕРНОГО ЗОНДИРОВАНИЯ
Рассмотрен принцип лазерного зондирования вредных выбросов ТЭС.
Ключевые слова: лазерное зондирование, воздушный бассейн, выбросы ТЭС.
С.С Dambiev, D.Sc. Engineering, Prof.
T.V. Yeremina, D.Sc. Engineering, Prof.
A.F. Munkuev
MONITORING OF AIR BASIN THERMAL POWER PLANTS USING LASER SENSING
The principle of laser probing of thermal power plant emissions.
Key words: leser sensing, air basin, emissions of TPP.
Сам по себе принцип лазерного зондирования аэрозолей в воздушном бассейне ТЭС чрезвычайно прост. В атмосферу направляется световой импульс малой длительности т. Он высвечивает область длиной ст, где с - скорость света. Остальное пространство является не освещенным. Каждая освещенная пылинка рассеивает во все стороны попавшее на нее лазерное излучение, часть которого попадает на приемный телескоп и детектируется. Чем больше концентрация пылинок, тем больше зарегистрированный сигнал. Фактически ситуация сильно осложняется дифракционными процессами, так как размеры частиц соизмеримы с длиной волны лазера. Такие процессы подробно изложены в целом ряде работ [1, 2].
Для практического использования в поставленной задаче вполне применимо уравнение лазерного зондирования (УЛЗ) в приближении однократного рассеяния света на аэрозолях. Задача облегчается тем, что для примененного лазера на 2-й гармонике иттрий - алюминиевого граната (зеленый свет) практически полностью отсутствует поглощение света в атмосферных газах и аэрозолей. То есть все ослабление света зависит только от светорассеивающих свойств зондируемой среды. При таких условиях УЛЗ выглядит следующим образом:
р е x
P(x) = 0 'h • • s(x)exp[-2ís(Z)dZ] (1)
2 •x 0
Здесь P(x) - мощность эхо-сигнала, пришедшего на детектор лидара из точки, находящейся на расстоянии х от лидара;
Р0- излученная лазером импульсная мощность;
S0 - площадь приемного телескопа;
П - кпд всей оптической системы;
о(х) - распределенное вдоль трассы зондирования значение объемного коэффициента рассеяния излучения;
уп - значение индикатрисы рассеяния.
Экспоненциальный сомножитель представляет собой пропускание атмосферы по трассе от 0 до х и назад от нее до точки 0, т.е. до лидара. Если атмосфера не слишком мутная (т.е. это еще не облако и не дымовой шлейф), то этот сомножитель близок к единице.
Теоретически и экспериментально получено, что связь между массовой концентрацией гомогенного аэрозоля и коэффициентом рассеяния излучения т имеет характер линейной регрессии:
M [ мг / м3] = m[ мг • км / м3] •а[км(2)
Коэффициент пропорциональности ц сложным образом зависит от вида функции распределения частиц по размерам. Дело в том, что величина о связана со среднеквадратичным размером частиц, а М -со среднекубическим. Если заранее известен спектр размеров частиц, то коэффициент ц может быть рассчитан достаточно точно. Так, в [1] приведено измеренное в атмосфере промышленного центра (но не в
дымах труб) ц=(0,52 ± 0,34) [мг-км/м3]. Высокое значение стандартного отклонения (58%) обусловлено широким набором аэрозолей различного происхождения. Для аэрозоля с фиксированными свойствами, например дыма от источника, питаемого углем постоянного состава, этот разброс значительно уменьшается, и достоверность измерений массовой концентрации повышается.
Таким образом, в УЛЗ полезную информацию представляет объемный коэффициент рассеяния о. Он содержится в уравнении как непосредственно, так и через множитель экспоненциального ослабления сигнала. На малых расстояниях до дымового шлейфа (сотни метров при отсутствии тумана) этот множитель слабо отличается от единицы. И если ближняя к лидару граница шлейфа находится на расстоянии х1, то коэффициент рассеяния дыма будет
* = *(*.) = 2 14 ^--(3)
P0 •с•т • S0
Алгоритм прост, но требует калибровки лидара, т.е. знания величин Р0, т, 80, П- Необходимо также априорное знание индикатрисы обратного рассеяния уп. Это его уязвимые позиции.
Специально для тонких дымовых слоев может быть предложен еще один несложный алгоритм. Пренебрегая потерями излучения в атмосфере до шлейфа, запишем мощность сигнала от его ближней к лидару границы:
Р(X,) = Ро "С"Т " 50 П -ГрМX,) (4)
2 • X:
Мощность сигнала от дальней границы облака при тех же условиях
Р • с• т • 5 • л хс
Р(х2) = 0 20 1 -уя •о(х2)"ехр-2" ¡«Я" (5)
2 • Х„
2 X,
2 » 2
Поскольку толщина слоя дыма невелика, то X, » Х2 , и среднее значение коэффициента рассеяния по шлейфу между точками х2 и х:
- 1 , Р(X])
- =--1п-— (6)
2(x2 - X,) Р(x2)
В тех ситуациях, когда слой дыма достаточно толст по пространству и плотен и выполняется условие , применим асимптотический алгоритм определения пространственного распре-
деления о(х):
-( X) = --(7)
2 J P(0 Z2 dZ
Для интерпретации результатов лазерного зондирования полезным является использование поляризационных эффектов. Из теории рассеяния света следует, что ансамбль частиц, состоящих из идеальных сфер, при рассеянии строго назад должен сохранять состояние поляризации, присущее зондирующему пучку. Например, если лазерное излучение линейно поляризовано в какой-то плоскости, то и однократно рассеянное в направлении 180° оно поляризовано в этой же плоскости. Поэтому наблюдающаяся в экспериментах деполяризация однократно рассеянного излучения обусловлена отклонением формы частиц от сферической. Так будет при зондировании дымов вблизи.
На более значительных дистанциях на приемник начнет попадать излучение, фотоны которого многократно испытали рассеяние на частицах дыма, капельках и тому подобное. В этом случае величина деполяризации эхосигнала будет возрастать по мере углубления в облака. Для светлого дыма, тумана над градирней деполяризация будет возрастать достаточно быстро. Для темного дыма этот процесс будет происходить медленнее, так как в темном веществе поглощается больше фотонов и снижается кратность излучения, отраженного в приемный телескоп.
Рис.1. Принципиальная схема лидара
Принципиальная схема лидара приведена на рисунке 1. Здесь 1 - лазерный излучатель типа ИЛТИ - 405, работающий на кристалле иттрий - алюминиевого граната; 2 - фотометрическая сфера для усреднения излучения; линзы 3 и 4 образуют телескоп коллиматора; линзы 5 и 6 - приемный телескоп; 6 -сменная диафрагма; 8 - интерференционный светофильтр; 9 - поворотный поляроидный фильтр; 10 -стеклянное кольцо для подачи излучения от световода 18 на катод фотоумножителя (ФЭУ) 11; 12 - блок высоковольтного питания фотоумножителя; 13 - аналогово-цифровой преобразователь (АЦП) с телевизионным дисплеем 14; 15 - интерфейс связи с ЭВМ 16; 17 - блок обеспечения работы лазера.
Проходящий в атмосфере вдоль оптической оси лазера импульс света частично рассеивается на встречающихся на его пути частицах дыма, капельках воды и иных аэрозолях. Чем выше счетная концентрация аэрозолей и их размер (т. е. их массовая концентрация), тем больше мощность излучения, собираемая приемным телескопом 5-7.
АЦП преобразует непрерывный сигнал с ФЭУ 11, промодулированный атмосферными неодно-родностями, в дискретный, цифровой вид, удобный для обработки на ЭВМ.
На телемониторе в реальном масштабе времени и в виде, близком к обычному телевизионному стандарту, представлены 256 эхосигналов лидара. При каждой последующей вспышке лазера один новый эхосигнал появляется на экране, а последний на данный момент исчезает с него. При этом можно вести зондирование дымов над срезами труб ТЭС как по вертикали, так и вдоль горизонтальной трассы в предположении, что атмосфера близка к однородной.
Калибровка лидара ставит в соответствие мощность принятого сигнала Р(х) мощности излученного импульса Р0. Последняя составляет миллионы ватт в импульсе и, естественно, не может быть прямо определена фотоприемником. Наиболее приемлема калибровка по свободной атмосфере, когда оценивается оператором лидара метеорологическая дальность видимости 8М. Она связана с коэффициентом рас-
3.9
сеяния соотношением о =- , когда
„
°м
Р(х) ■ х2 (8)
о(х) = -
%■ То2 - 2|Р(С) ■С'йС
Здесь х - расстояние, от которого начинается регистрация принимаемого сигнала; Т02 = ехр[-2ох] -пропускание атмосферы до этой дистанции; % - калибровочная постоянная, мультипликативно объединяющая все входящие в (1) константы, включая индикатрису атмосферы уп. Если проводить зондирование по горизонтали в условиях, когда можно положить а(х)=сои81;, записать решение (8) для двух расстояний х1и х2 и приравнять их друг к другу, то получим
С=—--X-. (9)
Т2 Нх)Х -ЯхгХ
При выполнении такого рода калибровки по результатам зондирования атмосферы с крыши административного здания Улан-Удэнской ТЭЦ-1 в восточном направлении было получено, что с =2,1-109. Это значит, что величина 0,5-Р0-с-т-80-П7л=2.1-109 единиц кода аналогово-цифрового преобразователя
(АЦП) (в среднем). Мощность лазера эквивалентна 7,47-1010 единиц кода АЦП.
В качестве эксперимента было проведено зондирование труб Улан-Удэнской ТЭЦ-1. Расстояние по трассе с наклоном 240 до границы дымового выброса равнялась 630 м. Луч лазера проходил в шлейфе на высоте 20-25 м над устьем трубы, высота которой составляет 180 м. Лидар располагался на высоте
около 25 м над уровнем земли.
Рис.2. Результаты зондирования шлейфа
На рисунке 2 представлены результаты зондирования шлейфа, который визуально и по видеозаписи стабильно представлял собой довольно прозрачное образование.
Библиография
1. Дамбиев Ц.Ц. Термодинамическая модель устойчивого развития региона: экологический и энергосберегающий аспекты. - Улан-Удэ: ВСГТУ, 2001. - 400 с.
2. Дамбиев Ц.Ц., Шаманаев В.С., Пеннер И.Э. и др. О влияние режимов работы ТЭЦ на характеристики ли-дарных сигналов // В кн. Материалы XI международного симпозиума по лазерному и акустическому зондированию атмосферы. - Томск: Институт оптики атмосферы СО РАН, 1992. - С. 12-16.
Bibliography
1. Dambiev Ts. Ts. Termodinamic model sleady evolution region: ecological and powereconomy aspect. - Ulan-Ude the East Siberian Sea state technological university, 2001. - 400 p.
2. Dambiev Ts. Ts.,Shamanaev V.S., Penner 1.Э., et al. Of influence regime work TES on description leader signal. In book materials Х1 international symposium the laser and acoustical probe atmosphere. - Tomsc: institute optics atmosphere СОРАН, 1992. - P. 12-16.
