Tổng quan nghiên cứu

Kiểm soát và định lượng vật liệu hạt nhân đặc biệt như Plutonium (ký hiệu hóa học Pu, số nguyên tử 94) là nhiệm vụ cốt lõi trong công tác thanh sát an ninh hạt nhân toàn cầu theo tiêu chuẩn của Cơ quan Năng lượng Nguyên tử Quốc tế IAEA. Trong thực tế, chỉ cần 6,2 kg Plutonium đã có thể tạo ra sức công phá tương đương 20.000 tấn thuốc nổ TNT, đồng thời khối lượng tới hạn của đồng vị 239Pu chỉ ở mức 10 kg khi có lớp phản xạ neutron. Tại nhiều quốc gia, số lượng nguồn phóng xạ kín PuBe tồn tại rất lớn trong công nghiệp và nghiên cứu (điển hình như Hungary từng quản lý tới 152 nguồn PuBe với tỷ lệ hỗn hợp chuẩn Pu/Be = 1/13).

Vấn đề kỹ thuật nan giải phát sinh khi phần lớn các nguồn phóng xạ này đã bị thất lạc hồ sơ kỹ thuật, không rõ kích thước hình học bên trong, hoặc có độ cháy cao dẫn đến hàm lượng đồng vị 239Pu giảm xuống dưới 75%. Trong các trường hợp này, những gói phần mềm thương mại phổ biến như MGA++ hoàn toàn mất khả năng tự động hội tụ và tính toán chính xác.

Mục tiêu cụ thể của nghiên cứu là xây dựng phương pháp phổ kế gamma độ phân giải cao kết hợp quy trình xử lý chuyên sâu trên nền tảng phần mềm Gamma Vision để xác định đồng thời tỷ lệ thành phần các đồng vị (238Pu, 239Pu, 240Pu, 241Pu, 242Pu, 241Am) và khối lượng tổng của Plutonium mà không cần phá hủy mẫu. Nghiên cứu được thực nghiệm và kiểm chứng trên các tập dữ liệu nguồn phóng xạ đo vào năm 2002 và năm 2009, giúp hạ sai số xác định khối lượng xuống dưới mức 6% đối với nguồn đã biết kích thước và dưới 10% đến 15% đối với các nguồn mất hoàn toàn thông số vỏ bọ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 vận dụng lý thuyết tương tác bức xạ gamma với vật chất, bao gồm hiệu ứng quang điện, tán xạ Compton và cơ chế tự hấp thụ bức xạ trong nguồn khối dày (self-absorption). Cơ sở toán học nền tảng dựa trên phương trình biểu diễn mối liên hệ giữa tốc độ đếm đỉnh quang điện với số lượng hạt nhân phóng xạ, hằng số phân rã, hệ số phân nhánh bức xạ, góc khối và hệ số truyền qua của môi trường.

Để giải quyết bài toán suy giảm cường độ gamma qua các lớp vật liệu, mô hình tính toán tích hợp thừa số hiệu chỉnh tự hấp thụ trong mẫu $F_{1E}$ thông qua nghiệm giải tích của hàm Bessel hiệu chỉnh bậc 1 ($I_1$) kết hợp hàm Struve hiệu chỉnh bậc 1 ($L_1$), cùng với thừa số hấp thụ của lớp vỏ thép bảo vệ $F_{2E}$. Các khái niệm định lượng cốt lõi được sử dụng xuyên suốt gồm: hoạt độ riêng (đơn vị Bq/g), độ cháy nhiên liệu hạt nhân (Burnup), đường cong hiệu suất ghi tương đối (Relative Efficiency Curve), và tỷ số hình học $R = H/D$ giữa chiều cao và đường kính khối phóng xạ.

Phương pháp nghiên cứu

Nghiên cứu sử dụng cỡ mẫu gồm toàn bộ phổ gamma của các nguồn chuẩn PuBe và các nguồn thực tế được thu thập qua hai đợt đo đạc kiểm chuẩn vào năm 2002 và năm 2009. Phương pháp chọn mẫu là chọn mẫu có chủ đích bao phủ nhiều dải hàm lượng 239Pu từ thấp (dưới 75%) đến cao (trên 95%) nhằm đánh giá toàn diện độ tin cậy của thuật toán.

