Institut Agronomique et Vétérinaire Hassan II

Département du Génie Rural

Initiation à GMS
MODFLOW – L’approche par maillage
Deux approches peuvent être utilisées pour construire des simulations Modflow dans GMS : l’approche
par maillage et l’approche conceptuel. L’approche par maillage implique le travail avec des mailles 3D
en appliquant les termes puits/sources et d’autres modèles bloc par bloc. L’approche conceptuelle
implique l’utilisation d’un outil SIG pour élaborer une carte pour développer le modèle conceptuel du
site à modéliser. Les données sont copiées après dans la grille du maillage.

L’approche par maillage est décrite dans le présent chapitre. Dans la plus part des cas, l’approche
conceptuelle est plus efficace que l’approche par maillage. Cependant, cette dernière est très pratique
pour les problèmes simples ou pour des exercices académiques où l’édition des données par maille est
nécessaire.

1. Description du problème
Le problème à résoudre dans cet exercice est donné dans la figure 1. Trois aquifères vont être simulés
en utilisant trois couches dans une grille conceptuelle de blocs. La grille couvre une région carrée de 75
km sur 75 km. Le maillage consiste en 15 lignes et 15 colonnes, les mailles sont carrées et mesurent
5000 m de coté dans le plan horizontal. Pour un souci de simplicité la topographie de chaque couche est
considéré plate. Les valeurs de conductivité hydraulique sont données dans la direction horizontale. Pour
la direction verticale on utilisera une fraction de la conductivité hydraulique horizontale.

Figure 1: Un problème simple à résoudre

Pour chaque couche les données sont comme suit :

- Couche 1: K = 15 m/j, Cote du niveau du sol = 60 m, cote de substratum = -45 m;

- Couche 2: K = 0.9 m/j, Cote du toit = -45 m, cote de mur = -120 m;

- Couche 3: K = 2 m/j, Cote du toit = -120 m, cote de mur = -215.

L’écoulement dans le système est dû à l’infiltration des précipitations qui sera défini comme une
recharge de la nappe. Les prélèvements du système sont ceux qui se font vers un drain, un puits de
pompage (non visible sur la figure 1) et vers un lac représenté par une limite à charge hydraulique
constante à gauche de la figure. Les charges initiales sont considérées égales à 0 et on cherche à trouver
la solution du régime permanent.

2. Mise en route
 1. Si nécessaire, lancer GMS. Si GMS est déjà lancé, sélectionner File | New dans le menu pour
assurer que la configuration du programme revient à son état initial.

3. Unités
A ce stade, on peut définir les unités utilisées dans le modèle. Les unités choisies seront modifiées en
éditant les champs de l’interface GMS.

1. Sélectionner la commande du menu Edit | Units.

2. Pour Length, enter m (pour mètre). Pour Time, entrer d (pour jours). On ignore pour l’instant
les autres unités (elles ne sont pas utilisées dans les simulations de l’écoulement).
3. Cliquer sur le bouton OK.

4. Création du maillage
La première étape de résolution du problème est de créer un maillage de différences finies en 3D.

1. Dans l’espace Explorateur de Projet (Project Explorer) cliquer par le bouton droit de souris dans
une zone vide puis dans le menu flotant choisir New | 3D Grid.
2. Dans la section intitulée X-dimension, entrer la valeur 22860 pour la longueur (Length value)
et 15 pour le numéro de mailles (Number cells) ;
3. Dans la section intitulée Y-dimension, entrer la valeur 22860 pour la longueur (Length value)
et 15 pour le numéro de mailles (Number cells) ;
4. Dans la section intitulée Z-dimension, entrer 3 pour le nombre de mailles (Number cells) ;

Plus tard, nous allons entrer les valeurs des cotes du toit et mur de chaque couche. Ainsi, l’épaisseur
des mailles dans la direction z introduit à ce niveau ne vont pas affecter les calculs du modèle
MODFLOW.

5. Cliquer sur le bouton OK.

Le maillage va apparaître sur l’écran de l’ordinateur. Une représentation simplifiée du maillage doit
également apparaître dans Mini-Grid Plot.

5. Création d’une Simulation MODFLOW

L’étape suivante dans la mise en œuvre du modèle est d’initialiser la simulation de MODFLOW.

1. Dans l’espace Explorateur de Projet (Project Explorer) cliquer par le bouton droit sur le
dossier 3D Grid Data puis sélectionner New MODFLOW dans le menu flottant.

5.1. Le Package global

Les données d’entrée de MODFLOW sont subdivisées en packages. Certains package sont optionnels
et d’autre sont obligatoires. L’un des packages obligatoire est le Package Global par lequel nous allons
commencer.

En premier lieu nous allons choisir le package.

 2. Cliquer sur le bouton Packages.

Les packages

La boîte de dialogue packages est utilisés pour spécifier quel package va être utilisé pour la mise en
œuvre du modèle. Le package de base est déjà utilisé et ne peut pas être arrêté. Pour choisir d’autres
packages il faut passer par les étapes suivantes:

3. Dans la section Boundary Conditions cocher les options Drain (DRN1) et pour la prise en
compte du drainage, Well (WEL1) pour la prise en compte des pompages.
4. Pour la prise en compte de la recharge choisir l’option Recharge (RCH1).
5. Dans la section Solver, choisir le package Stongly Implicite. Procedure (SIP1)’ pour choisir
la méthode de résolution des équations de différences finies.
6. Cliquer sur OK pour quitter la boîte de dialogue Packages.

Le tableau IBOUND

Dans l’étape suivante nous allons spécifier le tableau IBOUND (variable qui prend des valeurs <0, =0
ou >1 en fonction du type de la condition aux limites). Quand IBOUND>0 la maille est dite active,
quand IBOUND=0 la maille est inactive et quand IBOUND<0 la maille est à charge imposée. Pour notre
problème, toutes les mailles seront actives excepté pour la colonne de gauche des deux premières
couches dans les mailles où la charge hydraulique est imposée.

7. Cliquer sur le bouton IBOUND ;

La boîte de dialogue IBOUND affiche les valeurs du tableau IBOUND pour chaque couche à la fois. Le
champ Layer dans le coin haut à gauche est utilisé pour changer la couche en cours. Pour notre problème,
nous allons spécifier que toutes les valeurs du tableau sont supérieures à zéro, excepté pour la colonne
gauche des couches 1 et 2 où les valeurs sont inférieures à 0. Par défaut, les valeurs du tableau sont
supérieures à zéro. Ainsi, ce dont on a besoin est de changer les valeurs des mailles où la charge
hydraulique est imposée. Ceci peut être accompli en introduisant la valeur -1 pour chaque maille à
charge imposée. Cependant, il existe une autre façon pour éditer le tableau IBOUND qui est plus simple
dans ce cas. cette méthode va être décrite plus tard. Pour l’instant on va garder toutes les mailles active
(IBOUND=1).

8. Cliquer sur bouton OK pour quitter la boîte de dialogue IBOUND.

Charges initiales

Dans l'étape suivante on va spécifier les valeurs des charges initiales.

9. Cliquer sur le bouton Strating Heads.

Le tableau des charges hydrauliques est utilisé pour établir une valeur initiale en exécutant une
simulation en régime transitoire. Puisque nous allons exécuter une simulation en régime permanent, les
charges initiales ne vont pas changer la solution. Cependant, plus les charges initiales sont près des
valeurs de la solution plus MODFLOW converge rapidement. En outre, pour certains types de couches,
si les valeurs principales de départ sont très faible MODEFLOW peut l’interpréter par des mailles
deviennent sèches. Pour le problème qu’on est entrain de traiter des valeurs de la charge hydraulique
égales à 0 sont suffisantes.

