Radiochemia świeżego wypalonego paliwa: Np-239, I-131 i krótkotrwałe radionuklidy

Streszczenie

Wypalone paliwo nie jest jedną substancją o stałej "radioaktywności". Tuż po wyjęciu z rdzenia jest mieszaniną setek radionuklidów, z których część znika w minutach i godzinach, część w dniach, a część zostaje na dekady, tysiąclecia albo dłużej. Dlatego świeże wypalone paliwo ma inny podpis radiochemiczny niż paliwo chłodzone przez miesiące lub lata.^1,2

Ten artykuł opisuje krótkotrwałe radionuklidy jako zegary napromieniania i chłodzenia. Szczególną rolę pełnią Np-239, powstający po wychwycie neutronu przez U-238, oraz produkty rozszczepienia takie jak I-131, Xe-133, Ba-140/La-140, Zr-95/Nb-95 i inne nuklidy, które pozwalają czytać historię paliwa bez wchodzenia w procedury przerobu chemicznego.

Rozszerzenie tematu

Świeże paliwo wypalone jako zapis czasu

W reaktorze paliwo jest jednocześnie źródłem energii, celem napromieniania i archiwum historii neutronowej. Każde rozszczepienie tworzy produkty rozszczepienia. Każdy wychwyt neutronu przez U-238 może rozpocząć drogę do transuranowców. Każdy dzień pracy rdzenia zmienia proporcje radionuklidów.

Po wyłączeniu reaktora reakcja łańcuchowa ustaje, ale inwentarz promieniotwórczy nie znika. Produkty rozszczepienia i aktynowce rozpadają się według własnych czasów połowicznego zaniku. To właśnie dlatego czas chłodzenia jest jednym z najważniejszych parametrów interpretacji wypalonego paliwa.

Świeże wypalone paliwo, mierzone w skali godzin lub dni po zakończeniu napromieniania, zawiera silny wkład krótkotrwałych nuklidów. Po miesiącach dominują inne grupy. Po latach coraz ważniejsze stają się Cs-137, Sr-90, aktynowce i produkty aktywacji. Jedna próbka, ale wiele różnych "wieków" radiochemicznych.

Co znaczy "zegar radionuklidowy"

Zegar radionuklidowy to nie zegarek z jedną wskazówką. To para albo grupa nuklidów, których proporcje zmieniają się w przewidywalny sposób. Jeśli znamy sposób produkcji i czasy połowicznego zaniku, możemy wnioskować o czasie od napromieniania albo od rozładowania paliwa.

Najprostszy przykład to pojedynczy nuklid po zakończeniu produkcji:

A(t) = A0 * exp(-lambda*t).

Po jednym czasie połowicznego zaniku zostaje połowa aktywności. Po dwóch zostaje jedna czwarta. Po dziesięciu zostaje około jedna tysięczna. Dlatego nuklid o czasie połowicznego zaniku 8 dni jest ważny w pierwszych tygodniach, ale po roku staje się śladem bardzo małym w porównaniu z długowiecznymi składnikami.

W paliwie reaktorowym dochodzi produkcja podczas napromieniania. Wtedy proste "liczenie od zera" nie wystarcza. Trzeba uwzględnić, że nuklid jednocześnie powstawał i zanikał. Mimo to podstawowa intuicja pozostaje: krótkotrwałe nuklidy mówią o świeżej historii, długotrwałe o całkowitym inwentarzu i długim chłodzeniu.

Np-239: pomost od U-238 do plutonu

Jednym z najważniejszych łańcuchów w historii techniki jądrowej jest:

U-238 + n -> U-239 -> Np-239 -> Pu-239.

W materiałach Nuclear Weapons FAQ oraz w artykule o Hanford pokazano ten łańcuch jako podstawę produkcji plutonu w reaktorze: U-239 rozpada się beta z czasem połowicznego zaniku około 23,5 min, a Np-239 rozpada się beta do Pu-239 z czasem około 2,35 dnia.^3,4

Z punktu widzenia radiochemii Np-239 jest znakomitym wskaźnikiem świeżości napromieniania. Jeśli paliwo albo próbka uranu była niedawno w strumieniu neutronów, Np-239 może być widoczny. Jeżeli minęło wiele tygodni bez dalszej produkcji, większość Np-239 przejdzie już do Pu-239.

To nie znaczy, że Np-239 sam "daje pluton do ręki". Oznacza, że jego obecność w widmie i bilansie jest sygnałem niedawnego wychwytu neutronów przez U-238. Dla dydaktyki jest to świetny przykład zegara pośredniego: nuklid jest ważny nie jako produkt końcowy, lecz jako krótkożyciowe ogniwo łańcucha.

Dlaczego Np-239 jest ważny historycznie

W programie Manhattan kluczowe było nie tylko odkrycie plutonu, lecz także zrozumienie, że można go produkować w reaktorze przez napromienianie uranu. Artykuł o Chicago Pile-1 przypomina, że Seaborg i współpracownicy obserwowali powstawanie Np-239, a następnie jego rozpad do Pu-239, zanim reaktory produkcyjne umożliwiły skalę kilogramową.⁵

W Hanford ten sam łańcuch przestał być ciekawostką laboratoryjną i stał się osią przemysłową. Reaktor produkował pluton przez wychwyt neutronów w U-238, a niskie wypalenie ograniczało dalsze wychwyty prowadzące do Pu-240.⁴

Dla tego artykułu ważny jest jednak aspekt metrologiczny: Np-239 pokazuje, że skład materiału po napromienianiu nie jest statyczny. Nawet po wyjęciu z reaktora zachodzą rozpady beta, które zmieniają izotopowy obraz próbki w skali dni.

I-131: zegar tygodniowy i problem biologiczny

I-131 jest produktem rozszczepienia o czasie połowicznego zaniku około 8 dni. W świeżym wypalonym paliwie i świeżym opadzie promieniotwórczym może mieć duże znaczenie radiologiczne. Po kilku miesiącach jego aktywność spada o wiele rzędów wielkości.

