Skip to content
BY 4.0 license Open Access Published by De Gruyter Open Access March 4, 2024

Evaluation of the operational performance of Karbala waste water treatment plant under variable flow using GPS-X model

  • Noor Abdulhussein Khudhair EMAIL logo , Basim K. Nile and Jabbar Hammoud Al-Baidani
From the journal Open Engineering

Abstract

As Karbala City is a religious tourism destination, millions of tourists visit the city annually, so there is a high fluctuation of flow in the wastewater between the plants during the year. The aim of this study is to evaluate the performance of the wastewater treatment plant (WWTP) in Karbala, Iraq, in removing pollutants for different flow rates in five scenarios using the GPS-X model. The most important phase in modeling, which greatly impacts simulation accuracy, is characterizing the influent composition to meet the mass balance. As a result, the influent wastewater was initially described and thoroughly examined. The model has been calibrated, followed by the collected data’s validation. The sensitivity of different stoichiometric and kinetic factors in the GPS-X was examined and screened to calibrate the model. To demonstrate the consistency between the simulated and measured data, the route mean square error was used in this instance. The result showed that the Karbala WWTP complies with Iraqi environmental regulations for water discharged to surface water or water for other uses and has an appropriate efficiency of wastewater treatment even if the flow entering the plant reaches the peak flow rate of 180,000/day. An improvement in orthophosphate removal efficiency was observed as the flow rate increased because of the contact time in the anaerobic basins, which gradually decreased as the flow increased until it reached the appropriate time for PO 4 3 removal. The outcomes of the present study provide an impression to the operators of the treatment plant of the impact of fluctuating flow on the treatment plant. The developed model can also be used for future studies.

1 Introduction

The volume of wastewater entering the treatment facilities has increased due to urbanization and industrial growth. Wastewater contains substances that are harmful to both people and the environment. As a result, dumping untreated wastewater in natural water bodies creates serious sustainability issues, like endangering aquatic and terrestrial life and raising the expense of treating the polluted water [1]. Consistency with discharge characteristic criteria, harm inhibition on the environment and human health, feasible water reuse, and recovery concerns connected to energy and materials (e.g., nutrients) make proper wastewater management a significant global task [2]. To adhere to environmental regulations, enterprises that produce pollutants must effectively treat their wastewater before releasing it to maintain the long-term viability of their water resources [3]. Activated sludge-based systems are currently the most popular choice since they are the most cost-effective and efficient of the several biological wastewater treatment methods used today. Severe weather might bring either dry or rainy spells, and the population expansion will create new difficulties for operating a wastewater treatment facility. The wastewater treatment process (WWTP) can be simulated under various intake circumstances to prepare for these issues. In harsh climatic circumstances, the feed stream flow rates may exceed the intended maximum values, and the compositions may fluctuate greatly [4].

It is crucial to employ simulation to make the most of the model’s predictive skills, which may be used to quickly and effectively weed out the optimal design option and cut down on the time and expense of laboratory testing. Modeling is the process of simplification of actual representation. Mathematical operations and equations used time-dependent variables and parameters to define the model [5]. Recent times have seen a rise in dynamic modeling and simulation in wastewater treatment. Many models are created to enhance activated sludge. Several software programs, such as SIMBA, GPS-X, AQUASIM, BioWin, STOAT, FOR, and WEST that promote dynamic modeling of wastewater treatment facilities use these models. Tools for simulation and modeling are used to evaluate process procedures, optimize designs, and analyze costs [6]. The calibration and validation steps are crucial for these models to give accurate results. To precisely manage and calibrate the plant’s operation, the stoichiometry of the reactions inside the reactor controls the materials produced and consumed. To investigate the microbial transformation process in removing organic matter and nutrients, kinetics and stoichiometric measurements in the titration and verification procedures were considered necessary [7].

