Влияние тропосферной рефракции на распространение радиоволн

жүктеу/скачать 152.92 Kb.

Дата	03.03.2016
өлшемі	152.92 Kb.
	#36381

Рек. МСЭ-R P.834-6

РЕКОМЕНДАЦИЯ МСЭ-R P.834-6^{^*}

Влияние тропосферной рефракции на распространение радиоволн

(Вопрос МСЭ-R 201/3)

(1992-1994-1997-1999-2003-2005-2007)

Сфера применения

В Рекомендации МСЭ-R P.834 предоставляются методы расчета крупномасштабных рефракционных явлений в атмосфере, включая отклонение луча, волноводные слои, эквивалентный радиус Земли, видимый угол места и углы направления максимума излучения на трассах Земля-космос и эффективную длину радиотрассы.

Ассамблея радиосвязи МСЭ,

учитывая,

a) что для надлежащего планирования наземных линий связи и линий связи Земля-космос необходимы соответствующие процедуры расчетов для оценки влияния рефракции на радиосигналы;

b) что были разработаны процедуры, позволяющие количественно оценить некоторые рефракционные эффекты при распространении радиосигналов по наземным линиям связи и линиям связи Земля-космос,

рекомендует,

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

Приложение 1

1 Отклонение луча

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

, (1)

где:

 : радиус кривизны траектории луча;

n : индекс рефракции в атмосфере;

dn/dh : вертикальный градиент индекса рефракции;

h : высота точки над поверхностью Земли;

 : угол траектории луча относительно горизонтальной плоскости в рассматриваемой точке.

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

2 Эквивалентный радиус Земли

Если трасса является почти горизонтальной, величина  близка к нулю. Однако так как n очень близок к единице, уравнение (1) упрощается следующим образом:

. (2)

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

Если высотный профиль изменения рефракции линеен, т. е. градиент рефракции остается постоянным вдоль траектории луча, возможно преобразование, которое позволяет рассматривать распространение как прямолинейное. Это преобразование предполагает, что распространение происходит над гипотетической Землей с эквивалентным радиусом R_e  ka при:

, (3)

где a – фактический радиус Земли, а k – коэффициент эквивалентного радиуса Земли (k-фактор). В результате такого геометрического преобразования траектории лучей становятся линейными независимо от угла места.

Строго говоря, градиент рефракции остается постоянным, только если трасса горизонтальна. На практике для высот менее 1000 м экспоненциальную модель усредненного профиля индекса рефракции (см. Рекомендацию МСЭ-R P.453) можно заменить линейной. Соответствующий k фактор k  4/3.

3 Модифицированный индекс рефракции

Для некоторых применений, например при построении траектории луча, используется модифицированный индекс рефракции или модуль рефракции, который определен в Рекомендации МСЭ-R P.310. Модуль рефракции, M, определяется уравнением:

, (4)

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

4 Видимый угол направления максимума излучения на наклонных трассах

4.1 Введение

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

4.2 Видимость космической станции

Как описано выше в п. 1, радиолуч, излучаемый станцией, находящейся на поверхности Земли (на высоте h (км) и имеющей угол места q (градусов)), будет под влиянием атмосферной рефракции изгибаться в сторону Земли. Поправка на рефракцию, t (в градусах), вычисляется с помощью следующего интеграла:

, (5)

где j определяется в полярных координатах на основе закона Снеллиуса следующим образом:

(6)

(7)

r : радиус Земли (6370 км);

x : высота (км).

Поскольку отклонение луча в значительной степени определяется нижней частью атмосферы, для обычной атмосферы индекс рефракции на высоте x можно определить по формуле:

, (8)

где:

a  0,000315;

b  0,1361.

Эта модель основана на использовании экспоненциальной модели атмосферы для распространения по наземным трассам, которая приведена в Рекомендации МСЭ-R P.453. Кроме того, n' (x) является производной от n(x), т. е. n' (x)  –ab exp (–bx).

Значения t (h, ) (в градусах) рассчитывались для условий эталонной атмосферы, при этом оказалось, что следующая числовая формула обеспечивает хорошую аппроксимацию:

t (h, q)  1/[1,314  0,6437 q  0,02869 q²  h (0,2305  0,09428 q  0,01096 q²)  0,008583 h²]. (9)

