Дифференциальные спутниковые измерения и методы их обработки

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

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

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

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

Перспективы достижения качественно новых результатов решения задачи определения гравитационных полей Земли и планет в настоящее время связываются с так называемыми дифференциальными методами, использующими искусственный спутник (или систему спутников) в качестве пробной гравиметрической массы, слежение за движением которой доставляет необходимые детальные сведения о структуре внешнего гравитационного поля исследуемого небесного тела. Речь идет о методе межспутникового слежения (Satellite-to-Satellite Tracking - SST), предполагающего выполнение бортовых измерений характеристик относительного движения в системе двух или более аппаратов, и методе спутниковой градиентометрии (Satellite Gravity Gradiometry - SGG), основанного на бортовых измерениях вторых производных гравитационного потенциала.

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

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

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

Заметим в этой связи, что в современной литературе принято выделять два подхода к оцениванию параметров гравитационного потенциала исследуемого небесного тела по спутниковым данным: пространственный – space-wise(SW) и временной – time-wise(TW).

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

потенциала однородной сферической Земли и возмущающего потенциала , моделируемого суммой ряда по объемным сферическим функциям до степени :

Согласно идеологии SW-подхода запишем результат некоторого измерения, выполненного в околоземном пространстве в виде функционала

зависящего от – геоцентрического вектора положения и скорости одного (или нескольких) ИСЗ, участвующих в измерении , и – возмущающего потенциала, описываемого набором гармонических коэффициентов ряда (2). Вектор состояния , где – эвклидовое пространство размерности p, возмущающий потенциал , в реальной ситуации описываемый конечным числом гармонических коэффициентов, , где – эвклидовое пространство размерности q. Тогда функционал отображает произведение пространств и в пространство действительных чисел (т.е. на числовую ось):

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

Линеаризуя (3), приходим к уравнению поправок:

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

В рамках TW-подхода, который вернее было бы назвать орбитальным, возмущающий потенциал (2) или непосредственно измеряемая его трансформанта представляются функцией времени , играющего в данном контексте роль углового расстояния ИСЗ (в плоскости спутниковой орбиты) от его положения в начальный момент времени :

где

,

- функция наклонения орбиты ИСЗ, – долгота ее восходящего узла, – аргумент широты.

Коэффициенты ряда (5), получившие название «связанных» (lumped) коэффициентов, представляют собой линейную комбинацию гармонических коэффициентов модели геопотенциала того же порядка , но различных степеней . При вычислении коэффициентов по формулам (5) суммирование по индексу при начинается со значения , поскольку для ряда (2) имеем . Кроме того, суммирование по производится при этом с шагом равным 2, то есть так, чтобы четность значения индекса всегда совпадала с четностью значения индекса (если четное, суммирование по ведется только по четным индексам и наоборот). Параметры при фиксированном значении наклонения орбиты являются лишь функциями времени (в силу ) и иногда называются коэффициентами чувствительности. Для различных трансформант гравитационного потенциала разложение (5) сохраняет свой вид, меняется только вид коэффициентов чувствительности.

Разложение (5) есть ни что иное, как ряд Фурье с коэффициентами , представляющих собой амплитуды тригонометрических функций угловых параметров. Числовые значения «связанных» коэффициентов могут быть найдены по данным бортовых измерений с помощью соответствующих стандартных вычислительных процедур. На следующем этапе полученные значения коэффициентов рассматриваются как результат псевдоизмерений, позволяющих сформировать вектор свободных членов уравнений поправок, линейных относительно искомых поправок в параметры модели гравитационного поля Земли.

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

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

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