The GPS-X model has been used in a wide range of research. Faris et al. [5] investigated the medium volume, dissolve oxygen content, and sludge return to the moving bed biofilm reactor by using sensitivity analysis. The results demonstrated that the nitrification process is impacted when the dissolved oxygen content exceeds 3 mg/L in the ammonia-rich side streams generated from the rejected water of the anaerobic digester. Alwardy et al. [8] evaluated Al-muamirah WWTP in Hilla City, the results revealed that the facility has an acceptable level of wastewater treatment efficiency and produces water that complies with Iraqi environmental standards for water discharged to surface water or water for other uses. Kobeyeb et al. [9] Presented a study to design and model an on-site greywater treatment system for a hotel building in Los Angeles, California. Many greywater treatment plant options were rigorously examined and modeled using the GPS-X software. The MBR plant was the best option for considering location and standards and was advised for use in hotel structures. Hammed et al. [10] used the GPS-X model to investigate how to lower nutrient concentrations in a full-scale plant. The findings demonstrated that the lack of rbCOD in significant quantities to aid in reducing nitrates rendered having a proportion of internal recycle (IR) of 3% worthless. Cao et al. [7] introduced an innovative way using GPS-X and response surface methodology to improve the removal of total nitrogen (TN) in WWTPs. 61 parameters’ sensitivities were checked and examined. The findings demonstrated that the denitrification rate was significantly impacted by the DO concentration that diffused into various biological compartments. SRT and TN elimination go hand in hand. Key parameter relevance and optimization orders were examined.

Due to the city of Karbala’s popularity as a destination for religious travelers, there is a significant variation in wastewater flow into the plant. The present study aims to study the effect of variation in flow on the operational efficiency of the Karbala WWTP by using five scenarios with different flow rates (minimum flow to peak flow) in GPS-X software to model the plant. Additionally, the optimal operation of the WWTP was suggested for each scenario of different flow rates depending on the model results after calibration and validation of the model.

2 Materials and methods

2.1 Karbala WWTP

The WWTP in the city of Karbala serves 2.5 million people. The conventional activated sludge method is vital to Karbala’s WWTPs. Four sewage treatment facilities utilizing the conventional activated sludge system of type A2/O are used in an integrated wastewater treatment project in Karbala City. Each facility has a daily discharge capacity of 100,000 m3. This plant’s geographical coordinates are 32.525590° North and 44.074909° East. Figure 1 shows a satellite.00000000000e image of the Karbala WWTP.

Figure 1 
                  Google Earth satellite view of the Karbala WWTP.
Figure 1

Google Earth satellite view of the Karbala WWTP.

The plant involves five steps of treatment. The first step is preliminary (physical) treatment, which includes four units, which are coarse and fine screens used for removing coarse and fine solid materials. A grit and oil removal chamber is used to remove sand and oil by a physical machine. The next step is primary treatment, carried out by four primary sedimentation basins. After that, there is the secondary treatment, which involves four tanks. There are two anaerobic reactors with a volume of 8,736 m3 for both of them, two anoxic reactors with a volume of 14,112 m3, eight aeration reactor tanks with a volume of 54,054 m3, and eight sedimentation tanks with a surface area of 6,432 m2. The next step is tertiary treatment, which includes chemical disinfection through the chlorination tank with a surface area of 3,000 m3. The last stage is sludge treatment, which involves two gravity thickeners with a surface area of 400 m2, three mechanical thickeners with a surface area of 60 m2, four anaerobic digesters, and 60 cells of drying beds with a surface area of 50,000 m2.

2.2 Collection of samples

The purpose of the present study was to study the effect of variable flow on the operational performance of Karbala WWTP. Influent sewage has been analyzed for its properties to measure contaminants’ concentrations prior to treatment. Moreover, concentrations of effluent contaminants following treatment have been evaluated and compared with the Iraqi standard limits. Four tests for each pollutant were conducted monthly for 1 year, from January 2022 to December 2022. According to the Standard Procedures for the Examination of Water and Wastewater, the tests were conducted in the laboratory of the Karbala Sewage Department (APHA 2017). The biological oxygen demand (BOD5), chemical oxygen demand (COD), total suspended solids (TSS), ammonia-nitrogen ( NH 3 + N ), orthophosphate ( PO 4 3 ), nitrate ( NO 3 -N ), and nitrite ( NO 2 -N ) are measured. Table 1 displays the concentrations of contaminants in the influent and effluent wastewater.

Table 1

Concentrations of influent and effluent sewage of Karbala WWTP

