MODELISATION HYDROGEOLOGIQUE

DU SYSTEME AQUIFERE
DE LA PLAINE DES GALETS


















FEVRIER 2012



Date Etablissement Nombre de pages
Fvrier 2012 Office de lEau 39
TABLE DES MATIERES


I CONTEXTE ET OBJECTIFS 2
II DESCRIPTION DE LA ZONE DETUDE 3
II.1 LOCALISATION GEOGRAPHIQUE 3
II.2 CONTEXTE CLIMATIQUE ET PLUVIOMETRIQUE 3
II.3 CONTEXTE HYDROLOGIQUE 5
II.4 CONTEXTE GEOLOGIQUE 6
II.5 CONTEXTE HYDROGEOLOGIQUE 10
III MODELISATION DES NAPPES DE LA PLAINE DES GALETS 15
III.1 DEFINITION DU MODELE CONCEPTUEL 15
III.2 MODELISATION DES ECOULEMENTS EN REGIME PERMANENT 21
III.3 MODELISATION DES INTRUSIONS SALINES 25
IV CONCLUSIONS 25
IV.1 SYNTHESE 25
IV.2 PROBLEMES RENCONTRES 25
IV.3 PISTES DAMELIORATION DU MODELE 26

ANNEXES 27
ANNEXE 1 : ZONES DE GEOMETRIE 28
ANNEXE 2 : ZONES D'ALIMENTATION 37









Liste des figures

Figure 1 : Localisation de la zone dtude _____________________________________ 3
Figure 2 : Rpartition de la pluviomtrie de lanne 2006 ________________________ 4
Figure 3 : Prcipitations journalires de la station Dos dAne (2002 2010)__________ 4
Figure 4 : Prcipitations journalires de la station du Port (2002 2010) ____________ 5
Figure 5 : Carte de situation des biefs dinfiltration de la Rivire des Galets __________ 6
Figure 6 : Chronique de dbit reconstitue de la Rivire des Galets (1988 2006) __ 6
Figure 7 : Carte gologique 1/50000me BRGM________________________________ 7
Figure 8 : Coupe gologique interprtative A-B ________________________________ 7
Figure 9 : Photographie arienne, campagne IGN 1966 __________________________ 8
Figure 10 : Synthse des donnes gophysiques disponibles sur la zone du Port ______ 9
Figure 11 : Schma de principe des nappes de la Plaine des Galets (ANTEA) ________ 10
Figure 12 : Coupe NW-SE passant proche de la centrale EDF (BRGM 1998) _________ 11
Figure 13 : Coupe NNE-SSW passant vers le pont de la RN1 (BRGM 1998) __________ 11
Figure 14 : Extension des nappes identifies sur la Plaine des Galets ______________ 14
Figure 15 : Carte pizomtrique gnrale de la Plaine des Galets (1992 1996) _____ 14
Figure 16 : Extension du polygone de modlisation sur fond gologique ____________ 16
Figure 17 : Reprsentation du maillage et des diffrents types de limites du modle __ 17
Figure 18 : Schma conceptuel du modle (ANTEA 2010)______________________ 18
Figure 19 : Niveaux pizomtriques simuls dans la nappe suprieure _____________ 22
Figure 20 : Niveaux pizomtriques simuls dans la nappe moyenne ______________ 23
Figure 21 : Niveaux pizomtriques simuls dans la nappe infrieure ______________ 24



Liste des tableaux


Tableau 1 : Permabilits des diffrentes formations gologiques _________________ 19
Tableau 2 : Porosits efficaces des diffrentes formations gologiques _____________ 20
Tableau 3 : Alimentation souterraine en limite Est _____________________________ 20
Tableau 4 : Alimentation par infiltration le long de la Rivire des Galets ____________ 20
Tableau 5 : Forages dexploitation et dbits moyen annuel exploit (20002008) ____ 21
Tableau 6 : Rsultats du calage pour la nappe suprieure _______________________ 22
Tableau 7 : Rsultats du calage pour la nappe moyenne ________________________ 24
Tableau 8 : Rsultats du calage pour la nappe infrieure ________________________ 25



1

Prambule



Une modlisation mathmatique des aquifres de la plaine des Galets a t
ralise en 2010 la demande du Conseil Gnral de La Runion par le bureau d'tude
ANTEA. Cette modlisation a tenu compte des phnomnes de salinisation marine sur le
secteur d'tude.

L'objet de ce travail est de reprendre le modle hydrodynamique de 2010 avec le
logiciel FEFLOW afin de pouvoir le faire voluer en fonction du temps, des nouvelles
connaissances et des nouveaux enjeux territoriaux.

Pour cela, le Conseil Gnral a transmis l'ensemble des donnes qui ont permis la
ralisation du modle mathmatique. Ces donnes ont t utilises pour concevoir le
modle avec l'outil de modlisation disponible l'Office de l'eau Runion. Cette premire
tape reprend strictement les choix de conception du bureau dtude que ce soit au
niveau :

des donnes sur la gomtrie des formations et le modle conceptuel

des valeurs des paramtres hydrodynamiques

des donnes dalimentation et dexploitation du systme







2
I Contexte et objectifs

Lle de La Runion est globalement caractrise par labondance de ses
ressources en eau grce une trs forte pluviomtrie. Ce constat doit nanmoins tre
nuanc lchelle des microrgions du fait dune grande variabilit spatiale et temporelle
des prcipitations.

A titre dexemple, la cte au vent est caractrise par des records mondiaux de
pluviomtrie alors que la cte sous le vent et prcisment la commune du Port est une
zone considre comme sche avec moins de 1000 mm de prcipitations par an
concentrs sur quelques vnements pluviomtriques intenses.

