В работе представлены описание численной модели движения искусственных спутников Луны (ИСЛ), разработанной авторами, и результаты исследования особенностей динамики объектов в окололунном пространстве. Было проведено численное моделирование движения 5180 объектов, равномерно распределенных в окололунном пространстве. Показано, что в ряде областей окололунного пространства время жизни ИСЛ на орбитах весьма мало в зависимости от выбора начального значения большой полуоси и наклонения, и выяснены причины этого явления. Данные результаты могут быть использованы при выборе стратегии размещения спутниковых систем долговременного использования в окололунном орбитальном пространстве.