Aérage

J.E. KISSAI

1 10/17/2021
 2ème partie
Calcul des réseaux d’aérage

2 10/17/2021
 Réseau d’aérage
Quelques définitions:
Tête de charge (H) : hauteur de la colonne d'eau équivalent à la
pression exercée par l'air à cet endroit; utilisée pour mesurer les
différences causant un écoulement d'air; Pa
Viscosité absolue (μ): trainée, ou résistance de l'air au déplacement; en
Pa.s:
μ = (17 + 0.045 x ºC) x 10-6
Viscosité cinématique (v): rapport de la viscosité absolue à la masse
volumique de l'air; m2/s

3 10/17/2021
 Réseau d’aérage
Quelques définitions:
La loi des gaz parfaits:
La loi d’Avogadro stipule que pour des volumes égaux de gaz idéaux à la
même température contient le même nombre de moles.

La masse volumique de l’air :

En plus :
mR = constante ou MR = constante
Pour un gaz universel Ru = 8314.36 J /K.

Et donc :
4 10/17/2021
 Réseau d’aérage
Quelques définitions:

Exercice n°2 :
Calculez le volume de 100kg de méthane à la pression de 75kPa et à
la température de 42°C.

5 10/17/2021
 Réseau d’aérage
Quelques définitions:

L’équation de Bernoulli : cas d’un écoulement sans

échange de travail,

6 10/17/2021
 Réseau d’aérage
Quelques définitions:
Equation de Bernoulli:

7 10/17/2021
 Réseau d’aérage
Quelques définitions:
L’équation de Bernoulli : cas d’un écoulement avec
échange de travail,

Pnet est la puissance nette transmise à l’air

qm et le débit massique de l’air

8 10/17/2021
 Réseau d’aérage
Quelques définitions:
La loi de Bernoulli entre deux points (1) et (2) :

9 10/17/2021
 Réseau d’aérage
Ventilation naturelle:

Deux éléments peuvent contribuer à la générer :

• la différence de densité de l’air entre les accès d'entrée et
de sortie
•la différence des températures entre l'entrée et la sortie
Plus les différences d'élévations et de températures sont
prononcées, plus la pression de ventilation naturelle est
importante.
Cet apport d'énergie n'est pas fiable, puisque le courant
s'inverse dans le temps et que son intensité varie d'une journée
à l'autre.
10 10/17/2021
 Réseau d’aérage
Ventilation naturelle:

Pour un circuit simple, il faut considérer 2 colonnes d'air de même hauteur. La

procédure suivante est alors conseillée pour définir la direction
d'écoulement de l'air :
a) les colonnes d'air doivent partir du point le plus profond de la mine et se
terminer au point le plus haut;
b) la température de l'air à l'extérieur de la mine est supérieure à celle de
l'air de la mine en été et vice-versa en hiver;
c) la colonne d'air la plus lourde (la plus froide) a tendance à tomber et
déplacer l'air chaud, plus léger;
d) la direction d'écoulement dans la mine se fera donc de la colonne la plus
pesante vers la colonne la plus légère

11 10/17/2021
 Réseau d’aérage
Ventilation naturelle:

12 10/17/2021
 Réseau d’aérage
Ventilation naturelle:

La pression peut être estimée de la différence du poids

volumique des deux colonnes d’air dans les puits. Elle peut être
déterminée des différences de températures. Il s'agit de la
méthode "la plus précise", puisqu'elle tient compte de l'humidité
de l'air. On a:

Dans ce cas, on a les mêmes profondeurs L.

13 10/17/2021
 Réseau d’aérage
Ventilation naturelle:

Quand on n’a pas la même hauteur, on calcule le poids du puits

en ajoutant une colonne au dessus du puits le plus petit afin
d’avoir la même cote:

14 10/17/2021
 Réseau d’aérage
Ventilation naturelle:

En utilisant les températures à l’entrée et à la sortie, on aura:

NVP pression naturelle (en Pa)

g l’accélération
D la profondeur de la mine (en m)
B la pression barométrique moyenne (en kPa)
T température (en °K)

15 10/17/2021
 Réseau d’aérage
Ventilation naturelle:

Exemple d’application:

Calculez la valeur du NVP dans une mine avec les