Le tableau des charges hydrauliques initiales est aussi utilisé pour spécifier les valeurs de la charge
hydraulique des mailles à charge imposées. Pour notre problème, les charges hydrauliques imposées
sont égales à 0. On a n’a pas donc besoin de changer ces valeurs puisqu’elles sont déjà égale à 0.

10. Cliquer sur OK pour quitter la boîte de dialogue Starting Heads.

Les cotes du toit et du mur de chaque couche

Dans l’étape suivante nous allons spécifier les côtes du toit et du mur de chaque couche.

11. Cliquer sur le bouton Top Elevation.

12. S’assurer que la couche en cours en 1.
13. Cliquer sur le bouton ConstantàLayer.
14. Entrer la valeur 60 and cliquer sur OK.
15. Cliquer sur OK pour quitter la boîte de dialogue Top Elevations.

GMS force les cotes du toit de la couche d’être à la même cote que le mur de la couche en dessus.
Ainsi, nous avons juste besoin d’entrer les cotes des murs des autres couches et les toits vont être
spécifiés automatiquement.

16. Cliquer sur le bouton Bottom Elevation.

17. S’assurer que la couche en cours est 1.
18. Cliquer sur le bouton ConstantàLayer.
19. Entrer la valeur -45 et clique sur OK.
20. Changer la coche Layer à 2.
21. Cliquer sur le bouton ConstantàLayer.
22. Entrer la valeur -120 et cliquer sur OK.
23. Change the Layer to 3.
24. Cliquet sur le bouton ConstantàLayer.
25. Entrer la valeur -215 et clique sur OK.
26. Cliquer sur OK pour sortir de la boîte de dialogue Bottom Elevation.
27. Cliquer sur OK pour sortir de la boîte de dialogue MODFLOW Global Package.

6. Affectation des valeurs à IBOUND directement aux mailles

Comme mentionner plus haut, les valeurs de IBOUND peuvent être entrées à travers la boîte de dialogue
du Tableau IBOUND. Dans certains cas, il est plus simple d’affecter les valeurs directement aux mailles.
Ceci peut être accompli en utilisant la commande propriétés de la maille (Cell Properties). Avant
d’utiliser la commande, nous devons d’abord sélectionner les mailles de la colonne gauche des deux
couches de dessus.

6.1. Affichage de la colonne gauche

Pour simplifier la sélection des mailles, nous allons changer l’affichage de telle manière à voir la couche
gauche du maillage.
1. Cliquer sur le bouton Side View .

Le maillage apparait très fin. Pour faciliter les choses, nous allons amplifier la valeur de Z de telle
manière à faire apparaître les mailles dans la direction Z.
2. Clique sur le bouton Display .
3. Changer Z magnification à 15.
4. Cliquer sur OK.
6.2. Sélection des mailles

Pour sélectionner une maille

1. Cliquer sur l’outil Select Cells .
2. Changer le numéro de colonne à 1 dans l’affichage Mini-Grid.

Il faut noter que nous affichons la colonne 1 à gauche du domaine.

3. Dessiner un rectangle autour de toutes les mailles des deux premières colonnes en haut du
maillage.

6.3. Changement des valeurs de IBOUND

Pour éditer les valeurs de IBOUND:

1. Cliquer par le bouton droit sur l’une des mailles sélectionnées.
2. Choisir Properties dans le menu flottant.
3. Changer l’option IBOUND à Specified head.
4. Cliquer sur OK pour sortir de la boîte de dialogue 3D Cell Properties.
5. Cliquer Plan View .

Noter qu’un symbole est affiché dans les mailles que nous avons éditées, indiquant que ces mailles
sont des mailles à charge hydraulique imposée.

6.4. Vérification des valeurs

Pour s’assurer que les valeurs de IBOUND ont été entrées correctement:
1. Exécuter la commande du menu MODFLOW | Global Options.
2. Cliquer sur le bouton IBOUND.
3. Changer la valeur de la couche (Layer) en haut à gauche de la boîte de dialogue.

Noter que les mailles de la colonne gauche ont toutes une valeur de IBAOUND=-1 dans les deux couches
de dessus. Toutes les données d’entrée de MODFLOW peuvent être éditées dans GMS en utilisant à la
fois le tableau de données ou par sélection d’une série de maille directement.

4. Cliquer sur OK pour quitter la boîte de dialogue du tableau IBOUND.

5. Cliquer sur OK pour quitter la boîte de dialogue MODFLOW Global Package.

7. Le package LPF
La prochaine étape dans la mise en œuvre du modèle est d’entrer les données sur le package des
propriétés de l’écoulement dans chaque couches (Layer Property Flow ou LPF). Le package LPF
calcule la conductance entre les mailles et configure les équations de différences finies décrivant
l’écoulement d’une maille à l’autre. Pour entrer les données LPF:

1. Lancer le menu MODFLOW | LPF Package.

7.1. Types de couches

Les options dans la section Layer Data de la boîte de dialogue sont utilisées pour définir le type de la
couche et la conductivité hydraulique pour chaque couche. Pour notre problème, nous avons trois
couches. La couche supérieure est une nappe phréatique les autres couches sont des nappes captives. La
couche par défaut dans GMS est ‘convertible’ c’est à dire elle peut être captive ou libre. Ainsi nous
n’avons pas besoin de changer le type de la couche.
7.2. Les paramètres des couches

Les boutons dans la section Layer Data sont utilisés pour entrer les paramètres nécessaires pour calculer
les conductances d’une cellule à l’autre. MODFLOW requière un jeu de paramètres pour chaque couche
dépendant du type de la couche.

7.3. La couche supérieure

En premier lieu, nous allons entrer les données pour la couche supérieure:
1. Cliquer sur le bouton Horizontal Hydraulic Conductivity.
2. Clique sur le bouton ConstantàLayer.
3. Entrer la valeur 15.
4. Cliquer sur OK.
5. Cliquer sur OK pour sortir de la boîte de dialogue Horizontal Hydraulic Conductivity.
6. Répéter ce processus pour entrer la valeur 10 pour l’anisotropie verticale.

7.4 Couche du milieu

Dans la suite nous allons introduire les données pour la couche du milieu:
1. Changer la couche à 2 pour le champ layer pour pouvoir éditer les données de la couche 2.
2. Entrer les valeurs suivantes pour la couche 2:

Paramètre Value
Conductivité Hydraulique Horizontal 0.9 m/j
Anisotropie verticale 5

8.5 Couche inférieure

Nous introduisant maintenant les données pour la couche inférieure:

1. Changer la couche layer à 3 et entrer les valeurs suivantes:

Paramètre Value
Conductivité Hydraulique Horizontal 2 m/j
Anisotropie verticale 5

2. Cliquer sur OK pour quitter boîte de dialogue LPF Package.

8. Le package recharge
Dans la suite nous allons entrer les données sur la recharge en utilisant le package recharge. Le package
recharge est utilize pour simuler la recharge d’un aquifèredue à la pluie et l’infiltration. Pour entrer les
donnée sur la recharge:
1. Exécuter le menu MODFLOW | Source/Sink Packages | Recharge Package.
2. Cliquer sur ConstantàArray.
3. Entrer la valeur 0.0009 et cliquer sur OK.
4. Cliquer sur OK pour sortir de la boîte de dialogue Recharge Package.

9. Le package drain
Nous allons maintenant définir la ligne du drain dans la couche supérieure du modèle. Pour définir le
drain, il faut tout d’abord sélectionner les mailles où le drain est localisé et exécuter après la commande
Point Sources/Sinks.

