Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

FORECASTING OF THE DEADLINE «ENVIRONMENTAL LIFE» OF THE NEAR-EARTH ORBITS IN TERMS OF PROGRESSIVE PROPAGATION OF SPACE DEBRIS

Shestopalova O.L. 1 Shestopalov R.P. 2 Darbuzova K.O. 2
1 A Branch «Voskhod» of the Moscow aviation institute (national research university)
2 Lyceum «International space school named after. V.N. Chelomey»
The article deals with the issues of forecasting of the deadline «environmental life» of the near-earth orbits in terms of progressive propagation of space debris. We have identified three stages of development of the propagation of space debris: the stage of accumulation of critical density, the stage of cascade propagation of space debris and the stage of stabilization of the space debris density. Forecasting is based on logistic curve for growth of space debris density. For the first time are obtained the analytical expression, that allows to predict the deadline «environmental life» of the near-earth orbits as a solution of fuzzy equations, including fuzzy value of maximum allowable density of space debris. The solution method of this equation is described. Method deals with the triangular shape of the membership function of a fuzzy set values of maximum permissible density of space debris. Shows how to convert a fuzzy result of the prediction for interval form.
the near-earth space
space debris
space debris propagation
the near-earth orbits deadline «environmental life»
forecasting
fuzzy set
1. Veniaminov S.S., Chervonov A.M. Kosmicheskiy musor ugroza chelovechestvu [Space Debris – a Threat to Humanity]. Moskow, IKI RAN, 2012. 189 p.
2. Zade L.A. Ponyatie lingvisticheskoy peremennoy i ego primenenie k prinyatiyu priblizhennykh resheniy [The Concept of Linguistic Variable and its Application to Decision-Making]. Moskow, Mir, 1976. 165 p.
3. Chetyrkin E.M. Statisticheskie metody prognozirovaniya [Statistical Forecasting Methods]. Moskow, Statistika, 1977. 200 p.
4. Mironov A.N., Mironov E.A., Shestopalova O.L., Platonov S.A. Informatsiya i kosmos – Information and Space, 2015, no 3, pp. 167–171.
5. Mironov A.N., Mironov E.A., Shestopalova O.L., Platonov S.A. Fundamentalnye issledovaniya [Fundamental Research]. 2015, no. 2–13, pp. 2815–2818.
6. Mironov A.N., Mironov E.A., Shestopalova O.L., Platonov S.A. Fundamentalnye issledovaniya [Fundamental Research]. 2015, no. 2–13, pp. 2819–2823.
7. Kessler D., Cour-Palais B. Collisional Frequency of Artificial Satellites: The creation of a debris belt // J. Geophysical Research. 1978. Vol. 83. N. A6. 1978.

Одной из негативных тенденций, осложняющей дальнейшее освоение и использование околоземного космического пространства (ОКП) в интересах населения Земли, является тенденция его загрязнения техногенным мусором или так называемым, космическим мусором (КМ) [1]. На протяжении всей космической эры, начиная с 1957 года по настоящее время, степень засоренности ОКП неуклонно возрастала и будет продолжать расти, даже если полностью прекратить запуски космических аппаратов (КА). Этому в определенной мере способствует «подключение» нового механизма формирования КМ, связанного с «размножением» КМ за счет разрушения крупных фрагментов мусора при их неупругом столкновении или столкновении с действующими КА. Этот эффект (впервые описанный Дональдом Кесслером в 1978 году [7]) получил название эффекта каскадного размножения Кесслера. Механизм каскадного размножения Кесслера запускается при достижении плотностью КМ критического значения, после чего количество фрагментов увеличивается по экспоненте. Со временем процесс размножения КМ замедляется по причине уменьшения числа относительно крупных фрагментов мусора и связанного с этим преобладания упругого столкновения над неупругим.

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

Можно выделить три стадии формирования космического мусора: I – стадию накопления критической плотности КМ, II – стадию каскадного размножения КМ и III – стадию стабилизации плотности КМ.

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

В зависимости от стадии скорость Vρ(t) возрастания плотности ОКМ ρ(t) во времени изменяется нелинейно, т.е. сначала возрастает, затем уменьшается до нуля. В то же время для удельной скорости shestopal01.wmf, равной shestopal02.wmf, характерен линейный вид зависимости. Удельная скорость роста плотности ОКМ, равная в начальный момент r, по мере роста плотности ОКМ снижается, и при ρ(t) = P становится равной нулю, т.е. процесс размножения ОКМ останавливается, а плотность ОКМ стабилизируется, асимптотически приближаясь к числу P. Эта зависимость математически может быть записана в виде уравнения

shestopal03.wmf (1)

где shestopal04.wmf, где r – начальное значение удельной скорости возрастания плотности ОКМ, фиксируемое на момент начала первой стадии; P – максимальное установившееся значение плотности ОКМ по завершении третьей стадии

Приняв во внимание, что

shestopal05.wmf,

можно записать

shestopal06.wmf (2)

или

shestopal07.wmf (3)

Уравнение (3) носит название уравнения Ферхюльста, впервые примененного бельгийским математиком П. Ферхюльстом в 1838 году для описания роста численности населения.

Решив дифференциальное уравнение (3), получаем логистическую кривую [3].

shestopal08.wmf (4)

где P0 – начальное значение плотности ОКМ, фиксируемое на момент t0 начала первой стадии.

Можно найти координату точки перегиба функции (4):

shestopal09.wmf

и значение функции (4) в этой точке shestopal10.wmf.

