Next Article in Journal
Hydrogeochemistry, Water Quality, and Health Risk Analysis of Phreatic Groundwater in the Urban Area of Yibin City, Southwestern China
Next Article in Special Issue
The Hydrological Impacts of Retrofitted Detention Ponds for Urban Managed Aquifer Recharge in the Cape Flats, South Africa
Previous Article in Journal
A Numerical Study on Impact of Coal Mining Activity and Mine Water Drainage on Flow and Transport Behavior in Groundwater
Previous Article in Special Issue
A Practical, Adaptive, and Scalable Real-Time Control Approach for Stormwater Storage Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Probabilistic Model for Predicting the Performance of a Stormwater Overflow Structure as Part of a Stormwater Treatment Plant

1
Faculty of Environmental Engineering, Geomatics and Renewable Energy, Kielce University of Technology, Al. Tysiąclecia Państwa Polskiego 7, 25-314 Kielce, Poland
2
Institute of Geography and Environmental Sciences, Jan Kochanowski University, ul. Uniwersytecka 7, 25-406 Kielce, Poland
*
Author to whom correspondence should be addressed.
Water 2024, 16(24), 3598; https://doi.org/10.3390/w16243598
Submission received: 25 October 2024 / Revised: 6 December 2024 / Accepted: 11 December 2024 / Published: 13 December 2024
(This article belongs to the Special Issue Urban Stormwater Control, Utilization, and Treatment)

Abstract

:
The purpose of this study was to attempt to develop a stochastic model that describes the operation of the stormwater overflow located in the stormwater sewerage system. The model built for this study makes it possible to simulate the annual volume of the stormwater discharge, the maximum volume of the overflow discharge in a precipitation event, and the share of the latter in the total amount of stormwater conveyed directly, without pre-treatment, to the receiver. The dependence obtained with the linear regression method was employed to identify the occurrence of stormwater discharge. The prediction of the synthetic annual rainfall series was made using the Monte Carlo method. This was performed based on the determined log-normal distribution, the parameters of which were specified using 13-year rainfall series. Additionally, simulation of the stormwater overflow operation was performed with the use of a calibrated hydrodynamic model of the catchment. The model was developed using the Storm Water Management Model (SWMM). The results of the hydrodynamic simulations of the volume and number of discharges were within the scope of the probabilistic solution, which confirms the applicative character of the method presented in this study, intended to assess the operation of stormwater overflow.

1. Introduction

Linear drainage systems collect stormwater from the area surface and convey it to the stormwater treatment plant (STP). The technological facilities of the treatment plant (settling tanks, separators, hydrocyclones, and others) reduce the wastewater pollutant load. In intense precipitation events, when the amount of stormwater flowing from the catchment increases rapidly, exceeding the capacity of stormwater treatment plants, excess wastewater is directly discharged into the river via the stormwater overflow or may be stored for future use. This can lead to the disturbance of the biological and chemical equilibrium in the aquatic environment [1,2,3,4], the washout of the banks and bottom of the watercourse [5,6,7,8], and an increase in the peak flow below the wastewater discharge section [9]. As a result, computational methods have been developed to assess the effect of both the quantity and quality of the stormwater introduced into the watercourse on the conditions prevailing in it [10,11,12,13,14]. In quality-based methods, the amount of dissolved oxygen and concentrations of heavy metals and ammonia in the receiving body are analyzed [15], or the saprobity index is determined [16]. In quantity-based methods, the assessment is based on the number of wastewater discharges, or on the maximum instantaneous flow. The review of the literature [13,17,18] indicates that many studies focus on predictions of the annual number of discharge events through the overflow, and the methods applied can be divided into empirical, deterministic, and stochastic methods [19]. As regards the first category, the methods rely on multiyear measurements of precipitation amounts, and on the multiplicity of the overflow operation. This makes it possible to determine regression dependences or appropriate nomograms [9,17,20]. Deterministic models include the Kuiper method [21,22], in which the calculations of the number of stormwater discharge events are based on the rainfall time series recorded over multiyear periods. However, because the empirical and deterministic methods mentioned above disregard the random character of precipitation events, the dynamics of the event over time, and the catchment retention, the description of the runoff is greatly simplified. Consequently, these methods provide only estimates, and the results may be marred by major errors [19].
For an assessment of the performance of sewerage systems, the hydrodynamic model can be used [9,23,24]. One of the most commonly used is the Storm Water Management Model (SWMM) software [25,26,27], which results from the fact that the source code, written in the C++ programming language, can be modified to fit individual needs. It should be emphasized that to assess the performance of the stormwater overflow, retention reservoir, or other facility in the sewerage system, it is necessary to perform continuous simulations [28,29,30] using rainfall time series collected for many years [31]. The computations [9,32,33] indicate that simulation results can provide the basis to specify the annual number and volume of stormwater discharges through overflows in an assumed return period. However, because of the high costs of taking measurements of precipitation and flows, the development of the model is not always a viable option. Additionally, the complexity of the runoff phenomenon makes it necessary to introduce a number of empirical parameters to describe the catchment and precipitation characteristics, which can only be established by calibration. However, a large number of these and the interactions between the individual parameters [34,35,36] can cause problems related to the unambiguous determination of the numerical values of these parameters, which is a condition for obtaining a correct solution.
In view of the above remarks, stochastic models have been developed to assess the performance of stormwater overflow structures. These are based, among others, on the First-Order Reliability Method (FORM), the authors of which [37] accounted for the uncertainty of the estimated parameters in the hydrodynamic model. For the method to be used, it is necessary to employ complex numerical algorithms, which may limit the application of the method in engineering practice. In his hydrodynamic model, which also accounts for parameter uncertainty, Thorndahl [38] used the results of simulations performed on the basis of rainfall time series collected over a 50-year period. In their models, neither Thorndahl [38] nor other researchers [9,29] explained the dynamics of the variation in the discharge volume over the yearly cycle, which is a valuable piece of information with respect to the hydraulic conditions in the rainwater receiver.
The primary objective of this study was to develop a stochastic model to evaluate the performance of a stormwater overflow. The model makes it possible to determine the annual quantity and volume of stormwater discharges, taking into account the stochastic nature of rainfall events. The method developed by the authors also makes it possible to identify extreme events in which the maximum volume of wastewater discharge is found in the precipitation event, and also to compute the share of that volume in the total yearly discharge. For the identification of the event and the computation of the volume of the overflow discharge, a regression dependence was used. It was obtained on the basis of continuous simulations performed by the SWMM program. To predict yearly precipitation data, the Monte Carlo method was applied. Estimation of the statistical distribution parameters describing the dynamics of the rainfall profile in individual precipitation events was performed on the basis of observation data collected in the years 2008–2020.

