Einzugsgebiet: Unterschied zwischen den Versionen

Aus TALSIM Docs
(Diese Seite wurde zum Übersetzen freigegeben)
 
(26 dazwischenliegende Versionen von 3 Benutzern werden nicht angezeigt)
Zeile 1: Zeile 1:
<languages/>
<translate>
<!--T:1-->
{{Navigation|vorher=|hoch=Beschreibung der Systemelemente|nachher=Einleitung}}
[[Datei:Systemelement001.png|50px|none|Symbol Systemelement Einzugsgebiet]]Die Simulation natürlicher Einzugsgebiete verlangt die Bestimmung der Belastungsbildung, Abflussaufteilung und der Abflusskonzentration. Im Folgenden sind die zugrundeliegenden Berechnungsansätze aufgeführt.
[[Datei:Systemelement001.png|50px|none|Symbol Systemelement Einzugsgebiet]]Die Simulation natürlicher Einzugsgebiete verlangt die Bestimmung der Belastungsbildung, Abflussaufteilung und der Abflusskonzentration. Im Folgenden sind die zugrundeliegenden Berechnungsansätze aufgeführt.


==Belastungsbildung==
 
==Belastungsbildung== <!--T:2-->
 
<!--T:3-->
Die Belastungsbildung beschreibt die Ermittlung des Gebietsniederschlags für das betrachtete Einzugsgebiet. Pro Einzugsgebiet wird nur ein Niederschlag benutzt. Liegen mehrere Niederschlagsstationen im Einzugsgebiet vor, so ist es zweckmäßig das Gebiet in mehrere Systemelemente 'Einzugsgebiet' zu unterteilen, bis jedem Element wieder nur ein Niederschlag zugeordnet werden kann.
Die Belastungsbildung beschreibt die Ermittlung des Gebietsniederschlags für das betrachtete Einzugsgebiet. Pro Einzugsgebiet wird nur ein Niederschlag benutzt. Liegen mehrere Niederschlagsstationen im Einzugsgebiet vor, so ist es zweckmäßig das Gebiet in mehrere Systemelemente 'Einzugsgebiet' zu unterteilen, bis jedem Element wieder nur ein Niederschlag zugeordnet werden kann.


==Abflussbildung befestigter/unbefestigter Flächen==
 
==Abflussbildung befestigter/unbefestigter Flächen== <!--T:4-->
 
<!--T:5-->
Die Abflussbildung ermittelt aus dem gefallenen Niederschlag den Effektivniederschlag und daraus abgeleitet die Komponenten Oberflächenabfluss, Infiltration, Verdunstung und Interflow. Eine Schneeberechnung wird durchgeführt bei Temperaturen unter Null °C und erfolgt anhand des Snow-Compaction-Verfahrens. Bezüglich der Algorithmen des Verfahrens wird auf die einschlägige Literatur verwiesen.  
Die Abflussbildung ermittelt aus dem gefallenen Niederschlag den Effektivniederschlag und daraus abgeleitet die Komponenten Oberflächenabfluss, Infiltration, Verdunstung und Interflow. Eine Schneeberechnung wird durchgeführt bei Temperaturen unter Null °C und erfolgt anhand des Snow-Compaction-Verfahrens. Bezüglich der Algorithmen des Verfahrens wird auf die einschlägige Literatur verwiesen.  
Der natürlich ablaufende Prozess vom Niederschlag zum Abfluss wird für die mathematische Simulation in einzelne Phasen untergliedert. In der Abflussbildungsphase wird die Aufteilung des Niederschlages (Systembelastung) in den direkt zum Abfluss gelangenden "wirksamen Niederschlag" und die abflussunwirksamen Verluste (Benetzungs-, Mulden-, Verdunstungs- und Versickerungsverlust) vorgenommen. Dementsprechend wird diese Phase auch mit Belastungsaufteilung bezeichnet. Die resultierende mathematische Gleichung für die momentane Belastungsaufteilung schreibt sich wie folgt:
Der natürlich ablaufende Prozess vom Niederschlag zum Abfluss wird für die mathematische Simulation in einzelne Phasen untergliedert. In der Abflussbildungsphase wird die Aufteilung des Niederschlages (Systembelastung) in den direkt zum Abfluss gelangenden "wirksamen Niederschlag" und die abflussunwirksamen Verluste (Benetzungs-, Mulden-, Verdunstungs- und Versickerungsverlust) vorgenommen. Dementsprechend wird diese Phase auch mit Belastungsaufteilung bezeichnet. Die resultierende mathematische Gleichung für die momentane Belastungsaufteilung schreibt sich wie folgt:


<!--T:6-->
<math>N_W(t) =N(t) -VP(t) -I(t) - \frac{dO}{dt} - \frac{dS}{dt}</math>
<math>N_W(t) =N(t) -VP(t) -I(t) - \frac{dO}{dt} - \frac{dS}{dt}</math>


<!--T:7-->
mit:
mit:
{|style="margin-left: 40px;"
{|style="margin-left: 40px;"
Zeile 26: Zeile 39:
|}
|}


<!--T:8-->
Nachfolgend werden die in der Gleichung verwendeten Terme und deren Berechnung im Einzelnen erläutert.
Nachfolgend werden die in der Gleichung verwendeten Terme und deren Berechnung im Einzelnen erläutert.


===Niederschlag <math>N(t)</math>===
 
===Niederschlag N(t)=== <!--T:9-->
 
<!--T:10-->
Die Niederschlagsdaten müssen dem Simulationsmodell in Form von Regenreihen zur Verfügung gestellt werden. Hierbei ist es prinzipiell unerheblich, ob die Niederschlagsreihe ein Blockregen, ein Modellregen, ein gemessener natürlicher Regen, ein Regenspektrum oder eine langjährige Regenreihe ist. Je nach Zielsetzung der Simulationsrechnung ist die geeignete Belastung ausgewählt werden. Die Regenreihen stammen entweder aus der Zeitreihenverwaltung von Talsim-NG oder werden wie bei Anwendung einer Kurzfristprognose durch die Eingabe einer Regendauer, einer Niederschlagshöhe und der Wahl eines Modellregens direkt vor einer Simulation erzeugt.
Die Niederschlagsdaten müssen dem Simulationsmodell in Form von Regenreihen zur Verfügung gestellt werden. Hierbei ist es prinzipiell unerheblich, ob die Niederschlagsreihe ein Blockregen, ein Modellregen, ein gemessener natürlicher Regen, ein Regenspektrum oder eine langjährige Regenreihe ist. Je nach Zielsetzung der Simulationsrechnung ist die geeignete Belastung ausgewählt werden. Die Regenreihen stammen entweder aus der Zeitreihenverwaltung von Talsim-NG oder werden wie bei Anwendung einer Kurzfristprognose durch die Eingabe einer Regendauer, einer Niederschlagshöhe und der Wahl eines Modellregens direkt vor einer Simulation erzeugt.


===Verdunstung <math>VP(t)</math>===
[[Datei:Jahresgang_ETpot.png|thumb|Jahresgang der potentiellen Verdunstung nach /BRANDT, 1979/]]
[[Datei:Tagesgang_ETpot.png|thumb|Tagesgang der potentiellen Verdunstung als Vielfaches der mittleren Tagesverdunstung]]Die Verdunstung wirkt sich in zweifacher Weise auf die Abflussbildung aus. Zum einen sind die Anfangsbedingungen im Einzugsgebiet (Benetzung und Muldenfüllung auf der Oberfläche sowie eingeschränkt auch die Bodenfeuchte bei durchlässigen Flächen) ein Resultat des vor dem betrachteten Niederschlagsereignis stattfindenden Verdunstungsgeschehens, zum anderen wird der zu berechnende abflusswirksame Niederschlag um den Betrag der momentanen Verdunstungsrate geschmälert.
Die potentielle (energetisch mögliche) Verdunstung VP ist zeitlich und örtlich sehr unterschiedlich und einer genauen Berechnung nur schwer zugänglich. Aus ausgewerteten Messungen von 20 Stationen, deren Mittelwerte als Histogramm in dargestellt sind, wurde folgende Ausgleichsfunktion  ermittelt /BRANDT, 1979/.


