Kontrola jakości w radioanalizie: masa, pipetowanie, odzysk chemiczny i niepewność

Streszczenie

Radioanaliza kojarzy się z detektorami, widmami i nazwami radionuklidów, ale wynik często psuje się dużo wcześniej: przy ważeniu, odmierzaniu objętości, utracie części próbki, źle opisanej próbie ślepej albo pominiętym odzysku chemicznym. Spektrometr może działać idealnie, a wynik nadal będzie zły, jeśli do obliczeń trafi błędna masa lub nieznana część analitu zginie w przygotowaniu próbki.¹

Ten artykuł pokazuje radioanalizę od strony kontroli jakości. Najważniejsze pojęcia to tara i masa netto, typ pipety TC/TD, klasa szkła miarowego, powtarzalność, próba ślepa, odzysk chemiczny, wzorzec, próbka kontrolna, tło, wydajność pomiarowa i budżet niepewności.

Rozszerzenie tematu

Wynik zaczyna się przed detektorem

W pomiarze radiometrycznym łatwo skupić uwagę na aparaturze: napięciu detektora, widmie, kanale analizatora, piku gamma albo liczbie impulsów. To są ważne elementy, ale nie zastępują zwykłej metrologii chemicznej i masowej. Jeżeli próbka została źle zważona, roztwór źle odmierzony, a tło opisane nieporównywalnie, to nawet poprawne zliczenia dadzą niepoprawne Bq/kg.

Przewodnik KChRS zaczyna blok ćwiczeń radioanalitycznych od bardzo przyziemnego zadania: pipetowania roztworów i pomiaru masy na wadze analitycznej.¹ To nie jest dodatek dla początkujących. To fundament późniejszych obliczeń aktywności, wydajności i niepewności.

W radioanalizie łańcuch pomiarowy wygląda zwykle tak:

• pobrana próbka ma masę albo objętość,
• część próbki trafia do przygotowania,
• część analitu może zostać utracona albo rozcieńczona,
• próbka po przygotowaniu jest mierzona w określonej geometrii,
• detektor daje zliczenia,
• zliczenia są przeliczane na aktywność,
• aktywność jest dzielona przez masę, objętość albo inną wielkość odniesienia.

Błąd na początku tego łańcucha przechodzi do końcowego wyniku. Dlatego kontrola jakości nie jest formalnością po pomiarze. Jest warunkiem, żeby pomiar miał sens.

Masa: tara, netto i waga analityczna

Najprostszy przykład to ważenie. Wynik aktywności właściwej:

a_m = A / m

jest tak dobry, jak dobra jest masa m. Jeżeli masa próbki jest zawyżona, aktywność właściwa zostanie zaniżona. Jeżeli masa jest zaniżona, wynik w Bq/kg wzrośnie.

W ćwiczeniu KChRS podkreślono wagę analityczną o typowej dokładności rzędu 0,1 mg oraz konieczność rozdzielenia masy naczynka od masy dostarczonej cieczy.¹ W języku metrologii chodzi o tarę i masę netto:

m_netto = m_brutto - m_tara.

Ten wzór jest prosty, ale źródła błędu są praktyczne: naczynko może nie być suche, może mieć ślady z palców, masa może zmieniać się przez parowanie, a zapis może mieszać miligramy z gramami. W radioanalizie takie drobiazgi mogą mieć ten sam rząd znaczenia co statystyka zliczeń.

Objętość: pipeta nie jest tylko rurką

Pipety i szkło miarowe mają własne założenia kalibracyjne. KChRS rozróżnia pipety TC i TD: TC oznacza naczynie skalibrowane "to contain", a TD "to deliver". Dla pipety TD resztka cieczy pozostająca w końcówce jest częścią założenia kalibracji i nie powinna być traktowana tak samo jak przy pipetach TC.¹

To ma bezpośrednie znaczenie dla radioanalizy. Jeżeli student albo operator zmienia sposób użycia pipety, zmienia rzeczywistą objętość dostarczoną do próbki. Jeżeli objętość jest używana do obliczenia stężenia, rozcieńczenia albo odzysku, błąd przechodzi dalej.

Szkło miarowe ma też klasy dokładności. W przewodniku przypomniano, że klasa A i B różnią się dopuszczalnym błędem objętości, a wymagania klasy B są mniej rygorystyczne.¹ To nie jest kwestia estetyki sprzętu. To decyzja, czy niepewność objętości będzie pomijalna, czy stanie się istotnym składnikiem budżetu.

Powtarzalność pipetowania

Jedno poprawne pipetowanie nie mówi jeszcze, jak powtarzalna jest metoda. Dlatego w ćwiczeniu KChRS porównuje się kilka typów pipet i wykonuje powtórzenia, a potem oblicza średnią masę dostarczonej wody oraz odchylenie standardowe średniej.¹

W praktyce ten eksperyment uczy trzech rzeczy:

• błąd systematyczny: średnia objętość może być przesunięta względem wartości nominalnej,
• błąd przypadkowy: kolejne dostarczenia tej samej objętości mają rozrzut,
• zależność od operatora: technika pracy może być równie ważna jak klasa pipety.

Dla późniejszej radioanalizy ważne jest rozdzielenie tych składników. Powtarzalność mówi, jak bardzo rozchodzą się wyniki serii. Poprawność mówi, czy średnia trafia w wartość odniesienia. Precyzyjne, ale przesunięte pipetowanie nadal daje błędny wynik.

Próba ślepa

