Streszczenie

Kalibracja w laboratorium jądrowym nie jest formalnością wykonywaną przed „właściwym” pomiarem. To część definicji wyniku: bez niej liczba zliczeń, położenie piku albo wskazanie dawki są tylko sygnałem aparatury, a nie wielkością fizyczną możliwą do porównania z innymi laboratoriami.1,2

Niepewność pełni podobną rolę. Nie jest dodatkiem pesymistycznym, lecz informacją o jakości wyniku: ile pochodzi z losowych fluktuacji zliczeń, ile z tła, ile z geometrii, wydajności, masy próbki, kalibracji energii, czasu martwego i stabilności toru. Ten artykuł porządkuje te pojęcia na poziomie potrzebnym do czytania i projektowania ćwiczeń radiometrycznych.1,3

Rozszerzenie tematu

Wynik bez kalibracji jest tylko odczytem

Licznik pokazał 1250 impulsów. Spektrometr pokazał pik w kanale 417. Program wyświetlił 0,18 µSv/h. Każda z tych liczb wygląda konkretnie, ale żadna sama nie jest jeszcze w pełni zinterpretowanym wynikiem. Trzeba wiedzieć, jak przyrząd został skalibrowany, w jakiej geometrii wykonano pomiar, jakie było tło i jaka jest niepewność.

W radiometrii szczególnie łatwo pomylić odczyt z wielkością fizyczną, bo aparatura daje liczby automatycznie. Tymczasem przejście wygląda tak:

zliczenia -> częstość zliczeń -> częstość netto -> poprawki -> aktywność / dawka / stężenie -> niepewność.

Każda strzałka wymaga założenia albo kalibracji. Jeśli któregoś elementu brakuje, wynik może być użyteczny orientacyjnie, ale nie powinien być przedstawiany jako ilościowy pomiar laboratoryjny.

Kalibracja energii

W spektrometrii gamma i alfa analizator wielokanałowy zapisuje impulsy w kanałach. Kanał nie jest energią. Energią staje się dopiero po kalibracji, zwykle wykonanej na pikach znanych radionuklidów. W najprostszym przybliżeniu zależność ma postać:

E = a * k + b,

gdzie E to energia, k numer kanału, a a i b są parametrami dopasowania.

Przewodnik KChRS pokazuje taką logikę w ćwiczeniach spektrometrii alfa: źródło kalibracyjne z kilkoma znanymi energiami pozwala zbudować krzywą kalibracyjną, a następnie przypisywać energię pikom w widmie badanym.2 Ta sama idea wraca w spektrometrii gamma, tylko praktyczne problemy są inne: tło Comptona, nakładanie pików, wydajność zależna od energii i geometria próbki.

Kalibracja energii odpowiada na pytanie: gdzie na osi energii leży pik? Nie odpowiada jeszcze na pytanie: ile radionuklidu jest w próbce. Do tego potrzebna jest kalibracja wydajności.

Kalibracja wydajności

Wydajność detekcji mówi, jaka część emitowanych fotonów lub cząstek zostanie zarejestrowana w interesującym nas piku albo zakresie. To nie jest jedna stała przyrządu. Zależy od:

  • energii promieniowania,
  • typu i rozmiaru detektora,
  • odległości próbki od detektora,
  • kształtu i objętości próbki,
  • osłon i kolimatorów,
  • samopochłaniania w próbce,
  • prawdopodobieństwa emisji danej linii,
  • sposobu wyznaczania pola piku.

Dziunikowski i Kalita zwracają uwagę, że w pomiarach aktywności trzeba uwzględniać geometrię, rozmiary źródła, wydajność detekcji i inne poprawki.1 KChRS pokazuje to praktycznie w ćwiczeniach z żywnością, materiałami budowlanymi i próbkami o znanej aktywności: najpierw mierzy się wzorzec, potem próbkę, a wynik przelicza przez wydajność i tło.2

Uproszczony wzór na aktywność dla pojedynczej linii gamma można zapisać jako:

A = N_net / (t * epsilon * I_gamma * m),

jeśli wynik podajemy jako aktywność właściwą na jednostkę masy. N_net oznacza liczbę zliczeń netto w piku, t czas pomiaru, epsilon wydajność detekcji dla danej energii i geometrii, I_gamma prawdopodobieństwo emisji linii, a m masę próbki. Każdy z tych składników ma własną niepewność.

Tło: wynik, którego nie wolno pominąć

Tło nie jest „szumem do wyzerowania”. Jest mierzoną częścią środowiska i aparatury. W liczniku GM tło odejmujemy od częstości zliczeń próbki. W spektrometrii gamma odejmujemy widmo tła albo lokalne tło pod fotopikiem. W obu przypadkach tło ma własne fluktuacje statystyczne i własną niepewność.

Jeśli pomiar próbki trwa t_s, a pomiar tła t_b, można myśleć o częstościach:

R_s = N_s / t_s,

R_b = N_b / t_b,

R_net = R_s - R_b.

Niepewność R_net nie jest równa niepewności samego pomiaru próbki. Dla prostego przypadku zliczeń Poissona:

u(R_net) = sqrt(N_s / t_s^2 + N_b / t_b^2).

Wniosek praktyczny jest prosty: krótko zmierzone tło może ograniczać jakość całego wyniku. Długi pomiar próbki nie naprawi źle poznanego tła, jeśli odejmowana wartość jest porównywalna z sygnałem.

Statystyka zliczeń: Poisson przed Gaussem

Rozpady promieniotwórcze i rejestracja pojedynczych zdarzeń są dyskretne. Dziunikowski i Kalita omawiają błąd od fluktuacji statystycznych i rozkład Poissona jako podstawowy model liczby zliczeń.1 Najważniejsza intuicja:

u(N) ≈ sqrt(N).

Względna niepewność liczby zliczeń wynosi więc w przybliżeniu:

u(N) / N ≈ 1 / sqrt(N).

Jeśli mamy 100 zliczeń, statystyczna niepewność względna jest rzędu 10%. Jeśli mamy 10 000 zliczeń, spada do około 1%. To tłumaczy, dlaczego w pomiarach małych aktywności czas bywa równie ważny jak czułość przyrządu.

Dla dużych liczebności rozkład Poissona można często przybliżać rozkładem normalnym Gaussa. Dla małych liczebności trzeba uważać: symetryczne przedziały i automatyczne zaokrąglenia mogą dawać złudne poczucie precyzji. Ten temat zasługuje na osobny artykuł o statystyce zliczeń promieniotwórczych.

Propagacja niepewności

Wynik radiometryczny zwykle jest ilorazem i iloczynem kilku wielkości. Jeśli dla uproszczenia przyjmiemy:

A = N_net / (t * epsilon * I_gamma * m),

to względną niepewność można w pierwszym przybliżeniu złożyć kwadratowo:

u_r(A) = sqrt(u_r(N_net)^2 + u_r(t)^2 + u_r(epsilon)^2 + u_r(I_gamma)^2 + u_r(m)^2).

W praktyce czas jest zwykle znany bardzo dobrze, masa próbki może być znana dobrze, a największe wkłady często pochodzą z małej liczby zliczeń, tła, wydajności i geometrii. To jednak nie jest prawo natury. W radiochemii środowiskowej ważny może być odzysk chemiczny. W spektrometrii alfa bardzo ważna może być jakość cienkiego źródła. W pomiarach terenowych dominować może geometria i odpowiedź energetyczna przyrządu.

Najgorszy błąd redakcyjny polega na podaniu tylko niepewności statystycznej, gdy wynik w rzeczywistości ogranicza kalibracja lub przygotowanie próbki. Liczba ±1% wygląda elegancko, ale jeśli wydajność znana jest do 10%, jest po prostu nieuczciwa.

Niepewność typu A i typu B

W praktyce warto rozdzielać dwa źródła informacji o niepewności.

Niepewność typu A pochodzi z analizy danych: powtórzeń, odchylenia standardowego, statystyki zliczeń, rozrzutu wyników. Można ją policzyć z serii pomiarów albo z modelu Poissona.

Niepewność typu B pochodzi z innych informacji: certyfikatu wzorca, tolerancji wagi, dokładności pipety, dokumentacji źródła kalibracyjnego, charakterystyki detektora, stabilności napięcia, geometrii pozycjonowania, doświadczenia z daną metodą.

Laboratorium jądrowe potrzebuje obu. Sama seria powtórzeń pokaże rozrzut, ale nie wykryje źle skalibrowanego wzorca. Certyfikat wzorca nie pokaże, że operator pomylił geometrię albo źle odjął tło. Dlatego sensowny budżet niepewności jest listą składników, a nie jedną liczbą dopisaną na końcu.

Geometria jest częścią kalibracji

Geometria próbka-detektor często decyduje o wyniku. Jeśli wzorzec i próbka mają ten sam kształt, tę samą odległość, podobną gęstość i podobną matrycę, część problemów znika w porównaniu względnym. Jeśli geometria się zmienia, wydajność trzeba znać osobno.

KChRS wyraźnie wskazuje, że poprawka na samopochłanianie może nie być wymagana wtedy, gdy 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 praktyczna zasada: najlepszą poprawką bywa uniknięcie różnicy między wzorcem a próbką.

Ten temat rozwija osobny artykuł o geometrii źródło-detektor. Tutaj najważniejsze jest, żeby traktować geometrię jako parametr pomiaru, a nie ustawienie techniczne bez znaczenia.

Rozdzielczość, dryft i stabilność toru

Kalibracja wykonana rano nie musi być idealna wieczorem. Wzmocnienie fotopowielacza zależy od wysokiego napięcia. Elektronika może dryfować z temperaturą. Detektor półprzewodnikowy może wymagać stabilnych warunków chłodzenia. W spektrometrii gamma przesunięcie pików o kilka kanałów może pogorszyć dopasowanie energii, a w pomiarach aktywności zmiana rozdzielczości może utrudnić wyznaczenie pola piku.

Dlatego kontrola jakości powinna obejmować:

  • okresowe sprawdzanie położenia piku kontrolnego,
  • kontrolę szerokości piku,
  • pomiar tła,
  • pomiar wzorca kontrolnego,
  • zapis ustawień wysokiego napięcia i wzmocnienia,
  • kontrolę czasu martwego,
  • powtarzalność geometrii.

W dobrze prowadzonym laboratorium taka kontrola nie jest biurokracją. To sposób wykrywania, czy wynik dzisiejszy jest porównywalny z wczorajszym.

Błąd aparaturowy a rozrzut statystyczny

Dziunikowski i Kalita zwracają uwagę, że oprócz fluktuacji statystycznych mogą występować przypadkowe zakłócenia aparaturowe i błędy systematyczne.1 To bardzo ważne, bo samo zwiększenie czasu pomiaru zmniejsza tylko część statystyczną. Nie zmniejszy błędnej kalibracji, niestabilnej geometrii ani zaniżonej wydajności.

Przykład: mierzymy próbkę dwa razy dłużej. Liczba zliczeń rośnie, niepewność statystyczna względna spada. Jeśli jednak próbka jest w innym naczyniu niż wzorzec, błąd geometrii pozostaje. Jeśli fotopik jest źle wybrany, dłuższy pomiar tylko precyzyjniej policzy zły obszar.

