Tổng quan nghiên cứu

Việt Nam sở hữu đường bờ biển dài hơn 2.000 km và nằm trực diện trong tâm bão Tây Bắc Thái Bình Dương với tần suất bão đổ bộ trung bình khoảng 6 cơn/năm. Hệ thống đê kè biển đóng vai trò tuyến phòng thủ sinh tử bảo vệ tính mạng, tài sản và cơ sở hạ tầng cho hàng triệu người dân vùng duyên hải. Tuy nhiên, phần lớn các tuyến đê biển hiện nay có cao trình đỉnh đê tương đối hạn chế, phổ biến chỉ từ +4,0 m đến +5,5 m. Trong các đợt bão lớn kèm nước dâng và triều cường, lưu lượng sóng tràn qua đê có thể đạt hàng trăm lít trên giây trên mỗi mét dài. Đây chính là nguyên nhân trực tiếp gây xói lở mái trong và dẫn đến hơn 70% các vụ vỡ đê nghiêm trọng trong lịch sử, điển hình như sự cố đê biển do cơn bão số 7 năm 2005 tại khu vực Bắc Bộ và Bắc Trung Bộ.

Nhằm kiểm soát sóng tràn nhưng vẫn tiết kiệm kinh phí nâng cao toàn bộ thân đê, giải pháp xây dựng tường đỉnh chắn sóng đã được ứng dụng rộng rãi tại các địa phương như Hà Tĩnh (chiếm 50% chiều dài đê biển), Nam Định (chiếm 55%) và Hải Phòng. Mặc dù vậy, sự xuất hiện của tường đỉnh làm thay đổi căn bản chế độ thủy động lực học của dòng tràn, đồng thời làm phát sinh áp lực sóng cực lớn tác động trực tiếp lên kết cấu tường.

Luận văn thạc sĩ chuyên ngành Kỹ thuật Công trình tập trung giải quyết bài toán: Xác định áp lực sóng lớn nhất tác dụng lên tường đỉnh đê biển có mặt đứng và tường có mũi hắt sóng chắn tràn. Nghiên cứu giới hạn phạm vi địa lý tại dải ven biển Bắc Bộ và Bắc Trung Bộ, thực hiện tính toán ứng dụng trực tiếp cho tuyến đê biển Y Vích thuộc huyện Hậu Lộc, tỉnh Thanh Hóa với chiều dài 5,5 km, nơi trực tiếp bảo vệ an toàn cho 2.165 ha đất tự nhiên và hơn 60.000 người dân. Kết quả nghiên cứu giúp giảm thiểu nguy cơ mất ổn định trượt lật kết cấu tới 30-50%, tạo tiền đề kỹ thuật vững chắc để nâng cao năng lực ứng phó thiên tai biển trong bối cảnh biến đổi khí hậu và nước biển dâng toàn cầu.

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 thủy động lực học bờ biển và lý thuyết tương tác sóng - công trình ven bờ. Khung lý thuyết kế thừa các tiêu chuẩn kỹ thuật hiện đại gồm Sổ tay Kỹ thuật Bờ biển CEM-2006 của Quân đội Hoa Kỳ và Tiêu chuẩn Quốc gia TCVN 9901:2014 về yêu cầu thiết kế đê biển.

Cơ sở tính toán áp lực sóng dựa trên mô hình tải trọng sóng vỡ của Oumeraci (2001) cùng các nghiên cứu thực nghiệm máng sóng của Ramachandran (2012), Xuexue Chen (2014) và Koen Van Doorslaer (2015). Các khái niệm khoa học then chốt được vận hành trong nghiên cứu bao gồm:

  • Chiều cao sóng phổ mô men Hm0 và chu kỳ đỉnh phổ Tp biểu thị đặc trưng năng lượng sóng tới.
  • Độ cao lưu không Rc tính từ mực nước tính toán đến đỉnh tường.
  • Chỉ số sóng vỡ Iribarren xi_m xác định dạng sóng vỡ trên mái dốc đê biển.
  • Đại lượng áp lực sóng đặc trưng p1/250 (giá trị áp lực sóng trung bình của 1/250 số đợt tác động sóng lớn nhất).
  • Hệ số gia tăng áp lực sóng b_beta khi tường đỉnh tích hợp thêm mũi hắt sóng chắn tràn (recurve parapet).

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

