Modélisation Carte Piézo PDF

Modélisation en hydrogéologie: de la
carte piézométrique au modèle

numérique
Julio Gonçalvès
Maître de conférence UPMC
Extrait carte hydrogéologique de Grenoble Interprétation d’essai de puits (Trémosa et al., 2010)
Journées « modélisation en environnement », ECCOREV, 22-26/02/2010

Systèmes Carte Modèles
Piézométrique Conclusions
Hydrogéologiques numériques
Points abordés
•  Les systèmes hydrogéologiques modélisés

•  Un premier modèle « purement
mathématique » : la carte piézométrique
•  Modélisation numérique
•  Quelques réflexions et conclusions
Hydrodynamique souterraine
(Hyp.: ρ constant, milieu poreux)
h [m] : Charge hydraulique = potentiel hydraulique = cote piézométrique
h=P/ρg+z0 P : pression de l’eau (Pa)
Pa=0 ρ : masse volumique de l’eau kg/m3
(convention) g : accélération de la pesanteur (m2/s)
z0: cote de la mesure
Niveau piézométrique final h

Cote NGF de la surface libre
z
h constante sur une
colonne hydrostatique=altitude
de la surface libre
Tube crépiné
h ~ constante dans une nappe
Cote z0 de la mesure Entrée d’eau quelle que soit la cote de mesure
aquifère
aquitard
Hydrodynamique souterraine
Q : Débit (m3/s)
S : Section d’écoulement (m2)
Loi de DARCY (1803-1858) Δh : « Hauteur » de charge (m)
L : Longueur du milieu poreux (m)
Filtration d’eau à travers une carotte de milieu poreux…sous quel(s) effet(s)?

Loi de DARCY (1803-1858) Q : Débit (m3/s)

S : Section d’écoulement (m2)
Enseignement de l’expérience de Darcy? h : Hauteur de charge (m)
L : Longueur du milieu poreux (m)
K : Conductivité hydraulique
(m/s) (=perméabilité)
x = gradient hydraulique
L=Δx
Filtration d’eau à travers la carotte sous l’effet d’une différence de h = ’’Moteur de l’écoulement’’
Rq: que se passe-t-il si on coupe l’alimentation à gauche? ⇒U=Q/S
Notation différentielle (vecteur)
U=-K grad(h)
Signification de K et h…
Loi de DARCY (1802-1858) K : conductivité

hydraulique (m/s)
Roches meubles Sables fins et limons

Arènes granitiques Sables, graviers
Sables argileux
Argiles
10-10 10-9 10-8 10-7 10-6 10-5 10-4 10-3 10-2
Roches consolidées
Granites
Roches cristallines Craies
Schistes Calcaires, dolomies, grès, conglomérats
10-10 10-9 10-8 10-7 10-6 10-5 10-4 10-3 10-2

Très peu Peu perméable Perméable Très
perméable perméable
Conditions expérimentales : perméamètre horizontal

Débit Q « d’équilibre » Charge hydraulique h
hA
Mise en charge
hB h constante sur une
colonne hydrostatique=cote de la
L=10 cm Surface libre
A B
Conditions naturelles : écoulement dans un aquifère

z
h ~ constante dans une nappe
quelle que soit la cote de mesure
Écoulement
L=Δx=100 m
⇒ Écoulement depuis les charges élevées vers les charges faibles
• Aquifère : « corps (couche ou massif) de roches perméables comportant une zone saturée
et suffisamment conducteur d’eau pour permettre l’écoulement significatif d’une nappe
souterraine et le captage de quantités d’eau appréciables. »
• Nappe : « ensemble des eaux comprises dans la zone saturée d’un aquifère, dont toutes les
parties sont en liaison hydraulique »
• Perméabilité : capacité d’un milieu à conduire de l’eau
Aquifère dans sa partie Aquifère dans sa partie nappe captive

Nappe libre
F1
P F1
Lac SPnl SPnc
Toit de l’aquifère
aquifère
Mur de l’aquifère
Nappe libre
= contact direct avec l’atmosphère
Sol
Zone non saturée

Infiltration (matériau + eau + air)
Surface libre=surface pièzométrique
Zone saturée
Matériau = AQUIFERE
Eau = NAPPE SOUTERRAINE
Roche imperméable Substratum imperméable (mur)
La nappe libre « évacue » vers l’aval l’eau reçue par infiltration…

Exemples de nappes libres
- nappe de vallée en zone tempérée

