Wastewater treatment plant (WWTP) is a topic of pivotal importance in providing a sanitary water which is a vital resource for everyone. Protecting water sources and safeguarding the ecosystem is a long-standing challenge for the WWTP industry. WWTP are facing more stringent effluent standards which were formed for a safer ecosystem (Åmand and Carlsson, 2013; Han and Qiao, 2014; Olsson et al., 2014). The WWTP industry must come up with a solution that abide the stringent effluent requirements and is also economical.
Studies have shown that the energy consumption in the biological system such as activated sludge process, biological trickling filters, and membrane bioreactors can be curbed through good control of the aeration system. The issue of energy consumption has been investigated by various researchers and the findings suggest that the aeration section which is needed in the WWTP to detract nitrogen and natural or inorganic carbon in the biological process, contributes to 50–90% of the overall energy requirement of the WWTP (Åmand et al., 2013; Cristea et al., 2011; Ghoneim et al., 2016; Vrečko and Hvala, 2013).
In the last decade, there have been various studies investigating the effectiveness of various controller design utilizing dissolved oxygen (DO) control in lowering the aeration cost. This control configuration is the highlight during that time due to the availability of DO sensor probe that can continuously measure the DO concentration in the tank. The fundamental of using the DO sensor probe is to control the DO supply according to the oxygen demand of the microorganism in the tank. However, this solution has a weakness due to the difficulty in getting the exact value of the actual oxygen demand by the microorganism at a specific time; thus, most of the proposed DO control strategies implemented an elevated DO set point to avoid nitrification failure (Arnell, 2016; Ellis et al., 2009; Medinilla et al., 2020; Uprety et al., 2015). The DO control strategy has been extensively studied and many viable solutions have been developed and proposed, for example model predictive control (MPC) (Cristea et al., 2011; Han et al., 2012; Holenda et al., 2008; Martín et al., 2012; Yang et al., 2013), PID (Du et al., 2018; Husin et al., 2020b; Samsudin et al., 2013), fuzzy and neural network (NN) control (Han et al., 2020).
However, even with the DO control strategy, the aeration cost issues persist as DO control requires aerators and turbines which are operated by electrically powered motors that add extra cost to the system. This calls for a paradigm shift in the choice of methodology to solve the problems of energy consumption and cost of aeration control. This issue was explored in the publication by Åmand and Carlsson (2014) which has suggested that the aeration process can be regulated either using the aeration concentration control or tweaking the DO set point level corresponding to the ammonium (
ABAC is a control strategy which uses
Summary of recent research trend using ABAC control.
Author | Methods | Results | |
---|---|---|---|
PI ABAC | Uprety et al. (2015) | Feedback PID controller for ABAC to adjust DO in all aeration basins and zones | Decrease in supplemental carbon used for denitrification by 53% and overall decrease in energy consumption by 10% |
Várhelyi et al. (2018) | DO cascade, ABAC and combination of ABAC with the control of nitrate and return activated sludge recycles | ABAC combination is the most cost-saving methods (reduction of about 43%) | |
MPC ABAC | Santín et al. (2015a) | Fuzzy control and MPC (Feedforward ABAC) | Total Nitrogen ( |
Santín et al. (2016) | Risk detection of effluent violation using artificial NN, fuzzy controller to improve denitrification/nitrification and MPC to improve DO tracking |
According to the literature, most pilot or real plants use the PI control in their ABAC designs. The PI controllers utilized are configured in a decentralized configuration. This configuration is favorable because there is no need to deal with the coupling problem in a multi-input multi-output (MIMO) system. A PI controller, on the other hand, is known for being susceptible to disruptions and/or changes in the condition of the operation. Although advanced control schemes such as MPC have been shown to yield better results than PI controllers, MPC is known to be computationally complex and difficult to apply in a real plant. All the studies in the literature indicated that the MPC is implemented in simulation platform only, e.g., BSM1 and BSM2. Another observation into the recent research trend is the emphasis on aeration energy cost problem but less toward the pragmatic benefits brought by ABAC control strategy on effluent quality which has not been extensively explored by researchers.
Considering the advantages and disadvantages levied by these publications, an alternative control strategy that is more streamlined with lower complexity is desirable especially if the aim is to apply the controller in the real or pilot plant. Prior work in this area has been limited to the use of traditional control methods such as PI/PID control. Traditional approaches, such as PI/PID, have various drawbacks, including performance loss when applied to a highly nonlinear process (Aguilar-López et al., 2016) and the inability to meet the demand for high performance control as operating conditions vary (Samsudin et al., 2014). In light of these drawbacks, this research proposes a NN ABAC. The suggested NN ABAC’s main advantages are its simplicity and ability to effectively decouple. Interestingly, to the author’s knowledge and based on the literature research conducted, there have been no previous publications in the literature on the deployment of NN ABAC control to date. Unlike previous systems that relied on cascade configuration, this concept combines many cascaded controllers into a single NN ABAC controller. As a result, the objective of this study is to develop a NN ABAC controller which could be used in the BSM1 to improve the effluent quality especially the
The paper is divided into five sections: the first is an introduction, which includes some earlier and related publications. The simulation model is introduced in the second section. In addition to the typical PI ABAC as a comparison, the third section describes the design and implementation of the proposed NN ABAC. The fourth section summarizes the simulation findings and assesses the proposed NN ABAC’s performance, while the fifth section wraps up the work with some recommendations.
The proposed NN ABAC control is implemented using the MATLAB™/Simulink™ software, and the assessment for this study are established on BSM1. BSM1 is a benchmark simulation plant that consists of a model, a control system, a benchmarking procedure, and evaluation processes. Figure 1 illustrates the BSM1, which is based on the activated sludge model no. 1 (ASM1) and consists of five reactors and a clarifier. Many studies in this field have used this plant as a standard benchmark of comparison.
BSM1.
The simulation of BSM1 arrangement begins with initialization using 150 days of stabilization in a closed-loop state. Then, it follows by simulation using the dry weather file, and lastly, it progresses with weather files to be validated. The system turns out to be stable if the steady state is reached.
Various control procedures are measured corresponding to a guideline defines for the plant performance, which entails the Effluent Quality Index (EQI) and the Overall Cost Index (OCI) to evaluate the operating cost. The evaluation also includes the computation of the operating time that the intensity of the effluent is exceeding the limit, as shown in Table 2.
Concentration thresholds of pollutants in the effluent.
Variables | |||||
---|---|---|---|---|---|
Max. values | 18 | 100 | 4 | 30 | 10 |
For reference and comparison purposes, a feedback PI ABAC configuration is developed with the PI
Feedback PI ABAC configuration.
Both PI controllers for
In this study, a two-input single-output (TISO) NN architecture is used. NN is highly valuable and perform various critical tasks like classification (Sarabu and Santra, 2021), prediction (Fazelabdolabadi et al., 2021), clustering, and associations. A strong coupling problem might arise as the
Direct feedback NN ABAC.
The goal of this research is to create an NN ABAC controller that can be used in the BSM1 to improve effluent quality. As a result, the ABAC controller will be designed first, according to the flowchart shown in Figure 4. It is vital to identify the manipulated variables and control variables for the system during the designing phase. After that, the NN ABAC modelling takes place where the regression (R) and mean square error (MSE) are two criteria that are considered for the model. The implementation, which took place in the BSM1 simulation platform, took place only after the all the criteria were reached. This is a critical component because, even if the model is excellent, it does not always result in high plant performance due to overfitting problems.
Flowchart of the proposed NN ABAC methodology.
The NN model for the proposed TISO system is constructed using standard NN modelling procedures. Data pre-processing and data loading are the first two steps, with data pre-processing being crucial for better NN convergence. The model design came next, with the number of hidden neurons considered as a significant factor. Overfitting can occur when there are too many hidden neurons. Several researchers have proposed their methods to fix the number of hidden neurons as listed in Table 3.
Number of neurons suggested by the researcher and the corresponding MSE value
Researcher | Method | Number of hidden neurons | Mean square error (MSE) |
---|---|---|---|
Huang (2003) | 75 | 0.0113080 | |
Jinchuan and Xinzhe (2008) | 28 | 0.0052734 | |
Shibata and Ikeda (2009) | 1 | 0.0089480 |
In a two hidden layers feedforward network, Huang (2003) proposed a new approach for calculating the number of hidden neurons
The second method selected was proposed by Jinchuan and Xinzhe (2008) who stated that the optimum number of hidden neurons is determined by the network architecture’s complexity. As a result, they recommended that the number of hidden neurons be determined using three key components of network architecture:
The lowest MSE was found using the stated equation presented by Jinchuan and Xinzhe (2008) which calculated the number of hidden neurons based on three main components of network architecture, which are
Determining the training method is another crucial phase in the NN modelling process. After modelling, the NN model must be trained. The Backpropagation (BP) algorithm is used to train the model in this study. In NN training, however, a few BP algorithms are available. An examination of five different BP algorithms is performed before determining which method is best for the study. The following settings are used in this evaluation: the hidden layer has five hidden neurons, the tangent sigmoid transfer function (tansig) is used as the hidden layer’s activation function, and the linear transfer function (purelin) is used as the output layer’s activation function.
The flowchart of the NN training is shown in Figure 5. Assuming that the samples to be trained are
Flowchart of the NN training process.
The five selected BP algorithms are Levenberg–Marquardt (LM), Scaled Conjugate Gradient, Broyden–Fletcher–Goldfarb–Shanno (BFGS) Quasi-Newton, Batch Gradient Descent, and Batch Gradient Descent with Momentum. The LM algorithm combines the steepest gradient descent and the Gauss–Newton algorithms, inheriting the Gauss–Newton methodology’s speed and the steepest gradient descent method’s stability (Yu and Wilamowski, 2011). The Scaled Conjugate Gradient algorithm combines the LM algorithm with the Conjugate Gradient technique, which employs a step size scaling mechanism to skip a time-consuming line search learning cycle, making it faster than the BFGS algorithm (Møller, 1993).
Table 4 shows the comparison findings of the five BP algorithms. LM had the lowest mean square error (MSE) of 0.0057795 and the highest R of 0.99019 of all the BP algorithms tested. The Scaled Conjugate Gradient and the BFGS Quasi-Newton yield nearly identical results. The MSE values for Batch Gradient Descent and Batch Gradient Descent with momentum are both poor. As a result, the LM algorithm is the BP training technique of choice for the NN model.
Comparison of five backpropagation algorithms.
BP algorithm | Function | MSE | Epoch | R |
---|---|---|---|---|
Levenberg–Marquardt | trainlm | 0.0057795 | 23 | 0.99019 |
Scaled conjugate gradient | trainscg | 0.0073901 | 27 | 0.98264 |
BFGS quasi-Newton | trainbfg | 0.0074205 | 58 | 0.98849 |
Batch gradient descent | traingd | 0.0543580 | 1000 | 0.92262 |
Batch gradient descent with momentum | traingdm | 0.1869000 | 8 | 0.71436 |
Because the LM method was chosen to train the network in this study, a detailed explanation of the algorithm is required to fully comprehend it. LM algorithm is a combination of the steepest descent algorithm and the Gauss–Newton algorithm. The Gauss–Newton algorithm faces convergence problems like the Newton algorithm for complex error space optimization (Yu and Wilamowski, 2011). The problem can be interpreted as the matrix
The LM algorithm employs the approximation of the Hessian matrix as in Equation (6). During the training process, the LM algorithm shifts between the two approaches. When the combination coefficient is very tiny or almost zero, it employs the Gauss–Newton approach, which employs the approximate Hessian matrix as indicated in Equation (7). When the value of the combination coefficient is quite large, the steepest descent method as in Equation (6) is utilized. When reaching an error minimum, it is vital to switch to the Gauss–Newton strategy as soon as possible because it is substantially faster and more exact. The combination coefficient is reduced after each successful step and is only raised if a tentative step improves the performance function. As a result, the performance function decreases with each iteration of the algorithm:
According to the rule of the LM algorithm in Equation (6), if the error goes down, meaning it is smaller than the previous error, the quadratic approximation on the total error function is working, and the combination coefficient could be changed to a smaller value to reduce the influence of the gradient descent part. If, on the other hand, the error increases, indicating that it is higher than the previous error, it implies that it is necessary to follow the gradient more to find an appropriate curvature for quadratic approximation, and the value of the combination coefficient is increased. The flowchart of NN training utilizing the LM algorithm is shown in Figure 6 (Yu and Wilamowski, 2011).
Flowchart of the LM algorithm.
The training process using LM algorithm starts can be explained as follows: The total error is calculated using initial weights that are created at random. The computation of the Jacobian matrix is carried out. The weight is then updated using Equation (6). The total error is calculated using the new weight. If the current total error increases as a result of the update, the step is retracted, with the weight being reset to the prior value and the combination coefficient being increased by a factor of 10 or some other factor. The step then repeats from Step (1). If the update reduces the current total error, the step is approved, and the combination coefficient is reduced by a factor of 10 or some other factor. Step (2) is repeated with the new weight until the current total error is less than the required value.
The training parameters are given in Table 5. The NN training comes to an end when either of these conditions is met: the maximum number of epochs is reached, or the maximum amount of time is exceeded, or the performance is minimized to the goal, or the performance gradient falls below the minimum performance gradient, or the µu exceeds maximum µu, or validation performance has increased more than maximum validation failures times since the last time it decreased when using validation was used.
Parameter used for LM training algorithm.
Maximum number of Epochs to train | 1,000 |
Performance goal | 0 |
Maximum validation failures | 6 |
Minimum performance gradient | 1e–7 |
Initial µu | 0.001 |
µu decrease factor | 0.1 |
µu increase factor | 10 |
Maximum µu | 1e10 |
To summarize, 28 sigmoid hidden neurons consist of a two-layer NN model, and a linear output neuron is developed. The graphic representation of the NN model is shown in Figure 7, with DO and
The topological formation of the NN.
The proposed NN ABAC model is used as the controller to control the
The implementation of NN ABAC control architecture in BSM1.
The results obtained from the proposed NN ABAC control are compared with the results from PI ABAC, and benchmark PI to ensure its validity. The study focuses on improving effluent quality while maintaining and/or reducing the aeration costs. There are several methods to analyze the data outline in the BSM1. For this study, the NN ABAC control configuration is evaluated using the second level of performance assessment. In this level, the performances are evaluated using three separate categories.
The first category is the average effluent concentration obtained after a 7-day evaluation utilizing dry/rain/storm weather. The result obtained is compared to benchmark PI and PI ABAC which highlighting five key process variables (
The effluent quality limit.
Effluent average | |||||
---|---|---|---|---|---|
Dry | |||||
PI | 2.4783 | 13.0248 | 16.8908 | 48.2470 | 2.7587 |
PI ABAC | 2.5481 | 13.0244 | 15.8626 | 48.2736 | 2.7654 |
NN ABAC | 2.9118 | 13.0233 | 15.3519 | 48.2888 | 2.7689 |
Rain | |||||
PI | 3.1575 | 16.1970 | 14.7159 | 45.4587 | 3.4569 |
PI ABAC | 3.1299 | 16.197 | 14.1804 | 45.4702 | 3.459 |
NN ABAC | 3.2918 | 16.1958 | 13.9606 | 45.47 | 3.4581 |
Storm | |||||
PI | 2.9953 | 15.2935 | 15.8340 | 47.6875 | 3.2065 |
PI ABAC | 2.9965 | 15.2935 | 15.1311 | 47.7043 | 3.2103 |
NN ABAC | 3.2386 | 15.2923 | 14.8198 | 47.7119 | 3.2115 |
The number of times the effluent criteria were not exceeded over the latest 7 days simulation is the second category in the performance evaluation provided by the BSM1. The observation is carried out for
The effluent violations under dry, rain, and storm influent.
% of reduction | |||||
---|---|---|---|---|---|
PI | PI-ABAC | NN-ABAC | vs. PI | vs. PI-ABAC | |
Dry | |||||
17.86 | 11.90 | 11.61 | |||
7 | 5 | 5 | 0.00% | ||
16.82 | 16.52 | 16.67 | +0.91% | ||
5 | 5 | 5 | 0.00% | 0.00% | |
Rain | |||||
11.01 | 6.10 | 5.65 | |||
5 | 3 | 3 | 0.00% | ||
25.60 | 22.92 | 21.58 | |||
8 | 8 | 8 | 0.00% | 0.00% | |
Storm | |||||
15.48 | 10.86 | 10.71 | |||
7 | 5 | 5 | 0.00% | ||
26.34 | 25.15 | 25.15 | 0.00% | ||
7 | 7 | 7 | 0.00% | 0.00% | |
0.30 | 0.30 | 0.30 | 0.00% | 0.00% | |
2 | 2 | 2 | 0.00% | 0.00% |
The NN ABAC control configuration has demonstrated a marked improvement in the number of effluent violations for
To fully understand the results, the performances of
Performances of the last 7 days of simulation using dry weather with the PI (black line), PI-ABAC (red dotted line), and NN-ABAC (blue line).
Performances of the last 7 days of simulation using rain weather with the PI (black line), PI-ABAC (red dotted line), and NN-ABAC (blue line).
Performances of the last 7 days of simulation using storm weather with the PI (black line), PI-ABAC (red dotted line), and NN-ABAC (blue line).
It was found that only a small decrease of
It is well known that the contaminants
The last category in the performance assessment of the BSM1 is the evaluation of the effluent, aeration, and cost. The results are divided into average EQI, aeration energy cost index (AECI), and the total OCI as per Table 8. As can be seen from the results, the PI ABAC strategy produces the best EQI, while the proposed NN ABAC approach performs just slightly below it, with a deficit of 0.29–0.67%. However, the PI ABAC excess comes at a significant price, as evidenced by the AECI value, which is the highest among the other controllers. The proposed controller, on the other hand, has much lower AECI values, resulting in a reduction of 26.04%. Furthermore, the suggested NN ABAC outperforms the PI ABAC in terms of total OCI reduction by up to 5.15%. In comparison to the PI controller, the NN ABAC consistently outperforms the PI controller across all assessments.
The comparison of EQ, AECI, and Total OCI in dry/rain/storm weather.
% of reduction | |||||
---|---|---|---|---|---|
PI | PI ABAC | NN ABAC | vs. PI | vs. PI ABAC | |
Dry | |||||
EQI (kg poll.unit s/d) | 6,096.71 | 5,938.3021 | 5,978.3177 | +0.67% | |
AECI (kWh/day) | 3,697.57 | 3,769.517 | 2,835.2703 | ||
Total OCI | 16,366.30 | 16,500.995 | 15,689.4197 | ||
Rain | |||||
EQI (kg poll.unit s/d) | 8,146.75 | 8,005.5647 | 8,029.1791 | +0.29% | |
AECI (kWh/day) | 3,671.70 | 3,786.5543 | 2,832.47 | ||
Total OCI | 15,969.35 | 16,133.8675 | 15,302.504 | ||
Storm | |||||
EQI (kg poll.unit s/d) | 7,187.89 | 7,044.115 | 7,079.7043 | +0.51% | |
AECI (kWh/day) | 3,720.76 | 3,830.8403 | 2,833.1054 | ||
Total OCI | 17,328.67 | 17,403.9539 | 16,530.1204 |
Table 9 shows the comparison results in terms of AECI, EQI, OCI and percentage of time over limits of two highlighted effluents,
The comparison of AECI, EQI, OCI, and
Similar studies | Proposed NN ABAC | Husin et al. (2020b) | Husin et al. (2021b) |
---|---|---|---|
AECI (kWh/day) | 2,835.2703 | 3,641.69 | 3,749.24 |
EQI (kg poll.unit s/d) | 5,978.3177 | 6,081.46 | 5,975.75 |
Total OCI | 15,689.4197 | 16,366.30 | 16,435.9 |
11.61 | 15.77 | 13.8 | |
16.67 | 16.82 | 16.07 |
In summary, the study’s goal was accomplished satisfactorily. The NN ABAC controller was designed and applied effectively in the BSM1, yielding considerable results, particularly in terms of effluent violations. In dry and storm weather, the number of
The NN ABAC control configuration could be enhanced more in the future by taking into account the design constraints. Aside from that, injecting qEC in the first tank might be deemed to entirely suppress Ntot breaches below the authorized limit. However, the results produced using this method cannot be compared to those obtained without augmenting the qEC.