Nghiên cứu sử dụng phương pháp mô hình hóa vật lý kết hợp phân tích thứ nguyên và hồi quy thống kê đa biến.

Toàn bộ dữ liệu thực nghiệm được thu thập từ hệ thống máng sóng vật lý quy mô lớn với tổng cộng 324 kịch bản thí nghiệm độc lập. Cỡ mẫu 324 kịch bản được thiết kế theo phương pháp tổ hợp trực giao đầy đủ nhằm bao quát toàn diện các yếu tố thủy lực và hình học: chiều cao sóng Hm0 thay đổi từ 0,058 m đến 0,14 m; chu kỳ sóng Tp biến thiên từ 1,15 s đến 2,0 s; độ sâu nước trong máng D từ 0,4 m đến 0,6 m; chiều cao tường đỉnh W gồm các mức 6 cm, 9 cm, 12 cm; bề rộng thềm trước S; độ dốc mái đê m = 3 và m = 4; góc mũi hắt sóng beta từ 0 độ (tường đứng phẳng) đến 45 độ và 90 độ.

Dữ liệu áp lực thủy động được đo đạc liên tục bằng cụm 3 cảm biến áp điện chuyên dụng PDB-200kPa (đường kính 6,5 mm, độ dày 1,0 mm do Tokyo Sokki Kenkyujo Nhật Bản sản xuất), kết nối qua bộ chuyển đổi tín hiệu đa kênh SDA-380C truyền trực tiếp về máy tính. Thiết bị vận hành với tần số quét cực cao đạt 100 Hz (chu kỳ lấy mẫu 10 ms), đảm bảo ghi nhận chính xác các đỉnh áp lực tức thời mà không bị trễ pha.

Lý do lựa chọn định lý Buckingham Pi (phương pháp phân tích thứ nguyên) là nhằm chuẩn hóa mối quan hệ phức tạp giữa nhiều đại lượng vật lý thành các nhóm biến số phi thứ nguyên độc lập. Tiếp đó, phần mềm tính toán khoa học Matlab được sử dụng để phân tích hồi quy phi tuyến nhiều biến, giúp thiết lập các phương trình giải tích tường minh xác định áp lực sóng và hệ số phân bố áp lực Cp với hệ số tương quan đạt độ 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

Quá trình phân tích 324 kịch bản thí nghiệm trên mô hình vật lý đã chỉ ra những phát hiện khoa học mang tính quy luật rõ rệt:

Thứ nhất, về dạng xung lực sóng, tín hiệu áp lực đo được trên tường đỉnh chủ yếu xuất hiện dưới dạng sóng vỡ thoát khí có 1 đỉnh (dạng Wagner) hoặc 2 đỉnh kế tiếp nhau. Hoàn toàn không phát hiện dạng áp lực xung kích nén khí túi lớn kiểu Bagnold do năng lượng sóng đã bị tiêu hao một phần khi truyền qua mái dốc đê phía biển (m = 3 đến 4). Chu kỳ dao động của các đỉnh áp lực cực đại dao động trong khoảng từ 0,1 giây đến 1,5 giây.

