Tổng quan nghiên cứu

Việt Nam sở hữu dải bờ biển dài hơn 3.260 km trải dài từ Quảng Ninh đến Kiên Giang với 114 cửa sông lớn nhỏ. Khu vực ven biển bao gồm 28 trên tổng số 63 tỉnh, thành phố, đóng góp tới 42% diện tích tự nhiên và là nơi sinh sống của khoảng 15,5 triệu cư dân ven bờ cùng hơn 175 nghìn người dân tại hệ thống 3.000 đảo. Mật độ bờ biển trung bình đạt 1 km trên 100 km2 đất liền, cao gấp 6 lần mức trung bình toàn cầu. Tuy nhiên, dải ven duyên hải luôn đối mặt với các hiểm họa thiên tai cực đoan như nước dâng do bão, sóng thần và hoàn lưu dòng triều phức tạp, gây thiệt hại nghiêm trọng về kinh tế và con người.

Các mô hình thủy động lực học truyền thống thường dựa trên hệ phương trình nước nông với giả thiết phân bố áp suất thủy tĩnh. Giả thiết này đã triệt tiêu gia tốc thẳng đứng của chất lỏng, dẫn đến sai lệch lớn khi mô phỏng hiện tượng sóng biến dạng, sóng phân tán phi tuyến và sóng vỡ tại vùng nước ven bờ có đáy dốc thay đổi đột ngột.

Nghiên cứu tập trung giải quyết bài toán cốt lõi: Thiết lập mô hình số giải hệ phương trình nước nông phi thủy tĩnh (Non-hydrostatic Shallow Water Equations) nhằm tái hiện chính xác trường vận tốc ba chiều trung bình độ sâu, áp suất động và dao động mặt thoáng. Mục tiêu cụ thể là xây dựng thuật toán sai phân hữu hạn tối ưu, kiểm chứng qua mô hình máng sóng vật lý trong thời gian 40 giây và ứng dụng mô phỏng trường dòng triều thực tế tại Vịnh Đà Nẵng suốt chu kỳ 14 ngày. Kết quả nghiên cứu giúp nâng cao chỉ số hiệu quả Nash đạt từ 0,85 đến 0,99, giảm sai số biên độ sóng xuống dưới 2%, tạo công cụ tính toán tin cậy cho công tác cảnh báo thiên tai và bảo vệ công trình ven 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 thủy động lực học ba chiều xuất phát từ hệ phương trình Navier - Stokes đối với chất lỏng nhớt không nén được và phương trình liên tục. Để phản ánh chuyển động sóng trong vùng nước nông ven bờ, toàn bộ hệ phương trình được tích phân lấy trung bình theo phương thẳng đứng từ cao độ đáy z = -h đến cao độ mặt thoáng tự do z = ζ, với tổng độ sâu dòng chảy ký hiệu là D = ζ + h.

Khung lý thuyết phân tách trường áp suất toàn phần P thành hai thành phần riêng biệt: Áp suất thủy tĩnh tuyến tính theo độ sâu ρg(ζ - z) và áp suất động phi thủy tĩnh q(x, y) phát sinh do sự biến thiên gia tốc dòng chảy. Phương trình chuyển động theo phương ngang tích hợp hệ số ma sát đáy Chezy thông qua công thức Manning với hệ số nhám n = 0,022. Hệ thống phương trình khép kín bằng phương trình động lượng theo phương đứng mô tả vận tốc thẳng đứng W tại mặt thoáng và đáy biển, kết hợp phương trình Poisson đối với áp suất động phi thủy tĩnh. Mô hình còn tích hợp điều kiện biên dao động khô - ướt tự thích ứng dựa trên kiểm tra trạng thái tế bào lưới để theo dõi chính xác đường bờ di động khi sóng dâng hoặc triều rút.

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

Nguồn dữ liệu của nghiên cứu bao gồm hai tập dữ liệu chuẩn: Dữ liệu thực nghiệm vật lý trong máng sóng dài 37,7 m, sâu 0,8 m của nhóm tác giả Beji và Battjes (1993) với bãi ngầm hình thang cao 0,3 m và 8 điểm đo cao độ mực nước từ G1 đến G8; Dữ liệu thực địa số hóa địa hình đáy Vịnh Đà Nẵng cùng chuỗi 8 hằng số điều hòa thủy triều trạm Sơn Trà với mực nước trung bình 96,1 cm và sóng chính K1 có biên độ 18,33 cm.

