Giới thiệu dự án

Trong lĩnh vực quan trắc phóng xạ môi trường, y học hạt nhân và công nghiệp năng lượng nguyên tử, hệ phổ kế gamma sử dụng đầu dò nhấp nháy vô cơ NaI(Tl) (Sodium Iodide doped with Thallium) đóng vai trò then chốt nhờ hiệu suất ghi cao và khả năng vận hành ổn định ở nhiệt độ phòng. Tuy nhiên, việc xác định chính xác đường cong hiệu suất đỉnh năng lượng toàn phần (Full Energy Peak Efficiency - FEPE, $\varepsilon_{FEPE}$) đối với các hình thái mẫu và khoảng cách đo khác nhau bằng phương pháp thực nghiệm truyền thống đòi hỏi nguồn chuẩn phóng xạ đa dạng, chi phí tốn kém và làm tăng liều chiếu xạ cho kiểm định viên.

       +-------------------------------------------------------------+
       |                  NGUỒN PHÁT GAMMA (Eγ)                      |
       +-------------------------------------------------------------+
                                      |
                                      v (Góc khối Ω, khoảng cách D)
       +-------------------------------------------------------------+
       |                  ĐẦU DÒ NHẤP NHÁY NaI(Tl)                   |
       |  +-------------------------------------------------------+  |
       |  | Vỏ nhôm bảo vệ bên ngoài (Al, dày 0.05 - 0.16 cm)     |  |
       |  | +---------------------------------------------------+ |  |
       |  | | Đệm không khí / Lớp phản xạ Al2O3                 | |  |
       |  | | +-----------------------------------------------+ | |  |
       |  | | | Tinh thể NaI(Tl) (Kích thước: R x L)          | | |  |
       |  | | | Hấp thụ photon: Quang điện / Compton / Tạo cặp| | |  |
       |  | | +-----------------------------------------------+ | |  |
       |  | +---------------------------------------------------+ |  |
       |  +-------------------------------------------------------+  |
       +-------------------------------------------------------------+
                                      |
                                      v (Phát photon ánh sáng nhấp nháy)
       +-------------------------------------------------------------+
       |      ỐNG NHÂN QUANG ĐIỆN (PMT) & BỘ PHÂN TÍCH ĐA KÊNH (MCA) |
       +-------------------------------------------------------------+

Vấn đề nghiên cứu (Problem Statement)

Mô phỏng Monte Carlo bằng chương trình MCNP5 (Monte Carlo N-Particle version 5) cho phép tính toán hiệu suất ghi của đầu dò với độ linh hoạt cao. Tuy nhiên, khi sử dụng trực tiếp các thông số kỹ thuật do nhà sản xuất cung cấp (kích thước hình học tinh thể, mật độ và độ dày lớp phản xạ $Al_2O_3$, lớp vỏ bọc), kết quả mô phỏng thường sai lệch đáng kể so với thực nghiệm (độ lệch có thể lên tới 19.74% ở vùng năng lượng thấp như 31 keV). Nguyên nhân xuất phát từ sự không đồng nhất trong chế tạo, quá trình lão hóa của vật liệu và lớp chết (dead layer) không được định lượng chính xác.

Mục tiêu dự án

  1. Xây dựng phương pháp thực nghiệm kết hợp mô phỏng MCNP5 nhằm tách rời và tối ưu hóa độc lập 3 thông số kỹ thuật cốt lõi của đầu dò NaI(Tl): Mật độ lớp phản xạ $Al_2O_3$ ($\rho_{Al_2O_3}$), Bán kính tinh thể ($R$), và Chiều dài tinh thể ($L$).
  2. Thiết lập quy trình chuẩn trực và lựa chọn dải năng lượng photon đặc trưng (từ 31 keV đến 1408 keV) để cô lập ảnh hưởng của từng thông số hình học.
  3. Giảm độ lệch tương đối giữa hiệu suất mô phỏng và thực nghiệm xuống dưới mức 2.5% trên toàn bộ dải năng lượng khảo sát (31 keV – 1408 keV).

Phương pháp tiếp cận & Kết quả kỳ vọng

