Application of the Mathematical Simulation Methods for the Assessment of the Wastewater Treatment Plant Operation Work Reliability

The aim of the present work was the modeling of the wastewater treatment plant operation work using Monte Carlo method and different random variables probability distributions modeling. The analysis includes the following pollutants indicators; BOD5 (Biochemical Oxygen Demand), CODCr (Chemical Oxygen Demand), Total Suspended Solids (SSt), Total Nitrogen (TN), and Total Phosphorus (TP). The Anderson–Darling (A–D) test was used for the assessment of theoretical and empirical distributions compatibility. The selection of the best-fitting statistical distributions was performed using peak-weighted root mean square (PWRMSE) parameter. Based on the performed calculations, it was stated that pollutants indicators in treated sewage were characterized by a significant variability. Obtained results indicate that the best-fitting pollutants indicators statistical distribution is Gauss Mixed Model (GMM) function. The results of the Monte Carlo simulation method confirmed that some problems related to the organic and biogenic pollutants reduction may be observed in the Wastewater Treatment Plant, in Jaworzno.


Introduction
The dynamic development of water supply and sewage systems that has been observed in Poland for many years contributes to the discharge of a large amount of sewage.The sewage is defined as the community's spent water.Although it mainly consists of pure water (over 98%), it contains waste of almost every form and description, with the remainder being dirt.Approximately 25% of the pollutants in normal domestic sewage are in suspension and 75% in solution.Sewage contains many complex organic and mineral compounds.The organic portion of sewage is biochemically degradable and, as such, is responsible for the offensive characteristics usually associated with sewage.Furthermore, sewage contains large numbers of microorganisms, most of which are bacteria [1].Before the discharging of treated sewage into the natural sewage receivers, there is a necessity for reach the required level of sewage treatment, because of the natural environment [2][3][4].Because of this, it is important to make a general assessment each of the wastewater treatment plant, both in terms of the effectiveness of pollutants reduction and the reliability of the operation work of wastewater treatment plant facilities.Extensive analysis of the research results allows for predict the improprieties and for its elimination in the future.This is a significant point for attaining effective surface water protection against pollution [5].
Surface waters are very important source of water for large urban agglomerations.Continuous urban agglomerations development leads to an increase in human impact on the aquatic environment.Pollutants in rivers, lakes, and coastal waters pose a significant threat to the water reservoirs users, because of discharging the sewage after treatment processes into the natural environment.In extreme cases, pollutants that are introduced with sewage into the water can lead to the extinction all of the biological life forms in water and for purging the sewage natural receiver, e.g., surface water or ground water that have self-cleaning skills [6][7][8].In recent decades, some problems related to the deficiency of required quality water have intensified.Rapid urbanization and the increase in the number of residents contributed to the surface water quality deterioration.In addition, due to the still growing water demand and its limited amount, there is a necessity for water recover through the wastewater treatment plant processes [9][10][11].
In order to achieve the required level of pollutants indicators in treated sewage, it is necessary to control the operational work of the particular wastewater treatment plant facilities.It is connected with the permanent or periodic complex control of the sewage treatment processes.The aim of control of the wastewater treatment plant facilities is due to the necessity of its maintaining in permanent readiness to the proper work, which provides the required treated sewage quality at the outflow.The growth of the wastewater treatment plant technological reliability, affects the quick detection of devices failures and measuring apparatus failures.In addition, a technological control including daily observations of the treated sewage quality and quantity, the condition of activated sludge and other parameters, affects the sewage treatment process effects in a special way.
Generally, the assessment of the operation work of the wastewater treatment plant is performed based on some pollutants indicator concentrations in treated sewage, which are compared with legally required pollutants indicators values.These limitation values are determined by the obligatory legal acts, such as the water and legal permission acts, and other regulations concerning the quality of treated sewage.However, the method of assessing the operation work of the wastewater treatment plant that is presented above is a bit simplified.This is because of the basis only on the empirical sample of the analyzed random variables.The reliability of the wastewater treatment plant is mainly considered as the percentage share of time, during which the expected treated sewage content meets the requirements [12].A variability of the pollutants indicators values in sewage should be considered as a function of the probability distribution that determines the physicochemical parameters of treated sewage content.Thus, the mathematical simulations for the wastewater treatment plant operation work should be conducted.It allows for make more reliable assessment of the operation work of the wastewater treatment plant in comparison to the usually used methods that are based on short observations series for pollutants indicators.It should be noted that the using of a mathematical simulations for the optimization of sewage treatment processes is justified.Mathematical simulations give the possibility for elaboration the dynamic mathematical models of sewage treatment processes.Mathematical models allow making a decision about the changes in the sewage treatment processes.An additional advantage of using a mathematical models is the possibility of considering the variability of inflowing sewage quality (parameters).The wastewater treatment plant operator makes an assessment of the wastewater treatment plant operation work only based on the analysis of the pollutants indicators in treated sewage.Sometimes, the values of pollutants indicators in sewage at the inflow and sewage at the outflow from the wastewater treatment plant are compared with each other.Based on this comparison, it is possible to determine the reliability of the wastewater treatment plant operation work that is expressed by the reduction degree (%).Based on the conducted simulations, wastewater treatment plant operator is able to estimate the risk, i.e., the probability of occurrence in the future pollutants indicators values greater than its permissible limits.Thus, the higher risk value, the greater the need for the modernization or extension of the wastewater treatment plant facilities.A prediction of such values is even more important, because currently observed climate changes, demographic changes, and water consumption changes can affect significantly the changes of the wastewater treatment plant load, both in terms of the sewage volume (hydraulic load) and pollution sewage load.Finally, it can affect the treated sewage quality.Because of this, a prediction seems to be very important [13].Simulation issues of technological processes of sewage treatment are often considered research problems, both during designing process of the wastewater treatment plant facilities, its modernization, and exploitation too.Because the operation work of the wastewater treatment plant is basis on the complexity technology processes, many different models for describing its operation work may be found.Modeling of the wastewater treatment plant operation work can be performed by using Artificial Neural Networks (ANN), where the number of input variables can be separated by cluster analysis or by principal component analysis [14][15][16].Moreover, Activated Sludge Simulation Models (ASIM) are widely discussed in literature, e.g., Snip et al. [17], Wu et al. [18], Machado et al. [19], Alikhani et al. [20], Guo and Vanrolleghem [21].Modeling of the wastewater treatment plant operation work can be conducted using particular types of probability distributions.Bugajski et al. [22], Marzec [23], and Jóźwiakowski et al. [24] considered the possibility of using Weibull distribution for modeling the wastewater treatment plant operation work reliability.
The analysis of the operation work of the wastewater treatment plant using the probability distributions allows the description of the analyzed phenomena in a wider range.An important point of the presented method is the possibility of modeling complicated dependences between the considered variables and forecasting the pollutants indicators values in some period of time.Results obtained in this way may be helpful for the assessment of the risk of improper operation work of the wastewater treatment plant.However, it should be stressed that some assumptions about the unique form of pollutants indicators distribution may be incorrect.Thus, the stability of mechanisms having an impact on pollutants indicators reduction is suggested.In reality, there are many proposals basis on the same empirical sample.Because of this, for each empirical distribution of random variables, the best-fitting theoretical distribution should be determined.One of the methods used for selection of the best-fitting theoretical distribution is the Akaike (AIC), or Schwartz (BIC), information criteria [25].Alternative statistical tests also could be used for the assessment of theoretical distributions performance.One of them is an Anderson-Darling (A-D) test that compares the whole range, but it gives more weight to the upper tail [26].The main disadvantage of these goodness-of-fit metrics is that they not provide any information about the precision limit and the inclusion of untypical (outlying) values.Moreover, it is not possible to assess the values of test statistics.But on the other hand, some methods that are commonly used for the assessing of rainfall-runoff hydrological models, such as percentage error in peak flow, percentage error in volume, efficiency coefficient, peak-weighted root mean square error, sum of absolute residuals, or sum of squared residuals, should be considered.Also, for the verification of the best-fitting probability distributions, the other criterions, such as the Akaike criterion (AIC) or Bayesian information criterion (BIC), can be used.However these criterions also have some limitations, e.g., the necessity for the estimation of distributions parameters using only maximum likelihood method.
Based on the some information presented above, the aim of the study was determined.In this study, simulation of the wastewater treatment plant operation work using selected statistical distributions was performed.It will allow the determination of the wastewater treatment plant reliability coefficients, not only in the range of the empirical sample, but also based on the forecasting the values of pollutants indicators.The novelty of this study is the identification of the best-fitting theoretical function by adopting Peak-Weighted Root Mean Square Error (PWRMSE), which is commonly used for calibrating the parameters in hydrological models.Finally, the performed analysis aimed to determine the wastewater treatment plant operation work reliability coefficients.

