Cours COMPLET Hydrogeologie BTP

Filire dIngnieurs:
Btiment & Travaux Publics
Mme Malki Mouna

malki.mouna@gmail.com
Anne universitaire : 2016/2017

Hydrogologie M.Malki
I. Gnralit, dfinitions et notions de base
Aquifre:
Un aquifre est une formation hydrogologique permable permettant lcoulement

significatif dune nappe deau souterraine et le captage de quantits deau apprciables, par
des moyens conomiques. Cest une couche de terrain ou une roche, suffisamment poreuse
(qui peut stocker de leau) et permable (o leau circule librement), pour contenir un
rservoir deau souterraine.
Configuration de laquifre
La configuration de laquifre porte sur ses dimensions et les caractristiques de ses limites
gologiques et hydrodynamiques ou conditions aux limites.
La base de laquifre, appele substratum, est constitu par une formation hydrologique
impermable. Par contre sa limite suprieure est de trois types :
Hydrodynamique avec fluctuations libre : aquifre nappe libre
Gologique impermable : aquifre nappe captive
Gologique semi-permable : aquifre nappe semi-captive
Type daquifre :
Aquifre nappe libre
Formation gologique permable partiellement sature, limite par un plancher

impermable (aquitard infrieur) et o il existe une surface sur laquelle la pression de leau
est gale la pression atmosphrique. La surface pizomtrique constitue la limite suprieure
de laquifre. Cest une limite hydrodynamique. Cette limite peut sabaisser librement dans la
formation hydrogologique permable, do la dnomination aquifre nappe libre.
Aquifre nappe captive
Dans les aquifres plus profonds les eaux souterraines sont emprisonnes dans la formation
hydrogologique permable, entre deux formations impermables fixes : le substratum la
base et le toit au sommet.
Aquifre nappe semi-captive
Le toit ou le substratum (ou les deux) de laquifre sont souvent constitus par une formation
hydrogologique semi-permable. Celle-ci permet, dans certaines conditions
hydrodynamiques favorables (diffrence de charge) des changes deau (ou de pression) avec
laquifre superpos ou sous-jacent, appel drainance. Ce phnomne implique un aquifre
nappe semi-captive.
Anne universitaire 2016 / 2017 2

Aquiclude : Formation gologique qui ne laisse pas passer leau. Ex.: argiles
Aquitard : Formation qui ne peut pas alimenter un puits mais qui peut alimenter un aquifre
(Ex: des sables fins, le silt, les tills remanis
Structure de laquifre
Aquifre, complexe unique rservoir / eau souterraine
Complexe physico-chimique unique de deux constituants essentiels, ou phase, troitement lis

et en interactions :
Le rservoir: phase solide, milieu poreux ou fissur, constitue la trame de la structure,

squelette solide ou matrice (ex: grains de sable dune formation sableuse, roche fissure
de la craie, etc.)
Leau souterraine: phase liquide, dont la fraction mobilisable (eau gravitaire) constitue la
nappe souterraine alimentant les sources, rivires et captages.
Aquifre multicouche
Une combinaison de formations hydrogologiques semi-permables intercales entre des

formations permables, identifie un aquifre multicouche. Cest un systme hydrologique, car
chaque aquifre nappe semi-captive ne peut tre considr de manire indpendante.
Fonction du rservoir
Les rservoirs, suivant leurs configurations et leurs structures, peuvent remplir une ou
plusieurs de ces 3 grandes fonctions :
1 - Fonction capacitive : emmagasinement de leau, il stocke et libre de leau)
2 - Fonction conductrice : il permet la circulation de leau

3 - Fonction dchange : interactions physico-chimiques permanentes entre leau souterraine

et le rservoir (interaction eau / roche) et la roche (Ex: Chaleur, dissolution ou prcipitation
des sels, change dions, etc)
Zonalit :
Zone non sature :
Zone o la totalit de la porosit nest pas remplie deau. La circulation de leau se fait
principalement verticalement.
Zone sature :
Zone o la totalit de la porosit est remplie deau. La circulation de leau se fait

principalement horizontalement dans la direction o le gradient hydraulique est le plus fort.
Types deau dans un aquifre :
Eau de rtention (eau lie)
Englobe toute leau retenue dans les vides dun milieu poreux, satur ou non la surface des
grains ou des parois des microfissures. Elle est maintenue la surface du solide par de trs
grandes forces plus grandes que la gravit. Elle n'est donc pas mobilisable.
Eau adsorbe = hygroscopique : Lie la roche par attraction molculaire (adsorption)
Eau pelliculaire = eau dadhsion : Entoure les particules du sol et leau adsorbe dune
mince pellicule.
Eau capillaire
Eau maintenue dans un milieu poreux par des forces de capillarit: tension superficielle.
Eau gravitaire
Eau qui subit laction de la gravit. C'est l'eau mobilisable. Elle seule circule dans les
aquifres et peut tre exploite.
II. Identification hydrogologique du complexe eau/rservoir
Porosit
La plupart des roches contiennent des vides entres les grains qui les constituent, ces vides
peuvent tre occups par des fluides (eau et huile) ou des gaz. On parle de milieu poreux (Ex:
Sable, grs). Les fonctions conduite et rservoir, sont dtermines essentiellement par les
dimensions et les interconnections des vides. Ces dernires assurent la continuit du milieu
aquifre. Ltude morphologique des vides porte sur leur nature, leur forme et leurs

dimensions. Deux grands types de vides, pores et fissures, caractrisent respectivement le

milieu poreux et le milieu fissur.
Morphologie des pores et milieu poreux
Les pores sont des vides de forme plus ou moins sphrique, de petite dimensions (ordre de
grandeur millimtrique), mnags entre les particules solides ou grains, constituant le
rservoir.
Les dimensions des vides sont troitement lies celles des grains, dont la mesure est plus
directement accessible. Les diamtres des grains des roches meubles permables stalent
dans une gamme de 0,06 0,001 mm, soit dordre de grandeur micromtrique, dont les
argiles, milieu dit impermable.
Interconnexion des pores et milieu continu
Les pores communiquent entre eux, dans le sens de lcoulement de leau souterraine,
permettant le dplacement des particules deau (eau gravitaire). Celles-ci suivent des trajets
ou trajectoires, plus ou moins compliqus, identifiant les lignes de courant. Cet agencement
caractrise la continuit du milieu poreux qui est une des conditions de base pour la validit
des lois de lhydrodynamique souterraine.
Morphologie des fissures et milieu fissur
Les fissures sont des fentes de forme allonge, ouverture plus ou moins large. Leur
ensemble constitue la fissuration, phnomne naturel dont lorigine est essentiellement
mcanique. Les fissures sont classes suivant leur dimensions, en deux types : les
microfissures et les macro-fissures
Microfissures
Caractriss par une ouverture de quelques diximes de millimtres et une longueur de lordre
mtrique dcamtrique. Leur rle est comparable celui des pores interconnectes du milieu
poreux continu.
Macro-fissures
Ouverture suprieure quelques millimtres . Ce sont les zone de broyage et des failles. Le
milieu fissur est discontinu.

Dans les aquifres poreux, leau est contenue dans les pores ouverts de la roche et peut y
circuler librement (sables, craie, graviers, grs, scories volcaniques, etc.). La permabilit est
matricielle.
Dans les aquifres fissurs, leau est contenue et circule dans les failles, fissures ou diaclases
de la roche (calcaires , granites, coules volcaniques, etc.). La permabilit est fissurale.
Les aquifres karstiques sont des systmes complexes particuliers associant une zone
superficielle plus ou moins fissure et insature (en eau) servant de zone dinfiltration, et une
zone infrieure fissure, prsentant galement des conduits, grottes etc. Cette zone est sature
en dessous dun certain niveau et leau circule avec de grandes vitesses comparativement aux
systmes poreux.
Porosit totale (n)
La porosit totale (n) est la proprit d'un milieu poreux, ou fissur, de comporter des vides
interconnects ou non. C'est le rapport du volume des vides au volume total de la roche. Il est
exprim en %. Son inverse est la compacit.
n = (Volume des vides / Volume total )*100 = (Vv / Vt )*100
Porosit efficace ou effective
Ou encore porosit de drainage correspond leau mobilisable ou exploitable libre sous

laction de la gravit. C'est le rapport du volume de leau gravitaire au volume total de la
roche. Il est exprim en %.
ne = Volume de leau gravitaire / Volume total )*100 = (Vg / Vt )*100
Plus la taille des grains est rduite, plus la porosit efficace est rduite.
Indice des vides (e)
Il est dfini comme le rapport du Volume des vides au volume du solide. Il est exprim en %.
e = volume des vides / volume du solide = Ve / Vs
Capacit de Rtention ou coefficient de rtention (Cr)
Cest la quantit d'eau lie aux particules et/ou capillaire. La porosit n peut se dcomposer
en porosit efficace et porosit de rtention:
n = ne + Cr = Vr / Vt * 100 et Cr = (Vr-Vg) / Vt * 100

