Modèle de Gompertz modifié dans les tests BMP : retard, vitesse et potentiel

Le modèle de Gompertz modifié décrit la production cumulée de méthane par une courbe sigmoïde : B(t) = P·exp{−exp[(Rm·e/P)·(λ − t) + 1]}, où P est le potentiel final, Rm — la vitesse maximale de production (NmL/g de MOS/j), λ — le temps de retard (j), e — la base du logarithme népérien. On l’applique lorsque, après le début du test, la production démarre avec retard, puis s’accélère et s’éteint, c’est-à-dire lorsque le modèle du premier ordre ne s’ajuste pas aux données.

Quand la courbe BMP a une forme en S

• l’inoculum nécessite une adaptation au substrat (par ex. graisses, substrats contenant des inhibiteurs),
• le substrat nécessite une hydrolyse avant la production de méthane (matériaux lignocellulosiques),
• un I/S bas provoque une acidification temporaire et un retard,
• l’essai présente deux étapes de dégradation.

Paramètres

Exemple de calcul

Prenons P = 300, Rm = 25, λ = 2 j. Le coefficient Rm·e/P = 25 × 2,718 / 300 ≈ 0,2265 d⁻¹.

(Calcul pour t = 10 : exp(−0,812) ≈ 0,444, exp(−0,444) ≈ 0,641, B = 300 × 0,641 ≈ 192.) La vitesse maximale correspond au point d’inflexion, où B = P/e ≈ 110 NmL/g de MOS ; il découle de l’équation qu’il se situe aux alentours de t = λ + P/(Rm·e) ≈ 2 + 4,4 ≈ 6,4 j.

Ajustement pas à pas

1. Saisissez les données (t, B cumulé) dans un outil de régression non linéaire (tableur avec solveur, R, Python).
2. Fixez les valeurs initiales : P ≈ dernier point, Rm ≈ plus forte pente, λ ≈ point d’intersection de la tangente avec l’axe du temps.
3. Ajustez les paramètres ; vérifiez la convergence.
4. Évaluez le RMSE, les résidus, les intervalles de confiance.
5. Comparez avec le modèle du premier ordre (par ex. critère AIC) ; retenez le plus simple si l’écart est faible.

Pièges

• Excès de paramètres. Trois paramètres pour une courbe courte comportant peu de points donnent des estimations peu fiables.
• P comme BMP. Si le test a été arrêté prématurément, P sera une extrapolation ; indiquez-la avec un commentaire.
• Interprétation de λ. Le retard découle de l’inoculum et de l’I/S, ce n’est pas une caractéristique du matériau.
• Comparaison entre laboratoires. Les paramètres cinétiques présentent une dispersion plus grande que B₀/P.
• Absence de mécanisme. Le modèle est empirique ; il n’explique pas pourquoi un retard est apparu.

Utilisation en pratique

Le modèle de Gompertz est utile pour décrire des courbes dans des publications et des rapports, évaluer l’influence du prétraitement sur le retard, ainsi que pour détecter des problèmes : un λ allongé pour un substrat connu par rapport aux séries précédentes peut indiquer un inoculum plus faible ou un inhibiteur. Il ne doit pas être utilisé pour dimensionner directement un digesteur, car les conditions en continu diffèrent de celles du batch.

Documentation

Dans le rapport, indiquez : la forme de l’équation, le logiciel d’ajustement, les valeurs de P, Rm, λ avec incertitude, les mesures d’ajustement et la plage de données utilisée pour l’ajustement. Si vous utilisez un modèle biphasique ou un autre, nommez-le explicitement. Ainsi, le lecteur peut reproduire le résultat.

Choix du modèle — aide-mémoire

Exemple supplémentaire : effet du prétraitement

Supposons qu’avant traitement un substrat ait P = 280, Rm = 12, λ = 4 j, et après traitement P = 300, Rm = 22, λ = 1,5 j. L’augmentation de P d’env. 7 % est faible, mais la vitesse maximale a presque doublé et le retard a été raccourci de 2,5 j. Pour une installation à TSH court, c’est un changement important : une plus grande part du potentiel sera utilisée dans le digesteur. Si l’on évaluait le traitement uniquement d’après le BMP final, le bénéfice serait sous-estimé.

Informations sur les limites des données

L’ajustement de Gompertz exige un échantillonnage dense autour du point d’inflexion, c’est-à-dire lorsque la production croît le plus vite. Avec des mesures tous les quelques jours, le paramètre Rm sera estimé avec une forte incertitude. Dans un programme d’essais, il convient de prévoir des mesures plus fréquentes pendant les deux premières semaines.

Différences avec d’autres modèles sigmoïdes

Dans la littérature, on utilise aussi les modèles logistique et de Richards, qui ont une forme semblable, mais des paramétrages différents. Changer de modèle modifie les valeurs de λ et Rm ; la comparaison des paramètres n’est donc possible qu’au sein d’un même modèle. Pour choisir, fiez-vous à la qualité de l’ajustement et à la simplicité, et non au nombre de paramètres.