Research Object
The wastewater treatment plant that is considered in this study is located in Jaworzno, Silesian Voivodeship, Poland.The Wastewater treatment plant (WTP) in Jaworzno combines a mechanical and biological sewage treatment processes using an activated sludge method in a biological phase.Designed capacity of the object is 25 000 m 3 •d −1 .Total population equivalent (p.e.) for WTP in Jaworzno is 125 000.Mechanical facilities of the WTP in Jaworzno include bar screens, grit chamber, and primary radial settling tank.Biological sewage treatment is performed in bioreactor using the activated sludge processes.Before the treated sewage is discharged into the natural receiver, sewage is directed into the two secondary settling tanks in order to separate the activated sludge from treated sewage.After this, treated sewage is discharged into the Przemsza River (Figure 1).

Materials and Methods
The According to the Polish obligatory standard [27], for the wastewater treatment plant with over 100,000 p.e., the permissible value of BOD5 in treated sewage is 15 mg O2•dm −3 , CODCr -125 mg O2•dm −3 , SSt -35 mg•dm −3 , TN -10 mg TN•dm −3 , and TP -1 mg TP•dm −3 .Moreover, the Regulation of the Minister for Environment [27] determines the minimum amount of the collected sewage samples after treatment process per year.In this case, for the wastewater treatment plant that is characterized with over 100,000 p.e., there are 24 control sewage samples per year.As the cited Regulation (2014) indicates, at most, three of them do not have to meet the quality requirements.

