Spowalnianie neutronów i mierniki wodoru: od moderatora do pomiaru wilgotności

Streszczenie

Neutron szybki w materiale bogatym w wodór zachowuje się trochę jak kula bilardowa uderzająca w drugą kulę o prawie tej samej masie: może oddać bardzo dużą część energii w pojedynczym zderzeniu. Dlatego woda, parafina, polietylen i inne materiały wodorowe są tak skutecznymi moderatorami neutronów.¹

Ta sama fizyka pozwala budować neutronowe mierniki wodoru i wilgotności. Przyrząd nie „widzi wody” bezpośrednio; widzi zmianę widma i liczby neutronów termicznych lub epitermicznych po oddziaływaniu z materiałem, a dopiero kalibracja zamienia sygnał na gęstość wodoru albo zawartość wilgoci.¹

Rozszerzenie tematu

Od zderzenia do moderatora

Neutron nie traci energii przez jonizację tak jak cząstka naładowana. Jego główną drogą spowalniania w typowych moderatorach są zderzenia sprężyste z jądrami atomowymi. Jeżeli jądro jest ciężkie, neutron odbija się od niego jak lekka kulka od ciężkiej ściany i traci małą część energii. Jeżeli jądro ma masę podobną do neutronu, energia może zostać przekazana bardzo skutecznie.

Dziunikowski i Kalita podają dla zderzenia centralnego prosty wzór na stosunek energii neutronu po zderzeniu do energii przed zderzeniem:

E / E0 = ((M - m) / (M + m))^2,

gdzie m jest masą neutronu, a M masą jądra moderatora.¹ Dla węgla, którego masa jądra jest około dwanaście razy większa od masy neutronu, w zderzeniu czołowym neutron może stracić około 28% energii. Dla wodoru, którego jądro jest protonem o masie prawie równej masie neutronu, strata w zderzeniu czołowym może być bliska 100%.¹

To wyjaśnia, dlaczego wodór jest tak ważny w praktyce neutronowej. Nie chodzi tylko o to, że woda jest tania albo łatwo dostępna. Chodzi o dopasowanie mas: neutron i proton są prawie równymi partnerami w zderzeniu. Właśnie dlatego materiały zawierające dużo wodoru silnie zmieniają widmo neutronów.

Statystyczna utrata energii

Zderzenia w materiale nie są zawsze czołowe. Neutrony rozpraszają się pod różnymi kątami, więc utrata energii w pojedynczym zderzeniu jest zmienną losową. Dlatego wygodnie używać średniego logarytmicznego ubytku energii xi:

xi = < ln(E0 / E) >.

Jeżeli neutron ma spaść z energii E1 do E2, przybliżona średnia liczba zderzeń wynosi:

nu = ln(E1 / E2) / xi.

Dziunikowski i Kalita podają przykładową tabelę dla spowalniania neutronu od 2 MeV do energii termicznej 0,025 eV: dla wodoru potrzeba średnio około 18 zderzeń, dla deuteru około 25, dla węgla około 114, a dla uranu około 2150.¹ Tabela dobrze pokazuje różnicę między „materiałem, który może rozpraszać neutrony” a dobrym moderatorem.

Sama wartość xi nie wystarcza. Liczy się też liczba jąder w jednostce objętości i przekrój czynny na rozpraszanie. Dlatego wprowadza się zdolność spowalniania:

Z_s = xi * n0 * sigma_s,

gdzie n0 oznacza liczbę jąder w jednostce objętości, a sigma_s przekrój czynny na rozpraszanie.¹ Dobry moderator musi często zderzać neutrony i w każdym zderzeniu zabierać im znaczną część energii.

Moderator dobry nie tylko przez spowalnianie

W reaktorze albo detektorze nie wystarczy neutrony spowolnić. Trzeba jeszcze ich nie pochłonąć w niepożądanym miejscu. Dlatego obok zdolności spowalniania istotny jest stosunek spowalniania do absorpcji. Dziunikowski i Kalita definiują współczynnik spowalniania jako:

W_s = xi * Sigma_s / Sigma_a,

czyli w uproszczeniu stosunek skutecznego rozpraszania do pochłaniania.¹

To tłumaczy klasyczny kontrast między zwykłą wodą, ciężką wodą i grafitem. Zwykła woda bardzo skutecznie spowalnia neutrony, ale wodór ma niezerową absorpcję. Ciężka woda spowalnia gorzej na jedno zderzenie, bo deuter jest cięższy od protonu, ale ma bardzo mały przekrój absorpcji, więc w reaktorach może być znakomitym moderatorem. Grafit spowalnia wolniej, lecz przy wysokiej czystości może mieć małe straty absorpcyjne. Ten temat łączy się bezpośrednio z artykułami o ciężkiej wodzie, graficie reaktorowym i historycznym błędzie niemieckiego programu jądrowego.²

W miernikach wodoru interesuje nas jednak trochę inny aspekt. Nie szukamy najlepszego moderatora reaktorowego. Chcemy, aby sygnał przyrządu zmieniał się wyraźnie wraz z gęstością wodoru w badanym materiale.

Dlaczego wodór jest dobrym znacznikiem wilgotności

Woda zawiera dużo wodoru, a wiele suchych materiałów mineralnych zawiera go mało albo w znanej ilości. Jeżeli materiał suchy praktycznie nie ma wodoru, a dodana wilgoć jest głównym źródłem jąder wodoru, gęstość wodoru staje się dobrym wskaźnikiem wilgotności. To jest główna zasada neutronowych mierników wilgotności.

