Next Article in Journal
Ionic Storage Materials for Anodic Discoloration in Electrochromic Devices
Previous Article in Journal
Techno-Economic Analysis of a Highly Renewable and Electrified District Heating Network Operating in the Balancing Markets
Previous Article in Special Issue
Repurposing End-of-Life Coal Mines with Business Models Based on Renewable Energy and Circular Economy Technologies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Criteria Optimization of Energy and Water Consumption in Fruit- and Vegetable-Processing Plants in Poland

1
Institute of Mechanical Engineering, Warsaw University of Life Sciences, 02-787 Warszawa, Poland
2
National Fund for Environmental Protection and Water Management, 02-673 Warszawa, Poland
3
Faculty of Electrical and Computer Engineering, Cracow University of Technology, 31-155 Kraków, Poland
*
Author to whom correspondence should be addressed.
Submission received: 18 October 2023 / Revised: 4 December 2023 / Accepted: 14 December 2023 / Published: 17 December 2023

Abstract

:
Fruit and vegetable processing comes 6th in terms of energy consumption in the agri-food industry. At the same time, 88.4% of the industry’s final energy consumption structure is thermal energy, which depends heavily on electricity consumption. In addition, fruit and vegetable processing has a significant impact on the environment due to consumption of significant amounts of water. Reducing these three indicators simultaneously would increase the efficiency of the process while improving environmental protection. This paper proposes neural models of thermal energy, electricity and water consumption for selected major fruit- and vegetable-processing plants in Poland. These models were the basis for formulating a multi-criteria optimization task. Optimization of thermal energy, electricity and water consumption was carried out using genetic algorithms. The optimization results in the sense of Pareto can be the basis for the use of sustainable technology in selected fruit- and vegetable-processing plants.

1. Introduction

In Poland, fruit and vegetable processing comes 6th in terms of energy consumption in the agri-food industry, after the dairy, meat, fish, sugar and bakery sectors. At the same time, thermal energy accounts for 88.4% of the industry’s final energy consumption structure. Multi-directional processing of raw materials and their quality, variety of processes and operations, changeable operating conditions and non-simultaneous operation of equipment and seasonality of production contribute to high variability of consumption, especially of thermal energy [1,2,3,4,5]. Available publications [6,7,8,9,10] present selected cause–effect relationships in energy use but do not fully explain which factors influence energy and water consumption. Specific energy consumption rates for the same product may vary in individual cases. This is also due to differences in the technical equipment of the plants as well as the variety of methods used to establish these indicators. Detailed results of research on the determinants of fruit- and vegetable-processing energy intensity and water consumption are also discussed in publications by [11,12,13,14,15]. Gil et al. [16] also took into account aspects of production hygiene and related water consumption. The research is also justified due to the impact of fruit and vegetable industry plants on the environment [7,17] and the implementation of cleaner production principles [18]. The problem affects in particular African countries, where the demand for water, energy and food resources is increasing and communities have limited availability and affordability in these areas [19]. Optimization of water consumption in fruit and vegetable processing was achieved with the use of genetic algorithms [20]. The results obtained allow for the appropriate selection of the production structure to ensure the lowest water consumption. Attempts to simultaneously optimize the three criteria will therefore achieve even greater results. A publication [8,21] discusses attempts to describe the influence of various factors on the consumption of energy and water in agricultural production. However, their application is limited and they contribute only to a partial clarification of these poorly understood issues. Research in this area was also carried out by [22], and an analysis of energy and water consumption for sustainable paperboard production was undertaken by Man et al. [23]. The literature describes ranges of variation in specific heat consumption rates in agricultural processing companies but does not mention factors that may affect these values.
The aim of the work is to analyze the above issues leading to the construction of models of heat, electricity and water consumption in fruit- and vegetable-processing plants in Poland, which will become the basis for formulating an optimization task. The first part presents the current methodology for examining energy and water consumption in fruit- and vegetable-processing plants, which became the basis for creating a neural model. The complexity of the issue (many variables describing the three output variables, non-linear nature of the relationship) required a preliminary analysis of the data in order to select the most important predictors and the scope of research. A neural energy and water consumption model was used to describe the process. The use of neural networks is justified in such cases, as evidenced by similar studies [24,25]. The neural model was used to formulate the optimization task and then search for the optimal solution using genetic algorithms.
The following steps were carried out in the study.
  • On the basis of empirical studies in fruit- and vegetable-processing plants, a dataset was collected on the basis of which three regression models of heat, electricity and water consumption were developed.
  • The dataset was subjected to preliminary data analysis, excluding incomplete and outlier cases. The resulting 808 data cases were subjected to a non-hierarchical cluster analysis method, determining a set of 604 cases excluding cases with low values of independent variables (in such cases, the process is incomplete due to the limited volume of processed product and, in such a process, the use of equipment is inefficient).
  • The dataset of 604 cases was subjected to an analysis of variance (ANOVA) to isolate the variables (predictors) that have the most significant impact on all three dependent variables (heat, electricity and water consumption).
  • The extracted dataset and predictors were the basis for developing a neural model of the process describing heat, electricity and water consumption together.
  • The neural model of heat, electricity and water consumption was the basis for formulating a multi-criteria optimization task.
  • The solution of the optimization task was performed using genetic algorithms obtaining a set of optimal solutions in the Pareto sense.

2. Analysis of Fruit and Vegetable Processing

