Tổng quan nghiên cứu

Hệ thống sông ngòi và vùng ven biển Việt Nam với chiều dài đường bờ biển hơn 3.200 km đóng vai trò huyết mạch trong phát triển kinh tế - xã hội, song cũng là khu vực chịu tác động thủy văn phức tạp bậc nhất. Sự tương tác giữa dao động thủy triều biển khơi và chế độ dòng chảy ngọt từ thượng lưu tạo nên những biến động sâu sắc về mực nước, dòng chảy hai chiều và hiện tượng xâm nhập mặn. Thống kê thủy văn dọc bờ biển nước ta ghi nhận sự tồn tại của 8 vùng đặc trưng chế độ triều, nơi biên độ triều có sự chênh lệch lớn từ dưới 0,5 mét tại vùng biển miền Trung đến trên 3,5 mét tại vùng biển phía Nam và Bắc Bộ.

Vấn đề cốt lõi đặt ra là việc xác định các đặc trưng thủy văn thiết kế vùng cửa sông ảnh hưởng triều thường gặp nhiều sai số do tính phi tuyến của sóng triều và sự thiếu hụt tài liệu quan trắc đồng bộ. Đồng thời, sự phát triển mạnh mẽ của hệ thống hồ chứa nhân tạo ở thượng nguồn, điển hình như hồ chứa Hòa Bình với tổng dung tích 9,5 tỷ mét khối hay hồ chứa Sơn La với dung tích hiệu dụng 6,504 tỷ mét khối, đã làm thay đổi đáng kể chế độ dòng chảy tự nhiên và quá trình bồi lắng lòng dẫn.

Luận văn thạc sĩ này tập trung giải quyết bài toán tính toán thủy lực, thủy văn vùng sông chịu tác động của thủy triều và đánh giá quy luật điều tiết dòng chảy bằng hồ chứa. Phạm vi nghiên cứu bao quát các lưu vực sông ven biển trọng điểm, trạm thủy văn Trực Phương trên sông Ninh Cơ và các công trình hồ chứa quy mô lớn tại Việt Nam trong chuỗi số liệu quan trắc nhiều năm. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn, hỗ trợ giảm thiểu sai số tính toán mực nước thiết kế công trình xuống dưới 5%, tối ưu hóa hiệu quả cấp nước tưới tự chảy tăng 15% và đảm bảo an toàn xả lũ cho các vùng đồng bằng hạ du.

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 dựa trên sự tích hợp của hai nền tảng lý thuyết thủy văn và thủy lực công trình hiện đại:

Thứ nhất là lý thuyết động lực học thủy triều và tương tác thiên thể. Dao động thủy triều hình thành dưới tác động tương hỗ giữa trường hấp dẫn của Mặt trăng, Mặt trời và chuyển động tự quay của Trái đất xung quanh trục chung cách tâm Trái đất khoảng 0,73 bán kính Trái đất. Nghiên cứu áp dụng chỉ số hình dạng triều của Defant thông qua tỷ số giữa tổng các thành phần lực nhật triều và bán nhật triều: $F = (O_1 + K_1) / (S_2 + M_2)$. Khi chỉ số $F < 1$, chế độ bán nhật triều chiếm ưu thế; khi $F \ge 1$, chế độ nhật triều giữ vai trò chủ đạo.

Thứ hai là lý thuyết điều tiết dòng chảy và cân bằng bùn cát hồ chứa. Hồ chứa thực hiện việc khống chế, phân phối lại dòng chảy tự nhiên theo không gian và thời gian thông qua các thành phần dung tích được chuẩn hóa: dung tích chết $V_c$, dung tích hiệu dụng $V_h$, dung tích kết hợp $V_{kh}$ và dung tích siêu cao $V_{sc}$.

Khung lý thuyết định hình 5 khái niệm trọng tâm xuyên suốt: mực nước dâng bình thường ($H_{bt}$), nêm mặn và hiện tượng dòng chảy phân lớp hai chiều ($Q = Q^+ + Q^-$), điểm ngưng triều, mật độ bùn cát lơ lửng và quy luật bồi lắng lòng hồ.

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

Nguồn dữ liệu của luận văn được thu thập từ chuỗi quan trắc liên tục từ 25 đến 35 năm tại các trạm thủy văn cửa sông tiêu biểu như trạm Trực Phương và hệ thống hồ chứa lớn tại miền Bắc Việt Nam. Cỡ mẫu nghiên cứu bao gồm 30 tập mẫu số liệu cực trị mực nước đỉnh triều, chân triều hàng năm và hơn 300 đợt đo lưu lượng, bùn cát chi tiết qua các mùa lũ kiệt.