Degr de saturation (Sr) :
Il caractrise le pourcentage d'eau contenu dans les vides :
Sr = Volume deau contenu dans les vides / Volume totale des vides de la roche
Permabilit
Aptitude d'une roche conduire l'coulement de l'eau. Ce paramtre permet un classement des
roches en trois grandes catgories:- Formations permables ; impermable et semi-
permables .
Remarque : Ne pas confondre la porosit qui est la proprit du rservoir de stocker ou de

librer de l'eau souterraine et la permabilit qui est son aptitude conduire son coulement.
III. Ecoulements souterrains : Loi de Darcy
Loi de Darcy : applications
Enonc de la loi de Darcy
Darcy a montr que le volume deau, Q en m3/s, filtrant de haut en bas dans la colonne de
sable de hauteur l en m, travers la section totale, perpendiculaire la direction verticale
dcoulement, A en m, est fonction dun coefficient de proportionnalit, K en m/s,
caractristique du sable et de la perte de charge par unit de longueur du cylindre de sable h/l
sans dimension. Do lexpression de la loi de Darcy :
Q (m / s) = k(m/s) * A (m) * h/l
Le terme K est dfini comme un coefficient, dpendant de la permabilit de la couche ,

appel coefficient de permabilit.

Le quotient de charge h, par la longueur de la colonne de sable l, ou perte de charge par unit
de longueur h/l , est dfini comme le gradient hydraulique.
Do avec i= h/l, lexpression devient :
Q= K*A*i
Le dbit unitaire q est le dbit en m3/s traversant lunit de section, perpendiculaire la

direction dcoulement en milieu satur, dans lunit de temps en secondes. Cest aussi la
quantit deau traversant le milieu satur par unit de surface. tant le quotient dun dbit par
une surface, il a la dimension dune vitesse et sexprime en m/s.
Do, en combinant les deux dernires expressions, on obtient:
q (m/s) = K (m/s) * i
La vitesse de filtration V en m/s, rapporte la section totale A est :
V (m/s) = Q/A = K*i= q
La vitesse de filtration est ainsi la vitesse fictive dun flux deau en coulement uniforme,
travers un milieu aquifre satur, dduite du dbit dcoulement travers par ce flux. Elle est
gale au dbit unitaire.
Dispositif de laboratoire avec coulement latral
Un dispositif de laboratoire, avec coulement latral, reprsente mieux la circulation des

eaux souterraines dans laquifre. Les rsultats, obtenus par cette exprience, permettent
dcrire :
Q = - K *((h1 h2) / l)
(h1-h2)/l est le gradient hydraulique i ;
h1-h2 est la diffrence de charge h
Le signe ngatif est introduit devant le membre de droite de lquation car la charge dcroit
dans le sens de lcoulement et q ou V ne peuvent tre ngatif.
Mesure du niveau pizomtrique :
Il est calcul par diffrence entre la cote du sol, z et la profondeur de leau, Hp (H=z-Hp).
La profondeur de leau est mesure par les sondes: ficelles ou ruban avec flotteur, sondes
lectriques. Leur prcision est de lordre de plus ou moins 5 mm. Souvent les niveaux sont
enregistrs automatiquement, en continu, par les limnigraphes qui donnent des limnigrammes
pizomtriques.

Laltitude du sol est obtenue, soit par lecture de la carte topographique, soit lorsquune plus
grande prcision est recherche, par oprations de nivellement.
Gradient hydraulique :
Dans la pratique, le gradient hydraulique est calcul sur le terrain, laide des niveaux
pizomtriques mesurs dans deux ouvrages dobservation, aligns sur une ligne de courant.
Lun amont, H1, lautre aval, H2, spars par une distance L : i = (h1 h2) / l.
Mais la mthode recommande est celle de lutilisation des cartes pizomtriques. Les
valeurs des gradients hydrauliques, mesurs dans les conditions naturelles, sont faibles, de
0,001 0,00001
Conditions de validit de la loi de Darcy
La loi de Darcy est tablie par des expriences de laboratoire rpondant des conditions trs
strictes.
Quatre conditions doivent tre respectes: Continuit, Isotropie, Homognit du

rservoir et coulement laminaire.
Permabilit transmissivit diffusivit
Coefficient de permabilit
La permabilit est laptitude dun rservoir se laisser traverser par leau, sous leffet dun
gradient hydraulique. Elle exprime la rsistance du milieu lcoulement de leau qui le
traverse.
Le coefficient de permabilit, not K, est dfini par la loi de Darcy. Cest le volume deau
gravitaire en m3 traversant une unit de temps (seconde), sous leffet dune unit de gradient
hydraulique, une unit de section en m orthogonale la direction de lcoulement, dans les
conditions de validit de la loi de Darcy ( la temprature 20C). Il a la dimension dune
vitesse et sexprime en m/s.
La permabilit intrinsque, note k, est le volume de liquide en m3 dunit de viscosit

cinmatique traversant en une unit de temps (s), sous leffet dune unit de gradient
hydraulique, une unit de section (m) ou en darcy. Elle est parfois appele permabilit
gomtrique.
Facteurs du coefficient de permabilit
Daprs les deux expressions gnralises de la loi de Darcy, nous obtenons :
K = N d * /
Avec :

d10 , le diamtre efficace des grains, en cm ;
N, un facteur de forme sans dimensions, groupant les autres caractristiques

granulomtriques: forme et arrangement des grains.
, la viscosit dynamique qui exprime la rsistance du liquide l'coulement;
= g , le poids volumique, exprimant la force motrice (action de la force de gravit)
Dbit dune nappe
Le dbit dune nappe, Q , est le volume deau en m3 traversant par unit de temps(s) une
section transversale en m daquifre, sous leffet dun gradient hydraulique dtermin. Pour
simplifier les tudes sur le terrain, cest le plan vertical perpendiculaire la direction
dcoulement. Il est calcul par application des expressions de la loi de Darcy.
Ce dbit ne doit pas tre confondu avec le dbit de lcoulement souterrain Qw.
Qw reprsente les sorties du bassin hydrogologique ou de laquifre. Cest--dire son

drainage par les cours deau et lalimentation des sources du bassin hydrologique.
Calcul du dbit dune nappe
Il existe diffrentes mthodes pour calculer le dbit dune nappe.
Mthode de la carte pizomtrique
Le dbit peut tre calcul par interprtation de la carte pizomtrique.
Mthode de la section totale
Q=K*A*i
Mthode des sections lmentaires
La mthode prcdente est peu prcise. Cest pourquoi il est prfrable, lorsque les donnes
sont suffisantes, de procder au calcul des sections lmentaires. La section gnrale est
subdivise en section, de nombre gal celui des sondages dessai.

Lexpression de Darcy, applique par exemple la section 9, est :
Q = Km * bm * L * i
Km est la moyenne arithmtique des coefficients de permabilit mesurs dans les sondages
encadrant la section.
Exemple, pour la section 9:
Km = (K8 + K9)/2
bm est lpaisseur moyenne de laquifre dans les sondages 8 et 9.
bm = (b8 + b9)/2
i est le gradient hydraulique calcul sur une ligne de courant, trace sur la carte
pizomtrique.
Un calcul identique est effectu pour chaque section lmentaire. Le dbit de la nappe est
gal la somme des dbits traversant chaque section.
Le calcul du dbit peut galement tre effectu avec les transmissivit moyennes, en
appliquant lexpression :
Q = T*L*i
Vitesse de filtration et vitesse effective: Hydrodynamique souterraine
Lhydrodynamique souterraine, dont la base est la loi de Darcy, considre que lcoulement
travers un milieu, homogne et continu, seffectue selon des trajectoires thoriques

rectilignes, indpendantes de la structure microscopique du rservoir. Le trajet de la droite

moyenne ( ligne de courant) passe, indiffremment; travers les grains ou les pores. Cest
pourquoi la loi de Darcy nest valable que pour certaine grandeur de milieu, comprenant un
nombre suffisant de pores, donc dchelle macroscopique.
La vitesse de filtration V, calcule par la loi de Darcy, se rapporte la section totale A. Elle
na pas de ralit physique : Vitesse de filtration = Q / A
Vitesse efficace
La surface efficace dcoulement, ainsi rduite aux vides mnags par le corps solide ( grain +
eau de rtention), dpend de la porosit efficace, ne. Elle est gale A.ne. Lexpression de la
loi de Darcy, corrige, rapporte la section efficace, pour le calcul de la vitesse effective, Ve,
est donc :
Diffusivit
La diffusivit, note T/S, rgit la propagation dinfluences dans laquifre. Elle caractrise la
vitesse de raction d'un aquifre lors d'une perturbation (variation de niveau de la rivire, de la
nappe, pompage). Elle est gale au rapport entre la transmissivit et le coefficient
demmagasinement. Elle sexprime en m/s.
Transmissivit
La productivit dun captage dans un aquifre est fonction de son coefficient de permabilit,
K et de son paisseur, b. cest pourquoi un paramtre rcent, la transmissivit, note T, a t
cr. Il rgit le dbit deau qui scoule, par unit de largeur, L, dun aquifre, sous leffet
dune unit de gradient hydraulique, i . Il value la fonction conduite de laquifre.
La transmissivit est gale au produit du coefficient de permabilit K, par lpaisseur de

