water
Article Empirical Sewer Water Quality Model for Generating Influent Data for WWTP Modelling
Jeroen Langeveld 1,2,* , Petra Van Daal 2,3, Remy Schilperoort 1, Ingmar Nopens 4 , Tony Flameling 5 and Stefan Weijers 5
1 Partners4UrbanWater, Javastraat 104a, Nijmegen 6524 MJ, The Netherlands; [email protected] 2 Section Sanitary Engineering, Department of Water Management, Delft University of Technology, P.O. Box 5048, Delft 2628 CD, The Netherlands; [email protected] 3 Witteveen+Bos, P.O. Box 233, Deventer 7400 AE, The Netherlands 4 BIOMATH, Department of Mathematical Modelling, Statistics and Bioinformatics, Ghent University, Coupure Links 653, Gent 9000, Belgium; [email protected] 5 Waterschap De Dommel, P.O. Box 10.001, Boxtel 5280 DA, The Netherlands; tfl[email protected] (T.F.); [email protected] (S.W.) * Correspondence: [email protected]; Tel.: + 31-6-1897-6283
Received: 30 March 2017; Accepted: 30 June 2017; Published: 5 July 2017
Abstract: Wastewater treatment plants (WWTP) typically have a service life of several decades. During this service life, external factors, such as changes in the effluent standards or the loading of the WWTP may change, requiring WWTP performance to be optimized. WWTP modelling is widely accepted as a means to assess and optimize WWTP performance. One of the challenges for WWTP modelling remains the prediction of water quality at the inlet of a WWTP. Recent applications of water quality sensors have resulted in long time series of WWTP influent quality, containing valuable information on the response of influent quality to e.g., storm events. This allows the development of empirical models to predict influent quality. This paper proposes a new approach for water quality modelling, which uses the measured hydraulic dynamics of the WWTP influent to derive the influent water quality. The model can also be based on simulated influent hydraulics as input. Possible applications of the model are filling gaps in time series used as input for WWTP models or to assess the impact of measures such as real time control (RTC) on the performance of wastewater systems.
Keywords: sewer system; empirical model; influent generator; influent modelling; gap filling
1. Introduction Modelling of wastewater treatment plants (WWTPs) using activated sludge models (ASM) has become a standard in both industry and academia for a range of objectives, amongst others WWTP design [1,2], operation [3], and control [4,5]. Especially the latter objective, developing control strategies and assessing the performance of control strategies requires high frequency influent data [6]. In a comprehensive review, [7] discussed the available approaches for generating influent data. The approaches range from (1) data driven methods based on creating databases with monitoring and experimental data and derived models using the data; (2) very simple models based on harmonic functions; and (3) phenomenological models [8,9]. The data driven methods comprise two different approaches. The first method [10] uses pollutant release patterns derived from literature to generate dynamic influent data aggregating the punctual emissions from the database. The second method [11] interpolates available influent data at e.g., a daily timescale to e.g., hourly dynamics.
Water 2017, 9, 491; doi:10.3390/w9070491 www.mdpi.com/journal/water Water 2017, 9, 491 2 of 18
The simple models based on harmonic functions are very well suited for the analyses of dry weather flow (DWF) situations, but less so for wet weather flow (WWF) situations [7]. The phenomenological models are the most detailed influent models, that can give a phenomenological representation of dynamics of WWTP influent, including diurnal patterns, weekend, seasonal, and holiday variations as well as rain events [7,8]. Despite being labelled as ‘promising’ [7], todays phenomenological models cannot adequately reproduce the dynamics in WWTP influent during wet weather due to a relatively poor representation of the build-up and wash off of urban pollutants. This is also true for the most recently published influent generator by [12], who use a mix of statistical and conceptual modelling techniques for synthetic generation of influent time series. The limitations associated with the use of influent generators are not surprising given the state of the art knowledge on the physical-chemical, biological and transport processes occurring in sewer systems [13–17]. Sediment transport especially is not very well understood and not very successfully reproduced in deterministic sewer models. This is partly due to the fact that it is currently not possible to get enough data on the initial sewer sediment conditions throughout an entire sewer network. It is interesting to note that the developers of influent generators, showing many similarities with simplified or parsimonious sewer models [18,19], are facing the same issues as sewer modelers in the past, i.e., how to incorporate the contribution of in-sewer stocks during storm events to the outflow of sewers via either combined sewer overflow (CSO) or WWTP influent. In order to overcome the limitations of deterministic sewer models, regression models have been proposed, which are validated against monitoring data. A recent successful example of this approach is given by [20], who developed an empirical model for storm water total suspended solids (TSS), event mean concentrations with rainfall depth, and antecedent dry weather period as input variables. These empirical relations, that are valid at a CSO or storm sewer outfall (SSO), however, are not suitable for the prediction of WWTP influent quality, as these models do not predict the influent quality during DWF. For WWTP influent modelling, the empirical model described by [21], relates the influent concentration to the daily flows. A weak point of this approach is the impossibility to account for the dynamics during storm events. This is a major drawback, as storm events typically do not last a full day. Recent applications of water quality sensors have resulted in the availability of long time series of WWTP influent quantity and quality [22,23]. These time series contain a lot of information on the response of influent quality to storm events and the contribution of in-sewer stocks to WWTP influent [24]. In this study, time series analysis is used to understand the dynamics of WWF related variations in WWTP influent quality and to relate the variation in influent quality to influent hydraulics. This allowed the development of an empirical model based on understanding of the underlying physical processes. The paper is organized as follows: first, the available data set of WWTP Eindhoven is described. Second, the model development and calibration method are presented. Next, the calibration results, model results and the transferability of the concept are presented, discussed, and finally, some foreseen applications of the model are introduced.
2. Materials and Methods
2.1. System Description: The Dommel River IUWS The Dommel river is relatively small and sensitive to loadings from CSOs and WWTP effluent. The river flows through the city of Eindhoven (the Netherlands) from the Belgian border (south) into the River Meuse (north). The Dommel receives discharges from the 750,000 people equivalent (PE) WWTP of Eindhoven and from over 200 CSOs in 10 municipalities. In summer time, the base flow in the river just downstream of the WWTP comprises 50% of WWTP effluent, increasing to 90% during small storm events. The Dommel River does not yet meet the requirements of the European Water 2017, 9, 491 3 of 18 Water 2017, 9, 491 3 of 18
Unionoxygen Water (DO) Framework depletion, Directiveammonia (WFD)peaks [and25]. seasonal The water average quality nutrient issues to concentration be addressed arelevels dissolved [26,27]. oxygenEarlier (DO)research depletion, within ammoniathe KALLISTO peaks project and seasonal [18,25] averagedemonstrated nutrient that concentration the WWTP effluent levels [26 is,27 the]. Earliermain source research for within the toxic the KALLISTOammonia peaks project in [the18,25 Dommel] demonstrated river and that that the the WWTP ammonium effluent peaks is the in main the sourceWWTP for effluent the toxic can ammonia be significantly peaks in reduced the Dommel by applied river and integrated that the ammoniumreal time control peaks (RTC). in the WWTPIn [28], effluentthe use of can RTC be significantlyby activating reducedin-sewer bystorage applied volume integrated to reduce real and time delay control the hydraulic (RTC). In [peak28], theloading use ofof RTCthe byWWTP activating during in-sewer storm storageevents volumehas been to shown reduce andto be delay an effective the hydraulic measure. peak Reference loading of the[29] WWTPintroduced during a new storm RTC events concept: has been the smart shown buffer, to be anwh effectiveich minimizes measure. the Referencepeak load [ 29to] the introduced biology at a newWWTP RTC Eindhoven concept: the by smart applying buffer, the which aforementioned minimizes the RTC peak combined load to thewith biology using atonly WWTP one Eindhovenof the three byprimary applying clarifiers the aforementioned (PC) during dry RTC weather combined (DWF) with and using using only the one other of thetwo three PCs only primary during clarifiers storm (PC)events. during dry weather (DWF) and using the other two PCs only during storm events. TheThe 1010 municipalitiesmunicipalities contributingcontributing to to thethe WWTPWWTP influentinfluentare are divideddivided overover threethree catchmentcatchmentareas areas thatthat are are very very different different in in size size and and character, character, each each having having a separatea separate inflow inflow to theto the WWTP WWTP (see (see Figure Figure1). Wastewater1). Wastewater from from Eindhoven Eindhoven Stad Stad (ES, municipality(ES, municipality of Eindhoven) of Eindhoven) accounts accounts for approximately for approximately 50% 3 (in50% practice (in practice ranging ranging between between 14,000 14,000 and and 17,000 17,000 m /h) m3/h) of of the the hydraulic hydraulic capacity capacity and and is is discharged discharged directlydirectly toto thethe WWTP.WWTP. The The other other nine nine (much (much smaller) smaller) municipalities municipalities are are each each connected connected to toone oneof ofthe the twotwo wastewaterwastewater transporttransport mains, one one to to the the north north (N (Nuenen/Sonuenen/Son or or NS, NS, 7 km 7 km in inlength) length) and and one one to the to 3 thesouth south (Riool (Riool Zuid Zuid or RZ, or RZ, 32 32km km in inlength), length), accounting accounting for for respectively respectively 7% 7% (3000 (3000 m3 m/h)/h) and and 43% 43% (in 3 (inpractice practice ranging ranging from from 14,000 14,000 to to15,000 15,000 m3 m/h)/h) of the of thehydraulic hydraulic capacity. capacity. An Anelaborate elaborate description description of the of thestudied studied wastewater wastewater system system can can be befound found in [23]. in [23 ].
Closing section catchment area: ps Son river Nuenen/Son ps Nederwetten Eindhoven Stad Riool-Zuid ps Gerwen ps Nuenen
wwtp cs rwzi Eindhoven ps Eeneind ps Wintelre ps Mierlo
sludge Veldhoven Geldrop proc. inst. ps Aalst ps Knegsel cs de Meren Aalst ps Heeze ps Waalre ps Steensel cs Valkenswaard ps Eersel ps Sterksel ps Duizel Valkenswaard ps Riethoven ps Leende
Bergeijk ps Westerhoven wwtp Eindhoven ps Weebosch free flow conduits ps Borkel pressure mains pumping stations control stations Luijksgestel sludge proc. inst. N 0 3 6 9 km
FigureFigure 1.1.Wastewater Wastewater system systemof of Eindhoven Eindhoven( Left(Left)) and and its its receiving receiving streams streams and and schematic schematic lay lay out out of of thethe wastewater wastewater system; system; ( Right(Right)) Figure Figure reproduced reproduced with with permission permission from from [ 23[23].].
2.2. Monitoring Network and Data Validation 2.2. Monitoring Network and Data Validation At each of the three inflows into the WWTP (locations ‘A’ in Figure 1 at the right) on-line At each of the three inflows into the WWTP (locations ‘A’ in Figure1 at the right) on-line spectroscopy sensors (UV-VIS) have been installed that measure equivalent concentration values of spectroscopy sensors (UV-VIS) have been installed that measure equivalent concentration values of wastewater quality parameters: total suspended solids (TSS), chemical oxygen demand (COD), and wastewater quality parameters: total suspended solids (TSS), chemical oxygen demand (COD), and filtered COD (CODf), i.e., the dissolved fraction of COD, at an interval of 2 min. In addition, flow is filtered COD (CODf), i.e., the dissolved fraction of COD, at an interval of 2 min. In addition, flow is recorded every minute at these locations and ammonium (NH4, using Hach Lange Amtax sensors) at recorded every minute at these locations and ammonium (NH4, using Hach Lange Amtax sensors) the Eindhoven Stad and Riool Zuid catchments. In this study, monitoring data for the year 2012 was at the Eindhoven Stad and Riool Zuid catchments. In this study, monitoring data for the year 2012 used. was used. The monitoring data were validated manually, focusing on obtaining reliable data for calibration The monitoring data were validated manually, focusing on obtaining reliable data for calibration of of wet weather flow processes. Figure 2 shows an example of data and their evaluation. After wet weather flow processes. Figure2 shows an example of data and their evaluation. After validation, validation, only 38.5% of the data was considered to have an acceptable quality during the condition required. This percentage of data perceived ‘good enough’ after validation may seem relatively low.
Water 2017, 9, 491 4 of 18
Water 2017, 9, 491 4 of 18 only 38.5% of the data was considered to have an acceptable quality during the condition required. ThisDuring percentage earlier research of data projects perceived (2007–2008) ‘good enough’ at WWTP after Eindhoven validation on may UV-VIS seem sensors, relatively the low. percentage During ofearlier ‘good research enough’ projects data after (2007–2008) data validation at WWTP ranged Eindhoven between on UV-VIS 50% and sensors, 75%, thedespite percentage very intensive of ‘good maintenanceenough’ data afterand datasurveillance validation [23] ranged and betweenwithout 50%restrictions and 75%, on despite the influent very intensive conditions. maintenance WWTP influentand surveillance has shown [23 ]to and be withouta very difficult restrictions medium on the for influent water conditions. quality monitoring. WWTP influent The dataset has shown after validationto be a very comprises difficult mediumapproximately for water 30 storm quality events monitoring. with good The data dataset for each after calibration validation performed. comprises Inapproximately the model calibration, 30 storm events these with events, good including data for eachthe antecedent calibration performed.dry day and In several the model following calibration, dry days,these events,were used. including In Figure the antecedent2, and all other dry dayFigures, and severalthe data following used for calibration dry days, were is represented used. In Figure by the2, darkand all grey other bullets, Figures, data the not data used used for calibration for calibration with islight represented grey bullet bys.the An darkassessment grey bullets, of routine data 24 not h waterused for quality calibration samples with of WWTP light grey influent bullets. showed An assessment that the DWF of routine does not 24 hshow water a noticeable quality samples seasonal of WWTPpattern. influent showed that the DWF does not show a noticeable seasonal pattern.
Figure 2. Example of quality evaluation of monitoring data.
2.3. Data Analysis In earlier work [16], [16], a part of this datadata setset was used to studystudy thethe dynamicsdynamics of wastewater composition. ThisThis resulted resulted in wellin well described described typical typical diurnal diurnal patterns patterns during DWFduring and DWF typical and dynamics typical dynamicsduring WWF during (Figure WWF3). (Figure 3).
Figure 3. Wet weather flow (WWF) dynamics during the stages of a storm event on 12 June 2008 in Figure 3. Wet weather flowflow (WWF) dynamics during the stages of a storm event on 12 June 2008 in Eindhoven Stad catchment. Eindhoven Stad catchment.
WaterWater2017 2017, 9, ,9 491, 491 55 of of 18 18
For WWF, it has been observed that the concentration levels of the wastewater show a typical patternFor WWF,during it a has storm been event: observed a short that period the concentration called ‘onset’ levels of the of storm the wastewater event, with show an aincreased typical patternconcentration during level a storm for event:particulate a short matter period but called not ‘onset’for dissolved of the stormmatter, event, a longer with anperiod increased called concentration‘dilution’, where level dilution for particulate of both matter dissolved but not and for dissolvedparticulate matter, matter a longertakes place, period and called ‘recovery’, ‘dilution’, a whereperiod dilution where dissolved of both dissolved and particulate and particulate matter slowly matter return takes to place, DWF and levels. ‘recovery’, a period where dissolved and particulate matter slowly return to DWF levels. 2.4. Model Development 2.4. Model Development The general idea behind the model development is that the measured hydraulic influent data, i.e., Theflow general and water idea level behind in the the influent model pumping development station, is that can the be used measured to make hydraulic a distinction influent between data, i.e.,the flowfour andpatterns: water DWF, level onset in the of influent WWF, pumpingdilution du station,ring WWF can beand used recovery to make after a distinction WWF. Each between of these thepatterns four patterns: is denoted DWF, as a onsetsystem of state, WWF, during dilution which during a certain WWF relation and recovery between after flow WWF. and Eachconcentration of these patternslevel applies. is denoted This allows as a system the incorporation state, during of which the contribution a certain relation of in-sewer between stocks flow on and top concentration of the mixing levelprocess applies. between This allowswastewater the incorporation and stormwater. of the The contribution latter is ofa common in-sewer stocksfeature on of top influent of the mixingmodels processapplied between to simulate wastewater both dry and and stormwater. wet periods, The while latter explicitly is a common accounting feature of for influent the contribution models applied of in- tosewer simulate stocks both circumventing dry and wet the periods, relatively while limi explicitlyted knowledge accounting associated for the in contributionsewer processes. of in-sewer stocksThe circumventing water quality the average relatively dry limited weather knowledge diurnal pattern associated is the in core sewer of processes.the model. As long as the systemThe state water is quality‘DWF’, average the average dry weatherdry weather diurnal diurnal pattern pattern is the based core of on the monitoring model. As data long is as used, the systemtogether state with is the ‘DWF’, measured the average flow data. dry weatherThe aver diurnalage dry pattern weather based diurnal on monitoringpattern has databeen isderived used, togetherfrom flow with monitoring the measured data flowby averaging data. The the average monitori dryng weather data of diurnal 10 dry patterndays over has 5 been min derivedintervals from with flowthe same monitoring timestamp. data by averaging the monitoring data of 10 dry days over 5 min intervals with the same timestamp.During wet weather, the model superimposes a number of processes on the DWF pattern for waterDuring quality wet to weather,mimic onset, the modeldilution, superimposes and recovery. a numberThe type of of processes parameter on (NH the4 DWF or COD) pattern and forthe watertype of quality event to(small, mimic medium onset,dilution, or large) anddetermine recovery. wh Theich of type these of processes parameter is (NH to be4 orapplied. COD) The and type the typeof event of event is used (small, as characteristic medium or large) of a storm determine event, which as it was of these found processes that the isrelation to be applied.between The flow type and ofconcentration event is used levels as characteristic differs very of amuch storm event,between as itsmall, was foundmedium, that theand relation large storm between events. flow andThe concentrationmeasured hydraulics, levels differs in this very case much the betweenflow and small, water medium, level in the and influent large storm pumping events. station, The measured are used hydraulics,to determine in thiswhich case ofthe the flow described and water processes level inshould the influent be activated pumping in the station, mode arel, using used the to determine scheme of whichFigure of 4. the described processes should be activated in the model, using the scheme of Figure4.
yes