УДК 621.396.96
Е.В. Масалов, Н.Н. Кривин, К.В. Кокоулин
Анализ влияния гидрометеоров на величины дифференциальной радиолокационной отражаемости и линейного деполяризационного отношения радиолокационных сигналов
На основе анализа компонент вектора Джонса в собственном базисе среды распространения получены расчетные соотношения для оценки величин дифференциальной радиолокационной отражаемости и линейного деполяризационного отношения радиолокационных сигналов. При этом анализируются изменения радиолокационного сигнала, возникающие в процессе его распространения в среде, заполненной гидрометеорами. Под средой распространения понимается воздушное пространство, заполненное частицами осадков: каплями дождя, частицами града, снежинками. Приводятся результаты расчетов этих величин в точке приема радиолокационного сигнала для осадков в виде дождя различной интенсивности. Полученные результаты анализируются. Ключевые слова: поляризация, угол эллиптичности, угол ориентации эллипса поляризации, дифференциальное ослабление, дифференциальный фазовый сдвиг, дифференциальная радиолокационная отражаемость, линейное деполяризационное отношение. doi: 10.21293/1818-0442-2017-20-2-28-32
Дистанционное зондирование метеообразований с использованием поляризационной структуры сигналов, с одной стороны, позволяет значительно увеличить количество информации о метеообразованиях, с другой стороны, сама поляризационная структура радиолокационных сигналов оказывается в значительной степени подвержена искажающему воздействию факторов среды [1, 2]: дифференциальному ослаблению Да [дБ/км] и дифференциальному фазовому сдвигу ДФ [град/км].
Кроме того, ещё одним негативным фактором, влияющим на изменчивость дифференциальной радиолокационной отражаемости и линейного деполяризационного отношения (ЬБК), является рассогласование углов ориентации собственного поляризационного базиса среды распространения и измерительного базиса, связанного с РЛС.
Постановка задачи
Указанные выше факторы являются причиной неадекватности результатов оценки интенсивности осадков и алгоритмов распознавания опасных зон метеообразований на основе поляризационных измерений с применением двух ортогональных линейных поляризаций [3]. Перечисленные обстоятельства обусловливают необходимость поиска путей, позволяющих оценить значимость влияния указанных факторов на получаемые результаты измерения и определить характер изменчивости измеряемых параметров и ЬБЯ^), где г - длина трассы распространения радиолокационного сигнала в среде, заполненной гидрометеорами.
Методика решения
В данной работе величина ЬБК^) определяется соотношением
1БЯ( г) = 20 • ЦЁху (г )| / \Ёхх (г)|), (1)
где |Ё^ (г)| - амплитуда сигнала, излученного с горизонтальной поляризацией и принятого с вертикальной; \ЁХХ (г)| - амплитуда сигнала, излученного
и принятого с горизонтальной поляризацией.
Величина дифференциальной радиолокационной отражаемости г) на выходе приёмника с логарифмической характеристикой может быть определена в виде [1]
(г) = 20 • ЩЁхх (41\Ётт (г)), (2)
где |Ёхх (г) - амплитуда сигнала, излученного с горизонтальной поляризацией и принятого той же антенной; |Ёуу (г) - амплитуда сигнала, излученного с вертикальной поляризацией и принятого той же антенной.
Для нахождения указанных амплитуд использована методика определения компонент вектора Джонса в собственном базисе среды распространения [4] с последующим переходом в собственный базис РЛС и определением сигналов на выходе прямоугольных волноводов, принимающих сигналы горизонтальной и вертикальной поляризации соответственно.
Тогда комплексные амплитуды, необходимые для определения ЬБЯ^), могут быть определены из приведенных ниже выражений:
exx (z) =
Exy (z) =
1 0 0 0
[я_1фн)][ф( z,ßH))]
[^_1(ßH)][^(ß( z,ßH))]
cosa( z) j sin a(z)
cos a( z) j sin a( z)
(3)
. (4)
Для оценки параметра помимо указанной
выше комплексной амплитуды Ёхх (г), необходимо
определить величину Ёуу (г), для чего достаточно
учесть, что 7-компонента излучённого сигнала повёрнута относительно Х-компоненты на 90°, поэтому к углу рН следует добавить 90°. Также при расчетах коэффициента эллиптичности а(.г) и в угле ориентации в(2) для этой компоненты необходимо будет использовать начальный угол наклона плоскости поляризации, равный рН + 90°. Тогда
éyy (z) =
0 0 0 1
cos^ sin(
(Ри + 90o) -sin ( + 90o) in(pH + 90o) cos(pH + 90o)
p(z,pH + 90o) sin p(z,pH + 90o) inp(z ,pH + 90o) cosp(z,pH + 90o)
cos a (z, pH + 90°) j sin a(z, pH + 90o) В формулах (3), (4) и (5) использованы сле
cosp -sinp
(5)
дующие обозначения:
cos a( z) j sina(z)
- вектор Джонса
эллиптически поляризованной волны в базисе сре-
ды; [[z,ph ))] =
cosP(z,pH) sin Р( z ,ph ) -sinP(z,pH) cosP(z,pH)
- опе-
ратор перехода в базис среды из базиса локатора;
cospH -sinpH sinpH cospH
R_1(PH)
- оператор перехода в
базис локатора из базиса среды;
1 0 0 0
матрица
Джонса прямоугольного волновода, принимающего
"0 0"
сигнал с горизонтальной поляризацией;
0 1
матрица Джонса прямоугольного волновода, принимающего сигнал с вертикальной поляризацией; а и р^) - характеристики эллипса поляризации, представляющие собой его угол эллиптичности и его угол наклона соответственно; рн - угол рассогласования поляризационных базисов РЛС и среды распространения зондирующего сигнала.
Осуществим перемножение в (2) и (3) и определим модули полученных произведений Ехх (2)| и
\Ёху (2)|, тогда получим
(6)
(7)
\ÉXY (z)| = [cos2 a(z)• sin2 ( - P(z,Рн)) +
+sin2a(z)cos2(Ph -P(z,Ph))]0,5, \ÉXX (z)| = [cos2 a(z) • cos2 ( - P(z, Ph )) + +sin2 a(z) sin2 (Ph-P(z,Ph ))]0,5.
Для амплитуды |Éyy (z)|, исходя из формулы (5), получим
|ÉYY (z) = ^cos2 a(z,pH + 90o) • sin2 ((pH + 90o) --p(z ,pH + 90o)) + sin2 a(z ,pH + 90o )x
<cos'
((ph + 90o)-p((H + 90o))
0,5
(8)
Для расчета зависимостей а^) и р^) в формулах (3)-(5) воспользуемся известными по [1, 2] соотношениями:
a( z) =—arcsin 2
(2 • tgPH J00,05Aaz • sin(AOz)Л
1 + i00,1Aaz • tg2 PH
(9)
P( z) = 2 • arctg
2 40
0,05 Aaz
tgPH ^cos(A®^ z)
Л
1-1001Aaz • tg2 PH
+™Л10)
где рН - угол ориентации измерительного поляризационного базиса (связанного с РЛС) относительно собственного базиса среды.
На рис. 1 и 2 приведены результаты расчетов а^) и р^) соответственно, полученные для интенсивности осадков Я = 12,5 мм/ч, при которой величины Да и ДФ составляют 0,02 дБ/км и 1 град/км соответственно [2].
30 40 50
Длина трассы, z, км Рис. 1. Зависимость изменения угла эллиптичности от длины трассы
3 loo
Длина трассы, г. км Рис. 2. Зависимость изменения угла ориентации от длины трассы
Расчеты приведены для значения рН: 1 - рН = 0°; 2 - рН = 11,25°; 3 - рН = 22,5°; 4 - рН = 33,75°; 5 - Рн = 45°; 6 - РН = 56,25°; 7 - РН = 67,5°; 8 - рН = 78,75°; 9 - рН = 90°.
При интенсивности Я = 12,5 мм/ч характер изменения характеристик а^) и в(г) носит в целом близкий к линейному характер изменения на трассе протяженностью 50 км, величина а^) не превышает значения 23 град; Р^) изменяется от 0 до 90 град. Угол эллиптичности а^) в случае начального угла
ориентации эллипса поляризации рН = 45° по мере увеличения расстояния от точки излучения линейно возрастает и достигает максимального значения, равного 23 град. Угол ориентации р^) в зависимости от каждого начального угла рН ориентации эллипса поляризации изменяется слабо.
Определяя величину дифференциальной радиолокационной отражаемости г) на выходе приёмника с логарифмической характеристикой в виде (2), после подстановки формул (7) и (8) получим
(г) = 10• Ц|Ёху (г)2/\Ёгт (г)2) . (11)
Анализ полученных результатов
На рис. рисунке 3 приведены результаты расчетов ЬБЯ^) для различных значений интенсивности дождя Я. Значения рН соответствуют приведенным выше для рис. 1 и 2.
а
Длина трассы, г, км
б
Длина трассы, 2. км
в
Рис. 3. Зависимость линейного деполяризационного отношения от длины трассы при интенсивности осадков: а - 12,5 мм/ч; б - 50 мм/ч; в - 150 мм/ч
Результаты расчетов по формуле (11) приведены на рис. 4-9.
Расчеты проводились для однородной среды и линейно поляризованных излучаемых волн трехсантиметрового диапазона, для которых величины Я, Да, ДФ составят следующие значения [2]: Да = 0,02 дБ/км; ДФ = 1 град/км - при Я = 12,5 мм/ч; Да = 0,1 дБ/км; ДФ = 4 град/км - при Я = 50 мм/ч; Да = 0,8 дБ/км; ДФ = 14 град/км - при Я = 150 мм/ч.
Из рис. 3 видно, что для случая Я = 12,5 мм/ч величина ЬБЯ^) линейно возрастает и достигает значения порядка -7 дБ при рН = 45° и г = 50 км.
Как следует из расчетных зависимостей величины ЬБЯ^), для интенсивности дождя Я = 50 мм/ч (см. рис. 3, б) также имеет место монотонность зависимостей ЬБЯ^) от длины трассы г. Однако в этом случае появляются участки г, на которых возрастающий характер изменяется на убывающий. В данном случае это проявляется при начальных углах ориентации рН = 22,5; 33,75; 45 и 56,25°. Появляющийся при таком характере изменений перепад составляет величину порядка 15 дБ.
Еще более выраженной является изменчивость характера монотонной зависимости для случая интенсивности дождя Я = 150 мм/ч, как видно из результатов расчета ЬБЯ^), приведенных на рис. 3, в. В этом случае на трассе длиной г = 50 км имеют место три участка с таким же, как и в предыдущем случае, характером изменения направленности. При этом в случае рН = 22,5; 33,75 и 45° указанные особенности проявляются уже на интервале г от 0 до 15 км. Максимальный перепад в оценке ЬБЯ^) достигает порядка 20 дБ.
Далее рассмотрим результаты расчетов дифференциальной радиолокационной отражаемости 1ОЯ (г), приведенные на рис. 4-9.
Дальность, г, км
Рис. 4. Зависимости дифференциальной радиолокационной отражаемости (г) при различных значениях рН для интенсивности дождя Я = 12,5 мм/ч
Угол рН изменялся от 0 до 90° с дискретом 11,25°. На рис. 4, 6, 8 начальный угол наклона плоскости поляризации излучаемых колебаний рН принимает значения: 1 - 0°; 2 - 11,25°; 3 - 22,5°; 4 -33,75°; 5 - 56,25°; 6 - 67,5°; 7 - 78,75°.
10 20 30 40 50 Дальность, 2, км
а
б Дальность, 2, км Рис. 5. Зависимость дифференциальной радиолокационной отражаемости ZDR (2) для интенсивности дождя Я = 12,5 мм/ч при рН = 45° (а) и 90° (б)
0,1
-0,41
-0,92
-1,43
30 35 40 45
Рис. 6. Зависимости дифференциальной радиолокационной отражаемости ZDR (2) при различных значениях рН для интенсивности дождя Я = 50 мм/ч
20
из
ч 0
£
-20
-40
-60
-6
-6,5
из
ч
-7
7,5
I 10 20 30 40 50
а Дальность, 2, км
б
40 50
Дальность, 2, км
Рис. 7. Зависимость дифференциальной радиолокационной отражаемости ZDR (2) для интенсивности дождя Я = 50 мм/ч при рН = 45° (а) и 90° (б)
0,1 -0,1
-0,9
-1,9
45 50
0 5 10 15 20 25
Дальность, 2, км Рис. 8. Зависимости дифференциальной радиолокационной отражаемости ZDR (2) при различных значениях рН для интенсивности дождя Я = 150 мм/ч
20
0
1-3
ч -20
к
-40
-60
0
10 20
30 40 50 " Дальность, 2, км
-6,5
из
ч
-7,5
10
20
б
30 40 50 Дальность, 2, км
Рис. 9. Зависимость дифференциальной радиолокационной отражаемости ZDR (2) для интенсивности дождя Я = 150 мм/ч при рН = 45° (а) и 90° (б)
Как видно из рис. 4, при интенсивности осадков 12,5 мм/ч изменения дифференциальной радиолокационной отражаемости имеют убывающий характер в зависимости от начального угла наклона рН плоскости поляризации излучаемых колебаний и от длины трассы. Оценка ZDR(г) изменяется от 0 до -0,959 дБ. Однако при рН = 45 и 90° зависимость ZDR (2) резко меняет характер и пределы изменения. При рН = 45° дифференциальная радиолокационная отражаемость увеличивается с увеличением длины трассы 2 в пределах от -60 до -5 дБ. При рН = 90° дифференциальная радиолокационная отражаемость уже уменьшается в пределах от -6 до -6,75 дБ. Последнее иллюстрируют рис. 5, а, б.
По мере увеличения интенсивности осадков до 50 мм/ч (см. рис. 6) дифференциальная радиолокационная отражаемость Z2) имеет также убывающий характер изменения в пределах от 0 до
0
а
0
-1,452 дБ. Так же, как и при интенсивности осадков 12,5 мм/ч, при рН = 45 и рН = 90° зависимость
(г) резко меняет характер и предел изменения.
При рН = 45° (см. рис. 7, а) для указанной выше интенсивности осадков значение дифференциальной радиолокационной отражаемости растёт с увеличением длины трассы г в пределах от -40 до 10 дБ. При рН = 90° (см. рис. 7, б) дифференциальная радиолокационная отражаемость уже уменьшается в пределах от -6 до -7 дБ. Из рис. 6 следует, что величина (г) зависит, помимо длины трассы г, как от ориентации собственного базиса среды распространения рН, так и от дифференциальных характеристик среды Да и ДФ.
Из рис. 8 для максимальной интенсивности осадков 150 мм/ч следует, что оценка г), начиная с дальности г = 25 км, приобретает ярко выраженную зависимость от угла ориентации собственного базиса среды распространения относительно измерительного базиса. Также можно показать, что при рН = 56,25° оценка г) на 26 км достигает минимального значения -0,8 дБ. В случае рН = 33,75° и рН = 56,25° имеют место участки с изменением характера зависимости г) с убывающего на возрастающий и наоборот. Максимальный размах изменчивости достигает 1,5 дБ.
При рН = 45° (см. рис. 9, а) дифференциальная радиолокационная отражаемость также увеличивается с увеличением длины трассы. Нижний и верхний пределы изменения в этом случае составляют -30 и 0 дБ соответственно. При рН = 90° (см. рис. 9, б) дифференциальная радиолокационная отражаемость уменьшается в пределах от -6 до -7,405 дБ на участке трассы длиной от 0 до 39 км, далее происходит монотонное возрастание до -6,8 дБ.
При увеличении интенсивности осадков с 12,5 до 50 мм/ч оценка (г) уменьшается до -1,452 дБ при длине трассы 50 км, а при интенсивности 150 мм/ч уже достигает -2,3 дБ. Это говорит о том, что чем больше степень влияния осадков на оценку (г), тем сильнее она подвержена резким изменениям и большим затуханиям в конце трассы.
Заключение
С увеличением интенсивности осадков влияние дифференциального фазового сдвига и дифференциального ослабления увеличивается, существенно увеличивая изменчивость ЬБЯ^). Для более точной интерпретации полученных данных учёт этого влияния представляет собой значимый фактор и требует дальнейших исследований, поскольку такой характер изменчивости ЬБЯ^) может быть использован в качестве дополнительного информационного признака.
В случае использования линейных ортогональных поляризаций имеет место выраженный характер рассогласования углов ориентации собственного поляризационного базиса среды распространения и измерительного базиса, связанного с РЛС.
Перечисленные выше обстоятельства позволяют утверждать, что в случае несовпадения базиса
среды и измерительного базиса при реализации метода, основанного на использовании линейной поляризации, принимаемый сигнал обладает поляризационными характеристиками, зависящими в большей степени от характеристик процесса распространения, а не от характеристик рассеивающего метеообъёма.
С увеличением интенсивности осадков изменчивость дифференциальной радиолокационной отражаемости носит всё более выраженный характер.
Литература
1. Масалов Е.В. Точностные характеристики поляризационного радиолокатора при зондировании среды, заполненной гидрометеорами / Е.В. Масалов, С.В. Янов // Доклады Том. гос. ун-та систем управления и радиоэлектроники. - 2011. - № 2(24), ч. 3. - С. 17-19.
2. Масалов Е.В. Трансформация линейно поляризованных электромагнитных волн в средах, содержащих гидрометеоры // Матер. междунар. конф. «Актуальные проблемы электронного приборостроения». - Новосибирск: НГТУ, 2010. - С. 77-79.
3. Теория и практика поляризационных измерений в метеорологической радиолокации [Электронный ресурс] / Б.М. Вовшин, И.С. Вылегжанин, В.Ю. Жуков и др. // Матер. V Всерос. конф. «Вторые Всероссийские Армандов-ские чтения» (Муром). - 2012. - С. 49-54. - Режим доступа: http://www.mivlgu.ru/conf/armand2012/pdf/P5.pdf, свободный (дата обращения: 11.12.2016).
4. Татаринов В.Н. Введение в современную теорию поляризации радиолокационных сигналов. Поляризация плоских электромагнитных волн и её преобразования: учеб. пособие / В.Н. Татаринов, Л.П. Лигтхарт, С.В. Татаринов. - Томск: ТУСУР, 2012. - Т. 1. - 380 с.
Масалов Евгений Викторович
Д-р техн. наук, профессор каф. конструирования и производства радиоаппаратуры (КИПР) ТУСУРа Тел.: +7 (913) 846-42-13 Эл. почта: [email protected]
Кривин Николай Николаевич
Канд. техн. наук, ст. преп. каф. КИПР ТУСУРа
Тел.: +7 (913) 882-91-87
Эл. почта: [email protected]
Кокоулин Константин Владимирович
Студент каф. КИПР ТУСУРа Тел.: +7 (999) 495-04-94 Эл. почта: [email protected]
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
Based on the analysis of the Jones vector components in their own basis of the propagation medium, the calculated ratios for estimating the value of the linear depolarization ratio and the differential reflectivity are obtained. The results of calculations of this value at the point of reception of the radar signal are given. The results are analyzed.
Keywords: polarization, ellipticity angle, the angle of orientation of the polarization ellipse, the differential attenuation, differential phase shift, linear depolarization ratio, differential reflectivity.