Определение нормального поля с использованием условия бровара
На правах рукописи
Романовский Сергей Иванович
Определение нормального поля
с использованием условия Бровара
Специальность 25.00.32 — Геодезия
Автореферат диссертации
на соискание ученой степени
кандидата технических наук
Москва 2008
Работа выполнена в Московском государственном университете геодезии и картографии (МИИГАиК) на кафедре высшей геодезии
| Научный руководитель: | кандидат технических наук, доцент Огородова Людмила Валентиновна |
| Официальные оппоненты: | доктор технических наук, профессор Баранов Владимир Николаевич |
| кандидат технических наук, доцент Крылов Виктор Иванович | |
| Ведущая организация: | Государственный астрономический институт им. П. К. Штернберга (ГАИШ МГУ) |
Защита состоится: 20 ноября 2008 г. в 12 ч. на заседании диссертационного совета Д212.143.03 при Московском государственном университете геодезии и картографии (МИИГАиК) по адресу: 105064, г. Москва, Гороховский пер., 4, МИИГАиК, зал заседания Ученого совета
С диссертацией можно ознакомиться в библиотеке МИИГАиК
Автореферат разослан: « ____ » ________________ 2008 г.
| Ученый секретарь диссертационного совета | Ю. М. Климков |
Общая характеристика работы
Актуальность темы исследования. Определение параметров общего земного эллипсоида, т. е. такого «… эллипсоида вращения, который наилучшим образом представлял бы многие астрономо-геодезические измерения, распределенные по земной поверхности…» (Гельмерт), относят к одной из основных задач геодезии. Разным аспектам решения этой задачи посвящены работы многих геодезистов. Ныне общий земной эллипсоид определяют как уровенный эллипсоид, внешнее поле которого определяют четыре фундаментальные постоянные. Три из них — геоцентрическая гравитационная постоянная 
, разность 
полярного 
и среднего экваториального 
моментов инерции и угловая скорость 
вращения — существуют у реальной Земли и в настоящее время известны с высокой точностью. Однако существует известный произвол в выборе и определении четвертого параметра, определяющего размеры уровенного эллипсоида. Постоянные, задающие размеры эллипсоида, не относятся к стоксовым постоянным Земли и не могут быть получены непосредственно по измерениям. В связи тем, что общий земной эллипсоид объективно не существует, его определение, по мнению М. И. Юркиной, является не только геодезической, а и философской проблемой.
В 1995 г. В. В. Бровар вместо прямых поисков оптимальной величины четвертого параметра уровенного эллипсоида предложил использовать общее требование, при котором «… возмущающая сила во внешнем пространстве достигает наименьшего значения…». В этом случае модель можно считать оптимальной, а общим земным эллипсоидом будет уровенный эллипсоид оптимального нормального поля. Требование Бровара записывают в виде условия, накладываемого на аномальный потенциал 
как гармоническую функцию,
, (1)
где 
— поверхность Земли, 
— вешняя нормаль к 
. Такой подход объединяет требования к нормальному полю, выдвигаемые при решении как геометрических, так и физических задач, и открывает возможность создания оптимальной универсальной модели. Однако до настоящего времени предложение Бровара не реализовано на практике.
Принятое в геодезии поле уровенного эллипсоида не является единственной моделью нормального поля. Наряду с ним в смежных областях естествознания используются и иные модели, более удобные для решения специфических задач различных наук. В настоящей работе предложение Бровара реализовано в широкой постановке и использовано не только для поисков общего земного эллипсоида, а и для построения оптимальных в смысле условия Бровара моделей нормального поля при разных формах его представления.
Степень разработанности проблемы. Попытки использования предложения Бровара предпринимались международным коллективом авторов (М. Бурша, М. И. Юркина, Г. В. Демьянов, Я. Коуба и др.). Однако их работы основаны не на условии, предложенном Броваром, а на видоизмененных условиях, соответствующих аномальному потенциалу с иными аналитическими свойствами. Кроме того, решение построено с учетом измерений только на поверхности океана; в этом случае условие (1) Бровара неприменимо. Задача определения нормального поля, реализующего условие Бровара, до настоящего времени не решена.
Цель исследования: построение оптимального в смысле условия Бровара нормального поля и определение параметров общего земного эллипсоида, гарантирующего минимум интеграла Дирихле во внешнем относительно поверхности Земли пространстве.
Задачи исследования. Для достижения цели исследования были поставлены задачи:
- получить удобные для минимизации выражения нормального потенциала и интеграла
для разных моделей нормального поля; - исследовать геометрию исходного поля уровенного эллипсоида на произвольных расстояниях от Земли и установить область существования эллипсоидальных уровенных поверхностей;
- найти удобную форму представления информации о рельефе поверхности Земли.
Объект исследования: поверхность и внешнее гравитационное поле Земли.
Предмет исследования: поле притяжения и поле силы тяжести эллипсоида вращения и диполя.
Теоретическая, методологическая и информационная база исследования. Теоретическая база — теория ньютоновского потенциала, методология — определение потенциала притяжения из решения краевой задачи. Задача определения уровенного эллипсоида была поставлена не только как задача определения общего земного эллипсоида, а в более широком смысле: было снято требование близости эллипсоидальной уровенной поверхности к поверхности геоида. Другими словами, уровенная эллипсоидальная поверхность не отождествляется с поверхностью нормальной Земли, а ищется модель нормального поля, наиболее близкого к действительному, в котором единственная возможная эллипсоидальная уровенная поверхность занимает произвольное положение в пространстве. Информационная база — модель EGM96 гравитационного поля Земли и цифровая модель GLOBE земного рельефа.
Научная новизна. Новые научные результаты, полученные в диссертации, состоят в следующем:
- доказана теорема, обобщающая теорему Маклорена – Лапласа;
- выяснена геометрия уровенных поверхностей на больших расстояниях от Земли;
- установлены точные соотношения между параметрами уровенного эллипсоида;
- определены параметры уровенного эллипсоида, гарантирующего минимум интеграла Дирихле;
- предложен новый способ определения уровенного эллипсоида по стоксовым постоянным Земли;
- получены формулы, удобные для вычисления потенциала эллипсоида и диполя;
- выполнено разложение рельефа до 360-й степени и проведен анализ этого разложения.
Теоретическая и практическая значимость работы. Работа носит исследовательский характер. Полученные результаты дополняют теории уровенного эллипсоида и нормального поля. Найденные выражения для нормального потенциала и силы тяжести предлагаемого нормального поля могут найти применение при вычислениях нормальной высоты и действительного потенциала Земли. Модель земного рельефа целесообразно использовать при глобальных геофизических и геодезических исследованиях.
Реализация результатов работы. Результаты исследования внедрены в учебный процесс МИИГАиК. Модель рельефа Земли в виде разложения по сферическим функциям включена в учебник по высшей геодезии.
Апробация результатов исследования. Результаты исследования докладывались на научных конференциях молодых ученых и аспирантов МИИГАиК (1998, 2001 и 2008 гг.) и на Сагитовских чтениях в ГАИШ МГУ (1997 г.).
Публикации. Основное содержание диссертации опубликовано в 1998 –2008 гг. в пяти статьях в журнале «Известия вузов. Геодезия и аэрофотосъемка» и в сети Интернет.
Структура и объем работы. Диссертация состоит из введения, пяти глав, заключения и списка литературы, включающего 66 наименований, содержит 19 таблиц и 19 рисунков. Общий объем диссертации — 148 страниц.
Основные результаты, выносимые на защиту:
- модели оптимального в смысле условия Бровара нормального поля при разных формах его аналитического представления;
- новый принцип определения уровенного эллипсоида;
- параметры эллипсоида, аппроксимирующего нормальный сфероид
оптимального поля; - обобщенная теорема Маклорена-Лапласа о потенциале притяжения симметричного эллипсоида;
- обобщенное условие Клеро, устанавливающее предельные значения параметров уровенного эллипсоида;
- модель рельефа Земли в виде разложения по сферическим функциям.
Основное содержание работы
Во введении сформулированы цель и задачи исследования — практическое применение предложения Бровара для определения нормального поля.
Первая глава содержит описание аналитических форм задания нормального поля и методов его выбора. Отмечается значение условия (1) Бровара и его преимущество перед использованными прежде условиями, накладываемыми на элементы нормального поля.
Во второй главе для потенциала 
притяжения эллипсоида вращения на основании решения краевой задачи с краевым условием
найдено выражение
, (2)
где 
и 
— эллипсоидальные координаты, 
и 
— полиномы Лежандра второй степени первого и второго рода соответственно, 
— линейный эксцентриситет (половина фокального расстояния), 
— потенциал на поверхности эллипсоида, 
и 
— коэффициенты. Потенциал притяжения определен тремя постоянными — 
, 
и 
— и не зависит от размеров эллипсоида. Этот результат сформулирован в виде теоремы, обобщающей теорему Маклорена-Лапласа: Софокусные эллипсоиды вращения равной массы с равными разностями моментов инерции и симметричным относительно оси вращения и плоскости экватора распределением плотности создают во внешнем относительно наибольшего из них пространстве одинаковое гравитационное поле.
Формула (2) позволяет получить как частные случаи потенциал притяжения однородного эллипсоида вращения и эллипсоида, уровенного для потенциала силы тяжести.
В третьей главе установлено, что для потенциала притяжения эллипсоида, уровенного для потенциала силы тяжести, коэффициенты 
и 
имеют значения
, 
,
где 
— малая полуось эллипсоида, 
— потенциал притяжения на его полюсе, а потенциал 
силы тяжести имеет вид
Потенциал силы тяжести определен постоянными 
, 
, 
и 
и зависит от размеров эллипсоида. При 
коэффициент при 
равен нулю, поэтому потенциал 
на поверхности эллипсоида постоянен. При любом другом значении 
потенциал зависит от широты 
, и поэтому только одна уровенная поверхность потенциала силы тяжести является эллипсоидом вращения.
Эксцентриситет 
связан с разностью 
моментов инерции и распределением массы внутри эллипсоида. Установлены соотношения:
(3)
где 
. Это точное соотношение, устанавливающее пределы изменения эксцентриситета 
в зависимости от распределения плотности. Оно эквивалентно известному соотношению Клеро. Нижний предел в соотношениях (3) соответствует однородному шару, верхний — однородному эллипсоиду.
При 
величины силы тяжести на оси вращения уровенного эллипсоида и в экваториальной плоскости связаны соотношением
. (4)
При 
эта формула переходит в формулу Молоденского, устанавливающую связь между силой тяжести на экваторе и полюсе уровенного эллипсоида.
Для общего земного уровенного эллипсоида, заданного первичными постоянными
= 3,986004418 · 1014 м3 с–2,
= 1,0826359 · 10–3, (5)
= 7,292115 · 10–5 с–1
и параметром 
= 6 378 137 м, производные постоянные принимают следующие значения:
= 6 356 752,26 м,
= 521 854,72 м,
= 62 636 852 м2 с–2
и экваториальная сила тяжести равна нулю при 
= 42 161,466 км, т. е. на расстоянии 
= 42 164,695 км от оси вращения. Потенциал в этой точке минимален и равен 
= 14 180 418 м2 с–2. Поверхность 
является последней (наиболее удаленной от Земли) замкнутой уровенной поверхностью. Она пересекает ось вращения на расстоянии 
= 28 107,650 км. В силу свойств потенциала уровенного эллипсоида эта поверхность не может быть эллипсоидом вращения.
Уровенный эллипсоид, сила тяжести на экваторе которого равна нулю, имеет параметры:
= 40 277,338 км,
= 27 030,407 км, (6)
= 29 860,024 км,
а потенциал на его поверхности 
= 14 023 117 м2 с–2. Как видно, значение эксцентриситета превышает значение малой полуоси.
Сила тяжести 
на полюсе эллипсоида (6) согласно формуле (4) при 
равна
= 0,43 м с–2 = 43 гал.
Разность силы тяжести на предельном эллипсоиде (6) значительно (на порядок) больше разности силы тяжести на поверхности Земли.
Уровенные поверхности земного гравитационного поля и поверхность эллипсоида (6) приведены на рис. 1.
Таким образом, параметры внешней эллипсоидальной уровенной поверхности потенциала силы тяжести нормальной Земли могут изменяться в пределах, представленных в табл. 1.
Среди этих эллипсоидов находится тот, поле силы тяжести которого наиболее близко к земному гравитационному полю.
Основные результаты, полученные в третьей главе:
- определены требования, которым должна удовлетворять разность
моментов инерции для существования уровенной поверхности в виде эллипсоида вращения и сферы; - получено соотношение, связывающее силу тяжести на оси вращения и в экваториальной плоскости, обобщающее формулу Молоденского;
- установлено точное соотношение между параметрами уровенного эллипсоида, соответствующее приближенному неравенству Клеро;
- установлены пределы изменения параметров уровенного эллипсоида в зависимости от разности
и значения потенциала на его поверхности.
Рисунок 1 — Уровенные поверхности нормального поля
Таблица 1 — Изменение параметров уровенной поверхности в зависимости от потенциала
| 62 636 852 м2 с–2 | >  > | 14 023 117 м2 с–2 |
| 6 378,137 км | <  < | 40 277,338 км |
| 6 356,752 км | <  < | 27 030,407 км |
| 521,855 км | <  < | 29 860,024 км |
| 9,832185 м с–2 | >  > | 0,449651 м с–2 |
| 9,780325 м с–2 | >  > | 0 |
| 298,26 | >  > | 3,04 |
Четвертая глава посвящена выбору нормального поля под условием Бровара в сферическом приближении, а также определению общего земного эллипсоида и нормального поля, заданного стоксовыми постоянными.
Вывод оптимальных в смысле условия Бровара параметров нормального поля в сферическом приближении. Вычисление интеграла (1) для реальной Земли даже для безошибочных измерений связано со значительными трудностями. Это обусловлено, главным образом, сложностью формы поверхности 
Земли, вследствие чего в (1) входит производная аномального потенциала не по направлению радиус-вектора или отвесной линии, а по наклонному к ним и нерегулярно меняющемуся направлению. Еще более сложной является задача определения параметров эллипсоида с одновременным уравниванием результатов измерений. Выясним поэтому сначала только те возможности, которые предоставляет изменение параметров эллипсоида. Интеграл 
преобразован к виду
, (7)
где 
, 
, 
— геоцентрические координаты, 
— высота, 
, 
— поверхность единичной сферы. Для оценки первого члена
правой части (7) аномальный потенциал 
представлен суммой
, (8)
в которой
для зональных коэффициентов,
для всех остальных; 
— коэффициенты разложения нормального потенциала. Центробежные потенциалы действительной и нормальной Земли считаем одинаковыми, поэтому аномальный потенциал является потенциалом притяжения. В сферическом приближении нормаль 
к поверхности Земли совпадает с радиусом-вектором 
, производная 
становится чистой аномалией силы тяжести и интеграл 
принимает вид
, (9)
где
— (10)
степенная дисперсия аномального потенциала. Согласно (9) и (10) интеграл 
как сумма квадратов коэффициентов всегда положителен. Это его свойство вытекает из выражения для интеграла Дирихле и подтверждается практическими вычислениями, выполненными Бурша и др.: для всех использованных ими точек поверхности Мирового океана знаки аномалии высоты, т. е. аномального потенциала, и чистой аномалии силы тяжести, т. е. производной 
, одинаковы.
Будем считать массы и разность 
моментов инерции действительной и нормальной Земли совпадающими, начало координат поместим в центр масс Земли. Тогда члены нулевой и первой степени в (8) и (9) и зональный коэффициент 
второй степени равны нулю. При этих условиях вычислим интеграл 
для двух моделей: поля уровенного эллипсоида вращения и поля двух неподвижных центров, используя для нахождения дисперсий 
коэффициенты модели EGM96.
Для коэффициентов 
потенциала притяжения уровенного эллипсоида найдены выражения:
Результат вычислений интеграла 
с этими коэффициентами представлен на рис. 2.
Рисунок 2 — Зависимость интеграла 
от линейного эксцентриситета уровенного эллипсоида и график функции 
Точки экстремумов кривой 
найдены с помощью условия
, (11)
где
.
Выражение (11) сходится очень быстро, стабильный результат получается уже при 
= 8. При этом члены с 
= 3 и 4 мало влияют на результат. Производная (11) содержит независящий от степени 
множитель
,
который входит во все члены ряда (11). Поэтому условие 
выполняется при
, (12)
откуда 
.
При 
= 1 082,626684 · 10–6 и 
= 6 378 137 м (EGM96)
= ±469 265,05 м.
Это значение эксцентриситета соответствует локальному максимуму интеграла 
.
Минимума интеграл 
достигает при
= 298 527,97 м (
= – 5 024 643,07 м)
и
= 592 349,15 м (
= 6 946 842,74 м).
В обоих случаях на полюсе эллипсоида (5) потенциал практически одинаков и равен 62 636 805,46 ± 0,42 м2 с–2, поэтому для представления нормального поля можно взять любое значение 
. Однако, поскольку первое из них приходится на область отрицательных полуосей, для уровенного эллипсоида нормального поля выбран эллипсоид с параметрами
= 592 349,15 м,
= 6 946 842,74 м,
названный минимальным эллипсоидом или эллипсоидом Бровара. Его большая полуось 
= 6 972 051,46 м.
Полуоси эллипсоида Бровара — уровенного эллипсоида, соответствующего минимуму интеграла Дирихле, — оказались примерно на 600 км больше полуосей общего земного эллипсоида. На экваторе эллипсоид Бровара проходит на высоте 593,92 км, на полюсе на высоте 590,09 км над Землей. Сжатие эллипсоида Бровара равно 1:277 и значительно больше сжатия земного эллипсоида.
При дипольном представлении нормального поля коэффициенты 
имеют вид
,
где 
— параметр, определяющий массы и положение центров диполя. Интеграл 
достигает минимума при 
= –0,035535671. Для центрально-симметричного диполя, у которого массы равны и расположены симметрично относительно начала координат, 
.
Таблица 2 — Сравнение интеграла 
для различных моделей нормального поля
| уровенные эллипсоиды | диполи | |||
| WGS-84 ПЗ-90 | Бровара | центрально-симметричный | несимметричный | |
 , 1013 м5 с–4 | 3,389 | 3,380 | 3,383 | 3,267 |