10.1 Sélection des mailles

Le drain est localisé dans la couche supérieure (layer 1). Nous avons besoin de sélectionner les mailles
2 à 10 de la ligne 8. Pour cela :
1. Choisir l’outil sélection des mailles .
2. Noter que quand on déplace le curseur dans la grille, les indices i,j,k de la maille au-dessous du
curseur sont afficher dans la barre d’état en bas de l’écran comme s’est montré dans la figure 2.

Figure 2 : les indices i,j,k sous le curseur

3. Sélectionner les mailles à i=8, j=2, k=1.

4. Appuyer sur la touche Shift pour permettre le mode multi-sélection et sélectionner les mailles
dans les colonnes 3-10 dans la même ligne (Figure3).
 Figure 3 : Les mailles à sélectionner

10.2 Affectation des drains aux mailles

1. Cliquer par le bouton droit sur les mailles sélectionnées et exécuter la commande
Sources/Sinks à partir du menu flottant.
2. Cliquer sur l’onglet Drain.
3. Cliquer sur le bouton New. Ceci ajoute une nouvelle instance de drain à chacune des mailles
sélectionnée.

A ce point, nous devons entrer une cote et une conductance pour le mailles sélectionnée. Les mailles
drains ont toute la même conductance mais la cote n’est pas la même.

4. Entrer les valeurs suivantes pour la cote et pour la conductance des drains:

ID Elevation Conductance
107 0 7430
108 0 7430
109 3 7430
110 6 7430
111 9 7430
112 15 7430
113 20 7430
114 27 7430
115 30 7430

5. Cliquer sur OK.

6. Désélectionner les mailles en cliquant n’importe où sur le graphique.

10. Le package puits

Dans la suite nous allons définir plusieurs puits en sélectionnant les mailles où les puits sont localizes
en utilisant la commande Point Sources/Sinks.

10.1. Le puits de la couche supérieure

La plus parts des puits pompent dans la couche supérieure mais certains d’entre eux pompent dans la
couche du milieu ou la couche inférieure. On va tout d’abord définir les puits dans la couche
supérieure.
1. En maintenant le bouton de la souris et la touche Shift enfoncés, sélectionner les maille
 comme c’est montré dans la figure 4.

Figure 4 : mailles à sélectionner pour les puits

2. Cliquer par le bouton droit dans l’une des mailles sélectionnées et exécuter la commande
Sources/Sinks du menu flotant.
3. Cliquer sur l’onglet Well.
4. Cliquer sur le bouton New.
5. Entrer la valeur -12230 pour tous les puits (une valeur négative signifie que c’est un
pompage).
6. Cliquer sur le bouton OK.
7. Désélectionner les mailles en cliquant sur n’importe quel point de la grille.

10.2. La couche du milieu

1. Cliquer sur la flèche vers le bas dans Mini-Grid Plot.

Pour sélectionner les mailles:

2. Appuyer simultanément sur la touche Shift et la souris et sélectionner les mailles comme c’est
montré dans la figure 5.

Figure 5 : Maille à sélectionner dans la couche du milieu

3. Cliquer sur le bouton droit de la souris dans l’une des mailles sélectionnées et exécuter la
commande Sources/Sinks.
4. Cliquer sur l’onglet Well.
 5. Cliquer sur le bouton New.
6. Entrer la valeur -12230 pour les deux mailles.
7. Cliquer sur OK.
8. Déselectionner les mailles en cliquant sur un espace en dehors de ces deux mailles.

10.3 La couche inférieure

Finalement nous allons définir un seul puits dans la couche inférieure. Pour afficher la couche
inférieure:
1. Cliquer sur la fiche vers le bas de Mini-Grid Plot.
2. Cliquer sur la maille comme c’est montré dans la figure 6.

Figure 6 : Les mailles à sélectionner dans la couche inférieure

3. Cliquer sur le bouton droit dans la maille sélectionnée et execute la commande du menu
Sources/Sinks.
4. Cliquer sur l’onglet Well.
5. Cliquer sur New.
6. Entrer une valeur de -0.15.
7. Clique sur OK.
8. Désélectionner la maille.

A ce stade tous les puits ont été définis et on peut revenir à la première couche.

9. Clique sur la flèche vers le haut de Mini-Grid Plot.

11. Vérification de la Simulation

Nous avons maintenant complètement défini les données de MODFLOW et nous sommes près pour
lancer une simulation. Cependant, avant la sauvegarde de la simulation et l’exécution de MODFLOW,
nous devons exécuter le vérificateur du modèle MODFLOW et vérifier les erreurs. En raison de la
quantité importante de données nécessitées par les simulations de MODFLOW, il est souvent plus facile
d’oublier certaines données requises ou de trouver une inconsistance et une incompatibilité des options
ou des paramètres. De telles erreurs peuvent causer un plantage de MODFLOW ou la production de
résultats erronés. L’objectif du vérificateur du modèle (Model Checker) est d’analyser les données
d’entrée qui sont actuellement définies pour la simulation par MODFLOW et faire un rapport sur les
problèmes potentiels. Il faut noter que l’exécution avec succès du vérificateur du modèle ne garantit pas
que la solution soit correcte. Il sert simplement à une vérification initiale des données d’entrée pouvant
économiser un temps considérable.
Pour exécuter le vérificateur du modèle (Model Checker):
1. Exécuter la commande du menu MODFLOW | Check Simulation.
2. Cliquer sur le bouton Run Check.

Une liste de messages est affichée pour chaque package d’entrée de données de MODFLOW. Si tout a
été fait correctement, il n’y aura aucune erreur. Quand une erreur existe, si vous sélectionner une erreur,
GMS sélecte la maille ou la couche associé au problème.

3. Cliquer sur le bouton Done pour quitter le vérificateur du modèle.

12 Sauvegarde de la simulation
Now we are ready to save the simulation and run MODFLOW.
1. Exécuter la commande du menu File | Save.
2. Choisir le repertoire.
3. Sauver le projet sous le nom gridmod.gpr.

13 Exécution de MODFLOW
1. Exécuter la commande du menu MODFLOW | Run MODFLOW.

A ce stade MODFLOW est lancé dans une nouvelle fenêtre. Le fichier de donnée est passé à
MODFLOW comme une ligne de commande DOS. MODFLOW ouvre le fichier et commence la
simulation. Quand la solution commence, on voit du texte défiler indiquant que la simulation est en
progression.

2. Quand MODFLOW termine la simulation cliquer sur le bouton Close.

14 Affichage de la solution
GMS lit automatiquement la solution quand on ferme la fenêtre MODFLOW. A ce stade on peut voir
les courbes équipotentielles des charges hydrauliques pour la couche supérieure. On peut aussi voir que
quelques mailles contiennent un triangle bleu. Ces mailles sont “innodée”, c’est à dire que la charge
hydraulique simulée est au dessus du toit de la couche.

14.1. Changement de la couche

Pour afficher la solution de la couche du milieu:

1. Cliquer sur la flèche vers le bas dans Mini-Grid Plot.

Pour afficher la solution de la couche du milieu:

2. Cliquer sur la flèche vers le bas dans Mini-Grid Plot.

Pour retourner à la couche supérieure:

3. Cliquer deux fois sur la couche vers le haut.

14.2 Remplissage des courbes de niveau

On peut afficher les courbes de niveau en utilisant l’option couleur de remplissage (color fill).
1. Exécuter la commande du menu Data | Contour Options.
2. Changer Contour Method à Color Fill.
3. Clique sur OK.
14.3 Couleur de la légende