Parameter Influent, mg/l Effluent, mg/l Iraqi standard
COD 250 23 100
BOD5 115 4 40
TSS 140 11 60
NO3 0 16 8
NO2 0 0.2
NH 3 + 28 0.5 10
PO4 3− 4.5 2.5 3
H2S 30 0.3 3
DO 0 3.5
SO4 850 860 600

Note: all parameter units are in (mg/l).

2.3 Modeling of Karbala WWTP using GPS-X model

Mathematical modeling has emerged as an essential instrument for sustainable wastewater management, particularly in simulating intricate procedures inherent in activated sludge processes. Biological systems are very complicated. Temperature, flow fluctuation, the quantity of the wastewater pollutant, and the operational circumstances are just a few of the many variables significantly affecting treatment efficiency. The modeling technology company Hydromantis’ GPS-X program, version 8.0.1 (education license), was utilized in the present study. GPS-X is a multipurpose modeling platform designed to simulate WWTPs. This allows for the interactive and dynamic examination of the complicated relationships between the plant’s many unit processes [10]. It was constructed with integrated biological WWTP and many other processes that involve physical and chemical reactions. Data were obtained from the WWTP of Karbala and fed into the GPS-X program. A comprehensive carbon, nitrogen, and phosphorus pH (mantis2lib) library was selected. The dimensions of each unit and the input concentration were entered into the model. The assumptions used to create the model in this study are: the pH is stable and close to neutral, and the biological process takes place at a temperature of 20°C (Figure 2).

Figure 2 
                  Schematic diagram of Karbala WWTP.
Figure 2

Schematic diagram of Karbala WWTP.

To establish the best match between simulation results and actual plant effluent, the model was calibrated by adjusting stoichiometric, kinetic, and other important parameters to fit the simulation. After that, mathematical validation was used to check the model, and the following equations illustrate the results for root mean square error (RMSE) and correlation coefficient (R). The reasonable limits of these statistical criteria are 1 ≥ R > 0.8 and 0 ≤ RMSE < 1.5.

(1) R = ( C o C o ¯ ) ( C p C p ¯ ) ¯ σ c o σ c p ,

(2) RMSE = ( C o C p ) 2 ¯ C o ¯ C p ¯ ,

where C o is the actual data, C p is the modeled data, C o is the average of actual data, C p is the average of modeled data, and σ is the standard deviation over the dataset. Many simulation runs were carried out under different flow conditions to learn how varying flows may affect the plant’s efficiency. The steps listed above, illustrated in Figure 3, are used to simulate the modeling of the plant [11].

Figure 3 
                  Flow chart of the present study.
Figure 3

Flow chart of the present study.

3 Results and discussion

The average monthly flow entering the WWTP is 60,000 m3/day, under the design value of the plant, which has a value of 100,000 m3/day. There is little variation in flow during the year; the flow reaches its peak value during 2 months of the year because of the visitors coming to the city to revive religious rituals. The average monthly influent concentrations of TSS, COD, BOD5, NH4 +, and PO 4 3 are 140, 250, 115, 28, and 4.5, and there is a variation in the parameters mentioned earlier during the year. This fluctuation is due to seasonal variations in wastewater concentrations and industrial wastewater from industrial regions that enter the Karbala wastewater treatment system through the municipal wastewater system.

3.1 Modeling of Karbala WWTP by using GPS-X program

The Karbala WWTP was modeled using GPS-X by utilizing five different scenarios.

3.2 Model calibration

