Бесплатный автореферат и диссертация по наукам о земле на тему

Разработка и исследование методов, обеспечивающих повышение точности координатных определений в Социалистической Республике Вьетнам

ВАК РФ 25.00.32, Геодезия

Автореферат диссертации по теме "Разработка и исследование методов, обеспечивающих повышение точности координатных определений в Социалистической Республике Вьетнам"

005001888

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

НГУЕН ВАН ДОНГ

РАЗРАБОТКА И ИССЛЕДОВАНИЕ МЕТОДОВ, ОБЕСПЕЧИВАЮЩИХ ПОВЫШЕНИЕ ТОЧНОСТИ КООРДИНАТНЫХ ОПРЕДЕЛЕНИЙ В СОЦИАЛИСТИЧЕСКОЙ РЕСПУБЛИКЕ ВЬЕТНАМ

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

АВТОРЕФЕРАТ

диссертации на соискание ученой степени

кандидата технических наук 2 4 НОЯ 2011

МОСКВА-2011

005001888

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

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

Краснорылов Игорь Ильич

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

Татевян Сурия Керимовна I доктор технических наук, профессор

Шануров Геннадий Анатольевич

Ведущая организация: НИЦ (Топогеодезического и навигационного обеспечения) ФБУ «27 ЦНИИ Минобороны»

Защита диссертации состоится « б » с^оклЬря 2011 г. в «-/2» часов ш заседании диссертационного совета Д.212.143.03 в Московского государственном университете геодезии и картографии (МИИГАиК) по адресу 105064, Москва, Гороховский переулок, дом 4, МИИГАиК, зал заседани диссертационного совета.

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

Автореферат разослан « Ц- » Ноября 2011 г.

Ученый секретарь диссертационного совета

Климков

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ Актуальность темы диссертации. Современные требования к точности координатных определений очень высокие для решения новых задач, например, мониторинга смещений и деформаций сооружений, особенно имеющих большую высоту; изучения перемещений литосферных плит (блоков); кадастра и учета недвижимости; установки и эксплуатации буровых платформ на море; обеспечения безопасности различных видов транспорта (железнодорожное движение (особенно скоростное), воздушные полеты). Требуемые точности координатных определений уже менее .одного сантиметра. Однако с течением времени опорная геодезическая сеть постепенно «стареет», теряя свою первоначальную точность из-за влияния геодинамических эффектов. Результаты последних геодинамических исследований показывают, что изменение координат пунктов, находящихся на территории Социалистической Республики Вьетнам (СРВ), составляет около 50 мм/год, а изменение взаимного расположения пунктов, находящихся на севере и юге СРВ, составляет 4... 25 мм/год. Поэтому для повышения точности координатных определений, прежде всего, необходимо решить задачу по выявлению и учету геодинамических факторов.

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

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

влияние тропосферы. К настоящему времени для учета влияния тропосферы можно применять модели, позволяющие оценить величину тропосферных задержек (ТЗ): Хопфилд; Саастамойнена; Нейлла; MOPS (Minimum Operational Performance Standards); GCAT (GPS Code Analysis Tool); Университета Ныо Брунсвика (УНБ); Блэка. При построении этих моделей особенности условий Вьетнама не были учтены. Поэтому возникает вопрос о том, какая из моделей тропосферы дает наиболее точные ТЗ в соответствии с условиями Вьетнама. При учете влияния тропосферы также можно использовать значения зенитных тропосферных задержек (ЗТЗ), определенных методом точного позиционирования. В связи с этим для повышения точности координатных определений требуется определение наиболее эффективного способа учета влияния тропосферы в климатических условиях Вьетнама.

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

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

• Анализ современного состояния координатного обеспечения СРВ. На основе анализа особенностей координатного обеспечения СРВ сделан вывод о недостатках государственной геодезической сети (ГГС) СРВ и предложены основные направления для повышения точности координатных определений в СРВ.

* Разработка и исследование методов решения задач геодинамики.

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

оптимизация двух геодинамических сетей (Дьснбьсн-Лайчау и Сонла) в СРВ. Сформулированы рекомендации по практическому применению выполненных разработок.

• Анализ способов учета влняння тропосферы. Рассмотрены модели тропосферы: Хопфилд, Саастамойнена, Нейлла, MOPS, GCAT, УНБ, Блэка. Исследована возможность применения метода точного позиционирования для определения ЗТЗ. Рассмотрена эффективность использования измеренных метеопараметров при учете влияния тропосферы. Исследованы способы учета влияния тропосферы при спутниковых измерениях. Разработаны соответствующие рекомендации.

Научная новизна н результаты работы, выносимые на защиту:

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

• Метод оптимизации геодинамических сетей по критерию минимальной стоимости реализации проекта с использованием генетического алгоритма. Результаты оптимизации двух геодинамических сетей Дьенбьен-Лайчау и Сонла на территории СРВ с использованием этого метода.

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

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

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

Апробация работы. Основные результаты по теме диссертации докладывались на конференциях студентов, аспирантов и молодых ученых МИИГАиК 2009, 2010 и 2011 гг., на Международной научно-технической конференции «Геодезия, картография и кадастр - XXI век», посвященной 230-летию основания Московского государственного университета геодезии и картографии (25-27 мая 2009 г.).

Публикации. Основное содержание диссертации опубликовано в четырех работах в журнале Известия вузов «Геодезия и аэрофотосъемка», включенном в перечень ВАК.

Соответствие диссертации Паспорту научной специальности. Диссертация соответствует специальности 25.00.32 «Геодезия». Диссертационное исследование выполнено в соответствии с пунктом 3 паспорта специальности ВАК - «Геодезические (глобальные) навигационные спутниковые системы и технологии....»; пунктом 8 - «Геодезический мониторинг напряженно-деформированного состояния земной коры,...».

Структура и объем работы. Диссертация состоит из введения, трех глав, заключения, списка сокращенных слов и списка использованных источников. Основное содержание изложено на 171 странице текста. Работа содержит 40 рисунков и 14 таблиц. Список литературы включает в себя 76 наименований, в том числе 32 на русском языке, 44 на иностранных языках.

СОДЕРЖАНИЕ РАБОТЫ Введение содержит обоснование актуальности темы диссертации, формулировки цели, задач работы и ожидаемых научных результатов, которые автор предполагает получить при разработке диссертации.

В первой главе «Анализ современного состояния координатного обеспечения СРВ» приведены сведения из истории построения координатного обеспечения Вьетнама и его современное состояние. В результате выделены два направления исследования в соответствии с существующими требованиями геодезического обеспечения во Вьетнаме, возникающими при решении научно-практических задач: изучение геодинамических влияний; анализ способов учета влияния тропосферы при спутниковых наблюдениях в условиях СРВ.

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

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

