Tổng quan nghiên cứu

Trong ngành công nghiệp năng lượng hiện đại, nhiên liệu hạt nhân đóng vai trò trụ cột nhằm bảo đảm an ninh năng lượng quốc gia trước sự cạn kiệt của các nguồn nhiên liệu hóa thạch truyền thống. Trong tự nhiên, nguyên tố Urani tồn tại chủ yếu dưới ba đồng vị phóng xạ gồm 238U chiếm 99,2745%, 235U chiếm 0,720% và 234U chỉ chiếm 0,0055% với hàm lượng trung bình trong vỏ Trái Đất khoảng 10^-4% (tương đương vài ppm). Để phục vụ cho các lò phản ứng năng lượng thương mại, Urani tự nhiên bắt buộc phải trải qua quá trình làm giàu đồng vị 235U lên mức 3% đến 5%. Cứ từ 8,05 tấn Urani tự nhiên chứa 0,72% 235U, quy trình công nghệ sẽ sản xuất được 1 tấn Urani làm giàu đạt 3,6% 235U và thải ra 7,05 tấn Urani nghèo chứa khoảng 0,3% 235U.

Vấn đề cốt lõi đặt ra trong công tác quản lý an toàn, kiểm kê và thanh sát hạt nhân là xác định chính xác độ giàu đồng vị 235U mà không làm tổn hại đến cấu trúc thanh nhiên liệu. Các phương pháp phân tích phá hủy mẫu truyền thống như khối phổ kế cảm ứng plasma (ICP-MS) đòi hỏi quy trình nghiền mẫu phức tạp, chi phí cao và tiềm ẩn rủi ro phát tán phóng xạ. Vì vậy, nghiên cứu này tập trung giải quyết bài toán xác định các đặc trưng cơ bản của thanh nhiên liệu hạt nhân và định lượng độ giàu Urani bằng phương pháp phân tích không phá hủy mẫu (NDA) sử dụng phổ kế bức xạ gamma.

Mục tiêu cụ thể của luận văn bao gồm: hệ thống hóa các đặc trưng bức xạ và cơ nhiệt của thanh nhiên liệu trong các lò phản ứng PWR và VVER-1000; xây dựng cơ sở phương pháp phổ kế gamma dải năng lượng thấp dưới 300 keV kết hợp kỹ thuật chuẩn trong (nội chuẩn); và tiến hành đo đạc thực nghiệm trên các mẫu nhiên liệu làm giàu thấp nhằm đánh giá độ chính xác và sai số của phương pháp. Nghiên cứu được thực hiện tại Bộ môn Vật lý Hạt nhân, Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội và Viện Khoa học Đồng vị Phóng xạ Hungary trong giai đoạn 2015-2016. Kết quả nghiên cứu mang lại ý nghĩa khoa học và thực tiễn sâu sắc, cung cấp giải pháp đo lường nhanh với sai số tương đối dưới 1,0%, hỗ trợ đắc lực cho công tác kiểm định chất lượng và thanh sát an toàn hạt nhân.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nền tảng vật lý của đề tài dựa trên lý thuyết phân rã phóng xạ tự nhiên của hai chuỗi Urani chính: chuỗi 238U - 206Pb (có số khối xác định theo biểu thức A = 4n + 2) và chuỗi 235U - 207Pb (có số khối xác định theo biểu thức A = 4n + 3). Chu kỳ bán rã của 238U là khoảng 4,47 tỷ năm, trong khi chu kỳ bán rã của 235U là khoảng 704 triệu năm và 234U là khoảng 2,46 x 10^5 năm. Do tuổi của các mẫu quặng và vật liệu Urani rất lớn so với chu kỳ bán rã của các hạt nhân con cháu, hệ thống đạt đến trạng thái cân bằng phóng xạ thế kỷ. Tại trạng thái này, hoạt độ phóng xạ của các đồng vị trong cùng một chuỗi đều bằng nhau theo phương trình cân bằng liên tiếp.

