Streszczenie

Samopochłanianie oznacza, że promieniowanie gamma powstaje wewnątrz próbki, ale nie wszystkie fotony opuszczają ją z takim samym prawdopodobieństwem. Część energii zostaje rozproszona albo osłabiona w samej próbce, zanim foton dotrze do detektora, więc aktywność obliczona z fotopiku może być systematycznie zaniżona.1,2

Ten artykuł pokazuje samopochłanianie jako korektę metrologiczną: zależną od energii gamma, masy powierzchniowej, geometrii naczynia, gęstości, składu matrycy i tego, czy wzorzec kalibracyjny naprawdę przypomina próbkę badaną.

Rozszerzenie tematu

Próbka nie jest punktem

W najprostszym rysunku źródło promieniowania jest punktem przed detektorem. W rzeczywistym pomiarze środowiskowym próbka bywa objętością: pojemnikiem Marinelli, cylindrem z proszkiem, słoikiem z glebą, próbką materiału budowlanego albo naczyniem z solą. Radionuklid nie leży wtedy na powierzchni. Jest rozmieszczony w objętości.

Foton gamma wyemitowany tuż przy ściance skierowanej do detektora ma krótką drogę przez próbkę. Foton wyemitowany głębiej musi przejść przez większą masę materiału. Po drodze może zostać rozproszony komptonowsko albo pochłonięty fotoelektrycznie. Nawet jeśli nie znika całkowicie, może opuścić fotopik i trafić do tła ciągłego.

Dlatego aktywność z fotopiku nie zależy tylko od liczby rozpadów. Zależy też od tego, ile fotonów danej energii zdołało wyjść z próbki w kierunku detektora i zachować energię pełnej absorpcji.

Różnica między absorpcją zewnętrzną a samopochłanianiem

W klasycznym ćwiczeniu z osłabiania promieniowania mamy zewnętrzną wiązkę i warstwę absorbentu. Model jest prosty:

I = I0 * exp(-µx).

W samopochłanianiu źródło znajduje się w absorbencie. Foton może powstać w różnych miejscach próbki, pod różnymi kątami i z różną drogą do powierzchni. To bardziej złożona geometria. Mimo to ta sama fizyka osłabiania jest obecna: decyduje współczynnik osłabienia i masa materiału na drodze fotonu.3

Dla uproszczonej warstwy, w której emisja jest równomierna, można otrzymać współczynnik korekcyjny:

f = (1 - exp(-µx)) / (µx).

W przewodniku KChRS taki wzór jest użyty w ćwiczeniu o samopochłanianiu gamma w próbce KCl, gdzie µ jest masowym współczynnikiem osłabienia, a x masą powierzchniową próbki.1

Co oznacza współczynnik f

Współczynnik f opisuje, jaka część emisji jest widoczna w porównaniu z przypadkiem braku samopochłaniania w danym uproszczonym modelu. Jeżeli µx jest bardzo małe, rozwinięcie wykładnicze daje:

f ≈ 1 - µx/2.

To znaczy, że cienka albo słabo pochłaniająca próbka prawie nie wymaga poprawki. Jeżeli µx rośnie, f maleje: coraz większa część fotonów nie wychodzi z próbki w postaci użytecznej dla fotopiku.

W rachunku aktywności pojawia się to jako poprawka wydajności. Jeśli pole piku jest mniejsze z powodu samopochłaniania, a my tego nie uwzględnimy, uznamy, że próbka ma mniejszą aktywność niż w rzeczywistości. To błąd systematyczny, nie fluktuacja statystyczna.

Dlaczego przykład KCl jest dydaktycznie dobry

KCl zawiera naturalny potas, a naturalny potas zawiera niewielką frakcję K-40. K-40 emituje linię gamma około 1461 keV, użyteczną w spektrometrii gamma. Przewodnik KChRS wykorzystuje ten przykład do pokazania samopochłaniania w objętościowej próbce soli.1

Ten przykład jest wartościowy z trzech powodów:

  • radionuklid jest naturalny i równomiernie związany z masą próbki,
  • energia 1461 keV jest typową energią gamma używaną w analizie środowiskowej,
  • można pokazać różnicę między aktywnością obliczoną z masy i składu a aktywnością widzianą przez detektor.

W artykule nie chodzi jednak o instrukcję przygotowania stanowiska. Chodzi o pojęcie: nawet dobrze znana próbka może dać mniejszy fotopik, jeśli część fotonów jest tracona w samej próbce.

Masa powierzchniowa zamiast samej masy

W samopochłanianiu ważna jest nie tylko masa próbki, ale masa na drodze promieniowania. Dlatego pojawia się pojęcie masy powierzchniowej:

x = m / S,

gdzie m jest masą próbki, a S efektywną powierzchnią, przez którą promieniowanie „wychodzi” w danym uproszczeniu. Jednostką jest zwykle g/cm^2.

Dwie próbki mogą mieć tę samą masę, ale inną geometrię. Cienka, szeroka warstwa ma małą masę powierzchniową. Wysoki, wąski cylinder może mieć większą drogę przez materiał. W spektrometrii gamma naczynie Marinelli jest popularne właśnie dlatego, że poprawia geometrię względem detektora, ale nadal nie usuwa problemu samopochłaniania dla każdej matrycy.

Masa powierzchniowa jest pomostem między intuicją „ile materiału mam?” a fizyką „ile materiału widzi foton po drodze?”.

Energia gamma i skład matrycy

