BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP VÀ PTNT TRƯỜNG ĐẠI HỌC THUỶ LỢI LÊ XUÂN KHÁNH NGHIÊN CỨU DỰ BÁO LŨ PHỤC VỤ VẬN HÀNH HỒ CHỨA VÀ MÔ PHỎNG NGẬP LỤT VÙNG HẠ LƯU SÔNG DINH TỈNH QUẢNG BÌNH LUẬN VĂN THẠC SĨ Hà Nội – 2019 BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP VÀ PTNT TRƯỜNG ĐẠI HỌC THUỶ LỢI LÊ XUÂN KHÁNH NGHIÊN CỨU DỰ BÁO LŨ PHỤC VỤ VẬN HÀNH HỒ CHỨA VÀ MÔ PHỎNG NGẬP LỤT VÙNG HẠ DU SÔNG DINH TỈNH QUẢNG BÌNH Chuyên ngành: Thủy văn học Mã số: 181800148 LUẬN VĂN THẠC SĨ NGƯỜI HƯỚNG DẪN KHOA HỌC: TS Nguyễn Tiến Thành PGS.TS Trần Kim Châu Hà Nội, NĂM 2019 LÊ XUÂN KHÁNH LUẬN VĂN THẠC SĨ HÀ NỘI, NĂM 2019 LỜI CAM ĐOAN Tôi xin cam đoan cơng trình nghiên cứu riêng tơi Các số liệu, kết trình bày luận án trung thực chưa cơng bố cơng trình khác Tác giả luận án Lê Xuân Khánh i LỜI CẢM ƠN Luận văn hồn thành Bộ mơn Mơ hình tốn Dự báo khí tựợng thủy văn hướng dẫn khoa học PGS.TS Trần Kim Châu TS Nguyễn Tiến Thành Tác giả xin bày tỏ lòng biết ơn sâu sắc tới hai Thầy tận tình giúp đỡ hướng dẫn bước nghiên cứu hoàn thiện luận văn Tác giả luận án xin cảm ơn Trường Đại học Thủy lợi, Đài Khí tượng Thủy văn Khu vực Trung Trung Bộ, Đài Khí tượng thủy văn tỉnh Quảng Bình tận tình giúp đỡ, tạo điều kiện thuận lợi cho tác giả học tập nghiên cứu hoàn thành luận án Nhân dịp này, tác giả xin bày tỏ lời cảm ơn sâu sắc tới nhà khoa học, thầy giáo, Ban giám hiệu Trường Đại học Thủy Lợi, bạn bè đồng nghiệp tận tình giúp đỡ, chia sẻ động viên tác giả suốt trình học tập nghiên cứu thực luận án Cuối cùng, tác giả xin cảm ơn đến gia đình người thân gia đình, đặc biệt người vợ hiền động viên, khích lệ, tạo điều kiện suốt q trình học tập nghiên cứu hồn thành luận văn Tác giả luận văn Lê Xuân Khánh ii MỤC LỤC LỜI CAM ĐOAN……………….……………………………………………… i LỜI CẢM ƠN…… ……………………………………………………………………… … ii MỤC LỤC ……………………………………………………………………………… … iii DANH MỤC HÌNH VẼ …………………………………………………………….……… vi DANH MỤC BẢNG BIỂU……… ………………………………………………….………vii MỞ ĐẦU CHƯƠNG I: TỔNG QUAN NGHIÊN CỨU 1.1 Tổng quan nghiên cứu dự báo mưa lũ 1.1.1 Các nghiên cứu thế giới 1.1.1.1 Dự báo mưa 1.1.2.2 Dự báo lũ 11 1.2 Tổng quan nghiên cứu mô ngập lụt vận hành hồ chứa 12 1.2.1 Các nghiên cứu thế giới 12 1.3 Tổng quan lưu vực nghiên cứu 15 1.3.1 Đặc điểm tự nhiên 15 1.3.2 Đặc điểm khí tượng thủy văn 16 1.3.3 Mạng lưới sông ngòi 18 1.3.4 Đặc điểm địa hình, địa mạo 19 1.4 Những tồn tại phương pháp tiếp cận 20 CHƯƠNG II: CƠ SỞ KHOA HỌC DỰ BÁO LŨ VÀ MÔ PHỎNG NGẬP LỤT 22 2.1 Đặc điểm hình thế thời tiết gây mưa lũ lưu vực sông Dinh 22 2.2 Nghiên cứu sản phẩm mưa dự báo cho lưu vực sông Dinh 23 2.2.1 Các mô hình dự báo mưa số trị thế giới Việt Nam 23 2.2.2 Lựa chọn sản phẩm mưa dự báo phù hợp cho lưu vực sông Dinh 25 2.3 Đánh giá sản phẩm mưa dự báo cho lưu vực sơng Dinh tỉnh Quảng Bình 29 2.3.1 Đánh giá sản phẩm mưa dự báo 29 iii 2.3.1.1 Sai số trung bình (ME - Mean Error) 30 2.3.1.2 Sai số tuyệt đối trung bình (MAE-Mean Absolute Error) 30 2.3.1.3 Sai số bình phương trung bình MSE (Mean Square Error) 30 2.3.1.4 Sai số bình phương trung bình quân phương 31 2.3.1.5 Hệ số tương quan (Correlation coefficient) 31 3.1.1.6 Phương pháp đánh giá đối với dự báo pha 31 2.3.1.7 Chỉ số FBI (BS hay FBI - Bias score): 32 2.3.1.8 Xác suất phát hiện (Probability oƒ Detection - POD) 33 2.3.1.9 Tỷ phần phát hiện sai (False Alarms Ratio - FAR) 33 2.3.1.10 Điểm số thành công 33 2.3.1.11 Điểm số thành công hợp lý (Equitable Threat Score — ETS) 34 2.3.1.12 Điểm số so le (Odds Ratio) 34 2.3.1.13 Điểm số Hanssen & Kuipers (HK hay TSS) 34 3.1.2 Hiệu chỉnh sản phẩm mưa dự báo 37 2.3 Phương pháp dự báo lũ cho lưu vực sông Dinh 38 2.3.1 Thực trạng dự báo lũ cho lưu vực sông Dinh 38 2.3.2 Lựa chọn phương pháp dự báo lũ 39 2.3.2.1 Một số phương pháp dự báo lũ 39 1.3.2.2 Lựa chọn mơ hình mô 41 1.3.2.3 Giới thiệu mơ hình Mike NAM 41 1.3.2.4 Giới thiệu mơ hình Mike 11 46 1.3.2.5 Giới thiệu mơ hình Mike 21FM 50 1.3.2.6 Giới thiệu mơ hình Mike Flood 51 2.4 Phương pháp mô ngập lụt cho lưu vực sông Dinh 52 CHƯƠNG III: DỰ BÁO LŨ PHỤC VỤ VẬN HÀNH HỒ CHỨA VÀ ĐÁNH GIÁ ẢNH HƯỞNG CỦA HỒ CHỨA ĐẾN NGẬP LỤT VÙNG HẠ LƯU SƠNG DINH TỈNH QUẢNG BÌNH 53 3.1 Xây dựng phương án dự báo dòng chảy lũ 53 3.1.1 Mạng lưới trạm khí tượng thủy văn 53 iv 3.1.1.1 Trạm truyền thống 53 3.2.1.2 Trạm tự động 55 3.2.2 Thiết lập mô hình thủy văn dự báo dòng chảy lũ 56 3.2.2.1 Thiết lập mô hình thủy văn cho lưu vực Đồng Tâm 56 3.2.2.2 Thiết lập mô hình thủy văn cho lưu vực sông Dinh 58 3.2.3 Hiệu chỉnh kiểm định mô hình 58 3.3 Dự báo thử nghiệm 61 3.3.2 Lựa chọn trận lũ dự báo 61 3.3.3 Dự báo dòng chảy tới hồ phương án vận hành an toàn hồ chứa 62 3.4 Mô ngập lụt vùng hạ lưu sông Dinh 63 3.4.1 Thiết lập mô hình thủy lực mô ngập lụt hạ lưu sông Dinh 63 3.4.1.1 Thiết lập mô hình thủy lực Mike11 cho sông Dinh 63 3.4.1.2 Thiết lập mô hình 2D Mike 21 64 3.4.2 Hiệu chỉnh kiểm định mơ hình 64 3.4.2.1 Hiệu chỉnh mô hình 2D 65 3.4.2.2 Kiểm định mô hình 2D 67 3.5 Đánh giá ảnh hưởng của hồ chứa đến ngập lụt vùng hạ lưu sông Dinh 69 3.5.1 Khi chưa có hồ chứa 69 3.5.2 Khi hồ chứa được vận hành 69 3.5.2.1 Lựa chọn kịch bản 69 3.5.2.3 Kết quả tính toán kịch bản 71 KẾT LUẬN VÀ KIẾN NGHỊ 73 Kết luận 73 Kiến nghị 73 TÀI LIỆU THAM KHẢO 75 PHỤ LỤC 78 v DANH MỤC HÌNH VẼ Hình Sơ đồ khới nghiên cứu luận văn Hình 1.1: Bản đồ địa hình lưu vực sông Dinh 20 Hình 2.1: Sơ đồ biểu diễn số phục vụ đánh giá chất lượng dự báo mưa 32 Hình 2.2: Sơ đồ minh họa phương pháp nội suy song tuyến tính 35 Hình 2.3: Giao diện trang web dự báo mưa cho hồ Thác Ch́i 39 Hình 2.4: Cấu trúc mô hình NAM 42 Hình 2.5: Sơ đồ sai phân hữu hạn điểm ẩn Abbott 48 Hình 2.6: Sơ đồ sai phân điểm ẩn Abbott mặt phẳng x ~ t 48 Hình 3.1: Bản đồ mạng lưới trạm KTTV tỉnh Quảng Bình 53 Hình 3.2: Bản đồ mạng lưới trạm khí tượng thủy văn lưu vực sơng Dinh 56 Hình 3.3: Bản đồ lưu vực Đồng Tâm 57 Hình 3.4: Sớ hóa tiểu lưu vực cho mô hình thủy văn lưu vực Đồng Tâm 57 Hình 3.5: Sớ hóa tiểu lưu vực cho mô hình thủy văn lưu vực sông Dinh 58 Hình 3.6: So sánh đường trình lũ tính tốn thực đo 15-17/10/2013 59 Hình.3.7: So sánh đường trình lũ tính toán thực đo 13-15/10/2016 59 Hình 3.8: So sánh đường trình lũ tính toán thực đo 31/10-02/11/2016 60 Hình 3.9: So sánh đường trình lũ tính toán thực đo 14-16/9/2017 60 Hình 3.10: Biểu đồ dòng chảy lượng mưa thực đo hồ Thác Chuối 62 Hình 3.11: Kết quả dự báo thử dòng chảy hồ mô hình mike NAM 62 Hình 3.12: Sơ đồ thủy lực lưu vực sông Dinh 63 Hình 3.13: Tạo địa hình cho mô hình 2D 64 Hình 3.14: Minh họa thiết lập Mike flood 64 Hình 3.15: Bản đồ ngập lụt vùng hạ du sông Dinh 16-18/10/2010 67 Hình 3.16: Bản đồ ngập lụt vùng hạ du sông Dinh 01-02/10/1985 68 Hình 3.17: Bản đồ ngập lụt hạ du sông Dinh theo kịch bản xả lũ thiết kế 70 Hình 3.18: Bản đồ ngập lụt hệ thống sông Dinh ứng với kịch bản vỡ đập mưa lớn 72 Phụ lục hình 1: Đường tần suất mưa ngày max trạm Đồng Hới 78 vi Bảng 3.13: So sánh kết quả tính tốn mơ thực đo lũ 01-02/10/1985 Kí hiệu Tọa độ Cao trình vết lũ (mét) Độ ngập điều tra (mét) Độ ngập tính tốn (mét) Chênh lệch (mét) 3.45 1.60 1.98 -0.38 7.86 2.66 2.87 -0.21 VL_01 VL_02 X 667692.80 664739.99 Y 1941406.33 1941406.24 VL_03 666553.98 1941567.15 6.48 1.55 1.76 -0.21 VL_04 VL_05 664951.71 664528.32 1940411.40 1941057.01 7.52 1.84 2.81 -0.97 7.64 2.09 2.12 -0.03 VL_06 663628.82 1940083.36 VL_07 VL_08 662581.00 661862.08 1939660.08 1938189.26 7.81 9.69 1.91 2.31 2.17 2.45 -0.26 -0.14 VL_09 662327.03 1937151.82 10.5 11.35 1.94 1.91 1.98 2.15 -0.04 -0.24 VL_10 661787.29 1936876.61 11.84 3.21 3.45 -0.24 VL_11 VL_12 663173.65 659306.14 1940263.33 1935058.22 8.3 20.37 3.34 2.21 3.55 2.63 -0.21 -0.42 VL_13 VL_14 VL_15 659140.73 659549.59 659734.83 1932842.34 1932472.71 1932856.42 23.86 3.64 3.79 -0.15 23.9 23.49 3.13 2.07 3.54 2.34 -0.41 -0.27 Hình 3.16: Bản đồ ngập lụt vùng hạ du sông Dinh 01-02/10/1985 Từ kết mô ngập lụt số liệu vết lũ thu thập tính hệ số 68 Nash mơ hình 0.70, đối chiếu với tiêu xếp loại chất lượng mơ hình bảng 3.9 Nhận thấy kết mô ngập lụt cho trận lũ tháng 10 năm 1985 mức 3.5 Đánh giá ảnh hưởng hồ chứa đến ngập lụt vùng hạ lưu sông Dinh Để đánh giá ảnh hưởng hồ chứa đến ngập lụt vùng hạ du, cách sử dụng mơ hình mơ ngập lụt trường hợp chưa có hồ điều tiết có hồ tham vào q trình điều tiết 3.5.1 Khi chưa có hồ chứa Khi lưu vực chưa có hồ chứa có mưa lũ thượng nguồn sơng Dinh ảnh hưởng trực tiếp đến trình lũ lụt hạ du Cụ thể kết mô trận lụt 0102/10/2010 gây ngập lụt diện tích 1173.5ha, năm 1985 2249.8ha cho xã vùng hạ du bảng 3.10 bảng 3.13 Vậy chưa có hồ chứa lũ lụt xẩy tương đối nghiêm trọng cho xã phía hạ du hồ Thác Chuối 3.5.2 Khi hồ chứa vận hành Khi có hồ chứa đưa vào vận hành điều tiết lũ hồ chứa có tác dụng cắt lũ cho vùng hạ du, vấn đề ngập lụt vùng hạ du chắn sẽ giảm nhẹ khả ngập sâu diện tích ngập lụt Bên cạnh cơng tác điều tiết lũ hồ chứa gây lên mối nguy hiểm thường trực cho vùng hạ lưu hồ như: tạo trận lũ nhân tạo; gây lên tượng lũ trồng lũ đặc biệt tồi tệ hồ chứa bị vỡ thảm họa vỡ đập khơng thể tính tốn lường trước Để có sở đánh giá phần ảnh hưởng hồ chứa đến ngập lụt vùng hạ du sông Dinh, luận văn nghiên cứu số kịch cố hồ chứa Thác Chuối xảy 3.5.2.1 Lựa chọn kịch bản Đối với việc xây dựng quy trình vận hành hồ chứa quan tâm đến khả xả hồ mà chưa tính đến việc xả có gây ngập lụt hạ lưu hay không Do cần thiết phải xây dựng đồ ngập lụt trường hợp Trong bối cảnh biến đổi khí hậu xuất trận mưa, lũ đặc biệt lớn 69 nhiều khả xảy vận hành hồ chứa gặp cố cửa van Luận văn đưa kịch để tính tốn với tình khác để từ có biện pháp phòng tránh giảm nhẹ tới mức tối thiểu thiệt hại ngập lụt gây Kịch 1: Lũ đến hồ tương ứng với lũ thiết kế (1%), kết hợp lượng mưa lưu vực ứng với tần suất 1% lưu lượng ứng với lưu lượng xả thiết kế từ có lũ Tổng lượng xả thiết kế 1083.4m3/s Kịch 2: Lũ đến hồ tương ứng với lũ kiểm tra (0.2%), kết hợp với lượng mưa đệm khu vực cịn lại sơng Dinh 1%, tất cửa van bị kẹt hồ bị vỡ tượng tràn qua đỉnh đập 3.5.2.2 Kết quả tính toán kịch bản Khi lũ kiểm tra đến hồ, để trì mực nước hồ thấp 59.46m theo quy trình, hồ phải xả lũ với lưu lượng xả thiết kế (1083.4 m3/s) Trên sông Dinh ngồi dịng chảy xả lũ từ hồ Thác Chuối, cịn có lượng dịng chảy nhập lưu khu sinh mưa lớn lưu vực kết hợp, làm cho hạ lưu sông Dinh xuất lũ lớn, gây ngập lụt Hình 3.17: Bản đồ ngập lụt hạ du sông Dinh theo kịch bản xả lũ thiết kế 70 Bảng 3.14: Thống kê xã ngập ứng với kịch bản Diện tích ngập (ha) Ngập lụt Xã 0.5-1 m Nhân Trach Lý trạch 1-3 m 3-6 m Tổng >6 m 5.06 5.06 46.49 94 15.79 4.26 160.54 Đại Trạch 167.48 89.4 24.43 18.2 299.51 Nam trạch 117.13 174.31 63.03 15.93 370.4 Việt Trung 165.58 249.08 89.6 12.7 516.96 3.5.2.3 Kết quả tính toán kịch bản Khi lũ kiểm tra đến hồ chứa, nguyên nhân đó, điện, cành từ thượng lưu lũ kéo lấp van làm tràn xả lũ khơng cịn có tác dụng Mực nước hồ dâng tràn qua đỉnh đập gây nên tượng xói hình thành vết vỡ đỉnh đập Lưu lượng đỉnh lũ tràn qua vết vỡ lên đến 4000 m3/s (tham khảo kết tính tốn vỡ đập từ dự án “Xây dựng phương án phòng chống lũ lụt cho hạ du hồ Thác Chuối – Huyện Bố Trạch” Ban quản lý dự án đầu tư xây dựng cơng trình ngành Nơng nghiệp PTNT tỉnh Quảng Bình làm chủ đầu tư nghiệm thu) Tương tự kịch 1, dòng chảy vỡ đập kết hợp với dòng chảy lũ mưa tiểu lưu vực gây trận lụt lớn Đây trường hợp nguy hiểm vết vỡ lớn sinh lượng nước lớn tràn xuống hạ lưu, kết hợp với mưa lớn hạ lưu làm mực nước sông dâng cao gây ngập lụt diện rộng Trong trường hợp này, vùng đất thấp dọc sông Dinh bị ngâp Nặng Nơng Trường Việt Trung Đại Trạch có đến 1000ha bị ảnh hưởng Diện tích độ sâu ngập lụt mô xã thống kê bảng Bảng 3.15: Thống kê diện tích ngập của xã ứng với kịch bản vỡ đập mưa lớn Ngập lụt Xã Nhân Trach Lý trạch Đại Trạch Nam trạch Việt Trung 0.5-1 m 9.37 20.99 464.75 171.82 238.99 Diện tích ngập (ha) 1-3 m 3-6 m >6 m 18.41 18.73 84.29 514.81 5.81 69.12 21.32 185.38 202.82 136.94 7.11 71 Tổng 27.78 39.72 1069.66 447.64 585.86 Hình 3.18: Bản đồ ngập lụt hệ thớng sơng Dinh ứng với kịch bản vỡ đập mưa lớn Như với kết mô ngập lụt vùng hạ du sơng Dinh có hồ chứa Thác Chuối tác động, nhận thấy với trường hợp hồ chứa vận hành bình thường cửa xả hoạt động bình thường hồ có tác dụng cắt lũ cho hạ du tốt, trường hợp van xả lũ hồ gặp cố (kịch 2) gây vỡ đập thảm họa xảy đến với vùng hạ du sông Dinh, gây ngập lụt nghiêm trọng cho xã có địa hình trũng thấp xã Đại Trạch 72 KẾT LUẬN VÀ KIẾN NGHỊ Kết luận Luận văn nghiên cứu khái quát điều kiện khí hậu tự nhiên Quảng Bình nói chung lưu vực sơng Dinh nói riêng, nghiên cứu hình gây mưa lũ cho lưu vực sông Dinh Ban đầu luận văn tiếp cận, tìm hiểu, nghiên cứu đánh giá sơ phương pháp dự báo mưa, lũ mô ngập lụt giới nước Luận văn lựa chọn cải thiện sản phẩm dự báo mưa cho lưu vực sông Dinh, thông qua việc sử dụng số sai số thống kê Kết cho thấy sản phẩm dự báo mưa ban đầu thường cho sai số thiên nhỏ Tuy nhiên, sau cải thiện phương pháp tổ hợp có trọng số, chất lượng dự báo mưa tốt so với ban đầu Luận văn bước đầu tiếp cận, nghiên cứu xây dựng công cụ dự báo dòng chảy tới hồ Thác Chuối (Mike NAM) để phục vụ công tác vận hành hồ chứa Luận văn mô ngập lụt với kịch mưa lớn lưu vực kết hợp với xả lũ với lưu lượng thiết kế kịch mưa lớn lưu vực gây vỡ hồ Thác Chuối, từ 02 kịch thấy tác động hồ chứa tới hạ du vận hành hồ xảy cố Luận văn đạt mục đích nghiên cứu ban đầu đề ra: Nghiên cứu đặc điểm hình gây mưa lũ lưu vực sơng Dinh; xây dựng phương án dự báo lũ phục vụ vận hành an toàn hồ chứa Thác Chuối lưu vực sông Dinh; Đánh giá ảnh hư ng hồ chứa tới ngập lụt vùng hạ lưu sông Dinh Kiến nghị Trong khuân khổ luận văn thực số nghiên cứu định hướng (dự báo dòng chảy đến hồ, dự báo mưa lưu vực từ sản phẩm mơ hình số trị), mơ ngập lụt vùng hạ du tác động hồ chứa Tuy nhiên thời nghiên cứu luận văn nhiều hạn chế chưa thể khắc phục: Số lượng trạm đo mưa lưu vực ít, số liệu dung để tính tốn tiêu sai số mơ hình số liệu tự động, vậy nhiều có sai số, tương lai 73 có số liệu thực đo chi tiết cần đánh giá lại mức độ tin cậy sản phẩm dự báo mưa Phương án dự báo dự báo dòng chảy đến hồ mơ hình Mike NAM xây dựng, chưa có tài liệu thực đo dồng chảy hồ để đánh giá xác mức độ tin cậy mơ hình Trong tương lai hồ có tài liệu quan trắc dịng chảy hồ cần có nghiên cứu đánh giá lại Kết tính tốn dự báo dòng chảy phụ thuộc nhiều vào dự báo mưa, thời gian chảy truyền tỉ trọng gia nhập khu phụ thuộc vào tâm mưa Ngoài cần tham khảo số liệu mưa dự báo từ nguồn khác để có đánh giá độ tin cậy dự báo cần hiệu chỉnh mưa dự báo dự báo tác nghiệp Cần tiếp tục hiệu chỉnh thông số mơ hình, đảm bảo độ ổn định Dự báo thử nghiệm cho trận lũ Luận văn xây dựng 02 kịch tác động hồ chứa tới khả ngập lụt hạ du Do vậy, cơng trình nghiên cứu khác cần xây dựng thêm kịch ngập lụt vùng hạ du sông Dinh 74 TÀI LIỆU THAM KHẢO [1] Kalnay, "Ensemble forecasting at NCEP: the breeding method Mon Wea Rev.," American Meteorological Society, vol 125, no 1997, pp 3297-3318., 1995 [2] Zhu, "Initial perturbations based on the ensemble transform (ET) technique in the NCEP global operational forecast system," Dynamic Meteorology and Oceanography, vol 60, no 2016, pp 62-79, 2010 [3] Palmer et al, "Ensemble prediction systems have been part of the operational forecasting suite at ECMWF," in ECMWF, 1992 [4] Philippe Bougeault et al, “The THORPEX Interactive Grand Global Ensemble,” American Meteorological Society, số 2009, 2009 [5] Zhu, “NCEP Ensemble Forecast System,” 2009, Environmental Modeling Center, 2009 [6] Molteni et al, “A strategy for high‐resolution ensemble prediction,” Royal Meteological society, tập 127, số 576, pp 2069-2094, 2001 [7] Marsigli et al, “A strategy for high‐resolution ensemble prediction,” Royal Meteorological Society, tập 127, số 576, pp 2095-2115, 2001 [8] Montani et al, "Ensemble activities at ARPA-SIM: the COSMO-LEPS and COSMO-SREPS systems," 2003, Bologna, Italy, 2003 [9] C.Marsigli et al, "The COSMO-LEPS mesoscale ensemble system: validation of the methodology and verification," in https://hal.archives-ouvertes.fr/hal00302614, 2005 [10] B.T Gouweleeuw, "Flood Forecasting using medium-range probabilistic weather prediction," Hydrology and Earth system Sciences, vol 9(4), pp 365-380, 2005 [11] F.Pappenberger et al, "Medium range multi model weather forecast ensembles in flood forecasting," in ECMWF, 2008 [12] C.J.White, "Using subseasonal to seasonal (s2s) extreme rainfall forecasts for extended-range flood prediction in Australia," IAHS, vol 370, p 229–234, 2015 [13] Nguyễn Ngọc Thục, Lương Tuấn Minh, “Các hình synop gây mưa lớn miền Bắc Việt Nam,” Tuyển tập báo cáo khoa học tại Hội nghị Khoa học dự báo KTTV lần thứ III 1986-1990, 1990 [14] Nguyễn Ngọc Thục cộng sự, "Những nhân tố gây mưa lớn, lũ tỉnh miền Trung Tây Nguyên," Thuộc chương trình lũ miền Trung, Cục Dự báo 1994, 1994 [15] Lê Văn Ánh, "Mưa, bão, lũ thiệt hại chúng gây khu vực Bắc Trung Bộ," Tuyển tập Báo cáo NCKH năm 2004 Trung Tâm KTTV Quốc Gia., Hà Nội, 2004 [16] Nguyễn Đức Hậu, "Nghiên cứu loại hình thời tiết, khí tượng thủy văn, thiên tai bất thường miền Trung," Báo cáo tổng kết đề tài nhánh thuộc đề tài mã số ĐTĐL 2009/01, Hà Nội, 2011 75 [17] Lê Bắc Huỳnh, "Báo cáo thiên tai lũ - Dự án UNDP VIE/97/002 -," Báo cáo thiên tai lũ - Dự án UNDP VIE/97/002, Hà Nội, 2000 [18] Tổng cục KTTV, "Dự án Trung tâm Khí tượng Thủy văn Quốc gia," Tổng cục KTTV, Hà Nội, 2013 [19] Viện khoa học KTTV&BĐKH, "Điều tra, khảo sát, xây dựng đồ phân vùng nguy lũ quét khu vực Miền Trung, Tây Nguyên, xây dựng hệ thống thí điểm phục vụ cảnh báo cho địa phương có nguy cao xảy lũ quét phục vụ công tác quy hoạch, đạo điều hành phòng tránh thiên t," Viện khoa học KTTV&BĐKH, Hà Nội, 2015 [20] Nguyễn Văn Thắng, "Báo cáo tổng kết đề tài cấp Bộ Nghiên cứu công nghệ xác định lượng mưa, kết hợp số liệu đa, vệ tinh với số liệu đo mưa trạm phục vụ dự báo KTTV," Bộ TN&MT, Hà Nội, 2013 [21] Trần Đình Trọng, "Báo cáo tổng kết đề tài cấp BộNghiên cứu sở khoa học, lựa chọn áp dụng phƣ ơng pháp dự báo thời tiết hạn cực ngắn Việt Nam," Bộ TN&MT, Hà Nội, 2011 [22] Đặng Ngọc Tĩnh, "Báo cáo tổng kết đề tài cấp BộNghiên cứu ứng dụng số liệu vệ tinh, mƣa dự báo số trị kết hợp số liệu bề mặt dự báo hệ thống sông Hồng," Bộ TN&MT, Hà Nội, 2011 [23] Nguyễn Lan Châu, "Báo cáo tổng kết đề tài cấp Bộ Nghiên cứu xây dựng công nghệ dự báo lũ sơng Đà phục vụ điều tiết hồ Hồ Bình cơng tác phịng chống lũ lụt," Bộ TN&MT, Hà Nội, 2004 [24] Đ T Mai, "Báo cáo tổng kết đề tài cấp Bộ 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," Bộ TN&MT, Hà Nội, 2005 [25] Đăng Ngọc Tĩnh nnk, "Báo cáo tổng kết đề tài cấp Bộ Nghiên cứu ứng dụng số liệu vệ tinh, viễn thám, mưa số trị kết hợp số liệu bề mặt dự báo thủy văn," Bộ TN&MT, Hà Nội, 2010 [26] Đặng Thanh Mai nnk, "Báo cáo tổng kết đề tài cấp Bộ Nghiên cứu phát triển hệ thống phân tích, giám sát, cảnh báo dự báo thiên tai thủy văn cho hệ thống sông Ba," Bộ TN&MT, Hà Nội, 2011 [27] Vũ Đức Long, "Báo cáo tổng kết đề tài cấp Bộ Nghiên cứu xây dựng công nghệ cảnh báo, dự báo lũ cho sơng Quảng Bình, Quảng Trị," Bộ TN&MT, Hà Nội, 2012 [28] Bùi Đình Lập, "Báo cáo tổng kết đề tài Nghiên cứu xây dựng cơng nghệ dự báo dịng chảy lũ đến hồ chứa lớn hệ thống sông Hồng," Bộ TN&MT, Hà Nội, 2014 [29] Đặng Thanh Mai, "Nghiên cứu xây dựng công nghệ dự báo lũ cảnh báo ngập lụt cho sơng tỉnh Bình Định Khánh Hòa," Bộ TN&MT, Hà Nội, 2014 [30] Nguyễn Lan Châu, "Báo cáo tổng kết đề tài Nghiên cứu áp dụng để dự báo dịng chảy sơng Hồng sơng Thái Bình," Bộ TN&MT, Hà Nội, 2007 [31] Nguyễn Lan Châu, "Báo cáo tổng kết đề tài cấp Bộ Nghiên cứu xây dựng công nghệ dự báo lũ sông Đà phục vụ điều tiết hồ Hồ Bình cơng tác 76 [32] [33] [34] [35] phòng chống lũ lụt," Bộ TN&MT, Hà Nội, 2004 Nguyễn Viết Thi, "Xây dựng công nghệ dự báo trước ngày dòng chảy đến hồ chứa lớn hệ thống sông Đà sông Lô," Bộ TN&MT, Hà Nội, 2006 Vũ Minh Cát, "Báo cáo tổng kết Hợp tác nghiên cứu xây dựng công nghệ dự báo lũ trung hạn kết nối với công nghệ điều hành hệ thống cơng trình phịng chống lũ cho đồng sơng Hồng – sơng Thái Bình," Bộ NN&PTNT, Hà Nội, 2007 Hoàng Thanh Tùng, "Luận án Tiến sĩ Nghiên cứu dự báo mưa, lũ trung hạn cho vận hành hệ thống hồ chứa phòng lũ ứng dụng cho lưu vực sông Cả," Đại học Thủy lợi, Hà Nội, 2015 Nguyễn Anh Tuấn, "Báo cáo tổng kết đề tài Nghiên cứu xây dựng hệ thống xác định khách quan hình gây mƣa lớn điển hình cho khu vực Việt Nam," Bộ TN&MT, Hà Nội, 2015 77 PHỤ LỤC FFC 2008 © Nghiem Tien Lam ĐƯỜNG TẦN SUẤT MƯA NGÀY MAX TRẠM ĐỒNG HỚI 890 840 X1 max Đồng Hới TB=250.05, Cv=0.37, Cs=1.15 790 TB=250.05, Cv=0.37, Cs=1.15 740 690 640 Lượng mưa, X(mm) 590 540 490 440 390 340 290 240 190 140 90 0.01 0.1 10 20 30 40 50 60 70 80 90 99 99.9 Tần suất, P(%) © FFC 2008 Phụ lục hình 1: Đường tần suất mưa ngày max trạm Đồng Hới Phụ Lục 1: Bảng kết quả tính toán tần suất mưa ngày max trạm Đồng Hới Thứ tự Tần suất P(%) X mm Thời gian lặp lại (năm) 0.01 844.25 10000.000 0.10 693.26 1000.000 0.20 647.27 500.000 0.33 613.80 303.030 0.50 585.82 200.000 1.00 538.63 100.000 1.50 510.65 66.667 2.00 490.58 50.000 3.00 461.94 33.333 10 5.00 425.08 20.000 11 10.00 373.19 10.000 12 20.00 317.81 5.000 13 25.00 298.80 4.000 14 30.00 282.66 3.333 15 40.00 255.64 2.500 16 50.00 232.79 2.000 17 60.00 212.10 1.667 18 70.00 192.22 1.429 19 75.00 182.17 1.333 78 99.99 20 80.00 171.76 1.250 21 85.00 160.63 1.176 22 90.00 148.06 1.111 23 95.00 132.22 1.053 24 97.00 123.58 1.031 25 99.00 110.43 1.010 26 99.90 96.03 1.001 27 99.99 90.11 1.000 Phụ lục 2: Kết quả mô hình NAM lưu vực Dinh 16-18/10/2010 Thời gian 16/10/2010 07:00 16/10/2010 08:00 16/10/2010 09:00 16/10/2010 10:00 16/10/2010 11:00 16/10/2010 12:00 16/10/2010 13:00 16/10/2010 14:00 16/10/2010 15:00 16/10/2010 16:00 16/10/2010 17:00 16/10/2010 18:00 16/10/2010 19:00 16/10/2010 20:00 16/10/2010 21:00 16/10/2010 22:00 16/10/2010 23:00 17/10/2010 00:00 17/10/2010 01:00 17/10/2010 02:00 17/10/2010 03:00 17/10/2010 04:00 17/10/2010 05:00 17/10/2010 06:00 17/10/2010 07:00 17/10/2010 08:00 17/10/2010 09:00 17/10/2010 10:00 17/10/2010 11:00 17/10/2010 12:00 17/10/2010 13:00 BASIN KG1 45.07 44.314 45.126 42.884 43.69 62.13 93.642 111.587 124.944 137.355 143.209 143.533 156.328 186.544 229.298 294.076 395.869 448.864 406.445 341.284 279.644 228.104 186.874 154.533 129.161 109.119 93.144 80.285 69.833 61.256 54.155 BASIN KG2 7.521 7.394 7.529 7.155 7.288 10.35 15.576 18.552 20.77 22.834 23.813 23.875 26.002 31.013 38.1 48.831 65.686 74.501 67.54 56.786 46.584 38.037 31.188 25.808 21.583 18.243 15.578 13.432 11.687 10.254 9.067 79 BASIN BASIN BASIN KG3 CUARA THACCHUOI 13.82 12.428 67.381 13.587 12.219 66.151 13.835 12.442 68.085 13.147 11.823 64.019 13.392 12.043 66.016 19.037 17.119 104.433 28.675 25.784 172.022 34.161 30.717 208.994 38.245 34.391 233.046 42.041 37.806 253.276 43.834 39.421 258.277 43.937 39.516 252.431 47.846 43.033 273.845 57.068 51.328 331.128 70.11 63.06 412.61 89.855 80.822 536.177 120.859 108.714 728.453 137.033 123.273 801.809 124.171 111.714 685.976 104.352 93.892 544.296 85.571 76.999 424.581 69.846 62.853 332.267 57.252 51.522 262.837 47.364 42.625 210.963 39.601 35.639 171.815 33.465 30.117 141.848 28.571 25.713 118.57 24.631 22.167 100.23 21.426 19.283 85.59 18.796 16.915 73.76 16.618 14.954 64.096 17/10/2010 14:00 17/10/2010 15:00 17/10/2010 16:00 17/10/2010 17:00 17/10/2010 18:00 17/10/2010 19:00 17/10/2010 20:00 17/10/2010 21:00 17/10/2010 22:00 17/10/2010 23:00 18/10/2010 00:00 18/10/2010 01:00 18/10/2010 02:00 18/10/2010 03:00 18/10/2010 04:00 18/10/2010 05:00 18/10/2010 06:00 18/10/2010 07:00 18/10/2010 08:00 18/10/2010 09:00 18/10/2010 10:00 18/10/2010 11:00 18/10/2010 12:00 18/10/2010 13:00 18/10/2010 14:00 18/10/2010 15:00 18/10/2010 16:00 18/10/2010 17:00 18/10/2010 18:00 18/10/2010 19:00 48.227 43.238 39.01 35.401 32.302 29.625 27.299 25.269 23.487 21.917 20.527 19.292 18.202 17.61 17.042 16.497 15.974 15.474 14.995 14.536 14.125 13.763 13.464 13.186 12.912 12.642 12.377 12.134 12.81 16.041 8.076 7.242 6.535 5.931 5.412 4.964 4.575 4.235 3.937 3.674 3.441 3.234 3.052 2.953 2.857 2.766 2.679 2.595 2.514 2.437 2.368 2.308 2.257 2.211 2.164 2.119 2.075 2.034 2.143 2.677 14.798 13.267 11.969 10.862 9.91 9.088 8.373 7.75 7.202 6.72 6.293 5.913 5.58 5.397 5.223 5.055 4.895 4.741 4.593 4.452 4.326 4.214 4.123 4.037 3.952 3.869 3.788 3.713 3.915 4.897 13.317 11.939 10.771 9.773 8.917 8.177 7.534 6.972 6.479 6.045 5.661 5.319 5.019 4.855 4.697 4.547 4.402 4.264 4.131 4.004 3.89 3.79 3.707 3.63 3.554 3.479 3.406 3.338 3.52 4.403 56.12 49.476 43.895 39.17 35.141 31.682 28.696 26.102 23.837 21.85 20.099 18.549 17.172 16.332 15.635 14.969 14.331 13.721 13.137 12.58 12.084 11.663 11.333 11.037 10.743 10.45 10.159 9.894 11.788 18.36 Phụ lục 3: Kết quả mô hình NAM lưu vực Dinh 01-03/10/1985 Date Time 01/10/1985 13:00 01/10/1985 14:00 01/10/1985 15:00 01/10/1985 16:00 01/10/1985 17:00 01/10/1985 18:00 01/10/1985 19:00 01/10/1985 20:00 01/10/1985 21:00 01/10/1985 22:00 01/10/1985 23:00 02/10/1985 00:00 BASIN KG1 3.605 19.254 64.912 148.693 251.602 346.913 421.956 472.319 498.44 522.008 556.005 631.732 BASIN KG2 0.747 4.11 13.421 29.335 47.431 63.021 74.67 82.202 85.9 89.473 94.998 107.72 80 BASIN BASIN BASIN KG3 CUARA THACCHUOI 1.55 1.489 10.617 8.682 8.404 65.934 27.84 26.765 206.35 59.216 56.109 407.651 93.229 86.829 589.428 121.217 111.181 713.234 141.457 128.4 793.683 154.233 139.144 840.737 160.249 144.101 857.804 166.375 149.584 883.887 176.302 158.752 933.794 199.657 180.498 1060.585 02/10/1985 01:00 02/10/1985 02:00 02/10/1985 03:00 02/10/1985 04:00 02/10/1985 05:00 02/10/1985 06:00 02/10/1985 07:00 02/10/1985 08:00 02/10/1985 09:00 02/10/1985 10:00 02/10/1985 11:00 02/10/1985 12:00 02/10/1985 13:00 02/10/1985 14:00 02/10/1985 15:00 02/10/1985 16:00 02/10/1985 17:00 02/10/1985 18:00 02/10/1985 19:00 02/10/1985 20:00 02/10/1985 21:00 02/10/1985 22:00 02/10/1985 23:00 03/10/1985 00:00 03/10/1985 01:00 03/10/1985 02:00 03/10/1985 03:00 03/10/1985 04:00 03/10/1985 05:00 03/10/1985 06:00 03/10/1985 07:00 03/10/1985 08:00 03/10/1985 09:00 03/10/1985 10:00 03/10/1985 11:00 03/10/1985 12:00 03/10/1985 13:00 721.573 798.063 734.533 630.659 525.051 451.456 378.738 274.624 184.558 125.678 89.379 66.305 50.975 40.37 32.776 27.18 22.954 19.697 17.139 17.625 22.792 25.402 24.592 22.859 20.936 19.066 17.339 15.785 14.406 13.333 12.726 12.128 11.546 10.982 10.439 9.92 9.425 122.783 135.509 124.413 106.647 88.691 76.222 63.911 46.294 31.072 21.133 15.013 11.125 8.544 6.759 5.481 4.54 3.829 3.282 2.852 2.936 3.813 4.257 4.121 3.829 3.504 3.187 2.895 2.633 2.399 2.215 2.113 2.012 1.914 1.819 1.728 1.64 1.557 81 227.26 250.453 229.56 196.564 163.348 140.339 117.632 85.148 57.102 38.807 27.549 20.402 15.659 12.38 10.034 8.305 7.001 5.995 5.206 5.365 6.984 7.806 7.558 7.019 6.42 5.837 5.299 4.814 4.385 4.041 3.853 3.669 3.488 3.314 3.146 2.986 2.833 205.709 226.49 205.746 175.27 145.079 124.735 104.366 74.835 49.592 33.373 23.528 17.344 13.271 10.472 8.478 7.014 5.912 5.065 4.402 4.601 6.121 6.93 6.733 6.252 5.708 5.177 4.688 4.249 3.862 3.585 3.412 3.243 3.078 2.918 2.766 2.62 2.481 1194.917 1291.954 1136.238 949.063 774.439 663.217 550.557 386.553 249.652 164.108 113.402 82.109 61.767 47.922 38.134 30.995 25.649 21.557 18.364 19.615 28.397 33.373 32.54 29.969 27 24.097 21.431 19.057 16.975 15.165 13.703 12.927 12.171 11.44 10.738 10.066 9.427 82 ... Phương pháp mô ngập lụt cho lưu vực sông Dinh 52 CHƯƠNG III: DỰ BÁO LŨ PHỤC VỤ VẬN HÀNH HỒ CHỨA VÀ ĐÁNH GIÁ ẢNH HƯỞNG CỦA HỒ CHỨA ĐẾN NGẬP LỤT VÙNG HẠ LƯU SÔNG DINH TỈNH QUẢNG BÌNH ... học dự báo lũ mô ngập lụt Chương 3: Dự báo lũ phục vụ vận hành hồ chứa đánh giá ảnh hưởng hồ chứa đến ngập lụt vùng hạ du lưu vực sông Dinh CHƯƠNG I: TỔNG QUAN NGHIÊN CỨU 1.1 Tổng quan nghiên cứu. .. vùng hạ lưu sông Dinh Cách tiếp cận phương pháp nghiên cứu: Hình Sơ đồ khối nghiên cứu luận văn Đối tượng nghiên cứu: Dòng chảy lũ đến hồ chứa, vận hành hồ chứa, ngập lụt hạ lưu hồ chứa vận hành