Analyse de sensibilité paramétrique d'un outil de modélisation des conséquences de scénarios d'accidents. Application à la dispersion atmosphérique de rejets avec le logiciel Phast.

Christian Delvosalle

Description

Niveau: Supérieur, Doctorat, Bac+8

• mémoire

THÈSE En vue de l'obtention du DOCTORAT DE L'UNIVERSITÉ DE TOULOUSE Délivré par l'Institut National Polytechnique de Toulouse Discipline ou spécialité : Génie des Procédés et de l'Environnement JURY Christian DELVOSALLE, Professeur, Faculté Polytechnique de Mons (Rapporteur) Pascal FLOQUET, Professeur, INPT/ENSIACET Nadine GABAS, Professeur, INPT/ENSIACET (Directrice de thèse) André LAURENT, Professeur Emérite, INPL/ENSIC (Rapporteur) Eric MARSDEN, Docteur en Informatique, ICSI Yvon MOUILLEAU, Ingénieur, Société Technip Ecole doctorale : MEGeP Unité de recherche : Laboratoire de Génie Chimique (UMR 5503) Directeur(s) de Thèse : Nadine GABAS Présentée et soutenue par Nishant PANDYA Le 1er décembre 2009 Titre : Analyse de sensibilité paramétrique d'un outil de modélisation des conséquences de scénarios d'accidents. Application à la dispersion atmosphérique de rejets avec le logiciel Phast.

• outil de modélisation des conséquences de scénarios d'accidents


• dispersion atmosphérique de rejets avec le logiciel phast

Sujets

• Mémoire
• Laurent
• Deleuze
• Floquet
• Marsden
• Bertrand

Informations

• Publié par: profil-feym-2012
• Publié le: 01 décembre 2009
• Nombre de lectures: 76
• Langue: Français
• Poids de l'ouvrage: 3 Mo

Extrait

THÈSE


En vue de l'obtention du

DOCTORAT DE L’UNIVERSITÉ DE TOULOUSE DOCTORAT DE L’UNIVERSITÉ DE TOULOUSE

Délivré par l'Institut National Polytechnique de Toulouse
Discipline ou spécialité : Génie des Procédés et de l'Environnement


Présentée et soutenue par Nishant PANDYA
erLe 1 décembre 2009

Titre : Analyse de sensibilité paramétrique d’un outil de modélisation
des conséquences de scénarios d’accidents.
Application à la dispersion atmosphérique de rejets avec le logiciel Phast.

JURY
Christian DELVOSALLE, Professeur, Faculté Polytechnique de Mons (Rapporteur)
Pascal FLOQUET, Professeur, INPT/ENSIACET
Nadine GABAS, Professeur, INPT/ENSIACET (Directrice de thèse)
André LAURENT, Professeur Emérite, INPL/ENSIC (Rapporteur)
Eric MARSDEN, Docteur en Informatique, ICSI
Yvon MOUILLEAU, Ingénieur, Société Technip

Ecole doctorale : MEGeP
Unité de recherche : Laboratoire de Génie Chimique (UMR 5503)
Directeur(s) de Thèse : Nadine GABAS

















à Yujie ( 宇 杰).
i ii

Exprimer les remerciements reste une tache très difficile pour moi, que ce soit à
l’oral ou à l’écrit. En attendant de nouveaux développement dans la communication par
télépathie, je vais essayer d’exprimer ma gratitude en quelques mots à toutes les
personnes m’ayant -- directement ou indirectement, scientifiquement et/ou
moralement -- aidé pendant mon travail de thèse. On m’a appris que les remerciements
les plus longs sont parfois les moins sincères. Je serais donc bref.

Je remercie ma directrice de thèse, Professeur Nadine GABAS, pour ses conseils
avisés, sa grande disponibilité et ses encouragements tout au long de mes travaux.
J’adresse mes remerciements à Monsieur René DELEUZE, Directeur Général de
l’Institut pour une Culture de Sécurité Industrielle (ICSI) pour m’avoir accueilli au sein de
l’Institut, et Monsieur Gilles VACHER, Secrétaire Général de l’ICSI, pour son dynamisme et
ses encouragements.
Je tiens à également remercier Professeur Joël BERTRAND, Directeur du Laboratoire
de Génie Chimique, de m’avoir accueilli au sein du laboratoire.
Je tiens à remercier les partenaires industriels de ce projet pour avoir financé ma
thèse et leurs correspondants pour leur soutien et la richesse des échanges scientifiques
pendant nos réunions. J’aimerais également témoigner de mon appréciation pour
l’implication de Monsieur Henk WITLOX, de la société DNV Software, qui n’a pas hésité à
partager son expertise sur le logiciel Phast.
Un grand merci à Professeur André LAURENT et Professeur Christian DELVOSALLE
d’avoir accepté la charge de rapporter sur ce mémoire et pour leurs commentaires
constructifs. Je remercie également Professeur Pascal FLOQUET de me faire l’honneur de
présider mon jury de thèse, ainsi que pour son implication pendant la première phase de
thèse. Je remercie Monsieur Yvon MOUILLEAU d’avoir accepté de participer au jury.
Je ne peux pas oublier Monsieur Eric MARSDEN, ICSI, qui m’a beaucoup aidé pendant
ces trois années, particulièrement par son expertise en informatique et par sa bonne
humeur ; sans lui la forme des résultats serait moins attrayante.
Maintenant j’en profite pour remercier mes collègues de l’ICSI, Pauline, Dominique,
Zina, Caroline, Gilles, …, pour leur sympathie et l’ambiance amicale.
J’ai une pensée pour mes parents, mon frère, mes proches ainsi que ma belle
famille.
Enfin, toutes mes pensés vont à Yujie, ma femme, pour son dévouement éternel,
ses sacrifices et son attention sans limite depuis le début de mon séjour en France. Elle est
toujours à mes côtés pendant les moments difficiles de ma vie. Je lui dédie cette thèse.




