Original Article

Int. J. Chem. Sci, Volume: 16( 2) DOI: 10.21767/0972-768X.1000266

Modeling Using the Response Surface Methodology the Activation Process of Moroccan Oil Shale for Removal of Cd (II)

Nawal Elhammoudi1*, Mina Oumam1, Abdelkrim Abourriche2, Said Mansouri1, Az-iddin Chham1 and Hassan Hannache1,3 1Laboratory of Engineering and Materials LIMAT, Faculty of Science Ben M?sik, University Hassan II de Casablanca, Casablanca 20000, Morocco 2Laboratory of Materials, Processes, Environment and Quality, Cady Ayyed University, National School of Applied Sciences, Morocco 3Center of Advanced Materials, (EMINES), University Mohammed VI Polytechnic, Ben Guerir, Morocco

*Correspondence: Nawal Elhammoudi , Department of Chemistry, Laboratory of Engineering and Materials LIMAT, University Hassan II de Casablanca, Casablanca 20000, Morocco, Tel:+212660365877; E-mail: [email protected]

Received: April 09, 2018; Accepted: May 12, 2018; Published: May 15, 2018 Citation: Elhammoudi N, Oumam M, Mansouri S, et al. Modeling Using the Response Surface Methodology the Activation Process of Moroccan Oil Shale for Removal of Cd (II). Int J Chem Sci. 2018;16(2):266

Abstract

This study for the first time designed to modeling and optimize of activation process for removal of cadmium (II) ions by adsorption on Moroccan oil shale using the response surface methodology. The optimal activation parameters selection was carried out on the basis of three factors by central composite design (CCD). The factors that have been identified are: activation temperature, activation time and massic ratio (mH2SO4/mYC). The capacity of adsorption of the Cd (II), was chosen as a response. The measured parameters were found in accordance with the predicted values with a coefficient of determination (R2) of 0.92. The model has been confirmed by experimental tests to validate the results obtained. Based on the equation of the model and the desirability function, the activation temperature and activation time significantly effect of the adsorption capacity of Cd (II) ions. According to optimal activation conditions obtained using prediction profiler, the powder of oil shale mixed with H2SO4 at 213.8°C during 30 min in a ratio of 1.76, the adsorption capacity of Cd (II) ions is 68.38 (mg/g), the specific BET surface area of activated oil shale equal to 181.35 (m2/g).

Keywords

Modeling; Activation process; Moroccan oil shale; Adsorption; Cd (II); Response surface methodology

Introduction

Industrial wastewater containing cadmium is a threat to public health, due to its accumulation in aquatic life. It is a heavy metal of considerable toxicity with a destructive effect on most organ systems, this pollutant originating from the metallurgical alloying, industrial operations. The principal sources of contamination are welding contaminated food and beverages [1].

Experimental

The conventional processes for removing cadmium from aqueous solution include chemical precipitation, ion exchange, membranes and adsorption technologies [2]. Among these different physicochemical processes, adsorption has shown to be the best prospects owing to its economic feasibility [3]. Today the oil shale has received increasing interest in the removal of heavy metals. The oil shale is a sedimentary rock contained the organic and the mineral matter [4]; it is essentially an immature form of petroleum. Morocco has a great resource of the oil shale. The deposits that have been explored are the Timahdit and Tarfaya deposits, the Timahdit deposit is located about 250 kilometers southeast of rabat. The average sulfur content is two percent [5]. The deposits that have been strongly explored are the Timahdit and Tarfaya deposits. The total content from Moroccan oil shale is evaluated 86 billion tons [6]. The objective of the present study is to optimize and modeling of activation conditions of oil shale for removal of cadmium (II) ions using response surface methodology. The response surface methodology is one of statistical method to optimize of the parameters and reducing the number of experimental trials [7]. The initial step in planning the response surface design was to determine the activation parameters that have an influence on the removal of cadmium (II) ions [8]. The second step was the modeling of the response; it consists of selecting the appropriate plan and validating the model. For this reason, we employed a central composite design, to optimize of activation conditions of oil shale for removal of cadmium (II) ions. The final step was the optimization of the conditions for the working of the activation process [9]. The factors studied including, the activation temperature (X 1 ), activation time (X 2 ) and massic ratio (X 3 ). The response measured, it is the adsorption capacity of Cd (II) ions: Qads(mg/g). The oil shale is composed of the organic and the mineral material essentially formed by calcite, dolomite, silicate and clay. The choice of the Moroccan oil shale for this application was eluted by two factors, the estimated reserve and its composition rich in organic matter constituted a carbon source.

