Tổng quan nghiên cứu

Chi phí đầu tư xây dựng các hạng mục công trình tháo lũ thường chiếm từ 25% đến 40% tổng nguồn vốn xây lắp của một dự án đầu mối hồ chứa thủy lợi. Đặc thù của đường tràn xả lũ là dòng chảy luôn ở trạng thái chảy xiết với số Froude lớn hơn 1, thậm chí vượt ngưỡng 2.5 và vận tốc dòng chảy có thể đạt trên 30 m/s. Nhằm tối ưu hóa bài toán kinh tế kỹ thuật, giảm khối lượng đào đắp đất đá và bê tông thân dốc, giải pháp bố trí đoạn thu hẹp chuyển tiếp sau ngưỡng tràn được áp dụng phổ biến tại nhiều công trình như hồ Gò Miếu dài 125m, hồ Yên Lịch II dài 29m hay hồ Lộc Đại với chiều dài chuyển tiếp 30m. Tuy nhiên, việc thu hẹp lòng dẫn trên dòng xiết tạo ra các hiện tượng thủy lực phi tuyến rất phức tạp như sóng xiên, sóng phản xạ, nước nhảy và hiện tượng hàm khí. Phương pháp giải tích một chiều truyền thống bộc lộ nhược điểm lớn khi không thể mô phỏng được sự biến dạng mặt nước không gian hai chiều, dẫn đến sai số tính toán chiều cao tường biên từ 15% đến 25%. Nghiên cứu được thực hiện theo Quyết định số 2039/QĐ-ĐHTL ngày 27 tháng 10 năm 2015 của Hiệu trưởng Trường Đại học Thủy lợi, tập trung giải quyết bài toán mô phỏng dòng xiết trên dốc nước có đoạn thu hẹp đáy phẳng tường biên thẳng. Mục tiêu trọng tâm là ứng dụng phương pháp thể tích hữu hạn giải hệ phương trình nước nông hai chiều để mô phỏng chính xác đường mặt nước và trường lưu tốc, áp dụng tính toán thực tế cho công trình tràn xả lũ hồ chứa Hốc Xoài tại Quảng Ngãi và hồ chứa Khe Gia tại Quảng Ninh. Đề tài mang lại giá trị thực tiễn cao, giúp nâng cao độ an toàn công trình đầu mối, tiết kiệm từ 12% đến 18% chi phí gia cố kết cấu và ngăn ngừa nguy cơ xói lở nguy hiểm ở hạ lư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ủa hệ phương trình nước nông hai chiều (2D Shallow Water Equations - SWE) dạng hyperbolic phi tuyến kết hợp cùng lý thuyết sóng gián đoạn trong cơ học chất lỏng. Hệ phương trình nước nông 2 chiều phản ánh đầy đủ quy luật bảo toàn khối lượng và bảo toàn động lượng theo hai trục tọa độ không gian x, y cùng biến thời gian t, có xét đến các số hạng nguồn gồm độ dốc đáy $S_{0x}, S_{0y}$ và lực cản ma sát đáy qua hệ số nhám Manning $n$. Khung lý thuyết đi sâu phân tích bài toán Riemann tổng quát về hiện tượng vỡ đập cùng hai nghiệm giải tích chuẩn mực là bài toán Stoker cho điều kiện hạ lưu ướt và bài toán Ritter cho điều kiện hạ lưu khô. Các khái niệm thủy lực then chốt được lượng hóa bao gồm: số Froude ($Fr = V / \sqrt{gh}$ với trạng thái chảy xiết khi $Fr > 1$), tốc độ lan truyền sóng trọng trường nhỏ ($c = \sqrt{gh}$), tỷ năng mặt cắt ($\ni = h + \alpha V^2 / 2g$), độ sâu phân giới ($h_k$) cùng điều kiện sốc Rankine-Hugoniot và điều kiện entropy Lax nhằm nhận diện chính xác sóng nén (shock) và sóng dãn (rarefaction).

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