х,.(0 = х('+х1м,

. у;.(г) = у(° + у1д/, (1)

2,(/) = г,°+2,.Д/,

где Г" 2,°) , У. - координаты и составляющие скорости

движения пункта /' в начальную эпоху /п; Д/ = / - /0 - время, отсчитываемое от начальной эпохи. Выражение компонентов векторов базовых линий между пунктами / и / в эпоху I имеет вид

дх, (О = Х){1)-ХХЧ

дг,(/)=гу(0-^(0, (2)

[дг,со = ¿,(0-2,(0-

Подставляя (1) в (2), получим

=х) + - -ДУД*) = Г/ + У,Д/ - У,0 - (3)

дг. (о = г)+¿, д/ - г,0 - ¿,д/.

На основе выражении (3), приведем уравнение поправок в виде

АХ + Ь = У, Рь, (4)

где А - матрица коэффициентов. Ее фрагмент для одной линии имеет вид:

"-10 О -А( 0 0 1 0 0 Дг 0 0

А= 0 -1 0 0 -А/ О 0 1 0 0 ДГ 0 ; (5)

0 0 -1 0 0 -Дг 0 0 1 о о м X - неизвестные, поправки к координатам пунктов и составляющим скоростей их движения

Х = (АХ* Д У? Д2,° АХ: А У, АЙ, АХ° А У/ Ы) АХ1. ДУ,. Д2;)Г; (6) Ь - вектор свободных членов

L =

AXUO-AXit)

(7)

ДУ,°(0-ДУ,(0 ^(O-ÄZ,«,

где АХЦ(/), АУ°(/), AZ°(t) - вычисленные (счислимые) компоненты вектора

базовых линий на эпоху 1\ АХЦ(1),АУ^(1),А2Ц(1.) - «измеренные» компоненты

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

К недостатку предлагаемого метода можно отнести существенное увеличение размерности матриц, требующихся для его реализации. Для преодоления этого недостатка мы предлагаем использовать особые матричные преобразования. Поскольку у матрицы коэффициентов уравнения поправок А и матрицы весов уравнений поправок Рг существует большое число нулевых элементов, для уменьшения объема вычислительных операций нами предлагается преобразование матриц к виду векторов. Составляем один вектор для сохранения значений ненулевых элементов и два вектора для индекса (номер строки и столбца) соответствующих элементов. Необходимое время для совместного уравнивания координат и составляющих скоростей движений пунктов при разном количестве измеренных величин с использованием преобразования матриц к виду векторов и без его использования изображено на рис. 1. Расчеты выполнены на персональном компьютере с характеристиками: процессор Intel Core2 Duo Т6600 2,2 ГГц; объем памяти 4096 Мб; операционная система Windows 7 Домашняя расширенная.

500

— 300 я

% 200 S

100 о

С преобразованием Без преобразования

0 5000 10000 15000 20000

Число измеренных величин

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

: : тт^тМ; '

£3 «зжаатай »«кг

ршо 8М?шАМ «»г»»»«

Рис. 2. Движения станций, находящихся вблизи территории Вьетнама

оч&жтхмл' Мтяг

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

2. Оптимизация геодинамических СР8/ГЛОНАСС сетей. Рассматривался вопрос об оптимизации существующих сетей по критерию минимальной затраты для реализации сети. Целью оптимизации, прежде всего, является оптимизация весов измерений. Это приводит к многокритериальной задаче

ф= 2]с,Р<

а. <а.

<о<

Р, ь Р,.

1 = 1.../, 1=1...г.

(8)

где Ф - целевая функция; с, - коэффициенты относительной выгодности, которые определяются по формуле

затраты на 1-ое измерение затраты на 1 -ое измерение

(9)

р, - оптимизируемые параметры (веса измерений); р1м„„, р^лшк - минимально и максимально достижимое значение весов измерений базовых линий; а„ <т, мак -средние квадратические ошибки деформационных параметров и их предельное значение; / - число наблюдений; г - число деформационных параметров.

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

где С, - стоимость измерения у-го пункта; <т„ а1ек - средние квадратические ошибки деформационных параметров до и после включения пункта с номером Г, X], У], X/, мин, У], мии, мак, Ум мак - координаты пункта и их предел.

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

Нами было разработано программно-алгоритмическое обеспечение для оптимального проектирования геодинамических сетей. С его использованием мы провели оптимизацию двух геодинамических сетей во Вьетнаме: Дьенбьен-Лайчау и Сонла. При оптимизации было принято: средняя квадратическая ошибка единицы веса ао = 1; коэффициенты выгодности всех измерений при оптимизации с точки зрения весов измерений с, = 1, V/; стоимость измерения пункта одинакова при улучшении сети, т. е. Су = С2 = ... = С„, где и - число пунктов проектируемой сети; точность определения плановых компонент

(10)

У),мш< ~ У у ~ У\

сгах1 = стд г, =o'i=a + bLi=(3 + Q, )мм;

(п)

максимальные и минимальные веса

р:™ = Р": = рг=О, 02)

(7,

Геодинамическая сеть Дьенбьен-Лайчау построена для исследования деформаций земной коры, вблизи разлома Дьенбьен-Лайчау и состоит из 7 пунктов (рис. 3). Рассмотрим задачу оптимизации, в которой модель деформации имеет вид

к=0, у,=0, (13)

К=яо. ^¡¡=к

где м„ V, - составляющие движений пунктов, находящихся на блоке /'; а0, Ь„ -параметры, выражающие двюкения блока II относительно блока I.

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

Таблица 1. Результаты оптимизации весов измерений геодинамической сети Дьенбьен-Лайчау

Базовые линии Предел ошибок деформационных параметров а"™ = а"°к =

3,0 им/год 2,5 мм/год 2,0 мм/год 1,5 мм/год

КСА1-В01Ч1 0,90 0,84 0,96 1

НАМ1-ЬЕМ1 0.84 0,86 0.88 1

00141-НАМ 1 0,90 0,93 1

ТАШ-ТАТ1 0,98 1

ТАШ-ТАТ1 0,63 1

ЬЕМ1-ТАШ 1

ЬЕМЬТАШ 1

ША1-ЬЕМ1 1

БОМЬТАШ 0,99

ООЖ -ТАШ 0,88

Для исследования деформаций внутри литосферных блоков в модель, имеющую вид (13), добавлены три деформационных параметра. Модель деформаций района Дьенбьен-Лайчау имеет вид

Г», = еЛ + гхуу1, V, = е^х, + еуу,, }

\и„ =а0 + ехх„ + 8„,у„, \>„ =Ь0 + е„хп + е):у„, где «„ V, и X/, у, - составляющие движений и координаты пунктов, находящихся

на блоке /; а0, Ь0 - параметры, выражающие движения блока II относительно блока I; ех, еу, е„ - параметры, выражающие деформации внутри блоков. Проведенные расчеты показывают, что в этом случае построенная геодинамическая сеть обеспечивает определение деформационных параметров ао, Ьо, £х, £у, е-у с точностью 1,7 мм/год, 5,4 мм/год, 2,8-10"8 год"1, 31,5-10"8 год"1, 5,8-10"8 год"1 соответственно. Для определения деформационных параметров с более высокой точностью требуется включение в состав сети дополнительных пунктов. Пусть допустимая область изображена на рис. 3. С использованием разработанного программно-алгоритмического обеспечения, мы определили положения восьми наиболее информативных дополнительных пунктов РОТ1, POI2,..., POI8 (рис. 3). Ошибки деформационных параметров, полученные при включении в сеть дополнительных пунктов, показаны в таблице 2.

22« 30' -

УОГ5. Р»0|РО15 PÔI4 Is'GAl

DON]

«Ml

¿Ш.М1

•РОИ

21* 10'

№2» 55'

■ im-

pôt:

Рис. 3. Наиболее информативные дополнительные пункты геодинамической сети Дьенбьен-Лайчау

РОВ

Обозначения -—- Разлом

-Граница допустимой области

^ Сушесгвугагшк* тта: » Дополнительные пункт

при

103'05' ¡03* 15' 103« 25' Таблица 2. Априорные ошибки деформационных параметров включении в состав сети Дьенбьен-Лайчау дополнительных пунктов

Число БК. допол. пунктов с„о 1 оьо вех 1 Osv 1 Сот Число вк. допол. пунктов 1 Обо a«r 1 о£1, | ает

мм/год 10"8 год"1 мм/год 10"8год"'

1 1,4 2,2 2.3 8,6 2,3 5 0.9 1,1 1,5 4,4 1,4

2 1.2 1,6 2.1 6.5 2.1 6 0.8 0,9 1,4 4.2 1,3

3 1,2 1,5 1.9 4,9 1.8 7 0.7 0,8 1.3 4.0 1,2

4 1,0 1,2 1.7 4,6 1,6 8 0.7 0,8 1,2 3,6 1,2

Геодинамическая сеть Сонла создана для мониторинга деформаций земной коры из-за действий разломов Сонла и Черной реки и состоит из 7 пунктов (рис. 4). Рассмотрим задачу оптимизации, в которой модель деформаций имеет вид:

и, = а0 + ЕхХ, + Е у1, ш ~

У11 =5^11 +еуУ,П

(15)

"ш =с»+£ а + «л, = ¿о + ад» + еуУш. где V,, и X/,- составляющие двгокений и координаты пунктов, находящихся

на блоке /; яо, ¿'о, <-'о, "о - параметры, выражающие движения блоков I и III

относительно блока II; еЛ £у, сху - параметры, выражающие деформации внутри

блоков. Была проведена оптимизация для 9 вариантов пределов ошибок

деформационных параметров с принятием о™' = сг"°' = а"" = а^" = и

а"" = о"°к = = а"""'. В таблице 3 представлена оптимальная доля

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

Таблица 3. Результаты оптимизации весов измерений геодинамической сети Сонла

