Next Article in Journal
A Cooperative Phase-Steering Technique with On-Off Power Control for Spectrum Sharing-Based Wireless Sensor Networks
Next Article in Special Issue
A Novel Method for Estimating Monocular Depth Using Cycle GAN and Segmentation
Previous Article in Journal
The Design and Implementation of a Custom Platform for the Experimental Tuning of a Quadcopter Controller
Previous Article in Special Issue
Novel Method of Semantic Segmentation Applicable to Augmented Reality
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Soft Sensor Application in Identification of the Activated Sludge Bulking Considering the Technological and Economical Aspects of Smart Systems Functioning

1
Faculty of Environmental, Geomatic and Energy Engineering, Kielce University of Technology, Tysiąclecia Państwa Polskiego 7, 25-314 Kielce, Poland
2
Faculty of Civil and Environmental Engineering, Gdansk University of Technology, Narutowicza 11/12, 80-233 Gdansk, Poland
3
Faculty of Environmental Engineering, Lublin University of Technology, Nadbystrzycka 40B, 20-618 Lublin, Poland
4
Faculty of Fundamentals of Technology, Lublin University of Technology, Nadbystrzycka 38, 20-618 Lublin, Poland
5
Faculty of Environmental Engineering and Land Surveying, University of Agriculture in Krakow, Mickiewicza 24/28, 30-059 Kraków, Poland
6
Department of Science and Engineering of Materials, Environment and Urban Planning-SIMAU, Polytechnic University of Marche Ancona, 60121 Ancona, Italy
*
Authors to whom correspondence should be addressed.
Sensors 2020, 20(7), 1941; https://doi.org/10.3390/s20071941
Submission received: 29 February 2020 / Revised: 24 March 2020 / Accepted: 25 March 2020 / Published: 30 March 2020
(This article belongs to the Special Issue Smart City and Smart Infrastructure)

Abstract

:
The paper presented the methodology for the construction of a soft sensor used for activated sludge bulking identification. Devising such solutions fits within the current trends and development of a smart system and infrastructure within smart cities. In order to optimize the selection of the data-mining method depending on the data collected within a wastewater treatment plant (WWTP), a number of methods were considered, including: artificial neural networks, support vector machines, random forests, boosted trees, and logistic regression. The analysis conducted sought the combinations of independent variables for which the devised soft sensor is characterized with high accuracy and at a relatively low cost of determination. With the measurement results pertaining to the quantity and quality of wastewater as well as the temperature in the activated sludge chambers, a good fit can be achieved with the boosted trees method. In order to simplify the selection of an optimal method for the identification of activated sludge bulking depending on the model requirements and the data collected within the WWTP, an original system of weight estimation was proposed, enabling a reduction in the number of independent variables in a model—quantity and quality of wastewater, operational parameters, and the cost of conducting measurements.

1. Introduction

The processes occurring in the environment, including urban areas, are very difficult to describe. These changes are dynamic in nature and are governed by a number of external factors, that are random (weather conditions), anthropogenic (traffic volume, amount of water supplied to the network) and local. Taking into account that many of them affect the operational costs of infrastructure, urban development, and the comfort of city residents, it is necessary to predict, control and optimize these factors [1,2]. This approach results in the development of the so-called smart system and infrastructure within smart cities [3,4,5]. Due to the complex process mechanisms and the influence of numerous factors, great quantities of data have to be collected using sensors. This enables us to identify their dynamics and create so-called soft sensors based on the devised the mathematical models [6,7,8]. In numerous cases, this solution enables us to reduce the amount of collected data as well as identify the phenomena, the time course of which is difficult to measure and, for example, generates high operational costs of the measurement system. The approach involving the creation of soft sensors supporting infrastructure is increasingly popular, especially in large cities which require their on-line control and optimization [9], as well as in such technological objects as wastewater treatment plants (WWTPs). The concept and implementation of soft sensors at WWTPs is well established, because they are employed for simulating the complex processes occurring in bioreactors [10]. Soft sensors were used for predicting total nitrogen at the WWTP outlet on the basis of the influent wastewater quality and operational parameters; this enables on-line process control [11]. Another soft sensor concept was also devised for this purpose, where the content of nitrogen compounds was modeled at the outlet based on the primary settling tank operation, internal recirculation, and wastewater temperature parameters [8]. A much more complex soft sensor concept was presented by Canete et al. [12], who predicted the chemical oxygen demand (COD), total suspended solids (TSS), and total nitrogen (TN) values at the WWTP outlet using artificial neural networks. The results of calculations and the settings were verified through the simulations conducted by means of a mechanistic calibrated WWTP model. Luccarni et al. [13] presented the concept of a soft sensor used for N-NH4 prediction in bioreactor chambers based on pH, oxidation-reduction potential (ORP), and dissolved oxygen concentration in the bioreactor (DO), which obtained satisfactory calculation results. Hernández-del-Olm et al. [14], taking into account the high cost of calculating important quality markers at the plant inlet (COD, ammonium nitrogen (N–NH4)) proposed a soft sensor for predicting their value for dry-weather and wet-weather as well as intense rainfall events. The concept of soft sensor for BOD prediction at the WWTP outlet was presented by Yu et al. [15], who employed the learning machine model for its creation, in which the values of the weights obtained were optimized by means of the cuckoo search algorithm. When assessing the processes conducted in treatment plant devices and the quality of treated wastewater, soft sensors are used to analyze multidimensional data sets obtained with gas sensor arrays [16,17,18,19]. At present, attention is also increasingly drawn to the energy aspects at the stage of optimization and control of technological processes, which is reflected in numerous papers on this topic [20,21,22].
Activated sludge sedimentation is one of the essential processes affecting the operation of biological reactors, specifically their efficiency. In WWTPs, it is assessed in terms of the sludge volumetric index (SVI). The parameter enables control and optimization of a bioreactor operation, thus affecting the management of the facility. In the WWTPs with integrated removal of carbon, nitrogen and phosphorus compounds, the recommended maximum SVI is 150 cm3/g [23]. If the value is exceeded, problems may occur with sludge dewatering and with the deteriorated quality of the discharged water. In order to keep the value of SVI within the optimal range, a mathematical model which enables the dynamic control of the SVI by proper selection of settings and a balanced management of the WWTP, must be constructed. On the one hand, the mathematical model must accurately reflect the course of a phenomenon, and on the other, its implementation should not be very complicated. Moreover, the number of the exogenous variables should be as small as possible.
The literature survey [24,25,26] indicates that a predominant majority of the models used for the SVI simulations are those predicting continuous values. These are usually data-mining methods, specifically artificial neural network methods and their complex modifications. Building complex model structures improves predictive power in the presently-known models and reduces the number of independent variables in models [26,27,28], but it also may lead to management problems. In models based on conventional methods, the number of independent variables is usually higher than that in complex models. The literature review [29,30,31,32,33,34,35,36,37,38] indicates that the classification models are built using a number of statistical methods, ranging from analytical (logistic regression) through more modified models employing the concept of regression trees (BT—boosted trees, RF—random forests), to neural networks. In the case of classification models, the basis for the identification of the bulking process involves the limit SVI values indicating problems with sedimentation.
At the stage of model creation, it is necessary to find a balance between the model complexity and its accuracy. The conducted analyses [39,40,41,42,43,44,45,46,47,48] showed that in numerous cases, modification of the RF and BT method leads to an improvement in the prediction capability of the regression tree model. This confirms the essential effect of implementing the boosting method in BT and substituting a single tree with a forest in the RF method. In the case of unsatisfactory simulation results, more complex methods than BT and RF can be employed, i.e., artificial neural networks. The model structure is many times more complex than in the methods mentioned above.
According to [12,49,50,51,52,53] the multi-layer perceptron (MLP) is one of the most frequently used methods. It comprises three layers: the input, the hidden, and the output layers. In this method, the so-called learning step identifies the values of what is called “weight”, which connect the neurons in the consecutive layers. In that method, the initial weight factor values have an essential effect on the course of the learning process. This may affect the predictive power of a model, and the results obtained using RF and BT may even be better than those provided by MLP. The shortcomings of MLP were eliminated in the support vector achine (SVM) method [54,55]. In the case of a non-linear relationship between the model output (y) and explanatory variables (x1, x2, …, xn) transformation of n-dimensional space to K-dimensional (linear) space of variable features is performed using the kernel function. Vapnik [55] has developed a special learning algorithm in the support vector method. Vapnik identified the scale values by solving the problem of square programming, which guarantees the existence of a single minimum. Consequently, at the operating and management levels (intended to obtain the optimum SVI by selecting the correct settings to limit problems with the continuity of processes), using that model for the monitoring and correcting sludge sedimentation is a source of technical problems. Constructing the models with a significant number of independent variables requires a constant monitoring of the values of operating parameters and numerous indicators of the wastewater quality; at the operating level, this creates problems with the process control and management. Obtaining large amounts of data with high resolution under working conditions is difficult for a number of reasons. In order to provide reliable input data for a model, the analyzing equipment must be calibrated on an ongoing basis. Moreover, the hardware sensors are prone to failure or sustaining mechanical damage under real conditions. Although it is still feasible to eliminate such damage or to calibrate the equipment as frequently as necessary, these tasks may be quite expensive and the costs of their implementation under operating and management conditions in a real facility may be disproportionately high. In the facilities operated with no hardware sensors, the indicators of quality and the bioreactor operating parameters are measured with analytical methods. From the technical point of view, obtaining data in a continuous system in that case is rather difficult. Therefore, the implementation of a calculation model in a WWTP may be limited, thus complicating the management of the facility. Despite the many problems with data acquisition, with their reliability and with the costs of measurements, researchers [56,57] still strive to improve the predictive power of models. The approach calls for some criticism, not being fully consistent with the latest trends, as these tend to focus on cost minimization [58,59,60,61]. On the other hand, systems that enable an on-line process control, taking into account the potential failures of measurement schemes, are required. In order to solve this problem, researchers will create new standards for the construction of mathematical models for controlling technological processes in wastewater treatment plants.
Considering the aforementioned observations, the different conditions and data acquisition standards applicable to WWTPs, it seems desirable to seek systems enabling process simulation with high accuracy, taking into account the real-time differences in the operation of WWTP facilities. Therefore, it seems reasonable to aim at such solutions for the simulation of the wastewater treatment processes that enable the processes to be controlled and operation of the facility to be managed regardless of the measured data (including the wastewater quality, operating parameters and their combinations). On the other hand, a process engineer working in a WWTP will be interested in selecting such independent variables that make the model accurate, even though the number of variables included in the model is limited: less measured data will generate lower operating costs but will still enable control of the wastewater treatment process. These assumptions can be implemented under operating conditions and the choice of the simulation method is of major importance. On the one hand, the calculation model must provide a compromise between complexity and accuracy. On the other hand, the question arises whether a model must be created to simulate continuous values. By introducing binary input/output variables, in many cases it is possible to obtain the models in the form of analytical dependencies, which can be applied without the need to implement complex computational algorithms.
As part of the analyses, the concept of a soft sensor system expert system has been proposed which can be applied in selecting a method for the analysis of activated sludge bulking and which takes into account the economic aspects of the wastewater treatment technology. These include the measured indicators of the wastewater quality, the operating parameters of the bioreactor, the duration of the test period, the influence of measurement errors of several independent variables on the results of a simulation, and the complexity of the method used in building the model for soft sensor construction.

