- Chức năng của mơ hình: Mô phỏng ảnh hưởng của việc quản lý sử dụng đất trong hệ thống lưu vực sông đến nguồn nước, bùn cát, hàm lượng chất hữu cơ trên
Q trình hình thành dịng chảy mặt
3.1.3.1. Tiến hành chạy mơ hình
Hình 3-5: Giao diện của mơ hình SWAT
- Mơ hình SWAT với cơng cụ tự động khoanh lưu vực (Watershed Delineation) sẽ giới hạn lưu vực cần tính tốn và mạng lưới sơng suối trên lưu vực đó sau khi vào dữ liệu bản đồ địa hình dưới dạng DEM, mơ hình sẽ tính tốn các thơng số thứ tự như độ dốc, chiều rộng, diện tích… của các lưu vực con trong lưu vực tính tốn (xem hình 3-6).
- Sau khi phân lưu vực thành các lưu vực con tiến hành chồng ghép bản đồ sử dụng đất và bản đồ phân loại đất để tính tốn % từng diện tích sử dụng đất và % từng loại đất trên các lưu vực con đó.
Sau đó tiến hành chạy thử nghiệm mơ hình với thời gian bắt đầu 1/1/1972, thời gian kết thúc mô phỏng là 31/12/1985, thời đoạn tính tốn theo tháng, với kết quả có thể xem dưới dạng biểu đồ Hình 3.9.
ình 3-7: Bản đồ mạng lưới sơng và các lưu vực con
Hình 3-9: kết quả dịng chảy của từng lưu vực con
3.1.3.2.Lựa chọn và hiệu chỉnh bộ thơng số cho mơ hình a. Các thơng số của mơ hình
* Các thơng số thuộc modun lưu vực
(1) Các thơng số tính lượng mưa hiệu quả
Phương pháp đường cong SCS (1972) thấm tính lượng mưa hiệu quả
CN2 : Chỉ số CN ứng với điều kiện ẩm II (trong file *.mgt)
Phương pháp thấm Green & Ampt tính tổng lượng dịng chảy
SOL_K : Ksat : Độ dẫn thuỷ lực ở trường hợp bão hòa (trong file *.sol) CN2 : Chỉ số CN ứng với điều kiện ẩm II (trong file *.mgt)
SOL_BD : ρb : Mật độ khối của đất (mg/m3) (trong file *.sol) CLAY : mc : % đất sét (trong file *.sol)
SAND : ms : % đất cát (trong file *.sol)
(2) Các thơng số tính lưu lượng đỉnh lũ
CH_N(1) : n : Hệ số nhám kênh dẫn (trong file *.sub)
(3) Các thơng số tính hệ số trễ dòng chảy mặt
SURLAG : surlag : hệ số trễ dòng chảy mặt (trong file *.bsn)
(4) Các thơng số tính các tổn thất - Tổn thất dọc đường
CH_K(1) : KRchR: Độ dẫn thuỷ lực của kênh dẫn (trong file *.sub)
- Tổn thất do bốc hơi
CANMX : canRmxR: Lượng trữ lớn nhất của vòm cây (trong file *.hru) ESCO : esco : Hệ số bốc hơi của đất (trong file *.sub)
(5) Các thơng số tính tốn dịng chảy ngầm
GWQMN : aqRshthr,qR: Ngưỡng sinh dòng chảy ngầm (mm) (trong file *.gw) ALPHA_BF : αRgwR: Hệ số triết giảm dòng chảy ngầm (trong file *.gw) REVAPMN : aqRsthr,rvpR : Ngưỡng sinh dòng thấm xuống tầng ngậm nước sâu (mm) (trong file *.gw)
(6) Các thơng số tính xói mịn và chuyển tải các chất dinh dưỡng
USLE_K: KRUSLER: Hệ số xói mịn đất trong phương trình USLE (trong file *.sol) USLE_C: CRUSLER : Hệ số ảnh hưởng của độ che phủ thực vật đến xói mịn đất trong phương trình USLE (trong file *.sol)
USLE_P: PRUSLER : Hệ số bảo vệ đất trước xói mịn đất trong phương trình USLE (trong file *.sol)
SURLAG: surlag : Hệ số trễ dòng chảy mặt (trong file *.bsn)
LAT_SED: concRsedR: Lượng bùn cát trong dòng chảy ngầm (trong file *.hru) NPERCO: βRNO3R: Hệ số thấm của Nitơ (trong file *.bsn)
PHOSKD: kd,surf: Hệ số phân chia hàm lượng phốt pho trong đất (trong file*.bsn) LAT_TTIME : Thời gian chảy truyền (ngày) (trong file *.hru)
SOL_CBN : orgCly : Hàm lượng Các bon hữu cơ trong đất (trong file *.sol) PERCOP : Hệ số thấm của thuốc trừ sâu (trong file *.bsn)
(7) Các thơng số tính chất lượng nước: SOL_CBN: orgCly : Phần trăm hàm lượng các bon trong 10 mm lớp đất đâu (%) (trong file *.sol)
(8) Các thơng số tính lớp phủ thực vật, cây trồng: T_BASE: Nhiệt độ tối thiểu cho
sự phát triển của cây trồng (P o
P
C) (trong file crop.dat)
* Các thông số modun diễn toán trong kênh:
(1) Các thơng số diễn tốn dịng chảy trong kênh
CH_N(2): n: Hệ số nhám của kênh chính (trong file *.rte)
MSK_X: X: Hệ số trọng số trong phương pháp Muskingum (trong file *.bsn) MSK_CO1: coefR1R: Hệ số C1 trong phương pháp Muskingum (trong file *.bsn) MSK_CO2: coefR2R: Hệ số C2 trong phương pháp Muskingum (trong file *.bsn) CH_K(2): Độ dẫn thuỷ lực của kênh chính (mm/giờ) (trong file *.bsn)
EVRCH:coefev : Hệ số hiệu chỉnh bốc hơi của kênh chính (trong file *.bsn) GW_REVAP : βrev : Hệ số Revap (trong file *.gw)
ALPHA_BNK : Hệ số tỷ lệ bờ kênh (trong file *.rte)
(2) Các thơng số diễn tốn dịng chảy bùn cát trong kênh
SPCON: cRspR: Hệ số trong phương trình chuyển tải bùn cát (trong file *.bsn) SPEXP: spexp : Hệ số trong phương trình chuyển tải bùn cát (trong file *.bsn) CH_COV: CRCHR : Hệ số che phủ kênh (trong file *.rte)
(3) Các thơng số diễn tốn các chất dinh dưỡng và thuốc bảo vệ thực vật
NPERCO: Hệ số tràn thấm của Nitơ (trong file *.bsn) PPERCO: Hệ số tràn thấm của Phốt Pho (trong file *.bsn)
PHOSKD: Hệ số chia lượng Phốt Pho trong đất (trong file *.bsn)
CHPST_KOC: Kd : Hệ số gia nhập chất bảo vệ thực vật (trong file *.swq)