Dziunikowski i Kalita opisują warunek bardzo jasno: pomiar wilgotności metodą neutronową jest sensowny wtedy, gdy związek między całkowitą zawartością wodoru i zawartością wody w materiale jest jednoznaczny.¹ Jeżeli sucha masa zawiera wodór w nieznanej ilości, miernik neutronowy nadal mierzy wodór, ale przeliczenie na wodę staje się niepewne albo wręcz błędne.

To prowadzi do pierwszej ważnej lekcji metrologicznej: neutronowy miernik wilgotności nie jest uniwersalnym „miernikiem wody”. Jest miernikiem efektu neutronowego zależnego od wodoru. W piasku, kruszywie, rudzie albo masie formierskiej może działać dobrze, jeżeli matryca jest znana. W materiale organicznym, tworzywie sztucznym albo mieszaninie z wodorem związanym chemicznie może zawyżać wilgotność, jeżeli kalibracja nie uwzględnia składu suchej masy.

Zasada działania miernika wodoru

W modelu pojęciowym miernik wodoru składa się z pola neutronowego, próbki, moderatora, detektora neutronów termicznych i toru zliczającego. Neutrony wchodzą w materiał, zderzają się z jądrami, a wodór szczególnie silnie je spowalnia. Im więcej wodoru w objętości pomiarowej, tym bardziej zmienia się liczba neutronów, które wracają do detektora w zakresie energii, na który detektor jest czuły.¹

Dziunikowski i Kalita opisują, że w liczniku BF3 reakcja (n, alpha) ma przekrój czynny zależny od energii zgodnie z regułą 1/v, więc neutrony spowolnione do małych energii są łatwiej rejestrowane.¹ W uproszczeniu próbka bogatsza w wodór mocniej termalizuje neutrony, a to zmienia częstość zliczeń. Konkretna zależność zależy od konstrukcji przyrządu, widma wejściowego, geometrii, materiałów osłonowych i detektora, dlatego nie jest to czysta stała fizyczna. Jest to krzywa kalibracyjna.

W bezpiecznym artykule dydaktycznym nie trzeba opisywać operacyjnej konstrukcji miernika ani pracy z realnym źródłem neutronów. Wystarczy zrozumieć tor przyczynowy:

• wodór silnie spowalnia neutrony,
• detektor termiczny reaguje inaczej na neutrony spowolnione niż szybkie,
• częstość zliczeń koreluje z gęstością wodoru,
• krzywa kalibracyjna zamienia zliczenia na wilgotność tylko dla określonej klasy materiałów.

Krzywa kalibracyjna

Najprostsza kalibracja może mieć postać liniową:

J = a0 + a1 * rho_H,

gdzie J jest częstością zliczeń, rho_H gęstością wodoru, a a0 i a1 współczynnikami wyznaczonymi na wzorcach.¹ W praktyce zakres liniowości może być ograniczony. Przy dużej wilgotności, dużej absorpcji albo skomplikowanej geometrii odpowiedź może stać się nieliniowa.

Kalibracja musi odpowiadać materiałowi. Wzorce dla kruszywa nie muszą być dobre dla rudy metalu, a wzorce dla piasku nie muszą działać dla materiału z domieszką związków organicznych. Jeżeli sucha matryca zawiera wodór, trzeba odjąć albo uwzględnić wkład wodoru niezwiązanego z wodą.

Rachunkowo przeliczenie może wyglądać tak:

rho_H = (J - a0) / a1.

Jeżeli znamy gęstość nasypową materiału i udział wodoru w wodzie, można przejść od gęstości wodoru do wilgotności masowej. Dla wody H2O masa wodoru stanowi około 2/18, czyli około 11,1% masy wody. Zatem masa wody odpowiadająca masie wodoru jest około 9 razy większa. To proste przeliczenie jest poprawne tylko wtedy, gdy wodór mierzony przez przyrząd rzeczywiście pochodzi z wody.

Przykład modelowy bez źródła neutronów

Załóżmy syntetyczną kalibrację:

J = 120 + 85 * w,

gdzie w oznacza wilgotność w procentach masowych w określonej klasie suchego materiału. Jeżeli zmierzona częstość zliczeń netto wynosi 545 cps, model daje:

w = (545 - 120) / 85 = 5,0%.

To nie jest uniwersalny wzór na wilgotność. To wyłącznie przykład obliczeniowy pokazujący, jak działa kalibracja. W realnej metrologii trzeba znać materiał, objętość pomiarową, stabilność aparatu, tło, czas zliczania, niepewność wzorców i wpływ gęstości nasypowej.

W ćwiczeniu syntetycznym można dodać szum Poissona. Jeżeli pomiar trwa t = 60 s, a oczekiwana częstość wynosi 545 cps, oczekiwana liczba zliczeń to 32700, a statystyczne odchylenie standardowe to około sqrt(32700) = 181. Niepewność częstości wynosi więc około 3,0 cps, co przekłada się na niepewność wilgotności około 3,0 / 85 = 0,035% w samym składniku statystycznym. W praktyce błędy kalibracji i matrycy mogą być większe.

Efekty matrycy

Najpoważniejsze błędy nie muszą pochodzić ze zliczeń. W miernikach wodoru szczególnie ważne są efekty matrycy:

• sucha masa może zawierać wodór związany chemicznie,
• materiał może zawierać silne pochłaniacze neutronów,
• gęstość nasypowa może zmieniać liczbę jąder w objętości pomiarowej,
• skład ziarnowy może wpływać na porowatość i rozmieszczenie wody,
• wilgoć może być nierównomiernie rozłożona,
• geometria próbki może nie odpowiadać geometrii wzorcowania,
• temperatura i skład chemiczny mogą zmieniać odpowiedź detektora lub materiałów pomocniczych.