Dự án áp dụng nguyên lý suy giảm bức xạ truyền qua và phân bố năng lượng theo độ sâu tương tác để tách rời biến số:

  • Khảo sát mật độ $\rho_{Al_2O_3}$: Dùng chùm tia chuẩn trực hẹp năng lượng thấp 31 keV từ nguồn $^{133}\text{Ba}$ chiếu vào mặt trước và mặt bên.
  • Khảo sát bán kính $R$: Dùng nguồn điểm năng lượng thấp (31 – 121 keV) không chuẩn trực đặt ở khoảng cách $D = 40\text{ cm}$.
  • Khảo sát chiều dài $L$: Dùng nguồn điểm năng lượng cao (662 – 1408 keV) có độ xuyên sâu lớn.

Phạm vi nghiên cứu giới hạn trên hệ phổ kế NaI(Tl) kích thước danh định $76.2 \times 76.2\text{ mm}$ ($3 \times 3\text{ inch}$) kết nối bộ phân tích đa kênh MCA, nguồn bức xạ chuẩn $^{133}\text{Ba}, ^{137}\text{Cs}, ^{241}\text{Am}, ^{152}\text{Eu}, ^{60}\text{Co}, ^{22}\text{Na}$.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Phương pháp Ưu điểm Nhược điểm Độ lệch thực nghiệm
Chuẩn trực thực nghiệm toàn phần Độ tin cậy cao, trực quan Tốn kém nguồn chuẩn, mất thời gian, liều xạ cao Chuẩn cơ sở (0%)
Mô phỏng MCNP với thông số NSX Nhanh, không cần nguồn chuẩn Sai số rất lớn ở vùng năng lượng thấp và mép tinh thể 15% – 20%
Hiệu chỉnh đồng thời đa biến (Fit tổng) Khớp tốt một số đỉnh Dễ rơi vào cực tiểu cục bộ, thông số không có ý nghĩa vật lý 3% – 5%
Quy trình tách rời thông số (Đề xuất) Chính xác, thông số vật lý rõ ràng, áp dụng được cho mọi dải năng lượng Cần thiết lập bố trí đo chuyên biệt (ống chuẩn trực) < 2.16%

Thiết kế hệ thống mô phỏng và tham số đầu vào

Cấu trúc hình học đầu dò NaI(Tl) và hệ đo được mô hình hóa trong MCNP5 thông qua hệ thống thẻ ô mạng (Cell Cards), mặt (Surface Cards) và dữ liệu (Data Cards).

          Trục đối xứng Oz
                 |
                 +========================+ --- Mặt trước đầu dò
                 |  Vỏ nhôm (Al: 0.05cm)  |
                 |  +------------------+  | --- Lớp đệm không khí (0.16cm)
                 |  | Phản xạ Al2O3     |  | --- Lớp phản xạ (d_Al2O3)
                 |  | +--------------+ |  |
                 |  | |              | |  |
                 |  | | Tinh thể     | |  |
                 |  | | NaI(Tl)      | |  | Bán kính R
                 |  | |              | |  | Chiều dài L
                 |  | |              | |  |
                 |  | +--------------+ |  |
                 |  +------------------+  |
                 +========================+

C =================================================================
C FILE INPUT MCNP5: MO HINH HOA DAU DO NaI(Tl) TOI UU
C =================================================================
MESSAGE: Detector NaI(Tl) FEPE Simulation
-- Title Card: Optimization of NaI(Tl) detector parameters
C --- CELL CARDS ---
1    1 -3.667   -1 2 -3            $ Cell 1: Tinh the NaI(Tl)
2    2 -2.020   (-2:1:3) -4 5 -6   $ Cell 2: Lop phan xa Al2O3 toi uu
3    3 -0.0012  (-5:4:6) -7 8 -9   $ Cell 3: Lop dem khong khi
4    4 -2.700   (-8:7:9) -10 11 -12 $ Cell 4: Vo nhom (Al) ngoai
5    0          #1 #2 #3 #4 -99    $ Cell 5: Moi truong khong khi ngoai
6    0          99                 $ Cell 6: Ngoai the gioi (Graveyard)