En vue de favoriser le dveloppement agricole, le Dpartement de La Runion a
lanc un projet de transfert des eaux dest en ouest. La solution retenue consiste
raliser des captages dans les rivires du cirque de Mafate (Bras Sainte Suzanne et
Rivire des Galets) et du cirque de Salazie (Rivire Fleurs Jaunes et Rivire du Mt) puis
de les relier par un rseau de 30 Km de galeries souterraines afin de distribuer leau
rcupre dans l'Ouest sur huit antennes d'irrigation rparties entre les communes du
Port et de Saint-Leu.

La corrlation dune gestion durable de la ressource en eau en lien avec les enjeux
du bon tat des masses deau en 2015 (Directive cadre sur leau) milite pour une
amlioration des connaissances du fonctionnement des aquifres souterrains de la Plaine
des Galets. Le SDAGE 2010-2015 a ainsi retenu plusieurs dispositions afin dassurer
lquilibre ressources/besoins pour les diffrents usages en prservant le milieu
aquatique (dispositions 1.1.3, 1.1.4 et 1.5.1).

La modlisation mathmatique reprsente un trs bon outil daide la dcision.
Elle permet d'apprhender les mcanismes, le fonctionnement et les relations entre
diffrentes zones aquifres. Plusieurs modles ont t raliss sur le secteur de la Plaine
des Galets successivement par ARMINES (1995) et lObservatoire Runionnais de l'Eau
(1999). Cependant, le basculement des eaux de lest vers louest, oprationnel dans
toute sa dimension au cours de lanne 2014 (parties Salazie et Mafate), est une
situation substantielle non reprsente par ces anciens modles.

Dans ce contexte, lOffice de l'eau Runion reprend la modlisation
hydrogologique ralise par ANTEA en 2010 sur ce secteur. Lutilisation du logiciel
FEFLOW-FM3 doit permettre d'valuer :
lquilibre global de la masse deau souterraine entre les prlvements et la
recharge naturelle,
les impacts potentiels quantitatifs et qualitatifs,
les implications sur les phnomnes d'intrusion salines au sein du systme
aquifre.

Dans un premier temps, cette note prsentera le contexte hydrogologique et
pluviomtrique de la zone dtude. Dans un deuxime temps, une description du modle
conceptuel dtaillera les choix de construction (gomtrie, conditions aux limites) et les
diffrents paramtres hydrodynamiques introduits (permabilit, porosit). Dans un
troisime temps, les rsultats de la modlisation seront prsents ainsi que les
perspectives d'volution court terme de l'outil.






3
II Description de la zone dtude
II.1 Localisation gographique

La zone dtude correspond aux aquifres associs la Rivire des Galets. La
Plaine des Galets se situe au nord-ouest de lle de La Runion sur les communes du Port
et de Saint-Paul. Ses limites correspondent sensiblement la bordure de la planze
Sainte Thrse et au niveau des dbouchs des diffrentes ravines lest, la Rivire des
Galets au sud et lOcan Indien au nord et louest.


Figure 1 : Localisation de la zone dtude

II.2 Contexte climatique et pluviomtrique

Situ proximit immdiate du tropique du Capricorne dans le sud-ouest de
lOcan Indien lle de La Runion est caractrise par un climat tropical. Cependant, plus
de 300 microclimats sont recenss sur lle. La commune du Port, situe sur la cte sous
le vent, dtient les plus faibles cumuls de prcipitations annuelles avec le record de
temprature enregistr 40C.

Figure 2 prsente les cumuls pluviomtriques enregistrs en 2006. Malgr le
passage de deux temptes tropicales, la carte illustre bien l'ingale rpartition observe
entre lest et louest.

4



Figure 2 : Rpartition de la pluviomtrie de lanne 2006

Lanalyse des prcipitations au niveau de la Plaine des Galets porte sur deux
stations du rseau Mto France : Le Port (altitude : 11 mNGR) et Dos dAne (altitude :
930 mNGR). En moyenne, les pluies enregistres la station de Dos dAne sont deux fois
plus importantes que celles de la station du Port, consquence logique de lexistence dun
gradient altimtrique des prcipitations la Runion.

0
100
200
300
400
500
600
s
e
p
t
.
-
0
2
o
c
t
.
-
0
3
n
o
v
.
-
0
4
d

c
.
-
0
5
j
a
n
v
.
-
0
7
f

v
r
.
-
0
8
m
a
r
s
-
0
9
m
a
i
-
1
0
j
u
i
n
-
1
1
(
m
m
)

Figure 3 : Chronique de prcipitations journalires de la station Dos dAne
(2002 2010)
ZONE
D'ETUDE
Source : METEO FRANCE

5

0
20
40
60
80
100
120
140
s
e
p
t
.
-
0
2
j
a
n
v
.
-
0
4
m
a
i
-
0
5
o
c
t
.
-
0
6
f

v
r
.
-
0
8
j
u
i
l
.
-
0
9
n
o
v
.
-
1
0
(
m
m
)

Figure 4 : Chronique de prcipitations journalires de la station du Port
(2002 2010)

II.3 Contexte hydrologique

La Rivire des Galets est lexutoire principal du cirque de Mafate. Son bassin
versant a une superficie de 117 km son embouchure.

A partir de 1989, une station denregistrement du dbit install au niveau du Cap
Lebot, laval de la confluence avec le Bras Sainte-Suzanne, a permis de suivre en
continu lvolution du dbit de la Rivire des Galets. A partir de 2002, cette station est
remplace par celle de lIlet Malid, lgrement plus laval, qui a cess de fonctionner
en fvrier 2006 suite aux consquences de la tempte tropicale Diwa . Entre les deux
stations une lgre infiltration de lordre de quelques dizaines de litres par seconde est
observe.

Le dbit de la Rivire des Galets augmente de lamont du bassin versant jusquau
Cap Lebot. Ensuite, le dbit de la rivire diminue la faveur de trois zones dinfiltration
(cf. Figure 5) :

Le bief amont : il stend du Cap Lebot la prise du canal Lemarchand. Le
dbit maximum susceptible de sy infiltrer est de 0.07 m
3
/s,

