Tổng quan nghiên cứu

Dải bờ biển Việt Nam kéo dài 3.260 km, giữ vị thế chiến lược trong phát triển kinh tế biển nhưng đồng thời cũng là nơi chịu tác động nặng nề của thiên tai khí tượng thủy văn. Thống kê hàng năm cho thấy có khoảng 5 đến 6 cơn bão đổ bộ trực tiếp vào đất liền nước ta, gây ra hiện tượng nước dâng nghiêm trọng. Điển hình trong năm 2005, hai cơn bão Washi và Damrey đổ bộ trùng thời kỳ triều cường đã tàn phá hàng loạt đê biển tại Hải Phòng và Nam Định. Trên thế giới, thảm họa nước dâng do bão từng ghi nhận độ cao lên tới 7,0 m tại Bangladesh năm 1970 cướp đi hơn 300.000 sinh mạng, hay bão Katrina năm 2005 tại Mỹ gây nước dâng 7,6 m làm ngập 80% diện tích thành phố New Orleans.

Vấn đề nghiên cứu then chốt xuất phát từ tính chất phức tạp của địa hình đáy biển, cấu trúc đường bờ quanh co và đặc biệt là độ bất định lớn trong dự báo quỹ đạo bão hạn 24 giờ. Mục tiêu cụ thể của luận văn là nghiên cứu tính toán đường bao cực đại của nước dâng và xác định xác suất xuất hiện mực nước dâng do bão tại dải ven bờ Biển Đông Việt Nam. Phạm vi nghiên cứu bao phủ toàn bộ vùng Biển Đông từ kinh độ 98°E đến 118°E, vĩ độ 5°N đến 24°N, tập trung đánh giá chi tiết ba khu vực Bắc Bộ, Trung Bộ và Nam Bộ qua ba cơn bão tiêu biểu: bão Côn Sơn (năm 2010), bão Xangsane (năm 2006) và bão Durian (năm 2006).

Ý nghĩa thực tiễn của công trình thể hiện qua việc cung cấp bộ công cụ tính toán thủy động lực có độ tin cậy cao, hỗ trợ giảm thiểu sai số dự báo ngập lụt ven bờ xuống dưới 0,39 m. Kết quả nghiên cứu đóng vai trò là cơ sở khoa học trực tiếp phục vụ công tác quy hoạch, củng cố đê biển, tối ưu hóa các phương án sơ tán dân cư và bảo vệ an toàn cho hơn 20 triệu người dân sinh sống tại các đô thị, làng chài duyên hải.

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 dựa trên hai nền tảng lý thuyết thủy động lực học và khí tượng biển hiện đại:

  1. Lý thuyết dòng chảy nước nông hai chiều (2D Shallow Water Equations): Mô hình hóa sự chuyển động của khối nước dưới tác động của lực Coriolis, áp suất khí quyển, ma sát đáy và ứng suất gió bề mặt. Phương trình liên tục và phương trình động lượng tích phân theo chiều sâu mô tả chính xác sự lan truyền của sóng dài nước dâng.
  2. Lý thuyết trường gió và khí áp bão Young – Sobey: Xác định phân bố tốc độ gió quay và áp suất tiếp cận tâm bão dựa trên bán kính gió cực đại và áp suất trung tâm bão.

Ba khái niệm học thuật cốt lõi định hình toàn bộ nghiên cứu bao gồm:

  • Nước dâng do bão (Storm surge): Sự dâng cao bất thường của mực nước biển so với mực nước triều thiên văn do ứng suất gió bão và sự sụt giảm khí áp gây ra (theo quy luật khí áp giảm 1 hPa làm mực nước biển dâng khoảng 1 cm).
  • Mực nước tổng cộng trong bão (Storm tide): Giá trị mực nước thực tế tại bờ biển được tổng hợp từ ba thành phần chính gồm thủy triều, nước dâng do bão và nước dâng do sóng vỡ.
  • Đường bao nước dâng cực đại (Maximum Envelope of Water - MEOW): Đường phân bố không gian thể hiện giá trị mực nước dâng cao nhất ghi nhận được tại từng vị trí ven bờ trong toàn bộ thời gian bão di chuyển và đổ bộ.

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

Nghiên cứu sử dụng nguồn dữ liệu địa hình đáy biển toàn cầu GEBCO với độ phân giải 1/32 độ, kết hợp dữ liệu khí tượng bão tái phân tích từ Viện Tin học Quốc gia Nhật Bản (NII) và chuỗi số liệu mực nước thực đo tại 5 trạm hải văn ven bờ: Hòn Dáu, Hòn Ngư, Sơn Trà, Quy Nhơn, Vũng Tàu.

