Tổng quan nghiên cứu

Thảm rừng ngập mặn ven biển đóng vai trò như một bức tường xanh tự nhiên bảo vệ dải bờ trước tác động tàn phá của bão lũ và nước dâng. Thực tế ghi nhận vào tháng 7 năm 1996, khi cơn bão số 2 với sức gió đạt 103 đến 117 km/h đổ bộ vào khu vực Thái Thụy thuộc tỉnh Thái Bình, toàn bộ tuyến đê biển và đầm nuôi trồng thủy sản vẫn được giữ vững nhờ dải rừng ngập mặn che chắn. Ngược lại, tại huyện Tiền Hải liền kề, do dải rừng bị suy thoái và chặt phá nên hệ thống bờ đầm bị sóng đánh phá vỡ và xói lở nghiêm trọng. Đến năm 2005, cơn bão số 7 tiếp tục gây sóng lớn tại cửa sông Trà Lý, làm sạt lở hơn 650 m đê biển tại khu vực không có rừng ngập mặn ở thôn Tân Bồi, xã Thái Đô. Trong khi đó, các đoạn đê được thảm thực vật che chở đều an toàn tuyệt đối nhờ khả năng làm giảm 50% đến 75% chiều cao sóng tới.

Luận văn thạc sĩ khoa học chuyên ngành Hải dương học tập trung vào mục tiêu xây dựng công cụ số trị nhằm mô phỏng và lượng hóa chi tiết quá trình lan truyền cũng như tiêu tán năng lượng sóng khi đi qua thảm cây ngập mặn. Nghiên cứu được triển khai áp dụng cụ thể tại vùng ven biển cửa sông Trà Lý thuộc tỉnh Thái Bình trên bãi bùn có độ dốc khoảng 1:2000. Phạm vi dữ liệu phục vụ nghiên cứu bao gồm chuỗi số liệu quan trắc khí tượng hải văn nhiều năm từ năm 1960 đến năm 1994 kết hợp dữ liệu đo đạc thực địa chi tiết. Kết quả của công trình mang lại giá trị thực tiễn to lớn, hỗ trợ các nhà quản lý tính toán chính xác biên độ sóng giảm thiểu tới hơn 60% khi tấn công vào chân công trình kè đê, từ đó tối ưu hóa chi phí đầu tư xây dựng các công trình chỉnh trị bờ 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 và động lực học biển hiện đại, tích hợp các lý thuyết sóng cơ bản cùng mô hình số trị chuyên sâu:

  • Lý thuyết sóng tuyến tính: Áp dụng hệ phương trình sóng Airy (1845) và phương trình thế vận tốc hai chiều theo nghiên cứu của Young (1999), mô tả chuyển động tuần hoàn của hạt nước trong môi trường chất lỏng không nén được và bỏ qua độ nhớt ở vùng nước có độ sâu thay đổi.
  • Lý thuyết tiêu tán năng lượng do lực cản công trình và thực vật: Dựa trên phương trình lực Morison (1950) tính toán lực kéo và lực quán tính tác dụng lên thân cây hình trụ. Lý thuyết này được kế thừa và phát triển bởi Dalrymple (1984) và hiệu chỉnh mở rộng bởi Mendez và Losada (2004) để xác định sự suy giảm năng lượng sóng ngẫu nhiên qua các vùng bãi triều có cây ngập nước.
  • Mô hình phổ sóng thế hệ thứ ba SWAN: Mô hình Simulating Waves Nearshore giải phương trình cân bằng tác động phổ hai chiều trên miền không gian địa lý và phổ tần số. Năng lượng sóng được bảo toàn và biến đổi thông qua hàm nguồn tổng quát bao gồm: năng lượng truyền từ gió, tiêu tán do sóng bạc đầu, ma sát đáy theo công thức Collins (1972) hoặc Van Rijn (1989), sóng vỡ do địa hình nông và tương tác phi tuyến sóng ba bậc, bốn bậc.
  • Các khái niệm then chốt: Khái niệm chiều cao sóng hữu hiệu, chu kỳ sóng đỉnh phổ, mật độ thực vật trên một đơn vị diện tích và hệ số lực cản biểu kiến của thân cây ngập mặn.

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

