Corps de l’article

1. Introduction

Les effluents de l’industrie agroalimentaire sont une des causes principales de la pollution et de la contamination des eaux de surface et souterraines. Leur traitement, avant rejet, par les méthodes conventionnelles ne permet pas la minéralisation totale des polluants contenus dans l’effluent. Leur rejet dans la nature n’est pas sans conséquence sur l’environnement. De nouvelles technologies sont nécessaires pour éliminer et décomposer complètement de tels polluants organiques et minimiser leurs effets sur la nature. La technologie la plus attrayante et la plus prometteuse semble être la photocatalyse. L’utilisation de la photocatalyse hétérogène pour l’oxydation des polluants organiques et inorganiques dans l’eau et l’air a été intensément étudiée durant ces vingt dernières années (OLLIS et AL-EKABI, 1993; SCHIAVELLO, 1997). Ce procédé, faisant partie de la famille des procédés d’oxydation avancée (POA), utilise des semi-conducteurs, tels que TiO2 et ZnO, et peut mener à la minéralisation totale des composés organiques en CO2, eau et acides minéraux.

Bien que le TiO2 soit le photocatalyseur le plus utilisé pour la dégradation d’un large éventail de composés organiques, le ZnO peut représenter une alternative intéressante. En effet, d’après PIRKANNIEMI et SILLANPAA (2002), le mécanisme de photodégradation est semblable à celui du TiO2. De plus, le ZnO peut s’avérer plus efficace pour la dégradation de certains composés organiques tels que le 2-phénylphénol (KHODJA et al., 2001) et le phénol (MARCI et al., 2001). Cependant, certaines études signalent que le ZnO pourrait perdre de son efficacité en raison d’un phénomène de photocorrosion (YASSITEPE et al., 2008). Cependant, ceci ne semble pas rédhibitoire pour des utilisations de courtes durées. Le ZnO retrouve ses capacités initiales après un court séjour dans une solution aqueuse isolée du rayonnement lumineux.

L’objectif de cette étude consiste à étudier la photodégradation du 4-méthylphénol en présence d’UV/ZnO. L’influence de la concentration en polluant, de la quantité de catalyseur et du débit de circulation de la solution sur les performances du procédé sera appréhendée. En utilisant la méthodologie des surfaces de réponse (RSM) et un plan central composite (CCD) (Montgomery, 2001), une modélisation des effets de ces différents paramètres sur le rendement d’élimination et sur l’abattement de la demande chimique en oxygène (DCO) sera réalisée. Cette méthodologie, appelée aussi analyse multivariable, a été employée dans la conception expérimentale de la réaction de Fenton (Herreraet al., 2000; Oliveroset al., 1997) et dans le traitement photocatalytique des effluents industriels (Perezet al., 2001) (Kumar et Kanmani, 2006; CHO et ZOH, 2007).

Notre démarche pourra être résumée en trois étapes essentielles. Une première, durant laquelle le nombre et les niveaux des paramètres à tester sont choisis. Un modèle polynomiale sera proposé et sa validité discutée. Une deuxième étape, basée sur l’utilisation des graphes des effets des facteurs, permettra d’évaluer les effets des différents paramètres sur les performances du procédé. Enfin, en dernière étape, une optimisation et une identification des interactions seront réalisées grâce aux graphes de contours.

En s’appuyant sur cette méthode, on peut faire ressortir l’effet de la concentration initiale en 4-méthylphénol, de la teneur en catalyseur et du débit de circulation de la solution sur les rendements de dégradation et sur l’abattement en DCO. On peut également, déterminer les interactions des paramètres étudiés et les conditions optimales de fonctionnement du procédé. Grâce à cette approche, le nombre d’essais nécessaires est réduit au minimum.

2. Matériels et méthode

2.1 Matériel

Le photoréacteur utilisé se présente sous la forme d’un cylindre à double parois en verre (Figure 1). Le volume utile du réacteur est de 1 200 mL. La source d’irradiation est une lampe UV de marque Phillips référencée PL-L 24W/10/4P. Son spectre d’émission est compris entre 350 et 390 nm avec un maximum à 365 nm. Elle est placée au centre du photoréacteur. La circulation de la solution est assurée par une pompe reliée à un débitmètre qui permet d’ajuster le débit. Ce dernier varie de 0,5 mL•s‑1 à 4,03mL•s‑1.

