Materiał przygotowany z pomocą AI

Treść tego wpisu została wygenerowana z wykorzystaniem sztucznej inteligencji i nie została jeszcze zweryfikowana merytorycznie.

Ten materiał jest częścią cyklu Modelowanie odpływu ze zlewni

Hydrologiczne modele zlewni

Wprowadzenie

Zlewnia jako system hydrologiczny to obszar, w którym zachodzi dynamiczny proces transformacji zasilania (opadu) w odpływ. Wymaga to traktowania zlewni jako systemu złożonego z różnych elementów składowych, takich jak przepływ wody oraz procesy dynamiczne odpowiedzialne za transformację opadu w odpływ.

Każdy system przyrodniczy ma swoje otoczenie, z którym wchodzi w interakcje. W przypadku systemów hydrologicznych kluczowym elementem otoczenia jest atmosfera, która ma wpływ na obieg wody. Zakwalifikowanie konkretnego obiektu do systemu lub otoczenia zależy od przedmiotu i celu badań.

Elementy systemu hydrologicznego

Każdy element systemu opisuje się za pomocą zbioru cech (atrybutów), które mogą być stałe lub zmienne w czasie. Zależności między cechami poszczególnych elementów determinują:

  • strukturę wewnętrzną systemu
  • interakcje z otoczeniem

Współrzędne stanu systemu — wielkości charakteryzujące zachowanie systemu.

Stan systemu — najmniejszy liczebnie zbiór współrzędnych stanu, określonych dla danej chwili.


Modele hydrologiczne — definicja i zastosowanie

Modele hydrologiczne stanowią formę prezentacji obiektu, procesu lub systemu (zlewnia). Wykorzystywane są do opisu obiegu wody w zlewniach poprzez zbiór równań matematycznych, służących opisowi wyróżnionych cech badanego obiektu.

Zastosowania modeli hydrologicznych:

  • pozyskiwanie informacji hydrologicznej na potrzeby projektowania
  • narzędzie stosowane w prognozach hydrologicznych
  • głównie stosowane w zlewniach niekontrolowanych

Równanie ogólne modelu hydrologicznego

Jeżeli przez wyrazi się zmienne wejściowe do modelu hydrologicznego, a przez zmienne wyjściowe, wówczas modele hydrologiczne można opisać za pomocą równania:

gdzie:

  • — funkcja o założonej postaci
  • — parametry modelu (zlewni) pozyskane z pomiaru oryginału
  • — błąd modelu

Klasyfikacja modeli hydrologicznych

Przeprowadzenie klasyfikacji modeli hydrologicznych jest uzależnione od podejścia, wytypowania odpowiednich cech i kryteriów klasyfikacyjnych (por. klasyfikacja modeli). Modele można podzielić według następujących kryteriów:

  1. Ogólna struktura modelu oraz przedmiot modelowania
  2. Rola czynnika czasu
  3. Jednoznaczność opisu i uzyskiwanych wyników
  4. Stosowane podejścia i metody rozwiązania
  5. Właściwości funkcji operatora

1. Klasyfikacja według struktury modelu i przedmiotu modelowania

Modele pojedynczych procesów — odnoszą się do jednorodnych lub prawie jednorodnych zjawisk, które można wyrazić za pomocą jednego procesu:

  • parowanie
  • infiltracja
  • spływ powierzchniowy

Modele komponentów zlewni — opisujące zespół procesów zachodzących w poszczególnych środowiskach występowania wody w zlewni:

  • modele strefy aeracji
  • modele strefy saturacji
  • modele pokrywy śnieżnej

Modele integralne — obejmujące wszystkie komponenty zlewni, stanowiące logiczny system wzajemnych powiązań i oddziaływań wszystkich procesów hydrologicznych zachodzących w zlewni. Ważnym elementem jest prawidłowość struktury chronologicznej komponentów i procesów. Model integralny jest docelowym w badaniach hydrologicznych.

Modele globalne — opisują zlewnię jako całość; za pomocą jednego operatora (funkcji) symuluje się łączny efekt wszystkich procesów.


2. Klasyfikacja według roli czynnika czasu

Modele statyczne — opisują zjawiska i powiązania między nimi w pewnym momencie lub okresie przyjętym za jednostkę. Przykład: zależności regresyjne, wielowymiarowe zależności przyczynowo-skutkowe.

