Next Article in Journal
Household-Reported Availability of Drinking Water in Africa: A Systematic Review
Next Article in Special Issue
CFD Modeling of a Stirred Anaerobic Digestion Tank for Evaluating Energy Consumption through Mixing
Previous Article in Journal
Optimal Resilience Enhancement of Water Distribution Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Parameter Estimation in Activated Sludge Process Based Wastewater Treatment Practice

1
College of Electrical and Information Engineering, Lanzhou University of Technology, Lanzhou 730050, China
2
School of Engineering, Royal Melbourne Institute of Technology (RMIT University), Melbourne 3000, Australia
3
Key Laboratory of Gansu Advanced Control for Industrial Processes, Lanzhou University of Technology, Lanzhou 730050, China
4
National Demonstration Center for Experimental Electrical and Control Engineering Education, Lanzhou University of Technology, Lanzhou 730050, China
*
Author to whom correspondence should be addressed.
Water 2020, 12(9), 2604; https://doi.org/10.3390/w12092604
Submission received: 3 July 2020 / Revised: 8 September 2020 / Accepted: 12 September 2020 / Published: 17 September 2020

Abstract

:
Activated sludge models (ASMs) are often used in the simulation of the wastewater treatment process to evaluate whether the effluent quality parameters of a wastewater treatment plant meet the standards. The premise of successful simulation is to choose appropriate dynamic parameters for the model. A niche based adaptive invasive weed optimization (NAIWO) algorithm is proposed in this paper to find the appropriate kinetic parameters of activated sludge model 1 (ASM1). The niche idea is used to improve the possibility of convergence to the global optimal solution. In addition, the adaptive mechanism and periodic operator are introduced to improve the convergence speed and accuracy of the algorithm. Finally, NAIWO is used to optimize the parameters of ASM1. Comparison with other intelligent algorithms such as invasive weed optimization (IWO), genetic algorithm (GA), and bat algorithm (BA) showed the higher convergence accuracy and faster convergence speed of NAIWO. The results showed that the ASM1 model results agreed with measured data with smaller errors.

1. Introduction

Activated sludge process is a complex biochemical reaction process, which mainly uses the metabolism of microorganisms present in the activated sludge to remove the organic pollutants that are present in the wastewater. Therefore, it will be affected by external conditions such as the environment and the temperature of the wastewater. Activated sludge process is one of the most commonly used wastewater treatment methods, especially for large-scale urban wastewater treatment plants. Activated sludge model 1 (ASM1) was proposed by the International Water Association (IWA) in 1987. The ASM1 is a complex mathematical model with 65 differential equations and 19 kinetic or stoichiometric parameters [1]. ASM1 is the most commonly used and most researched simulation model in the field of activated sludge process [2,3,4,5]. However, in actual wastewater treatment plant applications, due to changes in operating conditions and differences in the external environments, the dynamics or stoichiometric parameters of ASM1 will have to be determined. If the reference parameters provided by IWA are used, they will cause a serious deviation between the model output and the actual results. Therefore, correcting the parameter values of ASM1 for practical applications is the key to the success of the application of the model [6].
In the early days, the methods of parameter identification of environmental models in environmental engineering mainly included trial and error techniques [7], optimization methods based on gradient descent [8], and statistical methods based on random sampling [9,10,11]. However, these methods are inadequate in improving parameter accuracy or convergence speed and it is difficult to obtain satisfactory results with minimum time and effort. Gradient-based optimization methods rely too much on the selection of initial values and tend to converge to a local optimum. The search mechanism of statistical methods is random. When the number of parameters increases, the number of calculation steps will increase exponentially, becoming inefficient in solving complex problems. With the extensive application of genetic algorithms (GA) in various fields [12], meta-heuristic search algorithms are gradually introduced into parameter estimation of the activated sludge process models. Compared with the above methods, a meta-heuristic search algorithm can converge to the global optimal solution with higher probability at a faster speed [13,14]. Kim et al. [15] used GA to estimate the sensitive parameters in ASM1. The results showed that using the optimized parameters to run ASM1 can reduce the error between the predicted and the actual values, and confirms the applicability of a meta-heuristic optimization algorithm in parameter optimization of a nonlinear system model. Du et al. [16] introduced the adaptive step control algorithm to improve the cuckoo search (CS), targeting the shortcomings of the optimization mechanism of the CS algorithm and the slow convergence speed in the later stages. Outcomes show that the improved CS is more effective in the parameter estimation of ASM1.
Invasive weed optimization (IWO) algorithm, proposed by Mehrabian et al. [17] in 2006, is a new type of meta-heuristic optimization algorithm. IWO mimics the strong colonial dominance power of wild grasses by imitating the growth, reproduction, diffusion, and competition of field weeds. It has been widely used in engineering optimization [18,19], fault diagnosis [20], and hybrid algorithm optimization [21]. In the IWO algorithm, as the number of iterations increases, the spatial distribution of the next generation of seeds will gradually narrow, which can ensure that the algorithm has a strong global search ability in the early stages and a strong local search ability in the later stages. However, it also leads to the lack of local search ability in the early stages and the lack of population diversity in the later stages. In order to get better optimization results, it is necessary to improve and fuse the IWO algorithm. Cuevas et al. [21] proposed a hybrid evolutionary method, which combines the search ability of IWO and probability model of an estimated distribution algorithm, making a hybrid method with higher accuracy, efficiency, and robustness. Liu et al. [19] introduced simplified quadratic approximation into the IWO algorithm and the application in the directional pattern synthesis of array antennas showed the effectiveness of the improved IWO algorithm.
In response to the above problems, the niche idea and adaptive mechanism are introduced into the IWO algorithm in this paper and a niche-based adaptive invasive weed optimization (NAIWO) algorithm is proposed. Niche is used to increase the population diversity of the algorithm. Then, a periodic operator and adaptive algorithm are introduced into the adaptive mechanism, so that the standard deviation of the spatial diffusion of individual weeds not only changes with the number of iterations, but also can dynamically change according to the parameters of the periodic operator and the fitness value of the individual. Finally, verified through parameter optimization of ASM1, the effectiveness of the NAIWO algorithm in global convergence ability and convergence speed is proven.
Figure 1 shows the schematic diagram of the process of ASM1 parameter estimation using the NAIWO algorithm. The input data of the wastewater treatment plant is taken as the input data of ASM1. The predicted data of effluent is obtained from the ASM1 model and the error is calculated by comparing the measured and the predicted data. Then, the error is sent to the NAIWO algorithm for parameter optimization. Parameter optimization is repeated until the error is minimized to a desired level.

2. Materials and Methods

