Научная статья на тему 'Прозрачномеры-газоанализаторы на двухволновом полупроводниковом лазере'

Прозрачномеры-газоанализаторы на двухволновом полупроводниковом лазере Текст научной статьи по специальности «Физика»

CC BY
129
42
i Надоели баннеры? Вы всегда можете отключить рекламу.
Область наук
Ключевые слова
ИЗМЕРИТЕЛЬ ОПТИЧЕСКИХ ХАРАКТЕРИСТИК / ДВУХВОЛНОВОЙ ПОЛУПРОВОДНИКОВЫЙ ЛАЗЕР / БАЗИСНЫЙ И НЕФЕЛОМЕТРИЧЕСКИЙ МЕТОДЫ

Аннотация научной статьи по физике, автор научной работы — Козлов В. Л., Кугейко М. М.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по физике , автор научной работы — Козлов В. Л., Кугейко М. М.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Transparency meters and gas analyzers based on two-wavelength semiconductor lasers

Construction techniques of precision measuring instruments of optical characteristics on the basis of two-wave lasers with use of basic and nephelometer methods are presented. System stability to changes of hardware constants, influence of an environment, pollution of optics is shown. The system automatically takes into account changes of a controllable line length that expands functionalities of a measuring instrument.

Текст научной работы на тему «Прозрачномеры-газоанализаторы на двухволновом полупроводниковом лазере»

СРЕДСТВА ИЗМЕРЕНИЙ

УДК 621.375.826

ПРОЗРАЧНОМЕРЫ-ГАЗОАНАЛИЗАТОРЫ НА ДВУХВОЛНОВОМ ПОЛУПРОВОДНИКОВОМ ЛАЗЕРЕ

Козлов В.Л., Кугейко М.М.

Белорусский государственный университет, г. Минск, Республика Беларусь

Приводятся методики построения прецизионных измерителей оптических характеристик с использованием двухволновых лазеров на основе базисного и нефелометрического методов. Показана устойчивость системы к изменениям аппаратурных констант, влиянию окружающей среды, загрязнению оптики. Система автоматически учитывает изменение длины контролируемой трассы, что расширяет функциональные возможности измерителя. (E-mail: KozlovVL@bsu.by)

Ключевые слова: измеритель оптических характеристик, двухволновой полупроводниковый лазер, базисный и нефелометрический методы.

Введение

Для повышения точностных и эксплуатационных характеристик информационно-измерительных и диагностических систем различного функционального назначения необходимо использование новых подходов, принципов и функциональных элементов. Одним из путей решения этой задачи является использование новых типов источников зондирующего излучения, к которым можно отнести двухволно-вые полупроводниковые лазерные диоды [1, 2]. Двухволновые лазеры имеют ряд существенных преимуществ перед дискретным набором лазерных источников, заключающихся в более высокой стабильности разности длин волн генерации, одинаковых флуктуаци-ях параметров излучения генерируемых сигналов, отсутствии проблем согласования оптических путей зондирующих сигналов. Эти преимущества обеспечивают более высокую эффективность работы измерительных систем на основе двухволновых лазеров перед аналогичными одноволновыми лазерными измерителями и перед многочастотными измерительными системами, в которых используются не-

сколько различных лазерных источников зондирующего сигнала.

Другим путем улучшения параметров лазерных измерительных систем является использования принципов концепции «безаприорно-сти» [3, 4], заключающейся в максимальном исключении методических погрешностей, обусловленных нестабильностью аппаратурных констант, приемо-передающего тракта, окружающей среды и исследуемых физических процессов. Данный подход снижает проблему априорной неопределенности, облегчает калибровочные измерения и используется для построения систем измерения оптических характеристик неоднородных рассеивающих сред на основе базисного и нефелометрического методов [4-6] для измерения оптической толщины, прозрачности, коэффициента обратного рассеяния, индикатрисы рассеяния подстилающей поверхности. Основная причина ограничения точности измерений в таких системах обусловлена неточным совмещением пучков зондирующего излучения, флуктуациями разностной частоты излучений и неточным знанием длины контролируемого участка среды. Использование в качестве источников зондирующего из-