Nghiên cứu thu thập và xử lý toàn bộ hồ sơ kỹ thuật, thông số địa hình và thủy văn từ 02 công trình tràn xả lũ thực tế gồm hồ Hốc Xoài (Quảng Ngãi) và hồ Khe Gia (Quảng Ninh), kết hợp với chuỗi dữ liệu thử nghiệm mô phỏng trên kênh lăng trụ đáy bằng dài 2000m. Để đánh giá tính ổn định và độ tin cậy của thuật toán, tác giả thiết lập cỡ mẫu kiểm chứng gồm 3 kịch bản mực nước hạ lưu khác nhau trong bài toán Stoker ($h_0 = 1.0\text{m}$, $h_0 = 5.0\text{m}$, $h_0 = 0.1\text{m}$ trên nền mực nước thượng lưu $h_1 = 10.0\text{m}$) và 3 mốc thời gian tính toán ($t = 40\text{s}$, $50\text{s}$, $80\text{s}$). Phương pháp thể tích hữu hạn (Finite Volume Method - FVM) với sơ đồ xấp xỉ Godunov được lựa chọn nhờ tính năng "bắt sốc" vượt trội và khả năng bảo toàn cơ học tuyệt đối tại các mặt cắt gián đoạn, khắc phục hoàn toàn hiện tượng mất ổn định số thường gặp ở phương pháp sai phân hữu hạn (FDM) hoặc phần tử hữu hạn (FEM). Miền tính toán được rời rạc hóa linh hoạt với các cấp ô lưới từ thô $\Delta x = \Delta y = 5.0\text{m}$ đến lưới mịn $\Delta x = \Delta y = 0.1\text{m}$. Toàn bộ quy trình tính toán, so sánh đối chuẩn giữa phương pháp số và phương pháp cộng trực tiếp giải tích 1 chiều được thực hiện trong giai đoạn từ cuối năm 2015 đến năm 2016 tại Bộ môn Thủy công, Trường Đại học Thủy lợi.

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

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

Thứ nhất, kết quả kiểm chứng mô hình trên kênh tiêu chuẩn dài 2000m chứng minh phương pháp thể tích hữu hạn FVM cho kết quả trùng khớp hoàn hảo với nghiệm giải tích chính xác của bài toán Ritter và Stoker, với sai số mực nước tại thời điểm $t = 50\text{s}$ nhỏ hơn 0.5% và lưu lượng đơn vị cực đại tại vị trí tim đập đạt mức $24.7\text{ m}^3/\text{s/m}$. Thứ hai, tại tràn xả lũ hồ Hốc Xoài và hồ Khe Gia, phương pháp số 2 chiều đã mô phỏng rõ nét chuỗi sóng xiên hình thành từ vị trí tường biên chuyển tiếp, gây ra hiện tượng mực nước dâng cao không đối xứng giữa hai bên thành dốc với độ chênh cao trình mặt nước cục bộ từ 1.2m đến 1.8m. Thứ ba, độ mịn của lưới tính toán có ảnh hưởng trực tiếp đến chất lượng kết quả: ô lưới thô $\Delta x = \Delta y = 5.0\text{m}$ chỉ phản ánh được xu hướng chung, trong khi ô lưới mịn $\Delta x = \Delta y = 0.1\text{m}$ bắt trọn vẹn từng đỉnh sóng phản xạ cục bộ và bộc lộ vận tốc dòng chảy thực tế vượt ngưỡng 28 m/s, cao hơn 15% so với kết quả tính theo sơ đồ 1 chiều.

Thảo luận kết quả

Sự xuất hiện của các dải sóng xiên có biên độ lớn trên đoạn dốc thu hẹp được giải thích là do đặc tính của dòng siêu xiết ($Fr > 2.5$), khi vận tốc dòng chảy lớn hơn vận tốc truyền sóng ($V > c$), các nhiễu động hình học từ tường biên không thể lan truyền ngược dòng mà bị dồn nén, va đập liên tục vào thành bên tạo thành chuỗi sóng phản xạ giao thoa kéo dài. Phương pháp cộng trực tiếp giải tích 1 chiều truyền thống mặc định đường mặt nước nằm ngang trên mỗi mặt cắt, bỏ qua hoàn toàn tương tác động lượng theo phương ngang y, do đó đã đánh giá thấp đáng kể chiều cao tường biên thực tế từ 18% đến 25%. Các kết quả nghiên cứu được trực quan hóa chi tiết thông qua các biểu đồ phân bố đường mặt nước không gian, bản đồ vector lưu tốc 2D và bảng số liệu thủy lực tại các mặt cắt đặc trưng dọc theo thân dốc. So sánh với các công bố của LeVeque và Toro, mô hình số sử dụng trong luận văn thể hiện tính ưu việt khi giải quyết mượt mà các điều kiện biên đáy khô - ướt mà không làm phát sinh dao động số giả tạo, mang lại độ tin cậy vượt trội cho công tác thiết kế hình học công trình tháo nước.

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

