Simulation of water quality parameters from the treatment of wastewater using stabilization ponds ( Case study : Tehran wastewater stabilization pond )

Today, discharging the wastewaters into the river is the main factor of water pollution and associated with environmental impacts. There are several methods for wastewater treatment. The natural ways like wastewater stabilization pond (WSP) is apply in a lot of countries especially in under developed countries because no needed to high technology. The WSP’s are suitable for wastewater treatment in small societies and some special industries like slaughterhouse, dairy products and meat products and etc. The advantages of this way are high efficiency in treatment of pathogens, toxic and organic materials. These ponds treat the wastewater naturally by alga and bacterial reactions. Ecological condition and wastewater properties affect the efficiency of WSP. Furthermore the efficiency of WSPs depends to several factors like the temperature, light, alga and bacterial reactions and Etc and must be analyzed before construction. In this regard the quality parameters in Tehran WSP simulated with CE-QUAL-W2. The model has been calibrated and evaluated on the basis of collected data from Tehran WSP and the ability of the model analyzed to optimum the efficiency. Furthermore the concentration trend of parameters analyzes during the day in different hours to determine the effects of solar radiation and the others parameters on them. The results shows BOD concentration increase at noon and decrease at night that shows it depends to solar radiation. Alga also is a function of solar radiation and the concentration increase in noon and decrease in afternoon. The trend of phosphate and nitrate is the same and since these two parameters are used for growth process of algae so in noon decrease and in night increased. Because at night the growth process of algae is limited and no need for phosphate and nitrate exist so the concentration increased. The results of modeling with CE-QUAL-W2 show good agreement between the experiments and the simulation data and CE-QUAL-W2 is a suitable model for analyzing the quality parameters in WSP and can be used for design of new stabilization ponds or for improving the effluent quality of existing ponds. @JASEM

Today, pollution of water resources with wastewater is one of the main environment impacts of the wastewater.The human activities have different effects on environment and to minimum the impacts the new methods and management works needs in all programs (Khodadadi et al. 2005).It is estimated that 1 cubic meter domestic wastewater can pollute about 40 to 60 cubic meter fresh water (Farzad kia et al. 2003).Furthermore the wastewater stabilization ponds can provide satisfactory sewage treatment especially in rural areas (Pirooz et al. 2010).High efficiency in removing the pathogens and organic materials are the advantages of WSP.Furthermore, wastewater treatment points are not very sensitive to shock loading (Mirbagheri et al. 2004).Also the experiments shows numerical methods are a suitable ways to simulate the pollution parameters (Abaspour et al. 1996).

MATERIALS AND METHODOLOGY
Case Study: The study area is situated in the south of Shahr rey City within Tehran province.The average elevation is 980m in the area (Figure 1).

Figure 1-The wastewater plan
Meteorological conditions:The efficiency of wastewater stabilization ponds is in relate with meteorological conditions.The study area situate in semi-arid district with hot summer and cold winter.The maximum temperature is more than 42 C and the minimum temperature is less than 10 C. The mean temperature is 10 C. the relative humidity is about 15% in summer and 75% in winter.The wind direction is from west and southwest to east in whole of the year and the average is 20 Km per hour (TPWW, 1998).
Background on numerical models: CE-QUAL-W2 is a laterally averaged, 2-D hydrodynamic and transport model designed for long-term, time-varying water quality simulations of lakes, reservoirs, and estuaries (Cole at al. 1995).Formulations of the fundamental equations, the structure of the computations, and summaries of applications are described in (Buchak et al. 1988 andCole 1994).CE-QUAL-W2 model is based on the assumption that flow and transport phenomena in water body with a distinct flow direction and a regular bathymetry can be computed from the laterally integrated Navier-Stokes equations (Bloss et al . 1988 andEdinger et al. 1978).CE-QUAL-W2 uses the laterally averaged equations of continuity, momentum and transport.The main formulations involve a vertically varying, longitudinal momentum balance, vertical momentum in the form of the hydrostatic approximation and local continuity, and longitudinal and vertical transport of the different constituents.The continuity equation used in CE-QUAL-W2 is of the form ( ) where x and z the horizontal and vertical coordinates, b the width of the water body, u the laterally averaged longitudinal velocity, w the laterally averaged vertical velocity, and q the lateral inflow or outflow per unit area.In the CE-QUAL-W2 model, variations in width with depth are taken into consideration.Eqs.(2a) and (2b) represent the horizontal momentum equation in the CE-QUAL-W2 model, respectively.
where z N is the vertical eddy viscosity, p the pressure, ρ the density of water, t the time, x N the longitudinal momentum dispersion coefficient and ) is the shear stress per unit mass resulting from the vertical gradient of the horizontal velocity, u The second term on the RHS of Eq. (2b) of CE-QUAL-W2 is the horizontal dispersion of momentum.In the CE-QUAL-W2 model, the effect of stratification on vertical eddy diffusivity is not included.The longitudinal pressure gradient used in model results from the hydrostatic pressure assumption (vertical momentum equation): Where g is gravitational acceleration and η the free water surface elevation.The transport equation in CE-QUAL-W2 is given by Where where sal c is the salinity.The density may also be adjusted for the effect of suspended solids Where ss c is the suspended solids concentration  Numerical solution:The CE-QUAL-W2 model uses a space-staggered, finite difference grid, but the CFL criterion has been avoided by the following numerical procedure to solve the unknowns in Eqs.
(1), ( 2b), ( 3), ( 4) and ( 5).At each time step, the numerical procedure is to first compute water surface elevations.With the new surface elevations, horizontal and then vertical velocities are determined so that the new constituent concentrations can be computed.Using the updated horizontal and vertical velocities, the water surface elevation (Eq.( 5)) can be solved for η simultaneously.The solution for η is thus spatially implicit and eliminates the surface gravity wave speed criterion, which can seriously limit time steps in deep water bodies (Cole et al. 1995).A major problem with upwind differencing is the introduction of numerical diffusion where there is longitudinal advection.In many cases, numerical diffusion can overwhelm physical diffusion and produce inaccurate results when strong gradients are present.In the CE-QUAL-W2 model, numerical diffusion has been reduced by implementing QUICKEST, an explicit, third-order accurate horizontal/vertical transport scheme (Fletcher et al. 1991).