лучения двухволновых лазеров позволяет устранить указанные недостатки. В литературных источниках отсутствуют сведения об использовании двухволновых лазеров в измерительных системах, построенных на основе концепции «безаприорности». Двухволновые полупроводниковые лазерные диоды с асимметричной квантово-размерной гетероструктурой [2] обеспечивают генерацию на двух различных оптических длинах волн. Переключение длины волны излучения в импульсе с на Х2 происходит при скачкообразном изменении амплитуды тока накачки в импульсе с 11 на 12. Длительность электрических импульсов и, соответственно, импульсов излучаемого света на разных длинах волн может быть сделана достаточно малой, менее единицы наносекунд. Разность длин волн генерации АХ = X] - Х2 достигает значений 10-90 нм. Если использовать терморегулятор и стабилизировать амплитуду тока инжекции, то достигается высокая стабильность разности длин волн генерации. Ниже приводятся принципы использования двухвол-новых лазеров в прецизионных измерителях оптических характеристик, построенных на основе концепции «безаприорности».

Нефелометрический измеритель прозрачности рассеивающей среды

Рассмотрим измеритель оптической толщины газовой струи с быстропеременной концентрацией в условиях ограниченных пространств. В данном случае использование локационного принципа, основанного на регистрации рассеянного в обратном направлении света, не приемлемо. В таких условиях рассеяние от близлежащих областей вызывает насыщение фотодетектора, что резко снижает точность регистрации или даже приводит к его повреждению. Использование же квазибистатической схемы измерений для исключения указанного выше фактора насыщения невозможно из-за образования теневой зоны, которая может даже превышать ширину исследуемой струи. Для базовых систем характерно влияние на результат измерений состояния среды, окружающей струю, загрязнения или запотевания оптики. Предлагаемая методика это устраняет, а также не требует калибровки системы [4]. Функциональная схема нефелометрического измерителя прозрачности рассеивающей среды на двухволновом лазере представлена на рисунке 1.

Рисунок 1 - Функциональная схема нефелометрического измерителя прозрачности рассеивающей среды на двухволновом лазере

Система содержит: двухволновой лазер, блок питания лазера БП, зеркало З, первый приемник излучения П1, второй приемник излучения П2, блок Процессор, блок рециркуляции РЦ.

Через исследуемую среду двухволновым лазером посылается импульсное зондирующее излучение на длине волны Х\. Излучение, рассеянное в первой и второй граничных точках и Я2 контролируемого участка среды под углом ф, величины которых «^(^оД) и ^«ДоД), рс-гистрируется приемниками П1 и П2 соответственно, и поступают в процессор. Затем зондирующее излучение отражается от зеркала и осуществляет обратный проход через исследуемую среду.

На приемниках П1 и П2 регистрируются значения сигналов «^(ДД) и кото-

рые записываются в память процессора. Из-за разных длин трасс рассеянное излучение будет поступать на приемники в разные моменты времени, поэтому система цифровой обработки процессора позволяет эффективно разрешить и идентифицировать эти сигналы. Для величин сигналов при прямом проходе среды на длине волны рассеянных под углом ф к направлению посылки в точках Д и Я2, восстановленных на квадрат расстояния, можно записать следующие выражения:

=А\ /', стф(/^,)'/(/^„./^,)'/■ ](/^,

И2) =А2Р! ^(К2)Ти(К0Л)Ти(К,Д2)Ти(К2Д5),

где Аь А2 - аппаратурные константы приемников П1 и П2 соответственно; 1'\ - мощность излучения на длине волны X] при прямом проходе, стф - коэффициент рассеяния под углом ф, Д, — координаты местонахождения приемников П1 и П2 соответственно; 11,. Д - координаты местонахождения соответственно лазера и зеркала; Ти(Яг, Р,) = ехр{-е-прозрачности участков [Д, |. 7,7 = 0...5 на длине волны

