Tổng quan nghiên cứu

Vịnh Bắc Bộ là một trong những khu vực chịu ảnh hưởng nặng nề nhất bởi bão và áp thấp nhiệt đới tại Việt Nam, với tần suất trung bình từ 4 đến 5 cơn bão quét qua mỗi năm và vận tốc gió giật có thể vượt mốc 50 m/s. Theo các thống kê lịch sử hải văn ven bờ, cứ 100 cơn bão đổ bộ vào dải bờ biển nước ta thì có tới 10 cơn gây nước dâng vượt 2,0 m và 30 cơn gây nước dâng trên 1,5 m. Hiện tượng nước dâng do bão khi kết hợp với triều cường tại vùng biển phía Bắc thường duy trì đỉnh cao từ 2 đến 3 giờ, gây ngập lụt nghiêm trọng, phá vỡ hệ thống đê điều và cuốn trôi mùa màng. Cơn bão số 7 vào năm 2005 là minh chứng điển hình khi tàn phá nặng nề các tỉnh ven biển, gây thiệt hại kinh tế hơn 747 tỷ đồng tại Thanh Hóa và 517 tỷ đồng tại Nam Định.

Vấn đề nghiên cứu trọng tâm của đề tài là làm sáng tỏ bản chất động lực học và thiết lập mối quan hệ thực nghiệm định lượng giữa dao động mực nước ven bờ với các tham số đặc trưng của bão. Mục tiêu cụ thể là xây dựng các công thức hồi quy thực nghiệm đáng tin cậy, giúp xác định nhanh trị số nước dâng và nước rút cực đại phục vụ công tác cảnh báo sớm.

Phạm vi nghiên cứu bao quát toàn bộ vùng biển ven bờ Vịnh Bắc Bộ trải dài từ vĩ độ 17 độ Bắc đến 21 độ Bắc, tập trung chuyên sâu tại hai trạm quan trắc hải văn chiến lược gồm trạm Hòn Dấu tại Hải Phòng và trạm Hòn Ngư tại Nghệ An trên chuỗi số liệu quan trắc liên tục 32 năm từ năm 1980 đến năm 2012. Ý nghĩa thực tiễn của công trình thể hiện qua việc cung cấp công cụ tính toán nhanh với thời gian xử lý dưới 30 phút, hỗ trợ bảo vệ an toàn cho hơn 763 km đường bờ biển Vịnh Bắc Bộ và hàng triệu cư dân sinh sống tại các vùng trũng thấp 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

Luận văn vận dụng lý thuyết phân tích điều hòa thủy triều dựa trên nguyên lý chồng chập sóng điều hòa của Darwin và Doodson để phân tách dao động mực nước biển thành hai thành phần độc lập. Thành phần thứ nhất là dao động tuần hoàn do lực hấp dẫn thiên văn từ Mặt Trăng và Mặt Trời gây ra. Thành phần thứ hai là dao động phi tuần hoàn hay còn gọi là mực nước dư, đại diện cho hiện tượng nước dâng hoặc nước rút do áp suất khí quyển và trường gió bão tác động lên mặt biển.

Mô hình nghiên cứu xây dựng dựa trên lý thuyết hồi quy thống kê tuyến tính và phi tuyến, kết hợp ma trận tương quan Pearson để xác định mức độ phụ thuộc giữa mực nước dư cực trị với tập hợp các tham số bão. Bốn khái niệm khoa học cốt lõi được định nghĩa xuyên suốt bao gồm: hằng số điều hòa biên độ và pha thiên văn, mực nước dư cực đại, áp suất khí quyển cực tiểu tại tâm bão, và khoảng cách động lực học từ tâm bão đến trạm quan trắc ven bờ.

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

Nguồn dữ liệu quan trắc mực nước từng giờ được thu thập đồng bộ tại trạm Hòn Dấu trong giai đoạn từ năm 1980 đến 2012 với tổng số 289.296 giờ quan trắc và trạm Hòn Ngư trong giai đoạn từ năm 1980 đến 2008. Cơ sở dữ liệu khí tượng bão được khai thác từ nguồn dữ liệu quỹ đạo chuẩn Best Track của cơ quan khí tượng hải quân quốc tế với chuỗi thời gian 50 năm từ năm 1960 đến 2009.