Приведенная выше формула была получена в качестве аппроксимации для 0 £ h £ 3 км и
q_m £ q £ 10, где _m – угол, под которым радиолуч пересекается с земной поверхностью
и который определяется как:

(10)

или, приблизительно, (градусы).

Уравнение (9) дает также вполне приемлемую аппроксимацию в диапазоне 10  q £ 90.

Пусть угол места космической станции в условиях распространения в свободном пространстве будет равен ₀ (градусы), а минимальный угол места от станции на поверхности Земли, для которого радиолуч не пересекается с земной поверхностью, – _m. Поправка на рефракцию, соответствующая _m, будет равна t (h, _m). Следовательно, космическая станция будет видима лишь в том случае, если выполняется неравенство:

. (11)

4.3 Определение видимого угла места

Если неравенство в уравнении (11) выполняется, то видимый угол места,  (градусы), можно определить с учетом атмосферной рефракции, решив следующее уравнение:

. (12)

Решение уравнения (12) имеет вид:

, (13)

где значения t_s (h, q₀) идентичны значениям t (h, q), но выражены в виде функции q₀.

Функция t_s (h, q₀) (градусы) хорошо аппроксимируется следующей числовой формулой:

t_s (h, q₀)  1/[1,728  0,5411 q₀  0,03723 q₀²  h (0,1815  0,06272 q₀

 0,01380 q₀²)  h² (0,01727  0,008288 q₀)]. (14)

Значение , рассчитанное с помощью уравнения (13), и будет видимым углом места.

4.4 Основные этапы вычислений

Этап 1: Угол места космической станции в условиях распространения в свободном пространстве обозначается как q₀.

Этап 2: С помощью уравнений (9) и (10) проверить, выполняется ли уравнение (11). Если оно не выполняется, то спутник невидим, и, следовательно, дальнейшие вычисления не требуются.

Этап 3: Если ответ на этапе 2 положительный, то с помощью уравнений (13) и (14) рассчитывается q.

4.5 Результаты измерения видимого угла направления максимума излучения

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

ТАБЛИЦА 1

Величины угловой девиации при распространении радиоволн через атмосферу в целом

Угол места,  (градусы)	Средняя общая угловая девиация,  (градусы)
Угол места,  (градусы)	Полярный континентальный воздух	Умеренный континентальный воздух	Умеренный морской воздух	Тропический морской воздух
1 2 4 10 20 30	0,45 0,32 0,21 0,10	– 0,36 0,25 0,11 0,05 0,03	– 0,38 0,26 0,12 0,06 0,04	0,65 0,47 0,27 0,14
	Колебания величины  ото дня ко дню (только для столбцов 1 и 4)
1 10	0,1 среднеквадратичное значение 0,007 среднеквадратичное значение

5 Фокусировка и расфокусировка волны при распространении через атмосферу

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

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

где:

₀: угол места линии, соединяющей передающий и приемный пункты, в градусах (₀< 10);

h: высота самой низкой точки над уровнем моря, в км (h < 3 км);

b: изменение уровня сигнала для волны, проходящей через атмосферу, по сравнению с условиями свободного пространства, в децибелах;

знак b в уравнении будет отрицательным ("–") для источника передачи, расположенного вблизи поверхности Земли, и положительным ("+") для источника, расположенного за пределами атмосферы.

6 Эффективная длина радиотрассы и ее изменения

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

, (15)

где:

s : расстояние вдоль трассы;

n : индекс рефракции;

A и B : концы трассы.

Уравнение (15) можно использовать лишь в том случае, когда известны изменения индекса рефракции n вдоль трассы.

Если температура, T, атмосферное давление, P, и относительная влажность, H, на уровне поверхности Земли известны, то увеличение длины трассы L будет вычисляться с помощью полуэмпирического метода, о котором речь пойдет ниже и который был разработан с использованием профилей, полученных с помощью радиозондирования атмосферы, на основе данных, собранных на 500 метеорологических станциях в течение одногодичной кампании 1979 года. В основе этого метода лежит следующее уравнение для определения увеличения длины трассы L:

, (16)

где:

₀ :угол места в точке наблюдения;

