Аннотация результатов, полученных в 2019 году
Околоземное космическое пространство (ОКП) все интенсивнее используется в интересах хозяйственной деятельности Человечества. Создаются и разворачиваются в космосе все новые и новые спутниковые системы. Отработавшие объекты этих систем, как правило, остаются на околоземных орбитах и превращаются в космический мусор. В настоящее время в околоземном пространстве по оценкам Европейского космического агентства находится свыше 30 тысяч объектов космического мусора размером более 10 см, порядка 900 тысяч объектов размером в диапазоне от 1 до 10 см, порядка 128 миллионов — около 1 мм и это число постоянно увеличивается. На орбитальную динамику этих объектов оказывают влияние множество различных факторов, которые способны значительно изменять как положение орбиты в пространстве, так и саму орбиту. К таким факторам относятся вековые и орбитальные резонансы. Знание резонансной структуры орбитального пространства позволяет с большой достоверностью прогнозировать возможные катастрофические события, оптимизировать выбор областей для размещения новых спутниковых систем и паркинга отработавших объектов. Исследование динамической структуры околоземного орбитального пространства требует проведения обширных численных экспериментов, сопровождаемых анализом больших объемов данных. В связи с этим важными задачами проекта являются внедрение автоматизации выполнения ряда этапов исследований, требующих большой обработки данных, в том числе с помощью методов машинного обучения и искусственных нейронных сетей (ИНС), а также разработка и усовершенствование методик и программного обеспечения для исследования динамики околоземных объектов искусственного происхождения. Что и было осуществлено в первую очередь в ходе выполнения первого этапа проекта в соответствии с календарным планом. В ходе выполнения первого этапа проекта были проведены исследования различных методов применения ИНС для выявления наличия резонансов в ОКП и определения их типа (устойчивые/неустойчивые). В результате была получена обученная ИНС, которая показала точность 99.3% при определении наличия резонанса и его классификации. Был применен алгоритм автоматической кластеризации данных HDBSCAN, который позволил быстро локализовать объекты, демонстрирующие нетипичную и интересную динамику. Разработано программно-алгоритмическое обеспечение для выявления эффекта Лидова–Козаи и автоматизирован отбор объектов с хаотической динамикой по параметру MEGNO. Проведена модификация программного комплекса «Численная модель движения систем ИСЗ», реализованного на кластере «СКИФ Cyberia» Томского Государственного Университета. Производительность программного комплекса повышена благодаря замене интегратора на новый и более эффективный интегратор, разработанный в НИИ ПММ ТГУ, и смене алгоритма распараллеливания. Данная модификация численной модели позволила повысить быстродействие расчетов почти в 3 раза при сохранении точности вычислений. К важным результатам первого этапа по выполнению проекта можно отнести модификацию методики выявления и исследования вековых резонансов в динамике околоземных объектов. Особенность методики в том, что в отличие от традиционно используемых приближенных аналитических методик все частоты в движении спутника, возникающие под действием возмущающих факторов, определяются численно по точным формулам. Такой подход дает более достоверные результаты, чем чисто аналитическая методика при исследовании резонансной структуры движения по орбитам с большим или быстро растущим эксцентриситетом. С помощью разработанного и модифицированного в ходе выполнения проекта программного обеспечения, в том числе с использованием ИНС, построены карты динамической структуры зон низкоорбитальной части ОКП (где уже не так сильно действие атмосферы, и имеет смысл рассматривать влияние вековых резонансов) и части зоны среднеорбитального пространства до 25500 км по большой полуоси (до зоны функционирования системы ГЛОНАСС). Выявлены области, где на движение объектов оказывают влияние только вековые резонансы. Было показано, что особенности орбитальной эволюции объектов в зоне, не содержащей орбитальных резонансов, определяются, прежде всего, действием вековых апсидально-нодальных резонансов первого порядка: апсидального резонанса Лидова–Козаи и нодального резонанса, связывающего скорости изменения долгот восходящего узла спутника и возмущающих тел. Выявлено, что нодальный резонанс первого порядка проходит через всю рассматриваемую область и проявляет устойчивость при наклонениях 0, 90 и 180 градусов. Результаты исследования показали, что в рассмотренной области орбитального пространства впервые становится заметным действие механизма Лидова–Козаи. Особенности его влияния на динамику объектов зависят от взаимодействия этого механизма с влиянием сжатия Земли. Выявлено, что механизм Лидова–Козаи начинает проявляться на орбитах с большими полуосями около 10000 км и затем проявляется во всех областях зоны МЕО. Кроме того, показано, что совместное действие нескольких вековых резонансов встречается нечасто в динамике объектов рассматриваемой зоны. Наложение нескольких устойчивых вековых резонансов разных типов не приводит к возникновению хаотичности в движении объектов, а при совместном действии устойчивых и неустойчивых вековых резонансов наблюдается рост хаотичности движения. Проведено сопоставление особенностей эволюции объектов нерезонансной зоны (не содержащей орбитальных резонансов) с эволюцией объектов, движущихся в зонах орбитального резонанса. Выявлено, что движение в резонансных зонах отличается большей хаотичностью. Если в нерезонансной зоне явление хаотичности движения наблюдается редко, то для резонансных зон оно является характерным свойством. Причем, определяющим фактором в возникновении хаотичности в движении объектов является наличие в динамике неустойчивых компонент орбитального резонанса. Вклад вековых резонансов в возникновение хаотичности в резонансных зонах второстепенен, однако они оказывают влияние на орбитальную эволюцию, которое проявляется ростом амплитуд долгопериодических колебаний позиционных переменных. Это относится, прежде всего, к вековым резонансам первого порядка, о которых говорилось выше. Наряду с космическим мусором опасность для Земли могут представлять и объекты естественного происхождения, а именно астероиды. Астероиды, сближающиеся с Землей, обладают сложным движением вследствие наличия сближений с большими планетами и возможных неустойчивых орбитальных и вековых резонансов. Исследование динамики таких объектов осуществимо только численными методами, поскольку требует учета влияния множества сил, в частности негравитационных эффектов, таких как световое давление и эффект Ярковского. Изучение вероятностной орбитальной эволюции требует применения многопроцессорных вычислительных систем и обработки больших объемов информации. В ходе выполнения первого этапа проекта был осуществлен переход к автоматизированной обработке результатов вычислений. В частности, автоматизирован процесс построения графиков с использованием скриптов и их анализ с помощью методов машинного обучения. Модифицирован программный комплекс «ИДА», предназначенный для исследования динамики астероидов, путем включения алгоритмов определения коэффициента трансверсального ускорения А2, создаваемого эффектом Ярковского в движении астероидов. Реализовано и проанализировано два метода получения этого параметра из условия минимума среднеквадратической ошибки представления наблюдений. В первом методе перебираются различные значения A2 в пределах заданного интервала с определенным шагом. При каждом значении параметра решается задача наименьших квадратов и определяется среднеквадратическая ошибка представления наблюдений. За искомое значение параметра трансверсального ускорения А2 берется такое, которое соответствует минимальному значению среднеквадратичной ошибки. Данная методика имеет ряд недостатков: она ограничивает исследователя интервалом перебора параметра А2 и дискретным шагом, и требует решения задачи наименьших квадратов на каждом шаге, что существенно увеличивает время работы программы. Во втором методе параметр включается в число оцениваемых параметров при решении задачи наименьших квадратов. Здесь задача решается однократно, и позволяет, кроме того, получить ошибку определения параметра А2. Разработанные алгоритмы были протестированы на нескольких астероидах, имеющих малые перигелийные расстояния, так как эффект Ярковского может оказывать значительное влияние на их динамику из-за регулярного прохождения объектов вблизи Солнца. В частности, проведено сравнение методик при определении трансверсального ускорения А2 для астероидов, имеющих малые перигелийные расстояния и различный состав наблюдений. В связи с тем, что часть астероидов имеет радарные наблюдения, в программный комплекс «ИДА» была включена возможность их обработки. На примере астероидов 425755 2011 CP4 и 3200 Phaethon (следует отметить, что эти два астероида имеют радарные наблюдения) показана возможность получать требуемые параметры с помощью разработанного программного обеспечения, с точностью сопоставимой с точностью NASA. Таким образом, на основе проведенных исследований можно сделать вывод о большей эффективности второй методики. Далее была детально исследована орбитальная эволюция астероидов 504181 2006 TC и 2007 PR10 c учетом и без учета эффекта Ярковского. Выявлены сближения 2006 TC и 2007 PR10 со всеми планетами внутренней группы (Меркурий, Венера, Земля, Марс), рассмотрена вековая резонансная динамика (выявлен ряд апсидально-нодальных резонансов) и получена оценка времени предсказуемости движения астероидов с помощью параметра OMEGNO. Результаты исследований показали, что эффект Ярковского оказывает существенное влияние на большую полуось орбит астероидов, что приводит к изменениям в числе и расстоянии сближений. Что касается оценки времени предсказуемости с помощью параметра OMEGNO, то движение астероида 504181 2006 TC можно считать регулярным лишь на интервале времени (1440, 2270) гг., а 2007 PR10 – на интервале (1390, 4310) гг. За пределами этих интервалов у орбит проявляются признаки хаотичности, к которым, скорее всего, приводят многократные и тесные сближения астероидов с большими планетами. Кроме того, в программном комплексе «ИДА» реализована возможность определения параметра светового давления из наблюдений. Методика протестирована на примере астероида 3200 Phaethon. Полученные нами результаты согласуются с данными NASA. В качестве подготовки к работе по второму этапу построена упрощенная модель численного движения спутников Луны, которая включает в себя основной потенциал (LP150Q) и точечную модель влияния внешних тел (Земли и Солнца). Усложнение модели и доведение ее до уровня высокоточной планируется на втором этапе проекта.