Preparation of the activated oil shale

The oil shale (Y) is crushed then ground until we get a fine powder [10]. At a mass of the oil shale was leached with hydrochloric acid, to obtain the concentrated oil shale to be named (YC). The precursor (YC) was mixed with different ratio (mH 2 SO 4 /mYC) of the sulfuric acid. The mixture was heated in electric furnace under oxidant atmosphere gaze at different temperatures, after preprocessing at 140°C in oven under air. The product obtained (YCS) is thoroughly washed with distilled water in a soxhlet apparatus, to remove excess sulfuric acid (H 2 SO 4 ) and to extract any soluble material, and then to dry at 80°C.

Experimental response

The measured response in this study was the adsorption capacity of Cd (II) ions. The removal of cadmium (II) ions on the activated oil shale (YCS) has been realized on the synthetic solution from cadmium chloride CdCl 2 . A solution of Cd (II) ions of concentration 20 mg/L was prepared by dissolving 32.62 mg of cadmium chloride CdCl 2 in a 1000 mL. The pH of solution from cadmium chloride CdCl 2 was adjusted with 0.1 N HCl or 0.1 N NaOH. In each adsorption test, 200 ml of cadmium solution of 20 mg/L was mixed with a 20 mg mass of activated oil shale (YCS), and then the mixture is stirred for 24 hours. The equilibrium concentration in Cd (II) ions was determined by atomic absorption spectroscopy (AAS), the quantity of Cd (II) ions adsorbed per unit mass of the activated oil shale (YCS) was calculated using the following form of eqn. (1):

Q ads ((mg/g))=((C 0 -C e )*V)/mYCS (1)

Where Q ads is the adsorption capacity of Cd (II) ions (mg/g), V: volume of the solution of Cd (II) ions (L), C 0 : initial concentration of Cd (II) ions (mg/L), C e : equilibrium concentration of Cd (II) ions (mg/L) and mYCS: mass of the activated oil shale (g).

Response surface methodology

Response surface methodology is a combination of mathematical and statistical techniques used to develop and optimize processes and to determine the relative significance and interaction of all variables, response surface methodology step followed by determination of the optimal region [11]. One of the most popular response surface methods is the Central Composite Design (CCD). The CCD was used to examine the influence of the activation parameters on the response [12,13]. In the central composite design step, experiments were randomly planned to minimize the effect of uncontrolled variables as shown in TABLE 1.

Factors Levels Low (-1) Central (0) High (+1) (X 1 ) 160�C 210�C 260�C (X 2 ) 15 min 30 min 45 min (X 3 ) 1.5 2 2.5

Table 1: Experimental factors and levels in the central composite design.

A central composite design was used to determine the pertinence of the three factors: the activation temperature (X 1 ), activation time (X 2 ) and massic ratio (X 3 ). The first step of this study is the choice of appropriate plan that suits us ideas to optimize the experimental conditions [14]. We have chosen to work with the composite plan on three factors with three levels. The chosen allows studying the influence of these factors on response. The plan consists of three parts: factorial design helps in optimizing the activation parameters [15], whose factors are two levels -1 and+1, experimental point in the center of the study area 0 and axial points on the axes of each factor. For three variables (n=3) and two levels (low (-) and high (+)), the total number of experiments was 16 determined by the expression: (2n (23TABLE 2=8: full factor points)+2n (2*3=6: axial points)+2 (center points)) (TABLE 2). The response measured was the adsorption capacity of Cd (II) ions, the design consisted of 16 experiments point carried out in random order. The suggested model for the quadratic model of the response variable, according to the independent variables [16], the mathematical relationship between the independent variables can be approximated by the second order polynomial model (eqn. 2):

