Vấn đề: lưu lượng xả thay đổi khi hồ tích nước
Một hồ điều hòa trong khu nhà máy làm chậm và giảm đỉnh dòng xả bằng cách giữ lại một phần nước tạm thời. Muốn tính mực nước cao nhất, cần cả thủy đồ vào lẫn quan hệ giữa cao độ, dung tích và lưu lượng ra. Chỉ lấy tổng lượng mưa trừ một lưu lượng xả cố định có thể bỏ qua việc cửa xả làm việc mạnh hơn khi nước trong hồ dâng lên.
Bài này xây dựng mô hình hồ có mặt nước ngang, không có thấm hoặc bốc hơi trong thời đoạn tính. USACE gọi mô hình này là level-pool routing, thường giải bằng phương pháp Modified Puls. [1] Ví dụ dùng một quan hệ cửa xả tuyến tính giả định để mọi bước có thể tính tay; không coi quan hệ ấy là công thức vật lý chung của lỗ xả, cống hay tràn.
Ba quan hệ cần biết trước khi lập bảng
Ký hiệu h là chiều cao nước so với mốc dung tích hoạt động, S(h) là dung tích m³ và O(h) là lưu lượng xả m³/s. Thủy đồ I(t) là lưu lượng vào m³/s theo thời gian. Ngoài ba quan hệ này, phải chọn trạng thái ban đầu h₀ hoặc S₀; hồ có sẵn nước sẽ khác hồ bắt đầu rỗng.
Trong ví dụ, chọn diện tích mặt bằng không đổi A=1.000 m² và S=Ah. Cửa xả giả định có O=kh, với k=0,50 m²/s. Đơn vị của k xuất phát từ m³/s chia cho m; nó không phải hệ số không thứ nguyên. Mốc h=0 đồng thời là trạng thái không còn dung tích hoạt động và O=0 của mô hình.
Với hồ thực có mái nghiêng, S(h) thường không tuyến tính. Lưu lượng cửa xả còn có thể thay đổi theo chế độ chảy và mức nước hạ lưu. HEC-HMS cho phép mô tả quan hệ cao độ–dung tích và ảnh hưởng ngập hạ lưu; các cao độ phải dùng cùng mốc. [2] Không thể lấy đường O(h) của cửa xả tự do rồi giữ nguyên khi cửa bị ngập đáng kể.
Cân bằng theo một bước thời gian
Từ bảo toàn thể tích, dS/dt=I−O. Dùng trung bình hai đầu bước thời gian Δt cho cả dòng vào và dòng ra:
S₂−S₁ = Δt·[(I₁+I₂)−(O₁+O₂)]/2
Chỉ số 1 là đầu bước, 2 là cuối bước; Δt tính bằng s. Phương trình này là dạng cân bằng rời rạc trong tài liệu USACE. [1] Điểm quan trọng là O₂ chưa biết vì nó phụ thuộc S₂ hoặc h₂. Bỏ O₂ rồi dùng O₁ cho cả bước tạo thành một phương pháp khác với đặc tính sai số khác.
Với S=Ah và O=kh, thay vào rồi chuyển vế cho công thức truy hồi:
h₂ = [(A−kΔt/2)h₁ + Δt(I₁+I₂)/2]/(A+kΔt/2)
Tử số có đơn vị m³, mẫu số m² nên h₂ ra m. Mỗi dòng bảng dùng h cuối dòng trước làm h đầu dòng sau; không khởi động lại từ h₀ ở mọi khoảng. Khi S(h) và O(h) phi tuyến, có thể lập quan hệ 2S/Δt+O theo h và nội suy hoặc giải nghiệm phù hợp.
Ví dụ giả định: một thủy đồ mưa kéo dài một giờ
Chọn Δt=600 s, h₀=0. Tại các thời điểm 0, 10, 20, 30, 40, 50 và 60 phút, lưu lượng vào lần lượt là 0; 0,50; 1,00; 0,50; 0,20; 0; 0 m³/s. Giữa các mốc, I được giả định biến thiên tuyến tính. Các giá trị không đại diện cho trận mưa tại một tỉnh hoặc chu kỳ lặp cụ thể.
Ở bước đầu, A−kΔt/2=850 m² và A+kΔt/2=1.150 m². h₁₀=[850×0+300×(0+0,50)]/1.150=0,130435 m. Dung tích cuối bước là 130,435 m³ và lưu lượng xả O₁₀=0,065217 m³/s.
Ở bước 10–20 phút, h₂₀=[850×0,130435+300×(0,50+1,00)]/1.150≈0,487713 m. Lặp cùng phép tính cho các bước còn lại tạo chuỗi kết quả sau; số liệu đã làm tròn để đọc, còn tính toán giữ thêm chữ số.
• 0 phút: h=0; S=0; O=0.
• 10 phút: h=0,130435 m; S=130,435 m³; O=0,065217 m³/s.
• 20 phút: h=0,487713 m; S=487,713 m³; O=0,243856 m³/s.
• 30 phút: h=0,751788 m; S=751,788 m³; O=0,375894 m³/s.
• 40 phút: h=0,738278 m; S=738,278 m³; O=0,369139 m³/s.
• 50 phút: h=0,597858 m; S=597,858 m³; O=0,298929 m³/s.
• 60 phút: h=0,441895 m; S=441,895 m³; O=0,220947 m³/s.
Trong bảng 10 phút, mực nước lớn nhất xuất hiện tại 30 phút, sau đỉnh dòng vào ở 20 phút. Đây là đỉnh tại các mốc tính, chưa chắc là đỉnh liên tục giữa hai mốc. Độ trễ xuất hiện vì hồ vẫn tiếp tục tích nước khi I lớn hơn O, dù I đã bắt đầu giảm.
Kiểm tra khép kín thể tích
Tổng dòng vào bằng tổng diện tích các hình thang dưới thủy đồ: 150+450+450+210+60+0=1.320 m³. Tổng dòng ra tính bằng cùng quy tắc hình thang là khoảng 878,105 m³. Dung tích cuối là 441,895 m³, nên S₆₀−S₀+Vra−Vvào≈0 trong sai số làm tròn.
Nếu lấy O cuối mỗi bước nhân cả 600 s rồi cộng, tổng ra sẽ khác vì đã thay quy tắc tích phân. Kiểm tra khối lượng cần dùng đúng quy tắc đã dùng để giải h. Tuy nhiên, khép kín khối lượng không tự chứng minh dự báo đỉnh chính xác: một mô hình với bước quá lớn vẫn có thể bảo toàn khối lượng rời rạc rất tốt.
Sau 60 phút, hồ chưa trở về trạng thái ban đầu. Nếu dừng bảng ở đây rồi coi trận tiếp theo bắt đầu rỗng, sẽ bỏ qua 441,895 m³ còn lưu lại. Cần kéo dài mô phỏng giai đoạn rút nước hoặc dùng lượng nước còn lại làm điều kiện đầu của sự kiện tiếp theo.
Độ nhạy bước thời gian và hình học
Chia bước xuống 300 s, nội suy tuyến tính cùng thủy đồ vào và giữ nguyên A, k. Đỉnh trên lưới mới khoảng 0,764187 m tại 35 phút, cao hơn đỉnh bảng 600 s khoảng 1,65%. Riêng tại 30 phút, lưới nhỏ cho 0,750271 m, rất gần 0,751788 m. Sai khác đỉnh chủ yếu liên quan việc lưới 10 phút không quan sát được thời điểm 35 phút, cùng với sai số tích phân.
So sánh này là một kiểm tra hội tụ bước đầu, không chứng minh 300 s đủ cho mọi hồ. Khi cửa tràn bắt đầu hoạt động, bơm đổi trạng thái hoặc thủy đồ vào nhọn, cần phân giải được quá trình chuyển tiếp. Hướng dẫn HEC-HMS mô tả việc giảm bước nội bộ khi mực nước, dung tích hoặc dòng ra thay đổi nhanh. [2]
Nếu tăng A từ 1.000 lên 1.500 m², giữ k và bước 600 s, đỉnh trên lưới giảm còn khoảng 0,582774 m tại 40 phút. Nhưng dung tích tương ứng là 874,162 m³, lớn hơn 751,788 m³ ở phương án đầu. Hồ rộng hơn làm mực nước thấp hơn đồng thời có thể giữ lại nhiều nước hơn; không nên so sánh chỉ một chỉ tiêu chiều sâu.
Nếu giữ A=1.000 m² nhưng giảm k từ 0,50 xuống 0,35 m²/s, đỉnh trên lưới tăng lên 0,859180 m tại 40 phút. Lưu lượng xả tại đó khoảng 0,300713 m³/s. Hạn chế cửa xả có thể giúp giảm tải hạ lưu nhưng làm yêu cầu trữ nước trong hồ lớn hơn.
Giới hạn và cách dùng kết quả
Mô hình giả định mặt nước ngang và một quan hệ đơn trị giữa trạng thái trữ với xả. USACE chỉ ra quan hệ trữ–xả có thể có vòng trễ đáng kể khi áp dụng cho dòng chảy trong lòng dẫn rất thoải. [3] Vì vậy, không nên dùng bài toán hồ đơn giản để thay thế mô hình truyền lũ trên mọi tuyến mương hoặc vùng ngập dài.
Trong mô hình tuyến tính, khi dòng vào bằng không, hệ số truy hồi là (A−kΔt/2)/(A+kΔt/2). Nếu chọn Δt>2A/k, hệ số trở thành âm và có thể sinh mức nước luân phiên âm–dương dù quá trình rút nước thực không như vậy. Kiểm tra dấu và miền vật lý bổ sung cho kiểm tra khối lượng; bước thỏa điều kiện dấu vẫn cần được xét về độ chính xác.
Dung tích tính từ thủy đồ chưa bao gồm mọi phần của hồ xây dựng. Phần dự trữ, bùn lắng, hình học mái, cao độ tràn và khoảng trống phải được xác định theo yêu cầu công trình. Một báo cáo có thể kiểm tra nên lưu cả thủy đồ, đường S(h), O(h), mốc cao độ, trạng thái đầu và bước thời gian. Các dữ liệu ấy cho biết vì sao hồ dâng đến một mức, thay vì chỉ đưa ra một thể tích tổng không thể tái dựng.
Nguồn tham khảo
[1] HEC-HMS Technical Reference Manual — Reservoir Modeling Concepts and Equations — U.S. Army Corps of Engineers, Hydrologic Engineering Center. Truy cập ngày 09/11/2026.
[2] HEC-HMS User’s Manual — Outflow Structures Routing — U.S. Army Corps of Engineers, Hydrologic Engineering Center. Truy cập ngày 09/11/2026.
[3] River Routing with HEC-1 and HEC-2, Training Document TD-30 — U.S. Army Corps of Engineers, Hydrologic Engineering Center — Vernon R. Bonner. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Một hồ điều hòa trong khu nhà máy làm chậm và giảm đỉnh dòng xả bằng cách giữ lại một phần nước tạm thời. Muốn tính mực nước cao nhất, cần cả thủy đồ vào lẫn quan hệ giữa cao độ, dung tích và lưu lượng ra. Chỉ lấy tổng lượng mưa trừ một lưu lượng xả cố định có thể bỏ qua việc cửa xả làm việc mạnh hơn khi nước trong hồ dâng lên.
Bài này xây dựng mô hình hồ có mặt nước ngang, không có thấm hoặc bốc hơi trong thời đoạn tính. USACE gọi mô hình này là level-pool routing, thường giải bằng phương pháp Modified Puls. [1] Ví dụ dùng một quan hệ cửa xả tuyến tính giả định để mọi bước có thể tính tay; không coi quan hệ ấy là công thức vật lý chung của lỗ xả, cống hay tràn.
Ba quan hệ cần biết trước khi lập bảng
Ký hiệu h là chiều cao nước so với mốc dung tích hoạt động, S(h) là dung tích m³ và O(h) là lưu lượng xả m³/s. Thủy đồ I(t) là lưu lượng vào m³/s theo thời gian. Ngoài ba quan hệ này, phải chọn trạng thái ban đầu h₀ hoặc S₀; hồ có sẵn nước sẽ khác hồ bắt đầu rỗng.
Trong ví dụ, chọn diện tích mặt bằng không đổi A=1.000 m² và S=Ah. Cửa xả giả định có O=kh, với k=0,50 m²/s. Đơn vị của k xuất phát từ m³/s chia cho m; nó không phải hệ số không thứ nguyên. Mốc h=0 đồng thời là trạng thái không còn dung tích hoạt động và O=0 của mô hình.
Với hồ thực có mái nghiêng, S(h) thường không tuyến tính. Lưu lượng cửa xả còn có thể thay đổi theo chế độ chảy và mức nước hạ lưu. HEC-HMS cho phép mô tả quan hệ cao độ–dung tích và ảnh hưởng ngập hạ lưu; các cao độ phải dùng cùng mốc. [2] Không thể lấy đường O(h) của cửa xả tự do rồi giữ nguyên khi cửa bị ngập đáng kể.
Cân bằng theo một bước thời gian
Từ bảo toàn thể tích, dS/dt=I−O. Dùng trung bình hai đầu bước thời gian Δt cho cả dòng vào và dòng ra:
S₂−S₁ = Δt·[(I₁+I₂)−(O₁+O₂)]/2
Chỉ số 1 là đầu bước, 2 là cuối bước; Δt tính bằng s. Phương trình này là dạng cân bằng rời rạc trong tài liệu USACE. [1] Điểm quan trọng là O₂ chưa biết vì nó phụ thuộc S₂ hoặc h₂. Bỏ O₂ rồi dùng O₁ cho cả bước tạo thành một phương pháp khác với đặc tính sai số khác.
Với S=Ah và O=kh, thay vào rồi chuyển vế cho công thức truy hồi:
h₂ = [(A−kΔt/2)h₁ + Δt(I₁+I₂)/2]/(A+kΔt/2)
Tử số có đơn vị m³, mẫu số m² nên h₂ ra m. Mỗi dòng bảng dùng h cuối dòng trước làm h đầu dòng sau; không khởi động lại từ h₀ ở mọi khoảng. Khi S(h) và O(h) phi tuyến, có thể lập quan hệ 2S/Δt+O theo h và nội suy hoặc giải nghiệm phù hợp.
Ví dụ giả định: một thủy đồ mưa kéo dài một giờ
Chọn Δt=600 s, h₀=0. Tại các thời điểm 0, 10, 20, 30, 40, 50 và 60 phút, lưu lượng vào lần lượt là 0; 0,50; 1,00; 0,50; 0,20; 0; 0 m³/s. Giữa các mốc, I được giả định biến thiên tuyến tính. Các giá trị không đại diện cho trận mưa tại một tỉnh hoặc chu kỳ lặp cụ thể.
Ở bước đầu, A−kΔt/2=850 m² và A+kΔt/2=1.150 m². h₁₀=[850×0+300×(0+0,50)]/1.150=0,130435 m. Dung tích cuối bước là 130,435 m³ và lưu lượng xả O₁₀=0,065217 m³/s.
Ở bước 10–20 phút, h₂₀=[850×0,130435+300×(0,50+1,00)]/1.150≈0,487713 m. Lặp cùng phép tính cho các bước còn lại tạo chuỗi kết quả sau; số liệu đã làm tròn để đọc, còn tính toán giữ thêm chữ số.
• 0 phút: h=0; S=0; O=0.
• 10 phút: h=0,130435 m; S=130,435 m³; O=0,065217 m³/s.
• 20 phút: h=0,487713 m; S=487,713 m³; O=0,243856 m³/s.
• 30 phút: h=0,751788 m; S=751,788 m³; O=0,375894 m³/s.
• 40 phút: h=0,738278 m; S=738,278 m³; O=0,369139 m³/s.
• 50 phút: h=0,597858 m; S=597,858 m³; O=0,298929 m³/s.
• 60 phút: h=0,441895 m; S=441,895 m³; O=0,220947 m³/s.
Trong bảng 10 phút, mực nước lớn nhất xuất hiện tại 30 phút, sau đỉnh dòng vào ở 20 phút. Đây là đỉnh tại các mốc tính, chưa chắc là đỉnh liên tục giữa hai mốc. Độ trễ xuất hiện vì hồ vẫn tiếp tục tích nước khi I lớn hơn O, dù I đã bắt đầu giảm.
Kiểm tra khép kín thể tích
Tổng dòng vào bằng tổng diện tích các hình thang dưới thủy đồ: 150+450+450+210+60+0=1.320 m³. Tổng dòng ra tính bằng cùng quy tắc hình thang là khoảng 878,105 m³. Dung tích cuối là 441,895 m³, nên S₆₀−S₀+Vra−Vvào≈0 trong sai số làm tròn.
Nếu lấy O cuối mỗi bước nhân cả 600 s rồi cộng, tổng ra sẽ khác vì đã thay quy tắc tích phân. Kiểm tra khối lượng cần dùng đúng quy tắc đã dùng để giải h. Tuy nhiên, khép kín khối lượng không tự chứng minh dự báo đỉnh chính xác: một mô hình với bước quá lớn vẫn có thể bảo toàn khối lượng rời rạc rất tốt.
Sau 60 phút, hồ chưa trở về trạng thái ban đầu. Nếu dừng bảng ở đây rồi coi trận tiếp theo bắt đầu rỗng, sẽ bỏ qua 441,895 m³ còn lưu lại. Cần kéo dài mô phỏng giai đoạn rút nước hoặc dùng lượng nước còn lại làm điều kiện đầu của sự kiện tiếp theo.
Độ nhạy bước thời gian và hình học
Chia bước xuống 300 s, nội suy tuyến tính cùng thủy đồ vào và giữ nguyên A, k. Đỉnh trên lưới mới khoảng 0,764187 m tại 35 phút, cao hơn đỉnh bảng 600 s khoảng 1,65%. Riêng tại 30 phút, lưới nhỏ cho 0,750271 m, rất gần 0,751788 m. Sai khác đỉnh chủ yếu liên quan việc lưới 10 phút không quan sát được thời điểm 35 phút, cùng với sai số tích phân.
So sánh này là một kiểm tra hội tụ bước đầu, không chứng minh 300 s đủ cho mọi hồ. Khi cửa tràn bắt đầu hoạt động, bơm đổi trạng thái hoặc thủy đồ vào nhọn, cần phân giải được quá trình chuyển tiếp. Hướng dẫn HEC-HMS mô tả việc giảm bước nội bộ khi mực nước, dung tích hoặc dòng ra thay đổi nhanh. [2]
Nếu tăng A từ 1.000 lên 1.500 m², giữ k và bước 600 s, đỉnh trên lưới giảm còn khoảng 0,582774 m tại 40 phút. Nhưng dung tích tương ứng là 874,162 m³, lớn hơn 751,788 m³ ở phương án đầu. Hồ rộng hơn làm mực nước thấp hơn đồng thời có thể giữ lại nhiều nước hơn; không nên so sánh chỉ một chỉ tiêu chiều sâu.
Nếu giữ A=1.000 m² nhưng giảm k từ 0,50 xuống 0,35 m²/s, đỉnh trên lưới tăng lên 0,859180 m tại 40 phút. Lưu lượng xả tại đó khoảng 0,300713 m³/s. Hạn chế cửa xả có thể giúp giảm tải hạ lưu nhưng làm yêu cầu trữ nước trong hồ lớn hơn.
Giới hạn và cách dùng kết quả
Mô hình giả định mặt nước ngang và một quan hệ đơn trị giữa trạng thái trữ với xả. USACE chỉ ra quan hệ trữ–xả có thể có vòng trễ đáng kể khi áp dụng cho dòng chảy trong lòng dẫn rất thoải. [3] Vì vậy, không nên dùng bài toán hồ đơn giản để thay thế mô hình truyền lũ trên mọi tuyến mương hoặc vùng ngập dài.
Trong mô hình tuyến tính, khi dòng vào bằng không, hệ số truy hồi là (A−kΔt/2)/(A+kΔt/2). Nếu chọn Δt>2A/k, hệ số trở thành âm và có thể sinh mức nước luân phiên âm–dương dù quá trình rút nước thực không như vậy. Kiểm tra dấu và miền vật lý bổ sung cho kiểm tra khối lượng; bước thỏa điều kiện dấu vẫn cần được xét về độ chính xác.
Dung tích tính từ thủy đồ chưa bao gồm mọi phần của hồ xây dựng. Phần dự trữ, bùn lắng, hình học mái, cao độ tràn và khoảng trống phải được xác định theo yêu cầu công trình. Một báo cáo có thể kiểm tra nên lưu cả thủy đồ, đường S(h), O(h), mốc cao độ, trạng thái đầu và bước thời gian. Các dữ liệu ấy cho biết vì sao hồ dâng đến một mức, thay vì chỉ đưa ra một thể tích tổng không thể tái dựng.
Nguồn tham khảo
[1] HEC-HMS Technical Reference Manual — Reservoir Modeling Concepts and Equations — U.S. Army Corps of Engineers, Hydrologic Engineering Center. Truy cập ngày 09/11/2026.
[2] HEC-HMS User’s Manual — Outflow Structures Routing — U.S. Army Corps of Engineers, Hydrologic Engineering Center. Truy cập ngày 09/11/2026.
[3] River Routing with HEC-1 and HEC-2, Training Document TD-30 — U.S. Army Corps of Engineers, Hydrologic Engineering Center — Vernon R. Bonner. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
