Tổng quan nghiên cứu

Cửa sông Châu Giang là một trong những vùng kinh tế hàng hải sôi động bậc nhất khu vực Đông Á với diện tích mặt nước rộng hơn 2.000 km2 và chiều dài luồng lạch khoảng 65 km từ Hổ Môn ra Biển Đông. Hàng năm, hệ thống sông này tiếp nhận tổng lượng dòng chảy lên tới 179,3 tỷ m3 nước ngọt và hơn 30,8 triệu tấn phù sa lơ lửng từ 4 cửa sông chính gồm Hổ Môn, Tiêu Môn, Hồng Kỳ Môn và Hoành Môn. Sự tương tác phức tạp giữa dòng chảy sông ngòi, dòng triều bán nhật triều không đều có biên độ trung bình 2,0 m (cực đại 3,4 m) và tác động của sóng gió đã tạo nên hiện tượng vùng cực đại độ đục. Vùng cực đại độ đục gây ra sự bồi lắng nghiêm trọng tại các luồng hàng hải trọng yếu, đặc biệt là Luồng Đông và Luồng Tây có độ sâu từ 5 m đến 20 m. Mục tiêu của nghiên cứu là xây dựng mô hình số trị ba chiều tích hợp tương tác sóng và dòng triều nhằm làm sáng tỏ cơ chế hình thành, phát triển và tiêu tán của vùng cực đại độ đục. Nghiên cứu thực hiện đánh giá toàn diện trên cơ sở dữ liệu quan trắc thủy văn đa thời kỳ từ năm 1978 đến năm 1998 tại vùng vịnh Linh Đinh, Hong Kong và Macau. Kết quả nghiên cứu có ý nghĩa quyết định trong việc tối ưu hóa chi phí nạo vét luồng tàu, giảm thiểu rủi ro sa bồi luồng hàng hải và nâng cao hiệu quả quản trị vùng ven biển với mục tiêu cắt giảm ít nhất 15% khối lượng nạo vét duy tu thường niê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 lý thuyết động lực học chất lỏng địa vật lý và lý thuyết vận chuyển bùn cát hạt mịn kết dính trong môi trường cửa sông chịu ảnh hưởng của triều. Ba khái niệm trung tâm bao gồm hoàn lưu trọng lực cửa sông do chênh lệch mật độ muối, hiện tượng bơm triều do tính bất đối xứng giữa triều lên và triều rút, cùng quá trình tái lơ lửng bùn cát đáy do ứng suất cắt đáy biển vượt ngưỡng ma sát tới hạn. Khung lý thuyết kết hợp phương trình Navier-Stokes ba chiều phi tĩnh học với hệ tọa độ chuẩn hóa sigma và mô hình khép kín rối bậc 2,5 của Mellor và Yamada để mô phỏng chính xác sự dập tắt nhiễu loạn do phân tầng mật độ nước. Về động lực học bùn cát, lý thuyết khuếch tán thẳng đứng biên dạng Rouse được ứng dụng đồng thời với mô hình bảo toàn tác động sóng nhằm tính toán ứng suất cắt tổng hợp đáy biển dưới tác động đồng thời của sóng và dòng chảy. Độ mặn được mô hình hóa biến thiên từ 0 ppt tại thượng lưu đến 33 ppt ngoài khơi, trong đó vùng keo tụ hạt mịn phát triển mạnh nhất tại dải độ mặn từ 3 ppt đến 10 ppt, giải thích sự gia tăng tốc độ lắng đọng bùn cát có đường kính trung vị khoảng 21,6 micromet.

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