2. Materials and Methods

2.1. Study Area

This study presents an analysis of the catchment of a Si9 sewer, with a total area of 62 ha, located in the central–eastern part of the city of Kielce in Poland (Figure 1). The main sewer, 1569 m in length, and with diameter ranging from 600 to 1250 mm, collects precipitation wastewater from 17 laterals. The total length of the system is 5583 m. The sewerage system is buried 1.59–3.55 m in the ground. The highest catchment point is located at 271.20 m above the mean sea level (amsl), and the lowest-lying point is at 260.00 m amsl. The building development of the area includes housing estates, institutional buildings, trunk roads, and side streets. In total, roofs constitute 14.3%, pavements 8.4%, roads 17.7%, car parks 11.2%, greenery 47.2%, and school playgrounds 1.3% of the catchment area [39].
Wastewater transported from the catchment by the main sewer, Si9, is delivered to the stormwater treatment plant via the separation chamber (DC in Figure 1). The route along which the wastewater flows through the treatment facility includes an oblong settling tank (SET), coalescence separator (SEP), and control manhole (CS) (Figure 1). When the water level in the chamber (DC) does not exceed 0.42 m, stormwater flows through four ϕ400 mm conduits into the settling tank, and when the level is higher, stormwater is delivered, via the stormwater overflow (OV) (Figure 2), into the discharge (outlet) channel, from where it travels to the receiving water body (i.e., the Silnica River).
Two ultrasonic flow meters (ISCO AV 2150, Teledyne ISCO, Lincoln, NE, USA) were installed in the Si9 collector. The first is located 3.0 m from the inlet to the separation chamber, and the other is at a distance of 5.0 m below the storm overflow. Devices indirectly measure the flow rate based on the direct measurement of the average velocity and the water level (hydrostatic pressure measurement). The water level measurement ranges from 0.01 to 3.0 m, with an accuracy of ± 0.003 m. Furthermore, a probe for measuring the filling, which triggers the automatic rainwater sampling device, was found in the separation chamber (sampler type 6712, Teledyne ISCO, Lincoln, NE, USA). The tipping bucket-type rain gauge RG50 from SEBA Hydrometrie GmbH (Kaufbeuren, Germany) was used to measure the depth of the rain. The measurement frequency ranged from 2 to 5 min, with a resolution of 0.1 mm.

2.2. Precipitation Event and Antecedent Period

The literature review [38,40,41,42,43] indicates that the rainfall definition is unclear because, for an event to be considered relevant, the rainfall depth should range 1–10 mm. As the present study was intended to predict the number of wastewater discharges through the overflow within a one-year cycle, the assumed minimal rainfall depth selected as a basis to show the statistical distribution should make it possible to identify cases when such discharge does not occur. Therefore, the rainfall depth cited in Arbeitsblatt DWA-A 118 [44], equal to Ptot = 10 mm, which defines an intense event, seems to be too high. Conversely, the minimum value of Ptot = 1 mm is not sufficient because, as shown in the investigations conducted by Kotowski et al. [45], it could be marred by a considerable measurement error related to the equipment used. In relation to the above, in this study, a precipitation event was assumed to mean a rainfall depth of a minimum of 2 mm. In many studies [40,41,42,43], this value provides a basis for determining multidimensional distributions of extreme values used to evaluate the operation of the sewerage system. The time interval that separates individual precipitation events was assumed to be 4 h [44].

2.3. Precipitation Data

Precipitation data were obtained from precipitation stations located within the study catchment area and about 2.0 km from its borders. These posts have been continuously collecting atmospheric data since 2008, and precipitation is continuously monitored with a resolution of 1–5 min. Data from 2008–2020 were used for the analysis, with the total precipitation (Ptot) ranging from 2.0 to 45.2 mm during the period. The highest observed daily precipitation totals were in 2010, when days with more than 20 mm of precipitation occurred eight times, including two days during which 43.8 mm and 45.2 mm fell, which, using the formula given by Suligowski [46], classifies them as precipitation with a probability of occurrence close to p = 5% (once in 20 years). The number of precipitation events analyzed during a given year ranged from 52 to 75. Finally, 480 rainfall events were used to determine the synthetic rainfall generator for the Pt10, Pt15 and Pt30 forecasts. The duration of rainfall of the observed events (td) varied in the range of 20–2366 min, and the rain-free period (tap) was 0.17–60 days. The maximum 10, 15, and 30 min values (Pt10, Pt15, Pt30, respectively) were adopted, according to the literature [47,48], as authoritative for describing the intensity of the course of the considered phenomenon in a rainfall event. In the analyzed catchment, they are Pt10 = 0.2–14.4 mm, Pt15 = 0.2–19.3 mm, and Pt30 = 0.3–21.2 mm.
For the observed precipitation events described by means of tap, Ptot, Pt10, Pt15, and Pt30, empirical statistical distributions were stated, to which one of the following distributions, namely, Weibull, chi-square, exponential, GEV, Fisher–Tippett, gamma, log-normal, or beta, was fitted. The congruence between the empirical and theoretical distributions was made based on the results of the Kolmogorov–Smirnov and chi-square tests.

2.4. Hydrodynamic Model

