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

shimbirev_b_p_teoriya_figury_zemli

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

В последнем интеграле взят знак минус по той причине, что внешняя нормаль п направлена на плоскости Оху в сторону отрицательных 2.

Поскольку потенциал точечной массы на поверхности сферы

т

,

йт

 

1

 

Я

'

получим

 

 

 

дТ

 

,

йт

дЯ

~~ ~

'

Если провести вспомогательную сферу Я = 1, то

6о = Я2(1(й,

где ш — часть поверхности вспомогательной сферы, соответствующая части сферы а (Я), лежащей ниже плоскости 2 = 0. Тогда будем иметь

п (Н)

о (й)

а

При Я -*- ьо, со -*• 2л, а интеграл по ограниченной части плоскости 2 = 0 переходит в интеграл по всей плоскости 2 = 0.

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

+ СО +оо

2

|

|

йх йу.

(XIV. 1

Если в точке С находится не точечная масса, а некоторое тело, то, сум-

мируя по всем элементам массы вт, тела и пользуясь свойством аддитивности потенциала, найдем

+ 0 0 4-СЮ

 

2я/Ма = |

^ ^ Л х й у ,

(Х1У.2'

— оо

- о о

 

где под Т следует понимать потенциал притяжения, создаваемый аномальное массой Ма. Поскольку, как мы условились, под аномалией силы тяжести пони-

мается вертикальная составляющая притяжения, то, очевидно, что

АД Т

& = ~дГ*

т. е. под интегралом в (XIV.2) стоит аномалия силы тяжести, создаваемая искомым аномальным телом массы Ма. Определим координаты центра инерция

тела (%0, г]0, ?„). Рассмотрим интеграл

+СО + О Э

СО - С О

Поскольку на плоскости Оху потенциал Т точечной массы и его производ-

ная по % симметричны относительно проекции точки С на плоскость, то

+со 4-<х>

390.

откуда

-Ьоэ +со 4-со -рсо

* 1

$

 

1 1

х^-йхйу.

— со —оо

 

—00 - с о

 

 

Учитывая (XIV.!), получим

+СО со

 

 

 

 

 

 

 

 

— со —со

 

 

Переходя от элемента массы йт к массе тела Ма,

находим

 

 

+ СО

+ С О

 

 

 

2 Л / М А | 0 = | | Х ^ Й Х Й У ,

( Х 1 У . З )

 

 

-оо -оо

 

 

так как на основании (11.29)

 

 

 

 

Аналогично получим 2я/Мат]0

= +^СО ^+СОУ ^ й х й у ,

(Х1У.4)

 

 

-00 -00

 

 

Глубина залегания центра массы

определяется по формуле Г. А. Гам-

бурцева, которую приводим без вывода,

 

 

 

+со +оо

 

 

 

6 Л ! М А 1 Й = \ ^ ( Т + 2 Т * ^ ) Й Х Д У .

( Х 1 У . 5 )

 

-со —оо

 

 

 

При практическом применении формул (Х1У.2), (Х1У.З), (XIV.4) и (Х1У.5) приходится ограничиваться интегрированием в некоторых конечных пределах. Входящие в формулу (Х1У.5) величины Т и д^Т/дт?, как правило, непосредственно не измеряют, а вычисляют через дТ/дг по интегральным формулам.

Применение столь сложных приемов определения при малой точности полученных результатов делается нецелесообразным. Для оценки глубины залегания аномального тела может быть применен метод, идея которого была высказана Ф. А. Слудским.

Так как направление вектора притяжения всегда проходит через притягивающее тело, то вполне возможна пространственная засечка положения тела. Правда, в общем случае направления векторов не пересекаются в одной точке, но этот метод дает возможность ориентировочно оценить положение притягивающего тела.

Для получения направления вектора притяжения следует вычислить проекции этого вектора на координатные оси, учитывая, что дТ/дг = Ад. дТ/дх и дТ/ду вычисляют по формулам для составляющих уклонений отвеса,

поскольку на основании (VI.7) имеет место зависимость

= — ?

дх ~

дТ

1 Г — ™

391.

Используя формулы для

плоскости

(IX.56) и (IX.57) и полагая в нп:

с1о = г (1г с1А, получим

оо 2 Я

 

дТ

 

5

^-соаАйгбА,

дх

о о

 

 

 

 

со 2Я

 

ду

1 Г Г М . * Ъ А Й Г Й А .

2д ^ ^

г

Для улучшения сходимости интегралов при малых г целесообразно Л.

уменьшить на постоянную величину Ад0, т. е. на значение аномалии в Т' : точке, для которой вычисляются дТ /дх и дТ /ду. Получаем известные формул-:

Б. В. Нумерова

дТ

дх

дТ_ ^ | —А§-0

ду - о о

Дальнейшее построение засечки можно осуществлять графически п.~: аналитически.

§ 88. ВТОРЫЕ ПРОИЗВОДНЫЕ ВОЗМУЩАЮЩЕГО ПОТЕНЦИАЛА

Эти производные, как и производные более высоких порядков от потекциала силы тяжести, играют важную роль при интерпретации, посколыг они обладают гораздо большей чувствительностью по отношению к притягивающим массам, нежели сами аномалии силы тяжести. Это объясняется те.", что чем выше порядок производной, тем сильнее ее зависимость от расстояния Так, например, если потенциал притягивающей точки У = / (т/г) есть функция, убывающая с расстоянием как 1/г, то первая производная потенциала —

сила притяжения

Р =

йУ/Лг =

—/ (т,/г2), — убывает уже как 1/г2,

а вторе.-:

производная дРУ/йг2 =

2/ (т/г3),

как 13

и т. д. На рис. 77 приведены г; -

фики изменения

аномалии силы тяжести

А§, ее вертикального

градиент

д Ад/дг и второй

вертикальной производной аномалии силы тяжести дг А§/о:-

от двух круговых цилиндров бесконечного протяжения. В приведенном примег

кривая Д^ не показывает наличия двух источников гравитационных аномалп;: оно слабо проявляется и на поведении кривой д Ад/дг и только кривая д2

ясно обнаруживает наличие двух аномальных масс.

Потенциальная функция силы тяжести УУ имеет шесть вторых производ-

ных, которые принято обозначать как

дх>±

Ухх'

ду2

д22

 

дхду

 

дхдг

А'2' ду д2

У2'

Эти производные потенциала имеют размерность 1/с2. Единицей изме; ния является этвеш, составляющий 10"9 (1/с2).

392

Производные д2УУ!дх2, д21У/ду2 и д2Ц>г2/дх ду, как следует из формулы

(1У.15), определяют кривизну уровенной поверхности силы тяжести, производная д2УУ/дг2 является вертикальным градиентом силы тяжести а д2Ш/дх дг и д2Ш/ду дг — горизонтальными градиентами силы тяжести.

Величины д2№/дх ду, д2Ц//дх дг, д2\У/ду дг и д2Ш/ду2 — д2Ш/дх2 =

получают из наблюдений с вариометрами, однако для целей гравиметрической разведки необходимо получить аномальные значения этих производных.

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

Приведем нормальные значения вторых производных для уровенного эллипсоида

° х г - дх дг дх М

 

 

 

ау

 

 

 

иУ*

дудг

ду

'

 

 

VХУ = дх^ду

= 0'

 

 

тт

дЮ

дт

Уее2

,

В,

 

= —^7,

=

 

СОЗ2

где М — радиус кривизны меридианного эллипса; а — большая полуось эллипсоида; е — его эксцентриситет. Подставив значения постоянных, получим

в этвешах

Vхг = 8,11 з т 2В, 11уг = 0, ЬТху = 0, ПА = 10,24 соз2 В.

393.

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

Влияние рельефа учитывается при помощи формул, выведенных различными авторами (Этвеш, Швейдар, Никифоров, и др.), а также графическими методами.

 

Аномальное значение второй производной, например

величины IVху

получается следующим образом:

 

где

(\Уху)а — аномальное значение второй производной;

ху — наблюден-

ное;

Лху — нормальное значение; {№ху)т — поправка за рельеф.

 

Однако в настоящее время вариометры в гравитационной разведке приме-

няются редко и их все в большей степени вытесняют высокоточные гравиметры.

Используя

результаты

наблюдений с гравиметрами,

можно

получить

горизонтальные

градиенты аномалии силы тяжести: д

Ад/дх

и д А§/ду и вто-

 

 

Т а б л и ц а

43

рую вертикальную производную ано-

 

 

малии силы тяжести д2 Ад/дг2.

 

 

 

 

 

 

Учитывая

связь между аномалией

в

ил

иху

Пхг

иУг

силы тяжести Д^ и возмущающим по-

 

 

 

 

 

 

тенциалом,

можно

написать

 

0

10,2

0

0,0

 

0

дАе

д

 

 

 

32т

 

10

9,9

0

2,8

 

0

дх

дх

 

 

 

 

 

20

9,0

0

5,2

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

30

7,7

0

7,7

 

0

дАё

д

 

 

 

 

 

40

6,0

0

8,0

 

0

 

 

 

 

 

 

ду

ду

 

 

 

 

 

50

4,2

0

8,0

 

0

 

 

 

 

 

60

2,6

0

7,0

 

0

 

92

 

 

 

 

 

70

1,2

0

5,2

 

0

дг2

(

дТ

\

_ 83т

,

80

0,3

0

2,8

 

0

дг2

V

дг

 

а2з

"

90

0

0

0

 

0

По предложению

М. У.

Сагитова

 

 

 

 

 

 

 

 

 

 

 

 

[11] эти величины можно получить через

аномалии

силы

тяжести.

Возьмем

прямоугольную

систему

координат хОу

с началом в точке, где предполагается вычислить вторую вертикальную производную Тггг; построим сеть квадратов со стороной 21. Координаты вершин

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

 

( 2 '

 

2 ) '

(

2 '

2 ) ' (

2 '

2 ) '

(

2 '

2~)

будем иметь

 

 

 

 

 

 

 

 

 

 

 

 

: (

2 '

2 )

 

1_ I \ _ А? (I, 1) + Ае

(I, 0)—Л^ (0, Р — Ац (0, 0)

 

 

 

 

21

 

 

 

 

 

 

(

2 '

2 )

 

_ А у ( 0 . Г) + А * ( 0 , 0 ) — А ^ ( — I, 0 ) — А § (—I, I)

 

Т х г

 

 

 

21

 

 

 

 

(Х1У.6)

т

(1_

 

I \

 

А$ (I, 0) + Ае (I, -0—А*

 

(0, 0)—А* (0, -I)

 

 

 

 

 

Х 2 \ 2 '

 

2 )

~

 

 

21

 

 

 

 

 

« Ч

2 '

 

2 )

 

А?(0,

0 ) + А ^ ( 0 ,

-1)-Аё(~1,

0 ) — Д ? ( — г ,

-I)

 

 

 

 

21

 

 

 

 

394.

Аналогично

получим

выражения

для

горизонтальных

градиентов Туг

в тех же точках:

 

 

 

 

 

 

 

 

тУг (г

2г

А У ( * ,

1) + А§

(0,

2/!

0 ) - Д У ( 0 , 0)

 

 

Туг(

 

г

г •

д*(0,

г)+д*(-*,

г>—а^г(о, о ) - д г ( - г ,

0)

 

Т '

 

 

 

21

 

(XIV. 7)

 

г

 

 

 

 

 

 

1 '

 

о)+д?(0, 0)—

-о—Ду(о,

- о

 

 

 

 

 

2

, —

 

 

 

21

 

 

 

г

 

1

А * ( 0 ,

 

<—Г,

0 ) — Д ^ ( О , - Г ) — Д у ( - Г . - г )

 

 

"2~*

2

 

 

 

 

 

 

 

г - ^

( П )

(~2>г)

М до

т )

М

0,-1)

Рис. 78

Имея значения вторых производных Тхг и Туг в выбранных точках, полу-

чим выражения для вторых производных аномалии силы тяжести в центре палетки (0, 0)

 

 

Тххг(

0,0) =

 

 

тХг ( т ,

+

( 4 - ' —т)~Тхг

( ~ Т ' т)~Тхг

( ~ ~ Т '

~ т )

 

 

 

22

 

 

 

 

 

0) =

 

 

—• т )

+ Г г / 2 (

2"' т ) ~ Т у г { т > ~ т ) ~ т

у г { ~ т >

~ т )

>

(Х1У.8)

Известно, что потенциал притяжения вне притягивающих масс удовлетворяет уравнению Лапласа

ТXX Туу-С Т22 = 0.

Продифференцируем все члены последнего уравнения по 2. Выразим вторую вертикальную производную Тггг аномалии силы тяжести через производные Тххг и Тууг в точке (0, 0)

7 ^ ( 0 , 0) = -Гх х г (0, 0 ) - Г д а ( 0 , 0).

395.

Принимая во внимание соотношения (Х1У.6), (XIV.7) и (XIV-8), в нос.-

нем равенстве выразим Тххг и Тууг через аномалии силы тяжести

Ад.

Т- (0

0 ) — О — А * (Л -0—Ау(-г. 1)-&В(-1, -I)

(XIV-

Таким образом, вторая вертикальная производная Тггг аномалии сг:

тяжести в некоторой точке (начале координат) с точностью до постоянна" множителя (1/2) равна учетверенному значению аномалии силы тяжести в с- точке минус аномалии силы тяжести в вершинах квадратов.

§ 89. РЕШЕНИЕ ПРЯМОЙ ЗАДАЧИ ГРАВИМЕТРИЧЕСКОЙ РАЗВЕДКИ

Как уже говорилось выше, нахождение путем вычислений значений ну

малий силы тяжести или значений вторых производных возмущающего пот» - циала от находящегося в недрах Земли тела любой формы и с любым распре.- лением плотности в нем для точек, расположенных на земной поверхност носит название прямой задачи гравиметрической разведки.

При этом поверхность Земли принимается за плоскость, а оси коордиг принято располагать следующим образом: оси х и у помещаются в плоско!:

горизонта, а ось г направляется вертикально вниз.

Прямая задача, в противоположность обратной, всегда имеет однозначн' решение, получение которого не связано с какими-либо принципиальны: трудностями. Однако математические вычисления, которые необходимо п; этом производить, отличаются большой сложностью. Сложность эта связг. с тем, что всякое аномальное тело является трехмерным и потому при оп; делении параметров гравитационного поля приходится выполнять трехкратнинтегрирование, что и приводит при сложной форме тела к чрезвычайно громо; ким выражениям. Вследствие этого в гравиметрической разведке нашли шп' кое применение так называемые двухмерные тела, определение гравиметриских характеристик которых значительно проще, чем трехмерных.

Под двухмерным телом понимается тело бесконечного простирания в однг направлении. Форма двухмерных тел целиком определяется формой их пои речного сечения. Двухмерное тело можно определить как цилиндрическ полученное при движении обра&ующей по контуру поперечного сечек: В том случае, если контур имеет форму многоугольника, мы получим прпг' тическое двухмерное тело. Горизонтальный однородный круглый цилинд притяжение которого равно притяжению эквивалентной массы, расположена по оси цилиндра в форме материальной линии, является простейшим случг - двухмерных цилиндрических тел.

Погрешности, вызванные заменой трехмерных тел двухмерными, обыч> невелики и для вытянутых по простиранию геологических тел, длина котор: в 5 раз и более превосходит их поперечник, можно считать, что гравитационгполе над серединой их такое же, как если бы эти тела простирались в бесконе- ность.

Только для очень ограниченного числа тел, имеющих простейшую форм~ удается получить удобные для практического использования формулы. Конечк подобные тела в природе встречаются очень редко. Обычно формулами для т-. простой геометрической формы пользуются лишь для первоначальной прпбгг зительной оценки основных параметров тела, после чего переходят к постеш ному подбору таких параметров, которые наилучшим образом соответствов;

396.

бы данным, полученным из наблюдений. Этот подбор ведется при помощи методов численного интегрирования.

Для расчета кривых аномалий силы тяжести, обусловленных простейшими по своей форме телами, применяются формулы:

а) цилиндрическое кольцо

 

 

 

 

А^ - 2л/ Аб [ / 1 \ + Я* -

\ / +

4 - 1к+1 - 1к],

что следует из (VI 1.12) при п =

1;

 

 

б) цилиндр

 

 

 

 

Аг = 2 * / А б Я ( 1 - 4 - ) ;

в) горизонтальный бесконечный слой

 

 

Д^ = 2я/Лб#;

 

г) конус

 

 

 

 

А^ = 2я/АбЯ(1 — з т I),

где угол I определяется как г =

агс

Н/а.

 

Во всех этих формулах под

Аб понимается аномальная плотность Дб =

= б2 — б-р

где б2 — плотность

данного тела;

бх — плотность окружающих

данное тело

пород.

 

 

 

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

производных.

 

Потенциал материальной точки с массой йт

в соответствии с (1.15) будет

7 =

(ХГОО)

где г определяется по формуле (1.5). Первые производные потенциала в соот-

ветствии

(1.16)

будут

 

 

 

 

 

 

 

 

 

 

 

 

 

 

• ^

=

 

 

 

 

=

 

 

 

 

=

 

 

 

Йт. (XIV. 11)

 

дх

 

1

г3

 

 

с*у

 

1 гЗ

 

 

да

 

1

/3

 

\

>

При

вычислении

вторых производных

потенциала

будем

учитывать,

что

 

 

 

 

дг

 

х—с,

 

дг

у

д г

 

2—5

 

 

 

 

 

 

 

 

 

дх

 

г

'

ду

г

 

'

Зг

г

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

д^У

д

(

дУ

\

 

|

дг

,

,

Дт

„ ,

х

,

,

 

Ат

 

 

 

 

I

 

х — I

(1т

(х—|)г

 

 

гя

 

дх%

дх

\

дх

}

'

г4

дх

 

-

г з

'

 

 

 

 

1

 

92V

д

/ дУ

\

п ,

у —Г1

дг

,

,

(1т

„ ,

(и—Г))2

,

 

,

(

 

ду*

ду

\

ду

)

'

г*

ду

 

'

г»

 

гб

 

 

 

'

гз

 

Получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

^

= •

=

3/

 

 

 

 

Л*.

 

 

 

(Х1У.12)

3 9 7

Далее

=

(XIV,1:

( X I V , -

Поскольку шар притягивает по тому же закону, что и материальная точка, воспользуемся только что выведенными формулами для получения соответствующих величин для однородного шара, положив Лт = 4/3 (яАбД3) = М

и понимая, как и прежде, под Дб аномальную плотность и под М а аномальнук массу.

О

Р(х,0,0)

Будем считать, что центр шара находится на оси 2 на некоторой глубине Вычислим возмущающий потенциал, аномалию силы тяжести и вторые производные возмущающего потенциала для точки М, расположенной на осп +

Полагая в формулах (ХГУ.10), (ХГУ.11), (Х1У.12), (Х1У.13), (ХГУ.14) и (XIV.!.:

^ = т] =

г/ = 2 = 0, и считая,

что

г

=

]/ж2 +

понимая под г

расстояние

от центра

шара до точки М,

получим

 

 

 

 

 

 

 

 

,

Ма

'

 

 

(XIV. 1-

 

 

 

~Т

 

Г

 

 

 

 

 

__ зт

 

 

 

 

 

(XIV. 17

 

 

 

 

 

 

 

 

 

 

 

дхду

=

Тху

= 0,

 

 

(XIV. 181

 

 

д2 Т

 

 

 

 

 

 

(XIV. V?

 

 

ду дг.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дх дг

х

 

 

 

 

 

 

(Х1У.2С,

 

дчт

д*Т

 

 

 

 

 

 

(Х1У.21

 

ду2

дх*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

На рис. 79 представлены графики кривых Д^ и Тхг.

Теперь рассмотрим

различные случаи для цилиндрических тел.

 

 

 

1. Г о р и з о н т а л ь н а я

в е щ е с т в е н н а я

п р я м а я

(горизон-

тальный круговой цилиндр).

Пусть

на

рис. 80

изображена вещественная

398.

горизонтальная

линия,

лежащая

в

 

плоскости уг на глубине

Обозначим

линейную плотность через Я, так что с1т =

Я

Расстояние от элементарной

массы йт до точки Р

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

г = Т / ^ + л' + ^о2-

 

 

 

Используя

формулу

(XIV.И),

получим

 

 

 

 

^

д г - / ?

 

 

^ ^

 

 

 

 

/я:

?'

^

 

Обозначим постоянную величину х2

+

^о через с2. Представим последний:

интеграл в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Г

дз

'

 

 

 

 

 

где

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

К =

1 / У + с 2 .

 

 

 

 

 

 

 

 

 

 

 

Отсюда получаем

 

=

-г|с?т]/Т/гг/2 —|- с2.

Поскольку

Я' =

т1/1/т]2 + с2 =

= "П/Д, имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ЯЯ' =

г|,

 

 

 

 

(XIV.22),

 

 

 

 

 

 

=

Л' йт).

 

 

 

(XIV.23).

Дифференцируем (Х1У.22)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ЯйЯ'

+ Я' йЯ = <1цч

 

 

Отсюда с учетом соотношения (XIV.23) будем иметь

 

 

 

 

 

 

 

<1П'

 

йт)

 

 

 

 

На основании

равенства

1 - й ' 2

 

Я

 

 

 

 

 

 

 

 

 

 

 

 

 

 

учитывая (Х1У.22),

получим

Я2 — г\2 = с2,

 

 

 

1

 

1—Я'2

 

 

 

 

 

 

 

 

 

 

 

 

 

(XIV.25).

 

 

 

 

 

Д2

 

С2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом,, подынтегральное выражение с учетом (XIV.24) и (XIV.2.5)?

можно представить как

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ЛИ^ЛI

д

' дг в

С2

 

 

 

 

 

 

л»

_

 

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

с ^

 

1

и н '

~

я

'

 

 

п

 

 

 

 

О Я3

Са

 

 

 

С2 -

 

с2|Л^)172 '

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

399 >

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