The materials and results of the measurements came from 16 fruit- and vegetable-processing plants researched during the summer. Fifty daily periods were analyzed at each site to obtain the necessary datasets. Figure 1 shows a generalized diagram of a fruit- and vegetable-processing plant.
Indicators We, Wc and Ww stand for specific consumption rates for electricity, heat and water, respectively, while Ae, Ac and Aw for the daily consumption of electricity, heat and water, while Z is the throughput of raw materials per day. The determinants of thermal energy consumption at the surveyed facilities can be divided into four groups. Group 1 is the general characteristics of the production facilities studied, described by the total installed capacity of the production facility’s electrical equipment P and the total volume of the facility’s premises V2. Group 2 is the structure of installed electrical capacity described by P1, P2, P3 and P4, which respectively mean the installed capacity of: electrical equipment in the plant boiler room, hydrophore plant and water treatment plant, beverage and juice production lines, equipment used in storage, freezing and air conditioning (including ammonia compressors) and electrical equipment in the administrative and amenity buildings and plant lighting. Group 3 characterizes the structure of daily raw material processing or production and is described by Z1, Z2, Z3, Z4, Z5 and Z6 respectively denoting the daily production volume of: fruit concentrates, vegetable concentrates, beverages, frozen fruit, frozen vegetables and juices. Group 4 is the K2 indicator, which determines the level of technical and technological equipment, organization of production processes and spatial development. Energy losses are described by the Qi index. Other variables adopted in the study that were found to be insignificant were not included. The correlation and strengths of the adopted independent variables on the selected dependent variables (daily consumption of energy carriers Ac and specific consumption rates for thermal energy Wc and water Ww) were established. Previous studies had analyzed cause-and-effect relationships that are a function of factors of low significance or that express trends assumed at the plant design stage. To explain the dependence of y on a number of independent variables (which are actual parameters observed in practice or functions of them), Formula (1) was adopted:
y = b0 + b1x1 + b2x2 + … + bkxk
in which: y—explanatory variable (Ac or Wc), x—explanatory variables (e.g., P1, P2, K2, V1, Z1, Z2, Z3). Application of the resulting empirical formulas when conditions (2) are met:
b1x1 + b2x2 + … +bkxk ≥ b0 and xi ≥ 0 for i = 1....k.
allows the problem in question to be explained to a large extent in the analyzed fruit and vegetable production facilities.

2.1. Thermal Energy Consumption

On the basis of the analyses, equations expressing the influence of the factors covered by the four adopted groups on thermal energy consumption are shown in Table 1. Only those regression equations for which the correlation coefficient R ˃ 0.75 were included. The average specific heat consumption of the analyzed plants, Wc, for the daily period varied considerably and was (0.2–45.7) × 10−3 [GJ/kg].
The lowest Wc rates were found in plants with a high share of refrigeration and freezing in the production technology. It can be seen from Table 1 that 61.8% of the variability in daily thermal energy consumption is attributed to Group 2 factors (installed capacity in the plant boiler house, hydrophore plant and water treatment plant and juice and beverage line). This is because the installed capacity of the electrical equipment in the boiler plant covered by P1 is correlated with the size of the thermal equipment used to generate and use heat. Values Z1 and Z2 (production of fruit and vegetable concentrates) and Z3 (production of beverages) included in Group 3 were responsible for 59.6% of the impact on daily thermal energy consumption. In this context, based on the work of Singh [26], we can analyze the example of a technological specific heat consumption rate (WT) for the production of canned vegetables, which is 5.187 × 10−3 [GJ/kg]. Within the structure of this indicator, nearly 50% is accounted for by heat in steam and hot water consumed in the process. It should be noted that about 23% of the heat energy was losses, indicating the need for closed-loop circuits and waste heat recovery (Qn in Figure 1). This is confirmed by detailed studies by Cuéllar and Webber [1]. Research on the energy intensity of production was also conducted by Gasparino et al. [27] and Sogut et al. [28]. Alvarez et al. [29] presented opportunities for implementing innovations in energy-efficient juice production. The production of frozen citrus concentrates required a thermal energy input of 8.234 × 10−3 [GJ/kg]. In this case, steam and hot water heat demand and direct heating oil consumption accounted for 83%. It should be noted that this rate was more than double the WT technological specific electricity consumption rate. The application of the four groups of factors provides information on the combined influence of technical and technological factors, degree of mechanization of production operations, organizational and production factors and spatial development on thermal management.
The results obtained confirm the trend that the cubic capacity of both production and non-production rooms has a strong influence on heat consumption Ac and the level of specific thermal energy consumption Wc. The K2 indicator in the equation (Table 1) is a function of the total volume of the plant and the daily throughput of raw materials. Due to the high degree of correlation (R = 0.916), it can be concluded that the equation obtained expresses the cause of the changes in specific thermal energy consumption in fruit- and vegetable-processing plants. Indeed, studies have shown that more than 84% of the variability in specific heat consumption is attributed to the K2 indicator. In practice, this formula has been shown to be useful when K2 is contained within the limits of (4000–30,000) × 10−3 [m3/kg] under established and fault-free operating conditions. The work [30] and WS Atkins International (1998) [31] show that in fruit- and vegetable-processing plants, the average specific heat energy consumption Wc for an annual period was 8.33 × 10−3 [GJ/kg] of processed raw materials, with a maximum value of 32.40 × 10−3 [GJ/kg] of raw material. Research [32] shows, for example, the technological specific thermal energy consumption rate WT for apple concentrate was 8.93 × 10−3 [GJ/kg] of product. For tomato concentrate, the corresponding value was 4.69 × 10−3 [GJ/kg]. The specific electricity consumption rate We was in the range of (22–1450) × 10−3 [kWh/kg].

2.2. Electric Energy Usage

On the basis of the analyses, empirical formulas expressing the influence of the factors covered by the four adopted groups on electricity consumption are shown in Table 2. When Group 1 factors were used, only the impact of the total installed capacity was significant. More than 63% of the variability in daily electricity consumption Ae was attributed to the mentioned factor P. Group 2 factors are used to study the impact of installed capacity structure on electricity consumption. These factors show that 70.8% of the variation in daily electricity consumption is attributed to the installed capacity of ammonia compressors and equipment used in storage, freezing and air conditioning and, to a lesser extent, to consumers in administrative and amenities buildings. It is a starting point to clarifying the poorly recognized issues of electricity management at these plants. The two formulas obtained should be analyzed together due to the physical nature of the independent variables. Values Z4 and Z5 (production of frozen fruit and vegetables), Z1 (production of fruit concentrates) and Z6 (production of drinking juices) included in Group 3 factors were responsible for 63% of the impact on daily electricity consumption. In this context, based on the work of Singh [26], an example can be given of the technological specific consumption rate WT for the production of canned vegetables, which is 200 [kJ/kg]. The production of frozen citrus concentrates required an electricity input of 4047 [kJ/kg]. The application of the four groups of factors provides information on the combined influence of technical and technological factors, degree of mechanization of production operations and organizational and production factors on electricity consumption.
The research shows that, again, the installed capacity P is significant, together with the utilization rate of the production facilities, as expressed by the Km indicator. This indicator is also dependent on the volume of daily throughput of raw materials. Due to the vast range of variability of the Km indicator, it can be considered that the resulting formula expresses the reason for changes in the specific consumption of electricity. Indeed, studies have shown that more than 94% of the variability in specific electricity consumption is attributed to the Km indicator. In practice, this formula has been shown to be useful when Km is less than 400 [kW/Mg], i.e., under conditions of established and fault-free operation and when the plant has more than 300 employees in the production area. Maximum throughput means the minimum Km indicator. The work of Kubicki (1998) [29] and WS Atkins International (1998) [30] shows that in the fruit and vegetable industry, the average specific electricity consumption for an annual period was 720 × 10−3 [kWh/kg] of raw materials processed and at some plants producing apple concentrate this ratio may be three times lower. The sources mentioned also state that refrigeration is the most energy intensive, and that the plants analyzed did not pay attention to the need to minimize electricity consumption. The issue of manufacturing innovations leading to a reduction in electricity consumption was also addressed in publications [33,34,35,36].

