1. Introduction
The cement industry is a significant contributor to global energy consumption and CO
2 emissions [
1], responsible for an estimated 7–8% of worldwide CO
2 emissions [
2]. With a current global production of 4.4 billion tons [
3], both production and demand are expected to grow substantially by 2030 [
1,
3]. The immense energy demand of cement production, typically consuming 3.4 GJ per ton [
3], highlights the urgency for innovative and sustainable alternatives to traditional cement production methods. Essentially, the conventional production of cement begins with mining, grinding raw materials into fine powders and heating the fine mixtures to form clinker [
3,
4]. The clinker formation process, which is the main stage in the production of cement, consumes over 70% of the total energy required [
2,
3,
4]. Reducing energy consumption and carbon emissions during cement production is critical as it affects both the cost of production and the impact on greenhouse gas emissions. A large portion of CO
2 emissions in the cement industry comes from the clinker production process, which involves heating limestone raw materials and burning fossil fuels. Addressing both energy consumption and carbon emissions in cement production is essential for safeguarding the environment, ensuring economic sustainability, and improving public health [
5].
Clinker formation involves heating to calcinate a raw mix of calcium carbonate (CaCO
3) and alumina–silicate materials such as clay in a rotating kiln to approximately 1450 °C [
6,
7,
8,
9,
10]. This process causes the thermal breakdown of limestone, releasing carbon dioxide (CO
2) and producing calcium oxide (CaO) [
6]. Concurrently, the remaining oxides in the mixture undergo solid-state sintering reactions [
7]. This high-temperature fusion, which inhibits full melting, promotes the formation of the main phases of tricalcium silicate (Ca
3SiO
5) and minor phases of dicalcium silicate (2CaO.SiO
2) alongside tricalcium aluminate (3CaO.Al
2O
3) [
8]. Essentially, the energy used to heat the raw mix to start chemical reactions usually radiates from the kilns through the kiln shell [
9], amounting to 15–16% of the total energy input [
10]. The hot gases also leave the kiln and carry heat energy in the form of CO
2 and other combustion products. Unfortunately, around 35–39% of the total energy input is wasted during the heating process [
11].
Energy consumption during clinker formation is highly dependent on the combustion system of the kiln [
12], and optimising the performance of the kiln entails reducing energy losses during clinker formation [
11,
13]. The applicability and effectiveness of optimising the performance of the kiln vary between different kilns depending on several factors, such as kiln design and type, technology and equipment, operational practices, availability of alternative fuels, and the cost of modifications relative to the potential savings from reduced energy consumption [
13]. One of the effective strategies for reducing energy losses for kilns is predictive maintenance, essentially done by adding extra sealing and insulation to the kilns and using superior refractory materials in the kiln lining [
14]. Equally, other strategies for reducing energy losses for real-time kiln operations involve integrating complex control systems, such as implementing waste heat recovery systems to capture and reuse heat from the kiln’s exhaust or using advanced process control (APC) systems, which predict the performance of the kiln and make adjustments in real time to optimise fuel consumption and clinker quality [
15,
16]. Furthermore, energy improvement strategies for the kilns also involve employing model-based algorithms that can learn and predict changes to processes and materials to enhance energy optimisation [
17,
18].
Optimisation of energy consumption for cement production, using a model-based approach, can improve kiln performance and reduce energy losses [
19]. The data-based approach usually involves the use of data analytics, machine learning, and other computational methods that establish process predictions and improvement scenarios for various processes involved in clinker production [
19,
20,
21]. Some studies [
20,
21,
22,
23,
24] have also suggested that model-based optimisation techniques offer a holistic and data-driven methodology for analysing energy consumption and improving energy efficiency in cement production [
20]. Further, integrating process modelling with advanced controls, and predictive algorithms is one of the most feasible techniques to optimise energy-intensive operations, such as kiln heating, grinding, and clinker formation [
21,
22]. Additionally, leveraging real-time data and sensor technologies enables the development of dynamic models that adapt to evolving production conditions and variables [
23]. These models can provide insights into the energy performance of the system, identify areas of inefficiency, and suggest optimal operating conditions [
24].
Recent advances in research to improve energy efficiency have focused on integrating model-based approaches with sustainable methods to strike a balance between improving plant performance regarding energy consumption and lessening environmental effects [
25]. These sustainable approaches include the use of alternative fuels such as biomass to potentially lower carbon emissions [
26], together with advanced technologies that involve optimising equipment performance [
26,
27]. Other sustainable methods that have been studied include the use of cement blends in raw materials and feeds to reduce the carbon footprint of cement manufacturing [
28]. Furthermore, the utilisation of waste materials as fuel sources, such as municipal solid waste and tire-derived gasoline, is being studied for its potential to reduce environmental impact [
29].
Mathematical modelling has also played an essential role in optimising production and energy consumption for the cement industry for several decades, with regression analysis mostly used in model formulation due to its capability to identify correlations between multiple variables [
30]. A fundamental regression equation usually connects a dependent variable, such as energy consumption, to one or more independent variables relevant to production [
31,
32]. This type of analysis has proven beneficial in predictive analytics and process optimisation, enabling producers to turn large amounts of operational data into usable information [
33]. Regression models also help in understanding how numerous parameters, such as kiln temperature and raw material properties, influence energy consumption, and this allows for focused interventions to improve efficiency. The ability to predict outcomes based on current and historical data using regression modelling is one of its unique advantages that promotes successful planning and decision-making, which is imperative for improving energy efficiency [
34].
The use of other modern control techniques, such as PID controllers, in cement production helps in improving energy efficiency [
35]. PID controllers use feedback techniques to reduce the difference between desired and actual output values. Control signals are generated efficiently by taking into account past, present, and anticipated errors [
36]. Modified PID controllers have shown great potential for optimising cement plant processes, specifically improving energy efficiency [
37]. The advantages of using advanced approaches with optimised and modified PID controllers include managing the nonlinear, multi-variable dynamics and process restrictions in cement plants [
37,
38]. Model simulation of dynamic systems is also one of the advanced approaches that have been used in energy control systems to minimise energy use, boost efficiency, and reduce environmental pollutants while maintaining production standards [
36,
39]. Implementing an integrated systems approach to energy optimisation in cement plants can result in significant cost reductions [
39,
40].
Despite extensive research that has been conducted [
12,
13,
14,
15,
16,
17,
18,
19,
20,
21,
22,
23,
24], on the use of model-based approaches to optimise energy consumption in cement production, there is still a considerable gap in knowledge on the systematic establishment of an energy consumption model in the context of clinker production, which explicitly links the optimisation of energy consumption to specific improvement strategies. Furthermore, while some studies [
25,
26,
27,
28,
29,
30,
31,
32,
33,
34,
35,
36,
37,
38,
39,
40] have also explained the potential of leveraging advanced analytics and machine learning for process optimisation, there remains a gap in developing models that continuously learn and adapt to new data, providing real-time, actionable insights and improvement strategies. The importance and urgency of the study in the current field of cement production is to address the ever-increasing challenges of climate change, the energy crisis, cost of production, and the need to reduce greenhouse gas emissions [
36]. The objective of this study is to develop a model that explains energy efficiency improvements during clinker production by incorporating energy-intensive parameters of the cement kiln system. The model quantifies the heat loss caused by each parameter and its impact on the overall kiln system.
Two semi-automated cement production facilities in Zambia, known as Plant A and Plant B, provided the production and process data needed for a complete analysis under a range of production schedules and operating conditions. The strategy was to first determine the important energy-intensive parameters as a function of energy consumption from two semi-automated cement plants. The relationship between clinker production and energy consumption was determined using the regression model. A PID-controlled kiln system was established using Simulink software, producing a dynamic model that can adjust to production conditions in real time, offer useful insights, and suggest the best operational plans.
3. Results
This section presents the findings, with an emphasis on the energy loss variables and predictive models created for two cement production plants, referred to as Plant A and Plant B. Pareto charts and regression analysis are used to analyse and present the factors that influence energy consumption during clinker production. Furthermore, the relationship between clinker production and energy consumption is investigated, offering insight into potential solutions for increasing energy efficiency in these plants. The findings presented in this part provide useful insights into the primary causes of energy inefficiency in cement production, as well as actionable recommendations for improving overall energy performance.
3.1. Selection of Energy Intensive Parameters
Figure 2 shows a Pareto chart analysing the energy loss details in Plant A. The graph identifies and prioritises the factors that contribute to energy loss. The graphic indicates that the total energy loss in Plant A is 1889.07 MJ per year. The graph also shows that the most significant sources of energy loss are the flue gas temperature leaving the kiln (
FgTK), clinker discharge temperature (
CDT), flue gas temperature at the recovery point (
FgTR), and overload flue gas temperature (
OFgT). These four parameters alone account for 96% of total energy losses, as seen in Equation (11).
Figure 2 also highlights essential opportunities where solutions to improving energy efficiency can be most successfully applied by identifying these key factors. For example, high temperatures in the flue gases leaving the kiln indicate large heat losses, meaning that enhancing heat recovery systems or optimising kiln operation might significantly minimise these losses. Similarly, high clinker discharge temperatures suggest inefficient heat transfer processes in the kiln and cooling systems. Enhancing the efficiency of clinker cooling can result in significant energy savings. High temperatures at the flue gas recovery point indicate possible inefficiencies in heat recovery units, implying that improvements could include greater insulation, enhanced heat exchanger design, or sophisticated control systems to maximise heat recovery. The factors should be considered for model design with a direct correlation to overall energy losses for Plant A.
Figure 3 shows a Pareto chart for energy loss details in Plant B. The figure provides a clear visual representation of the numerous aspects contributing to energy losses, allowing for prioritising energy improvement measures.
Figure 3 indicates that the total energy loss for Plant B is 1839.87 MJ annually. The graph shows the most significant sources of energy loss, with the tertiary air duct damper temperature (
TaDT), cooler air exhaust temperature (
CaET), flue gas temperature at the recovery point (
FgTR), and overload flue gas temperature (
OFgT) being the main contributors. These components account for 95% of total energy losses, as seen in Equation (12).
The substantial proportion of energy loss (95%) attributed to these specific factors suggests that these are important factors for designing the model that can optimise energy consumption for Plant B.
3.2. Correlation of Energy Consumption and Factors Affecting Energy Loss
Figure 4 shows the regression plot for clinker production for Plant A, which is based on a regression model that links the amount of energy consumed to clinker production. The graph compares the observed data points to their expected values. The figure shows that the R
2 value of 0.749 indicates a strong linear relationship, implying that clinker production volumes account for 74.9% of the variability in energy consumption.
Equation (13) describes the linear relationship between energy consumption and clinker production. According to this model, assuming no clinker is produced (
= 0 tons), the projected energy consumption per ton is 1850 MJ. This baseline energy consumption covers all of the plant’s fixed energy requirements.
where
is the energy consumed, and
represents clicker production in tons.
The regression model predicts that as clinker production increases, energy consumption per ton will decrease. This is supported by the regression equation’s negative slope coefficient (−0.016). The negative coefficient indicates that with each new ton of clinker produced, the energy consumption per ton drops significantly. This suggests that the plant becomes more energy efficient as production increases. Several factors contribute to the reduction in energy consumption with increasing tons of clinker produced. One major factor is the increased utilisation of the fixed energy requirement for the plant. When production is low, the energy required to keep the plant running (baseline energy) is distributed among fewer units, resulting in increased energy consumption per unit. As production increases, the fixed energy cost can be spread out over a larger number of units, lowering energy consumption per ton.
Figure 5 shows the regression plot for clinker production for Plant B, which is based on a regression model that links the amount of energy consumed to clinker production.
The graph compares the observed data points to their expected values according to a normal distribution. The linear relation is expressed by Equation (14), which describes the relationship between energy consumption and clinker production.
where
is the energy consumed, and
represents clicker production in tons.
The positive slope coefficient indicates that for every one-unit increase in clinker production, the energy consumed increased by 0.0045 units. In this context, if clinker production increases, the model predicts a slight decrease in energy consumption per unit of clinker produced. This suggests a trend towards greater energy efficiency with increasing production. The intercept value (583.7) suggests the baseline energy consumption that occurs regardless of production. As clinker production increases, the model predicts a decrease in energy consumption per unit of clinker produced, as indicated by the negative slope coefficient. This relationship implies some level of energy efficiency improvement or optimisation with increasing production levels. The R2 for model fit for Plant B is 0.83.
Figure 6 shows an energy consumption order plot, which indicates the normalised residuals of energy consumption data vs. the observation order for Plant A. The line plot joins each observation’s standardised residuals, allowing for visual detection of patterns, trends, or deviations. In a well-fitted regression model, residuals should be evenly distributed around the zero line, showing non-systematic tendency. However,
Figure 6 shows large variances, especially in the 4th and 12th observations, where standardised residuals exceed 2. These points could be regarded as outliers, suggesting considerable differences between the model’s predictions and actual values.
Figure 7 shows an energy consumption order plot, which indicates the normalised residuals of energy consumption data vs. the observation order for Plant B.
As shown in
Figure 7, there is a large variance, particularly at the 4th and 10th observations, where standardised residuals exceed 2. These points may be considered outliers, indicating significant discrepancies between the model’s predictions and actual values.
In
Figure 8, the logarithmic energy consumption curve, along with confidence intervals for Plant A, illustrates how energy efficiency improves with higher production volumes. It also highlights how energy consumption varies across different production levels.
According to the graph, clinker production for Plant A decreases the amount of energy consumed, indicating an inverse relationship. Solid lines indicate the predicted energy consumption based on clinker production, while dotted lines represent confidence intervals, which indicate the range of actual energy consumption values that are anticipated to decrease within a specific possibility. The shaded background indicates the area covered by the confidence intervals, illustrating the range of variability around the main trend. At higher volumes, the accuracy of energy consumption projections is enhanced as the confidence intervals diminish, demonstrating improved process stability and control. In particular, the lower end of the energy consumption range at peak output emphasises the importance of high-volume production for improving efficiency.
Figure 9 also shows that energy consumption for Plant B increases as clinker production increases. Initially, the energy consumption is relatively low for lower levels of clinker production and increases as production rises. This suggests that as more clinker is produced, the process becomes more energy-efficient. The confidence intervals indicate that there is variability in energy usage, which could be caused by numerous operational or external factors, such as fluctuations in raw material quality or weather conditions. Despite this variability, the overall trend shows improvement in energy efficiency with increased clinker production.
3.3. Kiln System Simulation
The average temperature of the gas entering the burner was 26.7 °C, the volume was 25 m
3, and the density of air was 1.2 kg/m
3. The system response using the equation can be solved as follows:
The transfer function of the kiln for Plant A, calculated using the actual volume of the kiln and inlet temperature, is as follows:
For Plant B, the transfer function of the kiln is as follows:
The transfer function of the heat sensor is calculated using Equation (8) and indicates the transfer function results as follows:
The transfer function of the PID is calculated using Equation (4):
The transfer function for the heat loss is calculated using Equation (2):
where
is defined as
, and
is the heat loss to the environment.
Therefore, the closed-loop equation for the kiln system model is as follows:
The energy efficiency of the system can be defined as follows:
The improvement is established as shown in Equation (23).
Figure 10 shows a Simulink schematic of a PID-controlled system that maintains the cement kiln temperature at 1450 °C. The system determines the difference between the desired and actual temperatures, which are subsequently controlled by a PID controller with proportional, integral, and derivative components. This PID controller controls the damper to manage the airflow and fuel mixture to the burner, with a transfer function indicating the dynamic response. The thermal dynamics of the kiln are also represented by a transfer function, which specifies how it responds to heat input. In addition, the system accounts for heat loss using another transfer function. A complete heat detector measures both generated and lost heat, providing feedback to maintain accurate temperature control.
Figure 11 shows the dynamic response model for the burner in the cement kiln system, which is modelled using the transfer function in Equation (13) to control inputs over time. This function describes a first-order system, in which the temperature response increases exponentially towards a steady-state value. The response curve indicates that the burner output can be modified to attain the ideal kiln temperature of 1450 °C, which is required for clinker formation. The constant of time τ, obtained from the denominator of the transfer function, is 60 s, showing the time required for the system to achieve approximately 63.2% of its final value. The time constant shows small values, which indicate how fast the system responds.
Figure 11 also indicates that the system stabilises after around 100 s without exceeding the specified temperature, proving that the PID controller is effectively tuned. The absence of overshoot and oscillations in the graph indicates that the control system provides stability and has the ability to efficiently handle process disruptions, emphasising the significance of an optimised burner response for energy efficiency.
Figure 12 presents the response function of a heat sensor in a cement kiln system. The graph indicates that the temperature immediately rises to roughly 1400 °C before stabilising at this level, suggesting that the system quickly achieves and maintains the desired temperature. The system stabilises after 2 s, which indicates that the control mechanism, likely a PID controller, is well tuned to respond swiftly and effectively. The absence of temperature overshoot in the curve suggests that the control system is highly effective in ensuring that the temperature does not exceed the desired level.
Figure 13 shows the temperature response of a kiln system over time, controlled by a PID control system. The figure shows the temperature rises rapidly and exceeds the set point, but then gradually decreases and settles around the ideal temperature of 1400 °C. The system settles in approximately 40 s, indicating how long it takes the PID controller to adapt and maintain the desired temperature without oscillations. This type of response is characteristic of PID-controlled systems, designed to reach the setpoint efficiently while reducing overshoot and ensuring rapid stability. The graph shows how the control system manages kiln temperature, responds fast to changes, and achieves steady-state conditions. This results in high efficiency and consistent product quality, demonstrating a powerful control system that meets contemporary cement production standards.
The graph in
Figure 14 presents the temperature response for a PID-controlled and non-controlled system in Plant B. The PID-controlled system initially reaches a peak temperature before settling to a stable temperature of approximately 1200 °C. This mirrors the actions of the PID controller, which responds quickly to reduce error and stabilise the system. The system without control has a significant overrun, reaching temperatures above the desired setpoint. It exhibits more noticeable oscillations, indicating a lack of control and difficulty in stabilising quickly. The graph validates the model’s behaviour by showing that the model for Plant A matches the effective control dynamics provided by the PID system in Plant B, proving its precision and applicability.
4. Discussion
According to the data collected and analysed, the annual energy loss for Plant A accumulates to 1889.07 MJ. The four major factors that contribute to this energy loss are overload flue gas temperature (OFgT), clinker discharge temperature (CDT), recovery point flue gas temperature (FgTR), and flue gas temperature (FgTK). Pareto analysis indicates that these four variables account for 96% of all energy losses. A significant reduction in energy loss could be achieved by implementing improvement strategies that are related to these four variables. The improvement strategies must include insulation and heat recovery. In recent studies [
41,
42,
43] and other studies [
44,
45,
46,
47], it has been shown that optimising heat recovery systems and improving insulation are effective strategies for improving energy efficiency in industrial processes. Therefore, by focusing on these areas, Plant A can significantly reduce its energy consumption and align with best practices in cement manufacturing energy management. Regression analysis for Plant A indicates that there is a significant relationship (R
2 = 0.749) between clinker production and energy consumption. The analysis indicates that energy consumption per ton reduces with increasing production volume. This is because fixed energy costs are distributed across higher production volumes, indicating substantial efficiencies in production.
Plant B, on the other hand, has an annual energy loss of 1856.5 MJ, with a different set of causes. The tertiary air duct damper temperature (TaDT), cooler air exhaust temperature (CaET), flue gas temperature at the recovery point (FgTR), and overload flue gas temperature (OFgT) all contribute significantly to energy losses, accounting for 95% of the total. Plant B has different production and energy consumption trends compared to Plant A. The regression model shows a stronger association (R2 = 0.83) between energy usage and clinker production. This implies that Plant B benefits from a lower drop in energy consumption per unit produced, indicating different operational constraints or efficiency compared to Plant A. The variations in operational efficiency could be attributed to differences in equipment, process technology, or operational methods. Using Pareto analysis, Plants A and B can identify solutions tailored to their specific energy challenges and prioritise them accordingly.
The dynamic simulations of the cement kiln system using Simulink indicate that the use of a PID-controlled system is effective at managing kiln temperature and maintaining process stability. The burner and kiln transfer functions define how well the system responds to input changes, resulting in a quick temperature rise to the desired setpoint with little overshoot. This rapid stabilisation within a short timeframe indicates that the PID controller is precise and effective in minimising fluctuations and maintaining steady-state conditions, as supported by studies [
48,
49,
50]. The PID controller helps optimise energy usage by minimising temperature variations, which has been highlighted in the literature as an important component in improving system sustainability [
50]. The rapid stabilisation is consistent with research findings [
51,
52], which indicate that well-tuned PID systems can greatly enhance energy efficiency, adding to the sustainability goals emphasised in current research. Furthermore, the absence of significant fluctuations and overshoot in the response of the system demonstrates precision in achieving desired operating conditions. This level of control indicates that the model is both energy-efficient and capable of delivering consistent controls required for product quality. These findings are consistent with previous research [
48,
49,
50,
51], which emphasises the importance of rigorous control strategies in cement manufacturing to preserve product quality while optimising energy use.
The energy efficiency equation expressed in Equation (21) provides a metric for determining the energy efficiency of a cement kiln system over a certain time. This equation emphasises the magnitude of the closed-loop transfer function, which describes the dynamic reaction of the system to inputs and disturbances while ensuring stable functioning. By integrating this response over time, the equation gives useful information on the effectiveness of the PID controllers, specifically its capacity to reduce temperature fluctuations and stabilise the system. For Plant A, where considerable energy losses are caused by factors such as high flue gas and clinker discharge temperatures, reaching a higher energy efficiency measure requires effective control over these parameters. This can lead to less energy waste and more operational efficiency. Similarly, in Plant B, where air duct and cooler exhaust temperatures are the key sources of energy loss, an elevated measure of energy efficiency would indicate first-rate management of these inefficiencies. This can ensure the kiln operates under optimal circumstances for clinker formation. By comparing the dynamic behaviour of Plant A and Plant B, the model for Plant A can be validated against the observed performance in Plant B, indicating that PID control effectively captures the required control dynamics. This comparison demonstrates that the model accurately reproduces real-world behaviour, demonstrating its precision and suitability for predicting and optimising kiln operations under identical conditions.
In both plants, optimising the PID control system to improve energy efficiency not only contributes to energy saving but also supports sustainable cement production by lowering the carbon footprint and operational expenses associated with excessive energy consumption. As a result, the established energy efficiency equation incorporates the identified areas for improvement as a function of heat loss and ensures that both plants waste as little energy as possible, in line with industry sustainability and cost-effectiveness standards.
5. Conclusions
Factors affecting energy efficiency were identified from the collected production data from two semi-automated plants. The analysis of Plants A and B revealed that high energy losses at Plant A are related to high flue gas temperatures and clinker discharge temperatures, which accounted for 96% of energy inefficiencies. The energy losses for Plant B were related to issues such as the tertiary air duct damper temperature, cooler air exhaust temperature, and other related factors, which account for 95% of total losses.
The regression analysis for Plant A established a strong correlation between energy consumption and clinker production, emphasising the advantages of higher production efficiency because of lower energy consumption per unit of clinker produced. The regression model for Plant B established an even stronger relationship between energy consumption and production, indicating that further improvements in techniques can greatly increase energy efficiency.
The dynamic behaviour of the PID-controlled system showed a rapid stabilisation of kiln temperature, quickly reaching the desired setpoint with minimal overshoot. The system effectively minimised fluctuations, maintaining steady-state conditions with precise control, contributing to improving energy efficiency. The absence of significant overshoot further confirmed the accuracy of the system in achieving the desired operating conditions. This stable and efficient performance ensured consistent process stability, essential for maintaining product quality in cement manufacturing.
The energy efficiency equation was developed as a mathematical model that incorporates heat loss to the environment as a system disturbance. This model captures the dynamic response of the system by focusing on the closed-loop transfer function, which is critical for maintaining steady operation. This approach simulates real-world energy losses and disturbances and indicates that the established model adapts to the dynamic operational requirements of processes like cement production, where maintaining stability and efficiency despite varying conditions is essential.