Kinetic studies on the removal of phenol by MBBR from saline wastewater.

BACKGROUND: Phenols are chemical compounds which are included in the high priority of pollutants by environmental protection agency (USEPA). The presence of high concentrations of phenols in wastewaters like oil refineries, petrochemical plants, olive oil, pesticide production and oil field operations contain high soluble solids (TDS) and in an olive oil plant, wastewater is acidic, high salty and phenol concentrations are in the range of 0.1- 1%.
METHODS: Kinetic parameters were calculated according to Monod, Modified Stover- Kincannon, Hamoda and Haldane models. The influence of different initial phenol concentrations on the biodegradation rate was performed. The concentrations of phenol varied from 0 to 500 mg/l.
RESULTS: The value of Ki in saline phenolic wastewater in attached growth systems was higher than suspended growth systems that represented a higher phenol inhibition in suspended growth systems. It was obvious that the best model fitting the obtained data are Hamoda model and the Modified Stover-Kincannon model, having highest R2 values of 0.991 and 1, respectively. The value of Ki in saline phenolic wastewater in attached growth system was higher than suspended growth systems which represented a higher phenol inhibition in suspended growth systems.
CONCLUSIONS: Hamoda model and the Modified Stover-Kincannon model having highest R2 value of 0.991 and 1, respectively, and also predicting reasonable kinetic coefficient values.

Background

Phenols are chemical compounds which have been categorized in the high priority of pollutants by environmental protection agency (USEPA) [1]. Phenols are toxic compound at very low concentrations [2, 3]. So, phenol should be effectively removed from phenolic wastewater before discharge into water bodies [4]. The presence of phenols in wastewater such as oil refineries, petrochemical plants, olive oil, pesticide production and oil field operations contain high soluble solids (TDS) [5, 6] and in a olive oil plant, wastewater is acidic, high salty [7]. The concentrations of phenol in refinery and petrochemical industries, are in the ranges of 6–500 and 2.8-1220 mg/ml, respectively [8]. Various methods such as physical and chemical ones exist for the removal of phenol from saline wastewater, but they often have some problem including high costs and also production of hazardous products [2, 5, 6, 9]. The environmental treatment of saline wastewater by means of physic-chemical approaches like reverse osmosis, ion exchange, electrocoagulation etc. is not cost-effective [10, 11].

Kinetic models which related to culture growth and substrate utilization in presence of some inhibitory materials are showed in Table 1: [12]. In this situation, Haldane model is applied for representing the growth in both lower and higher concentrations of inhibitory materials.

Table 1

Kinetic models for growth on the presence of inhibitory substrates (Rozich, Gaudy et al. [11])

μ = μmax. S ÷ [S+ Ks + (S2/Ki)]	(1)
μ = μmax. S [1 + (S/K)] ÷ [S+ Ks + (S2/Ki)]	(2)
μ = μmax. S ÷ [S+ Ks + (S2/Ki)] [1 + (S/Ki)]	(3)
μ = μmax. S exp. (−S/Ki) ÷ (S + Ks)	(4)
μ = μmax [exp (−S/Ki)] – [exp (−S/Ks)]	(5)

Results of Gaudy study are shown with the above given equation [13]. Hamoda proposed a mathematical model based on Monod model making it possible to determine the kinetic coefficients of aerobic submerged attached growth reactor. Kinetic coefficients of a submerged attached growth systems designing the removal of organic carbon are determined by this model (equation 1) [14].

Primitive equation (equation 2) of Modified Stover- Kincannon model designed by Stover for RBC reactor in 1982, has been reported according to the following reaction [15]:

If we replace the area disc with the effective volume of the reactor, in that case, we will have the equation below (equation 3) [15]:

This equation can also be written as a linear response as follows (equation 4):

Hussain et al. performed some experiments to study the kinetic study of aerobic treatment of phenolic wastewater. Kinetics of phenol degradation has also been studied using Haldane model [16]. In addition, Sahariah and Chakraborty were studied the kinetic analysis of phenol, thiocyanate and ammonia-nitrogen removal in an anaerobic-anoxic-aerobic moving bed bioreactor system [17].

