Tổng quan nghiên cứu

Dự án Nhà máy Điện hạt nhân Ninh Thuận 2 được quy hoạch với quy mô công suất 4x1000 MW tại xã Vĩnh Hải, huyện Ninh Hải, tỉnh Ninh Thuận, đóng vai trò chiến lược trong việc đảm bảo an ninh năng lượng quốc gia và giảm phát thải khí nhà kính. Tuy nhiên, hoạt động vận hành của nhà máy đòi hỏi hệ thống làm mát công suất lớn với lưu lượng xả thải lên tới 330.000 m3/h. Dòng nước làm mát sau quá trình trao đổi nhiệt có mức nhiệt độ chênh lệch cao hơn môi trường nước biển tự nhiên khoảng 7°C, tạo ra nguồn phát thải nhiệt liên tục ra vùng biển ven bờ.

Vấn đề nghiên cứu trọng tâm là xác định định lượng quy luật lan truyền và khuếch tán nhiệt từ dòng nước thải làm mát, từ đó đánh giá mức độ rủi ro đối với môi trường nước và hệ sinh thái biển nhạy cảm. Vùng biển Vĩnh Hải (tọa độ 109°10’E và 11°38’N) nằm trong vùng đệm của Vườn Quốc gia Núi Chúa và Khu bảo tồn biển, nơi sở hữu quần thể 333 loài san hô tạo rạn cùng bãi đẻ tự nhiên của rùa biển quý hiếm.

Mục tiêu cụ thể của luận văn là ứng dụng mô hình số trị thủy động lực học để mô phỏng trường nhiệt độ dư theo các kịch bản khí tượng thủy văn đặc trưng trong năm, xác định ranh giới ảnh hưởng nhiệt tại các tầng nước, và đề xuất các giải pháp kỹ thuật nhằm giảm thiểu tác động sinh thái. Phạm vi nghiên cứu bao trùm toàn bộ diện tích 180 ha mặt bằng nhà máy và vùng biển phụ cận, dựa trên chuỗi số liệu quan trắc thực địa liên tục từ năm 2011 đến năm 2012. Kết quả nghiên cứu mang ý nghĩa thực tiễn quan trọng, cung cấp cơ sở dữ liệu định lượng giúp các nhà quản lý tối ưu hóa thiết kế cửa xả ngầm ở độ sâu -14m, kiểm soát mức tăng nhiệt độ không vượt quá ngưỡng sinh thái 29°C của các loài sinh vật biể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 được xây dựng trên nền tảng cơ học chất lỏng tính toán và động lực học biển ven bờ, trọng tâm là hệ phương trình Navier-Stokes hai chiều tích phân theo độ sâu cho chất lỏng không nén được dưới giả thiết thủy tĩnh và giả thiết Boussinesq. Quá trình vận chuyển vật chất và năng lượng trong môi trường biển được mô tả thông qua hệ phương trình bảo toàn khối lượng, phương trình động lượng theo hai phương nằm ngang, kết hợp chặt chẽ với phương trình truyền tải - bình lưu - khuếch tán nhiệt và độ muối.

Khung lý thuyết tích hợp các khái niệm then chốt gồm: trường nhiệt dư (chênh lệch nhiệt độ giữa dòng thải và môi trường nền), hệ số xoáy rối theo phương ngang xác định theo mô hình Smagorinsky, ứng suất ma sát đáy biển tính toán qua hệ số nhám Manning, và cơ chế trao đổi nhiệt bức xạ qua bề mặt thoáng với khí quyển. Mô hình tiếp cận dựa trên lý thuyết tương tác nhiệt - sinh thái, trong đó tốc độ suy giảm nhiệt lượng phụ thuộc trực tiếp vào các điều kiện động lực học dòng triều, tốc độ gió và độ sâu cột nước.

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

