Simulation des transferts de polluants

dans les sols

Etat de l’art

RECORD 95-0504/1A février 1997 RECORD

 ETUDE N° 95-0504/1A

ETAT DE L’ART DE LA SIMULATION DES TRANSFERTS

DE POLLUANTS DANS LES SOLS

RAPPORT FINAL

février 1997

J. BOURGOIS - École des Mines de St Etienne

P. MOSZKOWICZ - LAEPSI - INSA de Lyon

www.record-net.org
Créée en 1989 à l’initiative du Ministère en charge de l’Environnement, l’association RECORD –
REseau COopératif de Recherche sur les Déchets et l’Environnement – est le fruit d’une triple
coopération entre industriels, pouvoirs publics et chercheurs. L’objectif principal de RECORD est le
financement et la réalisation d’études et de recherches dans le domaine des déchets et des pollutions
industrielles.
Les membres de ce réseau (groupes industriels et organismes publics) définissent collégialement des
programmes d’études et de recherche adaptés à leurs besoins. Ces programmes sont ensuite confiés
à des laboratoires publics ou privés.

Avertissement :
Les rapports ont été établis au vu des données scientifiques et techniques et d'un cadre réglementaire
et normatif en vigueur à la date de l'édition des documents.

Ces documents comprennent des propositions ou des recommandations qui n'engagent que leurs
auteurs. Sauf mention contraire, ils n'ont pas vocation à représenter l'avis des membres de RECORD.

Pour toute reprise d’informations contenues dans ce document, l’utilisateur aura l’obligation de
citer le rapport sous la référence :
RECORD, Etat de l’art de la simulation des transferts de polluants dans les sols, 1997, 129 p,
n°95-0504/1A.

Ces travaux ont reçu le soutien de l’ADEME (Agence de l'Environnement et de la Maîtrise de

l'Energie)
www.ademe.fr

© RECORD, 1997

Etude RECORD n°95-0504/1A

 Sommaire

ETAT DE L’ART DE LA SIMULATION DES TRANSFERTS DE

POLLUANTS DANS LES SOLS I
.
1 PRESENTATION DE LA MISSION ............................................................................................... 3
1.1. Introduction ....................................................................................................................................... 3
1.2. Objet de la mission ............................................................................................................................ 3
1.3. Contenu de la mission ....................................................................................................................... 4
I.3.1. Analyse de la demande ............................................................................................................... 4
1.3.2. Analyse de 1’oj5-e........................................................................................................................ 4
1.3.3. Elaboration de propositions ....................................................................................................... 4
.
II UN PEU DE GEOLOGIE ................................................................................................................ 5
II .1. Les roches......................................................................................................................................... 5
II.2. L’eau dans les roches ....................................................................................................................... 6
1’2.1. Comportement des roches vis-à-vis de l’infiltration ................................................................. 6
11.2.2. Les aquijères.............................................................................................................................. 6
11.2.3. Écoulement et alimentation des eaux souterraines ................................................................... 7
II.3. Etude d’un site hydrogéologique...................................................................................................... 8
II.3.1. La connaissance de 1’aquijèresur le terrain ............................................................................ 8
11.3.2. Où trouver les données hydrogéologiques ? ............................................................................. 9
II.4. Les sols et les sites industriels français............................................................................................ 9
II.4.1. Les grands types de sols ............................................................................................................ 9
II.4.2. Situation géographique des grandes zones industnellesfi.ançaises....................................... 10
.
IIl GENERALïïES SUR LES POLLUANTS DANS LES SOLS .................................................. 11
III.1. Le sol : un milieu poreux composé de trois phases ...................................................................... 11
III.2. Les relations entre le sol et l’eau .................................................................................................. 11
III.3. Les polluants ................................................................................................................................. 12
III.4. L’eau potable : gestion et normes................................................................................................. 14
.
IV TRANSFERT DE L’EAU ET DES POLLUANTS DANS LES SOLS ..................................... 15
IV.1.Fonctionnement général ............................................................................................................... 15
N.2. Transport de l’eau en milieu poreux............................................................................................. 16
N . 2 . I . La zone saturée ....................................................................................................................... 16
N . 2 . 2 . La zone non saturée................................................................................................................ 18
IV.3. Transport des solutés dans les milieux poreux ............................................................................. 19
N.3.1. La convection ......................................................................................................................... 20
N .3.2. La difusion ............................................................................................................................. 20
N.3.3. L equation générale de transport...........................................................................................
> I
21
IV.4. Les interactions entre le sol et le soluté........................................................................................ 22
N . 4 . 1 , Les mécanismes physico chimiques........................................................................................ 22
N . 4 . 2 . Les réactions microbiologiques ............................................................................................. 24
N.5. La notion de coefficient de répartition ......................................................................................... 25
IV.6. Etablir l’équation de transport ...................................................................................................... 26
.
V MODELISER POUR SIMULER .................................................................................................. 28
V .1. Qu’est-ce qu’un modèle ? .............................................................................................................. 28
V.2. Les problèmes d’échelle de travail ................................................................................................. 28
V.3. Les différentes approches de la modélisation ................................................................................ 29
K 3 . 1 . L ’approchedéterministe .......................................................................................................... 29
V.3.2. L ’approchestochastique.......................................................................................................... 31

-1-
 V.3.3. Les fonctions de transfert ........................................................................................................ 32
V.3.4. Les modèles type «boîte noire» ............................................................................................... 33
V.4. Les modèles en hydrologie ............................................................................................................ 33
V.5. Les modèles géochimiques ............................................................................................................ 34
V.6. Les différentes approches rencontrées au cours de la modélisation .............................................. 35
V.7. Conclusions.................................................................................................................................... 37
VI.SYNTHESE DE LA DEMANDE ................................................................................................. 38
VI.1.Motivations de la demande ........................................................................................................... 38
VI.2. La demande générale .................................................................................................................... 39
VI.2.1.Les principaux besoins ........................................................................................................... 39
VI.2.2. Comptes-rendus d’entretiens concernant la demande........................................................... 40
.
VII ANALYSE DE L’OFFRE ........................................................................................................... 43
VI1.1. Les modèles les plus utilisés et leur champ d’application .......................................................... 43
VII.1.1. Quelques modèles ................................................................................................................. 43
VII.2. Autres modèles existants sur le marché ...................................................................................... 47
VII.3. Les programmes en cours, groupes de travail, organismes de recherche et bureaux d’études ...49
VII.3.1. Quelques programmes en cours ............................................................................................ 49
VII.3.2. Groupe de réflexion .............................................................................................................. 51
VII.3.3.Laboratoires de recherches et écoles ................................................................................... 51
V11.3.4.Organismes de recherche. grands groupes et sociétés ......................................................... 57
.
VITI CONCLUSIONS ET PERSPECTIVES ................................................................................... 62
BIBLIOGRAPHIE .............................................................................................................................. 64
LISTE DES ANNEXES ....................................................................................................................... 66
LISTE DES FIGURES ........................................................................................................................ 67
LEXIQUE ............................................................................................................................................. 68
ANNEXES

-2-
 Présentation de la mission

1. PRESENTATION DE LA MISSION

1.1. Introduction

Les sociétés industrielles utilisent et produisent un très grand nombre de substances

qui sont susceptibles d’engendrer une pollution des milieux naturels. Une conséquence
particulièrement importante peut survenir par infiltration dans le sol de polluants solubles
dans l’eau, pouvant contaminer la nappe phréatique et menacer les ressources en eau potable.
La réglementation évolue pour garantir la préservation de cette ressource.
Il est donc nécessaire d’estimer l’impact d’une pollution accidentelle ou chronique,
réelle ou potentielle, en terme de risque. La modélisation, en permettant la simulation du
devenir des polluants dans les sols peut constituer un outil d’aide à la décision, soit en vue de
la prévention de la pollution de la ressource ou de la réparation de nuisances déjà existantes.
Dans un souci d’anticiper la demande réglementaire et en se conformant au principe de
précaution, l’association RJ3.CO.R.D. a initié cette étude pour faire un état de l’art des outils
permettant de simuler les transferts de polluants dans les sols. Des propositions de
programmes de recherches seront faites en vue d’améliorer les possibilités actuelles d’analyse
de risques pour les sols pollués.

1.2. Obiet de la mission

La mission a pour but de récolter les informations pertinentes sur les méthodes
existantes de simulation et de prévision de la migration des polluants dans les sols et plus
particulièrement pour les zones non saturées.
En effet, avant de lancer de nouvelles études concernant la modélisation du transfert
des polluants dans les sols il est nécessaire de disposer d’informations précises sur l’état
actuel des connaissances, les performances des outils existants et les travaux déjà engagés par
ailleurs (entreprises, bureaux d’études, laboratoires de recherche, programmes européens en
cours,...). De plus, il est nécessaire de faire apparaitre la nature des besoins concrets exprimés
par les membres de RE.CO.R.D.

Cette étude a été réalisée sous la responsabilité de:

Pierre Moszkowicz : Maître de conférence à 1’INSA de LYON : Laboratoire

d’Analyse Environnementale des Procédés et Systèmes
Industriels (LAEPSI).

Jacques Bourgois : Maître de recherche à 1’Ecole des Mines de Saint-Etienne

Département Ingénierie de l’Environnement.

-3-
 Présentation de la mission

avec les concours de:

Hervé Vaillant : Ingénieur de recherche à 1’Ecole des Mines de Saint-Etienne

Département Ingénierie de 1’environnement.

Ibrahim Alimi Ichola : Maître de Conférence à 1’INSA de LYON : Laboratoire de

Géotechnique.

Suzanne Foret : Stagiaire Ecole Nationale des Travaux hblics de 1’Etat de

Vaulx en Velin

1.3. Contenu de la mission

La mission est articulée en trois volets correspondant à l’analyse de la demande, de

l’offre, et à la formulation de propositions de programmes de recherche.

1.3.1. Analvse de la demande

Les experts sont chargés de rencontrer les interlocuteurs proposés par les différents
membres de l’association RE.CO.R.D. De ces entretiens, doit se dégager une formulation des
besoins exprimés par les membres de RE.CO.R.D.

1.3.2. Analyse de 1’off;e

L’analyse de l’offre doit permettre d’identifier les outils de calcul existants, en mettant
l’accent sur ceux qui visent à caractériser le transfert de l’eau et celui des polluants solubles
dans les milieux saturés et non saturés, en prenant en compte le couplage des phénomènes
physiques de transport (diffusion / convection) avec les interactions chimiques sols/polluants.
Plusieurs contacts avec des laboratoires spécialisés ont été conduits afin d’avoir la
vision la plus large possible sur les programmes de recherche et les préoccupations
scientifiques actuelles concernant la thématique de la modélisation des transferts de polluants
dans les sols.

1.3.3. Elaboration de vrovositions

L’ensemble de l’analyse de la demande et de l’offre doit permettre de dégager une liste

préliminaire des acteurs de la recherche. Les limitations empêchant aujourd’hui l’utilisation
satisfaisante des outils de simulation doivent être mises en évidence .
Ces résultats serviront à proposer des orientations qui seront discutées dans le cadre de
la programmation 1997 des programmes de recherche de RE.CO.R.D.

-4-
II. UN PEU DE GEOLOGIE

11.1. Les roches

Les roches constituant la croûte terrestre sont faites d’un assemblage varié de
minéraux. On peut classiquement diviser les roches en deux grands groupes :
- les roches qui se sont formées à la surface du globe : ce sont les roches sédimentaires,
- les roches qui ont pris naissance en profondeur dans des conditions de pression et de
température très différentes de celles qui règnent à la surface de la Terre. A l’intérieur de ces
roches, on sépare deux familles : les roches éruptives ou magmatiques et les roches
métamorphiques.

Les roches sédimentaires.

Ce sont les produits de sédiments meubles qui sont formés à la surface du globe, sur
les continents ou dans les mers. On les subdivise en trois grands groupes en fonction de leur
origine détritique, chimique ou biologique.
Les roches détritiques sont constituées de matériaux arrachés par l’érosion aux continents et
redéposés plus ou moins loin de leur source d’origine, dans des milieux variés sous des
formes diverses (sable, argile, conglomérats...).
Les roches d ’origine chimique résultent principalement de la précipitation physico-chimique
de divers sels dont l’eau de mer constitue le stock principal. C’est le cas des évaporites (gypse,
sel gemme, potasse...) et de certains calcaires et dolomites, de phosphates...
Les roches d’origine biologique, dont la génèse est liée indirectement ou directement à
l’intervention de la vie et des processus biologiques. Certaines se sont formées par
accumulation de coquilles ou de squelettes d’organismes qui ont puisé dans le milieu (la mer
le plus souvent) les éléments nécessaires à leur constitution (calcaire, silice, phosphates).
D’autres sont constituées par de la matière vivante (dite organique) élaborée par les
organismes. Ce sont les roches carbonées (charbons et hydrocarbures).

Les roches éruptives ou ignées.

Ces roches se forment par cristallisation d’un magma (d’où le terme de roche
magmatique) visqueux, constitué d’un mélange de gaz, de liquide et de cristaux, porté à haute
température (700 à 1200 OC). Le magma peut monter plus ou moins rapidement à la surface
du globe. Lorsque la montée est brusque, le magma se refroidit rapidement, et se trouve
comme figée lors de son arrivée en surface ou près de celle-ci. Ainsi se forment les roches
volcaniques comme les basaltes.
Lorsque le magma migre plus lentement et demeure à l’intérieur de l’écorce, il cristallise
progressivement et forme les roches plutoniques tel que les granites. Ces roches n’apparaîtront
à la surface qu’à la suite d’une érosion profonde.

Les roches métamorphiques.

Dans certaines conditions tectoniques et/ou d’enfouissement, les roches sédimentaires
ou ignées peuvent se transformer en d’autres roches sous l’effet d’une élévation de la
température et de la pression. Elles sont alors dites métamorphiques.
Les variations de température et pression engendrent de nouvelles structures (organisation en
feuillet) et conduisent à l’apparition de nouveaux assemblages minéraux dans les roches
(schistes, mica-schistes, gneiss...).

-5-
 Un peu de géologie

11.2. L’eau dans les roches

L’eau ne peut pénétrer dans les roches que dans la mesure où il existe des vides, et elle
ne peut circuler que si ces vides sont en communication les uns avec les autres. La porosité
exprime le volume de vides ; la perméabilité mesure l’aptitude d’une roche à se laisser
traverser par l’eau.

11.2.1. Comportement des roches vis-à-vis de 1l nfiltration

La pénétration de l’eau dans le sous-sol est conditionnée par la perméabilité des roches
qu’elle y rencontre. En fonction de ce critère, on sépare les roches en trois catégories:

- Roches poreuses mais imperméables. Ce sont des roches largement ou totalement

argileuses, où l’eau est solidement fixée par adsorption dans des pores de très petites
dimensions : la circulation y est pratiquement nulle.

- Roches poreuses, perméables en petit. Ces roches laissent circuler l’eau à travers
les pores interconnectés de la roche. Ce sont elles qui renferment les réserves d’eau
mobilisables les plus fréquentes. Certaines sont meubles (sables, graviers, cailloutis, scories
volcaniques...) mais on rencontre également dans cette catégorie des roches qui présentent, en
plus de la porosité interstitielle, une porosité de fissure (grès, craie, calcaires récifaux...)

- Roches non poreuses, perméables en prand. Pratiquement imperméables à l’échelle

de l’échantillon par suite de leur compacité, ces roches offrent à l’eau des voies de pénétration
localisées dans les fissures : diaclases, joint de cisaillement, de stratification, de schstosité,
fentes de retraits, etc. Toutes les roches dont la cohésion est élevée peuvent être rangées dans
ce groupe (granite, gneiss, calcaires, laves....). Certaines voient l’ouverture initiale de leurs
discontinuités s’agrandir par dissolution (roches carbonatées, gypse). Les roches carbonatées
sont très caractéristiques de ce type de dissolution qui conduit à une morphologie karstique
(cas des Causses, du Jura, de la Provence...).

11.2.2. Les aquifères

Les aquifêres sont des formations géologiques perméables permettant l’écoulement

d’une nappe souterraine ou le captage de quantité appréciable d’eau. Les relations de l’eau et
des terrains encaissants sont assimilables à celles d‘un liquide en mouvement dans un récipient
occupé par une roche poreuse.
Comme dans un récipient plus ou moins rempli, l’eau de la nappe contenue dans l’aquifêre
atteint un certain niveau appelé surface piézométrique ; cette surface ou ce niveau d’eau
fluctue dans le sens vertical en fonction de l’alimentationet de la vidange. On distingue quatre
grands types d’aquifêres :

1- Les nappes libres où l’eau occupe, à partir d’un substratum, tout l’espace vertical
sus-jacent, sans rencontrer d’obstacle imperméable limitant son développement vers la
surface.

-6-
 Un peu de géologie

2- Les aquifères fissurés et karstiques. Dans les roches peu poreuses perméables en
grand, une nappe d'eau peut se mettre en place, dans la mesure où les fissures et les fractures
sont en communication et où les fentes sont suffisamment nombreuses et larges. En général,
seules les zones de grandes fractures contiennent de l'eau en quantité notable. Dans les
aquifères karstiques, l'eau se concentre dans des vides agrandis par la dissolution de la matrice
calcaire par les eaux chargées de gaz carbonique acide. Ce type d'aquifère est caractérisé par
des variations considérables de la surface piézométrique et des débits.

3- Les nappes captives. Dans les aquifères plus profonds, les eaux souterraines sont
emprisonnées dans les formations hydrogéologiques perméables entre deux formations
imperméables fixes : le substratum en base et le toit au sommet. Cette formation est donc sous
pression et lorsque que l'on perce le toit de l'aquifère, l'eau s'en trouve expulsée. Son niveau se
stabilise à une altitude qui représente le niveau piézométrique. La pression peut être telle que
le niveau piézométrique peut se retrouver au-dessus de la surface (cas des puits artésiens).

4- Les aquifères superposés : La localisation des aquifères en profondeur est

déterminée par la succession des formations hydrogéologiques. Les formations imperméables
délimitent les aquifères. Ceux-ci peuvent ne comporter qu'une seule couche perméable
(aquifères monocouches) ou plusieurs couches (aquifères multicouches). Chaque aquifère est
caractérisé par une surface piézométrique particulière. Dans les aquifères multicouches, les
écoulements à travers les formations semi-perméables sont possibles grâce au mécanisme de
drainance.

11-2.3. Écoulement et alimentation des eaux souterraines

L'écoulement des nappes d'eau souterraines se fait depuis les zones hautes vers les
zones basses où se trouvent les exutoires. Les aires d'alimentation sont Caractérisées par un
bombement de la surface piézométrique. L'alimentation principale se fait par les eaux de pluie
qui s'infiltrent, mais une nappe peut être alimentée par une rivière ou par une autre nappe. De
façon générale, les nappes sont soumises à des écoulements liés à la différence de charge
hydraulique, entre les zones d'alimentation (zone haute) et les zones d'exploitation ou
d'émergence (zone basse).

L'alimentation par les précipitations se fait naturellement par les eaux météoriques et
par les eaux de fonte du manteau neigeux. Les fluctuations piézométriques des nappes ou du
débit des sources sont en liaison avec la lame d'eau d'infiltration réelle issue de la pluie après
prélèvement de l'évapotranspiration et du ruissellement.

Des relations hydrauliques existent entre les rivières et les nappes. En dehors des
périodes de pluie qui montrent une relation directe entre les précipitations et les débits
observés dans les rivières, on constate que dans notre pays, les rivières sont encore alimentées
en étiage (période de basses eaux). Ce débit de base provient des nappes souterraines qui ont
accumulé des eaux d'infiltration et les restituent progressivement aux rivières. En été
l'écoulement se fait donc des nappes vers les rivières; et en hiver le sens de l'écoulement
s'inverse lorsque la rivière en crue est plus haute que la nappe, elle alimente alors le réservoir
aquifère.

-7-
 Un peu de géologie

11.3. Etude d’un site hvdroeéologique

Beaucoup de pollutions affectent aujourd’hui les eaux souterraines. Ces pollutions sont
d’autant plus préocupantes que dans la plupart des cas, l’homme exploite ces ressources
souterraines par l’intermédiaire des captages d’eau potable. Une bonne connaissance de la
géologie et de l’hydrogéologie permet une meilleure prévention des risques de pollution et une
efficacité accrue en ce qui concerne les traitements nécessaires. Les caractéristiques d’un
aquifère (débit d’exploitation, vitesse d’écoulement, qualité des eaux, type d’alimentation,
géologie et structure du terrain...) sont les données essentielles pour comprendre et décrire les
transferts des polluants transportés par l’eau dans les sols.

II.3.1. La connaissance de I ’aquifèresur le terrain

Les études hydrogéologiques réalisées dans la plupart des cas à des fins d’exploitation
de captages, constituent une source de renseignements très précieuse pour la connaissance de
la géologie, de l’aquifêre, et des relations qui les lient.
L’approche hydrogéologique d’un site passe par la connaissance géologique des sols et
des roches concernés (identification des formations, type de la formation productive, structure
globale de la région, système de failles). L’agencement et la nature géologique du terrain
donne en général, les grandes notions de l’hydrogéologie du site (emplacement de la nappe,
type d’aquifère, type d’alimentation, sensibilité...).
La connaissance précise de l’aquifere nécessite la mise en oeuvre de techniques plus
lourdes et l’installation de sondages et puits d’exploitation par une société de forage, sollicitée
en général par le bureau d’étude ou l’hydrogéologue qui réalise le projet.
Ce sont les piézomètres, forages de très faible diamètre (une dizaine de cm en
moyenne), qui sont réalisés en premier lieu. Ils permettent l’investigation du sol et de
l’aquifère par l’intermédiaire de différentes sondes qui sont descendues à l’intérieur du forage.
Ces sondes, suivant leur fonctionnement (émission de rayon gama, résistance électrique...)
donnent accès a des informations qualitatives sur les roches en profondeur (perméabilité,
nature de la roche, compacité...) et le niveau de la nappe dont les variations sont suivies
directement dans le piézomètre. Toutes ces informations, indispensables pour l’implantation
et l’exploitation d’un captage AEP (eau potable), constituent une bonne base de données
facilement consultable.

C’est l’exploitation d’un tel ouvrage, qui passe par une série d’essais de pompage et
permet d’obtenir les paramètres quantitatifs et qualitatifs spécifiques de l’aquifère et de
visualiser et comprendre son comportement global. On connaît alors de facon très précise le
débit exploitable, la réponse de la nappe, l’influence de l’alimentation, la qualité des eaux, et
on peut donc évaluer le temps de transmission d’une pollution accidentelle déversée en amont.

-8-
 Un peu de géologie

113.2. Où trouver les données hvdropéolopiques ?

L’implantation d’un forage AEP nécessite des autorisations départementales et la

participation de plusieurs acteurs (privés, communaux, départementaux...). L’étude préalable
hydrogéologique est en général faite par un hydrogéologue agréé, ou par un bureau d’études.
Les hydrogéologues agréés dépendent d’une administration départementale (souvent la
DDASS) et sont seuls habilités a définir les périmètres de protection autour du captage. Ces
personnes sont sur une région, une source de renseignements de qualité. C’est auprès d’eux
que l’on peut obtenir des renseignements sur la géologie, sur la structure géologique de la
région et sur l’aquifêre.
Les captages sont en général gérés soit par la commune, soit par une filiale d’un grand
groupe (CGE, Lyonnaise...). La DDASS en association avec le ((fermien) (personne du groupe
qui s’occupe des captages), mesure la qualité de l’eau de façon mensuelle, veille au bon
fonctionnement des forages et au respect des périmètres de protection. On peut donc obtenir
auprès du fermier, de la DDASS, des hydrogéologues, tous les renseignements nécessaires sur
les caractéristiques de l’aquifére, utiles pour une modélisation de celui-ci.
A ces sources de renseignements on peut ajouter également les bureaux d’études et les
sociétés de sondages, qui en général travaillent depuis longtemps sur une même région et qui
en connaissent donc bien les principales spécificités ; les résultats sont cependant facilement
accessibles. Enfin, il faut ajouter la base de données du sous-sol du BRGM.

11.4. Les sols et les sites industriels français.

11.4.1.Les arands types de sols

La France se divise en cinq grandes zones métamorphiques et plutoniques qui sont le

Massif Armoricain, le Massif Central, les Alpes, les Pyrénées et les Vosges.
Ces zones sont composées de roches très anciennes (plusieurs centaines de millions
d’années) autour desquelles, au fil des millions d’années, les variations du niveau marin, les
glaciations, la succession des climats, et la tectonique, ont généré toute une série de roches
sédimentaires. Ce sont ces roches qui forment nos grands bassins (Bassin Parisien et
Aquitain). Des zones d’effondrement à l’intérieur des vieux massifs ont créé des fossés tels
que le fossé rhénan et le sillon rhodanien dans lesquels s’écoule le Rhin et le Rhône. Ces
mêmes variations du niveau marin ont comblé ces fossés avec l’aide de la sédimentation
fluviatile et de l’érosion glaciaire et forment ainsi, des grandes zones sédimentaires telle que
la Provence ou les Causses.
Toutes ces régions sédimentaires, métamorphiques et plutoniques se sont trouvées plus
ou moins fracturées, bousculées suivant l’endroit, la force et la durée des événements
tectoniques. La France géologique se couvre donc ainsi d’un réseau complexe de failles à
toute échelle.

Les grands ensembles sédimentaires vont à la fois présenter des roches calcaires,
manieuses, sableuses... qui proviennent d’une accumulation plus ou moins importante en
domaine océanique et en même temps des dépots fluviatiles et glaciaires (galets,
conglomérats, alluvions...p our les fleuves, moraines et autres pour les glaciers).
Les cratons sont composés de la vieille croûte, que l’on retrouve d’ailleurs sous les
zones sédimentaires à plus ou moins grande profondeur. Ces zones d’émergence de vieux

-9-
 Un peu de géologie

terrains sont un assemblage de roches métamorphiques, de protusions granitiques mise à nue

par l’érosion et de roches volcaniques plus ou moins récentes.

L’organisation structurale de ces massifs anciens est en général très complexe car elle
découle d’une succession de mouvement tectono-métamorphique générant différentes
conditions de pression et température au cours du temps. Ces zones sont pour la plupart très
peu perméables. Les ressources en eau souterraine sont faibles : on rencontre de petites nappes
locales dans les zones d’altération ou dans les dépots alluviaux de surface ou encore dans les
zones fissurées. Pour avoir de gros débits, il faut faire appel aux eaux de surface et ce, en
période hors étiage qui peuvent être sévères.
Dans les zones sédimentaires, il en est tout autrement puisque l’on s’adresse alors à
des roches perméables à différentes échelles, chaque type de dépots étant caractérisé par un
certain nombre de paramètres variables au cours du temps et dans l’espace, ainsi qu’un
cortège de réactions chimiques et biologiques encore mal déterminées. L’agencement des
nappes va dépendre du milieu sédimentaire auquel on s’adresse.
+ Dans les régions de plateaux calcaires (karst), c’est le domaine de l’écoulement dans les
fissures élargies, le réseau de fissures peut jouer un rôle de réservoir ou nappe en liaison avec
les eaux de surface et la rivière souterraine. Ces réservoirs sont en général situés en dessous de
la côte d’altitude des exutoires.
+ Les régions de bassins sédimentaires sont en général constituées par une alternance de
roches perméables et peu perméables. C’est là que l’on trouve les plus grands aquifêres.
+ Les plaines alluviales vont renfermer les nappes alluviales en relation avec les fleuves. Ces
nappes sont souvent les plus exploitées et souvent en danger d’être contaminées par des
pollutions accidentelles ou venues du fleuve qu’elles côtoient.

II.4.2. Situation aeopravhiaue des mandes zones industrielles fianeaises.

Les grandes régions industrielles sont généralement localisées sur les zones
sédimentaires où s’écoulent les grands fleuves, voies de communication. Le faible relief de
ces régions permet l’installation des voies routières et ferroviaires à moindre coût, à l’opposé
des zones de montagnes (vieux cratons en général) qui nécessitent tunnels, terrassements, et
mesures de protection géotechniques accrues contre les effondrements ou les avalanches).

Globalement, on retrouve donc les concentrations urbaines et industrielles, soit sur les
zones de faibles reliefs, présentant des voies de communication aisées, soit dans les régions de
dépots sédimentaires qui sont la plupart du temps, assez sensibles à cause de leur forte
perméabilité et de leur fonction de stockage d’eau. On peut remarquer quelques exceptions
telles que les vallées de la Maurienne et de la Tarentaise qui sont en zone de montagne mais
néanmoins en fond de vallée et donc sur des dépots glaciaires et fluviatiles de faible épaisseur.

-10-
 Généralités sur les polluants dans les sols

III. GENERALITES SUR LES POLLUANTS DANS LES SOLS

111.1. Le sol : un milieu poreux composé de trois phases

La charte européenne des sols établie par le conseil de l’Europe le 30 mai 1972, définit
le sol comme suit :
R Le sol est un milieu vivant et dynamique qui permet l’existence de la vie
végétale et animale. Il est essentiel à la vie de 1’homme en tant que source de nourriture et de
matière première. Il est un élément fondamental de la biosphère et contribue, avec la
végétation et le climat, à régler le cycle hydrologique et à influencer la qualité des eaux.»

Suivant le spécialiste, la définition d’un sol varie. Pour l’ingénieur, le sol est avant tout
une roche meuble, un support, un milieu poreux. Pour l’agriculteur, c’est la couche de terre
travaillée par opposition au sous-sol qui ne l’est pas. Pour le pédologue, au plein sens du
terme, qui étudie l’ensemble des matériaux et l’activité biologique, c’est l’aspect biologique et
géochimique qui domine. Le géologue travaille à l’échelle du globe (tectonique des plaques),
de la croûte terrestre (histoire et condition de formation des chaînes de montagnes) comme à
celle des premiers mètres de terrain (dépôts sédimentaires superficiels) ou encore de la
molécule (étude de la spéciation moléculaire à haute pression). L’hydrogéologue travaille sur
un ensemble restreint à l’échelle d’une région géographique que l’on appelle le bassin versant
et va considérer les relations eau - sol dans cette zone (quantification, mais aussi mode de
circulation, type de terrain, qualité des eaux et protection, possibilité d’exploitation. Enfin
l’hydrologue va considérer les phénomènes physico-chimiques concernant l’eau et ses
interactions avec le sol dans la tranche superficielle (quelques dizaines de mètres par rapport à
l’approche du géologue).

D’un point de vue physico-chimique, qui est l’approche de la modélisation concernant

les transferts de polluants, le sol peut être considéré comme un système poreux à trois phases :
la matrice du sol, une phase liquide composée d’eau et de solutés, une phase gazeuse
composée d’air et de vapeur d’eau ainsi que de composés organiques volatils. Chacune de ces
phases possède des caractéristiques physiques et mécaniques propres, utilisées pour la
modélisation afin de les définir en quantité et en qualité.

111.2. Les relations entre le sol et l’eau

Le sol se compose :
- d’une zone non saturée composée elle-même de trois zones : la zone d’évaporation
et d’évapotranspiration (la plus proche de l’air libre), la zone de transition où il y a un transfert
d’eau par infiltration et la partie non saturée de la fiange capillaire (partie au-dessus de la
nappe où se trouve l’eau capillaire, seulement retenue par les forces capillaires) où la pression
de l’eau est inférieure à la pression atmosphérique. L’eau remplit incomplètement les
interstices, le reste du volume étant occupé par l’air,
- d’une zone saturée où l’eau occupe cette fois complètement la porosité du sol. L‘eau
de cette zone est libre en ce sens qu’elle peut être mise en mouvement par l’effet de la
gravité.ou de charges hydrauliques Cette notion regroupe à la fois ce que l’on appelle la

-1 1-
 Généralités sur les polluants dans les sols

nappe souterraine et en même temps une partie de la frange capillaire qui peut également être
saturée.
L’eau en présence dans les sols non saturés peut avoir différents comportements. On
distingue l’eau de constitution du sol qui appartient à la structure même des composants et
qui est donc ni échangeable, ni mobile, l’eau liée, l’eau capillaire et l’eau gravitaire.

L’eau liée ne peut pas se déplacer par gravité et ne transmet pas la pression
hydrostatique. C’est une eau qui ne peut être extraite que par séchage (thermique ou
mécanique). Elle est liée à la surface des grains par le jeu des forces d’attraction moléculaires
très intenses.

L ’eau capilzaire s’élève au-dessus de la nappe pour former la zone capillaire, cette eau
est insensible à l’action de la gravité et est retenue par les seules forces capillaires du sol. Ces
forces sont uniquement présentes dans les milieux poreux non saturés et sont des forces
motrices pour le mouvement de l’eau. Afin d’illustrer de façon simple la notion de capillarité,
il suffit de considérer la montée de l’eau dans un mur à partir d’un horizon mouillé en base ;
la hauteur de montée d’eau est proportionnelle à la pression, ou succion qu’exercent les forces
capillaires du milieu.

L ’eau gravitaire ou eau libre remplit l’espace resté libre entre les pores, interstices ou
fissures. Cette eau de percolation (proprement dite) peut se déplacer sous l’effet de la gravité
ou des gradients de pression. Elle est en dehors du champ d’attraction des particules solides
C’est la partie vraiment active des eaux souterraines.

111.3. Les polluants

La pollution des eaux souterraines et des eaux continentales en général s’effectue soit
par infiltration d’eaux météoriques chargée en polluants, soit par déversement direct de
polluants liquides à la surface du sol ou encore par contamination indirecte par les eaux de
ruissellement qui sont chargées en produits plus ou moins toxiques et qui peuvent atteindre la
nappe en traversant la zone non saturée ou par l’intermédiaire des cours d’eau. Après avoir
atteint la nappe le polluant est entraîné par l’écoulement de celle-ci vers l’aval et c’est au
cours de ce trajet qu’il voit sa concentration varier en fonction de la circulation des eaux
souterraines et des différents facteurs biophysicochimiques pouvant affecter sa mobilité et sa
toxicité.
On peut lister les principales familles de polluants : radioactifs, chimiques
(hydrocarbures et produits de combustion, matières plastiques, pesticides, détersifs, dérivés du
soufre, nitrates, phosphates, métaux lourds, fluorures, composés organiques,...), biologiques
(matières organiques mortes, micro-organismes pathogènes).
Tous ces polluants peuvent être directement ou indirectement très toxiques. Ils peuvent
agir suivant différents modes de propagation (sol, air, eau), par concentration progressive dans
les organismes vivants et ainsi altérer tout ou partie d’une chaîne alimentaire (phénomènes de
bioamplification, biomagnification).
Les principales sources de polluants proviennent de la production et de l’utilisation des
diverses sources d’énergie, des activités industrielles et de l’agriculture. Le souci de la
pollution des sols apparaît dans la législation Concernant les déchets dès 1975 puisque la loi
du 15 juillet 1975 vise les déchets «de nature à produire des effets nocifs sur le sol, la flore et
la faune, à dégrader les sites ou les paysages, à polluer l’air ou les eaux)).

-12-
 Généralités sur les polluants dans les sols

L’urbanisation accélérée des pays industrialisés a eu pour effet une multiplication des
sources de contamination de l’environnement. L’expansion considérable de la production
industrielle a engendré des masses énormes de déchets minéraux ou organiques,
biodégradables ou non, parfois très toxiques. La pollution des sols va généralement provenir
de xénobiotiques alors que pour la pollution de la ressource en eau, on a une forte composante
microbiologique couplée aux xénobiotiques.

La production et 1 ’utilisation d ’énergie constitue une source importante de pollution

de la biosphère en général et des sols en particulier. Dans cette catégorie, vont rentrer les
produits pétroliers, mais aussi le charbon et le gaz. L’énergie nucléaire constitue un cas
particulier: la gestion des risques de dissémination des radioéléments porte essentiellement sur
les mesures de sécurité du fonctionnement des réacteurs et sur la gestion à long terme des
déchets. De très lourdes études sont actuellement engagées sur le stockage profond des
déchets à longue période. L’extraction et la combustion des produits pétroliers
s’accompagnent d’innombrables pollutions directes ou indirectes tel que les marées noires ou
les accidents de transports qui contaminent et polluent les eaux continentales et les sols, ainsi
que les vidanges sauvages et autres usages dispersifs des hydrocarbures. Les goudrons
résultant de l’exploitation du gaz à partir du charbon, sont à l’origine de pollutions très
importantes qui concernent plusieurs centaines de sites.
En définitive, la boulimie d’énergie, Caractéristique des pays industrialisés,
s’accompagne d’un risque de contamination sans cesse accru des eaux continentales et des
sols par les innombrables substances polluantes produites par les combustions et usages
annexes.

L ’industrie chimiaue, mais aussi la métallurgie, voire 1 ’électronique, mettent en