Próba ślepa jest odpowiedzią na pytanie: co zobaczyłaby metoda, gdyby nie było analizowanego radionuklidu w próbce? W spektrometrii gamma może to być pojemnik z wodą dejonizowaną albo materiałem matrycowym bez oczekiwanego radionuklidu. W radiochemii może to być pełny przebieg przygotowania bez próbki właściwej.

KChRS używa prób ślepych w ćwiczeniach spektrometrii żywności i materiałów budowlanych, gdzie trzeba wyznaczyć zliczenia pod fotopikami w konfiguracji porównywalnej z próbką.² Sens jest ogólny: odejmujemy nie tylko "tło sali", lecz tło całej konfiguracji pomiarowej.

Próba ślepa wykrywa kilka problemów:

• zanieczyszczenie odczynników lub naczyń,
• radionuklidy obecne w materiale pojemnika,
• tło aparaturowe w danej geometrii,
• błędne przypisanie piku,
• przeniesienie śladów między próbkami.

Jeżeli próba ślepa nie jest mierzona albo nie jest porównywalna, wynik netto może wyglądać dokładnie, ale być przesunięty.

Wzorzec i próbka kontrolna

Wzorzec jest materiałem o znanej wartości odniesienia. Może służyć do kalibracji wydajności, sprawdzenia poprawności metody albo kontroli stabilności aparatury. Próbka kontrolna jest mierzona nie po to, żeby odkryć nową wartość, lecz żeby sprawdzić, czy metoda nadal daje wynik zgodny z oczekiwaniem.

W ćwiczeniu z licznikiem Geigera-Müllera KChRS pokazuje zasadę wyznaczania wydajności pomiarowej na próbce o znanym stężeniu aktywnościowym i składzie identycznym jak próbki nieznane.³ To zdanie jest bardzo ważne: wydajność pomiarowa nie jest uniwersalną stałą licznika. Zależy od rodzaju próbki i geometrii.

Ten sam sens występuje w spektrometrii gamma i alfa. Detektor, pojemnik, odległość, energia promieniowania i matryca tworzą jeden układ pomiarowy. Wzorzec powinien sprawdzać właśnie ten układ, a nie abstrakcyjną zdolność detektora do rejestrowania promieniowania.

Odzysk chemiczny

Odzysk chemiczny mówi, jaka część analitu przeszła przez przygotowanie próbki i znalazła się w preparacie mierzonym. Jeśli w procesie przygotowania tracimy 20% radionuklidu, a wynik policzymy tak, jakby nic nie zginęło, aktywność w próbce wyjściowej będzie zaniżona.

Ogólnie:

R = ilość_odzyskana / ilość_początkowa.

Jeżeli mierzony preparat odpowiada tylko części analitu, wynik trzeba skorygować:

A_początkowa = A_zmierzona / R.

W radioanalizie środowiskowej odzysk bywa wyznaczany przez znacznik, wzorzec dodany do próbki albo próbkę wzbogaconą. W tym artykule ważna jest sama zasada, nie procedura chemiczna: odzysk jest składnikiem równania pomiarowego i ma własną niepewność. Nie wolno traktować go jako "szczegółu przygotowania".

Odzysk a wydajność detektora

Odzysk chemiczny i wydajność detektora to różne rzeczy. Odzysk mówi, jaka część radionuklidu dotarła do preparatu. Wydajność detektora mówi, jaka część emisji z preparatu została zarejestrowana jako użyteczny sygnał.

Schematycznie:

R_net = A_próbki * R_chem * epsilon * p * inne_poprawki.

Stąd:

A_próbki = R_net / (R_chem * epsilon * p * inne_poprawki).

Jeśli R_chem = 0,75, a pominiemy ten czynnik, wynik będzie o około jedną czwartą za mały. Jeśli pomylimy odzysk chemiczny z wydajnością detekcji, możemy skorygować wynik dwa razy albo nie skorygować go wcale.

Błąd przypadkowy i systematyczny

Kontrola jakości wymaga rozróżnienia dwóch rodzin błędów. Błąd przypadkowy powoduje rozrzut wyników wokół średniej. Błąd systematyczny przesuwa całą serię w jedną stronę.

Przykłady błędów przypadkowych:

• statystyczne fluktuacje zliczeń,
• drobne różnice kolejnych pipetowań,
• szum tła,
• niewielki rozrzut mas przy powtarzanym ważeniu.

Przykłady błędów systematycznych:

• błędnie skalibrowana pipeta,
• zła masa tary,
• wydajność detektora z innej geometrii,
• pominięty odzysk chemiczny,
• użycie wzorca o innej matrycy niż próbka,
• stałe zanieczyszczenie próby ślepej.

Dziunikowski i Kalita przypominają, że powtarzane pomiary nigdy nie są idealnie powtarzalne, a wartość prawdziwa pozostaje nieznana i musi być szacowana.⁴ To dobry język dla radioanalizy: wynik jest estymatą wraz z niepewnością, nie objawioną wartością.

Średnia, odchylenie i mała liczba powtórzeń

Jeżeli wykonano n powtórzeń, podstawowym oszacowaniem wyniku jest średnia:

x_srednie = (x1 + x2 + ... + xn) / n.

Rozrzut opisuje odchylenie standardowe. Dla małej liczby powtórzeń nie należy udawać, że średnia jest znana tak dobrze jak przy dużej próbie. Dziunikowski i Kalita omawiają użycie rozkładu Studenta do szacowania przedziałów ufności przy kilku lub kilkunastu danych pomiarowych.⁴

Dla artykułów serwisu praktyczna lekcja jest prosta: jeśli wynik opiera się na trzech powtórzeniach, należy mówić o ograniczonej precyzji. Sama średnia bez rozrzutu jest informacją niepełną.

