Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

shimbirev_b_p_teoriya_figury_zemli

.pdf
Скачиваний:
102
Добавлен:
04.06.2015
Размер:
14.99 Mб
Скачать

 

 

 

Значение радиусов концентрических зон

 

Т а б л и ц а ®

 

 

 

 

 

№ зоны

Радиусы зон, км

 

 

№ зоны

Радиусы зон, км

№ зоны

Радиусы зон. ш

 

5,0

 

 

 

 

XI

198,6

 

 

X X I

 

1000,0

I

 

 

 

 

XII

 

 

 

X X I I

 

7,3

 

 

 

 

 

246,7

 

 

 

1163,7

II

 

 

 

 

XIII

 

 

 

XXIII

 

10,7

 

 

 

 

 

305,4

 

 

 

1345,5

III

 

 

 

 

XIV

 

 

 

X X I V

 

15,7

 

 

 

 

 

357,8

 

 

 

1545,6

IV

 

 

 

 

XV

 

 

 

X X V

 

22,8

 

 

 

 

 

418,1

 

 

 

1763,9

V

 

 

 

 

XVI

 

 

 

X X V I

 

33,3

 

 

 

 

 

487,4

 

 

 

2000,0

VI

 

 

 

 

XVII

 

 

 

 

 

48,5

 

 

 

 

 

566,6

 

 

 

 

 

VII

 

 

 

 

XVIII

 

 

 

 

 

 

70,6

 

 

 

 

 

656,6

 

 

 

 

 

VIII

 

 

 

 

X I X

 

 

 

 

 

 

102,6

 

 

 

 

 

758,0

 

 

 

 

 

IX

 

 

 

 

X X

 

 

 

 

 

 

128,0

 

 

 

 

 

872,0

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

159,6

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ак =

— 0,005" соз Ак,

 

Ък =

— 0,002" соз

 

 

 

с* = — 0,0015" соз Ак,

 

 

 

 

 

 

ак=

-0,000871" соз Ак,

 

 

 

 

 

 

VIII

16

 

XIII

 

24

 

 

 

 

Л5-2000 = 2

2

к=1

— У ) 1 к + 2

 

2

ъ '

к ( § — у)1 к

+

 

 

 

 

1=1

 

 

/=1Х Ь=1

 

 

 

 

 

 

XXI

24

 

 

 

XXVI

 

48

 

 

 

 

 

+

 

2

2

 

 

2

<=ХХН

2

 

 

 

(х.:-:-.

где

 

/ - X I V й=1

 

к=1

 

 

 

ак == — 0,005" з т

 

=

 

— 0,002" з т Ак<

 

 

 

 

 

 

с к = -

0,0015" з т Ак,

й'к =

-

 

0,000871" з т

Ак.

 

Для вычисления по этим формулам используют палетки, предложенные

В. Ф. Еремеевым: палетка

(рис. 61)

для учета

аномалий в

первой областж.

палетка (рис. 62) для учета аномалий во второй и в третьей областях и палетк± для учета аномалий в четвертой области.

Палетки изготавливают на прозрачной основе в масштабе гравиметрической карты.

Как указывалось выше, для целей астрономо-гравиметрического нивели-

рования достаточно принять сферический радиус гр0

2—3°, что соответствует

учету XIII зоны палеток Еремеева (305,4 км).

Г

В настоящее время появилась возможность вычислять влияние зон VII — XIII, т. е. от 48,5 до 305,4 км, на ЭВМ по аномалиям в свободном воздухе, осредненным по площадям размером 10 X 15'.

300'

Влияние зоны от 0 до 48,5 км на уклонения отвеса пока что вычисляют вручную, путем численного интегрирования с использованием палеток, для чего палетку накладывают на гравиметрическую карту таким образом, чтобы центр палетки совпал с определяемым пунктом, а линия N5 — с меридианом. Затем в каждой трапеции определяют среднее значение аномалии у), умножают на соответствующие коэффициенты и подсчитывают суммы полученных произведений.

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

Рис. 61

Рис. 62

спомощью гипсометрических карт. Переход от аномалий Буге к аномалиям

всвободном воздухе осуществляется по формуле (VII.25)

где 2я/б — коэффициент Буге, равный 0,1119, при плотности б, принятой равной 2,67 г/см3 и 0,0964 при б = 2,3 г/см3.