laquifre, b. elle sexprime en m/s.
T= K * b
Lexpression de la loi de Darcy , Q= K*A*i , devient avec A=b*L :
Q=T*L*i
Incluant lpaisseur de laquifre, la transmissivit permet de reprsenter sur des cartes, les
zones de productivit. Elle est la base de la discrtisation du calcul par mailles des modles
mathmatiques. Elle est mesure, sur le terrain, par les pompages dessai.
Coefficient demmagasinement- Dbit dune nappe- vitesse dcoulement conditions

aux limites
Coefficient demmagasinement

Le coefficient demmagasinement S est le rapport du volume deau libr par unit de

surface, sous une charge hydraulique unitaire, et dtermine la fonction capacitive du
rservoir. Il exprime, en nappe libre, la capacit de vidange des pores de la roche et peut
donc tre assimil sa porosit efficace. En nappe captive, il rsulte de la dcompression de
la roche et de leau contenue dans ses pores. Sa valeur varie de 0,2 0,01 pour les nappes
libres et de 0,001 0,0001 pour les nappes captives.
Vitesse de dplacement. Hydrocinmatique souterraine. Dispersion
Lhydrocinematique, branche de la cinmatique, considre les dplacements rels des

particules deau dans les vides continus. Elle tudie les trajectoires relles dans les vides du
milieu lchelle microscopique. Elle introduit le concept de dispersion.
Mise en vidence de la dispersion : Traages
Lexprience de Darcy est effectue sur une colonne de sable, verticale, en introduisant au
sommet de lappareil, un traceur une concentration Co (poids de traceur par unit de volume
de solution). Cette opration, appele traage, permet de mesurer , sur le terrain, la vitesse de
dplacement , la direction relle de lcoulement et les paramtres de la dispersion.
Deux mthode dintroduction des traceur sont utilises: injection massive, ou bouffe, de
courte dure ou continu concentration constante de longue dure.
La concentration C est mesure en bas de la colonne, des intervalles de temps chelonnes

ou en continu, par des enregistreurs automatiques. Le temps coul entre lintroduction du
traceur et sa dtection la sortie, est appel temps de sjour.
Les donnes obtenues sont portes sur un graphique. En ordonnes, les rapports des
concentrations C/Co. En abscisses les temps de sjour. La courbe obtenue est la courbe de
restitution du traceur. Elle dtermine les temps de sjour et la vitesse de dplacement.
Explication du phnomne de dispersion
La courbe de restitution montre que les particules de traceur, donc les particules de leau,
injectes un instant donn, au point de dpart, narrivent pas simultanment en bas de
colonne. la sortie elles sont tales dans le temps et dans un volume plus ou moins grand.
Ce fait nest pas conforme la loi de Darcy. Ce phnomne est appel la dispersion.
Ltude lchelle microscopique montre que les particules deau se dplacent dans les vides
continus, aligns selon la direction moyenne, gnrale, dcoulement. Elles dcrivent des
trajets compliqus. Cest la tortuosit des trajectoires. Au cours de ces trajets, les
caractristiques physiques du milieu entrainent des variations de la vitesse des molcules,
causes de la dispersion mcanique.
Cette action mcanique nest pas la seule qui intervient. La dispersion est due trois groupes
de facteurs:

La structure physique du rservoir, cause de la dispersion mcanique tudie

prcdemment, laquelle est prpondrante.
La structure du fluide dont lagitation thermique des molcules provoque la diffusion

molculaire.
Les interactions eau/roche lorigine de ladsorption et de la dsorption.
Dtermination de la vitesse de dplacement sur le terrain
Ltude de lcoulement de leau souterraine, vhicule de transport de toutes substances

minrales ou organiques, ncessaire pour la prvention contre la pollution de lespace
souterrain, doit considrer les trajectoires relles. Do la mise en uvre de mthodes de
mesures, sur le terrain, des vitesses de dplacement et des paramtres de la dispersion. Elles
reposent sur la technique des traages.
Porosit cinmatique
Cest le rapport de la vitesse de dplacement la vitesse de filtration. Il quivaut au rapport

du volume des vides rellement parcouru par leau gravitaire au volume total du milieu (
satur ou non). Cest la teneur en eau mobile. Ce concept est proche de la porosit efficace,
dfini comme un rapport de volumes.
Conditions aux limites
Types de conditions aux limites gologique et hydrodynamiques
Lidentification du comportement de laquifre repose sur une dfinition rigoureuse des

conditions aux limites. Leur tude porte sur leur situation et leurs conditions. Elles sont
dtermines par des points (sources) des lignes (berges des rivires) ou des surfaces (aire
dalimentation).
Ltude de la configuration de laquifre a dgag deux grands types de limites:
Limites gologiques, closes, position fixe ou limites tanches ( flux nul), imposes
par les structures hydrogologiques. Ce sont : le substratum, le toit, les passages
latraux de facis, les biseaux de transgression et les failles.
Limites hydrodynamiques, ouvertes, position variable dans lespace et dans le

temps, imposes par les conditions extrieures (environnement de laquifre).
Elles sont identifies un instant donn (calage des modles mathmatiques) ou au cours
dune dure moyenne ( prvisions).
Elles sont classes en trois types :
Limites flux impos ou conditions de dbit. Les dbits peuvent tre nuls, entrant ou
sortant. Les dbits nuls sont imposs par les limites gologiques tanches. Les dbits
entrants ou affluant, sont les nappes affluentes, les aires dalimentation par infiltration des

prcipitation efficaces, les rivires infiltrantes, etc. les dbits sortants sont les sources et
lignes dmergences, les cours deau drainants, etc.
Limites potentiel impos ou conditions de potentiel. Elles sont identifies par une
courbe quipotentielle ou hydroisohypes de la surface pizomtrique. Ce sont les lignes
de sources, les plans deau de surface (rives des lacs et rivires), les lignes de rivage.
La surface pizometrique ou surface libre, laquelle rpond deux conditions

particulires: pression gale la pression atmosphrique sur toute sa surface et flux nul.
La drainance est impose par des conditions de flux (limite gologique semi-permable)
et de potentiel.
Ligne de courant et lignes quipotentielles: Rseau dcoulement
La schmatisation de lcoulement de leau dans un aquifre est ncessaire pour dresser des
coupes et des cartes hydrogologiques. En hydrodynamique souterraine, lcoulement est
considr comme dplacement de particules le long de trajectoires thoriques. Ces trajectoires
sont matrialises par les lignes de courant. Elles sont synonymes de lignes de flux et filets
liquides. Une ligne de courant est donc un ligne idale qui reprsente la trajectoire dune
particule deau en mouvement dans un aquifre.
La perpendiculaire aux ligne de courant successives est la ligne quipotentielle. Cest une
ligne dgale potentiel hydraulique assimile une ligne dgale charge et dgal niveau
pizomtrique.
Les modalits de lcoulement souterrain, dans un aquifre, sont schmatises, en coupe ou en

plan, par un quadrillage de lignes de courant et de ligne quipotentielles. Celui-ci constitue un
rseau dcoulement
Schmatisation de lcoulement de leau souterraine dans un aquifre de subsurface:

systme de flux
Dans un aquifre de subsurface, lcoulement de leau souterraine seffectue des zones de

recharge vers les zones de dcharge.
Les zones de recharge ou aires dalimentations, sont constitues par les reliefs,
vritables chteaux deau. Les eaux souterraines y sont alimentes par linfiltration des
prcipitation efficaces.
Les zones de dcharge ou drainage, se localisent dans les valles gnralement occupes
par des cours deau, les surfaces deau libre ( lacs, mers, ocans) et les dpressions
endoriques des zones arides.
Le moteur de lcoulement est la diffrence daltitude H( diffrence de charge) engendrant

des gradients hydrauliques. Les limites hydrodynamiques des systmes sont constitus par les
lignes de partage des eaux souterraines sous les reliefs et les zones ou axes de drainage de
surface. Lunit, ainsi dfinie, constitue un systme de flux. Elle correspond un ou plusieurs
aquifres.
Deux grands types de systme de flux sont identifis:
Flux latral des zones de recharge (alimentation) vers celles de dcharge (drainage et
exutoires) : il comporte trois branches : descendantes subverticale, latrale et ascendante. Il
est identifi par lanalyse de la surface pizomtrique qui permet de reconnaitre les zones
dalimentation et dexutoires et les grands axes dcoulement.
Flux de drainance de direction subverticale et de sens ascendant ou descendant. Il est

