shimbirev_b_p_teoriya_figury_zemli
.pdfность Н = 1Ц = сопвЪ и учесть их вертикальный градиент, то выражения
(VIII.115), (УШ.116) и (VIII.117) будут более точными.
Л. П. Пеллинен вывел формулы, которые с точностью порядка [(^ — у)Н\ /
позволяют получить тождественный результат с тем, как если бы при вычислении в поле полных топографических аномалий применять формулы М. С. М-:- лоденского в первом приближении. Формулы Л. П. Пеллинена имеют виг
С =
2Л Л
И |
т- + 2 л / б Я * + V ) 5(!>) Лв + А?р, |
(VIII. 121 |
(О |
|
|
* : = - 4 огоИ ~ у ) н - т - р - + 2 я / б я т + ® с о 8 ( у п т - 1 2 2
2Л Л |
|
|
|
|
|
|
|
|
|
Т1 = ~ 4 Г 1 |
^ К г - Т к т . р . + |
гя/бД' + |
б ^ ^ ^ в т Л Й Л ^ |
+ Алр, |
(VIII.1231 |
||||
о о |
|
|
|
|
|
|
|
|
|
выражения |
в квадратных |
скобках — поправка за вертикальный |
градиент |
||||||
аномалий силы тяжести, которая вычисляется по формуле |
|
|
|||||||
|
к |
Д(У-У)н.Т.р. _ |
Ж* |
Г Г |
Уп. Т. р. |
(ё |
У)п. Т. Р. |
|
|
|
|
АН |
|
2зт ^ ^ |
/*р |
|
|
||
Здесь |
к = |
Ну' — Ни — разность |
высот |
текущей |
и |
исследуемой точек, |
|||
в которой определяют I,, Е и т|. Индексом 0 отмечены величины, относящиеся |
|||||||||
к точкам, в которых определяют |
По ряду соображений формулы (VIII. 121), |
(VIII.122) и (VIII.123) удобнее, нежели те, в которых используются поправки 6^1. При получении члена 8§" не требуется обращаться к карте высот физической поверхности Земли, достаточно знать лишь высоты этих точек; отпадает необходимость в вычислении поправки за наклон в уклонение отвеса.
§ 49. УРАВНЕНИЯ ГРАДУСНЫХ ИЗМЕРЕНИИ
Высоты квазигеоида С и составляющие & и г) уклонений отвеса относительно принятого референц-эллипсоида вычисляют по формулам (VIII.90) и (VIII.91), которые содержат четыре постоянные величины: 1У0 — 170, х0,
Ув) На-
значит, может быть поставлена обратная задача: определить эти четыре постоянные в пунктах, в которых из наблюдений получены \ и т). Градусные измерения совместно с гравиметрической съемкой позволяют решить эту задачу. В астрономических пунктах высоты квазигеоида I, определяются методами астрономического или астрономо-гравиметрического нивелирования, а составляющие уклонений отвеса вычисляются по формулам (VI.5)
Т] = (А,— Ц созВ.
Высоты квазигеоида и составляющие уклонения отвеса, вычисленные через аномалии силы тяжести по формулам первого приближения (см. § 48), обозначим: I,, г).
220'
Согласно |
(VIII 90) устанавливаем |
|
|
|
|
ъ = |
|
+[зг0со5Дсо5^ + у сое 5 81пХ + 20 |
8тБ]. |
(VIII. 124) |
|
или |
|
|
|
|
|
I _ |
2= |
+ аг0 с08 Я С08 X + |
г/0 соз Я з т Ь + г0 |
з т В. |
(VIII. 125) |
В уравнении |
(VIII.125) левая часть, |
т. е. величина |
I, — |
известна. |
|
В правую часть входят неизвестные |
|
|
|
||
|
|
ТРо-^о. хо. |
го- |
|
|
Уравнения (VIII.125) можно составить для всех пунктов градусных изме-
рений, в которых известны разности ^ — При достаточном диапазоне изменения широт и долгот все четыре неизве-
стные определяются из решения по способу наименьших квадратов. Условия
для определения ]У0 — 110, х0, |
у0, |
г0 можно |
получить по разностям составля- |
|||
ющих уклонения отвеса. |
|
|
|
|
|
|
Дифференцируя (VIII. 124) по широте В, найдем |
|
|
||||
-Ц-=-Ц— |
1УВ —Хъ&тВсо&Ь |
— у0а\пВ&тЬ |
+ гсо8 |
В. |
||
Высота квазигеоида |
и составляющая уклонения отвеса в меридиане % |
|||||
связаны соотношением |
|
|
|
|
|
|
|
ЭС . |
% |
(М+НУ)дВ' |
|
|
|
е |
дх |
|
р0дВ |
|
|
|
где М — радиус кривизны меридиана. |
|
|
|
|||
Таким образом, |
|
|
|
|
|
|
следовательно, |
|
|
|
|
|
|
(М + №) (1-1) = х0а1пВсовЬ |
+ у0 вшД 8шХ - з 0 со8Д + |
|
||||
|
|
|
|
|
|
(VIII.126)- |
Уравнение (VIII. 126) |
используют только для определения |
координат |
||||
центра инерции Земли — х0, |
у0, |
г0, так как величина |
— !У0 входит в это |
|||
уравнение с малым коэффициентом (порядка сжатия Земли). |
|
Дифференцируя (VIII. 124) по долготе Ь, получим условие для определения неизвестных по разностям составляющей уклонения отвеса в первом вертикале.
Имеем |
|
|
|
|
= |
-Ц |
х0 соз В 8 т Ь + у0 соз В соз Ь, |
||
Ол-1 иЛ/ |
|
|
|
|
но |
_ |
|
ас |
а? |
а? |
|
|||
"Л — —ду ~ |
р0соз ВдЬ |
(Ы+НУ)совВдЬ, ' |
||
где N — радиус кривизны первого вертикала. |
||||
Отсюда |
|
|
|
|
|
0 |
= |
+ |
созЯ-т! |
|
|
|
|
221
и |
аналогично |
|
|
|
|
|
|
= — |
+ |
соз В л. |
|
|
Таким образом, получим |
|
|
|
|
|
(И + |
Ш) (Г) - |
л) = Х0 |
з т Ь — у0со&Ь. |
(VIII. 12" |
|
Неизвестные х0 и у0 |
целесообразно |
определить из уравнений |
(VIII. 127 |
|
г0 |
— из уравнений (VIII.126) и из уравнений (VII 1.125) — величину |
И7,, — С |
|||
Возможны также другие комбинации исходных данных. |
|
||||
|
Знание постоянных |
— Н0, |
х, у0, г0 дает возможность определить форму |
||
Земли и ее внешнее гравитационное поле. |
|
|
|||
|
§ 50. НОВЫЕ СПОСОБЫ РЕШЕНИЯ ЗАДАЧИ М. С. МОЛОДЕНСКОГО |
||||
|
Трудности, с которыми встретились |
при практическом применении фор- |
мул, полученных М. С. Молоденским, заставляют искать новые способы определения возмущающего потенциала Земли. При этом наметилось два направления: одни исследователи пытаются дать точное решение этой проблемы, другие — приближенное, сколь угодно близкое к точному.
Точные способы решения задачи Молоденского характеризуются следующими основными особенностями: 1) граничные значения предполагаются заданными непрерывной функцией на физической поверхности Земли; 2) решение представляется в виде разложения различных выражений для возмущающего потенциала по малому параметру.
М. С. Молоденский при составлении интегрального уравнения представил возмущающий потенциал Земли Т в виде потенциала простого слоя (УТИ.бЗ).
Это — не единственный способ представления возмущающего потенциала. Можно, например, представить Т в виде потенциала двойного слоя или в виде
линейной комбинации потенциалов простого и двойного слоев, подобрав соответствующее соотношение между их плотностями.
Для поисков новых интегральных уравнений, решающих поставленную краевую задачу, можно использовать более простые краевые задачи, чем данная, а именно краевые задачи с условиями на границе
= А;
где II ж V — вспомогательные, гармонические вне поверхности Е функции. В. В. Бровар показал, что решение краевых задач можно свести к решению одного интегрального уравнения для определения плотности простого слоя 11. Если искать решение внешней задачи Дирихле в виде производной потен-
циала простого слоя по р, умноженной на р |
|
ор л2д г |
(VIII. 128) |
|
•222
то для определения (х используется краевое условие, которому на поверхности Е должна удовлетворять функция II
= А = -2яцсоз(7г, р ) 4 - | |
(УШ.129) |
Р |
2 |
|
Найдем решение внешней краевой задачи с условием на границе
|
|
дУ |
|
|
|
Эр 2. = / » |
(VIII.130) |
где V (р, В, Ь) — искомая гармоническая функция. |
|
||
Положим, что вне поверхности |
2 |
|
|
= |
= р |
[с учетом |
(VIII. 128)], |
отсюда |
2 |
|
|
|
|
|
|
|
|
|
(VIII. 131) |
|
|
2 |
|
Подставляя (VIII.131) в (VIII.130), найдем |
|
||
|
|
д±_ |
|
/2 = |
— 2я|л соз (п, р) + |
(VIII.132) |
|
|
|
2 |
|
Формулы (VIII.132) и (VIII.131) решают поставленную задачу. Отметим, что интегральное уравнение (VIII.132) совпадает с уравнением (VIII.129).
Пусть, наконец, нужно решить внешнюю краевую задачу с условием на границе
2Т |
. дТ |
= Л , |
(VIII.133) |
Р |
др 2 |
где Т (р, В, Ь) — искомое решение задачи.
Положим, что функция Т связана с функцией II (VIII.128) условием
Т |
< р т > = : + Р - з ™ и -=р Ж И ^ |
( У Ш Ш ) |
|
2 |
|
Подставив в краевое условие уравнение (VIII.133), получим интегральноеуравнение для [1
Р |
= —2я|*сов(л, Р) + | |
(VIII.135) |
2 |
|
|
|
|
Для определения функции Т (р, Ва Ь) воспользуемся формулой (VIII.134)
1
д
дР
223'
или после интегрирования по р |
|
|
|
|
|
|
Г(р, в, ^Н-^г |
|
|
|
|
ь) • |
т п л з б ! |
|
2 |
|
|
|
|
|
Рассмотрим способ определения возмущающего потенциала Земли, пред- |
||||||
ложенный В. В. Броваром. |
|
|
|
|
|
|
Пусть произвольная точка К земной поверхности |
2 (см. рис. 51) имеет |
|||||
радиус-вектор р' = Я + Ну', |
а внешняя точка М, в |
которой |
вычисляется |
|||
возмущающий потенциал Т, радиус-вектор р = В |
Н. Через В', |
Ь' обозначим |
||||
Координаты текущей точки К, |
а через В, Ь — точки |
М. |
|
|
||
Расстояние г' между точками К и М |
определяется соотношением |
|||||
г'г |
= р2 + |
р'2 - |
2рр" соз 1)). |
|
|
(VIII. 137) |
Когда внешняя фиксированная точка сливается с точкой |
то ее радиус- |
|||||
-вектор р0 = В + 'Щ, а расстояние |
г между точками N ж К |
|
||||
г2 = р2 + |
р'2-2р'р0со81(3. |
|
|
(VIII.138» |
||
Проекцию г на сферу радиуса В обозначим г0 |
|
|
|
|||
|
г0 = 2Взт^-. |
|
|
(VIII.139) |
||
Рассмотрим решение внешней |
задачи Дирихле |
для некоторой функции |
||||
7 ( Р , Я , |
|
|
+ |
|
|
( у П Ш 0 ) |
|
2 |
|
|
2 |
|
|
Первое слагаемое правой части представляет собой потенциал простого слоя с плотностью (1/4я) р.'; второе — произведение р на производную потенциала простого слоя с плотностью (1/2:гс) (х'. Оба слоя распределены по поверхности 2. Функция V (р, В, Ь) является гармонической во внешнем пространстве
ирегулярной на бесконечности.
Всоответствии со свойствами производной от потенциала простого слоя [см. формулу (1.85)] второе слагаемое на поверхности 2 терпит разрыв, вели-
чина которого будет равна
— 2я |
ц.' соз (п, Ро) р0 = — р.' р0 соз (ге, |
р0), |
|||
где п — направление внешней к поверхности |
2 нормали. |
||||
Формуле (VIII.140) придадим иной вид, используя соотношение (VIII.137). |
|||||
Поскольку |
|
|
|
|
|
д г'— _ |
1 |
дг' |
дг' |
р —р'созф |
|
д р |
г'г |
др ' |
др |
г' |
' |
легко доказать, что |
|
д±_ |
|
|
|
|
|
|
|
|
|
|
± + |
|
|
|
(УШ.14П |
224'
Подставляя полученное соотношение в (VIII.140), получим |
|
||||
|
|
7(р, В, |
= |
|
(VIII.142) |
|
|
|
|
2 |
|
На поверхности |
2 (при р -> р0) |
функция V (р, В, Ь) должна |
удовлетво- |
||
рять уравнению |
|
|
|
|
|
Г(р0, |
В, |
Ь) = ^ ' с о з К |
Ро)Ро + 4 г 1 ^ ' |
(VIII. 143) |
|
|
|
|
|
' 2 |
|
Интегральное |
уравнение |
(VIII. 143) позволяет определить функцию ц/, |
|||
если на поверхности 2 задано |
V (р0, В, Ь). Подставляя найденное значение р/ |
||||
в (VIII. 142), получим решение внешней задачи Дирихле. |
|
||||
Отметим, что |
поскольку |
р' = р0, соз (п, р0) = 1, уравнение |
(VIII.143) |
||
хтя сферы принимает простой вид |
|
|
|||
|
|
У(р0, В, |
/,) = - р # ' . |
(VIII.144) |
|
Применяя тот же метод, рассмотрим задачу определения возмущающего |
|||||
потенциала Т. |
|
|
|
|
|
Введем функцию |
|
|
|
||
|
|
V (Р, |
= |
Г Л = 2Т + Р ^ - . |
(VIII.145) |
Эта функция гармоническая, регулярная на бесконечности и на поверхности 2 в соответствии с (VIII. 10) принимает значение
V (Ро, В, П) = -р0 б?0 .
Подставляя последнее равенство в (VIII.143), получим интегральное уравнение для определения [х'
|
|
соз (гс, |
= |
+ |
(VIII.146) |
|
|
|
|
|
2 |
Для нахождения Т подставим (VIII.145) в (VIII.140). Получим |
|||||
|
Р 7(Р, В, |
= |
|
= |
+ |
|
|
|
|
|
2 |
|
Интегрируя по р, найдем |
|
|
|
|
Г(р, |
4п |
и ' |
|
|
^ 2 + Г х { В : Ь) • (VIII.147) |
|
2 |
|
|
|
|
Внутренний интеграл в формуле (VIII.147) можно представить
15 Заказ 1379 |
225 |
При вычислении второго интеграла применим формулу иптегрированн* по частям
| и йи = ил) — | V Ли,
положив |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ц = 2р2; |
|
= |
|
|
|
|
|
|
||
|
|
|
|
йи = 4рф; |
V = |
|
|
|
|
|
|||
получим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
д |
± |
|
|
|
г |
|
|
|
|
|
|
|
|
|
|
|
|
|
и |
^ . |
|
|
(VIII.14- |
Таким образом, |
|
|
|
|
|
|
|
V |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
± г( . |
|
|
|
|
|
|
|
|
|
|
|
|
|
Р2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
Обозначив |
|
|
|
|
|
|
|
|
|
|
|
|
|
получим |
7-(р, |
В, |
Ц = |
|
|
|
|
|
|
|
|
( V I I I . I V |
|
Для |
определения |
функции Е' |
воспользуемся известными формулами ин |
||||||||||
тегрального исчисления (VIII.33) и (VIII.34). Полагая |
|
|
|
||||||||||
получим |
|
у = р; |
а — 1; |
Ь = — гр'соз'ф, |
с — р'*, |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
^ _р_йр = г , + |
р , с о 8 |
^ ^ |
= |
г , + |
р , с о з ^ 1 п |
+ р _ |
р/ с о 8 |
^ |
(VIII.15. |
||||
Следовательно, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Е ' |
= |
^ + |
^ |
^ |
1п (г' + |
Р - |
р'.соз о])) |
|
|
||
и |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
-у,— |
ЗЕ' — -р |
|
•— ^ р Г ^ |
1п(г- + р - р ' с о з ^ ) . |
(VIII. 15: |
|||||||
Полученную В. В. Броваром функцию удобно выразить через обобщенну |
|||||||||||||
функцию |
Стокса. Предварительно |
представим обобщенную |
функцию Стоке |
||||||||||
в виде |
|
|
|
00 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
п=2 |
— 1 |
(."+1 |
|
|
/•' |
1 р |
р2 |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
5р'совг[) |
Зр'соа^ |
^ т-' + р —р' собг|з |
|
(VIII I"-' |
||||||
|
|
|
|
р2 |
|
р2 |
|
|
|
2р |
' |
|
|
Функция (VIII.152) при р' = К совпадает с обычной обобщенной функппе> Стокса.
226'
Сравнивая (VIII.151) с (VIII.152), найдем
А _ З Я ' = * ( Р , |
+ |
Ы 2 Р . . |
Чтобы решение было возможно, необходимо исключить негармонические члены в выражении для возмущающего потенциала, т. е. удовлетворить условие
- З у ^ ц ' р ' с о з г Н ^ О
2
•ли
|ц'р'созг|)й2=0. (VIII.153)
2
Следовательно, в формуле (VIII. 149) сферические функции первой степени будут отсутствовать.
Поэтому окончательно
Т(Р, В, |
= |
[ я ( р , |
р ' ) - ± ] ^ + Г 1 ( р Д 2 ' ь) . |
(VIII.154) |
|
|
|
* 2 |
|
|
|
Введем новую плотность [г, определяемую из соотношения |
|
||||
|
|
соз (п', |
р') |
д2 |
|
|
|
— • |
|
||
|
|
|
|
р |
|
Обозначив через йа проекцию элемента |
йЕ на сферу радиуса |
В, будем |
|||
аметь |
|
<*Зсо8 {п', |
р') Д2 |
|
|
|
|
|
|||
|
|
р'* |
|
|
|
Элемент массы слоя
р/ йЕ = (г д,а
• основные формулы (VIII.146), (VIII.153) и (VIII.154) в этом случае примут вид
^ - [ Л С О З = |
+ |
(VIII.155) |
|
|
|
о |
|
I |
/ (Яр'соз Лр йа = О, |
(VIII.156) |
|
а |
|
|
|
Т(Р, В, |
^ С * 5 ^ |
|
(VIII.157) |
о |
|
|
|
Применяя метод Молоденского, Бровар получил соотношения для определения величин (х и Т
о
И , = Ц „ |
( * , Р„) + |
П ^ |
+ |
I |
& |
|
|
а |
|
|
а |
15* |
227 |
1 С Г — Ну
а
3 |
ГР |
(Я^'-ЯУ)3 |
|
1 |
Г Г |
(ЯУ- |
Я*)' |
|
- ^ г З |
о |
3 |
а |
а + ^ я |
3 3 111 |
|
<*а, |
|
|
|
|
|
|
а |
|
|
|
|
|
|
= |
(*) - 1 1 |
+ |
|
|
|
|
|
|
в |
|
|
|
|
|
|
|
|
о |
|
|
|
|
|
г 2 = |
|
1 |
|
- |
И |
^ |
( Я Т '^Я °, ) 2 |
|
4яЯ Я М * < + ) - 1 1 |
||||||||
|
|
|
а |
|
|
|
|
|
= |
4Й- И «*» ^ <+> - 1 1 ^ |
- 5 Г |
И ^ |
" ^ |
Г ^ |
Эти выражения проще аналогичных формул Молоденского, поскольку интегральное уравнение решается при помощи одних лишь функций р„, тогд; как в решении Молоденского приходится вычислять две функции 8@„ и х- В. В. Броваром рассмотрен также случай, когда возмущающий потенциал представляется суммой интегралов с ядрами Е и ЛЕ/йп'. Здесь п' — внешняя нормаль к поверхности 2 в текущей точке, а Е — фундаментальная гармони-
ческая функция, имеющая вид
|
Е { р, % |
|
= |
|
|
г ' + Р - Р ' с о з ^ |
|
|
(VIII.158) |
||
Не делая подробного вывода, опишем его основные этапы. Найдем решение |
|||||||||||
внешней задачи Дирихле для земной поверхности 2 в виде |
|
|
|
|
|||||||
|
|
|
|
|
а — |
|
|
|
|
|
|
|
и (Р, В, |
Ь) = Л |
V |
Й2 + Л |
|
й 2, |
|
(VIII.159 |
|||
|
|
|
2 |
|
|
2 |
|
|
|
|
|
где V — неизвестная плотность |
двойного |
слоя, |
распределенного |
на 2; |
р. = |
||||||
= V сов (п', р')/2р' — плотность |
простого |
слоя, |
зависящая |
от |
V, |
положения |
|||||
и наклона элемента поверхности в текущей точке. |
|
|
|
|
|||||||
Все обозначения даны в соответствии с рис. 51. |
|
|
|
|
|||||||
Если исследуемая точка стремится к поверхности 2, то |
р ->• р0, г' |
-»- г. |
|||||||||
При р = |
р0 потенциал |
двойного |
слоя испытывает разрыв, |
равный 2пу |
[см. |
||||||
(1.89)], а |
Нт II (р, В, |
Ь) |
равен заданной функции / (В, Ь). |
Сделав указан- |
|||||||
ный предельный переход, |
получим интегральное уравнение для плотности \ |
||||||||||
|
|
|
|
|
1 |
|
|
|
|
|
|
|
НВ, Ь) = |
|
П |
|
|
|
|
|
|
(VIII.1601 |
228'
После некоторых преобразований получим
О, сгг
|
|
т |
Ц ^ у |
+ ^ |
- |
^ |
Ъ |
+ |
^ ^ |
Я |
^ Ъ . |
(VIII. 162) |
|
|
|
|
СГ! |
|
|
|
01 |
|
|
|
|
|
Здесь |
йах — проекция |
элемента |
й2 |
на сферу радиуса р', йах |
= сов (п', |
||||||
|
<22, а дН'/дЗ — наклон земной поверхности в текущей точке. |
|
||||||||||
| |
Искомая |
плотность V |
находится в |
|
результате |
решения интегрального |
||||||
I сравнения |
(VIII.162). |
Решение |
внешней |
задачи |
Дирихле дается |
формулой |
||||||
' |
VIII.161). |
|
|
|
|
|
|
|
|
|
|
|
|
Воспользуемся данным методом для решения задачи Молоденского. Поло- |
'рНМ, что
~е Т — возмущающий потенциал Земли, удовлетворяющий на земной поверх-
ности 2 краевому условию (VII 1.62).
Функция Л на поверхности 2 принимает значения |
|
V |2 = -Р0 в*о. |
(У1И.163) |
Поскольку во внешнем пространстве функция II определяется по формуле VIII. 164), получим
Т Ж |
= И |
|
|
|
й |
|
Р Ш ) |
|
01 |
К |
|
|
О, |
|
|
Плотность V определяется из уравнения |
|
|
|||||
- р ^ о |
^ ^ - Ь |
^ * ! |
^ |
^ |
^ ДО |
V ^ |
(VIII.165) |
|
01 |
|
|
|
С , |
|
|
На основании (VIII.164) напишем |
|
|
|
|
|||
|
01 |
^ |
|
|
01 |
|
|
я-'и после интегрирования по р |
|
|
|
|
|
||
|
|
|
|
|
|
1 |
|
+ |
5$ •V 7 Г |
( р ' з 1 п |
* |
$ |
й р) |
+ С (р\ 115). |
(VIII. 166) |
|
01 |
|
|
|
|
|
|
Введем обозначения:
В{Р, гр, Р') = Р'8Ш^ Г - 4 - ф .