Samopochłanianie silnie zależy od energii. Dla wyższych energii gamma próbka bywa bardziej przezroczysta. Dla niższych energii rośnie znaczenie absorpcji fotoelektrycznej, zwłaszcza w matrycach zawierających pierwiastki o większej liczbie atomowej.

To oznacza, że poprawka dla K-40 przy 1461 keV nie może być automatycznie przeniesiona na niskoenergetyczną linię innego radionuklidu. W tej samej próbce jedna linia może być osłabiana słabo, a druga silnie. Dlatego krzywa wydajności detektora dla geometrii objętościowej powinna uwzględniać energię.

Skład chemiczny też ma znaczenie. Woda, sól, gleba, beton, popiół i ruda nie mają tych samych współczynników osłabienia. Gęstość objętościowa, wilgotność i uziarnienie mogą zmienić efektywną matrycę pomiarową.

Kiedy osobna poprawka nie jest potrzebna

Przewodnik KChRS podaje ważną praktyczną zasadę: poprawka na samopochłanianie nie jest wymagana, jeśli kalibracja wydajnościowa detektora była wykonana dla próbek w tej samej konfiguracji geometrycznej, o podobnej gęstości i podobnym współczynniku osłabienia promieniowania gamma.1

To zdanie warto zapamiętać. Najlepszą poprawką jest często dobra kalibracja matrycowa. Jeżeli wzorzec i próbka mają:

  • ten sam typ naczynia,
  • podobną objętość wypełnienia,
  • podobną gęstość,
  • podobny skład osłabiający,
  • podobną energię analizowanej linii,

to wydajność wyznaczona ze wzorca może już zawierać efekt samopochłaniania. Wtedy dokładanie osobnej poprawki mogłoby nawet podwójnie skorygować ten sam efekt.

Kiedy poprawka staje się krytyczna

Samopochłanianie staje się szczególnie ważne, gdy:

  • próbka ma dużą objętość,
  • próbka jest gęsta,
  • analizowana energia gamma jest niska,
  • matryca zawiera pierwiastki silnie pochłaniające,
  • wzorzec ma inną geometrię lub gęstość niż próbka,
  • próbka jest niejednorodna,
  • wynik ma być ilościowy, a nie tylko jakościowy.

W pomiarach środowiskowych dotyczy to gleb, popiołów, materiałów budowlanych, rud, osadów i próbek o zmiennej wilgotności. W takich przypadkach sam fotopik nie wystarcza. Trzeba wiedzieć, ile aktywności mogło „schować się” w objętości próbki.

Jak samopochłanianie wpływa na fotopik

Samopochłanianie może działać na dwa sposoby. Po pierwsze, foton może zostać pochłonięty i nie dotrzeć do detektora. Po drugie, może zostać rozproszony i dotrzeć z inną energią albo z innym kierunkiem. W obu przypadkach liczba zliczeń w fotopiku pełnej absorpcji maleje.

To ważne, bo aktywność z gamma liczy się zwykle z pola fotopiku:

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

Jeżeli epsilon nie zawiera samopochłaniania, trzeba uwzględnić f:

A = N_net / (t * epsilon_0 * f * I_gamma * m).

epsilon_0 oznacza tu wydajność bez samopochłaniania w próbce, a f korektę widoczności fotonów. W praktyce laboratorium może stosować inne konwencje zapisu, ale sens jest ten sam: im mniejsze f, tym większa aktywność odpowiada temu samemu polu piku.

Samopochłanianie a tło

Samopochłanianie nie jest tym samym co tło. Tło jest sygnałem niezwiązanym z analizowanym fotopikiem próbki albo pochodzącym z otoczenia i aparatury. Samopochłanianie jest stratą części sygnału próbki.

Oba efekty spotykają się w analizie. Jeżeli pik jest słaby, a tło duże, pole netto ma większą niepewność. Jeżeli dodatkowo samopochłanianie jest silne, pole netto jest mniejsze, więc względna niepewność rośnie. Wtedy wydłużenie czasu pomiaru może zmniejszyć niepewność statystyczną, ale nie usunie błędu systematycznego źle dobranej korekty.

Dlatego raport powinien oddzielać:

  • odjęcie tła,
  • korektę wydajności,
  • korektę samopochłaniania,
  • niepewność statystyczną,
  • niepewność modelu lub wzorca matrycowego.

Naczynie Marinelli i geometria objętościowa

Naczynie Marinelli otacza detektor próbką, zwiększając kąt bryłowy i wydajność dla słabych aktywności. Jest bardzo użyteczne w spektrometrii środowiskowej. Jednocześnie wprowadza geometrię objętościową: fotony powstają wokół detektora w różnych odległościach i mają różne drogi przez próbkę.

To oznacza, że prosty wzór warstwowy f = (1 - exp(-µx))/(µx) jest dydaktycznym modelem, nie pełnym opisem każdego naczynia. W profesjonalnej analizie można używać wzorców w tej samej geometrii, symulacji Monte Carlo albo empirycznych poprawek zależnych od energii i matrycy.

Najważniejsze dla czytelnika jest rozpoznanie pułapki: większa próbka zwykle daje więcej aktywności całkowitej, ale niekoniecznie proporcjonalnie większy sygnał w fotopiku, bo rośnie droga samopochłaniania.

Błąd typowy: jedna wydajność dla wszystkich próbek

Częsty błąd polega na przyjęciu jednej wydajności detektora dla wszystkich próbek, bo „mierzono tym samym spektrometrem”. Tymczasem wydajność w fotopiku zależy od energii i geometrii, a przy próbkach objętościowych także od matrycy.