Аннотация результатов, полученных в 2020 году
В течение второго этапа работы по проекту проведено исследование области околоземного космического пространства (ОКП) вплоть до сферы действия Луны по отношению к Земле (от 35000 до 315000 км по большой полуоси). Подробно изучена резонансная структура данного околоземного космического пространства. Показано, что значительная часть объектов рассматриваемой области имеет малое время жизни, что вполне объясняется ростом эксцентриситетов орбит. Большой и быстрый рост эксцентриситета и малое время жизни имеют все приполярные объекты, орбиты которых имеют наклонение к экватору, близкое к 90 градусам. Начиная с больших полуосей, равных 235000 км, значительная часть объектов имеет малое время жизни для большинства наклонений, а для больших полуосей, превосходящих 275000 км, все объекты имеют быстро растущие эксцентриситеты и малое время жизни независимо от наклонения. Это связано с прямым действием Луны, причем влияние Луны тем больше, чем ближе орбита объекта к экватору. Характер эволюции объектов, не подверженных прямому действию Луны, объясняется действием многочисленных вековых резонансов, как устойчивых, так и неустойчивых. Были построены карты распределения вековых апсидально-нодальных резонансов, связанных с Луной и с Солнцем в рассматриваемой области ОКП. Данные карты, а также, составленные на их основе карты наложений резонансов различных типов для Луны и для Солнца были получены с использованием искусственных нейронных сетей (ИНС). Анализ всех построенных карт позволил сделать следующие выводы: все устойчивые апсидально-нодальные резонансы, содержащие скорость изменения долготы перицентра, группируются в диапазоне наклонений 60–120 градусов, а устойчивые нодальные резонансы, содержащие скорость изменения долготы восходящего узла, группируются в окрестности трех наклонений 0, 90 и 180 градусов, причем большая часть в окрестности 90 градусов. Это объясняет, почему все объекты приполярной области имеют быстро растущий эксцентриситет и малое время жизни. Для выявления в данной области ОКП отличительных особенностей в движении околоземных объектов в областях, подверженных влиянию орбитальных резонансов, и вне их действия рассмотрены несколько областей орбитальных резонансов. Ранее при исследовании областей действия орбитальных резонансов низких порядков мы приходили к выводу, что области хаотичности движения определяются наложением прежде всего устойчивых и неустойчивых компонент орбитального резонанса, а наложение вековых резонансов только усиливает хаотизацию. Однако при рассмотрении областей резонансов высоких порядков, например 2:1 со скоростью вращения Земли, мы увидели ситуация прямо противоположную той, которую мы наблюдаем в области орбитальных резонансов низких порядков. Сопоставляя данные по всем видам резонансов, действующих в области резонанса 2:1 со скоростью вращения Земли, мы обнаружили, что орбитальный резонанс проявляется лишь незначительным ростом хаотичности движения объектов, попадающих в этот резонанс, в то время как основным фактором оказывающим влияние на орбитальную эволюцию являются вековые резонансы. Таким образом, в данной части ОКП динамика области, содержащей орбитальные резонансы высоких порядков, мало чем отличается по своим особенностям от областей, не содержащих данных резонансов. Разработана усовершенствованная за счет уточнения модели сил и использования более эффективного интегратора численная модель движения искусственных спутников Луны. Кроме того численная модель была дополнена возможностью вычисления осредненного параметра MEGNO, позволяющего судить о хаотичности движения космических объектов. Модель реализована в двух версиях: для персонального компьютера и в среде параллельных вычислений на кластере «СКИФ Cyberia» Томского Государственного Университета. Проведено всестороннее тестирование модели, которое подтвердило ее высокую точность. Проведено тестирование системы автоматической обработки данных с использованием искусственной нейронной сети путем подбора размеченных данных по окололунным объектам. Тестирование показало, что для корректного распознания наличия или отсутствия резонанса, необходимо помимо временных рядов, связанных с критическим аргументом, рассматривать в качестве сигналов еще два временных ряда, связанных с изменением эксцентриситета и наклонения орбиты. С применением разработанного программного обеспечения выполнен анализ влияния структуры селенопотенциала на движения окололунных объектов, который показал, что причиной роста эксцентриситета на низких орбитах является прямое действие селенопотенциала. Проведен анализ резонансной структуры окололунного орбитального пространства. По результатам исследования составлены карты действия на окололунные объекты вековых апсидально-нодалтьных резонансов 1-4 порядков и полувековых резонансов 2-5 порядков со средними движениями возмущающих тел Земли и Солнца. Показано, что главной причиной роста эксцентриситета и значительного сокращения времени жизни на орбите для объектов на средних и больших высотах следует считать влияние на орбитальную динамику ИСЛ векового резонанса типа Лидова–Козаи и близких к нему по типу резонансов. Полувековые резонансы со средними движениями Земли и Солнца локализованы в средних и верхних областях окололунного орбитального пространства, начиная с большой полуоси в 6000 км. Действие этих резонансов имеет неустойчивый характер, в изменении резонансных аргументов прослеживается перманентный переход от либрации к циркуляции и обратно. Поэтому значительного воздействия на орбитальную эволюцию они не оказывают. Класс астероидов, имеющих малые перигелийные расстояния (менее 0.15 а.е.), представляет большой интерес для ученых с точки зрения возможной опасности для Земли, исходящей от этих объектов. На динамику астероидов при прохождении их вблизи Солнца могут оказывать существенное влияние негравитационные возмущающие факторы, такие, как эффект Ярковского и световое давление. Не учет этих возмущений в модели сил в ряде случаев может привести к недостоверному прогнозу движения, в том числе к ошибочным оценкам вероятности столкновения. Основная проблема учета этих возмущений заключается в том, что для большинства астероидов в настоящее время плохо известны физические свойства, необходимые для определения параметров эффекта Ярковского и светового давления. В данной ситуации одним из способов определения параметров является метод получения их из наблюдений астероидов. В ходе выполнения второго этапа проекта с использованием модифицированного программного комплекса «ИДА» были получены параметры эффекта Ярковского А2, светового давления ks и их среднеквадратические ошибки для всех астероидов с малыми перигелийными расстояниями. По данным центра малых планет на январь 2021 года известно 50 таких астероидов. В процессе улучшения их орбит параметры А2 и ks включались в число оцениваемых, наряду с координатами и компонентами скорости. Наблюдения были взяты с сайта центра малых планет. Оценка точности определения параметра светового давления ks показала, что, для большинства астероидов среднеквадратическая ошибка ks совпадает по порядку с самим значением, а для 8 недавно открытых астероидов порядок ошибки превышает таковой для ks. Только для 3200 Phaethon погрешность на порядок меньше значения параметра. На примере астероидов с известными физическими параметрами (3200 Phaethon, 137924 2000 BD19, 276033 2002 AJ129, 386454 2008 XM, 394130 2006 HY51, 394392 2007 EP88) показано, что использование улучшенного значения ks несущественно влияет на значение среднеквадратической ошибки представления наблюдений, что, по-видимому, обусловлено малым влиянием светового давления на движение изучаемых астероидов. Оценка влияния светового давления на движение этих объектов позволила сделать вывод, что световое давление влияет меньше, чем сжатие Солнца и релятивистские эффекты от Солнца. Полученные значения параметра эффекта Ярковского и его среднеквадратической ошибки для астероидов с малыми перигелийными расстояниями позволили сделать следующие выводы. Подобно коэффициенту светового давления ks значения параметра А2 и его среднеквадратической ошибки для большинства объектов имеют один порядок. Анализ результатов определения параметра А2 для исследуемых астероидов показал, что с увеличением наблюдательного интервала точность определяемого параметра значительно улучшается. Данный вывод подтвердился и в случае с экспериментами с модельными и реальными наблюдениями: на примере астероидов 3200 Phaethon и 137924 2000 BD19 был проведен эксперимент по уменьшению числа реальных наблюдений от настоящего времени в прошлое, а на примере астероидов 2008 MG1 и 2020 BU13 – эксперимент по увеличению мерного интервала путем использования модельных наблюдений. Однако для каждого астероида есть предельный мерный интервал, при достижении которого точность определения параметра перестает улучшаться. Также анализ результатов показал, для объектов с длиной мерного интервала меньше 30 суток определение значения параметра эффекта Ярковского из наблюдений является недостаточно надежным и для его применения стоит дождаться новых наблюдений астероидов. С использованием программного комплекса «ИДА» построена вероятностная орбитальная эволюция астероидов 3200 Phaethon и 2007 PR10 с учетом эффекта Ярковского и без него. Для обоих астероидов выявлены многократные и тесные сближения с Меркурием, Венерой, Землей и Марсом и получена оценка хаотичности их орбит с помощью параметра OMEGNO (Orthogonal Mean Exponential Growth factor of Nearby Orbit). Учет эффекта Ярковского в обоих случаях привел к значительным изменениям в эволюции большой полуоси, что отразилось на количестве сближений астероидов с планетами и расстоянии до них. Вероятностные области движения астероидов в результате таких изменений ведут себя по-разному без учета эффекта и с ним. Поиск резонансов в движении астероидов показал, что 3200 Phaethon движется в окрестности орбитального резонанса 3/7 с Венерой, и без учета эффекта Ярковского астероид дольше находится в резонансе с планетой, что говорит о необходимости его учета для более правильной картины. Резонансов в динамике астероида 2007 PR10 выявлено не было.