C --- SURFACE CARDS ---
1    CZ  3.76                      $ Ban kinh tinh the toi uu (cm)
2    PZ  0.0                       $ Mat truoc tinh the
3    PZ  -7.55                     $ Mat sau tinh the
4    CZ  3.92                      $ Mat ngoai lop phan xa Al2O3
5    PZ  0.16                      $ Mat truoc lop phan xa
6    PZ  -7.71                     $ Mat sau lop phan xa
10   CZ  4.08                      $ Mat ngoai vo Al
11   PZ  0.24                      $ Mat truoc vo Al
12   PZ  -7.90                     $ Mat sau vo Al
99   SO  100.0                     $ Mat cau gioi han bai toan

C --- DATA CARDS ---
MODE P
SDEF POS=0 0 40.24 AXS=0 0 -1 RAD=D1 ERG=0.031 PAR=2
SI1  0 0.1
SP1 -21 1
C --- TALLY F8: Ghi nhan pho nang luong de lai trong Cell 1 ---
F8:P 1
E8   0.001 1000I 2.0
C --- MATERIAL CARDS ---
M1   11023 0.5 53127 0.5           $ NaI (ti le nguyen tu)
M2   13027 2.0 8016  3.0           $ Al2O3
M4   13027 1.0                     $ Al tinh khiet
NPS  10000000                      $ So hat chay: 10^7 photon

Phương pháp luận (Methodology)

Quy trình quản lý chất lượng dữ liệu mô phỏng dựa trên 10 bài kiểm tra thống kê tiêu chuẩn của MCNP5 (NPS $\ge 10^7$ hạt, đảm bảo sai số tương đối $R_{error} < 0.01$ tại các đỉnh phổ).

Implementation và kết quả

Cơ sở toán học và thuật toán xử lý

1. Định luật tương tác bức xạ và xấp xỉ suy giảm chùm tia hẹp

Cường độ chùm photon đơn năng truyền qua môi trường vật chất tuân theo định luật Beer-Lambert: $$I = I_0 e^{-\mu d} = I_0 e^{-\left(\frac{\mu}{\rho}\right) \rho d}$$

Đối với năng lượng thấp ($E_\gamma = 31\text{ keV}$ của $^{133}\text{Ba}$), hiệu ứng quang điện chiếm ưu thế tuyệt đối. Hiệu suất đỉnh toàn phần $\varepsilon_{FEPE}$ phụ thuộc vào mật độ lớp phản xạ theo hàm mũ: $$\varepsilon_{FEPE} = C \cdot e^{-\left(\frac{\mu}{\rho}\right){Al_2O_3} \cdot d{Al_2O_3} \cdot \rho_{Al_2O_3}}$$

Khai triển MacLaurin bậc nhất cho số mũ nhỏ ($\mu_{Al_2O_3} \cdot d_{Al_2O_3} \ll 1$): $$\varepsilon_{FEPE} \approx \alpha - \beta \cdot \rho_{Al_2O_3}$$

2. Hiệu suất đỉnh toàn phần thực nghiệm

$$\varepsilon_{FEPE} = \frac{N_{peak}}{A_0 \cdot e^{-\lambda t} \cdot I_\gamma \cdot t_{live}}$$

Trong đó:

  • $N_{peak}$: Diện tích đỉnh thuần (net peak area) sau khi trừ phông liên tục Compton.
  • $A_0$: Hoạt độ ban đầu của nguồn chuẩn tại ngày chứng nhận kiểm định.
  • $\lambda = \frac{\ln(2)}{T_{1/2}}$: Hằng số phân rã của đồng vị phóng xạ.
  • $I_\gamma$: Xác suất phát photon trên mỗi phân rã.
  • $t_{live}$: Thời gian đo sống (live time, giây).

3. Phương pháp trung bình trọng số xác định bán kính $R$ tối ưu

Giá trị bán kính tối ưu $\bar{R}$ và sai số lan truyền $\sigma_{\bar{R}}$ được tổng hợp từ $k$ đỉnh năng lượng độc lập: $$\bar{R} = \frac{\sum_{i=1}^{k} w_i R_i}{\sum_{i=1}^{k} w_i}, \quad w_i = \frac{1}{\sigma_i^2}, \quad \sigma_{\bar{R}} = \sqrt{\frac{1}{\sum_{i=1}^{k} w_i}}$$

Dữ liệu thực nghiệm và kết quả tối ưu hóa

Thí nghiệm 1: Xác định mật độ tối ưu lớp phản xạ $Al_2O_3$