Phương pháp phân tích dựa trên mô hình số sai phân hữu hạn (Finite Difference Method). Quá trình giải được phân tách thành hai bước tính toán độc lập: Bước tính thành phần thủy tĩnh sử dụng sơ đồ ngược dòng Upwind bảo toàn động lượng để xác định vận tốc tạm thời; Bước tính thành phần phi thủy tĩnh giải phương trình Poisson ma trận ngũ đường chéo nhằm tìm nghiệm áp suất động q, từ đó hiệu chỉnh vận tốc nằm ngang và vận tốc đứng.

Nghiên cứu khảo sát cỡ mẫu gồm 8 vị trí cảm biến đo đạc sóng trong phòng thí nghiệm và mạng lưới không gian với 6 mức tỷ lệ bước lưới từ L/10 đến L/200 (với bước sóng L = 4 m, chu kỳ sóng T = 2,02 s) cùng 5 cấp hệ số Courant từ 0,1 đến 0,5. Chuỗi mô phỏng thủy triều ven bờ được tính toán liên tục trong 14 ngày (từ ngày 01/05/2014 đến 14/05/2014). Phương pháp phân tách hai bước được lựa chọn vì đảm bảo tính ổn định tuyệt đối của nghiệm số, tiết kiệm tài nguyên bộ nhớ máy tính và mô tả mượt mà các dòng chảy không liên tục như sóng thứ cấp và sóng lan truyền qua bãi cạn.

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

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

Thứ nhất, mô hình phi thủy tĩnh đã tái hiện xuất sắc hiện tượng phân tán sóng phi tuyến và sự xuất hiện của các sóng thứ cấp tần số cao phía sau bãi ngầm hình thang cao 0,3 m. Tại các trạm đo từ G3 đến G8, đường cong mực nước tính toán bám sát tín hiệu thực nghiệm vật lý, trong khi mô hình thủy tĩnh truyền thống hoàn toàn làm phẳng và triệt tiêu các sóng thứ cấp này.

Thứ hai, đánh giá định lượng cho thấy mô hình đạt độ chuẩn xác vượt trội với hệ số hiệu quả Nash dao động từ 0,85 đến 0,99 tại hầu hết các trạm đo (trạm G1 đạt Nash 0,99; trạm G3 đạt 0,94). Sai số cực đại của mô hình phi thủy tĩnh giảm mạnh hơn 60% so với mô hình thủy tĩnh, khẳng định vai trò cốt yếu của thành phần áp suất động.

Thứ ba, phân tích độ nhạy bước lưới chỉ ra rằng khi kích thước lưới đạt mức L/50 so với bước sóng, sai số tương đối giảm xuống dưới 2%. Ngay cả khi sử dụng bước lưới thô ở mức L/10, sai số mô hình vẫn duy trì dưới ngưỡng 10%, chứng minh tính ổn định cao của sơ đồ sai phân đề xuất.

Thứ tư, khảo sát hệ số Courant trong dải từ Cr = 0,1 đến Cr = 0,5 cho thấy mức độ sai lệch kết quả giữa các trường hợp không vượt quá 0,05 cm. Việc thiết lập Cr = 0,5 cho phép tối đa hóa bước thời gian tính toán mà vẫn đảm bảo độ hội tụ, giúp giảm khoảng 40% tổng thời gian chạy mô phỏng.

Thứ năm, mô phỏng thủy triều tại Vịnh Đà Nẵng đã tái hiện chi tiết bức tranh thủy động lực học qua 4 pha triều chính (nước lớn, nước ròng, nước trung bình lên và xuống). Trường vận tốc dòng chảy được phân bố rõ nét tại các khu vực nhạy cảm quanh 77 điểm neo đậu tàu thuyền an toàn, mũi bán đảo Sơn Trà và hai cửa sông Hàn, sông Cu Đê.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự vượt trội của mô hình là việc giữ lại thành phần đạo hàm áp suất động theo phương thẳng đứng và phương ngang trong phương trình chuyển động. Khi sóng truyền từ vùng nước sâu 40 cm lên sườn dốc 1:20 của bãi ngầm, hiệu ứng phi tuyến tăng mạnh tạo ra các sóng điều hòa bậc cao. Mô hình thủy tĩnh do cố định áp suất theo trọng trường nên không thể sinh lực phân tán, dẫn đến biến dạng sóng sai lệch.