При обратном проходе среды излучением на длине волны отраженным от зеркала, величины сигналов, рассеянных под углом (180°-ф) в точках Д и Я2, имеют вид:

•4.1 Д,- Я\) = 1 / V» 1» т^ 1)'/ /. 1 (/¿2-1 ^) / /. 1Д Д) / /. 1Д -/<4)-

<$мДз, Я2) =Л2Р2 <31%сир(К2)Тх1(1{2,-ЯзДмДгД),

где Р2 - мощность излучения на длине волны X] при обратном проходе.

Отношение сигналов, рассеянных в точке ,равно:

1 ЗДД)

Последнее выражение можно записать следующим образом:

1п м = 1п в + 1п сх + тх1(/гь/г2),

(1)

где В=Р1ТХ1{Р0Д1)1Р2Ти{Р2Дъ)-

С,= ст0(Я,)/ст,^(Я,): 1МЯ2) = 8„(ДДНД-Д). Для сигналов, рассеянных в точке Я2:

м _^ д д)_ /дх/у^ддд.дд).

•4.1 Д Д ) а180-,р Д )ТЫ Д Дз )

Данное выражение можно записать в следующем виде:

(2)

где Ы2 = Л',^. ВДл Д, Д); С2 = Стф Д)/стш.ф(Д2).

Решение системы линейных уравнений (1) и (2) относительно т« равно:

т ж А)=1 ь^ДЬ^ДД). иК 15 21 зЛ)

+1п%дк

РД)

%Д)с>180 _Ф(Д2)

(3)

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

*„(»»)=I ы^ЖЛ2А1ДД) {■

(4)

Из выражения (4) легко перейти к коэффициенту ослабления:

ги(ЯМ = „ 1 „> У'^У'У'^У, ■ (5)

2а-я,) зД)

Затем через исследуемую среду двухвол-новым лазером посылается импульсное зондирующее излучение на длине волны Х2. Как и в предыдущих измерениях, при прямом и обратном прохождении т. R1, Я2 рассеянные излучения ВеЛИЧИНОЙ ^2(^0,^1) ¿ЪаДзД) И S/2(Rl,.R2). Л'/2 ДзД) записываются в процессор. Аналогично приведенным выше выкладкам, получаем коэффициент ослабления на длине волны '/.<

■<:/2( 1^.11) =

1

^.гДД) А2Д3Д2) . ^

2(R2-R1) Д,Д ) Д2 ДД )

Так как длины волн 'к \ и Х2 различаются незначительно, вычисление среднего значения коэффициента ослабления в соответствии с выражением:

В (ВД) =

8}лДД) + 8,2ДД)_ 1

Ш

(7)

хрп ^лДоД) А1ДД2) +

^лДД) ДлДзД) , 2 До Д1) А 2 Дз Д 2 ) j ^ДоДг) А2ДД)'

где AR= R2 - R1 - длина контролируемой трассы, обеспечит повышение точности измерений в y¡2 раз по сравнению с одноволновыми измерениями.

Одновременно с измерением интенсивности рассеянного излучения в системе осуществляется режим оптоэлектронной рециркуляции последовательно относительно первой и второй граничных точек. Зарегистрированный сигнал S,](R„, Ri), рассеянный первой граничной точкой R1, запускает блок рециркуляции. В результате замыкания петли оптической обратной связи в системе устанавливается процесс рециркуляции, период т которой определяется оптической задержкой излучения на дистанции при постоянной электрической задержке в блоке рециркуляции. Период (частота f1) рециркуляции относительно первой граничной точки R1 будет определяться следующим образом:

где, ¿ори = ^/с - время задержки излучения на трассе, Ь\ - длина трассы, состоящая из участка Д0..Д и участка Л1...Л4, с - скорость света в воздухе, ¿е - время электрической задержки. Время электрической задержки ¿е в блоке рециркуляции выбирается таким образом, чтобы оно было больше суммарной задержки излучения на всей контролируемой трассе. Следовательно, все импульсы рассеянного излучения при прямом и обратном проходе трассы будут зарегистрированы приемниками, и только после этого лазер посылает на трассу очередной оптический импульс.