Về mặt cơ nhiệt, mô hình truyền nhiệt trong thanh nhiên liệu dạng trụ được xác lập thông qua phương trình phân bố nhiệt độ hướng tâm phụ thuộc công suất sinh nhiệt thể tích. Trong quá trình vận hành ở công suất danh định 100%, độ chênh lệch nhiệt độ giữa tâm và bề mặt viên gốm nhiên liệu UO2 lên tới 1.200°C trên khoảng cách bán kính chỉ 3,8 mm. Nhiệt độ cực đại tại tâm viên gốm có thể đạt từ 1.600°C đến 1.940°C, trong khi bề mặt ngoài vỏ bọc hợp kim Zirconi duy trì ở mức khoảng 350°C.

Các khái niệm cốt lõi được áp dụng trong luận văn bao gồm:

  1. Độ giàu Urani: Tỷ lệ phần trăm khối lượng của đồng vị phân hạch 235U so với tổng khối lượng Urani trong mẫu vật liệu.
  2. Phân tích không phá hủy mẫu (NDA): Kỹ thuật nhận diện và định lượng đồng vị phóng xạ trực tiếp trên nguyên mẫu mà không làm thay đổi trạng thái lý - hóa của vật liệu.
  3. Kỹ thuật chuẩn trong (nội chuẩn): Phương pháp xác định hàm lượng đồng vị thông qua việc xây dựng đường cong hiệu suất ghi tương đối từ chính các vạch phổ gamma của mẫu đo, loại bỏ sự phụ thuộc vào dạng hình học và kích thước mẫu chuẩn bên ngoài.
  4. Quá trình rão và mỏi nhiệt: Hiện tượng biến dạng dẻo chậm và suy giảm độ bền cơ học của viên gốm UO2 cùng vỏ bọc Zircaloy/E110 dưới tác động liên tục của ứng suất nhiệt và mật độ thông lượng nơtron lên tới 10^15 n/cm2.s.

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

Nguồn dữ liệu thực nghiệm của luận văn được thu thập từ việc đo đạc 02 mẫu nhiên liệu Urani tiêu chuẩn: mẫu U0.7 là Urani tự nhiên tinh chế từ quặng do Viện Công nghệ Xạ hiếm (Viện Năng lượng Nguyên tử Việt Nam) cung cấp, và mẫu U4.46 là Urani làm giàu thấp do Cơ quan Năng lượng Nguyên tử Quốc tế (IAEA) cung cấp cho Viện Khoa học Đồng vị Phóng xạ Hungary. Mỗi mẫu có khối lượng chuẩn là 5,0 gam, được đóng kín trong hộp trụ kích thước đường kính 3,0 cm.

Phương pháp phân tích quang phổ gamma sử dụng hai hệ phổ kế bán dẫn siêu tinh khiết có độ phân giải cao:

  • Hệ phổ kế bán dẫn dải rộng BEGe (Broad Energy Germanium) Model BE5030 của hãng Canberra tại Đại học Khoa học Tự nhiên Hà Nội. Đầu dò có tinh thể bán dẫn đường kính 80,5 mm, độ dày 31 mm, diện tích hiệu dụng 5.000 mm2. Độ phân giải năng lượng (FWHM) tại các đỉnh 5,9 keV, 122 keV và 1.332 keV lần lượt đạt 0,5 keV; 0,75 keV và 2,2 keV. Đầu dò được làm lạnh liên tục bằng điện ở nhiệt độ 90 K và che chắn bằng buồng chì có tổng khối lượng 950 kg (gồm vỏ thép 9,5 mm, chì phông thấp 10 cm, lớp kẽm 1,0 mm và đồng tinh khiết 1,6 mm).
  • Hệ phổ kế bán dẫn HPGe bản mỏng Model GLP-10180/07 của hãng ORTEC tại Hungary, với diện tích hoạt động 80 mm2, chiều sâu tinh thể 7,0 mm, cửa sổ Beryli dày 0,13 mm và độ phân giải 180 eV tại đỉnh 5,9 keV.

