Бесплатный автореферат и диссертация по наукам о земле на тему
Разработка методов вычисления нормальных высот по результатам спутниковых измерений в инженерно-геодезических работах
ВАК РФ 25.00.32, Геодезия

Автореферат диссертации по теме "Разработка методов вычисления нормальных высот по результатам спутниковых измерений в инженерно-геодезических работах"

На правах рукописи

004615878

КРАВЧУК ИВАН МИХАИЛОВИЧ

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

Специальность: 25.00.32 - Геодезия

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

- 3 ЛЕК 2010

Москва 2010

004615878

Работа выполнена на кафедре прикладной геодезии Московского государст венного университета геодезии и картографии (МИИГАиК)

Научный руководитель: доктор технических наук, профессор

Клюшин Е.Б.

Официальные оппоненты: доктор технических наук, профессор

Матвеев С.И.

кандидат технических наук Куприянов А.О.

Ведущая организация: Государственный специализированный проектный

институт (ГСПИ)

Защита диссертации состоится «16» декабря 2010 г. в 14-00 час. на заседании диссертационного совета Д 212.143.03 при Московском государственном университете геодезии и картографии по адресу: 105064, Москва, К-64, Гороховский пер., 4, МИИГАиК, зал заседаний Ученого совета.

С диссертацией можно ознакомиться в библиотеке МИИГАиК.

2010 г.

Ученый секретарь

диссертационного совета

^¿щ

Яимков Ю.М.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы:

Применение методов определения местоположения по сигналам глобальных «вигационных спутниковых систем GPS и ГЛОНАСС (ГНСС) для геодезических делей началось в России в начале 90-х годов прошлого века. Выявились их существенные преимущества по сравнению с традиционными геодезическими методами. К ним относятся широкий диапазон точностей (от десятков метров до миллиметров), независимость от погоды, времени суток и года, отсутствие необходимости взаимной видимости между пунктами, высокая автоматизация и, как следствие, оперативность, возможность работы непрерывно и в движении. Эти качества обусловили высокую производительность и экономичность ГНСС, что особенно заметно в труднодоступных и малообжитых районах, которые занимают значительную часть нашей страны.

В настоящее время по результатам измерений двухчастотными спутниковыми приемниками достигнута средняя квадратическая ошибка вычисления приращений координат, равная 3 мм+1-10"бД где D - расстояние между спутниковыми приемниками. Но превышения между этими же пунктами могут быть получены со средней квадратической ошибкой 10-30 мм, которая значительно возрастает при увеличении расстояния D на десятки километров.

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

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

В основе спутниковых методов лежит геометрический принцип измерений,

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

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

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

Цель работы: Основной целью диссертационной работы является разработка методов повышения точности определения нормальных высот по результатам спутниковых измерений на основе тщательных точностных исследований.

Научная новизна работы:

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

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

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

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

5. Разработан алгоритм вычисления для повышения точности определения нормальных высот по данным спутниковых измерений для решения прикладных за-щч геодезии.

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

Апробация работы: результаты разработок обсуждены и одобрены на 5-й Международной научно-практической конференции «Геопространственные технологии и сферы их применения». 11-12 марта 2009 г.; на 64-й Юбилейной научно-технической конференции студентов, аспирантов и молодых ученых МИИГА-иК, посвященной 230-й годовщине со дня его основания, 7-8 апреля 2009 г.; на Международной научно-технической конференции, посвященной 230-летию основания МИИГАиК.

Публикации: по теме диссертации опубликованы 5 научных трудов, из них два в журналах, рекомендованных ВАК.

Структура и объем диссертации: диссертационная работа состоит из введения, пяти глав с подразделами, заключения и списка литературы. Общий объем работы - 116 стр. Диссертация содержит 4 таблицы и 19 рисунков. Список литературы составляет 66 наименований.

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

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

ГЛАВА 1. АНАЛИТИЧЕСКИЙ ОБЗОР НАУЧНОЙ И ТЕХНИЧЕСКОЙ ЛИТЕРАТУРЫ ПО СПУТНИКОВОМУ НИВЕЛИРОВАНИЮ

На основании изучения литературных источников выполнен анализ современного состояния развития спутниковых технологий в России, производство которых начиналось при полном отсутствии нормативных документов, регламентирующих технологию и точность выполнения работ. Кроме того, отсутствовали методы оценки точности спутниковых измерений. Всё это привело в 90-х годах прошлого века, с одной стороны, к неоправданному использованию высокоточных двухчастотных приёмников для создания съёмочного обоснования в низкоточных работах, а с другой - к созданию ответственных геодезических построе-

ний, например городских геодезических сетей, по упрощённой технологии.

Самыми последними работами в области нормативно-правовой документации в России по применению спутниковых методов являются «Спутниковая технология геодезических работ. Термины и определения» (РТМ 68-14-01, 2001) и «Инструкция по развитию съемочного обоснования и съемке ситуации и рельефа с применением глобальных навигационных спутниковых систем ГЛОНАСС и GPS» (ГКИНП (ОНТА)-02-262-02, 2002), в которой вопросам высокоточного спутникового нивелирования не уделяется достаточного внимания.

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

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

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

ГЛАВА 2. НОРМАЛЬНЫЕ И ГЕОДЕЗИЧЕСКИЕ ВЫСОТЫ

Геодезическая высота некоторой точки является минимальным расстоянием между поверхностью эллипсоида и определяемой точкой. Совместно с плановыми координатами В, L геодезическая высота Я определяет пространственное положение точки относительно заданного эллипсоида. Эти координаты называют эллипсоидальными. Геодезическая (эллипсоидальная) высота имеет только геометрический смысл и не может быть непосредственна измерена, поскольку в точке поверхности Земли неизвестны ни направление нормали к эллипсоиду, вдоль которой нужно измерять высоту, ни положение отсчетной точки на эллипсоиде, которая к тому же физически недоступна. Вычислить геодезические высоты или их

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

Спутниковые методы позволяют определить пространственные прямоугольные координаты X, У, 2 точки. Эти координаты связаны с геодезическими формулами перехода. Спутниковое определение координат может осуществляться двумя методами: абсолютным, как правило, с точностью до 5 м или относительным, точность которого колеблется от первых дециметров до нескольких метров при кодовых измерениях и 310 мм при фазовых методах измерений. В абсолютном способе по измеренным величинам возможно вычислить значения геоцентрических координат, а в относительном методе по результатам измерений можно вычислить лишь приращения АХ, ДУ, А7. пространственных координат. Наиболее эффективно проблема повышения точности решается за счет применения относительных методов спутниковых измерений.

