Tổng quan nghiên cứu

Khu vực cửa Đề Gi và đầm Nước Ngọt nằm tại ranh giới hai xã Mỹ Chánh (huyện Phù Mỹ) và Cát Minh (huyện Phù Cát), tỉnh Bình Định, giữ vai trò đặc biệt quan trọng trong phát triển kinh tế biển địa phương. Đây là cửa thoát lũ duy nhất của lưu vực sông La Tinh rộng 719 km², đồng thời kết nối không gian trao đổi nước cho đầm Đề Gi có diện tích mặt nước khoảng 15 km² tại tọa độ 14°09'N và 109°10'E. Vùng cửa biển này vừa là tuyến giao thông thủy huyết mạch phục vụ hàng trăm tàu cá cập cảng cá Đề Gi, vừa là nơi neo đậu trú tránh bão an toàn cho tàu thuyền khu vực Nam Trung Bộ.

Tuy nhiên, tình trạng bồi lấp luồng tàu diễn ra liên tục và ngày càng phức tạp. Mặc dù vào năm 2002, công trình đê chắn sóng dạng chữ I dài 480 m, cao trình 3,5 m với chân đế rộng 35 - 40 m đã được xây dựng vươn ra độ sâu gần 5 m, hiện tượng bồi tụ cát vẫn chưa được xử lý triệt để. Các doi cát ngầm dài 100 - 200 m liên tục hình thành, khiến độ sâu luồng tàu suy giảm mạnh từ mức thiết kế 4,0 - 5,0 m xuống dưới 1,5 m trong nhiều thời điểm, làm gia tăng nguy cơ mắc cạn cho phương tiện thủy.

Mục tiêu nghiên cứu là làm rõ nguyên nhân khiến công trình kè hiện hữu chưa phát huy tối đa hiệu quả, từ đó đề xuất các giải pháp công trình và phi công trình tối ưu nhằm chống bồi lấp cửa Đề Gi. Nghiên cứu khai thác chuỗi dữ liệu khí tượng hải văn 25 năm (1988–2012) và kết quả đo đạc thực địa chi tiết giai đoạn 2012–2013. Kết quả luận văn cung cấp cơ sở khoa học để khôi phục độ sâu luồng chạy tàu đạt trên 4,0 m, giảm thiểu hơn 75% khối lượng cát bồi lắng hàng năm và ổn định sinh kế cho hơn 10.000 hộ dân ven đầ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 ứng dụng hệ thống lý thuyết thủy thạch động lực học vùng cửa sông ven biển để mô tả sự tương tác đa chiều giữa sóng, dòng chảy và trường trầm tích. Mô hình thủy động lực dựa trên hệ phương trình Navier-Stokes ba chiều cho chất lỏng không nén trong điều kiện nước nông kết hợp phép xấp xỉ Boussinesq và phương trình cân bằng thủy tĩnh.

Mô hình sóng SWAN (Delft3D-Wave) thế hệ thứ ba giải phương trình cân bằng mật độ tác động phổ sóng (Action Density Wave spectrum), tích hợp các quá trình tăng trưởng sóng do gió, khúc xạ, nhiễu xạ và tiêu tán năng lượng do sóng đổ theo mô hình Battjes & Janssen (1978) kết hợp phổ JONSWAP với hệ số ma sát đáy 0,067.

Quá trình vận chuyển trầm tích lơ lửng được mô tả thông qua phương trình bình lưu - khuếch tán 3 chiều và công thức Partheniades - Krone. Lý thuyết tương tác sóng - dòng chảy của Soulsby được sử dụng để xác định ứng suất trượt đáy tổng hợp, làm rõ cơ chế xáo trộn hạt cát đáy có đường kính trung bình D50 = 250 µm và mật độ 2.650 kg/m³.

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

