Tổng quan nghiên cứu
Trong công tác bảo đảm an toàn bức xạ và an toàn hạt nhân, việc kiểm soát liều lượng bức xạ neutron đóng vai trò đặc biệt quan trọng do đặc tính ion hóa gián tiếp và hiệu ứng sinh học phức tạp của hạt neutron. Theo quy định tại Điều 24 Luật Năng lượng Nguyên tử Việt Nam năm 2008 và khuyến cáo của Cơ quan Năng lượng Nguyên tử Quốc tế (IAEA), 100% các thiết bị đo liều bức xạ cầm tay và liều kế cá nhân phải được hiệu chuẩn định kỳ từ 12 đến 14 tháng một lần. Tuy nhiên, tính đến năm 2014, Phòng Chuẩn cấp 2 về đo liều bức xạ ion hóa của Việt Nam (SSDL-VN) đặt tại Viện Khoa học và Kỹ thuật Hạt nhân (VKHKTHN) thuộc Viện Năng lượng Nguyên tử Việt Nam mới chỉ có năng lực chuẩn liều cho bức xạ photon (tia gamma và tia X), hoàn toàn thiếu vắng hệ thống chuẩn liều bức xạ neutron phục vụ nhu cầu thực tế.
Nghiên cứu của tác giả Lê Ngọc Thiệm thực hiện năm 2014 tại Hà Nội tập trung giải quyết bài toán cốt lõi: Xây dựng cơ sở khoa học và thiết lập trường chuẩn liều neutron sử dụng nguồn đồng vị phân hạch tự phát Californium-252 (252Cf). Mục tiêu cụ thể của luận văn là tính toán phổ thông lượng neutron, xác định tương đương liều môi trường $H^*(10)$, phân tách tường minh thành phần chùm neutron trực tiếp và tán xạ phòng, đồng thời thẩm định phương pháp trên thiết bị đo liều thực nghiệm Aloka TPS-451C. Kết quả nghiên cứu mang ý nghĩa tiên phong, đặt nền móng kỹ thuật vững chắc để mở rộng năng lực SSDL-VN đạt chuẩn quốc tế ISO 8529, kiểm soát sai số đo liều dưới 5%, trực tiếp phục vụ an toàn cho hơn 50 cơ sở y tế, công nghiệp và nghiên cứu hạt nhân trên cả nước.
Cơ sở lý thuyết và phương pháp nghiên cứu
Khung lý thuyết áp dụng
Nghiên cứu xây dựng trên nền tảng lý thuyết tương tác neutron với vật chất và hệ thống đo lường học bức xạ của Ủy ban Đo lường và Đơn vị Bức xạ Quốc tế (ICRU) cùng Ủy ban Quốc tế về An toàn Bức xạ (ICRP):
- Lý thuyết tương tác neutron vi mô và vĩ mô: Quá trình neutron truyền qua môi trường bao gồm tán xạ đàn hồi $(n,n)$, tán xạ không đàn hồi $(n,n')$, hấp thụ phóng xạ $(n,\gamma)$, và phân hạch hạt nhân $(n,f)$. Mối quan hệ giữa tiết diện vi mô $\sigma$ (đo bằng barn, $1\text{ b} = 10^{-24}\text{ cm}^2$) và tiết diện vĩ mô $\Sigma_t = N\sigma_t$ (đơn vị $\text{cm}^{-1}$) quyết định quãng đường chạy tự do trung bình $\lambda = 1/\Sigma_t$. Khả năng làm chậm neutron nhanh ($2\text{ MeV}$) xuống vùng nhiệt ($0.025\text{ eV}$) được đánh giá thông qua độ suy giảm năng lượng logarit trung bình $\xi$ và công suất làm chậm $\xi\Sigma_s$, trong đó Polyethylene giàu hydro (H) thể hiện hiệu suất cao với 27 lần va chạm trung bình.
- Hệ thống đại lượng liều lượng học chuẩn: Trường chuẩn được xác lập thông qua thông lượng neutron $\Phi$ ($\text{cm}^{-2}\cdot\text{s}^{-1}$), liều hấp thụ $D$ (Gy, J/kg), và tương đương liều môi trường $H^*(10) = D \cdot Q$ (đơn vị Sievert), với $Q$ là hệ số phẩm chất bức xạ phụ thuộc năng lượng hạt từ $0.025\text{ eV}$ đến $10\text{ MeV}$.
- Hệ số chuyển đổi thông lượng sang liều $h^*(10)$: Tỷ số giữa suất tương đương liều $H^*(10)$ và thông lượng neutron $\Phi$ tại điểm khảo sát, biểu diễn đặc trưng đáp ứng liều của trường bức xạ tự do và định hướng.
Phương pháp nghiên cứu
- Nguồn dữ liệu và cấu hình mẫu: Nghiên cứu khai thác thư viện số liệu hạt nhân chuẩn quốc tế kết hợp cấu hình hình học thực tế của phòng chuẩn SSDL-VN. Nguồn chuẩn sử dụng là đồng vị 252Cf phân hạch tự phát với phổ phân bố năng lượng Watt đặc trưng có năng lượng trung bình $\bar{E} = 2.13\text{ MeV}$.
- Phương pháp mô phỏng Monte Carlo: Phương pháp thống kê ngẫu nhiên mô phỏng quá trình truyền hạt của hàng triệu lịch sử neutron ($10^7$ hạt) qua các môi trường không khí, tấm chắn, sàn và tường bê tông nhằm giảm sai số thống kê xuống dưới 1.5%. Phương pháp này được chọn vì có khả năng giải quyết chính xác phương trình truyền xạ đa chiều trong không gian hình học phức tạp.
- Kỹ thuật tấm chắn hình nón cụt (Shadow Cone Technique): Phương pháp thực nghiệm - mô phỏng chuẩn tắc theo tiêu chuẩn ISO 8529-2 nhằm phân tách thông lượng tổng cộng ($n\text{-tot}$) thành thông lượng trực tiếp ($n\text{-dir}$) và thông lượng tán xạ phòng ($n\text{-sct}$). Tấm chắn thiết kế dạng nón cụt dài 30 cm, cấu tạo từ Polyethylene pha 5% Boron (PEB mật độ $1.0\text{ g/cm}^3$) kết hợp lớp chì 0.5 cm Pb phía trước để triệt tiêu chùm tia neutron trực tiếp tới đầu đo.
- Quy trình đo đạc thực nghiệm: Khảo sát tại 3 khoảng cách tiêu chuẩn: 75 cm, 125 cm và 150 cm tính từ tâm nguồn 252Cf đến điểm đo, sử dụng máy đo liều neutron chuyên dụng Aloka TPS-451C trong mốc thời gian nghiên cứu năm 2014.
Kết quả nghiên cứu và thảo luận
Những phát hiện chính
- Phân bố phổ thông lượng và sự phụ thuộc khoảng cách: Tại vị trí 75 cm tính từ nguồn 252Cf không có tấm chắn, thông lượng neutron trực tiếp $n\text{-dir}$ chiếm ưu thế tuyệt đối với 81.4% tổng thông lượng, trong khi thành phần tán xạ phòng $n\text{-sct}$ chỉ chiếm 18.6%. Khi tăng khoảng cách lên 125 cm và 150 cm, thông lượng trực tiếp giảm tỷ lệ nghịch với bình phương khoảng cách ($1/r^2$), khiến tỷ lệ đóng góp của tán xạ phòng tăng vọt lên tương ứng là 34.2% và 42.8% tổng thông lượng ghi nhận.
- Hiệu năng che chắn của cấu hình tấm chắn Shadow Cone: Cấu hình tấm chắn gồm 30 cm PEB kết hợp 0.5 cm Pb đạt hiệu suất suy giảm thông lượng trực tiếp trên 99.2%, tỷ lệ rò rỉ neutron trực tiếp xuyên qua tấm chắn ghi nhận được ở mức cực thấp dưới 0.8% trên toàn dải phổ. Điều này đảm bảo khi đặt tấm chắn, đầu đo chỉ ghi nhận thành phần tán xạ phòng thuần túy với độ tin cậy thống kê cao.
- Xác lập hệ số hiệu chuẩn cho máy đo liều Aloka TPS-451C: Sai số tương đối nội tại giữa giá trị mô phỏng tương đương liều trực tiếp $DE_{\text{dir}}$ và số liệu đo đạc thực tế từ thiết bị Aloka TPS-451C chỉ dao động trong khoảng từ 3.2% đến 4.8% tại khoảng cách 75 cm và 125 cm, hoàn toàn nằm trong giới hạn sai số cho phép ($\pm 15%$) của các chuẩn đo lường cấp 2 quốc tế.
Thảo luận kết quả
Thành phần tán xạ từ tường, trần và sàn phòng chuẩn đóng vai trò chi phối khi tăng khoảng cách nguồn - đầu đo. Hiện tượng này giải thích do neutron nhanh phát ra từ nguồn 252Cf va chạm đàn hồi với các hạt nhân Hydro, Oxy, Silic trong bê tông, bị nhiệt hóa và tán xạ ngược lại vùng đo. Dữ liệu phổ phân bố năng lượng neutron và đóng góp phần trăm được minh họa tường minh qua các đồ thị phổ thông lượng vi phân từ $10^{-8}\text{ MeV}$ đến $10\text{ MeV}$ và bảng đối chiếu tỷ lệ thành phần:
| Khoảng cách (cm) | Thông lượng trực tiếp $n\text{-dir}$ (%) | Thông lượng tán xạ $n\text{-sct}$ (%) | Suất tương đương liều trực tiếp (pSv/n) | Sai số thống kê (%) |
|---|---|---|---|---|
| 75 | 81.4 | 18.6 | $3.85 \times 10^{-1}$ | $\pm 0.6$ |
| 125 | 65.8 | 34.2 | $1.39 \times 10^{-1}$ | $\pm 0.9$ |
| 150 | 57.2 | 42.8 | $0.96 \times 10^{-1}$ | $\pm 1.2$ |
So sánh với các trường chuẩn quy chuẩn trong báo cáo kỹ thuật IAEA TRS-403, phương pháp phân tách thành phần bức xạ bằng kỹ thuật tấm chắn nón cụt tại SSDL-VN loại bỏ triệt để hiện tượng đánh giá vượt mức liều thực tế, khẳng định độ chính xác tương đương các phòng thí nghiệm hạt nhân tiên tiến trong khu vực.
Đề xuất và khuyến nghị
- Ban hành quy trình hiệu chuẩn chuẩn quốc gia (SOP): Cục An toàn Bức xạ và Hạt nhân phối hợp cùng VKHKTHN hoàn thiện quy trình kiểm định thiết bị đo liều neutron định kỳ 12 tháng/lần, kiểm soát sai số tổng thể dưới 5%, triển khai trong vòng 6 tháng tới.
- Nâng cấp vật liệu che chắn phòng đo SSDL-VN: Thi công lắp đặt thêm lớp vật liệu hấp thụ neutron nhiệt chứa Boron hoặc Cadmium dày 2 mm phủ kín bề mặt tường bê tông nhằm giảm tỷ trọng tán xạ phòng xuống dưới 15% tại khoảng cách 150 cm trước quý IV.
- Xây dựng chương trình phần mềm tự động hiệu chỉnh phân rã nguồn: Do nguồn 252Cf có chu kỳ bán rã ngắn ($T_{1/2} = 2.645\text{ năm}$), tỷ lệ phân rã hoạt độ khoảng 23.1%/năm, VKHKTHN cần tích hợp thuật toán tự động bù trừ hoạt độ và suất thông lượng neutron theo thời gian thực vào phần mềm điều khiển hiệu chuẩn hàng tháng.
- Thực hiện so sánh liên phòng cấp quốc tế: Đăng ký tham gia chương trình so sánh liên phòng chuẩn liều neutron với Cơ quan Cân đo Quốc tế (BIPM) và mạng lưới SSDL của IAEA/WHO định kỳ 2 năm một lần nhằm duy trì tính liên kết chuẩn đo lường hạt nhân.
Đối tượng nên tham khảo luận văn
- Cán bộ phụ trách an toàn bức xạ tại các cơ sở hạt nhân: Kỹ sư vận hành lò phản ứng nghiên cứu Đà Lạt, các cơ sở sản xuất đồng vị phóng xạ và máy gia tốc y tế cần tài liệu quy chuẩn để đánh giá suất liều neutron môi trường xung quanh khu vực làm việc.
- Chuyên viên phòng thí nghiệm kiểm định và hiệu chuẩn đo lường: Kỹ thuật viên tại các trung tâm tiêu chuẩn đo lường chất lượng khai thác quy trình thực hành kỹ thuật Shadow Cone để hiệu chuẩn chính xác máy đo liều cầm tay Aloka, Ludlum, Thermo Fisher.
- Giảng viên, học viên cao học và nghiên cứu sinh ngành Vật lý hạt nhân: Nhóm nghiên cứu học thuật sử dụng luận văn làm tài liệu chuyên khảo chuyên sâu về phương pháp tính toán tương tác neutron và kỹ thuật mô phỏng Monte Carlo trong vật lý che chắn.
- Kỹ sư thiết kế công trình che chắn bức xạ: Các chuyên gia thiết kế hầm máy gia tốc xạ trị proton/neutron và kho lưu giữ chất thải phóng xạ ứng dụng các công thức tính toán độ dày làm chậm và tiết diện truyền qua của vật liệu PEB và chì.
Câu hỏi thường gặp
Tại sao nguồn 252Cf được lựa chọn để xây dựng trường chuẩn liều neutron?
Nguồn 252Cf phát ra neutron thông qua cơ chế phân hạch tự phát với phổ năng lượng liên tục từ nhiệt đến trên $10\text{ MeV}$, năng lượng trung bình $2.13\text{ MeV}$, tương đồng với phổ neutron trong lò phản ứng và máy gia tốc. Suất phát neutron riêng của 252Cf rất cao ($2.314 \times 10^{12}\text{ n/s}\cdot\text{g}$), tạo ra trường bức xạ ổn định phục vụ hiệu chuẩn.
Kỹ thuật tấm chắn hình nón cụt (Shadow Cone) giải quyết vấn đề gì?
Trong phòng chuẩn thực tế, tường và sàn luôn gây ra bức xạ tán xạ ngược. Tấm chắn PEB 30 cm đặt giữa nguồn và đầu dò giúp che hoàn toàn chùm neutron trực tiếp, chỉ cho phép bức xạ tán xạ tới đầu đo. Bằng phép trừ thông lượng tổng và tán xạ, thông lượng trực tiếp thuần túy được xác định chính xác theo chuẩn ISO 8529.
Đại lượng tương đương liều môi trường $H^*(10)$ có ý nghĩa gì trong kiểm soát an toàn?
$H^*(10)$ là đại lượng đo lường quy ước được định nghĩa là tương đương liều tại độ sâu 10 mm trên bán kính quả cầu mô ICRU ($30\text{ cm}$ đường kính), ngược hướng chùm tia. Đại lượng này đại diện cho khả năng đâm xuyên mạnh của neutron vào cơ thể người, dùng làm tiêu chuẩn đánh giá an toàn khu vực làm việc.
Mức độ sai lệch giữa mô phỏng tính toán và đo đạc thực nghiệm là bao nhiêu?
Độ sai lệch tương đối giữa mô phỏng Monte Carlo và thực nghiệm bằng máy đo Aloka TPS-451C chỉ từ 3.2% đến 4.8% tại khoảng cách 75 cm và 125 cm. Kết quả này chứng minh mô hình tính toán hình học phòng và cấu tạo đầu đo đạt độ tin cậy vượt trội so với yêu cầu sai số $\pm 15%$ của phòng chuẩn cấp 2.
Chu kỳ hiệu chuẩn định kỳ bắt buộc của thiết bị đo bức xạ neutron là bao lâu?
Theo khuyến cáo của IAEA và Điều 24 Luật Năng lượng Nguyên tử, tất cả các thiết bị đo liều neutron cầm tay phải được hiệu chuẩn trước khi sử dụng lần đầu và tái hiệu chuẩn định kỳ sau mỗi 12 đến 14 tháng hoạt động nhằm đảm bảo độ chính xác của hệ số chuẩn $CF$.
Kết luận
- Luận văn đã xây dựng thành công cơ sở lý thuyết và mô hình tính toán trường chuẩn liều bức xạ neutron sử dụng nguồn đồng vị 252Cf đầu tiên tại Việt Nam.
- Thiết lập hoàn chỉnh phương pháp thực nghiệm tấm chắn hình nón cụt Shadow Cone (30 cm PEB + 0.5 cm Pb), bóc tách chính xác thành phần tán xạ phòng chiếm từ 18.6% đến 42.8% tổng thông lượng.
- Kiểm chứng thực nghiệm thành công trên thiết bị đo liều Aloka TPS-451C với độ sai lệch kiểm soát chặt chẽ dưới mức 5%.
- Đặt nền móng kỹ thuật trực tiếp để đưa Phòng chuẩn liều SSDL-VN tham gia toàn diện vào mạng lưới đo lường chuẩn bức xạ neutron quốc tế IAEA/WHO.
- Trong giai đoạn tiếp theo, quy trình cần được chuẩn hóa thành Tiêu chuẩn Quốc gia (TCVN) và ứng dụng kiểm định đại trà cho toàn bộ thiết bị đo bức xạ hạt nhân trên toàn quốc. Các đơn vị nghiên cứu và doanh nghiệp ứng dụng bức xạ có nhu cầu hiệu chuẩn có thể liên hệ trực tiếp Viện Khoa học và Kỹ thuật Hạt nhân để phối hợp triển khai.