INTRODUCTION A LA MODELISATION NUMERIQUE

DE L’ECOULEMENT DE L’EAU SOUTERRAINE

1. OBJECTIFS DE L’EXERCICE

Cet exercice est un didacticiel dont le but est de présenter l’interface Visual MODFLOW à
partir d’un modèle représentant les conditions hydrogéologiques d’un site contaminé réel.
Vous devrez suivre les instructions détaillées qui sont décrites dans cet exercice afin de
construire un modèle hydrogéologique simple à partir duquel vous visualiserez les voies de
migration des contaminants.

2. DESCRIPTION DU PROBLÈME

Le site est une chaufferie où un important déversement de diesel a eu lieu au printemps

2000. Ceci a causé une contamination importante des sols et de l’eau souterraine par des
hydrocarbures (diesel). Du diesel en phase flottante a été détecté sur une surface de 60 m
x 60 m.

Les dimensions du site sont de 200 m par 200 m. La stratigraphie consiste en un dépôt de
sable fin de 5 m d’épaisseur, suivi d’un sable moyen de 10 m d’épaisseur reposant sur le
roc. Ce dernier est considéré imperméable. Les conductivités hydrauliques, évalué à
partir d’essais de perméabilité (slug test), sont de 3x10-5 m/s pour le sable fin et de 1x10-4
m/s pour le sable moyen. L’écoulement de l’eau souterraine est du Nord-ouest vers le
Sud-est avec un gradient hydraulique de 0.01. Le niveau de la nappe se trouve entre 3 et
5.5 m de profondeur.
 2

3. GÉNÉRER UN MODÈLE D’ÉCOULEMENT DE L’EAU SOUTERRAINE

Étape 1 : Construire la grille du modèle

Pour commencer un nouveau modèle, cliquez sur :

• File (à partir de la barre de menu supérieure)

Et ensuite sur :

• New (à partir du menu déroulant)

La boîte de dialogue Create New Model sera affichée :

 3

C’est à partir de cette fenêtre que vous allez définir les dimensions de votre domaine, le
nombre de rangées, de colonnes et de couches ainsi que les unités des différentes
propriétés que vous allez entrer plus tard durant l’exercice.

Sous Model Domain, entrer 40 pour le nombre de colonnes (j) et de lignes (i). Entrez
ensuite 2 pour le nombre de couche (k), soit une pour le sable fin et une pour le sable
moyen. Entrez 0 pour Zmin et 15 pour Zmax.

À partir de cette fenêtre, vous pouvez également importer un fichier qui servira comme
arrière-plan du modèle. Pour ceci, cliquez sur :

• Import a site map

Ensuite, vous devez spécifier l'emplacement et le nom du fichier de l'arrière-plan (.bmp).

• Browse

Allez jusqu'à votre répertoire étudiant [K:\étudiant] et sélectionnez le fichier suivant:

• Plan_base.bmp

Cliquez sur :

• Create

La fenêtre suivante apparaîtra :

(107.9 , 97.9)

(44.7 , 47.9)

Cette fenêtre sert essentiellement à géoréférencer le plan de base à partir de points où les
coordonnées X-Y sont connues. Ces points sont les coins inférieur gauche (44.7 , 47.9)
et supérieur droit (107.9 , 97.9) de l’encadré gris (voir figure ci- haut). Pour géoréférencer
 4

le plan, cliquez sur le coin inférieur gauche de l’encadré gris. La fenêtre suivante
apparaîtra et vous pourrez définir les coordonnées de ce point :

Effectuez ensuite la même procédure pour le coin supérieur droit de l’encadré gris. Une
fois que vous aurez géoréférencé le plan de base, la fenêtre suivante apparaîtra :

Resize region

Fenêtre de
dimensionneme
nt

Cliquez sur :

• Resize region

Et déplacez la fenêtre de dimensionnement jusqu’au coin inférieur gauche de l’encadré

blanc. Vous aurez alors défini les dimensions du modèle. Cochez ensuite :

• Show Grid
 5

La fenêtre suivante apparaîtra :