Таким образом, в слагаемых формул (Х.25) и (X.26) появляются два члена. Их вычисляют с помощью палеток, первый член — по картам аномалий Буге, второй — по гипсометрическим картам, умножая результаты вычислений на коэффициент 2я/б.

При учете влияния дальних зон с помощью ЭВМ используются осредненные по трапециям аномалии Фая. В высокогорных районах учитывают также

нелинейную часть влияния на уклонения отвеса топографического массива — поправки А|р и Дг|р (см. формулы VIII.119 и VIII.120). Для учета влияния

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

301'

масштабов 1 : 100 ООО, 1 : 200 ООО, 1 : 300 ООО и 1 : 1 ООО ООО Все палеткп построены по одному принципу. Внутри 5-километровой зоны выделяются кольцевые зоны такого же влияния на уклонения отвеса, как зоны I—VI: они названы латинскими буквами А, В, С, Б, Е, Р. Радиусы этих зон даны

в табл. 20.

Счет зон ведется к центру палетки. Число зон различно в зависимости от масштаба карты.

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

вольно выбранного радиуса

р0

и

в

8 точках

на

окружности

радиуса ?-ср =

 

 

 

 

 

=

0

+

г0)/2.

 

 

 

 

 

 

 

 

Т а б л и ц а

20

 

 

При выводе формулы для вычисления влия-

 

 

 

 

 

ния на уклонения отвеса

зоны радиуса

г0 в ин-

 

Зона

Радиусы зон, км

тервале

от 0

до

р0

применен метод градиентов.

 

При этом полагают,

что аномалия

силы тяжестп

 

 

 

 

 

 

 

 

 

 

внутри окружности радиуса р0 может быть пред-

 

А

 

5000

 

ставлена

соотношением

 

 

 

 

 

3,419

 

 

 

 

 

 

 

 

 

 

 

 

 

В

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2,338

 

 

 

 

От - V) = От -

Т)о +

 

г,

(Х.27)

 

С

 

 

 

 

 

 

 

 

1,599

 

 

 

 

 

П

 

1,094

 

 

 

 

 

 

 

 

 

 

 

 

 

Е

 

0,748

 

где

(§ — у) о — значение

аномалии в

исследу-

 

Р

 

 

 

 

0,511

 

емой

точке;

д (§ — у)/дг — значение

горизон-

 

 

 

 

 

 

 

 

 

тального градиента аномалии вдоль данного

 

 

 

 

 

радиуса.

 

 

 

 

 

 

 

Влияние нулевой зоны на составляющую уклонения отвеса в меридиане

вычисляется по

формуле

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ро 2Я

 

 

 

 

 

 

 

 

 

Го-Ро=

 

 

 

о о

( ё - У ) т С 0 8 А й А й г >

 

<х-28>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

р" =

206 265".

 

 

 

 

 

 

 

 

 

 

 

 

После подстановки в (Х.28)

значения аномалии (Х.27) получим

 

 

 

 

 

 

 

 

 

 

 

Ро

 

 

 

 

 

 

 

 

 

 

 

 

О

 

О

 

 

 

 

 

 

2Л

 

Ро

 

 

 

 

 

 

 

 

 

 

 

 

о

 

о

 

 

 

 

 

 

 

 

 

 

где

=

0.

 

 

 

 

 

 

 

 

 

 

 

 

 

При вычислении / 2

будем считать градиент аномалии для данного азимута

величиной постоянной

вдоль

всего радиуса р0

 

 

 

 

 

 

 

 

 

д (8-—У) _

(8—у)—(е—- у)о

'

 

 

 

 

 

 

 

 

дг

 

 

 

ро

 

 

 

 

1

Рассчитанные

Е. М. Орловой.

 

 

 

 

 

 

 

 

 

302'

где под — у) будем понимать значение аномалии на расстоянии р0 от исследуемого пункта. Тогда

2 Я Р0

о

о

о

 

— |

— ?)о с о 8 А ( 1 А= ^ (8 —

у)созА<1А,

о

о

 

так как

 

2 Я

 

^ {8 — У)ос03 Ав,А = {§~- у)0 ^

соз АйА = 0.

о

о

 

Таким образом,

2 Я

о

При интегрировании по А воспользуемся формулой трапеций

ь

 

 

| / (г) йх =

[1 / К) + /(X,) + /(X,) + . . .

+/(**-!)+

а

положив: пределы интегрирования а = О, Ь = 2я, число ординат к = 8, подынтегральную функцию / (ж) = — у) соз А.

