O objetivo deste trabalho foi identificar entre os modelos Brody, Von Bertalanffy, Gompertz, France, Logístico, Logístico modificado e Logístico bicompartimental, aquele que apresenta maior qualidade de ajuste à curva de produção cumulativa de gases em silagens de girassol e milho. Os critérios adotados foram: coeficiente de determinação, quadrado médio do resíduo, análise gráfica das curvas observadas e estimadas, análise gráfica de dispersão dos resíduos estudentizados, erro percentual médio, eficiência relativa e número de iterações para atingir a convergência. Os modelos Brody, France e Logístico bicompartimental apresentaram os maiores valores de coeficiente de determinação em ambos os substratos, e a diferença entre eles pode ser considerada desprezível. Estes modelos apresentaram, também, os menores valores de quadrado médio do resíduo em silagens de girassol, e a diferença entre eles foi considerada desprezível. Os modelos Brody e France apresentaram menor quadrado médio do resíduo em silagens de milho. Todos os modelos apresentaram dispersão positiva dos resíduos em ambos os substratos após 144 horas de incubação. O modelo Brody apresentou menor erro percentual médio e número de iterações em ambos os substratos. Os modelos Logístico bicompartimental e France apresentaram maior eficiência relativa, respectivamente, em silagens de girassol e milho. Assim, o modelo Logístico bicompartimental apresenta maior qualidade de ajuste à curva de produção de gases em silagens de girassol e milho.