iii
ivRESUME
Dans le cadre de la réalisation des études de danger, la réglementation française
prévoit le calcul de périmètres de sécurité autour des établissements industriels à risques.
Ces valeurs conditionnent la mise en œuvre de mesures importantes en matière de
prévention des accidents et de protection des cibles (personnes, environnement et
installations). La libération accidentelle de gaz toxique et/ou inflammable constitue un
scénario d’accident majeur pouvant conduire à des distances d’effets importantes. Il en
résulte la propagation de nuages dans l’atmosphère selon des phénomènes extrêmement
complexes régis par de nombreuses variables. L’étude de la dispersion atmosphérique de
ces nuages fait appel à des logiciels de simulation comme le logiciel commercial Phast
(Process Hazard Analysis Software Tool, DNV) de type intégral, très souvent utilisé dans
les analyses de risques d’installations industrielles. Les utilisateurs de cet outil ont
constaté que les résultats de simulation sont dépendants d’un très grand nombre de
paramètres de modélisation ajustables. Or, même si l’utilisateur est expert, il ne dispose
pas toujours d’une bonne connaissance des valeurs pertinentes pour ces paramètres en
fonction des scénarios simulés, ni de leur impact sur les résultats du modèle.
L’objectif est de mener une étude approfondie sur l’analyse de la sensibilité
paramétrique du logiciel Phast. L’enjeu principal est l’acquisition de connaissances sur la
simulation de scénarios accidentels et de leurs impacts, permettant une meilleure
utilisation de cet outil de calcul.
La technique mise au point a consisté à coupler Phast et l’outil d’analyse de
sensibilité SimLab, ce qui permet d’exécuter automatiquement un grand nombre de
simulations en faisant varier l’ensemble des paramètres du modèle de façon simultanée.
La méthode d’analyse de sensibilité globale choisie, EFAST (Extended Fourier Amplitude
Sensitivity Test), est basée sur l’analyse de la variance des sorties du modèle pour le calcul
des indices de sensibilité.
Nous avons étudié des scénarios de rejet continus pour six produits différents
(monoxyde d’azote, ammoniac, chlore, azote, n-hexane et fluorure d’hydrogène),
sélectionnés pour couvrir une large gamme de caractéristiques physiques (gaz lourd, gaz
passif, gaz léger, gaz réactif avec l’atmosphère) et de conditions de stockage (gaz sous
pression, gaz liquéfié sous pression, liquide saturé, …) ainsi que pour leur importance
industrielle dans les analyses des conséquences d’accidents.
Nous avons tout d’abord fait l’analyse de deux termes sources de Phast (fuite sur
un réservoir et rupture totale d’une canalisation courte) en considérant les sorties telles
que le débit de rejet, la vitesse et la température finales du rejet. En ce qui concerne
l’analyse du modèle de dispersion de Phast appelé UDM (Unified Dispersion Model), la
démarche a été séparée en deux étapes. La première étape de « screening » a pour but de
comparer l’influence de l’ensemble des paramètres (y compris ceux liés aux conditions
météorologiques et au terme source). Elle nous a permis de vérifier que les paramètres les
plus influents sont ceux liés au terme source et aux conditions météorologiques, et non
les paramètres « internes », liés à la modélisation. La deuxième étape est l’analyse de
l’influence globale des paramètres de modélisation, autres que les paramètres météo et
du terme source, sur une plage large de valeurs (souvent les bornes minimales et
maximales permises par le logiciel de modélisation). Pour chaque produit, nous avons
décomposé les scénarios de base en sous-scénarios correspondant à des conditions de
rejet différentes (combinaisons de débit, hauteur de rejet, angle de rejet et conditions
météorologiques) dans le but de faciliter l’interprétation des résultats. Cette granularité
vfine dans les scénarios étudiés nous a permis d’interpréter les sensibilités des paramètres
sans la confusion qui résulterait du mélange de phénomènes physiques très différents.
Une étude complémentaire a consisté à effectuer une analyse de sensibilité locale
de ces paramètres autour de leur valeur par défaut. Ceci nous perm