Modelling and Simulation in Operations and Complex Supply ChainsView this Special Issue
Assessing the Impact of a Vinasse Pilot Plant Scale-Up on the Key Processes of the Ethanol Supply Chain
One of the byproducts generated in the cane sugar production is molasses, which is used for ethanol distillation. However, one of the problems of distilleries is vinasse. Vinasse is highly water pollutant and is dumped untreated in lakes or rivers and damages the environment. The company FALA developed a pilot plant that uses vinasse to produce a type of livestock feed called MD60. In this paper, the impact of the pilot plant’s scale-up in the key processes of the company’s supply chain is analyzed. With the help of a sensitivity analysis, this study finds the values that would allow the company to improve its order fulfillment indicator and to increase profits, assuming an expected demand by the introduction of this new product into the market. The results show that (1) the pilot plant fulfills 32% of the orders, (2) according to the current vinasse storage capacity, it is possible to fulfill up to 77% of the orders by scaling up the pilot plant, (3) to satisfy 100% of the orders, it is necessary to use all the vinasse generated, and (4) the highest profit is reached by processing all the vinasse and by considering the upper sale price.
The sugar industry is one of the most important in the world due to the products and byproducts generated from its production. These become raw material for many supply chains, including ethanol, biofuels, beverages, and pharmaceutical products. According to OECD/FAO , trends in sugar production indicate that, by 2023, 86% of the sugar globally produced will be obtained from sugar cane, and only a small portion will be extracted from sugar beet.
One of the products generated by the sugar industry is molasses, also called black treacle. The thick fluid is similar to honey but darker. According to experts, the darker it is, the more flavor and nutrients it has. For the international market, cane molasses have four main uses: livestock feed, due to their concentration of calcium, chloride, and magnesium; ethanol production; reprocessing for sugar production; and production of yeast, citric acid, and lysine, among others .
When molasses are used for ethanol production, the new byproduct or wastewater called vinasse is generated. Vinasse is a liquid residue remaining from the fermentation and distillation of alcoholic liquors. It has a dark, red color, ranging from crimson to almost black. Similarly, it has a strong odor, moderately acidic pH, temperature close to 85°C, and a biochemical oxygen demand (BOD) between 70 and 80 g/L. Vinasse is also composed of 93% water, 2% inorganic compounds (potassium, calcium, sulfates, chlorides, nitrogen, phosphorus, etc.), and 5% organic compounds that volatilize when heated at 650°C .
Vinasse from the ethanol industry is one of the major environmental pollutants. According to Moraes et al. , for every liter of produced ethanol, between 10 and 18 liters of vinasse are also generated, depending on the distilling process. Moreover, vinasse is usually discharged untreated into bodies of water, and it consequently pollutes surface water and groundwater . Furthermore, the growing ethanol production from cane juices and molasses to use it as automotive fuel  greatly increases the risk of further pollution if technologies that help reduce pollutants are not implemented.
According to CONADESUCA , Mexico produces 190 million liters of ethanol each year, and, considering the aforementioned ethanol-vinasse relationship, this equals 1,900 million liters of vinasse, also per year. Undoubtedly, this is a major issue, not only due to the progressive deterioration of the environment, but also because this is economically ineffective since costs for collection, transport, and disposal of vinasse are constantly increasing . However, if it is properly treated, vinasse can reduce its environmental impact, not to mention that it can become a continuous source of energy and raw materials with high added value .
For all these reasons, the company FALA (this name was given for confidentiality reasons) has designed a pilot plant for vinasse treatment, which, as a result, produces MD60, a type of livestock feed. Unfortunately, the plant does not have enough capacity to process 100% of the generated vinasse. Therefore, this paper analyzes the impact of the company’s pilot plant scale-up in the key processes of its supply chain. Systems Dynamics (SD) is used to carry out this analysis.
The remainder of this paper is organized as follows: Section 2 briefly discusses previous works reporting vinasse use and treatment, and it explains the concept of industrial scale-up and study cases. Then, Section 3 describes the methodology employed to analyze pilot plant scale-up in FALA, while Section 4 presents and discusses the obtained results. Finally, Section 5 includes conclusions and recommendations for future work.
2. Literature Review
2.1. Vinasse Treatment to Reduce Its Environmental Impact
As previously mentioned, the chemicals contained in vinasse make it highly hazardous for the environment. However, if properly treated, it can be used as an energy source, livestock feed, or fertilizer  or for many other purposes. For instance, vinasse has been widely used to generate biogas from anaerobic digestion [4, 9, 10]. Similarly, it has been employed as a fertilizer for sugar cane fields, although the main problem with this is that a portion of vinasse may infiltrate into the soil and damage aquifers [11–16]. In addition, Colin et al.  and Sankaran et al.  used bacteria to retrieve effluents generated from the ethanol industry, while Akram et al.  assessed the potential of vinasse for biofuels production. Another use of the vinasses is the production of livestock feed [5, 19–22].
2.2. Industrial Scale-Up
Industrial scale-up aims at implementing, at a full scale, innovations or new technologies already tested in laboratories or pilot plants in order to improve the performance of industrial processes .
The literature has reported several innovations in industrial processes, such as biogas generation from industrial waste and residues [24–26], and the utilization of vinasse generated from distilling processes . Similarly, other works have assessed the effectiveness of new technologies in improving the characteristics of heavy crude oil , reducing the environmental impact caused by the textile dyeing process , treating wastewater from a polyvinyl chloride production plant , and obtaining value-added products through biorefineries whose raw materials include biomass  and waste and residues [32, 33]. However, despite their success, most of these innovations have been tested only in laboratories or pilot plants. However, in some cases, their economic feasibility is analyzed and suggestions are provided for their implementation on a large scale.
Only few studies have implemented their proposed innovations and new technologies in the industry. Some of these implementations include a brand-new range of products from a steel company , performance improvement of activated sludge reactors , and the use of special equipment to control odors in a wastewater treatment plant . In addition to reporting results from both a pilot plant and a full-scale industrial plant, these studies have highlighted the problems that arise from implementing the technologies at the large scale. Some of these problems include installation issues, failures in operational stability, and excessive economic costs, among others. Unfortunately, none of the works found reports scale-up in pilot plants that use vinasse for the production of livestock feed.
2.3. Use of Systems Dynamics for Industrial Scale-Up
To facilitate industrial scale-up, certain tools must be employed to assess the impact on production processes, as well as the problems, costs, and benefits that would be brought by the implementation of innovations and new technologies into the plant. Simulation, especially System Dynamics (SD), is a tool that enables assessing these variables. SD is based on a holistic system analysis, which is why simulation models designed through it demonstrate the behavior that would emerge in the system if any of its components are modified [37–39]. This is the main reason why SD has been widely employed to analyze and test different alternatives that modify the performance of many industrial sectors, including the cement industry [40–42], the steel industry [43, 44], the oil industry , and the automotive industry, and supply chains [46, 47].
None of the articles found to support this work report the use of SD to simulate the large-scale implementation of innovations and technologies previously tested in a laboratory or in pilot plants. However, one particular study developed a simulation model that analyzed increased vinasse generation as the result of a surge in bioethanol production .
2.4. Complex Supply Chains
The structure of supply chains is complex due to the quantity of suppliers and customers involved [49, 50]. For this reason, different tools and techniques have been developed for their analysis, evaluation, and decision-making; one of these techniques is simulation. Considering its strengths, compared with mathematical programming methods or stochastic models, simulation especially through SD enables the user to find, analyze, and learn about the dynamic behaviors of the system studied, which is helpful in the decision-making process .
SD has been used to analyze the complexity of different types of supply chains [52, 53]. For example, Kumar and Kumar  study the effects of the lack of suppliers in healthcare supply chains for rural communities, caused by their geographical location, Lee and Chung  assess the coordination among the links of perishable food supply chains, since it is critical to guarantee the freshness of theirs products. Sagardi et al.  evaluate the feasibility of the biofuel supply chain in Mexico analyzing the ethanol production and three key variables were identified to develop this sector: the availability of raw materials, the production capacity, and the reduction of carbon dioxide emissions. Peng et al.  simulate the behavior of a postseismic relief supply chain, produced by information delay and damaged road infrastructure. Cannella et al.  tested the effect of inventory record inaccuracy in an information supply chain, finding that it causes instability among the links of the supply chain.
However, the analysis of the state of the art indicates that none of the research works found addresses the impact of industrial scale-up of vinasse treatment over a supply chain and uses SD to assess industrial scale-up. These are two aspects fully analyzed in this work.
To analyze the complexity of FALA’s supply chain, the use of a simulation tool as SD is necessary.
3. Materials and Methods
The company FALA was used to assess the impact of pilot plant scale-up in the key processes of an ethanol supply chain. This company is located in the South of Mexico, and it is a national pioneer in the use of vinasse to produce livestock feed in a pilot plant that was developed by its research department. Figure 1 shows the supply chain of the company.
The supply chain begins with the harvest of sugar cane, which is taken to the mills for sugar production. A residue that is generated in this process is molasses, and although they are commonly used as cattle feed, in this case they serve as raw material for FALA’s production process. That is, the company uses molasses to produce ethanol, which is later used in various products, such as makeup, drugs, beverages, and biofuel.
In its distilling process, FALA generates around 11 liters of vinasse for every liter of ethanol. This type of residue is usually dumped untreated into water bodies, thus damaging the flora and fauna, because it diminishes the oxygen concentration of these water bodies.
FALA perceived this problem and the serious environmental implications. It therefore decided to set up a pilot plant for vinasse treatment, obtaining MD60 as a result, which is a type of livestock feed. However, this pilot plant has the capacity to process only 2% of the daily-generated vinasse. The remainder is stored in holding tanks, but the storage generates additional costs. This paper therefore presents a SD model that assesses different scenarios in order to analyze the impact of the scale-up of the vinasse pilot plant in the key processes of the company’s supply chain. The best alternative will be chosen based on results of two key performance indicators (KPIs). The following section describes the use of SD to develop the evaluation model.
3.1. Causal Diagram
A causal diagram is a graphical tool that helps represent the behavior of a system and the different interactions among its variables. Causal diagrams consist of variables connected through arrows that indicate cause-effect relationships (. These relationships can be either positive or negative. If a relationship is positive, an increase in the cause will produce an increase in the effect, or a decrease in the cause will produce a decrease in the effect. Nevertheless, if the relationship is negative, any increase in the cause will produce a decrease in the effect, or any decrease in the cause will produce an increase in the effect.
Figure 2 depicts the causal diagram used to study industrial scale-up in the vinasse pilot plant. Some of the variables identified include master production schedule, molasses tank, ethanol production, ethanol production capacity, MD60 production capacity, and finished product inventory.
The feedback loops of the causal diagram are described as follows:(i)Loop B1: it has a positive impact on the amount of orders requested from the supplier. However, if the amount of order increases, the molasses inventory of the supplier decreases or vice versa.(ii)Loop B2: if the molasses inventory increases, the ethanol production may also increase, and when the ethanol production increases, the molasses inventory will decrease. Similarly, the ethanol production depends on the production capacity of the plant.(iii)Loop B3: when vinasse inventory increases, MD60 production also rises. An increase in MD60 production causes a decrease in vinasse inventory, which limits MD60 production.(iv)Loop B4: if the variable MD60 production increases, the finished product inventory will also increase. If there are high levels of finished product inventory, MD60 production decreases.
This section describes the equations from the simulation model that correspond to the key processes of the supply chain:
Orders from Supplier. FALA receives daily orders from the supplier (DOS) depending on the condition shown by (1). Thus, during harvest season, the company daily receives between 25 and 30 shipments with 30 tons of molasses. However, when it is not harvest season, the company receives maximum 8 shipments, although sometimes there is no daily supply:FALA uses the distillation method for ethanol production (EP). Therefore, the amount of ethanol obtained per ton of molasses can vary. To represent this variability, the RANDOM function is used, which generates a random value from a range of preestablished values. In this case, and according to historical data, the minimum and maximum values were 235,118 and 292,397 liters of ethanol:Vinasse is the raw material of FALA’s pilot plant, and it depends on EP. Every liter of produced ethanol generates around 11 liters of vinasse . This relationship is expressed by Due to the great amount of vinasse that is generated, FALA reuses it as its raw material. Its supply policy depends on the demand of the livestock feed, MD60. That is, if there is a high demand of MD60, vinasse is treated; otherwise, it is sent to holding tanks as expressed by In this equation, VHT stands for the vinasse that is sent to the holding tanks. It is the difference between the daily-obtained vinasse (DOV) and the vinasse that is temporarily stored in a pit (VP) for its further treatment. The amount of vinasse in is determined by the condition expressed by If the amount of daily-obtained vinasse (DOV) is smaller than the difference between the capacity of the pit (CVP) and the amount stored in (in time ), DOV is then stored in ; otherwise, the amount stored is only what results from this difference.
The pilot plant has a batch production system (BPS) for MD60 production. This system works as expressed in In this equation, PR is the percentage of vinasse retrieved and is the initial load added, while stands for four additional loads with different values which are added throughout the MD60 production process.
As regards industrial scale-up, MD60 production in the full-scale industrial plant (IP) follows the behavior described in where is the continuous flow per hour and represents the number of work hours in a day.
Both plants have PR that behaves as a normal distribution (see (8)) with a mean of 11.58% and a standard deviation of 1.8%:
The Company’s Inventories. Inventory of MD60 (IMD), inventory of vinasse (IV), and inventory of molasses (IM) behave as expressed in (8), (9), and (10), respectively:In these equations, represents the effective sales of MD60, is the amount of molasses to be processed, is the yield of molasses per ton of ethanol, is the amount of orders requested from the molasses supplier, and is the amount of molasses that is ordered. The company’s molasses inventory is considered in the model, since MD60 depends on ethanol production, which in turn could not be achieved with a molasses inventory that is low or inexistent.
Two performance indicators were considered to assess the impact of the pilot plant scale-up in the supply chain of the company FALA: order fulfillment of daily demand (OFDD) and profits gained in 365 days. On one hand, OFDD is calculated by where stands for the amount of days when there is product demand, represents the amount of this demand, and stands for the amount of effective sales.
On the other hand, profits (see (13)) are the difference between total earnings (TE) and total costs (TC):The next step in the methodology of SD is to develop the Forrester diagram considering the above equations.
4. Results and Discussion
The model simulation is run with software Stella® V10.0.4 during 52 weeks. , and hours is the time unit employed. Demand of MD60 is represented by a normal distribution with a mean of 164,000 kg and a standard deviation of 10,000 kg. The warehouse used to store the finished product has a storage capacity for 950,000 kg of MD60. It is assumed that the industrial plant would work 24 hours every day from Monday to Sunday. However, in these 52 weeks, there are two weeks of production stoppage due to the maintenance of the ethanol distillery plant.
Table 1 shows the parameters from both the pilot plant and the full-scale industrial plant. The parameters capacity vinasse pit, vinasse production for day, and capacity storage were chosen according to the technical specifications of the pilot and scale-up plant. Percentage of vinasse retrieve, order fulfillment, and profit were obtained by historical data. As can be observed, OFDD in the pilot plant reaches only 32% and equals 34,600 USD. An increase in both indicators is expected with the implementation of a full-scale industrial plant.
Due to the low percentage of OFDD obtained, a sensitivity analysis was performed to find the values that satisfy 100% of the orders.
4.1. Sensitivity Analysis
A low percentage in OFDD in the pilot plant reflects the company’s inability to satisfy a high demand of MD60. Therefore, to increase the values of the two performance indicators, a sensitivity analysis is carried out by modifying the values of the following variables: capacity vinasse pit (CVP), which represents the full capacity of the pit; continuous flow per hour (), which is obtained by a previous sensitivity analysis and represents the range where the flow can be continuous without stopping the production; initial raw material inventory (IIRM) and maintenance time (MT), which are defined through the experience of the company’s managers; and finally sale price (SP) per kilogram, which is the range of the molasses’ cost fluctuation, the biggest competitor of vinasses; all this is in order to improve performance indicators. Tables 2 and 3 present the results obtained of OFDD and . 91 different scenarios were tested, considering the following values: : 1,500,000. : 47,000, 48,000, and 49,000. : 500,000, 1,000,000, and 1,500,000. : 15, 20, and 25. : 0.11, 0.12, and 0.13.In Table 2, there is a summary of the development of 91 scenarios to achieve 100% of OFDD, but the maximum value of OFDD achieved is 77.90%, obtained with the next combination CVP = 1,500,000 L, = 49,000 L, IIRM = 1,000,000 L, MT = 15 days, and SP = 0.13 USD. In the case of profit, Table 3 shows the result obtained and the major profits achieved are 3,993,269.00 USD with CVP = 1,500,000 L, = 49,000 L, IIRM = 1,000,000 L, MT = 15 days, and SP = 0.13 USD.
Results of the first sensitivity analysis still show inability to fully satisfy the product demand. Thus, a second analysis is carried out, this time by considering that vinasse cannot be completely processed due to capacity constraints, which is why it is sent to holding tanks (Table 4). In this second case, the following values are analyzed: : 1,500,000. : 47,000, 48,000, 49,000. : 500,000, 1,000,000, and 1,500,000. : 15, 20, and 25. : 0.11, 0.12, and 0.13. = 15,000,000.As shown in Table 4, three different scenarios allow obtaining 100% of OFDD; the combinations to achieve are shown in Table 5. However, different scenarios are close to 100%, as 99.99%, 99.90%, and others. This is because of uncertainty and variability in the following variables and processes: demand, production process of ethanol, percentage of recuperation of vinasse, and inventory of raw material.
Table 6 shows the scenario to obtain the higher profits (4,65,820 USD). This scenario has the following combination: CVP: 1,500,000 ; : 49,000 ; IIRM: 1,000,000 ; MT: 15 (days); SP: 0.13 (USD); and holding tank = 15,000,000 . Although this scenario delivers the highest profit, it only achieves 99.92% of OFDD as shown in Table 4. This is due to the variability in the key process of the supply chain, which does not allow satisfying 100% of OFDD.
By using both the vinasse pit and the holding tank, the company could reach 100% of OFFD and obtain the highest profits. However, the use of the holding tank as storage will generate an extra cost that, given the improvement achieved for both OFFDD and P, is worth the investment.
Also, the sensitivity analysis shows a behavioral pattern among the set values, without missing the systems uncertainty and variability. In general, a larger MT decreases the OFDD and . On the other hand, a larger inventory and a larger increase the OFDD and . These relations prove the logic described in the causal diagram.
5. Conclusions and Future Work
This paper evaluates the impact of pilot plant scale-up in the key processes of a supply chain. The pilot plant reuses vinasse generated in ethanol production to produce livestock feed. The major contribution of this research is the identification of those key processes from the company’s supply chain that could be affected in the pilot plant scale-up process. A sensitivity analysis helps to find three scenarios to achieve 100% OFDD. One scenario allowed achieving the higher profits. The key processes of the supply chain of FALA are procurement of vinasses and production of ethanol because the major variability and uncertainty are present here. As future work, it is necessary to assess, by means of sustainability indicators, the environmental benefits of the full-scale industrial plant. Likewise, the relationship cost-benefit must be evaluated by considering time dedicated and investments made to set up the full-scale plant. Finally, although this research considered only two performance indicators, a multicriteria analysis is recommended in order to improve decision-making.
The authors declare that there is no conflict of interests regarding the publication of this paper.
The authors are grateful to National Technological of Mexico for supporting this research. This study was sponsored by the National Council of Science and Technology (CONACYT), as well as by the Public Education Secretary (SEP) through PRODEP.
OECD/FAO, OCDE-FAO Perspectivas Agrícolas 2014, OECD, 2014.
CONADESUCA, Cifras Definitivas Zafra 2014-2015, 2016, http://www.conadesuca.gob.mx/documentos%20de%20interes/CIFRAS%20DEFINITIVAS%20ZAFRA%202014-15.pdf.
SAGARPA, “Estudio de factibilidad para la instalación de una destileria de etanol en Atencingo,” http://www.sagarpa.gob.mx/Delegaciones/puebla/Documents/Evaluaciones%20Externas/2011/Estudio%20Fact%20Int%20Destileria%20Etanol%20Atencingo.pdf.View at: Google Scholar
G. Ortegón, F. Muñoz Arboleda, L. Candela, K. Tamoh, and J. Valdes-Abellan, “Vinasse application to sugar cane fields. Effect on the unsaturated zone and groundwater at Valle del Cauca,” Science of the Total Environment, vol. 539, pp. 410–419, 2016.View at: Google Scholar
M. Andalib, E. Elbeshbishy, N. Mustafa, H. Hafez, G. Nakhla, and J. Zhu, “Performance of an anaerobic fluidized bed bioreactor (AnFBR) for digestion of primary municipal wastewater treatment biosolids and bioethanol thin stillage,” Renewable Energy, vol. 71, pp. 276–285, 2014.View at: Publisher Site | Google Scholar
E. J. Olguín, G. Sánchez-Galván, R. E. González-Portela, and M. López-Vela, “Constructed wetland mesocosms for the treatment of diluted sugarcane molasses stillage from ethanol production using Pontederia sagittata,” Water Research, vol. 42, no. 14, pp. 3659–3666, 2008.View at: Publisher Site | Google Scholar
V. Parnaudeau, N. Condom, R. Oliver, P. Cazevieille, and S. Recous, “Vinasse organic matter quality and mineralization potential, as influenced by raw material, fermentation and concentration processes,” Bioresource Technology, vol. 99, no. 6, pp. 1553–1562, 2008.View at: Publisher Site | Google Scholar
K. Sankaran, M. Premalatha, M. Vijayasekaran, and V. T. Somasundaram, “DEPHY project: distillery wastewater treatment through anaerobic digestion and phycoremediation—a green industrial approach,” Renewable and Sustainable Energy Reviews, vol. 37, pp. 634–643, 2014.View at: Publisher Site | Google Scholar
A. Meneses-Jácome, R. Diaz-Chavez, H. I. Velásquez-Arredondo, D. L. Cárdenas-Chávez, R. Parra, and A. A. Ruiz-Colorado, “Sustainable energy from agro-industrial wastewaters in Latin-America,” Renewable and Sustainable Energy Reviews, vol. 56, pp. 1249–1262, 2016.View at: Publisher Site | Google Scholar
C. Sawatdeenarunat, D. Nguyen, K. C. Surendra et al., “Anaerobic biorefinery: current status, challenges, and opportunities,” Bioresource Technology, vol. 215, pp. 304–313, 2016.View at: Google Scholar
L. Zhang, S. Zhang, Y. Peng, X. Han, and Y. Gan, “Nitrogen removal performance and microbial distribution in pilot- and full-scale integrated fixed-biofilm activated sludge reactors based on nitritation-anammox process,” Bioresource Technology, vol. 196, pp. 448–453, 2015.View at: Publisher Site | Google Scholar
M. G. Cedillo and C. Sánchez, Análisis Dinámico de Sistemas Industriales, Trillas, Barcelona, Spain, 1st edition, 2008.
J. W. Forrester, “Industrial dynamics: a major breakthrough for decision makers,” Harvard Business Review, vol. 36, pp. 37–66, 1958.View at: Google Scholar
J. D. Sterman, Business Dynamics: Systems Thinking and Modeling for a Complex World, Management, Irwin McGraw-Hill, 2000.
S. Umeda and Y. T. Lee, “Design specifications of a generic supply chain simulator,” in Proceedings of the Winter Simulation Conference, R. G. Ingalls, M. D. Rossetti, J. S. Smith, and B. A. Peters, Eds., pp. 1158–1166, December 2004.View at: Google Scholar