Тогда получим

5о-Ро = - - 2 ^ ~

0 0 8 А ь

ь

и аналогичную формулу для составляющей уклонения отвеса в первом вертикале

 

о

Ло - Ро = — 2

( 2 — У ) к з 1 п А к -

 

к=1

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

Приведем окончательную формулу для вычисления влияния нулевой

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

 

к=1

Л=1

где А§к 0) — аномалия в свободном воздухе в точке с номером к на окруж-

ности радиуса р0;

А%к ср) — аномалия в аналогичной точке на окружности

радиуса гср

_ _

Ро+^о .

 

 

г°р ~~

2

Ак = пк/4, к =

2, 3, . . ., 8.

 

303'

Аномалии в свободном воздухе в формуле (Х.29) получают методом косвенной интерполяции, проводя численное интегрирование по картам аномалий Буге и гипсометрическим картам.

Рассмотрим случай, когда в центральной зоне (радиуса 5 км) палетки выделены 4 кольцевые зоны (А, В, С и ТУ) такого же влияния на уклонения отвеса, как зоны I —VI.

В этом случае радиус нулевой зоны г0

= 1,094 км, р0 = 0,300

км, гср

=

= 0,697

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

по формуле

 

 

 

8

8

 

 

 

{ 5 и

= -°'0 2 6 3 2 " 2 ^ Ц I Аи } -

°'0 2 "8 " 2

^ Сер) {

Ак } •

(Х'3

 

Й=1

к-Л

 

 

 

Если центральную зону разбить на иное число кольцевых зон (3, 2 или 1), то коэффициенты формулы (Х.ЗО) можно сохранить при соответствующем подборе значения р0. Радиус нулевой зоны г0 всегда равен внутреннему радиусу последней кольцевой зоны. Так, если центральная зона разбита только на три

зоны (4* В и С), то радиус нулевой зоны г0 =

1,599 км„ если на две зоны

и В), то г0 2,338 км и если — на одну (А),

то г0 = 3,419 км.

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

нения отвеса.

 

 

Формулу (IX.60) целесообразно

представить в этом случае в виде

Б а г - Б а г

=

+

где

 

 

= -0,00449 [(А?л + Д|в) АВ +

(Лги + Л*1в) сов Вт А/,],

Рабочие формулы для вычисления влияния на Д^ (2) зон I—VI и А — Р имеют вид

 

 

 

 

 

16

 

16

 

 

 

 

Л^ (2)1-^1 = 0,00014 2

Аек (I) +

0,000212 Д& (И) +

 

 

 

 

 

 

к=1

 

к=1

 

 

 

 

1 6

 

 

 

16

1 6

 

 

 

+

0,00030 2

Аёк (III) +

0,00044 2 А§к (IV) + 0,00062 2

Аек (V) +

 

 

к=Х

 

 

 

к=1

к=1

 

 

 

 

 

 

+

0,00085 2

(VI),

 

 

 

 

 

 

 

 

к=1

 

 

 

 

 

 

16

 

 

16

 

16

 

 

ДС(2)я_, = 0,000102 А ^ И ) + 0,000072

+ 0,00005 2

Аёк(С) +

 

 

 

к"1

 

 

к=1

 

к-1

 

 

 

+ 0,00003 2

 

+ 0 , 0 0 0 0 2 2 Ьёк{Е) + 0,000022

А8к(П

 

 

к=1

 

 

й-1

к-1 .

 

где

Д^б — аномалия

в

свободном

воздухе

в соответствующей

зоне, к = 1,

2,

3, . .

16.

 

 

 

 

 

 

 

304'

Влияние на А^ (2) центральной круговой области, где численное интегрирование проводится по точкам, зависит от ее радиуса г0, а также от выбора радиуса р0

 

 

 

дг; (2)0_Го = -На. ^ Дг* (Ро) +

2

М -

 

 

 

 

к-1

к=1

 

где

0) — аномалия в свободном воздухе

в точке с номером к, располо-

женной на окружности радиуса р0; А§ь (гср) т 0 же> н а

окружности радиуса

г с Р = ( Р о +

г о ) / 2 >

 

 

 