Cliquez finalement sur :

• OK

Visual MODFLOW construira alors une grille de différences finies comprenant 40 lignes,
40 colonnes et 2 couches avec un espacement uniforme des cellules selon les axes X et Y.

Étape 2 : Positionnement du contact stratigraphique entre le sable fin et le sable

moyen

Lors de la génération du maillage, MODFLOW a divisé uniformément les deux couches

du modèle par défaut. Cependant, la hauteur de ces deux couches doit être modifiée afin
de bien représenter l’épaisseur des deux unités stratigraphiques. Pour ceci, cliquez sur :

• View row (voir la portion supérieure gauche de l’interface)

Puis, sélectionnez n’importe quelle ligne du modèle. La fenêtre suivante devrait

apparaître :
 6

Cliquez sur :

• Edit Grid

Puis sur

• Edit Layers

La fenêtre suivante apparaîtra :

Cochez ensuite sur Move, sélectionnez ensuite la ligne qui sépare les deux couches du
modèle et déplacez- là jusqu’à 10 m d’élévation. Vous pouvez visualiser les coordonnées
X-Y- Z correspondant à la position du pointeur de la souris dans le coin inférieur droit de
l’interface.
 7

Coordonnées X,
Y et Z

Étape 3 : Définition des propriétés hydrogéologiques des unités stratigraphiques

Pour définir les propriétés hydrogéologiques des matériaux, cliquez sur :

• Properties

La fenêtre suivante apparaîtra alors :

Cette fenêtre sert à assigner les paramètres par défaut pour toutes les cellules du modèle.
Entrez les mêmes valeurs que celles de la fenêtre ci- haut. Ces paramètres correspondent
à celles du sable fin. Une fois que vous aurez entré les valeurs, cliquez sur :

• OK

Maintenant, vous devez changer les propriétés de la couche inférieure (sable moyen), car
elles sont différentes de celles par défaut. Cliquez sur :

• View Layer

Puis, afin d’aller sur la couche 2, cliquez sur :

 8

• Next (immédiatement à droite de «go to »)

Vous pouvez vérifier si vous êtes bel et bien sur la couche 2 en regardant vers le coin
inférieur gauche de l’interface.

Pour assigner des propriétés différentes pour cette couche, cliquez sur :

• Assign

Et ensuite sur :

• Window

Définissez ensuite un rectangle qui contient l’ensemble du maillage et ensuite appuyez

sur le bouton de gauche de votre souris. La fenêtre suivante apparaîtra :

Pour entrez les valeurs de conductivité hydraulique du sable moyen, cliquez sur :

• New

Entrez une valeur de 1x10-4 m/s pour Kx , Ky et Kx et cliquez finalement sur :

• OK
 9

Étape 4 : Définition des conditions aux limites

Cette section décrit les différentes étapes à suivre pour définir les conditions aux limites
du modèle. Trois types de conditions aux limites seront assignées, soit des charges
imposées, des flux nuls et une limite de recharge. Les limites de charges imposées seront
assignées aux frontières Nord et Sud du modèle de manière à bien représenter le gradient
hydraulique. Les limites à flux nul seront assignées sur les frontières Est et Ouest du
modèle qui correspondent chacune à la position d’une ligne d’écoulement de l’eau
souterraine. La base du modèle sera également une limite à flux nul, car il est supposé
que le socle rocheux est imperméable.

Limite de charge imposée = 12 m

Limite de flux nul

(ligne d’écoulement)

Limite de flux nul

(ligne d’écoulement)

Limite de charge imposée = 9.42 m

 10

Assignation des limites de charges imposées

Cette section décrit les étapes requises pour définir les conditions aux limites de charges
imposées qui doivent être assignées aux frontières Nord et Sud du modèle. Pour ceci,
allez dans:
• Boundaries (barre de menu supérieur)

Sélectionnez :

• Constant Head

Puis, allez dans :

• Assign (À partir de la barre de menu de côté)

Et finalement :

• Line