Le bief intermdiaire : il stend de la prise du canal Lemarchand au pont
de la RN1. Le dbit maximum susceptible de sy infiltrer est de 0.47 m
3
/s,

Le bief aval : il stend du pont de la RN1 locan Indien. Le dbit
maximum susceptible de sy infiltrer est de 0.45 m
3
/s.





6

Figure 5 : Carte de situation des biefs dinfiltration de la Rivire des
Galets

0
20000
40000
60000
80000
100000
120000
140000
160000
d

c
e
m
b
r
e
-
8
7
d

c
e
m
b
r
e
-
8
8
d

c
e
m
b
r
e
-
8
9
d

c
e
m
b
r
e
-
9
0
d

c
e
m
b
r
e
-
9
1
d

c
e
m
b
r
e
-
9
2
d

c
e
m
b
r
e
-
9
3
d

c
e
m
b
r
e
-
9
4
d

c
e
m
b
r
e
-
9
5
d

c
e
m
b
r
e
-
9
6
d

c
e
m
b
r
e
-
9
7
d

c
e
m
b
r
e
-
9
8
d

c
e
m
b
r
e
-
9
9
d

c
e
m
b
r
e
-
0
0
d

c
e
m
b
r
e
-
0
1
d

c
e
m
b
r
e
-
0
2
d

c
e
m
b
r
e
-
0
3
d

c
e
m
b
r
e
-
0
4
d

c
e
m
b
r
e
-
0
5
d

c
e
m
b
r
e
-
0
6
d

c
e
m
b
r
e
-
0
7
D

b
i
t

(
L
/
s
)

Figure 6 : Chronique de dbit reconstitue de la Rivire des Galets
(1988 2006)

La Figure 6 montre lamplitude des variations de dbit de la Rivire des Galets lors
dvnements pluvieux importants tels que les cyclones et temptes tropicales qui
touchent lle de la Runion. Lintensit du dbit peut monter jusqu' plus de 100 m3/s.
II.4 Contexte gologique
La Plaine des Galets rsulte de la mise en place dun cne de djection sur les
formations volcaniques du Piton des Neiges. La nature et la configuration des dpts
traduisent une capacit de transport de matriel sdimentaire exceptionnelle lie aux
priodes de crues. La granulomtrie est trs varie, des silts aux blocs plurimtriques.
Des niveaux discontinus plus ou moins limoneux peuvent sparer des alluvions plus
propres. Au sein des alluvions, les contrastes de permabilits crent diffrents niveaux
pizomtriques et impliquent la mise en jeu de phnomnes de drainance entre les
diffrentes units hydrogologiques.

Ilet Malid Cap Lebot

7
La Plaine des Galets se raccorde vers le sud-est aux planzes de Sainte-Thrse
et de la Plaine Saint-Paul qui forment les pentes du volcan bouclier. Le massif volcanique
est constitu par des coules anciennes des phases I et II du massif du Piton des Neiges.
Elles sont recouvertes au sud de la Ravine Marquet par des formations plus rcentes
des phases III et IV.

Une ancienne coule boueuse (lahar) s'est mise en place dans le lit de la Rivire
des Galets avant les coules de la phase III. Cette formation indure a une trs faible
permabilit. Elle contraint les coulements souterrains aliments par les infiltrations de
la rivire. Les flux se dirigent en grande partie vers la commune du Port sous l'effet de
cette formation.


Figure 7 : Carte gologique 1/50000me BRGM


Figure 8 : Coupe gologique interprtative A-B
A
B

8
La Figure 8 prsente une coupe selon l'axe AB (cf.Figure 7). Elle montre que les
coules basaltiques se prolongent sous le niveau de la mer et que les alluvions forment
une plateforme deltaique.

La position du lit majeur de la Rivire des Galets a varie au cours du temps. Sur
la photographie arienne (Figure 9), les palo-chenaux balayent toute la plaine alluviale.
Certains dentre eux sont mmes de direction sud-nord alors que les plus rcents sont de
direction est-ouest. On peut donc imaginer quil existe des palo-valles remplies
dalluvions ailleurs qu laplomb du lit actuel de la rivire.


Figure 9 : Photographie arienne, campagne IGN 1966

Cette hypothse a t vrifie par des campagnes de prospections gophysiques.
Une campagne sismique ralise en 1973 (Rapport BRGM n 73RME029AF) a mis en
vidence diffrentes palo-valles. La Figure 10 synthtise l'ensemble des informations
disponibles.

ANCIENS CHENAUX D'ECOULEMENTS
ORIENTES VERS LE NORD
LIT ACTUEL DE LA RIVIERE DES GALETS
9

Figure 10 : Synthse des donnes gophysiques disponibles sur la zone du Port

10
II.5 Contexte hydrogologique

Le schma hydrogologique de synthse bas sur l'analyse de nombreuses
donnes d'interprtations de travaux de forages et de suivis pizomtriques propose trois
nappes superposes et dlimites par des horizons impermables ou semi permables :

une nappe suprieure dans les alluvions grossires, dont les limites
sloignent trs peu du lit de la Rivire des Galets. Elle est libre, jusqu
65 m/sol de profondeur, son niveau oscille autour de +40 m NGR. Son
alimentation est essentiellement lie aux infiltrations des eaux de la Rivire
des Galets. Les flux d'infiltrations sont essentiellement verticaux ce qui
positionne la rivire bien au dessus du niveau de la nappe.

des alluvions limoneuses : aquiclude.

une nappe moyenne dans des alluvions grossires et anciennes. Elle est
gnralement sous-jacente la nappe suprieure, spare par un niveau
semi-permable. La nappe moyenne est semi captive entre 80 et 120 m de
profondeur, son niveau se situe entre +2 et +5 m NGR. Son alimentation
est semblable la nappe suprieure dont elle draine une partie des eaux.

des alluvions limoneuses : aquiclude

une nappe infrieure dans un mlange complexe de formations volcaniques
et dalluvions. Elle est captive partir denviron 125 m de profondeur et
son niveau se situe entre +2 et +2,5 m NGR. Son alimentation est
principalement reprsente par des coulements souterrains horizontaux
correspondant aux eaux infiltres travers les formations basaltiques de la
planze Sainte-Thrse, la Ravine Marquet et la Ravine des Lataniers. Un
complment par soutirage de la nappe moyenne est possible la faveur de
discontinuits des formations impermables qui les sparent.


Figure 11 : Schma de principe des nappes de la Plaine des Galets
(ANTEA)
Les Figure 12 et Figure 13 illustrent la complexit de cette plaine alluviale.
Lextension et lpaisseur des couches permables et des couches impermables peuvent
tre trs varies dans lespace.

11


Figure 12 : Coupe NW-SE passant proche de la centrale EDF (BRGM 1998)


Figure 13 : Coupe NNE-SSW passant vers le pont de la RN1 (BRGM 1998)
II.5.1 Secteur de la Rivire des Galets

Nappe suprieure

Elle est alimente par les infiltrations provenant de la Rivire des Galets qui
peuvent atteindre jusqu 600 l/s. Le gradient est denviron 3% lamont et infrieur
1% laval de la RN1.

Lanalyse des comportements pizomtriques vise rduire lextension de cette
nappe proximit de la rivire. La remonte de la coule boueuse proximit du centre
pnitencier induit un effet de seuil.

Les eaux de la nappe prennent une direction nord. Les volutions pizomtriques
sont comparables sur lensemble des ouvrages captant cette nappe. Lpaisseur de cette
nappe est dune dizaine de mtre.



12
Nappe moyenne

La nappe moyenne est alimente principalement par drainance de la nappe
suprieure et/ou directement partir des infiltrations de la rivire plus lamont. Une
autre source dalimentation proviendrait de limpact direct dpisodes de fortes
prcipitations.

On observe bien ce phnomne au niveau de lvolution pizomtrique.
Effectivement, les chroniques pizomtriques montrent une remonte rapide (quelques
jours) en dbut danne en cas dvnements pluvieux importants pouvant atteindre une
vingtaine de centimtres suivi dune seconde augmentation, plus lente et plus importante
(20 60 cm).

Les tudes du BRGM prennent en compte une nappe moyenne individualise dans
laxe de la Rivire des Galets et des nappes moyenne et suprieure confondues louest
et au nord. Lcoulement des eaux souterraines de cette nappe fortement exploite se
fait essentiellement vers le nord. Comme pour la nappe suprieure, cest probablement la
coule boueuse qui dtourne lcoulement de cette unit. Ses variations
pizomtriques sont retardes par rapport la nappe suprieure, ce qui illustre bien
leffet de drainance verticale entre les deux systmes.

Nappe infrieure

Une troisime nappe a t identifie sous les deux prcdentes. Ses eaux sont
saumtres, sa pizomtrie avoisine les 2,50 m NGR et ses variations annuelles sont trs
faibles. Elle prsente un comportement typique dune nappe alimente par les
infiltrations par les planzes lamont hydraulique.

Son maximum pizomtrique se situe juste aprs les fortes pluies de la saison
cyclonique. Elle semble en continuit hydraulique sur lensemble du secteur, dabord
dans les formations volcaniques de la planze, ensuite dans les alluvions. Les mesures de
conductivit via les ouvrages profonds montrent quelle est quasiment totalement
envahie par le sel.

II.5.2 La nappe libre de la Planze Sainte Thrse
Cette nappe existe la faveur des infiltrations dans les basaltes. Elle contribue
l'alimentation du systme de la Plaine des Galets. La pizomtrie est faible et avoisine les
2,5 m NGR sur lensemble du secteur. Elle prsente le comportement typique dune
nappe alimente de manire quasi continue par les infiltrations sur la planze avec un
soutien de ltiage en milieu danne et des alimentations rapides lors de la saison des
pluies. Des apports par les dbits de la Ravine Marquet sont visibles sur les
enregistrements pizomtriques.

L'amplitude de ses variations ne dpasse pas 50 cm. Ses caractristiques
gologiques sont identiques la nappe infrieure voque prcdemment mais les eaux
de cet aquifre sont douces. Du point de vue gologique, les donnes de forages sur ce
secteur rvlent qu'une trentaine de mtres de coules basaltiques de Phase IV reposent
sur des coules de Phase III ou directement sur des coules de Phase II.
II.5.3 Les nappes du Port Est

Dans la zone allant du centre pnitencier jusquau Port Est, plusieurs nappes
superposes ont t identifies sur les ouvrages FT6, P11, FT5 et P11. Linvasion par des
eaux dorigine marine est trs marque sur la nappe infrieure du FT6 ainsi que sur le
mlange des deux nappes au FT5 et proximit de la cote.


13
II.5.4 Le Secteur aval de la Ravine Marquet

Ce secteur a fait lobjet de plusieurs tudes menes par ANTEA dans le cadre de
lextension du Port Est. Une dizaine de pizomtres ont t implants pour amliorer les
connaissances gologiques sur cette zone. D'une profondeur moyenne de 50 mtres
(environ -40 m NGR), ils ont rencontr exclusivement des formations alluvionnaires
allant du sable limoneux aux graves puis aux blocs.

Laquifre rencontr au droit du site est en relation directe avec lOcan Indien. Il
est localis dans des formations fluvio-marines, alimentes par les basaltes directement
lest (amont hydraulique) et par les infiltrations directes de la Ravine Marquet.

La prsence dune coule boueuse ou de niveaux alluvionnaires fins plus ou moins
impermables implique galement lexistence de niveaux pizomtriques localiss et
perchs (Rapport ANTEA 36061/A novembre 2004). La Ravine Marquet semble
constituer une discontinuit du systme aquifre par la prsence suppose dune coule
boueuse qui oriente les coulements deau douce au dessus du niveau aquifre gnral.

On retiendra de cette zone :
Des niveaux de base observs entre 0,5 et 1 m NGR lest du secteur,
La prsence de niveaux perchs avec des cotes de lordre de 3,5
4 m NGR sur la partie proche du port Est,
Des amplitudes pizomtriques importantes lies la prsence de la
Ravine Marquet en cas dpisodes de fortes pluies (2,5 3,5 m NGR sur
la partie Est du secteur, 6 m NGR sur la partie ouest).
II.5.5 Les nappes aval de la Rivire des Galets
La partie aval de la Rivire des Galets prsente dans laxe de la rivire, des
niveaux saturs discontinus qui se rpartissent de part et d'autre de la rivire. Ces
niveaux discontinus existent la faveur des infiltrations du cours deau. Les variations y
sont importantes avec une diffrence marque entre la rive droite et la rive gauche.

II.5.6 Conclusions
Le schma conceptuel du fonctionnement hydrogologique de la Plaine des Galets
peut se dfinir comme :
La superposition de trois nappes dans laxe de la rivire, plus ou moins
connectes entre elles par drainance verticale travers les formations qui les
sparent.
Un systme monocouche ( lexception des niveaux perchs) lorsque on se
rapproche du littoral. Ces nappes sont certainement le rsultat de la disparition
progressive des niveaux (cf. Figure 12 et Figure 13).
Une nappe libre sur la Planze Sainte-Thrse dans les formations volcaniques.
Elle se poursuit dans les alluvions dont la pizomtrie est fortement lie aux
prcipitations sur la planze.


14

Figure 14 : Extension des nappes identifies sur la Plaine alluviale de la
Rivire des Galets

La Figure 15 donne une bonne ide des tendances pizomtriques de chacune des
trois nappes de la Rivire des Galets ainsi que le sens dcoulement gnral.


Figure 15 : Carte pizomtrique gnrale de la Plaine des Galets
(synthse BRGM 1992 1996)


15
III Modlisation des nappes de la Plaine des Galets

Lobjectif de cette tude est dobtenir un outil permettant de simuler les
coulements souterrains du systme aquifre de la Plaine des Galets puis d'valuer la
morphologie de l'interface eau douce/eau sale.

Pour cela, il est ncessaire d'adapter le modle ralis par le bureau d'tude
l'outil de modlisation FEFLOW. L'enjeu est de pouvoir actualiser et faire voluer le
modle en fonction des nouvelles connaissances. Cet outil permettra d'valuer la
ressource potentiellement mobilisable et sa vulnrabilit l'chelle de l'unit aquifre
littorale de Saint-Paul/La Possession. Par contre, l'chelle d'observation et la prcision du
modle ne permettront pas d'estimer l'impact d'un forage en particulier.
III.1 Dfinition du modle conceptuel
III.1.1 Extension et construction du modle
Dlimitation, maillage et limites.

Lextension de la zone modliser est base sur des critres gologiques et sur la
connaissance des bassins versants hydrologiques du secteur. Les limites du modle sont
illustres sur la Figure 16.

La limite sud est positionne sur la limite gologique suppose impermable que
constitue la coule boueuse. Cette hypothse peut tre discute car on ne connat pas
prcisment la profondeur de cette coule (qui passe par endroit sur les alluvions) et il
n'est pas exclu qu'une partie de l'alimentation aille vers Cambaie au niveau du Grand
Pourpier.

Vers le Sud-Ouest, la coule disparat et la limite devient localement permable
vers la rive gauche. La limite Est reprsente la limite de reconnaissance par forage. Au
Nord et lOuest, les limites sont naturellement tablies au niveau du trait de cte. Par
ailleurs, le modle a t tendu vers le sud-est dans l'axe de la rivire au sud-est afin
d'englober la totalit de la zone dinfiltration souterraine.


16

Figure 16 : Extension du polygone de modlisation sur fond gologique

17


Figure 17 : Reprsentation du maillage et des diffrents types de limites
du modle
Le maillage en lments finis prsent sur la Figure 17 est affin le long de la
Rivire des Galets et lemplacement des ouvrages de rfrence connus sur le
secteur. Au total, sur toute lpaisseur du modle, on compte 228 390 lments et
124 848 nuds.

Couches et gomtrie du modle

Le nombre de couches a t fix 15 par le bureau d'tude. La gomtrie des
couches est base sur la diffrenciation lithologique des formations gologiques
prsentes au sein du cne alluvial de la Plaine des Galets. Les valeurs de conductivits
hydrauliques ont t attribues pour chaque couche dans des zones de gomtries en
fonction de la gologie rencontre. Ces zones de gomtries sont prsentes pour
chaque couche en annexe 1.

Les 15 couches correspondent aux formations gologiques suivantes de la surface
vers le bas :

les couches 1 3 reprsentent les alluvions de la nappe suprieure. Les
zones de gomtrie individualisent les lahars de la Rivire des Galets et de
la Ravine Marquet, les basaltes des phases III-IV et les alluvions de la
nappe suprieure.

Forages dexploitations
Limites potentiel impos
Limites flux impos

18
la couche 4 : aquiclude reprsentant la fois le mur de la nappe
suprieure et le toit de la nappe moyenne.

les couches 5 7 correspondent globalement la nappe moyenne. Elles
sont nanmoins subdivises en plusieurs zones de gomtries,
correspondant aux alluvions, au lahar et aux basaltes des phases III-IV.
Dans la partie basaltique des planzes, la couche 7 correspond au palosol
qui existe rgulirement entre les formations de phases II et III. Ce niveau
est trs peu permable et bloque ou ralentit les coulements qui pourraient
se faire vers le bas.

la couche 8 correspond laquiclude reprsentant le mur de la nappe
moyenne et le toit de la nappe infrieure. Elle recoupe le lahar au sud du
modle dans sa partie infrieure.

les couches 9 15 se situent entre la base de laquiclude sus-jacent et la
base du modle dfinit -400 m NGR. Les zones de gomtries
individualisent les alluvions de la nappe infrieure et les basaltes de phase
II.



