Water">
Chap 7 - Filtration 2023-1
Chap 7 - Filtration 2023-1
Chap 7 - Filtration 2023-1
1 2
M1 SIE – Traitements physicochimiques
Les filtres
Domaine d’emploi
Filtres ouverts à l’air libre
Eaux résiduaires :
Urbaines : Affinage eau traitée (cas particulier)
Industrielles : après coagulation-floculation – précipitation - décantation Vidéo 1
Vidéo 2
3 4
Nombre : 55 à 90 par m2 suivant la
vitesse de lavage souhaitée.
Qair ≈ 1 m3.h-1 par buselure
Largeur des fentes : 0,5 to 1,0 mm (à
adapter à la nature des matériaux
utilisés)
https://www.oasen.nl/Lists/Downloads/Filtratie%20ontwerp%20en%20inrichting_eng.pdf
5 6
Vidéo
Source : awi
7 8
Filtres fonctionnant sous pression
Filtration Lavage
Vidéo 1
Boues Boues
Filtrat Filtrat
L’eau ne pouvant augmenter de niveau, la mise sous pression est dans ce cas
indispensable pour vaincre les pertes de charge (assurer l’écoulement de l’eau
Vidéo 2 au travers du filtre).
Vidéo 3 9 10
http://www.aquatext.com/images/systemph/sandfilters1.htm
Date 21/09/2020
Référence
Prise d'essai 100 g
N° Tamis AFNOR Mailles (mm) Poids retenu (g) Tamisat (g) Tamisats (%)
37 4 0 98,62 100,00%
graviers sable anthracite 36
35
3,15
2,5
0
0
98,62
98,62
100,00%
100,00%
34 2 0 98,62 100,00%
Matériaux particuliers 33 1,6 4,98 93,64 94,95%
32 1,25 30,71 62,93 63,81%
31 1 45,57 17,36 17,60%
30 0,8 16,77 0,59 0,60%
29 0,63 0,59 0,00 0,00%
28 0,5 0 0,00 0,00%
27 0,4 0 0,00 0,00%
26 0,315 0 0,00 0,00%
Total = 98,62
Perte % = 1,38%
11 12
• La taille effective D10 doit être choisie en fonction de la mise en
œuvre. La tolérance admise par rapport à la valeur nominale annoncée
Courbe granulométrique
par le producteur est de ± 5 %.
100%
90%
80%
Application D10 (mm)
70%
Eaux de piscine 0,55
Tamisats (%)
60%
Filtres bicouche anthracite-sable 0,75
50%
Eau destinée à la consommation humaine 0,95
40%
30%
Traitement tertiaire des ERU 1,35
D ’après P. Corsin, G. Mauguin et N. Villain. Le sable de filtration,
20% L’eau, l’industrie, les Nuisances, N° 278.
10%
0%
0 0,5 1 1,5 2 2,5 3 3,5 4 4,5 • Le coefficient d’uniformité CU doit être inférieur à 1,5. Plus la valeur
D10 D60 Ouverture des mailles (mm)
de CU se rapproche de 1, meilleure est l’homogénéité obtenue pour le
D 1,23 lit filtrant, ce qui permet de réduire les pertes de charge et d’obtenir
Taille effective CU
D 0,91
1,35
une rétention en profondeur des matières en suspension.
taille de la maille du tamis laissant
13 14
passer 10 % de la masse de l’échantillon
• La masse de grains passant par le tamis de maille égale à 66% de celui de Porosité
la taille effective ne doit pas dépasser 1%, et celle ne passant pas par le
tamis de maille égale à 180% de celui de la taille effective ne doit pas Vides
V
dépasser 5%. ε = vides fraction de lit occupé par des vides
Courbe granulométrique
V Lit VLit
100%
≥ 95% Vsable
90%
1− ε = fraction de lit occupé par le matériau
80% VLit
70%
Tamisats (%)
60%
Mesure :
50%
40%
• Ajout d’un volume d’eau V0 donné dans une éprouvette
30%
• Ajout d’un volume de matériau
20%
VSurn
10% 180%.D10 = 1,64 mm • Lecture VSurn et VLit
VLit
≤ 1% 0%
V0 − VSurn
0 0,5 1 1,5 2 2,5 3 3,5 4 4,5
D10 = 0,91 mm
ε=
Ouverture des mailles (mm)
Mécanismes de capture
Mécanismes de filtration
Ils sont essentiellement de deux natures :
Trois mécanismes principaux interviennent successivement :
capture, fixation et détachement. Tamisage mécanique
Leur importance dépend des caractéristiques des particules à Rétention des particules plus grosses que la taille
des « pores » du filtre.
retenir et du matériau filtrant mis en œuvre.
Se produit pour les particules dont le diamètre est
supérieur à 0,15 x diamètre des particules.
18 20
Rétention dans les espaces intergranulaires
Mécanismes de fixation
C’est un mécanisme très important Des réactions avec l’établissement de liaisons covalentes peuvent
dans la filtration en profondeur. également avoir lieu à la surface de l’adsorbant conduisant à une
fixation irréversible du composé.
21 22
turbidité
Les dépôts déjà retenus peuvent alors se détacher partiellement et être tp temps
entraînés plus avant dans le matériau filtrant (progression du « front de
Au cours de cette phase, d’une durée allant de 10 min à deux heures,
filtration ») ou même dans le filtrat (« crevaison »).
l’eau est impropre à la potabilisation et doit être rejetée, ou renvoyée
en tête de station.
Après cette phase, la turbidité reste à une valeur très faible (< 0,1 NTU)
jusqu'à ce que la capacité du filtre à capter les matières soit dépassée,
23
conduisant à la percée du filtre (temps tp). 24
Filtration lente
Vitesse de filtration 0,2 m. h-1 ≤ U ≤ 2 m.h-1
http://slayden.com/falls-city-slow-sand-filter/
la filtration lente
la filtration rapide La filtration lente nécessite toutefois une
occupation au sol importante (procédé
limité à 1000 habitants), et n’est efficace que
pour les eaux peu chargées en MO et MES.
25
26
http://www.elaguapotable.com
Un faible nombre de filtres permet de réduire les coûts des accessoires. Problème 2
Un trop faible nombre de filtres peut toutefois conduire à de fortes Calculer le nombre et la surface des filtres nécessaires à la filtration de
variations de la vitesse de filtration lorsque l’un des filtres est en période 230 000 m3.j-1 avec une vitesse de 15 m.h-1.
lavage. Avec seulement deux filtres, la vitesse est ainsi doublée lorsque
l’un des deux est mis hors service pour lavage , ce qui peut conduire à On rappelle que la surface de filtration d’un filtre unitaire ne doit pas
une détérioration de la qualité de l’eau. dépasser 100 m2.
La plupart des stations ont pour cette raison un minimum de quatre Solution
filtres en fonctionnement.
29 30
Inconvénients
Batterie de filtres alimentés par un même canal largement dimensionné pour que • Possible destruction des flocs du fait de la chute d’eau importante à la remise en
les pertes de charge y soient négligeables, et que par suite le niveau de l’eau y route d’un filtre lavé,
soit sensiblement constant.
• Perturbation de la surface du lit : hauteur d’eau minimale préconisée au dessus du
Débit d’alimentation Q identique pour tous les filtres via des déversoirs lit filtrant : 1,8 m. Plaque déflectrice requise.
32
31
Avantages/inconvénients du
Filtres à vitesse variable – Charge constante
dispositif à charge constante Canal commun
d’alimentation
Niveau 3C Niveau 3F
Niveau 2F Charge disponible
Charge disponible
Niveau 2C
Avantage Niveau 1C Niveau 1F
Apport d’eau
de lavage
Orifice d’alimentation
• Pas de risque de rupture de flocs. La vanne de régulation est ouverte au fur submergé
L’admission de l’eau est réalisée par un même canal sous le niveau de l’eau dans
Inconvénient le filtre : les niveaux d’eau sont les mêmes dans l’ensemble de la batterie de
filtre. Les vitesses de filtration diffèrent en revanche suivant l’état de colmatage
Dispositif couteux. Risque de pannes. des filtres.
33 34
Pertes de charge
Au fur et à mesure du colmatage des filtres, le niveau de l’eau
augmente sur l’ensemble des filtres.
Filtre homogène (CU ≈ 1) alimenté en eau propre
Après lavage, la vanne en sortie de filtre est partiellement fermée pour
éviter que le filtre lavé ait une vitesse de filtration beaucoup plus Loi de Darcy
élevée que les autres. ν
∆H = U .Z
gB
Q
Avantages du dispositif
B : perméabilité de la couche (m2)
• Peu de risque de rupture des flocs par suite d’une chute excessive. ∆Η
• Une vitesse peu élevée sur les filtres les plus colmatés : moins de risque ν : viscosité cinématique du liquide (Pa.s)
d’entraînement de MES dans l’effluent comparé aux dispositifs à vitesse
constante. Z Sable
ν=
.
T : température (°C)
, ,
∆H
= 150ν
(1− ε )
2
U 1− ε U
+ 1,75 3
2
Z ε 3 gd 2
L
ε 3 gd = κV
ε3 gd 2 ε 3 gd
L
Z
Valable pour des particules de sable ou d’anthracite (non sphériques), en prenant pour
d la valeur de la taille effective du matériau et les valeurs suivantes d’ε, de κV et κL : Le modèle d’Ergun permet de calculer la perte de charge observée au début de
la filtration (filtre propre).
Source : Trussell, R. R., and Chang, M. (1999) ‘‘Review of Flow through Porous Media Il peut également permettre le calcul de la vitesse de filtration obtenue en
as Applied to Head Loss in Water Filters,’’ J . Environ. Eng., 125, 11, 998–1006.
début de filtration (filtre propre) lorsque le filtre est conçu pour fonctionner à
charge constante.
39 40
Problème 3
2) En déduire la hauteur d’eau présente en début de filtration au
dessus du filtre représenté à la figure ci-contre.
1) Calculer la perte de charge (en mCE) d’un lit de sable propre d’1 m de
hauteur, de taille effective 0,8 mm alimenté avec une charge superficielle On négligera les pertes de charge liées au plancher et à la conduite.
égale à 20 m3.m-2.h-1.
Canal
On rappelle :
∆H
= κV
(1 − ε ) ν U + κ 1 − ε U 2
2
d’alimentation
ε3 gd 2 ε 3 gd
L
Z Niveau max.
Z =1 m
Solution
41 42
Problème 4 Canal
0,1 m
Cas d’un filtre hétérogène
d’alimentation
1m
Une charge d’1 m est maintenue au dessus Les perméabilités des différentes couches
d’un filtre à sable d’1 m de profondeur. 1m d’épaisseur ∆Z ne sont pas nécessairement
égales.
Le sable a une taille effective de 0,4 mm.
La perte de charge sur l’ensemble du filtre
Calculer la vitesse de filtration (m3.m-2.h-1) en début de filtration pour une doit être calculée comme la somme des
température de l’eau de 20°C. pertes de charge de chacune des couches de
perméabilité donnée :
Viscosité dynamique de l’eau à 20°C: 10-3 Pa/s
On donne :
Masse volumique de l’eau à 20°C: 998 kg.m-3 ∆H = ∆H ∆Z Z : hauteur totale du milieu filtrant (m)
Z
∆H (1 − ε ) ν U + κ 1 − ε U 2
2 ν
On rappelle : = κV avec : ∆H ∆Z = U .∆Z
ε3 gd 2 ε 3 gd gB∆Z
L
Z
∆Η∆Z : perte de charge de la couche ∆Z (m2)
Solution
Β∆Z : perméabilité de la couche ∆Z (m2) 44
43
Problème 5 Cas de la filtration d’une eau chargée
Pertes de charge
Q[MES]0 du plancher
On considère un filtre bicouche constitué d’une couche d’1 m de sable de
taille effective 0,8 mm et d’1 m d’anthracite de taille effective 1,5 mm. Evolution de la perte de charge 1
∆HL,t + ∆Hp
Calculer la perte de charge en début de filtration, pour une charge ∆H L ,t = ∆H L , 0 + k HLσ t (1) 2
superficielle égale à 20 m3.m-2.h-1. Z
V
Constantes plancher
On rappelle :
∆H
= κV
(1 − ε ) ν U + κ 1 − ε U
2 2
Q[MES]ef
∆Z ε3 gd 2 ε 3 gd
L
∆HL,0 : perte de charge du lit propre (σ = 0), m
σt : dépôt spécifique au temps t (masse de matière piégée par unité de volume
du lit au temps t ), kg MES.m-3
alimentation
turbidité
la percée du filtre.
Une fois connue, par expérimentations pilotes ou par expérience, cette capacité
tp temps permet le calcul du temps de percée tp correspondant à une autre épaisseur de lit :
Impact de la taille effective du matériau sur la perte de charge Effet sur la perte de charge
et la capacité du filtre
U ([MES ]0 − [MES ]e )
∆H L ,max = ∆H L , 0 + k HL t HL
Effet sur la capacité Z
La capacité d’un matériau à retenir les particules augmente avec la Loi de type puissance décroissante
diminution de la taille effective, suivant une loi puissance ∆H L , 0
= κV
(1 − ε )2ν U 1− ε U 2
+κL 3
décroissante : Z ε 3 2
gD10 ε gD10 k HL = b2 D10− m2 Constantes
D’autant plus forte que la D10 est faible. La perte de charge augmente d’autant
σ p = b1 D10− m1 b1, m1 : Constantes plus rapidement que la D10 est faible.
En diminuant la taille effective, on peut ainsi augmenter le temps de Une taille effective faible, nécessaire à l’obtention d’une capacité élevée,
percée tp : conduit à une forte perte de charge initiale et un temps tHL plus court.
Les résultats obtenus sont donnés au tableau suivant. L’écart nécessaire entre niveau min et niveau max, pour que le filtre perce
au moment où l’eau atteint le niveau max (tp = tHL). Solution
Q[MES]0
Niveau max.
On admettra : ∆Hp ≈ 0,1 m
D10 Niveau min.
On rappelle :
∆H l ,t = ∆H l , 0 + k HLσ t
hp
53 σ p = b1 D10− m 1
k HL = b2 D10− m2
54
Q[MES]
Optimisation
L’optimisation consiste à trouver, pour une vitesse de filtration U, une Calculs sur
hauteur de lit Z et une charge disponible ∆HL,max choisies, la taille tableur
effective conduisant à un temps de filtration tel que : tp = tHL
U [MES ]0
Principe σp = tp
tp =
b1Z
U [MES ]0 D10m1
Z
σ p ≈ b1.D − m1
10
U ([MES ]0 − [MES ]e )
∆H L ,max = ∆H L , 0 + k HL t HL
Z
Temps
∆H L ,0
= κV
(1 − ε )2ν U 1− ε U 2
+κL 3 tHL
Z ε3 gD102 ε gD10
k HL = b2 D10− m2 tp
t HL =
(∆H L , max − ∆H L , 0 )Z
D10m2
Taille effective D10 (mm)
b2U ([MES ]0 − [MES ]e )
55 56
Fichier Excel
Interprétation
Avec cette taille effective de 0,87 mm l’augmentation de la charge disponible à une valeur
de 3 m n’apporte rien car il faudra de toute façon arrêter la filtration avant la percée,
57 c’est-à-dire un peu avant 47,5 heures. 58
59 60
Effets des autres paramètres Des relations empiriques ont été développées par certains bureaux d’études
pour apprécier les valeurs de tp et tHL en fonction des caractéristiques de
( ∆H L ,max − ∆H L , 0 ) Z ∆H L ,0 (1 − ε )2ν 1− ε U 2 construction et de fonctionnement.
t HL = D10m2 = κV
U
+κL 3
b2U ([MES ]0 − [MES ]e ) Z ε3 gD102 ε gD10
Pour une eau coagulée avec un sel métallique, suez-environnement
b1Z propose par exemple :
tp = https://www.suezwaterhandbook.fr/eau-et-generalites/processus-elementaires-du-
U [MES ]0 D10m1 genie-physico-chimique-en-traitement-de-l-eau/filtration/filtration-sur-lit-granulaire
63 64
Autres valeurs suggérées pour la profondeur et la taille Nombre de filtres
effective à utiliser
Q < 8000 m3.j-1 : minimum deux filtres
Q > 8000 m3.j-1 : minimum quatre filtres
• N = 0,0195.Q0.5
• N : nombre de filtres
• Q : débit de pointe (m3.j-1)
Géométrie B
65 Filtre 66
Pour que le modèle soit représentatif du milieu poreux, il faut : λ : coefficient de pertes de charge
une égalité des volumes des vides (= volume d’eau interne) : La perte de charge dans le lit est alors égale à celle d’un capillaire de longueur Lc
πd p2 et de diamètre apparent dc, où l’eau circule à une vitesse Uc, soit :
εSZ = N .Z p
4
une égalité des surfaces mouillées : Nπd c Lc = A p Lc U c2
Rayon hydraulique
∆H = λ avec dc = 4. Rh,c d’un capillaire
d c 2g
Ap : surface mouillée des particules 67 68
Le problème est alors de relier les différents termes de l’équation à des Nous allons donc utiliser une autre expression de Rh :
paramètres mesurables : U, vitesse en fût vide, ε : porosité du lit, Z : profondeur
du lit. Pour un capillaire de longueur Lc :
Lc
1) Expression de dc
Le volume d’eau interne est : section mouillée .Lc
Pour une conduite non circulaire (cas des capillaires) : dc = 4. Rh
La surface mouillée est : périmètre mouillé.Lc
section mouillée
Rh ,c =
périmètre mouillé Le rayon hydraulique dans le capillaire peut donc également s’exprimer par :
Pour un lit de sable, même modélisé par un ensemble de capillaires, cette
Volume d' eau interne d' un capillaire
définition de Rh n’est pas utilisable car on ne connait ni les sections, ni les Rh =
périmètres mouillés des capillaires. Surface mouillée du capillaire
Pour le lit de sable : Par substitution de Rh dans l’expression de la perte de charge dans le lit
filtrant, on obtient :
Volume d' eau interne Volume d' eau interne Ve
ε= = = Lc U c2 λ 1 − ε A p U c2
Volume lit Volume d' eau interne + Volume des particules Ve + V p ∆H = λ = Lc
4 Rh 2 g 4 ε V 2g
ε p
Ve = Vp
1− ε
2) Expression de Lc
Le rayon hydraulique s’exprime dans ces conditions par :
La longueur des capillaires Lc est plus longue que la hauteur du lit filtrant Ζ.
Pour en tenir compte Carman-Kozeny ont introduit un nouveau facteur
Volume d' eau interne Ve ε Vp correctif, noté T, appelé tortuosité :
Rh = = =
Surface mouillée Ap 1 − ε Ap
Longueur des pores du modèle L
avec : T= =
c
Z Lc
Hauteur du milieu filtrant Z
Ap : surface mouillée des particules
Vp : volume des particules
71 72
3) Expression de Uc On en tire finalement l’expression des pertes de charge sur le lit :
Le débit de l’eau au dessus du lit filtrant étant égal à celui au travers du
milieu poreux, on a : λ 1 − ε A p U c2 λ 3 1 − ε A p
U 2
∆H = Lc = T Z
4 ε V 2g 4 ε 3 V p 2g
p
πd c2
NU c = U .S S : section du filtre (m2)
4 64 16ν A
En régime laminaire : λ= = (1 − ε ) p 1/B
Re UT Vp
Débit au sein Débit au dessus
des pores 2 2
du filtre ∆H
= 2T 2ν
(1 − ε ) 2
Ap U
= kc
(1 − ε )2 Ap ν U
S πd 2 Z ε3 V g ε 3 V p g
p
On en tire : U c = U . ou, puisque ε SZ = N c .Lc
πd c2 4 avec kc = Cte de Carman-Kozeny
N Pour des particules sphériques de diamètre d :
4 L T
(égalité des volumes de vides) : U c = U . c = U . πd 3
εZ ε Vp = N Vp d
=
6
Ap 6 ∆H
= 180ν
(1 − ε ) U 2
Ap = Nπd 2
Z ε 3 gd 2
73 kc = 5 (Carman,1937) Formule de Carman-Kozeny 74