Восстановление вектора скорости ветра по данным лидарного зондирования методом сплайн-аппроксимации

Математические технологии

УДК 519.6

DOI:10.25729/ESI.2023.30.2.009

Восстановление вектора скорости ветра по данным лидарного зондирования методом сплайн-аппроксимации Баранов Николай Алексеевич

Вычислительный центр им. А.А. Дородницына ФИЦ ИУ РАН, Россия, Москва, baranov@ians.aero

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

Ключевые слова: восстановление вектора скорости ветра, сплайн-аппроксимация, лидарное сканирование

Цитирование: Баранов Н.А. Восстановление вектора скорости ветра по данным лидарного зондирования методом сплайн-аппроксимации / Н.А. Баранов // Информационные и математические технологии в науке и управлении. - 2023. - № 2(30). - С. 83-97. - DOI: 10.25729/ESI.2023.30.2.009.

Введение. Наземные лидарные системы дистанционного измерения скорости ветра находят все большее применение. В частности, они применяются в аэропортах для обнаружения сдвига ветра, а также мониторинга вихревых следов за самолетами [1]. Особенностью измерений ветровой обстановки лидарными системами является то, что они обеспечивают измерение проекции скорости ветра на направление сканирования. Пространственная картина ветрового поля получается за счет кругового сканирования с некоторым постоянным углом места, так называемый режим сканирования PPI - Plan Position Indicator. При сканировании с малым углом места можно считать, что вклад вертикальной составляющей ветра в измерения пренебрежимо мал и измеряется только горизонтальная составляющая скорости.

Поскольку непосредственно измеряется только радиальная составляющая скорости ветра (проекция на направление измерения) возникает задача восстановления двумерного вектора скорости по данным измерений в нескольких направлениях.

В работе [2] представлен краткий обзор методов, применяемых для решения задачи восстановления вектора скорости горизонтального ветра по данным лидарного сканирования в режиме PPI. В частности, рассмотрены вариационные методы восстановления и метод оптимальной интерполяции.

В методе оптимальной интерполяции [3] компоненты скорости ветра определяются путем минимизации функционала вида