Le 4-méthylphénol est fourni par la société FLUKA, il a une pureté de 99 %. Appelé aussi para-crésol ou p-crésol, sa formule chimique est CH3C6H4OH. Le catalyseur utilisé est l’oxyde de zinc (ZnO) en poudre de la société Merck. Les concentrations en 4-méthylphénol testées varient de 30 à 150 µL•L‑1 d’eau bidistillée. La concentration de ZnO dans le photoréacteur varie de 0,5 à 3 g•L‑1.

Figure 1

Schéma du dispositif expérimental.

Experimental design.

Schéma du dispositif expérimental.

-> Voir la liste des figures

Le mélange 4-méthylphénol/ZnO est agité pendant 60 minutes à l’obscurité de façon à atteindre l’équilibre d’adsorption. La lampe UV est ensuite allumée (t = 0). Simultanément, on fait barboter un débit d’air de 2,5 L•min‑1 dans la solution de façon à assurer un apport régulier d’oxygène. Les expériences ont été effectuées à pH = 6,8 et à une température de 25 °C. Le pH et la température sont contrôlés avec un appareil multi-paramètre Inolab.

Le suivi de la concentration en polluant est réalisé par prélèvements réguliers d’échantillons de 0,5 mL qui sont centrifugés, filtrés (Millipore 0,45 µm) et analysés à l’aide d’un spectrophotomètre UV-Vis Shimadzu. La longueur d’onde d’absorption maximale est de 277 nm. L’estimation de la DCO a aussi été évaluée par la méthode standard (American Public Health Association, 1998).

Le rendement d’élimination du polluant est ainsi déduit par l’expression suivante : equation: 038331are002n.png où Yi désigne la concentration initiale et Y120 la concentration atteinte au bout de 120 minutes de traitement. Le même type d’expression est utilisé dans le cas du calcul de l’élimination de la DCO.

2.2 Méthodologie des plans d’expérience

La méthodologie classiquement utilisée pour étudier l’influence de paramètres opératoires sur une variable d’intérêt consiste à modifier la valeur d’un paramètre tout en maintenant les autres fixes. L’exploitation des résultats et l’étude expérimentale peuvent être fortement simplifiées en utilisant la méthodologie des plans d’expérience. Cette technique permet de créer un modèle statistiquement significatif d’un phénomène qui intègre les interactions entre les variables tout en optimisant le nombre d’essais.

La méthodologie des surfaces de réponse (RSM) est une technique statistique empirique utilisée pour l’analyse de régression multiple des données quantitatives obtenues à partir des expériences statistiquement conçues en résolvant les équations multivariables simultanément. La représentation graphique de ces équations s’appelle surfaces de réponse, et permet de décrire l’effet individuel et cumulatif des variables d’essai sur la réponse et de déterminer l’interaction mutuelle entre les variables d’essai et leur effet sur la réponse (KHURI et CORNEL, 1987). L’objectif principal de la RSM est de déterminer les conditions opérationnelles optimales pour un système donné qui satisfasse les conditions spécifiques opératoires (MYERS et MONTGOMERY, 2002).

Le concept de surface de réponse modélise une variable dépendante Y, dite variable de réponse, en fonction d’un certain nombre de variables indépendantes (facteurs), x1 , x2 , …, xk, permettant d’analyser l’influence et l’interaction de ces dernières sur la réponse. On peut ainsi écrire le modèle pour une réponse donnée (Y) sous la forme suivante :

où Y est la réponse observée, β0, βi, βij, βii représentent respectivement le terme constant, les coefficients des termes linéaires, des termes représentant les interactions entre variables et des termes quadratiques, et ε représente l’erreur résiduelle.

La méthode des moindres carrés permet de déterminer l’estimateur de β, une matrice qui minimise la somme des carrés des erreurs résiduelles et ses coefficients sont donnés par l’équation (2).

où x est la matrice d’expérience, xt la transposée de x et (x t x)‑1 est l’inverse de la matrice (xt x);