Wzór Punzeta (1981) na jednostkowy średni niski odpływ wód z karpackiego dorzecza Wisły:

gdzie:

  • — średnia wysokość zlewni [m n.p.m]
  • — opad średni roczny w zlewni [mm]
  • — spadek podłużny cieku [m·km⁻¹]
  • — wskaźnik nieprzepuszczalności gleb Bołdakowa [%]

Modele dynamiczne — opisują zjawiska i występujące między nimi relacje w funkcji czasu.

Równanie Hortona opisujące infiltrację w podłożu glebowym:

gdzie:

  • — natężenie infiltracji [mm·h⁻¹]
  • — końcowe natężenie infiltracji [mm·h⁻¹]
  • — początkowa wartość infiltracji równa natężeniu opadu [mm·h⁻¹]
  • — stała recesji
  • — czas

3. Klasyfikacja według jednoznaczności opisu

Modele deterministyczne — opisują w sposób jednoznaczny system, na podstawie których można jednoznacznie (z prawdopodobieństwem ) określić sygnał wyjściowy.

Modele probabilistyczne — nie opisują badanego zjawiska w sposób jednoznaczny i na ich podstawie nie potrafimy jednoznacznie przewidzieć sygnału wyjściowego.

Uwaga: W rzeczywistości modele systemów hydrologicznych są zawsze modelami probabilistycznymi.


4. Klasyfikacja według stosowanych podejść i metod rozwiązania

W hydrologii inżynierskiej, przy konstruowaniu modeli, stosuje się trzy podejścia metodyczne:

a) Modele hydrodynamiczne “białej skrzynki”

Opracowywane przy pełnej informacji o działaniu i wewnętrznej budowie modelowanego systemu. Wszystkie procesy są opisane za pomocą równań fizyki matematycznej.

Charakterystyka:

  • bazują na równaniach fizyki matematycznej opisujących zasadę zachowania masy, pędu i energii
  • w praktyce często trudne do przyjęcia ze względu na:
    • niedostateczną znajomość fizycznych praw zachowania się systemu
    • trudności z określeniem warunków początkowych i brzegowych
    • niejednorodność zlewni
b) Modele systemowe “czarnej skrzynki”

Opracowywane przy całkowitym braku informacji o wewnętrznej budowie modelowanego systemu i jego działaniu. Informacja jest uzyskiwana na podstawie jego sygnałów wejściowych i wyjściowych.

Charakterystyka:

  • pomija się prawa fizyki rządzące ruchem wody w zlewni
  • wnioski dotyczące funkcjonowania systemu są wyciągane na podstawie sygnału wejściowego (opad) oraz wyjściowego (odpływ)
c) Modele koncepcyjne “szarej skrzynki”

Oparte na częściowej znajomości wewnętrznej budowy modelowanego systemu i jego działaniu.

Charakterystyka:

  • próba zapełniania luki pomiędzy dwoma powyższymi podejściami
  • znajdują szerokie zastosowanie w hydrologii inżynierskiej
  • nie odwzorowuje się fizycznych procesów transformacji opadu w odpływ, tylko na wyjściu uzyskuje się taki wynik, który jest efektem rzeczywistego odziaływania zlewni na transformacje opadu w odpływ
  • zazwyczaj opisywane przez kilka parametrów

5. Klasyfikacja według właściwości funkcji operatora

Funkcja operatora w modelach hydrologicznych to matematyczny sposób na opisanie i przekształcenie danych wejściowych (np. opadów) na dane wyjściowe (np. odpływy).

W kontekście hydrologii operator ten działa jako “czarne pudełko” przetwarzające wejścia systemu hydrologicznego na jego wyjścia. Główne elementy tej funkcji to opisanie relacji między opadem a odpływem na podstawie szeregu procesów hydrologicznych zachodzących w zlewni:

  • infiltracja
  • spływ powierzchniowy
  • retencja
  • parowanie

Funkcja operatora może reprezentować różne procesy fizyczne zachodzące w zlewni, które wpływają na końcowy wynik (np. przepływ w rzece). Może uwzględniać dynamikę procesów hydrologicznych poprzez równania różniczkowe, czasami w sposób deterministyczny, a czasem stochastyczny.