Về cỡ mẫu và phương pháp thu thập dữ liệu, nghiên cứu đã xây dựng ma trận địa hình đáy tổng hợp từ 5.317 điểm đo độ sâu chi tiết (XYZ). Dữ liệu được số hóa và đồng hóa theo hệ quy chiếu VN2000 từ bốn nguồn chính: Hải đồ Hải quân Mỹ năm 1967 tỷ lệ 1:50.000, Hải đồ Hải quân Nhân dân Việt Nam năm 1991 tỷ lệ 1:100.000, bản đồ số C-map năm 1999 tỷ lệ 1:50.000 và đợt khảo sát đo sâu hồi âm tháng 10/2009 của Viện Hải dương học Nha Trang. Chuỗi số liệu sóng và gió nước sâu 25 năm liên tục (1988–2012) được trích xuất từ mô hình WaveWatch III và dữ liệu tái phân tích NCEP tại tọa độ 14,5°N - 109,5°E với độ phân giải lưới 0,5 độ. Mẫu quan trắc mực nước dùng để hiệu chỉnh mô hình gồm 2.184 giá trị đo liên tục với tần suất 1 giờ/lần từ ngày 04/11/2012 đến 04/02/2013 tại trạm Mỹ Thành (14,15°N - 109,17°E).

Lý do lựa chọn phương pháp phân tích: Luận văn sử dụng bộ phần mềm Delft3D với tính năng liên kết đồng thời (online coupling) giữa ba mô-đun Delft3D-Flow, Delft3D-Wave và Delft3D-Sed. Miền tính toán rộng 2.000 km² được phủ bởi hệ lưới cong trực giao 76 x 116 ô, kích thước ô lưới mịn từ 42 m đến 124 m, bước thời gian tích phân 5 phút. Phương pháp này mô phỏng chân thực tính chất nhật triều không đều thông qua 8 hằng số điều hòa chính (M2, S2, K1, O1, N2, K2, P1, Q1) và phản ánh chính xác tương tác phi tuyến giữa dòng chảy triều và sóng biển.

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

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

Thứ nhất, chế độ sóng tại khu vực mang tính phân mùa rõ rệt, trong đó sóng gió mùa Đông Bắc (NE) chiếm ưu thế áp đảo với tần suất xuất hiện trên 50%, chiều cao sóng gây biến đổi đáy Hmor = 2,12 m và chu kỳ đỉnh phổ Tp = 8,6 s. Sóng hướng Đông Nam (SE) chiếm 19% với Hmor = 0,85 m và sóng hướng Đông (E) chiếm khoảng 15%, trong khi các hướng sóng khác chỉ chiếm xấp xỉ 15% tổng tần suất năm.

Thứ hai, công trình đê chắn sóng 480 m hiện hữu sau khi xây dựng đã làm biến dạng cục bộ trường dòng chảy. Vận tốc dòng chảy trong lòng eo cửa khi triều rút đạt mức cực đại 0,8 - 1,2 m/s, tạo xói sâu cục bộ ở luồng chính nhưng lại làm giảm vận tốc dòng chảy xuống dưới 0,2 m/s ở vùng khuất sau kè phía Nam, tạo điều kiện thuận lợi cho bùn cát bồi tụ thành các đụn cát dài 100 - 200 m và rộng 5 - 15 m.

Thứ ba, nồng độ trầm tích lơ lửng trong mùa gió Đông Bắc đạt giá trị rất cao, dao động từ 0,35 kg/m³ đến 0,50 kg/m³ tại khu vực trước cửa và giảm dần vào phía trong đầm ở mức 0,05 - 0,10 kg/m³. Mức nồng độ này cao hơn từ 150% đến 200% so với mùa gió Tây Nam (chỉ đạt 0,05 - 0,15 kg/m³), dẫn đến tốc độ bồi lắng đáy luồng tàu ước tính khoảng 0,3 - 0,6 m/năm.

Thứ tư, quá trình hiệu chỉnh và kiểm định mô hình thủy động lực cho thấy độ tin cậy rất cao. Hệ số tương quan R giữa chuỗi mực nước tính toán và thực đo tại trạm Mỹ Thành trong giai đoạn 05/11/2012 – 25/11/2012 đạt 0,94, với sai số biên độ mực nước trung bình chỉ ở mức 0,1 m.

Thảo luận kết quả

Cơ chế bồi lấp cửa Đề Gi được làm rõ dưới góc độ động lực học: Tuyến kè chữ I dài 480 m vươn ra độ sâu gần 5 m đã chặn được dòng bùn cát ven bờ nông, nhưng lại tạo ra hiệu ứng hội tụ sóng tại đầu kè. Dưới tác động của sóng Đông Bắc có năng lượng lớn, trường sóng bị khúc xạ bao quanh mũi kè, kết hợp với dòng triều dâng có vận tốc lớn cuốn lượng cát bị xáo trộn vòng qua đầu kè rồi đẩy thẳng vào lòng eo cửa.

