LỜI CẢM ƠN Với lòng chân thành nhất, xin chân thành cảm ơn tới: Tiến sỹ Nguyễn Tiền Giang, thầy hướng dẫn khoa học cho luận văn tôi, điều đạt luận văn kiến thức quý báu mà thầy tận tình dẫn thời gian qua Quý thầy cô trường Đại học Khoa học Tự nhiên, Đại học Quốc gia Hà Nội, đặc biệt quý thầy cô Khoa Khí tượng Thủy văn Hải dương học, quý thầy cô Phòng Đào tạo sau Đại học nhiệt tình hướng dẫn giúp đỡ hoàn thành chương trình Cao học luận văn tốt nghiệp Trong luận văn, tác giả kế thừa số liệu bình đồ đáy sông Đề tài “Giải pháp Khoa học Công nghệ chỉnh trị đoạn sông Hồng qua hà Nội”, (KC.08.14/06-10) GS.TS Lương Phương Hậu chủ nhiệm Do vậy, tác giả xin chân thành cảm ơn tới Thầy Gia đình, bạn bè giúp đỡ, ủng hộ sâu sắc thời gian qua Đặc biệt cảm ơn Mẹ tôi, người bên cạnh động viên để vững tâm phấn đấu học tập, nghiên cứu hoàn thành luận văn Lãnh đạo, bạn đồng nghiệp Viện Khoa học Khí tượng Thủy văn Môi trường tạo điều kiện cho trình học tập hoàn thành luận văn tốt nghiệp Trong trình làm luận văn giới hạn thời gian bước đầu nghiên cứu mô hình mã nguồn mở, nghiên cứu động lực bùn cát nên không tránh thiếu sót Vì vậy, tác giả mong cảm thông ý kiến đóng góp quý báu Thầy cô người quan tâm TÁC GIẢ i MỤC LỤC MỞ ĐẦU V CHƯƠNG 1: TỔNG QUAN VIII 1.1 GIỚI THIỆU TỔNG QUAN VỀ MÔ HÌNH TOÁN HAI CHIỀU viii 1.2 GIỚI THIỆU KHU VỰC NGHIÊN CỨU xiii 1.2.1 Vị trí địa lý xiii 1.2.2 Địa hình, địa mạo xiii 1.2.3 Địa chất, thổ nhưỡng xiv 1.2.4 Lớp phủ thực vật xiv 1.2.5 Khí hậu xv 1.2.6 Đặc điểm thủy văn hệ thống sông ngòi xv 1.2.7 Hiện trạng đoạn sông nghiên cứu xx 1.2.7 Diễn biến hình thái đoạn sông nghiên cứu [4] xxi 1.2.8 Diễn biến đáy sông qua thời kỳ [4] .xxvi 1.3 MỘT SỐ NGHIÊN CỨU ĐÃ CÓ xxvii 1.3.1 Các nghiên cứu Viện Khoa học Thuỷ lợi xxvii 1.3.2 Nghiên cứu quy hoạch tăng khả thoát lũ sông Hồng Cục Đê điều .xxviii 1.3.3 Dự án "Hệ thống kè cứng hoá bờ sông Hồng đường đỉnh kè thành phố Hà Nội", Tổng công ty Tư vấn Xây dựng Thuỷ lợi I (HEC-I) lập năm 2002 xxix 1.3.4 Các nghiên cứu Bộ GTVT .xxx 1.3.4 Nhận xét xxxi CHƯƠNG 2: GIỚI THIỆU CƠ SỞ LÝ THUYẾT MÔ HÌNH TREM XXXII 2.1 CÁC PHƯƠNG TRÌNH CƠ BẢN [11,12,13] .xxxii 2.1.1 Phương trình dòng chất lỏng xxxii 2.1.2 Phương trình liên tục bùn cát xxxiii 2.1.3 Các phương trình sức tải cát xxxiv 2.1.4 Biến đổi phương trình tải cát đáy .xxxiv 2.2 LỜI GIẢI SỐ xxxv 2.2.1 Quan niệm gián đoạn hoá thể tích hữu hạn FVM (finite volume method) xxxv 2.2.2 Lời giải ma trận hệ số đường chéo .xxxvii ii 2.2.3 Gián đoạn hoá phương trình Exner .xxxviii 2.3 ĐIỀU KIỆN BIÊN .xxxviii 2.3.1 Điều kiện biên thuỷ lực xxxviii 2.3.2.Điều kiện biên bùn cát .xxxix CHƯƠNG 3: MÔ PHỎNG THỦY ĐỘNG LỰC CHIỀU BẰNG MÔ HÌNH TREM XL 3.1 THIẾT LẬP ĐIỀU KIỆN BIÊN CHO MÔ HÌNH TREM xl 3.1.2 Thiết lập điều kiện biên mô hình MIKE 11 xl 3.1.3.Nhận xét .xlviii 3.2 XÂY DỰNG BIÊN VÀ LƯỚI TÍNH TOÁN CHO ĐOẠN SÔNG xlviii 3.2.1 Xây dựng biên xlviii 3.2.2 Xây dựng lưới tính toán xlix 3.3 MÔ PHỎNG TRƯỜNG VẬN TỐC CHO TRẬN LŨ THÁNG NĂM 1996 .li 3.4 ĐÁNH GIÁ ẢNH HƯỞNG CỦA CÁC PHƯƠNG ÁN CHỈNH TRỊ lx 3.4.1 Các phương án chỉnh trị [4] lx 3.4.2 Đánh giá ảnh hưởng phương án chỉnh trị đến trường vận tốc lxi 3.4.3 Đánh giá ảnh hưởng phương án chỉnh trị đến khả bồi xói lxvii LXXI LXXII KẾT LUẬN VÀ KIẾN NGHỊ lxxv TÀI LIỆU THAM KHẢO LXXVII iii DANH MỤC HÌNH DANH MỤC BẢNG iv MỞ ĐẦU I TÍNH CẤP THIẾT CỦA ĐỀ TÀI Xói lở bồi tụ trình tất yếu dòng chảy tự nhiên Theo góc độ nhìn nhận khác mà trình gây lợi, hại có tác động lớn đối phát triển kinh tế - xã hội Nếu nắm bắt quy luật dự báo có giải pháp tích cực phục vụ cho phát triển kinh tế - xã hội Trên hệ thông sông Hồng, tượng xói lở xảy thời gian gần gây thiệt hại không nhỏ làm đất nông nghiệp, đất nhà ở… Về tính chất xói lở theo phương thẳng đứng phương nằm ngang Có thể thấy nhiều nơi, đoạn từ Thượng Cát tới phà Khuyến Lương, xói lở theo hướng ngang gây thiệt hại lớn hai bên bờ sông nơi tập trung điểm dân cư sở kinh tế Đây nơi tương lai có dự án thành phố hai bên bờ sông Hồng với tổng chiều dài khoảng 40 km với khu đô thị đường giao thông Chính vấn đề kiểm xoát đưa giải pháp chỉnh trị khả thi có vai trò đặc biệt quan trọng khu vực Ở nước ta, xói lở nghiên cứu từ lâu phương pháp khác nhau, nhiên phương pháp mô hình toán đánh giá cao giải nhiều vấn đề thực tế đặt Lớp mô hình sử dụng bao gồm lớp mô hình chiều, chiều, chiều Các mô hình chiều kể đến HEC-6 (phiên tích hợp HEC-RAS), MIKE 11, ISiS…; số mô hình 2, chiều MIKE 21, EFDC, DEL3D… Các mô hình giải nhiều vấn đề thực tế đặt xây dựng giao điện cho người sử dụng Tuy nhiên mô hình thương mại giới hạn mô hình chiều, người sử dụng sử dụng phần mềm mà tham gia vào thay đổi chương trình phục vụ công tác nghiên cứu ứng dụng Mô hình diễn biến lòng dẫn chiều hệ tọa độ phi trực giao TREM (Twodimensional River bed Evolution Model) mô hình mã nguồn mở, cho phép phân tích tính toán xói lở lòng sông theo chiều, phản ánh gần đầy đủ nguyên nhân v trạng xói lở lòng sông Mô hình TS Nguyễn Tiền Giang phát triển ứng dụng cho sông Việt Nam Với tảng vậy, việc nghiên cứu ứng dụng mô hình tính toán diễn biến lòng sông cho khu vực nghiên cứu tạo điều kiện trau dồi kiến thức, cách tiếp cận cho Học viên tảng quan trọng cho Học viên trình nghiên cứu sau vấn đề bùn cát, mô hình toán mã nguồn mở II ĐỐI TƯỢNG VÀ PHẠM VI NGHIÊN CỨU Luận văn tập trung vào đối tượng phạm vi nghiên cứu sau: Nghiên cứu mô hình mã nguồn mở TREM, khả ứng dụng mô hình cho sông Hồng Phạm vi nghiên cứu ban đầu luận văn tập trung đánh giá diễn biến lòng dẫn cho đoạn sông Hồng từ cầu Chèm đến Khuyến Lương ứng với kịch chỉnh trị lòng sông Tuy nhiên trình nghiên cứu, cho phép Thầy hướng dẫn khoa học điều kiện khách quan phạm vi nghiên cứu thu hẹp lại cho đoạn sông Hồng từ Cầu Long Biên đến Khuyến Lương Một số lý thay đổi sau: - Hiện tại, đoạn sông nghiên cứu có trạm thủy văn Hà Nội nằm hạ lưu cầu Long Biên, trạm thủy văn cấp I nên có số liệu đo đạc tốc độ, lưu lượng, mực nước phục vụ cho việc hiệu chỉnh kiểm nghiệm mô hình Trong đoạn sông từ Chèm đến ngã ba Thượng Cát số liệu thực đo - Ba bốn đoạn sông phương án chỉnh trị tập trung khu vực hạ lưu cầu Long Biên - Thời gian cho mô dòng chảy thời gian thực ½ máy tính, thời gian mô phong bùn cát với tỷ lệ 1/1 nhiều thời gian để mô chưa xét đến thời gian lấy thể kết tính toán III PHƯƠNG PHÁP NGHIÊN CỨU vi Trong luận văn học viên áp dụng phương pháp kế thừa, phương pháp điều tra khảo sát thực địa phương pháp mô hình toán IV NỘI DUNG LUẬN VĂN Ngoài phần mở đầu, kết luận, kiến nghị tài liệu tham khảo, luận văn bao gồm chương: Chương I: Tác giả giới thiệu tổng quan số mô hình toán mô diễn biến lòng dẫn hai chiều, đồng thời giới thiệu sơ lược khu vực nghiên cứu số công trình nghiên cứu có liên quan Chương II: Tập trung giới thiệu sở lý thuyết mô hình TREM xây dựng phát triển TS Nguyễn Tiền Giang Chương III: Trình bày thành nghiên cứu học viên thời gian thực luận văn tốt nghiệp vii CHƯƠNG 1: TỔNG QUAN 1.1 GIỚI THIỆU TỔNG QUAN VỀ MÔ HÌNH TOÁN HAI CHIỀU Việc nghiên cứu xói lở lòng sông tiến hành nhiều nơi giới Nhiều phương pháp mô hình biến dạng lòng sông xây dựng, góp phần giải toán thực tế đặt Đã có tương đối nhiều nghiên cứu tính toán biến hình lòng sông biến động đường bờ, xây dựng công trình chống xói lở bờ sông bờ biển Trên sở hệ phương trình động lực dòng chảy cân bùn cát người ta xây dựng sơ đồ sai phân để tính toán diễn biến lòng sông Velikanov Grixanhin (Nga,1965) đưa phương pháp biến hình giúp cho việc giải toán xói lở dễ dàng xác Sau loạt mô hình đời, dùng riêng cho diễn biến lòng sông, dùng chung cho xói lở bờ sông bờ biển Đó mô hình GENESIS Trung tâm nghiên cứu công nghệ bờ biển Hải quân Mỹ(1989), UNIBEST Viện thuỷ công Hà Lan(1990), LITPACK Viện Thuỷ lực Đan Mạch(1991), SAND94 Viện Thuỷ công Ba Lan(1994), mô hình chiều biến dạng đáy sông Phòng thí nghiệm thuỷ lực Trường Đại học kỹ thuật Deft Hà Lan mô hình miêu tả hiệu ứng xói ngang MIKE21, mô hình SMS Anh Ở nước, Viện Khoa học thuỷ lợi nghiên cứu tình hình xói lở sông Hồng sông Đà sau đập thuỷ điện Hoà Bình, khảo sát nghiên cứu tình hình xói lở sông suối miền Trung biên giới Viện học có nghiên cứu xói lở bờ biển Việt Nam Một loạt dự án chống xói lở nghiên cứu triển khai xây dựng khu vực trọng điểm xói lở thượng lưu sông Hồng, sông Lô, khu vực bờ sông Hồng ven Thủ đô Hà Nội, số điểm xói lở sông suối miền Trung sông Cả, sông Thu bồn, sông Kinh Dinh v.v Các công trình dựa tính toán khoa học đề xuất giải pháp hữu hiệu chống xói lở Tuy nhiên công trình giải cho điểm cụ thể, chưa sâu vào việc sử dụng mô hình cho diễn biến lòng sông Trong thời gian gần sử dụng số mô viii HEC-6, MIKE11 để phân tích, tính toán xói lở Tuy nhiên mô hình giải toán chiều, xem xét xói đáy, với giả thiết chiều sâu xói lở toàn mặt ngang, chưa xem xét vận chuyển bùn cát xói lở không theo chiều ngang Trong thực tế sông thường bị xói lở theo chiều dọc chiều ngang mạnh mẽ chúng có tác động tương hỗ với Rất nhiều nơi xói lở theo hướng ngang có vai trò quan trọng nhất, hai bờ sông nơi tập trung điểm dân cư sở kinh tế Trên hệ thống sông Hồng tình hình xói lở diễn nghiêm trọng, sau hồ chứa Hoà Bình vận hành trình xói lở lại nghiêm trọng Sông suối biên giới phía Bắc bị xói lở nghiêm trọng, ảnh hưởng đến đường biên giới chung nước Một số kết nghiên cứu dự án trước làm sáng tỏ nguyên nhân giải pháp cần thiết để giảm thiểu ảnh hưởng Tuy nhiên việc xói lở diễn mạnh bờ, mà nghiên cứu mô hình trước không giải Vì cần thiết có mô hình chiều để giải toán Gần đây, Việt nam mô hình chiều sử dụng số công trình nghiên cứu biến dạng lòng dẫn MIKE 21C Viện Khoa học Thuỷ lợi sử dụng nghiên cứu lòng dẫn hệ thống sông Hồng (2006-2007) trước sau có hoạt động hồ chứa Hoà Bình Tuyên Quang MIKE 21C áp dụng nghiên cứu biến dạng vùng Cửa sông Hậu (2006) mô lũ lòng sông Cửu long (2005) Trường Đại học Khoa học Tự nhiên Viện Khoa học Thuỷ lợi áp dụng mô hình TREM SMS (2006) nghiên cứu biến dạng đáy số sông biên giới phía Bắc Công ty tư vấn công trình giao thông (TEDY) sử dụng MIKE 21C phân tích trường dòng chảy cửa sông phân lạch khu vực Trung Hà (sông Đà) (2007) tác dụng công trình chỉnh trị Các kết thu cho thấy hiệu mô hình chiều phân tích mô thuỷ lực lòng dẫn đoạn sông Nguyễn Thế Hùng (2004) nghiên cứu giải hệ phương trình sóng nước nông chiều phương pháp phần tử hữu hạn hai giai đoạn theo thuật giải phân rã giải theo thời gian với sơ đồ sai phân Crank-Nicolson Lương Tuấn Anh Nguyễn Hoàng Minh (2006) sử dụng thuật toán Runge- Kutta bậc phương pháp nội suy tuyến tính nối ix tiếp để giải hệ phương trình sóng nước nông chiều Từ thu biểu đồ phân bố vận tốc, độ sâu (mực nước) khu vực nghiên cứu Lương Tuấn Anh Nguyễn Thanh Sơn (2007) áp dụng phương pháp phần tử hữu hạn cho toán thuỷ lực sườn dốc tựa chiều, giải sườn dốc chạy từ đường phân nước tới lòng sông song song với Trên giải sườn dốc giải hệ phương trình chiều Tuy nhiên nghiên cứu dạng giải mặt thuỷ lực chưa đề cập đến biến dạng lòng dẫn Một cách chặt chẽ, hầu hết dòng chảy thực tế sông ngòi dòng chảy với đặc trưng thủy lực biến đổi theo chiều Tuy nhiên, độ rộng dòng chảy lớn nhiều so với độ sâu biến đổi đặc trưng theo chiều thẳng đứng trở nên không đáng kể có xáo trộn mạnh mẽ gây nên ứng suất ma sát đáy Trong trường hợp đó, mô hình dòng chảy chiều (2D) sử dụng hệ phương trình nước nông (shallow water equations) hay hệ phương trình trung bình theo độ sâu (depth-averaged equations) xem đủ để mô tả đặc trưng dòng chảy Hệ phương trình depth-averaged lần Kuipers Vreugdenhill giới thiệu năm 1973 viết hệ tọa độ Đề Tuy nhiên, áp dụng hệ phương trình để tính toán dòng chảy sông ngòi tự nhiên, khó khăn xuất sử dụng điều kiện biên cho đoạn sông cong có biên phức tạp Ban đầu giải pháp sử dụng xấp xỉ "stair stepped" thử nghiệm (Vreugdenhill Wijbenga, 1982), qua ứng dụng, bộc lộ hạn chế độ thô lưới tính gần biên, tăng thời gian tính toán nhớ máy tính muốn tăng độ phân giải lưới Thompson (1980) người giới thiệu hệ tọa độ cong phù hợp biên để giải hạn chế Theo đó, trước hết hệ phương trình depthaveraged viết hệ tọa độ Đề biến đổi sang hệ tọa độ với trục tọa độ ( ξ ,η ) Giải hệ phương trình lưới cong xây dựng cho phép ta có tranh chi tiết đặc trưng dòng chảy sông, đặc biệt đoạn sông cong có biên phức tạp x Hình 3.29: Kết tính toán trường vận tốc lớn trận lũ tháng 8/2000 – PA1 lxiv Hình 3.30 Đáy sông phân bố vận tốc theo phương ngang mặt cắt 46 Hình 3.31 Đáy sông phân bố vận tốc theo phương ngang mặt cắt 74 lxv Hình 3.32 Đáy sông phân bố vận tốc theo phương ngang mặt cắt 111 Hình 3.33 Đáy sông phân bố vận tốc theo phương ngang mặt cắt 146 lxvi Hình 3.34 Kết lưu tốc Phương án mép nước bờ tả, hữu dòng 3.4.3 Đánh giá ảnh hưởng phương án chỉnh trị đến khả bồi xói Để tiến hành đánh giá ảnh hưởng phương án chỉnh trị đến khả bồi xói cần có lưu lượng bùn cát biên thượng lưu Trong luận văn trình bày kết mô mạng sông Hồng từ Hoà Bình, Yên Bái, Vụ Quang đến Hưng Yên mô hình MIKE 11 (mô đun AD NCST) Kết cho quan hệ lưu lượng bùn cát lưu lượng nước câu Long Biên hình 3.39 lxvii Hình 3.35.Quan hệ Qs~Q năm 2000 2001 từ kết mô MIKE 11 [5] Chú ý thiết lập quan hệ tương quan lưu lượng nước lưu lượng bùn cát cần tiến hành chuyển đổi đơn vị sau: Qs=bQa Trong đó: Qs lưu lượng bùn cát tổng cộng (tấn/ngày); Q lưu lượng nước (m3/s); a,b hệ số tương quan Đối với trận lũ năm 2000, sau trình tính toán xác định hệ số a lấy 2,7 b lấy 0,00002 (sau trừ 33% bùn cát đáy) Trong điều kiện sử dụng lưu lượng tạo lòng để mô lưu lượng bùn cát lấy theo lưu lượng bùn cát mô theo quan hệ hình 3.39 Đường kính hạt lấy D50=0,1mm, D90=0.2mm Với điều kiện sử dụng lưu lượng thực đo năm 2000 ( PA 0) lưu lượng tạo lòng để mô 7650m3/s [4] (đối với PA 1) Mô hình chạy ổn định với bước thời gian mô bùn cát 0.5 giây Kết sau 126 mô (trong có h đầu tính thủy lực để mô hình ổn định) cho số nhận xét sau: lxviii Phương án trạng (PA0): Trong điều kiện trạng, tượng xói xuất chủ yếu chiều sâu xói lớn sau 120 mô đạt 1.0m Tại khu vực có dòng chủ lưu áp sát bờ tượng xói xuất với độ sâu xói lớn đạt 0.2m đến 0.4m/ngày Đồng thời, điều kiện trạng, tượng xói khu vực bãi xảy ra, nhiên mức độ lớn 0.01m sau 120 mô Hiện tượng bồi xảy số vị trí thuộc dòng chủ lưu nhiên mức độ không lớn mang tính chất cục Như qua mô thủy động lực thấy mô hình mô quy luật cho thấy tranh tổng thể trình bồi/xói điều kiện trạng Phương án chỉnh trị (PA1): Các điều kiện thủy lực bùn cát trạng sử dụng phương án trừ địa hình đáy sông cập nhật Để kết diễn biến lòng sông mô cách xác hơn, với điều kiện địa hình đáy thay đổi theo phương án chỉnh trị cho phép thu hẹp miền tính theo phương ngang nhằm giảm nhỏ kích thước ô lưới theo phương Lưới thể hình 3.40 Kết mô cho thấy phạm vi mức độ bồi/xói khu vực nghiên cứu Từ kết phân tích hình 3.42 đến 3.47 cho thấy: Hiện tượng bồi xói xuất xen kẽ với tốc độ bòi xói nhỏ Nếu bỏ 10 mặt cắt miền tính để loại ảnh hưởng điều kiện biên tốc độ bồi lớn khoảng 0.1 m/ngày, tốc độ xói lớn 0.05m/ngày từ mặt cắt 100 đến 110 Do có thay đổi dòng chủ lưu, địa hình nên khu vực bị bồi/xói thay đổi rõ rệt Nếu điều kiện trang thấy hầu hết khu vực xuất tình trạng xói điều kiện chỉnh trị tượng xói xuất số vị trí mức độ không lớn lxix Đồng thời khu vực xói không áp sát bờ, điều khẳng định phương án chỉnh trị phát huy hiệu việc giảm thiểu nguy xói đỉnh cong Hình 3.36 Lưới tính toán diễn biến lòng dẫn cho PA lxx Hình 3.37 Bồi xói sau 120 mô (lưu lượng thực đo) - PA0 lxxi Hình 3.38 Bồi xói 120 mô phỏng, phương án (lưu lượng tạo lòng) Viền biên miền tính cũ, gần trùng với bờ đê TW trạng lxxii Hình 3.39 Bồi xói 120 mô phỏng, phương án (Q tạo lòng) – mc 46 Hình 3.40 Bồi xói 120 mô phỏng, phương án (Q tạo lòng) – mc 74 lxxiii Hình 3.41 Bồi xói 120 mô phỏng, phương án (Q tạo lòng) – mc 111 Hình 3.42 Bồi xói 126 mô phỏng, phương án (Q tạo lòng) – mc 146 Hình 3.43 Thay đổi đường lạch sâu dọc đoạn sông nghiên cứu (Zsau –Z trước) sau 120 mô với lưu lượng tạo lòng, phương án lxxiv KẾT LUẬN VÀ KIẾN NGHỊ Kết luận Đã có nhiều nghiên cứu bồi/xói sông Hồng đặc biệt thời gian gần chủ đề thành phố hai bên bờ sông Hông nên việc sử dụng phương pháp khác để đánh giá bồi/xói sở khoa học quan trọng để đưa giải pháp ổn định lòng dẫn, phát triển kinh tế - xã hội Đồng thời việc sử dụng mô hình toán mô diễn biến lòng dẫn hai chiều luận văn sơ đánh giá định lượng hiệu phương án chỉnh trị đoạn sông nghiên cứu Luận văn trình bày tóm tắt kết ứng dụng mô hình TREM để mô trường tốc độ, bồi/xói cho khu vực nghiên cứu Qua rút số kết luận sau: 1- Mô hình MIKE 11 mô tốt trình mực nước, lưu lượng cho toàn mạng sông nghiên cứu thể qua kết tính toán mực nước, lưu lượng Đồng thời kế thừa kết nghiên cứu bùn cát mô hình chiều Đây sở quan trọng để khôi phục biên cho mô hình TREM số liệu thực đo 2- Bằng việc chia lưới tính hợp lý với trình hiệu chỉnh thông số mô hình kỳ công tốn nhiều thời gian kết tính toán cho thấy mô hình TREM có khả mô tốt trường tốc độ phân bố tốc độ theo phương ngang, thể qua kết tính toán thực đo tốc độ theo thủy trực trạm thủy văn Hà Nội 3- Đối với trường tốc độ toàn miền tính, hình dạng đoạn sông nghiên cứu có tác động mấu chốt đến phân bố trường tốc độ mô hình mô tốt cho trục động lực khu vực bãi bồi doi cát 4- Qua mô diễn biến lòng dẫn phương án địa hình trạng cho thấy kết mô phản ánh quy luật bồi/xói khu vực nghiên cứu lxxv 5- Kết mô phương án chỉnh trị cho thấy chế độ động lực dòng chảy, diễn biến bồi xói ổn định, có lợi mặt ổn định lòng sông lâu dài Với phương án dành nhiều quỹ đất cho phát triển kinh tế xã hội cho hai bên bờ sông 6- Với ý nghĩa đặc biệt đoạn sông, nơi có hàng loạt công trình kè, mỏ hàn, cảng, phà, nơi có vị trí quan trọng phát triển kinh tế, đô thị, kết nghiên cứu góp phần bước đánh giá định lượng chế độ thủy động lực khả ứng dụng mô hình toán nghiên cứu Kiến nghị 1- Luận văn tổng hợp, thu thập chỉnh lý hệ thống hóa khối lượng lớn số liệu, tài liệu… làm sở tính toán đặc trưng thủy động lực tài liệu tham khảo cho nghiên cứu khác 2- Kết nghiên cứu thể luận văn dựa mô hình toán để định lượng trình thủy động lực điều kiện khác sử dụng kết cho công tác chỉnh trị nghiên cứu khác liên quan 3- Với kết luận văn nghiên cứu trước mô hình TREM cho thấy mô hình có khả mô tốt trình thủy động lực bùn cát, đặc biệt cho đoạn sông cong theo tác giả sử dụng mô hình cho sông khác Việt Nam 4- Là mô hình mã nguồn mở cần tiếp tục nghiên cứu hoàn thiện mô hình đáp ứng yêu cầu thực tế độ xác, thời gian mô tính thân thiện với người sử dụng 5- Với thời gian có hạn với bước đầu nghiên cứu mô hình mã nguồn mở, nghiên cứu bùn cát nên thân học viên thấy cần phải đầu tư nhiều thời gian công sức để có kết tốt thời gian tới hạn chế luận văn lxxvi TÀI LIỆU THAM KHẢO Tài liệu tiếng việt Vũ Thanh Ca, Nguyễn Quốc Trinh, 2008 Nghiên cứu nguyên nhân xói lở bờ biển Nam Định Nguyễn Kiên Dũng, 2001, Nghiên cứu xây dựng sở khoa học tính toán bồi lắng cát bùn hồ chứa Hòa Bình, Sơn La Luận án Tiến sĩ Địa lý Nguyễn Tiền Giang Ngô Thanh Nga (2010) Đánh giá ảnh hưởng phương án chỉnh trị đến khả thoát lũ đoạn sông Hồng chảy qua Hà Nội cũ mô hình mô Tạp chí Khoa học Đại học Quốc gia Hà Nội, Khoa học Tự nhiên Công nghệ 26, Số 3S, 322-332 Lương Phương Hậu, 2010 Giải pháp Khoa học Công nghệ chỉnh trị đoạn sông Hồng qua Hà Nội Báo cáo tổng kết đề tài cấp Nhà nước Chương trình KC.08.14/06-10 Lương Văn Thanh cs., 2004 Đánh giá ảnh hưởng chất độc hoá học môi trường hồ Trị An - Đề xuất giải pháp khắc phục Báo cáo tổng kết đề tài khoa học cấp Nhà nước Chương trình 33 Trần Đức Thịnh (2009), Khai thác mô đun vận chuyển bùn cát mô hình MIKE 11 chạy thử nghiệm cho hệ thống sông Hồng (đến trạm hà Nội) Khoá luận tốt nghiệp Cử nhân Khoa học Trường Đại học Khoa học Tự nhiên, ĐHQG Hà Nội Trần Anh Tú & NNK, 2009 Một số kết nghiên cứu xói lở, bồi tụ vùng ven bờ khu vực Hải Phòng Tài liệu tiếng Anh Cheng, N.S., 1997, Simplified settling velocity formula for sediment particle J Hydr Engrg.,ASCE, 123(2):149-152 DHI MIKE, 2007, MIKE 21 & MIKE Flow model FM - Hydrodynamic and Transport Module Scientific Documentation lxxvii 10 DHI MIKE, 2007, MIKE 11 - Hydrodynamic and Transport Module Scientific Documentation 11 Nguyen Tien Giang, 2001, Two-dimentional Riverbed Evelution Model (TREM) constructed in the non-orthogonal curvilinear coordinate system Project Report Pp52, Japanese Government Project on Regional Education and Development Project AIT, Thailand 12 N.T Giang and N Izumi (2001) Application of an Integrated Morphological Model to Red River Network and Son Tay Curved Bend, Vietnam, (Ed S Ikeda), pp 295-304, 2001.The proceedings of International congress on River and Coastal Morphology (RCEM), September 2001 Hokaido, Japan 13 Nguyen Duy Khang, 2001 Two-Dimensional free-surface flow model based on boundary-fitted coordinate system Thesis of Master.AIT, Thailand 14 Gessler, J., 1967 The beginning of bedload movement of mixtures investigated as natural armoring in channels W.M Keck Laboratory of Hydraulics and Water Resources, California Institute of Technology, Translation T-5 15 Krone, R B., 1962 Flume studies of the transport of sediment in estuarial shoaling processes Final Report to San Francisco District, U.S Army Corps of Engineers, Washington, D.C 16 John M Hamrick, 1996 User's manual for the environmental fluid dynamics computer code Updated & Revised by Paul M Craig-Dynamic Solutions, LLC Web page: www.ds-intl.biz 17 Van Niekerk, A., Vogel, K., and Slingerland, R and Bridge, J., 1992 Routing of heterogeneous sediments over movable bed: model development J Hydr Engr., 118(2): 246-263 18 Van Rijn, L.C., 1984a Sediment transport, Part I: Bed load transport J Hydr Engr., 110(10): 1431-1456 19 Van Rijn, L.C., 1984b Sediment transport, Part II: Suspended load transport J Hydr Engrg., ll6(ll): 1612-1638 lxxviii [...]... ca bựn cỏt l lng da trờn lý thuyt ca Galappatty (1983), ch xột khuch tỏn theo chiu thng ng, v cú dng tng t nh trong mụ hỡnh TREM Gii c tin hnh cho tng nỳt li, gii ng thi cho c thu lc cht lng v hỡnh thỏi lũng dn, tng t mụ hỡnh TREM Li c xỏc nh theo h to cong cho phộp gii trờn sai phõn hu hn v phn t hu hn, cú mt hng song song vi hng dũng chy, cũn mt hng gn nh trc giao vi nú Li khụng c nh, sau khi xúi... phn thu lc cht lng v SED2D gii cho phn bựn cỏt v bin dng lũng dn H phng trỡnh s dng cng gm 2 phng trỡnh mụmen cho 2 chiu x v y, mt phng trỡnh liờn tc cht lng v mt phng trỡnh liờn tc ca bựn cỏt Cỏc phng trỡnh thu lc c gii bng phng phỏp phn t hu hn vi phng phỏp s d trng s Galerkin Gii bng s n y v phng phỏp lp Newton-Raphson Cỏc s hng ngun v ra tr thnh cỏc s hng ngun v ra cho mụ hỡnh ỏy tớnh toỏn bin... cng H Ni sau l 1984 ó ct c mt on uụi khỏ di so vi nm 1983 v lm cho dung tớch lũng sụng vựng cng di cao trỡnh +5,0 t 1,4 triu m3 nm 1982 tng lờn 5,9 triu m3 nm 1984 S phỏt trin nhanh chúng ca Bói T Liờn t sau nm 1984 lm cho phng hng dũng chy thay i phớa thng lu ca ung uụi bói T Liờn phỏt trin mnh v ngy cng thu hp ca vo nhỏnh Phỳc Xỏ lm cho dũng chy i vo nhỏnh ny gp khú khn hn Nh vy, nhỏnh Gia Lõm cú... nờn lng cung cp cho sụng ớt v vỡ vy dũng chy sụng Lụ mang rt ớt bựn cỏt Vựng thuc dóy nỳi Phan-Xi-Pan cú dn tớch rng, cao v a hỡnh cú s thay i ln, khng ch nhng vựng khớ hu, th nhng rt khỏc nhau t vựng ny c phỏt trin t cỏc loi ỏ gc nh dip thch tinh th, hoa cng, cng xung phớa Tõy Nam dip thch v ỏ vụi cng nhiu cũn phớa ụng Nam l dip thch v hoa cng õy l khu vc cung cp bựn cỏt quan trng cho sụng , gúp... cho mụ hỡnh ỏy tớnh toỏn bin i ỏy ti tng nỳt vi vic dựng trng s Crank-Nicholson ỏy cỏt c coi l phự hp vi mt b bựn cỏt dy hu hn, di nú l mt mt khụng xúi Bựn cỏt c thờm vo hoc ly i khi ỏy c xỏc nh bi s hng ngun v ra bc thi gian trc ỏy bựn-sột c x lý nh mt dóy cỏc lp, mi lp cú cỏc c tớnh ca nú v dy, mt , tui v cng xúi, ng sut tip cho xúi, h s hp nht Bựn-sột mi s bi lờn lp c v lm tng ỏp sut nộn v tui... 1.10 th hin Chỳ ý rng, tuyn ng thu vch cho lũng dn mc nc +6,0m v tuyn chnh tr vn gi trng thỏi phõn lch vựng Tm Xỏ v vựng T Liờn - Trung H Cỏc bin phỏp cụng trỡnh c phõn thnh 2 giai on hon thin dn: giai on 2010 v giai on 2020 3- Hin nay TEDI ang c B GTVT giao nhim v lp d ỏn NCKT v "Ci to giao thụng thu sụng Hng khu vc H Ni", ch yu ly kt qu nghiờn cu ca JICA lm cn c cho vic lp quy hoch tuyn lung v thit... cụng trỡnh n nh lung lch Nhng TEDI vch tuyn chnh tr cho lũng dn mc nc +9,0m xxx 1.3.4 Nhn xột a) Nhng vn ó c khng nh Lũng dn sụng Hng khu vc H Ni ang b nhiu tỏc ng bt li lm suy gim kh nng thoỏt l, nht l s khai thỏc vụ k hoch vựng bói Din bin on sụng nghiờn cu ch yu l dch chuyn qua li ca lch chớnh v lch ph trong on sụng phõn lch, th sụng thun li nht cho vic tho món cỏc yờu cu khai thỏc v phỏt trin kinh... S 2.2.1 Quan nim ca giỏn on hoỏ trong th tớch hu hn FVM (finite volume method) C s ca phng phỏp th tớch hu hn (FVM) cn c trờn nguyờn tc bo ton khi nc cho mt th tớch hu hn Nú thng dựng dng tớch phõn ca phng trỡnh cõn bng nh im khi u Phng trỡnh cõn bng cho mt lng vn chuyn vụ hng bi dũng chy cú dng tớch phõn: d + VndS = gradndS + q d t S (1) (2) (3) ( 4) xxxv (2.22) Trong phng trỡnh (2.22), v... Tõy Cụn Lnh v vựng nỳi Ba Vỡ Tam o cú lng ma t 2400 mm/nm n 2800 mm/nm c bit tõm ma Bc Quang cú nm t n 5499 mm/nm Nhng vựng khut sau nhng dóy nỳi chn giú nh thung lng Yờn Chõu, cao nguyờn Sn La, lũng cho Ngha L, vựng thng ngun sụng Gõm cú lng ma nh, ch khong 1200 mm/nm n 1600 mm/nm Vựng ng bng cú lng ma trung bỡnh 1700 mm/nm Lng ma bin i qua cỏc nm trong thi k nhiu nm khụng ln: lng ma ca nm ma nhiu... Nam-Trung Quc cao trờn 2000 m Sụng Thao l in hỡnh v hng ca mt con sụng do vn ng to sn Himalaya vch ra Cú th núi sụng Thao cú hng chy khỏ n nh: tr mt on ngn u ngun, on cũn li khỏ thng theo hng Tõy Bc-ụng Nam cho ti Vit Trỡ v ca sụng Tng lng nc bỡnh quõn nhiu nm ca sụng Thao ti Vit Trỡ l 28.4 km 3 tng ng vi lu lng bỡnh quõn l 500 m 3/s v mụ un dũng chy nm l 17.31 l/s.km2 xvi Ch dũng chy trờn sụng Thao ph thuc