In order to evaluate the effectiveness of the model, it is recommended to perform calibration procedures utilizing organic or default data. The model was verified using effluent COD, effluent TSS, BOD, orthophosphate, ammonia nitrogen, nitrite, and nitrate as indicators. This is compared with the calibration of the Albarrakiya WWTP by GPS-X software done by ref. [12], which uses the data of COD, BOD, and TSS to calibrate the model. Calibration is necessary due to the system’s limited understanding of real physical, chemical, and biological processes and the desire to balance model simplicity with model accuracy. The model ran with the default data inside; hence, the outcomes were wildly inaccurate [13]. The model was calibrated with data on the average monthly pollutant concentration for 1 year in 2022. Compared with the study by Alwardy et al. [8], which uses data from 2 months to calibrate the Al-muamirah WWTP model by GPS-X extended from 1/2021 to 2/2021. Several critical parameters that influenced the alteration of outcomes were adjusted. Characterizing the influent wastewater is seen as the most important phase in the modeling process, requiring rigorous investigation. The laboratory data pertain to the influent and effluent data for an influent flow of 60,000 m3/day. This represents the second scenario, which involves adjustments to achieve minimal differences between the simulated and observed pollutants in the effluent wastewater. The model used for influent characterization is the carbon, nitrogen, phosphorus, and pH (mantis3lib) library. In this regard, the GPS-X already contains the default settings for the COD fractions; however, these default values have been modified to ensure better model calibration [14]. In Table 2, the significant influent raw wastewater portions are mentioned.

Table 2

Default and calibrated values of influent stoichiometry composition

Influent stoichiometry composition
Parameter Symbol Default value Calibrated value
Soluble inert fraction of the total COD frsi 0.05 (raw) 0.08 (primary) 0.065
Readily biodegradable fraction of the total COD frss 0.2 (raw) 0.32 (primary) 0.09
Particulate inert fraction of the total COD frxi 0.13 (raw) 0.12 (primary) 0.168

In this instance, the effluent TSS readings were found to be more than 30 mg/L, which was far from the result of the actual influent flow of 60,000 m3/day. As a result, the GPS-X model’s sensitive parameters with regard to the secondary sedimentation basin were changed. Some delicate parameters have been adjusted, including the maximum Vesilind settling velocity, which has gone from 410 to 981.95 m/day; the maximum settling velocity, which has been adjusted from 274 to 356 m/day; and the feeding point from the bottom, which has been changed from the default 1 m to the real 4.3 m. Along with other operational characteristics of a different dimension, these three parameters were the most sensitive.

The effluent COD reading was higher than the result of the plant settling correction factor for xii, and the settling correction factor for xbai was adjusted from their default values 1, 1 to 0.07, 0.26. Acetate fraction of total COD was changed from the default value of 0 to 0.102 to adjust the value of PO4 3−, whose reading was 5; after the calibration, the reading was 3.122. The default nitrate values were shown to be substantially lower than the actual quantities, but the adjustment resulted in fitting the projected nitrate concentration with reality after modifying some sensitive parameters on the kinetics of ammonia and nitrates, including (ammonium fraction of soluble TKN) and (Aerobic heterotrophic yield on soluble substrate). The effluent stoichiometry composition is shown in Table 3.

Table 3

Effluent stoichiometry composition

Effluent stoichiometry composition
Calibrated fraction Default Calibrated
Maximum settling velocity 274 365
Maximum Vesilind settling velocity 410 981.95
Settling correction factor for xii 1 0.07
Settling correction factor for xbai 1 0.26
Acetate fraction of total COD 0 0.102
Aerobic heterotrophic yield on soluble substrate 0.6666 0.36963

Figure 4 shows the technique to handle the calibration process. After calibration, the model is ready to be applied to the scenarios.

Figure 4 
                  Illustration diagram of calibration process.
Figure 4

Illustration diagram of calibration process.

3.3 Statistical analysis

The next task after the calibration of the model for a whole year in 2022 is to examine the calibration process using statistical analysis. Figure 5 shows the actual and simulated concentrations of each pollutant. The statistical equations used to examine the calibration process are the correlation coefficient (R) and RMSE.

Figure 5 
                  The calibration of the actual and the simulated values.
Figure 5

The calibration of the actual and the simulated values.

The correlation coefficient for each actual and simulated parameter was above 0.8, which means that there is a positive linear relationship so that as one variable increases, the other variable also increases proportionally. The RMSE for each parameter was close to zero, and that indicates that the model has more accurate predictions and matches the data well. Table 4 shows the correlation and the RMSE for each parameter of actual and simulated values for average monthly data for a period of 12 months during the year 2022.

Table 4

Value of R and RMSE of calibration process

