Ensayos BMPExplicación
Modelo de Gompertz modificado en ensayos BMP: retraso, velocidad y potencial
Cuando la curva BMP presenta una fase de retraso y forma sigmoidal, se aplica el modelo de Gompertz modificado. Se explican los parámetros P, Rm y λ, se muestra el cálculo y las trampas comunes.
Otras versiones de idioma: Polski, English, Deutsch, Français, Italiano, Nederlands.
El modelo de Gompertz modificado describe la producción acumulada de metano mediante una curva sigmoidal: B(t) = P·exp{−exp[(Rm·e/P)·(λ − t) + 1]}, donde P es el potencial final, Rm la velocidad máxima de producción (NmL/g SV/d), λ el tiempo de retraso (d) y e la base del logaritmo natural. Se emplea cuando, tras iniciar el ensayo, la producción comienza con retraso, luego se acelera y finalmente se agota, es decir, cuando el modelo de primer orden no se ajusta a los datos.
Cuándo la curva BMP tiene forma de S
- El inóculo requiere adaptación al sustrato (p. ej., grasas, sustratos con inhibidores).
- El sustrato necesita hidrólisis previa a la producción de metano (materiales lignocelulósicos).
- Una baja relación I/S provoca acidificación temporal y retraso.
- En la muestra se observan dos etapas de descomposición.
Parámetros
| Parámetro | Unidad | Significado |
|---|---|---|
| P | NmL CH₄/g SV | potencial asintótico (similar a BMP) |
| Rm | NmL/(g SV·d) | velocidad máxima de producción en el punto de inflexión |
| λ | d | tiempo de retraso, determinado por la tangente en el punto de inflexión |
Ejemplo de cálculo
Supongamos P = 300, Rm = 25, λ = 2 d. El coeficiente Rm·e/P = 25 × 2,718 / 300 ≈ 0,2265 d⁻¹.
| t [d] | exponente = 0,2265·(2 − t) + 1 | B(t) |
|---|---|---|
| 5 | 0,32 | aprox. 76 NmL/g SV |
| 10 | −0,81 | aprox. 192 |
| 20 | −3,08 | aprox. 287 |
(Cálculo para t = 10: exp(−0,812) ≈ 0,444, exp(−0,444) ≈ 0,641, B = 300 × 0,641 ≈ 192.) La velocidad máxima se alcanza en el punto de inflexión, donde B = P/e ≈ 110 NmL/g SV; según la ecuación, esto ocurre alrededor de t = λ + P/(Rm·e) ≈ 2 + 4,4 ≈ 6,4 d.
Ajuste paso a paso
- Introduzca los datos (t, B acumulado) en una herramienta con regresión no lineal (hoja de cálculo con solver, R, Python).
- Establezca valores iniciales: P ≈ último punto, Rm ≈ mayor pendiente, λ ≈ punto de intersección de la tangente con el eje de tiempo.
- Ajuste los parámetros; verifique la convergencia.
- Evalúe el RMSE, los residuos y los intervalos de confianza.
- Compare con el modelo de primer orden (p. ej., criterio AIC); elija el más simple si la diferencia es pequeña.
Trampas comunes
- Exceso de parámetros. Tres parámetros en una curva corta con pocos puntos dan estimaciones poco fiables.
- P como BMP. Si el ensayo finalizó prematuramente, P será una extrapolación; indíquelo con un comentario.
- Interpretación de λ. El retraso depende del inóculo y la relación I/S, no es una característica del material.
- Comparación entre laboratorios. Los parámetros cinéticos muestran mayor dispersión que B₀/P.
- Falta de mecanismo. El modelo es empírico; no explica por qué aparece el retraso.
Utilización práctica
El modelo de Gompertz es útil para describir curvas en publicaciones y informes, evaluar el efecto del pretratamiento sobre el retraso y detectar problemas: un λ prolongado para un sustrato conocido, en comparación con series anteriores, puede indicar un inóculo más débil o la presencia de inhibidores. No debe emplearse directamente para diseñar un digestor, ya que las condiciones continuas difieren de las discontinuas.
Documentación
En el informe, indique: la forma de la ecuación, el programa de ajuste, los valores de P, Rm y λ con su incertidumbre, las medidas de ajuste y el rango de datos utilizados. Si emplea un modelo bifásico u otro, nómbrelo explícitamente. Así, el lector podrá reproducir el resultado.
Selección del modelo — guía rápida
| Observación | Modelo |
|---|---|
| Sin retraso, agotamiento exponencial | de primer orden |
| Retraso + sigmoide | Gompertz modificado |
| Dos fases de descomposición | suma de dos exponenciales |
| Inhibición y velocidad variable | modelo mecanicista (p. ej., basado en ADM1) |
Ejemplo adicional: efecto del pretratamiento
Supongamos que un sustrato antes del pretratamiento tiene P = 280, Rm = 12, λ = 4 d, y tras el pretratamiento, P = 300, Rm = 22, λ = 1,5 d. El aumento de P en un 7 % es moderado, pero la velocidad máxima casi se duplica y el retraso se reduce en 2,5 d. Para una instalación con un TRH corto, este cambio es significativo: se aprovechará una mayor parte del potencial en el digestor. Si se evaluara el pretratamiento solo por el BMP final, el beneficio estaría subestimado.
Limitaciones de los datos
El ajuste del modelo de Gompertz requiere un muestreo denso en la zona del punto de inflexión, es decir, cuando la producción crece más rápidamente. Con mediciones cada varios días, el parámetro Rm se estimará con gran incertidumbre. En el diseño del programa de ensayos, conviene prever mediciones más frecuentes en las dos primeras semanas.
Diferencias respecto a otros modelos sigmoidales
En la literatura también se emplean los modelos logístico y de Richards, que presentan una forma similar pero distintas parametrizaciones. Cambiar de modelo altera los valores de λ y Rm, por lo que solo es posible comparar parámetros dentro de un mismo modelo. A la hora de elegir, guíese por el ajuste y la simplicidad, no por el número de parámetros.
Fuentes
Preguntas y respuestas
¿Cuándo elegir el modelo de Gompertz en lugar del modelo de primer orden?
Cuando la curva presenta una fase de latencia clara y un punto de inflexión (forma de letra S). Cuando la producción comienza de inmediato y decrece de manera exponencial, es suficiente con un modelo de primer orden.
Artículos relacionados