Cỡ mẫu nghiên cứu gồm 02 mẫu đại diện cho hai mức độ giàu điển hình (tự nhiên 0,72% và làm giàu thấp 4,46%). Phương pháp chọn mẫu có chủ đích được áp dụng nhằm kiểm chứng tính đúng đắn của phương pháp chuẩn trong trên cả dải Urani tự nhiên lẫn Urani đã qua chế biến làm giàu. Lý do lựa chọn phương pháp phân tích phổ kế gamma bán dẫn dải năng lượng thấp dưới 300 keV là vì kỹ thuật này tận dụng trực tiếp các vạch phát xạ gamma tự nhiên của 235U và các đồng vị con cháu ngắn ngày của 238U (như 234Th và 234mPa), giúp giảm thiểu tối đa hiện tượng tự hấp thụ bức xạ cứng và không đòi hỏi nguồn kích hoạt ngoài. Quá trình thu thập dữ liệu phổ được tiến hành với thời gian đo thực tế là 171.650,5 giây đối với mẫu U0.7 và 74.281,46 giây đối với mẫu U4.46.

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

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

Thứ nhất, nghiên cứu đã xây dựng thành công hàm hiệu suất ghi tương đối trong vùng năng lượng thấp dưới 300 keV dựa trên các đỉnh phổ gamma đặc trưng của đồng vị 235U. Cụ thể, các vạch năng lượng chính được nhận diện và xử lý bao gồm 143,76 keV (cường độ 10,96%), 163,33 keV (cường độ 5,08%), 185,715 keV (cường độ phát xạ mạnh nhất đạt 57,2%) và 205,311 keV (cường độ 5,01%). Đối với 238U, vạch phổ 258,227 keV của 234mPa và các vạch 63,29 keV, 92,38 keV của 234Th được dùng làm cơ sở đối chuẩn tỷ số hoạt độ.

Thứ hai, việc xác định độ giàu đồng vị 235U trong hai mẫu nhiên liệu thực nghiệm đạt độ chính xác rất cao so với giá trị danh định chuẩn:

  • Đối với mẫu quặng tự nhiên U0.7, kết quả tính toán độ giàu 235U đạt mức xấp xỉ 0,718% đến 0,720%, sai lệch không đáng kể so với tỷ lệ chuẩn tự nhiên 0,720%.
  • Đối với mẫu nhiên liệu làm giàu U4.46, kết quả định lượng đạt giá trị 4,46% với sai số tương đối của phép đo được kiểm soát chặt chẽ ở mức dưới 1,0%.

Thứ ba, nghiên cứu đã làm rõ quy luật biến đổi vi cấu trúc và phân bố nhiệt cơ trong viên gốm nhiên liệu UO2 của lò phản ứng VVER-1000. Dưới tác động của thông lượng bức xạ gamma 10^13 photon/cm2.s và nơtron 10^15 n/cm2.s, độ chênh lệch nhiệt độ xuyên tâm 1.200°C dẫn đến hiện tượng tái cấu trúc tinh thể gốm trong vòng 24 giờ đầu vận hành. Quá trình ép nóng trong 200 giờ đầu tiên làm tiêu biến các lỗ xốp nhỏ dưới 1,0 μm, nâng mật độ vật liệu lên mức 98% đến 99% mật độ lý thuyết và làm giảm chiều cao cột nhiên liệu tới 7,5 cm so với chiều cao thiết kế ban đầu 3,66 m.

Thứ tư, hệ phổ kế gamma BEGe với buồng chì che chắn 950 kg đã thể hiện hiệu suất ghi vượt trội, giảm phông phóng xạ môi trường tới hơn 90% ở vùng năng lượng thấp. Độ phân giải năng lượng 0,75 keV tại vạch 122 keV cho phép tách biệt hoàn toàn các đỉnh năng lượng liền kề như 92,38 keV và 92,80 keV của 234Th với đỉnh tia X 93,356 keV của 231Th.