Khảo sát chùm photon 31 keV với mật độ $\rho_{Al_2O_3}$ quét từ $0.4\text{ g/cm}^3$ đến $3.6\text{ g/cm}^3$:

  • Phương trình hồi quy mặt trước: $\varepsilon = (6.841 \times 10^{-4}) - (0.598 \times 10^{-4}) \cdot \rho \quad (R^2 = 0.9995)$
  • Phương trình hồi quy mặt bên: $\varepsilon = (6.611 \times 10^{-4}) - (0.564 \times 10^{-4}) \cdot \rho \quad (R^2 = 0.9977)$
  • Giá trị nội suy thực nghiệm:
    • Mặt trước: $\rho = 2.05 \pm 0.31\text{ g/cm}^3$
    • Mặt bên: $\rho = 1.99 \pm 0.33\text{ g/cm}^3$
  • Mật độ tối ưu trung bình: $\bar{\rho}_{Al_2O_3} = 2.02 \pm 0.32\text{ g/cm}^3$ (thay vì $0.55\text{ g/cm}^3$ do NSX công bố).

Thí nghiệm 2: Xác định bán kính tối ưu tinh thể NaI(Tl)

Khảo sát sự thay đổi $\varepsilon_{FEPE}$ theo $R$ trong khoảng $3.72\text{ cm} \le R \le 3.82\text{ cm}$:

Đồng vị Năng lượng (keV) $\varepsilon_{Exp} \times 10^{-4}$ $\varepsilon_{Sim}^{NSX} \times 10^{-4}$ Sai lệch ban đầu $R$ tối ưu nội suy (cm) Sai số $\sigma_R$ (cm)
$^{133}\text{Ba}$ 31 16.40 19.61 +19.57% 3.74 0.07
$^{137}\text{Cs}$ 32 16.64 19.88 +19.47% 3.75 0.06
$^{241}\text{Am}$ 59 19.48 20.25 +3.95% 3.77 0.06
$^{133}\text{Ba}$ 81 17.52 18.22 +3.99% 3.74 0.06
$^{152}\text{Eu}$ 121 16.32 17.15 +5.08% 3.78 0.06

Bán kính tinh thể tối ưu tính theo trọng số: $\bar{R} = 3.76 \pm 0.03\text{ cm}$ (giảm so với danh định $R = 3.81\text{ cm}$).

Thí nghiệm 3 & Đánh giá mô hình hoàn chỉnh

Sau khi cập nhật cả 3 thông số tối ưu ($\rho = 2.02\text{ g/cm}^3, R = 3.76\text{ cm}, L = 7.55\text{ cm}$):

So sánh độ lệch hiệu suất FEPE (%) giữa Mô hình gốc và Mô hình tối ưu:

Năng lượng (keV) |  0%        5%        10%       15%       20%
-----------------+---------------------------------------------
31 keV (Gốc)     |  [===================================] 19.57%
31 keV (Tối ưu)  |  [==] 1.28%
-----------------+---------------------------------------------
81 keV (Gốc)     |  [=======] 3.99%
81 keV (Tối ưu)  |  [===] 1.65%
-----------------+---------------------------------------------
662 keV (Gốc)    |  [======] 3.42%
662 keV (Tối ưu) |  [==] 1.15%
-----------------+---------------------------------------------
1408 keV (Gốc)   |  [=====] 2.89%
1408 keV (Tối ưu)|  [=] 0.85%
Mức năng lượng $E_\gamma$ (keV) $\varepsilon_{Exp} \times 10^{-4}$ $\varepsilon_{Sim}^{Opt} \times 10^{-4}$ Sai số tương đối ban đầu Sai số tương đối mô hình tối ưu
31 ($^{133}\text{Ba}$) 16.40 16.61 +19.57% +1.28%
32 ($^{137}\text{Cs}$) 16.64 16.89 +19.47% +1.50%
59 ($^{241}\text{Am}$) 19.48 19.82 +3.95% +1.74%
81 ($^{133}\text{Ba}$) 17.52 17.81 +3.99% +1.65%
121 ($^{152}\text{Eu}$) 16.32 16.67 +5.08% +2.14%
662 ($^{137}\text{Cs}$) 5.82 5.89 +3.42% +1.15%
1173 ($^{60}\text{Co}$) 5.71 5.65 +2.98% -1.05%
1332 ($^{60}\text{Co}$) 5.12 5.07 +2.73% -0.97%
1408 ($^{152}\text{Eu}$) 4.88 4.92 +2.89% +0.85%

