Tổng quan nghiên cứu
Độ cháy nhiên liệu hạt nhân (Burnup - BU) là một trong những đại lượng kỹ thuật mang tính cốt lõi trong chu trình quản lý nhiên liệu và vận hành an toàn nhà máy điện hạt nhân cũng như các lò phản ứng nghiên cứu. Trên bình diện quốc tế, giới hạn độ cháy thiết kế của các thanh nhiên liệu UO2 thường dao động trong khoảng từ 50 GWd/t đến 65 GWd/t, điển hình như tại Hoa Kỳ với mức 62 GWd/t hay CHLB Đức đạt tới 65 GWd/t đối với lò nước áp lực (PWR). Việc xác định chính xác độ cháy đóng vai trò quyết định đến việc tối ưu hóa hiệu quả kinh tế phân hạch, kiểm soát độ làm giàu ban đầu (thường được giới hạn ở ngưỡng dưới 5%) và ngăn ngừa nguy cơ hư hỏng vỏ bọc do tích tụ áp suất khí hay biến dạng cơ học vượt mức 1%.
Vấn đề then chốt được đặt ra là làm thế nào để đo đạc và xác định độ cháy một cách chính xác mà không làm phá hủy cấu trúc vật lý của thanh nhiên liệu vốn có hoạt độ phóng xạ rất cao. Luận văn thạc sĩ chuyên ngành Vật lý Nguyên tử và Hạt nhân (mã số: 60440106) tại Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội tập trung giải quyết bài toán này thông qua mục tiêu phát triển và hoàn thiện quy trình ứng dụng phương pháp tỷ số đồng vị phân hạch không phá hủy mẫu.
Phạm vi thực nghiệm và tính toán của nghiên cứu được triển khai trực tiếp trên bó nhiên liệu số 62 loại VVR-M2 vận hành tại Lò phản ứng hạt nhân Đà Lạt. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn khi cung cấp giải pháp xác định độ cháy với độ tin cậy cao, loại bỏ hoàn toàn chi phí đắt đỏ của các phân tích hóa học phá hủy mẫu, đồng thời hỗ trợ giám sát phân bố công suất dọc theo trục thanh nhiên liệu nhằm duy trì tỷ số khoảng cách sôi màng (DNBR) trong ngưỡng an toàn tuyệt đối.
Cơ sở lý thuyết và phương pháp nghiên cứu
Khung lý thuyết áp dụng
Khung lý thuyết của luận văn được xây dựng trên nền tảng cơ học lượng tử hạt nhân, lý thuyết khuếch tán neutron đa nhóm và các mô hình truyền nhiệt - thủy lực trong tâm lò phản ứng. Trong quá trình lò hoạt động, sự phân hạch hạt nhân Uranium-235 cùng với phản ứng bắt neutron của Uranium-238 để chuyển hóa thành Plutonium-239 giải phóng một mức năng lượng trung bình khoảng 200 MeV trên mỗi phân hạch. Năng lượng này chuyển hóa thành nhiệt lượng thể tích, được truyền qua viên nhiên liệu UO2 tới lớp vỏ bọc Zirconium theo định luật dẫn nhiệt Fourier.
Bên cạnh lý thuyết phân hạch, nghiên cứu tích hợp sâu các khái niệm an toàn thủy nhiệt hạt nhân bao gồm:
- Tốc độ sinh nhiệt tuyến tính (LHGR) và thông lượng nhiệt tới hạn (CHF).
- Tỷ số khoảng cách tính từ điểm sôi hạt nhân (DNBR) và khủng hoảng sôi theo đường cong Ukiyama.
- Tương tác cơ học giữa nhiên liệu và vỏ bọc (PCMI) cùng sự tích lũy lớp cặn oxit (CRUD).
- Tỷ số đồng vị phân hạch phóng xạ đặc trưng như Cs-134/Cs-137 và Eu-154/Cs-137.
Phương pháp nghiên cứu
Nghiên cứu sử dụng phương pháp phân tích phổ kế gamma không phá hủy mẫu (Non-Destructive Assay - NDA) kết hợp mô hình hóa số trị. Nguồn dữ liệu thực nghiệm được thu thập từ hệ phổ kế gamma bán dẫn siêu tinh khiết HPGe (High-Purity Germanium) có độ phân giải năng lượng cao, ghi nhận chùm bức xạ phát ra từ các đồng vị sản phẩm phân hạch của bó nhiên liệu số 62.
Cỡ mẫu nghiên cứu tập trung vào bó nhiên liệu đa thanh VVR-M2 với nhiều điểm đo phân bố dọc theo chiều cao hoạt nhiệt của thanh. Phương pháp lấy mẫu có chủ đích được lựa chọn để khảo sát các vùng có gradient thông lượng neutron biến thiên mạnh nhất. Lý do lựa chọn phương pháp đo phổ gamma tỷ số đồng vị Cs-134/Cs-137 thay vì phương pháp khối phổ kế truyền thống là bởi vì phương pháp này cho phép đo trực tiếp tại bể lưu trữ nhiên liệu, giữ nguyên vẹn 100% kết cấu thanh nhiên liệu, rút ngắn thời gian phân tích từ vài tuần xuống còn vài giờ và tiết kiệm hàng chục nghìn USD chi phí xử lý chất thải phóng xạ hóa học. Toàn bộ thuật toán tính toán phân rã phóng xạ, tích lũy đồng vị và giải hệ phương trình vi phân biến thiên số hạt nhân được lập trình tự động hóa trên nền tảng phần mềm MATLAB.
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 đã làm rõ quy luật tích lũy phi tuyến của đồng vị Cs-134 và quy luật tích lũy tuyến tính của Cs-137 theo thông lượng neutron và thời gian chiếu xạ. Cụ thể, Cs-137 có chu kỳ bán rã dài 30.17 năm và hiệu suất phân hạch khoảng 6.2%, đóng vai trò là thước đo tỷ lệ thuận tuyệt đối với tổng số phản ứng phân hạch; trong khi Cs-134 sinh ra từ phản ứng bắt neutron của Cs-133 bền mang tính chất phụ thuộc bậc hai vào thông lượng. Do đó, tỷ số hoạt độ Cs-134/Cs-137 thể hiện mối tương quan tuyến tính rất chặt chẽ với độ cháy của bó nhiên liệu số 62.
Thứ hai, việc xây dựng đường cong hiệu suất ghi tương đối của detector HPGe trong dải năng lượng gamma từ 604.7 keV đến 1408 keV đã triệt tiêu hoàn toàn sự phụ thuộc vào hình học đo và hiện tượng tự hấp thụ bức xạ trong vật liệu. Điều này giúp nâng cao độ chính xác phép đo với sai số thống kê thực nghiệm được kiểm soát dưới 5%.
Thứ ba, kết quả xác định phân bố hoạt độ của Cs-134 và Cs-137 dọc theo chiều dài bó nhiên liệu số 62 đã phác họa chính xác biên dạng phân bố độ cháy dạng cosine lệch. Giá trị độ cháy cực đại ghi nhận tại khu vực trung tâm thanh nhiên liệu cao hơn khoảng 25% đến 30% so với hai đầu mút của bó nhiên liệu, hoàn toàn phù hợp với lý thuyết phân bố trường thông lượng neutron nhiệt trong tâm lò phản ứng nghiên cứu Đà Lạt.
Thảo luận kết quả
Cơ chế vật lý giải thích cho sự thành công của phương pháp nằm ở đặc tính hóa lý của các sản phẩm phân hạch. Khác với các nguyên tố khí phóng xạ như Xenon và Krypton dễ bị khuếch tán ra khoảng hở giữa viên nhiên liệu và vỏ bọc dưới nhiệt độ cao, đồng vị Cesium tương đối ổn định trong mạng tinh thể nhiên liệu ở mức nhiệt độ vận hành của lò phản ứng nghiên cứu (dưới 1000 độ C).
Các dữ liệu thực nghiệm thu được từ luận văn có thể được trình bày và trực quan hóa vô cùng hiệu quả thông qua hệ thống đồ thị 2D biểu diễn đường cong suy giảm công suất theo năng lượng tia gamma, đồ thị hàm số biến thiên số hạt nhân theo thời gian chiếu xạ, cùng các bảng tổng hợp so sánh giá trị tính toán độ cháy trung bình với số liệu vận hành thực tế tại Lò phản ứng Đà Lạt. Kết quả so sánh giữa độ cháy tính toán bằng tỷ số đồng vị Cs-134/Cs-137 và giá trị tham chiếu từ lịch sử chạy lò cho thấy mức độ sai lệch chỉ dao động trong khoảng từ 3% đến 6%, khẳng định tính khả thi vượt trội của phương pháp trong công tác giám sát an toàn hạt nhân định kỳ.
Đề xuất và khuyến nghị
Nhằm phát huy tối đa giá trị nghiên cứu và đảm bảo an toàn tuyệt đối cho chu trình nhiên liệu, luận văn đưa ra 4 nhóm giải pháp cụ thể:
- Chuẩn hóa quy trình đo phổ gamma không phá hủy định kỳ: Ban Quản lý Lò phản ứng hạt nhân Đà Lạt cần ban hành quy chuẩn kỹ thuật đo kiểm tra độ cháy cho 100% các bó nhiên liệu VVR-M2 sau mỗi chu kỳ đảo nhiên liệu từ 6 đến 12 tháng, đảm bảo độ cháy không vượt quá ngưỡng an toàn cho phép.
- Nâng cấp và mở rộng thuật toán mô phỏng trên MATLAB: Nhóm nghiên cứu kỹ thuật lò cần tích hợp thêm mô hình hiệu chỉnh độ cháy cho các đồng vị phụ trợ như Europium-154 (với đỉnh năng lượng 1274 keV) và Samarium-149 nhằm giảm sai số tính toán xuống dưới 3% trong lộ trình 18 tháng tới.
- Xây dựng cơ sở dữ liệu số hóa theo dõi lịch sử nhiên liệu: Cục An toàn Bức xạ và Hạt nhân phối hợp cùng các viện nghiên cứu thiết lập hệ thống lưu trữ số hóa thông số vận hành, tốc độ sinh nhiệt tuyến tính LHGR và tỷ số đồng vị của từng bó nhiên liệu từ nay đến năm 2030.
- Ứng dụng kiểm tra sự toàn vẹn của vỏ bọc nhiên liệu Zirconium: Đội ngũ kỹ sư an toàn bức xạ cần kết hợp kết quả đo độ cháy với các phép kiểm tra độ ăn mòn hydrua hóa và độ dày lớp cặn oxit CRUD nhằm ngăn chặn sớm nguy cơ rão vỏ bọc do quá áp khi độ cháy đạt mức cao.
Đối tượng nên tham khảo luận văn
- Kỹ sư vận hành lò phản ứng hạt nhân: Nắm vững phương pháp đo phổ gamma trực tiếp tại chỗ để xác định độ cháy thực tế của từng thanh nhiên liệu, từ đó lập kế hoạch đảo nhiên liệu và tối ưu hóa cấu hình tâm lò nhằm duy trì tỷ số DNBR an toàn.
- Nhà nghiên cứu và giảng viên chuyên ngành Vật lý Hạt nhân: Khai thác kho dữ liệu khoa học gồm 44 hình vẽ, đồ thị chi tiết và 16 bảng thông số phân hạch để phục vụ công tác giảng dạy lý thuyết trường neutron và tương tác bức xạ với vật chất.
- Học viên cao học và nghiên cứu sinh ngành Năng lượng Nguyên tử: Tham khảo phương pháp xây dựng thuật toán trên MATLAB, kỹ thuật hiệu chuẩn phổ bán dẫn HPGe và cách giải bài toán phân rã dây chuyền hạt nhân cho các đề tài nghiên cứu chuyên sâu.
- Chuyên viên an toàn và thanh sát hạt nhân: Áp dụng phương pháp phân tích không phá hủy NDA để kiểm chứng khối lượng vật liệu phân hạch Uranium-235 và Plutonium-239 còn lại trong thanh nhiên liệu, phục vụ công tác thanh sát theo tiêu chuẩn của Cơ quan Năng lượng Nguyên tử Quốc tế.
Câu hỏi thường gặp
Độ cháy nhiên liệu hạt nhân là gì và tại sao cần xác định chính xác? Độ cháy là tổng năng lượng nhiệt giải phóng trên một đơn vị khối lượng vật liệu phân hạch ban đầu, tính bằng MWd/kg hoặc GWd/t. Xác định chính xác độ cháy giúp đánh giá mức độ suy giảm nhiên liệu, ngăn chặn hiện tượng quá nhiệt làm nóng chảy vỏ bọc và tối ưu hóa chi phí thay thế thanh nhiên liệu trong tâm lò.
Tại sao tỷ số đồng vị Cs-134/Cs-137 lại được ưu tiên lựa chọn để đo độ cháy? Cs-137 phân rã với chu kỳ 30.17 năm phản ánh tổng số phân hạch, trong khi Cs-134 sinh ra do bắt neutron phản ánh thông lượng neutron tích lũy. Tỷ số giữa hai đồng vị này có quan hệ tuyến tính với độ cháy, phát xạ các tia gamma năng lượng cao dễ ghi nhận và giúp loại bỏ hoàn toàn các sai số do hình học đo.
Phương pháp đo không phá hủy bằng phổ gamma có ưu điểm gì so với phân tích hóa học? Phương pháp phổ gamma giữ nguyên vẹn 100% hình thái thanh nhiên liệu, không tạo ra chất thải lỏng phóng xạ nguy hại, có thể tiến hành trực tiếp dưới nước tại bể lò phản ứng và tiết kiệm hàng trăm giờ thao tác phức tạp trong phòng thí nghiệm hóa phóng xạ kín.
Giới hạn độ cháy an toàn của nhiên liệu hạt nhân trên thế giới hiện nay là bao nhiêu? Tùy thuộc vào thiết kế lò, giới hạn độ cháy trung bình thường từ 40 GWd/t đến 65 GWd/t. Ví dụ, tại Pháp nhiên liệu UO2 giới hạn ở mức 52 GWd/t, nhiên liệu MOX là 47 GWd/t; tại Hoa Kỳ là 62 GWd/t; và tại Đức cho lò PWR lên tới 65 GWd/t.
Detector bán dẫn HPGe đóng vai trò như thế nào trong phép đo tỷ số đồng vị? Detector HPGe sở hữu độ phân giải năng lượng cực cao (thường dưới 2.0 keV tại đỉnh 1332 keV của Co-60), cho phép tách bạch rõ ràng các đỉnh năng lượng gamma nằm sát nhau của Cs-134 (604.7 keV, 795.8 keV) và Cs-137 (661.6 keV), đảm bảo độ chính xác định lượng cao nhất.
Kết luận
- Luận văn đã chứng minh tính ưu việt và độ chính xác cao của phương pháp tỷ số đồng vị Cs-134/Cs-137 trong việc xác định độ cháy nhiên liệu hạt nhân bằng kỹ thuật không phá hủy mẫu.
- Xây dựng thành công quy trình thực nghiệm đo phổ gamma bán dẫn HPGe kết hợp thuật toán tính toán tự động hóa trên nền tảng phần mềm MATLAB.
- Xác định đầy đủ giá trị độ cháy trung bình và phân bố độ cháy theo chiều cao hình học của bó nhiên liệu số 62 loại VVR-M2 tại Lò phản ứng hạt nhân Đà Lạt.
- Cung cấp hệ thống tư liệu khoa học chuẩn mực với 16 bảng số liệu tra cứu và 44 đồ thị đặc trưng phục vụ công tác nghiên cứu vật lý lò phản ứng.
- Đóng góp luận cứ khoa học quan trọng cho việc nâng cao tiêu chuẩn an toàn bức xạ, kiểm soát khủng hoảng sôi và tối ưu hóa chu trình kinh tế nhiên liệu tại Việt Nam.
Trong giai đoạn 12 đến 36 tháng tới, các cơ sở nghiên cứu hạt nhân cần tiếp tục chuẩn hóa và mở rộng kỹ thuật này cho toàn bộ các cấu hình bó nhiên liệu làm giàu thấp mới. Hãy liên hệ và hợp tác với các nhóm nghiên cứu chuyên ngành Vật lý Hạt nhân để tiếp nhận chuyển giao thuật toán, ứng dụng hiệu quả công nghệ đo đạc tiên tiến vào việc quản lý an toàn lò phản ứng hạt nhân.