Propagacja niepewności w radioanalizie

Typowy wynik może mieć postać:

A = (N/t - B/t_B) / (epsilon * p * R_chem).

Jeżeli chcemy aktywność właściwą, dzielimy jeszcze przez masę:

a_m = A / m.

Każdy czynnik ma niepewność. Zliczenia mają składnik Poissona. Czas jest zwykle bardzo dobrze znany, ale nie zawsze pomijalny przy krótkich pomiarach automatycznych. Wydajność pochodzi z kalibracji. Intensywność linii pochodzi z danych jądrowych. Odzysk chemiczny pochodzi z kontroli procesu. Masa pochodzi z ważenia.

Dla iloczynów i ilorazów względne niepewności często dodaje się kwadratowo:

u_rel(a_m)^2 ≈ u_rel(R_net)^2 + u_rel(epsilon)^2 + u_rel(p)^2 + u_rel(R_chem)^2 + u_rel(m)^2.

To jest przybliżony zapis dydaktyczny. W prawdziwym budżecie trzeba sprawdzić korelacje, sposób wyznaczenia tła i model dopasowania piku. Ale nawet ten prosty wzór pokazuje, że wynik nie jest "statystyką zliczeń plus trochę marginesu".

Granica wykrywalności i fałszywa pewność

Kontrola jakości obejmuje też pytanie, czy metoda mogła wykryć dany radionuklid. Jeśli wynik jest poniżej granicy wykrywalności, nie należy pisać "zero". Poprawniej jest podać, że aktywność była poniżej określonej granicy detekcji w danych warunkach pomiaru.

Granica wykrywalności zależy od tła, czasu pomiaru, wydajności, masy próbki, linii gamma lub alfa oraz kryterium statystycznego. Dwie próbki o tym samym rzeczywistym stężeniu aktywności mogą dać różne komunikaty, jeśli jedna była mierzona krótko w dużym tle, a druga długo w dobrej osłonie.

Dla edukacyjnego serwisu ważne jest, aby nie mylić trzech zdań:

• "nie wykryto" oznacza brak statystycznie istotnego sygnału,
• "poniżej MDA" oznacza wynik mniejszy niż możliwości metody w tej konfiguracji,
• "równe zero" jest prawie nigdy nieuzasadnione w pomiarach śladowych.

Karta kontroli jakości

Dobry raport z radioanalizy powinien mieć metadane. Minimalna karta kontroli jakości może zawierać:

• identyfikator próbki,
• masę albo objętość,
• opis matrycy,
• typ naczynia i geometrię,
• datę przygotowania i pomiaru,
• czas pomiaru,
• wynik próby ślepej,
• wynik próbki kontrolnej,
• wydajność pomiarową,
• odzysk chemiczny, jeśli dotyczy,
• sposób odjęcia tła,
• niepewność i jej składniki,
• osobę lub system wykonujący analizę.

Bez takiej karty trudno później rozstrzygnąć, czy nietypowy wynik jest realny, czy wynika z błędu przygotowania, aparatury albo opisu.

Prosty przykład: wpływ odzysku

Załóżmy syntetyczny przykład: z pomiaru preparatu otrzymano aktywność A_zmierzona = 120 Bq. Przygotowanie próbki miało odzysk chemiczny R = 0,80. Aktywność w próbce początkowej wynosi:

A_początkowa = 120 / 0,80 = 150 Bq.

Jeżeli masa próbki początkowej wynosiła 0,50 kg, aktywność właściwa:

a_m = 150 / 0,50 = 300 Bq/kg.

Gdyby pominąć odzysk, wynik wyniósłby 240 Bq/kg. Różnica nie wynika z detektora. Wynika z chemii i masy.

Prosty przykład: wpływ masy

Drugi przykład pokazuje, dlaczego należy pilnować jednostek. Jeśli aktywność próbki wynosi 75 Bq, a masa to 25 g, nie wolno dzielić przez 25, jeśli wynik ma być w Bq/kg. Najpierw:

25 g = 0,025 kg.

Potem:

a_m = 75 / 0,025 = 3000 Bq/kg.

Ten błąd jest banalny, ale w praktyce bardzo groźny: pomylenie gramów i kilogramów daje błąd tysiąckrotny. Żaden zaawansowany model nie uratuje wyniku, jeśli jednostki są wpisane źle.

Kontrola jakości jako ochrona przed nadinterpretacją

W serwisie poświęconym atomistyce i historii technologii jądrowej kontrola jakości ma jeszcze jedną rolę: chroni przed atrakcyjnymi, ale fałszywie pewnymi wnioskami. Jeśli materiał "ma pik", nie znaczy jeszcze, że jego aktywność została dobrze policzona. Jeśli wynik różni się od innego laboratorium, nie znaczy od razu, że jedna strona się myli co do fizyki. Różnice mogą wynikać z masy, geometrii, tła, wydajności, odzysku i matrycy.

Dlatego artykuły o spektrometrii, radiochemii środowiskowej, aktywacji neutronowej czy analizie materiałów powinny odsyłać do kontroli jakości. To wspólny język dla wszystkich tych tematów.

Procedury spike-and-recovery: metodyka i obliczenia

Metoda spike-and-recovery jest standardową procedurą walidacyjną w radioanalizie. Polega na dodaniu do matrycy próbki (np. do gleby, wody, żywności) wzorca radiochemicznego o ściśle znanych właściwościach (znany radionuklid, znana aktywność, certyfikowany) przed wykonaniem pełnej procedury przygotowania próbki. Po przeprowadzeniu całego procesu analitycznego mierzy się aktywność dodanego radionuklidu w preparacie i porównuje z aktywnością spodziewaną.