2.3. Water Consumption

Group 2 and 3 factors did not have a significant impact on the explanation of specific water consumption. However, factors were responsible for the impact on daily water consumption Aw in the range of 47.6–54.3% (Table 3). Using Group 4 factors, it was found that 84.3% of the variation in specific water consumption Ww was explained by the influence of indicator K2 (total plant volume per 1000 [kg] of raw material processed per day). There are significant ranges for the K2 indicator. The first of these ranges is around approximately 30 [m3/kg]—to which the lower values of the specific water consumption Ww correspond. This refers to the operation of plants using near nominal production lines. The second range of variability of K2 above 30 [m3/kg] is distinguished by increased values of specific water consumption rates Ww. This may concern plants at times of underutilized capacity. This phenomenon occurs, for example, during the start-up phase of production lines, in the event of a reduction in the supply of raw materials or semi-finished products or in the event of a breakdown. It is also due to the seasonality of production of the industry’s plants and the diversion of highly contaminated raw materials for processing.
Observations at the plants analyzed showed that significant water saving opportunities are associated with increasing the use of condensate (from water obtained from the concentration of fruit juices). Cuéllar and Webber [1] described the potential for reducing water consumption in the production of canned meat and vegetables from 15 × 10−3 to 7.5 × 10−3 [m3/kg]. An overview of the possibilities for reducing water consumption when blanching vegetables was described by Derden et al. [12]. It should be added that the results of the modeling of water consumption in fruit- and vegetable-processing plants by Trajer et al. [20] showed a definite impact of the production structure on the consumption.

3. Data Analysis

Assuming that the K2 indicator is out of range of application, it can be considered that the most important factors influencing the variability of thermal energy and electricity consumption are: P, (x1)—the installed capacity of the plant and the production structure for selected products, Z4 and Z5, (x2)—production of frozen fruit and vegetables, Z1 and Z2, (x3)—production of fruit and vegetable concentrates, Z3 and Z6, (x4)—production of juices and beverages, Z7, (x6)—production of processed fruit and vegetables and Z8, (x5)—other products. The designations in brackets were adopted for further analysis. Similar factors are indicated by literature data from [6,7,8]. For the construction of the neural model of thermal energy, electricity and water consumption, the abovementioned independent variables were used, while the dependent variable Wc was the rate of specific heat consumption in the plant [GJ/kg], We—the rate of specific electricity consumption [kWh/kg] and Ww—the rate of specific water consumption in the plant [m3/kg].
The dataset was subjected to a preliminary data analysis omitting incomplete and outlier cases. The resulting 808 data instances were subjected to a non-hierarchical cluster analysis method using the expectation maximization (EM) algorithm [37], in order to identify similar groups of data characterizing the processing (Appendix A, Figure A1). Cluster 1, with 604 observations, refers to the cases with the largest values of the independent variables (except x6—fruit and vegetable processing) and balanced values of the dependent variables Ww, Wc and We, signifying a process close to the full range of processing capabilities of the given plant. Cluster 2, with 145 observations, refers to cases with low values of the independent variables (except x6—fruit and vegetable preparations) and high values of the dependent variables Ww, Wc and We. This is an incomplete process due to the limited volume of the processed product. In such a process, the use of equipment is inefficient, hence the high energy and water consumption. The analysis of variance (ANOVA) showed that the variable (x6) has no effect on the dependent variables, the p-value being 0.840459 (Appendix A, Table A1). This variable was therefore omitted from the processing optimization. It is expedient to try to model and optimize the process described in cluster 1 in terms of water, heat and electricity consumption.

4. Modelling of ANN Architecture