Thảo luận kết quả

Độ chính xác cao của phương pháp chuẩn trong bắt nguồn từ việc triệt tiêu được các yếu tố hình học phức tạp giữa mẫu và đầu dò. Khi tính toán tỷ số hoạt độ giữa 235U và 238U trên cùng một phổ đo, góc khối và hệ số truyền dẫn bức xạ bị triệt tiêu hoàn toàn trong biểu thức giải tích. Điều này khắc phục được nhược điểm cố hữu của phương pháp đo chuẩn ngoài vốn đòi hỏi mẫu chuẩn phải có cùng kích thước, mật độ và hình dạng hình học chính xác với mẫu cần phân tích.

Khi so sánh với phương pháp phân tích khí UF6 bằng chùm chiếu 241Am (năng lượng 60 keV) có sai số xấp xỉ 1% ở áp suất trên 700 torr, phương pháp phổ kế gamma BEGe thể hiện tính ưu việt hơn khi không phụ thuộc vào trạng thái áp suất và áp dụng được trực tiếp cho nhiên liệu dạng rắn. Kết quả phân tích cũng tương đương với độ tin cậy của phương pháp khối phổ kế plasma ICP-MS nhưng tiết kiệm hơn 60% thời gian chuẩn bị mẫu và loại bỏ hoàn toàn việc phát sinh chất thải phóng xạ lỏng độc hại.

Dữ liệu thực nghiệm của nghiên cứu có thể được biểu diễn trực quan thông qua đồ thị hồi quy phi tuyến mô tả sự phụ thuộc của hàm hiệu suất vào các mức năng lượng từ 50 keV đến 300 keV, kết hợp với bảng số liệu tổng hợp diện tích đỉnh phổ thực tế, tốc độ đếm và sai số thống kê tương ứng. Sự đồng nhất giữa đường cong thực nghiệm và mô hình lý thuyết khẳng định tính đúng đắn của thuật toán khớp phổ Gauss kết hợp hiệu chỉnh đuôi hàm mũ.

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

Nhằm nâng cao hiệu quả kiểm soát an toàn hạt nhân và ứng dụng thực tiễn các kết quả nghiên cứu, bốn nhóm giải pháp chiến lược được đề xuất như sau:

  1. Ban hành quy trình kiểm tra không phá hủy mẫu chuẩn quốc gia: Thiết lập quy chuẩn kỹ thuật đo độ giàu nhiên liệu Urani bằng phổ kế gamma bán dẫn dải năng lượng thấp (<300 keV) cho toàn bộ các cơ sở nghiên cứu và đào tạo hạt nhân. Mục tiêu đạt sai số phân tích dưới 0,5% đối với các lô nhiên liệu nhập khẩu, hoàn thành lộ trình ban hành trong vòng 12 tháng dưới sự chủ trì của Cục An toàn Bức xạ và Hạt nhân.
  2. Đầu tư hiện đại hóa hệ thống đo phổ gamma độ phân giải cao: Trang bị đồng bộ hệ thống đầu dò bán dẫn BEGe hoặc HPGe bản mỏng có diện tích hiệu dụng từ 5.000 mm2 trở lên, tích hợp buồng chì che chắn phông thấp dày tối thiểu 10 cm cho các phòng thí nghiệm phân tích hạt nhân trọng điểm. Tiến độ giải ngân và lắp đặt thực hiện trong giai đoạn 18 tháng do Viện Năng lượng Nguyên tử Việt Nam điều phối.
  3. Ứng dụng mô hình tính toán nhiệt cơ vào quản lý an toàn vùng hoạt: Tích hợp các thông số phân bố nhiệt độ (gradient nhiệt 1.200°C) và mức độ co ngót cột nhiên liệu (khoảng 7,5 cm) vào phần mềm mô phỏng vận hành lò phản ứng nghiên cứu, giúp kiểm soát nguy cơ rão mỏi vỏ bọc Zircaloy/E110 và giảm thiểu xác suất rò rỉ phóng xạ xuống dưới 0,01% trong chu kỳ vận hành 24 tháng.
  4. Xây dựng cơ sở dữ liệu phổ chuẩn đồng vị Urani quốc gia: Chuẩn hóa thư viện số liệu về tỷ số rẽ nhánh và các vạch bức xạ gamma đặc trưng của các hợp chất oxit Urani (UO2, U3O8, UF4) nhằm rút ngắn thời gian xử lý dữ liệu phổ xuống dưới 12 giờ mỗi mẫu, giao Bộ môn Vật lý Hạt nhân - Trường Đại học Khoa học Tự nhiên phối hợp cùng các viện nghiên cứu chuyên ngành hoàn thiện trong 6 tháng.

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

