Next Article in Journal
Degradation of Sodium Acetate by Catalytic Ozonation Coupled with a Mn-Functionalized Fly Ash: Reaction Parameters and Mechanism
Next Article in Special Issue
Evaluating the Contamination by Indoor Dust in Dubai
Previous Article in Journal
Synthesis and Biodegradation Test of a New Polyether Polyurethane Foam Produced from PEG 400, L-Lysine Ethyl Ester Diisocyanate (L-LDI) and Bis-hydroxymethyl Furan (BHMF)
Previous Article in Special Issue
A New Method for Ecological Risk Assessment of Combined Contaminated Soil
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling the Chlorine Series from the Treatment Plant of Drinking Water in Constanta, Romania

by
Alina Bărbulescu
1 and
Lucica Barbeș
2,3,*
1
Department of Civil Engineering, Transilvania University of Brașov, 5 Turnului Str., 500152 Brasov, Romania
2
Department of Chemistry and Chemical Engineering, Ovidius University of Constanța, 124 Mamaia Bd., 900152 Constanta, Romania
3
Doctoral School of Biotechnical Systems Engineering, Politehnica University of Bucharest, 313, Splaiul Independentei, 060042 Bucharest, Romania
*
Author to whom correspondence should be addressed.
Toxics 2023, 11(8), 699; https://doi.org/10.3390/toxics11080699
Submission received: 1 July 2023 / Revised: 10 August 2023 / Accepted: 11 August 2023 / Published: 13 August 2023

Abstract

:
Ensuring good drinking water quality, which does not damage the population’s health, should be a priority of decision factors. Therefore, water treatment must be carried out to remove the contaminants. Chlorination is one of the most used treatment procedures. Modeling the free chlorine residual concentration series in the water distribution network provides the water supply managers with a tool for predicting residual chlorine concentration in the networks. With regard to this idea, this article proposes alternative models for the monthly free chlorine residual concentration series collected at the Palas Constanta Water Treatment Plant, in Romania, from January 2013 to December 2018. The forecasts based on the determined models are provided, and the best results are highlighted.

1. Introduction

Drinking water quality is essential, given its impact on the population’s health [1]. Therefore, ensuring a sufficient quantity and adequate quality must be a priority of each state/community to improve the health indicators and the population’s well-being [2]. The urban population’s primary drinking water supply sources are surface water and groundwater, whereas wells are used in rural areas [3]. In an ideal scenario, a water supply system would operate continuously, without changes in flow rate or other special conditions for individual treatment processes, when the raw water quality and quantity are constant. In reality, ideal conditions are not always met [4,5]. Given that various contaminants can affect the drinking water quality, it is crucial to treat the water before its distribution for consumption [6,7].
Due to its effectiveness (in killing viruses, bacteria, etc.), environmental feasibility, and long-lasting effects, chlorine is the primary disinfectant used for drinking water treatment [8,9]. Hypochlorous and hydrochloric acids are produced by adding chlorine or its derivatives to the raw water [10]. The active element in the disinfection process (the hypochlorite ion) results from the dissociation of the hypochlorous acid. During the water treatment, chlorine oxidizes the mineral substances and then produces chloramines by reacting with ammonia. Supplementing the chlorine dose leads to chloramine oxidation, increasing the free chlorine residual level [11,12], which is crucial for effective disinfection. The laboratory analyses performed on water samples taken at the outlet of the water treatment station and the distribution network indicate the disinfection stages and the necessary chlorine doses for ensuring water quality [13,14]. A balance in the chlorine dosing must be kept to protect the population against contamination, on the one hand, and avoid the by-products’ formation and pipes’ corrosion, on the other hand [14,15]. In these conditions, models that accurately predict the free chlorine residual in the distribution system have been proposed as a first step for optimizing the water treatment plant functioning.
Ghang et al. [16] introduced a chlorine decay model based on potential chlorine decay mechanisms and evaluated its performances on four raw surface and alum-treated waters. The results prove that the proposed model accurately predicts free chlorine residuals (R2 = 0.98). Gómez-Coronel et al. [17] reported satisfactory results in the chlorine concentration at the input of a water distribution system simulated in EPANET, with a genetic algorithm implemented in MATLAB. The EPANET MSX software was used to model chlorine decay in Algarve’s drinking water supply systems [18]. García-Ávila et al. [19] employed the same tool with a built-in first-order equation for modeling chlorine decay for a case study from Ecuador. Nejjari et al. [20] proposed a methodology for efficiently calibrating the free chlorine decay models tested on the Barcelona water transport network. Zhang et al. [21] elaborated a model for integrating water quality and operation for forecasting water production (using a genetic algorithm-enhanced artificial neural network). In contrast, other authors focused on optimizing the chlorine dosing [22,23].
Quantifying chlorine residual, turbidity, standard plate count (SPC), coliforms, etc., was performed using statistical methods in a water distribution system from Pakistan [24]. The correlations between the coliforms’ presence in the water and the free chlorine content in the Parisian distribution system were also analyzed based on statistics and econometrics approaches [25]. For Romania, only a few studies provide results on drinking water treatment [13,26,27].
To summarize, most results on the chlorine concentration series in water distribution systems use differential equations and a few other methods, such as artificial intelligence. Despite the last period, econometrics and hybrid methods proved their efficiency for modeling and forecast time series in different research fields, like economics [28,29,30], signal analysis [31], hydro-meteorology, environmental pollution [32,33,34,35], and pharmaceutics [36], they were less utilized in modeling the chlorine series at the outlet of the water treatment plants and in the water distribution systems.
In the above context, this article proposes alternative models (econometrics not based on differential equations) for the free chlorine residual concentrations series collected in the water treatment plant Palas (Constanta, Romania) from January 2013 to December 2018. It also emphasizes the possibility of using them for the forecast. The proposed approaches are univariate, not multivariate, as in most of the above-cited literature. They do not require deep specific knowledge in the modeling field (as in the case of differential equations and artificial intelligence) and are easily understood and utilized. Another advantage is extending the research to an area less explored in Romania, for which only a limited number of studies were performed. The models are compared, and their weaknesses and advantages are highlighted.