circulation de nombreux composés minéraux edou organiques de toxicité souvent élevée et
peu dégradables ou indestructibles (mercure, cadmium, vanadium. La chimie organique
produit d’innombrables molécules de synthèse plus ou moins toxiques. On considère que tout
au plus, le tiers de ces substances ont fait l’objet d’une étude sérieuse sur leurs impacts
potentiels sur l’environnement et sur l’homme.

La civilisation moderne produit également des masses colossales de déchets solides

qui peuvent se classer de différentes manières. On pourra distinguer les déchets domestiques,
agricoles et industriels. Les deux premières catégories sont constituées en partie par des
matières organiques donc biodégradables, mais la gestion des déchets ultimes constitue
toujours une préoccupation vis-à-vis de la pollution des sols où sont situés les sites de
stockage.

L ’agriculture moderne représente une importante source de pollution des sols cultivés
et des ressources en eau potable. L’usage massif d’engrais chimiques et le recours
systématique aux pesticides ont malheureusement été accompagnés d’une pollution accrue des
sols et des eaux continentales ainsi que des productions végétales et animales par divers
contaminants minéraux et organiques. L’abus de fertilisants a été tel que dans de nombreux
pays la pollution des nappes phréatiques atteint localement des niveaux excédant les
concentrations admissibles en nitrates. Les métaux et métalloïdes toxiques (cadmium,
vanadium, chrome, cuivre, arsenic, etc.) contenus comme impuretés dans les superphosphates,
s’accumulent dans les sols et peuvent passer dans les plantes cultivées.
Les pesticides et les insecticides, qui ont connu une expansion spectaculaire, possèdent
également un fort pouvoir biocide et donc une haute toxicité pour une multitude d’organismes

-13-
 Généralités sur les polluants dans les sols

que nous consommons par la suite. Les masses de ces substances dispersées dans l’espace
rural sont considérables si l’on regarde la toxicité et/ou la persistance de certaines d’entre
elles. L’insertion de ces produits dans les chaînes alimentaires, par l’intermédiaire des sols et
des cultures n’est plus à démontrer.

111.4. L’eau potable :gestion et normes.

Pendant longtemps, on a considéré rôle épurateur du sol comme suffisant. La pollution

présente en surface ne semblait pas affecter les eaux souterraines captées ou les eaux
superficielles. C’est seulement avec l’avènement de l’ère industrielle et de ses inévitables
pollutions que la société a commencé à se soucier du lien entre l’activité agricole, industrielle
ou domestique et la pollution des eaux de nappe et de surface en terme de potabilité, ou
d’autre usage. On s’est alors rendu compte que le sol ne suffisait plus pour épurer nos rejets,
qu’il pouvait être lui-même contaminé et contaminer à son tour, après un certain temps de
transfert, la ressource en eau. Il devenait donc essentiel d’établir un moyen de contrôle de la
qualité de l’eau.

La préservation de cette qualité ainsi qu’une meilleure gestion du patrimoine eau font
partie des grandes préoccupations au même titre que l’évolution des réserves en matières
premières. L’OMS (Organisation Mondiale pour la Santé) a édicté des valeurs guide
internationales pour l’eau de boisson entre autre, qui comprennent un exposé général des
valeurs guide bactériologiques, physiques, chimiques, biologiques et radiologiques. La plupart
des pays ont établi des normes de qualité en référence à ces valeurs guide. En France, la loi
impose de façon stricte que «les eaux ne doivent présenter aucun signe de dégradation de
qualité et doivent satisfaire aux exigences de qualité définis dans le décret n”89-3 (3/01/89)»
(annexe 1).
Ces exigences ont été récemment réactualisées par la Loi sur l’Eau (03/01/92). Les
objectifs définis sont de renforcer les mesures de contrôle et de sanctions pour la protection de
la ressource en eau et des milieux aquatiques (entamé en 1964) et de mettre en place une
véritable gestion intégrée du patrimoine eau (cf Code Permanent Environnement et
Nuisances). Il faut souligner que les normes varient en fonction de l’usage : les paramètres
observés, et les exigences ne seront pas les mêmes suivant que l’eau est destinée à un usage
domestique ou industriel.
Le principe de gestion intégrée fait peu à peu son chemin dans les esprits et on assiste
maintenant à une collaboration de plus en plus fréquente entre les différents acteurs d’un
même bassin, souvent en conflit d’usage (agriculteurs, entreprises, communes, habitants, état)
en vue de la protection de la ressource en eau souterraine et superficielle.

-14-
 Transfert de 1*eauet des polluants dans les sols

IV. TRANSFERT DE L’EAU ET DES POLLUANTS DANS LES SOLS

Un soluté est une substance soluble dans l’eau et qui se meut avec elle par transport
convectifet dz#ksifmais qui peut aussi interagir avec le milieu dans lequel la solution (eau +
soluté) se déplace. La solution aqueuse est à la fois vecteur et milieu réactionnel. Les
transfm de solutés impliquent différents mécanismes bio-physico-chimiques fortement
couplés et très complexes.

IV.1. Fonctionnement général

Lorsque l’on parle de transfert de polluants dans un sol, on considère dans la plupart
des cas, les polluants solubles dans l’eau (cependant les Composés Organiques Volatils
peuvent aussi être transportés dans la phase gazeuse). L’eau est donc très importante en tant
que vecteur mais aussi en tant que cible potentielle ou réelle. Le schéma suivant rend quelque
peu compte de la complexité des phénomènes qui interagissent entre eux dans le cas
particulier des produits phytosanitaires :

SOL NON SA-! -1

T
-
(VADOSE)

Figure 1 :Cycle des produits phytosanitaires (d’après Vachaud and Vauclin, 1992)

Une fraction du produit répandu, est assimilé par le sol et les plantes. La fraction
restante, volatilisée à la surface du sol et du feuillage des plantes, retourne à l’atmosphère.
Une partie de ces produits sous forme vapeur, est susceptible de se dissoudre dans l’eau de
pluie et d’être réinfilitrée dans le sol, pouvant ainsi former des zones d’accumulation. Lorsque
la pluie est trop intense pour être entièrement infiltrée dans le sol ou que la zone non saturée
disparaît par suite de la remontée d.e l’horizon saturé, le ruissellement apparaît. Les polluants
accumulés dans les couches supérieures du sol peuvent alors être emportés par les eaux de
ruissellement, soit sous forme dissoute, soit sous forme adsorbée quand l’eau de ruissellement
transporte des particules fines sur lesquelles ces produits sont fixés [Guinot, 19951. La partie

-15-
 Transfert de l’eau et des polluants dans les sols

des polluants qui parvient à s’infiltrer est soumise à de nombreux phénomènes physiques,
chimiques et biologiques qui sont détaillés ci-après.

IV.2. Transport de l’eau en milieu poreux

La pénétration et le transport de l’eau dans un milieu poreux sont globalement régis

par deux grands phénomènes : la charge hydraulique induite par les forces de gravité,
exprimée mathématiquement par la loi de Darcy et les phénomènes de capillarité (succion).
En fait, ce mouvement dépend de l’énergie de l’eau dans le sol. La définition de cette énergie
doit être liée à l’énergie libre de Gibbs qui permet d’exprimer les conditions d’équilibre ou les
conditions de changement d’état. L’eau dans le sol étant soumis à divers champs de forces,
son énergie résulte de la somme des potentiels de ces forces. Ce potentiel résultant
comprendra aussi bien le potentiel des forces gravitaires que le potentiel des forces capillaires,
que le potentiel des pressions osmotiques. La loi de Darcy sera généralisée en utilisant le
gradient de ce potentiel. Il convient de séparer les transports dans la zone saturée de ceux
dans la zone non saturée.

N.2.I. La zone saturée

I r 2.1.1.La loi de Darcy

La base fondamentale du calcul de quantité d’eau souterraine circulante (ou débit de

nappe), par l’hydrodynamique souterraine, est la loi empirique de Darcy établie en 1856. Avec
un dispositif expérimental, Darcy a montré que le débit d’eau Q en m3.s-l, filtrant de haut en
bas dans une colonne de sable de 1 m de hauteur, à travers une section totale perpendiculaire à
la direction de l’écoulement, A en m2, est fonction d’un coeficient de proportionnalité, K en
m.s-’ , caractéristique du milieu perméable et de la perte de charge par unité de longueur du
cylindre de sable, h/l (sans dimension).
On définit alors le débit Q fonction de :

Q(m3.s-l)= K(m.s-’).A(m 2).-h

1

Cette loi, établie sur des dispositifs particuliers, a été vérifiée expérimentalement au
laboratoire, dans toutes les conditions possibles. Ces études ont permis d’établir une
expression généralisée applicable dans tous les cas, quelle que soit la direction x
d’écoulement. Soit :

où * q représente le volume d’eau écoulé par unité de surface et par unité de temps soit le
débit unitaire en m.8’.
* K le facteur de proportionnalité généralement appelé conductivité hydraulique
(m.s-‘) ou coefficient de perméabilité.
* !!! est le gradient de charge hydraulique (sans dimension).
dx

-16-
 Transfert de I’eau et des poIIuants dans les sols

IV.2.1.2.Définitions

La conductivité hydraulique ou coefficient de perméabilité K permet de quantifier

la perméabilité qui est l’aptitude d’un réservoir à se laisser traverser par l’eau, sous l’effet
d’un gradient de charge hydraulique. Le coefficient de perméabilité définit le volume d’eau
gravitaire (en m3) traversant en une unité de temps (s), sous l’effet d’une unité de gradient
hydraulique, une unité de section (en m2) orthogonale à la direction de l’écoulement. Ce
paramètre a la dimension d’une vitesse. Des méthodes de mesure de la conductivité existent
en laboratoire sur des échantillons mais aussi in situ.

D’autre part, on distingue plusieurs grandeurs à ne pas confondre :

-+La charge hydraulique H qui représente le poids de la colonne d’eau, exprimée en
m de hauteur d’eau au dessus du niveau de référence.
+ Le gradient de charge hydraulique est défini comme le quotient de la charge H
appliquée par la longueur de la colonne de sable 1 (expérience de Darcy). Ceci représente la
perte de charge par Unité de longueur (sans dimension).

N.2.1.3.Les limites d’application de la loi de Darcy

La loi de Darcy n’est pas toujours vérifiée en écoulement saturé. Cette relation n’est
plus linéaire aux grandes vitesses d’écoulements pour lesquelles les turbulences ne sont plus
négligeables. La loi de Darcy est donc applicable seulement pour les écoulements laminaires
(c’est à dire qu’il y a glissement non turbulent de couches fluides superposées et parallèles) et
dans les cas où l’interaction sol-eau n’introduit pas de changements dans la fluidité ou de
perméabilité en fonction du gradient. L’écoulement laminaire est prépondérant dans les
limons et les matériaux plus fins, mais dans les sables grossiers et les graviers de porosité plus
importante, on peut avoir des conditions non laminaires et la loi de Darcy n’est alors pas
toujours applicable.
Le critère quantitatif permettant l’évaluation du domaine turbulent est le cnombre de
Reynolds)) :
Re=-dup
CI

où d est le diamètre effectif des pores, la vitesse moyenne d’écoulement, p la densité du

liquide et 1.1 sa viscosité dynamique.
Dans les milieux poreux, le flux reste linéaire par rapport au gradient hydraulique aussi
longtemps que Re est inférieur ci l’unité. Au fur et à mesure que la vitesse d’écoulement
augmente, particulièrement dans les systèmes composés de grands pores, l’apparition de
tourbillons turbulents contribue à la dissipation d’énergie effective. On observe alors des
phénomènes qui s’écartent de la loi de Darcy. Il faut néanmoins souligner que ces
phénomènes sont généralement d’une faible importance dans la pratique, et la loi de Darcy
peut être utilisée dans la grande majorité des cas en ce qui concerne l’écoulement de l’eau
dans un sol saturé [Hillel, 19881.

-17-
 Transfertde l’eau et des uolluants dans les sols

IV.2.2. La zone non saturée

IV.2.2.1.L’équationde Darcy appliquée au milieu non saturé

Les équations du milieu non saturé découlent de celles de la zone saturée. La

différence réside dans le fait que certains paramètres deviennent variables au cours du temps
et dans l’espace. C’est le cas de la conductivité hydraulique K qui, en zone non saturée,
dépend à la fois des caractéristiques du sol (sols plus ou moins compacts, graveleux, sableux,
dimension des pores, état hydrique du sol....) mais aussi de celles du fluide (densité, viscosité).
L’état insaturé du sol entraîne l’apparition du phénomène de capillarité qui se traduit par la
succion. On peut imager la succion par la hauteur d’eau qui va monter dans un milieu poreux à
partir d’un horizon saturé. On l’exprime en m d’eau.
En tenant compte de ce phénomène, on peut considérer K comme variable par rapport
à Y (potentiel de succion).

dH
Dans la zone saturée, le flux est induit par un gradient de charge hydraulique - qui
dx
représente un potentiel hydraulique. Dans la zone non saturée, en absence de charge
hydraulique, on peut avoir un flux qui cette fois résulte d’un gradient de concentration en eau
(ou d’un gradient de potentiel de succion).

Richards en 1931, étend donc la loi de Darcy conçue pour l’écoulement en milieu
saturé, à la zone non saturée, en stipulant que la conductivité est une fonction de la charge de
succion. La relation de Darcy devient alors l’équation de Richards écrite ici pour un
écoulement unidirectionnel:

où 9
dx
est le gradient de succion
On peut écrire l’équation de Richards, lorsque l’on est en présence à la fois de la
succion et d’une charge hydraulique, sous la forme suivante :

où H représente la charge hydraulique qui peut inclure les termes de gravité et de succion.
On peut bien écrire ces relations vectoriellement pour les trois directions de l’espace.

IV.2.2.2.L’approche difliionnelle en zone non saturée

On peut décrire le flux d’eau en milieu non saturé par l’intermédiaire de la diffusion au
lieu de la succion. En effet les gradients de teneur en eau peuvent être considérés comme les
potentiels induisant des transports diffusifs.

-18-
 Transfert de l’eau et des polluants dans les sols

On exprime mathématiquement la diffusion par une loi phénoménologique, la loi de

Fick, qui considère que le flux diffusif q est proportionnel au gradient de concentration du
produit :

q = -D,Grad 8

Où 8 est la teneur en eau

De, le coefficient de difision (m2.s-’)ou la diffusivité de l’eau.

Le coefficient de difision est alors lié à la conductivité hydraulique par le fait que 8
et Y sont eux-mêmes liés. Le flux peut donc s’écrire de deux manières suivant que l’on
considère la succion ou le mouvement par difision :

8 la teneur en eau volumique peut être exprimée en pourcentage; c’est la quantité

d‘eau (gravitaire, liée) contenue dans un reservoir, saturé ou non, rapporté au volume total,
quel que soit le type d‘eau considérée. Dans la zone non saturée, on parle d’humidité définit
par :
volume deau
0=
volume total

On peut alors relier les coefficients de difision et de perméabilité :

De(€))= K(O)$

Remaraue :
La loi de Darcy définie à l’aide du potentiel de succion n’a de sens que si le sens de
l’évolution hydrique du sol est connue. En effet, le phénomène d’hystérèse de la relation entre
la succion et la teneur en eau, oblige à préciser le sens de variation de la teneur en eau pour
fixer cette relation. La fonction K(w) doit être définie soit pour le mouillage (infiltration), soit
pour le drainage. Ceci complique sa détermination. L’hypothèse faite par Richards (193 1), en
posant la conductivité hydraulique fonction de la teneur en eau, ne peut être traduite par la
relation ~ ( w )que si la relation y@) est univoque. En précisant l’histoire de l’évolution de
l’humidité dans le sol, la relation (IV.4) a alors un sens. Sinon, l’équivalence formelle des
deux écritures de la loi de Darcy ne peut être justifiée.

W.3. Transport des solutés dans les milieux poreux

Le transport des solutés présents dans l’eau est régi par les lois physiques de transport
de l’eau que nous venons de voir. Les solutés sont ainsi véhiculés de façon «passive»c’est à

-19-
 Transfert de 1 ’eau et des polluants dans les sols

dire sans interactions d’ordre chimique, par des mécanismes de convection. En outre les
différences de concentration d’un soluté dans l’eau induisent un transport par diffusion.

IV.3.I . La convection

Ce phénomène est le plus simple que l’on puisse considérer et l’hypothèse implicite est
la non-interaction entre les solutés (le polluant) et le sol. Le soluté suit donc le mouvement
général de l’eau sans perturbation d’origine chimique. Au mouvement général du soluté,
s’ajoute le mécanisme de dispersion cinématique. Ce phénomène dépend de la vitesse
d’écoulement d’un fluide dans un milieu poreux. La dispersion est liée à l’hétérogénéité des
vitesses microscopiques lorsqu’un fluide s’écoule à travers un milieu poreux; sa distribution
de vitesse n’est pas uniforme. Cela est dû :
* au profil de vitesse qui, dans les capillaires intergrains est parabolique,
* à la dimension des pores,
* aux lignes de courant qui fluctuent par rapport à la direction principale de
l’écoulement en raison des obstacles constitués par les grains ou matériaux plus grossiers
[Dupont, 19941.

IV.3.2. La diffusion

La diffusion recouvre l’ensemble des phénomènes qui tendent à homogénéiser les

concentrations dans un fluide où existent des gradients de concentration. A la dispersion
cinématique décrite précédemment, résultant de l’écoulement du fluide , on ajoute la
diffusion moléculaire, phénomène physique qui provient de l’agitation moléculaire. On
utilise souvent le terme de diffusionpour parler à la fois de la dispersion cinématique et de la
diffusion moléculaire. La diffusion moléculaire préexiste dans un fluide au repos et provient
de la distribution non homogène du soluté dans la phase fluide. Le soluté va migrer des zones
plus concentrées vers les zones moins concentrées afin d’établir une distribution de
concentration uniforme. Ce mouvement est généré par la pression osmotique, commandé par
le gradient de concentration qui varie d’une zone à l’autre du milieu poreux.
On modélise le phénomène de diffusion (dispersion cinématique et diffusion
moléculaire) par la loi de Fick explicitée précédemment où on considère que le coefficient de
diffusion D (coefficient de diffusion apparent), est la somme du coefficient de dispersion
cinématique Dcin et du coefficient de diffusion moléculaire Dmol.

La dispersion cinématique étant liée à la vitesse d’écoulement, Dc, devient nul quand
la vitesse de l’écoulement tend vers O. Le coefficient de dispersion cinématique est lié à la
vitesse de l’écoulement par une loi semi-empirique :

où V est la vitesse de l’écoulement.

Le coefficient h est appelé la dispersivité et a la dimension d’une longueur. Sa valeur
déterminée au laboratoire est de quelques centimètres; in situ elle de l’ordre du mètre voire

-20-
 Transfertde 1 ’eau et des Dolluants dans les sols

des centaines de mètres suivant l’hétérogénéité de la formation (Lallemand-Barrès et al. cités

par de Marsily). La dispersivité peut varier avec la direction de l’écoulement.
L’étude du transport dispersif nécessite la mesure du coefficient de diffusion
moléculaire et de la dispersivité. Le coefficient de diffusion du soluté dans le sol est différent
du coefficient de diffusion des ions dans l’eau appelé coefficient de diffusion libre. La valeur
du coefficient de diffusion moléculaire est réduite par la présence des pores partiellement
remplis et par la tortuosité du sol. Pour tenir compte des différents facteurs influençant la
diffusion des ions dans la sol, le coefficient de diffusion moléculaire est multiplié par
différents coefficients réducteurs. On obtient alors le coefficient effectif de diffusion. Le
tableau ci-dessous donne une idée des coefficients correcteurs proposés. Les techniques de
détermination du cofficient de diffusion du soluté permettent l’obtention du coefficent effectif
de diffusion. On peut alors obtenir le cofficient correcteur global appelé facteur d’impédance,
si on connaît par ailleurs le coefficient de diffusion libre.

Définition Référence
D* = ZDmol Gilham et al. (1984)
D* = TCtDmol Li et Grégory (1974)
D* = zeDmoi Berner (1974), Drever (1982)
D* = zayeDmol Kemper et al. (1964); Olse et Kemper (1968); Nye (1979)
D* = V D m o l Poter et al.
D* = zCtûDmo, van Shalk et Kemper (1966)

z facteur de tortuosité
c1 facteur de mobilité
y facteur dû au compactage
8 teneur en eau volumique
L’utilisation de la loi de Fick appliquée à la diffusion du soluté dans le sol permet la
détermination de D*.

IK3.3. L ’éauationaénérale de transvort

Si on combine la convection et la diffusion, on aboutit à l’équation de convection-

diffusion, écrite ici pour un écoulement unidirectionnel dans un milieu poreux saturé :

-
où c la concentration du soluté,
Dij le coefficient de diffusion de l’élément i dans la directionj,
-
et V j la vitesse d’écoulement dans la direction j.
Le premier terme (a) représente la diffusion tandis que le second (b) traduit la
convection. <D représente le terme source (apparition ou disparition de masse de soluté
durant le transport).
Le champ de vitesse Vj s’obtient par la résolution des équations de Darcy ou de
Richards.

-21-
 Transfert de 1’eau et des polluants dans les sols

La convection et la dispersion peuvent être considérées comme conservatives lorsque:

1) il n’existe pas de phénomène de sorption-désorption qui modifie la concentration en soluté
en un point de l’écoulement par les échanges avec la matrice poreuse.
2) tout l’eau se déplace; il n’y a pas de couche de molécules d’eau adsorbées ou retenues par
capillarité, qui constituerait une phase immobile susceptible d’influencer la migration du
soluté.
3) il n’existe pas de phénomènes de dissolutiodprécipitation qui modifieraient la
concentration de soluté en un point du milieu poreux.

IV.4. Les interactions entre le sol et le soluté

Le transport d’un soluté dans un milieu poreux est fonction de facteurs spécifiques au
fluide, à la matrice solide, au soluté lui-même ou encore aux différentes interactions existant
entre les différentes phases. Les phénomènes à prendre en compte sont la convection, la
diffusion qui englobe la dispersion cinématique et la difision moléculaire pour les
mécanismes de transport du liquide, auxquels s’ajoutent les phénomènes d’adsorption et
d’échange d’ions, les réactions chimiques et les processus biologiques (liste non exhaustive)
pour les mécanismes influençant la mobilité des espèces.
Dans ce type de processus, la loi de conservation de la matière n’est pas respectée dans
la phase aqueuse, c’est pourquoi ces phénomènes sont «non conservatifs)).Ils se superposent,
en terme de «source» et de <cpuits»aux équations de transport par convection - dispersion.
Ainsi, on identifie plusieurs processus pouvant conduire à l’immobilisation des solutés
(adsorption, précipitation chimique...) ou à leur mobilisation diférée.

IV.4.1. Les mécanismes uhvsico chimiaues

La rétention