Przykład: kalibracja wykonana na wodzie w naczyniu Marinelli nie musi idealnie pasować do suchej gleby, popiołu albo soli. Nawet jeśli objętość jest taka sama, gęstość i współczynnik osłabienia mogą być inne. Przy niskoenergetycznych liniach różnice będą większe niż przy 1461 keV.

Wynik liczbowy bez informacji o matrycy i geometrii powinien być traktowany ostrożnie. Może być dobry jako screening, ale nie jako dokładna aktywność właściwa.

Minimalny opis korekty

Jeżeli artykuł albo raport podaje aktywność z próbki objętościowej, powinien powiedzieć:

  • jaki był typ detektora,
  • jaka energia gamma była analizowana,
  • jaki był kształt i rozmiar próbki,
  • jaka była masa, gęstość albo masa powierzchniowa,
  • czy wzorzec miał tę samą matrycę,
  • czy wydajność zawierała samopochłanianie,
  • czy stosowano osobny współczynnik f,
  • jaka niepewność wynika z tej poprawki.

Bez tego czytelnik nie wie, czy pole piku zostało przeliczone na aktywność w sposób porównywalny z innymi pomiarami.

Modele obliczeniowe samopochłaniania

Uproszczony wzór f = (1 − e^(−µx))/(µx) działa dobrze dla cienkiej, jednorodnej warstwy emitującej prostopadle do detektora. W praktyce laboratoryjnej próbka jest trójwymiarowa (cylindryczna, Marinelli), detektor ma skończoną geometrię, a fotony wychodzą pod różnymi kątami. W takich przypadkach konieczne są bardziej rozbudowane modele.

Symulacja Monte Carlo jest metodą referencyjną dla złożonych geometrii. Fotony generowane są losowo w całej objętości próbki z kierunkami i energiami zgodnymi ze spektrum emisji. Każdy foton jest śledzony przez materiał próbki, pojemnik i detector — z uwzględnieniem wszystkich procesów oddziaływania (fotoelektryczny, Compton, produkcja par). Kody: MCNP (Los Alamos), GEANT4 (CERN), EGSnrc (NRC Kanada), FLUKA. Precyzja: 1–3% przy dostatecznej ilości historii. Wadą jest konieczność precyzyjnego zdefiniowania geometrii i składu chemicznego.

LabSOCS (Laboratory Sourceless Object Calibration Software, Mirion/Canberra) jest programem opartym na parametryzowanych symulacjach Monte Carlo, który oblicza krzywe wydajności dla standardowych geometrii próbek (cylindry, Marinelli, ampułki Petri, całe ciało) bez potrzeby posiadania wzorców. Użytkownik podaje geometrię próbki, skład chemiczny i gęstość. LabSOCS generuje krzywą ε(E) z uwzględnieniem samopochłaniania. Weryfikacja LabSOCS przez wzorce certyfikowane jest wymagana w akredytowanych laboratoriach.

EFFTRAN (EFFiciency TRANsfer code) jest alternatywą akademicką obliczającą przeskalowanie wydajności między różnymi geometriami na podstawie współczynników osłabienia. Metodologia znana jako „efficiency transfer" jest opisana w IAEA TECDOC-1360 (2003) i jest podstawą akcelerowanej kalibracji bez wzorców fizycznych dla każdej nowej geometrii.

Empiryczna korekcja przez transmisję (transmission method): próbka mierzona jest bez dodatkowego radionuklidu i z punktowym źródłem wzorcowym umieszczonym za próbką w osi detektora. Stosunek sygnału ze źródłem do sygnału bez próbki daje transmisję T = e^(−µx). Stąd µx jest znane, a f obliczane z wzoru warstwowego. Metoda prosta, ale zakłada jednorodny model.

Współczynniki masowego osłabienia dla typowych matryc

Znajomość µ (masowy współczynnik osłabienia, cm²/g) dla danej energii i składu matrycy jest podstawą korekty. Wartości z bazy danych NIST XCOM (online: physics.nist.gov/xcom):

Matryca Gęstość (g/cm³) µ przy 60 keV (cm²/g) µ przy 662 keV (cm²/g) µ przy 1461 keV (cm²/g)
Woda (H₂O) 1,00 0,200 0,0857 0,0559
SiO₂ (kwarc) 2,65 0,268 0,0861 0,0561
Gleba (typowa) 1,3–1,6 ~0,27 ~0,087 ~0,056
CaCO₃ (wapień) 2,71 0,332 0,0883 0,0563
NaCl (sól) 2,16 0,456 0,0967 0,0581
Popiół lotny 0,7–1,2 ~0,26 ~0,086 ~0,056
Fe₂O₃ (ruda) 5,24 0,814 0,0985 0,0591
Beton 2,3 0,26 0,086 0,056
Biota (tkanki) ~1,0 ~0,20 ~0,086 ~0,056

Kluczowa obserwacja: µ przy 662 keV i 1461 keV różni się między wodą a rudą żelaza mniej niż 15%, co oznacza, że dla wysokoenergetycznych linii efekt matrycy (samopochłanianie) jest względnie mały i kalibracja wodna daje akceptowalną dokładność dla lekkiej matrycy. Przy 60 keV różnica µ między wodą a rudą żelaza wynosi czynnik 4 — w takim przypadku kalibracja musi być matrycowo dopasowana.

Samopochłanianie w zakresie niskich energii