R=a 0 +a 1 X 1 +a 2 X 2 +a 3 X 3 +a 11 X 1 2+a 22 X 2 2+a 33 X 3 2+a 12 X 1 X 2 +a 13 X 1 X 3 +a 23 X 2 X 3 (2)

Matrix of the coded factors Matrix of the reels factors Comment Experimental

value of Qads tr> 1 -1 -1 -1 160 15 1.5 Full factorial 75.1 2 1 -1 -1 260 15 1.5 Full factorial 63.12 3 -1 1 -1 160 45 1.5 Full factorial 50.61 4 1 1 -1 260 45 1.5 Full factorial 62.74 5 -1 -1 1 160 15 2.5 Full factorial 76.99 6 1 -1 1 260 15 2.5 Full factorial 72.69 7 -1 1 1 160 45 2.5 Full factorial 62.74 8 1 1 1 260 45 2.5 Full factorial 16.52 9 0 0 -1 210 30 1.5 Axial point 62.73 10 1 0 0 260 30 2 Axial point 62.74 11 0 1 0 210 45 2 Axial point 45.23 12 0 0 1 210 30 2.5 Axial point 38.68 13 0 -1 0 210 15 2 Axial point 90.91 14 -1 0 0 160 30 2 Axial point 82.74 15 0 0 0 210 30 2 Center point 98.61 16 0 0 0 210 30 2 Center point 87.41

Table 2: Experimental design in coded and reels variables for central composite design.

Where R is the predicted response: adsorption capacity of Cd (II) ions. Xis are the independent variables, (Activation temperature (X 1 ), activation time (X 2 ) and massic ratio (X 3 )), the parameter a0 is the model constant, a 1 , a 2 , a 3 are the linear coefficients, a 11 , a 22 and a 33 are the quadratic coefficients and a 12 , a 13 and a 23 are the interaction effect [17]. The modeling is performed adjusting polynomial equation obtained by analysis of variance (ANOVA) [18].

Each experiment was repeated three times.

Results and Discussion

Acid attack of oil shale

The powder of oil shale (Y-layer) was added, with stirring a solution of hydrochloric acid HCl (1 mol/L) until the complete release of carbon dioxide, according to the following equations.

Activation process of oil shale

Activation processes can be classified into physical and chemical activation. During physical activation, the oil shale may undergo gasification produced during carbonization [19,20]. On the other hand, chemical activation consists of impregnating with chemicals such as ZnCl 2 , H 3 PO 4 , HNO 3 or H 2 SO 4 [21,22]. The concentrated oil shale (YC) was activated by sulfuric acid H 2 SO 4 at the temperature 300°C for 1 hour. The product obtained, referenced by (YCS) where "S" denotes the activation with sulfuric acid H 2 SO 4 . The activation with sulfuric acid H 2 SO 4 has become a widespread method for the manufacture of activated carbon. The use of sulfuric acid (H 2 SO 4 ) has some environmental advantages such as easy recovery, low energy cost, plays two roles during the preparation of activated oil shale (YCS), acts as an acid catalyst to promote bond cleavage, dehydration, and may function as a template because the volume occupied by sulfuric acid in the interior of the activated is coincident with the micropore volume of the activated oil shale obtained [23].

Optimization using the response surface methodology

The variables studied are: the activation temperature (X 1 ), activation time (X 2 ) and massic ratio (X 3 ) and the measured response was the adsorption capacity of Cd (II) ions (Qads(mg/g)). After choosing the response and the variables, we have sought to optimize and modeling of the activation process for removal of Cd (II) ions by oil shale using the response surface methodology. The central composite design to three levels was conducted to evaluate the influence of experimental factors on the adsorption capacity of Cd (II) ions. We realized an adsorption test of the 16 experiments. The experimental data have been reproduced three times to evaluate the statistical quality of the results (TABLE 2).

