1. Introduction
Over the last decade, a constantly growing trend of the popularization of self-service automated stores has been observed. Vending machines have been expanded into fully automated stores, the offer of which is comparable to small, conventional stores. The application of technical solutions in the form of manipulators and belt feeders had a positive impact on the better space management, which, in turn, resulted in a higher number of goods available in stores. Among the advantages of the discussed solution, the most important are contactless sales, ensuring the safety and hygiene of purchases for customers, and 24/7 access. Apart from meeting consumer needs, one of the basic reasons for the popularization of modern automated stores is the reduction in a store’s energy consumption while ensuring a comparable range of products offered. Research into possibilities of reducing greenhouse gases emission is important in terms of the environment and climate protection [
1].
The climate and energy strategy adopted by the European Union assumes that by 2030 greenhouse gases emission will be reduces by 40% (to the level from 1990), while the energetic efficiency will be improved by at least 32.5% [
2,
3].
The construction and service sectors are of key importance for achieving these goals [
1,
2].
The carried-out review of the scientific literature regarding research into modelling the electrical energy consumption in buildings [
2,
3], including schools [
4] or office buildings [
5,
6], confirmed that the topic is up-to-date, but no publications concerning models of automated stores have been found.
Automated stores, due to the relatively limited usable space and the lack of open coolers or store refrigerators intended for storing the items, should allow for a reduction in electrical energy consumption and CO2 emission in comparison to the conventional stores.
Delfin company from Nowa Wies, Poland, which deals with the design and construction of automated stores, performed a preliminary estimation of the annual energy consumption of an automated store and compared the obtained results to the energy consumption measured for a conventional store. The average annual electricity consumption evaluated for a grocery store of an area of 70 m2 amounted to approximately 38.4 MWh, while for an automated store of an area of 9 m2 and a comparable product range, the electricity consumption was approximately 10.1 MWh, i.e., 74% smaller.
However, it should be stressed that the amount of electrical energy consumed by an automated store depends on many factors, such as geographical location and changes in the climatic conditions, including the store’s external temperature. The external temperature of the store determines the amount of energy consumed by the heating, ventilation, and air conditioning (HVAC) or heating and air conditioning (HAC) units, which ensures the appropriate conditions for storing the products. Another important factor is the number of products purchased by users, which translates into the amount of energy supplied to the electrical drives of the feeding system, which is activated when a store’s customer completes their shopping list.
The estimated energy demand for the store is the information of key importance for the potential investors and the manufacturer of the automatic store. Therefore, the research presented in the paper concerns the development of the model for determining the electricity consumption, operating costs and CO2 emission of the automated store designed and developed by Delfin company.
In the designed automated store, three basic systems responsible for the electrical energy consumption were distinguished: the feeding system, the HAC unit and other electrical energy consumption points (marked in the paper as “others”), e.g., the operator panel with the automation system, lighting, outdoor advertising screens, etc.
In many scientific works concerning the development of the energy consumption models, machine learning algorithms are used. The most commonly used algorithms are LSTM (long short-term memory) [
7,
8,
9,
10], MLP (multilayer perceptron) [
11,
12], MLR (multiple linear regression) [
13,
14], ANN (artificial neural network) [
15,
16,
17,
18,
19], SVR (support vector regression) [
20,
21], BPNN (back propagation neural network) [
22], SVM (support vector machine) [
20,
23], XGB (extreme gradient boosting) [
24,
25] and RF (random forest) [
26,
27]. Authors of these scientific papers often point out the possible fields of application of such solutions, e.g., controlling HVAC systems. Practical applications of ML algorithms involve the use of extensive databases, usually containing many types of measured values. As an example, [
5] can be taken, in which the authors compared the MLR, JPR (joinpoint regression), BP (back propagation), RF and JP-MLR (joinpoint−multiple linear regression) models, for which they used eight variables for analysis: five continuous variables (average outdoor air temperature, average relative humidity, daily global radiation intensity, average wind speed and daily temperature amplitude) and three decisive variables (gender of inhabitants, holiday index and sunshine index during the day).
Another approach consists of the usage of thermal models of buildings. In this case, three methodologies are most frequently used [
28,
29]: the temperature response method, e.g., [
30], the finite difference method, e.g., [
31], and the lumped heat capacity method, e.g., [
32]. The application of the lumped heat capacity method allows for the relatively easy determination of a mathematical model of temperature changes [
28,
29,
33], while the physical parameters can be determined based on the parameters of materials used in the construction of the building of interest (such as the heat capacity or thermal resistance of the constructional elements) [
28].
The novelty of the carried-out research results consists of the development of the first known energy consumption model for an automated store with a container structure. It should be stressed that the developed model provides results with sufficient accuracy, which has been validated experimentally for a single location.
In order to determine the automated store’s energy consumption, the temperature data from a nearby meteorological station were sufficient; over the course of the research, it has been proven that it does not have to be measured at the exact location of the store.
Forecasting the value of the automated store’s energy consumption in the various environmental conditions based on the data from the publicly available sources, such as meteorological stations, is convenient from a practical point of view. It enables quick access to a set of information concerning climatic conditions for a given location—the place where the store is installed. However, the information from meteorological stations, e.g., temperature, is often averaged for a given area, with a relatively low number of samples (low sampling frequency). Therefore, for the selected model, it was necessary to carry out the tests to determine the impact of the input data in the form of the measurements from the meteorological station on the accuracy of the calculations and, consequently, on the output values of the model.
In order to calculate the energy consumption of the store’s feeding system, it was necessary to solve the problem of finding a correlation between the energy consumption of the electrical drives of the feeding system and the arrangement of products on the store shelves, which are made up of belt feeders. Depending on the type and quantity of products selected by the customer, the movement trajectory of the feeding system is generated. The total permissible capacity of the goods that the manipulator can transport in one work cycle is limited by the volume of its container.
For a given customer’s shopping list, by knowing the approximate volume of each type of product, the store management system can optimize the length of manipulator trajectory necessary to complete the order. This is beneficial in terms of the order processing time and the amount of electrical energy necessary to power the engines during their operation. However, from the point of view of the created model, it results in the high complexity of the mathematical notation.
Additionally, the set of input parameters of such a model is a shopping list containing the exact number of products of different types, which is difficult to define at the stage of preliminary research regarding the interest of the consumer market. In the carried-out research, the general guidelines for the operation of the store model assumed the calculation of the average value of energy consumption per product, which ensured a simple analysis for various sale variants.
Project assumptions regarding the minimum number of types of data sets necessary to estimate the energy consumption, CO2 emission and costs over the assumed average operating period of one month (approximately 30 days) resulted in the final choice of the model description for the HAC unit based on the lumped method heat capacity, and a linear relationship between the energy consumption coefficient per product for the feeding system.
3. Model of Energy Consumption
In the automated store model (
Figure 4), the three main sources of energy consumption has been taken into account: HAC unit, feeding system and other electrical devices. The input data for the model are the outdoor temperature and the number of products sold. Based on the input data, the model calculates an estimated value of electricity consumption, the mass of CO
2 emitted as a result of generating this energy and the estimated cost of production.
The energy consumption of the considered automated store can be described as follows:
where
ETotal—overall consumption of energy;
EHAC,
EFE,
EOTH—electrical energy consumed by the HAC unit, feeding system and other electronic components;
PHAC(
t),
PFE(
t), and
POTH(
t)—instantaneous power consumed by the HAC unit, feeding system and other electrical components, respectively.
3.1. Model of HAC Unit
The accurate modelling of the air temperature changes inside an automated store is not possible due to the complexity of the occurring physical phenomena. The temperature inside the store is influenced by a number of phenomena, including non-deterministic ones, such as: sunlight, ambient temperature, wind strength, number and type of products in the store assortment, frequency of purchases, frequency of deliveries, etc. Therefore, it is only feasible to develop the approximate mathematical models of temperature changes inside the store and electrical energy consumption. The equations on the basis of which the model proposed in this research has been developed are presented below.
where
Tint,
Tout—temperature inside and outside the store, respectively;
Theat,
Tcool—temperatures of warm and cold air from the cooling and heating unit, respectively;
Mair,
Meq_rest—mass of air inside the store and the equivalent mass of other elements (except air) present inside the store (including food products, equipment elements, feeders);
Mdot—air mass supplied from the HAC unit per unit of time;
cair,
ceq_rest—heat capacity of air at the constant pressure and the equivalent heat capacity of other elements (apart from air) inside the store (including food products, equipment elements, feeders), respectively;
Pheat,
Pcool,
Ploss—heating and cooling power of the HAC unit and the power loss, respectively;
Peq_eq—equivalent power emitted by elements inside the facility (including power supplies, electronic systems and electric drives);
PHAC_Standby,
PHAC_Heat,
PHAC_Cool,
PAC—electrical power consumed by the heating and cooling unit in the standby, heating and cooling states, and the overall electrical power consumed by the heating and cooling unit;
uh(
t),
uc(
t)—time histories of on/off signals for heating and cooling, respectively.
Equations (3) and (4) describe heat flow in the cooling and heating unit working in the heating and cooling modes, respectively. The parameter Mdot represents the mass of air injected by the HAC unit. It was estimated on the basis of the air density and flow provided in the HAC catalogue note and then corrected in a way that allows for the best tuning of the model to the real data. Equation (5) describes thermal power losses resulting from the usage of materials with the equivalent thermal resistance Req in the construction of the store. It was assumed that the equivalent thermal resistance is equal to the thermal resistance of the PUR (polyurethane) foam used for insulation, and its value was read from the manufacturer’s documentation. Equation (6) describes the variability of the temperature inside the store, taking into account the power supplied by the HAC unit, heat losses resulting from the finite thermal resistance of the structural elements and the equivalent power dissipated in the form of heat by the store equipment elements located inside, i.e., power supplies, electronic systems, electrical drives. The latter was estimated based on the measurements of electrical energy consumed by the store and the efficiency of individual components. In order to tune the model, it was also necessary to estimate the mass of air and the mass of other elements inside the store (mainly food products and equipment), as well as the equivalent heat capacity of these elements.
These values were initially estimated on the basis of the average amount of products available inside the automated store, and then adjusted to obtain the best fit of the model to the real object. Equations (7) and (8) make it possible to determine the electrical energy consumption based on the time histories, thereby determining the moments the thermostat was turned on and off by the heating and cooling unit.
A diagram of the electrical energy consumption model developed in the MATLAB/Simulink environment is presented in
Figure 5. The main elements of the developed model are blocks representing a thermostat, a cooling and heating (HAC) unit, and a store container.
The first experiments on the HAC unit model concerned the identification of the model parameters. The model input includes a set of measured outdoor temperature values (
Figure 6,
TOUT Real) obtained from the automated store monitoring system (
Section 2.1).
The measurement time amounted to one day (24 h), and the sampling frequency was equal to 0.033 Hz. A comparison of the recorded temperature inside the store (
Figure 6,
TINT Real) and the temperature obtained at the model output (
Figure 6,
TINT Mod) allowed for the determination of the values of the model parameters, including the equivalent heat capacity (
ceq_rest, Equation (6), page 7) of objects inside the store (food products, equipment elements, manipulator and feeders).
Approximately 80% of the total capacity of the goods in the store consisted of liquids, which, in the model, was assumed as 40% of the store volume with the heat capacity of water. In the model, the equivalent power (Peq_eq, Equation (6), page 7) generated by the elements inside the store (including power supplies, electronic systems and electric drives) was also determined.
The total electrical energy consumption of the store HAC unit for the measurement day amounted to 22.49 kWh, while for the model, it amounted to 23.10 kWh. The absolute error was equal to 0.61 kWh.
The third test that was carried out over the course of the research concerned/aimed at investigation on the model for a set of external temperature values obtained from the meteorological station [
36] for the Krakow region, southern Poland. The data set covered a measurement period of 72 h, and the sampling period was 60 min. The outside temperature varied from 12 to 19 °C.
In
Figure 7, time histories of the external (
TOUT Real) and internal (
TIN Real) temperatures registered by the store measurement system, external temperature
TOUT WS registered by the meteorological station, and the temperature inside the store
TIN MOD WS estimated from the model are presented.
The last step of model validation consisted of a comparison of the measured and estimated energy consumption of the HAC unit. A comparison of time histories of the HAC unit’s electrical power, measured (
PAC_Real) and estimated for the assumed model (
PAC Mod K), is presented in
Figure 8.
The values of total energy consumption of the store (
EAC Real) and the total energy consumption (
EAC Mod R and
EAC Mod WB) calculated by the HAC model for two sets of external temperatures (
TOUT Real,
TOUT WS) are listed in
Table 1. Obtained results made it possible to calculate the absolute error (
AE) and the percentage absolute error (
APE). For the energy consumption
EAC Mod R, these errors amounted to 1.64 kWh (2.4%), while for
EACMod WB, it amounted to 0.49 kWh (0.7%).
The proposed model was also tested in the conditions of switching on the HAC unit heating mode. In order to verify the model correctness, the measurements gathered during 3 consecutive days (from 10 January to 12 January 2021) were used. The recorded external (outdoor) temperature varied in the range from 12 °C to −8 °C, and the average value was approximately 1 °C.
Analogously to the previously considered case, the errors for the HAC unit energy consumption calculated from the model did not exceed 2% for the external temperature measured by the store measurement system nor 1% for the external temperature recorded by the meteorological stations.
3.2. Model of Feeding System
In the first developed model of energy consumption by the store’s feeding system, on the basis of carried-out experiments, a linear energy consumption coefficient per product was determined.
where:
Over the course of the carried-out test, the first model was denoted by Mod 0. The value of the energy consumption coefficient per product was determined on the basis of measurements of the number of products (
Figure 9) and power (
Figure 10) carried out by the store’s monitoring system in a 24 h period. In total, 151 products were collected.
The results of the conducted experiment are presented in
Figure 11. The maximum value of energy consumption per product (
Figure 11, marked as Max.) amounted to 202.91 W/product, while the minimum value was 54.31 W/product. The obtained average value of energy consumption per product (
Figure 11, Average) was equal to 147.47 W/product, and in further research, it will be used as the energy consumption coefficient per product ((9),
fpp).
For the obtained energy consumption coefficient per product, a preliminary test was performed, which made it possible to determine the error between the measured value of power
PF Real (
Figure 10) and the power value
PF MOD 0 (
Figure 12) estimated from the model.
In the case of three products being purchased, the maximum error regarding power estimation was 131.1 (W). For the considered experiment, which was carried out for a period of 24 h, the total energy consumption of the feeding system amounted to 0.372 kWh, while the value calculated by the model was 0.371 kWh, and thus, the absolute percentage error (APE) was equal to 0.3%.
The second proposed model of energy consumption of the store’s feeding system was described by Equation (10) and denoted by Mod 1. Based on the orientation of the global coordinate system (
Figure 3), the following notation was assumed (10):
where
nx,
ny,
nz—number of drives in individual axes;
Px,
Py,
Pz—the nominal power of the drives; Δ
tx, Δ
ty, Δ
tz—drives activation time;
nop—number of products.
The average activation times of the drives in individual axes (Δtx, Δty and Δtz) were determined experimentally. For this purpose the following drive controllers were used: EL7211-9014 for AM8112 drives and AX8206 for AM8023 drives manufactured by Beckhoff.
These systems allow for the recording of the position, speed and acceleration, taking into account mechanical transmission ratios, with a maximum sampling frequency of 16 kHz, and the measuring of the motor current with a maximum sampling frequency of 32 kHz. The carried-out experiments aimed at determining the average activation (operation) time of drives in the consecutive axes, which is necessary to download the consecutive products form the store. The obtained experimental results are presented in the third column of
Table 2. The scheme of the considered model implemented in the MATLAB/Simulink R2016a environment is presented in
Figure 13.
The difference between the measured power and power calculated from the model is presented in
Figure 14. The maximal value of error, for three products downloaded, amounted to 136.5 W.
In the experiment carried out for a period of 24 h, the total energy consumption for the store’s feeding system amounted to 0.372 kWh, while the result obtained from the model amounted to 0.367 kWh. Therefore, the error of total energy consumption calculated for the store’s feeding system was equal to 5 W (1.4%).
In the next test, the values of energy consumption per product were calculated for a time period of 72 h. The distribution of the number of products downloaded from the store is shown in
Figure 15. The total number of products was 426.
In
Figure 15, the results of the power measurements carried out for the considered system,
PReal, power computed by the model,
PF Mod 0, described by Equation (9) and power calculated by the model
PF Mod 1 are presented and are described by Equation (10).
The difference between the power measured from the real object and power calculated by the second version model is presented in
Figure 16. The maximal error value for a time period of 72 h was equal to 279.7 W for the first version of the model (Mod 0) and 283.3 W for the second version of the model (
Figure 16).
For the purposes of comparing the results obtained from the models, the following error measures values were used: mean absolute error—
MAE (11), mean absolute percentage error—
MAPE (12), and mean squared error—
RMSE—and the variability coefficient—
CV (14).
where
ymod—sample calculated by the model;
yreal—measured sample;
n—total number of samples;
—average value of all measured samples.
In order to calculate the factors affecting data, the real object power measurement data
PAC Real was used. The
PAC Real, which constitutes a set of samples (
yreal), was compared to both the power data sets
PF Mod 0 obtained for the first model and
PF Mod 1 calculated by the second model. The samples of the data sets
PF Mod 0 and
PF Mod 1 were denoted by
ymod. The total number of samples
n corresponded to the number of products taken from the store. The obtained results are presented in
Table 3 and
Table 4.
The total energy consumption for the real object and models is presented in
Table 5. The error of total energy consumption for the first version model (Mod 0) was equal to 0.54 kWh (0.9%), and for the second version model (Mod 1), it was 1.31 kWh (2%).
In the subsequent tests, it was decided that the second version of the model, Mod 1, was to be used. The justification for the choice is presented in
Section 5.
5. Discussion
The first tests that were carried out over the course of the research presented in this paper allowed for the identification of the parameters of the thermal model of the automated store, including the equipment, assortment as well as the HAC unit. Both external temperature measurement data at the ZINU Shop location (recorded by the store’s measurement system) and generalized data for the given region obtained from a meteorological station were used as input data for the model.
The error of the HAC unit’s electrical energy consumption estimated from the data recorded over a three-day period with respect to the results obtained from the developed model (calculated on the basis of the input data in the form of the external temperature of the store) amounted to approximately 2.4%. However, for the outdoor temperature data obtained from a publicly accessible meteorological station, the error amounted to approximately 0.72%. The difference in the values of the above errors results from the different external-temperature time histories used as the model input.
One of the project assumptions concerned calculating the energy consumption of electrical drives per product. Two linear models were proposed: Equation (9),
PF MOD 0; and Equation (10),
PF MOD 1. The first one was based on the power factor per product determined experimentally. The second model took into account the actuation times of the drives, their number and rated power (10). The results of preliminary tests for three measurement days allowed for a comparison of the models. The values of calculated errors and error measures, presented in
Table 3 and
Table 4, showed similar values of deviation in power consumption per product. Taking into account that the motor parameters in the model make it possible to easily introduce the changes without the need to perform further experiments, it is possible to replace the drive in the feeding system with a drive of a different nominal power value since it does not affect the activation time. Therefore, for the further tests, the second model described by formula (10) was selected.
The second stage of the research concerned the tests carried out over the period four months. The carried-out investigation aimed at establishing the influence of the usage of the average daily outdoor temperatures obtained from a meteorological station on the accuracy of the HAC unit model’s calculations. The obtained values of energy consumption (
Figure 18) and the calculated absolute error (
Figure 19) showed that for the months from January to March, the error deviation values for two considered input datasets were comparable (
Table 6 and
Table 7). In April, this difference increased significantly. The largest spread of the absolute error (
Figure 19) occurred for the 97th measurement day (7 April), in which the average temperature measured was 3.4 °C (max.: 7.6 °C; min.: 0.9 °C), and for the meteorological station, it was 1.9 °C (max.: 3.1 °C; min.: 0.5 °C) (
Figure 17). The real temperature during the day remained within the temperature range set by the HAC unit thermostat: 4–7 °C (
Figure 17,
Th,
Tc parameters). This resulted in a relatively low value for the real energy consumption per day of 7.67 kWh (
Figure 18). Therefore, the temperature absolute error for that day for the calculations carried out based on the real temperature amounted to 1253 (16.32%), while for the data from the meteorological station, it amounted to 11,655 (151%).
Based on the results of the HAC unit’s total energy consumption (
Table 8), the largest error occurred for the month of January for the model based on the measured temperature (77.74 kWh, 15%,
Table 9), while for the model based on the temperature from the meteorological station, the error amounted to 68.81 kWh (11.5%,
Table 9). The error value in this month was influenced by both the average temperature, which oscillated between 4 and 7 °C (
Figure 17), and much greater daily temperature fluctuations, e.g., max. 18.37 °C and min. 1 °C on 1 July. The inaccuracies of the model based on the heat capacity for such a case resulted in a larger error than in case of days with lower dynamics of temperature changes.
The maximum absolute error of the total energy consumption of the automated store feeding system (
Table 10) was observed in April (0.34 kWh, 4.7%,
Table 11), in which a total of 3146 products were purchased. For the entire measurement period, the absolute error was 0.09 kWh (0.1%), which is fully acceptable.
In the last stage of work on the model, the energy consumption values (
Figure 4) of the electrical elements marked as “others” were taken into account. The values of total energy consumption in the individual months and for the entire measurement period (
Table 12) were calculated on the basis of the average energy consumption per day, which was equal to 13.5 kWh.
Table 13 shows the total energy consumption of the automated store over a period of 4 months. The maximum error for the model was equal to 59.15 kWh (5.38%)—based on the outdoor temperature measured in April (
Table 14)—and 65.78 kWh (5.69%)—based on the temperature from the meteorological station in March (
Table 14). The absolute error for the entire period was equal to 68.01 kWh (1.62%) and 90.71 kWh (2.15%), respectively.
The obtained results of the electrical energy consumption evaluated for the HAC unit (
Table 8), the store feeding system (
Table 10) and other elements (
Table 12) allowed for the analysis of the percentage contribution of the model elements in the overall energy consumption over a period of four months (
Figure 20). In January and February, the HAC unit consumed about 58% of the total energy, while in March and April, it was about 62%. The percentage contribution of the feeding system to the energy consumption in the entire period was at the level of approximately 1%.
On the basis of the average CO
2 emission for Poland [
37], 742 g/kWh, the values of mass generated for the production of electrical energy were calculated for the individual months of the measurement period (
Table 15). The total weight for four months was equal to 3057 kg, while the total cost of generating the energy used by the automated store was EUR 603.
6. Conclusions
The results presented in this paper should be treated as the first to be obtained on the basis of an energy consumption model developed for an automated store with a container structure, taking into account information on the external temperature and the number of products sold, for a single store location. According to the authors’ knowledge, no similar results have been published so far.
The hermetic container structure of the automated store together with the efficient HAC unit with the “defrost” function allowed for the application of the model based on the thermal capacity, on the basis of which electrical energy consumption was calculated.
In case of the considered automated store, contrary to the models of buildings like, e.g., schools, office buildings or houses, it is unnecessary to take into consideration additional factors such as humidity, wind, number of windows, frequency of door opening, etc.
The modelling of the store’s feeding system by means of a simple linear model proved to be sufficient to calculate the value of energy consumption with an accuracy of over 95% (
Table 12).
Data from the meteorological station (the store external temperature) allows for the estimation of monthly electrical energy consumption for the HAC unit with an accuracy of 89.5%.The accuracy obtained is sufficient to use the model to forecast energy consumption, CO2 emission and costs for a potential store locations.
The HAC unit (
Figure 20) has the largest contribution to the electrical energy consumption in the analysed time period of four months—approximately 60%. On the other hand, the feeding system, which was switched on when the products are taken, has the smallest contribution of approximately 1%.
To sum up, the main advantages of the proposed solution are a small number of the required model input data, which allows for the quick estimation of energy consumption for any location, easy access to the measurement data, and model implementation that does not require high computational power.
The main disadvantage of the developed solution consists of the fact that the model has currently been verified for one location only for the automated store and, therefore, tests for other locations (latitudes) are planned. Currently, the developed model finds applications in forecasting the costs of energy consumption and gas emissions for a precisely defined type of building structure with an efficient HAC unit.
Future work on the model will concern expanding it to include the application of solar panels as an alternative energy source for the automated store.