Tổng quan nghiên cứu

Trong bối cảnh nguồn nhiên liệu hóa thạch truyền thống như than đá và dầu mỏ đang suy giảm nhanh chóng, năng lượng hạt nhân đã trở thành giải pháp công nghệ quan trọng nhằm đảm bảo an ninh năng lượng toàn cầu. Tại Việt Nam, định hướng phát triển điện hạt nhân đã được xác lập cụ thể thông qua chủ trương đầu tư của Quốc hội vào ngày 25/11/2009 và hiệp định hợp tác quốc tế ký kết ngày 31/10/2010. Do tính chất công suất lớn của các tổ máy thương mại đạt từ 1000 MW đến 1200 MW, yêu cầu an toàn thủy nhiệt trong mọi điều kiện vận hành và sự cố luôn được đặt lên hàng đầu.

Vấn đề nghiên cứu trọng tâm của luận văn là phân tích an toàn thủy nhiệt của lò phản ứng nước áp lực (PWR) khi xảy ra sự cố mất nước tải nhiệt vết nứt nhỏ (SB-LOCA với diện tích nứt nhỏ hơn hoặc bằng 465 cm²). Khác với sự cố vết nứt lớn (LB-LOCA) vốn kích hoạt ngay lập tức toàn bộ hệ thống làm mát khẩn cấp, sự cố SB-LOCA có tiến trình suy giảm áp suất chậm, tiềm ẩn nguy cơ làm trễ sự can thiệp của hệ thống an toàn và dẫn đến phơi trần vùng hoạt.

Mục tiêu cụ thể của nghiên cứu là xây dựng mô hình mô phỏng thủy nhiệt chi tiết cho lò phản ứng nước áp lực Zion (Mỹ) bằng chương trình tính toán RELAP5/Mod3.3. Luận văn tiến hành khảo sát các hiện tượng chuyển tiếp thủy động phức tạp, đánh giá hiệu quả của hệ thống làm mát vùng hoạt khẩn cấp (ECCS) và đối chiếu với các tiêu chuẩn an toàn nghiêm ngặt của Ủy ban Pháp quy Hạt nhân Hoa Kỳ (US NRC). Kết quả nghiên cứu cung cấp cơ sở dữ liệu quan trọng giúp kiểm soát nhiệt độ vỏ thanh nhiên liệu dưới ngưỡng 1204°C và khống chế mức độ oxy hóa kim loại dưới 17%, đóng góp thiết thực cho quá trình chuẩn bị nhân lực và tiếp nhận công nghệ điện hạt nhân công suất lớn.

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 tổ hợp các lý thuyết thủy nhiệt và vật lý lò phản ứng hạt nhân hiện đại nhằm mô tả chính xác trạng thái chuyển tiếp của chất làm mát:

  • Mô hình thủy động hỗn hợp hai pha phi cân bằng (Two-Fluid Hydrodynamic Model): Sử dụng hệ 6 phương trình vi phân đạo hàm riêng để bảo toàn khối lượng, động lượng và năng lượng độc lập cho từng pha lỏng và pha hơi. Mô hình tích hợp phương trình vận chuyển nồng độ Boron hòa tan và phương trình trạng thái cho khí không ngưng tụ (như Nitơ hoặc Hydro sinh ra từ phản ứng Zirconium - nước).
  • Mô hình dẫn nhiệt một chiều không dừng: Ứng dụng phương trình Fourier truyền nhiệt qua cấu trúc chất rắn nhiều lớp (vỏ thùng lò, thành ống sinh hơi và các thanh nhiên liệu Zircaloy) kết hợp với các điều kiện biên truyền nhiệt đối lưu tự nhiên, sôi bọt và sôi màng.
  • Mô hình động học lò điểm (Point Kinetics): Tính toán biến thiên mật độ nơtron và công suất phân rã phóng xạ sau khi dập lò thông qua 6 nhóm nơtron trễ, với thời gian sống nơtron từ 10⁻⁸ đến 10⁻³ giây. Nhiệt phân rã được xác định chiếm khoảng 5% đến 6% công suất định mức trong giai đoạn đầu sau dập lò.
  • Các khái niệm cốt lõi: Luận văn làm rõ cơ chế của hệ thống làm mát vùng hoạt khẩn cấp ECCS (bao gồm bơm cao áp HPSI, bình tích áp ACC dung tích 28.3 - 42.5 m³ và bơm hạ áp LPSI lưu lượng 190 kg/s), hiện tượng giới hạn dòng thuận nghịch (CCFL) và hiệu ứng nghẽn thủy lực (Loop Seal Clearing).

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