Bor, kadm, gadolin i niektóre inne nuklidy mogą bardzo silnie pochłaniać neutrony termiczne. Jeżeli występują w materiale albo w pobliżu objętości pomiarowej, mogą zmienić sygnał w sposób, którego prosta kalibracja „zliczenia = wilgotność” nie obejmuje. To jest analogiczne do efektów matrycy w rentgenowskiej analizie fluorescencyjnej: przyrząd mierzy odpowiedź fizyczną całego układu, nie samą interesującą wielkość chemiczną.

Głębokość i objętość pomiarowa

Miernik neutronowy nie mierzy punktu matematycznego. Mierzy pewną objętość materiału, której rozmiar zależy od energii neutronów, przekrojów czynnych, geometrii i składu próbki. Gdy wilgotność jest duża, neutrony szybciej tracą energię i przestrzeń skutecznego pomiaru może się zmieniać. Gdy materiał jest suchy albo mało wodorowy, neutrony mogą penetrować dalej.

Dla interpretacji wyniku oznacza to, że lokalna niejednorodność ma znaczenie. Jeżeli w objętości pomiarowej są warstwy o różnej wilgotności, pustki, duże ziarna albo obce wtrącenia, wynik jest średnią ważoną odpowiedzią neutronową, a nie prostą średnią arytmetyczną wilgotności.

W dobrym raporcie trzeba więc zapisać nie tylko wartość wilgotności, ale także typ materiału, metodę kalibracji, czas zliczania, liczbę powtórzeń, gęstość nasypową i informację, czy próbka była jednorodna w skali objętości pomiarowej.

Związek z reaktorami i bronią jądrową

Ten sam mechanizm spowalniania neutronów przez lekkie jądra pojawia się w kilku działach serwisu. W reaktorach moderator decyduje, czy neutrony rozszczepieniowe zostaną spowolnione do zakresu, w którym przekrój czynny rozszczepienia U-235 jest duży. W historii programu niemieckiego wybór między grafitem a ciężką wodą był jednym z kluczowych ograniczeń technologicznych.²

W broni jądrowej moderator nie pełni tej samej roli co w reaktorze energetycznym, ale materiały bogate w wodór nadal wpływają na transport neutronów, osłony, diagnostykę i środowisko testowe. Dla dydaktyki najważniejsze jest rozumienie ogólnej fizyki: neutrony szybkie, pośrednie i termiczne nie są „tym samym strumieniem” z inną nazwą. Materiał może je rozpraszać, pochłaniać, termalizować albo odbijać, a skutki zależą od energii.

Czego nie mówi sam wynik zliczeń

Jeżeli miernik pokazuje większą częstość zliczeń, łatwo powiedzieć: „jest więcej wody”. To może być prawda tylko w dobrze skalibrowanym i jednorodnym przypadku. Sam wynik zliczeń nie mówi:

• czy wodór pochodzi z wody, czy z materiału organicznego,
• czy matryca nie zawiera pochłaniaczy neutronów,
• czy gęstość próbki jest zgodna ze wzorcami,
• czy pole neutronowe i detektor są stabilne,
• czy odpowiedź jest nadal liniowa,
• czy mierzona objętość odpowiada rzeczywistemu pytaniu technologicznemu.

Właśnie dlatego mierniki neutronowe są narzędziami metrologicznymi, a nie magicznymi czujnikami wilgoci. Ich siła polega na tym, że mogą badać materiał nieniszcząco i objętościowo. Ich słabość polega na tym, że wynik jest wiarygodny tylko w ramach poprawnej kalibracji.

Minimalny model kalkulatora

Przyszły kalkulator dla tego artykułu powinien mieć dwa tryby. Pierwszy tryb powinien liczyć spowalnianie: użytkownik wybiera energię początkową, energię końcową i materiał moderatora, a narzędzie oblicza przybliżoną liczbę zderzeń nu = ln(E1/E2)/xi. Dla wodoru, deuteru, węgla i tlenu można pokazać, jak zmienia się liczba zderzeń.

Drugi tryb powinien liczyć kalibrację miernika wodoru. Użytkownik wpisuje punkty wzorcowe rho_H, J, program dopasowuje J = a0 + a1 rho_H, pokazuje reszty, niepewność nachylenia i przelicza nowy pomiar na gęstość wodoru. W wariancie rozszerzonym można dodać przeliczenie na wilgotność przy założeniu, że cały dodatkowy wodór pochodzi z wody.

Najkrótsze podsumowanie: wodór jest wyjątkowy w metrologii neutronowej, bo ma masę prawie taką jak neutron. Miernik wilgotności wykorzystuje ten fakt, ale wynik jest tak dobry, jak kalibracja, znajomość matrycy i kontrola geometrii pomiaru.

Źródła neutronów stosowane w miernikach wilgotności — aspekty regulacyjne

Każdy neutronowy miernik wilgotności zawiera radioizotopowe źródło neutronów — najczęściej Am-241/Be (amerykowo-berylowe). Regulacje prawne dotyczące posiadania i użytkowania tych źródeł są istotnym elementem praktycznym, o którym podręczniki rzadko wspominają.

