Tại các lưu vực sông có trạm đo thủy văn với thời gian đủ dài, chế độ thủy văn được xác định theo số liệu thực đo. Tại các lưu vực sông không có trạm đo, hoặc có trạm nhưng thời gian đo đạc ngắn thì việc đánh giá chế độ dòng chảy phải dựa vào các phương pháp khác. Các phương pháp được sử dụng tùy thuộc vào nhân tố cần xác định. Ví dụ để xác định chuẩn dòng chảy năm, phân phối dòng chảy có thể dùng các phương pháp sau:
+ Lưu vực tương tự;
+ Bản đồ đẳng trị;
+ Nội suy;
+ Phương trình cân bằng nước;
+ Quan hệ giữa các thông số phân phối với các nhân tố ảnh hưởng (xây dựng cho từng vùng);
Ngày nay, với sự phát triển của khoa học kỹ thuật, các lưu vực không có trạm đo, hoặc có ít số liệu quan trắc, chuỗi số liệu dòng chảy thường được tính toán, khôi phục lại bằng các mô hình thủy văn, dựa trên có nhân tố hình thành và ảnh hưởng như mưa, bốc hơi, tính chất lưu vực, ...
Một trong những mô hình được sử dụng rộng rãi hiện nay cho các lưu vực sông vừa và nhỏ là mô hình NAM. Thông tin về mô hình này như sau:
Mô hình NAM được xây dựng tại Khoa Thuỷ văn Viện Kỹ thuật Thuỷ Động lực và Thuỷ lực thuộc Đại học Kỹ thuật Đan Mạch năm 1982. Mô hình Nam đã được sử dụng rộng rãi ở Đan Mạch và một số nước nằm trong nhiều vùng khí hậu khác nhau như Srilanca, Thailand, Ấn Độ, vv... và Việt Nam. Trong mô hình NAM, mỗi lưu vực được xem là một đơn vị xử lý. Do đó, các thông số và các biến là đại diện cho các giá trị được trung bình hóa trên toàn lưu vực. Mô hình tính quá trình mưa - dòng chảy theo cách tính liên tục hàm lượng ẩm trong năm bể chứa riêng biệt có tương tác lẫn nhau.
Cấu trúc mô hình NAM được xây dựng trên nguyên tắc các hồ chứa theo chiều thẳng đứng và các hồ chứa tuyến tính, gồm có 5 bể chứa theo chiều thẳng đứng như hình 10.1.
+ Bể chứa tuyết tan được kiểm soát bằng các điều kiện nhiệt độ. Đối với điều kiện khí hậu nhiệt đới ở nước ta thì không xét đến bể chứa này.
+ Bể chứa mặt: lượng nước ở bể chứa này bao gồm lượng nước mưa do lớp phủ thực vật chặn lại, lượng nước đọng lại trong các chỗ trũng và lượng nước trong tầng sát mặt. Giới hạn trên của bể chứa này được ký hiệu bằng Umax.
+ Bể chứa tầng dưới: là vùng đất có rễ cây nên cây cối có thể hút nước cho bốc, thoát hơi. Giới hạn trên của lượng nước trong bể chứa này được ký hiệu bằng Lmax, lượng nước hiện tại được ký hiệu là L và tỷ số L/Lmax biểu thị trạng thái ẩm của bể chứa.
+ Bể chứa nước ngầm tầng trên.
+ Bể chứa nước ngầm tầng dưới.
Mưa hoặc tuyết tan đều đi vào bể chứa mặt. Lượng nước (U) trong bể chứa mặt liên tục cung cấp cho bốc hơi và thấm ngang thành dòng chảy sát mặt. Khi U đạt đến Umax, lượng nước thừa là dòng chảy tràn trực tiếp ra sông và một phần còn lại sẽ thấm xuống các bể chứa tầng dưới và bể chứa ngầm.
Nước trong bể chứa tầng dưới liên tục cung cấp cho bốc thoát hơi và thấm xuống bể chứa ngầm. Lượng cấp nước ngầm được phân chia thành hai bể chứa: tầng trên và tầng dưới, hoạt động như các hồ chứa tuyến tính với các hằng số thời gian khác nhau. Hai bể chứa này liên tục chảy ra sông tạo thành dòng chảy gốc.
Dòng chảy tràn và dòng chảy sát mặt được diễn toán qua một hồ chứa tuyến tính thứ nhất, sau đó các thành phần dòng chảy được cộng lại và diễn toán qua hồ chứa tuyến tính thứ hai. Cuối cùng cũng thu được dòng chảy tổng cộng tại cửa ra.
![](/Portals/0/2020-huy/hinh10.1.jpg)
Hình 10.1. Cấu trúc của mô hình NAM
Thành phần cơ bản của mô hình
Lượng trữ bề mặt:
Lượng ẩm bị giữ lại bởi thực vật cũng như được trữ trong các chỗ trũng trên tầng trên cùng của bề mặt đất được coi là lượng trữ bề mặt. Umax biểu thị giới hạn trên của tổng lượng nước trong lượng trữ bề mặt. Tổng lượng nước U trong lượng trữ bề mặt liên tục bị giảm do bốc hơi cũng như do thấm ngang. Khi lượng trữ bề mặt đạt đến mức tối đa, một lượng nước thừa PN sẽ gia nhập vào sông với vai trò là dòng chảy tràn trong khi lượng còn lại sẽ thấm vào tầng thấp bên dưới và tầng ngầm.
Lượng trữ tầng thấp hay lượng trữ tầng rễ cây:
Độ ẩm trong tầng rễ cây, lớp đất bên dưới bề mặt đất, tại đó thực vật có thể hút nước để bốc thoát hơi đặc trưng cho lượng trữ tầng thấp. Lmax biểu thị giới hạn trên của tổng lượng nước trữ trong tầng này. Độ ẩm trong lượng trữ tầng thấp cung cấp cho bốc thoát hơi thực vật. Độ ẩm trong tầng này điều chỉnh tổng lượng nước gia nhập vào lượng trữ tầng ngầm, thành phần dòng chảy mặt, dòng sát mặt và lượng gia nhập lại.
Bốc thoát hơi nước:
Nhu cầu bốc thoát hơi đầu tiên được thoả mãn từ lượng trữ bề mặt với tốc độ tiềm năng. Nếu lượng ẩm U trong lượng trữ bề mặt nhỏ hơn yêu cầu (U < Ep) thì phần còn thiếu được coi rằng là do các hoạt động của rễ cây rút ra từ lượng trữ tầng thấp theo tốc độ thực tế Ea. Ea tương ứng với lượng bốc hơi tiềm năng và biến đổi tuyến tính theo quan hệ lượng trữ ẩm trong đất, L/L max, của lượng trữ ẩm tầng thấp.
![](/Portals/0/2020-huy/CT_1.png)
Dòng chảy mặt:
Khi lượng trữ bề mặt đã tràn, U > Umax, thì lượng nước thừa PN sẽ gia nhập vào thành phần dòng chảy mặt. Thông số QOF đặc trưng cho phần nước thừa PN đóng góp vào dòng chảy mặt. Nó được giả thiết là tương ứng với PN và biến đổi tuyến tính theo quan hệ lượng trữ ẩm đất, L/Lmax, của lượng trữ ẩm tầng thấp.
![](/Portals/0/2020-huy/CT_2.png)
(2) trong đó: CQOF = hệ số dòng chảy tràn trên mặt đất (0 ≤ CQOF ≤ 1),
TOF = giá trị ngưỡng của dòng chảy tràn (0 ≤ TOF ≤ 1).
Phần lượng nước thừa PN không tham gia vào thành phần dòng chảy tràn sẽ thấm xuống lượng trữ tầng thấp. Một phần trong đó, ΔL, của nước có sẵn cho thấm, (PN-QOF), được giả thiết sẽ làm tăng lượng ẩm L trong lượng trữ ẩm tầng thấp. Lượng ẩm còn lại, G, được giả thiết sẽ thấm sâu hơn và gia nhập lại vào lượng trữ tầng ngầm.
Dòng chảy sát mặt:
Sự đóng góp của dòng chảy sát mặt, QIF, được giả thiết là tương ứng với U và biến đổi tuyến tính theo quan hệ lượng chứa ẩm của lượng trữ tầng thấp.
![](/Portals/0/2020-huy/CT_3.png)
(3) trong đó CKIF là hằng số thời gian dòng chảy sát mặt và TIF là giá trị ngưỡng tầng rễ cây của dòng sát mặt (0 ≤ TIF ≤ 1).
Diễn toán dòng chảy mặt và dòng sát mặt:
Dòng sát mặt được diễn toán qua chuỗi hai hồ chứa tuyến tính với cùng một hằng số thời gian CK12. Diễn toán dòng chảy mặt cũng dựa trên khái niệm hồ chứa tuyến tính nhưng với hằng số thời gian có thể biến đổi.
![](/Portals/0/2020-huy/CT_4.png)
(4) trong đó OF là dòng chảy tràn (mm/hr) OFmin là giới hạn trên của diễn toán tuyến tính (= 0,4 mm/giờ), và β = 0,4. Hằng số β= 0,4 tương ứng với việc sử dụng công thức Manning để mô phỏng dòng chảy mặt.
Theo phương trình trên, diễn toán dòng chảy mặt được tính bằng phương pháp sóng động học, và dòng chảy sát mặt được tính theo mô hình NAM như dòng chảy mặt (trong lưu vực không có thành phần dòng chảy mặt) được diễn toán như một hồ chứa tuyến tính.
Lượng gia nhập nước ngầm:
Tổng lượng nước thấm G gia nhập vào lượng trữ nước ngầm phụ thuộc vào độ ẩm chứa trong đất trong tầng rễ cây.
![](/Portals/0/2020-huy/CT_5.png)
(5) trong đó TG là giá trị ngưỡng tầng rễ cây đối với lượng gia nhập nước ngầm (0 ≤ TG ≤ 1).
Độ ẩm chứa trong đất:
Lượng trữ tầng thấp biểu thị lượng nước chứa trong tầng rễ cây. Sau khi phân chia mưa giữa dòng chảy mặt và dòng thấm xuống tầng ngầm, lượng nước mưa còn lại sẽ đóng góp vào lượng chứa ẩm (L) trong lượng trữ tầng thấp một lượng ΔL.
ΔL = PN-QOF-G (6)
Dòng chảy cơ bản
Dòng chảy cơ bản BF từ lượng trữ tầng ngầm được tính toán như dòng chảy ra từ một hồ chứa tuyến tính với hằng số thời gian CKBF.
Các thông số cơ bản của mô hình NAM gồm:
+ CQOF: hệ số dòng chảy tràn không có thứ nguyên, có phạm vi biến đổi từ 0,0 đến 0,9. Nó phản ánh điều kiện thấm và cấp nước ngầm. Vì vậy nó ảnh hưởng nhiều đến tổng lượng dòng chảy và đoạn cuối của đường rút. Thông số này rất quan trọng vì nó quyết định phần nước dư thừa để tạo thành dòng chảy tràn và lượng nước thấm. Các lưu vực có địa hình bằng phẳng, cấu tạo bởi cát thô thì giá trị CQOF tương đối nhỏ, ở những lưu vực mà tính thấm nước của thổ nhưỡng kém như sét, đá tảng thì giá trị của nó sẽ rất lớn.
+ CQIF: hệ số dòng chảy sát mặt, có thứ nguyên là thời gian (giờ)-1. Nó chính là phần của lượng nước trong bể chứa mặt (U) chảy sinh ra dòng chảy sát mặt trong một đơn vị thời gian. Thông số này ảnh hưởng không lớn đến tổng lượng lũ, đường rút nước.
+ CBL: là thông số dòng chảy ngầm, được dùng để chia dòng chảy ngầm ra làm hai thành phần: BFU và BFL. Trường hợp dòng chảy ngầm không quan trọng thì có thể chỉ dùng một trong 2 bể chứa nước ngầm, khi đó chỉ cần CBFL=0 - tức là lượng cấp nước ngầm đều đi vào bể chứa ngầm tầng trên.
+ CLOF, CLIF: các ngưỡng dưới của các bể chứa để sinh dòng chảy tràn, dòng chảy sát mặt và dòng chảy ngầm, không có thứ nguyên và có giá trị nhỏ hơn 1. Chúng có liên quan đến độ ẩm trong đất. Khi các giá trị của ngưỡng này nhỏ hơn L/Lmax thì sẽ không có dòng chảy tràn, dòng chảy sát mặt và dòng chảy ngầm. Về ý nghĩa vật lý, các thông số này phản ánh mức độ biến đổi trong không gian của các đặc trưng lưu vực sông. Do vậy, giá trị các ngưỡng của lưu vực nhỏ thường lớn so với lưu vực lớn.
+ Umax, Lmax: thông số khả năng chứa tối đa của các bể chứa tầng trên và tầng dưới. Do vậy, Umax và Lmax chính là lượng tổn thất ban đầu lớn nhất, phụ thuộc và điều kiện mặt đệm của lưu vực. Một đặc điểm của mô hình là lượng chứa Umax phải nằm trong sức chứa tối đa trước khi có lượng mưa vượt quá, PN xuất hiện, tức là U< Umax. Do đó trong thời kỳ khô hạn, tổn thất của lượng mưa trước khi có dòng chảy tràn xuất hiện có thể được lấy làm Umax ban đầu.
+ CK1,2, CKBF: là các hằng số thời gian về thời gian tập trung nước. Chúng là các thông số rất quan trọng, ảnh hưởng đến dạng đường quá trình và đỉnh.
Điều kiện ban đầu của mô hình:
1. U là lượng nước chứa trong bể chứa mặt (mm),
2. L là lượng nước chứa trong bể chứa tầng dưới (mm),
3. QOF - cường suất dòng chảy mặt sau khi diễn toán qua bể chứa tuyến tính (mm/h),
4. QIF - cường suất dòng chảy sát mặt khi qua bể chứa tuyến tính (mm/h),
5. BF - cường suất dòng chảy ngầm (mm/h).
Hiệu chỉnh thông số mô hình:
Hiệu chỉnh thông số mô hình nhằm xác định các thông số của mô hình để cho đường quá trình tính toán phù hợp nhất với đường quá trình thực đo. Việc hiệu chỉnh các thông số mô hình có thể được tiến hành bằng 2 phương pháp: phương pháp thử sai hoặc phương pháp tối ưu.