Phương pháp chọn mẫu áp dụng nguyên tắc mẫu cực trị hàng năm kết hợp chọn mẫu theo thời đoạn tính toán điển hình (chọn 1 giá trị mực nước chân triều thấp nhất hoặc đỉnh triều cao nhất trong các tháng then chốt như tháng 1 hoặc tháng 3 hàng năm). Phương pháp này bảo đảm tính độc lập thống kê, tính đồng nhất và tính đại diện của chuỗi biến ngẫu nhiên thủy văn.

Lý do lựa chọn phương pháp phân tích:

  • Phương pháp thống kê xác suất kết hợp kỹ thuật tịnh tiến mốc cốt địa hình: Khi vẽ đường tần suất mực nước, việc cộng hoặc trừ một hằng số cao độ $a$ hợp lý giúp ổn định hệ số phân tán $C_v$, loại bỏ sai số quy chiếu mốc cao đạc, từ đó xác định chính xác mực nước thiết kế $Z_p$ theo phân phối xác suất Pearson loại III.
  • Phương pháp mô hình toán thủy lực một chiều: Đối với các đoạn sông không có trạm đo thực nghiệm, nghiên cứu giải hệ phương trình Saint-Venant (gồm phương trình liên tục và phương trình động lực) bằng phương pháp sai phân hữu hạn thông qua các mô hình thủy lực chuyên dụng như MIKE 11, VRSAP. Phương pháp này giải quyết triệt để tính chất dòng chảy không ổn định và tương tác nêm mặn tại mạng lưới sông cửa bãi.

Toàn bộ quy trình thu thập dữ liệu, hiệu chỉnh thông số mô hình và phân tích thống kê được triển khai hoàn chỉnh trong timeline nghiên cứu kéo dài 18 tháng.

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

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

Nghiên cứu đã lượng hóa chi tiết tính phân hóa không gian của chế độ thủy triều dọc bờ biển Việt Nam. Tại khu vực từ Quảng Ninh đến Thanh Hóa, chế độ nhật triều đều chiếm ưu thế tuyệt đối với biên độ triều cường đạt từ 2,6 đến 3,2 mét. Ngược lại, khu vực cửa Thuận An ghi nhận biên độ triều thấp nhất cả nước, chỉ dao động từ 0,4 đến 0,5 mét trước khi gia tăng mạnh mẽ trở lại tại vùng biển Hàm Tân - Cà Mau với biên độ bán nhật triều không đều đạt 2,0 đến 3,5 mét.

Tại vùng sông ảnh hưởng triều, nghiên cứu làm rõ cơ chế truyền triều và hiện tượng phân lớp dòng chảy. Khi sóng triều truyền sâu vào lòng dẫn, do lực cản địa hình và sự đối kháng của dòng chảy ngọt, biên độ sóng triều giảm dần về 0 tại vị trí giới hạn triều. Trong các thời kỳ lũ lớn (như đợt quan trắc tháng 8 năm 1969 tại trạm Trực Phương), mực nước đỉnh triều bị đẩy cao vượt mức 300 cm, làm biến dạng đường quá trình mực nước hình sin tự nhiên. Hiện tượng phân lớp nêm mặn hình thành rõ nét khi tỷ trọng nước biển cao hơn nước ngọt, dẫn đến lớp nước ngọt tầng mặt tiếp tục chảy xuôi ra biển trong khi lớp nước mặn đáy chảy ngược lên thượng lưu.

Về phương diện hồ chứa, kết quả phân tích cân bằng bùn cát cho thấy lượng ngậm cát lơ lửng chiếm tỷ trọng áp đảo, trong khi lưu lượng bùn cát di đáy tại khu vực đồi núi chiếm từ 10% đến 30% tổng lượng bùn cát lơ lửng. Việc xác định dung tích chết $V_c \ge V_b \times T$ dựa trên công thức thực nghiệm Poliakov giúp kiểm soát chính xác tốc độ bồi lấp lòng hồ qua vòng đời khai thác 50 đến 100 năm.

Bố chế dung tích phòng lũ hạ du tại các công trình đầu mối lớn được xác lập với các chỉ số tối ưu: hồ Hòa Bình bố trí 4,69 tỷ mét khối (sau khi có hồ Sơn La được điều chỉnh về 3,0 tỷ mét khối), hồ Sơn La đảm nhận 4,0 tỷ mét khối phòng lũ với mực nước trước lũ duy trì ở cao trình 198,1 mét, và hồ Tuyên Quang duy trì dung tích phòng lũ 1,0 tỷ mét khối tại mực nước trước lũ 105,22 mét.

Thảo luận kết quả

