Борис Шустов - Астероидно-кометная опасность - вчера, сегодня, завтра

После сближения Апофиса с Землей в 2029 г. точность прогнозирования его движения резко падает. Если же для части возможных траекторий после 2029 г. будет иметь место тесное сближение (не соударение!) в 2036 г., то дальнейшее движение по этим траекториям становится практически недетерминированным. Другими словами, современной точности начальных данных недостаточно для исследования методом численного интегрирования возможных соударений в период, последующий за 2036 г., если в этом году будет иметь место тесное сближение астероида с Землей. По мере получения новых наблюдений точность орбиты будет, конечно, расти, однако в ближайшие несколько лет вряд ли удастся достаточно точно предвычислить траектории тесных сближений в 2036 г., чтобы надежно проанализировать все возможные сближения и возможные столкновения в последующие годы.

Интуитивно ясно, что опасными траекториями, на которые может перейти астероид в результате тесного сближения с Землей, являются резонансные траектории невысокого порядка. Например, если после сближения астероид окажется на орбите, которая имеет соизмеримость средних движений 1:1 с Землей, то ровно через год астероид и Земля окажутся в том же месте и возможно их новое тесное сближение или соударение. Из-за расширения трубки возможных траекторий в результате взаимодействия с Землей существует больше возможностей попадания в такие опасные резонансы. Соответствующие области начальных условий называют зонами резонансного возврата. Так, в работе [Chesley, 2005] рассмотрены возможные резонансные возвраты астероида Апофис после сближения в апреле 2029 г. Ко времени написания работы Чизли орбита астероида имела еще сравнительно низкую точность, и в результате рассеяния трубки возможных траекторий после этого сближения астероид мог оказаться на различных резонансных орбитах, указанных на рис. 10.13 в главе 10. Эти орбиты приводили к опасным сближениям Апофиса с Землей в 2034, 2035, 2036, 2037, 2038, 2046 и 2048 гг. Однако по мере увеличения точности орбиты число возможных резонансных возвратов быстро уменьшалось, и уже к августу 2005 г. остался лишь один возможный резонансный возврат в апреле 2036 г. Существенно в этой связи, что в ближайшие несколько лет до 2012 г. точность орбиты астероида Апофис не будет расти в результате поступления новых наблюдений, поскольку в это время наблюдать его с Земли практически невозможно. После появления новых наблюдений и уточнения орбиты Апофиса в 2012–2013 гг. совершенно аналогично уменьшится число возможных траекторий с опасными сближениями в период после 2036 г., которые детально обсуждаются в работе [Соколов и др., 2008].

В первом приближении резонансные возвраты можно исследовать, используя прием, называемый методом точечных гравитационных сфер (ТГС). Для нескольких последовательных сближений с планетой этот метод позволяет с умеренной трудоемкостью получить картину возможных преобразований орбит. При этом, если сближения достаточно тесные и относительная скорость достаточно велика, метод дает приемлемую точность результатов. Точность определяется из сравнения с результатами численного интегрирования неупрощенных уравнений движения. Метод ТГС позволяет в ряде важных случаев рассматривать вместо сложной задачи трех тел последовательность простых задач двух тел. Этот прием в различных модификациях часто применяется в небесной механике и механике космического полета, в частности для проектирования траекторий космических аппаратов с несколькими гравитационными маневрами. Некоторые подробности, а также литературные ссылки можно найти, например, в работах [Соколов и др., 2008; Елькин и др., 2003].

В методе ТГС исследуемая траектория аппроксимируется последовательностью кеплеровых гелиоцентрических орбит соударения с планетами. При «соударении» происходит мгновенное преобразование скорости — поворот вектора планетоцентрической скорости на угол между асимптотами планетоцентрической гиперболы. Модуль планетоцентрической скорости остается при этом постоянным. Элементы кеплеровой орбиты преобразуются в этом случае по формулам классической задачи двух тел. Как известно, по координатам и скоростям в данный момент времени элементы орбиты определяются однозначно. Угол поворота вектора планетоцентрической скорости зависит от прицельного расстояния и в рамках рассматриваемой аппроксимации (сжатия гравитационной сферы в точку, обнуления прицельного расстояния) не определен, т. е. может быть произвольным в некоторых пределах.