===Verdunstung VP(t)=== <!--T:11-->
<!--T:12-->
Die Verdunstung wirkt sich in zweifacher Weise auf die Abflussbildung aus. Zum einen sind die Anfangsbedingungen im Einzugsgebiet (Benetzung und Muldenfüllung auf der Oberfläche sowie eingeschränkt auch die Bodenfeuchte bei durchlässigen Flächen) ein Resultat des vor dem betrachteten Niederschlagsereignis stattfindenden Verdunstungsgeschehens, zum anderen wird der zu berechnende abflusswirksame Niederschlag um den Betrag der momentanen Verdunstungsrate geschmälert.
<!--T:57-->
Die potentielle (energetisch mögliche) Verdunstung VP ist zeitlich und örtlich sehr unterschiedlich und einer genauen Berechnung nur schwer zugänglich.
<!--T:58-->
Es stehen verschiedene Optionen zur Eingabe/Berechnung der potentiellen Verdunstung zur Verfügung:
* Vorgabe einer externen Zeitreihe mit Werten für die potentielle Verdunstung
* Interne Berechnung der potentiellen Verdunstung: Die potentielle Verdunstung wird auf Basis von Temperatur und je nach Verfahren weiteren Parametern intern berechnet. Es stehen die folgenden Verfahren zur Verfügung:
** Penman: Erfordert die Angabe mittlerer Jahresgänge oder von Zeitreihen für Sonnenscheindauer, Windgeschwindigkeit und relative Luftfeuchte
** Haude: Erfordert die Angabe eines mittleren Jahresgangs oder einer Zeitreihe relative Luftfeuchte
** Turc: Erfordert die Angabe mittlerer Jahresgänge oder von Zeitreihen für Sonnenscheindauer und relative Luftfeuchte
** Blaney-Criddle: Vereinfachtes Verfahren, erfordert lediglich die Angabe des Breitengrads, woraus die Sonnenscheindauer abgeleitet wird.
* Angabe einer fixen Verdungstungshöhe pro Jahr: Bei dieser Option wird die eingegebene Jahressumme mit der Ausgleichsfunktion nach Brandt (siehe unten) innerjährlich verteilt.
====Ausgleichsfunktionen==== <!--T:59-->
Die folgenden Ausgleichsfunktionen nach Brandt werden verwendet, um eingegebene Jahressummen der Verdunstung über das Jahr zu verteilen bzw. um Tageswerte auf Stundenwerte umzurechnen.
<!--T:60-->
Aus ausgewerteten Messungen von 20 Stationen, deren Mittelwerte als Histogramm in dargestellt sind, wurde folgende Ausgleichsfunktion  ermittelt /BRANDT, 1979/.
<!--T:13-->
<math>VP=(0.96+0.0033 \cdot i) \cdot \sin\frac{2\pi}{365}(i-148)+158</math>
<math>VP=(0.96+0.0033 \cdot i) \cdot \sin\frac{2\pi}{365}(i-148)+158</math>


<!--T:14-->
mit:
mit:
{|style="margin-left: 40px;"
{|style="margin-left: 40px;"
Zeile 46: Zeile 85:
|}
|}


Die jährliche potentielle Gesamtverdunstungshöhe beträgt 642 mm. Liegen keine gemessenen Verdunstungswerte vor, kann optional dieser normierte Jahresgang der potentiellen Verdunstung für die Berechnung der aktuellen Verdunstung herangezogen werden. Ist das gewählte Berechnungszeitintervall kleiner als ein Tag, wird mittels des dargestellten Tagesgangs letztendlich die potentielle Verdunstung für jedes Berechnungszeitintervall ermittelt. Ist das Berechnungsintervall &le; 1 Tag entfällt die Berücksichtigung des Tagesganges.
<!--T:15-->
Die jährliche potentielle Gesamtverdunstungshöhe beträgt 642 mm. Liegen keine gemessenen Verdunstungswerte vor, kann optional dieser normierte Jahresgang der potentiellen Verdunstung für die Berechnung der aktuellen Verdunstung herangezogen werden. Ist das gewählte Berechnungszeitintervall kleiner als ein Tag, wird mittels des dargestellten Tagesgangs letztendlich die potentielle Verdunstung für jedes Berechnungszeitintervall ermittelt. Ist das Berechnungsintervall &ge; 1 Tag entfällt die Berücksichtigung des Tagesganges.
<gallery mode="packed" heights=300px">
Datei:Jahresgang_ETpot.png|Jahresgang der potentiellen Verdunstung nach /BRANDT, 1979/
Datei:Tagesgang_ETpot.png|Tagesgang der potentiellen Verdunstung als Vielfaches der mittleren Tagesverdunstung
</gallery>


===Oberflächenwasservorrat (versiegelter Flächenanteil) <math>O</math>===
===Oberflächenwasservorrat (versiegelter Flächenanteil) O=== <!--T:16-->


<!--T:17-->
Bei den versiegelten Flächenanteilen kann neben dem Schneevorrat auch die Infiltration vernachlässigt werden, so dass sich die Bilanzgleichung wie folgt vereinfacht:
Bei den versiegelten Flächenanteilen kann neben dem Schneevorrat auch die Infiltration vernachlässigt werden, so dass sich die Bilanzgleichung wie folgt vereinfacht:


<!--T:18-->
<math>N_W(t)=N(t)-VP(t)-\frac{dO}{dt}</math>
<math>N_W(t)=N(t)-VP(t)-\frac{dO}{dt}</math>


<!--T:19-->
wobei die Oberflächenwasservorratsänderung <math>dO/dt</math> die Benetzung der Oberfläche sowie die Auffüllung und Entleerung (durch Verdunstung) der Mulden repräsentiert.[[Datei:Schema_Modellansatz_Benetzungs-_und_Muldenverlust.png|thumb|Schema der Modellansätze Benetzungs- und Muldenverlust]]
wobei die Oberflächenwasservorratsänderung <math>dO/dt</math> die Benetzung der Oberfläche sowie die Auffüllung und Entleerung (durch Verdunstung) der Mulden repräsentiert.[[Datei:Schema_Modellansatz_Benetzungs-_und_Muldenverlust.png|thumb|Schema der Modellansätze Benetzungs- und Muldenverlust]]
Als Benetzungsverlust <math>BV</math> für versiegelte Flächen wird folgender Standardwert angesetzt.
Als Benetzungsverlust <math>BV</math> für versiegelte Flächen wird folgender Standardwert angesetzt.


<!--T:20-->
<math>BV = 0.5 \mbox{ mm}</math>
<math>BV = 0.5 \mbox{ mm}</math>


<!--T:21-->
Der Muldenverlust MV wird durch den Anwender vorgegeben. Der Standard- und gleichzeitig Maximalwert im Modell beträgt 4 mm  
Der Muldenverlust MV wird durch den Anwender vorgegeben. Der Standard- und gleichzeitig Maximalwert im Modell beträgt 4 mm  
Der Muldenverlust stellt den Mittelwert für eine geneigte Oberfläche dar. Da die Mulden jedoch nicht gleichmäßig verteilt sind und erfahrungsgemäß bereits ein Abfluss einsetzt, bevor überall die komplette Muldenauffüllung erreicht ist, wird unterstellt, dass jeweils
Der Muldenverlust stellt den Mittelwert für eine geneigte Oberfläche dar. Da die Mulden jedoch nicht gleichmäßig verteilt sind und erfahrungsgemäß bereits ein Abfluss einsetzt, bevor überall die komplette Muldenauffüllung erreicht ist, wird unterstellt, dass jeweils
Zeile 66: Zeile 115:
aufweist. Somit kommt es bereits zum Abfluss, wenn der um die Verdunstungsrate verminderte Niederschlag den Benetzungsverlust und 1/3 des Muldenverlustes übersteigt (bei trockener Vorgeschichte). Die o.g. Annahmen sind in der folgenden Abbildung schematisch skizziert.
aufweist. Somit kommt es bereits zum Abfluss, wenn der um die Verdunstungsrate verminderte Niederschlag den Benetzungsverlust und 1/3 des Muldenverlustes übersteigt (bei trockener Vorgeschichte). Die o.g. Annahmen sind in der folgenden Abbildung schematisch skizziert.


