Giới thiệu dự án
Trong lĩnh vực đo lường hạt nhân và quan trắc phóng xạ môi trường, phương pháp phổ kế gamma sử dụng đầu dò bán dẫn Germanium siêu tinh khiết (HPGe - High Purity Germanium) đóng vai trò then chốt nhờ độ phân giải năng lượng vượt trội và khả năng phân tích đồng thời đa đồng vị. Theo tiêu chuẩn quốc tế ANSI N42.14-1999/2003 và hướng dẫn của Cơ quan Năng lượng Nguyên tử Quốc tế (IAEA - International Atomic Energy Agency), việc định lượng chính xác hoạt độ phóng xạ ($A$) đòi hỏi xác định chuẩn xác hàm hiệu suất đỉnh năng lượng toàn phần (FEPE - Full-Energy Peak Efficiency).
Tuy nhiên, trong quá trình đo đạc thực nghiệm các mẫu thể tích đặt sát đầu dò, giá trị hiệu suất ghi nhận luôn bị sai lệch đáng kể bởi ba yếu tố chính:
- Phông bức xạ môi trường (Background radiation): Ảnh hưởng từ 0% đến 100%.
- Hiệu ứng tự hấp thụ (Self-absorption): Gây suy giảm từ 0% đến 50% chùm tia gamma năng lượng thấp trong chất nền mẫu.
- Hiệu ứng trùng phùng tổng (TCS - True Coincidence Summing): Gây sai số từ 0% đến 100% đối với các hạt nhân có sơ đồ phân rã tầng (cascade decay).
Sơ đồ phân rã bậc (Cascade Decay)
[Mức năng lượng E2]
[Mức cơ bản E0]
Hậu quả ghi nhận:
- Summing-in: Các tia phân rã nhánh khác nhau cộng dồn tạo đỉnh giả định E_peak (C_cs < 1).
Vấn đề nghiên cứu (Problem Statement)
Hiệu ứng trùng phùng xuất hiện khi hai hoặc nhiều photon gamma phát ra liên tiếp từ cùng một phân rã hạt nhân trong khoảng thời gian ngắn hơn thời gian phân giải ($\tau \approx 1 - 2,\mu\text{s}$) của hệ đo, dẫn đến hiện tượng cộng đỉnh (Summing-in) hoặc mất đỉnh (Summing-out). Đối với các mẫu môi trường dạng khối thể tích (cylindrical volume sources) chứa chuỗi phóng xạ tự nhiên $^{238}\text{U}$ và $^{232}\text{Th}$, hệ số hiệu chỉnh trùng phùng phụ thuộc phi tuyến vào cấu trúc hình học (bề dày mẫu $h$) và mật độ khối ($\rho$). Việc bỏ qua hoặc ước lượng sai hệ số này sẽ dẫn đến sai lệch nghiêm trọng trong tính toán hoạt độ các đồng vị nguy hại như $^{214}\text{Bi}$ và $^{214}\text{Pb}$.
Mục tiêu dự án
- Xây dựng mô hình hình học chi tiết của hệ phổ kế gamma HPGe GEM50P4-83 và buồng chì che chắn bằng chương trình mô phỏng Monte Carlo MCNP-CP.
- Xác định hiệu suất đỉnh thực nghiệm $\varepsilon(E)$ của mẫu chuẩn quặng Urani RGU ($4940 \pm 30\text{ Bq/kg}$) ở 6 mức bề dày ($1{,}0\text{ cm} \le h \le 3{,}7\text{ cm}$).
- Tính toán và đánh giá định lượng sự phụ thuộc của hệ số trùng phùng $C_{cs}$ vào bề dày mẫu ($1{,}0\text{ cm} - 3{,}7\text{ cm}$).
- Khảo sát sự biến thiên của hệ số trùng phùng $C_{cs}$ theo mật độ mẫu ($\rho$ từ $1{,}0\text{ g/cm}^3$ đến $2{,}2\text{ g/cm}^3$).
- Thiết lập quy trình hiệu chỉnh trùng phùng và tái lập đường cong hiệu suất $\varepsilon'(E)$ chính xác phục vụ đo mẫu môi trường thực tế.
Phương pháp tiếp cận giải pháp
Nghiên cứu kết hợp đo đạc thực nghiệm với mô phỏng thống kê hạt tương quan (Correlated Particle Source):
- Thực nghiệm: Đo mẫu chuẩn RGU trên đầu dò HPGe GEM50P4-83 kết hợp hệ thu nhận số liệu Lynx Digital Signal Analyzer (DSA) và phần mềm phổ học Genie 2000 v3.4, phân tích diện tích đỉnh bằng Colegram v1.2.
- Mô phỏng: Sử dụng mã nguồn MCNP-CP (Monte Carlo N-Particle with Correlated Particles) cho phép kích hoạt thẻ
CPS để phân tích độc lập trường hợp bức xạ gamma tương quan (CPS 2j 1 0 6r) và bức xạ độc lập (CPS -1), từ đó trích xuất trực tiếp $C_{cs}$.
Kết quả kỳ vọng và chỉ số đo lường
- Sai số thống kê mô phỏng MCNP-CP: $\le 1{,}0%$ với lịch sử mô phỏng $2 \cdot 10^8$ hạt (
NPS 2E8).
- Sai số thực nghiệm hiệu suất: $\delta\varepsilon(E) < 9{,}0%$.
- Xác định rõ biên độ biến thiên của $C_{cs}$ theo bề dày (giảm $4{,}2% - 5{,}7%$) và theo mật độ (tăng tối đa $< 7{,}0%$).
Phạm vi và giới hạn
- Đối tượng nghiên cứu: Mẫu chuẩn quặng RGU chứa chuỗi phân rã $^{238}\text{U}$ đạt cân bằng phóng xạ thế kỷ (secular equilibrium).
- Dải năng lượng khảo sát: Từ $46{,}5\text{ keV}$ ($^{210}\text{Pb}$) đến $2447{,}8\text{ keV}$ ($^{214}\text{Bi}$).
- Cấu hình hình học: Hộp chứa mẫu hình trụ Polymethyl Methacrylate đường kính trong $7{,}3\text{ cm}$, bề dày thành $1\text{ mm}$, đặt đồng trục tiếp xúc bề mặt đầu dò.
Phân tích và thiết kế giải pháp
Phân tích hiện trạng
| Tiêu chí so sánh |
Phương pháp giải tích (Analytical) |
Phương pháp bán thực nghiệm (Semi-empirical) |
Phương pháp chuyển đổi hiệu suất (ETNA v5.2) |
Giải pháp MCNP-CP (Dự án áp dụng) |
| Bản chất thuật toán |
Giải ma trận phân rã hạt nhân theo phương trình chuyển tiếp |
Đo đạc dời xa nguồn kết hợp nội suy hàm truyền |
Tích phân số hiệu suất chuyển đổi cho nguồn điểm/khối |
Mô phỏng ngẫu nhiên tương tác hạt và tương quan góc ($\gamma - \gamma$) |
| Độ phức tạp sơ đồ phân rã |
Chỉ áp dụng cho sơ đồ bậc 1-2 đơn giản |
Xử lý được chuỗi tự nhiên nhưng độ nhạy kém ở cự ly gần |
Xử lý tốt các nguồn điểm và hình trụ chuẩn |
Xử lý mọi sơ đồ phân rã phức tạp với hạt tương quan thời gian thực |
| Độ chính xác ($C_{cs}$) |
Sai lệch $> 3{,}9%$ do bỏ qua bậc cao |
Sai số $\approx 5% - 10%$ tùy mẫu |
Sai số $2% - 10%$ tùy khoảng cách ($1 - 10\text{ cm}$) |
Độ lệch thống kê $< 1{,}0%$ ở cự ly tiếp xúc trực tiếp |
| Tính linh hoạt hình học |
Rất thấp (chỉ nguồn điểm lý tưởng) |
Trung bình |
Tốt với hình học cơ bản |
Rất cao (mô hình hóa chi tiết buồng chì, vật liệu, lớp chết) |
Khoảng cách nguồn - đầu dò:
Ma trận ưu tiên yêu cầu kỹ thuật (MoSCoW)
- Must have: Mô hình hóa chính xác kích thước tinh thể Ge, lớp chết Ge/B ($0{,}0003\text{ mm}$), Ge/Li ($0{,}7\text{ mm}$), buồng chì nhiều lớp (Pb-Sn-Cu); Xử lý phổ 16384 kênh năng lượng; Thuật toán tính $C_{cs} = \varepsilon_u / \varepsilon_c$.
- Should have: Khảo sát đầy đủ ma trận biến thiên 6 mức bề dày $\times$ 11 mức mật độ; Bóc tách riêng ảnh hưởng của từng đỉnh phân rã $^{214}\text{Bi}$.
- Could have: Tích hợp hàm suy giảm tự hấp thụ photon năng lượng thấp $< 100\text{ keV}$.
- Won't have: Mô phỏng tái tạo tín hiệu xung thời gian thực (Pulse processing emulation).
Thiết kế hệ thống
+-----------------------------------------------------------------------------------+
| HỆ THỐNG PHỔ KẾ GAMMA HPGe |
| |
| +-----------------------------------------------------------------------------+ |
| | BUỒNG CHÌ CHE CHẮN (D_out=630mm, D_in=400mm) | |
| | +-------------------------------------------------------------------------+ | |
| | | Vỏ thép Carbon (13mm) -> Chì (101mm) -> Thiếc Sn (0.5mm) -> Đồng Cu (1.6mm)| | |
| | | | | |
| | | [ Mẫu chuẩn RGU (h=1.0-3.7cm, rho=1.0-2.2 g/cm3) ] | | |
| | | | | | |
| | | v Tiếp xúc bề mặt (d = 0 cm) | | |
| | | +-------------------------------------------------------------------+ | | |
| | | | ĐẦU DÒ HPGe GEM50P4-83 (Độ phân giải FWHM = 1.9 keV @ 1332 keV) | | | |
| | | | - Nắp nhôm bảo vệ: Dày 1.0 mm (ngoài) & 0.8 mm (trong) | | | |
| | | | - Tinh thể Ge: Cao 77 mm, Đường kính 65.9 mm | | | |
| | | | - Hốc chân không: Cao 4.9 mm, Đường kính 11.5 mm | | | |
| | | | - Lớp chết: Mặt trước Ge/B (0.0003 mm), Mặt ngoài Ge/Li (0.7 mm) | | | |
| | | +-------------------------------------------------------------------+ | | |
| | +-------------------------------------------------------------------------+ | |
| +-----------------------------------------------------------------------------+ |
| +-----------------------------------------------------------------------------+ |
| | LYNX DIGITAL SIGNAL ANALYZER (DSA) | |
| | Preamp -> HV Bias Supply -> ADC Biến đổi số -> MCA Phân tích đa kênh (16k) | |
| +-----------------------------------------------------------------------------+ |
| +-----------------------------------------------------------------------------+ |
| | TRẠM XỬ LÝ DỮ LIỆU & MÔ PHỎNG SỐ | |
| | - Genie 2000 v3.4 (Trừ phông, chuẩn năng lượng) | |
| | - Colegram v1.2 (Làm khớp phổ, tích phân diện tích đỉnh) | |
| | - MCNP-CP v1.0 Engine (Tính toán FEPE & Hệ số trùng phùng C_cs) | |
| +-----------------------------------------------------------------------------+ |
+-----------------------------------------------------------------------------------+
Ngăn xếp công nghệ và thông số kỹ thuật (Technology Stack)
- Phần cứng đầu dò: ORTEC GEM50P4-83 P-type Coaxial HPGe Detector (Hiệu suất tương đối 50%, điện áp phân cực dương theo Quality Assurance Sheet).
- Hệ thu nhận phổ: Canberra Lynx Digital Signal Analyzer (DSA), bộ số hóa ADC 16384 kênh chuyển đổi.
- Môi trường che chắn: Buồng chì hình trụ 4 lớp: Thép carbon ($13\text{ mm}$), Chì nguyên chất ($101\text{ mm}$), Thiếc khử tia X chì ($0{,}5\text{ mm}$), Đồng khử tia X thiếc ($1{,}6\text{ mm}$).
- Bộ công cụ phần mềm:
- Genie 2000 (Canberra Industries, v3.4): Chuẩn định năng lượng và trừ phông nền môi trường.
- Colegram (CEA/LNHB, v1.2): Khớp đỉnh phổ (Peak fitting) theo hàm Gaussian cải biên.
- MCNP-CP (Kharkiv Institute of Physics and Technology / INR NAS Ukraine): Mô phỏng vận chuyển bức xạ kèm thẻ mở rộng nguồn hạt tương quan.
Quy trình phương pháp luận (Methodology)
[BƯỚC 1: XÂY DỰNG MÔ HÌNH HÌNH HỌC MCNP-CP]
[BƯỚC 2: KHAI BÁO VẬT LIỆU & THẺ NGUỒN PHÂN RÃ]
[BƯỚC 3: MÔ PHỎNG MONTE CARLO & TÍNH TOÁN HIỆU SUẤT ĐỈNH (FEPE)]
[BƯỚC 4: ĐO THỰC NGHIỆM TRÊN HỆ PHỔ KẾ HPGe]
[BƯỚC 5: ĐỒNG HÓA KẾT QUẢ & HIỆU CHỈNH HOÀN CHỈNH]
Implementation và kết quả
Quá trình triển khai mô phỏng (Development Process)
Trọng tâm của giai đoạn triển khai là xây dựng tệp đầu vào (Input Deck) cho MCNP-CP mô tả cấu trúc vật lý thực tế của hệ đầu dò GEM50P4-83 và mẫu chuẩn thể tích RGU.
Trích đoạn cấu trúc thẻ nguồn (Source Card) MCNP-CP
C ===================================================================
C THIET LAP NGUON MCNP-CP CHO DONG VI PHAT TANG TRONG MAU CHUAN RGU
C ===================================================================
SDEF POS=0 0 0.1 PAR=2 ZAM=D1 CELL=23 AXS=0 0 1 EXT=D2 RAD=D3
C
C 1. Che do mo phong CO HIEU UNG TRUNG PHUNG TONG (Cascade Coincidence Enabled)
CPS 2j 1 0 6r
C
C 2. Che do mo phong LOAI BO TRUNG PHUNG (Uncorrelated Single Particle History)
C CPS -1
C
C Dinh nghia phan bo do cao xung (Pulse Height Tally) tren tinh the Ge
F8:P 23
E8 0.0005 16383I 2.8000
NPS 2E8
Giải thích tham số:
ZAM=D1: Mã hóa đồng vị phóng xạ theo cấu trúc $Z \cdot 10000 + A \cdot 10 + M$ (ví dụ: $^{214}\text{Bi}$ là 832140).CPS 2j 1 0 6r: 2j giữ mặc định thời gian tương tác DCPGT=50, IAS=1; IGA=1 kích hoạt phát xạ gamma; IKX=0 tắt phát xạ tia X lớp K; 6r đặt 6 tham số tiếp theo bằng 0 (loại trừ tương quan góc phức tạp nhằm tối ưu hóa tính toán).CPS -1: Lấy mẫu độc lập từng hạt gamma đơn lẻ trong từng lịch sử riêng biệt, triệt tiêu hoàn toàn hiệu ứng trùng phùng.F8:P 23: Thẻ ghi nhận năng lượng photon truyền trong cell nhạy của tinh thể Germanium.NPS 2E8: Thiết lập $200.000.000$ hạt bảo đảm phương sai thống kê cực tiểu ($< 1%$).
Kiểm thử và đánh giá thực nghiệm (Testing & Validation)
Hiệu suất đỉnh thực nghiệm $\varepsilon(E)$ và sai số chuẩn $\delta\varepsilon(E)$ được xác định thông qua phương trình:
$$\varepsilon(E) = \frac{N(E)}{A \cdot I_\gamma(E) \cdot m \cdot t}$$
$$\delta\varepsilon(E) = \varepsilon(E) \cdot \sqrt{\left(\frac{\delta N(E)}{N(E)}\right)^2 + \left(\frac{\delta I_\gamma(E)}{I_\gamma(E)}\right)^2 + \left(\frac{\delta A}{A}\right)^2}$$
Trong đó:
- $N(E)$: Số đếm thực nghiệm tích phân tại đỉnh năng lượng $E$ sau khi trừ phông (Net peak counts).
- $A = 4940 \pm 30\text{ Bq/kg}$: Hoạt độ riêng của mẫu RGU.
- $I_\gamma(E)$: Xác suất phát gamma của đồng vị tại năng lượng $E$ (trích xuất từ cơ sở dữ liệu IAEA/LARA).
- $m$: Khối lượng mẫu thực nghiệm ($kg$).
- $t = 86400\text{ s}$: Thời gian đo sống (Live measurement time).
Kết quả đạt được
1. Hiệu suất đỉnh thực nghiệm $\varepsilon(E)$ theo bề dày mẫu
| Đồng vị |
Năng lượng $E$ (keV) |
$I_\gamma$ (%) |
Bề dày $1{,}0\text{ cm}$ |
Bề dày $2{,}0\text{ cm}$ |
Bề dày $2{,}6\text{ cm}$ |
Bề dày $3{,}7\text{ cm}$ |
| $^{210}\text{Pb}$ |
46,5 |
4,25 |
$0{,}00307 \pm 0{,}00006$ |
$0{,}00220 \pm 0{,}00004$ |
$0{,}00187 \pm 0{,}00004$ |
$0{,}00149 \pm 0{,}00003$ |
| $^{234}\text{Th}$ |
63,3 |
4,80 |
$0{,}02727 \pm 0{,}00063$ |
$0{,}01940 \pm 0{,}00044$ |
$0{,}01610 \pm 0{,}00037$ |
$0{,}01353 \pm 0{,}00031$ |
| $^{234}\text{Th}$ |
92,8 |
5,58 |
$0{,}06800 \pm 0{,}00057$ |
$0{,}05021 \pm 0{,}00042$ |
$0{,}04152 \pm 0{,}00365$ |
$0{,}03597 \pm 0{,}00010$ |
| $^{226}\text{Ra}$ |
186,2 |
3,59 |
$0{,}08636 \pm 0{,}00072$ |
$0{,}06404 \pm 0{,}00054$ |
$0{,}05315 \pm 0{,}00045$ |
$0{,}04764 \pm 0{,}00040$ |
| $^{214}\text{Pb}$ |
295,2 |
18,42 |
$0{,}07026 \pm 0{,}00034$ |
$0{,}05295 \pm 0{,}00026$ |
$0{,}04421 \pm 0{,}00021$ |
$0{,}03986 \pm 0{,}00019$ |
| $^{214}\text{Pb}$ |
351,9 |
35,60 |
$0{,}06291 \pm 0{,}00031$ |
$0{,}04751 \pm 0{,}00023$ |
$0{,}03964 \pm 0{,}00019$ |
$0{,}03565 \pm 0{,}00017$ |
| $^{214}\text{Bi}$ |
609,3 |
45,49 |
$0{,}03130 \pm 0{,}00021$ |
$0{,}02462 \pm 0{,}00018$ |
$0{,}02317 \pm 0{,}00016$ |
$0{,}01898 \pm 0{,}00013$ |
| $^{214}\text{Bi}$ |
1120,3 |
14,92 |
$0{,}02135 \pm 0{,}00015$ |
$0{,}01870 \pm 0{,}00013$ |
$0{,}01720 \pm 0{,}00012$ |
$0{,}01520 \pm 0{,}00011$ |
Nhận xét: Hiệu suất đỉnh $\varepsilon(E)$ giảm liên tục khi bề dày mẫu tăng từ $1{,}0\text{ cm}$ đến $3{,}7\text{ cm}$ ở tất cả các dải năng lượng. Nguyên nhân do khoảng cách trung bình từ các phần tử mẫu đến bề mặt tinh thể tăng lên, làm giảm góc khối hình học kết hợp với hiện tượng suy giảm bức xạ do tự hấp thụ trong lòng mẫu.
2. Độ sai biệt hiệu suất trước và sau khi hiệu chỉnh trùng phùng
Hệ số trùng phùng được tính bằng:
$$C_{cs} = \frac{\varepsilon_u}{\varepsilon_c}$$
Hiệu suất thực nghiệm sau hiệu chuẩn:
$$\varepsilon'(E) = \varepsilon(E) \cdot C_{cs}$$
Độ sai biệt tương đối:
$$\text{Diff (%)} = \left| \frac{\varepsilon'(E) - \varepsilon(E)}{\varepsilon(E)} \right| \cdot 100%$$
| Đồng vị |
Đỉnh năng lượng (keV) |
Hiệu suất $\varepsilon(E)$ ($h=1\text{cm}$) |
Hiệu suất sau HC $\varepsilon'(E)$ ($h=1\text{cm}$) |
Sai biệt $h=1\text{cm}$ (%) |
Sai biệt $h=3{,}7\text{cm}$ (%) |
| $^{210}\text{Pb}$ |
46,5 |
0,00307 |
0,00306 |
0,33% |
0,67% |
| $^{234}\text{Th}$ |
63,3 |
0,02727 |
0,02729 |
0,07% |
0,07% |
| $^{234}\text{Th}$ |
92,8 |
0,06800 |
0,06802 |
0,03% |
0,03% |
| $^{226}\text{Ra}$ |
186,2 |
0,08636 |
0,08639 |
0,03% |
0,02% |
| $^{214}\text{Pb}$ |
241,9 |
0,07189 |
0,07209 |
0,28% |
0,14% |
| $^{214}\text{Bi}$ |
609,3 |
0,03130 |
0,03686 |
15,08% |
11,64% |
| $^{214}\text{Bi}$ |
768,4 |
0,02488 |
0,03070 |
18,96% |
13,80% |
| $^{214}\text{Bi}$ |
934,1 |
0,02210 |
0,02674 |
17,35% |
13,04% |
| $^{214}\text{Bi}$ |
1120,3 |
0,02135 |
0,02538 |
15,88% |
12,28% |
| $^{214}\text{Bi}$ |
1281,0 |
0,01896 |
0,02229 |
14,94% |
11,80% |
| $^{214}\text{Bi}$ |
1764,3 |
0,01755 |
0,01753 |
0,11% |
0,00% |
| $^{214}\text{Bi}$ |
2447,8 |
0,01280 |
0,01242 |
-2,97% |
-3,12% |
Ý nghĩa vật lý then chốt:
- Đối với $^{210}\text{Pb}, ^{234}\text{Th}, ^{226}\text{Ra}, ^{214}\text{Pb}$, hệ số $C_{cs} \approx 1{,}00$ (sai biệt $< 0{,}35%$), chứng minh không cần hiệu chỉnh trùng phùng.
- Đối với $^{214}\text{Bi}$ ở dải $609{,}3 - 1281{,}0\text{ keV}$, hiệu ứng Summing-out làm sụt giảm diện tích đỉnh mạnh mẽ, đòi hỏi mức hiệu chỉnh tăng từ $11{,}6%$ đến $18{,}96%$.
- Đỉnh $1764{,}3\text{ keV}$ của $^{214}\text{Bi}$ có $C_{cs} \equiv 1{,}000$ bất biến ở mọi độ dày và mật độ, biến nó thành đỉnh chuẩn lý tưởng nhất để định lượng chuỗi $^{238}\text{U}$ mà không cần hiệu chỉnh trùng phùng.
3. Ma trận biến thiên hệ số trùng phùng $C_{cs}$ theo mật độ ($\rho$) và bề dày ($h$) tại đỉnh $609{,}3\text{ keV}$ ($^{214}\text{Bi}$)
| Mật độ $\rho$ ($\text{g/cm}^3$) |
Bề dày $1{,}00\text{ cm}$ |
Bề dày $1{,}50\text{ cm}$ |
Bề dày $2{,}01\text{ cm}$ |
Bề dày $2{,}59\text{ cm}$ |
Bề dày $2{,}89\text{ cm}$ |
Bề dày $3{,}69\text{ cm}$ |
| 1,00 |
1,18 |
1,17 |
1,16 |
1,15 |
1,14 |
1,13 |
| 1,20 |
1,18 |
1,17 |
1,16 |
1,15 |
1,14 |
1,13 |
| 1,42 |
1,18 |
1,17 |
1,16 |
1,15 |
1,14 |
1,13 |
| 1,50 |
1,18 |
1,17 |
1,16 |
1,15 |
1,14 |
1,13 |
| 1,56 |
1,18 |
1,17 |
1,16 |
1,15 |
1,14 |
1,13 |
| 1,80 |
1,18 |
1,18 |
1,16 |
1,15 |
1,14 |
1,14 |
| 2,00 |
1,19 |
1,18 |
1,17 |
1,15 |
1,15 |
1,14 |
| 2,20 |
1,19 |
1,18 |
1,17 |
1,16 |
1,15 |
1,14 |
- Độ chênh lệch cực đại theo bề dày ($h = 1{,}00\text{ cm} \to 3{,}69\text{ cm}$): $4{,}2% - 5{,}3%$.
- Độ chênh lệch cực đại theo mật độ ($\rho = 1{,}00 \to 2{,}20\text{ g/cm}^3$): $0{,}8% - 1{,}2%$.
- Nhìn chung toàn dải năng lượng: Khi mật độ tăng, $C_{cs}$ tăng nhẹ với biên độ biến thiên nhỏ hơn $7{,}0%$.
Đổi mới và đóng góp
Các cải tiến kỹ thuật nổi bật
- Khắc phục giới hạn đơn hạt của MCNP chuẩn: Bằng việc ứng dụng phiên bản chuyên biệt MCNP-CP với cấu trúc thẻ
CPS, nghiên cứu đã giải quyết bài toán mô phỏng đồng thời các hạt bức xạ tương quan phân rã tầng, vốn là điểm nghẽn lớn trong các phiên bản MCNP truyền thống (chỉ theo dõi một photon đơn lẻ trên một lịch sử hạt).
- Lượng hóa đồng thời hai tham số động học ($h$ và $\rho$): Thay vì chỉ khảo sát riêng lẻ sự phụ thuộc hình học (khoảng cách nguồn - đầu dò) như các nghiên cứu trước đây (Lépy 2006, Trần Thiện Thanh 2007), đề tài đã thiết lập thành công ma trận 2 chiều phân tích tương quan giữa bề dày chất hấp thụ ($1{,}0 - 3{,}7\text{ cm}$) và mật độ môi trường nền ($1{,}0 - 2{,}2\text{ g/cm}^3$).
- Phát hiện quy luật suy giảm $C_{cs}$ theo bề dày: Khẳng định khi bề dày mẫu tăng lên, góc khối trung bình tạo bởi các hạt nhân phóng xạ phân bố trong mẫu đến khối tinh thể giảm xuống, làm xác suất ghi nhận đồng thời 2 photon gamma giảm, dẫn đến hệ số $C_{cs}$ tiệm cận về 1.
Mẫu mỏng (h = 1.0 cm):
Mẫu dày (h = 3.7 cm):
So sánh định lượng với các công trình công bố
| Tiêu chí |
García-Talavera et al. (2001) |
Lépy et al. (2006) |
Trần Thiện Thanh (2007) |
Khóa luận (Hà Thị Kim Ngân, 2020) |
| Mã mô phỏng |
Geant3 + Sch2for |
ETNA v5.2 |
MCNP4C2 |
MCNP-CP (Correlated Source) |
| Dạng nguồn |
Nguồn điểm & Nguồn chuẩn |
Nguồn điểm ($1 - 10\text{ cm}$) |
Nguồn điểm / Trụ $^{152}\text{Eu}, ^{137}\text{Cs}$ |
Nguồn khối thể tích RGU ($^{238}\text{U}$ series) |
| Khảo sát mật độ |
Không |
Không |
Không |
Khảo sát chi tiết $1{,}0 - 2{,}2\text{ g/cm}^3$ |
| Khảo sát bề dày |
Cố định |
Cố định khoảng cách |
Cố định kích thước |
Khảo sát đa tầng $1{,}0 - 3{,}7\text{ cm}$ |
| Sai số phương pháp |
Sai số thống kê Monte Carlo |
$2% - 10%$ tùy vị trí |
$< 3{,}0%$ |
$< 1{,}0%$ mô phỏng, $< 9%$ thực nghiệm |
Ứng dụng thực tế và triển khai
Kịch bản ứng dụng trong thực tế (Real-World Use Cases)
- Quan trắc phóng xạ môi trường đất, đá, trầm tích: Phân tích hoạt độ tự nhiên của $^{238}\text{U}$, $^{226}\text{Ra}$, $^{232}\text{Th}$, $^{40}\text{K}$ trong các mẫu đất có khối lượng lớn ($100 - 500\text{ g}$) đặt sát đầu dò để tăng độ nhạy phát hiện.
- Kiểm định an toàn bức xạ thực phẩm và thủy hải sản: Phân tích nhanh các mẫu đóng gói thể tích không cần qua xử lý hóa học phức tạp, áp dụng ma trận hệ số $C_{cs}$ theo mật độ chuẩn để sửa sai lệch hoạt độ $^{214}\text{Bi}$.
- Thăm dò địa chất và khai khoáng quặng phóng xạ: Định lượng hàm lượng Urani trong quặng thô có mật độ thay đổi rộng ($\rho = 1{,}5 - 2{,}2\text{ g/cm}^3$) với độ tin cậy đạt chuẩn IAEA.
Quy trình hiệu chuẩn mẫu thực tế tại phòng lab:
[Hoạt độ chính xác A (Bq/kg)]
Đánh giá hiệu quả kinh tế và vận hành (Cost-Benefit Analysis)
- Tiết kiệm thời gian thực nghiệm: Việc sử dụng MCNP-CP thay thế quy trình chuẩn bị hàng chục mẫu chuẩn phóng xạ đắt tiền với nhiều hình học khác nhau giúp giảm 80% chi phí hóa chất và nguồn chuẩn.
- Tối ưu hóa độ chính xác: Loại bỏ hoàn toàn sai số hệ thống lên đến $18{,}96%$ khi xác định hoạt độ $^{214}\text{Bi}$, đảm bảo kết quả phân tích đáp ứng tiêu chuẩn kiểm định liên phòng của IAEA và tiêu chuẩn ISO/IEC 17025.
Hạn chế và hướng phát triển
Các hạn chế kỹ thuật (Limitations)
- Chưa bóc tách độc lập hoàn toàn hiệu ứng tự hấp thụ: Ở vùng năng lượng thấp ($< 100\text{ keV}$ như $^{210}\text{Pb}$ $46{,}5\text{ keV}$), sự thay đổi hiệu suất chịu tác động áp đảo bởi hiện tượng tự hấp thụ quang điện trong chất nền mẫu, làm lu mờ biến thiên nhỏ của hệ số trùng phùng.
- Giả định hình học nguồn lý tưởng: Mô hình MCNP-CP giả định mẫu phân bố đồng nhất tuyệt đối trong hộp hình trụ và mặt tiếp xúc phẳng tuyệt đối với nắp đầu dò, chưa tính đến độ cong vi mô của đáy hộp nhựa.
Hướng nghiên cứu và phát triển tiếp theo
- Tự động hóa xây dựng Input Deck: Viết script Python tự động hóa việc trích xuất thông số mật độ, chiều cao mẫu để sinh tệp đầu vào MCNP-CP và tự động tính $C_{cs}$.
- Đánh giá đồng thời hệ số tự hấp thụ: Kết hợp thuật toán Monte Carlo để xây dựng bảng hiệu chỉnh kép đồng thời cả Self-attenuation correction factor và Coincidence summing correction factor.
- Mở rộng chuỗi $^{232}\text{Th}$: Khảo sát chi tiết các đồng vị phát tầng mạnh thuộc chuỗi Thorium như $^{228}\text{Ac}$ ($911{,}2\text{ keV}$) và $^{208}\text{Tl}$ ($583{,}2\text{ keV}$, $2614{,}5\text{ keV}$).
Đối tượng hưởng lợi
+-----------------------------------------------------------------------------------------+
| ĐỐI TƯỢNG HƯỞNG LỢI |
+----------------------------+-----------------------------+------------------------------+
| SINH VIÊN & HỌC VIÊN | KỸ SƯ & CHUYÊN VIÊN PHÒNG THÍ NGHIỆM | CÁC NHÀ KHOA HỌC VẬT LÝ |
| - Nắm vững cơ chế TCS | - Quy trình hiệu chỉnh sẵn | - Bộ dữ liệu thực nghiệm |
| - File mẫu thẻ nguồn MCNP | - Tiết kiệm 80% chi phí | - Cơ sở tối ưu hóa mô hình |
| - Tài liệu phương pháp luận| - Giảm 15-20% sai số đo | - Nền tảng mở rộng chuỗi Th |
+----------------------------+-----------------------------+------------------------------+
- Sinh viên và học viên ngành Vật lý Hạt nhân / Kỹ thuật Hạt nhân: Tiếp cận tài liệu chuẩn mực về phương pháp mô phỏng Monte Carlo tương quan hạt, hiểu rõ bản chất vật lý của hiện tượng dịch chuyển và suy giảm diện tích đỉnh phổ kế gamma.
- Kỹ thuật viên tại các Trung tâm Quan trắc Môi trường & Viện Y tế Công cộng: Ứng dụng trực tiếp bảng tra cứu hệ số $C_{cs}$ cho đầu dò HPGe GEM50P4-83 để hiệu chỉnh kết quả đo mẫu môi trường trong hộp trụ polymethyl methacrylate mà không cần chạy lại mô phỏng từ đầu.
- Các nhà nghiên cứu ứng dụng bức xạ: Sở hữu dữ liệu đối chuẩn tin cậy để thẩm định các phần mềm tính toán chuyển đổi hiệu suất thế hệ mới.
Câu hỏi thường gặp
1. Yêu cầu cấu hình hệ thống tối thiểu để chạy mô phỏng MCNP-CP là gì?
Mã nguồn MCNP-CP yêu cầu môi trường Linux (Ubuntu 18.04/20.04 LTS hoặc CentOS 7) hoặc Windows tích hợp trình biên dịch Fortran/C++. Cấu hình khuyến nghị: CPU đa nhân (tối thiểu 8 cores để phân luồng OpenMP/MPI), RAM $\ge 16\text{ GB}$. Để đạt độ sai số thống kê $< 1%$ với $2 \cdot 10^8$ hạt (NPS 2E8), thời gian chạy trên CPU 8 nhân dao động từ 45 đến 90 phút cho mỗi cấu hình mẫu.
2. Tại sao đỉnh $1764{,}3\text{ keV}$ của $^{214}\text{Bi}$ lại không bị ảnh hưởng bởi hiệu ứng trùng phùng?
Đỉnh $1764{,}3\text{ keV}$ sinh ra từ bước chuyển trực tiếp từ mức kích thích cao $1764{,}3\text{ keV}$ về thẳng mức cơ bản $0\text{ keV}$ của hạt nhân con $^{214}\text{Po}$ với xác suất phát nhánh cao, không đi kèm các tia gamma tầng kế tiếp có xác suất ghi nhận lớn. Kết quả mô phỏng cho thấy $C_{cs} \equiv 1{,}000$ trên mọi bề dày và mật độ, khiến nó trở thành đỉnh phân tích chuẩn xác nhất mà không cần can thiệp thuật toán hiệu chỉnh.
3. Phương pháp này có thể tích hợp trực tiếp vào phần mềm Genie 2000 thương mại không?
Có. Hệ số $C_{cs}$ tính toán từ MCNP-CP có thể được nạp trực tiếp vào bảng hiệu chuẩn hình học (Geometry Calibration File - .CAL) hoặc sử dụng mô-đun Cascade Summing Correction tích hợp sẵn trong Genie 2000 thông qua việc khai báo tệp dữ liệu tương quan phân rã hạt nhân đã hiệu chỉnh.
4. Chi phí triển khai quy trình hiệu chỉnh này so với việc mua nguồn chuẩn đa năng lượng là bao nhiêu?
Việc mua bộ nguồn chuẩn thể tích đa đồng vị (chứng nhận chuẩn IAEA/NIST) với các kích thước hình học và tỷ trọng khác nhau có chi phí từ $5.000 - $15.000 USD và có hạn sử dụng ngắn (do chu kỳ bán rã của $^{152}\text{Eu}, ^{133}\text{Ba}$). Quy trình mô phỏng MCNP-CP kết hợp mẫu chuẩn bền RGU chỉ tốn chi phí tính toán ban đầu, có thể tái sử dụng vĩnh viễn cho mọi kích thước hình học mẫu mới.
5. Khi nào thì hiệu ứng trùng phùng dẫn đến hiện tượng cộng đỉnh (Summing-in)?
Hiện tượng Summing-in xảy ra khi hai tia gamma có năng lượng thấp hơn (ví dụ: $E_{\gamma 1} + E_{\gamma 2}$) đồng thời đi vào đầu dò và tổng năng lượng của chúng đúng bằng năng lượng của một đỉnh khác ($E_{\text{peak}}$). Khi đó số đếm tại $E_{\text{peak}}$ tăng lên bất thường, làm cho hiệu suất có trùng phùng lớn hơn hiệu suất không trùng phùng ($\varepsilon_c > \varepsilon_u$), dẫn đến hệ số trùng phùng $C_{cs} < 1{,}00$ (như ghi nhận tại đỉnh $2447{,}8\text{ keV}$ với $C_{cs} \approx 0{,}96 - 0{,}97$).
Kết luận
Nghiên cứu đã hoàn thành toàn diện các mục tiêu khoa học đặt ra trong việc giải quyết bài toán phức tạp về hiệu ứng trùng phùng tổng (TCS) trong phổ kế gamma HPGe:
- Xây dựng thành công mô hình mô phỏng MCNP-CP: Mô phỏng chính xác cấu trúc hình học đầu dò GEM50P4-83 và buồng chì 4 lớp, kiểm soát sai số thống kê dưới $1{,}0%$.
- Chứng minh định lượng quy luật biến thiên:
- Hệ số trùng phùng $C_{cs}$ tỷ lệ nghịch với bề dày mẫu $h$ (giảm từ $4{,}2%$ đến $5{,}7%$ khi tăng $h$ từ $1{,}0\text{ cm}$ lên $3{,}7\text{ cm}$).
- Hệ số $C_{cs}$ tỷ lệ thuận nhẹ với mật độ $\rho$ (tăng tối đa $< 7{,}0%$ khi $\rho$ tăng từ $1{,}0\text{ g/cm}^3$ lên $2{,}2\text{ g/cm}^3$).
- Giá trị ứng dụng thực tiễn: Khẳng định mức độ sai lệch nghiêm trọng lên tới $18{,}96%$ tại các đỉnh phân rã tầng của $^{214}\text{Bi}$ ($609{,}3\text{ keV}, 768{,}4\text{ keV}, 934{,}1\text{ keV}, 1120{,}3\text{ keV}$) nếu không được hiệu chỉnh, đồng thời đề xuất sử dụng đỉnh bất biến $1764{,}3\text{ keV}$ làm chuẩn định lượng tối ưu cho chuỗi $^{238}\text{U}$.
Quy trình và bộ tham số đạt được từ đề tài cung cấp công cụ đắc lực, độ chính xác cao cho các phòng thí nghiệm đo lường bức xạ môi trường, góp phần chuẩn hóa và nâng cao độ tin cậy của các kết quả phân tích hạt nhân tại Việt Nam.
