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

shimbirev_b_p_teoriya_figury_zemli

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

ность Н = = сопвЪ и учесть их вертикальный градиент, то выражения

(VIII.115), (УШ.116) и (VIII.117) будут более точными.

Л. П. Пеллинен вывел формулы, которые с точностью порядка [(^ — у)Н\ /

позволяют получить тождественный результат с тем, как если бы при вычислении в поле полных топографических аномалий применять формулы М. С. М-:- лоденского в первом приближении. Формулы Л. П. Пеллинена имеют виг

С =

2Л Л

И

т- + 2 л / б Я * + V ) 5(!>) Лв + А?р,

(VIII. 121

 

 

* : = - 4 огоИ ~ у ) н - т - р - + 2 я / б я т + ® с о 8 ( у п т - 1 2 2

Л

 

 

 

 

 

 

 

 

Т1 = ~ 4 Г 1

^ К г - Т к т . р . +

гя/бД' +

б ^ ^ ^ в т Л Й Л ^

+ Алр,

(VIII.1231

о о

 

 

 

 

 

 

 

 

выражения

в квадратных

скобках — поправка за вертикальный

градиент

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

 

 

 

к

Д(У-У)н.Т.р. _

Ж*

Г Г

Уп. Т. р.

У)п. Т. Р.

 

 

 

АН

 

2зт ^ ^

/*р

 

 

Здесь

к =

Ну' Ни — разность

высот

текущей

и

исследуемой точек,

в которой определяют I,, Е и т|. Индексом 0 отмечены величины, относящиеся

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

По ряду соображений формулы (VIII. 121),

(VIII.122) и (VIII.123) удобнее, нежели те, в которых используются поправки 6^1. При получении члена 8§" не требуется обращаться к карте высот физической поверхности Земли, достаточно знать лишь высоты этих точек; отпадает необходимость в вычислении поправки за наклон в уклонение отвеса.

§ 49. УРАВНЕНИЯ ГРАДУСНЫХ ИЗМЕРЕНИИ

Высоты квазигеоида С и составляющие & и г) уклонений отвеса относительно принятого референц-эллипсоида вычисляют по формулам (VIII.90) и (VIII.91), которые содержат четыре постоянные величины: 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

 

 

 

'

 

 

Функция (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 - ф .

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