II.
Dispersion atmosphrique
1/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Organisation de la journe
Les principaux modes de rejets
Terme source
Rejet
Evaporation de nappe
Dispersion du nuage
Transport / Diffusion
Conditions mtorologiques
Conditions environnementales
Modlisation de la dispersion
Effets toxiques
Notion de seuil toxique
Modlisation
2/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Organisation de la journe
Les principaux modes de rejets
Terme source
Rejet
Evaporation de nappe
Dispersion du nuage
Transport / Diffusion
Conditions mtorologiques
Conditions environnementales
Modlisation
Effets toxiques
Notion de seuil toxique
Modlisation
3/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Les diffrents modes de stockage
Etat du fluide
pression et
temprature
ambiante
Etat de stockage du fluide
liquide
liquide
gaz
gaz
liquide
pressuris
cryognique
4/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Les principaux modes de rejet dans latmosphre
Schma de principe
Fumes
dincendie
+
Rejet de gaz toxique et/ou inflammable
Instantan
continu
Exemple:
Ruine de rservoir
Stationnaire
Non stationnaire
Exemple:
Exemple:
Brche
Rupture guillotine
5/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Organisation de la journe
Les principaux modes de rejets
Terme source
Rejet dun fluide latmosphre
Evaporation de nappe
Dispersion du nuage
Transport / Diffusion
Conditions mtorologiques
Conditions environnementales
Modlisation
Effets toxiques
Notion de seuil toxique
Modlisation
6/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Phnomnologie (1)
Brche ou rupture de canalisation (gaz liqufi sous pression)
Phase d expansion
Flash
Phase d entranement
Arosols
+ vapeur
Brche
Rain-out
vaporation
7/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Zone passive
Phnomnologie (2)
Phase dexpansion
Flash thermodynamique dtente
Vaporisation
Fractionnement
Phase dentrainement
Condensation des molcules deau contenues dans lair
Phase passive
Vitesse du jet = vitesse du vent
8/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Phnomnologie (3)
Exemple: Rejets de butane en prsence dobstacles
Obstacle
Accumulation de glace et de butane :
- sur lobstacle
- sur les sondes
Produit pig dans
le bac de rtention
FILM
Bac de rtention
9/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Phnomnologie (4)
Mcanismes physiques mis en jeu selon la nature du produit stock
rainout
Rapport -19 de lINERIS
Dtermination des grandeurs
caractristiques du terme source
ncessaires lutilisation dun modle
de dispersion atmosphrique des rejets
accidentels
10/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Les principaux modes de rejet
Stockages pression atmosphrique
- substance liquide temprature ambiante
(pression de vapeur saturante < 1 atm Tamb et Tbul > Tamb - ex. : styrne, xylne,
benzne)
- gaz liqufi par rfrigration (ex. : NH3 -33C, GNL -162C)
REJETS LIQUIDES
Stockages sous pression
- gaz liqufi sous pression (ex : GPL) REJET
DIPHASIQUES
- gaz sous pression REJET GAZEUX
Incendies mettant des fumes
11/72
REJETS GAZEUX + suies
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Modlisation du terme source: rejet 100% liquide
Calcul du dbit de fuite
Pres
Patm
Bernoulli:
2.P
Q .Cd .S . 2.g .h
: masse volumique du liquide (kg/m3)
Cd : coefficient de dcharge (0.62 bord pais / 0.82 bord mince /0.96 bord arrondi)
S : section de l'orifice (m2)
g : 9,81m/s2
h : hauteur de liquide (m)
P : pression relative dans le rservoir (Pres - Patm) en Pascal (Pa)
12/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Modlisation du terme source: rejet 100% gazeux (1)
Calcul du dbit de fuite: coulement non choqu
Pres
Pa
coulement
choqu
1
coulement non
choqu
Pa : Pression atmosphrique (Pa)
P rservoir
Prs : Pression du rservoir (Pa)
P atmosphrique
g : Masse volumique dans le rservoir (kg/m3)
Cd : Coefficient de dcharge
Pa
Brche
13/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
S : Surface de lorifice (m2)
= CP/Cv
Modlisation du terme source: rejet 100% gazeux (2)
Calcul du dbit de fuite: coulement choqu
Pres
Pa
coulement
choqu
1
coulement non
choqu
Pa : Pression atmosphrique (Pa)
Prs : Pression du rservoir (Pa)
P rservoir
g : Masse volumique dans le rservoir (kg/m3)
Cd : Coefficient de dcharge
P atmosphrique
Pa
S : Surface de lorifice (m2)
= CP/Cv
Brche
14/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Modlisation du terme source: pertes de charge
Pertes en ligne
Effets de constriction
Les pertes en ligne dpendent
principalement:
du diamtre des canalisations
de la gomtrie du rseau (coudes,
changements de diamtre, jonctions,
vannes )
Effet de constriction
dpend de la forme de la brche
se traduit par un coefficient de
dcharge Cd
15/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Modlisation du terme source: rejet diphasique (1)
Phnomnologie
Apparition de phnomnes complexes :
Fragmentation mcanique
Fragmentation thermique
Fort cart de dbit par rapport un rejet
purement gazeux ou purement liquide
Formules de dtermination de dbit
prcdentes inutilisables
16/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Modlisation du terme source: rejet diphasique (2)
Calcul du dbit de fuite: plusieurs cas
Rejet essentiellement liquide
Fluide fortement sous-refroidi dans le rservoir (Psto >> Psat)
Formule de Lackm [1]
Rejet intermdiaire
Fluide entre le sous-refroidissement et la saturation (Psto > Psat)
Formule de Fauske [2]
Rejet lquilibre
Fluide proche de la saturation (Psto ~ Psat)
Psto = Pression de stockage (Pa)
Formule de Wheatley [3]
Psat = Pression de vapeur saturante (Pa)
Rf: [1] Thermodynamics of Critical Two-Phase Discharge from Long Pipes of initially Subcooled Water , Heat Transfer in Nuclear Reactors, pp 391-407 1980
[2] Source term consideration in connection with chemical accicents and vapour cloud modelling, J. Loss Prev, Process Ind, Vol 1, p 75, 1988
[3] Discharge of liquid ammonia to moist atmospheres Survey of experimental data and model for estimating initial conditions for dispersion calculations,
SRD/HSE/R410. United Kindom atomic energy 1987
17/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Modlisation du terme source : rejets particuliers
Cas des rejets confins
Paramtres supplmentaires
Dilution du nuage dans le local
Influence de lextraction
Fuites au niveau des parois
Effets retardateurs
obstacle
Cas des rejets en prsence dobstacle
Paramtres supplmentaires
Partie liquide pige
Modification de la concentration
18/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
vaporation
Brche
Capture
Exercice
La rupture guillotine dun bras de dchargement de wagon est envisage. Les 60 tonnes de propylne
ou propne (C3H6) prsentes dans le wagon sont relches travers le bras de dchargement. Aucun
systme de scurit ne fonctionnant, la dure de fuite est indtermine.
35m
2m
1m
Question 1: Dterminer lensemble des paramtres qui vous semblent ncessaires la ralisation de la
modlisation.
Question 2: Calculer le dbit de fuite en partie amont. Le logiciel donne un dbit de 34,2 kg/s. Critiquer
cette valeur obtenue par rapport la configuration.
Question 3: Calculer le dbit de fuite laide du logiciel ALOHA
Data:
Meteorological conditions : Wind speed: 2 m/s (10 meters high), wind direction: West, Clear cloud cover, Air
temperature: 10 C, Humidity: 70%,
Ground roughness: Open country
19/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Organisation de la journe
Les diffrents modes de rejets
Terme source
Rejet
Evaporation de nappe
Dispersion du nuage
Transport / Diffusion
Conditions mtorologiques
Conditions environnementales
Modlisation
Effets toxiques
Notion de seuil toxique
Modlisation
20/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Phnomnologie (1)
Vent
Flux
solaire
vaporation
Surface de la
nappe
Temprature et pression de vapeur saturante du liquide
changes thermiques avec le sol
Evaporation gouverne par 2 phnomnes:
Les transfert de masses
Les transferts de chaleur
21/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Phnomnologie (2)
Evaporation gouverne par les transferts de masse (Tbullition > Tsol)
Vent
Flux
solaire
vaporation
Surface de la
nappe
Temprature et pression de vapeur saturante du liquide
changes thermiques avec le sol
22/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Transferts de masse: modlisation du terme source (1)
Estimation du dbit dvaporation
Corrlation de Mackay et Matsugu
( Evaporation Rates of Liquid Hydrocarbon
Spills on Land and Water , Canadian Journal of Chemical Engineering vol.51, August 1973)
Corrlation du TNO 3.75 (Yellow book, third edition, 1997)
Qvap = C * U0.78 * r - 0.11 * Psat(T)*M/(R*T)
Qvap
Dbit dvaporation surfacique (kg/m/s)
Constante = 5,15.10-3 (-)
Vitesse du vent (m/s)
Rayon de la flaque (m)
Psat
Pression de vapeur saturante du fluide (Pa)
Masse molaire (kg/mol)
Constante des gaz parfaits = 8.314 [Link]-1
Temprature du fluide (K)
23/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Transferts de masse: modlisation du terme source (2)
Limites sur les modles
Dbits dvaporation trs variables dune corrlation
lautre (facteur 2, 3pour une mme configuration)
Ecarts lis aux conditions dessais diffrentes dune corrlation
lautre
Nature du produit test (peu ou trs vaporant)
Conditions extrieures (air libre ou confin)
Vitesse du vent (air libre ou confin)
Taille de la nappe utilise
24/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Phnomnologie (3)
Evaporation gouverne par les transferts de chaleur (Tbullition << Tsol)
Vent
Flux
solaire
vaporation
Surface de la
nappe
Temprature et pression de vapeur saturante du liquide
changes thermiques avec le
sol
25/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Transfert de chaleur: modlisation du terme source
Estimation du dbit dvaporation
Corrlation du TNO 3.74 (Yellow book, third edition, 1997)
Qvap
Qsol
Hvap
avec
Qsol
s (Ts Teb)
as t
Qvap
= Dbit dvaporation (kg/s)
Hvap
= Enthalpie de vaporisation du liquide (J/kg)
Qsol
= Chaleur conduite par le sol (W)
= Conductivit thermique du sol, (J/m.s.K)
Ts
= Temprature du sol (K)
Teb
= Temprature dbullition du liquide (K)
as
= Diffusivit du sol (m/s)
= Temps (s)
26/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Evaporation par transferts de masse
(vent) et par flux solaire ngligeable
Organisation de la journe
Les diffrents modes de rejets
Terme source
Rejet
Evaporation de nappe
Dispersion du nuage en champ lointain
Transport / Diffusion
Conditions mtorologiques
Conditions environnementales
Modlisation
Effets toxiques
Notion de seuil toxique
Modlisation
27/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Dispersion du nuage en champ lointain
volution du nuage dans latmosphre rgie par 2 phnomnes physiques :
Le dplacement du nuage : transport (= advection)
La dilution du nuage : diffusion
Paramtres influenant cette volution
Types de gaz : lourd, lger
Conditions mtorologiques : vitesse du vent, stabilit
atmosphrique
Conditions environnementales : obstacles, rugosit, relief,
occupation du sol
28/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Champ lointain : transport et diffusion dun nuage
Rappel diffusion
1)
2)
3)
29/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Champ lointain: concentration du nuage
Dure/type du rejet
Rejet instantane
c (kg/m3)
Rejet continu
A
B
B
t (s)
30/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Champ lointain : types de gaz
Gaz passif
Gaz lourd
Gaz lger
Gaz de mme densit que
lair
Gaz trs dilu dans lair
Pas deffet de jet
(vaporation)
Gaz de forte densit Tamb
(Cl2, COCl2 )
Rejet diphasique (GPL
pressuris )
Gaz froid
Gaz de densit plus faible
que lair (H2, CH4, NH3
gazeux )
Gaz chaud (fumes
dincendie )
Transport la vitesse du
vent
Dilution par la turbulence
atmosphrique
Effets de gravit
Effondrement
Dilution latrale accrue
Frottement la surface du
sol
Fortes concentrations la
surface du sol
Elvation
Dilution augmente par
cisaillement du vent
Concentrations plus faibles
la surface
31/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Champ lointain : gaz dense
32/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Organisation de la journe
Les diffrents modes de rejets
Terme source
Rejet
Evaporation de nappe
Dispersion du nuage
Transport / Diffusion
Conditions mtorologiques
Conditions environnementales
Modlisation
Effets toxiques
Notion de seuil toxique
Modlisation
33/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Conditions mtorologiques: la vitesse du vent
Impact du vent sur le panache
3 m/s
5 m/s
10 m/s
34/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Conditions mtorologiques: la stabilit atmosphrique (1)
Rappel
Patm diminue avec l altitude
Dplacement vertical d une masse d air
Elle se dilate
Elle se refroidit (PV = nRT)
Taux de refroidissement - Gradient thermique adiabatique :
Air sec = 0,98 C / 100 m
Air humide = 0,55 C / 100 m
35/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Conditions mtorologiques: la stabilit atmosphrique (2)
Graphe dvolution de la temprature dans latmosphre
h
100 m
Adiabatique sche
Atmosphre
stable
Atmosphre
instable
OMEGA 12
T
-1C
Ti
Tpol
36/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Ts
Conditions mtorologiques: la stabilit atmosphrique (3)
37/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Conditions mtorologiques: la stabilit atmosphrique (4)
Classification
38/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Conditions mtorologiques: la stabilit atmosphrique (5)
Rglementation franaise
Rejet horizontal au
niveau du sol
D - 5 m/s
F - 3 m/s
ABRejet vertical ou rejet C Dde gaz lger
EF -
3 m/s
3 m/s , B - 5 m/s
5 m/s, C - 10 m/s
5 m/s, D - 10 m/s
3 m/s
Circulaire du 10 mai
3 m/s
2010
39/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Phast v6.51
Organisation de la journe
Les diffrents modes de rejets
Terme source
Rejet
Evaporation de nappe
Dispersion du nuage
Transport / Diffusion
Conditions mtorologiques
Conditions environnementales
Modlisation
Effets toxiques
Notion de seuil toxique
Modlisation
40/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Obstacle
Influence dun obstacle sur un rejet atmosphrique
FILM
41/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Rugosit du sol
Influence de la rugosit du sol sur la vitesse du vent
Dsert
Cultures
Ville
Altitude en m
10 m
0m
Module de la vitesse du vent en m/s
42/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Centre urbain
Conditions orographiques
Influence du relief sur le panache
43/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Brise littorale
Influence de la mer sur les courants dair
Turbulences d origine thermique
Air froid descendant
Air chaud montant
Terre
Brise de mer
Brise de mer
Brise de terre
Air froid descendant
Air chaud montant
Terre
44/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Brise de terre
Organisation de la journe
Les diffrents modes de rejets
Terme source
Rejet
Evaporation de nappe
Dispersion du nuage
Transport / Diffusion
Conditions mtorologiques
Conditions environnementales
Modlisation de la dispersion
Effets toxiques
Notion de seuil toxique
Modlisation
45/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Modlisation de la dispersion: 3 approches
Modle gaussien
dispersion pilote uniquement par la mto
validit : 100 m < d < 10 km
Exemples : ALOHA (EPA/NOAA, USA), EPIcode
(Homann Associates, Inc.), INPUFF (EPA).
Modle de type intgral
rsolution simplifie de la mcanique des fluides
validit : 20 m < d < 10 km
Exemples : EFFECTS (TNO, Pays-Bas), SAFER
(Safer Systems, USA), PHAST (DNV Software,
UK).
Modle CFD (= Computational Fluid Dynamics)
Equation de la mcanique des fluides moins simplifies
validit : 1 cm < d < 2 km
46/72
Exemples
Fluidyn
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Fluent,
Phoenics,
Star-CD,
Flacs,
Modle Gaussien (1)
La concentration est calcule sur laxe du panache
Une loi Gaussienne (statistique) permet den dduire la concentration dans tout
le panache
47/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Modle Gaussien (2)
Rejet instantan
2
2
2
M
x u.t
y
z H
C ( x, y , z , t )
. exp .
2 . exp .
3/ 2
2
2
2
.
2
.
2
.
x
y
z
(2. ) . . .
x y z
M : masse totale rejete (kg)
u : vitesse du vent (m/s)
H : hauteur d'mission (m)
x, y, z carts-types de la distribution gaussienne
48/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Modle Gaussien (3)
Rejet continu
Un rejet prolong est assimil une succession de bouffes lmentaires
Les effets sadditionnent aux diffrents points de lespace, chaque instant
Concentration en un point = somme des concentrations instantanes en ce
point
49/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Modle Gaussien (4)
Principe
Solution de lquation de dispersion est une loi de distribution gaussienne dans lespace
Applications
Gaz neutres et stable (qui ne ragit pas chimiquement)
Situation durgence
Valables pour des distances comprises entre 100 m et 10 km
Hypothses simplificatrices
Relation simple entre la diffusion turbulente et le gradient moyen de concentration
Nglige la diffusion molculaire
Champs de vent uniforme (en vitesse et direction) dans le temps et l espace
Turbulence homogne et isotrope
Pas dobstacle ni de relief
Exemples de logiciels
Aloha - Cameo (EPA) : logiciel pour situation durgence
ADMS : logiciel de nouvelle gnration (tient compte de la topographie)
50/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Modle intgral
Principe
Rsolution simplifie des quations de la mcanique des fluides
Approche paramtrique prudente
Applications
Tout type de gaz (lger, lourd, neutre)
Limitations
Mto constante, ni reliefs, ni obstacles
Inadapt pour des rejets faible vitesse, vent faible car nglige la diffusion molculaire
Exemple de logiciels
EFFECTS (TNO, Pays-Bas)
Phast (DNV Software, UK) le plus courant en France
SAFER (SAFER System, USA)
51/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Phast v6.51
Modle CFD
Principe
Rsolution des quations de la mcanique des fluides
Utilisation de modles de turbulence
Applications / Hypothses
Prise en compte des fluctuations turbulentes, des obstacles et du relief
FDS modelling
Simulation denvironnements complexes
Si les enjeux ncessitent daffiner lestimation
Limitations
Nombreuses donnes dentre et paramtres numriques
Temps de calcul important
Prcision dpendant de la mthode de rsolution, du maillage utilis, du choix des modles...
Exemples
CFX, Fluent (Ansys) - Gnraliste
PHOENICS (Cham) - Gnraliste
Mercure (EdF) - Dispersion atmosphrique
FLACS (GexCon) - Dispersion et explosion
52/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Modle CFD: application dun modle (1)
Cas : Installation frigorifique - Rejet dammoniac dans la salle des machines
Objectif : Calcul des concentrations au niveau des maisons
53/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Modle CFD: application dun modle (2)
Utilisation du code Saturne
Maisons
Tunnel de liaison
Maillage
Btiment "Chambre
de conglation"
Domaine
Btiment
"Consommables"
Route Nationale n3"
54/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Modle CFD: application dun modle (3)
Concentrations en NH3 - t=160 s - Mto F3
Coupe horizontale
mi-hauteur des btiments
Coupe verticale
Axe de la chemine
55/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Modle CFD: cadre rglementaire
Les deux conditions pour lesquelles lutilisation de la CFD est accepte par
ladministration sont:
Scnarios prsentant un milieu compos de nombreux obstacles naturels ou
anthropiques
Scnarios aboutissant des distances deffets irrversibles (SEI) de plusieurs
milliers de mtres avec des modles intgraux ou gaussiens
(Des distances deffets irrversibles (SEI) infrieures 2000 m seront particulirement examines)
56/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Exercice (suite)
La rupture guillotine dun bras de dchargement de wagon est envisage. Les 60 tonnes de propylne
ou propne (C3H6) prsentes dans le wagon sont relches travers le bras de dchargement. Aucun
systme de scurit ne fonctionnant, la dure de fuite est indtermine.
1m
Question 1: En utilisant le logiciel ALOHA, modliser le nuage toxique hauteur dhomme 1000,
1500, 2000 ppm suite la dispersion du propylne (seule la contribution amont est tudie)
Donnes:
Dbit de fuite en amont: 34,2 kg
Dure de fuite: 30 min
Conditions mtorologiques: B/2, Temps clair, Temprature de lair: 15 C, Humidit: 70%, direction du vent: Ouest
Rugosit du sol: 3 cm
Question 2: Donner la concentration de propylne 200 m de la source hauteur dhomme dans laxe
du rejet
57/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Organisation de la journe
Les diffrents modes de rejets
Terme source
Rejet
Evaporation de nappe
Dispersion du nuage
Transport / Diffusion
Conditions mtorologiques
Conditions environnementales
Modlisation de la dispersion
Effets toxiques
Notion de seuil toxique
Modlisation
58/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Notion de seuil deffets toxiques (1)
Dfinition
Concentration dun gaz toxique dans lair ambiant partir de laquelle des
effets sur la sant des personnes exposes sont observs
En France, plusieurs seuils rglementaires
SEI: Seuil dEffets Irrversibles
SEL (CL 1%): Seuils des premiers Effets Ltaux (dangers graves)
SELS (CL 5%) : Seuils des Effets Ltaux Significatifs (dangers trs graves)
Dtermination dun seuil deffets toxiques associ un produit
A partir de diffrents rsultats dessais mens sur des rats / souris
Pour diffrentes dures dexposition et diffrents niveaux de concentration
Extrapolation de ces rsultats ltre humain, par des experts toxicologues
Dtermination des couples de donnes (dure dexposition / niveau de
concentration) capables de gnrer
des effets chez lhomme
59/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Notion de seuil deffets toxiques (2)
Exemple: SEI pour du Dichlore (Cl2) (Rf: fiche INERIS)
Dure d'exposition (min)
10
20
30
60
Concentration d'exposition (ppm)
110
41
30
25
19
Donnes accessible sur le site de lINERIS pour:
60 produits
Les 3 seuils deffets toxiques rglementaires (SEI, SEL, SELS)
Pour des dure dexposition comprises gnralement entre 1 et 60 min
60/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Relation de Haber (1)
Relation dose / effets
I = Cn(t) * dt
(Relation de Haber)
I : Intensit de leffet
C : Concentration
Hypothse
Dure dexposition gale dure dabsorption
n : Paramtre associ un produit et des effets toxiques recherchs
Dtermination du paramtre n
n dpend du produit et des effets recherchs
n peut tre calcul partir du tableau prcdent
t : Dure dexposition
61/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Relation de Haber (2)
Dtermination du paramtre n
Log (C)
Le produit toxique suit la
loi de Haber
Log (C1)
Droite de rgression linaire (a = -1/n)
Log (t1)
Log (t)
62/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Relation de Haber (3)
Dtermination du paramtre n
Log (C)
Log (C1)
Le produit toxique ne suit
pas la loi de Haber
Droite de rgression linaire
(a = -1/n2)
Rupture de pente
Droite de rgression linaire (a = -1/n1)
Log (t1)
Log (t)
63/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Liste des seuils rglementaires
Rappel des seuils toxiques rglementaires franais
SEI: Seuil dEffets Irrversibles
SEL (CL 1%): Seuils des premiers Effets Ltaux (dangers graves)
SELS (CL 5%) : Seuils des Effets Ltaux Significatifs (dangers trs graves)
Dpendent de la nature du produit et de la dure dexposition
A noter quen pratique, les seuils peuvent varier dune personne lautre suivant
les hypothses
Bb
Sportif
Personne ge
En labsence de seuils rglementaire franais, possibilit dutiliser dautres seuils
Guide pratique de choix des valeurs seuils de toxicit aigue en cas dabsence de
valeurs franaises (sur le site de lINERIS) 2008
Seuils amricains
64/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Exemple de seuils issus de la rglementation franaise
SEI 10 min
SEL 10 min
SELS 10 min
(ppm)
(ppm)
(ppm)
CO
2 600
7 000
ND
CO2
25 000
37 000
42 000
HCl
240
1 300
2 149
HF
600
1 123
1 705
HCN
62
110
174
NO2
60
100
118
SO2
129
1 148
1 358
Gaz toxique
Toxicit dun mlange despces (ex: fumes dincendie)
1
Seuileq
i 1
Ci
Seuili
Ci
= Gas i concentration (-)
Seuili
= Gas i toxic threshold (-)
Seuileq
= Equivalent threshold (-)
65/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Exemple de seuils issus de la rglementation amricaine (1)
IDLH 1987
ppm /
30 min
IDLH
1994
ppm /
30 min
CO2
50 000
CO
AEGL (ppm / 30 min)
ERPG (ppm)
40 000
1 500
1 200
*ND
150
200
350
600
NO
100
100
80
ND
ND
NO2
50
20
0,5
15
15
25
SO2
100
100
0,2
0,75
0,3
32
HF
30
30
34
20
62
HBr
50
30
43
250
HCl
100
50
1,8
43
20
210
HCN
50
50
2,5
10
10
21
Acrylonitrile
500
85
10
35
Ammoniac
500
300
30
220
25
150
1600
Formaldhyde
30
20
0,9
14
10
70
Phnol
250
250
19
29
10
50
ND
Benzne
3 000
500
73
1 100
50
150
5 600
Styrne
5 000
700
20
160
50
250
1 900
Tolune
2 000
500
200
570
50
300
4 200
H2S
300
100
0,6
32
0,1
30
59
Acide formique
30
30
50
CS2
500
500
200
50
600
0,03
0,18
0,1
0,5
2,5
10 000
2 000
45
340
10
200
1 100
Acroline
Actaldhyde
66/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
AEGL : Acute Exposure Guideline Level
ERPG : Emergency Response Planification Guideline
IDLH : Immediately Dangerous to Life or Health
Exemple de seuils issus de la rglementation amricaine (2)
AEGL : Acute Exposure Guideline Level
AEGL-1 is the airborne concentration, expressed as parts per million or milligrams per
cubic meter (ppm or mg/m3) of a substance above which it is predicted that the
general population, including susceptible individuals, could experience notable
discomfort, irritation, or certain asymptomatic nonsensory effects. However, the effects
are not disabling and are transient and reversible upon cessation of exposure.
AEGL-2 is the airborne concentration (expressed as ppm or mg/m3) of a substance
above which it is predicted that the general population, including susceptible
individuals, could experience irreversible or other serious, long-lasting adverse health
effects or an impaired ability to escape.
AEGL-3 is the airborne concentration (expressed as ppm or mg/m3) of a substance
above which it is predicted that the general population, including susceptible
individuals, could experience life-threatening health effects or death.
67/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
SEI
SEL
Exemple de seuils issus de la rglementation amricaine (3)
The Emergency Response Planning Guideline (ERPG) values are intended to
provide estimates of concentration ranges where one reasonably might anticipate
observing adverse effects as described in the definitions for ERPG-1, ERPG-2, and
ERPG-3 as a consequence of exposure to the specific substance.
The ERPG-1 is the maximum airborne concentration below which it is believed that
nearly all individuals could be exposed for up to 1 hr without experiencing other
than mild transient adverse health effects or perceiving a clearly defined,
objectionable odor.
The ERPG-2 is the maximum airborne concentration below which it is believed that
nearly all individuals could be exposed for up to 1 hr without experiencing or
developing irreversible or other serious health effects or symptoms which could
impair an individual's ability to take protective action.
The ERPG-3 is the maximum airborne concentration below which it is believed that
nearly all individuals could be exposed for up to 1 hour without experiencing or
developing life-threatening health effects.
68/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Seuil dEffets
Rversibles (60 min)
SEI (60 min)
SEL (60 min)
Exemple de seuils issus de la rglementation amricaine (4)
IDLH is an acronym for Immediately Dangerous to Life or Health, and is defined
by the US National Institute for Occupational Safety and Health (NIOSH) as
exposure to airborne contaminants that is "likely to cause death or immediate or
delayed permanent adverse health effects or prevent escape from such an
environment." Examples include smoke or other poisonous gases at sufficiently
high concentrations.
SEL
69/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Dmarches adoptes par la rglementation franaise
Dmarche antrieure
Limitation de dure de fuite 30 min
Etude des effets au sol uniquement
Nouvelle dmarche
Dure de fuite physiquement possible (jusqu 60 min)
Modlisation du panache toxique en hauteur
Distances deffets toxiques en hauteur
Maitriser lurbanisation: btiments grande hauteur aux abords des sites industriels
FILM
Phast v6.51
70/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Exercice
Question 1: Une fuite de Cl2 se produit dans une canalisation. Les diffrentes iso-concentration observe dans le nuage sont prsentes sur le
schma suivant:
200 ppm
100 ppm
50 ppm
20 ppm
1m
0
20
40
Distance (m)
Tracer directement sur le schma les distances deffets 1 m de hauteur au Seuil dEffet Irrversible (SEI) pour une dure de fuite de:
a)1 min
b)30 min
Question 2: En utilisant le logiciel ALOHA, calculer les distances deffets irrversibles (SEI) hauteur dhomme lies un rejet de CO dune
dure de:
a)10 min
b)30 min
Donnes:
Dbit de fuite de 20 kg/s
Conditions mtorologiques: D/5, couvert, temprature de lair: 5C, taux dhumidit: 50%, direction du vent: Est
Rugosit du sol: 20 cm
Hauteur de rejet: 1 m
SEI (CO 30 min) = 1 500 ppm
Question 3: Mme question avec une vitesse de vent de 10 m/s. Conclusion
71/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
Synthse sur la dispersion
Effets tudis
Toxiques (uniquement sur lhomme)
Seuils deffets issus de la rglementation franaise
SEI
SEL (CL1%)
SELS (CL 5%)
dpendent de la nature du produit et de la dure dexposition
Phnomnes susceptibles de gnrer un nuage de gaz toxiques
Brche dans une canalisation ou un rservoir
Rupture de canalisation
Evaporation de nappe
Ruine de rservoir
72/72
Ce document est la proprit de lINERIS. Il ne peut pas tre reproduit ou exploit sans son accord - 2012