Podział ze względu na właściwości funkcji operatora:

a) Modele liniowe

Modele dynamiczne, w których spełniona jest zasada superpozycji. Reakcja systemu na sumaryczne wejście jest równa sumie reakcji na składowe wejścia z uwzględnieniem translacji w czasie (np. modele hydrogramu jednostkowego).

Układ liniowy można zapisać jako:

gdzie:

  • — wyjście z modelu
  • — składowa wejścia
  • — operator transformujący składowe wejścia w zmienne wyjścia
  • — liczba elementarnych sygnałów wejściowych

b) Modele nieliniowe

Modele dynamiczne, w których nie jest spełniona zasada superpozycji (występują nieliniowe funkcje operatorów).

Uwaga: W rzeczywistości zlewnia nigdy nie reaguje liniowo na zasilanie opadem, jednak ze względu na duże trudności obliczeniowe traktuje się ją jako model liniowy.

c) Modele stacjonarne — niestacjonarne

  • Modele stacjonarne — operatory nie ulegają zmianie w czasie
  • Modele niestacjonarne — opisują systemy, w których właściwości są zmienne w czasie

Uwaga: Zbudowanie modeli stacjonarnych jest dużo prostsze, jednak zachodzące w hydrologii zmiany sezonowe zmuszają do stosowania modeli niestacjonarnych.

d) Modele o parametrach skupionych — rozłożonych

  • Modele o parametrach skupionych — nie uwzględniają przestrzennej zmienności charakterystyk badanego systemu. Parametry sprowadza się pozornie do właściwości, jakimi charakteryzuje się jeden fikcyjny punkt systemu.

  • Modele o parametrach rozłożonych — uwzględniają przestrzenną zmienność charakterystyk systemu.


Budowa modeli hydrologicznych

Proces tworzenia modelu hydrologicznego obejmuje trzy zasadnicze etapy:

  1. Specyfikacja modelu — sformułowanie problemu i zadań, wybór struktury modelu, wybór funkcji operatorów
  2. Identyfikacja modelu — wyznaczenie parametrów modelu
  3. Weryfikacja modelu — określenie jakości pracy modelu

Etap 1: Specyfikacja modelu

Jednym z ważniejszych elementów modelu jest jego struktura, czyli sposób, w jaki jego części są ze sobą powiązane i uzależnione. Najczęściej struktury opisywane są przez schematy blokowe, uzupełniane odpowiednimi algorytmami, które opisują transformację rozważanych procesów.

Głównym elementem decydującym o wyborze struktury modelu jest:

  • jakość informacji możliwej do uzyskania
  • cel i zadanie stawiane przed modelem
Trzy zasadnicze etapy specyfikacji modelu:

a) Zdefiniowanie celu

Określenie, co ma być osiągnięte przez model:

  • prognoza przepływów w rzekach
  • analiza skutków zmian klimatycznych
  • ocena wpływu zmian użytkowania terenu
  • prognozy powodziowe
  • optymalizacja zarządzania zasobami wodnymi

b) Dobór odpowiedniego typu modelu

  • Modele koncepcyjne — oparte na uproszczonych założeniach, często stosowane do symulacji w skali zlewni
  • Modele fizyczne — oparte na szczegółowych równaniach fizycznych, opisujących rzeczywiste procesy zachodzące w hydrologii
  • Modele statystyczne lub empiryczne — oparte na analizie historycznych danych i relacjach pomiędzy zmiennymi

c) Struktura modelu

  • Ustalenie procesów hydrologicznych — wybór procesów, które będą modelowane (odpływ powierzchniowy, parowanie, infiltracja, przepływ wód podziemnych)
  • Podział zlewni na jednostki — modele często dzielą zlewnię na mniejsze obszary (siatki, strefy, zlewnie cząstkowe), w których wykonywane są obliczenia

Etap 2: Identyfikacja modelu

Identyfikacja to warunkiem możliwości opisu danego systemu naturalnego za pomocą modeli matematycznych jest określenie wartości liczbowych parametrów modelu na podstawie informacji uzyskanych w wyniku badań empirycznych.

Parametry modelu to wartości stanowiące tzw. węzeł łączący model z oryginałem. Każdy system charakteryzuje się innym zestawem parametrów.