Figure 18 : Schma conceptuel du modle
(ANTEA Rapport n 55794/D-anne 2010)

Points de contrles

Les points de contrle correspondent aux ouvrages (forages dexploitations,
pizomtres) pour lesquels des informations sur les charges hydrauliques sont

19
disponibles. Ils permettent de comparer les niveaux simuls aux niveaux mesurs sur le
terrain afin d'valuer la finesse et la justesse de la modlisation.

III.1.2 Paramtres hydrodynamiques du modle et dbits prlevs

Les paramtres hydrodynamiques introduits dans le modle correspondent aux
paramtres directement lis la nature des formations gologiques composant la Plaine
des Galets mais aussi aux valeurs dalimentation et dexploitation du systme.

Ils sont dfinis par :
les couches et les zones de gomtries correspondantes,
les valeurs de flux entrant par les bassins dalimentation en amont
hydraulique du modle,
les infiltrations de la Rivire des Galets,
Les valeurs de flux sortant par exploitation.

Illustrs par la Figure 18, les flux dalimentations en limite est du modle sont
distribus sur les couches 9 15 et les flux dinfiltrations de la Rivires des Galets sont
distribus sur la couche 1. Des limites potentiels imposs gals 0 sont assignes sur
toutes les couches au niveau du littoral.

Les permabilits.

Les valeurs issues de linterprtation de pompages dessais et celles dtermines
pendant le calage du modle (ANTEA, Rapport n55794/D) restent dans les mmes
gammes de valeurs. Elles seront utilises comme valeurs de rfrences.


Tableau 1 : Permabilits des diffrentes formations gologiques
(ANTEA Rapport n 55794)








20
Les porosits.

La porosit efficace gouverne lcoulement de leau au sein dun massif. Les
valeurs introduites dans le modle apparaissent dans le tableau 2 ci-dessous.

Porosits efficaces
Aquifres 15%
Aquicludes 25%
Lahars 5%
Basaltes 5%
Palosol 30%

Tableau 2 : Porosits efficaces des diffrentes formations gologiques
(ANTEA Rapport n 55794)


Lalimentation du systme aquifre.

Trois sources principales dalimentation des nappes de la Plaine des Galets ont t
dfinies : une dorigine souterraine provenant des coulements au sein du massif
basaltique, une autre due aux infiltrations le long de la Rivire des Galets et une lie aux
apports par les prcipitations (130 mm/an). La rpartition de ces flux sur le modle est
prsente en annexe 2.


Alimentation souterraine travers le massif basaltique

Ces flux souterrains sont la consquence des infiltrations deau depuis la Ravine
des Lataniers, la Ravine Marquet et la Planze Sainte-Thrse. Ils sont appliqus en
limite Est du modle.

Dnomination Volume annuel (m3) Dbit / nuds (m3/s)
R. des Lataniers 4 953 675 9.03516 E-7
R. Marquet 3 991 196 2.43089 E-7
Planze Sainte-Thrse 5 173 544 1.181065 E-7

Tableau 3 : Alimentation souterraine en limite Est


Infiltration le long de la Rivire des Galets

On diffrencie trois biefs dinfiltrations diffrents, amont, intermdiaire et aval (cf.
Figure 5). Les volumes correspondent au calcul de la distribution des trois modules
dinfiltration le long de la rivire. Ces valeurs sont appliques sur la couche de surface du
modle et les volumes se dversent donc dans la nappe suprieure.

Bief dinfiltration Volume annuel (m3) Dbits / noeuds (m3/s)
Amont 2 516 573 1.13 E-6
Intermdiaire 16 022 811 2.9245 E-6
Aval 11 324 578 2.87 E-6

Tableau 4 : Alimentation par infiltration le long de la Rivire des Galets






21
Lexploitation du systme aquifre

Lexploitation du systme se fait la faveur de 11 ouvrages destins
lalimentation en eau potable de la population. La plupart dentre eux exploitent la nappe
moyenne, alors que dautres captent les eaux du massif basaltique plus lEst.

Les valeurs fournies dans le tableau 5 sont des moyennes des modules annuels
sur la priode 2000-2008 (rapport ANTEA n 55794/D).

Nom usuel code BSS Code OE Dbits (m3/s) Couche atteinte
Forage Balthazar 12263X0087/F 13087 0.0262 9
Forage Ravine Marquet 12263X0111/S2 13001-B 0.0495 9
Forage SAMY 12263X0050/P 13050 0.0629 9
Forage F6 12262X0175/F6 12175 0.0249 9
Forage Mounien 12262X0044/F1 0.0153 5
Forage F2 12262X0135/F2 12135 0.0114 5
Forage F3 12262X0136/F3 12136 0.0463 5
Forage F4 12262X0048/F4A 12048 0.0296 5
Forage F5 12262X0168/F5 12168 0.0268 5
Forage P11 12262X0070/P1 12070 0.0537 5
Forage EDF 12262X0046/F 12046 0.0347 5

Tableau 5 : Forages dexploitation et dbits moyen annuel exploit
(2000 2008)
III.2 Modlisation des coulements souterrains en rgime permanent
III.2.1 Mthode de calage du modle en rgime permanent
Cette tape consiste ajuster le schma conceptuel en fonction des rsultats des
niveaux pizomtriques obtenus par les premires modlisations. Des rsultats de calage
satisfaisants ont t obtenus aprs plus de 20 essais. La principale difficult a rsid
dans l'adaptation des donnes d'entres aux mthodes de calculs diffrentes entre les
deux modles.

Dans l'ensemble, les valeurs de permabilits initiales n'ont pas t modifies
l'exception de celles de la nappe moyenne qui a du tre diminue d'un facteur 10
(5,0.10
-5
m/s).
III.2.2 Rsultats de calage

Les Figure 19, Figure 20 et Figure 21 prsentent la pizomtrie obtenue pour la
nappe suprieure, la nappe moyenne et la nappe infrieure. Les points de contrles
(pizomtres et forages dexploitations) sont aussi reprsents.