Energie poniżej 100–200 keV są szczególnie trudne dla spektrometrii gamma objętościowej z powodu silnego samopochłaniania. Ważne radionuklidy w tym zakresie:

  • Pb-210 (E_gamma = 46,5 keV): używany jako marker wieku osadów; silna absorpcja w gęstych osadach morskich (µ/ρ ~1 cm²/g);
  • Th-234 (E_gamma = 63,3 keV, 92,6 keV): produkt rozpadu U-238, wskaźnik nierównowagi;
  • U-235 (E_gamma = 143,8 i 185,7 keV): kluczowe linie do oceny wzbogacenia;
  • Am-241 (E_gamma = 59,5 keV): wskaźnik skażeń transuranic, ale przy 59,5 keV silna absorpcja;
  • Ba-133 (E_gamma = 53,2, 80,9, 160,6 keV): wzorzec kalibracyjny niskoenergetyczny;
  • Eu-152 (E_gamma = 121,8, 344,3, 778,9 keV): wzorzec wielopunktowy z linią niskoenergetyczną.

Dla Pb-210 w próbce gleby o masie 50 g w naczyniu cylindrycznym 50 mL (typowy parametr laboratorium morskiego) samopochłanianie może wynosić 60–80%, zależnie od gęstości. Zignorowanie tego oznacza zaniżenie aktywności Pb-210 o współczynnik 3–5. To bezpośrednio błędna ocena tempa akumulacji osadów w badaniach datowania przez Pb-210 (metoda CRS/CIC).

Metoda kompensacji: wzorce kalibracyjne niskoenergetyczne (Ba-133, Am-241, Co-57) muszą być przygotowane w tej samej geometrii i identycznym składzie matrycy co próbka. Alternatywnie: zastosowanie symulacji Monte Carlo (LabSOCS, MCNP) z pełnym opisem składu matrycy.

Wpływ wilgotności próbki na samopochłanianie

Wilgotność gleby, popiołu lub osadu zmienia gęstość efektywną i skład matrycy. Suchą glebę i glebę wilgotną opisuje inny efektywny współczynnik osłabienia.

Dla gleby o wilgotności w%: gęstość efektywna ρ_eff ≈ ρ_sucha × (1 + w/100). Efektywny µ_eff jest ważoną sumą wkładów gleby i wody:

µ_eff = (µ_gleba × m_gleba + µ_woda × m_woda) / (m_gleba + m_woda)

Przy wilgotności 30% (typowe dla świeżej gleby) i energii 662 keV: różnica µ_eff między suchą (µ=0,0862) a wilgotną glebą (µ_eff ~0,0862 × 0,77 + 0,0857 × 0,30 = ok. 0,0861) jest mała — zaledwie 0,1% przy 662 keV. Natomiast przy 60 keV: µ_gleba ~0,27, µ_woda ~0,200; przy 30% wilgotności µ_eff ~0,251 vs. 0,270 dla suchej gleby — różnica 7%.

W praktyce: próbki glebowe wysuszają się w 105°C i waży się suchą masę. Aktywność podaje się na kg suchej masy (Bq/kg d.m.). Używanie masy mokrej bez normalizacji jest błędem systematycznym.

Samopochłanianie w próbkach biologicznych i żywności

Próbki żywności (ryby, mleko, zioła, mięso) i tkanki biologiczne mają gęstość zbliżoną do wody i lekki skład (C, H, O, N), co czyni je stosunkowo przezroczystymi dla energii gamma > 200 keV. Kalibracja wodna w tej samej geometrii jest zazwyczaj akceptowalna z błędem < 5%.

Trudniejsze przypadki w żywności:

  • Suche grzyby (niski współczynnik wypełnienia, puste komórki powietrzne) — gęstość 0,1–0,3 g/cm³;
  • Suchy herbata (podobny problem);
  • Sól spożywcza (NaCl) — gęstość 2,16 g/cm³ vs. 1,0 g/cm³ dla wody, przy 60 keV µ_NaCl ≈ 2,2 × µ_woda.

Dyrektywa UE o monitoringu żywności po Fukushimie (Rozporządzenie 2016/52) i normy ISO 18589 (gleba), ISO 10703 (woda) wymagają walidacji metody pomiarowej z wyraźną informacją o kalibracji matrycowej lub obliczeniowej korekcie samopochłaniania.

W Polsce CLOR prowadzi akredytowane oznaczenia aktywności Cs-137, Sr-90, Am-241 w żywności i środowisku. Ich procedury SOP dla spektrometrii gamma opisują kalibrację matrycową i dokumentowanie poprawki samopochłaniania.

Kalibracja wydajności a samopochłanianie — wzorce matrycowe

Wzorce matrycowe (matrix-matched standards) są podstawową metodą radzenia sobie z samopochłanianiem bez obliczeniowej korekty. Zasada jest prosta: wzorzec i próbka mają ten sam skład chemiczny i tę samą geometrię — wtedy fotony tracą identyczną część energii w obu przypadkach i wydajność εa z wzorca jest już poprawna dla próbki.

Wzorce matrycowe są wytwarzane na zamówienie lub kupowane jako certyfikowane materiały odniesienia (CRM). Ważne instytucje wydające CRM dla spektrometrii gamma środowiskowej:

  • NIST (USA): SRM 4350B (River Sediment), SRM 4353 (Rocky Flats soil), SRM 4354 (Freshwater lake sediment) — certyfikowane aktywności Cs-137, Pb-210, Ra-226, Am-241 w matrycach osadowych;
  • IRMM/JRC Geel (Belgia): IAEA-soil-6, IAEA-375, IAEA-447 — gleby i osady z certyfikowanymi radionuklidami;
  • IAEA: cykliczne programy ALMERA dostarczają próbek do PT z deklarowaną aktywnością — nie są to jednak CRM, lecz próbki PT (proficiency testing).

