Research Article  Open Access
Yasaman Sanayei, Naz Chaibakhsh, Ali Chaibakhsh, Ali Reza Pendashteh, Norli Ismail, Tjoon Tow Teng, "LongTerm Prediction of Biological Wastewater Treatment Process Behavior via WienerLaguerre Network Model", International Journal of Chemical Engineering, vol. 2014, Article ID 248450, 7 pages, 2014. https://doi.org/10.1155/2014/248450
LongTerm Prediction of Biological Wastewater Treatment Process Behavior via WienerLaguerre Network Model
Abstract
A WienerLaguerre model with artificial neural network (ANN) as its nonlinear static part was employed to describe the dynamic behavior of a sequencing batch reactor (SBR) used for the treatment of dyecontaining wastewater. The model was developed based on the experimental data obtained from the treatment of an effluent containing a reactive textile azo dye, Cibacron yellow FN2R, by Sphingomonas paucimobilis bacterium. The influent COD, MLVSS, and reaction time were selected as the process inputs and the effluent COD and BOD as the process outputs. The best possible result for the discrete pole parameter was . In order to adjust the parameters of ANN, the LevenbergMarquardt (LM) algorithm was employed. The results predicted by the model were compared to the experimental data and showed a high correlation with and a low mean absolute error (MAE). The results from this study reveal that the developed model is accurate and efficacious in predicting COD and BOD parameters of the dyecontaining wastewater treated by SBR. The proposed modeling approach can be applied to other industrial wastewater treatment systems to predict effluent characteristics.
1. Introduction
The reactive dyecontaining effluents from dye manufacturing and application industries can cause serious environment pollution due to the toxicity and slow degradation of dyes [1]. In addition, the presence of dyes in water is highly visible and affects water transparency and aesthetics even in low concentrations. Therefore, the effluents must be treated before being released into the environment.
In recent years, researchers have shown interests in biological treatment of wastewaters with high concentrations of dyes [2, 3]. Treatment of these polluted wastewaters requires high effectiveness and low cost dye removal processes [4]. Sequencing batch reactor (SBR) is a promising biological system for treating dyecontaining wastewaters [5, 6]. This system is cost efficient and flexible to handle different feed characteristics. Furthermore, its operation is easier than other biological methods [1]. However, the SBR process, like other biological processes, is highly nonlinear, time varying, and subject to significant disturbances [7]. Modeling the treatment process can provide better understanding, design, operation, and control of the process [8].
The ability of artificial neural networks (ANNs) in blackbox modeling of nonlinear systems with complicated structure has made them the most popular tool for modeling of biological processes [9]. In recent years, recurrent neural networks are developed based on common NARX (Nonlinear AutoRegressive eXogenous model) topology in order to characterize the behaviors of dynamic systems. In such models, the past terms of input and output signals are used as the network inputs [10]. In literature, different applications of NARX networks have been reported for monitoring [11, 12], controlling [13], optimization [14], and simulation [15, 16] of wastewater treatment process plants. In spite of good performances of NARX networks, they are still suffering from some drawbacks such as the need for prior terms information of input and output variables. Furthermore, the presence of output feedback may cause prediction errors to deteriorate the model performances [10].
Wienertype models are a class of blockoriented representation of nonlinear systems, where a dynamic linear part is cascaded by a static nonlinear part [17]. In recent years, Wienertype models have been successfully applied to identify nonlinear dynamics in biological processes [18–20]. Employing standard multilayer perceptron (MLP) networks as nonlinear mapping function in Wiener models structure is proposed as an alternative solution to develop nonlinear models [21]. The linear part of Wiener models is often chosen to be as formal autoregressive with exogenous input (ARX) models. In this case, in order to attain high performance, a large amount of training data for estimating the order of model, time delays, and parameters would be required. A big obstacle is the high cost of preparing training data, where special instruments and data acquisition equipment for online system monitoring are required.
A combination of orthogonal functions as the linear dynamic part and artificial neural networks as the memoryless nonlinear part of Wiener models was proposed to deal with this problem [22, 23]. Laguerre functions are the most common orthogonal functions used in Wiener model structure. Simple structure and capability to give accurate description of systems with rather small number of tunable parameters are the reasons they have been widely employed for linear and nonlinear system identification [24].
In this study, a WienerLaguerre model is employed to describe the nonlinear behavior of an SBR used for dye wastewater treatment. The proposed method has not been previously used for biological treatment modeling of dyecontaining wastewaters. An aerobic bacterium, Sphingomonas paucimobilis, was used for decolorizing of the effluent containing a reactive textile azo dye, Cibacron yellow FN2R. S. paucimobilis has been shown previously to be efficient for degradation of azo and triphenylmethane dyes [25, 26].
2. Material and Methods
2.1. Microorganism
Sphingomonas paucimobilis was isolated from a closed drainage system located in Perai industrial area, Perai, Penang, Malaysia, where several heavy industries were located. The microorganism was grown on slant agar at 35°C for 2 days under static conditions.
2.2. Media Composition
A synthetic dyecontaining wastewater was used in this study. Seed culture medium consisted of distilled water 1000 mL, powdered Cibacron yellow FN2R as the main carbon source, urea as the nitrogen source, and K_{2}HPO_{4} and KH_{2}PO_{4} as the phosphorous sources and the nutrients. The composition of the wastewater gave a C/N/P ratio of approximately 100/10/1 by adding NH_{4}Cl and KH_{2}PO_{4}. The initial pH was adjusted to 7 using NaOH.
2.3. Preparation of the Synthetic Dye Wastewater
Cibacron yellow FN2R is one of the main reactive dyes that are used in textile industries. Reactive Cibacron FN2R was purchased from SigmaAldrich. The dye solution of various concentrations 250, 500, 1000, and 1500 mg/L were treated in four reactors with successively increasing of dye concentration and were analyzed for COD, BOD, and MLVSS everyday according to the standard methods [27].
2.4. Experimental Setup
In this study, four cylindrical Plexiglass reactors with 14 cm diameter and 46 cm height were used. The working volume and influent flow rate were 1.6 L and 3.0 L/d, respectively. Four air pumps and four mixers were used for continuous aeration and mixing. The speed of impeller was adjusted to 300 rpm to maintain the dissolved oxygen concentration at DO > 2 ppm. After acclimatization period, different concentrations of the dye and 2 mL of bacterium (S. paucimobilis) were added to each reactor. The pH of the mixture was adjusted to 7 after stirring at a constant temperature (35°C) for a certain period.
2.5. InputOutput Selection
In order to develop accurate and reliable models, appropriate data sets should be collected from the variables, which are governing the behavior of the process. A set of variables including the influent COD, MLVSS, and reaction time as the process inputs and the effluent COD and BOD as the process output were prepared. Two sets of data including 147 data points (in 147 days) and 90 data points (in 90 days) are chosen for training and validation of the models, respectively. As a result, the input/output vectors would be as For increasing the performance of modeling process, the recorded data are normalized in the range of by dividing on the maximum value of each datum column.
3. Model Structure and System Identification
Based on the obtained experimental data, appropriate nonlinear models were developed to predict the concentrations of COD and BOD in the effluent.
Wienertype model is a wellknown blockoriented representation of nonlinear systems, in which a linear dynamic part is cascaded with static nonlinearities. As it is shown in Figure 1, the sequences of input variable are mapped into the intermediate variable through the transfer function , where the model output is estimated throughout the nonlinear mapping function . Many different options are available for the linear part and nonlinear mapping functions.
Here, the considered Wiener model is a combination of Laguerre basis filters as the linear part and a feedforward NN for nonlinear mapping. In this case, without loss of generality by considering a SISO (singleinput and singleoutput) system, the nonlinear dynamics can be represented as a discretetime nonlinear mapping on filtered input terms, where denotes the Laguerre filters. This structure is a NN realization of the WienerLaguerre model [21].
As it is presented in (1), three inputs are chosen as the dominant variables on the effluent COD and BOD. In this case, each process model can be considered as a multiinput singleoutput (MISO) model, as it is presented in Figure 2.
3.1. Laguerre Filters
Laguerre functions can be represented as a set of discretetime transfer functions in zdomain, where is the pole parameter. The dominant pole “” determines the rate of exponential decay for Laguerre functions responses, which can be captured through optimization or by experiments [28]. The first section of is a firstorder lowpass filter, which is followed by identical allpass filters. The orthonormal bases are formed by a set of Laguerre functions on [0,), which is defined on zdomain. Consider For physical systems, the limited order of Laguerre stages can be employed to approximate system dynamics. In this case, for a given truncated order “,” the transfer function for a SISO linear discretetime system is obtained as follows: in which and are input and output of the system and are the Laguerre coefficients [28]. By considering , will turn to regular delay operators and to the usual FIR (finite input response) model. Introducing as the output of Laguerre filters, Thus, for a linear system, the output of an “aorder” Laguerre network is obtained as the linear combination of the Laguerre filters output as The linear combination of the Laguerre filters’ output can be replaced by a multilayer perceptrons (MLP) as the nonlinear mapping function in order to deal with nonlinearities. In order to increase the performances of Laguerre networks, it is required that the pole parameters and the order of the Laguerre filters have to be estimated.
Developing a first order linear model for the process could be helpful to find appropriate process time constant. For this aim, the influent and the effluent CODs can be considered as the model input and output signals. We have in which and are the process gain and time constant. The parameters of linear model can be captured using the MATLAB System Identification Toolbox (Ver. 7.4.2). The best possible result for the linear system pole was /, as the model output fit on the real data was about 52.1%. The equivalent discrete pole parameter can be captured by estimating in the zdomain, where the sampling time is equal to 1. This results in .
The step responses of Laguerre filters (normalized) are shown in Figure 3. It is noted that by increasing the order of Laguerre filters, the structure of the NN part would be more complicated and therefore computational efforts for the nonlinear part training would increase. A major advantage of Laguerre network is its ability to describe the behavior of nonlinear systems by hightruncated order models. Choosing the order of filter is an appropriate tradeoff between the complexity and the accuracy of models. The pole parameters and filter order are chosen to be 0.44 and 9, respectively, for all input variables.
3.2. Neural Network
A multilayer feedforward neural network (FFNN) is considered as the nonlinear part of the proposed models. A FFNN is developed by a number of cascaded layers that are interconnected by weight coefficients to the neighboring layers nodes. Neurons are the fundamental part of NN and responsible for information processing, which are consisting of weighting coefficients, activation functions, and biases. The structure of FFNN is presented in Figure 4. It is noted that, in this work, no bias connection was considered.
Generally, a FFNN may have one or more hidden layers; however, a network with a single hidden layer could approximate many systems with an acceptable degree of accuracy. In many cases, choosing the appropriate number of neurons in the hidden layer and the type of the activation functions is known as the dominant parameters on the model performances. The number of neurons in the input/output layers is chosen with respect to the number of input/output variables of the investigated system, which are 27 and 1 for input and output layers, respectively. For a single hiddenlayer NN, it is recommended that number of hidden neurons () be chosen with respect to the geometric pyramid rule proposed as where and are the number of network inputs and outputs, respectively, and is multiplication factor. The value of should be selected in the range of depends on complexity of system [29]. As a result, by considering , the number of neurons in the hidden layer would be equal to 10.
The best activation functions may often be selected through trial and error. For the hidden layer, activation functions are chosen to be hyperbolic tangent sigmoid function and the output layer chosen as linear transfer function.
In order to adjust the parameters of NN, the secondorder derivativebased LevenbergMarquardt (LM) algorithm was employed. The LM is an accelerated neural networks training algorithm, which is rather used for adjusting the parameters of moderatesized MLP. The model weight parameters (W) are tuned through iteration by the following equation: where is the Jacobian matrix that is captured by first derivation of network errors with respect to its parameters and μ is some nonnegative value known as the learning parameter [30]. The objective function for training process was considered as mean square of error (MSE), where and are the target and actual outputs for the th pattern, respectively, and is the total number of training patterns. The error vector is calculated as follows: The MATLAB Neural Network Toolbox (ver. 7.0.1) was employed to perform the process of model training.
4. Simulations and Results
The proposed modeling approach was applied to the experimental data. The performance of the models was measured by the coefficient of determination () and mean absolute error (MAE) between the predicted values of the model and the experimental values, which were calculated by (13) and (14), respectively. Consider where is the average of over the data, and and are the th target and predicted responses, respectively.
In Figure 5(a), the responses of the developed model for prediction of BOD are presented. In addition, the correlation between the experimental and predicted values is presented in Figure 5(b). Figure 6(a) shows the responses of the developed model for prediction of COD, whereas the correlation between the experimental and predicted values of COD are depicted in Figure 6(b).
(a)
(b)
(a)
(b)
The results indicate the accuracy and ability of the proposed model to predict the selected parameters perfectly.
The prediction accuracy of the developed models was also evaluated by performing a comparison between the responses of the models and the actual values. For this propose, the error functions were calculated, where the error is defined as the difference between the predicted values by the models and the experimental data. Here, the upper bound error (Max()), lower bound error (Min()), MAE, and are calculated and presented in Table 1. The results show a small deviation between the models predicted values and the experimental data. A significantly high and a low MAE obtained indicate the superior data fitting and prediction capability of the developed model.