Thứ hai, mức độ tác động của các tham số thủy lực và hình học cho thấy áp lực sóng p1/250 tỷ lệ thuận với chiều cao sóng Hm0 và chu kỳ đỉnh phổ Tp (áp lực tăng bình quân từ 1,8 đến 2,5 lần khi chiều cao sóng tăng gấp đôi). Ngược lại, áp lực sóng tỷ lệ nghịch theo hàm mũ với độ cao lưu không tương đối Rc/Hm0 và bề rộng thềm trước tương đối S/Hm0. Việc bố trí thềm trước giúp giảm áp lực sóng tác động trực tiếp lên tường đỉnh từ 25% đến 40%.

Thứ ba, sự hiện diện của mũi hắt sóng (beta > 0 độ) làm gia tăng đáng kể áp lực sóng lên thân tường so với trường hợp tường thẳng đứng phẳng (beta = 0 độ). Áp lực tăng dần theo chiều cao mũi hắt tương đối hn/W và góc nghiêng beta (áp lực tại góc 90 độ tăng cao hơn từ 20% đến 35% so với góc 0 độ).

Thứ tư, áp lực sóng phân bố không đồng đều theo chiều cao thân tường. Giá trị áp lực đạt cực đại tại vị trí sát chân tường (độ cao Y = 1 cm tương ứng hệ số kinh nghiệm C1 = 3,82) và giảm dần đều khi lên cao (tại vị trí Y = 8,5 cm, hệ số phân bố C6 chỉ còn 1,12, tương đương mức suy giảm hơn 70% so với chân tường).

Thảo luận kết quả

Cơ chế vật lý của hiện tượng gia tăng áp lực khi có mũi hắt sóng bắt nguồn từ sự đổi hướng đột ngột của dòng tia nước tràn. Khi khối nước leo lên thân tường chạm vào mép mũi hắt có góc cong hoặc góc nghiêng, dòng tia bị chặn lại và phản hồi ngược ra phía biển, tạo ra một phản lực xung lượng bổ sung tác dụng ngược vào thành tường. Điều này giải thích vì sao hệ số b_beta luôn lớn hơn 1,0 khi góc beta lớn hơn 0 độ.

So sánh với các nghiên cứu tiền nhiệm, kết quả của luận văn hoàn toàn tương thích với quy luật suy giảm áp lực theo hàm mũ của Koen Van Doorslaer (2015), nhưng đã mở rộng và hoàn thiện hơn khi tích hợp đầy đủ ảnh hưởng của độ dốc mái đê m, bề rộng thềm trước S và hình học mũi hắt sóng, khắc phục triệt để hạn chế chỉ khảo sát một mặt cắt đơn lẻ trong thí nghiệm của Ramachandran (2012).

Các kết quả này được thể hiện trực quan thông qua các đường hồi quy phi tuyến giữa nhóm thông số không thứ nguyên p/(rhogHm0) với chỉ số Iribarren xi_m, kết hợp cùng biểu đồ bao phân bố áp lực và biểu đồ phân bố mô men uốn chân tường. Đây là tài liệu trực quan giúp các kỹ sư kết cấu dễ dàng tra cứu, tính toán tải trọng tác dụng và kiểm tra khả năng chịu lực của vật liệu bê tông cốt thép.

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