Christian et al. studied phenol removal (754 mg/l) by Sulfolobus Solfataricus and calculated the kinetic parameters [18]. In the study of modeling of phenol removal by using attached growth olive pulp bacteria, the kinetic coefficients of μmax, Ks and Ki were determined as 1.296 d−1, 19.23 mg/l,1571 mg/l, respectively [19]. In another research performed by Rozich et al. which phenol was treated by activated sludge, the range of μmax was equivalent to 1.92- 8.64 d−1 (average:4.56 d−1) and amount of Ks was 75 mg/l [20]. Moreover, Yalcin et al. carried out a study on the removal of phenol from wastewater. They reported that the kinetic coefficients of μmax and Ks were found 4.432 d−1 and 87.4 mg/l, respectively. Marques et al. investigated the attached growth biomass and substrate utilization rate in the moving-bed biofilm reactor. They found that no significant difference was observed between the obtained kinetic parameters with those reported in the literature. The organic load and the substrate utilization rates were around 0.8 kg biomass/kg inert carrier which is considered so high [21]. Juang et al. reported kinetics of degradation of phenol in saline solution through solvent extraction with degradation in a two-phase partitioning bioreactor. They found that phenol would potentially be degraded in saline solution and in the aqueous cell medium, although the initial phenol concentrations showed an increasing trend up to 1850 g/m3 [22]. Naghizadeh et al. studied the application of a Hollow- FiberMembrane in removal of COD, TNa and TP from wastewater and reported that the applied method can potentially be considered as an effective remediation method for different applications in wastewater effluent reuse [23]. Mahvi et al. were envestigated the Photo-oxidation of phenol in water media: Toxicity of intermediates [24]. Maleki et al. studied the degradation and toxicity reduction of phenol by ultrasound waves. In this study pH, kinetic constants and initial phenol concentration on the sonochemical degradation of phenol and toxicity assay were investigated [25]. Ahmadi et al. were used Kocuria turfanesis strain M7, Halomonas Alkaliphila strain R4, Pseudomonas balearica strain Z8 for treatment a high saline petrochemical wastewater, results of kinetic evaluation demonstrated that the yield (Y), endogenous decay coefficient (kd), maximum reaction rate (Kmax), maximum specific growth rate (mmax) and saturation constant (Ks) were 0.54 mg VSS mg COD−1, 0.014 day−1, 1.23 day−1, 0.66 day−1, and 1315 mg L−1, respectively [26]. This study was performed to help bridging this gap in knowledge, by evaluating phenol degradation by MBBR in saline solution and by determining corresponding kinetic rate parameters.

Methods

Reactor setup

Tables 2 and 3 show MBBR and settling tank characteristics used in this study. Synthetic wastewater was fed by a peristaltic pump and the reactor was aerated using a stone air diffuser situated at the bottom of the reactor feed by the compressor. Figure 1 shows a schematic of MBBR used in this study. MBBR reactor was filled with polyethylene media characterized in Table 4; also, to avoid leaving the media from the reactor, a plate screen was installed on top of the reactor.

Table 2

Description of the MBBR reactors used in the pilot- scale

Characteristics (units)	Values
External diameter (mm)	160
Internal diameter (mm)	135
Total height (mm)	650
Internal height (mm)	635
Effective height (mm)	440
Total volume (L)	9.5
The volume occupied by the media (L)	3.2
Effective volume (L)	6.3

Table 3

Description of the settling tank used in the pilot- scale

Characteristics (units)	Values
External diameter (mm)	165
Internal diameter (mm)	135
Total height (mm)	350
Internal height (mm)	325
Effective height (mm)	210
Total volume (L)	9
Effective volume (L)	3

Fig. 1

The schematic of the MBBR experimental set up

Table 4

Description of the media used in the pilot- scale

Characteristics	Values
Specific area (m2. Each)	18
Number in each m2	36,100
Specific surface area (m2.m−3)	650
Density (Kg.m−3)	140
Porosity (%)	87

Reactor operation