donnés suivants :
Profondeur de la mine: 300m
Pression barométrique en surface : 101.2 kPa
Pression barométrique en profondeur : 98.1 kPa
Température moyenne en haut du puits: 304 K
Température moyenne au bas du puits : 320 K
16 10/17/2021
 Réseau d’aérage
Ventilation naturelle:

La charge occasionnée par cette pression est donnée par:

en m

d est la densité moyenne de l’air

g est l’accélération de la pesanteur

17 10/17/2021
 Réseau d’aérage
Quelques définitions:
Les pressions:

– Pression barométrique Pb = Pvapeur + Pair

– Pression totale Pt = Pstatique + Pvitesse
– Pression absolue Pabs = Pt + Pb

18 10/17/2021
 Réseau d’aérage
Quelques définitions:
Le nombre de Reynolds, sans dimensions, est une variable dépendant
principalement des dimensions et forme de l’ouvrage et de la vitesse
d'écoulement du fluide :

Où D : diamètre équivalent de l’ouvrage, en m

V : vitesse d'écoulement du fluide, en m/s
μ : viscosité absolue

19 10/17/2021
 Réseau d’aérage
Quelques définitions:
Le type de courant influence le nombre de Reynolds :
– Laminaire NR < 2000
– Intermédiaire 2000 ≤ NR < 4000
– Turbulent NR ≥ 4000
Dans une conduite sans changement de direction ou de dimensions, la
vitesse au centre est maximale (Vmax); il existe une relation entre la
vitesse moyenne et la vitesse au centre : Vmoy ≅ 0,8 Vmax dans les
mines.

20 10/17/2021
 Réseau d’aérage
Quelques définitions:
Exercice :
A travers un puits d’aérage de 5m de diamètre passe un débit de
200m3/s. la masse volumique de l’air est de 1,2kg/m3 . Calculez la
valeur de Re et déterminez le régime d’écoulement d’air. La
température de l’air est de 20°C.

21 10/17/2021
 Réseau d’aérage
Pertes de charges dans la mine:

HL pertes de charge totales dans une galerie

Hf perte de charge par frottement
Hc pertes de charge par choc sur obstacles ou singularité

22 10/17/2021
 Réseau d’aérage
Pertes de charges dans la mine:
Pertes statiques dans une galerie : écoulement laminaire,

A partir de la loi de Hagen-Poiseuille (1799-1869) on tire :

μ viscosité dynamique de l’air

L longueur de l’ouvrage en m
R de l’ouvrage en m
Q le débit d’air en m3/s
23 10/17/2021
 Réseau d’aérage
Pertes de charges dans la mine:
Pertes statiques dans une galerie : écoulement laminaire,

Sachant que la perte de charge peut s’écrire:

Alors :

RL est la résistance laminaire.

24 10/17/2021
 Réseau d’aérage
Pertes de charges dans la mine:
Pertes statiques dans une galerie ou par frottement:
Ecoulement turbulent, Loi de Chezy Darcy,

f coefficient de frottement
L longueur de l’ouvrage en m
per périmètre en m
A la section en m²
ρ la densité
u la vitesse en m/s
25 10/17/2021
 Réseau d’aérage
Pertes de charges dans la mine:
Pertes statiques dans une galerie ou par frottement:

En posant: k = ρ f/2 et u = Q/A on obtient:

k facteur de frottement d’Atkinson

L longueur de l’ouvrage en m
per périmètre en m
A la section en m²
Q le débit d’air en m3/s
26 10/17/2021
 Réseau d’aérage
Pertes de charges dans la mine:
Pertes statiques dans une galerie ou par frottement:
Où:

R s’appelle la résistance aéraulique de la galerie.

R est aussi la résistance de turbulence.
Qui donne une perte de charge:

Le travail contre les forces de frottement:

27 10/17/2021
 Réseau d’aérage
Pertes de charges dans la mine:
Pertes statiques dans une galerie ou par frottement:
Détermination de la valeur f selon Nikuradse:

28 10/17/2021
 Réseau d’aérage
Pertes de charges dans la mine:
Pertes statiques dans une galerie ou par frottement:
Détermination de la valeur f selon Moody:

29 10/17/2021
 Réseau d’aérage
Pertes de charges dans la mine:
Pertes statiques dans une galerie ou par frottement:
La valeur de f dans le régime laminaire:

La valeur de f dans le régime avec de très faible turbulence,

