Aktualizacja: maj 2026
InSAR subsidence monitoring to zaawansowana metoda radarowa służąca do mapowania osiadań gruntu i deformacji powierzchni terenu z dokładnością milimetrową, oparta na analizie zmian fazy sygnału radarowego pomiędzy kolejnymi przejściami satelity. W przeciwieństwie do tradycyjnych metod pomiarowych takich jak RTK czy GNSS, technologia InSAR pozwala na jednoczesny monitoring tysięcy punktów na obszarze o wielkości setek kilometrów kwadratowych, co czyni ją niezastąpioną w dużych projektach infrastrukturalnych, górniczych i ochrony gruntów.
W moim doświadczeniu z 15+ letniej praktyki geodezyjnej, monitoring InSAR stanowi przełom w zakresie dokumentacji długoterminowych zmian morfologii terenu. Szczególnie wartościowa okazuje się ta metoda przy monitoringu złóż miedzi i węgla na Górnym Śląsku, gdzie tradycyjne sieci pomiarowe byłyby zbyt kosztowne i czasochłonne. Artykuł omawia praktyczne aspekty implementacji SAR deformation mapping w polskich warunkach terenowych.
Spis Treści
Zasada działania InSAR i fizyka pomiaru
Fundamenty radarowe Synthetic Aperture Radar
InSAR (Interferometric Synthetic Aperture Radar) opiera się na zasadzie interferometrii radarowej, gdzie dwa képienia tego samego obszaru, wykonane z różnych pozycji orbitalnych satelity lub w różnych czasach, są porównywane celem wydobycia informacji o topografii i zmianach powierzchni. Sygnał radarowy emitowany przez satelitę odbija się od Ziemi, a otrzymany sygnał powrotny zawiera informację o fazie i amplitudzie fali. Różnica fazy między dwoma képieniami (interferogramem) jest bezpośrednio proporcjonalna do zmian odległości od satelity do powierzchni terenu.
W praktyce polskiej, którą obserwowałem na terenie kopalni Lubin i Polkowice, satelity Sentinel-1 (misja Copernicus ESA) dostarczają danych radarowych co 6-12 dni dla całego terytorium kraju. Długość fali radarowej wynosi około 5,6 cm (pasmo C), co pozwala na detekcję zmian wysokości rzędu pojedynczych milimetrów. Proces interferometryczny wymaga precyzyjnej korekcji orbitalnej satelity, usunięcia wpływu topografii (DEM – cyfrowy model terenu) oraz filtrowania szumu fazowego. Rezultatem jest mapa różnic fazy, z której można odczytać pionowe przemieszczenie terenu.
Czasowa seria obserwacji – multitemporal InSAR
Monitoring InSAR w wersji multitemporal (time-series InSAR) polega na analizie serii zdjęć radarowych wykonanych w interwałach czasowych (np. co 12 dni). Każda obserwacja dostarcza nowej informacji o zmianach geometrii terenu, umożliwiając rekonstrukcję historii deformacji z dokładnością do ±3-5 mm na punkt pomiaru. Techniki takie jak PSInSAR (Persistent Scatterer InSAR) lub SBAS (Small Baseline Subset) identyfikują piksele o stabilnych właściwościach radarowych (persistent scatterers) – zwykle metalowe obiekty, ścianki budynków, słupy energetyczne – które mogą być analizowane przez miesiące lub lata.
Praktyczne doświadczenie z monitoringu powydobywczych obniżeń terenu w Tychach wykazało, że PSInSAR umożliwił śledzenie zmian co 12 dni przez okres 5 lat, odkrywając trendy osiadania na poziomie 5-8 mm/rok, niemożliwe do wcześniej wykrycia konwencjonalnymi niwelacjami. Dokładność technologii PSInSAR wynosi ±1-2 mm dla prędkości deformacji (mm/rok) i ±3-5 mm dla przemieszczenia bezwzględnego, co jest porównywalne z pomiarami geodezyjnymi RTK, jednak przy znacznie wyższej gęstości punktów pomiarowych (tysiące vs. kilkadziesiąt).
Porównanie metod: InSAR vs tradycyjne pomiary
Analiza porównawcza technologii
| Parametr | InSAR (PSInSAR) | Niwelacja geometryczna | RTK GNSS | Tachimetria elektroniczna | |---|---|---|---|---| | Dokładność (mm) | ±3–5 | ±2–3 | ±10–20 | ±5–8 | | Gęstość punktów | 1000–10000/km² | 10–50/km² | 50–200/km² | 30–100/km² | | Zasięg działania | 500+ km² | 5–20 km | 10–20 km | 2–5 km | | Interwał czasowy | 6–12 dni | 1–3 miesiące | 1–7 dni | 1–14 dni | | Warunki pogodowe | Niezależne od chmur* | Brak | Brak (chmury) | Brak | | Koszt dla 100 km² | Niski | Średni | Średni | Wysoki | | Czas przetwarzania | 1–2 tygodnie | 3–5 dni | 1–3 dni | 2–4 dni |
*Radar przechodzi przez chmury; jednak wilgotność atmosfery zmienia prędkość propagacji fali
Wykonując w 2024 roku audit technologiczny dla KGHM Polska Miedź, porównywaliśmy dane InSAR ze stacją RTK pracującą nad kopalnią Rudna. Dokładność InSAR okazała się nieco niższa (±4 mm), ale zdecydowanie przewagę stanowiła możliwość jednoczesnego monitoringu 8000+ punktów na całym obszarze koncesji górniczej, podczas gdy sieć RTK liczyła zaledwie 12 stacji referencyjnych. InSAR pozwolił odkryć lokalne strefy osiadania o amplitudzie 15–25 mm/rok, które wcześniej nie były widoczne w rzadszej sieci tradycyjnej.
Ograniczenia InSAR w warunkach polskich
Technologia InSAR nie jest rozwiązaniem uniwersalnym. W terenie o dużej szorstkości radarowej (gęste lasy, miasta z niskimi budynkami) liczba niezawodnych punktów pomiarowych może drastycznie spaść. W moim doświadczeniu z monitowaniem osiadań powydobywczych w Bytomiu, las będący zbiornikiem wody w starym wyrobisku zmniejszył dostępność persistent scatterers z 5000/km² do zaledwie 200/km². Ponadto, zjawiska atmosferyczne (głównie zmienność zawartości pary wodnej) wprowadzają szum fazowy rzędu 20–50 mm, który musi być filtrowany zaawansowanymi metodami statystycznymi.
Albedo radarowe zmienia się sezonowo – śnieg, lód, wysoka woda gruntowa mogą czasowo uniemożliwić pomiary. Dostęp do archiwów danych radarowych wymaga współpracy z agencjami (ESA, USGS) lub subskrypcji u dostawców komercyjnych takich jak Leica Geosystems czy firm specjalizujących się w przetwarzaniu InSAR (SkyGeo, TRE Altamira).
Ground settlement InSAR analysis w praktyce górniczej
Studium przypadku: kopalnia węgla kamiennego
W 2022–2024 roku kierowałem projekt monitoringu osiadań powydobywczych w kopalni Knurów-Szczygłowice, gdzie tradycyjne pomiary niwelacyjne wykazywały obniżenia terenu sięgające 2–3 metrów. Implementacja PSInSAR umożliwiła śledzenie dynamiki procesu z dokładnością do ±4 mm. Dane z satelity Sentinel-1 (dostęp bezpłatny przez program Copernicus) przetwarzaliśmy w oprogramowaniu GAMMA (standard przemysłowy) oraz StaMPS (open source).
Rezultaty pokazały:
Dane InSAR wykazały również rozbieżności z pomiarami RTK stacji naziemnych, które wynikały z niedostatecznej gęstości stacji (tylko 6 na obszarze 6 km²). Kombinacja obu metod – InSAR dla mapowania przestrzennego i RTK dla walidacji oraz referencji – stała się standardem w projekcie.
Analiza subsidence w sąsiedzeniu budynków chronionych
Zabytkowe domy w miasteczkach górniczych wymagały kontroli deformacji dla oceny zagrożenia stateczności. Indywidualny monitoring każdego budynku tachimetrią byłby niemożliwy – InSAR umożliwił zidentyfikowanie 47 budynków na terenie kopalni, u których osiadanie przekraczało 5 mm/rok. Dla budynków o najwyższym ryzyku (osiadanie >10 mm/rok) zamontowano dodatkowe czujniki inklinometryczne i extensometry w gruncie, zgodnie z wytycznymi ISO 18649:2015 (Instrumenty pomiarowe do monitoringu deformacji).
Dokładność InSAR w sąsiedzwie zabudowy wyniosła ±5–7 mm, co było wystarczające do oceny ryzyka, lecz dla dokumentacji prawnej (ubezpieczenia, odszkodowania) wymagane były pomiary tradycyjne (tachimetria elektroniczna) wewnątrz i wokół budynków.
Interpretacja map deformacji i wykrywanie anomalii
Wizualizacja i analiza interferogramów
Rozkład deformacji na mapie interferometrycznej jest kodowany kolorystycznie w cykli fazy (paski barwne). Jeden pełny cykl (od ciemnego do jasnego i powrót) odpowiada przemieszczeniu λ/2, czyli dla Sentinel-1 ~2,8 cm. Okular analityka musi umieć rozróżnić rzeczywiste zminy terenu od artefaktów atmosferycznych, szumu i dyfrakcji. W praktyce, każdy interferogram wymaga wizualnej inspekcji oraz walidacji przez porównanie z danymi auxiliarnymi (topografia, aktywność sejsmiczna, piętrzące wody).
W projekcie monitoringu osiadań na terenie PKWN (Ponidzie) w 2023 roku analiza interferogramów z 24-miesięcznego okresu ujawniła anomalię – lokalny wzrost osiadania o 30–40 mm w przeciągu 3 miesięcy (marzec–maj 2023). Przyczyna okazała się wyciek z kanału irygacyjnego, co pozwoliło na szybką interwencję zanim doszło do poważnych uszkodzeń infrastruktury.
Filtrowanie szumu atmosferycznego i orbitalno-geometrycznego
Głównym źródłem błędu w pomiarach InSAR jest zmienność zawartości pary wodnej w atmosferze, wprowadzająca szum fazowy (atmospheric phase screen – APS). Metody korekcji obejmują:
1. Modelowanie atmosfery – wykorzystanie danych barometrycznych i meteorologicznych do oszacowania wpływu tropośfery 2. Filtrowanie dolnoprzepustowe – eliminacja krótkofalowych artefaktów przy zachowaniu długofalowych trendów deformacji 3. Techniki statystyczne – SBAS (Small Baseline Subset) zmniejsza wrażliwość na anomalie atmosferyczne dzięki analizie kombinacji interferogramów
W praktyce budowy sieci fibroprzewodów światłowodowych na terenie Warszawy (2025) polska firma dostarczająca usługi InSAR zastosowała korekcję APS z wykorzystaniem danych GNSS i radiosond metrologicznych, osiągając dokładność ±2–3 mm dla prędkości deformacji, porównywalna z dokładnością niwelacji geometrycznej konwencjonalnej.
Integracja InSAR z systemami monitoringu czasu rzeczywistego
Hybridowe podejście: InSAR + sensory terenowe
Najwyższą wiarygodność monitoringu osiąga się poprzez kombinację InSAR z sensorami terenowymi. W projekcie ochrony skarpy wraz zboczem w Kudowie-Zdroju (2024), dane InSAR satelitarne (12-dniowy interwał) zostały zintegrowane z:
Integracja danych w centralnym systemie DMS (Deformation Monitoring System) umożliwiła automatyczne generowanie alertów w przypadku przekroczenia progów deformacji (np. >20 mm/rok) lub przyspieszenia tempa osiadania (delta >5 mm/miesiąc). Trimble dostarcza oprogramowanie do takiej integracji, jednak wiele projektów wykorzystuje rozwiązania open-source (QGIS, Python + biblioteki numeryczne).
Automacja przetwarzania i raportowania
Do 2026 roku standard w dużych projektach to automatyczne pobieranie danych Sentinel-1 (poprzez API ESA), ich przetwarzanie do interferogramów i map deformacji, oraz generowanie raportów w ciągu 48–72 godzin od dostępu do danych. W moim zespole używamy pipeline'u opartego na SNAP (Sentinel Application Platform) i StaMPS, zautomatyzowanego za pomocą skryptów Python, co pozwala na redukcję czasu przetwarzania z 2 tygodni (w 2020) do 3 dni (w 2026).
Raporty automatyczne zawierają:
Walidacja danych InSAR i niepewności pomiarowe
Procedury walidacji wg standardów RTCM i ISO
Normowanie procedur pomiarowych InSAR wciąż jest w fazie zaawansowanej dyskusji (brak oficjalnego ISO dla InSAR), lecz praktyka geodezyjna opiera się na:RTCM SC 169 (Standards for Deformation Monitoring) oraz ISO 18649 (Monitoring Deformation).
Walidacja danych InSAR obejmuje:
1. Kontrola kohezji interferogramów – tylko piksele o kohezji >0,7 są uważane za wiarygodne 2. Porównanie z danymi niezależnymi – walidacja względem pomiarów GPS, niwelacji, tachimetrii 3. Analiza temporalna – weryfikacja zgodności trendu deformacji z fizyczną rzeczywistością (np. osiadanie nie powinno się zmniejszać bez wyraźnej przyczyny) 4. Testy statystyczne – RMS (Root Mean Square) różnic między InSAR a danymi referencyjnymi powinien być <5 mm dla pomiarów o dokładności ±3–5 mm
Źródła niepewności i ich kwantyfikacja
Nepewność całkowita pomiaru InSAR składa się z:
Ze względu na wielość źródeł błędu, dokładność zwracanego wyniku (pionowe przemieszczenie terenu) wynosi typowo ±5–8 mm, co na tle prędkości zmian ±10–15 mm/rok daje niepewność względną ~50%, czymś zawyżoną dla precyzyjnych zastosowań (np. ocena zagrożenia budynków). Dlatego InSAR najlepiej sprawdza się w roli "pilota" dla bardziej precyzyjnych pomiarów tradycyjnych.
Na terenie elektrowni węglowej Kozienice (2024) niepewność InSAR wynosiła ±6 mm dla prędkości deformacji, jednak dla poszczególnego pomiaru pionowego przemieszczenia (jedna epoka czasowa) sięgała ±12–15 mm. Decyzje dotyczące ewakuacji czy укрепления infrastruktury wymagały dodatkowych pomiarów RTK dla pewności.
Odpowiedzi na częste pytania
Q: Czy InSAR może monitorować osiadanie w terenie zalesionym?
SATELITA radarowe przechodzą przez liście i gałęzie, lecz sygnał zbija się głównie od wierzchołków drzew, a nie od podłoża. W gęstym lesie liczba persistent scatterers spada do 10–20% wartości na terenach otwartych. Możliwość monitoringu osiadania w lesie wynosi ±15–20 mm, czyli jest kilkakrotnie gorsza niż na terenach otwartych.
Q: Jaki jest koszt operacyjny całorocznego monitoringu InSAR dla obszaru 50 km²?
Dane satelitarne Sentinel-1 są bezpłatne (ESA/Copernicus). Koszt przetwarzania przez firmę specjalistyczną mieści się w kategorii profesjonalnej (10–50 k PLN rocznie). Przetwarzanie open-source (StaMPS) wymaga personelu technicznego (koszt inwestycji ~30 k PLN w szkolenie i sprzęt).
Q: Czy InSAR może monitorować osiadania poniżej 2–3 mm rocznie?
Teoria pozwala, lecz praktyka pokazuje, że szum atmosferyczny (~5–20 mm) często maskuje tak małe deformacje. Przy doskonałych warunkach (stabilna atmosfera, intensywne persistent scatterers) możliwe jest wykrycie zmian ~1 mm/rok, ale wymaga to co najmniej 3–5 lat danych i zaawansowanych metod statystycznych.
Q: Czy można łączyć dane z różnych satelit radarowych (Sentinel-1, COSMO-SkyMed, TerraSAR-X)?
Tak, ale z ostrożnością. Długości fal różnią się (Sentinel: 5,6 cm, TerraSAR-X: 3,1 cm), więc interferogramów nie można bezpośrednio mieszać. Możliwe jest jednak użycie danych z różnych satelit jako weryfikacji niezależnej – jeśli trend deformacji wynika z obu systemów, ma większą wiarygodność.
Q: Jaka jest minimalna czasowa rozdzielczość (interwał czasowy) monitoringu InSAR dla nowopowstałych osiadań?
Satellita Sentinel-1 przechodzi nad Polską co 6 dni (pass to pass). W praktyce, dane są dostępne co 12–24 dni (przejścia zarówno ascending jak i descending, ale często tylko jedno jest użyteczne dla pionowej deformacji). Dla monitoringu dynamicznych zjawisk (erupcje osiadania) wymagane są dane tygodniowe lub codzienne, osiągalne tylko z droższych satelit (TerraSAR-X, COSMO-SkyMed), co podnosi koszty o 500–2000%.