Parameter R-value RMSE
TSS 0.92 0.011
BOD5 0.87 0.082
COD 0.85 0.021
NH 3 + -N 0.89 0.022
NO2–N 0.83 0.027
NO3–N 0.86 0.138
PO 4 3 0.81 0.011

3.4 Sensitivity analysis

In order to assess the potential effects of variation in flow on the output variables of BOD, COD, TSS, PO 4 3 , NH 3 + -N , NO 3 -N , and NO 2 -N of effluent, a basic sensitivity analysis was carried out as part of this study. This research is advantageous because it improves model prediction and reduces the number of model parameters that require calibration as per EPA standards [15].

Figures 6 and 7 demonstrate how the influent discharge affects the plant’s production. The hydraulic and organic loads on the plant’s reactors increase with the increase in the discharge. This has a negative effect on the treatment processes and increases the concentrations of BOD, COD, and TSS in the treated effluent wastewater. The BOD, COD, and TSS concentrations did not exceed the Iraqi standards even while the Karbala wastewater treatment facilities reached 180,000 m3/day, the plant’s peak flow. Reducing the levels of pollutants in the treated wastewater might be possible if the mass balance could be immediately altered in response to the quantity and quality realities. It has been demonstrated that phosphate removal efficiency increases as discharge increases. For the removal of phosphorus during rbCOD fermentation, acetate is required. For rbCOD fermentation, retention times of 0.25–1.0 h are sufficient [16].

Figure 6 
                  Effect of flow rate on effluent COD, BOD, and TSS concentrations.
Figure 6

Effect of flow rate on effluent COD, BOD, and TSS concentrations.

Figure 7 
                  Effect of flow rate on effluent concentrations of ortho-phosphate, ammonia nitrogen, nitrite, and nitrate.
Figure 7

Effect of flow rate on effluent concentrations of ortho-phosphate, ammonia nitrogen, nitrite, and nitrate.

3.5 Applied scenarios

Upon completion of the calibration processes of the model, the model is ready to be applied in the required scenario and estimate the plant’s response to the variable flow. Five scenarios were used that depended on variable flow. In the first scenario, the number of primary sedimentation basins was decreased from 4 to 2, which resulted in a reduction in energy consumption by 50% as well as maintenance and operating expenses. The number of secondary sedimentation basins was reduced from 8 to 3, which helped reduce energy consumption and maintenance costs by about 38%. In the second scenario, the number of primary and secondary clarifiers was reduced to half, which reduced 50% of the energy and cost of operation. While in the third scenario, the number of secondary clarifiers was reduced from 8 to 6. All the primary and secondary clarifiers were used in the fourth and fifth scenarios.

Table 5 demonstrates that the COD and BOD removal efficiencies for each scenario’s examined under different flows were over 90%, thus demonstrating that the organic matter was removed highly effectively. Microorganisms, particularly heterotrophic bacteria, assisted in removing organic materials once the appropriate dissolved oxygen concentrations were present and there was enough mixing inside the reactor to cause these compounds to break down and turn into fixed substances [17]. When the flow increases, COD and BOD readings increase; this increase could result from a shorter detention time period. The time wastewater is allowed to remain in the treatment system is known as the detention time. Shorter detention times occur when the influent flow is increased because the wastewater is treated more quickly. As a result, the microorganisms in charge of breaking down organic waste have less time to work, which raises the COD and BOD levels in the effluent, as indicated in Figures 8 and 9.

Table 5

Effect of fluctuation of flow on effluent COD and BOD concentrations

Flow rate, m3/day 40,000 60,000 100,000 140,000 180,000
COD mg/L 24.36 26.97 31.22 35.65 40.86
COD, removal efficiency 90% 89% 87% 85% 83%
BOD mg/L 3.807 4.155 5.434 6.748 7.829
BOD, removal efficiency 96% 96% 95% 94% 93%
Figure 8 
                  The variation in the effluent concentrations of COD with various flow rates.
Figure 8

The variation in the effluent concentrations of COD with various flow rates.

Figure 9 
                  The variation in the effluent concentrations of BOD with various flow rates.
Figure 9

The variation in the effluent concentrations of BOD with various flow rates.

