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

УДК 53.072; 681.3

ИМИТАЦИОННАЯ МОДЕЛЬ ОПТИКО-ЭЛЕКТРОННОГО ИЗМЕРИТЕЛЯ ВЫСОТЫ НИЖНЕЙ ГРАНИЦЫ ОБЛАКА

А.В. Демин, К.В. Константинов

В статье предложен метод имитационного моделирования оптико-электронного комплекса для измерения высоты нижней границы облаков способом просмотра активностей. Разработана и апробирована имитационная модель измерения высоты нижней границы облаков. Ключевые слова: имитационная модель, высота нижней границы облаков.

Введение

Одним из важных режимов обеспечения безопасности полетов является обеспечение безаварийной посадки летательных аппаратов (ЛА) на аэродром, которая, помимо технического состояния ЛА, определяется метеорологической обстановкой в зоне посадки (МОАр). Контроль за МОАр обеспечивается аэродромным оборудованием, важным компонентом которого является система измерения высоты нижней границы облаков (НГО) в зоне посадке. Качество и достоверность работы измерительного комплекса (ИК) во многом определяется выбранной моделью облака и, в частности, моделью НГО. Процесс измерения высоты НГО может быть представлен в виде следующей последовательности: ИК^атмосфера^НГО^атмосфера^ИК.

Наряду с известными методами определения высоты НГО, наиболее объективным является ИК, выполненный как оптико-электронный комплекс (ОЭКнш).

Постановка задачи

Представление процесса измерения высоты НГО математической моделью является одним из важных моментов как при проектировании, так и при последующей эксплуатации ОЭКНГО. Целью настоящей статьи является построение имитационной модели процесса измерения высоты НГО «атмосфера^НГО^атмосфера» для оценки ожидаемых показателей работы ОЭКНГО в процессе проектирования и их корректировки при эксплуатации. На рис. 1 представлена структурная схема ОЭКНГО, оптическая схема которого выполнена по коаксиальной схеме, т.е. объектив излучателя и объектив приемной системы выполнены как общий для обоих каналов. Лазерный излучатель с ^=0,9 мкм имеет дополнительную оптическую систему для изменения угловой расходимости и систему модуляции излучения. Перед фотоприемным устройством (ФПУ) установлен интерференционный фильтр с ^=(0,9±0,05) мкм, тем самым обеспечивается минимизация воздействия на ФПУ постороннего излучения. Контроллер выполняет функцию программного управления комплексом в соответствии с математической моделью функционирования ОЭКНГО в параметрической или аналитической форме. Блок обработки данных выполняет функцию расчета высоты НГО.

Метод моделирования

Из всех известных методов моделирования был выбран метод имитационного моделирования, основанный на моделировании непосредственно алгоритмами функционирования подсистем моделируемой системы, или имитационное моделирование системы способом просмотра активностей его подсистем (ИМАК). Метод позволяет реализовать моделирование любой системы, особенностью которой является то, что функциональные действия каждой из ее подсистем различны, причем для выполнения дей-

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

Рис. 1. Структурная схема оптико-электронного измерительного комплекса высоты

нижней границы облака

На рис. 2 представлена имитационная модель измерения высоты нижней границы (ИМНГО) в виде активируемых блоков, где приняты следующие символьные и функциональные обозначения: АБ\ - модель измерителя; АБ2 - модель атмосферы между НГО и ОЭКНГО, а именно, модельное или параметрическое представление (метеосводка, географическое место проведения ИМ ОЭКВНГО, время суток, месяц, год, высота светила); АБ3 - модель помехи и возмущения, а именно, описание влияния атмосферы и приемопередающей аппаратуры на оптический сигнал ОЭКНГО; АБ4 - модель НГО.

Рис. 2. Имитационная модель измерения высоты нижней границы облаков

Особенностью ИМАК является то, что все функциональные действия подсистем могут быть синхронизированы и рассмотрены в процессе моделирования независимо друг от друга и в реальном масштабе времени [1].