prpondrant sur le flux latral dans les aquifres profonds, par les vitesses et les changes.
Zonalit verticale des aquifres dans un bassin hydrogologique
Trois grandes zones de systmes de flux. Soit de haut en bas:
Zone des aquifres de subsurface systme de flux locaux courts. Cest le domaine des
aquifres nappe libre et des premiers aquifres multicouches qui y sont troitement lis. Ils
sont imposs par la topographie locale, le rseau hydrographique et les surfaces deau libre.
Le rle de la gologie structurale est faible. Linfluence des zones climatiques actuelles est
importante. Les vitesses effective sont de lordre du Km par an pour les nappes libres et de
100 mtres pour les nappes captives. Cette zone atteint la profondeur de 50 100 m en
fonction des structures hydrogologiques et de lalimentation par la surface.
Zone des aquifres intermdiaires systme de flux rgionaux, longs. Les aquifres sont du
type nappe captive. Ils sont imposs par la topographie rgionale, les grands axes
hydrographiques, les grands lacs, ocans et mers intrieures et les dpressions endoriques
des zones arides. Le rle de la gologie structurale est prpondrant. Les vitesses effectives
sont de lordre de 10 m/an. Sa profondeur atteint 200 300 m en fonction des structures
hydrogologiques et des changes avec les aquifres de la zone prcdente.
Zone daquifre profonds systme de flux globaux trs longs. Les changes verticaux
dominent sur les flux latraux. Le rle de la gologie est prpondrant. Les zones climatiques,
de faible influence, cdent la priorit aux paloclimats. La profondeur dexploitation des
aquifres, pour les usages humains est limite par la minralisation croissante de leau
souterraine . Les vitesses effectives sont de lordre du mtre par an.
Zone hydrogologiques des aquifres dans un bassin hydrogologique
Trois zones hydrogologiques peuvent tre reconnues en fonction de la distance aux

affleurements. Elles sont identifies par les caractristiques gologiques, hydrodynamiques et
hydrochimiques avec apport des isotopes du milieu.
Zone dalimentation ou de recharge qui correspond gnralement aux limites gologiques

qui peuvent tre constitu entre autres, des affleurements de laquifre, et des aires de
drainance descendante.

Zone dexutoire ou dcharge, pouvant tre soit artificielle due lexploitation intensive par
les pompages, ou naturelle , constituant le reflet , en profondeur, du grand axe de drainage
impos par les cours deau sur la surface pizomtrique des aquifres de subsurface. Les
diffrences de charge ainsi cres, provoquent une drainance ascendante active.
Zone de transition qui marque le passage progressif des zones dalimentation celles
exutoires. Le flux latral diminue progressivement pour devenir trs faible dans les zones
dexutoire. Le flux de drainance augmente, imposant des changes verticaux ascendants,
rarement descendants.
IV. Techniques de forage et pompage dessai
Dfinitions et concepts de base
Dfinitions et concepts de base
Les exprimentations sur le terrain, par des puits et sondages, sont excutes par des essais de
puits et des pompages dessais. Elles consistent mesurer laccroissement des rabattements
du niveau pizomtrique en relation avec le temps de pompage et leur remonte aprs arrt de
lopration.
Ce sont des tests portant sur les modifications du comportement hydrodynamique du

complexe aquifre/ouvrage de captage, en rponse une impulsion cre par un pompage
dbit constant. Leur excution doit tre conduite avec la mme rigueur scientifique quune
exprience de physique.
Les essais par pompage poursuivent quatre buts, dans lordre croissant de complexit:
Dtermination des caractristiques du complexe aquifre/ouvrage de captage. Cest

lessai de puits, destin lquipement technique de louvrage.
Mesure sur le terrain des paramtres hydrodynamiques de laquifre: transmissivit et

coefficient demmagasinement.

tude quantitative des caractristiques particulires de laquifre: test des conditions

aux limites, structure, htrognit, drainance, etc.
Observation directe, en vraie grandeur , de leffet de lexploitation sur laquifre.

Prvisions de lvolution des rabattements en fonction des dbits pomps. valuation
de la ressource en eau souterraine exploitable.
Les trois dernires oprations sont conduites par le pompage dessai sur station de pompage
comportant au moins un pizomtre.
quipements techniques des puits et sondage
La colonne ascensionnelle
Constitue dun tube unique ou dlments tlescopiques, soutient la paroi du trou. Lespace
annuaire, entre le tubage et le terrain, est obstru par une colonne de ciment. Celle-ci joue un
double rle: consolidation de louvrage et suppression des fuites et intercommunications entre
aquifres. sa base elle est ancre dans le toit de laquifre nappe captive ou dans le
substratum de la nappe libre. Son diamtre est calcul pour le logement de la pompe et en vue
de limiter la perte de charge quadratique.
La partie captante
Comporte une crpine et, ventuellement, un massif filtrant. La crpine est un tube perfor
douvertures de formes diverses, travers lesquelles leau pntre dans le sondage. Dans les
terrains meubles lespace annulaire , entre la crpine et le terrain, est rempli de gravier calibr.
Son rle est double: filtre retenant les lments fins et augmentant la permabilit au
voisinage du sondage et soutnement du terrain. Diamtres des ouvertures et granulomtrie du
gravier sont calculs par des expressions empiriques introduisant le diamtre caractristique
des grains du rservoir.
Lensemble , partie captante et aquifre au voisinage immdiat du sondage, constitue un

complexe aquifre/ouvrage de captage.
Effets du pompage sur laquifre. Cne de dpression
Le pompage dans un aquifre, dont la surface pizomtrique initiale est suppose horizontale,
cre une dpression en forme dentonnoir. Son axe concide avec celui de louvrage. Cest le
cne de dpression.
Dans laquifre nappe libre, il affecte le rservoir tandis quil est fictif dans celui nappe
captive. Lexprimentation a pour but de mesurer, dbit constant, les dimensions de ce cne
un instant donn et leur volution dans le temps. Elle porte galement sur son effacement
aprs arrt du pompage, appel remonte.

Les deux donnes gomtriques du cne de dpression, un instant donn, t, sont:
Le rabattement, not s, mesur par labaissement du niveau pizomtrique dans le

puits de pompage ou dans un pizomtre implant une distance x, de laxe de
louvrage. Le plan deau , dans louvrage est le niveau dynamique. La profondeur du
niveau dynamique, au dessous du niveau pizomtrique initial, en rgime non
influenc, est le rabattement, mesur au cours de la remonte, est appel rabattement
rsiduel sr .
Le rayon dinfluence, not R, est la distance de laxe du puits laquelle le rabattement

est nul ou ngligeable. Dans la pratique ou il nest plus mesurable.
La donne mesure, un instant donn, sur le terrain, est donc le rabattement ou la profondeur
du niveau deau, soit dans le puits, soit dans un pizomtre.
Le cne de dpression est reprsent en coupe par une courbe de dpression. En plan par des
courbes dgal rabattement, cercles concentriques laxe du puits. Ces courbes sont
assimilables des lignes quipotentielles. Les lignes de courant convergent vers laxe du
puits. Le rseau dcoulement, ainsi trac, caractrise une nappe radiale convergente.

Essai de puits par paliers de dbit de courtes dures
Cet essai value les caractristiques du complexe aquifre/ouvrage de captage. Ce sont : le

dbit critique, le dbit spcifique, le dbit spcifique relatif, les pertes de charge dans
louvrage et son environnement immdiat et le dbit maximum dexploitation ou productivit.
Il permet dtablir le programme dquipement technique de louvrage: tubage, crpine et
massif filtrant, puissance de la pompe, etc.
Conditions de base
Les conditions de base dapplication des expressions dhydrodynamique souterraine en

rgime transitoire, auxquelles doit satisfaire le complexe aquifre/ouvrage de captage sont:
Validit de la loi de Darcy ( coulement laminaire, isotropie, homognit et
continuit)
Puits complet: captant toute lpaisseur de laquifre, atteignant le substratum et
crpin sur toute sa hauteur.
Puits correctement dvelopp et quip
Surface pizomtrique subhorizontale
Dbit de pompage constant
Rayon du puits le plus petit possible
Excution de lessai de puits. Paliers de dbit
Lessaie de puits est effectu en ralisant des paliers de dbit, dbit constant pendant une
courte dure dtermine, 1 3 heures. Il mesure deux donnes: le rabattement et le dbit
constant. Chaque palier de dbit est suivi dun arrt de pompage dune dure gale,
permettant la remonte du niveau deau et la mesure du rabattement rsiduel. La premire
remonte doit tre poursuivi jusqu atteindre approximativement le niveau pizomtrique
initial. En gnral le temps de pompage du premier palier est suffisant.
Les dures gales de pompages et darrts sont courtes, une trois heures au maximum. Le
dbit initial est gal celui de la puissance minimum de la pompe. Ensuite, les dbits
croissent selon une progression de 2, 3,4 . Le nombre de paliers de dbit est au minimum de
quatre, le dernier tant dbit lev suprieur au dbit critique. Il est prfrable de choisir des
dures courtes et den augmenter le nombre jusqu un optimum de six. Pour les aquifres
nappe captive, trois paliers de dbit peuvent tre suffisants. Il doit tre suprieur pour les
aquifres nappe libre.