_мак . <7Ц , мм/год 3,5 3,5 3,5 3,0 3,0 3,0 2,7 2,7 2,7

10"7 год'1 2,00 1,50 1,25 2,00 1,50 1,25 2,00 1,50 1,25

ШТ1-ЫА02 1 1 1 1 1 1 1 1 1

мож-моп 1 1 1 1 1 1 1 1 1

ЫА02-ИАН1 1 1 1 1 1 1 1 1 1

]\Ю11-ЫАН1 1 1 1 1 1 1 1 1 1

ТРШ-ШТ1 1 1 1 1 1 1 1 1 1

ТРШ-<ЗТА2 1 1 1 1 1 1 1 1 1

моьп-ИАог 0,99 1 1 1 1 1 1 1 1

я МОШ-ТРШ 0,63 0.72 0,99 0,92 0.92 0.99 1 1 0,97

Е 0ТА2-КАО2 1 1 1 1 1 1 1 1

Ц ТРШ-ЫАШ 1 1 1 1 1 1 1 1

МОЖ-ШЛ 1 1 1 1 1 1 1

о а ьа N011 -ЫАР2 I 1 1 1 1 1 1

ЬОТ1-КАН1 1 1 1 1 1

Ь0Т1-М011 1 1 1 1 1

МОШ-ЫАН! 1 1 1 1 1

МОМ1-0ТА2 1 1 1 1 1

КОП-(}ТА2 1 1 1 1 1

три 1-N011 1 1 1 1 1

0ТА2-КАН1 1 1 1

ТРШ^А02 1 1 1

Расчет, проведенный нами, свидетельствует, что наблюдения на всех пунктах сети дает точность определения деформационных параметров а0, Ь0, с0, с1о на уровне 2,7 мм/год, а параметров ех, еу, 1,23 Ю"7 год"1. Для получения

более высокой точности необходимо перепроектировать сеть с добавлением новых пунктов. Пусть дополнительные пункты должны находиться в допустимой области, которая изображена на рис. 4. С использованием разработанного программно-алгоритмического обеспечения мы определили положения восьми наиболее информативных дополнительных пунктов РОИ, POI2,..., POI8 (рис. 4). Ошибки деформационных параметров, полученные при включении в сеть дополнительных пунктов, показаны в таблице 4.

Рис. 4. Наиболее информативные дополнительные пункты геодинамической сети Сонла Обозначения ---Разлом

-Транша допустимой области

Существующие пункты « Дополнительные пункты

И' 05'

104» 2V

104» Ю 10-1= 10' Долгота

Таблица 4. Априорные ошибки деформационных параметров при включении в состав сети Сонла дополнительных пунктов

Число вк. допол. пунктов мм/год мм/год о;». мм/год О^о. мм/год Ю-8 год"1 Ю-8 год"' 10"8год"'

1 1,3 1,4 1,2 1,2 3,4 4,9 2,8

2 1,1 1,2 1,1 1,2 3,0 4,4 2,5

3 0,9 0,9 1,0 1,0 2,7 3,6 2.2

4 0,8 0,8 0.9 0.9 2,6 3,3 2,1

5 0,8 0,8 0,9 0,9 2,2 3,0 1,8

6 0,7 0,7 0,8 0,8 2,1 2,7 1,7

7 0,7 0,7 0,7 0,7 2,0 2,5 1,6

8 0,6 0,6 0,7 0,7 1,8 2,3 1,5

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

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

1. Определение ЗТЗ методом точного позиционирования и анализ моделей тропосферы. Приведем уравнения наблюдений спутника к на пункте i при наблюдениях двухчастотными приемниками с использованием данных Международной ГНСС службы (МГС)

(]6)

К, =pf+КК, + КЛ„+cdt,+4К-

где Pl,Q>l - псевдодальность и фаза комбинационной волны Ls, выраженные в единицах расстояния; р* - геометрическое расстояние пункта /' до спутника к; с - скорость света в вакууме; dt, - ошибки часов приемника; Лз - длина комбинационной волны L3, N*,- число неоднозначности волны L3; Dzci,Dzejtl -сухая и влажная компоненты ЗТЗ; - функции отображения.

Линеаризуя выражения (16), получим уравнения поправок. Решая уравнения поправок методом наименьших квадратов, получим параметры тропосферы, ошибки часов приемника, и положения пункта.

Определяя ЗТЗ на пунктах, для оценки эффективности моделей тропосферы нами введён показатель остатка а, выраженный в процентах

а = 1^|хюо, (17)

D-.

где ДО./', Dj - остаточная и уточенная величины ЗТЗ, полученные из обработки результатов наблюдений методом точного позиционирования.

С использованием метода точного позиционирования, мы определили ЗТЗ по результатам измерений на двадцати станциях, в том числе три станции, находящиеся на территории Вьетнама; двенадцать станций - в сходных с Вьетнамом условиях; пять станций в США в условиях, отличающихся от Вьетнама. ЗТЗ на станциях определены с использованием модельных

метеопараметров. Для четырех станций ОАЯ\У, ТСМБ, ТЫМЦ и ТХУТИ, на которых производились метеорологические измерения, расчеты выполнены еще с использованием измеренных метеоданных по моделям Блэка, Хопфилд, УНБ и Саастамойнена, вычислительные алгоритмы которых явно зависят от метеопараметров. Мы получили ЗТЗ с точностью на уровне 6 мм, и определили значения показателей остатка а, которые приведены в таблицах 5-6. Таблица 5. Показатель остатка а (%) при использовании метеопараметров, определяющихся

по модели Международной стандартной атмосферы ISA

№ Станция (страна) Модель тропосферы

Блэка GCAT Хопф. MOPS Нейлла УНБ Саас.

1 BLVI (Вьетнам)1 6 6 6 2 6 9 6

2 CODA (Вьетнам) 6 7 7 1 7 10 7

3 DANA (Вьетнам) 8 8 8 2 9 11 9

4 ВАКО (Индонезия) 9 9 10 3 10 12 10

5 CUSV (Таиланд) 9 9 9 2 9 12 9

6 DARW (Австралия) 7 7 8 2 8 11 8

1 HAV2 (Индия) 7 8 8 1 8 11 8

8 HUTB (Индия) 7 8 8 1 8 11 8

9 KUNM (Китай) 10 5 10 5 6 11 10

10 NTUS (Сингапур) 8 8 9 2 9 12 9

11 PIMO (Филиппины) 10 10 10 3 10 13 10

12 SAMP (Индонезия) 8 9 9 2 9 12 9

13 TCMS (Тайвань) 7 7 7 1 8 11 8

14 TNML (Тайвань) 7 7 8 1 8 11 8

15 TWTF (Тайвань) 7 7 8 1 7 11 8

16 PI67 (США) 2 1 2 3 1 4 3

17 PI68 (США) 2 1 2 3 1 4 3

18 PI 69 (США) 2 1 2 3 1 4 3

19 PI 70 (США) 2 1 2 3 1 4 3

20 Р326 (США) 2 1 3 2 1 4 3

- наблюдения на станциях BLVI, CODA, DANA производились в период с 22 по 24 марта 2009 г., а наблюдения на остальных станциях в период с 04 по 10 июля 2010 г.

Таблица 6. Показатель остатка а (%) при использовании

измеренных метеопараметров в период с 04 по 10 июля 2010 года

№ Станция (страна) Модель тропосферы

Блэка Хопфилд УНБ Саастамойнена

1 DARW (Австралия) 1 2 7 2

2 TCMS (Тайвань) 2 1 6 1

3 TNML (Тайвань) 2 1 6 1

4 TWTF (Тайвань) 1 1 7 1

Наше исследование подтверждает, что для учета влияния тропосферы в

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

MOPS оказывается наиболее эффективной. Она позволяет учитывать в среднем 98% влияния тропосферы, а с использованием оставшихся моделей учтено в среднем 91% влияния тропосферы. Если производились метеорологические измерения, использование моделей Блэка, Хопфилд и Саастамойнена будет эффективным. При этом в среднем учитывается 99% влияния тропосферы. Также видно, что применение моделей тропосферы в условиях США, особенности территории которых учитывались при построении моделей, эффективнее чем в условиях Вьетнама.