Am-241/Be w Polsce: Am-241 jest zaliczany do kategorii materiałów jądrowych objętych nadzorem Państwowej Agencji Atomistyki (PAA). Użytkownik źródła musi posiadać zezwolenie PAA, źródło musi być ujęte w rejestrze i podlega corocznej inwentaryzacji. Przenośne mierniki z Am-241/Be używane na budowach lub w rolnictwie wymagają więc dodatkowych procedur transportowych zgodnych z ADR (umowa europejska o transporcie drogowym).

Alternatywy nieradioizotopowe: Generatory neutronów DD (deuteron-deuteron, E_n ≈ 2,5 MeV) i DT (deuteron-tryton, E_n = 14,1 MeV) działają wyłącznie przy zasilaniu — po wyłączeniu nie emitują promieniowania. Nie wymagają zezwolenia PAA (brak materiału jądrowego), ale wymagają rejestracji urządzenia wytwarzającego promieniowanie jonizujące (dyrektywa BSS). Są droższe w zakupie, lecz prostsze regulacyjnie. W zastosowaniach on-line na wytwórniach betonu coraz częściej zastępują Am-Be.

Bezpieczeństwo źródeł mobilnych: Źródło Am-241/Be niesione przez operatora na plecach lub w ręku w terenie stwarza ryzyko dawki promieniowania neutronowego i gammma. Operatorzy powinni korzystać z pleksiglasowych lub polietylenowych osłon i ograniczać czas ekspozycji. Szkolenie z ochrony radiologicznej jest obowiązkowe.

Moderatory reaktorowe — porównanie i kryteria wyboru

W reaktorach jądrowych rola moderatora jest podwójna: spowalniać neutrony szybkie (wyemitowane przez rozszczepienie z energią ~2 MeV) do energii termicznej (~0,025 eV), gdzie przekrój czynny rozszczepienia U-235 jest duży (σ ≈ 584 barn), jednocześnie absorpcja w moderatorze musi być na tyle mała, żeby reaktor mógł działać krytycznie.

Parametry kluczowe moderatora:

• Zdolność spowalniania: ξ × Σ_s (duże = dobre)
• Współczynnik spowalniania (moderating ratio): ξ × Σ_s / Σ_a (duże = mało traci neutronów)

Zestawienie moderatorów reaktorowych:

Moderator	ξ	Σ_s (1/cm)	Σ_a (1/cm)	Wsp. spowalniania	Uwagi
H₂O (zwykła woda)	0,920	1,475	0,022	~62	Tania, duże Σ_a; wymaga paliwa wzbogaconego
D₂O (ciężka woda)	0,509	0,353	0,000033	~5500	Bardzo niska absorpcja; kosztowna; CANDU, ATUCHA
C (grafit)	0,158	0,385	0,000388	~157	Tani, duży, wymaga high purity; Chernobyl RBMK
Be (beryl)	0,209	0,700	0,0012	~122	Toksyczny, drogi; reaktory badawcze
BeO	0,174	0,476	0,00084	~98	Jak Be, ceramika

H₂O dominuje w LWR (PWR i BWR) — reaktory pracują na paliwie 3–5% UOX lub MOX. Gdyby użyć paliwa naturalnego (0,72% U-235) z H₂O, reaktor byłby pod-krytyczny z powodu zbyt dużej absorpcji w wodorze.

D₂O jest kluczem do CANDU — naturalny uran + D₂O działa krytycznie. To daje strategiczną niezależność: Kanada (i Indie, Argentyna, Korea, Pakistan, Rumunia) nie muszą wzbogacać uranu.

Grafit RBMK (ZSRR/Rosja): największe reaktory (1000–1500 MWe), paliwo lekko wzbogacone (2,0–2,8% U-235), moderator grafitowy, chłodziwo H₂O. Niekorzystny współczynnik temperaturowy czynnika chłodzącego (wzrost temperatury → parowanie → zmniejszenie absorpcji → wzrost mocy) był jedną z przyczyn katastrofy w Czarnobylu.

Dyfuzja neutronów i zasięg migracji

Neutron termiczny w moderatorze nie stoi w miejscu — dyfunduje losowo, odbijając się od jąder. Opisuje to równanie dyfuzji Ficka:
∂φ/∂t = D∇²φ − Σ_a×φ + S,
gdzie φ — fluencja, D = 1/(3Σ_tr) — współczynnik dyfuzji, Σ_tr — makroskopowy przekrój czynny transportu, S — źródło.

Długość dyfuzji L (odpowiednik średniej drogi od termalizacji do absorpcji):
L² = D / Σ_a = 1 / (3 × Σ_tr × Σ_a).

Dla H₂O: L ≈ 2,73 cm. Neutron termiczny w wodzie dyfunduje przeciętnie 2,73 cm od miejsca termalizacji do miejsca absorpcji.
Dla D₂O: L ≈ 170 cm — dramatycznie większy zasięg dyfuzji, co wyjaśnia, czemu reaktory D₂O mogą mieć paliwo rozmieszczone rzadziej.
Dla grafitu: L ≈ 54 cm.

Praktyczne znaczenie dla mierników: Objętość pomiarowa neutronowej sondy gleby to nie kula o promieniu L, ale sfera efektywnego moderowania i dyfuzji. Im suchy materiał (mało H), tym dłuższa droga neutronów → większa objętość pomiarowa. Im wilgotniejszy → krótsza → mniejsza sfera.

Spowolnienie neutronów szybkich — prawo 1/v i rezonanse

Reguła 1/v (Breit-Wigner) mówi: dla neutronów termicznych σ ∝ 1/v ∝ E^(−½). Jest to dobry opis dla większości substancji pochłaniających w zakresie 0,001–1 eV.