Таким образом, самый распространенный на данный момент способ определения геодезических высот - это спутниковое нивелирование. До появления этой технологии, все высоты получали путем измерений на поверхности Земли. Полученные таким образом высоты носят название гипсометрических. Для связи гипсометрических и геодезических высот имеется соотношение Н" = //'""с + С, где Н' — геодезическая высота точки, Н'""с - гипсометрическая высота (высота в физическом понимании), £ - геоидальная составляющая высоты.

Нормальная система высот была предложена М.С. Молоденским и принята на территории бывшего СССР в 1945 г. как вспомогательная, необходимая для приближенного определения физической поверхности Земли, но далее стала широко применяться в геодезической практике. Нормальной высотой является высота над эллипсоидом в нормальном поле такой точки, для которой разность нормального потенциала относительно эллипсоида тождественно равна геопотенциальному числу точки поверхности Земли. Нормальные высоты практически не зависят от выбора эллипсоида и определяются по результатам измерений, выполненных на физической поверхности Земли. Однако они связаны с выбором начальной точки, от которой выполняется геометрическое нивелирование и ведется счет гео-потенциапьных чисел. В России отсчет высот принято вести от нуль-пункта Крон-штадского футштока. Нормальная система высот выгодно отличается от ортомет-рической тем, что не требует гипотез о распределении масс внутри Земли, поэто-

