Например, Бобцов

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

УДК 53.072; 681.3
ИМИТАЦИОННАЯ МОДЕЛЬ ОПТИКО-ЭЛЕКТРОННОГО ИЗМЕРИТЕЛЯ ВЫСОТЫ НИЖНЕЙ ГРАНИЦЫ ОБЛАКА
А.В. Демин, К.В. Константинов
В статье предложен метод имитационного моделирования оптико-электронного комплекса для измерения высоты нижней границы облаков способом просмотра активностей. Разработана и апробирована имитационная модель измерения высоты нижней границы облаков. Ключевые слова: имитационная модель, высота нижней границы облаков.
Введение
Одним из важных режимов обеспечения безопасности полетов является обеспечение безаварийной посадки летательных аппаратов (ЛА) на аэродром, которая, помимо технического состояния ЛА, определяется метеорологической обстановкой в зоне посадки (МОАр). Контроль за МОАр обеспечивается аэродромным оборудованием, важным компонентом которого является система измерения высоты нижней границы облаков (НГО) в зоне посадке. Качество и достоверность работы измерительного комплекса (ИК) во многом определяется выбранной моделью облака и, в частности, моделью НГО. Процесс измерения высоты НГО может быть представлен в виде следующей последовательности: ИК⇒атмосфера⇒НГО⇒атмосфера⇒ИК.
Наряду с известными методами определения высоты НГО, наиболее объективным является ИК, выполненный как оптико-электронный комплекс (ОЭКНГО).
Постановка задачи
Представление процесса измерения высоты НГО математической моделью является одним из важных моментов как при проектировании, так и при последующей эксплуатации ОЭКНГО. Целью настоящей статьи является построение имитационной модели процесса измерения высоты НГО «атмосфера⇒НГО⇒атмосфера» для оценки ожидаемых показателей работы ОЭКНГО в процессе проектирования и их корректировки при эксплуатации. На рис. 1 представлена структурная схема ОЭКНГО, оптическая схема которого выполнена по коаксиальной схеме, т.е. объектив излучателя и объектив приемной системы выполнены как общий для обоих каналов. Лазерный излучатель с λ=0,9 мкм имеет дополнительную оптическую систему для изменения угловой расходимости и систему модуляции излучения. Перед фотоприемным устройством (ФПУ) установлен интерференционный фильтр с λ=(0,9±0,05) мкм, тем самым обеспечивается минимизация воздействия на ФПУ постороннего излучения. Контроллер выполняет функцию программного управления комплексом в соответствии с математической моделью функционирования ОЭКНГО в параметрической или аналитической форме. Блок обработки данных выполняет функцию расчета высоты НГО.
Метод моделирования
Из всех известных методов моделирования был выбран метод имитационного моделирования, основанный на моделировании непосредственно алгоритмами функционирования подсистем моделируемой системы, или имитационное моделирование системы способом просмотра активностей его подсистем (ИМАК). Метод позволяет реализовать моделирование любой системы, особенностью которой является то, что функциональные действия каждой из ее подсистем различны, причем для выполнения дей-

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

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

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

модель – НГО

Исходная задача на моделирование
ОЭКНГО

результат

ЦВМ
АБ1 - ОЭКНГО

МДВ в вертикальном направлении

АБ2 – среда

АБ3 – помехи

АБ4 – НГО
Рис. 2. Имитационная модель измерения высоты нижней границы облаков

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

Имитационная модель нижней границы облака
(имитатор)
Имитационная модель среды
()
Имитационная модель ОЭКНГО
(имитатор)
Имитатор помех

Алг1
Алг2 Алг3 Алг4

Mt1 Mt2 Mt3 Mt4

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

ПО начала
имитации

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

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

ПО окончания
имитации

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

Известно, что совокупность атмосферных аэрозолей следует рассматривать как систему твердых или жидких дисперсных частиц в воздухе, размеры которых в зависимости от высоты и температуры в месте их локализации возможны в пределах 10-3–102 мкм. Совокупность атмосферных аэрозолей можно классифицировать на характерные типы – облака, туманы, дымки, осадки, ядра конденсации. При этом аэрозоль – одна из наиболее оптически активных и наименее изученных составляющих атмосферы – присутствует практически на всех высотах [2].
Выделяются следующие источники аэрозоля в атмосфере:
− морская, океаническая и земная поверхность;
− вторичная генерация в атмосфере за счет фотохимических реакций;
− продукты индустриальной и хозяйственной деятельности;
− продукты вулканической деятельности (газ, пыль), дым лесных пожаров;
− космическая пыль и продукты сгорания метеоров. Облака бывают капельножидкими и кристаллическими. Капельножидкие облака
состоят из сферических частиц, кристаллические облака имеют частицы неправильной формы – кристаллики льда.
Основные процессы, регулирующие распределения аэрозоля в атмосфере:
− турбулентная диффузия и конвективный перенос аэрозольных частиц;
− синоптические факторы – циклонные факторы, температурные инверсии, осадки;
− микромасштабные процессы – коагуляция, седиментация, абсорбция водяного пара. Решить задачу о вертикальном распределении аэрозольных частиц в атмосфере
можно, только зная источники и процессы, регулирующих их эволюцию. Одним из факторов обеспечения безопасности при взлете–посадке является высо-
та НГО, которые формируются из атмосферных аэрозолей размером от 0,5 мкм до 10 мкм [3]. Возможно следующее представление моделей НГО [3]:
− моделирование оптических свойств атмосферы путем задания моделей микроструктуры атмосферных аэрозолей;
− в соответствии с теорией Ми.
− представление НГО как объекта, обладающего определенными оптическими и микрофизическими свойствами. Учитывая, что облака представимы набором атмосферных аэрозолей в форме сфе-
рических частиц, построение модели НГО, под которой будем понимать упорядоченный по высоте и спектру частот числовой массив объемных коэффициентов взаимодействия и компонент матрицы рассеяния, невозможно осуществить без достоверной количественной информации о микрофизических свойствах ансамбля аэрозольных частиц, статистически обоснованного для заданной геофизической ситуации.
Концентрация и функция распределения аэрозольных частиц по размерам, форме частиц и их химическому составу, как правило, известны на основе экспериментальных измерений. Опираясь на экспериментальные данные по аэрозолям, модель НГО можно представить в виде целочисленной модели:

