- •ГЛАВА 1. ПОВЕРХНОСТЬ ЗЕМЛИ И СИЛА ТЯЖЕСТИ
- •§ 1. Сила тяжести и ее потенциал
- •§ 2. Физическая поверхность Земли и геоид
- •§ 4. Геодезические прямоугольные системы координат
- •§ 5. Геодезическая эллипсоидальная система координат
- •§ 6. Сферическая система координат
- •§ 7. Специальная система координат сжатого эллипсоида вращения
- •§ 8. Натуральная система координат
- •§ 9. Связь натуральной и геодезической систем координат
- •§ 10. Топоцентрические системы координат
- •§ 11. Влияние движения полюса на координаты
- •§ 12. Международная служба широты и Международное условное начало
- •§ 13. Международная служба вращения Земли
- •ГЛАВА 3. НОРМАЛЬНАЯ ЗЕМЛЯ И ФУНДАМЕНТАЛЬНЫЕ ГЕОДЕЗИЧЕСКИЕ ПОСТОЯННЫЕ
- •§ 14. Нормальный потенциал и нормальное поле. Способы выбора
- •§ 15. Внешний потенциал уровенного эллипсоида
- •§ 16. Представление потенциала уровенного эллипсоида в виде ряда
- •§ 17. Сила тяжести на поверхности уровенного эллипсоида
- •§ 19. Система координат в нормальном поле
- •§ 21. Фундаментальные геодезические постоянные
- •§ 22. Связь системы координат в нормальном поле с натуральной
- •§ 23. Связь элементов аномального поля с аномальным потенциалом
- •§ 24. Уклонения отвеса в геометрическом и физическом определениях
- •§ 25. Астрономо-геодезические и гравиметрические уклонения отвеса
- •§ 26. Топографические уклонения отвеса
- •§ 27. Топографо-изостатические уклонения отвеса
- •§ 28. Астрономо-геодезическая и гравиметрическая аномалии высоты
- •ГЛАВА 5. ПРИНЦИПЫ ОПРЕДЕЛЕНИЯ ФУНДАМЕНТАЛЬНЫХ ПОСТОЯННЫХ. ГЛОБАЛЬНЫЕ МОДЕЛИ ЗЕМЛИ. ОБЩЕЗЕМНЫЕ СИСТЕМЫ КООРДИНАТ
- •§ 29. Определение фундаментальных постоянных нулевого порядка
- •§ 31. Глобальные модели потенциала. Результаты определения фундаментальных постоянных. Современные модели нормального поля
- •§ 32. Глобальные модели рельефа
- •§ 33. Общеземные системы координат
- •ГЛАВА 6. РЕДУЦИРОВАНИЕ РЕЗУЛЬТАТОВ ГЕОДЕЗИЧЕСКИХ ИЗМЕРЕНИЙ
- •§ 34. Редукция угловых измерений
- •§ 35. Редукция линейных измерений
- •§ 36. Приведение линейно-угловой сети в единую систему координат
- •§ 37. Требования к точности геодезических координат для редуцирования
- •ГЛАВА 7. СИСТЕМЫ ВЫСОТ
- •§ 39. Геодезическая высота и методы ее определения
- •§ 41. Нормальная высота и аномалия высоты
- •§ 42. Связь геодезической высоты с нормальной высотой и аномалией высоты
- •§ 43. Нормально-ортометрическая высота и высота когеоида
- •§ 45. Определение разности нормальных высот
- •§ 46. Динамическая высота
- •§ 47. Связь уклонения отвеса и аномалии высоты
- •§ 49. Способы определения аномалии высоты
- •§ 50. Астрономическое нивелирование
- •§ 51. Астрономо-гравиметрическое нивелирование
- •§ 53. Связь приращений геодезической высоты, нормальной высоты и аномалии высоты
- •§ 54. Определение разности нормальных высот по спутниковым наблюдениям. (Астрономо-гравиметрическое нивелирование теллуроида)
- •§ 56. Вычисление гравиметрической аномалии высоты
- •§ 57. Вычисление аномального потенциала по дискретным измерениям силы тяжести
- •§ 58. Вычисление аномалии высоты и уклонения отвеса по дискретным измерениям силы тяжести
- •ГЛАВА 9. ГОСУДАРСТВЕННАЯ ГЕОДЕЗИЧЕСКАЯ СЕТЬ. ПОСТРОЕНИЕ И УРАВНИВАНИЕ
- •§ 59. Историческая справка о построении государственной геодезической сети России
- •§ 60. Точность измерений в государственной геодезической сети
- •§ 61. Определение эллипсоида Красовского. Система координат 1942 г.
- •§ 62. Уравнивание государственной геодезической сети
- •§ 63. Система координат 1995 г.
- •§ 64. Перспективы развития государственной геодезической сети России
- •§ 65. Начало счета геопотенциальных чисел и высот
- •§ 66. Водное нивелирование
- •§ 67. Океанографическое нивелирование
- •§ 68. Определение потенциала в начале счета высот
- •§ 69. Уравнивание нивелирной сети
- •§ 70. Необходимость учета геометрии поля силы тяжести в специальных геодезических работах
- •§ 71. Особенности редукционных вычислений в специальных геодезических работах
- •§ 72. Редуцирование результатов измерений в местную прямоугольную систему координат
- •§ 73. Высоты в локальной системе координат
- •§ 74. Определение уклонений отвеса в местной системе
- •СПИСОК ЛИТЕРАТУРЫ
- •АВТОРСКИЙ УКАЗАТЕЛЬ
- •ПРЕДМЕТНЫЙ УКАЗАТЕЛЬ
- •ИСПОЛЬЗОВАННЫЕ АББРЕВИАТУРЫ
- •ОГЛАВЛЕНИЕ
- •ЗАМЕЧЕННЫЕ ОПЕЧАТКИ
меры Δψ площади осреднения связаны со степенью N учитываемых членов соотношением Пеллинена (5.46).
Следовательно, карта на рис. 5.4 соответствует высотам гео ида, осредненным на площадках 10°х10°, и отражает общие гло бальные особенности геоида. Такой геоид принято называть пла нетарным.
На карте геоида можно выделить три простирающиеся в ме ридиональном направлении области. Большую часть Азии, вос точную часть Индийского океана и прилегающую часть Антарк тиды занимает область отрицательных высот. В этой области к югу от Индостана расположено максимальное понижение геоида, превышающее 100 м. Вторая область отрицательных высот зани мает Северную Америку, Мексиканский залив, восточную часть Тихого океана. В этой области в море Росса у побережья Антарк тиды расположен второй по величине минимум высоты, превы шающий 60 м. Европа, большая часть Африки, Индийский океан к югу от Африки и Атлантический океан находятся в области положительных высот. Максимальные высоты расположены в Северной Атлантике (свыше 60 м) и Индийском океане (свыше 40 м). Положительным высотам геоида соответствует также цент ральная часть Тихого океана, где в районе Гавайского архипе лага расположено повышение свыше 50 м. Практически на всей поверхности Антарктиды высоты геоида отрицательны. В общем поверхность геоида не повторяет форму физической поверхности Земли и сформирована распределением масс в глубоких слоях Земли.
Средняя квадратическая высота геоида, осредненная на площад ке 1°х1°, получаемая суммированием в выражении (5.47) до 180-й степени и порядка, равна 30,6 м.
§ 32. ГЛОБАЛЬНЫЕ МОДЕЛИ РЕЛЬЕФА
Для обобщенного описания поверхности Земли используют каталоги значений высот или разложение рельефа по сферическим функциям
N п
(5.48)
где <*«*;&* - коэффициенты разложения по нормированным фун кциям.
157
Разложение высот рельефа до 16-й степени и порядка первым выпол нил в 1922 г. австрийский геодезист А. Прей.1Степенные дисперсии Dnh разложения Прея приведены в таблице 5.6.
В геодезических приложениях часто используют разложение амери канского геофизика Каула.
Т а б л и ц а 5.6
Глобальные модели рельефа. Степенные дисперсии Dnh.
Дисперсия 404 275 м2
|
Прей |
Каула |
|
|
|
|
1922 г., |
1967 г., |
GLOBE-98, м2 |
Морская топо |
|
|
м2 |
м2 |
|
|
|
n |
|
|
графическая |
||
|
|
|
|
||
|
|
|
|
поверхность, см2 |
|
|
|
|
|
|
|
|
Точечные |
5° х 5° |
7,5'х 7,5' |
1° х 1° |
Рапп,1998 |
|
значения |
|
|||
|
|
|
|
|
|
0 |
50715 |
48400 |
53 670 |
53 685 |
- |
1 |
16144 |
18 857 |
10591 |
10588 |
654,34 |
2 |
13 624 |
17006 |
19506 |
19516 |
1624,09 |
3 |
23 809 |
30960 |
37 556 |
37 569 |
289,00 |
4 |
21796 |
39613 |
33719 |
33720 |
29,27 |
5 |
23198 |
34848 |
38952 |
38952 |
44,89 |
6 |
19821 |
24340 |
37145 |
37136 |
275,89 |
7 |
9336 |
17 357 |
19214 |
19202 |
158,76 |
8 |
9325 |
14318 |
13 830 |
13818 |
196,56 |
9 |
10309 |
|
20123 |
20097 |
104,65 |
10 |
7134 |
|
13788 |
13 759 |
79,21 |
11 |
7048 |
|
7989 |
7972 |
14,98 |
12 |
5325 |
|
9303 |
9271 |
13,47 |
13 |
3362 |
|
8319 |
8280 |
6,92 |
14 |
3915 |
|
3 782 |
3766 |
6,76 |
15 |
3721 |
|
4767 |
4739 |
6,92 |
16 |
3844 |
|
8296 |
8249 |
10,57 |
8 |
|
|
|
|
|
Σ"1 . |
137053 |
197 399 |
210512 |
210501 |
3272,797 |
1 |
|
|
|
|
|
16 |
181711 |
— |
286881 |
286634 |
3516,262 |
|
1Prey. A. Darstellung der Hohen und Tiefenverhaltnisse der Erde durch eine Entwiclung nach Kugelfuntionen bis zur 16. Ordnung. «Abhandlungen der Gesellschaft der Wissenschaften», V.ll, № 1, 1922, 1-29, Gettingen.
158
|
|
|
|
Продолжение табл. 5.6 |
|
|
Прей |
Каула |
|
|
|
|
1922 г., |
1967 г., |
GLOBE-98, м2 |
Морская топо |
|
|
м2 |
м2 |
|
|
|
п |
|
|
графическая |
||
|
|
|
|
||
|
|
|
|
|
поверхность, см2 |
|
Точечные |
5°х 5° |
7,5'х7,5' |
1°х 1° |
Рапп,1998 |
|
значения |
|
|||
|
|
|
|
|
|
30 |
|
|
|
|
|
Σ°· |
- |
- |
335 853 |
335 123 |
3563,556 |
1 |
|
|
|
|
|
π |
225 м |
220 м |
232 м |
232 м |
- |
Д а |
|
|
624 м |
623м |
0,60 м |
V/1=0 |
|
|
|
|
|
Наиболее полная современная глобальная цифровая модель высот земного рельефа- The Global Land One-kilometer Base Elevation (GLOBE-98) - создана международной группой специалистов в 1990-1998 гг. Эта модель содержит средние значения высот на рав новеликих по географической сетке трапециях размером 30" х 30", что соответствует 1 х 1 км на экваторе.
В таблице 5.6 приведены степенные дисперсии высот модели GLOBE-98 в двух вариантах разложения: при осреднении высот на трапециях 7,5'х 7,5' и 1° х 1°. Разложение выполнено до 360-й сте пени и порядка.
Средняя высота поверхности Земли над уровнем моря, равная составляет для модели GLOBE-98 232 м. Среднее квадратичес кое значение высоты суши над уровнем моря при учете гармоник до
п—30 составляет 623 м, при учете гармоник до п - 360 - 845 м.
Вразложении Прея и модели GLOBE-98 высоты морской то пографической поверхности равны нулю. Применение метода спут никовой альтиметрии позволяет найти высоты МТП над эллипсо идом, а после вычитания из них высот геоида, описываемых выра жением (5.47), получить высоты морской топографической поверхности над геоидом. Степенные дисперсии высот морской то пографической поверхности также приведены в таблице 5.6. Сред няя квадратическая высота морской топографической поверхнос ти над геоидом равна 0,60 м.
159
На рис. 5.5 приведены высоты морской топографической по верхности над геоидом модели EGM-96. Высоты меняются в преде лах от -2 до +1 м, наклоны МТП объясняют глобальные океани ческие течения. Так, Гольфстриму соответствует понижение высот морской поверхности от +0,3 м у восточного побережья Цент ральной Америки до -0,5 м в Северной Атлантике, а течению Куросио - понижение высоты вдоль побережья Японии примерно на 1 м. Значительно понижается поверхность океана к побережью Антарктиды.
0 |
30 |
60 |
90 |
120 |
150 |
180 |
210 |
240 |
270 |
300 |
330 |
0 |
Рис. 5.5. Морская топографическая поверхность относительно геоида EGM-96. Горизонтали проведены через 0,1 м.
§ 33. ОБЩЕЗЕМНЫЕ СИСТЕМЫ КООРДИНАТ
Общеземной системой (см. главу 2) называют систему, начало которой находится в центре масс Земли, а одна из осей совмещена с ее осью вращения. Из-за изменения оси вращения в теле Земли за ось Ζ прямоугольной геоцентрической системы координат вы бирают направление на Международное условное начало (МУН), ось X этой системы лежит в пересечении плоскости экватора в сис теме МУН и Гринвичского меридиана, ось У дополняет систему до правой.
160
Понятие земной системы координат имеет двойной смысл: теоретически это система, центр которой точно совмещен с выб ранной точкой, а одна из осей также точно совмещена с выбран ным направлением. Практически любая система координат реа лизована координатами опорных точек на поверхности Земли, которые получены с неизбежными ошибками1. В связи с этим воз можно существование нескольких реализаций одной и той же иде альной системы, т.е. нескольких общеземных систем координат.
Внастоящее время как международную общеземную исполь зуют систему WGS-84. В этой системе принято положение по люса на эпоху 1900-1905 гг., исправленное в 1984 г. МБВ (Между народное бюро времени) согласно новым формулам для нута ции, принятым Международным астрономическим союзом (MAC)
в1979 г. Такое положение МУН называют Условным земным полю сом (СТРConventional Terrestrial Pole). Оси системы WGS-84 пока заны на рис. 5.6
ВРоссийской Федерации используется, наряду с WGS-84 , гео
центрическая система ПЗ-90.
МУН [СТР (1984,0) МБВ]
z
WGS 84
Рис. 5.6. Общеземная система координат WGS-84
1В англоязычной литературе первую (идеальную) систему называют тер мином System, вторую (практическую) - Frame.
161
Системы WGS-84 и ПЗ-90 получены динамическим методом космической геодезии по разнородной информации с одновремен ным выводом моделей гравитационного поля Земли. При выводе этих систем использованы постоянные:
с= 299792458 мс-1,
GM = 398 600,44 MV 2,
ω= 7292115 · l C V .
Как уже упоминалось, в системе WGS-84 использован эллипсо ид ГРС-80. Параметры эллипсоида, определенные при выводе ПЗ-90, приведены в таблице 5.7; для сравнения здесь же даны пара метры эллипсоида WGS-84.
|
|
|
Т а б л и ц а 5.7 |
Параметры эллипсоида в системе ПЗ-90 и их сравнение с WGS-84 |
|||
Параметр |
ПЗ-90 |
WGS-84 |
|
GM, м с-2 |
398 600,44-109 |
398 600,5 · 109 |
|
а, м |
6378136 |
6 378 137 |
|
\:а |
298,257 839303 |
298,257 223 563 |
|
U0, |
MV 2 |
62636861,074 |
62636860,8497 |
Jr |
108 |
108262,57 |
108263 |
^20 ' 10 |
-484,164953 |
-484,16685 |
|
|
|
Системы ПЗ-90 и WGS-84 являются базовыми для ГЛОНАСС
(Глобальная Навигационная Спутниковая Система) и GPS (Global Positions System) соответственно.
Система WGS-84 реализована пятью постоянно действующи ми станциями GPS, достаточно равномерно распределенными по долготе. Кроме того в 1991 г. основана международная сеть GPS для геодинамики - International GPS Service fo r Geodynamics (IGS).
Несколько пунктов этой сети расположены на территории России. На территории СССР система WGS-84 закреплена 162 пунктами
Доплеровской геодезической сети (ДГС), созданной в 1984-1993 гг. с использованием ИСЗ доплеровской спутниковой системы TRANSIT Точность воспроизведения системы WGS-84 сетью ДГС по каждой из осей 1-1,5 м.
С 1989 г. проводятся наблюдения на пунктах геоцентрической координатной системы ITRS (International Terrestrial Reference System)
162
Международной службы вращения Земли IERS {International Earth Rotation Service). На пунктах ITRSIERS выполняют наблюдения спутников системы GPS, а также используют методы РСДБ, лазер ной локации ИСЗ, ДОРИС. Точность определения пунктов этой сети оценивают погрешностью менее 10 см. Так как координаты пунктов сети от одной эпохи к другой изменяются вследствие геодинамических эффектов и ошибок измерений, каталоги периоди чески обновляют с указанием года наблюдений, начиная с катало га ITRF-89 (YTRF-International Terrestrial Reference Frame).
Система ПЗ-90 реализована координатами следящих станций системы ГЛОНАСС и пунктами Космической геодезической сети (КГС), созданной в 1985 г. топографической службой Вооружен ных сил СССР. На территории Российской Федерации находятся около 30 пунктов этой сети при средних расстояниях между ними 1-3 тыс.км и ошибкой взаимного положения 20 см. Таким образом, существуют два варианта ПЗ-90 - ПЗ-90 (КГС) и ПЗ-90 (ГЛОНАСС). Точность реализации системы ПЗ-90 (средняя квадратическая ошиб ка отнесения к центру масс) около 2 м.
В таблице 5.8 указаны реализации общеземных систем координат.
Т а б л и ц а 5.8
Реализации общеземных систем координат WGS-84 и ПЗ-90
Система |
|
WGS-84 |
ПЗ-90 |
|||
Сеть |
PURS ДОРИС |
GPS |
IGS |
ETRS ДГС |
ГЛО |
КГС |
|
ERS |
|
|
|
НАСС |
|
Реали |
ITRF-89 и т.д. |
- |
- |
EUREF ДГС-93 ФАГС |
КГС-85 |
зация
Примечания. В таблице приняты сокращения:
ETRS - European Terrestrial Reference System - Европейская земная референцная система;
EUREF - European Reference Frame - Европейская референцная сеть; IGS - International GPS Geodynamics Service (Международная служба
изучения геодинамики с помощью GPS);
ITRF - IERS (земная геодезическая сеть МСВЗ);
ДОРИС - Доплеровская орбитографическая радиопозиционная ин тегрированная спутниковая система;
ФАГС - фундаментальная астрономо-геодезическая сеть.
В настоящее время точность реализации общеземной системы координат Международной службы вращения Земли составляет несколько сантиметров.
163
Переход от системы WGS-84 к ПЗ-90 выполняют по формулам (2.4), используя определенные по измерениям параметры перехо да. По состоянию на 2001 г. установлено:
х0 = 1,1 м; у0 = 0,3 м; ζ0 = 0,9 м; εχ = ε = 0; εζ= -0,82 · 10-6 радиан; т = -0,12 · ЮЛ
где х0, у0 ζ0 - координаты центра системы ПЗ-90 в системе WGS84; εζ - угол поворота системы ПЗ-90 вокруг оси Ζ.
Согласно этим данным центр системы WGS-84 смещен на 1,45 м относительно центра системы ПЗ-90, а разворот вокруг оси Ζ составляет 0,17" (рис. 5.7).
7
^ ПЗ-90
Так как системы WGS-84 и ПЗ-90 близки между собой, можно считать
ту0 = mz0 = 0; £ ^0 = ε ^ 0 = 0, |
(5.49) |
тогда формулы (2.4), (2.5) примут вид
( γ |
λ ( у |
Л |
г \ |
εΣ 0"( у |
Л ( х |
Л |
л W G S |
л п з |
|
|
л п з |
л о |
|
_
Y \V G S |
|
У п з |
Z W G S > |
\ |
Z |
|
1и ) |
- т - ε ζ 1 0 У п з |
_ Уо |
||
О |
О |
,Ζ Π3 , |
Л , |
|
|
164
% W G S * % П 3 = ~ X o |
+ |
т ^ П З |
+ |
П З’ |
|
У W GS - У П 3 ~ ~ У о |
+ |
т У п З |
~ |
Ε ζ ^ Π 3 > |
(5.51) |
Z\VGS ~ Z пз = -z0 + т гпзз.
Подставив значения параметров перехода, получаем формулы для перехода от ПЗ-90 к WGS-84:
X |
|
|
1 |
-0,82· 1O'6 |
|
0 |
X |
|
Ц - |
Y |
|
|
=(1-0,12· 10Г6) 0,82· КГ6 |
1 |
|
0 |
Y |
|
0,3 |
z |
W G S -84 |
0 |
0 |
|
1 |
z |
|
0,9 |
|
|
|
|
|
|
П З-90 |
|
|||
|
Для обратного перехода от WGS-84 к ПЗ-90 следует использо |
||||||||
вать транспонированную матрицу: |
|
|
|
|
|
||||
X |
|
|
1 |
0,821 (Г6 |
0 |
X |
|
|
U ‘ |
Y |
|
=(1+0,12-10^) -0,82· КГ6 |
1 |
0 |
Y |
|
+ 0,3 |
||
z |
П З -90 |
|
0 |
0 |
1 |
z |
W ^84 |
|
0,9 |
В этих формулах поправки в координаты выражены в метрах, углы - в радианах.
Изменение геодезических координат в системах WGS-84 и ПЗ-90 следует оценивать по формулам (2.18), используя получен ные по (5.50) или (5.51) разности прямоугольных координат. В ре зультате выражение (2.18) примет вид
Н ^ FIj-j'j —ηϊ[(X jи cos L + y s i n LS)cos В + sin B~\ —
- (x0cos L +y0sin L)cos В - z0sin В -
(5.54)
~ (a W G S ~ а П З ) + ^ \ am i e W G S ~ е П З ) + е П З (a W G S ~ а П З )]s'n2 B,
(M + H){BWGS Виз) ——т[_{Х[j3cos L + Yj13sin L)sm B Ζπ3 cos B] +
+ (x0cos L + y0sin L) sin В - z 0cos В +
+ ^ [ а П З ( e W G S ~ е П З ) + е П З ( a W G S ~ а П З )] S^n 2 Д
(N + H)cos B(LWGS - Lm ) = x0 sinL - yQcosL - £z(Xn3 cosL + Yn3sinL).
165
Подставив в эти формулы численные значения параметров пе рехода, получим
HWGS - Нпз = -0,12 · 10Л(*яз cosL + Ym sinL)cos В + Zn3 sin£] - - (1,1 cosL + 0,3sinL)cosZ? - 0,9sinZ? - 1 + 0,047sin2/?;
(M + H){BWGS - Bn3)= 0,12 · \0~*[(Xn3 cosL + Yn3 sinL)sin В -
- ΖΠ3 cosВ] + (1,1 cosL + 0,3sinL)sin£ - 0,9cosfi + |
0,047sin22?; |
(N + H)COSB(LWGS —Ln3)= |
(5.55) |
=l,lsinL - 0,3cosL + 0,82 · 10Γ6Ζ(ΧΠ3 cosL + Yn3 sinL).
Подсчет по этим формулам дает изменения координат для Пул ково (В = 59,772°; L = 30,328°; Н = 71 м):
HWGS ~ Нпз ~ -3,060 м; (М + H)(BWGS - Впз) —
=0,542 м; (N + H)cosB(LWGS - LU3) = 2,936 м.