Streszczenie

Aktywność źródła nie jest tym samym co liczba zliczeń w detektorze. Detektor widzi tylko tę część promieniowania, która trafia w jego objętość czynną, ma odpowiednią energię, nie została pochłonięta po drodze i zostanie zarejestrowana przez tor pomiarowy. Dlatego bez geometrii, wydajności i poprawek liczba impulsów nie daje jeszcze aktywności bezwzględnej.1,2

Ten artykuł wyjaśnia, dlaczego kąt bryłowy, odległość, rozmiar źródła, samopochłanianie, matryca próbki i geometria naczynia są częścią wyniku. To fundament dla spektrometrii gamma, pomiarów żywności i materiałów budowlanych, ćwiczeń z licznikami GM oraz każdego raportu, który z CPS albo pola fotopiku chce przejść do Bq, Bq/kg lub Bq/l.2,3

Rozszerzenie tematu

Detektor nie widzi całej emisji

Źródło promieniotwórcze może emitować promieniowanie w wielu kierunkach. Detektor zajmuje tylko pewien fragment przestrzeni. Jeśli źródło traktujemy jako punktowe, a emisję jako izotropową, intuicyjnie interesuje nas ułamek pełnego kąta bryłowego 4π, który „zasłania” detektor. Ten ułamek jest geometryczną częścią wydajności.

W prostym przybliżeniu dalekiego pola częstość zliczeń spada jak odwrotność kwadratu odległości. Ale to przybliżenie działa tylko wtedy, gdy źródło jest małe względem odległości, detektor też można potraktować prosto, a oś układu jest dobrze zdefiniowana. W wielu pomiarach laboratoryjnych warunki te nie są spełnione. Próbka ma objętość, detektor ma powierzchnię i grubość, a naczynie pomiarowe ma kształt.

Dziunikowski i Kalita poświęcają temu osobne ćwiczenie: wyznaczaniu poprawki na rozmiary źródła i geometrię pomiaru.1 To bardzo dobra wskazówka dydaktyczna. Geometria nie jest „ustawieniem do zdjęcia”, tylko wielkością metrologiczną.

Aktywność, emisja i zliczenia

Uproszczony łańcuch wygląda tak:

aktywność -> liczba rozpadów -> emisja danej linii -> promieniowanie dochodzące do detektora -> oddziaływanie w detektorze -> impuls -> zliczenie.

Każdy etap ma własny współczynnik. Dla pojedynczej linii gamma można myśleć o zależności:

R = A * I_gamma * epsilon_geom * epsilon_int * f_abs * f_dead * ...,

gdzie R to częstość zliczeń netto, A aktywność, I_gamma prawdopodobieństwo emisji danej linii, epsilon_geom część geometryczna, epsilon_int wewnętrzna wydajność detektora, f_abs poprawki absorpcji i samopochłaniania, a f_dead poprawka czasu martwego.

W praktyce wiele składników łączy się w jedną krzywą wydajności epsilon(E) dla konkretnej geometrii. To wygodne, ale nie wolno zapominać, że ta krzywa obowiązuje tylko dla warunków, w których została wyznaczona.

Kąt bryłowy

Kąt bryłowy jest trójwymiarowym odpowiednikiem kąta płaskiego. Pełna przestrzeń wokół punktu ma 4π sr. Detektor widziany ze źródła zajmuje tylko część tej przestrzeni. Im bliżej detektora znajduje się źródło i im większa powierzchnia czynna detektora, tym większy ułamek emisji może zostać zarejestrowany.

W idealizowanym układzie punktowego źródła na osi małego detektora łatwo narysować stożek. W rzeczywistym układzie trzeba uwzględnić:

  • rozmiar detektora,
  • rozmiar źródła,
  • odległość,
  • przesunięcie poza oś,
  • osłony i kolimatory,
  • martwe warstwy i obudowę detektora,
  • objętość próbki.

To właśnie dlatego pomiar punktowego źródła kalibracyjnego i dużej próbki środowiskowej w słoiku nie są automatycznie porównywalne.

Źródło punktowe, tarczowe i objętościowe

Źródło punktowe to wygodne przybliżenie, ale wiele próbek nim nie jest. Cienka plamka aktywności na podłożu, dysk, filtr aerozolowy, proszek w pudełku, ciecz w naczyniu Marinelli i próbka gleby to różne geometrie.

Źródło tarczowe ma punkty położone w różnych odległościach od detektora i pod różnymi kątami. Źródło objętościowe ma dodatkowo warstwy znajdujące się głębiej w próbce. Promieniowanie z głębszych warstw może być częściowo pochłaniane przez sam materiał próbki, zanim dotrze do detektora.

Dziunikowski i Kalita opisują podejście, w którym wyznacza się zależność częstości zliczeń od położenia punktowego źródła, a potem używa jej do oceny poprawki dla źródła o skończonych rozmiarach.1 Dydaktycznie najważniejszy wniosek brzmi: duże źródło nie jest punktem powiększonym skalą aktywności. Jego geometria zmienia odpowiedź układu.

Naczynie Marinelli i geometria powtarzalna

W spektrometrii gamma próbek środowiskowych często używa się naczyń o powtarzalnej geometrii, na przykład naczyń Marinelli. Ich sens polega na tym, aby próbka otaczała detektor w kontrolowany sposób i aby wzorce oraz próbki były mierzone w możliwie tej samej konfiguracji.