Hệ thiết bị thực nghiệm sử dụng detector bán dẫn siêu tinh khiết HPGe kiểu Planar model GL2020R (hãng Canberra) có diện tích tinh thể 2000 mm2, đạt độ phân giải năng lượng vượt trội 490 eV tại đỉnh 122 keV của 57Co. Tín hiệu được thu nhận qua bộ phân tích đa kênh MCA cấu hình 8196 kênh với độ rộng 0,075 keV/kênh. Khoảng cách đo được thiết lập linh hoạt từ 20 cm đến 200 cm tùy theo cường độ nguồn nhằm bảo vệ tinh thể bán dẫn khỏi hư hại bởi thông lượng neutron nhanh.

Lý do lựa chọn phương pháp phân tích phổ gamma bằng phần mềm Gamma Vision 32 kết hợp hiệu chỉnh bán tự công là vì công cụ này cho phép can thiệp trực tiếp vào từng vùng đỉnh năng lượng (ROI), xử lý tách chập các cụm đỉnh phức tạp (như 94-104 keV, 125 keV, 148 keV, 160 keV, 208 keV, 332 keV và 375 keV), từ đó vượt qua giới hạn thuật toán tự động của chương trình MGA++ truyền thống.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Thứ nhất, đối với phép xác định hàm lượng đồng vị 239Pu, phương pháp Gamma Vision sử dụng đỉnh 104,24 keV và đỉnh 160,0 keV cho kết quả tương thích rất cao với phần mềm chuẩn MGA++, với độ lệch tương đối được kiểm soát chặt chẽ dưới 10% trên toàn bộ các dải hàm lượng thử nghiệm.

Thứ hai, nghiên cứu đã giải quyết thành công bài toán tính toán hàm lượng đồng vị 240Pu ở các nguồn có độ cháy cao (nơi hàm lượng 239Pu nhỏ hơn 85%), đạt độ chênh lệch dưới 40% so với MGA++. Đây là cải tiến then chốt vì phần mềm tự động thường bị mất hội tụ tại vùng này do tỷ phần đóng góp của 240Pu trong đỉnh tổng 160 keV bị suy giảm mạnh từ 70% (ở mẫu chứa 6% 240Pu) xuống chỉ còn 25% (ở mẫu chứa 20% 240Pu).

Thứ ba, việc phát hiện quy luật sai số hệ thống của thừa số tự hấp thụ $F_{1E}$ tăng tỷ lệ thuận theo giá trị độ dày quang học $\mu_{lE} D$ (với độ lệch khối lượng ban đầu lên tới 12%) đã dẫn đến việc xây dựng công thức hiệu chỉnh thực nghiệm $f(\mu_{lE} D) = 1,095 \log_2(\mu_{lE} D)$. Nhờ thừa số hiệu chỉnh này, sai số xác định tổng khối lượng Plutonium thông qua hai đỉnh 375,04 keV và 413,71 keV đã giảm mạnh từ 12% xuống dưới 6% khi so sánh với phương pháp chuẩn nhiệt lượng kế $M_{cal}$.

Thứ tư, đối với các nguồn PuBe mất hoàn toàn nhãn mác kỹ thuật, giải thuật tối ưu hóa tỷ số hình học hiệu dụng $R_{eff}$ dựa trên đỉnh 129,29 keV kết hợp phương pháp năng lượng vô hạn đã ước lượng thành công khối lượng Plutonium với độ dày vỏ thép quy ước 3,5 mm (hoặc 7 mm cho vỏ kép), duy trì sai số tổng thể trong khoảng 10% đến 15%.

Thảo luận kết quả

Khi biểu diễn mối quan hệ giữa tỷ số $M_\gamma / M_{cal}$ theo thông số $\mu_{lE} D$ trên đồ thị hàm số, đường cong thực nghiệm phản ánh rõ bản chất phi tuyến của quá trình tự hấp thụ gamma trong nguồn khối có mật độ vật chất lớn ($\rho \approx 3,7\text{ g/cm}^3$). Bảng dữ liệu phổ kế tại các vùng năng lượng cho thấy tại dải 100 keV có tới 14 vạch phổ gamma và tia X đặc trưng của Uranium, Neptunium và Plutonium chồng lấn dày đặc.