Lors de la mise en route du pompage, les premiers volumes deau exhaurs correspondent la
vidange de louvrage, donc un coulement quadratique non linaire. Laquifre nest
sollicit quaprs un certain dlai. Cest leffet de capacit du puits.
Lexcution du premier palier de dbit doit tre prcde dune dure de pompage tc suffisante
pour leffacer. Elle est fonction de la transmissivit T de laquifre et du volume contenue
dans louvrage.
Signification du rabattement dans louvrage: pertes de charge

Le rabattement, mesur dans louvrage instant t, est la somme de deux composantes,
nommes pertes de charge, exprimes en mtres de hauteur deau, caractrisant le complexe
aquifre/ouvrage de captage:
Perte de charge linaire provoque par lcoulement laminaire de laquifre au
voisinage du puits note BQ.
Perte de charge quadratique, non linaire, provoque par lcoulement turbulent dans
louvrage, crpine et tubage, note CQ.
Le rabattement total s linstant t est ainsi donn par lexpression de C.E. Jacob (1946):
s = BQ + CQ
Cette expression, la plus utilise, tablie pour laquifre nappe captive, est tendue
laquifre nappe libre sous condition que le rabattement mesur soit infrieur 0,1*b
Perte de charge linaire

La perte de charge BQ rsulte de deux effets:

Influence de laquifre o lcoulement laminaire est de rgime transitoire. La perte de

charge linaire, impose par les paramtres hydrodynamiques de laquifre, au
voisinage du puits, crot avec le temps de pompage. Elle est caractristique de
laquifre.
Influence de la partie captante du puits: remaniement du rservoir, crpine et
ventuellement massif filtrant. Leffet est positif si ce dispositif apporte une
amlioration de lcoulement. Il est ngatif en prsence de colmatage.
Perte de charge quadratique

Note CQ , elle est fonction uniquement du dbit. Elle est caractristique de lquipement
technique de louvrage, diamtre du tubage et de la crpine.
Pour de faibles dbits, impliquant une faible vitesse, les pertes de charge dans le tubage et
dans la crpine sont linaires ou ngligeables. Lquation se simplifie:
s = BQ
Lorsque la vitesse effective dans laquifre est suprieure la vitesse critique, lquation
dvient:
s = CQ
Ce terme inclut alors toutes les pertes de charge dans laquifre et dans louvrage.
Interprtation graphique des donnes de lessai de puits

Relations dbit / temps et rabattement / temps
Sur un papier graphique coordonnes linaires sont ports, en ordonnes le dbit ou les
rabattements, et en abscisses les temps. Deux graphiques sont obtenus:
Graphique dbits/temps de pompage, figurant les dures et les dbits des paliers de
dbit, les dures et les arrts de remontes.
Courbe rabattements/temps de pompage et rabattement rsiduels/temps de remonte.
Ces graphiques permettent de contrler le bon droulement de lexprimentation dessais de
puits.
Courbe dbits/rabattements. Dbit critique

Le couple de donnes de chaque palier de dbit, dbit constant, en m3/h et rabattement
rsiduel en m, est port sur un papier graphique linaire .
Les points obtenus tracent la courbe
dbits/rabattements ou courbe caractristique, reprsentant la fonction
s= f(Q).
Sa forme apporte des informations sur le comportement hydrodynamique du complexe
aquifre/ouvrage de captage lorigine de la perte de charge quadratique:

Une droite : Perte de charge nulle ou ngligeable

Une courbe convexe : Perte de charge importante
Une courbe concave : Un essai de puits non valable
La courbe dbits/rabattements prsente deux parties distingues par le point critique A,
correspond en abscisse au dbit critique Qc.
Droite dbits spcifiques/rabattements : Dbit spcifique et dbit spcifique relatif.
Le dbit spcifique dun puits, not qs, est le dbit pomp Q rapport au rabattement s, dans le
puits, dans des conditions dessai de puits dfinies.
qs Q
=
sr
La relation dbits spcifiques/rabattements est reprsente par une droite incline,

reprsentative de la fonction, s= f (Q). Cette droite permet de calculer le dbit spcifique
relatif, ou dbit constant pomp dterminant un rabattement unitaire (1mtre) la fin dun
palier de dbit.
Sur la droite dbit spfiques/rabattements, labscisse correspondant lordonne 1 mtre,

donne la valeur du dbit spcifique relatif.

Droite dbits/rabattements spcifiques

Le rabattement spcifique s/Q, est la hauteur de rabattement mesure dans le puits rapporte
au dbit pomp dans des conditions dessais de puits dfinies (paliers de dbit). Il est exprim
en m/m3.h.
Lquation de Jacob peut scrire s/Q= B +CQ
Cest lquation dune droite, sous rserve que les mesures de chacun des dbits soient
effectues au mme instant. Cette droite met en vidence certaines formulations simples de la
relation dbits/rabattements. Quatre cas peuvent se prsenter:
1. Droite passant par lorigine indiquant que le rgime turbulent est fortement prdominant
dans laquifre et dans le puits s = CQ
2. Droite ne passant par lorigine s= BQ + CQ
3. Droite pente nulle, verticale parallle laxe des ordonnes, traduisant un coulement
laminaire, avec pertes de charge dans la crpine et le tubage nulles ou ngligeables s = BQ
4. Courbe concave vers le haut s = BQ + CQn ; avec n = 3 , 4, etc
Calcul des pertes de charges

La droite dbits/rabattements spcifiques, permet de dterminer les coefficients B et C de
lquation s/Q = B + CQ
Le coefficient B est obtenu par lintersection de la droite reprsentative avec laxe des
rabattements spcifiques.
Le coefficient C est gal la pente de la droite reprsentative.

Le rabattement, dtermin linstant t, imputable la perte de charge linaire,

consquence de lcoulement laminaire de laquifre BQ, impose par les paramtres
hydrodynamiques, est donne par lexpression: S = BQ = 0,01 * Q
Les valeurs calcules donnent la droite reprsentative. Pour chaque dbit, la diffrence
entre le point de cette droite et celui de la courbe de dbits/rabattements, situ au-
dessous, donne la valeur de perte de charge quadratique CQ.
Dtermination de la productivit dun puits. Dbit dexploitation maximum

La productivit dun puits Pr, est le dbit maximum qui peut tre pomp dans louvrage,
pendant une dure dfinie, sans que le rabattement induit par le pompage ne dpasse le
rabattement maximum admissible.
Le rabattement maximum admissible est impos par :
Des contraintes physiques et techniques du complexe aquifre/ouvrage de captage,

exprimes par le dbit critique Qc et le rabattement critique sc correspondant mesur
par les essais de puits.
Des contraintes socio-conomiques, dont la principale est le cot de production de

leau, imposant la profondeur du niveau dynamique.
Le rabattement maximum retenu doit donc tre gal au rabattement maximum mesur
sans dpasser le rabattement maximum admissible.
Pompages dessais
l'inverse des essais ponctuels qui ne caractrisent que des horizons permables trs limits,
les pompages dessai sintressent toute la hauteur du rservoir aquifre (puits complets)
ou une paisseur suffisamment reprsentative de la nappe (puits incomplets). Ils
fournissent des informations variables, en fonction de leurs conditions
de ralisation et de la nature du dispositif de suivi.
But du pompage dessai
Le pompage dessai poursuit trois buts principaux:
Mesure sur le terrain des paramtres hydrodynamiques : transmissivit et

coefficient demmagasinement;
tude quantitative des caractristiques particulires de laquifre: conditions aux

limites (confirmation de la distance du puits la limite, colmatage des berges dune
rivire), structure (htrognit, drainance);

Observation directe, en vraie grandeur, de leffet de lexploitation sur laquifre.

Prvision de lvolution du rabattement en fonction des dbits pomps. valuation de
la ressource en eau souterraine exploitable.
Droulement du pompage dessai

Les pompages dessai de longue dure sont excuts par un seul palier de dbit, dbit
constant, prolong durant au moins 42 heures; avec un optimum de 72 heures. La remonte
des niveaux doit tre observe pendant une dure gale. Lexcution et linterprtation des
donnes mesures, rabattements et temps, reposent sur lemploi des expressions
dhydrodynamique en rgime transitoire.
La mise en uvre dun pompage dessai ncessite linstallation pralable de dispositifs

adapts pour permettre un suivi des paramtres recherchs :
Dans louvrage asservi au pompage, un tube pizomtrique, pour contrler le niveau

de la nappe par sonde manuelle ou par pizographe ;
En sortie du tube de refoulement, une vanne, un compteur et un robinet,