selon Nikuradse et Theodore Von Kármán (1881-1963):

30 10/17/2021
 Réseau d’aérage
Pertes de charges dans la mine:
Pertes statiques dans une galerie ou par frottement:
La valeur de f selon Paul R.H. Blasius (1873-1970),
3000 < Re < 105 :

20 000 < Re < 107 :

31 10/17/2021
 Réseau d’aérage
Pertes de charges dans la mine:
Pertes statiques dans une galerie ou par frottement:
La valeur de f dans le cas d’une parois rugueuse:
Von karman,

d : 2 fois le rayon hydraulique (d = 4 A/Per)

e : est la plus grande aspérité

32 10/17/2021
 Réseau d’aérage
Pertes de charges dans la mine:
Pertes statiques dans une galerie ou par frottement:
La valeur de f dans le cas d’une parois rugueuse:
Équation de Colebrook White ,

33 10/17/2021
 Réseau d’aérage
Pertes de charges dans la mine:
Pertes statiques dans une galerie ou par frottement:
Sachant que: k = ρ f/2
On a d’après Von Karman:

Exemple:
Pour une galerie de longueur 550m, à section rectangulaire de
4×2,5 m². L’aspérité la plus grande est de 0,125m.
Calculez le facteur de frottement et le coefficient d’Atkinson pour
une densité de 1,2kg/m3.
En supposant qu’on est en régime turbulent, calculez la résistance
34 de la galerie et la perte de charge pour un débit de 250m310/17/2021
/s.
 Réseau d’aérage
Pertes de charges dans la mine:
Pertes statiques dans une galerie ou par frottement:
Valeurs de f en fonction de la qualité d’ouvrage, revêtement et
utilisation:

35 10/17/2021
 Réseau d’aérage
Pertes de charges dans la mine:
Pertes statiques dans une galerie ou par frottement:
Valeurs de f en fonction de la qualité d’ouvrage, revêtement et
utilisation:

36 10/17/2021
 Réseau d’aérage
Pertes de charges dans la mine:
Pertes statiques dans une galerie ou par frottement:
Valeurs de f en fonction de la qualité d’ouvrage, revêtement et
utilisation:

37 10/17/2021
 Réseau d’aérage
Pertes de charges dans la mine:
Pertes statiques dans une galerie ou par frottement:
On sait que la résistance dépend du facteur de forme FF = per/A1/2
Forme de la galerie Facteur de forme
Circulaire 1
Arqué avec des montants 1,08
droits 1,09
Arqué avec des montants
évasés
Carré 1,13
Rectangulaire (larg/haut)
1,5 1,15
2 1,20
38 3 1,30 10/17/2021
4 1,40
 Réseau d’aérage
Pertes de charges dans la mine:
Pertes de charges par choc

Ces pertes sont définies comme suit:

X facteur du choc

39 10/17/2021
 Réseau d’aérage
Pertes de charges dans la mine:
Pertes de charges par choc
Notion de longueur équivalente,

Qui donne après simplification:

40 10/17/2021
 Réseau d’aérage
Pertes de charges dans la mine:
Pertes de charges par choc

Exemples:

41 10/17/2021
 Réseau d’aérage
Pertes de charges dans la mine:
Pertes de charges par choc

Exemples:

42 10/17/2021
 Réseau d’aérage
Pertes de charges dans la mine:
Pertes de charges par choc
Exemples:

43 10/17/2021
 Réseau d’aérage
Pertes de charges dans la mine:
Pertes de charges par choc

Exemples:

44 10/17/2021
 Réseau d’aérage
Pertes de charges dans la mine:
Pertes de charges par choc

Exemples:

45 10/17/2021
 Réseau d’aérage
Pertes de charges dans la mine:
Calcul des réseaux: branches en série,

Q = Q12 = Q23 = Q34

HL,TOT = HL,1-2 + HL,2-3 + HL,3-4 = ΣHL
et comme H = RQ2, R = R1-2 + R2-3 + R3-4

46 10/17/2021
 Réseau d’aérage
Pertes de charges dans la mine:
Calcul des réseaux: branches en parallèle,
QT = QA + QB + QC = ΣQi
HL,1-2 = HL,A = HL,B = HL,C
et comme H = RQ2,
HL,A = RAQA2 = HL,B = RBQB2 = HL,C= RCQC2
La résistance équivalente d’un tel circuit est :