Phương pháp nghiên cứu kết hợp giữa phân tích số liệu thực địa và mô phỏng số trị 2 chiều tích phân theo độ sâu cùng mô hình 3 chiều phi tuyến. Nguồn dữ liệu quan trắc đồng bộ được thu thập qua các chiến dịch khảo sát thủy văn quy mô lớn vào mùa lũ tháng 7 năm 1978 và mùa cạn tháng 3 năm 1979, bổ sung dữ liệu kiểm chứng toàn diện các năm 1992, 1993 và năm 1998. Cỡ mẫu nghiên cứu bao gồm chuỗi đo đạc liên tục 25 giờ tại 16 trạm đo cố định bố trí dọc theo các luồng hàng hải chính và khu vực bãi cạn. Phương pháp chọn mẫu phân tầng theo vị trí địa lý và điều kiện thủy triều được áp dụng để đảm bảo tính đại diện cho cả chu kỳ triều cường và triều kém qua các mùa. Lý do lựa chọn mô hình số trị 3 chiều kết hợp mô hình sóng là vì vùng cửa sông Châu Giang có địa hình đáy phức tạp với 2 lạch sâu và 3 bãi cạn rộng lớn xen kẽ, khiến các dòng chảy phân tầng thẳng đứng và dòng dư hoàn lưu ngang không thể biểu diễn chính xác bằng các mô hình giản lược. Quá trình tính toán số trị sử dụng phương pháp chia tách toán tử và sai phân hữu hạn trên lưới phi cấu trúc với 10 lớp phân tầng thẳng đứng, đảm bảo độ ổn định và tối ưu hóa thời gian tính toá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 ba chiều đã tái hiện chính xác cấu trúc và quy luật dịch chuyển của vùng cực đại độ đục, với biên độ dịch chuyển dọc luồng đạt khoảng 22 km dọc theo các cồn cát ngầm chính giữa triều lên và triều rút. Hiện tượng cực đại độ đục biểu hiện rõ rệt nhất trong các kỳ triều cường và gần như biến mất trong các kỳ triều kém.

Thứ hai, phân tích phân bố không gian cho thấy nồng độ bùn cát lơ lửng tại các bãi cạn luôn cao hơn từ 40% đến 60% so với trong các lạch sâu lân cận. Trong mùa lũ, 4 cửa sông phía Tây xả ra từ 92% đến 96% tổng lượng bùn cát cả năm, nhưng có tới gần 80% lượng phù sa này bị giữ lại và lắng đọng ngay trong phạm vi vịnh Linh Đinh.

Thứ ba, trường dòng chảy dư Eulerian đóng vai trò chi phối với vận tốc dòng chảy đạt cực đại khoảng 0,3 m/s ở khu vực cửa sông, cao gấp hơn 3 lần so với thành phần trôi dạt Stokes vốn chỉ đạt dưới 0,1 m/s.

Thứ tư, tương tác sóng và dòng chảy làm gia tăng ứng suất cắt đáy tại các bãi cạn ven bờ từ 25% đến 50%, làm tăng độ dày của lõi bùn cát nồng độ cao thẳng đứng nhưng không làm thay đổi đáng kể vị trí trung tâm của vùng cực đại độ đục.

Thảo luận kết quả

Nguyên nhân cốt lõi tạo nên sự phân hóa vùng cực đại độ đục giữa Luồng Tây và Luồng Đông xuất phát từ sự khác biệt về động lực học thủy văn. Tại Luồng Tây, hoàn lưu trọng lực kết hợp với hiệu ứng bẫy triều là động lực chính tạo nên sự tích tụ bùn cát do tiếp nhận trực tiếp nguồn xả nước ngọt khổng lồ từ 3 cửa Tiêu Môn, Hồng Kỳ Môn và Hoành Môn. Ngược lại, tại Luồng Đông, hiện tượng cực đại độ đục chủ yếu bị chi phối bởi các chu kỳ lắng đọng và tái lơ lửng bùn cát cục bộ do ứng suất cắt đáy mạnh mẽ khi triều rút. So sánh với các nghiên cứu tại cửa sông Dương Tử hay cửa sông Gironde, cơ chế bẫy phù sa tại cửa sông Châu Giang phức tạp hơn nhiều do sự hiện diện của hệ thống 3 bãi cạn và các đảo chắn sóng.

Dữ liệu động lực học này có thể được trình bày trực quan thông qua bảng tổng hợp 8 thành phần phân rã thông lượng bùn cát tịnh theo thời gian thực và biểu đồ mặt cắt thẳng đứng thể hiện mối tương quan giữa nồng độ bùn cát lơ lửng trung bình 0,172 g/l với gradien độ mặn từ 0 đến 33 ppt. Việc đối chiếu các biểu đồ đồng mức nồng độ trầm tích mô phỏng với ảnh vệ tinh cho thấy sự tương thích rất cao về mặt hình thái không gian, chứng minh tính tin cậy tuyệt đối của phương pháp mô hình hóa ba chiều.

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