KChRS wykorzystuje w ćwiczeniach pomiary w naczyniach o określonej geometrii, gdy przechodzi od zliczeń fotopików do aktywności w żywności lub materiałach budowlanych.2 To nie jest detal organizacyjny. Jeśli wzorzec był w innej geometrii niż próbka, wydajność detektora może być inna nawet dla tej samej energii.

Najlepszą praktyczną poprawką jest powtarzalność: taka sama objętość, taka sama masa lub gęstość w zakresie kalibracji, taki sam pojemnik, takie samo położenie i ta sama procedura tła. Im mniej zmieniamy między wzorcem a próbką, tym mniej trzeba „ratować” obliczeniami.

Samopochłanianie

Samopochłanianie oznacza, że część promieniowania zostaje pochłonięta w samej próbce. Dla gamma efekt zależy od energii fotonów, składu i gęstości próbki oraz drogi przez materiał. Dla alfa i beta problem jest jeszcze ostrzejszy, bo zasięgi są krótkie i geometria źródła może całkowicie zdominować wynik.

KChRS wskazuje, że poprawka na samopochłanianie może być pomijalna, jeśli kalibracja wydajnościowa była wykonana dla próbek w tej samej konfiguracji geometrycznej, o podobnej gęstości i współczynniku osłabienia.2 To bardzo praktyczne zdanie: nie chodzi o to, żeby każdą poprawkę liczyć od zera, lecz żeby wiedzieć, kiedy wzorzec i próbka są do siebie dostatecznie podobne.

W przeciwnym razie trzeba uwzględnić matrycę. Próbka soli potasowej, gleby, betonu, wody, popiołu i rudy nie osłabia promieniowania tak samo. Jeśli badamy aktywność właściwą w Bq/kg, masa próbki jest potrzebna, ale nie wystarcza. Liczy się też rozmieszczenie aktywności w objętości.

Wydajność wewnętrzna i geometryczna

Wydajność geometryczna mówi, jaka część emitowanego promieniowania dociera do detektora. Wydajność wewnętrzna mówi, jaka część promieniowania, które dotarło, faktycznie wytworzy zarejestrowany sygnał. Dla gamma zależy to od energii i materiału detektora. Dla GM detekcja gamma może mieć małą wydajność, bo foton musi oddziałać w gazie lub w ściance tak, aby powstał elektron uruchamiający impuls.1

W spektrometrii gamma trzeba jeszcze rozróżnić wydajność całkowitą od wydajności w fotopiku. Foton może oddziałać w detektorze, ale oddać tylko część energii i trafić do kontinuum Comptona zamiast do piku pełnej absorpcji. Jeśli liczymy aktywność z fotopiku, interesuje nas wydajność fotopiku, a nie każdy zarejestrowany impuls.

To jest powód, dla którego widmo gamma zawiera tak dużo informacji, ale wymaga ostrożnej kalibracji. Pole piku jest użyteczne dopiero po powiązaniu z wydajnością dla tej energii i geometrii.

Odległość nie zawsze rozwiązuje problem

Można pomyśleć: wystarczy odsunąć źródło daleko, wtedy wszystko stanie się punktowe. Czasem to pomaga, ale ma cenę. Większa odległość zmniejsza częstość zliczeń, wydłuża potrzebny czas pomiaru i pogarsza statystykę. Dla słabych próbek środowiskowych odsuwanie detektora może być niepraktyczne.

Z kolei bardzo bliska geometria daje dużo zliczeń, ale jest bardziej wrażliwa na milimetrowe przesunięcia i rozmiar próbki. Mały błąd położenia może wtedy mocno zmienić kąt bryłowy. Właśnie dlatego powtarzalne uchwyty, pojemniki i procedury ustawiania są tak ważne.

Dobór geometrii jest więc kompromisem:

  • blisko: więcej zliczeń, większa czułość na ustawienie,
  • daleko: prostsza geometria, mniej zliczeń,
  • objętościowo: dobra wydajność dla słabych próbek, większy problem samopochłaniania i matrycy,
  • kolimator: lepsza definicja kierunku, mniejsza wydajność.

Kalibracja względna

W wielu praktycznych pomiarach nie trzeba wyznaczać wszystkich składników osobno. Jeśli wzorzec i próbka są mierzone w tej samej geometrii, stosunek aktywności może być wyznaczony ze stosunku częstości zliczeń po uwzględnieniu tła i rozpadu. Dziunikowski i Kalita opisują takie podejście: przy tej samej geometrii część poprawek skraca się w ilorazie.1

To jest ogromnie użyteczne. Kalibracja względna mówi jednak tylko tyle, ile podobieństwo warunków. Jeśli próbka jest grubsza, gęstsza, przesunięta, ma inną matrycę albo inną energię promieniowania, skrócenie poprawek przestaje być pełne.

W artykułach serwisu warto więc pisać nie tylko „zmierzono względem wzorca”, ale także: czy wzorzec i próbka miały tę samą geometrię i matrycę, czy zastosowano korektę wydajności, czy wynik dotyczy aktywności całkowitej, właściwej, powierzchniowej czy objętościowej.

Geometria alfa i beta

Dla promieniowania alfa geometria jest bezlitosna. Cząstki alfa mają bardzo mały zasięg w powietrzu i materii, więc grubość źródła, nierówność warstwy, odległość od detektora i okienko mogą zdecydować o kształcie widma. Dlatego spektrometria alfa wymaga cienkich, jednorodnych źródeł i dobrze powtarzalnej geometrii.2