Dữ liệu nghiên cứu được trực quan hóa sinh động qua hệ thống đồ thị và bảng số liệu chuyên dụng. Biểu đồ hoa sóng 25 năm và hoa gió NCEP minh chứng rõ nét hướng tác động cơ bản của trường lực khí tượng hải văn. Biểu đồ so sánh đường cong biến trình mực nước thực đo và mô phỏng qua 211 giờ kiểm định khẳng định độ chuẩn xác của mô hình toán. Đồng thời, các bản đồ vector trường vận tốc dòng chảy 2 chiều và bản đồ bình đồ phân bố nồng độ bùn cát lơ lửng đã bóc tách rõ nét sự khác biệt giữa hai kịch bản có công trình và không có công trình kè biển.

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

Dựa trên các kết quả tính toán mô phỏng, luận văn đề xuất bốn nhóm giải pháp kỹ thuật và quản lý cụ thể:

  1. Kéo dài và điều chỉnh hướng đê chắn sóng phía Bắc: Tiến hành thi công kéo dài thân đê hiện hữu thêm 120 - 150 m và uốn cong mũi đê về phía Đông Nam với góc nghiêng khoảng 15 - 20 độ, nhằm dẫn hướng dòng bùn cát ra xa cửa lạch và duy trì độ sâu đáy luồng ổn định trên 4,5 m trước năm 2028; đơn vị thực hiện: Sở Nông nghiệp và Phát triển Nông thôn tỉnh Bình Định.
  2. Xây dựng đê mỏ hàn phụ phía bờ Nam: Thiết lập tuyến kè mỏ hàn bằng khối bê tông tự chèn dài 80 - 100 m tại mũi phía Nam cửa Đề Gi nhằm ngăn chặn dòng cát trườn men theo bờ Nam eo cửa, giảm thiểu khoảng 60% lượng cát bồi lấp vào luồng lạch trong giai đoạn 2026–2027; đơn vị thực hiện: Ban Quản lý Dự án Nông nghiệp tỉnh Bình Định.
  3. Nạo vét luồng định kỳ kết hợp tái sử dụng bùn cát: Triển khai nạo vét thông luồng định kỳ hàng năm sau mùa mưa bão với khối lượng duy tu khoảng 50.000 - 80.000 m³/năm, đồng thời sử dụng toàn bộ lượng cát nạo vét đạt chuẩn để bơm nuôi bãi chống xói lở cho đoạn bờ biển xã Mỹ Thành; đơn vị thực hiện: Chi cục Thủy sản và Ban Quản lý Cảng cá Đề Gi.
  4. Thiết lập hệ sinh thái trạm quan trắc hải văn tự động: Lắp đặt 01 trạm quan trắc tự động đo mực nước, sóng và dòng chảy thời gian thực tại eo cửa Đề Gi với chu kỳ truyền tin 15 phút/lần, hoàn thành trước quý IV/2027 nhằm cung cấp số liệu phục vụ vận hành cảnh báo an toàn luồng lạch; đơn vị thực hiện: Đài Khí tượng Thủy văn khu vực Nam Trung Bộ phối hợp với UBND huyện Phù Cát.

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

  1. Cơ quan quản lý nhà nước về tài nguyên, nông nghiệp và giao thông thủy: Cung cấp cơ sở khoa học thực chứng để Sở Nông nghiệp và Phát triển Nông thôn, Sở Giao thông Vận tải tỉnh Bình Định phê duyệt các dự án nâng cấp hạ tầng cảng cá và quy hoạch chỉnh trị bờ biển.
  2. Kỹ sư tư vấn thiết kế công trình thủy công: Cung cấp tài liệu tham khảo chuyên môn về quy trình ứng dụng mô hình Delft3D liên kết sóng - dòng chảy - bùn cát trong thiết kế đê kè chắn sóng vùng cửa sông miền Trung.
  3. Nhà nghiên cứu và học viên cao học ngành Hải dương học, Thủy văn: Cung cấp phương pháp luận chi tiết về xử lý dữ liệu địa hình phức hợp 5.317 điểm đo và đồng hóa chuỗi sóng nước sâu 25 năm phục vụ mô phỏng biến động đáy.
  4. Ban quản lý cảng cá và nghiệp đoàn nghề cá địa phương: Giúp nắm bắt quy luật nông sâu luồng tàu theo từng mùa gió và pha triều để xây dựng lịch điều phối tàu thuyền ra vào cảng an toàn, giảm thiểu rủi ro tai nạn hàng hải.

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

