Introducción
En el ron fresco, como en la mayoría de las bebidas alcohólicas destiladas, el aroma recuerda a la materia prima utilizada. Este varía cuando se deja en reposo en recipientes de roble por un tiempo determinado, es decir durante el proceso de permanencia de los aguardientes en barriles de roble y al cual se le conoce comúnmente como “envejecimiento” o tiempo de añejamiento. Este tiempo permite que las bebidas sufran una evolución que repercute decisivamente en la mejoría de sus características y propiedades organolépticas(1).
Durante el proceso de añejamiento del ron se producen pérdidas de producto. El añejamiento del ron no cambia ni transforma la bebida, sino que desarrolla y sublima sus cualidades latentes. Es por ello que, en el contexto de excelencia en que compiten estas bebidas, surge el interés por estudiar la merma en el volumen de ron durante su añejamiento en relación a las condiciones ambientales existentes en la solera.
En el mundo el valor medio de pérdidas anuales que se considera para toneles de 180 litros es de 10 % del volumen almacenado. Sin embargo, la literatura señala para los países del Caribe 13 % de pérdidas anuales(2). En un trabajo realizado en las soleras de la ronera de Cárdenas, en los años 1985-1986, se corroboró una pérdida de 10,5 % anual. No se han encontrado reportes recientes del valor de estas pérdidas en Cuba. En Jamaica se reporta un 6% de pérdidas(3). La tecnología existente en las bodegas permitió llevar a cabo un estudio del comportamiento de las mermas durante 13 meses de añejamiento. Se realizaron mediciones del nivel del líquido de los toneles, grado alcohólico, temperatura y humedad. Toda esta memoria almacenada constituye una valiosa fuente de información que puede ser útil en la comprensión del presente y en la predicción del futuro(4).
Se tienen así, cálculos desactualizados y altos porcentajes de mermas en volumen de las bases en solera con desconocimiento de los factores que influyen.
Materiales y métodos
Verificación continuada de procesos
Desde hace varias décadas la validación ha sido uno de los aspectos de mayor escrutinio dentro de las buenas prácticas de fabricación como vía principal para el aseguramiento de los procesos. En el enfoque actual la validación del proceso consiste en: “recolectar y evaluar información del proceso desde su etapa de diseño hasta la producción, que establece evidencia científica de que el mismo es capaz de entregar un producto de calidad de forma consistente”(5).
El ciclo de vida de validación cuenta con tres etapas: diseño del proceso, calificación del proceso y verificación continuada(5). La verificación continuada del proceso es la tercera y última etapa dentro del ciclo de vida de validación. Es una alternativa al enfoque tradicional de validación de procesos, donde una vez concluida la calificación del proceso productivo, se monitorea y evalúa sistemáticamente su mantenimiento del estado de control. En el desarrollo de un nuevo producto o de uno ya existente es necesario dominar el proceso conociendo cómo los materiales, procesos y controles afectan el producto final, ya que se puede identificar los factores críticos de entrada al proceso. La variación en las respuestas del proceso debe ser entendida y controlada. Desde una perspectiva de fabricación, las materias primas y los parámetros del proceso afectan en mayor o menor medida los atributos críticos de calidad.
Para conocer cuáles parámetros del proceso son críticos y cuáles claves, conocer en qué dirección y magnitud estos parámetros afectan los atributos críticos de calidad es fundamental para la excelencia de fabricación. Esta comprensión puede ser utilizada para desarrollar procesos más robustos o añadir medidas de control adecuadas para ajustar los procesos y así obtener el producto objetivo(6).
Procedimiento de cálculo de mermas en toneles de añejamiento de ron
Para la realización de este estudio se instaló un dispositivo de medición de temperatura y humedad relativa ambiente marca HOBO, modelo H08-003-02 con registrador de datos en una solera. Estas variables físicas se midieron ininterrumpidamente durante catorce meses.
Los toneles se hallan almacenados en cinco ramblas que poseen una distribución vertical de cinco niveles cada una; en ellas se observa la identificación numérica de cada tonel.
Se escogieron tres puntos de la solera (centro y extremos) y se midieron los toneles del primer, tercer y quinto nivel, para evaluar las pérdidas en relación a la posición y la altura de los mismos. Se escogieron 45 toneles los cuales tienen almacenados diferentes bases.
El dispositivo de medición se instaló sobre las duelas de uno de los barriles del nivel más elevado de la rambla, en las cercanías del tapón utilizado para el llenado/vaciado de la barrica y del techo de la solera.
Los toneles contienen cinco (cinco) productos con diferentes concentraciones de etanol (formulaciones) y se encuentran distribuidos aleatoriamente en la solera; denominados: B Vi, B Va, B A, Ag, V V(7).
Análisis de Pareto
El diagrama de Pareto es una herramienta estadística que permite organizar por orden de relevancia los problemas o las causas que los generan. La utilización de este diagrama permite que cuando se quiera mejorar un proceso o atender sus problemas se establezcan prioridades y se enfoquen los esfuerzos donde puedan tener mayor impacto.
El viejo adagio dice que el 80 % de los problemas reportados se asocia con el 20 % de las diversas causas esenciales. Concentrando los esfuerzos en corregir este 20 %, se puede tener el mayor impacto inmediato en la calidad del producto. En esta investigación fue elaborado utilizando la herramienta Statgraphics Centurion v17.2.
Análisis de componentes principales (ACP)
El ACP es una técnica que facilita el análisis exploratorio de datos al explicar las fuentes de variabilidad, reducir dimensionalidad, eliminar datos discrepantes y extraer la información importante. En el trabajo se utilizó el software The Unscrambler X v10.4, que es una herramienta integral de análisis de datos para estadística exploratoria, análisis multivariado, clasificación, predicción y diseño de experimentos(4).
Resultados
Determinación de los atributos críticos de calidad
Se clasificaron los atributos críticos de calidad de acuerdo al impacto sobre la calidad del producto final y el riesgo potencial sobre la aceptación del lote. Para ello se empleó el cuadro 1 de asignación de índices, obteniéndose los resultados que se muestran en el cuadro 2(6),(8).
Para la determinación de los parámetros críticos se formó una comisión de seis expertos, integrada por másteres y profesionales en las especialidades de ciencia y tecnología de los procesos biotecnológicos, ingeniería de alimentos y química analítica; y se analizó un conjunto de variables que, según la experiencia de trabajo y el conocimiento del proceso, tenían mayor o menor influencia en los atributos de salida del producto final. Se tuvo en cuenta en este análisis solo la etapa de añejamiento.
Determinación de los parámetros críticos de calidad
En la valoración numérica para cada parámetro de proceso se empleó el cuadro 3 de asignación de índices y con los valores de criticidad asignados a cada atributo crítico de calidad se determinó la suma de los productos por cada parámetro en la fila de la matriz correspondiente y su porcentaje respecto al gran total.
En el cuadro 4 se muestra la matriz de valoración de los resultados obtenidos.
En la figura 1 se muestra el Diagrama de Pareto, representando además la línea ascendente del porcentaje acumulativo para cada parámetro.
En la figura 1 se observan las variables que tienen mayor influencia en los atributos críticos de calidad del producto final. Esto se justifica a partir del criterio de Pareto 80-20%, donde estas variables representan los parámetros que necesitan de un control estricto para así garantizar que los atributos críticos de calidad del ron, se mantengan en las especificaciones establecidas (los que las tengan). Estos son los llamados parámetros críticos del proceso (CPP), y son: tiempo de añejamiento, volumen, concentración de dioxígeno; temperatura y humedad. Las variables: grado alcohólico inicial, quemado e hinchado del tonel son los parámetros claves del proceso (KPP), ya que influirán de cierta forma en el desempeño del proceso, pero no en los atributos críticos de calidad.
En este estudio queda recogida la mayor parte de la información teórica y las experiencias productivas de los especialistas y operarios del área. Para el estudio posterior se decidió excluir la variable concentración de dioxígeno debido a que no se cuenta con la tecnología apropiada para medirla. Además, se agregó al análisis la variable grado alcohólico inicial ya que su valoración se realizó teniendo en consideración que puede ser ajustado para las especificaciones finales del proceso y no por el impacto específico a lo largo de este.
Determinación de ruido, limpieza y selección de los datos que se utilizaron
Se contó con un conjunto de datos integrado por 900 instancias y diez variables, de ellas cinco cualitativas: producto, rambla, fecha y posiciones horizontal y vertical; mientras que las restantes: tiempo de añejamiento, volumen, temperatura, humedad y grado alcohólico inicial tienen carácter cuantitativo. La matriz inicial de datos se redujo a 546 valores ya que presentaba mediciones incompletas con respecto a tres toneles (19, 21 y 159) y al período de tiempo que comprende de febrero a septiembre de 2014.
Análisis estadístico y transformación de los datos
El hecho de que todas las variables se hayan registrado por su propia naturaleza justifica que se hayan escalado por el rango y centrado por la media. El ACP se realizó partiendo de un fichero de datos con las variables producto, rambla, posiciones horizontal y vertical, fecha, volumen, grado alcohólico, temperatura, humedad y tiempo de añejamiento.
El gráfico de la influencia para este ACP se puede observar en la figura 2.
Se puede observar que no existen muestras que constituyan puntos discrepantes peligrosos, los que de existir aparecerían en el espacio representado por el círculo rojo ubicado en la sección superior derecha del gráfico.
En el gráfico anterior se puede observar que las muestras 209 y 210 (enmarcadas por el círculo azul) tienen alto valor residual por lo que constituyen puntos discrepantes pero no son peligrosos y pudieran contener información valiosa. Además, este punto se encuentra dentro de la elipse de Hotelling mostrada más abajo en la figura 3; por tanto, esta muestra no fue eliminada del análisis.
La muestra 390 (enmarcada por el círculo verde) presenta un elevado valor de distanciamiento por lo que constituyen un punto discrepante influyente, pero contiene información valiosa y por tanto no fue eliminada del análisis.
Para complementar el análisis de los puntos discrepantes y descubrir agrupamientos no perceptibles fácilmente se utilizó el gráfico de las puntuaciones con la elipse de Hotelling (figura 3). En este se puede observar cómo la mayoría de los puntos se encuentran dentro de la elipse, evidenciando que todos pertenecen al mismo proceso. Además, se muestra un agrupamiento fuera de la elipse; en la parte inferior. Estas puntuaciones discrepantes corresponden a muestras del producto Ag que presentan la mayor graduación alcohólica de la base de datos, que se encuentran ubicadas en la rambla tres y en la posición horizontal inferior (piso). Esto está directamente relacionado con los cambios abruptos experimentados, en estas instancias de tiempo, por la humedad ambiental, porque en la medida que esta sea mayor o menor se va a favorecer la evaporación de etanol o de agua.
Para analizar el comportamiento de la variable mermas de acuerdo al gráfico de las puntuaciones, se determinaron los valores de la media y la desviación estándar y se estableció que los valores de esta variable menores que el límite inferior eran ideales, los menores que la media eran normales y los mayores que el límite superior eran altos. De esta manera se pudo identificar la zona de mayores valores de mermas, ubicada a la izquierda del gráfico y en color verde, perteneciente al período más caluroso del año en especial a septiembre de 2015; alcanzándose valores extremos de temperatura para el proceso en el intervalo de (34-37) ºC.

