제6장: 궤도 예측 및 접근 평가

궤도 예측은 현재 궤도 정보로부터 미래 위치와 속도를 계산하는 과정이며, 충돌 회피의 핵심입니다. 정확한 궤도 전파(orbit propagation) 없이는 물체 간 접근을 예측할 수 없고, 충돌 위험을 평가할 수 없습니다. 접근 평가(conjunction assessment)는 두 물체가 언제, 어디서, 얼마나 가까이 지나가는지 결정하고, 충돌 확률을 계산하여 회피 기동의 필요성을 판단합니다.

WIA-SPACE-026 표준은 궤도 전파 방법, 접근 선별(conjunction screening) 프로세스, 충돌 확률 계산 알고리즘, 그리고 회피 기동 결정 기준을 정의합니다. NASA CARA(Conjunction Assessment Risk Analysis), ESA COSY(Collision Avoidance System), 18th SDS(Space Defense Squadron)와 같은 주요 기관들은 매일 수천 개의 위성에 대해 접근 분석을 수행하며, 위험한 경우 운영자에게 경고를 전달합니다. 본 장에서는 이러한 과정의 기술적 세부사항과 실제 운영 경험을 다룹니다.

궤도 전파 방법: SGP4 vs 수치 적분

궤도 전파는 크게 두 가지 접근법으로 나뉩니다: 해석적(analytical) 방법과 수치적(numerical) 방법입니다. 해석적 방법은 섭동 효과를 수학적으로 모델링하여 폐쇄 형태(closed-form)의 해를 제공하며, SGP4가 대표적입니다. 수치적 방법은 뉴턴의 운동 방정식을 미분 방정식 형태로 직접 적분하며, 더 정확하지만 계산 비용이 높습니다.

SGP4 (Simplified General Perturbations 4): SGP4는 TLE와 함께 사용되도록 설계된 전파 모델로, 1960년대에 개발되어 지금까지 사용되고 있습니다. SGP4는 지구 중력장의 비구형성(J2, J3, J4 항), 대기 항력, 그리고 태양 복사압의 단순화된 모델을 포함합니다. 주요 장점은 계산 속도입니다: 단일 epoch에서 임의 시각으로의 전파가 밀리초 이내에 수행됩니다. 이는 25,000개 이상의 물체 카탈로그를 처리하는 데 필수적입니다.

그러나 SGP4의 정확도는 제한적입니다. LEO 위성의 경우 하루 예측 오차가 1-2km, 일주일 예측 오차가 10-50km에 달할 수 있습니다. 이는 대기 항력 모델링의 단순화 때문입니다. SGP4는 BSTAR 항력 계수를 사용하여 항력을 근사하지만, 실제 대기 밀도는 태양 활동, 지자기 폭풍, 계절 변화 등에 따라 크게 변동합니다. 또한 SGP4는 달과 태양의 섭동, 태양 복사압, 지구 반사광(albedo) 등을 무시하거나 단순화합니다. 따라서 정밀 충돌 회피에는 더 정확한 방법이 필요합니다.

수치 적분 (Numerical Integration): 수치 적분은 궤도 운동 방정식을 직접 풉니다: d²r/dt² = F_total/m, 여기서 F_total은 모든 섭동력의 합입니다. 고차 중력장(70×70 또는 그 이상), 정밀 대기 모델(NRLMSISE-00, JB2008), 달/태양 섭동, 태양 복사압, 지구 반사광, 조석력 등을 포함할 수 있습니다. Runge-Kutta, Adams-Bashforth-Moulton, Dormand-Prince 등의 수치 적분 알고리즘이 사용되며, 적응형 스텝 크기로 효율성과 정확도를 균형잡습니다.

수치 적분의 정확도는 SGP4보다 훨씬 높습니다: LEO에서 단기 예측(1-3일)의 경우 수십 미터, 장기 예측의 경우 수백 미터 정확도를 달성할 수 있습니다. 그러나 계산 비용이 크며, 단일 물체를 7일 전파하는데 SGP4의 수천 배 시간이 소요될 수 있습니다. 따라서 실무에서는 하이브리드 접근법을 사용합니다: SGP4로 초기 접근 선별을 수행하고, 위험한 접근만 수치 적분으로 재분석합니다.

표 1: 궤도 전파 방법 비교