2. О точности определения координат пунктов в абсолютном методе в зависимости от способа учета тропосферных задержек. Нами определены координаты станций МГС, находящейся в сходных с СРВ условиях, по результатам спутниковых наблюдений 08 июля 2010 г. Расчеты проведены для случаев: использование моделей тропосферы для учета влияния тропосферы; без учета тропосферы; с уточнением тропосферы в процессе обработки. В качестве метеопараметров использованы метеопараметры, определяющиеся по модели Международной стандартной атмосферы. Также проведены расчеты в случае использования измеренных метеопараметров с использованием модели тропосферы Блэка, Хопфилд, УНБ и Саастамойнена для учета тропосферы. Для оценивания эффективности модели тропосферы в абсолютном методе приняты расхождения положений станций 5г, полученных нами и МГС. Значения расхождений положений в зависимости от способов учета влияния тропосферы приведены на рис. 5.

Из проведенных расчетов можно сделать следующие выводы. В абсолютном методе, при отсутствии измеренных метеопараметров учет влияния тропосферы по модели MOPS позволяет определить координаты с точностью лучше, чем при учете влияния тропосферы по остальным моделям. При присутствии измеренных метеопараметров, модели Блэка, Хопфилд, Саастамойнена наиболее эффективны для учета влияния тропосферы. При учете влияния тропосферы моделями использование измеренных метеопараметров повысит точность координатных определений в 2 -н 10 раз.

Станция

0,6

' • Блэка

Хопфилд — УНБ ' Саасгамойнена С уточнением

а' Рис. 5. Расхождения вычисленных положений с решением МГС в зависимости от способов учета влияния

тропосферы

а) при использовании модельных метеопараметров

б) при использовании измеренных метеопараметров

^ # / Станции

Для получения более точных результатов определения координат при обработке необходимо вводить ЗТЗ как неизвестные и уточнять их. При этом в данной работе расхождения положений станции 8г, полученные нами и МГС, составляют ~ 3 см. Такие расхождения также получены в обоих случаях (присутствия и отсутствия измеренных метеопараметров).

3. Учет влияния тропосферы при относительных наблюдениях. Нами определены координаты пунктов двух сетей (рис. 6) по результатам спутниковых наблюдений 04 июля 2010 г., с использованием моделей тропосферы без измеренных метеопараметров и значений ЗТЗ, определенных методом точного позиционирования для учета влияния тропосферы. Расчеты также проведены без учета влияния тропосферы. В таблицах 7 - 8 приведены средние квадратические ошибки (СКО) определения координат в зависимости от способов учета влияния тропосферы и СКО в положении, которые определяются по формуле

тг = ^т\ +т]+т\, (18)

где ту, ту, т2 - ошибки в координатах. Здесь пункт ТАУТР сети в Тайване и

пункт р 168 сети в США приняты за опорные пункты.

Долгота Долгота

а) Сеть в Тайване б) Сеть в США

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

Способ учета влияния тропосферы Пункты

TCMS TNML

тх ту mz т, ту тг /и/ тг

Блэка 8.1 7,0 3.1 10,8 8,0 6,8 3.0 10,6

GCAT 1,4 1,2 0.5 1,9 1,3 1,1 0,5 1,9

Хопфилд 8,1 7,0 J 3.1 10,8 8,0 6,8 3.0 10,7

MOPS 4.0 3,5 1.5 5.5 4.0 3,4 1.5 5.4

Нейлла 1.4 1,2 0,5 1,9 1.4 1,2 0,5 1.9

УНБ 4,1 3.6 1.6 5,6 4,1 3,5 1,5 5,5

Саастамойнена 8,2 7,1 3.2 11,0 8,1 6.9 3.0 10,8

Без учета 13,7 11,7 5,2 18,2 13,5 11,5 5Д 17,9

ЗТЗррр 1.2 1,1 0,5 1,8 1,2 1,1 0,5 1,8

Таблица 8. Ошибки определения координат пунктов в сети США

при использовании разных способов учета влияния тропосферы, мм

Способ учета влияния троп. Пункты

pi 67 р169 р 170 Р326

тх ту тг тг тх ту mz тг тх ту mz т, тх ту т7. тг

Блэка 7,9 6,2 4,1 10,9 10,2 8.0 6,1 14,4 8,8 7,0 5,2 12,3 8,3 7,1 4,8 11,9

GCAT 2,0 1,6 1Д 2,8 2,6 2,0 1,6 3,7 2,2 1,8 1,3 3,1 2,1 1,8 1,2 3,0

Хопф. 7,9 6,2 4,1 10,8 10,2 8,0 6,1 14,3 8,8 6,9 5,2 12,3 8,3 7,1 4,8 11,9

MOPS 1,5 1,2 0,8 2,0 1,9 1,5 1,2 2,7 1,7 1,3 1,0 2,3 1,6 1,3 0,9 2,2

Нейлла 2,0 1,6 1,1 2,8 2,6 2Д 1,6 3,7 2,3 1,8 1,3 3,2 2,1 1,8 1.2 3,1

УНБ 1,6 1,2 0,8 2,1 2,0 1,6 1,2 2,8 1,7 1,4 1,0 2,4 1,6 1,4 0,9 2,3

Саас. 7,9 6,2 4,1 10,8 10,2 8,0 6,1 14,3 8,8 7,0 5,2 12,3 8,3 7,1 4,8 11,9

Без учета 5,9 4,6 3,1 8,1 7,6 6,0 4.6 10,7 6,5 5,2 3,9 9,2 6,2 5,3 3,6 8,9

ЗТЗРРР 1,0 0,8 0,5 1,4 1,3 1,0 0,8 1,8 1,1 0,9 0,7 1,6 1.1 0,9 0,6 1,5

Из полученных результатов следует вывод: в рассматриваемых сетях без учета влияния тропосферы мы получили СКО по каждой координате 1,5 см, СКО в положении пунктов 2 см. При учете влияния тропосферы моделями СКО по каждой координате составляют 3...11 мм в зависимости от применяемой модели, а СКО в положении пунктов - 4... 15 мм. При учете влияния тропосферы с использованием ЗТЗ, определенных методом точного позиционирования, СКО по каждой координате достигают 1,5 мм, а СКО в положении пунктов - 2 мм.

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

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

Реализован в среде Visual С++ алгоритм метода совместного уравнивания координат и составляющих скоростей движений пунктов. Разработанная программа апробирована на сети действующих станциях, находящихся вблизи территории Вьетнама. Полученная точность определения координат и составляющих скоростей движений станций этой сети предлагаемым методом составляет 7 + 17 мм и 2,0 + 5,4 мм/год соответственно. Точность горизонтальных составляющих скоростей составляет 1,1 1,7 мм/год, а вертикальной составляющей 4,4 + 6,6 мм/год. Точность определения положения

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

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

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

3. Выполнено исследование определения ЗТЗ методом точного позиционирования. Реализован алгоритм этого метода в среде Visual С++ с использованием библиотеки с открытым кодом. Проведено уточнение ЗТЗ при спутниковых наблюдениях, которые заранее определены с использованием моделей Блэка, Хопфилд, GCAT, MOPS, Нейлла, УНБ и Саастамойнена. В результате получены ЗТЗ в сходных с Вьетнамом условиях с точностью на уровне 3 + 14 мм, а в среднем 6 мм.

4. Проведена оценка эффективности моделей тропосферы Блэка, Хопфилд, GCAT, MOPS, Нейлла, УНБ и Саастамойнена. Полученные результаты свидетельствуют о том, что для учета влияния тропосферы в сходных с Вьетнамом условиях при использовании модельных метеопараметров модель MOPS оказывается наиболее эффективной. Она позволяет учитывать в среднем 98% влияния тропосферы, причем с

использованием других моделей учтено в среднем 91% влияния тропосферы. Если производить метеорологические измерения в момент наблюдений, использование моделей Блэка, Хопфилд и Саастамойнена будет эффективным. С их использованием учтено 99% влияния тропосферы.

5. Выполнено исследование точности координатных определений абсолютным методом в зависимости от способа учета влияния тропосферы. Из полученных результатов следует, что при определении координат абсолютным методом необходимо вводить влияние тропосферы, как дополнительное неизвестное, и уточнять его в процессе обработки. При этом в настоящей работе расхождения положений пунктов, полученных нами, с решением МГС составляют 3 см при использовании и модельных, и измеренных метеопараметров для учета влияния тропосферы. В случае, когда нет возможности уточнения ЗТЗ, для получения наивысшей точности определения координат при отсутствии измеренных метеопараметров предлагается использование модели тропосферы MOPS, а при присутствии измеренных метеопараметров, следует использовать модели Блэка, Хопфилд или Саастамойнена. При учете влияния тропосферы моделями использование измеренных метеопараметров повысит точность координатных определений в 2 + 10 раз.