Définissez une ligne horizontale le long de la dernière ligne de cellule au bas du modèle.
Cette ligne correspond à une équipotentielle de 9.42 m. Déplacez le pointeur de la souris
au coin inférieur gauche de la grille. Cliquez sur le centre de la cellule à l'aide du bouton
gauche de la souris, déplacez le pointeur au coin inférieur droit de la grille et cliquez sur
le bouton droit de la souris au centre de la cellule. Une série de cellules deviendra rose et
la boîte de dialogue Assign Constant Head apparaîtra :

Entrez les valeurs suivantes:

• Code #: 0
• Assign to appropriate layer
• Stop time: 3650
• Constant Head Start Pt.: 9.42
• Constant Head End Pt.: 9.42

La fonction Assign to appropriate layer fait en sorte que les charges imposées ont été
 11

assignées seulement sur la couche #2, car la valeur de charge assignée est inférieure au
contact entre le sable fin et le sable moyen. Les cellules ayant une condition de charges
imposées sont de couleur bourgogne.

Maintenant, assignez une frontière de charge constante sur la bordure supérieure du

modèle :

• Constant Head

Ensuite, allez dans :

• Assign (à partir de la barre de menu de côté)

Sélectionnez

• Line

Définissez une ligne horizontale le long de la dernière ligne de cellule au bas du modèle.
Cette ligne correspond à une équipotentielle de 12.00 m. Déplacez par la suite le
pointeur de la souris au coin supérieur gauche de la grille. Cliquez sur le centre de la
cellule à l'aide du bouton gauche de la souris, déplacez le pointeur au coin supérieur droit
de la grille et cliquez sur le bouton droit de la souris au centre de la cellule. Entrez les
valeurs suivantes:

• Code #: 0
• Assign to appropriate layer
• Stop time: 3650
• Constant Head Start Pt.: 12.00
• Constant Head End Pt.: 12.00

La série de cellules roses deviendra alors bourgogne, indiquant qu'une valeur de charge
constante a été définie. Cependant, les charges constantes sont assignées seulement dans
la couche #1. Pour copier ces cellules à charges constantes dans la couche # 2, allez
dans :

• Copy (À partir de la barre de menu de côté)

Et ensuite sélectionnez :

• Layer
 12

La fenêtre suivante apparaîtra alors :

Cochez ensuite sur

• Constant Head et Layer 2

Et cliquez sur

• OK

Maintenant, vous allez examiner le modèle en coupe transversale pour vérifier si Visual
MODFLOW a bel et bien copié les charges constantes dans la couche inférieure. Pour
ceci, allez dans :

• View Column (à partir de la barre de menu de côté)

Déplacez le pointeur de la souris (barre rouge) dans le domaine du modèle et sélectionnez

n'importe quelle colonne en cliquant sur le bouton gauche de la souris. Vous pouvez
observer que la charge constante a été attribuée aux cellules de la deuxième couche
(cellules rouges). Pour voir plus en détails le modèle, vous pouvez appliquer un facteur
d'exagération verticale :

• F8 - Vert. Exag. (barre de menu inférieure)

• Vertical Exaggeration: 3.5

Assignation de la recharge

La recharge annuelle de la nappe phréatique de ce site a été estimée à 300 mm/an. Avant
d’assigner la recharge au modèle, assurez-vous d’être dans la couche 1 (vérifiez au-
 13

dessous du cube de navigation dans le coin inférieur gauche). Suivez les étapes suivantes
pour spécifier la recharge annuelle:

• Boundaries (à partir de la barre de menu supérieure)

Puis,

• Recharge

La boîte de dialogue Default Recharge apparaîtra. Entrez l’information suivante dans

les cases appropriées:

• Stop Time: 3650

• Recharge: 300

Cliquez sur :

• OK

Note : Par convention, MODFLOW assigne une limite de flux nul aux cellules qui sont
en bordures du modèle si aucune condition frontière (charge imposée, recharge etc.) n’est
assignée par l’utilisateur. Dans le cas de ce modèle, les frontières Est et Ouest ainsi que
la base du domaine agiront comme limite à flux nul.