Y est le vecteur des réponses observées, représentant le pourcentage de photodégradation, utilisant les variables indépendantes x1 (Concentration en para-crésol), x2 (Charge en catalyseur) et x3 (Débit de circulation) en valeurs codées.

Le modèle de prédiction est alors donné par l’équation (3) suivante :

Afin de tester la validité du modèle, l’analyse de variance (ANOVA) est utilisée afin d’examiner la signification et l’adéquation du modèle. Ce dernier permet de tracer les surfaces de réponse, d’estimer l’influence et l’interaction simultanées des facteurs sur le rendement photocatalytique.

Comme les facteurs sont en général exprimés dans des unités différentes, leurs effets ne sont comparables que s’ils sont codés. Les valeurs des variables codées sont définies de la manière suivante :

où xi et Xi représentent respectivement les valeurs codées et les valeurs naturelles des ie variables indépendantes, Xi0 la valeur naturelle de la ie variable indépendante au point central et δXi la valeur du pas de variation correspondant.

Dans le but d’optimiser les paramètres photocatalytiques de réaction, chaque paramètre dans la conception (RSM) a été étudié à cinq niveaux différents (-α, -1, 0, +1, +α). Le choix de cinq niveaux pour chaque variable est exigé par cette conception afin d’explorer la région de la surface de réponse à proximité de l’optimum. Le plan central composite (CCD) (BOX et WILSON, 1951; BOX et HUNTER, 1987) est la conception expérimentale la plus employée pour modéliser une surface de réponse de second ordre.

Un plan CCD se compose de Nf = 2k points cubiques (plan factoriel ou fractionnaire), 2k points axiaux (±α, 0, 0,…, 0), (0, ±α, 0,…, 0), (0, 0, ±α,…, 0),…, (0, 0,…,±α) et un nombre de réplications choisi au centre du domaine (0, 0, 0,…, 0) ou points centraux, avec k le nombre de variables indépendantes (facteurs). Un CCD peut être rendu rotatif en choisissant la valeur appropriée de α soit α = (Nf)1/4.

Une conception expérimentale rotative signifie que la variance de la réponse prédite à un point x est fonction seulement de la distance du point au centre de la conception et ne dépend pas de la direction.

Pour nos essais, trois variables indépendantes sont étudiées. Par conséquent, pour k = 3, il faudrait 8 points cubiques, 6 points axiaux et 6 points centraux. Ce qui fait un total de 20 expériences nécessaires pour un modèle polynomial de second ordre. En utilisant trois facteurs, la valeur de α est égale à 1,68.

Les paramètres choisis, dans le cadre de cette étude sont (en codé et naturel) : la concentration initiale en 4-méthylphénol (x1 , X1), la charge en catalyseur (x2 , X2) et le débit de circulation (x3 , X3). Les gammes et les niveaux de ces paramètres sont représentés dans le tableau 1. L’efficacité de la photodégradation constitue la variable de réponse.

Tableau 1

Gammes et niveaux expérimentaux des variables indépendantes d’essais.

Range and levels of experimental parameters.

Variables indépendantes

Gammes et niveaux

α

‑1

0

+1

Concentration 4-méthylphénol (µL•L‑1)

30

54,28

90

125,71

150

Charge en catalyseur ZnO (g•L‑1)

0,5

1,01

1,75

2,49

3

Débit de recirculation (mL•s‑1)

0,5

1,21

2,26

3,31

4,03

-> Voir la liste des tableaux

3. Résultats et discussion

La dégradation par photocatalyse est dépendante de plusieurs paramètres opératoires. Il peut être utile de connaître leur effet et de pouvoir le modéliser. L’utilisation de la méthodologie des plans d’expérience permettra de déterminer l’influence de plusieurs paramètres opératoires sur les performances de la photodégradation.

3.1 Résultats de la conception expérimentale centrale composite

Les résultats des expériences du plan central composite (CCD) sont représentés dans le tableau 2 qui montre les combinaisons statistiques des variables indépendantes x1 (concentration initiale en 4-méthylphénol), x2 (concentration en catalyseur), x3 (débit de circulation), avec les rendements (%) observés (mesurés) et prédits sur une durée de 120 minutes.