Wzór odzysku (recovery R):
R = A_znalezione / A_dodane × 100%

Gdzie A_znalezione to aktywność radionuklidu wzorca wykryta w preparacie końcowym (skorygowana o rozpad, geometrię i wydajność detekcji), a A_dodane to aktywność wzorca dodanego do próbki.

Spike-and-recovery spełnia kilka funkcji jednocześnie:

• Sprawdzanie strat na etapach radiochemicznych: wytrącanie, filtracja, przemywanie, odparowanie — każdy etap może powodować straty analitu. Spike umożliwia ich skwantyfikowanie.
• Walidacja czasu pracy operatora: indywidualna sprawność technologa radiochemicznego jest mierzona przez odzysk jego spiked samples.
• Weryfikacja odczynników i ich partii: zanieczyszczenia nowej partii odczynnika mogą blokować współstrącanie lub niszczyć kolumnę jonowymienną.
• Weryfikacja specyfiki metody dla matrycy: metoda opracowana dla wody morskiej może mieć inny odzysk w glebie bogatej w materię organiczną.

Akceptowalne zakresy R według różnych protokołów:

Metoda/Regulacja	Minimalny R	Maksymalny R
ISO 17025 (ogólnie)	70%	130%
IAEA ALMERA	75%	115%
US EPA Methods	60%	125%
Typowe wytyczne CLOR	80%	110%

Spike-and-recovery a korekta wyniku: jeżeli odzysk jest znany i powtarzalny (np. R = 0,82 ± 0,03), wynik może być skorygowany przez podzielenie przez R. Jeżeli odzysk jest zmienny (np. zakres 0,55–0,95 w tej samej serii), korekta wprowadza dużą niepewność i metoda wymaga poprawy — nie wystarczy samo przeliczenie.

Rodzaje blank samples i ich zastosowanie

Próba ślepa to materiał przygotowany równolegle z próbkami właściwymi, przez który przepuszcza się tę samą procedurę analityczną, ale który nie zawiera analitu (lub zawiera go w ściśle kontrolowanej, zerowej ilości). Termin "blank" kryje kilka różnych typów, które spełniają różne cele:

Method blank (ślepa metodyczna): pełna procedura analityczna (ważenie, mineralizacja, rozdzielanie, odparowanie, pomiar) wykonana na materiale matrycowym bez analitu (np. czysta piaskownica woda dejonizowana). Wykrywa zanieczyszczenia z odczynników, naczyń i środowiska.

Reagent blank (ślepa odczynnikowa): same odczynniki bez matrycy. Izoluje wkład odczynników od wkładu matrycy. Używana np. gdy podejrzewa się zanieczyszczenie konkretnej partii HNO₃ lub H₂SO₄.

Field blank (ślepa terenowa): pojemnik z czystą wodą dejonizowaną lub matrycą referencyjną, który jedzie z ekipą próbkującą, jest otwierany w miejscu poboru, szczelnie zamykany i przywożony do laboratorium. Wykrywa zanieczyszczenia z transportu, powietrza lub otoczenia poboru.

Processing blank (ślepa procesowa): wykonywana po serii silnie aktywnych próbek, by sprawdzić, czy nie doszło do przeniesienia krzyżowego (carryover) między próbkami.

Instrument blank (ślepa aparaturowa): czyste rozpuszczalnik lub woda dejonizowana mierzone bezpośrednio detektorem, bez przygotowania próbki. Wyznacza tło aparatury i geometrii pomiarowej.

W rzetelnym laboratorium wszystkie te typy blanków są stosowane naprzemiennie i planowo, a ich wyniki trafiają do systemu jakości. Pominięcie jakiegokolwiek z nich pozostawia lukę w kontroli nad procesem analitycznym.

Metody radiochemicznego rozdzielania i wyzwania dla odzysku

Właściwy pomiar aktywności wielu radionuklidów środowiskowych (Sr-90, Pu-239/240, Am-241, Ra-226, Po-210) wymaga chemicznego oddzielenia interesującego nas nuklidu od matrycy i interferujących izotopów przed pomiarem. Każda technika rozdzielania ma swój charakterystyczny odzysk i zmienność:

Współstrącanie (coprecipitation): Sr-90 strąca się jako szczawian strontu lub węglan baru-strontu. Pu(IV) strąca się przez redukcję do Pu(III) z La(OH)₃ lub TiO(OH)₂. Odzysk: 80–95%, wrażliwy na pH, temperaturę i czas trawienia.

Ekstrakcja ciekło-ciekło (solvent extraction): Am i Cm wyciągane z roztworów amoniakalnych do HDEHP (bis(2-ethylhexyl)phosphoric acid) w C₄ frakcjach hekzanu. Odzysk: 85–98%, wrażliwy na stosunki objętości faz i czas mieszania.

Chromatografia jonowymienien (ion exchange): mocna anionowymiennicza żywica (np. Dowex 1×8 w HNO₃) rozdziela Fe/Zr/Np/Pu/Am przez różne powinowactwa kompleksowych form azotanowych. Odzysk: 75–95%, zależny od stężenia HNO₃ i szybkości przepływu.

Ekstrakcja stała (solid phase extraction, SPE): kolumny UTEVA (U), TRU (transplutonics), DGA (ogólna), TEVA (Tc, U) od Eichrom Technologies. Pozwalają na szybsze i powtarzalne rozdzielanie w porównaniu z klasyczną chromatografią. Odzysk: 85–99% dla optymalnych warunków.

