위성 궤도(Orbital Error)

궤도 오차의 개념

GNSS 위성은 특정 궤도를 따라 지구를 공전하며, 수신기 측에서는 위성 궤도 정보를 포함한 천문력(에페머리스, Ephemeris)을 활용하여 위성의 위치를 추정한다. 위성의 실제 궤도와 예측된 궤도가 정확히 일치하지 않을 경우, 측정된 위성-수신기 거리(측정 신호의 전파 도달 시간에 기반)는 실제 거리와 편차가 발생할 수 있다. 이를 통틀어 “위성 궤도 오차”라고 한다.

GNSS에서 위성 궤도 오차는 일반적으로 수 미터 수준으로 나타날 수 있으며, 이 오차가 누적되거나 보정되지 않으면 측위 결과에 심각한 영향을 끼친다. 따라서 정확한 궤도 모델링과 정밀천문력(Precise Ephemeris) 데이터가 중요한 역할을 한다.

궤도 기술 방식과 천문력

GNSS에서 위성의 궤도는 기본적으로 케플러 궤도 이론에 근거하여 표현되며, 미세 오차와 섭동(perturbation)을 반영하기 위해 오비탈 요소를 시간에 따라 보정해 준다. 일반적으로 위성 궤도를 나타낼 때 케플러 요소를 직접 제공하기보다는, 천문력 메시지에 포함된 일련의 계수를 통해 수신기에서 자체적으로 궤도를 예측하도록 한다.

에페머리스 기반 궤도 계산

수신기는 방송천문력(Broadcast Ephemeris)에서 제공하는 다음과 같은 변수들을 이용하여 위성 위치를 계산한다.

• 궤도 점근선과 관련된 참조 시각

• 궤도 모양을 나타내는 궤도 이심률 $e$

• 반장축 길이(또는 궤도 주기를 유도할 수 있는 제곱근 준장반경)

• 궤도 경사각, 상승 교점 적경, 근점 편각 등 케플러 요소에 대응하는 보정값

• 제곱 항 및 선형 항 보정을 포함하는 각종 천문력 보정 계수

이러한 데이터를 통해 수신기는 특정 시각에서 위성의 케플러 궤도를 예측하고, 여기에 지구의 비구형 중력장, 태양풍, 달-태양 중력 영향, 지자기장 상호작용 등 여러 섭동 요소를 고려할 수 있도록 일부 추가 계수를 적용하여 위성의 위치를 근사적으로 산출한다.

케플러 궤도에서의 기본 식

케플러 궤도에서 위성은 다음과 같은 이차체(two-body) 운동 방정식을 만족한다. 지구 중심을 기준으로 할 때, 위성의 상태 벡터 $\mathbf{r}$과 가속도 $\ddot{\mathbf{r}}$ 사이의 관계는

$$\ddot{\mathbf{r}} = -\mu \frac{\mathbf{r}}{\|\mathbf{r}\|^3}$$

여기서 $\mu$는 지구의 표준 중력 변수(지구 중력 상수 $GM$)이며, $|\mathbf{r}|$은 벡터 $\mathbf{r}$의 유클리드 노름이다.

케플러 요소

위성 궤도는 일반적으로 다음 여섯 개의 케플러 요소(Keplerian Elements)로 정의할 수 있다.

1. 궤도 장반경 $a$ (또는 준장반경 $\sqrt{a}$)

2. 이심률 $e$

3. 궤도 경사각 $i$

4. 승교점 적경 $\Omega$ (Right Ascension of the Ascending Node)

5. 근점 편각 $\omega$

6. 평균 근점 편각에 대응하는 위상 요소(진앙, 평균근점통과시각, 편각 등)

실제로 GNSS 방송천문력에서는 이들 요소에 대해 시간에 따른 보정항(선형/이차 보정 등)과 추가 섭동 항목이 함께 제공된다.

위성 궤도 오차 발생 원인

GNSS 위성 궤도에 오차가 발생하는 요인은 크게 다음과 같이 구분할 수 있다.

