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

shimbirev_b_p_teoriya_figury_zemli

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

Координаты В ж Ь входят в неявном виде в функцию Стокса 8 (тр). Это

вытекает из формулы (111.17), заменив в которой 0 через 90° — В, получим

соз ар = з т В зхп В' + соз В соз В' соз (Ь' Ь).

(VIII.54

Следовательно, можно написать

 

 

 

 

1

 

 

 

 

 

4яу О Л уь

 

Щ

 

Г\ = .

1

ГГ/..

л ^

Зф

 

4лу

Цт-у)

йг|)

сое В дЬ (1(0.

 

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

81П яр з т А = соз В' з т {Ь" Е)

з т 1р соз А з т В' соз В — соз В' з т В соз (Ь' — Ь)\'

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

данной точки (В, Ь).

Дифференцируя (VIII.54), по В и Ь, получим

— з т гр = з т В" соз В — соз В' з т В соз (Ь" — Ь),

(VIII.55)

(В', Ь'), из

- з т а р ^ = СОЗ В' СОЗ В 81П (Ь' — 1).

Принимая во внимание формулы (VI 11.55), будем иметь

^

л

 

ТТГ =

соз А

 

дБ

 

(УШ.5б!

 

 

соз В дЬ

: 81И А

 

 

Используя эти соотношения, представим

| и г) в виде

 

дурр)

соз А (1(0

4117 И и г - т )

а\р

 

 

(VIII.57)

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

1

оо

Я2 Я

-4.туч81П- 1" и

* - ^ ®

8^плр^ з т А й\р <2А.

 

о о

 

 

200'

Введем обозначение

() = ~

1

Л8 (ф)

81П

 

(УШ.58)

8111 1"

|)

 

:-та функция называется функцией Венинг-Мейнеса.

 

 

Окончательно получим

 

 

 

 

 

о о

 

 

 

 

(УШ.59)

я 2 я

 

 

 

 

Дифференцируя функцию Стокса (VIII.38), найдем

 

 

г|)

 

 

 

 

 

рпе '

 

 

 

 

 

- 5 81П 1|3 — 3 С 0 8 у +

3 8тг|)1п (з1П

31П2

1 СОЗ 1|3 ^

СОВ

+ 81П ^ СОЗ - у )

 

 

Проведем алгебраические преобразования

й8 (Ч>) .

^ -

- созес - у — 6 81П

+

20 81П2 - у +

ау

81П1|5 = С 0 8 2

 

 

 

 

 

 

 

 

3 соз