Destylacja: Po-210 (lotny w zakresie 800°C), Cs (lotny w warunkach HCl+HNO₃), radon (lotny w formie gazowej). Odzysk: 70–95%, wrażliwy na temp. pieca i szczelność układu.

W każdej z tych technik odzysk chemiczny jest zmienny między seriami i operatorami. Właśnie dlatego walidacja metodą spike-and-recovery i regularne sprawdzenia CRM są nieodzowne dla wiarygodności wyników.

Traceability i referencyjne materiały matrycowe

Traceability (spójność pomiarowa) w radioanalizie oznacza, że każdy wynik ilościowy jest powiązany z krajowym lub międzynarodowym wzorcem przez udokumentowany łańcuch pomiarów o znanych niepewnościach na każdym stopniu.

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

• Wzorce pierwotne (POLATOM, NIST, PTB): aktywność certyfikowana metodami absolutnymi (4π-β-γ koincydencja, elektrometria, kalorimetria). Niepewność standardowa: 0,5–3%.
• Materiały referencyjne IAEA: IAEA RGU-1 (ruda uranu), IAEA SOIL-7 (gleba), IAEA W-4 (woda morska), IAEA 300 (ruchomy Cs-137 w suchej glebie). Certyfikowane wartości z niepewnościami z wielolaboratoryjnych porównań.
• Krajowe CRM (POLATOM, LGC, Eckert & Ziegler): punktowe lub matrycowe wzorce aktywności. Certyfikowane w odniesieniu do wzorców wyższego rzędu.
• In-house (wzorce wewnętrzne laboratoriów): wzorce przygotowane przez laboratorium z CRM wyższego rzędu, przechowywane i cyklicznie kalibrowane.

Materiały matrycowe CRM są kluczowe, gdy metoda wymaga walidacji dla specyficznej matrycy. CRM gleby (np. IAEA SOIL-7: Cs-137 = 103 ± 5 Bq/kg, Ra-226 = 16,0 ± 1,5 Bq/kg) pozwala sprawdzić, czy metoda ektraktuje radionuklidy równie efektywnie z prawdziwej gleby jak z wodnego wzorca punktowego.

Wyzwania z traceability w praktyce:

• Wiele laboratoriów kalibruje HPGe wzorcem punktowym Cs-137, ale mierzy gleby w Marinelli. Geometria jest inna, co wymaga osobnej kalibracji lub symulacji Monte Carlo.
• CRM matrycowe mają swoją datę ważności — po kilku półokresach składniki krótkotrwałe zanikają i CRM staje się niekompletnym sprawdzeniem.
• Import certyfikowanych materiałów (NIST, Eckert & Ziegler) wiąże się z kosztami i formalnościami ADR (transport substancji promieniotwórczych). Część małych laboratoriów polskich pomija matrycowe CRM, stosując tylko wzorce punktowe — zwiększa to ryzyko błędu matrycowego.

Trzy przykłady obliczeniowe

Przykład 1: Obliczenie odzysku chemicznego Sr-90 metodą spike-and-recovery

Do próbki gleby o masie 50 g dodano spikę: Am-243 o aktywności A_spike = 2,00 Bq. Po przeprowadzeniu pełnej procedury radiochemicznej (mikrofalowa mineralizacja, separacja TRU SPE) zmierzono aktywność Am-243 w preparacie: A_found = 1,64 Bq (skorygowane o czas od dodania spiki do pomiaru i ε_detekcji).

R = A_found / A_spike = 1,64 / 2,00 = 0,82 (82%).

Jednocześnie zmierzono aktywność Am-241 w preparacie: A_Am241_measured = 0,048 Bq.

Skorygowana aktywność Am-241 w próbce: A_Am241 = A_Am241_measured / R = 0,048 / 0,82 = 0,0585 Bq.

Aktywność właściwa: a = A_Am241 / m = 0,0585 / 0,050 = 1,17 Bq/kg.

Bez korekty odzysku: 0,048 / 0,050 = 0,96 Bq/kg — błąd ~18%.

Przykład 2: Budżet niepewności z uwzględnieniem odzysku chemicznego

Pomiar Sr-90 w glebie:

• Aktywność zmierzona (bez korekty): A_meas = 15,0 Bq, u_rel(A_meas) = 4,2% (z propagacji Poisson tła i sygnału)
• Odzysk: R = 0,85, u(R) = 0,04 (wyznaczone z 10 spike-and-recovery w tej samej matrycy, u_rel(R) = 4,7%)
• Masa próbki: m = 100 g = 0,100 kg, u(m) = 0,1 g → u_rel(m) = 0,1%
• Intensywność beta Y-90 (brana pod uwagę przez ε_β): u_rel(ε_β) = 2,5%

Aktywność właściwa: a = A_meas / (R × m) = 15,0 / (0,85 × 0,100) = 176 Bq/kg.

Złożona niepewność względna:
u_rel(a)² = 4,2² + 4,7² + 0,1² + 2,5² = 17,6 + 22,1 + 0,01 + 6,25 = 45,9
u_rel(a) = √45,9 = 6,78%

u(a) = 176 × 0,0678 = 11,9 Bq/kg

Wynik: a = (176 ± 12) Bq/kg Sr-90 (k=1, ~68% CI).

Rozszerzona: (176 ± 24) Bq/kg (k=2, ~95% CI).

Dominujący składnik to odzysk chemiczny (4,7%), nieznacznie przed statystyką zliczeń (4,2%). Poprawa odzysku (lepsza technika, zoptymalizowana procedura) byłaby najefektywniejszym sposobem redukcji niepewności.

Przykład 3: Wykrycie trendu z wykresu kontrolnego próbki kontrolnej

Laboratorium co tydzień mierzy CRM Cs-137 w Marinelli (certyfikowana aktywność = 125 Bq/kg, U = ±8 Bq/kg, k=2).

Wyniki przez 12 tygodni (Bq/kg):
127, 129, 124, 131, 128, 133, 136, 138, 140, 143, 147, 150

Analiza trendu: wyniki rosną systematycznie od ok. 127 do 150 Bq/kg. Różnica = 23 Bq/kg ≈ 18% od wartości nominalnej, i przekracza U = 8 Bq/kg.

Reguła WE nr 8: 8 lub więcej kolejnych punktów po tej samej stronie od średniej (125 Bq/kg). Taki trend jest flagą dla dryftu systematycznego.

Prawdopodobne przyczyny:

• Rosnące napięcie detektora skutkujące przesunięciem piku poza okno analityczne w oprogramowaniu (błąd obliczeniowy sumy w oknie)
• Powoli rozszczelniający się kriostat LN₂ powodujący utratę chłodzenia HPGe i pogorszoną rozdzielczość z reasocjacją tła
• Błąd w kalibracji czasu (timer) skutkujący zawyżaniem czasu pomiaru (zawyżony czas → zawyżona aktywność)

Działanie: przeprowadzić kalibrację energii i FWHM z Am-241/Co-57, sprawdzić poziom LN₂, sprawdzić stabilność timera. Wyniki próbek z tych 12 tygodni muszą być wstrzymane do czasu rozstrzygnięcia przyczyny trendu.

Planowanie eksperymentalne: optymalizacja czasu pomiaru próbki i tła

Jedną z najczęstszych decyzji praktycznych w laboratorium radiometrycznym jest: jak długo mierzyć próbkę, a jak długo tło? Przy ograniczonym czasie aparatury podział czasu między próbkę (t_s) i tło (t_b) ma bezpośredni wpływ na osiągalną niepewność.

Dla przypadku, gdy sygnał netto R_net = R_s − R_b, złożona niepewność:

u(R_net)² = R_s/t_s + R_b/t_b

Minimalizacja u(R_net) przy stałym czasie całkowitym T = t_s + t_b prowadzi do optymalnego podziału:

t_s/t_b = √(R_s/R_b)

Jeżeli tło jest znacznie mniejsze niż sygnał (R_b << R_s), optymalny czas pomiaru tła jest krótszy niż próbki: t_b = t_s × √(R_b/R_s). Natomiast gdy sygnał jest mały w porównaniu z tłem (pomiar niskiej aktywności w dużym tle gamma), optymalny podział czasu jest bardziej wyrównany.

Przykład liczbowy: R_s = 2,0 imp/s, R_b = 0,5 imp/s, T = 3600 s.

Optymalny t_s/t_b = √(2,0/0,5) = √4 = 2,0, więc t_s = 2400 s, t_b = 1200 s.

u(R_net)² = 2,0/2400 + 0,5/1200 = 0,000833 + 0,000417 = 0,00125

u(R_net) = 0,0354 imp/s; u_rel = 0,0354 / (2,0 − 0,5) = 2,36%.

Gdyby mierzono t_s = t_b = 1800 s: u(R_net)² = 2,0/1800 + 0,5/1800 = 0,00139; u_rel = 2,49% — nieznacznie gorszy.

Optymalizacja czasu pomiaru staje się krytyczna przy seriach próbek: 100 próbek razy 30 s różnicy w czasie = prawie godzina przepustowości aparatury na serię.

Systemy zarządzania danymi i śledzenie próbek

W akredytowanym laboratorium radioanalitycznym wynik pomiaru jest uzasadniony tylko wówczas, gdy pełny łańcuch danych od poboru próbki do raportu jest dokumentowany w systemie zarządzania danymi laboratoryjnymi (LIMS — Laboratory Information Management System).

Minimalne wymagania LIMS w laboratorium radiometrycznym:

• Identyfikacja próbki: unikalne ID, data poboru, osoba pobierająca, lokalizacja poboru, rodzaj matrycy.
• Historia przygotowania: daty każdego etapu radiochemicznego, operatorzy, partie odczynników (numer LOT), sprzęt używany (wagi, pipety, autoklawy, kolumny).
• Dane kalibracyjne powiązane z serią: identyfikatory CRM, numery certyfikatów, aktualność kalibracji detektora.
• Wyniki próbek kontrolnych i blankówI: powiązane ze swoją serią analityczną, aby można było prześledzić, czy seria była „pod kontrolą" w momencie pomiaru.
• Zatwierdzenia: podpis analityka i kierownika laboratorium potwierdzający sprawdzenie budżetu niepewności i wyników QC przed wydaniem raportu.

LIMS nie tylko dokumentuje — aktywnie wskazuje niespójności: jeśli CRM wyszedł poza granice wykrywalne, LIMS może zablokować wydanie wyników z danej serii do czasu wyjaśnienia.

W Polsce małe laboratoria WIOŚ często używają arkuszy kalkulacyjnych zamiast pełnego LIMS. To nie jest niezgodne z ISO 17025, jeśli arkusze są walidowane (chronione przed przypadkową edycją, z kontrolą wersji i auditem zmian), ale jest podatne na błędy ludzkie i trudniejsze do audytowania niż zintegrowany system.