Dlatego pytanie „ile mamy zliczeń?” jest potrzebne, ale niewystarczające. Trzeba pytać też „czy mierzymy właściwą rzecz w porównywalnych warunkach?”.

Spójność pomiarowa

Kalibracja ma sens dlatego, że łączy wynik z wzorcem. W prostym ćwiczeniu wzorcem może być próbka o znanej aktywności. W profesjonalnym laboratorium potrzebna jest identyfikowalność metrologiczna: wiadomo, skąd pochodzi wzorzec, jaka jest jego niepewność, dla jakiej daty podano aktywność i jak przeliczono rozpady w czasie.

W radiometrii data jest szczególnie ważna. Aktywność wzorca zmienia się zgodnie z okresem półtrwania. Jeśli używamy certyfikatu sprzed kilku lat, trzeba przeliczyć aktywność na dzień pomiaru. Dla długożyciowych radionuklidów zmiana może być mała, dla krótszych już istotna.

To samo dotyczy próbek aktywowanych, generatorów krótkotrwałych i pomiarów środowiskowych po zdarzeniach uwolnienia. Czas pobrania, czas przygotowania i czas pomiaru są częścią metadanych.

Minimalna karta wyniku

Dobry wynik laboratoryjny powinien dać się odtworzyć intelektualnie. Minimalna karta wyniku powinna zawierać:

  • nazwę metody i detektora,
  • datę i czas pomiaru,
  • czas rzeczywisty i czas żywy, jeśli dotyczy,
  • geometrię próbka-detektor,
  • masę lub objętość próbki,
  • tło i sposób jego odjęcia,
  • wzorce użyte do kalibracji,
  • równanie kalibracyjne energii,
  • wydajność lub krzywą wydajności,
  • poprawki na samoabsorpcję, czas martwy i rozpad,
  • wynik z niepewnością,
  • informację, czy niepewność obejmuje tylko statystykę, czy pełniejszy budżet.

Bez takiej karty wynik bywa efektowny, ale trudno go zweryfikować. W tym serwisie warto przyjąć zasadę: jeśli artykuł pokazuje liczbę z pomiaru, powinien też mówić, skąd bierze się jej niepewność.

GUM: ustandaryzowane podejście do niepewności pomiarów

ISO/IEC Guide 98-3:2008 — znany jako GUM (Guide to the expression of Uncertainty in Measurement) — jest fundamentalnym dokumentem definiującym sposób wyznaczania i wyrażania niepewności pomiarowych. Choć powstał jako ogólny standard metryczny, jest stosowany w laboratoriach radiometrycznych na całym świecie, w tym przez CLOR w Polsce.

GUM definiuje:

  • Standardową niepewność u(x_i): odchylenie standardowe rozkładu prawdopodobieństwa dla wielkości x_i.
  • Złożoną niepewność standardową u_c(y): propagacja niepewności z wszystkich wejściowych wielkości x_i na wynik y = f(x_1, x_2, ..., x_N).
  • Rozszerzona niepewność U = k · u_c: stosowana do deklarowania przedziałów ufności. Dla normalnego rozkładu k = 2 odpowiada ~95% przedziałowi ufności.

Wzór propagacji niepewności z GUM (prawo propagacji niepewności — LPU):

u_c²(y) = Σ [∂f/∂x_i]² · u²(x_i) + 2 · Σ Σ [∂f/∂x_i][∂f/∂x_j] · u(x_i, x_j)

Dla nieskorelowanych wielkości wejściowych człon z u(x_i, x_j) = 0 i wzór redukuje się do sumy kwadratów pochodnych cząstkowych (sensitivities) razy kwadratów niepewności składowych.