2. Materials and Methods

2.1. Data Series and Statistical Analysis

The Palas Constanţa treatment, storage, and pumping complex (PCTC) is located in the industrial area of Constanţa city on the Black Sea Littoral in Romania (Figure 1) and provides water to about 350,000 inhabitants.
The groundwater sources that feed the treatment plant are Cișmea I A, Cișmea I B, Cișmea I C, and Cișmea II. Cișmea I A + B+C are formed of 36 wells with depths from 50 to 120 m, except P35, with a depth of 300 m. They have a total supply capacity of 7657 m3/h. Cișmea II has 12 wells with depths between 90 and 150 m and a pumping capacity of 1940 m3/h. The Galeșu surface water source, with 13,050 m3/h catching capacity, is situated along the banks of Poarta Alba–Midia (on the Channel Danube–Black Sea). It has five intakes equipped with metal sieves for retaining the suspended particles.
This source was created to cope with the high water consumption during the summer and supplement Constanța city’s water supply when necessary. The water quality is good even before its treatment, according to [35,37]. After the treatment, the water must satisfy the Directives of the Council of the European Communities [38,39] and the Water Framework Directive [40]. The PCTC stores the water, which is distributed to Constanța and the Littoral water supply system. According to [41], in 2020, the total amount of water supplied to the inhabitants of Constanta was 42,150 m3 per day.
Generally, for the drinking water distribution networks, there is a risk of insufficient drinking water distributed to consumers caused by phenomena such as the clogging of water sources or the lowering of the surface water level due to drought and lack of precipitation [42,43]. To avoid such situations, there are four water storage stations in Constanța, each of 20,000 m3, one of 6.000 m3, and another of 10,000 m3. The Caragea Dermen groundwater source can also be accessed. It is formed by 18 wells with depths between 35 and 90 m and has a supply capacity of 3.549 m3/h. The water from different sources undergoes different chlorination processes. Only after chlorination are the streams of water mixed and introduced into the distribution network. The studied data series (Figure 2) is formed of the monthly free chlorine residual concentration collected at the outlet of PCTP during January 2013–December 2018.

2.2. Statistical Analysis

Basic statistics (mean, median, standard deviation—SD, variation coefficient—CV) were first computed for the monthly series. Then, the following hypotheses were tested: normality against the non-normality (by the Jarque-Bera [44], Shapiro–Wilk [45], and Anderson–Darling [46] tests), homoscedasticity against heteroskedasticity (by the Levene test) [47], the series stationarity vs. its nonstationarity in mean and variance (by the KPSS test) [48]. The null hypothesis that there is no time series trend was tested against the alternative that a monotonic trend exists via the Mann–Kendall and seasonal Mann–Kendall test [49,50,51]. When the null hypothesis is rejected, Sen’s procedure [52] can be used to determine the monotonic trend.

2.3. Mathematical Modeling

Since the preliminary statistical analysis revealed the series seasonality, different approaches have been adopted to model the data series.
In the first approach, the series ( y t ) was decomposed using an additive model, of which its components are the trend, the seasonal component, and the random variable. In this case, the steps were the following [53]:
  • Determine the trend using the linear trend computed via Sen’s method;
  • Calculate the detrended series by subtracting the trend from the data series;
  • Determine the seasonal component;
  • Determine the remainder (random or residual component) as the difference between the detrended series and the seasonal component.
