El presente trabajo está encaminado a brindar una metodología para la obtención de los modelos cinéticos de las principales reacciones complejas en la etapa de añejamiento de aguardientes. Para ello se analizaron las curvas de añejamiento de aguardiente, siguiendo el comportamiento de las variables principales del proceso. Se aplicó el método integral de análisis de datos cinéticos, obteniendo los modelos cinéticos conformados por un sistema de 4 ecuaciones diferenciales ordinarias no lineales, este sistema se resolvió mediante el software MatLab R2010a y los resultados se validaron estadísticamente.