Лазерная локация космических аппаратов в последнее время широко используется для высокоточного определения положений космических аппаратов в околоземном пространстве. Однако существующие математические модели обработки результатов лазерных измерений не учитывают воздействия на распространение световых сигналов вращения Земли и ее гравитационного поля. Поэтому при проведении лазерной локации космических аппаратов, находящихся на околоземных орбитах, и при геодезических измерениях неявно предполагается, что лазерная станция, покоящаяся на поверхности вращающейся Земли, является инерциальной системой отсчета. В результате световые лучи, входящие в оптические приборы и выходящие из них, при математическом моделировании считаются прямыми линиями и распространение световых импульсов по этим лучам предполагается равномерным.
Однако эти приборы находятся в слабонеинерциальной системе отсчета, связанной с вращающейся Землей, и их лучи подвержены действию слабого гравитационного поля Земли. Поэтому согласно общей теории относительности Эйнштейна [1] в гравитационном поле Земли и при действии сил инерции вращающейся системы отсчета лучи должны искривляться, а движение световых импульсов по этим лучам должно быть неравномерным. И хотя эти релятивистские эффекты в гравитационном поле вращающейся Земли малы, с определенного уровня точности измерений они будут искажать наблюдательные данные.
В настоящее время в стадии подготовки находятся несколько космических программ по проведению прецизионных астрометрических измерений с помощью интерферометров, установленных на космических аппаратах. Для баллистико — навигационного обеспечения этих программ локацию участвующих в экспериментах космических аппаратов необходимо производить с точностью нескольких сантиметров по дальности и порядка 1 угловой секунды по угловому положению, что доступно только лазерным станциям. На таком уровне точности уже необходимо производить учет воздействия гравитационного поля Земли и неинерциальности системы отсчета лазерной станции, находящейся на вращающейся Земле, на распространение световых импульсов при лазерной локации космических аппаратов. Таким образом, без математической модели, учитывающей уравнения и соотношения общей теории относительности, осуществлять с требуемой точностью прогнозирование и контроль орбит космических аппаратов в процессе выполнения прецизионных космических программ невозможно. Однако в математическом плане эти вопросы были неразработанными.
Более того, при проведении лазерной локации космических аппаратов был замечен довольно странный эффект [2]: на тех лазерных станциях, которые использовали маломощные лазеры с относительно узким лазерным пучком (как известно, лазерный пучек не является цилиндром, а в силу ряда физических причин [3] представляет собой расширяющийся конус), было очень трудно попасть электромагнитным импульсом в движущийся космический аппарат. Обычно космический аппарат, движущийся по низкой околоземной орбите, достаточно хорошо виден в телескоп лазерной станции и оператор, вводя упреждение вдоль вектора его скорости, должен был попасть лазерным импульсом в отражатели, установленные на его внешних поверхностях.
Однако было замечено, что если лазерный импульс направить в расчетное место его встречи с космическим аппаратом, то этот импульс не попадет на отражателиесли же ось лазера немного довора-чивать и направлять электромагнитный импульс не в расчетное место встречи, а сканировать окрестности этого места, то при некотором небольшом угловом отклонении (составляющем несколько угловых секунд) от расчетного места космический аппарат попадет в лазерный пучек и в телескопе лазерной станции появится отраженный от космического аппарата электромагнитный импульс. После этого дальнейшее сопровождение космического аппарата лазерным пучком переключают на компьютер, который автоматически вводит найденную экспериментально дополнительную угловую поправку, и процесс локации проходит, как правило, без вмешательства оператора.
В качестве первой причины такого небольшого рассогласования между расчетным угловым направлением на место встречи лазерного импульса с космическим аппаратом и фактическим угловым направлением на это место было выдвинуто [2] предположение о несовпадении оптической оси телескопа с осью лазерного луча на лазерной станции. Однако лазерный импульс посылается к космическому аппарату через телескоп лазерной станции и дополнительные исследования точности юстировки показали [2], что эти оси с требуемой точностью совпадают.
На мощных лазерных станциях, использующих широкие лазерные пучки, этот эффект не наблюдался, так как электромагнитный импульс в окрестности космического аппарата имел большие поперечные размеры и отражатели, установленные на космическом аппарате, попадали на периферию этого импульса.
Приступая к исследованию этого вопроса, мы предположили, что причиной такого рассогласования является недостаточно точная математическая модель, использовавшаяся для расчетов угловых координат места встречи лазерного импульса и космического аппарата, не учитывающая эффекты общей теории относительности.
Поэтому целью настоящей диссертации было решение задачи построения математической модели релятивистских поправок, вызываемых воздействием гравитационного поля Земли и полями сил инерции вращающейся Земли, на распространение световых импульсов при лазерной локации космических аппаратов и при геодезических измерениях, а также разработка алгоритмов вычисления искажений, вносимых этими поправками, в наблюдательные данные.
Первая глава носит вводный характер. В ней рассматриваются задачи геодезических измерений и лазерной локации космических аппаратов. Дается описание основных систем отсчета, используемых при построении математической модели лазерной локации космических аппаратов. В § 4, следуя нашей работе [4], изложены основы математической модели слабонеинерциальной системы отсчета в присутствии гравитационного поля.
Вторая глава настоящей диссертации посвящена построению математической модели вращательной аберрации при геодезических измерениях. На основе релятивистского уравнения для изотропной геодезической нами построены уравнения световых лучей во вращающейся системе отсчета и найдены их аналитические решения. Анализ этих решений позволил исследовать ряд релятивистских эффектов во вращающейся системе отсчета.
Первый из них — это то, что распространение электромагнитных сигналов из одной точки пространства в другую точку осуществляется по одному лучу, а обратно — по другому лучу. Выведена формула определяющая угол 5а между касательными к этим лучам в точке излучения импульса.
Второй эффект — это линейное смещение электромагнитного импульса после его отражения от ретрорефлектора (оптического устройства, отражающего импульсы вдоль касательной к лучу, по которому эти импульсы приходят). Выведена формула, определяющая расстояние 811 между точкой излучения электромагнитного импульса и лучем, по которому распространяется отраженный импульс.
На основе полученных формул в § 8 проведено вычисление вращательной аберрации при геодезических измерениях углов на местности. Впервые такие вычисления были проведены в наших работах [5−6]. Однако величина релятивистских эффектов, вызываемых вращением Земли, в условиях геодезических измерений мала: на расстоянии между источником электромагнитного излучения и отражателем, равным 100 км, полученные формулы дают 5а = Ю-2 угловой секунды, 6Я = 0,5 см.
Эти величины меньше требуемой в настоящее время точности геодезических измерений: 5а = 0,1 угловой секунды. Так как при оптических измерениях в геофизике расстояние между источником света и наблюдателем обычно не превышает 100 км, то значение 5а = Ю-2 угловой секунды можно принять в качестве оценки точности угловых измерений, начиная с которой необходимо учитывать угловую вращательную аберрацию.
Следует отметить, что эти релятивистские эффекты из-за их малой величины в эксперименте пока не измерялись. Однако с инструментальной точки зрения в настоящее время имеются технические возможности для измерения указанных выше релятивистских эффектов, вызываемых вращением Земли. Как показано нами в § 9 диссертации, в настоящее время существует экспериментальная возможность провести измерение релятивистских эффектов в земных условиях.
Поэтому проведенное во второй главе математическое моделирование распространения электромагнитных сигналов во вращающейся системе отсчета позволило установить пределы точности геодезических измерений, начиная с которой будет необходимо учитывать релятивистские эффекты, а также предложить эксперименты по проверке теоретических предсказаний о величине этих эффектов.
В третьей главе настоящей диссертации построена математическая модель релятивистских поправок при прецизионной лазерной локации космических аппаратов. Постановка задачи включала в себя следующие существенные моменты: лазерная станция предполагалась расположенной на поверхности вращающейся Земли, а космический аппарат с отражателем электромагнитных импульсов двигался по орбите, ориентированной произвольным образом относительно оси вращения Земли.
При проведении лазерной локации как импульс электромагнитного излучения, так и космический аппарат согласно общей теории относительности движутся в римановом пространстве-времени, метрический тензор которого зависит от гравитационного поля Земли, Луны, Солнца и других планет Солнечной системы и полей сил инерции неинер-циальной системы отсчета.
Безразмерный гравитационный потенциал Земли и очень мал [7]: II ~ Ю-9. Однако именно он дает основной вклад в гравитационное поле в окрестности Земливклад остальных тел Солнечной системы согласно работе [8] существенно меньше. Так как нам было необходимо построить математическую модель лазерной локации, обеспечивающую относительную точность на уровне ~ Ю-10, то вкладами остальных тел Солнечной системы мы пренебрегали.
При вычислении релятивистских эффектов мы также учитывали, что лазерная станция находится в слабонеинерциальной системе отсчета из-за вращения Земли вокруг своей оси с угловой скоростью П = 7,3−10″ 5 рад/сек.
Лазерную локацию космических аппаратов и обработку полученных данных удобно было проводить в топоцентрической системе отсчета, начало которой совпадает с положением лазерной станции на поверхности Земли, ось г направлена по местной вертикали, ось х — по касательной к меридиану, ось у — по касательной к параллели.
Как показал анализ, для достижения при лазерной локации относительной точности на уровне ~ Ю-10 движение космического аппарата можно было рассматривать в ньютоновском приближении, а движение электромагнитного импульса — с учетом постньютоновских слагаемых в метрическом тензоре. Так как в рассматриваемой задаче имелись малые параметры, то для удобства дальнейшего использования формул математической модели вычисление релятивистских поправок было проведено в аналитическом виде методом последовательных приближений.
Следует отметить, что в общем случае в математическую модель лазерной локации космических аппаратов необходимо включать поправки, связанные с влиянием на движение лазерного импульса нестабильной атмосферы, с влиянием на движение космического аппарата муль-типольных гравитационных моментов Земли, а также учитывать ряд аналогичных (см., например, [9−10]) возмущающих факторов. Но все эти поправки к настоящему времени хорошо изучены и включены в существующие математические модели локации космических аппаратов. Так как они создают лишь малые возмущения в движении космических аппаратов и лазерных импульсов, то их можно считать аддитивными поправками, не влияющими на малые величины 5Ь, Дт и 8а.
Все полученные в третьей главе настоящей диссертации результаты являются новыми и опубликованы в нашей работе [И]. В научной литературе аналогичные задачи ранее рассматривались, но либо в недостаточно общем виде [12−13] (подробнее см. конец § 14), либо относились к спутниковой конфигурации, когда лазерная станция находилась не на вращающейся Земле, а на каждом из космических аппаратов [1416], т. е. находилась в системе отсчета, метрический тензор которой имел локально инерциальный вид.
И, наконец, в четвертой главе настоящей диссертации на основе построенной математической модели проведен численный расчет величины релятивистских эффектов при лазерной локации космических аппаратов.
На основе формул, построенной в третьей главе математической модели, можно определять положение космического аппарата в интересующие нас моменты времени, уточнять элементы орбиты и решать некоторые другие задачи. В настоящей диссертации мы основное внимание сосредоточили на одной из этих задач: определили величину релятивистских поправок к дальности 5L, различие в темпах хода часов лазерной станции и бортовых часов космического аппарата Дт, обусловленную различием физических условий, в которых они находятся, а также разность 8а угловых положений направлений, по которым посылается к космическому аппарату и, соответственно, принимается после отражения от него лазерный импульс.
Как следует из полученных в диссертации формул, численные значения 8L, Дт и 8а должны существенно зависеть от элементов орбиты космического аппарата. Поэтому для численных расчетов были выбраны два типа наиболее распространенных орбит: низкая почти круговая орбита с высотой апогея ha = 600 км и высотой перигея hp = 400 км и высокоапогейная с высотой апогея ha = 2 • 105 км и высотой перигея hp = 250 км.
Первый тип орбит характерен для космических аппаратов различного назначения (дистанционное зондирование Земли, испытание космической техники и т. п.), влючая Международную космическую станцию и существовавшие ранее орбитальные станции. Второй тип орбит обычно используется для космических аппаратов научного назначения (преимущественно астрофизического), представителями которого являются, например, готовящиеся к запуску в России аппараты космических программ «Радиоастрон», «Спектр-УФ», «Спектр-РГ», «Озирис» [17−19] и другие.
Результаты численных расчетов 8Ь, Дт и 8а для указанных выше двух типов орбит показали, что для обоих типов орбит релятивистские поправки 8Ь и Ат достигают вполне измеримых величин и их необходимо учитывать при лазерной локации космических аппаратов астрометрического назначения.
Величина же 8 а для низких почти круговых орбит в настоящее время находится на пределе возможностей измерения и ею можно пренебречь. Для высокоапогейных космических аппаратов угол 6а на отдельных участках орбит может достигать значений более 10 угловых секунд. В этом случае существенно необходимо измерять величину 8а и использовать полученные данные для устранения релятивистских искажений при прецизионном определении местоположения космического аппарата.
Проведенный нами в четвертой главе диссертации численный эксперимент с использованием полученных формул (18.2)-(18.6) математической модели лазерной локации позволил обнаружить неизвестный ранее эффект — наличие не равного нулю угла 8а при локации с Земли экваториального геостационарного космического аппарата, т. е. космического аппарата, находящегося на расстоянии около 35 900 км от поверхности Земли, имеющего период обращения, совпадающий с земными сутками. В результате космический аппарат, находящийся на экваториальной геостационарной орбите, все время как бы «висит» над заданной точкой земной поверхности. Результаты численных расчетов показали, что в этом случае угол 8а равен 3,5 угловой секунды и может быть измерен современными лазерными станциями.
Этот результат впервые получен нами и опубликован в работе [11]. Так как во вращающейся системе отсчета Земли геостационарный спутник неподвижен, то его можно рассматривать как часть вращающейся системы отсчета Земли. Поэтому к нему применимы формулы, выведенные нами во второй главе настоящей диссертации. Подстановка в эти формулы параметров орбиты геостационарного спутника, впервые осуществленная в нашей работе [20], показывает, что 8а = 3,5 угловой секунды, а <Ш=0,6 км в полном соответствии с результатами численного расчета. Это совпадение результатов численных расчетов в физически выделенном случае геостационарного спутника, для которого существует аналитическое решение, дает основание для утверждения о достоверности построенной нами в третьей и четвертой главах диссертации математической модели лазерной локации космических аппаратов.
В заключении перечислены основные результаты, полученные в настоящей диссертации и выносимые на защиту.
Основные результаты диссертации докладывались на XII Российской гравитационной конференции (20−26 июня 2005 г., Казань, КГУ), Международной конференции «Тихонов и современная математика'1 (19−25 июня 2006 г., Москва, МГУ), Пятом Международном аэрокосмическом конгрессе (27−31 августа 2006 г., Москва, МГУ), на научных семинарах МАТИ-РГТУ им. К. Э. Циолковского и опубликованы в семи работах [4−6, 11, 20, 36, 73], три из которых являются статьями в рецензируемых журналах, включенных в список ВАК.
ЗАКЛЮЧЕНИЕ
.
Перечислим основные результаты, составляющие содержание настоящей диссертации и выносимые на защиту.
1. Построена математическая модель распространения световых сигналов во вращающейся системе отсчета. На основе этой модели исследованы основные оптические эффекты, вызываемые действием полей инерции вращающейся системы отсчета. Предложен эксперимент по проверке предсказания общей теории относительности о движении электромагнитных импульсов во вращающейся системе отсчета.
2. Получена формула для угла вращательной аберрации при геодезических измерениях. Установлены пределы точности геодезических измерений, начиная с которых необходимо учитывать влияние вращения Земли на результаты геодезических измерений.
3. Построена математическая модель лазерной локации космических аппаратов, учитывающая воздействие на электромагнитные импульсы гравитационного поля Земли и вращения лазерной станции вместе с Землей вокруг ее оси. Используя существенные уравнения этой модели, вычислены релятивистские поправки к дальности 8Ь. На основе численного анализа показано, что для низких почти круговых орбит 8Ь не превышает 10 см, а для высокоапогейных орбит 8Ь > 20 см.
4. Проведено моделирование разности между темпом хода часов, находящихся на лазерной станции, и темпом хода бортовых часов космического аппарата согласно уравнениям релятивистской теории. Показано, что для низких почти круговых орбит бортовые часы отстают от часов лазерной станции на несколько мкс за один оборот космического аппарата по орбите, а в случае высокоапогейных орбит бортовые часы уходят вперед по сравнению с часами лазерной станции на несколько сотен мкс за один оборот космического аппарата по орбите. Проведена оценка точности, с какой необходимо контролировать этот эффект при проведении космических астрометрических программ.
5. На основе построенной математической модели лазерной локации космических аппаратов проведено вычисление угла 5а между направлением, по которому электромагнитный импульс излучается лазерной станцией к космическому аппарату, и направлением, по которому отраженный от космического аппарата импульс регистрируется на лазерной станции. Показано, что для низких почти круговых орбит угол ¿-а незначительно превышает одну угловую секунду, а для высокоапогейных орбит — доходит до 14 угловых секунд.
6. Проведенный нами численный эксперимент с использованием построенной математической модели лазерной локации позволил обнаружить неизвестный ранее эффект — наличие не равного нулю угла 8а при локации с Земли экваториального геостационарного космического аппарата, т. е. космического аппарата, находящегося на расстоянии около 35 900 км от поверхности Земли, имеющего период обращения, совпадающий с земными сутками. В результате космический аппарат, находящийся на экваториальной геостационарной орбите, все время как бы «висит» над заданной точкой земной поверхности. Результаты численных расчетов показали, что в этом случае угол 8а равен 3,5 угловой секунды и может быть измерен современной лазерной станцией.
7. Проведено вычисление угла аберрации Sip, создаваемого полями сил инерции вращающейся системы отсчета и гравитационным полем Земли. Этот угол возникает между истинным направлением на космический аппарат в момент прихода к нему лазерного импульса и направлением, по которому на космический аппарат посылается электромагнитный импульс. В результате численных расчетов было установлено, что с относительной точностью ~ 10~5 угол аберрации Sip равен половине угла 5а между касательной к лучу, по которой электромагнитный импульс был отправлен к космическому аппарату, и касательной к лучу, по которому он возвратился на лазерную станцию.
Результаты диссертационной работы были использованы при выполнении НИОКР по договору № 4/2006 от 01.02.2006 г. между Федеральным государственным унитарным предприятием «Научно — производственное объединение имени С.А.Лавочкина» и Институтом прикладных и теоретических исследований, по договору № 1/2006 от 10.01.2007 г. между Федеральным государственным унитарным предприятием «ЦНИИМАШ» и Институтом прикладных и теоретических исследований, а также в учебном процессе при чтении лекций по курсу «Математическое моделирование» студентам факультета № 2 «Аэрокосмические конструкции и технологии» ГОУ ВПО «МАШ'-РГТУ им. К. Э. Циолковского, о чем имеются соответствующие акты внедрения.
В заключение автор выражает глубокую благодарность научному руководителю доктору технических наук, профессору Чебурахину Игорю Федоровичу за внимание к работе и полезные советы, а также всему коллективу кафедры «Кибернетика» за создание атмосферы, способствовавшей успешной работе над диссертацией.