GII THIU V SWAT
Soil and Water Assessment Tool Theoretical documentation, version 2005
by
S.L. Neitsch, J.G. Arnold, J.R. Kiniry, J.R. Williams.
2005 Grassland, Soil and Water Research Laboratory, Agricultural Research
Service, and Blackland Research Center, Texas Aagricultural Experiment Station,
Texas.U.S.A
Ngi dch:
Nguyn Duy Bình
i hc Nông nghip Hà ni
Hà ni, Tháng 8 nm 2009
MC LC
1. Quá trình phát trin ca SWAT 1
2. Gii thiu tng quan v SWAT 4
2.1 Vòng tun hoàn thy vn trên đt lin 6
2.1.1 Khí hu 7
2.1.2. Thy vn 9
2.1.3 Lp ph thc vt / Sinh trng ca cây ci 12
2.1.4 Xói mòn đt 13
2.1.5 Các cht dinh dng 13
2.1.6 Thuc tr sâu 15
2.1.7 Ch đ canh tác 16
2.2 Tính truyn dòng chy và cht trong kênh 16
2.2.1 Tính truyn dòng chy nc và cht trong kênh chính 17
2.2.2 Tính truyn dòng chy và cht trong h cha 18
3. Tài liu tham kho 20
DANH MC HÌNH 21
SWAT-Intro_vie.docx 1
GII THIU V SWAT
SWAT hay là Soil and Water Assessment Tool, mô hình dùng cho din tích mt lu
vc đc phát trin bi TS. Jeff Arnold theo đn đt hàng ca Trung tâm Nghiên cu
Nông nghip (Agricultural Research Service, ARS, thuc B Nông nghip Hoa k).
SWAT đc phát trin nhm mc đích d báo tác đng ca nhng hot đng qun lý và
khai thác đt đai lên ch đ dòng chy, lên tài nguyên nc, ch đ cht l lng và d
báo tàn d ca vic áp dng các cht hóa hc trong nông nghip
đi vi nhng lu vc
sông ln và phc tp vi đc đim th nhng, s dng và qun lý đt đai trong điu kin
bin đng trong mt thi gian dài. có th đt đc nhng yêu cu đó, mô hình SWAT
phi tha mãn các điu kin sau.
Mô hình phi đc phát trin da trên c s các quá trình vt lý. Thay bng vic s
dng nhng công thc tng quan din t quan h gia d liu đu vào và đu ra, SWAT
cn phi có đc thông tin v khí hu, tính cht th nhng, d liu đa hình, thm thc
vt, và các hình thái canh tác s dng đt đai trong lu vc cn mô phng. SWAT s
dng nhng thông tin đu vào đó đ mô phng nhng quá trình vt lý liên quan đn
chuyn vn ca nc, c
a cht bùn cát l lng trong nc (phù sa), phát trin ca cây
trng, vòng tun hoàn ca các cht dinh dng, v.v.
SWAT s dng nhng d liu đu vào đã có sn. Trong khi SWAT có th đc s
dng đ nghiên cu nhng vn đ đc thù nh chuyn vn ca vi khun, thì đi vi
nhng vn đ dòng chy, nó vn có th vn hành đc vi nhng d li
u đã có sn trong
các c quan ca chính ph.
Các phng pháp tính toán trong SWAT cn đt đc tc đ ti u. Vic mô
phng nhng lu vc có din tích ln vi nhng chin lc qun lý phc tp phi thc
hin đc mà không cn đu t tài chính ln và thi gian dài.
SWAT phi cho phép đánh giá tác đng lâu dài. Rt nhiu vn đ cn gii quyt
hin nay bao gm quá trình tích ly chm rãi nhng cht ô nhim vi tác đng lên khi
nc dòng h lu. có th nghiên cu nhng vn đ nh vy, ngi ta cn nhng kt
qu mô phng kéo dài vài chc nm.
SWAT là mt mô hình mô phng theo chui thi gian liên tc, loi mô hình cho kt
qu thi gian dài. Mô hình không đc thit k mô phng cho nhng quá trình xy ra
trong mt thi gian ng
n nh quá trình truyn l do mt trn ma.
HP: u đim ca phng pháp tip cn này là
Có th mô phng cho nhng lu vc mà đó cha có d liu quan trc thc t.
Tác đng tng đi ca nhng b d liu đu vào khác nhau (ví d nh thay đi v
phng thc qun lý tài nguyên, bin đng khí hu, thm th
c vt, v.v.).
1. Quá trình phát trin ca SWAT
SWAT kt hp các đc trng ca nhiu mô hình ca ARS và là kt qu tha k trc
tip ca mô hình SWRRB
1
(Simulator for Water Resources in Rural Basins) (Williams
et al., 1985; Arnold et al., 1990). Nhng mô hình đã có đóng góp đáng k vào s phát
1
SWRRB là mô hình đng thái đc phát trin đ mô phng tác đng ca các ngun phát tán din
rng.
Gii thiu v SWAT (Soil and Water Assessment Tools)
trin ca SWAT bao gm: CREAMS
2
( Chemicals, Runoff, và Erosion from Agricutural
Management System) (Knisel, 1980), GLEAMS
3
(Groundwater Loading Effects on
Agicultural Management Systems) ) (Leonard et al., 1987), và EPIC
4
( Erosion –
Productivity Impact Calculator) (Williams et al., 1984).
Quá trình phát trin ca SWRRB bt đu vi vic sa đi mô hình thy vn v lng
ma ngày ca mô hình CREAMS. Các thay đi chính so vi mô hình thu vn CREAMS
bao gm: a) mô hình đã đc m rng đ cho phép đng thi tính toán trên nhiu lu vc
con đ d đoán dòng chy t nc ma; b) b sung thêm mô hình v nc ngm, hay mô
hình v dòng chy hi lu; c) b sung mô dun h cha nh
m tính toán tác đng ca ao h
nông tri và h cha nc đn ch đ dòng chy và lu lng bùn lng; d) thêm mô hình
mô phng thi tit cha đng các d liu cho lng ma, bc x mt tri, và nhit đ
nhm to điu kin thun li cho các mô phng lâu dài và cung cp bin đng thi tit
theo thi gian và không gian; e) phng pháp d báo giá tr cc đ
i ca tc đ dòng chy
mt đt đã đc ci thin; f) mô hình phát trin cây trng EPIC đã đc b sung nhm
mô phng tác đng ca phát trin cây trng hàng nm; g) b sung mô dun tính truyn l;
h) b sung mô đun vn chuyn bùn cát nhm mô phng chuyn vn ca phù sa bùn cát
qua các ao, h cha, dòng sông sui và thung lng; và i) các tính toán v tn tht truyn
dn.
Mô hình đc chú trng s dng trong cu
i nhng nm 1980 ch yu nhm đánh giá
cht lng nc và s phát trin ca SWRRB đã chng t điu này. Nhng ci tin đáng
k ca mô hình SWRRB vào thi gian đó bao gm s kt hp vi: a) thành phn mô
phng s phn ca thuc tr sâu trong mô hình GLEAMS; b) phng pháp SCS đ tính
toán giá tr cc đi ca tc đ dòng chy mt đ
t; và c) phng trình mi đuc phát trin
v bi lng. Nhng sa đi, b sung này đã tng cng kh nng s dng ca mô hình
trong vic gii quyt nhiu vn đ v qun lý lu vc.
Vào cui nhng nm 1980, Cc các vn đ v ngi da đ (the Bureau of Indian
Affairs) cn mt mô hình đ c lng tác đng vào dòng chy h lu c
a công tác qun
lý ngun nc trong phm vi lu vc khu đt giành cho ngi da đ Arizona và New
Mê-hi-cô. Trong khi SWRRB đã đc s dng mt cách d dàng cho các sông có din
tích lên đn vài trm km vuông, Cc các vn đ v ngi da đ mun mô phng dòng
chy cho lu vc rng hàng nghìn km vuông. i vi khu vc rng ln nh vy, lu vc
đc nghiên cu cn thit phi b chia thành vài trm l
u vc con. SWRRB ch cho phép
chia lu vc thành mi lu vc con và mô hình tính truyn vn chuyn dòng nc và
bùn lng ra khi lu vc con trc tip đn đim ra ca lu vc dòng sông. Nhng hn ch
này đã dn đn s phát trin ca mt mô hình có tên gi là ROTO (Routing output to
outlet) (Arnold et al., 1995), trong đó kt qu t nhiu ln chy mô hình SWRRB cho các
2
Tuân th lut Bo v ngun nc, vào đu nhng nm 1970 c quan ARS thành lp mt nhóm gm
các nhà khoa hc lien ngành đn t khp nc M, nhm phát trin mt mô hình mô phng tác đng ca
các ngun phát tán din rng da trên nguyên lý phân tích tng quá trình thành phn. Kt qu là mô hình
CREAMS đã đc to nên. CREAMS là mt mô hình dành cho din tích khu đng rung, đc thit k đ
mô phng tác đng ca các hot đng làm đt và qun lý đt đai đn dòng nc, lng cht l lng và bùn
lng, cht dinh dng và lng d tha thuc tr sâu thoát ra khi khu đng rung. Mt s nhng mô hình
khác do ARS đã phát trin nh GLEAMS, EPIC, SWRRB và AGNPS đu có ngun gc t mô hình
CREAMS
3
GLEAMS là mô hình mô phng tác đng ca các ngun phát tán din rng nhng chú trng tác đng
ca thuc tr sâu và cht dinh dng lên ngun nc ngm.
4
EPIC vn đc phát trin đ mô phng tác đng ca xói mòn đt lên nng sut cây trng và hin nay
đã phát trin thành mt mô hình toàn din v qun lý nông nghip trên din tích đng rung, mô phng tác
đng ca các ngun phát tán din rng.
Gii thiu v SWAT (Soil and Water Assessment Tools)
lu vc con đc chuyn theo dòng chy trong các kênh và h cha. ROTO cung cp
mt phng pháp tip cn tính toán theo tng đon sông và khc phc đc nhc đim
ca SWRRB v gii hn s lu vc con bng cách "liên kt" nhiu ln chy mô hình
SWRRB li vi nhau. Mc dù phng pháp tip cn này rt hiu qu, nhng d liu đu
vào và d liu
đu ra ca nhiu ln chy SWRRB tr nên cng knh và cn kh nng lu
tr ln trong máy tính. Ngoài ra, tt c các ln chy mô hình SWRRB phi đc thc
hin đc lp và sau đó nhp kt qu vào mô hình ROTO đ thc hin bc tính truyn
theo các kênh và h cha. khc phc đc nhng rc ri này, mô hình SWRRB và
Roto đã đc kt hp thành mt mô hình duy nht, có tên g
i là SWAT. Trong khi
SWAT cho phép mô phng khu vc rt rng ln, nó gi li tt c các tính nng đã làm
cho mô hình SWRRB có giá tr nh mt mô hình mô phng.
T khi SWAT đc to ra vào đu nhng nm 1990, nó đã liên tc tri qua nhiu ln
đc xem xét, đánh giá và ci tin nhm m rng kh nng mô phng. Nhng ci tin
đáng k nht ca các mô hình theo các phiên bn khác nhau bao gm:
SWAT94.2: B sung khái nim đn v đng nht v phn ng thu vn (HRUs).
SWAT96.2: Phng án t đng bón phân và ti nc đc thêm vào nh là nhng
qun lý tùy chn; tính toán lng nc do tán lá cây lu tr; thành phn mô phng CO
2
trong mô hình tng trng cây trng phc v các nghiên cu v bin đi khí hu; b sung
phng trình Penman-Monteith v bc thoát nc tim nng; dòng chy theo chiu ngang
trong đt da trên mô hình lu tr đng thái; b sung phng trình cht lng nc v
thành phn dinh dng ca dòng chy t mô hình QUAL2E; tính truyn vn chuyn
thuc tr sâu trong dòng chy sông sui.
SWAT98.1: Ci tin chng trình con v mô phng l
ng tuyt tan; ci thin tính
toán cht lng nc trong dòng sông sui; m rng tính truyn vòng tun hoàn cht dinh
dng; tác đng chn th đng c, tác đng cách thc áp dng phân bón, và thêm phng
án tiêu nc s dng cày sâu nh là mt phng thc qun lý, sa đi mô hình đ có th
áp dng khu vc nam bán cu.
SWAT99.2: Ci tin tính truyn vòng tun hoàn cht dinh dng, ci tin tính toán
ru
ng lúa/đm ly, b sung phn c tính lng tn tht cht dinh dng do quá trình bi
lng trong h cha / ao / đm ly; b sung lng nc cha đc do b sông, b sung
tính truyn kim loi theo th t các đon sông sui; tt c các nm tài liu tham kho
trong mô hình đã thay đi biu th t 2 ch s thành 4 ch s, b sung phng trình nh
hng các khu
đô th lên dòng chy t mô hình SWMM theo phng trình quan h ca
C quan Thm dò a cht (USGS).
SWAT2000: B sung tính truyn vn chuyn vi khun trong dòng chy; b sung
phng trình thm Green & Ampt, ci thin mô hình mô phng thi tit, cho phép đc
vào hoc mô phng d liu bc x mt tri hàng ngày, đ m tng đi, và tc đ gió;
cho phép đc vào hoc c tính các giá tr bc thoát nc tim nng ET cho l
u vc; xem
xét li tt c các phng pháp c tính ET tim nng; ci thin quá trình liên quan đn đ
cao ca b sông; cho phép mô phng không gii hn s lng h cha; b sung phng
pháp tính truyn Muskingum; sa đi tính toán mô phng trng thái ngng hot đng
sng (ng đông) cho phù hp vi các khu vc nhit đi.
SWAT2005: Ci thin tính truyn vn chuyn vi khun trong dòng chy; thêm kch
bn d báo thi tit; b sung phn mô phng lng ma ri; thông s lu tr nc trong
tính toán giá tr CN hàng ngày có th là là hàm s ca lng nc trong đt (đ m đt)
hay ca lng bc thoát hi nc t cây ci.
Gii thiu v SWAT (Soil and Water Assessment Tools)
Ngoài nhng thay đi đã đc lit kê trên, giao din cho các mô hình đã đc phát
trin cho môi trng h điu hành Windows (Visual Basic), GRASS, và ArcView. Mô
hình SWAT cng đã đc kim chng cht ch.
2. Gii thiu tng quan v SWAT
SWAT cho phép mô phng mt s các quá trình vt lý khác nhau xy ra trong mt
lu vc sông. Nhng quá trình này s đc xem xét và tng kt li trong phn này ca tài
liu. Bn đc cn phi tham kho thêm các mc chuyên đ nu mun tìm hiu sâu hn v
tng quá trình.
Hình 1 Bn đ lu vc h Fork đông bc Texas vi phân b s dng đt và mng li
sông ngòi.
Vi m
c đích mô phng, mt lu vc có th đc chia ra thành nhiu lu vc con.
Vic s dng khái nim lu vc con đc bit tin li khi mt lu vc bao gm nhiu vùng
có đc tính s dng đt khác nhau hay có loi đt khác nhau tác đng đn ch đ thy vn
khác nhau. Khi chia lu vc thành các lu vc con, ngi s dng có th phân bit và ch
ra các vùng din tích khác nhau theo v
trí không gian. Hình 2 th hin kt qu chia nh
lu vc trong Hình 1.
Gii thiu v SWAT (Soil and Water Assessment Tools)
Hình 2 Chia nh lu vc h Fork.
D liu đu vào cho mi mt lu vc con đc phân loi hay t chc li theo nhng
nhóm sau : khí hu, đn v đng nht v tác đng thy vn hay còn gi là HRU; ao h /
vùng ngp nc; nc ngm; và dòng kênh chính, hoc các khúc sông, đng tiêu nc
ca lu vc con; đn v đng nht v tác đng thy vn (HRU) là vùng con trong mt l
u
vc có cùng mt loi s dng đt, loi th nhng, và loi kt hp các bin pháp qun lý.
Trong bt k ng dng mô hình nào vi SWAT thì cân bng nc là đng lc cho
mi quá trình xy ra trong lu vc. có th mô phng chính xác chuyn vn ca thuc
tr sâu, bùn lng hay cht dinh dng, thì vòng tun hoàn thy vn đc mô phng phi
th hin đc b
n tính ca quá trình đó đang thc s xy ra trong mt lu vc trong th
gii thc. Nhng mô phng ch đ thy vn trong lu vc có th đc phân thành hai
loi. Vic phân chia đu tiên là phn lc đa ca vòng tun hoàn thy vn nh đc biu
th trong Hình 3. Nhng quá trình thy vn trong vòng tun hoàn xy ra trên đt lin kim
soát lng dòng chy, bùn lng và lng thu
c tr sâu đi vào trong chuyn vn ca nc
trong dòng sông chính ca mi mt lu vc con. Loi phân chia th hai là quá trình lng
nc hay các cht trong nc trong vòng tun hoàn thy vn, có th đc đnh ngha nh
là vn chuyn ca nc, bùn cát, v.v. theo mng li dòng sông và cui cùng đu đi ra
qua đim ra ca lu vc.
Gii thiu v SWAT (Soil and Water Assessment Tools)
Hình 3 S đ vòng tun hoàn thy vn.
2.1 Vòng tun hoàn thy vn trên đt lin
Vòng tun hoàn thy vn trong phn mm SWAT đc thành lp da trên phng
trình sau:
trong đó SW
t
là đ m đt (mm H
2
O) cui cùng, SW
0
là đ m đt ngày đu tiên (mm
H
2
O), t là ch s ngày (ngày), R
day
là lng ma/tuyt ri trong ngày i (mm H
2
O), Q
surf
là
lng dòng chy tràn mt đt trong ngày i (mm H
2
O), E
a
là lng bc thoát nc trong
ngày i (mm H
2
O), w
seep
là lng nc thm vào tng không bão hòa trong t tng r cây
trong ngày i (mm H
2
O), and Q
gw
là lng dòng chy hi lu trong ngày i (mm H
2
O).
Vic phân chia lu vc nghiên cu thành các lu vc con cho phép mô hình th hin
đc s khác nhau v lng bc thoát nc đi vi các loi cây trng và loi đt khác
nhau. Dòng chy tràn trên mt đt (runoff) đc mô phng riêng cho tng đn v đng
nht thy vn (HRU) và tính truyn l đ thu đc tng dòng chy tràn mt đt cho toàn
b lu vc. iu này làm tng đ chính xác ca mô hình và bi
u th tt hn phng trình
cân bng nc v mt vt lý.
Gii thiu v SWAT (Soil and Water Assessment Tools)
Hình 4 Vòng tính toán cho HRU/Lu vc con.
Hình 4 Th hin trình t quá trình mà SWAT mô hình hóa vòng tun hoàn thy vn
trên lc đa đt lin. Nhng d liu đu vào và các quá trình liên quan đn vòng tun hoàn
này đc trình bày c th hn phn tip sau đây.
2.1.1 Khí hu
Khí hu ca mt lu vc cung cp ngun nc và nng lng cho vòng tun hoàn
nc và khng ch tm quan trng tng thành phn trong vòng tun hoàn th
y vn.
Nhng bin s cn thit cho đu vào ca SWAT bao gm d liu hàng ngày v lng
ma/tuyt, nhit đ không khí cc đi và cc tiu, bc x mt tri, tc đ gió và đ m
Gii thiu v SWAT (Soil and Water Assessment Tools)
tng đi ca không khí. Nhng d liu này có th là d liu thu đc t các trm khí
tng hoc d liu mô phng.
D liu khí tng mô phng. Phn mm mô phng giá tr trung bình ngày ca các
thông s khí tng t các giá tr trung bình tháng. Sau đó mi lu vc con s đc gán
mt b d liu khí tng đã đc to ra trong quá trình mô phng khí hu. Các giá tr cho
các l
u vc con này đc to ra mt cách đc lp và không h có quan h không gian nào
gia chúng vi nhau.
Mô phng d liu ma/tuyt. SWAT s dng mô hình ca Nicks (1974) đ mô phng
d liu ma cho nhng lu vc con không có d liu quan trc. Mô hình ma/tuyt này
cng đc s dng đ to nhng d liu còn thiu trong chui d liu quan trc. Mô hình
mô phng ma s dng chui Markov bc nht đ xác đnh mt ngày s có ma hay
không ma bng cách so sánh mt s đc to ra ngu nhiên (0,0-1,0) vi xác sut xy ra
ma hàng tháng mà ngi s dng phi đa vào. Nu ngày đó có ma thì lng ma s
đc to nên t phân b lch chun hay phân b hàm m ci tin.
Bin đng cng đ trn ma trong ngày. N
u cn thit phi có phân b lng ma
theo gi trong ngày thì phn mm s dng mt hàm s m kép nh là hàm phân b ca
cng đ ma trong mt trn ma. Khi đó cng đ ma s tng lên theo hàm s m
theo thi gian đn mt giá t cc đi, hay là đnh cng đ ma. Sau đó cng đ ma s
gim d
n cho đn thi đim kt thúc trn ma.
Mô phng d liu nhit đ không khí và bc x mt tri Nhit đ cc đi và cc
tiu, cng nh bc x mt tri đc to nên t phân b xác sut chun. Phn mm s
dng phng trình liên tc đ biu th cho tác đng ca thi tit ma n
ng đi vi bin
đng ca nhit đ và bc x mt tri. Nhit đ cc đi và bc x mt tri b gim xung
khi có ma và tng lên khi tri không ma. Nhng điu chnh đó đc thc hin sao cho
v lâu dài thì các giá tr trung bình tháng phi trùng hp vi giá tr ca d liu đu vào.
Mô phng d liu tc
đ gió. Phn mm s dng phng trình hàm m đ to ra d
liu tc đ gió ngày t d liu tc đ gió trung bình tháng.
Mô phng d liu đ m không khí tng đi. Mô hình đ m tng đi hàng ngày
s dng phân b tam giác t giá tr trung bình tháng. Ging nh đi vi nhit đ và bc
x mt tri, giá tr trung bình ngày ca đ
m tng đi cng đc điu chnh đ biu th
tác đng ca ngày có ma và không có ma.
Tuyt. SWAT xác đnh ma hay tuyt trên c s nhit đ không khí trung bình ngày
xy ra hin tng ma/tuyt.
Lp tuyt ph. Thành phn ca SWAT chuyên tính toán đ dày tuyt ri ph xung
mt đt đã đc ci tin t mt mô hình đn gin mô ph
ng đ dày tuyt ph đng đu
thành mt mô hình phc tp hn, cho phép tính toán bin đng ca đ dày tuyt ph do
bóng râm, trt tuyt, điu kin đa hình và thm thc vt gây ra. Ngi s dng phi đa
d liu đ dày ln nht cho phép đi vi toàn b din tích nghiên cu mà nu tuyt ri
dày hn thì mt lng tuyt s tràn ra các di
n tích khác. Khi tng lng tuyt ri nh
hn đ dày cho phép này thì đ dày tuyt che ph s gim theo hàm s phi tuyn trên c
s đng cong gim din tích.
Lng tuyt tan. Lng tuyt tan ph thuc vào nhit đ không khí, nhit đ ca
khi tuyt, tc đ tuyt tan, và din tích lp tuyt che ph. nhng din tích có tuyt
ph, tuyt s tan vào nhng ngày có nhi
t đ ln nht vt quá 0°C theo mt hàm s
tuyn tính quan h gia nhit đ khi tuyt trung bình và nhit đ ti hn gây tuyt tan.
Tuyt tan cng đc mô phng tng t nh đi vi lng mua, bao gm c dòng chy
Gii thiu v SWAT (Soil and Water Assessment Tools)
mt đt và lng nc thm xung đt. i vi tuyt tan, nng lng ca ht nc đc
coi là bng không và phn mm cn c vào đó đ tính toán giá tr đnh ca lu lng
dòng chy mt đt vi gi thuyt tuyt tan đu trong vòng 24 ting đng h.
Di phân cách theo đ cao. Mô hình cho phép giá tr đ cao trong mi lu vc con
có th đ
c phân loi thành 10 khong đ cao. dày tuyt ph và lng tuyt tan đc
mô phng cho riêng tng di giá tr đ cao. Nh vy, mô hình có kh nng đánh giá đc
s khác nhau v đ dày tuyt ph và lng tuyt tan do bin đng không gian v đa hình,
lng ma tuyt và nhit đ không khí.
Mô phng nhit đ đt. Nhit đ đt có tác đng lên dòng chuyn đng ca nc
trong đt và tc đ phân hy ca nhng phn d tha có trong đt. Nhit đ trung bình
ngày đc mô hình tính toán đi vi v trí mt đt và đim gia tâm mi tng đt. Nhit
đ b mt mt đt là hàm s ca lng tuyt ph, thm thc vt và đ dày lp ph d
tha, nhit đ lp đt không có che ph
, và nhit đ đt ngày trc đó. Nhit đ ca
mt lp đt là hàm s ca nhit đ mt đt, nhit đ không khí trung bình hàng nm và đ
sâu mà đó nhit đ không b nh hng bi bin đng nhit đ không khí trên mt đt.
sâu này, còn có tên gi là ‘the damping depth’, ph thuc vào trng lng riêng ca
đt và đ m đt.
2.1.2. Thy v
n
Mt khi ri gn xung mt đt, nc ma có th b tán lá cây gi li mt phn và
phn còn li s tip tc ri xung mt đt. Khi đn b mt đt, nc s thm vào trong đt
hoc chy tràn trên mt đt to nên dòng chy mt. Dòng chy mt di chuyn nhanh v
phía các rãnh kênh và là yu t chính to nên hiu ng dòng ch
y ngn hn. Lng nc
thm xung đt có th nm li trong các lp đt đ sau đó bin thành hi nc hay tip
tc chuyn đng theo nhng khe dn đn h thng các tng nc ngm. Hình 5 trình bày
s đ đng vn chuyn ca nc nh đc mô phng bi mô hình SWAT trong mt
HRU (đn v đng nht thy vn).
N
c lu tr trong tán lá cây. Nc lu tr trong tán lá cây là lng nc đc gi
li bi lp thm thc vt b mt trái đt và có kh nng b bc hi sau đó. Trong
phng pháp đng cong ca SCS (The US Soil Conservation Service) đ tính toán dòng
chy mt, lng nc lu tr trong tán lá cây cng đc tính đn trong tính toán dòng
chy tràn trên b mt đt. Tuy nhiên, nu nhng phng pháp nh ca Green và Ampt
đ
c s dng đ tính toán dòng thm thì lng nc lu tr trong tán lá cây phi có mô
hình mô phng riêng bit. SWAT cho phép ngi dùng nhp vào d liu v lng nc
ti đa có th đc tán lá cây lu tr khi ch s din tích lá đt giá tr ti đa. Giá tr này và
ch s din tích lá đc s dng bi các mô hình đ tính toán lng nc lu tr ti đa
bt c
thi đim nào trong chu k tng trng ca cây trng/ lp ph thc vt. Khi tính
toán lng nc bc hi, nc lu tr trong các tán lá cây là lng nc đu tiên b tác
đng.
Dòng thm vào đt. Dòng thm đ cp đn hin tng dòng nc chuyn đng
xuyên qua b mt đt đ thm vào trong lòng đt. Khi dòng thm này xy ra liên tc, đt
ngày càng tr nên m
t, thì tc đ thm s gim bt theo thi gian cho đn khi đt đn
mt giá tr n đnh. Tc đ thm ban đu ph thuc vào đ m ca đt khi nc bt đu
thm vào đt. Tc đ thm ti hn cui cùng có giá tr tng đng vi h s dn nc
thu lc ca đt bão hoà. Bi vì ph
ng pháp đng cong SCS đc s dng đ tính toán
dòng chy mt có bc thi gian tính toán bng mt ngày, nên nó không th trc tip tính
toán mô hình thm. Lng nc thm vào trong đt đc tính toán nh là khong chênh
gia
trc
t
nc
ma
tron
g
phân
ra kh
b
ão
h
mt
l
B
mà n
nc
trn
g
thc
v
ch s
WAT
lng ma
v
t
ip mô hìn
h
trong lòng
đ
â
g
.
M
b c
qua tng l
i
h
òa .
.
l
p
.
B
c l
g
, và
bng và
v
t mt các
h
din tích l
á
(Soil and Wat
e
v
à l
ng dò
n
h
c
nh hn.
. Phâ
n
đ
t sau khi
) b m
t
â
M
t khi
a SWAT
p đt
lu l
phân b
0
°
5
.
ng hoc r
n
.
; bc
h
tuyt
h
riêng bit
n
đc c
á
(din
e
r Assessment
n
g chy b
m
nhng nó đ
ò
n
b
t
đt
tro
n
phân
ng n
b
°
C hoc
nc là
m
n
ti b mt
t
h
i t l
á
.
M
n
h đc m
ô
tính nh là
m
tích
Tools)
m
t. Phn
g
ò
i hi
đ cp đ
.
s khá
c
n
g sut toàn
b
.
nh hng b
,
m
t thut ng
t
rái đt hoc
bao gm
á
cây (qua q
u
M
ô hình tín
h
ô
t bi Rit
c
m
t
lá
g
pháp thm
G
d l
i
n chuyn đ
g
c
bit v
b
n đ d bá
h
các lp
i nhit đ
đ
S
W
chung
gn b m
b
u
h
bc
h
c
hie (1972).
ca
H
G
reen và A
m
i
u lng m
ng tc
. T
h
o dòng chu
y
m sâu xu
n
di
đ
t. Nu nhi
W
AT
tt c
c
t trái đt, t
r
c hi t sô
n
) và t
h
i t mt đ
bc
h
H
RU
m
pt không
a
ca
(ln
g
h
ành phn
y
n
gx
y
ch
a
t đ trong
không
.
c
ác quá t
r
ìn
h
thành hi
n
g h, đt
t và
h
i nc
v
g
y
a
h
v
à
Gii thiu v SWAT (Soil and Water Assessment Tools)
ca HRU đó). Lng bc thoát nc thc t đc c tính s dng mt hàm s m ca
chiu sâu đt vi đ m đt. Lng nc mt đi qua quá trình hô hp ca cây đc mô
phng nh là mt hàm s tuyn tính gia lng bc thoát nc ln nht và ch s din
tích lá.
Lng bc thoát nc ln nht. B
c thoát nc ln nht là tc đ bc hi nc có
giá tr khi nó xy ra t mt khu vc rng ln và đng đu đc ph bng thm thc vt
đang phát trin vi ngun cung cp nc không gii hn t đt. Tc đ này đc gi đnh
là không b nh hng bi các quá trình vi khí hu nh chuyn đng khi hoc tích tr
nhit lng. Mô hình cung cp ba la chn đ c tính bc thoát nc: phng pháp
Hargreaves (Hargreaves et al, 1985.), Priestley-Taylor (Priestley và Taylor, 1972), và
Penman-Monteith (Monteith, 1965).
Dòng chy ngm sát mt đt. Dòng chy ngm theo chiu ngang sát mt đt, hay còn
gi là ‘interflow’, là dòng chy vào sông ngòi bt ngun t tng bên di sát b mt đt,
nhng trên vùng có tng đá bão hòa nc. Dòng chy ngm sát mt đt (đ sâu 0-2 m)
đc tính toán đng thi vi quá trình phân b đ m đt nh
đã đc din t trên. Mt
mô hình đng thái lu tr nc đc s dng đ d báo dòng chy bên trong mi lp đt.
Mô hình tính đn nhng bin đng giá tr ca h s dn nc, đ dc đa hình và đ m
đt.
Dòng chy mt. Dòng chy mt, hay là dòng chy tràn trên mt đt, là dòng chy xy
ra theo chiu b mt sn dc. Mô hình SWAT s dng lng ma ngày hay ma gi đ
mô phng khi lng dòng chy mt và giá tr đnh ca tc đ dòng chy tràn trên mt
đt cho tng din tích HRU.
Tng lng dòng chy mt đc tính toán da trên mt sa đi ca phng pháp
đng cong SCS (Cc bo v đt đai, SCS, thuc B Nông nghip Hoa k, USDA, 1972)
hay s dng ph
ng pháp thm Green & Ampt (Green và Ampt, 1911). Trong phng
pháp s đng cong SCS, giá tr ch s đng cong bin đi phi tuyn tính vi đ m đt.
Giá tr tr s đng cong gim xung khi đ m đt có giá tr gn bng đ m cây héo và
tng đn gn 100 khi đ m đt đt đc giá tr gn bng đ m bão hòa. Phng pháp
Green & Ampt đòi hi phi có đc d li
u lng ma phân b theo bc thi gian nh
hn mt ngày và tính toán lng nc thm theo mt hàm s gia th nng gi nc v
trí b mt gii hn di ca khi nc thm và h s dn nc thy lc thc t. Lng
nc ma còn li (không thm đc xung đt) tr thành lng nc gây dòng ch
y tràn
mt đt. SWAT bao gm c mt module c lng dòng chy t đt đông lnh khi nhit
đ trong lp đt gn sát mt đt có giá tr bé hn 0°C. Mô hình nh vy làm tng dòng
chy tràn trên đt đông lnh nhng vn cho phép có giá tr đáng k v dòng thm khi các
loi đt khô b đông lnh.
Lu lng đnh ca dòng chy mt đc d
báo vi vic áp dng mt sa đi ca
phng pháp biu đ thy vn hp lý. Nu tóm tt, thì phng pháp biu đ thy vn hp
lý có nn tng cho rng nu mt lng ma bt đu ngay lp tc vi cng đ i và tip
tc vi cng đ i đó vô thi hn, thì lu lng dòng chy s liên tc tng lên cho đn
khi thi gian ma bng thi gian tp trung nc, t
c
, khi mà lng ma ri trên tt c các
lu vc con cùng có đóng góp vào dòng chy ti đim đu ra ca lu vc. Trong công
thc hp lý đã sa đi, lu lng đnh ca dòng chy mt là mt hàm s t l thun gia
lng ma hàng ngày trong lu vc con, tng lng dòng chy tràn hàng ngày, và thi
gian tp trung nc ca lu vc con, t
c
. Phn lng ma xy ra trong lu vc con đc
c tính nh là mt hàm s ca tng lng ma hàng ngày bng cách s dng k thut
tính toán xác sut. Thi gian tp trung nc ca mt lu vc con đc c tính bng cách
Gii thiu v SWAT (Soil and Water Assessment Tools)
s dng công thc Manning's bao gm c vic xem xét c dòng chy tràn và dòng chy
trong kênh rch.
Ao vng nc. Ao nc là các cu trúc đa hình lu tr nc nm trong mt lu vc
con và ct ngang dòng chy tràn trên mt đt. Din tích lu vc ca ao đc coi là mt
phn nh trong tng din tích lu vc con. Ao đc gi đnh có v trí nm cách bit vi
các kênh chính trong mt lu v
c con và s không bao gi nhn đc nc t các lu
vc con phía thng lu. Lng nc lu tr trong ao h là mt hàm s ca th tích ao,
lng dòng chy hàng ngày vào ao và thoát khi ao, lng nc thm qua b và đáy ao
và lng bc hi t mt nc ao. Yêu cu đu vào cho mô hình là dung lng lu tr và
din tích b mt ca ao khi ao cha đy nc. Din tích b m
t ao khi ao không đy nc
đc c tính nh là mt hàm s phi tuyn tính ca th tích lu tr.
Sông nhánh. Trong mt lu vc con ngi ta phân ra hai loi kênh: kênh chính và
kênh nhánh. Kênh nhánh là mng li các kênh ta ra t kênh chính trong phm vi lu
vc con. Mi kênh nhánh trong phm vi mt lu vc con ch tp trung nc ma/tuyt t
mt phn ca lu vc con và không nhn đc s đóng góp ca nc ngm vào dòng
ch
y ca kênh. Tt c các dòng chy trong các kênh nhánh đu kt thúc trong dòng ca
kênh chính và đc tính truyn dc theo kênh chính ca lu vc con. SWAT s dng các
thuc tính ca các kênh nhánh đ xác đnh thi gian tp trung nc cho các lu vc con.
Tn hao nc trong kênh rãnh. Tn tht nc trong kênh là tn tht nc ca dòng
chy trên mt đt do hin tng thm nc qua b mt đáy và b kênh. Loi tn tht này
xy ra các dòng sui có nc liên tc hay không liên tc, ni mà nc ngm ch cung
cp nc cho dòng chy trong sui vào mt s thi đim nht đnh trong nm, hoc
không h cho mt chút nc nào. SWAT s dng phng pháp Lane đc mô t trong
Chng 19 ca S tay Thy vn SCS (Cc bo v đt đai, B Nông nghip Hoa k,
1983) đ c tính tn tht n
c trong sông sui. Tht thoát nc trong kênh là mt hàm
s ca chiu rng và chiu dài kênh và thi gian xy ra dòng chy. C tng lng dòng
chy tràn và đnh lu lng dòng chy đc điu chnh khi xy ra tn tht dòng chy
trong các kênh nhánh.
Dòng chy hi lu. Dòng chy hi lu, hoc dòng chy c bn, là tng lng dòng
chy có ngun gc t nc ngm. SWAT chia n
c ngm thành hai h thng tng nc
ngm: mt tng nc ngm đ sâu nông, có mt nc t do, không có áp, và mt tng
nc ngm sâu di đt, b ngn phía trên bi mt tng đa cht không thm nc, và
có đóng góp cho dòng chy ca dòng sông sui bên ngoài lu vc (Arnold et al., 1993).
Nc chuyn đng thm qua vùng r cây hot đng đc phân thành hai dng nc, mà
mi dng cung cp nc cho mt trong hai loi tng nc ngm. Ngoài dòng chy hi
lu, nc đc lu gi trong các tng nc ngm nông có th b sung đ m cho tng đt
trong điu kin đt tr nên rt khô hay trc tip đc r cây s dng. Nc trong tng
nc ngm nông hoc sâu đu có th đc con ngi khai thác bng bm n
c.
2.1.3 Lp ph thc vt / Sinh trng ca cây ci
SWAT s dng mô hình tng trng cây ci đ mô phng lp ph thc vt. Mô hình
có kh nng phân bit gia các cây hàng nm và cây lâu nm. Thi gian sinh trng ca
cây hàng nm là t ngày trng cây đn ngày thu hoch hoc cho đn khi tng tích ly các
đn v nhit bng con s tim nng các đn v nhit ca loi cây tr
ng đó. Cây lâu nm
duy trì h thng r ca chúng trong sut c nm, chúng không hot đng trong nhng
tháng mùa đông. Chúng tr li tip tc phát trin khi nhit đ không khí trung bình hàng
ngày vt quá giá tr ti thiu. Các mô hình sinh trng cây trng đc s dng đ tính
toán lng nc và các cht dinh dng đc b r cây chuyn khi vùng r cây hot
Gii thiu v SWAT (Soil and Water Assessment Tools)
đng, thoát đi qua quá trình hô hp ca cây, và tích ly, chuyn hóa thành sinh khi đ
hóa thành sn lng cây trng.
Tim nng phát trin ca cây trng. Kh nng tng trng sinh khi trong khong
thi gian mt ngày đc đnh ngha nh là lng gia tng sinh khi ca cây trong điu
kin sinh trng lý tng. Lng sinh khi gia tng tim nng trong trong mt ngày là
mt hàm s ca lng n
ng lng đc hp th và hiu sut chuyn đi ca cây t nng
lng thành sinh khi. Nng lng tích ly đc ca cây đc c tính nh là mt hàm
s ca bc x mt tri và ch s din tích lá ca cây trng.
Lng nc tn tht do hô hp tim nng và thc t. Quá trình tính toán tn tht
nc tim nng ca cây trng đã đc mô t trong phn bc thoát nc trên. Lng
thoát nc thc t do hô hp ca cây là mt hàm s ca bc thoát nc tim nng và
lng nc mà đt có th cung cp.
Dinh dng hp th. Lng nit và pht pho mà cây trng s dng, đc c tính
trên c s quan nim v lng cht cung cp đc và nhu cu ca cây, trong đó nhu c
u
hàng ngày v nit và photpho đc tính nh là lng khác bit gia nng đ thc t trong
các phn t cu thành cây trng và nng đ ti u. Nng đ ti u ca các cht dinh
dng bin đng theo các giai đon phát trin ca cây nh đc Jones (1983) mô t.
Nhng yu t hn ch cây trng phát trin. Tim nng tng trng và nng su
t cây
trng cao nht thng không đt đc do nhng hn ch áp đt bi môi trng. Mô hình
có th d báo nhng nh hng gây ra bi s thiu nc, thiu ht cht dinh dng và
bin đng ca nhit đ không khí đi vi sinh trng và phát trin ca cây trng.
2.1.4 Xói mòn đt
Xói mòn và dòng bùn cát đc mô hình c tính cho mi HRU s dng công thc
tng hp xói mòn đt đã thay đi (Modified Universal Soil Loss Equation, MUSLE;
Williams, 1975). Trong khi công thc USLE s dng lng ma nh là mt ch báo ca
nng lng gây xói mòn, thì MUSLE s dng tng lng nc tham gia dòng chy tràn
trên mt đt đ mô phng xói mòn và lu lng dòng bùn cát. Vic ci tin nh vy đa
đn cho ta mt s ích li: đ chính xác ca d báo xói mòn đt tt hn, không cn đn
thông s t l bào mòn, và có th tính toán lng đt xói mòn trong mi tr
n ma. Mô
hình thy vn cung cp các c tính th tích dòng chy và đnh lu lng dòng chy, mà
cùng vi din tích lu vc con, đc s dng đ tính toán bin s v nng lng xói mo欝n
ca dòng chy. Yu t ch đ canh tác đc tính toán li cho tng ngày mà dòng chy
tràn xy ra. ây là mt hàm s ca sinh khi trên mt đt, d tha trên b mt đt, và y
u
t C ti thiu cho loi cây trng. Các yu t khác ca phng trình xói mòn đc tính
toán và c lng nh mô t ca Wischmeier và Smith (1978).
2.1.5 Các cht dinh dng
SWAT mô phng chuyn đng và bin đi ca nhiu loi phân đm và pht pho
trong mt lu vc. Hình 6 din t vòng tun hoàn bin chuyn nit t các dng khác nhau
trong đt. Quá trình chuyn đng và bin đi ca pht pho trong đt c
ng đc trình bày
di dng vòng chu trình nh trong Hình 7. Nhng cht dinh dng này có th xâm nhp
vào dòng chy trên các kênh chính và đc vn chuyn xung vùng h lu theo chuyn
đng ca dòng chy mt và dòng nc ngm.
NIT. Hình 6 trình bày nguyên lý mô phng ca SWAT v các quá trình khác nhau
xy ra trong HRUs và các vùng khác nhau mà nit có th đc gi li trong đt. Lng
nit do cây ci hp th đc c tính bng cách s dng nguyên tc cân bng gia nhu
cu
nng
đi
c
l
tron
g
P
WAT
c
cây tr
n
chy
đ trung bì
n
đc tín
h
Williams
l
à
g
đt.
P
HT PHO
.
(Soil and Wat
e
c
ung cp
n
g). Ngoài
và
n
h ca nitra
t
h
toán vi m
và Hann (1
9
lp đ
à
nng đ c
6
.
Hình
a
e
r Assessment
mô t tr
o
t
h
. L
n
t
trong lp
t d
o
9
78)
.
t t
r
ên cùng,
a N hu c
Tools)
o
ng phn
h
c vt
n
g NO
3
-
N c
c tính
.
o
McElroy
e
áp d
n
c t
l
ng
trong
Nitrogen
, nitrat v
à
trong dòn
g
e
t al. (1976)
p
n
g cho
í
nh tn tht
,
v
c
h
SWA
T
sinh tr
à
nit hu c
g
chy,
d
ca
u c
p
hát trin v
à
v
à
t
l
h
ia cho n
T
.
N hu
ng
có th
d
òng chy
nc và
à
sa
. T
.
c.
Gii thiu v SWAT (Soil and Water Assessment Tools)
Hình 7 Quá trình bin chuyn Phosphorous nh đc SWAT mô phng.
2.1.6 Thuc tr sâu
Mc dù SWAT không mô phng tác đng lên sinh trng và phát trin ca cây trng
ca c di, côn trùng hi cây, và các loi côn trùng khác, nhng vic s dng thuc tr
sâu trên din tích mt HRU đc mô phng đ nghiên cu chuyn đng ca các cht hóa
hc nông nghip trong mt lu vc. SWAT mô phng xâm nhp ca thuc tr sâu vào
dòng chy sông ngòi qua dòng ch
y tràn trên mt đt (di dng hòa tan và dng cht b
dòng bùn cát hp th trc khi trôi theo dòng chy mt), và chuyn đng vào trong tng
sâu phu din đt và nc ngm (di dng hòa tan). Phng trình ch đo cho quá trình
chuyn đng ca thuc tr sâu trong vòng tun hoàn thy vn đã đc ly t phn mm
GLEAMS (Leonard et al., 1987). Vn chuyn ca thuc tr sâu b nh hng ca m
c đ
hòa tan, thi gian gim thiu 50%, và h s hp th carbon hu c. Lng thuc tr sâu
nm trên thân lá cây và trong đt s t tiêu hy theo hàm s m vi c ch thi gian gim
thiu 50%. Lng thuc tr sâu vn chuyn theo nc và bùn cát đc tính toán cho tng
đt chy tràn và lng thuc tr sâu thm xung đt đc c tính cho mi lp đt mi
khi hin tng thm sâu xy ra.
Gii thiu v SWAT (Soil and Water Assessment Tools)
Hình 8 Bin đi và chuyn đng ca thuc tr sâu nh đc SWAT mô phng.
2.1.7 Ch đ canh tác
SWAT cho phép ngi s dng đa d liu v các ch đ canh tác áp dng cho tng
HRU. Ngi s dng có th xác đnh ngày bt đu và kt thúc ca mi mùa v cây trng,
thi đim áp dng và lng phân bón, thuc tr sâu và nc ti cng nh ch đ làm
đt. Khi mt thi v cây trng kt thúc, toàn b lng sinh khi ca cây trng có th
đc đa ra khi các HRU di dng sn phm nông nghip hay dng d tha hu c.
Ngoài các bin pháp canh tác c bn nh trên, thì mô hình cng cho phép đa vào tác
đng ca nhng bin pháp khác nh khai thác đng c cho chn nuôi, bón phân và ti
nc t đng, hay bt c bin pháp đc bit nào khác. Kh nng x lý tác đng s dng
đt mi nht ca mô hình là v tính toán dòng bùn cát và cht dinh dng có ngun gc
t các đim dân c.
Luân canh. Theo các t đin thì luân canh có ngha là trng mt s lng các loi
cây trng khác nhau liên tip nhau cùng trên mt khu rung, thng là theo mt chu k
lp li. Khái nim luân canh cây trng trong SWAT đc hiu là vic thay đi ch đ
canh tác theo thi gian nm, ví d nm này có ch đ canh tác khác vi nm ti. Không
có gii hn v s nm mà mt ch đ canh tác đc lp li trong mt vòng luân canh.
SWAT cng cho phép không hn ch s lng loi cây trng trng trong mt nm ti mt
HRU. Tuy vy, ch có mt loi cây trng có mt trên mt ô rung trong mt thi đim
nào đó.
S dng nc. Hai loi hình s
dng nc ph bin nht là nc cho nông nghip và
nc cung cp cho sinh hot các đim dân c. SWAT cho phép lng nc áp dng
vào mt HRU có ngun gc t trong lòng hay t ngoài lu vc nghiên cu. Nc cng có
th đc vn chuyn gia các h cha, các khúc sông và các lu vc con hay có th đc
chuyn vn ra khi lu vc nghiên cu.
2.2 Tính truyn dòng chy và cht trong kênh
Mt khi SWAT xác đnh đc ngun ti v nc, bùn cát, dinh dng và thuc tr
sâu vào trong dòng kênh chính, thì ti lng nhng cht này đc tính truyn dc theo
mng li song ngòi ca lu vc, vi vic s dng mt cu trúc tính toán tng t nh
ca mô hình HYMO (Williams and Hann, 1972). Ngoài dòng chy khi lng trong
kênh, SWAT còn mô phng bin chuyn ca các cht hóa hc trong dòng nc và trong
lp bùn lng đáy kênh. Hình 9 mô t các quá trình khác nhau xy ra trong kênh đ
cc
SWAT mô phng.
Gii thiu v SWAT (Soil and Water Assessment Tools)
Hình 9 Các quá trình bin đi và vn chuyn cht trong lòng sông ngòi nh đc SWAT
mô phng
2.2.1 Tính truyn dòng chy nc và cht trong kênh chính
Tính truyn dòng chy và cht trong kênh chính có th đc chia làm bn thành
phn: nc, bùn cát, cht dinh dng và cht hóa hc hu c.
Tính truyn l. Khi dòng nc chy v phía h lu, mt lng nc có th b tn tht
do hin tng bay hi, và thm qua đáy kênh. Nhng hao t
n khác có th do nhng s
dng nc cho nông nghip, công nghip và sinh hot. Bên cnh đó dòng chy trong
kênh có th đc b sung bi nhng trn ma ri thng xung mt nc kênh hay/và bi
nhng ngun kênh nhánh. Lng dòng chy đc tính toán truyn theo chiu dài kênh
vi phng pháp h s tr nc bin đi ca Williams (1969) hoc phng pháp tính
truyn Muskingum.
Tính truyn dòng cht l lng
. Chuyn vn bùn cát trong kênh b hai quá trình cùng
xy ra đng thi chi phi: lng đng và xói l. Nhng phiên bn trc đây ca SWAT s
dng quan nim sc mnh dòng chy đ c tính lng lng đng và xói l trong kênh
(Arnold et al, 1995). Bagnold (1977) đnh ngha sc mnh dòng chy nh là tích s ca
mt đ nc, lu lng và đ dc mt nc. Williams (1980) s dng đnh ngha c
a
Bagnold đ phát trin phng pháp xác đnh xói l nh là hàm s ca đ dc đáy kênh và
vn tc nc.
Trong phiên bn này ca SWAT, phng trình trên đã đc đn gin hóa và lng
bùn cát ln nht mà dòng nc trong mt khúc sông có th vn chuyn đc ra khi khúc
sông đó là mt hàm s ca vn tc dòng chy ln nht. Nng lng ca dòng chy đc
s d
ng đ làm bùn cát tan ra và lng đng xung cho đn khi tt c bùn cát đc chuyn
đi. Nng lng còn li ca dòng chy đc s dng to nên xói mòn đáy kênh. Lng xói
đáy kênh đc tính toán chnh sa theo kh nng b xói ca vt cht đáy sông.
Gii thiu v SWAT (Soil and Water Assessment Tools)
Tính truyn vn chuyn cht dinh dng. Bin đi dng ca các cht dinh dng
trong dòng kênh đc x lý bi module cht lng nc trong kênh ca mô hình. SWAT
s dng c ch tính truyn cht trong kênh t mô hình QUAL2E (Brown and Barnwell,
1987). Mô hình mô phng cht dinh dng hòa tan trong nc kênh và b hp th bi bùn
cát l lng. Cht dinh dng hòa tan đc chuyn vn theo dòng nc trong khi lng
cht b bùn cát hp th đc l
ng đng cùng vi bùn cát xung đáy kênh.
Tính truyn thuc tr sâu. Thc t thì có vô vàn loi thuc tr sâu đc áp dng vào
din tích các HRUs nhng SWAT ch tính truyn đc mt loi thuc theo mng li
sông ngòi do các quá trình liên quan quá phc tp. Cng tng t nh đi vi cht dinh
dng, tng lng thuc tr sâu khi vào trong nc kênh đc chia thành hai dng: hòa
tan và hp th vào bùn cát. Trong khi lng thuc tr sâu hòa tan chuy
n vn theo dòng
nc thì lng hp th b chi phi bi hin tng lng đng và chuyn vn ca bùn cát l
lng. Quá trình bin đi cht ca thuc tr sâu trong hai dng hòa tan và hp th b chi
phi bi quan h gim thiu bc nht. Nhng quá trình đáng k mà mô hình mô phng
đc bao gm lng đng, chôn vùi, l lng, bay hi, khuyt tán và bin đi cht nh phn
ng hóa hc.
2.2.2 Tính truyn dòng chy và cht trong h cha
Cân bng nc trong h cha bao gm dòng chy vào, chy ra, lng ma ri trên
mt h, bc hi nc, thm xung qua mt đáy h và dòng dn nc ra khi h.
Dòng chy ra khi h cha. Mô hình cho phép s dng ba chn la đ c tính dòng
chy ra khi h cha. Ch
n la th nht cho phép ngi s dng nhp d liu quan trc
v lu lng dòng chy ra khi h. Chn la th hai dành cho h cha nh, không th
kim soát dòng chy ra, bt buc ngi s dng nhp d liu v lu lng dòng chy ra.
Khi th tích h cha vt khi th tích cha nc chính thì mô hình s s dng l
u lng
trên đ tiêu nc ra khi h. Phn th tích vt quá kh nng đp tràn an toàn s đc
đa ra khi h trong khong thi gian mt ngày. Chn la th ba, đc thit k cho các
h cha ln, thng đc xây dng vi các công trình kim soát dòng chy, cho phép
ngi s dng nhp d liu v th tích dòng ra tng tháng.
Tính truyn vn chuyn bùn cát. Dòng vn chuy
n bùn cát vào h có th có ngun
gc t nhng đon sông ngòi thng lu hay t dòng chy tràn mt đt trong lu vc
con. Nng đ bùn cát trong nc h đcc c tính vi mt phng trình đn gin da
trên nguyên lý cân bng th tích và nng đ trong dòng chy vào, chy ra khi h, và
lng nc tích tr trong h. Quá trình bùn cát lng đng trong h cha b chi phi bi
n
ng đ bùn cát cân bng và đng kính ht trung bình ca bùn cát. Lng bùn cát l
lng trong dòng chy ra khi h là tích s ca th tích nc chy ra và nng đ bùn cát l
lng trong nc h vào thi đim dòng chy ra khi h.
Cht dinh dng trong h cha. Mt mô hình đn gin da trên nguyên lý cân bng
khi lng nitrogen và phosphorus có ngun gc t Chapra (1997). Mô hình này có các
gi thit sau: 1) h nc đ
c trn ln hoàn toàn; 2) phosphorus là cht dinh dng
khng ch; và, 3) tng s phosphorus là s đo ca mc đ phù dng ca h. Gi thuyt
đu tiên đã b qua tính cht tng lp ca h nc và cng đ phát trin ca phù du
(phytoplankton) trong lp nc nông (epilimnon). Gi thuyt th hai thng là phù hp
vi thc t khi ngun ô nhim phân tán có vai trò ch đo và gi thuyt th ba ngh
a là có
tn ti quan h gia tng lng phosphorus và sinh khi. Phng trình cân bng khi
lng phosphorus bao gm nng đ trong nc h, dòng chy vào, dòng chy ra, và tc
đ thp thoát nc.
Gii thiu v SWAT (Soil and Water Assessment Tools)
Thuc tr sâu trong h cha. Mô hình cân bng thuc tr sâu trong h có gc gác t
Chapra (1997) và gi thuyt có tn ti điu kin trn ln hoàn toàn. H thng nc h
đc chia ra thành tng nc gn b mt trn ln hoàn toàn, và tng phía di trn ln
bùn cát hoàn toàn. Lng cht thuc tr sâu đc chia làm phn hòa tan vào nc và
phn cht rn l lng hay cht b
hp th vào ht bùn cát. Mô hình mô phng các quá
trình chính nh ngun cht, dòng chuyn vn ra khi h, bin đi dng cht, bc hi, lng
đng, khuyt tán, vn chuyn l lng, và chôn vùi.
SWAT-Intro_vie.docx 20
3. Tài liu tham kho
Arnold, J.G., P.M. Allen, and G. Bernhardt. 1993. A comprehensive surfacegroundwater
flow model. J. Hydrol. 142:47-69.
Arnold, J.G., J.R. Williams, A.D. Nicks, and N.B. Sammons. 1990. SWRRB: A basin
scale simulation model for soil and water resources management. Texas A&M Univ.
Press, College Station, TX.
Arnold, J.G., J.R. Williams and D.R. Maidment. 1995. Continuous-time water and
sediment-routing model for large basins. Journal of Hydraulic Engineering
121(2):171-183.
Bagnold, R.A. 1977. Bedload transport in natural rivers. Water Resources Res. 13(2):303-
312.
Brown, L.C. and T.O. Barnwell, Jr. 1987. The enhanced water quality models QUAL2E
and QUAL2E-UNCAS documentation and user manual. EPA document EPA/600/3-
87/007. USEPA, Athens, GA.
Chapra, S.C. 1997. Surface water-quality modeling. McGraw-Hill, Boston.
Green, W.H. and G.A. Ampt. 1911. Studies on soil physics, 1. The flow of air and water
through soils. Journal of Agricultural Sciences 4:11-24.
Hargreaves, G.L., G.H. Hargreaves, and J.P. Riley. 1985. Agricultural benefits for
Senegal River Basin. J. Irrig. and Drain. Engr. 111(2):113-124.
Jones, C.A. 1983. A survey of the variability in tissue nitrogen and phosphorus
concentrations in maize and grain sorghum. Field Crops Res. 6:133-147.
Knisel, W.G. 1980. CREAMS, a field scale model for chemicals, runoff and erosion from
agricultural management systems. USDA Conservation Research Rept. No. 26.
Leonard, R.A. and R.D. Wauchope. 1980. Chapter 5: The pesticide submodel. p. 88-112.
In Knisel, W.G. (ed). CREAMS: A field-scale model for chemicals, runoff, and
erosion from agricultural management systems.
U.S. Department of Agriculture, Conservation research report no. 26.
Leonard, R.A., W.G. Knisel, and D.A. Still. 1987. GLEAMS: Groundwater loading
effects on agricultural management systems. Trans. ASAE 30(5):1403-1428.
McElroy, A.D., S.Y. Chiu, J.W. Nebgen, A. Aleti, and F.W. Bennett. 1976. Loading
functions for assessment of water pollution from nonpoint sources. EPA document
EPA 600/2-76-151. USEPA, Athens, GA.
Monteith, J.L. 1965. Evaporation and the environment. p. 205-234. In The state and
movement of water in living organisms. 19th Symposia of the Society for
Experimental Biology. Cambridge Univ. Press, London, U.K.
Nicks, A.D. 1974. Stochastic generation of the occurrence, pattern and location of
maximum amount of daily rainfall. p. 154-171. In Proc. Symp. Statistical Hydrology,
Tucson, AZ. Aug Sept. 1971. USDA Misc. Publ. 1275. U.S. Gov. Print. Office,
Washington, DC.
Priestley, C.H.B. and R.J. Taylor. 1972. On the assessment of surface heat flux and
evaporation using large-scale parameters. Mon. Weather Rev. 100:81-92.
Gii thiu v SWAT (Soil and Water Assessment Tools)
Ritchie, J.T. 1972. A model for predicting evaporation from a row crop with incomplete
cover. Water Resour. Res. 8:1204-1213.
USDA Soil Conservation Service. 1983. National Engineering Handbook Section 4
Hydrology, Chapter 19.
USDA Soil Conservation Service. 1972. National Engineering Handbook Section 4
Hydrology, Chapters 4-10.
Williams, J.R. 1980. SPNM, a model for predicting sediment, phosphorus, and nitrogen
yields from agricultural basins. Water Resour. Bull. 16(5):843-848.
Williams, J.R. 1975. Sediment routing for agricultural watersheds. Water Resour. Bull.
11(5):965-974.
Williams, J.R. 1969. Flood routing with variable travel time or variable storage
coefficients. Trans. ASAE 12(1):100-103.
Williams, J.R., A.D. Nicks, and J.G. Arnold. 1985. Simulator for water resources in rural
basins. Journal of Hydraulic Engineering 111(6): 970-986.
Williams, J.R. and R.W. Hann. 1978. Optimal operation of large agricultural watersheds
with water quality constraints. Texas Water Resources Institute, Texas A&M Univ.,
Tech. Rept. No. 96.
Williams, J.R. and R.W. Hann. 1972. HYMO, a problem-oriented computer language for
building hydrologic models. Water Resour. Res. 8(1):79-85.
Williams, J.R., C.A. Jones and P.T. Dyke. 1984. A modeling approach to determining the
relationship between erosion and soil productivity. Trans. ASAE 27(1):129-144.
Wischmeier, W.H., and D.D. Smith. 1978. Predicting rainfall losses: A guide to
conservation planning. USDA Agricultural Handbook No. 537. U.S. Gov. Print.
Office, Washington, D. C.
DANH MC HÌNH
Hình 1 Bn đ lu vc h Fork đông bc Texas vi phân b s dng đt và mng
li sông ngòi. 4
Hình 2 Chia nh lu vc h Fork. 5
Hình 3 S đ vòng tun hoàn thy vn. 6
Hình 4 Vòng tính toán cho HRU/Lu vc con. 7
Hình 5 S đ đng chuyn đng ca nc đc SWAT s dng. 10
Hình 6 Quá trình bin chuyn Nitrogen đc SWAT mô phng. 14
Hình 7 Quá trình bin chuyn Phosphorous nh đc SWAT mô phng. 15
Hình 8 Bin đi và chuyn đng ca thuc tr sâu nh đc SWAT mô phng. 16
Hình 9 Các quá trình bin đi và vn chuyn cht trong lòng sông ngòi nh đc
SWAT mô phng 17