permettant de rgler le dbit dexhaure, de mesurer le volume pomp et de procder
des prlvements deau ;
En priphrie et sur des distances variables en fonction du contexte et des objectifs

de lessai, des quipements destins contrler les effets loigns du pompage
(pizomtres, jaugeage de source, chelle limnimtrique sur vallon).
Dans certains cas, il savre ncessaire dvacuer les dbits exhaurs dans un rseau
spcifique, soit pour viter les rinjections deau dans la mme nappe, soit pour protger
des fonds infrieurs sensibles (risques de glissement de terrain, par exemple) ou trs
urbaniss.

Il est enfin indispensable que les oprations de pompage ne soient pas influences par la
pluviomtrie ou par dautres prlvements dans la mme nappe sur des sites proches,
afin de valoriser linterprtation des rsultats obtenus.
Interprtation graphique des pompages dessai
Relation entre les rabattements et les temps : droite reprsentative
la rsolution des expressions dapproximation logarithmique est obtenue par le trac et

linterprtation de la droite reprsentative rabattement/logarithmes des temps de remonte.
Le premier terme dans ces expressions est une constante avec Q et T constants. Le second
terme, seul le temps varie.
Les rabattements croissent en fonction du logarithme du temps de pompage. Cette condition

est conforme au concept de rgime transitoire.
Les donnes de pompages sont reportes sur un papier graphique semi-logarithmique.
Le niveau pizomtrique initial est indiqu en haut du graphique.
Les points obtenus tracent la droite moyenne reprsentative.
La courbe observe, au dbut du pompage, traduit leffet de capacit de louvrage,

provoquant un coulement turbulent non linaire.
Le point dintersection de la droite reprsentative avec le niveau pizomtrique

initiale, mesure le temps fictif lorigine, not t0.
Niveau
Module
Effet de puits logarithmique
coulement
turbulent
Temps pompage

Calcul des paramtres hydrodynamiques
La transmissivit est calcule par la pente de la droite reprsentative. Lchelle des

coordonnes ntant pas homogne, la pente est dtermine par laccroissement des
rabattements (ou des profondeurs des niveaux deau), au cours dun
module logarithmique, not c.
La transmissivit est calcule par lexpression:
Le coefficient demmagasinement est obtenu par calcul numrique dans le deuxime terme de
lexpression
ou plus simplement, avec t0, lorsquil peut tre dtermin, par

lexpression:
Relation entre les rabattements rsiduels et les temps de remonte de niveaux
Les donnes de la remonte des niveaux, aprs arrt du pompage, permettent dtablir le
graphique traant la droite reprsentative de lexpression de Jacob.
Les rabattements rsiduels, exprimes en m (ou les niveaux deau), sont ports en ordonnes
linaires, la valeur (t+t)/t en abscisses logarithmiques. Seule la transmissivit peut tre
calcule avec lexpression :
Types hydrodynamiques daquifres et application de lexpression dapproximation

logarithmique
Lapplication de lexpression dapproximation logarithmique de Jacob est limite au type

hydrodynamique daquifre nappe captive illimit et pontes (substratum et toit)
impermables (absence de drainance). Toutefois, lapplication peut tre tendue, sous
certaines conditions, laquifre nappe libre illimit.
Interprtation des pompages dessai : conclusion gnrale
Pour interprter un pompage dessai, le couple de donnes recueilli est port sur un papier
graphique semi-logarithmique:
Les rabattements ou les profondeurs de leau en m en ordonnes linaires.
Les temps de pompage ou (t + t)/t pour la remonte, en abscisses logarithmiques.

Lessai est valable sil est possible de tracer une droite moyenne reprsentative. Trois cas
types peuvent tre obtenus. Ils sont interprtables si le type hydrodynamique daquifre est
identifi par ltude hydrogologique.
1. Droite reprsentative pente constante: aquifre nappe captive ou libre (s infrieur

0,1b), illimit pontes impermables.
2. Droite reprsentative brise avec une pente double aquifre nappe captive ou libre,
limit latralement par une limite tanche: passage latral de facis ou faille. Dans ce cas,
La distance thorique d du puits de pompage la limite dalimentation est calcule par
lexpression:
Avec, x tant la distance du pizomtre laxe du puits et ti le temps dintersection
3. Droite reprsentative brise avec palier de stabilisation, indice de dbit entrant aux limites.
Deux cas:
Aquifres, nappe captive ou libre, limits latralement par une limite potentiel
impos (limite dalimentation): cas du systme global aquifre/rivire.
Aquifre nappe semi-captive, impliquant la drainance par une ou deux pontes.
Dans tous les cas, le type hydrodynamique daquifre tant identifi, le premier segement de
droite donne, lexception de la drainance, la solution numrique de lexpression
dapproximation logarithmiques de Jacob. La transmissivit et le coefficient
demmagasinement ( descente seulement) peuvent tre calculs.

V. Cartographie de laquifre et cartes pizomtriques

La cartographie de laquifre a pour but de reprsenter sa configuration, sa structure et de
schmatiser les fonctions du rservoir et son comportement hydrodynamique.
Les cartes sont de deux types : structurales (morphologie, positions des surfaces limites,
paisseurs) et pizomtriques (tude hydrogologique).
La synthse et l'interpolation spatiale des donnes sur les paramtres du rservoir, sont
prsentes par les cartes structurales. Celles des niveaux pizomtriques et des conditions aux
limites par les cartes pizomtriques. Elles schmatisent les fonctions capacitives et
conductrices du rservoir et le comportement hydrodynamique de l'aquifre. Elles sont
indispensables l'tablissement des modles physiques et mathmatiques.
CARTES STRUCTURALES DE L'AQUIFERE
Les cartes structurales de l'aquifre reprsentent sa configuration et sa structure. Elles sont

tablies par synthse des donnes sur la gologie, les conditions aux limites et les
paramtres physiques et hydrodynamiques des aquifres.
Cartes de la configuration de l'aquifre: Dimensions de l'aquifre
La cartographie de la configuration, ou l'enveloppe de l'aquifre, reprsente les limites

gologiques et hydrodynamiques. Ces cartes permettent de dterminer les dimensions et le
volume de l'aquifre.
L'interpolation spatiale des donnes ponctuelles est figure par trois types de cartes en
courbes d'isovaleurs :
Cartes en courbes isohypses ou d'gale altitude, figurant la morphologie de la

surface considre, au mme titre que les cartes topographique en courbes de niveau,
reprsentent la surface du sol ; ex : Cartes de la surface du substratum
Cartes en courbes isobathes ou d'gale profondeur, par rfrence la surface du sol,

situant dans le sous-sol la surface reprsente ; ex : Cartes de la limite suprieure de
l'aquifre
Cartes en courbe isopaches ou d'gale paisseur de l'aquifre, bases du calcul du

volume du rservoir ; ex : Cartes de l'paisseur de l'aquifre
Les conditions aux limites latrales gologiques de l'aquifre sont portes sur les cartes
structurales. Celles aux limites hydrodynamiques figurent sur les cartes pizomtriques.
Leur cartographie exacte, avec indication du type et de donnes numriques, est la base de
l'tude du comportement hydrodynamique de l'aquifre.

CARTES PIEZOMTRIQUES
Les cartes de la surface pizomtrique, dites pizomtriques, tablies avec les donnes sur les
niveaux pizomtriques, reprsentent, une date donne, la distribution spatiale des
charges et des potentiels hydrauliques. Elles figurent galement les conditions aux limites
hydrodynamiques. Des cartes des fluctuations de la surface pizomtrique des nappes libres,
dans l'espace et dans le temps, sont galement tablies. Les cartes pizomtriques sont les
documents de base de l'analyse et de la schmatisation des fonctions capacitive et
conductrice du rservoir et du comportement hydrodynamique de l'aquifre. C'est la
synthse la plus importante d'une tude hydrogologique.
Etablissement des cartes pizomtriques
Il repose sur :
1. Mesure des niveaux pizomtriques
Elles doivent tre effectues dans des conditions de stabilisation et pour l'ensemble de la
rgion cartographie au cours d'une priode la plus courte possible. En effet ce document a
une valeur de rfrence une date donne. En cas de variations importantes au cours de la
campagne de relevs, il faut effectuer des corrections en rapportant les rsultats une cote de
rfrence d'un (ou de plusieurs) ouvrage reprsentatif en observation continue par
enregistrement limnigraphique.
2. Report des niveaux pizomtriques: Echelle de la carte
Les points d'eau, affects de leur code de rfrence et de leur niveau pizomtrique, sont
reports sur une carte topographique en courbes de niveau grande chelle, en gnral 1/50
000. L'chelle de la carte est choisie en tenant compte de la densit des points de mesure et
des fonds topographiques existants. La priorit est donne la prcision du nivellement. La
date, laquelle ont t effectues les mesures, est porte sur la carte.
Trac des courbes hydroisohypses
La surface pizomtrique est, comme la surface du sol, reprsente par des courbes d'gale
altitude de niveau d'eau, soit d'gal niveau pizomtrique, dites courbes hydroisohypses. Le
dessin de ces courbes comporte successivement le choix de leur quidistance et la technique
de leur trac
a. Choix de l'quidistance des courbes hydroisohypses
L'quidistance des courbes hydroisohypses est la distance constante entre des plans
horizontaux d'gal niveau pizomtrique. Leur intersection avec la surface pizomtrique est
l'espacement qui se mesure sur la carte pizomtrique entre deux courbes successives.
L'quidistance dpend de la prcision et de la densit des mesures, des valeurs du gradient
hydraulique, de l'chelle de la carte et de la qualit du nivellement. En gnral elle est de
l'ordre du mtre (0,5, 1 ou 2 m) pour les cartes 1/10 000 et 1/20 000 de 5 ou 10 m pour
celles l/50 000 et 1/100 000.
b. Technique de trac des courbes hydroisohypses
Il est effectu par trois mthodes d'interpolations, adaptes la prcision et la densit des
donnes disponibles. Une carte dresse avec des ajustements trop nombreux ne peut tre
utilise pour des tudes srieuses car les conclusions dcoulent des hypothses de travail.
L'interpolation approximative des niveaux pizomtriques est effectue par une