1. Exécuter la commande du menu Data | Color Ramp Options.

2. Cocher l’option Legend.
3. Clique sur OK.

15 Zone du bilan
La zone du bilan (Zone Budget) est un programme développé par l’USGR (Harbaugh 1990) qui est
utilisé pour calculer un bilan hydrique sous-régional pour MODFLOW. GMS a incorporé un modèle de
bilan similaire. Dans GMS, l’utilisateur définit des zones en affectant à leurs mailles un identifiant (Zone
Budget ID). Une fois ces zones définies, un rapport peut être généré et qui affiche le bilan dans la zone.
Le rapport inclut également une composante qui montre les entrées et les sorties dans les zones
adjacentes.

15.1 Délimitation des zones de bilan

Dans ce modèle nous allons considérer que chaque couche est une zone.
1. Cliquer sur l’outil de sélection des mailles .
2. Si nécessaire changer l’affichage en Plan View .
3. S’assurer qu’on est entrain d’afficher la deuxième couche. Ajuster la couche dans Mini-Grid
Plot sur 2.
4. Dessiner un rectangle autour de toutes les mailles de la couche 2.
5. Cliquer par le bouton droit sur les mailles sélectionnées.
6. Sélectionner la commande Properties dans le menu flottant.
7. Entrer 2 pour Zone budget ID et cliquer sur OK.
8. Changer la couche à 3 dans Mini-Grid Plot.
9. Répéter les étapes de 4 à 7 sauf qu’il faut introduire 3 pour Zone budget ID.

15.2 Affichage du rapport de la zone du bilan

1. Exécuter la commande du menu Data | Flow Budget.

2. Cliquer sur l’onglet Zones.

On est entrain d’afficher le rapport de la première zone (couche supérieure). Le rapport est divisé en
deux section : écoulements à l’intérieur de la zone et écoulements vers l’extérieur de la zone.Chaque
terme puits/sources pris en compte dans le modèle est listé dans le rapport. Les échanges de flux entre
zones sont également affichés.

1. Dans la Zone drop-down box choisir 2.

On peut aussi afficher le rapport de la zone 3. Après avoir fini clique sur OK.
MODFLOW – Interpolation des données
d’une Couche

Pour des sites avec une stratigraphie complexe et avec un écoulement tridimensionnel, un modèle multi-
couches peut être plus précis qu’un modèle à une couche bidimensionnel. Pendant la création d’un
modèle multi-couches, la définition des données du modèle est un vrai défi. GMS contient une suite
d’outils pour interpoler et manipuler les cotes topographiques des différentes couches. Ce chapitre décrit
l’utilisation des points dispersés, de l’interpolation et de MODFLOW.

1. Modules/Interfaces Requis

On aura besoin des composants suivants pour compléter le présent excercice :

- Grid
- Geostatistics
- Map
On peut savoir si ces modules sont disponibles en exécutant la commande du menu File | Register.

2. Mise en route
1. Si nécessaire lancer GMS. Si GMS est déjà lancé select the File | New pour assurer que la
configuration du programme soit à son état par défaut.

3 Interpolation dans une couche MODFLOW

Une des façons de créer les cotes du mur et du toit de chaque couche de MODFLOW est de les interpoler
à partir d’un certain nombre de points dispersés. L’interpolation peut être réalisée directement en
utilisant la commande du menu Interpolation | Interpolate _ 2D Grid lorsque on a de points dispersés
en 2D. Les points peuvent être importés à partir d’une table de valeur dans un fichier. Pour des régions
ayant une stratigraphie complexe les points doivent être définis ace plus d’attention. L’interpolation peut
être vérifiée par le vérificateur de modèle (Model Checker).

4 Exemple de problèmes
Pour illustrerle processus d’interpolation des cotes et de correction des erreurs, on verra une série
d’exemples de problèmes. Chaque exemple illustre un problème différent et décrit une approche simple
pour une modélisation correcte de stratigraphie.

5 1er cas – Couche Complète

Le premier cas qu’on va examiner est montré dans la Figure 7. C’est le cas plus simple dans ce chapitre.
Le site à trois couches et toutes les trois couches s’étendent sur tout le domaine du modèle.

Figure 7 : Coupe Est/Ouest de 1er cas

5.1 Importing the Scatter Point Sets

La première étape dans la définition des données d’une couche est de créer une grille. L’étape suivante
est de créer une série de points à différentes localisations x,y dans le modèle. Chaque point a une cote
pour le toit et le mur pour chaque couche. Dans un problème réel ces données peuvent provenir des
forages d’exploration. Les données sont entrée une un tableau en fichier texte, en utilisant par exemple
Excel. On import, ensuite la table de données dans GMS. Quatre différentes séries de points dispersés
ont été préparées préalablement en utilisant la technique décrite plus haut. Le fichier a été déjà importé
dans GMS et sauver dans un fichier projet GMS. Pour lire ce projet :
1. Exécuter la commande Open .
2. Localiser et ouvrir le dossier C:\SAED\GMS\.
3. Choisir le ficher nommé points.gpr puis clique sur Open.

5.2 Changement en vue de face (Front View)

Avant d’interpoler, on va changer la vue de telle manière à avoir une section transversale de la grille.
De cette façon on va immédiatement voir le résultat de l’interpolation. Avant de change de vue, on va
sélectionner une maille à l’intérieur du domaine. Si une maille est sélectionner quand on chage de vue,
le nouvelle vue changera à la ligne ou à la colonne passant à travers la maille sélectionner.
1. Cliquer sur l’outil de sélection des mailles .
2. Sélectionner une maille près du centre de la grille.
3. Clique sur le bouton Front View .

5.3 Interpolating the Elevation Values

Les points dispersés ont chacun 4 données sur la cote topographique, top1, bot1, bot2, et bot3. L’étape
suivante consiste à interpoler chacune de ces donnée dans la couche appropriée dans le maillage de
MODFLOW. Avant d’interpoler ces points il faut initialiser les données de MODFLOW.
1. Exécuter la commande du menu MODFLOW | New Simulation.
2. Cliquer sur OK.

Pour interpoler les valeurs de des cotes:

3. Pour afficher les séries de données, dans Project Explorer, étendre Case 1 scatter point si
nécessaire en cliquant sur le signe +.
4. Cliquer par le bouton droit sur Cas 1 dans Project Explorer et exécuter Interpolate To |
MODFLOW layers dans le menu flottant.

La boîte de dialogue qui apparaît est utilisée pour définir quelles données des points dispersés sont
interpolées. Les données à interpoler sont listées en haut à gauche de la boîte de dialogue, et les variables
de MODFLOW dans la partie haute droite. Il faut donc chercher à faire correspondre les données au
variables en utilisant les boutons Map et Unmap. GMS essaie de faire ces correspondance d’une manière
automatique en se basant sur les noms des données. Dans le présent cas, toutes les relations ont été
correctement affectées et on peut continuer.

5. Cliquer sur OK.

Noter que les couches interpolées représentent correctement la section transversale de la Figure 7. On
peut voir l’ensemble des sections verticales en appuyant sur les couches Mini-Grid Display.

6 2ème Cas – Couche encastrée

The next case we will examine is illustrated in Figure 6-1. In this case, the middle layer is an embedded
seam that only exists on the left (West) side of the model. This type of layer is more difficult to model
with MODFLOW since a full array of K values must be defined for each layer. Making the cells in the
second layer inactive on the right side of the model will not work since this will result in a no-flow
boundary between the first and third layers. To model the seam correctly, we will use a three layer model
and the second layer will extend throughout the entire model domain. However, we will adjust the layer
thickness of the second layer so that the cells will be very small on the right side of the model, thus
simulating the pinchout.
MODFLOW – Approche du modèle
conceptual
1. Introduction