Nguồn dữ liệu đầu vào của luận văn được tổng hợp từ bản đồ địa hình đáy biển khu vực ven bờ Thái Bình, bình đồ độ sâu chi tiết và chuỗi quan trắc khí tượng thủy văn trạm Hòn Dáu trong suốt 34 năm (giai đoạn 1960 - 1994). Số liệu gió bão cực đoan được trích xuất từ các cơn bão lịch sử đổ bộ trực tiếp vào khu vực đồng bằng Bắc Bộ.

Phương pháp chọn mẫu và thiết lập miền tính được thực hiện theo cấu trúc lưới lồng hai cấp độ nhằm tối ưu hóa độ chính xác:

  • Lưới tính thô: Bao phủ toàn bộ vùng biển ven bờ với kích thước 32,56 km dọc theo bờ biển và 30,8 km theo phương vuông góc ra khơi. Miền tính thô bao gồm 297 x 281 điểm nút mạng với kích thước mỗi ô lưới là 110 m x 110 m (tương đương quy mô mẫu 83.457 điểm tính).
  • Lưới tính tinh: Tập trung chi tiết vào khu vực rừng ngập mặn cửa Trà Lý với kích thước 5,72 km dọc bờ và 5,368 km ra khơi, chia thành 245 x 261 điểm nút với kích thước ô lưới mịn 22 m x 22 m (tổng cộng 63.945 điểm tính).

Lý do lựa chọn phương pháp mô hình hóa số trị phân lớp là vì rừng ngập mặn có cấu trúc hình học phức tạp thay đổi rõ rệt theo phương thẳng đứng (từ hệ rễ chống, thân chính đến cành lá). Việc chia cắt thân cây thành nhiều lớp độ sâu và tích phân tham số cản vào mô hình SWAN giúp khắc phục triệt để nhược điểm của các phương pháp ma sát đáy tương đương truyền thống, mang lại độ tin cậy vượt trội.

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

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

Nghiên cứu đã mang lại các kết quả định lượng rõ nét về đặc trưng lan truyền sóng tại khu vực cửa Trà Lý:

  • Khả năng tiêu tán chiều cao sóng hữu hiệu: Qua mô phỏng số trị, dải rừng ngập mặn tại xã Thái Đô với chiều rộng trung bình từ 400 m đến 600 m có khả năng làm suy giảm từ 45% đến 70% chiều cao sóng tới. Trong điều kiện sóng ngoài khơi cao 0,8 m truyền qua dải cây có mật độ trung bình khoảng 1000 thân cây nhỏ trên một mét vuông, chiều cao sóng khi vào gần chân đê chỉ còn dao động từ 0,2 m đến 0,28 m.
  • Quy luật phân bố sóng theo mùa: Dữ liệu phân tích nhiều năm cho thấy trong mùa đông (từ tháng 10 đến tháng 3), sóng hướng Đông chiếm ưu thế với tần suất xuất hiện dao động từ 30% đến 56,3%. Vào các tháng mùa hè (tháng 6 và tháng 7), sóng hướng Đông Nam và Nam chiếm tỷ lệ từ 38% đến 70%, tạo ra các đợt sóng trung bình từ 0,6 m đến 0,8 m trực tiếp uy hiếp vùng cửa sông.
  • Tác động của cấu trúc phân lớp thẳng đứng: So sánh giữa mô phỏng thân cây đồng nhất và mô phỏng phân tầng rễ - thân cho thấy việc đưa cấu trúc thẳng đứng vào tính toán giúp tăng độ chính xác của trường tiêu tán năng lượng lên thêm 18,5% so với mô hình ma sát đáy đơn lẻ.

Thảo luận kết quả

Các phát hiện của luận văn chứng minh rằng cơ chế tiêu tán năng lượng sóng chủ đạo trong rừng ngập mặn bắt nguồn từ lực cản thủy động lực do thân và rễ cây tạo ra, vượt trội gấp 5 đến 8 lần so với lực ma sát đáy thuần túy. Hiện tượng này hoàn toàn phù hợp với các quan sát thực địa của các nhà nghiên cứu quốc tế tại khu vực biển Đông Nam Á sau thảm họa sóng thần năm 2004.

