Гравиметрическое обеспечение прилегающей к Вьетнаму акватории по результатам спутниковой альтиметрии

Нгуен Ван Шанг

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

Автореферат диссертации по теме "Гравиметрическое обеспечение прилегающей к Вьетнаму акватории по результатам спутниковой альтиметрии"

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

НГУЕН ВАН ШАНГ

ГРАВИМЕТРИЧЕСКОЕ ОБЕСПЕЧЕНИЕ ПРИЛЕГАЮЩЕЙ К ВЬЕТНАМУ АКВАТОРИИ ПО РЕЗУЛЬТАТАМ СПУТНИКОВОЙ АЛЬТИМЕТРИИ

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

Автореферат

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

- 8 НОЯ 2012

Москва-2012

005054496

Диссертация выполнена на кафедре Высшей математики Московского государственного университета геодезии и картографии

(МИИГАиК)

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

Нейман Юрий Михайлович

Официальные оппоненты:

Лебедев Святослав Владимирович, доктор технических наук, доцент, Московский государственный университет геодезии и картографии (МИИГАиК), профессор кафедры Высшей геодезии.

Плешаков Дмитрий Иванович, кандидат технических наук, ФБУ «27 ЦНИИ МО РФ», начальник Научно-исследовательского управления Научно-исследовательского центра (топогеодезического и навигационного обеспечения). Ведущая организация:

Центральный ордена «Знак почёта» научно-исследовательский институт геодезии, аэросъёмки и картографии имени Ф. Н. Красовского (ФГУП «ЦНИИГАиК»),

¡с- и

Защита состоится « ... » ноября 2012 г. в часов на заседании диссертационного

совета Д.212.143.03 в Московском государственном университете геодезии и картографии (МИИГАиК) по адресу: 105064 Москва, Гороховский пер., 4, зал заседания Ученого совета.

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

Автореферат разослан « 11» октября 2012 г.

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

диссертационного совета —-Климков Юрий Михайлович

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

Обоснование актуальности темы диссертации

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

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

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

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

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

2) адаптация программного пакета BRAT (Basic Radar Altimetry Toolbox) для решения конкретных задач и вычисление с его помощью высот морской поверхности на акватории вблизи Вьетнама;

3) выделение динамических возмущений морской поверхности путём уравнивания пересечений спутниковых треков;

4) исследование степени гладкости гравитационного поля на изучаемой территории;

5) практическое вычисление аномалий силы тяжести на акватории, прилегающей к Вьетнаму, по данным альтиметрии ENVISAT (Environmental Satellite);

6) исследование точности полученных результатов и исследование возможностей их уточнения по доступным результатам гравиметрической съёмки. Научная новизна работы

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

- учёт реальной специфики акватории,

- учёт реальной гладкости соответствующего гравитационного поля,

- ослабление влияния динамической топографии,

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

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

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

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

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

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

4 - Регулярная сетка с шагом 10' значений аномалии силы тяжести, полученная на прилегающей к Вьетнаму акватории по результатам альтиметрии ENVISAT. Вклад автора в проведенное исследование

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

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

2 - На прилегающей к Вьетнаму акватории по результатам альтиметрии ENVISAT получены конкретные значения аномалии силы тяжести на регулярной сетке с шагом 10'.

3 - Даны рекомендации по работе с современными форматами космической геодезии при использовании программного пакета BRAT.

Апробация результатов работы

Основные результаты по теме диссертации докладывались на конференциях студентов, аспирантов и молодых ученых МИИГАиК 6 апреля 2011 и 3 апреля 2012 гг.

Публикации

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

Диссертация состоит из введения, трёх глав основного текста и заключения. Список литературы насчитывает 46 наименований на русском, вьетнамском и английском языках. Кроме того, имеется список использованных сокращений, список принятых обозначений и 6 приложений. Общий объем работы составляет 157 страниц машинописного текста и содержит 71 рисунок и 66 таблиц.

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

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

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

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

1. Принцип «удаления - восстановления» при обработке данных альтиметрии

Высота морской поверхности SSH (Sea Surface Height) представляется формулой

SSH = Negm + AN+ hd. (!)

Согласно принципу «удаления - восстановления», до вычисления аномалии силы тяжести, высоты геоида ^относительно EGM удаляются. В результате роль исходных данных играют только остаточные значения AN и динамическая высота моря hd. Такое вычитание длинноволновой части поля в определённом смысле равносильно удалению влияния "дальних зон" и уменьшению радиуса корреляции гравитационного поля. Полезным следствием является возможность выполнять последующие вычисления в меньших регионах, учитывая остаточное влияние только "ближних зон". Понятно, что в конце вычислений низкочастотную часть поля необходимо восстановить. Высоты геоида относительно EGM вычислены по формуле:

