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

МВ до вик контрольної р (1)

.pdf
Скачиваний:
4
Добавлен:
28.06.2022
Размер:
186.37 Кб
Скачать

МІНІСТЕРСТВО ОСВІТИ І НАУКИ, МОЛОДІ ТА СПОРТУ УКРАЇНИ НАЦІОНАЛЬНИЙ УНІВЕРСИТЕТ “ЛЬВІВСЬКА ПОЛІТЕХНІКА

ОСНОВИ ГРАВІМЕТРІЇЇ

МЕТОДИЧНІ ВКАЗІВКИ

до виконання контрольної роботи для студентів 3-го курсу базового напряму

“Геодезія, картографія та землеустрій” заочної форми навчання

Затверджено на засіданні кафедри вищої геодезії

та астрономії.

Протокол № 8-10/11 від 11.02.2011 р.

Львів – 2011

Основи гравіметрії: Методичні вказівки до виконання контрольної роботи для студентів 3-го курсу базового напряму “Геодезія, картографія та землеустрій” заочної форми навчання / Укл. П.Д. Двуліт. – Львів: Видавництво Львівської політехніки, 2011. – 20 с.

Укладач

Двуліт П. Д., д-р техн. наук, проф.

Відповідальний за випуск Заблоцький Ф.Д., д-р техн. наук, проф.

Рецензент

Савчук С. Г., д-р техн. наук, проф.

2

1. Загальні зауваження

Відповідно до навчального плану дисципліну „Основи гравіметрії” вивчають на третьому курсі. Студіюючи дану дисципліну, студенти заочної форми навчання повинні знати основні характеристики гравітаційного поля Землі, вміти виконувати вимірювання прискорення вільного падіння сучасними маятниковими приладами та статичними гравіметрами під час створення опорних гравіметричних мереж, регіональних і детальних гравіметричних знімань та вміти опрацьовувати результати вимірювань, включаючи складання гравіметричних карт.

Самостійне вивчення дисципліни завершується виконанням контрольної роботи. Під час лабораторно-екзаменаційної сесії студенти заслуховують установчу і оглядові цикли лекцій, що триває 6 годин, і виконують 2 лабораторні роботи.

2. Програма основних розділів дисципліни

Вступ. Зміст предмету „Основи гравіметрії”. Гравітаційне поле Землі та їх характеристики.

Абсолютні вимірювання прискорення вільного падіння.

Маятниковий метод абсолютних визначень прискорення вільного падіння. Балістичний метод абсолютних визначень прискорення вільного падіння. Прилади для абсолютних вимірювань. Точність вимірювань. Результати абсолютних визначень.

Відносні вимірювання прискорення вільного падіння.

Відносні маятникові вимірювання та їх результати. Статичний метод відносних вимірювань прискорення вільного падіння. Прилади для вимірювань. Точність вимірювань. Методика польових вимірювань із сучасними маятниковими приладами і статичними гравіметрами та опрацювання результатів вимірювань.

Вимірювання других похідних потенціалу прискорення вільного падіння.

Гравітаційні варіометри і градієнтометри. Методика роботи і опрацювання результатів спостережень. Геодезичне використання вимірювань градієнтів прискорення вільного падіння і різниць кривин рівневої поверхні.

3

Опорні мережі та гравіметричні знімання.

Світова гравіметрична мережа. Національна і польова гравіметричні мережі. Топографо-геодезичне забезпечення гравіметричного знімання. Аномалії прискорення вільного падіння та їх види. Методика складання гравіметричних карт.

3. Методичні вказівки до виконання контрольної роботи

Контрольна робота складається із двох завдань і включає основні теми розділів програми дисципліни. Виконання контрольної роботи необхідно починати після вивчення тем основних розділів програми. Для виконання контрольної роботи студенту необхідно вибрати із вихідних даних, які поміщені в даних методичних вказівках свій номер варіанту, використовуючи дві останні цифри шифру.

Завдання № 1

Метою завдання є ознайомлення із методами обчислень аномалій прискорення вільного падіння за спостереженнями в різних умовах, а також довідатися про порядок величин цих аномалій на фізичній поверхні Землі. Для цього необхідно виконати обчислення аномалій у вільному повітрі й аномалій Буге для точок, розміщених:

1)на поверхні суходолу Землі у рівнинному районі;

2)на поверхні суходолу Землі у гірському районі;

3)на поверхні моря;

4)на дні моря;