1. 중력장 모델 오차 지구의 중력장은 완전한 구형이 아니며, 지구 내부 질량 분포의 불균질, 적도 팽창 등의 영향으로 인해 위성 공전에 추가적인 섭동이 발생한다. 이 효과는 일반적으로 지구 중력장 고차 항(자전, J2, J3 항 등)을 통해 모델링되지만, 모델 불일치로 오차가 남을 수 있다.

2. 태양-달 등 3체 섭동 태양 및 달의 중력 또한 무시할 수 없으며, 특히 고궤도 위성의 경우 태양-달의 위치에 따라 상당한 섭동이 발생한다. GNSS 위성은 MEO(Medium Earth Orbit)에 위치하여 하루 기준 수차~수십 cm 정도의 궤도 섭동이 누적될 수 있다.

3. 태양 복사압(Solar Radiation Pressure) 위성에 도달하는 태양 광자(Photon)에 의한 복사압이 위성에 미세한 힘을 가하여 궤도 변화를 일으킨다. 위성의 태양 전지판 각도나 위성 외형에 따라 복사압 모델이 달라지며, 완벽히 모델링하기 어렵기 때문에 궤도 오차의 원인이 될 수 있다.

4. 지자기장 상호작용 위성에 탑재된 여러 장비 및 추진계에는 지자기장과의 상호 작용에 따른 미세한 힘과 토크가 작용하며, 이 역시 완전히 보정하기 어렵다.

5. 추진계 오차 및 기타 잔류 요인 위성 자세 제어나 궤도 유지(Station Keeping) 과정에서 추진계의 동작 오차, 연료의 질량 변화, 위성 자세 제어 모드에 따른 불연속적 요인 등이 발생하여, 이상적 궤도에서 미세하게 벗어나는 원인이 된다.

오차 전파(Propagation)와 영향

위성 궤도 오차는 단순한 위치 오차로 그치지 않고, 측정 신호의 전파 경로와 위상에도 영향을 끼칠 수 있다. 일반적으로 GNSS의 측정 신호는 위성의 시간 기준에 따라 발사되므로, 위성 시각과 궤도 정보가 동기화되지 못할 경우 거리 측정에 추가 오차가 발생한다.

위성 위치 오차가 거리 측정에 미치는 영향

위성 위치 오차 $\Delta \mathbf{r}$가 있을 때, 측정된 거리 오차 $\Delta \rho$는 다음과 같이 근사할 수 있다.

$$\Delta \rho \approx \frac{\mathbf{r}_{\text{sat}} - \mathbf{r}_{\text{rcv}}}{\|\mathbf{r}_{\text{sat}} - \mathbf{r}_{\text{rcv}}\|} \cdot \Delta \mathbf{r}$$

여기서 $\mathbf{r}{\text{sat}}$은 실제 위성 위치 벡터, $\mathbf{r}{\text{rcv}}$는 수신기 위치 벡터이고, $\Delta \mathbf{r}$는 실제와 예측 위성 위치 차이(오차)이다. 결국 위성 위치 오차가 측정 거리로 직결되어, GNSS 최종 위치 해석에 오차를 유발한다.

간단한 궤도 결정(Dynamic Orbit Determination) 모델

궤도 결정(Orbit Determination)은 위성의 위치와 속도를 추정하기 위해, 동역학 방정식과 관측 데이터를 결합하여 해를 구하는 과정을 말한다. 가장 기본적인 2체 운동 방정식에 보정 항을 추가한 형태로 모델링할 수 있다.

$$\ddot{\mathbf{r}} = -\mu \frac{\mathbf{r}}{\|\mathbf{r}\|^3} + \mathbf{a}_{\text{pert}}$$

$\mathbf{a}_{\text{pert}}$는 앞서 언급한 고차 중력장, 태양 복사압, 3체 섭동, 지자기장 효과, 추진계 교정 등 여러 물리량을 합산한 섭동 가속도이다. 측정 방정식은

$$\rho_{\text{obs}} = \|\mathbf{r}_{\text{sat}} - \mathbf{r}_{\text{rcv}}\| + \epsilon_{\rho}$$