Jednak wiele nuklidów ma rezonanse absorpcji w zakresie epitermicznym (1 eV – 100 keV). Są to wąskie piki σ(E) opisane wzorem Breit-Wignera:
σ(E) = σ₀ × (E_r/E)^½ × Γ²/[(E−E_r)² + Γ²/4],
gdzie E_r — energia rezonansowa, Γ — szerokość rezonansowa.

Rezonanse kluczowe dla praktyki:

• Au-197: rezonans przy 4,9 eV, σ_peak ≈ 30 000 barn — klucz metody foliowej;
• In-115: rezonans przy 1,46 eV, σ_peak ≈ 30 000 barn;
• U-238: silne rezonanse w zakresie 6–200 eV (np. 6,67 eV, 20,9 eV) — pochłania neutrony epitermiczne w paliwie reaktorowym (efekt Dopplera stabilizuje reaktor przy wzroście temperatury).

Efekt Dopplera: Przy wzroście temperatury paliwa jądra U-238 drgają szybciej → rezonanse „rozszerzają się" energetycznie → więcej neutronów epitermicznych pochłoniętych w U-238 → liczba neutronów do rozszczepienia maleje → moc spada. To ujemny współczynnik reaktywności temperaturowej — ważny mechanizm bezpieczeństwa PWR/BWR.

Praktyczny wniosek dla pomiarów: W miernikach neutronowych rezonanse absorpcji materiałów w gruncie (np. In, Au, Fe w rudach) mogą selektywnie pochłaniać neutronowe o pewnych energiach i deformować widmo powracające do detektora. To kolejny powód, dla którego kalibracja na materiale referencyjnym odpowiadającym rzeczywistej próbce jest niepomijalnym elementem dobrej praktyki pomiarowej. Ignorowanie składu mineralogicznego gruntu może prowadzić do błędu wilgotności rzędu kilku procent — szczególnie w terenach o bogatej mineralizacji żelaza, boru lub litu.

Typy detektorów i geometrie mierników wodoru

Podstawą większości przenośnych mierników wodoru jest licznik proporcjonalny BF₃: cylindryczna rura z trifluorkiem boru pod ciśnieniem 0,5–1 atm wzbogaconym w izotop B-10 (naturalny bor: 19,9% B-10, wzbogacony: 96–99% B-10). Reakcja ¹⁰B(n,α)⁷Li ma przekrój czynny σ ≈ 3835 barn dla neutronów termicznych i silnie rośnie zgodnie z regułą 1/v przy obniżaniu energii.

Licznik He-3: Alternatywa dla BF₃, wyższy przekrój czynny (5327 barn), mniej toksyczny (gaz niejadowity), lecz od ok. 2008 r. He-3 jest trudno dostępny z powodu ograniczonej produkcji (pochodzi z rozpadu trytu z głowic jądrowych). Stosowany w wysokiej klasy instrumentach badawczych i safeguardowych.

Scyntylatory LiI(Eu) i Li-glass: Kryształ LiI(Eu) lub szkło litowe są wrażliwe na neutrony termiczne przez reakcję ⁶Li(n,α)³H (σ = 940 barn). Produkują błysk świetlny mierzony przez PMT lub SiPM. Niższy przekrój niż BF₃/He-3, ale dobre parametry energetyczne i brak toksycznych gazów.

Detektory PIPS z folią ⁶LiF: Plastikowa lub metalowa matryca z ⁶LiF przed detektorem PIPS. Cząstki α i ³H z reakcji wchodzą do detektora i są zliczane. Metoda stosowana w specjalistycznych przenośnych sondach.

Geometria sondy gleby: Klasyczna sonda do pomiaru wilgotności gleby (soil moisture probe) ma cylindryczne obudowanie zawierające źródło neutronów (Am-241/Be) i detektor. Spuszczana jest w pionowy otwór w gruncie i zlicza neutrony powracające z otaczającego gruntu. Objętość pomiarowa to sfera o promieniu 10–30 cm w zależności od wilgotności (więcej wody → krótsza droga → mniejsza sfera).

Geometria miernika przenośnego (powierzchniowego): Przyrząd kładziony na badanej powierzchni. Jedno gniazdo zawiera źródło, drugie detektor. Neutrony wchodzą w materiał, spowalniają się i wracają. Stosowany na budowach (kontrola wilgotności betonów i gruntów), w kopalniach (wilgotność rudy) i w rolnictwie (monitoring pola).

Mierniki przepływowe (on-line): W przemyśle budowlanym i spożywczym — głowice zamontowane stacjonarnie ponad taśmociągiem. Wynik ciągły w czasie rzeczywistym. Stosowane do kontroli wilgotności kruszywa, ziaren, makaronu, papieru.

Szczegółowe omówienie efektów matrycy z tabelą

Efekty matrycy w miernikach neutronowych są analogiczne do efektów matrycy w XRF (gdzie absorpcja i wzmocnienie fluorescencji są funkcją składu). Tabela głównych efektów i sposoby ich kontroli:

Efekt	Przyczyna	Wpływ na zliczenia	Korekcja
Wodór w suchej masie	CH, OH w związkach organicznych	Zawyżenie wilgotności	Kalibracja na właściwy materiał, analiza LOI
Pochłanianie przez Cd, Gd, B	Duże σ dla neutronów termicznych	Zaniżenie zliczeń	Dodatkowy pomiar transmisji, uwzględnienie w kalibracji
Gęstość nasypowa	Zmiana liczby jąder w objętości pomiarowej	Zmiana stałej kalibracyjnej	Pomiar gęstości, kalibracja wielokryterialna
Zróżnicowanie ziarnowe	Niejednorodna porowatość, rozkład wody	Niejednorodna odpowiedź	Próbki reprezentatywne, wiele pomiarów
Temperatura	Zmiana gęstości i ruchliwości wody	Zmiana kształtu widma neutronów	Korekcja temperaturowa kalibracji
Zawartość chloru	Cl absorbuje neutrony (σ ≈ 43 barn)	Zaniżenie zliczeń	Uwzględnienie Cl w obliczeniach