5)над Землею у повітрі;

6)у свердловині або в шахті.

Під аномалією прискорення вільного падіння розуміють різницю між виміряним значенням прискорення вільного падіння і його нормальним або теоретичним. Виміряне значення g отримують із абсолютних або відносних спостережень, а теоретичне γ розраховують за формулами нормального розподілу прискорення вільного падіння геодезичної референцної системи. Якщо виміряне і нормальне значення прискорення вільного падіння віднесені до однієї точки, то таку аномалію називають „чистою”. У разі, коли ці значення належать до різних точок простору, отримують “змішану” аномалію прискорення вільного падіння. При обчисленнях аномалій вводять різні поправки у виміряне або теоретичне значення прискорення вільного падіння.

4

Залежно від того, які вводять поправки (редукції), розрізняють різні види аномалій. Найбільше практичне застосування у гравіметрії та фізичній геодезії мають аномалії у вільному повітрі й аномалії Буге.

Спочатку приведемо необхідні формули для обчислень „змішаних” аномалій у вільному повітрі і аномалій Буге залежно від розташування гравіметричного пункту. Очевидно, для таких пунктів спостереження необхідно обчислити аномалії, але, насамперед, потрібно їх знати на фізичній поверхні Землі. Зауважимо, що у разі обчислення „змішаних” аномалій, замість геодезичних використовують нормальні висоти.

4. Обчислення аномалій у вільному повітрі

Обчислюючи аномалії, необхідно враховувати розміщення пункту спостереження. Так, для точок 1 і 2, розміщених на суходолі поверхні Землі, аномалію у вільному повітрі обчислюють за формулою

(g - γ )в.п. = g - (γ 0 + Dγ ),

де g – виміряне значення прискорення вільного падіння,

γ0 – нормальне значення прискорення вільного падіння на поверхні рівневого нормального еліпсоїда,

γ – редукція у вільному повітрі.

Значення γ 0 обчислюють за формулою 1980 року

γ0 = 978032,7(1+ 0,0053024sin2 B - 0,0000058sin2 2B)

або вибирають із таблиць. Величину γ , яку називають поправкою за висоту в нормальне значення прискорення вільного падіння або редукцією у вільному повітрі, обчислюють за формулою

Dγ = -0,30855(1+ 0,00071cos2B) H + 0,0723× H 2 ×10−6 ,

де B – геодезична широта пункту спостереження, H – геодезична висота.

Зауважимо, що якщо геодезичну висоту виражати у метрах, то тоді поправка γ буде виражатися в мілігалах.

Для точки 3, яка розміщена на поверхні моря, поправка γ дорівнює нулю, оскільки в цьому разі нормальна висота H γ точки спостереження дорівнює нулю. Точка 4 розташована нижче рівня моря і її нормальна висота буде від’ємною. Тут під висотою H у формулі розуміють глибину моря. Тоді аномалія у вільному повітрі для точки 3 дорівнюватиме

5

 

(g − γ )в.п. = g − γ0 ,

а для точки 4 аномалія на поверхні моря

 

(g − γ )в.п. = g (γ 0 + γ ) + 1g ,

 

D1g = 0.08387δм.в.,h = 0.0864h мГал ,

де D1g – притягання шару морської води

h

глибина моря в метрах,

 

δм.в.

– густина морської води, δм.в. = 1.03 г / см3.

Аномалію у вільному повітрі в точці 5 обчислюють за формулою

аналогічною як для точок 1 і 2, але у разі обчислення поправки γ замість H

необхідно

використати значення

H + h, бо висота точки спостереження

складається із нормальної висоти

H γ і висоти h літака (гелікоптера) над

фізичною поверхнею Землі.

 

Точка 6 розміщена у свердловині, а її висота над відліковою поверхнею

дорівнює

H γ h , а h – глибина свердловини. Щоб обчислити аномалію у

вільному повітрі на фізичній поверхні Землі, а не в точці спостереження, необхідно врахувати притягання 1 g шару товщиною h вище точки спостереження

D1g = 0,08387δ × h ,

де δ – густина гірських порід, h – глибина свердловини.

Тоді аномалію у вільному повітрі для точки 6 обчислюють за формулою

(g − γ )в.п. = g (γ 0 + γ ) + 1g .

5. Обчислення аномалій Буге

Для обчислення аномалій Буге необхідно враховувати місцерозташування точки спостереження: суходіл чи море. Для суходолу (точки 1, 2, 5 і 6) аномалію Буге обчислюють за формулою

