ООО “НПП “Лазерные системы”
Методы и алгоритмы обработки данных ветрового когерентного доплеровского лидарного профилометра с коническим сканированием
В.Р. Ахметьянов, Д.Н. Васильев, М.А. Коняев,
О.А. Мишина, М.С. Пенкин, Г.А. Петров,
Я.А. Тезадов, И.В. Шаталов, И.Ф. Ширяев
Аннотации
В работе рассмотрены вопросы, связанные с программно-алгоритмическим обеспечением ветрового когерентного доплеровского лидарного профилометра, предназначенного для оперативного мониторинга ветровых характеристик атмосферы и, в частности, для определения такого метеоявления, как сдвиг ветра в районе взлетно-посадочной полосы аэропорта. Описан принцип работы лидарного профилометра с непрерывным коническим сканированием. Приведены методы и алгоритмы обработки лидарного сигнала, составляющие основу программно-алгоритмического обеспечения профилометра.
Ключевые слова: скорость и направление ветра, доплеровский лидар, методы и алгоритмы обработки сигнала.
The paper deals with questions related to the program-algorithmic software of the wind coherent doppler lidar profiler, which is intended for operational monitoring of the wind characteristics of the atmosphere and, in particular, to determine such meteorological phenomenon as wind shear in the area of the airport's runway. The principle of operation of the lidar profiler with continuous conical scanning is described. Methods and algorithms of processing of lidar signal, as the basis of the profiler program-algorithmic support, are considered.
Key words: wind speed and direction, doppler lidar, methods and algorithms for signal processing.
1. Введение
В самых различных сферах жизнедеятельности человека требуется информация о ветровых характеристиках атмосферы. В результате в настоящее время для определения таких параметров ветра, как его скорость и направление, все более широкое распространение получают дистанционные методы зондирования атмосферы. Благодаря развитию технологической базы и бурному росту вычислительных мощностей компьютерной техники активно развиваются лидарные методы, в основе работы которых лежит эффект Доплера. В соответствии с достигнутым техническим уровнем в нормативных документах ICAO [1] в настоящее время рекомендуется внедрение ветровых когерентных доплеровских лидаров (КДЛ) в систему метеорологического обеспечения авиационной безопасности [2].
Работа КДЛ построена на допущении, что лазерное излучение на частоте f по мере распространения в атмосфере рассеивается на аэрозольных частицах. С учетом скорости движения этих частиц частота обратно рассеянной волны изменяется согласно эффекту Доплера.
В рассматриваемом типе когерентных доплеровских ладарах приемник рассеянного аэрозольными частицами сигнала расположен в том же месте, где и передатчик. Следовательно, выражение для частотного сдвига может быть записано следующим образом:
,
где Vr - радиальная составляющая скорости ветра или иначе проекция мгновенного вектора скорости V={Vx, Vy, Vz} на направление зондирования. Таким образом
лидарный профилометр доплеровский когерентный
,
где - непосредственно измеряемый доплеровский сдвиг частоты.
Необходимо отметить, что в существующих и разрабатываемых КДЛ может использоваться как импульсное, так и непрерывное лазерное излучение. КДЛ с импульсным излучением, как правило, позволяет проводить зондирование атмосферы на расстояние от сотни метров до десятков километров. В то же время для обеспечения мер авиационной безопасности в приземном слое атмосферы, особенно в районе взлетно-посадочной полосы (ВПП) аэропорта необходимо проводить оперативный мониторинг ветровых характеристик атмосферы на расстоянии от единиц метров до нескольких сотен метров.
Для получения полной информации о скорости ветра необходимо проводить измерения, как минимум, при трех различных положениях зондирующего пучка в пространстве. Для этого используются различные методы сканирования атмосферы. В соответствии с известными методами КДЛ может работать в таких режимах сканирования, как RHI, PPI, VAD, либо их различная комбинация.
1. Range Height Indicator (RHI) - режим, при котором азимутальный угол остается постоянным, а угол места изменяется в определенном диапазоне.
2. Plan Position Indicator (PPI) - режим, при котором угол места остается фиксированным, как правило менее 10°, а азимутальный угол меняется либо во всем возможном диапазоне, либо в выбранном секторе.
3. Velocity Azimuth Display (VAD) - режим измерения производятся на нескольких азимутальных углах (обычно не менее 3) при фиксированном угле места (обычно 20° - 30° от вертикали).
В данной статье рассматриваются в первую очередь методы и алгоритмы обработки сигнала ветрового КДЛ с непрерывным излучением, в котором применяется коническое сканирование приземного слоя атмосферы.
2. Принцип конического сканирования
Как было сказано выше, для получения полной информации о трех компонентах вектора скорости ветра в декартовой системе координат с помощью КДЛ необходимо провести измерения, как минимум, при трех различных положениях зондирующего пучка в пространстве. В соответствии с рисунком 1, поясняющим принцип конического сканирования атмосферы, на заданной высоте h, определяемой необходимостью измерения, радиальную составляющую скорости ветра можно записать следующим образом:
,
где и - угол азимута лазерного луча, который изменяется в интервале [0; 2р), а ц - угол места лазерного луча, изменяющийся в диапазоне [0; р/2].
Таким образом, при измерении радиальной скорости ветра в одной точке записывается только одно уравнение при трех неизвестных. При измерении радиальной составляющей скорости ветра в трех независимых точках в пространстве на заданной высоте для трех неизвестных параметров записывается система из трех независимых уравнений. Решение данной системы уравнений позволяет определить полный вектор скорости ветра. При непрерывном вращении лазерного луча по углу азимута путем отбора его дискретных положений получается переопределенная система уравнений, которая решается методом наименьших квадратов.
Отметим, что измерение полного вектора скорости ветра при коническом сканировании воздушного пространства лазерным лучом на различных высотах достигается путем соответствующего изменения фокусного расстояния оптической системы КДЛ. Ключевым моментом принципа конического сканирования является предположение о “замороженности" ветровых характеристик атмосферы за время сканирования на заданной высоте. Поэтому в разработанном алгоритме измерение на одной высоте происходит в течение одной секунды.
Рисунок 1. Принцип конического сканирования атмосферы лазерным лучом.
В случае, когда ветровые характеристики атмосферы за период полного сканирования лазерным лучом по углу азимута на заданной высоте остаются постоянными, радиальная составляющая вектора скорости ветра изменяется по синусоидальному закону, как это показано на рисунке 2.
Рисунок 2. Закон изменения радиальной составляющей скорости ветра при постоянстве ветровых характеристик атмосферы при коническом сканировании.
Так как практически осуществляется измерение абсолютных значений доплеровского сдвига частоты, то и полученные результаты должны аппроксимироваться следующей функцией:
,
где a, b - аппроксимирующие коэффициенты.
Направление ветра по азимуту определяется как угол , соответствующий максимальной положительной скорости. Величина горизонтальной составляющей скорости определяется амплитудой гармонической функции. Вертикальным смещением полученного графика определяется величина вертикальной составляющей вектора скорости.
3. Методы и алгоритмы определения скорости и направления ветра
Основными задачами программно-алгоритмического обеспечения КДЛ являются:
1. Определение спектральных характеристик гетеродинного сигнала на выходе приемного устройства.
2. Обнаружение полезного сигнала в спектральных данных и формирование набора частот, характеризующего скоростные параметры воздушных масс.
3. Подготовка спектральных данных для аппроксимации путем удаления выбросов из набора частот с целью уменьшения ошибочных вычислений параметров ветра.
4. Проведение операции аппроксимации с целью определения искомых параметров функции приближения.
5. Окончательный расчёт параметров ветра.
Блок-схема алгоритма вычисления параметров ветра приведена на рисунке 3.
Рисунок 3. - Блок-схема алгоритма вычисления параметров ветра.
Алгоритм вычисления параметров ветра рассмотрим на примере обработки сигнала с профилометра лидарного ветрового ПЛВ-300 [3], разработанного ООО “Лазерные системы" в Санкт-Петербурге.
Сигнал с выхода приемника обратно рассеянного излучения содержит переменную составляющую с частотой до 50 МГц. Обрабатываемая полоса частот непосредственно задается диапазоном скоростей, определяемых лидаром. Аналого-цифровое преобразование осуществляется на частоте 100 МГц. Объем одиночной выборки содержит 512 отсчетов временного сигнала. Длительность выборки по времени - 5.12 мкс. Полученный временной сигнал поступает на блок быстрого преобразования Фурье (БПФ). На выходе блока БПФ вычисляется квадрат модуля комплексного Фурье-спектра.
Квадраты модулей спектров усредняются. Количество спектров в серии - 4096. Общая продолжительность одной серии съема данных - ~21 мс. На рисунке 4 представлена схема обработки сигнала.
Рисунок 4. - Схема обработки сигнала в единичном измерении.
На рисунке 5 показан вид результирующего спектрального сигнала.
Рисунок 5. - Результирующий спектральный сигнал.
Спектр содержит искомую дискретную компоненту в виде колоколообразного импульса. Отношение сигнал/шум находится в диапазоне 3 - 10 для нормальных атмосферных условий с МДВ в диапазоне 7 - 20 км. Ширина сигнала, определяемая по уровню 0.7 от максимального значения, составляет 1 - 5 МГц. Задача алгоритма заключается в поиске дискретной составляющей полезного спектра, определении ее максимального значения и вычислении соответствующей ему частоты доплеровского сдвига:
,
где N0 - положение, соответствующее нулевой скорости ветра.
После этого вычисляется абсолютное значение радиальной составляющей скорости ветра
,
где м - рабочая длина волны КДЛ.
Для выделения полезного сигнала несущего информацию о доплеровском сдвиге частоты, необходимо провести обработку спектра, усредненного за 4096 реализаций. Основные этапы обработки спектра представлены ниже:
1. Нормировка спектра S.
Данная операция позволяет выделить изолинию спектра. Характерный вид спектра с неравномерной фоновой характеристикой показан на рисунке 6. Неравномерность полученного спектра обусловлена шумовыми характеристиками лазера, приемника и трака АЦП.
a) С помощью медианного фильтра с фиксированным окном спектр сглаживается, из него удаляются резко выделяющиеся значения (в том числе полезный колоколообразный сигнал). На выходе фильтра получается спектр . На рисунке 7 показан сглаженный спектр, исходный вид которого приведен на рисунке 6.
b) Производится нормировка: .
Рисунок 6. - Усреднённый спектральный сигнал.
Рисунок 7. - Изолиния спектрального сигнала.
2. Определение частоты доплеровского сдвига в спектре полезного сигнала.
a) Выделение полезного колоколообразного сигнала на нормированном спектре. Для этого производится сглаживание “скользящим" или “бегущим” средним с маленьким окном, чтобы удалить случайные выбросы. В результате на выходе получается спектр , показанный на рисунке 8.
b) Расчет математического ожидания m и среднеквадратического отклонения s спектра .
Рисунок 8. - Спектральный сигнал после нормировки и сглаживания.
Для точек нормированного фильтрованного спектра, удовлетворяющих условию спектра >, рассчитывается центральная частота колоколообразного импульса , множитель k - задает порог чувствительности и лежит в пределах от 3 до 4, s - стандартное отклонение, m - математическое ожидание. Положение центральной частоты рассчитывается в соответствии с известным центроидным методом по формуле центра масс:
,
где - частота i-й точки спектра,
- масса i-й точки (значение по амплитуде).
Кроме приведенного выше центроидного метода могут быть использованы такие методы, как метод гауссовой аппроксимации [4] либо метод с использованием порядковых статистик [5].
На рисунке 9 показан уровень максимума спектрального сигнала и рассчитанный уровень порога .
3. Пункты 1 и 2 повторяются для окна медианного фильтра шириной в 2 раза большей. Форма сигнала (ширина колоколообразного импульса) зависит от текущих метеоусловий и может меняться. Алгоритм повторяется с большим окном медианного фильтра, чтобы не пропустить импульсы полезного сигнала различной ширины. Это позволяет использовать алгоритм независимо от текущей формы сигнала.