Đổi mới và đóng góp

  1. Phương pháp tách biến độc lập (Parameter Decoupling Strategy): Khác với các nghiên cứu trước đây (như nghiên cứu của Hoàng Đức Tâm và cộng sự [5] chỉ hiệu chỉnh độ dày lớp phản xạ $d_{Al_2O_3}$ nhưng vẫn giữ nguyên mật độ lý thuyết), nghiên cứu này tách rời hoàn toàn ảnh hưởng của mật độ vật liệu hấp thụ bằng kỹ thuật chuẩn trực năng lượng thấp, triệt tiêu sự phụ thuộc vào kích thước dài và bán kính.
  2. Định lượng lại thông số thực tế của đầu dò: Chứng minh thực nghiệm rằng mật độ lớp phản xạ $Al_2O_3$ trong đầu dò thực tế đạt $2.02\text{ g/cm}^3$ (thay vì $0.55\text{ g/cm}^3$ dạng xốp lỏng lẻo như công bố của NSX), và bán kính hiệu dụng của tinh thể bị suy giảm $0.05\text{ cm}$ do lớp suy giảm quang học ở bề mặt tiếp xúc.
  3. Độ chính xác vượt trội: Hạ thấp sai số tương đối từ mức xấp xỉ 20% xuống dưới 2.16% trên toàn bộ dải năng lượng từ 31 keV đến 1408 keV.

Ứng dụng thực tế và triển khai

+-------------------------------------------------------------------------+
|                  HỆ THỐNG MÔ PHỎNG VÀ HIỆU CHUẨN ĐẦU DÒ                 |
+-------------------------------------------------------------------------+
       |                                              |
       v                                              v
+-----------------------------+                +--------------------------+
|  QUAN TRẮC MÔI TRƯỜNG       |                |  Y HỌC & ĐỒNG VỊ PHÓNG XẠ|
|  - Đo hoạt độ mẫu đất, nước |                |  - Đo liều I-131, Tc-99m |
|  - Mẫu thể tích Marinelli   |                |  - Kiểm chuẩn buồng đo   |
|  - Tính tự hấp thụ chính xác|                |  - Giảm liều xạ nhân viên|
+-----------------------------+                +--------------------------+
  • Xây dựng đường cong hiệu suất cho các hình học mẫu phức tạp: Sử dụng mô hình MCNP5 đã tối ưu để tính toán hiệu suất cho mẫu hình trụ, mẫu hộp Marinelli chứa đất, bùn, nước mà không cần chế tạo mẫu chuẩn phóng xạ đắt tiền.
  • Tối ưu hóa an toàn bức xạ: Giảm 85% thời gian thao tác trực tiếp của nghiên cứu viên với nguồn phóng xạ mạnh trong phòng thí nghiệm.
  • Cấu hình phần cứng yêu cầu:
    • Hệ thống máy tính: CPU 4 nhân trở lên (Intel Core i5/i7/Xeon), RAM 8GB.
    • Phần mềm: MCNP5 v1.40/v1.60, MCNPX, trình biên dịch ANSI Fortran 90, công cụ trực quan hóa Vised X_22S.

Hạn chế và hướng phát triển

  • Hạn chế kỹ thuật: Chưa khảo sát chi tiết hiện tượng trôi đỉnh do nhiệt độ làm thay đổi độ lợi quang điện của ống nhân quang điện (PMT); sai số nội suy phụ thuộc vào độ chính xác của số liệu hệ số suy giảm khối trong thư viện photon NIST.
  • Hướng phát triển tiếp theo:
    1. Ứng dụng thuật toán tối ưu hóa thông minh (Genetic Algorithm - GA hoặc Particle Swarm Optimization - PSO) tự động hóa quá trình tìm kiếm không gian tham số 3D ($\rho, R, L$).
    2. Mở rộng mô hình cho các loại đầu dò bán dẫn HPGe (High-Purity Germanium) và đầu dò nhấp nháy thế hệ mới $LaBr_3(Ce)$.