Nобл

=

⎤ ⎥ Kσ ⎦⎥

×



⎛ × ⎜⎜⎝

π×r2 ×h Vаэр

⎞ ⎟⎟⎠

⎡ ⎢ ⎣⎢

+

1,

(1)

где Nобл – число аэрозолей в НГО, необходимое для регистрации в ОЭКНГО; r – радиус

светового пучка на высоте H (радиус облучаемого участка облака); Kσ – коэффициент,

определяющий рассеяние излучения в НГО; Kα – коэффициент, определяющий ослабление излучения по трассе; Vаэр – объем аэрозоля; h – глубина облачного слоя.

Поскольку для регистрации высоты НГО в измерителе необходимо выдерживать

определенное отношение сигнал/шум при заданной мощности и расходимости излуче-

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


∫σ(λ) = N πr 2 ⋅ K p (r, m, λ) f (r)dr ,

(2)

0

где N – число аэрозольных частиц в единице объема; K p (r, m, λ) – фактор эффектив-

ного рассеяния излучения с длиной волны λ на частице радиусом r, имеющей ком-

плексный показатель преломления m; f(r) – функция распределения аэрозольных час-

тиц по размерам.

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



=

N

(h)

×

⎛ ⎜⎜⎝

⎤ ⎥ ⎦

dаэр dл

⎡ ⎢ ⎣

+1⎠⎟⎟⎞ ,

(3)

где dаэр – среднее значение диаметра аэрозольной частицы; dл – диаметр пучка отра-

женных от частицы лучей, возвращающихся в ОЭКНГО (не более 0,25 dаэр); N (h) – ко-

личество частиц аэрозоля в НГО в зависимости от высоты (рассчитанное из экспери-

ментальных данных).

Под коэффициентом аэрозольного ослабления понимается величина αa (λ) , зави-

сящая от длины волны λ и определяемая из закона Бугера в дифференциальной форме.

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

вертикальном направлении на момент измерения высоты НГО.

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

max

{

Nизл

}

=

⎤ ⎥ ⎥⎦

(

Dвых.зр. (λ)2

)

2

⎡ ⎢ ⎣⎢

+

1,

(4)

где Dвых.зр. – диаметр выходного зрачка объектива ОЭКНГО. Для обеспечения необходи-

мого отношения сигнал/шум η число «единичных лучей» требуемой мощности (Wу.е. – отношение мощности вышедшего излучения к минимальной мощности, определяемой

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

( )NизлW

=

η

×

⎛ ⎜ ⎝⎜

⎤ ⎥ ⎥⎦

K

σ

×

K

−1 α

×π Vаэр

×

r

2

×

dаэр

⎡ ⎢ ⎣⎢

+

⎞ 1⎟
⎟⎠

×

Wу.е.

,

(5)

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

В соответствии с соотношением (4) значение Dвых.зр. как оценка будет определять-

ся выражением

Dвых.зр. = λ × NизлW .

(6)

Таким образом, из соотношения (6) получаем диаметр выходного зрачка Dвых.зр. . Однако, как следует из соотношений (3), (5), этот диаметр зависит от числа аэрозолей в момент измерения, т.е. от метеоусловий (типа оптической погоды). Поэтому расчет Dвых.зр. производится по средним условиям, а нужное количество излучаемых лучей NизлW обеспечивается с помощью подбора числа зондирований. В дальнейших исследованиях будет представлена целочисленная модель излучающей и приемной части оптико-электронного комплекса измерителя высоты нижней границы облаков.
Заключение
Предложен метод имитационного моделирования способом просмотра активностей. Разработана целочисленная модель «атмосфера⇒НГО⇒атмосфера», на основании которой получены соотношения для определения параметров оптико-электронного комплекса в общей модели процесса измерения высоты НГО.
Литература
1. Максимей И.В. Имитационное моделирование на ЭВМ. – М.: Радио и связь, 1988. 2. Зуев В.Е., Зуев В.В. Дистанционное оптическое зондирование атмосферы. – СПб:
Гидрометеоиздат, 1992. 3. Зуев В.Е., Кабанов М.В. Перенос оптических сигналов в земной атмосфере (в усло-
виях помех). – М.: Сов. радио, 1977.

Демин Анатолий Владимирович – Санкт-Петербургский государственный университет ин-

формационных технологий, механики и оптики, доктор

технических наук, профессор, dav_60@mail.ru

Константинов Константин – Санкт-Петербургский государственный университет ин-

Владимирович

формационных технологий, механики и оптики, аспирант,

const_cv@mail.ru