Artificial neural networks (ANNs) were used to model the electricity consumption rate We, the water consumption rate Ww and the heat consumption rate Wc. The task of ANN was to map five input decision variables: total power (x1), frozen products (x2), concentrates (x3), juices and beverages (x4) and other products (x5) for three output variables: We, Ww and Wc to obtain the smallest mean squared error (MSE) and the highest correlation coefficient R. Input and output parameter values were normalized from 0 to 1 (dividing by their maximum values: 14,237 × 103, 282 × 103, 772.98 × 103, 773.312 × 103, 37 × 103 for the decision variables x1, x2, x3, x4, x5, respectively, and 1413.23 × 10−3, 88.87 × 10−3, 17.59 × 10−3 for the indicators We, Ww and Wc, respectively. The input data (474 cases) were randomly divided into sets of learning (80%), testing (10%) and validation (10%) cases. The Neural Networks Toolbox R2018a [38] located in Matlab was used for learning data. The learning algorithm for the artificial neural networks was the Levenberg–Marquardt algorithm. In order to find the best relationships between the input and output parameters, different activation functions and the number of neurons in the hidden layers were tested (Appendix A, Table A2). Finally, for optimization, an ANN architecture multi-layer perceptron (MLP) 5:14:3 (Appendix A, Figure A3) was selected with five neurons in the input layer (x1, x2, x3, x4 and x5), fourteen neurons in the hidden layer and three neurons in the output layer (We, Ww and Wc) with a logarithmic–sigmoidal activation function for the hidden and output layers (ID 6 in Appendix A, Table A2). The highest correlation coefficient R was 0.92 and the lowest mean squared error (MSE) 0.00489. The best fit was obtained after 48 epochs, for which the smallest MSE was 0.0048939 (Appendix A, Figure A2a). The correlation coefficients for the learning, validation and test data were, respectively: 0.92141, 0.92363 and 0.91156 (Appendix A, Figure A2b). The analysis of the sensitivity of the neural model shows that the production of juices and beverages and the power of devices installed in the plant have a more than twice greater impact on energy and water consumption than other independent variables.

5. Multi-Criteria Optimization of Fruit and Vegetable Processing

Optimization algorithms using linear and non-linear programming sometimes have difficulty finding global optima or, in the case of multi-objective optimization (MOO), a Pareto front. In multi-objective optimization, the definition of solution quality is much more complex than in single-objective optimization (SOO) problems. The main challenges in the MOO environment are: getting closer to the Pareto-optimal front and keeping the solution set as diverse as possible. The first task ensures that the resulting set of solutions is close to the optimum, while the second task ensures that a wide range of compromise solutions are obtained. Today, many engineering multi-objective optimization problems are solved using genetic algorithms (GAs). GAs are stochastic optimization methods, which are inspired by natural evolution [39]. Crossover and mutation are the key operators of the genetic algorithm. Crossover involves the random selection of genes from the chromosomes of a pair of parents. The probability of crossover is usually taken from a range of 0.5 to 0.8. A mutation parameter converts a random gene in a chromosome from 0 to 1 or vice versa. Mutation prevents premature convergence of the algorithm and the loss of valuable genetic information from the population. The probability of mutation is assumed to be between 0.005 and 0.50.
The objective function plays a main role in the genetic algorithm steps. The function should be well formulated because the main genetic operators perform their tasks based on the evaluation of the objective function. The function should be well formulated. Optimization of a multi-criteria genetic algorithm involves simultaneously minimizing or maximizing multiple objective functions (quality criterion) using constraints [40,41]. The basic steps of a multi-criteria genetic algorithm are described in the work [42,43]. In the present work, the objective functions We, Wc and Ww (Appendix A, Formulas (A1)–(A17)) were simultaneously minimized according to the constraints imposed on the decision variables: x1, x2, x3, x4 and x5 (see limitations (3)).
min x = min   W e = [ k W h / k g ]   min   W w = [ m 3 / k g ] min   W c = [ G J / k g ] 412.5   x 1 14,237 [ k W ]   0   x 2 282 × 1000   [ kg / day ]   0   x 3 773 × 1000   [ kg / day ]   0   x 4 312 × 1000   [ kg / day ]   0   x 5 37 × 1000   [ kg / day ]  
Multi-objective optimization was performed using a non-dominated sorting genetic algorithm (NSGA II), implemented in the Global Optimization Toolbox in Matlab R2018a. The following genetic parameters were used for optimization: crossover function was indirect, crossover probability of 0.8, migration was forward, mutation function was adapted, mutation probability of 0.15, number of generations of 300, Pareto front population fraction of 0.5, population size of 80 and selection function was tournament.

5.1. Pareto Solutions Using a Multi-Criteria Genetic Algorithm

Table 4 shows the sixteen potential solutions (ID1–ID16) of the optimal set in the Pareto sense. It can be seen from Table 4 that daily electricity consumption is inversely proportional to daily water and heat consumption. Minimum solutions for the indicators We, Ww and Wc were found for parameters in the following ranges: from 3433 to 4750 [kW] for total power, from 151 × 103 to 192 × 103 [kg/day] for frozen products, from 191 × 103 to 231 × 103 [kg/day] for concentrates, from 173 × 103 to 192 × 103 [kg/day] for 282 × 103 [kg/day] for juices and beverages and from 20 × 103 to 28 × 103 [kg/day] for other products. Figure 2 shows the Pareto curves created using a multi-criteria genetic algorithm.
The set of solutions consists of sixteen points forming a Pareto curve, whose boundaries are defined by the extreme points ID1 and ID16 (Table 4). Point ID16 in Figure 2 is the highest point on the Pareto curve with the lowest electricity consumption We = 1.274 × 103 [kWh/kg], the highest heat consumption Wc = 0.036 × 10−3 [GJ/kg] and the highest water consumption Ww = 2.579 × 10−3 [m3/kg] (Table 4. ID = 16). For point ID16, the total installed capacity of the electrical equipment is x1 = 3514 [kW], frozen products are x2 = 151 × 103 [kg/day], concentrates x3 = 199 × 103 [kg/day], juices and beverages x4 = 173 × 103 [kg/day] and other products x5 = 28 × 103 [kg/day]. The lowest point on the Pareto curve is ID1. For point ID = 1, the highest electricity consumption is We = 3.693 × 10−3 [kWh/kg], while the lowest water and heat consumption are, respectively: Ww = 0.567 × 10−3 [m3/kg] and Wc = 0.002 × 10−3 [GJ/kg]. For point ID = 1, the genetic algorithm found the following solutions: total power (x1 = 3470 kW), frozen products (x2 = 192 × 103 [kg/day]), concentrates (x3 = 204 × 103 [kg/day]), juices and beverages (x4 = 282 × 103 [kg/day]) and other products (x5 = 28 × 103 [kg/day]). It can be seen from Figure 2 that electricity consumption is inversely proportional to water and heat consumption.

5.2. The Optimization Results

The points located on the curves (Figure 2a–c) are non-dominated solutions, since an improvement in the first objective function (We index) causes a simultaneous deterioration in the other objective functions (Ww and Wc indices) and vice versa. No single solution was found in which all the considered functions (minimum We, Ww and Wc) would simultaneously reach optimal values. An improvement in the We indicator (its minimization) causes a deterioration in the other two indicators We and Wc (their maximization). Thus, the We indicator is in constant conflict with the Ww and Wc parameters. In this case, the solution is the set of non-dominated solutions in the Pareto sense (Table 4). The choice of a particular solution depends on the requirements of the chosen plant and the preferences concerning the processing conditions.
The difference in values between the extreme points ID1 and ID6 of the solutions in the case of electricity consumption is three times, in water consumption it is about five times and the greatest difference in heat consumption is as much as eighteen times. There is no attempt in the literature at multi-criteria optimization of energy and water consumption in fruit and vegetable processing, which makes it impossible to compare it with other results. The method described is limited to the technology used. The procedure for a different technology would be analogous in all steps.

6. Conclusions

The knowledge obtained from this study allows us to understand the impact of the most important factors on energy and water consumption in fruit and vegetable processing in Poland and indicates the possibilities of using sustainable technologies in food production. To sum up, it can be said that:
  • The power of the equipment installed in the plant and the production structure of selected products have the greatest impact on the consumption of energy and water in the processing of fruit and vegetables.
  • The optimization results showed that with similar power of the installed devices, it is possible to use up to five times less thermal energy and several times less water consumption, but the processing structure needs to be selected appropriately.
  • Sensitivity analysis of the neural model shows that the production of juices and beverages and the power of equipment installed in the plant have more than twice the impact on energy and water consumption than other independent variables.
  • The results of Pareto optimization can be the basis for the use of sustainable technology in selected fruit- and vegetable-processing plants.
  • Electricity consumption is inversely proportional to water and heat consumption in fruit and vegetable processing.
Further research should focus on analyzing the impact of the technology used on reducing energy and water consumption in processing. Focus should also be put on developing more efficient equipment and technology. The effort also needs to focus on solving the limitation in the availability of products for processing.

Author Contributions

Supervision, Neural modeling, Conceptualization, Writing—review and editing, Project administration, J.T.; Methodology, Software, Optimization, Writing, R.W.; Data processing, Formal analysis, Investigation, B.D.; Original draft, Supervision, Formal analysis, Writing, J.W.; Data analysis, Neural modeling, Writing—review and editing, R.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author. The data is not publicly available due to business confidentiality.

Conflicts of Interest

The authors declare that they have no known competing financial interest or personal relationships that could have appeared to influence the work reported in this paper.

Appendix A

Figure A1. Plot of mean scale variables from EM cluster analysis.
Figure A1. Plot of mean scale variables from EM cluster analysis.
Energies 16 08118 g0a1
Figure A2. (a) Course of the MSE, (b) function and correlation coefficients for individual sets.
Figure A2. (a) Course of the MSE, (b) function and correlation coefficients for individual sets.
Energies 16 08118 g0a2
Figure A3. ANN architecture of energy and water consumption.
Figure A3. ANN architecture of energy and water consumption.
Energies 16 08118 g0a3
Table A1. Analysis of variance (ANOVA) for continuous variables, total number of training cases: 808.
Table A1. Analysis of variance (ANOVA) for continuous variables, total number of training cases: 808.
ANOVA for Continuous Variables, Number of Clusters: 2, Total Number of Training Cases: 808
Outgroup
SS
dfIntragroup
SS
dfFp-Value
x21.943046 × 10511.837805 × 10680685.21550.000000
x34.484888 × 10517.026927 × 10680651.44240.000000
x41.628664 × 10512.308415 × 10680656.86600.000000
x52.088467 × 10318.127049 × 10480620.71240.000006
Ww6.805427 × 10515.072564 × 106806108.13410.000000
Wc2.218745 × 10411.817085 × 10580698.41630.000000
x11.323167 × 10911.400941 × 101080676.12550.000000
We1.440855 × 10811.963493 × 10980659.14610.000000
x63.956767 × 10117.864680 × 1058060.04060.840459
Table A2. Characteristics of the tested artificial neural networks.
Table A2. Characteristics of the tested artificial neural networks.
IDActivate Function in the Hidden LayerNumber of Neurons in the Hidden LayerActivate Function in the Output LayerStatistical Analysis
MSER
1log-sigmoid4 0.01680400.75844
29pureline0.01698100.70719
314 0.17320000.73911
44 0.02718200.63388
59log-sigmoid0.00726540.88576
614 0.00489390.92363
7tan-sigmoid4 0.01601800.71191
89pureline0.01234800.74598
914 0.00624690.88769
104 0.01102100.77659
119log-sigmoid0.00868620.81746
1214 0.00598450.90975
m i n   W e = 1 1 + e ( 23 F 1 0.6 F 2 + 1 F 3 + 3.7 F 4 + 1.3 F 5 1.1 F 6 1.3 F 7 + 0.8 F 8 0.8 F 9 + 0.8 F 10 3.4 F 11 1.3 F 12 + 0.2 F 13 0.5 F 14 + 20.5 )
m i n   W w = ( 1 1 + e ( 1 F 1 0.7 F 2 + 0.8 F 3 0.9 F 4 0.2 F 5 + 1.9 F 6 + 1.9 F 7 0.4 F 8 1.9 F 9 0.09 F 10 + 1.4 F 11 + 0.25 F 12 + 0.7 F 13 2 F 14 1.9 ) )
m i n   W c = ( 1 1 + e ( 9.1 F 1 + 8.3 F 2 1.0 F 3 2.9 F 4 18.2 F 5 + 3.5 F 6 + 2.8 F 7 + 3.8 F 8 28 F 9 + 2.2 F 10 + 6 F 11 + 3.2 F 12 0.6 F 13 + 4.2 F 14 11.4 ) )
F 1 = ( 1 1 + e ( 28 x 1 0.3 x 2 0.07 x 3 + 91 x 4 0.09 x 5 + 27 ) )
F 2 = ( 1 1 + e 5.08 x 1 4.75 x 2 + 10.98 x 3 3.27 x 4 + 4.7 x 5 + 0.5 )
F 3 = ( 1 1 + e ( 23 x 1 14 x 2 + 0.17 x 3 5.8 x 4 22.8 x 5 + 23 ) )
F 4 = ( 1 1 + e ( 11.9 x 1 + 13.4 x 2 5.1 x 3 + 13.2 x 4 + 3.6 x 5 4.3 ) )
F 5 = ( 1 1 + e ( 16.4 x 1 + 6.3 x 2 20.7 x 3 7.7 x 4 43 x 5 + 18 ) )
F 6 = ( 1 1 + e ( 1.4 x 1 + 7 x 2 38.6 x 3 1.5 x 4 + 4 x 5 + 6 ) )
F 7 = ( 1 1 + e ( 15 x 1 24 x 2 + 36 x 3 9.3 x 4 + 4.4 x 5 + 12 ) )
F 8 = ( 1 1 + e ( 86.8 x 1 + 27.7 x 2 60 x 3 + 1.9 x 4 3.4 x 5 64 ) )
F 9 = ( 1 1 + e ( 55.7 x 1 + 23 x 2 + 19 x 3 2.3 x 4 + 8.5 x 5 + 3.5 ) )
F 10 = ( 1 1 + e ( 1.7 x 1 6.4 x 2 31.7 x 3 9.7 x 4 2.1 x 5 + 5.3 ) )
F 11 = ( 1 1 + e ( 9.9 x 1 + 9.5 x 2 + 5 x 3 + 7 x 4 6 x 5 5.4 ) )
F 12 = ( 1 1 + e ( 41 x 1 + 6 x 2 1.8 x 3 + 18 x 4 + 17 x 5 + 12.6 ) )
F 13 = ( 1 1 + e ( 9.8 x 1 + 6.9 x 2 28.6 x 3 13.8 x 4 0.5 x 5 + 9.7 ) )
F 14 = ( 1 1 + e ( 29 x 1 + 10 x 2 + 19.8 x 3 + 6.2 x 4 14 x 5 + 16.2 ) )

References

  1. Cuéllar, A.D.; Webber, M.E. Wasted food, wasted energy: The embedded energy in food waste in the United States. Environ. Sci. Technol. 2010, 44, 6464–6469. [Google Scholar] [CrossRef]
  2. Hitchens, D.M.; Trainor, M.; Clausen, J.; Thankappan, S.; de Marchi, B. Fruit & Vegetable Processing Industry. In Small and Medium Sized Companies in Europe; Environmental Performance, Competitiveness and Management. International EU Case Studies; Springer: Berlin/Heidelberg, Germany; New York, NY, USA, 2003; pp. 117–159. [Google Scholar]
  3. Jongen, W. (Ed.) Fruit and Vegetable Processing: Improving Quality; Elsevier: Amsterdam, The Netherlands, 2002; ISBN 0-8493-1547-1. [Google Scholar]
  4. Lozano, J. Overview of the Fruit Processing Industry. In Fruit Manufacturing; Springer: Boston, MA, USA, 2006; p. 27. [Google Scholar] [CrossRef]
  5. Lu, J.; Vigneault, C.; Charles, M.T.; Raghavan, G.V. Heat treatment application to increase fruit and vegetable quality. Stewart Postharvest Rev. 2007, 3, 4. [Google Scholar]
  6. Classen, J.; Siebenmorgen, T.J.; Griffis, C.L.; Young, J.C. Processing energy requirements for several vegetables. Trans. ASEA 1992, 35, 973–974. [Google Scholar] [CrossRef]
  7. Grzybek, A. Impact of selected technologies on the environment and energy consumption of fruit and vegetable processing (In Polish: Wpływ wybranych technologii na środowisko i energochłonność przetwórstwa owocowo-warzywnego). Rozprawy habilitacyjne nr 12. Inżynieria Rol. 2003, 2. Available online: http://ir.ptir.org/ (accessed on 13 December 2023).
  8. Wojdalski, J.; Dróżdż, B.; Lubach, M. Conditions of heat management in fruit and vegetable processing plants. (In Polish: Uwarunkowania gospodarki cieplnej w zakładach przetwórstwa owocowo-warzywnego). MOTROL Motoryz. Energetyka Rol. 2006, 8A, 286–293. [Google Scholar]
  9. Wojdalski, J.; Drożdż, B.; Lubach, M. Factors influencing electrical energy consumption in fruit and vegetables processing plants. Pol. J. Food Nutr. Sci. 2007, 57, 195–199. [Google Scholar]
  10. Wojdalski, J.; Dróżdż, B.; Lubach, M. Factors influencing energy consumption in fruit and vegetable processing plants. MOTROL Motoryz. I Energetyka Rol. 2007, 7, 277–285. [Google Scholar]
  11. Asgharnejad, H.; Khorshidi Nazloo, E.; Madani Larijani, M.; Hajinajaf, N.; Rashidi, H. Comprehensive review of water management and wastewater treatment in food processing industries in the framework of water-food-environment nexus. Compr. Rev. Food Sci. Food Saf. 2021, 20, 4779–4815. [Google Scholar] [CrossRef]
  12. Derden, A.; Vercaemst, P.; Dijkmans, R. Best available techniques (BAT) for the fruit and vegetable processing industry. Resour. Conserv. Recycl. 2002, 34, 261–271. [Google Scholar] [CrossRef]
  13. Hackett, B.; Chow, R. Energy efficiency opportunities in fresh fruit and vegetable processing/cold storage facilities. In Proceedings of the 2005 ACEEE Summer Study on Energy Efficiency in Industry, Washington, DC, USA, 19 July 2005; American Council for an Energy-Efficient: Washington, DC, USA, 2005. [Google Scholar]
  14. Masanet, E. Energy efficiency improvement and cost saving opportunities for the fruit and vegetable processing industry. In An Energy Star Guide for Energy and Plant Managers; Lawrence Berkeley National Laboratory: Berkeley, CA, USA, 2008; pp. 68–70. Available online: http://escholarship.org/uc/item/8h25n5pr (accessed on 13 December 2023).
  15. Strzelczyk, M.; Steinhoff-Wrześniewska, A.; Rajmund, A. Indicators of water consumption and the quantity of wastewater formed in selected branches of food industry. Pol. J. Chem. Technol. 2010, 12, 6–10. [Google Scholar] [CrossRef]
  16. Gil, M.I.; Selma, M.V.; López-Gálvez, F.; Allende, A. Fresh-cut product sanitation and wash water disinfection: Problems and solutions. Int. J. Food Microbiol. 2009, 134, 37–45. [Google Scholar] [CrossRef]
  17. Wróbel-Jędrzejewska, M.; Stęplewska, U.; Kuleta, P.; Przybysz, Ł.; Polak, E. Carbon footprint of fruit paste technology. Pol. J. Environ. Stud. 2016, 25, 341–349. [Google Scholar] [CrossRef]
  18. Rahim, R.; Ramann, A.A.A. Cleaner production implementation in a fruit juice production plant. J. Clean. Prod. 2015, 101, 215–221. [Google Scholar] [CrossRef]
  19. Muhirwa, F.; Shen, L.; Elshkaki, A.; Chiaka, J.C.; Zhong, S.; Bönecke, E.; Hirwa, H.; Seka, A.M.; Habiyakare, T.; Tuyishimire, A.; et al. Alert in the dynamics of water-energy-food production in African countries from a nexus perspective. Resour. Conserv. Recycl. 2023, 194, 106990. [Google Scholar] [CrossRef]
  20. Trajer, J.; Winiczenko, R.; Dróżdż, B. Analysis of Water Consumption in Fruit and Vegetable Processing Plants with the Use of Artificial Intelligence. Appl. Sci. 2021, 11, 10167. [Google Scholar] [CrossRef]
  21. Wojdalski, J.; Dróżdż, B.; Lubach, M. Factors influencing water consumption in fruit and vegetable processing plants. (In Polish: Czynniki wpływające na zużycie wody w zakładach przetwórstwa owocowo-warzywnego). Postępy Tech. Przetwórstwa Spożywczego 2005, 15, 39–43. [Google Scholar]
  22. Gorbenko, E.; Cheban, A. Optimizations of Technological Process of Processing of Fruit and Vegetable Raw Material on Juice. Motrol 2008, 10b, 220–229. [Google Scholar]
  23. Man, Y.; Han, Y.; Liu, Y.; Lin, R.; Ren, J. Multi-criteria decision making for sustainability assessment of boxboard production: A life cycle perspective considering water consumption, energy consumption, GHG emissions, and internal costs. J. Environ. Manag. 2020, 255, 109860. [Google Scholar] [CrossRef]
  24. Mitroczuk, I.J. Energy Transformation: Challenges and Opportunities—The Polish Case. Environ. Prot. Nat. Res. 2022, 33, 21–34. [Google Scholar] [CrossRef]
  25. Lemes, D.L.; Jacques, M.M.; Sousa, N.B.; Bernardon, D.P.; Sperandio, M.; Silva, J.A.; Chiara, L.M. Estimation of Electrical Energy Consumption in Irrigated Rice Crops in Southern Brazil. Energies 2023, 16, 6742. [Google Scholar] [CrossRef]
  26. Singh, R.P. Energy Accounting of Food Processing Operations (in Energy in Food Processing); Elsevier: Amsterdam, The Netherlands; Oxford, UK; New York, NY, USA; Tokyo, Japan, 1986; Volume 33, p. 36. [Google Scholar]
  27. Gasparino, F.J.; Vitali, A.A.; Viegas, F.C.P.; Rao, M.A. Energy consumption in a concentrated orange juice plant. J. Food Process. Energy 1984, 7, 77–89. [Google Scholar] [CrossRef]
  28. Sogut, Z.; Ilten, N.; Oktay, Z. Energetic and exergetic performance evaluation of the quadruple-effect evaporator unit in tomato paste production. Energy 2010, 9, 3821–3826. [Google Scholar] [CrossRef]
  29. Alvarez, S.; Riera, F.A.; Alvarez, R.; Coca, J.; Cuperus, F.P.; Bouwer, S.T.; Panyor, L. A new integrated membrane process for producing clarified apple juice and apple juice aroma concent. J. Food Eng. 2000, 46, 109–125. [Google Scholar] [CrossRef]
  30. Kubicki, M. (Ed.) Environmental Protection in the Fruit and Vegetable Industry. (In Polish: Ochrona Środowiska w Przemyśle Owocowo-Warzywnym); Wyd. FAPA: Warszawa, Poland, 1998; pp. 30–36, 38–43. [Google Scholar]
  31. WS Atkins International. Ochrona Środowiska w Przemyśle Rolno-Spożywczym; Standardy środowiskowe; Fapa: Warszawa, Poland, 1998; pp. 66–69, 107. [Google Scholar]
  32. Grzybek, A.; Rogulska, M. Material and Energy Inputs and Cumulative Energy Consumption in the Production of Fruit and Vegetable Preserves (In Polish: Nakłady Materiałowo-Energetyczne i Energochłonność Skumulowana Produkcji Przetworów Owocowo-Warzywnych); Zeszyt IBMER: Warszawa, Poland, 1992. [Google Scholar]
  33. McDonald, K.; Sun, D.W. Vacuum cooling technology for the food processing industry: A review. J. Food Eng. 2000, 45, 55–65. [Google Scholar] [CrossRef]
  34. Morris, C.; Brody, A.L.; Wicker, L. Non-thermal food processing/preservation technologies: A review with packaging implications. Packag. Technol. Sci. 2007, 20, 275–286. [Google Scholar] [CrossRef]
  35. Rico, D.; Martin-Diana, A.B.; Barat, J.M.; Barry-Ryan, C. Extending and measuring the quality of fresh-cut fruit and vegetables: A review. Trends Food Sci. Technol. 2007, 18, 373–386. [Google Scholar] [CrossRef]
  36. Szymańska, E.J.; Mroczek, R. Energy Intensity of Food Industry Production in Poland in the Process of Energy Transformation. Energies 2023, 16, 1843. [Google Scholar] [CrossRef]
  37. Witten, I.H.; Frank, E. Data Mining: Practical Machine Learning Tools and Techniques; Morgan Kaufmann: Burlington, NY, USA, 2000. [Google Scholar]
  38. MATLAB 7.6 R2008a. Documentation R.; MathWorks Inc.: Natick, MA, USA, 2008.
  39. Deb, K. Optimization for Engineering Design: Algorithms and Examples; Prentice-Hall: New Delhi, India, 1996. [Google Scholar]
  40. Goldberg, D.E. Genetic Algorithms in Search, Optimization and Machine Learning; Addison-Wesley: Boston, MA, USA, 1989. [Google Scholar]
  41. Sahu, M.; Bhattacharjee, B.; Kaushik, S.C. Thermal design of air conditioned building for tropical climate using admittance method and genetic algorithm. Energy Build 2012, 5, 1–6. [Google Scholar] [CrossRef]
  42. Dyer, J.D.; Hartfield, R.J.; Dozier, G.V.; Burkhalter, J.E. Aerospace design optimization using a steady state real-coded genetic algorithm. Appl. Math Comput. 2012, 218, 4710–4730. [Google Scholar] [CrossRef]
  43. Winiczenko, R.; Górnicki, K.; Kaleta, A.; Martynenko, A.; Janaszek-Mańkowska, M.; Trajer, J. Multi-Objective Optimization of Convective Drying of Apple Cubes. Comput. Electron. Agric. 2018, 145, 341–348. [Google Scholar] [CrossRef]
Figure 1. Scopes of application of indicators of specific consumption of energy carriers in processing plants.
Figure 1. Scopes of application of indicators of specific consumption of energy carriers in processing plants.
Energies 16 08118 g001
Figure 2. The Pareto fronts solutions: (a) Wc–We. (b) Ww-We. (c) Wc–Ww.
Figure 2. The Pareto fronts solutions: (a) Wc–We. (b) Ww-We. (c) Wc–Ww.
Energies 16 08118 g002aEnergies 16 08118 g002b
Table 1. Factors influencing the variability of heat energy consumption in the surveyed plants (SE—standard error).
Table 1. Factors influencing the variability of heat energy consumption in the surveyed plants (SE—standard error).
Group of
Independent Variables
Regression EquationsR2SEIndependent Variables
Designation, DimensionNumerical Range
IIAc = −260.8 + 0.74⋅P2 + 25.15⋅ P 1 0.618269.0P1 [kW]
P2 [kW]
41–1715
25–932
IIIAc = 200.5 + 3.6Z3 + 211.9⋅logZ1 + 443.0⋅logZ2 0.596276.7Z1 [kg]
Z2 [kg]
Z3 [kg]
(0.585–772.980) × 103
(2.920–8.640) × 103
(0.765–191.094) × 103
IVWc = −0.39 + 0.0014K20.8458.9K2 [m3/kg](307–307.692) × 10−3
Table 2. Factors affecting the volatility of electricity consumption in the surveyed plants (SE—standard error).
Table 2. Factors affecting the volatility of electricity consumption in the surveyed plants (SE—standard error).
Group of
Independent Variables
Regression EquationsR2SEIndependent Variables
Designation, DimensionNumerical Range
IAe = 6806.04 + 0.0006⋅P20.63523,360P [kW]413–14,237
IIAe = −45,896.0 + 0.0013⋅ P 3 2 + 29,020.5⋅log P40.70820,920P1 [kW]
P2 [kW]
81–6566
35–3588
IIIAe = 8356.4 + 736.6⋅Z5 + 3468.1⋅ Z 4 +
13,703.4⋅logZ1 + 1.35⋅ Z 6 2
0.63023,560Z1 [kg]
Z4 [kg]
Z5 [kg]
Z6 [kg]
(0.6–773.0) × 103
(0.1–282.0) × 103
(0.7–155.6) × 103
(0.5–312.3) × 103
IVWe = 46.7 + 4.12 Km0.942420Km [kW/kg](9–7929) × 10−3
Table 3. Factors affecting the variability of water consumption in fruit and vegetable industry plants (SE—standard error).
Table 3. Factors affecting the variability of water consumption in fruit and vegetable industry plants (SE—standard error).
Group of
Independent Variables
Regression EquationsR2SEIndependent Variables
Designation,
Dimension
Numerical Range
IIAw = 408.4 + 2.30⋅P10.5431029P1 [kW]41–1715
IIIAw = 2180.0 + 66.6⋅logZ1 + 140.50⋅ Z 3 − 1420.0/Z5 0.4761003Z1 [kg]
Z3 [kg]
Z5 [kg]
(64.0–773) × 103
(11.1–191.1) × 103
(3.8–105.0) × 103
IVWw = 1.4 + 0.005K20.843133.7K2 [m3/kg](563–307,692) × 10−3
Table 4. The Pareto optimal solution.
Table 4. The Pareto optimal solution.
Inputs Outputs
IDx1
[kW]
x2 × 103
[kg/Day]
x3 × 103
[kg/Day]
x4 × 103
[kg/Day]
x5 × 103
[kg/Day]
We × 10−3
[kWh/kg]
Ww × 10−3
[m3/kg]
Wc × 10−3
[GJ/kg]
13470192204282283.6930.5670.002
24422160231180202.4050.7770.005
34750165191193242.9920.6880.005
43928162211182212.0451.1200.008
53541163219192221.8941.2240.009
64334155216177231.7321.3550.013
73786155217186241.6211.5100.015
83627159214196251.5881.6100.017
93859153216177251.4771.7520.020
104113152210180261.4331.8730.024
113600163216177271.4001.9810.025
123742154202179261.3562.1510.028
133433156208186281.3242.2300.029
143650151203175271.2882.3920.033
153510153201177281.2832.4560.034
163514151199173281.2742.5790.036
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Trajer, J.; Winiczenko, R.; Dróżdż, B.; Wojdalski, J.; Sałat, R. Multi-Criteria Optimization of Energy and Water Consumption in Fruit- and Vegetable-Processing Plants in Poland. Energies 2023, 16, 8118. https://0-doi-org.brum.beds.ac.uk/10.3390/en16248118

AMA Style

Trajer J, Winiczenko R, Dróżdż B, Wojdalski J, Sałat R. Multi-Criteria Optimization of Energy and Water Consumption in Fruit- and Vegetable-Processing Plants in Poland. Energies. 2023; 16(24):8118. https://0-doi-org.brum.beds.ac.uk/10.3390/en16248118

Chicago/Turabian Style

Trajer, Jędrzej, Radosław Winiczenko, Bogdan Dróżdż, Janusz Wojdalski, and Robert Sałat. 2023. "Multi-Criteria Optimization of Energy and Water Consumption in Fruit- and Vegetable-Processing Plants in Poland" Energies 16, no. 24: 8118. https://0-doi-org.brum.beds.ac.uk/10.3390/en16248118

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop