BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP VÀ PTNT TRƢỜNG ĐẠI HỌC THỦY LỢI DOÃN THỊ NỘI NGHIÊN CỨU SỰ BIẾN ĐỘNG CỦA MƢA LŨ VÀ ĐỀ XUẤT CƠ SỞ KHOA HỌC TÍNH LŨ CHO CÔNG TRÌNH GIAO THÔNG VÙNG NÚI ĐÔNG BẮC–VIỆT NAM LUẬN ÁN TIẾN SĨ KỸ THUẬT HÀ NỘI, NĂM 2016 BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP VÀ PTNT TRƢỜNG ĐẠI HỌC THỦY LỢI NGHIÊN CỨU SỰ BIẾN ĐỘNG CỦA MƢA LŨ VÀ ĐỀ XUẤT CƠ SỞ KHOA HỌC TÍNH LŨ CHO CÔNG TRÌNH GIAO THÔNG VÙNG NÚI ĐÔNG BẮC–VIỆT NAM Chuyên ngành: Thủy văn học Mã số: 62-44-02-24 NGƢỜI HƢỚNG DẪN KHOA HỌC PGS.TS Ngô Lê Long PGS.TS Hoàng Thanh Tùng HÀ NỘI, NĂM 2016 LỜI CAM ĐOAN Tác giả xin cam đoan công trình nghiên cứu thân tác giả Các kết nghiên cứu kết luận luận án trung thực, không chép từ nguồn dƣới hình thức nào.Việc tham khảo nguồn tài liệu (nếu có) đƣợc thực trích dẫn ghi nguồn tài liệu tham khảo quy định Tác giả luận án Doãn Thị Nội i LỜI CẢM ƠN Lời tác giả xin đƣợc bày tỏ lòng cảm ơn sâu sắc đến PGS.TS Ngô Lê Long, PGS.TS Hoàng Thanh Tùng tận tình hƣớng dẫn tác giả suốt thời gian nghiên cứu thực luận án Tác giả xin bày tỏ lòng biết ơn đến Ban Giám Hiệu, Phòng Đào tạo ĐH&SĐH, Tập thể Thầy cô giáo khoa Thủy văn Tài nguyên nƣớc, Phòng Khoa học Công nghệ, Trƣờng Đại Học Thủy Lợi - Hà Nội, giúp đỡ tạo điều kiện để tác giả hoàn thành luận án Tác giả xin chân thành cảm ơn Trƣờng Đại Học Giao Thông Vận Tải, Khoa Công Trình, Bộ môn Thủy Lực - Thủy Văn, nơi tác giả công tác, tạo điều kiện thời gian công việc giúp tác giả hoàn thành luận án Tác giả xin bày tỏ lòng biết ơn đến gia đình, bạn bè sát cánh động viên tác giả vƣợt qua khó khăn để thực luận án Tác giả luận án ii MỤC LỤC LỜI CAM ĐOAN .i LỜI CẢM ƠN ii MỤC LỤC iii DANH MỤC CÁC HÌNH ẢNH v DANH MỤC BẢNG BIỂU viii DANH MỤC CÁC TỪ VIẾT TẮT x MỞ ĐẦU 1 Tính cấp thiết đề tài .1 Mục tiêu nghiên cứu .2 Đối tƣợng phạm vi nghiên cứu Cách tiếp cận phƣơng pháp nghiên cứu Ý nghĩa khoa học thực tiễn Những đóng góp luận án Cấu trúc luận án .4 CHƢƠNG TỔNG QUAN VỀ NGHIÊN CỨU TÍNH LŨ THIẾT KẾ CHO CÔNG TRÌNH GIAO THÔNG 1.1 Tổng quan tính lũ thiết kế 1.1.1 Các nghiên cứu tính lũ thiết kế giới 1.1.2 Các nghiên cứu tính lũ thiết kế Việt Nam 1.2 Tổng quan tính lũ thiết kế cho công trình giao thông .8 1.2.1 Tính lũ thiết kế nƣớc Nhật, phƣơng Tây Mỹ 1.2.2 Tính lũ thiết kế nƣớc Đông Âu Nga 14 1.2.3 Tính lũ thiết kế Việt Nam 15 1.3 Những hạn chế tính lũ thiết kế cho giao thông Việt Nam 20 1.4 Đề xuất hƣớng tiếp cận phƣơng pháp nghiên cứu .22 1.5 Tổng quan khu vực nghiên cứu 24 1.5.1 Đặc điểm tự nhiên 24 1.5.2 Đặc điểm khí tƣợng thủy văn 26 1.5.3 Tình trạng giao thông cố công trình mùa mƣa lũ 30 1.5.4 Tình hình tài liệu nghiên cứu 36 iii 1.6 Kết luận chƣơng I 39 CHƢƠNG XÂY DỰNG CƠ SỞ KHOA HỌC TÍNH LŨ THIẾT KẾ CHO CÔNG TRÌNH GIAO THÔNG KHU VỰC NGHIÊN CỨU 41 2.1 Cơ sở lý thuyết phƣơng pháp tính lũ thiết kế 41 2.1.1 Phƣơng pháp Cơ quan bảo vệ thổ nhƣỡng Hoa Kỳ (SCS - CN) 41 2.1.2 Phƣơng pháp mô hình quan hệ 47 2.1.3 Phƣơng trình hồi quy .51 2.2 Cơ sở liệu phƣơng pháp tính lũ thiết kế .53 2.2.1 Xây dựng sở liệu mƣa 53 2.2.2 Xây dựng sở liệu mặt đệm 82 2.3 Kết luận chƣơng II 96 CHƢƠNG TÍNH TOÁN THỬ NGHIỆM VÀ ĐỀ XUẤT PHƢƠNG PHÁP TÍNH LŨ THIẾT KẾ CHO CÁC CÔNG TRÌNH GIAO THÔNG .98 3.1 Cơ sở phân nhóm công trình thoát nƣớc tính lũ thiết kế .98 3.2 Tính thử nghiệm theo phƣơng pháp khác 100 3.2.1 Thông số tính toán lƣu vực cầu tính thử nghiệm .100 3.2.2 Tính lũ thiết kế theo phƣơng pháp SCS - CN .108 3.2.3 Tính lũ thiết kế theo mô hình quan hệ 111 3.2.4 Tính lũ thiết kế theo phƣơng trình hồi quy tƣơng quan 112 3.2.5 Tính lũ theo pp Xokolopsky CĐGH (TCVN 9845:2013) 115 3.2.6 Đánh giá kết tính theo phƣơng pháp 117 3.3 Đề xuất phƣơng pháp tính phù hợp 119 3.4 Xây dựng chƣơng trình tính lũ thiết kế cho công trình giao thông 120 3.4.1 Giới thiệu chung chƣơng trình tính 120 3.4.2 Cấu trúc chƣơng trình tính 121 3.4.3 Hƣớng dẫn sử dụng chƣơng trình tính 123 3.5 Kêt luận chƣơng III .125 KẾT LUẬN VÀ KIẾN NGHỊ .126 DANH MỤC CÔNG TRÌNH ĐÃ CÔNG BỐ 130 TÀI LIỆU THAM KHẢO 131 PHỤ LỤC 136 iv DANH MỤC CÁC HÌNH ẢNH Hình 1.1 Các phƣơng pháp tính lũ cho giao thông giới .8 Hình 1.2 Các phƣơng pháp tính lũ cho giao thông Việt Nam 17 Hình 1.3 Sơ đồ tiếp cận nghiên cứu 23 Hình 1.4 Bản đồ khu vực nghiên cứu (vùng Đông Bắc) .24 Hình 1.5 Bản đồ địa hình khu vực nghiên cứu (Tỉnh Bắc Kạn Lạng Sơn) 25 Hình 1.6 Bản đồ tuyến đƣờng khu vực nghiên cứu 31 Hình 1.7 Bản đồ bố trí công trình thoát nƣớc khu vực nghiên cứu 32 Hình 1.8 Nƣớc chảy gây xói mái ta luy dƣơng rãnh dọc tuyến [39] 34 Hình 1.9 Tuyến đƣờng nội tỉnh Lạng Sơn bị ngập năm 2013 xói ta luy âm 34 Hình 1.10 Nƣớc lũ tràn mặt đƣờng 34 Hình 1.11 Cầu Sam Lang lúc khánh thành bị lũ trôi (sau hai tháng sử dụng) 35 Hình 1.12 Bản đồ trạm khí tƣợng khu vực nc (Tỉnh Bắc Kạn Lạng Sơn) 36 Hình 2.1 Đƣờng trình lũ đơn vị theo phƣơng pháp SCS 44 Hình 2.2 Sơ đồ tính lũ thiết kế theo phƣơng pháp SCS-CN 46 Hình 2.3 Sơ đồ tính lƣu lƣợng thiết kế theo phƣơng pháp mô hình quan hệ 51 Hình 2.4 Đƣờng trận bão năm 2012 năm 2013 .57 Hình 2.5 Đƣờng trận bão Utor năm 2013 57 Hình 2.6 Thành Phố Lạng Sơn ngập nƣớc lũ năm 2014 59 Hình 2.7 Đƣờng trận bão năm 2014 2015 .59 Hình 2.8 Biến đổi lƣợng mƣa ngày lớn trạm Bắc Kạn Bắc Sơn 62 Hình 2.9 Biến đổi lƣợng mƣa ngày lớn trạm Hữu Lũng Lạng Sơn 63 Hình 2.10 Biến đổi lƣợng mƣa ngày lớn trạm Ngân Sơn Thất Khê 63 Hình 2.11 Bản đồ hệ số biến thiên lƣợng mƣa ngày max (CV) vùng Đông Bắc 64 Hình 2.12 Sơ đồ xây dựng IDF theo hàm mũ .67 Hình 2.13 Tƣơng quan X1ng max-X1h max; X1ng max-X3hmax Bắc Kạn-Lạng Sơn (T=100) 68 Hình 2.14 Tƣơng quan X1ng max-X6h max; X1ng max-X12h max Bắc Kạn-Lạng Sơn(T=100) 69 Hình 2.15 Tƣơng quan X1ng max- X24 h max tỉnh Bắc Kạn-Lạng Sơn (T=100) 69 Hình 2.16 Tƣơng quan X1ng max-X1h max; X1ngay max-X3hmax Bắc Kạn-Lạng Sơn (T=50) 69 Hình 2.17 Tƣơng quan X1ng max-X6h max; X1ng max-X12hmax Bắc Kạn-Lạng Sơn (T=50) 70 Hình 2.18 Tƣơng quan X1ng max-X24h max Bắc Kạn-Lạng Sơn (T=50) 70 Hình 2.19 Sơ đồ bƣớc xây dựng đƣờng cong IDF (khu vực nghiên cứu) 75 Hình 2.20 Bộ đƣờng cong IDF (Lạng Sơn, T = năm T = 10 năm phút) .76 Hình 2.21 Bộ đƣờng cong IDF (Lạng Sơn, T = 25 năm T = 50 năm, phút) 76 Hình 2.22 Bộ đƣờng cong IDF (Lạng Sơn, T = 100 T = 200 năm, phút) 76 Hình 2.23 Bộ đƣờng cong IDF trạm Lạng Sơn Đình Lập 77 Hình 2.24 Bản đồ đẳng trị cƣờng độ mƣa (I-1-100) (Bắc Kạn Lạng Sơn) 78 Hình 2.25 Phân bố lũy tích mƣa 24h trạm Bắc Kạn Bắc Sơn 80 Hình 2.26 Phân bố lũy tích mƣa 24h trạm Chợ Rã Đình Lập 80 v Hình 2.27 Phân bố lũy tích mƣa 24h trạm Lạng Sơn Ngân Sơn 81 Hình 2.28 Phân bố lũy tích mƣa 24h trạm Thất Khê 81 Hình 2.29 Sơ đồ bƣớc xây dựng đồ số CN .84 Hình 2.30 Bản đồ loại đất tỉnh Bắc Kạn Lạng Sơn 86 Hình 2.31 Bản đồ trạng sử dụng đất tỉnh Bắc Kạn Lạng Sơn 88 Hình 2.32 Bản đồ số CN tỉnh Bắc Kạn Lạng Sơn 90 Hình 2.33 Sơ đồ bƣớc xây dựng đồ hệ số dòng chảy C 91 Hình 2.34 Bản đồ hệ số dòng chảy tỉnh Bắc Kạn Lạng Sơn (S > 6%) .92 Hình 2.35 Sơ đồ bƣớc xây dựng đồ hệ số nhám Manning n 93 Hình 2.36 Bản đồ hệ số nhám Manning tỉnh Bắc Kạn Lạng Sơn 94 Hình 2.37 Bản đồ lƣu vực cầu Bắc Khƣơng 95 Hình 3.1 Bản đồ lƣu vực cầu Can .100 Hình 3.2 Bản đồ lƣu vực cầu Bản Chắt .101 Hình 3.3 Bản đồ lƣu vực cầu Pắc Vằng 102 Hình 3.4 Bản đồ lƣu vực cầu Kỳ Lừa 103 Hình 3.5 Bản đồ trạng sử dụng đất hệ số nhám lƣu vực cầu Can 104 Hình 3.6 Bản đồ hệ số dòng chảy số CN lƣu vực cầu Can .104 Hình 3.7 Bản đồ trạng số CN lƣu vực cầu Bản Chắt 105 Hình 3.8 Bản đồ hệ số dòng chảy hệ số nhám lƣu vực cầu Bản Chắt 105 Hình 3.9 Bản đồ trạng số CN lƣu vực cầu Pắc Vằng .106 Hình 3.10 Bản đồ hệ số nhám hệ số dòng chảy C lƣu vực cầu Pắc Vằng 106 Hình 3.11 Bản đồ trạng sử dụng đất số CN lƣu vực cầu Kỳ Lừa .107 Hình 3.12 Bản đồ hệ số nhám hệ số dòng chảy lƣu vực cầu Kỳ Lừa 107 Hình 3.13 Quá trình lũ thực đo tính toán Trạm Lạng Sơn (2008 2013) .109 Hình 3.14 Quá trình mƣa lũ thiết kế lƣu vực cầu Kỳ Lừa cầu Pắc Vằng 110 Hình 3.15 Quá trình mƣa lũ thiết kế lƣu vực cầu Bản Chắt cầu Can 111 Hình 3.16 Các bƣớc tính theo phƣơng pháp CIA 112 Hình 3.17 Tƣơng quan Q100 ~A Q50~A Bắc Kạn - Lạng Sơn 114 Hình 3.18 Tƣơng quan Q25~A Q10 ~A Bắc Kạn - Lạng Sơn 114 Hình 3.19 Giao diện ban đầu ảnh vệ tinh Google map 121 Hình 3.20 Sơ đồ khối xây dựng chƣơng trình tính 122 Hình 3.21 Giao diện ban đầu chƣơng trình tính đồ 123 Hình 3.22 Kết tính lũ cầu Can theo pp CIA (Trên ảnh vệ tinh) 124 Hình Bộ đƣờng cong IDF trạm T = năm T = 10 năm 136 Hình Bộ đƣờng cong IDF trạm T = 25 năm T = 50 năm 136 Hình Bộ đƣờng cong IDF trạm T = 100 năm T = 200 năm 136 Hình Bộ đƣờng cong IDF trạm Thác Giềng Chợ Mới 137 Hình Bộ đƣờng cong IDF trạm Chợ Đồn An Tịnh 137 Hình Bộ đƣờng cong IDF trạm Bằng Khẩu Bằng Lũng 137 Hình Bộ đƣờng cong IDF trạm Bằng Phúc Cốc Đán 137 vi Hình Bộ đƣờng cong IDF trạm Côn Minh Đông Lạc 137 Hình Bộ đƣờng cong IDF trạm Dƣơng Phong Hảo Nghĩa 137 Hình 10 Bộ đƣờng cong IDF trạm Liên Thụy Nà Pạc 137 Hình 11 Bộ đƣờng cong IDF trạm Thuận Mang Xuân Dƣơng 137 Hình 12 Bộ đƣờng cong IDF trạm Xuân Lạc Yên Hán 137 Hình 13 Bộ đƣờng cong IDF trạm Yên Nhuận Yên Tịnh 137 Hình 14 Bộ đƣờng cong IDF trạm Vu Loan .137 Hình 15 Bản đồ lƣu vực cầu Khuổi Lu .137 Hình 16 Bản đồ phân vùng mƣa theo Đa giác Theissen cầu Can Bản Chắt 137 Hình 17 Bản đồ phân vùng mƣa theo Đa giác Theissen cầu Pắc Vằng Kỳ Lừa 137 Hình 18 Biến đổi lƣợng mƣa 1h 3h lớn trạm Thất Khê 137 Hình 19 Biến đổi lƣợng mƣa 6h 12h lớn trạm Thất Khê 137 Hình 20 Biến đổi lƣợng mƣa 1h 3h lớn trạm Ngân Sơn 137 Hình 21 Biến đổi lƣợng mƣa 6h 12h lớn trạm Ngân Sơn .137 Hình 22 Biến đổi lƣợng mƣa 1h 3h lớn trạm Lạng Sơn 137 Hình 23 Biến đổi lƣợng mƣa 6h 12h lớn trạm Lạng Sơn .137 Hình 24 Biến đổi lƣợng mƣa 1h 3h lớn trạm Hữu Lũng 137 Hình 25 Biến đổi lƣợng mƣa 6h 12h lớn trạm Hữu Lũng 137 Hình 26 Biến đổi lƣợng mƣa 1h 3h lớn trạm Đình Lập 137 Hình 27 Biến đổi lƣợng mƣa 6h 12h lớn trạm Đình Lập .137 Hình 28 Biến đổi lƣợng mƣa 1h 3h lớn trạm Chợ Rã 137 Hình 29 Biến đổi lƣợng mƣa 6h 12h lớn trạm Chợ Rã 137 Hình 30 Biến đổi lƣợng mƣa 1h 3h lớn trạm Bắc Sơn 137 Hình 31 Biến đổi lƣợng mƣa 6h 12h lớn trạm Bắc Sơn 137 Hình 32 Biến đổi lƣợng mƣa 1h 3h lớn trạm Bắc Kạn 137 Hình 33 Biến đổi lƣợng mƣa 6h 12h lớn trạm Bắc Kạn 137 Hình 0.34 Kết tính lũ cầu Pắc Vằng theo pp SCS (Trên ảnh vệ tinh) 137 Hình 0.35 Kết tính lũ cầu Pắc Vằng theo pp SCS (Trên ảnh vệ tinh) 137 vii DANH MỤC BẢNG BIỂU Bảng 1.1 Tóm tắt phƣơng pháp tính lũ thiết kế cho giao thông Mỹ 10 Bảng 1.2 Bảng thống kê phƣơng pháp tính lũ thiết kế Đông Âu Nga 15 Bảng 1.3 Quy định tần suất lũ 16 Bảng 1.4 Các đặc trƣng khí tƣợng trung bình nhiều năm .27 Bảng 1.5 Dòng chảy trung bình nhiều năm khu vực nghiên cứu 29 Bảng 1.6 Bảng phân loại cầu theo chiều dài khu vực nghiên cứu .33 Bảng 1.7 Các cố công trình giao thông khu vực nghiên cứu 35 Bảng 1.8 Các trạm đo mƣa thuộc tỉnh Bắc Kạn Lạng Sơn .37 Bảng 1.9 Các trạm quan trắc thủy văn thuộc tỉnh Bắc Kạn Lạng Sơn .38 Bảng 2.1 Xác suất bão, ATNĐ có lƣợng mƣa 24 cực đại cấp (%) 54 Bảng 2.2 Xác suất bão, ATNĐ có tổng lƣợng mƣa theo cấp (%) 54 Bảng 2.3 Kết kiểm định xu lƣợng mƣa ngày lớn 63 Bảng 2.4 Bảng hệ số CV lƣợng mƣa ngày lớn vùng Đông Bắc 65 Bảng 2.5 Bảng thông số trạm, tỉnh Bắc Kạn Lạng Sơn (T=100 năm) 68 Bảng 2.6 Bảng hệ số phƣơng trình đƣờng cong IDF, I = a.Dn .77 Bảng 2.7 Bảng hệ số phƣơng trình đƣờng cong IDF, I = a.Dn .77 Bảng 2.8 Kết tính cƣờng độ mƣa theo hàm mũ cho khu vực nghiên cứu 78 Bảng 2.9 Bảng so sánh kết tính theo luận án TCVN 9845:2013 .81 Bảng 2.10 Bảng phần trăm sai số kết tính theo luận án TCVN 82 Bảng 2.11 Bảng phân loại nhóm đất khu vực nghiên cứu 86 Bảng 2.12 Bảng thống kê trạng sử dụng đất khu vực nghiên cứu 87 Bảng 2.13 Bảng giá trị CN tình hình sử dụng đất nhóm đất 89 Bảng 2.14 Bảng tra hệ số dòng chảy C (nhóm A B) 91 Bảng 2.15 Bảng tra hệ số dòng chảy C (nhóm C D) 92 Bảng 2.16 Bảng tra hệ số nhám Manning .94 Bảng 2.17 Kết tính đặc trƣng (lƣu vực cầu Bắc Khƣơng) .95 Bảng 3.1 Các đặc trƣng lƣu vực (tính toán từ GIS) cầu Can 100 Bảng 3.2 Các đặc trƣng trƣng lƣu vực (tính toán từ GIS) cầu Bản Chắt .101 Bảng 3.3 Các đặc trƣng trƣng lƣu vực (tính toán từ GIS) Cầu Pắc Vằng 102 Bảng 3.4 Các đặc trƣng trƣng lƣu vực (tính toán từ GIS) Kỳ Lừa 103 Bảng 3.5 Các thông số CN, hệ số nhám hệ số dòng chảy cầu Can 104 Bảng 3.6 Các thông số CN, hệ số nhám hệ số dòng chảy cầu Bản Chắt 105 Bảng 3.7 Các thông số CN, hệ số nhám hệ số dòng chảy cầu Pắc Vằng .106 Bảng 3.8 Các thông số CN, hệ số nhám hệ số dòng chảy cầu Kỳ Lừa 108 Bảng 3.9 Kết tính lƣu lƣợng thiết kế theo phƣơng pháp SCS-CN .110 Bảng 3.10 Kết tính lƣu lƣợng thiết kế theo mô hình quan hệ .112 Bảng 3.11 Kết tính lƣu lƣợng thiết kế theo phƣơng trình hồi quy .115 Bảng 3.12 Kết tính lƣu lƣợng thiết kế Xokolopsky 116 viii Các lƣu vực thoát nƣớc tính thử nghiệm lƣu vực nghiên cứu Hình 15 Bản đồ lƣu vực cầu Khuổi Lu Hình 16 Bản đồ phân vùng mƣa theo Đa giác Theissen cầu Can Bản Chắt Hình 17 Bản đồ phân vùng mƣa theo Đa giác Theissen cầu Pắc Vằng Kỳ Lừa 142 Bảng Các công trình cầu đƣợc tính thử nghiệm luận án TT Tên Cầu A (km2) Quốc lộ Tỉnh Can 3.16 4B Lạng Sơn Dạ Tân 5.0 4B Lạng Sơn Khổng Thó 8.28 4B Lạng Sơn Khuổi Pẩu 6.29 4B Lạng Sơn Pắc Hán 1.21 4B Lạng Sơn Lâm Nghiệp 6.97 4B Lạng Sơn Khe Vuồng 4.2 31 Lạng Sơn Pò Mất 6.7 31 Lạng Sơn Còng Khƣa 7.6 31 Lạng Sơn 10 Pàn Mò 8.7 31 Lạng Sơn 11 Pàn Mò 9.8 31 Lạng Sơn 12 Nà Khƣơng 13.5 31 Lạng Sơn 13 Nà Phạ 3.3 31 Lạng Sơn 14 Nà Loong 42.6 31 Lạng Sơn 15 Pò Háng 294 31 Lạng Sơn 16 Toong soong 19,5 31 Lạng Sơn 17 Nà Cái 21,4 31 Lạng Sơn 18 Bản Chắt 25,24 31 Lạng Sơn 19 Khe Cháy 42,1 31 Lạng Sơn 20 Păc Lang 76,7 31 Lạng Sơn 21 Pò Tấu 116,7 31 Lạng Sơn 22 Treo Noi Thi 116,9 279 Bắc Kạn 23 Sau Hai 118,8 279 Bắc Kạn 24 Treo Na Luông 123,7 279 Lạng Sơn 25 Pắc Vằng 129,3 279 Lạng Sơn 26 Nà Lƣơng 174,6 279 Lạng Sơn 27 Păc Cù 261,9 279 Lạng Sơn 28 Hạo Nghĩa 281,3 279 Bắc Kạn 29 Nội Thị 331,2 279 Bắc Kạn Lạng Sơn 30 Na Loong 350,7 279 31 Chợ Mới 427,0 279 Bắc Kạn 32 Treo Hảo Nghĩa 453,6 279 Bắc Kạn 33 Khuổi Lu 505,2 279 Bắc Kạn 34 Bắc Khƣơng 506,5 279 Bắc Kạn 35 Xuất Hóa 776,4 279 Bắc Kạn 36 Cƣ Lễ 1041,2 279 Bắc Kạn 37 Kỳ Lừa 1559.9 T 238 Lạng Sơn 38 Bắc Sơn 52.5 31 Lạng Sơn 39 Đình Lập 31 Lạng Sơn 40 Đình Lập 31 Lạng Sơn 143 Kết đánh giá biến đổi lƣợng mƣa thời đoạn ngắn Hình 18 Biến đổi lƣợng mƣa 1h 3h lớn trạm Thất Khê Hình 19 Biến đổi lƣợng mƣa 6h 12h lớn trạm Thất Khê Hình 20 Biến đổi lƣợng mƣa 1h 3h lớn trạm Ngân Sơn 144 Hình 21 Biến đổi lƣợng mƣa 6h 12h lớn trạm Ngân Sơn Hình 22 Biến đổi lƣợng mƣa 1h 3h lớn trạm Lạng Sơn Hình 23 Biến đổi lƣợng mƣa 6h 12h lớn trạm Lạng Sơn 145 Hình 24 Biến đổi lƣợng mƣa 1h 3h lớn trạm Hữu Lũng Hình 25 Biến đổi lƣợng mƣa 6h 12h lớn trạm Hữu Lũng Hình 26 Biến đổi lƣợng mƣa 1h 3h lớn trạm Đình Lập 146 Hình 27 Biến đổi lƣợng mƣa 6h 12h lớn trạm Đình Lập Hình 28 Biến đổi lƣợng mƣa 1h 3h lớn trạm Chợ Rã Hình 29 Biến đổi lƣợng mƣa 6h 12h lớn trạm Chợ Rã 147 Hình 30 Biến đổi lƣợng mƣa 1h 3h lớn trạm Bắc Sơn Hình 31 Biến đổi lƣợng mƣa 6h 12h lớn trạm Bắc Sơn Hình 32 Biến đổi lƣợng mƣa 1h 3h lớn trạm Bắc Kạn 148 Hình 33 Biến đổi lƣợng mƣa 6h 12h lớn trạm Bắc Kạn Kết kiểm tra tiêu χ2 theo pp Gumbel Bảng Kết kiểm tra tiêu χ2 trạmLạng Sơn ( Gumbel) Khoảng cấp Xác suất P Ef Of Of-Ef (Of-Ef)2 (Of-Ef)2/Of Xn,max ≤ 0,10 5.5 ≥ 145 1.5 2.25 0.41 0,10