방법	계산 속도	정확도 (LEO)	모델링 충실도	사용 사례	입력 형식
SGP4	매우 빠름 (~1ms)	1-5km (하루)	낮음	카탈로그 규모 선별	TLE
SDP4 (심우주)	빠름 (~5ms)	변동적	중간	HEO/GEO 전파	TLE
수치 적분 (중간)	중간 (~1s)	10-100m (하루)	높음	정밀 접근 분석	상태 벡터
수치 적분 (고정밀)	느림 (~10s)	1-10m (하루)	매우 높음	회피 기동 계획	상태 벡터
반해석적 (DSST)	중간 (~100ms)	100m-1km (하루)	중간-높음	장기 예측	상태 벡터/요소

접근 선별 (Conjunction Screening) 프로세스

접근 선별은 카탈로그의 모든 물체 쌍을 검사하여 가까운 접근을 식별하는 과정입니다. 27,000개 물체 카탈로그의 경우 약 3억 6천만 개의 쌍이 존재하며, 매일 이들을 모두 검사해야 합니다. 따라서 효율적인 알고리즘이 필수적입니다. 접근 선별은 여러 단계의 필터링을 사용하여 계산 부담을 줄입니다.

1단계: 궤도 영역 필터링: 궤도 영역(LEO, MEO, GEO)이 크게 다른 물체 쌍은 접근할 가능성이 없으므로 제외합니다. 예를 들어, 400km LEO 물체와 20,000km MEO 물체는 결코 가까이 지나갈 수 없습니다. 반장축, 이심률, 경사각의 간단한 기준으로 불가능한 쌍의 대부분을 제거할 수 있습니다.

2단계: 시공간 게이팅: 각 물체의 궤도를 전파하여 특정 시간 창(예: 향후 7일) 동안의 위치 범위를 계산합니다. 두 물체의 시공간 볼륨이 겹치지 않으면 접근이 불가능합니다. 이는 효율적인 공간 데이터 구조(octree, k-d tree)를 사용하여 가속화됩니다.

3단계: 상세 접근 분석: 게이팅을 통과한 쌍에 대해 상세한 분석을 수행합니다. 두 궤도를 고해상도로 전파하고, 최소 접근 거리(miss distance)와 접근 시각(TCA, Time of Closest Approach)을 계산합니다. 최소 거리가 임계값(예: 5km) 이하이면 "접근 사건(conjunction event)"으로 분류합니다.

4단계: 충돌 확률 계산: 접근 사건에 대해 충돌 확률을 계산합니다. 이는 두 물체의 궤도 불확실성(공분산)을 고려하며, 다음 섹션에서 자세히 다룹니다. 충돌 확률이 임계값(예: 1/10,000) 이상이면 운영자에게 경고를 전달합니다.

18th SDS는 매일 약 20,000개의 접근 사건을 식별하며, 이 중 수백 개가 "관심 사건(events of interest)"으로 분류되어 운영자에게 통보됩니다. 대부분의 접근은 안전한 거리를 유지하지만, 일부는 회피 기동을 필요로 합니다. ISS는 평균적으로 연간 1-3회 회피 기동을 수행하며, 상업 위성 constellation은 더 빈번합니다(예: Starlink는 수천 회/년).

충돌 확률 계산 방법론

충돌 확률(Probability of Collision, Pc)은 두 물체가 충돌할 가능성을 정량화합니다. 이는 궤도 불확실성을 고려한 통계적 계산이며, 단순한 최소 거리보다 훨씬 많은 정보를 제공합니다. 두 물체가 1km 거리로 지나가더라도 불확실성이 크면 충돌 확률이 높을 수 있고, 100m 거리라도 불확실성이 작으면 확률이 낮을 수 있습니다.

기본 개념: 충돌 확률은 두 물체의 상대 위치 확률 분포가 "충돌 단면(combined hard-body region)"과 겹치는 정도를 나타냅니다. 각 물체의 위치는 공분산 행렬로 표현되는 불확실성을 가지며, 일반적으로 다변량 정규 분포로 가정됩니다. 상대 위치 분포는 두 공분산의 합으로 계산됩니다(독립 가정 하에).

2D 근사법 (Chan, Foster, Alfriend): 가장 널리 사용되는 방법은 3D 문제를 2D로 단순화하는 것입니다. TCA에서 두 물체의 상대 속도 방향을 따라 매우 빠르게 지나가므로, 이 방향의 불확실성은 충돌 확률에 거의 기여하지 않습니다. 따라서 상대 속도에 수직인 평면에서의 2D 문제로 근사합니다. 이 평면에서 상대 위치 분포는 2D 정규 분포이고, 충돌 단면은 원(두 물체 반지름의 합)입니다.