Phương pháp mô hình toán là công cụ chủ đạo thông qua việc ứng dụng bộ phần mềm MIKE 21 Flow Model FM (MIKE 21FM) do Viện Thủy lực Đan Mạch phát triển. Nghiên cứu sử dụng kỹ thuật lưới tam giác phi cấu trúc linh hoạt (Flexible Mesh), cho phép chia nhỏ mắt lưới chi tiết tại khu vực công trình cửa lấy nước, đê chắn sóng dài khoảng 300m và cửa xả ngầm, đồng thời mở rộng lưới ra ngoài khơi để tiết kiệm tài nguyên tính toán mà vẫn đảm bảo độ chính xác cao bám sát đường bờ.

Cỡ mẫu dữ liệu quan trắc phục vụ mô hình bao gồm chuỗi đo đạc nhiệt độ nước biển liên tục với tần suất 10 phút một lần suốt 377 ngày (từ ngày 26/6/2011 đến ngày 6/7/2012) tại 2 trạm chuyên dụng TS1 và TS2, đo đồng thời ở 5 tầng độ sâu từ -1m đến -15m với hơn 50.000 điểm dữ liệu thực đo. Phương pháp chọn mẫu là quan trắc cố định tự động kết hợp khảo sát mặt rộng định kỳ về thủy triều, sóng biển (đo được chiều cao sóng cực đại 7,50m) và đa dạng sinh học. Lý do lựa chọn mô hình MIKE 21FM là khả năng xử lý bài toán vùng nước nông có địa hình phức tạp, giải quyết tốt bài toán khô ướt bãi triều và mô phỏng chính xác sự xáo trộn nhiệt độ dưới tác động của hoàn lưu gió mùa. Mô hình được hiệu chỉnh và kiểm định độc lập cho hai mùa gió đặc trưng (mùa hè từ ngày 4 đến 26/8/2011 và mùa đông từ ngày 3 đến 26/12/2011), chất lượng mô phỏng được đánh giá chuẩn hóa qua chỉ số hiệu quả Nash-Sutcliffe đạt giá trị trên 0,85.

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

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

Kết quả mô phỏng số trị đã làm sáng tỏ cơ chế lan truyền nhiệt của dòng thải nước làm mát qua 4 phát hiện trọng tâm:

Thứ nhất, trường nhiệt độ dư có sự phân hóa không gian rõ rệt theo hai mùa gió chính. Vào mùa đông, dưới tác động của gió mùa Đông Bắc chiếm tần suất 39,9% đến 50,3% và tốc độ gió trung bình 2,5 m/s đến 3,8 m/s, vệt nước ấm bị đẩy dọc theo đường bờ trôi về phía Tây Nam. Ngược lại, vào mùa hè, gió Tây Nam chiếm tần suất 22,1% đến 31,4% kết hợp dòng nhật triều đẩy dòng nhiệt lan tỏa rộng ra vùng biển khơi phía Đông Bắc.

Thứ hai, cấu trúc cửa xả ngầm dạng 4 vòi phun hướng lên ở độ sâu -14m giúp thúc đẩy quá trình xáo trộn mật độ ban đầu rất mạnh mẽ. Dòng thải ban đầu có độ chênh nhiệt độ +7°C (đạt 32°C đến 38°C tại vòi phun) nhanh chóng giảm nhiệt lượng khi hòa trộn với khối nước biển tầng sâu. Độ chênh lệch nhiệt độ tại tầng mặt giảm xuống dưới 1,5°C ngoài bán kính 500m tính từ tâm cửa xả.

Thứ ba, mức độ ảnh hưởng đến 32 ha rạn san hô sống tại bãi rạn Thái An (Bãi Thịt rộng 210 ha) có sự biến thiên lớn theo mùa. Trong mùa đông, nhiệt độ nước biển tại bãi rạn được duy trì an toàn dưới 26°C. Tuy nhiên, vào mùa hè, nhiệt độ nền tự nhiên vốn đã đạt mức cao từ 27,2°C đến 31,4°C, khi cộng thêm nhiệt độ dư từ 0,5°C đến 1,0°C của dòng thải làm cho nhiệt độ vùng nước rạn san hô tiếp cận hoặc vượt ngưỡng chịu đựng cực hạn 29°C trong các thời điểm triều kiệt.

