Cinetica del primo ordine nel test BMP: costante k, B₀ e tempo t90

Il modello del primo ordine descrive il rendimento cumulato di metano mediante l'equazione B(t) = B₀·(1 − e^(−k·t)), dove B₀ rappresenta il potenziale finale (NmL CH₄/g SV), e k è la costante di velocità (1/d). Dal modello si ricava il tempo necessario per raggiungere il 90% del potenziale (t90 = ln 10 / k ≈ 2,3/k). Si tratta del metodo più semplice per descrivere la curva BMP e della base di molti calcoli progettuali, ma presenta dei limiti: non considera la fase lag e assume una degradazione omogenea e monofasica.

Origine del modello

Ipotesi: la velocità di produzione del metano è proporzionale alla quantità di materiale non ancora degradato: dB/dt = k·(B₀ − B). La soluzione fornisce la suddetta funzione esponenziale. Nei modelli di processo (ad es. ADM1), l'idrolisi dei polimeri è anch'essa descritta da una cinetica del primo ordine, pertanto la costante k del BMP rappresenta una "costante di idrolisi" approssimativa nelle condizioni del test.

Interpretazione dei parametri

I valori di k riportati in letteratura per i substrati variano in un ampio intervallo (da frazioni a diversi 1/d) e dipendono dal metodo; non devono essere considerati come costanti materiali.

Adattamento del modello

1. Regressione non lineare (minimi quadrati): adattamento di B₀ e k ai punti (t, B) — metodo preferenziale.
2. Linearizzazione: ln(1 − B/B₀) = −k·t; richiede la conoscenza di B₀ ed è sensibile agli errori quando B è vicino a B₀.
3. Valutazione dell'adattamento: R², RMSE, analisi dei residui (verifica se i residui presentano un andamento sistematico).
4. Incertezza dei parametri: intervalli di confidenza derivati dall'adattamento o dal bootstrap.

Esempio di calcolo

Si assuma B₀ = 320 NmL/g SV e 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.

Per k = 0,05 d⁻¹, t90 ≈ 46 d. Se l'HRT della camera è di 40 giorni, nel primo caso, con un tempo di ritenzione rigido (come nel test batch), verrà sfruttato circa il 98% del potenziale, mentre nel secondo caso circa l'86% (1 − e^(−2) ≈ 0,865). Questa è una semplificazione, poiché una camera a miscelazione completa presenta una distribuzione del tempo di ritenzione, non un tempo rigido; in una camera CSTR, la frazione degradata è k·HRT/(1 + k·HRT): per k = 0,10 e HRT di 40 d, si ottiene 4/5 = 80%, per k = 0,05, 2/3 ≈ 67%. La differenza tra l'approccio "in singola bottiglia" e il CSTR è significativa e spiega parte delle discrepanze tra BMP e scala tecnica.

Limitazioni del modello

• Fase lag. Quando la produzione iniziale è ritardata, il modello del primo ordine si adatta male; si utilizzano modelli con ritardo o quello di Gompertz.
• Curve bifasiche. I substrati contenenti una frazione facilmente e una difficilmente degradabile (ad es. insilato con fibra) richiedono un modello bifrazionale (somma di due esponenziali).
• Inibizione. Gli effetti di tossicità e acidificazione non derivano dalla cinetica del primo ordine.
• Dipendenza dalle condizioni. k varia con il rapporto I/S, la temperatura e la dimensione delle particelle.

Buone pratiche

1. Adatti il modello all'intera curva, non solo al punto finale.
2. Indichi B₀, k e le relative incertezze, oltre al metodo di adattamento.
3. Non confronti k tra serie diverse senza verificare che l'inoculo e il rapporto I/S siano comparabili.
4. Utilizzi k per una scelta approssimativa dell'HRT, ma confermi con test continui.

Applicazione nella pianificazione

Per substrati con B₀ simile, un k più elevato indica la possibilità di un HRT più breve senza perdita di rendimento. Questo è un argomento a favore del pretrattamento (frammentazione, idrolisi), che solitamente aumenta k, ma non necessariamente B₀. Se dopo il pretrattamento k raddoppia a parità di B₀, un HRT più breve alla stessa efficienza può ridurre il volume della camera o aumentare la capacità di trattamento.

Rapida checklist

• [ ] I dati cumulati, non quelli giornalieri, sono l'input per l'adattamento.
• [ ] Sono stati considerati tutti i punti, anche quelli all'inizio del test.
• [ ] È stata verificata l'assenza di una fase lag.
• [ ] Sono stati indicati B₀, k, le incertezze e R².
• [ ] È stato registrato lo strumento utilizzato.