Tổng quan nghiên cứu

Dự án phát triển nguồn điện hạt nhân với công suất dự kiến 2.000 MW tại xã Phước Dinh, huyện Thuận Nam, tỉnh Ninh Thuận từng được hoạch định với tổng mức đầu tư ước tính khoảng 200.000 tỷ đồng nhằm đáp ứng nhu cầu năng lượng quốc gia. Tuy nhiên, an toàn bức xạ và kiểm soát rủi ro môi trường luôn là điều kiện tiên quyết hàng đầu trong toàn bộ vòng đời của cơ sở hạt nhân. Khi xảy ra các kịch bản sự cố nghiêm trọng giả định, các sản phẩm phân hạch dạng khí và hạt nhân phóng xạ như các đồng vị phóng xạ iốt sẽ phát tán vào tầng đối lưu, gây nguy cơ chiếu xạ toàn thân và chiếu xạ tuyến giáp đối với cộng đồng dân cư xung quanh.

Nghiên cứu này tập trung vào bài toán đánh giá định lượng khả năng vận chuyển và khuếch tán chất phóng xạ trong môi trường khí quyển tại khu vực dự kiến xây dựng Nhà máy điện hạt nhân Ninh Thuận 1. Mục tiêu cốt lõi là thiết lập mô hình tính toán nồng độ phóng xạ pha loãng tương đối trong không khí, từ đó xác lập cơ sở khoa học chính xác cho việc phân định phạm vi vùng cấm dân cư và vùng hạn chế dân cư theo các quy chuẩn an toàn hạt nhân quốc gia và quốc tế.

Phạm vi nghiên cứu bao quát không gian địa lý thuộc thôn Vĩnh Trường, xã Phước Dinh trong bán kính 20 km với tổng số 310.561 người cư trú, tại tọa độ 11°25’54” đến 11°27’17” vĩ độ Bắc và 108°59’43” đến 109°01’00” kinh độ Đông. Dữ liệu khí tượng thu thập trong giai đoạn 2 năm liên tiếp (2008 - 2009) phục vụ mô phỏng động lực học khí quyển. Kết quả nghiên cứu cung cấp các chỉ số khoảng cách an toàn cụ thể theo 16 hướng gió, đóng góp trực tiếp vào công tác quy hoạch mặt bằng nhà máy, quản lý sử dụng đất đai và xây dựng phương án ứng phó sự cố bức xạ khẩn cấp cho địa phương.

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 hệ thống lý thuyết lan truyền chất ô nhiễm trong khí quyển dựa trên định luật khuếch tán Fick, lý thuyết vận chuyển gradient tại điểm cố định theo hệ quy chiếu Euler và lý thuyết thống kê nhiễu loạn Lagrange. Mô hình luồng khí Gauss phản xạ toàn phần trên bề mặt đất được lựa chọn làm công cụ tính toán nồng độ chất phóng xạ theo phương trình phân bố nồng độ ba chiều trong không gian:

$$\chi(x, y, z, H) = \frac{Q}{2\pi u \sigma_y \sigma_z} \exp\left(-\frac{y^2}{2\sigma_y^2}\right) \left[ \exp\left(-\frac{(z-H)^2}{2\sigma_z^2}\right) + \exp\left(-\frac{(z+H)^2}{2\sigma_z^2}\right) \right]$$

Trong đó, $\chi$ là nồng độ phóng xạ trong không khí (Bq/m³), $Q$ là tốc độ phát thải nguồn điểm (Bq/s), $u$ là vận tốc gió trung bình ở độ cao hiệu dụng (m/s), $H$ là độ cao hiệu dụng của ống khói (m), $\sigma_y$ và $\sigma_z$ là các hệ số phát tán theo phương ngang và phương thẳng đứng (m).

Trạng thái xáo trộn của khí quyển lớp biên được chuẩn hóa qua 6 cấp độ ổn định Pasquill-Gifford từ loại A (rất không ổn định) đến loại F (rất ổn định) và mở rộng cấp G (cực kỳ ổn định). Các hệ số khuếch tán $\sigma_y$ và $\sigma_z$ được tính toán thông qua quan hệ hàm mũ kinh nghiệm của Caraway và Klug tương ứng với từng dải khoảng cách từ 100 m đến trên 10.000 m. Sự biến thiên của tốc độ gió theo chiều cao lớp biên khí quyển từ mặt đất 10 m lên độ cao phát thải được hiệu chỉnh thông qua quy luật hàm mũ với chỉ số điều chỉnh $p$ dao động từ 0,07 đến 0,55 tùy theo cấp độ ổn định địa hình nông thôn hoặc ven biển.

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

Nguồn dữ liệu khí tượng phục vụ nghiên cứu được trích xuất từ Hệ thống xử lý dữ liệu toàn cầu GDAS thuộc Trung tâm Tài nguyên Khí quyển (ARL) của Cơ quan Quản lý Khí quyển và Đại dương Quốc gia Hoa Kỳ (NOAA). Bộ dữ liệu có độ phân giải cao bao gồm chuỗi quan trắc liên tục với cỡ mẫu chuẩn 2.920 mốc thời gian mỗi năm (chu kỳ quan trắc 3 giờ một lần trong suốt 365 ngày của hai năm 2008 và 2009).

Phương pháp chọn mẫu áp dụng là thu thập toàn bộ chuỗi chu kỳ thời gian đồng nhất tại tọa độ thực tế của nhà máy để loại trừ sai số mùa vụ. Lý do lựa chọn nguồn dữ liệu NOAA GDAS là khả năng đồng hóa dữ liệu vệ tinh kết hợp trạm mặt đất, phản ánh chính xác các thông số tốc độ gió, hướng gió, bức xạ nhiệt và độ che phủ mây tại khu vực ven biển nhiệt đới vốn thiếu trạm quan trắc cao không chuyên sâu.

Phương pháp xử lý và phân tích số liệu kết hợp phần mềm chuyên dụng WRPLOT View để lập biểu đồ hoa gió 16 hướng và phần mềm Google Earth để chồng lớp bản đồ không gian. Để xác định ranh giới an toàn bức xạ theo hướng dẫn kỹ thuật số 1.145 của Cơ quan Pháp quy Hạt nhân Hoa Kỳ (US NRC) và Thông tư 13/2009/TT-BKHCN, nghiên cứu áp dụng thuật toán thống kê phân bố tần suất tích lũy. Khoảng cách an toàn được xác định tại ngưỡng xác suất vượt quá 0,5% trên tổng số giờ quan sát trong năm, kết hợp phương pháp hồi quy tuyến tính với hệ số xác định $R^2$ đạt trên 0,95 để nội suy khoảng cách chính xác cho từng hướng gió.

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

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

Phân tích dữ liệu khí tượng 2 năm 2008 - 2009 tại khu vực Phước Dinh cho thấy trường gió mang tính phân hóa theo mùa sâu sắc nhưng có tính ổn định chu kỳ qua các năm:

  • Hướng gió thịnh hành tập trung chủ yếu ở hai trục chính: hướng Đông Bắc chiếm tỷ trọng trên 30% tổng thời gian quan sát với cường độ gió mạnh, và hướng Tây Nam chiếm trên 15% tổng thời gian quan sát với cường độ trung bình. Dải tốc độ gió từ 3,6 m/s đến 7,7 m/s (tương đương 7 đến 17 Knots) chiếm tỷ lệ áp đảo hơn 60% tổng lượng gió cả năm.
  • Phân loại độ ổn định khí quyển ghi nhận lớp D (trạng thái trung hòa) chiếm tỷ lệ cao nhất với hơn 40% thời gian trong năm. Các trạng thái kém ổn định (loại A, B, C) xuất hiện chủ yếu vào ban ngày khi có bức xạ mặt trời mạnh, trong khi trạng thái ổn định (loại E, F) chiếm tỷ lệ thấp hơn vào ban đêm.
  • Với giới hạn liều chiếu xạ tuyến giáp do các đồng vị iốt phóng xạ (I-131, I-132, I-133, I-134, I-135) không vượt quá 3 Sv trong 2 giờ đầu sự cố, giá trị nồng độ tương đối cho phép tối đa được xác định là $3,09 \times 10^{-4}\text{ s/m}^3$. Bán kính vùng cấm dân cư tương ứng tại ngưỡng xác suất 0,5% dao động từ 800 m đến 1.239 m tùy theo từng hướng gió.
  • Đối với kịch bản phơi nhiễm dài hạn trong 30 ngày với nồng độ tương đối giới hạn là $5,55 \times 10^{-6}\text{ s/m}^3$, bán kính vùng hạn chế dân cư mở rộng từ 12.154 m đến 16.576 m, kéo dài rõ rệt nhất theo hai hướng Đông Bắc và Tây Nam.