Thứ tư, khu vực bãi đẻ của rùa biển tại bãi Móng Tay và bãi Ngang ở phía Bắc nhà máy nằm ngoài vùng lõi chịu ảnh hưởng nhiệt trực tiếp ở tầng đáy. Tuy nhiên, mức tăng nhiệt bề mặt ven bờ khoảng 0,3°C đến 0,6°C trong mùa hè có thể làm thay đổi nhẹ nhiệt độ cát bãi đẻ, tác động đến tỷ lệ đực - cái của rùa non khi tỷ lệ sống sót tự nhiên chỉ đạt khoảng 1/10.000 cá thể.

Thảo luận kết quả

Nguyên nhân cốt lõi dẫn đến sự suy giảm nhiệt nhanh chóng ngoài khơi là nhờ chế độ động lực biển mạnh mẽ tại vùng biển Ninh Thuận với biên độ nhật triều từ 1,5m đến 2,0m kết hợp hệ thống sóng biển có chiều cao trung bình trên 1,0m. So sánh với các nghiên cứu trung tâm nhiệt điện ven biển miền Trung sử dụng cửa xả kênh hở trên mặt, phương án xả ngầm sâu -14m của Nhà máy Điện hạt nhân Ninh Thuận 2 cho thấy hiệu quả hòa trộn nhiệt cao hơn khoảng 45%, giúp diện tích vùng nước bị tăng nhiệt trên 2°C giảm đi đáng kể.

Trong phân tích học thuật, các dữ liệu lan truyền nhiệt được biểu diễn tối ưu qua hệ thống biểu đồ mặt cắt nhiệt độ thẳng đứng tương quan giữa tầng mặt và độ sâu -15m theo thời gian, kết hợp với các bản đồ đẳng nhiệt hai chiều phủ lớp màu (contour maps) tương ứng với từng kịch bản vận hành mùa hè (S1, S2, FS1, FS2) và mùa đông (W1, W2, FW1, FW2). Bảng thống kê diện tích vùng ảnh hưởng nhiệt độ cho thấy giải pháp kỹ thuật xả sâu đã khoanh vùng được phạm vi quá nhiệt 29°C trong phạm vi bán kính an toàn, hạn chế tối đa nguy cơ tẩy trắng trên diện rộng đối với 333 loài san hô thuộc Vườn Quốc gia Núi Chúa.

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

Dựa trên các kết quả tính toán định lượng, nghiên cứu đề xuất 4 nhóm giải pháp kỹ thuật và quản lý môi trường mang tính khả thi cao:

Một là, tối ưu hóa góc nghiêng và vận tốc phun của hệ thống cửa xả ngầm. Ban Quản lý dự án phối hợp với đơn vị tư vấn thiết kế điều chỉnh 4 vòi phun đường kính 2,4m nghiêng một góc 45 độ so với phương thẳng đứng tại độ sâu -14m, nhằm tối đa hóa quãng đường dòng tia chuyển động trong khối nước biển, nâng hiệu suất xáo trộn ban đầu thêm 40%, hoàn thành trong giai đoạn lập thiết kế kỹ thuật thi công (thời gian thực hiện từ 12 đến 18 tháng).

Hai là, thiết lập mạng lưới trạm quan trắc tự động nhiệt độ và chất lượng nước theo thời gian thực. Sở Tài nguyên và Môi trường tỉnh Ninh Thuận chủ trì lắp đặt 5 trạm phao quan trắc online tại khu vực cửa lấy nước, tâm cửa xả, bãi rạn Bãi Thịt và bãi đẻ rùa biển, mục tiêu duy trì sai số kiểm soát nhiệt độ dưới 0,2°C và phát cảnh báo sớm khi nhiệt độ vượt 28,5°C (hoàn thành trước thời điểm chạy thử nghiệm tổ máy số 1).