Evaluation of the quality of the model

The coefficient of adjusted regression R2 for the adsorption capacity was equal to 0.92. This indicated that more than 92% of the variation observed was explained by the direct effects of the factors. Therefore, the choice of quadratic model to modeling of activation process was best. The result of this test is shown in FIG. 1 and TABLE 3.

Figure 1: Graphical representation of observed values as a function of predicted values.

Variable studied Adsorption capacity of Cd (II) Regression analysis R square (R2) 0.92 Adjusted R square (R A 2) 0.91 ANOVA test F value 50.82 p-value <.0001* Validation model F value 3.101 p-value <.0001*

* Significant at the level 95%.

Table 3: Quality of quadratic model chosen to model the activation conditions for removal of Cd (II) ions.

Test of variance (ANOVA)

The validation quadratic model for the removal of Cd (II) ions on activated oil shale (YCS) was performed using the model analysis of variance (ANOVA) and the method of the absence of adjustment analysis, also called analysis bias. ANOVA was used to assess the variance of the model established in relation to the variance of the residue, using the test "Fisher Snedecor", the result was considered significant if (F exp >> F α , ν mod , ν res ), where α=0.05. According to the result of ANOVA test, analysis of experimental give a factor F exp =50.82. The theoretical values determined according to the table of Fisher Snedecor (for ν model =9, ν residue =38 and confidence level=95%) was F α , ν mod , ν res =F 0.05, 9, 38=2.19. This theoretical factor was much lower than the experimental factor, F exp =50.82 >> F theo =2.19. The results of the lack of adjustment of analysis indicate that the p-value was largely experimental factor less than 0.05, these analyses were determined using the JMP12 software, shown in TABLE 3. The results of ANOVA test and the lack of adjustment showed that the quadratic model established was validated.

Equation of the model

The significance of each coefficient was determined using the F-test and p-value given by JMP 12 software. A p-value less than 0.05 indicate the significance of an effect at 95% confidence level [24]. The significant effects and their coefficients are shown in TABLE 4. The effects coefficients of model showed that the adsorption capacity of Cd (II) ions was significantly affected by the quadratic term of activation temperature (X 1 ) and massic ratio (X 3 ), the antagonistic effect of activation time and massic ratio (X 23 ), was not significantly, p-value larger than 0.05, we can neglect this coefficient from the equation of the adsorption capacity of Cd (II) ions. The fitted quadratic model for removal of Cd (II) ions given by eqn. (3).

Q ads (mg/g)=76.15–15.95*X 1 –4.76*X 2 +10.19*X 3 –5.69 *X 12 +4.51*X 13 –12.07*X 1 2+5.95*X 2 2 – 6.50*X 3 2 (3)

Desirability function

After validation of the quadratic model, the optimal conditions of the activation process for removal of Cd (II) ions by oil shale were determined using the desirability function. The function of desirability established to discover the optimal conditions based globally on the function of the desirability of Derringer [25]. The main advantage of the function of desirability is the ability to obtain answers by a simple processing and rapid response.

First and foremost, the response is converted to a function of a particular desirability in the range of 0-1. The value equal to 0 shows a minimum applicability [26]. The prediction plot showing the effects of various parameters, the factors that appear to be more influential are the activation temperature and activation time, allows to conclude that its factors affect the response Q ads (mg/g). It is clear that an increase of the activation temperature and activation time causes an increase in the adsorption capacity of the activated oil shale (YCS).

The optimal condition determined using JMP 12 Software was represented in FIG. 2. On the other hand, the result indicated the highest adsorption capacity of Cd (II) ions of 68.381 mg/g was obtained to 213.8°C, during 30 min with 1.76 of massic ratio. Depending on the results presented in FIG. 2, we can conclude that the desirability function permitted to reduce activation temperature, activation time and massic ratio (m H 2 SO 4 /m YC ) for a high adsorption capacity.