W zależności od posiadanej informacji oraz przyjętej struktury modelu wyróżnia się dwa podejścia przy identyfikacji parametrów:

a) Bezpośredni pomiar w terenie

Stosowane w modelach białej skrzynki. Polega na zbieraniu rzeczywistych danych fizycznych z naturalnego środowiska w celu dokładniejszego odwzorowania procesów hydrologicznych w modelu.

Mierzone zmienne:

  • przepływ wód
  • opady
  • właściwości gleby
  • poziom wód gruntowych

Dane te są wykorzystywane do kalibracji modelu hydrologicznego, co pozwala na lepsze zrozumienie dynamiki wody w danym obszarze oraz na tworzenie bardziej precyzyjnych prognoz hydrologicznych.

Cel: Dopasowanie modelu do rzeczywistych warunków środowiskowych, co zwiększa jego dokładność i użyteczność w praktyce (np. w prognozowaniu powodzi, zarządzaniu zasobami wodnymi, ocenie skutków zmian klimatycznych).

b) Metody optymalizacyjne

Stosuje się je dla modeli “czarnej” bądź “szarej skrzynki” zawierające parametry będące bliżej nieokreślonymi współczynnikami liczbowymi.

Identyfikację taką prowadzi się na podstawie zarejestrowanych sygnałów wejściowych i wyjściowych systemu naturalnego. Proces optymalizacji parametrów polega na uzyskaniu maksymalnej zgodności pomiędzy elementami wyjścia z modelu a wartościami pomierzonymi.

Kryterium zgodności można zapisać jako sumę kwadratów odchyłek:

gdzie:

  • — odchyłki wartości obliczonych i pomierzonych
  • — liczebność ciągu

Uwaga: Zadania optymalizacyjne są tym prostsze, czym model zawiera mniej parametrów.


Etap 3: Weryfikacja modelu

Weryfikacja modelu hydrologicznego to proces sprawdzania, jak dobrze model odzwierciedla rzeczywiste warunki i procesy hydrologiczne. Polega na porównaniu wyników uzyskanych z modelu z danymi rzeczywistymi, które nie były używane podczas jego kalibracji.

Cel weryfikacji: ocena, czy model jest w stanie niezawodnie prognozować zjawiska hydrologiczne w różnych sytuacjach i warunkach.

Weryfikacja modelu przeprowadzana jest w dwóch etapach:

Etap I: Weryfikacja na tle zależnego materiału badawczego

Polega na przeprowadzeniu testów zgodności wartości obliczonych modelem z wartościami pomierzonymi. Obliczenia te wykonuje się na tzw. materiale zależnym (ten sam, który brał udział w procesie optymalizacji).

W zależności od przeznaczenia modelu zakładane są różne wartości błędów dopuszczalnych, decydujących o jakości modelu, a następnie decydujących o jego dopuszczeniu do eksploatacji.

Jeżeli na tym etapie nie uzyska się zadawalających wyników, należy dokonać korekty funkcji operatorów, a nawet jego struktury, aż do uzyskania zadawalających wyników.

Etap II: Weryfikacja na tle niezależnego materiału badawczego

Polega na sprawdzeniu kryterium zgodności na tzw. materiale niezależnym (taki który nie był wykorzystywany w procesach optymalizacyjnych).

Dalszy sposób postępowania taki sam jak w przypadku etapu I. Dopiero po uzyskaniu zadawalających wyników weryfikacji model można przeznaczyć do eksploatacji.


Miary jakości modelu

Do weryfikacji modeli hydrologicznych stosuje się różne miary jakościowe, które pozwalają ocenić, jak dobrze model odzwierciedla rzeczywiste dane:

1. Współczynnik Nash-Sutcliffe’a (NSE)

Mierzy efektywność predykcji modelu w porównaniu do średniej wartości obserwowanych danych.

gdzie:

  • — przepływ obserwowany [m³·s⁻¹]
  • — przepływ obliczony za pomocą modelu [m³·s⁻¹]
  • — wartość średnia z przepływów obserwowanych [m³·s⁻¹]

Zakres: od minus nieskończoności do 1, gdzie 1 oznacza doskonałą zgodność modelu z danymi obserwowanymi.