로 표현할 수 있는데, 여기서 $\epsilon_{\rho}$는 각종 측정 잡음 및 오차 항을 포함한다. 이 두 식(동역학 방정식과 측정 방정식)을 동시에 풀어 궤도와 오차 추정값을 결정한다. 실제로는 EKF(Extended Kalman Filter), SRIF(Square Root Information Filter), Batch Least Squares 등 다양한 추정 기법이 활용된다.

궤도 보정

GNSS 측면에서 궤도 보정은 크게 두 가지 형태로 이뤄진다.

• 방송천문력(Broadcast Ephemeris) 보정 위성 운용국에서 정기적으로 추정한 궤도 정보(에페머리스)를 업링크하여, 수신기에 최신 정보를 제공한다.

• 정밀천문력(Precise Ephemeris) 보정 IGS(International GNSS Service) 등을 통해 전 세계 관측소 데이터를 기반으로 정밀한 궤도를 추정하고, 이를 별도로 제공한다. PPP(Precise Point Positioning)나 기타 정밀 측위에서는 이 정밀천문력 데이터를 활용하여 궤도 오차를 획기적으로 줄인다.

위와 같은 여러 요인이 복합적으로 작용하여 발생하는 위성 궤도 오차는 다른 GNSS 오차 요소(시계 오차, 전리층 오차 등)와 함께 최종 측위를 결정한다. 이후 실제 측정 환경에서 이러한 각종 오차를 어떻게 제거(또는 줄일) 것인지는 보정 기법 파트에서 다룬다.

지구 중력장 고차 항 모델

GNSS 위성 궤도를 결정하고 유지하기 위해서는 지구 중력장을 정확히 모델링해야 한다. 하지만 지구는 완전한 구형이 아니며, 질량 분포가 균일하지도 않다. 이를 반영하기 위해 지구 중력장은 일반적으로 지구 중심 좌표계(예: ECEF 좌표계)에서 조화함수(Spherical Harmonics)를 이용하여 확장한다. 기본적인 두 객체(2-body) 운동 방정식에 지구 중력장 고차 항(대표적으로 $J_2, J_3, \dots$)을 포함시키면 다음과 같은 형태를 띤다.

$$\ddot{\mathbf{r}} = -\mu \frac{\mathbf{r}}{\|\mathbf{r}\|^3} \;+\; \mathbf{a}_{J_2} \;+\; \mathbf{a}_{J_3} \;+\;\dots \;+\; \mathbf{a}_{\text{기타 섭동}}$$

여기서 $\mathbf{a}{J_2}$, $\mathbf{a}{J_3}$ 등은 각 차수별 지구 중력장 계수를 이용해 계산되는 섭동 가속도 항이다.

$J_2$ 항(지구의 적도 팽창) 예시

$J_2$는 지구 적도 부분이 팽창되어 있는 효과(Oblateness)를 나타내는 계수로, 가장 큰 고차 항 섭동 요인이다. 그 공식적 형태는 다음과 같은 추가 가속도로 표현된다.

$$\mathbf{a}_{J_2} = -\frac{3}{2} J_2 \frac{\mu R_{\mathrm{E}}^2}{\|\mathbf{r}\|^4} \begin{bmatrix} x \left(1 - 5 \frac{z^2}{\|\mathbf{r}\|^2}\right) \\ y \left(1 - 5 \frac{z^2}{\|\mathbf{r}\|^2}\right) \\ z \left(3 - 5 \frac{z^2}{\|\mathbf{r}\|^2}\right) \end{bmatrix}$$

• $R_{\mathrm{E}}$: 지구 평균 반지름

• $x, y, z$: 지구 중심 좌표계에서의 위성 위치 좌표

• $|\mathbf{r}| = \sqrt{x^2 + y^2 + z^2}$

• $J_2$: 지구 중력장 계수 중 2차 항

이 항은 지구의 적도 돌출로 인해 위성 궤도에 추가 가속도가 발생함을 의미한다. 이를 포함하지 않으면 고도 수백~수천 km 단위에서 점차적 궤도 편차가 누적될 수 있다.

태양 복사압(Solar Radiation Pressure) 모델