Cỡ mẫu và phương pháp chọn mẫu được thiết lập khoa học: tác giả chọn 3 cơn bão đại diện cho ba phân vùng khí hậu hải văn đặc trưng của Việt Nam (Côn Sơn ở miền Bắc, Xangsane ở miền Trung, Durian ở miền Nam) và 2 cơn bão kiểm định (Damrey năm 2005, Haima năm 2011). Để giải quyết độ bất định trong dự báo khí tượng, nghiên cứu thiết lập tập hợp 13 kịch bản bão giả định cho mỗi cơn bão, gồm 1 quỹ đạo tái phân tích và 12 quỹ đạo tịnh tiến song song có điểm đổ bộ cách nhau 9,6 km, phân bố đều trong phạm vi sai số dự báo vị trí tâm bão 24 giờ là 115,2 km. Lý do lựa chọn phương pháp này nhằm bao quát toàn diện rủi ro nước dâng khi bão lệch hướng so với dự báo chính thức.

Phương pháp phân tích số trị sử dụng mô hình MIKE 21 HD của Viện Thủy lực Đan Mạch (DHI) trên hệ lưới phần tử hữu hạn phi cấu trúc gồm 16.797 nút lưới. Kết quả mô phỏng mực nước sau đó được nội suy không gian bằng phương pháp Kriging (phương sai nhỏ nhất) và đánh giá độ chính xác thông qua chỉ số hiệu quả Nash-Sutcliffe. Timeline nghiên cứu và kiểm định bao quát chuỗi dữ liệu thực đo từ năm 2005 đến năm 2014.

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

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

Quá trình hiệu chỉnh và mô phỏng số trị đã mang lại 4 phát hiện khoa học quan trọng:

Thứ nhất, mô hình MIKE 21 HD đạt độ chính xác xuất sắc trong việc tái hiện dao động mực nước dọc bờ biển Việt Nam. Chỉ số Nash-Sutcliffe tại 5 trạm hải văn hiệu chỉnh đều đạt mức rất cao, dao động từ 0,93 đến 0,98. Khi kiểm định với bão Damrey (2005) tại trạm Hòn Dáu, mực nước cực đại tính toán đạt 2,71 m so với giá trị thực đo 2,32 m (hệ số tương quan 0,91, sai số 0,39 m). Với bão Haima (2011), mô hình đạt hệ số tương quan 0,93 (thực đo 1,68 m, tính toán 1,43 m, sai số 0,25 m).

Thứ hai, đối với bão Côn Sơn (2010) tại Bắc Bộ, vùng chịu ảnh hưởng nước dâng kéo dài từ Quảng Ninh đến Hà Tĩnh. Tại điểm ven bờ Quảng Ninh, xác suất xuất hiện nước dâng từ 2,0 m đến 2,5 m đạt 69,77%, xác suất nước dâng từ 1,5 m đến 2,0 m chiếm 30,23%. Tại điểm ven bờ Hải Phòng, xác suất xuất hiện nước dâng từ 1,5 m đến 2,0 m đạt tới 98,61%, và tại điểm tiếp giáp đạt tuyệt đối 100%.

Thứ ba, đối với bão Xangsane (2006) tại Trung Bộ, dải bờ biển từ Quảng Trị đến Đà Nẵng chịu tác động lớn. Kết quả phân tích xác suất từ 13 kịch bản giả định cho thấy 100% khả năng xuất hiện mực nước dâng trong khoảng từ 0,5 m đến 1,0 m tại các trạm kiểm tra trọng điểm thuộc Thừa Thiên Huế và Đà Nẵng.

Thứ tư, đối với bão Durian (2006) tại Nam Bộ, đường bao nước dâng mở rộng từ Bình Thuận đến Cà Mau. Mực nước dâng cực đại ghi nhận có sự biến thiên lớn giữa các vùng ven biển mở và các cửa sông do tác động giao thoa giữa cấu trúc địa hình đáy nông và chế độ bán nhật triều không đều.

Thảo luận kết quả

Nguyên nhân vật lý của hiện tượng nước dâng cực đại tập trung lệch về phía bên phải của hướng bão di chuyển là do sự kết hợp của hoàn lưu gió xoáy nghịch chiều kim đồng hồ ở Bắc bán cầu với vận tốc tịnh tiến của bão. Khi bão tiếp cận vùng thềm lục địa nông dạng lòng chảo như Vịnh Bắc Bộ (độ sâu trung bình chỉ 40 - 50 m), hiệu ứng dồn nước do gió và sóng vỡ tăng vọt, khiến nước dâng tại bờ biển Quảng Ninh - Hải Phòng cao hơn đáng kể so với khu vực Trung Bộ nơi có sườn lục địa dốc đứng.