Figura 3 Gráfico de las puntuaciones con elipse de Hotelling resultado del ACP al fichero con 10 variables y 546 instancias
Tres componentes principales logran explicar el 88,32 % de la varianza de los datos iniciales. En la figura 4 se observa el gráfico de la varianza explicada. El primer componente explica el 46,43 % de la varianza el segundo el 27,34 % y el tercero el 14,55 %. Si se toma en consideración un mayor número de componentes se lograría explicar el 95,90 % de la varianza, pero el modelo sería más complejo.

Figura 4 Gráfico de la varianza explicada resultado del ACP al fichero con 10 variables y 546 instancias
El gráfico de los pesos correlacionados muestra una mejor visualización del comportamiento de las variables. En la figura 5 se presenta el gráfico de los pesos correlacionados.
Como se puede observar en la figura 5 con excepción del volumen, todas las variables se encuentran entre las dos elipses, por lo que todas son significativas ya que tienen un elevado aporte al primer y segundo componente que logran explicar el 74,64 % de la varianza total inicial. Este resultado justifica la inclusión de todas estas variables para la modelación posterior. La variable volumen al presentar el menor coeficiente variación (7,63 %) justifica que no se incluya entre las dos elipses que explican entre el 50 y 100 % de la varianza total de los datos.

Figura 5 Gráfico de los pesos correlacionados resultado del ACP al fichero con 10 variables y 546 instancias
Además, el gráfico explica el fenómeno físico que se esperaba encontrar, ya que la variable temperatura, ubicada a la izquierda a lo largo del componente principal 1 (PC1) se encuentra relacionada de manera inversa con la humedad; lo cual es lógico ya que este componente es el que mayor influencia tiene en la variable de salida que es el % de mermas en proceso. También se muestra la variable grado alcohólico que influye de manera negativa en el componente principal 2 (PC2), la cual tiene una relación inversa con la variable tiempo de añejamiento.
En el cuadro 5 se encuentran los pesos de cada una de las variables en los tres primeros componentes principales. Se subrayaron además las variables que mayor contribución tienen a cada componente principal; de esta manera el primer CP resultó estar relacionado a las variables medioambientales humedad y temperatura. El segundo CP se vincula al efecto del tiempo de añejamiento en el estudio y su repercusión en el grado alcohólico. Mientras que el tercer CP está mayormente influenciado por la humedad y el volumen, que son de cierta manera aleatorias en las muestras.
Análisis económico
Para la evaluación económica de los beneficios de este trabajo se calculan los porcentajes de pérdidas por mermas y se comparan con el 6 % admisible en la Norma General No. 03-2012 de República Dominicana(9) ya que no existe una norma cubana aprobada.
Los porcentajes de mermas se determinan teniendo en cuenta la diferencia entre los volúmenes final e inicial de los 45 toneles que se estudian. El precio de venta de ron es de 30 $/botella y cada botella contiene 700 mL, quedando determinadas las pérdidas monetarias por concepto de merma en bodegas.
Las mermas en bodega durante los 13 meses de estudio fueron 763,43 L, lo cual equivale a un 10,5 % del volumen inicial existente, que en comparación con lo establecido en la norma dominicana(9), para períodos de envejecimiento entre dos y cinco años (lo cual se ajusta a la mayoría de los datos existentes), supera el 6 % permitido. Además, representa una pérdida de $32 718.
El valor de mermas admisible correspondiente al 6 % establecido en la norma dominicana para la bodega, en este período de tiempo, es 434,72 L, que representa unas pérdidas monetarias de $18 631 y difiere del valor de pérdidas reales $14 087.
Conclusiones
Se determinaron cinco parámetros críticos que influyen en la calidad del proceso de añejamiento mediante un modelo de riesgo basado en matriz de entradas y salidas, estos fueron: humedad, temperatura, tiempo de añejamiento, concentración de dioxígeno disuelto y volumen medido.
Se realizó un análisis exploratorio de los datos que se miden en CERALBE, permitiendo reducir la dimensionalidad de los datos a tres componentes principales que explican el 89,11 % de la varianza, seleccionando las variables temperatura, humedad, tiempo de añejamiento, grado alcohólico y volumen para la modelación posterior.
Se calcularon las mermas en bodega durante los 13 meses de estudio, que resultaron ser 763,43 L, lo cual equivale a un 10,5 % del volumen inicial existente; superando el 6 % permitido en la norma de República Dominicana y que representa una pérdida de $32 718. Se determinó que el 6 % de mermas con respecto al volumen inicial de la solera es 434,72 L, que representa unas pérdidas de $18 631.