Kinetik erster Ordnung im BMP-Test: Konstante k, B₀ und Zeit t90

Das Modell erster Ordnung beschreibt die kumulierte Methanausbeute mit der Gleichung B(t) = B₀·(1 − e^(−k·t)), wobei B₀ das Endpotenzial (NmL CH₄/g oTS) und k die Geschwindigkeitskonstante (1/d) ist. Aus dem Modell lässt sich die Zeit ablesen, die zum Erreichen von 90% des Potenzials nötig ist (t90 = ln 10 / k ≈ 2,3/k). Es ist die einfachste Art, die BMP-Kurve zu beschreiben, und die Grundlage vieler Auslegungsrechnungen, hat aber Grenzen: Es berücksichtigt keine Lag-Phase und setzt einen homogenen, einstufigen Abbau voraus.

Herkunft des Modells

Annahme: Die Geschwindigkeit der Methanproduktion ist proportional zur Menge des noch nicht abgebauten Materials: dB/dt = k·(B₀ − B). Die Lösung ergibt die oben genannte Exponentialfunktion. In Prozessmodellen (z. B. ADM1) wird die Hydrolyse von Polymeren ebenfalls mit einer Kinetik erster Ordnung beschrieben, daher ist die Konstante k aus dem BMP eine angenäherte „Hydrolysekonstante“ unter den jeweiligen Testbedingungen.

Interpretation der Parameter

Die Werte von k für Substrate liegen in der Literatur in einem weiten Bereich (von Bruchteilen bis zu einigen 1/d) und hängen von der Methode ab; sie sind nicht als Materialkonstanten zu behandeln.

Anpassung

1. Nichtlineare Regression (kleinste Quadrate): Anpassung von B₀ und k an die Punkte (t, B) — bevorzugt.
2. Linearisierung: ln(1 − B/B₀) = −k·t; erfordert ein bekanntes B₀ und ist bei B nahe B₀ empfindlich gegenüber Fehlern.
3. Beurteilung der Anpassung: R², RMSE, Residuenanalyse (ob die Residuen eine systematische Form haben).
4. Unsicherheit der Parameter: Konfidenzintervalle aus der Anpassung oder dem Bootstrap.

Rechenbeispiel

Nehmen wir B₀ = 320 NmL/g oTS und k = 0,10 d⁻¹ an.

• t = 5 d: B = 320·(1 − e^(−0,5)) = 320 × 0,3935 ≈ 126 NmL/g,
• t = 10 d: B = 320·(1 − e^(−1)) = 320 × 0,632 ≈ 202,
• t = 20 d: B = 320·(1 − e^(−2)) = 320 × 0,8647 ≈ 277,
• t90 = 2,303 / 0,10 ≈ 23 d.

Für k = 0,05 d⁻¹ ist t90 ≈ 46 d. Beträgt die HRT des Fermenters 40 Tage, wird im ersten Fall bei starrer Verweilzeit (wie im Batch-Test) etwa 98% des Potenzials ausgenutzt, im zweiten — etwa 86% (1 − e^(−2) ≈ 0,865). Das ist eine Vereinfachung, da ein vollständig durchmischter Fermenter eine Verweilzeitverteilung hat und keine starre Zeit; in einem CSTR-Fermenter beträgt der abgebaute Anteil k·HRT/(1 + k·HRT): für k = 0,10 und HRT 40 d sind es 4/5 = 80%, für k = 0,05 – 2/3 ≈ 67%. Der Unterschied zwischen dem Ansatz „in der einzelnen Flasche“ und dem CSTR ist erheblich und erklärt einen Teil der Abweichungen zwischen BMP und technischem Maßstab.

Grenzen des Modells

• Lag-Phase. Ist die Produktion zu Beginn verzögert, passt sich das Modell erster Ordnung schlecht an; es wird ein Modell mit Verzögerung oder das Gompertz-Modell verwendet.
• Zweistufige Kurven. Substrate mit einer leicht und einer schwer abbaubaren Fraktion (z. B. Silage mit Fasern) erfordern ein Zweifraktionsmodell (Summe zweier Exponentialfunktionen).
• Inhibition. Toxizitäts- und Versauerungseffekte ergeben sich nicht aus der Kinetik erster Ordnung.
• Abhängigkeit von den Bedingungen. k ändert sich mit I/S, Temperatur und Partikelgröße.

Gute Praxis

1. Passen Sie das Modell an die gesamte Kurve an, nicht nur an den Endpunkt.
2. Geben Sie B₀, k und deren Unsicherheiten sowie die Anpassungsmethode an.
3. Vergleichen Sie k nicht zwischen Serien, ohne zu prüfen, dass Inokulum und I/S vergleichbar waren.
4. Nutzen Sie k zur näherungsweisen Wahl der HRT, bestätigen Sie sie aber mit kontinuierlichen Tests.

Anwendung in der Planung

Bei Substraten mit ähnlichem B₀ bedeutet ein schnelleres k die Möglichkeit einer kürzeren HRT ohne Ausbeuteverlust. Das ist ein Argument für eine Vorbehandlung (Zerkleinerung, Hydrolyse), die k in der Regel erhöht, nicht unbedingt B₀. Steigt k nach der Behandlung bei gleichem B₀ auf das Doppelte, kann eine kürzere HRT bei gleicher Leistung Fermentervolumen freisetzen oder den Durchsatz erhöhen.

Schnelle Checkliste

• [ ] Kumulierte, nicht tägliche Daten sind die Eingabe für die Anpassung.
• [ ] Alle Punkte wurden berücksichtigt, auch die vom Testbeginn.
• [ ] Geprüft, ob keine Lag-Phase vorliegt.
• [ ] B₀, k, Unsicherheiten und R² angegeben.
• [ ] Dokumentiert, welches Werkzeug verwendet wurde.