Значение минимума интеграла 
для несимметричного диполя меньше, чем при использовании эллипсоида Бровара и центрально-симметричного диполя (см. табл. 2). Это объясняется тем, что в первом случае в нормальный потенциал включены и нечетные зональные гармоники, т. е. учтена асимметрия северного и южного полушарий. Поле центрально-симметричного диполя и эллипсоида Бровара дают практически одинаковое значение минимума интеграла 
, увеличенное по сравнению с несимметричным диполем на 3,3%. Для эллипсоида системы ПЗ-90 интеграл 
больше по сравнению с несимметричным диполем на 3,6%.
Модели нормального поля, основанные на стоксовых постоянных. Поскольку коэффициенты 
разложения по шаровым функциям потенциалов моделей очень быстро убывают с возрастанием степени 
, для выбора нормального поля можно не определять минимум интеграла Дирихле, а обеспечить совпадение начальных коэффициентов модели и реального земного поля. Зональный коэффициент второй степени в разложении аномального потенциала равен нулю. Поэтому для уровенного эллипсоида изменение интеграла 
вызваны, в основном, изменением коэффициента 
четвертой степени
(13)
поскольку коэффициент 
не превышает 10–2 коэффициента 
и его изменение не может заметно повлиять на аномальный потенциал и интеграл 
. Зависимость квадрата коэффициента 
от эксцентриситета представлена на рис. 2. Как видно, характер изменения интеграла 
и квадрата коэффициента 
одинаков. Согласно (12) все коэффициенты достигают максимума при одном и том же значении эксцентриситета. Точки минимума найдены как точки пересечения кривой 
с осью абсцисс; в этих точках 
и из (13) получаем уравнение, связывающее линейный эксцентриситет 
эллипсоида с коэффициентами 
и 
, решение которого имеет вид
,
откуда
= ±298 291,89 м,
= ±592 824,86 м.
Эти значения близки к значениям, найденным из условия минимума интеграла Дирихле. Использование в качестве четвертого исходного параметра коэффициента 
снимает неопределенность в выборе уровенного эллипсоида, который вместо исходных постоянных 
, 
, 
, 
, или 
, 
, 
, 
, основан в данном случае только на реальных стоксовых постоянных Земли и угловой скорости ее вращения
, 
, 
, 
, (14)
где 
— зональная стоксова постоянная четвертой степени в разложении потенциала земного притяжения. Уровенный эллипсоид, основанный на четырех первичных постоянных (14), определяется однозначно; неточность его параметров явится следствием только ошибок этих постоянных.
Коэффициенты разложения потенциала притяжения центрально-симметриченого диполя не зависят от параметра 
и содержит только два параметра — 
и 
(или 
). Его нормированные коэффициенты определены только параметром 
.
Несимметричный диполь содержит уже три параметра; его коэффициенты включают параметр 
, который можно найти из условия равенства коэффициентов третьей степени диполя и Земли
= 0,957254173792 · 10–6,
отсюда 
= –0,035571550, что несколько отличается от значения, найденного из условия Бровара.
Определение общего земного эллипсоида. Найдем общий земной эллипсоид как эллипсоид, аппроксимирующий уровенную поверхность
(15)
поля эллипсоида Бровара, проходящую вблизи Земли, выбрав его под условием
, (16)
— радиус вектор уровенной поверхности (15), 
— радиус-вектор поверхности общего земного эллипсоида. Потенциал 
на поверхности 
не постоянен. На полюсе и экваторе потенциал эллипсоида Бровара меньше, а на широте 45° больше, чем потенциал 
на поверхности уровенного эллипсоида (5). В линейной мере разность 
составляет +4,7 м на полюсе, –1.9 м на широте 45° и +1,8 м на экваторе. Если же считать потенциал на поверхности эллипсоида (5) постоянным, то поле такого эллипсоида будет далеко от оптимального. Следовательно, соблюдение одновременно обоих требований к общему земному эллипсоиду, а именно
и 
для эллипсоида, находящегося вблизи Земли, невозможно. Эти условия выполнены только для эллипсоида Бровара, все остальные уровенные поверхности оптимального в смысле условия (1) поля эллипсоидальными не являются.
Подберем эллипсоид вращения к уровенной поверхности поля эллипсоида Бровара
= 62 636 857 м2 с–2 (17)
под условием (16). Экваториальный и полярный радиус-векторы уровенной поверхности (17) соответственно имеют значения
= 6 378 134,65 м,
= 6 356 747,05 м
и меньше полуосей эллипсоида (5). Параметры эллипсоида практически одинаковы при аппроксимации как уровенной поверхности оптимального поля, так и уровенной поверхности реальной Земли. Во всех случаях полуось 
совпадает с обычно используемой, а сжатие больше.
Результат аппроксимации уровенной поверхности 
= 62 636 857 м2 с–2 эллипсоидом вращения для разных моделей потенциала дан в табл. 3. На рис. 3 показано взаимное положение эллипсоида Бровара и общего земного эллипсоида, а также отступления от последнего уровенных поверхностей 
= 62 636 857 м2 с–2 поля эллипсоида Бровара (сфероп Бровара) и зональной части модели EGM96 (сфероид Мещерякова). Высоты этих двух уровенных поверхностей над общим земным эллипсоидом представлены в сильно увеличенном вертикальном масштабе.
Рисунок 3 — Взаимное положение общего земного эллипсоида, эллипсоида Бровара, а также его сферопа и сфероида по Мещерякову вблизи поверхности Земли
Таблица 3 — Параметры общих земных эллипсоидов, аппроксимирующих уровенные поверхности различных моделей поля силы тяжести
| Модель поля силы тяжести, уровенную поверхность которой аппроксимирует эллипсоид | Параметры аппроксимирующего общего земного эллипсоида | ||
 , м |  , м |  | |