Dla beta problemem jest rozpraszanie, samopochłanianie i widmo ciągłe energii. Detektor GM może zliczać cząstki beta, ale przejście od zliczeń do aktywności wymaga znajomości geometrii, okienka, absorpcji i energii. Prosty pomiar bez poprawek jest wskaźnikiem, nie pełnym oznaczeniem.

Dla gamma geometria bywa bardziej „wybaczająca”, bo promieniowanie jest przenikliwe, ale nie znika. Przy niskich energiach i dużych próbkach samopochłanianie może być bardzo ważne.

Co zapisać w raporcie

Jeśli wynik ma być ilościowy, opis geometrii powinien zawierać:

  • typ detektora i jego powierzchnię lub geometrię wejściową,
  • odległość i położenie próbki względem detektora,
  • kształt i rozmiar próbki,
  • masę, objętość i gęstość, jeśli są istotne,
  • typ pojemnika,
  • informację o kolimatorze lub osłonie,
  • geometrię wzorca,
  • zgodność wzorca i próbki,
  • poprawkę samopochłaniania lub uzasadnienie jej pominięcia,
  • krzywą wydajności używaną dla tej geometrii.

Bez tych danych aktywność bezwzględna może wyglądać precyzyjnie, ale nie będzie dobrze udokumentowana. Sama liczba Bq/kg nie mówi, czy porównanie z innym laboratorium ma sens.

Analityczny kąt bryłowy: od wzoru do intuicji

Dla punktowego źródła S na osi symetrii okrągłego detektora o promieniu r, w odległości d, kąt bryłowy można wyrazić wzorem analitycznym:

Ω = 2π × [1 − d / √(d² + r²)]

Ułamek kąta bryłowego: ε_geom = Ω / (4π) = (1/2) × [1 − d / √(d² + r²)].

