shimbirev_b_p_teoriya_figury_zemli
.pdfКоординаты В ж Ь входят в неявном виде в функцию Стокса 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'
Введем обозначение
() = ~ 2у |
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 р о ( Я * ' - Я г ) + р ' ( Я Г - Я У ) - р о ( Я * ' - Я У ) |
|
|||||
2ро/"3 |
|
2р0гЗ |
|
|
2р0гЗ |
|
— |
|
_ Я ^ - Я У |
р ' - р о ) |
|
Я ^ - Я ? , |
|
^ Т Т Т Й - |
|
||
— |
* г» |
"" |
2^75 |
— |
Г |
" |
^ |
' |
Используя обозначение
Я
перепишем уравнение (VIII.66) в более простом виде
( У Ш . 6 8 ,
Введем вместо плотности слоя <р вспомогательную функцию %, связанную
с ф соотношением,
__ Я» Ф — % соз ос.
Элемент Й2 поверхности 2 заменим элементом телесного угла йо)
асо = Л 2 соз а .
Р'2
Элемент массы слоя будет равен
Ф<1'^ = П2%<1(о,
а интегральное уравнение (VIII.68) примет вид
2л% соз2 ос = |
8§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 - |
2р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 — к21§2а + ... (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) + Л ДО {§ - у + 6ег) [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? |