태양 복사압은 GNSS 위성에 작용하는 주요 섭동 요인 중 하나로, 태양에서 방출된 광자가 위성 표면과 상호 작용할 때 생기는 힘을 말한다. 위성 본체와 태양 전지판에 반사·흡수되는 양이 정확히 모델링되지 않으면, 궤도 예측 오차가 발생한다.

단순 플럭스 모델

태양 복사압을 단순화하여 위성 단면적 $A$, 반사율(흡수율) $\alpha$, 태양 방사 플럭스 $S$ 등을 이용하면 다음과 같은 형태로 표현할 수 있다.

$$\mathbf{a}_{\text{SRP}} = \frac{S \, A \, (1 + \alpha)}{m c} \hat{\mathbf{n}}$$

• $m$: 위성 질량

• $c$: 광속

• $\hat{\mathbf{n}}$: 태양 방향으로의 단위 벡터

• $S$: 태양 복사 플럭스 (약 $1367,\text{W/m}^2$ @ 1 AU, 지구 근처)

실제로는 위성의 기하학적 형상, 태양 패널의 회전 각도, 위성 표면 특성(흡수/확산/정반사 비율)에 따라 훨씬 복잡한 모델이 필요하다.

박스-윙(Box-Wing) 모델

실제 위성 외형을 “박스와 윙(태양 전지판)”으로 단순화하여, 각 면에 대해 태양 복사압을 산정하고 합산하는 방식이 널리 사용된다. 이때는 각 면의 질량 중심 오프셋, 포인터(ADCS) 상태 변화, 태양각에 따른 동적 반사율 변화 등을 고려해야 한다. GNSS 위성 운영 기관(예: GPS, GLONASS, Galileo 등)은 위성의 공칭 모델 및 다양한 측정 자료를 토대로 복사압 계수를 추정하여 정기적으로 갱신하고, 정밀 궤도 결정에 반영한다.

태양/달 3체 섭동

지구-위성 간의 2체 문제만 고려해도 1차적인 궤도 해석은 가능하지만, 실제로는 태양이나 달의 중력장이 무시하기 힘들 만큼 영향을 준다. GNSS 위성 고도(약 20,200 km 전후)에서 태양 및 달의 중력에 의한 섭동은 일시적으로 수십 cm 이상의 오차로 작용할 수 있다.

3체 문제 접근

태양-지구-위성(또는 달-지구-위성) 3체 문제를 완전히 해석하는 것은 이론적으로 복잡하다. 대신, 천체력(예: 천문력에 기반한 달, 태양의 위치 정보를 이용)으로부터 태양/달의 위치 벡터 $\mathbf{r}{\odot}, \mathbf{r}{\mathrm{moon}}$를 얻고, 다음과 같이 추가 가속도를 근사적으로 계산한다:

$$\mathbf{a}_{\text{sun}} = \mu_{\odot} \left( \frac{\mathbf{r}_{\odot} - \mathbf{r}}{\|\mathbf{r}_{\odot} - \mathbf{r}\|^3} - \frac{\mathbf{r}_{\odot}}{\|\mathbf{r}_{\odot}\|^3} \right) \\ \mathbf{a}_{\text{moon}} = \mu_{\mathrm{moon}} \left( \frac{\mathbf{r}_{\mathrm{moon}} - \mathbf{r}}{\|\mathbf{r}_{\mathrm{moon}} - \mathbf{r}\|^3} - \frac{\mathbf{r}_{\mathrm{moon}}}{\|\mathbf{r}_{\mathrm{moon}}\|^3} \right)$$

• $\mathbf{r}$: 지구 중심에서 본 위성 위치 벡터

• $\mathbf{r}{\odot}$, $\mathbf{r}{\mathrm{moon}}$: 지구 중심에서 본 태양, 달 위치 벡터

• $\mu_{\odot}$, $\mu_{\mathrm{moon}}$: 태양, 달에 대한 중력상수

위 식에서 보이는 것처럼, “3체 섭동”은 위성-태양(또는 달) 간 중력과, 지구-태양(또는 달) 간 중력 차이를 통해 계산된다. 이 차이가 위성에 추가 가속도를 부여한다.