Przygotowanie wzorca matrycowego in-house: do pustej matrycy (odpopielona lub zsyntetyzowana gleba, kreda, sól) dodaje się certyfikowany roztwór radionuklidu, miesza, suszy i homogenizuje. Wynik musi być walidowany przez porównanie z zewnętrznym CRM lub PT. Niejednorodność wzorca własnego jest głównym ryzykiem.

Korekcja przez build-up factor

Przy grubych próbkach i wyższych energiach nie cały rozproszony foton jest stracony. Część fotonów Comptona, które zmieniły kierunek, może nadal trafić do detektora i zdeponować energię poza fotopikiem (do tła Comptonowskiego). Dla bardzo grubej geometrii może się zdarzyć, że promieniowanie rozproszone wielokrotnie i w końcu pochłonięte wewnątrz próbki „nie dotarło" do detektora — model wykładniczy e^(−µx) daje wtedy za małą transmisję przez pomijanie promieniowania odbitego od sąsiednich obszarów.

Fizycznie poprawny opis wymaga uwzględnienia buildup factor B, który z definicji jest stosunkiem całego strumienia fotonów w punkcie do strumienia promieniowania nierozproszonego:

I_total = I₀ × B × e^(−µx)

Dla spektrometrii gamma środowiskowej i typowych geometrii (masy 100–500 g, energie 200–1500 keV) buildup factor jest zazwyczaj nieistotny — korektę samopochłaniania obejmuje wzór f bez buildup. Buildup staje się ważny dla grubych ekranów i obliczeń osłonowych, nie dla standardowych naczyń laboratoryjnych.

Samopochłanianie a certyfikacja laboratoriów

Wymagania normy ISO/IEC 17025:2017 nie wymieniają wprost korekty samopochłaniania, ale §7.2.1 wymaga walidacji metody, §7.6.1 nakazuje identyfikację składowych niepewności, a §7.2.2 wymaga sprawdzania spójności metod. Korekta samopochłaniania jest składową niepewności i musi być opisana w procedurze SOP (standard operating procedure).

Audytor techniczny PCA podczas inspekcji laboratorium spektrometrii gamma sprawdza: czy wzorzec kalibracyjny ma ten sam rodzaj geometrii i porównywalną gęstość/matrycę co typowe próbki, lub czy laboratorium ma zatwierdzony model obliczeniowy (LabSOCS, Monte Carlo) z walidacją przez CRM. Brak dokumentacji korekty samopochłaniania jest niezgodnością podczas audytu.

Zakres akredytacji PCA dla spektrometrii gamma zwykle obejmuje: matryce gleba, osad, biota, żywność, woda; energie 100–2000 keV; nuklidy Cs-137, Co-60, K-40, Ra-226, Pb-210, Am-241. Dla każdej kombinacji matryca×energia laboratorium powinno wykazać, że samopochłanianie jest albo nieistotne (małe próbki, wysokie energie), albo skorigowane z udokumentowaną niepewnością.

Tabela: kiedy samopochłanianie jest krytyczne

Poniżej ocena konieczności korekty samopochłaniania dla typowych kombinacji matryca × energia:

Matryca / Energia 46 keV (Pb-210) 59 keV (Am-241) 185 keV (U-235) 662 keV (Cs-137) 1461 keV (K-40)
Woda (100 mL) Tak (15–30%) Tak (10–25%) Umiarkowana (5%) Mała (2–3%) Mała (1%)
Gleba (200 g) Krytyczna (>40%) Krytyczna (>30%) Tak (8–12%) Umiarkowana (4–8%) Mała (2–4%)
NaCl (100 g) Krytyczna (>50%) Krytyczna (>40%) Tak (10–15%) Umiarkowana (5–8%) Mała (3–5%)
Biota (100 g) Tak (15–30%) Tak (10–20%) Umiarkowana (5%) Mała (2–4%) Mała (1–2%)
Popiół lotny (50 g) Krytyczna (>35%) Krytyczna (>25%) Tak (6–10%) Mała (3–5%) Mała (1–3%)

Ocena „Krytyczna" oznacza, że pominięcie korekty prowadzi do błędu >10% na aktywności — co jest zazwyczaj nieakceptowalne w akredytowanych pomiarach. „Umiarkowana" = 3–10%, „Mała" = <3%.

Metoda transmisji — praktyczne wyznaczanie korekty

Metoda transmisji pozwala zmierzyć eksperymentalnie µx próbki bez modelu obliczeniowego. Procedura dydaktyczna:

  1. Mierzy się tło spektrometru (czas t₀);
  2. Umieszcza się źródło wzorcowe punktowe (np. Cs-137) na osi detektora w odległości d, mierzy N₀ zliczeń przez czas t;
  3. Umieszcza się próbkę między źródłem a detektorem (ta sama geometria), mierzy N_t zliczeń przez czas t;
  4. Transmisja: T = (N_t − tło) / (N₀ − tło);
  5. Z T = e^(−µx): µx = −ln(T);
  6. Korekta f = (1 − e^(−µx)) / (µx) = (1 − T) / (−ln T).

Ograniczenia: metoda zakłada jednorodność próbki, działa najlepiej dla geometrii warstwowej, zakłada brak rozproszenia do detektora (tzw. build-up). Dla niskich energii i grubych próbek foton rozproszony może fałszywie zwiększyć transmisję. W takich przypadkach należy stosować wąskie kolimatory i ołowiane ograniczniki geometrii.

Samopochłanianie w kalibracji bez wzorców fizycznych