Figure 2: Prediction plot show the effect of parameters: activation temperature (X1), activation time (X2) and massic ratio (X3) on the adsorption capacity of Cd (II) ions.

Experimental validation

To confirm the quality of the results determined by experimental data, experimental tests have been realized to validate the results obtained. The values of the confirmation experiment trials present the same values compared to the optimal conditions for removal of Cd (II) ions on activated oil shale (YCS) obtained using the response surface methodology, no significant difference was observed between the experimental and predicted values, the results of this test are presented in the TABLE 5.

�Response Activation conditions Experimental results Result of (RSM) Adsorption capacity of Cd (II) ions (mg g-1) (X 1 )=213.8°C; (X 2 )=30 min; (X 3 )=1.76 68.73 ± 0.96 68.38 (X 1 )=160°C; (X 2 )=25 min; (X 3 )=1.5 66.58 ± 1.77 67.71 (X 1 )=200°C; (X 2 )=40 min ; (X 3 )=1.7 68.37 ± 0.24 70.94

Table 5: Confirmation experiment of optimal activation process for removal of Cd (II) ions.

N 2 adsorption-desorption

The surface area of the oil shale (Y), concentrated oil shale (YC) and activated oil shale (YCS) are obtained by the BET method (N 2 adsorption in gaseous phase). The Nitrogen adsorption and desorption isotherms of the precursors are given in FIG. 3 and TABLE 6 respectively.

Figure 3: N 2 adsorption-desorption isotherms of oil shale (Y), concentrated oil shale (YC) and activated oil shale (YCS).

Sample Y YC YCS SBET(m2/g) 7.14 23.92 181.35

Table 6: Specific surface obtained by the BET method.

The BET surface area of the (Y), (YC) and (YCS) materials were calculated and found to be 7.14 m2/g, 23.92 m2/g and 181.35 m2/g respectively. The results of analyzes specific surface SBET are summarized in TABLE 6. This table shows that the total surface area of the activated oil shale (YCS) is higher than that obtained by the two precursors (YC) and (Y) which allows that chemical activation with sulfuric acid causes an increase in the surface.

SEM analysis

The scanning electron microscopy (SEM) was carried out for the activated oil shale (YCS) before and after activation in order to evaluate changes on their microstructures. The SEM analysis was executed in the Centre National de Recherche Scientifique et Technique (CNRST). As can be observed in FIG. 4, the rock of oil shale (Y) untreated and the concentrated oil shale (YC) present a rigid and very compact morphology, without any apparent porous structure, while the activated oil shale (YCS) has a porosity much developed compared to the rock of oil shale (Y) of departure where the grains are not consistent with the virtual absence of porosity. The SEM image of activated oil shale (YCS) clearly indicates surface alterations when compared to the rock of oil shale (Y) and the concentrated oil shale (YC), which may be due to oxidation of the organic matter existing in the (YC). The SEM photograph of the activated oil shale (YCS), demonstrate the catalytic role of the sulfuric acid in the chemical activation of oil shale and the development of microstructure.

Figure 4: (A) SEM images of oil shale (Y), (B) SEM image of concentrated oil shale (YC), (C) SEM image of the activated oil shale (YCS).

Conclusion

The activation conditions for removal of Cd (II) ions by oil shale were modeling and optimized using the response surface methodology under central composite design. The activated oil shale (YCS) was successfully prepared using a simple, low cost and highly efficient activation process. Under optimal activation conditions, the oil shale sample mixed with H 2 SO 4 at 213.8°C during 30 min in a massic ratio of 1.76. The adsorption capacity of Cd (II) ions was 68.38 mg/g.

Acknowledgments

I have taken efforts in this project. However, it would not have been possible without the kind support and help of many individuals of the Laboratory of Engineering and Materials (LIMAT) and organizations of the Agence Universitaire de la Francophonie (AUF) for support of this work. I would like to extend my sincere thanks to all of them.

References