Dans les sols non saturés, la proportion de l’eau liée plus ou moins fortement au
squelette peut être du même ordre de grandeur que celle de l’eau libre. Cette partition de la
phase liquide en phase mobile et phase immobile, influence le transport de soluté. La
différence de concentration dans la phase mobile et la phase immobile provoque un échange
entre les deux phases qu’on appelle (<échangelatéral ». Pour tenir compte de cet aspect du
transport, on introduit le coefficient de partage exprimé en pourcentage d’eau mobile. Soitf=
8,/6, où 8, est la teneur en eau mobile et 8 est la teneur en eau volumique total du milieu
poreux.
A chaque phase est alors associée la concentration appropriée: Cm pour la phase
mobile et Ci, pour la phase immobile. L’équation (IV.6) de convection-dispersion devient :

Cette équation suppose constant, le pourcentage d’eau mobile dans le poreux.. Cette
équation comporte deux inconnues ; il faut donc lui adjoindre une relation suplémentaire pour
la résolution.

-22-
 Transfert de l’eau et des polluants dans les sols

L ’adsomtion

L’adsorption de molécules de soluté à la surface du solide provoque leur

immobilisation. Ce phénomène résulte des forces d’attraction ou de répulsion existant entre
les ions ou les molécules des différentes phases.
On distingue :
- l’adsorption physique qui est le résultat de réactions impliquant de liaisons physico-
chimiques de faible énergie (échange d’ions et forces de Van der Waals),
- l’adsorption chimique qui concerne les liaisons chimiques de forte énergie (liaisons
covalentes).
Pour décrire les processus d’adsorptioddésorption, des hypothèses d’équilibre local
sont souvent utilisées :
- l ’hypothèse d’équilibre local (ZEA) : on considère que les processus d’adsorption sont
instantanés par rapport à la convection et la dispersion. La sorption et la désorption sont
décrites par des relations entre la concentration massique de soluté adsorbé S et la
concentration de soluté C de la solution d’équilibre. Ces relations sont les isothermes de
sorption-désorption. Les isothermes S(C) peuvent être de différents types (explicités ci-
après). Dans les conditions naturelles, les phénomènes de sorption ne sont réversibles que
lorsque le sol et la solution contiennent sufisamment de cations de remplacement pour
permettre la désorption. On peut tenir compte de ce phénomène d’hystérèse en utilisant des
coefficients différents pour les isothermes suivant qu’on est en sorption ou en désorption.
Cette approche conduit à la définition d’un coefficient Kd quand la relation entre les
concentrations dans les deux phases est linéaire.
Lorsqu’on considère que les processus d’absorptioddésorption ont une cinétique qui ne peut
être négligée, il faut prendre en compte la cinétique de sorptioddésorption.

Les réactions d’échange ioniaue

Ce processus s’ajoute à l’adsorption. L’existence de charges électriques non
neutralisées à la surface et / ou à l’intérieur des minéraux en est à l’origine. Des ions de
charges opposées viennent se fixer sur ces charges non neutralisées et créent ainsi une double
couche électrique. On distingue les échanges de cations des échanges d’anions suivant la
charge électrique des éléments en présence. Ainsi par exemple, les charges négatives du
réseau cristallin des argiles ne sont pas totalement compensées par les cations appartenant à ce
réseau. Des cations ((compensateurs)),interchangeables, peuvent être attirés et retenus à la
surface des minéraux argileux.

argileCa + 2K a argileK + Ca2’

Les réactions de solubilisation-DréciDitation

Ces réactions sont aussi importantes que les phénomènes d’adsorption ou d’échange
ionique. Elles se produisent dans l’eau des pores ou en surface des particules du sol. Ce
processus a pour action de précipiter les solutés dissous dans l’eau [Dupont 941.

A+ + X- <====== > A x
(forme ionique) (forme moléculaire solide)

L’évolution de l’équilibre réactionnel est fonction de nombreux paramètres

(concentrations respectives dans la phase aqueuse, composition du sol, température, pH,

-23-
 Transfert de 1’eau et des volluants dans les sols

pression,....). Ces réactions peuvent être modélisées à l’aide des constantes d’équilibre ou des
produits de solubilité que l’on peut trouver dans les bases de données thermodynamiques.

N.4.2. Les réactions microbioloni4ues

Dans les sols, tous les types de micro-organismes existent (champignons, algues,
bactéries, protozoaires...). Depuis longtemps on s’est aperçu qu’ils jouent un rôle fondamental
dans maintes transformations chimiques et mêmes physiques. Les approches modélisatrices
ont souvent négligé d’aborder ce sujet car les interconnections existants entre les micro-
organismes, le sol, les plantes et les composantes chimiques et physiques sont extrêmement
complexes et encore assez méconnues. La connaissance en ce domaine est très incomplète,
empirique ou complètement théorique. La prise en compte de l’action des micro-organismes
sur un polluant dans la modélisation est donc assez récente et liée à l’utilisation de plus en
plus courante de micro-organismes pour la dépollution des sols ou bioremédiation.
La biodégradation d’un composé organique en condition aérobie (02présent) ou
anaérobie (absence d’oxygène) peut être définie comme une réaction de réduction catalysée
biologiquement. Elle conduit fréquemment à la conversion d’un composé organique en un
produit inorganique. Une telle conversion est appelée minéralisation. Toutes les réactions
microbiennes conduisent au stade ultime, à un produit inorganique par l’intermédiaire d’une
minéralisation. Dans ce type de réaction, le carbone, l’azote, le phosphore, le souffre et les
autres éléments en présence, sont convertis en CO2 et autres formes inorganiques. Une partie
cependant reste sous forme organique et est incorporée dans les cellules de la population
active pour assurer le fonctionnement. Une autre partie peut apparaître comme un produit
typique du métabolisme de l’organisme.
L’action des micro-organismes est souvent perçue comme positive dans le sens où elle
transforme une substance active en un produit inactif. Mais cette transformation peut induire
la production de molécules toxiques. [Alexander, 19941.
L’évaluation de l’action, de l’efficacité des micro-organismes ainsi que la qualité et la
quantité de leur production paraît donc extrêmement importante. Cette évaluation passe
immanquablement par l’aspect modélisation. Plusieurs types de modèles sont utilisés afin
d’évaluer l’impact de l’activité microbienne et rendre compte des différents processus
chimiques qui interagissent entre eux . Ces modèles sont le plus souvent utilisés dans le cadre
de dépollution de sols in situ.
Les modèles de réaction biologique vont représenter la manière dont les organismes
transforment les polluants et ce, plus ou moins rapidement. Le taux de biodégradation dépend
de la somme de biomasse présente, de la concentration en polluant, du type d’accepteurs
d’électrons, des sources de nourriture disponibles, et de paramètres décrivant le taux de
transport des polluants suivant les organismes en présence et enfin, du rendement des
réactions enzymatiques en jeu.
Les limitations de ces modèles sont en général liées au caractère très spécifique de ceux-
ci. La validation s’effectue sur un site et les résultats disponibles sont difficilement
extrapolables d’un site à un autre. [ In situ bioremediation, 19931.
Il est évident que le développement des micro-organismes et l’activité qui en découle est
particulièrement présent dans les zones non saturée. En particulier, les phénomènes de
battement de nappe sont particulièrement propice à ces phénomènes.

-24-
 Transfertde l’eau et des polluants dans les sols

W.5. La notion de coefficient de répartition

Le transport de polluants dans un sol est le résultat de compétition entre les

phénomènes conservatifs qui tendent a véhiculer le polluant, et les phénomenes non
conservatifs qui tendent a l’immobiliser (adsorption, précipitation). Afin de rendre compte des
différents processus de sorption-désorption, le concept de Kd a été introduit. Le Kd ou
coefficent de répartition intervient dans toutes les réactions citées précédemment et regroupe
en quelque sorte les différentes interactions entre les polluants et le sol. Il exprime la
répartition de l’élément polluant considéré entre la phase liquide où il est transporté et la
phase solide où il est adsorbé. Ce coefficient a pour avantage d’introduire dans les fonctions
de transfert, un facteur de retard.

concentration totale dans la phase solide (par unité de masse) S

Kd =
concentration totale en solution (par unité de volume liquide) =

La répartition d’un élément ou d’un composé entre une fraction immobilisée et une
fiaction dissoute peut être exprimée suivant différentes relations mathématiques. Ces relations
mathématiques représentent les isothermes qui décrivent les équilibres thermodynamiques de
répartition de concentration entre la phase solide et la phase liquide On distingue
classiquement trois grandes classes d’équilibre dont les isothermes ont les formes suivantes :

Freundiich Langmuir i inear

c C C

Figure 2 : Isothermes d’adsorption les plus fréquents

Ainsi pour les métaux, l’isotherme d’adsorption est souvent de type Freundlich
(parabole) ou Langmuir (hyperbole). Pour les polluants organiques hydrophobes, 1’isotherme
d’adsorption est souvent linéaire, et par conséquent le Kd correspondant est indépendant de la
concentration totale du composé. Ii est généralement admis que pour les métaux, à des faibles
concentrations, le coefficient de répartition peut être considéré comme constant.

Il est possible d’estimer le & par l’expérimentation en comparant la vitesse de

transfert d’un élément par rapport à celle de l’eau. L’expression à laquelle on aboutit, est
considérée comme un facteur de retard qui peut être utilisé pour estimer la vitesse de
migration du polluant. Cette expression n’est valable que pour une adsorption linéaire,
réversible et instantanée.

-25-
 Transfertde 1’eau et des Dolluants dans les sols

Le facteur de retard est souvent exprimé dans les modèles déterministes par une
équation faisant intervenir la densité et la porosité du sol et le Kd de la substance. Le Kd pour
une fonction linéaire est calculé lui même par le Koc et la fraction organique du sol. Le
facteur de retard est déterminé par l’équation :

1+ (t* KCI) (cas de DEDALE 3D)

Il ne faut pas omettre d’autre part que l’adsorption peut rentrer en compétition avec les
réactions d’hydrolyse et de complexation - dissolution. Par exemple la migration du mercure à
travers une colonne de sable est fortement dépendante de la composition de l’eau.

Le coefficient de répartition permet de simplifier la modélisation des transferts de

polluants, et dans certains cas, cette notion permet d’obtenir des prédictions tout à fait
raisonnables. Cependant l’utilisation du & implique plusieurs contraintes à respecter si l’on
veut obtenir des résultats valides. D’une part la réaction d’adsorption, doit être suffisamment
rapide par rapport à l’écoulement de l’eau et réversible. Le & est soumis au comportement de
l’élément ou du composé que l’on considère (spéciations préférentielles, comportement
chimique,...). De plus, il dépend également de la stabilité de la chimie de l’eau environnante.
On peut noter aussi, l’influence de la présence de grosses molécules organiques dissoutes (ou
sous forme colloïdale) comme les substances humiques qui sont susceptibles de diminuer
l’adsorption des molécules organiques hydrophobes.

IV.6. Etablir l’équation de transport

La vitesse de l’écoulement du fluide est donnée par la loi de Darcy. Lorsque le sol est
non saturé, cette vitesse dépend en plus de la teneur en eau volumique 8.
Les équations de transfert de soluté doivent tenir compte de l’état hydrique du sol et de
toutes les possibilités de modification de la concentration. En tenant aussi compte de la
disparition ou de l’apparition de l’élément chimique, par application du principe de la
conservation de la masse, l’équation générale de transport isotherme d’une espèce chimique
est la suivante (transfert unidirectionnel) :

Ci est la concentration de l’espèce chimique dans la phase liquide.

Si est la concentration de l’espèce chimique dans la phase solide.
Di est le coefficient apparent de diffusion de l’espèce chimique, dépendant de la teneur en eau
volumique et du flux de Darcy donné par l’équation (3).
Le premier terme (a) de l’équation (IV.8) décrit le taux de passage de l’espèce (i) de la phase
liquide à la phase solide et vice versa; il inclue le processus de sorptioddésorption et
d’échange ionique.

-26-
 Transfert de l’eau et des polluants dans les sols

Le deuxième terme (b) est le taux d’accumulation de l’espèce (i) dans la phase liquide.
Le terme (c) décrit le transfert diffusif associant la difision moléculaire ou ionique à la
dispersion cinématique.
Le terme (d) représente le transfert convectif, l’espèce chimique (i) étant transporté par le flux
de liquide.
Le terme (e) prend en compte les sources et consommations de l’espèce (i).

L’utilisation du coefficient de partition & permet l’association des deux

concentrations Si et Ci (phase solide et phase liquide). L’utilisation de la relation entre Si et Ci
permet l’introduction dans l’équation de transfert de la fonction retard R(C). On peut alors
réécrire l’équation de transport sous la forme ( transport unidirectionnel dans un sol saturé):

avec R(C) = 1+ -K* si la relation entre s et c est linéaire soit : s = &C

n
yd est la densité sèche du sol, n la porosité.

Lorsque le transport a lieu dans un sol non saturé, les paramètres 8 (teneur en eau
volumique) et q (flux de Darcy, lié à la vitesse de pore) évoluent. On doit alors associer à
l’équation de transport de soluté l’équation de Richards en régime transitoire, soit :

at az

ou
avec
c(y) =-
ae est la capacité de rétention du sol; cette fonction rend compte de l’aptitude d’un
av
milieu poreux à libérer ou à emmagasiner de l’eau sous l’effet d’une variation de pression.
Elle correspond au coefficient d’emmagasinement d’une aquifëre lorsqu’on atteint la
saturation.
De(0) est la diffusivité de l’eau et w la succion. L’écriture de ces correspond aux deux
écritures de la loi de Darcy, soit à partir du gradient de succion (forme la plus fréquemment
utilisée), soit à partir du gradient de teneur en eau en introduisant la difisivité de l’eau dans
un sol.

Les équations de transport peuvent être établies pour un écoulement bidimensionnel et

tridimensionnel. Si le milieu est anisotrope, les paramètres de ces équations doivent être
évalués pour chacune des trois directions. La résolution d’un tel système est très lourde et très
compliquée. Une appréciation des différents phénomènes est nécessaire avant toute recherche
de solution d’un cas réel.

-27-
 Modéliser pour simuler

V. MODELISER POUR SIMULER

V.l. Qu’est-ce qu’un modèle ?

Comme il a été souligné précédemment, la compréhension des phénomènes de

transfert est complexe et la modélisation doit prendre en compte le fait que les processus
étudiés sont en général couplés les uns aux autres, qu’ils sont non linéaires et dans un milieu
polyphasé (eau-air-sol). Ce milieu est susceptible de se déformer et les concentrations en
fluide varient dans le temps et l’espace sous l’action conjuguée de variation d’energie
mécanique et thermique avec la possibilité de changements de phase et d’échanges physico-
chimiques entre les phases. Un modèle doit donc être considéré cornme une image simplifiée
d’une réalité complexe qui permet de simuler et d’extrapoler cette réalité.
On peut représenter grossièrement un modèle sous une forme générale telle que :

Y = F(X,A) + E

X est le vecteur des variables d’entrée et Y, le vecteur des variables de sorties. A

représente le vecteur des paramètres à prendre en compte et E, l’erreur commise en assimilant
l’objet de l’étude à son image modélisée mathématiquement. Tout l’art de la modélisation
réside donc dans l’appréciation de la validité des résultats en fonction des erreurs. Un modèle
doit être précis, robuste, adaptable et simple d’usage dans la mesure du possible.
Globalement on distingue trois grandes classes de modèles :
- lorsque F résulte de l’expérimentation, il s’agit de modèles empiriques. Cette appellation est
plutôt employée pour désigner le mode d’obtention de ceriaines lois (comme la Loi de Darcy,
tirée directement de l’expérimentation) utilisées dans les codes de calcul des modèles de
transfert,
- les modèles déterministes correspondent à une fonction F fondée sur une certaine
connaissance mathématique et physique des processus impliqués,
- si les variables d’entrée ou les paramètres sont considérés comme des fonctions aléatoires,
les modèles sont dit stochastiques,
- les modèles type (boîte noire» : ce type de modèles est basé sur la détermination de relations
entre les données d’entrée et la sortie après analyse du système et ceci, sans que la structure
interne de la boîte ne soit connue.

V.2. Les problèmes d’échelle de travail

La modélisation implique de travailler à des échelles différentes, et le passage d’une

échelle à une autre est en général très délicat. Si l’on veut faire une analyse fine des
phénomènes de migration de l’eau et des solutés, on doit distingue trois niveaux d’approche :
l’échelle microscopique, l’échelle du volume représentatif et l’échelle du terrain.
L’échelle microscopique prend en compte les phénomènes au niveau de la phase
aqueuse contenue dans les pores et à celui des interfaces du liquide avec les parois des pores.
On définit alors des paramètres locaux qui décrivent une quantité physique en un point, c’est-
à-dire un élément de volume infiniment petit. Les géométries de la porosité définissent les
domaines de continuité. Cette approche présenterait l’avantage de décrire les phénomènes

-28-
 Modéliser pour simuler

avec l’approche la plus rigoureuse mais elle n’est pas réaliste pour des milieux poreux aussi
complexes et hétérogènes que le sont des sols.
L’échelle spécifique des milieux poreux est très généralement une échelle plus
macroscopique. C’est à son niveau que le mélange polyphasique (eau air sol) peut-être
considéré comme un milieu continu. On définit alors un Element de Volume Représentatif
(EVR) [Gaudet, 19741, autour d’un point, choisi à l’intérieur d’un pore. Les grandeurs
(porosité, concentrations, pression, vitesses d’écoulement,. ..) peuvent alors être considérées
comme des variables continues. Le domaine d’étude est alors le milieu poreux dans les limites
où ces grandeurs peuvent être définies et continues. Il faut noter le domaine de recherche qui
vise à (< réconcilier N les deux approches précédantes: l’homogénéisation vise à trouver les
paramètres macroscopiques à partir de la description microscopique des phénomènes.
Le passage de l’échelle du laboratoire à celle du terrain pose évidemment de très gros
problèmes dans le cas des transferts de polluants dans les sols car les caractéristiques du
terrain ne peuvent être approchées que de manière très schématique ou incomplète à l’échelle
du laboratoire ou même du pilote. C’est certainement là l’une des limitations pour une
utilisation pertinente des modèles de simulation.

V.3. Les différentes approches de la modélisation

V.3.1. L approche déterministe

V.3.1.1.Objectifs et Iimites de I ’approchedéterministe

L’approche déterministe conduit à expliciter par des équations mathématiques

l’ensemble de tous les phénomènes qui régissent le transport des polluants, comme on l’a vu
au chapitre précédant. A chaque cas étudié correspond une hiérarchisation des phénomènes à
prendre en compte:
- nature et caractéristiques du milieu poreux,
- interactions des solutés avec le sol,
- conditions du scénario considéré...
Les lois de Darcy et/ou de Richards permettent théoriquement de calculer les
écoulements de la phase aqueuse et les diverses réactions bio-physico-chimiques peuvent être
traduites par des équations correspondantes.
Cette approche atteint cependant un certain nombre de limites et difficultés:
- pour modéliser le transfert de polluants dans un sol, les différentes connaissances à metire en
jeu proviennent de différents champs disciplinaires et les formalismes des équations
mathématiques qui en résultent ne sont pas toujours faciles à intégrer dans une approche
unifiée,
- le problème d’échelle de travail est crucial pour l’approche déterministe car il est très
difficile, voire impossible, de garantir l’identité des phénomènes au laboratoire et sur le
terrain,
- la connaissance précise des nombreux paramètres à prendre en compte est souvent très
lourde ou même impossible,
- la résolution numérique d’un modèle déterministe détaillé peut conduire ell-même à de
grosses difficultés, hors de proportion avec la qualité espérée des résultats de simulation,
- en général, une description déterministe détaillée menée sur un cas d’étude perd la
possibilité d’être généralisable.

-29-
 Modéliser pour simuler

Il ne faut cependant pas condamner sans appel l’approche déterministe. Au laboratoire,

elle permet, en tout cas, de mettre en évidence et de quantifier les phénomènes qui
déterminent le fonctionnement du système étudié, permettant ainsi d’étayer une approche plus
pragmatique. Elle constitue donc un excellent outil de recherche pour étudier en détail les
interactions entre l’eau, le sol et le soluté. Néanmoins, la détermination des paramètres, la
validation des modèles ainsi que l’utilisation de ceux-çi en condition naturelle, sont autant de
limites à considérer.

V.3.I .2.Les modèles déterministes

Ces modèles sont fondés sur la connaissance plus ou moins établie des processus mis
en jeu lors de la migration des polluants dans les sols. On peut distinguer les modèles
essentiellement fondés sur la notion de vitesse et de flux de matière et qui utilise les lois de
Darcy, de Fick. Ces modèles qui sont utilisés la plupart du temps, sont dit déterministes
mécanistes.
Ces outils déterministes mécanistes tiennent compte, dans l’absolu, des phénomènes
physiques suivants:
- les phénomènes de transport représentés sont la convection et la diffusion,
- les phénomènes physico-chimiques pris en compte sont l’adsorption du soluté sur la
matrice poreuse, la volatilisation dans la fraction gazeuse et la dégradation des composés par
les micro-organismes contenus dans le sol [Guinot, 19921.
En parallèle les modèles fonctionnels sont de type capacitifs : le sol est assimilé à un
réservoir. Ces modèles ne s’attachent pas au calcul du transport dans un seul milieu physique
particulier. Ils sont utilisés dans le cas de simulation à l’échelle du bassin versant. Ces
approches effectuent également une modélisation de l’hydrologie et du transport fondée sur
des hypothèses ou des descriptions simplificatrices des processus physiques.
Certains modèles négligent le transport par diffusion et dispersion, seuls sont pris en
compte les transports par convection. Ce type de modèle a été utilisé avec un certain succès
pour simuler le devenir de concentrations appliquées à la surface d’un sol, mais ils présentent
de sérieuses limitations lorsque le sol est hétérogène ou contient déjà des polluants.
D’autres modèles adoptent une approche chromatographique pour simuler le lessivage
des sols en considérant une (( discrétisation D du sol en compartiments disposés en cascade.
Cette approche apparaît pertinente dans le cas de transport de plusieurs espèces chimiques ou
de substances fortement interactives avec la phase solide. Ces modèles ont l’avantage de
présenter une structure mathématique simple et ils n’utilisent qu’un nombre réduit de
paramètres et de données d’entrée ce qui est un gros avantage. L’utilisation de ce type
d’approche est là aussi limitée lorsque le sol présente de fortes hétérogénéités ou un
géomorphisme complexe.
Les approches déterministes sont utilisées dans la plupart des logiciels de calcul que
l’on trouve sur le marché.
Les logiciels MARTHE et TAFFETAS développés au BRGM utilisent la méthode aux
volumes finis pour la résolution du modèle. Ils résolvent simultanément les équations de
conservation de la masse d’eau et de la masse de soluté. Une loi de variation de la densité est
prise en compte lors de la discrétisation de l’équation de conservation de la masse de soluté.
Les logiciels SIMUSCOP de lIFP, RESCOPP de l’EN1 (Italie), BURGEAP et 1’IFP
resolvent les équations de transport en lD, 2D et 3D en utilisant la méthode des éléments finis
ou des volumes finis. La géochimie est prise en compte par introduction du facteur retard. Ces

-30-
 Modéliser pour simuler

logiciels se sont pricipalement intéressés aux hydrocarbures et paramètres hydrodynamiques

nécessaires sont donnés par des tables (bases de données).
Le logiciel DEDALE 3-D utilise des équations déterministes plus détaillées. La
variation de la densité est prise en compte aussi bien dans l’écriture des équations de
conservation de la masse d’eau et de la masse de soluté.
Ces logiciels permettent d’établir la répartition spatio-temporelle de la charge
hydraulique ou de la succion et celle de la concentration des éléments chimiques. Le nombre
de paramètres à déterminer pour leur utilisation sont d’autant plus nombreux que le modèle
mise en oeuvre est complexe.

V.3.2. L ’amroche stochastique

Les propriétés physico-chimiques d’un sol présente de grandes fluctuations dans

l’espace et au cours du temps. Les paramètres utilisés pour la description du transport de
polluant n’échappent pas à ces fluctuations. Ces paramètres étant obtenus après
échantillonnage ou après mesures in situ en nombre limité, le problème de la représentativité
se pose toujours. Pour limiter ces inconvenients, on envisage une démarche statistique.
A partir de l’ensemble des mesures disponibles, on effectue une analyse statistique des
propriétés de l’échantillon, sa moyenne, sa variance et sa loi de distrbution. Cette approche à
été abondamment décrite dans la littérature (voir Wanick et Nielsen, Mac Bratney et Webster
1983 et Vauclin 1983). Elle conduit à considérer les phénomènes de transfert comme des
processus erratiques susceptibles d’être quantifiés par des modèles stochastiques.
Les propriétés hydrodynamiques sont réprésentées par des fonctions aléatoires. qu’on
introduit dans les modèles mécanistes. La solution d’un problème de transport est présentés
sous forme de probabilité ou en terme de risques. Par exemple, si on admet la loi de Darcy, la
conductivité hydraulique K(y) du milieu naturel sera considérée comme un processus
stochastique. De cette loi et de l’équation de Darcy résultent le fait que la vitesse de pore, du
cofficient de dispersion cinématique sont des processus stochastiques.
La résolution des équations de base après introduction des différents processus
stochastiques, doit donner en théorie, la distribution spatio-temporelle de la concentration en
soluté pour toute les réalisations possibles de K
L’équation stochastique de transport a été très peu étudiée. On peut citer le modèle de
Bresler et Dagan (1979, 1980 et 1983) qui reste à l’heure actuelle le plus représentatif. Les
résultats actuellement disponibles aussi bien pour l’hydrodynamique que pour le transport,
démontre la non existence de propriétés effectives capable de caractériser un milieu
“ “

homogène équivalent à un milieu réel pour lequel la solution du problème hydro-dispersif

dans un sol non saturé serait identique à la moyenne des solutions du problème dans le milieu
réel hétérogène (Vauclin 1993).
Cette approche des paramètres hydrodynamiques par les fonctions aléatoires suppose
la formation géologique étudiée bien définie (grès, argile, roche cristaline). La description de
la formation elle-même peut être fortement erratique. Une description stochastique du sol
s’impose. On définit alors les propriétés moyennes du milieu poreux. Cette description est
liée à l’échelle utilisée (microscopique ou macroscopique). Le passage de l’échelle
microscopique à l’échelle macroscopique nécessite un outil mathématique (de Marsily 1992).
Une fois les lois statistiques des propriétés du milieu poreux établie, elles sont alors étendues
à la loi de Darcy et aux équations de conservation de la masse pour une description
macroscopique de l’écoulement. Les champs de vitesses et de pression sont représentés par
leur moyenne d’après la loi de probabilité utilisée.

-3 1-
 Modéliser pour simuler

Une dernière approche stochastique de l’écoulement est la représentation du milieu

naturel par des milieux aléatoires. Des techniques de représentation de milieux poreux ont été
développées pour tirer meilleur parti des résultats des campagnes de reconnaissance
géologique ou pédologique et rendre compte de la variabilité des faciès. Il s’agit de trouver un
modèle de site d’après les informations données par les sondages. Des exemples d’application
telles que la dépollution, l’estimation de la ressource en eau d’un aquifêre demandent que
soient connus des modèles de la géomètrie du milieu. La technique des milieux aléatoires
“ “

propose de décrire le milieu par un nombre de “blocs élémentaires” pouvant atteindre

plusieurs millions. Ces nouvelles méthodes sont très liées à l’utilisation de méthodes
numériques. Les méthodes développées dépendent du domaine d’application.
La méthode booléenne a été développée dans le domaine pétrolier. Elle consiste à
générer dans l’espace, de façon aléatoire, des “objets” dont la position et les dimensions sont
tirées au hasard, tandis que la forme est à priori connue. Il suffit alors de définir les paramètres
hydrodynamique dans l’objet et dans la matrice pour calculer un écoluement. Le choix du type
d’objets élémentaires est infini. Plusieurs types d’objets peuvent être superposés, représentant
plusieurs faciès (sables grossiers, sables fins, limons, argiles, par exemple).
La demière méthode de description de la géométrie d’un milieu poreux est la
simulation de processus sédimentaires qui génèrent des situations géomorphiques d’un site.
La génération des faciès est faite non pas en simulant les processus statistiques ayant une
ressemblance avec les distributions observées, mais comme le résultat d’un processus
physique de sédimentation. On va donc engendrer des successions de crues dans les rivières
ou delta, en simulant le débit du fleuve, sa charge en sédiments et la sédimentation de cette
charge, et éventuellement la reprise partielle ou totale de cette charge par l’érosion. On suit
ainsi le mécanisme de formation des dépôts au cours du temps, en fonction du climat, de la
pente du sol, de l’érosion et du transport solide. Cette méthode n’est pas conditionnée par les
observations; il n’existe aucun moyen de forcer la simulation à respecter les faciès observés au
forage. Contrairement aux deux méthodes précédentes, une seule simulation est faite pour la
description du milieu poreux. Vue l’aspect aléatoire que comporte la succession des crues, il
serait possible d’engendrer plusieurs simulations mais la lourdeur des calculs l’interdit pour
l’instant.

K3.3. Les fonctions de transfert

Il s’agit, de manière schématique, d’établir une relation stochastique entre un signal et

sa transformée par le milieu traversé. Par exemple, si l’on désigne par Qe le débit d’entrée et
Qs le débit de sortie, cette relation peut se mettre sous la forme :

La fonction g(t1t’) est la densité probabilité des temps de séjour. Elle décrit la
probabilité qu’aura une molécule entrant dans le volume de contrôle ( par exemple la surface
du sol ) au temps t’ de sortie (par exemple à une profondeur définie) entre les instants t et t’.
Elle représente l’effet de l’ensemble des processus et de la nature des conditions initiales et
aux limites sur le temps de séjour d’une substance dans le sol.
Ce type d’approche a été utilisé avec succès pour simuler des expériences in situ
réalisées avec différentes substances (Br, C1, Nitrates, pesticides, Jury et al. 1986; White et al.
1986). Cette approche peut constituer un modèle de gestion de site interessant dans la mesure

-32-
 Modéliser pour simuler

où il ne requiert il ne requiert qu’un nombre réduit de données. Ces types de modèles très
souvent utilisés en hydrologie, nécessitent un étalonnage spécifique aux conditions
d’utilisation.
Toute modification du milieu, des contions aux limites exige la recherche d’une
nouvelle fonction. Ce type d’approche a conduit à des modèles de types «« boîte noire ».

V.3.4. Les modèles t w e «boîte noire»

Dans l’approche «« boîte noire », le modélisateur détermine les relations existant entre
les entrées et les sorties d’un système dynamique
Il faut signaler une technique d’identification assez récente des systèmes basée sur
l’utilisation des réseaux de neurones. En utilisant des séries d’entrées et de sorties observées,
le réseau de neurones peut déterminer une «« pondération >> entre les données d’entrées qui
puisse foumir la ou les sorties. Cette phase de calage est la phase d’apprentissage. La relation
qui existe entre les entrées et la sortie se trouve automatiquement calée par le modèle sans que
le modélisateur n’ait à faire aucun choix. La modélisation est réduite à sa plus simple
expression : il n’y a pas d’équations phénoménologiques à écrire, pas de géométrie, pas de
paramètres ni de conditions aux limites et pas de vérification d’un algorithme. Après le calage
(apprentissage du réseau), on peut utiliser le modèle en prévision et fournir la sortie pour un
nouveau jeu de données d’entrée. Il n’y a pas à notre connaissance d’application de cette
technique à la modélisation des transferts de polluants dans les sols. Une limitation évidente
est le grand nombre de données expérimentales nécessaires pour le calage du modèle.

V.4. Les modèles en hvdrologie

Les modèles en hydrologie sont utilisés afin de prédire l’évolution des débits à
l’exutoire d’un bassin versant, suivant différents paramètres (pluie, évapo-transpiration, état
hydrique du bassin versant et caractéristiques de celui-ci...). L’hydrologie s’attachant à des
phénomènes aléatoires, les approches s’appuient toutes sur une interprétation ou une
modélisation statistique des données hydrologiques disponibles, c’est-à-dire sur des mesures
de pluie et de débits. Les catégories de modèles sont identiques à celle inventoriées
précédement pour les modèles de transfert. Soit on progresse avec une approche déterministe
basée sur les connaissances physiques des phénomènes, soit avec une approche statistique ou
stochastique, soit sous la forme des modèles «réservoirs»où les bassins versants sont abordés
comme une suite de réservoirs.
On trouve classiquement des modèles déterministes couplés avec des modèles
stochastiques qui utilisent les équations avec des fonctions aléatoires à la place de paramètres
fixes mesurés sur le terrain ou par l’expérimentation. Cette approche courante en hydrologie
l’est moins en ce qui concerne les transferts en milieux poreux.
Dans les modèles déterministes hydrologiques, on retrouve le SHE (Système
Hydrologique Européen) qui est utilisé à l’intérieur de certains code de calcul pour la
détermination de la convection-dispersion de polluants dissous, ou pour la modélisation de
processus géochimiques (exemple de MIKE SHE).
Dans les modèles statistiques purs, on va avoir des modèles très spécifiques tels que
les modèles AR (auto regressif), ARX (AR exogène), ou ARMB( AR à moyenne mobile)...
Les modèles réservoirs sont parmi les plus utilisés en hydrologie, et suivant le type de

-33-
 Modéliser DoUr simuler

résolution, on va distinguer les modèles réservoirs linéaires, modèles de Muskingam, modèle

réservoir non linéaire, ou cascade de Nash. Ces modèles sont basés sur l’équation de
continuité suivante :

--
dS(t) - I(t) - Q(t)
dt

ainsi que sur les lois de comportement du bassin versant reliant le stock d’eau S à l’infiltration
1et au débit Q :

Si on s’attache plus en détail au terme d’infiltration, on peut faire appel aux différents
modèles d’infiltration, comme le modèle de Gree-Ampt (que l’on retrouve souvent dans les
codes de calculs pour les transferts de polluants dans les sols), de Philips ou de Horton. Ces
modèles sont basés soit sur Darcy, soit sur Richards suivant la zone qu’ils considèrent.

De la même manière que pour l’infiltration, les hydrologues ont également établi des
modèles pour l’évapo-transpiration qui prennent en compte les paramètres atmosphériques
(vitesse du vent, humidité de l’air, température, rayonnement solaire..) alliés aux parmètres
caractéristiques des sols.
L’approche hydrologique va etre utilisée lorsque l’on va vouloir modélisé l’impact d’une forte
pluie sur le relargage de polluants dans les sols.

V.5. Les modèles dochimiaues

La géochimie des sols pose encore beaucoup de problèmes quant à la modélisation.

Les logiciels détaillés n’en sont encore qu’à leur début et sont peu nombreux. On considère
souvent la chimie sous l’angle simplificateur du Kd qui joue un rôle de boîte noire. Les
modèles de géochimie sont des logiciels de calculs thermodynamiques des équilibres en
solution en contact avec les roches.
Comme exemple, on peut citer le modèle Geochimist’s Workbench, developpé par
Bethke depuis 1992, qui aborde le problème des équilibres chimiques de manière très
complète. Ce type de modèle utilise plusieurs codes de calculs qui déterminent les équilibres
thermodynamiques et les chemins réactionnels possibles pour les systèmes fluide + minéraux
+ gaz rencontrés en géochimie.
La minéralogie des formations du site ainsi que les compositions et caractères physico-
chimiques des eaux sont utilisés pour l’évaluation des réactions possibles, des spéciations des
solutions et des chemins réactionnels. La modélisation permet, à partir des interactions fluide/
roche simulées, de retrouver des solutions ayant des caractéristiques proches de celles des
eaux analysées.
La confrontation entre la composition des solutions simulées et des solutions analysées
aide à la compréhension de l’évolution de la solution polluante dans l’encaissant. La
modélisation doit permettre d’évaluer la réaction de l’encaissant vis à vis des éléments
impliqués, les types de réactions, le rôle des différentes formations, les éléments mobiles et
immobiles...
On peut citer quelques cas d’étude: la contamination de sols et des eaux sous-terraines
par des polluants minéraux, l’évaluation du devenir de produits industriels fluorés, le draînage

-34-
 Modéliser vour simuler

acide de rejets miniers sulfurés, l’injection de gaz dans un aquifère profond et l’impact
hydrochimique de la remise en eau d’une mine métallique.

V.6. Les différentes approches rencontrées au cours de la

modélisation

Lors de la modélisation il est indispensable de réfléchir au préalable sur (par

exemple) :
- l’approche dimensionnelle : 1-D, 2-D ou 3-D,
- le type de résolution spatiale de la zone à étudier : différences finies, éléments fuis,
volumes finis,
- le schéma de résolution en régime transitoire : implicite, explicite, implicite-
explicite.
Dans les différents outils de simulation de transfert de polluants dans les sols,
présentés par la suite, ces différentes approches sont abordées de manières différentes. Se pose
alors le problème du choix ou du bien fondé de ces différentes approches pour obtenir des
simulations les plus proches de la réalité. Nous ne donnerons ci dessous que les éléments de
réflexion pour l’approche dimensionnelle.

Les écoulements souterrains sont de façon générale de nature tridimensionnelle. Dans

certains cas, une hypothèse de «bidimensionnalité» peut cependant être faite : c’est
l’hypothèse de Dupuit. Dans cette hypothèse, pour de faibles pentes de la surface libre, les
lignes de courant sont considérées comme horizontales et sont associées à des vitesses
proportionnelles à la pente de surface libre et indépendantes de la profondeur (ce qui implique
que les lignes équipotentielles sont verticales). Les modèles 2-D s’appuient sur cette
hypothèse et considèrent en un point de l’espace (x,y) une vitesse unique de l’écoulement
souterrain, quelque soit la profondeur.

On peut rencontrer trois types d’approches :

1- approche 2-D ZS

Figure 3 :approche 2-DZS

-35-
 Modéliser pour simuler

2- approche 3-D Z N S / ZS

Y
P

Figure 4 :approche 3-D Z N S / Z S

3- approche 3-D complète (DEDALE 3-D)

a \
e

Figure 5 :approche 3-D complète (DEDALE 3-D)

Il est important de noter que le devenir d’un polluant est complètement différent suivant les
approches considérées, comme le montrent les deux figures suivantes.

1- Une approche bidimensionnelle implique d’appliquer le polluant non pas à la

surface du sol, mais directement dans l’aquifère et sur toute l’épaisseur de celui-ci, de façon
homogène. Par conséquent, même si l’hypothèse de Dupuit est valable, l’utilisation de ce type
de modèle pour des problèmes de transfert de polluants s’avère délicate.

Zone non sururée

Figure 6 : approche bidimensionnelle

-36-
 Modéliser pour simuler

2- Dans une approche tridimensionnelle, le polluant peut effectivement être appliquée

à la Surface du sol. Sa propagation se fait dans la zone non saturée puis ensuite dans la partie
supérieure de la zone saturée. Sa concentration selon une verticale est hétérogène. Dans le cas
d’une approche complète, a la différence d’une approche couplées Z N S / ZS, le polluant peut
migrer latéralement au sein de la zone non saturée.

Figure 7 :approche tridimensionnelle

V.7. Conclusions

Différentes approches modélisatrices sont utilisées pour décrire les transferts de

polluants dans les sols. L’enjeu principal est la prise en compte réaliste de la complexité des
phénomènes réels.
Une première difficulté réside dans le fait que la description macroscopique du
transfert de polluants dans un sol doit être menée sur de grandes échelles d’espace et de temps
alors que la connaissance précise des mécanismes mis en jeu n’est possible qu’à l’échelle du
laboratoire. De plus, les connaissances proviennent de champs disciplinaires très différents.
Un modèle global est a priori difficile à construire, ne serait-ce qu’à cause des formalismes
provenant des différentes approches.

L’approche déterministe, qui a pour ambition de coller au plus près des mécanismes
physico-chimiques réels, atteint ses limites lorsque la complexité des phénomènes réels
empêche de garantir la pertinence des choix nécessairement simplificateurs du modèle ou la
fiabilité des données nécessaires.

L’approche stochastique, séduisante a priori pour une approche des situations

compexes et les possibilités de changement d’échelles, ne s’impose qu’au prix d’une approche
très rigoureuse des données statistiques manipulées et de validations expérimentales poussées.

-37-
 Synthèse de la demande

VI. SYNTHESE DE LA DEMANDE

VI.1. Motivations de la demande

Une analyse de risques dus à la pollution d’un site passe par l’évaluation de la pollution
actuelle et de son devenir. Cette évaluation peut être menée à l’aide d’outils de modélisation
et de prévision de la migration de polluants, et plus particulièrement la migration de ceux-ci
dans la zone non saturée en direction de la nappe.
Aujourd’hui, les membres de FEC0.R.D. s’intéressent à la problématique de
l’évaluation du risque “sols pollués” essentiellement pour trois raisons :

+ La nouvelle réglementation :
Suite à la circulaire de Barnier (circulaire du 3 décembre 1993, Annexe 2) des
consignes sur l’étude des risques pour les sols pollués arrivent aux DFüRE. Cette circulaire
annonce aux préfets la mise en place d’une véritable politique de traitement des sites et sols
pollués.
Les études de risques impliquent :
* la caractérisation de la source de pollution,
* la caractérisation de la cible (attention toute particulière à l’eau potable),
* la caractérisation du transfert des polluants jusqu’à la nappe. C’est ici qu’intervient
l’utilité des modèles de migration de polluants dans les sols et notamment en zone non
saturée.
Une grille d’évaluation ainsi que les niveaux d’exigence de décontamination sont
joints à la circulaire.

+ La cessation et l’achat d’activités :

Les arrêtés d’autorisation prévoient les conditions de remise en état des lieux à
l’expiration d’une exploitation. L’article 34, alinéa 2 du décret du 21 septembre 1977
(Annexe 3) impose à l’exploitant d’informer le préfet dans le mois de la cessation d’activité
sous peine d’une amende contraventionnelle. Il doit remettre le site dans un état tel qu ’il n ’y
ait plus d’inconvénients ou de danger pour la santé ou 1 ’environnement.A défaut du respect
de cette obligation générale le préfet pourra mettre en oeuvre les sanctions administratives
prévues à l’article 23 de la loi du 19 juillet 1976 (Annexe 4). Dans le cas d’un rachat,
l’industriel sera responsable du terrain acheté, et par conséquent de la dépollution nécessaire si
elle n’a pas été faite. Le coût de rachat peut alors devenir exorbitant. Des litiges et des
poursuites judiciaires apparaissent déjà entre des acheteurs de terrains qui nécessitaient une
dépollution et qui n’en ont pas bénéficié pour cause (d’omissions volontaires ou non).
Dans le cas d’une cessation d’activité ou d’un rachat, la problématique de l’industriel
est simple : d partir de quand doit-il dépolluer et jusqu ’où ? Si le site est pollué, l’industriel
doit savoir si cette pollution va évoluer ou non (migration des polluants), et par conséquent,
quelle est l’urgence d’une mise en place d’une dépollution.

-38-
 Synthèse de la demande

+ Répondre à des problèmes internes de pollution :

Les industriels, suivant leurs activités présentes ou passées, peuvent avoir la
responsabilité de sites pollués. Ces sites sont connus, répertoriés et attendent d’être traités.
Mais lorsque l’on est en présence de plusieurs sites, il devient alors nécessaire d’avoir un outil
d’évaluation afin de hiérarchiser les sites à dépolluer. La modélisation de l’évolution de la
pollution dans le temps et dans l’espace peut permettre de prévoir quels sont les sites qui vont
demander un traitement rapide et ceux qui peuvent rester en attente (combien de temps ?) sans
conséquences néfastes pour l’environnement.

VI.2. La demande générale

v1.2.1.Les principaux besoins

La demande générale des membres de RJ3.CO.R.D. dans ce domaine est relativement

unanime. Les besoins s’orientent essentiellement vers des outils simples Ù utiliser, et donnant
la meilleure prévision possible avec une économie de données. Un reproche souvent rencontré
à propos des outils existants concerne leur lourdeur d’utilisation et la nécessité de posséder de
nombreuses données théoriques et surtout expérimentales (sur les polluants : répartition du
champ de concentrations, solubilités, coefficients de partage, constantes de Henry, pressions
de vapeur..., et sur les sols : perméabilité, porosité, teneur en carbone organique,
température,...) Or l’acquisition de ces données est très coûteuse, de plus elles sont
difficilement représentatives de la réalité complexe de sols en général très hétérogènes.

Les principaux besoins et les principales questions relatives à ces outils sont résumés
ci-dessous :
* étude des interactions polluants/sols, élaboration d’une base de données,
* problème de la connaissance des données physico-chimiques des polluants (solubilité,
coefficient de partage...), d’où le besoin d’un carnet d’adresse de laboratoires pouvant
effectuer ces analyses,
* demande d’outils simples et utilisables dans des cas concrets, et/ou domaine d’application
plus large,
* connaissance du domaine d’application des différents modèles : 2D, 3D, zone saturée, zone
non saturée, analytique, numérique...
* quels sont les moyens nécessaires pour utiliser ces outils de simulation : un béotien équipé
d’un micro ou alors une station de calcul avec une équipe de spécialistes ?
* problème de la “gourmandise” des logiciels actuels en données (implication sur les coûts),
* connaissance des paramètres nécessaires pour faire tourner le modèle avec une bonne
précision, et implication des données d’entrée sur les résultats de la simulation,
* compréhension des phénomènes,
* nécessité d’outils solidement validés sur des cas réels,
* nécessité de méthodologies d’acquisition de données sur le terrain qui soient pertinentes.
De manière générale, les demandes convergent vers la nécessité d’un outil fiable
d’aide à la décision, permettant, avec une économie de moyens mis en oeuvre, de caractériser
l’évolution future des polluants présents dans le sol.

-39-
 Synthèse de la demande

V1.2.2. Comutes-rendus d ’entretiensconcernant la demande.

Avertissement :
Ne figurent dans ce paragraphe que les éléments jugés significatifs des entretiens qui
ont été menés avec les représentants des membres de RE.CO.R.D. Ces comptes-rendus
pourront être complétés si nécessaire à la demande des tuteurs qui suivent l’étude.

ELF
Contact : M. PEREZ

Une préoccupation très actuelle est l’étude de risques concernant les sites pollués. Les
outils méthodologiques sont insuffisants, surtout en ce qui concerne la mobilité des polluants.
Pour évaluer les conséquences de la pollution d’un sol, un maillon manquant reste la
modélisation du terme de transfert des polluants de la source d’émission à la cible. Il faudrait
dépasser une approche qui reste relativement arbitraire pour caractériser la dangerosité des
polluants relativement à leur mobilité. En effet, la responsabilité des entreprises est engagée et
il faut un outil d’évaluation agréé par les différents acteurs.
Une approche critique des modèles existants est nécessaire pour déterminer leur
domaines d’application et la validité de leur pouvoir prédictif.
ELF s’est déjà engagé sur plusieurs programmes de recherche qui concernent les sols
pollués, notamment avec 1’IFP (projet EUREKA/RESCOPP).

GDF
Contact : M. VERCAMER

GDF est propriétaire de nombreux sites historiquement contaminés. Le problème

réside dans l’hétérogénéité extrêmement importante des sites. L’un des problèmes essentiels
va être de hiérarchiser la dépollution des sites. Au vu du nombre de sites, la dépollution va
être répartie au cours du temps. Mais dans quel ordre?
Les polluants sont souvent des hydrocarbures lourds résultant des anciens procédés de
production de gaz, mais les produits ont évolué au cours du temps et les sols ont connu des
remaniements. La modélisation des mécanismes de transferts est très complexe et il serait bon
d’avoir une approche d’expert pour avoir les éléments d’une bonne appréhension des
principaux phénomènes.
La demande est ciblée sur deux points essentiels :
+ * nécessité de mettre au point un inventaire des principaux couples polluants/sols afin
de pouvoir établir une méthodologie de travail pour un gain de temps et un meilleur
traitement.
+ * besoin de faire des évaluations de sols plus pertinentes afin de mieux gérer les sols
pollués, impliquant peut être une révision de la législation quand à la destination des sols
dépollués (décharge, incinération). La modélisation sert dans ce cas à prédire la mobilsation
ou non des polluants et d’en prédire l’évolution future.

-40-
 Synthèse de la demande

EDF
Contact : Melle DUMONT

Il n’y a pas de recherche sur la migration de polluants dans les sols chez EDF (sauf
pour le domaine particulier du stockage à long terme des déchets nucléaires pour lequel des
études sont engagées par I’ANDRA et le CEA).
EDF souhaiterait connaître l’offre commerciale existante, tant du point de vue des
logiciels que des laboratoires et bureaux d’études experts dans le domaine et capables d’offrir
des prestations. Un domaine intéressant EDF est le cas des PCB.

GEMMESmREDI
Contact : Mrs. HILERMANN, STRAUSS (Société GEMMES) et DURECU (Société TREDI)
Destruction des déchets industriels (TREDI)
Réhabilitation de sites (GEMMES)

Les techniques de réhabilitation de sites utilisées par la société GEMMES sont

diverses :
+ désorption par traitement thermique,
+ traitement par lavage,
+ traitement biologique (en cours de développement),
+ traitement par CO2 supercritique (en R&D).
La société GEMMES réalise des études prospectives sur les sites pollués. et possède des outils
de modélisation de migration de polluants dans les sols tel que le logiciel SUTRA (1984,
révisé en 1989) qu’elle utilise notamment pour la simulation de la migration des sels dans les
mines de potasse. Le principal reproche envers ce logiciel concerne sa lourdeur d’utilisation.
GEMMES recherche un logiciel plus simple d’utilisation, plus performant et avec un domaine
d’application plus large.

RHONE-POULENC
Contact : M. LEMAREC

La problématique “sols pollués” est essentielle pour 2 raisons :

+ - la réglementation nouvelle : directives européennes, Ministère de l’Environnement.. .
+ - cessation et achat d’activités.
RP a déjà effectué récemment un état de l’art sur les logiciels disponibles.
La préoccupation actuelle de Rp est de comprendre ce qu’il y a dans la “%boîtenoire”
de la modélisation. Pour cela RP mène des actions de recherche avec d’autres laboratoires et
les thèmes de travail sont les suivants:
+ - comportement des polluants dans les sols,
+ - modélisation de la migration des polluants :modèle d’écoulement et de transport de
polluants,
+ - analyse des risques avec une composante écotoxicologique.
L’intérêt est plutôt orienté vers les modèles numériques plutôt qu’analytiques. Les
modèles stochastiques ne semblent encore pas assez efficaces et fiables.

-41-
 Svnthèse de la demande

A l’heure actuelle, RP utilise plusieurs logiciels existant sur le marché :

+ DEDALE 3D (LHF - SOGREAH)
+ MARTHE 3D (ANTEA - BRGM)
+ MODFLOW (écoulement) et MT3D (transfert de masse)

SOCOTEC
Contact : M. MOREL

SOCOTEC Environnement est un cabinet d’études effectuant notamment le diagnostic

de sites pollués (historique, caractérisation géologique et hydrogéologique, sondages).
Le diagnostic donne un état du site (à un moment donné) qui permet une “analyse des
risques” puis d’arriver à une décision (dépollution ou non) mais le problème rencontré avec
les logiciels de calculs est la gourmandise en données qui exigent de nombreuses
investigations SUT les sites, et impliquent donc une augmentation sensible du coût pour le
bureau d’études qui utilise le logiciel ainsi que pour le client, sans garantir une prestation
meilleure que celle d’un expert.

SOLVAY
Contact : M. MICHIELS

Les modèles déjà utilisés sont MARTHE et MODFLOW (hydrodynamique), AQUA

3D (hydrodynamique et également transport des polluants en zone saturée).
Les polluants rencontrés chez SOLVAY :
+ organochlorés : dichlorométhane, di et trichloroéthylène ...,
+ pesticides : hexachlorocyclohexane, lindane et dérivés, pentachlorophénol ...,
+ hydrocarbures et solvants (cétones ...),
+ métaux lourds (métalliques) : mercure, chrome...
La demande de SOLVAY peut se résumer en quelques points:
+ tester AQUA 3D en zone saturée (à court terme),
+ une étude bibliographique en ZNS des différentes configurations selon le type de
polluants et le type de sols (déterminer la capacité de stockage de la ZNS),
+ étude de quelques polluants en Z N S (les plus sensibles comme par exemple le
mercure),
+ améliorer les techniques d’échantillonnage.

-42-
 Analyse de l’ofie

VII. ANALYSE DE L’OFFRE

VII.1.Les modèles les plus utilisés et leur champ d’application

Beaucoup des logiciels de calcul que l’on trouve maintenant sur le marché proviennent
du marché américain dont le fonctionnement permet de consulter gratuitement tous les
modèles réalisés dans les instituts d’état comme les Universités ou 1’EPA. Le monde des
modélisateurs fiançais s’est donc, dans beaucoup de cas, contenté d’améliorer les interfaces
qui sont par contre, rudimementaires et rendent les logiciels extrêmement compliqués à
l’utilisation à l’origine. Ces modèles sont tous basés sur l’approche déterministe et se
différencient par leur domaine d’application et la complexité des phénomènes pris en compte.

VII.1.1.ûuelques modèles

MARTHE
Modélisation d’Aquifère par un maillage Rectangulaire en régime Transitoire pour le calcul
Hydrodynamique des Ecoulements

Conçu et réalisée au BRGM @. Thiéry) pour la modélisation hydrodynamique et

hydrodispersive des écoulements souterrains en milieu poreux, le logiciel MARTHE
comprend :
- un module de calcul hydrodynamique et hydrodispersif (MARTHE proprement dit),
- quatre modules de préparation, gestion et représentation graphique des données et
résultats de modélisation
MARTHE permet de résoudre les problèmes d’écoulements rencontrés dans différents
domaines :
+ gestion des aquifères,
+ géniecivil,
+ environnement et plus particulièrement modélisation de la migration d’un polluant
dans une nappe,
+ stockage souterrain de déchets industriels,
+ exploitation minière.
MARTHE a été étendu pour lui permettre, en plus des écoulements classiques en aquifères, de
modéliser les écoulements complexes en zone non saturée et à densité variable.
Pour modéliser les écoulements en zone non saturée, on fait les hypothèses suivantes :
- écoulement isotherme,
- sol non gonflant (la porosité est constante au cours du temps),
- fluide incompressible,
- perméabilité à l’air infinie, d’où une pression de l’air égale à la pression
atmosphérique.

Ce modèle a le gros avantage d’avoir été testé dans divers contextes naturels qui ne sont pas
forcément ceux pour lesquels il a été conçu. On peut citer l’exemple de l’installation de
MARTHE chez un industriel du Val de Saône qui a investi dans l’achat et le calage de ce
logiciel sur des conditions hydrogéologiques précises, pour la prévision de pollution sur son
site.

-43-
 Analvse de l’offre

Il faut ajouter à MARTHE, qui simule les écoulements en milieu surtout sédimentaire, les
autres modèles présentés par le BRGM et ANTEA :
ROCKFLOW (simule des écoulements et transports en milieu à double
porosité),
CHANNET (simule les écoulements et transports dans les milieux fracturés et
imperméables).

DEDALE 3D
Le Laboratoire d’Hydraulique de France (LHF), le Centre de Recherche pour l’Hydraulique et
les Structures (CRIS) de l’Ente Nazionale per 1’Energia (ENEL) et SOGREAH développent
conjointement un système de modélisation tridimensionnel opérationnel pour les eaux
souterraines : DEDALE 3D
DEDALE 3D est un outil de simulation des écoulements, du transport et de la dispersion dans
les eaux souterraines. Parce qu’il modélise à la fois la zone saturée et la zone non saturée,
qu’il tient compte de l’influence sur l’écoulement des variations de la densité de fluide, qu’il
permet d’intégrer les hétérogénéités et anisotropies du sol et les singularités hydrauliques
(puits tranchées, drains, surfaces de suintement...), le système DEDALE est adapté à toute une
gamme de problèmes posés par la gestion ou la protection de la ressource en eau souterraine
ou des problèmes de génie civil.

DEDALE est applicable aux études de pollution des eaux souterraines comme :
- la contamination provenant de décharge,
- fuite de stockage souterrain profond ou de surface,
- pollution agricole par pesticides et fertilisants,
- pollution ponctuelle et diffuse,
- périmètre de protection,
- pénétration de solution salée...etc
ainsi qu’à des problèmes de génie civil :
- excavation drainage,
- tunnels,
- construction de barrage...
Ce système de calcul, dont le développement a été initié en 1992, a déjà fait l’objet de
nombreuses validation de terrain et continue à évoluer. DEDALE peut tourner sur des micro-
ordinateurs ou sur station.

MIKE SHE
Assez souvent cité, ce logiciel est commercialisé par le LHF en collaboration avec le Danish
Hydraulic Institute qui l’a créé. MIKE SHE est un système de modélisation des processus
majeurs de la phase terrestre du cycle hydrologique. Il simule le mouvement de l’eau et le
transfert de polluants dans les bassins hydrographiques en partant de la pluie (ou fonte des
neiges) via l’interception-évapotranspiration,le ruissellement, l’écoulement dans les cours
d’eau, l’infiltration à travers la zone non saturée, l’écoulement dans les nappes et les échanges
nappedrivières. MIKE SHE est construit autour d’un module principal MIKE SHE WM qui
simule le mouvement de l’eau auquel peuvent s’ajouter des modules spécifiques
complémentaires :
MIKE SHE AD (convection-dispersionde polluants dissous),

-44-
 Analyse de 1 ’0ffi-e

MIKE SHE GC (processus géochimiques),

MIKE SHE SE (érosion du sol),
MIKE SHE IR (irrigation).
Le module principal WM est lui même modulable suivant que l’on considère la fonte des
neige, les écoulements de surface, les écoulements en nappe ou les écoulements en zone non
saturée.
MIKE SHE a été testé et validé sur un grand nombre de bassins hydrographiques dans des
conditions hydrologiques et climatiques très diverses; il peut être installé sur des stations de
travail UNIX.

MODFLOW
MODFLOW est un logiciel de simulation pour Windows (PMWIN) pour décrire et estimer les
écoulements et le transport pour les eaux souterraines en 3D qui utilise les résolutions
d’équations par différences finies. A ce module MODFLOW, on peut ajouter différents
modules qui gèrent le transport de particules (PMPATH ou MODPATH), le transport de
soluté (MT3D) et l’estimation de paramètres (PEST). MODFLOW peut simuler l’influence de
puits, de rivières, de drains, l’impact de conditions limites, les phénomènes de recharge et
d ’évapotranspiration.
Le module MT3D qui permet de simuler le transport de solutés utilise des résolutions
hybrides Eulerienne et Lagrangienne pour résoudre l’équation de convection-dispersion en
3D. MT3D peut être utilisé pour modéliser le changement de concentration d’une espèce
miscible de polluant dans les eaux souterraines en considérant les phénomènes de convection,
dispersion et des réactions chimiques. Les réactions chimiques inclues dans ce code sont
limitées aux équilibres de sorption, aux réactions irréversibles de dégradation du premier
ordre ou la biodégradation.

SIMUSCOPP
Le logiciel SIMUSCOPP est développé dans le cadre du projet EUREKA RESCOPP par
l’Institut Français du Pétrole (développement du simulateur réalisé par Beicip-Franlab,
validation), BURGEAP (spéciations, validation) et des sociétés italiennes du Groupe ENI
(développement des interfaces).
Ce modèle numérique a pour but la simulation :
+ de la migration d’un polluant hydrocarbure ou solvant chloré de la surface d’un sol
au travers de la zone insaturée et saturée, sa dissolution puis son transport dans la nappe (aide
à l’analyse des risques),
+ d’un certain nombre de procédés de rémédiation : pompage d’huile mobile, venting
(froid ou chaud), biodégradation, injection de solutions tensioactives, injection de vapeur.
Ce logiciel peut travailler en lD, 2D ou 3D en géométrie cartésienne ou radiale, sous
maillages locaux, en milieu poreux homogène ou hétérogène (porosités et conductivités
hydrauliques). Les phénomènes physiques pris en compte sont les suivants :
- 1,2 ou trois phases présentes dans le système eau-phase organique (huile-air),
- écoulements mono, di ou tri phasiques (loi de Darcy généralisée),
- transfert de constituants entre les phases : dissolution dans l’eau, évaporation en
phase gazeuse,
- transport par convection et diffusioddispersion hydrodynamique,
- facteur retard : rétention de polluants,
- calcul des propriétés thermodynamique et d’équilibre des phases à l’aide d’équations
d’état ou de valeurs tabulées (base de données),

-45-
 Analyse de I’offi-e

- biodégradation des produits dissous et en phase résiduelle (en cours de

développement).
Les validations qui ont été effectuées se sont faites en laboratoires sur colonne, à l’échelle
d’un site pilote (SCERES) ainsi qu’in situ sur des problèmes de pompage d’huile mobile,
d’opérations de venting ou de transferts d’hydrocarbures dans le cadre d’accident.
Le code écrit en FORTRAN 77 fonctionne sur un grand nombre d’ordinateurs : PC, Stations
de travail UND(, ordinateurs centraux vectoriels ou non (CRAY, VAX, CONVEX...) .

LEACHMP
- LEACHM de New York State College of Agnculture and Life Sciences ComellUniversity.
Le LEACHM [Wagenet R.J. and Huston J.L., 19891 est un modèle d’écoulement
unidimensionnel. Il permet d’estimer l’impact due à la mobilisation et la migration d’un
polluant organique adsorbé, localisé dans la zone non saturée au dessus de la nappe.
LEACHM calcule les équilibres de distribution du polluant entre la phase liquide, gazeuse et
adsorbée sur la matrice solide.
Il utilise l’équation de Richards pour la description des écoulements. La résolution des
équations de l’hydrodynamique et du transport est effectuée suivant la méthode des
différences finies. Il prend en compte l’extraction racinaire de l’eau et des produits
phytosanitaires par les plantes, ainsi que la volatilisation. Il est possible de simuler l’évolution
simultanée de plusieurs produits et de leur métabolites. Le pas de temps de calcul est variable,
déterminé en fonction des données météorologiques disponibles. Mais ce modèle ne prend pas
en compte les écoulements dérivant des différences de densité.

Il est composé des modules :

- LEACHP pour simuler le mouvement des pesticides et leur dégradation,
- LEACHC pour décrire le mouvement des ions inorganiques (Ca, Mg, Na, K, SO4,
Cl, CO39 HC031,
- LEACHB pour décrire le développement microbien,
- LEACHW pour décrire le mouvement de l’eau seule.
Ces différents modèles s’appliquent aussi bien aux essais de laboratoire qu’aux essais in situ.
L’absorption de l’eau et du soluté par les racines des plantes est décrite par les différents
modèles. Il est possible d’obtenir le profil de températures à l’aide de LEACHP. On peut alors
modifier point par point les coefficients qui dépendent de la température.
Ce modèle a été développé pour la migration des pesticides avec des isothermes d’adsorption à
un ou deux sites.

SUTRA
SUTRA(Saturated-Unsaturated Transport) est un logiciel qui simule le mouvement de fluide
et le transport d’énergie ou de substance dissoute dans le sol. C’est un modèle bi-
dimensionnel utilisant a la fois la méthode des éléments finis et des différences finies pour
résoudre les équations qui gouvernent les processus simulés, soit :
- la densité de flux dépendante de l’écoulement d’eau dans la zone saturée et non
saturée,

-46-
 Analyse de 1 ’ofle

- le transport de solutés en eau souterraine où le soluté est soumis au phénomène

d’adsorption sur la matrice poreuse et des réactions de production et de dégradation,
- transport de chaleur dans l’eau souterraine et dans la matrice solide de la nappe.
Les écoulements souterrains sont simulés a travers une résolution numérique des équations de
bilans de masse. Le système hydrauiique souterrain peut être saturé ,partiellement insaturé ou
complètement insaturé et la densité de fluide peut être constante ou varier avec le temps
comme une fonction de la concentration du soluté ou de la température du fluide.
SUTRA peut être employé en analyse d’écoulement hydrauiique Uni ou bidimensionnel
d’aquifère ou pour la simulation du transport naturel, ou induit par l’homme, de solutés,
impliquant les processus de sorption, de production et de dégradation. L’option modélisant
les écoulements avec les différences de densités peut être utilisée dans le cas d’intrusion d’eau
salée dans un aquifère.
Ce code a été validé sur le terrain de nombreuses fois par la Société Hydroexpert entre autre et
son efficacité est réelle. Son point faible reste la convivialité : SUTRA est lourd d’utilisation
et assez complexe. Il fonctionne sur PC (Pentium nécessaire) et sur station UNIX.

Une description détaillée de ces principaux modèles se trouve dans l’annexe 5.

VII.2.Autres modèles existants sur le marché

Les modèles cités précédemment sont ceux qui ont été cité le plus souvent dans le
discours des industriels, bureaux d’études et laboratoires contactés. Vauclin en 1993 avait déjà
précédemment établi une liste non exhaustive des modèles disponibles, on peut toutefois
amener quelques’compléments (fig. 8,9)

Figure 8 :Liste non exhaustive de modèles disponibles (Vauclin, 1993)

-47-
 Analyse de 1 'ofle

Figure 9 :Liste complémentaire de modèles

-48-
 Analyse de 1 ’ofle

VII.3.Les programmes en cours, groupes de travail, organismes de

recherche et bureaux d’études
v11.3.1.Quelques growammes en cours

Projet Eureka ISMAP (Integrated System for the Management of Agricultural Pollutions).
Ce projet a pour souci de concilier 1 ’emploides produits phytosanitaires de protection
des cultures et la qualité de la ressource en eau. Cet objectif peut être atteint grâce à la
réalisation d’un système d’information et d’aide à la décision pour la gestion des pollutions
diffuses d’origine agricole. Dans ce système, la composante modélisation est appelée à jouer
un rôle important.

Dans le cadre de ce programme, une thèse, en collaboration avec le Laboratoire Hydraulique

de France (LHF-Grenoble), et le Laboratoire d’étude des Transferts en Hydrologie et
Environnement (LTHE-Grenoble), a déjà été réalisée par Vincent Guinot, sur (La
modélisation mécaniste du devenir des produits phytosanitaires dans l’environnement
souterrain. Application à la protection des captages en aquifere)). Cette thèse a fait l’objet
d’une validation du code MIKE SHE sur un site-test du projet ISMAP (Bassin de Grand
Morin, région parisienne) au cours de laquelle une réflexion a été menée sur la capacité de
l’outil numérique à modéliser les processus qui se produisent dans le sol. Le but de la thèse a
été de mettre en place une méthodologie qui permettra de déterminer les quantités maximales
de produits phytosanitaires que l’on peut tolérer à l’épandage, sans risquer d’entraîner, au
niveau des captages d’eau potable, un dépassement des normes de potabilité.

Projet Eureka RESCOPP

Le projet Rescopp est un projet européen sur la restauration des sols pollués, qui a
débuté en 1993 et qui s’achève fin 1996. Les objectifs de ce programme sont de répondre à la
demande en analyse de risques et en prédiction.
Dans le cadre de ce projet, l’IFP, BURGEAP et des sociétés italiennes du groupe ENI ont
développé le logiciel SIMUSCOPP. Le logiciel s’attache à la simulation :
-+ de la migration d’un polluant hydrocarbure ou solvant chloré de la surface d’un sol
au travers de la zone insaturée et saturée, sa dissolution puis son transport dans la nappe (aide
à l’analyse des risques),
-+ d’un certain nombre de procédés de rémédiation : pompage d’huile mobile, venting
(froid ou chaud), biodégradation, injection de solutions tensioactives, injection de vapeur.

Le programme PNETOX
L’écotoxicologie est un élément incontournable et indispensable dans l’évaluation des
risques liés à la présence de contaminants dans l’environnement. Les besoins en
écotoxicologie, de la part des industriels, des pouvoirs publics et des directions
opérationnelles du ministère ont donc conduit à la mise en place en 1995 d’un plan d’action et
d’un programme sur l’écotoxicologie, mené par le club CRIN-environnement avec
l’assistance de 1’ADEME.

Le programme PNETOX s’appuiera sur les programmes en cours dans les autres organismes
( S U E , ADEME, INERIS, IFREMER, INRA,CNRS, CEA, CEMAGREF). Il jouera un rôle
de coordination et d’animation par rapport à ces programmes et de structuration de

-49-
 Analyse de I’ofie

l’écotoxicologie. Il développera et soutiendra les recherches complémentaires à celles déjà

entreprises en même temps que de nouvelles recherches pluridisciplinaires prenant en compte
l’évaluation du risque et celles pouvant aboutir à une valorisation.
Les recherches portent sur :
-+ la précision de la capacité d’acceptation des milieux à recevoir des charges toxiques
(contamination par l’eau, l’air ou le sol), tout en maintenant une qualité satisfaisante de
l’environnement,
+ l’évaluation correcte du risque lié à la mise sur le marché de substances ou à des
pollutions installées ou susceptibles de se produire. Les tests existants sont à améliorer et de
nouveaux sont à mettre en place. Des modèles d’exposition sont à développer.
Les thèmes de recherche retenus concernent l’effet biologique des polluants ainsi que la
notion de biodisponibilité d’un contaminant.

Contrat de Plan Etat-Région

Dans le cadre du Xerne Contrat de Plan Etat-Région, le Programme Fédérateur Environnement
a fmancé des crédits au Conseil Régional Rhône-Alpes ainsi que des organismes d’Etat. Un
tableau récapitulatif des moyens attribués à l’ensemble des projets retenus est présenté en
annexe 6 .

Caracas
Depuis 1994, il a été lancé au sein de l’Europe plusieurs programmes concernant
l’environnement, comme le EC R&D Programme Environnement and Climate qui s’achèvera
en 1998. Ce programme a pour sujet de concertation le recyclage, le traitement des déchets et
le traitement de site contaminé. A cette occasion, différents projets ont démarré comme le
projet CARACAS (Concerted Action on Risk Assessement for Contarninated land qui a
débuté le 1“ février 1996 et d’une durée de trente mois.
Les objectifs principaux sont :
-+ d’identifier, de compiler, d’évaluer et de réajuster les différents projets et approche
scientifique pour l’évaluation des risques, actuellement développés dans les états membres de
l’Union européenne,
-+ de proposer des priorités scientifiques pour les futurs projets et programmes de la
commission européenne et des états membres de l’Europe,
+ d’élaborer des directives et des recommandations pour l’évaluation du risque sur les
sites contaminés.
Cette action concertée est menée par 14 participants des états membres de l’union européenne:
Belgique, Autriche, Danemark, Finlande, France, Allemagne, Grèce, Irlande, Italie, Pays-Bas,
Portugal, Espagne, Suède et Grande-Bretagne. Des pays tel que la Norvège, la Suisse, les
USA et l’Australie apportent leur contribution par leur expérience et leur connaissance sur ce
projet.
Les thèmes de travail sont les suivants :
+ toxicologie humaine,
+ évaluation du risque écologique,
+ devenir et transport des contaminants,
+ analyse et étude de site,
-+ modélisation,
+ méthode pour l’évaluation des risques,
+ bilan sur les normes : inventaire et directives.

-50-
 Analyse de l’ofie

Le projet CARACAS progresse en collaboration avec d’autres initiatives internationales dans

le domaine des sols contaminés, comme le projet NICOLE (Network industrially
Contaminated Land in Europe), ISOKEN-Working Group, EEA-Topic Center Soil, Ad Hoc-
working Group sur les sols contaminés.
Plus particulièrement, sur le thème du devenir et du transport des contaminants, le groupe doit
passer en revue la connaissance actuelle sur les processus biologiques, chimiques, physiques
qui interviennent dans le devenir et le transport des contaminants organiques et inorganiques
pour les sols et les eaux. L’accent est m i s sur le fait que ce sujet devrait être, à l’avenir, intégré
à la pratique d’évaluation de risques. Ce thème est abordé par l’analyse et la discussion autour
des différents processus tel que la sorptioddésorption, la biodégradation, la dispersion, la
diffusion.. .

VII.3.2.Groupede réflexion

Le club CRIN, déjà cité à l’occasion du programme PNETOX, regroupe des

laboratoires et centres de recherche, l’administration et des entreprises. La collaboration de
tous ces membres a débouché sur un travail en réseaux, où la fonction du club a pris un
caractère volontairement complémentaire des structures fédératives qui existent dans chacun
de ces milieux. Le club CRIN désire accentuer ce fonctionnement en réseaux organisés en
développant une coopération avec certaines associations (Entreprises pour l’Environnement,
Chimie et Ecologie, RE.CO.R.D....) et constituer des relais régionaux vers les collectivités
territoriales, acteurs d’importance croissante. Un de leur thème de travail s’articule autour de
la problématique des sols en considérant différentes approches :
+ la bioaccessibilité,
-+ le devenir des xénobiotiques,
+ lesmétaux.
Sur l’approche <unétam),le club CRIN s’est donné un programme de travail comportant
quatre volets, à savoir, les teneurs du sol en éléments traces, les friches et sites pollués, la
pollution diffuse des sol et les transferts de polluants entre le sol et les autres milieux.

VII.3.3.Laboratoiresde recherches et écoles

Le Laboratoire d’études des Transferts en Hydrologie et Environnement

Contact : Mrs GAUDET et VAUCLIN
Suite à l’évolution de l’Institut de Mécanique de Grenoble vers 4 laboratoires différents, l’ex-
groupe Hydrologie est devenu en Janvier 1992 le LTHE, unité biappartenante à l’Institut
National Polytechnique de Grenoble et l’université Joseph Fourier (Grenoble).
L’orientation générale des activités de recherche est focalisée sur des études de transferts
polyphasiques couplés de masse et de chaleur dans les matériaux poreux et des étude de la
phase continentale du cycle des eaux en terme de quantité et de qualité.
Les recherches sont réparties entre 4 équipes travaillant en étroite collaboration :
+ l’équipe Transferts en Matériaux Poreux (TMP) (Resp. J.P. Gaudet),
+ l’équipe Hydrologie de la Zone Non Saturés (HLNS) (Resp. G. Vachaud),
+ l’équipe Hydrométéorologie (HMET) (Resp. J.D. Creutin),
+ l’équipe hydrologie de surface (HSURF) (Resp. Ch. Obled).

-5 1-
 Analvse de l’ofie

Les perspectives de recherches vont s’orienter, entre autre, vers la modélisation des transferts
couplés d’eau, d’énergie et de solutés avec la prise en compte des équilibres géochimiques
dans le milieux poreux (sols, roches, matériaux de construction). Cette optique inclue plus
précisément :
--+ le traitement des non linéarités dans les processus de transfertltransport,
--+ le développement et l’introduction de procédures d’estimation des paramètres par
méthode inverse,
--+ les techniques de changement d’échelles et de macro (méga) scopisation des processus
en milieu hétérogène.
Pour les équipes TMP et HLNS, les buts à atteindre en ce qui concerne la modélisation des
transferts de polluants dans les sols seraient d’identifier les couples Polluants/Sols les plus
couramment rencontrés afin d’établir une méthodologie de travail efficace.
La caractérisation de la nature de la variation spatiale des paramètres du sol, des fluides et des
solutés constitue un thème de travail important à développer.

Rq : un groupe de travail s’est constitué en collaboration avec P. Verbelin (ENSEPI) dans le

cadre d’un programme Environnement TTEC sur l’écotoxicologie.

Laboratoire d’Analyse Environnementale des Procédés et Systèmes Industriels

INSA-lyon
Contact : Mrs NAVARRO et MOSZKOWICZ
Le LAEPSI a pour thèmes de recherche 6 axes principaux :
--+ analyse des cycles de vie (ACV),
+ organisation du système déchets et analyse des procédés,
--+ inertage des déchets ultimes,
+ décontamination des sites pollués,
+ traitements des effluents,
--+ génie des procédés appliqués à l’environnement.
Quelques thèse en cours :
* Etude des transferts de solutés dans les matériaux poreux soumis à lixiviation (Florence
Sanchez).
* Etude de la réversibilité de l’adsorption des polluants organiques aromatiques dans les sols
pollués (Rémy Bayard).
* Etude du transfert de polluants volatils en phase gazeuse pour le diagnostic rapide de sols
pollués par des essence ou gasoil (Pascale Le Goff).

Laboratoire des Sciences de l’Environnement

Contact : M. T. Winiarski, Melles Y. Bouvet et C. Delolme

L’équipe interdisciplinaire de recherche du Laboratoire des Sciences de

l’Environnement (LSE) travaille sur l’analyseet la modélisation de l’impact des activités
humaines sur l’environnementnaturel.
L’impact des centres d’enfouissement techniques sur la qualité des eaux souterraines et
Z‘évaluation des risques pour l’environnement (pollution des eaux et des sols, santé humaine)
liés à la qualité des lixivats produits par les décharges constituent le premier thème. Les effets
des aménagements (barrages) sur la capacité de prise en charge des effluents industriels et
urbains par les fleuves représentent le second champ de recherche. Le troisième thème

-52-
développe une méthodologie de diagnostic des sols contaminés friches industrielles) et de la
qualité des eaux.
Les outils développés pour traiter ces thèmes sont: un laboratoire de chimie analytique, des
modèles physiques et biologiques, des stations d’acquisition de données.
La problématique des sites contaminés s’oriente plus précisement vers :
-+ l’étude de la microflore du sol (C Delolme),
+ l’étude de la matrice solide du sol (T. Winiarski),
+ l’amélioration, la validation et l’utilisation de modèles (C Delolme et T. Winiarski)
Cette double approche permet de tester les modèles numériques sur modèle physique.
Exemples de modèles physiques utilisés :
. fond de décharge-zone non saturée-aquifêre à l’échelle 1:100,
. bassin d’infiltration d’eaux pluviales-zone non saturée à l’échelle 1.
Exemples de modèles numériques utilisés :
POLLUT2D dans le cas d’un fond de décharge,
LEACH dans le cas d’un fond de décharge et d’un bassin d’infiltration.

Quelques thèse en cours :

* Contribution des microorganismes dans le transfert des métaux lourds dans le sol (J.
Crosnier).
* Modélisation du devenir des métaux lourds dans la zone non -saturée du sol. Prise en
compte du compartiment microbien (C. Hébrard).
* Relation entre les fluctuations de l’activité microbienne et l’état de perturbation d’un sol (B.
Brohon).

Laboratoire des sciences du nénie chimiaue

Contact : M. SARDIN
M. SARDIN utilise le code IMPACT (modèle de réacteurs en cascade).
M. SARDIN a beaucoup travaillé sur les pollutions minérales (métaux lourds) et s’intéresse
depuis 2-3 ans aux HAP et COV. Son laboatoire a surtout travaillé jusqu’à maintenant en zone
saturée. Pour la zone insaturée, M. SARDIN soulève le problème de l’aspect triphasique et le
problème de la modélisation hydrodynamique dans cette zone.

Laboratoire de Géologie Amliquée-Université Paris VI

Contact : M. De MARSILY
Le Laboratoire de Géologie Appliquée organise ses recherches autour des thèmes suivants :
-+ fonctionnement des systèmes fluviaux anthropisés,
+ écoulements souterrains et interactions fluides-solides,
+ métallogénie,
-+ géophysique non linéaire et analyse des précipitations.

L’Institut de Mécanique des Fluides

Contact : M. ACKERER
L’IMF constitue 1’URA CNRS 854, à l’université Louis Pasteur de Strasbourg. Les pôles de
travail s’articulent autour de 4 axes en modélisation à I’IMF :
+ transport et écoulement en milieux poreux saturés,
+ méthode inverse,

-53-
+ modèle de transport des substances réactives (phénomènes de spéciation avec les
métaux lourds),
+ aspect multiphasique
Ils s’intéressent par exemple a la dissolution des différents constituants d’un gasoil dans l’eau.
Les méthodes classiques des éléments finis n’étant pas suffisantes (surtout en milieu
hétérogène), I’IMF développe une méthode d’éléments finis mixtes hybrides et discontinus.
Intérêt : temps de calcul plus court.
Contrainte : méthode explicite (problème des pas de temps).
Développé en 2D et 3D, en régime permanent ou transitoire.
Rq : L’IMF possède Swanflow (partie écoulement d’hydrocarbure).

Afin de lutter contre toute forme de propagation souterraine d’un contaminant, il est important
de connaître les mécanismes relatifs au cheminement d’un polluant. C’est dans ce but que le
Laboratoire d’Hydrodynamique des Milieux Poreux (LHMP) a mis en place le modèle
MARCEAUS (Modèle d’Analyse et de Recherche sur la Contamination des Eaux
Souterraines). Ce dispositif permet de reconstituer, dans une cuve de laboratoire, une portion
de nappe phréatique artificielle. Pour représenter au mieux la réalité physique, la cuve
expérimentale est constituée d’un milieu hétérogène. La cuve expérimentale est remplie d’un
milieux poreux qui occupe un volume de 500x100x95 cm3.

L’Institut Franco-Allemand de Recherche sur l’Environnement

Contact : M. RAZAKARISOA
L’objectif de 1’IFARE est la protection des aquifères contre la pollution par hydrocarbures et
dérivés.
Cette protection passe par :
+ la compréhension des mécanismes,
+ les méthodes de décontamination,
+ les modèles mathématiques.
L’IFARE travaille sur un grand dispositif expérimental baptisé SCERES (Site Contrôlé
Expérimental de Recherche pour la Réhabilitation des Eaux et des Sols) dans le but
d’améliorer la maîtrise des phénomènes de contamination d’aquifère alluvial par infiltration
d’hydrocarbures. L’outil SCERES, constitué par un bassin étanche de 25x12x3m, implanté
sur le campus de Strasbourg-Cronenbourg, permet de faire progresser les connaissances sur la
compréhension des mécanismes de migration du polluant et la détermination quantitative du
domaine contaminé. Le projet SCERES associe des équipes universitaires (CNRS, universités
fiançaises et allemandes), et des partenaires du secteur industriel et de la technologie avec
Anjou Recherche de la Générale des Eaux, le Burgéap et 1’IFP.

Laboratoire de Géomécaniaue - ENS de Géologie de Nancy

Contact : M. BUES
Le Laboratoire de Géomécanique travaille sur la migration de polluants en zone saturée et
s’attache plus précisément à décrire:
+ le transfert de micropolluants réactionnels, aspect bio-physico-chimique (code calcul
bio-physico-chimique),
+ la biodégradation des pesticides dans les sols (phényl acétate de mercure : accident de
Sandoz en 1986) : transport et interaction avec le sol,
+ le transport de produits de différentes densités (viscosité élevée) (solutions salines par
exemples),

-54-
 Analyse de 1’ofle

+ la modélisation en 1D pour interpréter les essais sur colonnes (convection, dispersion,

cinétique de dégradation),
+ la modélisation en 3D : viscosité et densité différentes.

L’Ecole des mines de Paris

Contact : M. LEDOUX
Le Centre d’Informatique Géologique consacre son activité à l’étude du cycle de l’eau et plus
particulièrement des transferts de fluides et de matières dans le sous-sol et à l’interface sol et
sous-sol. L’activité de recherche du centre est orientée autour de trois thèmes principaux
couvrant une large gamme de problèmes environnementaux liés à la partie continentale du
cycle de l’eau :
+ modélisation des systèmes hydrogéologiques et hydrologiques,
+ modélisation des transferts de matière et de chaleur et des interactions fluide-roche
dans le milieu souterrain,
+ étude des précipitations et des écosystèmes lacustres et fluviaux.
Modélisation des systèmes hydrologiques et hydrogéologiques
Ce thème de recherche traditionnel du laboratoire a donné lieu à des développements
théoriques portant essentiellement sur :
+ la modélisation hydro-thermo- mécanique des bassins sédimentaires,
-+ la modélisation tridimensionnelle de la structure hydrogéologique du Bassin de Paris,
entreprise pour le compte du CEA,
+ le changement d’échelle sur les perméabilités et les paramètres de transport dans les
réservoirs très hétérogènes (application directe pour les réservoirs pétroliers et sites de
stockage profond de déchets radioactifs),
+ la modélisation du transport en milieu fissuré (équipement hydroélectrique souterrain
pour EDF, stockage de déchets nucléaires).

Modélisation du transfert de matière et de chaleur et des interactionsfluides-roches.

Cette activité s’articule autour des problèmes de géochimie continentale et plus
particulièrement des transferts de matières minérales sous l’effet de circulation de fluides. ces
travaux trouvent leur application dans le domaine d’altération et accumulations superficielles,
dans la gestion de l’environnement des sites miniers et industriels, dans la cimentation des
réservoirs pétroliers et dans les expertises géochimiques sur des sites de friches industrielles.
Ce thème vise à l’identification, à la mise en oeuvre et à la validation de nouveaux concepts
dans les modèles. Elle s’appuie sur une expérimentation en laboratoire ou in-situ menée en
collaboration avec des organismes de recherches extérieurs.

Les perspectives
Les problèmes environnementaux se rapportant au sous-sol, tels que les pollutions
industrielles ou l’élimination des déchets, constituent actuellement un thème de recherche
porteur. De plus, la problématique scientifique du Centre d’Informatique Géologique, alliant
les moyens de caractérisation des milieux au potentiel de quantification grâce à la
modélisation, vont lui permettre de développer et promouvoir des méthodes et de réaliser des
expertises.

-55-
 Analyse de 1’offi.e

L’école des mines de Saint Etienne

Contact : Mrs BOURGOIS et VAILLANT
Le département “Ingénierie de l’Environnement’’ a été créé à l’Ecole des Mines de Saint-
Etienne en 1991. Les principaux domaines de recherche abordés, liés à l’environnement
industriel, sont:

la maîtise spatiale des risques (pollutions, inondations),

la modélisation des phénomènes de propagation de pollution (ressources en eau)
le diagnostic et l’aide à la décision pour la réhabilitation des sols pollués,
l’analyse des risques techniques et organisationnels en conception et maintenance des
réseaux,
la phention du BLEVE en cas d’incendie,
la conception d’outils d’aide à la décision pour le traitement de déchets industriels,
l’analyse des dysfonctionnements d’installations automatisées et l’aide à la conduite à
partir de systèmes
experts temps réel.

Les thèses en cours dans le laboratoire Pollution-Eau-Sol sont :

Evaluation du risque potentiel sur l’environnement lors de la réhabiliation d’un site
contaminé (V.L. Duong).
Modélisation des propagations de pollutions biphasiques en milieihétérogène (zone
saturée ou non) -application au cas de l’accident ferroviaire de Chavanay (Loire) (W.
Anker).

L’école des mines de St Etienne contibue au Centre National d’Etudes et de Recherche sur
les Sites et Sols Pollués de Douai via un groupement entre EMSE, EMP et ARMINES.

L’Ecole des mines de Douai

Contact : M. GUILLERMO et Melle SOUTO
L’école des mines de Douai a mis en place un Centre National d’Etudes et de Recherches
sur les Sites et Sois Pollués.

Actuellement le département Chimie et Environnement ne travaille pas sur la modélisation de

migration de polluants dans les sols. Toutefois, il y a une thèse en cours (depuis environ 1 an
et demi) sur “Etude d’un site pollué par les métaux - Evaluation du degré de contamination
des sols et prévision de son évolution” (c’est surtout l’aspect interaction sol-polluants qui est
traité).

Centre National d’Etudes et de Recherches sur les Sites et Sols Pollués

La première pierre du Centre a été posée début avril par Madame le Ministre de
l’Environnement.
Les statuts de ce centre prévoient deux types de partenaires :
noyau dur : partenaires mettant à disposition du centre 2 cadres à temps plein + 1
bourse de thèse + 125 KF/an. Le noyau dur décide des thèmes de recherche à effectuer dans le
centre.
membres associés : partenaire mettant à disposition du centre 250 KF/an. Les membres
associés ne sembleraient avoir qu’une voix consultative sur les thèmes de recherche à effectuer
dans le centre.

-56-
 Analyse de 1 ’ofle

Partenaires actuels : ARMINES-EMP-EMSE, BRGM, CDF, CEA, CGE, EM Douai,

GDF, INEFUS, Institut Pasteur, ISA, Rhône-Poulenc Chimie, Usinor, Université du Nord.

Les axes généraux des recherches qui seront réalisés dans ce centre sont actuellement définis
comme ci-dessous :
caractérisation des sites pollués :
. spéciation,
. devenir des polluants,
. mobilité/transport,
. bilan pollutiodindustrie,
. échantillonnage,
analyse de risques,
faisabilité de procédés.

VI7.3.4.Organismes de recherche. grands grouves et sociétés

Ce sont en général des grosses structures qui ont une connaissance poussée en la
matière. Du fait de leur diversification, on y trouve souvent un pôle de recherche, en relation
avec leur production, qui peut mener des études fort intéressantes et qui peut être
particulièrement avancé dans un domaine, suivant leur ancienneté. C’est le cas de I’IFP et de
Rhône-Poulenc.

L’IFP
Contact :M. Le THIEZ
Ce grand groupe possède une activité de recherche très développée et collabore avec de
nombreux laboratoires (cf. collaboration avec I’IMF, et I’IFARE de Strasbourg). L’IFP
développe actuellement un logiciel de migration de polluants organiques dans les sols et sous-
sols (2s et ZNS, di ou triphasique, migration et transfert du milieu poreux vers la nappe) dans
le cadre du programme européen RESCOPP : SIMUSCOPP.
Un des pôles d’activité de I’IFP est orienté vers la modélisation de procédés de remédiation :
pompage, venting, biodégradation (ZNS et ZS), injection de solution de tensioactif.

Pour l’avenir, I’IFP compte axer ces efforts sur :

+ l’introduction de stratégie d’échantillonnage,
-+ la géostatistique pour la prise en compte de l’hétérogénéité des milieux poreux,
-+ les procédures de calage automatique,
-+ l’amélioration de la connaissance de la biodégradation et des transfert physique.

ANTEA-BRGM
Contact : M. ROYER ( ANTEA) et M. THIERY (BRGM)
De l’étroite coopération de ces deux structures, il en découle une volonté d’être à l’écoute des
industriels en particulier en tant que maître d’oeuvre par l’intermédiaire d’ANTEA tout en
progressant en recherche appliquée et fondamentale via le BRGM et ses collaborations avec
de multiples laboratoires de géologie, géophysique, géochimie...en France et à l’étranger.
ANTEA travaille avec les industriels sur le stockage de déchets et également sur le diagnostic
des sols pollués. Ils sont maître d’œuvre pour les traitements de type land farming, traitement
de lixiviat (incinération ou mise en décharge). Le BRGM et ANTEA développe et
commercialise deux logiciels :

-57-
+ MARTHE : Modélisation d’AquiEre par un maillage Rectangulaire en régime
Transitoire pour le calcul Hydrodynamique des Ecoulements : simulation hydrodynamique en
nappe en milieu sédimentaire, milieu poreux continu, différences finies.
+ ROCKFLOW : simulation écoulement en zone de socle ou milieu fracturé, utilisé pour
site de stockage sous granite, éléments finis, nécessite beaucoup de données, milieu à double
porosité, approche stochastique de la fracturation.
Un autre logiciel est également développé : le logiciel CHANNET pour les milieux fracturés
(matrice imperméable).

Il faut souligner la parution d’un Guide Méthodologique, demandé au BRGM par le

Ministère de l’Environnement : ((Gestion des sites (potentiellement pollués)))
Contenu du mide :
Cet outil a été mis en place dans le souci d’une approche nationale cohérente et homogène du
problème des sites et sols pollués. L’élaboration de cet ouvrage est menée par le Ministère de
l’Environnement avec des organismes experts, le BRGM et I’INERIS et dans le cadre de
groupes de travail nationaux. Il se compose de plusieurs petits guides méthodologiques ainsi
qu’un grand nombre d’annexes dans lesquelles sont répertoriées toutes les informations
législatives et techniques relatives au problème.
Un premier guide pose les directives à suivre pour un prédiagnostic sur les sites en activités et
les sites abandonnés. Il est suivi par un guide sur 1’Etude des sols, un autre sur 1’Evaluation
simplifiée des risques et un dernier sur 1’Evaluation simplifiée des risques et classement des
sites. Dans les annexes, sont regroupés les lois, les normes pour l’eau, les sols, des
recommandations pour la collecte des échantillons préalable au diagnostic, des techniques, les
sources d’informations, les notes de danger, un glossaire des sites et sols pollués....
Présenté de façon claire et concise, cet ouvrage constitue un outil d’aide à la décision très
utile. Toutefois, il convient de considérer ces guides comme des outils <provisoires» qui
seront régulièrement modifiés et actualisés pour prendre en compte d’une part l’évolution des
connaissances et d’autre part le retour d’expérience de leur utilisation.

GEOSTOCK
Contact : Mrs. COLIN et GENESTE
Géostock est maître d’œuvre en dépollution. Il essaye de promouvoir les techniques in situ
(lessivage, pompage, écrémage, technique bio en excavé).
Géostock utilise depuis longtemps des modèles appliqués au stockage soutterain (modèle
hydraulique sans transport, modèle 2D hydraulique (SILMULOG) (éléments fuis)). Ce sont
des outils d’origine pétrolière.
Ils possédent également un modèle récent pour les stockages en résevoir : Workbench 3D
(avec module transport et module multiphasique).
En migration de polluants dans les sols, ils possèdent :
. Porflow : fait à l’origine pour les radioéléments, ZNS, volumes finis, écoulement et
transport, complet en physique, multiphasique, tourne sur une station X (apprentissage
délicat). Logiciel distribué par ACRI (Sofia Antipolis),
. AQUA 2D : éléments finis, écoulement et transport, uniquement nappe, tous types de
polluants, données : coefficient de rétention et constante de dégradation, paramètres
dimensionnels, paramètres hydrauliques, conditions limites spatiales et temporelles, porosité,
transmittivité...tourne sur PC, outil convivial,
.un petit logiciel américain : SPILLUOL.

-58-
 Analyse de I’ofie

HYDROEXPERT
Contact : Mrs BONNET et OLTEAN
HYDROEXPERT, possède plusieurs antennes dans la région parisienne, ainsi que dans le sud
de la France. Cette société d’ingénieurs conseil, indépendante de toute entreprise ou groupe
financier, possède entre autre SUTRA (plus MODFLOW, TALISMAN..) dont ils ont
beaucoup développés l’interface (pré-processeurs et post-processeurs). Leurs services sont
axés sur l’évaluation de projets d’hydraulique rurale, les problèmes de sédimentologie marine
et fluviale, d’hydrologie quantitative, d’impact et de sûreté.
Un pôle recherche met en oeuvre pour ses besoins propres ou à la demande de ses clients, des
procédés, des équipements ou des méthodes de mesures innovants comme la fiacturation
hydraulique, des tensioperméamètres ou des carottiers légers.

SOGREAH Ingénierie
Contact : Mrs Le BIDEAU et JARDIN
Située à Grenoble, filiale du groupe Sogelerg-Sogreah, Sogreah Ingenierie est une société
d’Ingénieurs Conseil axée sur les techniques de l’hydraulique et de l’environnement
principalement en vue de l’aménagement du territoire.
Leurs activités sont Concentrées sur l’environnement urbain et industriel, l’Ingénierie fluviale
et hydroagricole, les barrages et aménagements et les ports et les côtes.
Cette société collabore étroitement avec la filiale de recherche, le LHF et participe de ce fait a
de nombreux projet. Les modèles, m i s au point ou améliorés par le LHF sont commercialisés
par ces deux partenaires et peuvent être installés et calés sur une situation, chez un industriel,
dans le but d’une prévision des risques de pollution.

Le Laboratoire d’Hvdraulisue de France

Contact : Mrs CUNGE et BOISSON
Le LHF a été créé en 1988, il constitue une des premières sociétés privées fiançaises que
l’industrie gère en collaboration avec la recherche. Leurs compétences sont axées sur le
domaine de l’eau, de l’environnement et de l’informatique. La collaboration de leur division
numérique et physique avec la participation de 1’INPG et de I’UJF de Grenoble a permis le
développement de plusieurs codes de calculs dont les plus importants, en ce qui concerne le
sujet de l’étude, sont DEDALE 3D et MIKE SHE.
Le développement d’autres codes de calculs ne constitue plus un pôle d’intérêt, ou une
difficulté pour eux. Leurs efforts se veulent concentrés sur l’amélioration et la promotion de
leurs produits actuels par le biais de validation de perfectionnement sur site pilote (à trouver)
et de modification de programme incluant l’approche géochimique, d’une part, et de
l’influence de la microbiologie d’autre part.

TRANSOFT
Contact : Mme RENARD et M. TRIPATHI
TRANSOFT offre une gamme de logiciels provenant du marché de la recherche américaine,
destinés à étudier les phénomènes de mécaniques des fluides et de structure en utilisant les
technologies actuelles. Leurs logiciels sont en général disponibles sur la plupart des
ordinateurs. TRANSOFT développe beaucoup les interfaces des logiciels récupérés afin de les
rendre utilisable par le grand public. Il sont demandeurs de sites-test afin de valider leur
modèle.

-59-
 Analyse de I ’offre

Faculté Polytechnique de Mons - Laboratoire de Géologie

Contact : M. SQUERENS
Le Laboratoire de Géologie utilise des modèles de migration de polluants en zone saturée. Ils
utilisent le logiciel AQUA3D (modèle islandais) commercialisé par la société Natural Systems
Software. Ce laboratoire est intéressé par des modèles prenant en compte la zone non saturée.
Il illustre cette demande par un exemple de pollution : cas d’une citerne qui s’est renversée sur
le remblais d’une autoroute, la pollution n’est toujours pas apparue au niveau de la nappe !
(alors que seulement 3 O00 tonnes de carburant, sur environ 17 O00 perdues, ont été pompées).
La pollution est certainement retenue par les craies.
Mr. Squerens estime qu’un modèle chimique est pratiquement impossible à réaliser du fait du
grand nombre de réactions pouvant intervenir.

University of Illinois WSA)

Un certain nombre d’industriels (Amoco, Chevron, Hewlett-Packard, Mobil, Texaco.. .) ont
constitué en 1986 un consortium pour coopérer sur le programme “Hydrogelogy” (Annexe 7).
L’objectif est le développement d’outils de calculs pour la simulation des processus
hydrogéologiques et géochimiques. Les recherches sont effectuées par un groupe de
chercheurs de l’université de l’Illinois. Des logiciels sont élaborés et commercialisés. Le
logiciel Géochemist’s Workbench permet de modéliser les interactions géochimiques de
nombreux composés avec les roches et les sols (voir chapitre V.5).

OUESTOR
Contact : M. STEWART
Queen’s University Sciences and Technology Research Centre (Dublin, Irlande).
QUESTOR est une structure de mutualisation de crédits de recherche pour un groupe
d’industriels dans le domaine de l’environnement. A la différence de RE.CO.R.D.,
QUESTOR possède ses propres laboratoires, situés dans une université.
La thématique de la modélisation des transferts de polluants dans le sol a été abordée par les
chercheurs de QUESTOR et se développe. Deux projets sont actuellement en cours dans ce
domaine (annexe 8) :
. modélisation numérique de la dispersion des polluants dans la zone non saturée : il
s’agit de résoudre les problèmes numériques de résolution des modèles rencontrés surtout
dans les cas de la prise en compte des réactions chimiques. Une validation des résultats sur
sites est également prévue,
. modélisation biogéochimique du devenir de métaux dans les milieux saturés et non
saturés : les interactions de ces métaux avec les composés organiques dans les sols doivent
être spécifiquement étudiées et modélisées pour être prises en compte.

Il est évident que les membres industriels de QUESTOR considèrent comme prioritaire les
recherches en matière de caractérisation de la mobilité des polluants dans les sols. Les études
sont menées par des équipes très pluridisciplinaires associant chimistes, biochimistes et
numériciens.

QUESTOR est très ouvert a une collaboration avec RE.CO.R.D.

-60-
 Analyse de l’ofie

Il est important de souligner que faire un inventaire exhaustif de tous les laboratoires,
écoles, entreprises et organismes qui travaillent de près ou de loin sur la pollution des sols est
impossible. La revue ci-dessus se veut uniquement indicative et non pas directive quant aux
prises de contacts qui pourront s’établir par le futur. Il est bien entendu, qu’il existe en fait, un
grand nombre de travaux en cours sur le sujet. Entre autres, ne sont pas cités tous les acteurs
de la recherche en microbiologie, en géochimie, thèmes qui doivent être pris en compte dans
la modélisation d’aujourd’hui.
Un certain nombre de renseignements, concernant les thèmes de travail et les
perspectives des labo, organismes et entreprises sont disponibles sur les serveurs d’htemet.

-61-
 Conclusions et DersDectives

VIII. CONCLUSIONS ET PERSPECTIVES

A l’issue de cette étude, il est possible de dégager les principales questions et les
problèmes qui restent posés concernant la modélisation du transport des polluants dans les
sols.
La problèmatique nécessite une approche très pluridisciplinaire, comme toutes les
sciences de l’Environnement en général. Si de nombreux spécialistes de l’hydrogéologie, de la
géochimie et de la bio-physico-chimie maîtrisent bien des aspects particuliers du domaine
d’étude, par contre leur association dans une approche commune et globale ne se fait que très
difficilement. La collaboration avec des numériciens ou des informaticiens est certainement
encore plus problématique.
Comme pour toute modèlisation, les données d’entrée restent LE problème. En
effet, le sol et en particulier sa zone non saturée sont très hétérogènes et se pose alors le
problème des mesures et de leur acquisition. Les données sont-elles représentatives du milieu
à étudier ? Sinon, quelle est la signification des résultats obtenus ?
Comme pour toute simulation, il apparait aussi la nécessité de la fiabilité et de la
validité des prévisions obtenues par le calcul.
Si, a priori, il ne semble pas qu’il y ait des avancées importantes à attendre de
l’amélioration des schémas numériques, il existe sûrement de nombreux cas où certains
phénomènes physiques, physico-chimiques ou biologiques restent mal modélisés.

En se rappelant que modéliser c’est d’abord hiérarchiser les phénomènes régissant une
réalité complexe (examen des connaissances existantes, écritures des phénomènes sous forme
d’équations mathématiques, résolution numérique et validation), nous pouvons faire les
propositions d’études suivantes :

* l’évaluation de la pollution des sols étant une approche relativement nouvelle, il nous
parait important de rester pragmatique dans un premier temps et de recenser les connaissances
déjà acquises. Ce recensement pourrait faire l’objet de la création d’une base de données
concemant :
- les interactions polluants/sols afin de répondre aux questions : quel est le
devenir d’un polluant dans un sol donné ? Au bout de combien de temps le polluant examiné
sera-t-il dégradé naturellement ? Le sol peut-il être réutilisable et quand ?
- les polluants types générés par les différentes catégories d’industrie, dans le
but de cerner au mieux les recherches de polluants lors d’une évaluation de site.
Cette base de données ne saurait voir le jour sans une participation active des
diférents acteurs du domaine,

* mise en oeuvre de calculs couplant les différentes approches (géologiques,

chimiques, physiques, biologiques, ...), ce qui entraîne la participation d’équipes
pluridisciplinaires ou un rapprochement de différentes équipes de recherche. Cette intégration
favoriserait une vision globale des phénomènes afin de se rapprocher le plus possible de la
réalité,

-62-
 Conclusions et perspectives

* avant d'utiliser tout modèle, il est nécessaire d'améliorer la méthodologie

d'acquisition des données afin de représenter au mieux le site à étudier, cette amélioration
passant par un recensement des connaissances sur les sols et les sous sols après discussion
avec des géologues, hydrogéologues, utilisateurs des sites, ... Ce recensement pourrait aboutir
à la création de cartes en 3D. La gestion des données par un Système d'Information
Géographique servirait à paramétrer les modèles et à analyser dans l'espace les résultats de
simulation,

* les paramètres d'entrée avec quelle précision ? La réponse à cette question pourrait
être faite en étudiant la sensibilité des modèles aux différents paramètres d'entrée et la
hiérarchisation des paramètres clés.
En effet, d'une part, il n'est pas toujours nécessaire de connaître avec une grande
précision certains paramètres dont la sensibilité est faible vis-à-vis des simulations obtenues
par le modèle. Et d'autre part, la précision nécessaire pour certains paramètres ne peut être
obtenue par la mesure sur le terrrain, et, par conséquent, la simulation ne pourra donner que
des résultats approchés voire erronés,

* en ce qui concerne la validation des modèles, deux voies peuvent être pressenties :
- sur un site Darticulier très instrumenté,
- sur des sites anciens ou existants déià recensés et ayant déjà fait l'objet d'une étude.
La première voie implique le choix d'un site, son instrumentation et son suivi, la
seconde pourrait faire l'objet de l'écriture d'une base de données (paramètres nécessaires à la
validation) en collaboration avec un organisme tel que l'Ademe,

* les modèles de transport existants sur le marché sont tous du type déterministe, une
approche alternative se développe actuellement : "lesmodèles stochastiques" où les résultats
ne sont plus donnés comme des valeurs proprement dites mais en terme de probabilité (les
paramètres d'entrée pouvant être également probabilistes ou non numériques du type un peu,
moyen, beaucoup). Quelle sont la complexité et la validité de cette approche ?

* on peut aussi définir un sujet dont l'enjeu parait moindre mais qui mérite cependant
l'attention : les logiciels existant actuellement sur le marché semblent manquer de
convivialité. Il existe une demande dans ce domaine notamment pour la mise au point
d'interfaces et aussi pour une adaptation "locale" simplifiée de logiciels commerciaux,
paramétrés et adaptés à des cas restreints ou localisés (par exemple surveillance d'un site
sensible),

* une étude comparative des qualités et limites de quelques logiciels commerciaux

peut être utile à condition de pouvoir obtenir les éléments objectifs de la validation des
simulations,

* des études particulières, en laboratoire, des mécanismes fondamentaux régissant les

interactions des polluants avec les sols, sont déjà effectuées par des équipes spécialisées. Il
serait souhaitable de recenser les thèmes déjà étudiés afin de cerner les domaines encore peu
connus, et éventuellement, de lancer de nouvelles recherches.

-63-
 Biblioarar,hie

BIBLIOGRAPHIE

M.B. ABBOTT, J.C. BATHURST, J.A. CUNGE, P.E. O’CONNELL, J. RASMUSSEN

(1986) An introduction to The European Hydrological Sysytem - Système Hydrologique
Européen a SHE », 1 : History and Philosophy of a physically-based distributed modelling
system. 2 : Structure of a physically-based distributed modelling system. J. of Hydrology, 87,
p. 45-59 and 61-77
ALEXANDER M. (1994). Biodegradation and bioremediation. Academic Press, 302p.
BOURG A.C.M. & SAUTY J.P. (1987) . Can the distribution coefficient Kd concept
account successfully for the adsorption of polluants in hydrodynamic mass transfer model?
B.R.G.M. Orléans. Vulnérability of soil and groundwater to polluants :proceeding of the
International Conference on March 30, at Noordwijk aa ZEE.
CASTAGNY G. (1982). Principes et méthodes de l’hydrogéologie. Dunod Université236p.
DUPONT A. (1994). Validation d’un modèle mathémathique de transferts de polluants dans
les sols à partir de données d’un modèle réduit. TFE de I’ENTPE, V/A Environnement et
Hydraulique, 7Op.
ENCYCLOPAEDIA UNIVERSALIS. Solutions en équilibre (292-302) et Sols (272-285).
Enseigner la géologie. Collège - Lycée. Nathan, 286p, 196-223
DE MARSILY G. (1983) Hydrologie quantitative. Col1 Sciences de la Terre. Ed Masson, 2 15
P.
GAUDET J.P (1974). Ecoulement transitoire d’eau et de sel dans la zone non saturée. Thèse
3ème cycle mécanique des fluides. Université scientifique et médicales de Grenoble. 136p.
GOUNOT A.M. (199 1). Ecologie microbienne des eaux et sédiments souterrains. hydrologie,
n03, 239-248.
GUINOT V. (1995) . Modélisation mécaniste du devenir des produits phytosanitaires dans
l’environnement souterrain, apllication à la protection des captages en aquifere. Thèse, 2 16p.
H1LLEL.D. (1988). L’eau et le sol. 2ème Ed. Sciences et Techniques.288 p.
HISTOIRE DE LA TERRE NOTRE PLANETE Société Géologique de France.
IN SITU BIOREMEDIATION. When does it work (1993). National Research Council,
207p, 80-9 1.
KATYAL, A.K., KALUARACHCHIi, J.J. and PARKER, J.C. : MOFAT, a two
dimensionnal finite element program for multiphase flow and multicomponent transport,
program documentation and user’s guide, 1991.
LEDOUX E. (1993) . Modélisation du transport des polluants par les eaux souterraines :
Séminaire de programmation dfodéllisation du comportement des polluants dans les
hydrosystèmes».
G. de MARSILY, Soi1 and Rock Flow Properties, Stochastic Description, Encyclopedia of
Physical Science and Technology, vol. 15
G. de MARSILY, Milieux Aléatoires et Phénomènes Hydrologiques, C.R. Acad. Agric. Fr,
1994,80, no 6 , pp. 3-20
G. de MARSILY Quelques réflexions sur l’utilisation des modèles en hydrologie. Revue des
sciences de l’eau, 7 (1994) pp. 219-234
.
MARTELAT A. (1993). Estimation des flux hydriques et nitriques dans les sols agricoles.
Approches spatiale à plusieurs échelles dans la plaine du Rhin Supérieur. Thèse de 3ème
cycle. 18Op.
MOLINARI J.,ROCHON.J. (1976) . Mesure des paramètres de transport de l’eau et des
substances en solution en zone saturée. La houille Blanche, no 3 14,233-242.

-64-
 Bibliographie

NE-ZHENG S. (1996). Mathematical modelling of groundwater pollution, Geological

Publishing House réedité chez Springer, 377, 187-199.
PHILIP J.R (1957) The theory of infiltration : Sorptivity and algebraic infiltration equations.
Soil Science 84.257-264p.
WMBAULT G. (1986). Diffisivité et conductivité hydrauliques de matériaux ou sols non
saturés en eau. Bull. Liason Labo. P. et Ch. 145 Sept-Oct. Réf 3 113.
VAUCLIN M. (1993) Modélisation du transport de solutés dans la zone non saturée du sol :
revue et état de l’art. Séminaire de programmation ccA4odéllisation du comportement des
polluants dans les hydrosystèmes».
M. VAUCLIN, Modélisation du Transport de Soluté dans les Zones no Sauturées du Sol :
Revue et Etat de l’Art Laboratoire d’Etude des Transferts en Hydrologie et Environnement.
VESTED, H.J. et al. : Advection-dispersion modeling in three dimensions. Appl. Math.
Modeling, 16 : 506-519, 1992.
WAGENET RJ. AND HUSTON J.L. (1989) : LEACHM. A process-based mode1 of water
and solute movement, transformations, plant uptake and chemical reactions in the
unsaturated zone. Departement of agronomy, Comelle University-Ithaca, NY 14853
WHITE R. E., S. S. DYSON, R. A., HAIGH, W. A. JURY and G. SPOSITO 1986. A
Transfert Fonction Mode1 of Solute Transport through Soil. 2 - Illustrative applications. Water
Ressour. Res. 22 : 248-254
YONG R.N., MOHAMED A.M.O., WARKENTIN B.P (Ed 1992). Principles of
contaminant transport in soils. Developments in Geotechnical Engineering 73,327, 143-157.

-65-
 Liste des annexes

LISTE DES ANNEXES

Annexe 1 : Décret n”89-3 du 3 janvier 1989.

Annexe 2 : Circulaire de 3 décembre 1993

Annexe3 : Article 34, décret du 21 septembre 1977

Annexe 4 : Article 23 de la loi du 19 juillet 1976

Annexe 5 : Fiches détaillées de divers modèles

Annexe 6 : Tableau récapitulatif des projet retenus pour le XIèmeContrat Plan Etat-Région

Annexe 7 : Programme “Hydrology” de l’université de l’Illinois

Annexe 8 : Programmes de recherche de QUESTOR

-66-
 LISTE DES FIGURES

figure 1 : Cycle des produits phytosanitaires (Vachaud and Vauclin, 1992).

f i m e 2 : Isothermes d’adsorption les plus fréquents.

figure 3 : approche 2-D Z S

figure 4 : approche 3-D ZNS i ZS

figure 5 : approche 3-D complète (DEDALE 3-D)

figure 6 : approche bidimensionnelle

figure 7 : approche tridimensionnelle

figure 8 : Liste non exhaustive de modèles disponibles.

figure 9 : Liste complémentaire des modèles

-67-
 LEXIQUE

Argile
Terme désignant soit un minéral (=minéral argileux), soit une roche composée pour l’essentiel
de ces minéraux.
- Minéraux argileux - Phyllosilicates hydratés, se présentant en très petit cristaux
(quelques pm, en plaquettes hexagonales ou parfois en fibres). Ces minéraux sont très
nombreux, avec principalement : la kaolinite, l’illite, les smectites, les interstratifiés, les
minéraux fibreux (auxquels on ajoute les chlorites et les micas lorsqu’ils se présentent en très
petit cristaux).
- Roches argileuses - Roches sédimentaires ou résiduelles à grain très fin (classe des
lutites), cmtenant au moins 50% de minéraux argileux, auxquels peuvent s’ajouter d’autres
minéraux très divers, détritiques ou non, d’où des compositions très variées (argiles
calcareuses, argiles sableuses, argiles micacées.. .).

Calcaire
Roche sédimentaire carbonatée contenant au 50% de calcite CaCO3, pouvant être
accompagnée d’un peu de dolomite, d’aragonite, de sidérite.

Coefficient d’emrnagasinement S sans dimension

Ce coefficient permet de mesurer en place et sur un volume important, par le biais d’études et
expérimentations, la capacité du reservoir à emmagasiner de l’eau. C’est le rapport du volume
d’eau libérée ou emmagasinée, par unité de surface de l’aquifêre (1 m2) sous l’effet d’une
variation de charge hydraulique correspondante, Ah.

Conductivité hvdrauliaue K en m.s-’ ou coefficient de Derméabilité pour les hvdrogéologues

La conductivité définit le volume d’eau gravitaire en m3 traversant en une unité de temps (s),
sous l’effet d’une unité de gradient hydraulique, une unité de section en m2 orthogonale à la
direction de l’écoulement. Ce paramètre a la dimension d’une vitesse.

Diffusion moléculaire
C’est un phénomène physique lié à l’agitation moléculaire : dans un fluide au repos, le
mouvement brownien envoie des particules dans toutes les directions de l’espace. La
résultante de cette agitation moléculaire sera donc un transfert de particules des zones à
concentration élevée vers les zones à concentration plus faibles.

Difisivité en m2.s-l
Le rapport T/S est appelé la difisivité de l’aquifêre. Elle régit la propagation d’influences
dans l’aquifêre.

Dispersion
C’est un phénomène de mélange essentiellement lié à l’hétérogénéité des vitesses
microscopiques au sein du milieux poreux (cf paragraphe sur la dispersion cinématique, II-
3.2).

-68-
Gravier
Elément de quelques millimètres dans les roches sédimentaires détritiques (classe des rudites).
à 15 ou 30 mm, parfois plus.
Selon les classifications, les limites vont de 1 ou 2 ~llfn

Indice de vides noté e

C’est encore un indice du volume relatif des pores mais il se rapporte au volume solide et non
au volume total de sol. Cet indice est généralement préféré par les ingénieurs et les
mécaniciens des sols, tandis que la porosité est plus utilisée en physique du sol appliquée à
l’agriculture. En général, e varie de 0,3 à 2,O :
e = (vol d’air + vol d’eau) / vol solide
Limon
Dépot détritique meuble, argileux et silteux, à grain très fin (classe des lutites), continental et
d’origine fluviatile, lagunaire ou encore éolienne (limon des plateaux, loess).

Perméabilité
C’est l’aptitude d’un réservoir à conduire l’écoulement de l’eau dans des conditions
hydrodynamiques imposées, soit sous l’effet d’un gradient de charge hydraulique. On la
mesure par deux coefficients : le coefficient de perméabilité ou conductivité hydraulique et la
perméabilité intrinsèque.

Perméabilité intrinsèaue k en m2 ou en Darcy

C’est le volume de liquide en m3 d’unité de viscosité cinématique traversant en une unité de
temps (en s), sous l’effet d’une unité de gradient hydraulique, une unité de section (en m2)
orthogonale à la direction de l’écoulement. Elle s’exprime en m2 ou en Darcy. Le Darcy est la
perméabilité d’un milieu débitant 1 cm3.s-l à travers une surface de 1 cm2 sous un gradient
hydraulique, normal à cette surface, de 1 atmosphère par centimètre.

Permébilité relative
La perméabilité intrinsèque ki est fonction de la saturation du milieu en fluide i : plus la
portion du milieu poreux, occupée par le fluide i, sera grande, plus la perméabilité liée à ce
fluide sera grande. On définit alors la perméabilité relative kfi telle que :
k” = ki / k
k étant la perméabilité intrinsèque (saturée) du milieu poreux.

Porosité cinématique
C’est le rapport de la vitesse de déplacement à la vitesse de filtration. il équivaut au rapport du
volume des vides réellement parcouru par l’eau gravitaire au volume du milieu (saturé ou
non). C’est la teneur en eau mobile. Ce concept est assez proche de la porosité efficace.

Porosité efficace, notée ne en %

C’est le rapport du volume d’eau gravitaire, Ve que le réservoir peut contenir à l’état saturé,
puis libérer sous l’effet d’un égouttage complet, à son volume total Vt.
Porosité efficace, ne = Vol d’eau gravitaire (Ve)/vol totale (Vt).

Porosité totale O ou n en %
La plupart des roches et des sols contiennent naturellement un certain pourcentage de vides
qui peuvent être occupés par de l’eau ou des fluides. C’est ce que l’on appelle la porosité ou
porosité efficace suivant la littérature.
Porosité totale = volume de vides/volume total de roche.

-69-
Sable
Au sens courant, matériau meuble formé de grains de quartz (grains de sable) tel celui des
plages ou des dunes. Plus précisément, sédiment détritique meuble dont les grains sont en
majorité compris entre 1/16 mm (62,5 pm) et 2 mm (classe des arénites).

Teneur en eau 8 en %
8 la teneur en eau volumique exprimée en pourcentage, est la quantité d’eau contenue dans
un reservoir, saturé ou non, rapporté au volume total, quel que soit le type d’eau considérée.
Dans la zone non saturée, on parle d’humidité définit par:
volume d’eau
O=
volume total

définit le volume d’eau gravitaire en m3 traversant en une unité de temps (s), sous l’effet
d’une unité de gradient hydraulique, une unité de section en m2 orthogonale à la direction de
l’écoulement. Ce paramètre a la dimension d’une vitesse.

Transmissivité T en (m2.s-’)
La productivité d’un captage dans un aquifère est fonction de son coefficient de perméabilité,
K et de son épaisseur, b.
On définit la transmissivité comme le débit d’eau qui s’écoule, par unité de largeur, L, d’un
aquifère, sous l’effet d’une unité de gradient hydraulique.
T(m2.s-’) = K(m.s-’).b(m)
On la mesure sur le terrain grâce aux pompages d’essais.

Vitesses d’écoulement en m.s-’ ou m.an-’

L’écoulement des eaux souterraines peut-être considéré comme le déplacement de particules
d’eau dans l’espace et dans le temps, le long de trajectoires, appelées lignes de flux ou lignes
de courant.
Deux méthodes de détermination des vitesses d’écoulement, correspondant à deux
conceptions de leur mouvement, sont utilisées :
- application de la loi de Darcy, avec une correction introduisant la porosité efficace, calculant
la vitesse effective, notée Ve. C’est le concept d’hydrodynamique souterraine,
- opérations de traçages, sur le terrain, mesurant la vitesse de déplacement, notée Vd. c’est le
concept de l’hydrocinématique souterraine.
La vitesse de filtration V ,calculée par la loi de Darcy, se rapporte à la section totale (section
avec trajectoires théoriques qui ne tiennent pas compte des grains et des pores A). Elle n’a pas
de réalité physique. Par exemple pour un débit de nappe Q=lm3.s-’ , traversant une section
totale de 200 OOOm2 :
V = Q / A = 5.10-6m.s-’ = 150dan
La vitesse effective, Ve : la loi de Darcy est comgée par rapport à la section efficace (section
où l’on considère que les lignes de courant suivent un trajet complexe à travers les pores et les
grains) :
Ve = V / ne (de l’ordre de 10-5m.s-’)
La section efficace est plus petite que la section totale, donc pour un débit de nappe constant,
la vitesse effective est plus grande que la vitesse de filtration.
La vitesse de déplacement introduit le concept de dispersion et étudie les trajectoires réelles
dans les vides du milieu à l’échelle microscopique. L’opération de mesure pour se faire en
laboratoire sur colonne ou sur le terrain. On obtient la direction réelle de l’écoulement et les
paramètres de dispersion. On injecte un traceur à un temps t et on mesure sa concentration au

-70-
 Lexiaue

cours du temps. L’intervalle de temps écoulé entre l’entrée et la sortie est appelé temps de
séjour. La courbe obtenue (rapport des concentrations en fonction des temps de séjours)
permet d’obtenir la vitesse de déplacement.

-71-
Eau TEXTES
DéCret no89-3du 3 janvier 1989

ANNEXE 1-3
Exigences de quallté des eaux douceo superficieiles uüll8ées ou d e s t i h s à être UaliSees pour la production d'eau
destinée à la consommation humaine
(Dé- no 91-257 du 7mars 1991, art 2-III)
G :voleur gwde. 1 ; valeur lin& imp4rative.

Al A2 A3
GROUPES DE PARAMÈTRES PARAMhES
G I G 1 7
par am^^ organoleptiques. Coloration (après filtration simple) mg/l
(échelle Pt) ............................................ 10 1 1 20 50 100 150 2ûO

Paramètres physico-chimiques lib à la COnduaMtt p s l h 20 "C .................. 1OOO 1OOO

,amturc natureiie des eaux.
T - P ç F - E pc)..................................... 2 2 2 5 22 25

.......................................
pH (unités pH)

t2hmzqsSm) ................................
- -. . -.

s - -d h - ~ ; .....~..........-..........
.. .

Ma~èreSeas&pension (mg/i)
....-
...............
ZkEk
65-85

25

1
Demande biochimique en oxygène
@BO,) à 20 "C saus nitrification (m@l
O,) ............................................................ c5
Demande chimique en oxygène @CO)
( m a 03 _..............................................
~

I l 30

Paramètres concernant les substances

Taux de saturation en oxygène dissous
(% 09

Nitrates (m@l NO3

....................................................
................................
I l
'70
25 50
> 50

50
> 30

indésirables.
Ammoniaque ( m a m)........... 0,OS 1 15

Azote Kjeldhal (NO,) excepté (maN)

Hydrocarbures dissous ou émulsionnés

après extraction par éther de pétrole
(ma) ...................................................... 0.2
~~~ ~

Phénols (indice phénol) para-nitraniline

4 aminoantipyrine (rngll C,H,OH) ........ 0.001 0.001 0.005 0,Ol 1 0.1

Agents de surface réagissant au bleu de

méthylène (mg/i lauryl-sulfate) ............. 0.2 0.2 0.5
~~~~~ ~ ~

Substances extractibles au chloroforme

( m a ) ......................................................... 0,l 0.2 0.5

Fer dissous (mg/l Fe) ............................. 0.1 0.3 1 2

Manganèse (mg/l Mn) .................

.....__._. ~ 0.05 0.1

Cuivre (mgll Cu) ..................................... 0.02 0.05 0.05

Zinc (mg4 Zn) ............................................ 0,5 3 1 5

Phosphore (mgll P, O,) ......................... 0.4 0,7

~ ~~~~~~

~~~~
Fluor (mg/l ............................................

Bore (mgll B)...........................................

Baryum (mg/i Ba)

~~ ~

...................................
0.7-1

1
15

0.1
0.7-1.7

1
4-
0.7- 1.7

2858 Feuillets 98 (15 juin 1995) Code Permanent Environnement et Nuisances

Pour une information plus récente reportez-vous a la table d'actualisation, en tête du Code
TEXTES Eau
D&ret no89-3 du 3 janvier 1989

A l A2 A3
GROUPES M PARAInÈTRES PARAYmES ,

G I G I G I

paramttrcs CoIlaRmUt les subatanas

toxiques.

1 Arsenic (pg/l As) ....................................

(ps/i cd)................................

cyanmcsoig/t CN-)...............................

Plomb (pg/ï Pb) .....................................

..................................
Hydrocarbura polycydiques m-ti-
q~&tal6substan<sviséeSpardLcrrt
du 3janvier 1989 (&.-
l).......................
Pesticides. Totai :parathion,HCH, dieldrine (&).

Paamttres microbioiogiqucs. A_ Fsdisei;slestotaux 37 *c(100 mi) ..........

thcrmotolérants (100 mi) .....
.-
--strepuicoques fécaux (100 mi) ...............
salmonelles ............................................

ANNEXE II
(Dé-. no 90-330 du IO avr. 1990, art. 12-1)
Programmesd'analyse der écbnülion8 d'eau

La présente annexe fixe la programmes d'analyse d'échantillons,pour l a eaux disuibuées par un réseau collectif public ou pnVi (1), les eaux
conditionncts et la giaa alimentaire (II) et les eaux utilisées dans l a industries alimentaim (III) ainsi que la modalités d'adaptation du
programme d'analyse 0.

1. - EAUX DISTRIBUÉES PAR UN RÉSEAU COLLECTIF PUBLIC OU PRIVÉ

TABLEAU 1
Anai.vses bacrériologiques

ANALYSES B A ~ R I O L O G I Q U E S

Réduite (B 1) !Sommaire (B 2) Complbte (B 3)

Streptocoques fécaux. Streptocoques fécaux. Streptocoqua fécaux.

ColiforneS.

I ria aérobies revivifiables à

22 OC et 37 OC. I ries aérobies rcvivifiables à
22 OC et 37 OC.
Spores de bactéries anaéro-
bies suifitoréductrica.

Code Permanent Environnement et Nuisances Feuillets 81 (15 juin 1991) 2859

Pour une information plus récente reportez-vous à la Table Récapitulative du bulletin
 Eau TEXTES
,
Decret no893 du 3 janvier 1989

B. - FrCqiwre da prCWlcmtiib d'mi vements d'échantillons d'eau à effectuer analyses a effectuer doit être au moins égai à
i d p w chaque année dans la rcs~oura(RP. RS) et celui correspondant à la popuiatim des
dans I'cau avant son refoulement dans le coormunes dcsseMes par le ristau sans arc
(D&.n"91-257Chc 7-s 1991.~t.4-1) riSeau de distribution (P 1. P2 P. P 2 S, P 3) inférieur au nombre dcs communes desscr-
selon le débit journaiia de I'eau. vies ;
l0 Le tableau no 3 précise le type d'analyses
a effectuer selon que les échantillons d'eau Le tabkau 5 indique la fréqucna de ces 20 Lorsque la valeur de artaias paramétres
ont été préicvés dans la res~ourcc(R),dans prélèvements dans Peau distribuée aux figurant dans Ics analycs-typcs R ct P est
I'eau, tmitéc ou non, avant son refoulement consommateurs selon la population dcsscMe susxptibie d'erre modif'kcm cours de distri-
dans le riseau de distribution (P) ou dans le par le réseauet selon que I'eau est dcsinfecta bution, les prélèvements d'échantillons d'eau
réseau 0). ou non. Lorsqu'un k u de distribution p e ü m t êUe rtalisés h le rcsCaü de distri-
u Le tableau 4 indiquela friqucna des préik- dessert piusieurs commu~s, le nombre des bution. N

C4r Clb C4C

~~ ~~~ ~ ~

Paramcm or- - AsPtct(qd- -Aspaa (quali- - AWxt (q=-

titatif): odeur,
plloleptiques tatrf): odeur, tatif): odeur,
saveur, cou- saveur, cou- saveur, cou-
leur. leur. leur.
- Turbidité. - Turbidite. - Turbidité.

Paramètres - pH. - Température. - Température. Matières totales

physico- - Conductivité. - pH. - pH. en suspengon
:himiqueS. - Conductivité. - Conductivité. (M.E.S.).
- chlom. M a d e chi-
Structure na- - Sulfates. mique en oxy-
turelle des - silice. giae (D.C.O.).
zux. - calcium. Demande bio-
- Magnésium. chimique en
- Sodium. oxyginedissous
- Potassium. @.B.OF)
- Aluminium.
- Rtsidus KCS.
- Oxygtne dis-
sous.
- Anhydride
carbonique
libre (essai au
marbre) ou
calcul de 1%-
quilibre calco-
carbonique.
- Carbonates.
- Hydrogéno-
carbonates.

2860 Feuillets 81 (15 juin 1991) Code Permanent Environnementet Nuisances

Pour une information plus récente reportez-vous A la Table Récapitulative du bulletin
TEXTES Eau
Décret no 8 9 3 du 3 janvier 1989

ANALYSE!3 PHYSICO-CHIMIQUES
I I
Anillu- -Pw-
ddnilqrn- dJmipin-
(a)
ai
?aramètres (ficr.n"91-257du - Nitrates. - Azote.
mncemant les ~ W S m. 3- --
1991, Nitrites. Kjeldhal.
;ubstanca in- P) - Ammonium. - Hydr~~ar-
iéyrables. - Nitrates. - Oxydabilité bures dis-
- Nitrites. a u KMnO,, sous.
- AmmOniUm.~ à chaud, CII - Agents de
-DeuxParamé- milieuacide surfaa.
trcs parmi les - Hydrogéne - indice.
pilramètrrssui- sulfuré. Phénol.
~ a n t s :chIoni- - Fa.
res, sulfates; - Cuivre.
o x y d a b i l i t i ~ - zinc.
K M n 0 y - o ~ - Manganèse.
carbone mga- - Phosphore.
nique t o t a l , - Ruor.
titre alcaiimé-
vique complet
OU dUrC3 t ~ t a -
____- - Ir
-
€%
- &F
y-r --- - Chlore rési-
+- - -+-*d+
a ou-tout
w- duci ou tout
pa€~-è€ïere autre parame-
9 ~ r e s ë ë n i presentatifT trc représen-
tatifdu traite- traitewrrt-àt tatifdu traitc-
ment de dé- désinfecton. ment de de-
sinfection. sinîcction.
'aramètres - cadmium. - Arsmic.
onamant l a - Plomb. - cyanures.
ubstances - H.P.A. - Chrome.
oxiqnes. - MIXCU~C.
- Sélénium.
4utres pa- - Pesticides.
amètres. - Composés
organoha-
logénés
volatils.

TABLEAU 3
Anaiyses !yps
iDécr. no 91-257 du 7 mors 1991.art. 4-II)

EMPLACEMENT RESSOURC€ PRODUCTION DlSTRlüUfION

Au point de puisage, Aprds traiternent et avant refoulement

avant traitement ou au point de puisage en i'absenw de traitement

Eaux
(R.P.) (R.S.) (P 1) souterraines
et superilcielles
Eaux Eaux Eaux Eaux (Dl
Eaux Eaux Bouterraines souterraines superficielles souterraines
Origine de l'eau
routanaines supeifkielies et eaux et eaux
supciwcielles (p 2 s) supemcieiia

ANA LYSES
TYPES

Code Permanent Environnement et Nuisances Feuillets 81 (15juin 1991) 2861

Pour une information plus rkente reportez-vous a la Table Récapitulative du bulletin
Eau TEXTES
W r e t no89-3 du 3 janvier 1989

TABLEAU 4
Fréquences annuelies d'analyse (èchaniülonrprélevés à Io reJsoi<rccer en usine)
(Dècr. no 91-257 du 7 mars 1991, mi. 4-III)

&BIT JOURNAUER FREOUENCES ANNUELLES DlÉCHANllUONNAGE

(mV-
R.P. R.S. P1 P2P P2S P3

Inférieur à 100 ....................................... - - 1 - - -

De 100 à 399 .......................................... 112 2 2.5 1P 1 115
De400à999 .......................................... 112 2 23 1/2 1 115
De 1 O00 à 1 999.................................... 1P 2 33 112 1 115
De2O00à5999 .................................... 1 3 7 1 1 1P
Dc6ûûOP9999 .................................... -2 6 8 1 1 1P
De 1OOOO à 19 999 .............-...-.-...-..... - 2, 6 14 2 2 1
De 20000 à 29 999 ................................. 4 12 22 3 3 1
D e 3 0 0 0 0 à 5 9 9 9 9 ................................. 4 12 42 6 6 1
- .- -- - .
De60ûûOà99999 -...---.---. z 12 70 10 10 1 '
suptmur ou égal à 100 O00 ................-. e-- J2 140 20 20 1
i

TABLEAU 5
Fréquences annuelles d'analyses
(Échantillonsprkltvés en distribution)
(Dècr. n" 91-257 du 7mars 1991. an. 4-IV)

D
POPULAVON DESSERVIE

~ ~ _ _ _ _

-
~~ ~~~

500 habitants ...................................................................................................... 7 4

2 O00 habitants ...................................................................................................... 6

5 O00 habitants...................................................................................................... 12
10 O00 habitants ......... ............. 24 24
30 O00 habitants ...................................................................................................... 60
50 O00 habitants .......................... 90
100 O00 habitants .......................... ........................ 150 240
150 O00 habitants ...................................................................................................... 210
300 O00 habitants ...................................................................................................... 390 720

Pour la populations inférieures a 500 habitants. le nombre d'analyses D est égal a 2 dans le cas d'eaux non désinfectéeset a 4 dans le cas d'eaux disinfectées.
Pour l a populations supérieures a 500 habitants. le nombre d'analyses a effectuer est obtenu par interpolation linCaire entre l a chiffres ii's dans les
colonnes D, le chiffre étant arrondi a la valeur entiére la plus proche.
Pour les populations supérieures a 300 000 habitants. le nombre d'analyses a effectuer est obtenu par extrapolation linéaire. le chiffre étant arrondi a la valeur
entière la plus proche.

2862 Feuillets 81 (15 juin 1991) Code Permanent Environnement et Nuisances

Pour une information plus recente reportez-vous à la Table Flécapitulative du bulletin
TEXTES Eau
Décret no893 du 3 janvier î9û9

II. - EAUX CONDITIONNÉES ET GLACE ALIMENTAIRE

JAS Irialvpn et fréquences annuelles d'échantillonnage doivent respecter les prrsCriptiOnS
d t f i n i c s au tabiean no 6 ci-après :

F
ANALYSES TYPES

Rasoiirct (1). 3

Analyse bactériologique 3
COmplitc (B 3)
Anaiysc physico-cibique
CcJmpEtC (C 3).
__
E A n a l y physicctchimi- 1
A ques spcciaies (C4a, C 4 b
et c 4 c).

Après conditionncwnt Analyse bactériologiquc 6

ou après mbaiiagc. cornpletc (B3).
Analyse physico-chjmique
sommaire (C 2).
1

Dans le cas ou i'eau est traitée avant conditiomcment, le programme riglementairc

d'analyses est adapté a ia nature et aux caracttristiques du traitement effectué.

M Flacons vides aprèâ Analyse ùactériologique 6

A iavagc, désinfection et sommairc.(B2).
T rinçage.
l? Surfacc des bouchons,
R capsuies ou joints en
1 contact avec i'eau.
E
L

Eau de rinçage (2).

(1) Les analyses sont à effectuer lorsque i'cau ne provient pas d'une distribution publique.
(2) Lorsque ie rinçage n'est pas pratique avec ïeau a conditionner ou avec une eau de distribution
pubtique.

Code Permanent Environnement et Nuisances Feuillets 82 (15 septembre 1991) 2863

Pour une information plus r k e n t e reportez-vous a la Table Récapitulative du bulletin
Eau TEXTES
Décret no 89-3 du 3 janvier 1989

I
Analyses types. Analyse bactériolo@que Analyse physico-chimique
cornplite @ 3). compiite (C 3).
Analyse physico-chimique Anaiyses physiw-cùimiques
sommaire (C 2) particulières ( C 4 k C 4 b ,
c 4 c).
Fréquenas annueiics Les fréquences annueiics d'ichantitlonnagc sont identiques à
d'échantillOMl3gc. ceiies p r h pour les anai- effectuées sur les ressources
aiimentant lcs réseaux d'adduction (tabieau 4) ;clics ne peuvent
être inférie- a cdlcs indiquées cidessous :

IV. - ADAPTA'ITON 6. Plomb ............................ 50 PeJi (Pb)

DU PROGRAMME 7. Sélénium ........................ 10 Ppii (W
Pour l'appiication de l'article 9 du présent Pour ks strbrt_eoas suivantes, valeurs des 8.Pesticides pour le
k e t , le programme d'analyse peut être concentrations supérieures aux valeun ci- total des trois SU&-
modifié dans les conditions suivantes : aprts : tanas suivantes
(parathion. HCH,
1. A coût égai de réalisation du programme
d'analyses, la nature et la fréquence des para- dieidrine) ........................... 5 Pg/i
mitres analysés peuvent être modifiées.
1. Nitrates ................ 9. Hydrocarbures po-
2. a) Pour les eaux d'origine superfcielle, les lycycliques aroma-
fréquences d'içhantillonnage peuvent étre tiques :
réduites d'un facteur 2 à i'exccption de celles Pour le total des six
concernant les analyses bactériologiques ; substanm suivan-
b ) Pour les eaux souterraines ou profondes, 2. Ammonium .......... tes ..................................... 1 P&'I
les fréquences d'échantiiionnage peuvent être
réduites d'un facteur 4.
3. Oxydabilité Ruoranthine ;
(KMnO,) en milieu Bemo (3.4) fluoranthme ;
3. Le programme peut comprendre des acide ........................ 10 mp/i (02) Benzo (11.12) fluoranthme;
analyses supplhentairrs dont le coût ne Benu> (3.4) pyrène ;
conduit pas à un dépassement suptrieur à 4. Phénols (indice
phénols) para nitrani- Benzo (1,12) pérylène;
20p. 100 du programme d'analyse établi Indéno (1,2,3-cd) pyréne.
selon les modalités prévues aux tableaux 1 et iine et 4 aminoantipy-
2 de la présente annexe. nne ..........................
5. - Parrmétres microbiologiqoes
ANNEXE III 5.Agents de surface Eau contenant pius de 20000 coliformes
(réagissant au bleu de thermotolérants et plus de 10 000 streptoco-
IDécr. no 90-330 du 10 avr. 1990. méthylène ................ 0.5 mg/1 ques fécaux par 100 millilitres d'eau prelevte.
art. 13.1 et II et Décr. no 9 1 - 3 7 (lauryl-sulfate).
du 7 m r s 1991. art. 5 ) 6. Hydrocarbures dis-
sous émulsionnés
CIRCULAIRE DU 16 MAI 1989
Limites de quaüté des eaux brutes après extraction au relative & la présence
utilûées pour In productioa d'eau tétrachlorure de car- de composés organohalagénis volatils
destinée a In coisommation humaine bone ......................... 1 mgjl dans les eaux destin-
& la consommation humaine
7. Zinc ...................... 5 mgfl (Zn)
1. Paramètre orgaaoleptiqw (J.O. du 28 juin 1989)
8. Baryum ............. 1 mg.1 (Ba) pour Destinataires :.Uessirurs les préfets de région
Coloration après filtration dépassant 200 les eaux superfi- direction régionale des affaires sanitaires et
mg 1 de platine en référence à l'échelle plati- cielles ». sociales) et à Madame et Messieurs les préfets
ne. cobalt. de département (direction départementale des
4. Paramètres concernant affaires sanitaires et sociales).
2. Paramitres en relation des substances toxiques Il y a une quinzaine d'anntes. a été mise en
avec la striichue naturelle des eaux évidence la présence de certains composés
Pour les substances suivantes. valeurs de organohalogénes dans l a eaux traitées desti-
1. Température de l'eau supérieure à 25 OC. concentrations supérieures aux valeurs ci- nées à la consommation humaine. Peu aprtk.
2. Pour les substances suivantes. valeurs des après : du trichlorethylène était identifie dans des
concentrations supérieures aux valeurs ci- eaux de nappes.
après : 1. Arsenic .......................... Depuis lors. de nombreux travaux ont été
100 ~g 1 (As) menés dans différents pays. notamment en
a) Chlorures .................... (< 200 n mg/l (Ci)
2. Cadmium ...................... 5 pg,l (Cd) France. pour mieux connaître ces produits.
b ) Sulfates ........................ 250 mgjï (SO,) 3. Cyanures ....................... 50 pg;l (CN) les conditions de leur présence dans les eaux
3. Pour les eaux superficielles. pourcentage et de leur évolution. et pour apprécier les
4. Chrome total ................. 50 pg,l (Cr) risques sanitaires qu'ils représentent.
d'oxygéne dissous inférieur à 30 p. 100 de la
valeur de saturation. 5. Mercure ......................... 1 pg 1 (Hg) La directive C.E.E. no 80-778 du 15 juillet

2864 Feuillets 82 (15septembre 1991) Code Permanent Environnement et Nuisances

Pour une information plus récente reportez-vous A la Table Récapitulative du bulletin
 Le Ministre de l'Environnement

Messieurs les Préfets de Région

>lesdames et Messieun les Préfets de Dépanernent

Qbiet : Politique ae réhabilitation et de traitement des sites et sols poilues.

La prévenzion des risques indunriels et de la poilution des e a u ainsi que la

réduaion des rejets atmosphériques font l'objet depuis de nombreuses années d'actions
cohërentes et concenees. Celles-ci ont déjà conduit dans le domaine indusuiel a des résultats
notables. L'impact des activités humaines sur l'environnement se réduit unsi-Maïs les traces
des pollutions passées en sonr malheureusement d'autant plus perceptibles- L'histoire nous a en
efiet Iegue des pollutions hinoriques, constituées d'anciens dépôts dedechets, de-sols e~ d'eaux
souterraines pollués dont la mémoire collective a parfois deja oublié l'existence - er .la
localisarion précise.

Les c q de pollution les plus,flqants Ont fait l'objet par le passé d'adons de
résorption efficaces mais ponctueiles, Vos services etablisent pinodiquement lci\.enraire des
<._.. .
cas COMUS et la circulaire du 9 janvier 1989 détaille, pour le cas ou aucun respo&blesolkxble
- - .*-.--L:.-
n'est identifié, les modaiités d'exécution d'office de travaux par I'asence de'ïemironnement et
de la m-uise de ïenerge-

Des mesures de prévention dans l'aménagement et I'exploitation des indlationr

potentidement a l'orignede telles pollutions ont paraliëlement été mises en oeuwe depuis de-
nombreuses années, Les stockages de substances dansereuses sont par*-exemple-
systématiquement réa!isés sur.des. aires d e rétention et les cuves souterraines-comportent
cou&ment des enveloppes dédoublées. Les conditions d'élimination des déchets se soni'a;ssi
considthblement améliorées depuis v i n s ans.
 Annexe 7 - Aspects-réglementaires

La loi no 93-646 du 13 juillet 1992 rdative a I'elimination des dëchets et aux

installations ciassées pour la protection de I'environnement et I'arr&ë ministériel du ler mars
relatif aux prelevements et a la consommation d'eau ainsi qu'aux rejets de toute nature des
installations ciassées soumises a autorisation ont récemment encore renforcé ces mesures
préventives.

ïi &ut cependant définir aujourd'hui une véritable politique de traitement des sites
et sois pollués. L'exïaence de ces sites doit être admise et nous devons apprendre a les traiter
dans la transparence et avec sérénité. A l'époque ou seuls les cas les plus graves attiraient notre
attention, l'appellation de "points noin" etait peut-ètre fondée. Aujourd'hui. demx les
quelques milliers de sites potentieilement pollués et présentant des risques ués miables que
compte sans doute notre pays, il convient d'en parier ouvertement mais sans catastrophisme.

1) Les principes d'une politique réaliste de traitement des sites et sols

polluési

Une poiitique efficace mais raisonnée de traitement des sites er sois pollues doit
s'appuyer tout a la fois sur une recherche qstematique et orsanisée des sites potenllellement
concernés et permettre la définition concertée de priorités. Le traitement de chaque site doit
dépendre de son impact ef€ecùf sur I'entironnement et de I'usage auquel il est destiné. Il faut
.donc commencer par ékaiuer pfëasément ies conséquences actuelles et potentielles des
pollutions constatées. C'est l'appréciation de I'impacr sur l'environnement et la santé humaine
qui doit ensuite conduire au choix de techniques de traiternent et d'objectifs de dépollution
adaptés au devenir du site.

Cette démarche générale doit concerner de tres nombreux sites : des secteurs
industriels entiers ont pu e u e a l'origine de pollutiodplus ou moins étendues des sois: Lés
.--.
d-ts
~
anciens d'hydrocarbures ou de produits chiirniques ont en gënéral conduit a des
pollutions des sols. Mëme les stations senice. voire les cuves de fioul domestique si
communes, peuvent avoir conduit à une Iégke et trts locale contamination des sois par les
hydrocarbures. Les décharges municipales de nos villes n'ont sans doute pas non plus toujours
été aménagées et exploitées sans incidences pour I'entironnement-

Vous comprenez donc qulil s'agt d'une action de longue haieine. à i'ëcheiie d'un
siècle et demi d'histoire indusvielle dans noue pays, La mise en oeuvre de cette politique ne
pourra être que progressive et fonction des moyens publics-et privés qui1 sera possible de
mobiiiser. ii convient cependant de s'y consacrer d& maintenant avec un souci naturel
d'identifier des priorités.

36 Gestion des sites (potentiellement) poilu& -

Version O
0 BRGM Ediuons - Décembre 7995
 Annexe 1 - Aspects réglementaires

II) Les moyens et les outils nécessaires ii h mise en oeuvre de cette

poiitique.

J'entends alors progessivement mettre en place les éiements administratifs et

techniques qui vous seront a cet effet nécessaires :

ï) La recherche des sites et sois poilues.

L'identification des sites pollués repose. d'une pan. sur les constats de Ihspection
des installations c l a s s k ou sur des campagnes de détection (étude deS.sols) sur les sites
actuellement en activité et, d'autre pan, sur la coilecte de nombreuses informations relatives
aux activitês industrielles passées d'une région (archives, témoiyases etc.).

.Mes services préparent une grille d'orientation des recherches précisant les activités
ayant potentiellement pu conduire a une pollution des sols ainsi qu'un guide simplifié de
réaiisation d'études de sols. A I'imase des etudes-déchets. ces études porteront pour les
secteurs les plus concernés sur une appréciation des caracténniques des sols et sous-sols des
principaux sites.

Parallelement. des inventaires dits "historiques", fondis sur I'e..amen d'archives tant
privées que publiques et le recueil de témoigazes pourront contribuer a localiser les sites
potentiellement pollues par des activiites aujourd'hui arrêtées.

3 ) L'énluation des risques et de la xulnerabilirë de chaque site.

a) Ebaluation initiale simplifiee et hiérarchisation.

,
En complément a la définition d'un pide methodoIosique d'étude des sois, je
compte pouvoir bientôt vous diffuser une @lie simplifiée d'kduation et de hiérarchisation des
sires permettant d'apprécier rapidement l'ursence d'investi~ationscomplémentaires. Cene @Ie.
d'ëvaiuation simplifiée comprendra des informations de référence permettant d'apprécier pour
cenaines nibnances des teneurs jugées anormalement élevées dans ies sols. Ces seuils
poneront notamment sur les principaux poIIuants rencontrés et sur les niveaux de
contamination qui rendent des investigations complémentaires souhaitables.

Les démarches de recherche de site e: d'ébduation devront être étroitement liées.

L'identification d'un site potentiellement pollué devra s'accompagner Ie plus rapidemenr
possible de i'éwduation simpiifiéc des Nques qu5i ensendre. Cette première évaiuation devra
permettre d'apprécier l'opportunité et I'urgence d'invenigations complémentaires préaiables a
un traitement, Les rkques présentés par le site seront a cet efftt comparés it une échelIe de
gravité qui facilitera au niveau national l'établissement de priorités d'intervention Cette
démarche opérationnelle aboutira a une hiérarchisation de la pollution des sires.
..Y..

Gestion des sites (potentiellementlpollués - Version O 37

BRGM Editions - Décembre 7995
 Annexe 7 - Aspects réglementaires

ü n site sera considéri comme pollue lorsque cette première ésaluarion simplifiee
des risques conduira a prescrire une étude plus approfondie.

b) Etudes d'impact et objectifde rihabilitation.

Des études d'impact approfondies devront être réalisées sur les sites identifies
comme pollués afin de définir les travaux a mener ainsi que les objectifs de contamination
résiduelle a retenir en fonction de l ' u q e ulterieur du site. Dans certains cas. les interaitions
ainsi difinies pourront rendre le site de nouveau propre a tout usase. Dans d'autres cas. lei
travaux permettront simplement de limiter a son plus faible niveau techniquement réalisable et
économiquement acceptable l'impact cr& par le site sur son environnement. Il conviendra alors
de veiller durablement a son affeaation

Un guide methodologque à venir dcfinira les grandes lignes de ces études

approfondies et iilunrera ieur utilité pour &.er des seuils de décontamination a travers
quelques scénarii rypes. Parallèlement, des travaux complémentaires seront engazes afin
d'étudier I'opportunité d'un inventaire géochimique national sur la teneur habituelle des sols et
des roches naturelies en substances minérales. En liaison avec I'évolution de la résiementation
de l'urbanisme, j e proposerai enfin a mes colltipes charges de la santé, de I'équipement et de
l'agriculture, d'envisager la détermination de seuils d'usage banaiise qui fixeraient au cas par w
et pour des scenarii de réfërence des objectifs de dépollution en fonction de i'usage ulrérieur
souhaité.

Tous ces outils administratifs et techniques de\-ront être disponibles dans ieur
version initiale aLant la fin 1994 afin de lancer dés que possible la réalisation d'études
historiques résonales er d'ëtudes des sols dans les secteurs d'actiiites a prion potenrielicmtnt a
I'orisine de telles poliutions des sols.Je compte pour ce-faire organiser une larse concenarion
nationale avec l'ensemble des aaeurs concernes et crëer un groupe de trabail national sur ce
-- sujet. Cene structure Iteere-etouverte participera à la définition technique de tous les éléments
~

d ecette nouvelle poiitique.

3) La création de 'Rructures adaptées d'information et de concertation.

La d a o n à I'édide nationale de: méthodologies et d'outils d'appréciation de

référence CS indispuisabie, La-notionmême de site poiiué est cependant aujourd'hui encore
tellement subjective a p . a d m i s e qu'il convient de î'txptciia localement et de faciiiter tous
les échanges driormation sur ce sujetscnsiile- Lest pourquoi, il est nécessaire d'instituer sur
ce sujet un échange de vues et une concenation à I'èchelle régionaie.
.../...

38 Gestion des sites (potentiellementl poilués - Version O

.-
. Edirjons - Oécembfe 1995
BeGM
 Annexe i -Aspects réglementaires

Cette réflexion pourrait ëtre menée au niveau regional en liaison avec I'ëlaboration
du plan résionai d'élimination des dkhets industriels speciaux- La complexité du travail à
mener dans ce dernier cadre. nocarriment pour ouwir un nouveau centre de stockase de
déchets industriels ultimes dewait cependant conduire a privilé+r la création d'une
commission distincte pour les sites et sols industriels polIués. Ctrtains secrétariats permanents
pour la prévention des pollutions industrielles (S.P.P.P.I.) pourraient dors s'averes des
structures d'accueil appropriées.

Aussi je demande aux préfets de résion de bien vouloir consulter des aujourd'hui
leurs collëgues de dépanement et I'ensembie des panenaires concernés sur la stmcture qui leur
semble la plus appropriée pour traiter de ces questions. ïinteriocuteur naturel de Etat dans ce
domaine doit être la r+on et je demande ainsi aux préfets de regon de conwiter de manière
privilé3ite les présidents des conseils régionaux sur les modalités de création d'une structure de
concerration et d'information régionale sur les sites et sols pollués. Il appanient aux prife;s de
résion de coordonner les auions de lZcat dans ce domine.

Ces stmctures permettront denzaser dès a présent les consultations préalables au

lancement d'erudes hisroriques résionales. Celles-ci doivent en effet reposer sur des actions
concenées.des industriels et de leurs fédérations professionnelles, des élus. de l'administration
et des associations. La reaiisation d'une telle étude pourra donc utilement regrouper l'ensemble
des partenaires locaux concernes autour de la direction rtigionale de l'industrie. de la recherche
et d e i'emironnement ex d e ta delégation régionale de I'agence de I'environnement et de la
maitnse de I'ener@e. Cette même stmcture sera enfin panidierement appropriée. pour
consijner les informations issues de I'ensemble- des sources de - renseignement disponibles
(inspection des innailations classées. industriels. témoignages, etc.).

Je demande aux préfets de région de m'informer akant le ler mars 1994 des
dispositions qu'iis auront retenues pour instiruer cette nécessaire concenation a l'échelle
résionale sur le traitement des sites et des sols pollués. Il conviendra notamment d'étudier le
plus rapidement possible, et dans le cadre d'action collective. ia réalisation d'itudes historiques
régionales sur les sites potentiellement pollués par les acütités aujourd'hui disparues,

4) Le fichier national des sites et sols industriels p01lu&~

Le reculsement des sites potentiellement- pollués devra rester du ressort de

l'inspection des instdlations classées à i'échdon l d . il me s&bIe cep&dmi indispensable de
disposer d'un véritable inventaire national des sites effktivement pollués, où des actions de
dépollution doivent être entreprises. éventueIlement avec Ië concours d e la puissance publiqÜe
si aucun responsable solvable ne peut être identifié. L'existence de-certains sites dewa par
ailIeurs être consewëe en mémoire. notamment si leur état conduit a une restriction de I'usase
ultérieur du sol et du sous-sol.
.../...

Gesrion des sires (porentiellementl pollués - Version O 39

BRGM Edirions - Décembre 7 995
 Annexe I - Aspecrs réglementaires

La direction de la prévention des pollutions et des risques conscituera donc un

fichier national des sites et sols pollu& dont un inventaire.actualisé sera pinodiquement publie-
Aprk un recensement initial des sites COMUS, toute nouvdle découverte d'un site
effectivement pollué. c'est-a-dire dont I'Ctaluation simplifiée des risques conduira a prescrire
une etude plus approfondie, donnera lieu a une inscription sur ce ficfiier. Lorsque les actions de
d&pollutions pratiquées permettront un usase banalisé du sitec il sera raye de l'inventaire
nati0~1.

Iïï) Qualification initiale de l'inventaire des sites connus.

Dans la phase initiale de préparation des éléments techniques et administratifs de
cette nouvelle politique. il ne sera pas immédiatement possible de rialiser une évaluation'
simplifiée et une hiérarchisation des sites connus ou qui seront découvens prochainement. La
mise au point d'une methodoIogie nationale harmonisée d'évaluation simplifiée des risques et de
hiérarchisation des sites ne pourra en effet pas aboutir avant mi 1991.Il convient en effet de
tester la pertinence d'une telle méthode avant de la eeneraiiser.

11 est cependant dés a présent souhaitable de rassembler les informations objectives

disponibles sur chaquë site déjà connu. Je vous adresserai donc tres prochainement un
formulaire de recensement des sites pollues connus et de qualification des informations
disponibles sur chacun d'entre eux Il sera accompagne d'un mode d'emploi destine a
harmoniser le recueil des informations aujourd'hui disponibies sur ces sites, II conkiendra
d'étudier avec soin les sites recenses par le passe et qui ont deja pu fare l'objet de traitemënt
pius ou moins poussés. rattacherai en effet une imponance toute-paAicuiière ce que
ïexistence de ces sites et les restrictions d'usage qui en découlent éventuelIement ne soient pas
prooressivement oubiiêes.

Le fichier national a venir sera donc- initialement renseigne par ce recueil de

données sur les sites déjà connus et complété ultérieurement. Le reFoupement en une meme
base de donnies des informations objectives précises sur les sites CO MUS.^^ d& resultats des
campapzs eqérimentaies de hiérarchisation pennetrra en effet d'affiner la méthode utilisée et
d'établir une éthelie de Zravite de la pollution des sites. Les sites dëcouvens ensuite, a I h u e
de recherches historiques ou de la réalisation d'études des sols, seront -enreginris au fichier
national avec une classification immédiate sur l'échelle de gravité dessites pollués.

Vous le comprenez : i'objeaif est ambitieux et la tâche importante. L a polhtion

des sols et des eaux souterraines peut conduire à des risques réels Po& la-santéhumaine. C m
aujourd'hui bien souvent ia polIution d'une nappe exploitée voire d'un captage qui conduit %la
découverte d'un site pollue. Les risques a les coûts de teiles découvertes tardives sont
considirables :nouvelle source d'aiimentation en eau potable, dérivation de cours d'eau, gel de
milliers d'hectares, etc. L'identification et le traitement des cas les plus menaçants constitueront
ainsi de véritables actions de prévention.
.../...

40 Gestion des sires (potenrielfement] pollués - Version 0

-
?.BRGM Editions Décembre 1995
 Annexe 1 - Aspects réglementaires

Je ne doute pas des effons que vos services et votre inspection des installations
classées entreprendront dans le cadre de la mise en place progressive de cette nouveile
politique de réhabilitation des sites et sols pollués que je souhaite initier.

Je vous demande dans un premier temps de bien vouloir informer le conseil

départemental d'hygiène de votre département de ces orientations jénérales et d'entreprendre
sans tarder les consultations nécessaires à la constitution de.ia Structure d'information et de
concertation sur les sites et sols pollués adaptée au contexte local de votre région. Comme je
vous i'ai précisé plus haut je vous demanderai ensuite de consigier les informations disponibles
sur les sites actuellement connus. Chaque étape de la mise en place de cette nouvelle politique
fera enfin l'objet d'instructions spécifiques détailIëes. Vous voudrez bien me-tenir précisément
informé des conditions de leur mise en oeuvre et des difficultés éventuelles que vous poumez
rencontrer dans leur application.

3Iichel BARNER

Gestion des sites (potentiellement)polluk - Version O 41

-
BRGM Editions Décembre 7995
ANNEXE 3
 Annexe 7 - Aspects réglementaires

DECRET No 77-1133DU 21 SEPTEMBRE 1977

pris pour l'application de la loi no 7 - du 19 juiiiet 1976

d a t i v e auiin-3tiüations d a s s t e s - ~ - a r l a ~ i t ~ ï o ü ' d ë
l'environnement

Sur le rapport du garde des sceaux, minisuc de ia justice, du minian de l'intérieur, du minime de ia défense, &
ministre de ia cuiture et de l'emiironntmem, du minisat déltgué à l'économie et aw finances, ctu minisae de
l'équpement et de l'aménagement du ttrritoire. du minisa de l'agridture, du miaisue de l'industrie, &
mmmmetde I'aflsanaf du minis&e&travaiï et du minimede iasautéetdeia~tésociaie,

Vu la directive 85B37ICEE du conseil des communautés earopocnneS en date du 27 juin 1985 conarnant
i'évaiuarion des incidences de certains projets publics et privés SUT l'environnement;

Vu le code pénai, et notamment son article R 25 ;

Vu la loi no 76-663 du 19juiiiet 1976 reiative aux imaiiations classéts pour ia proteaion de i'environnement ;

Vu la loi no 61-842 du 2 août 1961 relative a la lune contre les pollutions aûnosphériques et les odeurs :

Vu ia loi no 75-633 du 15juillG 1975 reiafive à l'éiiminauon des déchets et à la récupération des matériaux ;

Vu la loi no 76029 du 10 juiiiet 1976 dative à la protection de la nature, et notamment son axticle 2 ;

Vu la loi no 83430 du 12 juiUet 1983 moM& reiative a la démocratisation dw enquêtes pubiiques et a la

protection de I'errvironnemerxt ;

Vu la loi no 92-3 du 3 janvier 1992 sur l'eau, notamment ses amcles 10 et 11 ;

Vu la loi no 92-1443 du 31 décembre 1992 portant réforme du régime pétroiia ;

Vu la loi no 93-24 du 8 jauvier 1 9 3 sur la proteaion et la miSe en valeur des paysages et mt- CeRaines
dispositions légisiatives en matière d'enquêtes pûbiiques, notamment son article 2 1 ;
/

Vu le décret du ler falrier 1925 institUaatla commission interminisérielie des depots d'hydmahres ;

Vu le décret no 53-578 du 20 mai 1953 modifié portant rtgiement dadministraa'on publique pour I'appücation des
articles 5 et 7 de la loi du 19 décembre 1917 modinée reMm aux étaMissemems dang- insalubres ou
inwmmodw;

Vu le décret no 72-1240 du 29 décembre 1972 fixant les modalitcs de recowrement de la redevance annuelie
appiicable a certains étaMissements ciassés comme dan- insaiubres ou incommodes et le décret no 75-1370
du31 dkembre 1975 h a n t la liste des aaivit& soumisa à la m o n de ia redEvance annuelle appiicabie à
certains établissements ciassés wmme dangerew, iruaiubresou incommodts
__ ;

Vu le décret no 73-361 du 23 mars 1973 fxxant les modalités de recouMement de ia taxe Unique appkable aux
établissements classes comme dangereux insaiubres ou incommodes ;

Le Conseil dZtat (sectiondes travaux pubiics) entendq

Décrète :

page- 1 -
Gestion des sites (potentiellement/ pollues - Version O 11
Oi BRGM Editions - Dëcembre 1995
 Annexe ? - Aspects réglementaires

Artide 29

LesarretéSpris en application de l'alinéa précédem ahsi que ceux qui sontprévllsaux articles 10 (3è alinéa) et 11
de la loi du ISjuiUet 1976 sont pris SUT le rapport de i ' i o n des installations ciassées et aprk avis & consefi
d'hygiène. Ils font l'objet dcs mesures de pubiicitéprévues à l'article 27.

Le déchanta la Eiailti de se fave entendre par le conseil ou de désigner à cet &et zlll mandataire. n'doit être
informé au moins huit joms à l'avance de la date et du iieu de ia réunion du conseil et reçoit simuitanément un
exemplaire des propositions de l'iPspeaion des installagons dassées.

Le projet d'arrêté est porté par le préfet à la wnnaisanœ du dédarant, auquel un délai de quinzejours est accordé
pour préçcatcrévenaacllement ses obscrvatonspar écrit au préfes direaemaitapar-.
Artide 31

Toute modiücation apportée par le &clarant à l ' w t i q à son mode d'exploitation ou à son voisinage,
enmbant un changement notable des éIéments du dossier de dédamion initiale doit être prtée avant sa
réakaion à la c
0
-œ du préfg qui peut exiger une d e dédaration

Tout transf" d'une ïmallaion soumise à déclaration sur un autre emplacement nécessite une nouvelle
dedaration
Les déciarations privues a m aiinéas précédents sont soumises aux mêmes formalités que les déclarations
plrimitiveS.
- ._
__._. .. - .

Articie 32.

...
-Ladéclaration cesse depmiuire enet lorsque l'installation n'a pas été miseen service dans le cikhi de mis ans ou
iorsqae rwoitatioi a été interrompue pendant pius de deux années CO- saufle cas de force majeure.

-DISPOSïTIONS COMMUNES
A TOUTES INSTALLAïïONS CLASSEES

Artide 33

-
('DCcrctnO9U8.4àn 9juin ï994, ort;29 errtrcrcr orl.iguervsirmoisaprèsiapubIirntinn dù déCrel)
"Le direcîenr régional d e I'hdus&rie,-.delarrcherche et de i'en&ronnement est chargé, sous I'antorite du
préfet du département, de l'organisation de l'inspection des instpllatio
* nsdassées,

Les inspectaus des kîailations c i a s s h sont des a d & techniques désignés par le préfet sur proposition du
direaaV régionai de l'industrie, de la h e r c h e et de i'mvimnnement et relevant :

page - 17 -

Gestion des sites (potentiellement) pollues - Version O 27

-
BRGM Editions Décembre 1995
 Annexe 7 - Aspects réglementaires
a) De h diredion dgionaie de I'indristrie, de la h u c h e et de l'environnement ;
b) Des services VttCrinaks de la dirrction dtpartematt.de de I'agricuitnre et de la fo& ;
c) DeladirectIon dcpnrtancntliedes rdfputs sadaims etsocinles

Artide 34

(Déad no 94-484 & 9juin 1994,art 30)

" S a d dam le cas prévu à l'article -2, lorsqu'aae iastdhon ciasde change d'qloitant, le nouvel
exploitnnt ea fait la déciPration an pdfet dams le mois qmi sait la prise en charge de I'exploitdm Cette
-
déciantion mentionne, s'a s'agit d'une personne physique * les nom, prénoms et domicile da nouvel
espioitamt et, s'ii s'agit d'une pusonme morale, sa dénomination ou sa rnWn soci?k, sa forme juridique,
!'adresse de son siège social absi que b qualité du s@ataire de in déduntioa Il est délivré na récepissé
sans frais de cette dédaration"

kticie 34-1

(Décreî no 94-484du 9juin 1994,art 31)

-
"1 Lorsqu'une installatiori classée esi mise à I'arrêt d é f i i , son exploitant remeî son site dans un état tel
qu'a ne s'y mnnaeSte anam des dangers ou inconvénientsmentionnés à l'article ler de la loi da 19 juillet
1976 sumisée

Le préfet p t à tout moment imposer à l ' e x p l o i t les prescriptions relatives à la remise en état da site,par
&pris dans les formes prénies à l'article 1% c i - d n

-
ï ï L'ëxploitant qui met B l'arrêt définitason installoti'ounotifie aa préfet k date de Wadt an moins an
mois avant d e c i

ToutefoY, dans le cas des inst?llations natorisées poar une durée limitée Wiiu à I'artide 17-1, cette
aotifwtion est ad& au pdfet sis mois an moins avaut h dafe d'erpiration de I'antorisaâion,

ïU-Dans le cas des installati*oas Joamyes à autorisation, il est joint à la natakation un*dossKrcomprenant
l e p h à joar des terrains d'emprise de l'inseallnti-0% ainsi qn'mr.mimoirr s ~ l fl'état du site Le mémoire
précise les mesms prises ou p h n i e s pour assurer ta protectioa des intérêts visés à I'articIe lu de la loi du
19juilietWï6 susviskeîpowratcomporterwtunmentr

10 L'évparah'on ou l'&nination des pmduits dangerair,ainsi que des déchets préseatS sur le site ;
20 La dépouatiOn des sois et des eaax soutemaines iveatadlementpoihées ;
-
30 I'inserîion du site de I'instnllabon dans SOLI environnement :
4 O En cas de besoin, la suweiiiane B exercer de l'impacî de 1'hstaI.h'on sur son environnement

- -
page 18

28 Gestion des sites ipotentiellementl pollués - Version O

BRGM Editions - Décembre 1995
 Annexe I - Aspects réglementaires
Le préict consulte le maire de In commiine concernk En l'absence d'obsuvatons dans le délai d'un mois,
mn nvisest répritcflvorable

o&d no 94484 du 9 j h 1994, art 3 ' "Pour les hsîaUxm

2 * eristaipks faisant l'objet des dispositions de
l'article 16 de ka loi du 19 j d e t 1976, I'cxploitUrt doit f d r an préfet ks indications Saivamcs :"

l0 Sii s'agit diane personne physique, ses nom, prénoms et domicile ; s'il s'agit d h e personne morale, sa
dénomination ou sa raison sociale, sa forme juridique a l'adresse de son siège social ainsi que la quaiité du
si- de la déciaration ;
2" L'empkement de i'-on ;
3" La nature et le volume des activités exercées ainsi que La ou les rubriques de la nomendanire dans lesqueiics
l'insiallan'on doit être rangée.

Article 36 (abrogépm Surticle 33 dtc déCrci no 9#W du 9 j n b I994)

Artide 37

@écret no 94-484 du 9juU, 1994, e . 3 4 "Dans

) le cas prévu à l'uticle 35, le préfet peut exiger le production
des pikcs meation.aées anx artides 3 00.25 du préSent dkreî."
Le préfet pcur prescrire, dans les conditions prémus aux anicks 18 et 30 cidessus, les mesures propres à
sawegardtr les intéritSmentionnés a i'artide ler de la loi chi 19 jUiiiet 1976.

Ces mesures ne pewent entraîner ckmodifcations importantes touchant lé gros oeuvre de I v e o n ou des
changements consiciérablesdans son mode d'exploitation

Les Clispostions de i'aiinéa précédent cessent d'êû-e appiicables. si i'exploiîation a été interrompue pendant deux
années-oc saufle casde force majeure, ou sil'imtabtion se trowe dans les cas prévus aux articies 2 0 ,
31 ou 39 du présentdke!;

Arîicie 38

L'expioitant dîme installation SoIILZLise à autorisation ou a déclaratonest tenu à dédarer @&ztno 86-1289 chi 19
dhmbre 1986, a~ 5) "dans les meilleurs délais" à Ikspection des instaliationsciassks les accidents ou incidents
sunenus &fait du fonctionaementdecenektaiiation qui sont de naâire àporter atteinte aux intérês mentionnés
i l'article ler de la loi h19jUillet 1976.
Article 39

Le préfet peut &der que.,k e-tn-eœ- d'une-instaiJation momentanément hors d'usage par suite d'un
incendie, d'une explosion ou detait autre &dent résultant de !'exploiiation sera subordonnée, selon le cas,à une
noweLle autorisation ou a une nowelle déclaration

Artide 40

~ a arrêté
r pris après avis ctu conseil SupaieUr cies w o n s classées, le minisne chargé àes installations
classées peut procéder à i'agrémcnt de labolatoires ou d'organismes en vue de la réaisarion des anaiyses et
Contdles qui p m n t êîre prescrits en appkation & présent décret et mis à la charge des exploitants.
page - 19 -

Gestion des sites [potentiellement)pollués - Version O 29

BRGM Editions - Décembre 7 995
 Annexe 1 - Aspecrs réglementaires

@Lcrctno 94484 du 9jxin 1994, aat 35)

"DepnrrêtéS pris àam les mêmes formes firmt l
a conditions dans 1esqaeUeS doivent être p d s e n t é e leS
demnndes d ' e t "

A dtona pour l'expbhnt de prendre les dispogtions n k s a k s , il pouna ëtre f i t application &s procidures
privuc=S à I'articfe 23 de ia loi du 19j\Srra1976.
(Décret no 80-813 drr 15 octobre 1980. art 8.) "En cas d'appIication de l'article 15 de la loi du 19 juiïiet 1976 à
une installaton publique ou privée traMiIlaat pour les arwcS, le projet de &ret p&u audit article est soamis
pour avis auministre de la défense, avant sonewmen parle Conscii Sripinaudes indiations classées."

Article 42

Lorsquimeinstallati'on doit être Mplantée SUT le temtoire de plusians départernenu, la demaude on la déclaration
prévue au présent décm est ackssée aux préfas de ces dépmments, qui prodent à I'insuuaion dans les
conditions prévues au présent décret ;les dicisons sont pNts pararrêté amjoint de ces w et
s ,saufdans le cas
pr&m à l'article 16.

Artide 43

(Decret no 86-1289 du 19 décembre 1986, art. 6.) "Serapuni de la peine d'amende privUe pour les contraventions
de la sème classe:

"l0Quiconque aura exploité une installas'on soumise à dédaration sans avoir fait la déciaration prévue à l'article 3
de la loi du 19@et 1976 ;

"2O Quiconsue n'aura pas pris les mesures imposétS en vertu de Ilartide 26 de la loi du 19 juillet 1976 sans qu'ait
été pris, en raison de iingence, l'avis du maire ou de la commissioadépartenatntaleconsultativecompétente ;

"3" Quiamque aura exploité une kiaiiation soumise à auton9tion ~ a r r rs&f&~ aux prescript~ns à
I'article 7 de la loi du 19juillet 1976 et aux articles 17et 18 du prisent décret ;"

4" Quiamque aura exploité une insmllaa-onsoumise à déciamion sans satisbk auxprescnptions générales ou
partiaiiièresprRrues a m articles 28,29 et30 du prrSent dkm

