Création d’un maillage

4 minutes de lecture

Projet maillage

La première étape consiste à créer un projet maillage, lequel sera utilisé pour le calcul d’écoulement et transport 2D saturé.

Si THYRSIS est activé, un menu Thyrsis doit apparaître dans la barre de menu principale de QGIS. Pour créer un projet maillage, cliquer sur : Thyrsis > Nouveau > Maillage... Figure 2 :

../../_images/tutorial_newmesh.png

Figure 90 Création d’un nouveau maillage

Note

THYRSIS crée le projet maillage à partir de la projection courante. Changer celle-ci dans QGIS pour être certain d’utiliser la bonne projection. Dans ce tutoriel, toutes les données sont projetées en EPSG: 27572 (Lambert II Etendu) :::

Le chemin du répertoire stockant la base de données maillage doit être renseigné. Le répertoire ~/.thyrsis/mesh est celui dans lequel devra être placée la base de données créée pour pouvoir être utilisée dans THYRSIS ultérieurement, mais il est conseillé de choisir un autre emplacement pour la mise au point du modèle. La base de données doit être nommée SITE_MODELE.mesh.sqlite, où :

  • SITE désigne le nom du site - créé et sauvegardé dans la base sites.sqlite s’il n’existe pas déjà,

  • MODEL désigne le nom du modèle - un modèle étant la combinaison d’un maillage et de champs ou paramètres hydrogéologiques.

Ici la base est dénommée SACLAY_TUTORIAL, pour indiquer que le maillage est basé sur le site de SACLAY et qu’il s’agit du modèle TUTORIEL. THYRSIS crée également le projet QGIS (fichier .qgs). A l’issue de cette étape le répertoire de maillage doit ressembler à la Figure 3.

../../_images/tutorial_meshdirectory.png

Figure 91 Répertoire maillage

Avoir une couche de fond peut être utile pour localiser l’environnement voisin du project. Le fond de carte OpenStreetMap peut être ajouté grâce au menu d’exploration des couches de QGIS Figure 4.

../../_images/tutorial_osmlayer.png

Figure 92 Canevas QGIS avec couche OSM

Il est possible d’ajouter une couche de rivière, qui sera utile pour délimiter la zone du projet. Pour le tutoriel, un fichier river.gpkg est fourni. Cliquer sur Couche>Ajouter une couche>Ajouter une couche vecteur et ajouter la couche de rivière Figure 5.

../../_images/tutorial_addVectorLayer.png

Figure 93 Ajouter une couche vecteur

Pour finir, deux rasters sont nécessaires pour la création du maillage :

  • un modèle numérique de terrain (MNT), nommé mnt.tif

  • un raster fournissant l’altitude du mur de l’aquifère, dénommé mur.tif

Pour les ajouter, cliquer sur Couche> Ajouter une couche> Ajouter une couche raster.

Il est recommandé de placer le fond cartographique et les couches raster sous les autres couches et d’ajouter de la transparence dans les propriétés des couches raster. Le canevas de la carte doit ressembler à :

../../_images/tutorial_mntlayer.png

Figure 94 Canevas QGIS avec couche MNT

Ne pas oublier de sauvegarder le projet (CTRL + S ou cliquer sur Projet > Sauvegarder)


Maillage

L’étape suivante est la création du maillage 2D.

La couche contours doit être définie en premier. C’est une couche de polylignes avec les attributs suivants :

  • groupe : nom du groupe auquel le contour appartient,

  • nom : nom du contour,

  • longueur_element : longueur caractéristique qui définit la densité du maillage dans le voisinage du contour,

  • potentiel_impose : champ égal à 1 si le contour considéré est à potentiel imposé, égal à NULL ou 0 sinon.

Warning

Pour créer les contours, activer l’accrochage dans QGIS et sélectionner l’option couche active seulement. Ceci est nécessaire afin de créer des contours formant des ensembles clos.

../../_images/tutorial_snapping.png

Figure 95 Activer l’accrochage

Ensuite sélectionner la couche de contours, cliquer sur le bouton d’édition,

../../_images/tutorial_editing.png

Figure 96 Édition de la couche contours

et activer l’outil Ajouter une entité linéaire (Figure 9).

../../_images/tutorial_addLineFeature.png

Figure 97 Ajouter une polyligne dans la couche contours

Pour le tutoriel, nous voulons créer un maillage sur le Plateau de Saclay. Nous recommandons de créer 4 contours :

  • un contour au nord le long de la rivière,

  • un contour à l’est entre les deux rivières,

  • un contour au sud le long de la rivière,

  • un contour à l’ouest perpendiculaire aux rivières nord et sud.

ous les contours doivent se trouver dans l’étendue des couches du MNT et du mur. À chaque validation d’un contour, la fenêtre de latutorial_contoursAttributes s’ouvre.

../../_images/tutorial_contoursAttributes.png

Figure 98 Attributs de la couche contours

Le champ longueur_element doit être renseigné et le champ potentiel_impose doit être égal à 1 si le contour correspond à une condition limite à potentiel imposé.

Des contours imbriqués peuvent être créés de manière à définir des zones avec différentes densités de maillage.

QGIS ne permet pas l’accrochage d’une entité sur elle même lors de la création de l’entité. Pour définir un contour fermé d’un seul tenant, créer un contour presque fermé, renseigner la longueur élémentaire puis utiliser l’outil de noeud pour déplacer le dernier point sur le premier (avec accrochage).

Enregistrer la couche contours. THYRSIS remplit alors la couche domaines avec les polygones définis par la couche contours.

L’absence de polygone créé dans la couche domaines est généralement lié à des erreurs topologiques, aucun contour fermé n’ayant été trouvé.

Le canevas devrait ressembler à l’image de la Figure 11, avec des attributs de contours similaires.

../../_images/tutorial_contoursCanvas.png

Figure 99 Canevas QGIS après validation de la couche contours

Pour créer le maillage, cliquer sur le bouton Maillage (Figure 12) dans la barre d’outils de THYRSIS.

../../_images/tutorial_meshButton.png

Figure 100 Bouton maillage

La fenêtre de la Figure 13 s’ouvre.

../../_images/tutorial_meshDialog.png

Figure 101 Dialogue de maillage

Verifier les couches et lancer la création du maillage. Les couches noeuds and mailles sont actualisées et les couches altitude et altitude_mur sont ajoutées.

Le canevas QGIS de la Figure 14 doit être obtenu.

../../_images/tutorial_meshResult.png

Figure 102 Canevas QGIS avec le résultat du maillage