시슬루너 공간을 위한 최소 가정 초기·궤도 결정과 확률적 추적 프레임워크

본 논문은 시슬루너(지구‑달 사이) 공간에서의 목표물(Resident Space Object, RSO) 추적을 위한 새로운 초기 궤도 결정(IOD) 및 궤도 결정(OD) 프레임워크를 제시한다. 전통적인 가우스 방법은 두 몸(Kepler) 역학과 짧은 관측 간격을 전제로 하지만, 시슬루너 영역에서는 달의 중력이 크게 작용해 이러한 가정이 깨진다. 따라서 저자들은 최소 가정(minimal‑assumption) 접근법을 채택한다. Ⅰ. 배경 및 필요성 시슬루너 영역은 4π 입체적 부피가 GEO 이하 궤도보다 약 1000배 넓으며, 최근 인도·중국·미국 등 여러 국가의 달 탐사 임무가 증가함에 따라 공간 상황 인식(SSA)과 충돌 회피가 중요해졌다. 기존의 IOD 방법으로는 가우스, 제한된 허용 영역(CAR), 확률적 허용 영역(PAR) 등이 있으나, 이들 모두 두 몸 역학에 의존하거나 다수의 사전 통계 가정을 필요로 한다. 특히 PAR은 범위 통계와 Jacobi 상수 등을 사전에 설정해야 하며, 복잡한 동역학 모델을 사용한다. Ⅱ. 제안된 IOD 방법 저자들은 시슬루너 목표물이 한 관측소의 시야에 10~20시간 동안 머무를 수 있다는 점을 활용한다. 관측소는 고정된 광학 센서로 방위(AZ)와 고도(EL)를 초당 여러 번 측정하고, 제한적인 거리(ρ) 정보를 추가한다. 이렇게 얻은 시계열 데이터를 다항식 스플라인으로 피팅하고, 시간에 대한 1차·2차 미분을 통해 위치와 속도를 추정한다. 관측 노이즈(각도 1.5 arcsec, 거리 5~10% 수준)를 반영해 여러 번 샘플링하면, 각 샘플은 초기 상태 후보 입자(particle)가 된다. 이 입자 집합을 Gaussian Mixture Model(GMM)로 클러스터링해 확률 밀도 함수(PDF)를 만든다. Ⅲ. 동역학 및 측정 모델 시슬루너 영역의 동역학은 원형 제한 삼체 문제(CR3BP)로 모델링한다. 질량 비율 μ≈0.01215, 비차원화된 거리와 시간 단위를 사용해 방정식(1a‑1c)을 제시한다. 측정 모델은 지구 중심 관측소(O)와 목표물(T) 사이의 상대 위치를 barycentric 좌표계에서 topocentric NEU 좌표계로 변환하고, 방위·고도 각을 계산한다. 거리 정보는 TDOA, 일방향 도플러(OWD) 등으로 추정하지만, 본 연구에서는 거리 노이즈가 크다는 가정 하에 각도 전용 측정만으로도 충분히 초기 상태를 생성할 수 있음을 보인다. Ⅳ. OD 단계 – Particle Gaussian Mixture(PGM) 필터 초기 입자‑GMM PDF를 기반으로 PGM 필터를 적용한다. PGM은 (1) 파티클 전파, (2) 클러스터링(새 GMM 생성), (3) 측정 업데이트, (4) 리샘플링의 네 단계로 구성된다. 파티클 전파 단계에서는 CR3BP 방정식을 수치 적분하고, 클러스터링 단계에서는 EM 알고리즘을 이용해 최적의 가우시안 컴포넌트를 찾는다. 측정 업데이트는 각 컴포넌트별로 Kalman‑like 업데이트를 수행하고, 리샘플링은 입자 고갈을 방지한다. 이 구조는 비선형·혼돈성이 큰 삼체 시스템에서도 가중치와 컴포넌트 수를 동적으로 조절해 차원 저주를 완화한다. Ⅴ. 시뮬레이션 및 성능 평가 다양한 궤도 유형에 대해 실험을 수행한다. 1) L1·L2 근처의 불안정 궤도: 급격한 궤도 변동과 민감한 초기 조건 때문에 EKF/UKF가 발산하지만, PGM은 안정적인 추적을 유지한다. 2) L4·L5 주변의 안정 궤도: 모든 필터가 수렴하지만, PGM이 가장 작은 RMS 오차를 기록한다. 3) 장거리 원궤도(수십만 km): 각도 전용 측정만으로도 초기 입자 구름을 형성하고, 장기간(수일) 동안 오차가 1 km 이하로 유지된다. 비교 대상에는 EKF, UKF, EnKF, 적응형 GMM 필터가 포함되었으며, PGM은 특히 비선형 구간에서 평균 오차를 30~50% 감소시켰다. Ⅵ. 거리 정보 활용 방안 초기 상태 추정에 필요한 최소 거리 정보를 얻는 세 가지 방법을 제시한다. (a) 직접 레이더 거리 측정, (b) 수동 RF 기반 TDOA, (c) 일방향 도플러를 이용한 상대 속도 추정. 각 방법은 측정 노이즈 모델링을 통해 입자 생성 단계에 반영된다. Ⅶ. 한계 및 향후 연구 - 관측이 완전히 각도 전용일 경우 거리 추정이 부정확해질 수 있다. - 입자 수가 많아지면 계산 비용이 급증하므로 GPU 병렬화가 필요하다. - 현재는 단일 목표물 추적에 초점을 맞추었으며, 다중 목표물(MTT) 확장은 추가 연구가 필요하다. 결론적으로, 이 논문은 시슬루너 공간이라는 특수한 환경에서 기존 두 몸 기반 IOD가 갖는 근본적인 한계를 극복하고, 확률적 입자‑GMM 혼합 필터를 통해 장기 궤도 결정을 실현하는 실용적인 방법론을 제공한다.