Một là, Cục Quản lý Đê điều và Phòng chống Thiên tai cùng các cơ quan thẩm định chuyên ngành cần ban hành quy định bắt buộc ứng dụng mô hình số 2 chiều (2D SWE) trong việc kiểm tra, thẩm tra thiết kế thủy lực cho các dốc nước tràn xả lũ có vận tốc dòng chảy $V > 15\text{ m/s}$, nhằm loại bỏ hoàn toàn nguy cơ tràn nước qua đỉnh tường biên, lộ trình hoàn thành trong vòng 12 tháng.

Hai là, các đơn vị tư vấn thiết kế công trình thủy lợi cần tối ưu hóa hình học đoạn thu hẹp bằng cách áp dụng tường biên cong liên hợp hoặc đường cong không nhiễu thay cho tường gãy khúc đột ngột, giúp giảm biên độ sóng xiên từ 20% đến 35%, qua đó tiết kiệm khoảng 10% khối lượng bê tông tường biên ngay từ giai đoạn lập dự án đầu tư.

Ba là, chủ đầu tư và nhà thầu xây lắp cần triển khai lắp đặt hệ thống mố nhám gia cường hoặc mũi phun tiêu năng tại các vị trí dốc nước có độ dốc đáy $i \ge 10%$, nhằm chủ động tiêu hao từ 40% đến 55% động năng dư thừa trước khi dòng chảy vào bể tiêu năng, đảm bảo tuổi thọ công trình an toàn trên 50 năm và chống xói lở bờ hạ lưu.

Bốn là, các viện nghiên cứu và trường đại học kỹ thuật cần chuẩn hóa phần mềm tính toán thủy lực chuyên dụng với kích thước lưới tính tối ưu $\Delta x \le 0.1\text{m}$, giúp rút ngắn thời gian mô phỏng phương án từ 10 ngày xuống dưới 24 giờ, hoàn thành việc chuyển giao công nghệ cho các kỹ sư thủy lợi trong thời gian 18 tháng tới.

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

Nhóm kỹ sư thiết kế công trình thủy lợi và thủy điện: Tiếp cận quy trình tính toán thủy lực 2 chiều tiên tiến, hỗ trợ xác định chính xác cao trình tường biên và tối ưu hóa kết cấu nối tiếp tiêu năng cho các đường tràn có đoạn thu hẹp phức tạp.

Nhóm học viên cao học, nghiên cứu sinh chuyên ngành Xây dựng Công trình Thủy: Sử dụng tài liệu như một cẩm nang học thuật chuyên sâu về phương pháp thể tích hữu hạn, thuật toán Godunov và kỹ thuật mô phỏng các bài toán sóng vỡ đập Riemann.

Nhóm cán bộ quản lý, giám sát và vận hành hồ chứa tại các Ban Quản lý Dự án: Nắm bắt các nguy cơ thủy lực tiềm ẩn như sóng xiên, nước nhảy và khí thực để xây dựng phương án duy tu, gia cố các tràn xả lũ hiện hữu trước mùa mưa bão.

Nhóm giảng viên và sinh viên ngành Kỹ thuật Tài nguyên Nước: Sử dụng làm tài liệu tham khảo giảng dạy trực quan cho các môn học Thủy lực công trình và Thiết kế đập tràn, kết nối giữa lý thuyết giải tích cổ điển và mô hình số hiện đại.

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

Tại sao phương pháp giải tích 1 chiều truyền thống không còn phù hợp cho dốc nước có đoạn thu hẹp? Phương pháp giải tích 1 chiều chỉ tính toán dựa trên giả thiết dòng chảy ổn định biến đổi dần với mặt nước nằm ngang. Trong thực tế dòng xiết ($Fr > 2.5$) qua đoạn thu hẹp xuất hiện sóng xiên và sóng va đập 2 chiều khiến mực nước dâng cao tại tường biên, phương pháp 1 chiều sẽ bỏ sót hiện tượng này và gây thiếu hụt chiều cao tường từ 18% đến 25%.

