NGHIÊN CỨU & TRAO ĐỔI NGHIÊN CỨU ỨNG DỤNG MÔ HÌNH TỐN PHỤC VỤ DỰ BÁO LŨ LỚN TRÊN LƯU VỰC SƠNG LẠI GIANG TỈNH BÌNH ĐỊNH Trịnh Xn Mạnh, Lê Thị Thường Đại học Tài nguyên Môi trường Hà Nội iện lũ lụt nguyên nhân gây nhiều tác động tiêu cực đến đời sống người dân khu vực Miền Trung nói chung lưu vực sơng Lại Giang, tỉnh Bình Định nói riêng Bài báo tập trung vào việc nghiên cứu thiết lập mơ hình tốn gồm mơ hình mưa dịng chảy MIKE NAM, thuỷ lực MIKE 11 xây dựng phương án dự báo với thời gian dự kiến 24h lưu vực sông Kết nghiên cứu cho thấy mơ hình tốn áp dụng cho kết khả quan, với độ tin cậy mơ hình 75% Ngồi dự báo thử nghiệm cho hai trận lũ lớn năm 1990 1992 cho mức đảm bảo phương án 80%, đồng thời phương án dự báo đánh giá đạt Như vậy, số lượng trạm khí tượng thuỷ văn cịn nguồn số liệu cịn thiếu song thấy mơ hình tốn có khả áp dụng tốt vào công tác dự báo thực tế lưu vực sơng Lại Giang Từ khố: Lũ lớn, Lại Giang, Dự báo lũ, MIKE NAM, MIKE 11 H Mở đầu Lũ lụt thiên tai liên quan đến dòng chảy gây nên thiệt hại vô to lớn người tài sản, đồng thời gây tác động xấu tới môi trường tự nhiên Đặc biệt, tỉnh miền Trung Việt Nam lại chịu ảnh hưởng mạnh mẽ rõ rệt tác động cộng hưởng điều kiện địa hình, lịng dẫn sơng đặc điểm thời tiết khí hậu Bình Định tỉnh thành nằm duyên hải Nam Trung Bộ nên hàng năm chịu ảnh hưởng nặng nề lũ lụt Trên lưu vực sơng Lại Giang nói riêng tỉnh Bình Định nói chung, lũ lụt gây nhiều thiệt hại người tài sản Điển năm 1998, tỉnh Bình Định có 62 người chết, 10 người bị thương, hư hại nặng 800 nhà người dân, sập đổ 14 trường học phòng y tế; tổng diện tích lúa bị ngập nước 23,378 ha, ước tính thiệt hại lũ gây năm khoảng 150 tỉ đồng [4] Trước tình hình lũ lụt diễn biến phức tạp, cơng tác cảnh báo, dự báo lũ đặc 16 TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 10 - 2016 biệt dự báo lũ lớn đóng vai trị vơ quan trọng Tuy nhiên số lượng trạm quan trắc, đo đạc khí tượng thuỷ văn lưu vực lại cịn q ít, chưa đảm bảo mật độ Bên cạnh đó, thời gian quan trắc, đo đạc cịn ngắn, dẫn đến nguồn số liệu đo đạc thiếu chưa đồng Đây coi thách thức cơng tác dự báo lũ Như biết, phòng tránh lũ lụt vấn đề ưu tiên hàng đầu nhằm hạn chế lũ lụt thiệt hại lũ lụt gây Trong quan trọng vấn đề cảnh báo, dự báo lũ sớm nhằm tránh tổn thất to lớn lũ gây nên Trải qua nhiều thời kỳ phát triển, với lớn mạnh không ngừng khoa học công nghệ thông tin nên công tác có nhiều phát triển Trong thủy văn học có nhiều phương pháp dự báo từ truyền thống đến đại Trong nghiên cứu tập trung vào ứng dụng phương pháp mơ hình tốn dự báo lũ cho lưu vực sơng Lại Giang Trong có sử dụng đến hai mơ hình tốn gồm mơ hình mưa dịng chảy MIKE NAM NGHIÊN CỨU & TRAO ĐỔI có nhiệm vụ tính tốn đầu vào cho mơ hình thuỷ lực MIKE 11 HD để diễn tốn dịng chảy xuống hạ du, đồng thời hai mơ hình kết nối với tạo thành mơ hình số hồn chỉnh phục vụ công tác dự báo lũ Mục tiêu nghiên cứu thiết lập hai mơ hình tốn với độ tin cậy mơ hình cao thơng qua q trình hiệu chỉnh kiểm định mơ hình Bên cạnh đó, nghiên cứu tập trung vào xây dựng phương án dự báo với thời gian dự kiến 24h tiến hành dự báo thử nghiệm để đánh giá chất lượng phương án dự báo đề Tổng quan lưu vực phương pháp nghiên cứu 2.1 Tổng quan lưu vực nghiên cứu Hình Bản đồ mạng lưới sơng trạm khí tượng thuỷ văn lưu vực Lại Giang Sông Lại Giang có vị trí địa lý nằm khoảng 14030’ tới 14040’ vĩ độ bắc 108050’ tới 108060’ kinh độ đơng Lưu vực nghiên cứu thuộc tỉnh Bình Định - tỉnh ven biển miền Trung Việt Nam Lưu vực Lại Giang nằm phía đơng dãy Trường Sơn trải dài đến Biển Đơng, tỉnh Bình Định Phía Đơng giáp với Biển Đơng lưu vực sơng Xương Phía Tây giáp với huyện Vĩnh Thạnh lưu vực sông Kơn Phía Nam giáp với huyện Phù Mỹ sơng Châu Trúc cuối phía Bắc giáp với tỉnh Quảng Ngãi [1] Lại Giang sông lớn thứ hai tỉnh Bình Định, sơng hình thành từ hợp hai dịng sơng An Lão Kim Sơn Lũ dạng lũ điển hình lưu vực sơng loại trung bình thuộc vùng đồi núi miền Trung nước ta Theo đó, lũ lên nhanh sau mưa lớn, có cường suất lớn thời gian kéo dài không lâu thường từ đến ngày tuần [3] 2.2 Phương pháp mơ hình tốn 2.2.1 Mơ hình mưa – dịng chảy MIKE NAM Mơ hình NAM xây dựng Khoa Thuỷ văn Viện Kỹ thuật Thuỷ động lực Thuỷ lực thuộc Đại học Kỹ thuật Đan Mạch năm 1982 NAM chữ viết tắt cụm từ tiếng Đan Mạch “Nedbør - Afstrømnings - Models” có nghĩa mơ hình mưa rào dịng chảy Mơ hình NAM với ngun lý dùng hàm phi tuyến để mơ tả dịng chảy vào sơng, mơ hình thường sử dụng để dự báo lũ mưa rào tuyết tan Trong mô hình NAM, lưu vực xem đơn vị xử lý Do đó, thơng số biến đại diện cho giá trị trung bình hóa tồn lưu vực Mơ hình tính q trình mưa – dịng chảy theo cách tính liên tục hàm lượng ẩm bể chứa riêng biệt có tương tác lẫn Cấu trúc mơ hình NAM xây dựng nguyên tắc hồ chứa theo chiều thẳng đứng hồ chứa tuyến tính, gồm bể chứa theo chiều thẳng đứng gồm bể chứa tuyết tan, bể chứa mặt, bể chứa tầng dưới, bể chứa ngầm tầng bể chứa ngầm tầng Trong mô hình NAM dịng chảy mặt xác định lượng trữ bề mặt tràn, U > Umax, lượng nước thừa PN gia nhập vào thành phần dòng chảy mặt Thông số QOF đặc trưng cho phần nước thừa PN đóng góp vào dịng chảy mặt Nó giả thiết tương ứng với PN biến đổi tuyến tính theo quan hệ lượng trữ ẩm đất, L/Lmax, lượng trữ ẩm tầng thấp QOF L / Lmax TOF PN nÕu L / Lmax ! TOF °CQOF TOF ® °0 nÕu L / Lmax d TOF ¯ TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 10 - 2016 17 NGHIÊN CỨU & TRAO ĐỔI Trong CQOF= hệ số dịng chảy tràn mặt đất (0 ≤ CQOF ≤ 1), TOF = giá trị ngưỡng dòng chảy tràn (0 ≤ TOF ≤ 1) 2.2.2 Mơ hình thuỷ lực MIKE 11 HD MIKE 11 phần mềm kỹ thuật chuyên dụng mô lưu lượng, chất lượng nước vận chuyển bùn cát sông, hệ thống tưới, kênh dẫn hệ thống dẫn nước khác Với môi trường đặc biệt thân thiện với người sử dụng, linh hoạt tốc độ̣, MIKE 11 cung cấp môi trường thiết kế hữu hiệu kỹ thuật cơng trình, tài ngun nước, quản lý chất lượng nước ứng dụng quy hoạch MIKE 11 mơ hình động lực chiều, mơ đun mơ hình thủy động lực (HD) phần trung tâm hệ thớng mơ hình MIKE 11 hình thành sở cho hầu hết mơ đun bao gờm tính tốn, dự báo lũ, tải khuếch tán, chất lượng nước mô đun vận chuyển bùn cát Mô đun MIKE 11 HD giải phương trình tổng hợp theo phương đứng để đảm bảo tính liên tục bảo tồn động lượng Phương trình mơ hình để tính tốn cho trường hợp dịng khơng ổn định phương trình liên tục phương trình động lượng (hệ phương trình Saint - Venant) với giả thiết chất lỏng không nén đồng nhất, dòng chảy chủ yếu chiều, độ dốc đáy nhỏ, thông số mặt cắt ngang biến động theo chiều dọc, phân bố áp suất thủy tĩnh cho phương trình liên tục (bảo toàn khối lượng) phương trình động lượng (bảo tồn động lượng) - Phương trình liên tục: wQ wA q wx wt - Phương trình động lượng: QQ wQ w § Q · wh wt wx ăă D A áá gA wx g C AR â Trong ú: h - Mực nước thời đoạn tính tốn (m); t - Thời gian tính tốn (s); Q - Lưu lượng dịng chảy qua mặt cắt (m3/s); x - Không gian (dọc theo chiều dịng chảy) (m); A - Diện tích mặt cắt ướt (m2); q - Lưu lượng đơn vị gia nhập 18 TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 10 - 2016 theo chiều dài dọc sông (m2/s); C - Hệ số Chezy; n - Hệ số nhám; R - Bán kính thủy lực (m); y Hệ số, theo Manning y = 1/6; g - Gia tốc trọng trường, g = 9,81 m/s2; α - Hệ số động lượng Thiết lập mơ hình tốn dự báo lũ sơng Lại Giang 3.1 Thiết lập mơ hình mưa - dịng chảy MIKE NAM Mơ hình MIKE NAM thiết lập nhằm xác định lượng dòng chảy đầu vào dịng chảy khu cho mơ hình thuỷ lực MIKE 11 Để thiết lập MIKE NAM cho lưu vực Lại Giang tiểu lưu vực phân chia dựa mơ hình số độ cao DEM kết hợp với đồ địa hình tỷ lệ lớn [2] Kết phân chia lưu vực sông Lại Giang xác định tiểu lưu vực, tiểu lưu vực nhỏ có diện tích 27 km2 lớn 571 km2 (hình1) Số liệu mưa trạm An Hoà, Bồng Sơn Hoài Ân dùng hiệu chỉnh, kiểm định mơ hình tính tốn biên trên, biên nhập lưu khu cho mơ hình thuỷ lực Thời gian cho hiệu chỉnh mơ hình áp dụng cho mùa lũ năm 1987 1990, năm có xuất lũ lớn có số liệu đo đạc đầy đủ, đồng trạm Trạm An Hoà sử dụng làm trạm kiểm tra xác định thơng số mơ hình trạm có đo lưu lượng dòng chảy nhiều năm lưu vực Việc hiệu chỉnh thơng số mơ hình chủ yếu tiến hành phương pháp thử sai Kết hiệu chỉnh (bảng 1) cho thấy tính tốn đường thực đo tương đối phù hợp với sai số lệch đỉnh giá trị thời gian không nhiều Hình thể kết hiệu chỉnh năm 1987 cho thấy hai đường q trình tính tốn thực đo tương đối bám sát pha dao động giá trị đỉnh Ngoài số NASH tương đối tốt, lớn 0,8 hệ số tương quan cao 0,9 Trong điều kiện hạn chế số liệu kết hồn tồn có độ tin cậy cao mơ hình ứng dụng tính tốn cho bước NGHIÊN CỨU & TRAO ĐỔI Bảng Thống kê kết đánh giá sai số hiệu chỉnh kiểm định MIKE NAM năm 1987, 1990 1992 trạm thuỷ văn An Hồ (TT: tính tốn; TĐ: thực đo) Các yu t Q max (m3/s) ăQ (m3/s) Sai s thi gian xuҩt hiӋn ÿӍnh HӋ sӕ NASH HӋ sӕ tѭѫng quan TT TĈ 1987 919 1000 81 1h 1h 0,84 0,97 TT TĈ 1990 670 683 13 1h 1h 0,85 0,95 Sau bước hiệu chỉnh mơ hình cho kết tốt, mơ hình MIKE NAM tiến hành kiểm định cho mùa lũ năm 1992 Kết kiểm định trình bày bảng hình Kết kiểm định thuỷ văn cho lưu vực An Hoà tương đối khả quan, nhận thấy đường tính tốn đường thực đo tương đồng Chênh lệch lưu lượng lớn tính tốn giá trị thực đo khơng đáng kể Sai số lệch đỉnh trạm kiểm tra nằm phạm vi cho phép Kết tính tốn hệ số NASH tương đối tốt (0,82), với kết thơng số mơ hình TT TĈ 1992 1291 1440 149 1h 1h 0,82 0,9 có độ tin cậy cao áp dụng vào dự báo Như vậy, thông qua hai bước hiệu chỉnh kiểm định, nghiên cứu xác đinh thông số tối ưu cho lưu vực lớn Lại Giang với giá trị thơng số CQOF = 0,79; CK1,2 = 18,20 ; Lmax = 80; Umax = 10 Tuy nhiên, thông số mùa lũ nên áp dụng tính tốn cho trận lũ, nhóm tác giả tiến hành hiệu chỉnh số thông số nhạy CQOF, CK1,2, TOF TIF cách tăng giảm thông số tuỳ thuộc vào đặc điểm trận lũ Hình Q trình lũ tính tốn thực đo hiệu chỉnh mơ hình MIKE NAM năm 1987 kiểm định năm 1992 3.2 Thiết lập mơ hình thuỷ lực MIKE 11 HD Mạng lưới sông nghiên cứu số hố mơ hình với 100 điểm số hố (hình 3) Tổng số mặt cắt sử dụng 21 mặt cắt, với tổng chiều dài 21,8 km, mặt cắt xã Ân Thạnh, huyện Hoài Ân mặt cắt cuối xã Hoài Sơn, huyện Hồi Nhơn Trong mơ hình thuỷ lực biên biên nhập lưu khu tính tốn thơng qua mơ hình thuỷ văn Hình Mạng thuỷ lực sơng Lại Giang số hố MIKE NAM, biên q trình mực nước mơ hình MIKE 11 HD triều thực đo có thời gian TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 10 - 2016 19 NGHIÊN CỨU & TRAO ĐỔI Thời gian dùng để hiệu chỉnh mơ hình mùa lũ năm 1987 1990 Ngoài ra, mùa lũ năm 1992 dùng để kiểm định lại thơng số mơ hình MIKE 11 Biên mực nước triều đo đạc trạm Quy Nhơn có thời gian Đặc biệt trạm thủy văn Bồng Sơn sử dụng làm nút kiểm tra nhằm dị tìm thơng số tối ưu cho mơ hình tính hệ số NASH, hệ số tương quan để đánh giá độ tin cậy mô hình Bảng Thống kê kết đánh giá sai số hiệu chỉnh kiểm định MIKE 11 năm 1987, 1990 1992 trạm thuỷ văn Bồng Sơn (TT: tính tốn; TĐ: thực đo) Các yӃu tӕ H max (m /s) ăH (m3/s) H s NASH H s tѭѫng quan TT TĈ 1987 7,6 8,6 1,0 0,79 0,83 TT TĈ 1990 7,13 7,25 0,12 0,77 0,80 TT TĈ 1992 7,47 7,38 0,90 0,75 0, 86 Hình Q trình lũ thực đo tính tốn trạm Bồng Sơn hiệu chỉnh năm 1987 1990 Kết đánh giá sai số trận lũ hiệu chỉnh 1987, 1990 kiểm định năm 1992 trạm Bồng Sơn thể bảng hình 4, cho thấy tính tốn thực đo tương đối sát Thời gian giá trị xuất đỉnh lũ thực đo tính tốn trùng khớp Tuy nhiên, lũ năm 1990 chênh lệch mực nước lớn tính tốn giá trị thực đo phần chân lũ lớn (1 m) sai số tính tốn thực đo tính 0,79 Tương tự, lũ năm 1992, chênh lệch mực nước lớn tính tốn giá trị thực đo 0,12 m sai số tính tốn thực đo (NASH) 0,75 Ngun nhân dẫn đến sai số trình hiệu chỉnh kiểm định mơ hình liệu đầu vào cho toán thuỷ lực xác định từ mơ hình mưa – dịng chảy, bên cạnh số lượng mặt cắt dùng cho tính tốn thuỷ lực cịn chưa đảm bảo nên khơng tránh khỏi có sai số Tuy 20 TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 10 - 2016 nhiên đánh giá độ tin cậy mơ hình thơng qua hệ số NASH hệ số tương quan mức đạt (> 0,7) Trong điều kiện lưu vực thiếu số liệu quan trắc việc dự báo đỉnh lũ đóng vai trị quan trọng kết chấp nhận để nghiên cứu dự báo Như vậy, qua q trình dị tìm thơng số tối ưu cho mơ hình MIKE 11, nghiên cứu xác định thơng số nhám (manning) cho đoạn sông hạ lưu sông Lại Giang giao động từ 0,025 - 0,032 Lập phương án dự báo dự báo thử nghiệm lưu vực sông Lại Giang 4.1 Lập phương án dự báo Dựa việc phân tích đặc điểm tự nhiên lưu vực, đặc điểm lũ, tình hình số liệu mơ hình tốn thiết lập, phương án dự báo lũ đưa nghiên cứu có thời gian dự kiến 24h vị trí dự báo lựa chọn trạm thuỷ văn Bồng Sơn Phương án dự báo NGHIÊN CỨU & TRAO ĐỔI sau, mơ hình MIKE NAM sử dụng số liệu giả định lại sau bước dự báo có thêm kết thiết mưa dự báo tính tốn lưu lượng dịng chảy dự báo thực đo sau làm đầu vào nhập lưu cho mơ hình thuỷ lực MIKE 11 diễn tốn dịng chảy đến vị trí dự báo Cụ thể, thời điểm bắt đầu làm dự báo, dựa vào mưa q trình dịng chảy thời đoạn trước tiến hành dự báo lũ đến biên trên, nhập lưu khu cho 24h sau 24h coi mưa tồn lưu vực mưa dự báo Ở nhóm tác giả tiến hành giả thiết lượng mưa dự báo Sau có giá trị Q dự báo đầu tiên, tiếp tục diễn toán thủy lực tới trạm Bồng Sơn để Hình Quá trình lũ thực đo tính tốn trạm Bồng Sơn kiểm định năm 1992 kết dự báo mực nước Tương 4.2 Dự báo thử nghiệm đánh giá sai số tự, thời điểm 24h tiếp theo, sau biết Trong phần này, phương án dự báo đưa mưa xảy thời đoạn trước tiếp hành cập nhật lại số liệu mơ hình tiếp tục dự báo cho dự báo thử nghiệm để đánh giá mức bảo 24h sau với giả thiết mưa dự báo khác đảm phương án Trong đó, hai trận lũ năm 1990 trận lũ năm 1992 sử dụng để dự báo thử tiếp tục bước kết thúc trận lũ Các mơ hình dự báo có sai số nhiều nghiệm Trận lũ năm 1990 xảy ngày từ ngun nhân như: sai số mơ hình, sai số ngày 11/11/1990 - 16/11/1990 Căn vào mưa liệu đầu vào, sai số thông số… Nghiên cứu tình hình số liệu, tiến hành dự báo lũ từ sử dụng công cụ thống kê để hiệu chỉnh 7:00:00 sáng ngày 11/11/1990 - 7:00:00 sáng sai số dự báo dựa việc xây dựng ngày 16/11/1990 Tương tự, trận lũ tháng phương trình hồi quy tuyến tính kết mơ 10/1992 có thời gian kéo dài ngày từ ngày thực đo Sau có kết dự báo, 22/10/1992 - 30/10/1992 Theo đó, kết dự số liệu dự báo hiệu chỉnh phương báo thể hình hình trình hồi quy phương trình xác Bảng Đánh giá phương án dự báo thử nghiệm trận lũ năm 1992 Mӭc ÿҧm HӋ sӕ bҧo tѭѫng phѭѫng án quan dӵ báo Tên trҥm Dӵ báo chѭa hiӋu Bӗng chӍnh Sѫn Dӵ báo ÿã hiӋu chӍnh Ĉánh giá phѭѫng án 0,67 0,74 0,67 Kém 0,45 0,90 0,89 Ĉҥt Từ kết dự báo thấy sau áp dụng Mô đun hiệu chỉnh sai số, kết cải thiện đáng kể với mức đảm bảo phương án đạt 0,83 phương án đánh giá mức đạt cho trận lũ dự báo năm 1990 (hình 6) Ngoài ra, trận lũ năm 1992 mức đảm bảo phương án 0,89 phương án đánh giá đạt (hình 7) Từ kết dự báo thấy giá trị mực nước dự báo sau hiệu chỉnh xấp xỉ mực nước thực đo so với giá trị mực nước dự báo chưa hiệu chỉnh Bên cạnh đó, giá trị dự báo thực đo có thời gian xuất đỉnh, đồng thời khơng có biến động lớn đường nước lên nước xuống chúng tương đối đồng dạng TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 10 - 2016 21 NGHIÊN CỨU & TRAO ĐỔI 6.80 6.30 5.80 5.30 4.80 4.30 3.80 3.30 2.80 2.30 1.80 H(m) Thӡi gian XI/11/1990 XI/12/1990 XI/13/1990 XI/14/1990 XI/15/1990 XI/16/1990 7:00:00 AM 7:00:00 AM 7:00:00 AM 7:00:00 AM 7:00:00 AM 7:00:00 AM H dӵ báo chѭa hiӋu chӍnh H dӵ báo ÿã hiӋu chӍnh H thӵc ÿo Hình Đường trình mực nước thực đo dự báo trạm Bồng Sơn trận lũ tháng 11/ 1990 6.80 6.30 5.80 5.30 4.80 4.30 3.80 3.30 2.80 2.30 1.80 H (m) Thӡi gian X/22/1992 X/23/1992 X/24/1992 X/25/1992 X/26/1992 X/27/1992 7:00:00 AM 7:00:00 AM 7:00:00 AM 7:00:00 AM 7:00:00 AM 7:00:00 AM H dӵ báo chѭa hiӋu chӍnh H dӵ báo ÿã hiӋu chӍnh H thӵc ÿo Hình Đường trình mực nước thực đo dự báo trạm Bồng Sơn trận lũ tháng 10/1992 Kết luận Nghiên cứu xây dựng phương án dự báo lũ lớn với thời gian dự kiến 24h cho lưu vực sông Lại Giang điều kiện hạn chế số liệu thuỷ văn địa hình sơng Ngồi ra, nghiên cứu này, phương pháp mơ hình tốn thuỷ văn sử dụng linh hoạt với kết hợp mơ hình mưa dịng chảy (MIKE NAM) mơ hình thuỷ lực diễn tốn dịng chảy (MIKE 11) hiệu chỉnh kiểm định tốt với hệ số NASH đạt yêu cầu, dao động từ 0,75 - 0,85 Kết dự báo thử nghiệm cho thấy phương án dự báo đề với thời gian dự kiến 24 có mức bảo đảm cao (89%) phương án đánh giá đạt Tài liệu tham khảo Bạch Huy Hoàng (2012), Nghiên cứu xây dựng đồ ngập lụt đánh giá rủi ro miền Trung Việt Nam-Đối tượng nghiên cứu: Lưu vực sơng Lại Giang, tỉnh Bình Định, Việt Nam Luận văn Thạc sĩ, Viện Khoa học Thuỷ lợi Việt Nam Nguyễn Văn Hiếu (2013),Ứng dụng mơ hình Mike-Flood tính tốn ngập lụt mơ phương án xây dựng đê cho sơng Lại Giang, Tạp chí Khoa học ĐHQGHN, Khoa học Tự nhiên Công nghệ 29, Số 1S (2013) 72-79 Nguyễn Ân Niên, Vũ Văn Nghị, Nguyễn Anh Đức (2002), Lũ lưu vực Lại Giang sơ đồ tính, Tuyển tập KHCN năm 2002, Viện Khoa học Thuỷ Lợi Việt Nam http://www.baobinhdinh.com.vn/binhdinhquaanh/2009/5/75957/ 22 TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 10 - 2016