BỘ TÀI NGUN VÀ MƠI TRƯỜNG TRUNG TÂM KHÍ TƯỢNG THUỶ VĂN QUỐC GIA Số – Đặng Thái Thân – Hoàn Kiếm – Hà Nội *** - BÁO CÁO TỔNG KẾT CÁC CHUYÊN ĐỀ NGHIÊN CỨU THUỘC ĐỀ TÀI NGHIÊN CỨU KHOA HỌC VÀ CÔNG NGHỆ CẤP BỘ NGHIÊN CỨU ỨNG DỤNG MÔ HÌNH WETSPA VÀ HECRAS MƠ PHỎNG, DỰ BÁOQ TRÌNH LŨ TRÊN HỆ THỐNG SÔNG THU BỒN – VU GIA CHỦ NHIỆM ĐỀ TÀI: Ths Đặng Thanh Mai Trung tâm Dự báo Khí tượng Thủy văn Trung ương 7517-1 15/10/2009 Hà Nội, – 2009 BỘ TÀI NGUYÊN VÀ MÔI TRƯỜNG TRUNG TÂM KHÍ TƯỢNG THỦY VĂN QUỐC GIA Số - Đặng Thái Thân – Hoàn Kiếm – Hà Nội BÁO CÁO TỔNG KẾT ĐỀ TÀI NGHIÊN CỨU KHOA HỌC Tên đề tài: NGHIÊN CỨU ỨNG DỤNG MƠ HÌNH WETSPA VÀ HECRAS MƠ PHỎNG, DỰ BÁO Q TRÌNH LŨ TRÊN HỆ THỐNG SÔNG THU BỒN – VU GIA Chỉ số phân loại: Chỉ số đăng ký: Chỉ số lưu trữ: Các quan cộng tác viên tham gia thực đề tài: Ks Bùi Đức Long Ths Vũ Đức Long Ths Nguyễn Thị Thu Trang Ths Võ Văn Hòa Ts Nguyễn Thanh Long Ts Lê Quốc Hùng Ths Phạm Văn Chiến Ks Nguyễn Hữu Thiêm Ks Phùng Hồng Long Ngày tháng năm 2009 CHỦ NHIỆM ĐỀ TÀI Ngày tháng năm 2009 CƠ QUAN THỰC HIỆN ĐỀ TÀI Bùi Minh Tăng Ths Đặng Thanh Mai Hà Nội, ngày tháng Trung tâm Dự báo KTTV Trung ương Trung tâm Dự báo KTTV Trung ương Trung tâm Dự báo KTTV Trung ương Trung tâm Dự báo KTTV Trung ương Viện Nghiên cứu Địa chất khoáng sản Viện Nghiên cứu Địa chất khoáng sản Đài KTTV khu vực Trung Trung Đài KTTV khu vực Trung Trung Đài KTTV khu vực Trung Trung năm 2009 HỘI ĐỒNG ĐÁNH GIÁ CHÍNH THỨC CHỦ TỊCH HỘI ĐỒNG Hà Nội, ngày Ngày tháng năm 2009 CƠ QUAN CHỦ TRÌ ĐỀ TÀI Trần Văn Sáp tháng năm 2009 CƠ QUAN QUẢN LÝ ĐỀ TÀI TL BỘ TRƯỞNG KT.VỤ TRƯỞNG VỤ KHOA HỌC CƠNG NGHỆ PHĨ VỤ TRƯỞNG Đề tài: Nghiên cứu ứng dụng mơ hình WETSPA HECRAS mơ phỏng, dự báo q trình lũ hệ thống sơng Thu Bồn – Vu Gia Chun đề 1: Nghiên cứu, tính tốn xử lý sản phẩm mưa 2A12 từ ảnh vệ tinh TRMM kết hợp với số liệu mưa mặt đất trận lũ lớn từ 1998-2005 lưu vực hệ thống sông để phù hợp với yêu cầu đầu vào mơ hình Trong nghiệp vụ dự báo lũ, vấn đề xử lý lượng mưa xảy lưu vực hình thời tiết gây phức tạp, mưa phân bố không theo không gian thời gian Đối với chế độ thủy văn hệ thống sông Thu Bồn – Vu Gia, tổng lượng mưa, cường độ mưa, thời gian mưa phân bố mưa định q trình hình thành dịng chảy lũ lưu vực Mưa lớn tập trung vùng núi sinh lũ lớn vùng đồng mưa lớn hay mưa nhỏ Nhưng mưa lớn tập trung vùng đồng bằng, vùng núi mưa nhỏ khơng gây lũ lớn sơng mà gây úng ngập cục số khu vực Trong đợt mưa lớn mà lượng mưa tập trung khắp lưu vực vùng núi vùng đồng nguy hiểm - vùng núi tốc độ dòng chảy lũ lớn vùng đồng xảy ngập lụt nghiêm trọng Việc xác định cách định lượng xác phân bố mưa hình thời tiết lưu vực sơng khó điều kiện mưa yếu tố ngẫu nhiên chịu tác động tổ hợp nhiều trình phức tạp mang tính ngẫu nhiên Tuy nhiên, việc xác định cách định tính phân bố mưa theo khơng gian thời gian hình thời tiết gây mưa thực phần dựa vào việc nghiên cứu, phân tích đợt mưa sinh lũ xảy khứ để đưa đặc điểm chung Đây phần quan trọng trình cảnh báo dự báo lũ lưu vực Mặt khác, yêu cầu đầu vào mơ u cầu cải tiến chất lượng mơ hình dự báo nghiệp vụ, cần thiết phải xét đến khác biệt phân bố mưa trận lũ hình thời tiết gây Để giải vấn đề Đề tài tập trung Nghiên cứu, phân tách, chuyển đổi sản phẩm mưa 2A12 ảnh vệ tinh TRMM thành dạng đồ RASTER Đặc điểm mưa hình thời tiết gây lũ lớn 1.1 Các hình thời tiết gây mưa lớn Việc phân tích ngun nhân hình thành trận lũ lớn xảy hệ thống sông Vu Gia – Thu Bồn từ 1976 đến 2007 (bảng 5.2) cho thấy, hình thời tiết gây mưa lũ lưu vực là: Bão, ATNĐ, KKL, dải HTNĐ, đới gió đơng tổ hợp chúng a/ Bão, áp thấp nhiệt đới (ATNĐ) Theo số liệu thống kê nhiều năm, khu vực Quảng Nam – Đà Nẵng vùng chịu ảnh hưởng nhiều bão ATNĐ, chiếm tới 20% số bão, ATNĐ ảnh hưởng đến Việt Nam Trung bình hàng năm có bão ATNĐ ảnh hưởng đến khu vực Thông thường, bão đổ vào khu vực Quảng Nam – Đà Nẵng nhiều tháng IX - XI (chiếm 70%) Đáng lưu ý, bão không đổ trực tiếp vào lưu vực, song chịu ảnh hưởng gián tiếp gây thiệt hại nặng nề mưa, lũ gây Có trường hợp bão mạnh bão đổ liên tiếp thời gian ngắn Như tháng IX/1964, IX/1978, khoảng 15 ngày, bão mạnh đổ liên tiếp vào miền Trung gây lũ lụt lịch sử lũ lụt đặc biệt lớn tồn miền, có khu vực Quảng Nam - Đà Nẵng, gây thiệt hại lớn Trường hợp bão bão ATNĐ đổ liên tiếp vào khu vực gây thiên tai nghiêm trọng Các loại hình thời tiết ảnh hưởng trực tiếp gián tiếp gây mưa lớn diện mưa rộng, thời gian mưa kéo dài, gây lũ lớn khắp sơng hệ thống sơng b/ Khơng khí lạnh (KKL) Áp cao lạnh ảnh hưởng gây mưa lớn tỉnh miền Trung bao gồm hình áp cao lạnh tác động đơn áp cao lạnh kết hợp với hình khác rãnh áp thấp bị nén, tín phong nhiễu động nhiệt đới khác dải hội tụ nhiệt đới, xoáy thuận nhiệt đới (XTNĐ) hay nhiễu động sóng đới gió đơng cao Trong loại hình synop đặc trưng gây mưa lớn tỉnh miền Trung loại hình áp cao lạnh kết hợp với dải hội tụ nhiệt đới XTNĐ gây mưa mãnh liệt Thực tế cho thấy, áp cao lạnh ảnh hưởng thường mưa miền Trung khơng phải đợt gây mưa lớn mà hồn tồn phụ thuộc hoàn lưu cao lớp tầng đối lưu tiêu biểu mực 1500m Thời gian bắt đầu kết thúc mưa, cường độ mưa, thời gian trì mưa lớn, phạm vi mưa khơng đồng đợt xâm nhập KKL Nó khơng phụ thuộc vào tác động rìa phía nam hay tây nam áp cao lạnh khu vực đơng nam Trung Quốc mà cịn phụ thuộc vào điều kiện mặt đệm, địa hình Thời gian hoạt động KKL gây mưa cho lưu vực thường vào tháng X, XI, XII c/ Ảnh hưởng kết hợp loại hình gây mưa lớn Nhìn chung, mưa - lũ lưu vực nhiều loại hình thời tiết gây trận lũ lớn thường nhiều loại hình thời tiết hoạt động kết hợp với Diễn biến mưa lũ thường phức tạp, gây nhiều khó khăn cho cơng tác dự báo phịng chống Tùy theo vị trí đổ bão ATNĐ, tổ hợp chúng với KKL mà mưa lưu vực khác Mưa lưu vực thường lớn, cường độ cao vùng mưa lớn nước Đặc điểm riêng biệt trận lũ lụt lớn, lũ lụt lịch sử xảy có tác động tổ hợp hồn lưu bão, ATNĐ, đới gió đơng dải HTNĐ với KKL Lũ lụt lớn mưa bão, ATNĐ bão, ATNĐ kết hợp với KKL chiếm 61%; KKL kết hợp với dạng hoàn lưu khác - 39% tổng số trận Đặc biệt, điều kiện địa hình, quy luật chung khí hậu chế hồn lưu khí quyển, tháng nhiều bão miền Trung khu vực Quảng Nam – Đà Nẵng trùng với thời kỳ đầu hoạt động đợt khơng khí lạnh từ phía bắc tràn xuống Bão, ATNĐ dải HTNĐ có tác động KKL thường gây mưa to, diện rộng, cường độ mưa lớn, thời gian mưa kéo dài hơn, gây lũ lụt nghiêm trọng phạm vi rộng Bảng 5.2: Thống kê trận lũ từ mức báo động II trở lên (1976-2007) sông Thu Bồn - Trạm Câu Lâu Bắt đầu Kết thúc Hđỉnh DH T N G T N G cm cm 1976 11 13 08 11 15 17 314 237 KKL 102 1980 11 14 12 11 17 21 454 382 Bão vào Khánh Hoà (16/11) + KKL 369 1981 10 13 06 10 15 12 390 304 Bão vào Khánh Hoà (20/09) + KKL 356 TT Năm Hình thời tiết Xtblv mm - 10 26 13 10 29 07 383 274 ATNĐ + KKL 536 - 10 31 18 11 01 13 371 128 Bão vào Phú Yên (01/11) 173 - 11 07 13 11 09 12 365 273 KKL 201 1983 10 28 12 10 31 03 436 425 Bão vào Khánh Hoà (30/10) + KKL 503 11 14 13 11 17 05 325 285 KKL 275 1984 10 12 13 10 14 02 397 372 Bão vào Thuận Hải (12/10) + KKL 314 - Năm Kết thúc Hđỉnh DH T N G TT Bắt đầu T N G cm cm Hình thời tiết Xtblv mm 10 - 11 01 13 11 04 14 315 288 Bão vào Khánh Hoà (02/11) 292 11 - 11 07 19 11 08 20 316 219 Bão vào Bình Định (07/11) 114 12 1985 11 30 05 12 01 08 384 273 KKL 292 13 1986 10 20 07 10 23 05 330 272 Bão vào Quảng Ngãi (22/10) 126 12 02 07 12 04 06 419 330 Bão vào Thuận Hải (03/12) 296 1988 10 09 07 10 11 08 374 323 Bão vào Bình Định (24/07) 237 10 15 19 10 16 19 324 197 Bão vào Phú Yên (15/10) 123 1990 10 13 14 10 16 17 439 431 ATNĐ vào Phú Yên - K Hoà + KKL 320 11 12 01 11 13 21 437 342 Bão vào NhaTrang-Tuy Hoà (12/11) 272 19 1991 10 20 13 10 25 06 335 307 KKL 441 20 1992 10 21 13 10 25 06 395 363 Bão vào Tuy Hoà (23/10) 339 10 28 07 10 29 13 358 182 Bão vào Qui Nhơn (28/11) 144 22 1993 10 23 19 10 25 02 329 193 HTNĐ + Gió mùa đơng bắc 225 23 1994 10 18 07 10 21 12 355 367 KKL 248 24 1995 10 05 13 10 09 01 345 368 HTNĐ + KKL 423 14 15 16 17 18 21 - - - - 25 - 10 26 07 10 27 16 386 283 Bão vào Tp Qui Nhơn ( 26/10 ) 246 26 - 10 28 07 10 30 07 342 65 KKL 125 27 - 11 01 13 11 02 11 372 207 Bão vào Quảng Ngãi (01/11) 179 28 - 11 08 07 11 11 09 368 279 KKL 170 1996 10 21 16 10 23 09 363 294 ATNĐ vào Phú Khánh (22/10) 155 29 30 - 10 24 13 10 26 07 361 36 HTNĐ + KKL 163 31 - 10 27 19 10 30 07 420 168 HTNĐ + KKL 295 32 - 11 16 07 11 19 03 401 304 HTNĐ + KKL 406 33 - 11 30 13 12 01 20 351 181 ATNĐ + KKL 198 1997 09 20 07 09 23 01 381 389 HTNĐ + KKL 356 34 35 - 09 24 24 09 26 11 343 233 Bão vào Quảng Nam - Đà Nẵng (25/09) 183 36 - 11 02 07 11 04 17 366 298 HTNĐ + KKL 264 37 Kết thúc Hđỉnh DH T N G TT Bắt đầu T N G cm cm 1998 11 13 13 11 15 15 339 301 Bão vào N.Thuận - B.Thuận (14/11)+KKL 216 Năm Hình thời tiết Xtblv mm 38 - 11 19 07 11 21 05 509 416 Bão vào P.Yên - K.Hoà (20/11) + KKL 248 39 - 11 26 06 11 27 03 386 180 Bão vào Bình Định (26/11) 109 1999 11 01 12 11 03 13 523 452 KKL + HTNĐ 587 12 01 13 12 04 16 454 344 ATNĐ + KKL 488 2000 08 21 07 08 23 08 310 306 Bão vào Quảng Nam - Đà Nẵng (22/8) 179 40 41 42 - 43 - 10 08 16 10 11 04 370 383 Hai ATNĐ di chuyển dọc bờ biển 449 44 - 10 27 19 10 29 10 378 301 KKL + Đới gió đơng 157 45 - 11 16 07 11 18 04 383 202 KKL + Đới gió đơng 322 46 2001 10 21 15 10 22 18 415 174 ATNĐ vào Phú Khánh ( 21/10 ) 276 47 2003 10 15 07 10 16 18 324 362 HTNĐ + KKL 235 11 12 07 11 14 20 378 346 HTNĐ + KKL 207 49 2004 11 25 07 11 27 23 459 438 Bão + KKL 423 50 2005 10 10 11 10 11 17 345 292 Áp cao lục địa + HTND 166 48 - 51 - 10 23 10 10 25 11 356 334 HTNĐ + KKL 282 52 - 10 31 16 11 02 09 315 219 Bão dọc bờ biển Quảng Nam-Hà Tĩnh 207 2006 09 30 16 10 02 06 365 331 Bão vào Đà Nẵng (1/10) 234 12 05 07 12 06 19 364 317 Bão + KKL 250 2007 10 14 07 10 17 05 325 254 Bão vào đèo Ngang 272 53 54 55 - 56 - 10 30 10 10 31 07 373 323 ATNĐ vào Bình Định + KKL 186 57 - 11 01 19 11 02 24 382 228 KKL 115 58 - 11 10 07 11 12 11 539 476 Bão vào Khánh Hòa 530 59 - 11 17 07 11 19 12 356 209 Áp cao lục địa + đới gió đơng cao 145 1.2 Đặc điểm phân bố mưa hình thời tiết gây mưa lớn a) Bão Đây loại hình bão, ATNĐ ảnh hưởng gây mưa lớn mà khơng có đóng góp trực tiếp hệ thống synop khác Với loại hình phạm vi vùng mưa lớn thường hẹp phân bố tương đối đối xứng phụ thuộc vào cấu trúc XTNĐ mà thể phân bố trường mây ATNĐ đổ vào miền Trung mưa lớn thường ác liệt bão, ATNĐ đổ vào nơi khác, đồng thời địa hình chia cắt tạo nên phân hoá lớn phân bố mưa địa phương Đối với lưu vực sông Vu Gia - Thu Bồn, điều kiện hướng bờ biển địa hình dãy Trường Sơn mà tình hình mưa xảy khốc liệt Mưa bão thường xảy khoảng 2-4 ngày tùy thuộc tốc độ bão di chuyển nhanh hay chậm, hướng di chuyển, hệ thống synop xung quanh, khả tồn xoáy khu vực ảnh hưởng Cường độ XTNĐ có liên quan đến mưa trình phát triển đối lưu mạnh hay yếu nên khơng thể nói mưa bão lớn mưa ATNĐ Nhìn chung dạng hình thời tiết gây mưa không kéo dài Tổng lượng mưa bão, ATNĐ thường lớn tập trung - ngày, trung bình 300 - 400mm Lượng mưa bão 24h mức 150mm - 400mm, lượng mưa 200 - 300mm phổ biến chiếm gần 1/2 số bão, ATNĐ ảnh hưởng Rất có bão cho tổng lượng mưa 500mm 150mm Mưa sau bão, ATNĐ đổ thường kéo dài có cường độ lớn trước đổ Phân bố lượng mưa trận trung bình trạm lưu vực sông Vu Gia – Thu Bồn bão gây trình bày hình 5.2 Biểu đồ phân bố mưa theo không gian thời gian số trận mưa bão gây trình bày hình 5.1 - 5.8 - Phụ lục chương X (mm) Lượng mưa trung bình ảnh hưởng mưa ATNĐ bão 600 500 400 300 200 100 Hiên Khâm Thành Hội Đức Mỹ Khách Ái Trà My Nghĩa Hiệp Đức Tiên Phước Nông Sơn Giao Thủy Câu Lâu Hội An Cẩm Lệ Tam Kỳ Đà Trạm Nẵng Hình 5.2: Lượng mưa trung bình trạm hệ thống sơng Thu Bồn – Vu Gia bão gây ATNĐ nằm dải hội tụ nhiệt đới thường bền vững XTNĐ hoạt động đơn lẻ Đôi khi, xuất XTNĐ ảnh hưởng liên tiếp thời gian ngắn khu vực hẹp làm kéo dài thời gian mưa tổng lượng mưa, gây hậu nghiêm trọng Bão bão ATNĐ đổ liên tiếp vào miền Trung thường gặp loại thiên tai nguy hiểm thường gây lũ lụt lớn Từ năm 1964 đến nay, có 15 đợt bão bão ATNĐ đổ liên tiếp vào miền Trung, có lưu vực sơng Vu Gia–Thu Bồn Đặc điểm dạng thiên tai là: Mưa thường to, cường độ lớn, diện rộng, kéo dài nhiều ngày (7 – 10 ngày), tâm mưa đạt từ 800 - 1.000mm, 1.500mm có tác động KKL b) Khơng Khí lạnh ảnh hưởng độc lập Hình trình xâm nhập áp cao lạnh xuống phía nam khơng kết hợp với hệ thống thời tiết khác Mưa hình thời tiết thường khơng lớn chủ yếu tác động KKL tác động địa hình cịn ngun nhân khác nguyên nhân phụ Một đợt mưa lớn gió mùa đông bắc đơn gây thường kéo dài 23 ngày, có đợt kéo dài 4-5 ngày Tổng lượng mưa đợt trung bình từ 100200mm (hình 5.3) Hình 5.9 phụ lục chương thể phân bố lượng mưa theo khơng gian khơng khí lạnh đơn gây Như lượng mưa hình thời .Type = chChartTypeLine End With End If If mBD2 And mBD2 < max Then Set oSeries = oChart.SeriesCollection.Add With oSeries Caption = "BD2" SetData chDimCategories, chDataLiteral, X 'xvalues SetData chDimValues, chDataLiteral, BD2 'yvalues2 Type = chChartTypeLine End With End If If mBD3 And mBD3 < max Then Set oSeries = oChart.SeriesCollection.Add With oSeries Caption = "BD3" SetData chDimCategories, chDataLiteral, X 'xvalues SetData chDimValues, chDataLiteral, BD3 'yvalues2 Type = chChartTypeLine End With End If oChart.PlotArea.Interior.color = RGB(255, 255, 255) oChart.Axes(chAxisPositionLeft).MajorGridlines.Line.Weight = oChart.Axes(chAxisPositionLeft).MajorGridlines.Line.DashStyle = oChart.Axes(chAxisPositionLeft).MajorGridlines.Line.color = RGB(120, 120, 120) oChart.Axes(chAxisPositionBottom).MajorGridlines.Line.Weight = oChart.Axes(chAxisPositionBottom).MajorGridlines.Line.DashStyle = oChart.Axes(chAxisPositionBottom).MajorGridlines.Line.color = RGB(120, 120, 120) oChart.Axes(chAxisPositionBottom).HasMajorGridlines = True oChart.Axes(chAxisPositionLeft).HasMajorGridlines = True oChart.Axes(chAxisPositionLeft).Line.Weight = oChart.Axes(chAxisPositionBottom).Line.Weight = oChart.Axes(chAxisPositionLeft).NumberFormat = "#,##0" If Round((max - min) / 100) * 10 = Then oChart.Axes(chAxisPositionLeft).MajorUnit = Else oChart.Axes(chAxisPositionLeft).MajorUnit = Round((max - min) / 100) * 10 End If oChart.Axes(chAxisPositionLeft).Scaling.Minimum = - oChart.Axes(chAxisPositionLeft).MajorUnit oChart.Axes(chAxisPositionLeft).Scaling.Maximum = max + oChart.Axes(chAxisPositionLeft).MajorUnit oChart.Axes(chAxisPositionBottom).MajorUnit = Round(N / 10) oChart.Axes(chAxisPositionBottom).TickLabelSpacing = Round(N / 10) oChart.HasLegend = True oChart.legend.Position = chLegendPositionBottom oChart.legend.Font.name = ".vntime" oChart.legend.Font.Size = 12 oChart.Axes(chAxisPositionLeft).HasTitle = True If Left(Trim(dong1), 9) = "L-u l-ỵng" Then oChart.Axes(chAxisPositionLeft).Title.Caption = "L-u L-ỵng (m3/s)" Else oChart.Axes(chAxisPositionLeft).Title.Caption = "Mùc n-íc(cm)" End If oChart.Axes(chAxisPositionLeft).Title.Font.name = ".vntime" oChart.Axes(chAxisPositionLeft).Title.Font.Size = 12 oChart.Axes(chAxisPositionBottom).HasTitle = True oChart.Axes(chAxisPositionBottom).Title.Caption = "Thêi gian" oChart.Axes(chAxisPositionBottom).Title.Font.name = ".vntime" oChart.Axes(chAxisPositionBottom).Title.Font.Size = 12 ChartSpace1.Refresh End Sub Private Sub kanman_Click() Dim goi As Integer Open (App.path + "\Kalman\kalmaninput.txt") For Output As #1 For i = To Nod Print #1, Htd1(i), HDb1(i) Next i Close #1 Open (App.path + "\Kalman\ParaList.dat") For Output As #1 Print #1, "&IOConfig" Print #1, "InPath = '" + App.path + "\Kalman\'," Print #1, "OutPath = '" + App.path + "\Kalman\'," Print #1, " InFile = 'kalmaninput.txt'," Print #1, " OutFile = 'kalmankq.txt'" Print #1, "/" Print #1, "&KFConfig" Print #1, " Kal = 1," Print #1, " Choice = 0," Print #1, " id = 0" Print #1, "/" Close #1 Open (App.path + "\kalman.bat") For Output As #1 Print #1, "@ECHO OFF" Print #1, "ECHO." Print #1, "ECHO This is a batch file" Print #1, "ECHO." Print #1, "Cls" Print #1, "CD " + App.path + "\kalman" Print #1, "start KalmanFilter.exe" 'Print #1, "PAUSE" Print #1, "Cls" Print #1, "EXIT" Close Dim mofile As Integer Dim SS() As Single ReDim SS(Nod) As Single Dim tam As Single goi = Shell(App.path + "\kalman.bat", vbHide) 'mofile = Shell(App.path + "\Kalman\KalmanFilter.exe") MsgBox "Ket thuc tinh" Open (App.path + "\Kalman\kalmankq.txt") For Input As #1 For i = To Nod Input #1, tam, tam, SS(i), HDb1(i), tam Next i Close For i = Nod + To Nof HDb1(i) = HDb1(i) + SS(Nod) Next i For i = To Nof HDb(i) = HDb1(i) Next i Call vedothi Call danhgia End Sub Private Sub Picture1_Click() For i = To N HDb(i) = HDb(i) + Val(Text1.text) HDb1(i) = HDb(i) Next i Call vedothi Call danhgia End Sub Private Sub Picture2_Click() For i = To N HDb(i) = HDb(i) - Val(Text1.text) Next i Call vedothi Call danhgia End Sub Public Sub Hienthidata() 'Dim oSps As Spreadsheet Spreadsheet1.Caption = "Quá trình Thực đo Dù b¸o" Spreadsheet1.Windows(1).DisplayWorkbookTabs = False 'Spreadsheet1.Sheets(3).Delete 'Spreadsheet1.Sheets(2).Delete Spreadsheet1.ActiveSheet.Cells.Clear Spreadsheet1.ActiveSheet.Cells(1, 1) = "Thực đo" Spreadsheet1.ActiveSheet.Cells(1, 2) = "Dự báo" Spreadsheet1.ActiveSheet.Cells.Font.name = ".vnArial" Spreadsheet1.ActiveSheet.Cells.Font.Size = 10 Spreadsheet1.ActiveSheet.Cells.Font.color = RGB(0, 0, 255) Spreadsheet1.Windows(1).ViewableRange = "$a1:$b" & N & "" 'Spreadsheet1.Windows(1).DisplayHorizontalScrollBar = False For i = To N Spreadsheet1.ActiveSheet.Cells(i + 1, 1) = htd(i) Spreadsheet1.ActiveSheet.Cells(i + 1, 2) = HDb(i) Next i Spreadsheet1.ActiveSheet.Protection.AllowDeletingColumns = False Spreadsheet1.ActiveSheet.Protection.AllowDeletingRows = False Spreadsheet1.ActiveSheet.Protection.AllowFiltering = True Spreadsheet1.ActiveSheet.Refresh End Sub Private Sub Form_Load() Call DocSL l=0 Call thembien Call tinhchung(r()) End Sub Public Sub DocSL() Open App.path + "\data\" + Trim(Str(Myhydro_stnObj(ibasin).Code)) + ".HQSL" For Input As #1 Line Input #1, T1 Line Input #1, T2 Line Input #1, DVX Line Input #1, DVY FILE2 = App.path + "\data\" + "HOIQUY." + Trim(T2) FILE33 = App.path + "\data\" + "HS." + Trim(T2) Open FILE2 For Output As #2 Input #1, N, M ReDim X(N, M), X1(N, M), TB(M + 1), r(2 * M, * M), B(M + 1), XR(M, M) ReDim SSS(2 * M, * M), V(M + 1), C(M + 1), T(N), Y1(N) ReDim Y(N), GOD(N) ReDim TEN1(20) For j = To M Line Input #1, TEN1(j) For i = To N Input #1, X(i, j): X1(i, j) = X(i, j) ' If j = And X(i, j) Then MsgBox OK Next i ' Print " X(N+1,"; j; ") = "; X(n, j) Next j For j = To M 'TEN1(j) = Right$(" " + TEN1(j), 9) Next j Line Input #1, TEN1(0) For i = To N Input #1, X(i, 0) Y(i) = X(i, 0) Next i Close For i = To N GOD(i) = i Next i For j = To M s=0 For i = To N s = s + X(i, j) Next i TB(j) = s / N Next j For j = To M s=0 For i = To N X(i, j) = X(i, j) - TB(j) s = s + X(i, j) * X(i, j) Next i If j = Then SY = s V(j) = Sqr(s) For i = To N X1(i, j) = X(i, j) / V(j) Next i Next j Print #2, " So bien gia thiet M = "; M Print #2, " Do dai liet tinh toan N = "; N Print #2, "O -O" Print #2, "| Bien | trung binh | dltc |" Print #2, "| | | |" labelKQHQ.Caption = " So bien gia thiet M = " + Str(M) labelKQHQ.Caption = labelKQHQ.Caption + Chr(13) + " Do dai liet tinh toan N = " + Str(N) labelKQHQ.Caption = labelKQHQ.Caption + Chr(13) + "O -O" labelKQHQ.Caption = labelKQHQ.Caption + Chr(13) + "| Bien | trung binh | dltc |" labelKQHQ.Caption = labelKQHQ.Caption + Chr(13) + "| - | | |" labelKQHQ.Caption = labelKQHQ.Caption + Chr(13) For j = To M DR = V(j) / Sqr(N) Print #2, " X ("; j; ")"; labelKQHQ.Caption = labelKQHQ.Caption + " X (" + Str(j) + ")" Print #2, Tab(20); TB(j); Tab(35); V(j) / Sqr(N) labelKQHQ.Caption = labelKQHQ.Caption + " " + Str(TB(j)) + " " + Str(Round((V(j) / Sqr(N)), 2)) labelKQHQ.Caption = labelKQHQ.Caption + Chr(13) Next j Print #2, " Y"; Print #2, Tab(20); TB(0); Tab(35); V(0) / Sqr(N) Print #2, labelKQHQ.Caption = labelKQHQ.Caption + Chr(13) + " Y" labelKQHQ.Caption = labelKQHQ.Caption + " " + Str(TB(0)) + " Sqr(N)), 2)) labelKQHQ.Caption = labelKQHQ.Caption + Chr(13) so = Int((0.674 * V(0) / Sqr(N - 1)) * 100 + 0.5) / 100 For j = To M For H = To (j - 1) For i = To N r(j, H) = r(j, H) + X1(i, j) * X1(i, H) Next i Rr = 0.01 If Abs(r(j, H)) >= 0.1 Then Rr = r(j, H) r(H, j) = r(j, H) Rrtext = Left(Str(Int(Rr * 100 + 0.5) / 100) + " ", 6) 'CH = j: CY = j ' Da sua den day If (H = 0) And Abs(rm) < Abs(r(j, H)) Then rm = r(j, H): CH = j: CY = j End If Next H r(j, j) = Next j r(0, 0) = Print #2, " MA TRAN TUONG QUAN " Print #2, " - " " + Str(Round((V(0) / labelKQHQ.Caption = labelKQHQ.Caption + Chr(13) + " MA TRAN TUONG QUAN " labelKQHQ.Caption = labelKQHQ.Caption + Chr(13) + " " labelKQHQ.Caption = labelKQHQ.Caption + Chr(13) + " " + TEN1(0) 'Left(TEN1(0), 9) For j = To M labelKQHQ.Caption = labelKQHQ.Caption + " " + TEN1(j) 'Left(TEN1(j), 9) Next j labelKQHQ.Caption = labelKQHQ.Caption + Chr(13) Print #2, TEN1(0); 'Left(TEN1(0), 9); For j = To M Print #2, " " + TEN1(j); ' Left(TEN1(j), 9); Next j Print #2, Print #2, For j = To M labelKQHQ.Caption = labelKQHQ.Caption + TEN1(j) ' Left(TEN1(j), 9) Print #2, TEN1(j) 'Left(TEN1(j), 9); For H = To M Print #2, " " + Str(Round(r(H, j), 2)); labelKQHQ.Caption = labelKQHQ.Caption + " " + Left(Str(Round(r(H, j), 2)) + " 10) Next H Print #2, labelKQHQ.Caption = labelKQHQ.Caption + Chr(13) Next j ' Tao ma tran mo rong ' -r(0, 0) = For j = To M ", For H = M + To * M r(j, H) = 0: If H = j + M Then r(j, H) = Next H Next j For j = M + To * M For H = To * M r(j, H) = 0: If j = H + M Then r(j, H) = -1 Next H Next j ' Thu tu bien ' -For j = To M C(j) = j Next j End Sub Public Sub tinhchung(r() As Single) s = r(0, 0) * SY / (N - l - 1) s = Sqr(s) For j = To l jj = C(j) B(j) = r(jj, 0) * V(0) / V(jj) Next Bo = For i = To N s=0 For j = To l jj = C(j): s = s + B(j) * (X(i, jj) + TB(jj)) Next j Y1(i) = s Bo = Bo + (Y(i) - s) / N Next i NO = For i = To N Y1(i) = Y1(i) + Bo T(i) = Y(i) - Y1(i) If Abs(T(i)) - so Then r(0, 0) = R1 = Sqr(1 - r(0, 0)) Print #2, "He so tuong quan boi : "; Left(Str(R1), 5) Print #2, "H e s o x a c d i n h : "; Left(Str(R1 * R1), 5) Print #2, Print #2, "O -O" Print #2, "| PTHQ | DLTC | THONG KE T | TEN BIEN Print #2, "| | -| -| |" Print #2, " Y = "; Bo; "+" |" labelKQHQ.Caption = labelKQHQ.Caption + Chr(13) labelKQHQ.Caption = labelKQHQ.Caption + Chr(13) + "He so tuong quan boi : " + Left(Str(R1), 5) labelKQHQ.Caption = labelKQHQ.Caption + Chr(13) + "H e s o x a c d i n h : " + Left(Str(R1 * R1), 5) labelKQHQ.Caption = labelKQHQ.Caption + Chr(13) labelKQHQ.Caption = labelKQHQ.Caption + Chr(13) + "O -O" labelKQHQ.Caption = labelKQHQ.Caption + Chr(13) + "O PTHQ | DLTC | THONG KE T | TEN BIEN O" labelKQHQ.Caption = labelKQHQ.Caption + Chr(13) + "O -| -| -| O" labelKQHQ.Caption = labelKQHQ.Caption + Chr(13) + " Y = " + Str(Bo) + "+" labelKQHQ.Caption = labelKQHQ.Caption + Chr(13) Dim Z As String Z=" " For j = To l jj = C(j): B(j) = Int(B(j) * 100000 + 0.5) / 100000 Print #2, B(j); "X("; jj; ")"; labelKQHQ.Caption = labelKQHQ.Caption + Chr(13) + Str(B(j)) + "X(" + Str(jj) + ")" DL = Sqr(r(jj + M, jj + M)) * s / V(jj) If DL < 0.001 Then DL = 0.001 TT = B(j) / DL: DL = Int(DL * 1000 + 0.5) / 1000 labelKQHQ.Caption = labelKQHQ.Caption + " " + Str(DL) + " " + Left(Str(TT) + " 12) labelKQHQ.Caption = labelKQHQ.Caption + " " + Left(TEN1(jj) + Z, 20) Print #2, Tab(18); DL; Tab(20); Left(Str(TT) + " ", 12); Print #2, " " + Left(TEN1(jj) + Z, 20) labelKQHQ.Caption = labelKQHQ.Caption + Chr(13) Next j End Sub Private Sub Text1_Click() Command1.Visible = False Command2.Visible = False End Sub Private Sub Text1_KeyPress(KeyAscii As Integer) If (KeyAscii = 13) Then If Not IsNumeric(Me.Text1.text) Then MsgBox "not is numric, try againt" Exit Sub Else Command1.Visible = True Command2.Visible = True CH = CLng(CVar(Text1.text)) End If End If Me.labelKQHQ.Caption = " " Me.labelKQHQ.Caption = Me.labelKQHQ.Caption + Chr(13) End Sub ", Đề tài Nghiên cứu ứng dụng mơ hình WETSPA HECRAS mơ phỏng, dự báo q trình lũ hệ thống sông Thu Bồn – Vu Gia Chuyên đề 32: Ứng dụng công nghệ để dự báo thử nghiệm q trình lũ năm 2008 lưu vực sơng Vũ Gia - Thu Bồn Sau xây dựng cơng nghệ dự báo dịng chảy lũ cho hệ thống sông Thu Bồn – Vu Gia đề tài tiến hành dự báo thử nghiệm trình lũ cho vị trí hệ thống sơng Thu Bồn-Vu Gia với thời gian dự kiến 24 cảnh báo lũ đến 36 Phòng Dự báo Thuỷ văn Trung bộ, Tây nguyên Nam - Trung tâm Dự báo KTTVTƯ song song với Phòng dự báo - Đài Khí tượng thủy văn khu vực Trung Trung Bộ Thời gian thử nghiệm từ 1/10/2008 đến 30/11/2008 Trình tự dự báo Dự báo thử nghiệm tiến hành theo trình tự sau: - Cập nhật số liệu mưa, mực nước lưu lượng thực đo đến h hàng ngày, kiểm tra số liệu thực đo - Cập nhật mưa số trị phát báo lúc h hàng ngày – Chọn phương án dự báo mưa Phân tích số liệu mưa thực đo dự báo, so sánh điều chỉnh mưa dự báo Hiệu chỉnh kết dự báo mưa theo kết dự báo mưa synop kết dự báo mưa tham khảo từ mơ hình Mỹ Nhật - Tính tốn dự báo cho vị trí trạm sơng theo hai mơ hình WETSPA HECRAS theo quy trình mục 6.3 - Hiệu chỉnh kết dự báo theo phương pháp - Lặp lại bước cho vị trí tính tốn dự báo từ thượng nguồn đến hạ lưu theo dẫn chương trình - Ra tin dự báo, kết xuất kết dự báo Kết dự báo thử nghiệm 24 Dự báo thử nghiệm q trình dịng chảy hệ thống sông Thu Bồn – Vu gia tiến hành hàng ngày vào lúc 10h, sử dụng kết dự báo mưa số trị 24 với thời đoạn mưa 6h phát báo vào lúc 7h ngày ngày, số liệu mưa dự báo cung cấp phòng Dự báo thời tiết số liệu dự báo mưa trang web Mỹ Nhật Năm 2008 năm có nhiều lũ hệ thống sơng Thu Bồn – Vu Gia Từ 1h/1/X/2008 đến 19h/30/XI có đợt lũ với mực nước đỉnh lũ lớn Câu lâu đạt 3,94m; Ái Nghĩa đạt 9,1m lớn mức báo động Đặc trưng trận lũ năm 2008 trình bày bảng Bảng 1: Đặc trưng trận lũ lớn từ 1h/1/X/2008 đến 19h/30/XI/2008 Từ ngày 1/X đến ngày 30/XI/2008 thực dự báo lũ hệ thống sông điều kiện tác nghiệp theo hai phương án Kết dự báo thử nghiệm gồm 61 tin Các tin lưu trữ thư mục kết chương trình (D:\wetspathubon\ketqua\bantin\bantinwetspa) Kết dự báo cho trạm phát báo trình bày hình đến Các tiêu chất lượng thống kê bảng đến Bảng 2: Kết dự báo tác nghiệp đỉnh lũ đến Ái Nghĩa mô hình WETSPA WETSPA +HECRAS (Từ 1/X đến 30/XI/ 2008) Hmax thực đo (m) 7.18 9.1 7.16 7.75 7.35 WETSPA WETSPA + HECRAS Thời gian xuất Hmax (m) ∆H (m) ∆τ(h) Hmax (m) ∆H (m) ∆τ(h) 5.12.X 19.17.X 10.28.X 8.20.XI 23.26.XI 6.87 8.72 6.82 7.56 7.22 -0.31 -0.38 -0.34 -0.19 -0.13 -6 -6 7.14 8.83 7.39 7.92 7.45 0.04 0.27 -0.23 -0.17 -0.1 -6 -6 Bảng 3: Kết dự báo tác nghiệp q trình lũ đến Ái Nghĩa mơ hình WETSPA WETSPA +HECRAS (Từ 1/X đến 30/XI/ 2008) Thời gian dự kiến 6h 12h 18h 24h Mức sai số WETSPA cho Số điểm P% phép(cm) đúng/Tổng số 23 53/61 86.9 83.6 37 51/61 50 50/61 82.0 78.7 64 48/61 WETSPA + HECRAS Số điểm P% đúng/Tổng số 52/61 85.2 82 50/61 49/61 80.3 47/61 77 Bảng 4: Kết dự báo tác nghiệp đỉnh lũ đến Giao Thủy mô hình WETSPA WETSPA +HECRAS (Từ 1/X đến 30/XI/ 2008) Hmax thực đo (m) 7.05 8.41 6.50 7.41 6.90 WETSPA WETSPA + HECRAS Thời gian xuất Hmax (m) ∆H (m) ∆τ(h) Hmax (m) ∆H (m) ∆τ(h) 5.12.X 18.17.X 08.28.X 09.20.XI 11.26.XI 6.81 8.2 6.38 7.61 6.69 0.24 0.21 0.12 -0.2 0.21 -6 -6 6.72 8.72 6.21 7.28 6.67 -0.33 0.31 -0.29 -0.13 -0.23 -6 -6 Bảng 5: Kết dự báo tác nghiệp q trình lũ đến Giao Thủy mơ hình WETSPA WETSPA +HECRAS (Từ 1/X đến 30/XI/ 2008) Mức sai số WETSPA Thời gian dự cho phép Số điểm kiến (cm) đúng/Tổng số 19 53/61 6h 25 50/61 12h 44 49/61 18h 60 47/61 24h P% 86.2 82 80.3 77 WETSPA + HECRAS Số điểm P% đúng/Tổng số 88.5 54/61 83.6 51/61 85.2 52/61 80.3 49/61 Bảng 6: Kết dự báo tác nghiệp đỉnh lũ đến Câu Lâu mơ hình WETSPA WETSPA +HECRAS (Từ 1/X đến 30/XI/ 2008) Hmax thực đo (m) 2.65 3.94 2.38 3.13 2.84 WETSPA WETSPA + HECRAS Thời gian xuất Hmax (m) ∆H (m) ∆τ(h) Hmax (m) ∆H (m) ∆τ(h) 9.12.X 01.18.X 13.28.X 13.20.XI 15.26.XI 2.53 3.61 2.2 2.9 2.66 0.12 0.33 0.18 0.23 0.18 -6 6 2.37 4.12 2.5 3.18 2.71 -0.28 0.18 0.12 0.05 -0.13 0 6 Bảng 7: Kết dự báo tác nghiệp q trình lũ đến Câu Lâu mơ hình WETSPA WETSPA +HECRAS (Từ 1/X đến 30/XI/ 2008) WETSPA Mức sai số Thời gian dự cho phép Số điểm kiến (cm) đúng/Tổng số 29 51/61 6h 36 52/61 12h 50 48/61 18h 64 46/61 24h P% 83.6 85.2 78.7 75.4 WETSPA + HECRAS Số điểm P% đúng/Tổng số 85.2 52/61 82 50/61 80.3 49/61 77 47/61 Bản nhận xét đánh giá kết thử nghiệm đài KTTV khu vực Trung Trung trình bày phụ lục chương Nhận xét chung: Kết dự báo dịng chảy lũ hệ thống sơng Thu Bồn – Vu Gia từ ngày 1/X đến ngày 30/XI/2008 cho thấy dự báo dịng chảy lũ hệ thống sơng Thu Bồn – Vu Gia mơ hình WETSPA WETSPA + HECRAS có chất lượng "Đạt", sai số dự báo đỉnh lũ nằm sai số cho phép; tỷ số s/σ mức 0.35 - 0.60, tương ứng với hệ số tương quan tính tốn thực đo mức 0.80 Mức đảm bảo dự báo h đạt từ 83 - 86 %; 12h, 18 h mức từ 78 - 82 %; 24 giảm với mức đảm bảo từ 75 - 80 % Mơ hình mơ tốt q trình thực tế trạm phát báo Q trình tính tốn tương đối sát với trình thực đo Thời gian xuất hịên đỉnh lũ sát với thực tế Một số trận lũ có thời gian xuất đỉnh chậm từ 3–6 Quá trình nước lên nước xuống phù hợp Đường nước lên có sai số lớn đường nước xuống - Kết dự báo đỉnh lũ tương đối tốt thường lệch phải (xuất muộn khoảng 4-6 h) Sai số lớn phương án 110cm (24 giờ), mức đảm bảo trung bình đạt 78,3 % Sai số lớn phương án 80cm (24 giờ), mức đảm bảo trung bình đạt 82,8 % Các trạm vùng hạ lưu hệ thống sông bị ảnh hưởng thủy triều Câu Lâu, Cẩm Lệ, Hội An kết mô phần chân lũ nước thấp chịu ảnh hưởng thủy triều chưa tốt Biên độ dao động mực nước ảnh hưởng triều so với thực đo cịn nhỏ Với kết này, cơng nghệ dự báo dịng chảy lũ hệ thống sơng Thu Bồn – Vu Gia sử dụng phương án dự báo Trung tâm Dự báo KTTV TƯ đài KTTV khu vực Trung trung Kết cảnh báo 36 Từ ngày 1/X đến ngày 30/XI/2008 thử nghiệm cảnh báo lũ hệ thống sông điều kiện tác nghiệp Quá trình thử nghiệm cảnh báo lũ hệ thống sông Thu Bồn – Vu gia tiến hành hàng ngày vào lúc 10h Kết dự báo thử nghiệm gồm 61 tin Do mơ hình số trị cung cấp định lượng mưa cho vị trí lưu vực 24 với thời đoạn mưa 6h phát báo vào lúc 7h Nên số liệu mưa dùng để cảnh báo lũ hệ thống sông chủ yếu dựa vào số liệu mưa dự báo cung cấp phòng Dự báo thời tiết 24 đến 48 giờ, xu hướng diễn biến mưa hình thời tiết Một nguồn dự báo mưa quan trọng tham khảo dự báo định lượng mưa trang web Mỹ Nhật Tuy nhiên, số liệu dự báo mưa đưa vài điểm, với lượng mưa dự báo lũy tiến bước thời gian rộng nên khó sử dụng q trình cảnh báo lũ cho hệ thống sông Qua thời gian theo dõi chất lượng dự báo mưa trang Web Mỹ Nhật cho thấy dự báo xu có đợt mưa lớn diễn tiến xu mưa chấp nhận Tuy nhiên, dự báo lượng mưa sai khác lớn Điều dẫn đến chất lượng cảnh báo lũ 36 cho hệ thống sơng có chất lượng khơng cao Bảng 8: Kết cảnh báo lũ 36 mơ hình WETSPA WETSPA +HECRAS (Từ 1/X đến 30/XI/ 2008) cho trạm phát báo hệ thống sơng Thu Bồn – Vu Gia Vị trí Ái Nghĩa Giao Thủy Câu Lâu WETSPA Mức sai số cho phép Số điểm P% (cm) đúng/Tổng số 80 80 80 32/61 37/61 34/61 52.5 60.6 55.7 WETSPA + HECRAS Số điểm đúng/Tổng số 32/61 37/61 34/61 P% 52.5 60.6 55.7 Kết cảnh báo trước 36 h cho trận lũ năm 2007 tương đối đạt yêu cầu Mực đảm bảo cảnh báo 36 cho Ái Nghĩa 52.5 %, Giao Thủy 60.6%, Câu Lâu 55.7 % Trong trận lũ lớn cảnh báo xuất lũ, mức báo động mực nước lũ 36 mức báo động đỉnh lũ vị trí hạ lưu Trong 36 giờ, sai số cảnh báo đỉnh lũ thường lớn, nhiều lên tới 40 – 50 % mức đảm bảo cảnh báo thường mức 40 - 50 % Kết cảnh báo với thời gian 36 phụ thuộc nhiều vào độ xác dự báo mưa kinh nghiệm dự báo viên -100 28/11/2008 26/11/2008 24/11/2008 22/11/2008 20/11/2008 18/11/2008 16/11/2008 400 14/11/2008 700 12/11/2008 800 10/11/2008 28/11/2008 26/11/2008 24/11/2008 21/11/2008 19/11/2008 17/11/2008 15/11/2008 12/11/2008 10/11/2008 08/11/2008 06/11/2008 03/11/2008 01/11/2008 30/10/2008 28/10/2008 25/10/2008 23/10/2008 21/10/2008 19/10/2008 16/10/2008 14/10/2008 12/10/2008 10/10/2008 07/10/2008 05/10/2008 03/10/2008 01/10/2008 H (cm) 800 08/11/2008 01 /1 03 /20 /1 08 0/ 05 20 /1 08 07 /20 /1 08 0/ 10 20 /1 08 12 /20 /1 08 0/ 14 20 /1 08 16 /20 /1 08 19 /20 /1 08 0/ 21 20 /1 08 23 /20 /1 08 0/ 25 20 /1 08 28 /20 /1 08 0/ 30 20 /1 08 01 /20 /1 08 03 /20 /1 08 06 /20 /1 08 08 /20 /1 08 1/ 10 20 /1 08 12 /20 /1 08 15 /20 /1 08 1/ 17 20 /1 08 19 /20 /1 08 1/ 21 20 /1 08 24 /20 /1 08 1/ 26 20 /1 08 28 /20 /1 08 1/ 20 08 H (cm) 900 06/11/2008 04/11/2008 02/11/2008 31/10/2008 29/10/2008 27/10/2008 25/10/2008 23/10/2008 21/10/2008 19/10/2008 17/10/2008 15/10/2008 13/10/2008 11/10/2008 09/10/2008 07/10/2008 05/10/2008 03/10/2008 01/10/2008 H (cm) 1000 Hình 6.1:Đường qúa trình mực nước thực đo dự báo trạm Ái Nghĩa Hthực đo Hdự báo 700 600 500 400 300 t Hình 6.2: Đường qúa trình mực nước thực đo dự báo trạm Giao Thuỷ 900 Hthực đo Hdự báo 600 500 400 300 200 100 t Hình 6.3: Đường qúa trình mực nước thực đo dự báo trạm Câu Lâu 500 Hthực đo Hdự báo 300 200 100 t -50 -100 t 28/11/2008 26/11/2008 24/11/2008 21/11/2008 19/11/2008 17/11/2008 300 15/11/2008 350 12/11/2008 28/11/2008 26/11/2008 24/11/2008 21/11/2008 19/11/2008 17/11/2008 15/11/2008 12/11/2008 10/11/2008 650 10/11/2008 750 08/11/2008 06/11/2008 03/11/2008 01/11/2008 30/10/2008 28/10/2008 25/10/2008 23/10/2008 21/10/2008 19/10/2008 16/10/2008 14/10/2008 12/10/2008 10/10/2008 25/11/2008 20/11/2008 15/11/2008 10/11/2008 05/11/2008 31/10/2008 26/10/2008 21/10/2008 16/10/2008 11/10/2008 06/10/2008 01/10/2008 H (cm) 800 08/11/2008 06/11/2008 03/11/2008 01/11/2008 30/10/2008 28/10/2008 25/10/2008 23/10/2008 21/10/2008 19/10/2008 16/10/2008 14/10/2008 12/10/2008 10/10/2008 07/10/2008 05/10/2008 03/10/2008 01/10/2008 H (cm) 300 07/10/2008 05/10/2008 03/10/2008 01/10/2008 H (cm) 1000 Hình 6.4:Đường qúa trình mực nước thực đo dự báo trạm Ái Nghĩa 900 Hthực đo Hdự báo 700 600 500 400 400 Hình 6.6: Đường qúa trình mực nước thực đo dự báo trạm Câu Lâu Hthực đo Hdự báo 250 200 150 100 50 t Hình 6.5: Đường qúa trình mực nước thực đo dự báo trạm Giao Thuỷ 950 850 Hthực đo Hdự báo 550 450 350 250 150 t ... thảm phủ lưu vực hệ thống sông Thu Bồn - Vu Gia Đề tài Nghiên cứu ứng dụng mơ hình WETSPA HECRAS mơ phỏng, dự báo q trình lũ hệ thống sơng Thu Bồn – Vu Gia Chuyên đề 4: Nghiên cứu, chuyển đổi... để chạy chương trình khác Đề tài Nghiên cứu ứng dụng mơ hình WETSPA HECRAS mơ phỏng, dự báo q trình lũ hệ thống sông Thu Bồn – Vu Gia Chuyên đề 3: Nghiên cứu, chuyển đổi đồ số sử dụng đất thành... tinh TRMM (Trận lũ từ ngày 17/11/2007 đến 19/11/2007) Đề tài Nghiên cứu ứng dụng mơ hình WETSPA HECRAS mơ phỏng, dự báo q trình lũ hệ thống sông Thu Bồn – Vu Gia Chuyên đề 2: Nghiên cứu, xử lý chuyển