Preliminary Data Analysis
Preliminary analysis of the pollutants indicators series includes the calculations of the descriptive statistics.Besides minimum (Min.),mean (Mean), and maximum (Max.)values, the measures of dispersion, such as standard deviation (s) and coefficient of variation (Cs) were calculated.

Statistical distributions Fitting
Based on the observational series of pollutants indicators in treated sewage, the best-fitting theoretical distribution of random variables to its empirical distribution was considered.The analysis includes General Extreme Values (GEV), Gaussian Mixture Model (GMM), Log-normal, Normal, Pareto, Rayleigh, Triangle, and Weibull theoretical distributions.Each of the f(x) functions is described as follows [28][29][30][31][32][33].

Materials and Methods
The analysis was made based on the data including the values of the following pollutants indicators; According to the Polish obligatory standard [27], for the wastewater treatment plant with over 100,000 p.e., the permissible value of BOD 5 in treated sewage is 15 mg O 2 •dm −3 , COD Cr -125 mg O 2 •dm −3 , SS t -35 mg•dm −3 , TN -10 mg TN•dm −3 , and TP -1 mg TP•dm −3 .Moreover, the Regulation of the Minister for Environment [27] determines the minimum amount of the collected sewage samples after treatment process per year.In this case, for the wastewater treatment plant that is characterized with over 100,000 p.e., there are 24 control sewage samples per year.As the cited Regulation (2014) indicates, at most, three of them do not have to meet the quality requirements.

Preliminary Data Analysis
Preliminary analysis of the pollutants indicators series includes the calculations of the descriptive statistics.Besides minimum (Min.),mean (Mean), and maximum (Max.)values, the measures of dispersion, such as standard deviation (s) and coefficient of variation (C s ) were calculated.