L’approche du modèle conceptuel consiste en l’utilisation du SIG pour développer un model conceptuel
sur le site à modéliser. La localisation des puits/sources, les paramètres des couches comme la
conductivité hydraulique, les limites du domaine et toute autre donnée nécessaire pour la simulation
peuvent être définis au niveau du modèle conceptuel. Une fois ce modèle est complété, le maillage est
généré et le modèle conceptuel est converti à modèle par maillage et les toutes les affectations sont faites
automatique maille par maille.

2 Description du Problème
Le problème qu’on va résoudre est illustrer dans Figure 8.

Figure 8 : Site modélisé. (a) Vue en pla ; (b) coupe Nord-Sud

On va modéliser l’écoulement souterrain dans les sédiments de vallée limitée par des collines au Nord
et de deux rivières au Sud. Une coupe Nord-Sud du site est représentée dans la Figure 8. Le substratum
imperméable (mur) est composé de calcaire qui affleure au Nord. Le site présente deux couches de
sédiments. La couche supérieure va être modélisée comme une nappe phréatique et la couche inférieure
sera considérée comme un aquifère captif.

La limite du Nord sera considérée à flux nul et les autres limites seront considérées à potentiel imposé
correspondant à la cote du lit des rivières. On supposera également la recharge du système se fait par
infiltration des eaux de pluie. Les ruisseaux à l’intérieur du domaine sont parfois secs mais de temps en
temps ils s’écoulent en raison du drainage de la nappe libre. On représentera ces ruisseaux par des drains.
Il existe aussi deux puits de pompage qu’on va introduire dans le modèle.

3 Mise en route
Si nécessaire lancer GMS. Si GMS est déjà lance, exécuter la commande du menu File | New pour
assurer que GMS restaure sa configuration initiale.

4 Importing the Background Image

La première étape dans mise en œuvre de la simulation est d’importer une image numérique du site à
modéliser. L’image a déjà été importée dans GMS et enregistrée dans un fichier de projet GMS. Pour
lire l’image on va ouvrir le fichier projet. Une fois l’image est importée dans GMS, elle peut être affichée
en arrière plan pour aidé à la conception du modèle.

4.1 Reading the Image

Pour importer l’image:

1. Cliquer sur le bouton Open de la barre d’outils.
2. Localiser et ouvrir le répertoire C:\SAED\GMS.
3. Ouvrir le fichier start.gpr.

Tous les autres objets de GMS sont dessinés sur l’image. L’image apparaît uniquement en vue plan.

5. Enregister le projet
Avant de faire des changements, il faut enregistrer le projet sous un autre nom.
1. Exécuter la commande du menu File | Save As.
2. Enregistrer le projet sous le nom conceptuel.gpr.

On peut maintenant cliquer à chaque fois sur le bouton Save ou le racourcis clavier ctrl s pour
enregistrer le projet à chaque instant.

6 Définition des unités

At this point, we can also define the units used in the conceptual model. The units we choose will be
applied to edit fields in the GMS interface to remind us of the proper units for each parameter.
1. Exécuter la commande du menu Edit | Units.
2. Pour Length choisir m (pour mètre). Pour Time, choisir d (pour jour). On ignore les autres unités
puisqu’elles ne sont pas utilisées dans cet exercice.
3. Cliquer sur OK.

7 Defining the Boundary

La première étape consiste à définir la limite extérieure du domaine. On fera cela par la création d’un
arc qui va former une limite fermée autour du site.

7.1 Create the Coverage

1. Dans Project Explorer cliquer par le bouton droit dans un espace vide puis dans le menu flottant
exécuter le menu New | Conceptual Model.
2. Pour Name, entrer ExempleMC. Pour le modèle choisir MODFLOW.
3. Cliquer sur OK.
4. Clique par le bouton droit sur le modèle conceptuel ExempleMC et exécuter New
Coverage (nouvelle couverture) dans le menu flottant.
 5. Changer Coverage name à Limites. Changer Default elevation à 213. Changer Default layer
range pour aller de 1 à 2.
6. Cliquer sur OK.
7.2. Création de l’Arc
1. Cliquer sur l’outil Create Arc .
2. Commencer le tracer d’une polyligne à gauche (ouest) du domaine dans la localisation montrée
dans Figure 9.

Figure 9 : Création de la limite externe du domaine

3. Créer une polyligne en allant vet l’Est et en commençant par le Sud. Ne pas faire attention à
l’espacement ou à la localisation exacte des points. Il faut juste un nombre suffisant de points
pour définir une localisation approximative de la limite. Les limites Sud et Est coïncide avec les
rivières. La limite Nord coïncide avec la limite calcaire au Nord du Site.
4. Pour terminer la polyligne, cliquer sur le point ou on a commencé.

8 Construction d’une couverture local Puits/Sources

L’étape suivante consiste en la construction d’une couverture locale puits/sources. Cette couverture
définit la limite de la région modélisée ainsi que les termes puits/sources comprenant les puits, les
rivières, les drains et les charges hydrauliques générales.

Les proprieties qui peuvent être assignees aux caractéristiques des objets dans la couverture dependent
du modèle et des options qu’on peut configurer dans la boîte de dialogue Coverage Setup. Avant de
créer les caractéristiques des objets (feature objects), on va changer les options dans la boîte de dialogue
Coverage Setup.
1. Cliquer par le bouton droit sur la couverture Limite et choisir la commande Duplicate dans le
menu flottant. Changer le nom de la nouvelle couverture en Puits & Sources (cela peut
être fait en cliquant par le bouton droit de la souris sur le nouvelle couverture et en choisissant
Rename dans le menu flottant).
2. Cliquer par le bouton droit sur la couverture Puits & Sources et exécuter la commande
Coverage Setup dans le menu flottant.
3. Dans la liste Sources/Sinks/BCs, cocher les options qu’on va utiliser dans cette exemple :
- Layer range (Etendue de la couche).
- Wells (Puits).
- Refine points (affiner les points).
- Specific head (charge hydraulique imposée)
- Drain
 4. S’assurer que la case Use to define model boundary (active area) est cochée.
5. Cliquer sur OK.

8.1 Définition de la limite à potentiel imposé

L’étape suivante consiste à définir la limite à charge imposée (specific head) le lond des limites Sud et
Est du modèle. Avant cela, on doit tout d’abord couper la polyligne en trois autres polylignes. Une
polyligne va définir la limite à flux nul au Nord et les deux autres limites vont définir les deux rivières
auSud et à l’Est. Une polyligne est coupée en sélectionnant un ou plusieurs vertex dans la plolyligne et
en convertissant la les vertex en nœuds.
1. Clique sur l’outil Select Vertices .
2. Clique sur deux vertex comme c’est indiqué dans la Figure 10. Le vertex 1 est localisé à la
jonction entre les deux rivières. Vertex 2 est localisé en haut de la rivière Est. Pour sélectionner
les deux vertex en une seule fois, sélectionner le premier puis appuyer sur la touche Majuscule
(Shift) tout en cliquant sur l’autre vertex.

Figure 10 : Conversion d’un vertex en un nœud.

3. Cliquer par le bouton droit sur l’une des vertex et exécuter la commande Vertex --> Node dans
le menu flottant.

Maintenant que nous avons les trois polyligne, nous allons spécifier les deux polyligne des rivières
comme des limites à charge imposée.