Nghiên cứu sử dụng phương pháp mô phỏng số thông qua phần mềm tiêu chuẩn công nghiệp RELAP5/Mod3.3 do Phòng thí nghiệm Kỹ thuật Quốc gia Idaho (INEEL) và US NRC phát triển. Lý do lựa chọn công cụ này là khả năng mô phỏng chính xác các quá trình nhiệt thủy động hai pha một chiều và các chế độ làm mát đối lưu tự nhiên trong trạng thái chuyển tiếp khắc nghiệt.

  • Cỡ mẫu và đối tượng mô phỏng: Mô hình hóa toàn diện nhà máy điện hạt nhân Zion với sơ đồ nút mạng thủy lực gồm hơn 100 thể tích điều khiển (volumes), các nhánh nối (junctions) và các cấu trúc nhiệt (heat structures). Cỡ mẫu khảo sát bao gồm 5 kịch bản kích thước vết nứt từ 1% đến 4% diện tích tiết diện ống dẫn (tương đương đường kính lỗ nứt từ 2 inch đến 4 inch) tại các vị trí kênh nóng và kênh lạnh.
  • Phương pháp chọn mẫu: Lựa chọn kịch bản dựa trên phương pháp phân tích bất định và đánh giá rủi ro định lượng đối với các vị trí dễ bị tổn thương nhất trong hệ thống áp lực 15.5 MPa.
  • Tiến trình nghiên cứu: Quá trình tính toán được thực hiện theo 2 bước: thiết lập trạng thái dừng ổn định (steady-state) trong 100 giây đầu để chuẩn hóa các thông số vận hành (áp suất 15.5 MPa, nhiệt độ nước đầu vào khoảng 560 K), tiếp theo là mô phỏng trạng thái chuyển tiếp sự cố SB-LOCA trong khoảng thời gian từ 0 đến 3600 giây.

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

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

Quá trình mô phỏng số trên phần mềm RELAP5/Mod3.3 đã mang lại các phát hiện khoa học quan trọng về động học sự cố SB-LOCA trên lò PWR:

  1. Quy luật suy giảm áp suất theo kích thước vết nứt: Tốc độ giảm áp suất tỷ lệ thuận với diện tích vết nứt. Với vết nứt 2% (đường kính khoảng 4 inch), áp suất vòng sơ cấp giảm từ 15.5 MPa xuống dưới 4.14 MPa sau khoảng 450 giây, kích hoạt bình tích áp ACC xả nước làm mát khẩn cấp. Đối với vết nứt 3%, thời gian kích hoạt ACC rút ngắn còn khoảng 280 giây, giúp hệ thống phục hồi mức nước nhanh hơn 35% so với kịch bản vết nứt 1%.
  2. Hiện tượng phơi trần vùng hoạt do Loop Seal Clearing: Tại thời điểm khoảng 300 đến 600 giây của quá trình chuyển tiếp, sự chênh lệch áp suất thủy tĩnh giữa các nhánh của ống chữ U trong bình sinh hơi đẩy chất lỏng ra khỏi vùng hoạt. Mức nước hai pha trong thùng lò tụt giảm cục bộ xuống còn khoảng 30% đến 40% chiều cao vùng hoạt, khiến nhiệt độ đỉnh vỏ thanh nhiên liệu (PCT) tăng vọt từ 320°C lên mức cực đại khoảng 650°C đến 850°C.
  3. Ảnh hưởng của vị trí vết nứt và dòng Bypass: Vết nứt xuất hiện tại kênh lạnh (Cold Leg) gây thất thoát lưu lượng nước lớn hơn so với kênh nóng do hiện tượng dòng nước ECCS đi tắt ra ngoài vết nứt, làm giảm khoảng 20% lượng nước nạp thực tế vào khoang lưu hồi thùng lò.
  4. Tác động của độ trễ hệ thống cấp nước phụ (AFW): Khi hệ thống AFW khởi động trễ trên 15 giây, khả năng tản nhiệt qua bình sinh hơi thứ cấp suy giảm nghiêm trọng, duy trì áp suất vòng sơ cấp ở mức cao trên 8.7 MPa lâu hơn, làm chậm thời điểm xả của bình ACC thêm khoảng 120 giây.

Thảo luận kết quả

Các kết quả tính toán có thể được hệ thống hóa trực quan thông qua các biểu đồ biến thiên áp suất theo thời gian (đường cong P-t), biểu đồ mức nước thùng lò và bảng đối chiếu nhiệt độ đỉnh vỏ thanh nhiên liệu. Dữ liệu từ mô phỏng chỉ ra rằng, dù vùng hoạt bị phơi trần tạm thời do hiệu ứng Loop Seal Clearing, nhiệt độ cực đại của vỏ thanh nhiên liệu trong mọi kịch bản khảo sát (dao động từ 650°C đến 850°C) vẫn thấp hơn đáng kể so với giới hạn an toàn 1204°C do US NRC quy định. Mức độ oxy hóa Zirconium bề mặt được duy trì dưới 3%, thấp hơn nhiều so với ngưỡng giới hạn 17%.

