CC BY
53
УДК 621.396
Н. В. Кхыонг
Московский физико-технический институт (национальный исследовательский университет)
Оценка влияния метеобразования на распространение радиоволн в Х-диапазоне
В статье проанализировано ослабление электромагнитных волн (ЭМВ) сантиметрового диапазона в облаках, осадках и других метеорологических явлениях. Приведены аналитические выражения для расчета коэффициентов ослабления, а также примеры численного расчета коэффициентов для данных, полученных на макете посадочного радиолокатора (ПРЛ) Х-дпапазона.
Ключевые слова: ослабление, поглощение, рассеяние, посадочный радиолокатор, Х-днапазон, снтиметровый диапазон, водность, плотность, облако, осадка, газ, коэффициент ослабления, тропосфера, атмосфера.
N. V. Khuong
Moscow Institute of Physics and Technology
Evaluation of the influence of meteorology on the propagation of radio waves in X-bands
The papers analyzes the attenuation of electromagnetic waves (EMW) of the centimeter range in clouds, precipitation and other meteorological phenomena. The analytical expressions for calculating the attenuation coefficients are given, as well as examples of the numerical calculation of coefficients for the data obtained by the model of the X-band landing radar (PRL).
Key words: attenuation, absorption, scattering, landing radar, X-band, meter range, water content, density, cloud, precipitation, gas, attenuation coefficient, troposphere atmosphere.
1. Введение
Одной из причин, влияющей на выбор рабочего диапазона частот в радиолокации, является наличие ослабления на трассе распространения радиосигналов из-за осадков, облаков и других метеорологических явлений. Величина потерь различна для разных длин волн, причем на высоких частотах потери будут значительно больше.
Свободное пространство принято рассматривать как однородную, безграничную и непо-глощающую среду, у которой относительная диэлектрическая проницаемость равна единице. Таким образом, не принимается во внимание множество факторов, влияющих на распространение радиоволн. В атмосфере Земли всегда присутствуют атмосферные газы, пары воды, гидрометеоры (облако, дождь, снег, град, туман), а также частицы пыли, поднятые с поверхности Земли, которые вызывают ослабление мощности принимаемого сигнала. Принято рассматривать два вида ослабления - поглощение и рассеяние [6]. При поглощении часть энергии радиоволн превращается в тепловую энергию, а при рассеянии - перераспределяется в пространстве в направлениях, отличных от заданного. Ослабление в тропосфере сигналов сантиметрового диапазона имеет заметную величину.
© Н. В. Кхыонг, 2020
(с) Федеральное государственное автономное образовательное учреждение высшего образования
«Московский физико-технический институт (национальный исследовательский университет)», 2020
д
_ _ -2 [ а(К)<т
Р = Рс10 0 , (1)
где К - расстояние от радиолокатора до метеорологической цели, Р - средняя мощность на входе приемника без учета ослабления, Ро - средняя мощность на входе приемника с учетом ослабления,
а(К) - характеризует закон изменения ослабляющих свойств среды в направлении радиолуча.
В случае однородной среды распространения, когда коэффициент ослабления не изменяется вдоль линии распространения радиоволн, мощность принятого сигнала выражается как
Р = Р~010-2аК, (2)
где а - коэффициент ослабления (дБ/км).
При использовании логарифмического масштаба уравнение (2) будет иметь вид
101д= = -2аК. (3)
Для того чтобы охарактеризовать полные потери на трассе распространения, вводят коэффициент Г. Величина а измеряется в децибелах (дБ) и представляет собой сумму коэффициентов ослабления аг в газе и водяном паре, в облаках а0бл и осадках аос:
Г = 2(аг Дгг + аоблД^обл + «ос Д^ос )• (4)
Рассмотрим подробнее названные выше коэффициенты затухания каждой среды распространения ЭМВ в атмосфере.
Ослабление сантиметровых радиоволн в газах
Ослабление обусловлено рассеянием и поглощением энергии сигнала за счет преобразования части энергии в тепловую частицами газа (кислорода и азота) и водяным паром. При распространении радиоволн с длиной волны меньше Л = 3, 2 см ослабление в тропосфере происходит в основном за счет поглощения, эффект рассеяния незначительный, и им можно пренебречь. Тогда общий коэффициент ослабления аТ представляет собой сумму коэффициентов ослабления газа и водяного пара:
аг = а0 + ,
(5)
где а0 и ат - коэффициенты ослабления (дБ/км) газа и паров воды соответственно.
При частоте / меньше 54 ГГц, соответствующей длине волны больше 0,5 см коэффициент ослабления а0(дБ/км) в газе определяется выражением:
а„ =
0, 72г,
2,8
+
0,062{з
100 + 0, 34г2г1'6 (54 - 10)1,16?1 + 0, 836
6 = г\
0,0717 -1,8132 [0,0156(1-п)-1,6515(1-п)]
е
Р
Л2'
(6) (7)
6 = <5146г,
0,5146^-4,6368е[-0,1921(1-г4)-5,7416(1-г4)]
£з = г0,3414г-6,5851е[0,213(1-г4)-8,5854(1-г4)],
где / - часто та (ГГц),
I _ температура воздуха (°С), р - давление сухого воздуха (кПа),
(8) (9)
2
— р
1,013
И П
288 273+4 •
Для паров воды коэффициент ослабления ат (дБ/км) определяется как
= Ь 10 000^2,0^2^ 22) + л0 ^2 752) +
+
(10 - 22, 235)2 + 9, 42^2'
11, 96^1е0'7(1-Г4) (10 - 183, 31)2 + 11,14^2
(10 - 752)2 17,4^1е1'45(1-Г4) 25, 37^1е1'09(1-Г4) + (- 448)2 + (- 380)2 +
+-,ю 9( А, 557) + пп„ д(Х, 1780) +
+
(10 - 557)2 0, 081ще0,1(1-Г()
+
(I0 - 1780)2 3,66^1е1'6(1-п)
(Ю - 321, 226)2 + 6, 29^2 (10 - 325,153)2 + 9, 22^2
т)'
2,5
Ъ Р-
Значения ^1,^2 и функция д(А, /¿) равны:
(10)
= 0, 955грг0'68 + 0, 006р
^2 = 0, 735грг0'5 + 0, 0353г4р
5(А,/г) = 1 +
( 10 - А/г \
^10 + А/^
(11) (12)
(13)
где р - плотность паров воды [г/м3], /г - частота линии.
На рис. 1 приведен характер изменения коэффициента ослабления для водяного пара и газа при различных значениях длины волны Л и стандартных метеорологических условиях.
Длина волны А, см
Рис. 1. Зависимость коэффициентов ослабления сигнала от длины волны
Синяя кривая изображает зависимость коэффициента ослабления для водяного пара 7, 5 / 3 70
ослабления имеет максимум, равный 0, 34 дБ/км, при длине волны Л = 0, 45 см. Далее он значительно уменьшается при увеличении длины волны. Красная кривая соответствует газу при общем давлении 1001 гПа. Коэффициент ослабления для нее также убывает с
0, 1 10
0, 0055
г
р
0,13 [дБ/км] и имеет максимум при длине волны Л = 0, 45 см. Для длины волны рассматриваемого ПРЛ Л = 3, 2 см общий коэффициент ослабления сигнала в газах равен 0,011 дБ/км. На максимальной дальности ПРЛ (ДМаК = 40 км) потери энергии составляют 0,88 дБ.
Ослабление сантиметровых радиоволн в облаках
Жидкие и кристаллические частицы облаков без осадков имеют малый размер по сравнению с длиной волны Л (их диаметр не превышает 200-400 мкм), как правило, нас ПРЛ их не обнаруживает. Но они могут вызывать существенное ослабление радиолокационных сигналов по эффектам: поглощение и рассеяние. Данные эффекты не трудно описать с помощью формул:
для рассеяния:
для поглощения:
2^3 _ ( m2 - 1\
** = -*ГIm\-• (15)
Из сравнения формул (14) и (15) видно, что при малых размерах частиц эффективная площадь поглощения (ЭПП) гораздо больше эффективной площади рассеяния (ЭПР). Следовательно, общая эффективная площадь ослабления практически равна ЭПП жидкой частицы, так как ЭПР можно пренебречь.
Для расчета коэффициента ослабления воспользуемся понятием водность ш [г/м3]:
ш = 1 £ d3> (16)
i
где ро - плотность воды или льда.
Тогда общий коэффициент ослабления сигнала в жидких и кристаллических облаках выражается следующим образом:
-■к / m2 — 1 \
аобл = 0,4343—/т--2— (17)
A \ m2 + 2 J
Для приближенных расчетов при температуре 18 °С коэффициент ослабления в таких облаках может быть вычислен по формуле: для водного облака:
aw = 0,4343-^2 • (18)
А2
Коэффициент ослабления в ледяных или кристаллических облаках:
aic = 0, 024[дБ/км], (19)
А
где q - удельная влажность облака, равна массе ледяных частиц в граммах на 1 кг влажного воздуха.
При сравнении формул (18) и (19) показано, что коэффициент ослабления в кристаллических облаках примерно на два порядка меньше, чем в водяных, при одной влажности. При наблюдении метеорологических целей ослаблением сигнала в таких облаках можно пренебречь.
Таким образом, общий коэффициент ослабления в облаках равен коэффициенту ослабления в водяных облаках:
аобл = 0, 4343^, [дБ/км]. (20)
На рис. 2 представлены графики зависимости коэффициента ослабления от длины волны для облаков с различной средней водностью:
^pi =
Ослабления в облаке
10"' 10° 101
Длина волны А, см
Рис. 2. Зависимость коэффициента ослабления сигнала от длины волны для различных видов облаков
Из рис. 2 видно, что зависимости имеют вид параллельных линий, и коэффициент ослабления обратно пропорционален длине волны. При увеличении длины волны в десять раз от 1 до 10 см ослабление сигнала уменьшается в 100 раз. При длине волны Л = 3, 2 см и среднем значении водности ш = 1+2,5 г/м3 величина коэффициентов ослабления составляет более 0,042 дБ/км.
Таким образом, общий коэффициент ослабления для Х-диапазона в случае облаков без осадков составляет а0бл = 0, 005+0,14 дБ/км.
Ослабление сантиметровых радиоволн в виде дождя
Дождь состоит из частиц воды с различными размерами, формами и плотностями. В связи с тем, что для различных размеров частиц закон распределения но размерам отличается, а расчет параметров плотности распределения представляет собой довольно трудную задачу, наиболее подходящим методом расчета коэффициента ослабления является эмпирический метод, при котором считается, что зависимость коэффициента ослабления от интенсивности осадков выражена экспоненциальным законом [5|:
«обл = №/км], (21)
где параметры а и Ь эмпирические величины, которые зависят от длины волны и места расположения радиолокатора; I- интенсивность осадков [мм/ч].
Значения параметров а и Ь при температуре 18 °С приведены в табл. 1.
Результаты расчетов коэффициентов ослабления радиолокационного сигнала в дожде, в соответствии с формулой (21) при значениях длины волны 3,2 см, 5,6 см и 10 см представлены на рис. 3.
Графики, приведенные на рис. 3, демонстрируют, что коэффициент ослабления в дожде линейно возрастает с увеличением интенсивности. Коэффициент ослабления имеет значительную величину при большой интенсивности (при интенсивности I = 50 мм/ч для длины волны Л = 3, 2 см коэффициент оставления с оставляет аос = 1, 244 дБ, ир и Л = 5, 6 см коэффициент ослабления равен аос = 0, 214 дБ и ир и Л = 10 см коэффициент ослабления имеет значение аос = 0,15 дБ).
Таким образом, при использовании сигнала с длиной волны Л = 3, 2 см и интенсивности дождя больше 10 мм/ч потери энергии на дальности 1 км составляют более 3,16 дБ.
Т а б .л и ц а 1
Параметры а и Ь
А(см) 0,9 3,2 10
а 0,22 0,0074 0,003
Ь 1 1,31 1
Рис. 3. Зависимость коэффициента ослабления сигнала в дожде от интенсивности при различной длине волны
Ослабление сантиметровых радиоволн в виде снега
Снежинки, составляющие снегопад, имеют различную форму, ориентацию, плотность, размеры и содержание влаги. Большинство снежинок имеет гексагональную плоскую форму и диаметр в диапазоне от 0,5 до 3 мм. Показатель преломления представляет собой суммарный показатель для смеси воздуха, воды и льда. Учитывая вышесказанное, расчет коэффициента ослабления сигнала для снега является непростой задачей. При предположении, что снежинка имеет сферическую форму и ее масса равна массе такой же водяной сферы, можно рассчитать коэффициент ослабления жидкого снега через интенсивность снегопада I по формуле [1|:
а, = 0, 094(22) Л
Для сухого снегопада при 0 °С коэффициент ослабления описывается формулой Баттана [7]:
I2 I
ап = 0, 035— + 0, 0022-. (23)
Л Л
На рис. 4 представлены графики зависимости коэффициента ослабления от интенсивности при различных типах снега для нескольких длин волн:
Сравнение графиков (рис. 4) для интенсивности, изменяющейся в диапазоне от 0,1 до 10 мм/ч, показывает, что коэффициент ослабления для жидкого снега значительно больше коэффициента ослабления для сухого снега, примерно в 40 раз. Средняя интенсивность сухого снегопада не превышает нескольких мм/ч, поэтому в сантиметровом диапазоне при метеорологическом наблюдении ослаблением сигнала в сухом снеге можно пренебречь.
Интенсивность!, мм/ч Интенсивность!, мм/ч
Рис. 4. Зависимость коэффициентов ослабления сигнала дожди от интенсивности
Из графика видно, что коэффициент ослабления для жидкого снегопада в Х-диапазоне линейно зависит от интенсивности I. Для определения коэффициента ослабления используем формулу (22) или график, показанный на рисунке 4.
На действующем макете ПРЛ Х-диапазона с параллельным обзором был проведен ряд экспериментов но работе но метеобразованием.
Рис. 5 демонстрирует, что на расстоянии 11 км от ПРЛ было обнаружено дождевое облако протяженностью ¿0бл = 4 км [5].
а)
б)
Рис. 5. а) Зависимость нормированного уровня принимаемого сигнала от расстояния, б) карта Погода Vontusky.com
На рис. 5а представлены результаты первичной обработки, показан уровень принимаемого сигнала в зависимости от расстояния. А на рисунке 56 отображена карта метеообразований «Ventusky.com» на момент измерений.
Рассмотрим данный пример, в котором интмрирование но формуле (4) заменим простым суммированием. Для примера рассчитаем на длине волны 3,2 см двухстороннее общее ослабление энергии радиоволн для расстояния 15 км, вызванное их поглощенном и рассеянием, если траектория радиоволн на отрезке 4 км проходит в водяных облаках с интенсивностью и температурой 0 °С.
Ослабление в газах атмосферы:
Г = 2агДгг = 2(0, 0072 + 0, 0039)15 = 0, 33[дБ]. (24)
Ослабление в облаках с осадками:
Гд = 2адДгд = 2.0, 045.4 = 0, 36[дБ]. (25)
В результате общее двухстороннее ослабление энергии на трассе при наблюдениях с метеорологическими целями в 11 часов 11-го июля 2019 года
~Р
10^= = -(0, 33 + 0, 36) = -0, 69[дБ], (26)
Ро
На рисунках 6а и 66 продемонстрировано, что ПРЛ расположен в зоне енмхшада протяженностью 4 км.
Рис. 6. Зависимость нормированного уровеня принимаемого сигнала от расстояния а. карта «Яндекс Погода» б
В данном случае имеем следующие параметры для расчета: дальность 4 км, протяженность 4 км, енмхшад со сродней интенсивностью 1,9 мм/ч при температуре 5 °С.
Ослабление в газах атмосферы:
Г = 2агДгг = 2(0, 0072 + 0, 0039)4 = 0, 088[дБ]. (27)
Ослабление в облаках с осадками:
Гсне = 2аСнеДгсне = 2.0, 055.4 = 0, 44[дБ]. (28)
Таким образом, общее ослабление энергии при наблюдениях с метеорологическими целями в 16 часов 6-1'о февраля 2019 года составляет:
101д= = -(0, 088 + 0, 44) = -0, 528[дБ]. (29)
Ро
На рассмотренных численных примерах продемонстрировано, что при работе на малой дальности при небольшой интенсивности метеорологических явлений ослабленном сигналов можно пренебречь.
2. Заключение
Основные выводы и результаты, представленные в статье:
- величина коэффициента ослабления в миллиметровом диапазоне значительно больше, чем в сантиметровом диапазоне;
- при на работе на малых дальностях ослаблением в газах и облаках без осадков можно пренебречь:
- для рассматриваемого ПРЛ сантиметрового диапазона (рабочая длина волны Л = 3, 2 см) необходимо учитывать затухание сигнала только при сильном дожде, граде и снегопаде. Например, при сильном дожде интенсивностью I = 50 мм/ч на дальности 10 км сигнал ослабляется на 3,16 dB, а при сильном снегопаде I = 10 мм/ч на дальности 10 км потери энергии составляют 5,8 dB.
В данной работе представлена методика расчета общих потерь энергии на трассе распространения радиоволн. Для основных факторов, оказывающих влияние на ослабление сигнала, приведены аналитические выражения для расчета соответствующих коэффициентов ослабления.
Литература
1. Брылёв Г.Б., Гтшта С.Б., Низдойминого, Г.Л. Радиолокационные характеристики облаков и осадков. Ленинград : Гидрометеоиздат, 1986. 234 с.
2. Грудипская Г.П.Распространение радиоволн: учеб. пособие. Москва : Высшая школа, 1975. 280 с.
3. Довиак Р., Зрпич Д. Доплеровские радиолокаторы и метеорологические наблюдения. Ленинград : Гидрометеоиздат, 1988. 512 с.
4. Кхыонг Н.В. Экспериментальная проверка способности обнаружения метеообразований новым посадочным радиолокатором // Труды МФТИ. 2020. Т. 12, № 1. С. 44-52.
5. Пожидаев В.Н. Ослабление и обратное рассеяние миллиметровых радиоволн в тумане, дожде, снегопаде // IV Всероссийская конференция «Радиолокация и радиосвязь». ИРЭ РАН, 29 ноября-3 декабря 2010 г.
6. Родос Л.Я. Электродинамика и распространение радиоволн (распространение радиоволн). Санкт-Петербург : Изд-во СЗТУ, 2007. 89 с.
7. Степаненко В.Д. Радилокация в метеорологии. Ленинград : Гидрометеоиздат, 1966. 350 с.
8. Черный Ф.Б. Распространение радиоволн. Москва : Советское радио, 1972. 463 с. References
1. Brylyov G.B., Gashina S.B., Nizdoyminoga G.L. Radar characteristics of clouds and precipitation. Leningrad : Gidrometeoizdat, 1986. 234 p.
2. Grudinskaya G.P. Propagation of radio waves: textbook, allowance. Moscow : Higher School, 1975. 280 p.
3. Doviak R., Zrnich D. Doppler radars and meteorological observations. Leningrad : Gidrometeoizdat, 1988. 512 p.
4. Khuong N. V. An experimental test of the ability to detect meteorological conditions with a new landing radar. Proceedings of MIPT. 2020. V. 12, N 1. P. 44-52.
5. Pozhidaev V.N. Attenuation and backscattering of millimeter-wave radio waves in fog, rain, snowfall. IV Ail-Russian Conference «Radar and Radio Communication». IRE RAS, November 29-December 3, 2010.
6. Rhodes L.I. Electrodynamics and radio wave propagation (radio wave propagation). Sankt-Peterburg : Izd-vo SZTU, 2007. 89 p.
7. Stepanenko V.D. Radar in meteorology. Leningrad : Gidrometeoizdat, 1966. 350 p.
8. Cherny F.B. Propagation of radio waves. Moscow : Soviet Radio, 1972. 463 p.
Поступим в редакцию 24-06.2020
