đồ án tốt nghiệp kĩ sư thủy văn nghiên cứu dự báo mực nước đỉnh lũ tại trạm thủy văn diêu trì trên sông hà thanh tỉnh bình định bằng phương pháp hồi quy tuyến tính bội

đồ án nghiên cứu dự báo mực nước đỉnh lũ tại trạm thủy văn Diêu Trì trên sông Hà Thanh tỉnh Bình Định bằng phương pháp hồi quy tuyến tính bội. Đồ án tốt nghiệp kĩ sư thủy văn trường đại học tài nguyên và môi trường TP.HCM

LỜI CAM ĐOAN Tôi xin cam đoan công trình nghiên cứu riêng tơi hướng dẫn khoa học Ths Lê Văn Phùng Các nội dung nghiên cứu kết đồ án tốt nghiệp trung thực Những số liệu bảng biểu phục vụ cho phân tích, tính tốn đánh giá thu thập từ nguồn có ghi rõ phần tài liệu tham khảo Trong đồ án có sử dụng số đánh giá số liệu quan tổ chức khác có trích dẫn thích nguồn gốc Nếu phát có gian lận nào, tơi xin hồn tồn chịu trách nhiệm nội dung đồ án tốt nghiệp TP Hồ Chí Minh, ngày tháng năm 2018 Sinh viên thực i LỜI CẢM ƠN Việc hoàn thành đồ án tốt nghiệp sau năm học tập trường Đại học Tài nguyên Môi trường TP Hồ Chí Minh đánh dấu móc quan trọng đời học tập, nghiên cứu, thực hành cá nhân em Em xin bày tỏ lòng cảm ơn chân thành biết ơn sâu sắc tới Ths Lê Văn Phùng tận tình hướng dẫn giúp đỡ cho em suốt trình nghiên cứu thực đồ án tốt nghiệp Xin cảm ơn thầy nhiều góp ý bảo tận tình thời gian qua Em xin gửi lời cảm ơn chân thành đến Ban lãnh đạo tập thể cán Đài KTTV tỉnh Bình Định tạo điều kiện thuận lợi cho em có hội thực tập, học hỏi giúp đỡ em hoàn thành đồ án tốt nghiệp Trong trình làm đồ án tốt nghiệp, cố gắng kiến thức kinh nghiệm hạn hẹp nên đồ án khơng tránh khỏi thiếu sót tính tốn phân tích Kính mong q thầy đóng góp ý kiến, giúp em làm tốt hoàn thiện đồ án tốt nghiệp công việc sau trường Cuối cùng, xin dành lời cảm ơn đặc biệt đến q thầy giáo Khoa Khí Tượng Thủy Văn nhiệt tình truyền thụ cho em kiến thức quý giá bổ ích suốt thời gian học tập nghiên cứu trường Em xin chân thành cảm ơn! TP Hồ Chí Minh, ngày tháng năm 2018 Sinh viên thực ii MỤC LỤC DANH MỤC CHỮ VIẾT TẮT .vi DANH MỤC HÌNH .vii DANH MỤC BẢNG viii MỞ ĐẦU 1 TÍNH CẤP THIẾT CỦA ĐỀ TÀI TÌNH HÌNH NGHIÊN CỨU 2.1 Tình hình nghiên cứu dự báo lũ giới 2.2 Tình hình nghiên cứu dự báo lũ nước MỤC ĐÍCH NGHIÊN CỨU PHƯƠNG PHÁP NGHIÊN CỨU 5 DỰ KIẾN KẾT QUẢ NGHIÊN CỨU NỘI DUNG NGHIÊN CỨU, CHƯƠNG TỔNG QUAN VỀ LƯU VỰC SÔNG KÔN – HÀ THANH 1.1 ĐẶC ĐIỂM ĐIỀU KIỆN ĐỊA LÝ TỰ NHIÊN CỦA LƯU VỰC 1.1.1 Vị trí địa lý 1.1.2 Điều kiện địa hình 1.2 ĐẶC ĐIỂM KHÍ TƯỢNG THỦY VĂN 1.2.1 Số nắng 1.2.2 Nhiệt độ 12 1.2.3 Chế độ ẩm 13 1.2.4 Gió 13 1.2.5 Bốc 13 1.2.6 Mưa 15 1.3 SƠNG NGỊI VÀ CHẾ ĐỘ CHẢY 18 1.3.1 Mạng lưới sơng ngịi lưu vực 18 1.3.2 Mạng lưới trạm đo đạc 18 1.4 TÌNH HÌNH KINH TẾ XÃ HỘI TRONG LƯU VỰC 19 iii 1.5 TÌNH TRẠNG THIÊN TAI CỦA KHU VỰC 20 1.5.1 Tình trạng thiên tai lũ lụt 20 1.5.2 Tình trạng sa bồi thủy phá 21 1.6 TÌNH HÌNH SỐ LIỆU 22 1.7 ĐẶC ĐIỂM DÒNG CHẢY NĂM 22 1.8 PHÂN PHỐI DÒNG CHẢY NĂM 22 1.8.1 Phân mùa dòng chảy 23 1.8.2 Phân phối dòng chảy năm 24 1.9 CHẾ ĐỘ DÒNG CHẢY MÙA LŨ 26 1.9.1 Nguyên nhân gây mưa sinh lũ 27 1.9.2 Chế độ dòng chảy lũ 28 1.9.2.1 Mực nước đỉnh lũ 29 1.9.2.2 Lưu lượng đỉnh lũ 30 1.9.3 Các trận lũ điển hình 31 1.10 CHẾ ĐỘ DÒNG CHẢY MÙA CẠN 32 CHƯƠNG 35 TỔNG QUAN VỀ CÁC PHƯƠNG PHÁP DỰ BÁO LŨ VÀ LỰA CHỌN PHƯƠNG PHÁP DỰ BÁO 35 2.1 CÁC PHƯƠNG PHÁP DỰ BÁO LŨ THƯỜNG DÙNG 35 2.1.1 Phương pháp xu 35 2.1.2 Phương pháp lưu lượng mực nước tương ứng 35 2.1.3 Phương pháp lượng trữ 36 2.1.4 Phương pháp phân tích thống kê 37 2.1.5 Phương pháp hồi quy tuyến tính bội 37 2.1.6 Phương pháp mạng Nowrowrron thần kinh - Mạng trí tuệ nhân tạo ANN (Artificial Neural Network) 37 2.2 PHÂN TÍCH LỰA CHỌN PHƯƠNG PHÁP HỒI QUY TUYẾN TÍNH BỘI 38 2.2.1 Cơ sở phương pháp 38 2.2.2 Các tài liệu kế thừa 38 2.2.3 Phân tích lựa chọn 38 iv CHƯƠNG 40 XÂY DỰNG PHƯƠNG ÁN DỰ BÁO MỰC NƯỚC ĐỈNH LŨ TẠI TRẠM THỦY VĂN DIÊU TRÌ TRÊN SƠNG HÀ THANH – TỈNH BÌNH ĐỊNH BẰNG PHƯƠNG PHÁP HỒI QUY TUYẾN TÍNH BỘI 40 3.1 TÌNH HÌNH CƠNG TÁC DỰ BÁO THỦY VĂN TRÊN LƯU VỰC TRƯỚC ĐÂY 40 3.2 TÌNH HÌNH CƠNG TÁC DỰ BÁO THỦY VĂN TRÊN LƯU VỰC HIỆN NAY 41 3.3 XÂY DỰNG PHƯƠNG ÁN DỰ BÁO MỰC NƯỚC ĐỈNH LŨ TẠI TRẠM THỦY VĂN DIÊU TRÌ TRÊN SƠNG HÀ THANH – TỈNH BÌNH ĐỊNH BẰNG PHƯƠNG PHÁP HỒI QUY TUYẾN TÍNH BỘI 41 3.3.1 Cơ sở lí luận phương pháp 41 3.3.2 Dự báo thử nghiệm đánh giá kết 50 Thuyết minh kết phương án 50 Xác định thời gian chảy truyền 51 Cơ sở phương pháp 51 Các bước tính tốn xây dựng phương án dự báo 52 KẾT LUẬN – KIẾN NGHỊ 54 PHỤ LỤC TÍNH TỐN 55 TÀI LIỆU THAM KHẢO 64 v DANH MỤC CHỮ VIẾT TẮT Số thứ tự Danh mục chữ viết tắt Ý nghĩa B Hệ số phân phối dòng chảy H Mực nước Q Lưu lượng X Lượng mưa KTTV Khí tượng thủy văn KCN Khu công nghiệp ODA Dự án hỗ trợ phát triển thức BCH PCLB Ban huy phòng chống lụt bão TT Thứ tự 10 UBND Ủy ban nhân dân vi DANH MỤC HÌNH Hình 1.1 Bản đồ lưu vực sơng Kơn – Hà Thanh Hình 1.2 Mơ hình phân phối số nắng năm trung bình nhiều năm Trạm Khí tượng Quy Nhơn 11 Hình 1.3 Mơ hình phân phối số nắng năm trung bình nhiều năm Trạm Khí tượng An Nhơn 11 Hình1 Mơ hình phân phối số nắng năm trung bình nhiều năm Trạm Khí tượng Hồi Nhơn 12 Hình 1.5 Mơ hình phân phối bốc năm trung bình nhiều năm Trạm khí tượng Quy Nhơn 14 Hình1 Mơ hình phân phối bốc năm trung bình nhiều năm Trạm Khí tượng Hồi Nhơn 15 Hình 1.7 Mơ hình phân phối mưa năm trung bình nhiều năm Trạm Khí tượng Quy Nhơn 17 Hình1 Mơ hình phân phối mưa năm trung bình nhiều năm Trạm Khí tượng Hồi Nhơn 17 Hình 1.9 Phân phối dịng chảy năm trạm thủy văn Bình Tường 25 Hình 3.1 Quan hệ mực nước Vân Canh thời gian chảy truyền 51 vii DANH MỤC BẢNG Bảng 1.1 Tổng số nắng trung bình tháng năm 10 Bảng 1.2 Số ngày nắng trung bình tháng năm 10 Bảng 1.3 Nhiệt độ tháng Bình Định 12 Bảng 1.4 Độ ẩm tương đối trung bình tháng năm 13 Bảng 1.5 Tốc độ gió trung bình nhiều năm tháng năm 13 Bảng 1.6 Lượng bốc tháng Bình Định 14 Bảng 1.7 Lượng mưa tháng, năm trung bình nhiều năm trạm Quy Nhơn 16 Bảng 1.8 Lượng mưa tháng, năm trung bình nhiều năm trạm Hồi Nhơn 16 Bảng 1.9 Mạng lưới trạm yếu tố quan trắc 19 Bảng 1.10 Số liệu đặc trưng trận lũ ngày 18- 19/XI/1987 20 Bảng 1.11 Số liệu đặc trưng trận lũ ngày 01- 07/XII/1999 21 Bảng 1.12 Các đặc trưng dịng chảy lưu vực sơng Kơn – Hà Thanh 22 Bảng 1.13 Tần suất xuất dịng chảy bình qn tháng ≥ 8,3% dòng chảy năm (%) 24 Bảng 1.14 Phân phối dòng chảy theo mùa 24 Bảng 1.15 Phân phối dòng chảy theo tháng năm 25 Bảng 1.16 Trung bình số trận lũ xuất năm 28 Bảng 1.17 Tỷ lệ % lũ xuất lớn năm 29 Bảng 1.18 Mực nước cao trạm lưu vực sông Kôn – Hà Thanh 29 Bảng 1.19 Đặc trưng cường suất, biên độ lũ trạm 30 Bảng 1.20 Lưu lượng lớn trạm 30 Bảng 1.21 Tần suất lưu lượng đỉnh lũ, môđun đỉnh lũ 30 Bảng 1.22 Số liệu đặc trưng trận lũ ngày 18- 19/XI/1987 31 Bảng 1.23 Số liệu đặc trưng trận lũ ngày 01- 07/XII/1999 32 Bảng 1.24 Khả xảy kiệt ngày nhỏ năm trạm Bình Tường (%) 33 Bảng 1.25 Dòng chảy nhỏ trạm lân cận lưu vực 34 viii Biểu diễn hệ phương trình theo cách nhân ma trận ta được: n  x1 n i =1 n  x1 i =1 n  ( x1.x1) i =1 n  x2 …… i =1 n  ( x1.x2) …… i =1 n  xm i =1  xm i =1 n  ( xm.x1) i =1 n  ( xm.x2) …… i =1 n  ( x1.xm)  ( y.x1) i =1 n  ( xm.xm) i =1 n y i =1 n …………………………………………………  n a a i =1 = … …… n  ( y.xm) i =1 A  B=C Hệ phương trình có nghiệm detA  0, tức ma trân A không suy biến Sử dụng biến đổi ma trận tương đương để biến đổi ma trân A ma trận tam giác trên, ma trận B C biến đổi tương đương theo A Giải ngược hệ theo phương pháp để tìm nghiệm hệ.(Nguồn: PGS.TS Nguyễn Thông – Kinh tế lượng Hà Nội – 2000) 3.3.2 Dự báo thử nghiệm đánh giá kết Dùng phương pháp hồi quy với nhân tố ảnh hưởng mưa lưu vực để xây dựng phương án dự báo mực nước đỉnh lũ cho trạm Diêu Trì Thuyết minh kết phương án Vùng đồng hạ lưu sơng Hà Thanh vùng trọng điểm phía nam tỉnh Bình Định hoạt động kinh tế, trị, xã hội Vào mùa bão lũ vùng hạ lưu sông Hà Thanh bị thiệt hại nặng nề, nhiều xã ngập hoàn toàn nước, nhiều đoạn đường giao thông liên xã, liên huyện bị ách tắc Dự báo lũ cho sơng Hà Thanh góp phần khơng nhỏ việc giảm nhẹ thiệt hại thiên tai bão lũ gây mục tiêu quan trọng giúp cho việc chủ động điều hành hoạt động kinh tế, trị, xã hội cho tỉnh nhà 50 Để lập phương án dự báo chọn tài liệu từ năm 1996 – 2014 trạm thủy văn Vân Canh Diêu Trì, tài liệu chọn sau (Xem phụ lục bảng tỉnh toán) Xác định thời gian chảy truyền Áp dụng phương pháp điểm đặc trưng tương ứng để xác định  Vì tài liệu trạm đọc thưa nên phần chân lũ khó xác định thời gian xuất mực nước chân lũ điểm chuyển tiếp xác, xác định  theo đỉnh lũ tương ứng hai trạm Thời gian truyền lũ  chênh lệch thời gian xuất điểm đặc trưng tương ứng đường trình hai trạm, kết tính lập thành bảng (3.1) Qua bảng (3.1) hình 3.1 thấy khó xác định rõ quan hệ  = f(Hvc), lấy giá trị trung bình  cho cấp mực nước, kết quả: tb = /n = 372/63 = 5.90 = 6h Quan hệ mực nước Vân Canh thời gian chảy truyền 4900 4800 Hvc (cm) 4700 4600 4500 4400 4300 4200 4100 10 12 14 Thời gian chảy truyền  (h) Hình 3.1 Quan hệ mực nước Vân Canh thời gian chảy truyền Cơ sở phương pháp Hồi quy tuyến tính bội dựa số liệu quan trắc khứ để thiết lập tuyến tính mơ tả mối quan hệ yếu tố dự báo nhân tố ảnh hưởng Trong trường hợp tổng quát biểu diễn dạng sau: 51 Ht +  = C + ∑ m i=1 iXi,t Trong đó: Ht + : yếu tố cần dự báo Xi,t: nhân tố ảnh hưởng m: số lượng nhân tố ảnh hưởng i: hệ số hồi quy C: số t, : thời điểm phát báo thời điểm dự kiến Các bước tính tốn xây dựng phương án dự báo Xuất phát từ đặc điểm lưu vực sơng Hà Thanh tình hình điện báo có, đề xuất phương án dự báo mực nước lũ Diêu Trì dạng quan hệ tương quan lấy quan hệ mưa, mực nước tương ứng Vân Canh Diêu Trì vì: Trạm Vân Canh nằm thượng nguồn sơng Hà Thanh cách tâm lưu vực khống chế khoảng 10km Khu vực có độ dốc lịng sơng lớn nên thời gian tập trung nước đến trạm Vân Canh khoảng HDiêu Trì (t+) = f(HVân Canh (t), HDiêu Trì (t), XVân Canh (t)) Trong đó: : thời gian truyền lũ từ Vân Canh đến Diêu Trì HDiêu Trì (t+): mực nước trạm Diêu Trì thời điểm t+ HVân Canh (t): mực nước trạm Vân Canh thời điểm t HDiêu Trì (t): mực nước trạm Diêu Trì thời điểm t XVân Canh (t): lượng mưa thời đoạn trạm Vân Canh đến thời điểm t * Xây dựng phương án dự báo: Dùng tài liệu mưa trạm Vân Canh với thời đoạn 6, 12 tính đến thời điểm tính tốn, mực nước đỉnh trạm Vân Canh mực nước tương ứng trạm Diêu Trì thời gian chảy truyền  = 6h (đã xác định) lập thành bảng để xây dựng đường hồi quy dự báo Kết thu phương trình tương quan hồi quy dạng sau: HDiêu Trì (t+) = 0.161 HVân Canh (t) + 0.513 HDiêu Trì (t) + 0.619 XVân Canh (6h) – 495 52 Trong đó: : thời gian truyền lũ từ Vân Canh đến Diêu Trì lấy HDiêu Trì (t+): mực nước dự báo trạm Diêu Trì HVân Canh (t): mực nước trạm Vân Canh thời điểm t HDiêu Trì (t): mực nước trạm Diêu Trì thời điểm t XVân Canh (6h): lượng mưa trạm Vân Canh thời đoạn đến thời điểm t * Đánh giá sai số phương án + Sai số cho phép: Scp = 65cm, với: ∑(Hd−Hbq)2 =√ n−1 ∑(Hdb−Hd)2 ‘= √ 𝑛−3 = 97cm = 42cm ‘/  = 42/97 = 0.43 σ′  = √1 − ( )2 = 0.90 σ Hd: mực nước thực đo Hbq: mực nước bình quân thực đo Hbd: mực nước dự báo tính theo quan hệ + Mức đảm bảo phương án dự báo P% = m/n = 55/63 = 87% m: tổng số điểm có ΔH < Scp n: tổng số điểm lập quan hệ Qua kết thấy phương án dự báo mực nước đỉnh lũ cho trạm Diêu Trì dùng tốt Xem kết đánh giá sai số phương án bảng (3.4) 53 KẾT LUẬN – KIẾN NGHỊ Kết luận Tính tốn phân tích lập phương án dự báo mực nước đỉnh lũ cho trạm Diêu Trì thực theo bước sau: - Thu thập số liệu mực nước đỉnh lũ trạm Vân Canh Diêu Trì để tính thời gian chảy truyền theo đỉnh lũ tương ứng hai trạm, lấy trung bình - Dựa vào số liệu lượng mưa Vân Canh lúc 12 mực nước Hmax Vân Canh Diêu Trì để tính H (t), H (t+6) trạm Diêu Trì để xây dựng quan hệ dự báo - Mức đảm bảo phương án dự báo mực nước đỉnh lũ trạm Diêu Trì phương pháp hồi quy tuyến tính bội đạt 87% Dự báo đỉnh lũ trạm Diêu Trì phương pháp đạt kết tốt Kiến nghị Để xây dựng phương án dự báo đỉnh lũ cho trạm Diêu Trì phương pháp hồi quy tuyến tính bội đạt kết xác cao cần phải tối ưu hóa yếu tố sau - Khi tính tốn cần xử lý sai số q trình tính toán cho kết dự báo tốt - Ghi chép, nhập số, liệu tính tốn phải chỉnh xác Việc đo đạc, thu thập số liệu phải tuân theo quy định, quy phạm ngành Phương pháp dự báo mực nước đỉnh lũ cho trạm thủy văn Diêu Trì theo phương pháp hồi quy tuyến tính bội đạt kết tốt, sử dụng phương pháp để dự báo cho trạm thủy văn khác vùng hạ lưu sông nhằm phục vụ cho quy hoạch, kinh tế Nâng cao chất lượng dự báo nhằm hạn chế lũ gây thiệt hại nặng sống kinh tế người dân nơi 54 PHỤ LỤC TÍNH TỐN Bảng 3.1 Thời gian chảy truyền  Vân Canh – Diêu Trì  Năm 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 1996 1997 1998 1999 2000 2001 2002 2003 Vân Canh H(cm) Tgian 4207 07h 17/X 4305 01h 28/X 4377 04h 03/XI 4478 10h 16/XI 4479 01h 18/XI 4482 12h 01/XII 4254 20h 30/X 4449 09h 02/XI 4390 08h 04/XI 4421 13h 21/X 4517 19h 22/X 4370 07h 07/XI 4428 15h 13/XI 4534 05h 20/XI 4487 06h 26/XI 4375 21h 10/XII 4424 10h 14/XII 4280 06h 18/X 4353 16h 23/X 4421 10h 24/X 4359 16h 02/XI 4442 11h 05/XI 4623 20h 21/XI 4447 11h 02/XII 4420 07h 03/XII 4387 19h 05/XII 4377 07h 16/XII 4180 07h 10/X 4244 07h 18/X 4303 13h 13/XI 4292 13h 17/XI 4170 03h 22/X 4246 05h 12/XI 4158 07h 23/IX 4287 13h 11/XI 4157 16h 04/X H(cm) 308 307 470 454 540 497 287 487 434 440 550 429 455 554 571 460 500 271 356 386 383 462 527 436 462 396 453 293 418 432 397 277 344 341 481 349 Diêu Trì Tgian 11h 17/X 07h 28/X 07h 03/XI 18h 16/XI 07h 18/XI 18h 01/XII 02h 31/X 16h 02/XI 13h 04/XI 20h 21/X 23h 22/X 11h 07/XI 23h 13/XI 11h 20/XI 09h 26/XI 01h 11/XII 15h 14/XII 17h 18/X 03h 24/X 18h 24/X 01h 03/XI 19h 05/XI 01h 22/XI 17h 02/XII 11h 03/XII 01h 06/XII 13h 16/XII 13h 10/X 13h 18/X 21h 13/XI 22h 17/XI 09h 22/X 12h 12/XI 11h 23/IX 18h 11/XI 19h 04/X  (giờ) 6 7 4 11 11 6 6 55 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 Tổng 2005 2007 2008 2009 2010 2011 2013 2014 4319 4364 4296 4248 4292 4314 4161 4467 4272 4233 4360 4368 4241 4142 4254 4805 4329 4296 4219 4277 4336 4221 4190 4434 4243 4342 4176 372 07h 19/X 07h 20/X 13h 25/X 07h 15/XII 12h 17/X 3h 27/X 7h 30/X 18h 04/XI 16h 23/XI 07h 17/X 10h 19/X 02h 25/XI 07h 31/XII 01h 29/IX 16h 17/X 22h 02/XI 22h 01/XI 15h 08/XI 08h 19/XI 14h 29/XI 07h 20/X 07h 07/XI 10h 26/XI 10h 15/XI 15h 17/XI 04h 30/XI 09h 13/XII 510 386 512 527 462 547 310 608 488 320 520 526 372 180 414 733 450 531 455 537 412 452 438 656 400 430 453 12h 19/X 14h 20/X 19h 25/X 07h 15/XII 21h 17/X 6h 27/X 11h 30/X 23h 04/XI 19h 23/XI 10h 17/X 14h 19/X 14h 25/XI 1/1/2009 07h 29/IX 22h 17/X 01h 03/XI 9h 02/XI 19h 08/XI 13h 19/XI 15h 29/XI 11h 20/X 13h 07/XI 13h 26/XI 15h 15/XI 19h 17/XI 09h 30/XI 10h 13/XII 3 12 24 6 11 5 56 Bảng 3.2 Số liệu xây dựng quan hệ dự báo TT Năm 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 1996 1997 1998 1999 2000 2001 2002 2003 Tháng Ngày X X XI XI XI XII X XI XI X X XI XI XI XI XII XII X X X XI XI XI XII XII XII XII X X XI XI X XI IX XI X X 17 28 16 18 30 21 22 13 20 26 10 14 18 23 24 21 16 10 18 13 17 22 12 23 11 19 Giờ X Vân Canh (6h) (mm) X Vân Canh (12h) (mm) 10 12 20 13 19 15 21 10 16 10 16 11 20 11 19 7 13 13 13 16 12 32 59 68 90 70 18 123 60 72 130 70 41 67 75 44 99 32 52 35 30 50 100 27 63 28 40 13 116 113 62 19 52 67 86 46 52 87 37 106 91 95 118 33 129 75 133 178 90 80 100 148 86 125 81 101 54 47 63 125 54 87 30 77 63 116 169 103 51 117 99 120 81 84 Hmax Vân Canh (cm) 4207 4305 4377 4478 4479 4482 4254 4449 4390 4421 4517 4370 4428 4534 4487 4375 4424 4280 4353 4421 4359 4442 4623 4447 4420 4387 4377 4180 4244 4303 4292 4170 4246 4158 4287 4157 4319 Hmax Diêu Trì (cm) 308 307 470 454 540 497 287 487 434 440 550 429 455 554 571 460 500 271 356 386 383 462 527 436 462 396 453 293 418 432 397 277 344 341 481 349 510 H () Diêu Trì (cm) 264 296 441 360 482 436 278 308 397 335 470 383 360 473 522 449 418 192 239 330 332 334 280 401 447 346 360 267 267 267 347 240 230 250 438 311 496 H (t+6) Diêu Trì (cm) 308 307 428 431 540 497 287 484 428 438 456 423 450 554 555 429 499 244 327 369 369 438 526 436 456 396 453 293 418 427 395 277 340 338 477 345 510 57 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 2005 2007 2008 2009 2010 2011 2013 2014 XI X X X X X XI XI X XI XI XII IX X XI XI XI XI XI X XI XI XI XI XI XII 13 20 25 17 27 30 23 17 19 25 31 29 17 19 29 20 26 15 17 30 13 13 12 18 16 10 12 16 22 22 15 14 7 10 10 15 62 111 84 111 87 54 131 18 80 96 70 44 65 356 40 58 31 67 63 80 36 200 48 42 18 120 118 124 142 135 57 204 53 88 151 75 85 45 120 474 63 100 55 133 71 111 63 258 56 106 90 4326 4364 4296 4292 4314 4161 4467 4272 4233 4360 4368 4241 4142 4254 4805 4329 4296 4219 4277 4336 4221 4190 4434 4243 4342 4176 539 386 512 462 547 310 608 488 320 520 526 372 180 414 733 450 531 455 537 412 452 438 656 400 430 453 526 166 471 159 496 240 543 467 200 440 453 330 175 361 686 330 490 389 532 283 400 390 456 390 253 433 520 375 512 462 547 310 608 456 310 511 501 335 180 410 733 435 531 455 512 402 452 403 646 380 429 427 58 Bảng 3.3 Bảng tính hệ số phương trình tương quan tuyến tính bội STT Dự báo Diêu Trì 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 308 307 428 431 540 497 287 484 428 438 456 423 450 554 555 429 499 244 327 369 369 438 526 436 456 396 453 293 418 427 395 277 340 338 477 345 510 520 X1 Vân Canh 4207 4305 4377 4478 4479 4482 4254 4449 4390 4421 4517 4370 4428 4534 4487 4375 4424 4280 4353 4421 4359 4442 4623 4447 4420 4387 4377 4180 4244 4303 4292 4170 4246 4158 4287 4157 4319 4326 X2 Diêu Trì 264 296 441 360 482 436 278 308 397 335 470 383 360 473 522 449 418 192 239 330 332 334 280 401 447 346 360 267 267 267 347 240 230 250 438 311 496 526 X3 Mưa (mm) 12 32 59 68 90 70 18 123 60 72 130 70 41 67 75 44 99 32 52 35 30 50 100 27 63 28 40 13 116 113 62 19 52 67 86 46 52 62 59 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 375 512 462 547 310 608 456 310 511 501 335 180 410 733 435 531 455 512 402 452 403 646 380 429 427 4364 4296 4292 4314 4161 4467 4272 4233 4360 4368 4241 4142 4254 4805 4329 4296 4219 4277 4336 4221 4190 4434 4243 4342 4176 166 471 159 496 240 543 467 200 440 453 330 175 361 686 330 490 389 532 283 400 390 456 390 253 433 111 84 111 87 54 131 18 80 96 70 44 65 356 40 58 31 67 63 80 36 200 48 42 18 60 SUMMARY OUTPUT Regression Statistics Multiple R 0.905951 R Square 0.820748 Adjusted R Square 0.811633 Standard Error 42.15277 Observations 63 ANOVA Regression Residual Total df 59 62 SS 480008.8 104834.5 584843.3 Intercept X Variable X Variable X Variable Standard Coefficients Error -495.426 244.2876 0.160599 0.058838 0.512817 0.057221 0.618947 0.134958 Significance MS F F 160002.9 90.04833 5.35E-22 1776.856 t Stat -2.02804 2.729524 8.962055 4.58622 P-value 0.047079 0.008347 1.32E-12 2.4E-05 Lower 95% -984.244 0.042865 0.398318 0.348897 Upper 95% -6.60778 0.278333 0.627316 0.888997 Lower 95.0% -984.244 0.042865 0.398318 0.348897 Upper 95.0% -6.60778 0.278333 0.627316 0.888997 Các số màu đỏ hệ số cần tìm H(V Canh) = H (Dieu Tri) = Xbqlv = DB H(Dieu Tri) = 4300 300 100 411 tren bao dong I 61 Bảng 3.4 Bảng tính sai số phương án dự báo Năm Tgian XVCanh (6h) Hmax Vân Canh (cm) H (t) Diêu Trì (cm) H (t+6) Diêu Trì (cm) HDiêu Trì Dự báo (cm) ΔH1= Hđo Htb (ΔH1)2 ΔH2= Hđo Htb (ΔH2)2 1996 07h 17/X 12 4207 264 308 325 -124 15376 -17 295 01h 28/X 32 4305 296 307 370 -125 15625 -63 3939 04h 03/XI 59 4377 441 428 472 -4 16 -44 1976 10h 16/XI 68 4478 360 431 453 -1 -22 472 01h 18/XI 90 4479 482 540 529 108 11664 11 119 12h 01/XII 70 4482 436 497 494 65 4225 12 20h 30/X 18 4254 278 287 344 -145 21025 -57 3209 09h 02/XI 123 4449 308 484 455 52 2704 29 816 08h 04/XI 60 4390 397 428 453 -4 16 -25 605 13h 21/X 72 4421 335 438 433 36 23 19h 22/X 130 4517 470 456 554 24 576 -98 9568 07h 07/XI 70 4370 383 423 448 -9 81 -25 644 15h 13/XI 41 4428 360 450 428 18 324 22 485 05h 20/XI 67 4534 473 554 519 122 144884 35 1218 06h 26/XI 75 4487 522 555 542 123 15129 13 179 21h 10/XII 44 4375 449 429 467 -3 -38 1440 10h 14/XII 99 4424 418 499 493 67 4489 36 06h 18/X 32 4280 192 244 312 -188 35344 -68 4676 16h 23/X 52 4353 239 327 361 -105 11025 -34 1131 1997 1998 1999 2000 2001 2002 2003 10h 24/X 35 4421 330 369 408 -63 3969 -39 1500 16h 02/XI 30 4359 332 369 396 -63 3969 -27 712 11h 05/XI 50 4442 334 438 422 36 16 242 20h 21/XI 100 4623 280 526 455 94 8836 71 5063 11h 02/XII 27 4447 401 436 443 16 -7 55 07h 03/XII 63 4420 447 456 485 24 576 -29 837 19h 05/XII 28 4387 346 396 406 -36 1296 -10 103 07h 16/XII 40 4377 360 453 419 21 441 34 1147 07h 10/X 13 4180 267 293 323 -139 19321 -30 900 07h 18/X 116 4244 267 418 397 -14 196 31 439 13h 13/XI 113 4303 267 427 405 -5 25 22 497 13h 17/XI 62 4292 347 395 412 -37 1369 -17 303 03h 22/X 19 4170 240 277 311 -155 24025 -34 1173 05h 12/XI 52 4246 230 340 339 -92 8464 1 07h 23/IX 67 4158 250 338 344 -94 8836 -6 38 13h 11/XI 86 4287 438 477 473 45 2025 15 16h 04/X 46 4157 311 345 362 -87 7569 -17 299 07h 19/X 52 4319 496 510 487 78 6084 23 529 03h 13/XI 62 4326 526 520 510 88 7744 10 106 Đánh giá không đạt không đạt không đạt 62 2005 2007 2008 2009 2010 2011 2013 2014 07h 20/X 111 4364 166 375 361 -57 3249 14 183 13h 25/X 84 4296 471 512 490 80 6400 22 472 12h 17/X 111 4292 159 462 346 30 900 116 13389 3h 27/X 87 4314 496 547 508 115 13225 39 1532 7h 30/X 54 4161 240 310 331 -122 14884 -21 461 18h 04/XI 131 4467 543 608 584 176 30976 24 584 16h 23/XI 18 4272 467 456 444 24 576 12 156 07h 17/X 80 4233 200 310 339 -122 14884 -29 820 10h 19/X 96 4360 440 511 492 79 6241 19 357 02h 25/XI 70 4368 453 501 484 69 4761 17 290 07h 31/XII 44 4241 330 335 384 -97 9409 -49 2433 01h 29/IX 4142 175 180 267 -252 63504 -87 7605 16h 17/X 65 4254 361 410 415 -22 484 -5 28 22h 02/XI 356 4805 686 733 851 301 90601 -118 13897 22h 01/XI 40 4329 330 435 396 39 1520 15h 08/XI 58 4296 490 531 484 99 9801 47 2216 08h 19/XI 31 4219 389 455 403 23 529 52 2703 14h 29/XI 67 4277 532 512 508 80 6400 16 07h 20/X 63 4336 283 402 387 -30 900 15 217 07h 07/XI 80 4221 400 452 439 20 400 13 161 10h 26/XI 36 4190 390 403 402 -29 841 1 10h 15/XI 200 4434 456 646 577 214 45796 69 4816 15h 17/XI 48 4243 390 380 418 -52 2704 -38 1437 04h 30/XI 42 4342 253 429 360 -3 69 4782 09h 13/XII 18 4176 433 427 411 -5 25 16 Tổng cộng 27190 Trung Bình 432 584854 khơng đạt không đạt không đạt không đạt không đạt 269 105150 63 TÀI LIỆU THAM KHẢO PGS.TS Đặng Văn Bảng - Giáo trình Dự báo thủy văn Nhà xuất nông nghiệp (Hà Nội – 2000) PGS.TS Lê Văn Nghinh – Giáo trình nguyên lý thủy văn Nhà xuất Nông Nghiệp (Hà Nội – 2000) ThS Nguyễn Tấn Hương cộng tác viên: KS Thiệu Quang Tân, KS Trần Sỹ Dũng, KS Thân Văn Đón, KS Lê Văn Dũng, KS Nguyễn Ngọc Quỳnh, KS Nguyễn Trung Thiếp, KS Võ Anh Kiệt, KS Nguyễn Văn Lý, KS Lương Ngọc Lũy – Báo cáo kết nghiên cứu khoa học đề tài Đặc Điểm Khí Hậu Thủy Văn tỉnh Bình Định (2004) Phạm Việt Hùng - Tai biến thiên nhiên tỉnh Bình Định giải pháp thích ứng, giảm nhẹ Trung tâm nghiên cứu tài nguyên môi trường Đại học quốc gia Hà Nội Quy hoạch phát triển kinh tế - xã hội tỉnh Bình Định đến năm 2020 Tác giả: Nguyễn Văn Tuần, Đoàn Quyết Trung, Bùi Văn Đức – Giáo trình Dự báo thủy văn Nhà xuất Đại học quốc gia Hà Nội (Hà Nội – 2003) PGS.TS Nguyễn Thông - Kinh tế lượng (Hà Nội – 2000) 64 ... PHƯƠNG ÁN DỰ BÁO MỰC NƯỚC ĐỈNH LŨ TẠI TRẠM THỦY VĂN DIÊU TRÌ TRÊN SƠNG HÀ THANH – TỈNH BÌNH ĐỊNH BẰNG PHƯƠNG PHÁP HỒI QUY TUYẾN TÍNH BỘI 40 3.1 TÌNH HÌNH CƠNG TÁC DỰ BÁO THỦY VĂN TRÊN... TÁC DỰ BÁO THỦY VĂN TRÊN LƯU VỰC HIỆN NAY 41 3.3 XÂY DỰNG PHƯƠNG ÁN DỰ BÁO MỰC NƯỚC ĐỈNH LŨ TẠI TRẠM THỦY VĂN DIÊU TRÌ TRÊN SƠNG HÀ THANH – TỈNH BÌNH ĐỊNH BẰNG PHƯƠNG PHÁP HỒI QUY. .. sông Kôn – Hà Thanh Chương : Tổng quan phương pháp dự báo lũ lựa chọn phương pháp dự báo Chương : Xây dựng phương án dự báo mực nước đỉnh lũ Trạm thủy văn Diêu Trì sơng Hà Thanh – Tỉnh Bình Định