LEO 위성 단일관측으로 방송형 GNSS 스푸퍼 위치추정

초록

본 논문은 저궤도(LEO) 위성에 탑재된 단일 GNSS 수신기가 한 번의 궤도 통과 동안 방송형 GNSS 스푸퍼의 위치를 추정하는 방법을 제시한다. 스푸퍼가 발생시키는 비정상적인 주파수 변동을 수신기의 시계 드리프트 추정값에서 제거함으로써, 수신기‑스푸퍼 사이의 레인지‑레이트(time‑history)를 얻고 이를 비선형 최소제곱으로 처리해 지리적 위치를 복원한다. 실험에서는 Spire Global과 협력해 알려진 지상 스푸퍼를 비GNSS 대역에서 송신하고, LEO 수신기가 이를 포착해 10 m 이내의 오차로 위치를 재현하였다.

상세 분석

이 연구는 LEO 기반 GNSS 수신기가 제공하는 연속적인 궤도·속도 정보를 활용해, 방송형 스푸퍼가 방출하는 다중 신호의 공통 도플러 성분을 추출하는 새로운 알고리즘을 제안한다. 기존의 다중 위성·다중 수신기 기반 T/FDOA 방식은 동시 캡처와 시간 동기화가 필요해 운영 비용이 크게 증가하지만, 본 방법은 단일 위성·단일 패스만으로도 충분히 정확한 위치 추정이 가능하도록 설계되었다. 핵심 아이디어는 스푸퍼가 각 위성 신호에 삽입하는 ‘가짜’ 도플러 ˜fₙ(t) 를 무시하고, 모든 신호에 공통으로 나타나는 시계 드리프트 γ(t) 를 GNSS 수신기의 내부 PVT 추정기가 자동으로 추정한다는 점이다. γ(t) 에는 수신기‑스푸퍼 사이의 레인지‑레이트 ˆrᵀvᵣ, 수신기 자체 시계 드리프트 δ̇tᵣ, 송신기(스푸퍼) 시계 드리프트 δ̇tₜ, 그리고 스푸퍼가 가짜 위치에 대해 가정한 시계 드리프트 δ̇t̃ᵣ이 포함된다. 이 중 ˆrᵀvᵣ가 가장 큰 기여를 하므로, γ(t) 를 시간에 따라 측정하면 스푸퍼와 수신기 사이의 거리 변화율을 직접 얻을 수 있다.

제안된 비선형 최소제곱 추정기는 γ(t) 를 입력으로 받아 스푸퍼의 3차원 좌표를 최적화한다. 수식적으로는 γ(t)=−(1/λ)·ˆrᵀvᵣ+δ̇tᵣ−δ̇tₜ+δ̇t̃ᵣ 로 표현되며, 여기서 λ는 캐리어 파장이다. 알고리즘은 먼저 γ(t) 를 고정된 시간 간격으로 샘플링하고, LEO 위성의 알려진 궤도·속도 정보를 이용해 ˆrᵀvᵣ 를 계산한다. 이후 δ̇tᵣ, δ̇tₜ, δ̇t̃ᵣ 은 각각 위성 자체 시계 모델, 스푸퍼의 알려진(또는 추정 가능한) 시계 특성, 그리고 스푸퍼가 가짜 위치에 대해 설정한 시계 드리프트로 보정한다. 보정 후 남은 항목이 순수한 레인지‑레이트가 되며, 이를 기반으로 최적화 문제를 풀어 스푸퍼의 위치를 도출한다.

오차 분석에서는 송신기 시계 불안정성(위상 잡음, 주파수 드리프트)과 스푸퍼가 의도적으로 삽입할 수 있는 최악의 주파수 변조를 고려하였다. 수학적으로는 시계 불안정성을 σ²·Δt 형태의 잡음 항으로 모델링하고, 최악의 경우 γ(t) 에 0.5 Hz 수준의 추가 편향이 발생할 수 있음을 보였다. 시뮬레이션 결과, 이러한 편향이 존재하더라도 1 km 이내의 위치 오차로 제한되며, 실제 실험에서는 10 m 수준의 정밀도를 달성하였다.

실험은 Spire Global의 LEO 위성에 탑재된 GNSS 수신기를 이용해, 비GNSS 대역(예: 2.4 GHz)에서 알려진 좌표의 지상 스푸퍼를 송신한 뒤 데이터를 수집하였다. 수신기는 스푸퍼 신호를 기존 GNSS 신호와 구분해 추적하고, 내부 PVT 추정기로부터 γ(t) 를 추출하였다. 추정된 γ(t) 를 이용해 비선형 최소제곱을 수행한 결과, 실제 스푸퍼 위치와 평균 8 m, 최대 12 m 차이의 오차를 보였다. 이는 기존 다중 위성 기반 방법보다 비용·복잡도는 낮지만, 실용적인 정확도를 제공함을 입증한다.

본 연구는 LEO 기반 스펙트럼 감시 서비스가 방송형 스푸퍼 탐지·위치추정에 바로 적용될 수 있음을 시사한다. 특히, 저비용 단일 위성 수집만으로도 전 세계적인 실시간 감시가 가능해, 항공·해운 분야의 GNSS 보안 강화에 큰 기여를 할 것으로 기대된다. 향후 연구에서는 이동형(비정지) 스푸퍼, 다중 주파수 동시 캡처, 그리고 시계 불안정성에 대한 실시간 보정 알고리즘을 확장할 계획이다.