(g − γ )Б = (g − γ )в.п. + 2 g ,

де (g − γ )в.п. – обчислене значення аномалії у вільному повітрі для відпо-

відної точки,

D2g – притягання топографічних мас, які розміщені між рівнем моря і фізичною поверхнею Землі, яку приймають за площину, або поправка за вплив проміжного шару.

Поправку D2g обчислюють за формулою

D2g = -0.0419δ × H ,

або вибирають із таблиць.

6

У разі, коли точка розміщена на поверхні моря (точки 3 і 4), під час обчислення аномалії Буге враховують дефект мас, обумовлений тим, що густина морської води δм.в. є меншою, ніж середня густина земної кори δ .

Притягання 2g мас, які доповнюють густину морської води до середньої густини порід, обчислюють за формулою

2g = 0.0419(2.67 −1.03)h = 0.0687h .

Ця поправка завжди додатна. Її можна обчислити за формулою або вибрати із таблиць. Тоді аномалію Буге для точок 3 і 4 обчислюють за формулою

(g − γ )Б = (g − γ )в.п. + 2g .

Числовий приклад обчислення аномалій у вільному повітрі й аномалій Буге

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B º '

Hγ, м

 

 

h, м

 

Hγ +h,

 

g, мГал

 

 

γ0, мГал

точок

 

 

 

 

 

 

 

м

 

 

 

 

 

 

 

 

1

 

52º13'

5

 

0

 

5

 

 

981274.8

 

981266.6

2

 

36º48'

384

 

0

 

384

 

 

979851.0

 

979888.3

3

 

4º 22'

0

 

3820

 

0

 

 

978072.8

 

978062.6

4

 

25º45'

0

 

125

 

125

 

 

979069.3

 

979008.0

5

 

67º17'

143

 

500

 

643

 

 

982192.7

 

982442.3

6

 

48º50'

125

 

40

 

85

 

 

980924.7

 

980965.9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Δγ,

 

γ ,

 

1g, мГал

 

(g-γ)в.п.,

2g, мГал

 

(g-γ)Б,

точок

 

мГал

 

мГал

 

 

 

мГал

 

мГал

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1.5

 

981265.1

 

 

 

 

 

 

+9.7

 

-0.6

 

+9.1

2

 

118.5

 

979769.8

 

 

 

 

 

 

+81.2

 

-43.0

 

+38.2

3

 

0

 

 

 

 

 

 

 

 

 

+10.2

+262.4

 

+272.6

4

 

38.6

 

979046.6

 

10.8

 

 

 

+33.5

 

+8.6

 

+42.1

5

 

198.4

 

982243.9

 

 

 

 

 

 

-51.2

 

-16.0

 

-67.2

6

 

26.2

 

980939.7

 

9.0

 

 

 

-6.0

 

-14.0

 

-20.0

Завдання № 2

Визначити коефіцієнти формули розподілу нормального значення прискорення вільного падіння і стиснення рівневої поверхні.

Метою завдання є практичне знайомство з методом визначення стиснення рівневої поверхні за виміряними значеннями прискорення вільного падіння.

7

в.п.

Відомо, що зміну прискорення вільного падіння на поверхні рівневого еліпсоїда можна описати виразом

γ0 = γ e (1+ β sin2 B - β1sin2 2B),

вякому коефіцієнти γe , β , β1 називають фундаментальними геодезичними сталими. Для їх визначення необхідно мати значення γ0 . Крім цього, коефіцієнт β1 належить до сталих вищих порядків і будемо вважати його таким, що дорівнює 0,000006. Отже, для визначення коефіцієнтів γe і β

необхідно мати значення γ0 у двох точках. Але γ0 неможливо отримати із результатів вимірювань, тоді для знаходження γe і β використовують вимірювання істинного прискорення вільного падіння в точках фізичної поверхні Землі. Нормальне значення прискорення вільного падіння γ 0 можна виразити через результати вимірювань g із співвідношення

γ0 = g - Dγ + D1g - (g - γ )

Тоді можна записати, що

γe + γeβ ×sin2 B - (g - Dγ + D1g + γ eβ1sin2 2B) = -(g - γ )в.п.

Якщо ввести позначення

x = γ e − 978000 , y = γ eβ ,

b= sin2 B ,

υ= -(g - γ )в.п. ,

