Определение углов полной рефракции при больших зенитных расстояниях Текст научной статьи по специальности «Физика»

Определение углов полной рефракции при больших зенитных расстояниях

В.И. Куштин, Д.Л. Дробязко Ростовский государственный строительный университет, г. Ростов-на-Дону

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

,^п

&с = -р —Щ р (1)

п

где ф - угол между радиусом сферы и направлением оптического луча в текущей точке; п

- коэффициент преломления газовой среды или воздуха, р” = 206265” [1].

Недостатком этого уравнение является то, что при ф = 900 dф = да, хотя известно, что в горизонте угол астрономической рефракции равен не бесконечности, а примерно 35’. Поэтому выражение (1) используют обычно при углах ф, отличающихся от 900. Кроме того, при статическом состоянии атмосферы изодиоптрические поверхности обычно совпадают с уровенными, положение которых в пространстве более точно аппроксимируется поверхностями, параллельными поверхности Земного эллипсоида. Поэтому при определении углов рефракции необходимо учитывать и это обстоятельство.

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

В инварианте Снеллиуса

п sinр = п' Бтр' (2)

п ± п' р±р , * , *

примем —-— = пс , —-— = рс , п - п = Ап , р - р = Ар .

Тогда

Ап , Ап Ар , Ар

п = пс ±— , п = пс-------------, р = рс----------, р =рс ±---------.

с 2 с 2 2 2

С учетом этих значений инвариант Снеллиуса примет вид

(пс ±Т )

После преобразований имеем

Ап | . ( Ар | ( Ап | . ( Ар

пс +-2 ) Ч рс ^ 1 = 1 пс -) Ч рс + —

Ар Ап

= ~ Шс . (3)

2 2пс

Полученная формула является строгой и ее можно использовать при любых ал

во внимание п = 1 + N -10-6, где N - индекс преломления,

Ар Ар

значениях ф. Учитывая малое значение Аф, можно принять ~ . Тогда, принимая

Ар = р"-10-6 щ«^«+1

па 2

Угол полной рефракции

= ^Ар = Р-10-6 р1-р^^ (4)

1 п • 2

Г

где N = (n - l)l06 - индекс преломления.

Для определения входящих в формулу (4) углов Ф , фi+1 в сферической модели атмосферы используют известный инвариант an sin р = c = a1n1 sin р1 , откуда

an sin ^i

sin р = —1—1--L, (5)

an

где а 1 - расстояние от центра Земли до начальной точки траектории оптических волн, а ¡

- до текущей точки i; n 1, n - коэффициенты преломления в начальной и текущей точках

соответственно; ф 1 , ф - углы между нормалью к изодиоптрической поверхности и

траекторией луча в начальной и текущей точках.

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

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

N

1 + e’2 cos2 B cos2 A ’

Ra = i , „.г*»* (6)

где радиус кривизны первого вертикала

N = , \ = ■ С . (7)

V1 - e2 sin2 B V1 + e'2 cosB

В формулах (6), (7) Є и е’2 - квадраты первого и второго эксцентриситетов эллипса; а, с - радиус кривизны экватора и полярный радиус кривизны; В - геодезическая широта; А - азимут нормального сечения.

Значения Ra , вычисленные по формуле (6) для различных значений широты В и азимута А, показали, что максимальная разность Ra достигает 64,146 км в плоскости меридиана. Выясним, на каком расстоянии от пункта наблюдения оптический луч будет выходить из атмосферы. Известно, что на высоте 30 км в стандартной атмосфере ГОСТ 4401-81 давление р = 1,93 гПа, Т= 226,509 К. На уровне моря р 0 = 1013,25 гПа, Т0 =

NT Р

288,15 К, N0 = 278,24. С учетом этих значений N30 = —■ — = 4,17 . На высоте 40 кмр =

Ро T

0,287 гПа, Т = 250,350 К, а значение N40 = 0,091.

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

атмосферы с радиусом кривизны Ra на высоте 30 км. Для слоев одинакового

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

(Rg + H h sin <pg =(r + H1h sin (p1, для отрезка 1 - 2 (R1 + H1 )n1 sin (p1 = (R1 + Hг )пг sin (рг, для отрезка 2 - А (R2 + Hг )пг sin (рг = (R2 + Hn )na sin <pa.

^ (R2 + Hn )na sin (p

Подставляя из последнего выражения пг sin (рг =-------------------- во второе, получим

R2 + H г

( \ • (R + н2 )(R2 + н)

(R1 + н 1 )ni sin (px =----- ---- -------na Sin (pa .

R2 + H 2

Подставляя из полученного выражения n1 sin (р1 в первое, находим

к + Н

Я + Н Хл, + н 2 Хл + Н) .

8Ш <Ре =----------------------------77.-77777-^-па в1П сра .

(Я, + Н, )Я + Н 2 )

Или (я^ + Н ^ БІП

1 -

Я1- Я V

V Яі + н іу

1 - Я2 - Я1

Л (

Я2 + Н 2 У

1 -

Я - Я

Я, + я

а У

(Яа + На )па ^П

После небольших преобразований с учетом величин первого порядка малости имеем

(Я + Н П БІПФ =

1-

Яі - Яг

Я2 - Я1

Я - Я.

Я1 + Н1 Я2 + Н 2

Я + Н

(Яа + На )Па ^П

а У

Учитывая бесконечно малые отрезки трассы электромагнитных волн О - 1, 1 - 2, ...,п -А, получим

+ Нё)пё 81ПФ =

а

і-І

V Я

dЯ

Я + Н

(Яа + На )па *ІП Фа .

Используя при интегрировании теорему о среднем значении, находим

(

Я + Н*г-г

БІПФ

Ґ Я„

1 -V

Я

(Яа + На )Па §ІП

Фа .

(8)

(9)

(Я ± Н )т )

Значения N для изотермической атмосферы, в которой температура с высотой изменяется по линейному закону, вычисляют по формуле

N = N

( т ^

т

Vі0 У

Яв

Сравнение значений гс, вычисленных по формуле (4) для фё = 890 и фё = 900 и условий стандартной атмосферы ГОСТ 4401-81 численным интегрированием с шагом по высоте 0,25 км и путем удвоения шага (0,50 км) показало, что разности значений гс не превышают для фё = 890 0,15”, а для фё = 900 при Н через 0,1 км и для Н через 0,2 км, практически дают одинаковые результаты (разности не превышают 0,07”), что свидетельствует о высокой точности полученной формулы и ее возможности определять углы гс при любых значениях фё, включая и фё = 900. т.е. точность определения углов гс является достаточно высокой.

Литература

1. В.И. Куштин. Учет влияния атмосферы на результаты измерения длин радиоэлектронными системами. Монография. Москва, МИИГАиК, 2003, 179с.

2. Атмосфера стандартная. Параметры ГОСТ 4401-81.-М.: Издательство стандартов, 1981, 280с.