| эллипсоид Бровара | 6 378 137,24 | 6 356 749,62 | 298,216269 |
| зональная часть EGM96 (только четные гармоники) | 6 378 137,45 | 6 356 749,11 | 298,206382 |
| полная зональная часть EGM96 | 6 378 137,47 | 6 356 749,10 | 298,205978 |
| зональная часть GEM-10C (только четные гармоники) | 6 378 137,00 | 6 356 748,71 | 298,207 |
Таблица 4 — Параметры уровенных эллипсоидов
| Параметр | Модель потенциала EGM96 | Уровенные эллипсоиды | ||
| WGS-84 | ПЗ-90 | Бровара | ||
 , 1014 м3 с–2 | 3,986004418 | 3,986005 | 3,9860044 | 3,986004418 |
 , м | — | 521 854,01 | 521 853,39 | 592 349,15 |
 , м | — | 6 378 137,00 | 6 378 136,00 | 6 972 051,46 |
 , 10–6 | –484,165372 | –484,166855 | –484,164953 | –484,165372 |
 , 10–6 | 0,539874 | 0,790304 | 0,790296 | 0,542423 |
 , 10–6 | –0,149958 | –0,001687 | –0,001687 | 0,000993 |
 , м2 с–2 | — | 62 636 861 | 62 636 861 | 57 326 349 |
 , м2 с–2 | — | 62 636 861 | 62 636 851 | 62 636 806 |
 , м2 с–2 | — | 62 636 861 | 62 636 852 | 62 636 871 |
 , м2 с–2 | — | 62 636 861 | 62 636 852 | 62 636 834 |