Sự biến động biên độ triều giữa các vùng biển ven bờ bắt nguồn từ sự cộng hưởng địa hình vịnh nông và sự giao thoa sóng triều từ đại dương. Chiều cao sóng triều ngoài khơi chỉ đạt khoảng 1,0 mét nhưng khi tiến sâu vào các vịnh hẹp có thể bị biến dạng và tăng vọt lên 3,0 đến 5,0 mét. Trên các biểu đồ phân bố lưu tốc theo chiều sâu và biểu đồ quá trình lưu lượng $Q(t)$ lồng ghép mực nước $Z(t)$, sự tồn tại của các thời điểm ngưng triều ($Q = 0$) nhưng vẫn có dòng chảy phân lớp đã giải thích rõ nguyên nhân bùn cát mịn và độ mặn lắng đọng sâu vào nội đồng.

So sánh với các nghiên cứu thủy lực quốc tế, việc ứng dụng hệ phương trình Saint-Venant kết hợp phương trình cân bằng bùn cát một chiều đạt mức độ tương thích trên 92% so với số liệu đo thực tế sau 20 năm vận hành tại hồ Thác Bà (dung tích chết 780 triệu mét khối) và hồ Tuyên Quang (dung tích chết 561 triệu mét khối). Điều này chứng minh việc lựa chọn mực nước trước lũ hợp lý không chỉ giảm áp lực ngập lụt thượng lưu mà còn bảo đảm 100% dung tích trữ nước đầy hồ phục vụ tưới tiêu và phát điện vào cuối mùa lũ.

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

Nhằm nâng cao hiệu quả quản trị tài nguyên nước vùng cửa sông và khai thác bền vững các hồ chứa, luận văn đề xuất 4 nhóm giải pháp chiến lược:

Thứ nhất, nâng cấp mạng lưới quan trắc thủy văn tự động và kiểm soát triều thông minh. Bộ Nông nghiệp và Môi trường cần chủ trì lắp đặt trạm đo mực nước và độ mặn tự động theo thời gian thực tại 100% các cửa sông chính, hướng tới mục tiêu giảm 25% diện tích ngập úng và thiệt hại do xâm nhập mặn nội đồng trong giai đoạn 2026-2030.

Thứ hai, tối ưu hóa quy trình vận hành liên hồ chứa đa mục tiêu. Tập đoàn Điện lực Việt Nam và Ban quản lý lưu vực sông cần áp dụng thuật toán điều tiết thời gian thực, duy trì linh hoạt mực nước trước lũ tại hồ Sơn La (cao trình 198,1 mét) và hồ Hòa Bình (cao trình 88 đến 98,8 mét) nhằm cắt giảm từ 15% đến 20% lưu lượng đỉnh lũ xả về đồng bằng hạ du trước năm 2028.

Thứ ba, triển khai giải pháp kỹ thuật kiểm soát bồi lắng và xả cát định kỳ lòng hồ. Các đơn vị quản lý hồ chứa thủy điện lớn cần thực hiện đo đạc trắc dọc lòng hồ định kỳ 3 năm một lần, ứng dụng công nghệ xả cát đáy qua công trình tháo sâu để duy trì dung tích hiệu dụng, khống chế tỷ lệ bồi lấp dưới 0,8% tổng dung tích mỗi thập kỷ.

Thứ tư, hoàn thiện quy hoạch công trình lấy nước tưới nông nghiệp và đê biển. Chính quyền các tỉnh ven biển cần điều chỉnh lịch đóng mở cống tự chảy dựa trên dự báo chu kỳ triều kém và chân triều thấp nhất, bảo đảm nguồn nước tưới tiêu an toàn cho hơn 150.000 hecta đất canh tác nông nghiệp ven biển đến năm 2030.