Kết quả này hoàn toàn tương đồng với các công bố học thuật quốc tế của Stelling, Zijlema và Yamazaki khi nghiên cứu về động lực học sóng ven bờ. Dữ liệu nghiên cứu được biểu diễn trực quan qua các bảng thống kê sai số bình phương trung bình, biểu đồ chuỗi thời gian so sánh dao động mực nước thực đo và tính toán tại 8 trạm kiểm chứng, cùng hệ thống bản đồ vector trường dòng chảy 2D tại các mốc thời gian 0h, 3h, 10h, 17h qua từng ngày triều. Sự phù hợp giữa phân tích hằng số điều hòa và kết quả mô phỏng tại trạm Sơn Trà đã khẳng định tính khả thi cao của mô hình khi áp dụng vào các thủy vực tự nhiên có địa hình đường bờ khúc khuỷu phức tạp.

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

Thứ nhất, chuẩn hóa tỷ lệ phân chia bước lưới tính toán không gian ở ngưỡng tối thiểu L/50 đối với các bài toán mô phỏng sóng ven bờ có địa hình biến đổi mạnh. Giải pháp này giúp kiểm soát sai số tương đối dưới 2% và cần được các kỹ sư công trình thủy áp dụng ngay trong giai đoạn thiết kế tính toán sóng tác động.

Thứ hai, thiết lập hệ số Courant cố định ở mức giới hạn trên Cr = 0,5 trong các phần mềm tính toán thủy động lực học vùng ven biển. Việc này giúp các cơ quan tư vấn tối ưu hóa thời gian chạy mô hình trên máy tính từ 30% đến 40% mà không gây mất ổn định nghiệm số, triển khai áp dụng ngay trong quý tiếp theo.

Thứ ba, mở rộng tích hợp module phi thủy tĩnh vào hệ thống dự báo thiên tai quốc gia nhằm mô phỏng bài toán tương tác giữa nước dâng do bão và sóng triều. Tổng cục Khí tượng Thủy văn cùng các viện nghiên cứu chuyên ngành nên chủ trì triển khai giải pháp này trong lộ trình 2025 - 2027 nhằm nâng cao năng lực cảnh báo ngập lụt sớm trước 24 đến 48 giờ cho 28 tỉnh ven biển.

Thứ tư, hoàn thiện điều kiện biên hấp thụ sóng phản xạ tại các biên mở nhân tạo. Các nhóm phát triển phần mềm mô phỏng cần nâng cấp thuật toán điều kiện biên phát xạ tự do trong vòng 6 tháng tới để khử triệt để sóng dội ngược, nâng cao độ chính xác tại các trạm đo nằm gần biên thoát của miền tính.

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

Nhóm học viên cao học, nghiên cứu sinh chuyên ngành Cơ học chất lỏng, Thủy văn và Toán tin ứng dụng có thể sử dụng tài liệu này như một cẩm nang chi tiết về kỹ thuật sai phân hữu hạn giải hệ phương trình đạo hàm riêng phi tuyến và phương pháp tách bước áp suất Poisson.

Nhóm kỹ sư thiết kế công trình thủy lợi, cảng biển và công trình ngoài khơi có thể khai thác các công thức tính toán tải trọng sóng, trường vận tốc và mực nước dâng để thiết kế kết cấu đê chắn sóng, bến cập tàu và các khu neo đậu an toàn ven bờ.

Nhóm chuyên gia tại các trung tâm dự báo khí tượng thủy văn và cơ quan phòng chống lụt bão có thể ứng dụng mô hình số trị để thiết lập các kịch bản dự báo ngập lụt ven duyên hải, mô phỏng lan truyền sóng thần và nước dâng trong bão với độ chính xác cao.

Nhóm cơ quan quản lý nhà nước về tài nguyên, môi trường biển đảo và quy hoạch đô thị ven biển tại các địa phương như Đà Nẵng, Quảng Nam, Hải Phòng có thể tham khảo trường hoàn lưu dòng chảy để đánh giá khả năng lan truyền chất ô nhiễm từ cửa sông và xây dựng chiến lược phát triển kinh tế biển bền vững.

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