4. Clique sur l’outil Select Arcs .

5. Cliquer sur les polylignes Sud et Est en sélectionnant la première polyligne et en cliquant sur la
deuxième tout en appuyant sur la touche Majuscule (Shift) du clavier.
6. Cliquer par le bouton droit sur l’une des polyligne sélectionnées puis exécuter Attribute Table
dans le menu flottant.
7. Chercher la ligne All et la colonne Type du tableau. Dans cette cellule choisir spec. Head.
Ceci va changer le type pour les deux polyligne.
8. Cliquer sur OK.
9. Clique quelque part sur l’image pour désélectionner les polylignes.

Noter que la couleur des deux polyligne a change indiquant le type de condition aux limites.
L’étape suivante consiste à définer les charges hydrauliques aux bouts des deux polylignes. La charge
le long des deux polylignes sera supposée variée linéairement le long de la polyligne.
10. Cliquer sur l’outil Select Points/Nodes .
11. Double cliquer sur le noeud à l’ouest du domaine au bout de la limite Sud.
12. Entrer la valeur 212 pour Head-Stage.
13. Cliquer sur OK.
14. D’une manière similairen affecter la valeur 208 au noeud de la jonction des deux rivières et la
valeur 214 au neoud en haut de la limite Est.

8.2 Defining the Drain Arcs

A ce stade, on va entrer les polylignes définissant les ruisseaux à l’intérieur du domaine constituant des
drains.
1. Cliquer sur l’outil Create Arc .
2. Créer la plolilygne nommée arc 1 dans Figure 11. Commencer par cliquer sur la polyligne Sud,
créer la polyligne en cliquant sur des points le long du ruisseau et terminer la polyligne en double
cliquant sur la polyligne Nord.

Figure 11 : Les polylignes représentant les drains

Noter que lorsqu’on clique sur les environs d’un vertex dans une polyligne existante, GMS fait
intercepter automatiquement la polyligne crée à la polyligne existante et crée un nœud à la jonction des
deux polylignes.

3. Créer les polylignes arc 1 et arc 2 comme c’est montré dans Figure 11 et de la même manière
que la plyligne arc 1.

Dans la suite, on va définir les polylignes comme des drains et on va leur affecter des conductances et
des cotes topographiques.
4. Cliquer sur l’outil Select Arcs .
5. Sélectionner l’ensemble des polylignes correspondant aux drains en cliquant dessus tous en
maintenant la touche Shift enfoncé.
6. Cliquer par le bouton droit sur l’une des polylignes sélectionnées et exécuter la commande du
menu Attribute Table dans le menu flottant.
7. Dans la ligne All et la colonne Type choisir l’option drain.
 8. Entrer une conductance de 555 dans la ligne All et la colonne Cond.. Ceci représente une
conductance par unité de logueur. GMS calcule automatiquement calcul automatiquement la
conductance dans les mailles du domaine quant la les drains sont affectés à au maillage du
modèle.
9. Changer les propriétés From layer et To layer à 1 pour chaque polygne. Cela signifie sont
uniquement affectés à la couche 1 du modèle.
10. Cliquer sur OK.

Les cotes des drains sont spécifiées dans les extrémités des polylignes. Elles sont supposées variées
linéairement le long de chaque polyligne.

11. Cliquer sur l’outil Select Points/Nodes .

12. Double cliquer sur le noeud 2 in dans la Figure 12. Noter que ce nœud a 2 propriétés qui lui sont
associées puisqu’il est attacher à deux limites.

Figure 12 : Cotes des nœuds² des drains

13. Entrer 212 pour Bot. elev. De la propriété drain. Ne rien changer dans la propriété spec. head.
Cliquer OK.
14. Répéter cette procédure pour affecter les cotes des noeuds du bout de chaque drain en utilisant
les valeurs de la Figure 12.

8.3 Building the polygons

Avec le type de couverture puits/sources locale, l’ensemble de la région à modéliser doit être couverte
par polygones qui ne se chevauchent pas. Ceci définir la région active du maillage du modèle. Dans la
plus parts des cas, tous les polygones vont être des polygones à charge variable. Cependant, d’autres
polygones peuvent être utilisés. Par exemple, pour modéliser un la, un polygone de charge général peut
être utilisé. La façon la plus simple de définir les polygones est de créer en premier lieu toutes les
polylignes utilisées dans la couverture et utiliser ensuite la commande construction de polygones (Build
polygons). Cette commande cherche à travers les polylignes et crée un polygone pour chaque boucle
fermée définie par les polylignes. Ces polygones sont par défaut de type « NONE » mai peuvent être
converti à d’autres types en sélectionnant le polygone et en utilisant la commande propriétés
(Properties).

Maintenant que toutes les polylignes ont été crées dans la couverture, on est prsè pour construire les
polygons. Tous les polygones auront une charge hydraulique variable.

1. Dans le menu principal, exécuter la commande Feature Objects | Build Polygons.

Noter que le polygone n’est pas rempli. On peut changer l’affichage du polygone si l’on veut en
exécutant la commande Display | Display Options.

8.4 Création de puits

L’étape finale dans la création de la couverture du terme Puits/Sources local est de définir les puits de
pompage. Les puits sont définis comme des types d’objete point. Deux puits vont être créés.
1. Cliquer sur l’outil Create Point .
2. Déplacer le curseur la localisation approximative Puits 1 montré dans la Figure 8 et cliquer par
la souri pour créer le point.
3. Quand le nouveau point est sélectionner introduire les coordonnées (835, 1425) dans les champs
d’édition X et Y en haut de la fenêtre GMS.
4. Cliquer sur le bouton Properties .
5. Pour Type, choisir l’option well.
6. Pour Flow rate, entrer la valeur constante -680.
7. Changer les proprieties From layer et To layer pour être égales à 1. Cela signifie que le puits 1
pompe dans la couche 1 uniquement.
8. Clique sur OK.
9. D’une manière similaire, créer l’autre puits à la localisation (3220, 1000) et lui affecter un débit
de pompage de –2830. Cependant, pour ce puits, changer From layer et To layer à 2 pour que
ce puits soit affecté à la couche 2.

Affinement du maillage

Un puits de pompage représente un point de convergence des écoulements souterrains et cause un

gradient hydraulique très raide à son voisinage. Dans le but modéliser d’une manière précise
l’écoulement près du uits, le maillage est affiné au voisinage du puits. Ce type d’affinage peut être réalisé
automatiquement par GMS en affectant des données d’affinage directement au puits dans le modèle
conceptuel.
1. Cliquer sur l’outil Select Points/Nodes .
2. Sélectionner les deux puits e cliquant sur la souris pendant que la touche Shift du clavier est
enfoncée.
3. cliquer sur le bouton Properties .
4. Chercher la colonne Refine et dans la ligne All cocher la case. Ceci permettra au maillage d’être
affiner pour au voisinage des deux puits.
5. Changer la taille de base (Base size) à 25, Bias à 1.1 et Max size à 150 pour les deux points.
6. Cliquer sur OK.

9 Délimitation des zones de recharge

L’étape suivante dans la conception du modèle conceptuel est de construire une couverture qui définie
les zones de recharge. On va supposer que la recharge sur l’ensemble de la région à modéliser est
uniforme excepter pour une décharge situé à l’intérieur de cette région. La recharge dans la surface de
la décharge sera réduite en raison d’un système d’étanchéisation.

9.1 Copie de la limite