Đố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 giá trị ứng dụng cao cho 4 nhóm đối tượng chuyên môn:

  1. Kỹ sư thiết kế và chuyên gia quy hoạch công trình thủy lợi, thủy điện: Tài liệu cung cấp quy trình giải tích chuẩn xác để tính toán mực nước triều thiết kế, xác định cao trình ngưỡng cống tiêu úng, cống ngăn mặn và phân bổ tối ưu các thành phần dung tích hồ chứa.
  2. Cơ quan quản lý nhà nước và ban chỉ huy phòng chống thiên tai: Cung cấp cơ sở khoa học để xây dựng kịch bản ngập lụt, lập bản đồ phân vùng xâm nhập mặn và ban hành quy trình vận hành liên hồ chứa giảm lũ an toàn cho hạ du.
  3. Ban quản lý khai thác công trình thủy điện và cấp nước: Hỗ trợ đánh giá chính xác tốc độ bồi lắng dung tích chết, dự báo tuổi thọ công trình và tối ưu hóa biểu đồ tích nước phát điện trong mùa kiệt.
  4. Giảng viên, học viên cao học và nghiên cứu sinh chuyên ngành Kỹ thuật tài nguyên nước, Thủy văn học: Luận văn là tài liệu tham khảo học thuật chuẩn mực, kết hợp nhuần nhuyễn giữa phương pháp giải tích lý thuyết, mô hình hóa thủy lực số và chuỗi dữ liệu thực nghiệm tại Việt Nam.

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

  1. Kỹ thuật tịnh tiến mốc cốt địa hình có vai trò gì trong việc tính toán mực nước triều thiết kế? Kỹ thuật này giải quyết hiện tượng hệ số phân tán $C_v$ bị sai lệch khi giá trị số đo mực nước thực tế quá lớn hoặc âm. Bằng cách cộng hoặc trừ một hằng số cao độ cố định $a$, chuỗi số liệu được quy về mốc chuẩn mới giúp giữ nguyên hệ số thiên lệch $C_s$, giảm thiểu sai số của đường tần suất Pearson loại III về mức dưới 3%.

  2. Sự khác biệt căn bản giữa triều phân lớp và triều không phân lớp tại vùng cửa sông là gì? Triều phân lớp xảy ra khi năng lượng sóng triều không đủ mạnh để xáo trộn hoàn toàn với dòng chảy ngọt thượng lưu, tạo ra nêm mặn hình tam giác ở đáy và hình thành dòng chảy hai chiều ngược hướng. Triều không phân lớp xuất hiện khi năng lượng triều rất lớn làm xáo trộn đồng nhất hai khối nước thành nêm mặn hình thang chuyển động cùng chiều.

  3. Tại sao dung tích kết hợp lại đóng vai trò quyết định trong việc giảm ngập lụt vùng thượng lưu hồ? Dung tích kết hợp được tạo ra bằng cách hạ mực nước hồ về mực nước trước lũ trong mùa mưa bão. Phần dung tích này giúp trữ tạm thời lượng nước lũ đột biến mà không cần nâng cao trình đập hoặc mực nước siêu cao, từ đó hạn chế tối đa diện tích đất đai và tài sản bị ngập lụt tại vùng lòng hồ thượng lưu.

  4. Phương pháp nào được ưu tiên để xác định mực nước thiết kế khi tuyến công trình không có số liệu quan trắc? Phương pháp mô hình toán một chiều giải hệ phương trình vi phân Saint-Venant bằng thuật toán sai phân hữu hạn là lựa chọn tối ưu. Mô hình sử dụng điều kiện biên trên là quá trình lưu lượng thiết kế $Q_p(t)$ và biên dưới là quá trình mực nước triều thiết kế $Z_p(t)$ để mô phỏng chính xác mực nước tại mọi mặt cắt.

  5. Tỷ lệ lưu lượng bùn cát di đáy được xác định như thế nào trong bài toán tính bồi lắng lòng hồ? Do việc đo đạc bùn cát di đáy rất phức tạp, lưu lượng này thường được xác định thông qua hệ số tỷ lệ kinh nghiệm so với bùn cát lơ lửng. Tại các vùng sông đồng bằng, bùn cát di đáy chiếm khoảng 5% đến 10% lượng cát lơ lửng, trong khi tại các lưu vực sông miền núi độ dốc lớn, tỷ lệ này dao động từ 10% đến 30%.

Kết luận

  • Luận văn đã hệ thống hóa và làm sáng tỏ quy luật phân hóa không gian của 8 vùng chế độ thủy triều dọc hơn 3.200 km bờ biển Việt Nam.
  • Thiết lập hoàn chỉnh phương pháp luận tính toán mực nước triều thiết kế thông qua kết hợp đường tần suất lý luận Pearson loại III và kỹ thuật xử lý mốc cốt địa hình.
  • Ứng dụng thành công hệ phương trình thủy lực Saint-Venant và phân tích động lực học nêm mặn để mô phỏng chính xác hiện tượng dòng chảy phân lớp hai chiều tại cửa sông.
  • Chuẩn hóa các nguyên tắc lựa chọn thành phần dung tích hồ chứa, dung tích kết hợp phòng lũ và quy trình tính toán bồi lắng bùn cát theo vòng đời công trình.
  • Cung cấp bộ thông số thực nghiệm thủy lực giá trị cao từ các đại công trình thủy điện Hòa Bình, Sơn La, Tuyên Quang và Thác Bà phục vụ phát triển hạ tầng nguồn nước.

Trong giai đoạn 2026-2030, các nghiên cứu tiếp theo cần phát triển mô hình toán thủy động lực học 2 chiều và 3 chiều tích hợp kịch bản nước biển dâng do biến đổi khí hậu. Các cơ quan quản lý và đơn vị tư vấn thiết kế hãy áp dụng ngay các phương pháp tính toán và giải pháp điều tiết trong luận văn để nâng cao tính chuẩn xác của các dự án công trình thủy lợi ven biển, bảo đảm an toàn dân sinh và phát triển kinh tế bền vững.