( 1 +

2

+ 1 2 з ш 2 | - 1 П ( з т | . + 8 т 2 ^ ) -

1 + 81П

 

 

 

 

:реобразуем последний член, стоящий в скобке,

Зсозг|> ( 1 + 28Ш

3

б 8 ш - 1 ^ 1 _ 8 т - | - — 2 з ш 2 - | - )

\

1 + 8111-

1 + ЗШ-|-

1

 

 

 

 

в й п | - ( 1 - 2 й п 4

. ) =

+

6 81П

12 81П2

1 + 81П

 

1 + 8т^|-

 

 

 

 

окончательно

 

 

 

 

 

 

81П

— -СОЗ2 ± [созес - 1 + 12 8Ш ±

-

32 81П2

+

 

+ — — - — ^ —128ш2 -у 1п (зш| - + 8 т 2 4 - ) ] .

 

(У1И.60)

201

Таким образом, функция Венинг-Мейнеса

 

V м = -^ТПГГ с о з 2 1 [ с о з е с т

+ з т

- 32 эт* ± +

+

Ц : г - 1 2 8 ш ^ 1 П

 

(УШ.61)

 

1+8111-5-

4

' л

 

и

 

 

Числовые значения функции (} (гр) приведены в табл. И.

 

 

 

7 = 979,77

гал*

 

Т аб л и ц а

11

 

 

 

 

 

 

 

г|)°

 

Я

 

Я

 

 

Я

 

0

 

с\о

60

+ 0 , 2 2

130

 

- 0 , 3 4

 

1

 

+ 1 2 , 3 5

70

+ 0 , 0 3

140

 

- 0 , 2 4

 

10

 

+ 1 , 5 9

80

- 0 , 1 5

150

 

— 0 , 1 6

 

20

 

+ 1 , 0 2

90

- 0 , 2 9

160

 

— 0 , 0 8

 

30

 

+ 0 , 7 9

100

— 0 , 3 8

170

 

— 0 , 0 2

 

4 0

 

+ 0 , 6 1

110

— 0 , 4 1

180

 

0,00

 

50

 

+ 0 , 4 3

120

- 0 , 4 0

 

 

 

 

* Для

V взято

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

Гельмерта

1901 — 1909 гг.

 

 

 

 

 

 

Функция

(гр) непрерывна во всей области 0°

гр ^ 180°,

кроме точки

гр = 0°. С увеличением

значения функции () (гр) быстро убывают и при гр

^

70° (? (гр) =

0. Затем

функция (} (гр)

принимает отрицательные

значения,

убывая до —0,41 при гр=110° и далее возрастает до 0 при гр = 180°.

Эти осо-

бенности изменения функции () (гр) учитываются при вычислении

| и т] мето-

дами численного интегрирования.

 

 

 

 

 

Формулы Стокса (VIII.53) и Венинг-Мейнеса (VIII.59) следует рассматривать как первое приближение к определению высоты квазигеоида и составляющих уклонения отвеса в точках физической поверхности Земли.

§46. ОПРЕДЕЛЕНИЕ ВОЗМУЩАЮЩЕГО ПОТЕНЦИАЛА ЗЕМЛИ НА ЗЕМНОЙ ПОВЕРХНОСТИ (РЕШЕНИЕ МОЛОДЕНСКОГО)

Краевое условие для возмущающего потенциала, полученное М. С. Молоденским для поверхности 2 (VIII.8) с учетом (VIII.6), имеет вид

где индексы 2 и 8 означают, что соответствующие величины относятся либо

к поверхности Земли первого приближения 2, либо к физической поверхности 8.

В качестве исходной уровенной поверхности нормального потенциала берется сфера радиуса В и считается, что радиус-вектор точки на поверхности 2

равен

Рх = В+Н1.

202'

Проблема состоит в определении функции Т, которая на поверхности 2 удовлетворяет условию (VIII.62), вне Б является гармонической функцией координат и на бесконечности регулярна. М. С. Молоденский предложил рассматривать возмущающий потенциал во внешнем пространстве как потенциал простого слоя плотности ф, распределенного на поверхности 2,

- И * <22. (УШ.бЗ)

В таком случае задача сводится к определению функции ф. Внешнюю1 производную потенциала (VIII.63) на слое 2 по направлению Н получим ана-

логично (VIII.41)

'(ЭТ.)

 

 

V дН

)п->1Г<

 

 

 

 

 

 

 

и Ф

СОв(г,

НУ)

IV

 

0_

л

 

 

 

 

 

,

 

 

— 2яф соз а,

 

 

 

 

 

 

 

 

 

 

(VIII.64)

 

 

 

где а — угол

между направле-

 

 

 

нием Щ и нормалью п к поверх-

 

 

 

ности 2

 

(рис.

51).

 

 

 

 

 

 

 

 

Обозначая радиус-вектор дан-

 

 

 

ной на поверхности 2 точки N,

 

 

 

определяемой координатами Б,

Ь

 

 

 

н Н1 через рп, а радиус-вектор

 

 

 

текущей

 

точки

К,

определяемой

 

 

 

координатами

 

В',

Ь',

Я7 '

че-

 

 

 

рез р',

получим

 

 

 

 

 

 

 

 

Отсюда

 

 

 

 

р'2 = р02 + гг-2р0лсоз(;\

р0).

 

 

 

 

 

 

 

 