So với các công thức kinh nghiệm cổ điển như phương pháp Saffir-Simpson hay công thức thực nghiệm của Cơ quan Khí tượng Nhật Bản (JMA) vốn chỉ tính toán tĩnh và bỏ qua biến động địa hình đáy, cách tiếp cận mô hình số trị động lực kết hợp xác suất 13 kịch bản quỹ đạo giả định đã khắc phục triệt để các sai số cục bộ. Phương pháp này cho phép nhận diện rõ ràng độ bất định 115,2 km trong dự báo bão hạn 24 giờ.

Về mặt biểu diễn dữ liệu, các kết quả mô phỏng được thể hiện trực quan qua bản đồ phân bố đường bao nước dâng cực đại (MEOW) và bảng ma trận xác suất phân tầng độ cao (theo các dải 0,0 - 0,5 m; 0,5 - 1,0 m; 1,0 - 1,5 m; 1,5 - 2,0 m; 2,0 - 2,5 m). Cách trình bày này giúp các nhà quản lý rủi ro thiên tai dễ dàng đọc và chuyển đổi số liệu thành các kịch bản ngập lụt thực địa.

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

Dựa trên các kết quả tính toán đường bao cực đại và xác suất xuất hiện nước dâng, luận văn đề xuất 4 nhóm giải pháp thực thi cấp thiết:

  1. Nâng cấp và kiên cố hóa hệ thống công trình đê biển: Bộ Nông nghiệp và Phát triển nông thôn phối hợp với Ủy ban nhân dân các tỉnh ven biển từ Quảng Ninh đến Nam Định rà soát, nâng cao trình đỉnh đê thêm từ 0,5 m đến 1,0 m tại các đoạn xung yếu. Mục tiêu đảm bảo an toàn đê điều trước các cơn bão mạnh cấp 11-12 kết hợp triều cường, triển khai hoàn thành trong giai đoạn 2026 - 2030.
  2. Tích hợp bộ công cụ mô phỏng xác suất vào quy trình dự báo bão nghiệp vụ: Trung tâm Dự báo Khí tượng Thủy văn Quốc gia đưa mô hình MIKE 21 HD kết hợp 13 kịch bản bão giả định vào vận hành thường trực, phát hành bản đồ cảnh báo nước dâng dạng xác suất trước 24 giờ khi bão đổ bộ, phấn đấu kiểm soát sai số dự báo mực nước dưới 0,3 m từ năm 2027.
  3. Hiện đại hóa và dày đặc hóa mạng lưới quan trắc hải văn: Tổng cục Khí tượng Thủy văn đầu tư lắp đặt thêm các trạm đo mực nước tự động truyền số liệu thời gian thực (tần suất 10 - 15 phút/lần), rút ngắn khoảng cách trung bình giữa các trạm đo dọc bờ biển từ 250 km hiện nay xuống dưới 80 km vào năm 2028.
  4. Cập nhật bản đồ rủi ro ngập lụt và phương án sơ tán dân cư ven biển: Ban Chỉ huy Phòng chống thiên tai và Tìm kiếm cứu nạn tại 28 tỉnh, thành phố ven biển cập nhật kịch bản ứng phó khẩn cấp, thiết lập kế hoạch di dời dân cư chi tiết tại các vùng có xác suất nước dâng trên 1,5 m vượt quá 50%, thực hiện diễn tập định kỳ hàng năm.

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