To czyni I-131 zegarem tygodniowym. Jeśli w próbce środowiskowej widzimy istotny I-131, zwykle pytamy o niedawne uwolnienie, niedawny proces rozszczepienia albo źródło medyczne. Jeśli próbka pochodzi z dawno chłodzonego paliwa, I-131 nie powinien być jednym z głównych składników.

I-131 ma też znaczenie biologiczne, ponieważ jod jest wychwytywany przez tarczycę. Artykuły o skutkach opadu promieniotwórczego zwracają uwagę na drogę przez mleko i szczególną wrażliwość dzieci.⁶ W tekście o wypalonym paliwie ten kontekst pomaga zrozumieć, dlaczego krótkożyciowy nuklid może być bardzo ważny mimo tego, że po roku prawie znika.

Nuklidy godzinowe i dzienne

Najświeższy podpis po rozszczepieniu tworzą nuklidy o czasach połowicznego zaniku liczonych w godzinach i dniach. W lokalnym kalkulatorze inwentarza wypalonego paliwa znajdują się między innymi Kr-88, Sr-91, Zr-97, I-132, I-133, I-135, Xe-133, Xe-135, Ba-141, La-141 i Ce-143.⁷

Tak krótkie czasy życia sprawiają, że po dłuższym chłodzeniu nuklidy te przestają dominować. Ich obecność jest natomiast cenna w analizie bardzo świeżego materiału, awarii, testu lub eksperymentu aktywacyjnego.

Ważne jest rozróżnienie dwóch zjawisk:

• nuklid może być krótkożyciowy, ale stale produkowany podczas pracy reaktora,
• po zakończeniu produkcji jego aktywność szybko maleje.

Dlatego próbka pobrana w czasie pracy albo tuż po wyłączeniu wygląda inaczej niż próbka badana po miesiącu.

Pary rodzic-córka

W wypalonym paliwie wiele sygnałów pochodzi z par lub łańcuchów rodzic-córka. Przykłady dydaktyczne to Te-132/I-132, Ba-140/La-140, Sr-90/Y-90 i Cs-137/Ba-137m. Córka może mieć krótki czas połowicznego zaniku, ale być stale zasilana przez dłużej żyjącego rodzica.

To prowadzi do równowagi promieniotwórczej. Jeśli rodzic żyje znacznie dłużej niż córka, aktywność córki może po pewnym czasie zbliżyć się do aktywności rodzica. Wtedy obecność krótkotrwałej córki nie oznacza koniecznie świeżego rozszczepienia. Może oznaczać, że istnieje dłużej żyjący rodzic.

Przykład Cs-137/Ba-137m jest znany ze spektrometrii gamma: krótkotrwały stan metastabilny Ba-137m daje linię około 662 keV, ale źródłem długiego trwania sygnału jest Cs-137. W analizie wypalonego paliwa podobna logika jest konieczna dla wielu par.

Skala miesięcy: Zr-95, Nb-95, Ru-103, Ce-144

Po pierwszych dniach i tygodniach nadal pozostaje wiele nuklidów o czasach połowicznego zaniku liczonych w tygodniach lub miesiącach. W dydaktycznym inwentarzu kalkulatora pojawiają się między innymi Zr-95, Nb-95, Ru-103, Ba-140, Ce-141, Ce-144, Pm-148m i inne produkty rozszczepienia.⁷

Te nuklidy są użyteczne jako pośrednie zegary chłodzenia. Nie znikają tak szybko jak I-135, ale też nie są dominującym problemem po dekadach. Ich proporcje pomagają odróżnić materiał chłodzony krótko od materiału, który przeszedł już przez pierwszy etap rozpadu krótkotrwałych produktów.

W praktyce pełna interpretacja wymaga modelu inwentarza: wypalenia, historii mocy, udziału rozszczepień U-235, Pu-239 i Pu-241, czasu chłodzenia oraz danych jądrowych. Dydaktycznie wystarczy uchwycić kierunek: im dłuższe chłodzenie, tym bardziej przesuwa się znaczenie od nuklidów godzinowych i dniowych ku miesięcznym, rocznym i długowiecznym.

Skala lat i dekad: Cs-137 i Sr-90

Po latach chłodzenia wciąż bardzo ważne są Cs-137 i Sr-90, oba o czasach połowicznego zaniku rzędu trzydziestu lat. To one często dominują w dyskusjach o odpadach i skażeniach długoterminowych w skali ludzkiego życia.

Cs-137 jest łatwy do śledzenia spektrometrią gamma dzięki linii potomnego Ba-137m. Sr-90 jest trudniejszy w bezpośredniej spektrometrii gamma i wymaga innego podejścia analitycznego. Ta różnica pokazuje, że "ważny radionuklid" nie zawsze znaczy "łatwy do zmierzenia tym samym detektorem".

W kontekście świeżego wypalonego paliwa Cs-137 i Sr-90 nie są jedynymi bohaterami, ale tworzą tło długowieczne. Gdy krótkotrwałe nuklidy znikają, te izotopy nadal pozostają.

Aktynowce: mniej aktywności chwilowej, więcej długiego ogona

Produkty rozszczepienia często dominują aktywność i ciepło rozpadu tuż po wyłączeniu. Aktynowce są inne. Pluton, neptun, ameryk i kiur tworzą długi problem materiałowy i radiotoksyczny. Ich znaczenie rośnie w dyskusji o wieloletnim chłodzeniu, przerobie, MOX, transmutacji i składowaniu.

Np-239 jest krótkotrwałym ogniwem prowadzącym do Pu-239. Np-237 jest zupełnie innym problemem: długożyciowym aktynowcem mniejszym. Pu-241 rozpada się do Am-241, co zmienia ciepło, aktywność alfa i sygnaturę materiału w czasie. Dlatego pojęcie "neptun" albo "pluton" bez numeru masowego jest w radiochemii zbyt ogólne.