6. Выполнено исследование точности координатных определений относительным методом в зависимости от способа учета влияния тропосферы. Разработана подпрограмма, включение которой в программу обработки результатов спутниковых измерений DDBase позволит учитывать влияние тропосферы на относительные наблюдения путем использования ЗТЗ, полученных методом точного позиционирования. Это позволяет получить наивысшую точность координатных определений. Согласно полученным в работе данным, при расстоянии между пунктами 10 40 км и использовании ЗТЗ, полученных методом точного позиционирования, ошибки в положении пунктов достигают 2 мм. Если в относительном методе не учтено влияние тропосферы, эти ошибки составляют 2 см, а с использованием моделей для

учета влияния тропосферы эти ошибки составляют 4... 15 мм в зависимости от применяемой модели.

Согласно имеющимся у нас данным есть основания полагать, что полученные в диссертации результаты будут учтены при модернизации координатного обеспечения СРВ. В частности, будет осуществлена привязка ГТС СРВ к геодинамической сети МГС. Планируется также продолжение нами проведенных в диссертации исследований. Оно будет касаться обеспечения геодинамических исследований во Вьетнаме точной эфемеридной информации спутников GPS/ГЛОНАСС. Основы для проведения такой работы заложены в нашей статье [1].

Основные результаты исследований по теме диссертации опубликованы в журнале Известия вузов «Геодезия и аэрофотосъемка»:

1. Краснорылов И.И., Нгуен Ван Донг, Сорокин Н. А. Исследование возмущенного движения космического аппарата в случае его запуска на орбиту Земли в треугольную точку либрации L4 системы Солнце-Земля // Известия вузов. «Геодезия и аэрофотосъемка». 2009. № 5. С. 69-73.

2. Нгуен Ван Донг. Повышение точности Государственной геодезической сети Вьетнама // Известия вузов. «Геодезия и аэрофотосъемка». - 2010. № 3. С. 19-21.

3. Нгуен Ван Донг. Определение координат пунктов и их изменений из наблюдений спутников, образующих глобальные навигационные спутниковые системы (на примере одного из районов Юго-Восточной Азии) // Известия вузов. «Геодезия и аэрофотосъемка». 2011. № 4. С. 30-35.

4. Нгуен Ван Донг. Анализ моделей тропосферы по результатам спутниковых наблюдений в регионе Юго-Восточной Азии // Известия вузов. «Геодезия и аэрофотосъемка». 2011. № 5. С. 21 - 25.

Подписано в печать 02.11.2011. Гарнитура Тайме Формат 607 90/16. Бумага офсетная. Печать офсетная.

Объем 1,5 усл. печ. л.

Тираж 80 экз. Заказ №227 Цена договорная Издательство МИИГАиК 105064, Москва, Гороховский пер .,4

Содержание диссертации, кандидата технических наук, Нгуен Ван Донг

ВВЕДЕНИЕ.

1. АНАЛИЗ СОВРЕМЕННОГО' СОСТОЯНИЯ КООРДИНАТНОГО ОБЕСПЕЧЕНИЯ СОЦИАЛИСТИЧЕСКОЙ РЕСПУБЛИКИ ВЬЕТНАМ.

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

1.2'. Современные геодезические построения в СРВ.

1.3. О соответствии с существующими требованиями геодезического обеспечения в СРВ, возникающими при решении научных и производственных задач.:.

1.3.1. Изучение геодинамических процессов.

1.3.2. Анализ способов учета влияния тропосферы'и-в перспективе построение моделей тропосферы, соответствующих условиям Вьетнама.

1.4. Выводы по главе 1.

2. РАЗРАБОТКА И ИССЛЕДОВАНИЕ МЕТОДОВ РЕШЕНИЯ ЗАДАЧ ГЕОДИНАМИКИ.29'

2.1. Постановка задачи исследований.

2.2. Особенности геологической структуры СРВ.

2.3. Совершенствование теоретических основ методов изучения геодинамических влияний.

2.3.Ь. Алгоритм метода совместного уравнивания координат и составляющих скоростей движений пунктов.

2.3.2. Определение матрицы весов уравнений поправок.

2.3.3. Преобразование матриц к виду векторов.

2.3.4. Определение скорости движения пунктов, находящихся на сопредельной с Вьетнамом территории.

2.4. Оптимизация геодинамических ОР8/ГЛОНАСС сетей.

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

2.4.2. Методы оптимизации геодинамических сетей.

2.4.2.1. Априорная оценка точности деформационных параметров проектируемой сети.

2.4.2.2. Генетический алгоритм оптимизации.

2.4.2.3. К определению расположения пунктов относительно литосферного блока.

2.4.2.4. Алгоритм оптимизации весов измерений.

2.4.2.5. Алгоритм улучшения геодинамических сетей.

2.4.3. Примеры оптимального проектирования геодинамических GPS/ГЛОНАСС сетей.

2.4.3.1. Оптимизация геодинамической сети Дьенбьен-Лайчау.

2.4.3.2. Оптимизация геодинамической сети Сонла.

2.5. Выводы по главе 2.

3. О ВЛИЯНИИ ТРОПОСФЕРЫ ПРИ СОЗДАНИИ СПУТНИКОВЫХ ПОСТРОЕНИЙ В СРВ.

З.Г. Постановка задачи и общие замечания.

3.2. Принятые способы учета влияния тропосферы.

3.2.1. Модель Хопфилд.

3.2.2. Модель Саастамойнена.

3.2.3. Модель Нейлла.

3.2.4. Модель MOPS (Minimum Operational Performance Standards)

3.2.5. Модель GCAT (GPS Code Analysis Tool).

3.2.6. Модель Университета Нью Брунсвика.

3.2.7. Модель Блэка.

3.3. Особенности метеорологических условий в СРВ.

3.4. Априорная оценка тропосферных задержек на спутниковые измерения в СРВ.

3.5. Зависимость ошибки определения тропосферной задержки от входящих в формулу параметров.

3.6. Алгоритм метода точного позиционирования.

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

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

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

3.7.3. Об использовании измеренных значений метеопараметров при определении зенитных тропосферных задержек.

3.7.4. О точности определения координат пунктов в абсолютном методе в зависимости от способа учета влияния тропосферы.

3.7.5. Учет влияния тропосферы при относительных наблюдениях с использованием моделей и значений ЗТЗ, определенных методом точного позиционирования.

3.8. Выводы по главе 3.

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

Современные требования к точности координатных определений очень высокие для решения новых задач: мониторинга смещений и деформаций зданий и сооружений, особенно имеющих большую высоту; изучения перемещений литосферных плит (блоков); мониторинга деформаций и смещений земной поверхности в районах активных разломов - главных зонах генерации сильных землетрясений; прогнозирования землетрясений наблюдениям спутников; модернизации аэропортов; прокладки нефте- и газопроводов; кадастра и учета недвижимости; установки и эксплуатации буровых платформ на море; обеспечения безопасности различных видов транспорта (железнодорожное движение (особенно скоростное), воздушные полеты,.). Например, для координатного обеспечения кадастра точность определения координат в системе координат ITRF должна характеризоваться средними квадратическими ошибками менее 5 см по каждой из координат, а взаимное положение пунктов должно быть определено со средними квадратическими ошибками 1 см по каждой координате [6]. Взаимное положение пунктов, участвующих в геодинамических исследованиях, мониторинге смещений сооружений должно быть определено со средними квадратическими ошибками на уровне первых единиц миллиметров.