Statistical Distributions Fitting
Based on the observational series of pollutants indicators in treated sewage, the best-fitting theoretical distribution of random variables to its empirical distribution was considered.The analysis includes General Extreme Values (GEV), Gaussian Mixture Model (GMM), Log-normal, Normal, Pareto, Rayleigh, Triangle, and Weibull theoretical distributions.Each of the f (x) functions is described as follows [28][29][30][31][32][33].
GMM distribution: where, Log-normal distribution: where, σ Y ; µ Y -parameters of the log-normal distribution: variance and mean value, respectively.
Weibull distribution: The assessment of theoretical distributions compatibility with pollutants indicators empirical distributions was made using Anderson-Darling (A-D) test.This type of the statistical test is sensitive in the whole distribution function range.In this case, the possibility of detection of some differences between the distribution functions is greater.Thus, in comparison to other tests, the A-D test provides better assessment of the distribution compatibility.It test includes two opposed hypothesis: H 0 : data follow a specified distribution and H 1 : data not follow a specified distribution.The calculations were performed for the statistical significance level of α = 0.05.The Anderson-Darling (A-D) statistic is described as follows [34][35][36] A where, n-number of observations; F-theoretical cumulative distribution; X i -ordered data.

The Best Statistical Distributions Fitting
The best-fitting of theoretical function to the empirical random variables distribution was identified based on the metric used in the assessment, e.g., in hydrological modeling: peak-weighted root mean square error (PWRMSE).The studies concerning the possibilities of its using for the assessment of the best-fitting statistical distributions were considered by Mły ński et al. [37] preliminary.The analyzed goodness-of-fit metric is described by the following formula [38,39].
where, n-size of the observation series; O i -observed value; P i -predicted value; Also, the Anderson-Darling statistic was used to compare the fit of analyzed distributions to determine which one was the best.The distribution with the largest p-value usually has the closest fit to the data.

Pollutants Indicators Modeling
Modeling of the pollutants indicators values in treated sewage was performed using Monte Carlo method.This method is defined as the using of the sequence of random numbers for creating the sample from the hypothetical population.Based on this population, there is a possibility of determining a statistical parameter estimator of the searched solution [40][41][42].Essentially, this is a creation of the stochastic model that describes real phenomena.Next, this model is running many times with using randomly generated variables (according to the assumed probability distribution).Finally, statistical analysis of the obtained results is performed.As some papers show [43], the Monte Carlo method is commonly used for the analysis of water and sewage management facilities operation work.It should be emphasized that the Monte Carlo method has some advantages and disadvantages.This method is very flexible and the empirical distributions of random variable can be handled.Moreover, the method can generally be easily extended and developed as required.Finally, this method is easily understood by non-mathematicians.Regarding the disadvantages, use of the Monte Carlo method usually requires specialist software.Next, the calculations can take much longer than analytical models.Moreover, the solutions are not exact, but depend on the number of repeated runs used to produce the output statistics.That is, all outputs are estimates.
It should be emphasized that Monte Carlo method has already been used for assessment of the wastewater treatment plant operation work.However, previous research has focused on the assumption of a specific form of statistical distribution for pollutants indicators (e.g., normal distribution), values of distribution parameters and the prediction of the pollutants indicators values based on the assumed information, such as at work [44].However, it should be stressed that some assumptions about the unique form of pollutants indicators distribution may be incorrect.In reality, there are many functions which may describe the empirical sample.Hence in this paper some novelties related to the Monte Carlo approach were included.Namely the specific form of distribution was not assumed.In performed analysis a few functions were considered.Next the best-fitting statistical distribution describing the empirical distributions of pollutants indicators was identified.Next for this function a simulation, using Monte Carlo method, was performed.One simulation series consist of 365 created random values of pollutants indicators.For each of the generated data set, 24-element's sample was generated in a random way.In each of the series, the number of simulation observations that not meet the technological efficiency condition was determined, i.e., if simulated values were greater than the permissible levels.Next, the result of the control of the wastewater treatment plant operation work was determined as follows; CF = 1, in case of the negative assessment result, i.e., if the number of samples that not meet the requirements exceeded the maximum amount of these type samples; CF = 0, in otherwise.Then, each simulation series was repeated 100 times.

