Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống
1
/ 18 trang
THÔNG TIN TÀI LIỆU
Thông tin cơ bản
Định dạng
Số trang
18
Dung lượng
367,14 KB
Nội dung
1 V các mô hình toán ca dòng chy TS Tô Vn Trng I. KHÁI NIÊM V MÔ HÌNH 1. Các lo"i mô hình Ngoài các loi mô hình toán hc, trong các sách và tài liu tham kho ta còn gp nhng t khác nh!: mô hình t# l hay mô hình v%t lý, mô hình t!'ng t(, mô hình nh%n th)c, mô hình t*t +,nh, mô hình ng-u nhiên, i) Mô hình t$ l% (mô hình v(t lý): Thay cho làm th(c nghiêm trên mô hình nguyên m-u ngoài th(c +,a, ng!2i ta ti3n hành thí nghim trên các mô hình thu nh5 nh!ng gi nguyên t# l gia các chi6u (dài, r9ng, cao, ). T*t nhiên, phi tuân th; m9t s< nguyên t=c và tiêu chu>n. Ph!'ng pháp này th!2ng r*t +=t, +òi h5i nhi6u th2i gian +@ xây mô hình. Ph!'ng pháp này ch; y3u +!Bc áp dCng +@ nghiên c)u chi ti3t khi th%t cDn thi3t +@ thi3t k3 công trình nh! c<ng +%p. ii) Mô hình t+,ng t-: Vì các ph!'ng trình mô t dòng chy ngDm (th*m) t!'ng t( nh! ph!'ng trình mô t dòng +in, cho nên thay vì nghiên c)u mng dòng chy ph)c tp ng!2i ta l%p m9t mng +in và suy các k3t qu t mng +in sang mng dòng chy. iii) Mô hình nh(n th/c hay khái ni%m: Tr!Gc khi nghiên c)u m9t quá trình ng!2i ta phi xem quá trình +ó có bao nhiêu thành phDn, cách th)c liên h và nh h!Kng l-n nhau c;a các thành phDn +ó, sau +ó xem xét chi ti3t t ng thành phDn. iv) Mô hình t1t 23nh: Là m9t ph!'ng pháp nghiên c)u khi bi3t các thành phDn tham gia vào quá trình và cách th)c nh h!Kng gia các thành phDn. 2. Mô hình v(t lý và mô hình toán h6c Môi tr!2ng th(c r*t +a +ng và ph)c tp, các m<i quan h gia các y3u t< +an xen chOng ch,t nh h!Kng l-n nhau. P@ kho sát hoc nghiên c)u các m<i quan h +ó +ã t lâu con ng!2i phi sR dCng mô hình hoá nh! là m9t công cC, có nghSa là phi +'n gin hoá b)c tranh th(c hoc môi tr!2ng th(c. Mô hình không bao gi2 ch)a +!Bc t*t c các +c +i@m c;a môi tr!2ng th(c mà chT gi li các +c +i@m chính các m<i quan h chính c;a h th<ng th(c mà chúng có th@ +c tr!ng cho h th<ng +ó. Ví dC: Khi xem xét ch*t l!Bng n!Gc ng!2i ta chT xem xét m9t vài chT tiêu ch; y3u nh! nh! +9 pH, nhu cDu ô xy sinh hóa BOD, +9 c)ng, Eli-Coliform. Khi xem xét n!Gc bi@n và n!Gc sông ta chT cDn xem xét +9 mn. T!'ng t( nh! vây, khi thi3t k3 m9t con tDu ng!2i ta th!2ng làm các mô hình v%t lý +@ xem xét các hình dng nào có s)c cn nh5 nh*t ch) ch!a cDn chú ý tGi vic b< trí ca bin, hDm tDu. Nh! v%y, quá trình mô hình hoà là quá trình xem xét +@ chT cDn gi li các +c +i@m chính +c tr!ng cho môi tr!2ng hoc v%t nào +ó cDn phi nghiên c)u. Quá trình làm m-u m9t con tDu thu# có kích th!Gc theo m9t t# l nào +ó r`i cho vào thR trong n!Gc vGi m9t s< +i6u kin v6 sóng gió, +!Bc xem là mô hình v%t lý. M9t 2 loi mô hình khác th!2ng có tên là mô hình toán sb +!Bc gii thích kc trong các phDn d!Gi, nh!ng có th@ hi@u nôm na là, các m<i quan h gia các hin t!Bng, s( v%t hay các y3u t< c;a môi tr!2ng bao gi2 cdng có th@ bi@u dien bOng các quan h (hay ph!'ng trình) toán hc. BOng cách nghiên c)u hoc gii các ph!'ng trình toán +ó ng!2i ta có th@ phát hin +!Bc các tính ch*t c;a các hin t!Bng cDn quan tâm xem xét. Mô hình hoá bOng các mô hình toán +!Bc phát tri@n r*t nhanh trong các th%p niên gDn +ây, bKi vì: • S( phát tri@n nh! vd bão c;a công ngh máy tính và công cC tin hc, +`ng th2i các công cC mGi v6 toán hc cdng phát tri@n. Hai y3u t< này giúp cho con ng!2i có th@ gii quy3t r*t nhanh các bài toán ph)c tp v6 mt toán hc. • Mt khác, yêu cDu v6 phát tri@n kinh t3, xã h9i và dân s< d-n +3n vic suy thoái môi tr!2ng, +c bit vic ô nhiem môi tr!2ng n!Gc d-n +3n +e do s( s<ng trên hành tinh trong t!'ng lai gDn. Vì v%y, xu*t hin các bài toán ph)c tp v6 mt môi tr!2ng mà chT có công cC mô hình hoá mGi có th@ d( báo +!Bc bi3n +gi có th@ x>y ra. 3. Mô hình là công c8 qun lý. S( +ô th, hoá, s( bùng ng dân s< và s( phát tri@n công ngh ngày càng gia ting áp l(c và tác +9ng lên môi tr!2ng ta +ang s<ng. Các ch*t ô nhiem thi vào h sinh thái, +c bit là h sinh thái n!Gc, +ang làm gia ting các ch*t +9c hi c v6 hóa, lý, sinh, có th@ hu# dit các loài s<ng trong +ó hoc phá vj c*u trúc c;a h sinh thái. H sinh thái ngày nay r*t ph)c tp, nhim vC c;a chúng ta là d( +oán s( bi3n +gi c;a môi tr!2ng d!Gi tác +9ng c;a các y3u t< khác nhau. Trong b<i cnh +ó mô hình hoá sb cung c*p m9t b)c tranh vGi các +áp )ng khác nhau, t +ó có nhng bin pháp qun lý và thích )ng cdng nh! l(a chn các gii pháp công ngh hay pháp lý thích h'p trong xR lý cdng nh! qun lý. 4. Mô hình hoá là m<t công c8 khoa h6c Mô hình là loi công cC +!Bc sR dCng r9ng rãi trong khoa hc. Các nhà khoa hc tr!Gc +ây sR dCng r9ng rãi các mô hình v%t lý +@ ti3n hành các thí nghim ngoài hin tr!2ng cdng nh! trong phòng thí nghim +@ nghiên c)u các m<i quan h chính mà ng!2i ta quan tâm. Ngày nay, do s( phát tri@n c;a công ngh máy tính và công ngh thông tin xu th3 phg bi3n là làm các thí nghim trên máy tính tr!Gc khi ti3n hành b=t bu9c m9t s< thí nghim v%t lý nhOm ki@m +,nh các k3t qu t máy tính, và do +ó mô hình toán +!Bc sR dCng r*t r9ng rãi. P,nh lu%t Newton mà ta quen thu9c K m9t phm vi nào +ó là m9t mô hình toán hc v6 s( nh h!Kng c;a l(c trng tr!2ng lên v%t th@ khi b5 qua l(c ma sát và nh h!Kng c;a gió. Do tính ph)c tp c;a môi tr!2ng và h sinh thái vic mô hình hoá là b=t bu9c +@ khám phá ra m<i liên h gida các y3u t< và t!'ng tác gia các y3u t<. Chlng hn khi xem xét s( ô nhiem n!Gc sông ch,u nh h!Kng c;a thu# tri6u. Khi thi ch*t b>n vào dòng chy, d!Gi tác +9ng c;a thu# tri6u ch*t thi b>n lan to +i các h!Gng khác nhau và cdng gim dDn n`ng +9 b>n do quá trình t( làm sch. P@ tính toán +!Bc phm vi nh 3 h!Kng c;a các ngu`n ô nhiem cdng nh! n`ng +9 ti t ng th2i +i@m thì chT có mô hình hoá mGi gii quy3t +!Bc. Có th@ tóm l!Bc m9t s< !u +i@m c;a công cC mô hình hoá nh! sau: - Là công cC hu ích và không th@ thi3u trong kho sát các h sinh thái ph)c tp - SR dCng mô hình có th@ khám phá ra các tính ch*t c;a h th<ng. - Nh2 mô hình có th@ hoàn thin s( hi@u bi3t v6 ki3n th)c môi tr!2ng, sinh thái. - Mô hình là m9t công cC +@ thR nghim các gi thuy3t v6 khoa hc và so sánh gia b)c tranh th(c và b)c tranh c;a môi tr!2ng +ã +!Bc +'n gin hoá. 5. Th? nào là mô hình toán: Trong nhi6u lSnh v(c hot +9ng hàng ngày ta phi th(c hin các tính toán t +'n gin tGi ph)c tp. n tr!2ng phg thông phi th(c hin các phép c9ng tr nhân chia, r`i cao h'n là các phép +o hàm, vi tích phân. R`i trong các tr ong Pi hc phi hc các ph!'ng pháp s< nh! sai phân hu hn, phDn tR hu hn. Nói chung, ta có th@ gi chung là các công cC toán hc và sR dCng chúng +@ gii quy3t các bài toán trong th(c t3 hàng ngày, t +'n gin +3n phúc tp. VGi s( phát tri@n r*t nhanh c;a công ngh thông tin, kS thu%t máy tính và các công cC toán hc hin +i, mô hình toán hc +ã +!Bc phát tri@n r*t nhanh và +ã trK thành công cC nhanh mnh, không th@ thi3u +<i vGi nhng ng!2i làm công tác qui hoch và ra quy3t +,nh. V%y mô hình toán hc là gì ? P@ gii thích ta xét m9t ví dC sau +ây: t1 t i t2 Gi sR ta phi tính khong cách S c;a ôtô chy trong th2i gian T t th2i +i@m t1 tGi t2 vGi v%n t<c t)c th2i v(t). Cách tính +'n gin sb nh! sau: Ta chia khong th2i gian T thành n khong nh5 h'n vGi b!Gc th2i gian là t i , có nghSa là : = =++++= n i in tttttT 1 321 Trong m9t b!Gc th2i gian t i ta xem v%n t<c v(t) gDn nh! không +gi vGi giá tr, v i , và khong cách S mà xe chy trong khong th2i gian T có th@ +!Bc tính x*p xT nh! sau: 11 2 2 1 . . . . (2.1) n nn ii i Svt vt vt vt = + ++ = Dùng công th)c (2.1) bài toán có th@ xem nh! +!Bc gii quy3t (m9t cách gDn +úng). N3u t i +; nh5 thì t ki3n th)c hc trong Pi hc ta có : 4 2 1 0 1 lim . . (2.2) t n ii t i t Svtvdt = = = Nh! v%y vGi các cán b9 +ã t<t nghip +i hc bài toán tính quãng +!2ng c;a xe chy +!Bc tính toán theo các b!Gc nh! sau: i) Xu*t phát t công th)c (2.2) : 2 1 . (2.2) t t Svdt= ii) Ti3n hành r2i rc hoá (2) theo dng (1): 11 2 2 1 . . . . (2.1 ) n nn ii i Svt vt vt vt a = + ++ = Trong +ó v%n t<c v i c xem là hOng s< trong khong t i và bOng giá tr, trung bình c;a v%n t<c t)c th2i v i trong b!Gc th2i gian t i . iii) N3u ta gp khó khin khi tính bOng tay (chlng hn +<i vGi các bài toán lGn ph)c tp) ta có th@ l%p trình trên máy tính +@ tính toán. N3u ta cm th*y k3t qu tính toán nói trên ch!a +; chính xác ta có th@ ti3p tCc làm t<t h'n theo cách sau:. iv) Làm nh5 h'n t i (t)c là ting n) và tính li công th)c trên bOng cách dùng cùng m9t ch!'ng trình máy tính +ã vi3t. Quá trình này +!Bc lp li (t)c là ting dDn n) cho +3n khi k3t qu thu +!Bc có th@ xem là +t yêu cDu. 6. Các b+Bc trong xây d-ng m<t mô hình toán: Quá trình tính toán vGi các b!Gc +'n gin nh! nêu trong ví dC chi3c ô tô K trên +!Bc xem nh! các b!Gc xây d(ng môt mô hình toán (tr!2ng hBp r*t +'n gin). Trên th(c t3 ta gp r*t nhi6u bài toán kS thu%t ph)c tp, nhi6u khi không th@ tính toán bOng tay +!Bc, vì th3 vic xây d(ng m9t quá trình tính toán (hoc xây d(ng môt mô hình toán hc) sb r*t ph)c tp, t<n công, nh!ng nói chung sb g`m các b!Gc sau +ây: i) B+Bc 1: L(a chn các ph!'ng trình toán hc c' bn mô t các quá trình v%t lý (hoc bài toán ta phi gii quy3t. Nói chung hDu h3t các bài toán th(c t3 +6u có th@ mô t bOng các ph!'ng trình toán hc). VGi bài toán tính khong cách c;a ô tô nh! nêu K trên thì +ó là ph!'ng trình (2). P<i vGi các bài toán kS thu%t thì các ph!'ng trình c' bn (hoc h ph!'ng trình) +6u là các ph!'ng trình vi phân, tích phân hay +o hàm riêng mà +@ gii chúng cDn phi có s( giúp +j c;a ph!'ng pháp s< và máy tính. Thông th!2ng +@ có +!Bc các ph!'ng trình c' bn mô t m9t quá trình v%t lý nào +ó ta th!2ng áp dCng các lu%t bo toàn cho các quá trình v%t lý +ó nh! bo toàn kh<i l!Bng, bo toàn mô men +9ng l!Bng hay bo toàn ning l!Bng. Nguyên lý bo toàn chung cho m9t +i l!Bng b*t ký (th@ tích, n`ng +9, kh<i l!Bng,…) sb nh! sau: S bin i theo thi gian ca mt i lng bt k trong mt th tích V s$ b%ng tng lng vào th tích tr' i tng lng ra kh(i th tích cng (ho*c tr') v,i lng phát 5 sinh (ho*c mt i) do các nguyên nhân khác nhau trong chính th tích V ó. Ch6ng hn 7i v,i n,c trong th tích V Q1 Q2 2 121 SSQQ dt dV += Hình 2: S' +` cân bOng n!Gc cho th@ tích V VGi Q1, Q2 là l!u l!Bng vào ra ti 2 mt th@ tích; S1 là ngu`n n!Gc bg xung (x n!Gc vào) còn S2 là ngu`n n!Gc b, l*y +i. Pây là nguyên lý +!Bc sR dCng khi thi3t l%p ph!'ng trình liên tCc c;a ph!'ng trình Saint-Venant. VGi BOD nguyên lý cân bOng trên +!Bc vi3t nh! sau BOD vào th tích V + BOD s=n sinh trong V – (BOD ra kh(i V + BOD b? chuy n hoá) = s thay i BOD trong th tích V trong kho=ng thi gian t. VGi B là n`ng +9 BOD, Q là l!u l!Bng ti mt c=t, g là t<c +9 sn sinh, f là t<c +9 m*t +i, +an phát bi@u trên +!Bc tóan hc hóa bOng bi@u th)c sau: . . BVB QB gV Q B x fV xt + + += Hình 3: S' +` cân bOng n`ng +9 BOD trong th@ tích V Chlng hn +@ mô t chuy@n +9ng c;a n!Gc và +9 mn trên kênh sông ng!2i ta th!2ng dùng h ph!'ng trình Saint-Venant m9t chi6u cho dòng chy và ph!'ng trình ti khu3ch tán cho +9 mn nh! d!Gi +ây: Ph!'ng trình liên tCc cho n!Gc (bo toàn th@ tích n!Gc): q, S q (2.3) HQ Bq tx += Ph!'ng trình +9ng l!'ng (bo toàn mô men +9ng l!Bng) : H A Hình : q, S q gV-fV Q.B V. B/ t [] B QB x x + 6 2 2 0 (2.4) gQQ QQ H gA txA x AC R +++ = 0 Hình 4: Mt c=t ngang sông Ph!'ng trình liên tCc cho +9 mn (bo toàn kh<i l!Bng ): 2 2 () .(2.5) oiq AS QAS S EA q S q S tx x += + Trong +ó : H = M(c n!Gc so vGi cao +9 chu>n (m); Q = l!u l!Bng (m 3 /s); B = +9 r9ng mt n!Gc ti m9t mt c=t ngang sông bao g`m c phDn tr (m); A = din tích mt c=t ngang (m 2 ) C = H s< cn Chezy ; g = gia t<c trng tr!2ng (m/s 2 ); R = bán kính thu# l(c (m); q =q i -q o : dòng gia nh%p dc dòng chy (q i ) hoc m*t +i (q o ) trên m9t +'n v, +9 dài c;a dòng chy (m 2 /s) t = th2i gian (s) x = khong cách dc dòng chy (m) S(x,t) : P9 mn (hay n`ng +9 ch*t) trung bình trên mt c=t ngang (g/L) E : H s< phân tán dc (dispersion) Sq: P9 mn (hay n`ng +9 ch*t) trong dòng gia nh%p ii) B+Bc 2: P<i vGi các bài toán ph)c tp mô t bKi các ph!'ng trình +o hàm riêng thì +@ gii +!Bc cDn phi cho +i6u kin biên, +i6u kin +Du, các tham s< và các h s<. iii) B+Bc 3: Nói chung các ph!'ng trình mô t các quá trình vât lý hDu nh! không có nghim gii tích hoc nghim chính xác (theo nghSa toán hc) vì th3 phi dùng các ph!'ng pháp s< +@ gii gDn +úng. Pi6u +ó có nghSa rOng bài toán chT +!Bc gii gDn +úng và k3t qu thu +!Bc cdng là k3t qu gDn +úng ch) không phi k3t qu chính xác. S( khác nhau gia k3t qu chính xác và k3t qu gDn +úng phC thu9c vào ph!'ng pháp s< +!Bc sR dCng. Có r*t nhi6u ph!'ng pháp s<, vic l(a chn ph!'ng pháp nào phC thu9c vào trình +9 và ki3n th)c c;a ng!2i l%p mô hình. Ví dC, K trên là m9t thu%t toán s< r*t +'n gin +@ tính quãng +!2ng chy c;a ôtô. P@ ting +9 chính xác giá tr, v%n t<c trung bình v i có th@ l*y theo các cách khác nhau, chlng hn l*y giá tr, trung bình +Du +on và cu<i +on, hoc l*y giá tr, gia +on. iv) B+Bc 4: VGi s( phát tri@n r*t nhanh c;a kc thu%t máy tính hDu h3t các ph!'ng pháp s< +6u có th@ th(c hin trên máy tính vGi +i6u kin thu%t toán s< t!'ng )ng +ã +!Bc l%p trình và chy thông +!Bc trên máy tính. Pây là b!Gc không th@ thi3u +!Bc khi xây d(ng m9t mô hình toán hc 7 v) B+Bc 5: ThR tính +úng +=n c;a k3t qu qua m9t s< bài toán m-u +@ bo +m rOng k3t qu phn ánh t!'ng +<i chính xác các qui lu%t v%t lí (vì ta chT tính gDn +úng) c;a quá trình +!Bc mô ph5ng. Chlng hn tính bo toàn kh<i l!Bng, nh! cân bOng n!Gc, hoc +9 mn không th@ âm. N3u +9 mn tính ra b, âm thì có +i6u gí +ó sai trong thuât toán +!Bc sR dCng. M9t ví dC khác trong thR nghim là tính +<i x)ng. N3u t*t c các +i6u kin c;a bài toán (mi6n, biên, +i6u kin +Du, ) là +<i x)ng thì nghim s< c;a nó cdng +<i x)ng. N3u th*y k3t qu tính ra có sai sót thì li phi xem xét t b!Gc 1. VGi 5 b!Gc chính nh! +ã nêu K trên ta có m9t mô hình toán hc. P9 chính xác c;a k3t qu tính toán phC thu9c vào chính mô hình và ch*t l!Bng c;a s< liu +Du vào. Vì v%y khi sR dCng b*t c) m9t mô hình nào (phDn m6m máy tính) +@ gii quy3t m9t bài toán th(c tien hai b!Gc nêu d!Gi +ây cDn phi +!Bc th(c hin: • Hiu chTnh mô hình: Khi m9t mô hình +ã +!Bc xây d(ng nó có th@ sR dCng cho b*t kì m9t bài toán kc thu%t nào nh!ng ch*t l!Bng c;a k3t qu phC thu9c vào ch*t l!Bng s< liu +Du vào và giá tr, c;a các tham s< c;a mô hình. D(a trên m9t s< s< liu +Du vào +!Bc +o +c xác +,nh và hiu chTnh các tham s< (chlng hn h s< nhám trong mô hình thu# l(c) +@ có +!Bc l2i gii t<t nh*t. Quá trình này gi là hiu chTnh mô hình. • Ki@m +,nh mô hình: M9t khi mô hình +ã d!Bc hiu chTnh t<t cDn phi thR vGi m9t t%p s< liu khác +@ ki@m tra xem vGi các tham s< mô hình +ã +!Bc xác +,nh liu có +úng vGi tr!2ng hBp khác không, n3u k3t qu tính không sai nhi6u vGi k3t qu th(c +o thì mô hình có th@ coi là +!Bc ki@m +,nh và dùng +@ tính toán các k,ch bn khác nhau II. MÔ HÌNH TÓAN CHO DÒNG CHMY VÀ CHPT LRSNG NRTC TRÊN HV THWNG KÊNH SÔNG – MÔ HÌNH DELTA (VRSAP-SAL) 2.1 M] 2^u Hin ti, +@ tính dòng chy ld kit, xâm nh%p mn, trng thái ô nhiem hu c', trên các h th<ng kênh sông c;a Vit nam, ch; y3u là P`ng bOng sông CRu Long (PBSCL), h th<ng sông Sài gòn-P`ng Nai-Th, vi, P`ng bOng sông H`ng, sông H!'ng các kc s! và cán b9 kc thu%t th!2ng dùng m9t s< phDn m6m máy tính c;a n!Gc ngoài và trong n!Gc. V6 mt hc thu%t các mô hình tính dòng chy và ch*t l!Bng n!Gc trong sông +6u xu*t phát t h ph!'ng trình Saint-Ve nant 1 chi6u (K các dng khác nhau) và ph!'ng trình lan truy6n ch*t m9t chi6u. Tuy nhiên, s' +` và thu%t tóan gii các h ph!'ng trình này li khác nhau tùy thu9c tác gi c;a t ng mô hình, t +ó +9 chính xác c;a k3t qu cdng nh! th2i gian tính trên máy có khác nhau. 2.2 Các Mô hình t_ n+Bc ngoài +!Bc du nh%p vào Vit nam theo con +!2ng các d( án (trong +ó các phDn m6m kèm theo +!Bc tính vào ti6n d( án, t)c là phi mua phDn m6m) hoc bOng con +!2ng c;a du hc sinh hoc hBp tác song ph!'ng. 2.2.1 Nhóm mô hình th+,ng m"i: Pây là nhóm mô hình mua tr(c ti3p hoc tính thành ti6n thông qua các d( án song ph!'ng hoc +a ph!'ng: 8 A. Nhang mô hình dòng chy và ch1t l+bng n+Bc có tính th+,ng m"i trên th? giBi phi kf 2?n h6 mô hình MIKE, trong 2ó MIKE11 (vGi mô+un thu# l(c HD, mô +un tính mn, ch*t l!Bng n!Gc AD, ECOLAB, ) Pây là b9 phDn m6m c;a Vin DHI Pan Mch, +!Bc )ng dCng, nghiên c)u cho d( án quy hoch và qun lý tài nguyên n!Gc và phòng ch<ng thiên tai ti nhi6u n!Gc trên th3 giGi nh! Nh%t Bn, Thái Lan, Bangdales Trong khuôn khg c;a D( án ting c!2ng ning l(c các Vin Ngành n!Gc K Vit Nam, DHI +ã +ào to và chuy@n giao bn quy6n cho m9t s< c' quan ngành n!Gc thu9c B9 NN&PTNT. M9t s< phDn m6m h MIKE khác nh! MIKEBASIN (dùng cho tính cân bOng n!Gc), MIKE FLOOD dùng cho mô ph5ng ld, MIKE21 dùng cho bài toán n<i 1 chi6u và 2 chi6u trong m9t vùng nh5, V6 b9 MIKE 11: MIKE 11 là phDn m6m thu9c h MIKE vGi modun tính dòng chy HD và modun AD dùng cho tính lan truy6n ch*t (mn,…), +@ tính lan truy6n ch*t ô nhiem phi dùng ECOLAB vGi các y3u t< lan truy6n ch*t t th*p tGi cao. P@ tính dòng chy trong sông kênh MIKE 11 cdng sR dCng h ph!'ng trình Saint-Venant m9t chi6u và sR dCng s' +` sai phân 6 +i@m xen kb Q, H c;a Abbott và Ionescu; tài liu +,a hình +!Bc cho ti các mt c=t tính H; v%n t<c u +!Bc tính ti +i@m Q; H ph!'ng trình sai phân +!Bc gii tr(c ti3p và bOng ph!'ng pháp lp, vì v%y t<c +9 tính ch%m và cDn có kinh nghim xR lý khi to +i6u kin ban +Du (hotstart file). Trong MIKE 11 +ã xét các công trình c<ng +%p phg bi3n, tuy nhiên +ôi khi khi gp tr!2ng hBp không gn +,nh khi phi v%n hành công trình. Trong modun AD +ã sR dCng ph!'ng pháp sai phân hu hn cho ph!'ng trình lan truy6n ch*t m9t chi6u vì th3 th!2ng gp hin t!Bng khu3ch tán s< nh h!Kng +3n +9 chính xác c;a k3t qu nh! n`ng +9 có khi b, âm, hoc khi không có ngu`n sinh v%t ch*t trong mi6n mà n`ng +9 trong mi6n cao h'n giá tr, K biên,… P@ sR dCng công cC GIS, trong MIKE 11 +ã dùng k3t hBp vGi b9 ArcView/Arcview GIS +@ tg ch)c c' sK d liu và bi@u dien k3t qu (thông qua các script bOng ngôn ng Avenue). Nhìn chung nhng !u nh!Bc +i@m c;a b9 MIKE 11 (+!Bc sR dCng nhi6u K Vit nam ch; y3u qua d( án ting c!2ng ning l(c cho các Vin ngành n!Gc) nh! sau: + †u +i@m: - Là phDn m6m th!'ng mi nên phDn giao din r*t mnh, hu hiu. - PhDn n<i k3t vGi công cC GIS r*t mnh k@ c to Database (Mc dù phi cDn thêm các phDn m6m GIS nh! ArcView hay ArcGIS, ) - Các tin ích +Dy +;, de cho ng!2i sR dCng. - Thu%n tin cho vic gii quy3t các bài toán v a và nh5. + Nh!Bc +i@m: - Không bi3t +!Bc phDn lõi (phDn thu%t toán, tg ch)c ch!'ng trình, ) nên ng!2i sR dCng không th@ ci biên, c%p nh%t mà phi qua n'i bán, khi +ó phi tr thêm ti6n và m*t th2i gian ch2 +Bi, ) - Khi phi tính cho bài toán lGn nh! PBSCL trong m9t th2i gian dài (mô ph5ng c m9t nim cho ld và cn) MIKE 11 +òi h5i nhi6u th2i gian tính trên máy không 9 thu%n tin cho giai +on chy hiu chTnh vì phi chy r*t nhi6u lDn mGi hiu chTnh +!Bc m9t tham s< nên t<n th2i gian chy trên máy. H'n na, +@ to +i6u kin ban +Du (hotstart file) +òi h5i nhi6u kinh nghim và th!2ng phi xu*t phát t b!Gc th2i gian nh5 - P9 chính xác c;a k3t qu tính, +c bit cho các bài toán lan truy6n ch*t (mn, BOD, DO, ) nhi6u khi không +m bo do bn ch*t thu%t toán +!Bc sR dCng (khu3ch tán s< d-n +3n n`ng +9 âm hoc n`ng +9 sát biên lGn h'n biên khi không có ngu`n trong mi6n) - Vì là phDn m6m th!'ng mi nên giá thành r*t +=t (MIKE11+ECOLAB giá 18000EU, cj 400 triu +`ng Vit nam cho m9t license) mŠi license, dng khoá c)ng, chT dùng +!Bc cho m9t máy tính, hoc cdng có phiên bn chy n<i k3t máy tính trên mng nh!ng giá thành cao h'n nhi6u. - Nhi6u nghiên c)u trong n!Gc +ã sR dCng mô hình MIKE11 +@ làm công cC tính toán thu# l(c và ch*t l!Bng n!Gc. Nh!ng sau khi hoàn thành d( án không chuy@n giao công ngh +!Bc vì các c' quan h!Kng lBi t d( án không có bn quy6n sR dCng MIKE11 và d( án cdng th!2ng không có +; kinh phí +@ mua phDn m6m chuy@n giao. B. ISIS: B9 phDn m6m này c;a Công ty Halcrow và tr!2ng Wallingford ph<i hBp xây d(ng, +!Bc sR dCng trong ch!'ng trình sR dCng n!Gc (WUP) c;a U# H9i sông Mê Công. MŠi m9t n!Gc thành viên có +!Bc 2-3 license. Tuy phDn m6m này, +<i vGi Vit Nam, ch!a th!'ng mi hoá nh! MIKE , nh!ng du nh%p vào Vit nam thông qua các d( án có th@ chuy@n giao công ngh nh! Ch!'ng trình WUP nói trên +ây. Gi<ng nh! b9 MIKE11, phDn m6m ISIS cdng sR dCng h ph!'ng trình Saint- Venant m9t chi6u cho dòng chy và ph!'ng trình lan truy6n ch*t m9t chi6u cho mn. Khác vGi MIKE 11, trong ISIS sR dCng s' +` sai phân Preissmann cho dòng chy và lan truy6n mn. Cdng nh! MIKE 11 phDn m6m ISIS ch!a có kh ning tính mn trong +`ng. Vì là phDn m6m th!'ng mi, ISIS cdng có phDn giao din khá +Œp và tin dCng, tny nhiên cdng b9c l9 m9t s< y3u +i@m và khó khin khi gii quy3t bài toán trên phm vi r9ng, nhi6u liên k3t nh! PBSCL. Pc bit các lŠi v6 +9 chính xác c;a k3t qu tính, +c bit v6 mn. Trong khuôn khg c;a ch!'ng trình WUP, phDn m6m ISIS +ã +!Bc sR dCng cho PBSCL K dng mng kênh sông +!Bc +'n gin hoá r*t nhi6u (b5 mng kênh c*p 2, ch; y3u gi li dòng chính), nh!ng ch!a cho k3t qu có th@ sR dCng +!Bc, +c bit là phDn tính mn. Mô +un ch*t l!Bng n!Gc v-n ch!a +!Bc thR nghim nên ch!a có k3t qu +ánh giá cC th@. T<c +9 tính tóan c;a ISIS cdng r*t ch%m và cdng k3t hBp vGi ArcView +@ n<i k3t vGi GIS và Database. Nh! v%y khi mua MIKE 11 hoc ISIS phi tr c ti6n bn quy6n c;a ArcView. 2.2.2. Nhóm mô hình phi th+,ng m"i (theo nghia Vi%t nam ch+a phi mua mà có 2+bc qua các con 2+jng khác nhau nh+ d- án hk trb song ph+,ng homc 2ào t"o) Các b9 phDn m6m khác nh! Duflow, Sobek/Wendy,Telemax, Qual2-E, Wasp6 vv +!Bc du nh%p qua các các con +!2ng c;a du hc sinh hoc các d( án nh5 song ph!'ng. P<i vGi các d( án qu<c t3 thì +ây cdng là các b9 phDn m6m th!'ng mi, phi mua bn quy6n nên khi sR dCng th!2ng +!Bc c' quan c*p phDn m6m khuy3n cáo rOng 10 có th@ ch*p nh%n m9t s< r;i ro gây thit hi do không +!Bc +ào to, t%p hu*n và không hi@u bi3t nhng hn ch3 c;a mô hình nên khi áp dCng gây lŠi. Vì không có mã ngu`n nên không hi@u +!Bc h3t phDn lõi bên trong xR lý ra sao (nh! thu%t tóan, các xR lý +c bit, ) và ch!a +!Bc áp dCng cho các bài toán lGn và ph)c tp nh! PBSCL. Các phDn m6m này có ngu`n g<c t châu Âu (hoc Mc) vGi +i6u kin sông ngòi khác hln +i6u kin Vit nam (chlng hn K Vit nam mng kênh sông có dng mch vòng ph)c tp, ch,u nh h!Kng c;a th;y tri6u, ) cho nên không phi khi nào cdng sR dCng +!Bc các phDn m6m nêu trên. Có th@ xét qua phDn m6m Sobek, Duflow và Qual2-E: SOBEK: PhDn m6m này do Delft,Hà lan, phát tri@n, g`m phDn dòng chy và tính tóan ô nhiem 1,2 chi6u, +ã n<i k3t vGi công cC GIS. Pã sR dCng h ph!'ng trình Saint- Venant 1 chi6u cho dòng chy trong kênh sông (trong ph!'ng trình có k@ s< hng gió và nh h!Kng c;a góc nh%p l!u). SOBEK cdng sR dCng l!Bc +` sai phân xen kb gi<ng nh! MIKE11, có +i@m H và +i@m Q; +,a hình +!Bc cho ti các +i@m tính H. Các y3u t< ô nhiem +!Bc mô ph5ng bOng ph!'ng trình lan truy6n ch*t 1 chi6u có k@ tGi quá trình bi3n +gi sinh hóa c;a các ch*t ô nhiem. Ph!'ng trình lan truy6n ch*t m9t chi6u +!Bc gii bOng ph!'ng pháp sai phân, mc dù có các l(a chn các s' +`, nh!ng do bn ch*t c;a l!Bc +` sai phân, k3t qu tính v-n b, nh h!Kng bKi hin t!Bng khu3ch tán s<. Qual2-E: PhDn m6m này do c' quan bo v môi tr!2ng c;a Mc (EPA) phát tri@n và +ã +!Bc sR dCng r9ng rãi K Mc và m9t s< n!Gc châu Âu. Qual2-E +ã +!Bc du nh%p vào Vit Nam qua m9t s< d( án. Qual2-E cdng sR dCng h ph!'ng trình Saint-Venant và lan truy6n ch*t m9t chi6u và gii bOng ph!'ng pháp sai phân và có th@ sR dCng cho nhi6u y3u t< ô nhiem (BOD, DO,To, Nit', Ph<t pho, ). Nh!Bc +i@m c;a Qual2-E là chT áp dCng cho mng sông +'n gin có dng hình cây (không áp dCng cho mng sông dng mch vòng); thi3t din kênh sông phi +6u dng hình thang, hay hình ch nh%t và không ch,u nh h!Kng c;a th;y tri6u Duflow: Pây là phDn m6m +!Bc phát tri@n bKi Vin th;y l(c (IHE) c;a Hà lan, Pi hc công ngh Delft, STOWA và tr!2ng Pi hc nông nghip Wageningen. Duflow +!Bc thi3t k3 +@ sR dCng cho nhi6u mCc tiêu (tính tri6u, ld, sR dCng n!Gc, ). Duflow cdng gii quy3t các bài tóan lan truy6n ch*t trong kênh sông có các công trình. S' +` sai phân 4 +i@m c;a Preissmann +ã +!Bc sR dCng cho bài tóan th;y l(c. Duflow có giao din +` ha tin dCng. Vì +ây là phDn m6m thi3t k3 ch; y3u cho ging d%y và +ào to, cho nên khi sR dCng cho các bài tóan lGn cDn có ci biên. 2.3 Mô hình trong n+Bc: Do các yêu cDu c;a th(c tien quy hach và sR dCng tài nguyên n!Gc, nhi6u chuyên gia trong n!Gc phi t( xây d(ng các b9 phDn m6m, +@ khi cDn thi3t, có th@ t( sRa +gi và c%p nh%t thu%t tóan, mã ngu`n (code) +@ có th@ +áp )ng +!Bc các yêu cDu tính tóan cC th@. Các b9 phDn m6m do các cán b9 trong n!Gc +!Bc nh=c tên và áp dCng nhi6u cho các d( án trên 2 P`ng bOng g`m: VRSAP, +ây là b9 phDn m6m +!Bc xem là +Du tiên cho tính tóan th;y l(c mng kênh sông, do c< PGS Nguyen nh! Khuê phát tri@n sau +Bt th(c t%p ti Hà Lan vào nim 1978. VRSAP +ã +!Bc Phân vin Kho sát Quy hoch Th;y lBi Nam b9 (Nay là Vin [...]... (GIS) V6 c' b n ch!'ng trình VRSAP +ã +áp )ng +!Bc các yêu cDu tính toán, tuy nhiên do nhu cDu phát tri@n, kích cj c;a các bài toán qui ho ch cdng ting dDn, không chT K m)c +9 P`ng bOng c;a Vi t nam mà K m)c +9 châu thg (chlng h n c Vi t Nam và Cim pu chia) và ph i mô t vGi th2i gian dài và vGi các k,ch b n ph)c t p v6 v%n hành các h th . nh!: mô hình t# l hay mô hình v%t lý, mô hình t!'ng t(, mô hình nh%n th)c, mô hình t*t +,nh, mô hình ng-u nhiên, i) Mô hình t$ l% (mô hình v(t lý): Thay cho làm th(c nghiêm trên mô hình. 1 V các mô hình toán ca dòng chy TS Tô Vn Trng I. KHÁI NIÊM V MÔ HÌNH 1. Các lo"i mô hình Ngoài các loi mô hình toán hc, trong các sách và tài liu tham. tham gia vào quá trình và cách th)c nh h!Kng gia các thành phDn. 2. Mô hình v(t lý và mô hình toán h6c Môi tr!2ng th(c r*t +a +ng và ph)c tp, các m<i quan h gia các y3u t< +an xen