Однако с течением времени опорная геодезическая сеть постепенно «стареет», теряя свою первоначальную точность из-за влияния геодинамических эффектов. Координатным обеспечением на территории Социалистической Республики Вьетнам (СРВ) в настоящее время служит государственная геодезическая сеть (ГГС) 0, 1,2, 3-го классов и морская сеть, в том числе сети 0, 3-го классов и морская сеть, построенные спутниковым методом, и сети 1 и 2-го классов — триангуляционным методом. Прошло почти 15 лет с момента последнего построения и обновления ГТС О, 1 и 2-го классов и 10 лет с момента построения ГГС 3-го класса. Результаты последних геодинамических исследований показывают, что изменение координат пунктов в системе ГГИ7, находящихся на территории Вьетнама, составляет около 50 мм/год, а изменение взаимного расположения пунктов, находящихся на севере и юге СРВ, составляет 4-^-25 мм/год. Причем при построении ГГС СРВ не были учтены движения и возможные деформации литосферных блоков.

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

Для решения геодинамических задач в настоящее время широко применяются спутниковые технологии. Совершенствование и развитие спутниковых технологий с использованием глобальных навигационных спутниковых систем (ГНСС) позволят обеспечить построение однородной по точности геодезической сети на совершенно новых принципах при существенном сокращении временных затрат на построение сети. При использовании спутниковых технологий привязка исследуемой сети к Международной геодинамической сети также не представляет принципиальных трудностей.

Традиционно после выполнения наблюдений на пунктах не менее, чем за два периода, для определения геодинамических эффектов вычисляются, уравниваются координаты пунктов по каждому периоду. В результате будут получены средние положения пунктов в данный период. После этого, определяются изменения координат за интервал времени между периодами, и получается скорость движения пунктов. При реализации традиционного метода существуют следующие недостатки. Во-первых, во времени наблюдений каждого периода пункты также движутся. Если по каждому периоду вычисляются, уравниваются координаты пунктов, то не учитывается изменение координат пунктов во времени наблюдений каждого периода. Во-вторых, при геодинамических исследованиях более выигрышным экономически оказался подход, когда в общее уравнивание стараются включить максимальное число разных измерений, хотя бы они относились к разным эпохам. Эти измерения могут быть выполнены по разным программам, следовательно, в некоторые эпохи имеются данные наблюдений только на определяемых пунктах, а данные наблюдений на опорных пунктах отсутствуют. При реализации традиционного метода такие наблюдения не могут быть включены в уравнивание. Следовательно, не всегда имеется возможность включить в уравнивание максимальное число измерений. В-третьих, при реализации традиционного метода необходимо иметь данные о положении опорных пунктов на каждую t эпоху. Во избежание таких недостатков можно провести совместное уравнивание координат и составляющих скоростей движений пунктов. Этот метод позволяет учитывать движения пунктов во времени наблюдения каждого периода и также позволяет включить в обработку данные измерений в эпохи, когда данные наблюдений на опорных пунктах отсутствуют. Благодаря этому, точность полученных результатов будет выше. Однако недостаток метода совместного уравнивания координат и составляющих скоростей движений пунктов состоит в том, что размерность матриц, требующихся при реализации метода, большая. В настоящей диссертационной работе все расчеты проведены на персональном компьютере с характеристиками: процессор Intel Core2 Duo Т6600 2,2 ГГц; объем памяти 4096 Мб. Из-за ограничения объема памяти ЭВМ реализация этого метода при большом числе наблюдений невозможна. Способу решения этой проблемы будет посвящена настоящая работа.

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

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

Известно, что к настоящему времени разработано значительное количество- разных моделей, позволяющих оценить величину тропосферных задержек (ТЗ), влияющую на результаты спутниковых наблюдений: модель Хопфилд, модель Саастамойнена, модель Нейлла, модель MOPS (Minimum Operational Performance Standards), модель GCAT (GPS Code Analysis Tool), модель Университета Нью Брунсвика (УНБ) и модель Блэка. При учете влияния тропосферы по моделям всегда стоит вопрос о том, какая из них дает наиболее точные ТЗ в соответствии с местными метеорологическими условиями. Причем следует иметь в виду, что из-за муссонов территория Вьетнама отличается от других стран, находящихся на такой же широте: немного прохладнее зимой и не так жарко летом. Это приводит к необходимости анализа эффективности моделей тропосферы при спутниковых измерениях во Вьетнаме. Отсюда можно узнать, какая из моделей дает наиболее точные ТЗ в условиях Вьетнама, и нужно ли в перспективе строить модель, соответствующую метеорологическим условиям данной территории. Однако работ, посвященных исследованию этой проблемы в геодезической литературе, нами не обнаружено.

В последние годы для обработки результатов наблюдений спутников ГНСС, развит пакет программного обеспечения DDBase с открытым кодом, опубликованной в интернете [44]. Он позволяет вычислить компоненты базовых линий с учетом влияния тропосферы по моделям. Причем, наиболее точные по сравнению с модельным учетом значения задержек тропосферы находят уточнением значений зенитных тропосферных задержек, заранее определяемых по моделям, в процессе обработки результатов измерений с использованием метода точного позиционирования (Precise Point Positioning — РРР). В связи с этим возникает задача о разработке подпрограммы, включение которой в. программу DDBase позволит учитывать влияние тропосферы путем использования ЗТЗ, заранее полученных методом точного позиционирования, и исследовании эффективности этого способа над применением модельных учетов влияния-тропосферы.

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

Заключение Диссертация по теме "Геодезия", Нгуен Ван Донг

Основные результаты исследований по теме диссертации^ опубликованы в журнале Известия вузов «Геодезия и аэрофотосъемка»:

1. Краснорылов И.И., Нгуен Ван Донг, Сорокин Н. А. Исследование возмущенного движения космического аппарата в случае его запуска на орбиту Земли в треугольную точку либрации L4 системы Солнце-Земля // Изв. вузов. «Геодезия и аэрофотосъемка». 2009. Jte 5. С. 69 — 73.

2. Нгуен Ван Донг. Повышение точности Государственной геодезической сети Вьетнама // Изв. вузов. «Геодезия и аэрофотосъемка». - 2010. - № 3. - стр. 19-21.

3. Нгуен Ван Донг. Определение координат пунктов и их изменений из наблюдений, спутников, образующих глобальные навигационные спутниковые системы (на примере одного из районов Юго-Восточной Азии) // Изв. вузов. «Геодезия и аэрофотосъемка». - 2011. - № 4. - стр. 30 - 35.

4. Нгуен Ван Донг. Анализ моделей тропосферы по результатам спутниковых наблюдений в регионе Юго-Восточной Азии // Изв. вузов. «Геодезия и аэрофотосъемка». 2011. № 5. С. 21 - 25.

163

ЗАКЛЮЧЕНИЕ

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

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

Реализован в среде Visual С++ алгоритм метода совместного уравнивания координат и составляющих скоростей движений пунктов. Разработанная программа апробирована на сети действующих станций, находящихся вблизи территории Вьетнама. Полученная точность определения координат и составляющих скоростей движений станций этой сети предлагаемым методом составляет 7 + 17 мм и 2,0 5,4 мм/год соответственно. Точность горизонтальных составляющих скоростей составляет 1,1 1,7 мм/год, а вертикальной составляющей 4,4 6,6 мм/год. Точность определения положения станций и скорости предлагаемым методом в 2 3 раза выше, чем при использовании обычного метода. Из этого следует, что предлагаемый метод является эффективным для определения движений станций по результатам спутниковых наблюдений и может применяться для геодинамических исследований.

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

Реализован алгоритм метода оптимизации геодинамических сетей в среде Visual С++. С использованием разработанного программного обеспечения, проведена оптимизация двух геодинамических сетей (Дьенбьен-Лайчау и Сонла) во Вьетнаме. Практическое использование составленного программного обеспечения позволит в ряде случаев достичь существенной экономии финансовых и людских ресурсов за счет уменьшения количества подлежащих измерению величин и пунктов сети. И также предложены варианты повышения точности определения деформационных параметров включением в состав сети дополнительных пунктов. Разработанное программное обеспечение может быть внедрено при проектировании сетей, создаваемых для исследования движений и деформаций литосферных блоков; мониторинга смещений грунтов и сооружений; прогнозирования землетрясений по наблюдениям спутников.

3. Выполнено исследование определения ЗТЗ методом точного позиционирования, который позволяет определить координаты пунктов по результатам абсолютных наблюдений и уточнить поправку к измеренным дальностям за влияние тропосферы. Реализован алгоритм этого метода в среде Visual С++ с помощью библиотеки с открытым кодом [44].

