Eerste-orde kinetiek in de BMP-test: constante k, B₀ en tijd t90

Een eerste-orde model beschrijft het cumulatieve methaanopbrengst met de vergelijking B(t) = B₀·(1 − e^(−k·t)), waarbij B₀ het eindpotentieel is (NmL CH₄/g VS) en k de snelheidsconstante (1/d). Uit het model wordt de tijd afgeleid die nodig is om 90% van het potentieel te bereiken (t90 = ln 10 / k ≈ 2,3/k). Dit is de eenvoudigste manier om de BMP-curve te beschrijven en vormt de basis voor veel ontwerpberekeningen, maar kent beperkingen: het houdt geen rekening met de lag-fase en gaat uit van een homogene, éénstaps afbraak.

Oorsprong van het model

Aanname: de snelheid van methaanproductie is evenredig met de hoeveelheid nog niet afgebroken materiaal: dB/dt = k·(B₀ − B). De oplossing hiervan levert de bovenstaande exponentiële functie. In procesmodellen (bijv. ADM1) wordt de hydrolyse van polymeren eveneens beschreven met eerste-orde kinetiek, waardoor de constante k uit de BMP een benadering is van de "hydrolyseconstante" onder de gegeven testomstandigheden.

Interpretatie van parameters

De waarden van k voor substraten in de literatuur variëren binnen een breed bereik (van fracties tot enkele 1/d) en zijn afhankelijk van de methode; ze mogen niet worden beschouwd als materiaalconstanten.

Aanpassing

1. Niet-lineaire regressie (kleinste kwadraten): aanpassing van B₀ en k aan de punten (t, B) — voorkeursmethode.
2. Linearisatie: ln(1 − B/B₀) = −k·t; vereist een bekende B₀ en is gevoelig voor fouten bij B dicht bij B₀.
3. Evaluatie van de aanpassing: R², RMSE, residuanalyse (of residuen een systematisch patroon vertonen).
4. Onzekerheid van parameters: betrouwbaarheidsintervallen uit de aanpassing of bootstrap.

Rekenvoorbeeld

Stel B₀ = 320 NmL/g VS en 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.

Voor k = 0,05 d⁻¹ bedraagt t90 ≈ 46 d. Als de HRT van de vergister 40 dagen is, wordt in het eerste geval bij een vaste verblijftijd (zoals in een batchtest) ongeveer 98% van het potentieel benut, in het tweede geval ongeveer 86% (1 − e^(−2) ≈ 0,865). Dit is een vereenvoudiging, omdat een volledig gemengde vergister een verdeling van de verblijftijd heeft, en geen vaste tijd; in een CSTR-vergister bedraagt de afgebroken fractie k·HRT/(1 + k·HRT): voor k = 0,10 en HRT 40 d is dit 4/5 = 80%, voor k = 0,05 is dit 2/3 ≈ 67%. Het verschil tussen de benadering "in een enkele fles" en CSTR is aanzienlijk en verklaart een deel van de discrepanties tussen BMP en technische schaal.

Beperkingen van het model

• Lag-fase. Als de productie aan het begin vertraagd is, past het eerste-orde model slecht; er wordt dan een model met vertraging of het Gompertz-model toegepast.
• Tweefasige curves. Substraten die een fractie makkelijk en moeilijk afbreekbaar materiaal bevatten (bijv. kuilvoer met vezels) vereisen een tweefractiemodel (som van twee exponentiële functies).
• Inhibitie. Effecten van toxiciteit en verzuring vloeien niet voort uit eerste-orde kinetiek.
• Afhankelijkheid van omstandigheden. k verandert met I/S, temperatuur en deeltjesgrootte.

Goede praktijken

1. Pas het model aan op de volledige curve, niet alleen op het eindpunt.
2. Geef B₀, k en hun onzekerheden aan, evenals de aanpassingsmethode.
3. Vergelijk k niet tussen series zonder te controleren of het inoculum en I/S vergelijkbaar waren.
4. Gebruik k voor een benaderende keuze van HRT, maar bevestig dit met continue tests.

Toepassing in de planning

Voor substraten met een vergelijkbare B₀ betekent een hogere k de mogelijkheid van een kortere HRT zonder verlies van opbrengst. Dit is een argument voor voorbehandeling (verkleining, hydrolyse), die meestal k verhoogt, maar niet noodzakelijkerwijs B₀. Als k na voorbehandeling verdubbelt bij dezelfde B₀, kan een kortere HRT bij dezelfde efficiëntie het vergistervolume verkleinen of de verwerkingscapaciteit verhogen.

Snelle checklist

• [ ] Cumulatieve data, niet dagelijkse, vormen de input voor de aanpassing.
• [ ] Alle meetpunten zijn meegenomen, ook die van het begin van de test.
• [ ] Er is gecontroleerd of er geen lag-fase aanwezig is.
• [ ] B₀, k, onzekerheden en R² zijn vermeld.
• [ ] Er is genoteerd welk hulpmiddel is gebruikt.