•Ро

 

 

 

 

соз (г,

р0) = соз(/\

Щ) =

 

 

 

 

2рог

2ро

 

 

 

 

 

 

 

 

 

 

 

Подставив значение соз (г, Ну) в (VIII.64), получим

 

 

 

 

 

 

1

я2

ф-

-Р8

Й2 — у - Л ^ - й 2 - 2 л ф с о з а . (УШ.65)

 

 

 

 

 

 

 

 

\дН

)Н^-Н'1

 

2ро

гЗ

 

2

 

Подставим (VIII.65) и (VIII.63) в краевое условие (VIII.62)

_лф соз а = д — у •

.

з ("Сш

^ се п'2-п?

Ро

2ро

 

 

 

это и есть основное интегральное уравнение Молоденского.

Существует несколько способов решения этого уравнения. Изложим способ, предложенный М. С. Молоденским, с учетом уточнений, сделанных впоследствии В. В. Броваром. Приведем (VIII.66) к более простому виду.

203'

Учитывая, что р' = В + , р0 = Д + Щ, а р' - р0 = Нг - Щ. выражение, стоящее в последнем члене уравнения (VIII.66), можно представить как

р'2 -р§

_ (яу '-Я0 у )(р'+Ро)

2 р о ( Я * ' - Я г ) + р ' ( Я Г - Я У ) - р о ( Я * ' - Я У )

 

о/"3

 

0гЗ

 

 

0гЗ

 

 

_ Я ^ - Я У

р ' - р о )

 

Я ^ - Я ? ,

 

^ Т Т Т Й -

 

* г»

""

2^75

Г

"

^

'

Используя обозначение

Я

перепишем уравнение (VIII.66) в более простом виде

( У Ш . 6 8 ,

Введем вместо плотности слоя <р вспомогательную функцию %, связанную

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

__ Я» Ф — % соз ос.

Элемент Й2 поверхности 2 заменим элементом телесного угла йо)

асо = Л 2 соз а .

Р'2

Элемент массы слоя будет равен

Ф<1'^ = П2%<1(о,

а интегральное уравнение (VIII.68) примет вид

2л% соз2 ос =

0 +

^

йсо +

 

 

о

 

Ь>

0)

 

 

Выразим возмущающий потенциал Т через функцию %

(VIII.70)

М. С. Молоденский предложил следующий метод решения. Поверхность 2 преобразуется в поверхность 2 так, что угловые координаты в полярной системе координат сохраняют свое первоначальное значение, а радиус-вектор (см. рис. 51) изменяется в соответствии с выражением

р = В + к(р~В) = В+Ш7,

(VIII.71)

204'

где р. — радиус-вектор преобразованной поверхности 2; к — постоянный коэф-

4пцпент. Очевидно, что при к =

1 поверхности 2 и 2 совпадают, а при к = О

преобразованная

поверхность

2

совпадает с поверхностью сферы. В

общем

случае 0 ^ к

1.

 

 

 

При переходе от 2 к 2

элемент йсо не меняется и аномалии 8§0

будем

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

Г =

со

 

 

(VIII.72)

 

г

 

 

где г — расстояние ЫК между проекцией данной точки N с радиус-вектором

р0 = В + Щ и проекцией текущей точки К с радиус-вектором р' =

Е -+- Я г

на поверхности 2 .

 

 

 

 

Новая плотность % на поверхности 2

определится уравнением,

аналогич-

ным (VIII. 69)

 

 

 

 

2Яу„ соз2а = 6 * + ДО

I А» + Д»ДО*

* » +

 

СО

Г

0)

г

 

+

 

 

 

(УШ.73)

со

 

г

 

 

где а — угол между нормалью к поверхности 2 и радиус-вектором р0 в точке N.

Величины г и ос являются функциями к и могут быть представлены в виде

рядов по степеням к.

Для поверхности сферы (см. рис. 51) имеем

г0 = 2Е з т у ,

где — угол между р_о_и р'.

Для поверхности 2 можно написать

 

г2 = р« + Р72 -

0р' соз 41) = (Е + кН1)2

+

 

+

(В + кНг)2 -

2 (В + кЩ) (В + кНг) соз г}>.

 

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

 

 

 

 