2D 충돌 확률은 다음과 같이 계산됩니다: Pc = ∫∫_A p(x,y) dx dy, 여기서 A는 충돌 단면 영역, p(x,y)는 2D 정규 분포입니다. 이 적분은 일반적으로 수치적으로 계산되지만, 특수한 경우(원형 공분산) 해석적 해가 있습니다. Chan과 Foster는 정확한 수치 적분 방법을 개발했고, Alfriend는 빠른 근사 공식을 제안했습니다. 이들은 각각 장단점이 있으며, 실무에서는 정확도와 속도를 고려하여 선택합니다.

Monte Carlo 시뮬레이션: 대체 방법으로 Monte Carlo 시뮬레이션이 있습니다. 공분산 행렬에서 수천 또는 수만 개의 샘플 궤도를 생성하고, 각각을 전파하여 충돌 여부를 확인합니다. 충돌한 샘플의 비율이 충돌 확률의 추정치입니다. Monte Carlo는 개념적으로 간단하고 복잡한 상황(비정규 분포, 기동 중인 물체)에도 적용 가능하지만, 계산 비용이 높습니다. 낮은 확률(예: 10^-5) 추정에는 수백만 샘플이 필요할 수 있습니다.

표 2: 충돌 확률 임계값 및 대응 조치

Pc 범위	위험 수준	대응 조치	통보 대상	추가 분석
< 1×10⁻⁶	무시 가능	없음	없음	없음
1×10⁻⁶ ~ 1×10⁻⁵	낮음	모니터링	운영팀 정보 공유	궤도 업데이트 시 재계산
1×10⁻⁵ ~ 1×10⁻⁴	중간	집중 모니터링	운영자 공식 통보	고정밀 전파, 추가 관측 요청
1×10⁻⁴ ~ 1×10⁻³	높음	기동 검토	운영자 경고 + 권고	회피 기동 계획 수립
> 1×10⁻³	위험	기동 권장/필수	긴급 경고	즉각 기동 실행 (ISS 기준)

NASA CARA: ISS와 NASA 위성 보호

NASA의 CARA(Conjunction Assessment Risk Analysis) 팀은 ISS와 NASA의 모든 운영 위성에 대한 충돌 위험 평가를 담당합니다. CARA는 18th SDS로부터 TLE와 접근 데이터를 받아, 자체적으로 고정밀 분석을 수행합니다. ISS의 경우 매일 수백 개의 잠재적 접근이 평가되며, Pc > 1/10,000 (10^-4) 임계값을 사용합니다.

ISS 충돌 회피 프로토콜: ISS는 유인 우주선이므로 특별히 보수적인 기준을 적용합니다. Pc > 1/10,000이고 접근 거리가 1km 이내인 경우 "Red threshold"로 분류되어 회피 기동이 고려됩니다. 기동 결정은 접근 시각, 기동 비용(연료), 크루 일정, 기타 운영 제약을 종합적으로 고려합니다. 일반적으로 TCA 24-48시간 전에 최종 결정이 내려지며, Progress, ATV, Cygnus 등의 화물선 또는 ISS 자체 엔진을 사용하여 기동합니다.

ISS 기동은 일반적으로 수백 미터의 고도 상승 또는 하강으로 수행되며, 이는 수 m/s의 delta-V를 필요로 합니다. 기동 후에는 궤도를 재평가하고, 다른 접근 사건에 미치는 영향을 분석합니다. 때로는 기동이 다른 접근을 악화시킬 수 있으므로, 여러 사건을 동시에 고려하는 최적화가 필요합니다. CARA는 JSpOC(현 18th SDS)와 긴밀히 협력하여 최신 추적 데이터를 공유하고 분석을 조율합니다.

NASA 과학 위성: Hubble, Terra, Aqua, Aura 등 NASA의 과학 위성도 CARA의 보호를 받습니다. 이들은 ISS만큼 빈번한 기동은 하지 않지만(연료 제약), 임무에 중요한 경우 회피 기동을 수행합니다. 예를 들어, Hubble은 2020년대 초반까지 수차례 회피 기동을 수행했습니다. CARA는 각 위성의 운영팀과 협력하여 기동 가능 여부, 비용, 임무 영향을 평가합니다.

ESA COSY와 18th SDS의 역할

ESA(유럽우주국)의 COSY(Collision Avoidance System)는 유럽 위성들의 충돌 위험을 평가하고 회피 기동을 지원합니다. COSY는 18th SDS, LeoLabs, 그리고 유럽 SST(Space Surveillance and Tracking) 네트워크로부터 데이터를 받아 독립적인 분석을 수행합니다. Sentinel-1, Sentinel-2, Sentinel-3 등 Copernicus constellation과 Swarm, Aeolus 등 과학 위성이 주요 보호 대상입니다.