The MBBR used in this work was operated with HRT of 6 h for 233 days. Reactor operation was started with 3 L of the municipal activated sludge and 3 L of synthetic wastewater containing acetate and essential nutrients. Municipal sludge characterization was evaluated in terms of mixed liquor suspended Solids (MLSS), sludge volume index (SVI) and microscopic assay; furthermore, microscopic assessment was done during experiments. The samples were also collected from influent and sampling port of each reactor. The applied synthetic wastewater consisted of acetate phenol as carbon source and nutrient containing dipotassium mono hydrogen phosphate, potassium dihydrogen phosphate, sodium bicarbonate, ammonium chloride, magnesium sulfate, calcium chloride, boric acid, ferric chloride, copper sulfate, Iodide potassium, chloride, manganese, zinc, cobalt chloride, and magnesium thiosulfate (Merck company, Germany) resulting in C/N/P as 100/5/1. In order to prevent decreasing pH, the alkalinity was enhanced through addition of sodium bicarbonate, pH value range of aerobic biological process was maintained and dissolved oxygen in the reactor was maintained at a range of 1–2 mg. In adaptation stage, phenol concentration showed an increasing during 120 days and simultaneously acetate decreased stepwise until phenol increased to 100% of the target amount and acetate reached to 0 mg/l. Effect of phenol concentration (0–500 mg/l) on the MBBR performance was also carried out in three phases 0–100, 100–300 and 300–500 mg/l, where phenol concentrations increased gradually during 105 days of reactor operation with TDS = 0. A total of 105 days of reactor operation comprised 15, 23 and 21 days for the three steps and 46 days between three steps respectively were attributed to evaluation of the effect of phenol concentration and effluent stability. Moreover, the effect of TDS concentration (0-3%) on the MBBR performance was evaluated during 114 days with a gradual increase in TDS in three phases 0–1, 1–2 and 2–3% with a constant concentration of 500 mg/l phenol. A total of 114 days of reactor operation comprised 18, 18 and 13 days of long operation for the three steps and 65 days between three steps respectively were attributed to evaluate the effect of TDS concentration and effluent stability. Seven days of constant phenol concentration in the effluent were defined as stability of the reactor.

Microorganisms and culture acclimatization condition

The Moving-bed Biofilm Reactor in this study for the growth studies was performed with HRT of 6 h over a period of 233 days. Reactor was operated with activated sludge obtained from municipal wastewater treatment. The applied synthetic medium in the reactor for growth studies is reported in Table 5.

Table 5

Stock solutions used in growth medium solutions used

Compound	mg/l
K2HPo4	58
KH2Po4	25
CaC12	50
MgSo4	75
Na2Co3	200
NH4Cl	191
H3Bo3	1500
Fecl3 .6H2o	150
CuSO4 .5H2O	30
MnCL2 .4H2O	30
ZnSO4 .7H2O	120
COCl2 .6H2O	60
Na2MnSO4 .2H2O	150

In adaptation stage, phenol concentration increased during 120 days and simultaneously acetate decreased stepwise so that phenol evaluated to 100% of the target amount and acetate reached 0 mg/l. Effect of phenol concentration (0–500 mg/l) over the MBBR performance conducted in three separate phases including 0–100, 100–300 and 300–500 mg/l. Phenol concentration witnessed a gradual increasing trend during 105 days of reactor operation with TDS = 0. A total 105 days of reactor operation consisted of 15, 23 and 21 days operation for the three steps and 46 days between three steps, respectively. They were attributed to the evaluation of effect of phenol concentration and effluent stability. Moreover, the effect of TDS concentration (0-3%) on the MBBR performance was evaluated during the 114 days with which showed a gradual increase in TDS in three phases 0–1, 1–2 and 2–3% with a constant phenol concentration of 500 mg/l. A total of 114 days of reactor operation consisted of 18, 18 and 13 days of long operation for the three steps and 65 days between three steps respectively were ascribed to the evaluation of effect of TDS concentration and effluent stability. Seven days of constant phenol concentration in the effluent were defined as stability of the reactor.

Analytical procedure

During the period of continuous operation, some samples were obtained from the influent and effluent of the MBBR every day. The concentration of phenol was determined using spectrophotometry and by the colorimetric 4-aminoantipyren procedure as detailed in the Standard Methods (0.002 to 0.200 mg/L) [27] using a Spectrophotometer DR 5000 UV/VIS (HACH company, USA). TDS was also measured by means of a TDS meter (HACH company, USA) as well.

Modeling the kinetics of phenol biodegradation