Alternatywą dla wzorców CRM jest tzw. kalibracja obliczeniowa (sourceless calibration), w której wydajność detektora ε(E) jest obliczana numerycznie bez wykonywania pomiarów na fizycznych wzorcach radionuklidowych. Metody:

LabSOCS / ISOCS (Mirion Technologies): oprogramowanie parametryzuje odpowiedź detektora HPGe na podstawie charakterystyk fabrycznych (współrzędne geometryczne kryształu HPGe) zmierzonych przez producenta. Użytkownik definiuje geometrię próbki (typ naczynia, wymiary, skład i gęstość matrycy) i oprogramowanie oblicza ε(E). Niepewność metody deklarowana przez producenta: 2–5% dla geometrii Marinelli, 5–10% dla złożonych geometrii. Walidacja przez CRM jest wymagana przy wdrożeniu i przy każdej zmianie geometrii.

ANGLE (Angle Geometry calculation for radioactivity measurements, IAEA/Institute Vinča): kod podobny do LabSOCS, tańszy, dostępny dla krajów rozwijających się przez IAEA/RCA. Stosowany w polskich uczelniach i instytutach.

Kalibracja obliczeniowa jest szczególnie użyteczna gdy: (a) masz wiele geometrii i nie możesz zakupić CRM dla każdej; (b) geometria próbki jest niestandardowa (fragmenty skały, obiekt biologiczny); (c) chcesz szybko ocenić, czy samopochłanianie wymaga pełnej korekcji. Nie zastępuje całkowicie pomiarów sprawdzających przez CRM.

Połączenie LabSOCS z GENIE 2000 (oprogramowanie analizy widm Mirion) umożliwia automatyczne raportowanie aktywności z korekcją samopochłaniania, zintegrowaną krzywą wydajności i kompletnym budżetem niepewności — w formacie wymaganym przez normy ISO i akredytację PCA. W Polsce uczelnie techniczne (AGH, PW) coraz częściej wdrażają LabSOCS lub ANGLE jako podstawę kalibracji laboratoriów pomiarowych, co eliminuje konieczność utrzymywania kosztownego zbioru fizycznych wzorców kalibracyjnych w każdej geometrii naczynia Marinelli i cylindra pomiarowego.

Budżet niepewności korekty samopochłaniania

Niepewność korekty samopochłaniania pochodzi z kilku źródeł:

Niepewność µ: wartości z tabel NIST mają typowo 1–3% dla energii powyżej 100 keV i 3–5% poniżej. Dla matryc złożonych niepewność składu chemicznego (zwłaszcza zawartości żelaza, wapnia, baru) przekłada się na niepewność µ_eff.

Niepewność x (masy powierzchniowej): pochodzi z ważenia próbki, znajomości geometrii pojemnika i jednorodności rozkładu próbki. Przy niepewności masy 0,5% i geometrii 1% — łączna niepewność x może wynosić 1–2%.

Niejednorodność próbki: frakcjonowanie ziarnowe (grube i drobne frakcje mają inne gęstości), segregacja składników (np. gęste minerały na dnie), pęcherze powietrzne. Niejednorodność jest trudna do kwantyfikacji i może wprowadzać niesymetryczne błędy. Standardem jest jednorodne mielenie i mieszanie próbki przed ważeniem.

Model geometryczny: wzór f dla warstwy jest przybliżeniem. Dla geometrii Marinelli, cylindra lub ampułki Petri wzór może być nieadekwatny. Niepewność modelu: 2–10% zależnie od energii i geometrii.

Łączna niepewność korekty samopochłaniania: 3–15%, dominująca przy niskich energiach i złożonych geometriach. Dlatego dla energii poniżej 100 keV zalecane jest stosowanie Monte Carlo lub kalibracji matrycowej, a nie prostego wzoru f.

Normy i wytyczne MAEA dotyczące korekty samopochłaniania

IAEA TECDOC-1360 (2003) „Standardization of Radioactive Waste Packages" i IAEA Safety Series 295 (2019) „Radiation Protection and Safety in Medical Uses of Ionizing Radiation" zawierają szczegółowe wytyczne dotyczące korekcji samopochłaniania.

Dla spektrometrii gamma próbek środowiskowych kluczowe dokumenty MAEA:

  • IAEA TECDOC-619 (1991) „Measurement of Radionuclides in Food and the Environment": rozdział 3 opisuje kalibrację wydajności i korekcję samopochłaniania dla próbek objętościowych;
  • IAEA TECDOC-1215 (2001) „Analytical Quality in Nuclear Applications": wymagania dotyczące kalibracji matrycowej i dokumentacji korekty;
  • ISO 20042:2019 „Measurement of radioactivity in the environment — gamma-ray spectrometry": definiuje wymagania dla wydajności w fotopiku, korekcji samopochłaniania i niepewności metody.

W Polsce wymagania wynikają z akredytacji PCA (według EA-4/16 i ISO 17025). Laboratoria muszą walidować metodę korekcji samopochłaniania i dokumentować źródła niepewności. PCA przeprowadza audyty techniczne obejmujące ocenę tego aspektu.

Trzy przykłady rachunkowe

Przykład 1 — korekta f dla próbki KCl

Próbka: 100 g KCl w cylindrycznym naczyniu o przekroju S = 50 cm². Masa powierzchniowa: x = m/S = 100 g / 50 cm² = 2 g/cm². Skład KCl: K (39,1%), Cl (60,9%). Masowy współczynnik osłabienia µ dla KCl przy E_gamma = 1461 keV (linia K-40): µ ≈ 0,0581 cm²/g (z NIST XCOM). Parametr: µx = 0,0581 × 2,0 = 0,1162.