Dựa trên kết quả nghiên cứu thực nghiệm và tính toán kiểm chứng, luận văn đưa ra 4 nhóm đề xuất hành động cụ thể:

  1. Chuẩn hóa quy trình tính toán thiết kế kết cấu tường đỉnh đê biển: Bộ Nông nghiệp và Phát triển Nông thôn cùng các Viện Nghiên cứu Chuyên ngành cần tích hợp công thức xác định áp lực p1/250 và hệ số phân bố áp lực Cp vào sổ tay hướng dẫn tính toán thiết kế đê biển Việt Nam trong giai đoạn 2021-2023, giúp giảm thiểu sai số tính toán tải trọng động lực học từ mức 35% xuống dưới 10%.

  2. Áp dụng kết cấu tường đỉnh có thềm giảm sóng và mũi hắt sóng tối ưu: Các Ban Quản lý Dự án Công trình Nông nghiệp và Giao thông ven biển tại Thanh Hóa, Nam Định, Hà Tĩnh nên ưu tiên lựa chọn giải pháp tường đỉnh có thềm trước rộng từ 0,8 m đến 1,5 m kết hợp mũi hắt sóng góc 45 độ đến 60 độ. Giải pháp này giúp giảm thiểu 40-60% lượng sóng tràn trong khi vẫn kiểm soát được mức gia tăng mô men uốn chân tường dưới 25%.

  3. Cải tạo và nâng cấp đồng bộ tuyến đê biển xung yếu Y Vích (Hậu Lộc, Thanh Hóa): Sở Nông nghiệp và Phát triển Nông thôn tỉnh Thanh Hóa cần sớm triển khai dự án kiên cố hóa 5,5 km đoạn đê Y Vích với cao trình đỉnh tường +5,5 m, chiều cao tường W = 1,0 m theo biểu đồ tải trọng sóng mới, hoàn thành trong vòng 12 đến 18 tháng nhằm đảm bảo an toàn tuyệt đối cho 60.000 dân cư thuộc 5 xã ven biển.

  4. Thiết lập mạng lưới trạm quan trắc áp lực sóng tự động ven bờ: Tổng cục Phòng chống Thiên tai phối hợp với các cơ sở đào tạo đại học lắp đặt thử nghiệm hệ thống cảm biến đo áp lực sóng tần số cao (tối thiểu 100 Hz) trên 3 tuyến đê biển trọng điểm tại Hải Phòng, Nam Định và Hà Tĩnh trước năm 2025 nhằm thu thập dữ liệu sóng thực tế phục vụ kiểm chứng mô hình.

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