Z ~ Ё(с«." C0SOT¿+ *»,», sin тЯ)^, (sin .

/V -ЯМ.

цел I

у.г

(2)

Модельные аномалии силы тяжести относительно ЕвМ вычисляются по формуле:

Д?£СИ = —

Z " ("-')]£ (с».ш c°s тЯ + s„ ,„ sin тЛ)р„ „, (sin ф)

(3)

2. Вычисление высоты средней динамической топографии моря

Динамическая высота моря hd складывается из высоты средней динамической топографии моря (Mean Dynamic Topography)/,„„,. и динамической высоты моря, изменяющейся со временем И,. Высота средней динамической топографии моря является разностью между поверхностью геоида и средней морской поверхностью Её значения изменяются в пределах ±1,8 м и во многом обусловлены морскими течениями. Высоту средней динамической топографии моря можно вычислить с помощью моделей геоида (EGM) и моделей средней морской поверхности (Mean Sea Surface - MSS) по формуле

hMDT = hilss -Ntm/ (4)

или с помощью моделей средней динамической топографии моря (Mean Dynamic Topography - MDT).

В данной работе использованы модели MSS_GNES CLS 10 и

GO_CONS_EGM_GOC_2I. Высота средней динамической топографии моря на акватории, прилегающей к Вьетнаму, вычислена с помощью программы GUT. В результате мы получили 15 053 значений высоты средней динамической топографии моря на регулярной сетке с размером ячеек 6 х 6'.

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

тиметрии, то неизбежно возникает необходимость интерполяции. В нашей работе эта задача решается методом оптимального линейного прогнозирования с помощью нашей программы «/п(ег_етр», «1п!ег_/И» и «/ШегроШюп».

3. Уравнивание пересечений треков

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

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

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

- классификация нисходящих треков и восходящих треков.

- сравнение каждого восходящего трека со всеми нисходящими треками с целью поиска пересечения.

- определение приближенного места пересечения.

- определение уточнённого места пересечения треков.

102 104 106 108 110 и: 114 Рис. 1. Пересечения треков спутника ENVISAT

- вычисление разности высоты на пересечении.

Уравнивание результатов пересечения треков

Для интересующей нас акватории Вьетнама длины треков достигают 1500 км, динамические высоты моря каждого трека могут быть смоделированы в предположении постоянства смещения а и уклона Ъ трека, то есть

еШ„ = (а, + Ь,^) - (а, + Ь^,) + , (5)

где Ь,, Ь, - неизвестные параметры уклона, которые относятся к восходящему и нисходящему треку; //,, ц] - относительные долготы ;'-го иу'-го треков;

= (6)

где относительная долгота точки Р на /'-ом треке; ?!,- долгота точки Р на /-ом треке; Х>- средняя долгота /'-го трека.

В матричном виде эти уравнения связи принимают вид:

й = Ах + \>, (7)

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

*Г=к Ьп а,2 Ьп ... ат Ь„„ в/| Ь,х ап Ър ... аА (8)

А - матрица коэффициентов уравнений связи;

й - матрица разностей высот морской поверхности в точках пересечения тре-

ков

är=(dH, dH2 ... dHn); (9)

v - матрица остатков

vr=(v, v2 ... v„). (Ю)

Матрица коэффициентов уравнений связи (7) имеет ранг, на четыре меньший полного. Уравнения связи можно решать или путём предварительной фиксации каких-нибудь двух треков, или совместным уравниванием пересечения всех треков под условием хТС;[х = min.

Недостатком обоих этих подходов является то, что после уравнивания результатов пересечений треков разность между поверхностью, созданной из уравненных точек пересечения треков, и поверхностью геоида иногда остаётся довольно большой. Для уменьшения этой разности желательно приближать каждый трек к геоиду. В связи с этим высоты морской поверхности каждой точки на А-ом треке рекомендуется смоделировать в предположении постоянства смещения (а°) и уклона (6°) трека, то есть

(И)

Здесь SSHtii - высота морской поверхности после удаления высоты геоида относительно модели гравитационного поля Земли NCOil и высоты средней динамической

топографии моря hilDT; - относительная долгота точки /'ь Vh - остаток высоты в точке /].

Параметры а° и Ь°. оцениваются методом наименьших квадратов под условием jV,* = min, щ - количество измеренных точек на к-ом треке.

Приближение каждого трека к поверхности геоида и уравнивание пересечений треков можно выполнять одновременно (т.е. <я° = а и b° = Ь). Тогда система уравнений связи принимает вид:

¿Я, = (а, + b,n¡tI - (а, + b¡M¡) + v..

SSH^=(ak+bif¡tl¡) + Vtl¡ ■ { >

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

£v¿ + v£X=min. (13)

Величина w является относительным весом. Если w мало, то разность в пересечениях уменьшается достаточно эффективно, но средняя морская поверхность не приближена к поверхности геоида. Наоборот, если w велико, то средняя морская поверхность хорошо приближена к поверхности геоида, но разность в пересечениях всё ещё остаётся значительной. В расчётах относительный вес принят w = 0,1.

После того как параметры смещения уклона оценены, остаточные высоты геоида вычисляются по этим параметрам:

AiV = SSH"-(a + b/j). (14)

4. Определение аномалии силы тяжести по данным альтиметрии методом кол-локации

Трактуем величины остаточных высот геоида AÑt,AÑ2,...,AÑn как измеренные значения функционалов ¿„ i = 1,2, ..., п на соответствующем возмущающем потенциале. Требуется вычислить значение другого функционала F - остаточную аномалию силы тяжести Agp в точке заданной Р.

Согласно локальной задаче коллокации,

F = Kr(L,F).(K(L,L) + CJ-4. (15)

Здесь F - оценка функционала F, К(\ •) - воспроизводящее ядро используемого гильбертова пространства (ковариационная функция), Сд - ковариационная матрица ошибок измерений, \ результаты измерений функционалов L¡, i =1,2,..., п.

Таким образом, аномалия силы тяжести AgP в точке Р вычисляется по следующей формуле

&gP = Kl„(AN,Ag,,){K(AN,AN) + C^AV,,,,. (16)

Соответствующая дисперсия оценивается величиной . = км( ) - (ш, лёр дло+.

Здесь Kт(AN,Agp) = [K(&Nl,AgP) K(AN2,Agp)

~К(Щ,Щ) К(Щ,АКг) К(АМ2,Щ) К(АЫ2,Щ)

К(АМ„,Аёр)]; К(Щ,Щ)

С,

ДУУ =

(17)

(18)

(19)

(20)

Ковариационная функция остаточных высот геоида записывается в формуле

м У,У,

я2

Р,{ ССКу/) +

Р,( СОЧ/). (21)

ьи (/-2)(/ + 6) ^

Взаимная ковариационная функция аномалий силы тяжести и остаточной высоты геоида имеет вид:

и-Ы Я

У, ы

1-1

+ г,,

Р,(«в«0, (22)

А ковариационная функция аномалий силы тяжести имеет вид:

ДА?,,Л?,) =

(/-1)г

'¡-О

Р,(соя//)+ £

(1-1)

Р,( соя?), (23)

где Р^соя!//) - полиномы Лежандра степени у/ - сферическое расстояние между точками ¡' и у; г, и - рассстояния этих точек от начала координат; с!, - степенные дисперсии ошибок коэффициентов низкочастотной части потенциала до N включительно; Я - средний радиус Земли; Ъ - число, которое выбирается равным 4, но иногда, чтобы максимально приблизиться к низким степенным дисперсиям, достигает 24;. Степень М, масштабный коэффициент а, константа Л в единицах (м/с)4 и радиус целиком расположенной в теле Земли сферы Бьерхаммара Яя являются свободными параметрами и подлежат определению по результатам ковариационного анализа остаточных значений высот геоида. В нашем случае среднее расстояние Ау/ между точками с данными приближённо равно 10 дуговых минут.

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

д? = Дг + Л^,г;л( , (24)

где ~ аномалия силы тяжести относительно модели гравитационного поля Земли вычисляется по формуле (3).

5. Уточнение аномалий силы тяжести по результатам гравиметрической съём-

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

В самом деле, пусть на изучаемой акватории имеются к значений аномалий силы тяжести Д.£,'Г"", Л£Г'"> вычисленных по данным альтиметрии, и т зна-

чений аномалий силы тяжести Д^"", 7%"™,..., полученных по результатам гравиметрической съёмки на море или на ближайшем побережье. Тогда уточненное значение аномалия силы тяжести &gr в точке Р можно получить по формуле:

""Мр) г

ЦА8" Кт{

-1 V""

М".

(25)

Соответствующая дисперсия оценивается величиной

К( Ag^'•■',Agp)

^ = К^&вг)

-1 'К( Дg"• ".аг,)"

К(

) + Сл~д_ АГ^—.Дв"")

Здесь К(\ •) - ковариационная функция аномалии силы тяжести, Сд ная матрица ошибок измерений.

Ковариационная функция аномалии силы тяжести вычислена аналогично по формуле (23). Параметры

а, <Л/, /V, А, и Ид подлежат определению по результатам корреляционного анализа эмпирических данных аномалий силы тяжести.

(26)

Гдд - ковариацион-

6. Общая схема определения аномалий силы тяжести по данным альтиметрии ЕтавАт

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

В схеме: ББН - высота морской поверхности; Ицам - высота геоида относительно определённой модели гравитационного поля Земли (ГПЗ); Д/У - остаточная высоты геоида; Имот- высоты средней динамической топографии моря; И,- динамическая высота моря, изменяющаяся со временем; д^ - остаточная аномалия силы тяжести; ^есм - аномалия силы тяжести относительно модели (ГПЗ); Дg - полная аномалия силы тяжести; Аёусред - усреднённая аномалия силы тяжести по данным многих циклов; Ag¡ - аномалия силы тяжести ¡-го цикла; п - число циклов.

Рис. 2. Схема определения аномалий силы тяжести по данным спутниковой альтиметрии ЕКУГЭАТ.

В третьей главе «Практические результаты вычисления аномалий силы тяжести на акватории, прилегающей к Вьетнаму, по данным альтиметрии ЕЫУ^АТ» приведены экспериментальные результаты вычисления аномалии силы тяжести на акватории, прилегающей к Вьетнаму, по данным альтиметрии ЕЫУГБАТ, результаты сравнения с каталогом О^СШСЯАУ, результаты сравнения с результатами гравиметрической съёмки на побережье и результаты совместного уравнивания аномалии силы тяжести, полученные по результатам гравиметрической съёмки на побережье и по данным альтиметрии.

1. Результаты вычисления аномалий силы тяжести при использовании ЕСМ96 для процедуры «удалення-восстановлення»

В этом разделе описаны результаты вычисления аномалии силы тяжести на акватории, прилегающей к Вьетнаму (ср: 8° - 22°, к 105° - 114°), в соответствии с изложенной выше теорией. Исходными данными для вычислений послужили высоты морской поверхности, полученные со спутника ЕМУ15АТ на десяти циклах и приведенные в первой главе. Для процедуры «удаления - восстановления» использована известная модель глобального гравитационного поля Земли ЕСМ96. Аномалии силы тяжести вычислены для каждого цикла и затем усреднены по десяти циклам.

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

Здесь Аё>'с<*'> - усредненная аномалия силы тяжести для десяти циклов в /-й точке; &g¡J- аномалия силы тяжести у'-го цикла в г'-й точке.

В результате для изучаемого района получена регулярная сетка 10' х 10' усредненных аномалий силы тяжести, общее количество точек 4675, максимальное значение равно 70,8 мгал, минимальное значение равно -62,3 мгал, среднее значение равно -3,8 мгал, среднеквадратическое значение равно 16,5 мгал.

Таблица 1 содержит краткие статистики сравнений аномалий силы тяжести для каждого цикла с усреднёнными результатами для всех десяти циклов.

Здесь п - общее количество точек;

I 10

(27)

<%/ = -

,усред .

/ = 1,2,...,«; у = 81, 83, 84, 85, 86, 87, 88, 89, 90, 91; - максимальное отклонение аномалий силы тяжести;

- минимальное отклонение аномалий силы тяжести; ¿%ср - среднее отклонение аномалий силы тяжести

п ы

а^ - среднеквадратическое отклонение

Таблица 1. Краткие статистики сравнений каждого цикла с усреднёнными результатами

Цикл п ¿£ио> (мгал) <%„„>гал) ^(мгал) (мгал)

81 -ср 4675 11,6 -17,3 -0,3 2,5

83-ср 4675 14,5 -13,9 -0,0 2,9

84-ср 4675 11,1 -13,6 0,0 2,4

85-ср 4675 23,8 -16,1 -0,2 3,5

86-ср 4675 14,7 -24,2 -0,1 2,9

87-ср 4675 31,2 -4,9 0,3 3,0

88-ср 4675 19,9 -20,9 0,1 3,2

89-ср 4675 24,0 -13,1 -0,0 3,2

90-ср 4675 40,4 -56,9 0,4 7,9

91 - ср 4675 11,0 -10,2 -0,2 2,6

Из таблицы 1 видно, что результаты 90-го цикла сильно отличаются от усредненных значений, а другие циклы от усредненных значений отличаются мало. Поэтому мы удалили результаты 90-го цикла и ещё раз усреднили остаточные девять циклов.

Усредненные аномалии силы тяжести девяти циклов вычислены по формуле:

(31)

где - усредненная аномалия силы тяжести для девяти циклов в /'-й точке; &glJ -аномалия силы тяжестиу'-го цикла в /-й точке.

Краткая статистика усредненных результатов: максимальное значение = 70,7 мгал, минимальное значение = -62,5 мгал, среднее значение = -3,9 мгал, среднеквадратическое значение = 16,4 мгал. Усредненные аномалии силы тяжести для девяти циклов показаны на рисунке За.

(30)

(а) при использовании ЕОМ96 для проце- (б) при использовании ЕОМ2008 для про-дуры «удаления-восстановления». Высота цедуры «удаления-восстановления». Высо-сечения 5 мгал та сечения 10 мгал

Рис. 3. Усредненные аномалии силы тяжести для девяти циклов

2. Результаты вычисления аномалий силы тяжести при использовании ЕСМ2008для процедуры «удаления-восстановления»

В этом разделе описаны результаты вычисления аномалии силы тяжести на акватории, прилегающей к Вьетнаму. Исходными данными для вычислений послужили высоты морской поверхности, полученные со спутника ЕТчМБАТ на десяти циклах и приведенные в первой главе. Для процедуры «удаления - восстановления» использована известная модель глобального гравитационного поля Земли ЕСМ2008. Аномалии силы тяжести вычислены для каждого цикла и затем усреднены по десяти циклам.

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

1 10

Л^'= -!>,. (32)

Здесь - усредненная аномалия силы тяжести для десяти циклов в /-й точке;

Л?,у — аномалия силы тяжестиу'-го цикла в /-й точке.

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

значение равно 137,8 мгап, минимальное значение равно -86,0 мгал, среднее значение равно -4,2 мгал, среднеквадратическое значение равно 21,2 мгал.

Таблица 3 содержит краткие статистики сравнений аномалий силы тяжести для каждого цикла с усреднёнными результатами для всех десяти циклов.

Таблица 2. Краткие статистики сравнений каждого цикла с усредненными результатами

Цикл п <%,0>гал) ¿¡^ (мгал) <%ср (мгал) а^ (м.гал)

81 -Ср 4675 12,3 -16,6 -0,1 2,4

83-Ср 4675 15,1 -18,9 0,1 2,8

84-Ср 4675 10,0 -11,8 -0,0 2,1

85-Ср 4675 40,8 -22,4 -0,1 3,5

86-Ср 4675 16,8 -26,5 -0,2 2,8

87-Ср 4675 35,5 -19,5 0,1 3,1

88-Ср 4675 25,4 -25,6 0,0 3,3

89-Ср 4675 21,6 -13,3 0,1 3,0

90-Ср 4675 68,7 -80,4 0,1 9,3

91-Ср 4675 9,2 -10,6 -0,1 2,4

Из таблицы 3 видно, что результаты 90-го цикла сильно отличаются от усредненных значений, а другие циклы от усредненных значений отличаются мало. Поэтому мы удалили результаты 90-го цикла и ещё раз усреднили значения остаточных девяти циклов.

Усредненные аномалии силы тяжести девяти циклов вычислены по формуле

= (зз)

где Д^"1 - усредненная аномалия силы тяжести девяти циклов в г-й точке; аномалия силы тяжестиу'-го цикла в ;'-й точке.

Краткая статистика усредненных результатов: максимальное значение = 137,8 мгал, минимальное значение = -86,1 мгал, среднее значение = -4,2 мгал, среднеквадратическое значение = 21,2 мгал. Усредненные аномалии силы тяжести для девяти циклов показаны на рисунке 36.

3. Сравнение аномалий силы тяжести, полученных при использовании для процедуры «удаления - восстановления» моделей ЕвМ96 и ЕСМ2008

Вычисленные результаты аномалий силы тяжести, полученные при использовании для процедуры «удаления - восстановления» моделей ЕОМ96 и ЕСМ2008, сравнены между собой. 16

Результаты сравнений на 4675 точках представлены горизонталями на рис. 5 с высотой сечения 10 мгал. На рис. 4 показана гистограмма распределения частоты отклонений. Краткая статистика результатов сравнений: среднее отклонение равно 0,3 мгал, и среднеквадратическое отклонение равно 14,1 мгал.

■110 -100 -90 -80 -70 -60 -50 -40 -30 -20 -10 0 10 20 30 40 50 60 70 80

Разность (мгал)

Рис. 4. Гистограмма распределения частоты отклонений

Рис. 5. Разность аномалий силы тяжести, полученных при использовании ЕСМ96 и ЕСМ2008 для процедуры «удаления - восстановления» на акватории, прилегающей к Вьетнаму

Из рис. 5 видно, что на суше и островах, где нет данных альтиметрии, разности большие, а на море, где есть данные альтиметрии, разности маленькие. Если рассмотреть только морской район (широта: 08° н- 18°, долгота: 109°30' - 114°), окруженный синей линией, то среднее отклонение равно 1,0 мгал, и среднеквадратическое отклонение равно 10,6 мгал.

4. Сравнение вычисленных аномалий силы тяжести с каталогом Б^С08СКА

ОЫ5С08СЯА является одним из наиболее подробных современных глобальных каталогов аномалии силы тяжести, созданный на сетке с размером ячеек 1' х 1' датским национальным космическим центром в 2008 году. Среднее отличие данных этого каталога от 321 400 результатов морской гравиметрии составило 0,39 мгал, а среднеквадратическое отклонение равно 4,36 мгал.

Сравнение с данными каталога 0№С0801А аномалий силы тяжести, вычисленных при использовании для процедуры «удаления — восстановления» модели ЕСМ96

Вычисленные аномалии силы тяжести девяти циклов, полученные при использовании ЕвМ96 для процедуры «удаления - восстановления», сравнены с каталогом ВЫ8С08С11А. Результаты сравнения представлены горизонталями на рис. 7а. Гистограмма распределения частоты отклонений показана на рис. 6. Краткая статистика отклонений аномалий силы тяжести: среднее отклонение равно 0,4 мгал; среднеквадратическое отклонение равно 14,2 мгал.

Рис. 6. Гистограмма распределения частоты отклонений, сравненных с каталогом ОЫ5С08011А

Из рисунка 7а видно, что на суше и островах, где нет данных альтиметрии, отклонения аномалий силы тяжести большие, а на море, где есть данные альтиметрии, отклонения маленькие. Если рассмотреть только морской район (широта: 08° ^ 18°, долгота: 109°30' 114°, окруженный синей линией), то среднее отклонение равно всего 0,9 мгал, но среднеквадратическое отклонение равно 11,1 мгал.

Ь Г *

V Л

(а) при использовании ЕСМ96 для процедуры «удаления - восстановления». Высота сечения 10 мгал

(б) при использовании ЕСМ2008 для процедуры «удаления - восстановления». Высота сечения 2 мгал

Рис. 7. Отклонения аномалий силы тяжести, сравненные с каталогом ОЫБСЖОКА

Сравнение с данными каталога ОЮСОЗСЛА аномалий силы тяжести, вычисленных при использовании для процедуры «удаления - восстановления» модели ЕСМ2008

Аномалии силы тяжести, вычисленные с использованием в процедуре «удаления - восстановления» модели ЕСМ2008, сравнены с каталогом ОН8С08СЯА. Отклонения аномалий силы тяжести показаны на рис. 76 горизонталями с высотой сечения 2 мгал. На рис. 8 показана гистограмма распределения частоты отклонений. Краткая статистика отклонения аномалий силы тяжести: среднее отклонение равно 0,1 мгал, среднеквадратическое отклонение равно 3,1 мгал.

Из рисунка 76 видно, что на суше и на островах, где нет данных альтиметрии, отклонения аномалий силы тяжести по-прежнему больше, чем на море, где есть данные альтиметрии. Если рассмотреть только морской район (широта: 08° - 18°, долгота: 109°30' + 114°, на рисунке ограничен синей линией), то максимальное отклонение равно 8,9 мгал, минимальное отклонение равно -11,3 мгал, среднее отклонение равно -0,1 мгал, и среднеквадратическое отклонение равно 1,8 мгал.

Рис. 8. Гистограмма распределения частоты отклонений, полученных при сравнении с каталогом О^СШ

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

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

О результатах гравиметрической съёмки па побережье Вьетнама Наибольший интерес, конечно, представляет сравнение полученных нами аномалий силы тяжести с непосредственно измеренными значениями на побережье. Мы располагали такими данными в 133 хаотично расположенных точках (см. рис. 9). Из них 123 точки удалось получить в BGI (International Gravimetric Bureau, Франция) и 10 точек получены из работы ЦНИИГАиК в 2011 г. по созданию государственной гравиметрической сети Вьетнама. Все эти точки находятся на островах и на побережье не далее 30 км от моря.

Сравнение с результатами гравиметрической съёмки аномалий силы тяжести, полученных по данным альтиметрии при использовании модели EGM96 Аномалии силы тяжести на 133 точках интерполированы методом коллокации из усреднённых аномалий силы тяжести, полученных по данным альтиметрии девяти циклов спутника ENVISAT при использовании EGM96 и сравнены с аномалиями силы тяжести в свободном воздухе, полученными по результатам гравиметрической съёмки на побережье Вьетнама.

Краткая статистика результатов сравнения: среднее отклонение равно -12,0 мгал; среднеквадратическое отклонение равно 15,2 мгал. На рис. 10 показана гистограмма распределения частоты отклонений.

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

Разность (мгал)

Рис. 10. Гистограмма распределения частоты отклонений

Сравнение с результатами гравиметрической съёмки аномалий силы тяжести, полученных по данным альтиметрии при использовании модели ЕС,М2008

Аномалии силы тяжести на 133 точках снова интерполированы методом кол-локации из усреднённых аномалий силы тяжести, полученных по данным альтиметрии девяти циклов спутника ЕМУ^АТ при использовании ЕСМ2008, и сравнены с аномалиями силы тяжести в свободном воздухе, полученными по результатам гравиметрической съёмки на побережье. Краткая статистика результатов сравнения: среднее отклонение равно -3,3 мгал; среднеквадратическое отклонение равно 8,7 мгал. На рис. 11 показана гистограмма распределения частоты отклонений.

-30 -20 -10 0 10 20 30

Разность (мгал)

Рис. 11. Гистограмма распределения частоты отклонений

Из сравнения аномалий силы тяжести, полученных по данным альтиметрии, с результатами гравиметрической съёмки видно, что на акватории, прилегающей к Вьетнаму, точность аномалий силы тяжести, полученных по данным альтиметрии при использовании ЕОМ2008 гораздо лучше, чем при использовании ЕвМ96.

6. Результаты совместного уравнивания аномалий силы тяжести, полученных по результатам гравиметрической съёмки на побережье и по результатам альтиметрии на акватории

Аномалии силы тяжести, полученные по результатам гравиметрической съёмки, имеют, конечно, более высокую точность, чем вычисленные по данным альтиметрии. Поэтому целесообразно выполнить совместное уравнивание аномалий силы тяжести, полученных по результатам гравиметрической съёмки на побережье и по результатам альтиметрии на акватории. Для расчётов использованы аномалии силы тяжести, полученные по данным альтиметрии спутника Е1МУ18АТ девяти циклов при использовании ЕОМ2008 для процедуры «удаления - восстановления», и аномалии силы тяжести в свободном воздухе, полученные по результатам гравиметрической съёмки на побережье в 133 точках. Предварительно у всех аномалий силы тяжести удалена низкочастотная часть поля с помощью модели ЕОМ96. Процесс уравнивания выполнен методом среднеквадратической коллокации.

В результате, после восстановления низкочастотной части поля получены уточненные аномалии силы тяжести на регулярной сетке 10' х 10'. Общее количество точек 4675, максимальная аномалия силы тяжести равна 132,2 мгал, минимальное значение равно -84,5 мгал, среднее значение равно -4,2 мгал, среднеквадратиче-ское значение равно 21,0 мгал. Полученные значения аномалий силы тяжести представлены на рис. 12 горизонталями с высотой сечения 10 мгал.

Сравнение аномалий силы тяжести до и после уравнивания

Сравнены 4675 значений аномалии силы тяжести до и после уравнивания. Максимальная разность оказалась равной 13,9 мгал, минимальная -13.1 мгап, сред-22

неквадратическое отличие равно 2,0 мгал, но среднее изменение равно нулю. Эти изменения показаны на рис. 13 горизонталями с высотой сечения 1 мгал.

Рис. 12. Уравненные аномалии силы Рис. 13. Разности между значениями

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

Вьетнаму нивания

Из рис. 13 видно, что после уравнивания аномалии силы тяжести изменились только в районе, где есть результаты гравиметрической съёмки, и в ближайших окрестностях съёмки. И это понятно, так как чем дальше от района съёмки, тем меньше сказывается её влияние. Количественной мерой при этом служит радиус корреляции поля, который в данном эксперименте равен 0.4 градуса.

ЗАКЛЮЧЕНИЕ

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

1 - собран необходимый материал в виде десяти циклов альтиметрии ЕЫУ18АТ (за период с 20 июля 2009 г по 9 августа 2010 г) и выполнен отбор данных, необходимых для работы на акватории, прилегающей к Вьетнаму;

2 - разработан алгоритм определения аномалий силы тяжести по данным спутниковой альтиметрии с учётом формы акватории Вьетнама и наличия островов;

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

4 - оптимизирована методика ковариационного анализа аномалии силы тяжести на акватории и показано преимущество использования метода коллокации и программного пакета GRAVSOFT при обработке альтиметрии, в частности, для совместного уравнивания результатов обработки с доступными данными гравиметрии;

5 - разработанные алгоритмы и фортран-программы автоматизируют весь процесс обработки результатов спутниковой альтиметрии ENVISAT и доступны для оперативного повторения с другими исходными данными;

6 - сформирована регулярная сетка с шагом 10' значений аномалии силы тяжести, полученная на прилегающей к Вьетнаму акватории по результатам альтиметрии ENVISAT;

7 - даны рекомендации по работе с современными форматами космической геодезии с помощью программного пакета BRAT.

ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ

1. Нейман Ю.М, Нгуен Ван Шанг. Определение аномалии силы тяжести на акватории Вьетнама по данным альтиметрии ENVISAT // Изв. вузов. «Геодезия и Аэрофотосъемка». -2011. -№ 5. -С.15 - 21.

2. Нгуен Ван Шанг. Уравнивание пересечений треков спутниковой альтиметрии на акватории, прилегающей к Вьетнаму // Изв. вузов. «Геодезия и Аэрофотосъемка».-2012.-№ 3-С.8 - 12.

3. Нгуен Ван Шанг. Определение высот морской поверхности по данным ENVISAT на акватории Вьетнама // Научно-технический журнал ХГГУ. Ханойский горногеологический университет. -2011. -№ 35. -С.81 - 85.

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

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

СОДЕРЖАНИЕ.

СПИСОК ОБОЗНАЧЕНИЙ.

ВВЕДЕНИЕ.

ГЛАВА 1. ОБЩИЕ СВЕДЕНИЯ ОБ АЛЬТИМЕТРИИ И ОПРЕДЕЛЕНИИ ВЫСОТ МОРСКОЙ ПОВЕРХНОСТИ ПО РЕЗУЛЬТАТАМ ПРОЕКТА ENVISAT.

1.1. КРАТКИЕ СВЕДЕНИЯ ОБ АЛЬТИМЕТРИИ.

1.1.1. Историческое развитее альтиметрии.

1.1.2. Принцип альтиметрии.

1.1.3. Типы спутников альтиметрии и их параметры.

1.1.4. Некоторые области применения альтиметрии.

1.2. ПРОЕКТ ENVISAT.

1.2.1. О спутнике ENVISAT.

1.2.2. Инструменты и приборы спутника ENVISAT.

1.3.0 РАБОТЕ С ПРОГРАММНЫМ КОМПЛЕКСОМ BRAT.

1.3.1. Структура пакета BRAT.

1.3.2. Чтение данных.

1.4. ОСНОВНЫЕ ФУНКЦИИ ПАКЕТА BRAT ПРИМЕНИТЕЛЬНО К АЛЬТИМЕТРИИ НА АКВАТОРИИ, ПРИЛЕГАЮЩЕЙ К ВЬЕТНАМУ.

1.4.1. Представление данных спутниковой альтиметрии.

1.4.2. Вычисление высоты морской поверхности.

1.4.3. Выбор необходимых данных и вычисление аномалии уровня океана.

1.4.4. Представление полного цикла.

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

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

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

1.5.1. Описание результатов.

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

ГЛАВА 2. ОПРЕДЕЛЕНИЕ АНОМАЛИИ СИЛЫ ТЯЖЕСТИ ПО ДАННЫМ АЛТИМЕТРИИ ЕМУ18АТ.

2.1. ПРИНЦИП «УДАЛЕНИЯ - ВОССТАНОВЛЕНИЯ» ПРИ ОБРАБОТКЕ СПУТНИКОВОЙ АЛЬТИМЕТРИИ.

2.1.1. Принцип «удаления - восстановления».

2.1.2. Модель ЕвМ96.

2.1.3. Модель ЕвМ2008.

2.2. ДИНАМИЧЕСКАЯ ВЫСОТА МОРЯ. ОПРЕДЕЛЕНИЕ ВЫСОТЫ СРЕДНЕЙ ДИНАМИЧЕСКОЙ ТОПОГРАФИИ МОРЯ.

2.2.1. Динамическая высота моря.

2.2.2. Определение высот средней динамической топографии моря

2.3. УРАВНИВАНИЕ ПЕРЕСЕЧЕНИЙ ТРЕКОВ СПУТНИКОВОЙ АЛЬТИМЕТРИИ В ЛОКАЛЬНОМ РАЙОНЕ.

2.3.1. Определение местоположения и разности высот пересечения.

2.3.2. Уравнивание пересечений треков.

2.4. ВЫЧИСЛЕНИЕ АНОМАЛИИ СИЛЫ ТЯЖЕСТИ ПО ДАННЫМ АЛЬТИМЕТРИИ МЕТОДОМ КОЛЛОКАЦИИ.

2.4.1. Вычисление аномалии силы тяжести по данным альтиметрии как локальная задача коллокации.

2.4.2. Определение эмпирической ковариационной функции.

2.4.3. Решение уравнений коллокации.

2.5. УТОЧНЕНИЕ АНОМАЛИЙ СИЛЫ ТЯЖЕСТИ ПО РЕЗУЛЬТАТАМ ГРАМЕТРИЧЕСКОЙ СЪЁМКИ.

2.6. ОБЩАЯ СХЕМА ОПРЕДЕЛЕНИЯ АНОМАЛИЙ СИЛЫ ТЯЖЕСТИ ПО ДАННЫМ АЛЬТИМЕТРИИ ЕЫУ^АТ.

ГЛАВА 3. ПРАКТИЧЕСКИЕ РЕЗУЛЬТАТЫ ВЫЧИСЛЕНИЯ АНОМАЛИЙ СИЛЫ ТЯЖЕСТИ НА АКВАТОРИИ, ПРИЛЕГАЮЩЕЙ К ВЬЕТНАМУ, ПО ДАННЫМ АЛЬТИМЕТРИИ ЕЫУ^АТ.

3.1. РЕЗУЛЬТАТЫ ВЫЧИСЛЕНИЯ АНОМАЛИЙ СИЛЫ ТЯЖЕСТИ ПРИ ИСПОЛЬЗОВАНИИ ЕСМ96 ДЛЯ ПРОЦЕДУРЫ «УДАЛЕНИЯ-ВОССТАНОВЛЕНИЯ»

3.1.1. Результаты удаления высот геоида, соответствующих ЕОМ96.

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

3.1.3. Результаты уравнивания пересечений треков.

3.1.4. Результат определения эмпирической и аналитической ковариационной функции.

3.1.5. Результаты оценки остаточных аномалий силы тяжести и восстановления низкочастотной части поля.

3.1.6. Усреднение аномалий силы тяжести для десяти циклов.

3.1.7. Усреднение аномалий силы тяжести для девяти циклов.

3.2. РЕЗУЛЬТАТЫ ВЫЧИСЛЕНИЯ АНОМАЛИЙ СИЛЫ ТЯЖЕСТИ ПРИ ИСПОЛЬЗОВАНИИ ЕОМ2008 ДЛЯ ПРОЦЕДУРЫ «УДАЛЕНИЯ - ВОССТАНОВЛЕНИЯ».

3.2.1. Удаление высот геоида, соответствующих модели ЕОМ2008.

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

3.2.3. Результат уравнивания пересечений треков.

3.2.4. Результаты определения эмпирической и аналитической ковариационной функции.

3.2.5. Результаты оценки остаточных аномалий силы тяжести и восстановления низкочастотной части поля.

3.2.6. Усреднение аномалий силы тяжести для десяти циклов.

3.2.7. Усреднение аномалий силы тяжести для девяти циклов.

3.2.8. Сравнение аномалий силы тяжести, полученных при использовании для процедуры «удаления - восстановления» моделей ЕвМ96 и ЕОМ2008.

3.3. СРАВНЕНИЕ ВЫЧИСЛЕННЫХ АНОМАЛИЙ СИЛЫ ТЯЖЕСТИ С КАТАЛОГОМ ВЫ8С080ЫА.

3.3.1. О каталоге Б^ОШНА.

3.3.2. Сравнение с данными каталога ОК8С080ЫА аномалий силы тяжести, вычисленных при использовании для процедуры «удаления -восстановления» модели ЕвМ96.

3.3.3. Сравнение с данными каталога ОН8С080ЯА аномалий силы тяжести, вычисленных при использовании для процедуры «удаления -восстановления» модели ЕвМ2008.

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

3.4.1. О результатах гравиметрической съёмки на побережье.

3.4.2. Сравнение с результатами гравиметрической съёмки аномалий силы тяжести, полученных по данным альтиметрии при использовании модели ЕвМ96.

3.5. РЕЗУЛЬТАТЫ СОВМЕСТНОГО УРАВНИВАНИЯ АНОМАЛИЙ СИЛЫ ТЯЖЕСТИ, ПОЛУЧЕННЫХ ПО РЕЗУЛЬТАТАМ ГРАВИМЕТРИЧЕСКОЙ СЪЁМКИ НА ПОБЕРЕЖЬЕ И ПО РЕЗУЛЬТАТАМ АЛЬТИМЕТРИИ НА АКВАТОРИИ.

3.5.1. Ковариационный анализ аномалии силы тяжести.

3.5.2. Результаты уравнивания.

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

Обоснование актуальности темы диссертации

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

Цель и основные задачи исследований

3) выделение динамических возмущений морской поверхности путём уравнивания пересечений спутниковых треков;

4) исследование степени гладкости гравитационного поля на изучаемой территории;

6) исследование точности полученных результатов и исследование возможностей их уточнения по доступным результатам гравиметрической съёмки.

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

- учёт реальной специфики акватории,

- учёт реальной гладкости соответствующего гравитационного поля,

- ослабление влияния динамической топографии,

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

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

Положения, выносимые на защиту

1 - Общий алгоритм определения аномалии силы тяжести по результатам альтиметрии спутника ЕМУГБАТ применительно к акватории Вьетнама и с учётом её специфики.

4 - Регулярная сетка с шагом 10' значений аномалии силы тяжести, полученная на прилегающей к Вьетнаму акватории по результатам альтиметрии ЕМУ18АТ.

Вклад автора в проведенное исследование

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

Практическая значимость работы

1 - Разработанные алгоритмы и соответствующее программное обеспечение на языке Фортран автоматизируют весь процесс обработки результатов спутниковой альтиметрии ЕМУГБАТ и доступны для оперативного повторения с другими исходными данными.

2 - На прилегающей к Вьетнаму акватории по результатам альтиметрии ЕЫУ18АТ получены конкретные значения аномалии силы тяжести на регулярной сетке с шагом 10'.

3 - Даны рекомендации по работе с современными форматами космической геодезии при использовании программного пакета BRAT. Структура и объем диссертации

Диссертация состоит из введения, трёх глав основного текста и заключения. Список литературы насчитывает 46 наименований на русском, вьетнамском и английском языках. Кроме того, имеется список использованных сокращений, список принятых обозначений и 06 приложений. Общий объем работы составляет 157 страницы машинописного текста и содержит 71 рисунок и 66 таблиц.

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

ЗАКЛЮЧЕНИЕ

1 - собран необходимый материал в виде десяти циклов альтиметрии ENVISAT (за период с 20 июля 2009 г по 9 августа 2010 г) и выполнен отбор данных, необходимых для работы на акватории, прилегающей к Вьетнаму;

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

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

2. Бернхард Гофман-Велленгхоф, Гельмут Мориц. Физическая геодезия.: Перевод с английского Ю. М. Нейман, Л. С. Сугаипова/Под редакцией Ю.М. Неймана. -М.: Изд-во МИИГАиК, 2007. 426 с.

3. Коллектив авторов ответственный редактор Бровар Б. В. Гравиметрия и геодезия. М.: Научный мир, 2010. - 572 с.

4. Нгуен Ван Шанг. Уравнивание пересечений треков спутниковой альтиметрии на акватории, прилегающей к Вьетнаму // Изв. вузов. «Геодезия и Аэрофотосъемка». -2012. -№ 3. -С.8 12.

5. Нейман Ю.М, Нгуен Ван Шанг. Определение аномалии силы тяжести на акватории Вьетнама по данным альтиметрии ЕгмБа! // Изв. вузов. «Геодезия и Аэрофотосьемка». -2011. -№ 5. -С.15 21.

6. Нейман Ю. М., Вергасов В. А., Журкин И. Г., Красикова М. В., Смирном С. А. Вычислительная математика. М.: Недра, 1976. - 230 с.

7. Нейман. Ю. М. Вариационный метод физической геодезии и коллока-ция. М.: Научный мир, 2010. - 50 с.

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

9. Огородова Л. В., Шимбирев Б. П., Юзефонвич А. П. Гравиметрия. М.: Недра, 1978.-325 с.

10. Ю.Пигин А. П., Березина С. В. Глобальная модель геоида ЕвМ2008 предварительный анализ. Г. Минск, 2009. 6 с.

11. СканЭкс. Сайт http://www.spacestudy.ru

12. Шимбирев Б. П. Теория фигуры Земли. М.: Недра, 1975. - 432 с.

13. Юзефонвич А. П., Огородова Л. В. Гравиметрия: Учебник для вузов. -М.: Недра, 1980.-320 с.

14. Яблонский JI. Н. Научно-технический отчет с результатами измерений по контракту № 3/7 09 // Разработка проекта и производство гравиметрических измерений на территории Вьетнама. Москва, 2011. - 42 с.

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

16. Alexander Braun and Shum С.К. Satellite Altimetry Database and Processing. Byrd Polar Reseach Center and laboratory for Space Geodesy and Remote Sensing, The Ohio State University, Columbus, OH, USA. 2004. - 6 P

17. Andersen О. B, Knudsen P, Philippa A. M. Berry. The DNSC08GRA global marine gravity field from double retracked satellite altimetry // Journal of Geodesy. 2010. -№ 3. -P.191 - 199.

18. Andersen O.B. Marine Gravity and Geoid from Satellite Altimetry. Geodetic Department, DTU Space, Juliane Maries Vej 30, DK - 2100, Denmark. 2010.-50 p.

19. Andersen. B, Knudsen. P. The DNSC08MSS global Mean Sea Surface. Danish national Space Center, Denmark. 2008. 16p.

20. Andersen. О. B, Knudsen. P, Berry. P, Kenyon. S. The DNSC08 ocean-wide altimetry derived gravity anomaly field. Danish national Space Center, Denmark. 2008. 14p.

21. AVISO. DT CorSSH and DT SLA Product Handbook. Toulouse France. 2010.-22 p.

22. AVISO. Website http://www.aviso.oceanobs.com

23. Cheinway Hwang, Barry Parsons. Gravity anomalies derived from Seasat, Geosat, ERS-1 and T/P altimetry and ship gravity: a case study over the Reykjanes Ridge // Geophys. J. Int. 1995. - № 122. - P. 551 - 568.

24. ESA. GUT Installation Guide, February 2009. 5 p.

25. ESA. GUT Tutorial, March 2011. 80 p.

26. ESA. GUT User Guide and Algorithm Descriptions, February 2009. 97 p.

27. Jing Xie. Implementation of Parallel Least-Squares Alorithms for Gravity Field Estimation. Dep. of Civil and Environmental Engineering and Geodetic Science, The Ohio State University, Columbus, Ohio, Report № 474, 2005. -65 p.

28. Lee-Lueng Fu, Anny Cazenave. Satellite Altimetry and Earth Sciences. ACADEMIC PRESS, San Diego San Francisco - New York - Boston -London - Sydney -Tokyo. 2001. - 463 p.

29. Moritz H. Stepwise and Sequential Collocation. Dep. of geodetic science, OSU, Report. № 203, 1973.

30. NIMA. Department of Defense World Geodetic System 1984. National Imagery and Mapping Agency, American. 2000.

31. Proefschrift A Decade of ERS Satellite Orbits and Altimetry. Delft University Press, Netherlands. 2002. 219 p.

32. Rene Forsberg, C.C. Tscherning. Geodetic Gravity Field Modelling Programs. National Space Institute and Niels Bohr Institute, University of Copenhagen, Denmark, 2008. 75p.

33. Rosmorduc V. Basic Radar Altimetry Toolbook practical. Bergen, Norway. 2009. 32 p.

34. Sylvain Bonvalot. Annual Report 2007-2008. International Gravimetic Bureau, Toulouse, France, April 2008. 5p.

35. Tcherning С. C. A FORTRAN program for the determination of the anomalous potential using Stepwise least squares collocation. Dep. of geodetic science, OSU, Report № 212,1974.

36. Tscherning С. C. Gravity field modelling with GRAVSOFT least-squares collocation // Lecture Notes, International School for the Determination and Use of the Geoid, Milano, Oct. 10-15. International Geoid Service. 1994. - P. 101-134.

37. Veronique Amans, Henri Laur. Access to Envisat data. European Space Agency. 2007. 64 p.

Информация о работе

Нгуен Ван Шанг
кандидата технических наук
Москва, 2012
ВАК 25.00.32

Диссертация

Гравиметрическое обеспечение прилегающей к Вьетнаму акватории по результатам спутниковой альтиметрии - тема диссертации по наукам о земле, скачайте бесплатно

Автореферат

Похожие работы