Проведено уточнение ЗТЗ при спутниковых наблюдениях, которые заранее определены с использованием моделей Блэка, Хопфилд, GCAT, MOPS, Нейлла, УНБ и Саастамойнена. В результате получены ЗТЗ в сходных с Вьетнамом условиях с точностью на уровне 3-14 мм, а в среднем 6 мм.

4. Проведена оценка эффективности моделей тропосферы Блэка, Хопфилд, GCAT, MOPS, Нейлла, УНБ и Саастамойнена. Полученные результаты свидетельствуют о том, что для учета влияния тропосферы в сходных с Вьетнамом условиях при использовании модельных метеопараметров модель MOPS оказывается наиболее эффективной. Она позволяет учитывать не менее 92% (в среднем 98%) влияния тропосферы, причем с использованием моделей Хопфилд, Саастамойнена, Нейлла, GCAT, УНБ и Блэка учтено только 85% -ь 96% (в среднем 91%) влияния тропосферы. Если производить метеорологические измерения в момент наблюдений, использование моделей Блэка, Хопфилд и Саастамойнена. будет эффективным. Их использование дает возможность определять не менее 95% (в среднем 99%) влияния тропосферы.

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

• при определении координат абсолютным методом необходимо вводить поправку за влияние тропосферы,. как дополнительное неизвестное, и уточнять его в процессе обработки. В настоящей работе расхождения- положений пунктов, полученных нами, с решением МГС абсолютным методом составляют 3 см в случаях использования и модельных, и измеренных метеонараметров для учета влияния^ тропосферы. Также следует, что при определении координат методом точного позиционирования с уточнением ЗТЗ в процессе обработки не требуются измеренные метеопараметры.

• как показали наши вычисления, в случае, когда нет возможности уточнения: ЗТЗ и отсутствуют измеренные метеопараметры, модель MOPS оказывается наиболее эффективной для учета влияния: тропосферы при спутниковых наблюдениях в: сходных, с СРВ условиях. Согласно полученным в диссертационной работе данным, с использованием модели MOPS для учета влияния тропосферы отклонения положений, полученных нами' и МГС, составляют 0,04 0,25 м, а с использованием моделей Хопфилд, Саастамойнена, Нейлла, GCAT, УНБ и Блэка соответственно - 0,27 1,03 м.

• в случае, когда нет возможности уточнения ЗТЗ и присутствуют измеренные метеопараметры, модели Блэка, Хопфилд, Саастамойнена будут наиболее эффективными. При использовании этих моделей мы получили отклонения положений пунктов, полученных нами, с решением МТС 0,01 — 0.15 м. При использовании модели УНБ, отклонения составляют 0,45 + 0,58 м.

Из полученных результатов также следует, что учет влияния тропосферы по моделям при использовании измеренных метеопараметров дает возможность определять координаты с точностью в 2 - 10 раз выше, чем использование модельных метеопараметров.

Согласно результатам исследования предлагается, что при абсолютных измерениях мы< уточним влияние тропосферы в обработке результатов измерений. Если нет возможности уточнения ЗТЗ- для учета влияния тропосферы, предлагаем использование модели MOPS в случае отсутствии измеренных метеопараметров, а в противном случае, когда производились метеорологические измерения, предлагаем- использование модели Хопфилд, Саастамойнена или Блэка.

6. Выполнено исследование точности координатных определений относительным методом, в зависимости от способа учета влияния тропосферы. Разработана подпрограмма, включение которой в программу обработки результатов спутниковых измерений DDBase [44] позволит учитывать влияние тропосферы на относительные наблюдения путем использования ЗТЗ, полученных методом точного позиционирования. Из полученных результатов следует, что учет влияния тропосферы на относительные наблюдения путем использования ЗТЗ, полученных методом точного позиционирования позволяет получить точность координатных определений выше, чем при использовании моделей- тропосферы. Согласно полученным в работе данным, при расстоянии между пунктами 10 40 км и использовании ЗТЗ, полученных методом точного позиционирования, ошибки по каждой координате составляют 1,5 мм, а ошибки в положении — 2 мм. С использованием моделей для учета влияния тропосферы ошибки по каждой координате составляют 3-^-11 мм в зависимости от применяемой модели, а ошибки в положении - 4 + 15 мм. Если в относительном методе не учтено влияние тропосферы, средние квадратические ошибки по каждой координате достигают 1,5 см, а ошибки в положении 2 см.

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

Согласно имеющимся у нас данным есть основания полагать, что' полученные в диссертации результаты! будут учтены при модернизации координатного обеспечения СРВ. В частности, будет осуществлена привязка ГГС СРВ к геодинамической сети МТС (IGS).

Планируется также продолжение* нами- проведенных в диссертации исследований. Оно будет касаться обеспечения геодинамических исследований во Вьетнаме точной эфемеридной информации спутников^ GPS/TJIOHAGC. Основы для проведения такой работы заложены в нашей статье [1].

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

1. Абалакин В. К., Краснорылов И. И., Плахов Ю. В. Геодезическая астрономия и астрометрия. Москва : Картгеоцентр - Геодезиздат, 1996. - стр. 435.

2. Антонович К. М. Использование спутниковых радионавигационных систем в геодезии. Москва : Картоцентр, Новосибирск: Наука, 2005. - Т. 1 : стр. 334.

3. Антонович K.M. Использование спутниковых радионавигационных систем в геодезии: Москва : Картоцентр, Новосибирск: Наука; 2006. - Т. 2 : стр. 360.

4. Аплонов C.B. Геодинамика. СПб : СПбГУ, 2001. - стр. 360.

5. Баранов В.Н., Бойко Е.Г., Краснорылов И.И., Машимов М.М., Плахов Ю.ВЧ, Урмаев М.С., Яшкин С.Н. Космическая геодезия. Москва : Недра, 1986. - стр. 407.

6. Бойков A.B. Теоретические основы и практическая реализация координатного, обеспечения спутниковой системы межевания земель: проект «Москва». Дис. канд. техн. наук / МИИГАИК. Москва, 2008. - стр. 184.

7. Большаков В.Д., Маркузе Ю.И., Голубев В.В. Уравнивание геодезических построений: Справочное пособие. Москва : Недра, 1989. - стр. 413.

8. Валландер С. В. Лекции по гидроаэромеханике. Л.: Изд. ЛГУ, 1978. -стр. 296.

9. Войтенко А. В. Применение калибровочных районов при топографо-геодезических работах со спутниковыми приемниками GPS // Геопрофи. Москва. 2006. № 1. С. 25 27.

10. Генике А. А., Побединский Г. Г. Глобальные спутникове системы определения местоположения и их применение в геодезии. Москва : Картгеоцентр, 2004. - стр. 355.

11. П.Герасименко М. Д. Оптимальное проектирование и уравнивание геодезических сетей. Москва : Наука, 1992. - стр. 160.

12. Герасименко М. Д., Касахара М. К вопросу об учете физической корреляции при оптимальном проектировании геодезических измерений // Изв. вузов «Геодезия и аэрофотосъемка». 2005. № 1. С. 17 21.

13. Герасименко М. Д., Коломец А.Г. Исследование математических проблем при* построении и обработке космических геодезических сетей для целей геодинамики' // Дальневосточный математический журнал. 2008. Том 8. № 1.С. 31-45.

14. Демьянов Г. В. Государственные геодезические сети: современное состояние и перспективы развития // Информационный бюллетень, ГИС-Ассоциации. 2009. № 5 (72). С. 41 42.

15. Крылов В. И. Космическая геодезия. Москва : МИИГАиК, 2002. -стр. 175.

16. Кузин' С.П., Махматгазиев Б. Международная геодинамическая служба-IGS // Сборник научных трудов^«Космическая4 геодезия и современная геодинамика» / ред. А.Г. Масевича. Москва : МИД, 1996. - стр. 253.

17. Маркузе Ю.И., Голубев- В.В. Теория математической обработки геодезических измерений: Учебное пособие для вузов. Москва: Академический Проект; Альма Матер, 2010. - стр. 247.

18. Машимов М.М. Уравнивание геодезических сетей. Москва : Недра, 1979. - стр. 367.

19. Михайлов Ю.В., Морозов В.Н., Татаринов В.Н. Глобальные навигационные спутниковые системы в обеспечении геодинамической безопасности разработки рудных месторождений. Москва : МГОУ, 2008.