Współczynnik korekcyjny:
f = (1 − e^(−0,1162)) / 0,1162 = (1 − 0,8902) / 0,1162 = 0,1098 / 0,1162 = 0,945

Interpretacja: 94,5% emisji jest widoczne w fotopiku. Samopochłanianie wynosi 5,5%. Jeśli pominięto by tę korektę, aktywność K-40 byłaby zaniżona o 5,5%. Dla K-40 w glebach i materiałach budowlanych typowy x = 1–5 g/cm² i µ = 0,056 cm²/g daje f = 0,97–0,87 — zakres korekcji 3–13%.

Przykład 2 — samopochłanianie niskiej energii (Pb-210)

Próbka osadu morskiego: 20 g suchego materiału, gęstość ρ = 1,4 g/cm³, cylinder o średnicy 6 cm (S = 28,3 cm²). x = 20/28,3 = 0,707 g/cm². Skład: krzemionka z pyłem węglowym. µ przy E = 46,5 keV (Pb-210): szacujemy µ ≈ 0,50 cm²/g (przybliżenie dla lekkiej matrycy). µx = 0,50 × 0,707 = 0,354.

f = (1 − e^(−0,354)) / 0,354 = (1 − 0,702) / 0,354 = 0,298 / 0,354 = 0,842

Tylko 84,2% fotonów Pb-210 opuszcza próbkę. Pominięcie korekty zawyżyłoby obliczone tło i zaniżyło aktywność Pb-210 o 15,8%. W badaniach datowania przez Pb-210 taki błąd przekłada się bezpośrednio na błędną ocenę tempa akumulacji osadów. Przy x = 2 g/cm² (próbka grubsza):
µx = 1,00, f = (1 − e^(−1,00)) / 1,00 = 0,632 — korekta 36,8%.

Przykład 3 — efekt wilgotności na samopochłanianie przy 662 keV

Gleba sucha: ρ = 1,3 g/cm³, µ = 0,0862 cm²/g. Naczynie cylindryczne, masa 200 g, S = 70 cm², x = 2,86 g/cm². µx = 0,0862 × 2,86 = 0,246. f_sucha = (1 − e^(−0,246)) / 0,246 = (1 − 0,782) / 0,246 = 0,218 / 0,246 = 0,886.

Gleba wilgotna (30% wody): masa wzrasta do 260 g (200 g gleby + 60 g wody). Efektywne µ_eff = (0,0862 × 200 + 0,0857 × 60) / 260 = (17,24 + 5,14) / 260 = 0,0861 cm²/g. x_eff = 260/70 = 3,71 g/cm². µx_eff = 0,0861 × 3,71 = 0,319. f_mokra = (1 − e^(−0,319)) / 0,319 = (1 − 0,727) / 0,319 = 0,273 / 0,319 = 0,856.

Różnica f: f_sucha = 0,886, f_mokra = 0,856. Różnica korekty: 3,4% — nie duża przy 662 keV. Jeśli aktywność w Bq/kg podawana jest na suchą masę, ale pik mierzony ze świeżą (mokrą) próbką i kalibracja z suchą — różnica wynosi ~3% przez samą zmianę matrycy plus musi być uwzględniona normalizacja masy. Przy niskiej energii (60 keV) podobne obliczenie dałoby >10% efekt wilgotności.

Pytania otwarte

  1. Dlaczego wzór f = (1 − e^(−µx)) / (µx) nie jest adekwatny dla próbki w naczyniu Marinelli, a jedynie dla cienkiej jednorodnej warstwy? Jakie przybliżenia zostały poczynione i kiedy są zbyt duże?

  2. W tabeli współczynników osłabienia różnica µ między wodą a NaCl jest duża przy 60 keV (czynnik ~2,3), ale mała przy 662 keV (czynnik ~1,1). Jak to wpływa na dobór matrycy wzorcowej dla laboratoriów oznaczających Cs-137 i Pb-210 w próbkach wody morskiej (matryca NaCl)?

  3. Metoda transmisji zakłada jednorodność próbki. Jak wpłynęłaby na wynik niejednorodna próbka gleby z gęstymi skupieniami minerałów? Czy metoda transmisji dałaby wynik zawyżony czy zaniżony µx, i dlaczego?

  4. LabSOCS oblicza krzywe wydajności bez wzorców fizycznych. Jakie dane wejściowe są konieczne i jakie są wymagania akredytacyjne dotyczące walidacji tej metody? Co różni walidację LabSOCS od zwykłej kalibracji matrycowej wzorcami CRM?

  5. Dla próbki żywności mierzonej metodą Marinelli: masz do dyspozycji wzorzec Cs-137 w naczyniu Marinelli wypełnionym wodą. Próbka to suszony herbata (gęstość 0,2 g/cm³). Oszacuj jakościowo, czy kalibracja wodna zawyży czy zaniży aktywność Cs-137 w herbacie, i dlaczego.

  6. Pb-210 jest używany do datowania osadów jeziora metodą CRS (constant rate of supply). Błąd 20% w korekcie samopochłaniania jak propaguje do błędu datowania i tempa akumulacji osadu? Jakie konsekwencje ma to dla interpretacji historii środowiskowej jeziora?

  7. Rozporządzenie EU 2016/52 definiuje limity aktywności Cs-137 w żywności po Fukushimie. Czy laboratorium mierząc Cs-137 metodą Marinelli bez korekcji samopochłaniania może fałszywie stwierdzić zgodność próbki z limitem? Przy jakich parametrach próbki (gęstość, zawartość soli) błąd byłby największy?

  8. IAEA zaleca stosowanie ISO 20042 i TECDOC dokumentów dla spektrometrii gamma. Które wymagania dotyczące korekty samopochłaniania musi spełnić laboratorium, żeby uzyskać akredytację PCA dla oznaczania Cs-137 w glebie metodą Marinelli?