Preliminary Analysis of Pollutants Indicators
At the beginning, the values of the descriptive statistics for tested pollutants indicators observational series were determined and presented in Table 1.In turn, Figure 2 shows, noted over the whole investigation period, pollutant indicators values in treated sewage with its their permissible levels.The results presented in Figure 2 and Table 1 indicate that the range between the minimum and maximum BOD 5 values was 23 mg O 2 •dm −3 .For the whole investigation period, 17% of the observations (eight samples) exceeded the limitation value for BOD 5 of 15 mg O 2 •dm −3 .In turn, the statistical range of COD Cr limit values was 67 mg O 2 •dm −3 .In the case of COD Cr , the exceedance of its permissible level of 125 mg O 2 •dm −3 was not observed.For total suspended solids (SS t ), a difference between the minimum and maximum concentration was 44 mg•dm −3 .Only one treated sewage sample did not meet the SS t quality requirements.The range of between the minimum and maximum total nitrogen concentration (TN) was 22.3 mg TN•dm −3 .Between 2015 and 2016, 19 samples of treated sewage (38% of all) exceeded the permissible level of total nitrogen concentration (10 mg TN•dm −3 ).In turn, the difference between the lowest and the highest noted total phosphorus (TP) concentration was 1.3 mg TP•dm −3 .Three samples of treated sewage were characterized by the exceedance of the limitation value for TP concentration of 1mg TP•dm −3 .The exceedance of the permissible values of pollutants indicators in treated sewage may be caused by different factors, such as sewage inflow irregularity or sewage pollutants content irregularity and sewage temperature in bioreactor.Both sewage quantity and sewage quality is strongly dependent on the city households' character.Inhabitants without an access to sewage systems and treatment plants are accustomed to save the water.This is the reason for greater pollutants concentration in smaller water volume.Moreover, the water usage trends in households cause noticeable sewage drainage heterogeneousness.It leads to sewage treatment systems operation work instability [43].Taking into consideration a coefficient of variation C s (Table 1), it can be stated that pollutants indicators values were characterized by a considerable variation.Such irregularity is typical for household sewage.Research results presented by [45][46][47][48][49] also indicate on the irregularity of pollutants concentration in household sewage.

The Results of Theoretical and Empirical Distributions Fitting
Theoretical distributions fitting to empirical pollutants indicator distributions was performed using the following functions; GEV, GMM, Log-normal, Normal, Pareto, Rayleigh, Triangle, and Weibull distribution.The assessment of the compatibility of distributions fitting was made using Anderson-Darling (A-D) test for the significance level of α = 0.05.Analysis results are presented in Table 2. Based on the results presented in Table 2, it can be stated that BOD 5 observational series can be described using all of the analyzed probability distributions.This is evidenced by the p-values for each of the calculated A-D parameters.Its values were greater than assumed significance level of α = 0.05.In the case of COD Cr , A-D test showed that for Pareto distribution, p-value is less than 5%.Therefore, in this case, H 0 hypothesis about the compatibility of theoretical distribution with empirical distribution was rejected.As the results indicate (Table 2), for SS t , there is no possibility for empirical distribution using Pareto and Triangular functions.Statistical analysis of TN showed the lack of the compatibility of theoretical and empirical distributions in the case of Pareto and Triangular distributions.Finally, H 0 hypothesis for TP concerning the distributions compatibility was rejected in the case of Triangular and Rayleigh functions.
The performed analysis showed that almost in every case, Triangular distribution cannot be used for the statistical describing of pollutants indicators observational series.Only in the case of BOD 5 , there is a possibility of using a Triangular distribution.However, it must be noted that the calculated p-value for BOD 5 of Triangular distribution (p-value = 0.054) was close to the H 0 hypothesis acceptability limit for Anderson-Darling test.The lack of the possibility of using a Triangular distribution is the result of the function characteristic.According to the Stein and Keblis [50] and Pereira et al. [51], Triangular distribution is a homogeneous function that is described by the most probable minimum and maximum values.Sewage treatment is a complex process.There are many time-variable factors affecting the sewage treatment processes that are the reason for the pollutants indicators variability in treated sewage.In consequence, prediction of the most repeated minimum and maximum values for Triangular distribution is difficult.In addition, it must be mentioned, that not much random variables types may be described using a Triangular distribution.