Réservoir de sub-surface, l’infiltration par la pluie alimente la nappe ⇒ exutoires sont les
points bas de la topographie (source, rivières). Le drainage se fait par les vallées
suffisamment profonde (≠ vallées sèches).
- nappe alluviale
Nappe libre contenue dans les alluvions d’un fleuve
Niveau en équilibre avec celui du fleuve.
Nappe captive
= pas de contact direct avec l’atmosphère
Nappe captive Nappe libre

Infiltration possible
Artésianisme Infiltration nulle
Roche
imperméable Forage
Roche imperméable
Différences nappe libre / nappe captive
Nappe libre (eau non confinée)

•  Contact direct avec l’atmosphère
•  Existence d’une zone non saturée
•  Recharge (alimentation) rapide
•  Pas de protection naturelle face aux pollutions de surface
Nappe captive (eau confinée)
•  Pas de contact direct avec l’atmosphère
•  Niveau piézométrique> toit de l’aquifère
•  Recharge (alimentation) lente
•  Protection naturelle face aux pollutions de surface
•  Eau généralement fortement minéralisée
Réalisation d’une carte piézométrique
Réalisation d’un premier forage Réalisation d’un second forage
H = 94 m H = 94 m
H = 92 m
94 m
?
92 m
88 m Points
102,5 m 82 m d’observation :
94 m puits, forage,
piézomètre,
88 m source
99 m ? 90 m
Cote
92 m piézométrique
97 m 90 m 85 m interpolée
92,5 m
84 m
79 m
94 m
Géostatistique (statistique des variables régionalisées, Matheron 50-60)
Concept de base : la distribution d’une variable naturelle (%Fe, [O3], e…)

régie par une loi statistique ⇒ distribution est UNE réalisation d’un processus
aléatoire
milieu ⇔ infinité de variables aléatoires Z(X(x,y,z)) spatialement corrélées
Déterminer la structuration spatiale : corrélation entre Z(X) et Z(X+D)
45°: corrélation parfaite sinon écart

γ(D)

Si on fait varier D ⇒ fonction variogramme expérimental

γ(D) Puis calage d’une fonction théorique continue
Seuil C(0)=Var(Z(X))
Portée (distance au delà de laquelle plus de structure)
⇒ Combler les vides d’échantillonnage

par interpolation en tenant compte de la structure
•  Définition de la population à moyenner (voisinage)

•  Construction du variogramme expérimental à partir de tous les
couples de points
•  Choix d’un variogramme modèle ⇔ variogramme expérimental
•  Utilisation de ce modèle dans une démarche de comblement de
vide d’échantillonnage
Points
d’observation :
94 m puits, forage,
piézomètre,
source
?
Isopièze
92 m
Vocabulaire spécifique aux cartes piézométriques
Cote
piézométrique
Δx Ligne de
courant
Isopièze
(hydro-isohypse) Gradient hydraulique
i = ΔH / ΔX
Points
d’observation :
94 m puits, forage,
piézomètre,
source
? Sens
d’écoulement de l’eau
92 m
Diverses configurations d’écoulement
Zone divergente
Pompage
Zone convergente
Barrière
imperméable
Zone d’alimentation
60
40
35 55
30 50
25 45
Relation Nappe/Rivière : nappe drainée par la rivière
Rivière
puits
plateau
infiltration
ou forage
125
de la pluie
124
vallée
rivière
123
surface
de la nappe
nappe écoulement (niveau de l’eau
dans les puits)
roche poreuse
couche imperméable formant la base de la nappe

Relation Nappe/Rivière : nappe alimentée par la rivière
puits
ou forage
Rivière
ruissellement 125
oued
en crue 124
infiltration 123
surface
de la nappe
nappe (niveau de l’eau écoulement
dans les puits) roche poreuse
base de la nappe
Équation de conservation (continuité) + loi d’écoulement (Darcy)
(ω: porosité)
Δh?
Sol
Infiltration Roche imperméable
Roche imperméable
Roche imperméable
Un modèle numérique pour quoi faire?
•  Contexte trop complexe (milieu hétérogène, non linéarité

des processus…) pour avoir des solutions simples
•  Tester des hypothèses de travail : fonctionnement, la
réactivité, l’importance relative de processus
•  Une fois validé, il peut constituer un outil prédictif: outils
pour la gestion des ressources et/ou de la qualité des
eaux
Étapes de la construction du modèle
•  Description géométrique du domaine (milieu géologique)

•  Construction du maillage (méthodes numériques)
•  Estimation de l’hétérogénéité des propriétés : habillage
du maillage
•  Conditions initiales et limites
•  Calage et Validation
Forages ⇔ bases de données

Exemple : Bassin de Paris
Log de Forage
Géométrie 2D, 3D
Très onéreux et rare…Géophysique