Thảo luận kết quả

Đặc thù địa hình Ninh Thuận với 63,2% diện tích là đồi núi bao bọc ba mặt tạo nên hiệu ứng kênh dẫn gió ven biển, khiến luồng khí phát thải khi gặp sự cố có xu hướng phát tán dọc theo dải bờ biển hẹp thay vì tỏa đều tròn đồng tâm. Biểu đồ hoa gió 16 hướng phản ánh rõ nét sự bất đối xứng này.

Khi trình bày dữ liệu phát tán trên bản đồ số học không gian qua các lớp đẳng nồng độ nơron phóng xạ, vùng bảo vệ an toàn bức xạ thực tế có dạng elip kéo dài theo trục Đông Bắc - Tây Nam. Điều này lý giải tại sao cự ly an toàn ở hướng Đông Bắc (16.576 m) lại lớn hơn đáng kể so với các hướng Tây hoặc Tây Bắc (chỉ khoảng 7.000 m đến 9.000 m).

So sánh giữa hai hệ số thực nghiệm Caraway và Klug cho thấy hệ số Caraway cho kết quả phân hóa chi tiết và sát với thực tế địa hình ven biển hơn ở cự ly trên 5.000 m, trong khi hệ số Klug có xu hướng đánh giá đồng nhất hơn ở phạm vi gần. Sự phù hợp giữa mô hình luồng khí Gauss ngoại suy và số liệu khí tượng NOAA khẳng định độ tin cậy của phương pháp thống kê xác suất 0,5% trong việc bảo vệ an toàn cho 310.561 người dân trong bán kính 20 km quanh tâm lò phản ứng.

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

Dựa trên các kết quả tính toán định lượng về phạm vi phát tán phóng xạ, nghiên cứu đề xuất 4 nhóm giải pháp chiến lược nhằm tối ưu hóa an toàn bức xạ:

  1. Thiết lập vành đai kiểm soát quy hoạch dân cư bất đối xứng: Ủy ban nhân dân tỉnh Ninh Thuận phối hợp cùng Bộ Khoa học và Công nghệ ban hành quy chế sử dụng đất nghiêm ngặt. Xác lập vùng cấm dân cư cố định bán kính tối thiểu 1.200 m tính từ tâm lò phản ứng (đảm bảo mật độ dân số 0 người/km²) và vùng hạn chế dân cư mở rộng theo hình elip với cự ly 16,5 km về hướng Đông Bắc và Tây Nam trước khi triển khai thi công xây dựng dự án.
  2. Xây dựng trạm quan trắc khí tượng và đo bức xạ tự động đa tầng: Viện Năng lượng nguyên tử Việt Nam chủ trì lắp đặt tháp quan trắc khí tượng cao 100 m ngay tại địa điểm Phước Dinh với chu kỳ cập nhật dữ liệu tự động 10 phút/lần. Hệ thống phải tích hợp cảm biến đo trực tiếp građien nhiệt độ thẳng đứng và độ lệch chuẩn dao động gió ngang nhằm liên tục hiệu chỉnh hệ số Pasquill theo thời gian thực.
  3. Chuẩn hóa quy trình vận hành kịch bản ứng phó sự cố khẩn cấp: Cục An toàn bức xạ và hạt nhân xây dựng hệ thống phần mềm dự báo phát tán phóng xạ thời gian thực tích hợp GIS, kết nối dữ liệu trực tuyến với cơ quan khí tượng quốc gia. Tổ chức diễn tập sơ tán định kỳ 12 tháng/lần cho cộng đồng dân cư trong dải ảnh hưởng ưu tiên 15 km thuộc các xã Phước Dinh và Vĩnh Hải.
  4. Phát triển mô hình mô phỏng thủy động lực học chuyên sâu: Các cơ sở nghiên cứu và trường đại học tiếp tục triển khai các mô hình phát tán khí động học phi tĩnh (như CALPUFF hoặc WRF-Chem) để đánh giá chi tiết ảnh hưởng của vi khí hậu ven biển và tương tác gió đất - gió biển tại Ninh Thuận trong giai đoạn chuẩn bị dự án tiếp theo.

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