Cỡ mẫu nghiên cứu bao gồm 909 cơn bão và áp thấp nhiệt đới hoạt động trên Biển Đông trong 30 năm từ 1980 đến 2009, trong đó tác giả chọn lọc 63 cơn bão mạnh có vận tốc gió cực đại đạt từ 64 hải lý/giờ (tương đương 32,9 m/s) trở lên. Phương pháp chọn mẫu tập trung vào các cơn bão có quỹ đạo đổ bộ vào phía Nam của trạm quan trắc như khu vực Nam Định, Thanh Hóa để khảo sát tác động dâng nước mạnh nhất ở bán khuyên nguy hiểm bên phải đường đi của bão.

Phương pháp phân tích điều hòa bằng thuật toán bình phương tối thiểu được sử dụng để trích xuất gần 30 sóng triều chính, giúp khử bỏ triệt để thành phần thủy triều khỏi chuỗi số liệu thực đo. Toàn bộ các thuật toán xử lý dữ liệu, đồng hóa tọa độ không gian và phân tích hồi quy đa biến được lập trình tự động hóa bằng ngôn ngữ Fortran. Lý do lựa chọn phương pháp thống kê thực nghiệm xuất phát từ tính ưu việt trong việc phản ánh trung thực đặc thù địa hình vịnh nông có độ sâu trung bình từ 40 m đến 50 m, khắc phục triệt để các sai số do lưới tính toán quá thưa trong các mô hình số trị thủy động truyền thống.

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

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

Thứ nhất, kết quả thống kê tần suất cho thấy sự cân bằng giữa các đợt nước dâng và nước rút tại trạm Hòn Dấu với 12.175 đợt dâng và 12.177 đợt rút trong 32 năm. Tuy nhiên, các đợt nước dâng đáng kể có độ cao trên 50 cm chỉ xuất hiện 142 lần, chiếm tỷ lệ 1,17%. Tại trạm Hòn Ngư, số đợt nước dâng trên 50 cm lên tới 352 lần, chiếm tỷ lệ 4,49% trên tổng số 7.844 đợt dâng, chứng tỏ khu vực bờ biển Nghệ An chịu biên độ dao động nước dâng mạnh hơn đáng kể so với khu vực Hải Phòng.

Thứ hai, nghiên cứu đã xác lập các mốc giá trị cực trị lịch sử của nước dâng và nước rút do bão. Tại trạm Hòn Dấu, mực nước dâng lớn nhất ghi nhận đạt 161,1 cm trong cơn bão ngày 23 đến 25 tháng 8 năm 1996. Tại trạm Hòn Ngư, đỉnh nước dâng kỷ lục đạt 198,2 cm xuất hiện vào ngày 29 tháng 8 năm 1990. Về hiện tượng nước rút, giá trị âm sâu nhất tại Hòn Dấu đạt âm 91,5 cm vào ngày 10 tháng 11 năm 1990 và tại Hòn Ngư đạt âm 143,5 cm vào ngày 18 tháng 5 năm 2007.

Thứ ba, mối quan hệ giữa độ lớn nước dâng với khoảng cách từ tâm bão đến trạm quan trắc tuân theo quy luật phi tuyến rõ rệt. Điển hình trong cơn bão Charlotte năm 1962, quan hệ bậc 6 trên toàn dải khoảng cách đạt hệ số tương quan R bình phương bằng 0,8282. Khi phân tích riêng phạm vi cận bờ dưới 400 km, mối quan hệ chuyển sang dạng hàm đa thức bậc hai với hệ số tương quan R bình phương đạt 0,704. Quy luật thực nghiệm chỉ ra rằng khi bão cách bờ khoảng 400 km, mực nước có xu hướng rút nhẹ, nhưng khi tâm bão tiến sát bờ dưới 100 km, nước dâng tăng vọt theo cấp số nhân.

Thứ tư, phân tích hồi quy đa biến khẳng định mực nước dâng cực đại phụ thuộc chặt chẽ vào sự suy giảm khí áp tâm bão và vận tốc gió cực đại. Khí áp tâm bão khi giảm sâu xuống mức 974 mb so với chuẩn 1013 mb sẽ làm nước biển dâng thêm khoảng 40 cm thuần túy do hiệu ứng áp kế nghịch.

Thảo luận kết quả

Nguyên nhân cốt lõi tạo nên các đợt nước dâng cực trị bắt nguồn từ sự cộng hưởng giữa hiệu ứng hút nước của tâm áp thấp và ứng suất gió bão dồn nước vào lòng chảo Vịnh Bắc Bộ. Địa hình vịnh tương đối nông với hơn 60% diện tích có độ sâu dưới 30 m đã khuếch đại biên độ sóng lùi và sóng dâng khi bão áp sát bờ.