Имитационная модель нижней границы облака (имитатор)

Имитационная модель среды

Алг1

Я

Алг2

С

Имитационная

модель ОЭКНГО (имитатор)

Имитатор помех

Алг,

Алг4

шг1

ши

ши

шгл

II

Управляющая программа имитационного моделирования

7Т

7ч

/ ПО сбора статистики; критерий оценки

ПО начала имитации Алгоритмы: преобразования; анализа; регистратора; определения характеристик объекта

Рис. 3. Блок-схема имитационного моделирования оптико-электронного комплекса измерения высоты нижней границе облака способом просмотра активностей

В соответствии с ИМНго на рис. 3 представлена блок схема имитационного моделирования ОЭКНго, где приняты обозначения: Mt - модельное время имитируемого блока; АЛГ - алгоритм имитации (математическая модель) имитируемого блока; ПО -программное обеспечение имитационного моделирования по принятым аналитическим соотношениям. Начало имитации обусловлено установкой значений начальных состояний компонентов модели. Проверка выполнимости условия инициализации активностей выполняется либо путем определения параметров модели, либо путем вычисления момента активизации блока или параметров модели, либо путем вычисления момента активизации блока. Поскольку все функциональные действия в подсистемах ИМИУС происходят независимо друг от друга, то и модельное время Mt для каждого АБг разное.

Построение математической модели тракта «атмосфера^НГО^атмосфера»

Предполагаем, что функционирование ОЭКНГО можно описать линейным оператором. Рассмотрим основы для построения математической модели тракта «атмосфе-ра^НГО^атмосфера» как наиболее важного звена в общей модели процесса измерения НГО при определении параметров комплекса.

Известно, что совокупность атмосферных аэрозолей следует рассматривать как систему твердых или жидких дисперсных частиц в воздухе, размеры которых в зависимости от высоты и температуры в месте их локализации возможны в пределах 10-3-102 мкм. Совокупность атмосферных аэрозолей можно классифицировать на характерные типы

- облака, туманы, дымки, осадки, ядра конденсации. При этом аэрозоль - одна из наиболее оптически активных и наименее изученных составляющих атмосферы - присутствует практически на всех высотах [2].

Выделяются следующие источники аэрозоля в атмосфере:

- морская, океаническая и земная поверхность;

- вторичная генерация в атмосфере за счет фотохимических реакций;

- продукты индустриальной и хозяйственной деятельности;

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

- космическая пыль и продукты сгорания метеоров.

Облака бывают капельножидкими и кристаллическими. Капельножидкие облака состоят из сферических частиц, кристаллические облака имеют частицы неправильной формы - кристаллики льда.

Основные процессы, регулирующие распределения аэрозоля в атмосфере:

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

- синоптические факторы - циклонные факторы, температурные инверсии, осадки;

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

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

Одним из факторов обеспечения безопасности при взлете-посадке является высота НГО, которые формируются из атмосферных аэрозолей размером от 0,5 мкм до 10 мкм [3]. Возможно следующее представление моделей НГО [3]:

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

- в соответствии с теорией Ми.

- представление НГО как объекта, обладающего определенными оптическими и микрофизическими свойствами.

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

Концентрация и функция распределения аэрозольных частиц по размерам, форме частиц и их химическому составу, как правило, известны на основе экспериментальных измерений. Опираясь на экспериментальные данные по аэрозолям, модель НГО можно представить в виде целочисленной модели:

ПХ^ ХИ +1, (1)

N =

обл

КX Ках

V

аэр

где ^бл - число аэрозолей в НГО, необходимое для регистрации в ОЭКнго; г - радиус светового пучка на высоте Н (радиус облучаемого участка облака); К - коэффициент, определяющий рассеяние излучения в НГО; Ка - коэффициент, определяющий ослабление излучения по трассе; ^эр - объем аэрозоля; И - глубина облачного слоя.