к =

1 ,

2 ,

3 , . . . ,

8 .

 

 

Для палетки с зонами А — Р имеем

 

 

Гд =

0,511

км, ро =

0,140 км, гср = 0,326

км,

 

 

 

 

 

8

8

 

 

 

А^(2)0_Го = 0,000022 А^(Р0) + 0 , 0 0 0 0 4 2

Д ^ М -

 

 

 

 

к=1

к=1

 

§64. СПЕЦИАЛЬНЫЕ ПАЛЕТКИ ДЛЯ ВЫЧИСЛЕНИЯ ГРАВИМЕТРИЧЕСКОЙ ПОПРАВКИ В РЕЗУЛЬТАТ

АСТРОНОМИЧЕСКОГО НИВЕЛИРОВАНИЯ

Для определения гравиметрических поправок в результаты астрономического нивелирования применяют специальные палетки, предложенные различными авторами: М. С. Молоденским (1937 г.), А. К. Маловичко (1952 г.), Фан-Цзюном (1958 г.), И. М. Тироном (1963 г.) и др. Все палетки, построенные для вычисления гравиметрических поправок, можно разделить на два вида. К первому виду относятся палетки, область интегрирования на которых представлена рядом точек с вычисленными для них коэффициентами влияния аномалий силы тяжести на превышение квазигеоида между конечными пунктами А и В нивелирования. Ко второму виду относятся палетки, область интегрирования на которых разбита координатными линиями на элементарные площадки. Аномалии силы тяжести в элементарной площадке считаются постоянными во всех ее точках. Каждой площадке соответствует свой коэффициент влияния аномалий силы тяжести на превышение квазигеоида между астрономо-геодези- ческими пунктами А и В.

К первому виду относится палетка, построенная Молоденским в системе биполярных криволинейных координат с фокусами в точках А и В.

 

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

сов

(с полуосями: а0 = 1,1547, с0 = 0,5744, а = 3,00, с = 2,83, а = 4,00,

с =

3,87; за единицу принято I — половина расстояния между астрономо-гео-

дезическими пунктами А и В) разбивается на 10 участков, обозначенных на рис.63 через 8г, 8%, 8й, 54, 8Ъ 86, 87, 8а, 89, 810.

Для выделения околофокальных участков и (на рис. 63 они заштри-

хованы) вокруг пунктов А ж В, как фокусов, проводятся софокусные эллипсы (с полуосями а0 ж с0) ж гипербола (с полуосями Ь0 ж й0).

Исходя из условия, чтобы внутри фокальных участков аномалии силы тяжести изменялись линейно, Молоденский остановился на следующей сис-

теме значений (принимая I = 1):

 

 

00 = 1.1547,

с0 =

0,5744,

&0 = 0,9041,

й0 =

0,4274.

2 0 заказ 1379

 

305

Каждый околофокальный участок оказался заключенным внутри пря1*

угольника по оси х со сторонами 0,2506 (а0 Ъ0 = 0,2506)

и 0,4936 по

у (2у = 2с0й0 = 0,4936). Таким образом, при расстоянии АВ,

равном 100 :-

сторона прямоугольника по оси х достигает 12,5 км, а по оси у — 25 км.

Влияние околофокальных участков вычисляется особо, а влияние остальных участков — по формулам численного интегрирования Гаусса. Со< ветственно этому поправка разбивается на три поправки: Д^д, Д 1,в и Д^

Первые две

и Добусловлены влиянием околофокальных участков, огр

ниченных координатными линиями а0

и Ь0, третья Д^с — влиянием остальн:.

части области 2 . Таким образом,

 

 

Формула для учета влияния околофокальных участков

имеет вид

Л2Дш1'В =

+0-02755" [{& - у)й0 -

- у)Ьо+ {(* - V)».

У)«.}в 1

Точки а0 и Ь0 около пунктов А и В, в которых должны быть определен"-

значения аномалий — у), показаны на рис. 63.

Для вычисления поправки Д^с Молоденским получена формула

I /

Коэффициенты Ац в этом выражении вычисляют по формуле численно: интегрирования Гаусса.

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

имеет вид

 

 

 

 

 

= -

у

+ <

г )

+ 0,02755" [{(* - у)а,

- ( в - 7 ) * о } А +

+