4. SIMULATION DE L’ÉCOULEMENT DE L’EAU SOUTERRAINE AVEC

MODFLOW

Vous êtes maintenant prêt à simuler l’écoulement de l’eau souterraine avec MODFLOW.
Retournez dans le menu principal de Visual MODFLOW en cliquant sur :

• F10 – Main Menu (barre d’outils inférieure)

Allez ensuite dans le menu :

• Run (menu supérieur)

La fenêtre suivante apparaîtra alors :

 14

Gardez l’option Steady State Flow, car vous simulerez l’écoulement en régime
permanent. Cliquez sur :

• OK

Et ensuite recliquez sur :

• Run (menu supérieur)

La fenêtre suivante apparaîtra :

Cochez seulement :

• MODFLOW 2000

À cette étape, il vous fa ut seulement simuler l’écoulement de l’eau souterraine.

Puis cliquez sur :

• Translate & Run

L’interface Visual MODFLOW va par la suite créer les fichiers nécessaires à l’exécution
du modèle MODFLOW. Après cela, la boîte de dialogue suivante sera affichée :
 15

Des crochets bleus apparaîtront à gauche de Intro et MODFLOW 2000 lorsque la

simulation sera terminée. Après l’apparition de ces crochets bleus, fermez la fenêtre et
procédez par la suite à la visualisation des résultats de la simulation.
 16

5. VISUALISATION DES RÉSULTATS

Allez dans le menu :

• Output (menu supérieur)

La piézométrie résultant de votre simulation apparaîtra alors et devrait ressembler à ceci:

Pour changer les options de contour des équipotentielles, allez dans le menu :

• Option (menu de gauche)

Vous pouvez alors changer, par exemple, l’intervalle des contours et le nombre de
décimales pour les étiquettes de contours. Espacez les iso-contours de 0.25 m et
définissez 2 décimales pour les étiquettes. Cliquez ensuite sur :

• OK

pour sortir de ce menu.

 17

La zone de couleur olive qui est présente au bas du modèle est une zone de cellule sèche.
Ces cellules sont en fait inactivées par MODFLOW, car le niveau de la nappe phréatique
est inférieur à la base de la cellule.

Cellule active Cellule inactive

ou sèche

Cellule active Cellule active

Visualisation des vecteurs de vitesses

Pour visualiser les vecteurs de vitesses, allez dans :

• Map (menu supérieur)

Et ensuite :

• Velocities

Si vous visualisez les vecteurs de vitesses sur la couche #1 et ensuite ceux de la couche
#2, vous observerez que la longueur des flèches est différente. En effet, Visual
MODFLOW dimensionne les vecteurs de vitesses selon leur magnitude.

Vous pouvez également observer les vitesses selon leur direction, ce qui est utile pour
visualiser des vecteurs de vitesse dont les valeurs sont extrêmement faibles. Une telle
situation peut se produire lorsque l’écoulement à travers une couche de faible
conductivité résulte en de petits vecteurs de vitesse qui peuvent être difficiles à voir.
Pour afficher les flèches d’écoulement selon leur direction, allez dans le menu :
• Direction (à partir de la barre de menu de côté)

Vous pouvez également visualiser les résultats en coupe. Pour cela, sélectionnez une
colonne du modèle.
 18

Exagération verticale de 5

6. SIMULATION D’UN TRAÇAGE DE PARTICULES

Le traçage de particules est un outil généralement utilisé pour définir et visualiser les
voies de migration d'une contamination et pour déterminer l’aire de captage d’un puits de
pompage. Pour cet exercice, vous allez larguer des particules à partir de la zone où le
déversement a un lieu pour

1) Définir les voies de migration du contaminant.

2) Estimer dans combien de temps le contaminant va se rendre à la frontière Sud du

modèle qui correspond à la limite de la propriété.

Étape 1 : Ajouter des particules

Vous devez placer des particules à la limite sud de la zone où la phase flottante en
hydrocarbure a été observée :
 19

Position initiale des particules