mthode visuelle. Les courbes sont traces en tenant compte, implicitement, des lois
gnrales de la morphologie de la surface pizomtrique. Dans la plupart des cas
cette mthode donne des rsultats satisfaisants. Mais elle doit tre utilise avec
prudence car il faut viter, lors de l'analyse ultrieure, de confondre les hypothses
de travail avec les faits observs.
Pour le trac par la mthode d'interpolation du triangle, les donnes sont groupes
par trois aux sommets de triangles. Les cts du triangle sont tracs et diviss en
segments proportionnels. L'quidistance retenue dans cet exemple est de 0,5 m. Les
courbes hydroisohypses sont obtenues en joignant, par des segments de droite, les
points d'gal niveau. Les tracs sont lisss pour obtenir des courbes rgulires.
Dans les secteurs dpourvus de points de mesure, les courbes en tirets sont ajustes
sur celles qui les encadrent en amont et en aval. Les courbes maitresses, multiples de
5, 10, 50 ou 100,sont soulignes en traits forts. Cette mthode donne d'excellents
rsultats lors que les points de mesures sont suffisants.
3. Interprtation des cartes pizomtriques
L'interprtation des cartes pizomtriques, appuye sur les cartes structurales du rservoir,
aboutit cinq oprations :
Analyse morphologique de la surface pizomtrique;
Etude de la structure de l'aquifre, des anomalies structurales du rservoir et de la

distribution spatiale des paramtres hydrodynamiques;

Etude des fonctions du rservoir distribution spatiale des stocks d'eau et rgime de
l'coulement de l'eau souterraine;
Etude du comportement hydrodynamique de l'aquifre dbits imposs entrant et

sortant, potentiels imposs;
Analyse des fluctuations de la surface pizomtrique des aquifres nappe libre.

Prvision de l'volution des niveaux pizomtriques.
Les donnes obtenues sont valables une date dtermine, celle du recueil des informations
portes sur la carte. Des cartes de donnes moyennes sont galement tablies : tiage moyen
annuel, surface pizomtrique moyenne annuelle, etc.
Habillage de la carte pizomtrique
Pour faciliter l'interprtation d'une carte en premier lieu, il est utile de procder deux
oprations d'habillage :
Trac des lignes de courant de la surface pizomtrique, donc superficielles. Elles

matrialisent la direction moyenne de l'coulement. Ce sont les droites de plus grande
pente, donc la perpendiculaire leve sur chaque courbe hydroisohypse.
Flchage des lignes de courant indiquant le sens de l'coulement dduit des niveaux
pizomtriques.
Les axes principaux du flux, correspondant aux trajets les plus courts et les plus
simples, sont souligns. Ce sont galement les rayons de courbure des arcs
lmentaires successifs.
Les droites brises, ainsi traces, perpendiculairement chaque courbe successive

qu'elles recoupent, sont lisses. Il est alors possible de schmatiser les caractristiques
principales de l'coulement, donc de la fonction conduite du rservoir.
Les types d'aquifres sont identifis. Les lignes de partage des eaux souterraines
dlimitent les bassins hydrogologiques.
La courbure dun arc de cercle est identifie par son orientation et son rayon
L'orientation de la concavit, par rapport au sens de l'coulement, distingue deux types
d'arcs de cercle:
L'orientation aval avec concavit ouverte vers l'aval. Les lignes de courant
convergentes identifient un aquifre nappe convergente.
L'orientation amont avec concavit ouverte vers l'amont. Les lignes de courant
divergentes identifient un aquifre nappe divergente.

La gomtrie des courbes hydroisohypses conduit des reprsentations parfois

complexes, qui rsultent de la combinaison de quelques formes lmentaires simples :
Continuit linaire avec des lignes de courant rectilignes et parallles, ce qui

traduit un coulement relativement uniforme ;
Courbure avec concavit tourne vers laval o les lignes de courant
convergent vers un axe de drainage privilgi ;
Courbure avec concavit tourne vers lamont et des lignes de courant
divergentes, ce qui matrialise une crte pizomtrique et caractrise souvent une
zone dapport par infiltration
Courbes fermes avec des lignes de courant convergentes. Ces dpressions
pizomtriques indiquent des fonctions dans la nappe par pompages ou par fuites
vers un aquifre sous-jacent
Courbes fermes avec des lignes de courant divergentes. Il sagit alors de
dmes pizomtriques qui correspondent des aires privilgies dinfiltration.
Le module d'espacement des courbes, dans le sens de l'coulement est constant,

dcroissant ou croissant. Il caractrise, sur une coupe verticale passant par une ligne
de courant le profil pizomtrique. Il est numris par le gradient hydraulique. Le
module d'espacement identifie deux types d'coulement rpartis en trois grands types
d'aquifres lmentaires :

Ecoulements uniforme
L'coulement uniforme est caractris par des dbits unitaires et une direction constante en
tous points du domaine aquifre. Il se traduit, en coupe, par un profil pizomtrique
linaire et, en plan, par un module d'espacement constant. Ce rgime exige un aquifre
homogne paisseur constante.
Ecoulements non uniforme
L'coulement non uniforme prsente des dbits unitaires et une direction variable selon les
points du domaine aquifre. Il se traduit par deux types de profils pizomtriques
paraboliques avec d'espacement dcroissant et hyperbolique avec module d'espacement
croissant. Dans le premier type, les autres paramtres tant constants, le dbit de la nappe
dcroit dans le sens de lcoulement. Dans le second il croit.
Grands types d'aquifre lmentaires
La courbure et le module d'espacement des courbes hydroisohypses identifient trois grands

types d'aquifres lmentaires, et chaque type est caractris, en plan, par le trac des courbes
hydroisohypses et, en coupe, par son profil pizomtrique :
Les aquifres nappe plate;
Les aquifres nappe plate prsentent une surface pizomtrique plane, incline dans le
sens de lcoulement. Ce type apparait donc, sur la carte pizomtrique, par des droites
hydroisohypises parallles et module d'espacement constant. Les lignes de courant sont
rectilignes et parallles. Le profil pizomtrique est linaire.
Ce type de nappe, le seul caractriser un coulement uniforme, est restreint aux aquifres
homognes paisseur constante, donc pratiquement aux aquifres nappe captive.
Les aquifres nappe cylindrique;
Les aquifres nappe cylindrique sont caractriss par une surface pizomtrique
cylindrique. C'est--dire dont les gnratrices horizontales, orthogonales aux lignes de
courant, se confondent avec les droites hydroisohypses. Le rgime d'coulement est non
uniforme. Il apparait, en plan, par des droites parallles au module d'espacement variable.
Deux modules d'espacement caractrisent deux profils: dcroissant pour le profil parabolique
et cCroissant pour le profil hyperbolique nappe cylindrique est rare.
Les aquifres nappe radiale.
Les aquifres nappe radiale sont les plus frquents. En gnral la surface pizomtrique
une forme conique bombe (surface convexe) ou dprime (surface concave). Il en rsulte que
les lignes de courant sont des rayons divergents ou convergents. D'o deux types :

Aquifre nappe radiale divergente avec arcs de cercle orientation amont et