Bor-10 jest szczególnie problematyczny: σ = 3835 barn (termiczne), więc już kilka procent boru w minerale (tourmalin, boraks, sassolite) może dramatycznie zmienić odpowiedź. Podobnie gadolin (σ = 49 000 barn dla termicznych — najwyższy znany) — jest obecny w rzadkich ziemiach i w odpuszczalnej rudzie.

Calibration — modele kalibracyjne

Model liniowy (1-parametrowy): J = a₀ + a₁ × θ, gdzie θ — wilgotność masowa lub objętościowa. Wystarczający dla materiałów jednorodnych w wąskim zakresie θ.

Model potęgowy: J = A × θ^B. Lepiej opisuje krzywoliniowość przy dużych wilgotnościach, gdzie efekt moderatora saturuje.

Model wielokryterialny (MLR): J = a₀ + a₁×θ + a₂×ρ + a₃×[B], gdzie ρ — gęstość, [B] — stężenie boru lub innego pochłaniacza. Stosowany gdy matryca jest znana, ale zmienia się w czasie (kontrola procesów przemysłowych).

Kalibracja wzorcem TDR: Porównanie neutronowego miernika wilgotności z Time Domain Reflectometry (TDR) lub grawimetrycznym pomiarem wilgotności. Kilkanaście wzorców w różnych warunkach. Wzorzec grawimetryczny: próbkę gleby ważymy, suszymy w 105°C do stałej masy, ważymy ponownie. θ_grawim = (masa_wilg − masa_sucha) / masa_sucha.

Neutronowe sondy gleby w nauce i praktyce

Hydrologia i agronomia: Pomiar wilgotności gleby jest kluczowy do modelowania bilansu wodnego, planowania nawadniania i kalibracji modeli hydrologicznych. Neutronowe sondy gleby (np. CPN 503DR, Hydroprobe) były pierwszymi metodami nieinwazyjnego pomiaru wilgotności objętościowej w głębi profilu glebowego. Dziś konkurują z nimi TDR (Time Domain Reflectometry), FDR (Frequency Domain Reflectometry) i czujniki pojemnościowe — ale sondy neutronowe nadal mają przewagę w zakresie objętości pomiarowej i braku czułości na zasolenie.

Monitoring gruntów składowisk: Na składowiskach odpadów komunalnych i przemysłowych monitoruje się migrację odcieków w gruncie. Sonda neutronowa wykrywa zmianę zawartości wody w pakiecie uszczelniającym. Brak potrzeby wykopywania i instalacji długich kabli.

Przemysł cementowy: Kontrola wilgotności kruszywa (piasek, żwir) wchodzącego do wytwórni betonów gotowych. Zbyt mokre kruszywo zmniejsza wskaźnik w/c (woda/cement) bez wiedzy operatora, co zaburza wytrzymałość betonu. Mierniki montowane nad taśmociągiem pracują ciągłe online z pomiarami co kilka sekund.

Przemysł spożywczy i papierniczy: Wilgotność ziaren zboża, makaronu, pieczywa suchego, celulozy. Regulacje (normy ISO, EN) wymagają, by produkt miał wilgotność w określonym przedziale. Mierniki neutronowe są tu certyfikowane jako narzędzie produkcyjne.

Detekcja min i IED: Specjalna aplikacja wojskowa — miną lub IED zasypany gruntem zmienia lokalną gęstość wodoru (lub węgla). Systemy bazujące na rozpraszaniu neutronów (PFNA — pulsed fast neutron analysis) mogą wykrywać ładunki przez pokład pojazdu. To złożone zagadnienie — ważny kontekst historyczny i technologiczny, wykracza poza typowe „pomiar wilgotności".

Polska infrastruktura

W Polsce neutronowe pomiary środowiskowe i przemysłowe prowadzone są przez:

• Instytut Meteorologii i Gospodarki Wodnej (IMGW): stacje pomiarowe z sondami gleby, kalibracje do modelowania hydrologicznego;
• IUNG (Instytut Uprawy, Nawożenia i Gleboznawstwa, Puławy): badania wilgotności gleby w kontekście nawadniania i bilansu wodnego;
• AGH Kraków (Wydział Wiertnictwa, Nafty i Gazu): mierniki profilowe do odwiertów — gamma-gamma gęstościomierze i neutronowe mierniki porowatości stosowane w geofizyce wiertniczej;
• NCBJ Świerk: wzorcowania przyrządów pomiarowych, badania reaktorowe;
• Przemysł cementowy (Lafarge, CRH, Holcim — zakłady w Polsce): mierniki wilgotności kruszywa na wytwórniach betonu towarowego.

Geofizyczna metoda neutronowa stosowana w wiertnictwie naftowym (podciśnienie naftowe NMR i gamma/neutron log) jest pokrewna, lecz używa głowic w odwiercie. AGH i PIG (Państwowy Instytut Geologiczny) korzystają z tych metod do interpretacji profilów litologicznych.

Trzy przykłady rachunkowe

Przykład 1 — Liczba zderzeń do termalizacji i dominacja wodoru