The hydrodynamic model of the catchment was developed using SWMM 5.0 software. The model consisted of 92 subcatchments, with an area of 0.12–2.10 ha, 200 manholes, and 72 sections of the sewers. Stormwater from the roofs of buildings is discharged directly into the sewer system and surface runoff is discharged into street inlets. The Manning roughness coefficients for impervious and pervious areas of the catchment were 0.015 and 0.15 m−1/3∙s, respectively, and the roughness coefficient of sewer walls was 0.015 m−1/3∙s. The amounts of ground retention of the impervious areas (52.8% of the catchment) and pervious ones (47.2%) were 1.5 and 6.0 mm, respectively, and the weighted average retention (davg,ret) determined based on the share of each type of development was equal to davg,ret = 0.528∙1.5 + 0.472∙6.0 = 3.62 mm.
The mathematical model of the STP was developed based on field measurements and the design documentation of the facility [49]. Details of the construction and calibration of the model can be found in the works of Szeląg and Bąk [50]. To compute the amount of wastewater flowing into the treatment plant and that of the wastewater discharged through the stormwater overflow (Vd), simulations performed with the SWMM 5.0 program [51,52,53] were used. To calibrate the developed model of the separation chamber with the stormwater overflow, flow measurement sequences from 2008 to 2020 were used.

2.5. Calibration of the Hydrodynamic Model

The following parameters were used to assess the consistency of the measured and simulated fillings of the separation chamber:
-
The ratio of the maximum measured to simulated filling:
δ = h K R ( m ) h K R ( s i m )
where hKR(m) is the maximum measured filling of the separation chamber (m), and hKR(sim) is the maximum calculated filling of the separation chamber in the SWMM program (m);
-
The correlation coefficient (R):
R = i = 1 n ( h o i h s r , o ) · ( h s i h s r , s ) i = 1 n ( h o i h s r , o ) 2 · i = 1 n ( h s i h s r , s ) 2
where ho(i) is the rainwater distribution chamber fill values observed with a one-minute time step (m); hs(i) is the distribution of the chamber fill values obtained from simulations with a one-minute time step (m); hsr,o is the average measured distribution chamber fill (m); hsr,s is the average distribution chamber fill calculated by numerical simulations (m); n is the total number of observations.

2.6. Analysis of the Effect of the Precipitation Characteristics on the Overflow Performance

Due to the fact that hydrodynamic simulations are time-consuming, and it is not always possible to calibrate the model in the optimal way, for practical purposes, only those measurement results are used on the basis of which the regression dependences are derived. The later application of these dependences is much simpler and faster than that of the hydrodynamic model. The results obtained using these dependences are close to those produced through numerical simulations. To determine the discharge volume or the size of the pollutant load, multiple regression or the least-squares method is usually used [38,54,55].
To determine the effect of the individual parameters describing the non-uniformity of the rainfall depth during the event (Ptot, Pt10, Pt15, Pt30), the event duration (td), and time interval separating individual events (tap) on the volume of the overflow discharge (Vd), a Pearson correlation matrix was generated. Also, the relations that were statistically relevant at the adopted confidence level of p = 0.05 were found. Using computation results to avoid the necessity of applying the hydrodynamic model to determine the number and volume of discharges, the multiple regression equation was employed:
V d = i = 1 j α i · P i + β · t d + χ · t a p θ
where the characteristics of precipitation events during which wastewater discharge through overflow occurred were identified on the basis of the dependence:
i = 1 j α i · P i + β · t d + χ · t a p = θ
where αi, β, χ, and θ are the empirical parameters found with the least-squares method, and Pi, tap, and td are the parameters that describe the precipitation.
Equation (4) describes the multidimensional space of precipitation events during which stormwater discharge through the overflow will not occur. After appropriate transformations, dependence (4) offers the possibility of determining the limit value of the rainfall depth, the exceeding of which leads to discharge.

2.7. Prediction of the Annual Number and Volume of Discharge Events

Because of the relatively short period of precipitation observations, it was assumed in the computations that the number of rainy days was the same as that in individual years, which is a fairly common assumption made by researchers [56]. Computational analyses performed by Szeląg et al. [47], which took into account a uniform distribution of the number of rainfall days per year, showed that the average numbers of overflow discharges obtained in the variant taking this distribution into account and in the variant excluding it are practically identical. In the probabilistic approach, to determine an annual series (n) of precipitation events, described by means of the td, td(t−1), tap, Pi, and Pi(t−1), it is necessary to run a simulation of a single series at least 1000 times using the Monte Carlo method [56,57].
Using the multiple regression model developed, the annual number and volume of discharges through the stormwater overflow were obtained for the generated precipitation series. On this basis, it was possible to find empirical distribution functions that describe the probability that, in the assumed time period, the discharge number and volume will not exceed a certain value. Additionally, the occurrence of events with the maximum stormwater discharge through the overflow was analyzed. Also, the share of such events in the overall annual discharge of wastewater into the receiving water body was defined.

3. Results and Discussion