<!--T:22-->
Der Abflussbeiwert der versiegelten Flächen (nach Abdeckung der Anfangsverluste) wird mit <math>\Psi = 1</math> angesetzt. Bei der Festlegung des versiegelten Flächenanteils in einem Teileinzugsgebiet ist zu beachten, dass nicht alle befestigten oder versiegelten Flächen tatsächlich in eine Kanalisation entwässern.
Der Abflussbeiwert der versiegelten Flächen (nach Abdeckung der Anfangsverluste) wird mit <math>\Psi = 1</math> angesetzt. Bei der Festlegung des versiegelten Flächenanteils in einem Teileinzugsgebiet ist zu beachten, dass nicht alle befestigten oder versiegelten Flächen tatsächlich in eine Kanalisation entwässern.
Die kontinuierliche Bereitstellung der Benetzungs- und Muldenverluste erfolgt über die laufende Bilanzierung dieser Speicher und der Verdunstung.
Die kontinuierliche Bereitstellung der Benetzungs- und Muldenverluste erfolgt über die laufende Bilanzierung dieser Speicher und der Verdunstung.


===Oberflächenwasservorrat (unversiegelter Flächenanteil) <math>O</math>===
 
Der Oberflächenwasservorrat wird über die Bilanzierung eines Verlustspeichers in Abhängigkeit des gewählten Abflussbildungsansatzes berechnet. Einzelheiten dazu finden sich in den folgenden Abschnitten zur Berechnung der Infiltration bzw. abflusswirksamer Niederschlag.  
===Oberflächenwasservorrat (unversiegelter Flächenanteil) O=== <!--T:23-->
===Infiltration bzw. abflusswirksamer Niederschlag <math>I(t)</math>, <math>N_W(t)</math>===
 
<!--T:24-->
Der Oberflächenwasservorrat wird über die Bilanzierung eines Verlustspeichers in Abhängigkeit des gewählten Abflussbildungsansatzes berechnet. Einzelheiten dazu finden sich in den folgenden Abschnitten zur Berechnung der Infiltration bzw. abflusswirksamer Niederschlag.
 
 
===Infiltration bzw. abflusswirksamer Niederschlag I(t), N<sub>W</sub>(t)=== <!--T:25-->
 
<!--T:26-->
Bei den durchlässigen Flächen kann die Infiltration in den Boden nicht vernachlässigt werden, da diese das Abflussgeschehen entscheidend prägt. Für die Berechnung wurden drei Ansätze im Modell implementiert:
Bei den durchlässigen Flächen kann die Infiltration in den Boden nicht vernachlässigt werden, da diese das Abflussgeschehen entscheidend prägt. Für die Berechnung wurden drei Ansätze im Modell implementiert:
# Konstanter Abflussbeiwert <math>\Psi</math>
# Konstanter Abflussbeiwert <math>\Psi</math>
Zeile 77: Zeile 134:
# Bodenfeuchtesimulation
# Bodenfeuchtesimulation


====Konstanter Abflussbeiwert <math>\Psi</math>====


====Konstanter Abflussbeiwert &Psi;==== <!--T:27-->
<!--T:28-->
Bei Angabe eines <math>\Psi</math>-Wertes kommt nach Abdeckung der Anfangsverluste (Benetzungs- und Muldenverlust) der übrige Anteil des Niederschlages im Verhältnis des Abflussbeiwertes <math>\Psi</math> zum Abfluss und zwar unabhängig von der Vorgeschichte und den Merkmalen des Niederschlages (Höhe, Intensität, Dauer). Auf diesen Ansatz sollte nach Möglichkeit verzichtet werden, da hier der Prozess der Abflussbildung nur grob vereinfachend beschrieben wird.
Bei Angabe eines <math>\Psi</math>-Wertes kommt nach Abdeckung der Anfangsverluste (Benetzungs- und Muldenverlust) der übrige Anteil des Niederschlages im Verhältnis des Abflussbeiwertes <math>\Psi</math> zum Abfluss und zwar unabhängig von der Vorgeschichte und den Merkmalen des Niederschlages (Höhe, Intensität, Dauer). Auf diesen Ansatz sollte nach Möglichkeit verzichtet werden, da hier der Prozess der Abflussbildung nur grob vereinfachend beschrieben wird.


====Ereignisspezifischer Abflussbeiwert in Anlehnung an das Verfahren des Soil-Conservation-Service (SCS)====
 
====Ereignisspezifischer Abflussbeiwert in Anlehnung an das Verfahren des Soil-Conservation-Service (SCS)==== <!--T:29-->
 
<!--T:30-->
Bei Angabe eines vom Bodentyp und der Bodennutzung abhängigen CN-Wertes (siehe /DVWK, 1991/) lässt sich ein vorgeschichtsabhängiger Anfangsverlust sowie eine vorgeschichtsabhängige Beziehung des Abflussbeiwertes von der bis zum betrachteten Zeitpunkt akkumulierten Niederschlagshöhe formulieren /Zaiss, 1987/; d.h. der Abflussbeiwert wächst mit zunehmendem Niederschlag im Verlauf des Ereignisses an.
Bei Angabe eines vom Bodentyp und der Bodennutzung abhängigen CN-Wertes (siehe /DVWK, 1991/) lässt sich ein vorgeschichtsabhängiger Anfangsverlust sowie eine vorgeschichtsabhängige Beziehung des Abflussbeiwertes von der bis zum betrachteten Zeitpunkt akkumulierten Niederschlagshöhe formulieren /Zaiss, 1987/; d.h. der Abflussbeiwert wächst mit zunehmendem Niederschlag im Verlauf des Ereignisses an.
Die Quantifizierung der Vorgeschichte erfolgt über den 21-Tage-Vorregenindex <math>VN</math>
Die Quantifizierung der Vorgeschichte erfolgt über den 21-Tage-Vorregenindex <math>VN</math>
[[Datei:Abhängigkeit_des_Abflussbeiwertes_von_der_Vorgeschichte.png|thumb|Abhängigkeit des Abflussbeiwertes von der Vorgeschichte]][[Datei:Abhängigkeit_des_Abflussbeiwertes_von_der_kumulierten_Niederschlagssumme.png|thumb|Abhängigkeit des Abflussbeiwertes von der kumulierten Niederschlagssumme]]
 
<!--T:31-->
<math>V_N=\sum_{j=1}^21 C(j)^j \cdot hN_j</math>
<math>V_N=\sum_{j=1}^21 C(j)^j \cdot hN_j</math>


<!--T:32-->
mit:
mit:
{|style="margin-left: 40px;"
{|style="margin-left: 40px;"
Zeile 95: Zeile 159:
|}
|}


<!--T:33-->
Der Einfluss der Jahreszeit wird durch einen Jahresgang des Faktors C wiedergegeben.
Der Einfluss der Jahreszeit wird durch einen Jahresgang des Faktors C wiedergegeben.


<!--T:34-->
<math>C=0.05 \cdot \sin\frac{2\pi}{365}(i+0.75)+0.85</math>
<math>C=0.05 \cdot \sin\frac{2\pi}{365}(i+0.75)+0.85</math>


<!--T:35-->
mit:
mit:
{|style="margin-left: 40px;"
{|style="margin-left: 40px;"
Zeile 105: Zeile 172:
|}
|}


