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

shimbirev_b_p_teoriya_figury_zemli

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

Тогда

 

 

 

 

 

4 п у дБ

 

 

 

 

1

 

д

п С 0 8 В СОЗ Ь

Уо СОЗ В 81П Ь +

2 0 8111 В].

 

 

Ро

дБ

 

 

 

 

 

 

=

 

4 я у вес Я Д О

7) Я (г|>)Жв-

 

 

 

вес В

 

д

0 соз В соз Ь +

у0 соз В з т Ь + 2 0

з т 5].

 

 

ро

дЬ

 

 

 

 

Используя формулы (VIII.56), получим

 

 

1

ГГ , _

..ч й5(г|))

• соз А <2со +

 

 

4 л у

ДООГ-7)

 

 

_1_

0 81П 5 СОЗ Ь + у0 81П В 81П Ь — 20 СОЗ Б]

(У1Н.91)

ро

 

 

 

 

 

 

 

 

Ло:

4лу

 

 

 

 

81П А ЙСО + А - [х0 81П Ь — у0 СОЗ

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

МШ,Н)

поверхности

действительного

потенциала

относительно уровенной поверхности нор-

 

 

мального

потенциала,

заданного

произ-

 

вольно.

 

 

 

 

 

 

 

 

Для определения поправок к форму-

 

лам Венинг-Мейнеса (VIII.91)

необходимо

 

учесть слагаемые Тг, Т2

в возмуща-

 

ющем потенциале Т (VIII.86).

При этом

 

приходится вычислять производные

дТе0дВ

 

и зес В/р0-дТе/дЬ,

где индексом е обозна-

 

чены значения внешних предельных произ-

 

водных

возмущающего потенциала

Т.

 

При

вычислении

производных

следует

 

учитывать, что возмущающий потенциал есть

 

функция не

только

координат

В

и Ь, но

 

и высоты Н.

 

 

 

 

 

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

внешнего

пространства,

радиус-

вектор которой р = В + Н, будем обозначать Т (В, Ь,

Н).

 

 

 

Для точки на поверхности 2, радиус-вектор которой р0

= В +

Щ, потен-

циал будет Т (В, Ь, Щ). В свою очередь нормальную высоту Щ

следует рас-

сматривать как функцию координат В и Ь,

т. е. считать,

что Щ

=

/ (В, Ь).

Вид этой функции нам неизвестен. Потенциал Т на самой поверхности 2 в общем виде обозначим Т [В, Ь, Щ (В, Ь)].

При вычислении производных от потенциала йТ [В, Ь,

(В, Ь)}

 

ав

ЗЮ

