Source: http://docplayer.fi/2171116-Pc-cosyman-kayttoonottoja-soveltaminen.html
Timestamp: 2017-10-19 14:34:15+00:00
Document Index: 15143716

Matched Legal Cases: ['kko ', 'kko ', 'kko ', 'kko ', 'kko ', 'kko ', 'kko ', 'kko ', 'kko ', 'kko ']

PC COSYMAn käyttöönottoja soveltaminen - PDF
PC COSYMAn käyttöönottoja soveltaminen
Download "PC COSYMAn käyttöönottoja soveltaminen"
1 STUK-YTO-TR 96 FI PC COSYMAn käyttöönottoja soveltaminen Jukka Rossi JOULUKUU 1995 SÄTEILYTURVAKESKUS Strölsäkerhetscentralen Finnish Cenfre for Radiation and // QJUKUlhi Nuclear Safety
2 STUK-YTO-TR 96 JOULUKUU 1995 PC COSYMAn käyttöönottoja soveltaminen Jukka Rossi VTT Energia Ydinenergia Tutkimuksen johto Säteilyturvakeskuksessa Risto Paltemaa, Arto Isolankila Tutkimus on tehty Säteilyturvakeskuksen tilauksesta. SÄTEILYTURVAKESKUS PL 14, HELSINKI (90)
3 ISBN ISSN Painatuskeskus Oy Helsinki 1995
4 STUK-YTO-TR 96 SÄTEILYTURVAKESKUS ROSSI, Jukka (VTTEnergia). PC COSYMAn käyttöönottoja soveltaminen. STUK-YTO-TR 96. Helsinki s. + liitteet 11 s. ISBN ISSN Avainsanat: PSA3, reaktorionnettomuuksien seurausvaikutusten arviointi, radionuklidien dispersio TIIVISTELMÄ COSYMA on EU-maiden tutkimuslaitosten, pääasiassa NRPB:n ja KfK:n, 1980-luvulla kehittämä reaktorionnettomuuksien ympäristövaikutusten laskentamalli. Vuonna 1993 ohjelmasta valmistui mikrotietokoneella toimiva yksinkertaistettu versio, PC COSYMA. Tämä versio on tarkoitettu käytettäväksi tutkimuskäytössä sekä determinististen että todennäköisyyspohjaisten ympäristövaikutusten arvioimiseen. Deterministinen laskenta sisältää vain yhden leviämistilanteen ja todennäköisyyspohjainen laskenta kaikki mahdolliset leviämistilanteet ja niiden esiintymistodennäköisyydet sääaineistoon perustuen. Ohjelmassa on merkkipohjainen käyttöliittymä sekä lähtötietojen syöttämistä että tulosten katselua varten, joten tuloksista ei suoraan saada graafista esitystä näytölle. Laitoskohtaisten ympäristö- ja säätietojen valmistaminen tapahtuu kuitenkin tiedostojen välityksellä. Ohjelman laskenta-ajat ovat suhteellisen lyhyitä. Käyttöä tukee selkeä opaskirja. Kaikilla lähtötiedoilla on oletusarvot, joita käyttäjä voi muuttaa sovellustensa ja tarpeidensa mukaan. Sovelluslaskuina käsiteltiin eräitä pitoisuus-, annos-, terveysvaikutus-, saastuneen alueen pinta-alasekä vastatoimenpidelaskuja. Osassa laskuja tuloksia verrattiin ARANO-ohjelmalla saatuihin tuloksiin, jotta voitiin varmistua siitä, että laskettiin todella sitä, mitä oli ajateltu laskettavan. Vertailu osoitti tulosten hyvää yhteensopivuutta näiden kahden mallin välillä pitoisuus- ja annoslaskennassa. PC COSYMA-ohjelman käyttöönotto ja soveltaminen osoittivat, että ohjelma sopii käytettäväksi laitoskohtaisissa ympäristövaikutusarvioinneissa. Käyttäjän olisi tunnistettava ne ohjelmiston monipuolisista piirteistä johtuvat ominaisuudet, jotka ratkaisevasti vaikuttavat laskentatuloksiin.
5 SÄTEILYTURVAKESKUS STUK-YTO-TR 96 ROSSI, Jukka (VTT Energy). Installation and implementation of PC COSYMA. STUK-YTO-TR 96. Helsinki pp. + Appendices 11 pp. ISBN ISSN Keywords: PSA3, reactor accident consequence assessment, dispersion of radionuclides ABSTRACT COSYMA is a probabilistic reactor accident consequence assessment model developed in the 1980's by the research institutes of the EC countries, primarily by NRPB and KfK. In 1993 a simplified version of the program was released to be implemented on a PC. This version is purposed to be used in research usage for deterministic as well as probabilistic offsite consequence assessments. Deterministic approach includes only a single set of atmospheric conditions whereas probabilistic approach takes into account all possible atmospheric conditions and their respective frequency of occurence on the basis of the weather data. The program has an alfanumeric menu driven user interface for both input and output, so no graphical results are available on the display. The site specific gridded environmental and weather data has to be prepared through input files. The execution times of the program are relatively short. Operation is supported by a comprehensive user's guide. All the input parameters have default values, which can be modified according to applications and needs of the user. In this study the PC COSYMA program was applied to calculate as the end points nuclide concentrations, doses received, numbers of health effects, areas affected by countermeasures and other results related to countermeasure strategies. Some results are compared with the results obtained by the ARANO program developed by VTT Energy to quarantee proper interpretation of the application. This comparison indicates rather good accordance of the results in the case of concentrations and doses received as end points. The installation and implementation of PC COSYMA has proved out that the program is a flexible tool for site specific offsite consequence assessments. The user should be rather careful regarding some properties of the program which are due to the extensive features of the system and which decisively affect the results obtained.
6 STUK-YTO-TR 96 SÄTEILYTURVAKESKUS SISÄLLYSLUETTELO TIIVISTELMÄ Sivu ABSTRACT ESIPUHE 6 1 JOHDANTO 7 2 OHJELMISTON MALLINNUKSEN YLEISKUVAUS Päästön määrittely Dispersio ilmassa ja depositio Annoslaskenta Ulkoisen annoksen laskenta Sisäisen annoksen laskenta Terveysvaikutukset Vastatoimenpiteet Pelkkä sisällesuojautuminen Sisällesuojautuminen ja sitä seuraava evakuointi Joditabletti Väestön pitempiaikainen poissiirto Elintarvikkeiden käyttökielto 16 3 OHJELMISTON KÄYTTÖKELPOISUUS 17 4 OPERATIIVISET TYÖVAIHEET Ohj elman käyttöönotto j a lähtötiedot Sovelluslaskut Yksittäisen leviämistilanteen pitoisuus- ja annoslaskut TVO:n poistokaasujärjestelmän vaurioon liittyvät annoslaskut Vakavaan reaktorionnettomuuteen liittyvän päästön raja-arvo Vakavaan onnettomuuteen liittyvät vastatoimenpidelaskelmat 26 5 YHTEENVETO 31 VIITTEET 32 MUITA COSYMA-RAPORTTEJA 32 LIITE 1 Katsaus ohjelmassa olevaan laitoskohtaiseen ympäristötietoon ja COSYMA-ohjelman vaatimukset 33 LIITE 2 PC COSYMAn asennus ja käyttöönotto 38 LIITE 3 Experiences with PC COSYMA 1.0 (full version) 40
7 SÄTEILYTURVAKESKUS STUK-YTO-TR 96 ESIPUHE Työ on tehty Säteilyturvakeskuksen tilauksesta Valtion teknillisen tutkimuskeskuksen (VTT Energia) Ydinenergian tutkimusalueella. Yhteyshenkilöinä ovat toimineet ylitarkastaja Risto Paltemaa ja tarkastaja Arto Isolankila (STUK) sekä tutkija Jukka Rossi VTT:stä.
8 STUK-YTO-TR 96 SÄTEILYTURVAKESKUS 1 JOHDANTO Tämän projektin tarkoituksena oli ottaa käyttöön Valtion teknillisen tutkimuskeskuksen (VTT) Ydinenergian tutkimusalueella EY: n MARIA-tutkimusohjelman lopputuotteena syntyneen COSYMA-ohjelman (Code System from MARIA) IM mikrotietokoneen DOS-käyttöjärjestelmässä toimiva ohjelmaversio (PC COSY- MA121) ja avustaa Säteilyturvakeskusta (STUK) tämän ohjelman käyttöönotossa. Lisäksi ohjelmalla tuli suorittaa eräitä sovelluslaskuja, joista vaativimmat koskisivat vastatoimenpiteiden tehokkuutta suomalaisessa ydinlaitosympäristössä. Samalla oli tarkoitus kerätä, päivittää ja saattaa COSYMA:n vaatimaan muotoon tarpeellinen ympäristöä koskeva lähtötieto sekä sääaineisto. Vertailulaskujen tekoon VTT:n ARA- NO-ohjelmalla varauduttiin. Projektin tarkoitus oli arvioida COSYMAm käyttöä ja soveltuvuutta suomalaisessa ympäristössä sekä parantaa valmiutta osallistua EY:n projekteihin. COSYMA on tutkimuskäyttöön tarkoitettu ydinvoimalaitosonnettomuuksien ympäristön seurausvaikutusten todennäköisyysperustainen arviointimalli. COSYMA-mallin pääkehittäjä on ollut Kernforschungszentrum Karlsruhen (KfK) ja National Radiological Protection Boardin (NRPB) tutkijaryhmä, jota on avustanut joukko eurooppalaisia yhteistyökumppaneita. COSY- MA-ohjelmaa on kehitetty toistaiseksi siten, että sen kahdesta laiteriippuvasta versiosta on keskuskoneympäristössä toimivaa täydellistä ohjelmaversiota kehitetty ensisijaisesti ja vain harkittu osa sen ominaisuuksista on siirretty PCversioon, joka on pysynyt aina suppeampana laskentamahdollisuuksiltaan. PC-versiosta on saatavilla ns. laajempi ja suppeampi versio, joista viimemainitussa on vähäisin parametrien vaihtelumahdollisuus.
9 SÄTEILYTURVAKESKUS STUK-YTO-TR 96 2 OHJELMISTON MALLINNUKSEN YLEISKUVAUS Seuraavassa on esitelty pääpiirteissään ne mallinnuksen ominaisuudet, jotka vaikuttavat oleellisesti ohjelman käyttökelpoisuutta ja sopivuutta arvioitaessa. Mallin kuvaus perustuu vahvasti viitteeseen 2, joka antaa hyvän yleiskuvan PC COSYMA.sta. 2.1 Päästön määrittely Päästön suuruus voidaan määritellä kahdella vaihtoehtoisella tavalla: käyttämällä lyhyttä nuklidiluetteloa tai lähdetermiohjelmaa, jossa on käytettävissä 145 nuklidia. Päästössä voi olla mukana korkeintaan 50 nuklidia. Päästön suuruus määrätään antamalla inventaari ja päästöosuus. Tässä yhteydessä voi määritellä kaikki olennaiset päästöön liittyvät parametrit: päästön alkuviive ja kesto, päästökorkeus, vapautuva lämpöenergia ja reaktorirakennuksen mitat. Päästöjaksoja voi olla enintään kolme, niiden ei tarvitse olla peräkkäisiä ja niissä on radioaktiivinen hajoaminen aina mukana. Kukin jakso on tunnin pituinen. Jodi voidaan jakaa alkuainemuotoiseen, orgaaniseen ja aerosolimuotoon. Koska lähdetermiohjelmassa on mahdollista antaa yli 50 nuklidia, on käyttäjän joko annettava annosreiteittäin katkaisuraja siten, että nuklidien määrä on korkeintaan 50 tai muuten supistettava nuklidien määrää. Tämän rajoituksen tarpeellisuus ilmenee vasta ohjelmaa ajettaessa virheilmoituksena, minkä seurauksena voi joutua työlääseen iterointiin nuklidimäärän rajoittamiseksi. Käyttäjällä saattaa myös olla tarvetta laskea sellaisia päästöjä, joissa tietty aktiivisuusmäärä nuklidia vapautuu, mikä ei suoraan ole mahdollista tällä ohjelmalla, koska radioaktiivinen hajoaminen on aina mukana. Lähdetermin voi tallettaa tiedostoksi, jota voi siten jatkossa käyttää sujuvasti muissa ajoissa. 2.2 Dispersio ilmassa ja depositio Ilmassa leviämisen laskenta jakaantuu olennaisesti kolmeen osaan: säätilanteen poiminta, dispersio ja depositio. Parametrien asettelumahdollisuudet riippuvat paljon siitä lasketaanko determinististä vai probabilistista leviämistä. Yksittäistä leviämistilannetta laskettaessa voidaan joko lukea alkamishetken perusteella leviämistilanne käyttäjän antamasta meteorologisesta tiedostosta tai antaa leviämiseen liittyvät parametrit suoraan valikosta. Jälkimmäisessä tapauksessa sama leviämistilanne säilyy koko leviämisen ajan, mutta ensimmäisessä se voi muuttua datatiedostossa olevan mukaisesti. Yksittäistä leviämistilannetta voidaan muokata monipuolisesti sopivilla parametriasetuksilla. Alkuvaiheen pilvennousema sekä rakennusten vaikutus käsitellään Briggsin ja Hannan periaatteilla /3,4/. Alkuperäisen ohjelman viidestä dispersiomallista on tässä mukana vain yksi. Dispersio ilmassa käsitellään segmentoidulla gaussisella mallilla (MUSEMET), joka perustuu säätiedon käsittelyyn tunnin aika-askelella siten, että tuulensuunta ja -nopeus, stabiilius sekä sateen intensiteetti voivat vaihdella leviämisen aikana tällä taajuudella. Leviämissuunta voidaan asettaa kolmella vaihtoehtoisella tavalla: jokainen jakso kulkee päästön alkaessa vallinneeseen suuntaan, jokainen jakso kulkee niiden alkuhetkellä vallinneeseen suuntaan, jokainen jakso kulkee tunneittain päästölähteessä vallitsevan suunnan mukaisesti. Meteorologiseen datatiedostoon voi sisältyä yhden tai kahden vuoden sääaineisto. Todennäköisyyspohjaisessa laskennassa voidaan käyttää joko syklistä tai aineiston jaotteluun (stratified sampling) perustuvaa poimintamenettelyä. Edellisessä yhden näytteen todennäköisyys on l/n, 8
10 STUK-YTO-TR 96 SÄTEILYTURVAKESKUS kun n on näytteiden kokonaismäärä. Syklisessä poiminnassa eivät välttämättä tule esille harvinaisimmat leviämistilanteet. Siinä myöskin aikavälin valinnalla voi olla merkitystä lopputulokseen (esim. 24 h). Nyt on mahdollista valita 144 näytettä, mikä tarkoittaa 61 tunnin sykliä vuoden aineistosta. Toinen mahdollisuus on jakaa sääaineisto ryhmiin siten, että kukin ryhmään sisältyvä näyte tuottaa suunnilleen samansuuruisen tuloksen. Mikäli ryhmien määrittely on onnistunut, saadaan aineisto jaoteltua siten, että näytteet edustavat laskentatulosten komplementaarisen kumulatiivisen todennäköisyysjakauman (CCDF = Complementary Cumulative Distribution Function) koko arvoaluetta. Kustakin ryhmästä poimitaan edustaja tai useampi varsinaiseen laskentaan ja todennäköisyys saadaan ryhmässä olevien näytteiden kokonaismäärän ja valitun näytemäärän avulla. Näytteet voidaan poimia satunnaisesti ryhmän sisältä. Tässä versiossa kuitenkin vain yksi näyte kustakin ryhmästä voidaan poimia. Leviämisen aikana lasketaan radioaktiivinen hajoaminen ja tytärnuklidien muodostuminen. Deposition käsittely tapahtuu neljässä eri ryhmässä: aerosolit, orgaaniset aineet sekä jodin elementaarinen ja orgaaninen (esim. metyylijodi) muoto. Käyttäjä voi antaa depositionopeuden ja pesukertoimeen liittyvät parametrit ryhmittäin. Pesukerroin lasketaan sateen intensiteettiin perustuen. 2.3 Annoslaskenta Kun mitään vastatoimenpiteitä ei käytetä, ovat COSYMA:n altistus- ja annosreitit: ulkoinen gammasäteily pilvestä, pilvessä olevien radioaktiivisten aineiden hengittämisestä aiheutuva altistus, ulkoinen gammasäteily laskeumasta, ulkoinen beetasäteily iholle tai vaatteisiin deposoituneista aineista, laskeumasta resuspension välityksellä hengityselimistöön tulevat radionuklidit, laskeuma-alueella tuotettujen ravintoaineiden nauttimisesta aiheutuva annos. Annokset lasketaan ulkona olevalle keskimääräiselle aikuiselle 20-vuotiaalle henkilölle eikä mitään esimerkiksi lapsia koskevia tuloksia lasketa. Annokset lasketaan elinkohtaisina ja halutulle ajanjaksolle lyhyen ajan annosten osalta (1 365 vuorokautta), mutta pitkän ajan annoksissa altistusaikaa ei voi antaa vaan se on aina 50 vuotta. Altistuksen alkuhetkeä ei voi antaa. Käyttäjä voi valita valikoissa lyhyen ja pitkän ajan annosreitit. Altistusreittikohtaisesti elimet ja siten myös terveysvaikutukset vaihtelevat. Varsinaisten yksilö- ja väestöannosten ohella näitä arvoja käytetään jatkossa terveysvaikutusten riskin laskennassa. Ravintoaineannoksia ei käytetä arvioitaessa varhaisvaikutuksia eikä niiden aiheuttamaa yksilöannosta lasketa lainkaan tässä versiossa. Annos lasketaan napakoordinaatiston pisteessä ja on siis segmentin keskiakselin arvo. Väestöannos saadaan kertomalla annos napakoordinaatiston segmentin väestömäärällä Ulkoisen annoksen laskenta Pilvestä tulevan ulkoisen gamma-annoksen laskenta perustuu puoliäärettömän avaruuden aproksimaatioon ja maanpintapitoisuuden arvoon. Korjauskertoimella otetaan huomioon pilven todellinen geometria ja altistujan etäisyys keski - akselilta. Laskeumasta aiheutuva ulkoinen annos lasketaan nuklidikohtaisen kirjaston avulla. Se sisältää annoksen elintä ja yksikkölaskeumaa kohden eripituisille ajanjaksoille ottaen huomioon tytärnuklidien aiheuttama vaikutus annokseen laskeuman jälkeen. Ohjelmassa on kaksi eri mallia nuklidien merkittävyyden mukaan jaoteltuna. Tärkeiksi katsotuille nuklideille (Ru-103, Ru-106, 1-131, Te-132, Cs-134, Cs-137 ja Ba- 140) kehitetty malli kuvaa taajaan asuttua aluetta ja sisältää tarkennetun laskennan mainitulle nuklideille. Se sisältää laskeuman käsittelyn erityyppisille pinnoille (katto, seinä, nurmi, maaperä) sekä siirtymistä ja migraatiota niissä. Mallin kertoimet on laskettu Iso-Britannian kaupunkiympäristöjen väestömäärillä painotetuiksi arvoiksi. Muiden nuklidien osalta käytetään yksinkertaisempaa avoimen alueen olosuhteita kuvaavaa mallia, jossa on mukana vain maaperän migraatio.
11 SÄTEILYTURVAKESKUS STUK-YTO-TR 96 Iholle ja vaatteille deposoituneiden nuklidien aiheuttama beeta-annos iholle on verrannollinen maahan kuivadeposoituneeseen määrään samassa paikassa. Käyttäjä voi antaa verrannollisuuskertoimen, altistusajan sekä tehollisen puoliintumisajan iholla. Vaatteiden pinnalla oleva kontaminaatio aiheuttaa paljon pienemmän annoksen kuin paljaalle iholle deposoitunut materiaali Sisäisen annoksen laskenta Pilvessä olevan pitoisuuden aiheuttama hengitysannos lasketaan integroidun ilmapitoisuuden, hengitysnopeuden ja annos/yksikköpitoisuus - tekijöiden tulona. Hengitysnopeuden tulisi olla aikuisen arvo, koska annosmuunnoskertoimet ovat aikuiselle. Resuspension aiheuttama annos lasketaan samoista tekijöistä, mutta pitoisuus arvioidaan laskeumasta riippuvalla kaavalla, jonka kertoimet käyttäjä voi antaa. Maaperän pitoisuudessa otetaan huomioon radioaktiivinen hajoaminen sekä migraatio syvemmälle maahan. Koska ravinnosta ei lasketa yksilöannoksia, perustuu altistusreitin väestöannoksen laskenta laskeuman, ravintoainepitoisuuden yksikkölaskeumaa kohden, tuotetun ravintoainemäärän ja ravintoaineessa olevan yksikköpitoisuuden aiheuttaman annoksen tuloon. Altistus alkaa päästön alkaessa. Eräät ravintoaineet voidaan jakaa tuoreeseen, jonka kulutusviive on lyhyt sekä jalostettuun osaan, jonka kulutusviive on pitempi. Nämä viiveet käyttäjä voi antaa. Ravintoainemallin parametrien arvioidaan olevan sopivia keski- ja pohjoiseurooppalaisessa viljelyssä. Ravintoaineet ovat: lehmän maito, naudan liha, naudan maksa, lampaanliha, lampaan maksa, vihannekset, juurekset ja vilja. 2.4 Terveysvaikutukset Terveysvaikutusten laskenta jaetaan varhais- ja myöhäisvaikutusten laskentaan. Varhaisvaikutusten (fatal ja non-fatal) laskenta perustuu NRC:n ja NRPB:n suosittelemiin riskimalleihin (hazard function)/2,5,6/. Mallin elimet ja niiden vaikutukset parametreineen (muotokerroin, kynnysarvo ja se annoksen arvo, jolla puolet altistuvasta väestöstä saa terveysvaikutuksen) käyvät ilmi taulukosta I. Malleihin sisältyy myös annosnopeusriippuvuuksia siten, että pienen annosnopeuden aiheuttama riski on myös alhaisempi kuin suuren annosnopeuden. Tämän laskemiseksi pitää annokset laskea aikaväleittäin. Luuytimelle ja keuhkoille on käytössä hieman tarkempi malli kuin muille elimille (taulukko II). Myöskin usean eri elimen annoksen aiheuttaman riskin yhteisvaikutus otetaan huomioon laskemalla yhdistetty riski. Myöhäis vaikutuksina käsitellään 11 syöpälajia sekä perinnöllisiä seurauksia. Luuytimen altistuksesta aiheutuva leukemia sekä luunpinnan altistuksesta aiheutuva syöpä lasketaan absoluuttisen riskin periaatteella, jolloin ei oleteta riippuvuutta normaalista taustalla vallitsevasta syöpäriskistä. Muiden syöpälajien riskikertoimet perustuvat suhteelliseen malliin, jossa otetaan huomioon väestön taustariski sekä ikäjakauma altistuksen aikana. Myöskin tässä lasketaan sekä kuolemaanjohtavat että ei-fataalit tapaukset lineaarisen kertoimen, jolla ei ole kynnysarvoa, avulla. Myöhäisvaikutusten laskenta on kokonaisuudessaan komplisoitu tehtävä, koska väestö koostuu eri ikäisistä yksilöistä ja ilmenevät syövät voivat esiintyä vasta vaihtelevan latenssiajan jälkeen ja lisäksi altistusta voi tapahtua pitkän ajan kuluessa, joka tulisi ottaa huomioon riskiä laskettaessa. Tämän takia tässä ohjelmaversiossa riski lasketaan yksinkertaisemmin annoksen ja riskikertoimen tulona. Käyttäjä voi valita joko GSF:n tai ICRP-60:n riskikertoimet, jotka on annettu taulukossa III. Näin napakoordinaatin pisteessä lasketun yksilöriskin avulla saadaan väestön riski kertomalla se väestömäärällä. Ravintoaineannosten riskiä laskettaessa oletetaan, että tuotanto ja kulutus tapahtuvat samassa paikassa. 2.5 Vastatoimenpiteet Vastatoimenpiteillä pyritään rajoittamaan joko akuuttia säteilyaltistusta onnettomuuden alkuvaiheessa ja pian sen jälkeen tai myöhempää altistusta laskeuman johdosta. Vastatoimenpiteiden mallinnus on pyritty tekemään mahdollisimman joustavaksi, jotta eri käyttäjät pystyisivät soveltamaan vapaasti haluamiaan toimenpiteitään ja kriteereitään. 10
12 STUK-YTO-TR 96 SÄTEILYTURVAKESKUS Taulukko I. PC COSYMA:ssa käytetyt oletusarvot terveydellisten varhaisvaikutusten mallinnuksessa /2/. dose in organ/tissue external internal irradiation irradiation effect slope parameter S median dose D50 (Sv) (1) upper end of exposure interval (days) threshold (Sv) mortality lung lung pulmonary syndrome 7.0 See table 2b 5.0 red bone marrow remainder (2 red bone marrow - haematopoietic syndrome gastrointestinal syndrome See table 2b skin - skin bums (5%) ovaries uterus pre- and neonatal death (3) lung lung morbidity lung function impairment 7.0 See table 2b 2.3 thyroid thyroid hypothyroidism external exposure iodine inhalation 10.0 skin - skin bums skin - cataracts ovaries uterus mental retardation (3) Notes (1) In accidents involving releases of a-emitters, the absorbed dose should be multiplied by a factor to account for the relative biological effectiveness of acute a- irradiation. (2) Since the dose conversion factors for the Gl-tract were not included In the GSF-data, the calculations are based on the dose to the remainder organs. (3) Risks predicted using these parameters have to be multiplied by 0.01 to be applicable to the general population, to allow for the fraction of that population who are pregnant at any one time Mahdolliset vastatoimenpiteet ovat: sisällesuojautuminen, evakuointi (suojaväistö), joditabletit, maaperän ja rakennusten puhdistus, väestön pitempiaikainen poissiirto (relokaatio) ja elintarvikkeiden käyttökielto. Näiden kumoaminen tapahtuu määritellyn kriteerin pohjalta. Sekä varhais- että myöhäisvaikutuksiin voi vaikuttaa paikallissuojautumisella, evakuoinnilla ja joditableteilla. Maaperän ja rakennusten puhdistus, väestön pitempiaikainen poissiirto ja elintarvikkeiden käyttökielto vaikuttavat vain myöhäisvaikutusten määrään. Aluksi lasketaan vastatoimenpiteen vaikutus annokseen ja terveysvaikutusten riskiin. Tämän pelkän annostulostuksen lisäksi on vastatoimenpiteiden vaikutusta ympäristössä havainnollistettu siten, että tulostetaan vastatoimenpiteen todennäköisyys etäisyyspisteittäin ja sektoreittain sekä henkilömäärä ja Taulukko II. PC COSYMA:n oletusarvot terveydellisille varhaisvaikutuksille keuhko-ja luuydinannoksen laskennassa /2/. Effect Mortality Pulmonary syndrome Haematopoietic syndrome Morbidity Lung function impainment Note D o GyVhour (1) For these effects the D50in any time period k given by D_ + D0/dose rate. alueen koko, joita vastatoimenpide koskettaa sekä elintarvikkeiden osalta hylätyn elintarvikkeen määrä. 11
13 SÄTEILYTURVAKESKUS STUK-YTO-TR 96 Taulukko III. PC COSYMA.n oletusarvot myöhäisvaikutusten riskikertoimille /2/. Organ Effect Fraction of cancers which are fatal Based on GSF Risk per Sv(1) Based on ICRP-60 Bone marrow Leukaemia E E-04 Bone surface Cancer E E-04 Breast Cancer E E-04 Lung Cancer E E-04 Stomach Cancer E E-02 Colon Cancer E E-04 Liver Cancer E E-03 Pancreas Cancer E E-03 (2) Thyroid Cancer ,77E E-04 Skin Cancer E E-04 Remainder Cancer E E-04 (3) Total 5.05E E-02 Gonads Hereditary 2.00E E-02 Notes (1) For leukaemia and cancer, these give the risk of a fatal effect. For hereditary effects, it gives the risk of severe effects in all future generations. (2) ICRP does not consider this organ. The value is obtained from the GSF value using a DDREF of 2.0. (3) This covers all organs not identified above, and is not the same as the ICRP remainder. Taulukko IV. PC COSYMA.ssa käytössä olevat vastatoimenpiteet, niihin liittyvät altistusreitit sekä niiden oletusarvoiset integrointiajat /2/. Action invoked Dose combination considered Default integration period Sheltering only (12) Sheltering and evacuation Relocation Stable iodine tablets Return from evacuation or relocation Food bans Notes (1 2) Cloud ydose Deposited y dose Inhalation dose Cloud ydose Deposited y dose Inhalation dose Cloud ydose Deposited y dose Inhalation dose Inhalation of resuspended material Automatic only Cloud ydose Deposited ydose Inhalation dose Inhalation of resuspended material Concentration in food n.a< 3» 7 days 50 years n.a (3) 7 days 50 years n.c' 4» 1 year n.c <4 > n.c (4 > n.c< 4) 1 year n.c (4) n.c 141 instantaneous (1) These actions may also be invoked automatically. (2) It ttiese actions are both invoked on the basis of dose, then the same set of organs and integration times must be used for both. (3) An integration period is not applicable for this route of exposure. The dose from the whole of the plume passage is considered. (4) This route of exposure can be considered in the criterion, but is not included by default 12
14 STUK-YTO-TR 96 SÄTEILYTURVAKESKUS Taulukko V. PC COSYMA.ssa käytössä olevat vastatoimenpiteet sekä niiden elinkohtaisiin annoksiin liittyvät oletusarvoiset toimenpiderajat /2/. Action invoked Sheltering (1) Sheltering and evacuation (1> Relocation Return from relocation or evacuation Organ considered Bone marrow Lung Gl Tract Thyroid Effective dose Bone marrow Lung Gl Tract Thyroid Effective dose Effective dose Effective dose Default dose intervention level Sv n.c < 2 > 0,05 n.c (2) 0,05 0,005 n.c (2) 0,5 n.c < 2 > 0,5 0,05 0,05 0,025 Notes (1) It both these actions are considered, the same set of organs must be used, but the intervention levels may be different. (2) This organ can be considered, but h Is not included by default. Taulukko VI. PC COSYMA.ssa käytössä olevat elintarvikkeiden vastatoimenpiteiden nuklidiryhmät ja niiden oletusarvot /2/. Nuclide groups Strontium Iodine a-emitters Caesium For milk Bq/I Intervention level (1) For other foods Bq/kg Notes (1) Different levels can be used for imposing and withdrawing food bans. The default values are the same hr both. The levels are those recommended by CEC (reference 42). Taulukossa IV on esitetty vastatoimenpiteet, niihin liittyvät annosreitit ja mahdollisten integrointiaikojen oletusajat. Taulukossa V on esitetty elinkohtaisiin annoksiin perustuvat kriteerit ja taulukossa VI elintarvikkeiden käyttökiellon ja poistamisen kriteerit oletusarvoina pitoisuuksille. Joditabletin käyttö suoritetaan aina automaattisesti, sisällesuojautuminen ja evakuointi voidaan suorittaa joko automaattisesti tietyllä alueella tai annoskriteerin perusteella ja muiden toimenpiteiden käyttöönotto perustuu aina annoskriteeriin paitsi elintarvikkeiden osalta kriteeri ilmoitetaan pitoisuutena tuotteessa. Automaattisten toimenpiteiden käsittelyä havainnol- 13
15 SÄTEILYTURVAKESKUS STUK-YTO-TR 96 listaa kuva 1. Sisällesuojautumiseen ja evakuointiin liittyvät aikatekijät on esitetty kuvassa 2, jonka mukaan on mahdollista käsitellä määrätty viive vastatoimenpiteen aloittamiselle, sisällesuojautumisen kesto, ajoaika alueelta pois, viive alueen A evakuoinnin loppumisesta alueen B evakuoinnin alkamiseen sekä aikaviive ihon puhdistamiseen alueilta A ja B evakuoinnin jälkeen. Vastatoimenpide voi olla joko yksittäinen tai useamman toimenpiteen yhdistelmä. Myös annoskriteeri voi perustua eri annosreittien summaan. Evakuointi tapahtuu vuorokauden kuluessa ja väestön pitempiaikainen poissiirto muutaman viikon kuluessa päästön alusta lukien. Lyhin toimenpiteiden aikaskaala on yksi tunti, mutta usein myös yksi vuorokausi. Suojauskertoimet voidaan määrittää annosreiteittäin normaaleissa ja sisällesuojautumisen olosuhteissa sekä evakuoinnin aikana. Alueet A, B ja S määritellään siten, että kukin napakoordinaatiston piste kuuluu vain yhteen näistä alueista tai sitten kokonaan niiden ulkopuolelle Pelkkä sisällesuojautuminen Kuvassa 1 alue S kuvaa sisällesuojautumista, kun mitään muuta vastatoimenpidettä ei suoriteta. Se aloitetaan annoskriteerillä, joka perustuu joko pilvestä tai laskeumasta tulevaan ulkoiseen annokseen tai hengityksen kautta tulevaan sisäiseen annokseen tai näiden yhdistelmiin. Automaattinen sisällesuojautuminen on mahdollista vain pitkän ajan vaikutuksia laskettaessa käyttäjän määrittelemän ympyrän sisällä. Molemmissa tapauksissa on annettava suojautumisen kestoaika. Annoskriteeriin perustuvassa vaihtoehdossa ohjelma automaattisesti kasvattaa suojautumisaikaa, mikäli pilvi ei ole vielä kulkenut ohi. Tämä voi kasvattaa sisällesuojautumisaikaa olennaisesti alkuperäiseen asetukseen verrattuna Sisällesuojautuminen ja sitä seuraava evakuointi Evakuointi voidaan suorittaa joko automaattisesti kuvan 1 alueella A, jonka määräävät pienempi ympyrän säde r sekä suuremman säteen R ja kulman a rajaama alue tai sitten annokseen perustuen alueella B. Toimenpidealueeksi voidaan valita jompikumpi näistä tai molemmat. Alueen B käsittely on verrannollinen edellisen kohdan alueen S käsittelylle. Evakuointia edeltää aina sisällesuojautuminen. Käyttäjä voi antaa säteille r ja R yhtä suuret arvot, jolloin evakuointialueeksi tulee ympyrä tai antaa r:n arvoksi nollan, jolloin evakuointialueeksi saadaan sektori. Ohjelma valitsee sektorin päästön alussa vallitsevan tuulensuunnan perusteella. Annoskriteeri voi perustua vapaasti valittuun annosyhdistelmään halutulta ajanjaksolta. Jopa eri elimille voidaan asettaa erilaiset annoskriteerit. Kuitenkin samat annosreitit, elimet, integrointiajat ja suojauskertoimet ovat käytössä alueilla B ja S, vaikkakin toimenpidetasot voivat olla erisuuruiset. Evakuointi tapahtuu, jos minkä tahansa elimen annoskriteeri ylittyy. Ohjelma tarjoaa mahdollisuuden kasvattaa annoskriteeriä suuremmilla etäisyyksillä, jotta laajoilla alueilla ei jouduttaisi suurten väestömäärien evakuointiin. Koska toimenpiteet voivat tapahtua nopeammin alueella A, voidaan viiveille antaa erisuuret arvot alueella A ja B. Evakuoinnin alkamisen viiveeksi voidaan antaa myös nolla, joka kuvaa spontaania evakuointia. Evakuoinnin aikana kertyvä annos lasketaan lähtöpaikan annosnopeuden ja vaaditun matkustusajan perusteella. Alueen B evakuointi voi aikaisintaan alkaa vasta kun alueen A evakuointi on loppunut. Ajoaika oletetaan samaksi sekä automaattisessa että annoskriteeriin perustuvassa evakuoinnissa. Sisällesuojautuvan väestön suojaustekijät voidaan antaa kolmessa aikajaksossa: ennen suojautumista, sen aikana sekä evakuoinnin aikana, joka tarkoittaa kulkuvälineiden suojauskerrointa Joditabletti Joditabletin käyttö määräytyy automaattisesti käyttäjän määrittelemän ympyrän sisällä, jonka alan ei tarvitse riippua muiden toimenpiteiden alueista. Joditabletilla rajoitetaan radioaktiivisen jodin kilpirauhaselle aiheuttamaa annosta, joka riippuu altistuksen alkamisen ja tabletin nautinnan välisestä viiveestä. PC COSYMA:ssa tabletin vaikutuksen käsittely on yksinkertais- 14
16 STUK-YTO-TR 96 SÄTEILYTURVAKES KUS Area A: actions: definition: Area B: actions: definition: Area S: actions: definition: sheltering and/or evacuation keyhole shaped area with inner radius r, outer radius R and sector angle a (e.g., r = 2.4 km, R = 5.6 km, a = 60 sheltering and/or evacuation dose to lung, bone marrow Gl-tract, thyroid or effective dose above organ specific thresholds dose to lung, bone marrow Gl-tract, thyroid or effective dose above organ specific thresholds Kuva 1. Välittömien vastatoimenpiteiden mallinnus PC COSYMA:ssa /2/. i,tacc t UAH. Area A TINA TSHA TDR TSKIN TDELB Area B TINB TSHB TDR TSKIN Area A - adions taken automatically Area B - actions taken on the basis of predicted doses TACC = TINA, TINB = TSHA, TSHB = TDELB = TDR = TSKIN = time of accident (end of chain readion) initial delay in areas A and B duration of sheltering in areas A and B delay time between end of evacuation in area A and begin of evacuation in area B driving time to leave area A,B time between leaving area A.B and skin decontamination Kuva 2. Sisällesuojautumisen ja evakuoinnin aikamääreet PC COSYMA.ssa fi./. 15
17 SÄTEILYTURVAKESKUS STUK-YTO-TR 96 tettu siten, että mallille annetaan haluttu kerroin ( 1), jolla kilpirauhasen annos kerrotaan määrätyllä alueella paikasta ja altistuksen ajasta riippumatta Väestön pitempiaikainen poissiirto Väestön pitempiaikainen poissiirto tapahtuu muutaman päivän kuluttua onnettomuudesta ja siten se vaikuttaa myöhäisvaikutusten riskiin. Toimenpide perustuu efektiivisen annoksen kriteeriin täysin itsenäisenä vastatoimenpiteenä, eikä sitä voida yhdistää aiempaan sisällesuojautumiseen. Käytössä olevat annosreitit ovat ulkoinen annos pilvestä ja laskeumasta sekä sisäinen annos hengityksen kautta pilven ja resuspension aiheuttamista maanpintapitoisuuksista. Resuspension ja laskeuman osalta voidaan antaa annoksen integrointiajat (< 1 vuosi), joiden tulisi olla riittävän lyhyet, jotta kriteerit vastaisivat hetkellisiä annosnopeuksia. Hengitysannokset integroidaan vain määrättyyn ajanhetkeen saakka eikä annoskertymän tapaan. Sisäänhengitetystä radioaktiivisesta aineesta aiheutuva annos voi kertyä suurimmaksi osaksi jo lyhyen ajan kuluessa reaktorionnettomuuksien tapauksissa, mutta laskeumasta voi annosta kertyä oleellisesti vielä myöhemminkin. Sen väestön osan, jota ei ole evakuoitu ennen väestönsiirtoa, annokset lasketaan suojauksiltaan normaaliolosuhteiden mukaisesti. Poissiirron kuluessa väestö ei saa annosta. Väestön paluu evakuoinnin tai pitempiaikaisen poissiirron jälkeen vanhalle asuinalueelleen tapahtuu saman annoskriteerin pohjalta. Tällöin tarkastellaan annosta vuoden aika-askelella ja jos se alittaa kriteerin, väestö palaa alueelle. Efektiivinen annos voi olla laskettuna joko ulkoinen annos laskeumasta tai hengitysannos pilvestä tai resuspensiosta tai molempien summa. Uudelleen asuttaminen ja puhdistus käsitellään yhdessä. Puhdistamisen ansiosta takaisinpaluuta voidaan aikaistaa. Ensin verrataan tilannetta suurimman aika-askelen kohdalla, onko annos pienempi kuin kriteeri. Jos näin on, väestö palaa alueelle eikä puhdistusta tarvita. Jos kriteeri ei täyty sovelletaan puhdistusta ja katsotaan täyttyykö kriteeri sitten. Käyttäjä voi määritellä puhdistuksen aikajaksot ja yhden puhdistuskertoimen. Paluun jälkeen sovelletaan normaaleja suojausolosuhteita ja mikäli maaperän puhdistus on suoritettu niin laskeuman aiheuttama annos jaetaan puhdistuskertoimella Elintarvikkeiden käyttökielto Maataloustuotteiden toimenpidetasot annetaan tuotteen pitoisuusarvona nuklidikohtaisesti kahdelle elintarvikeryhmälle: maidolle ja ei-maitotuotteille. Tämä menettely vastaa EU:n toimenpidetasomenettelyä. Tasot kuvaavat elintarvikkeiden pitoisuusarvoja tuoreessa tuotteessa eli ilman mitään viivettä tuotannon ja kulutuksen välillä. Tuotteelle annetaan sekä toimenpideraja että toimenpiteen poistoraja, joiden ei tarvitse välttämättä olla samoja. Nuklidiryhmät ovat: kesium, strontium, jodi ja oc-säteilijät. 16
18 STUK-YTO-TR 96 SÄTEILYTURVAKESKUS 3 OHJELMISTON KÄYTTÖKELPOISUUS Seuraavassa on esitelty käyttäjän kannalta oleellisia havaintoja PC COSYMA:n käytettävyydestä. Ohjelmaversiomme on laajempi 1.0- versio, jonka laskentamahdollisuudet ovat koottuna: 1) ilmapitoisuus ja laskeuma napakoordinaatin pisteessä, 2) annos määritetyltä ajanjaksolta napakoordinaatin pisteessä tai keskimääräisenä annoksena etäisyyspisteessä, 3) väestöannos, 4) varhais- ja myöhäisvaikutusten (fatal ja nonfatal) yksilöriski etäisyyden funktiona, 5) varhais- ja myöhäisvaikutusten (fatal ja nonfatal) lukumäärä väestöjakauman pohjalta, 6) vastatoimenpiteisiin liittyen: vastatoimenpiteen kesto napakoordinaatin pisteessä, henkilöiden ja alueiden määrä ja niiden aikaintegraalit, hylätyn ruoan määrä. Merkittävimmät rajoitukset verrattuna ohjelman pääversioon ovat: 1) leviämismallina on suoraviivainen MUSEMET-malli, 2) vastatoimenpiteiden aiheuttamia taloudellisia seurauksia ei sisälly laskentaan, 3) ravintoaineista ei saada yksilöannoksia, 4) päästöjaksoja on enintään kolme, joiden kesto on kiinteä yksi tunti, 5) nuklidien lukumäärä on 50, 6) vastatoimenpiteet ovat rajoitetumpia mm. seuraavasti: väestön osien suojaus eri ajanjaksoina on yksinkertaistettu, evakuoinnin ja väestönsiirron yksinkertaistettu matkustusmalli, elintarvikkeiden käyttökieltokriteerit vain pitoisuuteen perustuen. Ohjelmalla ei voi laskea suoranaisesti vastatoimenpiteellä säästettyä annosta. Annoskriteerit perustuvat vain annokseen eikä esimerkiksi ulkoiseen annosnopeuteen. Lisäksi vasta-toimenpiteistä aiheutuvien taloudellisten vaikutusten laskenta puuttuu kokonaan. Liitteessä 1 on tarkasteltu taloudellisten vaikutusten laskennan periaatteita ja siihen liittyvien parametrien määritystä COS YMA:ssa. Sovellusosassa on ARANO-ohjelmalla laskettu eräitä tilanteita vertailtavuuden vuoksi ja siten myös havaittuja eroavuuksia näiden kahden ohjelman välillä on pyritty tuomaan esiin. Dispersion osalta ehkä merkittävin eroavuus on se, että PC COSYMA:n MUSEMET-mallissa käytetään pystysuuntaisen sekoittumisen kuvaamiseen perinteistä Gaussin jakaumaa, mutta ARANO:ssa diffuusioyhtälön ratkaisuun perustuvaa ns. K z -mallia, jossa kuivadepositio vaikuttaa pystysuuntaisen profiilin muotoon. Taulukoissa VII ja VIII on esitetty kaikki aktiivisessa käytössä olevat kotimaiset ilmaleviämistä ja ympäristön seurausvaikutuksia laskevat tietokoneohjelmat tärkeimpien ominaisuuksiensa mukaan eriteltynä. Kuten taulukosta käy ilmi on ARANO ainoa PC COSYMA:an verrattavissa oleva kotimainen ohjelma. Kansainvälisellä tasolla eri maiden käyttämiä ohjelmia on vertailtu OECD/NEA:n ja CEC:n järjestämässä vertailussa, jossa COSYMA.n pääversio ja ARANOohjelma olivat mukana neljän muun ohjelman kanssa 11,%/. PC COSYMA.n käyttöliittymä perustuu hiirellä ja näppäimistöllä ohjattaviin valikoihin ja siten sen käyttö on visuaalisesti helppoa. Sekä lähtötietojen syöttö että tulosten katselu tapahtuvat valikkonäyttöjen kautta. Valikoiden käytön oppii ehkä parhaiten "istu ja kokeile" -menetel- 17
19 SÄTEILYTURVAKESKUS STUK-YTO-TR 96 Taulukko VII. Kotimaiset radioaktiivisten päästöjen ilmeleviämis- ja annoslaskentamallit. Malli Kehittäjä Käyttäjä Leviämismalli. ulottuvuus Annoslaskenta Vastatoimenpiteet Tervey svai kut u kset Taloudelliset vaikutukset Deterministiset, todennäk. laskut sääaineistolla Käyttötarkoitus ARANO IL, VTT VTT suoraviivainen K- lähialue, mesoskaala pilvi, laskeuma, hengitys, ravintoaineet suojautuminen, evakuointi, joditabletit, väestönsiirto, maaperän puhdistus varh.+myöh. vastatoimien kustannukset det + tod tutkimus OIVA STUK, IL STUK trajektori K 2, lähialue, mesoskaala pilvi, laskeuma, hengitys - det valmiusohjelma, tutkimus ROSA VTT IVO, TVO, VTT suoraviivainen K, lähialue pilvi, laskeuma, hengitys - det reaaliaikainen valmiusohjelma harjoittelu TRADOS IL, VTT IL, VTT trajektori K 7, mesoskaala, kaukokulkeutuminen pilvi, laskeuma, hengitys, ravintoaineet - det + tod tutkimus, valmuisohjelma TUULET IVO IVO, TVO suoraviivainen Gauss, lähialue, mesoskaala pilvi, laskeuma, hengitys, ravintoaineet - det+tod tutkimus Taulukko VIII. Mallien ohjelma- ja laiteominaisuudet. Malli ARANO OIVA ROSA TRADOS TUULET Ohjelmointikieli Fortran Fortran Fortran, C Fortran, C Fortran Käyttöympäristö PC/DOS/OS2 VAX/VMS PC/OS2, VAX/VMS Työasema/ UNIX VAX/VMS, PC/DOS Käyttöliittymä tiedostot VMS virtual + grafiikka GKS/CGI C/Motif tiedostot 18
20 STUK-YTO-TR 96 SÄTEILYTURVAKESKUS mällä, joka on yleistynyt kaikissa kaupallisissa PC-ohjelmissa. Tämän menetelmän periaate on, että varsinaisten ohjelman käyttöön opastavien manuaalien käyttö ei ole tarpeen vaan käyttöliittymä tarjoaa riittävästi opastusta. Ohjelman käyttäjille on järjestetty myös maksullisia noin viikon kestäviä kursseja. Kaikilla parametreillä on oletusarvot ja niitä muuttamalla käyttäjä voi suorittaa melko sujuvasti herkkyystarkasteluja. Kuitenkin parametrien sallitut arvot puuttuvat, mikä voi aiheuttaa turhaa epätietoisuutta niiden vaihtelualueen laajuudesta. Annetut lähtötiedot voi tallettaa tiedostoksi (*.SAV), jota voi käyttää jatkossa pohjana uusia syöttötiedostoja muodostettaessa. Sen sijaan laskentatuloksina säilyvät vain viimeisimmän ajon tulokset. Ohjelman käyttöön liittyvä manuaali on myös hyvin käyttäjäystävällinen 191. Kuitenkin ohjelmaa käytettäessä voi törmätä kysymyksiin, joihin vastausta joutuu hakemaan KfK:n tai NRPB:n COSYMAraporteista. Niiden mallien kuvaukset koskevat kuitenkin ohjelman pääversiota, joten aina ei niistäkään saata löytää vastausta kysymyksiin. COSYMAm käyttäjille tarkoitettu lehti "COSY- MA Newsletter" ilmestyy myös tarpeen mukaan muutaman kerran vuodessa /10/. Lehdessä esitetään ajankohtaisia asioita COSYMA-ohjelmien kehityksestä ja sovelluksista. On myöskin perustettu ohjelman käyttäjien ryhmä "COSYMA Users Group", johon kokoonnutaan noin kerran vuodessa vaihtamaan mielipiteitä mallin kehittäjien ja käyttäjien välillä 711,12/. 19
SÄTEILYTURVAKESKUS STUK-YTO-TR 64. Jukka Rossi HELMIKUU 1994. Strälsäkerhetscentralen r j-inry/ Finnish Centre for Radiation and MIIAI Nuclear Safety
Loppuraportti osallistumisesta OECD/NEA:n ja CEC:n organisoimaan kansainväliseen reaktorionnettomuuksien ympäristöseurauksia laskevien ohjelmien vertailuun Jukka Rossi HELMIKUU 1994 SÄTEILYTURVAKESKUS