Khác với phương pháp đo nhiệt lượng kế vốn đòi hỏi thiết bị phức tạp, thời gian chờ cân bằng nhiệt kéo dài hàng giờ và hoàn toàn không thể xác định tỷ lệ từng đồng vị, phương pháp phổ kế gamma cải tiến trong luận văn vừa cung cấp phổ đồng vị chi tiết vừa tính toán chính xác tổng khối lượng hạt nhân trong thời gian đo từ 30 đến 60 phút.

Đề xuất và khuyến nghị

  1. Nâng cấp phần mềm phân tích phổ gamma hạt nhân tại các cơ sở nghiên cứu và kiểm định bức xạ: Viện Năng lượng Nguyên tử Việt Nam cùng các phòng thí nghiệm chuyên ngành cần tích hợp thuật toán hiệu chỉnh tự hấp thụ phi tuyến và giải chập cụm đỉnh 160 keV trong vòng 12 tháng tới, nhằm mục tiêu đưa sai số phân tích đồng vị của các nguồn phóng xạ độ cháy cao về dưới mức 5%.

  2. Chuẩn hóa quy trình thanh sát và kiểm kê nguồn PuBe mất nhãn mác lưu kho: Cục An toàn Bức xạ và Hạt nhân cần ban hành hướng dẫn kỹ thuật áp dụng mô hình ước lượng tỷ số hình học hiệu dụng $R_{eff}$ qua đỉnh 129 keV vào quý 3 năm tới, bảo đảm 100% các nguồn phóng xạ không rõ nguồn gốc trong diện quản lý quốc gia được định lượng chính xác.

  3. Áp dụng hệ phin lọc hấp thụ gamma năng lượng thấp (tấm Cadmium dày 0,15 mm hoặc tấm đồng/chì mỏng) trong các phép đo phổ: Đội ngũ kỹ thuật viên đo lường bức xạ cần thực hiện ngay trong vòng 6 tháng tới để triệt tiêu hiện tượng chồng phổ do đỉnh 59,5 keV của 241Am gây ra, tăng tỷ số tín hiệu trên phông (S/N) tại vùng phổ 100-208 keV lên hơn 25%.

  4. Xây dựng cơ sở dữ liệu phổ chuẩn đồng vị Actinide và siêu Uranium: Bộ Khoa học và Công nghệ chủ trì phối hợp với các trường đại học triển khai dự án nghiên cứu 24 tháng, hoàn thiện thư viện chuẩn trên 500 mẫu phổ phục vụ đào tạo thực hành và hỗ trợ công tác thanh sát an ninh hạt nhân theo hiệp định IAEA.

Đối tượng nên tham khảo luận văn

  1. Thanh sát viên an toàn bức xạ và cán bộ quản lý vật liệu hạt nhân: Nắm vững phương pháp phân tích không phá hủy (NDA) để giám định, lập hồ sơ kiểm kê và phát hiện các hành vi vận chuyển trái phép nguồn phóng xạ chứa Plutonium.

  2. Kỹ sư vận hành lò phản ứng và chuyên gia tái chế nhiên liệu hạt nhân: Ứng dụng quy trình tách chập phổ để kiểm tra độ cháy, tỷ lệ làm giàu và hàm lượng đồng vị trong thanh nhiên liệu ôxít hỗn hợp MOX hoặc chất thải phóng xạ siêu Uranium.

  3. Học viên cao học và nghiên cứu sinh ngành Vật lý hạt nhân, Kỹ thuật hạt nhân: Sử dụng tài liệu như một cẩm nang mẫu mực về kỹ thuật phổ kế bán dẫn HPGe, phương pháp hiệu chuẩn hiệu suất ghi tương đối và kỹ năng xử lý phổ phức hợp trên phần mềm Gamma Vision và MGA++.

  4. Chuyên gia đo lường và ứng phó sự cố bức xạ: Tiếp cận giải thuật ước lượng nhanh khối lượng hạt nhân của các nguồn phóng xạ bị biến dạng hình học hoặc mất thông tin kỹ thuật phục vụ đánh giá rủi ro an ninh phóng xạ khẩn cấp.

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