COSY는 자동화된 접근 선별과 위험 평가를 수행하며, 위험 수준에 따라 운영팀에 통보합니다. Pc > 1×10^-4 (ESA의 "red" 임계값)인 경우 상세 분석과 기동 검토가 시작됩니다. ESA는 또한 예방적 기동(preventive maneuvers)을 수행하기도 하는데, 이는 확률이 임계값 이하이지만 불확실성이 크거나 다른 위험 요소가 있을 때 안전 마진을 확보하기 위한 것입니다.

최근 ESA는 자율 충돌 회피 기술을 개발하고 있습니다. 이는 위성이 스스로 접근을 감지하고 기동을 결정하는 시스템으로, 지상 통제 없이도 작동합니다. 이는 특히 constellation 운영에 중요하며, 수십 또는 수백 개의 위성을 인간이 모두 관리하기 어렵기 때문입니다. 자율 시스템은 기계 학습을 사용하여 접근 패턴을 학습하고, 안전한 기동 전략을 최적화합니다. 2023년 ESA는 실험적 자율 회피를 성공적으로 시연했습니다.

미국 우주군의 18th SDS(구 JSpOC, Joint Space Operations Center)는 전 세계 우주 활동을 위한 SSA 서비스를 제공합니다. 18th SDS는 US SSN을 운영하며, 모든 추적 물체에 대한 카탈로그를 유지하고, 접근 평가를 수행하며, 위험한 접근에 대해 위성 운영자에게 경고를 전달합니다. 이 서비스는 미국 위성뿐 아니라 국제 파트너와 상업 운영자에게도 제공됩니다. 18th SDS는 매일 모든 활성 위성에 대해 접근 선별을 수행하며, 위험한 접근이 식별되면 운영자에게 CDM(Conjunction Data Message)을 전송합니다.

회피 기동 전략과 불확실성 관리

회피 기동은 충돌 확률을 낮추기 위해 위성의 궤도를 변경하는 것입니다. 기동은 연료를 소비하므로 위성 수명을 단축시키며, 임무 운영을 방해할 수 있습니다. 따라서 기동 결정은 위험과 비용의 트레이드오프를 신중히 고려해야 합니다. 최적 기동은 최소 연료로 충돌 확률을 충분히 낮추는 것입니다.

기동 타이밍: 일반적으로 TCA 24-72시간 전에 기동합니다. 너무 일찍 기동하면 예측 불확실성이 크고, 너무 늦으면 준비 시간이 부족합니다. 최적 타이밍은 추적 데이터 업데이트, 기동 계획 및 검증, 지상국 패스 가용성 등을 고려하여 결정됩니다. 긴급 상황(Pc가 빠르게 증가)에서는 12시간 전 또는 그 이내에도 기동할 수 있습니다.

기동 방향과 크기: 기동 방향은 접근 기하학에 따라 결정됩니다. 일반적인 전략은 상대 속도에 수직인 방향으로 miss distance를 증가시키는 것입니다. 접선 방향(along-track) 또는 법선 방향(cross-track) 기동이 주로 사용되며, 둘의 조합도 가능합니다. 기동 크기는 수십 cm/s에서 수 m/s이며, 원하는 Pc 감소량에 따라 계산됩니다. Monte Carlo 시뮬레이션을 사용하여 기동 후 Pc를 예측하고 최적 delta-V를 찾습니다.

불확실성 관리: 충돌 확률 계산의 정확성은 궤도 불확실성 추정에 크게 의존합니다. 불확실성이 과소평가되면 Pc가 과소평가되어 위험한 접근을 놓칠 수 있고, 과대평가되면 불필요한 기동으로 연료를 낭비할 수 있습니다. 따라서 신뢰할 수 있는 공분산 추정과 검증이 중요합니다. 공분산 현실성 검증, 민감도 분석, 앙상블 방법을 통해 불확실성을 관리하며, 이는 신뢰할 수 있는 충돌 확률 계산과 강건한 의사 결정에 필수적입니다.

표 3: 주요 충돌 회피 사례 (2020-2024)

날짜	위성	위협 물체	초기 Pc	기동 후 Pc	Delta-V
2020-09-22	ISS	우주 쓰레기	1.4×10⁻⁴	<10⁻⁶	0.5 m/s
2021-11-11	Sentinel-1A	Cosmos 파편	8×10⁻⁵	2×10⁻⁶	0.3 m/s
2022-07-02	Starlink-1234	OneWeb-0089	5×10⁻⁴	조율	-
2023-03-15	ISS	러시아 ASAT 파편	2×10⁻⁴	<10⁻⁶	0.7 m/s
2024-01-08	Terra	Fengyun 파편	1.5×10⁻⁴	3×10⁻⁶	0.4 m/s