Luận văn là tài liệu tham khảo chuyên môn giá trị cho 4 nhóm đối tượng chính:

  1. Kỹ sư tư vấn thiết kế công trình thủy lợi, thủy công ven biển: Ứng dụng trực tiếp hệ thống công thức thực nghiệm và hệ số phân bố Cp để tính toán ổn định trượt, lật, kiểm tra nứt và bố trí cốt thép chịu uốn cho kết cấu tường đỉnh chắn sóng trên đê.

  2. Cán bộ quản lý đê điều và cơ quan quản lý nhà nước: Các Sở Nông nghiệp và Phát triển Nông thôn, Chi cục Thủy lợi và Ban Chỉ huy Phòng chống Thiên tai cấp tỉnh sử dụng luận văn làm cơ sở khoa học để thẩm định hồ sơ thiết kế, đánh giá mức độ an toàn của hệ thống đê hiện hữu trước nguy cơ bão cấp 12 đến cấp 14.

  3. Giảng viên, nhà nghiên cứu và học viên cao học chuyên ngành Công trình Thủy: Khai thác phương pháp luận thiết kế thí nghiệm máng sóng vật lý, kỹ thuật phân tích thứ nguyên Buckingham Pi và phương pháp xử lý số liệu áp lực thủy động lực học tần số cao trong các nghiên cứu chuyên sâu tiếp theo.

  4. Doanh nghiệp xây dựng và nhà thầu thi công kè biển: Nắm vững đặc tính phân bố tải trọng động lực học của sóng để lập biện pháp thi công tối ưu, lựa chọn mác bê tông phù hợp và bố trí mạch ngừng thi công hợp lý nhằm ngăn ngừa hiện tượng nứt gãy kết cấu.

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

  1. Vì sao áp lực sóng p1/250 được lựa chọn làm đại lượng tính toán đại diện trong nghiên cứu? Áp lực p1/250 đại diện cho giá trị trung bình của 1/250 số đợt sóng tác động mạnh nhất. Thống kê này loại bỏ các giá trị dị biệt ngẫu nhiên, vừa phản ánh đúng bản chất tải trọng cực trị nguy hiểm nhất trong suốt thời gian bão duy trì đỉnh triều, vừa đảm bảo tính kinh tế và an toàn kỹ thuật theo khuyến nghị quốc tế của Van Doorslaer (2015).

  2. Việc tích hợp mũi hắt sóng mang lại ưu điểm và thách thức gì cho kết cấu tường đỉnh? Mũi hắt sóng chuyển hướng dòng nước tràn ngược ra biển, giúp giảm đáng kể lượng nước tràn qua mặt đê. Tuy nhiên, lực phản hồi của dòng nước làm gia tăng áp lực sóng lên thân tường từ 20% đến 35%, đòi hỏi kết cấu thân tường và móng đê phải được gia cường khả năng chịu uốn và chống lật lớn hơn.

  3. Phương pháp phân tích thứ nguyên Buckingham Pi đóng vai trò gì trong việc xử lý số liệu? Định lý Buckingham Pi giúp quy đổi các biến số thứ nguyên rời rạc như chiều cao sóng, chu kỳ, chiều cao tường, bề rộng thềm thành các nhóm đại lượng phi thứ nguyên có ý nghĩa vật lý. Điều này cho phép khái quát hóa kết quả từ mô hình máng sóng thí nghiệm thành công thức tổng quát ứng dụng cho mọi công trình đê biển thực tế.

  4. Vì sao thí nghiệm không ghi nhận áp lực sóng xung kích dạng nén khí Bagnold? Áp lực sóng xung kích Bagnold chỉ hình thành khi sóng đập trực diện vào tường đứng ở vùng nước sâu. Tại đê biển có mái dốc thoải (m = 3 đến 4), sóng biển đã bị vỡ và tiêu tán năng lượng trên mái dốc trước khi chạm vào tường, khiến khối nước va đập ở dạng sóng vỡ thoát khí có đỉnh phẳng và êm hơn.

  5. Tuyến đê biển Y Vích tại Thanh Hóa đạt được lợi ích gì khi áp dụng công thức của đề tài? Khi áp dụng công thức mới cho 5,5 km đê Y Vích, kỹ sư xác định chính xác biểu đồ áp lực sóng và biểu đồ mô men uốn dọc thân tường với chiều cao 1,0 m, loại bỏ tình trạng thiết kế thiếu an toàn hoặc lãng phí vật liệu, bảo vệ vững chắc cho hơn 60.000 cư dân trước thiên tai.

Kết luận

Luận văn đã hoàn thành xuất sắc các mục tiêu nghiên cứu với những đóng góp học thuật và thực tiễn nổi bật:

  • Thiết lập thành công cơ sở dữ liệu thực nghiệm toàn diện từ 324 kịch bản máng sóng vật lý sử dụng cảm biến thủy động tần số cao 100 Hz.
  • Làm sáng tỏ quy luật tương tác thủy động và chứng minh áp lực sóng trên tường đỉnh đê biển dốc chủ yếu tồn tại ở dạng sóng vỡ thoát khí một hoặc hai đỉnh.
  • Xây dựng bộ công thức thực nghiệm hoàn chỉnh xác định áp lực sóng phi thứ nguyên có xét đến đầy đủ các yếu tố hình học và thủy lực.
  • Định lượng hóa chính xác hệ số gia tăng áp lực do mũi hắt sóng b_beta và quy luật phân bố áp lực Cp theo chiều cao thân tường.
  • Ứng dụng thành công vào bài toán thiết kế thực tế cho tuyến đê biển Y Vích thuộc huyện Hậu Lộc, tỉnh Thanh Hóa.

Trong lộ trình 2 đến 3 năm tới, hướng nghiên cứu cần tiếp tục mở rộng mô phỏng với tác động của sóng ngẫu nhiên 3 chiều và sóng xiên góc. Kết quả nghiên cứu là tài liệu tham khảo không thể thiếu cho các kỹ sư và nhà quản lý đê điều trong công cuộc bảo vệ bờ biển Việt Nam.