자세(Attitude) 및 추진계 오차

GNSS 위성은 추력기(Thruster)를 통해 궤도를 유지하고, 위성 몸체와 태양 패널을 태양 방향 등에 적절히 조향하기 위해 자세 제어(ADCS: Attitude Determination and Control System)를 수행한다. 이 과정에서 발생하는 불완전한 추력, 자세 편차 등은 미세하지만, 장기적으로 누적되어 궤도 예측에 오차를 일으킨다.

• 추진계 불균형: 추력 벡터가 위성 질량 중심에서 어긋나거나, 추력기의 성능 편차가 존재할 수 있음.

• 연료 소모에 따른 질량 변화: 궤도 유지용 연료가 소모되면서 위성 질량이 작아지고, 복사압/추력 등이 달라짐.

• 비선형 자세 제어 모드: 태양 전지판이 태양을 추적하는 과정에서 비연속적인 회전이 발생할 때, 순간적으로 복사압 모델에 예기치 못한 변화를 유발할 수 있음.

지자기장 상호작용

고궤도 위성이라 할지라도 지자기장의 영향이 완전히 사라지지는 않는다. 특히 위성 내부 전자장비와의 상호 작용(가령, 잔류 자기 모멘트, 전류 루프 등)은 미약한 자기력, 로런츠 힘 등을 일으킬 수 있으며, 그 결과 오비탈 섭동이 나타난다. 이를 정밀 모델링하기 위해서는 위성 자체의 자기 특성과 지자기장 모형(IGRF 등)을 함께 고려해야 한다.

위에서 살펴본 다양한 섭동 요소들은 서로 상호작용하여 종합적인 궤도 편차를 유발한다. GNSS 시스템에서는 이를 최소화하기 위해 정밀 궤도 결정 과정을 수행하며, 그 결과가 수신기에게 전달되어야만 한다.

위성 궤도 오차 구분: Radial, Along-track, Cross-track

GNSS 위성 궤도 오차는 일반적으로 위성의 공전 궤도 면을 기준으로 다음 세 방향 성분으로 나누어 해석한다.

• Radial(반경방향) 성분: 위성 궤도 중심(지구 중심)에 대해 위성의 거리 오차가 발생하는 방향

• Along-track(궤도 진행방향) 성분: 위성의 공전 진행 방향에 대한 오차

• Cross-track(궤도 수직방향) 성분: 궤도 면에 수직하게 발생하는 오차

이렇게 구분하면 궤도 오차가 어느 방향으로 크게 발생하는지 쉽게 파악할 수 있고, 보정 알고리즘을 설계할 때 해당 성분별 계수를 통해 편리하게 적용할 수 있다.

Radial 오차

Radial 성분은 위성-지구 중심 간의 거리 예측값과 실제 거리 간 불일치로 발생한다. 예를 들어 중력장 모델이 부정확하거나 태양 복사압 모델이 부정확하면, 궤도 고도가 과대/과소로 예측되어 Radial 오차를 유발한다. Radial 오차는 직접적으로 GNSS 측정 신호 거리에 가장 큰 영향을 주기 때문에, 궤도 결정에서 매우 중요한 지표로 취급된다.

Along-track 오차

Along-track 성분은 위성의 공전 운동에 따라 축적되는 위성의 위치 오차를 나타낸다. “예상 시간 대비 위성이 더 앞서 있거나, 뒤처져 있다”라는 형태로 해석할 수 있다. 방송천문력에서 시간을 기점으로 예측하는 위성의 위치에 비해, 실제 위성이 공전 주기나 섭동 요소로 인해 더 빨리(또는 더 늦게) 지나가는 경우가 여기에 해당된다.

Cross-track 오차

Cross-track 성분은 궤도면을 기준으로 하여, 실제 위성이 궤도면 밖으로 벗어나서 발생하는 오차다. 이는 궤도 경사각 $i$ 또는 승교점 적경 $\Omega$ 등의 변동(또는 섭동)이 충분히 정확히 예측되지 못한 경우 발생한다. Cross-track 오차는 레이의 투영 거리에 의한 영향이 Radial이나 Along-track만큼 직접적으로 크지는 않을 수 있으나, 누적되면 측정 범위나 위성가시도 예측 등에 영향을 줄 수 있다.

