Tổng quan nghiên cứu

Việt Nam sở hữu đường bờ biển dài hơn 3.260 km cùng vùng đặc quyền kinh tế rộng lớn hơn 1 triệu km², tạo điều kiện thuận lợi để phát triển mạng lưới hơn 34 cảng biển huyết mạch. Trong bối cảnh biến đổi khí hậu và tần suất bão biển gia tăng, các công trình bảo vệ bờ như đê chắn sóng dạng tường đứng và đê chắn sóng dạng mái nghiêng phải thường xuyên chịu tác động của những đợt sóng bão có chiều cao sóng có nghĩa vượt trên 5,0 m đến 8,5 m. Việc xác định sai lệch tải trọng sóng trong giai đoạn thiết kế có thể làm tăng chi phí xây dựng từ 15% đến 30% do dư thừa vật liệu, hoặc nghiêm trọng hơn là gây phá hủy cục bộ thân đê khi xảy ra bão lớn, dẫn đến thiệt hại kinh tế hàng trăm tỷ đồng.

Mục tiêu cốt lõi của nghiên cứu là hệ thống hóa, khảo sát và đánh giá so sánh chuyên sâu 8 phương pháp tính toán tải trọng sóng phổ biến trên thế giới, bao gồm lý thuyết sóng tuyến tính, sóng phi tuyến bậc cao, giải pháp giải tích phương trình Laplace bậc bốn và mô hình số Navier-Stokes hai chiều. Trên cơ sở đó, tác giả kiểm chứng độ tin cậy của từng phương pháp với số liệu thực nghiệm chuẩn quốc tế để đề xuất quy trình tính toán tối ưu, phù hợp nhất với điều kiện thủy hải văn ven bờ Việt Nam.

Nghiên cứu tập trung vào hai cơ chế thủy động lực học điển hình: sóng không vỡ (sóng đứng) và sóng vỡ tác động lên đê tường đứng (khối xếp, thùng chìm caisson) cũng như đê mái nghiêng đá đổ. Kết quả của luận văn giúp nâng cao độ chính xác trong dự báo áp lực sóng lên thân đê thêm 12% đến 25% so với các công thức kinh nghiệm cổ điển, đồng thời giảm hệ số phân tán rủi ro trong tính toán ổn định trượt lật xuống dưới mức 5,0%.

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 thủy động lực học công trình biển phi tuyến, tích hợp các hệ phương trình cơ bản về bảo toàn khối lượng và động lượng trong môi trường chất lỏng lý tưởng và chất lỏng thực. Hai nhánh lý thuyết chính được áp dụng bao gồm:

Thứ nhất là lý thuyết thế vận tốc cho bài toán giá trị biên phi tuyến dựa trên phương trình Laplace đối với dòng chảy không xoáy. Hiện tượng sóng đứng toàn phần trước công trình phản xạ tuyệt đối được mô tả bằng phép khai triển chuỗi xấp xỉ liên tục bậc bốn của Tadjbaksh, Keller (1960) và Goda, Kakizaki (1966). Mô hình này xác định trường áp lực động lực học và độ dâng mực nước đỉnh sóng dựa trên các điều kiện biên động học và động lực học tại mặt thoáng tự do.

Thứ hai là hệ phương trình vi phân Navier-Stokes 2D hai chiều mô tả dòng chảy rối thực có xét đến độ nhớt động học của nước biển, kết hợp phương trình liên tục và phương trình biến thiên mực nước tự do. Tác giả vận dụng mô hình số cải tiến của Duy (1996) với phép biến đổi tọa độ Jacobian nhằm giải bài toán mặt thoáng di động và địa hình đáy biến đổi.

Các khái niệm chuyên ngành cốt lõi được chuẩn hóa gồm:

  • Áp lực sóng không vỡ (áp lực sóng đứng hình thành do sóng phản xạ với hệ số phản xạ dao động từ 0,90 đến 1,00).
  • Áp lực đẩy nổi đáy công trình (chiếm khoảng 30% đến 50% tổng áp lực tác động lên cấu kiện chìm).
  • Tham số phi thứ nguyên Ursell (phân định miền áp dụng lý thuyết sóng Stokes bậc 5 khi tham số nhỏ hơn 25 và lý thuyết sóng Cnoidal bậc 3 khi tham số lớn hơn hoặc bằng 25).

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