Artykuł o MOX podkreśla, że przerób paliwa nie cofa historii izotopowej plutonu. Skład plutonu jest wynikiem wypalenia i chłodzenia, a kolejne izotopy oraz produkty rozpadu zmieniają właściwości materiału.⁸

Ciepło powyłączeniowe

Krótkotrwałe produkty rozszczepienia są ważne nie tylko dlatego, że są aktywne. Ich rozpady beta i gamma wytwarzają ciepło powyłączeniowe. Reaktor po wyłączeniu nie produkuje już mocy z kontrolowanej reakcji łańcuchowej, ale paliwo nadal grzeje się od rozpadu promieniotwórczego.

To jest powód, dla którego świeże wypalone paliwo wymaga chłodzenia i osłon. W pierwszych godzinach i dniach malejący, ale nadal znaczący wkład krótkotrwałych nuklidów decyduje o wymaganiach cieplnych i radiologicznych.

Kalkulator burnup w serwisie modeluje wypalanie paliwa, ciepło powyłączeniowe po SCRAM oraz uproszczony inwentarz ORIGEN-lite z rankingami aktywności, ciepła i wskaźników gamma.⁹ To dobry kierunek edukacyjny: zamiast podawać jedną liczbę "radioaktywności", pokazuje się zmienność w czasie.

Spektrometria jako narzędzie odczytu historii

Widmo gamma świeżego materiału po rozszczepieniu może być bardzo złożone. Linie wielu nuklidów nakładają się, tło komptonowskie jest wysokie, a aktywności zmieniają się w czasie. Analiza nie polega więc na prostym odczytaniu jednego piku.

Trzeba rozdzielić:

• identyfikację linii gamma,
• odejmowanie tła,
• wydajność detektora,
• czas od zakończenia napromieniania,
• czas pomiaru,
• rozpady podczas pomiaru,
• równowagi rodzic-córka,
• samopochłanianie i geometrię próbki,
• historię produkcji nuklidów.

Właśnie dlatego temat pasuje do działu metrologii. Radiochemia świeżego paliwa jest mniej "katalogiem izotopów", a bardziej problemem odwrotnym: z widma i modelu próbujemy odtworzyć historię.

Co wolno pokazywać dydaktycznie, a czego unikać

Ten temat leży blisko technologii przerobu paliwa, więc trzeba wyraźnie oddzielić warstwę edukacyjną od operacyjnej. W artykule wolno i warto omawiać:

• rozpady i czasy połowicznego zaniku,
• produkcję radionuklidów podczas napromieniania,
• jakościową interpretację widm,
• zależność inwentarza od czasu chłodzenia,
• różnicę między produktami rozszczepienia i aktynowcami,
• rolę krótkotrwałych nuklidów w cieple powyłączeniowym.

Nie należy natomiast publikować instrukcji separacji chemicznej, receptur przerobu, sekwencji odzysku plutonu, parametrów procesowych gorących komór ani wskazówek umożliwiających praktyczne obchodzenie się ze świeżym wypalonym paliwem. To nie jest potrzebne do zrozumienia fizyki zegarów radionuklidowych.

Prosty model: nuklid po zakończeniu produkcji

Załóżmy syntetyczny nuklid X o czasie połowicznego zaniku T1/2 = 8 dni. Po zakończeniu produkcji jego aktywność maleje:

A(t) = A0 * 2^(-t/T1/2).

Po 8 dniach zostaje 50%, po 16 dniach 25%, po 80 dniach około 0,1%. To pokazuje, dlaczego nuklid podobny czasowo do I-131 jest znakiem świeżości, a nie długoterminowym składnikiem po latach.

Jeżeli drugi nuklid ma T1/2 = 30 lat, jego aktywność po 80 dniach prawie się nie zmieni. Porównanie tych dwóch krzywych jest pierwszym krokiem do zrozumienia, dlaczego widmo paliwa starzeje się radiochemicznie.

Prosty model: Np-239 po napromienianiu

Po zakończeniu napromieniania U-239 szybko przechodzi do Np-239, a Np-239 dalej do Pu-239. Jeżeli pominiemy dalszą produkcję, aktywność Np-239 maleje w skali kilku dni. Po jednym czasie połowicznego zaniku zostaje połowa, po pięciu zostaje około 3%.

To daje intuicję: Np-239 jest znaczący bardzo krótko po napromienianiu, natomiast Pu-239 pozostaje jako produkt długowieczny. W analizie historycznej i safeguards obecność lub brak Np-239 może więc mówić o czasie, a skład plutonu o historii wypalenia.

Ograniczenia prostych zegarów

Żaden pojedynczy radionuklid nie rozwiązuje całej historii paliwa. Proste zegary zawodzą, jeśli:

• historia mocy była zmienna,
• paliwo miało niejednorodne wypalenie,
• część nuklidów była lotna i mogła migrować,
• próbka nie reprezentuje całego elementu paliwowego,
• doszło do separacji chemicznej,
• pomiar wykonano po czasie porównywalnym z rozpadem analizowanego nuklidu,
• użyto niepełnej biblioteki danych jądrowych.

Dlatego profesjonalna analiza korzysta z modeli inwentarza, bibliotek przekrojów i danych o historii rdzenia. W serwisie edukacyjnym można jednak zacząć od prostych wykresów: aktywność w funkcji czasu dla kilku nuklidów o różnych czasach połowicznego zaniku.

Metody pomiaru krótkożyciowych radionuklidów

Pomiar inwentarza świeżego wypalonego paliwa jest trudny eksperymentalnie. Aktywność tuż po wyłączeniu może wynosić miliardy bekereli na gram, a szybkie zmiany składu wymagają pomiaru w precyzyjnie określonym czasie. Stosuje się kilka komplementarnych technik.

