Bilanzieren

Die klassische Physik kennt neben der Energie sieben weitere bilanzierfähige Mengen. In dieser Lektion wird anhand der Grössen Masse und Volumen erklärt, was eine Bilanz ist und wie man eine Bilanz formuliert. Bilanzieren (System abzugrenzen, Einflussgrössen identifizieren, Bilanz formulieren) gehört zu den grundlegenden Fähigkeiten eines gebildeten Menschen. Der praktische Nutzen des Bilanzierens reicht weit über die Physik hinaus. Ökonomen, Ingenieure, Naturwissenschafter, Ärzte und Politiker haben fast täglich mit Bilanzen zu tun. Ob Geld-, Material- oder Informationsflüsse quantitativ zu erfassen sind oder ein Medikament verabreicht werden muss, macht keinen grossen Unterschied. Ausgangspunkt dieser Überlegung ist immer das Bilanzieren. Parallel zu dieser Vorlesung werden die Studierenden mit einer Methode vertraut gemacht, welche das Bilanzieren massiv erleichtert.

Lernziele

In dieser Vorlesung werden sie lernen

wie eine Bilanz bezüglich eines Systems im allgemeinen Fall zu formulieren ist
wie eine einfache Volumenbilanz oder Massenbilanz aussieht
wie aus dem Volumenstrom-Zeit-Diagramm das geflossene Volumen herausgelesen werden kann
wie die Stärke eines Volumenstromes mit der mittleren Strömungsgeschwindigkeit zusammenhängt
wie die Änderungsrate aus dem Inhalt berechnet wird
wie aus dem Füllhöhen-Zeit-Diagramm die Volumenänderungsrate bestimmt wird

Problemstellung

10. August 2007, 19:18; Letzte Aktualisierung: 20:43; Quelle: sf drs

Berner liessen die Aargauer untergehen/Schleusen im Bielersee zu stark offen

Der Kanton Aargau wurde besonders stark von den Überschwemmungen getroffen. Mitverantwortlich ist der Kanton Bern, der die Schleusen beim Bielersee zu stark geöffnet hatte. Das Bundesamt für Umwelt nimmt die Berner in Schutz.

Der Aargauer Gewässerkenner Pierre-Yves Christen rieb sich am Donnerstag die Augen, als er die unglaublichen Zahlen sah: Die Aare floss mit knapp 1300 Kubikmetern pro Sekunde durch Murgenthal (AG) – viel mehr als erlaubt. Denn laut «Murgenthaler Bedingung» darf die Aare bei Murgenthal nicht mehr als 850 Kubikmeter Wasser pro Sekunde führen. Offensichtlich hatte der Kanton Bern die Schleusen am bedrohlich steigenden Bielersee zu stark geöffnet, wie die «Aargauer Zeitung» mutmasste.

Tatsächlich waren die Schleusen in Biel nach Informationen von SF Tagesschau für lange Zeit zu weit offen. Der Kanton Bern hatte sich auf eine Prognose des Bundes gestützt, die den Höchststand des Wassers erst für Donnerstag Vormittag vorsah. Die Wassermassen waren jedoch schneller und die Schleusen konnten nicht mehr rechtzeitig geschlossen werden.

Das Bundesamt für Umwelt übt Selbstkritik: «Jede Prognose ist mit einer Unsicherheit verbunden; wir haben gedacht, die Spitze käme später» , sagt Vizedirektor Andreas Götz. «Die Berner haben sich klar an die Regeln gehalten und sich auf die Prognose verlassen.»

Kann der Bielersee den Aargau vor dem Hochwasser retten?
Wie stark ist die "Murgenthaler Bedingung" verletzt worden?
Wieso braucht es überhaupt eine Prognose?

Bilanz bezüglich des Bielersees

Der Bielersee weist ein Fläche von 39.3 km² auf, besitzt drei grössere Zuflüsse (Hagneck-Kanal, Schüss und Twannbach), ist mit dem Neuenburgersee über den Zihlkanal leitend verbunden und entwässert in den Nidau-Büren-Kanal. Gemessen werden die Volumenströme im Hagneck-, im Zihl- und im Nidau-Büren-Kanal sowie der Schüss (Suze) bei Sonceboz.