Thứ nhất, tối ưu hóa quy hoạch tuyến luồng và lịch trình nạo vét duy tu luồng tàu hàng hải. Cục Hàng hải cùng các Ban Quản lý Cảng biển cần thực hiện nạo vét tập trung vào giai đoạn triều kém của mùa khô với mốc thời gian từ tháng 11 đến tháng 3 hàng năm, giúp giảm thiểu 20% hiện tượng tái bồi lắng nhanh và tiết kiệm từ 15% đến 25% ngân sách duy tu thường niên đối với các đoạn luồng sâu trên 12,5 m.

Thứ hai, thiết lập hệ thống trạm quan trắc số trị và cảnh báo bồi lắng tự động thời gian thực. Viện Nghiên cứu Khoa học Thủy lợi và Thủy điện cần phối hợp với Trung tâm Khí tượng Thủy văn triển khai mạng lưới 10 trạm đo cảm biến quang học và vận tốc dòng chảy trong khung thời gian 12 tháng, hướng tới mục tiêu nâng cao độ chính xác dự báo bồi tích đạt trên 85%.

Thứ ba, xây dựng công trình chỉnh trị giảm sóng và hướng dòng tại các cồn cát ngầm trọng điểm. Ban Quản lý Dự án Công trình Giao thông Thủy cần triển khai thi công đê ngầm giảm sóng tại khu vực bãi cạn phía Tây trong lộ trình 24 tháng, nhằm kiểm soát độ dày lõi độ đục và hạn chế 30% lượng bùn cát mịn tràn từ bãi cạn vào luồng tàu chính.

Thứ tư, hoàn thiện quy chế phối hợp liên ngành trong việc điều tiết xả lũ và bùn cát từ các hồ chứa thượng nguồn. Ủy ban Quản lý Lưu vực sông cần ban hành hướng dẫn điều tiết lưu lượng dòng chảy qua 4 cửa sông chính trong giai đoạn 6 tháng đầu năm, đảm bảo duy trì lưu lượng tối thiểu để hạn chế mặn xâm nhập sâu quá 25 km vào nội địa.

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

Thứ nhất, kỹ sư công trình thủy và chuyên gia nạo vét hàng hải: Tiếp cận phương pháp dự báo chính xác vị trí tích tụ bùn cát và quy luật dịch chuyển 22 km của vùng độ đục cực đại để thiết kế tuyến luồng và lập kế hoạch nạo vét luồng hàng hải đạt chuẩn độ sâu từ 12,5 m đến 15 m với chi phí thấp nhất.

Thứ hai, nhà quản lý quy hoạch và tài nguyên môi trường ven biển: Sử dụng bức tranh tổng thể về sự phân bố của 30,8 triệu tấn bùn cát hàng năm và cơ chế bẫy trầm tích để khoanh vùng bảo tồn bãi cạn sinh thái và quy hoạch không gian phát triển cảng biển bền vững trong tầm nhìn 10 đến 20 năm.

Thứ ba, nhà nghiên cứu động lực học thủy văn và mô hình hóa ven biển: Ứng dụng thuật toán chia tách toán tử, hệ tọa độ sigma và mô hình khép kín rối bậc 2,5 để phát triển các mô hình 3 chiều mô phỏng lan truyền sóng, dòng triều và độ mặn cho các vùng cửa sông ven biển phức tạp.