Dla kilku wartości odległości (zakładając r = 3 cm, typowy detektor NaI 3"):

Odległość d [cm] ε_geom [%]
1 cm 31,4%
3 cm 14,6%
5 cm 7,3%
10 cm 2,2%
20 cm 0,57%
50 cm 0,09%

Widać, że dla odległości 5× promień detektora (d = 15 cm przy r = 3 cm) ułamek kąta bryłowego wynosi zaledwie ~0,3%, co odpowiada około 0,3% wszystkich emitowanych fotonów trafiających w area detektora. Z tego powodu źródła o małej aktywności wymagają pomiarów w bliskiej geometrii (d~1–3 cm) lub długich czasów pomiaru.

Wzór zakłada źródło punktowe i płaski detektor bez efektów brzegowych. Dla detektorów HPGe o walcowatej geometrii „wejście" efektywnego detektora jest bardziej złożone, bo foton musi dotrzeć nie tylko do powierzchni, ale i do strefy aktywnej (po przejściu przez nieaktywną warstwę wierzchnią, dead layer).

Korekcja na martwą warstwę detektora (dead layer)

Nowoczesne detektory HPGe mają cienką nieaktywną warstwę zewnętrzną, gdzie pochłonięte fotony nie generują zarejestrowanego sygnału. Grubość tej warstwy (dead layer) wynosi typowo 0,3–1,0 mm dla HPGe p-type i 0,3 mm dla HPGe n-type.

Dla fotonów o niskich energiach (poniżej ~100 keV) martwa warstwa ma istotny wpływ na wydajność rejestracji — foton jest pochłaniany przed dotarciem do strefy aktywnej. Dla wyższych energii (>200 keV) fotony przechodzą przez martwą warstwę z wysokim prawdopodobieństwem.

Korekcja na martwą warstwę jest szczególnie ważna przy:

  • Pomiarach nuklidów emitujących nisko-energetyczne promieniowanie X lub gamma (Pb-210, Am-241)
  • Kalibracji wydajności w zakresie <100 keV
  • Obliczeniach Monte Carlo geometrii detektora

Producenci dostarczają specyfikacje grubości martwej warstwy dla konkretnego detektora, ale warstwa może rosnąć ze starzeniem kryształu i wymagać ponownej charakteryzacji.

Krzywa wydajności i jej kalibracja

Wydajność fotopiku ε_peak(E) jest funkcją energii promieniowania gamma. Dla typowych detektorów HPGe i NaI(Tl) ma charakterystyczny kształt:

  • Wzrost od ~50 keV do maksimum (~100–200 keV)
  • Plateau lub łagodny spadek dla wyższych energii
  • Monotoniczny spadek powyżej ~200 keV dla HPGe (spada z energią), choć wolniejszy niż dla NaI

Kalibracja krzywej wydajności wymaga zestawu źródeł certyfikowanych o znanych aktywnościach i różnych energiach gamma. Standardowe nuklidy kalibracyjne:

Nukleid Linia gamma [keV] T₁/₂ Zastosowanie
Am-241 59,5 432,6 lat Niskie energie
Cd-109 88,0 462 dni Niskie energie
Co-57 122,1 / 136,5 271,8 dni Środkowy zakres
Cr-51 320,1 27,7 dni Środkowy zakres
Cs-137 661,7 30,17 lat Środkowy zakres
Mn-54 834,8 312 dni Wysoki zakres
Co-60 1173,2 / 1332,5 5,27 lat Wysoki zakres
Y-88 898,0 / 1836,1 106,6 dni Bardzo wysoki zakres

Punkty z tych pomiarów dopasowuje się krzywą empiryczną (np. wielomian stopnia 3–4 do ln ε vs. ln E) lub krzywą segmentową. Niepewność kalibracji wydajności (typowo 2–5% dla dobrze wzorcowanej geometrii) stanowi często dominującą składową niepewności w analizie próbek środowiskowych.

Koincydencje sumacyjne (true coincidence summing)

Gdy nukleid emituje w jednym rozpadzie dwa lub więcej fotonów gamma praktycznie jednocześnie (np. Co-60 emituje 1173 keV i 1332 keV w kaskadzie), detektor o bliskiej geometrii może zarejestrować oba fotony jako jeden impuls o energii 1173+1332 = 2505 keV. Ten efekt nazywany jest koincydencją sumacyjną (TCS — True Coincidence Summing).

TCS powoduje:

  • Ubytek zliczeń w obu pikach (1173 keV i 1332 keV)
  • Pojawienie się pozornych pików na wyższych energiach
  • Systematyczne zaniżenie wyznaczonej aktywności, jeżeli nie jest stosowana korekcja

TCS zależy TYLKO od geometrii (im bliżej detektora, tym większy efekt — więcej fotonów z kaskady trafia do detektora jednocześnie) i od schemy rozpadu (liczba i intensywność kaskad gamma). Nie zależy od aktywności źródła.

Korekcja TCS wymaga znajomości pełnego schematu poziomów energetycznych nuklidu oraz wydajności sumy (Σε) dla danej geometrii. Dla dokładnych pomiarów Co-60, Eu-152, Ba-133 i innych nuklidów z kaskadami gamma korekcja TCS jest obowiązkowa przy bliskiej geometrii.

Praktyczna zasada: unikanie TCS jest możliwe przez odsunięcie źródła na odległość ~10–15 cm od detektora, gdzie wydajność sumy jest pomijalnie mała. Kosztem jest znacznie gorsza statystyka.

Korekcja na macierz: absorpcja w samej próbce

Próbka środowiskowa (gleba, żywność, woda, materiał budowlany) nie jest przezroczysta na promieniowanie gamma. Im niższa energia fotonów i im gęstsza próbka, tym więcej promieniowania ulega pochłoniętym w materiału próbki przed dotarciem do detektora.

Korekcja na samopochłanianie próbki (matrix correction) polega na wyznaczeniu współczynnika transmisji T_mat dla danej próbki i energii:

T_mat = (1 − exp(−μ_mat × d_eff)) / (μ_mat × d_eff)

gdzie μ_mat jest masowym współczynnikiem osłabienia materiału próbki (cm²/g), a d_eff = ρ × t to grubość masowa próbki (g/cm²).

W praktyce stosuje się dwie metody:

  1. Wzorcowanie matrycowe: Wzorzec ma skład chemiczny i gęstość podobną do próbki. Efekty samopochłaniania wchodzą w krzywą wydajności wzorca i automatycznie korygują próbkę.
  2. Transmisja źródła zewnętrznego: Przez próbkę przepuszcza się wiązkę z zewnętrznego źródła gamma i mierzy tłumienie. Znając T_mat, oblicza się korekcję.

Dla próbek o niskiej gęstości i wysokich energiach gamma (>500 keV) korekcja matrycowa może być pomijalnie mała (<2%). Dla próbki betonu (ρ ≈ 2,3 g/cm³) i niskich energii (<100 keV) może być istotna.

Monte Carlo jako narzędzie do kalibracji geometrii

Kiedy eksperymentalna kalibracja geometrii jest trudna (brakuje odpowiednich wzorców, próbka ma niestandardowy kształt), można stosować obliczenia Monte Carlo do symulacji wydajności detekcji.

Programy takie jak GEANT4, MCNP6 i EffTran modelują:

  • Geometrię 3D detektora (kryształ, martwa warstwa, obudowa)
  • Kształt i skład próbki
  • Rodzaje oddziaływań fotonów (fotoelektryczne, Compton, parowanie)
  • Deponowaną energię w objętości aktywnej

Na podstawie symulacji oblicza się „oczekiwaną" wydajność detekcji dla danej geometrii i energii. Wyniki są następnie walidowane przez porównanie z mierzalnymi wzorcami o zbliżonej geometrii.

Metoda Monte Carlo jest standardem w laboratoriach referencyjnych (NIST, PTB, NPL) i jest stosowana przez CLOR do wyznaczania wydajności w niestandardowych geometriach dla próbek środowiskowych.

Kluczowa zasada: symulacja MC jest tylko tak dobra, jak model geometryczny i parametry materiałowe. Niedokładna grubość martwej warstwy detektora lub błędna gęstość próbki przekłada się na błąd symulowanej wydajności.

Trzy przykłady obliczeniowe

Przykład 1: Kąt bryłowy i wydajność geometryczna

Detektor NaI(Tl) o promieniu r = 3,81 cm (3") jest umieszczony w odległości d = 10 cm od punktowego źródła Cs-137.

Wydajność geometryczna:
ε_geom = (1/2) × [1 − d/√(d² + r²)] = (1/2) × [1 − 10/√(100 + 14,52)] = (1/2) × [1 − 10/10,70] = (1/2) × [1 − 0,9346] = (1/2) × 0,0654 = 0,0327 = 3,27%.

Wewnętrzna wydajność NaI(Tl) przy 662 keV (dla kryształu 3"×3"): ε_int ≈ 25% (typowa wartość dla fotonów faktycznie trafiających w detektor).

Wydajność fotopiku: ε_peak ≈ ε_geom × ε_int × ε_fotoabs, gdzie ε_fotoabs to udział fotonów oddziałujących fotoelektrycznie vs. Compton. Dla NaI(Tl) przy 662 keV: ε_peak ≈ 3,27% × 25% = 0,82% (szacunkowe, bez uwzględnienia szczegółów przekrojów czynnych).

Wnikliwa kalibracja daje ε_peak (certyfikowana) ≈ 0,8–1,0% dla tej geometrii. W praktyce używa się zmierzonej krzywej wydajności, nie przybliżonego wzoru.

Przykład 2: Wpływ odległości na wyznaczoną aktywność przy błędnej geometrii

Zmierzono 500 impulsów/min w fotopiku Cs-137 przy odległości d = 5 cm. Kalibracja wydajności była wykonana dla d = 10 cm z wydajnością ε = 0,015 (1,5%).

Jeżeli błędnie zastosujemy ε = 0,015 do pomiaru przy d = 5 cm, wyznaczona aktywność:

A_błędna = R / ε = 500/(60 s) / 0,015 = 8,333 / 0,015 ≈ 556 Bq.

Rzeczywista wydajność przy d = 5 cm (szacunek ze wzoru, r = 3,81 cm):
ε_geom(5cm) = (1/2)×[1 − 5/√(25+14,52)] = (1/2)×[1 − 5/6,286] = (1/2)×0,2046 = 0,1023
ε_geom(10cm) = (1/2)×[1 − 10/√(100+14,52)] = (1/2)×0,0654 = 0,0327

Stosunek wydajności: ε(5cm)/ε(10cm) = 0,1023/0,0327 ≈ 3,13.

Rzeczywista aktywność:
A_rzeczywista = R / ε(5cm) = 500/60 / (0,015 × 3,13) = 8,333 / 0,04695 ≈ 177 Bq.

Błąd przy użyciu niewłaściwej wydajności: (556 − 177)/177 × 100% ≈ 214%.

Wniosek: zmiana odległości z 10 cm na 5 cm bez przeliczenia wydajności daje błąd ponad 200%. Geometria nie jest drobiazgiem.

Przykład 3: Korekta na samopochłanianie w próbce gleby

Próbka gleby o gęstości ρ = 1,6 g/cm³ i grubości t = 3 cm mierzona jest przy energii 661 keV. Masowy współczynnik osłabienia dla gleby (SiO₂ dominuje): μ/ρ ≈ 0,0828 cm²/g.

Masowa grubość: d_eff = ρ × t = 1,6 × 3 = 4,8 g/cm².

Współczynnik osłabienia: µ_lin = (μ/ρ) × ρ = 0,0828 × 1,6 = 0,1325 cm⁻¹.

Transmisja (dla promieniowania gamma przechodzącego przez środek próbki, d = t/2 = 1,5 cm):
T = exp(−0,1325 × 1,5) = exp(−0,1988) = 0,8197.

Korekcja samopochłaniania (dla geometrii objętościowej, czynnik całkowania):
f_SA = (1 − exp(−µ × t)) / (µ × t) = (1 − exp(−0,1325×3)) / (0,1325 × 3) = (1−0,672) / 0,3975 = 0,328/0,3975 = 0,825.

Wyznaczona aktywność musi być podzielona przez f_SA, aby skorygować o samopochłanianie:
A_koryg = A_mierz / f_SA = A_mierz / 0,825 = 1,21 × A_mierz.

Wniosek: zaniedbanie samopochłaniania dla 662 keV (Cs-137) w glebie daje zaniżenie aktywności o ~18%. Dla niższych energii (np. 60 keV) efekt byłby wielokrotnie silniejszy.

Pytania otwarte

  1. Jak wyznaczenie wydajności detekcji zmienia się, gdy zamiast punktowego źródła używamy kalibrowanego roztworu ciekłego w naczyniu Marinelli? Jakie nowe czynniki wchodzą do rachunku wydajności?

  2. Dlaczego efekt true coincidence summing (TCS) rośnie z bliskością geometrii, a nie z aktywnością źródła? Jak modelować TCS w obliczeniach Monte Carlo?

  3. Dla jakich energii gamma samopochłanianie w próbkach środowiskowych (gleba, żywność) jest najbardziej istotne? Gdzie przebiega praktyczna granica, poniżej której korekcja matrycowa jest konieczna?

  4. Jak weryfikuje się, że krzywą wydajności wyznaczona metodą MC zgadza się z pomiarami eksperymentalnymi? Jakie tolerancje rozbieżności stosują laboratoria referencyjne?

  5. Jaką rolę odgrywa grubość martwej warstwy detektora HPGe w kalibracji wydajności dla niskich energii (<100 keV) i jak zmienia się ona ze starzeniem detektora?

  6. Jak sprawdzić, czy pomiary żywności (np. Cs-137 w mięsie) wymagają korekcji matrycowej, jeśli żywność ma niską gęstość i wysoką zawartość wody, a Cs-137 emituje 662 keV?

  7. Jak geometria wpływa na pomiary aktywności alfa i beta w porównaniu z gamma? Dlaczego pomiary alfa mają zwykle dużo ostrzejsze wymagania geometryczne niż gamma?

  8. W jaki sposób zastosowanie kolimatorów zmieniając pole widzenia detektora, wpływa na absolutną aktywność wyznaczoną z pomiaru? Czy kolimator zawsze redukuje aktywność pozorną?

Podsumowanie dydaktyczne

  1. Detektor mierzy zdarzenia, nie aktywność — przejście od zliczeń (CPS) do aktywności (Bq) zawsze wymaga co najmniej czterech składników: kąta bryłowego, wewnętrznej wydajności detekcji, intensywności linii gamma i poprawek na absorpcję. Pominięcie któregokolwiek prowadzi do wyniku pozornie liczbowego, ale fizycznie bezsensownego.

  2. Geometria jest parametrem pomiaru — kąt bryłowy, odległość źródło-detektor i rozmiar próbki są mierzalnymi wielkościami, nie tylko „ustawieniami". Zmiana odległości z 10 cm na 5 cm trzykrotnie zwiększa ε_geom; zmiana bez aktualizacji kalibracji wydajności prowadzi do błędu ponad 200%.

  3. Kalibracja względna skraca ścieżkę — jeżeli wzorzec i próbka mają tę samą geometrię i matrycę, wiele poprawek skraca się w ilorazie aktywności. To jest dlatego geometria Marinelli i procedury powtarzalności (te same pojemniki, ta sama objętość) są tak cenne w rutynowej analizie środowiskowej.

  4. Samopochłanianie rośnie z gęstością i maleje z energią — dla niskich energii gamma i gęstych próbek (beton, gleba) korekcja matrycowa może być powyżej 20%. Dla 662 keV i gleby efekt jest rzędu 18%; dla 60 keV — wielokrotnie większy. Energetyczna dependencja samopochłaniania jest kolejnym powodem, dla którego krzywa wydajności musi obejmować odpowiedni zakres energii.

  5. True coincidence summing jest efektem geometrycznym — koincydencje sumacyjne (TCS) dla nuklidów emitujących kaskady gamma rosną z bliskością geometrii, nie z aktywnością. Przy bliskiej geometrii pomiar Co-60, Eu-152 lub Ba-133 bez korekcji TCS daje systematycznie zaniżoną aktywność. Korekcja wymaga znajomości pełnego schematu poziomów jądrowych.

  6. Monte Carlo uzupełnia eksperyment — gdy brakuje certyfikowanych wzorców dla danej geometrii, obliczenia MC (GEANT4, MCNP6, EffTran) dostarczają krzywą wydajności. Jakość wyniku MC zależy od jakości modelu geometrycznego i danych o martwej warstwie detektora — nieznane parametry detektora są dominującym źródłem błędu.

  7. Opis geometrii jest częścią raportu — wynik Bq/kg bez opisu pojemnika, odległości, wzorca, metody kalibracji i poprawek matrycowych nie jest w pełni audytowalny i nie może być bezpośrednio porównany z wynikami z innych laboratoriów. Dokumentacja geometrii jest metrological traceability — element walidujący porównywalność wyników.

  8. Alfa wymaga jeszcze ostrzejszej geometrii — promieniowanie alfa ma zasięg ~4 cm w powietrzu i submilimetrowy w materii. Każda warstwa między źródłem a detektorem (powietrze, folia, cienka warstwa matrycy) silnie tłumi i zniekształca widmo. Dlatego spektrometria alfa wymaga źródeł cienkich (<1 µg/cm²), próżni lub atmosfery inertnej i pomiarów w odległości ~1 cm.

Walidacja geometrii przez pomiar w kilku odległościach

Prostym i eleganckim sposobem weryfikacji, czy dana geometria jest poprawnie skalibrowana, jest pomiar tego samego źródła w kilku odległościach i sprawdzenie, czy wyznaczona aktywność pozostaje zgodna w granicach niepewności.

Jeżeli kalibracja wydajności jest prawidłowa dla każdej z odległości, pomiar źródła certyfikowanego w odległości d₁, d₂ i d₃ powinien dawać ten sam wynik A [Bq] w granicach niepewności kombinowanej (pomiarowej + kalibracyjnej).

Rozbieżność wskazuje na błąd jednej z wydajności — np. zbyt wygładzoną krzywą wydajności w pobliżu detektora, gdzie TCS zaczyna mieć znaczenie, lub efekty objętościowe źródła nieuwzględnione w modelu punktowym.

Praktyczny protokół:

  1. Zmierzyć certyfikowane źródło przy d = 3 cm, 7 cm i 15 cm.
  2. Dla każdej odległości wyznaczył ε_peak z kalibracji lub modelu.
  3. Wyznaczone A₁, A₂, A₃ powinny leżeć w granicach ±2σ od wartości certyfikowanej.
  4. Jeżeli A₁ (bliskie) jest systematycznie niższe od A₂ i A₃, wskazuje to na TCS lub efekty objętościowe przy bliskiej geometrii.
  5. Jeżeli A₃ (dalekie) ma dużą niepewność, pomiar wymaga dłuższego czasu dla lepszej statystyki.

Specyfika próbek ciekłych i naczynia Marinelli

Naczynie Marinelli (beaker Marinelli lub Marinelli geometry) jest naczyniowym standardem stosowanym w spektrometrii gamma próbek środowiskowych. Jego kształt — cylindryczna próbka z wewnętrznym otworem dopasowanym do głowicy detektora — maximalizuje wydajność detekcji przez otoczenie detektora próbką z góry, z boku i od dołu.

Typowe objętości naczyń Marinelli: 0,45 L, 1 L i 2 L. Norma ANSI/IEEE 325-1996 i zalecenienia IAEA definiują standardowe geometrie Marinelli. Wydajności dla tych standardowych geometrii są tabelaryzowane i dostępne od producentów źródeł kalibracyjnych (Eckert & Ziegler, Polatom, Amersham).

Dla naczyń Marinelli kluczowe znaczenie mają:

  • Masa objętościowa próbki (density) — musi być bliska gęstości wzorca kalibracyjnego. Zmiana gęstości o 10% może dawać kilkuprocentową zmianę wydajności dla niskich energii.
  • Wysokość słupa cieczy — musi odpowiadać wzorcowi.
  • Skład chemiczny — szczególnie zawartość pierwiastków o dużych przekrojach czynnych na absorbcję promieniowania.

Technika tzw. „spike & recovery" polega na dodaniu do próbki wzorca certyfikowanego, pomiarze i sprawdzeniu, czy wyznaczona aktywność wzorca odpowiada certyfikowanej. Jest to standardowy test poprawności geometrii i kalibracji.

Znaczenie powtarzalności i dokumentacji w laboratorium

W laboratoriach rutynowej analizy środowiskowej (WIOŚ, CLOR, SANEPID) setki próbek żywności, wody i gleby są mierzone rocznie. Utrzymanie jakości geometrycznych warunków pomiaru wymaga protokołów powtarzalności:

Standardowe naczynie dla każdego typu matrycy: Gleba w pudełku Petri, żywność w naczyniu polietylenowym o stałej objętości, woda w naczyniu Marinelli. Każde naczynie jest ważone przed i po załadowaniu próbki — masa musi mieścić się w zdefiniowanym przedziale (np. 500 ± 10 g).

Regularne testy geometrii: Raz na miesiąc pomiar certyfikowanego źródła w pełnym zestawie geometrii stosowanych w laboratorium. Wynik jest logowany i porównywany z wynikami historycznymi — trendy wskazują na dryfowanie wydajności (starzenie detektora, zmiana grubości martwej warstwy).

Akredytacja laboratorium: Laboratoria analityczne akredytowane przez PCA (Polskie Centrum Akredytacji) muszą wykazywać zdolność do wyznaczania aktywności z niepewnością poniżej określonego limitu (np. ±10% dla Cs-137 w glebie przy aktywności powyżej MDA). Walidacja geometryczna jest częścią tego wymagania.

Śledzenie niepewności kalibracji: Niepewność certyfikowanej aktywności wzorca (np. ±3% dla standardów POLATOM/Eckert & Ziegler) przenosi się na niepewność całą krzywą wydajności. Budżet niepewności kalibracyjnej jest formalnie dokumentowany zgodnie z GUM (Guide to the expression of Uncertainty in Measurement, ISO/IEC Guide 98-3).

Geometria pomiarów w żywności: przykłady z polskich przepisów

Polska ustawa o prawie atomowym i rozporządzenia wykonawcze określają maksymalne dopuszczalne poziomy aktywności (MAD) radionuklidów w żywności i wodzie pitnej. Pomiary zgodności wymagają wiarygodnych oznaczeń aktywności, w tym poprawnie skalibrowanych geometrii.

Typowe geometrie stosowane w polskich laboratoriach monitoringu żywności:

  • Zboże i mąka: 500 g w pojemniku walcowym ∅80 mm × h75 mm, HPGe lub NaI(Tl), pomiar 3600 s.
  • Mleko i napoje ciekłe: 1 L w naczyniu Marinelli, HPGe, pomiar 3600–7200 s.
  • Mięso i ryby: 200 g w pojemniku walcowym, HPGe.
  • Gleba i osady: 150–200 g w pojemniku walcowym, HPGe.

Każda z tych geometrii wymaga oddzielnej krzywej wydajności lub certyfikowanego wzorca w identycznym pojemniku o identycznej gęstości. Brak tej dokumentacji stanowi niezgodność z wymaganiami normy PN-EN ISO/IEC 17025.

Porównanie z innymi technikami: niezależność wyznaczenia aktywności

W spektrometrii ciekłoscyntylacyjnej (LSC) i proporcjonalnej geometria ma inne znaczenie: próbka jest najczęściej mieszana ze scyntylatorem i umieszczana wewnątrz pary fotokatod, co daje prawie 4π geometrię detekcji. Efektywność jest wysoka dla emiterów beta (>90% dla trytu, >95% dla C-14), ale zależy od tzw. quenching — gaszenia scyntylacji przez składniki chemiczne próbki.

Porównanie:

Technika Geometria Typowa wydajność beta Samopochłanianie
GM (beta, zewnętrzny) 2π 30–50% Duże (okienko)
LSC ~4π 50–95% Chemiczne (quench)
HPGe (gamma) ε_peak(E) < 10% N/A Matrycowe (niskie)
NaI(Tl) (gamma) ε_peak(E) 5–30% N/A Matrycowe
Detektor Si (alfa) ~2π N/A Krytyczne (żaden)

Ta tabela pokazuje, że „problem geometrii" jest ogólny, ale manifestuje się różnie w zależności od rodzaju promieniowania i techniki pomiarowej. Każda technika ma swój specyficzny model kalibracji i odpowiednie wymagania na warunki pomiarowe.

Granica wykrywalności a geometria

Granica wykrywalności (MDA — Minimum Detectable Activity lub MDC — Minimum Detectable Concentration) zależy bezpośrednio od geometrii przez wydajność detekcji. Im wyższa wydajność dla danej geometrii (bliżej, większy detektor, lepsze ε_peak), tym niższe MDA można osiągnąć przy tym samym czasie pomiaru.

Wzór Curie-Currie'go na MDA dla spektrometrii gamma:

MDA = (4,65 × √N_bg + 2,71) / (ε_peak × I_γ × t_pom)

gdzie N_bg to liczba zliczeń tła w oknie fotopiku, ε_peak to wydajność fotopiku, I_γ to intensywność emisji linii gamma, a t_pom to czas pomiaru.

Widać wyraźnie, że MDA jest odwrotnie proporcjonalne do ε_peak. Podwojenie wydajności (np. przez zbliżenie próbki o połowę na odległość lub użycie większego detektora) obniża MDA o połowę. Jest to bezpośredni argument za stosowaniem naczyń Marinelli i bliskich geometrii dla próbek o niskiej aktywności.

Z drugiej strony, bliższa geometria zwiększa ryzyko TCS, wymaga dokładniejszej kalibracji matrycowej i jest bardziej wrażliwa na zmianę geometrycznego ułożenia próbki. Dlatego wybór geometrii jest zawsze kompromisem między MDA a jakością/stabilnością kalibracji.

Dla typowego HPGe (∅70 mm × 80 mm) i próbki 1 L w naczyniu Marinelli MDA dla Cs-137 (przy pomiarze 3600 s) wynosi zwykle 0,3–1 Bq/kg w zależności od poziomu tła. Dla pobliżu reaktora jądrowego lub po incydencie ta sama geometria może wykryć setki lub tysiące Bq/kg.

Polskie perspektywy: normalizacja i harmonizacja wyników

W Polsce pomiary aktywności próbek środowiskowych są koordynowane przez CLOR i Państwową Agencję Atomistyki (PAA). Wyniki zgłaszane do EURDEP (sieć europejska) muszą być wyrażone w ustandaryzowanych jednostkach i geometriach, aby były porównywalne z danymi z innych krajów.

CLOR uczestniczy w programach porównawczych (ring tests) organizowanych przez IAEA ALMERA i EC NCRP, gdzie laboratoria z różnych krajów mierzą te same certyfikowane próbki. Wyniki porównania ujawniają rozbieżności między laboratoriami, które często wynikają właśnie z różnic w geometrii pomiaru, stosowanych wzorcach kalibracyjnych lub modelach korekcji matrycowej.

Harmonizacja metod pomiarowych — wybór tych samych pojemników, tych samych wzorców, tych samych protokołów korekcji — jest jednym z priorytetów europejskiej metrologii jądrowej (EURAMET). Jest to praktyczny przykład tego, że „opis geometrii" nie jest biurokratycznym wymogiem, lecz warunkiem porównywalności wyników między niezależnymi laboratoriami.

Polskie wyniki monitoringu aktywności żywności po awarii jądrowej w Czarnobylu (1986) i w związku z incydentem Ru-106 (2017) były raportowane do sieci europejskich właśnie przez standardowe geometrie Marinelli i cylindryczne. Możliwość porównania polskich wyników z wynikami z Finlandii, Szwecji czy Niemiec wynikała bezpośrednio z tej harmonizacji. Geometria pomiaru okazała się fundamentem, który pozwolił skoordynować odpowiedź regulacyjną na skalę kontynentu — nie tylko naukową ciekawostką. To dlatego znajomość zasad kalibracji geometrycznej i poprawek korekcyjnych nie jest wiedzą wyłącznie laboratoryjną, lecz kompetencją niezbędną dla każdego, kto interpretuje wyniki monitoringu radioekologicznego lub podejmuje decyzje ochronne na ich podstawie.

Dodatkowe materiały multimedialne

Warto przygotować interaktywny model geometryczny: punktowe źródło, tarcza i próbka objętościowa nad detektorem. Użytkownik zmieniałby odległość, promień próbki, grubość i współczynnik osłabienia, a model pokazywałby zmianę wydajności geometrycznej i samopochłaniania.

Najkrótsze podsumowanie: detektor nie mierzy aktywności źródła bezpośrednio. Mierzy zdarzenia w konkretnej geometrii, a aktywność pojawia się dopiero po uwzględnieniu kąta bryłowego, wydajności, tła, samopochłaniania i podobieństwa próbki do wzorca.

Powiązane kalkulatory i narzędzia

Ćwiczenia praktyczne

Pierwsze ćwiczenie geometryczne: dla punktowego źródła i okrągłego detektora użyć uproszczonego modelu 1/r^2, aby porównać względną częstość zliczeń przy odległościach 5, 10 i 20 cm. Następnie zapisać, dlaczego ten model przestaje być dobry, gdy źródło jest dużym dyskiem blisko detektora.

Drugie ćwiczenie względne: student dostaje częstość zliczeń wzorca i próbki zmierzonych w tej samej geometrii. Ma policzyć aktywność próbki z proporcji, odjąć tło i wskazać, które poprawki skróciły się dzięki identycznej geometrii.

Trzecie ćwiczenie o samopochłanianiu: porównać dwie próbki o tej samej aktywności, ale różnej gęstości i grubości. Na danych syntetycznych pokazać, że mniejsza częstość zliczeń nie musi oznaczać mniejszej aktywności, lecz większe osłabienie w matrycy.

Czwarte ćwiczenie raportowe: przygotować opis geometrii pomiaru próbki w naczyniu Marinelli. Wypisać masę, objętość, pojemnik, położenie względem detektora, wzorzec użyty do kalibracji i informację, czy próbka oraz wzorzec mają porównywalną gęstość.

Piąte ćwiczenie krytyczne: porównać dwa wyniki Bq/kg z różnych opisów. Jeden zawiera geometrię, kalibrację wydajności i tło, drugi tylko końcową liczbę. Student ma ocenić, który wynik można włączyć do tabeli porównawczej i jakie pytania trzeba zadać autorowi drugiego.

Przejdź do ćwiczenia interaktywnego

Powiązane artykuły