Аннотация результатов, полученных в 2021 году
Исследование динамической структуры околоземного и окололунного орбитальных пространств является весьма актуальной задачей, поскольку знание динамических особенностей различных областей космического пространства необходимо и при разработке новых спутниковых систем, и при определении областей захоронения отработавших объектов. В рамках настоящего проекта было проведено исследование динамической структуры околоземного космического пространства (ОЗКП), начиная от 8000 км по большой полуоси и до сферы действия Луны, а также окололунного космического пространства (ОЛКП) от 1911.8 км до 26070 км по большой полуоси. В обоих случаях наклонения варьировались от 0 до 180 град. Моделирование орбитальной эволюции объектов данных областей проводилось при помощи усовершенствованного и разработанного в ходе выполнения предыдущих этапов проекта программного обеспечения: численных моделей движения систем ИСЗ и ИСЛ, а также с помощью методов машинного анализа, в том числе обученной искусственной нейронной сети. Были построены карты распределения вековых апсидально-нодальных резонансов, связанных с Луной и с Солнцем в ОЗКП и с Землей и Солнцем в ОЛКП. Анализ орбитальной эволюции всей совокупности рассмотренных объектов (более 100 тыс.) и карт резонансов позволил сделать ряд общих выводов о характере движения объектов изучаемых областей. Значительная часть объектов имеет малое время жизни, что связано с ростом эксцентриситетов орбит. Большой и быстрый рост эксцентриситета и малое время жизни имеют все приполярные объекты, орбиты которых имеют наклонение к экватору, близкое к 90 град. Это объясняется действием вековых резонансов и справедливо и для ОЗКП, и для ОЛКП. Все устойчивые апсидально-нодальные резонансы, содержащие скорость изменения долготы перицентра, группируются в диапазоне наклонений 60–120 град. В ОЗКП устойчивые нодальные резонансы, содержащие скорость изменения долготы восходящего узла, группируются в окрестности трех наклонений 0, 90 и 180 град., причем большая часть в окрестности 90 град. В ОЛКП данный тип резонансов в основном сосредоточен в диапазоне от 70 до 110 град. по наклонению. Наложение большого количества апсидально-нодальных резонансов приводит к стремительному росту эксцентриситетов орбит объектов, вследствие чего объекты прекращают свое существование либо сталкиваясь с Землей/Луной, либо уходя из сферы их тяготения. Резонансы, связанные со средним движением возмущающих тел Луны и Солнца, появляются только в нижней части рассматриваемой области ОЗКП. Вековые резонансы со средним движением Солнца распределены в диапазоне больших полуосей от 8000 км до 21000 км по обе стороны от наклонения 90 град., а вековые резонансы со средним движением Луны проявляются только на орбитах, большие полуоси которых менее 16000 км и только в динамике объектов с обратным движением. Резонансы, связанные со средним движением возмущающих тел Земли и Солнца (2-5 порядка) в ОЛКП не обнаружены. Помимо влияния резонансов на движение ИСЗ и ИСЛ действуют и иные факторы, приводящие к уменьшению времени жизни. Так в ОЗКП, начиная с больших полуосей, равных 235000 км, значительная часть объектов имеет малое время жизни для большинства наклонений, а для больших полуосей, превосходящих 275000 км, все объекты имеют быстро растущие эксцентриситеты и малое время жизни независимо от наклонения. Это связано с прямым действием Луны, причем влияние Луны тем больше, чем ближе орбита объекта к экватору. В ОЛКП при рассмотрении динамики низколетящих ИСЛ исследования показали, что главным источником роста эксцентриситета на низких орбитах следует считать прямое действие сложного гравитационного поля Луны. С использованием для ИСЗ идей, предложенных Т. Галардо (2006) для астероидов, была получена сравнительная оценка сил всех действующих в этом пространстве орбитальных резонансов. Эта оценка показала, что чем ниже объект, тем сила резонанса больше и наоборот, чем он выше, тем сила меньше. Анализ всей совокупности полученных данных позволяет сделать вывод, что на действие практически всех устойчивых компонент орбитального резонанса налагаются действия либо неустойчивых компоненты этого же резонанса, либо вековой резонанс, что приводит к хаотичности. Причем зоны усиления хаотичности совпадают с наложением вековых резонансов. В рамках данного этапа проведено также исследование динамических структур ОЗКП и ОЛКП с учетом влияния светового давления (СД) для различных значений парусности (отношением площади миделева сечения к массе). Параметры численного моделирования и способ распределения объектов аналогичен эксперименту, проведенному без СД. Сопоставление карт резонансов, построенных без учета СД, и с учетом СД показывает, что в большинстве случаев резонансные зоны сохраняются, хотя для ряда резонансов при росте парусности зоны могут, как становиться шире, так и совсем исчезать. Кроме того возрастание величины парусности приводит к увеличению скорости роста эксцентриситетов орбит и, как следствие, к уменьшению времени жизни. Исследование резонансных движений астероидов, сближающихся с Землей, является одной из важных задач на пути к решению проблемы астероидной опасности. Поскольку устойчивый орбитальный резонанс может служить защитным механизмом от тесных сближений с планетами, а неустойчивый повышает риск сближения астероида с планетой, что в свою очередь может привести к значительным изменениям элементов орбит объектов и возможным сближениям с Землей. Особенно актуально исследование резонансов средних движений для астероидов, имеющих малые перигелийные расстояния, так как орбиты этих объектов имеют вытянутую форму и потенциально способны взаимодействовать с большинством планет Солнечной системы. В процессе исследования выявлены все астероиды с малыми перигелийными расстояниями, движущиеся в окрестности устойчивых и неустойчивых резонансов средних движений с большими планетами. Так как эффект Ярковского (ЭЯ) может оказывать существенное влияние на движение этих астероидов при прохождении их вблизи Солнца, рассмотрено влияние эффекта на поведение резонансных характеристик со временем. СД также может менять значение большой полуоси, поэтому отдельно было исследовано его влияние на резонансные характеристики. Оценено влияние ЭЯ на резонансные характеристики: резонансную щель и критический аргумент. Результаты исследования показали, что влияние ЭЯ на устойчивые резонансные соотношения незначительно: оно может несущественно менять амплитуду либрации, но не приводит к разрушению резонанса. В случае же неустойчивого резонанса изменяется число проходов через точную соизмеримость, а для некоторых объектов резонанс становится даже более устойчивым. Оценка влияния СД на поведение характеристик орбитального резонанса астероидов показала, что СД меньше влияет на движение изучаемых астероидов, чем ЭЯ, особенно при диаметре больше 1 км. Для устойчивого резонанса влияние первого имеет сходный характер с влиянием ЭЯ, а при неустойчивом резонансе учет СД приводит к изменению числа прохождений через точную соизмеримость, в некоторых случаях частично разрушая участки устойчивого резонансного движения или наоборот увеличивая интервалы резонансного движения. Проведено подробное исследование орбитальной эволюции астероидов 137924 2000 BD19 и 394130 2006 HY51 без учета ЭЯ и СД и отдельно с каждым из них. Результаты исследования показали, что учет СД и ЭЯ практически не оказывает влияния на эволюцию орбитальных элементов обоих астероидов, лишь незначительно изменяя поведение большой полуоси. Такой результат влияния отражается на взаимном расположении астероидов и планет в процессе эволюции и приводит к изменению количества сближений с большими планетами. Так как исследуемые возмущения не оказывают влияния на долготу восходящего узла и аргумент перицентра орбит астероидов, на основе которых вычисляются характеристики вековых резонансов, изменений в поведении вековых резонансов выявлено не было. Эволюция параметра OMEGNO с учетом ЭЯ и СД показала, что данные возмущения не влияют на интервал предсказуемости движения астероидов 137924 2000 BD19 и 394130 2006 HY51.