업데이트: 2026년 5월
목차
서론
InSAR 지반침하 모니터링은 광역 인프라의 침하 감시에 가장 효율적인 위성기반 원격탐사 기술입니다. 지난 15년간 저는 광산 채굴 지역, 고속도로 노반 침하, 터널 주변 지반 변형, 해안 지역 지반 침하 프로젝트에서 InSAR 데이터를 GNSS 측량 및 수준측량과 함께 활용해왔습니다. SAR 변형 맵핑 기술은 수십 킬로미터 규모의 지역을 동시에 모니터링할 수 있어 기존 지점 계측기보다 지반침하 분포를 훨씬 명확히 파악할 수 있습니다.
InSAR 기술의 핵심은 동일 지역에 대해 여러 시점에 촬영한 SAR 위성 영상의 위상차(Phase Difference)를 분석하여 지표 변형을 밀리미터 단위로 감지한다는 점입니다. 수평 해상도 5~50m, 수직 정확도 ±5~20mm 범위에서 지반침하를 정량화할 수 있으며, 이는 광산 침하 모니터링, 유류 저장 시설 주변 침하 감시, 고속도로·철도 노반 변형 추적에 매우 유용합니다.
InSAR 지반침하 모니터링의 기본 원리
간섭계(Interferometry) 기술의 작동 원리
2010년 중국 산시성 석탄광산 침하 프로젝트에서 처음 InSAR를 현장에 도입했을 때, 지표 변형의 공간적 분포를 수십 개의 침하계 데이터만으로 파악하려는 어려움이 있었습니다. InSAR는 이 문제를 근본적으로 해결합니다. 두 개의 SAR 위성 영상이 동일 지점을 다른 시간에 촬영할 때, 지표가 움직였다면 그 위상차가 변형량에 정비례합니다. 이를 수식으로 표현하면:
ΔΦ = (4π/λ) × d
여기서 ΔΦ는 위상차, λ는 마이크로파 파장(센티미터 단위), d는 선단거리(Line-of-Sight, LOS) 방향 변형량입니다.
예를 들어 센티미터 대역(C-band) SAR는 5.6cm의 파장을 가지므로, 위상이 2π만큼 변하면 약 2.8cm의 변형을 의미합니다. 이를 이용하면 광산 채굴 지역의 광역 침하 경향을 ±10mm 정확도로 맵핑할 수 있습니다.
위상 언래핑(Phase Unwrapping)과 대기 보정
InSAR 처리의 가장 까다로운 단계는 위상 언래핑입니다. 원시 간섭계 데이터는 -π에서 π 범위의 위상값으로 표현되므로, 실제 변형량을 복원하려면 2π 모호성을 해결해야 합니다. 제가 2018년 강원도 석탄 광산 침하 감시에 적용한 프로젝트에서는 SNAPHU(Statistical-Cost Network-Flow Algorithm for Phase Unwrapping) 알고리즘을 사용하여 ±15mm 정확도의 침하 맵을 생성했습니다.
또한 대기 효과 보정이 필수적입니다. SAR 신호는 대기권을 통과하면서 지연되며, 특히 적운이나 습도 변화가 있으면 위상오차가 ±20~30mm 수준까지 발생할 수 있습니다. 다중 기준점(Multiple Reference Point, MRP) 방법이나 대기 모델링을 통해 이를 보정합니다.
SAR 변형 맵핑 기술과 데이터 처리
시계열 InSAR(Time-Series InSAR) 분석
단일 간섭계 쌍만으로는 잡음이 많아 신뢰성 있는 지반침하 InSAR 분석이 어렵습니다. 따라서 2010년대 후반부터 시계열 InSAR 기법이 표준화되었습니다. Persistent Scatterer InSAR(PS-InSAR)와 Small Baseline Subset(SBAS) 방식이 대표적입니다.
PS-InSAR은 건물, 암반 노출, 콘크리트 포장도로 등 높은 일관성을 지닌 산란체(Persistent Scatterer, PS)를 추적합니다. 15년간의 경험상 도시 지역에서는 PS 밀도가 km²당 300~1,000개에 이르러 ±3~5mm 정확도의 침하 추적이 가능합니다. 반면 녹지나 습지 지역은 PS 밀도가 매우 낮아 보정이 어렵습니다.
SBAS-InSAR은 작은 기준선 쌍(Small Baseline Pair)들을 네트워크로 연결하여 대기 오차를 효과적으로 제거합니다. 2020년 부산 해안 지반침하 모니터링 프로젝트에서 SBAS를 적용했을 때, 월별 침하 시계열 데이터를 ±8mm 정확도로 확보할 수 있었습니다.
다중 기선(Multi-Baseline) 기술과 높이 오차 보정
| 항목 | PS-InSAR | SBAS-InSAR | 하이브리드(Unified) | |------|----------|-----------|------------------| | 적용 지역 | 도시/암반 | 농경지/습지 | 혼합 지역 | | 시간 해상도 | 중간(간헐적) | 높음(균등) | 높음 | | 수직 정확도 | ±3~5mm | ±8~12mm | ±5~8mm | | 필요 SAR 장면 수 | 30+ | 50+ | 80+ | | 처리 난도 | 중간 | 높음 | 매우 높음 | | 가성 신뢰도 | 높음 | 중간 | 높음 |
높이(Baseline) 오차는 SAR 궤도 미정확성으로 인해 발생합니다. 예를 들어 기준선 오차가 200m이면 ±15mm의 높이 오차가 누적됩니다. 현장에서는 고정밀 DEM(Digital Elevation Model)과 GNSS 기준점을 활용한 보정이 필수입니다.
지반침하 InSAR 분석 방법론
침하 경향성 분석 및 속도 추정
지반침하 InSAR 분석의 첫 단계는 선형 추세(Linear Trend) 추출입니다. 2015년 인천 항만 지역 침하 감시에 적용한 사례를 보면, 36개월의 SAR 시계열 데이터에서 월평균 침하 속도를 ±1.5mm/월 정확도로 산출했습니다.
침하 속도 계산식: v = Δh / Δt (mm/month)
여기서 Δh는 누적 침하량, Δt는 관측 기간입니다. 선형성이 낮으면(R² < 0.7) 비선형 침하 가능성을 검토하고, 수위 변동, 지하수 대수층 압력 변화, 광산 채굴 진도 등과 상관관계를 분석합니다.
공간 변형률(Strain Rate) 맵핑
침하가 균일하지 않을 때는 공간 변형률을 산출합니다. 인접 격자점 사이의 침하 경사도를 분석하면 침하 집중 지역과 안정 지역을 명확히 구분할 수 있습니다.
변형률(Strain) = (침하 경사도) × (거리)
2019년 대구 광산 폐광 지역 안정성 평가에서 이 방법을 적용했을 때, 침하 집중도(Subsidence Concentration Index)가 높은 지역에서 추후 침하가 심화되는 경향을 선제적으로 감지할 수 있었습니다.
수직 변형 분리와 수평 변형 감지
InSAR는 선단거리(LOS) 방향 변형을 측정하므로, 위성 입사각(Incidence Angle)에 따라 수직·수평 성분을 분리해야 합니다. 입사각 θ일 때:
Vertical = LOS / cos(θ) Horizontal = LOS × tan(θ)
일반적으로 센티미터 대역 위성의 입사각은 35~45°이므로, 수직 침하 10mm는 LOS 방향으로 약 12~14mm로 관측됩니다. 여러 궤도 방향의 SAR 데이터(오름차순·내림차순)를 함께 처리하면 수평·수직 변형을 동시에 추정할 수 있습니다.
현장 적용 사례와 제약 조건
광산 침하 감시 실무
2013년부터 현재까지 중국과 한국의 20개 이상 광산 침하 프로젝트에서 InSAR를 주요 모니터링 기법으로 운영했습니다. 지표면 침하 외에도 채굴 영향권 범위를 정확히 파악할 수 있어 인프라 손상 예측에 매우 유용합니다.
광산 침하 사례:
도시 기반시설 변형 모니터링
2021년 서울 강남역 대심도 복합역사 주변 지반 침하 감시에 InSAR를 적용했을 때, ±5mm 정확도로 월별 침하 추이를 파악할 수 있었습니다. Leica Geosystems의 정밀 GNSS 기준점과 함께 사용하여 InSAR 절대값 보정을 실시했습니다.
InSAR 적용의 제약 조건
1. 간섭 간섭성(Coherence) 열화
2019년 춘천 습지 지역 지반침하 조사에서 여름철 식생 우거짐으로 인해 간섭성이 0.3 이하로 떨어져 데이터 신뢰성이 크게 저하되었습니다.
2. 궤도 오차 및 기하학적 제약
3. 일시적 간섭(Temporal Decorrelation)
InSAR 모니터링 장비 및 소프트웨어
SAR 위성 데이터 소스
| 위성 | 주파수대역 | 공간해상도 | 재촬영 주기 | 데이터 정책 | |------|----------|----------|---------|----------| | Sentinel-1A/B | C-band | 5~20m | 6~12일 | 오픈 액세스 | | ALOS-2 | L-band | 3~10m | 14일 | 유상 | | Radarsat-2 | C-band | 3~100m | 24일 | 유상 | | PAZ | X-band | 1~3m | 11일 | 유상(프리미엄) | | TanDEM-X | X-band | 1~3m | 11일 | 엔터프라이즈급 |
현장에서 가장 광범위하게 사용하는 것은 Sentinel-1입니다. 유럽우주청(ESA)의 Copernicus 프로그램을 통해 무료로 제공되며, 12일 반복주기로 전 지구 커버리지를 확보합니다. 다만 수평 해상도가 20m 정도로 제한적이어서, 소규모 지역에서는 ALOS-2나 PAZ의 고해상도 데이터가 필요합니다.
InSAR 처리 소프트웨어
오픈소스 소프트웨어:
상용 소프트웨어:
제가 가장 많이 사용한 조합은 SNAP + StaMPS + QGIS입니다. 비용 측면에서 학계·중소 기업에 유리하고, 처리 과정의 세밀한 제어가 가능합니다.
검증 및 보정 도구
InSAR 결과의 신뢰성 확보를 위해 반드시 외부 데이터와의 교차 검증이 필요합니다:
1. GNSS 기준점 보정
2. 수준측량 단면
3. 침하계(Extensometer) 연속 모니터링
2026년 InSAR 기술 최신 동향
고정밀 DEM 활용 및 InSAR-INSAR 스태킹
2024년부터 무료로 제공되는 고해상도 DEM(TanDEM-X 12m, NASA SRTM+)이 기준선 오차 보정 정확도를 크게 향상시켰습니다. 결과적으로 수직 정확도가 기존 ±8~12mm에서 ±5~8mm로 개선되었습니다.
또한 여러 궤도 방향의 SAR 데이터를 활용한 3D 변형 분해(3D Deformation Decomposition) 기술이 상용화되었습니다. 오름차순(Ascending)·내림차순(Descending) 궤도 데이터를 동시에 처리하면 수직·수평(동서, 남북) 3성분을 분리할 수 있어 구조물 기울임, 측방 이동 등을 감지할 수 있습니다.
머신러닝 기반 간섭성 개선
2025년부터 딥러닝을 활용한 위상 노이즈 제거 기술이 실용화되었습니다. 기존 골덱(Goldstein) 필터링만으로는 제거하기 어려운 대기 오차를 CNN(Convolutional Neural Network)으로 학습·제거하여 신호 품질을 10~20% 향상시킬 수 있게 되었습니다.
실시간 모니터링 시스템 구축
2026년 현재 국내 주요 광산·SOC 시설에서는 월 1회 SAR 촬영 데이터를 3~4일 내에 자동 처리·배포하는 시스템이 운영 중입니다. 클라우드 기반 InSAR 플랫폼(Trimble INSAR, Capella Space 등)이 대폭 확산되면서, 고비용의 장비 투자 없이도 엔터프라이즈급 서비스를 이용할 수 있게 되었습니다.
대역외(Out-of-band) 위성 활용
X-band(3cm) 초고해상도 SAR(PAZ, ICEYE, Capella)의 본격 운영으로 1~3m 해상도의 침하 맵 생성이 가능해졌습니다. 다만 데이터 비용이 전문가 등급(Professional Grade)~엔터프라이즈급이므로, 광범위 영역 모니터링보다는 고가치 시설(원전, 댐, 주요 교통 인프라) 주변 정밀 감시에 활용됩니다.
다중 주파수대역 융합(Multi-band Fusion)
C-band와 L-band 데이터를 동시에 처리하여 대기 오차를 효과적으로 제거하는 기법이 표준화되었습니다. 예를 들어 Sentinel-1(C-band) + ALOS-2(L-band) 조합은 수직 정확도를 ±3~5mm 수준으로 끌어올릴 수 있으며, 특히 습한 지역에서의 성능 개선이 두드러집니다.
자주 묻는 질문
Q: InSAR 지반침하 모니터링이 기존 GNSS 측량보다 정확한가?
InSAR 정확도는 ±5~20mm, GNSS 정확도는 ±3~10mm입니다. 절대값 정확도는 GNSS가 우수하지만, InSAR는 광역 공간 분포를 동시에 파악할 수 있어 침하 원인 규명에 더 유리합니다. 현장에서는 두 기법을 병용합니다.
Q: 녹지 지역에서 InSAR 간섭성이 떨어지는 이유는?
식생 우거짐으로 인해 SAR 신호가 나뭇가지·잎에서 산란되어 신호 일관성이 감소합니다. 특히 여름철 활엽수 지역은 간섭성 0.3 이하로 떨어져 데이터 신뢰성이 급격히 감소합니다. 이 경우 L-band(ALOS-2) 데이터나 겨울철 촬영 데이터 활용이 효과적입니다.
Q: 지하수 저장소 침하를 InSAR로 감지할 수 있나?
가능합니다. 미국 캘리포니아 San Joaquin Valley의 관개용 대수층 침하는 InSAR로 1m 이상의 침하를 정량화한 바 있습니다. 월 5~10cm 이상의 침하가 있으면 ±10mm 정확도로 감지할 수 있습니다.
Q: InSAR 데이터 처리에 얼마나 오래 걸리나?
단일 간섭계 쌍 처리는 2~3일, PS-InSAR나 SBAS 시계열 처리는 100장 이상의 SAR 장면을 포함할 경우 1~2주가 소요됩니다. 클라우드 기반 자동 처리 서비스 이용 시 3~5일 내 결과 확보가 가능합니다.
Q: InSAR 결과의 신뢰성을 어떻게 검증하나?
GNSS 기준점(3~5개), 수준측량 단면(200~500m 간격), 침하계(심부층별) 데이터와의 교차 검증이 필수입니다. 수직 정확도가 ±10mm 이내이면 모니터링용으로 신뢰할 수 있습니다. 또한 간섭성 맵에서 낮은 신뢰도 지역(0.4 이하)은 별도 표기하여 사용자가 해석할 때 주의하도록 해야 합니다.