5O Quiconque aura omis de procéder aux notificatons prévues aux articles 20 (leraiin&) et 3 1 ( l a aiin&) du
présmt&et

(Z)écrrtno 94-484 du 9j& 1994, mt 36)

" 6 O Quiconque aura omis de faire la désiaratioo ou la o&ication préwe aorartides-34434-1 du préSent
décret ;"

@émet no 94-484 du 9juin 1994, art 37

"70 Quiconque, après cessation d e I'apioitatioo: d ' u n e - i n s t 2 l l ? t i o n ~ ' - ~ ~ ~ n 'pas
w r arespeCté les
prescriptions de 1'an-W pdfecîorai pris en appücatioo de i'articie 34-1 du présent décret"

S0 Quiconque aura omis de fournir les informations prévues aux articles35 et 36 du Prcçent décret;

9O Quiconque aura omis d'adresser la déciamion prévue à l'article 38 ùuprSent déact ;

- -
page 20

30 Gestion des sites (porenrieiiernent)pollués - Version O

-
BRGM Editions Décembre 1995
 Annexe 7 -Aspects réglemenraires

@&et no 9 U 8 4 hr 9 j n h 1994, -38)

"100 Qlriconque aura mis en ocllyrt d u substances, d u prodpits, des organismes ou des procédés de
fabricrtioù s a m i s Q fgrbwnt en- du qoaîrhe aJinéa de I'utidel de la loi du 19 jnïiict 1976 psvUee
maus avollobîam I'pgréwnt oi1 saas avoirnspc&les c o n d i t i o l w p parcet
~ agdme"

du 9juin 1994, ort 39)

. ..(Décret no94-484
"mquwKmsdltMs niuopérritionssoimwesBI'igrémem
institriépar le quatrieme aibéa de I ' d d e 4 de la loi du
19 juilia 1976 modifwt s a s v k k

If-Unarrêtéduaiinistrr charge des établissements classés fue la CompOsition du dossier à fournir à l'appui
de la demande.

L'agrément est déiiVré par l'autorité compétente pour délivrer I'antorisation ou donner récépisst

11 ne pest être accordé que de&&= expresse.Le délai mPrimnm de déüvrancc del'agrément, qui c01ut à
partir de la date de notification de l'accusé de que l'arrtorité compétente adresse à l'erploitant
lorsqw le dossier est compieq est de trois mois. Ce délai pait itr+ p d o n g é parune decision motivée lorsqne
desconsultpti'om sont néassaireS.

Les délais préwxs ci-dessns cournt h partir de l'accusé de &cePtion que le préfa adresse à I'agloitmt
lorsque le dossier est complet.

IU - La demande d'agrément pour la mise en oernre d'organismes génétiquement modifiés dans me

in.+9iiltian ciasde est transmise à la commission de génie généticirie, notvnment pour déterminer le
ctvsemcnt des organismes mis en oe~lyreCet avis pent ne pas êîm demandé dans le cas'-& opérations
dev;rntlinmipisbe de la defense L'autorité compétente dispose d'nu-délai de huit jours à compter de la
date de l'aansé de réception du dossier complet pour transmettre la demande d'avis à ia cornmuSion d e
génie géuétiqne

La cOmmiSsjon dispose d'un déiai de quaranî&q jours poar formuler son avis Si elle ne s'est pas
prwonde dans ce déiai, son avis est réputé favorable.

-
IV La commSon de génie génétique est consnltée par le minisSe c h q é des installltion~clusées SU^ lu-
Rgles généraies applicables 201iastlllations cl?ssées figurant à la nibrique 26ûû d e ia nomen+$re, f&
en applicatiori des artides 7 et 10-1 de la loi du 19 juillet 1976. Elle dispose d'un déiai de 'dam mois +r
fomuier soo avis. Si eUe ne s'est pas prononcéle dans ce délai, son avis est r@té fawrabl&&

DISPOSITIONS TRANSITOIRES

Article 44

t a b ~ insalubres ou î l x o m résultant du dkret chi

A t i t s e ~ t o ~ l a n o s n e n c i a ~ d e s édangerrux,
20 mai 1953 m&é constitue la nomencisure des installations classées pour ia pioteCtion de i'enviromement
prévue à i'anicie 2 de ia loi du 19juiliet 1976.

page - 21 -

Gestion des sites (potenriellemenr)poilues - Version O 31

BRGM Editions - Décembre 1995
ANNEXE 4
Installations classées TEXTES
Circulaire du 24 mai 1976

conscil supt?rieur des installations darwrrs. LI doit renouveler sa demande d'autorisation

Ce d k t soumet les installations P autorisa- o u sa dalaration soit en cas de mnsferi. soit
tion ou a déclaration suivant la 3ravite d a en BJ dexunsion ou dc transformation de
dangers ou d a inconveninits que peut ses installations. ou de changenunt dans xs
p r Q m t a leur exploitation. p r o d i s de fabricatioioa. muainant d a
ArL3,- Sont s o u m P a a autorisation prifec- dangers ou inconvénients mentionnks a i:ar-
?O
.& les installations qui présentent de ticle lu-
gnivet - d a n g m o r r i n ~ u - p o r p o u rt c r ---
inririts visés a raniide lu. IL d 924S4 du 13/yüI. 1992. art. 7-11 La
miSe en CCUVTC. dans catrina af&.ories
L'autorisation ne peut are aocorda que si d'iristaliations classk. de substanas, de
ces dangm ou imonwininits peuvent ètre produits. d ' o r g a n b ou de pro&& de
prévenus par des mesuresque spéCifEl'arrëte fabrication peut, pour I'application de dira-
prefeuoral. tiva communautairrs datives à ia protcc-
tion de l'environnement. d u c subordonnée P .
IL d'93-3 du 4 j a n v . 1993. art.2-II) La un agrimcnt. Un décret en Conseil d'tut
düivrance de ïaurorisation. pour a s i-1- rue l a conditions Capplicruon du prisent
lationr peut éue subordonnée notamment a alinka et nouaunent la conditions de deli-
leur eloignement d a habitations. immeubks v r a e de I'apriment aiCrsi que l a d e k i dans
habitucllcment oacupks par d a tiers. etabtis- lesquels il est accorde. ou a ïexpiration
scments TeCCvant du pubiic COUK d'eau. desquels il eFt rbuze aaordé.
voies de communication. captages d'eau, ou
d a z o n a destin& à ïhabitation par d a IL. ri" 92446 du 13 j$l. 1992. ari.-6-lI) Un
documents d'urbanisme opposables aux d k r e t en Conseil d'Etat definit l a cas ct
tien. G Elle prend en compte l a capcitts conditions dans lesquels le changementd'a-
techniqua et financieres dont dispose le ploitant est soumis a une autorisation prtfec-
demandeur, a mème de lui permettre de totale déliVrte en considération d a capaatk
conduire son projet dans le respect da inte- techniqua et finan* aixoSaires pour
Ru vis& a I'anicle 1". >B metire en auMe I'activite ou remettre en état
le site diins le respect de la protection d a
Sont soumises à didaration l a installations interets vis& a ïanicle IQ.
qui. ne p k n t a n t .pas de tels dangers ou
inconvinients. doivent n h m o i n s respecter
les prescriptions générala c d i par le AR 4-1 (L. n"92446 ah 13 juiil. 1992.
prkfet en vue d'assurer dans le dkpanement la -
arr. 4-1IIj. Lcs dispositions prisa en appli-
protection d a interits VisCs a I'anicle Io. cation de la prisenu loi doivent. Iorsqu'etla
intërcsrnt l a décàets. prendre en compte la
Art. 3-1 (L. nu 91-13881 du 3 O d k 1991. objectifs de la loi no 75-633 du 15 juillet 1975
-
art. 2). Le srockage souterrain en coucha pdtée.
géologiques profonda de produits dange-
reux. de quelque ~ t l l t equ'ils soient. at Les dWions relatives aux instaüations d'tli-
soumis a autorisation adminkmtive. Cette mination des déchets prisa au titre de la
autorisation ne peut iue accordée ou prolon- prSsente loi doivent comporter les mesures
g&e que pour une dur& limitée et peut en pevues aux mida 7 et 7-1 de la loi
conskquence prévoir l a conditions de réni- no75633 du 15 juillet 1975 précittc.
sibdite du stocka- Inproduits doivent ètrc
retiris à ïexpintion de i'autorisation.
Art. 4-2 (L. n"93-3 du. 4jmv. 1993.
LOI No 76463 (Lnu 92-646 du 13 juiil.. 1992. CVI. 6-1) u A
I'issue d'une. piriode d.e CoaCtionnaneaL
art. 2-11, - La mix CD acuvitk. tant aprcs
SaÜiorisation initiak qu'aprù une autonsa-
DU 19 JUlLLET 1976 auto&.de vingtCinq ans au moins. l'autc- tion de changement d'exploitant des inctal-
relative aux inshuations Jarséts risation peut ètre prolongée pour une dur& lations définies par décret en Conseil d'Etat
pour la prottcoonde l'environnement iiiimitke, sur la base.d'un-bilan Ccologique praCnmt d a WUCS i m p o ~ de
n ~pollu-
comprenant une étude d'impactet i'utposi tion ou d'accident, des can5ks et da u-stai-
(JOdu 20 juillet 1976) d a soiutionr.altmiativa au maintien. d u lations de stockage de dkhcts est subordon-
, s t ~ ~ i u p ede tl c ~ Cn O ~ - U C O ~LZCCUOU-
S née a la constitution de garanties finanairrs.
veiIemmt s'accompagne d'une;. nouvelle
TITRE PREMIER .évaluationd a garanties f i è r e s prorUa à ~a -garanties' sont desinéCs à assum.
olsPoslrlolvs GÉN€RALES-- I'artide 7-1 d e la loi no 75633 du ISjUiiiet suivant la n a t u r e d a dangers ou inconvk-
.'1975 relative a I'éiiminationdes déchaset à .nienrsde chaque catégo&âïmallationr la
la. rinrpiration'do-matiriauxou à ïarti- survallancedu site et Le maintien en Sinirite
Articie p r e m i e r ( l d 93-3 &.4jmr,1993. cleS3 de ia loi .no 87-565 du 22 juiilier t987 de 17ixiuliation.l a inrmrntioat tvrntueUa
~WK- -
I**-I) -Sont soumis a u x dispositions d e
la priscntrlohIct- rrsina;:atdiers=,dépàtr
'relative à I'organisatioa de ia sccUrit~'avile;à - encasd'accident avanrouapràlafemieture.
e t ia remiseen Mt apris f a m u u s Elles ne
chaukctd'uncrmniirrginCraklesinstal- a indamisaiions dues par
couvrent p ~ r t l
lations exploitées ou: dircnuec-par toute ïcxploitant aux tien qui pourraient subir un
pasonne physique ou morale. publique ou Pour l a S t O C k a g e s ~soutaTairrs .de&ets pn5judie par fu't de pollution ou dàc9dent
pri* qui peuvent pr&scnter des dangers ou t ultimes, L'Agena d e ltenvironnanentet de la
d a inconvCPicnts soit pouda cornmodili du mU& de l'énergie peutconclureavsra-.
voisinagq soit pour.k arikia--té. k- :ploitant. avant ïoctroi deltautorïsatioa visie
salubrité pubüqria. soit pour ïaagricUlturc- :au premier a ü n s une convention quïdéter-- p a md a * g a m t i t set les.@CS de firion
soit pour la protection d e la naturcet d e ;mine ks conditions techniqua et financi& de4euc montani:
ïenvironnemmt soit pour k conservation de I ' n i g a g ~et~de ~ la~ poursuite
~t del'utploi-; -*-%-&Y.. _- . ,=..
,.an- . .
des sites et d a monumentL -- - 2-
:tatios q>mpte,tcnu d e ïivcntualité-du refus ,Sanspréjudice-de. k.procid~&d'amende
l
-c .-.- 'de.-= 'prolongation., Cette. convaition-est , administrative ptinie à ïanicic-7-1 de la loi
(Lna 93-3 r6r 4 jcmR 1993. art. l'YI) k soumise pour avis a u représentant de FÉUL n? 75-633 du-ISjuiilu 1975 relative à ïkürni-
dispositions de la prtsCnte loi sont égalanent __ .. . .. . .
a p p l i a b k aux exploitations de ~ r r i t r au
sens d a &da lu et 4 du Code minier.
a 'L;adispositions da de& alin& prtçadnik
nes'appliquent pas au stockagesouterrain de
nation.d& dëchetiet a h rkup6nlion d a
matériaux pour l a i d l a t i o n s de stockage
dédieu radioactifs. D .. - de déchets. l a manquements aux obligations
-
Art. L Les installations Vism a Iànicle 1- de garantis TcnanMcs donnent lieu- a Sap-
sont dtCinies dans k aomtnclature des insral- -
Ah 4. L'exploirant est tenu d ' a d k r sa plication de la procédure & consignation
lations c k î s k s établie par h i en Conseil demande d'autorisation ou sa didaration en privUe a I'anicle 23 de la prixnte loi. inde-
d'État, pris sur le rapport du m m s ~ r charge
e mMe temps que sa demande de permis de pendammeni d a poursuites pknaia qui
d a Installations c1ass.k. apres avis du construire. peuvent itrc exercks.

4128 F e u i l l e t s 88 (15février 1993) C o d e PermaGent Environnement et Nuisances

Pour une information plus r k e n t e reooflez-vous a la Table Rkapitulative du bulleltn
5
Installations classées TEXTES
Loi no 76-663 d u 19 juillet ;976

tuer une panie du prix ;il peut aussi deman- s'imposent de picin droit aux installations autorisations d'exploitation de carri& pour
der la remise en etat du sicc aux frais du nouvella. Ils pririsent. ap& avis des o r g - ksqudla:iedéiaiderecouautf~àsllmois
vendeur. lorsque le coiit.de a t t e remise en nisations profasionnelles interessées, l u a compter de l'achévancnt d a focmalit& de
itat ne pzrait p ~ sdisproportionni par dckiset la conditions dans Iesqwk ils s ' a p publiate de la déclaration de déùat d'exploi-
n p p o n au prix de vente. pliquent aux installrtionsaistantes. tation trapsmise par ïaploirant au m e n -
Arf 9 (Lno 90-557 du 2juiü. 1990. art. 6). - It précisent tgrianent les conditions dans tant de I'Etat dans le dtpurement. s
Dans l a communes compomnt une aire de ksqucias ccs pracriptions peuvent être iatiasquin'ontacquisoupriràbzild~
-produaion -de vins - d - a p p d î a t i o n - t f o t i ~ - . a d a p t k ~ ~ a t o n k a n x - c-immeubles i ~ ou-n'onc-tlc*i-derwon-
l'avis du ministre de I'Agridture doit tut UMaloCrhs. dans le V O ~ S ~ M d'~w iastaiiah daztte
demandé cï1 vue de i.autor9Ition prévue au Art 11 (L n.93-3 du 4janv. 1993. q w pontiieurrment i i'atliigc ou a la
pmnier alinéa de I'artide4 a-âcsus Cet
avis ut donni apr& coasuiratio~le gs
art. 34I). - Si 1s inténits mentiorin& à publication de l'arrètë autorisant ïouvaturc
de cette insraüation ou atténuant ks p r r s r i p
r&eiu & la p r t u n u loi ne sont pas
&c!kant. de I ï i t u t nationai d u appella- -ris par r ~ t i o d n u prescriptions tions primitives ne sont pas rsPcMblez a
tions d'origine. ghénks COQLÏCks inconvsphnu inhércnuà dXércr M i t a m t é à la jundiaion admiais-
Le ministre de I'Agriculture ut en outre rupioitation d'une instaüation soumise à tntiw.
consulté. sur sa demande. lorsqu'un etablis- dickrrtion. le préfet. tventueIlancnt a la (Lf f 76-1285 du 31 a&. 1976. art. 69) n Le
sement soumis à ïautocisation Wie adessus danande da tiers i n W et ap& avis a de permis de construire et ïactc de venu, a da
doit être OUVM dans une commune liitro- h cornminion départanentale consultative tiers. de biens fonciers et immobiliers
phc d'une commune Comportant une aire de compilcntcs. peut imposer. par arrêté. doivent. le cas échéant. mentionner cxplicite-
production de vins d'appeilation d'origine. mutes &utions suécialu nécessaires. mcnt lu servitudu affërmtes innitu& en
n Le ministre de I'Am-iculture et de la Forët

qu'un etablissement soumis à l'autorisation p&cr. k rcprëscncant de l'Êta[ peut pres- Art. 15 (L n"92654 du 13 hiil. 1992.
visée cidessus doit être ouvert dans une crire la rkkation d u évaluations et la mise art. 7 - I I I ) . -S'il appam-t qu'une installation
commune ou dans une commune limitrophe en œuvre d a rrmtda que rendent n k s s a i r u classa priscna. pour les intérëtsmentionnts
d'une commune comportant une aire de soit les conséqucnœs d'un aaident ou inci- à l'article 1'' cidessur des dangers ou des
production d'un produit d'appellation d'ori- dent survenu dam l'installation. soit lu inconvknienrs qui n'étaient pas connu lors
@nec o n t d t e autre que le via * conskquenrrz enminées par l'inobservation de son autorisation ou de sa dédaration, le
des conditions i m w en application de la ministrechar@ d a installationsckrva peut
Le ministre de I'Agriculture dispose d'un prismte loi. Ces mesura sont prucrita par ordonner la suspension de son exploitation
délai de trois mois pour donner son avis. Ce des &tés pris, saufcas d'urgence. ap& avis pendant le delai nécessaire la miSe en œuvre
delai court a partir de la date à laqueüe il a été N de la commission dipartemmtale consulta- des mesures propres a faire disparaitre ces
saisi par le préfet du dossier auquel at joins tive compétente ». dangers ou inconvenimts. Sauf cas durgen-
son avis.
-
Art 12. La installations qui, soumisu à cc la suspension intervient aprks avis des
organes consultatifs cornpicents et apks que
diciaration en vertu de la prisente loi. biné-
TITRE Ill f i t d'une autorisation riguliere avant la i'exploitanr a etc mis a mtme de prtsenter ses
date d'en& en vigueur de la loi du 19 de- observations.
DlSPOSITIONS APPUCABLES
AUX INSTALLATIONS m b r e 1917. sont dispaks de toute décla- (+II" 92#6& 13j+iI. 1992.mt.4-VIl) Un
SOUMISES A O&URATlON
ration :dla sont soumises aux dispositions demet en Conseil d'Etat pris aprS avis du
d u arrjcla 10 et LI. conscil supérieur da installations daaia.
peut ordonner n la suppnsion ou la ferme-
Ah 10 (L. n'9û-85 du Ujanv. 1990. TITRE IV ture * de toute installation,figurant ou non a
art. 11-11 et L no 93-3 du 4jamv.. 1993. la nomenclature. qui présente. pour l a inu-
art.3-il). -
La presCriptio2 gênédes DISPOSlTIONS APPUCAôLES
A TOUTES LES INSTALLATIONS rèts mentionn6sà I'anide 1 . da dangers ou
prévues à t'anicie 3. dernier aiinéa. sont ëdic- inconvénienrs tels que l u mesures prévues
tm par arri<a préïectonrq p N a p r à am CLAss€ES
par la prauitc loi ne puisent l a faire dizpa-
n de la connumon diparimimtale consulta-
uve COaip6tcIIUB a.pour k atciius hors Art; 13. -
La pasonus c h a r e de I'inspec- raitrc. ..
htt 16 (L n' 93-3 du Q@v. 1993.
sol. de la co&on dipartwcntale d a tion des iriarllarionsda;,ou d'apeniws
structuits agricolu. Eiics s'appliquent aucw sont-assamenck et astreintes au secret -
art. 7). La installauons qui. a p r à avoir été
matiqucment à route installation nouvclle ou profasionad dans la conditions et sous l a réguiiércment mises en servi- sont soumi-
soumise à nouveiie déclaration. sanctions privws à ranide 3 7 s . d ~GAC ses,e n w n u d'un décret relatif à la Nomen-
Les modifications ulcérieurrsde ces p r u a i g ptiieL.ktucllement. aux arriciu . .
70 et daturc d u installations classés, i autoksa-
suxvants du même codr .. tion ou à dédaration peuvent continuer à
tions giniaia peuvmtétrc renduesappiica-
bles. au* insullations erristantn don l u ~ p c w c n t ~ t c r à t o u t m o m a t i a i i s ~ - fonctionner sanscetcc autorisation ou dida-
latiom-soumises& h r surveillance;=? .i . ration a la seule condition que rcrploitant se
modalités et d o n ks dilais prénn dans soit déjà fait connaive du r e p b t a n t de
i'arrett préfeuorai qui f i igakment l a A& 14:. (c-992-654 du 13 '~%iC".1992, I'État; dans le. dépanement OU. sc fasse
conditions dans laqucues I c s plucriptioar- 'are.-7 - F Z L"'n' 93L3 . du 4 j a n n ' 1993, connaitredelui dans I'aonœ suivant la publi-
@&aies peuvent être a d a p c k aux aiu~us-
:ta-locala- .
- - .wr-9)5 Lerdéüsionsprises eu application
.derarrida3;a4;4-2», 6.11; 1 2 16;U;24 e e
cation-du d k e t .y

---.II-

26 de- k: prtrcnte loi- sonvsoumises;à un'' fa 'rr;;cigiiemeàu q u e i~uploiLant~

doit
L& &blissanents sou& a décta;alion&us VMmetuc au reprisn+icic ritardansle
je k g h c de la loi d u 1 9 d é a m b r e 1917 ec contentieux -de pleine juridictioa ..EIIes.
dipancment ainsi que Iesmaurrs que œlui-
ayant, ob-4 m vcltll de- l'artick19;. peuvent int dX* a la juridiction admi-.
a peur imposa afin de sauvegarder l
u inté-
alinëa 1" ou 4 de ladite loi;la supprasiou réu mmtionnks I l'article 1" sont précisés
ou i'attinuation d'une O: piusicurs prescrip- Io Par kr dermedeun ou exploitanu, dans par décret m Coowil &iht
tiok-"raultanL d ' y m r s pféfeaonius. .uudéiai.dednirmoisquicommcnocicourir
'cowrvent le W i œ de CCT dkogations- K .dupuc.où k u i i t s , kur ~ été notifiés :
~ ~ ont
peuttoutcfoisy ttrcmnfin paearrëttptcfœ- T[TREN BIS
:2? P a r Z l & e pm0kcs:physiqUa- ou. - .
torai pris aprks a v i s - a d e k cornmimon morrlcs; kcommunes intksskes ou leurs
'dtpanementale consui~compttcnte,.. ;grciupanc+k M o ndesinconvtnieatsou-
:. - DlSPOSITIONS APPUCABLES .
--
-. AUX EXPLOlTATiONS DE CARRi&ES
selon l .
amodaiik*.---..,
et dans k-
.- -_..,
düai
.
fiés par '.desdan-
-.,.----. .-.*-
- c . -
q u e k fonaionnaacnt. de I ' i (Lno 93-3 du 4 janv. 1993, art. 8 )
. .. -A.,-. .
~

- c a l h i ~ p r 9 t n e p O u k-
P intiritr.*à. ._
&. lû-1 --(i f f 92454 . du' I3-'@ ïaraclc l";,danf.Un dtti-dequatre-aris à
-
1992.sart. 7-Ir). Pour la protcuion des
intérëtsmentionnis à I'artide leadessus. le
.eoasp~~~deL:pubücation
desdiu -ad&
oude-i'alrichage
érant. le cas éckéant
An 16-1
art-8).
(Ln'93-3 du 4janw. 1993.
- La exploitations de carriira sont
ministre chargé d a inscaiiacionsdassm peut proloo+jusqu'i k ira d'une përiodede deux soumises à l'autorisation adminimative
fua par arrêté. après consuication d u mi& années suivant k misc en activité de ïinual- prinre à I'anKle 3. ,
tres inté& et du Conseil suprieur d a latioa .--- L ' a u t o h t i o n administrative vistc à l'alinéa
installationsciassiia. les prrscnptions +lé- (L. n'93-3 du 4jmv. 1993. ari.6) u La p M e n t ne peut u a i d c r trente ans.
ales applicables a ccrca~ncseatigonu d'ins- dispositionsde I'alinia précédent ne sont pas Cette auronsatton ne peut exctder quinze ans
tallations soumtses a d a h r a t i o n Ca arrttes applicabh aux dkisions concemant l a pour les terrains dont le défrichement cst

4130 F e u i l l e t s 88 (15 février 1993) Code P e r m a n e n t E n v i r o n n e m e n t e t Nuisances

Pour une informabon plus récente reportez-vous a la Table R k a p i h i l a t i v e du bulletin
7
 Insta1lations classées TEXTES
Loi no 76663 du 19 juillet 1976

puni d'une peine d'emprisonnement de ... (1) La p e i n a encourua par l a personna mora- l'exploitation de l'installation jusqu'au d e t
a....(l)etd'uneamendede...(1)à...(l) F.ou l a sont : de la dalaration o u jusqu'a la d e i o n d a -
de l'une de c4 deux peiner 1' L'amendc suivant les modaiitks prevua tive i la demonde d'autorisauon.
II. Quiconque pounuit l'exploitation d'une par l'article 131-38 d u Code +rial ; Si ïupioimnt ne défère pas a la mise en
installation clauic sans Y confomer a duneurc de régulariser sa situation ou si sa
ï&ti de mise en demeure d'avoir a respa- 2- La peina mmtionnk aux 2'. 3'. 4'. Y. demande d ' a u t o h u o n es rrjetcS le préfet
ta. au urrne d'un d i h i firt. les prescriptions 6'. 8' et 9' de ïanicle 13 1-39 du mame c o de peu^ en cas de n k s s i l t . ordonaa la ferme-
techniqua détamin& en appïkation des L ' i n d i o n mentionnœ au 2' de l'ar- ture ou la suppression d e ïinstallatioa. Si
-anidet3-,-6.7. thel 1saapuni-d'unepcine-- - tick 13 1-39d u mime code ponc sur I'activitt -Ibcpioitunt n'a pas obtanpiré-clrns-kdani
d h n p r i s o ~ ~de n t._. ~ ~a ...(1) et d'une
~(1) dans ïexcrcice ou a ïoccasion de ïuercice de fixe. le préfet peut faire appüation d a
amendede ...( I ) a ...( 1)F.ouderunedeacs laqucile l'infraction a été commise (2). procidures prtvua a l'anicle 23 (3= et 4=aii-
deux peina. ncIis).
(L. t f 93-3 & 4 + ~ . 1993. ~ f 3-11) . Sem
puni d a mânes peines quiconque poursuit
niw VII Le prifet peut faire procider. par un agent de
la foroc publique. a l'apposition des a l l a
ïuploitation d'une instaiiation sans se SANCTIONS ADMINISTRATIVES sur une installation qui est mainunue en
confomer a un arrête de mise en demeure fonctionnement soit en infraction a une
pris en a p p ï i t i o n de l'anide 26 par le repri- ArL 23 (L n*92-654 du 13juiIi. 1992. mesure de suppression, de fermeture ou de
sentant de l'Eut dans k département sur avis
du maire et Y de la commission départanen-
-
arf. 7-Y). Indépendamment d a pounuita suspension prise en application de ïarci-
de 15. de l'article23 ou d a deux premiers
p é M k qui peWCXlt C X U C k et lOlSqU'Un
tak consuitativc compkcentc ». inspecteur d a installations cl& ou un alinéas du p-t article. b i t en dépit d'un
(L.tf92-646 du 13 j d l . 1992. arf. 12-11) cxpcn désigné par le ministre chargi d a amite de refus dautori9tion.
III. Quiconque ne Y conforme pas à i'arrite
de mise en demeure de prend% dans un düai
Installations cias& a constate I'inobsena-
tion des conditions impos&s a l'exploitant
-
AR 2% Pendant la du& de la suspension
détcnnint. l a m a u m de surVgliancc ou de de fonctionnement prono& en a p p ï i t i o n
d'une installation claaie. le prifet met en de l'article 23 ou de l'article 24a-daais. ïex-
tunise en etat &une installation ou d e son demeure cc dernier d e satisfaire a ces sondi-
sitc prcsaita en a p p ï i t i o n des anicks 6.7. ploitant a t tenu d'assurer a son personnel le
tions dans un délai détermine. paiement d a salaires. indemnités et rémune-
10. 1 1. 15.24 ou 26 lorsque I'activité a Cast
sera puni d'une peine d'emprisonnement Si. a l'expiration du délai fut pour Iéxtcu- ntionsde toute natunauxqueis il avait droit
___
de (1) a ._. (1) et d'une amende de ... (1) tion. l'exploitant n'a pas obumpiri a cette jusqu'alon.
a _.. (1) F. ou de l'une de ccs deux peina. injonction. le préfet peut :
TITRE Vlll
AR 2 i IL. no 85-661 du 3JkiiL 1985. art. 4 ) . (L. no 92-654 oh 13 juül. 1992. an. 7-Y)
u a ) Obliger l'exploitant a consigw entre les DISPOSITIONS DIVERSES
Quiconque obstacle a I'e~ercied a fonc-
tions d a personna chargks de I'ins&on mains d'un comptable public une somme Ah 26 IL. n"9J-3 du 4 janv. 1993.
ou de l'expertise d a insrahions ciassies
sen puni d'une peine d'emprisonnement '
ripondant du montant d a travaux a rkliser.
laquelle sen restituée a i'urpioitant au fur et a
urf. 3-11). - Lorsque l'exploitation d'une
installation non comprise dans la nomencia-
de ... (1) a .__ (1) et d'une amende de ... (1) maure de l'exécution d a mesures prcsCntes : ture d a installations classies prisenle des
a _.. (1) F, ou de l'une de m deux peina. il a t proadt au r u o u m m m t de m e somme
dangers ou des inconvcnienrs graves pour l a
comme mari& de créances etrrngks i
Arf 22. - La infractions sonc constatirs par intertu mentionnks a ranick 1" de la
les procis-Verbaux d a offiers de poiice
I'imp6t et au do-inc. Pour le recouwemmt
-
présente loi. le prefet. a p e avis sauf cas
judiciaire a d a inspecteurs d a installations
de une somme. l'Etat k c f i e &un priviltge
de même r a q que celui priw a ïanicle 1920 -
d'urgence du maire et cc de la commission
ciassk. Ca pro&-verbaux sont dresSa en du Code générai d a impôts : départementale consultative compétente»,
double exemplaire dont i'un a t ad- au met l'expioitant en demeure de prendre l a
préfet et l'autre au procureur de la Républi- b) Faire proctder d'office. aux frais de l'ex- mesures ntccStaim pour faire d i s p a r a b ks
que. Ils font foi jusqu,'a preuve contraire. ploitant, a l'exhtion d a mesures prcvrita : dangers ou l a inconvénients dümcnt consta-
c ) (i.. n'93-3 du 4janv. 1993. art.3-II) t&. Faute par l'exploitant de Y conformer a
AR 22-î (ILno 92-1336 du 16 dëc. I992.
-
u f . 307). & tribunai peut ordonner l'afli- Suspendre par arrite. apris avis u d e la
commission dépariementale consultative
cette injonction dans le délai imparti, ii p u t
itrc fait application d a maufcs prévues a
chap ou la diffusion intégraie ou partielk de
la daiPon prononcée dans les conditions compétente », le fonctionnement de I'instai- l'article -3 ci-dessus,
prtwcS par l'anicle 131-35 du Code pkd. latioo jusqu'à e x h t i o n d a condkions A h 26-1 (L. no 92446 du IJ J X L 1992.
An 22-2 (L n'8%61 du 3Jkiii. 198s.
imposstS et prendre l
res nksaircs.
a dispositions provisoi- -
m. 4- Vlll). Lorsque l a personnes moiala
mi.6). - Toute asochiion Rguiièrancnt
Les somma consignks en application d a
de droit public interviennent. matmellemmt
ou financitmncnt. pour attinuer ladomma-
déclarce depuis au moins cinq ans à la date
des faiu Y proposant par ses statuts la dispositions d u u peuvent ètre utilistes pour g a rhultant d'un incident ou d'un accident
sauvegarde de tout ou panie des interèts vis& régier l a dépenses enuainées par l'exécution cauY par une installation mentionnée a i'ar-
à ïuricle premier de la p k t e loi. p u t d'ofice d a m a u ~ c ps R v u a aux 6 et c. u ticle2 ou pour évita l'aggravation de cb.
umrr l a droits reconnusà la pame civile en dommages. ella ont droit au tunboursc-
(L. n"92-654 du 13 juili. 1992. urt. 7-VI) ment,par l a pusonna rrfponsabks d e t'in-
--qui concerne l a faiu constituant une Lorsque l'état exécutoire p N en application
infraction aux dispositions de la p k c n t e loi adent ou de Saccideni d a frais qu'eiies ont
d'une meave de consignation ordonnœ par engag&. sans prijudice de t'indemnisation-
ou d a r&glcmentset arrêt& pris pour son l'autorité administrative fait l'objet d'une
application et ponant un pRjudicc direct ou d a autres ,donunaga subis. A oc tiuc. d a
opposition devant le juge administratif. le peuvent se constitua partie avile devant l e '
indirect aux inttrirs c o u d s qu'elle a pour président du tribunal adminisrratif ou le
objet de défendre. juridictions pénales saisies. der pou&-.
magistrat qu'il délige. statuant en référé. c o n k u t i v a a l'incident ou a l'accident .. .
Art 2?4 (L n'85-661 du 3 W l - 1985. peut. nonobstant cette opposi!ioe a la Cette action s'excrœ sans préjudia des droits
-
art- 7)- Pendant la du& d e t'interdiction demande du reprkntant de l'Etat ou de
toute pcrsonne int- d e d e r que le ou- par l'article 22-2 de la p h u loi
d ' u t i i l'insdlatioa pronon& en applica- aux associations ripondant aux conditions
non de l'article 19 adessus ïexploitant est m u r s ne seza pas suspensif. d& Ion que l a de cet article. . . ... .
tenu d'assurer a son pcrro~elle paiement moyens avan& par ïcxploitant ne lui
des salairrs.indemnités et rémunérations de paraisoCntpas - . Le prisident du tribu- AR P (L. n"87-565 d u - 2 2 . S 1987.
toute nature auxqueis celui-ci avait droit nal statue dans ks quinze jours de sa saisine. -
art. 25). En u qui wnCQnel a instaü+ons
jusqu'aion. -
Ari 24. Lorsqu'une installation clas& est
appartenant a- servias et organrmier
dêpcndant d e l'Etat qui Seront inscrits sur
Ari Z-4 (Ln'92-1336 du 16 dic 1992. exploit& sans avoir fait l'objet de la déclara-
-
art. 308). La personnes morales peuvent
&c dkiarées responsables dans l a condi-
tion ou de l'autorisation requise par la
prrSmte loi. le préfet met iicxpioitant en
une liste etablie par diCrci. les pouvoirs anri-
buis au préfet par la priWnte loi seront
e x e soit par le ministre c h a r e des Éta-
tions prémes par l'article121-2 d u Code demeure de rtguhriscr sa situation dans un
düai dé&& en déposmt suivant le cas biiients c b soit par le ministre
p k d d a infnaions définia aux artides 18 chargi de h Wf- pour l a installations
et 20 de la p k t c loi. une dédaration ou une demande d'auto&-
tion. 11 peuf par arrêté motive. suspendre qui r e l k n t de son d ~ n e m m t
u La dispositions d a articles 7-1 à 7 4 de la
(1) NDLR : Pour le aux de ces peines. v. au p r k n t e loi ne sont pas applicables à d e s de
comnunwùe: Tabieau &s sanc~ionrpimles. (2) .&tir& en vipeur uù nouveau Cide pinal Ir ces installations qui relévent du minisVe de la
I f W. I" mars 1994 ( Ln' 93-913 du 19juill. 19931. Défenw. m

4132 Feuillets 93 (15 février 1994) C o d e F' e r m a n e n t Environnement et N u i s a n c e s

Pour une inforrnation plus récente reponez-vous a la table d'actuaiis.ation. e n tete du Code
9
TEXTES Insta!laticns classées
Circrrlairs du 24 aoüt i976

La +nalit& p & v u a a u titre VI sont appli-.

a b l e s aux justiciables des juridictions mifi-
taires d a foras a m & confonniment au
Code de justia militaire et notamment cn ses
anicla L 56 et 100.
-
Ait ZS- La modaliris d'appiication de la
présente l+ seront fw-pair dicms m
ConaiI+E6f1r
.
Ces d k u d i m i n e r o n t QI outre:
1" Pour ks insullations vis& 6 I'anide27
a-dersrn. l a p r o c t d u u &enquête et d'auto-
risation. ainsi que les conditions & s u d -
lana a de conudc ;
P Pour les au- services d e I'Etat. ainsi que
pour les coiicctivittr locales et l a établisse-
nunu publics a caraukre adsninistratif:
a ) La conditions d'application d a mesures
prbues aux arcides 19.23.34.25 et 26 ;
6)La pusoanes qui seront rrgardascomme
pinaiem«it rcsponsabies des infrauions
cornmisa.
-
A r t 29. Les dispositions de la prbente loi
entrent en vigueur le 1"janGa 19n. A cette
date, sont a b r o e la loi rnodifik du 19 de-
cembre 1917 relative aux ctablissmients
dangereux. insalubres o u incommoda, le
d k t - l o i valide du 1" avril 1939 insraurant
une procidure d'urgena pour I'instruction
des demanda de construction de dépots
d'hydrocarbures. et les dispositions applica-
bles aux insrallations soumises la p h n i e
loi et qui lui sont contraires.
La référence a la p r k n u loi est substitute a
la reFeience 1 la loi du 19 décembre 1917 dans
cous les textescontenant une telle disposition.
présente loi sera exicuttc comme loi de
I'Etat.

Code P e r m a n e n t Environnement et N u i s a n c e s Feuillets 88 (15 f é v r i e r 1993) 4 133

Pour une inforrna!ion plus recente reporter-vous a la Table R h o i t u i a u v p d u bulletin
- 10
ANNEXE 5
 Nom :MARTHE
Société :B.R.G.M
Type de modèle
Modélisation hydrodynamique et hydrodispersive des écoulements souterrains
en milieu poreux.
Modélisation des écoulements en 2D ou en 3D anisotropes horizontal-vertical
Aquifères monocouches ou multicouches
Ecoulement en régime permanent ou en régime transitoire
Ecoulement en zone saturée ou non saturée
Migration dispersive d'un effluent dans la nappe
"Transport inverse" pour la localisation des sources de pollution
Couplage du transport avec des interactions physico-chimiques entre eau,
effluent et matrice poreuse
Ecoulements diphasiques, eau-gaz, eau-huile