L_V :увеличение длины трассы по вертикали;

k и  (₀, L_V) : поправочные коэффициенты, при вычислении которых используется экспоненциальная модель атмосферы.

Изменения угла места на трассе учитываются с помощью коэффициента k. Член  (₀, L_V) отражает влияние рефракции (трасса не является прямой линией). Этот член всегда очень мал по величине, за исключением тех случаев, когда угол места очень мал, а потому им можно пренебречь в вычислениях; это приводит к погрешности всего лишь 3,5 см для угла ₀ порядка 10 и 0,1 мм для угла при ₀ порядка 45. Кроме того, можно заметить, что при очень низких углах места, когда нельзя пренебречь членом , допущение относительно плоской стратифицированной атмосферы, лежащее в основе всех методов расчета увеличения длины трассы, более недействительно.

Увеличение длины трассы по вертикали (м) описывается уравнением:

L_V  0,00227 P  f (T ) H. (17)

В первом члене правой части уравнения (17) P – атмосферное давление (гПа) в точке наблюдения.

Во втором члене, полученным эмпирическим путем, H – относительная влажность (%); функция температуры f (T ) зависит от географического местоположения и имеет следующий вид:

f (T )  a 10^bT , (18)

где:

T выражена в C;

a выражена в м/% относительной влажности;

b выражена в C^–1.

Параметры a и b в зависимости от географического местоположения приведены в таблице 2.

ТАБЛИЦА 2

Местоположение	a (м/%)	b (C^–1)
Прибрежные районы (острова или пункты, удаленные менее чем на 10 км от берега моря)	5,5  10^–4	2,91  10^–2
Неприбрежные экваториальные районы	6,5  10^–4	2,73  10^–2
Все прочие районы	7,3  10^–4	2,35  10^–2

При вычислении поправочного коэффициента k, входящего в уравнение (16), делается предположение об экспоненциальном законе изменения атмосферной рефракции N в зависимости от высоты h:

N(h)  N_s exp (– h / h₀), (19)

где N_s – среднее значение рефракции у поверхности Земли (см. Рекомендацию МСЭ-R P.453), а h₀ определяется как:

(20)

Тогда k вычисляется с помощью следующего выражения:

, (21)

где n_s и n (h₀) – значения индекса рефракции вблизи поверхности Земли и на высоте h₀ (определяемые из уравнения (20)), соответственно, а r_sи r (h₀) – соответствующие расстояния до центра Земли.

Для трасс Земля спутник с углом места q более 10° увеличение длины трассы в тропосфере (м) можно представить в виде суммы "сухой" и "влажной" составляющих:

м, (22)

где:

p_s, e_s: общее давление воздуха и давление водяного пара у поверхности Земли (гПа);

T_ms: средняя температура водяного пара над поверхностью (K);

λ: коэффициент уменьшения давления пара;

R_d: R/M_d = 287,0 (Дж/кг K);

R: молярная постоянная газа = 8,314 (Дж/мол K);

M_d: молярная масса сухого воздуха = 28,9644 (г/мол);

k₁ = 77,604 (K/гПа);

k₂ = 377600 (K²/гПа);

g_m = g(lat,h_s) = 9,784 · (1 – 0,00266 · cos (2 · lat) – 0,00028 · h_s) = ускорение силы тяжести (м/с²);

lat: широта местоположения (град.);

h_s: высота местоположения над уровнем моря (км).

Для приемников, расположенных на высоте над средним уровнем моря, h (км), отличном от высоты поверхности, h_s, значения входных метеорологических параметров могут быть экстраполированы из значений поверхности, T_ms, e_s и p_s, используя следующие уравнения:

м (23a)

K (23b)

гПа (23c)

гПа, (23d)

где:

a_m : вертикальный градиент средней температуры водяного пара (K/км).

=вертикальный градиент температуры воздуха K/км (23e)

= R_d /1 000 = 0,287 Дж/г K (23f)

K (23g)

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

, (24)

где:

X_i: p_s, e_s, T_ms, l или am;

a1_i: среднее значение параметра;

a2_i: сезонная флуктуация параметра;

a3_i: день минимального значения параметра;