Nội dung và kết quả thực nghiệm của luận văn mang lại giá trị thực tiễn cho 4 nhóm đối tượng chính:

  1. Cơ quan quản lý nhà nước và thanh sát an toàn hạt nhân: Cung cấp cơ sở khoa học và quy trình kỹ thuật để thanh tra, kiểm kê độ giàu Urani (từ 0,72% đến 5,0%), phục vụ công tác giám sát chống phổ biến vũ khí hạt nhân và quản lý an toàn phóng xạ quốc gia.
  2. Kỹ sư vận hành và chuyên gia thiết kế lò phản ứng: Cung cấp tập dữ liệu chi tiết về đặc trưng bức xạ, nhiệt trở khe hở khí Heli và cơ chế mỏi rão của vỏ bọc nhiên liệu E110/Zircaloy dưới thông lượng nơtron 10^15 n/cm2.s, phục vụ tối ưu hóa cấu hình vùng hoạt và kéo dài tuổi thọ thanh nhiên liệu.
  3. Giảng viên, nghiên cứu sinh và học viên ngành Vật lý hạt nhân: Làm tài liệu tham khảo chuyên sâu về phương pháp giải phổ gamma, thuật toán khớp đỉnh Gauss phi tuyến và kỹ thuật chuẩn trong ứng dụng trên hệ phổ kế bán dẫn Ge siêu tinh khiết.
  4. Chuyên viên phân tích tại các phòng thí nghiệm đo lường phóng xạ: Hướng dẫn chi tiết quy trình vận hành đầu dò BEGe Model BE5030, kỹ thuật thiết lập hệ điện tử xử lý xung đa kênh và phương pháp làm lạnh đầu dò về 90 K nhằm tối ưu hóa độ phân giải FWHM trong phân tích mẫu phóng xạ thấp.

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

  1. Phương pháp phổ kế gamma chuẩn trong có ưu điểm gì vượt trội so với phương pháp chuẩn ngoài truyền thống? Phương pháp chuẩn trong xây dựng đường cong hiệu suất trực tiếp từ các đỉnh gamma của chính mẫu đo, giúp loại bỏ hoàn toàn yêu cầu về mẫu chuẩn có cùng hình học và mật độ. Khi đo hai vạch gamma trong cùng dải hẹp dưới 300 keV, hệ số góc khối và độ truyền dẫn qua vỏ bọc bị triệt tiêu, giúp phép đo đạt sai số dưới 1,0% ngay cả với các mẫu có hình dạng bất kỳ.

  2. Tại sao dải năng lượng dưới 300 keV lại được lựa chọn để xác định độ giàu Urani trong nghiên cứu này? Vùng năng lượng dưới 300 keV tập trung hầu hết các vạch phát xạ gamma đặc trưng nhất của 235U như 143,76 keV, 163,33 keV và vạch 185,715 keV có cường độ phát xạ lớn tới 57,2%. Đồng thời, các đồng vị con cân bằng của 238U như 234Th (63,29 keV) và 234mPa (258,227 keV) cũng phát photon trong dải này, tạo điều kiện thuận lợi để so sánh tỷ số hoạt độ trực tiếp giữa hai chuỗi phân rã.

  3. Sự biến đổi hình dạng và mật độ viên gốm UO2 trong lò phản ứng ảnh hưởng như thế nào đến an toàn hạt nhân? Gradient nhiệt độ lên tới 1.200°C giữa tâm (1.600°C) và bề mặt viên gốm gây ra ứng suất nhiệt cực lớn, dẫn đến nứt nẻ và làm giảm chiều cao cột nhiên liệu tới 7,5 cm trên tổng chiều cao 3,66 m. Hiện tượng này làm thay đổi mật độ công suất cục bộ và gia tăng áp lực khí phân hạch lên vỏ hợp kim E110, đòi hỏi phải tính toán chính xác khe hở để tránh vỡ vỏ bọc.

  4. Đầu dò bán dẫn BEGe có điểm gì khác biệt so với đầu dò HPGe đồng trục thông thường? Đầu dò BEGe (Model BE5030) có cấu trúc tinh thể dạng bản mỏng to bản với diện tích hiệu dụng 5.000 mm2, tối ưu hóa hiệu suất ghi cho photon năng lượng thấp dưới 1 MeV. Nhờ độ phân giải vượt trội 0,75 keV tại đỉnh 122 keV và phông nền thấp hơn hẳn đầu dò đồng trục, BEGe giúp tách rõ các đỉnh phổ năng lượng sát nhau trong dải 50 keV đến 300 keV.

  5. Mối liên hệ giữa chu kỳ bán rã và khả năng phân hạch của 235U và 238U được ứng dụng như thế nào? Đồng vị 235U có chu kỳ bán rã 704 triệu năm, sở hữu tiết diện phản ứng phân hạch nơtron nhiệt rất lớn, đóng vai trò tạo năng lượng chính trong lò phản ứng thương mại ở độ giàu 3% đến 5%. Ngược lại, 238U có chu kỳ bán rã 4,47 tỷ năm, chỉ phân hạch với nơtron nhanh nhưng có khả năng hấp thụ nơtron để chuyển hóa thành 239Pu, tạo nguồn nhiên liệu thứ cấp quan trọng.

