Tổng quan nghiên cứu
Dải ven biển châu thổ sông Hồng có tổng chiều dài bờ biển 135 km, trải dài qua địa bàn 4 tỉnh gồm Hải Phòng, Thái Bình, Nam Định và Ninh Bình. Với diện tích tự nhiên khoảng 16.000 km², đây là khu vực có mật độ dân cư tập trung đông đúc bậc nhất cả nước, đạt xấp xỉ 1.700 người/km². Do đặc điểm tự nhiên là vùng đất thấp trũng, toàn bộ đời sống kinh tế và sinh kế của hàng triệu người dân phụ thuộc trực tiếp vào mức độ an toàn của hệ thống đê biển. Tuy nhiên, châu thổ sông Hồng là khu vực chịu tần suất bão đổ bộ rất cao, chiếm tới 28% tổng số các cơn bão tại Việt Nam.
Trong lịch sử phát triển, vùng châu thổ sông Hồng tồn tại đồng thời hai xu thế đối lập: xu thế bồi tụ lấn biển tại các khu vực cửa Ba Lạt, cửa Đáy với tốc độ bồi đắp trên 500 ha mỗi năm, song hành với hiện tượng xói lở bãi biển và hạ thấp nền địa hình cực kỳ nghiêm trọng tại dải bờ biển huyện Hải Hậu nằm giữa hai cửa sông. Số liệu đo đạc nhiều năm chỉ ra rằng cao trình bãi biển khu vực này bị bóc mòn liên tục từ 2 cm đến 4 cm mỗi năm do hiện tượng xói lở cấu trúc. Trước khi có hệ thống đê biển kiên cố, tốc độ xâm thực lấn sâu vào đất liền đạt từ 5 m đến 10 m mỗi năm, phá hủy hoàn toàn ít nhất 2 tuyến đê biển kể từ năm 1960.
Mục tiêu cốt lõi của nghiên cứu là xây dựng và hiệu chỉnh hệ thống mô hình số động lực SBEACH trong bộ phần mềm công nghệ ven bờ CEDAS, nhằm tính toán chính xác trường sóng ngẫu nhiên, cơ chế vận chuyển bùn cát ngang bờ và mức độ xói lở bãi biển, chân đê trong điều kiện bão cực đoan tại thị trấn Hải Thịnh. Luận văn tập trung đánh giá mức độ an toàn của hệ thống đê biển hiện hữu được thiết kế ở cao trình đỉnh từ +5,0 m đến +5,5 m dưới tác động của các cơn bão có chu kỳ lặp 50 năm và 100 năm. Kết quả nghiên cứu mang lại ý nghĩa thực tiễn to lớn, cung cấp cơ sở khoa học định lượng giúp các cơ quan quản lý giảm thiểu 30% đến 50% nguy cơ sạt lở đê biển trong mùa mưa bão và tối ưu hóa hàng trăm tỷ đồng ngân sách đầu tư công trình bảo vệ bờ.
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 lý thuyết mặt cắt cân bằng động (Equilibrium Beach Profile - EBP) và lý thuyết sóng ngẫu nhiên phi tuyến trong đới sóng nông. Khung lý thuyết trọng tâm bao gồm:
- Lý thuyết mặt cắt cân bằng động (EBP): Được xây dựng dựa trên nguyên lý tiêu tán năng lượng sóng của Dean và phát triển mở rộng bởi Larson và Kraus. Khi tác động của trường sóng và mực nước duy trì liên tục, mặt cắt đáy biển sẽ tự điều chỉnh về trạng thái cân bằng để tiêu tán tối đa năng lượng sóng. Vận chuyển bùn cát ngang bờ xuất hiện do sự chênh lệch giữa mức tiêu tán năng lượng sóng thực tế và mức tiêu tán năng lượng sóng cân bằng trên một đơn vị thể tích nước.
- Lý thuyết trường sóng ngẫu nhiên phân bố Rayleigh: Mô tả sự lan truyền sóng từ vùng nước sâu vào vùng nước nông, tính toán sự suy giảm thông lượng năng lượng sóng, sự biến dạng sóng và quá trình sóng vỡ trong đới ven bờ.
- Cơ chế vận chuyển bùn cát ngang bờ và sạt lở (Avalanching): Phân định ranh giới giữa hai quá trình bùn cát di chuyển ra khơi gây xói lở trong bão và bùn cát bồi tụ vào bờ trong điều kiện bình thường thông qua tham số thời gian sa lắng. Thuật toán sạt lở tự động tái phân bố bùn cát giữa các ô tính lân cận khi độ dốc cục bộ vượt quá góc dốc tới hạn của vật liệu đáy.
- Các khái niệm chuyên ngành then chốt: Sóng đổ bổ nhào (plunging breaker), sóng tràn (overwash), sóng leo trên mái dốc (wave runup), ứng suất bức xạ sóng vuông góc với đường bờ và gradient lưu lượng bùn cát dọc bờ.
Phương pháp nghiên cứu
Luận văn kết hợp chặt chẽ giữa phương pháp mô phỏng số trị trên máy tính và phương pháp phân tích, kiểm định qua số liệu thực nghiệm đo đạc hiện trường:
- Cỡ mẫu và nguồn dữ liệu nghiên cứu: Bộ dữ liệu bao gồm 3 chuỗi thí nghiệm quy mô tự nhiên từ bể sóng siêu lớn SUPERTANK tại Đại học Bang Oregon (bể có chiều dài 104 m, rộng 3,7 m, sâu 4,6 m với 600 m³ cát mẫu đường kính hạt trung bình 0,22 mm); số liệu đo đạc địa hình mặt cắt bãi biển vuông góc với bờ tại trạm Hải Thịnh trước cơn bão ngày 20/9/2005 và sau cơn bão ngày 29/9/2005; chuỗi dữ liệu mực nước triều tại trạm Hòn Dấu và trường sóng biển sâu từ mô hình WAVEWATCH-III với độ cao sóng hữu hiệu 3,5 m và chu kỳ sóng 10 s trong bão Damrey.
- Phương pháp phân tích và mô phỏng: Ứng dụng mô hình SBEACH tích hợp trong hệ thống CEDAS. Lưới tính toán không gian được thiết lập linh hoạt với bước lưới chi tiết từ 0,5 m đến 1,5 m ở vùng sát bờ, 1,5 m đến 3,0 m ở đới sóng vỡ và mở rộng từ 5,0 m đến 15,0 m ở vùng ngoài khơi, tương ứng bước thời gian tích phân từ 5 đến 10 phút. Lý do lựa chọn mô hình SBEACH là khả năng mô phỏng chi tiết các quá trình động lực học phi điều hòa, sóng tràn và xói lở chân đê trong thời gian ngắn của cơn bão mà các mô hình cân bằng tĩnh thông thường không xử lý được.
- Phương pháp đánh giá độ tin cậy: Sử dụng hai chỉ số thống kê định lượng gồm sai số trung bình bình phương (RMS) và chỉ số số dư tương đối (Res) để tối ưu hóa bộ tham số bùn cát thực nghiệm, hệ số vận chuyển do độ dốc bãi và độ sâu biên đới sóng tràn.
Kết quả nghiên cứu và thảo luận
Những phát hiện chính
Mô hình hóa số trị và phân tích số liệu đo đạc đã mang lại những phát hiện quan trọng về động lực học ven bờ tại khu vực Hải Thịnh:
- Hiệu chuẩn thành công bộ tham số mô hình SBEACH: Quá trình kiểm định qua 3 bài toán mẫu SUPERTANK và cơn bão thực tế Damrey (9/2005) cho thấy bộ tham số hiệu chỉnh đạt độ chính xác rất cao. Sai số trung bình bình phương RMS giảm mạnh xuống dưới 0,15 m và hệ số số dư Res đạt giá trị dưới 0,25 (tương đương mức độ sai lệch giảm hơn 40% so với việc sử dụng bộ tham số mặc định của mô hình), phản ánh chính xác quy luật hình thành và dịch chuyển của các bar ngầm chắn sóng.
- Xác định mức độ xói lở bãi biển trong bão Damrey: Trong cơn bão Damrey với sức gió mạnh cấp 10 - 11 (vận tốc gió 28 m/s), bãi biển Hải Thịnh bị xói sâu từ 0,8 m đến 1,2 m ngay tại chân đê. Lượng bùn cát bị dòng chảy ngang bờ bóc tách và vận chuyển ra phía ngoài khơi, tích tụ tạo thành cồn cát ngầm tại vị trí cách chân đê từ 120 m đến 180 m, làm thay đổi hoàn toàn hình thái mặt cắt ban đầu.
- Dự báo xói lở cực đoan theo chu kỳ bão 50 năm và 100 năm: Với kịch bản bão chu kỳ 50 năm (mực nước tổng cộng dâng 1,41 m, sóng biển sâu cao 9,37 m, chu kỳ 11,5 s), cao trình bãi biển bị hạ thấp thêm từ 1,5 m đến 1,8 m. Đối với kịch bản bão chu kỳ 100 năm (mực nước dâng 2,24 m, sóng biển sâu cao 9,69 m, chu kỳ 12,0 s), độ sâu xói lở tại chân đê vượt mức 2,1 m, đẩy vùng sóng vỡ trực diện áp sát vào thân đê.
Thảo luận kết quả
Kết quả nghiên cứu chứng minh rằng cơ chế gây xói lở bãi biển và phá hủy đê điều trong bão hoàn toàn khác biệt với quá trình xói lở cấu trúc theo mùa. Trong điều kiện bão, sự kết hợp giữa nước dâng do gió (chiếm 80% đến 90% tổng lượng nước dâng), áp suất khí quyển giảm và trường sóng có độ dốc lớn làm cho thông lượng năng lượng sóng tiêu tán vượt xa ngưỡng cân bằng động. Điều này kích hoạt dòng vận chuyển bùn cát cực mạnh hướng từ bờ ra khơi.
Dữ liệu tính toán khi được trực quan hóa qua biểu đồ biến đổi mặt cắt địa hình đáy trước và sau bão (Profile Cross-section Plot) cho thấy sự sụt giảm thể tích bãi tập trung mạnh nhất ở đới sóng tràn và đới sóng vỡ. Điểm trũng xói lở sâu nhất xuất hiện ngay tại chân đê, đồng thời bảng so sánh định lượng các phương án hiệu chỉnh tham số chỉ ra rằng nếu không xét đến tham số sạt lở góc dốc thì sai số tính toán có thể tăng lên trên 35%.
So sánh với các nghiên cứu tại vùng bờ biển miền Trung cho thấy bãi biển châu thổ sông Hồng có độ dốc thoải hơn nhiều, nhưng do vật liệu đáy là cát mịn pha bùn nên tốc độ xói lở và mất ổn định chân đê diễn ra nhanh hơn gấp 1,5 lần. Mặc dù hệ thống đê biển hiện nay đã được nâng cấp lên cao trình +5,0 m đến +5,5 m đủ khả năng ngăn sóng tràn đỉnh trong bão cấp 10, nhưng nguy cơ vỡ đê vẫn ở mức báo động do chân đê bị xói rỗng chân, gây trượt mái và sập đổ toàn bộ khối thân đê.
Đề xuất và khuyến nghị
Dựa trên kết quả tính toán định lượng về trường sóng và mức độ xói lở bãi biển, luận văn đề xuất 4 nhóm giải pháp kỹ thuật công trình và quản lý cụ thể:
- Gia cố khẩn cấp kết cấu hộ chân đê biển: Thi công tăng cường hệ thống thảm đá trọng lực bọc lưới thép mạ kẽm kết hợp xếp khối bê tông dị hình Tetrapod tiêu sóng 3 chiều, hạ sâu cao trình đáy hộ chân xuống mức âm 2,0 m đến âm 2,5 m so với mực nước chuẩn lục địa. Mục tiêu đảm bảo an toàn tuyệt đối cho chân đê trước nguy cơ xói sâu 2,1 m trong bão chu kỳ 100 năm, hoàn thành trước quý 3 năm 2027 do Ủy ban nhân dân tỉnh Nam Định và Sở Nông nghiệp và Phát triển Nông thôn làm chủ đầu tư.
- Xây dựng tuyến đê mỏ hàn chữ T giảm sóng xa bờ: Thiết lập hệ thống kè mỏ hàn ngầm chữ T (T-groins) cách mép bờ từ 150 m đến 200 m nhằm triệt tiêu 60% năng lượng sóng biển trước khi truyền vào bờ bãi, đồng thời giữ lại lượng bùn cát bồi tụ dọc bờ. Tiến độ triển khai trong giai đoạn 2026 - 2028 dưới sự chủ trì kỹ thuật của Viện Khoa học Thủy lợi Việt Nam.
- Thực hiện giải pháp nuôi bãi nhân tạo định kỳ: Bơm bổ sung khối lượng cát thạch anh khoảng 150.000 m³ đến 200.000 m³ với kích thước hạt hiệu dụng từ 0,2 mm đến 0,3 mm tại các vị trí xung yếu thuộc xã Hải Lý và thị trấn Hải Thịnh trước tháng 5 hàng năm nhằm tạo vùng đệm tiêu tán sóng tự nhiên, do Ban Quản lý Dự án Đê điều Nam Định điều hành.
- Nâng cấp hệ thống quan trắc và cảnh báo tự động: Lắp đặt mạng lưới trạm đo tự động ghi nhận liên tục độ cao sóng, mực nước triều và định kỳ quét mặt cắt đáy biển bằng công nghệ đo sâu Sonar đa chùm tia 6 tháng một lần, nâng độ chính xác của công tác dự báo rủi ro xói lở lên trên 95% từ năm 2026 do Trung tâm Dự báo Khí tượng Thủy văn Quốc gia quản lý.
Đối tượng nên tham khảo luận văn
Công trình nghiên cứu mang giá trị ứng dụng cao cho nhiều nhóm chuyên gia và cơ quan ban ngành:
- Kỹ sư công trình thủy công và tư vấn thiết kế ven biển: Cung cấp bộ quy trình tính toán thủy thạch động lực học chuẩn hóa trên phần mềm SBEACH/CEDAS, hỗ trợ trực tiếp việc tính toán kết cấu chân khay, mái đê và kè bảo vệ bờ biển chịu được tải trọng bão cấp 11 đến cấp 12.
- Cán bộ quản lý đê điều và phòng chống thiên tai cấp tỉnh: Cung cấp các kịch bản dự báo ngập lụt, xói lở chi tiết cho các chu kỳ lặp 50 năm và 100 năm, làm căn cứ khoa học vững chắc để phê duyệt quy hoạch đê biển và bố trí ngân sách tu bổ hàng năm.
- Giảng viên, nghiên cứu sinh ngành Cơ kỹ thuật và Hải văn học: Nguồn tài liệu tham khảo chuyên sâu về mô hình số tính toán tương tác sóng ngẫu nhiên, cơ chế vận chuyển bùn cát vùng nước nông và phương pháp hiệu chỉnh tham số động lực học ven bờ.
- Các tổ chức quốc tế và cơ quan tài trợ dự án thích ứng biến đổi khí hậu: Sử dụng nguồn dữ liệu thực nghiệm và kết quả mô phỏng tại dải ven bờ châu thổ sông Hồng để đánh giá hiệu quả kinh tế - kỹ thuật của các dự án đê biển kết hợp thích ứng thiên tai.
Câu hỏi thường gặp
Mô hình SBEACH có ưu điểm gì vượt trội trong tính toán biến động bãi biển do bão?
Mô hình SBEACH mô phỏng chính xác quá trình vận chuyển bùn cát ngang bờ dưới tác động đồng thời của trường sóng ngẫu nhiên phân bố Rayleigh, nước dâng bão và triều cường. Mô hình tích hợp cơ chế sạt lở góc dốc và xói lở đụn cát với bước lưới linh hoạt từ 0,5 m đến 15 m, giúp kiểm soát sai số trung bình RMS dưới 0,15 m trong các bài toán thực địa.
Tại sao bờ biển huyện Hải Hậu lại bị xói lở nặng nề dù nằm trong châu thổ sông Hồng?
Bờ biển Hải Hậu nằm ở trung tâm của hai cửa sông lớn nhưng chịu tác động của hiện tượng xói lở cấu trúc do gradient vận chuyển bùn cát dọc bờ bị suy giảm, khiến đáy bãi hạ thấp từ 2 cm đến 4 cm mỗi năm. Khi bão lớn đổ bộ, dòng bùn cát ngang bờ cuốn trôi lớp cát mặt ra khơi làm chân đê bị xói sâu từ 0,8 m đến 1,2 m và sụp đổ nhanh chóng.
Dữ liệu từ bể thí nghiệm SUPERTANK đóng vai trò gì trong việc kiểm định mô hình?
Bể sóng siêu lớn SUPERTANK với chiều dài 104 m và 600 m³ cát mẫu cung cấp dữ liệu thủy động lực học quy mô thực tế không bị ảnh hưởng bởi hiệu ứng tỷ lệ. Việc thử nghiệm thành công trên 3 dạng mặt cắt đại diện (xói đụn cát, biến động bar ngầm, tương tác công trình) đã xác lập độ tin cậy khoa học tuyệt đối cho mô hình trước khi áp dụng tính toán tại bờ biển Hải Thịnh.
Đê biển được nâng cấp lên cao trình +5,0 m đến +5,5 m có đảm bảo an toàn trong bão chu kỳ 100 năm không?
Cao trình đỉnh đê từ +5,0 m đến +5,5 m đủ khả năng ngăn chặn hiện tượng sóng tràn đỉnh trước mực nước dâng tổng cộng 2,24 m, nhưng chưa thể đảm bảo an toàn tuyệt đối. Trường sóng cực đại với độ cao 9,69 m sẽ gây xói sâu chân đê trên 2,1 m, dẫn đến nguy cơ mất ổn định trượt cung tròn thân đê nếu không có công trình hộ chân tăng cường.
Phương pháp tính toán trong luận văn có thể chuyển giao áp dụng cho các vùng biển khác không?
Phương pháp tính toán hoàn toàn có thể chuyển giao cho mọi dải bờ biển tại Việt Nam. Người sử dụng chỉ cần thay thế chuỗi số liệu đầu vào gồm: bình đồ địa hình mặt cắt đo đạc thực tế, kích thước cấp hạt bùn cát đáy và chuỗi số liệu khí tượng hải văn (sóng, triều, gió bão), sau đó tiến hành hiệu chuẩn lại hệ số vận chuyển bùn cát theo dữ liệu đo đạc thực tế tại địa phương.
Kết luận
- Xây dựng và thiết lập thành công hệ thống mô hình số SBEACH trong bộ phần mềm CEDAS, cho phép tính toán chính xác trường sóng ngẫu nhiên và biến động bãi biển trong bão.
- Làm rõ cơ chế xói lở kép tại dải bờ biển châu thổ sông Hồng: xói lở cấu trúc làm mất bãi lâu dài từ 2 cm đến 4 cm mỗi năm kết hợp với xói lở tức thời do dòng bùn cát ngang bờ trong bão.
- Tái hiện sát thực tế cơn bão Damrey năm 2005 với độ sâu xói lở chân đê đạt 1,2 m và xác định vị trí cồn cát ngầm tích tụ cách bờ từ 120 m đến 180 m.
- Dự báo định lượng mức độ hạ thấp cao trình bãi biển từ 1,5 m đến 2,1 m dưới tác động của các cơn bão cực đoan có chu kỳ lặp 50 năm và 100 năm.
- Đề xuất hệ giải pháp công trình tích hợp gồm gia cố chân khay hộ chân đê bằng khối Tetrapod, xây dựng mỏ hàn ngầm chữ T và nuôi bãi cát nhân tạo định kỳ.
Đóng góp lớn nhất của luận văn là cung cấp bộ tham số thực nghiệm đã được kiểm định chuẩn xác cho bài toán động lực học bùn cát ven bờ tại Việt Nam. Định hướng nghiên cứu tiếp theo trong giai đoạn 2026 - 2030 là mở rộng mô hình tính toán không gian 3 chiều và tích hợp các kịch bản nước biển dâng do biến đổi khí hậu toàn cầu. Hãy tải toàn văn tài liệu và bộ tham số mô phỏng để ứng dụng ngay vào các dự án thiết kế, nâng cấp công trình bảo vệ đê biển ven bờ!