Spektrometria gamma jest główną metodą nieinwazyjną. HPGe (germano wysoko-czysty) pozwala rozdzielić linie gamma nakładające się na widmo złożone. Silne promieniowanie gamma wymaga jednak grubych osłon ołowianych lub wolframowych i zazwyczaj znacznego rozcieńczenia lub geometrii dalekiej (standoff). Aktywność zmienia się podczas pomiaru, gdy T₁/₂ jest krótki — konieczna korekta na rozpad integracyjny. Normalizacja do czasu pomiaru stosuje wzór Bateman — dla nuklidu zanikającego podczas czasu t₀ do t₀+Δt pole piku jest proporcjonalne do A₀ × (1 − e^(−λΔt)) / λ.

Pomiary beta (liczniki proporcjonalne, detektory plakietkowe) mierzą strumień elektronów i pozytronów. Są mniej selektywne niż gamma, ale czulsze dla nuklidów emitujących wyłącznie beta (np. Sr-90, Y-90, Pm-147). W aktywnym rdzeniu spektrometria beta jest praktycznie niemożliwa z powodu tła gamma; stosuje się ją po separacji radiochemicznej.

Techniki radiochemiczne — mineralizacja, separacja na kolumnach (TEVA, TRU, Eichrom SR) i pomiar spektrometryczny poszczególnych frakcji — są pracochłonne, ale dają selektywność niedostępną spektrometrycznie. Środki ostrożności przy pracy z materiałem świeżo napromieniowanym są znaczne (gorące komory, manipulatory zdalne). Akademickie opisy tych technik są dopuszczalne na poziomie zasady i nie zawierają operacyjnych receptur.

Neutrony (detekcja przez He-3, BF₃, komory szczelinowe) mierzą strumień neutronów z opóźnionego rozszczepienia lub reakcji (γ,n) w jądrach ciężkich. Technika stosowana w kontroli procesu i pomiarach kriticality.

Produkcja radionuklidów podczas napromieniania — równania Batemana

Dla nuklidów produkowanych neutronami (aktywacja lub rozszczepienie) obowiązuje równanie różniczkowe łączące produkcję P i rozpad λ:

dN/dt = P − λN

Gdzie N jest liczbą atomów (lub aktywnością przez przeliczenie A = λN). W stanie stacjonarnym (tzw. saturacji): N_sat = P/λ, czyli A_sat = P. Po zakończeniu napromieniania (P = 0): A(t) = A₀ × e^(−λt).

Dla pary rodzic-córka (nuklid 1 → nuklid 2):

dN₂/dt = λ₁N₁ − λ₂N₂

Rozwiązanie analityczne to klasyczne równanie Batemana. Jeżeli T₁/₂(rodzic) >> T₁/₂(córka), po czasie wystarczająco długim aktywność córki zbliża się do aktywności rodzica (równowaga wiekowa). Jeżeli T₁/₂(rodzic) jest porównywalne z T₁/₂(córka), proporcja oscyluje w czasie.

Przykład Ba-140 (T₁/₂ = 12,75 d) → La-140 (T₁/₂ = 1,68 d): La-140 osiąga maksimum aktywności ok. 3 dni po izolacji Ba-140, następnie aktywność La-140 maleje razem z Ba-140. To klasyczny przykład równowagi przejściowej.

Kody obliczeniowe: ORIGEN i ORIGEN-S

Profesjonalna analiza inwentarza paliwa nie polega na prostym sumowaniu klonowalnych wzorów. Wymaga kodu numerycznego uwzględniającego setki łańcuchów rozpadów równoległych i sekwencyjnych, przekroje czynne zależne od widma neutronów oraz historię mocy reaktora.

ORIGEN (Oak Ridge Isotope Generation and Depletion Code) opracowany w ORNL jest od lat 1970. standardem przemysłowym. Wersja ORIGEN-S jest składnikiem pakietu SCALE (Standardized Computer Analysis for Licensing Evaluation). Podstawowe dane wejściowe: wypalenie (MWd/tHM), skład paliwa, wzbogacenie, geometria, widmo neutronów. Wynik: aktywność, ciepło, toksyczność i emisja neutronów dla każdego nuklidu w zadanym czasie chłodzenia.

FISPACT (EAF, EAF-2010, FISPACT-II) jest alternatywą europejską, stosowaną przy aktywacji materiałów strukturalnych reaktorów fuzyjnych i fisyjnych. Dane jądrowe z bibliotek TENDL i JEFF.

Z dydaktycznego punktu widzenia ważne jest zrozumienie, że kody te nie „wymyślają" inwentarza — korzystają z danych jądrowych (przekroje, stałe rozpadu, gałęzie rozpadów) zbieranych przez dekady w bazach ENDF, JEFF, JENDL. Niepewności tych danych propagują się do niepewności obliczonego inwentarza.

Radiochemia paliwa w programach safeguards MAEA

Międzynarodowa Agencja Energii Atomowej (MAEA) używa analizy składu izotopowego paliwa (tzw. destructive assay i non-destructive assay) jako instrumentu safeguards. Cel: potwierdzić, że materiał nie był przetaczany przez separację chemiczną i że jego historia jest zgodna z deklaracją.

Wypalenie paliwa można niezależnie wyznaczyć z kilku metod:

• stosunek Nd-148/U: Nd-148 jest produktem rozszczepienia akumulującym się z wypaleniem; jego masowe stężenie w stosunku do uranu rośnie monotonicznie;
• stosunek Cs-137/Nd-148: kombinacja pozwala na wzajemną weryfikację;
• Pu-238/(Pu-239+Pu-240): rośnie z wypaleniem, bo Pu-238 pochodzi z wielokrotnych wychwytów;
• stosunek Pu-240/Pu-239: rośnie z wypaleniem — stąd wniosek, że paliwo „reaktorowe" ma zupełnie inny skład plutonu niż paliwo „do broni" z krótkim napromienianiem.