Selection of the Best-Fitting Statistical Distribution
The performed analysis showed that pollutants indicators observational series may be described by using more than one theoretical distribution.Thus, the selection of the best-fitting statistical distribution using PWRMSE parameter was conducted.Calculations were made only for these functions, where the A-D test confirmed functions compatibility with random variables empirical distribution.The results of the analysis are presented in Table 3.In turn, Figure 3 shows quantile-quantile (Q-Q) graphs for the best-fitting statistical distributions.The results presented in Table 3 indicate that both for BOD 5 and COD Cr , the best-fitting theoretical distribution is the GMM function, which is evidenced by the PWRMSE values equal to 0.802 mg O 2 •dm −3 and 1.511 mg O 2 •dm −3 for BOD 5 and COD Cr , respectively.In the case of SS t , the best-fitting theoretical distribution is GEV function (PWRMSE SSt = 1.222 mg•dm −3 ).In turn, for biogenic compounds, the best-fitting statistical distribution is GMM.This is evidenced by the PWRMSE values for total nitrogen (PWRMSE Nt = 0.545 mg TN•dm −3 ) and for total phosphorus (PWRMSE Pt = 0.028 mg TP•dm −3 ).The analysis allowed us to indicate, unequivocally, that for both organic pollutants indicators (BOD 5 and COD Cr ) and for biogenic compounds (Total Nitrogen and Total Phosphorus) the best-fitting theoretical distribution is the GMM function.Although in the case of SSt, the best-fitting function is GEV distribution, nevertheless, there is also a possibility of using GMM distribution.This is because of the similarity of PWRMSE parameter for total suspended solids GEV and GMM distributions.The possibility of using Gaussian Mixture Model (GMM) for the assessment of the operation work of the wastewater treatment plant was also presented by Yu [52] and Bouzenad and Ramdani [53].The results of the analysis performed using PWRMSE parameter turned out to be compatible with the results of Anderson-Darling (A-D) test.Both for GMM distribution (in case of BOD 5 and COD Cr ) and for GEV distribution (in case of SS t ), A-D test achieved the lowest values, while the p-value was the greatest among the tested distributions.One of the main A-D test disadvantages is the necessity for determination the critical values for each statistical distribution separately [54].In addition, the use of the A-D test requires using advanced mathematical tools.Therefore, PWRMSE method can be an effective alternative tool for A-D test, because of its simplicity of using.This is especially important in the case of using the statistical methods in practical activities, especially by the wastewater treatment plants operators.The results presented in Table 3 indicate that both for BOD5 and CODCr, the best-fitting theoretical distribution is the GMM function, which is evidenced by the PWRMSE values equal to 0.802 mg O2•dm −3 and 1.511 mg O2•dm −3 for BOD5 and CODCr, respectively.In the case of SSt, the bestfitting theoretical distribution is GEV function (PWRMSESSt = 1.222 mg•dm −3 ).In turn, for biogenic

