MỞ ĐẦU Tính cấp thiết đề tài Giao thông đƣợc coi huyết mạch quốc gia, muốn phát triển đất nƣớc cần phải hoàn thiện đại hệ thống giao thông Hàng năm nhà nƣớc đầu tƣ nhiều kinh phí cho phát triển giao thông toàn quốc Các tuyến đƣờng đƣợc đầu tƣ xây mới, nâng cấp mở rộng để đảm bảo giao thƣơng kinh tế văn hóa vùng nƣớc quốc tế Với 2/3 diện tích đất nƣớc địa hình đồi núi, theo thống kê Bộ GTVT đƣờng miền núi chiếm 70% km tổng km chiều dài đƣờng nƣớc Đƣờng miền núi chủ yếu cấp V, số km cấp III IV công trình thoát nƣớc nhƣ cầu, cống rãnh thoát nƣớc chƣa đƣợc trọng thiết kế xây dựng, việc tính toán thủy văn thủy lực hạn chế Trong năm gần đây, việc thiết kế đƣờng ô tô chuyển từ tiêu chuẩn 22 - TCN - 4054 - 85 sang 22 - TCN - 4054 - 98, 22TCN - 4054 - 2005 tiêu chuẩn tính toán lũ thoát nƣớc chuyển từ 22TCN - 1995, 22TCN 273 - 01; 22TCN 273 - 05 sang TCVN 9845:2013 với yêu cầu kỹ thuật cao hơn, tuyến đƣờng cần cải tạo theo tiêu chuẩn để tăng mức độ an toàn chạy xe nhằm đáp ứng chiến lƣợc an toàn giao thông quốc gia nâng cao hiệu vận tải Việt Nam quốc gia nằm khu vực khí hậu nhiệt đới gió mùa với chế độ mƣa lũ khắc nghiệt dƣờng nhƣ mức độ ngày tăng lên biến đổi khí hậu toàn cầu làm cho công trình giao thông thƣờng bị hƣ hỏng nặng nề Các hƣ hỏng hệ thống đƣờng có nhiều nguyên nhân nhƣ: chế độ khí hậu có thay đổi đáng kể xu ngày ác liệt; nguyên nhân địa chất, móng công tác xây dựng, vận hành bảo dƣỡng Tuy nhiên nhóm nguyên nhân kể vấn đề ảnh hƣởng mƣa, lũ nguyên nhân gây tác động đáng kể tới hƣ hỏng công trình Các hƣ hỏng thƣờng kể đến nhƣ tƣợng sạt lở ta luy dƣơng, âm, trôi cầu cống, hỏng mố trụ nguyên nhân tính lũ thiết kế chƣa chƣa phù hợp Trong TCVN 9845:2013,''Tiêu chuẩn tính toán đặc trƣng dòng chảy lũ'' (đƣợc biên soạn dựa theo QP.TL C - - 77) giới thiệu số phƣơng pháp tính lũ thiết kế từ mƣa rào nhƣ Cƣờng độ giới hạn, Xokolopsky, phƣơng pháp trƣờng Đại học Xây dựng Các phƣơng pháp tác giả Liên Xô (cũ) xây dựng đƣợc đƣa vào sử dụng nƣớc ta Tuy nhiên, công thức có nhiều thông số đƣợc xác định điều kiện nƣớc Nga, đƣa vào quy phạm tính toán Việt Nam dù đƣợc hiệu chỉnh nhƣng điều kiện thiếu số liệu nên chƣa đƣợc chuẩn hóa; có thông số khó xác định, với phạm vi thay đổi lớn, dẫn tới kết tính toán có độ xác không cao, tùy thuộc vào quan điểm lựa chọn thông số ngƣời sử dụng hậu rủi ro hƣ hỏng công trình tăng lên Khó khăn tính toán lũ thiết kế cho công trình giao thông Việt Nam thiếu tài liệu từ đặc trƣng lƣu vực, đất đai, lớp phủ, đặc biệt mƣa dòng chảy thực đo thời đoạn ngắn thƣờng có tài liệu mƣa ngày, nhiên khối lƣợng hạn chế Vì vậy, luận án“Nghiên cứu biến động mưa lũ đề xuất 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” mà NCS lựa chọn cấp thiết có ý nghĩa khoa học Kết luận án sở khoa học đề xuất phƣơng pháp tính lũ thiết kế đơn giản với độ tin cậy mức độ ổn định cao hơn, phục vụ xây dựng công trình thoát nƣớc đƣờng, góp phần cập nhật xây dựng qui trình tính hợp lý cho công tác tu, bảo dƣỡng nâng cấp công trình có, nhƣ xây dựng công trình an toàn, giảm thiểu thiệt hại cho công trình giao thông Mục tiêu nghiên cứu Nghiên cứu biến động mƣa lũ, chi tiết hóa mƣa mặt đệm, xác lập sở khoa học tính lũ thiết kế cho công trình giao thông Đối tƣợng phạm vi nghiên cứu - Phạm vi nghiên cứu luận án gồm hai tỉnh Bắc Kạn Lạng Sơn nằm vùng núi Đông Bắc - Việt Nam; - Đối tƣợng nghiên cứu mƣa lũ thiết kế phục vụ xây dựng công trình thoát nƣớc nằm quốc lộ QL3, 3B, 279, 3, 4A, 4B, 1A thuộc khu vực nghiên cứu Cách tiếp cận phƣơng pháp nghiên cứu Để đạt đƣợc mục tiêu đề ra, tác giả thu thập số liệu, tài liệu cần thiết, tiến hành nghiên cứu tổng quan biến động mƣa lũ (mƣa sinh lũ), phƣơng pháp tính lũ thiết kế cho công trình giao thông nƣớc giới từ lựa chọn hƣớng tiếp cận phù hợp, vừa mang tính kế thừa vừa đảm bảo tính sáng tạo nghiên cứu Các phƣơng pháp đƣợc sử dụng luận án bao gồm: i) phƣơng pháp phân tích, thống kê, kế thừa có chọn lọc tài liệu có nhằm tập hợp, đánh giá biến động mƣa lũ khu vực nghiên cứu; ii) phƣơng pháp phân tích ảnh viễn thám, GIS phục vụ mô lƣu vực mô hình toán xây dựng đồ chyên đề làm sở khoa học cho phƣơng pháp tính lũ kiến nghị; iii) phƣơng pháp mô hình toán, tính toán thử nghiệm làm sở cho việc kiến nghị phƣơng pháp tính lũ thiết kế phù hợp cho công trình giao thông khu vực Đông Bắc - Việt Nam Ý nghĩa khoa học thực tiễn Quá trình phát triển dân sinh kinh tế Việt Nam đòi hỏi phải nâng cấp xây hàng loạt tuyến đƣờng giao thông huyết mạch, đặc biệt tuyến đƣờng giao thông miền núi kết nghiên cứu luận án có ý nghĩa thực tiễn cao giải khó khăn tính toán thủy văn, thủy lực Việc nghiên cứu biến động mƣa lũ cho khu vực Đông Bắc xác lập sở khoa học tính lũ thiết kế cho công trình giao thông khu vực có ý nghĩa khoa học việc tiếp cận với phƣơng pháp tính toán đại tiện lợi làm tiền đề cho việc xây dựng quy trình tính toán phù hợp với điều kiện Việt Nam tƣơng lai gần Những đóng góp luận án - Luận án hoàn thiện phƣơng pháp tính lũ thiết kế cho công trình giao thông có xét đến biến động mƣa lũ chi tiết hóa đặc điểm tự nhiên khu vực Đông Bắc, sở ứng dụng thành tựu công nghệ đại hệ thống thông tin địa lý (GIS); - Đã bƣớc đầu xây dựng đƣợc phần mềm hỗ trợ tính lũ cho công trình thoát nƣớc đƣờng giao thông Cấu trúc luận án Ngoài phần Mở đầu, Kết luận kiến nghị, luận án đƣợc trình bày chƣơng: Chƣơng 1: Tổng quan nghiên cứu tính lũ thiết kế cho công trình giao thông Chƣơng 2: Xây dựng 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 Chƣơng 3: Tính thử nghiệm đề xuất phƣơng pháp tính lũ cho công trình giao thông vùng núi Đông Bắc - Việt Nam 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 Dòng chảy lũ đặc trƣng quan trọng tính toán thiết kế công trình, tính toán lũ vấn đề đặc biệt đƣợc quan tâm nghiên cứu Lũ thiết kế trận lũ đƣợc sử dụng thiết kế công trình có trạng thái bất lợi, độ lớn phụ thuộc vào cấp công trình đƣợc quy định quốc gia Lũ thiết kế bao gồm đỉnh lũ, tổng lƣợng trình lũ Hiện nay, tính lũ thiết kế phân thành hai nhóm: phƣơng pháp trực tiếp (phân tích thống kê xác suất) phƣơng pháp gián tiếp (phân tích mƣa mặt đệm) Việc tính lũ thiết kế trải qua trình dài nghiên cứu mang tính kế thừa phát triển nhằm xác đại hóa phục vụ xây dựng công trình an toàn mùa mƣa lũ đặc biệt điều kiện biến đổi khí hậu Một số nghiên cứu lũ thiết kế điển hình kể đến nhƣ: Chow (1964), Shaw [1] sổ tay tính toán thủy văn có đề cập đến phƣơng pháp tính toán lũ thiết kế phụ thuộc vào diện tích lƣu vực tình trạng số liệu: lƣu vực lớn, đủ số liệu dùng phƣơng pháp ngẫu nhiên (thống kê xác suất), lƣu vực nhỏ dùng phƣơng pháp mô hình quan hệ, đường lũ đơn vị quan hệ lưu lượng với diện tích thời gian Chow, Maidment (1988) [2] tài liệu có đề cập đến tính toán thủy văn đặc trƣng thủy văn thiết kế nhƣ trình thu phóng, lựa chọn mƣa thiết kế xây dựng đƣờng cong IDF, biểu đồ mƣa thiết kế dạng đƣờng cong tích lũy 24h, ƣớc tính thời gian mƣa giới hạn, tính toán lƣợng mƣa lớn khả (PMF), đồ đẳng trị mƣa với thời gian mƣa, D = - 60 phút hay 30 phút - 24h cho thời kỳ lặp lại T = - 100 năm Các phƣơng pháp chuyển đổi mƣa hiệu xác định dòng chảy thiết kế gồm đỉnh lũ, tổng lƣợng trình lũ thiết kế dùng để thiết kế công trình thoát nƣớc, mô vùng ngập lụt, thiết kế hồ chứa, sử dụng quản lý tài nguyên nƣớc Đối với thoát nƣớc, Chow giới thiệu phương pháp tính lũ cho lưu vực vừa nhỏ theo mô hình quan hệ với A diện tích lưu vực, I cường độ mưa, C hệ số dòng chảy Ngoài ra, đường lũ đơn vị đề cập sử dụng cho lưu vực vừa nhỏ Vijay (2002) [3] trình bày mô hình toán ứng dụng để tính lũ cho lƣu vực lớn lƣu vực nhỏ Đối với lƣu vực nhỏ mô hình ứng dụng trình bày 15 mô hình đại diện toàn giới Về lý thuyết để xây dựng mô hình kiến thức ứng dụng từ tài liệu Chow hay Maidment Raghunath (2006) [4] tài liệu nguyên lý thủy văn, trình bày vấn đề tính thủy văn vùng Tapti, Ấn độ (miền trung Ấn độ) Phần tính lũ thiết kế gồm tổng lƣợng lũ, đỉnh lũ, tần suất lũ, xác suất rủi ro với phƣơng pháp đề xuất nhƣ: Đƣờng lũ đơn vị tức thời, mô hình Nash, mô hình Clark, đƣờng lũ đơn vị SCS, hồi quy tuyến tính, phân tích thống kê xác suất, mô hình toán, tính lũ vị trí số liệu quan trắc theo phƣơng pháp hồi quy đa biến Ngoài tài liệu nêu, có nhiều tài liệu nghiên cứu liên quan đề cập đến phƣơng pháp tính lũ thiết kế giới Về bản, lý thuyết tập trung dòng chảy hay phƣơng thức chuyển đổi mƣa hiệu nhƣ tài liệu trên, nhiên từ hai thập kỷ trở lại với phát triển vƣợt bậc công nghệ máy tính, kỹ thuật viễn thám GIS cho phép nhà khoa học phân tích thử nghiệm, cập nhật công nghệ đại nhằm xác hóa tham số mà phƣơng pháp trƣớc chƣa xây dựng đƣợc 1.1.2 Các nghiên cứu tính lũ thiết kế Việt Nam Một số tài liệu điển hình đề cập đến tính lũ thiết kế Việt Nam nhƣ: Quy phạm QP.TL C-6-77 (1977) [5] trình bày phƣơng pháp tính toán đặc trƣng thủy văn cần thiết cho việc thiết kế công trình thủy lợi sông không bị ảnh hƣởng thủy triều Việt Nam Các đặc trƣng thủy văn đƣợc hƣớng dẫn tính quy phạm bao gồm: lƣu lƣợng bình quân năm, lƣu lƣợng lớn nhất, lƣu lƣợng nhỏ nhất, phân phối dòng chảy năm, loại mực nƣớc thiết kế thông số khác Khi tính lũ phục vụ thiết kế công trình sông trƣờng hợp đủ số liệu tiến hành phân tích tần suất, trƣờng hợp số liệu sử dụng công thức kinh nghiệm nhƣ: công thức cƣờng độ giới hạn cho lƣu vực có diện tích nhỏ 100 km2, công thức Triết giảm Xokolopxky cho lƣu vực có diện tích 100 km2 Cho đến quy định phƣơng pháp tính lũ thiết kế chủ yếu dựa quy phạm nên có số bất cập nhƣ bảng tra không đƣợc cập nhập; việc xác định số thông số phụ thuộc vào kinh nghiệm chủ quan ngƣời tính toán Đỗ Cao Đàm nnk (1990) [6] xuất Thủy văn công trình, tài liệu có trình bày cách tính lũ thiết kế, phƣơng pháp chủ yếu phƣơng pháp đƣợc đề cập QP.TL C - - 77 Lê Đình Thành (1997) [7] nghiên cứu tìm khả điều kiện ứng dụng phƣơng pháp tính mƣa lớn khả (PMP) lũ lớn khả (PMF), từ kiến nghị tiêu chuẩn tính lũ thiết kế hợp lý cho điều kiện Việt Nam Kết nghiên cứu đề cập cách chi tiết đến phƣơng pháp nhƣ tính lũ liên quan đến lũ lớn khả năng, nhiên công trình giao thông mức độ tiêu chuẩn nhƣ tính chất công trình xét theo toán cần phải có nghiên cứu cụ thể tƣơng lai Lê Văn Nghinh (2000) [8] biên soạn Nguyên lý Thủy văn, tài liệu quan trọng đề cập đến tính toán đặc trƣng thiết kế nhƣ dòng chảy năm, tháng, lũ, kiệt mực nƣớc thiết kế Các phƣơng pháp tính lũ thiết kế bao gồm phƣơng pháp nằm QP.TL C - - 77 Tuy nhiên, tài liệu chủ yếu đề cập đến dòng chảy tháng năm thiết kế phục vụ cho xây dựng vận hành hồ chứa Bộ môn TV&TNN (2003) [9] biên soạn Thủy văn thiết kế, tài liệu quan trọng dùng để tính toán đặc trƣng thiết kế công trình Tuy nhiên phƣơng pháp cách tiếp cận dựa QP.TL C - - 77 Phạm Ngọc Quý nnk (2005) [10], đề tài nghiên cứu khoa học cấp Bộ: “Nghiên cứu cảnh báo dự báo lũ vƣợt thiết kế - Giải pháp tràn cố” tiến hành xây dựng phần mềm tính lũ thiết kế Phần mềm cho phép tính lũ theo tần suất thiết kế dựa vào công thức kinh nghiệm QP.TL C - - 77 nêu trên, phƣơng pháp tính lũ đơn vị SCS, tính lũ lớn khả PMF theo phƣơng pháp thống kê Hasfield Tuy nhiên phần mềm chƣa có cập nhập bảng tra Hà Văn Khối nnk (2012) [11] cập nhật cho tái giáo trình Thủy văn công trình (ấn phẩm đƣợc xuất năm 1993) gồm tập Tập trình bày phƣơng pháp tính toán lũ thiết kế Về phƣơng pháp tính toán theo QP.TL C - - 77, nhiên giáo trình có cập nhập giới thiệu thêm kỹ thuật sử dụng tính toán lũ thiết kế nhƣ mô hình toán thủy văn bao gồm mô hình thủy văn tất định tính toán dòng chảy từ mƣa, mô hình lũ đơn vị Ngô Lê Long nnk (2015) [12] đề tài nghiên cứu khoa học cấp Nhà nƣớc “Nghiên cứu sở khoa học đề xuất tiêu chuẩn thiết kế lũ, đê biển điều kiện biến đổi khí hậu, nƣớc biển dâng Việt Nam giải pháp phòng tránh, giảm nhẹ thiệt hại” tiến hành nghiên cứu đề xuất phƣơng pháp tính lũ thiết kế cho công trình hồ chứa có xét tới tác động biến đổi khí hậu Nghiên cứu xác lập đƣợc sở khoa học thực tiễn tiêu chuẩn thiết kế lũ đƣợc đề xuất điều kiện biến đổi khí hậu, nƣớc biển dâng nƣớc ta đảm bảo an toàn, an sinh xã hội 1.2 Tổng quan tính lũ thiết kế cho công trình giao thông Việc tính lũ phục vụ thiết kế công trình giao thông giới đƣợc nhiều nhà khoa học quan tâm nghiên cứu Qua tìm hiểu phân tích tài liệu liên quan, phân phƣơng pháp tính lũ thiết kế cho công trình thoát nƣớc giao thông nói chung thành hai nhóm chính: i) nhóm phương pháp sử dụng nước Nhật, phương Tây Mỹ ii) nhóm phương pháp sử dụng nước Đông Âu, Nga Việt Nam (hình 1.1) Hình 1.1 Các phƣơng pháp tính lũ cho giao thông giới 1.2.1 Tính lũ thiết kế nước Nhật, phương Tây Mỹ +) Các phƣơng pháp tính lũ thiết kế Nhật Tính lũ thiết kế Nhật đƣợc đề cập nhiều tài liệu tính toán thủy văn hƣớng dẫn, quy phạm, quy chuẩn thiết kế công trình; tài liệu mang tính pháp lý đƣợc sử dụng nhiều cho ngành giao thông là: Hƣớng dẫn tính thủy văn thủy lực - Hƣớng dẫn tiêu chuẩn kỹ thuật cho dự án thiết kế (phần kiểm soát lũ) [13], [14] Nội dung tài liệu đề cập đến việc sử dụng phƣơng pháp mô hình quan hệ cho lƣu vực có diện tích A < 20km2 (dùng đƣờng cong IDF cho vùng có số liệu mƣa, trƣờng hợp số liệu mƣa sử dụng từ vùng có diện tích A > 100km2) Đối với lƣu vực A > 20km2 phƣơng pháp mô hình quan hệ sử dụng đƣờng lũ đơn vị phƣơng trình lƣợng trữ Tại công trình có số liệu lũ thực đo sử dụng phƣơng pháp thống kê xác suất theo lý thuyết Bulletin (1982) với yêu cầu tối thiểu n ≥ 10 năm đo đạc +) Các phƣơng pháp tính lũ thiết kế Anh Tính lũ thiết kế đƣợc trình bày Hƣớng dẫn thiết kế cầu đƣờng Tiêu chuẩn kỹ thuật Cơ quan đƣờng quốc gia [15], [16] nhiều tài liệu khác, phƣơng pháp phân theo diện tích: diện tích lƣu vực nhỏ (A < 20km2) không đủ số liệu đo đạc lũ sử dụng công thức đơn giản từ mƣa (mƣa năm), diện tích lƣu vực số đất; trƣờng hợp nhiều số liệu lũ thực đo tính theo phƣơng pháp thống kê xác suất +) Các phƣơng pháp tính lũ thiết kế Mỹ Các phƣơng pháp dùng tính lũ thiết kế giao thông Mỹ thƣờng đƣợc đề cập tài liệu nhƣ: [17] Hƣớng dẫn tính thoát nƣớc đƣờng (AASHTO) [18]; Tài liệu giới thiệu mô hình toán thủy văn HEC; Thủy văn đƣờng (FHWA) [19], [20]; Hƣớng dẫn kỹ thuật (TR55) hay tiêu chuẩn thiết kế [21], [22], [23], [24]; Nhìn chung, phƣơng pháp đƣợc chia thành hai nhóm: i) Đối vùng rộng lớn, có số liệu thực đo sử dụng phƣơng pháp thống kê xác suất với phân bố nhƣ Log Normal, PIII, Gumbell; ii) Đối với vùng số liệu dựa vào đặc tính vùng để tính theo phƣơng pháp nhƣ: mô hình quan hệ; đƣờng lũ đơn vị tổng hợp SCS; phƣơng trình hồi quy vùng hồi quy theo USGS (dùng quy hoạch) Ngoài ra, phƣơng pháp TR55 (mô hình WinTR55) thƣờng dùng để tính lũ cho lƣu vực nhỏ A(F) < 65km2 cho kết tốt Do điều kiện số liệu đầy đủ chi tiết mƣa mặt đệm nên Mỹ xây dựng nhiều bảng tra, đồ cƣờng độ mƣa thiết kế toàn quốc, lƣợng mƣa thời đoạn dài (d > 1h) đƣợc chuyển đổi thành lƣợng mƣa thời đoạn ngắn (d < 1h) Các phƣơng pháp đƣợc tóm tắt bảng sau: Bảng 1.1 Tóm tắt phƣơng pháp tính lũ thiết kế cho giao thông Mỹ Phƣơng pháp 1)Mô hình quan hệ Tóm tắt yêu cầu phƣơng pháp * Lƣu vực nhỏ (A < 1,3 km2) * Thời gian tập trung dòng chảy Tc < 1h Số liệu yêu cầu * Tc thời gian tập trung dòng chảy (h); * A Diện tích lƣu vực (km2) * Mƣa phân bố theo không gian thời gian * C hệ số dòng chảy * Dòng chảy tràn bề mặt chủ yếu * I cƣờng độ mƣa (mm/h) * Lƣợng trữ kênh không đáng kể 2) Phƣơng trình hồi * Diện tích lƣu vực giới hạn vùng * Diện tích lƣu vực A (km2) quy vùng USGS * Giá trị đỉnh lũ điều kiện tự nhiên * Mƣa trung bình năm P(mm) không chịu tác động yếu tố * Độ cao (độ dốc lƣu vực) khác Thƣờng dùng quy hoạch 3) NRCS * Lƣu vực vừa nhỏ A < 65 km2 * Mƣa 24h 4) TR55 * Tc = 0,1–10h * Phân bố mƣa dạng I, II, IA, III * Dòng chảy gồm chảy tràn chảy * Chỉ số CN kênh * Tc (h) * Diễn toán kênh đơn giản * A (km2) 5) Đƣờng lũ đơn vị * A = 0,4–2500 km2 SCS * Thời gian mƣa cƣờng độ mƣa * Bản đồ đẳng trị, phân bố mƣa * Diện tích A (km2) * Quan hệ mƣa dòng chảy tuyến * Chiều dài lƣu vực L (km; m) tính * Chiều dài lƣu vực Lc (km; m) * Biểu đồ lũ đơn vị tổng hợp 6) Phƣơng pháp thống kê: Log PIII; Bullentin 17B * Lƣu vực vừa lớn có trạm đo đạc, số *n > 10 năm liệu đầy đủ * Số liệu H, Q 7) Chuyển đổi * Lƣu vực tƣơng tự * Q, A1 (km2) A2 (km2) +) Các phƣơng pháp tính lũ thiết kế Úc Phƣơng pháp tính lũ thiết kế đƣợc đề cập 10 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