Ba là, xây dựng quy trình điều tiết công suất làm mát linh hoạt theo mùa và pha triều. Đơn vị vận hành nhà máy áp dụng quy chế chủ động giảm tải lưu lượng làm mát từ 10% đến 15% vào các đợt triều ròng kết hợp nhiệt độ nền vượt 30°C trong các tháng 5, 6 và 7, nhằm bảo vệ an toàn cho 32 ha san hô sống (thực hiện định kỳ hằng năm trong suốt vòng đời vận hành nhà máy).

Bốn là, thực thi chương trình bảo tồn và phục hồi hệ sinh thái biển chuyên biệt. Ban Quản lý Vườn Quốc gia Núi Chúa triển khai các biện pháp nhân tạo hỗ trợ ấp trứng rùa biển tại bãi Móng Tay nhằm duy trì tỷ lệ nở thành công trên 85%, đồng thời tiến hành cấy ghép phục hồi 10 ha rạn san hô tại các khu vực chịu áp lực nhiệt trong khung thời gian 5 năm.

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

Luận văn là tài liệu tham khảo chuyên sâu và hữu ích cho 4 nhóm đối tượng chính:

Thứ nhất, các kỹ sư tư vấn thiết kế công trình ven biển và hạ tầng năng lượng. Luận văn cung cấp phương pháp luận chi tiết và các thông số thực nghiệm về bố trí cửa nhận nước sâu từ 0 đến -25m và đường ống xả ngầm dài 1.300m, hỗ trợ thiết kế các hệ thống làm mát tối ưu cho nhà máy nhiệt điện, điện hạt nhân hoặc điện khí LNG.

Thứ hai, các nhà khoa học, giảng viên và học viên cao học ngành Kỹ thuật Tài nguyên nước, Hải dương học và Khoa học Môi trường. Tài liệu mang lại case-study chuẩn mực về quy trình ứng dụng mô hình MIKE 21FM, từ xây dựng lưới phi cấu trúc, thiết lập điều kiện biên đến kỹ thuật hiệu chỉnh và kiểm định mô hình với chuỗi số liệu thực tế 377 ngày.

Thứ ba, cơ quan quản lý nhà nước về tài nguyên, môi trường và năng lượng. Báo cáo cung cấp cơ sở khoa học định lượng giúp các hội đồng thẩm định đánh giá chính xác Báo cáo đánh giá tác động môi trường (ĐTM) và cấp phép xả thải đối với các dự án ven biển quy mô lớn.

Thứ tư, ban quản lý các khu bảo tồn biển và vườn quốc gia. Cung cấp dữ liệu sinh thái chi tiết về ngưỡng chịu nhiệt 29°C của 333 loài san hô và tập tính làm tổ của rùa biển, phục vụ xây dựng kế hoạch quản lý môi trường thích ứng dài hạn.

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

Nước làm mát từ nhà máy điện hạt nhân tác động tiêu cực đến môi trường biển như thế nào? Dòng nước thải với lưu lượng lớn 330.000 m3/h có nhiệt độ cao hơn nước biển tự nhiên khoảng 7°C làm suy giảm lượng oxy hòa tan, tăng độ nhớt động học và gây sốc nhiệt cục bộ cho sinh vật phù du, đồng thời đe dọa các loài nhạy cảm nhiệt như san hô và cỏ biển.

Tại sao mô hình MIKE 21FM lại phù hợp cho bài toán khuếch tán nhiệt ven bờ? MIKE 21FM sử dụng phương pháp lưới phần tử hữu hạn tam giác phi cấu trúc, giúp mô phỏng chân thực đường bờ biển Vĩnh Hải khúc khuỷu, các công trình cảng biển và chế độ dòng chảy nhật triều với độ chính xác cao hơn nhiều so với lưới ô vuông truyền thống.

