Kinetyka pierwszego rzędu w teście BMP: stała k, B₀ i czas t90

Model pierwszego rzędu opisuje skumulowany uzysk metanu równaniem B(t) = B₀·(1 − e^(−k·t)), gdzie B₀ to potencjał końcowy (NmL CH₄/g VS), a k — stała tempa (1/d). Z modelu odczytuje się czas potrzebny na osiągnięcie 90% potencjału (t90 = ln 10 / k ≈ 2,3/k). Jest to najprostszy sposób opisu krzywej BMP i podstawa wielu obliczeń projektowych, ale ma ograniczenia: nie uwzględnia fazy lag i zakłada jednorodny, jednostopniowy rozkład.

Pochodzenie modelu

Założenie: tempo produkcji metanu jest proporcjonalne do ilości jeszcze nierozłożonego materiału: dB/dt = k·(B₀ − B). Rozwiązanie daje powyższą funkcję wykładniczą. W modelach procesu (np. ADM1) hydroliza polimerów jest również opisywana kinetyką pierwszego rzędu, więc stała k z BMP jest przybliżoną „stałą hydrolizy” w danych warunkach testu.

Interpretacja parametrów

Wartości k dla substratów w literaturze mieszczą się w szerokim zakresie (od ułamków do kilku 1/d) i zależą od metody; nie należy traktować ich jako stałych materiałowych.

Dopasowanie

1. Regresja nieliniowa (najmniejsze kwadraty): dopasowanie B₀ i k do punktów (t, B) — preferowane.
2. Linearyzacja: ln(1 − B/B₀) = −k·t; wymaga znanego B₀ i jest wrażliwa na błędy przy B bliskim B₀.
3. Ocena dopasowania: R², RMSE, analiza reszt (czy reszty mają systematyczny kształt).
4. Niepewność parametrów: przedziały ufności z dopasowania lub bootstrapu.

Przykład obliczeniowy

Załóżmy B₀ = 320 NmL/g VS i k = 0,10 d⁻¹.

• 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.

Dla k = 0,05 d⁻¹ t90 ≈ 46 d. Jeśli HRT komory wynosi 40 dni, w pierwszym przypadku przy idealnym mieszaniu zostanie wykorzystane około 98% potencjału, w drugim — około 86% (1 − e^(−2) ≈ 0,865). To uproszczenie, bo komora o pełnym wymieszaniu ma rozkład czasu przebywania, a nie sztywny czas; w komorze CSTR frakcja rozłożona wynosi k·HRT/(1 + k·HRT): dla k = 0,10 i HRT 40 d to 4/5 = 80%, dla k = 0,05 – 2/3 ≈ 67%. Różnica między podejściem „w pojedynczej butelce” a CSTR jest znaczna i tłumaczy część rozbieżności BMP vs skala techniczna.

Ograniczenia modelu

• Faza lag. Gdy na początku produkcja jest opóźniona, model pierwszego rzędu dopasuje się źle; stosuje się model z opóźnieniem lub Gompertza.
• Krzywe dwustopniowe. Substraty zawierające frakcję łatwo i trudno rozkładalną (np. kiszonka z włóknem) wymagają modelu dwufrakcyjnego (suma dwóch wykładniczych).
• Inhibicja. Efekty toksyczności i zakwaszenia nie wynikają z kinetyki pierwszego rzędu.
• Zależność od warunków. k zmienia się z I/S, temperaturą, wielkością cząstek.

Dobre praktyki

1. Dopasowuj model do całej krzywej, nie tylko do punktu końcowego.
2. Podawaj B₀, k i ich niepewności oraz metodę dopasowania.
3. Nie porównuj k między seriami bez sprawdzenia, że inokulum i I/S były porównywalne.
4. Wykorzystuj k do przybliżonego doboru HRT, ale potwierdzaj testami ciągłymi.

Zastosowanie w planowaniu

Dla substratów o zbliżonym B₀ szybsze k oznacza możliwość krótszego HRT bez straty uzysku. To argument za obróbką wstępną (rozdrobnienie, hydroliza), która zwykle zwiększa k, a niekoniecznie B₀. Jeśli po obróbce k rośnie dwukrotnie przy tym samym B₀, krótszy HRT przy tej samej wydajności może zwolnić objętość komory lub zwiększyć przerób.

Szybka lista kontrolna

• [ ] Dane skumulowane, nie dobowe, są wejściem do dopasowania.
• [ ] Uwzględniono wszystkie punkty, także te z początku testu.
• [ ] Sprawdzono, czy nie ma fazy lag.
• [ ] Podano B₀, k, niepewności i R².
• [ ] Zapisano, jakiego narzędzia użyto.