The increasing flow rate through the settling tanks leads to decreased settling efficiency. Table 6 demonstrates that TSS values also increase with the increase in flow rate. As a result, the amount of suspended particles in the effluent may increase, increasing the TSS values, which can be attributed to shorter retention times. Higher TSS levels in the effluent might result from shorter contact times between the wastewater and the treatment processes, which can prevent suspended particles from settling. For all scenarios, the secondary sedimentation basins achieved a removal efficiency of greater than 90% of the suspended solids. The removal of nutrients was impacted by the very low concentration of rbCOD in the wastewater of the city of Karbala [18]. Figure 10 shows the increase in TSS concentrations with the increase in the plant flow rates.

Table 6

Effect of fluctuation of flow on effluent TSS concentration

Flow rate, m3/day 40,000 60,000 100,000 140,000 180,000
TSS, mg/L 7.9 11.6 17.23 22.05 27.37
TSS removal efficiency 94% 91% 87% 84% 80%
Figure 10 
                  The variation in the effluent concentrations of TSS with the variation in flow rates.
Figure 10

The variation in the effluent concentrations of TSS with the variation in flow rates.

The NH 3 + -N values change for each scenario, as shown in Figure 11. The value of NH 3 + -N reaches 3.835 when the flow is 180,000 m3/day because the HRT decreases due to the increase in flow rate. With less time available for treatment, the nitrification process, which converts ammonium ( NH 3 + -N ) to nitrate ( NO 3 -N ), may be incomplete. As a result, higher ammonium levels can be observed in the effluent. Another reason for increasing the concentrations of NH 3 + -N is that when the influent flow increases, the contact time between the wastewater and the air or aeration mechanism tends to be reduced. This can limit the stripping of ammonia, leading to higher NH 3 + -N concentrations in the effluent. As the effluent concentration of ammonia nitrogen increases with the increase in the flow, the removal efficiency decreases. It was 98% when the flow was at its minimum value, which reduced the efficiency to 86% when it increased to 180,000 m3/day, as shown in Table 7. The nitrification process was functioning effectively at the Karbala WWTP, with ammonia being oxidized by over 95% for minimum and actual flow. The denitrification process at the Karbala WWTP is experiencing suboptimal performance due to an inadequate supply of appropriate organic substrate.

Figure 11 
                  The variation in effluent concentrations of 
                        
                           
                           
                              
                                 
                                    NH
                                 
                                 
                                    3
                                 
                                 
                                    +
                                 
                              
                              -N
                           
                           {\text{NH}}_{3}^{+}\text{-N}
                        
                      with different values of flow rates.
Figure 11

The variation in effluent concentrations of NH 3 + -N with different values of flow rates.

Table 7

Effect of fluctuation of flow on effluent NH 3 + -N concentrations

Flow rate, m3/day 40,000 60,000 100,000 140,000 180,000
NH 3 + -N , mg/L 0.329 0.41 0.7273 1.491 3.835
NH 3 + -N removal efficiency 98% 98% 97% 94% 86%

The potential impact of nitrogen and phosphorous removal on the high levels of dissolved oxygen discharged through IR and return activated sludge flow rates should also be considered [19]. The primary source of nutrition for phosphorous bacteria is rbCOD. Additionally, it serves as a source of feeding and functions as an electron donor in nitrate reduction. The competition for removing phosphates and nitrates by this compound may result in an uneven distribution of these nutrients. This parameter is highly responsive to the depletion of nutrients [20]. Table 8 shows the removal efficiency of phosphorous of the Karbala WWTP is considered low for minimum and actual flow, 31 and 29%, and this is in agreement with the study by Hammed et al. [10]. An increase in rbCOD concentration correlates with improved efficiency in phosphate removal and enhanced nitrification and denitrification processes. Due to the low concentrations present, it is necessary to incorporate an external carbon source to enhance the nutrient removal processes at the Karbala WWTP. The introduction or increase in rbCOD levels has the potential to impact the BOD and COD concentrations present in the plant’s effluent.

Table 8

Effect of fluctuation of flow rates on effluent PO 4 3 concentrations

Flowrate m3/day 40,000 60,000 100,000 140,000 180,000
PO 4 3 mg/L 3.03 3.122 1.543 0.932 0.66
PO 4 3 removal efficiency 31% 29% 67% 81% 86%