Nội dung và phương pháp luận của luận văn mang lại giá trị khoa học và ứng dụng thực tiễn cho 4 nhóm đối tượng chính:

  • Cơ quan quản lý nhà nước về an toàn bức xạ và hạt nhân: Các chuyên viên thẩm định tại Cục An toàn bức xạ và hạt nhân, Bộ Khoa học và Công nghệ có thể sử dụng bộ số liệu tính toán nồng độ tương đối và cự ly xác suất 0,5% làm tài liệu tham chiếu trong quá trình xây dựng các quy chuẩn kỹ thuật quốc gia về lựa chọn địa điểm nhà máy điện hạt nhân.
  • Kỹ sư quy hoạch và thiết kế nhà máy điện: Các chuyên gia thuộc Tập đoàn Điện lực Việt Nam và các đơn vị tư vấn năng lượng có thể ứng dụng trực tiếp thuật toán xử lý dữ liệu gió GDAS NOAA và mô hình luồng khí Gauss để quy hoạch bố trí tổng mặt bằng công trình, định vị ống khói phát thải và thiết kế hàng rào an ninh sinh học.
  • Giảng viên và nhà nghiên cứu học thuật: Các viện nghiên cứu vật lý hạt nhân, giảng viên và học viên cao học chuyên ngành Vật lý nguyên tử, Kỹ thuật hạt nhân, Khoa học môi trường có thể khai thác khung lý thuyết khuếch tán Fick, hệ số Caraway - Klug và quy trình xử lý dữ liệu hoa gió làm tài liệu giảng dạy và mở rộng đề tài nghiên cứu.
  • Ban chỉ huy phòng chống thiên tai và tìm kiếm cứu nạn cấp tỉnh: Lực lượng ứng phó tình huống khẩn cấp tại địa phương có thể sử dụng bản đồ phân vùng rủi ro 16 hướng gió để xây dựng các tuyến đường sơ tán dân cư tối ưu khi xảy ra sự cố rò rỉ hóa chất độc hại hoặc chất phóng xạ trên địa bàn ven biển.

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

1. Vì sao mô hình luồng khí Gauss được ưu tiên sử dụng trong đánh giá an toàn nhà máy điện hạt nhân?
Mô hình luồng khí Gauss được thừa nhận rộng rãi bởi Cơ quan Năng lượng Nguyên tử Quốc tế (IAEA) và US NRC nhờ nền tảng toán học giải tích chặt chẽ, tốc độ xử lý nhanh và độ chính xác phù hợp với dữ liệu thực nghiệm trong phạm vi bán kính phát tán dưới 20 km khi đánh giá các kịch bản phát thải liên tục hoặc tức thời.

