Rapid urbanization affecting demand patterns, coupled with potential water shortages due to supply side impacts of climatic changes, has led to the emergence of new technologies for water and wastewater reuse. Sewer mining (SM) is a novel decentralized option that could potentially provide non-potable water for urban uses, including for example the irrigation of urban green spaces, providing a mid-scale solution to effective wastewater reuse. SM is based on extracting wastewater from local sewers and treatment at the point of demand and entails in some cases the return of treatment residuals back to the sewer system. Several challenges are currently in the way of such applications in Europe, including public perception, inadequate regulatory frameworks and engineering issues. In this paper we consider some of these engineering challenges, looking at the sewer network as a system where multiple physical, biological and chemical processes take place. We argue that prior to implementing SM, the dynamics of the sewer system should be investigated in order to identify optimum ways of deploying SM without endangering the reliability of the system. Specifically, both wastewater extraction and sludge return could result in altering the biochemical process of the network, thus unintentionally leading to degradation of the sewer infrastructure. We propose a novel Monte-Carlo based method that takes into account both spatial properties and water demand characteristics of a given area of SM deployment while simultaneously accounting for the variability of sewer network dynamics in order to identify potential locations for SM implementation. The outcomes of this study suggest that the method can provide rational results and useful guidelines for upscale SM technologies at a city level.
Rapid urbanization and potential water shortages due to, inter alia, climatic variability have led to the emergence of new technologies in water and wastewater reuse aiming to provide alternative water sources for more resilient cities. Sewer mining (SM) is such a technology, which is based on extracting wastewater from local sewers for reuse applications (after treatment) and (often) returning treatment residuals (sludge) back to the system (e.g., Sydney Water 2008). Typical uses of this recycled water are industrial cleaning and cooling, as well as irrigation of urban green spaces (Hadzihalilovic 2009; Marleni et al. 2012). Literature classifies this technology as a decentralized option (Makropoulos et al. 2017) because it is applicable (and suitable) at a development level (for example, up to 5,000 households). This is also highlighted in Marleni et al. (2012) where it is argued that this practice is not intended for individual use (indoor appliance) rather than for collective/cluster scale developments. Furthermore, the latter authors remark that usually such systems are not managed by central water utilities (or governmental organizations) but by private establishments under some license agreements. As such, SM is a promising reuse option that lies in the interplay between reuse at household scale (e.g., grey water reuse; cf. Liu et al. (2010) and Makropoulos & Butler (2010)) and centralized reuse at the wastewater treatment plant (WWTP) level (Andreadakis et al. 2006). Current SM projects mostly involve park and sports fields' irrigation. Most of them are operating in Australia (Sydney Water 2009) where the climate is dry and water should be treated carefully. It is worth highlighting that in most cases the treated water is used for non-drinking purposes. Despite public perception, concerns and inadequate regulatory frameworks that may raise potential barriers for SM implementation, there are engineering issues that have to be addressed. A sewer network is a system where multiple physical, biological and chemical processes take place (Pomeroy 1990). Hence, prior to implementing SM, the dynamics of the system should be investigated in order to identify optimum ways of deploying SM without endangering the reliability of the system. Specifically, both wastewater extraction and sludge return could result in altering the biochemical process of the network, thus unintentionally leading to degradation of the infrastructure. In this paper we focus on addressing some of the engineering challenges linked to the potential deployment of such technologies at the city scale. Typical engineering issues associated with SM are odour and corrosion. Both of them are related to the production of hydrogen sulphide (H2S) in sewer pipes. This study focuses on identifying possible locations for SM placement subject to minimizing H2S build-up.
MATERIALS AND METHODS
While trying to address this issue, i.e., SM placement by taking also into consideration H2S production, we propose a methodology consisting of three steps: (I) a spatial data pre-processing step during which the spatial properties and water demand characteristics are being identified, (II) a Monte-Carlo simulation (MCS) step, which involves the simulation of the sewer network in order to account for the variability of sewage discharge into the network and, finally, (III) a post-processing step which comprises the definition of appropriate metrics that quantify the output of interest (III-a) and a multi-criteria analysis of the results (III-b). A schematic description of the proposed methodology is given in Figure 1. During the first step the available spatial information (i.e., sewage network topology and assets, topography, water and land uses) is imported into the procedure in order to identify land uses that will benefit from SM (in our case green areas and parks). It involves a procedure of locating neighbouring sewer network components (e.g., nodes) which are close to areas of interest. In more detail, this is done by delineating a wider area surrounding the original one (e.g., add 10 m offset to green areas) and subsequently identifying the nodes that lie in those wider areas. Finally, the paths from the identified nodes to an ‘exit’ node are identified and stored. The exit node could be a WWTP or a node that links the understudy network with a broader larger network. It is worth noticing that this path is unique for each node due to the ‘collective nature’ of sewer networks. The purpose of the second step is to propagate uncertainties related to the input parameters to the quantities of interest (e.g., 5-day biochemical oxygen demand (BOD5) concentration or flow of each pipe). Furthermore, the use of MCS allows the use of probabilistic functions and metrics, which in-turn provide uncertainty-aware outputs. Typical examples of uncertain parameters are the daily water consumption, daily and hourly variation coefficients of wastewater discharge and BOD5 loading (in terms of g/cap). Alternatively, one could use a similar scenario based approach to sample those parameters (or in conjunction with MCS) in order to investigate the effect of certain predefined scenarios (e.g., worst, base, favorable conditions).
The third and final step involves the definition and the use of metrics, i.e., utility functions or risk functions that quantify the output of interest, in our case H2S build-up, for a chain of pipes (the paths specified in step I). We remark that BOD5 can be directly associated with H2S through empirically derived relationships (e.g., Lahav et al. 2006; Marleni et al. 2015). Furthermore, as a final procedure, we use multi-criteria analysis which eventually leads to derivation of a Pareto front (based on conflicting criteria, e.g., suitability of location and green area water demand), which includes all the potential locations for SM.
The methodology is demonstrated in a sewer network designed for the city of Kalyvia Thorikou in Greece (Figure 2). The network has not been constructed yet, although it is foreseen to accommodate an area of 98 ha, from which 17 ha are green areas. It is part of a larger engineering project of Saronikos municipality (service 10–15 thousand people) which aims to extend the existing sewage network of the coastal zone. It consists of 1,030 pipes of total length ∼38 km, while their diameter varies from 0.2 m to 0.5 m. The pipe slope varies from 2‰ to 150‰, with an average slope of 35‰. The understudy area can be considered appropriate for testing the proposed methodology, since it consists of various network elements and has adequate number of green areas which could benefit from SM practices.
The design period of the network was assumed, T = 40 years, as in the original study of Hydroexigiantiki, the engineering firm which conducted the study of the above network. The design population (E) is adjusted using the compound rate formula , where Eois the current population, ɛ is the increase rate (assumed 1.5%) and n is the extrapolation year (n = 0, …, T). The value of n can be varied in order to assess the performance of the system at different time periods. In this study, q was assumed to be equal to 250 lpd, λL was assumed equal to 0.725 for year 0 and 0.85 for year 40. Similarly, λs was assumed equal to 0.625 for year 0 and 0.65 for year 40. The values of λL and λs for intermediate years can be calculated using linear interpolation. The value of λDWF was set equal to 0.2. Finally, we assumed λ1 and λ2 as uncertain parameters that follow uniform distribution; i.e., we assumed λ1 ∼ Uniform[0.7, 1.3] and λ2 ∼ Uniform[0.8, 1]. Concerning parameter n, we employed three scenarios, 0, 20 and 40 years. Also, the mass of BOD5 was varied using three scenarios 40, 50, and 65 g/(day cap). The maximum allowable number of simulation runs for the MCS step was set equal to 500. The desired quantile x (i.e., reliability level) for the calculation of Q[MZc]x was set to 75%.
RESULTS AND DISCUSSION
Figure 3 illustrates the final result of the post-processing step III in a form of a Pareto front, using as objectives the minimization of modified Z index and the maximization of green area. It is notable that one could also interpret those two objectives as the simultaneous maximization of suitability and benefit from SM practices respectively. The suggested procedure located three potential locations for SM units' placement that optimize both criteria simultaneously, while on the other hand it discarded other inferior locations. Additionally, the map depicted in Figure 4 provides a visual summary of all the green areas (light-grey polygons, green in online version) of the case study, as well as the three areas (dark-grey polygons, red in online version) identified by the proposed methodology since they were suitable for SM placement. Furthermore, in order to visually illustrate the concept of optimum path it presents the selected optimum path (light-grey line, magenta in online version) for the green area with ID3. This path has the lowest MZ value compared to all other alternative paths of ID3.
Figure 5 depicts the cross-section of optimum path of green area ID3 (light-grey/magenta line in Figure 4). The path starts from pipe C215 which is located close to the green area ID3 and ends at C122 which is linked with the ‘exit’ node of the understudy system. More specifically, the upper panel of Figure 5 shows the variability of the MZ across that path. Furthermore, the lower panel of Figure 5 shows the probability (P) of non-exceedance of the threshold values (Z < 7,500). It can be seen that until C171 the system demonstrates high non-exceedance probabilities (∼90%), i.e., high reliability. After that point the reliability decreases but it is still preserved within acceptable levels (70–80%).
In order to overcome the engineering challenges imposed by the multiple physical, biological and chemical processes that take place in a sewer network, we introduced a novel Monte-Carlo based method for the identification of potential locations for SM units. The proposed risk based approach allows the safe planning of SM deployments taking into due consideration system performance objectives regarding water quantity and quality. As such it can be used to enhance the decision making process with useful guidelines and insights. More specifically, the proposed method has been demonstrated through a case study (Kalyvia Thorikou, Greece) where we focused on identifying optimum locations for SM units subject to the generation (minimize) of hydrogen sulphide and water demand. The results showed that the proposed methodology was able to identify potential locations for SM units' placement while simultaneously taking into consideration the spatial properties of the area as well as the variability and hydraulic characteristics of the sewer network. Future work will focus on improving the proposed framework through the integration of a dynamic simulation model, such as SWMM 5.0, into the computational procedure.
The research leading to these results has received funding from the European Union Seventh Framework Programme under grant agreement no. 619039 (ENV.2013.WATER INNO&DEMO-1), for the research project DESSIN ‘Demonstrate Ecosystem Services Enabling Innovation in the Water’. The research and its conclusions reflect only the views of the authors and the European Union is not liable for any use that may be made of the information contained herein. Also, we would like to sincerely thank Hydroexigiantiki Consulting Engineers for providing the data of the understudy area (Kalyvia Thorikou, Greece). An earlier version of this paper was presented at the 13th IWA Specialized Conference on Small Water and Wastewater Systems (Tsoukalas et al. 2016).