г2 = г

\\ + к я?' + яо

+ к2

+ к2 {ВТ -

Но)2.

 

 

Я

'

Я 2

 

 

Если пренебречь малыми порядка

НУЧУ

 

 

— т о

 

 

 

 

 

Я

 

 

 

г2 = г% + к2

Г-НЪ)2 = г% (1 + к2*2),

(VIII.74)

где

и = н^'-Щ Го

Это упрощение вполне оправдано, так как само интегральное уравнение (VIII.73) имеет относительную погрешность порядка сжатия Земли а вследствие того, что за отсчетную поверхность принята сфера.

205'

Наклон поверхности 2 изменяется пропорционально к, т. е. будет име* место равенство

1,§а = к%%<%,

откуда

соз2а = (1 + к2а)"1 — 1 — к22а + ... (VIII.:'

Величины Т, % также могут быть представлены в виде рядов по степеням со

т = 2

кпт„,

(VIII.?

п=0

 

 

Х = 2

Ъп%п.

(У1П.7-

71=0

 

 

Неизвестные функции Тп и %„ полностью определяют искомые на поверх ности 2 величины Т и поскольку при к = 1 ряды (VIII.76) и (VIII.77) дак:

т = 2> тп,

(УШ.7Ч

71=0

 

со

(VIII.79

ЗС = 2 Ъп.

Чтобы определить функциональную зависимость между Тп и %п, подстави: (VIII.76) и (VIII.77) в формулу (VIII.72), кроме того, функцию 1/г представим в виде ряда по степеням к, использовав для этого соотношение (VIII.74)

4 = А- (! +

= — ( ' 1 — 1 к2*2 + 4 к Ы * - . . . ) .

'

" '

О

 

о

/

В результате

 

 

 

 

 

со

 

 

 

т =

2

ЪпТп2^(%1.

+ к%1 +

к*Ъ+...)

''о

 

Пе»0

 

 

 

 

 

 

Множители при к" слева и справа должны быть равны, т. е.

г . - л ' И - й - * » '

га

"

"

(VIII.80

йсо.

(УШ.81>

Следовательно, для определения возмущающего потенциала по формуле (VIII.78) нужно знать %п. Подставим в интегральное уравнение (VIII.73) полученные ряды (VIII.75), (VIII.77) и (VIII.80)

2я (Хо + Нх + к2Хг + • • •) (1 - к21ё2 а +...)

=

206'

2 • 55 (Хо+ кгх + /с2Ха + . . . ) &2х8 ( 1 - 1 /с2*2 + . . . ) Ао.

Перемножая ряды, найдем

 

 

 

 

 

 

 

2я%0 +

2лЛх1 +

2я/с2х2 + . . . — 2я/с2Хо а + . . . =

 

 

зя

 

/СХ1 I

кгЪ

 

27-0

о

. . . ^ й© +

 

 

 

 

''о

'"о

 

 

 

 

 

 

 

 

 

+ Д'

 

 

,

АЗхХг

•Ь •

3 АЗхЗхо

 

. . . )й<в +

 

2

2

г1

 

'

 

 

 

 

№хо , кЯ*2Ул |

 

2

 

 

 

 

 

 

го

 

г0

 

 

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

 

 

 

Л1

О

 

АО

I о .1.1

-щг

Хо Ло + 2я 1^2а Хо.

 

 

• м

 

0)

"

 

 

 

 

пли в общем виде

 

 

 

 

 

 

 

 

 

=

(VIII.82)

где

2 (Шр-Цр)

Хо" йсо,

(VIII.83)

- Т И Ь

Г3Я °) 2

+ 2яХо а,

207'

Интегральное уравнение (УШ.82) подробно рассмотрено в § 44. Аналогично (VIII.44) и (VIII.45) напишем общие решения уравнений (УШ.82) в следующем виде:

 

 

Я"

 

 

(VIII.84»

 

и

 

СО

со

(VIII.85!

 

 

Подставляя (VIII.85) в (VIII.81), получим уравнения для определения Т- всех степеней

= 1 г 11 6 ^

- 1 ] йсо +

 

 

со

 

 

 

 

Г . - ^ у о * ™ - ! ] * .

 

 

( у ш а д

 

®

 

Главную часть возмущающего потенциала Т на земной поверхности дает

первое уравнение, соответствующее к =

0, и, следовательно,

а = 0, р = В.

Подставляя значение

из (VIII.83),

получим

 

 

е>

что совпадает с формулой Стокса (VIII.47).

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

Наибольшую поправку к решению Стокса дает возмущающий потенциал Т

Если ограничиться только двумя членами Т0

и Тх, то

 

 

Т = 2 0 ~ И0) + Л ДО {§ - у + г) [5ДО_ 1] А» +

А ь М . .

( уШ .87)

 

ю

 

 

 

Для

вычисления поправочных членов

6 ^ ,

. . . служат

формулы

(VIII.83).

Входящие в правые части этих формул подынтегральные выражения

убывают

как 1/гВ или даже быстрее. Поэтому при вычислении

. . .

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

208'

Подставив значение Т (VIII.87) в (VI.19), получим формулу Молоденского. лля вычисления аномалии высоты

4Кя у $$ а - V •+ Ш [8 <ф) -1] Ао +

. (УШ.88)

Условием существования решения интегральных уравнений типа (VI 11.82)- является требование, чтобы свободный член не содержал сферической функции первой степени, т. е. применительно к формулам (VIII.87) и (VII 1.88) аномалии силы тяжести должны удовлетворять условию

 

Й

+

с03'фЖ» = 0.

 

(VIII.89)

Подставив значение функции У 1 (6, %) из

(VIII.48)

в (VIII.88)

и полагая

ввиду малости величин х0,

у0,

г0) }М/В2

у,

получим

 

 

^ _ И^о

Цо_

^

с о 8 ^ с о з

 

^ с о 8 % 8 } п

81П В] +

 

 