It was noticed that when the discharge increases, the nitrates and nitrites are only slightly influenced. Figures 12 and 13 show that, unlike NO 2 -N , the concentration of NO 3 -N decreases when the flow increases (Figure 14).

Figure 12 
                  The variation in the effluent concentrations of orthophosphate with different values of flow rates.
Figure 12

The variation in the effluent concentrations of orthophosphate with different values of flow rates.

Figure 13 
                  The variation in the effluent concentrations of nitrate with different values of flow rates.
Figure 13

The variation in the effluent concentrations of nitrate with different values of flow rates.

Figure 14 
                  The variation in the effluent concentrations of nitrite with different values flow rates.
Figure 14

The variation in the effluent concentrations of nitrite with different values flow rates.

4 Conclusion

The GPS-X program was successfully used in the present study to model the largest WWTP in Karbala, Iraq. The default settings of some model parameters, such as the kinetic and stoichiometric parameters and the influent wastewater characterization, have been calibrated. The developed model can be used for future studies because the model’s accuracy was acceptable, and the simulation result was close to the actual parameter concentrations of the plant. The results of the five scenarios with different flows show that the removal efficiency of the plant remains acceptable according to Iraqi standards when the influent flow is the minimum and actual flow for all pollutants except PO4 3−, which has a removal efficiency of 31% for minimum flow and 29% for basic flow. The removal efficiency is acceptable even if the influent flow reaches 180,000 m3/day, which is above the peak flow of the plant. In the future, it is suggested to decrease the number of primary clarifiers to two and secondary clarifiers to three when the influent flow is 40,000 m3/day or less and use two primary clarifiers and four secondary clarifiers when the influent flow to the treatment plant reaches 60,000 m3/day in summer or dry days of the year.

  1. Funding information: We declare that the manuscript was done depending on the personal effort of the author, and there is no funding effort from any side or organization.

  2. Conflict of interest: The authors state no conflict of interest.

  3. Data availability statement: Most datasets generated and analyzed in this study are in this submitted manuscript. The other datasets are available on reasonable request from the corresponding author with the attached information.

References

[1] Aghalari Z, Dahms HU, Sillanpää M, Sosa-Hernandez JE, Parra-Saldívar R. Effectiveness of wastewater treatment systems in removing microbial agents: a systematic review. Globalization Health. 2020;16:1–11.10.1186/s12992-020-0546-ySearch in Google Scholar PubMed PubMed Central

[2] Jafarinejad S, Jiang SC. Current technologies and future directions for treating petroleum refineries and petrochemical plants (PRPP) wastewaters. J Environ Chem Eng. 2019;7(5):103326.10.1016/j.jece.2019.103326Search in Google Scholar

[3] Ismail MS, Yahya MD, Auta M, Obayomi KS. Adsorption behavior of heavy metal ions by hybrid inulin-TEOS for water treatment. Civ Eng J. 2022;8(9):1787–98.10.28991/CEJ-2022-08-09-03Search in Google Scholar

[4] Lahdenperä E, Koiranen T. Effects of variable feeds in wastewater plants by means of process simulation. 11th ICARD, IMWA, MWD Conference – Risk to Opportunity; 2018.Search in Google Scholar

[5] Faris AM, Zwain HM, Hosseinzadeh M, Majdi HS, Siadatmousavi SM. Start-up and operation of novel EN-MBBR system for sidestreams treatment and sensitivity analysis modeling using GPS-X simulation. Alex Eng J. 2022;61(12):10805–18.10.1016/j.aej.2022.04.026Search in Google Scholar

[6] Hernandez-Sancho F, Sala-Garrido R. Cost modelling in waste water treatment processes: an empirical analysis for Spain. Dangerous pollutants (xenobiotics) in urban water cycle. Netherlands: Springer; 2008.Search in Google Scholar

[7] Cao J, Yang E, Xu C, Zhang T, Xu R, Fu B, et al. Model-based strategy for nitrogen removal enhancement in full-scale wastewater treatment plants by GPS-X integrated with response surface methodology. Sci Total Environ. 2021;769:144851.10.1016/j.scitotenv.2020.144851Search in Google Scholar PubMed