Sơ đồ thể tích hữu hạn Godunov mang lại ưu thế gì vượt trội so với các phương pháp số khác? Sơ đồ Godunov thuộc hệ phương pháp "bắt sốc", có khả năng bảo toàn cơ học tuyệt đối khối lượng và động lượng qua từng ô lưới. Thuật toán này xử lý chính xác các nghiệm gián đoạn như sóng shock và sóng rarefaction mà không gây hiện tượng phân kỳ hay nhiễu số như phương pháp sai phân hữu hạn thông thường.

Việc phân chia kích thước ô lưới ảnh hưởng như thế nào đến độ chính xác của bài toán dòng xiết? Kích thước ô lưới quyết định độ phân giải của các gợn sóng phản xạ. Khi chia lưới thô $\Delta x = \Delta y = 5.0\text{m}$, mô hình chỉ cho giá trị trung bình; nhưng khi giảm kích thước xuống $\Delta x = \Delta y = 0.1\text{m}$, mô hình nắm bắt trọn vẹn các đỉnh sóng cục bộ va vào thành bên, giúp sai số tính toán giảm xuống dưới 1%.

Ý nghĩa của việc giải bài toán Stoker và Ritter trước khi áp dụng vào công trình thực tế là gì? Bài toán Stoker (hạ lưu ướt) và Ritter (hạ lưu khô) là các trường hợp sóng gián đoạn hiếm hoi có nghiệm giải tích chính xác. Việc mô phỏng kênh dài 2000m tại các mốc thời gian 40s, 50s, 80s là bước đối chuẩn bắt buộc nhằm khẳng định độ tin cậy tuyệt đối của mã nguồn chương trình tính toán trước khi áp dụng cho địa hình thực địa.

Giải pháp thu hẹp lòng dẫn trên dốc nước mang lại hiệu quả kinh tế ra sao cho công trình? Việc thu hẹp lòng dẫn dốc nước giúp thu nhỏ bề rộng thân dốc, từ đó giảm từ 15% đến 30% khối lượng đào phá đất đá và vật liệu bê tông. Tuy nhiên, giải pháp này bắt buộc phải đi kèm tính toán mô phỏng 2 chiều chuẩn xác để đảm bảo dòng xiết không tràn qua tường gây xói lở thân đập.

Kết luận

  • Luận văn làm sáng tỏ các đặc tính thủy lực phi tuyến 2 chiều của dòng chảy xiết ($Fr > 2.5$) trên dốc nước có bố trí đoạn thu hẹp đáy phẳng tường thẳng.
  • Ứng dụng thành công phương pháp thể tích hữu hạn FVM giải hệ phương trình nước nông 2 chiều bảo toàn năng lượng và động lượng.
  • Kiểm chứng độ tin cậy của thuật toán đạt mức chính xác trên 99.5% thông qua các bài toán sóng gián đoạn kinh điển Stoker và Ritter trên kênh 2000m.
  • Tính toán thành công cho tràn xả lũ hồ chứa Hốc Xoài và Khe Gia, xác định chính xác cao trình dâng sóng va đập và khắc phục sai số từ 18% đến 25% của phương pháp 1 chiều.
  • Đề xuất các giải pháp công trình hữu hiệu như tối ưu hóa góc lệch tường biên và bố trí mố nhám tiêu năng để bảo vệ an toàn cho tuyến dốc.

Đóng góp then chốt của luận văn là đã thiết lập một công cụ tính toán số hiện đại, thay thế phương pháp giải tích truyền thống còn nhiều hạn chế trong ngành thủy công Việt Nam. Định hướng nghiên cứu tiếp theo trong giai đoạn 12 đến 24 tháng tới là tích hợp mô hình số 3 chiều nhằm phân tích sâu hiện tượng hàm khí và áp lực mạch động lên đáy dốc nước. Các kỹ sư và chuyên gia thủy lợi hãy áp dụng ngay phương pháp số này vào công tác tính toán thiết kế tràn xả lũ để đảm bảo an toàn tuyệt đối cho các công trình hồ chứa!