+ д

 

 

 

 

 

 

(УШ.ЭО)

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

для определения

аномалии высоты, помимо

аномалий:

силы тяжести, необходимо знать четыре постоянные величины

 

Уо, 2о»

 

 

Д

Л

(8 —Т + бй)Ао

 

 

 

Алу

 

§ 47. ОПРЕДЕЛЕНИЕ СОСТАВЛЯЮЩИХ УКЛОНЕНИЯ ОТВЕСА НА ПОВЕРХНОСТИ ЗЕМЛИ

Для определения составляющих уклонения отвеса \ и т) в точках физической поверхности Земли необходимо слагающие возмущающего потенциала Т — величины Т0, Тг, Г2, . . . подставлять в формулу (VI.8).

Слагаемые в | и т), соответствующие потенциалам Т0,

Тх,

Т2, . . ., обозна-

чим соответственно

через

^

|2> • • •> Ло» т1х> Лг» • • • Вычислим сначала,

величины

 

*

1

дТ0

 

 

 

 

 

 

 

 

5 0

ро?

дВ '

 

 

 

 

_

вес В

дТ0

 

 

 

1

1 0

Р ^ Г

ЭЬ '

 

 

Дифференцируя формулу Стокса (VIII.47) по переменной Ву получим

 

 

В

д

 

 

 

 

4яр0

•^•ЭДог-т)-^)*»-

 

 

Д2р0у

дВ [«о 008

ВСО&1 + У о СОЗ В 81П Ь + 20

81П

В\

Будем считать

 

 

 

 

 

 

 

р о - 1 '

 

 

 

14 Заказ 1379

209?

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