Tổng quan nghiên cứu

Vịnh Bắc Bộ có diện tích khoảng 150.000 km², chiều dài xấp xỉ 600 km, độ rộng dao động từ 200 đến 320 km và độ sâu trung bình khoảng 50 đến 60 m. Với hệ thống hơn 3.000 hòn đảo lớn nhỏ chiếm diện tích gần 3.000 km² tại vùng Đông Bắc cùng mạng lưới lạch biển phức tạp, đây là ngư trường trọng điểm, trung tâm giao thương hàng hải và du lịch biển lớn của Việt Nam. Tuy nhiên, dải ven bờ Tây vịnh Bắc Bộ từ Quảng Ninh đến Quảng Trị thường xuyên chịu tác động nặng nề của các hiện tượng thời tiết cực đoan như bão nhiệt đới, áp thấp và các đợt gió mùa hoạt động mạnh.

Vấn đề cốt lõi mà nghiên cứu giải quyết là hiện tượng dao động dâng, rút mực nước biển phi tuần hoàn do trường gió và biến thiên áp suất khí quyển gây ra. Sự kết hợp giữa các pha triều cường và nước dâng do gió mạnh có thể làm mực nước dâng cao đột biến trên 1,00 m đến 1,30 m, gây vỡ đê kè, ngập lụt diện rộng các đô thị ven biển và tê liệt hệ thống cảng nước sâu. Mục tiêu cụ thể của luận văn là lượng hóa ảnh hưởng của 6 hướng gió thịnh hành cùng sự suy giảm áp suất khí quyển đến chế độ dâng rút mực nước phi điều hòa, từ đó thiết lập các phương trình tương quan định lượng và bộ thông số thủy động lực học phục vụ dự báo.

Phạm vi nghiên cứu bao quát toàn bộ bờ Tây vịnh Bắc Bộ với 8 vị trí trọng điểm: Mũi Ngọc, Cửa Ông, Hòn Dáu, Ba Lạt, Lạch Trường, Diễn Châu, Vũng Áng và Cửa Tùng. Kết quả nghiên cứu có ý nghĩa thực tiễn to lớn trong việc cung cấp dữ liệu đầu vào chuẩn xác cho công tác quy hoạch đê điều, xây dựng bến cảng, giúp giảm thiểu ước tính khoảng 25% đến 35% thiệt hại kinh tế do tai biến ngập lụt ven bờ gây ra hàng năm.

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 xây dựng trên nền tảng lý thuyết thủy động lực học chất lỏng phi nén hai chiều và ba chiều, ứng dụng hệ phương trình Navier-Stokes kết hợp chặt chẽ với giả thiết Boussinesq và nguyên lý áp suất thủy tĩnh. Quá trình vận chuyển khối nước được chi phối bởi phương trình liên tục, phương trình bảo toàn động lượng theo các trục tọa độ Descartes, cùng các phương trình truyền tải nhiệt, độ muối và nồng độ chất vô hướng.

Nghiên cứu vận dụng 4 khái niệm và cơ chế vật lý nền tảng:

  • Hiệu ứng ma sát tiếp tuyến giữa gió và bề mặt biển: Ứng suất gió truyền động lượng trực tiếp xuống lớp nước mặt, dồn nước về phía bờ gây hiện tượng nước dâng hoặc đẩy nước ra khơi gây nước rút.
  • Hiệu ứng nghịch áp khí quyển (Inverse Barometer Effect): Quy luật thủy tĩnh xác định khi áp suất khí quyển giảm 1 mbar thì mực nước biển dâng cao tĩnh khoảng 10 mm (1 cm) và ngược lại.
  • Hiệu ứng bơm Ekman: Sự hội tụ khối nước tầng mặt dưới tác động của lực Coriolis và trường gió xoáy cục bộ, đẩy lớp nước ấm tích tụ ven bờ.
  • Lý thuyết dao động mực nước biển tổng hợp: Mực nước thực tế là sự chồng chập tuyến tính giữa dao động tuần hoàn (thủy triều thiên văn gồm 60 đến 68 sóng phân triều chính) và dao động phi tuần hoàn (nước dâng, rút khí tượng).

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

Nghiên cứu kết hợp hai phương pháp chủ đạo là mô hình hóa số trị và phân tích thống kê toán học:

Về phương pháp mô hình hóa, tác giả ứng dụng mô-đun thủy lực MIKE 21 FM dựa trên kỹ thuật thể tích hữu hạn trung tâm với lưới tính phi cấu trúc tam giác. Toàn bộ miền tính vịnh Bắc Bộ được rời rạc hóa thành 4.014 nút lưới và 6.830 phần tử tam giác, khoảng cách giữa các nút biến đổi linh hoạt từ 2,1 km tại dải bờ nông ven bờ đến 43 km tại vùng biển khơi sâu. Điều kiện biên lỏng được trích xuất từ mô hình thủy triều toàn cầu có độ phân giải 0,25 độ, kết hợp công cụ tạo trường gió bão chuyên dụng.

Về phương pháp thống kê, chuỗi số liệu quan trắc mực nước thực đo theo từng giờ trong khoảng thời gian 25 ngày (từ ngày 05/01/2005 đến ngày 29/01/2005) tại hai trạm hải văn chuẩn Hòn Dáu và Hòn Ngư được sử dụng để hiệu chỉnh và kiểm nghiệm mô hình. Tác giả áp dụng phương pháp phân tích điều hòa để tách mực nước thủy triều dự tính ra khỏi mực nước thực đo theo biểu thức hiệu số, từ đó cô lập chuỗi mực nước phi tuần hoàn. Mối liên hệ giữa tốc độ gió (thực nghiệm từ 5 m/s đến 25 m/s) với độ dâng rút được xác định qua phương trình hồi quy tuyến tính với tiêu chuẩn đánh giá nghiêm ngặt: hệ số tương quan đạt từ 0,80 trở lên và tỷ số giữa hệ số tương quan trên độ lệch xác suất đạt giá trị lớn hơn 10.

Lý do lựa chọn phương pháp hỗn hợp này là vì mô hình MIKE 21 FM phản ánh xuất sắc hình học đường bờ khúc khuỷu và địa hình đáy phức tạp, trong khi phương pháp thống kê giúp lượng hóa trực tiếp quy luật dâng rút thành các phương trình giải tích tường minh, thuận tiện cho việc tra cứu nhanh trong dự báo tác nghiệp.

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

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

Mô phỏng số trị và phân tích dữ liệu cho thấy trường gió là nhân tố quyết định độ lớn và xu thế dâng rút của mực nước phi tuần hoàn tại bờ Tây vịnh Bắc Bộ:

  • Gió hướng Đông Nam (SE) gây ra hiện tượng nước dâng mạnh nhất trên toàn dải ven bờ. Độ dâng trung bình toàn khu vực đạt 0,34 m. Mực nước dâng cực đại ghi nhận tại trạm Cửa Ông lên tới 1,30 m khi tốc độ gió đạt 25 m/s, tại trạm Mũi Ngọc đạt 1,19 m và tại Hòn Dáu đạt 1,19 m. Hệ số tương quan hồi quy tại các trạm phía Bắc đạt mức rất cao, dao động từ 0,95 đến 0,99.
  • Gió hướng Đông (E) tạo ra mức dâng nước đáng kể thứ hai. Độ dâng trung bình đạt khoảng 0,35 m. Khi gió đạt vận tốc 25 m/s, trạm Diễn Châu ghi nhận mức dâng cao nhất đạt 1,09 m, tiếp đến là Cửa Ông đạt 1,04 m và Lạch Trường đạt 0,94 m, trong khi các trạm phía Nam như Cửa Tùng chỉ dâng tối đa 0,46 m.
  • Gió hướng Đông Bắc (NE) mang đặc trưng của gió mùa mùa đông gây nước dâng trung bình khoảng 0,18 m, thấp hơn khoảng 47% so với gió hướng Đông Nam. Giá trị dâng cao nhất theo hướng này tập trung tại Diễn Châu với 0,74 m và Vũng Áng với 0,64 m ở cấp gió 25 m/s.
  • Các hướng gió Tây (W), Tây Bắc (NW) và Tây Nam (SW) đóng vai trò chủ đạo gây ra hiện tượng nước rút phi tuần hoàn. Khi gió Tây thổi mạnh 25 m/s, mực nước tại các vùng vịnh nông bị hạ thấp đáng kể, gây khó khăn cho luồng tàu ra vào cảng.

Thảo luận kết quả

Quy luật phân bố không gian chỉ ra rằng biên độ nước dâng có xu thế giảm dần rõ rệt từ Bắc xuống Nam dọc theo bờ Tây vịnh Bắc Bộ. Khu vực ven biển Quảng Ninh - Hải Phòng (Mũi Ngọc, Cửa Ông, Hòn Dáu) có độ dâng lớn nhất do hiệu ứng vịnh đóng, đáy biển nông thoải và nhiều cửa sông xen kẽ làm khối nước bị dồn ép mạnh. Ngược lại, khu vực Quảng Bình - Quảng Trị (tiêu biểu là trạm Cửa Tùng) có địa hình mở thông với biển sâu, độ sâu tăng nhanh nên biên độ dâng nước do gió chỉ đạt 0,19 m đến 0,46 m.