du maillage
Différentes méthodes numériques « classiques » : Différences finies, Volume

finis et Éléments finis
Inconnue: h(x,y)
Pas de solution analytique (K hétérogène) ⇒ solution approchée
i,j+1
Différences finies y
Application deux fois de l’approximation de Taylor :
x
Approximation de Taylor i-1,j i,j i+1,j
i,j-1
e(x,y)
Perméabilité à l’interface ≈
perméabilité équivalente
entre i et i+1 ou i-1⇒ pb de stabilité numérique
Quand milieu trop contrasté
Volumes finis, éléments finis…


du maillage

du maillage
Conditions de Dirichlet : limites à charge h imposée

• Contact nappe avec un plan d’eau libre (mer, lac, rivière, ….)
⇒La charge hydraulique est constante, imposée

par la côte de la rivière (f(x,t)
• Ligne de source
⇒ La charge hydraulique est constante, imposée par la côte de la source
• Cas de la nappe captive ayant un taux d’infiltration > au flux d’eau s’écoulant dans la nappe
⇒ h imposé = cote de l’affleurement

Conditions de Neumann : limites à flux imposé
• Contact nappe avec un imperméable (ex. faille)
60
55
50
⇒ Limite à flux nul
45
• Cas de la nappe captive ayant un taux d’infiltration < au flux d’eau s’écoulant dans la nappe
⇒ Limite à flux imposé qui est le taux d’infiltration de la pluie


du maillage
Thématiques de modélisation hydrogéologiques

actuelles
•  Ressources en eaux (stress hydrique…)

•  Sureté et impacts stockages (déchets nucléaires, gaz…)
•  Processus de couplages (écoulements « non purement
darciens » dans les argiles
•  Écoulements multiphasiques (Pétrole, Gaz, zone non-
saturé, pollutions)
•  Transport réactif (précipitations, dissolution)
•  Contexte côtier mélanges eaux douces/eaux salées
•  « Accommodation climatique » des aquifères/rivières…
Exple 1: Modélisation de l’écoulement dans

l’Albien, Raoult, 1999
Influence de l’exploitation sur la ressource et la

chimie de l’aquifère de l’Albien
E.1 Définition de la (des) géométrie(s)

E.2 Synthèse et acquisition

des données
Carte Piézométrique 1930

≈ R. Permanent
Carte Piézométrique 1995

Exploitant, BRGM,
Raoult(1995)
E.3 Modélisation « 3D»
Maillage et CL
Calage/Validation
Simulation ‘à l’actuel’
Scénario d’exploitation

actuelles

darciens » dans les argiles
STRATIGRAPHY LITHOLOGY
MAIN RESERVOIR BEDS
Albian Thickness (m)

Neocomian 10 -‐ 50
50 -‐ 100
100 -‐ 150
Lusitanian 150 -‐ 200
CALLOVO-‐OXFORDIAN 200 -‐ 250
250 -‐ 300
Dogger target beds
300 -‐ 350
TOARCIAN
Rhaetian
Keuper target beds
Ss (m-‐1)
(0.56 -‐ 1.13).10-‐4
(1.14 -‐ 1.23).10-‐4
(1.24 -‐ 1.80).10-‐4
(1.81 -‐ 5.14).10-‐4
(0.52 -‐ 2.46).10-‐3
(0.26 -‐ 1.38).10-‐2
Exple 2 : Impact
hydrogéologique
du stockage -Log[Kh(m/s)]
Anne JOST, Julio GONÇALVES,
Guillaume BOUQUARD and Hugo DAYAN
de CO2? <4
Deep Saline Aquifers for Geological Storage
4-5
5-6 of CO2 and Energy
6-7 27 - 29 May 2009, IFP, Rueil-Malmaison (France)
7-8
•  Modèle 4 couches extraites

d’un 3D complet maillage 5
km*5km avec ModFlow96
(PMWIN package)
•  Lias et Kim. Imperméables et
charges imposées à
l’affleurement
•  Simulation en permanent sans
injection
•  Simulation transitoire: Injection
« équivalent eau douce » à 5
Mt/an CO2⇔ 9Mm3/an sur un
puits pendant 50 ans.
Fin injection Fin injection
Δh (m)
< 0.2
0.2 -‐ 0.5
0.5 -‐ 1
1 -‐ 2
2 -‐ 3
3 -‐ 5
5 -‐ 10
10 -‐ 25
25 -‐ 50
50 -‐ 100
100 -‐ 150
+500 ans +100 ans

Impact sur la drainance flux

d’échange verticaux entre
deux couches en % par
rapport au régime
permanent (sans injection)

actuelles

darciens » dans les argiles)
•  Adapter le modèle et la méthode au problème posé

•  Tester des hypothèses et/ou prédiction
•  Implique test de sensibilité ⇒ erreur de la
modélisation
•  Réflexion avant de concevoir le modèle
« Ce n’est donc pas au modèle que l’on devrait faire
confiance mais à l’homme derrière le modèle, à son
expérience et à sa capacité à représenter par le
modèle ce qu’il a perçu de la réalité représentée » G.
de Marsily

Modélisation Carte Piézo PDF

Transféré par

Droits d'auteur :

Formats disponibles

Modélisation Carte Piézo PDF

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Modélisation Carte Piézo PDF

Transféré par

Droits d'auteur :

Formats disponibles

Modélisation en hydrogéologie: de la

carte piézométrique au modèle

Journées « modélisation en environnement », ECCOREV, 22-26/02/2010

• Les systèmes hydrogéologiques modélisés

Niveau piézométrique final h

Filtration d’eau à travers une carotte de milieu poreux…sous quel(s) effet(s)?

Loi de DARCY (1803-1858) Q : Débit (m3/s)

Loi de DARCY (1802-1858) K : conductivité

Roches meubles Sables fins et limons

10-10 10-9 10-8 10-7 10-6 10-5 10-4 10-3 10-2

Schistes Calcaires, dolomies, grès, conglomérats

10-10 10-9 10-8 10-7 10-6 10-5 10-4 10-3 10-2

Conditions expérimentales : perméamètre horizontal

Conditions naturelles : écoulement dans un aquifère

• Perméabilité : capacité d’un milieu à conduire de l’eau

Aquifère dans sa partie Aquifère dans sa partie nappe captive

Zone non saturée

Surface libre=surface pièzométrique

Roche imperméable Substratum imperméable (mur)

La nappe libre « évacue » vers l’aval l’eau reçue par infiltration…

- nappe de vallée en zone tempérée

Nappe captive Nappe libre

Différences nappe libre / nappe captive

Nappe libre (eau non confinée)

Réalisation d’une carte piézométrique

Réalisation d’un premier forage Réalisation d’un second forage

Réalisation d’une carte piézométrique

Réalisation d’une carte piézométrique

Géostatistique (statistique des variables régionalisées, Matheron 50-60)

Concept de base : la distribution d’une variable naturelle (%Fe, [O3], e…)

milieu ⇔ infinité de variables aléatoires Z(X(x,y,z)) spatialement corrélées

Déterminer la structuration spatiale : corrélation entre Z(X) et Z(X+D)

45°: corrélation parfaite sinon écart

Si on fait varier D ⇒ fonction variogramme expérimental

Portée (distance au delà de laquelle plus de structure)

⇒ Combler les vides d’échantillonnage

• Définition de la population à moyenner (voisinage)

Réalisation d’une carte piézométrique

Vocabulaire spécifique aux cartes piézométriques

Réalisation d’une carte piézométrique

Diverses configurations d’écoulement

Diverses configurations d’écoulement

Relation Nappe/Rivière : nappe drainée par la rivière

couche imperméable formant la base de la nappe

Diverses configurations d’écoulement

Relation Nappe/Rivière : nappe alimentée par la rivière

Équation de conservation (continuité) + loi d’écoulement (Darcy)

Infiltration Roche imperméable

Un modèle numérique pour quoi faire?

• Contexte trop complexe (milieu hétérogène, non linéarité

Étapes de la construction du modèle

• Description géométrique du domaine (milieu géologique)

Forages ⇔ bases de données

Très onéreux et rare…Géophysique

Étapes de la construction du modèle

• Description géométrique du domaine (milieu géologique)

Différentes méthodes numériques « classiques » : Différences finies, Volume

Volumes finis, éléments finis…

Étapes de la construction du modèle

•  Les systèmes hydrogéologiques modélisés

⇒ Combler les vides d’échantillonnage

•  Définition de la population à moyenner (voisinage)

•  Contexte trop complexe (milieu hétérogène, non linéarité

•  Description géométrique du domaine (milieu géologique)

•  Description géométrique du domaine (milieu géologique)

•  Description géométrique du domaine (milieu géologique)

•  Description géométrique du domaine (milieu géologique)

•  Description géométrique du domaine (milieu géologique)

•  Ressources en eaux (stress hydrique…)

•  Ressources en eaux (stress hydrique…)

•  Modèle 4 couches extraites

•  Ressources en eaux (stress hydrique…)

•  Adapter le modèle et la méthode au problème posé