22
Nappe suprieure


Figure 19 : Niveaux pizomtriques simuls dans la nappe suprieure
(en m NGR)

Les niveaux simuls sur la nappe alluviale suprieure sont en corrlation avec les
niveaux mesurs sur les points de rfrence. Lalimentation due aux infiltrations de la
Rivire des Gales et leffet de seuil cr par le lahar en limite sud sont bien visible. Cette
formation trs impermable a pour effet de dvier les coulements vers le nord. Certains
nuds prsentent des pizomtries abrrantes suprieures 150 mNGR. Ils se trouvent
proximit du lit de la Rivire des Galets o le modle doit combiner des flux importants
venant des infiltrations et la barrire hydraulique impermable qui bloque les
coulements (lahar). Ces artfacts numriques restent localiss et ils n'ont pas
d'influence relle sur les rsultats du modle.

Code BSS Code OE Mesure Simule
12263X0077/F5 13077 0.28 0.69
12262X0056/P3 12056-A 2.2 3.9
12262X0172/P1-2A 12172 41.73 46.04
12262X0057/P4-A 12057 67.38 69.77
12262X0065/P8-A 12065 27.22 26.93
12262X0052/SPZ1 12052 17 11.62
Pizomtre Pizomtrie (mNGR)
Nappe alluviale suprieure

Tableau 6 : Rsultats du calage pour la nappe alluviale suprieure

23

Comme le montre le tableau 6, les ordres de grandeur des niveaux
pizomtriques simuls sont en accord avec les niveaux rels. La prcision n'est pas trs
bonne avec des niveaux simuls suprieurs la ralit mais les tendances et les sens
d'coulement sont globalement bien reprsents.

Nappe moyenne


Figure 20 : Niveaux pizomtriques simuls dans la nappe moyenne
(en m NGR)
Au niveau de la nappe alluviale moyenne de la Rivire des Galets, la corrlation
est bonne entre des niveaux pizomtriques simuls et les niveaux mesurs sur les
ouvrages de rfrence. Seuls les niveaux pizomtriques simuls pour les forages 12171-
B et 12123 sont en dclage avec les niveaux observs.

Le secteur Sud-Est est marqu par une alimentation importante associe au bief
dinfiltration amont de la Rivire des Galets. Le modle reprsente clairement l'influence
des puits de pompage EDF et P11 (cercles bleus).



24
Code BSS Code OE Mesure Simule
12262X0134/GTM-2 12134 1.8 1.66
12262X0066/FT-6A 12066 2.97 2.24
12262X0064/FT-5A 12064 1.76 1.12
12262X0056/P3 12056-B 2.47 2.77
12262X0337/P2-B 12173-B 5.35 4.98
12262X0336/P1-1B 12171-B 3.22 0.29
12262X0123/SEP-2 12123 0.57 2.59
Nappe alluviale moyenne
Pizomtre Pizomtrie (mNGR)

Tableau 7 : Rsultats du calage pour la nappe alluviale moyenne
Nappe infrieure


Figure 21 : Niveaux pizomtriques simuls dans la Nappe infrieure
(en m NGR)