Pour ceci, retourner dans le menu Input (en cliquant sur F10-Main Menu et sur Input
dans la barre de menu supérieure). Puis, cliquez sur :

• Particles (menu supérieur)

Allez ensuite dans :

• Add (menu de gauche)

Puis, sur :

• Add Line

Avec votre souris, définissez une ligne le long de la limite inférieure de la zone où du
diesel flottant a été détecté. Lorsque la ligne sera définie, la fenêtre suivante apparaîtra :
 20

Cochez:

• Forward

pour effectuer un traçage de particules avant, c’est à dire de l’amont vers l’aval.

L’option Backward sert à larguer des particules qui se déplaceront dans la direction
opposée au gradient hydraulique, soit de l’aval vers l’amont. Cette option est surtout
utilisée pour définir l’aire de captage d’un puits de pompage. Fixez le nombre de
particules à 9 et ensuite cliquez sur :

• OK

Étape 2 : Exécuter MODPATH

Vous êtes maintenant prêts à simuler le traçage de particules avec MODPATH.

Retournez dans le menu principal de Visual MODFLOW en cliquant sur :

• F10 – Main Menu (barre d’outils inférieure)

Allez ensuite dans le menu :

• Run (menu supérieur)

La fenêtre suivante apparaîtra alors :

Gardez l’option Steady State Flow. Cliquez ensuite sur :

• OK

Cliquez à nouveau sur :

 21

• Run (menu supérieur)

La fenêtre suivante apparaîtra alors:

Cochez :

• MODFLOW 2000 et MODPATH

Ensuite cliquez sur :

• Translate & Run

 22

Visual MODFLOW va exécuter MODFLOW 2000 et ensuite MODPATH. Des crochets

bleus apparaîtront à gauche de Intro, MODFLOW 2000 et MODPATH lorsque la
simulation sera terminée. Après l’apparition de ces crochets bleus, fermez la fenêtre et
procédez par la suite à la visualisation des résultats de la simulation.

Étape 3 : Visualiser les résultats de la simulation de traçage de particules

Allez dans le menu :

• Output (menu supérieur)

Et ensuite dans :

• Maps

et

• Pathlines

La fenêtre suivante apparaîtra alors :

Trajet des
particules
 23

Le trajet des particules est illustré par les lignes rouges. Les flèches rouges indiquent la
distance parcourue durant un intervalle donné. Pour savoir quelle est la durée de cet
intervalle, allez dans :

• Option (menu de gauche)

L’intervalle entre les flèches est indiqué dans la section Time markers. Modifiez
l’intervalle à 30 jours et revenez à la fenêtre principale en cliquant sur :

• OK

Choisissez une colonne du modèle qui passe par un trajet de particules afin de visualiser
les résultats en coupe.

Point de départ des

particules

En comptant le nombre d’intervalle de 30 jours entre le point de départ des particules

(zone du déversement) et la frontière du modèle (limite de la propriété), vous pouvez
évaluer le temps auquel la contamination en diesel va atteindre la limite de la propriété.
Ce temps de calcul correspond au scénario du pire cas, car il ne tient pas compte des
phénomènes de sorption et de biodégradation qui peuvent respectivement ralentir la
migration de la contamination et réduire les dimensions de la zone de contaminée.
 24

7. IMPLANTATION D’UN PUITS DE POMPAGE

Suite aux résultats de la simulation précédente, il a été décidé d’implanter un puits de

pompage afin d’empêcher la contamination de sortir du site. Les étapes à suivre pour
simuler le pompage de l’eau souterraine sont présentées dans les sections suivantes

Étape 1 : Raffinement de la grille de différences finies

Le raffinement de la grille est nécessaire aux endroits où les gradients hydrauliques

deviennent importants localement, comme à proximité d’un puits de pompage.
L’utilisation d’un maillage grossier à pour simuler de telles conditions peut entraîner des
erreurs de calculs importantes. Vous allez donc raffiner le maillage dans le secteur où le
puits de pompage projeté sera implanté, soit à environ 60 m au sud du modèle.

Zone où le
maillage doit
être raffiné

Pour raffiner le maillage, cliquez sur :