Mô hình nước nông phi thủy tĩnh khác biệt như thế nào so với mô hình thủy tĩnh truyền thống? Mô hình phi thủy tĩnh tính thêm thành phần áp suất động q phát sinh từ gia tốc chuyển động thẳng đứng của chất lỏng. Điều này giúp mô hình mô phỏng chính xác hiện tượng sóng phân tán và sóng thứ cấp qua bãi ngầm, trong khi mô hình thủy tĩnh bỏ qua thành phần này khiến dạng sóng bị biến dạng phẳng và sai lệch biên độ.

Tại sao nghiên cứu phải chia quy trình giải số thành hai bước tách biệt? Phương pháp tách bước giúp cô lập thành phần thủy tĩnh phi tuyến giải bằng sơ đồ ngược dòng Upwind, sau đó giải phương trình Poisson để tìm áp suất động phi thủy tĩnh. Kỹ thuật này giúp chuyển bài toán động lượng phức tạp về hệ phương trình đại số tuyến tính dạng ma trận ngũ đường chéo, giúp tiết kiệm bộ nhớ và tăng tốc độ xử lý.

Độ dài bước lưới ảnh hưởng ra sao đến độ chính xác khi tính toán sóng nước nông? Độ dài bước lưới tỷ lệ nghịch với sai số tính toán. Khi bước lưới nhỏ hơn hoặc bằng L/50 so với chiều dài sóng, sai số tương đối đạt mức lý tưởng dưới 2%. Nếu bước lưới lớn hơn L/20 (như L/10 hay L/15), các dao động sóng thứ cấp tần số cao sẽ bị triệt tiêu, gây sai lệch lớn so với thực tế.

Thuật toán điều kiện biên khô - ướt xử lý biến động đường bờ theo nguyên lý nào? Thuật toán gán chỉ số trạng thái cho từng ô lưới dựa trên tổng độ sâu dòng chảy D. Ô có D lớn hơn 0 là ướt, ngược lại là khô. Khi nước dâng qua biên, ô khô được kiểm tra điều kiện cao độ mặt thoáng của ô liền kề để chuyển thành ô ướt và gán vận tốc tương ứng, đảm bảo tính liên tục của dòng chảy tràn bờ.

Mô hình đã chứng minh tính hiệu quả thực tế tại Vịnh Đà Nẵng như thế nào? Mô hình đã mô phỏng liên tục 14 ngày triều tại Vịnh Đà Nẵng dựa trên 8 hằng số điều hòa trạm Sơn Trà. Kết quả đã tái hiện chi tiết phân bố vận tốc dòng chảy qua 4 pha triều tại các cửa sông Hàn, sông Cu Đê và 77 điểm neo đậu tàu thuyền, chứng minh khả năng ứng dụng cho vùng bờ biển có địa hình phức tạp.

Kết luận

Nghiên cứu đã giải quyết trọn vẹn bài toán xây dựng và ứng dụng mô hình số hệ phương trình nước nông phi thủy tĩnh với 5 kết quả đóng góp cốt lõi:

  • Thiết lập thành công thuật toán số trị phân tách hai bước giải hệ phương trình nước nông phi thủy tĩnh kết hợp điều kiện biên khô - ướt linh hoạt.
  • Kiểm chứng hoàn hảo mô hình qua thí nghiệm vật lý của Beji và Battjes (1993), đạt chỉ số hiệu quả Nash xuất sắc từ 0,85 đến 0,99 và tái hiện chính xác sóng thứ cấp.
  • Xác lập các chuẩn thông số kỹ thuật tối ưu cho tính toán sóng gồm bước lưới L/50 và hệ số Courant Cr = 0,5 giúp cân bằng giữa độ chính xác cao và hiệu năng xử lý.
  • Ứng dụng mô phỏng thành công hoàn lưu thủy triều thực tế suốt 14 ngày tại khu vực Vịnh Đà Nẵng với hệ số nhám đáy Manning n = 0,022.
  • Mở ra hướng tiếp cận học thuật tiên tiến trong việc mô phỏng sóng phân tán ven bờ, thay thế hiệu quả cho các mô hình thủy tĩnh truyền thống vốn còn nhiều hạn chế.

Trong giai đoạn tiếp theo, mô hình cần được tiếp tục nâng cấp để tích hợp trường ứng suất gió bão và lan truyền sóng thần đáy biển. Các nhà khoa học, kỹ sư công trình và nhà quản lý tài nguyên biển hãy khai thác và phát triển khung mô hình này nhằm phục vụ hiệu quả công tác quy hoạch công trình cảng biển và chủ động ứng phó với biến đổi khí hậu ven bờ.