La répartition du débit est donnée par:

10/17/2021
47
 Réseau d’aérage
Avantages des galeries en parallèle
• Réduit la résistance à l'écoulement pour un débit
donné, donc réduit l'énergie requise pour faire circuler Q
dans le réseau ;
• Plusieurs points d'apport d'air frais.
Inconvénient principal : dans un réseau à séparation naturelle,
difficulté de contrôle des débits d'air.

48 10/17/2021
 Réseau d’aérage
Exemple:
Trois galeries avec la même résistance R, débit : Q.

• Si ces galeries sont montées en série quelles seront les pertes

de charge totale?
• Si ces galeries sont montées en parallèles quelles seront les
pertes de charge?

Conclusion.

49 10/17/2021
 Réseau d’aérage
Réseau non maillé:
Exemple: calculez les résistances des branches (régime
turbulent),
Branche K(kg/m ) L(m)
3
P(m) A(m ) R 2

AB 3,7 x 10-3 300 18,85 28,27  

BC 3,7 x 10-3 100 18,85 28,27  
CD 3,7 x 10-3 100 18,85 28,27  
DE 16,7 x 10-3 1000 14 12  
EF 18,6 x 10-3 100 14 11,25  
CF 16,7 x 10-3 800 11 7,5  
FG 18,6 x 10-3 100 14 11,25  
BG 16,7 x 10-3 600 11 7,5  
GH 18,6 x 10-3 300 14 11,25  

Le débit total étant de 75m3/s,

Calculez la résistance équivalente, le débit dans chaque branche et la perte
de charge totale.
Quelle est la branche qui occasion le plus de perte de charge?
50 10/17/2021
 Réseau d’aérage
Régulateurs du débit d’air,
Pour rediriger l’air à la direction qu’on veut, on ajoute des
résistances artificielles dans les galeries. C’est ce qu’on appelle des
régulateurs.
La résistance totale du réseau est augmentée, ce qui accroît la
consommation d'énergie par rapport à la séparation naturelle des
débits.

51 10/17/2021
 Réseau d’aérage
Régulateurs du débit d’air,
Exemples:
Supposons que dans l'exemple précédent on exige un débit fixe de
20 m3/s dans la branche BG (au lieu de 26 m3/s).
Quelle résistance artificielle devrait être ajoutée dans le réseau ?

52 10/17/2021
 Réseau d’aérage
Régulateurs du débit d’air,
Exemples:
Supposons que dans l'exemple précédent on exige un débit fixe de
20 m3/s dans la branche BG (au lieu de 26 m3/s).
Quelle charge de ventilation devrait être ajoutée dans le réseau ?

53 10/17/2021
 Réseau d’aérage
Régulateurs du débit d’air,
Sachant que :

On a aussi :

Selon Mc Elroy, pour un orifice circulaire ou carré C = 0.65 (coefficient de

contraction)

54 10/17/2021
 Réseau d’aérage
Régulateurs du débit d’air,
Exemple:
Pour une densité de l’air de 1,2kg/m3,
Calculez la résistance qu’il faut,
Calculez la surface Aartif qu’il faut.

AN: afin de créer une résistance artificielle

de 0,3 Ns2/m8, quelle est la surface du
régulateur qu’il faut ouvrir?

55 10/17/2021
 Réseau d’aérage
Réseaux maillés (lois de Gustav R. Kirchhoff (1824-87))
Première loi, loi des nœuds
La somme algébrique des débits dans les branches aboutissant à
un nœud est nulle (tout ce qui entre au nœud doit en sortir) :
ΣεiQi = 0 avec εi = ± 1
Deuxième loi, loi des mailles
Le long d'un contour fermé, la somme algébrique des différences
de pression entre les extrémités de chacune des branches du
contour est nulle/
(H1 - H2) - (H3 - H2) - (H1 - H3) = 0
Cette équation peut également être exprimée en fonction du
débit et de la résistance, en écrivant HL = R |Q| Q pour
56 10/17/2021
maintenir les conventions de signes.
 Réseau d’aérage
Réseaux maillés :
problématique
Exemple :
Un simple réseau est servi par deux
puits d’entrée et de sortie d’air.
Le débit pompé à l’intérieur de la mine
est de 100m3/s.
Dans la branche centrale un ventilateur
assure un débit de 40m3/s. Déterminez
la répartition de l’air et la pression
nécessaire du ventilateur de la branche
BC.
57 10/17/2021
 Réseau d’aérage
