Afdækning af muligheder for etablering af standardværktøjer og/eller -kriterier til vurdering af sundheds- og miljørisici i forbindelse med større uheld (gasudslip) på risikovirksomheder 5 Modeller til gasspredningsberegning
Fysiske modeller er den basis på hvilken spredningsberegningsværktøjer er bygget. Værktøjerne kan ikke forventes at fungere godt hvis der er problemer med modelleringen. Der er dog mange problemer selv om man har studeret gasspredning siden 1930’ern. Hovedgrunden hertil er, at luftens turbulens er meget kompleks. Man har dog opnået et niveau hvor man kan sige at modellerne er ”brugbare”. Der er stadigvæk nogle områder ved beregning af udslipskilderne som er stadigvæk ikke undersøgt, tilsyneladende på grund af manglende interesse eller forglemmelser. Status for modeludviklingen beskrives i dette kapitel. Usikkerhed i beregninger er vigtig fordi resultaterne anvendes til planlægning og har konsekvenser både for anvendelse på kommuneområdeplanlægning og på sikkerhed. Enkeltmodeller til gasspredning har en usikkerhed som er, for de bedste modeller, indenfor en faktor 30 % for simple tilfælde, men mange af modellerne der anvendes i dag har en afvigelse fra eksperiment resultater på en faktor 2. Dette kan betyde betydelig øgning i størrelsen af sikkerhedszoner, for eksempel fra 300 m til 600 m. i praktiske tilfælde. En sådan usikkerhed har stor betydning for planlægning. Der er derfor her lagt vægt på modellernes nøjagtighed. Usikkerhed i resultaterne afhænger ikke kun af spredningsberegninger. Uslipsberegninger der anvendes til at fodre spredningsberegninger med input data har også usikkerheder. Usikkerhederne lægges sammen. Det er derfor nødvendig at sikre så god kvalitet i modellering på alle stadier i udslippet og spredningen, hvis man vil opnå et rimeligt samlet resultat. 5.1 ModellerBeregning af gasspredning kræver at der er fysiske modeller til hver fase af udslippet og spredningen. Selve spredningen kan beskrives med ca. tolv modeller (jet og fane, for let og tung gas, fri, i bymiljø og indendørs), men der skal også være beregninger til frigivelse. Figur 5.1 viser de forskellige modeller og deres sammenhæng i fremstilling af rimelige nøjagtige og dækkende beregninger. I figur 5.1 anvendes grøn til at vise de steder hvor der er velafprøvede og validerede modeller som passer til de fleste praktiske tilfælde. Gul viser tilfælde hvor modelsættet er kun delvis udviklet. Rød viser beregninger hvor der er kun lidt eksperimentalinformation i dag. Bemærk at en grøn farve betyder ikke at en bestemt modeltype er totalt problemfri. Toksicitetsdata er for eksempel problemfri for ca. 100 af de vigtigste stoffer, men der er mange stoffer der anvendes i mindre målestok, hvor data mangler. Modeldetaljer beskrives i Bilag A og et oversigt over status ved udvikling af modeller gives i tabel 10.2. Figur 5.1 Modeller som anvendes til beregning af relevante udslipsscenarier i Danmark. Fremgangsmåden for beregninger vist i figur 5.1 er kun en del af den samlede riskikoanalyse- og risikovurderingsproces. En samlede risikoanalyse omfatter også: beregninger for brande og eksplosioner; beregninger af hyppigheder for udslip og pålidelighed af sikkerhedsforanstaltninger; plotning af risikokort og frekvens/konsekvens kurver; og sammenligning af resultaterne med akceptkriterier. 5.2 Usikkerheder i beregningerneAlle de modeller der anvendes til beregninger involverer usikkerhedsmomenter. Usikkerhederne stammer fra flere forskellige kilder:
Et problem med disse kilder til usikkerhed er, at de opstår i flere modeller, som skal anvendes en efter den anden. En usikkerhed på 50 % i modeller til udslip, til initialspredning i jet, og til gasfanespredning vil medføre en usikkerhed på ca. 150 % (se Bilag A), dvs., at der vil være en chance for, at resultatet er 150 % større end det er beregnet. Usikkerheder skal derfor bekæmpes i alle de forskellige delmodeller i en beregning. 5.3 Manglende scenarierDet er ønskeligt at et standardværktøj skal dække alle de scenarier, der vil kunne forekomme og give alvorlige konsekvenser. Det er dog usandsynligt at dette mål vil kunne opnås. En mere rimelig krav er måske at kunne beregne alle de uheldskonsekvenser der er opstået i praksis. De eksisterende programmer har betydelige svagheder, set fra dette synspunkt. Eksempler er modeller til spredning af gas fra sprøjtende væsker, eller fra væsker der løber fra tanke på grund af overfyldning. Sådanne mangler vises fra uheldet i Buncefield, England i 2006, hvor benzin løb fra toppen af en tank, fordampede, og hvor den resulterende dampfane var af meget større omfang end forventet på basis af pølfordampningsmodeller (ref. 30). Der er identificeret i alt 7 uheldssituationer i procesanlæg for hvilken der ingen publicerede modeller eksisterer. 5.4 Simple modeller og modelpakkerDe programmer, der anvendes i dag til risikoanalyse, er næsten udelukkende baseret på simple analytiske formler, dvs. formler der for det meste kan beregnes med en lommeregner, hvis man har tålmodighed nok. Nogle programmer implementerer kun enkelte modeller, for eksempel spredning af gas, men de programmer der avendes er oftest integrerede pakker, der inkluderer modeller for udslip, fordampning, spredning, brand og eksplosioner, og inkluderer databaser over stofegenskaber. Eksempler beskrives i det følgende kapitel. Simple modeller er ikke nødvendigvis dårligere end de mere avancerede og indviklede modeller der beskrives nedenfor. Den mest nøjagtige model til spredningsberegning i dag er Britter og MacQaid’s, ref. 32, der er en simpel formel tilpasning til eksperiment resultater. De mere indviklede modeller er tilpasset de sammen eksperimenter, som de simple. Fordelen med de mere avancerede modeller er, at de kan behandle situationer som f.eks. gasspredning omkring bygninger, hvor de simple modeller slet ikke passer. 5.5 Computational Fluid DynamicsDer er kommet metoder i de senere år der anvender avancerede numeriske metoder til at løse de fundamentale ligninger for gasspredning. Metoderne har hidtil været for krævende til almindelig anvendelse i sikkerhedsvurdering og især i risikoanalyse. Metoderne har typisk krævet timer eller dage til at udføre en enkelt beregning. Der er kommet nyere metoder i den senere tid som kræver meget mindre beregningstid, i visse tilfælde ned til 10 minutter per beregning. Der er dog ofte flere dages eller ugers indsats i opbygning af modellerne. Værktøj der anvender disse metoder beskrives i det følgende kapitel. Disse modeller anvendes i stort omfang til beregning af gasspredning på olie- og gasplatforme, hvor anlægsbudgetter berettiger omfattende beregninger. Disse modeller kan ikke være meget mere nøjagtige til beregning af simpelt tilfælde end de simplere modeller, idet alle skal justeres, således at de passer til de samme eksperimenter. CFD modellerne har den fordel, at de kan anvendes i langt flere situationer end simple modeller, og kan for eksempel behandle udslip indendørs og i tætte bebyggelser. De fleste CFD programmer kræver en del ekspertise til at opnå et korrekt resultat. (Det samme kan siges om simple modeller, men omfanget af den nødvendige ekspertise er større for CFD metoder). Det anbefales ofte, at man tager et kursus i programanvendelse før man anvender den. Nogle kommercielle programmer sælges kun med kursusafgift indbygget i prisen. Det vil være nødvendigt, for at få reproducerbar resultater ved anvendelse af disse programmer, at der bliver udarbejdet grundige vejledninger i, hvordan man opnår et pålideligt resultat. 5.6 Vigtigheden af modelundersøgelserDer er ingen naturlig teori, der kan forudsige spredningen af gasfaner, og som samtidig tillader praktiske beregninger med eksisterende computer. Alle de anvendte modeller baseres mere eller mindre på teori suppleret med observationer og eksperimenter. Derfor bør ingen model bør anvendes, med mindre dens resultater er blevet bekræftet med eksperimenter. Der er blevet foretaget mange undersøgelser af spredningsberegningsmodeller for at bestemme nøjagtigheden af beregningsresultater. Disse har for det meste beskæftiget sig med fanespredning og jet udslip, med nogle få sammenligninger af momentanudslip og af fordampning af kryogen væsker. Undersøgelser af selve gasspredning af tung gas er relativt godt undersøgt, med undersøgelser fra US EPA (ref. 33 ), og et undersøgelse af S.R. Hanna et al. (ref. 28). Nyere undersøgelser er SMEDIS, lavet af en samarbejdsgruppe for EU (resultater kun delvis publiceret (ref. 34). Dispersion modeller vurderes ofte i forhold til eksperimenter med parametrene ”geometric mean bias” og ”geometric variance”. ”Mean bias” måler graden med hvilken beregningerne systematisk undervurderer eller overvurderer resultater. ”Geometric variance” måler den gennemsnitlig afvigelse fra måleværdierne for simple tunggasspredning. CCPS (ref. 28) offentliggjorde resultater for 10 meget anvendte programmer vurderet mod 6 eksperimentserier (32 eksperimenter i alt, med 132 målepunkter). To af programmerne kunne udelukkes på grund af gennemgående dårlige resultater. Resten havde en ”mean bias” og ”geometric variance” af mindre end en faktor 2. Vi har gentaget disse vurderinger med sammenligninger med nogle nyere eksperimenter. Resultaterne gives i bilaget, sektion 9.9. Modeller for mere komplicerede situationer, som for eksempel når en gasfane rammer en væg (”impingement”), er mere sjældne, men der findes nyere eksperimentserier der kan anvendes (ref. 35 ). Der findes flere sæt eksperimenter om spredning i bymiljø, der blev lavet for at kunne vurdere frigivelse af giftstoffer ved terrorisme. Disse er relevant i det nuværende sammenhæng fordi de omfatter de samme fænomener som udslip ved uheld (ref. 36, 37). Der er også flere eksperimenter om let gas spredning, lavet primært for forureningsvurderinger, der kan anvendes til vurdering af uheldsudslip. Nogle eksperimenter er også blevet lavet til vurdering af tunggas udslip i industriel miljø, ved anvendelse af vindtuneller (f.eks. ref. 38). Der er mange nyere eksperimenter og vurderinger om tofase udslip, med henblik især at bestemme omfanget af ”rain out” dvs. nedfald ad dråber fra ammoniak og lignende udslip af fordråbede gasarter(se f.eks. ref. 39 til 44). Der er kun lidt data offentliggjort i nyere tid om mængden og hastigheden af tofaseudslip, men omfattende eksperimenter blev foretaget til afprøvning af de oprindelige tidlige publikationer (se f.eks. ref. 43) Eksperimenter og valideringer for udslip og af spredning af væske på jorden, og indendørs spredning, er sjældne. Det er ofte nødvendigt at søge oprindelige data, i visse tilfælde 60 år eller mere tilbage, for at finde eksperimental data. Data om fordampning blev indsamlet af MacKay og Matsugu (ref 47) og det er disse data der anvendes oftest. Tidligere data blev samlet af Sutton (ref 45) og Pasquill (ref.47). MacKay and Matsugu’s er stadigvæk den mest anvendte model men bør erstattes af en senere model af Kawamura og MacKay (ref.48). Modellerne passer enten til letflygtige væsker ((Kawamura og McKay) eller til sværfordampelige væsker(MacKay og Matsugu). Modellerne passer kun til fordampning under neutralt atmosfærisk forhold. Modellerne anvendes dog oftest ukritisk. Modellen af Brighton er tilpasset moderne forståelse af turbulens i atmosfæren, og passer meget bedre til eksperimentale resultater, men Brighton’s model anvendes sjældent. Der er nogle relevante områder hvor der næsten intet findes eller intet hvor der er ikke offentliggjorte måledata. Eksempler er væskespredning på jorden på realistisk industrielt jordforhold, udslip fra realistisk hulformer, og gasspredning i bygninger. Det er næsten meningsløst at tale om valideret modeller, når der er disse huller i det videnskabelig grundlag for beregninger. Nogle få eksperimenter er blevet fortaget for at ”lappe på hullerne” og disse afslører en relativ stor usikkerhed i de ”tommelfingerregler” som anvendes på nogle områder. Keller og Simmons har udført eksperimenter om pool spredning (ref. 50 til 52), med henblik på miljøbeskyttelse, men eksperimenter er igen næsten ukendt i risikoanalyse kredse. En af de største ”huller” er manglen på validerede modeller for fordampning fra plaskende og sprøjtende udslip. Rediphem som er et europæisk samarbejde om data spredning, med hjemsted på Risø, giver en omfattende samling af data om tung gas spredning (ref. 53). En oversigt og samling af referencer om brugbar eksperimentresultater gives i ref. 54. 5.7 Kriterier der kan stilles for spredningsberegningerDet er ikke overraskende, med de ovenstående usikkerheder, og med mange alternative modeller, at to forskellige analyser af det samme situation producerer forskellige resultater. Situationen er kendt fra mange ingeniørgrene. Der er, for eksempel, betydelige usikkerheder i beregning af styrken af en bygning eller en bro. Disse usikkerheder håndteres normalt ved anvendelse af standardberegningsmetoder og sikkerhedsfaktorer. Spørgsmålet er, om vi kan stille lignende krav til spredningsberegninger og de sikkerhedsanalyser der baseres på dem. Der er forskellige typer af krav der kan stilles til beregningerne:
5.8 A priori kontra tunede modellerModeller til udslipsberegning kan laves på basis af en grundlæggende teori. Beregninger af udslippet af væske fra et rundt hul kan for eksempel baseres på Bernoulli’s ligning, som følger direkte fra energibevarelse og fra løsning af flow langs strømlinier. Modellen baseres ikke på eksperimentale resultater. Modellen bekræftes igennem eksperimenter men anvender ingen eksperimentalt baserede konstanter eller tilpasninger i sin opbygning. Sådanne modeller hedder ”a priori” modeller. Der er kun få modeller der anvendes i spredningsberegninger som er rene a priori modeller. Men mange er i sin overordnede form baseret på et teoretisk grundlag, med delmodeller som er baseret på eksperimenter der beskriver enkelte fænomener. De fleste tunggas spredningsmodeller, for eksempel, er baseret på masse, energi og momentumbevarelse, men anvender resultater fra detaileksperimenter til bestemmelse af værdier af diffusionskonstanter, varmetransmissionskonstanter og konstanter for effektiviteten af omdannelse af potentiel til kinetisk energi. Modsætningen til ”a priori” modeller er rene empiriske modeller, der baseres på tilpasning af formler til eksperimentale resultater. Formlerne vælges frit for at opnå det bedste tilpasning til eksperimenterne uden nødvendigvis at tage hensyn til underliggende teori. Modellerne gøres rimeligt bredanvendelig ved at laver tilpasning på basis af dimensionsfrie variabler. En sådan model er Britter og McQuaid’s model for tunggas spredning. ”A priori” modeller betragtes som bedre end empiriske modeller fordi de forklarer mere og burde i princip være mere robust, idet de vil passe situationer der ligger udenfor området der er dækket af eksperimenter. Empiriske modeller vil altid viser sig at være de bedste hvis kvaliteten af modellerne vurderes alene udefra nøjagtigheden af overensstemmelse med eksperimentale resultater. Overensstemmelsen kan gøres perfekt, hvis modellen gøres tilstrækkelig kompliceret. Men en empirisk model vil måske vise sig at være unøjagtig hvis den sammenlignes med nye eksperimentale resultater. En mellemstadie imellem rene empiriske modeller og a priori modeller er ”tunede” modeller. Disse baseres på a priori modeller, men med en eller få konstanter for hvilken der er ingen teoretiske grundlag. Tilpasningen til eksperimentale resultater laves med at vælge værdier for disse konstanter. Muligheden for tuning gør valget af modeller mere vanskelig. Man bør ikke vælge modellen der passer eksperimenterne bedst hvis den model er baseret på rene empiri fordi man ved ikke hvornår der vil opstå situationer hvor modellen passer ikke. Der er til gengæld ingen grund til at foretrække en model med et godt teoretisk grundlag hvis de beregnede resultater passer dårligt til eksperimenterne. Et kompromis er oftest nødvendigt. Et rimeligt kriterium er at man bør vælge den a ”priori model” der passer eksperimentale resultater bedst såfremt den underlæggende teori er i orden og den ikke er væsentlig dårligere end den bedste empiriske model. Der er alligevel den beslutning at tage, om man bør tune modellerne. Ref. 25 beskriver vanskelighederne i at vælge den bedste model. 5.9 Fremgangsmåder til at opfylde kriterierneIngen af de nuværende vejledninger til spredningsberegninger af toksiske gasarter der anvendes i dag opfylder kravene fra den forrige sektion. Der er dog flere fremgangsmåder hvormed det vil være muligt at forbedre situationen:
5.10 Modeller for konsekvensreduktionKonsekvensreduktion er et krav i de bekendtgørelser der anvendes til styring af risici ved industrielle anlæg. Konsekvensreduktion er også ønskelig, fordi der er flere situationer hvor store konsekvenser kan give betydelige kommuneplanlægningsproblemer. Eksempler på konsekvensreducerende foranstaltninger der overvejet eller anvendt i praksis i Danmark er:
Det er ønskelig at have anerkendte modeller til dette formål, både til støtte for projektering af risikoreduktion, og for at kunne vurdere indvirkning af disse effekter på planlægning. Der findes modeller og eksperimentale data til alle disse formål. (f.eks. refs. 62 til 68) 5.11 Spredningsberegninger og risikoanalyserKvantitative risikoanalyser anvendes i dag til kommuneplanlægning, tillige med de kvalitative. Kvantitativ risikoanalyse foretages med beregning af mange scenarier og scenario varianter, ofte op til 90 scenarier per beholder. Det er specielt for kvantitative metoder, at der skal studeres et bredt spektrum af scenarier (i princip samtlige scenarier) og at sandsynligheder for disse scenarier skal beregnes. Kvalitative risikoanalyser baseret oftest på anvendelse af hazard og operabilitetsmetoden (HAZOP), der undersøger alle de forstyrrelser der teoretisk forstyrrelser, der kan opstå på et procesanlæg. Fremgangsmåderne er komplementære. Risikoacceptkriterier fra ”Kvantitative og kvalitative kriterier for risikoaccept”, Miljøprojekt 112 (ref. 69), har igennem ca. 15 år gjort det muligt at vurdere risikoen ved ”Seveso virksomheder” på en objektiv og konsistent måde under den forudsætning, at der blev taget rimeligt hensyn til usikkerhederne i beregningerne. Spredningsberegninger anvendes i dag i Danmark til både kvalitative og kvantitative risikoanalyser. Det vurderes at ca. halvdelen af de sikkerhedsrapporter der laves i Danmark i dag udføres med kvantitative metoder som et element af rapporten (forfatterens vurdering), og at anvendelsen af kvantitative metoder er så godt som 100 % for stort anlæg og anlæg med potentiale for stor konsekvenser. Risikoberegninger af denne type kræver data for hyppigheden af udslip og for pålidelighed af sikkerhedsforanstaltninger. Usikkerheden i risikoberegninger stammer til dels fra spredningsberegninger, og til dels fra hyppigheds data. I dag foregår der et ihærdigt arbejde på EU planen for at udvikle standard hyppighedsdata til anvendelse i sikkerhedsrapporter (ref. 70). Det forventes at dette arbejde vil medføre en betydelig forbedring i både reproducerbarhed og nøjagtighed af risikoanalyserne indenfor det kommende år. 5.12 OpsummeringBeregninger af spredning af toksisk gas involverer flere modeller til de forskellige stadier i et udslipsscenarier typisk mellem tre og fem modeller. Der er betydelige usikkerheder på op til 50 % for flere af modellerne og større for nogle få modeller til nogle scenarier. En gennemgang af tabel A.2 viser, at der kan foretages spredningsberegninger for mange af scenarierne med et minimum af usikkerhed, ca. en faktor 2 overalt. Dette kræver omhyggelig tilrettelæggelse af beregningsvejledninger, således at usikkerhederne der opstår fra datavalg og scenarievalg minimeres. Modeller der baseres direkte på eksperimentdata kan opnå en bedre nøjagtighed end en faktor 2. Der er større områder hvor usikkerheden i beregningerne er større og nogle områder hvor der ikke findes anerkendte metoder. Spredningsberegninger fortages fortrinsvis med simple modeller i dag og der er områder hvor disse modeller slet ikke passer. Eksempler er udslip i tæt bebyggelser og indendørs. Resultatet er, at beregninger ofte ikke vil kunne anvendes til vurdering af risikoen for medarbejderne. Der er muligheder for forbedring af situationen, for eksempel ved validering af modeller til nærfeltsberegninger og indendørsberegninger og ved en større anvendelse af CFD metoder. Der vil ligeledes være muligheder for reduktion af usikkerheder ved at udvikle bedre vejledninger til beregninger, der minimerer unødvendig variationer i beregningsresultater. 5.13 Krav til spredningsberegningsmodellerDet er muligt at udlede funktionelle krav til standardberegningsmetoder på basis af de ovenstående modelvurderinger. Andre krav der stammer fra anvendelsesbehov, specielt ved beredskabsberegninger, udledes i de senere kapitler.
|