Polska MAEA uczestniczy jako uczestnik Układu NPT i poddaje swoje instalacje regularnym inspekcjom MAEA. NCBJ Świerk jest objęty porozumieniem o zabezpieczeniach z MAEA.

CTBT i monitorowanie świeżych produktów rozszczepienia

System monitoringu CTBT (opisany osobno w artykule o IMS) używa ksenonu jako markera testu jądrowego. Ksenon jest chemicznie obojętny, nie absorbuje się w glebie i szybko migruje w atmosferze. Kluczowe izotopy:

Izotop	T₁/₂	Znaczenie
Xe-131m	11,9 d	Naturalny i reaktorowy; może maskować sygnał testu
Xe-133	5,24 d	Dominujący produkt rozszczepienia; główny marker
Xe-133m	2,19 d	Metastabilny; stosunek Xe-133m/Xe-133 pomaga datować zdarzenie
Xe-135	9,1 h	Bardzo świeży marker; zanika szybko

Świeży test daje charakterystyczny stosunek aktywności: Xe-135 >> Xe-133 >> Xe-131m w pierwszych dniach; po tygodniach hierarchia odwraca się. Na tej podstawie system CTBT może odróżnić wykrycie tuż po teście od wykrycia odległego w czasie.

I-131 (T₁/₂ = 8 d), Cs-137 i Ba-140/La-140 są też markerami, ale ich transport atmosferyczny jest bardziej skomplikowany (depozycja sucha i mokra, inne czasy życia w atmosferze). Ksenon zachowuje się jak gaz szlachetny i jest preferowanym analitem do wczesnego ostrzegania.

Polska perspektywa: doświadczenia po Czarnobylu i Fukushimie

Polska stacja CTBT (IFJ PAN Kraków, Belsk-IFJ) prowadziła ciągłe pomiary stężenia radionuklidów w powietrzu. W 1986 roku opady z Czarnobyla dotarły do Polski 28 kwietnia — I-131, Cs-137, Cs-134, Te-132/I-132, Ba-140/La-140 były wyraźnie widoczne w pomiarach CLOR i IMiGW. Proporcje aktywności tych nuklidów pozwoliły datować zdarzenie i odróżnić wkład różnych faz pożaru reaktora.

Po Fukushimie w 2011 roku I-131 dotarł do Europy w ultra-niskich stężeniach (rzędu mBq/m³), ale HPGe o niskim tle i długie czasy pomiaru pozwoliły na detekcję. Proporcja I-131/Cs-137 była markierem czasu transportu: krótka (dni) oznacza świeżą emisję, długa (tygodnie) oznacza transport z odległego źródła.

Polskie laboratoria (CLOR, IFJ PAN, NCBJ) publikowały swoje wyniki monitoringu i uczestniczyły w europejskich programach porównań (EURDEP, RODOS) po obu zdarzeniach.

Safeguards jądrowe a analiza składu izotopowego paliwa

Oprócz monitoringu środowiskowego analiza radiochemiczna wypalonego paliwa jest centralnym narzędziem systemu zabezpieczeń (safeguards) MAEA. Cel: potwierdzić, że deklarowana historia napromieniania jest prawdziwa i że nie doszło do nielegalnej separacji materiałów rozszczepialnych.

W tym celu stosuje się dwie klasy metod: destructive assay (DA) i non-destructive assay (NDA). W DA rozpuszcza się fragment paliwa i bada skład izotopowy Pu, U, Nd, Cs, Ba chemicznie i spektrometrycznie. W NDA mierzy się emisję gamma i neutronów bez roztwarzania. Obie metody mają swoje typowe niepewności i ograniczenia.

Konsystencja między deklarowanym wypaleniem a zmierzonymi stosunkami izotopowymi dostarcza niezależnego potwierdzenia. Na przykład: stosunek Pu-240/Pu-239 rośnie monotonicznie z wypaleniem. Jeśli deklarowane wypalenie wynosi 45 GWd/tHM, a zmierzony stosunek Pu-240/Pu-239 odpowiada tylko 10 GWd/tHM, istnieje rozbieżność wymagająca wyjaśnienia. Takie niespójności są sygnałem w systemie MAEA.

Pu-241 (T₁/₂ = 14,3 lat, beta) jest ważnym nuklid przekraczającym granicę między plutonowcami a americem. Przez rozpad Pu-241 → Am-241 skład wypalonego paliwa zmienia się nawet po długim chłodzeniu. Am-241 akumuluje się przez dekady. Laboratorium mierzące paliwo chłodzone 10 lat musi uwzględnić inkrementalny wzrost Am-241 i związany z nim wzrost aktywności alfa. Nieuwzględnienie tego powoduje zawyżenie stosunku Pu-241/Am-241 i potencjalnie zły wniosek o wieku materiału.

Forensyka jądrowa — od składu izotopowego do historii materiału

Forensyka jądrowa (nuclear forensics) stosuje analizę izotopową, pierwiastkową i chemiczną do odtworzenia historii próbki nieznanego pochodzenia. Dla materiałów wypalonego paliwa podstawowe pytania to:

1. Z jakiego rodzaju reaktora pochodzi materiał? (LWR, CANDU, RBMK, FBR — różne widma neutronów dają różny stosunek Pu-240/Pu-239)
2. Jakie było wypalenie? (Nd-148/U, Cs-137 jako sygnatury)
3. Kiedy paliwo było ostatnio w reaktorze? (czas chłodzenia z krótkożyciowych nuklidów i modeli Batemana)
4. Czy nastąpiła separacja chemiczna? (brak produktów rozszczepienia przy obecności aktynowców byłby anomalią)
5. Jaka jest historia wzbogacenia uranu i typ reaktora? (stosunek U-234/U-235 i U-236/U-235 wskazuje na reaktorowe czy wzbogacane paliwo)

