(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long(Luận văn thạc sĩ) Nghiên cứu sự biến động của bùn cát vùng đồng bằng sông Cửu Long
M CL C M U Ch ng 1: T NG QUAN 1.1 T ng quan vùng nghiên c u 1.1.1 V trí vƠ u ki n đ a lý t nhiên 1.1.2 c m đ a hình 1.1.3 a ch t, th nh ng 1.1.4 Ch đ khí h u 11 1.1.5 Ch đ th y v n 13 1.2 Các nghiên c u liên quan đ n bùn cát vùng BSCL 25 1.2.1 Các nghiên c u n 1.2.2 Các nghiên c u n 1.3 CH c 25 c 26 Tình hình thay đ i bùn cát vùng BSCL 27 NG PH NG PHÁP NGHIÊN C U 30 2.1 Ph ng pháp thu th p phân tích t ng h p thông tin, d li u 31 2.2 Ph ng pháp th ng kê 31 2.3 Ph ng pháp mơ hình tốn 33 2.3.1 Module thu đ ng l c 36 2.3.2 Module v n chuy n bùn cát l l ng 41 2.3.3 Thi t l p mơ hình 2D-FEM cho vùng nghiên c u 46 CH 3.1 NG K T QU TÍNH TỐN 55 Các tiêu chí đánh giá sai s 55 i 3.2 K t qu tính tốn hƠm l ng phù sa t t ng quan 56 3.2.1 K t qu xây d ng t 3.2.2 S thay đ i bùn cát l l ng theo th i gian 61 3.2.3 S thay đ i bùn cát l l ng theo không gian 69 3.3 ng quan t i tr m 56 K t qu mơ t mơ hình 2D-FEM 72 3.3.1 K t qu mô ph ng thu l c 72 3.3.2 K t qu mô ph ng hƠm l ng phù sa 78 K T LU N VÀ KI N NGH 84 K t lu n 84 Ki n ngh 85 TÀI LI U THAM KH O 87 ii DANH M C CÁC HỊNH NH Hình 1.1: B n đ vùng ng b ng sông C u Long khu v c ph c n Hình 1.2: B n đ th hi n đ a hình vùng tính tốn ảình 1.3: S ho đ d c d c sông Mê Kông t ngu n sông đ n c a Hình 1.4: Bi u đ l ng m a trung bình tháng t i m t s tr m 11 Hình 1.5: B n đ phân b l Hình 1.6: M ng l ng m a trung bình n m l u v c sông Mê Kông 12 i sông ngịi, kênh r ch vùng Hình 1.7: H th ng th y h vùng ng B ng sông C u Long 14 ng b ng sông c u long 15 Hình 1.8: nh v tinh th hi n khu v c ng p tr n l tháng 9/2000 (Vi n KHKTTV&MT) 17 Hình 1.9: B n đ ti m n ng tr l ng n c ng m vùng BSCL (Vi n KHKTTV&MT - 2010) 19 Hình 1.10: Phân b dịng ch y ki t tính tốn theo mơ hình tri u bán nh t (Hà lan, 1974) 20 Hình 1.11: M ng l i tr m th y v n l u v c 24 Hình 1.12: Di n bi n hàm l ng phù sa, ch t r n l l ng tr m Kratie 28 Hình 1.13: Di n bi n hàm l ng phù sa, ch t r n l l ng tr m Kratie bình quân tháng 28 ảình 1.14: Thay đ i hàm l ng phù sa Kratie tháng đ u mùa m a 29 Hình 2.1: S đ th hi n ph ng pháp n i dung nghiên c u lu n v n 30 Hình 2.2: S đ t ng qt c a mơ hình ph n t h u h n hai chi u 2D-FEM 35 iii Hình 2.3: S đ th hi n đ c tr ng dòng ch y h to đ ph ng trình đ c tr ng c a module thu đ ng l c 37 Hình 2.4: S đ q trình tính tốn c a module thu đ ng l c 40 Hình 2.5: S đ q trình tính tốn c a module v n chuy n bùn cát 45 Hình 2.6: S đ th hi n sông nhánh đ vào vùng bi n h Tonle Sap (Campuchia) 47 Hình 2.7: Ải i h n vùng tính tốn (bên ph i) ví d v ô l i tam giác (bên trái) dùng đ th hi n l i vùng tính tốn 47 ảình 2.8: L i tính tốn c a vùng đ ng b ng sơng h l u sông Mê Kông, g m 128,815 nút 255,996 ô l i tam giác 49 Hình 2.9: B n đ th hi n đ a hình vùng tính tốn 50 ảình 2.10: ng trình m c n c t i tr m Kratie n m 2000 51 ảình 2.11: ng q trình dịng ch y c a 12 nhánh nh p l u đ vào bi n h Tonle Sap n m 2000 52 ảình 2.12: ng trình dịng ch y c a nhánh nh p l u đ vào bi n h Tonle Sap n m 2011 53 ảình 2.13: ng trình m c n c t i tr m Kratie n m 2011 53 ảình 2.14: ng trình m c n c t i biên h l u n m 2000 54 ảình 2.15: ng trình m c n c t i biên h l u n m 2011 54 Hình 3.1: Quan h gi a bùn cát l l ng l u l ng t i tr m Kratie 58 Hình 3.2: Quan h gi a bùn cát l l ng l u l ng t i tr m Tân Châu 58 Hình 3.3: Quan h gi a bùn cát l l ng l u l ng t i tr m Châu Hình 3.4: Quan h gi a bùn cát l l ng l u l ng t i tr m Vàm Nao 59 Hình 3.5: Quan h gi a bùn cát l l ng l u l ng t i tr m M Thu n 60 c 59 iv Hình 3.6: Quan h gi a bùn cát l l ng l u l Hình 3.7: ng t i tr m C n Th 60 ng trình bùn cát l l ng (trên) l u l ng n c (d i) t i Kratie 62 Hình 3.8: ng trình bùn cát l l ng (trên) l u l ng n c (d i) t i Tân Châu 63 Hình 3.9: Châu ng trình bùn cát l l ng (trên) l u l ng n c (d i) t i c 64 Hình 3.10: ng trình bùn cát l l ng (trên) l u l ng n c (d i) t i Vàm Nao 65 Hình 3.11: ng trình bùn cát l l ng (trên) l u l ng n c (d i) t i tr m M Thu n 66 Hình 3.12: ng trình bùn cát l l ng (trên) l u l ng n c (d i) t i tr m C n Th 67 Hình 3.13: Bùn cát l l ng hay hàm l ng phù sa trung bình n m t i m t s v trí vùng nghiên c u [5] 68 Hình 3.14 : ảàm l ng phù sa trung bình tháng dịng sơng Mê Kơng 68 Hình 3.15: N ng đ bùn cát l l ng th i k khác n m 2000 70 Hình 3.16: N ng đ bùn cát l l ng th i k khác n m 2011 71 Hình 3.17: ng q trình m c n c tính tốn th c đo t i: a) Kompong Cham, b) Neak Luong, c) M Thu n d) C n Th 73 Hình 3.18: K t qu mơ ph ng đ sâu dịng ch y t i th i m 12 gi ngày 2509-2000 74 ảình 3.19: ng trình m c n c tính tốn th c đo t i C n Th cho ki m đ nh mơ hình 75 v Hình 3.20: ng trình m c n c tính tốn th c đo t i M Thu n cho ki m đ nh mơ hình 76 Hình 3.22: K t qu mơ ph ng phân b đ sâu v n t c dịng ch y vùng tính tốn t i th i m 12:00:00 ngày 01 tháng 10 n m 2011 77 Hình 3.23: (d ng trình bùn cát l l ng t i Tân Châu (trên) Châu c i) 81 vi DANH M C B NG BI U B ng 1.1: L ng m a trung bình tháng, n m t i m t s tr m [3] 13 B ng 1.3: T ng h p th ng kê chu i s li u th y v n tr m th y v n dòng 23 B ng 1.5: T i l ng phù sa trung bình nhi u n m dịng sơng Mê Kơng t i tr m Kratie [5] 27 B ng 2.1: Tên kí hi u chu i s li u m c n B ng 3.1: T c t i biên h l u 51 ng quan gi a bùn cát l l ng l u l ng n c t i m t s tr m vùng nghiên c u 56 B ng 3.2: Bùn cát l l ng (hàm l ng phù sa) trung bình tháng t i m t s v trí vùng nghiên c u [5] 68 vii M U Sông Mê Kông đ c x p hàng th h th ng sông l n nh t th gi i, có ngu n n c t ng đ i d i v i t ng l ng n c bình quơn hƠng n m kho ng 475 t m3 Ph n l n lãnh th c a Lào (97%) Campuchia (86%) n m l u v c Mê Kơng , 36% di n tích c a Thái Lan (vùng ông B c) n m l u v c Ph n di n tích c a Vi t Nam chi m kho ng 11%, có ph n vùng Tơy Nguyên vƠ vùng đ ng b ng sông C u Long ( BSCL), 82% t ng l ng dòng ch y đ c hình thành t b n n c h l u: LƠo 35%, Thái Lan 18%, Camuchia 18% vƠ Vi t Nam 11% Ph n dòng ch y đóng góp t hai n c th ng l u chi m kho ng 18% t ng l ng dòng ch y c n m, đ c đánh giá lƠ ph n đóng góp quan tr ng dịng ch y mùa ki t cho l u v c ng b ng sông C u Long vùng đ ng b ng quan tr ng nh t c a n c ta, v i di n tích t nhiên vào kho ng tri u (x p x 1/9 di n tích c a c n c), đư đóng góp h n 50% s n l ng l ng th c (là n n t ng an ninh l ng th c Qu c gia), h n 90% s n l ng g o xu t kh u, chi m h n 70% s n l ng nuôi tr ng th y s n 57% t ng s n l ng th y s n c a c n c Tuy m t đ ng b ng có giàu ti m n ng, nh ng v n nhi u v n đ t n t i liên quan đ n n c nh : l l t hƠng n m, xói l b sơng, ven bi n, xâm nh p m n, h n hán, đ t phèn vƠ n c phèn, ô nhi m ngu n n c BSCL n m h ngu n l u v c sông Mê Kông, tác đ ng phát tri n th ng l u nh gia t ng di n tích nơng nghi p, xây d ng th y n, đ c bi t phát tri n th y n dòng đ c xem lƠ nguy c đe , nh h ng đ n ngu n l i phù sa th y s n, có th lƠm gia t ng xói l b bi n đ i lòng d nầ đơy lƠ nguy c cho s phát tri n n đ nh c a vùng đ ng b ng sông C u Long Trong nh ng n m g n đơy, hi n t ng xói l phía h du vùng ng b ng sơng C u Long phát tri n m nh gây nh h ng x u đ n khu v c dơn c Có nhi u ngun nhân gây tình tr ng s t l nghiêm tr ng M t nh ng nguyên nhân gây nên hi n t ng xói l đ c cho lƠ l ng bùn cát đ bi n t th ng ngu n sông Mê Kông suy gi m Các nghiên c u cho th y, t n t i nhi u nguyên nhân d n đ n s suy gi m l ng bùn cát h th ng sông Lu n v n đ c đ t v i m c tiêu tìm hi u s bi n đ ng c a bùn cát l l ng vùng ng b ng sông C u Long M c dù đư có m t s nghiên c u v s bi n đ ng bùn cát l l ng h th ng sông Mê Kông vƠ đánh giá nguyên nhơn gơy s bi n đ ng Tuy nhiên, ph n l n nghiên c u ch y u t p trung vào s đánh giá s bi n đ ng bùn cát l l ng t i v trí c th vƠ xác đ nh mƠ ch a có nghiên c u t p trung vào vi c mô ph ng cho tồn vùng h l u sơng Mê Kơng c ng nh toƠn b vùng ng b ng sông C u Long s d ng mơ hình tốn ph c t p nh mơ hình hai chi u ho c ba chi u ph c v cho nhu c u nói trên, c n thi t ph i có m t cơng c h tr k thu t đ đánh giá đ c tình hình bi n đ ng c a bùn cát vùng đ ng b ng sơng C u Long Do đó, đ tài lu n v n “Nghiên c u s bi n đ ng c a bùn cát vùng ng b ng sông C u Long” đư đ a l a ch n C ng c n ph i nh n m nh r ng lu n v n s ch y u t p trung nghiên c u tìm hi u s bi n đ ng c a bùn cát l l ng h t m n có kích th c nh h n 62 m, c s phân tích s bi n đ i l u l ng bùn cát quan tr c đ c t i tr m th y v n (khi s d ng ph ng pháp th ng kê ậ quan h t ng quan) NgoƠi ph ng pháp th ng kê ra, lu n v n c ng m nh d n s d ng mơ hình ph n t h u h n hai chi u 2D-FEM đ mô ph ng thu l c vƠ bùn cát l l ng vùng nghiên c u Ngoài ph n m đ u ki t lu n, lu n v n đ c b c c nh sau Ch ng trình bày khái qt chung v ng b ng sơng C u Long Ch ng gi i thi u v ph ng pháp th ng kê dùng đ xây d ng t ng quan gi a l u l ng n c vƠ hƠm l ng phù sa, ph ng pháp mơ hình tốn c th mơ hình hai chi u 2D-FEM Ch ng trình bày c th k t qu tính tốn c a lu n v n CH NG 1: T NG QUAN 1.1 T ng quan vùng nghiên c u 1.1.1 V trí u ki n đ a lý t nhiên Sông Mê Kông lƠ m t nh ng sông l n nh t th gi i, b t ngu n t cao nguyên Tơy T ng (Trung Qu c) ch y qua n c (Trung Qu c, LƠo, Myanmar, Thái Lan, Campuchia vƠ Vi t Nam) tr c đ Bi n ông, v i chi u dƠi sơng vƠ di n tích l u v c l n l t lƠ 4.800 km 795.000 km2 [1, 2, 3, 4] L u v c sông Mê Kông tr i dài t 90 đ n 350 V B c 930 đ n 1070 Kinh ơng (Hình 1.1) L u v c sơng Mê Kơng có th chia thƠnh vùng lƠ th ng l u vƠ h l u (Hình 1.1) Vùng th ng l u gi i h n t ngu n sông kéo dƠi đ n biên gi i Trung Qu c, Myanmar vƠ LƠo, v i chi u dƠi sơng 3.000 km, chi m kho ng 24% di n tích toƠn l u v c vƠ đóng góp kho ng 18% l ng n c ch y vƠo sông Mê Kông Trong vùng th ng l u, l u v c sơng h p, dịng ch y m nh, lịng sơng h p vƠ sơu, nhi u gh nh thác vƠ ch y qua nhi u vùng núi cao có đ a hình ph c t p Xói mịn lƠ m t v n đ l n vùng nƠy c tính kho ng 50% tr m tích sơng Mê Kơng b t ngu n t vùng th ng l u L ng tr m tích nƠy theo dịng ch y di chuy n xu ng h l u, m t ph n b i l ng sông, hay bưi b i ven sơng q trình di chuy n, m t ph n tích t h Tonle Sap vƠ ph n l i b i l ng vùng thu c ng b ng sông C u Long Vùng h l u gi i h n t B c Vientiane (LƠo) đ n Bi n ơng, có th chia thƠnh vùng gi i h n t B c Vientiane (LƠo) đ n StungtrengậKratie (Campuchia) vùng đ ng b ng gi i h n t Kratie đ n Bi n ông - Vùng gi i h n t B c Vientiane (Lào) đ n Stungtreng–Kratie (Campuchia) có chi u dƠi sơng 750 km vƠ chi m kho ng 57% di n tích toƠn l u v c Trong vùng nƠy, sông ch y song song v i dưy Tr ng S n b ng qua m t cao nguyên sa th ch kh ng l v i t ng đ a ch t n m ngang t ng n, sông nh n n c t ph l u nh sông N m Re, N m U, N m Suông, N m Ng m, N m Th ng, Sê Bang Phai, Sê Bang Hiên, Sê Pôn h u ng n, sông nh n n c t ph l u N m Mum (bao trùm cao nguyên Carat) vƠ ph l u Mênam Xongkhram o n nƠy có hai thác r t l n lƠ thác Kemmarat (có d ng m t h m v c) vƠ thác Khone (r t hi m tr ) Trong vùng gi i h n t B c Vientiane (LƠo) đ n StungtrengậKratie (Campuchia) thu c lưnh th Campuchia, sông Mê Kông nh n n c t (i) ph l u nh sông SêKong, SêSan, Sêrêpok (t Tơy Nguyên Vi t Nam đ xu ng t ng n) vƠ (ii) h TonLe Sap ( Tơy B c Campuchia) đ vƠo TonLe Sap có ch đ sơng h , đơy t n t i m t h n c kh ng l gi a dòng TonLe Sap (v i chi u dƠi 150 km Hình 3.20: ng q trình m c n c tính tốn (Htt) th c đo (Htđ) t i M Thu n cho ki m đ nh mơ hình Hình 3.21: K t qu mô ph ng phân b đ sâu v n t c dịng ch y vùng tính tốn t i th i m 12:00:00 ngày 25 tháng n m 2011 Hình 3.21 th hi n phân b đ sâu v n t c dòng ch y t i th i m 12:00:00 ngƠy 25 tháng n m 2011, Hình 3.22 ví d minh ho k t qu mơ ph ng phân b đ sâu v n t c dịng ch y vùng tính tốn t i th i m 76 12:00:00 ngƠy 01 tháng 10 n m 2011 K t qu tính tốn th hi n r ng ph m vi kho ng 100 km h l u dịng sơng Mê Kơng gi i h n t Kraite đ n th đô Phnom Penh (C mpuchia), v n t c dòng ch y t ng đ i l n, v i biên đ t i m t s v trí lên đ n 2,0 m/s Trong v n t c dịng ch y nh xu t hi n vùng bi n h Tonle Sap Trong vùng ng b ng sông C u Long vùng bãi ng p d c theo dịng sơng Ti n sơng H u c ng nh sông nhánh c a sơng này, v n t c dịng ch y thay đ i kho ng t đ n 0,5 m/s Trong vùng c a sông c a ng b ng sơng C u Long, n i có s nh h ng m nh m c a thu tri u, v n t c dịng ch y dao đ ng m t d i l n h n, v i biên đ lên đ n 1,0 m/s Hình 3.22: K t qu mơ ph ng phân b đ sâu v n t c dòng ch y vùng tính tốn t i th i m 12:00:00 ngày 01 tháng 10 n m 2011 77 T ng t nh k t qu mô ph ng dịng ch y hi u ch nh thơng s mơ hình, k t qu mơ ph ng đ sơu vƠ tr ng phân b v n t c dòng ch y b c ki m đ nh c ng đư th hi n đ c di n tích khơ/ t khác vùng tính tốn i u có ngh a r ng mơ hình đư cho phép th hi n di n tích khô/ t khác c ng nh cho phép xem xét s thay đ i ph c t p c a đ a hình vùng tính tốn r ng l n c a vùng h l u sông Mê Kơng (tính t Kraite đ n vùng c a sơng ng b ng sông C u Long) Nh đư th hi n Hình 3.22, m c n c ngày mô ph ng th hi n s dao đ ng nh c a m c n c so v i th c đo S khác bi t v dao đ ng m c n c có th nhi u nguyên nhơn khác nh vi c s d ng h s nhám h ng s su t th i gian mô ph ng, s phơn chia vƠ nhơp l u dòng ch y t i ngã ba sông, nh h ng c a thu tri u c ng nh vi c b qua dòng ch y th c p tính tốn 3.3.2 K t qu mô ph ng hàm l ng phù sa 3.3.2.1 K t qu hi u ch nh mô ph ng bùn cát l l ng vùng tính tốn, giá tr c a ng su t t i h n cho h t bùn cát b t đ u b i l ng d = 0.05 N/m2 đ c s d ng ng th i, ng su t t i h n cho h t bùn cát b t đ u b xói đ c gi đ nh b ng ng su t t i h n cho h t bùn cát b t đ u b i l ng hay nói cách khác d = c H s t l xói M = 510-6 kg/s/m2 t c đ l ng chìm c a h t bùn cát w = 1.3 mm/s [26] Hình 3.23: N ng đ bùn cát l l ng th c đo tính tốn t i Tân Châu 78 Hình 3.24: N ng đ bùn cát l l ng th c đo tính tốn t i Châu ơc Hình 3.25: N ng đ bùn cát l l ng th c đo tính toán t i C n Th 79 K t qu tính tốn n ng đ bùn cát l l ng cho hi u ch nh thơng s mơ hình module v n chuy n bùn cát đ c th hi n hình t Hình 3.23 đ n Hình 3.25 D dàng nh n th y r ng mơ hình th hi n t t xu th thay đ i c a n ng đ bùn cát l l ng t i tr m: Tân Chơu, Chơu c vƠ Cơn Th th i k so sánh t ngày 27 tháng 05 đ n ngƠy tháng n m 2011 K t qu mô ph ng c ng th hi n đ c s thay đ i c a n ng đ bùn cát l l ng t i th i m khác th i k so sánh Tuy nhiên, k t qu tính tốn n ng đ bùn cát l l ng thiên l n th i k n ng đ bùn cát nh (t ngày 27-05-2011 đ n ngày 15-06-2011) th i k cịn l i k t qu mơ ph ng t mơ hình l i thiên nh Ngun nhân d n đ n s khác bi t xu th thay đ i có th vi c gi đ nh ng su t t i h n cho h t bùn cát b t đ u b xói b ng ng su t t i h n cho h t bùn cát, ho c có th vi c s d ng giá tr c a thơng s mơ hình moduel v n chuy n bùn cát ch a th t s xác 3.3.2.2 K t qu ki m đ nh Hình 3.26 Hình 3.27 l n l t th hi n k t qu mô ph ng n ng đ bùn cát l l ng t i Tơn Chơu vƠ Chơu c K t qu tính tốn n ng đ bùn cát l l ng t mơ hình hai chi u 2D-FEM đ c so sánh v i k t qu tính tốn theo t ng quan L u Ủ r ng s h n ch v s li u bùn cát vƠ hƠm l ng phù sa n m 2000 nên k t qu tính tốn t mơ hình m i ch d ng l i vi c so sánh v i k t qu c tính bùn cát l l ng t ph ng trình t ng quan theo ph ng pháp th ng kê D dàng nh n th y r ng k t qu th hi n s t ng đ ng v s thay đ i hƠm l ng phù sa theo th i gian gi a ph ng pháp VƠo th i gian có l u l ng dịng ch y l n hàm l ng phù sa tính tốn theo c ph ng pháp đ u cho k t qu l n Ng c l i vào th i k có l u l ng dịng ch y nh k t qu mơ ph ng t mơ hình tốn c ng cho hƠm l ng phù sa nh Hình 3.28 th hi n phân b n ng đ bùn cát l l ng vùng tính tốn t i th i th i m 12:00:00 ngƠy 25 tháng n m 2000 k t qu mô ph ng phân b n ng đ bùn cát l l ng vùng tính tốn t i th i điêm 12:00:00 ngƠy 01 tháng 10 n m 2000 đ c th hi n nh Hình 3.29 K t qu mô ph ng n ng đ bùn cát l l ng th hi n rõ s thay đ i theo khơng gian vùng tính tốn N ng đ bùn cát l l ng l n th ng xu t hi n t i v trí thu c dịng c a sơng Mê Kơng c ng nh t i v trí có t c đ dịng ch y l n, t i v trí bãi ng p thu c vùng ng b ng C mpuchia vƠ vùng ng b ng Sơng C u Long n ng đ bùn cát l l ng th ng nh b i t i v trí t c đ dịng ch y nh 80 Hình 3.26: ng trình bùn cát l l ng t i Tân Châu Hình 3.27: ng trình bùn cát l l ng t i Châu c 81 Hình 3.28: K t qu mơ ph ng phân b n ng đ bùn cát l l ng vùng tính tốn t i th i m 12:00:00 ngày 25 tháng n m 2000 82 Hình 3.29: K t qu mơ ph ng phân b n ng đ bùn cát l l ng vùng tính tốn t i th i m 12:00:00 ngày tháng 10 n m 2000 83 K T LU N VÀ KI N NGH K t lu n Bùn cát m t đ i l ng không th thi u chuy n đ ng c a dòng ch y vùng ng p n c nói chung trrong vùng ng b ng sông C u Long nói riêng Nghiên c u s thay đ i hay bi n đ ng c a bùn cát (nói chung bùn cát l l ng nói riêng) theo khơng gian th i gian có Ủ ngh a vô to l n th c s c n thi t B i vì, th nh t, s thay đ i c a bùn cát dòng ch y vùng ng b ng sông C u Long có s liên h ch t ch v i s màu m c a di n tích tr ng nh t di n tích đ t nông nghi p Th hai, nghiên c u chuy n đ ng c a bùn cát l l ng không ch cho phép xem xét đ c tr ng c a dòng ch y nh v n t c, đ sâu dòng ch y, l u l ng n c mƠ cịn có Ủ ngh a quan tr ng nghiên c u ch t l ng n c Th ba, nghiên c u s chuy n đ ng c a bùn cát cịn giúp ích cho vi c xem xét q trình b i xói sơng d c theo sông Th t , nghiên c u s thay đ i c a bùn cát có th giúp ích cho vi c đánh giá xu th tang gi m hàm l ng phù sa n c, qua có th cung c p thơng tin h u ích cho vi c đánh giá s t ng gi m di n tích đ t cho nơng nghi p ó lƠ lý đ tài nghiên c u s bi n đ ng c a bùn cát vùng ng b ng sông C u Long đ c l a ch n th c hi n lu n v n Ngoài n i dung t ng quan v u ki n đ a lý t nhiên, th nh ng, khí t ng khí h u, thu v n, m t s k t qu nghiên c u c a lu n v n có th tóm t t nh sau: HƠm l ng phù sa trung bình n m t i Kratie, Tơn Chơu vƠ Chơu c cho th i k t n m 1985 đ n 2016 có xu th gi m rõ r t, nh t k t n m 1998 Nguyên nhân d n đ n s suy gi m hƠm l ng phù xa s ho t đ ng c a h ch a thu n th ng ngu n sông Mê Kông vƠo khai thác v n hƠnh, nh đ p thu n Manwan (Trung Qu c) Hàm l ng phù sa trung bình tháng t i Kratie, Tơn Chơu, Chơu c, Vàm Nao, C n Th vƠ M Thu n có xu th hƠm l ng phù sa t ng cao tháng mùa l vƠ gi m vào tháng mùa ki t t i t t c tr m xem xét C th , hƠm l ng phù sa l n nh t x y vào tháng VIII (trùng v i th i k có dịng ch y l l n) hƠm l ng phù sa nh nh t xu t hi n vào tháng IV (trùng v i th i k có dịng ch y nh mùa ki t) T ng quan gi a hƠm l ng phù sa l l ng vƠ l u l ng n c đư đ c xây d ng s d ng s li u quan tr c vƠ đo đ c đ c tr ng n m 2011 K t qu th hi n r ng có s liên h ch t ch gi a hƠm l ng phù sa vƠ l u l ng n c, nh t t i v trí nh Kratie T i v trí n m vùng ng b ng Sông 84 C u Long ậ h l u c a dịng sơng Mê Kông, t ng quan gi a hƠm l ng phù sa vƠ l u l ng n c có s phân tán nh t đ nh i u có ngh a lƠ ngoƠi l u l ng dịng ch y hƠm l ng phù sa l l ng n c ph thu c vào y u t khác nh ngu n b sung bùn cát vào sông, nh h ng c a s phân chia dòng ch y nhánh sông thu tri u Tuy nhiên, y u t nêu l i không đ c xem xét m i t ng quan Do đó, đơy c ng lƠ lỦ mƠ t ng quan gi a hƠm l ng phù sa l l ng vƠ l u l ng n c t i v trí vùng ng b ng sơng C u Long có s phân tán D a vào s li u hƠm l ng phù sa thu th p vƠ c tính đ c, b n đ phân b hƠm l ng phù sa trung bình mùa l , mùa ki t c n m đư đ c xây d ng cho n m 2000 vƠ 2011, cho phép xem xét vƠ hình dung s b v s thay đ i c a hƠm l ng phù sa theo không gian vùng nghiên c u B n đ phân b n ng đ bùn cát theo không gian th i gian đ c tính tốn lu n v n m i ch mang tính ch t minh h a cho kh n ng ng d ng c a mô S d ng t ng quan gi a l u l ng n c vƠ hƠm l ng phù sa đư xơy d ng, hàm l ng phù sa t i v trí xem xét vùng nghiên c u (Kratie, Tân Châu, Châu c, Vàm Nao, M Thu n C n Th ) c ng đư đ c c tính cho n m 2011 Ngồi ra, s d ng t ng quan đư xơy d ng hƠm l ng phù sa t i 12 nhánh sông đ vào vùng bi n h Tonle Sap đư đ c xác đ nh s li u đ u vào mô ph ng s d ng mơ hình hai chi u 2D-FEM S d ng đ c tr ng dòng ch y vƠ hƠm l ng phù sa mùa l n m 2000 2011, mô ph ng hƠm l ng phù sa s d ng mơ hình hai chi u 2D-FEM đư đ c th c hi n K t qu mô ph ng hƠm l ng phù sa t mơ hình sau đ c so sánh v i s li u đo đ c c ng nh k t qu c tính hƠm l ng phù sa l l ng t quan h t ng quan, th hi n r ng có s t ng đ ng v s thay đ i hƠm l ng phù sa theo th i gian gi a ph ng pháp VƠo th i gian có l u l ng dịng ch y l n hƠm l ng phù sa tính tốn theo c ph ng pháp đ u cho k t qu l n Ng c l i vào th i k có l u l ng dịng ch y nh k t qu mơ ph ng t mơ hình tốn c ng cho hƠm l ng phù sa nh Ki n ngh Bên c nh m t s k t qu đư đ t đ c, c ng c n ph i nh n m nh r ng khu v c nghiên c u r ng bao g m toàn b vùng ng b ng sơng C u Long, vi c thu th p s li u, d li u nh s li u đ a hình, th y v n, h i v n (thu tri u), đ c tr ng dòng ch y (v n t c, đ sâu hay m c n c, l u l ng n c), bùn cát g p r t nhi u khó khan H u h t s li u thu th p đ c r t r i r c vƠ không đ ng nh t v m t th i gian ó c ng lƠ lỦ k t qu xây d ng t ng quan 85 gi a l u l ng n c vƠ hƠm l ng phù sa c ng nh k t qu mơ ph ng t mơ hình hai chi u 2D-FEM nhi u h n ch M t khác s không đ ng b v d li u s li u thu th p (ví d n m có đ c tr ng dịng ch y nh l u l ng, m c n c l i khơng có s li u v bùn cát ho c tr ng h p ng c l i có s li u v bùn cát l i khơng có s li u v thu tri u vƠ đ c tr ng dòng ch y Do vi c th c hi n mơ ph ng s d ng mơ hình hai chi u 2D-FEM g p r t nhi u khó kh n ơy c ng lƠ lỦ thơng s mơ hình đ c k th a t nghiên c u tr c đơy Mơ hình hai chi u 2D-FEM đ c phát tri n b i TS Ph m Thanh H i đư đ c m nh d n s d ng lu n v n thay mơ hình th ng m i nh MIKE 21 hay mơ hình khác nh TELEMAC Do mơ hình 2D-FEM mơ hình t phát tri n nên tài li u h ng d n v s d ng mơ hình c ng nh tƠi li u liên quan khác ch a nhi u th t s phong phú Do đó, vi c l n đ u tiên ti p c n s d ng mơ hình cho nghiên c u s bi n đ ng c a bùn cát vùng ng b ng sông C u Long lu n v n khơng tránh kh i nh ng khó kh n không ch vi c s d ng mô hình mà c vi c ti p c n hi u v c u trúc c ng nh code c a mơ hình x lý d li u đ u vƠo vƠ đ u c a mơ hình V i di n tích vơ r ng l n nh vùng ng b ng h l u sông Mê Kông, vi c s d ng mơ hình th ng m i ví d nh mơ hình MIKE 21 lƠ vi c g n nh khơng th b i s ph n t hay s ô l i dùng vi c th hi n l i vùng nghiên c u l n Cu i cùng, đ đ m b o đ xác c a k t qu tính tốn mơ ph ng t mơ hình, tài li u v đ a hình, thu v n, h i v n, bùn cát c n ph i đ c thu th p b sung c p nh t h n n a Sau đó, thơng s mơ hình s đ c hi u ch nh c n th n tr c mơ hình đ c ki m đ nh áp d ng cho tính tốn mô ph ng ng v i u ki n đ u vào (t i biên th ng l u vƠ h l u) khác nhau, nh m đ m b o t o c s t t h n cho vi c xác đ nh đ nh l ng chi ti t s bi n đ ng c a bùn cát theo c th i gian không gian vùng nghiên c u 86 TÀI LI U THAM KH O 10 11 12 13 14 15 16 Tr n Th c vƠ Lê Nguyên T ng (2010) Vi t Nam ng phó thích ng v i bi n đ i khí h u T p chí tài ngun mơi tr ng, s 3/2010 Báo cáo “Tác đ ng c a bi n đ i khí h u lên tài nguyên n c bi n pháp thích ng”, 11/2010, Vi n khoa h c khí t ng th y v n vƠ môi tr ng MeKong River Commission (2011) Climate Change Adaptation Initiative, 20112015 Programme Document MeKong River Commission (2011) Annual MeKong Flood Report Báo cáo đ tƠi “Nghiên c u đánh giá tác đ ng c a b c thang th y n dịng h l u sơng Mê Kơng đ n dịng ch y, mơi tr ng, kinh t xã h i vùng ng b ng sông C u Long đ xu t gi i pháp gi m thi u b t l i” ậ KC08-1311 Volume II: Supporting Report, Paper IV: Development of Hydro-Hydraulic Model for the Cambodian Floodplains; WUP-JICA, March 2004 Jansson M.B (1997) Comparison of sediment rating curves developed on load and on concentration Nordic Hydrology, 28(3), 189ậ200 Walling D.E (1977) Assessing the accuracy of suspended sediment rating curves for a small basin Water Resources Research, 13(3), 531ậ538 Cohn T.A (1995) Recent advances in statistical methods for the estimation of sediment and nutrient transport in rivers Reviews of Geophysics, 33(1), 1117ậ 1123 Cohn T.A., Delong L.L., Gilroy E.J., Hirsch R.M., Wells D.K (1989) Estimating constituent loads Water Resources Research, 25(5), 937ậ942 Nguy n T t c (2016), Bàn v m t s cách tính biên cho tốn th y l c c a QHTT-TL BSCL u ki n bi n đ i khí h u n c bi n dâng Nguy n T t c (2005) Mơ hình tốn cho dịng ch y ch t l ng n c h th ng kênh sông Nhà xu t b n Nông Nghi p, 240 Trang Nguy n T t c (2010) Ph n m m DELTA cho tính tốn dịng ch y ch t l ng n c h th ng kênh sông Nhà xu t b n Nông Nghi p, 124 Trang Papanicolaou, A.N., Elhakeem, M., Krallis, G., Prakash, S and Edinger, J (2008) Sediment transport modelling review - current and future developments Journal of Hydraulic Engineering, 134, 1-14 Pham Thanh Hai, T Masumoto, K Shimizu (2008) Development of a twodimensional finite element model for inundation processes in the Tonle Sap and its environs Hydrological Processes 22, 1329-1336 C Pham Van, Pham Thanh Hai (2016) Preliminary results of the flow in the Mekong river-lake-delta system APHydro2016, Hanoi, Vietnam, pp 151-158 87 17 Pham Thanh Hai (2010) Large-scale flooding analysis in the suburbs of Tokyo Metropolis caused by levee breach of the Tone River using a 2D hydrodynamic model Water Science and Technology, 62(8), 1859ậ1864 18 Pham Thanh Hai (2013) ng d ng mơ hình th y đ ng l c h c hai chi u phân tích q trình ng p l t v i quy mô l n gây v đê, v i tác đ ng c a bi n đ i khí h u ậ Tr ng h p sông Tône, Nh t B n T p chí H i C h c Th y khí, Trang 269-276 19 Leclerc, M, Bellemare, JF, Dumas, G, and Dhatt G(1990) A finite element model of estuarian and river flows with moving boundaries Advances in Water Resources, 4, 158ậ168 20 Kawachi T (1987) Finite element analysis of tidal flow in the Nakaumi Shinjiko estuarine basin Journal of the Japan Society of Irrigation, Drainage and Reclamation Engineering, 128, 93ậ101 21 Kawahara, M and Umetsu T(1986) Finite element method for moving boundary problems in river flow.International Journal for numerical methods in Fluids, 6, 365ậ386 22 Kawahara, M, Hirano, H, and Tsubota K (1982) Selective lumping finite element method for shallow water flow.International Journal for Numerical Methods in Fluids, 2, 89ậ112 23 Ariathuru R and Krone R.B (1976) Finite element model for cohesive sediment transport Journal of hydraulic Engineering, 102(3), 323-338 24 Krone R.B (1962) Flume studies of the transport of sediment in estuarial processes Hydraulic Engineering Laboratory and Sanitary 200 References Engineering Research Laboratory, University of California, Berkeley, California 25 Mehta A.J., Partheniades E (1975) An investigation of the depositional properties of flocculated fine sediments Journal of Hydraulic Research, 12(4), 361381 26 Manh N.V, Dung N.V, Hung N.N, Merz B, Apel H (2014) Large-scale suspended sediment transport and sediment deposition in the Mekong Delta Hydrol Earth Syst Sci 18, 3033-3053 27 Partheniades, E (1965) Erosion and deposition of cohesive soils Journal of the Hydraulic Division, 91, 105-139 28 Lang G., Shubert R., Markofsky M., Fanger H U., Grabemann I., Neumann H L., Riethmuller R (1989) Data interpretation and numerical modeling of the mud and suspended sediment experiment Journal of Geophysical Research, 94(C10), 14318-14393 29 C Pham Van, O Gourgue, M.G Sassi, A.F.J Hoitink, E Deleersnijder, S Soares-Frazão (2016) Modelling fine-grained sediment transport in the Mahakam land-sea continuum, Indonesia Journal of Hydro-Environmental Research, 13, 103-120 88 30 Neumeier U., Lucas C H., Collins M (2006) Erodibility and erosion patterns of mudflat sediments investigated using annular flume Aquatic Ecology, 40(4), 543-554 31 Le Hir P., Monbet Y., Orvain F (2007) Sediment erodability in sediment transport modelling: Can we account for biota effects? Continental Shelf Research, 27, 1116-1142 32 C Pham Van, B de Brye, E Deleersnijder, A.J.F Hoitink, M.G Sassi, B Spinewine, H Hidayat, S Soares-Frazão (2016) Simulations of the flow in the Mahakam river-lake-delta system, Indonesia Environmental Fluid Mechanics, 16, 603-633 33 Lu, XX and Siew, RY (2006) Water Discharge and Sediment Flux Changes over the Past Decades in the Lower Mekong River: Possible Impacts of the Chinese Dams Hydrology and Earth System Sciences, 10, 181-195 34 Fu, K D and He, D M.( 2007) Analysis and prediction of sediment trapping efficiencies of the reservoirs in the mainstream of the Lancang River, Chinese Sci Bull., 52, 134ậ140 35 Fu, K D., He, D M., and Lu, X X (2008) Sedimentation in the Manwan reservoir in the Upper Mekong and its downstream impacts, Quatern Int., 186, 91ậ99 36 Kummu, M and Varis, O (2007) Sediment-related impacts due to upstream reservoir trapping, the Lower Mekong River, Geomorphology, 85, 275ậ293 37 Kummu, M., Penny, D., Sarkkula, J., and Koponen, J.( 2008) Sediment: curse or blessing for Tonle Sap Lake?, Ambio, 37, 158ậ163 38 Kummu, M., Lu, X X., Wang, J J., and Varis, O (2010) Basinwide sediment trapping efficiency of emerging reservoirs along the Mekong, Geomorphology, 119, 181-197 39 Walling, D E (2008) The changing sediment load of the Mekong River, Ambio, 37, 150ậ157 40 Gupta, H., Kao, S.-J., and Dai, M (2012) The role of mega dams in reducing sediment fluxes: A case study of large Asian rivers, J Hydrol., 464ậ465 41 Liu, C., He, Y., Walling, E., and Wang, J.( 2013) Changes in the sediment load of the Lancang-Mekong River over the period 1965ậ2003, Sci China Technol Sci., 56, 843ậ852 42 Liu, X and He, D (2012) A new assessment method for comprehensive impact of h dropower development on runoff and sediment changes, J Geogr Sci., 22, 1034ậ1044 43 Kondolf, G M., Rubin, Z K., and Minear, J T (2014) Dams on the Mekong: Cumulative sediment starvation, Water Resour Res., 50, 1ậ12 44 Thuyen, L X., Tran, H N., Tuan, B D., and Bay, N T (2000) Transportation and de osition of fine sediment during flood season in Long Xuyen Quadrangle ậ Technical Report (Vietnamese), Minist Sci.Technol., HCM city, Vietnam, 1ậ70 89 45 Hung, N N., Delgado, J M., Tri, V K., Hung, L M., Merz, B., Bárdossy, A., and Apel, H.( 2012) Floodplain hydrology of the Mekong Delta, Vietnam, Hydrol Process., 26, 674ậ686 46 Hung, N N., Delgado, J M., Güntner, A., Merz, B., Bárdossy, A., and Apel, H (2014a) Sedimentation in the floodplains of the Mekong Delta, Vietnam Part I: suspended sediment dynamics, Hydrol Process., 28, 3132ậ3144 47 Hung, N N., Delgado, J M., Güntner, A., Merz, B., Bárdossy, A., and Apel, H (2014b) Sedimentation in the floodplains of the Mekong Delta, Vietnam Part II: deposition and erosion, Hydrol Process., 28, 3145ậ3160 48 Vien, D M., Binh, V V., Huong, H T., and Guong, V T (2010) Impact of sediment of Mekong river on rice yield and soil fertility, Proceedings of co ference “Sustainable Agricultural Development in the Mekong Delta”, Cantho University, Cantho, 1ậ11 49 Manh, N V., Merz, B., and Apel, H (2013) Sedimentation monitoring in cluding uncertainty analysis in complex floodplains: a case study in the Mekong Delta, Hydrol Earth Syst Sci., 17, 3039ậ3057 50 Dung, N V., Merz, B., Bárdossy, A., Thang, T D., and Apel, H (2011) Multiobjective automatic calibration of hydrodynamic models utilizing inundation maps and gauge data, Hydrol Earth Syst Sci., 15, 1339ậ1354 51 MRCS/WUP-FIN report: Research findings and recommendations, Mekong River Comm., available from: http://wacc.edu.vn/vi/ wp content/uploads/2013/06/wup-fin2_final-reportpart2.pdf (last access: 15 April 2014), 2007 90 ... c a bùn cát l l ng vùng ng b ng sông C u Long M c dù đư có m t s nghiên c u v s bi n đ ng bùn cát l l ng h th ng sông Mê Kông vƠ đánh giá nguyên nhơn gơy s bi n đ ng Tuy nhiên, ph n l n nghiên. .. đánh giá đ c tình hình bi n đ ng c a bùn cát vùng đ ng b ng sơng C u Long Do đó, đ tài lu n v n ? ?Nghiên c u s bi n đ ng c a bùn cát vùng ng b ng sông C u Long? ?? đư đ a l a ch n C ng c n ph i nh... qua vùng T giác Long Xuyên - Hà Tiên Sông Vàm Nao n i li n sông Ti n sông H u, chuy n m t l ng n c l n t sông Ti n qua sông H u lƠm n c sông H u t ng lên kho ng l n Sau sông Vàm Nao, l ng n c sông