Nguồn dữ liệu kiểm chứng của đề tài được trích xuất từ chuỗi thí nghiệm mô hình vật lý chính xác của Goda và Kakizaki (1966). Cỡ mẫu nghiên cứu bao gồm 5 bộ kịch bản thủy động lực học đại diện (từ trường hợp 1 đến trường hợp 5), với dải độ sâu nước tương đối biến thiên từ 0,12 đến 0,48 và độ dốc sóng dao động trong khoảng 0,018 đến 0,062.

Phương pháp chọn mẫu dựa trên nguyên tắc bao phủ toàn diện các trạng thái sóng từ vùng nước nông, nước trung bình đến nước sâu, đồng thời bao hàm cả dải sóng có độ dốc thấp lẫn sóng dốc cao tiệm cận giới hạn vỡ.

Phương pháp phân tích số trị sử dụng thuật toán sai phân hữu hạn (FDM) trên hệ lưới so le không đều. Lý do lựa chọn mô hình lưới so le kết hợp phép chuyển đổi tọa độ thực sang miền tính toán chữ nhật là nhằm đảm bảo bậc chính xác cấp 2 đối với mọi đạo hàm riêng phân, triệt tiêu hiện tượng dao động số không vật lý và tối ưu hóa thời gian tính toán ma trận tuyến tính thông qua thuật toán khử Gauss và truy đuổi lặp. Toàn bộ chu trình tính toán, lập trình và đối sánh thực nghiệm được thực hiện liên tục trong khung thời gian 12 tháng.

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

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

Phân tích định lượng từ 5 trường hợp kiểm chứng thực nghiệm đã đưa ra 4 phát hiện quan trọng về mức độ chính xác của các phương pháp tính toán tải trọng sóng:

Thứ nhất, mô hình giải tích bậc bốn của phương trình Laplace theo Goda - Kakizaki (1966) và mô hình số Navier-Stokes của Duy (1996) cho độ chính xác cao nhất. Sai số tương đối trung bình giữa giá trị tính toán và số liệu đo đạc thực nghiệm chỉ dao động trong khoảng 3,2% đến 4,8% trên toàn bộ chiều cao mặt cắt ướt của tường đứng.

Thứ hai, phương pháp cổ điển Sainflou (1928) bộc lộ sự sai lệch rất lớn khi áp dụng cho các đợt sóng có độ dốc cao (với tỷ số chiều cao trên chiều dài sóng vượt mức 0,040). Cụ thể, giá trị áp lực sóng tính toán theo Sainflou cao hơn từ 22,5% đến 41,2% so với giá trị đo đạc thực nghiệm, dẫn đến tình trạng tính toán quá dư thừa kích thước kết cấu.

Thứ ba, Tiêu chuẩn ngành 22 TCN 222-95 của Việt Nam cho kết quả tính toán áp lực sóng đứng rất tiệm cận với số liệu thực nghiệm tại vùng nước nông và trung bình, với mức sai số kiểm định chỉ từ 5,8% đến 8,6%. Điều này khẳng định biểu đồ hệ số kinh nghiệm trong tiêu chuẩn Việt Nam có độ an toàn và tin cậy cao.

Thứ tư, đối với đê chắn sóng mái nghiêng, áp lực sóng rút theo phương pháp Djunkovski (1949) và Tiêu chuẩn 22 TCN 222-95 cho thấy áp lực thủy động âm gây lực hút ngược lên lớp đá phủ đạt giá trị cực đại bằng khoảng 0,45 đến 0,65 lần áp lực thủy tĩnh tại chân sóng, đòi hỏi khối lượng viên phủ bảo vệ phải đạt tiêu chuẩn ổn định trọng lượng tối thiểu từ 3,5 tấn đến 12,0 tấn tùy thuộc độ dốc mái đê từ 1:1,5 đến 1:2,5.

Thảo luận kết quả

Nguyên nhân dẫn đến sự sai khác vượt trội của phương pháp Sainflou bắt nguồn từ giả thiết gần đúng bậc hai không hoàn toàn dựa trên sóng trocoide, bỏ qua các tương tác phi tuyến bậc cao của các thành phần điều hòa khi sóng tiến gần bờ. Ngược lại, phương pháp Miche-Rundgren (1958) với việc hiệu chỉnh hệ số phản xạ từ 0,90 đến 1,00 đã cải thiện đáng kể kết quả, đưa sai số về mức dưới 10,5% đối với sóng dốc.