Введем параметр PПред, равный предельному значению плотности ОКМ, при достижении которого использование околоземной орбиты становится невозможным по критериям опасности или нерентабельности затрат на защиту КА от ОКМ.

Предельным сроком «экологической жизни» TС.Ж околоземной орбиты будем называть интервал времени от начала ее использования до достижения предельного уровня PПред загрязнения орбиты КМ.

Составим и решим относительно TС.Ж следующее уравнение:

shestopal11.wmf (5)

Откуда

shestopal12.wmf (6)

Окончание использования околоземной орбиты зависит от возможностей КА противостоять воздействию потока высокоэнергетических частиц КМ. При этом надо учитывать, что средства защиты с течением времени развиваются и включают в себя как пассивные, так и активные средства, построенные на различных физических принципах. Стоимость таких средств защиты различна, эффективность также отличается. Все это приводит к тому, что трудно назначить устраивающее всех конкретное значение PПред. В связи с чем выходом из положения может быть отказ от назначения одного фиксированного значения предельного уровня PПред загрязнения орбиты в пользу рассмотрения некоторого спектра значений PПред, отличающихся по степени предпочтительности [4–6].

Математически такую идею можно реализовать, если использовать возможности теории нечетких множеств (ТНМ). С позиций ТНМ предельный уровень PПред загрязнения орбиты космическим мусором зададим в виде нечеткого множества (НМ):

shestopal13.wmf (7)

функция принадлежности (ФП) shestopal14.wmf которого имеет треугольную форму и описывается в аналитическом виде как

shestopal15.wmf (8)

С учетом (8) НМ shestopal16.wmf можно кратко записать в виде тройки параметров:

shestopal17.wmf (9)

Рассчитаем предельный срок «экологической жизни» околоземной орбиты при нечетком задании предельного уровня загрязнения орбиты космическим мусором. Для этого подставим (7) в формулу (6):

shestopal18.wmf (10)

В выражении (10) необходимо произвести ряд операций над нечеткой величиной shestopal19.wmf. Л. Заде в работе [2] предложен способ нахождения произвольных функций от нечетких аргументов, который называется «принцип обобщения Л. Заде». В нашем случае, когда должна быть вычислена функция вида (10) от одного нечеткого аргумента shestopal20.wmf, применяя принцип обобщения Л. Заде, получаем правило: результат нахождения функции вида (10) от нечеткого аргумента shestopal21.wmf есть нечеткая величина shestopal22.wmf, значения которой рассчитываются по формуле (6), а функция принадлежности совпадает с функцией принадлежности нечеткого аргумента shestopal23.wmf (8).

Математически это правило можно записать следующим образом:

shestopal24.wmf

shestopal25.wmf

shestopal26.wmf (11)

В случае с треугольной функцией принадлежности вида (8) по аналогии с (9) ФП НМ shestopal27.wmf можно записать в виде

shestopal28.wmf (12)

shestopal29.wmf

shestopal30.wmf

shestopal31.wmf

Наиболее наглядно изложенный выше способ прогнозирования предельного срока «экологической жизни» околоземной орбиты при нечетком задании предельного уровня загрязнения орбиты КМ показан на рис. 1.

Способ нечеткого прогнозирования предельного срока «экологической жизни» околоземной орбиты более адекватно учитывает вид фактической неопределенности информации о значении предельного уровня загрязнения орбиты космическим мусором. Однако нечеткость результата прогнозирования не всегда удобна для практического использования. В связи с этим исследуем вопрос возможности преобразования нечеткой оценки предельного срока «экологической жизни» околоземной орбиты к четкому виду.

В теории нечетких множеств существует операция нахождения альфа-сечения НМ, позволяющая преобразовать нечеткое множество к четкому.

Альфа-сечением Aα нечеткого множества shestopal32.wmf называется обычное, то есть четкое подмножество, определяемое формулой shestopal33.wmf α ∈ [0; 1]. Принцип нахождения альфа-сечения показан на рис. 2.

В зависимости от выбора конкретного значения α могут быть найдены различные множества Tα. Очевидно, что необходимо выбирать значение α не произвольно, а некоторым оптимальным способом.

pic_89.wmf

Рис. 1. Способ нечеткого прогнозирования предельного срока «экологической жизни» околоземной орбиты

pic_90.wmf

Рис. 2. Нахождение альфа-сечения Tα нечеткого множества предельного срока «экологической жизни» околоземной орбиты

Будем руководствоваться следующим правилом выбора. Чтобы наилучшим способом преобразовать НМ shestopal34.wmf к четкому множеству, необходимо выбрать такое значение α*, при котором четкое множество shestopal35.wmf наиболее близко к рассматриваемому НМ. В теории нечетких множеств доказано, что для того, чтобы наилучшим способом преобразовать НМ shestopal36.wmf к четкому множеству, необходимо выбрать значение α* = 0,5. Полученное описанным способом четкое множество T0,5 есть интервальная оценка момента окончания срока «экологической жизни» TС.Ж околоземной орбиты по критерию ее загрязнения космическим мусором.

Заключение

В работе получены аналитические выражения для прогнозирования срока «экологической жизни» орбит космических аппаратов с учетом прогрессирующей по закону Кесслера засоренности околоземного космического пространства техногенным мусором.

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

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

Рецензенты:

Козлов В.В., д.т.н., профессор кафедры, Военно-космическая академия имени А.Ф. Можайского, г. Санкт-Петербург;

Садин Д.В., д.т.н., профессор кафедры, Военно-космическая академия имени А.Ф. Можайского, г. Санкт-Петербург.