The principal aim of the present research was to define the growth and biodegradation kinetics of microorganisms at phenol concentration of 500 mg/l and salt concentration of 3% with HRT of 6, 12, 18, 24 and 36 h. For determination of the amount of biomass concentration at each stage, some number of media were removed randomly from the MBBR reactor which was followed by weighting. By measuring the weight difference of media with initial weight of media at each stage, biomass concentration was calculated. For conduction a comprehensive research, various kinetic models for degradation of phenol were applied. The kinetic parameters were also calculated according to Monod, Modified Stover- Kincannon, Hamoda and Haldane models. In Monod model, K and Ks were calculated by plotting of S-S0/X vs. 1/S, Y and Kd were calculated with plotting 1/S vs. S0-S/X. The Haldane equation is used (μ = μmax. S ÷ [S + Ks + (S2/Ki)]) for determination of the inhibitory coefficient (Ki). In Hamoda model, K and Ks were calculated by plotting 1/S vs. [A /Q (S0-S)]. The slope of regression line was equal to K/Ks and the intercept of regression line was equal to 1/Ks by plotting of [(S0-S)/X] vs. (A /QX), the slope of the regression line was equal to Kd/Y and intercept was equal to 1/Y. The Modified Stover-kincannon were calculated KB and Umax by plotting V/Q (S0 - Se) vs. V/QS0.

Results

The kinetic coefficients were determined at phenol concentration of 500 mg/l and salinity of 3%. A plot of 1/S vs. X/S0-S (Fig. 2) was also depicted for evaluation of K and Ks for the phenol biodegrading microorganisms, based on Monod equation. Therefore, the amounts of K and Ks were calculated as 16.47 d−1, 130.32 mg/l, respectively.

Fig. 2

Regression between X/S0-S and 1/S to determine Ks and K using Monod model

Assuming a specific decay rate, Kd of 0.0282 d−1 and Y of 0.4391 mg/mg were estimated by statistically fitting the observed data to an attached-growth Monod-type rate expression using linear regression techniques according to Fig. 3. The amount of Ki in Haldane equation, (μ = μmax. S ÷ [S + Ks + (S2/Ki)]) was calculated as 1641 mg/l. The kinetic coefficients of K, Ks, Y and Kd in Hamoda model, is shown in Figs. 4 and 5. The kinetic parameters of K, Ks, Y and Kd were calculated 1.1 d−1, 91.74 mg/l, 0.482 mg/mg and 0.011 d−1, respectively. The kinetic parameters of Umax and KB in modified Stover - Kincannon model, according to Fig. 6, were calculated 47.61 g/l.d and 13.47 g/l.d, respectively. This study was investigated the growth and biodegradation kinetics at 500 mg/l of phenol and salt concentration of 3% with HRT of 6, 12, 18, 24 and 36 h.

Fig. 3

Regression between S0-S/X θ and 1/θ to determine Y and Kd using Monod model

Fig. 4

Regression between [A /Q(S0-S)] and 1/θ to determine K and Ks using Hamoda model

Fig. 5

Regression between (A /QX) and [(S0-S)/X] to determine Y and Kd using Hamoda model

Fig. 6

Determination of Kinetic coefficients of the modified Stover - Kincannon model

Discussion

Kinetic studies are very important for industrial laboratory, in order to justify kinetic studies to generalize why using their results in the industrial scale [24]. In the study reported by Nadafi et al., the effect of lyophilized on the kinetic coefficients of the activated sludge was evaluated. The amounts of K, KS, Kd, Y and μmax were obtained as 3.2 d−1, 71.7 mg/l, 0.077 d−1, 0.48 mg/mg, 1.5 d−1, respectively [25]. Kinetic parameters such as K, Ks, Y and Kd according to Hamoda model were calculated as 1.1 d−1, 9.74 mg/l, 0.479 mg/mg, 0.01 d−1. The kinetic coefficients of Modified Stover – Kincannon model, Umax and KB were calculated, 47.61 g/l.d and 47.13 g/l.d, respectively. The amount of kinetic coefficients of K, Ks, Kd, Y and μmax in Monod model were calculated 47.16 d−1, 130.32 mg/l, 0.0282 d−1, 0.439 mg/mg and 7.23 d−1, respectively. The amount of the Ki according to Haldane equation, (μ = μmax. S ÷ [S + Ks + (S2/Ki)]) was calculated 1641 mg/l. The growth kinetics parameter values obtained from different models are shown in Table 6.

