BASES DE LA MODELISATION
HYDROGEOLOGIQUE
2.1. EQUATIONS
Les équations fondamentales permettant la
modélisation de l’écoulement de l’eau souterraine
proviennent de 2 sources :
- la loi de conservation du moment;
-l’équation de diffusivité qui exprime la
conservation de la masse.
2.1. EQUATIONS
✓ la loi de conservation du moment
Loi de Darcy pour un milieu poreux en est un exemple.
h (1)
q=K
x
Avec :
q : Flux volumique ;
K : Perméabilité de Darcy ;
h : Charge hydraulique
2.1. EQUATIONS
✓l’équation de diffusivité qui exprime la
conservation de la masse
Dans un Volume Élementaire Représentatif (VER)
fermé, la variation de la masse de fluide par unité
de temps est égale à la somme algébrique des flux
massiques traversant (entrant et sortant) la surface
du volume considéré.
2.1. EQUATIONS
✓l’équation de diffusivité qui exprime la
conservation de la masse
2.1. EQUATIONS
✓l’équation de diffusivité qui exprime la
conservation de la masse
qx + qy qz h
+ = Ss +W (x, y,t) (2)
x y z t
qx, qy, qz : Flux entrant ou sortant à travers la nappe ;
h : Charge hydraulique ;
Ss : Coefficient d’emmagasinement ;
t : Temps ;
W : Flux volumétrique par unité de volume représentant le terme source.
2.1. EQUATIONS
La combinaison des équations (1) et (2), avec
certaines hypothèses, conduit à l’établissement
d’une équation aux dérivées partielles représentant
le modèle conceptuel de l’écoulement en trois
dimensions dans un milieu poreux :
Equation fondamentale du mouvement:
2.1. EQUATIONS
h h h h
(K xx ) + (K yy ) + ( K zz ) −W = S s
x x y y z z t
Kxx, Kyy, Kzz : Perméabilités selon les trois directions de l’espace (x, y,z);
h : Charge hydraulique ;
Ss : Emmagasinement spécifique du milieu poreux;
t : Temps ;
W : Flux d’échanges verticaux par unité de volume
2.1. EQUATIONS
Cette équation définit entièrement l'écoulement
en permettant la détermination du champ de
charge hydraulique h.
C'est cette équation que les modèles
d'écoulement en milieu poreux s'efforcent de
résoudre.
2.2. REGIMES D’ECOULEMENT
Expression de la variation de la charge hydraulique
h
en fonction du temps ( t ): 2 régimes d’écoulement:
✓ Régime d’écoulement permanent :
tous les paramètres d’écoulement (charges, débits,
vitesses d’écoulement, etc…) sont indépendants du
temps (conditions fixes dans le temps).
h
= 0 (h est considérée constante)
t
2.2. REGIMES D’ECOULEMENT
✓ Régime d’écoulement transitoire :
Au moins l’un des paramètres d’écoulement est
fonction du temps (conditions variables dans le
temps).
h
0 (h varie en fonction du temps)
t
2.3. RESOLUTION DE L’EQUATION
✓Méthode analytique (cas simples): fournit une
expression mathématique explicite pour la charge et le
débit.
Résolution des équations différentielles de 2nd
ordre : n’est possible dans la plupart des cas que
par des méthodes numériques.
2.3. RESOLUTION DE L’EQUATION
✓Méthode numérique (cas complexes) :
résolution par intégration du système d’équations
différentielles.
Solution représente la distribution spatiale de la
charge hydraulique H(x; y; z) en régime permanent,
ou la distribution spatiale et temporelle H(x; y; z; t)
en régime transitoire.
2.3. RESOLUTION DE L’EQUATION
Les équations de base de l’écoulement ne possèdent
pas de solutions analytiques pour des applications
réelles.
✓ la géométrie complexe du domaine;
✓et des propriétés hétérogènes et
anisotropiques des matériaux l’aquifère.
Milieu isotrope : capacité de Milieu anisotrope : capacité de
transfert identique dans toutes les transfert différente dans toutes les
directions (Kx = Ky = Kz ) directions (Kx ≠ Ky ≠ Kz)
2.3. RESOLUTION DE L’EQUATION
Des méthodes et des techniques numériques sont
employées pour approcher la solution.
Ces méthodes permettent de transformer ces équations
en un système d’équations algébriques dont les
inconnues sont les charges hydrauliques.
2.3. RESOLUTION DE L’EQUATION
Ce système est résolu par l’utilisation de techniques
itératives ou d’autres méthodes directes.
Les modèles mathématiques numériques
d’écoulement, à base de différences finies ou
d’éléments finis, permettent d’approcher la
solution de l’équation de l’écoulement.
2.3. RESOLUTION DE L’EQUATION
Méthode des Différences Finies (MDF)
✓consiste à discrétiser le domaine d'étude en blocs
ou mailles (carrée ou rectangulaire), identifiés par
des points ou nœuds ;
✓remplacer l’équation différentielle en chaque nœud
par une série d'équations algébriques de différences
finies issues du développement en séries.
2.3. RESOLUTION DE L’EQUATION
Méthode des Élements Finis (MEF)
✓Consiste à subdiviser le domaine d'étude en sous
domaines ou éléments (triangulaire ou quadrilatère)
délimités par des nœuds (appelés éléments finis);
✓utiliser une approche intégrale (fonctions continues)
pour résoudre les équations différentielles.
2.3. RESOLUTION DE L’EQUATION
2.3. RESOLUTION DE L’EQUATION
Elément Finis
Différences Finies
2.3. RESOLUTION DE L’EQUATION
L’équation fondamentale est valable à l'intérieur du
domaine d’écoulement (domaine à modéliser)
Pour sa résolution, il faut connaître les conditions aux
limites du domaine :
✓conditions de Dirichlet portant sur la variable h
imposée ;
✓conditions de Neumann qui portent sur le flux
imposé h
n
✓ conditions de Fourrier ou Mixte qui portent sur h
et h : h + h imposé.
n n
2.3. RESOLUTION DE L’EQUATION
Parmi les codes de résolutions:
✓ Visual ModFlow
Modèle aux Différences Finies (maillage carré ou
rectanguliare)
✓ Feflow
Modèle aux Éléments Finis (maillage triangulaire)
✓ Etc.
FIN DE LA SEANCE