Trong khi đó, Tiêu chuẩn kỹ thuật công trình cảng Nhật Bản (2002) và Tiêu chuẩn 22 TCN 222-95 có ưu điểm là đơn giản hóa biểu đồ áp lực thành dạng hình thang và tuyến tính phân đoạn, vừa dễ ứng dụng trong thực hành kỹ thuật vừa đảm bảo hệ số dự phòng an toàn từ 1,15 đến 1,25 cho kết cấu.

Dữ liệu so sánh trong luận văn được trực quan hóa rất hiệu quả thông qua hệ thống bảng đối chiếu sai số tương đối và biểu đồ phân bố áp lực không thứ nguyên theo chiều sâu thẳng đứng. Trên biểu đồ mặt cắt từ cao trình đỉnh sóng leo (+1,25H) đến đáy biển (-d), các đường cong áp lực của mô hình Navier-Stokes và phương trình Laplace bậc bốn bám sát các điểm đo thực nghiệm, thể hiện rõ vùng tập trung áp lực cực đại tại vị trí mực nước tĩnh với giá trị đạt khoảng 1,5 đến 2,2 lần áp lực thủy tĩnh cơ bản.

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

Nhằm hoàn thiện quy trình thiết kế và nâng cao chất lượng xây dựng đê chắn sóng tại Việt Nam, luận văn đưa ra 4 nhóm giải pháp kỹ thuật cụ thể:

  1. Rà soát và cập nhật tiêu chuẩn thiết kế công trình thủy: Bộ Xây dựng và Bộ Giao thông Vận tải cần phối hợp bổ sung các thuật toán tính toán phi tuyến bậc cao vào Tiêu chuẩn Quốc gia TCVN trong lộ trình 24 tháng tới, giúp tối ưu hóa khối lượng bê tông đúc sẵn và đá hộc từ 8,0% đến 15,0% cho các dự án đê kè biển.
  2. Ứng dụng mô hình số hóa Navier-Stokes cho các dự án cảng nước sâu: Các viện nghiên cứu và doanh nghiệp tư vấn thiết kế công trình biển (trong giai đoạn 2026-2028) cần triển khai phần mềm thủy động lực học giải hệ phương trình Navier-Stokes 2D/3D đối với các tuyến đê chắn sóng caisson có độ sâu nước trên 15,0 m hoặc vùng có chiều cao sóng bão thiết kế vượt 6,0 m.
  3. Thiết lập mạng lưới quan trắc áp lực sóng thực địa: Ban quản lý các dự án cảng biển lớn như Lạch Huyện, Liên Chiểu, Cái Mép - Thị Vải cần lắp đặt hệ thống cảm biến áp lực sóng đa điểm tại 3 đến 5 vị trí xung yếu trên thân đê, theo dõi liên tục trong chu kỳ 36 tháng để thu thập số liệu bão thực tế nhằm hiệu chỉnh hệ số an toàn thiết kế nội địa.
  4. Kiểm soát ổn định mái dốc và lớp phủ đê chắn sóng: Đơn vị thiết kế cần kết hợp phương pháp Klein Breteler và 22 TCN 222-95 để xác định chính xác hệ số ổn định của khối phủ nhân tạo (như Tetrapod, Accropode), đảm bảo chỉ số ổn định đạt từ 7,0 đến 8,5, giúp giảm nguy cơ xói lở chân đê xuống dưới mức 2,0% sau 50 năm khai thác.

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

Nội dung và kết quả nghiên cứu của luận văn mang lại giá trị thực tiễn và học thuật cho 4 nhóm đối tượng chính:

  • Kỹ sư tư vấn thiết kế công trình thủy và thềm lục địa: Nắm bắt và làm chủ công thức tính toán áp lực sóng đứng, sóng vỡ và áp lực đẩy nổi, giúp giảm từ 10% đến 18% chi phí xây lắp nhờ chọn giải pháp kết cấu chính xác.
  • Chủ đầu tư và Ban quản lý dự án hạ tầng hàng hải: Có cơ sở khoa học độc lập để thẩm tra thiết kế cơ sở, kiểm soát chặt chẽ rủi ro mất ổn định công trình trước các kịch bản sóng bão chu kỳ 50 đến 100 năm.
  • Giảng viên, nghiên cứu sinh và học viên cao học: Sử dụng nguồn dữ liệu toán học về phương trình vi phân Laplace bậc bốn và mô hình sai phân Navier-Stokes làm tài liệu giảng dạy, phát triển các hướng nghiên cứu chuyên sâu về công trình biển.
  • Cơ quan quản lý nhà nước và hội đồng chuyên ngành: Sử dụng kết quả nghiên cứu làm cơ sở thực chứng để biên soạn, soát xét và ban hành các quy chuẩn kỹ thuật quốc gia về tải trọng tác động lên công trình thủy ven bờ.

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