BOD and COD are two major parameters for examining the quality of discharged wastewater. Their measurements need significant time and commitment to make proper adjustments in the wastewater treatment process. The model proposed has the ability to provide good generalization performance in capturing nonlinear relationships between parameters to estimate COD and BOD of the treated wastewater.
5. Conclusion
The present study aimed to explore the potential of the WienerLaguerre network model in estimation of BOD and COD of the output stream of a SBR system. A multilayer feedforward neural network was considered as the nonlinear part of the proposed model. The performance of the model was found to be reasonably good with a high (0.99) and a small deviation between the predicted and experimental values. The proposed model can be used as a flexible alternative to the firstorder models commonly used for longterm prediction of BOD and COD parameters.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
References
 O. J. Hao, H. Kim, and P. C. Chiang, “Decolorization of wastewater,” Critical Reviews in Environmental Science and Technology, vol. 30, no. 4, pp. 449–505, 2000. View at: Google Scholar
 B. Saba, A. Khalid, A. Nazir, H. Kanwal, and T. Mahmood, “Reactive black5 azo dye treatment in suspended and attach growth sequencing batch bioreactor using different cosubstrates,” International Biodeterioration & Biodegradation, vol. 85, pp. 556–562, 2013. View at: Publisher Site  Google Scholar
 I. Khouni, B. Marrot, and R. B. Amar, “Treatment of reconstituted textile wastewater containing a reactive dye in an aerobic sequencing batch reactor using a novel bacterial consortium,” Separation and Purification Technology, vol. 87, pp. 110–119, 2012. View at: Publisher Site  Google Scholar
 N. Koprivanac, J. JovanovićKolar, G. Bosanac, and J. Meixner, “Studies on wastewater decolorization by the precipitation/flocculation process,” Microchemical Journal, vol. 46, no. 3, pp. 379–384, 1992. View at: Google Scholar
 P. E. Lim and C. C. Er, “Treatment of dyecontaining wastewater by sequencing batch reactor with powdered activated carbon addition,” Toxicological and Environmental Chemistry, vol. 75, no. 12, pp. 75–87, 2000. View at: Google Scholar
 S. V. Mohan, N. C. Rao, and P. N. Sarma, “Simulated acid azo dye wastewater treatment using suspended growth configured sequencing batch reactor (SBR) under anoxicaerobicanoxic microenvironment,” Applied Ecology and Environmental Research, vol. 7, no. 1, pp. 25–34, 2009. View at: Google Scholar
 C. K. Yoo, K. Villez, I. B. Lee, C. Rosén, and P. A. Vanrolleghem, “Multimodel statistical process monitoring and diagnosis of a sequencing batch reactor,” Biotechnology and Bioengineering, vol. 96, no. 4, pp. 687–701, 2007. View at: Publisher Site  Google Scholar
 D. Güçlü and Ş. Dursun, “Artificial neural network modelling of a largescale wastewater treatment plant operation,” Bioprocess and Biosystems Engineering, vol. 33, no. 9, pp. 1051–1058, 2010. View at: Publisher Site  Google Scholar
 J. S. Almeida, “Predictive nonlinear modeling of complex data by artificial neural networks,” Current Opinion in Biotechnology, vol. 13, no. 1, pp. 72–76, 2002. View at: Publisher Site  Google Scholar
 J. M. P. Menezes Jr. and G. A. Barreto, “Longterm time series prediction with the NARX network: an empirical evaluation,” Neurocomputing, vol. 71, no. 16–18, pp. 3335–3343, 2008. View at: Publisher Site  Google Scholar
 S. H. Hong, M. W. Lee, D. S. Lee, and J. M. Park, “Monitoring of sequencing batch reactor for nitrogen and phosphorus removal using neural networks,” Biochemical Engineering Journal, vol. 35, no. 3, pp. 365–370, 2007. View at: Publisher Site  Google Scholar
 H. G. Han, Q. L. Chen, and J. F. Qiao, “An efficient selforganizing RBF neural network for water quality prediction,” Neural Networks, vol. 24, no. 7, pp. 717–725, 2011. View at: Publisher Site  Google Scholar
 G. M. Zeng, X. S. Qin, L. He, G. H. Huang, H. L. Liu, and Y. P. Lin, “A neural network predictive control system for paper mill wastewater treatment,” Engineering Applications of Artificial Intelligence, vol. 16, no. 2, pp. 121–129, 2003. View at: Publisher Site  Google Scholar
 Y. E. Hamzaoui, J. A. Hernández, S. SilvaMartínez, A. Bassam, A. Álvarez, and C. LizamaBahena, “Optimal performance of COD removal during aqueous treatment of alazine and gesaprim commercial herbicides by direct and inverse neural network,” Desalination, vol. 277, no. 1–3, pp. 325–337, 2011. View at: Publisher Site  Google Scholar
 J. C. Chen, N. B. Chang, and W. K. Shieh, “Assessing wastewater reclamation potential by neural network model,” Engineering Applications of Artificial Intelligence, vol. 16, no. 2, pp. 149–157, 2003. View at: Publisher Site  Google Scholar
 B. Ráduly, K. V. Gernaey, A. G. Capodaglio, P. S. Mikkelsen, and M. Henze, “Artificial neural networks for rapid WWTP performance evaluation: methodology and case study,” Environmental Modelling and Software, vol. 22, no. 8, pp. 1208–1216, 2007. View at: Publisher Site  Google Scholar
 S. I. Biagiola and J. L. Figueroa, “Wiener and Hammerstein uncertain models identification,” Mathematics and Computers in Simulation, vol. 79, no. 11, pp. 3296–3313, 2009. View at: Publisher Site  Google Scholar
 V. Quaglini, F. Previdi, R. Contro, and S. Bittanti, “A discretetime nonlinear Wiener model for the relaxation of soft biological tissues,” Medical Engineering and Physics, vol. 24, no. 1, pp. 9–19, 2002. View at: Publisher Site  Google Scholar
 A. D. Kalafatis, L. Wang, and W. R. Cluett, “Linearizing feedforwardfeedback control of pH processes based on the Wiener model,” Journal of Process Control, vol. 15, no. 1, pp. 103–112, 2005. View at: Publisher Site  Google Scholar
 I. H. Song and H. K. Rhee, “Application of Wiener type predictive controller to the continuous solution polymerization reactor,” Studies in Surface Science and Catalysis, vol. 159, pp. 861–864, 2006. View at: Google Scholar
 M. Ławryńczuk, “Computationally efficient nonlinear predictive control based on neural Wiener models,” Neurocomputing, vol. 74, no. 1–3, pp. 401–417, 2010. View at: Publisher Site  Google Scholar
 S. Mahmoodi, J. Poshtan, M. R. JahedMotlagh, and A. Montazeri, “Nonlinear model predictive control of a pH neutralization process based on WienerLaguerre model,” Chemical Engineering Journal, vol. 146, no. 3, pp. 328–337, 2009. View at: Publisher Site  Google Scholar
 V. Arto, P. Hannu, and A. Halme, “Modeling of chromatographic separation process with WienerMLP representation,” Journal of Process Control, vol. 11, no. 5, pp. 443–458, 2001. View at: Publisher Site  Google Scholar
 B. Wahlberg and P. M. Mäkilä, “On approximation of stable linear dynamical systems using Laguerre and Kautz functions,” Automatica, vol. 32, no. 5, pp. 693–708, 1996. View at: Publisher Site  Google Scholar
 Y. Sanayei, N. Ismail, T. T. Teng, and N. Morad, “Biological treatment of reactive dye (Cibacron yellow FN_2R) by sequencing batch reactor performance,” Australian Journal of Basic and Applied Sciences, vol. 3, no. 4, pp. 4071–4077, 2009. View at: Google Scholar
 L. Ayed, S. Achour, E. Khelifi, A. Cheref, and A. Bakhrouf, “Use of active consortia of constructed ternary bacterial cultures via mixture design for Congo Red decolorization enhancement,” Chemical Engineering Journal, vol. 162, no. 2, pp. 495–502, 2010. View at: Publisher Site  Google Scholar
 APHA, Standard Methods for the Examination of Water and wasteWater, American Public Health Association, New York, NY, USA, 2005.
 M. A. MasnadiShirazi and M. Aleshams, “Laguerre discretetime filter design,” Computers and Electrical Engineering, vol. 29, no. 1, pp. 173–192, 2003. View at: Publisher Site  Google Scholar
 A. K. Palit and D. Popovic, Computational Intelligence in Time Series Forecasting: Theory and Engineering Applications, Springer, London, UK, 2005.
 J. S. R. Jang, C. T. Sun, and E. Mitutani, NeuroFuzzy and Soft Computing: A Computational Approach to Learning and Machine Intelligence, PrenticeHall, Upper Saddle River, NJ, USA, 1st edition, 1997.
Copyright
Copyright © 2014 Yasaman Sanayei et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.