Tổng quan nghiên cứu
Trong quản lý chu trình nhiên liệu hạt nhân, việc xác định độ cháy (burnup) giữ vai trò then chốt nhằm bảo đảm an toàn vận hành và tối ưu hóa hiệu quả kinh tế. Theo các tiêu chuẩn quốc tế, giới hạn độ cháy thanh nhiên liệu trong các lò phản ứng năng lượng hiện dao động từ 40 đến 65 GWd/t (điển hình tại Mỹ là 62 GWd/t, Đức là 65 GWd/t và Pháp là 52 GWd/t). Tại Việt Nam, việc giám sát độ cháy nhiên liệu của Lò phản ứng hạt nhân Đà Lạt (công suất danh định 500 kW) có ý nghĩa sống còn đối với công tác nạp đảo nhiên liệu và an toàn bức xạ.
Vấn đề nghiên cứu trọng tâm xuất phát từ hạn chế của các kỹ thuật phân tích truyền thống. Các phương pháp phá hủy mẫu như phân tích hóa phóng xạ hoặc khối phổ kế đòi hỏi quy trình phức tạp, chi phí đắt đỏ, thời gian kéo dài và phá vỡ cấu trúc cơ học của thanh nhiên liệu. Nhằm khắc phục nhược điểm này, luận văn tập trung vào mục tiêu nghiên cứu và ứng dụng phương pháp tỷ số đồng vị phân hạch không phá hủy mẫu, sử dụng hệ phổ kế gamma bán dẫn độ phân giải cao.
Phạm vi nghiên cứu được thực hiện trực tiếp trên bó thanh nhiên liệu VVR-M2 (bó số 62) chiếu xạ tại Lò phản ứng hạt nhân Đà Lạt. Kết quả của luận văn mang lại giá trị thực tiễn to lớn khi giúp giảm hơn 80% chi phí thử nghiệm so với kỹ thuật phá hủy, rút ngắn thời gian phân tích từ vài tuần xuống dưới 3 giờ làm việc, đồng thời duy trì 100% tính toàn vẹn vật lý của thanh nhiên liệu, hỗ trợ kiểm soát biên độ an toàn nhiệt DNBR trên mức tiêu chuẩn 1.15.
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 ứng dụng đồng thời hai khung lý thuyết cốt lõi trong vật lý lò phản ứng:
- Lý thuyết khuếch tán neutron đa nhóm và trạng thái tới hạn: Sử dụng mô hình khuếch tán 2 nhóm năng lượng (nhóm neutron nhanh $E > 0.1 \text{ MeV}$ và nhóm neutron nhiệt $E < 0.625 \text{ eV}$) để mô tả phân bố không - thời gian của thông lượng neutron $\Phi(r,t)$ và xác định hệ số nhân neutron hiệu dụng $k = 1$ tại trạng thái cân bằng tới hạn.
- Lý thuyết biến đổi thành phần hạt nhân và động học sản phẩm phân hạch: Thiết lập hệ phương trình vi phân mô tả chuỗi chuyển hóa của các đồng vị nặng ($^{235}\text{U}$, $^{238}\text{U}$, $^{239}\text{Pu}$) và sự tích tụ của các mảnh phân hạch.
- Các khái niệm chính: Độ cháy nhiên liệu ($BU$), Tốc độ sinh nhiệt tuyến tính (LHGR), Thông lượng nhiệt tới hạn (CHF), Tỷ số khoảng cách sôi hạt nhân (DNBR) và Hiện tượng tương tác cơ học vỏ - nhiên liệu (PCMI). Phản ứng phân hạch một hạt nhân $^{235}\text{U}$ giải phóng trung bình 200 MeV (trong đó động năng mảnh phân hạch chiếm 169 MeV và năng lượng gamma phân rã chiếm 6 MeV).
Phương pháp nghiên cứu
Nguồn dữ liệu thực nghiệm được thu thập từ bó nhiên liệu số 62 (loại cấu trúc hình trụ đồng tâm VVR-M2 với độ giàu ban đầu 36% $^{235}\text{U}$) được chiếu xạ qua nhiều chu kỳ vận hành tại Lò phản ứng Đà Lạt. Phương pháp chọn mẫu chủ đích được áp dụng nhằm chọn bó nhiên liệu có lịch sử chiếu xạ ổn định và đại diện cao cho cấu hình tâm lò.
Hệ thống đo đạc bức xạ sử dụng detector bán dẫn siêu tinh khiết HPGe (High-Purity Germanium) kết hợp với khối phân tích đa kênh MCA. Lý do lựa chọn hệ đo HPGe là vì thiết bị có độ phân giải năng lượng tuyệt vời (FWHM đạt mức 1.8 đến 2.0 keV tại đỉnh năng lượng 1332 keV của $^{60}\text{Co}$), cho phép bóc tách rõ ràng các đỉnh gamma đặc trưng trên 500 keV gồm $^{137}\text{Cs}$ (661.6 keV) và $^{134}\text{Cs}$ (604.7 keV và 795.8 keV).
Nghiên cứu sử dụng phương pháp chuẩn nội để xây dựng đường cong hiệu suất ghi tương đối mà không cần nguồn chuẩn hoạt độ tuyệt đối. Toàn bộ thuật toán giải hệ phương trình tích tụ đồng vị và khớp phổ gamma được lập trình tự động trên nền tảng phần mềm MATLAB. Toàn bộ quy trình thực nghiệm và tính toán mô phỏng được hoàn thiện trong chu kỳ nghiên cứu năm 2014 - 2015.
Kết quả nghiên cứu và thảo luận
Những phát hiện chính
Quá trình phân tích thực nghiệm và mô phỏng trên bó nhiên liệu số 62 đã mang lại các phát hiện khoa học quan trọng:
- Tương quan động học đồng vị theo thông lượng neutron: Mật độ hạt nhân $^{137}\text{Cs}$ tỷ lệ tuyến tính bậc nhất với thông lượng neutron và thời gian chiếu xạ do được tạo thành trực tiếp từ chuỗi phân rã phân hạch. Ngược lại, mật độ hạt nhân $^{134}\text{Cs}$ biến thiên theo hàm bậc hai của thông lượng neutron vì được tạo thành qua phản ứng bắt neutron cấp hai $^{133}\text{Cs}(n,\gamma)^{134}\text{Cs}$.
- Tính tuyến tính của tỷ số đồng vị $^{134}\text{Cs}/^{137}\text{Cs}$: Tỷ số diện tích đỉnh phổ gamma của hai đồng vị này có quan hệ tuyến tính hoàn hảo với độ cháy của thanh nhiên liệu, triệt tiêu hoàn toàn các sai số do suy giảm bức xạ trong vỏ bọc và sai số hình học đo.
- Xác định chính xác độ cháy bó nhiên liệu số 62: Độ cháy trung bình của bó nhiên liệu số 62 được xác định là 32.5% lượng $^{235}\text{U}$ đã tiêu hao (tương đương khoảng 28.6 MWd/kgHM). Kết quả này có độ sai lệch tương đối chỉ 4.8% khi đối chiếu với số liệu tham chiếu tính toán từ hồ sơ vận hành lò.
- Phân bố độ cháy dọc trục nhiên liệu: Phân bố hoạt độ bức xạ dọc theo chiều cao 60 cm của thanh nhiên liệu tuân theo dạng đường cong cosine, trong đó độ cháy tại vị trí trung tâm tâm lò (cách đáy 25-30 cm) cao hơn khoảng 22% so với hai đầu mút thanh nhiên liệu.
Thảo luận kết quả
Nguyên nhân chính giúp phương pháp tỷ số đồng vị đạt độ tin cậy cao là do tính chất hóa lý của Cesium. Trong điều kiện nhiệt độ vận hành của lò phản ứng nghiên cứu kiểu bể hở (nhiệt độ chất làm mát dưới 60 độ C), hiện tượng di cư nhiệt của Cs bên trong ma trận nhiên liệu oxit nhôm uran diễn ra không đáng kể.
Khi so sánh với phương pháp khối phổ kế phân tích cặp đồng vị Zirconium ($^{93}\text{Zr}/^{90}\text{Zr}$) hoặc Neodymium vốn có sai số khoảng 5.0%, phương pháp tỷ số đồng vị gamma không phá hủy mẫu đạt độ chính xác tương đương (sai số dưới 4.8%) nhưng vượt trội hoàn toàn về tốc độ phân tích và mức độ an toàn bức xạ.
Dữ liệu nghiên cứu được trình bày trực quan thông qua 44 hình vẽ, đồ thị phổ đa kênh và 16 bảng biểu chi tiết. Sự biến thiên của tốc độ đếm xung theo năng lượng và phân bố công suất thể tích được minh họa rõ ràng qua các biểu đồ phân bố không gian, giúp các kỹ sư vận hành dễ dàng nhận diện vùng suy giảm nhiên liệu cục bộ.
Đề xuất và khuyến nghị
Nhằm chuẩn hóa quy trình giám sát độ cháy nhiên liệu và nâng cao năng lực an toàn hạt nhân, tác giả đưa ra các khuyến nghị cụ thể:
- Tích hợp module tính toán tự động vào hệ thống điều khiển lò: Ban Quản lý Lò phản ứng Đà Lạt cần chuyển giao và tích hợp thuật toán phân tích tỷ số đồng vị viết trên MATLAB vào hệ thống thu nhận dữ liệu SCADA, nhằm tự động hóa quy trình đánh giá độ cháy thanh nhiên liệu trong vòng 15 phút sau mỗi đợt đo, hoàn thành trước quý 4 năm 2026.
- Nâng cấp ống chuẩn trực và hệ thống đo phổ gamma: Viện Năng lượng Nguyên tử Việt Nam cần đầu tư hệ thống collimator chì đa khe đo (từ 0.5 mm đến 2.0 mm) tích hợp cơ cấu quét tự động theo trục đứng, nhằm nâng cao tỷ số tín hiệu trên nhiễu (SNR) thêm ít nhất 25% trong giai đoạn 2026 - 2027.
- Xây dựng cơ sở dữ liệu số hóa lịch sử vận hành: Bộ phận an toàn bức xạ cần số hóa 100% lịch sử nạp đảo vùng hoạt, cập nhật thông lượng neutron 2 nhóm theo từng chu kỳ vận hành 100 giờ để giảm sai số tính toán mật độ ban đầu xuống dưới 2.0%.
- Mở rộng nghiên cứu sang các cặp tỷ số đồng vị phụ trợ: Nhóm nghiên cứu tại các viện và trường đại học chuyên ngành nên mở rộng đánh giá tỷ số $^{154}\text{Eu}/^{137}\text{Cs}$ và $^{106}\text{Ru}/^{137}\text{Cs}$ đối với các thanh nhiên liệu có độ cháy cao vượt ngưỡng 50 GWd/t, hoàn thành báo cáo thử nghiệm trong vòng 18 tháng tới.
Đối tượng nên tham khảo luận văn
Nội dung luận văn cung cấp nền tảng lý thuyết và số liệu thực nghiệm giá trị cho 4 nhóm đối tượng chính:
- Kỹ sư vận hành và quản lý nhiên liệu lò phản ứng: Khai thác quy trình đo không phá hủy để kiểm soát chính xác độ cháy từng bó nhiên liệu, tối ưu hóa sơ đồ nạp đảo tâm lò và duy trì hệ số DNBR an toàn trên 1.15.
- Học viên cao học và nghiên cứu sinh ngành Vật lý Hạt nhân: Sử dụng làm tài liệu nghiên cứu chuyên sâu về lý thuyết khuếch tán neutron đa nhóm, động học biến đổi đồng vị trong nhiên liệu và kỹ thuật phổ kế gamma bán dẫn HPGe.
- Cơ quan quản lý và thanh sát an toàn hạt nhân (VARANS): Tham chiếu phương pháp luận và số liệu thực nghiệm để xây dựng quy chuẩn kỹ thuật quốc gia về kiểm toán vật liệu hạt nhân và giám sát thanh nhiên liệu đã qua sử dụng.
- Chuyên gia mô phỏng và phát triển phần mềm hạt nhân: Khai thác các hệ phương trình vi phân và thuật toán xử lý dữ liệu để xây dựng các phần mềm mô phỏng vật lý - nhiệt thủy lực tích hợp cho lò phản ứng.
Câu hỏi thường gặp
Tại sao phương pháp tỷ số đồng vị Cs-134/Cs-137 lại ưu việt hơn đo hoạt độ tuyệt đối?
Phương pháp tỷ số đồng vị tự triệt tiêu các yếu tố gây sai số hệ thống như sai số hình học đo, khoảng cách từ nguồn tới detector và độ suy giảm bức xạ qua lớp vỏ bọc kim loại. Nhờ việc lập tỷ số giữa hai đỉnh năng lượng lân cận (604.7 keV hoặc 795.8 keV của $^{134}\text{Cs}$ so với 661.6 keV của $^{137}\text{Cs}$), kết quả đo đạt độ tin cậy cao với sai số dưới 4.8% mà không cần xác định khối lượng mẫu tuyệt đối.
Độ cháy nhiên liệu ảnh hưởng trực tiếp đến an toàn lò phản ứng như thế nào?
Khi độ cháy tăng cao, độ dẫn nhiệt của viên nhiên liệu giảm từ khoảng 8 W/m·K xuống dưới 3 W/m·K do sự tích tụ khí phân hạch Xe và Kr. Điều này làm tăng nhiệt độ tâm nhiên liệu, gia tăng áp suất nội tại của thanh và đẩy hệ số DNBR lại gần ngưỡng tới hạn 1.15. Do đó, kiểm soát độ cháy trong giới hạn 50-65 GWd/t giúp ngăn ngừa nguy cơ nứt vỡ vỏ bọc trong trường hợp xảy ra sự cố mất chất tải nhiệt (LOCA).
Vai trò của đầu dò bán dẫn siêu tinh khiết HPGe trong phép đo là gì?
Detector HPGe sở hữu độ phân giải năng lượng vượt trội với độ rộng bán phần FWHM dưới 2.0 keV tại đỉnh 1332 keV, cao gấp hơn 25 lần so với detector nhấp nháy NaI(Tl). Nhờ đó, hệ thống phân tách rạch ròi các đỉnh năng lượng phức tạp của $^{134}\text{Cs}$, $^{137}\text{Cs}$ và $^{154}\text{Eu}$ trong trường bức xạ gamma hỗn hợp từ thanh nhiên liệu đã qua chiếu xạ.
Quy trình nghiên cứu trong luận văn có thể áp dụng cho nhà máy điện hạt nhân thương mại không?
Hoàn toàn có thể áp dụng. Dù thực nghiệm thực hiện trên bó nhiên liệu VVR-M2 của lò nghiên cứu Đà Lạt, hệ phương trình vi phân biến đổi đồng vị và nguyên lý đo tỷ số Cesium mang tính phổ quát, áp dụng trực tiếp cho các thanh nhiên liệu $UO_2$ và MOX trong các lò phản ứng năng lượng loại PWR, BWR hoặc VVER có độ cháy từ 40 đến 65 GWd/t.
Tại sao cần xây dựng đường cong hiệu suất ghi tương đối bằng phương pháp chuẩn nội?
Đường cong hiệu suất tương đối giúp hiệu chỉnh xác suất ghi nhận hạt photon gamma của detector ở các mức năng lượng khác nhau từ 500 keV đến 1500 keV mà không cần nguồn chuẩn phóng xạ đa năng lượng đắt tiền. Kỹ thuật này giúp tiết kiệm hơn 70% thời gian chuẩn bị và loại bỏ hoàn toàn sai số dịch chuyển vị trí nguồn đo.
Kết luận
- Luận văn đã hệ thống hóa toàn diện cơ sở lý thuyết về vật lý khuếch tán neutron đa nhóm và động học tích lũy sản phẩm phân hạch trong thanh nhiên liệu hạt nhân.
- Xây dựng thành công quy trình xác định độ cháy nhiên liệu bằng phương pháp tỷ số đồng vị $^{134}\text{Cs}/^{137}\text{Cs}$ không phá hủy mẫu, kết hợp phổ kế gamma bán dẫn HPGe.
- Ứng dụng thực nghiệm trên bó nhiên liệu số 62 tại Lò phản ứng hạt nhân Đà Lạt cho kết quả độ cháy đạt 32.5% $^{235}\text{U}$ với sai số thấp dưới 4.8%.
- Xác lập chính xác đường cong phân bố độ cháy và hoạt độ phóng xạ dọc trục thanh nhiên liệu theo hàm dạng cosine.
- Cung cấp giải pháp kỹ thuật có giá thành rẻ, tốc độ phân tích nhanh và bảo đảm tuyệt đối tính toàn vẹn cơ học của nhiên liệu hạt nhân.
Đóng góp lớn nhất của công trình là làm chủ công nghệ đo đạc hạt nhân tiên tiến, tạo tiền đề kỹ thuật vững chắc cho việc quản lý an toàn vùng hoạt lò phản ứng tại Việt Nam. Trong giai đoạn 2026 - 2028, các đơn vị vận hành cần tiếp tục tự động hóa phần mềm xử lý dữ liệu và mở rộng quy trình kiểm tra cho toàn bộ các bó nhiên liệu trong vùng hoạt. Quý độc giả và các nhà nghiên cứu quan tâm có thể khai thác tài liệu này để phục vụ công tác giảng dạy, nghiên cứu và ứng dụng thực tiễn trong ngành kỹ thuật hạt nhân.