Simulation and modeling
Qualitative analyses of wastewater stabilization ponds with CE-QUAL-W2: Due to analyze and simulate the quality parameters of WSP with CE-QUAL-W2 some parameters include: PH, BOD, TDS, Alga, Nitrate, Organic Nitrogen, Ammonia, Temperature and phosphate measured in 5 sections of the pond in a 15 days period (88/2/1 to 88/2/15) with 4 hours step (4, 10, 14 and 18).The distance between the samples is 28m.The pond divided to 6 layers in vertical.Two layers with o.5m from the surface and four 1m layers.First and last sections and layers take as the boundary condition.
Calibration and validation the model: For calibration and validation the data divided to two groups.65 percent of data used to calibration and 35 percent to validation.

Calibration of CE-QUAL-W2:
In this part at first sensibility analyses have been done and then the simulation data for calibration obtained.The mean of absolute error and performance coefficient compute to determine the accuracy of the model and are shown in table (1).The results show that Ammonia, TSS and BOD have the minimum mean-squared error and the best simulation by CE-QUAL-W2 model.

Validation of CE-QUAL-W2
After calibration, the validation of the model has been done with 5 days data (30 percent of data).As it can be seen in the figures BOD concentration increase at noon and decrease at night that shows it depends to solar radiation.Alga also is a function of solar radiation and the concentration increase in noon and decrease in afternoon.The trend of phosphate and nitrate is the same and since these two parameters are used for growth process of algae so in noon decrease and in night increased.Because at night the growth process of algae is limited and no need for phosphate and nitrate exist so the concentration increased.After validation the data, with analytical analysis the mean of absolute error and performance coefficient compute in this stage to determine the accuracy of the model and are shown in table (2).The data from validation compared with the field data and the results in the outlet (the fifth section) in 88/2/15 are shown in figures ( 6) to (9).The results show the good capability of the model to simulation the quality parameters of WSPs.

Conclusion:
The results of modeling with CE-QUAL-W2 show good agreement between the experiments and the simulation data.The CE-QUAL-W2 is a suitable model for analyzing the quality parameters in WSP and to optimum the efficiency of the operation.It is recommended to simulate the quality parameters in WSP to optimum the efficiency.As it can be seen in the figures BOD concentration increase at noon and decrease at night that shows it depends to solar radiation.Alga also is a function of solar radiation and the concentration increase in noon and decrease in afternoon.The trend of phosphate and nitrate is the same and since these two parameters are used for growth process of algae so in noon decrease and in night increased.Because at night the growth process of algae is limited and no need for phosphate and nitrate exist so the concentration increased.The results of modeling with CE-QUAL-W2 show good agreement between the experiments and the simulation data and CE-QUAL-W2 is a suitable model for analyzing the quality parameters in WSP and can be used for design of new stabilization ponds or for improving the effluent quality of existing ponds.

xDb
is the longitudinal temperature or constituent dispersion coefficient and n H the source or sink term for the particular constituent.The fourth term in Eq. (4) of the CE-QUAL-W2 model represents the horizontal diffusion term for transport of a constituent.In CE-QUAL-W2, there is an additional equation for free water surface elevation.is the width, which varies in time and space, h the total water depth and Q the lateral boundary inflow or outflow in the specific gravity of suspended solids, which if taken to be 2.65 yields.
based on equations for effects of latent and evaporative heat exchange and shortand long-wave radiation.

Figure ( 2
Figure (2) simulation of TSS and BOD in sections 1 to 5 at fifteenth day (mg/l)