Thứ tư, giảng viên và học viên cao học chuyên ngành kỹ thuật công trình biển: Khai thác bộ tư liệu chuẩn mực với hơn 100 biểu đồ phân bố thủy động lực học và dữ liệu quan trắc thực địa đa thời kỳ phục vụ giảng dạy, nghiên cứu chuyên sâu về cơ chế bồi tụ bờ biển.

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

  1. Hiện tượng cực đại độ đục tại cửa sông Châu Giang có đặc điểm dịch chuyển như thế nào?
    Vùng cực đại độ đục dịch chuyển tuần hoàn với biên độ khoảng 22 km dọc theo các luồng hàng hải chính giữa kỳ triều lên và triều rút. Hiện tượng này phát triển mạnh nhất trong kỳ triều cường mùa mưa khi lưu lượng dòng chảy đạt đỉnh, và tiêu tán gần như hoàn toàn trong các kỳ triều kém do năng lượng xáo trộn đáy suy giảm.

  2. Yếu tố nào đóng vai trò chính hình thành cực đại độ đục ở Luồng Tây và Luồng Đông?
    Tại Luồng Tây, hoàn lưu trọng lực do nước ngọt từ 3 cửa sông phía Tây gặp khối nước mặn 33 ppt cùng hiệu ứng bẫy triều là động lực chính. Trong khi đó, tại Luồng Đông, quá trình tái lơ lửng và lắng đọng bùn cát cục bộ do ứng suất dòng triều rút mạnh chiếm ưu thế vượt trội.

  3. Tương tác giữa sóng và dòng chảy ảnh hưởng thế nào đến nồng độ bùn cát?
    Tương tác sóng và dòng chảy làm tăng ứng suất cắt đáy tại các bãi cát ngầm từ 25% đến 50%, giúp mở rộng bề dày lõi bùn cát nồng độ cao thẳng đứng. Tuy nhiên, hiệu ứng sóng không làm thay đổi vị trí trung tâm hay biên độ dao động 22 km của cực đại độ đục.

  4. Nồng độ bùn cát lơ lửng trung bình và cấp hạt tại khu vực nghiên cứu là bao nhiêu?
    Nồng độ bùn cát lơ lửng trung bình toàn vùng đạt khoảng 0,172 g/l với hơn 92% khối lượng tập trung vào mùa lũ. Thành phần cấp hạt đáy dao động từ 2 đến 250 micromet, trong đó hạt lơ lửng có đường kính trung vị khoảng 21,6 micromet thuộc nhóm bùn sét kết dính dễ keo tụ.

  5. Vì sao cần áp dụng mô hình ba chiều thay vì mô hình hai chiều truyền thống?
    Mô hình 2 chiều tích phân độ sâu không thể mô phỏng hiện tượng phân tầng mật độ thẳng đứng và hoàn lưu trọng lực hai tầng với vận tốc dòng dư 0,3 m/s. Mô hình 3 chiều với 10 lớp sigma giải quyết triệt để cấu trúc trường vận tốc và nồng độ theo chiều sâu.

Kết luận

Luận văn đã mang lại những đóng góp khoa học và thực tiễn mang tính đột phá cho công tác nghiên cứu động lực học cửa sông:

  • Xây dựng thành công mô hình số trị ba chiều tích hợp tương tác sóng và dòng triều với lưới tính 10 lớp thẳng đứng, mô phỏng hoàn chỉnh quá trình động lực học bùn cát.
  • Xác định chính xác quy luật dịch chuyển 22 km của vùng cực đại độ đục và sự phân hóa cơ chế bồi lắng giữa Luồng Đông và Luồng Tây.
  • Định lượng được 80% trong tổng số 30,8 triệu tấn bùn cát hàng năm bị bẫy lại và bồi tụ ngay trong lòng vịnh Linh Đinh.
  • Làm rõ vai trò của trường dòng dư Eulerian với vận tốc cực đại 0,3 m/s và tác động gia tăng ứng suất cắt đáy do sóng tới 50%.
  • Cung cấp cơ sở khoa học tin cậy phục vụ tối ưu hóa kế hoạch nạo vét duy tu luồng hàng hải và quản trị môi trường ven biển.

Về lộ trình tiếp theo, các đơn vị quản lý cần tiến hành thử nghiệm mô hình cảnh báo bồi lắng thời gian thực trong 12 đến 24 tháng tới để kiểm chứng thực địa trên diện rộng. Hãy khai thác ngay bộ khung mô hình và các giải pháp kỹ thuật trong luận văn để nâng cao hiệu quả kinh tế và tính bền vững cho các dự án công trình biển.