Warto pamiętać, że audytorzy PCA podczas oceny na miejscu sprawdzają nie tylko wyniki, ale też integralność ścieżki dokumentacyjnej: czy możliwe jest odtworzenie każdego kroku, który doprowadził do konkretnego liczby w raporcie. Jeśli któregoś etapu brakuje w dokumentacji, wynik jest nominalnie wątpliwy nawet jeśli procedura analityczna była technicznie poprawna. Ta zasada integrity of audit trail odróżnia formalną akredytację od nieformalnego pomiaru "na własne potrzeby". Nowoczesne LIMS-y takie jak LabVantage, STARLIMS czy dedykowane systemy dla radiometrii (np. LabSolutions Radiation) automatyzują tę ścieżkę audytową: każda zmiana rekordu jest datowana i sygnowana, co praktycznie eliminuje ryzyko nieumyślnego skasowania historii analitycznej.

Niepewność próbkowania jako składnik budżetu

Większość artykułów o niepewności pomiaru skupia się na etapie laboratoryjnym. Tymczasem w praktyce środowiskowej istotną składową może być niepewność próbkowania — błąd wynikający z niejednorodności materiału w terenie i niereprezentacyjności pobranej próbki.

Niejednorodność Cs-137 w glebie po awarii Czernobyla jest dobrze udokumentowana: różnica aktywności między próbkami pobranymi w odległości 5 m od siebie może wynosić 20–50% na glebach o zróżnicowanej rzeźbie, roślinności i historii orki. W tym kontekście niepewność laboratoryjna rzędu ±5% jest mniej istotna niż niejednorodność przestrzenna.

Formalnie, całkowita niepewność wyznaczenia aktywności w danym miejscu:

u_total² = u_sampling² + u_preparation² + u_measurement²

EURACHEM/CITAC Guide CG4 ("Quantifying Uncertainty in Analytical Measurement") zaleca, aby dla celów decyzyjnych (np. deklarowanie zgodności z limitami regulacyjnymi) uwzględniać wszystkie trzy składniki i oceniać, który dominuje. Jeżeli u_sampling >> u_measurement, poprawianie procedury laboratoryjnej nie zmniejszy u_total — trzeba poprawić strategię poboru próbek.

W radioanalizie środowiskowej często stosuje się:

• Kompozytowe pobieranie próbek: łączenie kilku podpróbek z tego samego obszaru w jeden pojemnik, by uśrednić niejednorodność przestrzenną przed analizą.
• Duplikaty terenowe (field duplicates): dwie niezależne próbki pobrane z tej samej lokalizacji lub obszaru, które po analizie laboratoryjnej powinny dać podobne wyniki. Różnica między duplikatami charakteryzuje łączny efekt heterogeniczności terenu i powtarzalności poboru.
• Quartowanie (sample splitting): mechaniczne dzielenie dużej próbki (np. 2 kg gleby) na mniejsze porcje przy zachowaniu reprezentatywności frakcji ziarnowych i mineralogicznych.

Pytania otwarte

1. Jakie są granice stosowalności metody spike-and-recovery, gdy domyślamy się, że radionuklid w próbce środowiskowej jest powiązany chemicznie z frakcją ilastą gleby (sorbowany), a spike dodano w formie jonowej z wodnego roztworu? Jak zbadać, czy odzysk spika jest reprezentatywny dla odzysku radionuklidu natywnego?

2. Jak ocenić, czy różnica wyników między dwoma laboratoriami (np. laboratorium A: 210 Bq/kg, laboratorium B: 185 Bq/kg) wynika z błędu systematycznego jednej z metod, czy z rzeczywistego zróżnicowania niejednorodnej próbki? Jakie dodatkowe eksperymenty pozwoliłyby to rozstrzygnąć?

3. Dlaczego próba ślepa terenowa (field blank) jest konieczna przy próbkowaniu wody pitnej w sieci wodociągowej, a nie tylko próba ślepa metodyczna? Jakie zanieczyszczenia może wykryć field blank, a które umknęłyby method blank?

4. Jak interpretować wynik odzysku chemicznego R = 1,12 (112%)? Czy jest to możliwe fizycznie? Jakie błędy mogły spowodować odzysk powyżej 100% i jak je wykryć?

5. Jak poprawnie raportować wyniki z serii, w której część próbek dała wyniki poniżej MDA, a część powyżej? Jakie są zasady obliczania średniej aktywności dla takiej serii (np. metodą MLE dla ocenzurowanych danych)?

6. Jaką rolę odgrywa niejednorodność próbki w budżecie niepewności? Jak scharakteryzować niepewność pobierania próby z heterogenicznych materiałów (np. gleba z kamieniami, żywność z niejednorodnym rozmieszczeniem Cs-137)?

7. Dlaczego sam certyfikat CRM nie gwarantuje spójności pomiarowej laboratorium klienta? Co musi zrobić laboratorium, żeby rzeczywiście uzyskać spójność od BIPM przez NMI do własnego wyniku?

8. Jak formalnie rozróżnić w ISO 17025 termin „verification" i „validation" metody analitycznej? Kiedy wystarczy weryfikacja, a kiedy wymagana jest pełna walidacja, i co konkretnie musi ona obejmować dla metody radioanalitycznej?

Podsumowanie dydaktyczne

1. Wynik radioanalityczny to iloczyn — aktywność końcowa w Bq/kg jest wynikiem mnożenia sygnału detektora przez szereg czynników korygujących: czas, wydajność detekcji, intensywność linii, odzysk chemiczny i masę próbki. Błąd w którymkolwiek z tych czynników przekłada się bezpośrednio na wynik.