Determination of the Wastewater Treatment Plant Operation Work Reliability Coefficients
Determination of the pollutants indicators values was made based on the results of the best-fitting distributions selection using the Monte Carlo method.This method is commonly used in generating draws from a probability distribution.The conducted calculations relied on repeated random sampling to obtain numerical results.The essential idea of simulation was using randomness to solve problem that might be deterministic in principle as in the case of assessment of operation work reliability for WTP.GMM distribution was used for organic and biogenic pollutants indicators modeling.In turn, for SS t , simulation was conducted based on the GEV distribution.Obtained simulation results were used for the calculation of technological efficiency coefficient R, reliability coefficient CR and the risk of the "negative" control of the wastewater treatment plant operation work Re (Table 4).Based on the obtained results (Table 4), it can be stated that wastewater treatment plant in Jaworzno works properly in terms of the organic pollutants reduction.The values of technological efficiency coefficient R are close to or equal to 1.0 for BOD 5 and COD Cr , respectively.These values indicate that the number of treated sewage samples that not fulfill the quality requirements may be noted occasionally.This is also evidenced by the reliability coefficient CR values.In turn, Re parameter for BOD 5 showed that the number of treated sewage samples with exceeded BOD 5 limitation values may be greater than the permissible samples amount.As the results indicate, in the case of COD Cr , such risk is not predicted (Table 4).Obtained results showed that sewage treatment processes on the WTP in Jaworzno provide an effective SS t reduction.The values of R and CR parameters showed, that exceedings of the permissible SS t concentration at the wastewater treatment plant outflow may be observed.However, the number of these exceedations is not greater than the permissible limit for over 100,000 p.e. wastewater treatment plant.This is evidenced by the calculated value of Re parameter (Table 4).Performed simulation showed, that in the case of biogenic pollutants indicators, the requirements concerning to the TN and total TP concentration at the wastewater treatment plant outflow may be not fulfill.The values of technological efficiency coefficient R and reliability coefficient CR confirmed that biogenic compounds concentration in treated sewage may exceed its limit content.In addition, the values of Re parameter for biogenic compounds, suggest that the number of TN and TP exceedings may be greater than its limitation.Especially in the case of TN, the value of the risk of the "negative" control of the wastewater treatment plant operation work (Re=1.0)allowed for the statement, that, in the analyzed period, WTP in Jaworzno did not provide the required reduction of total nitrogen concentration.Based on the obtained simulation results, it can be stated that the sewage treatment processes on the WTP in Jaworzno were partially improper.Therefore, wastewater treatment plant operators should take actions to increase the organic and biogenic pollutants reduction effectiveness.The results of the analysis were based on the simulation (prediction) of the wastewater treatment plant reliability for pollutants indicators values in 49 sewage samples.These pollutants indicators indicate on the current wastewater treatment plant reliability.However, using this data, it is possible to forecast the sewage treatment processes reliability in the future.It is possible to forecast the wastewater treatment plant operation work in the future with the current hydraulic load and current pollutants content in sewage.It must be emphasized that presented simulation method can be used for modeling the sewage quality parameters after the next stages of sewage treatment process.It is possible to make an analysis of sewage quality parameters before biological sewage treatment (mechanical treatment).Thank to this, wastewater treatment plant operators are able to make a prediction about the occurrence of increased pollutants indicators values in sewage after mechanical treatment processes.Increased pollutant indicator values in sewage after mechanical treatment may cause some problems related to the achieving the required pollutants reduction after biological treatment.The performed simulation can give some information about the necessity for performing in the future the modernization of wastewater treatment plant, including, first, mechanical facilities modernization, without modernization of biological facilities.Similar to the above, calculations may be performed only for biological facilities.Obtained results can indicate to the risk of the negative operation work of some wastewater treatment plant facilities.It will allow us to indicate which of the facilities should be modernized in the future.
The assessment of the water and sewage infrastructure (facilities) reliability may be also carried out using the other methods.In the paper [55], the graph theory method was presented as a means to identify the most critical elements in a network with respect to malfunctioning of the system as a whole.The method is objective and independent on the type of storm event and requires limited computational effort.In the paper [56], it was found that the 'bathtub' type failure rate curve can describe failures for particular objects of sewage infrastructure.Registration of failure data is very important aspect too.In the paper of Korving and Langeveld [57], the uniform registration of failures was performed.This method enables cooperation of the different management authorities on the level of daily operations.Additionally, the comparison of the performance can be made between different water and sewage facilities and changes in the performance of system components can be identified more easily.
As it was already mentioned, the irregularity of pollutants content in sewage inflowing to the bioreactor is one of the main factor affected the disruption of sewage treatment processes.Both in the case of the WTP in Jaworzno and any other similar objects located on rural areas in Poland, sewage coming from septic tanks represent a substantial part of treated sewage.Sewage from septic tanks is characterized by high pollutants concentration; in addition, it is rotten sewage.Because of this, it is necessary to determine the permissible share of the sewage coming from septic tanks in total amount of treated sewage [58][59][60].In order to increase the reduction of sewage biogenic pollutants, it is suggested to consider the application of biofilters before the bioreactor with activated sludge in wastewater treatment plant facilities.Biofilters will be used for nitrification process.In turn, total nitrogen and total phosphorus reduction process will take place in anaerobic and anoxic activated sludge chambers.Separation of the nitrobacteries from heterotrophic organisms will result with its no competition and greater stability of nitrification population [61].Research related to determination of influence to the reduction of biogenic pollutants was also conducted by W ąsik et al. [62].It should be noted that the decisions about the optimization of the water and sewage management, cannot be make only based on the empirical data analysis.Negotiations between the involved managers of the water and sewage facilities and other similar activities should be also considered [63].It is observed very often, that some decisions about the sewage infrastructure facilities is partially making based on the intuition, without any reliable confirmation or arguments [64].