[8] Alwardy A, Alquzweeni SS, Al-Saadi RJ. Modelling and simulation of Al-muamirah wastewater treatment plant by GPS-X software. Kerbala J Eng Sci. 2021;1(2):113–30.10.1088/1755-1315/877/1/012027Search in Google Scholar

[9] Kobeyev S, Tokbolat S, Nazipov F, Satyanaga A. Design and modeling of an on-site greywater treatment system for a hotel building. Int J Build Pathol Adapt. 2023;41(1):201–24.10.1108/IJBPA-08-2021-0109Search in Google Scholar

[10] Hammed AN, Nile BK, Al-Baidani JH. Improving the performance of the Karbala wastewater treatment plant by using GPS-X Model. International Symposium on Multidisciplinary Studies and Innovative Technologies (ISMSIT); vol. 2022; 2022. IEEE.10.1109/ISMSIT56059.2022.9932706Search in Google Scholar

[11] Zwain HM, Nile BK, Faris AM, Vakili M, Dahlan I. Modelling of hydrogen sulfide fate and emissions in extended aeration sewage treatment plant using TOXCHEM simulations. Sci Rep. 2020;10(1):22209.10.1038/s41598-020-79395-8Search in Google Scholar PubMed PubMed Central

[12] Awad A, Jafar R, Ali K. Modeling Banias refinery treatment plant using the GPS-X. Tishreen Univ J Eng Sci Ser. 2019;41:6.Search in Google Scholar

[13] Mu’azu ND, Alagha O, Anil I. Systematic modeling of municipal wastewater activated sludge process and treatment plant capacity analysis using GPS-X. Sustainability. 2020;12(19):8182.10.3390/su12198182Search in Google Scholar

[14] Sadri Moghaddam S, Pirali MR. Modeling and calibration of a full-scale wastewater treatment plant using GPS-X model (A case study of Tehran). Numer Methods Civ Eng. 2021;5(4):67–76.10.52547/nmce.5.4.67Search in Google Scholar

[15] EPA. Environmental protection agency (EPA). Wastewater Technology Fact Sheet: Oxidation Ditches. Washington, DC, USA: US EPA; 2000.Search in Google Scholar

[16] Metcalf EAO, Mohammad B, Gregory B, Franklin LP, William S, H David Tchobanoglous, et al. Wastewater engineering: treatment and resource recovery. McGraw Hill Education; 2014.Search in Google Scholar

[17] Bankston E, Wang Q, Higgins BT. Algae support populations of heterotrophic, nitrifying, and phosphate-accumulating bacteria in the treatment of poultry litter anaerobic digestate. Chem Eng J. 2020;398:125550.10.1016/j.cej.2020.125550Search in Google Scholar

[18] Khursheed A, Gaur RZ, Sharma MK, Tyagi VK, Khan AA, Kazmi AA. Dependence of enhanced biological nitrogen removal on carbon to nitrogen and rbCOD to sbCOD ratios during sewage treatment in sequencing batch reactor. J Clean Prod. 2018;171:1244–54.10.1016/j.jclepro.2017.10.055Search in Google Scholar

[19] Stewart RD, Bashar R, Amstadt C, Uribe-Santos GA, McMahon KD, Seib M, et al. Pilot-scale comparison of biological nutrient removal (BNR) using intermittent and continuous ammonia-based low dissolved oxygen aeration control systems. Water Sci Technol. 2022;85(2):578–90.10.2166/wst.2021.630Search in Google Scholar PubMed

[20] Figdore BA, Stensel HD, Winkler M-KH. Bioaugmentation of sidestream nitrifying-denitrifying phosphorus-accumulating granules in a low-SRT activated sludge system at low temperature. Water Res. 2018;135:241–50.10.1016/j.watres.2018.02.035Search in Google Scholar PubMed

Received: 2023-09-21
Revised: 2023-11-04
Accepted: 2023-11-08
Published Online: 2024-03-04

© 2024 the author(s), published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Downloaded on 5.6.2024 from https://www.degruyter.com/document/doi/10.1515/eng-2022-0558/html
Scroll to top button