To address the need for robust positioning, navigation, and timing services in lunar environments, this paper proposes a novel fault detection framework for satellite constellations using inter-satellite ranging (ISR). Traditionally, navigation satellites can depend on a robust network of ground-based stations for fault monitoring. However, due to cost constraints, a comprehensive ground segment on the lunar surface is impractical for lunar constellations. Our approach leverages vertex redundantly rigid graphs to detect faults without relying on precise ephemeris. We model satellite constellations as graphs where satellites are vertices and inter-satellite links are edges. We identify faults through the singular values of the geometric-centered Euclidean distance matrix (GCEDM) of 2-vertex redundantly rigid sub-graphs. The proposed method is validated through simulations of constellations around the Moon, demonstrating its effectiveness in various configurations. This research contributes to the reliable operation of satellite constellations for future lunar exploration missions.
- ID статьи: 2406.09759
- Название: Autonomous Constellation Fault Monitoring with Inter-satellite Links: A Rigidity-Based Approach
- Авторы: Кейдай Иияма, Даниэль Неамати, Грейс Гао (Стэнфордский университет)
- Классификация: cs.RO (Робототехника)
- Время публикации/конференция: 2024 ION GNSS+ Conference
- Ссылка на статью: https://arxiv.org/abs/2406.09759
Для обеспечения надежных услуг позиционирования, навигации и синхронизации времени (PNT) в лунной среде в данной работе предлагается новая схема обнаружения неисправностей спутникового созвездия с использованием межспутниковых измерений расстояния (ISR). Традиционно навигационные спутники могут полагаться на надежную сеть наземных станций для мониторинга неисправностей. Однако из-за ограничений по стоимости создание комплексного наземного сегмента для лунного созвездия нецелесообразно. Предложенный метод использует графы с избыточной жесткостью вершин для обнаружения неисправностей без зависимости от точных эфемерид. Спутниковое созвездие моделируется как граф, где спутники являются вершинами, а межспутниковые каналы — ребрами. Неисправности выявляются по сингулярным значениям матрицы геометрически центрированных евклидовых расстояний (GCEDM) 2-вершинно-избыточно-жестких подграфов. Метод проверен на моделировании созвездия вокруг Луны и продемонстрирована его эффективность при различных конфигурациях.
Основная проблема, которую решает данное исследование, заключается в том, как реализовать автономное обнаружение неисправностей спутникового созвездия в лунной среде при отсутствии наземных станций мониторинга. В частности:
- Отсутствие наземных станций мониторинга: невозможно создать сеть станций, подобную земной системе SBAS, на поверхности Луны
- Сложность получения точных эфемерид: из-за ограниченного количества станций и низкой стабильности бортовых часов сложно получить точные эфемериды
- Проблема циклической зависимости: если измерения ISR используются для определения орбиты и синхронизации времени (ODTS) для генерации эфемерид, это создает циклическую зависимость между ODTS и обнаружением неисправностей
Данное исследование имеет важное значение для будущих лунных миссий:
- Поддержка LunaNet: обеспечение надежных услуг PNT для сети LunaNet НАСА
- Критически важные миссии: обеспечение надежной работы лунных навигационных спутников критически важно для безопасности
- Экономическая эффективность: предоставление экономически эффективного решения для мониторинга неисправностей в условиях ограниченных ресурсов на Луне
Традиционные методы мониторинга неисправностей спутников имеют следующие ограничения:
- Зависимость от точных эфемерид: требуют точной информации об орбите для расчета ожидаемых расстояний
- Сильные предположения: предполагают, что измерения ISR достаточно точны, хорошо откалиброваны и не содержат ошибок
- Зависимость от наземных станций: требуют сети наземных станций мониторинга
- Предложена схема онлайн-обнаружения неисправностей на основе теории жесткости: без необходимости в точных эфемеридах или наблюдениях с наземных станций
- Определены условия топологии графа, необходимые для обнаружения неисправностей: доказано, что граф должен быть 2-вершинно-избыточно-жестким для обнаружения неисправного спутника
- Предоставлена математическая теория рангов EDM и GCEDM: доказана математическая основа использования четвертого и пятого сингулярных значений GCEDM для обнаружения неисправностей
- Верификация моделирования лунного созвездия: продемонстрировано влияние гиперпараметров и амплитуды неисправностей на производительность обнаружения
Входные данные: двусторонние измерения расстояния между спутниками rijВыходные данные: идентификация и локализация неисправного спутника
Ограничения: отсутствие точных эфемерид, ограниченная поддержка наземных станций
Спутниковое созвездие моделируется как взвешенный граф G=⟨V,E,W⟩:
- Вершины V: спутники
- Ребра E: межспутниковые каналы
- Веса W: значения расстояний
\|x_i - x_j\| + w_{ij} + f_{ij} = \frac{c}{2}\tau_{ij} & (i \neq j) \\
0 & (i = j)
\end{cases}$$
где:
- $w_{ij} \sim N(0, \sigma_w)$: шум измерения
- $f_{ij} = f_i + f_j$: смещение неисправности
- $f_k = \begin{cases} \bar{f} & \text{спутник k неисправен} \\ 0 & \text{спутник k исправен} \end{cases}$
#### Основная теорема теории жесткости
**Теорема II.4**: Для взвешенного графа G граф является опровергаемым по неисправностям тогда и только тогда, когда G является k-вершинно-избыточно-жестким (k≥2) и тогда и только тогда, когда G является k-реберно-избыточно-жестким (k≥2).
### Матрица геометрически центрированных евклидовых расстояний (GCEDM)
#### Метод построения
1. Построение EDM: $D_{ij} = r_{ij}^2$
2. Геометрическое центрирование: $G = -\frac{1}{2}J^n D J^n$
3. где $J^n = I_n - \frac{1}{n}11^T$ — матрица геометрического центрирования
#### Ключевые свойства
**Предложение III.1**: Ранг EDM удовлетворяет $\text{rank}(D_{n,d,m}) \leq \min(d + 2 + 2m, n)$
**Предложение III.2**: Ранг GCEDM удовлетворяет $\text{rank}(G_{n,d,m}) \leq \min(d + 2m, n-1)$
где m — количество неисправных спутников, d=3 для трехмерного пространства.
### Алгоритм обнаружения неисправностей
#### Тестовая статистика
$$\gamma_{\text{test}} = \frac{\lambda_4 + \lambda_5}{\lambda_1}$$
где $\lambda_i$ — i-е сингулярное значение GCEDM.
#### Процедура алгоритма
1. **Обнаружение клик**: использование алгоритма Чибы-Нишизеки для поиска всех k-клик подграфов (k≥5)
2. **Онлайн-обнаружение**:
- Построение GCEDM для каждой k-клики
- Расчет тестовой статистики $\gamma_{\text{test}}$
- Если $\gamma_{\text{test}} > \bar{\gamma}_{\text{test}}$, подграф отмечается как содержащий неисправность
- Определение неисправного спутника по максимальному элементу четвертого сингулярного вектора $u_4$
3. **Подтверждение неисправности**: подсчет количества раз, когда каждый спутник отмечен как неисправный; спутник считается неисправным, если это количество превышает пороговое значение
#### Ключевые гиперпараметры
- **Размер клики k**: рекомендуется использовать k=6 для повышения точности идентификации неисправного спутника
- **Интервал времени обнаружения DI**: контролирует количество временных шагов, используемых для идентификации неисправности
- **Минимальное количество подграфов с неисправностями $\delta_{nf}$**: избегает определения неисправности на основе небольшого количества образцов
- **Минимальный коэффициент обнаружения неисправности $\delta_{rf}$**: должен удовлетворять $\delta_{rf} < \frac{1}{n_{fs,max}}$
- **Пороговое значение тестовой статистики $\bar{\gamma}_{\text{test}}$**: контролирует компромисс между TPR и FPR
## Экспериментальная установка
### Конфигурация созвездия
Используется эллиптическое лунное замороженное орбитальное (ELFO) созвездие:
- **Количество спутников**: 12 спутников, равномерно распределенные на 4 орбитальных плоскостях
- **Параметры орбиты**:
- Большая полуось: 6142,4 км
- Эксцентриситет: 0,6
- Наклонение: 57,7°
- Долгота восходящего узла: -90°, 0°, 90°, 180°
### Параметры моделирования
- **Шум измерения**: $\sigma_w = 1$ м
- **Количество неисправностей**: 1, 2, 3 неисправных спутника
- **Амплитуда неисправности**: $\bar{f} = 5, 8, 10, 15, 20$ м
- **Интервал времени обнаружения**: DI = 1, 2, 3, 5
- **Установка порогов**: 95%, 99%, 99,9% квантили
### Метрики оценки
1. **Коэффициент истинно положительных результатов (TPR)**: $\frac{TP}{TP + FN}$
2. **Коэффициент ложно положительных результатов (FPR)**: $\frac{FP}{FP + TN}$
3. **Показатель P4**: $\frac{4 \cdot TP \cdot TN}{4 \cdot TP \cdot TN + (TP + TN) \cdot (FP + FN)}$
## Результаты экспериментов
### Основные результаты
1. **Влияние амплитуды неисправности**: с увеличением амплитуды неисправности улучшаются значения TPR, FPR и P4
2. **Верхний предел TPR**: TPR составляет примерно 0,92, поскольку примерно 8,3% спутников находятся вблизи перилуния, что приводит к ограниченному количеству самосодержащихся подграфов из-за затенения
3. **Компромисс порогов**:
- Снижение порога $\bar{\gamma}_{\text{test}}$: увеличение TPR, увеличение FPR
- Увеличение длины обнаружения DI: увеличение TPR, увеличение FPR
### Ключевые выводы
1. **Геометрическое влияние**: обнаружение работает лучше, когда спутники распределены в трехмерном пространстве; обнаружение затруднено, когда неисправные спутники находятся в одной плоскости
2. **Влияние количества неисправностей**: чем меньше количество неисправных спутников, тем лучше производительность обнаружения
3. **Оптимизация параметров**:
- Малая амплитуда неисправности: высокий порог + длительное время обнаружения работают лучше
- Большая амплитуда неисправности: низкий порог + короткое время обнаружения работают лучше
### Показатели производительности
На примере одного неисправного спутника:
- **Амплитуда неисправности 20м, 99% порог, DI=1**: TPR=0,900, FPR=0,006, P4=0,952
- **Амплитуда неисправности 5м, 99,9% порог, DI=5**: TPR=0,124, FPR=0,005, P4=0,345
## Связанные работы
### Традиционные методы мониторинга неисправностей
1. **Системы SBAS**: зависят от сети наземных станций мониторинга
2. **Алгоритм RAIM**: автономный мониторинг целостности приемника
3. **Методы ISR**: работы Вольфа (2000) и Родригеса-Переса и др. (2011)
### Преимущества данной работы
1. **Отсутствие необходимости в точных эфемеридах**: решает основное ограничение лунной среды
2. **Прочная теоретическая основа**: обеспечивает математические гарантии на основе теории жесткости
3. **Адаптация к динамической топологии**: повышает надежность за счет обнаружения в нескольких подграфах
## Заключение и обсуждение
### Основные выводы
1. Предложена первая схема автономного обнаружения неисправностей спутникового созвездия на основе теории жесткости
2. Доказано, что 2-вершинная избыточная жесткость является необходимым и достаточным условием для обнаружения неисправностей
3. Верифицирована эффективность метода сингулярных значений GCEDM в лунном созвездии
### Ограничения
1. **Геометрические ограничения**: требуется распределение спутников в трехмерном пространстве; обнаружение затруднено при плоском распределении
2. **Обнаружение множественных неисправностей**: производительность жадного алгоритма снижается при наличии нескольких неисправностей
3. **Вычислительная сложность**: алгоритм обнаружения клик имеет высокую сложность
### Направления будущих исследований
1. **Распределенная реализация**: разработка распределенных алгоритмов обнаружения неисправностей
2. **Интеграция методов**: объединение методов остатков для повышения надежности
3. **Прогностические модели**: разработка моделей прогнозирования распределения статистики при отсутствии неисправностей
4. **Оптимизация выбора подграфов**: выбор эффективных подграфов для снижения вычислительной нагрузки
## Глубокая оценка
### Преимущества
1. **Высокая теоретическая новизна**: впервые применена теория жесткости к обнаружению неисправностей спутников, предоставлена полная математическая теоретическая база
2. **Высокая практическая ценность**: решает реальные инженерные проблемы в лунной среде, имеет важное значение для будущих глубокосмических исследований
3. **Полнота метода**: от теоретического доказательства к реализации алгоритма и верификации моделирования образует полную цепь исследований
4. **Математическая строгость**: предоставлены строгие доказательства рангов EDM и GCEDM, заполнены пробелы в соответствующей теории
### Недостатки
1. **Ограничения предположений**: предполагается, что неисправность представляет собой постоянное смещение, в то время как реальные неисправности могут быть более сложными
2. **Чувствительность к геометрии**: высокие требования к геометрическому распределению спутников, что может ограничить практическое применение
3. **Ограничения экспериментов**: верификация только в среде моделирования, отсутствие тестирования в реальных условиях
4. **Вычислительная эффективность**: алгоритм обнаружения клик имеет высокую сложность, что может повлиять на реальную производительность
### Влияние
1. **Академический вклад**: вводит новые теоретические инструменты в область обнаружения неисправностей спутников
2. **Инженерная ценность**: предоставляет важные рекомендации для проектирования будущих лунных навигационных систем
3. **Воспроизводимость**: подробное описание алгоритма облегчает воспроизведение и улучшение
### Применимые сценарии
1. **Глубокосмические исследования**: среда Луны, Марса и других мест, где отсутствует наземная инфраструктура
2. **Автономные системы**: спутниковые созвездия, требующие высокой степени автономии
3. **Среда с ограниченными ресурсами**: сценарии с ограниченными вычислительными и коммуникационными ресурсами
## Библиография
В статье цитируется 81 соответствующий источник, охватывающий важные работы в области теории жесткости, теории графов, спутниковой навигации, обнаружения неисправностей и других областей, обеспечивающих прочную теоретическую основу для исследования.
---
**Общая оценка**: Это отличная статья, имеющая важное значение как с точки зрения теоретической инновации, так и практического применения. Авторы умело применили теорию жесткости к области обнаружения неисправностей спутников, решив реальные инженерные проблемы в лунной среде. Хотя в некоторых аспектах существуют ограничения, ее новаторский теоретический вклад и практическая ценность делают ее важным прогрессом в этой области.