Table 6

Summary of growth kinetics parameter values obtained from different models during biodegradation of phenol

Model	μmax(d−1)	KS(mg/l)	K (d−1)	Kd(d−1)	Ki(mg/l)	Y (mg/mg)	Umaxg/l.d	KBg/l.d
Monod	7.23	130.32	47.16	0.0282	–	0.439	–	–
Haldane	–	–	–	–	1641	–	–	–
Hamoda		9.74	1.1	0.479	–	0.01	–	–
Stover-Kincannon	–		––	–	–	–	47.61	47.13

Tallaei et al. studied the kinetic parameters of crude oil decomposition by means of pseudomonas aeruginosa bacteria, The amount of K, KS, Kd and Y were also determined as 9.39 mg/l, 169.3 mg/l, 0.1071 d−1, 0.882 mg/mg, respectively [28]. For finding the kinetic parameters of phenol in saline wastewater in high concentration of 50 mg/l, μmax, K and Y were calculated as 5.28- 9.36 d−1, (33.4- 27.7) and 0.19- 0.28 mg/mg, respectively. By increasing phenol concentration up to 300 mg/l, the amount of μmax, Y, K were decreased to 5.28- 2.16, 0.2-0.29, 10.2- 18.2, respectively [29]. Results of this study were similar to present study. Rozich et al., reported that phenol was treated by activated sludge, the range of μmax and amount of Ks were 1.92-8.64 d−1 (by average: 4.56 d−1) and 75 mg/l in their investigation, respectively [19]. In the study by Yalcin et al., about phenol removal from wastewater, the amounts of μmax and Ks were calculated as 4.432 d−1 and, 87.4 mg/l, respectively. The findings of that research were different from present one [30]. In another study performed by Christian et al., the removal of phenol (754 mg/l) was found by Sulfolobus solfataricus, the quantities of μ, Ks, Ki and μmax were also calculated as 2.256 d−1, 77.7 mg/l, 319.4 mg/l and 1.128 d−1, respectively [17]. Moreover, in the study of modeling phenol by using attached growth olive pulp bacteria, μmax, Ks, Ki were calculated as 1.296 d−1, 19.23 mg/l, 1571 mg/l, respectively. The inhibitor constant, in this study, was found 1641 mg/l [18]. The best model was selected for the evaluation process by evaluation of the correlation coefficients (R ). According to Figs. 1 and 2, the correlation coefficients of Monod model were obtained 0.986 and 0.9771, respectively. According to Figs. 3 and 4, the correlation coefficients of Hamoda model were obtained 0.9886 and 0.9913, respectively. According to Fig. 5, the correlation coefficients of Modified Stover-Kincannon model were also equal to 1. Results revealed that the best models for fitting the experimental data of this study are Hamoda and Modified Stover- Kincannon models, having highest R value and predicting reasonable kinetic coefficient values. The literature reported that performance of both Hamoda and Modified Stover-Kincannon models were more suitable for designing the aerobic submerged attached growth biological reactors [13, 31, 32].

Conclusion

In this research, we studied the growth and biodegradation kinetics of microorganisms at 500 mg/l phenol concentration and salinity of 3%. According to Monod equation, K, Ks, Kd, Y and μmax were calculated 16.147 d−1, 130.32 mg/l, 0.0282 d−1, 0.4391 mg/mg and 7.23 d−1, respectively. Kinetic parameters like K, Ks, Y and Kd, according to Hamoda model, were calculated as 1.1 d−1, 91.47 mg/l, 0.482 mg/mg, 0.011 d−1, respectively. Kinetic parameters of the modified Stover- kincanon model (Umax and KB) were calculated as 47.61 g/l.d and 13.47 g/l.d, respectively. The value of Ki in saline phenolic wastewater in attached growth system was higher than suspended growth systems which represented a higher phenol inhibition in suspended growth systems. Hamoda model and the Modified Stover-Kincannon model having highest R value of 0.991 and 1, respectively, and also predicting reasonable kinetic coefficient values.