2. Materials and Methods

2.1. Experimental Data

Experimental data were obtained from the test facility the WWTP in Sitkówka-Nowiny, located near the city of Kielce in the south of Poland. The daily load of the test facility is Qn = 42,000 m3/d (equivalent to 2,750,000 PE) of municipal wastewater from Kielce and the adjacent area. The wastewater is first handled mechanically on bar screens, in a preliminary sedimentation tank basin and a sand trap. It is then sent downstream for biological treatment in a bioreactor, operated in a modified BARDENPHO system with a pre-denitrification chamber. The treated wastewater is separated from the activated sludge in the secondary sedimentation tank and is then discharged to the receiving body—the Silnica river.
Monitoring in the WWTP includes measurements of the quantity and quality of the influent wastewater and the bioreactor operating parameters. During the research period, a qualitative analysis of both the influent and effluent wastewater was performed once a week to determine its BOD, COD, TSS, TN, TP, N–NH4, and sedimentation parameters of activated sludge (SVI). The organic compounds were determined as COD in accordance with PN–ISO 6060: 2006 and as BOD with the method using the OXITOP, in accordance with PN–EN 1899–1:2002. TSS were determined by means of glass fiber filters with the methods set out in PN–EN 872: 2007. TP was determined in accordance with PN–EN ISO 6878: 2006. TN and nitrogen as N–NH4 were determined in accordance with PN–C–04576–14: 1973 and PN–ISO 5664: 2002, respectively. SVI was determined by means of the method set out in PN–EN 14702–1: 2008. Moreover, the quantity of the influent wastewater and the operating parameters of the reactor were measured, whereas the data were recorded at hourly intervals using SCADA. In this way, 250 datasets representing all the above mentioned variables were obtained. The monitored parameters included MLSS—mixed liquor suspended solids concentration of activated sludge, T—temperature in the bioreactor, DO—oxygen level, mPIX—dosage of PIX chemical coagulant, WAS—excess sludge, RAS—degree of recirculation [return activated sludge], and MLSSR—concentration of return sludge. The results of measurements in the various seasons of the years 2013–2015 are shown in Table 1. These data show that both the bioreactor operating parameters as well as the amount and quality of wastewater varied considerably in the different seasons of the study period, affecting the activated sludge sedimentation to a large extent. In the winter season, the average values of SVI were higher than those from the spring to fall. The relation obtained is confirmed by analysis of variance (ANOVA) test results. For the assumed value of p = 0.05, p-test value equal to p = 0.000024 and F = 42.257 were obtained. The SVI lability depending on temperature may indicate that the activated sludge bulking occurred in the period of interest due to a temperature drop, as confirmed in the analyses reported by Bayo et al. [62]. Other reasons could include a drop in the concentration of organic compounds in the influent wastewater as well as the changes in the oxygen concentration in the activated sludge chambers, as suggested by Comas et al. [29].
These results and the wide range of SVI indicates the requirement to build a mathematical model enabling the construction of soft sensor to control of the process of activated sludge sedimentation.
The literature [23,24,25,26] shows that the sludge bulking is a complex phenomenon. Its course and dynamics are influenced by temperature and pH, because they determine the growth of microorganisms. An increase in temperature raises the number of the floc-forming bacteria. The pH value should range between 7–7.5, since a drop below 6.0 results in the growth of fungi, which leads to sludge bulking. The MLSS and TSS values are the factors describing the amount of activated sludge in a bioreactor. The wastewater quality markers (BOD, COD, TN, TP, N-NH4) are a source of nutrients for microorganisms. Insufficient TN and TP lead to the creation of particles with high floc content and loss of sedimentation ability. The amount of coagulant usually leads to an improvement in the sedimentation of the activated sludge. Oxygen is used during obtaining energy for the biochemical processes conducted by microorganisms in the bioreactor.

2.2. Model Concept

This paper presents a concept of soft sensor for selecting the data-mining method for simulating the activated sludge bulking (or absence of bulking) which takes into account the following aspects: complexity of the data-mining method, time of analysis of independent variables in these simulations, different accessibilities of the recorded measurement data within the test facility, reliability of the measured data, usefulness of the model in process control and optimization (Figure 1). The approach presented below also indicates the possibility of extending the soft sensor system, to select the calculation method for various volumes and qualities of the influent wastewater and for various weather conditions.
In the approach presented, the model for assessing the activated sludge sedimentation was built based on the technological data recorded within the test facility. The various categories of technological data were determined by the type of the measurement system used by the test facility and were limited to the following: quantity and quality of the influent wastewater, the bioreactor operating parameters or their combinations (Step A, Figure 1). The recorded results enabled the construction of a classification model for simulating the activated sludge sedimentation based on the sludge volumetric index, SVI. The binary system was selected, which indicates whether sludge bulking does (0) or does not occur (1).
Sedimentation is a complex phenomenon and the amount of data recorded at the test facility was limited; therefore, the variables which have a statistically significant impact on the investigated phenomenon were identified before the model for soft sensor could be constructed. The Fischer–Snedecor test was used for that purpose (Step B, Figure 1). The variables with numerical values of test probability over p = 0.05 were omitted from the model.
In Step C (Figure 1) the combinations of independent variables (xi) which are considered a basis for the construction of classification models, were created. In these analyses, the combinations of the influent wastewater quality and the bioreactor operating parameters were adopted. This issue is of importance to the feasibility of the control and monitoring of the processes taking place in a complex facility such as a WWTP, considering the numerous technical problems in its operation. Moreover, a lot depends on the function that it is supposed to satisfy. Importantly, the number of independent variables included in the model should not entail high-cost determinations of the wastewater quality indicators or the operating parameters of a reactor. Instead, it should be optimal for the applicable requirements and enable management of the facility.
In the next step of the analyses (Step D, Figure 1), the selected combinations of the independent variables of interest were assigned numerical values in the range 0–3, the sum of which within the respective variables (quality of wastewater and operating parameters) is standardized. The values of the adopted weight factors and the proposed standardizing functions describe the following:
  • the number of independent variables xi relating to the wastewater quality, duration of their determinations, and the bioreactor operating parameters: δ, λ;
  • the cost of measurements of the wastewater quality indicators, bioreactor operating parameters (the value 0 relates to the lowest cost and 3 to the highest cost): f ( δ ) , f ( λ ) ;
  • the duration of measurement of the wastewater quality indicators (0 relates to those measurements in which the duration was lesser than one day, and 1 pertains to those lasting for more than one day): f ( t ) ;
  • the possibility of using the selected model for the management, control or optimization of the WWTP functioning (the weight factor of 1 indicates that the model is applicable in the control, optimization, and the weight factor of 0 means that the model cannot be used for improving the efficiency of the WWTP): F(S).
In the next step (Step E, Figure 1), the combination of independent variables is used for creating classification models using different data mining methods. In the present analyses, before selecting a method the authors considered its complexity, the number of estimated coefficients within its structure, and the possibility of interpreting the results obtained. On the basis of results of such calculations, it is determined how they fit to those of the measurements—SPEC (determines the correctness of data classification in the set of data, including events when a activated sludge bulking occurred), SENS (determines the correctness of data classification in the set of data constituting cases when no activated sludge bulking occurred) [63,64,65].Using these results, the model class is found and, depending on the weight factor value, the calculation method is selected.
The established standardized weight values expressing the number of the independent variables that describe such parameters as the wastewater quality δ = f(z1, z2, …, zq), bioreactor operating parameters λ = f(f1, f2, …, fg), bacterial flora (e1, e2, …, et), cost of analyses f(δ) and f(λ), duration of analyses f(t), and usefulness of the model for control and management of the bioreactor operation F(S), were used in determining the multi-dimensional weight factor vectors [δ f(δ) λ f(λ) f(t) F(S)]. They enabled the matrices to be found facilitating a fast and simple choice of the suitable data-mining method for identification of the activated sludge bulking (Step F, Figure 1).
When selecting a model, the proposed solution offers the possibility to introduce economic criteria to vary the independent variables that are included in the model (Step G, Figure 1). These independent variables relate to the number of data to be measured (volume and quality of the wastewater, and the reactor operating parameters) as well as to the possibility of using a specific model for the control, management and optimization of the WWTP operation.
The resulting soft sensor system for selecting the optimal data-mining method for the analysis of the activated sludge sedimentation can be applied in the daily operation of the WWTP and in its management, based on a variety of calculation methods, taking the wastewater volume and the bioreactor operating parameters into account.

2.3. Methods of Data Mining, Choice of Independent Variables and ConstructionCriteria for Soft Sensor