Damit sind das Bilanzgebiet und die grösseren Verbindungen mit der Umwelt festgelegt. Als bilanzierbare Menge nehmen wir die Masse des Wassers. Der Bielersee kann seinen Gehalt an Wasse aber auf drei verschiedene Arten ändern

durch Zu- und Abstrom in Flüssen und Bächen
durch Regen, Verdunstung und Austausch mit dem Grundwasser
durch Vernichtung (Wasser kann chemisch gebunden werden) und Erzeugung (Wasser wird beim Abbau von organischem Material erzeugt)

Den Austausch über Flüsse und Bäche messen wir mit Massenstromstärken I_m, den diffusen Austausch über die Oberfläche nennen wir Quelle und messen diese mit der Quellenstärke Σ_m, die Produktion von Wasser ist durch die Produktionsrate Π_m bestimmt. Ein wegfliessender Massenstrom, eine Senke (z.B. durch Verdunstung) oder eine Vernichtung von Wasser auf chemischem Weg wird mit einem negativen Vorzeichen versehen. Strom- und Quellenstärken sowie Produktionsrate ändern den Inhalt des Bielersees. Die Änderung des Inhalts pro Zeit heisst Änderungsrate. Eine Änderungsrate ist demnach gleich der Inhaltsänderung dividiert durch die dazu benötigte Zeit. Die Änderungsrate der Masse wird mit einem Punkt über dem Formelzeichen (m) abgekürzt.

Nun sind wir in der Lage, für den Bielersee eine saubere Massenbilanz zu formulieren: Die Summe über alle Stromstärken, die totale Quellenstärke und die totale Produktionsrate ergeben zu jedem Zeitpunkt die Änderungsrate des Systems. Dies lässt sich kompakt mit einer Formel schreiben

[math]\Sum_i I_{m_i} + \Sigma_m + \Pi_m = \dot m[/math]

Merken Sie sich diese Bilanzgleichung in Worten und mit Formelzeichen. Sie gehört zu den wichtigsten Formeln der Physik.

Statt der Masse kann auch das Volumen bilanziert werden. Der Zusammenhang zwischen der Massen- und der Volumenbilanz lässt sich am Beispiel der Stromstärke gut zeigen. Weil die Dichte gleich Masse durch Volumen ist, berechnet sich der Massenstrom als Dicht mal Volumenstrom

[math]I_m = \rho I_V[/math]

Die Volumenbilanz lautet damit ohne Quellenterm und Produktionsrate wie folgt

[math]\Sum_i I_{V_i} = \dot V[/math]

Die Summe über alle Volumenstromstärken ist gleich der Änderungsrate des Inhalts. Mit dieser Form der Volumenbilanz wollen wir nun weiter arbeiten. Alle Terme, die in einer Bilanzgleichung auftauchen, haben die Einheit Menge pro Zeit. In der Massenbilanz kommen demnach nur Grössen vor, die in Kilogramm pro Sekunde gemessen werden. Alle Grössen der Volumenbilanz sollten in m³/s angegeben werden.

Volumenstrom und tranportiertes Volumen

Die Graphik zeigt, wie viel Wasser während der starken Regenfälle im August 2008 durch die Suze bei Sonceboz geflossen sind. Etwas präziser formuliert, stellt die Graphik das Volumenstrom-Zeit-Diagramm der Suze bei Sonceboz dar. Wie viel Wasser ist nun aber in den drei Tagen vom Jura her in den Bielersee geflossen?