On the basis of the computations performed using SWMM software and the results of the measurements of the total rainfall depth (Ptot), rainfall duration (td), and stormwater flow rate, the stormwater overflow parameters were stated by means of calibration. This made it possible to compute the criteria of the fitting of the water levels in the separation chamber for seven selected precipitation events (Table 1). Events 1–5 involved intense rainfall, during which overflow discharge occurred. Conversely, during events 6–7, such discharge was not noted.
The analysis of the results of the fitting of the stormwater overflow model (δ) shows that the maximum instantaneous flows discharged through the overflow, determined with the SWMM model, are underestimated by, maximally, 12%. As part of the calculations, 20 rainfall events were selected from each year, for which calibration was performed (10 each with and without discharge). Table 1 lists seven sample precipitation events for which the fill of the separation chamber was determined and compared with the results of measurements determining the δ and R, varying in the ranges of 0.88–0.97 and 0.86–0.94, respectively.
The values of the correlation coefficients indicate the satisfactory fitting of the measurement results and hydrodynamic simulations with respect to variation in the water level of the separation chamber. In view of the above, it can be stated that the hydraulic model of the separation chamber with the overflow, developed by the authors, describes the chamber operation with satisfactory accuracy and can provide a basis for further analyses. Figure 3 shows an exemplary comparison of the filling level of the separation chamber, measured and obtained from simulations for precipitation event No. 4 (Table 1). The patterns of both curves show large similarities.
The shift in the theoretical curve relative to the curve obtained from measurements may be due to the lower average value of the catchment retention for real conditions than what was assumed in the simulation, which results in a later response of the drainage system to surface runoff. In the model, rainfall events were separated out when performing the simulation, ignoring the current state of the catchment retention between the rainfall events, which affects the results obtained. In addition, there are many factors, such as the amount of sediment in the sewers, the condition of the pipes, and sections with contraflows, that are difficult to determine precisely without very thorough field surveys.
The computations showed that the coefficient of the overflow is 0.36, and the values of the coefficients of the local resistances at the inlets and outlets of the channels conveying stormwater to the settling tank were in the ranges of 0.48–0.53 and 0.95–1.10, respectively. Simulations carried out to model the stormwater overflow performance made it possible to establish a likely number and volume of the overflow discharges in individual years (2008–2020). The values are as follows: 20 (discharge events) and 9234 m3 in 2008; 17 and 7911 m3 in 2009; 15 and 6274 m3 in 2010; 18 and 14,219 m3 in 2011; 21 and 13,259 m3 in 2012; 22 and 13,749 m3 in 2013; 30 and 29,553 m3 in 2014; 26 and 23,456 m3 in 2015; 24 and 21,432 m3 in 2016; 28 and 32,126 m3 in 2017; 19 and 18,754 m3 in 2018; 17 and 11,652 m3 in 2019; 14 and 9873 m3 in 2020, which gives 271 discharges altogether over 13 years.
Based on computations of the discharge volumes and measured rainfall parameters (Ptot, Pt10, Pt15, Pt30, td, tap), the correlation matrix was determined (Table 2). The matrix provided the basis for deriving, with the use of multiple regression, the dependence for the determination of the volume of the overflow discharge.
In Table 2, a very strong correlation (r = 0.92) can be observed between the volume of the overflow discharge (Vd) and the maximum 30 min rainfall depth (Pt30), and a strong correlation (r = 0.85 and r = 0.83) holds for the 10 and 15 min rainfall depths. A medium-strength correlation (r = 0.62) exists between the variables Vd − Ptot, whereas for the pair of variables Vd − tap, a weak correlation is shown.
Taking into account the results presented in Table 2, and also the strong correlations between the parameters considered, to determine the regression dependence, the method of multiple stepwise regression backwards was used. On the basis of the computations, the following dependence was obtained:
V d = 243.02 ± 18.86 · P t 30 1.85           ( r = 0.920 )
where the value ± 18.86 is the standard deviation.
The dependence (5) describes the amount of wastewater discharged through the overflow during a single precipitation event, which is confirmed by the correlation coefficient r = 0.920. It follows from Equation (5) that stormwater overflow discharge will occur when Pt30 exceeds the rainfall amount of 1.85 mm, whereas the weighted average value of the catchment retention (davg,ret) is much higher and amounts to 3.62 mm.
For predictions of the annual number and volume of overflow discharges, the statistical distribution of the maximum 30 min rainfall depth in the precipitation event was determined. For this purpose, values of the probability (p) of the Kolomogorov–Smirnov and chi-square tests were determined for the following distributions: Weibull, chi-square, exponential, GEV, Fisher–Tippett, log-normal, and beta. The analysis of the results of the computations (Table 3, Figure 4) indicates that for the measurement data from the years 2008–2020, the log-normal distribution is the best fitted one (Kolmogorov–Smirnov test: p = 0.3130; and chi-square: p = 0.2526 for α = 5%), having a standard deviation of σ = 0.970 mm and a mean value of μ = 0.403 mm.
Simulations conducted using the Monte Carlo method generated, 1000 times, 58 precipitation events in a year. The results of the computations of the number of overflow discharges, based on the generated series, were illustrated in the form of a cumulative probability density distribution (Figure 5). Next, from Formula (5), discharge volumes for 30 min synthetic rainfall depths were determined. This allowed for the determination of the annual volume (Figure 6) and the identification of extreme events that occurred in the year and during which the maximum discharge through the overflow occurred (Figure 7). Additionally, on the basis of the simulation results, the share of the stormwater volume in a single precipitation event during which the maximum wastewater discharge took place in the total yearly discharge was determined (Figure 8).
The graph in Figure 5 shows that the expected value (p = 0.50) of the annual number of discharges through the stormwater overflow in the urbanized catchment of concern is 22, whereas the percentile values p = 0.05 and p = 0.95 equal 17 and 28, respectively. The yearly number of the discharges obtained with the SWMM program ranges from 15 to 29, which corresponds to percentile values of p = 0.01 and p = 0.97.
The yearly discharge volume required for the percentile of p = 0.50 is 12450 m3, and for p = 0.05 and p = 0.95, the values are 7303 m3 and 19,903 m3, respectively (Figure 6). According to the simulations carried out using the hydrodynamic model, the lowest annual discharge volume is 7911 m3, which corresponds to p = 0.07, whereas the highest is 29,953 m3 (for p = 0.99).
On the basis of the computations performed using the stochastic model shown in Figure 7, it can be stated that the expected value of the maximum volume of the overflow discharge in a single precipitation event equals 2341 m3, and values of the percentiles p = 0.05 and p = 0.95 correspond to discharge volumes of 1184 m3 and 5840 m3.
The comparison of the results of the computations obtained using the hydrodynamic model with the stochastic solution shown in the graph (Figure 7) indicates that the maximum volume of the discharge in a single precipitation event based on the SWMM simulations amounts to 6120 m3. This value is close to that of the determined percentile p = 0.95 in the probabilistic model (empirical distribution function). In contrast, the minimal volume of discharge in the precipitation event, computed with the mathematical model of the catchment, is 2120 m3 and corresponds to p = 0.41.
The developed probabilistic model offers the possibility of determining the expected value that describes the share of the maximum volume of the stormwater discharge in a single precipitation event in the total annual discharge (Figure 8). This share amounts to 0.18 for the percentile p = 0.50, whereas the extreme values, i.e., for p = 0.05 and p = 0.95, are equal to 0.12 and 0.36, respectively.
The shares of the maximum discharges in the overall annual discharge, computed using the SWMM software, correspond to a percentile range of p = 0.07–0.97. This indicates considerable variation in the individual discharge values in the years 2008–2020 relative to the total annual discharge.
The results of the computations of the number and volume of overflow discharges obtained using the calibrated hydrodynamic model in the SWMM program and the stochastic model show high variation in the parameters describing the performance of the stormwater overflow. The annual number and volume of discharges obtained, and also the discharge volume in a single precipitation event obtained with the SWMM model, remain in the range of the probabilistic solution. This confirms the application capacity of the method presented in this study. In addition, a wide range of variation in the share of maximum discharge volumes in a single event in the total yearly discharge indicates substantial differences as regards the wastewater volumes in individual precipitation events, in which the wastewater is directly conveyed to the receiving water body.
To order to expand the applicability of the developed model, it is recommended to extend the study to other urban catchments with different land use patterns and different topologies of stormwater sewer systems. This approach can be implemented in parallel with computational experiments performed for the study catchment. Using experiment planning techniques, it is possible to design virtual catchments with different characteristics, for which simulations of the overflow performance can be performed using the developed SWMM model. The statistical models obtained in this way can be compared with predictions for other urban catchments.