Nội dung và kết quả của luận văn mang lại giá trị thực tiễn và học thuật cho 4 nhóm đối tượng chính:

  1. Cơ quan dự báo khí tượng thủy văn và cảnh báo thiên tai: Ứng dụng quy trình thiết lập 13 kịch bản bão giả định trên nền sai số dự báo 24 giờ để nâng cao chất lượng bản tin cảnh báo nước dâng và triều cường.
  2. Đơn vị tư vấn thiết kế công trình thủy và quy hoạch ven biển: Khai thác số liệu đường bao nước dâng cực đại và chu kỳ mực nước cực trị để tính toán khẩu độ cống hàn, cao trình đê biển, cảng biển và khu neo đậu tàu thuyền.
  3. Cơ quan quản lý nhà nước và hoạch định chính sách phòng chống thiên tai: Sử dụng bản đồ phân vùng xác suất ngập lụt làm căn cứ phê duyệt quy hoạch sử dụng đất duyên hải, xây dựng kế hoạch ứng phó biến đổi khí hậu và bảo hiểm rủi ro thiên tai.
  4. Giảng viên, nghiên cứu sinh và sinh viên chuyên ngành Hải dương học, Khí tượng học, Biến đổi khí hậu: Tiếp cận phương pháp kết hợp mô hình trường gió Young – Sobey, mô hình thủy lực 2 chiều MIKE 21 và kỹ thuật nội suy Kriging làm tài liệu tham khảo học thuật chuyên sâu.

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

  1. Nước dâng do bão khác với mực nước tổng cộng trong bão như thế nào? Nước dâng do bão là độ chênh lệch mực nước biển dâng cao bất thường do gió bão và giảm áp suất khí quyển. Mực nước tổng cộng trong bão là tổng đại số của nước dâng do bão, dao động thủy triều và nước dâng do sóng. Trong thực tế, bão đổ bộ vào thời kỳ triều cường sẽ đẩy mực nước tổng cộng lên rất cao, gây vỡ đê và ngập lụt nghiêm trọng như trường hợp bão Damrey năm 2005.

  2. Tại sao nghiên cứu lại chọn 13 kịch bản quỹ đạo bão giả định cách nhau 9,6 km? Theo báo cáo của ngành khí tượng, sai số dự báo vị trí tâm bão trung bình hạn 24 giờ của Việt Nam là 115,2 km. Việc chia khoảng cách này thành 12 đoạn đều nhau tạo ra 13 kịch bản bão giả định (mỗi kịch bản cách nhau 9,6 km) giúp bao quát toàn bộ vùng bất định, từ đó xác định chính xác xác suất xuất hiện nước dâng cực đại tại từng điểm bờ biển.

  3. Độ tin cậy của mô hình MIKE 21 HD trong tính toán nước dâng tại Việt Nam đạt mức nào? Mô hình MIKE 21 HD có độ tin cậy rất cao khi được hiệu chỉnh chuẩn xác. Kết quả kiểm định tại 5 trạm hải văn lớn cho chỉ số Nash-Sutcliffe từ 0,93 đến 0,98. Khi kiểm định trong bão Damrey và Haima tại trạm Hòn Dáu, hệ số tương quan đạt 0,91 - 0,93 và sai số tuyệt đối lớn nhất chỉ từ 0,25 m đến 0,39 m.

  4. Vì sao vùng ven biển Bắc Bộ có nguy cơ nước dâng do bão cao hơn miền Trung và miền Nam? Vịnh Bắc Bộ có thềm lục địa nông, độ sâu trung bình 40 - 50 m, diện tích vùng nước dưới 30 m chiếm tới 60%, lại có dạng hình lòng chảo kín ba bề. Khi gió bão thổi mạnh hướng vào bờ, khối nước bị dồn nén mãnh liệt tạo nên đỉnh nước dâng từ 2,0 m đến 2,5 m, cao hơn nhiều so với vùng bờ biển dốc đứng ở miền Trung.

  5. Làm thế nào để ứng dụng kết quả đường bao cực đại vào xây dựng đê biển? Đường bao cực đại cung cấp giá trị mực nước dâng cao nhất có thể xảy ra dọc tuyến bờ. Các kỹ sư công trình thủy sẽ cộng giá trị này với mực nước triều cao nhất và độ dâng do sóng để xác định cao trình an toàn đỉnh đê, bổ sung hệ số an toàn từ 0,5 m đến 1,0 m nhằm ngăn chặn hiện tượng sóng tràn qua đê.

Kết luận

  • Luận văn đã xây dựng thành công chuỗi công cụ tính toán tích hợp giữa mô hình trường gió Young – Sobey và mô hình thủy động lực MIKE 21 HD trên lưới phi cấu trúc 16.797 nút lưới bao phủ toàn bộ Biển Đông.
  • Hiệu chỉnh và kiểm định mô hình đạt độ chính xác cao với chỉ số Nash-Sutcliffe từ 0,93 đến 0,98 tại 5 trạm hải văn, sai số mô phỏng mực nước cực đại trong các cơn bão thực tế được khống chế dưới 0,39 m.
  • Xác lập phương pháp tính toán xác suất xuất hiện nước dâng cực đại dựa trên tập hợp 13 kịch bản bão giả định, phản ánh trung thực vùng bất định dự báo 115,2 km hạn trước 24 giờ.
  • Xây dựng chi tiết bản đồ đường bao nước dâng cực đại cho 3 cơn bão điển hình Côn Sơn (Bắc Bộ), Xangsane (Trung Bộ) và Durian (Nam Bộ), chỉ ra xác suất nước dâng trên 2,0 m tại Quảng Ninh lên tới 69,77%.
  • Đóng góp luận cứ khoa học quan trọng cho công tác quy hoạch đê biển, phân vùng rủi ro thiên tai và thiết lập hệ thống cảnh báo sớm ngập lụt ven bờ.

Về lộ trình phát triển tiếp theo trong giai đoạn 2026 - 2028, nghiên cứu cần mở rộng tích hợp tương tác sóng - dòng chảy - thủy triều phi tuyến 3 chiều và ứng dụng dữ liệu viễn thám độ phân giải siêu cao. Hãy liên hệ và khai thác ngay kết quả luận văn này để tối ưu hóa các giải pháp bảo vệ dải ven biển Việt Nam trước hiểm họa thiên tai biển.