Với đầu vào: st lượng trữ thực tế, qt độ sâu dòng chảy mặt thượng lưu, pt độ
sâu mưa, etđộ sâu bốc hơi
Lưu lượng vào được hiệu chỉnh theo công thức:
it = qt (a - at) + pt at - et at (1) trong đó: at là diện tích mặt hồ
a là diện tích lưu vực thực tế
Lượng tổn thất do rò rỉ, bốc hơi, ước lượng theo công thức:
lt (zt) = α zt3 + β zt2 + γ zt + δ (2) trong đó: zt là mực nước thực tế
α, β, γ, δ: là hằng số, thay đổi theo tháng (vì tổn thất có thể biến động theo mùa của mực nước ngầm).
Lượng nước thực, trước khi xả (i = 0), được cho bởi:
s0t = st + it - lt (3)
Ngoài ra, lượng xả từ hồ chứa được thực hiện theo phân cấp sử dụng nước (đối với những sử dụng i = 1, ..., n). Giả thiết, nhu cầu sử dụng nước thứ i, lượng xả
mong muốn tương ứng được xác định trên cơ sở mực nước thực tế và nhu cầu thực tế, nghĩa là:
rit* = λij dit, với λij = f (zt) (4) Lượng xả thực tế không được vượt quá dung tích hữu ích thực:
rit = min (rit*, si-1, t - smin) (5) Khi đó, lượng nước có trong hồ chứa (lượng trữ hồ chứa):
sit = si-1, t - rit (6) Sau đó, tính toán lại mực nước hồ chứa, tương ứng với các lượng trữ mới. Khi tất cả các nhu cầu được đáp ứng, lượng tràn ước tính bằng:
Bắt đầu t=1
Cho: St, qt, pt, ct, Lượng vào: It=qt(a-at)+ptat-ctat
Tổn thất: et(zt)=azt3+bzt2+gzt+d
Lượng trữ trước khi xả: S0t=St+it-et Nhu cầu nước di, Quy tắc điều hành:li,j=f(zi,j)
i,j=1
j=j+1
Lượng xả cho yêu cầu i, ứng với zj,j:
R*
i,j=li,jdi,j
(lượng xả thực): Rịt=min(R*
i,j, Si-1,t-Smin) Lượng trữ hồ cập nhât: Sit=Si-1,t-Rit Mực nước hồ: zit=f(Sit) i=i+1 Lượng xả tổng: Wt=max(Snt-Smax) Qt t=t+1 Kiểm tra số cặp (li,j zi,j) Kiểm tra số sử dụng i Kết quả Qxả=Qt =f(t) Kiểm tra số ngày tính t Kết thúc