Kết luận

  • Luận văn đã phân tích toàn diện các đặc trưng vật lý bức xạ, cơ nhiệt và hiện tượng rão mỏi của thanh nhiên liệu hạt nhân UO2 trong cấu trúc lò phản ứng VVER-1000 và PWR.
  • Thiết lập thành công cơ sở lý thuyết và quy trình thực nghiệm xác định độ giàu Urani bằng phương pháp phổ kế gamma không phá hủy mẫu dựa trên kỹ thuật chuẩn trong ở dải năng lượng dưới 300 keV.
  • Ứng dụng đầu dò bán dẫn BE5030 Canberra và HPGe GLP-10180/07 ORTEC đo đạc chính xác hai mẫu thực nghiệm U0.7 và U4.46 với thời gian đo lên tới 171.650,5 giây, đạt sai số tương đối dưới 1,0%.
  • Chứng minh tính khả thi của việc sử dụng các vạch gamma 185,715 keV (235U) và 258,227 keV (234mPa/238U) để xác định tỷ số hoạt độ mà không cần chế tạo mẫu chuẩn hình học tương đồng.
  • Đóng góp luận cứ khoa học quan trọng cho công tác kiểm định nhiên liệu hạt nhân, hỗ trợ đắc lực cho hoạt động nghiên cứu vật lý nguyên tử và đào tạo chuyên sâu tại Việt Nam.

Đóng góp chính của công trình là chuẩn hóa phương pháp phân tích gamma năng lượng thấp cho vật liệu hạt nhân trong nước. Kế hoạch tiếp theo trong vòng 12 tháng tới là mở rộng nghiên cứu xác định tuổi nhiên liệu Urani và đo đạc các mẫu nhiên liệu cháy có độ giàu cao hơn. Quý độc giả và các nhà nghiên cứu quan tâm vui lòng liên hệ Bộ môn Vật lý Hạt nhân, Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội để cùng hợp tác phát triển và ứng dụng công nghệ.