W typowym reaktorze LWR (lekkowodnym) widmo neutronów jest termiczne. Reaktory ciężkowodne (CANDU) mają wyższą konwersję U-238 → Pu-239. Reaktory na neutronach prędkich (FBR) produkują pluton o wyższym stosunku Pu-240/Pu-239 niż LWR przy tym samym wypaleniu. Te fingerprints pozwalają wnioskować o typie reaktora z składu izotopowego plutonu.

Polska nie prowadzi oddzielnego programu forensyki jądrowej, ale IFJ PAN i NCBJ posiadają kompetencje spektrometryczne i radiochemiczne niezbędne do wykonywania tego rodzaju analiz w razie potrzeby kryminalistycznej lub na zlecenie MAEA. Baza danych fingerprints reaktorów (IAEA Nuclear Forensics Library) gromadzi sygnatury izotopowe z różnych instalacji i jest dostępna dla organów śledczych państw członkowskich MAEA. Polska jako uczestnik Traktatu NPT i Konwencji o zabezpieczeniach fizycznych (CPPNM i jego zmienionego protokołu z 2005 r.) jest zobowiązana do współpracy w przypadku incydentów z materiałem jądrowym.

Tabela: dominujące grupy nuklidów a czas chłodzenia

Poniżej zestawienie, które grupy nuklidów dominują aktywność i ciepło wypalonego paliwa LWR w różnych fazach chłodzenia:

Czas chłodzenia	Dominujące nuklidy aktywności	Dominujące nuklidy ciepła	Uwaga
Sekundy–minuty	Kr-88, I-135, Xe-135, Ba-141, La-141	Te same krótkotrwałe + Sr-91, Y-91	Krytycznie ważne dla LOCA
Godziny–1 dzień	I-131, I-132, I-133, Xe-133, Np-239	I-131..133, Te-132, Ba-140, Np-239	Zegar Np-239 aktywny
1–7 dni	I-131, Ba-140, La-140, Ru-103, Te-132	Ba-140/La-140, Ru-103/Rh-103m, Nb-95	Równowaga Ba/La ważna
1–4 tygodnie	Zr-95/Nb-95, Ce-144/Pr-144, Ru-103	Zr-95, Ru-103, Ce-144, Ba-140 (zanika)	I-131 maleje 6-krotnie/tydzień
Miesiące	Ce-144, Ru-106, Pm-147, Cs-134, Cs-137	Cs-134, Ce-144, Pm-147, Ru-106	Sr-90/Y-90 ważne dla beta
Lata	Cs-137, Sr-90, Pm-147, Cs-134	Cs-137, Sr-90, Cm-244, Am-241	Aktynowce ważniejsze
Dekady	Cs-137, Sr-90, aktynowce	Am-241, Cs-137, Sr-90	Cs-137 i Sr-90 przez 300 lat
Tysiąclecia	Aktynowce: Pu-239, Np-237, Am-241	Pu-239, Am-241, Pa-231	Problem składowania długookresowego

Tabela jest poglądowa — rzeczywiste proporcje zależą od wypalenia, wzbogacenia i historii mocy. Służy do demonstracji, że nie ma jednego „dominującego radionuklidu" niezależnego od kontekstu czasowego.

Ciepło powyłączeniowe a bezpieczeństwo reaktora

Ciepło powyłączeniowe (decay heat) jest parametrem bezpieczeństwa klasy pierwszej. W projekcie reaktora wymagania dotyczące chłodzenia awaryjnego (ECCS — Emergency Core Cooling System) muszą uwzględniać moc ciepła powyłączeniowego jako funkcję czasu.

Regulacje wymaga, by projekt ECCS był wystarczający do chłodzenia paliwa z pełną mocą ciepła powyłączeniowego przez co najmniej 24 godziny bez dostarczania dodatkowej wody z zewnątrz (dotyczy różnych klas zdarzeń projektowych). Standard ANS-5.1-2014 „Decay Heat Power in Light Water Reactors" jest podstawą obliczeń w USA; norma IAEA SSG-2 obejmuje podobne wymagania.

Z punktu widzenia dydaktycznego: wzór Way-Wignera jest przybliżeniem empirycznym o dokładności ~20–30%. Dokładniejsze obliczenia wymagają kodu ORIGEN lub sumowania wkładów setek nuklidów. Dokładność ma znaczenie dla oceny bezpieczeństwa i niepewności obliczonych temperatur rdzenia.

Trzy przykłady rachunkowe

Przykład 1 — aktywność I-131 po kilku dniach chłodzenia

Aktywność I-131 w chwili t=0 (wyłączenie reaktora): A₀ = 10¹⁴ Bq na tonę ciężkiego metalu (typowa wartość przy wypaleniu ~40 GWd/tHM, szacunek z kodu ORIGEN). T₁/₂ = 8,02 dnia, λ = ln(2)/8,02 = 0,0864 d⁻¹.

Po t = 8 dniach:
A(8) = 10¹⁴ × e^(−0,0864 × 8) = 10¹⁴ × e^(−0,691) ≈ 10¹⁴ × 0,501 = 5,0 × 10¹³ Bq/tHM

Po t = 80 dniach:
A(80) = 10¹⁴ × e^(−0,0864 × 80) = 10¹⁴ × e^(−6,91) ≈ 10¹⁴ × 0,001 = 1,0 × 10¹¹ Bq/tHM

Po t = 160 dniach:
A(160) ≈ 10¹⁴ × 10⁻⁶ = 10⁸ Bq/tHM — o 6 rzędów wielkości mniej niż w chwili wyłączenia.

Wniosek: I-131 jest ważny w ciągu pierwszych kilku tygodni, po 5–6 miesiącach jego wkład jest marginalny. W tym czasie Cs-137 (T₁/₂ = 30 lat) traci zaledwie 0,004% swojej aktywności.

Przykład 2 — równowaga Ba-140/La-140

Ba-140 (T₁/₂ = 12,75 d, λ₁ = 0,0544 d⁻¹) produkuje La-140 (T₁/₂ = 1,68 d, λ₂ = 0,413 d⁻¹).