D_y: день года (1 ... 365,25),
1 = 1 января, 32 = 1 февраля, 60,25 = 1 марта.

Коэффициенты a1, a2 и a3 были вычислены путем подгонки гармонической модели к месячной статистике метеорологических параметров, взятой из базы данных ECMWF ERA15.

Значения параметров этой гармонической модели, a₁, a₂, a₃ для каждого метеорологического параметра, содержатся в цифровых картах, общим числом 15 карт, с разрешением 1,5 × 1,5 (град.). Коэффициенты гармонической модели содержатся в файлах под названием <met_param>_<mod_par>.dat. При этом <met_param> представляет собой press, vapr, tmpm, lamd или alfm для p_s, e_s, T_ms, l или am, соответственно, а <mod_par> – это a1, a2 или a3. Данные представлены в диапазоне от 0 до 360 по долготе и от +90 до –90 по широте, с разрешением 1,5° как по широте, так и по долготе. Дополнительно необходимые данные с помощью модели следующие: lat.dat = эталонные значения широты (90 ... –90), long.dat = эталонные значения долготы (0 ... 360), hreflev.dat = средняя высота опорного уровня по отношению к уровню моря (м).

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

Точность предлагаемой модели была проверена с использованием радиозонда и измерений в рамках системы GPS, и полученные результаты оказались в пределах 2 и 6 см по всему миру. Если необходима более высокая точность, то в качестве входных данных в модели могут использоваться совместные локальные измерения общего давления воздуха и давления водяных паров.

7 Распространение при наличии волноводных слоев

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

Существование волноводов является очень важным фактором, поскольку они могут способствовать возникновению аномальных условий распространения, особенно на наземных линиях связи или линиях связи Земля космос с очень небольшими углами места. Волноводы создают такой механизм распространения радиосигналов достаточно высокой частоты далеко за пределы их обычной прямой видимости, что является причиной возможных помех для других служб (см. Рекомендацию МСЭ R P.452). Кроме того, они играют важную роль с точки зрения возникновения многолучевой интерференции (см. Рекомендацию МСЭ-R P.530), хотя их наличие не является ни необходимым, ни достаточным условием для возникновения многолучевого распространения на отдельной линии.

7.1 Влияние угла места

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

, (25)

_{где d}_M_/d_h_{– вертикальный градиент модифицированной рефракции}

,_а__h_{– толщина}волновода, которая представляет собой высоту его верхней границы над передающей антенной.

На рисунке 2 показан максимальный угол места, при котором лучи могут быть захвачены волноводом. Максимальный угол захвата быстро увеличивается при уменьшении градиентов рефракции ниже уровня –157 N/км (т. е. при увеличении вертикальных градиентов) и при увеличении толщины волновода.

7.2 Минимальная частота захвата

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

Минимальную частоту захвата сигнала тропосферным волноводом можно рассчитать с помощью фазового интеграла. На рисунке 3 представлены минимальные значения частоты захвата для приземных волноводов (сплошные кривые), где допускается постоянный (отрицательный) градиент рефракции от поверхности Земли до заданной высоты и со стандартным профилем рефракции выше этой высоты. Для частот, используемых наземными системами (обычно 8–16 ГГц), минимальная толщина волноводного слоя должна быть приблизительно от 5 до 15 м, и в этих случаях минимальная частота захвата, f_min, в значительной мере зависит как от толщины волновода, так и от градиента индекса рефракции.

В случае приподнятых волноводов, даже когда рассматривается простой линейный профиль рефракции, вводится дополнительный параметр. Этот параметр описывает форму профиля индекса рефракции подлежащего слоя. На рисунке 3 штриховые линии показывают минимальную частоту захвата для волноводного слоя с постоянным градиентом, расположенного над поверхностным слоем со стандартным градиентом рефракции, равным –40 N/км.

РИСУНОК 2

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

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

РИСУНОК 3

Минимальная частота захвата сигнала в атмосферные радиоволны
с постоянными градиентами рефракции

___________________

*3 я Исследовательская комиссия по радиосвязи внесла редакционные поправки в настоящую Рекомендацию в 2009 году в соответствии с Резолюцией 1 МСЭ R.

жүктеу/скачать 152.92 Kb.

Достарыңызбен бөлісу: