3. M ụ c tiêu, n ộ i dung và ph ươ ng pháp nghiên c ứ u
3.3.2 C ấ u trúc mô hình FRASC
Đất, hay bất cứ môi trường xốp nào, có khả năng giữ một lượng nước nào đó trong bể chứa chống lại lực trọng lực. Điều này đôi khi được gọi “khả năng độ ẩm đồng ruộng”. Bằng định nghĩa, nước được trữ trong các bể chứa này không thể trở thành dòng chảy và bể chứa có thể bị rút bởi bốc thoát hơi. Do đó bốc thoát hơi nước trở
thành thông sốđiều khiển gây nên sự thiếu hụt độẩm của đất.
Trong những vùng ẩm ướt, các khả năng đồng khảo sát độ ẩm của đất chỉ ra rằng lượng ẩm của đất có thể thường đạt đồng ruộng trong phạm vi toàn bộ vỏđất. Điều này ám chỉ rằng sự bổ sung của lượng mưa trong thời kỳ kế tiếp sẽ bằng lượng bốc thoát hơi trong thời kỳ xảy ra đó. Điều này đề xuất một mối quan hệđơn giản giữa mưa và dòng chảy rằng trước khi sinh ra dòng chảy lượng mưa phải thỏa mãn độ
thiếu hụt bên dưới khả năng đồng ruộng mà lần lượt bị gây ra bởi bốc thoát hơi nước khi độ ẩm đồng ruộng xảy ra lần cuối và sau khi trận mưa tiếp theo đó sẽ chảy
đi hết vì đất không còn khả năng giữ chúng được nữa. Đây là khái niệm hình thành dòng chảy trong trạng thái đầy đủ của bể chứa và là khái niệm cơ bản của mô hình Xinanjiang.
Lưu vực sông được chia thành các tiểu lưu vực. Dòng chảy ra từ các tiểu lưu vực
đầu tiên được mô phỏng và rồi sẽ diễn toán xuống các sông suối tới cửa chính của lưu vực. Dựa vào khái niệm hình thành dòng chảy trong trạng thái đầy đủ của bể
chứa, mô phỏng dòng chảy từ mỗi tiểu lưu vực bao gồm bốn phần chính sau:
- Dòng chảy sinh ra theo lượng mưa và sự thiếu hụt lượng ẩm trong đất; - Dòng chảy được phân thành ba thành phần: dòng mặt, dòng sát mặt và dòng ngầm; - Diễn toán dòng chảy từ dòng chảy địa phương tới cửa ra của mỗi tiểu lưu vực tạo nên dòng chảy ra của tiểu lưu vực đó. Hình 3.2 Sơđồ cấu trúc mô hình FRASC[27][29]
Sơđồ mô hình FRASC được chỉ ra trong Hình 3.2. Tất cả các ký hiệu bên trong hộp là các biến bao gồm đầu vào, đầu ra, biến trạng thái và biến nội bộ, trong khi các ký hiệu bên ngoài các hộp là các thông số.
Đầu vào mô hình là dữ liệu GIS từ DEM, thảm phủ thực vật, mưa trung bình lưu vực (P) và bốc hơi nước tiềm năng theo phương pháp Penman – Monteith (EM).
Đầu ra là lưu lượng dòng chảy (TQ) từ toàn bộ lưu vực và từng vị trí ô lưới trong mô hình DEM.
Các biến trạng thái là bể chứa nước ứng suất trung bình lưu vực (W), và bể chứa nước tự do trung bình lưu vực (S). Nước ứng suất trung bình lưu vực (W) có ba thành phần WU, WL và WD tương ứng ở tầng trên, tầng dưới và tầng sâu. FR là thông số diện tích góp phần dòng chảy có quan hệ với W. Phần còn lại của các kí hiệu bên trong hộp là tất cả các biến nội bộ. RB là dòng chảy trực tiếp từ phần nhỏ
thành phần RS, RI, và RG tương ứng đó là dòng chảy mặt, sát mặt và dòng ngầm. Ba thành phần dòng chảy này sẽđược truyền và chuyển hóa thành QS, QI và QG và kết hợp với nhau tạo thành tổng dòng chảy vào hệ thống sông suối của tiểu lưu vực. Dòng chảy cửa ra của lưu vực là Q.
Mô hình bao gồm 15 thông sốđược diễn tả bên ngoài các hộp trong Hình 3.2. K là hệ số giữa lượng bốc thoát hơi tiềm năng và bốc hơi chậu nếu tài liệu lượng bốc hơi chậu thực đo được sử dụng trong mô hình. WM và B là hai thông số mô tả sự phân bố nước ứng suất. WM là sức chứa nước ứng suất trung bình lưu vực có các thành phần UM, LM và DM. B là hệ số mũ của đường cong phân bố sức chứa nước ứng suất. IM là hệ số diện tích không thấm. SM và EX tương tự như WM và B trong khi chúng diễn tả phân bố sức chức nước tự do. KI và KG là các hệ số liên quan đến RI và RG. CI, CG, L, CS, KE và XE là các thông số cho việc diễn toán dòng chảy truyền. Ý nghĩa của tất cả các thông số sẽ được thảo luận thêm ở những phần tiếp theo.
Bảng 3.2 Các thông số của mô hình FRASC
Category Parameter Description Các thông số
bốc thoát hơi nước
K [] Hệ số giữa lượng bốc thoát hơi tiềm năng và bốc hơi chậu
UM [mm] Sức chứa nước ứng suất trung bình lưu vực tầng trên
LM [mm] Sức chứa nước ứng suất trung bình lưu vực tầng dưới
C [] Hệ số bốc hơi tầng sâu Các thông số
hình thành dòng chảy
WM [mm] Sức chứa nước ứng suất trung bình lưu vực
B [] Hệ số mũ của đường cong phân bố sức chứa nước ứng suất
IM [] Hệ số diện tích không thấm Các thông số
phân bố dòng chảy
SM [mm] Sức chứa nước tự do trung bình lưu vực
EX [] Hệ số mũ của đường cong phân bố sức chứa nước tự do
KG [] Hệ số dòng ra của nước tự do thành dòng chảy ngầm KI [] Hệ số dòng ra của nước tự do thành dòng chảy sát mặt Các thông số tập trung dòng chảy CG [] Hằng số sút giảm của lượng trữ nước ngầm CI [] Hằng số sút giảm của lượng trữ nước sát mặt CS []
Hằng số sút giảm trong phương pháp “trễ và lộ trình” cho diễn toán dòng chảy trên hệ thống sông suối trong phạm vi mỗi tiểu lưu vực