Đối tượng hưởng lợi

  • Sinh viên & Học viên cao học ngành Vật lý hạt nhân/Kỹ thuật hạt nhân: Tài liệu tham khảo chuẩn mực về phương pháp xây dựng file input MCNP5 và xử lý dữ liệu phổ gamma.
  • Kỹ sư vận hành hệ đo bức xạ: Công cụ đắc lực để hiệu chuẩn thiết bị đo mà không phụ thuộc hoàn toàn vào nguồn chuẩn ngoại nhập.
  • Phòng thí nghiệm phân tích hạt nhân: Tiết kiệm hàng ngàn USD chi phí kiểm định và mua sắm nguồn chuẩn phóng xạ hàng năm.

Câu hỏi thường gặp

1. Tại sao mật độ lớp phản xạ $Al_2O_3$ tối ưu ($2.02\text{ g/cm}^3$) lại lớn hơn nhiều so với thông số NSX ($0.55\text{ g/cm}^3$)?

Thông số $0.55\text{ g/cm}^3$ từ NSX thường là mật độ biểu kiến khi đóng gói bột xốp tự do. Trong quá trình ép chân không và lắp ráp kín tinh thể, bột oxit nhôm bị nén chặt, làm mật độ thực tế tăng lên tiệm cận giá trị tinh thể đặc ($3.95\text{ g/cm}^3$).

2. Tại sao nguồn $^{133}\text{Ba}$ (31 keV) được chọn để xác định mật độ lớp phản xạ?

Ở năng lượng 31 keV, tiết diện phản ứng quang điện trên nhôm và iod rất lớn, quãng chạy của photon cực ngắn (chưa tới 1 mm trong tinh thể). Do đó, photon bị tán xạ và hấp thụ mạnh nhất tại lớp vỏ/lớp phản xạ trước khi chạm tới thể tích nhạy của tinh thể, biến 31 keV thành điểm khảo sát nhạy nhất với biến số mật độ $\rho_{Al_2O_3}$.

3. Có thể bỏ qua việc chuẩn trực trong phép đo xác định mật độ không?

Không. Nếu không chuẩn trực, chùm tia sẽ quét qua các góc tới khác nhau, làm thay đổi góc khối hình học và chịu ảnh hưởng đan xen bởi bán kính mặt cong tinh thể, khiến việc tách biến độc lập không thể thực hiện chính xác.

4. MCNP5 xử lý bài toán vận chuyển photon trong tinh thể NaI(Tl) bằng Tally nào?

MCNP5 sử dụng Tally F8 (Pulse Height Tally) cho photon (F8:P). Tally này ghi nhận toàn bộ năng lượng bị mất đi bởi các hạt thứ cấp trong thể tích cell tinh thể trên từng phân rã, tái hiện chính xác đáp ứng xung của detector.

5. Chi phí triển khai phương pháp tối ưu hóa này là bao nhiêu?

Chi phí thực tế gần như bằng không nếu phòng thí nghiệm đã có sẵn máy tính cá nhân và bộ phổ kế cơ bản, giúp giảm chi phí mua bộ nguồn chuẩn đa năng lượng từ 3.000 – 5.000 USD.

Kết luận

Đồ án đã giải quyết triệt để bài toán bất định về thông số kỹ thuật của đầu dò nhấp nháy NaI(Tl) trong mô phỏng bức xạ. Bằng quy trình tách rời 3 giai đoạn kết hợp thực nghiệm và mô phỏng MCNP5, nghiên cứu đã tìm ra bộ thông số tối ưu: mật độ lớp phản xạ $\rho_{Al_2O_3} = 2.02 \pm 0.32\text{ g/cm}^3$, bán kính tinh thể $R = 3.76 \pm 0.03\text{ cm}$, và chiều dài tinh thể $L = 7.55\text{ cm}$. Kết quả này đưa độ lệch hiệu suất giữa mô phỏng và thực nghiệm từ mức xấp xỉ 20% về dưới 2.16% trên toàn bộ dải năng lượng 31 – 1408 keV, cung cấp một giải pháp tin cậy, tiết kiệm và an toàn cho các ứng dụng đo đạc bức xạ hạt nhân hiện đại.