4. Conclusions

On the basis of the computations, it can be stated that the stochastic model presented to describe the performance of the stormwater overflow can be applied to determine the number and volume of overflow discharges, and also to identify the extreme events during which the maximum wastewater discharge occurred during the year. In the analyses presented above, the results of the numerical computations performed with the SWMM program allowed for establishing the regression dependence between the Pt30 and the volume of the overflow discharge. Additionally, as regards research paths and the number of parameters describing the catchment and the sewerage system characteristics in the SWMM model, and to account for the interactions between these parameters, it is advisable to conduct an analysis of the uncertainties of the estimated parameters in the hydrodynamic model. This makes it possible to identify reliable values of the parameters by means of suitable statistical methods, which allows for the physical interpretation of the parameters obtained in Equation (5).
The results of the hydrodynamic simulations of the volume and number of discharges were within the scope of the probabilistic solution, which confirms the applicative character of the method presented in this study, intended to evaluate the performance of stormwater overflows.
It is necessary to validate the parameters of the statistical distribution describing the variation in the maximum Pt30, and to account for the varied number of precipitation events in the yearly cycle by developing a statistical distribution to describe this parameter.
The method described in this paper can be applied to other small urban catchments for stormwater drainage networks with simple topologies, with the prerequisite of continuous high-resolution rainfall measurement series, as well as the identification of a regression relationship between the overflow discharge volume and the parameters describing the rainfall dynamics during individual rainfall events.
The adopted calibration method ignores interactions between the identified SWMM parameters. Further research directions should be focused on developing models for the Vd, td, and number of overflow discharge predictions that take into account the uncertainty of the SWMM parameters describing the catchment retention.

Author Contributions