Tableau 2

Résultats du plan central composite (CCD) du rendement de la photodégradation du para-crésol et de la DCO pour trois variables d’essais indépendantes dans les unités codées et naturelles avec les réponses observées et prédites.

CCD design matrix for three test variables in coded and natural units along with the observed and predicted responses

No exp.

Valeurs codées

Valeurs naturelles

Rendement (%)

x1

x2

x3

X1

X2

X3

Y(4‑méthyl phénol)

Ypred(4méthyl phénol)

Y (DCO)

Ypred (DCO)

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

‑1

1

‑1

1

‑1

1

‑1

1

‑α

0

0

0

0

0

0

0

0

0

0

‑1

‑1

1

1

‑1

‑1

1

1

0

0

‑α

0

0

0

0

0

0

0

0

‑1

‑1

‑1

‑1

1

1

1

1

0

0

0

0

‑α

0

0

0

0

0

0

54,28

125,71

54,28

125,71

54,28

125,71

54,28

125,71

30

150

90

90

90

90

90

90

90

90

90

90

1,006

1,006

2,49

2,49

1,006

1,006

2,49

2,49

1,75

1,75

0,5

3

1,75

1,75

1,75

1,75

1,75

1,75

1,75

1,75

1,215

1,215

1,215

1,215

3,315

3,315

3,315

3,315

2,26

2,26

2,26

2,26

0,5

4,03

2,26

2,26

2,26

2,26

2,26

2,26

85,05

69,75

98,50

77,39

69,85

52,25

91,04

63,65

96,30

59,85

61,03

86,50

84,12

61,68

68,52

68,01

68,45

68,65

67,10

68,03

84,20

69,35

99,36

76,70

70,05

50,91

90,96

64,01

95,99

60,85

62,22

85,99

84,53

68,11

68,11

68,11

68,11

68,11

68,11

68,11

69,99

63,10

95,60

65,60

58,33

48,25

81,20

52,94

85,00

49,30

51,07

75,03

77,30

53,00

58,95

58,04

59,45

59,05

57,65

57,97

70,72

61,59

95,23

65,46

57,39

47,54

81,62

51,13

84,58

51,26

52,00

75,64

77,55

54,29

58,47

58,47

58,47

58,47

58,47

58,47

-> Voir la liste des tableaux

Les résultats de la conception expérimentale ont été étudiés et interprétés par le logiciel MINITAB 14 (PA, USA) pour estimer la réponse de la variable dépendante dans l’ensemble des expériences.

La conception centrale composite (CCD) se compose de trois variables indépendantes citées plus haut. Les réponses sont l’abattement (%) du 4-méthylphénol et de la DCO. Elles sont représentées par les modèles de second ordre suivants :

En traçant les courbes de parité donnant les valeurs prédites en fonction des valeurs observées (Figures 2a et 2b), on constate que le modèle polynomial d’ordre deux corrèle bien les résultats expérimentaux.

Figure 2

Comparaison des valeurs prédites et observées pour les rendements (a) en 4‑méthylphénol et (b) en DCO.

Comparison of the observed values and the predicted values for elimination of (a) 4‑methylphenol and (b) COD.

(a)

(b)

-> Voir la liste des figures

L’analyse de la variance (ANOVA) du modèle de régression quadratique montre que les modèles sont hautement significatifs au seuil de 1 %. Ceci est justifié par le fait que la valeur du F calculée est beaucoup plus grande que la valeur de F tabulée soit :

equation: 038331are011n.png

equation: 038331are012n.png

Les probabilités (P-valeurs) indiquent la signification du modèle (Tableau 3). La P-valeur relative au modèle de régression est très faible (P-valeur = 0,000), ce qui montre la signification du modèle. De plus, la contribution de chacun des termes linéaires, quadratiques et interactions est hautement significative au seuil de 1 % (P-valeur < 0,01). Par conséquent, on peut conclure que les deux modèles de prédiction (Eq. 5 et 6) sont hautement significatifs et estiment de manière adéquate les données observées.

Tableau 3