Затем процессор, зарегистрировав сигнал S, ] (11,. R2), рассеянный второй граничной точкой R2, запускает режим рециркуляции относительно точки R2. Период (частота /2) рециркуляции относительно второй граничной точки R2 будет определяться следующим образом:

\lfl~ t0pt2 + U

(9)

где ¿ор12 = Ь2/с - время задержки излучения на трассе; Ь2 - длина трассы, состоящая из участка R0...RЬ участка R1...R2 и участка R2...R5. Из выражений (8, 9) можно найти длину контролируемой трассы AR по разности частот рециркуляции из выражения

(

AR = R2-Rl=c

J_

v/2

J_

А

(10)

1//Í - Xl - í0ptl + te,

(8)

где /1, /2 - частоты рециркуляции относительно первой и второй граничных точек, соответственно. Полученное значение длины контролируемой трассы AR учитывается в (5-7) при расчете оптических характеристик.

Если разность измеренных коэффициентов 8и(ДД2) и гк2(К1,К2) имеет большое значение, то наблюдается зависимость коэффициента ослабления от длины волны, что можно использовать для измерения концентрации газа. Ширина спектральных линий поглощения различных газов составляет величину в доли и единицы нанометров, поэтому для измерения концентрации газа, длины волн генерации двухволнового лазера выбираются таким образом, чтобы длина волны 'к \ находилась в центре полосы поглощения контролируемого газа, а длина волны /.2 - вне полосы поглощения газа. Представим коэффициенты ослабления е(л) на длинах волн Х\ и Х2 в виде:

г,1л=СМ\) + т1\ 8Х, =сда2) + (за2),

(11)

где К('к\) - коэффициент поглощения исследуемой газовой компоненты на длине волны К(Х2) - коэффициент поглощения исследуемой газовой компоненты на длине волны Х2, Сх -концентрация исследуемой газовой компоненты; Р(?ч), - суммарные коэффициенты поглощения и рассеяния на А,ь Х2 всеми другими компонентами, присутствующими в среде.

Решение системы уравнений (5, 6) с учетом выражений (11) относительно Сх будет иметь вид:

С = —

1

-[1п

1п ^ДЛ) ЛлДЛ^ 1

БМЛ) • ЯМА? [К(\)-К(Х2)]

Р^-РСМ

Так как длины волн генерации двухволно-вого лазера выбираются близкими по значению X] ~ Х2 с разницей в единицы нанометров, то Р(?ц)-Р(Л„) ->0- следовательно, искомая концентрация газа будет равна

С =

1

2 АКАЯ ЗД.ЛЫл^зЛ) $12(^0 А) -^2(^3^2)

(12)

где АК = КО,,) - К(Х2) ~ дифференциальный коэффициент поглощения; АЯ - длина контролируемой трассы.

Базисный измеритель прозрачности рассеивающей среды

Наиболее точные методы измерения прозрачности участков рассеивающих сред - базисные, суть которых в измерении проходящего излучения. При этом погрешности, связанные с нестабильностью аппаратурных констант при-емо-передающего тракта, с неточным совмещением пучков зондирующего излучения, присущи и этому методу. Использование в качестве источника излучения двухволнового лазера позволяет устранить упомянутые недостатки. Функциональная схема базисного измерителя прозрачности на двухволновом лазере представлена на рисунке 2.

Рисунок 2 - Функциональная схема базисного измерителя прозрачности на двухволновом лазере

Система содержит: двухволновой лазер Л, блок питания лазера БП, первый приемник излучения П1, второй приемник излучения П2, блок перемещения лазера БПЛ, блок питания лазера БП, блок Процессор.