방송천문력(Broadcast Ephemeris)와 예측 범위

GNSS 위성이 방송하는 천문력(Broadcast Ephemeris)은 미래 일정 시간(수 시간~하루 정도)에 대한 궤도 예측 값을 제공한다. 즉, 매 시각마다 새로 즉시 관측되는 위성 궤도가 아니라, 특정 알고리즘(궤도 역학 + 추정 기법)으로 미리 계산된 궤도 매개변수를 사용자가 받아보는 것이다.

• 업데이트 주기: GPS의 경우 보통 2시간 간격으로 업데이트되는 방송천문력을 수신기에 제공한다. 수신기는 이 방송천문력에 적힌 기준 시각(Ephemeris reference epoch) 이후 특정 시간 범위 안에서 위성 위치를 예측한다.

• 예측 오차: 방송천문력은 궤도 파라미터를 간략화하여 전송하므로, 예측 시간이 길어질수록 섭동 요인을 완벽히 반영하기 어렵게 되어, 오차가 누적된다.

IODE(정보 식별자)와 적용 구간

방송천문력 메시지에는 IODE(Issue of Data, Ephemeris)나 IODC(Issue of Data, Clock) 같은 식별자가 포함되어, 해당 궤도/시계 정보가 유효한 시간 범위를 명시한다. 수신기는 이를 참조하여 가장 최신의 궤도/시계 정보를 사용한다. 만약 IODE가 지난 오래된 천문력을 적용한다면, 오차가 커져 정확도가 급격히 떨어질 수 있다.

정밀천문력(Precise Ephemeris)과 개선

IGS(International GNSS Service)를 비롯한 다양한 기관에서는 전 지구에 분포된 GNSS 관측 자료를 수집하고, 매우 정밀한 궤도 및 시계 정보를 계산해 제공한다. 이를 “정밀천문력(Precise Ephemeris)”라고 부른다. 정밀천문력은 다음과 같은 과정으로 산출된다.

1. 전 세계 여러 지상 관측소에서 GNSS 관측 데이터를 수집한다.

2. 고차 중력장, 태양/달 중력, 복사압, 위성 안테나 페이즈 센터 등 세부적인 보정값을 반영하여 궤도와 시계 궤적을 추정한다.

3. 필터 알고리즘(베치(Batch)나 칼만 필터 계열)을 사용해 여러 시간 구간의 데이터를 동시 추정한다.

4. 수 시간~수 일 후, 최종적으로 정밀궤도가 완성되어 공개된다.

이 정밀천문력을 사용하면, 일반적인 방송천문력에 비해 궤도 오차를 수십 cm~1 m 이하 수준으로 줄일 수 있다. PPP나 정밀 측위 기법에서 이를 활용하여 고정밀 위치 결과를 얻을 수 있다.

다른 GNSS 간 비교

GPS뿐 아니라 GLONASS, Galileo, BeiDou, QZSS 등 다른 GNSS 시스템도 마찬가지로 위성 궤도 예측을 위한 방송천문력을 운영하고, IGS(또는 각 GNSS 서비스 기관)에서 정밀천문력을 제공한다. 다만, 각 시스템별로 사용되는 궤도 구조(예: GLONASS FDMA와 그 궤도 특징), 업데이트 주기, 전송 형식, 예측 정확도 등이 조금씩 다르며, 이는 결과적으로 위성 궤도 오차 규모에도 영향을 미친다.

위성 궤도 오차는 GNSS 오차 요소 중에서도 매우 중요한 축을 차지하며, 위성 시계 오차와 더불어 GNSS 측위 정확도를 좌우하는 핵심 요소이다. 여기까지는 “위성 궤도 오차”라는 주제에 대한 전반적인 설명이며, 이후 다른 보정 기법이나 활용 사례는 별도 섹션에서 다룰 수 있다.