Rạn san hô tại Vườn Quốc gia Núi Chúa chịu tổn thương ra sao khi nhiệt độ tăng? San hô cứng có ngưỡng chịu nhiệt giới hạn khoảng 29°C. Khi nhiệt độ vượt ngưỡng này do nhiệt thải kết hợp nắng nóng mùa hè, tảo cộng sinh Zooxanthellae bị trục xuất khiến san hô bị tẩy trắng và chết hàng loạt trên diện tích rạn 32 ha.

Dòng thải nhiệt có ảnh hưởng đến hoạt động sinh sản của rùa biển không? Rùa biển làm tổ tại các bãi cát ven biển từ bãi Móng Tay đến bãi Ngang với 90 đến 130 trứng mỗi ổ. Nhiệt độ nước và cát tăng cao sẽ rút ngắn thời gian ấp, làm lệch tỷ lệ giới tính sang rùa cái và cản trở đường bơi của rùa con ra đại dương.

Giải pháp xả ngầm ở độ sâu -14m mang lại ưu thế gì so với xả tràn bề mặt? Xả ngầm ở độ sâu -14m qua 4 cửa xả hướng lên trên giúp tận dụng thế năng và động năng của dòng tia để hòa trộn nhanh chóng nước thải nóng với khối nước biển lạnh tầng đáy, giảm độ chênh nhiệt bề mặt xuống dưới 1,5°C trong phạm vi 500m.

Kết luận

Nghiên cứu đã giải quyết trọn vẹn bài toán đánh giá tác động nhiệt từ dòng thải làm mát của Nhà máy Điện hạt nhân Ninh Thuận 2 thông qua 5 kết luận cốt lõi:

  • Thiết lập thành công mô hình toán thủy động lực học MIKE 21FM với lưới phi cấu trúc, phản ánh chính xác quy luật lan truyền nhiệt của dòng thải 330.000 m3/h trong môi trường biển Vĩnh Hải.
  • Xác định rõ sự khác biệt mùa vụ trong cơ chế khuếch tán nhiệt: mùa đông gió Đông Bắc đẩy dòng nhiệt trôi dạt dọc bờ về phía Tây Nam, mùa hè gió Tây Nam mở rộng trường nhiệt ra phía ngoài khơi.
  • Lượng hóa cụ thể rủi ro sinh thái đối với 333 loài san hô cứng tại bãi rạn Thái An, chỉ rõ nguy cơ tẩy trắng tập trung chủ yếu vào các tháng mùa hè khi nhiệt độ nước nền vượt 30°C.
  • Chứng minh hiệu quả vượt trội của giải pháp công trình xả ngầm sâu -14m qua đường ống 1.300m trong việc thúc đẩy quá trình hòa trộn nhiệt tầng sâu và thu hẹp diện tích ô nhiễm nhiệt bề mặt.
  • Đề xuất hệ thống giải pháp đồng bộ kết hợp tối ưu hóa cửa xả, quan trắc tự động 5 trạm thời gian thực và quy trình điều tiết lưu lượng vận hành linh hoạt.

Đóng góp chính của luận văn là xây dựng bộ dữ liệu khoa học định lượng có độ tin cậy cao, phục vụ trực tiếp cho quá trình hoàn thiện thiết kế kỹ thuật và thẩm định môi trường dự án điện hạt nhân tại Việt Nam. Bước đi tiếp theo cần mở rộng mô hình hóa không gian 3 chiều (MIKE 3) và tích hợp các kịch bản biến đổi khí hậu dài hạn. Các cơ quan quản lý và đơn vị phát triển dự án cần nhanh chóng ứng dụng các khuyến nghị này vào quy hoạch chi tiết để đảm bảo mục tiêu phát triển năng lượng bền vững song hành với bảo tồn hệ sinh thái biển.