Dữ liệu tính toán từ mô hình có thể được biểu diễn trực quan qua biểu đồ không gian phân bố chiều cao sóng hữu hiệu hai chiều và bảng thống kê suy giảm năng lượng sóng theo từng mặt cắt trắc dọc từ ngoài khơi vào bờ. Khi sóng truyền từ vùng nước sâu vào bãi triều, biểu đồ cho thấy đường năng lượng sóng giảm đột ngột tại ranh giới ngoài cùng của dải rừng ngập mặn. Điều này khẳng định rằng vành đai ngoài cùng của rừng đóng vai trò hấp thụ phần lớn động năng ban đầu của sóng biển.

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

Dựa trên kết quả nghiên cứu định lượng, luận văn đưa ra 4 nhóm giải pháp công trình và phi công trình cụ thể:

  1. Thiết lập vành đai rừng phòng hộ ven biển chuẩn hóa: Ủy ban nhân dân tỉnh Thái Bình và Sở Nông nghiệp và Phát triển nông thôn cần quy hoạch duy trì chiều rộng vành đai rừng ngập mặn tối thiểu từ 300 m đến 500 m phía trước các tuyến đê xung yếu, đặt mục tiêu duy trì tỷ lệ che phủ thảm thực vật trên 80% tại khu vực cửa Trà Lý trong giai đoạn 2026 - 2030.
  2. Ứng dụng mô hình SWAN cải tiến trong thiết kế đê biển: Đề nghị các viện thiết kế công trình thủy lợi tích hợp mô-đun tiêu tán sóng do thực vật vào quy trình tính toán tải trọng sóng tác động lên thân đê, giúp giảm kích thước tiết diện công trình và tiết kiệm từ 15% đến 20% chi phí xây dựng kè bảo vệ.
  3. Hiện đại hóa hệ thống quan trắc sóng ven bờ: Bộ Tài nguyên và Môi trường cần lắp đặt bổ sung 02 trạm đo sóng tự động tại vùng bãi triều cửa sông Thái Bình và Trà Lý trước năm 2028 nhằm liên tục cập nhật chuỗi số liệu hiệu chỉnh tham số lực cản biểu kiến cho mô hình số.
  4. Xây dựng cẩm nang kỹ thuật trồng cây chắn sóng đa tầng: Các cơ quan khuyến lâm địa phương cần triển khai kế hoạch phối hợp các loài cây như bần chua, trang, vẹt dù và sú theo từng dải độ sâu thích hợp, tạo nên cấu trúc rừng nhiều tầng tán nhằm tối đa hóa khả năng triệt tiêu sóng bão từ nay đến năm 2030.

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

Tài liệu này là nguồn tham khảo chuyên sâu và hữu ích cho các nhóm đối tượng sau:

  • Kỹ sư công trình thủy lợi và đê điều ven biển: Ứng dụng phương pháp tính toán tải trọng sóng thực tế đã qua tiêu tán để thiết kế tối ưu mái đê và chân kè biển.
  • Nhà nghiên cứu và giảng viên chuyên ngành Hải dương học: Sử dụng thuật toán tích phân lực cản thực vật trong mô hình SWAN làm tài liệu giảng dạy và mở rộng nghiên cứu mô phỏng vận chuyển trầm tích.
  • Cán bộ quản lý phòng chống thiên tai tại các tỉnh duyên hải: Tham khảo các kịch bản suy giảm sóng bão để hoạch định chính sách bảo tồn rừng ngập mặn và lập phương án ứng phó ngập lụt ven bờ.
  • Học viên cao học và sinh viên ngành Khoa học Trái đất: Tiếp cận phương pháp thiết lập lưới lồng số trị chi tiết và xử lý chuỗi số liệu khí tượng hải văn phức tạp trong đồ án tốt nghiệp.

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

1. Rừng ngập mặn có khả năng làm giảm bao nhiêu phần trăm chiều cao sóng biển?

Theo kết quả mô phỏng số trị tại cửa Trà Lý, một dải rừng ngập mặn trưởng thành có bề rộng từ 400 m đến 600 m có thể làm giảm từ 45% đến 70% chiều cao sóng bão khi truyền từ ngoài khơi vào bờ. Khi chiều cao sóng ngoài khơi đạt 0,8 m, sóng truyền tới chân đê phía sau rừng chỉ còn dưới 0,28 m, giúp triệt tiêu phần lớn năng lượng tàn phá của sóng biển.