Tại sao phần mềm MGA++ tiêu chuẩn lại không tính toán được khi hàm lượng 239Pu dưới 75%?
Khi hàm lượng 239Pu giảm dưới 75%, mật độ và hoạt độ của các đồng vị 240Pu, 241Pu và 241Am tăng cao làm thay đổi căn bản cấu trúc vùng phổ 100 keV. Sự chồng chập nghiêm trọng của hơn 16 vạch phổ khiến thuật toán lặp đa nhóm tự động của MGA++ bị phân kỳ, đòi hỏi phải can thiệp bằng phương pháp xử lý đỉnh độc lập trên Gamma Vision.

Phương pháp phổ gamma trong luận văn có ưu điểm gì vượt trội so với phương pháp đo nhiệt lượng kế?
Phương pháp nhiệt lượng kế cần từ 2 đến 4 giờ để đạt trạng thái cân bằng nhiệt và chỉ đo được tổng công suất nhiệt sinh ra mà không phân biệt được các đồng vị riêng lẻ. Phổ kế gamma vừa xác định chính xác tỷ lệ từng đồng vị (238Pu đến 242Pu), vừa xác định tổng khối lượng với thời gian đo ngắn và sai số dưới 6%.

Đỉnh năng lượng nào mang tính quyết định nhất khi tính toán hàm lượng đồng vị 240Pu?
Đỉnh 160,28 keV và đỉnh 104,24 keV là hai đỉnh then chốt. Trong đó, đỉnh 160,28 keV (kèm phép giải chập loại trừ đóng góp của 241Pu tại 159,59 keV và 239Pu tại 160,19 keV) là công cụ hữu hiệu nhất khi detector không ghi nhận tối ưu ở dải năng lượng dưới 100 keV.

Làm thế nào để xác định khối lượng nguồn PuBe khi hoàn toàn không có thông số về kích thước vỏ bọc?
Quy trình sử dụng phương pháp năng lượng vô hạn kết hợp so sánh số đếm đỉnh 129,29 keV với hai đỉnh 375,04 keV và 413,71 keV để tối ưu hóa tỷ số hình học $R_{eff}$, kết hợp độ dày vỏ thép giả định 3,5 mm giúp xác định khối lượng với sai số chỉ từ 10% đến 15%.

Detector Planar HPGe GL2020R đóng vai trò như thế nào trong sự thành công của nghiên cứu?
Với diện tích tinh thể lớn 2000 mm2 và độ phân giải siêu cao đạt 490 eV tại đỉnh 122 keV, detector GL2020R cho phép bóc tách rõ nét các vạch phổ đơn lẻ trong các dải năng lượng phức tạp từ 40 keV đến 375 keV, bảo đảm số liệu thống kê đạt độ chính xác cao.

Kết luận

  • Nghiên cứu đã xây dựng thành công phương pháp phổ kế gamma độ phân giải cao kết hợp phần mềm Gamma Vision để xác định đồng thời thành phần đồng vị và khối lượng Plutonium trong nguồn PuBe.
  • Khắc phục triệt để hạn chế của phần mềm MGA++, xử lý chính xác các mẫu có hàm lượng 239Pu dưới 75% với độ lệch tỷ phần 239Pu dưới 10%.
  • Hoàn thiện công thức hiệu chỉnh tự hấp thụ phi tuyến $f(\mu_{lE} D)$, giúp hạ thấp sai số xác định khối lượng tổng Plutonium so với phương pháp nhiệt lượng kế xuống dưới 6%.
  • Đề xuất quy trình ước lượng tỷ số hình học hiệu dụng $R_{eff}$ qua đỉnh 129 keV, cho phép giám định hiệu quả các nguồn phóng xạ mất nhãn mác với sai số kiểm soát từ 10% đến 15%.
  • Đóng góp cơ sở lý thuyết và quy trình thực nghiệm giá trị cho công tác thanh sát an toàn bức xạ quốc gia giai đoạn 2026-2028. Bạn đọc và các nhà nghiên cứu quan tâm có thể tham khảo toàn văn luận văn để ứng dụng chi tiết các mô hình giải thuật vào thực tiễn đo lường hạt nhân.