Taking into the account the issues with optimal selection of a method for the simulation of wastewater treatment plant processes, which was discussed in the introduction, the classification models were adopted for the analysis of the phenomenon considered. In order to determine the optimum between the accuracy and complexity of a model, as well as the number of independent variables, an analytical model—logical regression (LR)—was considered in the analyses. The model had the following form:
p = e x p ( j = 1 k β j · x j + β 0 ) 1 + e x p ( j = 1 k β j · x j + β 0 )
in which the established βj coefficients enable determining the impact of the input data (xj) on the probability of the occurrence of the phenomenon (p)—activated sludge bulking in this case.
Simultaneously, the application of models with increasing complexity was tested, starting from relatively simple ones (RF, BT) and finishing with a neural network methods like MLP and its subsequent modification, i.e., SVM. The mentioned methods, which were applied in this research work in creating the models for soft sensor predicting the activated sludge bulking, are highly efficient, as confirmed by a number of research works [32,49,62]. Taking into account the above, these models have been implemented in such statistical packages as R, STATISTICA, XLSTAT, SPSS, and can be used for building the models for predicting activated sludge bulking by many groups of users, including those responsible for the management of a WWTP.
In this elaboration, the statistical models were indicated by means of STATISTICA 12 software package. In order to devise models, the measurement data were divided into two sets: a learning set (75%) and a test set (25%). In the case of the logit model, based on the appropriate combinations of independent variables and on the assumed confidence interval of 0.05, the empirical coefficients were determined using formula (1). In the RF and BT model, assuming the number of trees up to 300, the structure of models was determined. In the MLP models, the number of neurons in the hidden layer was analyzed in the range (j–2·j + 1) (where: j—number of independent variables). In order to optimize the structure and weight values for the assumed number of neurons, different activation functions were analyzed, including: linear, exponential, sigmoidal, sinus, tangent—hyperbolic. For a complex number of neurons in the hidden layer and consecutively assumed activation functions, the fitting of calculation results to the measurements was determined. When the calculated values of SENS and SPEC for the assumed independent variables and the number of neurons as well as the activation function were minimal, the structure obtained was considered optimal. In the SVM method, Gaussian kernel was assumed, while the optimum values of C and γ were sought with the iteration method, substituting the values of above-mentioned parameters until the minimal SENS and SPEC values were achieved.
Given the fact that sludge bulking depends on the wastewater volume, its quality, operating parameters and the technology, the following general relationship can be formulated [11,62,65,66,67]:
{ S V I = f ( Q ,   z 1 ,   z 2 ,   z 3 ,   z q ,   f 1 , f 2 , f 3 , f g ,   e 1 ,   e 2 , e 3 , , e t ) S V I S V I l i m   a c t i v a t e d   s l u d g e   d o e s   n o t   b u l k      S V I > S V I l i m   a c t i v a t e d   s l u d g e   b u l k s  
where: z1,2,3,q—independent variables describing the quality of wastewater, f1,2,3,g—independent variables describing the bioreactor operating parameters, e1,2,3,t—independent variables describing the bacterial communities present in the activated sludge, and the technological solution of the bioreactor, q—the number of wastewater quality indicators, included in the mathematical model, g—the number of operating parameters of the reactor, included in the model, t—the number of independent variables, which describes the bacterial communities and has been included in the model.
Taking into account the complexity of the sedimentation process and the differences in the accessibility of the measured data within a WWTP, this paper provides detailed analyses to find the specific combination of data including the values of z1,2,3,q and f1,2,3,g for which the sludge bulking process can be identified. In the analyses, the variables e1,2,3,t (Formula 2) were omitted. The influence of micaceous organisms, including filamentous bacteria, on the results of calculations of volumetric index of activated sludge was analyzed by Bezak-Mazur et al. [68]. Consequently, the following models were developed for the analyses:
S V I q , g = f (   Q , z 1 ,   z 2 ,   z 3 ,   z q ,   f 1 , f 2 , f 3 , f g )
S V I q = f (   Q , z 1 ,   z 2 ,   z 3 ,   z q )
  S V I g = f (   f 1 , f 2 , f 3 , f g )
The model described by the relationship (4) includes the volume and indicators of the quality of the influent wastewater. Within the respective models, described by Equations (3)–(5), the authors sought such combinations of independent variables, including for instance the wastewater quality, of which the number is as small as possible and the results of simulation are consistent with the measured data:
{ 1 , 2 , 3 , q ,   1 , 2 , 3 , g } min   a n d   S P E C ( S E N S ) m a x
{ 1 , 2 , 3 , q } min   a n d   S P E C ( S E N S ) m a x S V I q = f (   Q , z 1 ,   z 2 ,   z 3 ,   z q )
{ 1 , 2 , 3 , g } min   a n d   S P E C ( S E N S ) m a x
On the basis of the work [69] it was shown that for the Siktówka-Nowiny WWTP there is a general relationship:
{ S V I ( t ) = f (   Q ( t 1 ) ,   B O D ( t 1 ) , T N ( t 1 ) , T P ( t 1 ) ,   M L S S ( t 1 ) , D O ( t 1 ) ,   m P I X ( t 1 ) , T ( t 1 ) )       S V I 150   a c t i v a t e d   s l u d g e   d o e s   n o t   b u l k          S V I > 150   a c t i v a t e d   s l u d g e   b u l k s
where: (t)—values of independent variables measured at moment t; (t–1)—values of independent variables measured at moment (t–1).
On the basis of the relationship (9), it can be stated that the load of organics, nitrogen and phosphorus in the influent wastewater has a statistically significant effect on the sludge sedimentation in the test facility. The operating parameters of the reactor, such as the concentration of activated sludge, oxygen level, dosage of PIX, as well as temperature in the sludge chambers strongly correlate with the seasons of the year [33,62,69]–are important as well. After substituting the obtained relationships in Equations (3)–(5), the following was obtained:
S V I ( t ) = f ( Q ( t 1 ) ,   B O D ( t 1 ) , T N ( t 1 ) , T P ( t 1 ) ,   M L S S ( t 1 ) , D O ( t 1 ) ,   m P I X ( t 1 ) , T ( t 1 ) )
S V I ( t ) = f ( Q ( t 1 ) ,   B O D ( t 1 ) , T N ( t 1 ) , T P ( t 1 ) , T ( t 1 ) )
S V I ( t ) = f (   M L S S ( t 1 ) , D O ( t 1 ) ,   m P I X ( t 1 ) , T ( t 1 ) )
In the case of interest, for the Sitkówka-Nowiny WWTP, operated as an integrated system for the removal of carbon, nitrogen and phosphorus compounds, it was assumed that SVIlim = 150 cm3/g, as reported in the literature [23].

2.4. Determining the Values of Weight Factors and Matrices of a Method Selection for Identification of Sludge Bulking

For a quantitative assessment of the differentiation of the adopted combinations of independent variables in predicting activated sludge sedimentation, as described above, the authors proposed an original system for determination of weight factors. The total effect of the aforementioned factors is expressed using the cumulative weight factor (Wtot) the minimum value of which, taking into account the model-building criteria (mentioned below), can be calculated from the formula:
W t o t , m i n = m i n { m i n { δ , f ( δ ) , f ( t ) } + m i n { λ , f ( λ ) } + F ( S ) }
The general relationship (14) takes the following form:
W t o t = { δ , f ( δ ) , f ( t ) } + { λ , f ( λ ) } + F ( S )
where: δ, f(δ), f(t), λ, F(λ), F(S)–according to the notation in the subsequent sections.
In practice, numerous variants of independent variables in models are encountered in the operation of a WWTP. This depends on the data obtained in measurements in the facility and on the requirements that a given model is expected to meet.
In Equation (13), the boundary conditions (π) facilitating the choice of independent variables for predicting SVI values were defined in parallel and in series:
  • the condition determining the costs of determination and the number of indicators of the wastewater quality which are included in the model and reducibility of their value:
    π 1 : m i n { δ , f ( δ ) , f ( t ) }
  • the condition determining the cost and number of measurements of the operating parameters which are included in the model and their reducibility:
    π 2 : m i n { λ , f ( λ ) }
  • the condition enabling minimization of the number and costs of measurements pertaining to the values of wastewater quality indicators and the reactor operating parameters:
    π 3 : m i n { δ , f ( δ ) , f ( t ) } + m i n { λ , f ( λ ) }
  • the condition minimizing the number and cost of measurements of the wastewater quality indicators and the reactor operating parameters so that the obtained model can be used in controlling the reactor operation:
    π 4 : m i n { m i n { δ , f ( δ ) , f ( t ) } + m i n { λ , f ( λ ) } + F ( S ) }
In spite of the defined boundary conditions, the choice of the suitable data-mining method for the identification of activated sludge bulking and for the management of the WWTP operation is possible if the operating parameters of the bioreactor are included; this corresponds to the following relationships (16), (17), (18).
The data-mining method for simulating the volumetric index of activated sludge is selected after establishing independent variables in the model, according to the aforementioned calculation algorithm.
The function δ in Formula (14) describes the number of the quality indicators which were adopted in simulating the values of SVI. It takes the following form:
δ = i = 1 q [ z i ] i = 1 q [ z i ] m a x
where: zi = 1, 2, …, q–number of the wastewater quality indicators included in the model; if the indicator is included in the model, then [zi] = 1, otherwise [zi] = 0, [zi]max—maximum weight of 1 indicating that the i-th indicator is included in the calculation model, in the case below q = Σ[zi]max = 3.
The symbol f(δ) is the function of costs showing the impact of the cost of determining the quality indicators; it is described by the equation:
f ( δ ) = i = 1 q [ K ( z i ) ] i = 1 q [ K ( z i ) ] m a x
where: K(zi)—weight of cost of measurement of the i-th indicator of the wastewater quality, included in the model; in this case, the number of indicators q = 3, and they are as follows: BOD, TN and TP. A literature survey [70,71] showed that the following relationship is true: K(TN) > K(TP) > K(BOD), where the value of TN is found by calculation Kim et al. [72]; therefore, the weight factors of K(TN) = 3 and K(TP) = 2 and K(BOD) = 1 were adopted in the calculations. For instance, if BOD and TP are included in the model, then the value of the function f(δ) = 1 + 2 + 0/(1 + 2 + 3) = 3/6. The proposed method takes into account the fact that the adopted weight factor values may be uncertain. In that case, it is acceptable to assume that the values K(zi) are stochastic and may be predicted with the Monte Carlo method.
The function f(t) in Formula (15) takes into account the duration of determination of the i-th indicator of the wastewater quality for simulating the values of SVI. It takes the values f(t,zi) = 1 for the time of determination longer than 24 h, otherwise, it is 0. In the model of interest, it was assumed that the values of the wastewater quality indicators in Equations (10), (11) may be modeled using statistical models [69,73,74,75]—for these cases, f(t) = 0.
The function λ describes the number of the bioreactor operating parameters adopted in calculating SVI. The function takes the following form:
λ = k = 1 g [ z k ] k = 1 g [ z k ] m a x
where: zk = 1, 2, …, g–number of the bioreactor operating parameters included in the model, if a parameter is included in the model, then [zk] = 1, otherwise, [zk] = 0, [zk]max—the maximum weight factor value of 1 indicates that the k-th operating parameter is included in the calculation model; in the case below, it is q = Σ[zk]max = 3 (namely, DO, MLSS, T).
The symbol f(λ) is the function of costs, describing the impact of the cost of determining the quality indicators. It can be described by the following equation:
f ( λ ) = k = 1 g [ K ( z k ) ] k = 1 g [ K ( z k ) ] m a x
where: K(zk)—the weight factor of the measurement of the k-th bioreactor operating parameter, included in the model; in this case, the number of the indicators included is g = 3, and these are: DO, MLSS and T. Literature [61,70] has shown that the following relationship is true: K(MLSS) > K(DO) > K(T). Therefore, the following weight factors were adopted in the calculations: K(MLSS) = 3, and K(DO) = 2, and K(T) = 1. For instance, when the values MLSS and T are included in the model, then the value of the function is f(λ) = 1 + 3 + 0/(1 + 2 + 3) = 4/6.
Under operating conditions, the measured values from several hardware sensors are exposed to error, potentially affecting all the devices. Therefore, in sensitivity analyses, the recommended approach is to analyze the impact of several factors on the results of calculation. This can be described by the following equation:
S ( x 1 , x 2 , , x j ) = Δ p p 0 ( x 1 , x 2 , , x j ) = p ( x 1 + p · Δ x 1 , x 2 + p 1 · Δ x 2 , , x j + p j · Δ x j ) p 0 ( x 1 , x 2 , , x j )
Therefore, the function F(S) takes into account the fact that the selected independent variables in a model enable the control and adjustment of the SVI values. This can be described as follows:
F ( S ) = f ( λ , δ )
If the variables adopted in the calculations enable the control of the SVI values, then the function F(S) = 0, otherwise, F(S) = 1.
If the essential criteria for the model creation, encompassing the measured values showing the wastewater volume and quality and the bioreactor operating parameters are adopted, then a matrix of cases can be written (Figure 2) and the weight factor values W can be found from Equation (13). For the established combinations of independent variables in the matrix above and for the weight factor values Wtot, the mathematical models for the identification of the activated sludge bulking are found using the selected data mining methods.
The obtained results of simulation are expressed in the form of a matrix for the selection of methods (Figure 3). Such a matrix includes the models with the optimum fitting between the results of calculation and the measured results among all the methods considered.
In order to supplement the aforementioned methods and the matrices, additional criteria can be included which govern the developed model (controllability—MΘ,F(S), duration of determination—MΘ,f(t)). A representative matrix for the selection of a method for the simulation of sludge bulking is shown in Figure 3. For MΘ = 0, the selected method cannot be applied for assessing the capacity of sedimentation because independent variables were included in the model or due to the duration of determination of the wastewater quality indicators, which are included in the model. For instance, for MΘ = MΘ,F(S) = MΘ,f(t) ≠ 0, the selected data-mining method can be used for predicting activated sludge bulking.

3. Results and Discussion

3.1. The Sludge Bulking Identification Method and the Choice of Independent Variables

On the basis of the results of measurements in the Sitkówka-Nowiny WWTP, the classification models for predicting SVI using the boosted tree (BT), random forest (RF), multilayer perceptron (MLP), support vector machines (SVM) and logistic regression (LR) methods were determined using the relationships (3)–(8), (11), (12), and (13). The variables A–G, which were statistically significant to the phenomenon of interest, were used in various combinations. Table 2 presents the values of measures of fit between the calculated and the measured results (SENS and SPEC).
The data in Table 2 indicate that the models found by means of the BT and RF methods are not over-learned. This is confirmed by the number of trees (N.T.) for the consecutive calculation variants (<300 trees). In the MLP method, the number of neurons is not greater than 2 · j + 1 (where: j –number of independent variables, included in the model). This seems to indicate that the obtained models are not over-learned [24]. In the hidden and output layers, usually the following activation functions are present: linear, exponential, tangent—hyperbolic, sinusoidal, and logistic.
The analysis of the data in Table 2 also indicates that the classification models for predicting SVI are characterized by different predictive powers. This is largely due to the choice of the independent variables which include the wastewater quality and the bioreactor operating parameters (or their combinations) as well as the selected methods (BT, RF, MLP, SVM, RL). The results of simulation have confirmed a significant effect of the content of biogenic compounds on the activated sludge sedimentation in the test facility. The parameter indicates that the phenomenon of activated sludge sedimentation is taking place [29].
The predictive power was observed to improve in the classification models after the inclusion of the various independent variables (BOD, TN, TP). Such a relationship can be perceived based on the model formed with the RF, BT and SVM methods. However, the approach where the wastewater quality indicators (A, B and C) and temperature in the activated sludge chambers were included, enables the management of the facility and application of the model for a soft sensor in optimizing the operation of the WWTP. These results were also confirmed in the studies of Pretorius and Rössle [76], who formulated the empirical relationships for predicting SVI from MLSS, T obtaining the value of R2 = 0.84. The above-mentioned relationship was also confirmed by the analyses of Szeląg and coauthors [69], who indicated the instrumental influence of the independent variables on sludge bulking by employing the logistic regression method. The possibility of employing logistic regression for sludge bulking identification was also shown by Bayo et al. [62], who presented the seasonal nature of the phenomenon, which confirms the influence of temperature on its course. An analysis of the examples of calculations for various combinations of independent variables A–G, indicates that the determination of a model with a high classification power does not require the inclusion of as many factors as the regression models. Such complex calculation algorithms are not necessary, either. Bagheri et al. [77] obtained a model of soft sensor characterized by a good fit with the measured data after including COD, TN, T, MLSS, DO, pH, TSS and modifying the MLP method in which the weight factors were corrected using a genetic algorithm. A similar approach to that of Bagheri et al. [77] was described by other researchers [26,56,57,78], focusing on the modification of an artificial neural network model in order to obtain better predictive power of their model. Interesting results of analyses were presented by Boztoprak et al. [25], who obtained reliable calculation results through observations with high-resolution cameras and the application of cellular neural networks (weights were optimized with genetic algorithms). Han and Qiao [78], based on the measurement results of only COD, TN, BOD, pH, DO and using modified MLP neural network models (hierarchical, hybrid– radial basis function (RBF) networks with clustering model (Kohonen maps) etc.), indicated the possibility of bulking simulation with a limited number of independent variables in comparison to the works of other authors [31,69,79].
Nonetheless, the use of such complex model structures may be problematic during the implementation in a WWTP facility. A comparison of the simulation results described by the aforementioned authors with the results of the studies indicates that there is no need to modify the methods used in the research work. The results of simulation, obtained with those methods for data relating only to the wastewater quality or only to the bioreactor operating parameters enable identification of the sludge bulking with high accuracy. This is confirmed by the obtained values of SPEC, SENS (Table 2). Looking at the obtained results in the aspect of reducibility of the number of independent variables, depending on the adopted method [24,69], one can conclude that using the classification model for the identification of sludge bulking and a rather limited number of independent variables, it is possible to obtain a good fit between the calculated and the measured results. This effect can be seen using non-modified methods (RF, BT, SVM, MLP, LR) and a considerable number of input data, including both the volume and quality of the wastewater and the bioreactor operating parameters. The aforementioned aspects were discussed in the work of Han and co-authors [26], and Han and Qiao [78], who compared the results of SVI value simulations obtained using modified models with those found by means of classic models (ARX, SVM, MLP), where the model structure was not corrected. It is also worth noting that the resulting soft sensor models can be used for the on-line control and adjustment of the activated sludge bulking. This allows improving the operation of the test facility and enables its real-time management. On the other hand, a compromise between the number of independent variables and accuracy of the soft sensor model is desirable. A large number of independent variables, as presented in the paper by Luo and Zhao [24]—Q, BOD, COD, TN, N-NH4, TP, TSS, T, MLSS, pH—or a higher number of measurable indicators of the wastewater quality, as described by Bagheri et al. [77], may lead to the problems with the continuity of independent variables and applicability of the soft sensor model in the management of a WWTP.

3.2. The Choice of a Sludge Bulking Identification Method for Various Independent Variables of Model

The data in Table 2 indicate that among the data-mining methods of interest in this paper, a high measure of fit (SPEC > 0.8 and SENS > 0.7) between the calculated and the measured results for SVI could be obtained based solely on the temperature in the activated sludge chambers, using a model built according to the RF and BT methods in the simulation. The results of calculations for the activated sludge simulation are improved by 10% after taking into account the values of wastewater quality indicators such as BOD, TN and TP in the model. The analysis results obtained are confirmed by the findings of Bagheri et al. [77] who indicated the influence of wastewater quality on the improvement of the prediction capabilities of the model (MLP) for the sludge bulking simulation. Similar dependences were also shown by Mirbagheri et al. [79], who investigated the influence of independent variables (wastewater quality, operational parameters) on the results of COD and TP simulations at the wastewater treatment plant discharge.
As an alternative to BT, the SVM algorithm can also be used for sludge bulking identification, except that this is one of the most complex models of all those considered in this research work. It is also worth noting that among all the employed methods, the BT method enables the prediction of SVI from indicators of the amount and quality of the influent wastewater. This is possible knowing just the value of total nitrogen. Inclusion of the other independent variables, such as TN, TP, one by one enables conducting the simulation calculations for the activated sludge bulking with high accuracy (SENS, SPEC > 0.80). Even though the sedimentation capacity of the sludge can be identified based on the measured temperature, amount and concentration of contaminants in the influent wastewater, these values cannot be used for adjusting the bioreactor operating parameters. In the aspect of operation and reliability of a WWTP, these models are applicable when the influent conditions do not indicate any possibility of activated sludge bulking. Similar dependences were obtained by Bayo et al. [62], who indicated the possibility of the sludge bulking process identification based on the knowledge of weather conditions, omitting the operational parameters of the bioreactor. Such a case and an analysis of the condition of activated sludge in the aspect of its bulking are only possible for a linear combination of independent variables, including the amount and quality of the wastewater, and weather conditions [62]. However, for the aforementioned data, the possibility of optimal adjustment of the bioreactor settings is limited, and this affects the operation of the facility.
From the simulation calculations performed, it can be observed that among the bioreactor operating parameters of interest (MLSS, DO), which enable adjustment of the SVI values, the lowest predicting errors are obtained for the model determined from MLSS. The aforementioned dependence is confirmed by the analyses conducted by Comas et al. [29], who devised an expert system for the assessment of the impact and investigating the interaction between the wastewater quality at the inlet and outlet, operational parameters and activated sludge bulking. In that case, a satisfactory fit (SPEC = 0.96 and SENS = 0.87) between the calculated and the measured values is obtained by means of the BT method. Lower values of error were obtained by MLP (SENS by 4% and SPEC by 22.67%) as well as SVM (SENS by 3.12% and SPEC by 17.95%). Reduced measurement errors (SENS, SPEC values) were obtained in the models, taking into account the sludge temperature. Therefore, for the variables MLSS, T, satisfactory results of calculations were obtained by means of BT, RF, SVM, MLP and LR. Under the operating conditions, the soft sensor model based on logistic regression is the simplest and easiest to implement. In the case of landslide identification, a good fit between the results of calculations obtained with a logit model was achieved, in relation to more complex statistical methods (SVM) [80]. Similar results were obtained through the analysis of methodical data, while comparing the results of the analyses conducted with the ANN method [81]. The calculation of SVI from DO alone provided satisfactory results with BT, RF, as confirmed by the values of SPEC >0.9 and SENS of at least 0.7. In the case of a combination of independent variables including DO and T, a good fit between the calculated and the measured data was obtained with BT, RF as well as by means of MLP and SVM, which are more complicated methods. A compilation of the aforementioned independent variables (which are the controlling variables since they enable adjustment of SVI, namely MLSS, DO and T), enables the activated sludge bulking identification using the methods based on neural networks (MLP, SVM), developed from the regression tree concept, that is BT and RF, as well as using them for monitoring and for the control of technological processes within a WWTP. The results obtained show that inclusion of the values of wastewater quality indicators (BOD, TN,TP) in the calculation model SVI = f(DO, T) or SVI = f(DO) provides lower errors since SENS, SPEC are higher. Among the methods considered, the lowest values of errors were obtained for the RF and BT methods in most cases. Using the SVM method, good measures of fit between the measured and the simulated data were obtained for the variables including the values of DO, T, indicators of amount and quality—BOD, TN. Among the classification models discussed in this paper (Table 2), the best predictive power was shown by those in which the values of wastewater quality indicators as well as the bioreactor operating parameters were included. For that combination of independent variables, the results of simulation are characterized by a good fit with the measured values for those methods in which the soft sensor model structure is not very complex, the number of calibrated parameters is small (RF, BT, LR) as well as in more complex methods (SVM, MLP). Significant similarity of the results pertaining to the calculations of technological wastewater treatment plant parameters (energy consumption, biogas production) obtained with the RF and MLP methods was confirmed by the analyses performed by Kusiak et al. [58,82]. Among the combinations considered in this paper, the best results of calculations were obtained for the combination involving Q, TN, TP, MLSS, DO, T, mPIX. On the one hand, these data contribute the information on the number of compounds that the microorganisms in the bioreactor are provided with and, on the other, they describe the dynamics of the processes taking place in the activated sludge. In this approach, the variables are supplemented with the coagulant dosage (it enables the adjustment of SVI and of bulking so that the values of MLSS and DO can be kept within the optimal range). The proposed approach offers the possibility of choosing the optimal settings in the operation of a WWTP and improving the sedimentation capacity of the activated sludge, even though the values of MLSS, DO may not guarantee SVI below 150 cm3/g.
Moreover, the data in Table 2 show that the activated sludge bulking can also be identified for a slightly reduced range of independent variables as well as for other combinations, and the results of calculations will also fit well to the measured values (BOD/TP, BOD/TN, DO, MLSS, T, mPIX). For instance, the logit model for a combination of independent variables relating to the amount and quality of wastewater and the bioreactor operating parameters is described by the following relationship:
X = 28.24 0.68 · T 2.66 · M L S S 1.82 · D O 0.56 · m P I X + 0.0008 · Q · T N + 0.0009 · Q · T P + 0.0001 · Q · T P
where: LTN—total nitrogen in the incoming wastewater, LTP—total phosphorus in the influent wastewater, LBOD—biochemical oxygen demand of the influent wastewater. The model above comprises the following variables having a random effect on the sludge bulking:
ω = 0.0008 · L T N + 0.0009 · L T P + 0.0001 · L B O D
and the controlling variables:
τ = 28.24 0.68 · T 2.66 · M L S S 1.82 · D O 0.56 · m P I X
where: τ, ω, LTP, LTN, LBOD, MLSS, DO, T, mPIX–independent variables, marked as above.

3.3. Impact of Uncertainty of Measured Data on Activated Sludge Bulking Identification

In order to demonstrate that it is necessary to continually calibrate the hardware sensors installed in the WWTP for measuring T, MLSS, DO and TN, the values of sensitivity coefficients (Sxj, %) were determined. Those calculations included the measurement errors for a single hardware sensor and for two and three hardware sensors. The measured values of the sensitivity coefficients, based on the relationship (23) are shown in Figure 4, Figure 5 and Figure 6. Moreover, the analyses included the determination of variability of the sensitivity coefficient vs. measurement error for the values MLSS, DO, TN; (Figure 7). The measurement errors for MLSS, DO were expressed as the coefficient τ, described by Equation (27).
The plotted curves in the figures below indicate that the accuracy of the measurement of MLSS, DO, TN has a significant effect on the values of the sensitivity coefficients—SMLSS,DO (Figure 4) as well as ST,DO (Figure 5), SMLSS,TN (Figure 6) and Sτ,TN (Figure 7). The curves in Figure 4, Figure 5, Figure 6 and Figure 7 indicate that a significant impact on the uncertainty of identification and on the monitoring of the sedimentation capacity of the activated sludge is particularly attributable to the errors resulting in the over-estimation of real values (theoretical ones, not affected by the measurement errors). In particular, this is visible if two or three hardware sensors fail at a time. It is also worth noting that the analyses described by other authors in previous reports only related to the single hardware sensors. The curves obtained (Figure 4, Figure 5 and Figure 6) show that the proposed approach enables the analyses of a bioreactor operation in a wider range than before, should the measured data be uncertain. Thus the results obtained are confirmed in the analyses performed by Comas et al. [29], who indicated a diversified influence of bioreactor operational parameters (MLSS, DO) and the quality of wastewater at the inlet (BOD, TN) on the activated sludge bulking by employing the fuzzy set method. The aforementioned relation is also confirmed by the analyses of Bagheria et al. [77], who identified the influence of the analyzed independent variables describing the wastewater quality (COD, TN, TSS) and operational parameters (T, MLSS, DO, pH) on the activated sludge bulking.
This is the condition of the optimal choice of settings for the bioreactor, to guarantee the desirable sedimentation capacities (SVI <150 cm3/g) and high reliability of the operation within the entire WWTP. These two factors are directly connected with the correct management of the WWTP operation, resulting in unhindered sludge dewatering and providing the optimal values of quality indicators of the effluent wastewater.

3.4. Matrices of the Choice of a Sludge Bulking Identification Method

From the weight factor values, established for the wastewater quality indicators (BOD—1, TN—3, TP—2) and for the bioreactor operating parameters (MLSS—3, DO—2, T—1), respectively, the numerical values for the functions f ( λ ) and f ( δ ) were determined. On this basis, the vectors [ f ( λ )   f ( δ ) ] were determined; then, it was established which of the considered methods should be used to make the soft sensor model useful for the purposes of calculations and predicting the sludge bulking, without taking into account any aspects concerning the process control and management (Table 3 and Table 4). The models for which the established values of the measures of fit for SPEC, SENS were the highest among all of those considered in this paper were regarded as the dedicated soft sensor models. Table 3 shows those methods for which the obtained values of the measures of fit were the maximum. Table 4 shows the applicability of the logistic regression LR (as the only one among the methods considered in which the chosen model is in the form of an empirical equation) for the sludge bulking simulation for the value [ f ( λ )   f ( δ ) ] . The same table comprises the cases for which the established values of SENS, SPEC >0.8 (a very good model—VG) and SENS, SPEC >0.9 (an excellent model). In the other cases of calculations, where one of the above-mentioned conditions was not met, the sludge bulking identification methods, referred to in Table 3 were adopted.
On the basis of the matrices, determined in Table 3, it is possible to choose the suitable method without having to perform the simulation calculations with various methods, taking into consideration the measured data that are available from the WWTP and that include the following: quality indicators (δ), their costs expressed as the function f(δ), the bioreactor operating parameters (λ) and their costs—f(λ). Moreover, to simplify the method selection with regard to the number of independent variables in the system: wastewater quality indicator–reactor operating parameters (1–1, 2–2, 3–3), the typical regions for such combinations were determined, which also helps choosing the bulking simulation method. Other variants (1–2, 2–1, 2–3, etc.) depend on the local conditions and data recorded within the specific WWTP. The proposed approach enables diagnosing and supervising the activated sludge bulking in the continuous system, which directly affects the input data (independent variables) included in the model. Moreover, the results in Table 3 and Table 4 indicate that—in many cases—logistic regression is an alternative to complex calculation models. The model obtained is the empirical equation in which the established coefficients can be used for determining the effect of selected independent variables on the activated sludge sedimentation process.
The other methods (BT, RF, MLP, SVM) require performing additional calculations. The matrices obtained and described above (Table 3, Table 4) enable the methods to be matched to simulation with respect to the number of the wastewater quality indicators, the bioreactor operating parameters, and the costs of their determination. Nonetheless, taking into consideration the operation of the wastewater treatment facilities, the soft sensor models which can be applied and which are potentially useful in the monitoring and supervision of the wastewater treatment processes are of the greatest importance. Therefore, using the developed algorithm and introducing a limit of the duration of determination indicators—f(t), a matrix of cases was obtained and is shown in Table 5.
The resulting matrix is an extended version of the matrix shown in Table 3; however, the variants included in it can be used for monitoring the operation of a WWTP, for its control, as well as for the adjustment of the bioreactor settings. The number of feasible variants results from the fact that the BOD determination is a time-consuming process, which affects the applicability of the model. However, assuming that the value of the parameter or that of the other selected wastewater quality indicators can be modeled by means of a statistical model (value of the function f(t) = 1), every case needs to be analyzed individually. The cost estimation is made even more complex by the fact that the costs of analysis pertaining to the wastewater quality indicators need to be assessed on a case-by-case basis. The reason is that it is very difficult to establish independent variables describing the values of selected wastewater quality indicators taking into account the local conditions.

4. Conclusions

The results of the analyses presented above show that the soft sensor model developed based on the selection of a data-mining method for sludge bulking simulation takes into account both the number of the wastewater quality indicators and the bioreactor operating parameters, as well as the costs of their determination. The reliability of the input data for the soft sensor model is an important criterion, as shown by a detailed analysis of sensitivity which enabled an assessment of the effect of one, two or three independent variable(s) on the results of calculation of the activated sludge bulking.
The simulations indicate that the procedure of soft sensor development described in this paper helps in choosing the calculation methods (RF, BT, SVM, MLP, LR) for the sludge bulking simulation, while minimizing the costs of measurements of such parameters as: wastewater quality indicators, the bioreactor operating parameters, duration of determinations, and the number of the independent variables included in the soft sensor model. The calculations show that a suitably selected simulation method enables the activated sludge bulking to be predicted with high accuracy from the load of contaminants arriving in the WWTP with the influent wastewater or the bioreactor operating parameters, and the combinations thereof.
The simulations and the matrices of the cases obtained and described above seem to indicate that the management of a WWTP in the aspects of its control, adjustment of settings and monitoring of activated sludge bulking can be performed by means of complex calculation methods such as neural networks and their modifications, having a simpler structure of mathematical algorithms such as logistic regression. The proposed soft sensor model enables on-line diagnosis in the operation of a WWTP because it includes various conditions potentially encountered in the course of a WWTP operation (failure of the wastewater quality hardware sensors, lack of continuity of measurements, errors in measuring the operating parameters, etc.).
The selection of the data-mining method depends on the user, although it is governed by the availability, the commercial attractiveness of the available software for modeling the operation of WWTP, and how the given software is operated. The soft sensor models for predicting the activated sludge bulking by means of the methods proposed in this paper (RF, BT, SVM, MLP, LR) can be built using the generally available statistical models, which are accessible to many groups of users.

Author Contributions

B.S. developed the concept of the article and wrote the outline of the paper draft; J.D. G.Ł, D.M., E.D., and F.F. analyzed information, and improved the paper draft. All authors of this article provided substantive comments. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by National Science Centre as a result of the research project no. 2017/26/D/ST8/00967 and Ministry of Science and Higher Education in Poland, within the statutory research of particular scientific units under subvention for a science program.

Acknowledgments

The participation of the Sitkówka-Nowiny WWTP staff is gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Castillo, A.; Porro, J.; Garrido-Baserba, M.; Rosso, D.; Renzi, D.; Fatone, F.; Gómez, V.; Comas, J.; Poch, M. Validation of a decision support tool for wastewater treatment selection. J. Environ. Manag. 2016, 184, 409–418. [Google Scholar] [CrossRef] [PubMed]
  2. Longo, S.; Mauricio-Iglesias, M.; Soares, A.; Campo, P.; Fatone, F.; Eusebi, A.L.; Akkersdijk, E.; Stefani, L.; Hospido, A. ENERWATER—A standard method for assessing and improving the energy efficiency of wastewater treatment plants. Appl. Energy 2019, 242, 897–910. [Google Scholar] [CrossRef] [Green Version]
  3. Hollands, R.G. Will the real smart city please stand up? City 2008, 12, 303–320. [Google Scholar] [CrossRef]
  4. Hashem, I.A.T.; Chang, V.; Anuar, N.B.; Adewole, K.; Yaqoob, I.; Gani, A.; Ahmed, E.; Chiroma, H. The role of big data in smart city. Int. J. Inf. Manag. 2016, 36, 748–758. [Google Scholar] [CrossRef] [Green Version]
  5. Visvizi, A.; Lytras, M. It’s Not a Fad: Smart Cities and Smart Villages Research in European and Global Contexts. Sustainability 2018, 10, 2727. [Google Scholar] [CrossRef] [Green Version]
  6. Chéruy, A. Software sensors in bioprocess engineering. J. Biotechnol. 1997, 52, 193–199. [Google Scholar] [CrossRef]
  7. Kadlec, P.; Gabrys, B.; Strandt, S. Data-driven Soft Sensors in the process industry. Comput. Chem. Eng. 2009, 33, 795–814. [Google Scholar] [CrossRef] [Green Version]
  8. Pisa, I.; Santín, I.; Vicario, J.; Morell, A.; Vilanova, R. ANN-Based Soft Sensor to Predict Effluent Violations in Wastewater Treatment Plants. Sensors 2019, 19, 1280. [Google Scholar] [CrossRef] [Green Version]
  9. Habibzadeh, H.; Boggio-Dandry, A.; Qin, Z.; Soyata, T.; Kantarci, B.; Mouftah, H.T. Soft Sensing in Smart Cities: Handling 3Vs Using Recommender Systems, Machine Intelligence, and Data Analytics. IEEE Commun. Mag. 2018, 56, 78–86. [Google Scholar] [CrossRef]
  10. Haimi, H.; Mulas, M.; Corona, F.; Vahala, R. Data-derived soft-sensors for biological wastewater treatment plants: An overview. Environ. Model. Softw. 2013, 47, 88–107. [Google Scholar] [CrossRef]
  11. Szeląg, B. Modelowanie Matematyczne Optymalizacja i Sterowanie Pracą Przepływowych Oczyszczalni Ścieków; Monografia IBS PAN: Warszawa, Poland, 2019. [Google Scholar]
  12. Fernandez de Canete, J.; Del Saz-Orozco, P.; Baratti, R.; Mulas, M.; Ruano, A.; Garcia-Cerezo, A. Soft-sensing estimation of plant effluent concentrations in a biological wastewater treatment plant using an optimal neural network. Expert Syst. Appl. 2016, 63, 8–19. [Google Scholar] [CrossRef]
  13. Luccarini, L.; Porrà, E.; Spagni, A.; Ratini, P.; Grilli, S.; Longhi, S.; Bortone, G. Soft sensors for control of nitrogen and phosphorus removal from wastewaters by neural networks. Water Sci. Technol. 2002, 45, 101–107. [Google Scholar] [CrossRef] [PubMed]
  14. Hernández-del-Olmo, F.; Gaudioso, E.; Duro, N.; Dormido, R. Machine Learning Weather Soft-Sensor for Advanced Control of Wastewater Treatment Plants. Sensors 2019, 19, 3139. [Google Scholar] [CrossRef] [Green Version]
  15. Yu, P.; Cao, J.; Jegatheesan, V.; Du, X. A Real-time BOD Estimation Method in Wastewater Treatment Process Based on an Optimized Extreme Learning Machine. Appl. Sci. 2019, 9, 523. [Google Scholar] [CrossRef] [Green Version]
  16. Łagód, G.; Duda, S.M.; Majerek, D.; Szutt, A.; Dołhańczuk-Śródka, A. Application of Electronic Nose for Evaluation of Wastewater Treatment Process Effects at Full-Scale WWTP. Processes 2019, 7, 251. [Google Scholar] [CrossRef] [Green Version]
  17. Guz, Ł.; Łagód, G.; Jaromin-Gleń, K.; Suchorab, Z.; Sobczuk, H.; Bieganowski, A. Application of Gas Sensor Arrays in Assessment of Wastewater Purification Effects. Sensors 2014, 15, 1–21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Szulczyński, B.; Wasilewski, T.; Wojnowski, W.; Majchrzak, T.; Dymerski, T.; Namieśnik, J.; Gębicki, J. Different Ways to Apply a Measurement Instrument of E-Nose Type to Evaluate Ambient Air Quality with Respect to Odour Nuisance in a Vicinity of Municipal Processing Plants. Sensors 2017, 17, 2671. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Gebicki, J.; Byliński, H.; Namieśnik, J. Measurement techniques for assessing the olfactory impact of municipal sewage treatment plants. Environ. Monit. Assess. 2016, 188, 32. [Google Scholar] [CrossRef] [PubMed]
  20. Geng, Z.; Chen, G.; Han, Y.; Lu, G.; Li, F. Semantic relation extraction using sequential and tree-structured LSTM with attention. Inf. Sci. 2020, 509, 183–192. [Google Scholar] [CrossRef]
  21. Geng, Z.; Zhang, Y.; Li, C.; Han, Y.; Cui, Y.; Yu, B. Energy optimization and prediction modeling of petrochemical industries: An improved convolutional neural network based on cross-feature. Energy 2020, 194, 116851. [Google Scholar] [CrossRef]
  22. Han, Y.; Zhou, R.; Geng, Z.; Bai, J.; Ma, B.; Fan, J. A novel data envelopment analysis cross-model integrating interpretative structural model and analytic hierarchy process for energy efficiency evaluation and optimization modeling: Application to ethylene industries. J. Clean. Prod. 2020, 246, 118965. [Google Scholar] [CrossRef]
  23. Henze, M.; Harremoes, P.; Arvin, E.; Lacour, J. Wastewater Treatment. Biological and Chemical Processes; Springer: Berlin/Heidelberg, Germany, 2002. [Google Scholar]
  24. Lou, I.; Zhao, Y. Sludge Bulking Prediction Using Principle Component Regression and Artificial Neural Network. Math. Probl. Eng. 2012, 2012, 1–17. [Google Scholar] [CrossRef] [Green Version]
  25. Boztoprak, H.; Özbay, Y.; Güçlü, D.; Küçükhemek, M. Prediction of sludge volume index bulking using image analysis and neural network at a full-scale activated sludge plant. Desalin. Water Treat. 2016, 57, 17195–17205. [Google Scholar] [CrossRef]
  26. Han, H.-G.; Li, Y.; Guo, Y.-N.; Qiao, J.-F. A soft computing method to predict sludge volume index based on a recurrent self-organizing neural network. Appl. Soft Comput. 2016, 38, 477–486. [Google Scholar] [CrossRef]
  27. Sutha, K.; Tamilselvi Jebamalar, J. Machine Learning Based Supervised Feature Selection Algorithm for Data Mining. Int. J. Innov. Technol. Explor. Eng. 2019, 8, 3396–3401. [Google Scholar]
  28. Mageswary, G.; Karthikeyan, M. Statistical based Feature Selection and Ensemble Model for Network Intrusion Detection using Data Mining Technique. Int. J. Recent Technol. Eng. 2019, 8, 858–864. [Google Scholar]
  29. Comas, J.; Rodríguez-Roda, I.; Gernaey, K.V.; Rosen, C.; Jeppsson, U.; Poch, M. Risk assessment modelling of microbiology-related solids separation problems in activated sludge systems. Environ. Model. Softw. 2008, 23, 1250–1261. [Google Scholar] [CrossRef]
  30. Aarnio, P.; Minkkinen, P. Application of partial least-squares modelling in the optimization of a waste-water treatment plant. Anal. Chim. Acta 1986, 191, 457–460. [Google Scholar] [CrossRef]
  31. Capodaglio, A. Sludge bulking analysis and forecasting: Application of system identification and artificial neural computing technologies. Water Res. 1991, 25, 1217–1224. [Google Scholar] [CrossRef]
  32. Cortés, U.; Martínez, M.; Comas, J.; Sànchez-Marrè, M.; Poch, M.; Rodríguez-Roda, I. A conceptual model to facilitate knowledge sharing for bulking solving in wastewater treatment plants. AI Commun. 2003, 16, 279–289. [Google Scholar]
  33. Flores-Alsina, X.; Comas, J.; Roda, I.R.; Poch, M.; Gernaey, K.V.; Jeppsson, U. Evaluation of plant-wide WWTP control strategies including the effects of filamentous bulking sludge. Water Sci. Technol. 2009, 60, 2093–2103. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Flores-Alsina, X.; Comas, J.; Rodriguez-Roda, I.; Gernaey, K.V.; Rosen, C. Including the effects of filamentous bulking sludge during the simulation of wastewater treatment plants using a risk assessment model. Water Res. 2009, 43, 4527–4538. [Google Scholar] [CrossRef] [PubMed]
  35. Deepnarain, N.; Nasr, M.; Kumari, S.; Stenström, T.A.; Reddy, P.; Pillay, K.; Bux, F. Decision tree for identification and prediction of filamentous bulking at full-scale activated sludge wastewater treatment plant. Process Saf. Environ. Prot. 2019, 126, 25–34. [Google Scholar] [CrossRef]
  36. Vitorino, D.; Coelho, S.T.; Santos, P.; Sheets, S.; Jurkovac, B.; Amado, C. A Random Forest Algorithm Applied to Condition-based Wastewater Deterioration Modeling and Forecasting. Procedia Eng. 2014, 89, 401–410. [Google Scholar] [CrossRef] [Green Version]
  37. Youssef, A.M.; Pourghasemi, H.R.; Pourtaghi, Z.S.; Al-Katheeri, M.M. Landslide susceptibility mapping using random forest, boosted regression tree, classification and regression tree, and general linear models and comparison of their performance at Wadi Tayyah Basin, Asir Region, Saudi Arabia. Landslides 2016, 13, 839–856. [Google Scholar] [CrossRef]
  38. Yurtsever, M.; Yurtsever, U. Use of a convolutional neural network for the classification of microbeads in urban wastewater. Chemosphere 2019, 216, 271–280. [Google Scholar] [CrossRef]
  39. Zhao, X.; Barber, S.; Taylor, C.C.; Milan, Z. Interval forecasts based on regression trees for streaming data. Adv. Data Anal. Classif. 2019. [Google Scholar] [CrossRef]
  40. Shih, Y.-S.; Liu, K.-H. Regression trees for detecting preference patterns from rank data. Adv. Data Anal. Classif. 2019, 13, 683–702. [Google Scholar] [CrossRef]
  41. Ma, X. Using Classification and Regression Trees: A Practical Primer; Information Age Publishing: Charlotte, NC, USA, 2018. [Google Scholar]
  42. Esmaily, H.; Tayefi, M.; Doosti, H.; Ghayour-Mobarhan, M.; Nezami, H.; Amirabadizadeh, A. A Comparison between Decision Tree and Random Forest in Determining the Risk Factors Associated with Type 2 Diabetes. J. Res. Health Sci. 2018, 18, 1–7. [Google Scholar]
  43. Barnard, D.M.; Germino, M.J.; Pilliod, D.S.; Arkle, R.S.; Applestein, C.; Davidson, B.E.; Fisk, M.R. Cannot see the random forest for the decision trees: Selecting predictive models for restoration ecology. Restor. Ecol. 2019, 27, 1053–1063. [Google Scholar] [CrossRef]
  44. Ahmad, M.W.; Reynolds, J.; Rezgui, Y. Predictive modelling for solar thermal energy systems: A comparison of support vector regression, random forest, extra trees and regression trees. J. Clean. Prod. 2018, 203, 810–821. [Google Scholar] [CrossRef]
  45. Breiman, L. Random Forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef] [Green Version]
  46. Friedman, J.H. Stochastic gradient boosting. Comput. Stat. Data Anal. 2002, 38, 367–378. [Google Scholar] [CrossRef]
  47. Paul, A.; Mukherjee, D.P.; Das, P.; Gangopadhyay, A.; Chintha, A.R.; Kundu, S. Improved Random Forest for Classification. IEEE Trans. Image Process. 2018, 27, 4012–4024. [Google Scholar] [CrossRef] [PubMed]
  48. Senagi, K.; Jouandeau, N. A Non-Deterministic Strategy for Searching Optimal Number of Trees Hyperparameter in Random Forest. In 2018 Federated Conference on Computer Science and Information Systems; IEEE: Piscataway, NJ, USA, 2018; pp. 73–80. [Google Scholar]
  49. Güçlü, D.; Dursun, S. Artificial neural network modelling of a large-scale wastewater treatment plant operation. Bioprocess Biosyst. Eng. 2010, 33, 1051–1058. [Google Scholar] [CrossRef]
  50. Hastie, T.; Tibshirani, R.; Friedman, J. The Elements of Statistical Learning; Springer Series in Statistics; Springer: New York, NY, USA, 2009. [Google Scholar]
  51. Rumelhart, D.E.; Hinton, G.E.; Williams, R.J. Learning Internal Representations by Error Propagation. In Parallel Distributed Processing; Rumelhart, D.E., Mcclelland, J.L., Eds.; MIT Press: Cambridge, MA, USA, 1986; Volume 1, pp. 318–362. [Google Scholar]
  52. Haykin, S. Neural Networks: A Comprehensive Foundation, 2nd ed.; Prentice Hall PTR: Upper Saddle River, NJ, USA, 1998. [Google Scholar]
  53. Hecht-Nielsen, R. Kolmogorov’s mapping neural network existence theorem. In Proceedings of the 1st IEEE International Joint Conference of Neural Networks; Academic Press: New York, NY, USA, 1987; Volume 3, pp. 11–14. [Google Scholar]
  54. Burges, C. A Tutorial on Support Vector Machines for Pattern Recognition. Data Min. Knowl. Discov. 1998, 2, 121–167. [Google Scholar]
  55. Vapnik, V. The Support Vector Method of Function Estimation. In Nonlinear Modeling; Springer: Boston, MA, USA, 1998; pp. 55–85. [Google Scholar]
  56. Han, H.; Qiao, J. Hierarchical Neural Network Modeling Approach to Predict Sludge Volume Index of Wastewater Treatment Process. IEEE Trans. Control Syst. Technol. 2013, 21, 2423–2431. [Google Scholar] [CrossRef]
  57. Honggui, H.; Ying, L.; Junfei, Q. A fuzzy neural network approach for online fault detection in waste water treatment process. Comput. Electr. Eng. 2014, 40, 2216–2226. [Google Scholar] [CrossRef]
  58. Kusiak, A.; Zeng, Y.; Zhang, Z. Modeling and analysis of pumps in a wastewater treatment plant: A data-mining approach. Eng. Appl. Artif. Intell. 2013, 26, 1643–1651. [Google Scholar] [CrossRef]
  59. Asadi, A.; Verma, A.; Yang, K.; Mejabi, B. Wastewater treatment aeration process optimization: A data mining approach. J. Environ. Manag. 2017, 203, 630–639. [Google Scholar] [CrossRef]
  60. Zaborowska, E.; Czerwionka, K.; Makinia, J. Strategies for achieving energy neutrality in biological nutrient removal systems—A case study of the Slupsk WWTP (northern Poland). Water Sci. Technol. 2017, 75, 727–740. [Google Scholar] [CrossRef]
  61. Flores-Alsina, X.; Arnell, M.; Amerlinck, Y.; Corominas, L.; Gernaey, K.V.; Guo, L.; Lindblom, E.; Nopens, I.; Porro, J.; Shaw, A.; et al. Balancing effluent quality, economic cost and greenhouse gas emissions during the evaluation of (plant-wide) control/operational strategies in WWTPs. Sci. Total Environ. 2014, 466, 616–624. [Google Scholar] [CrossRef] [PubMed]
  62. Bayo, J.; Angosto, J.M.; Serrano-Aniorte, J.; Cascales-Pujalte, J.A.; Fernández-López1, C.; López-Castellanos, J. Evaluation of physicochemical parameters influencing bulking episodes in a municipal wastewater treatment plant. In Water Pollution VIII: Modelling, Monitoring and Management; WIT Transactions on Ecology and the Environment, Volume 95; WIT Press: Southampton, UK, 2006; Volume 1, pp. 531–540. [Google Scholar]
  63. Feuerman, M.; Miller, A.R. Relationships between statistical measures of agreement: Sensitivity, specificity and kappa. J. Eval. Clin. Pract. 2008, 14, 930–933. [Google Scholar] [CrossRef] [PubMed]
  64. Procházka, A.; Vyšata, O.; Ťupa, O.; Yadollahi, M.; Vališ, M. Discrimination of axonal neuropathy using sensitivity and specificity statistical measures. Neural Comput. Appl. 2014, 25, 1349–1358. [Google Scholar] [CrossRef]
  65. Farmer, E.J. Detecting Leaks in Pipelines; International Society of Automation: Research Triangle Park, NC, USA, 2018. [Google Scholar]
  66. Martins, A.M.P.; Heijnen, J.J.; van Loosdrecht, M.C.M. Bulking sludge in biological nutrient removal systems. Biotechnol. Bioeng. 2004, 86, 125–135. [Google Scholar] [CrossRef] [PubMed]
  67. Kowalska, E.; Paturej, E.; Zielińska, M. Use of Lecane inermis for control of sludge bulking caused by the Haliscomenobacter genus. Desalin. Water Treat. 2016, 57, 10916–10923. [Google Scholar] [CrossRef]
  68. Bezak-Mazur, E.; Stoińska, R.; Szeląg, B. Evaluation of the impact of operational parameters and particular filamentous bacteria on activated sludge volume index—Case study. Ann. Set. Environ. Prot. 2016, 18, 487. [Google Scholar]
  69. Szeląg, B.; Barbusiński, K.; Studziński, J. Application of the model of sludge volume index forecasting to assess reliability and improvement of wastewater treatment plant operating conditions. Desalin. Water Treat. 2019, 140, 143–154. [Google Scholar] [CrossRef]
  70. Alex, J.; Benedetti, L.; Copp, J.; Gernaey, K.V.; Jeppsson, U.; Nopens, I.; Pons, M.-N.; Steyer, J.-P.; Vanrolleghem, P.A. Benchmark Simulation Model No. 1 (BSM1); Lund University: Lund, Sweden, 2008; Volume TEIE-7229. [Google Scholar]
  71. Corominas, L.; Larsen, H.F.; Flores-Alsina, X.; Vanrolleghem, P.A. Including Life Cycle Assessment for decision-making in controlling wastewater nutrient removal systems. J. Environ. Manag. 2013, 128, 759–767. [Google Scholar] [CrossRef] [Green Version]
  72. Kim, M.; Kim, Y.; Kim, H.; Piao, W.; Kim, C. Evaluation of the k-nearest neighbor method for forecasting the influent characteristics of wastewater treatment plant. Front. Environ. Sci. Eng. 2016, 10, 299–310. [Google Scholar] [CrossRef]
  73. Wilimowska, Z.; Borzemski, L.; Świątek, J. (Eds.) Information Systems Architecture and Technology: Proceedings of 38th International Conference on Information Systems Architecture and Technology—ISAT 2017. In Advances in Intelligent Systems and Computing; Springer International Publishing: Cham, Switzerland, 2018; Volume 657. [Google Scholar]
  74. Ahnert, M.; Marx, C.; Krebs, P.; Kuehn, V. A black-box model for generation of site-specific WWTP influent quality data based on plant routine data. Water Sci. Technol. 2016, 74, 2978–2986. [Google Scholar] [CrossRef] [Green Version]
  75. Jurik, L.; Kaletova, T.; Sedmakova, M.; Balazova, P.; Cervenanska, A. Comparison of service characteristics of two town’s WWTP. J. Ecol. Eng. 2017, 18, 61–67. [Google Scholar] [CrossRef] [Green Version]
  76. Rössle, W.H.; Pretorius, W.A. Batch and automated SVI measurements based on short-term temperature variations. Water SA 2019, 34, 237. [Google Scholar] [CrossRef] [Green Version]
  77. Bagheri, M.; Mirbagheri, S.A.; Bagheri, Z.; Kamarkhani, A.M. Modeling and optimization of activated sludge bulking for a real wastewater treatment plant using hybrid artificial neural networks-genetic algorithm approach. Process Saf. Environ. Prot. 2015, 95, 12–25. [Google Scholar] [CrossRef]
  78. Han, H.-G.; Qiao, J.-F. Prediction of activated sludge bulking based on a self-organizing RBF neural network. J. Process Control 2012, 22, 1103–1112. [Google Scholar] [CrossRef]
  79. Mirbagheri, S.A.; Bagheri, M.; Boudaghpour, S.; Ehteshami, M.; Bagheri, Z. Performance evaluation and modeling of a submerged membrane bioreactor treating combined municipal and industrial wastewater using radial basis function artificial neural networks. J. Environ. Health Sci. Eng. 2015, 13, 17. [Google Scholar] [CrossRef] [Green Version]
  80. Pham, B.T.; Pradhan, B.; Tien Bui, D.; Prakash, I.; Dholakia, M.B. A comparative study of different machine learning methods for landslide susceptibility assessment: A case study of Uttarakhand area (India). Environ. Model. Softw. 2016, 84, 240–250. [Google Scholar] [CrossRef]
  81. Dreiseitl, S.; Ohno-Machado, L. Logistic regression and artificial neural network classification models: A methodology review. J. Biomed. Inform. 2002, 35, 352–359. [Google Scholar] [CrossRef] [Green Version]
  82. Kusiak, A.; Wei, X. Prediction of methane production in wastewater treatment facility: A data-mining approach. Ann. Oper. Res. 2014, 216, 71–81. [Google Scholar] [CrossRef]
Figure 1. A schematic of the soft sensor system for selecting a data mining method for simulating sludge bulking.
Figure 1. A schematic of the soft sensor system for selecting a data mining method for simulating sludge bulking.
Sensors 20 01941 g001
Figure 2. Matrix of the cases of interest and the weight factor values W g q . where the following relationships are satisfied: λ = f(f1, f2, f3 …, fg) = f(λ), δ = f(z1, z2, z3 …, zq) = f(δ).
Figure 2. Matrix of the cases of interest and the weight factor values W g q . where the following relationships are satisfied: λ = f(f1, f2, f3 …, fg) = f(λ), δ = f(z1, z2, z3 …, zq) = f(δ).
Sensors 20 01941 g002
Figure 3. Matrices for the selection of methods for identification of sludge bulking.
Figure 3. Matrices for the selection of methods for identification of sludge bulking.
Sensors 20 01941 g003
Figure 4. The effect of variation of the MLSS and DO values on the sensitivity coefficient SMLSS,DO.
Figure 4. The effect of variation of the MLSS and DO values on the sensitivity coefficient SMLSS,DO.
Sensors 20 01941 g004
Figure 5. The effect of variation of the T and DO values on the sensitivity coefficient ST,DO.
Figure 5. The effect of variation of the T and DO values on the sensitivity coefficient ST,DO.
Sensors 20 01941 g005
Figure 6. The effect of variation of the MLSS and TN values on the sensitivity coefficient SMLSS,TN.
Figure 6. The effect of variation of the MLSS and TN values on the sensitivity coefficient SMLSS,TN.
Sensors 20 01941 g006
Figure 7. The effect of variation of the τ and TN values on the sensitivity coefficient Sτ, TN.
Figure 7. The effect of variation of the τ and TN values on the sensitivity coefficient Sτ, TN.
Sensors 20 01941 g007
Table 1. Results of measurements of the volume and quality of wastewater and the reactor operating parameters in the Sitkówka-Nowiny wastewater treatment plant (WWTP).
Table 1. Results of measurements of the volume and quality of wastewater and the reactor operating parameters in the Sitkówka-Nowiny wastewater treatment plant (WWTP).
IndicatorsUnitsWinterSpring-Fall
MinMeanMaxStandard DeviationMinMeanMaxStandard Deviation
Qm3/d29,95239,36488,986656330,12541,84294,7728559
BODmgO2/l15129048981.8313234055781.2
CODmg O2/l3847821183161.43428201703178.2
TSSmg/l13631547462.7611035057289.4
N-NH4mg/l2848.9625.682254.5266.97.13
TNmg/l56.282.0195.168.4239.995.15124.111.58
TPmg/l3.17.2212.11.443.57.8312.61.65
T°C1011.913.50.811.317.8233.1
DOmg/l1.82.853.250.81.512.23.250.65
MLSSmg/l2.854.956.540.842.154.115.280.95
MLSSRmg/l6.628.714.920.725.037.8111.860.1
WASt.s.m/d12.6915.3518.353.5110.0212.3517.253.77
RAS%85.2102.915216.2575.283.06120.524.4
SVIcm3/g154198291359013820037
mPIXm3/d00.811.750.2700.841.820.28
where: Q—inflow to wastewater treatment plant, BOD, COD—biochemical/chemical oxygen demand, TSS—total suspended solids, N–NH4—ammonium nitrogen, TN—total nitrogen, TP—total phosphorus, T—temperature in the activated sludge chambers, DO—oxygen concentration in the bioreactor, MLSS—concentration of activated sludge, MLSSR—concentration of return sludge, WAS—excess amount of sludge, RAS—degree of recirculation, mPIX—dosage of PIX, SVI—sludge volumetric index.
Table 2. List of values of the measures of fit between the calculated and the measured results (testing data) (SENS– sensitivity, SPEC–specifity) for SVI, obtained with various methods for a combination of independent variables (quantity and quality of wastewater, bioreactor operating parameters).
Table 2. List of values of the measures of fit between the calculated and the measured results (testing data) (SENS– sensitivity, SPEC–specifity) for SVI, obtained with various methods for a combination of independent variables (quantity and quality of wastewater, bioreactor operating parameters).
IndependentBTLRMLPSVMRF
SPSEN.T.SPSESPSEN:H:ESPSEC, γSPSEN.T.
A0.650.81190.270.880.420.875:exp,lin0.380.92900,0.350.460.88170
B0.880.81880.350.940.310.964: sin,lin0.310.98100,0.400.350.8860
C0.620.83190.270.920.460.883: lin,lin0.500.851000,0.250.460.8820
A,B0.850.83190.380.900.480.928:exp,tanh0.420.98800; 0.350.001.0025
A/B0.770.73110.270.850.350.943: tanh, lin0.540.85800;0.250.580.8140
A/C0.730.75380.350.850.350.923:sin,lin0.500.90800;0.350.690.7940
A,C0.650.75110.270.850.310.943:tanh,lin0.540.87800;0.800.580.7940
D0.810.881830.650.870.730.853:exp, sin0.620.88900;0.330.880.7730
D,A/B0.880.77100.680.870.650.874:lin,exp0.730.85800;0.350.880.8140
D,A/C0.880.881550.680.880.690.853:lin,lin0.730.87800;0.330.880.8030
D,B0.880.77100.620.880.650.907:tanh,lin0.730.87900;0.330.880.8210
D,B,C0.880.80240.620.880.680.906:tanh,lin0.760.88800;0.250.880.8630
D,A,B,C0.920.83140.650.920.690.949:lin,sin0.770.901000;0.250.850.8730
E0.880.751190.620.830.880.925:log,exp0.850.88900;0.400.500.98119
E,D0.950.96980.880.960.921.005:exp,lin1.000.98900;0.250.880.9230
E,F0.960.87980.650.900.880.965:lin,lin0.880.96900;0.200.461.00200
F0.920.811950.001.000.460.904:tanh,tanh0.580.85900;0.500.920.7130
D,F0.880.77150.620.880.730.857:log,lin0.650.88700;0.330.880.7730
D,F,B0.920.90150.650.900.690.907:exp, lin0.690.90700;0.250.850.8530
D,F,A/B0.900.8450.700.880.850.796:tanh,lin0.810.87900;0.330.880.7730
D,F,B,C0.880.81100.650.900.690.928:tanh,tanh0.810.88800;0.250.880.8530
F,B0.850.881140.350.920.500.873:tanh,tanh0.460.96700;0.200.880.77140
F,B,C0.920.921860.310.940.550.965:sin,lin0.540.85500;0.350.540.8520
F,C0.850.73150.230.900.500.883:lin,lin0.540.87600;0.400.500.7530
D,F,E1.000.931620.880.960.960.987:exp,lin0.960.98600;0.250.081.0030
A/C,F,E1.000.881190.690.900.880.947:lin,lin0.880.94600;0.350.851.0030
A/C,F,E,D1.000.961950.920.960.910.957: lin,sin0.931.00600;0.200.461.0060
A/C,A/B,F,E,D1.000.971000.940.960.921.003:tanh,exp0.940.98200;0.330.881.00100
B,C,E,F,D1.000.961190.960.940.951.008:tanh,lin0.880.98100;0.170.811.0070
B,C,E,F,D,G1.000.981950.960.980.991.004: lin,lin0.920.98100;0.250.880.9670
B,C,E,F0.960.881950.850.940.930.994: lin,lin0.690.96200;0.350.660.9170
B,E,F,D,G1.000.961950.850.940.951.006:lin,tanh0.920.9850;0.250.850.9230
B,E,F,D1.000.941950.960.960.971.005: tanh,log0.920.9570;0.200.731.0030
A/C,A/B,F,E,D,G1.000.981950.970.961.000.989:exp,lin0.980.99100;0.140.931.0030
A/C,A/B,F,E,G1.000.901190.730.880.810.964:lin,lin0.880.94900;0.250.850.96170
F,E,G0.960.881190.650.900.780.968:sin,exp0.810.94700;0.350.770.9430
F,E,G,D1.000.961990.960.960.960.967:exp,lin1.000.94700;0.20.811.0040
where: A—BOD, B—TN, C—TP, D—T, E—MLSS, F—DO, G—mPIX, SP—SPEC, SE—SENS.
Table 3. A matrix of cases showing relationships between the values of the functions f(δ), f(λ) and the classification models for predicting activated sludge bulking.
Table 3. A matrix of cases showing relationships between the values of the functions f(δ), f(λ) and the classification models for predicting activated sludge bulking.
δ/λ 01/31/31/3 (2/3)2/32/33/3
f(δ)/f(λ)01/62/63/64/65/66/6
00XBTBTBTBTBTBT
1/31/6RFBTRFRFRFRFBT
1/32/6BTMLPBTBTBTBTBT
1/3(2/3)3/6BTBTMLPBTBTBTMLP
2/34/6SVMBTMLPSVMBTMLPBT
2/35/6SVMMLPMLPBTMLPMLPBT
3/36/6SVMSVMMLPBTSVMBTBT
Table 4. A matrix of cases showing relationships between the values of the functions f(δ), f(λ) and the classification models, with a particular focus on the logit for predicting activated sludge bulking.
Table 4. A matrix of cases showing relationships between the values of the functions f(δ), f(λ) and the classification models, with a particular focus on the logit for predicting activated sludge bulking.
δ/λ 01/31/31/3 (2/3)2/32/33/3
f(δ)/f(λ)01/62/63/64/65/66/6
00X BTBTBTBTBTBT
1/31/6RFBTRFRFRFRFBT
1/32/6BTMLPBTBTBTBTBT
1/3(2/3)3/6LRVGLRVGLRVGLRVGLRVGLRVGLRVG
2/34/6LRVGLREXLREXLREXLREXLREXLREX
2/35/6SVMLRVGLRVGLRVGLRVGLRVGLREX
3/36/6LRVGSVMMLPBTSVMLREXLREX
where: VG—the logit model, where SPEC, SENS >0.8, EX—the logit model, where SPEC, SENS >0.9.
Table 5. A matrix of cases showing the relationships between the values of the functions f(δ), f(λ), f(t) and the classification soft sensor models for predicting the activated sludge bulking.
Table 5. A matrix of cases showing the relationships between the values of the functions f(δ), f(λ), f(t) and the classification soft sensor models for predicting the activated sludge bulking.
f(t) 0100101
δ/λ 01/31/31/3 (2/3)2/32/33/3
f(δ)/f(λ)01/62/63/64/65/66/6
00 X XBTBT XBT X
1/31/6RF XRFRF XRF X
1/32/6BT XBTBT XBT X
1/3(2/3)3/6BT XMLPBT XBT X
2/34/6SVM XMLPSVM XMLP X
2/35/6SVM XMLPBT XMLP X
3/36/6SVM XMLPBT XBT X

Share and Cite

MDPI and ACS Style

Szeląg, B.; Drewnowski, J.; Łagód, G.; Majerek, D.; Dacewicz, E.; Fatone, F. Soft Sensor Application in Identification of the Activated Sludge Bulking Considering the Technological and Economical Aspects of Smart Systems Functioning. Sensors 2020, 20, 1941. https://doi.org/10.3390/s20071941

AMA Style

Szeląg B, Drewnowski J, Łagód G, Majerek D, Dacewicz E, Fatone F. Soft Sensor Application in Identification of the Activated Sludge Bulking Considering the Technological and Economical Aspects of Smart Systems Functioning. Sensors. 2020; 20(7):1941. https://doi.org/10.3390/s20071941

Chicago/Turabian Style

Szeląg, Bartosz, Jakub Drewnowski, Grzegorz Łagód, Dariusz Majerek, Ewa Dacewicz, and Francesco Fatone. 2020. "Soft Sensor Application in Identification of the Activated Sludge Bulking Considering the Technological and Economical Aspects of Smart Systems Functioning" Sensors 20, no. 7: 1941. https://doi.org/10.3390/s20071941

APA Style

Szeląg, B., Drewnowski, J., Łagód, G., Majerek, D., Dacewicz, E., & Fatone, F. (2020). Soft Sensor Application in Identification of the Activated Sludge Bulking Considering the Technological and Economical Aspects of Smart Systems Functioning. Sensors, 20(7), 1941. https://doi.org/10.3390/s20071941

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