<!--T:36-->
Damit schwankt der Wert C zwischen 0.8 < C < 0.9. Hierdurch wird erreicht, dass bei gleichem Vorregen zu unterschiedlichen Jahreszeiten unterschiedliche Vorregenindizes berechnet und damit eine veränderte Abflussbereitschaft in Rechnung gestellt wird.
Damit schwankt der Wert C zwischen 0.8 < C < 0.9. Hierdurch wird erreicht, dass bei gleichem Vorregen zu unterschiedlichen Jahreszeiten unterschiedliche Vorregenindizes berechnet und damit eine veränderte Abflussbereitschaft in Rechnung gestellt wird.
In Abhängigkeit von der auf diese Weise quantifizierten Vorgeschichte kann unter Verwendung der gebietsspezifischen und für mittlere Vorfeuchteverhältnisse gültigen CN-Werte ein aktueller Abflussbeiwert berechnet werden. In der folgenden Abbildung ist für unterschiedliche CN-Werte dargestellt, wie sich der aktuelle Abflussbeiwert in Abhängigkeit von der Vorgeschichte verändert.
In Abhängigkeit von der auf diese Weise quantifizierten Vorgeschichte kann unter Verwendung der gebietsspezifischen und für mittlere Vorfeuchteverhältnisse gültigen CN-Werte ein aktueller Abflussbeiwert berechnet werden. In der folgenden Abbildung ist für unterschiedliche CN-Werte dargestellt, wie sich der aktuelle Abflussbeiwert in Abhängigkeit von der Vorgeschichte verändert.
Da sich im Verlaufe eines Regenereignisses durch die Durchfeuchtung des Bodens die Abflussbereitschaft eines Einzugsgebiets verändert, wird ebenfalls eine Anpassung des Abflussbeiwertes während eines Ereignisses als Funktion der kumulierten Niederschlagshöhe vorgenommen.
<gallery mode="packed" heights=300px>
Datei:Abhängigkeit_des_Abflussbeiwertes_von_der_Vorgeschichte.png|Abhängigkeit des Abflussbeiwertes von der Vorgeschichte
Datei:Abhängigkeit_des_Abflussbeiwertes_von_der_kumulierten_Niederschlagssumme.png|Abhängigkeit des Abflussbeiwertes von der kumulierten Niederschlagssumme
</gallery>
====Bodenfeuchtesimulation==== <!--T:37-->


Da sich im Verlaufe eines Regenereignisses durch die Durchfeuchtung des Bodens die Abflussbereitschaft eines Einzugsgebiets verändert, wird ebenfalls eine Anpassung des Abflussbeiwertes während eines Ereignisses als Funktion der kumulierten Niederschlagshöhe vorgenommen.


====Bodenfeuchtesimulation====
=====Landnutzung===== <!--T:38-->
=====Landnutzung=====


<!--T:39-->
Bei der Anwendung der Bodenfeuchtesimulation ist die Angabe von Landnutzungen notwendig. Aus den Angaben zur Landnutzung wird die Durchwurzelungstiefe benötigt, um die Dicke der Durchwurzelungsschicht zu ermitteln. Weitere Parameter der Landnutzung, die zur Berechnung der Interzeption und der Transpiration dienen, sind:
Bei der Anwendung der Bodenfeuchtesimulation ist die Angabe von Landnutzungen notwendig. Aus den Angaben zur Landnutzung wird die Durchwurzelungstiefe benötigt, um die Dicke der Durchwurzelungsschicht zu ermitteln. Weitere Parameter der Landnutzung, die zur Berechnung der Interzeption und der Transpiration dienen, sind:
* Wurzeltiefe
* Wurzeltiefe
Zeile 121: Zeile 196:
Die Angabe von Haude-Faktoren zur besseren Berücksichtigung der Verdunstung je Landnutzung sind über Eingabe von Jahresgängen beliebig möglich und können den gewünschten Landnutzungen zugeordnet werden.
Die Angabe von Haude-Faktoren zur besseren Berücksichtigung der Verdunstung je Landnutzung sind über Eingabe von Jahresgängen beliebig möglich und können den gewünschten Landnutzungen zugeordnet werden.


=====Bodentyp/ Bodenart=====
 
[[Datei:Programminterne_Zusammenfassung_Bodenschichten.png|thumb|Beispiel der Zusammenfassung der Bodenschichten zu einer programminternen Schicht anhand der Durchwurzelungsschicht]][[Datei:Schema_Bodenfeuchtesimulation.png|thumb|Mit der Bodenfeuchtesimulation berechnete Größen]]
=====Bodentyp/ Bodenart===== <!--T:40-->
 
<!--T:41-->
[[Datei:Programminterne_Zusammenfassung_Bodenschichten.png|thumb|Beispiel der Zusammenfassung der Bodenschichten zu einer programminternen Schicht anhand der Durchwurzelungsschicht]]
[[Datei:Schema_Bodenfeuchtesimulation.png|thumb|Mit der Bodenfeuchtesimulation berechnete Größen]]
Die Bodenfeuchtesimulation basiert auf einer nichtlinearen Berechnung der einzelnen Bodenhorizonte. Der Boden wird dabei in verschiedene Horizonte (Schichten) eingeteilt. Jede Schicht wird berechnet und mit den (falls vorhanden) darunter bzw. darüber liegenden Schichten abgeglichen. Als Parameter zur Bodenfeuchteberechnung dienen folgende bodenphysikalischen Größen:
Die Bodenfeuchtesimulation basiert auf einer nichtlinearen Berechnung der einzelnen Bodenhorizonte. Der Boden wird dabei in verschiedene Horizonte (Schichten) eingeteilt. Jede Schicht wird berechnet und mit den (falls vorhanden) darunter bzw. darüber liegenden Schichten abgeglichen. Als Parameter zur Bodenfeuchteberechnung dienen folgende bodenphysikalischen Größen:
* Welkepunkt (WP)
* Welkepunkt (WP)
Zeile 132: Zeile 211:
* Zuordnung zu einer Bodenart: Sand, Schluff, Ton
* Zuordnung zu einer Bodenart: Sand, Schluff, Ton


<!--T:42-->
Die mögliche Anzahl der Bodenschichten läuft von minimal einer bis maximal sechs. Die Erfahrung zeigte, dass die besten Ergebnisse mit einer Aufteilung in drei Schichten erzielt werden konnten. Aus diesem Grund werden die eingegebenen Schichten programmintern immer in drei Horizonte unterteilt.  
Die mögliche Anzahl der Bodenschichten läuft von minimal einer bis maximal sechs. Die Erfahrung zeigte, dass die besten Ergebnisse mit einer Aufteilung in drei Schichten erzielt werden konnten. Aus diesem Grund werden die eingegebenen Schichten programmintern immer in drei Horizonte unterteilt.  
* Infiltrationsschicht (Standarddicke [cm] = 20)
* Infiltrationsschicht (Standarddicke [cm] = 20)
Zeile 137: Zeile 217:
* Transportschicht (Mindestdicke [cm] = 5)
* Transportschicht (Mindestdicke [cm] = 5)


<!--T:43-->
Die Berechnung der neuen Bodenkennwerte für die programmintern verwendeten Schichten erfolgt durch eine Gewichtung entsprechend den vorgegebenen original Dicken der Schichten. Im Fall der gesättigten Leitfähigkeit läuft die Berechnung nach dem Prinzip der Erhaltung der Kontinuität der Strömung ab. Bei senkrechter Strömung soll aufgrund der Kontinuität der Strömung die Geschwindigkeit v bei gegebener Durchflussmenge in einer programminternen Schicht denselben Wert besitzen. Damit ist das hydraulische Gefälle nicht mehr konstant.
Die Berechnung der neuen Bodenkennwerte für die programmintern verwendeten Schichten erfolgt durch eine Gewichtung entsprechend den vorgegebenen original Dicken der Schichten. Im Fall der gesättigten Leitfähigkeit läuft die Berechnung nach dem Prinzip der Erhaltung der Kontinuität der Strömung ab. Bei senkrechter Strömung soll aufgrund der Kontinuität der Strömung die Geschwindigkeit v bei gegebener Durchflussmenge in einer programminternen Schicht denselben Wert besitzen. Damit ist das hydraulische Gefälle nicht mehr konstant.