20. Монахова М.А. Методы применения спутниковой системы (проект «Москва») для геодезического обеспечения кадастра объектов» недвижимости. Дис. канд. техн. наук / МИИГАИК. Москва, 2007. - стр. 163.

21. Нгуен Ван Донг Анализ моделей тропосферы по результатам спутниковых наблюдений в регионе Юго-Восточной Азии // Изв. вузов «Геодезия и аэрофотосъемка». 2011. № 5. С. 21 25.

22. Нгуен Ван Донг Повышение точности Государственной геодезической' сети Вьетнама // Изв. вузов «Геодезия^ и аэрофотосъемка». 2010. № 3. С. 19-21.26: Рульков Д.И. Навигация 1 и, лоция. Москва-: Транспорт, 1973.- стр. 232: ' ^

23. Самратов У.Д., Филатов В.Н: О современном состоянии и путях модернизации Государственной геодезической сети Российской Федерации// Информационный бюллетень. ГИС-Ассоциации. 2009. №>5 (72). С. 35-40.

24. Хаин В.Е. Тектоника континентов и океанов. Москва : Научный мир, 2001.-стр. 606.

25. Хаин В.Е., Ломизе М.Г. Геотектоника с основами геодинамики. -Москва : МГУ, 1995. стр. 480.

26. Шестаков Н.В. Оптимальное проектирование деформационных ОРБ/ГЛОНАСС сетей. Дис. канд. техн. наук/ Институт автоматики и процессов управления ДВО РАН. Владивосток, 2005. - стр. 159.1 ?

27. Шестаков Н.В., Герасименко М.Д. К вопросу об оптимальном проектировании деформационных геодезических GPS-сетей // Изв. вузов «Геодезия и аэрофотосъемка». 2009. № 5. С. 11 17.

28. Яковлев Н.В. Высшая геодезия: Учебник для вузов. Москва : Недра, 1989. - стр. 445.

29. А С++ Library of Genetic Algorithm Components. http://lancet.mit.edu/ga/.

30. Alfred Leick. GPS Satellite Surveying. 2004. -3rd : p. 464.

31. Black H. D. An easily implemented algorithm for the tropospheric range correction // Journal of Geophysical Research. 1978. Vol 83. P. 1825-1828.

32. Black H. D., Eisner A. Correcting Satellite Doppler Data for Tropospheric Effects //Journal of Geophysical Research. 1978. Vol, 89. P. 2616-2626.

33. Bo Tai nguyen vä Moi trirong Viet Nam1. http://www.monre.gov.vn.

34. Dennis D. McCarthy, Gerard Petit. IERS Technical1 Note No. 32. IERS Conventions / International-Earth Rotation and Reference Systems Service. 2003. -p. 127.

35. Duong Chi Cong, Hong-Sic Yun, Jae-Myoung Cho. GPS measurements of horizontal deformation across the Lai Chau -Dien Bien (Dien Bien Phu) fault, in Northwest of Vietnam 2002 2004 // J. Earth, Planets Space. 2006. Vol. 58. P. 523-528.

36. Frank Kleijer. Troposphere Modeling and Filtering for Precise GPS Leveling: Delft : Netherlands Geodetic Commission, 2004. - p. 282.

37. GNSS data // The Crustal Dynamics Data: Information System. -ftp://cddis.gsfc.nasa.gov.

38. GPS ToolKit Software Library.- http://www.gpstk.org.

39. Grafarend E.W. Optimization of Geodetic Networks // Geod. Sei. Äff. 1974. No. 4 : Vol 33. P. 351-406.

40. Grids & Datums Vietnam // Photogrammetric' engineering and-Remote Sensing. Journal of the American society for photogrammetry and- remote sensing. 2002. Vol. 68. , •

41. Kuang Shan-Long Optimization and Design of Deformation Monitoring Schemes. : Technical report. No 157/ Department of Geodesy and Geomatics Engineering ; University of New Brunswick. -1991.

42. Lai K. Y., Chen L. Chung, D. D. Lam Active tectonics and Quaternary basin formation along: the Dien Bien Phu fault zone, northwestern Vietnam // Ji Geology. Series B. Hanoi. 2006. No. 27. P. 49-57.

43. Mevlut Yetkin, Cevat Inal, Cemal Ozer Yigit Optimal Design of Deformation Monitoring Networks by Using PSO Algorithm // 13th FIG Symposium on Deformation Measurements and Analysis, 4th IAG Symposium on Geodesy for Geotechnical. Lisbon, Portugal, 2008.

44. Minimum; Operational Performance Standards for Global Positioning System/Wide Area Augmentation System Airborne Equipment. RTCA/DO-229C. -2004.

45. Nikolaidis R. Observation: of geodetic and seismic deformation with the Global Positioning System: Ph.D. Thesis / University of California. San Diego, 2002.

46. Pesci A., Loddo F., Cenni N., Teza G.,, Gasulai Gr Analyzing Virtual Reference Station for GPS surveying: experiments and applications in a test site of the northern Apennine (Italy) // Annals of Geophysics. 2008: No. 4 :. Volt 51. P. 633 6531

47. The Toolkit for GPS/GLONASS/Galileo/SBAS Data // University Navstar Gonsortium. http://facility.unavco.org/software/teqc/teqc.html.

48. Tran Dinh To» Processing results of GPS measurement? data on from the Son La and Song Da fault zones // Journal of Geology. Series B. Hanoi; 2006; No. 27. P: 115-122.

49. Transformation parameters // The International' Terrestrial*; Reference Frame. -http://itrf.ensg.ign:fr/transpara.php.

50. Trimble DocuShare // Trimble http://trl;trimble.com/dscgi/ds:py/Get/File-544869/SPEPOCH50UGlA.pdf.

51. Trimble DocuShare // Trimble. http://trl.trimble.com/docushare/dsweb/ Get/ Document-148711/022543-079J-RUTrimbleR8GNS SDS1109LR.pdf.

52. Zbigniew Michalewicz. Genetic Algorithms + Data Structures = Evolution Programs. New York : Springer, 1996. - p. 387.r r

53. Bô tài nguyên và môi trucmg Viêt Nam. Six dung hê thông tham sô tinhr r rchuyên giûa hê toa dô quôc tê WGS-84 và hê toa dô quôc gia Viêt Nam VN-2000. Quyét dinh sô 05/2007/QD-BTNMT ngày 27 thâng 2 nàm 2007. HàNôi. - 2007.

54. Cao Dinh Triêu, Nguyên Danh' Soan. Hê thông durt gây chinh lânh thô Viêt Nam trên ca sà phân tich kêt hop tài lieu trong lue, tù và ânh vê tinh // Tap chi Dia chât. Loat A. HàNôi. 1998. Sô 247. Tr. 17-27.> r r >

55. Cao Dinh Triêu, Thâi Anh Tuân, Cao Dinh Trong. Mot sô nét dac trung vêr r r rkiên tao dia chân khu vue Dông Nam A // Tap chi'Dia chât. Loat A. Hà Nôi. 2008. So 306. Tr. 3-13.w r r

56. Dô Ngoc Duàng. Giâo trinh trâc dia cao câp. Hà Nôi : Truorng Dai hoc Mo - Daij chât. - 2008i - tr. 1702r t r r

57. Ky yêu hôi nghi tông kêt, giôi thiêu luôi toa dô quôc gia hang III, bân dôr F •»dia hinh ty le 1:50000 phû trùm toàn quôc và câc tram GPS quôc gia. HàNôi : Tông eue Dia Chinh. - 2004.

58. Nguyên-Van Sang. Thuât toàn binh sai kêt hop lirai GPS nhiêu chu ky trong quan trâc dia dông. Dai hoc Mo Dia chât.1.r r *

59. Tông eue dia chinh. Bâo câo xây dung hê quy chiêu«và hê thông diêm toa. dô quoe gia. HàNôi. 1999, 250 tr.

60. Tông eue thông kê Viêt Nam. http://www.gso.gov.vn.

61. Vu Van Tich, Chu Van Ngoi, Luorng Thi Thu Hoài, Duong Thi Toan,r t r

62. Pham Khâc Hùng. Dac diêm biên dang dai«dut gay hoat dông Diên Biên Lai Châu và tiêm nâng dia nhiêt vùng U Va,,Tây Nam trung Diên Biên // Tap chi Dia chât. Loat A. HàNôi. 2009. So 313.