W chwili t = 0 izolujemy czyste Ba-140 z aktywnością A₀(Ba) = 1000 Bq. La-140 = 0.

Aktywność La-140 jako funkcja czasu (równanie Batemana):
A_La(t) = A₀(Ba) × [λ₂/(λ₂−λ₁)] × (e^(−λ₁t) − e^(−λ₂t))

Podstawiamy: λ₂/(λ₂−λ₁) = 0,413/(0,413−0,0544) = 0,413/0,359 = 1,150

Po t = 3 d:
A_La(3) = 1000 × 1,150 × (e^(−0,163) − e^(−1,239)) = 1000 × 1,150 × (0,850 − 0,290) = 1000 × 1,150 × 0,560 = 644 Bq

Maksimum La-140 wyznaczamy z warunku dA/dt = 0 → t_max = ln(λ₂/λ₁)/(λ₂−λ₁) = ln(7,59)/0,359 = 2,03/0,359 ≈ 5,65 d. To maksimum La-140 wypada po ok. 5,6 dniach od izolacji Ba-140 (ok. 0,65 Ba-140 nadal jest czynne).

Wniosek: La-140 osiąga maksimum po kilku dniach, nie natychmiast. To jest klasyczny przykład równowagi przejściowej — przydatny do zrozumienia wielu innych par rodzic-córka w wypalonym paliwie.

Przykład 3 — szacunek ciepła powyłączeniowego

Ciepło powyłączeniowe po wyłączeniu reaktora (t = 0) maleje przybliżoną regułą Way-Wignera:

P(t) ≈ 0,066 × P₀ × t⁻⁰·² (jednostki: t w sekundach, P₀ i P(t) w niespójnych jednostkach, tylko proporcja)

Precyzyjniej korzysta się z kodu ORIGEN lub danych ANS-5.1. Dla uproszczenia edukacyjnego:

Reaktor o mocy termicznej P₀ = 3000 MW_th (typowy PWR). Po t = 1 s:
P(1) ≈ 0,066 × 3000 = 198 MW_th

Po t = 1 min (60 s):
P(60) ≈ 0,066 × 3000 × 60⁻⁰·² = 198 × 60⁻⁰·² = 198 × 0,569 ≈ 113 MW_th

Po t = 1 h (3600 s):
P(3600) ≈ 198 × 3600⁻⁰·² = 198 × 0,288 ≈ 57 MW_th

Po t = 1 d (86400 s):
P(86400) ≈ 198 × 86400⁻⁰·² = 198 × 0,181 ≈ 36 MW_th

Po t = 7 d (604800 s):
P(604800) ≈ 198 × 604800⁻⁰·² = 198 × 0,127 ≈ 25 MW_th

Wniosek: godzinę po wyłączeniu reaktor nadal produkuje 57 MW ciepła — niemal 2% mocy nominalnej — wyłącznie z rozpadu promieniotwórczego. To uzasadnia, dlaczego aktywne chłodzenie jest konieczne przez wiele godzin i dni po wyłączeniu. Brak chłodzenia prowadzi do stopienia paliwa, co było przyczyną wypadku w Fukushimie Daiichi (utrata zasilania pomp chłodzących).

Pytania otwarte

1. Dlaczego sam czas połowicznego zaniku nuklidu nie wystarcza do oceny, czy jest on „świeżym" wskaźnikiem rozszczepienia? Jakie inne czynniki (historia napromieniania, równowagi rodzic-córka, transport chemiczny) komplikują interpretację?

2. Np-239 ma T₁/₂ = 2,35 dnia. Ile procent wyjściowej aktywności Np-239 pozostaje po 14 dniach, po 30 dniach i po 100 dniach? Co te wartości mówią o oknie czasowym, w którym Np-239 jest użytecznym markerem niedawnego wychwytu neutronów?

3. Porównaj trzy nuklidy: I-131 (T₁/₂ = 8 d), Ba-140 (T₁/₂ = 12,75 d) i Cs-137 (T₁/₂ = 30 lat). Jak wyglądałoby widmo gamma próbki tuż po wyłączeniu reaktora vs. 90 dni po wyłączeniu vs. 1 rok po wyłączeniu? Który nuklid dominowałby aktywność w każdym przypadku?

4. System CTBT używa stosunku aktywności Xe-133m/Xe-133 jako znacznika czasu od wybuchu. Opisz jakościowo, dlaczego ten stosunek maleje wraz z upływem czasu od zdarzenia.

5. Równowaga wiekowa definiuje stan, gdy aktywność córki równa się aktywności rodzica. Jakie warunki muszą być spełnione (proporcja T₁/₂, czas od izolacji)? Kiedy równowaga wiekowa nigdy nie zostaje osiągnięta?

6. Wartość wypalenia paliwa wyznaczona przez Nd-148 vs. stężenie Cs-137 powinna dawać spójny wynik. Jakie systematyczne odchylenia mogą wystąpić, jeśli część Cs-137 wypłynęła z paliwa przez pęknięty pręt paliwowy? Jak safeguards MAEA mogłyby to wykryć?

7. Ciepło powyłączeniowe maleje według reguły Way-Wignera w przybliżeniu jak t^(−0,2). Dlaczego zależność jest potęgowa, a nie wykładnicza? Jakie właściwości rozkładu czasów połowicznego zaniku produktów rozszczepienia powodują takie zachowanie zbiorcze?

8. Kalkulator eduakcyjny modeluje inwentarz z uproszczonym zestawem nuklidów. Jakie klasy błędu systematycznego może to wprowadzić w porównaniu z pełnym kodem ORIGEN? Jak student może ocenić, czy uproszczony model wystarczy do jego celu?

Podsumowanie dydaktyczne

1. Świeże wypalone paliwo jest dynamiczną mieszaniną setek radionuklidów o różnych czasach połowicznego zaniku. Pojęcie „radioaktywności paliwa" bez podania czasu chłodzenia jest niejednoznaczne — aktywność i skład radionuklidów zmieniają się o rzędy wielkości w ciągu dni, miesięcy i lat.