W praktyce radiometrii GUM wymaga:

  1. Zapisania modelu pomiarowego: y = f(x_1, ..., x_N)
  2. Identyfikacji wszystkich wielkości wejściowych x_i
  3. Wyznaczenia u(x_i) dla każdej składowej (A lub B)
  4. Obliczenia pochodnych cząstkowych c_i = ∂f/∂x_i (tzw. współczynniki czułości)
  5. Złożenia: u_c²(y) = Σ c_i² · u²(x_i)
  6. Wyznaczenia stopni swobody ν_eff (przybliżenie Welcha-Satterthwaite'a) dla właściwego k

ISO/IEC 17025: wymagania akredytacyjne dla laboratoriów

Norma PN-EN ISO/IEC 17025 definiuje wymagania dla laboratoriów badawczych i wzorcujących, które chcą wykazać kompetencje i dostarczyć wiarygodne wyniki. W Polsce akredytację laboratoriów przyznaje PCA (Polskie Centrum Akredytacji), a laboratoria radiometryczne akredytowane przez PCA muszą spełniać wymagania tej normy.

Kluczowe wymagania ISO/IEC 17025 istotne dla laboratoriów jądrowych:

  • Walidacja metod: każda stosowana metoda musi być zwalidowana — jej czułość, specyficzność, liniowość, powtarzalność i granica wykrywalności muszą być określone eksperymentalnie.
  • Nadzór nad sprzętem: kalibracja detektorów, wag, pipet, timer ów i mierników napięcia musi być dokumentowana z określoną częstotliwością. Sprzęt musi mieć identyfikatory i metadane kalibracyjne.
  • Spójność pomiarowa (traceability): wyniki muszą być powiązane z krajowymi lub międzynarodowymi wzorcami miary przez udokumentowany łańcuch porównań.
  • Budżet niepewności: dla każdego wydawanego wynika ilościowego musi być dostępny szacunek niepewności zgodny z GUM.
  • Kontrola jakości: regularne pomiary kontrolne (QC samples), certyfikowanych materiałów referencyjnych i uczestnictwo w testach biegłości (proficiency testing).

Certyfikowane materiały referencyjne (CRM) i ich rola

Certyfikowany materiał referencyjny (CRM) to materiał o prawnie ustalonych właściwościach, wyznaczonych z niepewnościami i udokumentowanym łańcuchem spójności pomiarowej. W radiometrii CRM to np. certyfikowany roztwór Cs-137 o aktywności (10 000 ± 200) Bq/mL na dzień odniesienia.

Producenci CRM dla radiometrii: POLATOM (Polska), Eckert & Ziegler (USA/Niemcy), Amersham Biosciences (UK), AEA Technology (UK), NIST (USA), PTB (Niemcy).

Zastosowania CRM w laboratoriach:

  • Kalibracja wydajności detektora: CRM o certyfikowanej aktywności pozwala wyznaczić ε(E) dla danej geometrii.
  • Sprawdzanie odzysku chemicznego: wzorzec CRM dodany do matrycy przed procedurą radiochemiczną pozwala ocenić, jaka część analitu zostaje odzyskana.
  • Codzienne sprawdzenie toru: pomiar CRM przed seriami analitycznymi potwierdza stabilność układu.
  • Walidacja metody: porównanie wyniku z wartością certyfikowaną daje bias (błąd systematyczny) metody.
  • Szkolenie personelu: CRM pozwala nowym analitykom ćwiczyć procedurę z materią o znanych właściwościach.

Ekspiracja CRM: po kilku półokresach aktywność wzorca spada do poziomu, gdzie certyfikowana aktywność (po przeliczeniu na dzień pomiaru) ma zbyt dużą relative uncertainty ze względu na niepewność aktywności nominalnej i czasu. Typowy termin użycia to 5–10 T₁/₂ od daty odniesienia.

Wykresy kontrolne (control charts) jako narzędzie stabilności

W laboratorium akredytowanym pomiar CRM nie jest wydarzeniem jednorazowym — jest regularną (tygodniową, miesięczną) procedurą, której wyniki są nanoszone na wykres kontrolny. Wykres Shewhart'a lub Levey-Jennings'a pozwala wizualnie obserwować trendy i odchylenia od wartości oczekiwanej.

Typowe reguły Western Electric (WE):

  • Reguła 1: Jeden punkt poza ±3σ (prawdopodobieństwo przypadkowe < 0,3%) → alarm, nie raportować do czasu zbadania przyczyny.
  • Reguła 2: Dwa kolejne punkty poza ±2σ po tej samej stronie → podejrzenie trendu.
  • Reguła 4: Cztery z pięciu kolejnych punktów poza ±1σ → podejrzenie dryftu.
  • Reguła 8: Osiem kolejnych punktów po tej samej stronie od średniej → podejrzenie systemowego błędu.

Dla detektora HPGe typowe reguły kontrolne:

  • Codzienne sprawdzenie FWHM piku Cs-137 (662 keV): musi być w granicach μ ± 3σ (typowo 1,8–2,0 keV).
  • Codzienne sprawdzenie centroid piku (pozycja w keV): musi stabilnie wskazywać na ±0,1 keV.
  • Tygodniowy pomiar CRM Cs-137 w geometrii Marinelli: wynik musi być w certyfikowanej wartości ±2U (95% CI).

Jeżeli wykres kontrolny sygnalizuje problem, laboratorium musi przeprowadzić dochodzenie przyczyny przed kontynuowaniem produkcji wyników.

Granica wykrywalności (MDA) i granica decyzyjna

Granica wykrywalności (MDA — Minimum Detectable Activity lub MDA w IAEA) i granica decyzyjna (LC — Critical Level, lub L_c) to dwa różne pojęcia definiujące, kiedy sygnał jest odróżnialny od tła. Są kluczowe dla raportowania wyników bliskich zeru.

Granica decyzyjna L_C (ang. Critical Level, Decision Threshold): wartość sygnału, powyżej której stwierdzamy, że aktywność nie jest zerowa. Oparta na błędzie I rodzaju (alfa): jeśli prawdziwa aktywność = 0, to z prawdopodobieństwem 1-α wynik NIE przekroczy L_C. Typowo α = 0,05.

L_C = k_α · u_0

gdzie u_0 jest niepewnością przy zerowej aktywności (wynikającą z tła).

Granica wykrywalności L_D (ang. Detection Limit, Minimum Detectable Activity): minimalna aktywność, którą można wykryć z zadaną mocą testu (beta, błąd II rodzaju). Przy L_D próbka o takiej aktywności będzie zidentyfikowana jako „powyżej zera" z prawdopodobieństwem 1-β. Typowo β = 0,05.

Przybliżony wzór Curie-Currie'go dla L_D przy α = β = 0,05:

L_D = 4,65 · √R_b / t_b + 2,71 / t_s

(dla pomiaru próbki czasu t_s i tła czasu t_b, z tłem R_b zliczeń na sekundę).

Przeliczenie na aktywność: MDA = L_D / (ε × I_γ × m).

Raportowanie „poniżej granicy wykrywalności" oznacza właśnie, że zliczony sygnał netto nie przekroczył L_C. Podanie wartości ujemnej jako wyniku jest błędem — powinna być podana wartość MDA lub L_D z informacją, że wynik jest niższy od granicy detekcji.

Testy biegłości (proficiency testing) i interlaboratoryjne porównania

Testy biegłości to procedura, w której laboratorium otrzymuje próbkę o nieznanej (dla laboratorium) zawartości i raportuje wynik, który jest następnie porównywany z wynikami innych uczestników lub z wartością przypisaną.

W radiometrii testy biegłości organizują m.in.:

  • IAEA ALMERA Network: regularne testy biegłości dla laboratoriów środowiskowych (gleba, woda, biota). Polska (CLOR, IFJ PAN) regularnie uczestniczy.
  • CETAMA (Francja), LGC (UK): europejskie interlaboratoryjne porównania dla różnych matryc.
  • EURAMET: europejskie porównania kluczowe dla wzorców aktywności (kalibratory aktywności, standardy kalibracyjne).

Wyniki testu biegłości są wyrażane jako z-score:

z = (x_lab − x_ref) / σ_PT

gdzie x_lab to wynik laboratorium, x_ref to wartość przypisana (reference value), σ_PT to odchylenie standardowe testu biegłości.

Interpretacja z-score:

  • |z| ≤ 2: wynik zadowalający
  • 2 < |z| ≤ 3: wynik wątpliwy (questionable), wymagana analiza przyczyn
  • |z| > 3: wynik niezadowalający, wymagane działanie korygujące

Laboratoria akredytowane muszą wykazywać uczestnictwo w testach biegłości jako część systemu kontroli jakości wymaganego przez ISO 17025.

Trzy przykłady obliczeniowe

Przykład 1: Pełny budżet niepewności pomiaru aktywności Cs-137 w glebie

Pomiar HPGe, geometria cylindryczna, masa próbki m = 250 g, czas t = 7200 s.

Dane zmierzone:

  • N_s = 2450 zliczeń w fotopiku 662 keV
  • N_b = 180 zliczeń tła w oknie 662 keV (zmierzone przez t_b = 3600 s)
  • ε_peak (z kalibracji wzorcem CRM) = 0,0220 (2,20%)
  • I_γ (intensywność emisji Ba-137m) = 0,8510 (85,10%)
  • m = 250,3 g (waga analityczna)

Obliczenia:
R_s = 2450 / 7200 = 0,3403 s⁻¹
R_b = 180 / 3600 = 0,0500 s⁻¹
R_net = 0,3403 − 0,0500 = 0,2903 s⁻¹

Aktywność: A = R_net / (ε_peak × I_γ) = 0,2903 / (0,0220 × 0,8510) = 0,2903 / 0,01872 = 15,51 Bq

Aktywność właściwa: a = A / m = 15,51 / 0,2503 = 61,97 Bq/kg

Budżet niepewności:

Składowa Wartość Niepewność u_rel [%]
N_s (Poisson) 2450 √2450 = 49,5 → u(R_net) —
N_b (Poisson) 180 √180 = 13,4 → u(R_bg) —
u_rel(R_net) — √(N_s/t_s² + N_b/t_b²)/R_net 3,52%
ε_peak (kalibracja) 0,0220 ±0,0007 (certyfikat) 3,18%
I_γ (jądrowe) 0,8510 ±0,0030 0,35%
Masa m 250,3 g ±0,1 g 0,04%

u_rel(a) = √(3,52² + 3,18² + 0,35² + 0,04²) = √(12,4 + 10,1 + 0,12 + 0) = √22,6 = 4,76%

u(a) = 61,97 × 0,0476 = 2,95 Bq/kg

Wynik: a = (62,0 ± 3,0) Bq/kg Cs-137 (k=1, 68% CI).
Rozszerzona niepewność (k=2): a = (62,0 ± 5,9) Bq/kg (95% CI).

Przykład 2: Obliczenie MDA dla pomiaru Sr-90

Sr-90 jest czystym emiterem beta, więc mierzy się go metodą LSC lub proporcjonalną po radiochemicznym oddzieleniu. Pomiar proporcjonalny:

Czas pomiaru próbki: t_s = 3600 s. Czas tła: t_b = 3600 s. Tło: R_b = 0,020 s⁻¹ (6 zliczeń/5 min, typowe dla komory proporcjonalnej).

Wydajność detekcji Sr-90 (przez Y-90, przy osiągnięciu równowagi): ε = 0,55 (55%).

MDA (wzór Curie-Currie): MDA = (4,65 × √(R_b × t_b) + 2,71) / (t_s × ε × m).

Podstawienie (masa m = 1 g, R_b × t_b = 72 zliczeń):
MDA = (4,65 × √72 + 2,71) / (3600 × 0,55 × 1) = (4,65 × 8,485 + 2,71) / 1980 = (39,45 + 2,71) / 1980 = 42,16 / 1980 = 0,0213 Bq/g = 21,3 mBq/g = 21,3 Bq/kg.

To jest MDA dla 1 g próbki. Dla 10 g: MDA = 2,1 Bq/kg. Dla 100 g: 0,21 Bq/kg — co pozwala osiągnąć czułość wymaganą do analizy środowiskowej tła.

Wniosek: zwiększenie masy próbki (do m = 10–100 g) jest kluczem do obniżenia MDA dla Sr-90, obok optymalizacji tła komory.

Przykład 3: Dryft kalibracji i jego wykrycie z wykresu kontrolnego

Laboratorium codziennie mierzy CRM Cs-137 (certyfikowana aktywność: 1000 ± 30 Bq). Przez 10 kolejnych dni wyniki (Bq):

985, 990, 988, 992, 1005, 1010, 1018, 1025, 1030, 1038

Średnia pierwszych 4 wyników: 988,75 Bq. Ostatnie 4 wyniki to 1018–1038 Bq. Różnica ~50 Bq, czyli ~5% od wartości nominalnej.

Reguła WE nr 4: cztery z pięciu kolejnych punktów poza +1σ (gdzie σ = 30 Bq, więc +1σ = 1030 Bq). Punkty 8–10 (1025, 1030, 1038) oraz punkt 7 (1018) — wszystkie powyżej 1000, a 3 z 4 są powyżej 1030. Reguła alarmuje.

Prawdopodobna przyczyna: dryfowanie wzmocnienia wzmacniacza lub zmiana napięcia zasilania fotopowielacza. Diagnoza: sprawdzenie napięcia HV, ponowna kalibracja energii, pomiar rozdzielczości FWHM. Jeżeli FWHM się zwiększyła, wskazuje na gorsze zbieranie ładunku w detektorze.

Działanie korygujące: przed raportowaniem wyników od dnia 5 należy zbadać, czy wyniki próbek z tego okresu były skorygowane na aktualną kalibrację, czy muszą być powtórzone.

Metody kalibracji wydajności dla różnych zakresów energii

Krzywa wydajności ε(E) detektora promieniowania gamma nie jest liniowa ani prosta — zazwyczaj ma kształt krzywej z maksimum w okolicach 100–300 keV, spadającą ku niskim energiom (z powodu pochłaniania w otulinie i oknie detektora) oraz malejącą ku wysokim energiom (z powodu malejącego przekroju czynnego fotoelektrycznego).

Do kalibracji wydajności stosuje się zestaw punktowych wzorców energetycznych obejmujących zakres roboczych energii. Standardowe izotopy kalibracyjne dla spektrometrii gamma HPGe:

Izotop Energia [keV] T₁/₂ Intensywność γ Uwagi
Am-241 59,54 432,2 lata 35,92% Niska energia — wrażliwy na okno
Cd-109 88,03 461,9 dni 3,61% Niska intensywność — dłuższy czas pomiaru
Co-57 122,06 271,7 dni 85,60% Podstawowy punkt niskiej energii
Cr-51 320,08 27,7 dni 9,91% Środkowy zakres, szybki rozpad
Cs-137 661,66 30,1 roku 85,10% Punkt referencyjny środkowej energii
Mn-54 834,85 312,1 dni 99,98% Wysoka energia, wysoka intensywność
Y-88 898,04 / 1836,06 106,6 dni 93,7% / 99,2% Dwa wysokoenergetyczne punkty
Co-60 1173,24 / 1332,50 5272 dni 99,85% / 99,98% Punkt wysokiej energii, TCS!

Uwaga dotycząca Co-60 i Ba-133: te izotopy mają kaskady gamma z koinbencjami wymagające korekcji sumowania prawdziwego (TCS — True Coincidence Summing). Jeżeli laboratorium używa Co-60 do kalibracji bez korekcji TCS, wyznaczona wydajność będzie zaniżona (bo część zliczeń odpada do kaskad sum).

Metody kalibracji wydajności dla warunków bez gotowego CRM wieloenergetycznego:

  • Mieszanina wzorców punktowych: kilka izotopów w jednym naczyniu (ale bez szans na sumowanie prawdziwe między nimi).
  • Wzorce matrycowe: CRM o certyfikowanym składzie matrycy i aktywności (np. certyfikowana gleba IAEA-375 lub woda IAEA W4).
  • Obliczenia Monte Carlo: gdy wzorzec w identycznej geometrii nie jest dostępny, symulacje GEANT4/MCNP6 pozwalają obliczyć ε(E) dla dowolnej geometrii; wymagają walidacji przez co najmniej jeden punkt doświadczalny.
  • Interpolacja krzywą wielomianową: dopasowanie ln(ε) = a + b·ln(E) + c·[ln(E)]² lub podobnego wyrażenia do punktów eksperymentalnych.

Specyfika kalibracji spektrometrów alfa

Dla spektrometrii alfa — stosowanej do pomiarów Pu, Am, U, Ra, Po — kalibracja wydajności i energii ma odmienną specyfikę niż dla spektrometrii gamma:

  1. Geometria 2π zamiast 4π: próbka jest cienką źródłem alfa napylonym elektrotermicznie lub strąconym po radiochemii na dysku. Detektor powierzchni barierowej lub PIPS umieszczony w odległości kilku mm od źródła widzi geometrycznie ~35–45% całkowitego emisji (z uwagi na kąt bryłowy). Wzorzec kalibracyjny musi być w identycznej geometrii.

  2. Rozdzielczość energetyczna: zależy od grubości i jednorodności źródła. Cienkie źródło (< 5 µg/cm²) daje FWHM ~20 keV; grubsze daje ogon energetyczny ku niższym energiom (samopochłanianie alfa).

  3. Izotopy kalibracyjne alfa:

    • Am-241: 5485,6 keV (85,2%), 5442,8 keV (12,8%), 5388,2 keV (1,4%) — standardowy punkt energii
    • Pu-239: 5156,6 keV (73,3%), 5144,3 keV (15,1%), 5105,6 keV (11,5%) — zbliżona energia do Am-241
    • Po-209 (spike izobaryczny dla Po-210 w próbce)
    • Np-237, Th-229 — stosowane jako trage (spike) w metodzie rozcieńczenia izotopowego alfa-spektrometrii

  4. Metoda wzorca wewnętrznego (spike): do próbki przed mineralizacją dodaje się znane ilości (aktywność znana z certyfikatu) tracera (np. Po-209 dla pomiaru Po-210, Am-243 dla pomiaru Am-241, Pu-242 dla Pu-239/240). Po radiochemii oba izotopy są w tej samej próbce, więc ich wydajność chemiczna i geometria są identyczne. Stosunek aktywności próbka/spike eliminuje wydajność jako czynnik kalibracyjny.

  5. Interferencje izobaryczne: Am-241 (5486 keV) może interferować z Cm-244 (5805 keV) lub Pu-238 (5499 keV). Dobra rozdzielczość energetyczna (<25 keV FWHM) pozwala oddzielić te piki.

Rola wzorców jądrowych w skali krajowej i międzynarodowej

Hierarchia spójności pomiarowej dla aktywności promieniowania:

BIPM (Biuro Miar i Wag): na szczycie piramidy. Utrzymuje wzorzec Système International (SI) poprzez komitet CCRI (Comité Consultatif pour les Rayonnements Ionisants). Porównania kluczowe między krajowymi laboratoriami wzorcowymi prowadzone są przez BIPM KCDB (Key Comparison Database).

Krajowe laboratoria wzorcowe (NMI): w Polsce rolę NMI dla aktywności radiologicznej pełni POLATOM (Narodowe Centrum Badań Jądrowych, Ośrodek Radioizotopów w Świerku). Posiada akredytację PCA jako laboratorium wzorcujące dla aktywności izotopów. PTB (Niemcy) i NPL (UK) są głównymi NMI europejskimi. NIST pełni tę rolę w USA.

Akredytowane laboratoria kalibracyjne: otrzymują wzorce od NMI z certyfikowanymi aktywnościami, następnie kalibrują sprzęt laboratoriów niższego rzędu lub certyfikują własne CRM.

Laboratoria badawcze (akredytowane): otrzymują wzorce od laboratoriów kalibracyjnych lub bezpośrednio od NMI. Używają wzorców do kalibracji detektorów i metod. Raportują wyniki z budżetem niepewności wynikającym z GUM.

Spójność dla żywotności: wzorce aktywności ulegają rozpadowi. Certyfikat podaje aktywność na datę odniesienia i niepewność tej aktywności. Użytkownik oblicza aktywność na dzień pomiaru: A(t) = A_0 × e^{−λt}. Niepewność aktywności aktualnej zawiera niepewność aktywności wzorca A_0 i niepewność czasu t (zazwyczaj pomijalna) oraz wartość T₁/₂ (z tabel DDEP lub IAEA) z jej własną niepewnością.

Praktyczne wyzwania w laboratoriach polskich

Laboratoria środowiskowe w Polsce (WIOŚ — Wojewódzkie Inspektoraty Ochrony Środowiska, CLOR, SANEPID) wykonują rutynowe pomiary gammaspectrometryczne i radiochemiczne. Ich specyfika obejmuje:

  • Serie próbek i wydajność: duże laboratoria analizują setki próbek miesięcznie. Czas pomiarów (od godziny do 24 godzin na próbkę dla HPGe) jest ograniczeniem przepustowości. Optymalizacja czasu pomiaru (t_s vs MDA vs koszty) jest realnym problemem praktycznym.
  • Certyfikacja metod przez GIOŚ/EEA: laboratoria biorące udział w ogólnopolskim monitoringu (Inspekcja Ochrony Środowiska) muszą stosować metody certyfikowane przez GIOŚ lub normalizowane (np. ISO 10703 dla aktywności gamma w wodzie, ISO 9696 dla Ra-226).
  • Sprawdzanie czujności przez POLATOM: POLATOM dostarcza wzorce kalibracyjne i certyfikowane próbki dla polskich laboratoriów; jest punktem kontaktowym dla kalibracji wzorców krajowych i dostawcą standardów aktywności.
  • EURDEP i raportowanie do Europy: przez platformę EURDEP (European Radiological Data Exchange Platform) polskie dane z sieci stacji monitoringowych (DOZ, sieci wczesnego ostrzegania) trafiają do europejskiej bazy danych. Wymaga to jednolite sformalizowanego formatu wyników z określoną niepewnością.
  • Moc obliczeniowa a czas pomiarów: nowoczesne laboratoria stosują oprogramowanie analizy widm (GammaVision, Genie 2000, Ortec Pro, Canberra) z automatycznym wyznaczaniem niepewności. Jednak użytkownik musi rozumieć, co oprogramowanie robi, by poprawnie interpretować raporty i wykrywać anomalie.

Pytania otwarte

  1. W jaki sposób korelacja między błędem tła a błędem piku (gdy oba są mierzone tym samym detektorem) wpływa na prawidłowe zastosowanie prawa propagacji niepewności GUM? Czy standardowy wzór kwadratowy jest właściwy dla tej sytuacji?

  2. Jakie minimum danych jest potrzebne do walidacji metody radiometrycznej zgodnie z ISO 17025? Czym różni się walidacja metody od sprawdzenia wewnętrznej zgodności wyników?

  3. Jak wyznaczenie MDA zmienia się, gdy stosuje się zamiast reguły Currie'go podejście bayesowskie? Kiedy ta różnica ma praktyczne znaczenie dla decyzji regulacyjnych?

  4. Jak skalibrować wydajność detekcji dla energii, dla których nie ma dostępnych certyfikowanych wzorców punktowych (np. zakres 3–30 keV)? Jakie metody interpolacji i jakie źródła informacji o przekrojach czynnych są stosowane?

  5. Jak minimalizować błąd systematyczny związany z samopochłanianiem w próbkach o niejednorodnym składzie (np. gleba z frakcją organiczną i mineralną)? Kiedy lepiej stosować korekcję obliczeniową, a kiedy przygotowanie wzorców matrycowych?

  6. Dlaczego samo uczestnictwo w testach biegłości bez oceny wyników i działań korygujących nie spełnia wymagań ISO 17025? Jakie procedury zarządzania niezgodnościami muszą być wdrożone?

  7. Jak efekt dryftu napięcia wysokiego (HV) fotopowielacza różni się od efektu dryftu wzmocnienia wzmacniacza w spektrometrze scyntylacyjnym? Jak każdy z nich manifestuje się na wykresie kontrolnym?

  8. Jak budżet niepewności dla aktywności próbki zmienia się, gdy zamiast metody bezwzględnej stosuje się metodę względną (porównanie z wzorcem w tej samej geometrii)? Jakie składowe się eliminują, a jakie pozostają?

Podsumowanie dydaktyczne

  1. Kalibracja = definicja wyniku — bez kalibracji energii i wydajności liczba zliczeń jest tylko odczytem aparatury, nie mierzalną wielkością fizyczną. Kalibracja dostarcza matematycznego modelu (funkcji kalibracyjnej), który zamienia sygnał w wynik fizyczny wyrażony w Bq, Bq/kg lub innej jednostce.

  2. Niepewność to struktura, nie pesymizm — budżet niepewności nie jest dopiskiem na końcu. Jest mapą tego, co w pomiarze wiemy dobrze (masa próbki, czas pomiaru), a co gorzej (kalibracja wydajności, odzysk chemiczny). Mówi, które składowe naprawić, aby poprawić wynik.

  3. Typy A i B są komplementarne — niepewność typu A (ze statystyki powtórzeń) i B (z certyfikatów, tolerancji, doświadczenia) opisują różne źródła nieoznaczoności. Samo powtórzenie pomiaru redukuje typ A, ale nie wykryje złego wzorca kalibracyjnego (błąd B).

  4. Tło jest pomiarem, nie zerowaniem — tło ma własne fluktuacje Poissonowskie i własną niepewność. Złe wyznaczenie tła (zbyt krótkie, w złej geometrii, bez powtórzeń) może dominować w budżecie niepewności pomiaru o niskiej aktywności. Czas pomiaru tła jest parametrem równie ważnym jak czas pomiaru próbki.

  5. MDA i L_C są decyzyjnymi narzędziami — granica wykrywalności (MDA) i granica decyzyjna (L_C) nie są właściwościami detektora, lecz funkcją tła, czasu, wydajności i masy próbki. Poprawne raportowanie „poniżej MDA" wymaga podania wartości MDA — samo „nie wykryto" jest informacyjnie niekompletne.

  6. Wykresy kontrolne chronią serię wyników — rutynowe pomiary CRM i wykresy kontrolne Shewhart'a/Levey-Jennings'a są systemem wczesnego wykrywania dryftu, niestabilności i awarii sprzętu. Bez nich systematyczny błąd może pozostać niewykryty przez tygodnie.

  7. Spójność pomiarowa (traceability) jest warunkiem porównywalności — wynik laboratorium A jest porównywalny z wynikiem laboratorium B tylko wówczas, gdy obydwa mają udokumentowany łańcuch spójności do tego samego wzorca wyższego rzędu (NIST, PTB, BIPM). Uczestnictwo w testach biegłości IAEA ALMERA jest praktycznym sprawdzeniem tej spójności.

  8. GUM, ISO 17025 i akredytacja = wspólny język — standardy GUM i ISO 17025 są wspólnym językiem laboratoriów radiometrycznych na całym świecie. Ich zastosowanie przez CLOR, IFJ PAN i laboratoryjne sieci WIOŚ pozwala Polsce raportować dane porównywalne z danymi z Finlandii, Niemiec czy Japonii. To jest dlaczego metrologiczna kultura w laboratorium ma znaczenie poza samym laboratorium.

Dodatkowe materiały multimedialne

Warto przygotować kalkulator budżetu niepewności dla prostego pomiaru gamma: użytkownik wpisuje zliczenia próbki, tła, czas, masę, wydajność i prawdopodobieństwo emisji. Narzędzie pokazuje aktywność, wkłady niepewności i ostrzeżenie, który składnik dominuje.

Najkrótsze podsumowanie: kalibracja mówi, jak odczyt aparatury zamienić na wielkość fizyczną, a niepewność mówi, jak bardzo można tej zamianie ufać. Bez obu elementów wynik jest tylko liczbą z urządzenia.

Ćwiczenia praktyczne

Pierwsze ćwiczenie z tłem: student dostaje N_s, t_s, N_b, t_b dla pomiaru próbki i tła. Ma obliczyć częstość netto oraz niepewność u(R_net). Następnie porównuje dwa warianty: długie tło i krótkie tło, aby zobaczyć, kiedy tło dominuje budżet.

Drugie ćwiczenie kalibracji energii: na podstawie trzech znanych pików dopasować prostą E = a * k + b, obliczyć energię piku nieznanego i sprawdzić, jak wynik zmienia się po usunięciu jednego punktu kalibracyjnego. Celem jest pokazanie, że kalibracja ma własną stabilność i niepewność.

Trzecie ćwiczenie wydajności: użyć uproszczonego wzoru A = N_net / (t * epsilon * I_gamma * m) dla danych syntetycznych. Następnie zwiększyć niepewność epsilon z 3% do 15% i sprawdzić, czy dalsze wydłużanie czasu pomiaru nadal istotnie poprawia wynik.

Czwarte ćwiczenie geometrii: porównać dwa raporty z pomiaru tej samej próbki. W jednym wzorzec i próbka mają tę samą geometrię, w drugim różną. Student ma wskazać, który raport wymaga dodatkowej poprawki wydajności i dlaczego.

Piąte ćwiczenie audytowe: przygotować „kartę wyniku” dla wymyślonego pomiaru K-40 w materiale budowlanym. Trzeba wypisać detektor, czas, tło, masę, geometrię, kalibrację energii, wydajność, wynik i niepewność. Jeśli któregoś elementu brakuje, wynik należy oznaczyć jako orientacyjny.

Przejdź do ćwiczenia interaktywnego

Powiązane artykuły