Objet et champ d'application

Résolution des problèmes d'écoulement dans plusieurs domaines
Gestion des aquifères
Génie civil : assèchement des fouilles - simulation de l'effet hydraulique des
parois étanches
Environnement : migration des polluants dans une nappe
Stockage souterrain de déchets : étude du confinement d'un stockage souterrain
Exploitation minière : calcul de débits d'exbaures miniers et des rabattements
associés

Méthodes et techniques utilisées

Résolution des équations d'écoulement par une méthode des différences finies
Calcul du transport par différentes méthodes selon l'importance du phénomène
de convection ou de diffusion (différences finies - méthode des déplacements aléatoires -
méthodes des caractéristiques (méthodes MOC)).

Type de transport
Transport de soluté
Transport polyphasique (eau-gaz, eau-huile)

Données requises
Paramètres hydrauliques du sol
Perméabilité : K = f(h) (paramètres décrivant la fonction )
Coefficient d'emmagasinement spécifique (courbe caractéristique de l'humidité
du sol)
Teneur en eau de saturation (porosité)
Paramètres de représentation du milieu
Données spatiales (cotes du substratum, topographie ...)
Paramètres de discrétisation et de maillage
Paramètres de transport
Coefficient de diffusion
Fonction retard (paramètres du phénomène d'adsorption-désorption)
Terme source
 Conditions initiales et aux limites (de l’écoulement et de transport)

Expériences acquises et validation

Modélisation des aquiferes
Ecoulement à travers la berge d’un cours d’eau
Réalimentation artificielle de nappe
Ecoulement dans un aquifere en contact avec la mer
Ecoulement en régime transitoire produit par convection naturelle

Limitations
Modèles des paramètres hydrodynamiques ( expression unique pour la fonction
perméabilité, courbe de rétention)
Validité de la cinétique d’adsorption-désorption
Difficilement utilisable dans les sols fins à faible perméabilité à cause du
phénomène de partition de l’eau interstitielle.
 Nom :TAFFETAS
Société :BRGM
Type de modèle
Simulation des écoulements et transport de masse (multi-composant) et de
chaleur dans un milieu poreux et/ou fracturé (régime permanent ou transitoire)
Représentation du milieu fracturé par le mélange des éléments finis à une
dimension (écoulement dans les canaux ou dans les cavités), les éléments à deux dimensions
Etude des écoulement en régime permanent et transitoire dans un domaine
bidimensionnel ou tridimensionnel
Etude du couplage hydrodynamique, du transport de masse et du transport de
chaleur grâce à la variation de la densité et de viscosité du liquide
Calage automatique des paramètres lors de la simulation par la méthode de
Khan-Newton

Objet et champs d'application

Calcul de l'écoulement d'un fluide incompressible monophase dans la matrice
poreuse (problème de filtration) et dans les fractures planes par utilisation de la loi de Darcy
linéaire dans le milieu poreux et de la loi linéaire ou non linéaire dans les fractures planes
Simulation du transport des polluants radioactifs et non radioactifs avec prise
en compte des réactions d'adsorption grâce à l'utilisation du facteur de retard
Application au calcul des transports dans les couches rocheuses fracturées
tridimensionnelles
Gestion des aquiferes
évaluation des termes du bilan hydraulique (entrées-sorties, stockage-
déstockage)
interaction entre forage
Environnement
modélisation de la migration d'un ou plusieurs polluants dans une nappe
Stockage souterrain des déchets industriels
étude du confinement d'un stockage souterrain
Géothermie
étude des doublets géothermiques

Méthodes et techniques utilisées

Modèle utilisant le code de ROCKFLOW développé par l'université de Hanovre et
mis en forme par le BRGM de 1989 à 1994.
Ecoulement
Utilisation des éléments à 1,2 et 3 dimensions pour la description de
l'écoulement de l'eau ou de la solution dans les fractures, la matrice poteuse et la matrice
poreuse connectées à des fractures
Ecoulement décrit par la loi de Darcy généralisée; utilisation des lois
d'écoulement non linéaires dans les fractures
Prise en compte des hétérogénéités et anisotropies dans toutes les directions
Transport de soluté
Description de transport de mono et poly-composants chimiques et/ou de
chaleur selon trois méthodes en fonction du problème posé :
Méthode de Taylor-Galerkine avec correction des oscillations liées au
phénomène de convection
 Méthode des déplacements aléatoires (Random Walk), appelée méthode
des particules
Méthode des caractéristiques (MOC), combinaison des deux méthodes
précédentes
Couplage de l’hydrodynamique, du transport de masse et du transfert de chaleur
Calage automatique des paramètres du modèle selon une variante de la
méthode de Khan-Newton

Données requises
Matrice des conductivités hydrauliques et des coefficients d’emmagasinement
Facteur de répartition des fractures, porosité
Coefficient de diffusion, coefficient de dispersion longitudinale et transversale
Fonction densité de la solution et fonction viscosité dynamique (nécessité d’un
modèle de fonction)
Conditions initiales - conditions aux limites sur :
Flux de liquide, flux de masse (de soluté), flux de chaleur dépendant du
temps
Concentration
Données géométriques du milieu et du maillage

Expériences acquises et validation

Modélisation d’un écoulement en régime transitoire et en régime permanent dans un
milieu poreux à double porosité (Aspo, Suède)
Etude de la migration des déchets radioactifs dans les roches granitiques (Grimsel,
Suisse ;Finnsjon, Suède ;Neuvy Boin, France)
Diffusion des matériaux de déchets contaminés en dépôt à Münchehagen, Allemagne ;
à Mol (Belgique)

Limitations et future évolution

Absence d’interactions physico-chimiques
Ecoulement multiphasique (fluides non miscibles)
Transport et échange de soluté entre fluides non miscibles
Peu adapté à l’étude des transferts dans les sols non saturés
 Nom :DEDALE 3D
Société :LHF
Le Laboratoire d’Hydraulique de France (LHF), le Centre de Recherche pour l’Hydraulique et
les Structures (CRIS) de l’Ente Nazionale per 1’Energia (ENEL) et SOGREAH développent
conjointement un système de modélisation tridimensionnel opérationnel pour les eaux
souterraines : DEDALE 3D

DEDALE 3D est un outil de simulation des écoulements, du transport et de la dispersion dans

les eaux souterraines. Parce qu’il modélise à la fois la zone saturée et la zone non saturée,
qu’il tient compte de l’influence sur l’écoulement des variations de la densité de fluide, qu’il
permet d’intégrer les hétérogénéités et anisotropies du sol et les singularités hydrauliques
(puis tranchées, drains, surfaces de suintement...), le système DEDALE est adapté à toute une
gamme de problèmes posés par la gestion ou la protection de la ressource en eau souterraine
ou des problèmes de génie civil.

DEDALE est applicable aux études de pollution des eaux souterraines comme :
- la contamination provenant de décharges,
- fuite de stockage souterrain profond ou de surface,
- pollution agricole par pesticides et fertilisants,
- pollution ponctuelle et diffuse,
- périmètre de protection,
- pénétration de solution salée...
ainsi qu’à des problèmes de génie civil :
- excavation drainage,
- tunnels,
- construction de barrage ...
Ce système de calcul, débuté en 1992, a déjà fait l’objet de nombreuses validations de terrain
et continue à évoluer. DEDALE peut tourner sur des micro-ordinateurs ou sur station.

Formulation des équations de bases

On suppose que la formulation de Darcy est applicable. Par conséquent les équations
générales sont donc les suivantes :

- Loi de Darcy avec densité variable du fluide:

q = -K(VH + p’y)

- Conservation de la masse de soluté avec une cinétique de dégradation du premier ordre :

R(- a ( P W - =O
- p e K ) + div[pCi - pûDgradC]
at

- Loi de variation de la densité du fluide en fonction de la teneur en soluté :

p = f l(C) (3)
- Loi de variation de la conductivité hydraulique en fonction de la teneur en eau:
K =tqe) (4)

- Loi de variation de la succion capillaire en fonction de la teneur en eau :

h+d(0)

avec :
q vecteur vitesse d’écoulement
H charge hydraulique
C concentration du soluté
p masse volumique de l’eau avec soluté
8 teneureneau
K conductivité hydraulique
z côte par rapport à un niveau de référence (croissant vers le bas)
D tenseur de d i f i i o n
R facteur de retard
h coefficient de dégradation
y vecteur unilatéral
p’densité relative : p’= -P -1
Po
avec po densité de l’eau à 10°C.
L’aspect biogéochimique est transcrit par un facteur de retard R et un coefficient de
dégradation h. Cette approche très classique se retrouve dans la plupart des modèles actuels.

Le schéma numérique

La méthode de résolution est fondée sur la méthode des différences finies. Cette méthode qui
présente deux niveaux d’itérations a été préférée à la résolution plus classique de l’équation
(2) exprimée sous la forme de l’équation de Richards car elle permet une meilleure
conservation de la masse d’eau notamment dans la zone saturée. L’équation de transport et de
dispersion de soluté est à son tour, discrétisée selon un schéma implicite aux différences
finies. Ce schéma a l’avantage de fiabiliser la résolution sans pour autant introduire une forte
diffusion numérique.

Les caractéristiques de sol et les conditions limites

Le calcul requiert la connaissance de la succion capillaire et de la conductivité hydraulique en

fonction de la teneur en eau. DEDALE permet de définir ces courbes caractéristiques en
chaque point du domaine ou dans des zones pour représenter les hétérogénéités du sol. Des
courbes différentes peuvent être obtenues dans les directions verticales et horizontales pour
tenir compte d’une éventuelle anisotropie.
En ce qui concerne les conditions limites, DEDALE gèrent toute une gamme de conditions
aux limites externes (flux ou charge imposée, concentration imposée ou sortie libre) ou
internes (surfaces de suintement, puits de pompage et d’injection, tranchées, drains, parois
minces peu imperméables ou imperméables).
 Nom :SIMüSCOPP
Société :IFP
Le logiciel SIMUSCOPP est développé dans le cadre du projet EUREKA RESCOPP par
l’Institut Français du Pétrole (développement du simulateur réalisé par Beicip-Franlab,
validation), BURGEAP (spéciations, validations) et des sociétés italiennes du Groupe ENI
(développement des interfaces).
Ce modèle numérique a pour but la simulation :
+ de la migration d’un polluant hydrocarbure ou solvant chlorés de la surface d’un sol
au travers de la zone insaturée et saturée, sa dissolution puis son transport dans la nappe (aide
à l’analyse des risques),
+ d’un certain nombre de procédés de remédiation :pompage d’huile mobile, venting
(froid ou chaud), biodégradation, injection de solutions tensioactives, injection de vapeur.

Ce logiciel peut travailler en lD, 2D ou 3D en géométrie cartésienne ou radiale, sous

maillages locaux, en milieu poreux homogène ou hétérogène (porosités et conductivités
hydrauliques).

Les phénomènes physiques pris en compte sont les suivants :

- 1,2 ou 3 phases présentes dans le système eau-phase organique (huile-air),
- écoulements mono, di ou tri phasiques (loi de Darcy généralisée),
- transfert de constituants entre les phases : dissolution dans l’eau, évaporation en
phase gazeuse,
- transport par convection et diffusioddispersionhydrodynamique,
- facteur retard : rétention de polluants,
- calcul des propriétés thermodynamiques et d’équilibre des phases à l’aide
d’équations d’état ou de valeurs tabulées (base de données),
- biodégradation des produits dissous et en phase résiduelle (en cours de
développement).

Les conditions initiales et aux limites prises en compte sont les suivantes :
- équilibres gravitaires et capillaires en zone insturée (eau-air ou eau-air-huile), avec
ou sans pendage de nappe,
- conditions aux limites :
. nappe : surfaces d’injection et productrices d’eau, recharges, pluies, puits
d’injectiodpompage,
. polluant : injection en surface, pompage en zone insaturée ou saturée,
. échanges avec la surface du sol : entréedsorties d‘air en fonction des gradients
de pression (naturels ou imposés) dans le milieu poreux.
Chaque condition aux limites peut fonctionner en pression ou en débit et supporte la présence
d’une ou plusieurs phases (par ex. le pompage d’huile et d’eau dans un puits de rabattement).

Aspects numériques :
- discrétisation spatiale et temporelle par volumes finis,
- schémas explicites ou implicites en temps, linéarisation par méthode de Newton,
- résolution des systèmes linéaires par méthodes de GAUSS ou par méthodes itératives
(gradient conjugué préconditionné ou non).
Les validations qui ont été effectuées se sont faites sur colonnes en laboratoires, à l'échelle
d'un site pilote (SCERES) ainsi qu'in situ sur des problèmes de pompage d'huile mobile,
d'opération de venting ou de transferts d'hydrocarbures dans le cadre d'accident.
Le code écrit en FORTRAN 77 fonctionne sur un grand nombre d'ordinateurs : PC, Stations
de travail UNIX, ordinateurs centraux vectoriels ou non (CRAY, VAX, CONVEX...) .

Formulation des équations de base :

La description mathématique du "flux multiphasique composite" en milieu poreux est basée

sur l'équation de conservation de masse pour chaque composé, et de conservation des
moments pour chaque phase. L'équation générale de conservation de masse pour un composé i
existant dans les trois phases fluides (p) et adsorbé sur la matrice r est :

avec :
g phase gazeuse
i composé
O phase huile
P phase fluide (g, O, w)
Y phase solide
W phase aqueuse

ip flux de diffusion-dispersion de i dans la phase p , mol.s-''

P pression de phase, Pa
9: fraction molaire de i adsorbé sur la phase solide
Qi terme source d'injectiodproduction du composé i, mol.s-'
Ri terme source de biodégradation du composé i, mo1.s-'
SP saturation de la phasep
t temps
- vélocité de Darcy de la phasep, m.s-'
uP
Xp fraction molaire de i dans la phase p
ap flux modificateur de i dans la phase i
4 porosité
PP densité molaire de la phase p , m01.m'~
 Nom :SUTRA
SUTRA (Saturated, Unsaturated, TRAnsport) est un logiciel qui permet de simuler les
mouvements d'eau et le transport soit d'énergie thermique soit d'éléments dissous, entraînés
par cette eau, dans un milieu poreux continu.
Le modèle est basé sur un méthode numérique hybride combinant les éléments finis et les
différences finies pour discrétiser les équations qui gouvernent les phénomènes
interdépendants qui sont simulés :

1) Le mouvement d'un fluide aqueux, de densité pouvant être variable, en zone saturée ou non
saturée.
2a) Le transport d'éléments dissous dans l'eau qui peuvent être en outre soumis au cours du
mouvement à des phénomènes d'adsorption sur la matière poreuse, et à des phénomènes de
croissance ou de production à la suite de réactions géochimiques locales.
2b) Le transfert d'énergie thermique par l'intermédiaire du fluide en circulation et de la matrice
poreuse (par conduction thermique).

SUTRA donne comme résultats fondamentaux des simulations : la pression du fluide en

mouvement et soit, selon l'option choisie, la concentration des éléments dissous, ou la
température, en fonction du temps et en tout point du domaine représenté.
SUTRA peut évidemment être utilisé pour simuler tout phénomène plus simple inclus dans
l'ensemble défini ci-dessus.
SUTRA est essentiellemnt bidimensionnel, ce qui veut dire qu'il peut être utilisé pour
représenter des coupes verticales dans le cas des phénomènes où intervient la gravité
(écoulement en zone non saturée, écoulement saturé à densité variable), ou pour représenter
des domaines en projection sur un plan horizontal dès lors que la gravité n'intervient pas. Dans
ce dernier cas, la 3emedimension est prise en compte par intégration, dans la valeur de certains
paramètres physiques qui sont fonction de l'épaisseur du milieu aquifêre (par exemple la
transmissivité) ou de distances verticales entre ce milieu et certaines conditions composées
(ex : la drainance).
SUTRA utilise, pour la discrétisation spatiale, la méthode des éléments fuiis. Ce qui veut dire
que la conception du maillage est très souple, à base d'éléments quadrilatéraux. En outre,
SUTRA autorise la création de "nœuds de pincement'' (pinch nodes) sur n'importe quel côté
des quadrilatères du maillage. Ce qui permet de resserrer très facilement le maillage dans tous
les endroits où la dynamique des phénomènes est importante (et réciproquement de le relacher
dans les endroits plus calmes).
Pour la discrétisation du temps, SUTRA utilise une approximation par différences finies
implicite avec prévision ou plus rarement un schéma explicite dans les débuts des périodes
transitoires.
Les conditions aux limites prises en compte peuvent être permanentes ou variables dans le
temps. Il peut s'agir de conditions de charge de concentration ou de température composées de
débits en eau evou d'éléments dissous ou de calories.
Les équations fondamentales qui sont traitées par le modèle sont :
1) les équations de conservation de la masse et de l'énergie,
2) l'équation de DARCY généralisée pour traduire l'équilibre hydrodynamique de la phase
fluide,
3) l'équation de dispersion des éléments dissous ou de la température au cours de leur
entraînement par l'eau dans le milieu poreux,
4) éventuellement une loi d'adsorption ou une loi de production (ou disparition) - par réaction
géochimique avec la matrice - de l'élément dissous dont la concentration est considérée.

Ces équations sont couplées dans la mesure où certains des paramètres physiques qui
interviennent dans chacune d'elles sont des fonctions de variables évolutives en fonction des
autres. Par exemple, la perméabilité à l'eau est une fonction de la température (par
l'intermédiaire de la viscosité) et de la densité ;la densité peut dépendre de la concentration en
éléments dissous, etc.
Les résultats des simulations sont donnés sous forme de tableaux de valeurs pour chaque
nœud du maillage ou sous forme de courbes isovaleurs tracées par un module de contourage
intégré dans SUTRA.
SUTRA est sans doute un modèle très complet du point de vue théorique. Compte tenue de sa
complexité, il doit être réservé à l'analyse de cas bien circonscrits, si possible par l'utilisation
préalable de modèles plus simplifiés.
ANNEXE 6
 v> f-
m r-
O
cn N

O
8
s 2
m

-L
8 8
0-
f-
(v
+
u>
O

+
œ
L
O I
n

8 8
a
m c
O
:
c
8
2
r-

8 8
N -
O'
O> -
el
O

8
-4
O
ANNEXE 7
consomum

- -
8
- _.-..
--:c:

dl ’ Research Consortium

The Hydrogeology Program at the University of Illinois is an active group that develops
and apphes computational methods to investigate hydrogeologic and geochemical
processes in the subsurface. The Program has more than ten years of histoy in cooperative
research with domestic and international research groups fiom the petroleum, minerals,
and environmental protection industries, as well as govemment research laboratones and
other universities.

An industrial consortium was created in 1986 to facilitate such cooperatioa The current
consortium members are:

1 Amoco Production Research Company

1 ARC0 Oil and Gas Company
Chevron Oil Field Research
1 Conoco,Inc.
1 Exxon Production Research Company
Hewlett-Packard
Japan National Oil Corporation
1 Lawrence Livermore National Laboratorv
1 Marathon Oil Company
1 Mobil Research and Development
1 Sandia National Laboratorv
1 SiliconGrarhics Computer Systems
1 Texaco,USA
1 U.S. Geological Survev (Reston)
iU.S. Geological Survev (Denver)
As OUT research sponsors, consortium members enjoy a number of benefits, including
licenses to software developed with the Hydrogeology Program.

-
hvdro hornmaoe - consomum software - short course- uubiications

Last updated 9/13/96:mrt

ICTML written using HTML Pro and testeci on Netscape Navigator 3.0 for Mac and Microsofi Intemet
Explorer 2.0 for Mac. Ail graphic images were manipulated using Adobe Photoshop. This page is maintained
on a semer at the University of Illuiois at Urbana-Champaign
URL: “http://www.hercules.geology.uiuc.edu/-bethke/consort.hrnil”

Lof1 09/27/96 17:jZ:Z 1

 _. . . ....
.
.-:,..-i .; . .
-- . . -+
6 .

- ,,
.* --:
f ~ . People in the Program

scraig Bethke, Professor, c-bethkeeuiuc.edu

-&mette Fud, Research Scientist,
[email protected]

3Tren Haselton, Research Programmer,

[email protected]

-2ilar Manzano, Research Programmer,

[email protected]

.@mg Zhao, Research Associate,

[email protected]

+mv Berger, Graduate student,

[email protected]

-elinda Tidrick, Graduate student,

[email protected]

;;*Jun&o Park,Graduate student,

[email protected]

+George Roadcap, Graduate student,

[email protected]

s h i y h g Cao, Graduate student,

[email protected]

~ T a m aShinar,
r Student programmer,
[email protected]

+Jean Daly, Business Manager,

jdalyeaquifer .geology.uiuc.edu

hvdro hommage - consortium - software - short course- publications

Lasî updated 919196:mrt

1 of2 09/27/96 175053

ANNEXE 8
QUESTOR http://web.qub.ac.uk/qc/questor.htn

The Oueen’s Universitv of Belfast

Questor Centre WWW Server

THE QUEEEPS u”ERsrry E”MENTALSCIENCE AND TECHNOLOGY RESEARCH CENTRE

Ftroductionl Fdustrv Membersl rCenter Oreanisationl Benefits to uIdu-1 rcment Research Proiectsl@&J
Advertisea JOUCATJ

+44 (0)1232 33557718 +44 (0)1232 661462

Emaii address: QUESTOR@,aub.ac.uk

INTRODUCTION

The Questor Centre was the first hdustry-miversitycooperative environmentai research centre outside the USA. It was set up
in May 1989 in The Queen’s University of Belfast with nine founder industry members and five participating University
departments. The Centre is physicaily located near the School of Chemistry which initiated the development.

-
The Centre aliows both partners to do that which they are best at the University to carry out basic research, the indusûy
partners to apply this research to minimise their environmental impact while generating wealth. This real partnership between
industry and the University is very advauîageous. It allows input fiom indusûy into the choice of research topics and helps to
focus the relevance of the research. The research students and post doctoral staffworking in the Centre are in regular contact
with senior industrialists Who can advise on the relevance of the work. The industry links help ensure that the students can
integrate more smoothly into the work environment.

INDUSTRY MEMûERS

The members are drawn fiom Ireland, North and South, and fiom GB. They are currently:

CI Beckman Instruments

D BNFL

0 BP Chemicals

0 Chemvite

C Coca Cola Bottiers (Ulster) Ltd

CI DoE Environment Service

3 DoE WaterService

Du Pont (UK)Ltd

3 Exxon Research & Engineering

1 of 3 10/17/96 11:41:4
QUESTOR http://web.qub.ac.uk/qc/questor.htm

Z Fisons Instruments

O Fruit of the Loom

[3 International Biochemicals

O ICI Chemicak & Polymers

O Northern Ireland Elemicity

0 Perkin-Elmer/Foss Electric

O Schering Plough

O Sevm Trent Water

O UnileverResemch

O Welsh Water

0 Yamanouchi Ireland

3 Yorkshire Water

CENTER ORGANISATION

C m n t l y the Departments of Agriculture, Chemical Engineeriug, Chemisîry, Civil Engineering, CornDuter Science, the
School of Biology and Biochemisûy, Division of Molecular Biology, and Psychology are participating in the Centre research
programme.

The smngth of the Centre lies in a carefuriy stnictured interaction between indutry and the University. One represenîative
h m each member organisation sits on the Centre Indusûy Advisos. Board which meets every six months. All proposed
research projects are presented to the Board for selection based on indusûy relevance, quality of the research and abiiity to
achieve the goals set out.

A mixture of research sîudents, research associates and research fellows are carrying out the research and they are selected
h m a large field of Young people Who are atûacted to apply to the Centre by a combination of environmental relevance and
industrial participation in the research. These Young people leave the Centre after two or îhree years and, in today's situation of
increasing demand and decreasing supply of science graduates, they are a valuable human resource. Members of the Centre
have the opportunity to assess theb performance and thus a m suitable candidates to work with them.

The interdisciplinary nature of the Centre, bringing various disciplines to bear on environmental problems, has proved to be
very important. This has, for example, led to a joint project between Chemical Engineering and Microbiology which has been
funded by EPSRC. This project has also been linked with a complementary one in the Hazardous Substance Management
Research Centre in New Jersev. USA. A M e r chemistrykhemical engineering project on flocculation has been fimded by
IRTU and this has been linked with a research project on the computational fluid dynamics of mixing in HSMRC.
Currentiy, the Centre is expandhg to cany out research into clean technology and build a demonstration facility in addition to
the current remit of clean-up and end-of-pipe technology. This expansion will include work on inorganic and organic clean
technology, some of which will be catalysis based, and also biocatalytic work. 'The staff of the Centre wiii number 36.

a SUMMARY OF BENEFITS TO INDUSTRY

Prior access to a waste management, clean technology and demonsiration research programme fimded at over
&pound;6m over the first six years.

m The opportunity to help select the research to be carried out and review it every six months.
The opportunity to assess and eventually offer employnent to the graduates and postgraduates working in the
Centre.

m Royalty-fiee access to al1 patents arising ftom the research.

2 of 3 10/17/96 11:41:5
QUESTOR http://web.qub.ac.uk/qc/questor.htm

1)Close contact with the Centre and thereby with The Queen'sUniversity of Belfast.

Involvement in an international environmentalresearchprogramme between USA, France, Germany, Russia and

Questor.

The opportunity for networlcing withh the Membership and Internationally.

For further information piease contact:

Prof. Jim Swindall
Director
The Questor Centre
The Q u m ' s University of Belfast
Belfast B'i9 5AG

Last upd&d 10 ûctober 1996 E-Mail [email protected]

3 of 3 10/17/96 11:41:!
 Pre-proposai for a Questor Centre Research Project
September 1994

1. TitLe
Numex?cal m o d e h g ofpoUutant dispersion in the unsatumed (surface) zone.
2. Principal UNestigator
John Stewart, QuestorKomputer Science.
Ph. (0232)245133 Ext 3375 Fax (0232) 33 1232
3. Oventiew
Many of Questois members may have pmbiems or potential probiems with contaminated
land. Modelling studïes are usuaUy necessary to determine the short and long term
dispersion of contaminanu and any ensuing risks and to decide on the need (and potentiai
strategies) for r e m a o n .
There are many commercial and research codes availabie to undercalce such studies and
m o d e h g of hydrology and contaminant aan~pofiin the (genedy) deeper saairated zone
has been wideiy undenaken. However, modelling of the unsaturated (nearsurface) zone,
and particuiarly the use of coupied modeis, which cope adequateiy with migration h m .
unsaturated ro saturated layes, is not so Weil advanced because of some wideiy reco-&ed
dBcuities.
These difnculties arise when attempting to soive the non-hear systems of quafions
needed to describe the compiex processes in the unsarurated zone. In parti&, the
presence of 'sharp fronts'associated with mi@aung contaminants causes probiems for
modeilers. These problems are d e s t e d in two main ways, either the mnrhems\n'Cal
soiver fails to converge to a stable solution or oscillahons about a solution OCCUII m e n
coupied with chemicai reacrion modeis, oscillations to negative contaminant concentrarionS.
will obviousiy iead tc p q p m fadure.
These difnculties were noted in a major comparafive smdy of gound water modefling
suategies carried out by the OECD (The Intemational HYDROCOIN Project) where none
ofthe six projet teams were able to obtain a converged soiution using the paxiimeters
suppiied with a test case. and have been connmied by researchers at British Nuciear Fueis
pic (A. Pinner, pers. CO-) Who are currently deveioping a modei for the unsaairated
zone.
This project seeks to combine the skius and experience in High Performance C o m & n g
and mathematical methods availabie at Queen's with the mo&Uing skills of xesearchers
from ttiose Queszor members Who are working in this fieid to overcome these diffïculties.
Some potential soiutions have been identified and these are deta3ed beiow.
4. IndustMl relèvame
Modelling dispersion (and retention) processes in the unsatmated zone is important
- for
severai reasons,
i) spiiis and accidental reieases tend to occut at the ground surface (Le. in the unsaturated
zone) and so retention in ancVor transport through the unsaturated zone is of prime
importance,
ii) predictions of where a contaminant wiü be retained are ais0 important when planning a
remediation strategy. If the contaminant is ansported to the sanirated zone (i.e. into the
aquifer) it may be possible to remove it by pumping h m weiis. If contaminants remain in
the unsaturated zone they may be difncult to remove without changing the hydrogedogical
characteristics of the system (Excavation is currently the most common strategy, but
others inciuding the use of leachhg soiutions and pumping Wells may be possr'bie.)

Queen's Univemty Environmentai Science and Technology Research Centre

iii) assessmentS of the long-tem risk associated with contaminanCs retained in the
u n s a m t e d zone will help to decide whether or not to undertake remedial action.
The widesprcad use of computer modeis for environmentai m o d e h g indicaies their c o s
effdveness. It would be unreaiistic to form a compiete chamcterisation of ground
contamixmion by physical examination. This panicular project is ainied at an important
-
area near-surface contaminaton - where many sites cannot be adequately modelied.

5. Route to exploitdon
The project wodd be carried out in cosperation with researchers h r n those Questor
members Who wish to take parf. The work shodd provide new modelling codes, with
acknowledged effectiveness, for use both in-house and commerclally. In order to estaùiish
the value of the mode1 detailed studies of conraminated sites would be made and it is
intended that Questor members would find direct benefit by co-operating in studies of their
own sites.
6. Metltodology
The project would proceed through a number of stages.
i) A review of available techniques and models (pdcularly those which have been
demonstrated, tested, calibrated and fully documented)
ii) Improving the performance of existing models. Some of the difficulties o u h e d above
may be ameiiorated by r e m g the three dimensional finite element mesh, for example
using frner meshes, moving or adapting meshes. for regions where flow and concentrarion
gradients change rapidiy. Mesh refiement techniques impiy the need for much increased
cornpuMg powcr and the use of Hi@ Performance Computing such as parallel or
mainfnme machines will be investigated
iii) An aitemative approach involves the investigation of different mathematical solvers for
the finite element scheme which mi@t circumvent some of the problems of non-
convergence or oscillations. It is possible that researchers in related fields (e.g
aeronaurics) have developed techniques which codd be adapted to overcome some of the
diffcuities.
iv) Having idenùfled andor produced candidate model(s) for unsaturated zone problems
these wiU be tested and caiiïrated at sites repsenting a range of conditions. Questor
members would be requested to take part in this study by tesMg codes against their own
contaminaed sites.
v) In Line with the aim of producing modelling codes which are acceptable to regdatory
authorities and others, details of the study will be documented for scrutiny and approval by
interested parties.

7. Objective
The objective of the project is to arrive at a consensus on which codes are usefui for
modelling in the unsatura~ed zone and codd be taken fomard for regulatory examination
and approval and subsequently be used both in-house and commercially.

Queen's University Environmental Science and Technology Research Centre

8. Timing

MONFHÇ

Review

tmpmved perfomian-
Alternative soivers
Testicalibrate
--
- 6 12

I
18 24 -c
.a& 36

Doamrent

9. Budget
This would be the main îask of the phcipal kvestigator but, in addition. a-mh
assistant, preferabiy with a back-mund in mathematicai computation wouk ais0 work on
the project.
Travei to base with co-operating rnembers would be necessary.
in addition, there would be coss associated with caiibration studies of cant;c,te codes -
these costs would not be part of the project budget,

Cost surmnary Staff 2 3 Senior Research FeUow

2/3 Postdoc. Research Assistant
Travei Bo00

Queen's University Environmental Science and Technology Research Centre

 Pre-proposal for a Questor Project

A CAST award studentship to s t u d y biogeochemical modelling of the fate of metai

contaminanis in saturated and unsaturated media

Pmject Title
Computer simulation of coupled transport and chemical interaction between organic
compowids and metais in the environment.

fincipal Investigators
This wodd be ajointly supervised CAST award studentship (a CAST award is the NI
equivalent of a CASE award in Britain). The Queen's University supervisors wouid be Dr.
RM. Kalin and Mr. J.R. Stewart. The Tndustrial partner would be British Nuclear Fuels.

Overview
The capability to model the movement of metai pollutants in saturated and unsaturated media
is limited by OUT ability to include ail the relevant physical and chemicai controls into
biogeochemical computer codes. In particuiar, the movement of metais is not oniy controlled
by reduction-oxidation reactions and the solubility of inorganic saits, but also interactions with
naturally occurring or anthropogenicaily derived organic compounds through chelation or
other bonding mechanisms. M e n bound to organic compounds their subsequent transport
may be retardeci, if the organic material is essentially immobile, or modified depending on the
transport of the organic.

This project will add to the ongoing groundwater modelling efforts at The Queen's University.
There is currently one Questor project on groundwater modelling. John Stewart 1s workmg on
the development of improved numerical methods to solve problems encountered wheb
modelling the transport of contaminants in the unsaturated zone, especially where 'sharp
concentration fronts' occur. Dr. Robert Kalin has several post-graduate students working on
various aspects of groundwater modelling and contamination. His new facility, the
Environmentai Engineering Research Centre VERC), was initiated to offer novel techmques
to vaiidate isotope-geochemical and biogeochemical modelling of contaminated land and will
work closely with the Questor centre. This project wilI form a usefid link between an
existing Questor research programme and h s new research group at Queen's.

Industrial Relevmce
The fate of metai contaminants in water is of considerable interest to both the subsequent
users of potentially contaminated water and fro those deemed responsible for the
contamination. Adequate computer simulation can provide a quantifiable projection of the
fate of metals in the environment. These models can ais0 be used as a predictive tool for the
potential success of remedial strategies.
Objectives
The objectives of thrs study are as follows:

1) Conduct a review of the literature on interactions between organic compounds (of specific
interest) and metais in the saturated and unsaîurated environment.

2) Study and develop the chemicai thermodynamic and kinetic relationships between
organics and m d s of interest needed for incorporation into existing computer modelling
code.

3) Design and Wnte appropriate computer code for incorporation into biogeochemical
modelling programs

Methodology
The course of the project will be dtiven by the availability of appropriate thennodynamic and
lunetic data for interaction between metals and organic materiais. We anticipate much o f dus
data currently exists due to research efforts during the past decade in the USA which
considered metal-organic interaction as a method of extracting metais from mine tailings and
waste sites. Some research at Queen's has been published on the formation of complexes
between copper ions and nat-ally occurring humic and fuivic acids. This work of Dr. Allen
from the Department of Chemical Engmeering may provide additional data. We anticipate
that the experimentai work will be undertaken jointly between EERC labs and the Questor
centre.

The project proceed viai the following stages

1) A literature survey will be made for details of reaction kinetics and equilibrium data on
relevant metal-orgauic compound interactions. The extent of the data available will determine
the need for experimentai work.

2) Experimental work will be undertake to sturiy. the interactions between specific meîals and
organics of interest. Emphasis will be given to *ose combinations of interest to Questor
-members.

3) Computer modelling code will be writîen to incorporate the effects of these chemicd
interactions into contaminant transport models.

4) The new code will. be incorporated into a current complete biogeochemical model, tested
and calibrated using both test and field data
Timing
The project will run for the normal three year period associated with doctoral research. It is
anticipated that the student would work for approxirnately half of the first two years at British
Nuclear Fuels research laboratory.

n i e proposed timing is as follows:

Months L ocation/w O&

1 to 3 -
BNFL familiarizalion wiîh biogeochemical modelling techniques 'and
computer codes. Literature survey on metal-organic compound
interaction

4 to 15 QUB - EERC IaboratoqdQuestor centre. Experhental work on reaction

b e t i c s and equilibria Initial design and coding of routines for
biogeochemical mode1 subroutuie.

16-2 1 -
BNFL further coding of computer routines and incorporation into
modelling program. Testing and calibration of new routines using
appropriate data from contaminated sites(s).

22-36 -
QUB remaining work on computer codes, program calibration and
writing üp thesis.

Dependhg on the avatlability of pubiished mformation on metal-organic compound reaction

kinetics and equilibria, there will be a greater or lesser need for experimental work in the
Iaboratory.

Budget
The student would be supported by a DENI CAST award. A supplement of €1000 (or BO00
for those with a first class honours degree) will be paid.
Résumé :
Les sociétés industrielles utilisent et produisent un très grand nombre de substances qui sont
susceptibles d’engendrer une pollution des milieux naturels et notamment du sol, sous-sol et nappes
phréatiques. Ii est donc nécessaire d’estimer l’impact d’une poilution accidentelle ou chronique,
réelle ou potentieiie, en terme de risque. La modélisation, en permettant la simulation du devenir
des polluants dans les sols peut constituer un outil d’aide à la décision, soit en vue de la prévention
de la pollution de la ressource ou de la répmtion de nuisances déjà existantes.
Dans ce rapport, nous présentons dans un premier temps l’état de l’art sur la modélisation de
ûansfert de poiiuats dans ies sols. Nous abordons ensuite les principaux outils existants SUT le
marché à l’heure actuelle, et les principaux organismes et laboratoires travaillant dans ce domaine,
en France tout particulierement.
A l’issue de cette étude, il est possible de dégager les principales questions et les problèmes qui
restent posés concernant la modélisation du transport des p o i i ~ t sdans les sols et qui pourrront
fêire l’objet d’études ultérieures. Quatre points se sont dégagés :
- la sensibilité des modèles aux différents paramètres d’entrée et la hiérarchisation des
paramètres clés,
- la réalisation d’me base de dom& faisant l’état de l’art des études modélisation rialisées
dans le cadre des sites pollués,
- la mise au point d’iidicateurs de mobilité d’un polluant dans un sol,
- le domaine d’application et la pertinence des modèles.

Summary :
Industriai wmpanies use and produce numerous substances which can induce a poilution of our
enwOnment and speciaiiy of soil and groundwatex. Thus, it’s necessary to estimate the risk of an
environmental impact h m an accidmîal or chronic, reai or potential pollution. Modeliing, which
ailow the simulation of poiiutant migration, can be used as a decision support system, either for the
pollution control and prevention of the resource, or for the monitoring of the remdation of
polluted sites. 6
in the first part of îhis study, we established a state of the art on modeliing of pollutant migration in
soils. In the second part, we f w e d on the main simulation tools currently available on the market,
and on the main agencies or laboratories working on this subject, especially in France.
At the end of this study, we àrew some conclusions concerning modelling of pollutant kgration in
soils and the mainpoints which will f m the subject of furîher studies :
- sensitivity analysis of mode1 to input parameters and rankuig of the main parameters,
- achievement of a database on the state of the art of the results on modelhg realized on
case studies,
- development of a mobility indicator of pollutant in soil,
- application field and relevance of the models.

Mots clés :sois pollués, poiiuants, migration, modélisation, simulation, logiciels.

Key words :poiiuted soiis, poliutants. mimition, modeiiing, Simulation, software.