му вычислить значение высоты можно достаточно точно, а качество полученных высот зависит в основном только от качества измерений. Нормальные высоты могут быть определены из результатов нивелирования, а также вычислены по геодезическим высотам, если есть информация о разностях (//' - //'), то есть высотах квазигеоида С,, поскольку сумма нормальной высоты и аномалии высоты дает геодезическую высоту точки поверхности Земли: Я" = 1Р+С,

Аномалия высоты возникает из-за того, что гравитационное поле Земли не совпадает с нормальным, поэтому нормальная высота Н' не равна геодезической Н\ Поскольку нормальный потенциал всегда выбирают близким к действительному, аномалия высоты будет сравнительно малой величиной. Если отложить нормальные высоты по нормали к эллипсоиду от каждой точки поверхности Земли с учетом знака высоты вниз или вверх, то геометрическое место точек концов этих отрезков образуют поверхность, которую называют квазигеоидом. Следовательно, высота квазигеоида над эллипсоидом равна аномалии высоты.

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

ГЛАВА 3. ОБЗОР МЕТОДОВ СПУТНИКОВОГО НИВЕЛИРОВАНИЯ

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

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

обеспечивает точность 3-5 м, по фазовым измерениям может достигать средней квадратической ошибки, равной 1 мм.

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

Одним из методов спутникового нивелирования является спутниковая альтиметрия, под которой понимается измерение высоты спутника относительно поверхности Земли по времени прохождения сигнала. Спутниковое нивелирование и спутниковую альтиметрию объединяет то, что в обоих методах по результатам измерений необходимо вычислять геодезические и нормальные высоты объекта наблюдения. Суть этого метода заключается в следующем. Радиовысотомер, установленный на борту ИСЗ, измеряет высоту мгновенного положения ИСЗ над средним уровнем океана. Отраженный сигнал в идеальных условиях прохождения тем же кратчайшим путем возвращается в приемник радиовысотомера. Если на моменты измерения высот спутниковым радиовысотомером вычислять геодезические координаты мгновенных положений ИСЗ, то возможно для некоторой площадки геоида определить высоту геоида над эллипсоидом. Отклонения (положительные или отрицательные) геоида относительно эллипсоида называют ондуляциями геоида.

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

ляющие судить о течениях и других процессах.

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

ГЛАВА 4. ОЦЕНКА ТОЧНОСТИ СПУТНИКОВОГО НИВЕЛИРОВАНИЯ

Целью исследований, выполненных в данной главе, является разработка формул оценки точности вычисления координат геодезической высоты и разности геодезических высот. Для выполнения этой задачи рассмотрено вычисление геодезической высоты по прямоугольным координатам. Формулы, связывающие криволинейные координаты с прямоугольными декартовыми координатами: X = (N + Н)eosBcosL; Y = (N + H)cosBsmL\Z = (N + tf )sinB-e2NúnB, (1)

где В - геодезическая широта пункта; L - геодезическая долгота пункта; Я - гео-

„2 1,2 2 а -о „

дезическая высота пункта; е =----первый эксцентриситет; N - радиус кри-

а

визны первого вертикала; N = . а ; а - большая полуось эллипсоида;

VI -е2 sin2 В

b - малая полуось эллипсоида.

Геодезическую высоту можно вычислить несколькими путями. Суммируя квадраты первых двух уравнений (1), нетрудно получить:

Н = ——N, rj\eR = \lx2 + Y1. (2)

eos В

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

Я = ——(l-e2W. (3)

sin В v > w

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

H = (N + Я)eos2 В + (N + Я^п2 В -e2Nsin2 B~n( 1-е2sin2 В).

Это тождество легко проверить, раскрыв скобки в правой части. Используя ранее приведенные формулы (1), нетрудно получить

Я = 4х2 + Y2 eos В + Z sin В - aVi - е2 sin2 В. (4)

Формула (4) свободна от недостатков формул (2) и (3). Особого внимания заслуживает формула Боуринга, которая позволяет получить достаточно точное значение широты, минуя итерационный процесс. Так как формула Боуринга публикуется в отечественной литературе с опечатками, в диссертации приведен полный вывод этой формулы. В окончательном виде формула Боуринга имеет вид:

Z+óe sin 0 ,, е2

■-:—г—, где е =— ,

й-ое eos 8 1-е2

tgВ=--——, где е =-- второй эксцентриситет. (5)

Из формулы (5) следует: tg9=—, (6)

W

Точность формулы (5) достаточна практически для всех геодезических задач.

Для северных широт из уравнения (6) следует

. 0 __2 1 й(1-е2)"2

51п0=-р^==-т====; со5б=-;-=—г= (7)

Ф+Щ2е ^г+Н2(\-е2) р2+Я2(\-е2)

С учетом (7) формулу Боуринга (5) можно представить в виде, удобном для вычисления геодезической высоты через декартовы координаты

-11/2

ff=[x2+y2+Z2(l-alfí2r~3)~2(]+blZ2r~3)2J -

е Z R

(8)

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

Следует особо подчеркнуть, что формула вычисления геодезической высоты (8) не обладает недостатками формул (2) и (3), она также универсальна, как и формула (4) и может с успехом применяться как на полюсе, так и на экваторе.

Для оценки точности вычисления геодезической высоты по результатам спутниковых измерений формула (8) представлена в следующем виде:

Н=Р-И, (9)

где N - радиус кривизны первого вертикала, выраженный через прямоугольные

координаты; Р=^Х1+У2+21(\-ахР2г'3)'2^Ь122г-3)2. (10)

Полный дифференциал выражения (8) имеет вид:

.„ дР дИ дР ЗЫ^дР- дИ _

йН-—йХ--(IX +—<1У--йУ+—¿2--с12. (11)

дХ дХ дУ дУ д2 82

Для упрощения дальнейших преобразований введены обозначения

Г,=1-а,Л2г"3; (2=\+Ь122г~3, следовательно, Р=^Х2+У2+22022 .

После нахождения частных производных окончательно полный дифференциал (11) имеет вид:

йН

X 3Ь,(1-е%Х24 2а/2Х22 За.(\-е2У2Х22 Я2 ае2Х2

I —12------- '-г----- |

ЯЙГ+

Р Р(У Р(У Р(У ^(14-11)312 нА

а(} хг1___2аа,ЦгХ Заа^ф-е2)ХЯ2 ЗаЬхф3Х22

У 36,(1 -е%У2А 2а{\У22 За,(\-е2УгУ22Я2 ае2У22 ,Р Р(У + Р1у р(у Г Я4 (12)

М,У22___2аа,у3У Заа^ф-е2)УКг 3 аЬ^уУ22

Г 2% ЗЬ^г5 За^г3 ае2г__Ш,2

\Р12 Р(У р(у ^А-ц)У1к2+С{ичг)У2Я2

с!У+

7-3

Заа^Я2 ЗаЩУ2:

ег22 22 ?

й2, где

Воспользовавшись тем обстоятельством, что коэффициенты и /2 отлича-

212

ются от единицы на малую величину, член формулы (12) можно представить

Ч

в следующем виде:

21\ ^(Ьа^У3)2 2 2а,2Я2 2Ьх23 Р(\+Ь122г'3)1~р Рг3 Рг3 " С учетом (13) формулу (12) можно записать как

. ЗЬ.(1-е%Х24 2аДХ22 За{\-е2)йХ22К2 ае2Х22 где А=--1-г-^-+ 1 \ .---!--+-^--

Р(У Р13У Р(зу ^4-/3)3/2я4

Ш2Х22 2аа^Х _+Заа^^{\-е2)хК2 ЗаЬ^(}Х22

(13)

(14)

С^Гя* ФПи-^г3 ф44,3((4-13)тг5 ФГ(/4-03/2^

12

_ З^О-е2)^4 2 а&У2г За,(1-е2)/2Хг2Л2 аегУ22 £> =--

Л,/-' Р'У Р(У /4и(/4-/, )3,2л4

а/зкг2 2аа^цУ Ъаа^{\-е2)УЯ2 3 аЬ^%У2г

/4/3^ *\3/2Р4 / I3'2-3 ,2,4/3/, , \3'2„5 ,3,4/3/, , \3/2 5

4 \ 4 3/ Л Уа (.4 3 / ' 2 4 (Л"^] ^ 2 4 "'з) Г

2а^Я2 _2Ь]21 За^г3 ае2г

рг Ргъ пу р{у

Ъаа^цХВ2 ЗаЬ^^г3

(17)

Переходя от дифференциалов к конечным приращениям, а от них к средним квадратическим ошибкам, из (14) следует:

< J < Ч <■ 8)

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

тн=Щ\\-=5-+—+~+~^г+А+в+с- О9)

Х2+У2+22 2АХ 2ВУ 2С2 -+-+-+-

Р2 Р Р Р

Анализ формулы (19) показал, что с ошибкой, не превышающей 4%, она

X2 л-У2+22

может быть представлена в виде: тн=тк, так как -—-«1, а остальные

члены малы.

Например, в средней полосе В = 45°, 1= 37, #=200 м; Х= 3 608 020 м ; У=

v-2.y2.y2 <у ду ?/?У

2 718 839 м; 2= 4 487 489 м; . =1,023; ^-=-0,0048; —=0,0089;

рТ- ,,р,,р,г

1Г7

-=0,0397; Л2 = 0,000018; Л"=1,0237; В2 = 0,000108; С2 = 0,000791.

На экваторе В =0°, ¿=37°, //=200 м; Х=5 093 965м; 7=3 838 578м; г=0;

*2+722+22 = 1; ^=0; ^=0; Л2 = 0; Я2 = 0; С2 = 0.

Р РРР

В Заполярье В =72°, 1=37°, Я=200 м; Х= 1 578 909 м; У= 1 189 793м;

р2 р 2С7

—^—=0,1163; Л2 = 0,0060; Я2 = 0,00009; С2=0,0037.

Несмотря на сложную и нелинейную зависимость геодезической высоты от координат пункта, определенных по результатам спутниковых измерений (8),

г=6 043 875м; --г-=1,0012; —^-=0,0027; -^-=0,0036; #=1,0593;

средняя квадратическая ошибка вычисления геодезической высоты практически не зависит от положения пункта при равноточных декартовых координатах.

Для оценки точности вычисления разности геодезических высот по результатам спутниковых измерений в соответствии с принятыми обозначениями и формулой (9), можно записать:

Ь=Нг-Нх=Р2-Р\+Их-Иг, (20)

где Р,=

р2=

(1 -а^У \

Х2+Г2+22

(1

ег2{ВГ{

<?АК

(1

{иъ^)2 Ч (1+№г?у

По аналогии с выводами формул оценки точности вычисления геодезической высоты (11), полный дифференциал разности геодезических высот имеет вид:

(21)

Координаты второго пункта (Х2, У2, 22) вычисляют через координаты первого пункта (Лг,, 2^ и приращения координат (ДА1, А У, А2), вычисленные из результатов фазовых измерений. Следовательно,

йХг =с1Х, +с!АХ- с1У2 =с1У, +ЛДГ; с12г +с!А2. (22)

С учетом (22) выражение (21) примет вид:

й7г=

Р Р

М 2

У У

р р \ Г1 ' 1

2 2 ^ Рг

с12. ■

X.

^+В2)с1АУ-{^-+С2 \clA2.

(23)

Переходя от дифференциалов к конечным приращениям, а от них к средним

квадратическим ошибкам, из выражения (23):

"\2

ОТ,

X, X.

Р Р

. Г1 2

1Г+Л-А

Р2 г) ■ и рг 2)

X2 . —+А,

Р

V 2

-±+В2 | тм +

(24)

—-С2

Р

УГ2

тл

Для случая равноточных определений координат пункта 1

(т2^ =тг^ =т2к} и приращений координат (т\х=т\г=т12=т1^) формула оценки точности разности геодезических высот (24) имеет вид:

Рх Рг

Р Р

'I 2

Р Р

тг

+

X,

М | +

Р -г) {Р2 2

(25)

Как показал анализ формулы оценки точности (19), с высокой степенью точности величина

И

+

| +

^2

4^2

1-с,

Следовательно, выражение (25) преобразуется к виду:

тк =

*21 2 + [У, 2 + Г*, 2гЛ

и Рг) и Рг) и Рг)

У, У-,

г, г,

+2{В-ВЛ +2(С-СЛ

К Рг) 7и Р2У

(26)

(27)

Анализ членов, составляющих формулу (27), выполнен при условии Р2=Р1+АРГ Такой подход обеспечивает минимальное значение поправки АР. После преобразований величина АР равна:

72

АП Х.АХ У.АУ 2Ш г22(5,-53) АХ2+АУ2+Аг2 АР=—-+-1—+—1—+ 24 2—-

Р, Р, Ъ Рх ^

а формула для оценки точности (27) преобразуется к виду:

(28)

где Дт=2 (А,-А2)

КРХ Рг) 11 2\/{ />,

+2(С,-С2)||-|^

+{А-А2)2+{В-В2)2+{С-С2)2+

2ХЛРАХ 2УАРАУ 2гЛРАИ -2-+-2-+-2-

+Б2АР2

(29)

(30)

Б =АХ +АУ +А2. .

(31)

Величина Ат является незначительной, и ею можно пренебречь. Например, при В = 45°, Ь = 37°, Я=200 м; Ат=0,000674.

15

Окончательную формулу для оценки точности разности геодезических высот можно представить в виде:

Щ=^гп1+т1, где R

см5

О 10 20 30 40 50 60 S, км Рисунок 1. График средней квадратической ошибки определения разности геодезических высот

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

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

НИВЕЛИРОВАНИЯ

При использовании спутниковых методов определение координат выполняется в пространственной прямоугольной системе координат X, Y, Z, начало которой совпадает с центром общего земного эллипсоида. В настоящее время в глобальной системе GPS-позиционирования применяется система координат WGS-84, а в ГЛОНАСС используют ПЗ-90 (параметры Земли 1990 г.). Системы WGS-84 и ПЗ-90 получены динамическим методом космической геодезии по разнородной информации с одновременным выводом моделей гравитационного поля Земли. В системе WGS-84 использован эллипсоид GRS-80 (Geodetic Reference System, 1980). В инженерной геодезии используют координаты в референцной системе, например в СК-42, СК-95 или местной системе. Возникает необходимость пересчета координат из одной системы в другую.

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

з - средний радиус Земли. (32)

Приведенный анализ точности показал, что при вычислении разности геодезических высот существенное влияние оказывает ошибка координат опорного пункта. Средняя квадратическая ошибка определения координат опорного пункта в среднем равна тк = 5 м, а средняя квадратическая ошибка определения приращений координат /ид = 5 мм.

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

С помощью спутникового нивелирования определяется нормальная высота как сумма геодезической высоты, которая получается с большой ошибкой, и высоты квазигеоида, сведения о которой в полном объеме неизвестны. Целью данной работы является решение вопроса об использовании результатов спутникового нивелирования для решения геодезических задач прикладного (инженерного) характера без привлечения дополнительных (гравиметрических) измерений.

Формулы, определяющие зависимость геодезической высоты от значений декартовых координат, минуя итерационный процесс, имеют вид (8) или:

^х'+у^-а^г-'У+г^+ь^г-3)2 Н—-;—;--

г

-(1 -е2)а

Для перехода из одной системы в другую необходимо знать следующие параметры преобразования: смещение начала координат одной системы относительно другой (АХ0, ДУ0, три параметра, характеризующие разворот осей X, У, 2 в двух системах относительно друг друга (углы Эйлера \|/, 0, ф или углы Кар-дано ех, су, е;); масштабный множитель (коэффициент) т (в настоящее время применяется редко).

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

Это можно доказать на основе формул (8) и (33) нахождения геодезической высоты, для которой необходимо вычислить радиус кривизны первого вертикала N. Причем погрешность определения величины N полностью войдет в вычисляемые нормальные высоты, для этого необходимо произвести исследование изменения ра-

1__е2г2К'2_

227Г2+(1-я,Л2/-3)2(1+^2/"3)~2

диуса кривизны первого вертикала в зависимости от геодезической широты В. Формула вычисления радиуса кривизны первого вертикала имеет вид:

Л^ = а(1-Лт25)"1/2. (34)

В результате дифференцирования этого выражения получается:

ёИ = ^аг\т2В{\ - е\т2Ву™с1В. (35)

Из формулы (35) следует, что наименьшая зависимость радиуса кривизны N проявляется вблизи экватора (б=0°=>зт25=0) и полюсов (Д=±90°=>5т2Д=0), а максимальная зависимость - при В ~ 45°. Зависимость дифференциала геодезической широты от дифференциала дуги меридиана йХ,, имеет вид:

(1-е25т2яГ

йВ = ±-?-тт—ЛХ .

2(1-в>)

(36)

Следовательно, дифференциал радиуса кривизны первого вертикала йЫ связан с дифференциалом дуги меридиана с!Хм следующим соотношением:

(Ш = £ , а обратная зависимость имеет видаX = -г——¿М.

2(1 -е2) " е'ьШВ

На рис. 2 представлен график расчетов, показывающий, насколько могут отличаться координаты пунктов в направлении север-юг, при этом ошибка вычисления радиуса кривизны первого вертикала составляет ¿N=10 мм.

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

Рисунок 2. График зависимости с1Хм от геодезической широты В при с1Ы = 10 мм

бежная операция, которой следует уделять особое внимание.

Для того чтобы выполнить спутниковое нивелирование, необходимо иметь как минимум три опорных пункта (рис. 3) в национальной системе координат: геодезические широты и долготы относительно эллипсоида Красовского и нормальные высоты из высокоточного наземного нивелирования. На одном из пунктов определяются координаты спутниковым методом по кодовым измерениям (X, У, 2Л\'ск-з4 с ошибкой в 5 м, и вычисляются координаты двух остальных пунктов относительно первого по приращениям ДА^, АУ1ч, Таким образом, взаимное положение этих пунктов уже будет определено более точно с ошибкой около 3-5 мм.

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

По спутниковым измерениям геодезическая высота получается с ошибкой ~5 м - это означает, что, отложив по нормали к эллипсоиду измеренную высоту, получается отрезок 1-1"' вместо искомого отрезка 1-1". Расхождение измеренной геодезической высоты и нормальной высоты по каталогу одного и того же пункта 1 составит величину Д,. Вычислив по приращениям АХ^2, Д^-2, Д^м-г координаты пункта 2, найденная геодезическая высота этого пункта будет представлять собой отрезок 2-2"' (рис. 4). Необходимо определить нормальную высоту любого /-го пункта только по результатам спутниковых измерений. Если /-й пункт расположен в створе, который образуют пункты 1 и 2, то можно поступить следующим образом. Путем интерполирования вычисляется отклонение Д/ по формуле

Д.-^А/.+Д,,

(37)

где /, - отрезок Г"-2"'; /, - отрезок Г"-/"'.

г

Гчфепеиц--,^^^

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

Затем следует отложить найденное значение А, по нормали от точки /"'. Таким образом получается искомая нормальная высота, так как

Я", = И'. - А/. (38)

Следует отметить, значение Я", будет отличаться от действительной высоты

/-го пункта относительно квазигеоида на минимально возможную величину. Это связано с тем, что поверхность квазигеоида имеет сложную форму, и для его математического представления необходимо вводить геометрические оформляющие в виде аппроксимирующей прямой, плоскости или поверхности, которые наиболее близки к изучаемой поверхности в данном районе. С другой стороны, максимальное расстояние между пунктами, как правило, не превышает 20-30 км. Высоты квазигеоида в таких пределах изменяются плавно относительно отсчетного эллипсоида.

Для лучшей сходимости вычисленных и действительных высот пунктов целесообразно рассматривать аппроксимирующую плоскость, поскольку /-й пункт не всегда находится в створе опорных пунктов. Величину отклонения геодезической высоты от нормальной высоты /-го пункта можно найти как

Д, = 5о +8, + 82(1;-10, (39)

где Х:,У. - плоские координаты в проекции Гаусса-Крюгера (удобнее использовать); 50, 5], 82-параметры оформляющей плоскости.

Согласно рис. 4 и формулам (38) и (39):

Д, = /Л - Л". = §о; Л2 = /Л - Н"2= 50 + 5, [Х2 - Х)+ 52()^ - ; д3 = я"3 _ //'3 = 50 + 5, (X, - X,) + 52 (?3 - ,

или

н; - н; = щ - н;+5, (х2 - +52 (72 -щ - н; = щ - н;+5, -)+82 (у; - ?)}''

Решая систему уравнений (40), можно вычислить параметры оформляющей плоскости §1 и б2. Теперь поправку Д,- к геодезической высоте, определенной по спутниковым наблюдениям на пункте с номером /, можно найти по формуле

4=//,-//', +Щ-?).

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

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

2. Определение координат одного из пунктов в системе \VGS-84 или ПЗ-90 при использовании спутниковых навигационных систем. Эти координаты определяются только один раз и остаются неизменными на весь цикл работ, и по ним вычисляются координаты определяемых пунктов, для которых необходимо вычислить нормальные высоты. Назовем этот пункт опорным.

3. По результатам фазовых измерений определяются приращения координат до пунктов с известными координатами в национальной или местной системе координат.

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

5. Вычисляются параметры преобразования координат из общеземной системы координат в национальную или местную систему координат.

6. Вычисляются параметры оформляющей плоскости (при наличии трех общих пунктов) или поверхности, если количество общих пунктов более трех.

7. Выполняются полевые измерения для вычисления приращений координат между опорным пунктом и определяемыми пунктами. Используя ранее опреде-

ленные координаты опорного пункта, вычисляются координаты определяемых пунктов.

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

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

ЗАКЛЮЧЕНИЕ

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

В диссертации доказано, что относительно невысокая точность вычисления разностей геодезических высот связана не с особенностями условий пространственной засечки, а в большей степени зависит от необходимости использовать абсолютные координаты пунктов, которые определяются с невысокой точностью (3-5 м).

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

В настоящее время спутниковое нивелирование находит широкое применение в различных областях геодезии, в том числе и в инженерно-геодезических ра-

ботах. Итоговой целью способа, разработанного в диссертации, является получение нормальных высот пунктов. По результатам спутниковых измерений вычисляется геодезическая высота пункта, и для определения нормальной высоты необходимо в районе работ наличие пунктов с известными нормальными высотами. После обработки спутниковых наблюдений получаются прямоугольные координаты в системе WGS-84 или ПЗ-90 (при использовании спутниковых навигационных систем NAVSTAR GPS или ГЛОНАСС), которые существенно отличаются от координат пунктов с известными нормальными высотами в системе координат СК-42, СК-63 либо СК-95. Такие расхождения не могут не сказаться на точности определения нормальных высот вновь определяемых пунктов. Однако ни в научной, ни в технической литературе нет рекомендаций по вопросу преобразования координат при выполнении спутникового нивелирования. Помимо этого жесткая конкурентная борьба между фирмами привела к тому, что все программные продукты стали закрытыми от специалистов и представляются потребителям как конечная разработка, которую он вынужден покупать, полагаясь на авторитет фирмы.

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

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

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

Публикации по теме диссертации:

1. Гайрабеков И.Г., Кравчук И.М. Оценка точности вычисления геодезической высоты по результатам спутниковых измерений. Геодезия и картография. № б, 2010, с. 5-7.

2. Юпошин Е.Б., Кравчук И.М. Спутниковое нивелирование. Сборник статей по итогам Международной научно-технической конференции, посвященной 230-летию основания МИИГАиК. Выпуск 2, часть II, МИИГАиК, М.: 2009.

3. Клюшин Е.Б., Кравчук И.М. Теоретические исследования точности спутникового нивелирования. Материалы 5-й Международной научно-практической конференции «Геопространственные технологии и сферы их применения». 1112 марта 2009 г. Материалы конференции. М.: Информационное агентство «Гром», 2009, 92 с.

4. Кравчук И.М. Особенности вычисления нормальных высот по результатам спутниковых измерений. Известия вузов. Геодезия и аэрофотосъемка. № 4, М.: МИИГАиК, 2010.

5. Кравчук И.М., Семеонова А.Е. Спутниковое нивелирование. Доклад на 64-й Юбилейной научно-технической конференции студентов, аспирантов и молодых ученых МИИГАиК, посвященной 230-й годовщине со дня его основания. М., 7-8 апреля 2009.

Подписано в печать 08.11.2010. Гарнитура Тайме Формат 60x90/16. Бумага офсетная. Печать офсетная. Объем 1,5 усл. печ. л. Тираж 80 экз. Заказ №261 Цена договорная Издательство МИИГАиК 105064, Москва, Гороховский пер., 4

Содержание диссертации, кандидата технических наук, Кравчук, Иван Михайлович

ВВЕДЕНИЕ.

1. АНАЛИТИЧЕСКИЙ ОБЗОР НАУЧНОЙ И ТЕХНИЧЕСКОЙ ЛИТЕРАТУРЫ ПО СПУТНИКОВОМУ НИВЕЛИРОВАНИЮ.

1.1. Нормативная документация.

1.2. Научно-техническая литература.

1.3. Примеры определения высот с помощью спутниковых технологий.

1.3.1. Исследование точности спутникового нивелирования при коротких векторах.

1.3.2. Исследование точности определения высотных отметок с помощью спутниковых приемников в хвойном лесу.

1.3.3. Измерение высотной деформации покрытия аэродрома Шереметьево-2 с применением спутниковой аппаратуры.

1.3.4. Мониторинг оседания нефтяного месторождения с помощью системы GPS.

2. НОРМАЛЬНЫЕ И ГЕОДЕЗИЧЕСКИЕ ВЫСОТЫ.

2.1. Физический смысл высоты. Типы гипсометрических высот.

3. ОБЗОР МЕТОДОВ СПУТНИКОВОГО НИВЕЛИРОВАНИЯ.

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

3.2. Применение спутниковой альтиметрии.

4. ОЦЕНКА ТОЧНОСТИ СПУТНИКОВОГО НИВЕЛИРОВАНИЯ.

4.1. Вычисление геодезической высоты по прямоугольным координатам.

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

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

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

5.1 Использование спутниковых методов для определения нормальных высот.

5.2. Особенности определения нормальных высот по результатам спутниковых измерений.

Введение Диссертация по наукам о земле, на тему "Разработка методов вычисления нормальных высот по результатам спутниковых измерений в инженерно-геодезических работах"

Применение методов определения местоположения по сигналам глобальных навигационных спутниковых систем GPS и ГЛОНАСС (ГНСС) для геодезических целей началось в России в начале 90-х годов прошлого века. Выявились их существенные преимущества по сравнению с традиционными геодезическими методами. К ним относятся широкий диапазон точностей (от десятков метров до миллиметров на значительных расстояниях), независимость от погоды, времени суток и года, отсутствие необходимости взаимной видимости между пунктами, высокая автоматизация и, как следствие, оперативность, возможность работы непрерывно и в движении. Эти качества обусловили высокую производительность и экономичность ГНСС, что особенно заметно в труднодоступных и малообжитых районах, которые занимают значительную часть нашей страны. Кроме этого, преимуществом явилась возможность совместного позиционирования в трехмерном пространстве.

В настоящее время по результатам измерений двухчастотными спутниковыми приемниками достигнута средняя квадратическая ошибка вычисления приращений координат, равная 3 мм+1-10"6Д где D - расстояние между спутниковыми приемниками. Но превышения между этими же пунктами могут быть получены со средней квадратической ошибкой 10-30 мм, которая значительно возрастает при увеличении расстояния D на десятки километров.

Существует мнение, что низкая точность вычисления превышений связана с неудачной засечкой высот антенн спутниковых приемников, вызванной неудачным расположением спутников в период наблюдения. С этим мнением трудно согласиться. Например, на экваторе с долготами, равными 0°, 90°, 180° и 270°, высоты пунктов совпадают с осями прямоугольных координат X и У, но и в этих местах как координаты пунктов, так и разности координат вычисляются с той же точностью, что и в других регионах Земли. Наиболее вероятно, что причина относительно невысокой точности спутникового нивелирования заключается в особенностях вычисления геодезических и нормальных высот, которая объясняется принципиальным различием между традиционными и спутниковыми методами геодезии, которое состоит в следующем.

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

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

Для успешной реализации результатов спутниковых измерений и получения неискаженных координат пунктов в государственной геодезической системе координат и, особенно, для спутникового нивелирования, необходима точная модель высот квазигеоида. В России последний вариант модели высот квазигеоида над эллипсоидом Красовского, созданный методом астро-номо-гравиметрического нивелирования, был получен при завершении общего уравнивания астрономо-геодезической сети в 90-х годах прошлого века. Для большей части России средние квадратические ошибки высот квазигеоида относительно исходного пункта в Пулково не превышают 1 м, а для наиболее удаленных регионов (северная часть полуострова Ямал и восточные части Чукотки и Камчатки) ошибки доходят до 1,5 м из-за ошибки передачи высот в эти районы астрономо-геодезическим нивелированием. Превышения высот квазигеоида между пунктами, находящимися в пределах одного листа карты масштаба 1:1 ООО ООО, характеризуются средними квадратическими ошибками не более 0,2 м. Для решения региональных задач для сравнительно небольших территорий создавались модели высот квазигеоида по гравиметрическим данным, а также модели, сочетающие определения высот квазигеоида по спутниковым и нивелирным измерениям с гравиметрическими определениями высот квазигеоида. Для обеспечения геодезических работ вдоль железной дороги Москва-С.-Петербург, а также для геодезических работ в г. Москве были созданы отдельные модели высот квазигеоида. Следует особо отметить, что в настоящее время нет модели высот квазигеоида с точностью, необходимой для успешной обработки результатов спутниковых измерений и, особенно, для спутникового нивелирования. В результате нормальная высота, найденная с использованием спутниковых данных, существенно уступает по точности высоте, найденной из высокоточного геометрического нивелирования, которое является одним из наиболее точных видов геодезических работ.

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

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

- рассмотреть достоинства и недостатки метода определения нормальных высот с точки зрения использования спутниковых технологий;

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

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

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

Перед автором диссертации поставлена научная задача: на основе тщательных точностных исследований разработать методы повышения точности определения нормальных высот по результатам спутниковых измерений.

Данная работа содержит 116 страниц машинописного текста, состоит из пяти разделов с подразделами, включающих в себя 4 таблицы и 19 рисунков.

Заключение Диссертация по теме "Геодезия", Кравчук, Иван Михайлович

ЗАКЛЮЧЕНИЕ

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

Впервые доказано, что относительно невысокая точность вычисления разностей геодезических высот связана не с особенностями условий пространственной засечки, а в большей степени зависит от необходимости использовать абсолютные координаты пунктов, которые определяются с невысокой точностью (3 — 5 м).

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

В настоящее время спутниковое нивелирование находит широкое применение в различных областях геодезии, в том числе и в инженерно-геодезических работах. Итоговой целью способа, разработанного в диссертации, является получение нормальных высот пунктов. По результатам спутниковых измерений вычисляется геодезическая высота пункта, и для определения нормальной высоты необходимо в районе выполняемых работ наличие пунктов с известными нормальными высотами. Следует помнить, что после обработки спутниковых наблюдений получаем прямоугольные координаты в системе WGS-84 или ПЗ-90 (при использовании спутниковой навигационной системы NAVSTAR GPS или ГЛОНАСС соответственно), которые существенно отличаются от координат пунктов с известными нормальными высотами в системе координат СК-42, СК-63 либо СК-95. Такие расхождения не могут не сказаться на точности определения нормальных высот вновь определяемых пунктов. Однако ни в научной, ни в технической литературе нет рекомендаций по вопросу преобразования координат при выполнении спутникового нивелирования. Помимо этого жесткая конкурентная борьба между фирмами привела к тому, что все программные продукты стали закрытыми от специалистов и представляются потребителям как конечная разработка, которую он вынужден покупать, полагаясь на авторитет фирмы.

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

В результате проведенных исследований на тему повышения точности определения нормальных высот по данным спутниковых измерений для решения прикладных задач геодезии предложена следующая методика:

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

2. Определение координат одного из пунктов в системе WGS-84 или ПЗ-90 при использовании спутниковой навигационной системы NAVSTAR GPS или ГЛОНАСС соответственно. Эти координаты определяются только один раз и остаются неизменными на весь цикл работ, и по ним вычисляются координаты определяемых пунктов, для которых необходимо вычислить нормальные высоты. Назовем этот пункт опорным.

3. По результатам фазовых измерений определяются приращения координат до пунктов с известными координатами в национальной или местной системе координат.

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

5. Вычисляются параметры преобразования координат из общеземной системы координат в национальную или местную систему координат.

6. Вычисляются параметры оформляющей плоскости (при наличии трех общих пунктов) или поверхности, если количество общих пунктов более трех.

7. Выполняются полевые измерения для вычисления приращений координат между опорным пунктом и определяемыми пунктами. Используя ранее определенные координаты опорного пункта, вычисляют координаты определяемых пунктов.

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

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

Библиография Диссертация по наукам о земле, кандидата технических наук, Кравчук, Иван Михайлович, Москва

1. Антонович K.M. Использование спутниковых радионавигационных систем в геодезии. Т. 2.- М.: ФГУП «КАРТГЕОЦЕНТР», 2006.- 360 С.

2. Бородко A.B., Макаренко Н.Л., Демьянов Г.В. Развитие системы геодезического обеспечения в современных условиях. Геодезия и картография, 2003, № 10, с. 7-13.

3. Васильев A.C. Лупачев Ю.В. Характерные закономерности поведения среднего уровня морей России в период потепления климата. В сб.: Исследования океанов и морей, выпуск 207, Гидрометеоиздат, 2000, с. 131-141.

4. Временная инструкция по обследованию и восстановлению пунктов и знаков государственной геодезической и нивелирной сетей СССР. РИО ВТС, М.,1970, с. 23.

5. Гайрабеков И.Г., Кравчук И.М. Оценка точности вычисления геодезической высоты по результатам спутниковых измерений. Геодезия и картография. № 6, 2010, с. 5-7.

6. Генике A.A., Лобазов В.Я., Ямбаев Х.К. Результаты исследований аппаратуры спутникового позиционирования GPS Wild-System 200. Геодезия и Картография.-1993.- №10. С.8-13.

7. Генике A.A., Побединский Г.Г. Глобальные спутниковые системы определения местоположения и их применение в геодезии. Изд. 2-е, перераб. и доп.- М.: КАРТГЕОЦЕНТР, 2004.- 355 С.

8. Гофман-Велленгоф Б., Мориц Г. Физическая геодезия. М.: 2005.

9. Дронов A.B., Помогаев О.Н. (Hl 111 «Навгеоком») Исследование точности определения с помощью GPS высотных отметок в хвойном лесу, http ://www.navgeocom.ru/proj ects/5 700forest/index.htm.

10. Еремеев В.Ф. Теория ортометрических, динамических и нормальных высот// Труды ЦНИИГАиК. 1951. Выпуск 86. с. 11- 52.

11. Жалковский Е.А, Демьянов Г.В., Зубинский В.И., Макаренко П.Л., Пьян-ков Г.А. О концепции и программе перехода топографо-геодезического производства на автономные методы спутниковых координатных опреде-лений//Геодезия и картография. 1998, № 5.

12. Жданов Н. Д., Макаренко Н. Л. О концепции перехода топографо-геодезического производства на автономные методы спутниковых координатных определений. — Геодезия и картография, 1998, № 3).

13. Инструкция по вычислению нивелировок. Утверждена ГУГК и ВТУ. М.: Недра, 1971.

14. Инструкция по нивелированию I, II, III и IV классов. М., Недра, 1990.

15. Инструкция по нивелированию I, И, III и IV классов ГКИНП (ГНТА)-03-010-03, МОСКВА, ЦНИИГАИК. 2004.

16. Инструкция по построению государственной геодезической спутниковой сети ГКИНП (ГНТА) 01-006-03. - М., Федеральная служба геодезии и картографии России, 2003 г.

17. Инструкция по развитию съемочного обоснования и съемке ситуации и рельефа с применением глобальных навигационных спутниковых систем ГЛОНАСС И GPS. ГКИНП (ОНТА)-02-262-02. Москва, ЦНИИГАиК, 2002.

18. Инструкция по составлению технических отчетов о геодезических, астрономических, гравиметрических и топографических работах. М.: Недра, 1971.

19. Кашин Л.А. Построение классической астрономо-геодезической сети России и СССР. М., Картгеоцентр-Геодезиздат, 1999, 192 с.

20. Клюшин Е.Б. Перспективы развития спутниковых методов измерений в геодезии// Геодезия и картография. 1997. № 3. С. 11-13.

21. Кравчук И.М., Семеонова А.Е. Спутниковое нивелирование. Доклад на 64-й Юбилейной научно-технической конференции студентов, аспирантов и молодых ученых МИИГАиК, посвященной 230-й годовщине со дня его основания. М., 7-8 апреля 2009.

22. Клюшин Е.Б., Кравчук И.М. Спутниковое нивелирование. Сборник статей по итогам Международной научно-технической конференции, посвященной 230-летию основания МИИГАиК. Выпуск 2, часть II, МИИГАиК, М.: 2009.

23. Кравчук И.М. Особенности вычисления нормальных высот по результатам спутниковых измерений. Известия вузов. Геодезия и аэрофотосъемка. № 4, М.: МИИГАиК, 2010.

24. Крылов В.И. Космическая геодезия: Учебное пособие М.: Изд. МИИГАиК УПП «Репрография», 2002.

25. Куприянов А. О., Бородко Е. А. Комплексные испытания интегрированной картографо-геодезической спутниковой аппаратуры// Геодезия и картография, 2006, №10. С.41-45.

26. Логинов С.А., Бойков В.В., Рогозин В.П., Монахова М.А. Исследование точности Государственной геодезической сети на территории Московской области с использованием Спутниковой системы (проект «Москва»). // Геодезия и картография, 2006, №6. С.23-28.

27. Майоров А.Н. Разработка технологии и создания модели высот квазигеоида с использованием спутниковых данных. Диссертация кандидата технических наук. М.: 2008.

28. Методические указания. Нивелирование морских уровенных постов. Вып. 9, ч. 1. Л., Гидрометеоиздат, 1980, с. 48.

29. Михедов Е.Ю. GPS-нивелирование при коротких базах.- Донецк, 2005. -129 с.

30. Михедов Е.Ю. Исследование точности GPS-нивелирования при коротких векторах, http://masters.donntu.edu.ua/2004/ggeo/mikhedov/diss/index.

31. Наставление гидрологическим станциям и постам. Вып. 9, ч. 1- Л., Гидрометеоиздат, 1968, с. 424.

32. Непоклонов В. Б., Чугунов И.П., Яковенко П.Э., Орлов В.В. Новые возможности развития сети нормальных высот на территории России// Геодезия и картография. 1996. № 7. с. 20-22.

33. Огородова JI.B. Высшая геодезия. Часть III. Теоретическая геодезия: Учебник для вузов.- М.: Геодезкартиздат, 2006.- 384 с.

34. Основные положения о государственной геодезической сети Российской Федерации ГКИНП (ГНТА) 01-006-03. - М., Федеральная служба геодезии и картографии России, 2003 г.

35. Остроумов В.З., Шануров Г.А., Масленников A.B. Повышение точности определения высот уровенных постов. Геодезия и картография, 2003, № 11, с. 25-29.

36. Остроумов В.З., Шануров Г.А., Епишин В.И. Высотная основа уровенных постов: геодезический аспект. Геодезия и картография, 2005, № 4, с. 20-29.

37. Перов А.И., Харисов В.Н. ГЛОНАСС. Принципы построения и функционирования. М.: 2005.

38. Правила закладки центров на пунктах геодезической и нивелирной сетей. М., КГЦ Геодезиздат, 1993 г.

39. Правила закрепления центров на пунктах спутниковой геодезической сети. Утверждены Роскартографией как дополнение к ГКИНП-07-016-91 от 07.05.2001 г.

40. Правила по технике безопасности на топографо-геодезических работах. -М., Недра, 1991.

41. Проект «МОРЯ». Монографический справочник. «Гидрометеорология и гидрохимия морей». Т. VI, КАСПИЙСКОЕ МОРЕ. С.-Петербург, Гид-рометеоиздат, 1992.

42. РТМ. Применение приёмников спутниковой системы WILD GPS System 200 фирмы Лейка (Швейцария) при создании и реконструкции городских геодезических сетей. Нижний Новгород, Верхневолжское аэрогеодезическое предприятие, 1995.

43. Руководство по созданию и реконструкции городских геодезических сетей с использованием спутниковых систем ГЛОНАСС/GPS. — М.: ЦНИИГАиК, 2003.

44. Руководящий технический материал. «Спутниковая технология геодезических работ. Термины и определения». РТМ 68-14-01. Москва, ЦНИИГАиК, 2001.

45. Руководящий технический материал. Высотная привязка уровенных постов. -ГКИНП-03-215-88, М., ЦНИИГАиК, 1988, с. 41.

46. Рутлетдж Д., Хэндерсон К., Коэрнер Р., Ремонди Б. Мониторинг оседания нефтяного месторождения с помощью системы GPS. Перевод статьи опубликованной в журнале GPS World в 2002. http://www.agp.ru/projects/oilmon/index.htm.

47. Серебрякова Л.И., Горобец В.П., Козлова Л.Ю., Сермягин Р.А. К вопросу оценки точности спутниковых определений, выполняемых на геодинамических полигонах. // Геодезия и картография, 2006. №6. - С.34-39.

48. Соловьев Ю.А. Спутниковая навигация и ее приложения. М.: 2003.

49. Технические условия на опытно-методические работы по использованию GPS при высокоточном нивелировании./ В.И. Забнев. Роскартография, 2004.

50. Федосеев Ю.Е. Особенности применения спутникового нивелирования для построения городских высотных сетей. Автоматизированные технологии изысканий и проектирования. №2(17), 2005, С. 18.

51. Филлипов М.В., Янкуш А.Ю. Сравнение GPS- и традиционных геодезических работ//Геодезия и Картография.-1995.- №9. С.15-19.

52. Центр «Геодинамика». Аэропорт «ШЕРЕМЕТЬЕВО» г. Москва, http ://avia. geodinamika.ru/rus/proj ects/aeroportsheremetevogmoskva/#.

53. Шануров Г.А., Мельников С.Р. Геотроника. МИИГАиК, НПП «Геокосмос». -М., 2001, с. 136.

54. Ярмоленко А.С., Богомья А.Н. Оценка точности GPS-измерений. Геодезия и Картография.-1996.- № . С. 14-16.

55. Яшкин С.Н. Некоторые аспекты спутниковой градиентометрии. ГЕОПРОФИ, №; 2007.

56. Alan Н. Phillips. Geometrical Determination of PDOP. Navigation: Journal of The Institute of Navigation. Vol. 31, No. 4, Winter 1984-85, p. 329-337.

57. Assessment of height variations by GPS at Mediterranean and Black Sea coast tide gauges from the SELF projects. Global and Planetary Change, 34 (2002), 5-35.

58. Modern adjustment of geodetic network considering the elimination both of the systematic errors and the blunders-Surveying Science in Finland, 2002, Volume 14, no.1-2, p. 47-55.