Am Dienstag hat es wohl nur wenig geregnet. Die Schüss hat im Tagesmittel etwa 7 m³/s Wasser geführt. Multipliziert man diese Stromstärke mit der Zeit (86'400 s) sind am Dienstag, den 7. August 2007, etwa 600'000 m³ Wasser von der Schüss durch Sonceboz transportiert worden. Das ist immerhin ein 60 Meter tiefer See mit einem Grundriss von 100 auf 100 Meter.

Wie berechnet man nun die geflossene Wassermenge für Mittwoch, den 8. August 2007? Dafür gibt es mehrere Verfahren. Liegen die Messwerte in Form von Zahlenreihen vor (Volumenstrom und zugehörige Zeitpunkte), multipliziert man jeden Messwert für den Volumenstrom mit dem Zeitabschnitt, der vom jeweiligen Zeitpunkt bis zu dem des nächsten Messwerts reicht. Danach zählt man alle so berechneten Teilvolumen zusammen. Dieses Verfahren kann mit Hilfe von Excel oder Matlab blitzschnell ausgeführt werden.

Zeichnet man nur einzelne Messwerte in das Volumenstrom-Zeit-Diagramm hinein und zieht von jedem Messpunkt eine horizontale Linie bis zum nächsten, erhält man eine Treppenstruktur. Man nimmt dann an, dass die Volumenstromstärke so lange konstant bleibt, bis man einen neuen Messwert kennt. Das Produkt aus konstant gehaltenem Volumenstrom und Zeitabschnitt entspricht dann der Fläche unter einer Stufe. Mit der Summe zählt man die Flächen unter allen Stufen zusammen. Liegt eine dichte Folge von Messwerten vor, sind die Stufen nicht mehr erkennbar. Die Fläche unter allen Stufen nähert sich dann der effektiven Fläche unter der Kurve an.

Nun sind wir in der Lage, die Frage nach der in drei Tagen geflossenen Menge näherungsweise zu beantworten. Wir müssen nur die Fläche unter der Kurve möglichst gut abschätzen. Dazu gibt es wiederum mehrere Verfahren. Am einfachsten unterteilt man die Kurve in ein paar wenige Stufen unterteilen. Dabei ist darauf zu achten, dass die Stufe eine möglichst exakte Flächenverwandlung der ursprünglichen Kurve darstellt. Danach multipliziert man die Höhe der Stufe (mittlere Volumenstromstärke) mal Länge der Stufe (Zeitabschnitt) und zählt alle so berechneten Teilvolumen zusammen. Die Suze hat demnach vom 8. August 0.00 Uhr bis zum 10. August 24.00 Uhr ..m³ Wasser durch Sonceboz transportiert.

Volumenstrom und Strömungsgeschwindigkeit

Der Hagneck-Kanal wurde als Teil der 1. Juragewässerkorrektion zwischen 1875 und 1878 gebaut, um die Aare in den Bielersee zu leiten. Mit den drei Seen als Ausgleichsbecken wollte man die Wassermenge im unteren Teil der Aare regeln und das Berner Seeland vor weiteren Überschwemmungen schützen. Bei der Mündung in den Bielersee ist 1900 ein Stauwehr errichtet worden, welches die 9 Meter Höhenunterschied zur Stromerzeugung ausnutzt.

Das Diagramm zeigt unter anderem die Fliessgeschwindigkeit und den Volumenstrom. Man kann deutlich erkennen, dass bei Hochwasser der Volumenstrom stärker ansteigt als die Fliessgeschwindigkeit. Wie aus dem Pegelstand abzulesen ist, fliesst dann das Wasser durch einen grösseren Querschitt. Das macht Sinn, weil bei hohem Pegel schon bei gleicher Fliessgeschwindigkeit mehr Wasser durchgeführt würde als bei tiefem Pegelstand.

Der Zusammenhang zwischen Geschwindigkeit und Volumenstrom ist ziemlich kompliziert, weil die Strömung nicht überall gleich stark ist. Beschränkt man sich auf die mittlere Strömungsgeschwindigkeit, kann der Zusammenhang mittels einer einfachen Überlegung gefunden werden. Stellt man sich ein Prisma vor, das den Kanal auf einer bestimmten Länge vollständig ausfüllt, ist das zugehörige Volumen (ΔV) gleich Querschnitt (A) mal Länge des Prismas (Δ s). Bewegt sich nun das Prisma durch den Kanal, gilt für den Volumenstrom

[math]I_V = \frac {\Delta V}{\Delta t} = \frac {A \Delta s}{\Delta t} = A v[/math]

Die Stärke des Volumenstromes kann als Querschnitt mal mittlere Strömungsgeschwindigkeit geschrieben werden. Oder umgekehrt, die mittlere Strömungsgeschwindigkeit ist definiert als Volumenstromstärke durch Querschnitt.

Vorausgesetz, die Fliessgeschwindigkeit ist die mittlere Strömungsgeschwindigkeit, kann der wasserführende Querschnitt des Hagneck-Kanals berechnet werden. Am Mittwoch morgen, den 8. August 2007, betrug dieser Querschnitt etwa 230 m². Am Donnerstag ist der Wasser führende Querschnitt auf etwa 270 m² angewachsen. Weil gleichzeitig der Pegel um 2.2 m gestiegen ist, kann die Breite des Hagneck-Kanals auf etwa 18 m geschätzt werden.

Inhalt und Änderungsrate

Die Zu- und Abflüsse verändern den Inhalt eines Systems: die Summe über die Stärken der Volumenströme und Quellen plus die Produktionsrate sind gleich der Änderungsrate des Volumens. Beim Volumen lässt man die Produktionsrate oft weg (nur geringe Kompression, Stoff bleibt erhalten) und unterscheidet nicht zwischen Strom und Quelle. Kennt man dann Volumenströme, ist die Änderungsrate über die Bilanzgleichung gegeben. Doch wie berechnet man die Volumenänderungsrate, wenn nicht alls Ströme gegeben sind? Die Änderungsrate ist als Inhaltsänderung pro Zeit definiert. Kennt man den Inhalt zu zwei verschiedenen Zeitpunkten, kann die mittlere Änderungsrate im Zeitabschnitt Δ t(Differenz der beiden Zeitpunkte) berechnet werden

[math]\dot V_{mittel} = \frac {V_2 - V_1}{t_2 - t_1} = \frac {\Delta V}{\Delta t}[/math]

Gemäss dieser Definition wird die Änderungsrate automatisch negativ, sobald der Inhalt abnimmt.

Weil beim Bielersee bloss der Pegelstand, die Höhe über Meer, gemessen werden, muss noch eine Zusatzüberlegung gemacht werden. Nimmt man an, dass sich die Fläche des Sees mit der Füllhöhe nicht stark ändert, kann der veränderliche Teil des Sees als Prisma modelliert werden. Die Volumenänderung ist dann gleich Seefläche A mal Höhe (Δ h). Daraus folgt

[math]\dot V_{mittel} = \frac {A \Delta h}{\Delta t} = A \dot h_{mittel}[/math]

Die Volumenänderungsrate ist gleich Grundfläche mal Höhenänderungsrate. Bei zylinderförmigen Gefässen gilt dieser Zusammenhang exakt und zu jedem Zeitpunkt. Nun ist die Höhenänderungsrate die Geschwindigkeit, mit welcher der Seespiegel steigt oder fällt. Also ist die Volumenänderungsrate bei prismatischen Gefässen gleich Grundfläche mal Geschwindigkeit des Flüssigkeitsspiegels. Diese Formel entspricht exakt dem Zusammenhang zwischen Volumenstromstärke und Strömungsgeschwindigkeit.

Liegt der Pegelstand in Funktion der Zeit nur als graphische Darstellung vor, liest man zwei Punkte heraus und bildet den Quotienten aus Höhenänderung und zugehörigem Zeitabschnitt. Dieser Quotient beschreibt - bildlich gesprochen - die Steigung einer Geraden durch die beiden ausgewählten Punkte (Sekante). Verschiebt man den einen Punkt gegen den andern, geht die Sekante in die Tangente des festgehaltenen Punktes über. Diese Tangentkonstruktion (Punkt auswählen, Tangente an die Kurve legen, Steigung bestimmen) erlaubt es uns, aus einem beliebigen Höhen-Zeit-Diagramm die Geschwindigkeit zu bestimmen. So kann aus jedem Orts-Zeit-Diagramm (fälschlicherweise oft Weg-Zeit-Diagramm genannt) die Geschwindigkeit herauszulesen.

Liegt direkt ein Volumen-Zeit-Diagramm vor, gilt der analoge Zusammenhang. Um die Änderungsrate zu einem bestimmten Zeitpunkt zu bestimmen, zeichnet man die Tangente an die Kurve. Die zugehörige Steigung entspricht dann der momentanene Volumenänderungsrate.

Volumenbilanz

Volumenstrom und Pegelstand im Zihlkanal

Die Graphik zeigt den Volumenstrom im Zihlkanal. In dieser Verbindung zwischen dem Neuenburger und Bielersee kann das Wasser je nach Pegelstände in beide Richtungen strömen. Fliesst eine Menge wechselseitig durch einen Leitungsabschnitt, muss das Vorzeichen der Stromstärke festgelegt werden. Dazu führt man einen Bezugspfeil längs der Leitung ein. Verbindet diese Leitung zwei Speicher, überträgt sich die Vorzeichenregel auf den einen Speicher. Bezüglich des andern ist die gleiche Stromstärke mit dem entgegengesetzten Vorzeichen zu bilanzieren. Beim Zihlkanal wird ein Strom, der vom Neuenburger zum Bielersee fliesst, positiv gezählt. Fliesst nun das Wasser wie am 10. August 2007 vom Bielersee in den Neuenburgersee, hat der Volumenstrom eine negative Stärke. Dieser Wert wird in der Volumenbilanz bezüglich des Bielersees direkt und in die Bilanz bezüglich des Neuenburgersees mit umgekehrtem Vorzeichen verrechnet. Die systemdynamischen Modellierungstools wenden die Vorzeichenregel genau so an. Zeigt eine Flussgrösse (Pipeline) in die Bestandesgrösse (Topf) hinein, wird eine gegebene Stromstärke direkt in die Bilanz eingefügt. Zeigt die Flussgrösse hinaus, wird in der Bilanz die zugehörige Stromstärke subtrahiert.

Nun soll die Bilanz bezüglich des Bielersees zum Zeitpunkt 8. August 24 Uhr bzw. 9. August 0 Uhr ausgewertet werden. Aus dem Hagneckkanal ergossen sich 1170 m³/s in den Bielersee, die Schüss brachte zusätzlich 45 m³/s und aus dem Ziehlkanal flossen 140 m³/s gegen den Neuenburgersee. Der Pegel des Bilersees stieg mit einer Geschwindigkeit von 1.9 10^-5 m/s, was auf den ganzen See hochgerechnet eine Änderungsrate von 760 m³/s ergibt. Folglich mussten m³/s 315 m³/s über den Nidau-Büren-Kanal abgeflossen sein. Die Messstelle in Brügg hat einen Wert von 280 m³/s festgestellt. Wie können wir die Differenz erklären?

Drei Gründe für diesen Unterschied zwischen theoretischer Vorhersagen (Bilanzgleichung)und der Messung können aufgeführt werden:

weitere Zuflüsse und Quellenterme (Dauerregen)
Grösse des Sees (kein homogenes System)
Laufzeit der Flutwelle

Weitere Zuflüsse abzuschätzen dürfte eine sehr aufwendige Arbeit sein, aber insgesamt kaum mehr als den Volumenstrom der Suze ausmachen. Ein kräftiger Regen kann 10 bis 50 mm Niederschlag pro Stunde bringen. Geht man von 20 mm aus, gibt das auf den ganzen Bielersee hochgerechnet gut 200 m³/s.

Der Pegel wird bei Ligerz gemessen. Weil die Seeoberfläche durch den Wind und das zufliessende Wasser bewegt wird, könnte man an einer andern Stelle des Sees andere Werte messen.

Ergiesst sich plötzlich ein starker Strom in den See, dauert es eine gewisse Zeit, bis sich dieses Wasser auf den ganzen See verteilt hat. Zudem müsste man noch wissen wie weit die Messstelle Brügg von der Mündung entfernt ist. Je weiter flussabwärts eine Messstelle gelegen ist, umso länger geht es, bis sich dort eine Änderung der Stromstärke bemerkbar macht.

Modellbildung

Durch den Nidau-Büren-Kanal sind in der kritischen Phase maximal 500 m³/s geflossen. Dieser Strom hätte die "Murgenthaler-Bedingung" nicht verletzen können. Ander Flüsse haben das ihre beigetragen. So hat in der fraglichen Nacht allein die Emme 600 m³/s in die Aare eingebracht. Die Idee wäre aber gewesen, dass man den Ausfluss aus dem Bielersee so weit drosselt, dass trotz der Emme und paar weiterer Flüssen die "Murgenthaler-Bedingung" erfüllt geblieben wäre. Weil alle diese Wassermengen eine gewisse Zeit brauchen, bis sie den Kanton Aargau erreicht haben, muss man das ganze hydraulische System der Aare und ihrer Zuflüsse in ein Simulationsmodell abbilden.