Задача авиационной гравиметрии с использованием градиентометрических измерений тема автореферата и диссертации по механике, 01.02.01 ВАК РФ

ЗМ^'.-- - ■ I

йа-:-- ■ |

Ж'

И-Ь-:- ■

ЗЙЙ^г&'С

аейя«'

ЭЙВИЙЯ: I

т0т

иЫЙ)®«-ЙЭаЗЙ?..

£>04

У

■{ а ■ к ,

д/

Московский ордена Ленина, ордена Октябрьской Революции и ордена Трудового Красного Знамени государственный университет им. М. В. Ломоносова

Механико-Математический факультет

¡Я»,.-...'

ЙЭ :

Шй'.-Г;---;

Й-^'-'г

онш

ИР®5* I

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

Папуша Ирина Анатольевна

Задача авиационной гравиметрии с использованием градиентометрических измерений

¿¡й!

Специальность 01.02.01 — "Теоретическая механика'

Диссертация на соискание ученой степени кандидата физико-математических наук

Научный руководитель —

доктор физико-математических наук

профессор Н. А. Парусников

ЙЁЬ-Ц.- . !

С^ге-..........|

* - : ■ ■ |

Ее^-' ■ I

ШШМ-, ~ ' ,-мзй:;

Москва — 1999

Оглавление

I

1г

Введение 4

Глава 1. Алгоритмы решения задачи авиагравиметрии

с градиентометром 13

1.1. Исходные соотношения градиентометрии....... 13

1.2. Формирование тензора гравитационного градиента по измерениям градиентометра, установленного на подвижном основании ................. 16

1.3. Гравитационные градиентометры ........... 18

1.4. Уравнения ошибок..................... 30

1.5. Алгоритмы решения задачи оценки гравитационных аномалий..............................33

Глава 2. Математические модели погрешностей канонического гравитационного градиентометра 46

2.1. Математическое описание акселерометра и его ошибок .................................48

2.2. Модель погрешностей канонического градиентометра ............................. 51

2.3. Стохастическая мера оцениваемости.......... 59

2.4. Исследование оцениваемости инструментальных погрешностей градиентометра .............. 61

Глава 3. Анализ точности решения задачи оценки гравитационных аномалий с использованием градиен-тометрических измерений 69

3.1. Процедуры предварительной обработка гравиизме-рений............................ 71

3.2. Результаты численного моделирования воздействия вибровозмущений на точность решения задачи авиационной гравиметрии .................. 74

Заключение 87

Библиография 89

Приложение 1 94

Приложение 2 98

а:

зЬяййй

Г

к;

■■■ I. глг,

Ж;

Й

Введение

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

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

Е'П-.' |

;; 1 Величина силы тяжести в основном зависит от распределения

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

можно выделить следующие:

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

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

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

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

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

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

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

ционной системы (ИНС). Однако помимо полезной информации о силе тяжести, показания гравиметра включают в себя составляющие, порожденные силами инерции, возникающими при движении ЛА. Эти силы могут быть исключены при помощи информации о параметрах движения объекта. Такую информацию получают от ИНС и дополнительных источников, например, высокоточной спутниковой навигационной системы (СНС), высотомеров, системы измерения вертикальных скоростей. Методы инерциаль-ной и спутниковой навигации хорошо известны [1] - [4].

Достаточно полное отражение современного состояния авиагравиметрии за рубежом содержится в [5], [6]. Из отечественных разработок следует отметить проводимую при участии лаборатории управления и навигации механико-математического факультета МГУ им. Ломоносова, МИЭА и Института физики Земли. Создан экспериментальный образец комплекса, в состав которого входят двухкомпонентная ИНС, блок гравиметров, два приемника (базовый и бортовой) СНС, работающей в дифференциальном режиме, система математической обработки гравиметрической информации. Образец прошел три серии летных испытаний, подтвердивших его вполне удовлетворительные качества [7].

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

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

8'зй:.

Ш

а, -л- ■

а.......

Ц

Officii-;.

' же время принцип действия другого прибора — градиентометра основан на измерении изменения гравитации по разности кажущихся ускорений в двух точках. Так как ускорение сил инерции

| ; в любой точке стабилизированнои системы координат одинаково,

grg-; I можно разделить инерционные и гравитационные эффекты, изме-

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

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

Интерес к вопросам градиентометрии возник еще в 70-е годы. Начальные исследования носили теоретический характер, в них си-стематезируется теория градиента гравитационного поля. Среди первых молено выделить работу [8], где приводится теоретическое решение задачи определения поля сил тяготения и сил инерции переносного движения в связанных с движущимся телом осях посредством инерциальных измерений. Выводятся общие соотношения и исследуется вопрос о составе измерений, необходимых для решения данной задачи в космосе.

В начале 80-х разработчики приступают к проектированию конструкций градиентометров. Создаются демонстрационные модели, отрабатывается элементная база. В статье [9] рассмотрены три разрабатываемых модели градиентометра на подвижном основании фирм Bell Aerospace, the Draper Laboratory и Hughes Aircraft. Описан принцип получения полезного сигнала для вращающихся конструкций и поплавковой камеры. Найдена оптимальная относительная ориентация чувствительных датчиков и количественный состав измерителей для определения всех компонент тензора гравитационного градиента. Эти три схемы приборных реализаций стали впоследствии классическими.

При переходе от теоретических исследований к практическим

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

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

Ш'

II

'}■:

¡Ь, дованы возможные источники ошибок таких приборов. В [13] на-

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

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

Некоторые варианты конструкций градиентометра предлагаются в [14], [15], также исследуется возможность их приборной реализации.

Работы [16], [17] затрагивают вопрос об использовании колебательной системы крутильного типа, установленной на подвесе, для создания гравитационного градиентометра на подвижном основании.

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

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

расчетов погрешностей калибровки чувствительности.

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

Алгоритм предварительной подготовки данных для фильтрации сигнала градиентометра предлагается в [21]. Предусмотрены процедуры перекодировки, упорядочения, отбраковки и усреднения по множеству реперных точек экспериментальных данных.

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

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

Появившийся в настоящее время большой интерес к авиагравиметрическим исследованиям, в том числе с использованием градиентометров, отчасти объясняется значительным повышением качества гравиизмерительной аппаратуры. Открылись перспективы применения принципиально новых технологий (эффекта сверхпроводимости при низких температурах). Наличие на борту ЛА при-

5 ~ ■ if-

¿¡--Г,

it'-

¡i ■:;.

r

3 •-••

Ш

fci: ¡¿ад-

SiisTi _

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

Из числа недавних публикаций по градиентометрии научно-популярная статья [24] содержит исторические сведения и современные концепции. Сделаны общие выводы по результатам морских испытаний в Мексиканском заливе системы с градиентометром фирмы Bell Aerospace.

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

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

'"'¿''¿г

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

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

Тр