Tại sao công trình kè chắn sóng xây dựng năm 2002 chưa thể ngăn chặn triệt để bồi lấp cửa Đề Gi?
Kè chữ I dài 480 m chỉ ngăn được dòng cát ven bờ sát mép nước nhưng chưa bao quát hết đới sóng đổ ngoài độ sâu 5 m. Dưới tác động của sóng Đông Bắc mạnh với Hmor = 2,12 m, bùn cát bị xáo trộn mạnh và cuốn vòng qua đầu kè, sau đó dòng triều dâng tải thẳng lượng cát này vào lòng eo cửa.

Mô hình Delft3D đóng vai trò gì trong việc giải quyết bài toán bồi lấp tại cửa Đề Gi?
Delft3D giải quyết đồng thời sự tương tác phức tạp giữa ba quá trình: thủy động lực học (Delft3D-Flow), trường sóng phổ SWAN (Delft3D-Wave) và vận chuyển bùn cát (Delft3D-Sed) trên lưới cong 76 x 116 ô, giúp mô phỏng chính xác sự thay đổi địa hình đáy theo từng kịch bản công trình.

Mùa gió nào gây ra hiện tượng bồi lắng bùn cát nghiêm trọng nhất tại khu vực?
Gió mùa Đông Bắc (tháng 11 đến tháng 2 năm sau) là thời kỳ bồi lắng mạnh nhất do sóng hướng Đông Bắc chiếm trên 50% tần suất xuất hiện, tạo nồng độ trầm tích lơ lửng cực đại lên đến 0,50 kg/m³, cao gấp 2 đến 3 lần so với mùa gió Tây Nam.

Độ tin cậy của các kết quả mô phỏng trong luận văn được đánh giá như thế nào?
Mô hình toán được hiệu chỉnh và kiểm định độc lập với chuỗi quan trắc mực nước thực đo 2.184 giờ tại trạm Mỹ Thành. Kết quả đạt hệ số tương quan R = 0,94 và sai số biên độ mực nước chỉ khoảng 0,1 m, khẳng định độ tin cậy rất cao của các kịch bản dự báo.

Giải pháp công trình nào được đánh giá là tối ưu nhất để ổn định luồng chạy tàu cửa Đề Gi?
Giải pháp tối ưu là kéo dài kè chắn sóng phía Bắc thêm 120 - 150 m và uốn cong về hướng Đông Nam, kết hợp xây dựng kè mỏ hàn phía Nam dài 80 - 100 m để hướng dòng chảy triều bóc tách bùn cát ra khỏi luồng chính, giữ độ sâu luồng tàu ổn định trên 4,5 m.

Kết luận

  • Luận văn đã xác định rõ nguyên nhân cốt lõi gây bồi lấp cửa Đề Gi là do tác động phối hợp giữa trường sóng Đông Bắc chiếm trên 50% tần suất và dòng triều dâng mang bùn cát vòng qua mũi kè 480 m vào lòng lạch.
  • Thiết lập và kiểm định thành công mô hình Delft3D 3 chiều với hệ số tương quan R = 0,94, tái hiện chính xác chế độ thủy thạch động lực trên miền tính rộng 2.000 km² của đầm Đề Gi.
  • Định lượng chi tiết nồng độ trầm tích lơ lửng cực đại đạt 0,50 kg/m³ và tốc độ bồi lắng luồng tàu từ 0,3 - 0,6 m/năm, làm cơ sở khoa học tin cậy cho công tác chỉnh trị luồng.
  • Đề xuất đồng bộ gói giải pháp kéo dài đê chắn sóng thêm 120 - 150 m uốn góc 15 - 20 độ, xây dựng kè mỏ hàn bờ Nam và nạo vét duy tu định kỳ 50.000 - 80.000 m³/năm.
  • Công trình đóng góp nền tảng kỹ thuật vững chắc để tỉnh Bình Định triển khai dự án nâng cấp luồng tàu và cảng cá Đề Gi giai đoạn 2026–2028.

Các cơ quan quản lý, viện nghiên cứu và doanh nghiệp quan tâm có thể ứng dụng toàn bộ khung tính toán và bộ thông số mô hình của luận văn để lập báo cáo nghiên cứu khả thi cho các dự án chỉnh trị cửa sông ven biển tương tự.