Поскольку для регистрации высоты НГО в измерителе необходимо выдерживать определенное отношение сигнал/шум при заданной мощности и расходимости излуче-

ния, то в зависимости от концентрации аэрозолей в облаке глубина облачного слоя, а также рассеяние и ослабление излучения могут быть различны. Тем самым возникает неопределенность в определении высоты НГО, которая может уменьшаться двумя путями - либо повышением мощности излучения, либо изменением размера зоны облученности. Безусловно, разумное сочетание мощности и размера зоны облучения, а также учет возможного значения коэффициентов рассеяния и ослабления определяют некий оптимум параметров ОЭКНГО.

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

ст(А) = N\пт2 • К р (г, т, А)/(г)йг ,

(2)

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

Представим соотношение (2) в следующем виде:

К с= N (И) х

г й - л

аэр +1

й - )

V- л

(3)

где ^аэр - среднее значение диаметра аэрозольной частицы; йл - диаметр пучка отраженных от частицы лучей, возвращающихся в ОЭКНГО (не более 0,25 ^аэр); N (И) - количество частиц аэрозоля в НГО в зависимости от высоты (рассчитанное из экспериментальных данных).

Под коэффициентом аэрозольного ослабления понимается величина а а (А), зависящая от длины волны А и определяемая из закона Бугера в дифференциальной форме. В этой связи коэффициент аэрозольного ослабления примем равным баллу видимости в вертикальном направлении на момент измерения высоты НГО.

Излучение на выходе ОЭКНГО также представим в виде набора дискретных лучей, предельное значение которого будет равно

тах

} =

(дДвых.зр. )

(А)2

+1,

(4)

где Двыхзр - диаметр выходного зрачка объектива ОЭКНГО. Для обеспечения необходимого отношения сигнал/шум п число «единичных лучей» требуемой мощности (Жуе. -отношение мощности вышедшего излучения к минимальной мощности, определяемой приемником) может быть определено следующим выражением:

Г

N

излW

= П х

VJ

КаХ К-1ах(пх Г 2 )х йаэр

V

аэр

1

хЖ„

где Ка - коэффициент аэрозольного ослабления.

В соответствии с соотношением (4) значение Двых.зр. ся выражением

Ах^

Д

вых.зр.

излЖ

(5)

как оценка будет определять-

(6)

Таким образом, из соотношения (6) получаем диаметр выходного зрачка А,ыхзр .

Однако, как следует из соотношений (3), (5), этот диаметр зависит от числа аэрозолей в момент измерения, т.е. от метеоусловий (типа оптической погоды). Поэтому расчет £ производится по средним условиям, а нужное количество излучаемых лучей

NИiJ¡W обеспечивается с помощью подбора числа зондирований. В дальнейших исследованиях будет представлена целочисленная модель излучающей и приемной части оптико-электронного комплекса измерителя высоты нижней границы облаков.

Заключение

Предложен метод имитационного моделирования способом просмотра активностей. Разработана целочисленная модель «атмосфера^НГО^атмосфера», на основании которой получены соотношения для определения параметров оптико-электронного комплекса в общей модели процесса измерения высоты НГО.

Литература

1. Максимей И.В. Имитационное моделирование на ЭВМ. - М.: Радио и связь, 1988.

2. Зуев В.Е., Зуев В.В. Дистанционное оптическое зондирование атмосферы. - СПб: Гидрометеоиздат, 1992.

3. Зуев В.Е., Кабанов М.В. Перенос оптических сигналов в земной атмосфере (в условиях помех). - М.: Сов. радио, 1977.

Демин Анатолий Владимирович

Константинов Константин Владимирович

Санкт-Петербургский государственный университет информационных технологий, механики и оптики, доктор технических наук, профессор, dav_60@mail.ru Санкт-Петербургский государственный университет информационных технологий, механики и оптики, аспирант, const_cv@mail.ru