Talk:SCS-Method: Difference between revisions

From BlueM
Jump to navigation Jump to search
(Die Seite wurde neu angelegt: ==SCS in BlueM== ===einmalig berechnete Parameter=== :<math>CNI = \frac{CN}{(2.3340 - 0.01334 \cdot CN)}</math> was soll das??? Maximaler Gebietsrückhalt (Speicherver...)
 
 
(23 intermediate revisions by 2 users not shown)
Line 1: Line 1:
==SCS in BlueM==
==SCS in BlueM==
===einmalig berechnete Parameter===
Eine Excel-Datei zur Berechnung des SCS-Verfahrens, wie es unter [[SCS-Verfahren]] beschrieben ist, findet sich unter {{file|xls|SCSVerfahren.xls|SCSVerfahren}}.
:<math>CNI = \frac{CN}{(2.3340 - 0.01334 \cdot CN)}</math> was soll das???


Maximaler Gebietsrückhalt (Speichervermögen) <code>S<sub>max</sub></code> [mm]:
==SCS in TALSIM==
:<math>S_{max} = \frac{25400}{CNI} - 254</math>
<div class="comment">In BlueM ist nur der variable Verlustansatz implementiert. Auch einen Endabflussbeiwert gibt es nicht.<br/>- [[Benutzer:Froehlich|Froehlich]] 13:56, 19. Nov. 2008 (UTC)</div>
<blockquote>
Bei der Abhängigkeit der Abflussbereitschaft zum kumulierten Niederschlag bietet TALSIM 2.2 zwei Möglichkeiten:
 
# Variabler Verlustansatz (default):<br/>Die Anpassung eines Verlustwertes für die Funktion des Abflussbeiwertes zum kumulierten Niederschlag erfolgt für jeden Zeitschritt neu.<br/>(liefert insgesamt höhere Abflussbeiwerte, so dass auf die Berücksichtigung eines Vorregens verzichtet werden kann)
# Konstanter Verlustansatz:<br/>Die Anpassung des Verlustwertes erfolgt nur zu Ereignisbeginn einmalig.<br/>(der Ansatz eines Vorregens ist in diesem Fall zweckmäßig)
 
Welcher Ansatz bessere Ergebnisse liefert geht nur aus einem Vergleich mit gemessenen Ganglinien hervor. Grundsätzlich ergeben sich mit dem variablen Verlustansatz höhere Abflussspitzen und Füllen bei gleichen Bedingungen.
 
Eine weitere Möglichkeit zur Beeinflussung der Abflussbildung besteht in der Option, einen Endabflussbeiwert festzulegen. Damit beschränkt man unabhängig vom gewählten Verlustansatz den maximalen Abflussbeiwert. Standardmäßig setzt TALSIM den Endabflussbeiwert auf 1.
</blockquote>
:''&ndash; {{:Literatur:TALSIM_2.2}}''
 
==Weiterentwicklung (Bug 23)==
;Motiv:
:Das SCS-Verfahren war ursprünglich für einzelne Regenereignisse konzipiert. Sobald in einem Zeitschritt kein Niederschlag mehr fällt, wird das vorhergehende Ereignis als abgeschlossen betrachtet und <code>h<sub>NE</sub></code> wird wieder auf 0 zurückgesetzt. Das führt dazu, dass auch nach nur kurzen Niederschlagspausen der Abflussbeiwert wieder herabgesetzt wird.
 
===Ansatz 1: Zwischenverlust===
[[Bild:SCS Zwischenverlust.png|thumb|Aufbau des Zwischenverlustes ZV in Regenpausen (t<sub>p</sub>){{:Literatur:Sartor_1999|}}]]
''(implementiert in Rev 547)''<br/>
Es wird ein Zwischenverlust nach {{:Literatur:Sartor_1999|Sartor}} <code>ZV</code> eingeführt, der am Anfang der Regenpause 0 beträgt, und der sich mit zunehmender Dauer der Regenpause asymptotisch dem Wert des Anfangsverlusts <code>h<sub>va</sub></code> annähert:
 
:<math>ZV_i = h_{va} - (h_{va} - ZV_{i-1}) \cdot e^{-K_p \cdot dt}</math>
:''aus {{:Literatur:Sartor_1999}}''


Anfangsverlust <code>I<sub>a</sub></code> [mm]:
:<math>I_a = a \cdot S_{max}</math>
:mit
:mit
:<code>a</code> = Konstante, ursprünglich mit <code>0,2</code> angenommen ({{:Literatur:USDA_1964}}), in BlueM für europäische Verhältnisse angepasst an <code>0,05</code> ({{:Literatur:DVWK_1991}})
:<code>K<sub>p</sub></code> = Regenerationskonstante [1/h]
:<code>dt</code> = Simulationszeitschritt [h]
:<code>i</code> = laufendes Intervall über die Regenpause [-]
 
:und
:<code>ZV<sub>i=0</sub> = 0</code>
 
:laut {{:Literatur:Sartor_1999}}:
<blockquote>
Die Größenordnung der zu kalibrierenden Regenerationskonstante K<sub>p</sub> liegt nach bisherigen Erfahrungen zwischen 0,01 und 0,2 h.
</blockquote>
 
Zusätzlich bietet es sich an, das Ende eines Regenereignisses (d.h. den Zeitpunkt an dem <code>h<sub>NE</sub></code> wieder zurückgesetzt wird) mithilfe des Zwischenverlustes zu redefinieren:  
 
:<code>h<sub>NE</sub></code> wird erst dann zurückgesetzt, wenn der Zwischenverlust 90% des Anfangsverlusts erreicht hat (<code>ZV >= 0.9 h<sub>va</sub></code>)


:<math>CVW = \frac{-100.}{\log(\frac{0.5}{I_a})}</math>
===Ansatz 2: Kontinuierlicher Verlauf der Regensumme===
Der Anfangsverlust <code>h<sub>va</sub></code> bleibt unverändert, dafür wird die Regensumme <code>h<sub>NE</sub></code> bei Regenende nicht auf 0 zurückgesetzt, sondern allmählich reduziert.


===Vorwert===
Dies kann entweder proportional zur Abnahme des Inhalts der Kaskadenspeicher geschehen (Rev 548):
Die Quantifizierung der Vorgeschichte erfolgt über den 21-Tage-Vorregenindex <code>VN</code>


:<math>V_N = \sum_{j=1}^{21} C(j)^j \cdot h_{N,j}</math>
:<math>h_{NE}(t) = h_{NE}(t-1) \cdot \left( 1. - \frac{S(t-1)}{S(t)} \right)</math>


:mit
:mit
:<code>h<sub>N,j</sub></code> = Niederschlagshöhe des j-ten Vortags
:<code>S</code> = Inhalt der Speicherkaskaden für unbefestigte Flächen
:<code>C(j)</code> = Faktor, der den Einfluss des j-ten Vortags beschreibt


Der Einfluss der Jahreszeit wird durch einen Jahresgang des Faktors C wiedergegeben.
(Dies birgt aber Probleme in der Implementierung, da der aktuelle Speicherinhalt abhängig von <code>h<sub>NE</sub></code> ist.)


:<math>C = 0.85 \cdot \sin\left(\frac{2 \pi}{365}\right) (i + 0.75 ) + 0.85</math>
Oder aber ganz einfach exponentiell (Rev 549):


:<div class="achtung">0.85 anstatt 0.05!</div>
:<math>h_{NE}(t) = h_{NE}(t-1) \cdot e^{-\frac{t}{K_R}}</math>
 
:mit
:<code>K<sub>R</sub></code> = Konstante für die Abnahme der Regensumme [h]


:mit
Ein guter Anhaltspunkt für <code>K<sub>R</sub></code> scheint der Wert von <code>K<sub>2</sub></code> sein (siehe http://bugs.bluemodel.org/attachment.cgi?id=47)
:<code>i</code> = lfd. Tag des Abflussjahres


===Verlust===
==Literaturangaben==
:<math>h_v = I_a \cdot \exp(-\frac{V_N}{CVW})</math>
<references/>

Latest revision as of 01:40, 5 March 2013

SCS in BlueM

Eine Excel-Datei zur Berechnung des SCS-Verfahrens, wie es unter SCS-Verfahren beschrieben ist, findet sich unter SCSVerfahren information.png.

SCS in TALSIM

In BlueM ist nur der variable Verlustansatz implementiert. Auch einen Endabflussbeiwert gibt es nicht.
- Froehlich 13:56, 19. Nov. 2008 (UTC)

Bei der Abhängigkeit der Abflussbereitschaft zum kumulierten Niederschlag bietet TALSIM 2.2 zwei Möglichkeiten:

  1. Variabler Verlustansatz (default):
    Die Anpassung eines Verlustwertes für die Funktion des Abflussbeiwertes zum kumulierten Niederschlag erfolgt für jeden Zeitschritt neu.
    (liefert insgesamt höhere Abflussbeiwerte, so dass auf die Berücksichtigung eines Vorregens verzichtet werden kann)
  2. Konstanter Verlustansatz:
    Die Anpassung des Verlustwertes erfolgt nur zu Ereignisbeginn einmalig.
    (der Ansatz eines Vorregens ist in diesem Fall zweckmäßig)

Welcher Ansatz bessere Ergebnisse liefert geht nur aus einem Vergleich mit gemessenen Ganglinien hervor. Grundsätzlich ergeben sich mit dem variablen Verlustansatz höhere Abflussspitzen und Füllen bei gleichen Bedingungen.

Eine weitere Möglichkeit zur Beeinflussung der Abflussbildung besteht in der Option, einen Endabflussbeiwert festzulegen. Damit beschränkt man unabhängig vom gewählten Verlustansatz den maximalen Abflussbeiwert. Standardmäßig setzt TALSIM den Endabflussbeiwert auf 1.

– SYDRO (2000)[1]

Weiterentwicklung (Bug 23)

Motiv
Das SCS-Verfahren war ursprünglich für einzelne Regenereignisse konzipiert. Sobald in einem Zeitschritt kein Niederschlag mehr fällt, wird das vorhergehende Ereignis als abgeschlossen betrachtet und hNE wird wieder auf 0 zurückgesetzt. Das führt dazu, dass auch nach nur kurzen Niederschlagspausen der Abflussbeiwert wieder herabgesetzt wird.

Ansatz 1: Zwischenverlust

Aufbau des Zwischenverlustes ZV in Regenpausen (tp)[2]

(implementiert in Rev 547)
Es wird ein Zwischenverlust nach Sartor[2] ZV eingeführt, der am Anfang der Regenpause 0 beträgt, und der sich mit zunehmender Dauer der Regenpause asymptotisch dem Wert des Anfangsverlusts hva annähert:

[math]\displaystyle{ ZV_i = h_{va} - (h_{va} - ZV_{i-1}) \cdot e^{-K_p \cdot dt} }[/math]
aus Sartor (1999)[2]
mit
Kp = Regenerationskonstante [1/h]
dt = Simulationszeitschritt [h]
i = laufendes Intervall über die Regenpause [-]
und
ZVi=0 = 0
laut Sartor (1999)[2]:

Die Größenordnung der zu kalibrierenden Regenerationskonstante Kp liegt nach bisherigen Erfahrungen zwischen 0,01 und 0,2 h.

Zusätzlich bietet es sich an, das Ende eines Regenereignisses (d.h. den Zeitpunkt an dem hNE wieder zurückgesetzt wird) mithilfe des Zwischenverlustes zu redefinieren:

hNE wird erst dann zurückgesetzt, wenn der Zwischenverlust 90% des Anfangsverlusts erreicht hat (ZV >= 0.9 hva)

Ansatz 2: Kontinuierlicher Verlauf der Regensumme

Der Anfangsverlust hva bleibt unverändert, dafür wird die Regensumme hNE bei Regenende nicht auf 0 zurückgesetzt, sondern allmählich reduziert.

Dies kann entweder proportional zur Abnahme des Inhalts der Kaskadenspeicher geschehen (Rev 548):

[math]\displaystyle{ h_{NE}(t) = h_{NE}(t-1) \cdot \left( 1. - \frac{S(t-1)}{S(t)} \right) }[/math]
mit
S = Inhalt der Speicherkaskaden für unbefestigte Flächen

(Dies birgt aber Probleme in der Implementierung, da der aktuelle Speicherinhalt abhängig von hNE ist.)

Oder aber ganz einfach exponentiell (Rev 549):

[math]\displaystyle{ h_{NE}(t) = h_{NE}(t-1) \cdot e^{-\frac{t}{K_R}} }[/math]
mit
KR = Konstante für die Abnahme der Regensumme [h]

Ein guter Anhaltspunkt für KR scheint der Wert von K2 sein (siehe http://bugs.bluemodel.org/attachment.cgi?id=47)

Literaturangaben

  1. ↑ SYDRO (2000): Dokumentation: TALSIM 2.2, Theoretische Grundlagen. http://www.sydro.de/down/talsim22/talsim_theorie.pdf
  2. ↑ 2.0 2.1 2.2 2.3 Sartor, J. (1999): Einsatz der Langzeit-Seriensimulation für kleine Einzugsgebiete, In: Berichte des Fachgebietes Wasserbau und Wasserwirtschaft der Universität Kaiserslautern, Heft 9 (PDF)
Retrieved from "https://wiki.bluemodel.org/index.php?title=Talk:SCS-Method&oldid=6963"