Aangepast Gompertz-model in BMP-tests: vertraging, tempo en potentieel

Het aangepaste Gompertz-model beschrijft de cumulatieve methaanproductie met een sigmoïdale curve: B(t) = P·exp{−exp[(Rm·e/P)·(λ − t) + 1]}, waarbij P het eindpotentieel is, Rm — de maximale productiesnelheid (NmL/g VS/d), λ — de vertragingstijd (d), en e — de basis van de natuurlijke logaritme. Het model wordt toegepast wanneer na het starten van de test de productie met vertraging op gang komt, vervolgens versnelt en uitdooft, dat wil zeggen wanneer een eerstegraadsmodel niet past bij de gegevens.

Wanneer de BMP-curve een S-vorm heeft

• het inoculum moet zich aanpassen aan het substraat (bijv. vetten, substraten met remstoffen),
• het substraat vereist hydrolyse voorafgaand aan methaanproductie (lignocellulosehoudende materialen),
• een lage I/S-verhouding veroorzaakt tijdelijke verzuring en vertraging,
• in de test treden twee afbraakfasen op.

Parameters

Rekenvoorbeeld

Stel P = 300, Rm = 25, λ = 2 d. De coëfficiënt Rm·e/P = 25 × 2,718 / 300 ≈ 0,2265 d⁻¹.

(Berekening voor t = 10: exp(-0,812) ≈ 0,444, exp(-0,444) ≈ 0,641, B = 300 × 0,641 ≈ 192.) De maximale snelheid valt op het buigpunt, waar B = P/e ≈ 110 NmL/g VS; uit de vergelijking volgt dat dit rond t = λ + P/(Rm·e) ≈ 2 + 4,4 ≈ 6,4 d plaatsvindt.

Aanpassing stap voor stap

1. Voer de gegevens (t, B cumulatief) in een hulpmiddel met niet-lineaire regressie in (rekenblad met solver, R, Python).
2. Stel beginwaarden in: P ≈ laatste punt, Rm ≈ grootste helling, λ ≈ snijpunt van de raaklijn met de tijdas.
3. Pas de parameters aan; controleer de convergentie.
4. Beoordeel RMSE, residuen, betrouwbaarheidsintervallen.
5. Vergelijk met het eerstegraadsmodel (bijv. AIC-criterium); kies het eenvoudigere model als het verschil klein is.

Valkuilen

• Overdaad aan parameters. Drie parameters bij een korte curve met weinig datapunten leveren onbetrouwbare schattingen op.
• P als BMP. Als de test voortijdig is beëindigd, is P een extrapolatie; vermeld dit met een toelichting.
• Interpretatie van λ. De vertraging is het gevolg van het inoculum en de I/S-verhouding, niet een eigenschap van het materiaal.
• Vergelijken tussen laboratoria. Kinetische parameters vertonen een grotere spreiding dan B₀/P.
• Geen mechanisme. Het model is empirisch; het verklaart niet waarom er vertraging optreedt.

Praktische toepassing

Het Gompertz-model is nuttig voor het beschrijven van curves in publicaties en rapporten, het beoordelen van de invloed van voorbehandeling op vertraging, en het opsporen van problemen: een verlengde λ voor een bekend substraat in vergelijking met eerdere series kan wijzen op een zwakker inoculum of een remstof. Het model mag niet rechtstreeks worden gebruikt voor het ontwerpen van een vergister, omdat continue omstandigheden verschillen van batchomstandigheden.

Documentatie

In het rapport vermeldt u: de vorm van de vergelijking, het aanpassingsprogramma, de waarden van P, Rm, λ met onzekerheid, aanpassingsmaten en het bereik van de gebruikte gegevens. Als u een tweefasenmodel of een ander model gebruikt, benoem dit dan expliciet. Zo kan de lezer het resultaat reproduceren.

Modelkeuze — korte gids

Aanvullend voorbeeld: effect van voorbehandeling

Stel dat een substraat vóór voorbehandeling P = 280, Rm = 12, λ = 4 d heeft, en na voorbehandeling P = 300, Rm = 22, λ = 1,5 d. De toename van P met ca. 7% is gering, maar de maximale snelheid is bijna verdubbeld en de vertraging is met 2,5 d afgenomen. Voor een installatie met een korte HRT is dit een significant verschil: een groter deel van het potentieel wordt in de vergister benut. Als de voorbehandeling alleen op basis van het uiteindelijke BMP zou worden beoordeeld, zou het voordeel worden onderschat.

Informatie over gegevensbeperkingen

De aanpassing van het Gompertz-model vereist een dichte bemonstering rond het buigpunt, dat wil zeggen wanneer de productie het snelst stijgt. Bij metingen om de paar dagen wordt de parameter Rm met grote onzekerheid geschat. In het onderzoeksprogramma is het raadzaam om in de eerste twee weken frequentere metingen te plannen.

Verschillen ten opzichte van andere sigmoïdale modellen

In de literatuur worden ook het logistische model en het Richards-model gebruikt, die een vergelijkbare vorm hebben maar andere parametriseringen. Het wijzigen van het model verandert de waarden van λ en Rm, dus het vergelijken van parameters is alleen mogelijk binnen één model. Laat u bij de keuze leiden door de aanpassing en eenvoud, niet door het aantal parameters.