{(г -

у)ь. - {

§ -

?)а.}в 1 + 2 Ац (в -

у)и,

 

 

 

 

а

 

306'

где у)а0 и у)ь0 — аномалии силы

тяжести, определяемые в точках

в0 и 60 по гравиметрической карте; у)ц

— аномалии силы тяжести в точ-

ках г]; Ац—коэффициенты влияния аномалий на превышение квазигеоида, относящиеся к точкам г/ и выраженные в 0,00001 долях секунды.

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

А ж В.

Эти палетки вычерчиваются на прозрачной основе, на них наносятся фокусы координатных линий — пункты А и В ж те точки г], для которых с гравиметрической карты должны быть сняты аномалии. Рядом с этими точками подписывают соответствующие коэффициенты А Э т а палетка изображена на рис. 64. Палетку накладывают на карту таким образом, чтобы точки А и В палетки совпали с астропунктами на гравиметрической карте. Далее, в отмеченных точках палетки определяют с карты значения аномалий — у),/ и умножают на соответствующие коэффициенты Ац.

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

17*

307

метрической карте. Высоты Н в точках палетки при этом умножаются 0,0963 и на коэффициенты Ац палетки. Сумма результатов вычислений по к. аномалий Буге и гипсометрической карте определяет гравиметрическую " правку.

Аналогично выполняется интерполирование аномалий в неполной топе фической редукции.

Ко второму виду палеток относятся палетки, предложенные Малови*1*" Фан-Цзюном и Тироном. Они построены в системе прямоугольных коорднв:

Например, А. К. Маловичко была рассчитана квадратная палетка. - область интегрирования разбита на квадраты со сторонами 0,41, 0,21, 0,051. Размеры квадратов увеличиваются с удалением их от астропунк* Для каждого квадрата палетки коэффициенты вычисляются по формуле

А„ —

4яу 8Щ 1"

х 1П ? + / ( * + ! ) 2

+ г / 2

+ у 1 п

.*+ 1 + 1

)2+г/2

«1-1

 

 

 

 

 

 

 

 

где хь

уI — прямоугольные координаты

углов

квадратов

палетки, выра—

ные в единицах половины расстояния между астропунктами I. Влияние г " ральных зон, расположенных вокруг астропунктов, учитывается с поые-ш

дополнительной

палетки.

 

 

 

Гравиметрические поправки с использованием квадратной палетки

деляются по формуле

 

 

 

2 г Д 8 ^ г

= {0,02949" [(Д& -

А ё а ) в + (А^а - Д ^ Ы + 2

Аё п Ап ,

где

и А^ь — значения аномалий

силы тяжести для точек

а и Ъ пале?™

Глава XI

ДВИЖЕНИЕ ИСКУССТВЕННЫХ СПУТНИКОВ ЗЕМЛИ В ПОЛЕ ТЯГОТЕНИЯ И ОСНОВНЫЕ УРАВНЕНИЯ НЕБЕСНОЙ МЕХАНИКИ

§ 65. ДВИЖЕНИЕ ТЕЛА В ЦЕНТРАЛЬНОМ ПОЛЕ. ВОЗМУЩАЮЩАЯ ФУНКЦИЯ

Законы движения небесных тел, открытые Кеплером в начале XVII века на основании обработки многочисленных наблюдений, являются по сути следствием закона всемирного тяготения Ньютона. Эти законы не получили в то время теоретического объяснения и их динамический смысл был выявлен Ньютоном после того, как он сформулировал свои законы движения. Ниже будет показано, что если тело совершает движение под действием силы ньютонианского тяготения, то орбита этого тела должна лежать в плоскости, занимающей неизменное направление в пространстве, и являться одним из конических сечений, в частности эллипсом. Он получил наименование Кеплерова эллипса.

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

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

случае / х

= 0

и

на основании (IV.34)

 

 

 

р

00

00

П

 

V = ^

1 -

2

(т )" 7»Р» <со8 0) + 2

2 (т)" ^ 008 кК + 81п

рпь (0) '

 

[_

п=2

п-2 к=1

^

 

 

 

 

 

(Х1.1)

где г — радиус-вектор орбиты спутника.

Главный член этого ряда Щ - соответствует случаю центрального поля,

а вся совокупность последующих членов представляет собой отличие потенциальной функции Земли от потенциальной функции шара.

309.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]