Summary
The aim of the paper was the modeling of the wastewater treatment plant operation work using Monte Carlo method.The analysis includes some pollutants indicators in treated sewage: BOD 5 , COD Cr , SS t , TN, and TP.A simulation was performed using the following probability distributions: GEV, GMM, Log-normal, Normal, Pareto, Rayleigh, Triangular, and Weibull functions.A-D test was used for the assessment of theoretical and empirical distributions compatibility.In turn, the selection of the best-fitting statistical distributions was performed using PWRMSE parameter.Obtained results indicate on a significant variability of tested pollutants indicators in sewage after treatment processes on the wastewater treatment plant in Jaworzno.As the results of A-D test showed, Triangular distribution should not be used for statistical description of pollutants indicators observational series.Based on the results obtained using PWRMSE parameter, it was stated that the best-fitting random variables statistical distribution is GMM function.Because there was observed a similarity of PWRMSE method results and A-D test results, it was found that PWRMSE parameter can be an effective alternative for using A-D test.The results of Monte Carlo simulation method confirmed that some problems related to the organic (mainly BOD 5 ) and biogenic pollutants indicators reduction may appear on the WTP in Jaworzno.Finally, it was found that Monte Carlo simulation method is a useful tool for modeling the wastewater treatment plant operation work reliability.However, it is important to indicate an appropriate theoretical distribution for random variables.With this, it is possible to create a new data, maintaining simultaneously a current correlation structures between the random variables.New values obtained in this way may be helpful for the assessment of the wastewater treatment plant operational work and for preparing different scenarios of its operation, which has a significant content in many aspects of the assessment of the wastewater treatment plant operation work.Simulation method can be used by the wastewater treatment plant operators for making an assessment of the risk of the failures the wastewater treatment plant facilities that can occur at specific random incidents.Presented methods allow for consider the risk of the failures in a wider scope.This "wider scope" includes especially a prediction of the changes of the amounts of inflowing sewage and the changes of pollutants content in sewage.Thanks to this, it is possible to speed up or to dismiss over time the need for the wastewater treatment plant modernization or the expansion of some wastewater treatment plant facilities.It is very important economical aspect.Such activities allow planning the expenses intended for modernization or expansion the wastewater treatment plant.

Water 2019 , 17 Figure 1 .
Figure 1.A simplified scheme of the wastewater treatment plant in Jaworzno.
analysis was made based on the data including the values of the following pollutants indicators; BOD5, CODCr, SSt, TN, and TP.Data was sourced from the period of 2015 to 2016 and included measured values of pollutants indicators for 48 samples of treated sewage.Performed analysis consists of a few elements: the initial analysis of the pollutants indicators series, theoretical distribution fitting for pollutants indicators, selection the best-fitting theoretical function to the empirical distribution, modeling of the pollutants indicators values, and determination of the wastewater treatment plant reliability coefficients values based on the performed simulation results.

Figure 1 .
Figure 1.A simplified scheme of the wastewater treatment plant in Jaworzno.

BOD 5 ,
COD Cr , SS t , TN, and TP.Data was sourced from the period of 2015 to 2016 and included measured values of pollutants indicators for 48 samples of treated sewage.Performed analysis consists of a few elements: the initial analysis of the pollutants indicators series, theoretical distribution fitting for pollutants indicators, selection the best-fitting theoretical function to the empirical distribution, modeling of the pollutants indicators values, and determination of the wastewater treatment plant reliability coefficients values based on the performed simulation results.

−O
-mean of observed values; − P-mean of predicted values.

Table 1 .
The values of the descriptive statistics for pollutants indicators observational series.

Table 1 .
The values of the descriptive statistics for pollutants indicators observational series.

Table 2 .
The results of theoretical distributions fitting to empirical random variables distributions.

Table 3 .
The results of the peak-weighted root mean square (PWRMSE) analysis for tested pollutants indicators.

Table 4 .
The values of the wastewater treatment plant operation work reliability coefficients determined based on the Monte Carlo simulation method.