cole de technologie suprieure
Dpartement de gnie de la construction
MGC-860 HYDRAULIQUE SOUTERRAINE
Introduction la modlisation numrique de lcoulement de leau souterraine avec
Visual MODFLOW 3.0. Exercice de rvision.
Par : Michel Mailloux, M.Sc., ing.stag.
Jean-Sbastien Dub, ing., Ph.D.
1.
OBJECTIFS DE LEXERCICE
Cet exercice est un didacticiel dont le but est de prsenter linterface Visual MODFLOW
partir dun modle reprsentant les conditions hydrogologiques dun site contamin
rel. Vous devrez suivre les instructions dtailles qui sont dcrites dans cet exercice afin
de construire un modle hydrogologique simple partir duquel vous visualiserez les
voies de migration des contaminants.
2.
DESCRIPTION DU PROBLME
Le site est une chaufferie o un important dversement de diesel a eu lieu au printemps
2000. Ceci a caus une contamination importante des sols et de leau souterraine par des
hydrocarbures (diesel). Du diesel en phase flottante a t dtect 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 dpt de
sable fin de 5 m dpaisseur, suivi dun sable moyen de 10 m dpaisseur reposant sur le
roc. Ce dernier est considr impermable. Les conductivits hydrauliques, valu
partir dessais de permabilit (slug test), sont de 3x10-5 m/s pour le sable fin et de 1x10-4
m/s pour le sable moyen. Lcoulement de leau 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.
3.
GNRER UN MODLE DCOULEMENT DE LEAU SOUTERRAINE
tape 1 : Construire la grille du modle
Pour commencer un nouveau modle, cliquez sur :
File
( partir de la barre de menu suprieure)
Et ensuite sur :
New
( partir du menu droulant)
La bote de dialogue Create New Model sera affiche :
Cest partir de cette fentre que vous allez dfinir les dimensions de votre domaine, le
nombre de ranges, de colonnes et de couches ainsi que les units des diffrentes
proprits que vous allez entrer plus tard durant lexercice.
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 fentre, vous pouvez galement importer un fichier qui servira comme
arrire-plan du modle. Pour ceci, cliquez sur :
Import a site map
Ensuite, vous devez spcifier l'emplacement et le nom du fichier de l'arrire-plan (.bmp).
Browse
Allez jusqu' votre rpertoire tudiant [K:\tudiant] et slectionnez le fichier suivant:
Plan_base.bmp
Cliquez sur :
Create
La fentre suivante apparatra :
(107.9 , 97.9)
(44.7 , 47.9)
Cette fentre sert essentiellement gorfrencer le plan de base partir de points o les
coordonnes X-Y sont connues. Ces points sont les coins infrieur gauche (44.7 , 47.9)
et suprieur droit (107.9 , 97.9) de lencadr gris (voir figure ci- haut). Pour gorfrencer
le plan, cliquez sur le coin infrieur gauche de lencadr gris.
apparatra et vous pourrez dfinir les coordonnes de ce point :
La fentre suivante
Effectuez ensuite la mme procdure pour le coin suprieur droit de lencadr gris. Une
fois que vous aurez gorfrenc le plan de base, la fentre suivante apparatra :
Resize region
Fentre de
dimensionneme
nt
Cliquez sur :
Resize region
Et dplacez la fentre de dimensionnement jusquau coin infrieur gauche de lencadr
blanc. Vous aurez alors dfini les dimensions du modle. Cochez ensuite :
Show Grid
La fentre suivante apparatra :
Cliquez finalement sur :
OK
Visual MODFLOW construira alors une grille de diffrences 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 gnration du maillage, MODFLOW a divis uniformment les deux couches
du modle par dfaut. Cependant, la hauteur de ces deux couches doit tre modifie afin
de bien reprsenter lpaisseur des deux units stratigraphiques. Pour ceci, cliquez sur :
View row (voir la portion suprieure gauche de linterface)
Puis, slectionnez nimporte quelle ligne du modle.
apparatre :
La fentre suivante devrait
Cliquez sur :
Edit Grid
Puis sur
Edit Layers
La fentre suivante apparatra :
Cochez ensuite sur Move, slectionnez ensuite la ligne qui spare les deux couches du
modle et dplacez- l jusqu 10 m dlvation. Vous pouvez visualiser les coordonnes
X-Y- Z correspondant la position du pointeur de la souris dans le coin infrieur droit de
linterface.
Coordonnes X,
Y et Z
tape 3 : Dfinition des proprits hydrogologiques des units stratigraphiques
Pour dfinir les proprits hydrogologiques des matriaux, cliquez sur :
Properties
La fentre suivante apparatra alors :
Cette fentre sert assigner les paramtres par dfaut pour toutes les cellules du modle.
Entrez les mmes valeurs que celles de la fentre ci- haut. Ces paramtres correspondent
celles du sable fin. Une fois que vous aurez entr les valeurs, cliquez sur :
OK
Maintenant, vous devez changer les proprits de la couche infrieure (sable moyen), car
elles sont diffrentes de celles par dfaut. Cliquez sur :
View Layer
Puis, afin daller sur la couche 2, cliquez sur :
Next (immdiatement droite de go to )
Vous pouvez vrifier si vous tes bel et bien sur la couche 2 en regardant vers le coin
infrieur gauche de linterface.
Pour assigner des proprits diffrentes pour cette couche, cliquez sur :
Assign
Et ensuite sur :
Window
Dfinissez ensuite un rectangle qui contient lensemble du maillage et ensuite appuyez
sur le bouton de gauche de votre souris. La fentre suivante apparatra :
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
tape 4 : Dfinition des conditions aux limites
Cette section dcrit les diffrentes tapes suivre pour dfinir les conditions aux limites
du modle. Trois types de conditions aux limites seront assignes, soit des charges
imposes, des flux nuls et une limite de recharge. Les limites de charges imposes seront
assignes aux frontires Nord et Sud du modle de manire bien reprsenter le gradient
hydraulique. Les limites flux nul seront assignes sur les frontires Est et Ouest du
modle qui correspondent chacune la position dune ligne dcoulement de leau
souterraine. La base du modle sera galement une limite flux nul, car il est suppos
que le socle rocheux est impermable.
Limite de charge impose = 12 m
Limite de flux nul
(ligne dcoulement)
Limite de flux nul
(ligne dcoulement)
Limite de charge impose = 9.42 m
10
Assignation des limites de charges imposes
Cette section dcrit les tapes requises pour dfinir les conditions aux limites de charges
imposes qui doivent tre assignes aux frontires Nord et Sud du modle. Pour ceci,
allez dans:
Boundaries (barre de menu suprieur)
Slectionnez :
Constant Head
Puis, allez dans :
Assign ( partir de la barre de menu de ct)
Et finalement :
Line
Dfinissez une ligne horizontale le long de la dernire ligne de cellule au bas du modle.
Cette ligne correspond une quipotentielle de 9.42 m. Dplacez le pointeur de la souris
au coin infrieur gauche de la grille. Cliquez sur le centre de la cellule l'aide du bouton
gauche de la souris, dplacez le pointeur au coin infrieur droit de la grille et cliquez sur
le bouton droit de la souris au centre de la cellule. Une srie de cellules deviendra rose et
la bote de dialogue Assign Constant Head apparatra :
Entrez les valeurs suivantes:
Code #: 0
Assign to appropriate layer
Stop time: 3650
Constant Head Start Pt.:
Constant Head End Pt.:
9.42
9.42
La fonction Assign to appropriate layer fait en sorte que les charges imposes ont t
11
assignes seulement sur la couche #2, car la valeur de charge assigne est infrieure au
contact entre le sable fin et le sable moyen. Les cellules ayant une condition de charges
imposes sont de couleur bourgogne.
Maintenant, assignez une frontire de charge constante sur la bordure suprieure du
modle :
Constant Head
Ensuite, allez dans :
Assign ( partir de la barre de menu de ct)
Slectionnez
Line
Dfinissez une ligne horizontale le long de la dernire ligne de cellule au bas du modle.
Cette ligne correspond une quipotentielle de 12.00 m.
Dplacez par la suite le
pointeur de la souris au coin suprieur gauche de la grille. Cliquez sur le centre de la
cellule l'aide du bouton gauche de la souris, dplacez le pointeur au coin suprieur 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.:
Constant Head End Pt.:
12.00
12.00
La srie de cellules roses deviendra alors bourgogne, indiquant qu'une valeur de charge
constante a t dfinie. Cependant, les charges constantes sont assignes 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 ct)
Et ensuite slectionnez :
Layer
12
La fentre suivante apparatra alors :
Cochez ensuite sur
Constant Head et Layer 2
Et cliquez sur
OK
Maintenant, vous allez examiner le modle en coupe transversale pour vrifier si Visual
MODFLOW a bel et bien copi les charges constantes dans la couche infrieure. Pour
ceci, allez dans :
View Column ( partir de la barre de menu de ct)
Dplacez le pointeur de la souris (barre rouge) dans le domaine du modle et slectionnez
n'importe quelle colonne en cliquant sur le bouton gauche de la souris. Vous pouvez
observer que la charge constante a t attribue aux cellules de la deuxime couche
(cellules rouges). Pour voir plus en dtails le modle, vous pouvez appliquer un facteur
d'exagration verticale :
F8 - Vert. Exag. (barre de menu infrieure)
Vertical Exaggeration: 3.5
Assignation de la recharge
La recharge annuelle de la nappe phratique de ce site a t estime 300 mm/an. Avant
dassigner la recharge au modle, assurez-vous dtre dans la couche 1 (vrifiez au-
13
dessous du cube de navigation dans le coin infrieur gauche). Suivez les tapes suivantes
pour spcifier la recharge annuelle:
Boundaries ( partir de la barre de menu suprieure)
Puis,
Recharge
La bote de dialogue Default Recharge apparatra. Entrez linformation suivante dans
les cases appropries:
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 modle si aucune condition frontire (charge impose, recharge etc.) nest
assigne par lutilisateur. Dans le cas de ce modle, les frontires Est et Ouest ainsi que
la base du domaine agiront comme limite flux nul.
4.
SIMULATION DE LCOULEMENT DE LEAU SOUTERRAINE AVEC
MODFLOW
Vous tes maintenant prt simuler lcoulement de leau souterraine avec MODFLOW.
Retournez dans le menu principal de Visual MODFLOW en cliquant sur :
F10 Main Menu (barre doutils infrieure)
Allez ensuite dans le menu :
Run (menu suprieur)
La fentre suivante apparatra alors :
14
Gardez loption Steady State Flow, car vous simulerez lcoulement en rgime
permanent. Cliquez sur :
OK
Et ensuite recliquez sur :
Run (menu suprieur)
La fentre suivante apparatra :
Cochez seulement :
MODFLOW 2000
cette tape, il vous fa ut seulement simuler lcoulement de leau souterraine.
Puis cliquez sur :
Translate & Run
Linterface Visual MODFLOW va par la suite crer les fichiers ncessaires lexcution
du modle MODFLOW. Aprs cela, la bote de dialogue suivante sera affiche :
15
Des crochets bleus apparatront gauche de Intro et MODFLOW 2000 lorsque la
simulation sera termine. Aprs lapparition de ces crochets bleus, fermez la fentre et
procdez par la suite la visualisation des rsultats de la simulation.
16
5.
VISUALISATION DES RSULTATS
Allez dans le menu :
Output (menu suprieur)
La pizomtrie rsultant de votre simulation apparatra 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, lintervalle des contours et le nombre de
dcimales pour les tiquettes de contours. Espacez les iso-contours de 0.25 m et
dfinissez 2 dcimales pour les tiquettes. Cliquez ensuite sur :
OK
pour sortir de ce menu.
17
La zone de couleur olive qui est prsente au bas du modle est une zone de cellule sche.
Ces cellules sont en fait inactives par MODFLOW, car le niveau de la nappe phratique
est infrieur la base de la cellule.
Cellule active
Cellule inactive
ou sche
Cellule active
Cellule active
Visualisation des vecteurs de vitesses
Pour visualiser les vecteurs de vitesses, allez dans :
Map (menu suprieur)
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 flches est diffrente. 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 extrmement faibles. Une telle
situation peut se produire lorsque lcoulement travers une couche de faible
conductivit rsulte en de petits vecteurs de vitesse qui peuvent tre difficiles voir.
Pour afficher les flches dcoulement selon leur direction, allez dans le menu :
Direction ( partir de la barre de menu de ct)
Vous pouvez galement visualiser les rsultats en coupe. Pour cela, slectionnez une
colonne du modle.
18
Exagration verticale de 5
6.
SIMULATION DUN TRAAGE DE PARTICULES
Le traage de particules est un outil gnralement utilis pour dfinir et visualiser les
voies de migration d'une contamination et pour dterminer laire de captage dun puits de
pompage. Pour cet exercice, vous allez larguer des particules partir de la zone o le
dversement a un lieu pour
1) Dfinir les voies de migration du contaminant.
2) Estimer dans combien de temps le contaminant va se rendre la frontire Sud du
modle qui correspond la limite de la proprit.
tape 1 :
Ajouter des particules
Vous devez placer des particules la limite sud de la zone o la phase flottante en
hydrocarbure a t observe :
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 suprieure). Puis, cliquez sur :
Particles (menu suprieur)
Allez ensuite dans :
Add (menu de gauche)
Puis, sur :
Add Line
Avec votre souris, dfinissez une ligne le long de la limite infrieure de la zone o du
diesel flottant a t dtect. Lorsque la ligne sera dfinie, la fentre suivante apparatra :
20
Cochez:
Forward
pour effectuer un traage de particules avant, cest dire de lamont vers laval.
Loption Backward sert larguer des particules qui se dplaceront dans la direction
oppose au gradient hydraulique, soit de laval vers lamont. Cette option est surtout
utilise pour dfinir laire de captage dun puits de pompage. Fixez le nombre de
particules 9 et ensuite cliquez sur :
OK
tape 2 :
Excuter MODPATH
Vous tes maintenant prts simuler le traage de particules avec MODPATH.
Retournez dans le menu principal de Visual MODFLOW en cliquant sur :
F10 Main Menu (barre doutils infrieure)
Allez ensuite dans le menu :
Run (menu suprieur)
La fentre suivante apparatra alors :
Gardez loption Steady State Flow. Cliquez ensuite sur :
OK
Cliquez nouveau sur :
21
Run (menu suprieur)
La fentre suivante apparatra alors:
Cochez :
MODFLOW 2000 et MODPATH
Ensuite cliquez sur :
Translate & Run
22
Visual MODFLOW va excuter MODFLOW 2000 et ensuite MODPATH. Des crochets
bleus apparatront gauche de Intro, MODFLOW 2000 et MODPATH lorsque la
simulation sera termine. Aprs lapparition de ces crochets bleus, fermez la fentre et
procdez par la suite la visualisation des rsultats de la simulation.
tape 3 :
Visualiser les rsultats de la simulation de traage de particules
Allez dans le menu :
Output (menu suprieur)
Et ensuite dans :
Maps
et
Pathlines
La fentre suivante apparatra alors :
Trajet des
particules
23
Le trajet des particules est illustr par les lignes rouges. Les flches rouges indiquent la
distance parcourue durant un intervalle donn. Pour savoir quelle est la dure de cet
intervalle, allez dans :
Option (menu de gauche)
Lintervalle entre les flches est indiqu dans la section Time markers.
lintervalle 30 jours et revenez la fentre principale en cliquant sur :
Modifiez
OK
Choisissez une colonne du modle qui passe par un trajet de particules afin de visualiser
les rsultats en coupe.
Point de dpart des
particules
En comptant le nombre dintervalle de 30 jours entre le point de dpart des particules
(zone du dversement) et la frontire du modle (limite de la proprit), vous pouvez
valuer le temps auquel la contamination en diesel va atteindre la limite de la proprit.
Ce temps de calcul correspond au scnario du pire cas, car il ne tient pas compte des
phnomnes de sorption et de biodgradation qui peuvent respectivement ralentir la
migration de la contamination et rduire les dimensions de la zone de contamine.
24
7.
IMPLANTATION DUN PUITS DE POMPAGE
Suite aux rsultats de la simulation prcdente, il a t dcid dimplanter un puits de
pompage afin dempcher la contamination de sortir du site. Les tapes suivre pour
simuler le pompage de leau souterraine sont prsentes dans les sections suivantes
tape 1 : Raffinement de la grille de diffrences finies
Le raffinement de la grille est ncessaire aux endroits o les gradients hydrauliques
deviennent importants localement, comme proximit dun puits de pompage.
Lutilisation dun maillage grossier pour simuler de telles conditions peut entraner 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 modle.
Zone o le
maillage doit
tre raffin
Pour raffiner le maillage, cliquez sur :
F10 Main Menu (menu infrieur)
Ensuite sur :
Input (menu suprieur)
Grid (menu suprieur)
Edit Grid (Menu de gauche)
25
Cliquez finalement sur :
Edit Row
pour raffiner lespacement des lignes. La fentre suivante apparatra :
Slectionnez :
Refine by
Et indiquer 3 pour raffiner le maillage par un facteur de 3. Dplacez ensuite le pointeur
de la souris sur la ligne 28 et ensuite sur la ligne 34 afin de dfinir la zone o
lespacement entre les lignes sera raffin.
Maintenant, vous allez raffiner lespacement entre les colonnes par un facteur de 3. Pour
ceci, cliquez sur :
Edit Grid (Menu de gauche)
Et ensuite sur :
Edit column
pour raffiner lespacement des lignes. La fentre suivante apparatra :
Slectionnez :
Refine by
26
et indiquez 3 pour raffiner le maillage par un facteur de 3. Dplacez ensuite le pointeur
de la souris sur la colonne 14 et ensuite sur la colonne 28 afin de dfinir la zone o
lespacement entre les colonnes sera raffin.
Le maillage rsultant devrait tre le suivant :
tape 2 : Positionnement du puits de pompage
Pour assigner un puits de pompage, allez dans :
Wells (menu suprieur)
Puis, sur :
Pumping Wells
Et dans :
Add Well (menu de gauche)
Positionnez par la suite le pointeur de votre souris dans la cellule correspondant
lintersection entre les lignes et colonnes 35 (35 , 35). Cette position est lemplacement
du puits de pompage. Cliquez sur le bouton de gauche de votre souris pour placer le
puits. La fentre suivante apparatra alors :
27
Cest dans cette fentre que vous pouvez entrer les proprits de votre puits comme le
dbit de pompage et la position de la crpine. Les paramtres 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 dindiquer un dbit ngatif votre puits afin que celui-ci pompe de leau.
Autrement, celui-ci agirait comme un puits dinjection. Cliquez ensuite sur :
OK
tape 3 :
Excuter MODFLOW et MODPATH
Vous tes maintenant prt simuler lcoulement de leau souterraine en condition de
pompage ainsi que le traage de particules avec. Retournez dans le menu principal de
Visual MODFLOW en cliquant sur :
F10 Main Menu (barre doutils infrieure)
Allez ensuite dans le menu :
Run (menu suprieur)
28
La fentre suivante apparatra alors :
Gardez loption Steady State Flow. Cliquez ensuite sur :
OK
Cliquez nouveau sur :
Run (menu suprieur)
La fentre suivante apparatra :
Cochez :
MODFLOW 2000 et MODPATH
Et cliquez sur :
Translate & Run
Aprs lapparition de ces crochets bleus, fermez la fentre et procdez la visualisation
des rsultats de la simulation.
29
tape 4 :
Visualisation des rsultats
Allez dans le menu
Output (menu suprieur)
Et ensuite dans :
Maps
et
Pathlines
La fentre suivante apparatra:
30
Les rsultats indiquent que lopration dun puits de pompage un dbit de 50 m3 /jour
nest pas suffisante pour contenir la contamination lintrieur des limites du site.
Dautres simulations seraient ncessaires 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 slectionnez :
Drawdown
Si vous tes dans la couche #2, la fentre suivante apparatra :
31
Comme pour les quipotentielles, vous pouvez modifier lintervalle des iso-contours dans
le menu Option.
