3-aprin.ppt [mode de compatibilité] - cedre
TRANSCRIPT
![Page 1: 3-aprin.ppt [Mode de compatibilité] - Cedre](https://reader030.vdocument.in/reader030/viewer/2022011723/61d2dd7655e15644f271126a/html5/thumbnails/1.jpg)
Paris,�27�mars�2013
Les�Journées�d’information�du�Cedre�
La�modélisation de�la�dispersion�des�produits chimiques
Laurent�Aprin
Institut�des�Sciences�des�Risques,�LGEI,�Ecole�des�Mines�d’Ales
![Page 2: 3-aprin.ppt [Mode de compatibilité] - Cedre](https://reader030.vdocument.in/reader030/viewer/2022011723/61d2dd7655e15644f271126a/html5/thumbnails/2.jpg)
Plan
1. Introduction• Contexte• Description phénoménologique
2. Les mécanismes de dispersion en milieu aquatique
3. Les stratégies de modélisation• Rejet de produit en surface• Rejet de produit en profondeur
4. Exemple de modélisation
2
![Page 3: 3-aprin.ppt [Mode de compatibilité] - Cedre](https://reader030.vdocument.in/reader030/viewer/2022011723/61d2dd7655e15644f271126a/html5/thumbnails/3.jpg)
MarineTraffic.com 3
Contexte
![Page 4: 3-aprin.ppt [Mode de compatibilité] - Cedre](https://reader030.vdocument.in/reader030/viewer/2022011723/61d2dd7655e15644f271126a/html5/thumbnails/4.jpg)
Accidents�maritimes significatifs impliquant des�Substances�Nocives�et�Potentiellement�Dangereuses�(SNDP)
Hazardand�
Noxious Substances(HNS)
Contexte
Source: European Maritime Safety Agency (EMSA)4
![Page 5: 3-aprin.ppt [Mode de compatibilité] - Cedre](https://reader030.vdocument.in/reader030/viewer/2022011723/61d2dd7655e15644f271126a/html5/thumbnails/5.jpg)
Evolution�des�accidents�maritimes impliquant des�SNPD
� Augmentation�constante du�transport�des�produits chimiques par�voies maritimes
Number of accidents involving HNS
Augmentation�de�la�probabilité�d’accidents�ayant�des�conséquences�graves�
5
Contexte
� Multiplication�des�échanges�par�3.5�en�20�ans
� Grande diversité des produits chimiques
� Comportements�physico�chimiques�différents�en�milieu�marin�� Code�SEBC
![Page 6: 3-aprin.ppt [Mode de compatibilité] - Cedre](https://reader030.vdocument.in/reader030/viewer/2022011723/61d2dd7655e15644f271126a/html5/thumbnails/6.jpg)
Grande�variété�des�produits�transportés
Différentes�stratégies�d’intervention�
Caractéristiques�du�polluant
� Pétroles� Produits chimiques� Produits agroalimentaires
Caractéristiques�des�zones�polluées
� Mer, océan, estuaires…� Parc marins, zones protégées� Zone d’aquacultures, de pêches ou de tourismes 6
Contexte
![Page 7: 3-aprin.ppt [Mode de compatibilité] - Cedre](https://reader030.vdocument.in/reader030/viewer/2022011723/61d2dd7655e15644f271126a/html5/thumbnails/7.jpg)
Exemple�d’accident
Ievoli�Sun�(2000)� Quantités�transportées
� Styrene :4000�T� Methyl Ethyl Ketone (MEK)�:1000�T� Iso�Propylique Alcohol (IPA)�:�1000�T
� Quantités�rejetées:�� Styrene :400�T�� Methyl Ethyl Ketone :100�T� Iso�Propylique Alcohol :�1000�T
7
Contexte
Entrée�d’eau�dans�le�compartiment�du�propulseur�d’étrave
Source: French Navy
![Page 8: 3-aprin.ppt [Mode de compatibilité] - Cedre](https://reader030.vdocument.in/reader030/viewer/2022011723/61d2dd7655e15644f271126a/html5/thumbnails/8.jpg)
Exemple�d’accident
Ievoli�Sun�(2000)
8
Contexte
SubstancesTransport�maritime
• Styrène• MEC• IPA
Phase S d Vp[%] [kg/m3] [kPa]
LLL
0,0326,379
906805786
0,66710,334,1
FEDED
� Formation�d’un�nuage�toxique�?� Risque�de�feu�ou�d’explosion�?
![Page 9: 3-aprin.ppt [Mode de compatibilité] - Cedre](https://reader030.vdocument.in/reader030/viewer/2022011723/61d2dd7655e15644f271126a/html5/thumbnails/9.jpg)
� Conditions�non�représentatives� Présence�possible�de�différents�produits
� Modification�de�la�solubilisation� Réaction�chimique
� Aucune�prise�en�compte�de�la�dynamique� Solubilisation� Compétitions�des�phénomènes
� Produit�pur� Petite�quantité� 20°C� Pression�atmosphérique� Eau�douce
Conditions�de�la�classification Conditions�réelles
� Pureté�non�garantie� Plusieurs�tonnes� 4°C� 10�bar� Eau�de�merLimites
9
Description�phénoménologique
![Page 10: 3-aprin.ppt [Mode de compatibilité] - Cedre](https://reader030.vdocument.in/reader030/viewer/2022011723/61d2dd7655e15644f271126a/html5/thumbnails/10.jpg)
Evaporant�(ex�:�Styrene)Coulant�(ex:�Toluene Diisocyanate)�
Flottant�(ex�:Dodecylbenzene) Soluble�(ex�:�butanol)
Exemple�de�comportements
10
Description�phénoménologique
![Page 11: 3-aprin.ppt [Mode de compatibilité] - Cedre](https://reader030.vdocument.in/reader030/viewer/2022011723/61d2dd7655e15644f271126a/html5/thumbnails/11.jpg)
11
Quantité
évaporée
Quantité
solubilisée
Etalement
Les�mécanismes�de�dispersion
Rejets de�polluants en�surface�1. Détermination�du�terme�source�(quantité�rejetée�ou�débit)�>�Bernoulli2. Bilan�de�matière�en�fonction�du�temps
Quelle�surface�de�produitQuelle�quantité�évaporée?Quelle�quantité�solubilisée�?Combien�de�temps�?
![Page 12: 3-aprin.ppt [Mode de compatibilité] - Cedre](https://reader030.vdocument.in/reader030/viewer/2022011723/61d2dd7655e15644f271126a/html5/thumbnails/12.jpg)
12N surfaces Sspillet Qevap et QsolSurface S
Modélisation
Stratégie de�ModélisationDivision�de�la�nappe�en�éléments�unitaires�(spillet)�pour�représenter�la�surface�totale�de�la�nappe
1. Calcul�de�l’évolution�de�chaque�spillet au�cours�du�temps�(Schéma�de�discrétisation�de�type�Runge�Kutta)�• Taux�d’étalement�de�chaque�spillet (Nihoul,�1983)• Quantité�évaporée�(Mackay�et�Matsugu,�1973)• Quantité�solubilisé�(Hayduk�et�Laudie,�1974)
2. Bilan�de�matière�pour�chaque�spillet
![Page 13: 3-aprin.ppt [Mode de compatibilité] - Cedre](https://reader030.vdocument.in/reader030/viewer/2022011723/61d2dd7655e15644f271126a/html5/thumbnails/13.jpg)
13
Modélisation
Stratégie de�ModélisationDivision�de�la�nappe�en�éléments�unitaires�(spillet)�pour�représenter�la�surface�totale�de�la�nappe
1. Calcul�de�l’évolution�de�chaque�spillet au�cours�du�temps�(Schéma�de�discrétisation�de�type�Runge�Kutta)�• Taux�d’étalement�de�chaque�spillet (Nihoul,�1983)• Quantité�évaporée�(Mackay�et�Matsugu,�1973)• Quantité�solubilisé�(Hayduk�et�Laudie,�1974)
2. Bilan�de�matière�pour�chaque�spillet
3. Couplage�avec�les�prévisions�hydrodynamiques1. Déplacement�des�spillets à�chaque�pas�de�temps
![Page 14: 3-aprin.ppt [Mode de compatibilité] - Cedre](https://reader030.vdocument.in/reader030/viewer/2022011723/61d2dd7655e15644f271126a/html5/thumbnails/14.jpg)
Rejet de�polluants en�profondeur
1. Détermination�du�terme�source�(quantité�rejetée�ou�débit)�>�Bernoulli2. Mode�de�dispersion�du�produit3. Hydrodynamique�de�remontée�(vitesse)4. Transfert�de�matière�(solubilisation)�
Apparition�du�produit�en�surfaceQuel�volume�en�surface�?Quand�?Combien�de�temps�?
Fraction
soluble
Fraction
flottante
Fraction
évaporée
14
Les�mécanismes�de�dispersion
![Page 15: 3-aprin.ppt [Mode de compatibilité] - Cedre](https://reader030.vdocument.in/reader030/viewer/2022011723/61d2dd7655e15644f271126a/html5/thumbnails/15.jpg)
Modélisation du�débit à�la�brèche
1. Deux�comportements�hydrodynamiques
2. Modélisation�par�un�bilan�énergétique�local�(Bernoulli)
15
Sortie�de�liquide
Entrée�d’eau
2�orifices
1�orifice
15
Modélisation
![Page 16: 3-aprin.ppt [Mode de compatibilité] - Cedre](https://reader030.vdocument.in/reader030/viewer/2022011723/61d2dd7655e15644f271126a/html5/thumbnails/16.jpg)
Modélisation du�débit à�la�brèche
16
Ecoulement�diphasique�à�l’orifice�de�sortie
Débit�Q�constant,Section�� linéaireDébit�Q�constant,Section�� linéaire
Ecoulement�monophasique�à�l’orifice�de�sortie
Modélisation
Débit�Q�linéaire,Section�� constante
Débit�Q�linéaire,Section�� constante
![Page 17: 3-aprin.ppt [Mode de compatibilité] - Cedre](https://reader030.vdocument.in/reader030/viewer/2022011723/61d2dd7655e15644f271126a/html5/thumbnails/17.jpg)
17
Bas
Haut
0 15 25 85 100 140Temps [s]
20 cm20 cm
2,2
m
Modélisation
Modélisation de�la�fragmentation�du�jet
Evaluation�de�la�distribution�granulométrique
![Page 18: 3-aprin.ppt [Mode de compatibilité] - Cedre](https://reader030.vdocument.in/reader030/viewer/2022011723/61d2dd7655e15644f271126a/html5/thumbnails/18.jpg)
Modélisation de�la�fragmentation�du�jet
� 2�types�de�distributions�granulométriques� Diamètre�maximal�des�gouttes�de�22�mm
18
Q < 6 L/minQ > 6 L/min
18
Modélisation
![Page 19: 3-aprin.ppt [Mode de compatibilité] - Cedre](https://reader030.vdocument.in/reader030/viewer/2022011723/61d2dd7655e15644f271126a/html5/thumbnails/19.jpg)
Modélisation de�la�fragmentation�du�jet
� Représentation�des�distributions�granulométrique�par�des�log�normale1. Rejet�monophasique�:�distribution�mono�modale2. Rejet�diphasique�:�distribution�bi�modale
19Rejet�monophasique Rejet�diphasique
Modélisation
![Page 20: 3-aprin.ppt [Mode de compatibilité] - Cedre](https://reader030.vdocument.in/reader030/viewer/2022011723/61d2dd7655e15644f271126a/html5/thumbnails/20.jpg)
Modélisation de�l’hydrodynamique�des�gouttes
1. Modélisation�de�la�vitesse�de�remontée�des�gouttes�avec�la�corrélation�de�Clift (1978)
2. Couplage�avec�les�prévisions�hydrodynamiques�pour�prendre�en�compte�le�déplacement�du�panache�de�bulles
20
Modélisation
CDOG�model,�Clarkson�Univeristy
![Page 21: 3-aprin.ppt [Mode de compatibilité] - Cedre](https://reader030.vdocument.in/reader030/viewer/2022011723/61d2dd7655e15644f271126a/html5/thumbnails/21.jpg)
Modélisation d’une�fuite�de�butanol�sur�un�chimiquier
1. Fuite�d’une�cuve�de�600�m32. Brèche�sur�un�piquage�de�10cm�de�diamètre
Source: BEAMerCEPPOL�:�Centre�d'Expertises�Pratiques�de�Lutte�Antipollution
21
Exemple
![Page 22: 3-aprin.ppt [Mode de compatibilité] - Cedre](https://reader030.vdocument.in/reader030/viewer/2022011723/61d2dd7655e15644f271126a/html5/thumbnails/22.jpg)
Modélisation d’une�fuite�de�butanol�sur�un�chimiquier
1. Modélisation�du�débit�de�fuite�:�vidange�en�1h442. Dispersion�dans�la�colonne�d’eau
1. Nombre�de�gouttes�par�classe�de�diamètre�pour�1�m3�de�produit�libéré2. Répartition�du�volume�par�classe�pour�1�m3�de�produit�libéré
22
Exemple
0 2 4 6 8 10 12 14 16 18 200
1
2
3
4
5
6
7x 105
Diamètre [mm]
Effe
ctif
0 2 4 6 8 10 12 14 16 18 200
0.01
0.02
0.03
0.04
0.05
0.06
0.07
0.08
Diamètre [mm]
Vol
ume
[m3]
![Page 23: 3-aprin.ppt [Mode de compatibilité] - Cedre](https://reader030.vdocument.in/reader030/viewer/2022011723/61d2dd7655e15644f271126a/html5/thumbnails/23.jpg)
Modélisation d’une�fuite�de�butanol�sur�un�chimiquier
1. Solubilisation�des�différentes�classes�de�gouttes• Evolution�du�volume�en�fonction�du�temps�(�t =�15 s)
23
Exemple
0 2 4 6 8 10 12 14 16 18 200
0.01
0.02
0.03
0.04
0.05
0.06
0.07
0.08
Diamètre [mm]
Vol
ume
[m3]
t�=�0
t�=�2�min�30�s
![Page 24: 3-aprin.ppt [Mode de compatibilité] - Cedre](https://reader030.vdocument.in/reader030/viewer/2022011723/61d2dd7655e15644f271126a/html5/thumbnails/24.jpg)
Modélisation�d’une�fuite�de�butanol�sur�un�chimiquier
Evolution�du�volume�en�surface�sur�le�volume�initial�en�fonction�de�la�hauteur
Après 10 m : 96% du volume total est solubiliséAprès 20 m : le butanol est totalement solubilisé
0 5 10 15 20 250
10
20
30
40
50
60
70
80
X: 10Y: 3.24
Hauteur parcourue h [m]
Vol
ume
à h
/ Vol
ume
libér
é [%
]
24
Exemple
![Page 25: 3-aprin.ppt [Mode de compatibilité] - Cedre](https://reader030.vdocument.in/reader030/viewer/2022011723/61d2dd7655e15644f271126a/html5/thumbnails/25.jpg)
Modélisation de�la�dispersion�des�produits�chimiques�en�mer
Prise�en�compte�d’autres�mécanismesSédimentationBio�accumulationVolatilisation
Prise�en�compte�des�paramètres�environnementauxTempérature�de�l’eauSalinitéPression�hydrostatique
25
Le�futur
![Page 26: 3-aprin.ppt [Mode de compatibilité] - Cedre](https://reader030.vdocument.in/reader030/viewer/2022011723/61d2dd7655e15644f271126a/html5/thumbnails/26.jpg)
Merci�de�votre�attention
26