Conceptualization, B.S. and J.G.; methodology, B.S. and Ł.B.; investigation, B.S., J.G., Ł.B. and A.Ś.; resources, A.Ś. and J.G.; writing—original draft preparation, B.S. and J.G.; writing—review and editing, Ł.B.; supervision, B.S. and Ł.B.; project administration, J.G. and Ł.B.; funding acquisition, J.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Science Centre of Poland under the Miniatura 3, grant number 2019/03/X/ST8/01446, entitled “Studies on the dynamics of stormwater quality changes during rainfall and snowmelt events in drainage systems in urbanized areas”.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of the data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Yuan, Q.; Guerra, H.B.; Kim, Y. An investigation of the relationships between rainfall conditions and pollutant wash-off from the paved road. Water 2017, 9, 232. [Google Scholar] [CrossRef]
  2. Szeląg, B.; Górski, J.; Bąk, Ł.; Górska, K. Modelling of stormwater quantity and quality on the example of urbanised catchment in Kielce. Ecol. Chem. Eng. A 2013, 20, 1305–1316. [Google Scholar]
  3. Llopart-Mascaró, A.; Farreny, R.; Gabarrell, X.; Rieradevall, J.; Gil, A.; Martínez, M.; Puertas, J.; Suárez, J.; Del Río, H.; Paraira, M. Storm tank against combined sewer overflow: Operation strategies to minimise discharges impact to receiving waters. Urban Water J. 2015, 12, 219–228. [Google Scholar] [CrossRef]
  4. Sałata, A.; Bąk, Ł.; Dąbek, L.; Ozimina, E. Assessment of the degree of pollution of sediments from the rainstorm sewer system in the urbanized catchment. Desalin. Water Treat. 2016, 57, 1478–1489. [Google Scholar] [CrossRef]
  5. Rauch, W.; Krejci, V.; Gujer, W. REBEKA—A software tool for planning urban drainage on the basis of predicted impacts on receiving waters. Urban Water 2002, 4, 355–361. [Google Scholar] [CrossRef]
  6. Lek, S.; Scardi, M.; Verdonschot, P.F.M.; Descy, J.P.; Park, Y.S. Modelling Community Structure in Freshwater Ecosystems; Springer: Berlin/Heidelberg, Germany, 2005. [Google Scholar]
  7. Bąk, Ł.; Michalik, A.; Tekielak, T. Appraise of consequences of erosive processes occurring in the Skawa river sector. Environ. Prot. Eng. 2011, 37, 83–92. [Google Scholar]
  8. Bąk, Ł.; Michalik, A.; Tekielak, T. The relationship between bank erosion, local aggradation and sediment transport in a small Carpathian stream. Geomorphology 2013, 191, 51–63. [Google Scholar] [CrossRef]
  9. Vaes, G.; Berlamont, J. Emission predictions with a multi–linear reservoir model. Water Sci. Technol. 1999, 39, 9–16. [Google Scholar] [CrossRef]
  10. US EPA. Combined Sewer Overflows. Guidance for Nine Minimum Controls; EPA 832-B-95-003; U.S. Environmental Protection Agency, Office of Water: Washington, DC, USA, 1995.
  11. Gamerith, V.; Bertrand-Krajewski, J.L.; Mourad, M.; Rauch, W. Implications of long-term stormwater quality modelling for design of combined sewer infrastructure. Urban Water J. 2011, 8, 155–166. [Google Scholar] [CrossRef]
  12. Mantegazza, S.A.; Gallina, A.; Mambretti, S.; Camylyn, L. Designing CSO storage tanks in Italy: A comparison between normative criteria and dynamic modelling methods. Urban Water J. 2010, 7, 211–216. [Google Scholar] [CrossRef]
  13. Zabel, T.; Milne, I.; Mckay, G. Approaches adopted by the European Union and selected Member States for the control of urban pollution. Urban Water 2001, 3, 25–32. [Google Scholar] [CrossRef]
  14. Rauch, W.; Henze, M.; Koncsos, L.; Reichert, P.; Shanahan, P.; SomlyóDy, L.; Vanrolleghem, P. River water quality modelling: I. state of the art. Water Sci. Technol. 1998, 38, 237–244. [Google Scholar] [CrossRef]
  15. US EPA. Methods for Measuring the Acute Toxicity of Effluents and Receiving Waters to Freshwater and Marine Organisms; EPA 821-R-02-012; U.S. Environmental Protection Agency, Office of Water: Washington, DC, USA, 2002.
  16. Schwoerbel, J.; Brendelberger, H. Einführung in die Limnologie. (Introduction to Limnology); Spektrum Akademie: Heidelberg, Germany, 2005. (In German) [Google Scholar]
  17. ATV-A 128E; Standards for the Dimensioning and Design of Stormwater Structures in Combined Sewers. GFA: Hennef, Germany, 1992.
  18. Fenz, R. Gewässerschutz bei Entlastungsbauwerken der Mischkanalisation; Wiener Mitteilungen, Band 174. Institut für Wassergüte und Abfallwirtschaft; Technische Universität: Wien, Austria, 2002. [Google Scholar]
  19. Dąbrowski, W. Estimating, computing and measuring multiplication factors for the working of storm overfalls. Gaz Woda Tech. Sanit. 2007, 11, 19–22. (In Polish) [Google Scholar]
  20. Fidala-Szope, M.; Sawicka-Siarkiewicz, H.; Koczyk, A. Protection of Surface Waters Against Stormwater Discharges from Combined Sewers: Decision Guide; Instytut Ochrony Środowiska: Warszawa, Poland, 1999. (In Polish) [Google Scholar]
  21. Vaes, G. The Influence of Rainfall and Model Simplification on Combined Sewer System Design. Ph.D. Thesis, Catholic University of Leuven, Leuven, Belgium, 1999. [Google Scholar]
  22. Butz, J. Stoffstrombilanzen für Phosphor und Sechs Schwermetalle am Beispiel des Oberen Kraichbachs; Universität Fridericiana zu Karlsruhe (TH): Karlsruhe, Germany, 2005. [Google Scholar]
  23. Leandro, J.; Chen, A.; Djordjević, S.; Savić, D. Comparison of 1D/1D and 1D/2D coupled (sewer/surface) hydraulic models for urban flood simulation. J. Hydraul. Eng. 2009, 135, 495–504. [Google Scholar] [CrossRef]
  24. Coutu, S.; Del Giudice, D.; Rossi, L.; Barry, D.A. Parsimonious hydrological modeling of urban sewer and river catchments. J. Hydrol. 2012, 464–465, 477–484. [Google Scholar] [CrossRef]
  25. Barco, J.; Wong, K.; Stenstrom, M. Automatic calibration of the U.S. EPA SWMM model for a large urban catchment. J. Hydraul. Eng. 2008, 134, 466–474. [Google Scholar] [CrossRef]
  26. Kim, J.; Park, J.; Cha, S.; Kwon, S. Applying low-impact development techniques for improved water management in urban areas. Water 2024, 16, 2837. [Google Scholar] [CrossRef]
  27. Ahmad, S.; Jia, H.; Ashraf, A.; Yin, D.; Chen, Z.; Ahmed, R.; Israr, M. A novel GIS-SWMM-ABM approach for flood risk assessment in data-scarce urban drainage systems. Water 2024, 16, 1464. [Google Scholar] [CrossRef]
  28. Andrés-Doménech, I.; Múnera, J.C.; Francés, F.; Marco, J.B. Coupling urban event-based and catchment continuous modelling for combined sewer overflow river impact assessment. Hydrol. Earth Syst. Sci. 2010, 14, 2057–2072. [Google Scholar] [CrossRef]
  29. Carleton, M.G. Comparison of overflows from separate and combined sewers—Quantity and quality. Water Sci. Technol. 1990, 22, 31–38. [Google Scholar] [CrossRef]
  30. Licznar, P. Stormwater reservoir dimensioning based on synthetic rainfall time series. Ochr. Srodowiska 2013, 35, 27–32. (In Polish) [Google Scholar]
  31. Arbeitsblatt DWA-A 110; Hydraulische Dimensionierung und Leistungsnachweis von Abwasserleitungen und—Kanälen. DWA: Hennef, Germany, 2006. Available online: https://shop.dwa.de/media/65/8d/ac/1726527742/a_110_08_2006.pdf (accessed on 3 December 2024).
  32. Kroll, S.; Blumensaat, F.; Dirckx, G.; Thoeye, C.; De Gueldre, G.; Van de Steene, B.; Tränckner, J. Assessment of CSO activity using simple volume balancing. In Proceedings of the 6th International Conference Novatech, Lyon, France, 25–28 June 2007. [Google Scholar]
  33. Cambez, M.J.; Pinho, J.; David, L.M. Using SWMM 5 in the continuous modelling of stormwater hydraulics and quality. In Proceedings of the 11th International Conference on Urban Drainage, Edinburgh, Scotland, 31 August–5 September 2008. [Google Scholar]
  34. Kiczko, A.; Romanowicz, R.J.; Osuch, M.; Karamuz, E. Maximising the usefulness of flood risk assessment for the River Vistula in Warsaw. Nat. Hazard. Earth Sys. 2013, 13, 3443–3455. [Google Scholar] [CrossRef]
  35. Kleidorfer, M.; Deletic, A.; Fletcher, T.D.; Rauch, W. Impact of input data uncertainties on urban stormwater model parameters. Water Sci. Technol. 2009, 60, 1545–1554. [Google Scholar] [CrossRef] [PubMed]
  36. Korving, H.; Van Gelder, P.; Van Noortwijk, J.; Clemens, F. Influence of model parameter uncertainties on decision-making for sewer system management. In Proceedings of the Fifth International Conference on Hydroinformatics, Cardiff, UK, 1–5 July 2002. [Google Scholar]
  37. Thorndahl, S.; Willems, P. Probabilistic modelling of overflow, surcharge and flooding in urban drainage using the first-order reliability method and parameterization of local rain series. Water Res. 2008, 42, 455–466. [Google Scholar] [CrossRef]
  38. Thorndahl, S. Stochastic long term modelling of a drainage system with estimation of return period uncertainty. In Proceedings of the 11th International Conference on Urban Drainage, Edinburgh, Scotland, 31 August–5 September 2008. [Google Scholar]
  39. Bąk, Ł.; Górski, J.; Górska, K.; Szeląg, B. Suspended solids and heavy metals content of selected rainwater waves in an urban catchment area: A case study. Ochr. Środ 2012, 34, 49–52. (In Polish) [Google Scholar]
  40. Fu, G.; Butler, D.; Khu, S.T.; Sun, S. Imprecise probabilistic evaluation of sewer flooding in urban drainage systems using random set theory. Water Resour. Res. 2011, 47, 1–13. [Google Scholar] [CrossRef]
  41. Fu, G.; Butler, D. Copula-based frequency analysis of overflow and flooding in urban drainage systems. J. Hydrol. 2014, 510, 49–58. [Google Scholar] [CrossRef]
  42. Fontanazza, C.M.; Freni, G.; La Loggia, G.; Notaro, V. Uncertainty evaluation of design rainfall for urban flood risk analysis. Water Sci. Technol. 2011, 63, 2641–2650. [Google Scholar] [CrossRef]
  43. Balistrocchi, M.; Bacchi, B. Modelling the statistical dependence of rainfall event variables through copula functions. Hydrol. Earth Syst. Sci. 2011, 15, 1959–1977. [Google Scholar] [CrossRef]
  44. Arbeitsblatt DWA-A 118; Hydraulische Bemessung und Nachweis von Entwässerungssystemen. DWA: Hennef, Germany, 2006.
  45. Kotowski, A.; Dancewicz, A.; Kaźmierczak, B. Accuracy of measurements of precipitation amount using standard and tipping bucket pluviographs in comparison to Hellmann rain gauges. Environ. Prot. Eng. 2011, 37, 23–34. [Google Scholar]
  46. Suligowski, R. Struktura Czasowa i Przestrzenna Opadów Atmosferycznych w Polsce. Próba Regionalizacji; Prace Instytutu Geografii Akademii Świętokrzyskiej w Kielcach: Kielce, Poland, 2004. [Google Scholar]
  47. Szeląg, B.; Kiczko, A.; Studzinski, J.; Dąbek, L. Hydrodynamic and probabilistic modelling of storm overflow discharges. J. Hydroinform. 2018, 20, 1100–1110. [Google Scholar] [CrossRef]
  48. Szeląg, B.; Suligowski, R.; Studzinski, J.; De Paola, F. Application of logistic regression to simulate the influence of rainfall genesis on storm overflow operations: A probabilistic approach. Hydrol. Earth Syst. Sci. 2020, 24, 595–614. [Google Scholar] [CrossRef]
  49. Butech Limited Liability Company. Construction of Stormwater Treatment on the Collector Si9, Building Project; Manuscript: Kielce, Poland, 2003. (In Polish) [Google Scholar]
  50. Szeląg, B.; Bąk, Ł. Probabilistic model for the annual number of storm overflow discharges in a stormwater drainage system. Urban Water J. 2017, 14, 604–611. [Google Scholar] [CrossRef]
  51. Huber, W.C.; Dickinson, R.E. Storm Water Management Model; Version 4: User’s Manual; EPA/600/3-88/001a; U.S. Environmental Protection Agency, Office of Research and Development: Athens, GA, USA, 1988.
  52. Rossman, L.A. Storm Water Management Model; User’s Manual. Version 5.0. EPA/600/R-05/040; U.S. Environmental Protection Agency, National Risk Management Research Laboratory, Office of Research and Development: Cincinnati, OH, USA, 2010.
  53. Fu, G.; Kapelan, Z. Flood analysis of urban drainage systems: Probabilistic dependence structure of rainfall characteristics and fuzzy model parameters. J. Hydroinform. 2013, 15, 687–699. [Google Scholar] [CrossRef]
  54. Brzezińska, A.; Zawilski, M.; Sakson, G. Assessment of pollutant load emission from combined sewer overflows based on the online monitoring. Environ. Monit. Assess 2016, 188, 502. [Google Scholar] [CrossRef]
  55. Yu, Y.; Kojima, K.; An, K.; Furumai, H. Cluster analysis for characterization of rainfalls and CSO behaviours in an urban drainage area of Tokyo. Water Sci. Technol. 2013, 68, 544–551. [Google Scholar] [CrossRef]
  56. Muhaisen, O.S.; Osorio, F.; García, P.A. Two-copula based simulation for detention basin design. Civ. Eng. Environ. Syst. 2009, 26, 355–366. [Google Scholar] [CrossRef]
  57. Osorio, F.; Muhaisen, O.; García, P.A. Copula-based simulation for the estimation of optimal volume for a detention basin. J. Hydrol. Eng. 2009, 14, 1378–1382. [Google Scholar] [CrossRef]