Tại sao phương pháp Sainflou truyền thống không còn phù hợp với sóng có độ dốc cao?
Phương pháp Sainflou dựa trên lý thuyết sóng trocoide gần đúng bậc hai không hoàn toàn. Khi độ dốc sóng lớn hơn 0,040, công thức này bỏ qua sự biến dạng phi tuyến mạnh của sóng, làm áp lực tính toán cao hơn thực tế từ 22% đến 41%, gây lãng phí kết cấu nghiêm trọng.

Tiêu chuẩn 22 TCN 222-95 của Việt Nam có độ chính xác như thế nào so với thực nghiệm?
Kiểm định với số liệu của Goda và Kakizaki cho thấy tiêu chuẩn 22 TCN 222-95 đạt độ chính xác cao tại vùng nước nông và trung bình, với độ lệch chỉ từ 5,8% đến 8,6%, hoàn toàn đáp ứng yêu cầu an toàn kỹ thuật trong điều kiện ven bờ Việt Nam.

Ưu điểm vượt trội của mô hình số Navier-Stokes (Duy, 1996) là gì?
Mô hình Navier-Stokes giải quyết trực tiếp dòng chảy rối thực trên hệ lưới so le biến đổi tọa độ Jacobian. Phương pháp này mô phỏng chính xác sự biến đổi tức thời của biên mặt thoáng di động, đạt độ chính xác bậc hai với sai số dưới 4,5% so với thực nghiệm.

Khi nào cần tính toán tải trọng sóng vỡ thay vì sóng không vỡ?
Khi tỷ số giữa chiều sâu nước và chiều cao sóng trước công trình nhỏ hơn 1,28, sóng sẽ bị biến dạng và vỡ ngay tại chân đê. Lúc này, áp lực va đập xung kích có thể tăng đột biến gấp 2,0 đến 3,5 lần so với sóng đứng, đòi hỏi phải áp dụng công thức Hiroi hoặc Minikin.

Hệ số phản xạ của công trình ảnh hưởng ra sao đến tải trọng sóng thiết kế?
Hệ số phản xạ đại diện cho mức độ phản xạ năng lượng sóng trước tường đứng. Đối với tường trơn nhẵn thẳng đứng, hệ số này tiệm cận 1,00; khi bề mặt có độ nhám hoặc đệm đá, hệ số giảm xuống khoảng 0,90, giúp giảm áp lực sóng tác động lên thân đê khoảng 5% đến 10%.

Kết luận

  • Luận văn đã hệ thống hóa toàn diện 8 phương pháp tính toán tải trọng sóng tác động lên đê chắn sóng tường đứng và đê mái nghiêng phổ biến trên thế giới và tại Việt Nam.
  • Khẳng định mô hình phương trình vi phân Laplace bậc bốn và mô hình số Navier-Stokes 2D có độ chính xác cao nhất với sai số dưới 4,8% so với dữ liệu thực nghiệm.
  • Chứng minh tính phù hợp và độ tin cậy của Tiêu chuẩn ngành 22 TCN 222-95 trong việc xác định tải trọng sóng đối với các công trình thủy ven bờ Việt Nam.
  • Chỉ rõ giới hạn áp dụng và mức độ sai lệch của phương pháp Sainflou cổ điển đối với các đợt sóng có độ dốc cao tiệm cận giới hạn vỡ.
  • Đề xuất các giải pháp kỹ thuật mang tính ứng dụng cao giúp tối ưu hóa từ 8% đến 15% khối lượng vật liệu xây dựng kết cấu đê chắn sóng.

Đóng góp lớn nhất của công trình nghiên cứu là cung cấp bức tranh đối sánh định lượng khoa học và chuẩn xác giữa các phương pháp lý thuyết và thực nghiệm, định hình hướng dẫn tính toán tối ưu cho kỹ sư công trình biển. Lộ trình tiếp theo trong 12 đến 24 tháng tới là mở rộng kiểm thử mô hình số Navier-Stokes 3D tích hợp tương tác giữa sóng và nền đất yếu ven biển. Kính mời các chuyên gia, kỹ sư và nhà nghiên cứu quan tâm tham khảo toàn văn tài liệu để áp dụng hiệu quả vào các dự án cảng biển thực tế.