Interpretacja: Wartości powyżej 0 wskazują na to, że model lepiej odwzorowuje dane niż średnia obserwacji, natomiast wartości poniżej 0 sugerują, że model jest mniej dokładny niż średnia wartość obserwacji.


2. Współczynnik determinacji (R²)

Mierzy stopień dopasowania modelu do danych obserwowanych, pokazując proporcję zmienności wyjaśnianą przez model.

gdzie:

  • — przepływ obserwowany [m³·s⁻¹]
  • — przepływ obliczony za pomocą modelu [m³·s⁻¹]
  • — wartość średnia z przepływów obserwowanych i modelowanych [m³·s⁻¹]

Zakres: od 0 do 1, gdzie 1 oznacza doskonałe dopasowanie.

Interpretacja: Czym wartość bliższa 1, tym lepsze dopasowanie modelu do danych obserwowanych.


3. Pierwiastek z błędu średniokwadratowego (RMSE)

Mierzy średni błąd modelu, wyrażony w tych samych jednostkach co dane.

gdzie:

  • — przepływ obserwowany [m³·s⁻¹]
  • — przepływ obliczony za pomocą modelu [m³·s⁻¹]
  • — liczba symulacji i obserwacji

Interpretacja: Czym mniejsza wartość RMSE, tym lepsze dopasowanie. Niższa wartość oznacza mniejszy błąd, czyli lepszą dokładność modelu.


4. Całkowity błąd średni (MAE)

Średnia różnica bezwzględna między wartościami obserwowanymi a modelowanymi. MAE jest prostszym wskaźnikiem niż RMSE, ponieważ unika kwadratów błędów.

gdzie:

  • — przepływ obserwowany [m³·s⁻¹]
  • — przepływ obliczony za pomocą modelu [m³·s⁻¹]
  • — liczba symulacji i obserwacji

Interpretacja: Czym mniejsza wartość MAE, tym lepsze dopasowanie. Niższa wartość oznacza mniejszy błąd, czyli lepszą dokładność modelu.


5. Współczynnik Kling-Gupta (KGE)

Kompleksowy współczynnik uwzględniający korelację, bias (odchylenie średnie) i współczynnik zmienności pomiędzy obserwacjami a wynikami modelu.

gdzie:

  • — współczynnik korelacji
  • — stosunek średnich wartości
  • — stosunek odchyleń standardowych

Zakres: od minus nieskończoności do 1, gdzie 1 oznacza doskonałe dopasowanie.

Interpretacja: Wartość bliższa 1 oznacza lepszą jakość modelu.


6. Błąd systematyczny (Bias)

Mierzy średnią różnicę między wartościami obserwowanymi a modelowanymi, pokazując, czy model systematycznie przeszacowuje lub niedoszacowuje wartości.

gdzie:

  • — przepływ obserwowany [m³·s⁻¹]
  • — przepływ obliczony za pomocą modelu [m³·s⁻¹]
  • — liczba symulacji i obserwacji

Interpretacja: Wartość dodatnia oznacza przeszacowanie, a ujemna niedoszacowanie przez model.


7. Procentowy błąd systematyczny (PBIAS)

Mierzy procentowy błąd systematyczny, wyrażając odchylenie modelu w odniesieniu do obserwowanych wartości.

gdzie:

  • — przepływ obserwowany [m³·s⁻¹]
  • — przepływ obliczony za pomocą modelu [m³·s⁻¹]

Interpretacja: Wartości dodatnie wskazują na niedoszacowanie przez model, a ujemne na przeszacowanie.


8. Błąd średniokwadratowy reszt (MSE)

Średnia wartość kwadratów różnic między obserwacjami a predykcjami modelu.

gdzie:

  • — przepływ obserwowany [m³·s⁻¹]
  • — przepływ obliczony za pomocą modelu [m³·s⁻¹]

Interpretacja: Im mniejszy MSE, tym lepsze dopasowanie modelu.


Zastosowanie modeli hydrologicznych

Modele hydrologiczne znajdują szerokie zastosowanie w rozwiązywaniu problemów inżynierskich związanych z wodą.

1. Prognozowanie powodzi (HEC-HMS)

Modele matematyczne opad-odpływ (np. HEC-HMS) są używane do symulacji zachowania rzek podczas intensywnych opadów deszczu. Modele te pomagają przewidywać, kiedy i gdzie mogą wystąpić powodzie oraz jakiej mogą być wielkości.