Podsumowanie dydaktyczne

  1. Samopochłanianie jest stratą sygnału z fotopiku, która powstaje, gdy fotony gamma są pochłaniane lub rozpraszane wewnątrz samej próbki na drodze do detektora. Efekt jest tym większy, im niższa energia gamma, im większa masa próbki i im gęstsza lub ciężka chemicznie matryca.

  2. Uproszczona korekcja f = (1 − e^(−µx)) / (µx) stosuje się dla cienkiej jednorodnej warstwy. Dla geometrii Marinelli, cylindra lub 3D konieczne są modele Monte Carlo (LabSOCS, MCNP) lub metoda transmisji — wzór warstwowy jest zbyt przybliżony.

  3. Masowy współczynnik osłabienia µ (cm²/g) jest silnie zależny od energii i składu chemicznego. Przy 662 keV różnice między typowymi matrycami środowiskowymi są małe (<15%); przy 60 keV mogą sięgać czynnik 4 (woda vs. ruda żelaza). Dlatego kalibracja matrycowa jest kluczowa przy niskich energiach.

  4. Najlepsza metoda eliminacji błędu samopochłaniania: kalibracja wzorcem w identycznej geometrii i identycznym składzie matrycy co próbka. Wtedy wydajność w fotopiku ε wyznaczona ze wzorca „zawiera" już samopochłanianie, bez potrzeby osobnej korekty f.

  5. Wilgotność próbki wpływa na samopochłanianie przez zmianę gęstości i składu efektywnego. Przy 662 keV efekt jest mały (<5%), przy 60 keV może wynosić >10%. Standard: suszyć próbki i podawać aktywność na suchą masę.

  6. Symulacje Monte Carlo (LabSOCS, MCNP, GEANT4) są metodą referencyjną dla złożonych geometrii. Wymagają dokładnego składu chemicznego i geometrycznego opisu próbki, ale pozwalają kalibrować bez fizycznych wzorców dla każdej geometrii.

  7. Normy ISO 20042, IAEA TECDOC-619 i wymagania akredytacyjne PCA/ISO 17025 nakładają obowiązek dokumentacji i walidacji korekty samopochłaniania. Wynik bez tej informacji jest niekompletny metrologicznie.

  8. Kluczowa reguła: aktywność gamma z próbki objętościowej bez opisu kalibracji matrycowej, gęstości, energii, geometrii i korekty samopochłaniania jest semi-ilościową oceną, nie pomiarem. Różnica między bezdekumentowanym wynikiem a w pełni opisanym pomiarem może wynosić 10–50% systematycznie.

Dodatkowe materiały multimedialne

Warto przygotować kalkulator samopochłaniania: użytkownik podaje energię, masowy współczynnik osłabienia, masę powierzchniową i wydajność bez samopochłaniania. Narzędzie oblicza f = (1 - exp(-µx))/(µx), pokazuje zmianę aktywności i ostrzega, gdy uproszczony model warstwowy jest zbyt słaby dla geometrii objętościowej.

Gotowe narzędzia: kalkulator samopochłaniania gamma, wizualizacja samopochłaniania próbki i model 3D pojemnika Marinelli.

Najkrótsze podsumowanie: samopochłanianie jest stratą sygnału wewnątrz samej próbki. Nie znika dlatego, że spektrometr jest dobrze skalibrowany, chyba że kalibracja została wykonana w tej samej geometrii i podobnej matrycy. Przy próbkach objętościowych aktywność z fotopiku wymaga myślenia o drodze fotonu przez materiał.

Ćwiczenia praktyczne

Pierwsze ćwiczenie rachunkowe: dla danych syntetycznych obliczyć f = (1 - exp(-µx))/(µx) dla kilku wartości µx: 0,01, 0,1, 0,5, 1,0 i 2,0. Student ma opisać, kiedy poprawka jest mała, a kiedy zaczyna dominować wynik.

Drugie ćwiczenie: porównać dwie próbki o tej samej aktywności właściwej, ale różnej masie powierzchniowej. Obliczyć, jak zmienia się pole fotopiku, jeśli epsilon_0 i I_gamma są takie same, a różni się tylko f.

Trzecie ćwiczenie: użyć syntetycznych danych dla K-40 w próbce soli. Student dostaje masę, masę molową, udział K-40, okres półtrwania, intensywność linii gamma i pole piku netto. Ma osobno policzyć aktywność teoretyczną z masy oraz aktywność widzianą przez spektrometr po uwzględnieniu f.

Czwarte ćwiczenie: przygotować tabelę decyzji, czy osobna poprawka samopochłaniania jest potrzebna. Warianty obejmują: wzorzec wodny i próbkę wodną, wzorzec wodny i suchą glebę, wzorzec matrycowy i próbkę o tej samej gęstości, niskoenergetyczną linię w ciężkiej matrycy oraz linię 1461 keV w lekkiej matrycy.

Piąte ćwiczenie audytowe: przeczytać opis pomiaru aktywności w Bq/kg i zaznaczyć, czy podano geometrię, masę, gęstość, energię, wzorzec, matrycę i sposób uwzględnienia samopochłaniania. Jeśli nie, wynik oznaczyć jako niepełny ilościowo.

Przejdź do ćwiczenia interaktywnego

Powiązane artykuły