Neutron startuje z energią E₀ = 2 MeV. Chcemy go zwolnić do energii termicznej E_term = 0,025 eV = 2,5×10⁻⁸ MeV.

Logarytmiczny ubytek energii: ln(E₀/E_term) = ln(2 / 2,5×10⁻⁸) = ln(8×10⁷) = 18,2.

Wartości ξ (średni logarytmiczny ubytek energii na zderzenie):

• Wodór (M=1): ξ_H = 1,000
• Deuter (M=2): ξ_D = 0,725
• Węgiel (M=12): ξ_C = 0,158
• Tlen (M=16): ξ_O = 0,120
• Uran (M=238): ξ_U ≈ 0,0084

Liczba zderzeń n = ln(E₀/E_term) / ξ:

• Wodór: n_H = 18,2 / 1,000 ≈ 18 zderzeń
• Węgiel: n_C = 18,2 / 0,158 ≈ 115 zderzeń
• Uran: n_U = 18,2 / 0,0084 ≈ 2167 zderzeń

Wniosek: Wodór jest 6× lepszym moderatorem niż węgiel pod względem liczby zderzeń. Uran praktycznie nie moderuje — neutrony szybkie w U-238 odbijają się bez utraty energii.

Przykład 2 — Przeliczenie zliczeń na wilgotność z uwzględnieniem tła

Miernik przenośny na gruncie. Kalibracja dla piasku kwarcowego (bez wodoru w suchej masie):
J = 180 + 75 × θ [cps, θ w % masowych]

Pomiar: czas t = 60 s, zliczone N_brutto = 21 780 impulsów. Tło bez materiału (powietrze): N_tło/60s = 300 impulsów.

Częstość netto: J_netto = (21 780 − 300) / 60 = 358,0 cps.

Wilgotność: θ = (J_netto − 180) / 75 = (358,0 − 180) / 75 = 2,37%.

Niepewność statystyczna: σ_N = √N_brutto = √21780 ≈ 147,6 → σ_J = 147,6/60 = 2,46 cps.
Niepewność θ = σ_J / 75 = 2,46 / 75 = 0,033%.

Wynik: θ = (2,37 ± 0,03)% — ale to tylko niepewność statystyczna. Jeśli kalibracja została wykonana na innym gatunku piasku (z 0,3% wodoru w suchej masie), błąd systematyczny wynosi 0,3% — 9-krotnie większy od statystycznego.

Przykład 3 — Wpływ pochłaniacza (bor) na wynik pomiaru

W gruncie odkrytym na terenie starych hut zawierającym 0,5% B₂O₃ masowo. Porównujemy czysty grunt (bez boru) i grunt z borem.

Bor-10 stanowi 20% naturalnego boru. W B₂O₃: % B-10 = 0,5% × (2×11 / 70) × 0,199 ≈ 0,5% × 0,314 × 0,199 = 0,031% B-10 masowo.

Makroskopowy przekrój czynny absorpcji dla neutronów termicznych: Σ_a = n × σ. Dla 0,031% B-10 w matrycy o gęstości 1,5 g/cm³: n_B10 = 0,031%/100 × 1,5 g/cm³ / (10 g/mol) × 6,022×10²³ = 2,8×10¹⁸ atomów/cm³. Σ_a = 2,8×10¹⁸ × 3835×10⁻²⁴ cm² = 10,7 cm⁻¹. To bardzo duże! Droga swobodna 1/Σ_a = 0,09 cm — neutrony termiczne prawie nie przejdą przez ten grunt.

Efekt: zliczenia drastycznie spadają (pozorny brak wilgoci), mimo że wilgotność rzeczywista jest np. 5%. Bez uwzględnienia boru kalibracja daje wynik bliski 0%. Konieczne: pomiar fluoresencji XRF na borze lub analiza chemiczna przed zastosowaniem miernika neutronowego.

Pytania otwarte

1. Dlaczego kalibracja neutronowego miernika wilgotności wykonana dla piasku nie jest ważna dla gliny zwałowej lub torfu? Jakie parametry materiału należy zmierzyć niezależnie, żeby przeliczyć wynik na rzeczywistą wilgotność?

2. Porównaj trzy metody pomiaru wilgotności in situ: neutronowa sonda gleby, TDR (Time Domain Reflectometry) i czujnik FDR. Dla każdej metody podaj jedną sytuację, w której ma zdecydowaną przewagę nad pozostałymi.

3. W neutronowym mierniku wilgotności kruszywa używanym na wytwórni betonów wynik jest zawyżony o 2% w stosunku do metody grawimetrycznej. Wymień trzy możliwe przyczyny rozbieżności i zaproponuj test, który pozwoli je odróżnić.

4. Reguła 1/v mówi, że przekrój czynny σ rośnie odwrotnie proporcjonalnie do prędkości neutronu. W temperaturze pokojowej neutrony termiczne mają v ≈ 2200 m/s. Jaki byłby przekrój czynny He-3 dla neutronów o energii 0,001 eV (ultra-zimne neutrony)?

5. Reaktor MARIA w Świerku używa zwykłej wody jako moderatora i pochłaniacza. W typowym reaktorze lekkowodnym woda pełni obie funkcje jednocześnie. Jakie jest kryterium równowagi między właściwościami moderatora a pochłaniacza, które decyduje, czy reaktor jest krytyczny?

6. Neutronowe mierniki porowatości stosowane w geofizyce wiertniczej mierzą „porowatość wodorową" (NPHI). Dlaczego w wapieniu, piaskowcu i dolomicie porowatość NPHI jest inna niż porowatość rzeczywista nawet przy tej samej ilości wody w porach?