On va créer une couverture de recharge en copiant la limite :
1. Cliquer par le bouton droit de la souris sur la couverture Limite et exécuter la commande
Duplicate dans le menu folottant.
2. Change le nom de la couverture de new coverage à Recharge.
3. Cliquer par le bouton droit sur la couverture Recharge et exécuter la commande Coverage Setup
 dans el menu flottant.
4. Dans la liste Areal Propertiest, cocher la propriété Recharge rate.
5. Cliquer sur OK.

9.2 Création de la limite de la décharge

Dans la suite on va créer une polyligne délimitant la décharge. Pour cela, on va tout d’abord créer une
boucle fermée de points sous forme de rectangle dans la zone de la décharge (Figure 8). On va ensuite
éditer les nœuds et les vertex de telle façon à ce que la polyligne coïcide exactement avec les limite de
la décharge.

1. Cliquer sur l’outil Create Arc tool .

2. Crée un polygone rectangulaire dans une zone proche de la localization de la décharge comme
c’est montré dans la Figure 13.

Figure 13 : délimitation de la décharge

Maintenant que la polyligne est créée dans une localization approximative, on va éditer les coordonnées
des vertex et des noeuds pour definer des coordonnées précises.
3. Cliquer sur l’outil Select Vertices .
4. Dessiner un rectangle autour de l’ensemble de la décharge ce qui permettra à tous lex vertex et
les nœuds d’être sélectionnés.
5. Clique par le bouton droit sur l’un des vertex sélectionné puis exécuter la commande Vertex ->
Node du menu flottant.
6. Clique sur l’outil Select Points/Nodes .
7. Cliquer sur l’un des noeuds des coins du rectangle.
8. Pendant que le nœud est sélectionner, entrer les coordonnées exactes du nœud dans les champs
d’édition X et Y en haut de la fenêtre GMS.
9. Répéter ce processus pour les 3 nœuds restant.

9.3 Construction de polygones

Maintenant que les polylignes sont définies on peut construire les polygones.
Now that the arcs are defined, we can build the polygons.
1. Dans le menu principal exécuter la commande Feature Objects | Build Polygons.

9.4 Affectation des valeurs de la recharge

Les zones de recharge étant définies, on peut affecter les valeurs de la recharge. On va affecter une
valeur de la recharge à la décharge et une autre valeur au reste du polygone.
1. Clique sur l’outil Select Polygons tool .
 2. Double cliquer sur le polygone de la décharge.
3. Changer Recharge rate 0.00006.

NB: Le taux de recharge rate est faible que celui affecté à d’autre polygone. La décharge sera revêtue et
aura un impact faible sur la recharge.

4. Clique sur OK.

5. Double cliquer sur l’autre polygone.
6. Changer Recharge rate à 0.00695.
7. Cliquer sur OK.

10 Definition de la conductivité Hydraulique

Dans la suite on va entrer la conductivité hydraulique pour chaque couche. Dans beaucoup de cas, on
souhaite définir plusieurs polygones correspondant à différentes zones de conductivité hydraulique.
Dans un souci de simplification on utilisera une seule valeur par couche.

10.1 Copying the Boundary

On va créer une autre couverture en copiant la couverture Limite.

1. Cliquer par le bouton droit sur la couverture Limite et exécuter la commande Duplicate du
menu dans le menu flottant.
2. Changer le nom new coverage à Couche 1.
3. Cliquer par le bouton droit sur la couverture Couche 1 et exécuter la commande Coverage
Setup dans le menu flottant.
4. Dans la liste Areal Properties, cocher les options suivantes :
- Horizontal K.
- Vertical anis.
5. Changer Default layer range de 1 à 1.
6. Cliquer OK.
7. Cliquer par le bouton droit sur la couverture Couche 1 puis exécuter la commande
Duplicate du menu flottant. Changer le nom de la couverture à Couche 2.
8. Cliquer par le bouton droit sur la couverture Couche 2 et exécuter la commande Coverage
Setup du menu flottant.
9. Changer Default layer range de 2 à 2.
10. Cliquer OK.

10.2 Couche supérieure

Tout d’abord, on va affecter les valeurs à K pour la couche supérieure.
11. Cliquer sur la couverture Couche 1 dans Project Explorer.
12. Dans le menu principal exécuter la commande Feature Objects | Build Polygons.
13. Cliquer sur l’outil Select Polygons , double cliquer sur le polygone.
14. Changer la valeur de Horizonal K à 5.5.
15. Change la valeur de Vertical anis. à 4.
16. Cliquer OK.

10.3 Couche inférieure

Pour la couche inférieure :

1. Cliquer sur la couverture Couche 2 dans Project Explorer.
2. Dans le menu principal exécuter la commande Feature Objects | Build Polygons.
3. Cliquer sur l’outil Select Polygons , double cliquer sur le polygone.
4. Changer la valeur de Horizonal K à 10.
 5. Change la valeur de Vertical anis. à 4.
6. Cliquer OK.

Ceci complète la définition des couvertures dans le modèle conceptuel. Avant de continuer la création
du maillage, on va faire de Puits & Sources la couverture active.

7. Cliquer sur la couverture Puits & Sources dans Project Explorer.

11 Localisation du cadre du maillage

L’ensemble des couvertures du modèle conceptuel étant comlétées, on est près pour créer un le maillage.
La premi²ère étape dans la création du maillage est de définir la localisation et l’orientation du maillage
en utilisant le cadre du maillage. Celui-ci représente la limite externe du maillage. Il peut être posisionner
au-dessus de carte géographique du site à modéliser.

1. Dans Project Explorer cliquer par le bouton droit sur un espace vide et exécuter la commande
New | Grid Frame dans le menu flottant.
2. Dans Project Explorer cliquer par le bouton droit sur Grid Frame et exécuter la commande Fit
to Active Coverage dans le menu flottant.
3. Double-cliquer sur Grid frame dans Project Explorer pour lancer la boîte de dialogue propriétés.
4. Changer Origin z: à 170 et Dimension z: à 60. Ceci fournit une série de valeurs initiales pour
les valeurs de cotes des couches du modèle MODFLOW. Plus tard, on va interpoler les valeurs
des cotes des différentes couches.
5. Clique sur OK pour sortie de la boîte de dialogue Grid Frame.

12 Creating the Grid

Maintenant que les couvertures et le cadre du maillage sont créés, on est près pour créer le maillage.
1. Dans le menu principale de GMS, executer la commande Feature Objects | Mapà 3D Grid.

Noter que le maillage est dimensionnée en utilisatn les données du cadre de maillage. Si le cadre du
maillage n’existe pas, le maillage va couvrir le modèle conceptuel. De plus, le nombre de maille dans
les directions x et y ne peut pas être changé du fait que le nombre de lignes et de colonnes et la
localisation des mailles limites vont être contrôlé par les données d’affinage entrer dans les puits.

2. Dans Z-Dimension changer Num cells à 2.

3. Clique sur OK.

13 Initializing the MODFLOW Data

Maintenant le maillage est construit et les zones actives/inactives sont délimitées. L’étape suivante
consistera à convertir le modèle conceptuel à modèle numérique basé sur le maillage. Avant cela, on
doit tout d’abord initialiser les données de MODFLOW :

1. Clique sur 3D Grid Data dans Project Explorer.

2. Dans le menu principal exécuter la commande MODFLOW | New Simulation.
3. Cliquer OK.

14 Defining the Active/Inactive Zones

Maintenant le maillage est créer, l’étape suivante est de definer les zones actives et inactives du modèle.
Ceci est accompli automatiquement en utilisant les informations sur la couverture des termes
Puits/Sources local.
 1. Clique sur Map Data dans Project Explorer.