В качестве источника излучения используется двухволновой полупроводниковый лазерный диод с асимметричной квантово-размерной гетероструктурой, описанный выше [2]. Двухволновой лазер имеет 2 полупрозрачных зеркала резонатора, поэтому оптическое излучение выводится в обе стороны. Через исследуемую среду двухволновым лазером посылается импульсное зондирующее излучение на длине волны . При одновременной посылке лазером световых потоков в противоположных направлениях, когда лазер находится на одной из границ участка (например, Е), для величин потоков, зарегистрированных вторым приемником излучения (измерительного ]> ) и первым приемником излучения (компенсационного р. )

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

можно записать следующие выражения:

(К г) = АРошТи (Я к + ЩТи (Я + М, г),

( п+ж Т

Т^Д + ДД^ехр]- | гп(г)с1г ;

где

ГХ1(Д + ДД,г) = ехр]- | еи(г)л1

Д+АЙ

где А, - аппаратурная константа приемника излучения П2; - величина светового по-

тока в направлении приемника П2; '/]_] - прозрачность участка среды на длине волны /ч; - коэффициент ослабления на длине волны

РУ/1(М,г]) = А2Р(}Г/1 ехр

л

-|е и{г)с1г

где А ^ - аппаратурная константа приемника излучения П1; - величина светового по-

тока в направлении приемника П1 .

Затем двухволновой лазер посылает через исследуемую среду импульсное зондирующее излучение на длине волны Х2. Как и в предыдущем измерении для величин потоков, зарегистрированных вторым приемником излучения

(измерительного р_ ) и первым приемником

излучения (компенсационного Р. ) можно за-

•2

писать следующие выражения:

Р1Х2 (К,г) = А1РШТХ2(К,К + тГг2(Р + М,г), РЗХ2(Я,г1) = А2Р02ХТХ2(Г1,Я).

С помощью блока перемещения лазер передвигается в точку Я+АЯ и посылает в исследуемую среду коллинеарно направлениям посылки из точки Е в противоположных направлениях импульсное зондирующее излучение на длине волны Для величин измерительного потока Р4„ , регистрируемого приемником П2,

и компенсационного Р2Х , регистрируемого

приемником П1, можно записать:

Р2Ч (К + АЛ, г) = АР0Щ ехр ] - I еХ1 (г)с/г \ =

Я+&Я

РЛ/1(Р + АР,г) = А2Р()2/1 ехр

л+дд

= АРпхТк +

Затем лазер посылает через среду импульсное излучение на длине волны Х2. Аналогично получаем выражения для Р4} и Р2} :

Р2Хг (Я + АД,г) = АР01ХТХ2(Я + ДД,г),

РА. (Я + АВД = А.Р^Т^ЩТ.МЛ + АД).

Отношения Р. к р и Р к Р соответ-

1 1/-1 1 3/ч 1 2 >4 1 4/4

ственно равны:

}Т_Рщ _ЛРощТи(КК + ЩТи(К + АЯ,г) 1 л,

^2 =

- 3^ р

АРтти(я+М1,г)

РаК А Р^ТхМЛШЯА + т

2Х1

Последние выражения можно представить в виде системы линейных уравнений:

к

In TV] = In 5 + SAR, In iV2 = In 5 - sAR,

(13)

где In В = In [A1P0mTu(R+AR,r)/A2P02UTu(riJi)].

Решение системы (13) относительно 8 xi имеет вид:

1 Р Р -In- 11

2АR P2XP3h

(14)

Прозрачность же на участке AR на длине волны равна

р Р

1 Щ1 4Â.1 р. р

1/2

(15)

Аналогично получаем коэффициент ослабления и прозрачность на участке АЛ на длине волны Х2

1 Р Р

Р^Ли

р р

1 Г/.21 4Â2

Р Р

1 2Х2 312

1/2

(16)

Очевидно, что вычисление среднего значения и '¿/ 2 обеспечит повышение точности по сравнению с одноволновыми измерениями. Если относительная разность и '¿/ 2 имеет

