Zmodyfikowany model Gompertza w testach BMP: opóźnienie, tempo i potencjał

Zmodyfikowany model Gompertza opisuje skumulowaną produkcję metanu krzywą sigmoidalną: B(t) = P·exp{−exp[(Rm·e/P)·(λ − t) + 1]}, gdzie P to potencjał końcowy, Rm — maksymalne tempo produkcji (NmL/g VS/d), λ — czas opóźnienia (d), e — podstawa logarytmu naturalnego. Stosuje się go, gdy po rozpoczęciu testu produkcja rusza z opóźnieniem, a potem przyspiesza i wygasa, czyli gdy model pierwszego rzędu nie pasuje do danych.

Kiedy krzywa BMP ma kształt S

• inokulum wymaga adaptacji do substratu (np. tłuszcze, substraty z inhibitorami),
• substrat wymaga hydrolizy przed wytwarzaniem metanu (materiały lignocelulozowe),
• niski I/S powoduje czasowe zakwaszenie i opóźnienie,
• w próbie występują dwa etapy rozkładu.

Parametry

Przykład rachunku

Przyjmijmy P = 300, Rm = 25, λ = 2 d. Współczynnik Rm·e/P = 25 × 2,718 / 300 ≈ 0,2265 d⁻¹.

(Obliczenie dla t = 10: exp(−0,812) ≈ 0,444, exp(−0,444) ≈ 0,641, B = 300 × 0,641 ≈ 192.) Maksymalne tempo przypada na punkt przegięcia, gdzie B = P/e ≈ 110 NmL/g VS; wynika z równania, że następuje ono w okolicy t = λ + P/(Rm·e) ≈ 2 + 4,4 ≈ 6,4 d.

Dopasowanie krok po kroku

1. Wprowadź dane (t, B skumulowane) do narzędzia z regresją nieliniową (arkusz z solverem, R, Python).
2. Ustal wartości początkowe: P ≈ ostatni punkt, Rm ≈ największe nachylenie, λ ≈ punkt przecięcia stycznej z osią czasu.
3. Dopasuj parametry; sprawdź zbieżność.
4. Oceń RMSE, reszty, przedziały ufności.
5. Porównaj z modelem pierwszego rzędu (np. kryterium AIC); wybierz prostszy, jeśli różnica jest mała.

Pułapki

• Nadmiar parametrów. Trzy parametry przy krótkiej krzywej o małej liczbie punktów dają niewiarygodne oszacowania.
• P jako BMP. Jeśli test zakończono przedwcześnie, P będzie ekstrapolacją; podawaj ją z komentarzem.
• Interpretacja λ. Opóźnienie wynika z inokulum i I/S, nie jest cechą materiału.
• Porównywanie między laboratoriami. Parametry kinetyczne wykazują większy rozrzut niż B₀/P.
• Brak mechanizmu. Model jest empiryczny; nie wyjaśnia, dlaczego pojawiło się opóźnienie.

Wykorzystanie w praktyce

Model Gompertza jest przydatny do opisu krzywych w publikacjach i raportach, oceny wpływu obróbki wstępnej na opóźnienie, a także do wykrywania problemów: wydłużone λ dla znanego substratu w porównaniu z poprzednimi seriami może wskazywać na słabsze inokulum lub inhibitor. Nie powinien być stosowany do projektowania komory wprost, bo warunki ciągłe różnią się od okresowych.

Dokumentowanie

W raporcie podaj: postać równania, program dopasowania, wartości P, Rm, λ z niepewnością, miary dopasowania i zakres danych wykorzystanych do dopasowania. Jeśli używasz modelu dwufazowego lub innego, nazwij go jawnie. Dzięki temu czytelnik może odtworzyć wynik.

Wybór modelu — krótka ściąga

Dodatkowy przykład: efekt obróbki wstępnej

Załóżmy, że substrat przed obróbką ma P = 280, Rm = 12, λ = 4 d, a po obróbce P = 300, Rm = 22, λ = 1,5 d. Wzrost P o ok. 7% jest niewielki, ale maksymalne tempo wzrosło niemal dwukrotnie, a opóźnienie skróciło się o 2,5 d. Dla instalacji o krótkim HRT to zmiana istotna: większa część potencjału zostanie wykorzystana w komorze. Gdyby oceniać obróbkę tylko po końcowym BMP, korzyść byłaby niedoszacowana.

Informacje o ograniczeniach danych

Dopasowanie Gompertza wymaga gęstego próbkowania w okolicy punktu przegięcia, czyli wtedy, gdy produkcja rośnie najszybciej. Przy pomiarach co kilka dni parametr Rm będzie szacowany z dużą niepewnością. W programie badań warto przewidzieć częstsze pomiary w pierwszych dwóch tygodniach.

Różnice względem innych modeli sigmoidalnych

W literaturze używa się także modeli logistycznego i Richardsa, które mają podobny kształt, lecz inne parametryzacje. Zmiana modelu zmienia wartości λ i Rm, więc porównywanie parametrów jest możliwe tylko w obrębie jednego modelu. Przy wyborze kieruj się dopasowaniem i prostotą, nie liczbą parametrów.