••Т [В, Ь, НУ ( В, Ь

следует

пользоваться

V

 

 

производной

от неявной

" —1

 

формулой,

С1Л-1

 

 

 

 

 

 

 

 

 

 

 

 

функции,

например

 

 

 

 

 

 

 

 

 

 

АТ [В, Ь,

Щ (в, Ь)

дТ [В, ь, Щ

(В, Ь)\

• дТ

[В,

Ь,

Щ (В,

Щ дЩ

 

тТто9ч

 

с1В

 

~

дВ

 

 

 

 

дЩ

дВ

 

Рассмотрим вопрос о значении внешних предельных производных возму-

щающего

потенциала

Т,

причем вывод будем делать только

для

вели-

чины дТе0дВ,

так как производная

дТе0дЬ получается аналогично.

 

Найдем производную

от Т [В, Ь, Щ

(В, Ь) ] по направлению

I, касатель-

ному к поверхности 2

и лежащему в плоскости меридиана. Получим

 

 

 

*ту,ь.тв.щ

=

[В, ь,тв,

щ

со8 (В>

^

(У1П

93)

так как

 

 

 

 

6,1 — р0 вес (В,

I),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где (В, () — угол между касательными к сфере и к поверхности 2 в точке N (рис. 52).

Производная от Т (В, Ь, Н) во внешней точке по тому же направлению Ъ- будет

ат (в ь, н) =

ат {в ь н)

 

г)

а

р ав

4

"

Если внешняя точка (В, Ь, Н) неограниченно приближается к поверхности 2, р -»- Ро, то в пределе мы получим значение внешней предельной произ-

водной возмущающего потенциала Т по касательному к 2 направлению I

, АТе (В, Ь, Н)

АТе{В, Ь, Н) д„„,р

V

Т(

^АВ

008

1) -

«)• (УШ.94)

Возмущающий потенциал Т определяется формулой (VIII.63). Поскольку потенциал простого слоя есть функция непрерывная во всем пространстве, то

Те (В, Ь, И) = Т {В, Ь, Щ).

(VIII.95)

При этом условии должны быть равны и производные потенциалов по направлению I, касательному к поверхности 2, т. е.

ат (в, ь, Щ) _

Ате (в, ь, н)

аь

~~

а

Сравнивая (VIII.93) и (VIII.94), получим

АТ (в, Ь, Н1) _

дте (В, Ь, Н)

дТе (В, Ь, Н)

дЩ

 

ро АВ

~

ро дВ +

дН

ро дВ '

 

отсюда

 

 

 

 

 

 

 

дТе (В, Ь, Н) =

АТ (В, Ь, Щ)

дТе (В, Ь, Н)

дЩ

(УШ.96>

ро дВ

ро АВ

дН

ро дВ '

 

Используя краевое условие (УШЛО), определим

 

 

 

дТе (В, Ь, Н) _

я

е (В, Ь, И)

 

 

 

 

дН

~

 

ро

'

 

 

14*

211:

или на основании

(VIII.95)

 

 

 

 

 

 

 

дТе

(В, Ь, Н) __

,

2Г (В,

Ь, НЪ)

 

 

 

Ш

 

 

 

 

й

 

 

 

от

 

 

 

 

 

Подставляя значение - щ - в (VIII.96), получим

 

дТе(В, Ь, Н)

АТ(В, Ь, НХ)

• Г»

,

2Т (В, Ь, НЪ) -| ддУ

Ро дВ

 

роав

 

 

 

ро

}дВ'

С учетом формулы (VIII.92), окончательно найдем

 

дТе,

Я)

_

дТ (В, ь, Щ) ,

 

,

2Г (В, Ь, НЦ)

 

ро дВ

 

родВ

•[«йг

I

р0

т

 

'

[_ 60

 

 

 

дТ(В, ь,

П)1_дн]_

 

 

/ у ш 97)

 

 

+

дЩ

] р одВ'

 

 

(VIII.У//

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

На основании формул (VIII.63) и (VIII.95) получаем

 

Т{В, Ь,

=

(УШ.98)

 

2

 

где

/ • 2 = Р о + Р ' 2 - 2 Р О Р ' С 0 8 1|>.

Следовательно,

но

Ро дв

 

 

^

^

Ро

 

 

 

 

 

 

 

 

 

 

 

д

{ т )

дг

_

 

1

ег

Ро дВ

~~

 

дг

ро дв

 

 

7-2 ро дВ '

ц согласно (VIII.56)

р'япф

аг|) _

 

р'зштр

.

дг

 

 

7

 

дВ~

 

 

г

С 0 8 Л -

Используя эти соотношения, получим

 

 

 

 

 

V Г

)

Р 31

 

соз А.

 

 

Ро дв

 

 

 

Г3

 

 

 

 

Введем сюда величину г0, связанную с г|з соотношением

 

 

г0 =

2Ввт±,

 

 

имеем

 

 

 

 

 

 

 

 

ЗШ

л

Ф

 

''О

 

Ф

= 2 31П

СОЗ -у =

 

соз -у .

212'

Следовательно^

 

г0

Ф

л р'

ЯП = - + 0 0 8 -^-СОЗ А^—

ро дВ

г®

2

Я

2

На основании формулы (VIII.65) получим прямое значение производной

дТ(В, Ь, Щ)

_

1

Г Г Ф ^у |

1

Г Г,. Р,2~Р§ лу

 

 

 

 

2ро

2

 

 

2

 

Используя точное

соотношение

(VIII.67)

и

полагая р'/Д

1, 1/2р0 ^ .

1 2Н, что вносит ошибку порядка Ну /В, найдем

 

 

дТ

(В,

Ь,

Ц

Ф^|-соз^-СО8 4ЙЕ,

(УШ.99)

 

Ро эв

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

 

 

 

 

 

 

7' - Я0У)2

 

( У Ш . Ю О )

 

 

 

'2

 

 

 

 

 

 

 

 

 

 

 

 

Подставим (VIII.98), (УШ.99) и (УШ.ЮО) в формулу (УШ.97)

дТе(В, Ь,

И)

 

 

 

РО ^

- И Ф -Й-С08 Т С03 ^

+

[б*о +

_1_

 

 

 

 

2 В

2

2

2

^

 

П'Хгле преобразований получим

РО дВ

Вместо поверхности 2 примем, как и раньше, поверхность 2, которая ••ределяется преобразованием радиус-векторов согласно (VIII.71). В этом сгучае необходимо иметь в виду соотношения

 

р'-р0 =

к(нг-нъ),

(Р'-Ро)2 2

( Б Г ' - Щ У

р ' 2 - р * _

к ( н ^ - Щ )

, к* ( я у ' - я у ) 2

0гЗ

гз

0гЗ

213'

Кроме того, введем новую плотность % и элемент телесного угла йсо. Элемент массы слоя будет

=ЯЧ<1со.

Для поверхности 2 будем иметь

 

 

 

 

 

Я ' * ' 1 ' Н )

=

г з

соз 1

соз А (1<о +

 

 

Ро дВ

 

2

1

 

 

 

г»

 

 

п с г у И я т ' - я 0 ? ) 2

ЙС0

X

 

 

 

гЗ

 

 

X А; Ро ОВ

Подставляя значение 1/г из формулы (VIII.80) в бесконечный ряд (VIII.77), получим

Ро ° В

,) «3

го

 

2

 

1

г о

 

Л» +

 

 

 

 

*

, ЗЯГГ

2* " Х »

1

 

г,

,

й(о +

 

 

1 1

 

п»

 

 

+ •

 

+ д . ^

гл - чс, (

^

Ж

|-4 _

3. ^

( я У - я Д У +

_ (1а -

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

дН*

(УШ.Ю1;

 

 

 

 

 

 

 

ро дВ

В соответствии с (VIII.76) можно представить

 

 

 

 

 

 

Я)

со

Г

Ь,

Я) -1

 

 

 

дТе(В,

Ь,

 

 

 

 

ро дВ

 

^

I

ро дВ

V

 

 

 

 

 

 

п=о

 

 

 

 

 

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

 

 

+

 

 

 

 

ро

дВ

 

Г

дТе (В, Ь, Я )

П

0 2

Г Г Х2

соз

^

л ?

(VIII. 102^

I

я р

_|2

= Н

\ \

2

соз А й(0 —

 

 

Ро 9В

 

Л ,)

 

 

 

 

 

| Хо

 

соз

соз Л Ло +

 

 

 

н>

 

 

 

 

 

 

 

+

214'

Заменим интегралы с обкладками %п на интегралы с обкладками

8%п, для

его используем соотношения (VIII.85).

 

 

 

 

Дифференцируя их по Б, получим

 

 

 

 

 

5

 

! д 2

 

 

 

 

 

 

Ро дВ

 

ф

1

 

©

 

 

ЯП

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

\ Го /

 

1

дг0

 

 

 

 

 

Ро дВ

 

ГI Ро дВ •

 

Поскольку 7*0 =

2В 81П -у-, то

 

 

 

 

 

 

 

 

 

дг о

В

\|)

 

 

 

 

 

 

р0дВ ~ро

 

С08ТЖ-

 

Полагая

= 1 и Ц

= —соз А, получим

 

 

 

 

 

 

дг0

 

Ф

соз

л

 

 

 

 

 

ТЕ» —=

—соз^ул о

А,

 

 

 

 

 

Ро дВ

 

2

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

— )

соз-^-соЗА

 

 

 

 

 

\ Г о /

 

2

 

 

 

 

 

 

1 дВ

 

 

 

 

 

Поэтому

 

 

 

 

 

 

 

 

 

 

д

 

В2

^ - ^ - й ш ! = Л2

Л

-§-соз А соя 4 А»;

 

 

Ро дВ

 

 

 

со

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В

 

 

 

 

 

 

 

 

Ро дВ ) 4л

со

 

 

)

 

со

 

 

 

 

 

со

 

 

 

 

 

Таким образом, получаем соотношение

 

 

 

В 2 1 5 ^ - с о з

соз А А» = -

 

^

6 * „ с о з А Ао,

(УШ.ЮЗ)

 

СО

 

 

 

 

 

0)

 

 

юторое используем для преобразования формул (VIII. 102).

Получим аналогично соотношение для случая п = 0. Используя соотно-

ваае

(VIII.85),

найдем

 

 

^ ^

ж соз 4

соз Л Ао = -

{ $ 6,0

соз 4 Ао + ^ Ш . (УШ. 104)

СО

 

0)

 

215'

Произведя в (VIII.102) замены в соответствии с соотношениями (VIII.! и (VIII.104), получим

[

дТе(В, Ь,

ш

 

 

Н)

^ г С08 Аё(0 + ж ж

 

Ро ев

•]„ *- - 4г И

 

 

 

 

 

 

дН]_

2в ) ро ад.

ЗД2

3 3 Хо"

^— С08-|-С08 4Л» +

 

дН]_

+ 0 +

5 1 ] Ро д В

Здесь

в соответствии с обозначениями, принятыми ранее.

Нулевое приближение совпадает с полученными ранее формулами ВенингМейнеса (VII 1.91). Поправки к нулевому решению будут иметь вид

у \р0дВ

Л

Ыу 3 3

61 ду

 

 

 

зт0

\ дн1

1г=~

 

 

<0

. (VIII. 106)

 

 

 

ЗД2

 

 

 

 

(//*'-

У/.Т)2

соз -у соз Ай(.о

Хо *

 

 

 

 

 

 

дН]_

Ро дВ

216'

Аналогично определяются поправки гц, х\г . . .

1 / зесВ дТ \

1

Г Г

 

 

 

У V

Ро дь

 

Л

 

о з

01

 

 

 

 

 

 

ЗГ0

вес В

дН%

 

 

 

 

 

Ро

 

дЬ

 

1 /

вес В

дТ \

1

Г Гс

(ф) . . , ,

 

 

 

 

 

 

 

 

(VIII.107)

ЗЛ2

Хо

— г — с о з

2

81П А с1(й -

 

 

 

 

 

г*

 

 

 

 

 

 

 

374 Л вес 5

(ЗЯ?

 

 

 

 

 

Ро

 

дЬ

 

Составляющие уклонения отвеса, таким образом, находятся по формулам

 

11 =11о+ "41 +

112 4

 

(VIII. 108)

 

 

 

 

Исследование на моделях, проведенное

В. Ф. Еремеевым, показало, что

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

 

и г^ уменьшает погрешности

|

п >1 от 1,5—2,2" до 0,3—0,2".

 

Следующие приближения снижают ошибку, но

менее существенно. Формулы

 

(VIII.106)

и (VIII.107) позволяют вычислять

|

п I] с ошибкой, меньшей 0,1",

независимо от аномальности района и сложности

его рельефа.

 

 

 

 

 

 

 

 

§ 48. ФОРМУЛЫ ПЕРВОГО ПРИБЛИЖЕНИЯ

Приведем формулы М. С. Молоденского первого приближения для общего земного эллипсоида. На основании (VIII.87) получим, что возмущающий потенциал на физической поверхности Земли будет вычисляться по формуле

(VIII. 109)

Поправка в аномалию силы тяжести, учитывающая влияние рельефа, в соответствии с (VIII.83) имеет вид

(УШ.ИО)

Вспомогательную функцию %0 на основании (VIII.84) вычисляют как

Х о

- V , 3

(VIII.111)

 

217

Учитывая вышеизложенное, а также формулы (УШ.88), (VIII.91 (VIII.106), (VIII. 107) и (VIII.108), при вычислении высот квазигеоида и сост: вляющих уклонения отвеса можно пользоваться формулами

 

Б = - Й Л $ От -

V + «й) Я (ф) Ж».

(VIIIЛ1-

 

а>

 

 

 

 

ш

 

 

 

 

 

1 С Г /

, , > й » ) .

, ,

1 /

.

\ вес В

ю

 

 

 

 

 

 

 

 

 

 

(VIII.114;

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

достаточно сложными ввиду большой зависимости поправок

от рельефа.

По расчетам Л. П. Пеллинена поправки

 

по своей величине могут превы-

шать аномалии силы тяжести. Исследования Е. М. Орловой показали, что для получения надежного результата эти поправки необходимо вычислять для всех экстремальных точек физической поверхности.

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

Поправки 8%2, . . . определяются в основном рельефом местности, непосредственно окружающем исследуемую точку, и при их вычислении не требуется производить интегрирование по всей поверхности Земли, однако

при переходе от 8%п _ г к 8$п необходимо предварительно

определить %„_1

по формулам (VIII.84), для чего возникает необходимость

вычислить 8§-„ _ 1

по поверхности всей Земли.

 

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

Так, например,, уклонения отвеса по формулам нулевого приближения в равнинных районах получаются со средней квадратической ошибкой ±0,3— 0,5", тогда как в горах — с ошибкой ±1,0—1,4" х. Чтобы получить в горных районах такую же точность, как и в равнинных, необходимо тем или иным способом исключить из аномалий силы тяжести влияние топографических масс и все вычисления (высот квазигеоида и уклонений отвеса) производить в поле остаточных аномалий. Затем необходимо восстановить влияние исключенных масс непосредственно на высоты квазигеоида и уклонения отвеса. Таким образом, высоты квазигеоида и уклонения отвеса определяются как сумма двух слагаемых — влияния топографических масс и аномальных масс, которое вычисляется по формулам Стокса и Венинг-Мейнеса.

1 По исследованиям Е. М. Орловой.

218'

Наиболее удобны для вычислений формулы, предложенные Л. П. Пел-

лпненом

 

 

 

 

^ =

К ^ - Т)н. X. р. +

2л/6НЦ 8 М

йсо + де„,

(VIII. 115)

 

 

 

 

Я

 

 

 

| =

|[(?-7)н.т.р.+2я/бЯ7 ](?(ф)с08Л^Л|; + Д5р,

(УШ.116)

о о

 

 

 

Я

 

 

 

Т| = - - | г 1

{ К е - Т к т . г . + ^ ^

^ ^ ^ ^

^ + Ч -

(VIII. 117)

о о

Стоящие в квадратных скобках выражения представляют собой аномалии Фая, вычисленные методом косвенной интерполяции через аномалии в неполной топографической редукции [см. формулу (VII.29)].

Высоты квазигеоида и составляющие уклонения отвеса, вычисленные с использованием аномалий Фая, соответствуют случаю, когда топографические

массы сконденсированы на поверхности Н = Щ = соп§1, проходящей

через

исследуемую точку.

 

 

Поправки в высоту квазигеоида и составляющие уклонения отвеса,

явля-

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

| и г] при их действи-

тельном расположении и конденсации на поверхности Н = Щ, = сопз1, вычи-

сляются по формулам

(VIII. 118)

со

(VIII.119)

со

(VIII. 120)

со

 

 

где б — плотность топографических масс;

к = Ну' — Н1\ г — расстояние

между исследуемой точкой и текущей точкой физической поверхности;

г0

расстояние между их проекциями на отсчетную поверхность.

 

По исследованиям Л. П. Пеллинена член

не превышает 0,5 м,

вели-

чины А Ер и Ат)р в горах могут достигать нескольких секунд. Для вычисления поправок А Ер и Ду]р имеются таблицы, составленные Е. М. Орловой

Формулы (VIII. 115), (VIII.116) и (VIII.117) получили распространение при вычислении уклонений отвеса в пунктах геодезической сети СССР; точность вычислений вполне соответствует современному состоянию гравиметрической съемки в горных районах. Однако после исключения влияния топографических масс не учитываются уклонения физической поверхности Земли от поверхности Н = Щ = сопзЪ, проходящей через исследуемую точку. Иными

словами, остаточные аномалии силы тяжести, получившиеся после удаления топографических масс, без всякого изменения переносят на поверхность Н —

= Щ = сопв!. Если перенести аномалии с физической поверхности на поверх-

1 Е. М. О р л о в а . «Тр. ЦНИИГАиК». М., «Недра», 1965, вып. 157.

219'

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