2. Cliquer sur l’outil Select Polygons .
3. Sélectionner un des polygones.
4. Cliquer sur le bouton Properties .
5. Confirmer que l’affectation des couches est 1 à 2 et cliquer sur OK.
6. Dans le menu principal executer la commande Feature Objects | Activate Cells in Coverage(s).

Chacune des mailles à l’intérieur de chaque polygone à l’intérieur de la couverture Puits/Sources est
désignée comme active et chaque maille à l’extérieur est désignée comme inactive. Noter que les mailles
sur la limite sont activées pour que la limite à flux nul au nord et la limite et la limite à potentiel impose
au Sud et à l’Est.

15 Interpolation des cotes des couches

Maintenant on a besoin de définir les côtes des couches et les charges initiales. Puisqu’on utilise le
package LPF, les cotes des toits et des murs sont définies pour chaque couche en considérant le type de
la couche. Pour un modèle à deux couches, on a besoin de définir les cotes du toit de la couche 1
(altitudes du sol), le mur de la couche 1 (substratum imperméable) et le mur de la couche 2. Il est supposé
que toit la couche 2 correspond au mur de la couche 1.

Pour définir les cotes, on importera des séries de points dispersés définissant les cotes et on interpolera
les cotes directement dans les mailles de chaque couche. Dans certains cas, cela est fait en utilisant une
série de points dispersés. Dans le présent cas, on va utiliser deux séries de points dispersés : la première
pour l’altitude de la surface du sol et la seconde pour les cotes du substratum de la couche 1 et le toit de
la couche 2.

15.1 Importing the Ground Surface Scatter Points

Les points disperses ont déjà été importés puisqu’ils ont été inclus dans le projet que nous avons lu au
départ. Ces points viennent d’un fichier texte. Les séries de points sont cachés qu’on peut faire apparaître
par :
1. Dans Project Explorer, cocher les cases terrain et elevs de 2D Scatter Data. Les points
apparaissent dans le modèle conceptuel

15.2 Interpolation des charges initiales et des cotes du terrain

Dans la suite on va interpoler les altitudes du sol et les charges initiales dans le maillage MODFLOW.
1. Cliquer par le bouton droit sur le série de point terrain et exécuter la commande
Interpolation To | MODFLOW Layers du menu flottant.

Cette boîte de dialogue informe GMS sur quelles données de MODFLOW les points seront interpolés.

2. Cliquer sur ground_elev et Starting Heads, et cliquer le bouton Map.

3. Cliquer sur ground_elev et Top Elevations Layer 1 puis cliquer le button Map.
4. Cliquer sur OK pour réaliser l’interpolation.

15.3 Interpolation des cotes des couches

Pour interpoler les cotes des couches:

1. Dans Project Explorer Cliquer sur elevs sous 2D scatter Data pour le rendre actif.
2. Cliquer par le bouton droit sur elevs et exécuter Interpolation To | MODFLOW Layers dans le
menu flottant.
GMS va automatiquement faire correspondre les valeurs Bottom 1 et Bottom 2 au mur de la couche 1 et
au mur de la couche 2 en se basant sur les noms des données.

3. Cliquer OK.

15.4 Ajustement de l’affichage

On peut maintenant cacher les series de points dispersés.

1. Décocher 2D scatter data dans Project Explorer.
2. Décocher Grid Frame dans Project Explorer.

15.5 Affichage d’une section transversal du modèle conceptuel

Pour vérifier l’interpolation on va afficher une section transversale.

1. Cliquer sur 3D Grid Data dans Project Explorer.

2. Clique sur une maille près du centre du maillage.
3. Cliquer sur Side View .

Pour avoir un meilleure affichage de la section, on va amplifier z.

To get a better view of the cross section, we will increase the z magnification.

4. Cliquer sur le bouton Display Options .

5. Entrer la valeur 5 pour Z magnification factor.
6. Cliquer OK.
7. Cliquer le bouton Frame .

Cliquer sur les flèches de Mini Grid pour se déplacer d’une colonne à l’autre et afficher toutes les
sections.

Noter que dans la partie droite de la section transversale, le mur de la couche 2 et au dessus de son toit.
Ceci doit être corrigé avant l’exécution du modèle.

15.6 Correction des cotes des couches

GMS fournit un certains nombre d’outils pratiques pour corriger les problèmes de cotes. Ces outils sont
inclus dans le vérificateur du modèle.
1. Dans le menu de GMS exécuter la commande MODFLOW | Check Simulation.
2. Cliquer sur le bouton Run Check.
3. Cliquer sur le bouton Fix Layer Errors.

Noter que plusieurs erreurs ont été trouvées. Il existe plusieurs méthode de corriger les erreurs des
couches. On choisira l’option de tronquer le substratum (Truncate to bedrock) option. Toutes les mailles
au dessous du substratum seront considérées inactives.

4. Cliquer sur l’option Truncate to bedrock.

5. Cliquer sur le bouton Fix Affected Layers.
6. Clique sur OK.
7. Cliquer sur le bouton Done.

Noter que les erreurs des couches on été corrigées. Les couches peuvent aussi être affichées en vue de
plan en cliquant sur le bouton Plan View .
16 Convertion du modèle Conceptual
1. Cliquer par le bouton droit sur le modèle conceptual ExempleMC et excuter la ommande Map
To | MODFLOW / MODPATH dans le menu folttant.
2. S’assurer que l’option All applicable coverages cocher puis cliquer sur OK.

Noter que les mailles interceptées par les drains, les puits et les limites à potential impose ont été
identifies par le modèle. Les charges hydrauliques et les cotes ont été déterminées par interpolation
linéaire. Les conductances des drains dans les mailles ont été déterminées en calculant la longueur du
drain qui intercepte la maille. De plus, les recharges et les conductivités hydrauliques ont été affectées
aux différentes mailles.

17 Vérification de la simulation

1. Cliquer sur 3D Grid Data dans Project Explorer.

2. Exécuter la commande MODFLOW | Check Simulation du menu principal.
3. Cliquer sur le buton Run Check.
4. Cliquer sur le bouton Done.

18 Enregistrement du projet

1. Cliquer sur le bouton Save .

19 Running MODFLOW
1. Exécuter la commande MODFLOW | Run MODFLOW du menu GMS.
2. Quand la solution est obtenue cliquer sur le bouton Close.

20 Affichage des lignes pizométriquesA set of contours should appear. To get

1. Exécuter la commande du menu Data | Contour Options.
2. Changer et teste différentes options.

21 Affichage du niveau de la nappe dans une vue latérale

1. Cliquer sur l’outil Select Cell .
2. Sélectionner une maille dans le modèle.
3. Cliquer sur le bouton Side View .
4. Cliquer sur le bouton Plan View .

22 Viewing the Flow Budget

1. Cliquer sur Map Data dans Project Explorer.

2. Cliquer sur l’outil Select Arcs .
3. Click on the rightmost drain arc.
1. Notice that the total flow through the arc is displayed in the strip at the bottom of the
2. window. Next, we will view the flow to the river.
3. Click on one of the specified head arcs at the bottom and view the flow.
4. Hold down the Shift key and select each of the specified head arcs.
5. Notice that the total flow is shown for all selected arcs. Flow for a set of selected cells
6. can be displayed as follows:
7. Select the 3D Grid Data Folder in the Project Explorer.
8. Select a group of cells by dragging a box around the cells.
9. Select the Data | Flow Budget command.
10. This dialog shows a comprehensive flow budget for the selected cells.
8. Select Done to exit the dialog.
9. Click anywhere outside the model to unselect the cells.