Luận án tiến sĩ kỹ thuật 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

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 Luận án tiến sĩ Kỹ Thuật 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 Hồng Thanh Tùng HÀ NỘI, NĂM 2016 Luận án tiến sĩ Kỹ Thuật 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 Luận án tiến sĩ Kỹ Thuật 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 Hồ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 giáo khoa Thủy văn Tài ngun 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ả hồn thành luận án Tác giả xin bày tỏ lịng biết ơn đến gia đình, bạn bè ln 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 Luận án tiến sĩ Kỹ Thuật 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 Luận án tiến sĩ Kỹ Thuật 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 TỐ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 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 tố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 Luận án tiến sĩ Kỹ Thuật 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 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 khơng có 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 Luận án tiến sĩ Kỹ Thuật 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 Q trình lũ thực đo tính tố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 Luận án tiến sĩ Kỹ Thuật 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 Luận án tiến sĩ Kỹ Thuật 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 tố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 tố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 Luận án tiến sĩ Kỹ Thuật Các lƣu vực 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 Luận án tiến sĩ Kỹ Thuật 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 Luận án tiến sĩ Kỹ Thuật 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 Luận án tiến sĩ Kỹ Thuật 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 Luận án tiến sĩ Kỹ Thuật 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 Luận án tiến sĩ Kỹ Thuật 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 Luận án tiến sĩ Kỹ Thuật 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 Luận án tiến sĩ Kỹ Thuật 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