2. Kalibracja sprzętu i walidacja metody to różne rzeczy — kalibracja mierzy cechy metryczne detektora (np. ε(E) w danej geometrii). Walidacja metody sprawdza, czy cały łańcuch analityczny (pobór próbki → przygotowanie → pomiar → obliczenia) daje wyniki z wymaganą dokładnością dla konkretnej matrycy i radionuklidu. Oba są konieczne.

3. Blank sample nie jest stratą czasu — każda próba ślepa mierzona równolegle z próbkami właściwymi dostarcza informacji o poziomie zanieczyszczenia tła procesu. Pominięcie blanków pozwala zaoszczędzić czas pomiarowy, ale usuwa mechanizm wczesnego ostrzegania przed błędami systematycznymi.

4. Odzysk chemiczny ≠ wydajność detekcji — to są dwa niezależne czynniki korekcyjne. Odzysk chemiczny opisuje, jaka część analitu przeszła przez procedurę radiochemiczną. Wydajność detekcji opisuje, jaka część promieniowania emitowanego przez preparat daje zliczenia użyteczne. Mylenie tych pojęć prowadzi do błędów korygowania wyniku.

5. MDA jest parametrem układu, nie stałą materiałową — granica wykrywalności dla danego radionuklidu zależy od tła detektora, geometrii, czasu pomiaru, masy próbki i wydajności. Można ją obniżyć zwiększając masę próbki, wydłużając czas pomiaru lub osłaniając detektor od tła. Jest błędem informacyjnym twierdzić, że dany izotop „nie jest wykrywalny" bez podania warunków pomiaru.

6. Wykresy kontrolne chronią całą serię — regularne pomiary CRM i śledzenie wyników na kartach Shewhart'a/Levey-Jennings'a pozwalają wykryć dryft lub skok systematyczny. Bez tych narzędzi seria kilkudziesięciu próbek analitycznych może być obciążona błędem, który ujawni się dopiero w teście biegłości, kiedy przekorygowanie jest kosztowne lub niemożliwe.

7. Spike-and-recovery jest narzędziem walidacji, nie tylko kontroli — regularne spike-and-recovery w każdej serii analitycznej to narzędzie walidacyjne: sprawdza powtarzalność i poprawność procedury w bieżących warunkach. Jednorazowy test spike przy wdrożeniu metody nie wystarczy, bo warunki laboratorium (partia odczynnika, skład matrycy, stan kolumny) zmieniają się z czasem.

8. Traceability to infrastruktura porównywalności — wyniki z polskich laboratoriów środowiskowych trafiają do europejskich baz danych (EURDEP, IAEA) i są porównywane z wynikami z Francji, Finlandii czy Japonii. Porównywalność ta jest możliwa tylko wtedy, gdy każde laboratorium utrzymuje udokumentowaną spójność pomiarową ze wspólnymi wzorcami wyższego rzędu. To jest metrologia jako infrastruktura naukowa — niewidoczna, ale konieczna.

Dodatkowe materiały multimedialne

Najkrótsze podsumowanie: wynik radioanalityczny jest łańcuchem. Detektor daje zliczenia, ale masa, objętość, tło, wydajność i odzysk decydują, czy zliczenia zamienią się w wiarygodne Bq/kg.

Powiązane kalkulatory i narzędzia

• Kalkulator: Odzysk chemiczny — Korekta aktywności o odzysk znacznika oraz prosta propagacja niepewności względnej.

Ćwiczenia praktyczne

Pierwsze ćwiczenie rachunkowe: student dostaje cztery syntetyczne pomiary masy dostarczonej wody dla pipety nominalnie 1,000 mL. Ma policzyć średnią masę, objętość przy zadanej gęstości, błąd względny i odchylenie standardowe średniej.

Drugie ćwiczenie: porównać pipetę jednomiarową i automatyczną na dwóch syntetycznych seriach. Jedna seria ma mały rozrzut, ale średnią przesuniętą; druga ma średnią bliską wartości nominalnej, ale większy rozrzut. Student ma wskazać, która ma problem poprawności, a która powtarzalności.

Trzecie ćwiczenie: z danych N_sample, t_sample, N_blank, t_blank, epsilon, p_gamma, m i R_chem obliczyć aktywność właściwą w Bq/kg. Następnie powtórzyć obliczenie przy R_chem = 1, aby zobaczyć błąd pominięcia odzysku.

Czwarte ćwiczenie: przygotować tabelę budżetu niepewności dla wyniku radioanalizy. Student ma podać względne niepewności zliczeń, wydajności, masy, intensywności linii i odzysku, a następnie złożyć je kwadratowo.

Piąte ćwiczenie audytowe: przeczytać opis wyniku Bq/kg i sprawdzić, czy zawiera masę próbki, czas pomiaru, tło, wydajność, próbę ślepą, odzysk chemiczny i niepewność. Brakujące pola należy oznaczyć jako ryzyko jakościowe.

Szóste ćwiczenie kontrolne: dla serii wyników próbki kontrolnej narysować prostą kartę kontroli. Student ma zaznaczyć średnią, granice ostrzegawcze i jeden punkt odstający, a potem opisać, czy problem wygląda na losowy, systematyczny czy związany z błędem opisu próbki.

Przejdź do ćwiczenia interaktywnego

Powiązane artykuły

• Kalibracja i niepewność pomiaru w laboratorium jądrowym
• Statystyka zliczeń promieniotwórczych: Poisson, Gauss, chi-kwadrat i błędy aparaturowe
• Spektrometria gamma w praktyce: kalibracja energii, rozdzielczość i wydajność detektora
• Naturalna promieniotwórczość żywności i materiałów budowlanych: od widma do Bq/kg
• Wielkości dozymetryczne: dawka pochłonięta, równoważna, efektywna i operacyjna