So với các công thức kinh nghiệm quốc tế như công thức của Cheng tại Hồng Kông hay Isozaki tại Nhật Bản, các hệ số hồi quy trong luận văn mang tính địa phương hóa sâu sắc, phù hợp hoàn hảo với đường bờ khúc khuỷu và thềm lục địa phía Bắc Việt Nam. Kết quả này vượt trội hơn các mô hình số trị trước đây vốn sử dụng lưới tính thô 27 km và giả định trên nền nước tĩnh.

Toàn bộ dữ liệu nghiên cứu được biểu diễn trực quan thông qua các bảng thống kê ma trận tương quan nhiều biến và hệ thống đồ thị chuỗi thời gian so sánh mực nước thực đo đường nét liền với mực nước tính toán đường gạch nối. Sự trùng khớp giữa hai đường đồ thị với hệ số tương quan trên 80% khẳng định độ tin cậy vượt trội của mô hình thực nghiệm.

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

Thứ nhất, nâng cấp mạng lưới trạm hải văn quan trắc tự động dọc dải bờ biển Vịnh Bắc Bộ. Bộ Tài nguyên và Môi trường cần chủ trì phân bổ ngân sách trong giai đoạn 2026 - 2028 để thay thế các thiết bị đo cơ học tại 5 trạm quan trắc trọng điểm bằng cảm biến sóng âm và áp suất điện tử, đặt mục tiêu giảm tỷ lệ khuyết số liệu xuống dưới 0,5% và truyền dữ liệu thời gian thực mỗi 10 phút một lần.

Thứ hai, tích hợp bộ phương trình hồi quy thực nghiệm vào phần mềm nghiệp vụ dự báo khí tượng hải văn. Trung tâm Dự báo Khí tượng Thủy văn Quốc gia cần triển khai ứng dụng các hàm tương quan bậc hai giữa khoảng cách bão, khí áp và mực nước dâng vào hệ thống hỗ trợ ra quyết định từ đầu năm 2026, nhằm nâng cao độ chính xác dự báo đỉnh nước dâng cực đại lên trên 85% trước khi bão đổ bộ từ 12 đến 24 giờ.

Thứ ba, điều chỉnh quy chuẩn thiết kế kỹ thuật công trình đê kè và cảng biển ven bờ. Bộ Nông nghiệp và Phát triển Nông thôn phối hợp cùng Bộ Giao thông Vận tải hoàn thành việc rà soát cao trình đỉnh đê tại các tỉnh Hải Phòng, Nam Định, Nghệ An trước năm 2027, bắt buộc bổ sung biên độ an toàn nước dâng tối thiểu từ 1,5 m đến 2,0 m kết hợp với triều cường đỉnh kỳ vào quy chuẩn xây dựng.

Thứ tư, xây dựng bản đồ phân vùng rủi ro ngập lụt do nước dâng bão chi tiết cấp xã phường. Ủy ban nhân dân các tỉnh ven biển Bắc Bộ và Bắc Trung Bộ cần hoàn thành số hóa bản đồ cảnh báo ngập lụt tỷ lệ 1:10.000 trong giai đoạn 2026 - 2030, giúp các địa phương chủ động phương án sơ tán dân cư, hướng tới mục tiêu giảm thiểu 30% thiệt hại về tài sản khi có bão mạnh đổ bộ.

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

Nhóm chuyên viên dự báo khí tượng thủy văn và hải văn có thể ứng dụng trực tiếp các phương trình hồi quy và bảng tra cứu ma trận tương quan để tính toán nhanh mực nước dâng bão trong các ca trực nghiệp vụ khẩn cấp chỉ với 15 phút xử lý thông tin.

Nhóm kỹ sư thiết kế công trình thủy công, cảng biển và đê điều có thể khai thác các số liệu thống kê mực nước dâng cực trị từ 161,1 cm đến 198,2 cm cùng chuỗi số liệu mực nước rút sâu âm 143,5 cm để xác định chuẩn xác tĩnh không cầu cảng, cao trình đỉnh kè và độ sâu an toàn cho luồng tàu hàng hải.

Nhóm cán bộ chỉ huy phòng chống thiên tai và tìm kiếm cứu nạn các cấp có thể tham khảo quy luật tương quan khoảng cách bão dưới 400 km để thiết lập các kịch bản sơ tán dân cư vùng trũng thấp trước 12 giờ khi bão đổ bộ vào đất liền.

Nhóm giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Hải dương học, Khí tượng học và Quản lý biển có thể sử dụng phương pháp tách triều bằng thuật toán bình phương tối thiểu và bộ mã nguồn Fortran làm tài liệu tham khảo chuẩn mực trong giảng dạy và phát triển các đề tài nghiên cứu mở rộng.

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