l = -(g - Dγ + D1g + γeβ1 sin2 2B)+ 978000 ,

то отримаємо рівняння похибок

x + by + l =υ .

Кількість таких рівнянь буде дорівнювати кількості точок із виміряними значеннями прискорення вільного падіння g . Невідомі x і y визначають з розв’язку системи нормальних рівнянь за умовою, що [υυ]= min .

nx + [b] y + [l] = 0, [b] x + [bb] y + [bl] = 0 ,

де n – кількість вимірювань.

Розв’язуючи систему нормальних рівнянь, знаходимо невідомі коефіцієнти

γe = x + 978000 ,

β= g e .

8

Використовуючи теорему Клеро, визначаємо стиснення α

α + β = 52 q ,

а звідси

α = 52 q - β .

Параметр q дорівнює відношенню відцентрової сили ω2a на екваторі до прискорення вільного падіння γe на екваторі. Якщо прийняти a = 6378.2 км , ω2 = 53,175×10−10 ( рад / с)2 і γe = 978000 МГал , то тоді q = 0.003468 .

Для розв’язування нормальних рівнянь складемо визначник D системи

D =

 

 

n

 

[b]

 

 

= n[bb]-[b]·[b] .

 

 

 

 

 

 

 

 

 

 

 

[b] [bb]

 

 

 

 

 

 

 

 

 

Невідомі x і y , а також їх вагові коефіцієнти знаходимо за формулами

Dx =

 

-[l]

 

[b]

 

 

= -[l][bb]+ [b]·[bb] ,

 

 

 

 

 

 

 

-[bb] [bb]

 

 

 

 

 

 

 

 

 

D

y

=

 

 

n

 

-[l]

 

= -n

[bl]+ [b]·[l] ,

 

 

 

 

 

 

 

 

[b] -[bl]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x =

 

D

x

, y =

 

Dy

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

 

 

 

 

 

 

 

 

 

 

 

D

 

 

 

 

 

 

 

 

 

 

 

β =

 

y

, α = 5 q - β .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γ e

 

 

 

2

 

 

 

 

 

Точність результатів оцінюють за відомими правилами математичного

опрацювання

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

μ =

 

 

[υυ]

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n - 2

 

 

 

 

 

 

 

D

Px

=

 

 

1

 

 

 

[b]

 

= [bb] ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

[bb]

 

 

 

 

1 = DPx ,

Px D

mγe = μPx ,

DPy =

n

0

= n,

[b]

1

9

1 = DPy ,

Py D mγeβ = μPy ,

mβ = mγeβ ;

γ e m= mαβ2 .

Числовий приклад обчислення коефіцієнтів формули розподілу нормального значення прискорення сили ваги

№ точок

B

g-

 

-Δγ+

 

 

1g,

 

γeβ1sin22B,

 

l,

γeβ∙b,

V,

b

978000,

 

мГал

 

 

мГал

 

мГал

мГал

мГал

 

 

мГал

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

52˚13′

3275

 

 

 

 

 

 

+2

 

 

 

+6

 

 

 

 

 

-3283

3189

-13

0.62463

2

36˚48′

1851

 

 

 

 

+118

 

 

 

+6

 

 

 

 

 

-1975

1832

-62

0.35883

3

4˚ 22′

73

 

 

 

 

 

 

0

 

 

 

0

 

 

 

 

 

-73

80

-38

0.00580

4

25˚45′

1069

 

 

 

 

 

 

-28

 

 

 

+4

 

 

 

 

 

-1045

964

0

0.18874

5

67˚17′

4193

 

 

 

 

+198

 

 

 

+4

 

 

 

 

 

-4395

4344

+30

0.85087

6

48˚50′

2925

 

 

 

 

 

 

+35

 

 

 

+7

 

 

 

 

 

-2967

2894

+8

0.56670

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-13738

 

+1

2.59557

 

 

 

D =

 

 

 

 

6

 

 

 

2.59557

 

 

= 2.86105,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.59557

1.59967

 

 

 

 

 

 

 

 

 

 

 

Dx

=

 

 

13738

2.59557

 

= 230.581,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8378

1.59967

 

 

 

 

 

 

 

 

Dy =

 

 

6

 

 

 

13738

 

= 14610.059 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.59557

8378

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x =

 

 

Dx

= 80.6

МГал ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y =

Dy

 

= 5106

МГал ,

 

 

 

 

 

 

 

 

 

 

 

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10