2. Mô hình SWAN mô phỏng quá trình tiêu tán sóng qua rừng ngập mặn bằng cách nào?

Mô hình SWAN bổ sung hàm tiêu tán năng lượng do thực vật dựa trên công thức lực cản Morison cải tiến của Dalrymple kết hợp Mendez và Losada. Thuật toán phân chia thân cây thành các tầng lớp thẳng đứng với các tham số hình học như đường kính thân, chiều cao cây và mật độ phân bố trên một mét vuông, sau đó tính toán công tiêu tán trực tiếp của hạt nước trong chu kỳ dao động sóng.

3. Vì sao cần chia cắt cấu trúc cây ngập mặn thành nhiều lớp theo phương thẳng đứng?

Trong tự nhiên, hình thái cây ngập mặn không đồng nhất từ gốc đến ngọn; tầng rễ chống phía dưới có mật độ rất dày đặc với đường kính nhỏ, trong khi thân chính và tán lá phía trên có kích thước lớn hơn nhưng thưa hơn. Việc mô hình hóa phân lớp giúp phản ánh chính xác phân bố vận tốc hạt nước theo độ sâu, làm tăng độ chuẩn xác của kết quả tính toán thêm 18,5% so với việc coi cây là một cọc trụ trơn đơn giản.

4. Chế độ sóng tại khu vực cửa sông Trà Lý - Thái Bình có đặc điểm gì nổi bật?

Khu vực cửa Trà Lý chịu chi phối mạnh mẽ bởi hai mùa gió chính: mùa đông có sóng hướng Đông chiếm tần suất từ 30% đến 56,3%, trong khi mùa hè xuất hiện sóng hướng Đông Nam, Nam và Tây Nam với tần suất từ 38% đến 70%. Chiều cao sóng trung bình mùa bão dao động từ 0,6 m đến 0,8 m kết hợp biên độ triều lớn, tạo ra áp lực xói lở đáng kể lên dải bờ biển nếu không có rừng che chở.

5. Nghiên cứu này đóng góp giải pháp gì cho bài toán kinh tế trong xây dựng đê biển?

Nghiên cứu chứng minh rằng việc bảo tồn và trồng mới dải rừng chắn sóng phía trước công trình đê biển cho phép giảm thiểu đáng kể chiều cao lưu không và độ dày lớp đá lát mái kè. Giải pháp công trình mềm kết hợp rừng ngập mặn này giúp các địa phương tiết kiệm từ 15% đến 20% tổng mức vốn đầu tư xây dựng đê kè, đồng thời nâng cao tuổi thọ công trình trước các biến cố thiên tai cực đoan.

Kết luận

  • Tổng kết năng lực bảo vệ: Rừng ngập mặn là công trình bảo vệ bờ biển tự nhiên đặc biệt hiệu quả, triệt tiêu đến 70% chiều cao sóng và bảo vệ an toàn hệ thống đê điều trước các cơn bão có sức gió trên 100 km/h.
  • Đóng góp học thuật: Luận văn đã tích hợp thành công mô-đun lực cản thực vật đa tầng vào mô hình sóng SWAN thế hệ thứ ba, giải quyết bài toán tương tác phức tạp giữa thủy động lực học sóng và sinh thái ven bờ.
  • Giá trị thực tiễn tại Thái Bình: Cung cấp bộ thông số tính toán tin cậy trên hệ lưới lồng chi tiết 22 m x 22 m cho vùng cửa sông Trà Lý, hỗ trợ trực tiếp công tác quy hoạch dải ven biển địa phương.
  • Kế hoạch triển khai tiếp theo: Tiếp tục mở rộng mô hình cho các vùng bờ biển Nam Bộ và tích hợp bài toán vận chuyển bùn cát gây bồi xói trong giai đoạn 2026 - 2028.
  • Lời kêu gọi hành động: Các cơ quan quản lý và chuyên gia quy hoạch hãy áp dụng ngay phương pháp tính toán số trị thủy văn kết hợp giải pháp bảo tồn hệ sinh thái rừng ngập mặn để kiến tạo dải bờ biển Việt Nam an toàn và phát triển bền vững.