2. Tiêu chuẩn xác định vùng cấm dân cư và vùng hạn chế dân cư dựa trên chỉ số liều bức xạ nào?
Theo Thông tư 13/2009/TT-BKHCN và tài liệu kỹ thuật TID-14844, vùng cấm dân cư được xác định sao cho liều chiếu toàn thân trong 2 giờ đầu sự cố không vượt quá 250 mSv và liều chiếu tuyến giáp do iốt phóng xạ không vượt quá 3 Sv (3.000 mSv). Vùng hạn chế dân cư áp dụng cùng định mức liều này cho toàn bộ thời gian diễn ra sự cố kéo dài 30 ngày.

3. Tại sao dữ liệu khí tượng NOAA GDAS lại phù hợp cho bài toán tính toán phát tán phóng xạ tại Ninh Thuận?
Hệ thống GDAS của NOAA cung cấp chuỗi dữ liệu khí tượng toàn cầu chuẩn hóa với mật độ quan trắc định kỳ 3 giờ/lần liên tục trong 365 ngày. Nguồn dữ liệu này phản ánh đầy đủ các thông số bức xạ mặt trời, nhiệt độ đối lưu và tốc độ gió đa tầng, bù đắp hiệu quả sự thiếu hụt trạm đo tháp cao tại khu vực bán sơn địa Phước Dinh.

4. Đặc điểm khí tượng tại Ninh Thuận ảnh hưởng như thế nào đến hình dạng vùng an toàn dân cư?
Do ảnh hưởng của chế độ gió mùa với hai hướng gió chủ đạo là Đông Bắc (chiếm trên 30%) và Tây Nam (chiếm trên 15%), luồng khí phát tán mạnh nhất theo trục này. Do đó, vùng bảo vệ an toàn không phải hình tròn đồng tâm mà có dạng hình elip kéo dài, với bán kính an toàn hướng Đông Bắc lớn gấp đôi các hướng khác.

5. Điểm khác biệt mấu chốt giữa hệ số phát tán Caraway và Klug trong nghiên cứu là gì?
Hệ số Caraway chia nhỏ các khoảng cách lan truyền thành nhiều đoạn cụ thể (dưới 500 m, 500 m đến 5.000 m và trên 5.000 m) với các cặp số mũ riêng biệt, giúp mô tả chính xác hơn sự suy giảm nồng độ ở cự ly xa trên 10 km so với công thức hệ số đơn của Klug.

Kết luận

  • Luận văn đã xây dựng thành công quy trình đánh giá định lượng khả năng phát tán phóng xạ trong môi trường khí quyển tại địa điểm xây dựng Nhà máy điện hạt nhân Ninh Thuận 1 bằng việc tích hợp mô hình luồng khí Gauss với chuỗi dữ liệu khí tượng NOAA GDAS 2008 - 2009 gồm 2.920 mốc quan trắc mỗi năm.
  • Xác định rõ nét đặc trưng khí tượng lớp biên tại xã Phước Dinh với độ ổn định khí quyển loại D chiếm ưu thế trên 40%, cùng chế độ gió phân hóa mạnh theo hai hướng Đông Bắc (trên 30%) và Tây Nam (trên 15%) với tốc độ gió phổ biến từ 3,6 m/s đến 7,7 m/s.
  • Tính toán chính xác các ngưỡng nồng độ pha loãng tương đối tới hạn gồm $3,09 \times 10^{-4}\text{ s/m}^3$ (liều 2 giờ) và $5,55 \times 10^{-6}\text{ s/m}^3$ (liều 30 ngày) nhằm đảm bảo liều chiếu tuyến giáp không vượt quá 3 Sv theo đúng quy chuẩn an toàn hạt nhân.
  • Xác lập cơ sở khoa học để khoanh vùng bán kính vùng cấm dân cư từ 800 m đến 1.239 m và vùng hạn chế dân cư kéo dài từ 12 km đến 16,5 km theo trục gió mùa thịnh hành, bảo vệ an toàn cho 310.561 cư dân địa phương.
  • Đóng góp luận cứ kỹ thuật quan trọng cho công tác quy hoạch không gian năng lượng, sẵn sàng tích hợp các công cụ dự báo khí động học hiện đại phục vụ chiến lược phát triển năng lượng bền vững và an toàn bức xạ quốc gia trong tương lai.