2. Np-239 (T₁/₂ = 2,35 d) jest krótkotrwałym ogniwem łańcucha U-238 → U-239 → Np-239 → Pu-239. Jego obecność w próbce wskazuje na niedawny wychwyt neutronów w U-238. Po kilku tygodniach Np-239 jest nieobecny, a Pu-239 pozostaje jako produkt długowieczny.

3. I-131 (T₁/₂ = 8 d) jest zegarem tygodniowym — istotnym zagrożeniem radiologicznym w pierwszych tygodniach, ale śladowym składnikiem po pół roku. Jego biologiczne znaczenie (wychwyt przez tarczycę) uzasadnia szczególne zainteresowanie mimo relatywnie krótkiego T₁/₂.

4. Równania Batemana opisują produkcję i rozpad par rodzic-córka. Równowaga przejściowa (Ba-140/La-140) i wiekowa (Cs-137/Ba-137m) są przykładami, gdzie krótkotrwała córka może być stale zasilana przez dłużej żyjącego rodzica — jej obecność nie musi wskazywać na świeże rozszczepienie.

5. Kody obliczeniowe (ORIGEN, FISPACT) są niezbędne do ilościowej analizy pełnego inwentarza. Proste modele jednomuklidowe dają intuicję fizyczną, ale nie zastąpią kodu przy obliczeniach inżynierskich lub porównaniach z danymi pomiarowymi.

6. System monitoringu CTBT używa ksenonu (Xe-133, Xe-133m, Xe-135, Xe-131m) jako markera jądrowego — ksenon jest gazem szlachetnym, nie deponuje się, i szybko migruje w atmosferze. Stosunki aktywności tych izotopów pozwalają datować zdarzenie i odróżnić reaktorowe tło od sygnału testu.

7. Ciepło powyłączeniowe po wyłączeniu reaktora maleje według przybliżonej reguły Way-Wignera. Godzinę po wyłączeniu reaktor PWR o mocy 3000 MW nadal oddaje ~57 MW tylko z rozpadów promieniotwórczych. Brak chłodzenia jest przyczyną wypadków takich jak Fukushima Daiichi.

8. Kluczowa metareguła: interpretacja widma gamma świeżego paliwa wymaga uwzględnienia czasu od wyłączenia (lub ekspozycji), historii mocy, równowag rodzic-córka, geometrii, wydajności detektora i korekty rozpadu podczas pomiaru. Pominięcie któregokolwiek z tych czynników prowadzi do systematycznych błędów w wyznaczanej aktywności.

Dodatkowe materiały multimedialne

Najkrótsze podsumowanie: świeże wypalone paliwo jest zegarem wielowskaźnikowym. Np-239 mówi o niedawnym wychwycie neutronów w U-238, I-131 o skali tygodni, a dłużej żyjące produkty rozszczepienia i aktynowce przejmują znaczenie po miesiącach i latach.

Powiązane kalkulatory i narzędzia

• Hodowla plutonu-239 — pokazuje profil U-239, Np-239, Pu-239 i Pu-240 oraz chłodzenie po napromienianiu.
• Burnup — pokazuje wypalanie paliwa, ciepło powyłączeniowe i uproszczony inwentarz produktów rozszczepienia.
• Inwentarz odpadów — pokazuje grupowy bilans krótkotrwałych produktów rozszczepienia, Cs/Sr, aktynowców i produktów aktywacji po chłodzeniu.
• Aktywność odpadów jądrowych — pokazuje aktywność, ciepło i dominujące nuklidy w uproszczonym inwentarzu odpadów.
• Kalkulator: Równowaga szeregu — Dwuczłonowy model Batemana: rodzic → córka i narastanie aktywności potomnej.

Ćwiczenia praktyczne

Pierwsze ćwiczenie rachunkowe: policzyć, jaka część aktywności nuklidu o T1/2 = 8 dni pozostaje po 8, 16, 40 i 80 dniach. Następnie porównać z nuklidem o T1/2 = 30 lat w tym samym okresie.

Drugie ćwiczenie: dla łańcucha U-239 -> Np-239 -> Pu-239 narysować jakościowy wykres po zakończeniu napromieniania. Student ma zaznaczyć, dlaczego Np-239 jest sygnałem krótkiego czasu, a Pu-239 produktem długowiecznym.

Trzecie ćwiczenie: użyć syntetycznej tabeli aktywności I-131, Ba-140, Zr-95, Cs-137 i Sr-90 dla trzech czasów chłodzenia. Student ma wskazać, która próbka jest najświeższa i które nuklidy o tym decydują.

Czwarte ćwiczenie: przeanalizować parę rodzic-córka, np. Cs-137/Ba-137m albo Ba-140/La-140. Student ma wyjaśnić, dlaczego krótki czas połowicznego zaniku córki nie zawsze oznacza świeże rozszczepienie.

Piąte ćwiczenie audytowe: przeczytać opis "widma świeżego paliwa" i sprawdzić, czy podano czas od wyłączenia, historię mocy, wypalenie, typ detektora, czas pomiaru, korektę rozpadu podczas pomiaru oraz model równowag rodzic-córka.

Szóste ćwiczenie modelowe: w arkuszu kalkulacyjnym porównać trzy krzywe A(t) dla czasów połowicznego zaniku 1 dzień, 8 dni i 30 lat. Wniosek ma opisać, dlaczego świeżość materiału najlepiej czyta się z kilku nuklidów naraz.

Przejdź do ćwiczenia interaktywnego

Powiązane artykuły

• Uran-238: rola w reaktorze i bombie
• Pluton-239: produkcja i właściwości
• Pluton-240: izotop trujący
• Głębokość wypalania paliwa (Burnup)
• Proces PUREX
• Paliwo MOX (mixed oxide)