클럭 에러(위성 시계, 수신기 시계)
개요
GNSS(Global Navigation Satellite System)에서 시간 동기화(Time Synchronization)는 매우 중요한 요소이다. 모든 위성은 동일한 기준시계를 유지하려고 하지만, 실제로는 완벽하게 일치하는 시계가 존재하기 어렵다. 또한 지상 수신기 역시 이상적인 기준시계를 가지고 있지 않다. 이로 인해 발생하는 것이 위성 시계 오차와 수신기 시계 오차이다. 이 장에서는 GNSS에서 고려되는 두 종류의 시계 오차가 어떻게 발생하고 모델링되는지 자세히 다룬다.
GNSS 수신에서 위성 시계 및 수신기 시계의 오차는 측정신호(코드 측정치, 위상 측정치)에 매우 직접적으로 영향을 미치므로, 이를 제대로 이해하고 보정하는 과정이 필요하다. 시계 오차는 기본적으로 거리 측정에 선형적으로 포함되며, 가장 기초적인 오차 요소에 속한다. 실제로 GNSS 측위 방정식을 풀기 위해서는 위성 시계와 수신기 시계의 상대적 오차를 모델링하고 보정해야 한다.
위성 시계 오차 이론
기본 개념
위성 시계 오차(Satellite Clock Error)는 GNSS 위성이 기준시계와 정확히 동기화되지 못하여 생기는 오차로서, 다음과 같은 원인들로 인해 발생한다.
원자시계 불안정성: 위성에 탑재된 원자시계(Cesium, Rubidium 등)가 완벽하게 안정적이지 못해 장기(Drift) 및 단기(Flicker) 에이징이 발생
발진기의 온도 특성
전원 공급 변화
위성에 대한 상대론적 효과 (지구 중심 기준좌표계와 위성 궤도 좌표계 간의 상대론적 보정이 필요)
GNSS 위성은 일반적으로 지상국(Ground Control Station)에서 지속적으로 시각 보정을 받으며, 각종 궤도 보정값과 함께 위성 시계 보정 항을 내보낸다. 이는 항시 정확한 UTC(협정 세계시)에 맞추려고 노력하지만, 실제로는 완벽 일치를 이루지 못한다. 그로 인해 위성 시계 오차는 측정치에서 제거해야 하는 중요한 항이 된다.
위성 시계 오차 모델링
위성 시계 오차를 이론적으로 서술하면, 한 위성의 실제 송신 시각과 참조 시각 간의 불일치로 나타난다. 이를 수학적으로 표현하면, 위성 $i$에 대해 다음과 같이 나타낼 수 있다.
여기서
$\Delta t_i$ : 위성 $i$에서 발생하는 시계 오차(초 단위)
$\Delta t_{i,0}$ : 기준 시각 $t_0$에서의 위성 시계 오차
$\dot{\Delta t}_{i}$ : 기준 시각에서 측정된 1차 시계 드리프트(Drift)
$\ddot{\Delta t}_{i}$ : 2차 드리프트(시계 가속도) 항
위 식은 테일러 전개 방식으로 위성 시계의 변화를 나타낸 것이다. GNSS에서는 특정 시점($t_0$)에 측정된 시계 오차와 그 드리프트 값을 예측 모형에 적용함으로써, 임의 시각 $t$에서의 시계 오차를 추정한다. 이는 대개 항을 2차항 이상까지 고려하기보다는, 1차항 정도까지를 주로 고려한다. 시각적 오차가 곧 거리 오차로 환산될 때, 광속 $c$를 곱한 형태로 위성 시계 오차가 측정식에 포함된다.
위성 시계 오차는 기본적으로 위성 시간 – 기준 시간 간의 차이를 의미하며, 실제 관측 방정식에는 $c \cdot \Delta t_i$로 들어가게 된다. 여기서 $c$는 진공 중 광속이다.
위성 시계 보정 신호
GNSS(예: GPS)는 내장된 Navigation Message를 통해 위성 시계 보정 파라미터를 제공한다. 다음은 GPS 신호에 포함되는 시계 보정 파라미터 예시이다.
Clock Bias 계수($a_{f0}$)
Clock Drift 계수($a_{f1}$)
Clock Drift Rate 계수($a_{f2}$)
이들은 단순화된 다항식을 통해 위성 시계 오차를 추정하기 위한 값이며, 수신기는 이를 이용해 다음과 같이 보정값을 계산한다.
위 공식에서 $(t - t_0)$는 메시지에 명시된 기준 에포크와 수신 시각 간의 차이다. 이 때의 결과 $\Delta t_{sat}$을 관측 방정식에서 빼주어 위성 시계 오차를 보정한다. 물론 여기에 상대론적 보정(Relativistic Correction) 항이 추가로 들어간다. GPS나 다른 GNSS에서도 궤도 이심률에 따른 상대론적 보정 항을 다음과 같이 준다.
여기서
$\mathbf{a}$ : 위성 궤도 장반경
$\mu$ : 지구 중력상수
$e$ : 궤도 이심률
$E$ : 궤도 타원에서의 이심각
$c$ : 광속
이와 같은 시계 보정 매커니즘을 통해 위성 시계 오차를 최소화한다.
수신기 시계 오차 이론
수신기 시계 오차의 정의
수신기도 이상적인 원자시계를 사용하지 않는 한, 내부 발진기의 오류 및 온도 변화, 안정도 등에 의해 측정 시각이 기준 시간과 불일치한다. 이를 **수신기 시계 오차(Receiver Clock Error)**라고 한다. 수신기 시계 오차는 위성으로부터 도달하는 신호를 수신할 때 신호의 도달 시간을 기록하는 과정에 직접 영향을 주므로, 측정된 거리(= 신호 도달 시간 × 광속)에 오차를 유발한다. 마찬가지로 위성 시계 오차와 동일하게, $c \cdot \Delta t_r$의 형태로 관측식에 들어가게 된다. 이때 $\Delta t_r$은 수신기의 시계 오차이다.
수신기 시계 오차의 발생 원인
발진기(OCXO, TCXO 등)의 온도 민감성
전원 불안정성
노화(Aging)에 따른 장기 드리프트
잡음성 진동(Phase Noise)
GNSS 수신기는 측정치(코드, 위상 등)에서 수신기 시계 보정 항을 별도로 추정해야 한다. 일반적으로 GNSS 측정에서 ‘시계 오차’라고 하면 위성 시계 오차와 함께 수신기 시계 오차도 동시에 고려된다.
수신기 시계 오차 모델링
수신기 시계 오차(Receiver Clock Error)는 수신기의 내부 발진기 특성을 나타내는 파라미터로, 일반적으로 GNSS 측정치에서 “추정해야 할 불확정 파라미터”로 포함된다. 수신기 측정치(코드 측정, 위상 측정)를 수학적으로 표현할 때, 수신기 시계 오차는 시간 바이어스와 그 드리프트(변화율)로 모델링될 수 있다. 예를 들어 수신기 시계 오차를 다음과 같은 다항식 형태로 나타낼 수 있다.
여기서
$\Delta t_r$: 수신기 시계 바이어스
$\Delta t_{r,0}$: 기준 시각 $t_0$에서의 수신기 시계 바이어스
$\dot{\Delta t}_{r}$: 기준 시각에서의 시계 드리프트(1차 항)
$\ddot{\Delta t}_{r}$: 시계 가속도(2차 항)
실제로는 수신기가 에포크마다 시계 오차를 재추정하기 때문에, GNSS 수신 소프트웨어에서 장기 드리프트나 가속도 항을 모두 반영하기보다는 매 에포크마다 독립적인 바이어스를 추정하는 방식을 많이 사용한다. 그러나 고정밀 응용에서의 시계 안정성 분석에는 위 다항식을 간단히 1차 혹은 2차까지 확장하여 사용할 수 있다.
수신기 시계 보정 기법
수신기 시계 오차를 측정식에서 제거(또는 추정)하기 위해서는 다음과 같은 방법들이 활용된다.
동시 추정(Estimation in the Navigation Filter) GNSS 측정 방정식을 세울 때, 수신기의 위치(3개 좌표)와 함께 수신기 시계 바이어스($\Delta t_r$)를 하나의 파라미터로 둔다. 예를 들어, 코드 측정식(유클리드 거리 + 위성 시계 보정 포함)을 단순화하여 표현하면 아래와 같다.
Pi=∥r−ri∥+cΔtr−cΔti+ϵi여기서
$P_i$: 위성 $i$에 대한 코드 유효측정치(피seudorange)
$\mathbf{r}$: 수신기 위치 벡터
$\mathbf{r}_i$: 위성 $i$의 위치 벡터
$c$: 광속
$\Delta t_r$: 수신기 시계 바이어스
$\Delta t_i$: 위성 $i$의 시계 바이어스(항상 보정 파라미터를 통해 처리됨)
$\epsilon_i$: 기타 오차(잡음, 전리층/대류층 오차 등)
이 때 $\Delta t_r$는 네 번째 미지수로 포함되어 에포크마다 추정된다. 이 방법을 통해 모든 위성으로부터의 측정치가 한꺼번에 얽혀 수신기 시계 바이어스를 최소제곱해 추정하게 된다.
차분(Differencing) 기법 다중 위성 신호를 동시에 수신할 때, 서로 다른 위성으로부터 온 측정치 간에 차분을 취하면 공통 수신기 시계 오차가 소거된다. 예를 들어, 단일차분(Single Difference)은 두 위성 $i, j$의 측정값을 서로 빼서 수신기 시계 항을 제거하거나, 여러 수신기 간의 비교(기준국 - 이동국)에서 시계 오차 일부를 상쇄시킬 수 있다. 다만, 차분 기법은 여러 위성/여러 기지국과의 측정값을 조합해야 하므로, 오차 전파 경로가 달라질 수 있고, 잡음이 증가할 수 있다는 단점이 있다. 하지만 고정밀 측위(예: RTK)에서는 수신기 시계 오차뿐 아니라 위성 시계 오차도 차분을 통해 상당 부분 제거가 가능하다는 이점이 있다.
외부 타이밍(External Timing) 참조 어떤 수신기는 시각 동기를 별도로 맞춰주는 안정화된 외부 기준(예: 루비듐 원자시계, 지상 시각표준국 신호 등)을 참조하여 시계 바이어스를 최소화한다. 이러한 장치는 전력, 유지보수 비용이 많이 드는 대신, 수신기 자체의 시계 안정도가 높아져 추가 보정 항을 단순화할 수 있다.
수신기 시계 오차 보정 기법은 수신기의 사용 목적(일반 내비게이션, 정밀 측량, 시간동기 등)과 가용 자원(고성능 발진기, 지상기준국 등)에 따라 선택된다.
측정 방정식에서의 시계 오차 처리
GNSS 기본 측정(코드, 반송파 위상)을 관측 방정식으로 세웠을 때, 시계 항은 아래와 같은 형태로 나타난다.
여기서
$P_i$: 코드 측정치(pseudorange)
$\Phi_i$: 반송파 위상 측정치(carrier phase)
$\mathbf{r}$, $\mathbf{r}_i$: 수신기/위성 위치 벡터
$\Delta t_r, \Delta t_i$: 수신기/위성 시계 오차
$I_i, I'_i$: 전리층 영향 (코드/위상에서의 영향이 다름)
$T_i$: 대류층 영향
$\epsilon_{P,i}, \epsilon_{\Phi,i}$: 잔여 잡음 및 기타 미지 오차
수신기는 위성 시계 보정 정보(내비게이션 메시지 등)로 $\Delta t_i$를 보정한 후, 남은 부분을 수신기 시계 항 $\Delta t_r$로서 추정한다. 이때 위성 시계 보정 정보는 보완적 오차(상대론 효과 포함)를 포함하므로, 모델 정확도에 따라 남는 오차가 달라진다.
시계 드리프트와 동적 모델
GNSS 필터링 기법(예: 칼만 필터, 확장 칼만 필터)을 적용할 때, 수신기 시계 오차를 동적 상태 변수로 포함하기도 한다. 예를 들어 수신기 시계 바이어스 $\Delta t_r$와 시계 드리프트 $\dot{\Delta t}_{r}$를 상태 변수로 두어, 시계 드리프트가 작은 일부 확률과정(예: 랜덤 워크, 아주 작은 분산의 백색 잡음)으로 가정한다. 이 경우, 시간(에포크) 간에
와 같은 이산형 상태 방정식으로 시계 항들을 예측하고, 측정 업데이트 과정에서 실제 GNSS 코드/위상 데이터를 이용하여 보정한다.
$\Delta T$: 필터의 에포크 간격
$w_1, w_2$: 잡음 항(필터 설계에 따라 달라짐)
이와 같이 동적 모델을 도입하면, 각 에포크마다 시계 오차를 단순히 독립 파라미터로 두는 것보다 시간적 상관성을 반영할 수 있어, 추정 성능을 개선할 수 있다.
상대론적 효과와 시계 오차
GNSS에서 시계 오차를 정확히 다루기 위해서는 고전역학적 시간 개념만으로는 부족하며, 상대론적 효과(Relativistic Effect)를 고려해야 한다. 위성은 지구 주위를 고속(약 3~4 km/s)으로 공전하고 있으며, 지구 표면의 중력 퍼텐셜과도 다소 차이가 있다. 특수 상대론 및 일반 상대론에 의하면, 빠르게 움직이는 물체나 다른 중력 퍼텐셜 영역에 있는 물체의 시계는 지상 기준시계와 다르게 흐른다. 이를 상대론적 시간 지연 또는 **상대론적 주파수 편이(Shift)**라고 한다.
특수 상대론 효과
위성이 고속 이동할 때 발생하는 시간 지연
위성에 탑재된 원자시계는 지상 기준시계에 비해 다소 느리게 흐른다.
위성 속도가 $v$, 광속을 $c$라고 할 때, 시계 비율은 대략 $\sqrt{1 - \frac{v^2}{c^2}} \approx 1 - \frac{v^2}{2c^2}$로 표현된다.
일반 상대론 효과
지구 중력 퍼텐셜 차이에서 오는 시간 지연
지구 중심부나 지표면보다 중력 퍼텐셜이 높은 곳(위성 고도)에서는 시계가 더 빨리 간다.
위성과 지상국 간 높이 차이에 의한 퍼텐셜 차이($\Delta U$)를 고려하여, 시계 비율이 $1 + \frac{\Delta U}{c^2}$ 꼴로 나타난다.
GPS 사례를 보면, 특수 상대론 효과(위성 운동)로 인해 위성 시계는 하루에 약 $-7,\mathrm{\mu s}$ 정도 느려지고, 일반 상대론 효과(중력 퍼텐셜)로 인해 약 $+45,\mathrm{\mu s}$ 빨라진다. 최종적으로는 $+38,\mathrm{\mu s}$ 수준이 되어, 시스템 설계 시 이를 미리 보정하거나 운영 중 보정 파라미터에 반영한다. 이러한 보정이 적용되지 않으면, 하루에 수십 미터 이상의 측위 오차가 누적될 수 있다.
GNSS 위성 시계 오차 보정 공식에는 이 상대론적 효과가 간단한 형태로 추가된다. 예를 들어, GPS의 경우 이심 궤도에서 발생하는 상대론적 보정 항(이심각 $E$ 기반)은 다음과 같이 나타난다.
이 외에도 위성의 공전 궤도, 지상국 모니터링, 운용 스케줄에 따라 세부적으로 조정된 값이 내비게이션 메시지에 반영된다.
GNSS 시계 동기화 절차
GNSS는 자체 시간 스케일을 유지하고, 이를 보정하기 위해 지상 추적국과 위성 간에 주기적인 동기화를 수행한다. 이 과정에서 위성 시계를 모니터링하여, 편차가 발생하면 그 보정값(Orbit/Clock correction)을 업링크(Uplink)하여 위성에 반영한다. 이어서 위성은 그 값을 내비게이션 메시지로 다시 방송한다. 여기서 주요한 점은 **“모든 GNSS 위성 시계를 최대한 동일한 기준 타임 스케일로 맞추는 것”**이다. 예를 들어,
GPS Time은 1980년 1월 6일을 주 시작 에포크(Week Epoch)로 하는 자체 시간 스케일
Galileo System Time(GST), GLONASS Time, BeiDou Time(BDT) 등도 각각의 시스템 고유 스케일
일반적으로는 UTC(협정 세계시)와의 오프셋(Leap Second 고려 등)이 수신기에서 함께 보정됨
따라서 GNSS 수신기는 내비게이션 메시지에 담긴 시계 파라미터를 사용해 먼저 GNSS Time을 복원하고, 필요시 UTC나 다른 시간 스케일로 변환할 수 있도록 한다.
GPS Time과 수신기 시계
GPS Time은 윤초(Leap Second)를 고려하지 않고 계속 증가하는 주 단위 카운트(Week)와 주 내초(SOW)로 표현된다. 반면, 실제 국제 표준시(UTC)에는 일정 주기로 윤초가 삽입된다. 수신기 내부에서는 “GPS 시간”을 기준으로 동작하다가, 시각 표시나 다른 시스템 연동(예: Network Time Protocol, 타이밍 어플리케이션) 시점에 UTC 변환 과정을 거친다. 이 때도 수신기 시계($\Delta t_r$)가 정확히 GPS Time과 동기되어 있어야 올바른 변환이 가능하다.
시계 오차의 추가 이슈: 지상국 동기와 배치
지상국(Tracking/Monitoring Station)의 시계 또한 원자시계에 기반하여 각 위성의 궤도 및 시계 정보를 모니터링한다. 만약 지상국 시계가 부정확하거나 관리가 소홀하면, 위성 시계 보정 파라미터에 영향을 준다. 따라서 글로벌 GNSS 운용체계는 다음과 같은 요소들을 관리한다.
전 지구적으로 배치된 모니터링 스테이션: 동시에 여러 위성을 추적, 수집된 측정치로 Orbit/Clock 상태를 추정
지상 마스터 컨트롤 스테이션: 전체 위성 궤도/시계 파라미터를 산출하고, 특정 시간 간격으로 업링크
원자시계(루비듐, 세슘 등)의 교정: 지상국과 위성 시계 간의 상호 검증 및 차이 추적
만약 한 모니터링 스테이션이 오프셋을 일으키면, 그 오프셋이 추정된 모든 위성 시계 값에 편향(Bias)을 줄 수 있다. 다만, GNSS는 여러 관측소 데이터를 종합하여 견고하게 추정하므로, 단일 지상국 문제만으로 전체 시스템이 크게 흔들리지는 않는다.
시각 동기 및 측정 원리 예시(mermaid)
GNSS 위성 시계 및 수신기 시계가 어떻게 동기되고, 보정 정보가 어떻게 순환되는지 간략히 나타내면 아래와 같은 흐름이 있다.
위성은 지상국에서 업데이트된 시계 보정 파라미터와 궤도 정보를 실시간 송신
수신기는 위성 신호를 받아 각 위성까지의 코드/위상 측정치와 함께 시계 파라미터를 이용해 수신기 시계를 추정
수신기 시계가 동기화되면, 정확한 측위와 표준 시간 산출이 가능해짐
지상국은 지속적으로 여러 위성을 모니터링하여 시계 및 궤도 오차가 누적되지 않도록 보정값을 생성, 위성으로 전송
고정밀 응용에서의 시계 오차 처리
일반 내비게이션 수준에서는 위성 시계 보정 정보(방송 궤도/시계)만으로 수 m 정도의 정확도를 얻을 수 있다. 그러나 고정밀 측위(예: 정밀 측량, 과학기술 응용 등)에서는 수 cm 이하의 정확도가 요구되므로, 훨씬 정교한 시계 오차 처리 기법이 필요하다.
Precise Ephemeris 및 정밀 시계 보정 국제 GNSS 서비스(IGS, International GNSS Service) 등에서는 방송 궤도/시계보다 훨씬 정밀한 후처리 제품(Precise Orbit and Clock)을 제공한다. 이들은 전 지구에 걸쳐 분산된 IGS 스테이션이 여러 GNSS 위성을 관측한 결과를 토대로, 매우 높은 정확도(수 cm 수준의 궤도, 수 ns 수준의 시계)를 확보하여 공개한다. 이를 사용하는 후처리(post-processing) 방식으로 측정하면, 위성 시계 오차를 크게 줄일 수 있다.
실시간 고정밀 서비스(RTK, PPP 등)
RTK(Real-Time Kinematic): 기준국과 이동국이 같은 위성을 관측하고, 위성 시계 오차 및 기지국 시계 오차가 동일하게 반영되는 부분을 차분(Differencing)으로 제거해 준다. 이때 기준국 자체도 고정밀 좌표와 안정된 시계를 가정하거나, 더 상위의 네트워크(RTK Network)에서 보정받는다.
PPP(Precise Point Positioning): Precise Orbit/Clock을 실시간(또는 근실시간)으로 제공받아, 단일 수신기에서도 고정밀 측위를 가능케 한다. 상용 위성 기반 보정 서비스나 인터넷 기반 보정 서비스를 통해 위성 시계 오차를 줄이고, 수신기 시계도 필터링 과정에서 정교하게 추정한다.
다중 GNSS 시스템 활용 GPS 외에도 GLONASS, Galileo, BeiDou 등 여러 GNSS를 동시에 관측하면, 훨씬 많은 위성/주파수를 이용하여 시계 오차 추정이 안정적이고 정교해진다. 다만 시스템별 기준 시계가 다를 수 있고(예: GPS Time vs. GLONASS Time), 위성 시계 오차 모델이나 메시지 포맷이 달라서 통합 처리가 복잡해진다. 일반적으로는 다음과 같은 **시스템 간 시간 오프셋(ISB, Inter-System Bias)**까지 고려해야 한다.
Δtr(GPS)−Δtr(GLONASS),Δtr(GPS)−Δtr(Galileo),…이와 같은 추가 편차까지 추정해야 다중 GNSS를 통합하여 위성/수신기 시계 오차를 제대로 보정할 수 있다.
시계 품질 평가 및 안정도 분석
시계의 품질(시계 안정도)은 일반적으로 **ADEV(Allan Deviation)**나 TDEV(Time Deviation) 등 통계 지표로 분석한다. 위성에 탑재된 원자시계, 지상국 원자시계, 수신기의 발진기(OCXO/TCXO/루비듐 등)가 모두 다른 안정도 특성을 가진다.
장기 안정도: 드리프트(Aging)나 온도 변화에 의해 느리게 변하는 오차
단기 안정도: 위상 잡음(Phase Noise), 순간적인 열적/전자적 변동에 의한 오차
GNSS 시계 처리에서는 이러한 시계 잡음을 모델링하여 필터 설계(예: 랜덤 워크, 루프 필터) 시 반영한다. 이를 통해 예측 단계에서 과도하게 드리프트를 추정하거나, 반대로 지나치게 강한 제약을 두어 오차가 누적되지 않도록 조정한다.
정밀 타이밍 응용에서의 고려사항
GNSS는 위치 측정뿐 아니라, 안정된 시각 동기가 필요한 많은 분야(통신기지국 동기, 전력망 동기 등)에도 활용된다. 이러한 정밀 타이밍 응용에서는 특히 다음을 신중히 고려한다.
운영 중 시간 오프셋 모니터링: $1,\mathrm{ns}$ 편차는 곧 $30,\mathrm{cm}$ 정도의 위치 오차에 해당하므로, 시간/주파수 표준으로는 작은 편차도 큰 영향이 될 수 있다.
루프 홀드오버(Hold-over) 기능: GNSS 신호가 일시적으로 끊겨도 내부 발진기로 일정 시간 정확도를 유지해야 한다. 이를 위해 고품질 발진기(OCXO, Rubidium 등)을 사용해 시계를 안정적으로 보강한다.
여러 GNSS 시스템/주파수 동시 사용: 장애나 스푸핑 공격 등에 대비해 이중·삼중 안전장치로서 여러 시스템을 병행 감시하며 시계 이상을 조기 탐지한다.
위성 시계와 수신기 시계의 예측 정밀도
위성 시계도 지상국 모니터링에서 업링크하기까지 일정 시간이 소요되어, “최신 추정값과 실제 시계 상태” 간의 편차가 발생한다. 통상 GPS 등에서는 15분2시간 주기로 시계 파라미터를 갱신하거나, 더 빈번히 갱신하는 확장 메시지를 사용하기도 한다. 수신기 시계는 필터링으로 순간순간 추정이 가능하지만, 실제로는 잡음이 존재하므로 짧은 시간(초 단위) 스케일에서는 비교적 정확하게 맞출 수 있어도, 장기(분시간 단위)에서는 발진기 특성에 따라 누적 편차가 커질 수 있다. 이러한 특성을 이해하고 모델링해야, 고성능 GNSS 측정이 가능하다.
Last updated