Urban Drainage Water Quality Modeling on the SWMM Software, Northeastern Sector, Santa Inés District, Tunja
Modelación de calidad de agua del drenaje urbano en el software SWMM, sector nororiental, Santa Inés, Tunja
Palabras clave:alcantarillado, calidad del agua, modelado, SWMM (es).
Palabras clave:sewerage, water quality, modeling, SWMM (en).
Amaya-Tequia, E. (2019). Modelización hidráulica de drenaje urbano. Aplicación sector nororiental distrito Santa Inés Tunja-Boyacá [Undergraduate theis, Unversidad Santo Tomás]. http://hdl.handle.net/11634/22086
Cortés, M., & Vera, E. (2019a). Durability considerations for mitigating corrosion caused by biogenic sulfuric acid in sewerage system. Journal of Physics: Conference Series, 1386, 012088. http://dx.doi.org/10.1088/1742-6596/1386/1/012088 DOI: https://doi.org/10.1088/1742-6596/1386/1/012088
Cortés, M., & Vera, E. (2019b). Methodology for simulating the process of microbiologically induced degradation of concrete by chemical exposure of speccimens under controlled laboratory conditions. Journal of Physics: Conference Series, 1386, 012089. https://doi.org/10.1088/1742-6596/1386/1/012089 DOI: https://doi.org/10.1088/1742-6596/1386/1/012089
Cortés-Zambrano, M., Amaya-Tequia, E., & Gamba-Fernández, D. (2020). Implementation of the hydraulic modelling of urban drainage in the northeast sector, Tunja, Boyacá. Revista Facultad de Ingeniería Universidad de Antioquia, 101, 74-83. https://doi.org/10.17533/udea.redin.20200798 DOI: https://doi.org/10.17533/udea.redin.20200798
Gironás, J., Roesner, L. A., & Davis, J. (2009). Storm Water Management Model Applications Manual. United States Environmental Protection Agency. http://nepis.epa.gov/Exe/ZyPURL.cgi?Dockey=P100EA7H.txt
Hogue, T. S., Gupta, H., & Sorooshian, S. (1988). A 'user-friendly' approach to parameter estimation in hydrologic models. Journal of Hydrology, 320(1-2), 202-217. https://doi.org/10.1016/j.jhydrol.2005.07.009 DOI: https://doi.org/10.1016/j.jhydrol.2005.07.009
Pérez-Carmona, R. (2013). Diseño y construccion de alcantarillados sanitario, pluvial y drenaje en carreteras. Ecoe Ediciones.
Rodríguez, J. (2005). Desarrollo y aplicación de un modelo de calidad del agua del interceptor Fucha - Tunjuelos - Canoas [Master's thesis, Universidad de los Andes]. https://repositorio.uniandes.edu.co/bitstream/handle/1992/8998/u270910.pdf?sequence=1
Romero-Rojas, J. A. (2004). Tratamiento de aguas residuales. Escuela Colombiana de Ingeniería Julio Garavito.
Seco, R., & Gómez-Valentín, M. (2011). Aplicación del modelo SWMM 5.0 para la estimación del arrastre de depósitos de sedimentos y contaminantes asociados en sistemas de alcantarillado unitario [Conference presentation]. IWA Spain National Young Water Professionals Conference, Madrid, Spain. http://hdl.handle.net/2117/18111
Temprano, J., & Tejero, I. (2002). Detention storage volume for combined sewer overflow into a river. Environmental Technology, 23(6), 663-675. https://doi.org/10.1080/09593332308618381 DOI: https://doi.org/10.1080/09593332308618381
Roca-Hernández, R. (2012). Estudio de la eliminacion y de la prevencion de la generacion de sulfuros en redes de saneamiento [Master's thesis, Universidad Politecnica de Cartagena]. https://www.upct.es/hidrom/publicaciones/Tesis_pfc/redes_de_saneamiento%20.pdf
Obermann, M., Rosenwinkel, K., & Tournoud, M. (2009). Investigation of first flushes in a medium-sized mediterranean catchment. Journal of Hydrology, 373(3), 405-415. http://dx.doi.org/10.1016/j.jhydrol.2009.04.038 DOI: https://doi.org/10.1016/j.jhydrol.2009.04.038
Marcor, J., & Pedraza, R. (2012). Efectos de la discretización en la simulación de la escorrentía urbana. Ingeniería del Agua, 13(1), 35-46. https://doi.org/10.4995/ia.2006.2882 DOI: https://doi.org/10.4995/ia.2006.2882
Ministerio de Ambiente y Desarrollo Sostenible (2015). Resolución número 0631. Por la cual se establecen los parámetros y los valores límites máximos permisibles en los vertimientos puntuales a cuerpos de aguas superficiales y a los sistemas de alcantarillado público y se dictan otras disposiciones. https://www.rds.org.co/es/recursos/resolucion-631-de-2015-parametros-vertimientos#:~:text=Por%20la%20cual%20se%20establecen,y%20se%20dictan%20otras%20disposiciones.
Ramos-Velandia, C. C. (2017). Evaluacion de un sistema de tratamiento de lodos activados para la remocion de carga organica en las aguas residuales domesticas generadas en el colegio San Viator [Undergraduate thesis, Universidad de América]. https://hdl.handle.net/20.500.11839/6578
U.S. Environmental Protection Agency. (1974). Process Design Manual for Sulfide Control in Sanitary Sewerage Systems. United States of America: EPA.
Recibido: 9 de febrero de 2021; Aceptado: 7 de febrero de 2022
This research aimed to design a drainage quality hydrodynamic model of the northeastern sector of Tunja (Santa Inés neighborhood), Boyacá, Colombia, which was validated and calibrated by means of the mean square error method, where a comparison between the values observed in the field and the simulation results in the model was determined in order to observe and analyze the effect of precipitation on pollutants before entering the treatment plant.
For this study, rainfall and water flow equipment were installed in the district. A validated and calibrated urban drainage model was created on the SWMM 5.1 software, which allowed a temporal hydrodynamic modeling of the system. In addition, samples were taken to determine the quality of the water arriving at the treatment plant and evaluate the transportation of pollutants, as well as the possibility of hydrogen sulfide generation in the study area.
The results obtained with the model show that all pollutants generate drag in their concentrations. An excess of these values is observed within the total flow, which demonstrates the relevance of accumulation, washing, and transport processes, in light of the pollutant remnants causing the flow to decrease before the pollutant load.
Our sampling indicates that it is necessary to analyze the coliforms present in wastewater discharges, as well as the sulfate, nitrite, and nitrate values in the case of a specific discharge. The entity in charge of discharge control must be then informed.
Universidad Santo Tomas - Tunja
Keywords:sewerage, water quality, modeling, SWMM.
Esta investigación pretendió diseñar un modelo hidrodinámico de calidad de drenaje urbano del sector nororiental de Tunja (barrio Santa Inés), Boyacá, Colombia, calibrado y validado por el método del error cuadrático medio, donde se determinó una comparación entre los valores observados en campo y los resultados de la simulación en el modelo en aras de observar y analizar el efecto de la precipitación sobre los contaminantes antes de ingresar a la planta de tratamiento.
Para el presente estudio se instalaron equipos de precipitación y flujo de agua en el distrito. Se creó un modelo de drenaje urbano validado y calibrado en el programa SWMM 5.1, que permitió obtener la modelización hidrodinámica temporal del sistema. Adicionalmente se tomaron muestras de agua en la red para conocer la calidad de agua que llega a la planta de tratamiento y evaluar el transporte de los contaminantes y la posible generación de sulfuro de hidrógeno en la red de estudio.
Los resultados obtenidos con el modelo muestran que todos los contaminantes generan arrastre de sus concentraciones. Se observa un exceso de estos valores frente al caudal máximo, lo cual demuestra la relevancia de los procesos de acumulación, lavado y transporte, en vista de que el remanente de contaminantes es lo que ocasiona que el caudal disminuya primero que la carga contaminante. Adicional a esto se presenta un análisis en los tramos de concreto en la red para así determinar la posible generación de sulfuro de hidrogeno, como se presenta en la Tabla 6 de este documento.
Nuestro muestreo indica que es necesario analizar los coliformes presentes en los vertimientos de aguas residuales, así como los valores de sulfatos, nitritos y nitratos en el caso un vertimiento puntual. Posteriormente se debe informar a la entidad encargada de control de vertimientos.
Universidad Santo Tomas, Tunja, Boyacá, Colombia.
Palabras clave:alcantarillado, calidad del agua, modelado, SWMM.
Nowadays, the city of Tunja is experiencing a phenomenon of occupation or alteration of the natural channels of rainwater or runoff that lead these flows to the receiving bodies, thus causing an impact on the drainage since the hydrology of the natural basins is not respected, which affects their drainage capacity. In addition, there is considerable and accelerated urban development in the city’s northeastern sector, which has caused concern due to this and the continuous substitution of drainage surfaces that allow infiltration by hard and impermeable areas, causing overflows and floods that affect the citizens (Amaya-Tequia, 2019).
In a city whose aqueduct and sewerage networks are growing, it is necessary to know the operating behavior of the system. Knowing the hydraulic operation of the networks in real time is a vital tool to identify areas of conflict, vulnerabilities, and risks. Thus, it is also possible to propose prevention measures, dictate alternative solutions, optimize the functioning of the urban drainage system, prevent damage and problems in the structures, prioritize resources, meet contractual targets, and improve the citizens’ quality of life (Gerencia de Planeación y Construcciones, 2017).
The over-dimensioning of the elements of a unitary urban sewerage system, even though it allows the transport of flows collected from rainfall runoff, provides a significant potential for sedimentation of wastewater solids during dry periods. This is because runoff speeds are, in many cases, inadequate to keep them in suspension. As a result, during periods without rainfall, the ducts act as reservoirs for said sediments, thus affecting water quality (Seco & Gómez-Valentín, 2011). Earlier, there was a tendency to not consider water quality aspects in the design and operation of drainage systems. This trend remains in the development of most projects (Rodríguez, 2005)
The management of pollution associated with urban runoff is a difficult problem to solve, especially considering the stochastic nature of rainfall and the hydrological regimes of some rivers (Obermann et al., 2009). The natural channels of some regions have flows whose seasonal differences are significant, which makes them more sensitive to discharges from unit systems at times when the water flow provided by the spillways may be of the same order or greater than the flow rate. The assessment of urban rainwater quality is of great importance for the current approach to integrated urban drainage management. The characteristics of wastewater can be determined in different ways depending on the specific purpose. A sampling software for water characterization and quality control involves careful analysis of the samples type, the number of samples, and the parameters to be studied (Romero-Rojas, 2004).
It is important to know where and how the flows are incorporated into the network in order to be able to analyze their hydraulic behavior. For this reason, together with the Veolia Aguas de Tunja S.A. E.S.P. company, a measurement and modeling of the contaminants present in the water was carried regarding BOD (biological oxygen demand), COD (chemical oxygen demand), lead, TSS (total suspended solids), and nitrogen, with the purpose of improving rainfall drainage to avoid economic, environmental, and public health consequences as a result of climate change and its impact on the design flow of urban drainage systems (Veolia Tunja, 2017)
In addition to the above, an analysis of the possible generation of sulfides is also carried out, since one of the main causes for the deterioration of drainage systems is microbiological corrosion. An example of this is Thiobacillus Thiooxidans, which is also known as Acidthiobacillus thiooxidans, a sulfur-oxidizing organism that has been reported in sewerage systems in several countries including Mexico, the United States, Japan, Belgium, and China. The restoration costs of the concrete elements affected by microbiological corrosion represent a great investment in some countries and cities such as Germany, where it accounts for 40% of the US$100 billion invested in wastewater infrastructure; Belgium, where it accounts for 10% of the total expenditure; and Los Angeles, USA, where approximately US$400 million are invested. It is also estimated that the United States spends around US$25 billion annually on the maintenance of sewerage systems (Cortés & Vera, 2019). The corrosion of concrete structures is very common in the world; it is a silent phenomenon that must be addressed because of the large amount of money that is lost in performing maintenance works (Cortés & Vera, 2019).
1. Area characterization and sampling
Topological base drainage network of the Santa Inés district (Amaya-Tequia, 2019)
Once the activities involving inspection, network cadastre, verification of the connectivity of the system, and identification of initial wells were performed, the delimitation of the Santa Inés district, as indicated in Figure 1, was carried out.
The points where the wastewater sampling was conducted are also indicated in Figure 1. These points were chosen given their location, as they are at the beginning, the middle (where a significant part of the system’s water is collected), and at the end of the network.
2. EPA modeling - SWMM (Gironás et al., 2009)
SWMM’s rainwater management model is a dynamic rain runoff simulation model that calculates quantity and quality mainly in urban areas. It operates in a collection of sub-catchment areas that receive precipitation and generate runoff and polluting loads. Routing transports the runoff through a pipe system, tracking its quality within each sub-uptake and the quality of the water in each pipe during a simulation period composed of multiple time steps.
This software shows how to simulate the accumulation and washing of pollutants in an urban basin. The influence of different land uses on pollutant accumulation is considered, and the average concentrations of events such as exponential functions are used to represent the washing process. The quality of surface runoff is an extremely important but very complex problem in the study of wet climate flows and their environmental impact. SWMM provides a flexible set of mathematical functions that can be calibrated to estimate the accumulation of pollutants on the Earth’s surface during dry weather periods, as well as their release from runoff during storm events.
3. Hydrogen sulfide (US Environmental Protection Agency, 1974)
Hydraulic characterization of sewage is necessary to determine the correlation between existing and predicted sulfur concentrations. It requires measuring the speed, flow depth, and slope of the sewer, by means of hydraulic relations. The flow rate and depth are obtained at the time of sampling of sulfides.
Speed can be measured through several methods (e.g., flotation speed). The slope of the culvert can be measured on the ground or taken in the form of constructed plans. The actual hydraulic roughness coefficient (n) can be calculated using the Manning Equation (1):
V: average speed ft/s
A: hydraulic radio (ft)
N: roughness coefficient
Parameters such as biochemical oxygen demand, pH, temperature, dissolved oxygen, and sulfate need to be controlled; their concentrations can be measured by individual or composite laboratory tests.
The daily mean sulfide concentration and pH values can be calculated by applying a correction factor (which is derived from the diurnal sulfide variation graph) to the sulfide concentrations measured in one day. Since speed is one of the determining factors in the generation of sulfides, it is feasible that all three conditions of possible sulfide generation exist in a pipe at some point within a period of 24 h. Therefore, for definition purposes, the three categories that define the sulfide generation characteristics of a pipeline are based on the average speed for the maximum flow period of 6 h.
The quantitative Pomeroy and Parkhurst method for the prediction of sulfides has proven its effectiveness in studies conducted in California, Louisiana, and Texas. It was developed and is applicable only for partially filled trunk sewage culverts when conditions are favorable for sulfide accumulation. Since misleading results may be obtained under other conditions, indicators are evaluated to determine whether there is a possibility of sulfide generation before using this method.
Sulfur generation indicators
The recommended method of analysis is to evaluate the formula Z and A/B curves. If any of the indicators shows that conditions are favorable for the generation of sulfides, the Pomeroy and Parkhurst method is used to determine if there is a real problem. If both indicators show no potential for sulfide generation, no further analysis is required.
The first equation to express the necessary conditions for sulfide generation in gravity nets was developed in 1946. This formula did not deal with sulfide levels, but merely with whether a build-up of sulfide could occur. In 1950, Davy presented a more complete formula that related the Reynolds number, the BOD5 effectiveness, the flow cross-sectional area, and the surface width. This work was later modified by Pomeroy to develop what is known as the Z formula:
Z: defined function
eBOD: effective BOD5, (mg/L)
Q: flow rate (ft3/s)
P: wet perimeter (ft)
b: width (ft)
Z < 5.000
Sulphur is rarely produced
5.000 ≤ Z ≤ 10.000
Secondary condition for sulfide accumulation
Z > 10.000
Sulfur accumulation is common
The sewerage system of the neighborhood is of combined nature since the pipe system collects and transports both waste and rainwater. The pipes are built from concrete (96%) and PVC (4%). The outflow well provides the right conditions for the installation of a flow meter, with a depth to level of 1,58 m. The main collector is built from concrete pipe with a diameter of 30 inches. A simultaneous monitoring of rainfall and flow during 4 months was carried out, and the representativeness of these events was executed by cross-referencing the data recorded in the measuring equipment, verifying the magnitude and time of occurrence between hydrograms and hietograms (Amaya-Tequia, 2019). Within the monitoring period (year 2018), six representative rainfall events were captured, the model of this research was implemented on July 16th because sampling was conducted on that day, with the results shown in Table 32 and obtaining a maximum flow of 161,13 l/s.
Total rain (mm)
Rainfall volume (m3)
< 2 years
Sampling was carried out for 24 h in dry and rainy periods at the established points of the network. During the sampling, in situ parameters such as pH, temperature, and percentage of dissolved oxygen were measured, and the corresponding samples were taken to the laboratory in order to analyze the other parameters. Tables 3, 4, and 5 present the results at each of the sampling points in the study area.
SM 5220 B
SM 4500 S2-F
SM 5220 B
SM 4500 S2-F
SM 5220 B
SM 4500 S2-F
2. SWMM model
To this date, there are several studies on modeling sewerage systems. It has been deduced that the spatial simplification scale influences the results of the SWMM simulation (Marcor & Pedraza, 2012). The total concentration time of the urban watershed of the study area was estimated via the Carter method at 21,60 minutes, considering the weighted average slope of sub-basins and the main length of the flow (Amaya-Tequia, 2019).
In 2012, a study was conducted with the purpose of determining specific quality patterns of the behavior of the sewerage system of a sector in Catalonia, where the unitary network was calibrated and validated with data of rainfall episodes. For ease of access, the SWMM software and the variables adopted in the simulation were tailored when performing the calibration and validation of the proposed model (Seco & Gomez-Valentín, 2011). This model treats each basin as a non-linear reservoir obtained by the continuity and the Manning equation for each sub-basin (Hogue, et al., 1988). Figure 2 presents the continuity of the simulated runoff in the model.
Figure 3 shows the runoff quality in the continuity balance throughout the study area. Input loads are expressed as "Initial Buildup" prior to the start of the simulation, "Surface Buildup" during the dry period, and "Wet Deposition" pollutants in the rain. The output loads include "Infiltration Loss", generated by direct rain; "Surface Runoff", the pollutant load which includes a portion of accumulation; and the continuity report, indicated by "Remaining Buildup".
Figure 4 shows the quality routing. Only runoff loads through the transport system are shown. Dry weather is not shown; it is only supplied by the user through external inputs. Therefore, the three variables represented are "Wet Weather Inflow", "External Outflow", and "Final Stored Mass". The entry for wet weather is the same.
Figure 5 compares the results of the measured flow with respect to the different pollutants established in the model at the point at the outlet of the basin. For all pollutants, the concentrations exceeded the maximum flow, thus demonstrating the initial analysis in which the importance of the accumulation, washing, and transport processes was mentioned, given that it is the remains of the pollutants generated that affect the flow.
Figure 6 shows the behavior of all the pollutants studied in a rainy period. It is evident that BOD5 has a lower concentration than the others and COD is the one with the highest concentration.
2. Hydrogen sulfide prediction
Table 6 shows the results of all sections of the study area. Green indicates that sulfide is rarely produced; yellow indicates a probability that it will occur; and red shows hydrogen sulfide accumulation.
Flow rate (ft3/s))
Oxygen influence on sulphur generation (Roca-Hernández, 2012)
In a sewer system that circulates by gravity, water is aerated. This process can be slow in large collectors due to the gentle slope and flow depth. The rate increases in smaller pipes, and this does not happen in pressure systems, so the oxygen is consumed in a shorter time, thus yielding a higher concentration of sulfides. Oxygen is consumed by microorganisms present in the body of water (in the biofilm of the pipes), and the rate of consumption and oxygen can vary depending on the distance that the wastewater has to travel due to the diffusion in the biologically active film of the pipe wall.
The structure of the biofilm is formed by several layers, with an aerobic and an anaerobic zone. If the former prevails, the conditions will be given for sulfate reduction to occur. The relationship between these zones is delimited by the concentration of organic matter. If the oxygen concentration in the current is close to zero, then not all the sulfide can be oxidized and passed into the current.
The behavior of the physicochemical and microbiological characteristics measured along the sewerage network at the three established points evidences that there is spatial uniformity in the analyzed parameters and that there is little variability in their behavior when the first wash occurs at each point. All this, considering the results indicated in Tables 4, 5, and 6.
The design of the hydrodynamic model of urban drainage quality of the northeastern sector of Tunja, Santa Inés neighborhood, was performed, calibrated, and validated via the SWMM 5.1 software based on data obtained in the field.
The implementation of the model allowed concluding that the sewerage of the city of Tunja, despite being combined, and considering the current dumping regulations in Colombia, namely Resolution 0631 of 2015 (Ministerio de Ambiente y Desarrollo Sostenible, 2015) which indicates the maximum allowable values for any type of dumping, requires prior treatment to avoid contamination in the dry period. As indicated by the rule, the maximum allowed value for suspended solids is 90 mg/L. This limit was exceeded in all three points. During the rainy period, although the value was not exceeded at the two points near the exit, it was very close to the limit: 85 and 86 mg/L, respectively.
It is recommended that treatment alternatives be considered in the activated sludge process, based on the relationship between BOD and COD, since waste can be degraded through a biological process (Ramos-Velandia, 2017).
Given that the Z formula has generally been successful in predicting the occurrence of sulfide problems in gravity sewers, for our network, specifically with 96% of the pipelines, approximately 48% show no likelihood of sulfur generation. For pipelines that are likely to generate hydrogen sulfide, it is recommended to ensure proper aeration in order to avoid gas accumulation.
Todos los textos incluidos en la Revista Tecnura están protegidos por derechos de autor. Conforme a la ley, está prohibida su reproducción por cualquier medio, mecánico o electrónico, sin permiso escrito del Comité Editorial. Los textos completos de los artículos son de acceso abierto, es decir, se pueden leer, descargar, copiar, distribuir, imprimir, buscar o vincular. Las opiniones expresadas en los artículos publicados son las de los autores y no coinciden necesariamente con las del Comité Editorial ni las de la administración de la Facultad.