Analyse (ANOVA) pour le rendement de photodégradation du 4‑méthylphénol (% rendement) et pour l’abattement de la DCO (% DCO).

Analysis of variance (ANOVA) for 4‑methylphenol degradation (%) and COD elimination (%COD).

Source

Degrés de liberté

F‑valeur

(% rendement)

P‑valeur

(% rendement)

F‑valeur

(% DCO)

P‑valeur

(% DCO)

Régression (modèle)

9

396,28

0,000

200,54

0,000

Linéaire

3

1067,32

0,000

510,23

0,000

Quadratique

3

100,02

0,000

50,57

0,000

Interactions

3

21,50

0,000

40,82

0,000

-> Voir la liste des tableaux

Les valeurs des coefficients de détermination respectivement du 4-méthylphénol et de la DCO sont R2 = 0,997 et R2 = 0,994. Ces valeurs élevées indiquent que 99,7 % et 99,4 % d’une simple variation du rendement photocatalytique du 4-méthylphénol et de la DCO sont expliquées par les variables indépendantes. Les valeurs du coefficient de détermination ajustées sont également très élevées, soit R2adj (4-méthylphénol) =  99,5 % et R2adj (DCO) =  99 %, suggérant une haute signification des modèles (Eq. 5 et 6) et un bon ajustement de ces derniers aux données expérimentales.

3.2 Graphes des effets des facteurs

Les effets individuels, quadratiques et interactions des différents facteurs ont été estimés (Tableau 4). Un coefficient avec un signe (+) signifie que le facteur a un effet synergétique. Par contre, un signe (-) met en évidence un effet antagoniste du facteur.

Tableau 4

Estimation des coefficients de régression pour le rendement d’élimination du 4-méthylphénol (Y1) et pour l’abattement en DCO (Y2) en utilisant les valeurs codées.

Factor effects and associated p-values for two responses.

Relations

 

Facteur

Réponses

Y1

Y2

Effet facteur

p‑valeur

Effet facteur

p‑valeur

Effets principaux

Linéaire

X1

‑41,382

0,000

‑27,729

0,000

Linéaire

X2

27,988

0,000

19,672

0,000

Linéaire

X3

‑26,574

0,000

‑19,672

0,000

Interactions

Quadratique

X12

14,833

0,000

9,601

0,000

X22

8,633

0,000

5,433

0,000

X32

7,389

0,000

7,568

0,000

 

Croisées

X1 x x2

‑5,911

0,000

‑22,059

0,000

X1 x x3

‑3,251

0,009

0,388

0,706

X2 x x3

4,357

0,001

‑0,147

0,866

-> Voir la liste des tableaux

L’effet individuel de la concentration en catalyseur est positif sur le rendement photocatalytique du fait que le coefficient correspondant est positif. À noter, nous n’avons pas observé d’effet de photocorrosion du ZnO lors de nos différents essais.

Par contre, la concentration initiale en 4-méthylphénol a un effet négatif important sur la photodégradation car le coefficient correspondant est affecté du signe négatif. La gamme de concentrations étudiée varie de 30 à 150 µL•L‑1 (Tableau 1). Ceci paraît évident car plus on augmente la concentration en polluant, plus le rendement de photodégradation décroît en raison de la quantité de produit à traiter. Le débit de recirculation influe aussi négativement sur les performances du procédé. Ceci est dû probablement à la diminution du temps de séjour dans le réacteur.

Sur les graphes des effets (Figure 3), on visualise l’évolution du rendement photocatalytique en fonction de chacun des paramètres opératoires étudiés.

Figure 3

Graphes des effets des facteurs sur le rendement de photodégradation à t = 120 min (x1 = concentration en 4-méthylphénol, x2 = concentration en ZnO, x3 = débit de recirculation).

Two-way interaction plots for 4-methylphenol removal at t = 120 min (x1 = concentration of 4-méthylphénol, x2 = concentration of ZnO, x3 = recirculating flow).

Graphes des effets des facteurs sur le rendement de photodégradation à t = 120 min (x1 = concentration en 4-méthylphénol, x2 = concentration en ZnO, x3 = débit de recirculation).

-> Voir la liste des figures