7. Szybkie spowalnianie neutronów przez wodór stwarza problem radiologiczny: w pobliżu źródła neutronów prędkich w polietylenie generowane są neutrony termiczne z wysoką fluencją. Jak ten efekt wpływa na projektowanie osłon biologicznych reaktorów i akceleratorów?

8. W technice PFNA (pulsed fast neutron analysis) do wykrywania min i materiałów wybuchowych stosuje się neutrony prędkie (14 MeV z generatora DT). Jakie mierzalne wielkości fizyczne (gamma-prompt, gamma-opóźnione, profil energii?) różnicują glębię zwykłą od terenu z ukrytym ładunkiem?

Podsumowanie dydaktyczne

1. Wodór jest dominującym moderatorem neutronów ze względu na dopasowanie masy neutronu i protonu (M ≈ m). Zderzenie czołowe przekazuje prawie 100% energii kinetycznej — to podstawa zarówno moderatorów reaktorowych, jak i neutronowych mierników wilgotności.

2. Miernik neutronowy mierzy efekty rozpraszania neutronów przez wodór — nie mierzy wody jako takiej. Przeliczenie na wilgotność jest możliwe tylko dla znanych materiałów z wystandaryzowaną kalibracją.

3. Kalibracja jest centralnym zagadnieniem metrologicznym: współczynniki a₀, a₁ są specyficzne dla materiału, geometrii i konfiguracji przyrządu. Użycie złej kalibracji daje systematyczny błąd niezależny od precyzji statystycznej.

4. Efekty matrycy (wodór chemicznie związany, pochłaniacze neutronów, niejednorodna gęstość) są głównymi źródłami błędu systematycznego. Szczególnie bor i gadolin mogą zmienić wynik o rząd wielkości nawet przy śladowych stężeniach.

5. Objętość pomiarowa sondy neutronowej nie jest punktem — to sfera o promieniu 10–30 cm, zależna od wilgotności. Gradient wilgotności lub niejednorodność materiału w tej sferze wpływa na wynik pomiaru.

6. Liczba zderzeń potrzebnych do termalizacji jest odwrotnie proporcjonalna do ξ (sredni logarytmiczny ubytek energii). Dla wodoru ξ = 1 → ok. 18 zderzeń. Dla węgla ξ = 0,158 → ok. 115 zderzeń. Dla uranu ξ ≈ 0,0084 → ponad 2000 zderzeń — uran jest praktycznie bez wartości jako moderator.

7. Polska infrastruktura miernikowa obejmuje zastosowania rolnicze (IUNG), hydrologiczne (IMGW), wiertnicze (AGH) i przemysłowe (cementownie). Akademickie kompetencje kalibracyjne utrzymuje NCBJ i uczelnie techniczne.

8. Kluczowe przesłanie dydaktyczne: żaden miernik neutronowy nie działa w izolacji od składu chemicznego materiału. Zrozumienie reakcji ¹H(n,n)¹H, reguły 1/v i efektów matrycy jest warunkiem poprawnej interpretacji wyników — niezależnie od tego, czy mierzymy wilgotność gleby w polu, kruszywo na wytwórni betonów, czy strumień neutronów w rdzeniu reaktora.

Dodatkowe materiały multimedialne

Gotowe narzędzia do ćwiczeń: kalkulator neutronowego miernika wilgotności oraz model 3D neutronowego miernika wilgotności.

Ćwiczenia praktyczne

Pierwsze ćwiczenie: obliczyć E/E0 dla zderzenia czołowego neutronu z jądrem wodoru, węgla i tlenu, używając wzoru ((M - m)/(M + m))^2. Porównać wynik z intuicją bilardową.

Drugie ćwiczenie: korzystając z wartości xi dla kilku moderatorów, policzyć średnią liczbę zderzeń potrzebną do spowolnienia neutronu od 2 MeV do 0,025 eV. Wyniki porównać z tabelą Dziunikowskiego i Kality.

Trzecie ćwiczenie: przygotować syntetyczne punkty kalibracyjne J = a0 + a1 rho_H z dodanym szumem Poissona. Dopasować prostą metodą najmniejszych kwadratów i obliczyć niepewność przewidywanej gęstości wodoru.

Czwarte ćwiczenie: przeliczyć gęstość wodoru na wilgotność masową dla materiału, w którym sucha masa nie zawiera wodoru. Następnie dodać stały wkład wodoru z suchej matrycy i pokazać, jak bardzo zawyża on wynik, jeżeli zostanie pominięty.

Piąte ćwiczenie: zasymulować efekt pochłaniacza neutronów w próbce przez zmniejszenie nachylenia krzywej kalibracyjnej. Student ma porównać wynik uzyskany z kalibracji właściwej i niewłaściwej dla danej matrycy.

Szóste ćwiczenie: zaprojektować minimalną kartę wyniku pomiaru wilgotności metodą neutronową. Powinna zawierać: typ materiału, model kalibracji, czas zliczania, liczbę powtórzeń, tło, gęstość nasypową, wynik, niepewność i ostrzeżenie o wodorze niezwiązanym z wodą.

Przejdź do ćwiczenia interaktywnego

Powiązane artykuły

• Spektrometria neutronowa: kule Bonnera, komory rekombinacyjne i widma neutronów
• Albedo neutronów jako praktyczne pojęcie dyfuzji neutronów
• Ciężka woda D2O jako moderator
• Grafit reaktorowy i zanieczyszczenia borem
• Kalibracja i niepewność pomiaru w laboratorium jądrowym