La simulation des niveaux pizomtriques de la nappe alluviale infrieure de la
Rivire des Galets est quand elle moins bien corrle avec les niveaux mesurs sur les
ouvrages de rfrences. En effet, les niveaux simuls dans les pizomtres de rfrence
sont plus bas que les niveaux mesurs (environ un mtre d'cart).

Il est trs probable que cet cart soit li la prsence de sels dans la nappe
infrieure. A ce stade du calage, l'apport de sels n'a pas t modlis. L'ajout de sels

25
dissous dans l'eau augmentera la densit du fluide ce qui entranera une hausse des
niveaux pizomtriques dans les ouvrages de rfrence.

De plus lalimentation en limite Est, cense alimenter prioritairement cette nappe
nest pas visible et donne limpression dune limite impermable. Cependant, la tendance
dune nappe avec un faible gradient et des niveaux pizomtriques bas (2 2.5 mNGR)
compar aux nappes qui lui sont suprieure est respecte.

Code BSS Code OE Mesure Simule
12262X0318/FT-6B 12066-B 2.7 1.1
12262X0174/P3 12174 2.37 1.3
12263X0036/S4 13036 2.59 0.75
12262X0053/S2 12053 2.8 1.1
Pizomtre Pizomtrie (mNGR)
Nappe alluviale infrieure

Tableau 8 : Rsultats du calage pour la nappe alluviale infrieure

III.3 Modlisation des intrusions salines

Cette tape n'a pas pu tre modlise car la structure complexe de la gomtrie
du modle entrane de nombreuses erreurs numriques qui bloquent le modle pendant
la phase de calcul.

La mthode ncessite des ajustements de la gomtrie du modle conceptuel afin
dobtenir des rsultats conformes la ralit. De nouveaux essais seront raliss lors
dune prochaine tude aprs une simplification de la structure du modle et une
extension de ses limites la rive gauche de la Plaine des Galets.

IV Conclusions
IV.1 Synthse
La conception du modle mathmatique reprsentatif de lcoulement souterrain
permet dapporter une aide la comprhension du systme aquifre complexe de la
plaine alluviale de la Rivire des Galets.

Les rsultats obtenus sont suffisamment reprsentatifs des diffrences de
comportement reconnues entre les trois units principales de ce systme. Les tendances
pizomtriques caractristiques des nappes suprieure, moyenne et infrieure indiquent
une image en corrlation avec la ralit. Lquilibre rel des circulations deau
souterraines entre la recharge, les pertes en mer et les prlvements du systme
aquifre est respect.

Cet outil informatique permet de simuler limpact de certains scnarii comme par
exemple une augmentation des prlvements ou une diminution de la recharge et ainsi
tre capable dalerter le public des consquences encourues.

IV.2 Problmes rencontrs
La reprise du modle hydrodynamique ralis par ANTEA, sous le logiciel FEFLOW
a t difficile.

Le logiciel de modlisation MARTHE utilis par ANTEA travaille en diffrences finies
alors que le logiciel FEFLOW fonctionne en lments finis. Cette diffrence demande un
certain nombre de rajustements des donnes dentre.

26

Dans un premier temps, une restructuration complte de ces fichiers dentre a
t ncessaire. Les fichiers concernant la topographie des diffrentes couches du modle
ont d tre arrangs de faon obtenir la bonne structure gomtrique. Ensuite, les
fichiers concernant les paramtres hydrodynamiques ont d tre, eux aussi, retravaills
afin de les introduire dans le logiciel de modlisation.

Dans un second temps, les flux dalimentation ont d tre recalculs pour les
assigner un certain nombre de nuds contrairement au logiciel MARTHE qui rparti les
flux sur des surfaces reprsentes par les mailles.

Toutes ces tapes compliquent la reprise dun modle partir des donnes
utilises par un autre logiciel de modlisation, qui plus est, lorsque lun travaille en
diffrences finies et lautre en lments finis.

Dautres difficults sont apparues lobtention des rsultats de modlisation. La
pizomtrie de certains nuds est incohrente avec lensemble des rsultats. Des
pizomtries extrmement leves se retrouvent des endroits prcis du modle et sur
des couches diffrentes. Lorigine de ces problmes provient de la dfinition du modle
conceptuel et plus particulirement du choix des limites. Les contrastes de faibles et
fortes permabilits aux contacts entre basaltes fissurs (0.02 m/s) et lahars (1,0E
-9
m/s) sont trop proches et provoquent des erreurs de calcul.

IV.3 Pistes damlioration du modle

Il parait essentiel dobtenir de meilleurs rsultats en rgime permanent,
notamment pour la nappe alluviale infrieure, avant de pousser l'analyse en rgime
transitoire et dtudier les intrusions salines dans ce systme aquifre. Ainsi, l'influence
de la prsence de sels dans la nappe infrieure pourra tre observe sur les niveaux
pizomtriques. .

Afin damliorer les rsultats du modle en rgime permanent, plusieurs pistes
sont envisages :
Extension de la zone de modlisation au niveau de la rive gauche de la
Rivire des Galets. Cette extension aura pour but dviter de mettre une
limite flux impos proximit dune limite flux nul ;
Rduction du nombre de couches pour faciliter le calcul ;
Attnuation des contacts brutaux entre des coules boueuses trs peu
permables et les alluvions plus permables.





27


















ANNEXES


28


















Annexe 1
Zones de gomtries

























29

Extension des diffrentes zones de gomtries Couche 1





Extension des diffrentes zones de gomtries Couche 2

30

Extension des diffrentes zones de gomtries Couche 3





Extension des diffrentes zones de gomtries Couche 4

31

Extension des diffrentes zones de gomtries Couche 5





Extension des diffrentes zones de gomtries Couche 6


32

Extension des diffrentes zones de gomtries Couche 7





Extension des diffrentes zones de gomtries Couche 8


33

Extension des diffrentes zones de gomtries Couche 9





Extension des diffrentes zones de gomtries Couche 10



34

Extension des diffrentes zones de gomtries Couche 11





Extension des diffrentes zones de gomtries Couche 12


35

Extension des diffrentes zones de gomtries Couche 13





Extension des diffrentes zones de gomtries Couche 14


36

Extension des diffrentes zones de gomtries Couche 15


































37

























Annexe 2
Zones dalimentation
























38

Extension des diffrentes zones infiltration au niveau de la Rivire des Galets






Extension des diffrentes zones infiltration au niveau de la limite Est du modle

39

Bibliographie

Rapports :


ANTEA (2009) Cne alluvial de la Rivire des Galets Elaboration dun modle
mathmatique dvaluation des impacts sur la ressource souterraine. Rapport de Phase I.
Juin 2009. Rapport n 54086/B.

ANTEA (2010) Cne alluvial de la Rivire des Galets Elaboration dun modle
mathmatique dvaluation des impacts sur la ressource souterraine. Rapport de Phase
II : Calage. Juillet 2010. Rapport n 55794/E.

ANTEA (2010) Cne alluvial de la Rivire des Galets Elaboration dun modle
mathmatique dvaluation des impacts sur la ressource souterraine. Rapport de Phase
III : Rsultats des simulations. Dcembre 2010. Rapport n60176/B.

Office de lEau Runion (1999) Modle Mathmatique de laquifre infrieur de la
Plaine des Galets. Etat des connaissances Analyse du modle

Office de lEau Runion (1996) Modlisation mathmatique du systme aquifre de
la Plaine des Galets. Rapport de stage de 3 anne. Juin 1996

Office de lEau Runion (2002) Redfinition dun modle de nappe multicouche pour
la prvision des niveaux pizomtriques dtiage des aquifres de la Plaine des Galets.
Modlisation sous GMS/Modflow en approche True Layer prliminaire au calage au
pas de temps mensuel de dcembre 1991 novembre 1997. Rapport de stage de 3
anne. Septembre 2002

ARMINES Etude des ressources en eau souterraine de la rgion de Saint-Paul. Ile de la
Runion. Construction et exploitation du modle mathmatique (texte).
LHM/RD/92/69

ARMINES Etude des ressources en eau souterraine de la rgion de Saint-Paul. Ile de la
Runion. Construction et exploitation du modle mathmatique (figures).
LHM/RD/92/69

CG974 Direction de lAgriculture et de la Fort (1995) Actualisation du modle
mathmatique du systme aquifre de la Plaine des Galets (texte)

CG974 Direction de lAgriculture et de la Fort (1995) Actualisation du modle
mathmatique du systme aquifre de la Plaine des Galets (figures)


Sites Internet :

http://www.eaureunion.fr/

http://infoterre.brgm.fr/