Zastosowanie:

  • zarządzanie ryzykiem powodziowym
  • projektowanie wałów przeciwpowodziowych
  • planowanie ewakuacji
  • ocena skutków zmian klimatycznych

2. Zarządzanie zasobami wodnymi (SWAT)

SWAT (Soil and Water Assessment Tool) służy do symulowania obiegu wody, jakości wód i transportu osadów w zlewni.

Zastosowanie:

  • zarządzanie zasobami wodnymi w rolnictwie
  • projektowanie zbiorników retencyjnych
  • ocena ilości wód gruntowych
  • optymalizacja zaopatrzenia w wodę dla regionów

3. Modelowanie przepływów rzecznych (HEC-RAS, MIKE 11)

Modele hydrauliczne, takie jak MIKE 11 czy HEC-RAS, symulują przepływy w rzekach, kanałach i innych ciekach wodnych. Modele te uwzględniają zarówno przepływy jednowymiarowe (1D), jak i wielowymiarowe (2D, 3D).

Zastosowanie:

  • projektowanie konstrukcji budownictwa wodnego (tamy, kanały nawadniające, przepusty, mosty)
  • analizy stabilności koryta rzecznego
  • prognozowanie skutków zmian w gospodarce wodnej

4. Ochrona przed erozją i transport osadów (WEPP)

WEPP (Water Erosion Prediction Project) przewiduje procesy erozji gleby oraz transport rumowiska. Modele te symulują wpływ różnych czynników na erozję:

  • opady
  • nachylenie terenu
  • rodzaj gleby
  • pokrycie terenu roślinnością

Zastosowanie:

  • ochrona gleby w rolnictwie
  • projektowanie systemów przeciwerozyjnych
  • ocena skutków działalności gospodarczej na erozję i transport osadów w ciekach

5. Ocena wpływu zmian klimatycznych (VIC)

Modele klimatyczne połączone z modelami hydrologicznymi pozwalają na symulację przyszłych warunków hydrologicznych w odpowiedzi na zmiany klimatyczne. Przykładem mogą być modele typu VIC (Variable Infiltration Capacity) czy PRECIS.

Zastosowanie:

  • prognozowanie zmian w dostępności zasobów wodnych
  • planowanie adaptacji infrastruktury wodnej
  • ocena ryzyka związanego z suszami, powodziami i innymi ekstremalnymi zjawiskami hydrologicznymi

6. Projektowanie systemów odwadniających (SWMM)

Storm Water Management Model (SWMM) symuluje spływ powierzchniowy w obszarach zurbanizowanych, uwzględniając sieci kanalizacyjne i systemy odwadniające.

Zastosowanie:

  • projektowanie kanalizacji
  • zarządzanie wodami opadowymi w miastach
  • minimalizacja ryzyka podtopień w terenach zabudowanych

Podsumowanie

Modele hydrologiczne zlewni stanowią fundamentalne narzędzie w hydrologii inżynierskiej, umożliwiające symulację i prognozowanie procesów hydrologicznych. Prawidłowy wybór typu modelu, dokładna identyfikacja parametrów oraz staranna weryfikacja są kluczowe dla uzyskania wiarygodnych wyników. Stosowanie odpowiednich miar jakości modelu (NSE, R², RMSE, MAE, KGE, Bias, PBIAS, MSE) pozwala obiektywnie ocenić przydatność modelu do zastosowań praktycznych.


Powiązane materiały

Źródła

  1. Byczkowski A. 1996. Hydrologia, tom I. Wyd. SGGW, Warszawa.
  2. Soczyńska U. 1997. Hydrologia dynamiczna. Wyd. PWN, Warszawa.
  3. https://www.hec.usace.army.mil/software/hec-hms/
  4. https://www.researchgate.net/figure/Overview-of-the-swat-model-model-input-output-parameters_fig2_259527294
  5. https://pl.wikipedia.org/wiki/HEC-RAS
  6. https://en.wikipedia.org/wiki/WEPP
  7. https://vic.readthedocs.io/en/master/Overview/ModelOverview/
  8. https://www.appsolutelydigital.com/ModelPrimer/chapter10_section2.html

Zobacz także