ỨNG DỤNG MÔ HÌNH MIKE11 NGHIÊN CỨU CHẾ ĐỘ THUỶ VĂN, THUỶ LỰC NGÃ BA SÔNG QUẢNG HUẾ - LƯU VỰC VŨ GIA - THU BỒN Sinh viên thực hiện: Vũ Minh Đức, 47B Giáo viên hướng dẫn: PGS TS Vũ Minh Cát Th.S Phạm Quang Chiến Tóm tắt: Hệ thống sông Vũ Gia - Thu Bồn thuộc tỉnh Quảng Nam với diện tích lưu vực tính tới cửa 10.350 km2 Hệ thống gồm sông Vũ Gia phía bắc Thu Bồn phía nam Hai sông có liên hệ thuỷ lực qua khu vực sông Quảng Huế, nơi dòng chảy tập trung lòng dẫn mùa kiệt chảy tràn qua bãi chủ yếu từ sông Vũ Gia sang sông Thu Bồn vào mùa lũ Khu vực hạ lưu vùng đồng thành phố Đà Nẵng, Hội An, nơi có mật độ dân số cao trung tâm kinh tế xã hội vùng nên nhu cầu nguồn nước mùa kiệt lớn, mùa lũ lại cần thiết chia nước cho nhánh Thu bồn, nhằm giảm nhẹ tình hình lũ lụt cho đồng Do chế độ thuỷ văn không đồng với mùa lũ ngắn tháng từ tháng IX đến tháng XII tổng lượng dòng chảy 80% năm làm khu vực ngã ba Quảng Huế thường xuyên nước tràn bờ, xói bãi tạo lòng mùa kiệt đến, để tự nhiên nước sông Vũ Gia chuyển hết sang Thu bồn gây tình trạng thiếu nước, xâm nhập mặn nghiêm trọng cho huyện phía bắc tỉnh thành phố Đà Nẵng Báo cáo trình bày kết nghiên cứu từ việc ứng dụng mô hình MIKE11 để tính toán đánh giá tác động dòng chảy mùa lũ tới khu vực nghiên cứu Kết thu từ MIKE11 sử dụng đầu vào mô hình chiều MIKE21 tính toán thông số để có tranh chi tiết diễn biến khu vực có lũ đề xuất giải pháp nhằm ổn định lâu dài, giảm thiểu xói lở mùa lũ kiểm soát, điều chỉnh dòng chảy mùa kiệt đáp ứng tỉ lệ phân chia nước phục Vũ phát triển kinh tế xã hội bền vững cho lưu vực sông Đặt vấn đề Sông Quảng Huế sông nhánh thuộc đoạn cuối sông Vũ Gia nối sang sông Thu Bồn Trong thời kỳ mùa kiệt, nước thường chảy lòng dẫn với lưu lượng thay đổi lòng sông bị bồi lắng Tuy nhiên, thời kỳ mùa lũ, nước sông thường tràn bờ tạo vùng giao thoa rộng lớn khu vực này, gây xói lở bãi lòng sông Quảng Huế, tạo phân bố lưu lượng không mong đợi mùa kiệt cho phần hạ lưu Vũ Gia - Thu Bồn, trực tiếp ảnh hưởng đến đời sống dân sinh kinh tế quanh khu vực sông Quảng Huế Trong báo này, sử dụng mô hình chiều MIKE11 để mô hệ thống sông Thu Bồn – Vũ Gia, đưa kết tính toán lũ đề xuất giải pháp công trình phù hợp Kết từ mô hình MỊKE11 số liệu đầu vào cho mô hình MIKE21 dùng để tính toán chi tiết phân bố không thời gian trường lưu lượng, lưu tốc nhằm đánh giá chi tiết tình hình xói bồi đề xuất giải pháp chỉnh trị ổn định lâu dài đoạn ngã ba sông Quảng Huế Tuy nhiên, hạn chế thời gian tài liệu đầu vào nên báo cáo trình bày kết MIKE11 Khái quát vùng nghiên cứu Lưu vực sông Thu Bồn thuộc địa phận tỉnh Quảng Nam Thành phố Đà Nẵng, giới hạn từ 14054 đến 16013 vĩ Bắc 107013 đến 108044 kinh Đông Phía Bắc giáp Thừa Thiên - Huế, phía Nam giáp Quảng Ngãi, phía Tây giáp Kon Tum nước cộng hoà dân chủ nhân dân Lào, phía Đông biển Diện tích tự nhiên lưu vực 10.350 km2 Địa hình lưu vực phức tạp, phần lớn đồi núi Đường chia nước qua đỉnh núi cao đỉnh A Tuất 500m, đỉnh Mang 1708m, đỉnh Ba Na 1483m, đỉnh Ngọc Lĩnh 2598m, đỉnh Hòn Ba 1358m, đỉnh Mũi Chúa 1362m Lưu vực bao bọc núi cao phía với độ dốc lớn tạo nên khả tập trung nước nhanh có lũ Dải đồng hẹp xen kẽ với đồi bát úp nơi nhận nước trực tiếp từ phần thượng lưu lưu vực Hình Bản đồ khu vực nghiên cứu Hệ thống gồm sông Vũ Gia phía bắc Thu Bồn phía nam Hai sông có liên hệ thuỷ lực qua khu vực sông Quảng Huế, nơi dòng chảy tập trung lòng dẫn mùa kiệt chảy tràn qua bãi chủ yếu từ sông Vũ Gia sang sông Thu Bồn vào mùa lũ Khu vực hạ lưu vùng đồng thành phố Đà Nẵng, Hội An, nơi có mật độ dân số cao trung tâm kinh tế xã hội vùng Khu vực có chế độ thuỷ văn không đồng với mùa lũ ngắn tháng từ tháng IX đến tháng XII, tổng lượng dòng chảy chiếm tới 80% năm làm khu vực ngã ba Quảng Huế thường xuyên nước tràn bờ, xói bãi tạo lòng mới, sông Lũ năm 1998, 1999, 2000 phá hủy đoạn cong cuối tạo dòng nối sang sông Thu bồn với chiều dài 1.1 km chiều rộng 80 - 100m Về mặt thủy lực việc cắt dòng hoàn toàn logic đường chảy dòng nước ngắn Nhưng sang mùa cạn phần lớn dòng chảy từ thượng nguồn Vũ Gia chuyển sang Thu Bồn gây tình trạng thiếu nước cho phần lưu vực phía bắc, nơi có 14,000 đất nông nghiệp nhu cầu nuớc cho sinh hoạt thành phố Đà Nẵng Đồng thời với trình tồn nhánh sông sông Quảng Huế cũ với chiều dài km bị bồi lấp dần bị chết vào mùa cạn Để khắc phục tình trạng đó, tháng III/2001 nhân dân địa phương xây dựng đập tạm chặn dòng mới, lưu lượng chảy sông Ái nghĩa không mong muốn Tháng VI/2001 Viện Khoa học Thủy lợi nghiên cứu, đề xuất giải pháp xây dựng đập chìm chặn cửa sông mở lũ tới cao trình 4.5m củng cố đập tạm dân đắp; củng cố bờ kênh phải trái tới cao trình tự nhiên 6.5m Bờ bảo vệ kéo dài 50m thượng lưu 100m hạ lưu phía bờ trái; 150m thượng lưu 50 m hạ lưu bờ phải Tuy nhiên sau trận lũ tháng X/2001, công trình xử lý khẩn cấp lại bị phá hướng dòng sông theo xu tạo chế độ thủy lực thích hợp Nếu để dòng sông phát triển theo qui luật tự nhiên làm đảo lộn hoạt động kinh tế xã hội theo hướng bất lợi Đó tình hình thiếu nước nghiêm trọng cho phần phía bắc dẫn tới tình hình nhiễm mặn nghiêm trọng hơn; nguồn nước cấp cho thành phố Đà Nẵng không đủ mùa cạn, mùa lũ, dòng chảy tập trung gần hết sang nhánh Thu Bồn làm cho tình hình lũ lụt phần phía Nam tỉnh nghiêm trọng Trước tình hình đó, cần thiết phải có nghiên cứu tổng thể nhằm tìm kiếm giải pháp hợp lý đáp ứng yêu cầu sau: • Tạo đường dẫn ổn định lâu dài năm thuộc nhánh sông hệ thống, phù hợp với chế độ thuỷ lực • Tạo phân phối dòng chảy hợp lý nhánh sông, đặc biệt mùa kiệt, tránh gây tình trạng thiếu nước cho phần phía bắc Việc ứng dụng phần mềm MIKE với mục tiêu tìm thông số thuỷ văn, thuỷ lực, bùn cát thực chất khu vực để có sở đề xuất giải pháp chỉnh trị cho đoạn sông nghiên cứu 3 Ứng dụng mô hình MIKE11 tính toán lũ hệ thống sông Thu Bồn – Vũ Gia 3.1.Giới thiệu mô hình MIKE11 MIKE 11, DHI Water & Environment phát triển, gói phần mềm dùng để mô dòng chảy, lưu lượng, chất lượng nước vận chuyển bùn cát cửa sông, sông, kênh tưới vật thể nước khác Mô hình xây dựng sở hệ phương trình Saint Venant Hệ phương trình gồm hai phương trình: Phương trình liên tục: ∂Q ∂A + = q ∂x ∂t Phương trinh động lượng: α QQ ∂Q ∂ Q2 ∂h + (β ) + gA + g = ∂t ∂x A ∂x C RA Giải hệ phương trình theo phương pháp sai phân hữu hạn điểm (AbbottIonescu 6-point) xác định giá trị lưu lượng, mực nước đoạn sông mặt cắt ngang mạng sông thời điềm khoảng thời gian tính toán 3.2.Ứng dụng mô hình MIKE11 tính toán lũ hệ thống sông Thu Bồn – Vũ Gia 3.2.1.Sơ đồ mạng lưới biên tính toán 3.2.1.1.Sơ đồ mạng lưới Mạng lưới sông mô trạm Thành Mỹ dòng Vũ Gia trạm Nông Sơn sông Thu Bồn Hệ thống sông kể từ điềm đầu vào mô gồm nhánh sông Vũ Gia, sông Thu Bồn, sông Bung, sông Kôn, sông Cẩm Lệ, sông Vĩnh Điện, sông Trường Giang Sông Quảng Huế đoạn sông nhỏ nối hai dòng sông Vũ Gia Thu Bồn Các sông mạng lưới có độ uốn khúc cao, hạ lưu sông Vũ Gia có chia dòng rõ ràng [meter] Standard - Run11.res11 Trạm Đà Nẵng 1784000.0 1782000.0 1780000.0 Trạm Cẩm Lệ 1778000.0 1776000.0 1774000.0 1772000.0 1770000.0 1768000.0 1766000.0 1764000.0 Trạm Ái Nghĩa 1762000.0 Trạm Hội An Trạm Câu Lâu 1760000.0 1758000.0 1756000.0 Trạm Hội Khách 1754000.0 1752000.0 1750000.0 1748000.0 Trạm Giao Thủy 1746000.0 1744000.0 1742000.0 1740000.0 Trạm Thành Mỹ 800000.0 Trạm Nông Sơn 810000.0 820000.0 830000.0 840000.0 850000.0 Trạm Tam Kỳ 860000.0 870000.0 [meter] Hình Sơ đồ thủy lực hệ thống sông Thu Bồn – Vũ Gia 3.2.1.2.Biên tính toán Các biên đầu cuối mô hình gồm có: - biên lưu lượng thực đo (Q ~ t) trạm đo thượng nguồn trạm Thành Mỹ sông Vũ Gia trạm Nông Sơn sông Thu Bồn Mặt cắt lòng sông trạm có đặc điểm hẹp, bờ sông dốc khống chế lưu lượng lũ lớn chảy qua Hai trạm có số liệu đo lưu lượng mực nước từ năm 1995 – 2006 - biên mực nước (H ~ t) vị trí Hội An, Đà Nẵng sông Trường Giang Tài liệu mực nước thực đo thời gian với biên lưu lượng có Đà Nẵng Hội An Mực nước biên sông Trường Giang, nội suy từ mực nước Hội An - biên nhập lưu (Q ~ t) dạng tập trung (Point source) gồm có: Sông Bung sông Kôn chảy vào sông Vũ Gia phía bờ trái Quá trình dòng chảy biên nhập lưu mô mô hình MIKE NAM, modul diễn toán thủy văn theo phương pháp mưa dòng chảy Diện tích lưu vực cho bảng sau: Lưu vực Diện tích (km2) Thu Bồn 3150 Vũ Gia 1850 Bung 2530 Kôn 627 - Các biên nhập lưu (Q ~ t) dạng phân tán (Distributed sources) trải dọc theo dòng hai sông Vũ Gia Thu Bồn Các biên nhập lưu phân tán mô mô hình MIKE NAM với nguyên tắc cho tổng lượng dòng chảy phù hợp theo quy luật diện tích so với sông Các biên nhập lưu dọc sông bao gồm: • biên sông Vũ Gia có diện tích 83 km2, 86 km2, 77 km2 • biên sông Thu Bồn có diện tích 168 km2 34 km2 Tổng cộng có biên nhập lưu toàn hệ thống 3.2.1.3.Trạm kiểm tra Tài liệu thủy văn trạm kiểm tra thu thập gồm tài liệu thực đo trạm sau: • Quá trình H ~ t trạm Hội Khách sông Vũ Gia (t = 1h) • Quá trình H ~ t trạm Ái Nghĩa sông Vũ Gia (t = 1h) • Quá trình H ~ t trạm Giao Thủy sông Thu Bồn (t = 1h) • Quá trình H ~ t trạm Câu Lâu sông Thu Bồn (t = 1h) • Quá trình H ~ t trạm Cẩm Lệ sông Cẩm Lệ (t = 1h) 3.2.2.Tài liệu phục vụ tính toán 3.2.2.1.Tài liệu địa hình Hệ thống mô gồm 64 mặt cắt Tài liệu mặt cắt đo đạc vào tháng 12/2007 có hệ cao độ thống đủ độ tin cậy cho tính toán Tuy nhiên cần lưu ý phần hạ lưu kể từ ngã ba Quảng Huế xuống biển phần đồng hẹp với địa hình phẳng nên lũ mức trung bình nước ngập toàn đồng bằng, khu vực trở thành khu chứa lũ Số lượng mặt cắt lại không đủ dày, vị trí mặt cắt không xác, cao độ điểm đầu cao mặt cắt không xác định khó khăn việc thành lập mạng lưới sông tính toán mô hình Tất yếu tố ảnh hưởng không nhỏ việc tính toán hiệu chỉnh mô hình Trên thực tế, ta vào số dấu hiệu file đo đạc kết hợp với GoogleMap để thiết lập mạng lưới sông 3.2.2.2.Tài liệu thủy văn Như trình bày trên, tài liệu thủy văn gồm trạm biên trạm kiểm tra 3.2.3.Hiệu chỉnh mô hình 3.2.3.1.Chọn thời đoạn tính toán Căn vào nhiệm vụ toán, số liệu thực đo thu thập được, để hiệu chỉnh thông số mô hình, sử dụng tài liệu thực đo trạm thủy văn nói từ 1h ngày 04/12/2006 đến 22h ngày 08/12/2006 Và để kiểm định thông số mô hình sử dụng tài liệu thực đo trạm cho thời đoạn khác năm 1995, 2002 2003 3.2.3.2.Hiệu chỉnh thông số mô hình Việc hiệu chỉnh thông số mô hình chủ yếu tiến hành cách thay đổi độ nhám lòng dẫn Phương pháp hiệu chỉnh thông số phương pháp thử dần Quá trình hiệu chỉnh tóm tắt thành bước sau đây: Bước 1: Giả thiết thông số (chủ yếu độ nhám), điều kiện ban đầu Bước 2: Sau có thông số giả thiết, tiến hành chạy mô hình Bước 3: So sánh kết tính toán với số liệu thực đo trạm có số liệu đo đạc lưu lượng mực nước Việc so sánh tiến hành trực quan (so sánh hai đường trình tính toán thực đo biểu đồ), đồng thời kết hợp tiêu Nash để kiểm tra Nash = - ∑ ( Xo, i − Xs, i ) ∑ ( Xo, i − Xo ) 2 Xo,i: Giá trị thực đo Xs,i: Giá trị tính toán mô Xo : Giá trị thực đo trung bình Bước 4: Nếu kết so sánh tốt dừng hiệu chỉnh lưu thông số Nếu kết không đạt, tiến hành phân tích đánh giá sai lệch, sau tiếp tục hiệu chỉnh lại thông số Chạy mô hình So sánh thực đo tính toán Đạt Dừng Không đạt Giả thiết thông số Thay đổi thông số Hình Sơ đồ trình hiệu chỉnh thông số mô hình 3.2.3.3.Kết hiệu chỉnh mô hình năm 2006 Ta chọn thời đoạn lũ điển hình với mực nước khống chế cho mặt cắt tượng chảy tràn, cụ thể sau: Mực nước Ái Nghĩa không vượt 8,5 m Mực nước Giao Thủy không vượt 8,0 m Mực nước Cẩm Lệ không vượt 2,1 m Mực nước Câu Lâu không vượt 4,6 m Căn vào tài liệu mực nước lũ thực đo ta định chọn thời đoạn sau phục Vũ chạy kiểm định mô hình: Từ 00:00 ngày 5/XII/2006 đến 00:00 ngày 9/XII/2006 Theo sơ đồ thủy lực hệ thống sông Thu Bồn – Vũ Gia hình vẽ trên, ta thấy sông Quảng Huế nối tiếp sông Vũ Gia sông Thu Bồn Có thể nói khu vực “nhạy cảm”, chế độ phân lưu thủy lực phức tạp Sự sơ sài tài liệu mặt cắt dẫn đến việc kiểm định trạm Ái Nghĩa Giao Thủy tương đối phức tạp Sau kết kiểm định mô hình cho thời đoạn chọn • Trạm Hội Khách (sông Vũ Gia) [meter] Time Series Water Level 15.0 14.5 14.0 13.5 13.0 12.5 12.0 11.5 11.0 10.5 10.0 9.5 9.0 8.5 12:00:00 4-12-2006 00:00:00 5-12-2006 12:00:00 00:00:00 6-12-2006 12:00:00 00:00:00 7-12-2006 12:00:00 00:00:00 8-12-2006 12:00:00 00:00:00 9-12-2006 Hình Mực nước trạm Hội Khách (đường tính toán màu đen, thực đo màu xanh) Quá trình đường mực nước tương đối giống cho đường tính toán thực đo Đây xu phù hợp với thực tế trạm Hội Khách gần với biên Thành Mỹ nên pha lũ truyền từ Thành Mỹ xuất tương đối giống Đỉnh lũ xuất vào khoảng 6h ngày 06/XII/2006 Tuy nhiên, để có phù hợp đỉnh lũ nói chung đường thực đo – tính toán nói riêng, ta có hiệu chỉnh độ nhám lòng sông Cụ thể khu vực lòng sông, hệ số nhám Manning n = 0.11 (lòng sông), n = 0.085 (giữa) n = 0.03 (bãi sông) Điều trái với quy luật thực tế thông thường phần bãi sông thường nhám so với phần lòng sông, nghĩa hệ số nhám lòng sông phải nhỏ Thêm vào đó, việc hiệu chỉnh hệ số nhám theo vùng (3 zone) có ảnh hưởng rõ với đường mực nước tính toán Sự thay đổi độ nhám đột ngột từ dẫn đến tượng đường mực nước tính toán không ổn định cao trình 12m Hệ số Nash = 87.79% • Trạm Ái Nghĩa (cuối sông Vũ Gia) [meter] Time Series Water Level 8.5 8.0 7.5 7.0 6.5 6.0 5.5 5.0 4.5 4.0 3.5 3.0 12:00:00 4-12-2006 00:00:00 5-12-2006 12:00:00 00:00:00 6-12-2006 12:00:00 00:00:00 7-12-2006 12:00:00 00:00:00 8-12-2006 12:00:00 00:00:00 9-12-2006 Hình Mực nước trạm Ái Nghĩa (đường tính toán màu đen, thực đo màu xanh) Từ chân lũ đến đỉnh lũ, đường tính toán thực đo tương đối giống Tuy nhiên, đỉnh lũ tính toán nhanh so với thực đo khoảng 1h Tuy nhiên, nửa sau lũ, xu hướng đường nước tính toán nhỏ thực đo khoảng 20cm Điều khắc phục cách nâng lưu lượng biên nhập lưu cục sông Bung sông Kon thời đoạn nửa cuối lũ Tuy nhiên, nói trên, chế độ thủy lực khu vực phức tạp nên việc thay đổi biên cục dẫn đến diễn biến thủy lực khác không theo mong muốn xảy Theo nhiều phương án thay đổi biên lưu lượng thử nghiệm phương án đem lại kết tốt Hệ số Nash = 92.59% • Trạm Giao Thủy (sông Thu Bồn) 10 [meter] Time Series Water Level 8.0 7.5 7.0 6.5 6.0 5.5 5.0 4.5 4.0 3.5 3.0 2.5 2.0 12:00:00 4-12-2006 00:00:00 5-12-2006 12:00:00 00:00:00 6-12-2006 12:00:00 00:00:00 7-12-2006 12:00:00 00:00:00 8-12-2006 12:00:00 00:00:00 9-12-2006 Hình Mực nước trạm Giao Thủy (đường tính toán màu đen, thực đo màu xanh) Như hình vẽ trên, ta thấy việc kiểm định trạm Giao Thủy tốt Tuy nhiên, đoạn cuối lũ, đường tính toán nhỏ so với thực đo lượng lớn Ta cần có phương án thích hợp để chỉnh trạm Ái Nghĩa Giao Thủy để có kết tốt Hệ số Nash = 82.71% • Trạm Cẩm Lệ (sông Cẩm Lệ) [meter] Time Series Water Level 1.3 1.2 1.1 1.0 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0.0 -0.1 -0.2 -0.3 12:00:00 4-12-2006 00:00:00 5-12-2006 12:00:00 00:00:00 6-12-2006 12:00:00 00:00:00 7-12-2006 11 12:00:00 00:00:00 8-12-2006 12:00:00 00:00:00 9-12-2006 12:00:00 00:00:00 10-12-2006 12:00:00 Hình Mực nước trạm Cẩm Lệ (đường tính toán màu đen, thực đo màu xanh) Hệ số Nash = 87.23% • Trạm Câu Lâu (cuối sông Thu Bồn) [meter] Time Series Water Level 4.0 3.8 3.6 3.4 3.2 3.0 2.8 2.6 2.4 2.2 2.0 1.8 1.6 1.4 1.2 1.0 0.8 0.6 0.4 0.2 0.0 12:00:00 4-12-2006 00:00:00 5-12-2006 12:00:00 00:00:00 6-12-2006 12:00:00 00:00:00 7-12-2006 12:00:00 00:00:00 8-12-2006 12:00:00 00:00:00 9-12-2006 12:00:00 00:00:00 10-12-2006 12:00:00 Hình Mực nước trạm Câu Lâu (đường tính toán màu đen, thực đo màu xanh) Hệ số Nash = 86.76% Trạm Câu Lâu Cẩm Lệ trạm thuộc cuối sông Do vậy, việc hiệu chỉnh kiểm định mô hình không quan trọng Chính lý nên khuôn khổ nghiên cứu này, ta không thực nhiều việc hiệu chỉnh trạm Cẩm Lệ Câu Lâu 3.2.4.Kiểm định mô hình Từ kết hiệu chỉnh mô hình ta có thông số mô hình, dùng thông số tiến hành chạy kiểm tra năm Công việc kiểm định mô hình nhằm mục đích xem xét tính ổn định thông số mô hình năm 2006 Kết kiểm định mô hình năm trình bày dạng tổng hợp sau 3.2.4.1.Kiểm định cho năm 1995 Thời đoạn chọn kiểm định năm 1995 từ 01:00 ngày 8/X/1995 đến 23:00 ngày 10/X/1995 • Đường trình H ~ t 12 [meter] Time Series Water Level 23.0 22.0 21.0 20.0 19.0 18.0 17.0 16.0 15.0 14.0 13.0 12.0 11.0 10.0 9.0 8.0 7.0 6.0 5.0 4.0 3.0 2.0 1.0 0.0 06:00:00 8-10-1995 12:00:00 18:00:00 00:00:00 9-10-1995 06:00:00 12:00:00 18:00:00 00:00:00 10-10-1995 06:00:00 12:00:00 18:00:00 Hình Đường trình H ~ t cho trạm (từ xuống dưới: Thành Mỹ, Hội Khách, Nông Sơn, Ái Nghĩa, Giao Thủy, Câu Lâu, Cẩm Lệ) • Đường trình Q ~ t [m^3/s] Time Series Discharge 4200.0 4000.0 3800.0 3600.0 3400.0 3200.0 3000.0 2800.0 2600.0 2400.0 2200.0 2000.0 1800.0 1600.0 1400.0 1200.0 1000.0 800.0 600.0 400.0 200.0 0.0 06:00:00 8-10-1995 12:00:00 18:00:00 00:00:00 9-10-1995 06:00:00 12:00:00 18:00:00 00:00:00 10-10-1995 06:00:00 12:00:00 18:00:00 Hình Đường trình Q ~ t cho số vị trí Trên hình vẽ thể đường trình lưu lượng số mặt cắt Đường màu đen lấy mặt cắt trước cửa vào sông Quảng Huế, đường màu xanh nước biển lấy gần trạm Ái Nghĩa, đường màu xanh mạ lấy sông Quảng Huế mới, 13 đường màu đỏ lấy sông Thu Bồn đoạn trước gặp sông Quảng Huế đường màu hồng gần trạm Giao Thủy Biểu đồ thể tương đối rõ tính liên tục lưu lượng hệ thống sông Đường màu đen tổng đường màu xanh nước biển xanh mạ, đường màu hồng tổng đường màu xanh mạ màu đỏ Tuy nhiên, tượng tràn bãi ảnh hưởng đến xác biểu đồ • Đường trình V ~ t [m/s] Time Series Flow velocity 2.0 1.9 1.8 1.7 1.6 1.5 1.4 1.3 1.2 1.1 1.0 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0.0 06:00:00 8-10-1995 12:00:00 18:00:00 00:00:00 9-10-1995 06:00:00 12:00:00 18:00:00 00:00:00 10-10-1995 06:00:00 12:00:00 18:00:00 Hình 10 Đường trình V ~ t cho số mặt cắt sông Quảng Huế Biểu đồ thể lưu tốc số mặt cắt sông Quảng Huế Ở đầu sông, lưu tốc nói chung tương đối lớn (đường màu đen) có xu hướng giảm dần phía cuối 3.2.4.2.Kiểm định cho năm 2002 Thời đoạn chọn kiểm định năm 2002 từ 01:00 ngày 24/IX/2002 đến 22:00 ngày 26/IX/2002 • Đường trình H ~ t 14 [meter] Time Series Water Level 16.0 15.0 14.0 13.0 12.0 11.0 10.0 9.0 8.0 7.0 6.0 5.0 4.0 3.0 2.0 1.0 0.0 06:00:00 24-9-2002 12:00:00 18:00:00 00:00:00 25-9-2002 06:00:00 12:00:00 18:00:00 00:00:00 26-9-2002 06:00:00 12:00:00 18:00:00 Hình 11 Đường trình H ~ t cho trạm (từ xuống dưới: Thành Mỹ, Hội Khách, Nông Sơn, Ái Nghĩa, Giao Thủy, Câu Lâu, Cẩm Lệ) • Đường trình Q ~ t [m^3/s] Time Series Discharge 4400.0 4200.0 4000.0 3800.0 3600.0 3400.0 3200.0 3000.0 2800.0 2600.0 2400.0 2200.0 2000.0 1800.0 1600.0 1400.0 1200.0 1000.0 800.0 600.0 400.0 200.0 0.0 06:00:00 24-9-2002 12:00:00 18:00:00 00:00:00 25-9-2002 06:00:00 12:00:00 18:00:00 00:00:00 26-9-2002 06:00:00 12:00:00 18:00:00 Hình 12 Đường trình Q ~ t cho số vị trí Trên hình vẽ thể đường trình lưu lượng số mặt cắt Đường màu đen lấy mặt cắt trước cửa vào sông Quảng Huế, đường màu xanh nước biển lấy gần trạm Ái Nghĩa, đường màu xanh mạ lấy sông Quảng Huế mới, 15 đường màu đỏ lấy sông Thu Bồn đoạn trước gặp sông Quảng Huế đường màu hồng gần trạm Giao Thủy Biểu đồ thể tương đối rõ tính liên tục lưu lượng hệ thống sông Đường màu đen tổng đường màu xanh nước biển xanh mạ, đường màu hồng tổng đường màu xanh mạ màu đỏ • Đường trình V ~ t [m/s] Time Series Flow velocity 0.85 0.80 0.75 0.70 0.65 0.60 0.55 0.50 0.45 0.40 0.35 0.30 0.25 0.20 0.15 06:00:00 24-9-2002 12:00:00 18:00:00 00:00:00 25-9-2002 06:00:00 12:00:00 18:00:00 00:00:00 26-9-2002 06:00:00 12:00:00 18:00:00 Hình 13 Đường trình V ~ t cho số mặt cắt sông Quảng Huế Biểu đồ thể lưu tốc số mặt cắt sông Quảng Huế Ở đầu sông, lưu tốc nói chung nhỏ (đường màu đen), sau cao phần bãi (đường màu xanh nước biển), nhỏ dần cuối đoạn sông Quảng Huế 3.2.4.3.Kiểm định cho năm 2003 Thời đoạn chọn kiểm định năm 2003 từ 01:00 ngày 16/X/2003 đến 23:00 ngày 19/X/2003 • Đường trình H ~ t 16 [meter] Time Series Water Level 19.0 18.0 17.0 16.0 15.0 14.0 13.0 12.0 11.0 10.0 9.0 8.0 7.0 6.0 5.0 4.0 3.0 2.0 1.0 0.0 06:00:00 12:00:00 16-10-2003 18:00:00 00:00:00 06:00:00 17-10-2003 12:00:00 18:00:00 00:00:00 06:00:00 18-10-2003 12:00:00 18:00:00 00:00:00 06:00:00 19-10-2003 12:00:00 18:00:00 Hình 14 Đường trình H ~ t cho trạm (từ xuống dưới: Thành Mỹ, Hội Khách, Nông Sơn, Ái Nghĩa, Giao Thủy, Câu Lâu, Cẩm Lệ) • Đường trình Q ~ t [m^3/s] Time Series Discharge 7000.0 6500.0 6000.0 5500.0 5000.0 4500.0 4000.0 3500.0 3000.0 2500.0 2000.0 1500.0 1000.0 500.0 0.0 -500.0 06:00:00 12:00:00 16-10-2003 18:00:00 00:00:00 06:00:00 17-10-2003 12:00:00 18:00:00 00:00:00 06:00:00 18-10-2003 12:00:00 18:00:00 00:00:00 06:00:00 19-10-2003 12:00:00 18:00:00 Hình 15 Đường trình Q ~ t cho số vị trí Trên hình vẽ thể đường trình lưu lượng số mặt cắt Đường màu đen lấy mặt cắt trước cửa vào sông Quảng Huế, đường màu xanh nước biển lấy gần trạm Ái Nghĩa, đường màu xanh mạ lấy sông Quảng Huế mới, 17 đường màu đỏ lấy sông Thu Bồn đoạn trước gặp sông Quảng Huế đường màu hồng gần trạm Giao Thủy Biểu đồ thể tương đối rõ tính liên tục lưu lượng hệ thống sông Đường màu đen tổng đường màu xanh nước biển xanh mạ, đường màu hồng tổng đường màu xanh mạ màu đỏ Từ biều vẽ, ta thấy xuất thời đoạn lưu lượng âm sông Quảng Huế Về thực chất, lưu lượng âm mà thời điểm đó, nước không chảy sông Quảng Huế, mô hình tính toán để đảm bảo tính liên tục tính toán • Đường trình V ~ t [m/s] Time Series Flow velocity 0.26 0.24 0.22 0.20 0.18 0.16 0.14 0.12 0.10 0.08 0.06 0.04 0.02 0.00 -0.02 -0.04 -0.06 -0.08 -0.10 -0.12 06:00:00 12:00:00 16-10-2003 18:00:00 00:00:00 06:00:00 17-10-2003 12:00:00 18:00:00 00:00:00 06:00:00 18-10-2003 12:00:00 18:00:00 00:00:00 06:00:00 19-10-2003 12:00:00 18:00:00 Hình 16 Đường trình V ~ t cho số mặt cắt sông Quảng Huế Biểu đồ thể lưu tốc số mặt cắt sông Quảng Huế Ở đầu sông, lưu tốc nói chung nhỏ (đường màu đen), sau cao phần bãi (đường màu xanh nước biển), nhỏ dần cuối đoạn sông Quảng Huế Trên biểu đồ, có thời đoạn đầu cuối xuất không ổn định Có thể nguyên nhân chưa hợp lý mặt cắt thông số ban đầu Kết luận • Bài toán thực với việc áp dụng mô hình MIKE11, mô hình thương mại sử dụng tính toán nhiều Việt Nam cho độ xác cao • Về số liệu mặt cắt địa hình: Khu vực có số liệu mặt cắt tương đối ít, không rõ ràng vị trí số mặt cắt gây khó khăn tính toán hiệu chỉnh mô hình 18 • Về số liệu khí tượng thủy văn: Chuỗi số liệu mưa năm 2006 số liệu triều Đà Nẵng Hội An sử dụng cho tính toán • Từ kết tính toán mô hình chiều, ta có số liệu đầu vào cho mô hình chiều tính toán cụ thể cho khu vực sông Quảng Huế Để tiếp tục phát triển mô hình đòi hỏi thời gian tương đối dài, phần tính toán chiều không đề cập đến báo cáo 19 ... dòng Vũ Gia trạm Nông Sơn sông Thu Bồn Hệ thống sông kể từ điềm đầu vào mô gồm nhánh sông Vũ Gia, sông Thu Bồn, sông Bung, sông Kôn, sông Cẩm Lệ, sông Vĩnh Điện, sông Trường Giang Sông Quảng Huế. .. trực tiếp từ phần thượng lưu lưu vực Hình Bản đồ khu vực nghiên cứu Hệ thống gồm sông Vũ Gia phía bắc Thu Bồn phía nam Hai sông có liên hệ thu lực qua khu vực sông Quảng Huế, nơi dòng chảy tập... thông số thu văn, thu lực, bùn cát thực chất khu vực để có sở đề xuất giải pháp chỉnh trị cho đoạn sông nghiên cứu 3 Ứng dụng mô hình MIKE11 tính toán lũ hệ thống sông Thu Bồn – Vũ Gia 3.1.Giới