2.1. Invasive Weeds Optimization (IWO)

The IWO algorithm imitates the basic process of weed diffusion, growth, propagation, and competitive survival in the field. “Weeds” represent a feasible solution to the problem and “populations” represent the set of all weeds.
The basic steps of IWO can be expressed as follows:
(1)
Population initialization. P0 of weeds (feasible solutions) were randomly generated in the solution space (D-dimension). Generally, the size of P0 can be adjusted according to the actual situation.
(2)
(Growth and reproduction. After seeds grow and bloom, they produce a new generation of seeds according to their own fitness. The number of seeds produced by the parents is related to the fitness of the parents. The specific relationship is shown in Equation (1):
S N = F N F min F max F min ( S max S min ) + S min ,
where FN and SN represent the fitness value of the Nth parent and the number of seeds it should produce, respectively; Fmax and Fmin represent the maximum and minimum fitness of the parents, respectively. Smax and Smin represent the maximum and minimum number of seeds produced by individual weeds, respectively. The number of seeds generated is a rounded down number of SN.
(3)
Space diffusion. The generated seeds are normally distributed in the D-dimensional space near their parents. The mean value of the normal distribution is 0 and the standard deviation is σiter. The variation of the standard deviation with the number of iterations is shown in Equation (2):
σ i t e r = ( i t e r max i t e r ) n i t e r max n × ( σ i n i t i a l σ f i n a l ) + σ f i n a l ,
where iter and itermax represent the current number of iterations and the maximum number of iterations; σinitial represents the initial standard deviation value; σfinal represents the final standard deviation value; and n represents the nonlinear harmonic coefficient. It is generally ensured that σinitial is greater than σfinal.
(4)
Competition exclusion. After several generations of reproduction, environmental resources will not be able to bear the number of offspring produced. The maximum population size is determined as the preset maximum population number Pmax. When Pmax is reached, firstly reproduce freely according to the previous steps. Then, based on the population upper limit requirement, the parents and children are eliminated together according to the adaptive value.
(5)
Repeat steps 2 to 4 until the maximum number of iterations is reached or the solution satisfying the required conditions is found.

2.2. Niche-Based Adaptive Invasive Weed Optimization (NAIWO)

It can be seen from Equations (1) and (2) that in the IWO algorithm, the individual weeds of each generation determine the number of next generations according to the degree of fitness, so as to ensure that the excellent individual genes are inherited to the next generation. As the iteration goes on, the diffusion range (i.e., the size of the standard deviation σiter) of the next generation seed is gradually reduced, which ensures that the algorithm has a strong global search ability in the early stages and local search ability in the later stages. However, this also leads to insufficient local search capability in the early stages of the algorithm and later searches only near seeds with higher fitness, resulting in a lack of diversity in the later stages of the algorithm. In order to solve these problems, a niche-based adaptive invasive weed optimization (NAIWO) algorithm is proposed by introducing the idea of a niche and adaptive mechanism.

2.2.1. Dynamic Adaptive Mechanism

In step 3 of the IWO algorithm, the distribution of the next-generation seeds generated by each parent follows the same normal distribution and the standard deviation of its search walking is σiter, which decreases with the number of iterations. Although the global search ability in the early stages and the local search ability in the later stages are considered in this distribution method, the population diversity in the later stage of the algorithm is insufficient, which makes the algorithm fall easily into a local optimum. In this paper, the dynamic adaptive mechanism is introduced into the spatial diffusion step of the IWO algorithm to balance the global and local optimization ability of the algorithm.
The spatial diffusion criterion is shown in Equations (3)–(5):
σ i t e r 1 = ( i t e r max i t e r ) n i t e r max n × ( σ i n i t i a l σ f i n a l ) + σ f i n a l
σ i t e r = σ i t e r 1 × ( K × | cos ( π × i t e r T ) | + 1 K )
σ j = σ i t e r × e ( 1 F max F j F max F min ) .
In the dynamic adaptive spatial diffusion mechanism, the spatial distribution σj of the children generated by each parent is shown in Equations (4) and (5). In Equation (4), a cosine periodic function is employed, where T is a periodic parameter and K is a scaling factor. By adjusting the values of K and T, the dynamic diffusion standard deviation σiter can be changed with a period T between [1/K, K]. For the jth parent‘s weed, the distribution of its offspring σj is also related to its fitness in the population, as shown in Equation (5).
In each iteration, Fmax and Fmin are the maximum and minimum fitness values in the parent population; Fj is the fitness value of the jth parent. The distribution of each generation’s population is not only related to the iteration number (iter), but also has a functional relationship with the fitness of the parent in the population. In the iterative process, the higher the adaptability of the weeds, the larger the number of the next generation seeds and the more centralized the distribution of the next generation, making the seeds continuously concentrate to high adaptability. However, for the parents with low adaptability, the number of seeds produced is less in a larger distribution range, so as to improve the possibility of finding the global optimal solution.

2.2.2. Niche Idea

The niche idea is derived from biology and refers to a living environment under a specific condition. In the evolution process of organisms, they generally always live with the same species and jointly reproduce offspring. Each generation is divided into several classes according to the fitness value and each class can represent a niche. The combination of a niche idea and intelligent algorithm shows strong utility [22,23]. In this paper, the niche idea is introduced into the competitive exclusion step of the IWO algorithm, and the characteristics of classification competition of niche is used to increase the diversity of the population and improve the overall optimization ability of the algorithm.
The determination of the radius R of the niche is based on Equation (6):
R = max ( d n i ) × [ k ( i t e r i t e r max ) 1 a + b k ] ,
where d n i represents the Euclidean distance from the ith weed to the most adaptable weed in the iter iteration; a, b, and k are adjustable parameters. By adjusting their values, we can adjust the rate of change and the start and end values of the radius R of the niche accordingly.
The process of dividing niches is described as follows:
(1)
Arrange the weed individuals in the population in descending order according to the degree of fitness. If the population number is greater than the maximum population number Pmax, take the first Pmax weed individuals as the parents for the next generation.
(2)
The center (H1) of the first niche is the position of the weed with the highest adaptability, and R is its radius. If the Euclidean distance (d1,i) of the other weeds in the population from the center of niche H1 is less than R, they would be included into the niche H1. Otherwise, it would be excluded.
(3)
The second niche, H2, is marked by the weed with the highest adaptability among the rest of the individuals not belonging to H1. Repeat the process (2) until all weeds in the population are tagged.
The steps of NAIWO can be described as follows:
(1)
Population and parameter initialization.
(2)
Calculate the fitness of each individual weed and arrange all the individual weeds according to the above-mentioned niche classification method.
(3)
According to Equation (1), the growth and propagation of weeds are carried out to produce seeds.
(4)
According to Equations (3)–(5), the adaptive spatial diffusion based on fitness is carried out.
(5)
Solution is considered optimal when the solution meets the requirements or the maximum number is reached.
(6)
If the current number of individual weeds Piter is greater than the maximum population number Pmax, go to the competition exclusion process in step 7; otherwise, return back to step 2.
(7)
Competition exclusion: Select a certain number of weed individuals from each niche for the next iteration as parents. The number of individuals selected in each niche is related to the ratio of the number of individuals in the niche to the total number of individuals in the population as shown in Equation (7):
N ( i ) = r o u n d ( P max × X ( i ) sum ( X ) ) ,
where Pmax is the maximum population; X(i) represents the number of individuals in the ith niche and sum(X) is the total number of individuals in all niches.
(8)
Repeat steps 2 to 7 until the optimal solution is found or the maximum number of iterations is reached.
Based on the above improvement ideas, a new algorithm (NAIWO) with excellent performance is established. To verify the effectiveness of the NAIWO algorithm, we have done a series of optimization performance tests. The results and more information could be found in Appendix A.

2.3. Application: ASM1 Parameter Optimization for Actual Wastewater Treatment Plants

During the operation of a biological wastewater treatment plant, its internal biochemical reaction mechanisms are extremely complex. Therefore, before the ASM1 was proposed, establishing the process model for biological wastewater treatment plants was always difficult and lacked accuracy; this led to difficulties in controlling the effluent water quality to meet increasingly strict standards using controllers.
The purpose of ASM1 is to describe the reaction mechanisms of the activated sludge process as accurately as possible. After selecting the appropriate reaction parameters for ASM1, the reaction processes of an activated sludge plant can be described more accurately, laying the foundation for a precise control of the effluent water quality parameters.
However, for different wastewater treatment plants located in different places, the environmental conditions and inflow conditions may be quite different. If the recommended parameters given by IWA are continued to be used, the model will lack accuracy when predicting the quality of the treated effluent.

2.3.1. Plants and Data Description

In order to obtain an accurate model of the wastewater treatment plant, the proposed NAIWO algorithm was used to estimate seven important parameters of ASM1, and the process models of two wastewater treatment plants, the Pingliang City Wastewater Treatment Plant (PC-WWTP) in Gansu Province, China and the Wushan County Wastewater Treatment Plant (WC-WWTP) in Tianshui City, Gansu Province, China were established. PC-WWTP is a large full-scale activated sludge process-based wastewater treatment plant with a designed wastewater treatment capacity of 50,000 m3/day and an average daily wastewater inflow of 20,291 m3/day. WC-WWTP is a small biological wastewater treatment plant with a designed wastewater treatment capacity of 8000 m3/day and an average daily wastewater inflow of 5742 m3/day. Figure 2 shows the basic processing of both plants, PC-WWTP and WC-WWTP, based on activated sludge process (ASP).
All the experimental data used in this work comes from the dry weather dynamic inflow and effluent water quality measurement data of these two wastewater treatment plants. The inflow data is collected as an “Influent” input, shown in Figure 2, and the effluent data is collected as an ‘Effluent’ output. The data used in the experiment were measured by the sensors in the two wastewater treatment plants. The inflow wastewater data includes the components of SI, SS, XI, XS, XB,H, XB,A, XP, SO, SNO, SNH, SND, XND, SALK, TSS, and Q0, where TSS represents total amount of solids (mg SS /L) and Q0 stands for influent flow rate (m3/day). The effluent concentration of four representative components, SNH, SNO, SS, and XS, was selected as the standard to judge the accuracy of parameter estimation. A detailed description of these components and reactions is given in Table A5. See more information in Appendix B.
The simulation data were collected in dry weather and the duration was 1 day. The sampling frequency was 15 min, so there were 97 sets of sampling data for computations.

2.3.2. ASM1 Parameter Estimation

The purpose of parameter estimation for ASM1 is to select a set of appropriate parameter values to minimize the errors between the model output and the observed values in the wastewater treatment plant. There are five stoichiometric parameters of biochemical reactions and 14 kinetic parameters involved in ASM1. The correctness of these 19 parameters ensures the accuracy of the ASM1 model in simulating the performance of actual wastewater treatment plants and the specific information of those 19 parameters is shown in Table 1.
Estimation of all the 19 parameters in ASM1 is very complicated. Therefore, seven parameters [24] that have a greater influence on the output results were selected in this section as the estimation objects and the remaining 12 insensitive parameters were used as recommended by the ASM1 model description.
In the process of parameter estimation of ASM1, the sum of squares of relative errors f ( t ) were used as the objective function as shown in Equation (8) in order to minimize the difference between model outputs and observed data:
f ( t ) = i = 1 p j = 1 q ( y i j y i j y i j ) 2 ,
where p is the total number of times the samples were taken and q is the number of effluent quality parameters considered; y i j is the model prediction of the jth effluent quality parameter of the ith sampling time and y i j is the measured value of the jth effluent quality parameter of the ith sampling time. The process of parameter optimization is to find the minimum value of the objective function. Figure 3 shows the flow chart of the NAIWO algorithm used in parameter estimation of the ASM1 model.
By bringing the optimal parameter values estimated by the proposed NAIWO algorithm to run ASM1, this will simulate the effluent quality from an activated sludge plant accurately. In this study, the effluent concentration of ammonia nitrogen (SNH), nitrate nitrogen (SNO), soluble rapidly biodegradable organic matter (SS), and insoluble slowly degradable organic matter (XS) were selected as the objects to verify the validity of the ASM1 parameter estimation.

3. Results and Discussion

In order to prove the advantages of the proposed NAIWO algorithm, seven parameters of ASM1 are estimated by IWO, GA, BA, and NAIWO, respectively. The optimization results are shown in Table 2 (PC-WWTP) and Table 3 (WC-WWTP).
The fitness curves of the four algorithms in the parameter estimation process of ASM1 for the two wastewater treatment plants are shown in Figure 4a,b, respectively. The calculation result of Equation (8) (i.e., the sum of squares of relative errors) is used as the fitness value. Therefore, the fitness value reflects the error between the predicted values by ASM1 and the actual measured values at the wastewater treatment plant. The smaller the error value is, the better the simulation performance.
The shape of fitness curve reflects the change in error during the process of parameter optimization carried out by an algorithm. The faster the fitness curve drops, the faster the convergence speed of the algorithm and the smaller the final fitness value, the higher the convergence accuracy of the algorithm.
It can be seen from the fitness curve shown in Figure 4a (PC-WWTP) that the convergence accuracies of GA and BA are significantly lower than that of IWO and NAIWO. For IWO, its convergence accuracy is only slightly less than NAIWO, but the convergence speed (declining speed of the curve) is significantly lower than NAIWO. In addition, in the process of estimating the parameters of WC-WWTP, the final fitness value of NAIWO is significantly smaller than the other three algorithms and the falling speed of the fitness curve is also significantly higher than that of the other three.
In summary, the NAIWO algorithm can achieve the minimum error with the fastest speed in comparison with GA, BA, and IWO, whether it is applied for large (PC-WWTP) or small (WC-WWTP) wastewater treatment plants.
In order to further verify the effect of parameter estimation, the parameters estimated by the four algorithms were introduced into ASM1 and the error between the predicted data of the ASM1 and the actual effluent data were compared. Ammonia nitrogen (SNH), nitrate nitrogen (SNO), soluble rapidly biodegradable organic matter (SS), and insoluble slowly degradable organic matter (XS) were used as comparative objects. The comparison curves of predicted values and actual values are shown in Figure 5 and Figure 6.
Figure 5a–d and Figure 6a–d are the effluent SNH, SNO, SS, and XS concentration curves of PC-WWTP and WC-WWTP. The red solid line in the figure represents the actual measured water quality parameters of the two wastewater treatment plants. The closer the other curves are to the red solid line, the smaller the error during the simulation of ASM1 and the higher the accuracy of the parameter estimation.
It can be seen that when the recommended parameter values given by IWA were used, the prediction errors of the model were larger. After optimal parameters were obtained by GA, BA, IWO, and NAIWO algorithms, the prediction accuracy was improved to a certain extent.
In Figure 5a, when the concentration of effluent SNH was small, the prediction errors of the four algorithms were small. When the concentration increased rapidly (the 50th–70th sampling points), the four algorithms had obvious differences. Only the parameters estimated by NAIWO could track the change of SNH concentration well. Meanwhile, it can be seen from Figure 5b–d that the NAIWO algorithm had the best performance with the highest tracking accuracy to the measured data of SNO, SS, and XS.
As can be seen from Figure 6, IWO had a similar prediction accuracy to that of NAIWO only when predicting the SS concentration of the effluent. The prediction errors of the NAIWO algorithm for the effluent concentration of SNH, SNO, and XS were much smaller than the GA, BA, and IWO algorithms.
According to the results simulated by ASM1 for the effluent quality of two wastewater treatment plants, it was found that the parameters estimated by GA, BA, and IWO can reduce the prediction error of ASM1 to a certain extent, while the parameters estimated by NAIWO can produce the minimum prediction error. Therefore, NAIWO is more effective in optimizing ASM1 parameters than GA, BA, and IWO.
The results illustrate that for both the Pingliang City Wastewater Treatment Plant (PC-WWTP) and Wushan County Wastewater Treatment Plant (WC-WWTP), the proposed NAIWO algorithm is effective for the optimization of the seven sensitive parameters of ASM1. The ASM1 model with optimized parameters can effectively predict the effluent quality of Pingliang City Wastewater Treatment Plant and Wushan County Wastewater Treatment Plant. These two wastewater treatment plants represent different scales with respect to treatment capacity and different water quality environments. The successful application of NAIWO in ASM1 parameter estimation for these two wastewater treatment plants gives confidence for NAIWO to be applied to other activated sludge process based wastewater treatment plants.

4. Conclusions

In this paper, a niche-based adaptive invasion weed optimization (NAIWO) algorithm was proposed to overcome the shortcomings of invasive weed optimization (IWO), such as insufficient local search ability in the early stages of the iteration and insufficient population diversity in the later stages. On the basis of high convergence accuracy, the NAIWO algorithm achieved the balance of global convergence ability and convergence speed and the stability of the algorithm was also improved.
IWA recommended model parameters cannot be used to run ASM1 for each and every activated sludge process based wastewater treatment plants as they are being operated under different environmental conditions. Thus, seven sensitive parameters of the ASM1 model were estimated by the NAIWO algorithm using the data obtained from a large and a small to medium scale treatment plant. In tracking the effluent parameters of those two actual wastewater treatment plants, the NAIWO algorithm achieved better prediction accuracy than the GA, BA, and IWO algorithms.

Author Contributions

X.D., Y.M. and X.W. conceived and designed the experiments; Y.M. performed the experiments; X.D., Y.M. and V.J. analyzed the data; X.D., Y.M., X.W. and V.J. wrote the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the National Natural Science Foundation of China (No. 61563032, No. 61763208, No. 61963025), the Natural Science Foundation of Gansu Province, China (No. 1506RJZA104, No. 2017GS10945), and the University Scientific Research Project of Gansu Province (No. 2015B-030).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In order to evaluate the performance of the proposed NAIWO algorithm, two simulation experiments were designed and the details of which are provided in this section. In the first part, nine well-known test functions were selected to verify the searching ability of the proposed NAIWO. Table A1 shows the basic information of the nine benchmark functions and the corresponding 2D perspectives are shown in Figure A1.
Table A1. Two-dimensional benchmark test functions considered in the simulations.
Table A1. Two-dimensional benchmark test functions considered in the simulations.
FunctionsNamesDimensionSolution RangeGlobal Minimum
F1Ackley2 x i [ 32.768 , 32.768 ] fmin = 0 at (0, 0)
F2Cross-in-Tray2 x i [ 20 , 20 ] fmin = −2.06261 at (±1.3491, ±1.3491)
F3Drop-wave2 x i [ 10 , 10 ] fmin = −1 at (0, 0)
F4Griewank2 x i [ 100 , 100 ] fmin = 0 at (0, 0)
F5Levy2 x i [ 10 , 10 ] fmin = 0 at (1, 1)
F6Rosenbrock2 x i [ 2.048 , 2.048 ] fmin = 0 at (1, 1)
F7Schaffer2 x i [ 20 , 20 ] fmin = 0 at (0, 0)
F8Schwefel2 x i [ 500 , 500 ] fmin = 0 at (420.9687, 420.9687)
F9Three-hump camel2 x i [ 10 , 10 ] fmin = 0 at (0, 0)
Figure A1. A two-dimensional perspective with contours of the benchmark test functions, including (a) Ackley function (F1), (b) Cross-in-Tray function (F2), (c) Drop-wave function (F3), (d) Griewank function (F4), (e) Levy function (F5), (f) Rosenbrock function (F6), (g) Schaffer function (F7), (h) Schwefel function (F8), and (i) Three-hump camel function (F9).
Figure A1. A two-dimensional perspective with contours of the benchmark test functions, including (a) Ackley function (F1), (b) Cross-in-Tray function (F2), (c) Drop-wave function (F3), (d) Griewank function (F4), (e) Levy function (F5), (f) Rosenbrock function (F6), (g) Schaffer function (F7), (h) Schwefel function (F8), and (i) Three-hump camel function (F9).
Water 12 02604 g0a1
Except for F6 and F9, all the other seven functions have a large number of local minimums. Even for F6 and F9, the global minimum of the function cannot be easily found when reaching the convergence point. These functions are widely used in the performance test of any newly proposed intelligent algorithm [25,26].
During the test, the initial search range ( σ initial ) of the IWO algorithm and the NAIWO algorithm is generally 1% of the definition domain and the final search range is σ final = 1 e - 5 . The initial population P0 = 25 and the maximum population P max = 100 ; Maximum number of iterations i t e r max = 300 . Where in the NAIWO algorithm, the determination parameters of the niche radius are a = 3, b = 1.05, and k = 0.6, Adaptive spatial diffusion parameter K = 5 and T = 10.
Figure A2 shows the fitness curves of the two algorithms. The maximum number of iterations selected in the test is 300, while the global search of the two algorithms is mainly reflected in the first 100 generations and the iterative process of the latter is mainly local search in order to obtain better convergence accuracy. Therefore, when generating the fitness curve, in order to better display the contrast, the first 100 generations of data were selected. The nine images in Figure A2 illustrate that the NAIWO algorithm proposed in this study has obvious advantages in both global convergence ability and convergence speed than the IWO.
In order to obtain more reliable results, each function was tested 30 times and the average value, minimum value, and variance are shown in Table A2. By comparing the minimum convergence values of the two algorithms in Table A2, it can be seen that in most cases, the NAIWO algorithm has comparable or higher convergence accuracy than the IWO algorithm. Meanwhile, the mean of convergence can better reflect the overall convergence of the algorithm. The analysis reflects that NAIWO has a higher probability of converging to the global minimum than IWO. The variance data shows that the stability of the NAIWO algorithm converging to the global optimal solution is much higher than IWO.
Figure A2. Comparison of the fitness curves between NAIWO and IWO. The maximum number of iterations selected in the test is 300. In order to better display the contrast, the first 100 generations are selected when generating the image. (a) Ackley function (F1), (b) Cross-in-Tray function (F2), (c) Drop-wave function (F3), (d) Griewank function (F4), (e) Levy function (F5), (f) Rosenbrock function (F6), (g) Schaffer function (F7), (h) Schwefel function (F8), and (i) Three-hump camel function (F9).
Figure A2. Comparison of the fitness curves between NAIWO and IWO. The maximum number of iterations selected in the test is 300. In order to better display the contrast, the first 100 generations are selected when generating the image. (a) Ackley function (F1), (b) Cross-in-Tray function (F2), (c) Drop-wave function (F3), (d) Griewank function (F4), (e) Levy function (F5), (f) Rosenbrock function (F6), (g) Schaffer function (F7), (h) Schwefel function (F8), and (i) Three-hump camel function (F9).
Water 12 02604 g0a2
Table A2. Comparison data of IWO and NAIWO algorithms.
Table A2. Comparison data of IWO and NAIWO algorithms.
FunctionsAlgorithmsMeanMinimumVariance
F1IWO0.960311.329 × 10−711.153
NAIWO1.702 × 10−61.4251× 10−75.9939 × 10−13
F2IWO−2.0511−2.06261.919× 10−3
NAIWO−2.0626−2.06261.2302 × 10−28
F3IWO−0.8402−11.983 × 10−2
NAIWO−1−17.125 × 10−24
F4IWO5.5380 × 10−24.6835 × 10−22.1 × 10−4
NAIWO1.9692 × 10−141.22125 × 10−155.3704 × 10−28
F5IWO1.8955.7394 × 10−135.95262
NAIWO2.125 × 10−121.1457 × 10−135.7893 × 10−24
F6IWO6.29 × 10−33.8735 × 10−144.1 × 10−4
NAIWO8.3774 × 10−137.91675 × 10−165.876 × 10−25
F7IWO6.474 × 10−305.0604 × 10−5
NAIWO6.6613 × 10−1702.0912 × 10−32
F8IWO68.4312.5455 × 10−57463.60
NAIWO7.89592.5455 × 10−5902.928
F9IWO3.98 × 10−21.66 × 10−151.066 × 10−2
NAIWO1.32 × 10−137.04 × 10−161.229 × 10−26
To further verify the optimization ability of the proposed NAIWO algorithm, two well-known optimization algorithms, genetic algorithm (GA) and bat algorithm (BA), are compared using four higher-dimensional benchmark functions in this section. Table A3 shows the details of the four test functions. The number of iterations is 500 and the selection methods of the value of other parameters are the same as above.
Table A3. High dimensional benchmark test functions considered in the simulations.
Table A3. High dimensional benchmark test functions considered in the simulations.
Functions NamesDimensionEquationsGlobal Minimum
F10Griewank3 f ( x ) = i = 1 n x i 2 4000 i = 1 n cos ( x i i ) + 1 , 600 x i 600 fmin = 0 at (0,0,0)
F11Rastirgin3 f ( x ) = i = 1 n [ x i 2 10 cos ( 2 π x i ) ] + 10 n , 5.12 x i 5.12 fmin = 0 at (0,0,0)
F12Colville4 f ( x ) = 100 ( x 1 2 x 2 ) 2 + ( x 1 1 ) 2 + ( x 3 1 ) 2 + 90 ( x 3 2 x 4 ) 2 + 10.1 ( ( x 2 1 ) 2 + ( x 4 1 ) 2 ) + 19.8 ( x 2 1 ) ( x 4 1 ) , 10 < x i < 10 fmin = 0 at (1,1,1,1)
F13Powell4 f ( x ) = ( x 1 + 10 x 2 ) 2 + 5 ( x 3 x 4 ) 2 + ( x 2 2 x 3 ) 4 + 10 ( x 1 x 4 ) 4 , 4 < x i < 5 fmin = 0 at (0,0,0,0)
Figure A3 shows the fitness curve of the iterative process of the four algorithms, which shows that in the process of finding the minimum value of F10 and F11 functions, the NAIWO algorithm has the most outstanding comprehensive convergence ability and fastest convergence speed. BA also has a faster convergence speed, but its poor convergence stability makes it possible to converge to a local minimum with a certain probability.
The statistical information of the 30 repeated simulations shown in Table A4 also confirms the conclusions made above. In the 30 runs of the two functions F10 and F11, the minimum convergence accuracy of the BA algorithm is the highest of the four algorithms, but its average value and variance value are higher than that of the NAIWO algorithm by one or more orders of magnitude, indicating that its convergence probability and stability are far inferior to the NAIWO algorithm. The comparison of the statistical data of F10–F13 shows that the NAIWO algorithm has greater advantages in stability and convergence accuracies compared with the other three algorithms, which shows the effectiveness of the proposed algorithm.
Figure A3. Comparison of the fitness curves between NAIWO, IWO, GA, and BA. The maximum number of iterations selected in the test is 500. (a) Griewank function (F10), (b) Rastirgin function (F11), (c) Colville function (F12), and (d) Powell function (F13).
Figure A3. Comparison of the fitness curves between NAIWO, IWO, GA, and BA. The maximum number of iterations selected in the test is 500. (a) Griewank function (F10), (b) Rastirgin function (F11), (c) Colville function (F12), and (d) Powell function (F13).
Water 12 02604 g0a3
Table A4. Statistical comparison data of IWO, GA, BA, and NAIWO algorithms.
Table A4. Statistical comparison data of IWO, GA, BA, and NAIWO algorithms.
FunctionsAlgorithmsMeanMinimumVariance
F10IWO1.459272.58225 × 10−101.27896
GA0.423610.023780.17285
BA0.431152.70894 × 10−120.25147
NAIWO2.3897 × 10−96.7921 × 10−113.5773 × 10−18
F11IWO0.188750.007400.01780
GA0.352220.081010.01243
BA0.010528.88178 × 10−140.00023
NAIWO0.002631.61982 × 10−131.71477 × 10−5
F12IWO0.256761.02182 × 10−100.39066
GA0.196000.011250.02123
BA1.779760.014743.29279
NAIWO0.000114.29174 × 10−111.66852 × 10−7
F13IWO1.73784 × 10−71.98073 × 10−123.33946 × 10−13
GA0.086090.002254.89659 × 10−3
BA8.79407 × 10−57.17666 × 10−65.63957 × 10−9
NAIWO1.63239 × 10−91.12409 × 10−121.03073 × 10−17

Appendix B

ASM1 contains 13 substrate components (seven soluble components and six insoluble components) as shown in Table A5. The influent concentrations of these 13 components reflects the degree of contamination of the influent to the wastewater treatment plant. All the 13 substrates have their reaction rates and the reaction rate equations are given in Table A6.
Table A5. The 13 components of ASM1.
Table A5. The 13 components of ASM1.
ComponentsDefinitionUnit
SISoluble inert organic matterg∙COD∙m−3
SSReadily biodegradable substrateg∙COD∙m−3
XIParticulate inert organic matterg∙COD∙m−3
XSSlowly biodegradable substrateg∙COD∙m−3
XB,HActive heterotrophic biomassg∙COD∙m−3
XB,AActive autotrophic biomassg∙COD∙m−3
XPParticulate product arising from biomass decayg∙COD∙m−3
SOOxygen (negative COD)g∙COD∙m−3
SNONitrate and nitrite nitrogeng∙N∙m−3
SNHNH4+ and NH3 nitrogeng∙N∙m−3
SNDSoluble biodegradable organic nitrogeng∙N∙m−3
XNDParticulate biodegradable organic nitrogeng∙N∙m−3
SALKAlkalinitymol∙m−3
Table A6. Reaction rate equations of the 13 components.
Table A6. Reaction rate equations of the 13 components.
ComponentsReactions
SI (i = 1) r 1 = 0
SS (i = 2) r 2 = 1 Y H ρ 1 1 Y H ρ 2 + ρ 7
XI (i = 3) r 3 = 0
XS (i = 4) r 4 = ( 1 f p ) ρ 4 + ( 1 f p ) ρ 5 ρ 7
XB,H (i = 5) r 5 = ρ 1 + ρ 2 ρ 4
XB,A (i = 6) r 6 = ρ 3 ρ 5
XP (i = 7) r 7 = ρ 4 + ρ 5
SO (i = 8) r 8 = 1 Y H Y H ρ 1 4.57 Y A Y A ρ 3
SNO (i = 9) r 9 = 1 Y H 2.86 Y H ρ 2 + 1 Y A ρ 3
SNH (i = 10) r 10 = i X B ρ 1 i X B ρ 2 ( i X B + 1 Y A ) ρ 3 + ρ 6
SND (i = 11) r 11 = ρ 6 + ρ 8
XND (i = 12) r 12 = ( i X B f p i X P ) ρ 4 + ( i X B f p i X P ) ρ 5 ρ 8
SALK (i = 13) r 13 = i X B 14 ρ 1 + ( 1 Y H 14 2.86 Y H i X B 14 ) ρ 2 + ( i X B 14 + 1 7 Y A ) ρ 3 + 1 14 ρ 8
where ρ 1 to ρ 8 are the growth, decay, and hydrolysis of microorganisms in wastewater listed in Table A7. These 21 reaction processes, shown in Table A6 and Table A7, represent the basic processes involved in the biological treatment of activated sludge-based wastewater treatment plants.
Table A7. Reaction rate of eight sub-reaction processes.
Table A7. Reaction rate of eight sub-reaction processes.
ProcessesDescriptions
ρ 1 μ H ( S S K S + S S ) ( S O K O , H + S O ) X B , H
ρ 2 μ H ( S S K S + S S ) ( K O , H K O , H + S O ) ( S N O K N O + S N O ) η g X B , H
ρ 3 μ A ( S N H K N H + S N H ) ( S O K O , A + S O ) X B , A
ρ 4 b H X B , H
ρ 5 b A X B , A
ρ 6 k a B N D X B , H
ρ 7 k h X S / X B , H K X + ( X S / X B , H ) [ η h ( S O K O , H + S O ) + η h ( K O , H K O , H + S O ) ( S N D K N O + S N O ) ] X B , H
ρ 8 ρ 7 ( X N D / X S )
There are five stoichiometric parameters of biochemical reactions and 14 kinetic parameters involved in Table A6 and Table A7, respectively. The correctness of these 19 parameters ensures the accuracy of the ASM1 model in simulating the performance of actual wastewater treatment plants.

References

  1. Henze, M.; Grady, L., Jr.; Gujer, W.; Marais, G.; Matsuo, T. Activated Sludge Model No. 1; IAWPRC: London, UK, 1987; Available online: https://www.researchgate.net/publication/243624144_Activated_Sludge_Model_No_1 (accessed on 15 September 2020).
  2. Al Madany, A.M.; El-Seddik, M.M.; Abdallah, K.Z. Extended Activated Sludge Model No. 1 with Floc and Biofilm Diffusion for Organic and Nutrient Removal. J. Environ. Eng. 2020, 146, 04020008. [Google Scholar] [CrossRef]
  3. Salles, N.A.; de Souza, T.S.O. Activated Sludge Model No. 1 (ASM1) applicability for simulation of sanitary sewage and landfill leachate co-treatment in aerated lagoons. Eng. Sanit. Ambient. 2020, 25, 293–301. [Google Scholar] [CrossRef] [Green Version]
  4. Hauduc, H.; Gillot, S.; Rieger, L.; Ohtsuki, T.; Shaw, A.; Takács, I.; Winkler, S. Activated sludge modelling in practice: An international survey. Water Sci. Technol. 2009, 60, 1943–1951. [Google Scholar] [CrossRef] [PubMed]
  5. Van Loosdrecht, M.C.M.; Lopez-Vazquez, C.M.; Meijer, S.; Hooijmans, C.M.; Brdjanovic, D. Twenty-five years of ASM1: Past, present and future of wastewater treatment modelling. J. Hydroinform. 2015, 17, 697–718. [Google Scholar] [CrossRef] [Green Version]
  6. Hauduc, H.; Rieger, L.; Oehmen, A.; Van Loosdrecht, M.C.M.; Comeau, Y.; Héduit, A.; Vanrolleghem, P.; Gillot, S. Critical review of activated sludge modeling: State of process knowledge, modeling concepts, and limitations. Biotechnol. Bioeng. 2012, 110, 24–46. [Google Scholar] [CrossRef]
  7. Bhuyan, S.; Koelliker, J.; Marzen, L.; Harrington, J. An integrated approach for water quality assessment of a Kansas watershed. Environ. Model. Softw. 2003, 18, 473–484. [Google Scholar] [CrossRef]
  8. Yeh, W.W.-G.; Yoon, Y.S.; Lee, K.S. Aquifer parameter identification with kriging and optimum parameterization. Water Resour. Res. 1983, 19, 225–233. [Google Scholar] [CrossRef]
  9. Ratto, M.; Tarantola, S.; Saltelli, A. Sensitivity analysis in model calibration: GSA-GLUE approach. Comput. Phys. Commun. 2001, 136, 212–224. [Google Scholar] [CrossRef]
  10. Omlin, M.; Reichert, P. A comparison of techniques for the estimation of model prediction uncertainty. Ecol. Model. 1999, 115, 45–59. [Google Scholar] [CrossRef]
  11. Beven, K.; Binley, A. The future of distributed models: Model calibration and uncertainty prediction. Hydrol. Process. 1992, 6, 279–298. [Google Scholar] [CrossRef]
  12. Cho, J.; Sung, K.S.; Ha, S.R. A river water quality management model for optimising regional wastewater treatment using a genetic algorithm. J. Environ. Manag. 2004, 73, 229–242. [Google Scholar] [CrossRef] [PubMed]
  13. Miró, A.; Pozo, C.; Guillén-Gosálbez, G.; Egea, J.A.; Jiménez, L. Deterministic global optimization algorithm based on outer approximation for the parameter estimation of nonlinear dynamic biological systems. BMC Bioinform. 2012, 13, 90. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Gandomi, A.H.; Yang, X.S.; Talatahari, S.; Alavi, A.H. (Eds.) Metaheuristic Applications in Structures and Infrastructures; Elsevier Science Publishers B.V.: Amsterdam, The Netherlands, 2013. [Google Scholar]
  15. Kim, S.; Lee, H.; Kim, J.; Kim, C.; Ko, J.; Woo, H.; Kim, S. Genetic algorithms for the application of Activated Sludge Model No. 1. Water Sci. Technol. 2002, 45, 405–411. [Google Scholar] [CrossRef] [PubMed]
  16. Du, X.; Wang, J.; Jegatheesan, V.; Shi, G. Parameter estimation of activated sludge process based on an improved cuckoo search algorithm. Bioresour. Technol. 2018, 249, 447–456. [Google Scholar] [CrossRef]
  17. Mehrabian, A.; Lucas, C. A novel numerical optimization algorithm inspired from weed colonization. Ecol. Inform. 2006, 1, 355–366. [Google Scholar] [CrossRef]
  18. Wu, H.; Liu, C.; Xie, X. Pattern synthesis of planar antenna arrays based on invasive weed optimization algorihm. J. Naval Univ. Eng. 2015, 27, 16–19. [Google Scholar]
  19. Liu, Y.; Jiao, Y.-C.; Zhang, Y.-M. A novel hybrid invasive weed optimization algorithm for pattern synthesis of array antennas. Int. J. RF Microw. Comput. Aided Eng. 2014, 25, 154–163. [Google Scholar] [CrossRef]
  20. Duan, M.; Yuan, C. Fault diagnosis of nuclear power plant based on invasive weed optimization algorithm. At. Energy Sci. Technol. 2015, 49, 719–724. [Google Scholar]
  21. Cuevas, E.; Rodríguez, A.; Valdivia, A.; Zaldívar, D.; Pérez, M. A hybrid evolutionary approach based on the invasive weed optimization and estimation distribution algorithms. Soft Comput. 2019, 23, 13627–13668. [Google Scholar] [CrossRef]
  22. Martins, T.M.; Neves, R.F. Applying genetic algorithms with speciation for optimization of grid template pattern detection in financial markets. Expert Syst. Appl. 2020, 147, 113191. [Google Scholar] [CrossRef]
  23. Zhou, S.J.; Li, H.M.; Gao, T.B. Combination of Ant Algorithm and niche Genetic Algorithm. In Proceedings of First International Conference of Modelling and Simulation, Volume IV: Modelling and Simulation in Business, Management, Economic and Finance; Yu, H.G., Jiang, Y., Eds.; World Academic Union-World Academic Press: Liverpool, UK, 2008; pp. 169–172. [Google Scholar]
  24. Mirjalili, S.; Mirjalili, S.M.; Lewis, A. Grey Wolf Optimizer. Adv. Eng. Softw. 2014, 69, 46–61. [Google Scholar] [CrossRef] [Green Version]
  25. Ju, X.; Peng, D. Kinetic parameters sensitivity analysis of model ASM1. Environ. Sci. Technol. 2010, 33, 312–314. [Google Scholar]
  26. Faramarzi, A.; Heidarinejad, M.; Stephens, B.; Mirjalili, S. Equilibrium optimizer: A novel optimization algorithm. Knowl. Based Syst. 2020, 191, 105190. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of activated sludge model 1 (ASM1) parameter estimation using the niche-based adaptive invasive weed optimization (NAIWO) algorithm.
Figure 1. Schematic diagram of activated sludge model 1 (ASM1) parameter estimation using the niche-based adaptive invasive weed optimization (NAIWO) algorithm.
Water 12 02604 g001
Figure 2. Basic processing blocks of the activated sludge process based wastewater treatment plants.
Figure 2. Basic processing blocks of the activated sludge process based wastewater treatment plants.
Water 12 02604 g002
Figure 3. Flow chart of the NAIWO algorithm used in ASM1 parameter estimation.
Figure 3. Flow chart of the NAIWO algorithm used in ASM1 parameter estimation.
Water 12 02604 g003
Figure 4. The fitness curves of four algorithms produced during the parameter estimation for (a) Pingliang City Wastewater Treatment Plant (PC-WWTP) and (b) WC-WWTP.
Figure 4. The fitness curves of four algorithms produced during the parameter estimation for (a) Pingliang City Wastewater Treatment Plant (PC-WWTP) and (b) WC-WWTP.
Water 12 02604 g004
Figure 5. Comparison of predicted and measured values of effluent water quality parameters for (a) Ammonia nitrogen (SNH), (b) nitrate nitrogen (SNO), (c) soluble rapidly biodegradable organic matter (SS), and (d) insoluble slowly degradable organic matter (XS) in PC-WWTP.
Figure 5. Comparison of predicted and measured values of effluent water quality parameters for (a) Ammonia nitrogen (SNH), (b) nitrate nitrogen (SNO), (c) soluble rapidly biodegradable organic matter (SS), and (d) insoluble slowly degradable organic matter (XS) in PC-WWTP.
Water 12 02604 g005
Figure 6. Comparison of predicted and measured values of effluent water quality parameters for (a) Ammonia nitrogen (SNH), (b) nitrate nitrogen (SNO), (c) soluble rapidly biodegradable organic matter (SS), and (d) insoluble slowly degradable organic matter (XS) in WC-WWTP.
Figure 6. Comparison of predicted and measured values of effluent water quality parameters for (a) Ammonia nitrogen (SNH), (b) nitrate nitrogen (SNO), (c) soluble rapidly biodegradable organic matter (SS), and (d) insoluble slowly degradable organic matter (XS) in WC-WWTP.
Water 12 02604 g006
Table 1. Nineteen parameters of ASM1.
Table 1. Nineteen parameters of ASM1.
SymbolDescriptionValue ++Unit
5 stoichiometric parameters
YAYield coefficient of autotrophic bacteria0.24g cell COD formed·(g COD oxidized)−1
YHHeterotrophic yield coefficient0.67g cell COD formed·(g N oxidized)−1
fpProportion of inert particles in microorganisms0.08dimensionless
ixbProportion of nitrogen content in microbial cells0.086g N·(g COD)−1 in biomass
ixpProportion of nitrogen content in microbial products0.06g N·(g COD)−1 in particulate products
14 kinetic parameters
μ H Maximum specific growth rate coefficient of heterotrophic bacteria6.0d−1
KSHalf saturation coefficient of heterotrophic bacteria20.0g COD·m−3
KO,HOxygen half-saturation coefficient of heterotrophic bacteria0.2g COD·m−3
KNONitrate half-saturation coefficient of denitrifying bacteria0.5g NO3-N·m−3
bHAttenuation coefficient of heterotrophic bacteria0.62d−1
η g Hypoxia correction factor for heterotrophic bacteria under hypoxia0.8dimensionless
η H Correction factor for hydrolysis under anoxic conditions0.4dimensionless
khMaximum specific hydrolysis rate3.0g slowly biodegradable COD·(g cell COD·d)−1
KXHalf saturation coefficient of slowly biodegradable substrate hydrolysis0.03g slowly biodegradable COD·(g cell COD)−1
μ A Maximum specific growth rate of autotrophic bacteria0.80d−1
KNHAmmonia half-saturation coefficient of autotrophic bacteria1.0g NH3-N·m−3
KO,AOxygen half saturation coefficient of autotrophic bacteria0.4d−1
kaDissolved organic ammonia nitrating rate0.08g COD·m−3
bAAutotrophic bacteria specific decay rate0.05m3·(g COD·d)−1
++ The values are for 20 °C, extracted from [1].
Table 2. Estimated results of seven parameters of ASM1 for PC-WWTP.
Table 2. Estimated results of seven parameters of ASM1 for PC-WWTP.
ParameterRangeInternational Water Association (IWA) Recommendation (20 °C) Values Estimated by Selected Algorithms
Genetic Algorithm (GA)Bat Algorithm (BA)Invasive Weed Optimization (IWO)NAIWO
YH0.4–0.90.670.66630.67200.66490.6690
bH0.3–1.20.620.34480.37280.44130.4521
µH2–1862.63882.89413.21014.9771
µA0.3–2.40.80.72930.92430.92250.6733
KOA0.13–1.20.40.21890.60111.10640.5139
KNH0.3–311.92322.08191.24791.0908
KS10–402012.950213.646313.101921.8191
Table 3. Estimated results of seven parameters of ASM1 for Wushan County Wastewater Treatment Plant (WC-WWTP).
Table 3. Estimated results of seven parameters of ASM1 for Wushan County Wastewater Treatment Plant (WC-WWTP).
ParameterRangeIWA Recommendation (20 °C) Values Estimated by Selected Algorithms
GABAIWONAIWO
YH0.4–0.90.670.76300.74470.72470.7510
bH0.3–1.20.620.49530.54800.54820.5891
µH2–1864.96085.74968.50816.2899
µA0.3–2.40.82.24811.44341.99180.9003
KOA0.13–1.20.41.12050.56001.18550.3851
KNH0.3–312.82311.93521.97321.0687
KS10–402018.202416.898628.136918.8017

Share and Cite

MDPI and ACS Style

Du, X.; Ma, Y.; Wei, X.; Jegatheesan, V. Optimal Parameter Estimation in Activated Sludge Process Based Wastewater Treatment Practice. Water 2020, 12, 2604. https://doi.org/10.3390/w12092604

AMA Style

Du X, Ma Y, Wei X, Jegatheesan V. Optimal Parameter Estimation in Activated Sludge Process Based Wastewater Treatment Practice. Water. 2020; 12(9):2604. https://doi.org/10.3390/w12092604

Chicago/Turabian Style

Du, Xianjun, Yue Ma, Xueqin Wei, and Veeriah Jegatheesan. 2020. "Optimal Parameter Estimation in Activated Sludge Process Based Wastewater Treatment Practice" Water 12, no. 9: 2604. https://doi.org/10.3390/w12092604

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