Kết quả này hoàn toàn tương đồng với các nghiên cứu trước đây về thủy triều và nước dâng do bão của các tác giả tại Viện Cơ học và Trung tâm Hải văn, khẳng định tính chuẩn xác của mô hình MIKE 21 FM. Về mặt trình bày trực quan, dữ liệu dâng rút phi tuần hoàn được biểu diễn hiệu quả thông qua hệ thống đồ thị tán xạ phân tán (scatter plot) đi kèm đường hồi quy tuyến tính, kết hợp cùng các bảng ma trận số liệu chiết xuất chi tiết theo 5 cấp vận tốc gió (5 m/s, 10 m/s, 15 m/s, 20 m/s và 25 m/s). Cách tiếp cận trực quan này cho phép người đọc nhận diện ngay điểm uốn gia tăng mực nước (bắt đầu tăng vọt khi gió vượt ngưỡng 5 m/s) và dễ dàng tra cứu trị số thiết kế cho từng đoạn bờ cụ thể.

Đề 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 công trình và phi công trình:

  • Hoàn thiện và nâng cấp mạng lưới quan trắc khí tượng hải văn ven bờ: Thiết lập bổ sung tối thiểu 5 trạm đo mực nước tự động độ chính xác cao tại các điểm xung yếu cửa sông (như Ba Lạt, Lạch Trường, Cửa Sót) trong giai đoạn 2026-2028, do Tổng cục Khí tượng Thủy văn và Cục Biển và Hải đảo Việt Nam chủ trì thực hiện.
  • Ứng dụng mô hình số trị MIKE 21 FM vào hệ thống cảnh báo sớm thời gian thực: Tích hợp module thủy động lực học vào quy trình nghiệp vụ dự báo bão, nhằm giảm sai số dự báo nước dâng xuống dưới mức 15 cm trước khi bão đổ bộ ít nhất 24 giờ, phục vụ chỉ đạo sơ tán dân cư an toàn.
  • Rà soát và cập nhật tiêu chuẩn kỹ thuật thiết kế đê biển và công trình cảng: Bộ Xây dựng cùng Bộ Nông nghiệp và Phát triển nông thôn cần ban hành quy chuẩn mới trước năm 2027, trong đó bắt buộc cộng thêm biên độ nước dâng phi tuần hoàn từ 1,00 m đến 1,30 m vào cao trình đỉnh đê và mặt cầu cảng từ Quảng Ninh đến Nghệ An để thích ứng với bão mạnh.
  • Xây dựng bản đồ phân vùng nguy cơ ngập lụt do nước dâng ven bờ: Ủy ban nhân dân các tỉnh ven biển phối hợp với các viện nghiên cứu chuyên ngành hoàn thành việc lập bản đồ rủi ro ngập lụt tỷ lệ 1:25.000 cho 8 phân khu kinh tế trọng điểm duyên hải trước năm 2028.

Đố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ị tham khảo thiết thực cho 4 nhóm đối tượng chuyên môn:

  • Kỹ sư thiết kế công trình thủy công, đê điều và cảng biển: Sử dụng trực tiếp các phương trình tương quan và bảng tra cứu mực nước dâng rút cực trị theo từng hướng gió để xác định chính xác cao trình đỉnh đê, độ sâu nạo vét luồng tàu và tải trọng sóng tác động lên thân công trình.
  • Chuyên viên dự báo khí tượng thủy văn và cảnh báo thiên tai: Khai thác các hệ số tương quan giữa tốc độ gió, áp suất khí quyển và mực nước phi tuần hoàn nhằm nâng cao độ tin cậy của các bản tin cảnh báo nước dâng do bão và triều cường.
  • Nhà quản lý quy hoạch đô thị và không gian biển duyên hải: Ứng dụng các kịch bản phân bố mực nước dâng để hoạch định ranh giới hành lang bảo vệ bờ biển, quy hoạch khu đô thị lấn biển và bố trí hệ thống thoát nước mưa chống ngập úng.
  • Học viên, nghiên cứu sinh ngành Hải dương học và Thủy khí động lực: Tiếp cận một tài liệu mẫu mực về quy trình kết hợp giữa phân tích điều hòa thủy triều, kỹ thuật xây dựng lưới phi cấu trúc trong mô hình MIKE 21 FM và xử lý số liệu quan trắc hải văn.

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