<!--T:44-->
<math>kf_V=\frac{\sum_{i=1}^n d_i}{\sum_{i=1}^n \frac{d_i}{k_i}}</math>
<math>kf_V=\frac{\sum_{i=1}^n d_i}{\sum_{i=1}^n \frac{d_i}{k_i}}</math>


<!--T:45-->
mit:
mit:
{|style="margin-left: 40px;"
{|style="margin-left: 40px;"
Zeile 154: Zeile 237:
Die zu lösende Gleichung ist:
Die zu lösende Gleichung ist:


 
<!--T:46-->
<math>\frac{dBF(t)}{dt}=Inf(t)-Perk(t)-Eva_{akt}(t)-Trans_{akt}(t)-Int(t)+Kap(t)</math>
<math>\frac{dBF(t)}{dt}=Inf(t)-Perk(t)-Eva_{akt}(t)-Trans_{akt}(t)-Int(t)+Kap(t)</math>


<!--T:47-->
mit:
mit:
{|style="margin-left: 40px;"
{|style="margin-left: 40px;"
Zeile 174: Zeile 258:
|-
|-
|}
|}
 
[[Datei:Bodenprozessfunktionen.png|thumb|Darstellung ausgewählter Bodenprozessfunktionen]]
Infiltration, Perkolation, Evaporation, Transpiration, Interflow und Kapillaraufstieg sind von der aktuellen Bodenfeuchte abhängig. In der Simulation wird diese Abhängigkeit durch folgende Funktionsverläufe beschrieben.
Infiltration, Perkolation, Evaporation, Transpiration, Interflow und Kapillaraufstieg sind von der aktuellen Bodenfeuchte abhängig. In der Simulation wird diese Abhängigkeit durch folgende Funktionsverläufe beschrieben.


<math>Inf(BF(t))=a_v \cdot \left(GPV-BF(t) \right)^{1.4}+k_f</math>
<!--T:48-->
<math>Inf(BF(t))=a_v \cdot \left(GPV-BF(t) \right)^{1.4}+k_f </math> &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;(Ansatz nach HOLTAN)


<!--T:49-->
<math>
<math>
Perk(BF(t)) =  
Perk(BF(t)) =  
\begin{cases}  
\begin{cases}  
0,  & BF(t)\le f_{PK} \cdot nFK + WP \\
0,  & BF(t)\le f_{PK} \cdot nFK + WP \\
k_f \cdot \left(\frac{BF(t)-(f_PK \cdot nFK +WP)}{GPV-(f_PK \cdot nFK +WP)}  \right)^exp,PK, & BF(t)> f_{PK} \cdot nFK + WP  
k_f \cdot \left(\frac{BF(t)-(f_{PK} \cdot nFK +WP)}{GPV-(f_{PK} \cdot nFK +WP)}  \right)^{exp,PK}, & BF(t)> f_{PK} \cdot nFK + WP
\end{cases}
</math>
:::::::::::::::::(mod. Ansatz nach /OSTROWSKI, 1992/)
<math>
Eva(BF(t)) =
\begin{cases}
0,  & BF(t)\le WP \\
f_{Eva} \cdot \left(\frac{BF(t)-WP}{GPV-WP}  \right)^{exp,PK}, & BF(t)> WP  
\end{cases}
\end{cases}
</math>
</math>


<math></math>
<!--T:50-->
 
<math>
<math></math>
Trans(BF(t)) =
\begin{cases}
0,  & BF(t)\le f_{Trans} \cdot nFK + WP \\
f_{Trans} \cdot \left(\frac{BF(t)-(f_{Trans} \cdot nFK +WP)}{GPV-(f_{Trans} \cdot nFK +WP)}  \right)^{exp,PK}, & BF(t)> f_{Trans} \cdot nFK + WP
\end{cases}
</math>


<!--T:51-->
mit:
mit:
{|style="margin-left: 40px;"
{|style="margin-left: 40px;"
|<math>d_i</math>: || anteilige Schichtdicke der jeweiligen Original-Schicht [mm]
|<math>a_v</math>: || Infiltrationsfaktor nach HOLTAN (in Talsim-NG <math>a_v=1</math>
|-
|<math>k_f</math>: || Durchlässigkeitsbeiwert des gesättigten Bodens
|-
|<math>nFK</math>: || nutzbare Feldkapazität (<math>nFK=FK-WP</math>)
|-
|<math>WP</math>: || Welkepunkt
|-
|<math>FK</math>: || Feldkapazität
|-
|<math>GPV</math>: || gesamtes Porenvolumen
|-
|<math>f_{PK}</math>: || bodenabhängiger Skalierungfaktor der Perkolationsfunktion
|-
|<math>exp,PK</math>: || bodenabhängiger Krümmungsparameter der Perkolationsfunktion
|-
|<math>f_{Eva}</math>: || bodenabhängiger Skalierungsfaktor der Evaporationsfunktion
|-
|-
|<math>k_i</math>: || gesättigte Leitfähigkeit der jeweiligen Original-Schicht [mm/h]
|<math>f_{Trans}</math>: || bodenabhängiger Skalierungsfaktor der Transpirationsfunktion
|-
|-
|<math>kf_V</math>: || gesättigte Leitfähigkeit der programmintern verwendeten Schicht [mm/h]
|<math>exp,Trans</math>: || Krümmungsparameter der Transpirationsfunktion
|-
|-
|}
|}


[[FORMELN]]
<!--T:52-->
Die Programmparameter werden intern berechnet. Der Anwender muss lediglich die Bodenkennwerte kf, WP, FK und GPV angeben.
Die Berechnung der Bodenprozesse erfolgt mit einem neuentwickelten [[Special:MyLanguage/Berechnungsschema von Speichern|Baustein zur Simulation von Speichern]].
 
 
=====Elementarflächen===== <!--T:53-->


=====Elementarflächen=====
<!--T:54-->
[[Datei:Aufteilung_EZG_in_Elementarflächen.png|thumb|Aufteilung eines Einzugsgebietselementes in Elementarflächen]][[Datei:Elementarflächen_Zuordnung_Bodentyp_Landnutzung.png|thumb|Zuordnung von Bodentyp und Landnutzung zu Elementarflächen]]
Wird mit der Bodenfeuchtesimulation die Abflussbildung berechnet, wird gleichzeitig das Elementarflächenkonzept angewandt. Ein Einzugsgebietselement wird dabei in beliebig viele hydrologisch homogene Flächen unterteilt, d.h. Flächen gleichen Bodentyps und gleicher Landnutzung.
Wird mit der Bodenfeuchtesimulation die Abflussbildung berechnet, wird gleichzeitig das Elementarflächenkonzept angewandt. Ein Einzugsgebietselement wird dabei in beliebig viele hydrologisch homogene Flächen unterteilt, d.h. Flächen gleichen Bodentyps und gleicher Landnutzung.
Für jede Elementarfläche gilt genau eine Zuordnung von Landnutzung und Bodentyp. Die aus einer Elementarfläche resultierende Wassermenge wird am Elementausgang angesetzt, d.h. alle Elementarflächen geben unabhängig ihrer Lage im Einzugsgebiet Wasser mit der gleichen zeitlichen Verzögerung ab.
Für jede Elementarfläche gilt genau eine Zuordnung von Landnutzung und Bodentyp. Die aus einer Elementarfläche resultierende Wassermenge wird am Elementausgang angesetzt, d.h. alle Elementarflächen geben unabhängig ihrer Lage im Einzugsgebiet Wasser mit der gleichen zeitlichen Verzögerung ab.
<gallery mode="packed" heights=200px>
Datei:Aufteilung_EZG_in_Elementarflächen.png|Aufteilung eines Einzugsgebietselementes in Elementarflächen
Datei:Elementarflächen_Zuordnung_Bodentyp_Landnutzung.png|Zuordnung von Bodentyp und Landnutzung zu Elementarflächen
</gallery>
==Abflusskonzentration== <!--T:55-->


==Abflusskonzentration==
<!--T:56-->
Die Abflusskonzentration bestimmt die Verzögerung des Oberflächenabflusses aus dem Einzugsgebiet. Es wird eine Parallelspeicherkaskade mit drei Speichern für unbefestigte und eine Kaskade für befestigte Flächen benutzt. Der Abfluss der Komponenten Interflow und Grundwasser wird über einen linearen Einzelspeicher verzögert an den Elementausgang abgegeben.
Die Abflusskonzentration bestimmt die Verzögerung des Oberflächenabflusses aus dem Einzugsgebiet. Es wird eine Parallelspeicherkaskade mit drei Speichern für unbefestigte und eine Kaskade für befestigte Flächen benutzt. Der Abfluss der Komponenten Interflow und Grundwasser wird über einen linearen Einzelspeicher verzögert an den Elementausgang abgegeben.
[[Datei:Abflusskonzentration.png|frame|none|Berechnung der Abflusskonzentration von Einzugsgebieten]]
[[Datei:Abflusskonzentration.png|frame|none|Berechnung der Abflusskonzentration von Einzugsgebieten]]
</translate>

Aktuelle Version vom 17. Dezember 2020, 10:32 Uhr

Sprachen:
Symbol Systemelement Einzugsgebiet

Die Simulation natürlicher Einzugsgebiete verlangt die Bestimmung der Belastungsbildung, Abflussaufteilung und der Abflusskonzentration. Im Folgenden sind die zugrundeliegenden Berechnungsansätze aufgeführt.


Belastungsbildung

Die Belastungsbildung beschreibt die Ermittlung des Gebietsniederschlags für das betrachtete Einzugsgebiet. Pro Einzugsgebiet wird nur ein Niederschlag benutzt. Liegen mehrere Niederschlagsstationen im Einzugsgebiet vor, so ist es zweckmäßig das Gebiet in mehrere Systemelemente 'Einzugsgebiet' zu unterteilen, bis jedem Element wieder nur ein Niederschlag zugeordnet werden kann.


Abflussbildung befestigter/unbefestigter Flächen

Die Abflussbildung ermittelt aus dem gefallenen Niederschlag den Effektivniederschlag und daraus abgeleitet die Komponenten Oberflächenabfluss, Infiltration, Verdunstung und Interflow. Eine Schneeberechnung wird durchgeführt bei Temperaturen unter Null °C und erfolgt anhand des Snow-Compaction-Verfahrens. Bezüglich der Algorithmen des Verfahrens wird auf die einschlägige Literatur verwiesen. Der natürlich ablaufende Prozess vom Niederschlag zum Abfluss wird für die mathematische Simulation in einzelne Phasen untergliedert. In der Abflussbildungsphase wird die Aufteilung des Niederschlages (Systembelastung) in den direkt zum Abfluss gelangenden "wirksamen Niederschlag" und die abflussunwirksamen Verluste (Benetzungs-, Mulden-, Verdunstungs- und Versickerungsverlust) vorgenommen. Dementsprechend wird diese Phase auch mit Belastungsaufteilung bezeichnet. Die resultierende mathematische Gleichung für die momentane Belastungsaufteilung schreibt sich wie folgt:

[math]\displaystyle{ N_W(t) =N(t) -VP(t) -I(t) - \frac{dO}{dt} - \frac{dS}{dt} }[/math]

mit:

[math]\displaystyle{ N_W }[/math]: abflusswirksamer Niederschlag
[math]\displaystyle{ N }[/math]: Niederschlag
[math]\displaystyle{ VP }[/math]: potentielle Verdunstung
[math]\displaystyle{ I }[/math]: Infiltration in den Bodenraum
[math]\displaystyle{ O }[/math]: Oberflächenwasservorrat
[math]\displaystyle{ S }[/math]: Schneevorrat

Nachfolgend werden die in der Gleichung verwendeten Terme und deren Berechnung im Einzelnen erläutert.


Niederschlag N(t)

Die Niederschlagsdaten müssen dem Simulationsmodell in Form von Regenreihen zur Verfügung gestellt werden. Hierbei ist es prinzipiell unerheblich, ob die Niederschlagsreihe ein Blockregen, ein Modellregen, ein gemessener natürlicher Regen, ein Regenspektrum oder eine langjährige Regenreihe ist. Je nach Zielsetzung der Simulationsrechnung ist die geeignete Belastung ausgewählt werden. Die Regenreihen stammen entweder aus der Zeitreihenverwaltung von Talsim-NG oder werden wie bei Anwendung einer Kurzfristprognose durch die Eingabe einer Regendauer, einer Niederschlagshöhe und der Wahl eines Modellregens direkt vor einer Simulation erzeugt.


Verdunstung VP(t)

Die Verdunstung wirkt sich in zweifacher Weise auf die Abflussbildung aus. Zum einen sind die Anfangsbedingungen im Einzugsgebiet (Benetzung und Muldenfüllung auf der Oberfläche sowie eingeschränkt auch die Bodenfeuchte bei durchlässigen Flächen) ein Resultat des vor dem betrachteten Niederschlagsereignis stattfindenden Verdunstungsgeschehens, zum anderen wird der zu berechnende abflusswirksame Niederschlag um den Betrag der momentanen Verdunstungsrate geschmälert.

Die potentielle (energetisch mögliche) Verdunstung VP ist zeitlich und örtlich sehr unterschiedlich und einer genauen Berechnung nur schwer zugänglich.

Es stehen verschiedene Optionen zur Eingabe/Berechnung der potentiellen Verdunstung zur Verfügung:

  • Vorgabe einer externen Zeitreihe mit Werten für die potentielle Verdunstung
  • Interne Berechnung der potentiellen Verdunstung: Die potentielle Verdunstung wird auf Basis von Temperatur und je nach Verfahren weiteren Parametern intern berechnet. Es stehen die folgenden Verfahren zur Verfügung:
    • Penman: Erfordert die Angabe mittlerer Jahresgänge oder von Zeitreihen für Sonnenscheindauer, Windgeschwindigkeit und relative Luftfeuchte
    • Haude: Erfordert die Angabe eines mittleren Jahresgangs oder einer Zeitreihe relative Luftfeuchte
    • Turc: Erfordert die Angabe mittlerer Jahresgänge oder von Zeitreihen für Sonnenscheindauer und relative Luftfeuchte
    • Blaney-Criddle: Vereinfachtes Verfahren, erfordert lediglich die Angabe des Breitengrads, woraus die Sonnenscheindauer abgeleitet wird.
  • Angabe einer fixen Verdungstungshöhe pro Jahr: Bei dieser Option wird die eingegebene Jahressumme mit der Ausgleichsfunktion nach Brandt (siehe unten) innerjährlich verteilt.

Ausgleichsfunktionen

Die folgenden Ausgleichsfunktionen nach Brandt werden verwendet, um eingegebene Jahressummen der Verdunstung über das Jahr zu verteilen bzw. um Tageswerte auf Stundenwerte umzurechnen.

Aus ausgewerteten Messungen von 20 Stationen, deren Mittelwerte als Histogramm in dargestellt sind, wurde folgende Ausgleichsfunktion ermittelt /BRANDT, 1979/.

[math]\displaystyle{ VP=(0.96+0.0033 \cdot i) \cdot \sin\frac{2\pi}{365}(i-148)+158 }[/math]

mit:

[math]\displaystyle{ i }[/math]: laufender Tag des Abflussjahres
[math]\displaystyle{ i=1 }[/math]: 1. November

Die jährliche potentielle Gesamtverdunstungshöhe beträgt 642 mm. Liegen keine gemessenen Verdunstungswerte vor, kann optional dieser normierte Jahresgang der potentiellen Verdunstung für die Berechnung der aktuellen Verdunstung herangezogen werden. Ist das gewählte Berechnungszeitintervall kleiner als ein Tag, wird mittels des dargestellten Tagesgangs letztendlich die potentielle Verdunstung für jedes Berechnungszeitintervall ermittelt. Ist das Berechnungsintervall ≥ 1 Tag entfällt die Berücksichtigung des Tagesganges.

Oberflächenwasservorrat (versiegelter Flächenanteil) O

Bei den versiegelten Flächenanteilen kann neben dem Schneevorrat auch die Infiltration vernachlässigt werden, so dass sich die Bilanzgleichung wie folgt vereinfacht:

[math]\displaystyle{ N_W(t)=N(t)-VP(t)-\frac{dO}{dt} }[/math]

wobei die Oberflächenwasservorratsänderung [math]\displaystyle{ dO/dt }[/math] die Benetzung der Oberfläche sowie die Auffüllung und Entleerung (durch Verdunstung) der Mulden repräsentiert.

Schema der Modellansätze Benetzungs- und Muldenverlust

Als Benetzungsverlust [math]\displaystyle{ BV }[/math] für versiegelte Flächen wird folgender Standardwert angesetzt.

[math]\displaystyle{ BV = 0.5 \mbox{ mm} }[/math]

Der Muldenverlust MV wird durch den Anwender vorgegeben. Der Standard- und gleichzeitig Maximalwert im Modell beträgt 4 mm Der Muldenverlust stellt den Mittelwert für eine geneigte Oberfläche dar. Da die Mulden jedoch nicht gleichmäßig verteilt sind und erfahrungsgemäß bereits ein Abfluss einsetzt, bevor überall die komplette Muldenauffüllung erreicht ist, wird unterstellt, dass jeweils

  • 1/3 der versiegelten Fläche einen verminderten Muldenverlust von 1/3·MV
  • 1/3 der versiegelten Fläche den mittleren Muldenverlust von 3/3·MV
  • 1/3 der versiegelten Fläche einen erhöhten Muldenverlust von 5/3·MV

aufweist. Somit kommt es bereits zum Abfluss, wenn der um die Verdunstungsrate verminderte Niederschlag den Benetzungsverlust und 1/3 des Muldenverlustes übersteigt (bei trockener Vorgeschichte). Die o.g. Annahmen sind in der folgenden Abbildung schematisch skizziert.

Der Abflussbeiwert der versiegelten Flächen (nach Abdeckung der Anfangsverluste) wird mit [math]\displaystyle{ \Psi = 1 }[/math] angesetzt. Bei der Festlegung des versiegelten Flächenanteils in einem Teileinzugsgebiet ist zu beachten, dass nicht alle befestigten oder versiegelten Flächen tatsächlich in eine Kanalisation entwässern. Die kontinuierliche Bereitstellung der Benetzungs- und Muldenverluste erfolgt über die laufende Bilanzierung dieser Speicher und der Verdunstung.


Oberflächenwasservorrat (unversiegelter Flächenanteil) O

Der Oberflächenwasservorrat wird über die Bilanzierung eines Verlustspeichers in Abhängigkeit des gewählten Abflussbildungsansatzes berechnet. Einzelheiten dazu finden sich in den folgenden Abschnitten zur Berechnung der Infiltration bzw. abflusswirksamer Niederschlag.


Infiltration bzw. abflusswirksamer Niederschlag I(t), NW(t)

Bei den durchlässigen Flächen kann die Infiltration in den Boden nicht vernachlässigt werden, da diese das Abflussgeschehen entscheidend prägt. Für die Berechnung wurden drei Ansätze im Modell implementiert:

  1. Konstanter Abflussbeiwert [math]\displaystyle{ \Psi }[/math]
  2. Ereignisspezifischer Abflussbeiwert in Anlehnung an das Verfahren des Soil-Conservation-Service (SCS)
  3. Bodenfeuchtesimulation


Konstanter Abflussbeiwert Ψ

Bei Angabe eines [math]\displaystyle{ \Psi }[/math]-Wertes kommt nach Abdeckung der Anfangsverluste (Benetzungs- und Muldenverlust) der übrige Anteil des Niederschlages im Verhältnis des Abflussbeiwertes [math]\displaystyle{ \Psi }[/math] zum Abfluss und zwar unabhängig von der Vorgeschichte und den Merkmalen des Niederschlages (Höhe, Intensität, Dauer). Auf diesen Ansatz sollte nach Möglichkeit verzichtet werden, da hier der Prozess der Abflussbildung nur grob vereinfachend beschrieben wird.


Ereignisspezifischer Abflussbeiwert in Anlehnung an das Verfahren des Soil-Conservation-Service (SCS)

Bei Angabe eines vom Bodentyp und der Bodennutzung abhängigen CN-Wertes (siehe /DVWK, 1991/) lässt sich ein vorgeschichtsabhängiger Anfangsverlust sowie eine vorgeschichtsabhängige Beziehung des Abflussbeiwertes von der bis zum betrachteten Zeitpunkt akkumulierten Niederschlagshöhe formulieren /Zaiss, 1987/; d.h. der Abflussbeiwert wächst mit zunehmendem Niederschlag im Verlauf des Ereignisses an. Die Quantifizierung der Vorgeschichte erfolgt über den 21-Tage-Vorregenindex [math]\displaystyle{ VN }[/math]

[math]\displaystyle{ V_N=\sum_{j=1}^21 C(j)^j \cdot hN_j }[/math]

mit:

[math]\displaystyle{ hN_j }[/math]: Niederschlagshöhe des j-ten Vortags
[math]\displaystyle{ C(j) }[/math]: Faktor, der den Einfluss des j-ten Vortags beschreibt

Der Einfluss der Jahreszeit wird durch einen Jahresgang des Faktors C wiedergegeben.

[math]\displaystyle{ C=0.05 \cdot \sin\frac{2\pi}{365}(i+0.75)+0.85 }[/math]

mit:

[math]\displaystyle{ i }[/math]: laufender Tag des Abflussjahres

Damit schwankt der Wert C zwischen 0.8 < C < 0.9. Hierdurch wird erreicht, dass bei gleichem Vorregen zu unterschiedlichen Jahreszeiten unterschiedliche Vorregenindizes berechnet und damit eine veränderte Abflussbereitschaft in Rechnung gestellt wird. In Abhängigkeit von der auf diese Weise quantifizierten Vorgeschichte kann unter Verwendung der gebietsspezifischen und für mittlere Vorfeuchteverhältnisse gültigen CN-Werte ein aktueller Abflussbeiwert berechnet werden. In der folgenden Abbildung ist für unterschiedliche CN-Werte dargestellt, wie sich der aktuelle Abflussbeiwert in Abhängigkeit von der Vorgeschichte verändert. Da sich im Verlaufe eines Regenereignisses durch die Durchfeuchtung des Bodens die Abflussbereitschaft eines Einzugsgebiets verändert, wird ebenfalls eine Anpassung des Abflussbeiwertes während eines Ereignisses als Funktion der kumulierten Niederschlagshöhe vorgenommen.


Bodenfeuchtesimulation

Landnutzung

Bei der Anwendung der Bodenfeuchtesimulation ist die Angabe von Landnutzungen notwendig. Aus den Angaben zur Landnutzung wird die Durchwurzelungstiefe benötigt, um die Dicke der Durchwurzelungsschicht zu ermitteln. Weitere Parameter der Landnutzung, die zur Berechnung der Interzeption und der Transpiration dienen, sind:

  • Wurzeltiefe
  • Bedeckungsgrad
  • Jahresgang des Bedeckungsgrades
  • Blattflächenindex
  • Jahresgang des Blattflächenindexes

Die Angabe von Haude-Faktoren zur besseren Berücksichtigung der Verdunstung je Landnutzung sind über Eingabe von Jahresgängen beliebig möglich und können den gewünschten Landnutzungen zugeordnet werden.


Bodentyp/ Bodenart
Beispiel der Zusammenfassung der Bodenschichten zu einer programminternen Schicht anhand der Durchwurzelungsschicht
Mit der Bodenfeuchtesimulation berechnete Größen

Die Bodenfeuchtesimulation basiert auf einer nichtlinearen Berechnung der einzelnen Bodenhorizonte. Der Boden wird dabei in verschiedene Horizonte (Schichten) eingeteilt. Jede Schicht wird berechnet und mit den (falls vorhanden) darunter bzw. darüber liegenden Schichten abgeglichen. Als Parameter zur Bodenfeuchteberechnung dienen folgende bodenphysikalischen Größen:

  • Welkepunkt (WP)
  • Feldkapazität (FK)
  • Gesamtporenvolumen (GPV)
  • Gesättigte Leitfähigkeit (kf-Wert)
  • Maximale Infiltrationskapazität (Max.Inf.)
  • Maximale Rate des Kapillaraufstiegs (Max.Kap.)
  • Zuordnung zu einer Bodenart: Sand, Schluff, Ton

Die mögliche Anzahl der Bodenschichten läuft von minimal einer bis maximal sechs. Die Erfahrung zeigte, dass die besten Ergebnisse mit einer Aufteilung in drei Schichten erzielt werden konnten. Aus diesem Grund werden die eingegebenen Schichten programmintern immer in drei Horizonte unterteilt.

  • Infiltrationsschicht (Standarddicke [cm] = 20)
  • Durchwurzelte Schicht (Mindestdicke [cm] = 5)
  • Transportschicht (Mindestdicke [cm] = 5)

Die Berechnung der neuen Bodenkennwerte für die programmintern verwendeten Schichten erfolgt durch eine Gewichtung entsprechend den vorgegebenen original Dicken der Schichten. Im Fall der gesättigten Leitfähigkeit läuft die Berechnung nach dem Prinzip der Erhaltung der Kontinuität der Strömung ab. Bei senkrechter Strömung soll aufgrund der Kontinuität der Strömung die Geschwindigkeit v bei gegebener Durchflussmenge in einer programminternen Schicht denselben Wert besitzen. Damit ist das hydraulische Gefälle nicht mehr konstant.

[math]\displaystyle{ kf_V=\frac{\sum_{i=1}^n d_i}{\sum_{i=1}^n \frac{d_i}{k_i}} }[/math]

mit:

[math]\displaystyle{ d_i }[/math]: anteilige Schichtdicke der jeweiligen Original-Schicht [mm]
[math]\displaystyle{ k_i }[/math]: gesättigte Leitfähigkeit der jeweiligen Original-Schicht [mm/h]
[math]\displaystyle{ kf_V }[/math]: gesättigte Leitfähigkeit der programmintern verwendeten Schicht [mm/h]
Schema der Ermittlung der aktuellen Verdunstung

Auf der Basis der bereichsweisen linearen Abbildung der die Bodenfeuchte beeinflussenden Prozessfunktionen Infiltration, aktuelle Verdunstung (Evaporation + Transpiration), Perkolation, Interflow und Kapillaraufstieg wird für eine Bodenschicht die Wasserbilanzgleichung gelöst. Die Eingangsgröße für die Evaporation und Transpiration ermittelt sich aus der potentiellen Verdunstung. Die zu lösende Gleichung ist:

[math]\displaystyle{ \frac{dBF(t)}{dt}=Inf(t)-Perk(t)-Eva_{akt}(t)-Trans_{akt}(t)-Int(t)+Kap(t) }[/math]

mit:

[math]\displaystyle{ BF(t) }[/math]: aktuelle Bodenfeuchte
[math]\displaystyle{ Inf(t) }[/math]: Infiltration in den Boden
[math]\displaystyle{ Perk(t) }[/math]: Perkolation (Durchsickerung)
[math]\displaystyle{ Eva_{akt}(t) }[/math]: aktuelle Evaporation
[math]\displaystyle{ Trans_{akt}(t) }[/math]: aktuelle Transpiration
[math]\displaystyle{ Int(t) }[/math]: Interflow
[math]\displaystyle{ Kap(t) }[/math]: Kapillaraufstieg
Darstellung ausgewählter Bodenprozessfunktionen

Infiltration, Perkolation, Evaporation, Transpiration, Interflow und Kapillaraufstieg sind von der aktuellen Bodenfeuchte abhängig. In der Simulation wird diese Abhängigkeit durch folgende Funktionsverläufe beschrieben.

[math]\displaystyle{ Inf(BF(t))=a_v \cdot \left(GPV-BF(t) \right)^{1.4}+k_f }[/math]           (Ansatz nach HOLTAN)

[math]\displaystyle{ Perk(BF(t)) = \begin{cases} 0, & BF(t)\le f_{PK} \cdot nFK + WP \\ k_f \cdot \left(\frac{BF(t)-(f_{PK} \cdot nFK +WP)}{GPV-(f_{PK} \cdot nFK +WP)} \right)^{exp,PK}, & BF(t)\gt f_{PK} \cdot nFK + WP \end{cases} }[/math]

(mod. Ansatz nach /OSTROWSKI, 1992/)

[math]\displaystyle{ Eva(BF(t)) = \begin{cases} 0, & BF(t)\le WP \\ f_{Eva} \cdot \left(\frac{BF(t)-WP}{GPV-WP} \right)^{exp,PK}, & BF(t)\gt WP \end{cases} }[/math]

[math]\displaystyle{ Trans(BF(t)) = \begin{cases} 0, & BF(t)\le f_{Trans} \cdot nFK + WP \\ f_{Trans} \cdot \left(\frac{BF(t)-(f_{Trans} \cdot nFK +WP)}{GPV-(f_{Trans} \cdot nFK +WP)} \right)^{exp,PK}, & BF(t)\gt f_{Trans} \cdot nFK + WP \end{cases} }[/math]

mit:

[math]\displaystyle{ a_v }[/math]: Infiltrationsfaktor nach HOLTAN (in Talsim-NG [math]\displaystyle{ a_v=1 }[/math]
[math]\displaystyle{ k_f }[/math]: Durchlässigkeitsbeiwert des gesättigten Bodens
[math]\displaystyle{ nFK }[/math]: nutzbare Feldkapazität ([math]\displaystyle{ nFK=FK-WP }[/math])
[math]\displaystyle{ WP }[/math]: Welkepunkt
[math]\displaystyle{ FK }[/math]: Feldkapazität
[math]\displaystyle{ GPV }[/math]: gesamtes Porenvolumen
[math]\displaystyle{ f_{PK} }[/math]: bodenabhängiger Skalierungfaktor der Perkolationsfunktion
[math]\displaystyle{ exp,PK }[/math]: bodenabhängiger Krümmungsparameter der Perkolationsfunktion
[math]\displaystyle{ f_{Eva} }[/math]: bodenabhängiger Skalierungsfaktor der Evaporationsfunktion
[math]\displaystyle{ f_{Trans} }[/math]: bodenabhängiger Skalierungsfaktor der Transpirationsfunktion
[math]\displaystyle{ exp,Trans }[/math]: Krümmungsparameter der Transpirationsfunktion

Die Programmparameter werden intern berechnet. Der Anwender muss lediglich die Bodenkennwerte kf, WP, FK und GPV angeben. Die Berechnung der Bodenprozesse erfolgt mit einem neuentwickelten Baustein zur Simulation von Speichern.


Elementarflächen

Wird mit der Bodenfeuchtesimulation die Abflussbildung berechnet, wird gleichzeitig das Elementarflächenkonzept angewandt. Ein Einzugsgebietselement wird dabei in beliebig viele hydrologisch homogene Flächen unterteilt, d.h. Flächen gleichen Bodentyps und gleicher Landnutzung. Für jede Elementarfläche gilt genau eine Zuordnung von Landnutzung und Bodentyp. Die aus einer Elementarfläche resultierende Wassermenge wird am Elementausgang angesetzt, d.h. alle Elementarflächen geben unabhängig ihrer Lage im Einzugsgebiet Wasser mit der gleichen zeitlichen Verzögerung ab.


Abflusskonzentration

Die Abflusskonzentration bestimmt die Verzögerung des Oberflächenabflusses aus dem Einzugsgebiet. Es wird eine Parallelspeicherkaskade mit drei Speichern für unbefestigte und eine Kaskade für befestigte Flächen benutzt. Der Abfluss der Komponenten Interflow und Grundwasser wird über einen linearen Einzelspeicher verzögert an den Elementausgang abgegeben.

Berechnung der Abflusskonzentration von Einzugsgebieten