On note bien l’effet antagoniste de la concentration initiale en polluant (x1) et du débit de recirculation (x3). La concentration en catalyseur (x2) a un effet positif dans l’intervalle étudié. On signalera que le même type de comportement est observé concernant l’abattement en DCO.

3.3 Optimisation des paramètres

L’étude des graphes de contours fournit une méthode simple d’optimisation du taux de traitement et d’identification des interactions entre les variables (voir Figures 4 et 5). Chaque courbe représente, dans notre cas, une infinité de combinaisons entre deux variables quand la troisième variable est maintenue à un niveau constant. Ces niveaux constants sont les niveaux centraux de chacune des variables.

Figure 4

Graphes de contours pour le rendement de photodégradation de 4 méthylphénol (%) pour chaque paire de variables en valeurs codées lorsque la troisième variable est fixée à la moyenne.

Contour plots of 4-methylphenol removal (%) for each pair of variables in coded values when the third variable is held constant at the mean value.

(a)

(b)

(c)

-> Voir la liste des figures

Figure 5

Graphes de contours du rendement de la demande chimique en oxygène DCO (%) pour chaque paire de variables en valeurs codées lorsque la troisième variable est fixée à la moyenne.

Contour plots of COD removal (%) for each pair of variables in coded values when the third variable is held constant at the mean value.

(a)

(b)

(c)

-> Voir la liste des figures

Les principales informations que l’on peut en tirer concernant les conditions optimales de fonctionnement du procédé sont :

  • L’efficacité du rendement de photodégradation est inversement affectée par la concentration initiale en 4-méthylphénol (Figure 4a). Ce qui est assez couramment observé en photocatalyse, l’augmentation de la cinétique de dégradation ne compense pas l’augmentation de la quantité de produit à dégrader.

  • L’augmentation de la concentration en ZnO entraîne une augmentation du rendement (Figure 4b). L’effet de la quantité de catalyseur est toujours positif dans la gamme de concentration choisie.

  • La diminution du débit de circulation entraîne l’accroissement du rendement de dégradation (Figure 4c). Ceci s’explique par l’augmentation du temps de séjour de la solution dans le réacteur.

  • L’effet le plus important est celui de la concentration initiale en 4-méthylphénol, suivie par la concentration en catalyseur et enfin par le débit de recirculation.

Les effets linéaires, quadratiques et les interactions sont significatifs pour les deux réponses, hormis les interactions (xx3) et (x2 x3) pour la DCO, où on remarque un faible effet (Figures 5b et 5c). Avec une concentration élevée en catalyseur et un faible débit de circulation de la solution, fixée à la moyenne dans l’intervalle du domaine expérimental choisi, une dégradation complète est atteinte pour une faible concentration initiale en 4-méthylphénol. Les valeurs optimales conduisant à un rendement de 100 % pour les deux réponses sont vérifiées dans notre cas et elles apparaissent sur les figures 4a et 5a. Les valeurs naturelles correspondantes sont respectivement : [4-méthylphénol] = (48,16 µL•L‑1; 43,76 µL•L‑1); [ZnO] = (2,34 mg•L‑1; 2,42 mg•L‑1) et débit de circulation = (1,24 mL•s‑1; 1,22 mL•s‑1). Une validation expérimentale de ces conditions optimales a été réalisée.

Conclusion

La dégradation photocatalytique du 4-méthylphénol en solution aqueuse est étudiée dans un réacteur en circulation. L’influence des trois paramètres, la concentration initiale en polluant, la charge en catalyseur et le débit de circulation, furent étudiés. Les résultats montrent que l’application de la méthodologie des surfaces de réponses (RSM) et le plan central composite (CCD) permettent de décrire et de modéliser d’une manière correcte l’influence de ces trois paramètres expérimentaux sur l’efficacité du traitement. Les valeurs optimales des paramètres donnant un rendement maximal (100 %) ont ainsi pu être déterminées. La fiabilité des deux modèles de prévision de second ordre fondée sur la régression multiple a été testée par la méthode d’analyse (ANOVA). Cette analyse a montré que les modèles sont hautement significatifs et en bonne adéquation avec les résultats expérimentaux.