Figure 1. Diagram of the urbanized catchment of the Si9 sewer in the SWMM program.
Figure 1. Diagram of the urbanized catchment of the Si9 sewer in the SWMM program.
Water 16 03598 g001
Figure 2. Diagram of the separation chamber (DC).
Figure 2. Diagram of the separation chamber (DC).
Water 16 03598 g002
Figure 3. Comparison of the results of measurements and computations of the separation chamber for precipitation event No. 4.
Figure 3. Comparison of the results of measurements and computations of the separation chamber for precipitation event No. 4.
Water 16 03598 g003
Figure 4. Empirical and theoretical distribution functions for the maximum Pt30 in the precipitation event.
Figure 4. Empirical and theoretical distribution functions for the maximum Pt30 in the precipitation event.
Water 16 03598 g004
Figure 5. Empirical distribution function of the probability of not exceeding the annual number of overflow discharges.
Figure 5. Empirical distribution function of the probability of not exceeding the annual number of overflow discharges.
Water 16 03598 g005
Figure 6. Empirical distribution function of the probability of not exceeding the annual volume of overflow discharges.
Figure 6. Empirical distribution function of the probability of not exceeding the annual volume of overflow discharges.
Water 16 03598 g006
Figure 7. Empirical distribution function describing the probability that the maximum volume of wastewater discharge through the overflow is not exceeded in the precipitation event.
Figure 7. Empirical distribution function describing the probability that the maximum volume of wastewater discharge through the overflow is not exceeded in the precipitation event.
Water 16 03598 g007
Figure 8. Empirical distribution function describing the probability that the share of the maximum volume of wastewater overflow discharge is not exceeded in the precipitation event in relation to the overall annual discharge.
Figure 8. Empirical distribution function describing the probability that the share of the maximum volume of wastewater overflow discharge is not exceeded in the precipitation event in relation to the overall annual discharge.
Water 16 03598 g008
Table 1. Characteristics of precipitation events and parameters of fitting of the separation chamber water level, measured and obtained from simulations.
Table 1. Characteristics of precipitation events and parameters of fitting of the separation chamber water level, measured and obtained from simulations.
No.Hydrograph of the Flow from the CatchmentSeparation Chamber
Ptot
mm
td
min
Vtot
m3
δ
-
R
-
18.66017330.880.90
29.228622210.950.94
312.510719080.970.91
416.527034150.940.86
54.22621330.960.88
65.4566840.910.87
73.6923270.920.90
Note: Ptot—total rainfall depth in the precipitation event; td—rainfall duration; Vtot—volume of hydrograph of flow from catchment; δ—ratio of the maximum measured water level to the maximum simulated water level; R—correlation coefficient.
Table 2. Values of the correlation coefficients between precipitation parameters and discharge volume.
Table 2. Values of the correlation coefficients between precipitation parameters and discharge volume.
PtotPt30Pt15Pt10tdtapVd
Ptot1.000.410.290.250.59−0.380.62
Pt30-1.000.950.75−0.26−0.130.92
Pt15--1.000.87−0.39−0.050.85
Pt10---1.00−0.43−0.020.83
td----1.00−0.27−0.05
tap-----1.00−0.21
Vd------1.00
Note: Ptot—total rainfall depth in the precipitation event; Pt30—maximum 30 min rainfall depth; Pt15—maximum 15 min rainfall depth; Pt10—maximum 10 min rainfall depth; td—rainfall duration; tap—antecedent period duration; Vd—volume of overflow discharge.
Table 3. Results of computations (p) of Kolmogorov–Smirnov and chi-square tests of the fitting of the empirical distribution function to the theoretical distribution function obtained with the use of statistical distributions of concern.
Table 3. Results of computations (p) of Kolmogorov–Smirnov and chi-square tests of the fitting of the empirical distribution function to the theoretical distribution function obtained with the use of statistical distributions of concern.
Distribution Kolmogorov–Smirnov TestChi-Square Test
Beta0.00010.0002
Chi-square0.00020.0003
Weibull0.02530.0312
Exponential0.00070.0005
GEV0.00020.0003
Fisher–Tippett0.00030.0004
Log-normal0.31300.2526
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

Górski, J.; Szeląg, B.; Bąk, Ł.; Świercz, A. A Probabilistic Model for Predicting the Performance of a Stormwater Overflow Structure as Part of a Stormwater Treatment Plant. Water 2024, 16, 3598. https://doi.org/10.3390/w16243598

AMA Style

Górski J, Szeląg B, Bąk Ł, Świercz A. A Probabilistic Model for Predicting the Performance of a Stormwater Overflow Structure as Part of a Stormwater Treatment Plant. Water. 2024; 16(24):3598. https://doi.org/10.3390/w16243598

Chicago/Turabian Style

Górski, Jarosław, Bartosz Szeląg, Łukasz Bąk, and Anna Świercz. 2024. "A Probabilistic Model for Predicting the Performance of a Stormwater Overflow Structure as Part of a Stormwater Treatment Plant" Water 16, no. 24: 3598. https://doi.org/10.3390/w16243598

APA Style

Górski, J., Szeląg, B., Bąk, Ł., & Świercz, A. (2024). A Probabilistic Model for Predicting the Performance of a Stormwater Overflow Structure as Part of a Stormwater Treatment Plant. Water, 16(24), 3598. https://doi.org/10.3390/w16243598

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