Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
fulltext-2011-Kuzmin_Gravimetria_2011.doc
Скачиваний:
0
Добавлен:
01.07.2025
Размер:
9.51 Mб
Скачать
  • графу Σ 3 .

    1. Суммируем значения сумм по секторам (по всей строке Σ 0 ) и значения

    сумм по зонам (по столбцу Σ 3 ). Суммы по строке и по столбцу должны быть равными:

    0  3 .

    (5.96)

    После этого контроля приступаем к обработке результатов.

    22. Выполняем действия, описанные в пунктах 4–9. Следует иметь ввиду, что в этой схеме каждой группе секторов соответствует свой коэффициент C i .

    В группе секторов 1, 7, 9, 15:

    ξ 

    cos 22 , 5

    0 , 0046

    

    C

    1

     0 , 005 

    .

    η

    sin 22 , 5

    0 , 0019 

    В группе секторов 2, 6, 10, 14:

    ξ 

    cos 45

    0 , 0035

    

    C

    2

     0 , 005 

    .

    η

    sin 45

    0 , 0035 

    В группе секторов 3, 5, 11, 13:

    ξ 

     0 , 005

    

    cos 67 , 5

    0 , 0019

    

    C

    3

    0 , 0046

    .

    η

    sin 67 , 5

    

    В группе секторов 16, 8 ,4, 12:

    ξ 

     0 , 005

    

    cos 360

    0 , 0050 

    C

    4

    0 , 0050 

    .

    η

    sin 90

    Вычисленные по

    группам секторов значения ξ i и ηi суммируем и

    записываем в строки ξ 5100 и η5100 :

    ξ 

    100

     ξ 

     ξ 

     ξ 

     ξ 

    5

    1

    2

    3

    4

    (5.97)

    η

     η

     η

     η

    η .

    100

    5

    1

    2

    3

    4

    23. Значения

    ξ 5100 и η5100 записываем в схему 5 (прил. 12), строка 5–

    100 км.

    24. Вычисляем значения ζ i (до 0,01 м). Сумму g i для каждой зоны умножаем на коэффициенты Ki 105 .

    1. Суммируем значения ζ i по всем строкам и записываем результат в строку Σ 4 (до 0,1 м).

    2. Вычисленное значение ζ 5100 записываем в схему № 6 (прил. 12),

    строка 5–100.

    27. Вычисляем суммарные значения составляющих ξ,η и ζ:

    ξ 

    100

     ξ 

     ξ 

    0

    05

    5100

    η

     η

    5

     η

    (5.98)

    0100

    0

    5100

    ζ 0

    100

     ζ 0

    ζ 05

    ζ 5100 (М).

    Требования к оформлению и защите лабораторной работы

    1. Работа выполняется на листах А4 с приложением результатов вычислений на бланках (прил. 11).

    1. Работа должна иметь титульный лист, на котором указывается наименование учебного заведения, номер лабораторной работы и ее название, номер группы и фамилия студента (на листе слева). Справа пишется фамилия преподавателя, его должность. Внизу страницы – год.

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

    1. Далее прилагаются вычислительные схемы с подписью вычислителя и

    датой.

    1. Записываются контрольные вопросы.

    КОНТРОЛЬНЫЕ ВОПРОСЫ

    1. Дать определение уклонения отвеса.

    1. Назвать виды уклонений отвеса и дать их определения.

    2. Величины уклонений отвеса и их размерность.

    3. Где и каким образом используются уклонения отвесной линии.

    1. Что означает высота квазигеоида, как она определяется, пределы ее изменения на земном шаре, размерность.

    2. Для чего нужно знать высоту квазигеоида?

    ЛАБОРАТОРНАЯ РАБОТА № 8

    ОПРЕДЕЛЕНИЕ ФУНДАМЕНТАЛЬНЫХ ГЕОДЕЗИЧЕСКИХ ПОСТОЯННЫХ ПО ГРАВИМЕТРИЧЕСКИМ ДАННЫМ

    ЦЕЛЬ РАБОТЫ: освоить метод определения сжатия уровенного эллипсоида по измерениям силы тяжести.

    МЕТОДИЧЕСКИЕ РЕКОМЕНДАЦИИ: к выполнению задания следует приступить после изучения работ [28, гл. VI, § 47, 48]. Общие теоретические сведения представлены в лабораторной работе № 5 данного издания. Исходные данные для выполнения лабораторной работы № 8 приведены в прил. 7.

    СОДЕРЖАНИЕ РАБОТЫ

    Определить коэффициенты формулы распределения нормальной силы тяжести.

    Вычислить значение сжатия уровенной поверхности α двумя методами.

    Динамический коэффициент

    J 2 0 , 00108 N ( N

    две последние цифры

    шифра).

    ПОРЯДОК ВЫПОЛНЕНИЯ РАБОТЫ

    1. Значение силы тяжести γ 0 на поверхности

    уровенного эллипсоида

    вычисляется по формуле:

    γ 0  γ e 1  βsin 2 B  β

    1 sin 2 2 B 14 ,мГал

    (5.99)

    Коэффициенты γ e , β , β

    1 относятся к так называемым фундаментальным

    геодезическим постоянным;

    γ e  978016, мГал – значение нормальной силы

    тяжести на экваторе; β  0 , 005302 – коэффициент порядка сжатия эллипсоида. Коэффициент β1 относится к постоянным высших порядков, которые уверенно определяются при задании параметров нулевого порядка и порядка сжатия. В формуле (8.1) коэффициент β1 будем считать равным 0,0000071. Таким образом, в формулу (8.1) входят два неизвестных коэффициента: γ e , β . Для их определения нужно иметь значения γ 0 в двух точках. Но так как нормальную силу тяжести нельзя получить по результатам измерений, то для нахождения неизвестных γ e и β используют измерения действительной силы тяжести в точках поверхности Земли.

    Выразим значение силы тяжести γ 0 на поверхности уровненного эллипсоида через результаты измерений и разрешим уравнение относительно

    gСВ:

     γ e 1  βsin 2 B  β 1 sin 2 2 B 14  g  δg 1

      g  γ  СВ 14

    γ 0

    (5.100)

    g СВ  γ e  γ e βsin 2 B   g 978 000 δg 1  γ e β 1

    sin 2 2 B 

    978 000 ,

    (5.101)

    где

    gСВ

    аномалия силы тяжести в редукции за свободный воздух;

    δg 1

    поправка за высоту точки над уровнем моря;

    g

    измеренное значение силы тяжести.

    Введем обозначения:

    a

    i

    1 ,

    b

     sin

    2

    B

    ,

    i

    i

    2

    2 B i ),

    l i  ( g i  978 000  δg 1 i  γ e β 1 sin

    ν i

    

    g СВ ,

    (5.102)

    x 1

     γ e

     978 000 ,

     γ e β .

    x 2

    Составим уравнения поправок в виде

    υ i

    a i

    x 1

    b i

    x 2

    l i ,

    (5.103)

    где i

    номер точки с измеренным значением силы тяжести, i = 1, …, n;

    n

    число точек с измеренным значением силы тяжести.

    Уравнения поправок запишем в матричном виде:

    Ax L υ,

    (5.104)

    где векторы (матрицы-столбцы) –

    υ1

    l

    υ 2

    1

    l2

    x1

    υ 

    , L

    (5.105)

    , x

    ,

    x 2

    υ n

    ln

    а матрица –

    a1

    b1

    a2

    b2

    A

    (5.106)

     .

    an

    bn

    2.

    Составить систему нормальных уравнений под условием υυ   min.

    Число нормальных уравнений равно числу неизвестных:

    Nx l 0,

    (5.107)

    где матрица

    N AT A

    aa

    ab

    ;

    (5.108)

    ab

    bb

    AT

    a

    a

    ...

    a

    (5.109)

    1

    2

    n

    ,

    b1

    b2 ...

    bn

    а вектор свободных членов

    l

    AT L

    al

    (5.110)

    .

    bl

    3. Определить искомые неизвестные, решив систему нормальных

    уравнений:

    x N 1 l Ql,

    (5.111)

    где Q – обратная по отношению к N матрица;

    Q11

    Q12

    (5.112)

    Q N 1

    .

    Q12 Q 22

    1. Определить коэффициенты формулы нормального распределения силы

    тяжести:

    γ e

     ( x 1  978 000), мГал;

    β 

    x 2

    .

    (5.113)

    γ e

    5. Определить сжатие уровенного эллипсоида в соответствии с теоремой

    Клеро:

    α 

    5

    q β,

    (5.114)

    2

    где

    q

    ω 2 a

    ;

    γ e

    ω – угловая скорость вращения Земли,  = 7,29  10-5 рад/с;

    a

    большая полуось эллипсоида Красовского a = 6 378 245 м.

    6. Определить сжатие уровенного эллипсоида через динамический

    коэффициент J 2 :

    α

    1

     3 J 2q .

    (5.115)

    2

    7. Оценить точность полученных результатов в соответствии с

    формулами:

    vv

    µ 

     µ Q11 ;

    ;

    m γ e

    n k

    m γ e β

    m β

    m γ e β

    (5.116)

    µ Q 22 ;

    ;

    m γ e

    m

    β

    m α

    ,

    α 2

    где Q11,

    Q 22 диагональные коэффициенты обратной матрицы N 1 Q .

    Пример выполнения задания

    Исходные данные (в соответствии с вариантом) приведены в табл. 5.15.

    Таблица 5.15. Исходные данные

    № п/п

    B ,о '

    ( g  978 000 ),мГал

    1

    5213'

    3 275

    2

    3648'

    1 851

    3

    422'

    73

    4

    2545'

    1 069

    5

    6717'

    4 193

    6

    4850'

    2 925

    1. Составим уравнения поправок.

    Величины δg 1 , вычисленные по формуле (5.66), выбираются из табл. 5.9 лабораторной работы № 5. Значения bi , l i , вычисленные по формулам (5.100–5.102), оформлены в виде таблицы (табл. 5.16).

    Таблица 5.16.Значения bi , l i

    B ,

    ( g  978 000),

     δg 1 ,

    γe β

    1 sin

    2

    2 B ,

    l ,

    γ e βb ,

    υ,

    b ,

    п/п

    

    мГал

    мГал

    мГал

    мГал

    мГал

    мГал

    мГал

    1

    52 13

    3 275

    1,5

    6

    –3 282

    3 189

    –13

    0,62463

    2

    36 48

    1 851

    118,5

    6

    –1 976

    1 832

    –62

    0,35883

    3

    4

    22

    73

    0

    0

    –73

    30

    38

    0,00580

    4

    25 45

    1 069

    –38,6

    4

    –1 034

    964

    0

    0,18874

    5

    67

    17

    4 193

    198,2

    4

    –4 395

    4 344

    30

    0,85087

    6

    48 50

    2 925

    26,2

    7

    –2 958

    2 894

    8

    0,56670

    Матрица A примет вид:

    1

    0 , 62463

    1. 0 , 35883

    1

    0 , 00580

    A

    1

    0 ,18874

    ,

    1. 0 , 85087

    1 0 , 56670

    вектор свободных членов –

     3 282

     1 976

     73

    L

     1 034

    .

     4 395

     2 958

    Уравнения поправок примут вид:

    x 1

     0,62463 x 2

     3 282  v 1

    x 1

     0,35883

    x 2

     1 976  v 2

    x 1

     0,00580 x 2

     73  v 3

    (5.117)

    x 1

     0,18874

    x 2

     1 045  v 4

    .

    x 1

     0,85087

    x 2

     4 395

    v 5

    x 1

     0,56670

    x 2

     2 958

    v 6

    2. Составим систему нормальных уравнений:

    N

    6

    2 , 59557

    ;

    2 , 59577

    1, 59967

     13 718

    l

     8 370

    ;

    6 x1  2 , 59557 x 2

    13 718

     0

    2 , 59557 x1 1, 59967 x 2

     8 370

     0

    .

    3. Вычислим неизвестные, решив систему нормальных уравнений. Решение нормальных уравнений возможно двумя методами:

    • Определение неизвестных с помощью программы MathCad;

    • Вычисление неизвестных с помощью определителей. Составим определитель D системы [27]:

    D

    n

    b

    n bb b  b 2 , 86105.

    b

    bb

    Неизвестные x1 и x 2 ,

    а также их

    весовые

    коэффициенты находят по

    формулам:

    x i

    D i

    ,

    D

    где D i

    определитель соответствующей величины.

    Имеем:

     l

    b

    D x 1

    13 738

    2 , 59557

    ;

      bl

    bb

    8 378

    1 , 59967

    x 1

    D x

    1

     76 , 7 мГал;

    D

     l

    6

    13 718

    D x 2

    n

    ;

    b

      bl

    2 , 59557

    8 370

    x 2

    D x

    2

     5 108 мГал;

    D

    b

    D P

    x 1

    1

      bb 1,59967 ;

    0

    bb

    1

    D P x 1

     0 , 5599 ;

    Px

    1

    D

    D P x

    n0

    n 6 ;

    2

    b

    1

    1

    D P x

    2

     2,097.

    Px

    2

    D

    Вычислим неизвестные с помощью программы MathCad:

    x  Ql

    0 , 559

     0 , 907

     13 718

    76 , 668

    ;

     0 , 907

    2 , 097

     8 370

    5 ,108 10 3

    x 76 , 7 мГал.

      1. 1 08

    1. Определим коэффициенты формулы нормального распределения силы тяжести:

    γ e

    x 1 978 000 978 077 мГал;

    β 

    x 2

     0 , 005222 .

    γ e

    5. Определим сжатия уровенной поверхности, приняв q  0 , 003468:

    • 5 q β, 2

    тогда α  0 , 003448, α  1 . 290

    6. Определим сжатие уровенного эллипсоида через динамический коэффициент J 2 :

    q 0 , 003468 ; α  0 , 0033623 ; α 1 . 297

    7. Оценим точность полученных результатов:

    µ 

     υυ 

     40 мГал,

    n k

    где n – количество измерений  n  6 ;

    k

    количество неизвестных  k  2 .

    µ

     µ

     30 мГал;

    m γ e

    Q 11

    P x 1

    µ

     µ

     58 мГал;

    m γ e β

    Q 22

    P x 2

    m β

    m γ e β

     0 , 000059;

    γ e

    m α

    m β

     5.

    α 2

    КОНТРОЛЬНЫЕ ВОПРОСЫ

    1.

    Сравните полученные значения γ e , β , α с современными значениями

    э тих величин. Объясните, почему полученные Вами коэффициенты отличаются от приведенных в справочной литературе [26, § 16, 46–48].

    1. Что относится к фундаментальным геодезическим постоянным?

    1. Какая существует связь между динамическими и геометрическими параметрами, характеризующими фигуру Земли?

    БИБЛИОГРАФИЧЕСКИЙ СПИСОК

    1. Астапов Ю.М. Земля и гравитация. – М.: Изд-во МГТУ им. Н.Э.

    Баумана, 1996. – 188 с.

    1. Астрономический ежегодник СССР на 1988 год. – Л.:Наука, 1988. – 690

    с.

    1. Боярский Э.А. Об оценке точности гравиметрических связей по формулам Буланже Ю.Д. // Изв. вузов, геодезия и аэрофотосъемка. – 1960. –

    • 6. – С. 45–51.

    1. Бровар В.В. Гравитационное поле в задачах инженерной геодезии. – М.:

    Недра, 1983. – 112 с.

    1. Буланже Ю.Д. Формулы для вычисления ошибок гравиметрической связи двух пунктов при многократных измерениях, выполненных группой

    гравиметров // Изв. АН СССР. Сер. Геофизика. – 1956. – № 7. – С. 755–764.

    1. Буланже Ю.Д., Арнаутов Г.П. Результаты первого международного сравнения абсолютных гравиметров, Севр, 1981 // Изв. АН СССР. Сер. Физика Земли. – 1983. – № 3. – С. 43–51.

    1. Веселов К.Е., Сагитов М.У. Гравиметрическая разведка. – М.: Недра,

    1968. – 512 с.

    8. Гладкий К.В. Гравиразведка и магниторазведка. – М.: Недра, 1967. –

    1. с.

      1. ГОСТ 13017–83. Гравиметры наземные. Общие технические условия. – М.: Изд-во стандартов, 1984. – 36 с.

      1. ГОСТ Р 52334–2005. Гравиразведка. Термины и определения. – М.: Стандартинформ, 2005. – 28 с.

      1. Гравиразведка: справочник геофизика / Под ред. Е.А. Мудрецовой, К.Е. Веселова. – 2- е изд., перераб. и доп. – М.: Недра, 1990. – 607 с.

      1. Грушинский Н.П. Теория фигуры Земли. – М.: Наука, 1976. – 512 с.

      1. Грушинский Н.П., Грушинский А.Н. В мире сил тяготения. – М.: Недра,

    1978. – 176 с.

    1. Грушинский Н.П., Сажина Н.Б. Гравитационная разведка: учеб. для техникумов. – 4- е изд. – М.: Недра, 1988. – 364 с.

    2. Еремеев В.Ф. Расчет палетки для вычисления высот квазигеоида

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