Примечание:  ,  и  — значения нормального потенциала на полюсе, широте 45° и экваторе эллипсоида (5) соответственно | ||||
В табл. 4 для сравнения приведены параметры поля эллипсоида Бровара, эллипсоидов основных моделей ГРС-80 и ПЗ-90 модели EGM-96. Как видно, уровенные эллипсоиды ГРС-80 и ПЗ-90 дают коэффициент 
в полтора раза больше, чем модель EGM96 и эллипсоид Бровара.
Несмотря на существенное различие полуосей эллипсоида Бровара и общего земного эллипсоида, сила тяжести оптимального нормального поля близка к силе тяжести на поверхности общего земного эллипсоида, рассматриваемого как уровенный. Разность значений нормальной силы тяжести рассмотренных моделей не превышает 4 мгл, то есть оптимальная модель дает близкие к земным значения силы тяжести.
Основные результаты, полученные в четвертой главе:
- выяснено, что наиболее близким к земному гравитационному полю является поле уровенного эллипсоида вращения, проходящего на высоте около 600 км над Землей;
- получено выражение для силы тяжести поля эллипсоида Бровара;
- показано, что возможности минимизации аномального поля при использовании нормального поля уровенного эллипсоида вращения ограничены. Меньшие значения интеграла
соответствуют дипольному представлению нормального поля. - получены параметры неуровенного общего земного эллипсоида, обеспечивающего наименьшее значение интеграла
.
В пятой главе выполнена приближенная оценка второго члена 
в выражении (7)
с использованием разложения высот 
физической поверхности Земли по сферическим функциям
.
Коэффициенты 
и 
найдены разложением по сферическим функциям до 360-й степени и порядка цифровой модели рельефа GLOBE.
Поправка 
представлена в виде 
,
,
где 
— максимальное значение аномального потенциала, 
— член, содержащий произведения полиномов Лежандра разных степеней и тригонометрических функций, обращающихся в нуль при нечетных значениях 
. Для 
принято значение 
= 103 м2 с–2, 
— нормальная сила тяжести, 
= 100 м — максимальная аномалия высоты, тогда
= 3,93 · 1011 м4 с–4 
= 8,31 · 1010 м5 с–4.
Значение 
составляет около 0,25% 
, поэтому можно ожидать, что учет влияния рельефа не приведет к качественному изменению полученных в диссертации глобальных оценок. Для более надежной оценки представляет интерес вычисление интеграла 
по наземным измерениям с выполнением интегрирования по реальной поверхности Земли.
Основные результаты пятой главы:
- получены ошибка представительства и дисперсия высоты рельефа Земли;
- выполнено разложение высот земного рельефа по сферическим функциям в четырех вариантах;
- найдены значения средней высоты земного рельефа, равное 232 м, и средней квадратической высоты — 636 м;
- выполнена приближенная оценка влияния рельефа на величину интеграла
.
Заключение
В диссертации представлен первый опыт применения предложения Бровара о выборе общего земного эллипсоида под условием минимума аномального гравитационного поля. Задача Бровара поставлена и решена в более широком смысле, как поиск не только параметров общего земного эллипсоида, но и определения нормального поля, заданного любой моделью. Для практического решения в качестве исходных данных о гравитационном поле Земли использованы коэффициенты модели EGM96, о поверхности Земли — цифровая модель GLOBE.
Основные результаты, полученные в диссертации, состоят в следующем:
- определен уровенный эллипсоид вращения, поле силы тяжести которого наиболее близко к земному. Этот эллипсоид расположен на высоте 600 км над Землей;
- получены параметры неуровенного общего земного эллипсоида, обеспечивающего наименьшее значение интеграла
: большая полуось 6 378 137,24 м, сжатие 1:298,216. Сжатие этого эллипсоида заметно больше сжатия эллипсоидов используемых ныне моделей; - предложен новый способ определения уровенного эллипсоида;
- показано, что возможности минимизации аномального поля при использовании нормального поля уровенного эллипсоида вращения ограничены. Меньшие значения интеграла
соответствуют дипольному представлению нормального поля; - доказана теорема, обобщающая теорему Маклорена-Лапласа;
- получено соотношение, связывающее силу тяжести на оси вращения и в экваториальной плоскости, обобщающее формулу Молоденского;
- установлено точное соотношение между параметрами уровенного эллипсоида, соответствующее приближенному неравенству Клеро;
- выполнено разложение высот земного рельефа по сферическим функциям в четырех вариантах.
Основное содержание диссертации опубликовано в статьях:
- Огородова Л. В., Романовский С. И. Применение интеграла Дирихле для выбора нормального поля. // Известия вузов. Геодезия и аэрофотосъемка, 1998, 6, с. 33 – 46.
- Романовский С. И. Сила тяжести оптимальной модели Земли. // Известия вузов. Геодезия и аэрофотосъемка, 2001, 1, с. 26 – 30.
- Романовский С. И. Некоторые характеристики цифровой модели GLOBE-98 рельефа Земли и ее спектральных представлений // Известия вузов. Геодезия и аэрофотосъемка, 2002, 1, с. 40 – 51.
- Огородова Л. В., Романовский С. И. Соотношения между параметрами уровенного эллипсоида. // Известия вузов. Геодезия и аэрофотосъемка, 2003, 4, с. 15 – 28.
- Огородова Л. В., Романовский С. И. Уровенный эллипсоид, основанный на стоксовых постоянных. // Известия вузов. Геодезия и аэрофотосъемка, 2008, 4, с. 3 – 6.
- Романовский С. И. Разложение высот цифровой модели рельефа Земли GLOBE в ряд по сферическим функциям // http://www.galaxy44.fatal.ru/files/globe98. rar.