So sánh với các báo cáo thực nghiệm quốc tế tại các cơ sở thí nghiệm bán quy mô như LOFT và ROSA-IV, mô hình tính toán bằng RELAP5 thể hiện độ tin cậy cao khi tái hiện chính xác 3 chế độ truyền nhiệt tuần hoàn tự nhiên: đối lưu đơn pha, dòng hai pha và ngưng tụ hồi lưu (reflux condensation). Sự hiện diện của khí không ngưng tụ (Nitơ từ bình ACC khi xả cạn) được xác định là nguyên nhân làm suy giảm khoảng 15% hệ số truyền nhiệt ngưng tụ trong các ống truyền nhiệt chữ U của bình sinh hơi.

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

Dựa trên các kết quả phân tích thủy nhiệt định lượng, nghiên cứu đưa ra 4 nhóm khuyến nghị kỹ thuật cụ thể nhằm nâng cao mức độ an toàn cho các tổ máy PWR:

  1. Chuẩn hóa quy trình vận hành khẩn cấp (EOP): Cơ quan quản lý an toàn hạt nhân và đơn vị vận hành nhà máy cần xây dựng kịch bản ứng phó sự cố SB-LOCA chuyên biệt. Cần quy định thời gian can thiệp của nhân viên vận hành trong khoảng 10 đến 15 phút đầu tiên, ưu tiên duy trì liên tục hệ thống cấp nước phụ AFW với lưu lượng tối thiểu 190 kg/s để hỗ trợ tản nhiệt thứ cấp.
  2. Tối ưu hóa thông số thiết kế hệ thống bình tích áp (ACC): Doanh nghiệp thiết kế và vận hành cần duy trì áp suất nạp khí Nitơ trong các bình ACC ở ngưỡng ổn định từ 2.0 MPa đến 4.14 MPa, đồng thời dung tích nước dự trữ phải đạt tối thiểu 35 m³ mỗi bình để đảm bảo khả năng ngập lại vùng hoạt ngay sau khi áp suất vòng sơ cấp suy giảm.
  3. Nâng cấp hệ thống giám sát mực nước thùng lò: Đơn vị công nghệ cần lắp đặt bổ sung các cảm biến đo mức hai pha độc lập trong thùng lò và khoang lưu hồi. Việc này giúp loại bỏ hoàn toàn hiện tượng chỉ thị sai lệch mức nước do ứ đọng dòng hơi tại ống nối bình điều áp (với lưu lượng hơi từ 2.27 đến 3.63 kg/s), ngăn chặn quyết định sai lầm của người vận hành như từng xảy ra trong sự cố Three Mile Island.
  4. Định kỳ diễn tập và tái đánh giá an toàn thủy nhiệt: Thực hiện đánh giá định kỳ 12 tháng một lần đối với toàn bộ hệ thống ECCS bằng các phần mềm mô phỏng nâng cao. Ban quản lý an toàn cần đưa các kịch bản hỏng hóc kép (như mất điện toàn nhà máy kết hợp vỡ ống 2% đến 4%) vào chương trình đào tạo trên buồng điều khiển mô phỏng toàn phần (Full-scope Simulator).

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

Luận văn là tài liệu chuyên khảo chuyên sâu, mang lại giá trị ứng dụng thiết thực cho 4 nhóm đối tượng:

  • Kỹ sư vận hành và chuyên viên an toàn nhà máy điện hạt nhân: Tiếp cận dữ liệu phân tích chi tiết về các pha chuyển tiếp nhiệt động để xây dựng hướng dẫn quản lý tai nạn nghiêm trọng (SAMG) và tối ưu hóa thao tác điều khiển ECCS trong điều kiện sự cố mất nước làm mát.
  • Cán bộ quản lý và chuyên viên thẩm định pháp quy hạt nhân: Sử dụng phương pháp luận và các kết quả định lượng trong luận văn làm cơ sở thẩm định Báo cáo Phân tích An toàn (SAR), đánh giá tính tương thích của hệ thống an toàn thụ động theo các tiêu chuẩn IAEA và US NRC.
  • Giảng viên, nghiên cứu sinh và học viên cao học ngành Năng lượng nguyên tử: Khảo sát khung lý thuyết thủy nhiệt hai pha, mô hình động học lò điểm và phương pháp lập trình mô phỏng bằng mã RELAP5/Mod3.3 phục vụ công tác giảng dạy và nghiên cứu mở rộng.
  • Kỹ sư thiết kế và tư vấn công nghệ năng lượng: Tham khảo các thông số thủy động học, cấu trúc bó nhiên liệu và cơ chế giải trừ nhiệt phân rã (5% công suất) để thiết kế tối ưu hệ thống trao đổi nhiệt và bình sinh hơi cho các lò phản ứng thế hệ mới.

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

Sự cố SB-LOCA là gì và khác biệt như thế nào so với LB-LOCA?

Sự cố SB-LOCA là tai nạn mất nước tải nhiệt do vết nứt có diện tích nhỏ hơn hoặc bằng 465 cm² trên đường ống áp lực của vòng sơ cấp. Trong khi sự cố vết nứt lớn (LB-LOCA) làm sụt áp tức thì và kích hoạt toàn bộ hệ thống an toàn trong vài giây, SB-LOCA diễn biến chậm hơn, khiến áp suất giảm từ từ và đòi hỏi sự phối hợp chính xác giữa các cấp bơm HPSI, ACC và LPSI.

Tại sao sự cố SB-LOCA lại tiềm ẩn rủi ro nóng chảy thanh nhiên liệu cao?

Do áp suất vòng sơ cấp giảm chậm, hệ thống bơm áp suất thấp và bình ACC không thể tự động xả nước ngay lập tức. Nếu hiện tượng Loop Seal Clearing hoặc Core Boil-off xuất hiện, nước trong vùng hoạt bị đẩy ra ngoài hoặc hóa hơi nhanh chóng, làm phơi trần các thanh nhiên liệu và khiến nhiệt độ vỏ Zircaloy tăng vượt 800°C nếu không được cấp nước kịp thời.

Chương trình tính toán thủy nhiệt RELAP5 đóng vai trò gì trong an toàn hạt nhân?

RELAP5 là bộ mã mô phỏng nhiệt thủy động tiêu chuẩn quốc tế, cho phép giải hệ phương trình bảo toàn hai pha phi cân bằng nhiệt và phi cân bằng vận tốc. Phần mềm mô phỏng chính xác các thông số áp suất, lưu lượng, nhiệt độ vỏ nhiên liệu và mật độ Boron theo thời gian thực, phục vụ công tác cấp phép và phân tích an toàn nhà máy điện hạt nhân.

Hiện tượng Loop Seal Clearing ảnh hưởng thế nào đến vùng hoạt?

Loop Seal Clearing xảy ra khi khối lượng hỗn hợp hơi - nước tích tụ tại đoạn ống chữ U giữa bình sinh hơi và bơm sơ cấp tạo ra chênh lệch áp suất thủy tĩnh lớn. Áp lực này đẩy ngược mức nước hai pha trong vùng hoạt tụt xuống chân thanh nhiên liệu, làm suy giảm đột ngột khả năng làm mát và tạo ra đỉnh nhiệt độ vỏ nhiên liệu tạm thời.

Các tiêu chuẩn an toàn nào của US NRC được áp dụng để đánh giá chấp nhận sự cố?

US NRC quy định 3 tiêu chuẩn cốt lõi: đỉnh nhiệt độ của lớp vỏ nhiên liệu (PCT) không được vượt quá 1204°C (2200°F), tổng chiều dày lớp oxy hóa Zirconium không vượt quá 17% độ dày vỏ ban đầu, và tổng lượng khí Hydro sinh ra từ phản ứng oxy hóa không vượt quá 1% lượng hydro tối đa giả định.

Kết luận

  • Luận văn đã mô hình hóa thành công các quá trình nhiệt thủy động chuyển tiếp phức tạp trong sự cố mất nước tải nhiệt vết nứt nhỏ (SB-LOCA) trên lò phản ứng nước áp lực PWR bằng chương trình tiêu chuẩn RELAP5/Mod3.3.
  • Kết quả mô phỏng trên lò Zion khẳng định đỉnh nhiệt độ vỏ thanh nhiên liệu trong các kịch bản nứt từ 1% đến 4% đạt mức cực đại khoảng 650°C đến 850°C, hoàn toàn nằm trong giới hạn an toàn 1204°C của US NRC.
  • Làm rõ cơ chế vật lý của hiện tượng phơi trần vùng hoạt do hiệu ứng Loop Seal Clearing và ảnh hưởng bất lợi của độ trễ hệ thống cấp nước phụ AFW (vượt quá 15 giây).
  • Kế hoạch tiếp theo cần tập trung mở rộng mô phỏng không gian 3 chiều đối với khoang lưu hồi và tích hợp các mô hình phân tích bất định (Uncertainty Analysis) cho các thế hệ lò phản ứng PWR tiên tiến (VVER-1200, AP1000).
  • Công trình là nguồn tài liệu học thuật và kỹ thuật giá trị cao; các tổ chức, kỹ sư và nhà nghiên cứu trong ngành năng lượng hạt nhân nên khai thác chi tiết các mô hình số và khuyến nghị vận hành của luận văn để nâng cao năng lực phân tích an toàn công nghệ.