SCS-Method: Difference between revisions

From BlueM
Jump to navigation Jump to search
(Vorsicht)
m (remove hierarchy extension template)
 
(9 intermediate revisions by 2 users not shown)
Line 1: Line 1:
{{HierarchieKopf}}
{{BlueM.Sim_nav}}
'''<big>Ereignisspezifischer Abflussbeiwert in Anlehnung an das Verfahren des Soil-Conservation-Service (SCS)</big>'''
{{BlueMTheory_nav}}
'''<big>Event specific run-off coefficient on the basis of the Curve-Number-Method (CN-Method)of the Soil-Conservation-Service (SCS)</big>'''
__TOC__
__TOC__
==Theorie==
==Theory==
Das in BlueM verwendete Verfahren ist eine Weiterentwicklung nach {{:Literatur:Zaiß_1989}} der ursprünglichen ''Curve Number Method'' ({{:Literatur:USDA_1964}}).
The Method applied in BlueM is a further development of the CN-Method ({{:Literatur:USDA_1964}}) by {{:Literatur:Zaiß_1989}}


[[Bild:Theorie_Abb36.gif|thumb|Abbildung 36: Abhängigkeit des Abflussbeiwertes von der Vorgeschichte]]
Bei Angabe eines vom Bodentyp und der Bodennutzung abhängigen CN-Wertes (siehe {{:Literatur: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; d.h. der Abflussbeiwert wächst mit zunehmendem Niederschlag im Verlauf des Ereignisses an.


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 [[:Bild:Theorie_Abb36.gif|Abbildung 36]] ist für unterschiedliche CN-Werte dargestellt, wie sich der aktuelle Abflussbeiwert in Abhängigkeit von der Vorgeschichte verändert.


[[Bild:Theorie_Abb37.gif|thumb|Abbildung 37: Abhängigkeit des Abflussbeiwertes von der kumulierten Niederschlagssumme]]
[[Bild:Theorie_Abb36.gif|thumb|Abbildung 36: Dependency of the run-off coefficient on previous history]]
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. In [[:Bild:Theorie_Abb37.gif|Abbildung 37]] ist dieser Zusammenhang für unterschiedliche CN-Werte dargestellt.
By supplying a soil type and land use dependent CN
(refer to{{:Literatur:DVWK_1991}}) an antecedent dependent initial loss as well as an antecedent dependent relationship of the run-off coefficient on the accumulated amount of rain up to a desired point in time can be determined. This results in an increasing run-off coefficient due to the accumulating rain amount in the course of a rainfall event.


<div class="TALSIM">
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)
A current run-off coefficient can be determined in dependency of the quantified previous history, as described above, by using the area specific and for average ''antecedent moisture conditions (AMC)'' ({{:Literatur:USDA_1964}}) valid CN.
# Konstanter Verlustansatz:<br/>Die Anpassung des Verlustwertes erfolgt nur zu Ereignisbeginn einmalig.<br/>(der Ansatz eines Vorregens ist in diesem Fall zweckmäßig)
In [[:Bild:Theorie_Abb36.gif|Abbildung 36]] the development of the run-off coefficient depending on previous history is depicted for different CN.


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.
[[Bild:Theorie_Abb37.gif|thumb|left|Abbildung 37: Dependency of the run-off coefficient on the accumulated amount of rain]]
As soil moisture increases during the course of a rainfall event the conditions for run-off development change which is why the run-off coefficient is additionally adjusted during a rainfall event as a function of accumulated rain. 
This correlation is depicted for different CN in [[:Bild:Theorie_Abb37.gif|Abbildung 37]].  


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.
;Attention:The CN-Method was developed for the simulation of solitary events on a daily basis. A further development is being undertaken for continuous simulation with smaller time steps.(refer to Bug 23 and [[Talk:SCS-Verfahren#Weiterentwicklung_.28Bug_23.29|Discussion]])<br clear="all"/>
</div>


;Vorsicht:Das SCS-Verfahren wurde für Simulationen von Einzelereignissen auf Tagesbasis konzipiert, für kontinuierliche Simulationen mit kleineren Zeitschritten ist eine Weiterentwicklung in Bearbeitung (siehe Bug 23 und [[Diskussion:SCS-Verfahren|Diskussion]])
==Calculation==
 
===one-time calculated parameters===
==Berechnung==
'''Input:''' <code>CN<sub>II</sub></code>
===Einmalig berechnete Parameter===
'''Eingangsgröße:''' <code>CN<sub>II</sub></code>


<div class="comment">
<div class="comment">
Die Umrechnung von <code>CN<sub>II</sub></code> in <code>CN<sub>I</sub></code> bedeutet, dass davon ausgegangen wird, dass das Gebiet zu Beginn der Simulation trocken ist?!
Die Umrechnung von <code>CN<sub>II</sub></code> in <code>CN<sub>I</sub></code> bedeutet, dass davon ausgegangen wird, dass das Gebiet zu Beginn der Simulation trocken ist?!
</div>
</div>
'''Umrechnung''' von <code>CN<sub>II</sub></code> in <code>CN<sub>I</sub></code>:
'''Conversion''' of <code>CN<sub>II</sub></code> to <code>CN<sub>I</sub></code>:
:<math>CN_I = \frac{CN_{II}}{(2.3340 - 0.01334 \cdot CN_{II})}</math>
:<math>CN_I = \frac{CN_{II}}{(2.3340 - 0.01334 \cdot CN_{II})}</math>


'''Maximaler Gebietsrückhalt''' (Speichervermögen) <code>S<sub>max</sub></code> [mm]:
'''Maximum retention capability of the area''' (storage capacity) <code>S<sub>max</sub></code> [mm]:
:<math>S_{max} = \frac{25400}{CN_I} - 254</math>
:<math>S_{max} = \frac{25400}{CN_I} - 254</math>


'''Gebietsspezifischer Anfangsverlust''' <code>I<sub>a</sub></code> [mm]:
'''area specific initial loss''' <code>I<sub>a</sub></code> [mm]:
:<math>I_a = a \cdot S_{max}</math>
:<math>I_a = a \cdot S_{max}</math>
:mit
:with
:<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>a</code> = constant, originally set to <code>0,2</code> {{:Literatur:USDA_1964|}}, adjusted to European conditions in BlueM as <code>0,05</code>{{:Literatur:DVWK_1991|}}


'''Krümmungsparameter''' <code>CVW</code>:
'''Krümmungsparameter''' <code>CVW</code>:
Line 55: Line 51:
</blockquote>
</blockquote>


===kontinuierlich berechnete Parameter===
===continuously calculated parameters===
====Vorgeschichte====
====previous history====
Die Quantifizierung der Vorgeschichte erfolgt über den 21-Tage-Vorregenindex <code>V<sub>N</sub></code>:
Previous history is quantified through the 21-day-antecedent rain index <code>V<sub>N</sub></code>:


:<math>V_N = \sum_{j=0}^{21} C(j)^j \cdot h_{N,j}</math>
:<math>V_N = \sum_{j=0}^{21} C(j)^j \cdot h_{N,j}</math>
:''Gl. 2.1 in {{:Literatur:Zaiß_1989}}''
:''Gl. 2.1 in {{:Literatur:Zaiß_1989}}''


:mit
:with
:<code>h<sub>N,j</sub></code> = Niederschlagshöhe des j-ten Vortags (<code>j = 0</code> ist aktueller Tag)
:<code>h<sub>N,j</sub></code> = rain height of the preceding day j (<code>j = 0</code> is the current day)
:<code>C(j)</code> = Faktor, der den Einfluss des j-ten Vortags beschreibt
:<code>C(j)</code> = factor, which describes the influence of the preceding day j


Der Einfluss der Jahreszeit wird durch einen Jahresgang des Faktors C wiedergegeben.  
Seasonal influence is considered through the annual pattern of factor c.  


:<math>C = 0.85 \cdot \sin\left(\frac{2 \pi}{365}\right) (i + 0.75 ) + 0.85</math>
:<math>C = 0.85 \cdot \sin\left(\frac{2 \pi}{365}\right) (i + 0.75 ) + 0.85</math>
:''Quelle? bei Zaiß finden sich nur so ähnliche Formeln (2.2 & 2.3)''
:''Quelle? bei Zaiß finden sich nur so ähnliche Formeln (2.2 & 2.3)''


:mit
:with
:<code>i</code> = lfd. Tag des Abflussjahres
:<code>i</code> = ongoing day of the hydrological year


Damit schwankt der Wert <code>C</code> zwischen <code>0,8 < C < 0,9</code>. Hierdurch wird erreicht, dass bei gleichem Vorregen zu unterschiedlichen Jahreszeiten unterschiedliche Vorregenindizes berechnet und damit eine veränderte Abflussbereitschaft in Rechnung gestellt wird.
<code>C</code> will alternate between <code>0,8 < C < 0,9</code>. This allows for different antecedent rain indices due to seasonal differences for same antecedent rain amounts and thereby leads to different conditions for run-off development.  


====Ereignisspezifischer Anfangsverlust====
====Event dependent initial loss====
:<math>h_{va} = I_a \cdot e^{-\frac{V_N}{CVW}}</math>
:<math>h_{va} = I_a \cdot e^{-\frac{V_N}{CVW}}</math>
:''Gl. 4.5b in {{:Literatur:Zaiß_1989}}''
:''Gl. 4.5b in {{:Literatur:Zaiß_1989}}''


====Abflussbeiwert====
====run-off coefficient====
[[Bild:SCS PSI.png|thumb|Abhängigkeit von <code>&psi;</code> zu <code>h<sub>va</sub></code> und <code>h<sub>NE</sub></code> (mit <code>A<sub>v</sub> = 0.05</code>)]]
[[Bild:SCS PSI.png|thumb|Dependency of <code>&psi;</code> on <code>h<sub>va</sub></code> and <code>h<sub>NE</sub></code> (with<code>A<sub>v</sub> = 0.05</code>)]]
:<math>\psi = \begin{cases}
:<math>\psi = \begin{cases}
0, & h_{va} \ge h_{NE} \\
0, & h_{va} \ge h_{NE} \\
Line 88: Line 84:
:''Gl. 4.4 in {{:Literatur:Zaiß_1989}}''
:''Gl. 4.4 in {{:Literatur:Zaiß_1989}}''


:mit
:with
: <code>A<sub>v</sub></code> = Verlustverhältnis = <code>0,05</code>
: <code>A<sub>v</sub></code> = loss ratio = <code>0,05</code>
: <code>h<sub>NE</sub></code> = ereignisbezogene Summe des Niederschlags [mm]
: <code>h<sub>NE</sub></code> = event-driven sum of rainfall [mm]


==Literaturangaben==
==Literature==
<references/>
<references/>
{{HierarchieFuss}}


[[Kategorie:BlueM Theorie]]
[[Kategorie:BlueM Theorie]]

Latest revision as of 08:14, 9 January 2015

BlueM_icon.png BlueM.Sim | Downloads | Application | Theory | Development

Event specific run-off coefficient on the basis of the Curve-Number-Method (CN-Method)of the Soil-Conservation-Service (SCS)

Theory

The Method applied in BlueM is a further development of the CN-Method (USDA (1964)[1]) by Zaiß (1989)[2]


Abbildung 36: Dependency of the run-off coefficient on previous history

By supplying a soil type and land use dependent CN (refer toDVWK (1991)[3]) an antecedent dependent initial loss as well as an antecedent dependent relationship of the run-off coefficient on the accumulated amount of rain up to a desired point in time can be determined. This results in an increasing run-off coefficient due to the accumulating rain amount in the course of a rainfall event.


A current run-off coefficient can be determined in dependency of the quantified previous history, as described above, by using the area specific and for average antecedent moisture conditions (AMC) (USDA (1964)[1]) valid CN. In Abbildung 36 the development of the run-off coefficient depending on previous history is depicted for different CN.

Abbildung 37: Dependency of the run-off coefficient on the accumulated amount of rain

As soil moisture increases during the course of a rainfall event the conditions for run-off development change which is why the run-off coefficient is additionally adjusted during a rainfall event as a function of accumulated rain. This correlation is depicted for different CN in Abbildung 37.

Attention
The CN-Method was developed for the simulation of solitary events on a daily basis. A further development is being undertaken for continuous simulation with smaller time steps.(refer to Bug 23 and Discussion)

Calculation

one-time calculated parameters

Input: CNII

Die Umrechnung von CNII in CNI bedeutet, dass davon ausgegangen wird, dass das Gebiet zu Beginn der Simulation trocken ist?!

Conversion of CNII to CNI:

[math]\displaystyle{ CN_I = \frac{CN_{II}}{(2.3340 - 0.01334 \cdot CN_{II})} }[/math]

Maximum retention capability of the area (storage capacity) Smax [mm]:

[math]\displaystyle{ S_{max} = \frac{25400}{CN_I} - 254 }[/math]

area specific initial loss Ia [mm]:

[math]\displaystyle{ I_a = a \cdot S_{max} }[/math]
with
a = constant, originally set to 0,2 [1], adjusted to European conditions in BlueM as 0,05[3]

Krümmungsparameter CVW:

Laut Sartor[4], der die selbe Gleichung verwendet, stammt dieser Ansatz aus der Dokumentation von SMUSI 3.0

[math]\displaystyle{ CVW = \frac{-100.}{\ln(\frac{0.5}{I_a})} }[/math]
entspricht b1 in Gl. 4.5b in Zaiß (1989)[5]
laut Zaiß:

Eine Abhängigkeit des "Krümmungsparameters" b1 von Gebietskenngrößen konnte im Rahmen dieser Arbeit nicht gefunden werden. Sie läßt sich nach den hier aufgeführten Zusammenhängen lediglich über Regressionsanalysen mehrerer N-A-Ereignisse für das jeweils betreffende Einzugsgebiet ermitteln.

continuously calculated parameters

previous history

Previous history is quantified through the 21-day-antecedent rain index VN:

[math]\displaystyle{ V_N = \sum_{j=0}^{21} C(j)^j \cdot h_{N,j} }[/math]
Gl. 2.1 in Zaiß (1989)[6]
with
hN,j = rain height of the preceding day j (j = 0 is the current day)
C(j) = factor, which describes the influence of the preceding day j

Seasonal influence is considered through the annual pattern of factor c.

[math]\displaystyle{ C = 0.85 \cdot \sin\left(\frac{2 \pi}{365}\right) (i + 0.75 ) + 0.85 }[/math]
Quelle? bei Zaiß finden sich nur so ähnliche Formeln (2.2 & 2.3)
with
i = ongoing day of the hydrological year

C will alternate between 0,8 < C < 0,9. This allows for different antecedent rain indices due to seasonal differences for same antecedent rain amounts and thereby leads to different conditions for run-off development.

Event dependent initial loss

[math]\displaystyle{ h_{va} = I_a \cdot e^{-\frac{V_N}{CVW}} }[/math]
Gl. 4.5b in Zaiß (1989)[7]

run-off coefficient

Dependency of ψ on hva and hNE (withAv = 0.05)
[math]\displaystyle{ \psi = \begin{cases} 0, & h_{va} \ge h_{NE} \\ 1 - \left(\frac{h_{va}}{A_v \cdot h_{NE} + (1 - A_v) \cdot h_{va}}\right)^2, & h_{va} \lt h_{NE} \end{cases} }[/math]
Gl. 4.4 in Zaiß (1989)[8]
with
Av = loss ratio = 0,05
hNE = event-driven sum of rainfall [mm]

Literature

  1. 1.0 1.1 1.2 U.S. Department of Agriculture, Soil Conservation Service (1964): National Engineering Handbook, Section 4 Hydrology, Washington
    (überarbeitete Fassung von 2004: NEH Part 630 Ch. 10)
  2. Zaiß, H. (1989): Simulation ereignisspezifischer Einflüsse des Niederschlag-Abfluß-Prozesses von Hochwasserereignissen kleiner Einzugsgebiete mit N-A-Modellen. Technischer Bericht des Instituts für Ingenieurhydrologie und Hydraulik, TH Darmstadt, Nr. 42
  3. 3.0 3.1 DVWK (1991): Beitrag zur Bestimmung des effektiven Niederschlags für Bemessungshochwasser aus Gebietskenngrößen. Ergebnis einer vergleichenden Untersuchung durch den DVWK-Fachausschuß "Niederschlag-Abfluß-Modelle", Materialien, Heft 2
  4. 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)
  5. Zaiß, H. (1989): Simulation ereignisspezifischer Einflüsse des Niederschlag-Abfluß-Prozesses von Hochwasserereignissen kleiner Einzugsgebiete mit N-A-Modellen. Technischer Bericht des Instituts für Ingenieurhydrologie und Hydraulik, TH Darmstadt, Nr. 42
  6. Zaiß, H. (1989): Simulation ereignisspezifischer Einflüsse des Niederschlag-Abfluß-Prozesses von Hochwasserereignissen kleiner Einzugsgebiete mit N-A-Modellen. Technischer Bericht des Instituts für Ingenieurhydrologie und Hydraulik, TH Darmstadt, Nr. 42
  7. Zaiß, H. (1989): Simulation ereignisspezifischer Einflüsse des Niederschlag-Abfluß-Prozesses von Hochwasserereignissen kleiner Einzugsgebiete mit N-A-Modellen. Technischer Bericht des Instituts für Ingenieurhydrologie und Hydraulik, TH Darmstadt, Nr. 42
  8. Zaiß, H. (1989): Simulation ereignisspezifischer Einflüsse des Niederschlag-Abfluß-Prozesses von Hochwasserereignissen kleiner Einzugsgebiete mit N-A-Modellen. Technischer Bericht des Instituts für Ingenieurhydrologie und Hydraulik, TH Darmstadt, Nr. 42