J = (w - wb Б-1 (w — wb) + (wG - H • w R—1 (w0 — H • w),

где w - вектор искомых параметров скорости ветра; wb - вектор априорной информации о скорости ветра в анализируемом пространстве; wG - вектор измерений скорости ветра; Б и R - соответственно ковариационные матрицы скорости ветра и ошибки измерений; H - матрица наблюдений. В методах вариационного типа в функционал качества добавляются еще функции штрафа, учитывающие ограничения на дивергенцию или завихренность векторного поля [4 - 6].

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

Для восстановления вектора скорости ветра по данным лидарного сканирования также активно используется метод наименьших квадратов (МНК) с весовыми коэффициентами [8 -10], согласно которому вычисление горизонтальных компонент и, V скорости ветра в точке г

с координатами (р0, г0) по измерениям V в точках г = (р,г ) из некоторой ее окрестности

основано на минимизации функционала вида

I

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

Проблема применения МНК связана с тем, что из-за наличия пропусков в данных измерений оценки компонент скорости ветра могут быть смещенными, что приводит к выбросам в решении. Эту проблему иллюстрирует пример, показанный на рис. 1. Рассматривается постоянное поле ветра с компонентами скорости и0 = 5 м/с, V = 6 м/с. Вектор, соответствующий истинному направлению и скорости ветра, показан на рис. 1 зеленым цветом.

Измерения представляют собой значение радиальной составляющей скорости в точках г, равномерно расположенных вдоль направлений измерений р. :

Щ] = ио • совР. + V • ,

где - ошибка измерений, которая является нормально распределенной случайной величиной с математическим ожиданием, равным нулю, и среднеквадратическим отклонением 0.25м/с. Предполагается, что некоторые измерения являются недостоверными - они показаны на рис. 1 выколотыми точками. Для вычисления скорости в точке р0 = 0°, г0 = 2200 м отбираются измерения в точках, удаленных от точки (р0, г0 ) не более, чем на Аг = 200 м. Соответствующая окрестность точки г показана пунктирной линией.

Можно видеть (рис. 1), что в окрестность |гу. - г0| < Аг попадает некоторое множество

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

Результатом решения задачи восстановления методом наименьших квадратов является вектор скорости ветра, показанный на рис. 1 красным цветом. При этом погрешность определения скорости составила »2м/с, что составляет 25% от истинного значения, а погрешность определения направления »20°.

Представленный пример показывает, что МНК довольно чувствителен к пропускам данных.

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

Рис. 1. Пример восстановления вектора скорости ветра методом наименьших квадратов

при наличии недостоверных данных В данной работе рассматривается возможность применения альтернативного подхода к восстановлению двумерного вектора горизонтальной скорости ветра по данным лидарного сканирования в режиме РР1, при котором компоненты скорости ветра вычисляются не в отдельной точке, а в заданном наборе узлов (г,...,гп) вдоль заданного направления (0. Предлагаемый алгоритм использует идею покоординатного вычисления компонент вектора: вначале вычисляется компонента скорости вдоль направления (0, а затем - ортогональная ей составляющая скорости ветра. Каждая компонента вычисляется путем сплайн-аппроксимации специальным образом сформированного набора данных измерений.

1. Постановка задачи. Будем предполагать, что измерения выполняются вдоль направ-

лений

( = (о + 7 Л( 7 = -п,..., П,

в контрольных точках, расположенных с постоянным шагом:

Г = г0 + г Лг, г = 0,..., т.

В результате имеется множество измерений V проекций скорости ветра на направления измерений ( на дальностях г. Для каждого измерения также определен признак его достоверности ^ :

Бу = 1, если измерений достоверно, £г = 0, в противном случае.

В качестве показателя достоверности может использоваться, например, соотношение сигнал-шум SNR, соответствующее измерению V [1, 11]. Если значение SNR превышает

некоторое заданное пороговое значение SNRlim:

SNRg> SNRm

то измерение считается достоверным S!} = 1. В противном случае измерение недостоверно и S„ = 0.

Расчет выполняется в повернутой декартовой системе координат, ось абсцисс которой совпадает с направлением % . Вектор скорости ветра имеет компоненты (и, v), где и - компонента скорости ветра вдоль направления %, v - компонента, ортогональная направлению %

Измерения V выражаются через компоненты скорости ветра (utj, Vy ) в этой точке зависимостью вида

V = и..  cos ( j А%) + Vу  sin ( j A%) + S tj, где S - погрешность измерений. Предполагается, что ошибки измерений независимы.

Предполагается, что поле ветра является плоскопараллельным, т.е. компоненты скорости являются некоторыми одномерными функциями дальности вдоль направления %, определенными на интервале [r0, rm ], и не зависят от азимутальной координаты. Гипотеза о пло-скопараллельности поля ветра означает, что скорость ветра в точке r. вдоль направления % равна скорости в точке г  cos (j А%) вдоль направления % :

и (Г ) = и (Г cos (j А%) %о), v(г , %) = v (Г cos (j А%) , %) •

Другими словами, согласно гипотезы о плоскопараллельности поля ветра компоненты скорости ветра являются функциями только координаты x (в декартовой системе координат, ось абсцисс которой совпадает с направлением % ) и не зависят от ординаты (рис. 2). Соответственно значение радиальной скорости V в точке r вдоль направления % является проекцией скорости ветра в точке г  cos (j А%) вдоль направления % на направление % .

Заметим, что при восстановлении скорости ветра для другого направления %% Ф % по

данным измерений радиальной скорости с учетом этой гипотезы будет приниматься допущение, что компоненты скорости ветра являются функциями только координаты x в декартовой системе координат, ось абсцисс которой совпадает с направлением %01, а не % .

Требуется найти оценки (u0 (r), V0 (r)) скорости ветра вдоль направления %, которые минимизируют функционал качества вида

m n 2

J

i0''i 0() = ZZ Sij (Vj u0 ( r - cos ( j A(ll- cos 1 j A( I > r - cos 1 j A(/|- Sin 1 j A()) +

i=0 j=-n

+(flV- <2>

r

Здесь Л > 0 - параметр сглаживания, контролирующий соотношение между точностью аппроксимации данных и гладкостью аппроксимирующей функции.

Заметим также, что хотя точки измерений вдоль каждого направления расположены равномерно с одинаковым шагом, аппроксимация компонент скорости ветра осуществляется на неравномерной сетке узлов {r • cos(j А^) j . Кроме того, часть узлов сетки может быть исключена из рассмотрения в силу недостоверности измерений, что определяется весовым коэффициентом .

Рис. 2. Геометрическая иллюстрация гипотезы плоскопараллельного поля ветра

2. Алгоритм вычислений. Для решения задачи определения компонент скорости ветра вдоль направления % будем использовать модифицированный функционал (2), учитывая

только измерения, находящиеся в малой 5-окрестности направлений %, % ± пА%, т.е. измерения вдоль направлений, отличающихся от направлений %, % ± пА% не более, чем на 5 (рис.

3):

%.-п -А%<5,

-% + п -А% < 5. При таком ограничении функционал (1) будет иметь вид

т 2

1

1=0 г]] А%<5, (п -)А%<5

+

1 р) #)

Выполним декомпозицию задачи, пользуясь тем, что проекция боковой составляющей скорости ветра на направление % равна нулю. Запишем функционал (3) в виде

1 ({^0> Яо }) * Л ({Яо }) +12 ({Яо, Яо })>

где

J1 ({йгo j) = E Z Sj (V - Uo (r • cos (j аф)У cos (j А^}) +Л|

i=0 j:\jAriiS

du0

dr

dr,

2

r

m

r

0

J2

i=0 j: (n-|j|)A(<S

+4 i f !2 dr

Рис. 3. Геометрическая иллюстрация отбора направлений измерений для формирования

модифицированного функционала В силу того, что параметр 5 мал, функционал ^ (|й0}) представим в виде

mm

Ji'№U)«Z Z Sj(Vj-щ(r)) +4

du

dr

dr.

1 =0 ]': ]Лф5 _

Особенностью функционала (4) является то, что в узлах г заданы несколько значений аппроксимируемой функции й0 (г): {у.. }.|.д | д с учетом признака достоверности измерений .

у

Поэтому задачу минимизации функционала (3) запишем в виде

min i Z sv  Vi - "о +д J

du0 dr

dr !

mm

j:| j A(<S

Коэффициент max Sr учитывает возможность, что на дистанции r могут отсутство-

вать достоверные измерения. Формулировка задачи построения функции u0 (r), аппроксимирующей измерения Vij на интервале [r0, rm] в диапазоне направлений измерений |% -%| < S в виде (5), (6) позволяет выполнить ее декомпозицию на две подзадачи:

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

- построение аппроксимирующей функции й0 (г) по заданным значениям Vi в узлах ri как решение задачи (5).

Для узлов г , для которых имеется хотя бы одно достоверное значение

Z s. >0,

j: I j A%<S

2

, г = о,...,т.

задача вида (6) имеет очевидное решение вида

у _ ]-\1 А%<5_

,= I *

] А%<5

Решение задачи (5) будем искать в классе кубических сплайнов, удовлетворяющих естественным граничным условиям [12]

ио'(Го) = иоЧГт) = о . В этом случае на каждом интервале [г , гм ] функция и0 (г) имеет вид

К2

и

( г) = ио. ( г) = (1 - 0 + - (1 - 0 [( 2 - +(1 + ^)пг+1 ]

где

К = Г-+1- Г, * =

г - г

К

Коэффициенты щ вычисляются путем решения системы уравнений [13, 14]

аопо + Ьоп1 + Соп2 = ёо ,

Ъ0п0+ап + Ьп+сп = &,

С-2п.-2 + Ъ-п- + ап + ЪЛ+1 + сЛ+2 = ё., I = 2,..., т - 2,

Ст-3пт-3 + Ът-2пт-2 + т-1пт-1 + Ът-1Ит = ёт-1 ,

с ,п , + Ъ ,п , + а п = ё ,

т-2 т-2 т-1 т-1 т т от'

где

а =1 (К-1 + К Р- +

Ь =1К -1

' 6 ' К

(Г

г-1

1 1

— + —

КК

V К-1 К /

Р + 7ГР+1, * = 1,...,т-1

Л ^ Р +

VV г-1 г/

1

К+К

1

К

л л

р+1

V "г "г +1 / /

, / = 2,..., т - 2,

с =

г КК

' Р+1, г' = 1,..., т - 3,

г+1

77 К-К

/' = !,...,1И-1,

ао = 1, Ъо = 1, со = ° ёо = о,

а = 1, Ъ ,= 1, с , = о, ё = о.

т 7 т-1 7 т-2 7 от

Здесь р > о - весовые коэффициенты. Представленная система уравнений для нахождения коэффициентов п имеет симметричную, положительно определенную матрицу с диагональным преобладанием. Известно [12, 13], что решение такой системы может быть осуществлено методом исключения Гаусса без выбора главных элементов, т. е. в данном случае хорошо апробированным методом пятиточечной прогонки [13 - 15]. Хорошая обусловленность матрицы системы уравнений и наличие высокоэффективных процедур решения этой системы и определили выбор метода построения аппроксимирующей функции с использованием кубических сплайнов. Размерность матрицы варьируется в зависимости от дистанции, на которой

осуществляется восстановление вектора скорости. Так, например, для дальности 2000м размерность матрицы составит 200*200.

При известных значениях коэффициентов щ коэффициенты zi вычисляются по формулам

zi=Vi~PiDi, i = l...,m,

где

D0 = — (Щ1 " Щ0 ) , Dm = —-( Пт " Щт-1 ) ,

h0 K-1

Di = У (Щ+1 -ЩЬт1"(Щ -Щ-1)> i = 1,...,m-1.

Таким образом, функция u0 (r) будет однозначным образом определена. После этого осуществляется восстановление боковой составляющей скорости ветра v0 (r) . Для ее вычисления преобразуем функционал J ({йг0, vi0}), учитывая малость параметра 5, к виду

J2({V,}) = Z Z S. Sin2 ( j A(

i=о j: (n-\ j|)A(<S

rm •cosf n A

f

V - u

sin

- V,

2

^ Í dr.

n A

dr

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

ветра на подзадачи:

min

S^^/.-^^o^-cos^A^)) нЧ J

-=о V 1 r0 •cos

rm•cosfnA() / \2 ^

dvо

dr

dr

min Z Sijsin2

j: (n-| j\)A(<S

V -u (r • cos (n A())' cos ( j A((

-Ж

sin^A^) y

Задача (В) имеет решение вида

£ sv sin UA(P) {vy - ûo (гг cos (n A(P)) cos (7 A(P)}

,i = 0,...,m.

Ж =

j: (n-| j|)A(<S

]': (и-|]\)Ар<5

а задача (7) решается аналогично задаче (4) для продольной составляющей.

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

й (X, у ) = йо +Ай  х,

у (X У ) = V) + \ X

где й0 = 5 м/с, V = 6 м/с, Ам = _0.004, Ау = 0.004. Будем предполагать, что ошибки измерений являются независимыми, нормально распределенными случайными величинами с математическим ожиданием 0 и среднеквадратической ошибкой 0.25м/с.

2

На рис. 4 представлены результаты расчета модуля скорости ветра для направления % — 0° при различных значениях параметра Л, регулирующего степень близости аппроксимирующей функции к заданным значениям в узлах. Можно видеть, что при значениях параметра Л~102...103 ошибка восстановления скорости составляется 10...15%, что является довольно высоким показателем. При значениях Л ~ 105 метод теряет чувствительность к вариации восстанавливаемой функции на области определения, фактически вырождаясь в линейную регрессию. Поэтому для дальнейших расчетов было принято значение Л —104.

131

и «

э-

Щ ™ 11

8 о

—•— Л-10е-02 —•— л=1.0е-03

—•— Л—I Ое • эталон -05

250

500

750

1000 1250

Дистанция, м

1500

1750

2000

Рис. 4. Влияние параметра сглаживания Л на восстановление модуля скорости ветра Также на этом модельном примере было рассмотрено влияние размера сектора данных, используемых для восстановления скорости вдоль заданного направления % — 5°, т.е. значения параметра пА%. Соответствующие результаты представлены на рис. 5. Можно видеть, что при малом размере сектора методу не хватает чувствительности (график зеленого цвета) в силу сглаживающих свойств аппроксимирующего сплайна, поскольку величина проекции боковой составляющей скорости ветра на направление пА% — 5° оказывается мала. При большом размере сектора пА% — 30° на качество восстановления начинает влиять пространственная изменчивость ветрового поля. Таким образом, целесообразно размер сектора данных выбирать в диапазоне 10°...20°.

250 500 750 1000 1250 1500 1750 2000

Дистанция, м

Рис. 5. Влияние размера пА% сектора данных на восстановление модуля скорости ветра

Апробацию представленного алгоритма расчета выполним на примере данных реальных измерений ветровых полей с помощью импульсного ветрового лидара Windex-5000 [16]. Ли-дар обеспечивает измерение радиальной составляющей скорости ветра с пространственным шагом Ar « 41 м. Измерения выполнялись с угловым разрешением Аф = 1° . Вид поля измеренных радиальных скоростей представлен на рис. 6. На рисунке показаны достоверные данные измерений.

-2000 -1500 -1000 -500 О 500 1000 1500 2000

Рис. 6. Данные измерений ветрового поля На рис. 7, 8 показаны результаты восстановления вектора скорости ветра для нескольких направлений измерений % е {70°, 80°, 90°, 100°, 110°} . Результаты представлены в виде значений модуля вектора скорости (рис. 7) и направления скорости ветра (рис. 8). Начальное направление % — 70° показано на рис. 2 линией черного цвета, что позволяет судить об изменчивости данных сканирования в анализируемой области. Для восстановления поля скоростей рассматривался сектор размером 20°, т.е. п=10. Величина 5 была принята равной 2°. Дистанция восстановления вектора скорости ветра составила 2000м.

250 500 750 1000 1250 1500 1750 2000

Дистанция, м

Рис. 7. Результаты восстановления скорости ветра в диапазоне направлений 70.110°

Для сравнения на этих же рисунках пунктирными линиями показаны результаты, полученные методом наименьших квадратов (1). На рис. 9, 10 аналогичные результаты представлены для направлений измерений % е {20°, 30°, 40°, 50°, 60°}.

1000 1250

Дистанция, м

Рис. 8. Результаты восстановления направления ветра в диапазоне направлений 70.110°

1000 1250

Дистанция, м

Рис. 9. Результаты восстановления скорости ветра в диапазоне направлений 20.60° Можно видеть, что применение описанного в данной работе алгоритма обеспечивает более гладкое восстановление характеристик скорости ветра по сравнению с методом наименьших квадратов. При применении предложенного алгоритма отсутствует резкие перепады параметров скорости ветра на короткой дистанции.

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

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

Рис. 10. Результаты восстановления направления ветра в диапазоне направлений 20.. .60°

При этом следует отметить, что на результаты восстановления существенным образом влияет размер сектора данных, используемых для решения задачи. Этот факт иллюстрируют рис. 11, 12, на которых представлены результаты расчета скорости и направления ветра при различных размерах сектора используемых данных: 10° (круглые маркеры), 20° (ромбовидные маркеры), 30° (треугольные маркеры).

Рис. 11. Влияние размера сектора измерений на результаты восстановления скорости ветра

Рис. 12. Влияние размера сектора измерений на результаты восстановления направления

ветра

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

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

Список источников

1. ISO 28902-2:2017. Air quality - Environmental meteorology - Part 2: Ground-based remote sensing of wind by heterodyne pulsed Doppler lidar.

2. Liu, Z., Barlow J.F., Chan P.-W., Fung J.C.H., Li Y., Ren C., Mak H.W.L., Ng E. A review of progress and applications of pulsed doppler Wind LiDARs. Remote Sens., 2019, 11, 2522.

3. Choukulkar A., Calhoun R., Billings B., Doyle J.D. A modified optimal interpolation technique for vector retrieval for coherent doppler LIDAR. IEEE Geosci. Remote Sens. Lett., 2012, v. 9, pp. 1132-1136.

4. Cherukuru N.W., Calhoun R., Krishnamurthy R., Benny S., Reuder J., Flügge M. 2D VAR single Doppler lidar vector retrieval and its application in offshore wind energy. Energy Procedia, 2017, v. 137, pp. 497-504.

5. Qiu C.J., Shao A.M., Liu S. Xu Q. A two-step variational method for three-dimensional wind retrieval from single Doppler radar. Meteorol. Atmos. Phys., 2006, v. 91, pp. 1-8.

6. Xu Q., Gong J. Background error covariance functions for Doppler radial-wind analysis. Q.J.R. Meteorol. Soc., 2003, v. 129, pp. 1703-1720.

7. Menke R., Vasiljevic N., Wagner J., Oncley S. P., Mann J. Multi-lidar wind resource mapping in complex terrain, Wind Energ. Sci., 2020, 5, 1059-1073.

8. Baranov N. Algorithms of 3D wind field reconstructing by lidar remote sensing data, numerical computations: Theory and Algorithms. NUMTA 2019. Lecture notes in Computer Science, 2020, v. 11974. Springer, Cham.

9. Newsom R.K., Brewer W.A., Wilczak J.M., Wolfe D.E., Oncley S.P., Lundquist J.K. Validating precision estimates in horizontal wind measurements from a Doppler lidar, Atmos. Meas. Tech., 2017, 10, 1229-1240.

10. Bodini N., Zardi D., Lundquist J.K.: Three-dimensional structure of wind turbine wakes as measured by scanning lidar, Atmos. Meas. Tech., 2017, 10, 2881-2896.

11. Banakh V.A., Smalikho I.N. Lidar studies of wind turbulence in the stable atmospheric boundary layer. Remote sens. 2018, 10, 1219.

12. Завьялов Ю.С. Методы сплайн-функций. / Ю.С. Завьялов, Б.И. Квасов, В.Л. Мирошниченко. - М.: Наука, 1980.

13. Квасов Б.И. Численные методы анализа и линейной алгебры / Б.И. Квасов. - Новосибирск, Новосиб. гос. ун-т, 2012.

14. Квасов Б.И. Методы изогеометрической аппроксимации сплайнами / Б.И. Квасов. - М.: Физматлит, 2006.

15. Самарский А.А. Методы решения сеточных уравнений / А.А. Самарский, Е.С. Николаев. - М.: Наука, 1978.

16. Kim A.A., Orlov A.V., Luginya V.S., Baranov N.A. About the methodology and tools for meteo LIDAR metro-logical support, Proc. SPIE, 2020, vol. 11560, 26th International symposium on atmospheric and ocean optics, Atmospheric physics, 1156022.

Баранов Николай Алексеевич. Доктор технических наук, главный научный сотрудник Вычислительного центра им. А.А. Дородницына Федерального исследовательского центра «Информатика и управление» Российской академии наук, AuthorID; 14242, SPIN: 5876-2490, ORCID: 0000-0002-6419-2051, baranov@ians.aero, 119333, Москва, ул. Вавилова, 40.

UDC 519.6

DOI:10.25729/ESI.2023.30.2.009

Recovering of the wind speed vector from lidar sensing data using the spline

approximation

Nikolay A. Baranov

Dorodnicyn Computing Centre, FRC CSC RAS, Russia, Moscow, baranov@ians.aero

Abstract. The paper considers the problem of recovering the horizontal wind speed from lidar scanning data, which provide measurements of the radial component. An approach is proposed in which the wind speed components are calculated not at a single point, but at a given set of nodes located along a given direction. The proposed algorithm uses the idea of coordinate-wise calculation of vector components. Each component is calculated by the spline approximation of a specially selected subset of measurements. The results of numerical calculations are presented.

Keywords: wind speed recovering, spline approximation, lidar scanning

References

1. ISO 28902-2:2017. Air quality - Environmental meteorology - Part 2: Ground-based remote sensing of wind by heterodyne pulsed Doppler lidar.

2. Liu, Z., Barlow J.F., Chan P.-W., Fung J.C.H., Li Y., Ren C., Mak H.W.L., Ng E. A review of progress and applications of pulsed doppler Wind LiDARs. Remote Sens. 2019, 11, 2522.

3. Choukulkar A., Calhoun R., Billings B., Doyle J.D. A modified optimal interpolation technique for vector retrieval for coherent doppler LIDAR. IEEE Geosci. Remote Sens. Lett., 2012, v. 9, pp. 1132-1136.

4. Cherukuru N.W., Calhoun R., Krishnamurthy R., Benny S., Reuder J., Flügge M. 2D VAR single Doppler lidar vector retrieval and its application in offshore wind energy. Energy Procedia, 2017, v. 137, pp. 497-504.

5. Qiu C.J., Shao A.M., Liu S. Xu Q. A two-step variational method for three-dimensional wind retrieval from single Doppler radar. Meteorol. Atmos. Phys. 2006, v. 91, pp. 1-8.

6. Xu Q., Gong J. Background error covariance functions for Doppler radial-wind analysis. Q.J.R. Meteorol. Soc. 2003, v. 129, pp. 1703-1720.

7. Menke R., Vasiljevic N., Wagner J., Oncley S. P., Mann J. Multi-lidar wind resource mapping in complex terrain, Wind Energ. Sci., 2020, 5, 1059-1073.

8. Baranov N. Algorithms of 3D wind field reconstructing by lidar remote sensing data, numerical computations: Theory and Algorithms. NUMTA 2019. Lecture notes in Computer Science, 2020, v. 11974. Springer, Cham.

9. Newsom R.K., Brewer W.A., Wilczak J.M., Wolfe D.E., Oncley S.P., Lundquist J.K. Validating precision estimates in horizontal wind measurements from a Doppler lidar, Atmos. Meas. Tech., 2017, 10, 1229-1240.

10. Bodini N., Zardi D., Lundquist J.K.: Three-dimensional structure of wind turbine wakes as measured by scanning lidar, Atmos. Meas. Tech., 2017, 10, 2881-2896.

11. Banakh V.A., Smalikho I.N. Lidar studies of wind turbulence in the stable atmospheric boundary layer. Remote sens. 2018, 10, 1219.

12. Zavyalov Yu.S., Kvasov B.I., Miroshnichenko V.L. Metody splayn-funktsiy [Methods of spline functions]. M. Nauka [Moscow. Scince], 1981.

13. Kvasov B.I. Chislennyye metody analiza i lineynoy algebry [Numerical methods of analysis and linear algebra]. Novosibirsk, Novosib. gos. un-t [Novosibirsk-Novosib. State university], 2012.

14. Kvasov B.I. izogeometricheskoy approksimatsii splaynami [Methods of isogeometric approximation by splines]. M. Fizmatlit [Moscow.Fizmatlit], 2006.

15. Samarsky A.A., Nikolaev E.S. Metody resheniya setochnykh uravneniy [Methods for solving grid equations]. M. Nauka [Moscow. Scince], 1978.

16. Kim A.A., Orlov A.V., Luginya V.S., Baranov N.A. About the methodology and tools for meteo LIDAR metro-logical support, Proc. SPIE, 2020, vol. 11560, 26th International symposium on atmospheric and ocean optics, atmospheric physics, 1156022.

Baranov Nikolay Alekseevich. Doctor of Science, Dorodnicyn Computing Centre, Federal Research Center "Computer Science and Control" of Russian Academy of Sciences, Chief Researcher, AuthorID: 14242, SPIN: 58762490, ORCID: 0000-0002-6419-2051, baranov@ians.aero, 119333, 40, Vavilov str., Moscow, Russia.

Статья поступила в редакцию 11.05.2023; одобрена после рецензирования 13.06.2023; принята к публикации 16.06.2023.

The article was submitted 05/11/2023; approved after reviewing 06/13/2023; accepted for publication 06/16/2023.