Réseaux maillés : problématique
Exemple :
Nombre des nœuds : 4 (N)
Nombre des branches : 5 (B)
Nombre de mailles : 5-(4-1) = 2
B-(N-1)

58 10/17/2021
 Réseau d’aérage
Réseaux maillés : problématique
Dans un réseau à b branches et n nœuds on peut écrire:

Dans chaque branche : H = f(Q) = R |Q| Q (b équations)

Dans chaque nœud : ΣεiQi = 0 avec εi = ± 1 (n-1 équations)
Dans chaque maille: Σ Hi= 0 (b-n+1 équations)

Ce qui fait 2b équations avec 2b inconnues.

Le problème est que les équations ne sont pas linéaires.

59 10/17/2021
 Réseau d’aérage
Méthode de Hardy Cross
Une fois le réseau établi et les résistances déterminées dans les
branches, il faut commencer par orienter l'air et déterminer une
procédure de rotation dans les mailles (horaire).
On définit les équations suivantes dans chaque maille :
ΣR |Q|Q et 2 Σ R|Q|
Puis on calcule la correction Δ:
Δ = - ΣR|Q|Q / 2 Σ R|Q|
Les Δ obtenus sont les corrections de débits à apporter dans
chaque branche de la maille et peuvent être positifs ou
négatifs par rapport à l'orientation choisie dans la maille.

60 10/17/2021
 Réseau d’aérage
Exemple de Hardy Cross
Soit 47 m3/s à faire circuler dans le réseau. Déterminez les débits et leur
orientation dans chaque branche, puis la perte de charge par branche.
Hypothèse initiale : Q2 = 30 m3/s, Q3 = 17m3/s, Q5 = 27m3/s Q6 = 20 m3/s et Q4 = 3
m3/s de B vers C

Combien de boucles on a? Calculez pour chaque boucle ∆ et faites les

corrections des débits.

61 10/17/2021
 Réseau d’aérage
1ère itération :
Boucle I :
ΣR | Q | Q = 4x10 – 1 x 302 + 10x10 – 1 x 32 – 6x10 – 1 x 172 = 195,6
2ΣR | Q | = 2 (4x10– 1 x 30 + 10x10 – 1 x 3 + 6x10– 1 x 17) = 50,4
Δi = – 195 / (50,4) = – 3,87 m3/s
Boucle II :
ΣR | Q | Q = 5x10 – 1 x 272 – 7x10 – 1 x 202 – 10x10 – 1 x 32 = 75,5
2ΣR | Q | = 2(5x10 – 1 x 27 + 7x10 – 1 x 20 + 10x10 – 1 x 3) = 61
ΔII = – 75,5 / 61 = – 1,237 m3/ s

1ère corrections :
Q2 =30 – 3,87 = 26,13 m3/s Q3 = 17 + 3,87 = 20,87 m3/s
Q4 = 3 – 3,87+1,237 = 0,37m3/s Q5 = 27 – 1,237 = 25,76 m3/s
Q6 = 20 + 1,237 = 21,24 m3/s

62 10/17/2021
 Réseau d’aérage
En fait, la 2ème loi de Kirchhoff s’écrit aussi:

p est la perte de pression par frottement

pf la charge offerte par le ventilateur
NVP pression de la ventilation naturelle

63 10/17/2021
 Réseau d’aérage
En général, la correction ∆Q est donnée par:

Sf et Snv sont les pentes (p, Q) pour le ventilateur et la

ventilation naturelle.
En pratique Snv est toujours pris égal à zéro.

64 10/17/2021
 Réseau d’aérage
La courbe caractéristique:
L’aérage dans la mine est modélisé par la relation:
Le changement du débit suit la loi de similitude:
Exemple:
Pour une mine avec un ventilateur qui pour une charge
statique de 497,7Pa et une charge totale de 749,5 Pa qui
donne un débit de 188,8m3/s. Complétez le tableau et tracez la
caractéristique. Débit dans la charge charge
mine statique totale
0
94,4
188,8 497 745
283,2
65 377,8 10/17/2021
 Fin de la deuxième partie

66 10/17/2021