Dao động mực nước phi tuần hoàn khác gì so với dao động thủy triều thông thường?
Dao động thủy triều mang tính tuần hoàn nghiêm ngặt do lực hấp dẫn thiên văn giữa Mặt Trăng, Mặt Trời và Trái Đất tạo ra với chu kỳ nửa ngày hoặc ngày. Ngược lại, dao động phi tuần hoàn hình thành do các yếu tố khí tượng bất định như ma sát gió, biến thiên áp suất và dòng chảy lũ, làm mực nước biển dâng cao hoặc hạ thấp bất thường ngoài dự tính thiên văn.

Tại sao gió hướng Đông Nam lại gây ra độ dâng nước lớn nhất tại khu vực bờ Tây vịnh Bắc Bộ?
Gió Đông Nam thổi vuông góc trực tiếp vào dải bờ biển cong hình cánh cung của vịnh Bắc Bộ, kết hợp với hiệu ứng đẩy dồn khối nước từ vùng biển sâu phía ngoài vào vùng nước nông ven bờ bị chặn bởi lục địa. Tại trạm Cửa Ông và Mũi Ngọc, tốc độ gió 25 m/s hướng Đông Nam đã tạo ra mức dâng cực đại lên tới 1,30 m.

Hiệu ứng nghịch áp khí quyển tác động đến mực nước biển ven bờ theo tỷ lệ nào?
Theo quy luật tĩnh học, khi áp suất khí quyển giảm 1 mbar (hPa) ở tâm bão hoặc áp thấp nhiệt đới, mực nước biển sẽ tự động dâng lên xấp xỉ 10 mm (1 cm). Khi một cơn bão có khí áp giảm sâu 40 mbar tràn qua, riêng thành phần nghịch áp tĩnh đã có thể đóng góp khoảng 40 cm vào tổng độ dâng của mực nước biển.

Mô hình MIKE 21 FM có ưu thế vượt trội nào trong tính toán thủy văn vịnh Bắc Bộ?
Mô hình MIKE 21 FM sử dụng phương pháp thể tích hữu hạn trên lưới phi cấu trúc tam giác, cho phép linh hoạt làm mịn mắt lưới tại các vùng bờ biển có địa hình phức tạp, nhiều đảo nhỏ như vịnh Hạ Long và Bái Tử Long, giúp tính toán chính xác dòng chảy và dao động mực nước mà không gặp lỗi số trị.

Làm thế nào để tách được thành phần mực nước phi tuần hoàn từ chuỗi dữ liệu thực đo?
Người nghiên cứu sử dụng phương pháp phân tích điều hòa để tính toán bộ hằng số điều hòa biên độ và pha của các sóng triều, sau đó tính toán mực nước thủy triều lý thuyết cho từng giờ. Mực nước phi tuần hoàn được xác định bằng cách lấy giá trị mực nước quan trắc thực tế trừ đi giá trị mực nước thủy triều tính toán tại cùng thời điểm.

Kết luận

  • Luận văn đã làm sáng tỏ cơ chế tác động của trường gió và biến thiên áp suất khí quyển đến dao động dâng, rút mực nước phi tuần hoàn tại bờ Tây vịnh Bắc Bộ.
  • Xác lập định lượng hướng gió Đông Nam và hướng Đông là hai nguyên nhân chính gây nước dâng nguy hiểm, với giá trị dâng cực đại lên tới 1,30 m tại Cửa Ông và 1,09 m tại Diễn Châu ở cấp gió 25 m/s.
  • Xây dựng thành công bộ 4.014 nút lưới và 6.830 phần tử tam giác trên mô hình MIKE 21 FM, kiểm định khớp pha và biên độ với chuỗi thực đo tại hai trạm chuẩn Hòn Dáu và Hòn Ngư.
  • Cung cấp hệ thống phương trình hồi quy tuyến tính có độ tin cậy cao với hệ số xác định đạt từ 0,95 đến 0,99, phục vụ trực tiếp công tác thiết kế công trình biển và quy hoạch đê điều.
  • Đề xuất kế hoạch hoàn thiện hệ thống quan trắc tự động và tích hợp mô hình cảnh báo sớm thời gian thực trong lộ trình 2026-2030 nhằm chủ động phòng chống thiên tai biển.

Hãy áp dụng ngay các giải pháp kỹ thuật và bộ dữ liệu thủy động lực học từ luận văn vào công tác tính toán thiết kế công trình thủy ven bờ để nâng cao tính bền vững và bảo vệ an toàn cho các dải đô thị duyên hải Việt Nam.