
книги из ГПНТБ / Проворов К.Л. Радиогеодезия учеб. пособие
.pdfX I I Генеральная ассамблея МГГС в 1960 г. рекомендовала для индекса показателя преломления воздуха следующие формулы:
1. Для светового диапазона формулу Баррела и Сирса
|
- Л к - |
Г 0,055 |
(63) |
N |
' i + at е, |
||
|
г р ~ l + at 760 |
|
где ІѴ?rrp— индекс преломления сухого воздуха находится по формуле
(53) с коэффициентами |
А, |
Ѵ5 и |
С, |
приведенными |
в табл. 6. |
2. Для диапазона |
ультракоротких радиоволн |
формулу Фрума |
|||
и Эссена |
|
|
|
|
|
103,49 |
(Р-е) |
+ |
86,27 |
(64) |
|
. 0 |
20 |
-20 |
0 20 |
W |
|
Температура, "С |
. Температура, "С |
|||
|
|
Р и с . |
20 |
|
|
Обе формулы |
можно |
получить |
из формул |
(59) и (61), положив |
|
в формуле (59) Т |
= 273,16(1 + at), |
где а = 2 |
? g ^ — коэффициент |
объемного температурного расширения воздуха.
На рис. 20, а приведен график зависимости показателя преломле ния воздуха от температуры и влажности, рассчитанный по формуле (62) при давлении 1013 мб для радиоволн. На рис. 20, б приведен аналогичный график для световых волн, рассчитанный по формуле (60) при Nrp= 304,5. Графики показывают значительно меньшую зависимость индекса преломления световых волн от влажности воздуха.
Так как плотность атмосферы, определяющая показатель прелом ления, в общем плавно убывает с высотой, то и показатель преломле ния с высотой также уменьшается. Зависимость индекса преломления
от высоты можно представить в виде |
|
N=A + BH+C№, |
(65) |
где коэффициенты А, В я С для некоторой территории можно |
найти |
по данным многолетних наблюдений давления, температуры |
и влаж |
ности воздуха на разных высотах. Имея такие наблюдения, |
можно |
в пределах некоторых равноотстоящих значений высот Hlf |
Н2, . . . |
50
Hп найти для каждой высоты соответствующую среднюю величину Ni, после чего составить уравнения погрешностей вида
Vt = A+BHt + CHt — Nt,
а величины А, В и С получить из соответствующих нормальных уравнений. Для стандартной атмосферы с постоянной относительной влажностью 70% индекс преломления для радиоволн достаточно х о р о ш о представляется формулой
ІѴ = 3 2 9 - 4 0 , 4 Я + 1 , 6 9 # 2 , |
(66) |
где //^выражено в километрах. На рис. 21, а показана зависимость индекса преломления от высоты, рассчитанная по формуле (66). На рис. 21, б показана аналогичная зависимость для светового диа пазона.
6
лп |
лп |
1 |
• • I |
I I |
I I I I I |
I |
0 2 Ч 6 8 10 0 ) 2 3 4 5 |
||||
• |
Высота, км |
Высота, км |
|
|
|
|
Рис . |
21 |
|
Более точно величину |
индекса |
преломления |
для высот до 9 км |
с ошибкой около 5 единиц можно получить по номограмме, составлен ной А. И. Соловьевым, Р. М. Гришечевой и Г. П. Гричуком (рис. 22).. По известному индексу преломления N0 в какой-либо точке и ее. высоте Н0 по номограмме находят соответствующую кривую, на кото рой и получают затем искомую величину N для высоты Н. Номограмма построена по данным многолетних наблюдений метеорологических элементов на 38 пунктах, равномерно расположенных на террито рии СССР.
Из формулы (66) найдем вертикальный градиент индекса |
прелом |
|||
ления |
|
|
|
|
|
gN = ~ |
= B + 2CH. |
(67) |
|
Для стандартной |
атмосферы |
из (67) |
получим |
|
^ |
= ( _ 40,4 + |
3,38Я) |
. Ю - 3 (на 1 м), |
(68) |
где H выражено в километрах. Отсюда видно, что модуль градиента убывает с высотой. В приземном слое при H = 0 градиент индекса
4* |
51 |
преломления в стандартной атмосфере составляет около — 4 - Ю - 2 на 1 м. Так как кривые на рис. 21 в общем мало уклоняются от пря
мых, то для приближенных расчетов можно |
принимать для некото |
||||||
рого слоя тропосферы |
между |
высотами |
Нг |
и Н2 |
градиент |
индекса |
|
преломления равным |
среднеинтегральному |
значению |
|
||||
gN- Я |
1 Я |
j ( в |
+ 2CH)dH |
= B + |
2СНт. |
(69) |
|
2 |
1 |
я , |
|
|
|
|
|
Н(м)
9000
8000
7000
6000
5000,
200 |
250 |
350 |
380 |
Р и с . 22
Для стандартной атмосферы на высоте от 0 до 1 км средний градиент индекса преломления, найденный по формуле (69), соста вляет
gN = ( - 40, 6 + 1 , 6 9 - 1 ) - Ю - 3 = |
- 3 , 9 - Ю - 2 (на 1 м). |
||
Если известны значения показателя преломления пх и п2 на высо |
|||
тах соответственно |
І7_ и Н2, то средний градиент для этого слоя можно |
||
найти по формуле |
|
|
|
|
gN- п2 — "1 |
(70) |
|
|
На |
HЛ |
|
Приближенной |
величиной показателя преломления, |
найденной |
|
по эмпирическим формулам или на основании различных |
допущений |
о градиенте показателя, можно пользоваться при радиогеодезических измерениях сравнительно невысокой точности, а также при общих исследованиях условий и точности измерений. Точные радио- и светодальномерные измерения выполняются в основном в приземном слое воздуха высотой до 100 м над подстилающей поверхностью, где про-
52
исходят наибольшие изменения и возмущения метеорологических элементов. В этих условиях показатель преломления также подвер жен наибольшим систематическим и случайным как суточным, так и сезонным колебаниям. Так, в течение суток индекс преломления в одном и том же месте может измениться на величину до 50 единиц, а флуктуации его в неблагоприятных условиях достигают 5 единиц. По результатам специальных измерений (Шонеманн и Билл), выпол ненных с применением рефрактометра, среднее за 42 дня значение флуктуации индекса преломления на высотах до 79 м получилось 0,53, а число флуктуации в минуту составило в среднем 4,8. В отдель ных случаях зафиксированы колебания индекса преломления вели чиной порядка 3—5. Амплитуды и частоты флуктуации индекса
преломления наблюдались |
максимальными в приземном слое (до |
1 м); с увеличением высоты |
они быстро убывают. |
Изменение индекса показателя преломления воздуха с высотой часто происходит скачками, а его градиент в слое высотой до 1 км колеблется от —0,04 до —0,30 на 1 м. При неблагоприятных мете орологических условиях (температурная инверсия, резкое умень шение влажности воздуха с высотой) на отдельных участках тропо сферы возникают аномалии градиента, который принимает нулевое значение и даже становится положительным. При этом показатель преломления вместо обычного уменьшения с высотой возрастает. Колебания показателя преломления воздуха в общем уменьшаются с увеличением высоты и на высоте около 9 км его значение практи чески постоянно в течение всего года.
Еще меньшая закономерность изменения показателя преломления воздуха наблюдается вдоль поверхности земли, в особенности в приземном слое атмосферы. Известны случаи изменения индекса показателя преломления на малых высотах на расстоянии несколь ких километров до 50 единиц. Характер изменений в большой степени зависит от подстилающей поверхности и рельефа местности. Так, над лесом величина вертикального градиента индекса преломления на блюдалась в 6—8 раз большая, чем над открытой местностью в одно и то же время. С увеличением высоты колебания показателя прелом ления в горизонтальном направлении сглаживаются.
Наибольшие изменения показателя преломления происходят в ясную теплую погоду и в особенности при заходе и восходе Солнца. При осадках, тумане, ветре, в пасмурную погоду изменения показа теля преломления выражены менее резко. Наименьшие изменения индекса преломления чаще бывают во вторую половину дня и наи большие — ночью. Наибольшие флуктуации обычно наблюдаются днем. На рис. 23 показан суточный ход градиента индекса преломле
ния по данным |
за март — сентябрь в Ленинградской |
области. |
На |
рис. 24 показан |
суточный ход градиента в течение |
полного |
года |
по данным трехлетних наблюдений в Англии (штриховкой на рис. 24 показано ночное время).
Из изложенного видно, что при точных радио- и светодальномерных измерениях показатель преломления следует находить по
53
значениям температуры, давления и влажности для момента изме рений, а также выбирать для измерения расстояний наиболее благо приятное время, избегая по возможности периодов с резкими коле баниями метеорологических элементов.
Влияние ошибок метеорологических факторов на точность опре деления показателя преломления нетрудно найти путем дифференци рования формул (60) и (62). Для светового
диапазона получим
dNrp = |
0.270 |
(N%p-ASe)dT + |
f2 |
Соответственно для радиоволн
dN = — 77,6 |
е |
9582e |
77,6 dp- |
|
' б |
"г" Т |
|
, |
77,6 / |
4791 |
1 \ , |
оарь |
Март |
Май |
Июль |
|
Сентябрь Ноябрь |
Февраль |
|
Апрель |
Июнь |
Август |
Октябрь Декабрь |
|
|
|
Рис . 24 |
|
|
Приняв для коэффициентов средние условия, т. е. при А^р— 305;
Т = 288° К; р = |
1013 мб и е — 10 мб, получим следующие |
выра |
||
жения: |
|
|
|
|
для |
световых |
волн |
|
|
для |
радиоволн |
dNrp = 0,3dp-l,0dT |
— 0,04de; |
(71) |
|
|
|
||
|
|
dN = 0,3dp-l,3dT |
+ 4,5de; |
(71)' |
54
здесь dp и de выражены в миллибарах и dT — в градусах. Соответ ствующие средние квадратические ошибки индекса преломления будут:
для световых волн
mN |
= V (0,3mp )2 + (1,0m.)2 + (0,04m,)2; |
(72) |
|
для радиоволн |
|
|
|
mN |
= V(0,3mp )2 + (l,3m,)2 |
+ (4,5me )2 . |
(72)' |
Экспериментальными исследованиями |
установлено, что ошибки |
измерения метеорологических элементов зависят от погоды; при этом
с повышением температуры и в особенности влажности |
воздуха |
||||||||||
ошибки mt и те |
возрастают. При наиболее |
благоприятных |
условиях |
||||||||
(при температуре около |
0° С ночью |
и в |
пасмурную |
погоду днем) |
|||||||
измеренные |
температура |
и |
влажность получаются, |
как |
правило, |
||||||
с наименьшими |
ошибками: |
mt — 0,1—0,2° С; те =0,2—0,3 мб. |
|||||||||
При высокой температуре и влажности, в особенности |
в ясную по |
||||||||||
году днем, |
ошибка измерения температуры |
достигает |
0,6—0,8° С, |
||||||||
а ошибка |
измерения |
влажности — 1,3—1,5 мб. Положив |
средние |
||||||||
квадратические |
ошибки |
метеорологических |
элементов |
равными: |
|||||||
тр = ± 0,7 мб; mt = |
± 0 , 4 ° С; те = |
± 0 , 7 мб, найдем по формулам |
(72) средние квадратические ошибки индекса преломления для сред них условий:
для |
световых волн |
|
|
|
mN |
= 1/0,04 + 0,16 + 0,00 = |
± 0 , 4 ; |
для |
радиоволн |
|
|
|
mN |
= 1/0,04 + 0,31 + 9,92 = |
± 3,2. |
При неблагоприятных условиях ошибки могут увеличиться в 2— 3-раза. Как видно из расчета, основное влияние на точность индекса преломления световых волн оказывает ошибка измерения темпера туры воздуха. Точность измерения влажности в этом случае почти
не |
имеет значения. Наоборот, для радиодиапазона главное влияние |
на |
точность индекса преломления оказывает ошибка определения |
влажности.
Существуют способы непосредственного определения показателя преломления воздуха с помощью рефрактометра. Один из способов основан на сравнении частот двух объемных резонаторов — эталон ного, вакуумированного и герметизированного и измерительного, сообщающегося с воздухом. Индекс преломления в этом случае
равен относительной разности |
частот резонаторов, |
т. е. |
_Ѵ = |
10е. |
(73) |
/изм Индекс преломления, пропорциональный разности частот из
мерительного и эталонного резонаторов, непрерывно отмечается
55
с помощью самописца в виде кривой, характеризующей ход его изме нения.
Способы непосредственного определения показателя преломле ния позволяют получить его значение в какой-либо точке или в ряде отдельных точек. Определение среднего показателя преломления вдоль измеряемой линии рассмотрено в § 7.
§ 6. ФОРМА П У Т И Э Л Е К Т Р О М А Г Н И Т Н Ы Х В О Л Н
При обработке результатов радиоэлектронных измерений суще ственное значение имеет форма пути, по которому распространяются электромагнитные волны. В свободном пространстве или в оптически однородной среде эти волны распространяются прямолинейно. В не однородной среде электромагнитные волны, согласно принципу
Ферма, распространяются по пути, оп тическая длина которого (произведение геометрической длины пути Dr на пока затель преломления п) наименьшая.
Вследствие неоднородности атмосферы электромагнитные волны распространя ются по некоторой сложной кривой. По казатель преломления воздуха в общем уменьшается с высотой, поэтому путь волны при переходе ее из нижележащего слоя в соседний будет удаляться от на правления отвесной линии (рис. 25) и соответствующая кривая будет обращена вогнутостью к Земле. Искривление пути волны в атмосфере называют атмосферной рефракцией.
Проекцию пути волны на поверхность земного эллипсоида счи тают прямым нормальным сечением, проведенным из начальной точки. Кроме того, считают, что ввиду малого времени прохождения волны в прямом и обратном направлениях, как это имеет место при радиоэлектронных измерениях расстояний, соответствующие пути волн совпадают.
Представим атмосферу разделенной на слои, параллельные уровенным поверхностям, настолько тонкие, что показатель преломле
ния |
в каждом слое |
можно |
принять |
постоянным (см. рис. 25). На |
||||
границе двух |
таких |
слоев, |
согласно |
закону |
Снеллиуса, |
|
||
|
|
|
п sin ô = щ sin ß t |
= . .. = |
const, |
(74) |
||
где |
ô — у г о л |
падения, a |
— угол |
преломления луча. |
|
|||
|
Обозначив радиусы соответствующих поверхностей через R, |
Ях, |
||||||
Л 2 и т. д., напишем |
|
|
|
|
|
|
||
|
|
|
|
Д |
sm ô |
|
(75) |
|
|
|
|
|
|
sin |
ß |
|
|
|
|
|
|
|
|
|
56
|
Исключая из (74) и (75) угол б, получим |
|
|
|
||||||||||
|
|
|
|
Rnsin ß = RxHx sin ßj = |
. . . = const. |
|
|
|
||||||
|
Или |
в общем |
виде |
Rn cos а = const, |
|
|
|
(76) |
||||||
|
|
|
|
|
|
|
|
|||||||
где а — 90° — ß — угол |
встречи |
волны |
с |
границей |
соответству |
|||||||||
ющего |
слоя (угол скольжения). |
|
|
|
|
|
|
|
||||||
|
Уравнение |
(76) определяет |
форму пути |
электромагнитной |
волны |
|||||||||
в |
атмосфере. |
Дифференцируя |
(74), |
получим |
|
|
|
|||||||
|
|
|
|
|
d ô = |
- t g ô — , |
|
|
|
|
|
|||
а |
из рис. 25 |
|
|
|
|
|
° |
п |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
7 |
dh |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
cos ô |
|
|
|
|
|
|
где dh — толщина |
слоя, |
a |
ds — линейный элемент |
пути |
волны. |
|||||||||
Следовательно, |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
db |
|
1 |
|
dn |
sin ô |
|
|
|
|
|
|
|
|
|
ds |
|
r |
|
dh |
n |
' |
|
|
|
|
откуда |
радиус кривизны |
пути |
волны |
в |
вертикальной |
плоскости |
||||||||
будет |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
I dn |
\ |
|
|
|
|
(77) |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
Ы ) С 0 8 |
а |
|
|
|
|
|
где а — угол |
встречи волны |
с |
уровенной |
поверхностью |
в |
данной |
точке. Вследствие небольшой кривизны пути волны этот угол можно принимать равным углу встречи волны с земной поверхностью в на
чальной |
точке. Величина этого угла при измерении расстояний |
обычно |
близка к 0й . Поэтому можно принять cos ix 1. Показатель |
преломления воздуха тадже близок к единице, поэтому можем на писать следующую приближенную формулу:
'• |
- - / - і у |
( 7 8 > |
|
\dh ) |
|
Таким образом, радиус кривизны пути электромагнитной волны
равен обратной величине градиента |
показателя |
преломления. |
Если бы градиент был постоянным, т. е. показатель |
преломления |
|
изменялся бы с высотой по линейному |
закону, то путь волны пред |
|
ставлял бы собой дугу окружности. |
В реальной атмосфере путь |
волны будет иметь переменную кривизну. С возрастанием высоты модуль градиента уменьшается, а следовательно, радиус кривизны пути волны увеличивается.
Если все члены выражений (71) разделить на dh, то получим формулы для градиента индекса преломления. Для средней атмосферы
57
на уровне моря (р = 1013 мб, Т = 288° К и е = |
10 мб) формулы |
|||
будут |
иметь вид: |
|
|
|
для |
световых волн |
|
|
|
|
^ - 0 , 2 8 ^ - 1 , 0 0 ^ - 0 , 0 4 ^ ; |
|
(79,а) |
|
для |
радиоволн |
|
|
|
|
gN |
= 0,27gp - 1,26g, + 4,48ge. |
|
(79,а)' |
Для |
этой же атмосферы на высоте 2 км (р = 795 мб, Т = 275° К |
|||
и е = 5 мб) получим следующие формулы: |
|
|
||
для |
световых волн |
|
|
|
|
ftv = 0,30g p - 0,86& - 0,04g,; |
|
(79,6) |
|
для |
радиоволн |
|
|
|
|
gN |
= 0,28gp -0,99g, + 4,92ge. |
|
(79,6)' |
Принимая градиенты метеорологических элементов равными (на |
||||
1 м) на уровне моря gp |
= —0,13 мб, gt = —0,006°, |
ge = |
—0,003 мб |
|
и на высоте 2 км gp = |
—0,10 мб, gt = —0,006°, |
ge = |
—0,002 мб, |
|
найдем следующие величины градиента индекса преломления: |
||||
для |
световых волн |
|
|
|
gN = —0,030 на 1 м, — на уровне моря, g w = — 0,025 на 1 м, — на высоте 2 км;
для радиоволн
gN= —0,041 на 1 м, —на уровне моря, g w = —0,032 на 1 м, — на высоте 2 км.
Отсюда видно, что в средней атмосфере радиус кривизны световой волны изменяется от 33 000 км на уровне моря до 40 000 км на вы соте 2 км. Соответствующие радиусы кривизны для радиоволн из меняются от 24 000 до 31 000 км. Полагая радиус Земли R равным 6371 км *, можно приближенно считать, что в стандартных условиях радиус кривизны пути световой волны в среднем равен 6R, а радиус кривизны радиоволны — 4і? **.
Радиус кривизны пути волны в горизонтальной плоскости, пред полагая, что изобарические поверхности совпадают с уровенными поверхностями, выражается формулой, аналогичной (77), если в ней
* За средний радиус Земли везде принят радиус сферы, площадь поверх ности которой одинакова с площадью поверхности земного эллипсоида. Д л я эллипсоида Красовского R = 6371,116 км .
** Величину A = -5- = —R |
называют коэффициентом рефракции |
светового луча или радиолуча. |
|
58
заменить вертикальный градиент показателя преломления гори зонтальным, т. е.
/ |
dn |
\ |
dn |
(80) |
\ |
dD |
) cos а |
ID |
|
Величина этого радиуса в 3000—4000 раз больше соответству ющего радиуса в вертикальной плоскости.
Из изложенного видно, что путь электромагнитной волны имеет
малую |
кривизну. Так, |
для световых |
|
волн при длине линии 60 км |
|||||||||
наибольшее отклонение |
траектории от прямой (отклонение по высоте) |
||||||||||||
составляет |
около |
12 |
м, |
а |
разность |
длин |
|||||||
кривой |
и |
стягивающей |
ее хорды |
меньше |
|||||||||
7 мм. |
Для |
радиоволн |
при |
|
расстоянии |
||||||||
200 км |
соответствующие |
величины |
будут |
||||||||||
180 м |
и 44 |
см. |
Поэтому |
при |
обработке |
||||||||
радиоэлектронных |
измерений |
допустимо |
|||||||||||
считать |
путь распространения |
волны |
ду |
||||||||||
гой окружности |
|
некоторого |
среднего |
ра |
|||||||||
диуса. Величину этого радиуса или при |
|||||||||||||
нимают |
равной |
указанным |
выше |
вели |
|||||||||
чинам 6І?, или |
iR, |
или |
же |
находят |
по |
||||||||
упрощенным формулам. Так, если известны |
|||||||||||||
показатели |
преломления |
пх |
|
и |
п 2 |
двух |
|||||||
точек с высотами |
Н1 |
и |
Н2, |
то на |
основа |
||||||||
нии (77) можно |
написать |
|
|
|
|
|
|
||||||
|
|
г |
= |
|
Н2—Hi |
|
|
|
|
|
|
(81) |
|
|
|
(ni — пг) cos |
a |
|
|
|
|
ѵгде a — угол встречи радиоволны с земной поверхностью в точке стояния антенны, величину которого можно найти по приближенной формуле
sin a = • |
D |
|
В случае, если известны показатели преломления для нескольких точек пути волны, весь путь можно представить как совокупность сопряженных дуг окружностей различных радиусов. Необходимые высоты точек траектории можно н айти на основании следующих рассуждений.
Пусть точка / (рис. 26), расположенная на пути электромагнит ной волны, высота которой Н( подлежит определению, удалена от точки А (начальной) на расстояние st. По рис. 26 напишем
|
Bf |
• кг |
-J- h2 |
где Н1 |
и # 2 — высоты точек А и В; |
s — расстояние между их проек |
|
циями |
на сфере; h1 и h2 —стрелки |
прогиба |
соответствующих дуг. |
59