УДК 621.396.41
Е.В. Масалов, Н.Н. Кривин, К.В. Кокоулин
Анализ изменчивости линейного деполяризационного отношения при радиолокационном зондировании неоднородной среды, заполненной гидрометеорами
Получены расчетные соотношения для оценки величины линейного деполяризационного отношения в точке приема радиолокационных сигналов. Это сделано с использованием компонент вектора Джонса собственного базиса однородного участка среды распространения и комплексного фазора рассеянного сигнала. При этом имеется в виду, что данный сигнал рассеивается вторым участком, который следует непосредственно за однородным участком; характеризуется анизотропными по поляризации свойствами и имеет угол ориентации собственного поляризационного базиса, отличный от угла ориентации собственного базиса однородного участка среды распространения. При этом анализируются изменения поляризационной структуры радиолокационного сигнала в процессе его распространения в указанной среде в зависимости от начального угла ориентации поляризационного базиса РЛС относительно собственного базиса однородного участка, величины степени анизотропии второго участка, угла ориентации собственного поляризационного базиса второго участка. Приводятся результаты расчетов величины линейного деполяризационного отношения в точке приема радиолокационного сигнала для осадков в виде дождя различной интенсивности. Полученные результаты анализируются. Ключевые слова: поляризация, угол эллиптичности, угол ориентации эллипса поляризации, дифференциальное ослабление, дифференциальный фазовый сдвиг, угол ориентации поляризационного базиса однородного участка, степень поляризационной анизотропии, угол ориентации собственного поляризационного базиса анизотропного участка, линейное деполяризационное отношение. doi: 10.21293/1818-0442-2018-21-3-7-13
Под неоднородной средой распространения понимается воздушное пространство, заполненное частицами осадков: каплями дождя, частицами града, снежинками и состоящее из участка однородной среды и следующего за ним второго участка, отличающегося анизотропными свойствами [1, 2]. Под однородным участком среды следует понимать участок, заполненный частицами метеообразований, имеющими одинаковую ориентацию и размеры, находящимися в одном и том же фазовом состоянии [1, 3, 4]. Базис такого участка среды является линейным [4].
В случае дистанционного зондирования метеообразований с использованием метеорадиолокатора поляризационная структура радиолокационных сигналов оказывается в значительной степени подвержена искажающему воздействию дифференциальных факторов среды (дифференциального ослабления Да и дифференциального фазового сдвига АФ) [3, 4]. Это обстоятельство существенно снижает эффективность идентификации анизотропного участка метеообразования, расположенного на периферии от передней границы метеообразования.
Кроме того, то обстоятельство, что угол ориентации собственного поляризационного базиса анизотропного участка среды, как правило, является отличным от ориентации базиса однородного участка среды, обусловливает многофакторность воздействия метеообразований на поляризационную структуру радиолокационного сигнала. Указанные отличия угла ориентации базиса 9 могут быть обусловлены как наличием поперечных (к трассе распространения) сдвигов ветра, так и наличием зон повышенного статического электричества (зон грозовой активности) [5].
Постановка задачи
Результаты оценки интенсивности осадков и алгоритмы распознавания опасных зон метеообразований на основе поляризационных измерений с применением двух ортогональных линейных поляризаций [6, 7] и использовании существующих методик не позволяют учесть указанные выше факторы и вызывают затруднение при идентификации опасных участков метеообразований (сдвиги ветра, зоны грозовой активности). Кроме того, весьма существенным оказывается многофакторное воздействие таких метеообразований на поляризационную структуру радиолокационного сигнала. Указанное воздействие влияет:
- на величину угла рн между измерительным базисом, связанным с РЛС, и собственным базисом первого (однородного) участка среды;
- величину угла 9 между ориентацией поляризационного базиса второго (анизотропного по поляризации) участка и ориентацией первого (однородного) участка среды;
- степень поляризационной анизотропии второго участка среды д.
Эти обстоятельства обусловливают необходимость поиска способов, позволяющих оценить значимость влияния указанных факторов на получаемые результаты измерения и определить характер изменчивости параметра ЬБК^), где г - длина трассы распространения радиолокационного сигнала в однородном участке среды, заполненной гидрометеорами.
Методика решения
Как известно, величина ЬБК^) определяется соотношением [8-11]
ЬБЯ(г) = 2О• (г)/\Evxx(г)), (1)
где [Еру (г)| - амплитуда сигнала, излученного с
горизонтальной поляризацией и принятого с вертикальной; |Е?рхх (г) - амплитуда сигнала, излученного
и принятого с горизонтальной поляризацией.
Для нахождения указанных амплитуд целесообразно использовать подход, описанный в [4]. Данный подход основан на использовании методики определения компонент вектора Джонса [12] в собственном базисе среды распространения с последующим переходом в собственный базис РЛС. При этом анализируемые сигналы определяются на выходе прямоугольных волноводов, принимающих сигналы горизонтальной и вертикальной поляризации соответственно. При этом согласно результатам работы [4] определение угла ориентации измерительного поляризационного базиса и определение угла ориентации собственного поляризационного базиса анизотропного участка целесообразно производить относительно поляризационного базиса однородного участка среды распространения, а оценку трансформации поляризационной структуры радиолокационного сигнала целесообразно проводить в собственном базисе однородного участка среды распространения.
Тогда комплексные амплитуды, необходимые для определения ЬВК(г) в точке приема (расположение РЛС), могут быть определены из приведенных ниже выражений:
"0 ОЦя-МИЗр(г,Ри
Ерхх(г) =
Ерху (г) =
О О О 1
[я_1(Рн)][(Рр (г,Ри))] [я_1(Рн)][(Рр (г,Ри))]
(2)
соБар( г) '1пар (г)
. (3)
В формулах (2) и (3) использованы следующие обозначения:
соБар (г) ' бш аР (г)
- вектор Джонса эллиптически по-
ляризованной волны в базисе анизотропного участка;
[[ (г,Рн ))] =
собРр (г ,Рн) Б1иРр (г,вн) -Б1ПРр (г ,Рн ) СОБРр (¿,Рн )
оператор перехода из базиса анизотропного участка в базис однородного участка среды;
^_1(Рн)] = [СОпрРН -ОЖ?] - оператор перехода в базис локатора из базиса однородного участка среды; [1 О
О О
- матрица Джонса прямоугольного вол-
новода, принимающего сигнал с горизонтальной поляризацией; "О О"
О 1
- матрица Джонса прямоугольного вол-
новода, принимающего сигнал с вертикальной поляризацией;
аР(г) и рР(г) представляют собой значения угла ориентации эллипса поляризации и угла эллиптич-
ности эллипса поляризации соответственно, изменяющиеся в процессе распространения сигнала, рассеянного анизотропным участком, в однородном участке среды.
Для учета указанных выше факторов воспользуемся моделью среды распространения, когда зондирующий сигнал распространяется в однородной среде, а его рассеяние происходит анизотропным участком, заполненным гидрометеорами. В рамках предложенной модели неоднородной среды, для характеристики анизотропного участка метеообразования, целесообразно использовать матрицу рассеяния следующего вида [13]:
5=О,5( +Х 2){
соб29 б1и29 б1и29 -СОБ29
(4)
где Я,1, Х2 - собственные числа матрицы рассеяния, Я,1 > ^2; 9 - угол ориентации собственного базиса метеообъекта относительно измерительного;
степень поляризационной анизотро-
Х1 -Х2 Х1 +Х2
пии объекта.
Поскольку, как отмечается в работе [14], метео-бразования являются объемными целями, образованными независимыми рассеивателями, то матрицу рассеяния (4), степень анизотропии д, собственные значения матрицы Я,1, ^2 и угол ориентации 9 собственного базиса анизотропного участка следует понимать как совокупную матрицу рассеивания объемной метеоцели и ее характеристики (д, Я,1, Я,2 , 9).
Выражения, определяющие аР(г) и рР(г), будут
иметь следующий вид:
ар(г) =1 х
хагсБШ
1О°,°5Даг • урН (г) Бт(ФрН (г) + ДФг)
хarctg
1+1О°,1-Даг • tg2 урн (г) Рр(г)=1 х
2 • Ш°,О5-Даг • tgурн (г) • соб(фрн (г) + ДФг)
1 -юОД Даг • tg2 урн (г)
(6)
где tgурн(г) - начальное значение модуля комплексного фазора рассеянной волны в зависимости от длины трассы распространения г; фрн (г) - начальное значение фазы комплексного фазора рассеянной волны.
начальное значение комплексного фазора определяется с учетом матрицы рассеяния, представленной выше в виде [4]:
Шурн (г) = |Еурн ( / |Ехрн (г)| = = 512 + 522 \Ёу1 (г)/ЕХ1 (г))
(5)
пп
5П + 512 •(((г)/ЁХ1 (г))
(7)
Поскольку, как было указано выше, оценку трансформации поляризации сигнала, облучающего анизотропный участок, целесообразно проводить в собственном базисе однородного участка, то комплексный фазор ЁЁу\ (г)/ Ёх\ (г) облучающей этот
участок волны может быть представлен согласно результатам работы [4] в виде
где
1^ = tg Yl( z ) • e-^ф( z ), ÈÈxl( z )
tgYl(z) = 100,05•Дaz • tgß^
(8)
(9)
а изменение фазы фазора ф(z) в процессе распространения составляет величину
ф(z) = ДФ- z . (10)
Подставив значения Sj из формулы (4) и выражения (9) и (10) в формулу (7) и проведя необходимые преобразования, получим, что
, ч r2|a-sin28-(l-|a-cos28)-tgy1 (z)cosДФz +
yPH(z)= arctg-------f--—-
2|-(1+|-cos26)-tgy1 (z)sin26-cos ДФz +
+|2 • sin2 29 + (1 -1 • cos29)2 • tg2 y1 (z)П°'5
. (11)
+(1 +1 • cos29)2 +12 • tg2 y1 (z) • sin2 29
Начальное значение фазы комплексного фазора рассеянной волны находится по формуле
9ph(z) = (z)-^(z) , (12)
где (z) = argÉyPH(z); ^2 (z) = argÉxPH (z) , и будет определяться следующим выражением: д • sin20 • tg y1 (z)
фpн(z) = arctg sin^z
cosДФz
- arctg
((1 + ц • cos29) + ц • sin29 • tg Y1 (z))
Г (1- ц • cos29 )tg Yl (z ) (ц • sin29 + (1 - ц • cos29 ) • tg Y1 (z ))
sinДФz
, (13)
Подставив значения уРН (г) и фРН (г), полученные по формулам (11) и (13), в формулы (5), (6) и используя полученные (после указанной подстановки) формулы (2), (3), получим выражения для амплитуды \Ёруу (сигнала, принятого с вертикальной поляризацией, и амплитуды Ёрхх (г )| сигнала,
принятого с горизонтальной поляризацией.
Осуществим перемножение в (2) и (3) и оп делим модули полученных произведений |ЕЁрХх (г)
реи
E
Рху (z)| , тогда получим
|£pxy (z) = [cos2 ар (z) • sin2 (Рн - Pp (z,Ph )) + +sin2ар(z)• cos2 (Рн -Pp(z,Ph))]0,5, (14)
\Épxx (z) =[cos2 аР (z)•cos2 (рН - pp (z, рН )) +
+sin2ap(z)• sin2 (Рн -Pp(z,Ph))]0,5. (15)
Анализ полученных результатов
Расчеты величины ЬБЯ(г) осуществлялись по формуле (1) с использованием значений амплитуд (14), (15).
На рис. 1 приведены результаты расчетов ЬБЯ(г) для различных значений степени поляризационной анизотропии д при интенсивности дождя Я = 12,5 мм/ч.
£ S "
-40
30 40 50 Длина трассы, z, км
£ S
-20
-40
10
20 б
30 40 50 Длина трассы, z, км
£ S
-40
30 40 50 Длина трассы, г, км
в
Рис. 1. Зависимость линейного деполяризационного отношения от длины трассы при интенсивности осадков 12,5 мм/ч, рН = 45° и степени поляризационной анизотропии д: а - 0,1; б - 0,2; в - 0,3
а
0
0
Расчеты приведены для следующих значений 9: 1 - 9 = О°; 2 - 9 = 11,25°; 3 - 9 = 22,5°; 4 - 9 = 33,75°; 5 - 9 = 45°; 6 - 9 = 56,25°; 7 - 9 = 67,5°; 8 - 9 = 78,75°.
Расчеты приведены для значения рн = 45°. С одной стороны, при этом значении рн происходят наиболее характерные изменения оценки ЬБЯ(г), с другой стороны, согласно техническому описанию ДМРЛ-С [6], излучение сигнала с линейной поляризацией с углом ориентации плоскости поляризации рн = 45° является одним из наиболее характерных режимов работы радиолокатора.
Расчеты проводились для линейно поляризованных излучаемых волн трехсантиметрового диапазона, для которых величины Я, Да, ДФ составляют следующие значения [4, 14]:
Да = О,О2 дБ/км; ДФ = 1 град/км при Я = 12,5 мм/ч; Да = О,1 дБ/км; ДФ = 4 град/км при Я = 5О мм/ч; Да = О,8 дБ/км; ДФ = 14 град/км при Я = 15О мм/ч.
При этом под дифференциальным фазовым сдвигом следует понимать фазовый сдвиг между ортогональными компонентами облучающей волны, возникающий при распространении радиолокационного сигнала в среде, заполненной гидрометеорами. наличие этого фазового сдвига установлено экспериментально в работах [14, 15].
Предел изменения значений степени анизотропии д при расчетах выбран с учетом их соответствия минимальным и максимальным размерам капель, дающим основной вклад в радиолокационную отражаемость и представляющими значимый интерес при решении задач метеорологической радиолокации [15, 16].
Как следует из расчетных зависимостей величины ЬБЯ(г), для случая интенсивности осадков 12,5 мм/ч характер изменения линейного деполяриза-ционного отношения подвержен изменению в зависимости от угла ориентации собственного поляризационного базиса анизотропного участка и степени поляризационной анизотропии этого участка.
на рис. 1, а величина ЬБЯ(г) принимает значение при длине трассы распространения г = О от -14 до -33 дБ в зависимости от угла ориентации собственного базиса 9. При длине трассы г = 5О км значение ЬБЯ(г) для различных 9 находится в диапазоне от -15 до -17 дБ. Подобное изменение расчетного линейного деполяризационного отношения также характерно для случаев, показанных на рис. 1, б и в. При этом величина ЬБЯ(г) принимает значение при г = О от -1О до -33 дБ для рис. 1, б и от -7 до -33 дБ для случая представленного на рис. 1, в, а при г = 5О км - от -12 до -14 и от -9 до -11 дБ соответственно.
Также, как видно из рис. 1, а, имеет место изменение характера ЬБЯ(г) с убывающего на возрастающий. наиболее выраженным случаем (см. рис. 1, а) является расчетная зависимость ЬБЯ(г) при 9 = 67,5° при длине трассы распространения от О до 1О км, значение которой изменяется от -33 (г = О км) до -38 дБ (г = 5 км).
Расчеты также показывают, что при интенсивности осадков Я = 12,5 мм/ч и значениях углов рн < 45° наблюдается изменение характера зависи-
мости ЬОЯ(г) по сравнению с зависимостями, представленными на рис. 1, а-в. При рн = О°, д = О,1 разброс значений ЬБЯ(г) составляет порядка 1О дБ (между значениями ЬБЯ(г) при 9 = О° и значениями ЬБЯ(г) при 9 = 33,75°) на всем протяжении длины трассы распространения (при г = 5О км), при этом изменения ЬБЯ(г) имеют возрастающий характер. Для значений углов рн > 45° характер изменения ЬБЯ(г) убывающий. В случае рн = 78,75°, д = О,1 имеет место разброс значений ЬБЯ(г) до 4 дБ (между значениями ЬБЯ(г) при 9 = 67,5° и значениями ЬБЯ(г) при 9 = 78,75°).
на рис. 2 представлены расчетные зависимости ЬБЯ(г) для случая интенсивности дождя 5О мм/ч и рн =45°. Значения 9 соответствуют значениям на рис. 1.
Изменение характера зависимости ЬВЯ(г) происходит на интервале г от О до 4 км, при этом имеет место перепад, составляющий порядка 6 дБ (см. рис. 2, а для случая 9 = 67,5°), тогда как в предыдущем случае интервал г соответствующий подобным характерным изменениям ЬБЯ(г), составляет от О до 1О км.
В отличие от предыдущего случая (Я = 12,5 мм/ч), при длине трассы г от 4О до 5О км появляется повторная смена характера изменения расчетной величины ЬВЯ(г): возрастающий характер меняется на убывающий с наличием характерных пиков (до 18 дБ для случая 9 = 78,75°, см. рис. 2, в) на графике ЬБЯ(г). С увеличением д эти пики становятся более выраженными.
на длине трассы распространения до 11 км весьма выражено влияние угла 9 на расчетную величину: разброс значений ЬБЯ(г) составляет порядка 2О дБ при изменении 9 от 67,5 до 78,75° (см. рис. 2, а) и от 45 до 78,75° (см. рис. 2, б и в).
Расчеты также показывают, что при интенсивности осадков Я = 5О мм/ч и значениях углов рн < 45° наблюдается уменьшение влияния значения 9 при всех значениях степени поляризационной анизотропии д. При д = О,3 и рн = О° разброс значений ЬВЯ(г) уменьшается до 1О дБ. Причем в области положительных значений ЬБЯ(г) указанный разброс также уменьшается до 1О дБ. Для случая, когда значения углов рн > 45°, также наблюдается существенное уменьшение влияния 9 на величину ЬБЯ(г) при всех д. При д = О,3 и рн =78,75° разброс значений ЬБЯ(г) уменьшается до 12 дБ, при этом ЬБЯ(г) находится в области отрицательных значений для любых значений 9.
Выраженный характер изменений имеют расчетные зависимости величины ЬБЯ(г) при интенсивности осадков 15О мм/ч и том же угле начальной ориентации измерительного базиса относительно базиса среды рн = 45° (рис. 3). Значения 9 соответствуют значениям на рис. 1.
Для всех расчетных зависимостей величины ЬБЯ(г), представленных на рис. 3, имеют место участки г со сменой характера изменений с возрастающего на убывающий. Появляющиеся при таком характере изменений перепады достигают величины порядка 12 дБ (см. рис. 3, в для 9 = 78,75° в диапазоне г от О до 8 км), при этом сама величина ЬБЯ(г) изменяется в пределах значений от -43 (см. рис. 3, б для 9 = 45°) до 4 дБ (см. рис. 3, б для всех значений 9).
20
0 £ 3 4 8 1 2
N ч
ч -20 " -у
-40 7
10
0 10 20 30 40 50
а Длина трассы, г, км
20
10 20 30 40 50 а Длина трассы, г, км
10
£
S к S
-20
-40
2
8 1 3
S'
10 20 30 40 50 Длина трассы, z, км
10
20 б
30 40 50 Длина трассы, z, км
10
£
-20
10
20
в
30 40 50 Длина трассы, z, км
0
10
20
Рис. 2. Зависимость линейного деполяризационного отношения от длины трассы при интенсивности осадков 50 мм/ч, рН = 45° и степени поляризационной анизотропии ц: а - 0,1; б - 0,2; в - 0,3
С увеличением степени поляризационной анизотропии влияние 9 на величину ЬБК^) увеличивается.
Расчеты также показывают, что при интенсивности осадков равной 150 мм/ч, в случае когда углы рН < 45°, зависимость ЬВК(г) будет аналогичной приведенной на рис. 3 а-в.
30 40 50 в Длина трассы, г, км
Рис. 3. Зависимость линейного деполяризационного отношения от длины трассы при интенсивности осадков 150 мм/ч, рН = 45° и степени поляризационной анизотропии ц: а - 0,1; б - 0,2; в - 0,3
Однако минимальные значения величины линейного деполяризационного отношения увеличиваются до значений порядка -20 дБ (рН = 0°, 9 = 0° и ц = 0,3). В случае значений углов рН > 45° для указанной интенсивности осадков также сохраняется характер зависимости ЬБК^), приведенной на рис. 3 а-в от -25 дБ при 9 = 67,5° до -5 дБ при 9 = 0 для случая Рн = 78,75°, ц = 0,1.
0
0
0
0
0
В указанных расчетных зависимостях угол 9 с физической точки зрения соответствует изменению ориентации поляризационного базиса анизотропного участка, возникающему вследствие изменения геометрической ориентации капель дождя под воздействием ветрового сдвига и (или) статического электричества [5, 15]. В приведенных расчетах указанное влияние анизотропного участка проявляется в использовании матрицы рассеяния S, содержащей тригонометрические функции угла 9, вследствие чего зависимости LDR(z) носят немонотонный характер.
Заключение
При дистанционном зондировании неоднородной вдоль трассы среды распространения, в отличие от однородной среды [8], появляется ряд характерных особенностей изменчивости линейного деполя-ризационного отношения, связанных как с воздействием дифференциальных факторов среды распространения, так и, в особенности, с изменением угла ориентации собственного поляризационного базиса анизотропного участка и изменением степени поляризационной анизотропии этого участка. Эти особенности проявляются в выраженном уменьшении величины линейного деполяризационного отношения на малых дальностях (до величин LDR(z) порядка 40 дБ) с последующим резким возрастанием и изменением знака величины LDR(z) при интенсивности осадков 50 и 150 мм/ч, причем в последнем случае изменение знака величины LDR(z) наблюдается при гораздо меньшей длине трассы распространения порядка 12 км.
Перечисленные особенности изменчивости линейного деполяризационного отношения позволяют обеспечить дистанционное получение информации для принятия решения о наличии в зондируемой области метеообразований
1) процессов, связанных с поперечными сдвигами ветра;
2) зон повышенного статического электричества (зон грозовой активности).
Приведенные в данной статье результаты расчетов справедливы, как было указано выше, для радиолокационных сигналов трехсантиметрового диапазона.
Поскольку указанные выше явления представляют собой опасные факторы для различных областей хозяйственной деятельности, например полетов воздушных судов гражданской авиации, результаты, полученные в данной работе, имеют практическую значимость.
Литература
1. Galletti M. Measurement and characterization of entropy and degree of polarization of weather radar targets / M. Galletti, D.H.O. Bebbington, M. Chandra, T. Boerner // IEEE T. Geosci. Remote. - 2008. - Vol. 46, No. 10. - P. 3196-3207.
2. Galletti M. Fully polarimetric analysis of weather radar signatures / M. Galletti, D.H.O. Bebbington, M. Chandra, T. Boerner // IEEE Radar Conf. - 2008. - P. 561-566.
3. Масалов Е.В. Точностные характеристики поляризационного радиолокатора при зондировании среды, заполненной гидрометеорами / Е.В. Масалов, С.В. Янов // Доклады ТУСУР. - 2011. - № 2(24), ч. 3. - С. 17-19.
4. Масалов Е.В. Трансформация линейно поляризованных электромагнитных волн в средах, содержащих гидрометеоры // Матер. междунар. конф. «Актуальные проблемы электронного приборостроения». - Новосибирск: НГТУ, 2010. - С. 77-79.
5. Масалов Е.В. Поляризационные измерения в задачах радиолокационной метеорологии / Е.В. Масалов, В.Н. Татаринов // Зарубежная радиоэлектроника. - 1987. -№ 4. - С. 44-52.
6. Вовшин Б.М. Теория и практика поляризационных измерений в метеорологической радиолокации [Электронный ресурс] / Б.М. Вовшин, И.С. Вылегжанин,
B.Ю. Жуков и др. // Матер. V Всерос. конф. «Вторые Всерос. Армандовские чтения» (Муром). - 2012. - С. 49-54. -Режим доступа: http://www.mivlgu.ru/conf/armand2012/pdf/ P5.pdf, свобод-ный (дата обращения: 25.05.2018).
7. Galletti M. Degree of polarization at simultaneous l transmit: Theory and applications for weather radar / M. Galletti, D.S. Zrnic, V.M. Melnikov, R.J. Doviak // IEEE T. Geosci. Remote. - 2012. - Vol. 9, No. 3. - P. 383-3871.
8. Масалов Е.В. Анализ влияния гидрометеоров на величины дифференциальной радиолокационной отражаемости и линейного деполяризационного отношения радиолокационных сигналов / Е.В. Масалов, Н.Н. Кривин, К.В. Кокоулин // Доклады ТУСУРа. - 2017. - № 2(20). -
C. 28-35.
9. Eigenvalue Signal Processing for Weather Radar Po-larimetry: Removing the Bias Induced by Antenna Coherent Cross-Channel Coupling / M. Galletti, D. S. Zrnic, F. Gekat, P. Goelz // IEEE T. Geosci. Remote. - 2014. - Vol. 52, No. 12. - P. 7695-7707.
10. Galletti M. Zenith/Nadir pointing mm-wave radars: linear or circular polarization? / M. Galletti, D. Huang, P. Kollias // IEEE T. Geosci. Remote. - 2014. - Vol. 52, No. 1. - P. 628-639.
11. Montopoli M. Meteorological radar systems / M. Montopoli, F.S. Marzano // Integrated ground-based observing systems. - Berlin: Springer, Heidelberg, 2010. - P. 33-57.
12. Татаринов В.Н. Введение в современную теорию поляризации радиолокационных сигналов. Поляризация плоских электромагнитных волн и её преобразования: учеб. пособие / В.Н. Татаринов, Л.П. Лигтхарт, С.В. Татаринов. - Томск: ТУСУР, 2012. - Т. 1. - 380 с.
13. Масалов Е. В. Представление матрицы рассеяния метеоцели / Е.В. Масалов, В.А. Потехин, В.Н. Татаринов // Изв. высш. учеб. завед. Физика. - Томск, 1983. - № 7. - С. 1-10.
14. Родимов А.П. Особенности использования поляризационных параметров ЭМВ в системах связи миллиметрового диапазона / А.П. Родимов, В.В. Поповский, В.И. Дмитриев // ЗРЭ. - 1980. - № 7. - C. 25-37.
15. Огути Т. Распространение и рассеяние электромагнитных волн в дожде и других гидрометеорах // ТИИЭР. - 1983. - Т. 71, № 9. - С. 6-65.
16. Экспериментальные исследования анизотропии рассеяния радиолокационных сигналов облаками и осадками / Н.Н. Бадулин, А.П. Бацула, Е.Б. Кульшенева и др. // Изв. АН СССР. Физика атмосферы и океана. - 1984. -Т. 20, № 6. - С. 505-510.
Масалов Евгений Викторович
Д-р техн. наук, профессор каф. конструирования и производства радиоаппаратуры (КИПР) Томского государственного университета систем управления и радиоэлектроники (ТУСУР) Ленина пр-т, д. 4О,г. Томск, Россия, 634О5О Тел.: +7-913-846-42-13 Эл. почта: e-v-masalov@yandex.ru
Кривин Николай Николаевич
Канд. техн. наук, доцент каф. КИПР ТУСУРа Ленина пр-т, д. 4О, г. Томск, Россия, 634О5О Тел.: +7-913-882-91-87 Эл. почта: freeman46@yandex.ru
Кокоулин Константин Владимирович
Студент каф. КИПР ТУСУРа
Ленина пр-т, д. 4О, г. Томск, Россия, 634О5О
Тел.: +7-913-845-15-О6
Эл. почта: jelikavar@gmail.com
Masalov E.V., Krivin N.N., Kokoulin K.V. Analysis of the linear depolarization ratio variability in the radiolocation sensing of inhomogeneous medium filled with hydrometeors
The theoretical results for estimating the linear depolarization ratio are obtained by using the Jones vector in eigenbasis of the homogeneous propagation medium and the complex phas-or of the radar signal scattered by the second medium region which follows directly after the first, homogeneous medium region. The second medium region is characterized by anisot-ropic polarization properties and orientation angle of its eigenbasis. Herewith changes in the polarization structure of the radar signal in the process of its propagation in this medium are analyzed. Such changes depend on initial orientation angle of the polarization basis of the radar relative to the eigenbasis of the homogeneous region, the degree of anisotropy of the second region and orientation angle of its eigenbasis. Linear depolarization ratio calculation results in the radar signal reception point for the different rain intensity are given. Obtained results are analyzed.
Keywords: polarization, ellipticity angle, polarization ellipse orientation angle, differential attenuation, differential phase shift, orientation angle of the anisotropic region's eigenbasis, polarization anisotropy degree, linear depolarization ratio, precipitation.
doi: 10.21293/1818-0442-2018-21-3-7-13
References
1. Galletti M., Bebbington D.H.O., Chandra M., Boerner T. Measurement and characterization of entropy and degree of polarization of weather radar targets. IEEE T. Geosci. Remote, 2008, vol. 46, no. 10, pp. 3196-207.
2. Galletti M., Bebbington D.H.O., Chandra M., Boerner T. Fully polarimetric analysis of weather radar signatures. IEEE Radar Conf., 2008, pp. 561-566.
3. Masalov E.V., Yanov S.V. Punctuality parameters of polarization radar at the hydrometeors remote sensing. Proceedings of TUSUR University, 2011, vol. 24, no. 2, pp. 17-19.
4. Masalov E.V. Transformation of linearly polarized electromagnetic waves in media containing hydrometeors. Proceedings of international conference «Actual problems of electronic instrument making», Novosibirsk, NGTU, 2010, pp. 77-79.
5. Masalov E.V., Tatarinov V.N. Polarization measurements in radar meteorology. Russian physics journal, 1987, no. 4, pp. 44-52.
6. Vovshin B.M., Vyligjanin I.S., Zhukov V.U., Pushkov A.A., Schukin G.G. Theory and practice of polarization measurements in meteorological radiolocation. Proceedings of 5th Russian conference «Second Russian Armand Readings» (Murom), 2012, pp. 49-54 (In Russ.). Available at: http://www.mivlgu.ru/conf/armand2012/pdf/P5.pdf (accessed: 25 May 2018).
7. Galletti M., Zrnic D.S., Melnikov V.M., Doviak R.J. Degree of polarization at simultaneous l transmit: Theory and applications for weather radar. IEEE T. Geosci. Remote, 2012, vol. 9, no. 3, pp. 383-3871.
8. Masalov E.V., Krivin N.N., Kokoulin K.V. An analysis of the influence of hydrometeors on the linear depolarization ratio and differential reflectivity of the radar signals. Proceedings of TUSUR University, 2017, vol. 20, no. 2, pp. 28-35.
9. Galletti M., Zrnic D. S., Gekat F., Goelz P. Eigenvalue Signal Processing for Weather Radar Polarimetry: Removing the Bias Induced by Antenna Coherent Cross-Channel Coupling. IEEE T. Geosci. Remote, 2014, vol. 52, no. 12, pp. 7695-7707.
10. Galletti M., Huang D., Kollias P. Zenith/Nadir pointing mm-wave radars: linear or circular polarization? IEEE T. Geosci. Remote, 2014, vol. 52, no. 1, pp. 628-639.
11. Montopoli M., Marzano F.S. Meteorological radar systems. Integrated ground-based observing systems. Springer, Berlin, Heidelberg, 2010, pp. 33-57.
12. Tatarinov V.N., Ligthart L.P., Tatarinov S.V. Vvedenie v sovremennuju teoriju polarizacii radiolocacionnyh signalov. Polarizacija ploskih electromagnitnih voln i eje preobra-zovanija. [An introduction to the modern polarization theory of radar signals. The polarization of a plane electromagnetic waves and its transformations]. Tomsk, TUSUR, 2012, vol. 1, 380 p. (In Russ.).
13. Masalov E.V., Potehin V.A., Tatarinov V.N. Representation of the meteorological objects scattering matrix. Russian Physics Journal, Tomsk, 1983, no. 7, pp. 1-10. (In Russ.)
14. Rodimov A.P., Popovsky V.V., Dmitriev V.I. Features of using the EMW polarization parameters in millimeter-wave communication systems, Russian physics journal, 1980, no. 7, pp. 25-37.
15. Oguti T. Rasprostranenie i rassejanie electromagnit-nih voln v dojde I drugih gidrometeorah.[Electromagnetic waves propagation and scattering in rain and other hydrome-teors], IEEE, 1983, vol. 71, no. 9, pp. 6-65 (In Russ.)
16. Badulin N.N., Bacula A.P., Kulsheneva E.B. Experi-mentalnie issledovanija anizotropii rassejanija radiolokacion-nyh signalov oblakami i osadkami [Experimental investigation of anisotropy of radar signal scattering by clouds and precipitation]. Izvestiya. Atmospheric and Oceanic Physics, 1984, vol. 20, no. 6, pp. 505-510.
Evgenii V. Masalov
Doctor of Technical Sciences, Professor,
Department of Radio Equipment Design and Production
Tomsk State University of Control Systems
and Radioelectronics (TUSUR)
40, Lenina pr., Tomsk, Russia, 634050
Phone: +7-913-846-42-13
Email: e-v-masalov@yandex.ru
Nikolai N. Krivin
PhD of Technical Sciences, Assistant Professor,
Department of Radio Equipment Design
and Production TUSUR
40, Lenina pr., Tomsk, Russia, 634050
Phone: +7-913-882-91-87
Email: freeman46@yandex.ru
Konstantin V. Kokoulin
Student, Department of Radio Equipment Design
and Production TUSUR
40, Lenina pr., Tomsk, Russia, 634050
Phone: +7-913-845-15-06
Email: jelikavar@gmail.com