Cinétique du premier ordre dans le test BMP : constante k, B₀ et temps t90

Le modèle du premier ordre décrit le rendement cumulé en méthane par l’équation B(t) = B₀·(1 − e^(−k·t)), où B₀ est le potentiel final (NmL de CH₄/g de MOS) et k — la constante de vitesse (1/j). Le modèle permet de lire le temps nécessaire pour atteindre 90 % du potentiel (t90 = ln 10 / k ≈ 2,3/k). C’est la façon la plus simple de décrire une courbe BMP et la base de nombreux calculs de dimensionnement, mais il a des limites : il ne tient pas compte de la phase de lag et suppose une dégradation homogène, en une seule étape.

Origine du modèle

Hypothèse : la vitesse de production de méthane est proportionnelle à la quantité de matière encore non dégradée : dB/dt = k·(B₀ − B). La solution donne la fonction exponentielle ci-dessus. Dans les modèles de procédé (par ex. ADM1), l’hydrolyse des polymères est également décrite par une cinétique du premier ordre ; la constante k issue du BMP est donc une « constante d’hydrolyse » approximative dans les conditions données du test.

Interprétation des paramètres

Les valeurs de k pour les substrats dans la littérature couvrent un large intervalle (de fractions à plusieurs 1/j) et dépendent de la méthode ; il ne faut pas les considérer comme des constantes de matériau.

Ajustement

1. Régression non linéaire (moindres carrés) : ajustement de B₀ et k aux points (t, B) — à privilégier.
2. Linéarisation : ln(1 − B/B₀) = −k·t ; exige un B₀ connu et est sensible aux erreurs lorsque B est proche de B₀.
3. Évaluation de l’ajustement : R², RMSE, analyse des résidus (les résidus ont-ils une forme systématique).
4. Incertitude des paramètres : intervalles de confiance issus de l’ajustement ou du bootstrap.

Exemple de calcul

Supposons B₀ = 320 NmL/g de MOS et k = 0,10 d⁻¹.

• t = 5 j : B = 320·(1 − e^(−0,5)) = 320 × 0,3935 ≈ 126 NmL/g,
• t = 10 j : B = 320·(1 − e^(−1)) = 320 × 0,632 ≈ 202,
• t = 20 j : B = 320·(1 − e^(−2)) = 320 × 0,8647 ≈ 277,
• t90 = 2,303 / 0,10 ≈ 23 j.

Pour k = 0,05 d⁻¹, t90 ≈ 46 j. Si le temps de séjour hydraulique (TSH) du digesteur est de 40 jours, dans le premier cas, avec un temps de séjour rigide (comme dans un test en batch), env. 98 % du potentiel seront utilisés, dans le second — env. 86 % (1 − e^(−2) ≈ 0,865). C’est une simplification, car un digesteur à mélange complet présente une distribution des temps de séjour, et non un temps rigide ; dans un digesteur CSTR, la fraction dégradée est k·TSH/(1 + k·TSH) : pour k = 0,10 et un TSH de 40 j, cela donne 4/5 = 80 %, pour k = 0,05 – 2/3 ≈ 67 %. L’écart entre l’approche « en flacon unique » et le CSTR est important et explique une partie des divergences entre BMP et échelle technique.

Limites du modèle

• Phase de lag. Lorsque la production est retardée au début, le modèle du premier ordre s’ajuste mal ; on utilise un modèle avec retard ou celui de Gompertz.
• Courbes à deux étapes. Les substrats contenant une fraction facilement et une fraction difficilement dégradable (par ex. un ensilage avec fibres) exigent un modèle à deux fractions (somme de deux exponentielles).
• Inhibition. Les effets de toxicité et d’acidification ne découlent pas de la cinétique du premier ordre.
• Dépendance aux conditions. k varie avec l’I/S, la température, la taille des particules.

Bonnes pratiques

1. Ajustez le modèle à l’ensemble de la courbe, et pas seulement au point final.
2. Indiquez B₀, k et leurs incertitudes ainsi que la méthode d’ajustement.
3. Ne comparez pas les k entre séries sans vérifier que l’inoculum et l’I/S étaient comparables.
4. Utilisez k pour un choix approximatif du TSH, mais confirmez par des tests en continu.

Application à la planification

Pour des substrats de B₀ voisin, un k plus rapide signifie la possibilité d’un TSH plus court sans perte de rendement. C’est un argument en faveur du prétraitement (broyage, hydrolyse), qui augmente généralement k, et pas nécessairement B₀. Si, après traitement, k double pour un même B₀, un TSH plus court à rendement égal peut libérer du volume de digesteur ou accroître le débit traité.

Liste de contrôle rapide

• [ ] Les données cumulées, et non journalières, sont l’entrée de l’ajustement.
• [ ] Tous les points ont été pris en compte, y compris ceux du début du test.
• [ ] On a vérifié l’absence de phase de lag.
• [ ] B₀, k, les incertitudes et R² sont indiqués.
• [ ] L’outil utilisé est consigné.