lignes de courant divergentes d'alimentation type de nappe convexe. Il
caractrise souvent les aires par infiltration des prcipitations efficaces;
Aquifre nappe radiale convergente avec arcs de cercle orientation aval et
lignes de courant convergentes. Habituellement il caractrise les zones de
drainage gnral par les cours deau.
Module despacement et paramtres hydrodynamiques
Dans les conditions dapplication de la loi de Darcy, une approximation satisfaisante au calcul
du dbit de la nappe:
Q=A.K.i
Avec K.b=T , la transmissivit
Q=T.L.i
Do,
i= Q / (A.K) = Q / (T.L)
La section sera considre comme une constante, ses variations pouvant tre identifies par
ltude hydrogologique. Les dbits Q1 et Q2, traversant deux sections constantes, A1 et A2,
successives de laquifre dans le sens de lcoulement, mais de coefficient de permabilit,
K1 et K2 ou de transmissivit, T1 et T2 diffrents, sont gaux.
Q1 = Q2 = A.K1.i1 = A.K2.i2
Do,
K1 / K2 = i2 / i1
Et avec T, L tant la largeur constante de la section :
Q1 = Q2 = L.T1.i2 = L.T2.i1
do,
T1 / T2 = i2 / i1
Le module despacement E, section constante, est fonction inverse du gradient hydraulique.

Do:
E1 / E2 = K1 / K2 = T1 / T2 = i2 / i1
Daprs lexpression i= Q / (A.K) = Q / (T.L) la pente du profil pizomtrique, mesure par

le gradient hydraulique, i est :
section constante fonction directe du dbit de la nappe et fonction inverse du

coefficient de permabilit.
largeur constante fonction directe du dbit de la nappe et fonction inverse de la
transmissivit.

Lexpression E1 / E2 = K1 / K2 = T1 / T2 = i2 / i1 montre que le module despacement est,

section constante, fonction directe des paramtres hydrodynamiques.
Etude de la structure de l'aquifre
Lanalyse morphologique de la surface pizomtrique, confronte avec les caractristiques

gologiques du rservoir, permet l'chelle rgionale :
d'identifier les htrognits et les anomalies structurales du rservoir .
d'valuer la distribution spatiale des paramtres hydrodynamiques coefficient de

permabilit ou transmissivit et dbit de la nappe.
Identification des anomalies structurales du rservoir
L'effet des variations latrales de facies tant identifies par la gologie, les principales
anomalies structurales du rservoir se traduisant dans la morphologie de la surface
pizomtrique, sont les variations de section, la surface du substratum et les accidents
tectoniques:
Les variations de la section portent sur l'paisseur ou la largeur sparment ou

conjointement. Elles se traduisent, les autres facteurs tant constants, directement par
une modification locale du module d'espacement : (1) Soit une diminution
(resserrement) accompagnant une rduction de section, (2) Soit un accroissement lors
d'une augmentation. Ces anomalies sont reconnues, avant ou aprs l'tude
hydrogologique, par la gologie structurale.
L'anomalie de la surface du substraturn la plus importante est sa dformation

provoquant surtout une variation de puissance, donc de section. Une dpression
augmente le module d'espacement et une protubrance le diminue. Les accidents
tectoniques, affectant le substratum, apparaissent dans la morphologie de la surface
pizomtrique.
Les accidents tectoniques du rservoir se traduisent par des anomalies de la surface

pizomtrique. C'est le cas de failles importantes. D'ailleurs ils sont parfois dcels par
l'hydrogeologie, laquelle apporte ainsi sa contribution la gologie structurale
Interprtation globale des carte pizomtrique , implantation douvrages de captages
L'interprtation globale des cartes structurales et pizomtriques aboutit l'identification des

zones privilgies pour l'implantation des stations des ouvrages de captages. Elle contribue
galement la prescription des mesures de protection de la qualit des eaux souterraines
captes pour l'alimentation humaine. Elle porte sur les variations locales des paramtres
tudis par rapport au contexte gnral. Elle dgage quatre grands types d'aquiferes d'intret
dcroissant :

type I. Aquifre nappe radiale convergente profil hyperbolique
type II. Aquifre nappe radiale divergente profil hyperbolique,
type Ill. Aquifere nappe radiale convergente profil parabolique;
type IV. Aquifre nappe radiale divergente profil parabolique
Analyse des fluctuations de la surface piezometrique
L'analyse de la surface pizomtrique des aquifres nappe libre aboutit l'tude de leur
comportement hydrodynamique dans l'espace la date des mesures de niveaux
pizomtriques. Celle des fluctuations introduit leur variabilit dans le temps.
Elle conduit trois ensembles de donnes :
Evolution du comportement hydrodynamique dans le pass (historique), base des

prvisions au cours d'tapes futures. Etude des variations dans le temps des conditions
aux limites hydrodynamiques et plus particulirement sur l'alimentation par infiltration
des prcipitations efficaces;
Evaluation de la variation de la rserve en eau souterraine ou rserve rgulatrice et

son volution dans le temps;
Prvision des niveaux pizomtriques d'tiage et des dbits d'tiage des cours d'eau
drainants.

VI. Conclusion
Les cartes structurales et pizomtriques de l'aquifre sont un moyen efficace d'tude de sa

configuration, de sa structure, de ses conditions aux limites, des fonctions de son rservoir
et de son comportement hydrodynamique. La carte pizomtrique, figurant les conditions
aux limites gologiques et hydrodynamiques, est le document essentiel de synthse d'une
tude hydrogologique. Elle est indispensable au calage des modles mathmatiques de
simulation hydrodynamique en rgime permanent.

Cours COMPLET Hydrogeologie BTP

Transféré par

Droits d'auteur :

Formats disponibles

Cours COMPLET Hydrogeologie BTP

Transféré par

Informations du document

Description originale:

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

Cours COMPLET Hydrogeologie BTP

Transféré par

Droits d'auteur :

Formats disponibles

Filire dIngnieurs:

Btiment & Travaux Publics

Mme Malki Mouna

Anne universitaire : 2016/2017

I. Gnralit, dfinitions et notions de base

Un aquifre est une formation hydrogologique permable permettant lcoulement

Hydrodynamique avec fluctuations libre : aquifre nappe libre

Gologique impermable : aquifre nappe captive

Gologique semi-permable : aquifre nappe semi-captive

Aquifre nappe libre

Formation gologique permable partiellement sature, limite par un plancher

Aquifre nappe captive

Aquifre nappe semi-captive

Anne universitaire 2016 / 2017 2

Aquifre, complexe unique rservoir / eau souterraine

Complexe physico-chimique unique de deux constituants essentiels, ou phase, troitement lis

Le rservoir: phase solide, milieu poreux ou fissur, constitue la trame de la structure,

Une combinaison de formations hydrogologiques semi-permables intercales entre des

1 - Fonction capacitive : emmagasinement de leau, il stocke et libre de leau)

2 - Fonction conductrice : il permet la circulation de leau

Anne universitaire 2016 / 2017 3

3 - Fonction dchange : interactions physico-chimiques permanentes entre leau souterraine

Zone non sature :

Zone o la totalit de la porosit est remplie deau. La circulation de leau se fait

Types deau dans un aquifre :

Eau de rtention (eau lie)

Eau adsorbe = hygroscopique : Lie la roche par attraction molculaire (adsorption)

II. Identification hydrogologique du complexe eau/rservoir

Anne universitaire 2016 / 2017 4

dimensions. Deux grands types de vides, pores et fissures, caractrisent respectivement le

Morphologie des pores et milieu poreux

Interconnexion des pores et milieu continu

Morphologie des fissures et milieu fissur

Anne universitaire 2016 / 2017 5

Porosit totale (n)

n = (Volume des vides / Volume total )*100 = (Vv / Vt )*100

Porosit efficace ou effective

Ou encore porosit de drainage correspond leau mobilisable ou exploitable libre sous

ne = Volume de leau gravitaire / Volume total )*100 = (Vg / Vt )*100

Indice des vides (e)

e = volume des vides / volume du solide = Ve / Vs

Capacit de Rtention ou coefficient de rtention (Cr)

n = ne + Cr = Vr / Vt * 100 et Cr = (Vr-Vg) / Vt * 100

Anne universitaire 2016 / 2017 6

Degr de saturation (Sr) :

Il caractrise le pourcentage d'eau contenu dans les vides :

Remarque : Ne pas confondre la porosit qui est la proprit du rservoir de stocker ou de

III. Ecoulements souterrains : Loi de Darcy

Loi de Darcy : applications

Enonc de la loi de Darcy

Q (m / s) = k(m/s) * A (m) * h/l

Le terme K est dfini comme un coefficient, dpendant de la permabilit de la couche ,

Anne universitaire 2016 / 2017 7

Do avec i= h/l, lexpression devient :

Le dbit unitaire q est le dbit en m3/s traversant lunit de section, perpendiculaire la

Do, en combinant les deux dernires expressions, on obtient:

La vitesse de filtration V en m/s, rapporte la section totale A est :

V (m/s) = Q/A = K*i= q

Dispositif de laboratoire avec coulement latral

n = (Volume des vides / Volume total )100 = (Vv / Vt )100

ne = Volume de leau gravitaire / Volume total )100 = (Vg / Vt )100

Lexpression de la loi de Darcy , Q= KAi , devient avec A=b*L :