In the multiplicative decomposition, the steps are similar, but the addition is replaced by multiplication and the subtraction by division in the second and fourth steps from the previous method.
In the second approach, the decomposition was conducted following a similar procedure, but the trend was determined using a moving average method of the 12th order.
The third approach was to use the Holt–Winters method, where the series was decomposed using Equations (1)–(4) in the additive model, with a seasonal period p = 12 as follows:
y ^ t + h = a t + h b t + s t 11 + h 1 m o d   12 ,
with
a t = α y t s t 12 + 1 α a t 1 + b t 1 ,
b t = β a t a t 1 + 1 β b t 1 ,
s t = γ y t a t + 1 γ s t 12 ,
In the multiplicative model, the equations are (5)–(8), which are expressed as follows:
y ^ t + h = ( a t + h b t ) s t 11 + h 1 m o d   12 ,
where
a t = α y t / s t 12 + 1 α a t 1 + b t 1 ,
b t = β a t a t 1 + 1 β b t 1 ,
s t = γ y t / a t + 1 γ s t 12 ,
in the hypothesis that a t and s t 12 are not zero.
In (1)–(8), α ,   β ,   γ are smoothing parameters that must be determined for the level, a t , trend, b t , and seasonal component, s t , respectively [54].
The fourth proposed model is a Seasonal Autoregressive Integrated Moving Average model, SARIMA. An ARIMA (p,d,q) process ( x t ) with a constant is defined by the following:
ϕ L 1 L d y t = c + θ ( L ) ε t ,
where L is the backward operator and
ϕ L y t = 1 i = 1 p ϕ i L i y t ,
θ L ε t = 1 + i = 1 q θ i L i ε t ,
where
y t = x t x t d .
p and q are the numbers of autoregressive and moving average terms, respectively, d is the differentiation degree, and ( ε t ) is white noise.
A SARIMA (p,d,q) × ( P , D , Q ) m (seasonal ARIMA model) is expressed as the following equation:
ϕ L Φ ( L m ) ( 1 L ) d ( 1 L m ) D y t = θ L Θ ( L m ) ε t ,
where
Φ L y t = 1 i = 1 P Φ i L i y t ,
Θ L ε t = 1 + i = 1 Q Θ i L i ε t ,
m, D, P, and Q represent the number of seasonal periods, the seasonal differencing, autoregressive, and seasonal moving average terms, respectively [55].
The residual independence was tested using the Box–Ljung test [56].
In all cases, apart from the residuals‘ analysis (normality, homoscedasticity, and randomness), the mean absolute deviation (MAD), mean standard deviation (MSD), and mean absolute percentage error (MAPE) were also computed to assess the models’ quality. Comparisons of the models, their advantages, and drawbacks are finally discussed.
The MINITAB 17, trial version (https://www.minitab.com/en-us/products/minitab/, accessed on 15 June 2023) and the R software, v.4.3.1 (https://www.r-project.org/, accessed on 15 June 2023) were utilized for testing the statistical hypotheses and mathematical modeling.

3. Results and Discussion

3.1. Results of the Statistical Analysis

The basic statistics of the data series are as follows: minimum = 0.200, maximum = 0.7400, mean = 0.4835, median =0.5000, standard deviation (SD) = 0.1181, coefficient of variance (CV%) = 24.42, skewness = −0.22, and kurtosis = −0.07. Thus, there is a small variation in the series values, and the distribution is left skewed.
Based on the above results, the computed value of the Jarque–Bera statistics was 0.4384, indicating that the normality hypothesis cannot be rejected at a significance level of 0.05. A similar result was obtained by applying the Shapiro–Wilk test. The p-value computed in the Levene test is 0.582 > 0.05, so the homoscedasticity hypothesis cannot be rejected. The statistics of the KPSS test for level (trend) stationarity is 0.59209 (0.03532), and the p-value is 0.02336 (0.1). So, the hypothesis of the level stationary is rejected, and that of the trend stationarity cannot be rejected at the significance level of 0.05. The Mann–Kendall test and its seasonal version rejected the null hypothesis. Therefore, based on Sen’s procedure, a linear trend, with the following Equation (16) can be fitted:
Y t = 0.001429 ( t 1 ) + 0.542143 ,
where Y t is the concentration in the month t.

3.2. Models

When using the first approach, the series decomposition via the additive model (denoted as DECA) is presented in Figure 3a. The recorded (Actual) and the computed (Fitted) values are represented in blue and brown, respectively, and the trend is in green. The violet curve represents the series forecast for the next 48 months. The residuals are normally distributed, according to the Q-Q plot (Figure 3b) and the results of the Shapiro–Wilk test. They are homoscedastic (the p-value of the Levene test is 0.582 > 0.05) and autocorrelated (the first-order correlation coefficient is −0.3195).
The highest seasonal index corresponds to November, and the lowest to June (Figure 4a). The highest variations of the detrended series (Figure 4b) are those from November and March and the lowest from October.
The highest percentage variations per season (Figure 4c) were in March and November. The highest variation in the residual component (therefore, the worst fitted value) was in March, and the lowest one was in October (Figure 4d).
A similar behavior is noticed in the case of the multiplicative decomposition model (denoted in the following as DECM). Figure 5 shows the original series, the detrended one, the seasonally adjusted series, and the residual one.
Removing the trend from the initial series increases the series range. The seasonally adjusted series presents a lower variance than the original one, indicating that seasonality is a significant component of the series. The multiplicative decomposition model with a linear trend is slightly worse than the additive one since the mean absolute deviation (MAD) of 0.0773, mean standard deviation (MSD) of 0.0114, and mean absolute percentage error (MAPE) of 18.642 are higher than those in the additive model (0.0767, 0.0098, and 18.4257, respectively). Still, the models do not provide significant differences between the seasonal components, percent variation per season, or residuals per season. The hypotheses of the residuals series normality and homoscedasticity could not be rejected, but the randomness could. Therefore, one should look for a model with uncorrelated residuals to avoid the errors’ propagation.
In the second approach (decomposition with a 12th-order moving average trend), the best model was the additive one (denoted as MAA12). Figure 6 shows the initial series (observed), its trend, the seasonal, and the random component (residual). Due to the moving average computation, the trend is not linear or monotonically decreasing.
The seasonal indices are, respectively, Jan = 0.04668, Feb = 0.03876, Mar = −0.01790, Apr = 0.00360, May = −0.01940, June = −0.08790, July = −0.05149, Aug = −0.078569, Sept = −0.05648, Oct = 0.10001, Nov = 0.08059, and Dec = 0.04210. In this case, the highest values of the seasonal component are recorded in October, followed by November, and the lowest in June. The highest seasonal values are correlated to the higher chlorination necessity (in November and December) after the high season and the precipitation absence in summer (to maintain the quality of the drinking water), respectively, to the lowest chlorination necessity in June after the spring season and the high precipitation period.
The random component’s analysis provides a p-value of 0.9195 in the Shapiro–Wilk test, so the normality hypothesis cannot be rejected. The correlogram (Figure 7) shows again a first-order autocorrelation of the random component’s values.
The hypothesis of the random component’s homoscedasticity could not be rejected. For all statistical tests, the significance level was kept at 0.05. In MAA12, which is better than the multiplicative model with a 12th-order moving average trend (denoted as MAM12), MAD = 0.07601, MSD = 0.00870, and MAPE = 18.6546. In terms of MSD and MAD, the MAA12 is the best, while with respect to MAPE, the best is DECA. In both situations, a first-order autocorrelation of the residual series is present, so a third approach, the Holt–Winters method, was proposed to describe the series evolution.
Figure 8a provides the series decomposition using the multiplicative Holt–Winters method (denoted as MHW). The smoothing parameters are α = 0.04697, β = 0.07233, and γ = 0.43818, and the initial parameters and seasonality indices are a = 0.38112 ,   b = 0.00191 ,   s 1 = 0.10025, s 2 = 0.03342, s 3 = 0.06512, s 4 = 0.03010, s 5 = 0.01835, s 6 = −0.0806, s 7 = 0.00810, s 8 = −0.09001, s 9 = −0.04201, s 10 = 0.11143, s 11 = 0.12542, and s 12 = 0.11005.
In MHW, the level decreases in time, the trend increases (but not monotonically), and the seasonal component is not constant, according to the regression Equation (4) (or (8) in the multiplicative model).
Adding up the values of the level with the corresponding ones of the trend will result in a decreasing series of values (a decreasing trend in the first approach). Similar results were obtained using the multiplicative Holt–Winters method.
The level compound’s shape in MHW is concordant with the time series non-stationarity in level. Among the seasonal components, the highest values are recorded in November and October, followed by December. The seasonal values of the chlorine introduced in water in the treatment station after the high season are higher than in other periods (do not forget that the treatment plant is situated on the Black Sea Littoral in a tourist area, and during summer, the pollution is higher than in the rest of the year) and depends as well on the precipitation record during summer (that can carry the pollutants affecting the source water quality). In the additive Holt–Winters model (denoted as AHW), MAD = 0.0803, MSD = 0.0130, and MAPE = 18.8673, whereas in MHW, the corresponding values are MAD = 0.0772, MSD = 0.0118, and MAPE = 18.2619.
The tests on residuals did not reject their normality (see the histogram in Figure 8b) and homoscedasticity, but the randomness (see the correlogram in Figure 8c).
Figure 9 illustrates the MHW model’s forecast for the next 48 months.
The series values are represented in blue, and the confidence intervals at 99% and 95% confidence levels are represented in two nuances of grey. The shape of the forecast curves is similar to that of the data series, confirming the modeling quality.
The advantage of this approach is that the level is considered, and the seasonal indices are updated at each step of the algorithm. The first two models incorporate the level and trend into a single component (trend), which does not reflect the series variation from the base.
The last model is of SARIMA(0,1,1)(0,1,1)12 type. For its validation, the residuals’ series analysis was performed. The Shapiro–Wilk test indicates that the hypothesis that the series in Gaussian cannot be rejected (p-value > 0.100 > 0.05; Figure 10a), the correlogram (Figure 10b) indicates the correlation absence, and the Levene test (Figure 10c) rejected the heteroskedasticity hypothesis. The p-value associated with the Box–Ljung test is p = 0.1137, indicating that the hypothesis of residuals’ series independence cannot be rejected. Moreover, MAD = 0.0695, MSD = 0.00868, and MAPE = 16.5426, showing that the SARIMA performs best among all the proposed models.
Figure 11 presents the series forecast based on the built SARIMA model in the blue curve and the confidence intervals at 99% and 95% confidence levels.
To emphasize the performances of the forecast obtained using the MHW and SARIMA, their output was compared with the series values in recorded 2019 (that were not used for modeling). Figure 12 shows that the predicted values obtained using SARIMA are closer to the recorded values via comparison to MHW. The worst forecast was obtained for July and the best one for December.
The goodness of fit indicators for SARIMA (MHW) are MAD = 0.01039 (0.02118), MSD = 0.00016 (0.00068), and MAPE = 2.8181 (5.5738), showing that the SARIMA model is better than MHW.
All the approaches gave good results in modeling the free residual chlorine series, but the best (and more complex one) is the SARIMA(0,1,1)(0,1,1)12.
As mentioned, the chlorine quantity decreases during the disinfection processes due to the reactions with different substances. Keeping its concentration within optimal limits can be done if this parameter is monitored over time. Traditionally, process-based models to forecast chlorine decay use generally first-order equations [18,23,57]. To build such models, advanced knowledge of the phenomena that appear in the pipes, and accurate and sufficient data on some water parameters in the distribution system are necessary (the last must be experimentally obtained). Often, the coefficients in such models depend on the loading conditions and are not practical for modeling purposes [57]. Therefore, other approaches are required [17,21].
The second approach involves utilizing data-driven statistical models; this means that the forecast of residual chlorine utilizes relationships between the response variable and some regressors. If the experimental data on some variables are difficult to obtain, imprecise, or unavailable, the data-driven models are excellent alternatives to the process-based models [58]. In such models, the knowledge of the processes from the system is less important [14]. Their main advantage is that a deep knowledge of the mathematics and chemistry laws governing chlorine behavior is not necessary [59]. Among the data-driven statistical methods, we mention the linear autoregressive models to predict chlorine concentration and its decay in distribution systems and storage [14,60,61]. This present study falls into this category. Based on our best knowledge, it proposed four models that were first employed for modeling free chlorine monthly series. Therefore, comparisons with the results of similar studies conducted on different series cannot be performed.

4. Conclusions

This article proposed four alternative approaches for modeling monthly free chlorine residual concentration series from PCTP using decomposition, Holt–Winters, and SARIMA models. The novelty of this approach is the use of univariate econometric models in engineering and extending the results of other studies on the water treatment plant in Romania (that previously presented only basic statistical analysis or models of chlorine decay).
In the first approach, the trend was built using a nonparametric Sen’s method, which has the advantage that no other restrictions are to be satisfied by the parameters of the linear trend. Another advantage of this method is its simplicity. The second method has the advantage that it can be applied even in a situation when the hypothesis that a monotonic trend exists is rejected. Nevertheless, the twelve values of the series cannot be estimated. In the Holt–Winters method, the seasonality factors and the trend are updated at each step, which gives a more realistic picture of the evolution of each component compared to the classical decomposition. While the first two approaches are simpler, the third one includes a fourth component, the level, as a base from which the series vary. The Holt–Winters model is in concordance with the stationary test results. The SARIMA(0,1,1)(0,1,1)12 model is more complex since it involves the first-order differentiation of the series and its seasonal components (to reach its stationarity), and considering the innovation process (by the presence of the moving average, one for both series and seasonality). While the last methodology provides the most accurate results, all the others may be used for modeling and forecast given the easiness and availability of their implementation in MINITAB and R.
In Romania, the studies in the above field are either experimental, present basic statistics of some water parameters series (without correlations to each other) or use the first-order chlorine decay model. Therefore, this article completes the very sparse research in the field. Since the chlorine concentration is regularly monitored, and exceeding the limits imposed by regulation may give birth to protests from the residents that acknowledge the smell and taste of the drinking water, the amount of chlorine must be dosed taking into account the input water quality, resulting from the analyses of chlorine concentrations and the necessity to conform the Romanian regulations.
Despite their performances, the models presented here should be used only for short-time prediction without updating the input given a decreasing trend from the level from which the series’ values vary. Updating the input of the models is recommended for improving the forecast. Automating the chlorine concentration monitoring will result in a better dosage and forecast.
Another note is that the models do not include the risk factors and the solution for the situation when the water quality decreases. Therefore, in a future study, these aspects will be considered because there is a need to constantly monitor the water resources and the quality in the water treatment process and to intervene to maintain it.

Author Contributions

Conceptualization, A.B. and L.B.; methodology, A.B. and L.B.; software, A.B.; validation, A.B.; formal analysis, L.B.; investigation, A.B. and L.B.; resources, A.B.; and data curation, A.B.; writing—original draft preparation, A.B. and L.B.; writing—review and editing, A.B. and L.B.; visualization, A.B.; supervision, A.B.; project administration, A.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data will be available upon request from the authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Brandt, M.J.; Johnson, K.M.; Elphinston, A.J.; Ratnayaka, D.D. (Eds.) The Demand for Potable Water. In Twort’s Water Supply, 7th ed.; Butterworth-Heinemann: Amsterdam, The Netherlands, 2017; Chapter 1; pp. 1–36. [Google Scholar]
  2. Brandt, M.J.; Johnson, K.M.; Elphinston, A.J.; Ratnayaka, D.D. (Eds.) Disinfection of Water. In Twort’s Water Supply, 7th ed.; Butterworth-Heinemann: Amsterdam, The Netherlands, 2017; Chapter 11; pp. 475–511. [Google Scholar]
  3. Moran, S. (Ed.) Clean water unit operation design: Biological processes. In An Applied Guide to Water and Effluent Treatment Plant Design; Butterworth-Heinemann: Amsterdam, The Netherlands, 2018; Chapter 9; pp. 111–116. [Google Scholar]
  4. Mosse, P.; Murray, B. Good Practice Guide to the Operation of Drinking Water Supply Systems for the Management of Microbial Risk; Final Report–WaterRA Project 1074; Water Research Australia Limited: Adelaide, Australia, 2015; p. 23. Available online: https://www.wsaa.asn.au/sites/default/files/publication/download/Good%20Practice%20Guide%20April%2015.pdf (accessed on 10 June 2023).
  5. Bărbulescu, A.; Barbeș, L. Statistical methods for assessing the water quality after the treatment on a Sequencing Batch Reactor. Sci. Total Environ. 2021, 752, 141991. [Google Scholar] [CrossRef] [PubMed]
  6. Ecological Risk Models and Tools. U.S. Environmental Protection Agency (EPA). 2021. Available online: https://www.epa.gov/risk/ecological-risk-models-and-tools (accessed on 10 June 2023).
  7. Alcayhuamán Guzmán, R.M.; Al-Emam, R.; Alhassan, H.; Ali, A.; Allély-Fermé, D.; Ampomah, B.; Anarna, M.S.S.; Bakir, H.; Bani-Khalaf, R.; Bartram, J.; et al. Water Safety Plan Manual-Step-by-Step Risk Management for Drinking Water Suppliers; World Health Organization: Geneva, Switzerland, 2009; Available online: www.who.int (accessed on 10 June 2023).
  8. WHO. Water Safety in Distribution Systems; World Health Organization Document Production Services: Geneva, Switzerland, 2014; Available online: www.who.int (accessed on 10 June 2023).
  9. Basic Information about Chloramines and Drinking Water Disinfection. Available online: https://www.epa.gov/dwreginfo/basic-information-about-chloramines-and-drinking-water-disinfection (accessed on 10 June 2023).
  10. Dubey, S.; Gusain, D.; Sharma, Y.C.; Bux, F. The occurrence of various types of disinfectant by-products (trihalomethanes, haloacetic acids, haloacetonitrile) in drinking water. In Disinfection By-products in Drinking Water; Priya, T., Mishra, B.K., Prasad, M.N.V., Eds.; Butterworth-Heinemann: Amsterdam, The Netherlands, 2020; Chapter 15; pp. 371–391. [Google Scholar]
  11. Free Chlorine Residual Definition. Available online: https://www.lawinsider.com/dictionary/free-chlorine-residual (accessed on 15 July 2023).
  12. Romanian Law no. 458/2002 about the Quality of the Drinking Water. Available online: http://legislatie.just.ro/Public/DetaliiDocument/37723 (accessed on 15 July 2023).
  13. Iordache, A.; Woinaroschy, A. Analysis of the efficiency of water treatment process with chlorine. Environ. Eng. Manag. J. 2020, 19, 1309–1313. [Google Scholar]
  14. Gibbs, M.S.; Morgan, N.; Maier, H.R.; Dandy, G.C.; Nixon, J.B.; Holmes, M. Investigation into the relationship between chlorine decay and water distribution parameters using data driven methods. Math. Comput. Modell. 2006, 44, 485–498. [Google Scholar] [CrossRef]
  15. Priya, T.; Mishra, B.K.; Prasad, M.N.V. (Eds.) Physico-chemical techniques for the removal of disinfection by-products precursors from water. In Disinfection By-Products in Drinking Water; Butterworth-Heinemann: Amsterdam, The Netherlands, 2020; Chapter 2; pp. 23–58. [Google Scholar]
  16. Gang, D.C.; Clevenger, T.E.; Banerji, S.K. Modeling Chlorine Decay in Surface Water. J. Environ. Inform. 2015, 1, 21–27. [Google Scholar] [CrossRef]
  17. Gómez-Coronel, L.; Delgado-Aguiñaga, J.A.; Santos-Ruiz, I.; Navarro-Díaz, A. Estimation of Chlorine Concentration in Water Distribution Systems Based on a Genetic Algorithm. Processes 2022, 11, 676. [Google Scholar] [CrossRef]
  18. Monteiro, L.; Figueiredo, D.; Dias, S.; Freitas, R.; Covas, D.; Menaia, J.; Coelho, S.T. Modeling of chlorine decay in drinking water supply systems using EPANET MSX. Procedia Eng. 2014, 70, 1192–1200. [Google Scholar]
  19. García-Ávila, F.; Avilés-Añazco, A.; Ordoñez-Jara, J.; Guanuchi-Quezada, C.; Flores del Pino, L.; Ramos-Fernández, L. Modeling of residual chlorine in a drinking water network in times of pandemic of the SARS-CoV-2 (COVID-19). Sustain. Environ. Res. 2021, 31, 12. [Google Scholar] [CrossRef]
  20. Nejjari, F.; Puig, V.; Pérez, R.; Quevedo, J.; Cugueró-Escofet, M.À.; Sanz, G.; Mirats, J. Chlorine Decay Model Calibration and Comparison: Application to a Real Water Network. Procedia Eng. 2014, 70, 1221–1230. [Google Scholar] [CrossRef] [Green Version]
  21. Zhang, Y.; Gao, X.; Smith, K.; Inial, G.; Liu, S.; Conil, L.B.; Pan, B. Integrating water quality and operation into prediction of water production in drinking water treatment plants by genetic algorithm enhanced artificial neural network. Water Res. 2019, 164, 114888. [Google Scholar] [CrossRef]
  22. Gámiz, J.; Grau, A.; Martínez, H.; Bolea, Y. Automated Chlorine Dosage in a Simulated Drinking Water Treatment Plant: A Real Case Study. Appl. Sci. 2020, 10, 4035. [Google Scholar] [CrossRef]
  23. Pérez, R.; Martínez Torrents, A.; Martínez, M.; Grau, S.; Vinardell, L.; Tomàs, R.; Martínez Lladó, X.; Jubany, I. Chlorine Concentration Modelling and Supervision in Water Distribution Systems. Sensors 2022, 22, 5578. [Google Scholar] [CrossRef] [PubMed]
  24. Farooq, S.; Hasmi, I.; Qazi, I.A.; Qaiser, A.; Rasheed, S. Monitoring of Coliforms and chlorine residual in water distribution network of Rawalpindi, Pakistan. Environ. Monit. Assess. 2008, 140, 339–347. [Google Scholar] [CrossRef] [PubMed]
  25. Cun, C.; Durand, M.; Leguyader, M.; Martin, J.; Vilagines, R. Statistical study of the relationship between free chlorine levels and bacteriological checks B on systems in the Paris area. Sci. Total Environ. 2002, 284, 49–59. [Google Scholar] [CrossRef]
  26. Paun, I.; Chiriac, F.L.; Iancu, V.I.; Pirvu, F.; Niculescu, M.; Vasilache, N. Disinfection by-products in drinking water distribution system of Bucharest City. Rom. J. Ecol. Environ. Chem. 2021, 3, 13–18. [Google Scholar] [CrossRef]
  27. Vîrlan, C.-M.; Toma, D.; Stătescu, F.; Marcoie, N.; Prăjanu, C.-C. Modeling the chlorine-conveying process within a drinking water distribution network. Environ. Eng. Manag. J. 2021, 20, 487–494. [Google Scholar]
  28. Aivaz, K.-A.; Florea Munteanu, I.; Stan, M.-I.; Chiriac, A. A Links Between Transport Noncompliance and Financial Uncertainty in Times of COVID-19 PandMultivariate Analysis on the emics and War. Sustainability 2022, 14, 10040. [Google Scholar] [CrossRef]
  29. Aivaz, K.-A. Correlations Between Infrastructure, Medical Staff and Financial Indicators of Companies Operating in the Field of Health and Social Care Services. The Case of Constanta County, Romania. In Under the Pressure of Digitalization: Challenges and Solutions at Organizational and Industrial Levels, 1st ed.; Edu, T., Schipor, G.-L., Vancea, D.P.C., Zaharia, R.M., Eds.; Filodiritto Publisher, Inforomatica SRL: Bologna, Italy, 2021; pp. 17–25. [Google Scholar]
  30. Vancea, D.P.C.; Aivaz, K.-A.; Duhnea, A. Political Uncertainty and Volatility on the Financial Markets-the Case of Romania. Transform. Bus. Econ. 2017, 16, 457–477. [Google Scholar]
  31. Bărbulescu, A.; Dumitriu, C.S. ARIMA and Wavelet-ARIMA models for the signal produced by ultrasound in diesel. In Proceedings of the 2021 25th International Conference on Systems, Theory, Control and Computing (ICSTCC), Iasi, Romania, 20–23 October 2021. [Google Scholar] [CrossRef]
  32. Bărbulescu, A.; Nazzal, Y.; Howari, F. Assessing the groundwater quality in the Liwa area, the United Arab Emirates. Water 2020, 12, 2816. [Google Scholar] [CrossRef]
  33. Bărbulescu, A.; Barbeș, L. Models for pollutants’ correlation in the Romanian littoral. Rom. Rep. Phys. 2014, 66, 1189–1199. [Google Scholar]
  34. Nazzal, Y.H.; Bărbulescu, A.; Howari, F.; Al-Taani, A.A.; Iqbal, J.; Xavier, C.M.; Sharma, M.; Dumitriu, C.S. Assessment of metals concentrations in soils of Abu Dhabi Emirate using pollution indices and multivariate statistics. Toxics 2021, 9, 95. [Google Scholar] [CrossRef]
  35. Bărbulescu, A.; Barbeș, L. Assessing the Danube River water quality of the Danube River (at Chiciu, Romania) by statistical methods. Environ. Earth Sci. 2020, 79, 122. [Google Scholar]
  36. Singh, I.; Juneja, P.; Kaur, B.; Kumar, P. Pharmaceutical Applications of Chemometric Techniques. Int. Scholarly Resear. Not. 2013, 13, 795178. [Google Scholar] [CrossRef] [Green Version]
  37. Frîncu, R.-M. Long-Term Trends in Water Quality Indices in the Lower Danube and Tributaries in Romania (1996–2017). Int. J. Environ. Res. Public Health 2021, 18, 1665. [Google Scholar] [CrossRef] [PubMed]
  38. Council Directive 80/778/EEC Relating to the Quality of Water Intended for Human Consumption. Available online: https://www.fao.org/faolex/results/details/en/c/LEX-FAOC037618/ (accessed on 15 July 2023).
  39. Council Directive 98/83/EC on the Quality of Water Intended for Human Consumption-Repealed. Available online: https://www.fao.org/faolex/results/details/en/c/LEX-FAOC018700 (accessed on 15 July 2023).
  40. Water Framework Directive. Directive 2000/60/EC of the European Parliament and of the Council. 2000. Available online: https://eur-lex.europa.eu/resource.html?uri=cellar:5c835afb-2ec6-4577-bdf8-756d3d694eeb.0004.02/DOC_1&format=PDF (accessed on 15 July 2023).
  41. County Report. Drinking Water Quality. Constanta. 2020. Available online: https://dspct.ro/wp-content/uploads/2021/08/RAPORT-JUDETEAN-APA-POTABILA-2020.pdf (accessed on 20 July 2023). (In Romanian).
  42. Xie, Y.; Zilberman, D. Theoretical implications of institutional, environmental, and technological changes for capacity choices of water projects. Water Resour. Econ. 2016, 13, 19–29. [Google Scholar] [CrossRef]
  43. Lanz, B.; Provins, A. The demand for tap water quality: Survey evidence on water hardness and aesthetic quality. Water Resour. Econ. 2016, 16, 52–63. [Google Scholar] [CrossRef] [Green Version]
  44. Gel, Y.R.; Gastwirth, J.L. A robust modification of the Jarque-Bera test of normality. Econ. Lett. 2008, 99, 30–32. [Google Scholar] [CrossRef]
  45. Shapiro, S.S.; Wilk, M.B. An analysis of variance test for normality (complete samples). Biometrika 1965, 52, 591–611. [Google Scholar] [CrossRef]
  46. Anderson, T.W.; Darling, D.A. A Test of Goodness-of-Fit. J. Am. Stat. Assoc. 1954, 49, 765–769. [Google Scholar] [CrossRef]
  47. Levene, H. Robust test for equality of variances. In Contributions to Probability and Statistics: Essays in Honor of Harold Hotelling; Olkin, I., Ed.; Stanford University Press: Stanford, CA, USA, 1960; pp. 278–292. [Google Scholar]
  48. Kwiatkowski, D.; Phillips, P.C.B.; Schmidt, P.; Shin, Y. Testing the null hypothesis of stationarity against the alternative of a unit root: How sure are we that economic time series have a unit root? J. Econ. 1992, 54, 159–178. [Google Scholar]
  49. Kendall, M.G. Rank Correlation Methods, 5th ed.; Oxford University Press: London, UK, 1990. [Google Scholar]
  50. Mann, H.B. Non-parametric tests against trend. Econometrica 1945, 13, 245–259. [Google Scholar] [CrossRef]
  51. Hipel, K.W.; McLeod, A.I. Time Series Modelling of Water Resources and Environmental Systems; Elsevier Science: New York, NY, USA, 1994. [Google Scholar]
  52. Sen, P.K. Estimates of the regression coefficient based on Kendall’s tau. J. Am. Stat. Assoc. 1968, 63, 1379–1389. [Google Scholar] [CrossRef]
  53. 3.2 Time Series Components. Available online: https://otexts.com/fpp3/components.html (accessed on 23 April 2023).
  54. 8.3 Methods with Seasonality. Available online: https://otexts.com/fpp3/holt-winters.html (accessed on 23 April 2023).
  55. SARIMA Models. Available online: https://real-statistics.com/time-series-analysis/seasonal-arima-sarima/sarima-models/ (accessed on 23 April 2023).
  56. Ljung, G.M.; Box, G.E.P. On a Measure of a Lack of Fit in Time Series Models. Biometrika 1978, 65, 297–303. [Google Scholar] [CrossRef]
  57. Jonkergouw, P.M.R.; Khu, S.-T.; Savic, D.A.; Zhong, D.; Hou, X.Q.; Zhao, H.-B. A Variable Rate Coefficient Chlorine Decay Model. Environ. Sci. Technol. 2009, 43, 408–414. [Google Scholar] [CrossRef] [PubMed]
  58. Rodriguez, M.J.; West, J.R.; Powell, J.; Serodes, J.B. Application of two approaches to model chlorine residuals in Severn Trent Water Ltd. (STW) distribution systems. Water Sci. Technol. 1997, 36, 317–324. [Google Scholar] [CrossRef]
  59. Serodes, J.B.; Rodriguez, M.J.; Ponton, A. Chlorcast ©: A methodology for developing decision-making tools for chlorine disinfection control. Environ. Modell. Softw. 2001, 16, 53–62. [Google Scholar] [CrossRef]
  60. Serodes, J.B.; Rodriguez, M.J. Predicting residual chlorine evolution in storage tanks within distribution systems: Application of a neural network approach. J. Water Supp. Resear. Techn.-Aqua 1996, 45, 57–66. [Google Scholar]
  61. Rodriguez, M.J.; Serodes, J.B. Assessing empirical linear and non-linear modelling of residual chlorine in urban drinking water systems. Environ. Modell. Softw. 1999, 14, 93–102. [Google Scholar] [CrossRef]
Figure 1. (a) Map of Romania; (b) the Palas Constanța water treatment complex (PCTC).
Figure 1. (a) Map of Romania; (b) the Palas Constanța water treatment complex (PCTC).
Toxics 11 00699 g001
Figure 2. The monthly series of free chlorine residuals from January 2013 to December 2018.
Figure 2. The monthly series of free chlorine residuals from January 2013 to December 2018.
Toxics 11 00699 g002
Figure 3. (a) Time series decomposition plot for the studied series. DECA; (b) the Q-Q plot of the random component. Mean is the average of the residual component’s values, StDev is the standard deviation of the residual component’s values, N is the number of the values, AD is the value of the Anderson–Darling statistics from the Anderson–Darling applied to the residual component, and P-value is the p-value computed in the Anderson–Darling test on the residual component.
Figure 3. (a) Time series decomposition plot for the studied series. DECA; (b) the Q-Q plot of the random component. Mean is the average of the residual component’s values, StDev is the standard deviation of the residual component’s values, N is the number of the values, AD is the value of the Anderson–Darling statistics from the Anderson–Darling applied to the residual component, and P-value is the p-value computed in the Anderson–Darling test on the residual component.
Toxics 11 00699 g003
Figure 4. DECA: (a) Seasonal indices (1 corresponds to January and 12 to December); (b) Detrended data by season; (c) Percent variation by season; (d) Residuals by season.
Figure 4. DECA: (a) Seasonal indices (1 corresponds to January and 12 to December); (b) Detrended data by season; (c) Percent variation by season; (d) Residuals by season.
Toxics 11 00699 g004
Figure 5. DECM: (a) Original Data, (b) Detrended series, (c) Seasonally adjusted series, (d) Seasonally adjusted and detrended series (residual component).
Figure 5. DECM: (a) Original Data, (b) Detrended series, (c) Seasonally adjusted series, (d) Seasonally adjusted and detrended series (residual component).
Toxics 11 00699 g005
Figure 6. MAA12: (a) The initial series. (Observed is the default name given to it by the software); (b) Trend; (c) Seasonal component; (d) Random component.
Figure 6. MAA12: (a) The initial series. (Observed is the default name given to it by the software); (b) Trend; (c) Seasonal component; (d) Random component.
Toxics 11 00699 g006
Figure 7. The random component’s correlogram in MAA12. The blue dotted line represents the limits of the confidence interval at a 95% confidence level.
Figure 7. The random component’s correlogram in MAA12. The blue dotted line represents the limits of the confidence interval at a 95% confidence level.
Toxics 11 00699 g007
Figure 8. (a) Holt–Winters multiplicative model; (b) Residuals’ histogram; (c) Residuals correlogram.
Figure 8. (a) Holt–Winters multiplicative model; (b) Residuals’ histogram; (c) Residuals correlogram.
Toxics 11 00699 g008
Figure 9. Forecast with the MHW model. The black curve is the series, the blue one is the forecast and the grey backgrounds are the confidence intervals at 95% and 99%, respectively.
Figure 9. Forecast with the MHW model. The black curve is the series, the blue one is the forecast and the grey backgrounds are the confidence intervals at 95% and 99%, respectively.
Toxics 11 00699 g009
Figure 10. SARIMA model. Residual series analysis (a) Results of the Shapiro–Wilk test; (b) Correlogram; (c) Results of the Levene test.
Figure 10. SARIMA model. Residual series analysis (a) Results of the Shapiro–Wilk test; (b) Correlogram; (c) Results of the Levene test.
Toxics 11 00699 g010
Figure 11. Forecast based on the SARIMA model. The black curve is the series, the blue one is the forecast and the grey backgrounds are the confidence intervals at 95% and 99%, respectively.
Figure 11. Forecast based on the SARIMA model. The black curve is the series, the blue one is the forecast and the grey backgrounds are the confidence intervals at 95% and 99%, respectively.
Toxics 11 00699 g011
Figure 12. The series recorded in 2019 and the forecast based on MHW and SARIMA models.
Figure 12. The series recorded in 2019 and the forecast based on MHW and SARIMA models.
Toxics 11 00699 g012
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bărbulescu, A.; Barbeș, L. Modeling the Chlorine Series from the Treatment Plant of Drinking Water in Constanta, Romania. Toxics 2023, 11, 699. https://doi.org/10.3390/toxics11080699

AMA Style

Bărbulescu A, Barbeș L. Modeling the Chlorine Series from the Treatment Plant of Drinking Water in Constanta, Romania. Toxics. 2023; 11(8):699. https://doi.org/10.3390/toxics11080699

Chicago/Turabian Style

Bărbulescu, Alina, and Lucica Barbeș. 2023. "Modeling the Chlorine Series from the Treatment Plant of Drinking Water in Constanta, Romania" Toxics 11, no. 8: 699. https://doi.org/10.3390/toxics11080699

APA Style

Bărbulescu, A., & Barbeș, L. (2023). Modeling the Chlorine Series from the Treatment Plant of Drinking Water in Constanta, Romania. Toxics, 11(8), 699. https://doi.org/10.3390/toxics11080699

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