• F10 – Main Menu (menu inférieur)

Ensuite sur :

• Input (menu supérieur)

• Grid (menu supérieur)

• Edit Grid (Menu de gauche)

 25

Cliquez finalement sur :

• Edit Row

pour raffiner l’espacement des lignes. La fenêtre suivante apparaîtra :

Sélectionnez :

• Refine by

Et indiquer 3 pour raffiner le maillage par un facteur de 3. Déplacez ensuite le pointeur

de la souris sur la ligne 28 et ensuite sur la ligne 34 afin de définir la zone où
l’espacement entre les lignes sera raffiné.

Maintenant, vous allez raffiner l’espacement entre les colonnes par un facteur de 3. Pour
ceci, cliquez sur :

• Edit Grid (Menu de gauche)

Et ensuite sur :

• Edit column

pour raffiner l’espacement des lignes. La fenêtre suivante apparaîtra :

Sélectionnez :

• Refine by
 26

et indiquez 3 pour raffiner le maillage par un facteur de 3. Déplacez ensuite le pointeur

de la souris sur la colonne 14 et ensuite sur la colonne 28 afin de définir la zone où
l’espacement entre les colonnes sera raffiné.

Le maillage résultant devrait être le suivant :

Étape 2 : Positionnement du puits de pompage

Pour assigner un puits de pompage, allez dans :

• Wells (menu supérieur)

Puis, sur :

• Pumping Wells

Et dans :

• Add Well (menu de gauche)

Positionnez par la suite le pointeur de votre souris dans la cellule correspondant à

l’intersection entre les lignes et colonnes 35 (35 , 35). Cette position est l’emplacement
du puits de pompage. Cliquez sur le bouton de gauche de votre souris pour placer le
puits. La fenêtre suivante apparaîtra alors :
 27

C’est dans cette fenêtre que vous pouvez entrer les propriétés de votre puits comme le
débit de pompage et la position de la crépine. Les paramètres sont les suivants :

• Screen Bottom = 5 m
• Screen Top = 10 m
• Start Time = 0 day
• End Time = 3650 days
• Rate = - 50 m3 /d

Prenez soin d’indiquer un débit négatif à votre puits afin que celui-ci pompe de l’eau.
Autrement, celui-ci agirait comme un puits d’injection. Cliquez ensuite sur :

• OK

Étape 3 : Exécuter MODFLOW et MODPATH

Vous êtes maintenant prêt à simuler l’écoulement de l’eau souterraine en condition de

pompage ainsi que le traçage de particules avec. Retournez dans le menu principal de
Visual MODFLOW en cliquant sur :

• F10 – Main Menu (barre d’outils inférieure)

Allez ensuite dans le menu :

• Run (menu supérieur)

 28

La fenêtre suivante apparaîtra alors :

Gardez l’option Steady State Flow. Cliquez ensuite sur :

• OK

Cliquez à nouveau sur :

• Run (menu supérieur)

La fenêtre suivante apparaîtra :

Cochez :

• MODFLOW 2000 et MODPATH

Et cliquez sur :

• Translate & Run

Après l’apparition de ces crochets bleus, fermez la fenêtre et procédez à la visualisation

des résultats de la simulation.
 29

Étape 4 : Visualisation des résultats

Allez dans le menu

• Output (menu supérieur)

Et ensuite dans :

• Maps

et

• Pathlines

La fenêtre suivante apparaîtra:

 30

Les résultats indiquent que l’opération d’un puits de pompage à un débit de 50 m3 /jour
n’est pas suffisante pour contenir la contamination à l’intérieur des limites du site.
D’autres simulations seraient nécessaires pour optimiser le design du puits.

Il est également possible de visualiser les iso-contours de rabattement. Pour cela, allez
dans :

• Maps

Ensuite :

• Contouring

Et sélectionnez :

• Drawdown

Si vous êtes dans la couche #2, la fenêtre suivante apparaîtra :

 31

Comme pour les équipotentielles, vous pouvez modifier l’intervalle des iso-contours dans
le menu Option.