большее значение, то это можно использовать для измерения концентрации газа аналогично нефелометрическому измерителю. Для измерения концентрации газа длины волн генерации двухволнового лазера выбираются таким образом, чтобы длина волны X] находилась в центре полосы поглощения контролируемого газа, а длина волны Х2 - вне полосы поглощения газа. Коэффициенты ослабления е(/.) на длинах волн X] и к2 могут быть представлены выражениями (11). Решение системы уравнений (14-16) с учетом выражений (11) относительно Сх будет иметь вид:

С = —

1

1

In Vk

Р Р

Р Р

1 2л/ АХг

[K{\)-KQ,2)]

ß(Äa)-ß(A,2) .

Так как длины волн генерации двухволнового лазера выбираются близкими по значению то Р(^)-Р(Х2) -»О, следовательно, искомая концентрация газа будет равна

е.. - -

1

2АKAR

P P P P

Jn 3>Ч 4Я-2

P P P Р

(17)

где АК = К(Х\) - К(Х2) - дифференциальный коэффициент поглощения.

Заключение

Предложены и обоснованы принципы использования двухволновых лазеров в прецизионных измерителях оптических характеристик, построенных на основе концепции «безаприор-ности». Как видно из приведенных выше выкладок, выражения для оптических характеристик (Т, s) и концентрации газа Сх как для базисного, так и для нефелометрического измерителей не содержат ни аппаратурных констант, ни энергии посылаемого излучения, ни параметров, отражающих влияние окружающей среды. Это означает, что отсутствует необходимость калибровки измерительной системы и установки аппаратурных константы ввиду их отсутствия в алгоритмах, которые получены без использования каких-либо допущений или пренебрежений этими константами. Устойчивость системы к изменениям аппаратурных констант означает и устойчивость к изменениям окружающей среды (появление пыли, влаги и др.), загрязнению оптики.

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

Возможность реализации двухчастотного нефелометрического метода в ограниченных по

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

Список использованных источников

1. Ikeda, S. Evidence of the wavelength switching caused by a blocked carrier transport in an asymmetric dual quantum well laser / S. Ikeda,

A. Shimizu // Appl. Phys. Lett. 1991. -V. 59. -P.504-506.

2. Афоненко, А.А. Полупроводниковый лазер / А.А. Афоненко, В.К. Кононенко, И.С. Ма-нак // Патент Республики Беларусь №1385. 1996.

3. Кугейко, М. М. Теория и методы оптико-физической диагностики неоднородных рассеивающих сред. / М. М. Кугейко, Д. М.Оношко - Минск : БГУ, 2003. - 185 с.

4. Кугейко, М.М. Лазерные системы (в условиях априорной неопределенности). / М.М. Кугейко. - Минск : БГУ, 1999. 196 с.

5. Сергеев, Н.М. Измерение прозрачности атмосферы с использованием двухлазеров / Н.М. Сергеев // Тез.докл. VI Всесоюзного симпозиума по лазерному и акустическому зондированию атмосферы. Томск. - Ч.1. -1980. - С.123-125.

6. А.с. СССР №1523974, МКИ G 01 N 21/47. Способ определения прозрачности участка рассеивающей среды / Б. Б. Виленчиц [и др.] Опубл. 1989 г. Бюл. №43.

Kozlov V.L., Kygeiko M.M.

Transparency meters and gas analyzers based on two-wavelength semiconductor lasers

Construction techniques of precision measuring instruments of optical characteristics on the basis of two-wave lasers with use of basic and nephelometer methods are presented. System stability to changes of hardware constants, influence of an environment, pollution of optics is shown. The system automatically takes into account changes of a controllable line length that expands functionalities of a measuring instrument. (E-mail: KozlovVL@bsu.by)

Key Words: two-wavelength semiconductor laser, optical characteristics measuring instrument, basic and nephelometer methods.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Поступила в редакцию 20.10.2011.

i Надоели баннеры? Вы всегда можете отключить рекламу.