Phương pháp loại triều theo dự tính thủy triều được thực hiện như thế nào? Phương pháp này áp dụng thuật toán bình phương tối thiểu để phân tích chuỗi số liệu mực nước thực đo liên tục nhiều năm nhằm xác định bộ hằng số điều hòa biên độ và pha của gần 30 sóng triều chính. Sau đó, giá trị triều tính toán được trừ khỏi mực nước tổng cộng để thu được chuỗi mực nước dư phi tuần hoàn chính xác phục vụ nghiên cứu bão.

Tại sao nước dâng cực đại thường xuất hiện ở phía bên phải đường đi của bão? Tại bán cầu Bắc, hoàn lưu bão xoáy ngược chiều kim đồng hồ kết hợp với vận tốc tịnh tiến của tâm bão làm cho trường gió ở nửa bên phải đường đi đạt cường độ mạnh nhất. Hướng gió thổi vuông góc từ biển vào bờ dồn ép khối nước lên vùng đáy nông, tạo ra đỉnh nước dâng cao nhất ở sườn phải bão.

Khoảng cách từ tâm bão đến trạm bờ ảnh hưởng thế nào đến mực nước? Khi tâm bão cách trạm quan trắc khoảng 400 km, trường gió ngoại vi có xu hướng đẩy nước ra xa bờ tạo hiện tượng nước rút nhẹ. Tuy nhiên, khi bão tiến vào cự ly dưới 100 km, mực nước dâng tăng vọt theo hàm phi tuyến bậc hai với hệ số tương quan đạt 0,704, đạt giá trị đỉnh điểm ngay trước khi tâm bão đổ bộ.

Mức độ biến động nước dâng giữa trạm Hòn Dấu và Hòn Ngư khác nhau ra sao? Trạm Hòn Ngư ghi nhận số đợt nước dâng trên 50 cm cao gấp 2,5 lần trạm Hòn Dấu với 352 đợt so với 142 đợt, đỉnh dâng cao nhất đạt 198,2 cm so với 161,1 cm tại Hòn Dấu. Nguyên nhân là do tần suất bão đổ bộ vào khu vực Thanh Nghệ Tĩnh cao hơn vùng ven biển Đông Bắc và đặc điểm đáy biển cục bộ dốc hơn.

Phương pháp thống kê thực nghiệm có ưu điểm gì so với mô hình số trị? Phương pháp thống kê thực nghiệm phản ánh trực tiếp các đặc thù địa phương ven bờ mà không phụ thuộc vào kích thước ô lưới hay các giả định nền nước tĩnh phức tạp của mô hình số trị. Phương pháp này cho phép tính toán tức thời với độ chính xác cao và hệ số tin cậy vượt 80% mà không đòi hỏi hệ thống máy chủ cấu hình lớn.

Kết luận

  • Xây dựng thành công bộ cơ sở dữ liệu hải văn và khí tượng chuẩn hóa gồm 289.296 giờ quan trắc mực nước và 909 cơn bão trong giai đoạn 1980 - 2012.
  • Xác lập chính xác các mốc cực trị nước dâng lịch sử đạt 161,1 cm tại trạm Hòn Dấu và 198,2 cm tại trạm Hòn Ngư.
  • Thiết lập hệ thống phương trình hồi quy phi tuyến bậc cao với hệ số tương quan R bình phương đạt từ 0,704 đến 0,8282.
  • Làm sáng tỏ quy luật dao động mực nước dâng rút theo cự ly tâm bão dưới 400 km và mức suy giảm khí áp xuống 974 mb.
  • Đóng góp luận cứ khoa học thực tiễn vững chắc phục vụ công tác cảnh báo thiên tai và thiết kế an toàn công trình ven biển Vịnh Bắc Bộ.

Đóng góp lớn nhất của luận văn là đã địa phương hóa thành công các công thức dự báo nước dâng bão bằng phương pháp thống kê thực nghiệm chính xác cao. Kế hoạch phát triển tiếp theo trong giai đoạn 2026 - 2030 là mở rộng thuật toán ra toàn bộ dải ven biển miền Trung và tích hợp trí tuệ nhân tạo để tự động hóa dự báo. Hãy ứng dụng ngay các phương trình thực nghiệm và cơ sở dữ liệu của luận văn vào quy hoạch đê điều, xây dựng cảng biển và nâng cao năng lực ứng phó thiên tai cho các địa phương ven biển Vịnh Bắc Bộ.