Ứng dụng mô hình phần tử hữu hạn sóng động học đánh giá tác động của việc sử dụng đất đến dòng chảy lưu vực sông Trà Khúc - trạm Sơn Giang

MỤC LỤC Trang Mở đầu………………………………………………………………………………………………… Chƣơng Tổng quan mơ hình tốn thuỷ văn…………………………………………… 1.1 Khái niệm mơ hình tốn thuỷ văn phân loại………………………………… 1.2 Giới thiệu số mơ hình tốn thuỷ văn………………………………………… 1.3 Một số mơ hình ứng dụng qui hoạch lƣu vực sơng…………………… 11 1.4 Mơ hình phần tử hữu hạn sóng động học đánh giá tác động việc sử dụng đất lƣu vực đến dòng chảy………………………………… Chƣơng Điều kiện địa lý tự nhiên lƣu vực sông Trà Khúc - trạm Sơn Giang 13 25 2.1 Vị trí địa lý………………………………………………………………………………… 25 2.2 Địa hình……………………………………………………………………………………… 25 2.3 Địa chất thổ nhƣỡng……………………………………………………………………… 29 2.4 Thảm thực vật……………………………………………………………………………… 29 2.5 Khí hậu……………………………………………………………………………………… 30 2.6 Mạng lƣới sơng suối tình hình nghiên cứu thuỷ văn……………………… 31 Chƣơng Ứng dụng mơ hình phần tử hữu hạn sóng động học đánh giá tác động việc sử dụng đất đến dòng chảy lƣu vực sông Trà Khúc - trạm Sơn Giang 35 3.1 Mô tả xử lý số liệu…………………………………………………………………… 35 3.2 Chƣơng trình tính………………………………………………………………………… 50 3.3 Kết tính tốn………………………………………………………………………… 52 Kết luận………………………………………………………………………………………………… 63 Tài liệu tham khảo………………………………………………………………………………… 65 Phụ lục………………………………………………………………………………………………… 66 MỞ ĐẦU Lũ pha chế độ nƣớc năm thƣờng gây tác hại lớn ngƣời việc tìm hiểu, khống chế đƣợc lũ vấn đề thời đƣợc quan tâm từ trƣớc tới Dự báo lũ biện pháp phòng chống hữu hiệu nhất, hạn chế nguồn tƣ liệu thơng tin khí tƣợng mặt đệm nên từ trƣớc tới dự báo lũ thƣờng xác mƣa xuất thành dòng chảy trạm thƣợng nguồn để áp dụng dự báo theo phƣơng pháp mực nƣớc tƣơng ứng Có nghĩa sử dụng lƣợng thơng tin lịng dẫn Để khắc phục nhƣợc điểm để tạo thời gian dự kiến dài thƣờng dùng mơ hình tốn Các mơ hình tốn mƣa - dòng chảy từ trƣớc tới thƣờng dùng mơ hình thơng số tập trung, mà nhƣợc điểm khơng tính đƣợc hết thay đổi đa dạng mặt đệm lƣu vực Nhằm khắc phục điều này, điều kiện công nghệ GIS phát triển, khố luận sử dụng phƣơng pháp mơ dựa phƣơng pháp SCS phƣơng pháp phần tử hữu hạn để sử dụng thông tin đa dạng mặt đệm với số liệu khí tƣợng thuỷ văn đồ số Khoá luận chọn lƣu vực sông Trà Khúc - trạm Sơn Giang làm đối tƣợng nghiên cứu để thực hai nhiệm vụ chính: - Xác định thơng số ổn định để mô lũ làm sở cho việc thiết lập phƣơng án, cảnh báo lũ phục vụ cho việc phòng chống thiên tai lũ lụt - Đánh giá đƣợc mức độ ảnh hƣởng việc sử dụng đất tới việc hình thành dịng chảy làm sở cho kiến nghị cho việc qui hoạch lƣu vực hợp lý Do hạn chế kiến thức phân tích tổng hợp, thời gian nghiên cứu nên kết cịn có nhiều khiếm khuyết mong góp ý tận tình thầy, bạn đồng nghiệp để vấn đề đƣợc giải ngày tốt Chƣơng TỔNG QUAN MƠ HÌNH TỐN THUỶ VĂN 1.1 KHÁI NIỆM MƠ HÌNH TỐN THỦY VĂN VÀ PHÂN LOẠI 1.1.1 Khái niệm mơ hình tốn thuỷ văn Mơ hình hệ thống thủy văn mơ hình vật lý hay tốn học Mơ hình vật lý mơ hình biểu thị hệ thống thật dƣới dạng thu nhỏ, nhƣ mơ hình thủy lực đập tràn Mơ hình tốn học miêu tả hệ thống dƣới dạng toán học, tập hợp phƣơng trình tốn, mệnh đề logic thể quan hệ biến thông số mơ hình để mơ hệ thống tự nhiên, hay nói cách khác mơ hình tốn học hệ thống biến đổi đầu vào (hình dạng, điều kiện biên, lực v.v ) thành đầu (tốc độ chảy, mực nƣớc, áp suất v.v ).[1, 8] 1.1.2 Phân loại mơ hình tốn thuỷ văn Có nhiều cách phân loại mơ hình tùy theo quan điểm ý tƣởng ngƣời phân loại Một cách phân loại dựa sở xem xét phân bố biến vào hệ thống trƣờng không gian, thời gian Đối với mơ hình, ngƣời ta xem xét định sau: Các biến mơ hình có ngẫu nhiên khơng? Chúng biến đổi theo không gian nhƣ nào? Chúng biến đổi theo thời gian sao? Tuỳ thuộc lựa chọn định trên, mơ hình tốn thuỷ văn phân loại theo “cây phân loại” nhƣ hình 1.1 Mơ hình tốnThuỷ văn Mơ hình tất định Mơ hình ngẫu nhiên Mơ hình thơng số tập trung Mơ hình hộp đen Mơ hình thơng số phân phối Mơ hình quan niệm Hình 1.1 Sơ đồ phân loại mơ hình tốn thuỷ văn Mơ hình vật lý - tốn Mơ hình động lực - ngẫu nhiên Trong sơ đồ phân loại hình 1.1, mơ hình tốn thuỷ văn đƣợc phân loại thành: mơ hình tất định mơ hình ngẫu nhiên Mơ hình ngẫu nhiên mơ trình dao động thân trình thủy văn mà không ý đến nhân tố đầu vào tác động hệ thống [8] Mơ hình tất định mơ hình mơ q trình biến đổi tƣợng thuỷ văn lƣu vực mà ta biết trƣớc Xét quan điểm hệ thống, mơ hình thuỷ văn tất định có thành phần sau: - Đầu vào hệ thống - Hệ thống - Đầu hệ thống Đầu vào (I) Đầu (Q) Hệ thống Một phân loại mơ hình thuỷ văn tất định khác đƣợc trỡnh by hỡnh 1.2 [1] Mô hình tất định (Deterministic models) Mô hình thuỷ động lực học (Hydro-dynamical models) Mô hình nhận thức (Conceptual models) Mô hình thông số dải (Distruibuted models) Phân phối theo đơn vị diện tích nhỏ (l-ới tính km2) Mô hình hộp đen (Black-box models) Mô hình thông số tập trung (Concentrated models) Phân phối theo đơn vị diện tích lớn (tiểu vùng thuỷ văn) Hỡnh 1.2 S phõn loi cỏc mụ hỡnh thuỷ văn tất định * Các mơ hình tiêu biểu nhƣ TANK, SSARR 1.2 GIỚI THIỆU MỘT SỐ MƠ HÌNH TỐN THUỶ VĂN (MƠ HÌNH MƯA - DỊNG CHẢY) 1.2.1 Mơ hình Trung tâm khí tƣợng thuỷ văn Liên Xơ (HMC) [1] Mơ hình mơ q trình tổn thất dịng chảy lƣu vực sau ứng dụng cách tiệm cận hệ thống để diễn tốn dịng chảy tới mặt cắt cửa lƣu vực Lƣợng mƣa hiệu sinh dòng chảy mặt P đƣợc tính tƣ phƣơng trình: P=h-E-I (1.1) đó: h - Cƣờng độ mƣa thời đoạn tính tốn (6h, 24h, ); E - Lƣợng bốc thoát nƣớc; I - Cƣờng độ thấm trung bình Mơ hình có tính đến lƣợng bốc mà lƣợng bốc lƣu vực thiếu nhiều, chủ yếu đƣợc ƣớc tính từ phƣơng trình hạn chế điều kiện đo đạc để xác định trực tiếp lƣợng bốc Ngoài cƣờng độ thấm trung bình thƣờng đƣợc lấy trung bình cho tồn lƣu vực với thời gian khơng xác định nên mơ hình cịn nhiều hạn chế 1.2.2 Mơ hình SSARR [1,8,10] Mơ hình SSARR Rockwood D xây dựng từ năm 1957, gồm thành phần bản: - Mơ hình lƣu vực - Mơ hình điều hồ hồ chứa - Mơ hình hệ thống sơng Trong mơ hình lƣu vực, phƣơng trình SSARR sử dụng để diễn tốn dịng chảy lƣu vực luật liên tục phƣơng pháp trữ nƣớc áp dụng cho hồ thiên nhiên sở phƣơng trình cân nƣớc:  I1  I   O1  O    t    t  S  S1     Phƣơng trình lƣợng trữ hồ chứa : (1.2) dS dQ  Ts dt dt (1.3) Mơ hình SARR cho phép diễn tốn tồn lƣu vực nhƣng bên cạnh mơ hình SSARR cịn hạn chế áp dụng đƣợc với lƣu vực không lớn, với lƣu vực có điều kiện ẩm khơng đồng tính tốn cho kết mơ khơng xác Mơ hình khơng thể sử dụng cách trực tiếp để điều tra (kiểm tra tác động thủy văn việc thay đổi đặc điểm lƣu vực sơng ví dụ nhƣ kiểu thảm thực vật, việc bảo vệ đất hoạt động quản lý đất tƣơng tự khác) 1.2.3 Mơ hình TANK [8, 10] Mơ hình TANK đƣợc phát triển Trung tâm Nghiên cứu Quốc gia phòng chống thiên tai Tokyo, Nhật Bản Theo mơ hình này, lƣu vực đƣợc mô chuỗi bể chứa (TANKS) theo tầng phù hợp với phẫu diện đất Nƣớc mƣa tuyết tan đƣợc quy bể chứa Mỗi bể chứa có cửa đáy hai cửa cuối thành bể, phía đáy Lƣợng nƣớc chảy khỏi bể chứa qua cửa đáy vào bể chứa tầng sau trừ bể chứa tầng cuối, bể lƣợng chảy xuống đƣợc xác định tổn thất hệ thống Lƣợng nƣớc qua cửa bên bể chứa trở thành lƣợng nhập lƣu cho hệ thống lòng dẫn Số lƣợng bể chứa, kích thƣớc nhƣ vị trí cửa thơng số mơ hình Hệ thức mơ hình Mƣa bình quân lƣu vực (P) n P n  Wi x1 /  Wi i 1 (1.4) i 1 đó: n - số điểm đo mƣa; Xi - lƣợng mƣa điểm thứ i; Wi - trọng số điểm mƣa thứ i Theo M.Sugawara Wi đƣợc chọn bốn số sau: 0,25; 0,5; 0,75; 1,0 Bốc lƣu vực (E) Khi XA  PS  E  Khi XA  PS  E   0,8EVT 0,75(0,8EVT  h f )  h f  E   0,6EVT  va XA  PS  H f  XA  PS (1.5) Cơ cấu truyền ẩm bể chứa đƣợc chia làm hai phần: dƣới, chúng xảy trao đổi ẩm Tốc độ truyền ẩm từ dƣới lên T xuống T2 đƣợc tính theo công thức: XA )TB PS XS T2  TC0  (1  )TC SS T1  TB0  (1  (1.6) (1.7) đó: XS, SS - lƣợng ẩm thực lƣợng ẩm bão hoà phần dƣới bể A; TBo,TB, TCo, TC - thông số truyền ẩm, theo MSugawar chúng nhân giá trị: TB = TB0 = mm/ngày đêm TC = 1mm/ngày đêm TC0 = 0,5mm/ngày đêm Dòng chảy từ bể A Lƣợng nƣớc vào bể A mƣa (P) Dòng chảy qua cửa bên(YA1, YA2) đáy (YA0) đƣợc xác định theo công thức sau: Hf XA + P - PS (1.8) YA0 = HfA0 (1.9)  (H f  HA1 ); H f  HA1 YA1   H f  HA1 0 (1.10) Trong mơ hình, tác dụng điều tiết sƣờn dốc tự động đƣợc xét thông qua bể chứa xếp theo chiều thẳng đứng Nhƣng hiệu tác động không đủ mạnh coi tổng dịng chảy qua cửa bên bể YA2+YA1+YB2+YC1+YD1 lớp cấp nƣớc điểm Đây hạn chế mô hình TANK 1.2.4 Mơ hình NAM [1] Mơ hình NAM đƣợc 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 Mơ hình dựa nguyên tắc bể chứa theo chiều thẳng đứng hồ chứa tuyến tính Trong mơ hình NAM, lƣu vực đƣợc xem đơn vị xử lý Do đó, thơng số biến đại diện cho giá trị đƣợc trung bình hố 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 năm bể chứa riêng biệt có tƣơng tác lẫn nhau: + Bể chứa tuyết đƣợc kiểm soát điều kiện nhiệt độ khơng khí + Bể chứa mặt bao gồm lƣợng ẩm bị chặn lớp phủ thực vật, lƣợng điền trũng lƣợng ẩm tầng sát mặt Umax giới hạn lƣợng nƣớc bể + Bể chứa tầng dƣới vùng dễ mà từ cối rút nƣớc cho bốc thoát Lmax giới hạn lƣợng nƣớc bể + Bể chứa nƣớc tầng ngầm bể chứa nƣớc tầng ngầm dƣới hai bể chứa sâu Dòng chảy tràn dòng chảy sát mặt đƣợc diễn toán qua hồ chứa tuyến tính thứ nhất, sau thành phần dịng chảy đƣợc cộng lại diễn toán qua hồ chứa tuyến tính thứ hai Cuối thu đƣợc dịng chảy tổng cộng cửa Phƣơng trình mơ hình: Dịng chảy sát mặt QIF: L   CLIF  L max CQIF U Víi   CLIF QIF     Khi   L  CLIF L max (1.11) L  CLIF L max đó: CQIF - hệ số dịng chảy sát mặt; CLIF - ngƣỡng dịng chảy; U, Lmax thơng số khả chứa Dòng chảy tràn QOF: 10 L   CLOF  L max CQOF PN   CLOF QOF       Víi Khi L  CLOF L max (1.12) L  CLOF L max đó: CQOF - hệ số dịng chảy tràn; CLOF - ngƣỡng dịng chảy Trong tính tốn giả thiết dịng chảy khỏi hồ tn theo quy luật đƣờng nƣớc rút: Q out  Q e out  t CK t    1  e Ck   Q in     (1.13) đó: Q dịng chảy tính thời điểm trƣớc; Qin dịng chảy vào thời out điểm tính; CK số thời gian hồ chứa Mơ hình NAM tính đƣợc dịng chảy sát mặt dịng chảy tràn, song bên cạnh thơng số biến đƣợc tính trung bình hố cho tồn lƣu vực Nên việc cụ thể hố tính tốn cho đơn vị nhỏ lƣu vực bị hạn chế 1.3 MỘT SỐ MƠ HÌNH ỨNG DỤNG TRONG QUI HOẠCH LƢU VỰC SƠNG 1.3.1 Mơ hình USDAHL [1,7] Mơ hình đƣợc cơng bố vào năm 70, mơ hình thơng số dải theo tiểu vùng thuỷ văn Mơ hình chia bề mặt lƣu vực thành tiểu vùng thuỷ văn với đặc trƣng nhƣ loại đất, sử dụng đất vùng, trình nhƣ mƣa, bốc hơi, thấm, điền trũng, dịng chảy đƣợc tính toán xử lý mối liên kết vùng với vùng khác Q trình hình thành dịng chảy đƣợc mơ nhƣ sau: Dịng chảy mặt bao gồm q trình thấm, trình trữ chảy tràn Quá trình thấm đƣợc mơ phƣơng trình Holtan: ft  A GI S1.4  f c at (1.14) đó: ft - cƣờng độ thấm; A - hệ số phụ thuộc vào độ rỗng đất, mật độ rễ cây; GI- số phát triển thực vật, phụ thuộc vào nhiệt độ khơng khí loại cây; fc 11 cƣờng độ thấm ổn định; Sat- độ thiếu hụt ẩm đất hàm số theo thời gian: Sat  Sat -1 - f t -1  f c Quá trình trữ, chảy tràn đƣợc thực dựa sở phƣơng trình cân nƣớc Q trình dịng chảy dƣới mặt đất đƣợc xem xét dựa sở phƣơng trình cân độ ẩm đất Dịng chảy lịng dẫn đƣợc diễn tốn theo mơ hình tuyến tính Sơ đồ cấu trúc mơ hình USDAHL đƣợc thể hình 1.3 Mơ hình có khả đánh giá tác động yếu tố lƣu vực quy mơ trung bình đến hình thành dịng chảy Mơ hình USDAHL xét đến tất thành phần phƣơng trình cân nƣớc, thành phần đƣợc xử lý xem xét dựa phƣơng trình Song việc xử lý lƣợng thấm, bốc hơi, điền trũng gặp nhiều khó khăn ngồi với lƣu vực lớn khả đánh giá tác động yếu tố lƣu vực đến hình thành dịng chảy Gi¸ng thủ Thấm trực tiếp Trữ mặt Dòng chảy mặt Bốc thoát Thấm xuống tầng d-ới Dòng chảy d-ới mặt Trữ ẩm đất Diễn toán Trữ n-ớc ngầm lòng dẫn Thấm tầng sâu Dòng chảy Hỡnh 1.3 S cu trúc mơ hình USDAHL 1.3 Mơ hình HEC-1 [8] Về ngun tắc mơ hình tiến hành giải thành phần: 12 Cứ tiếp tục hết dải sông I tiếp tục nhập sơng II tƣơng tự (Phụ lục3) 3.3 KẾT QUẢ TÍNH TỐN Từ file số liệu thơng số đƣợc xác lập qua việc tối ƣu hoá theo chƣơng trình cho trận lũ năm 1998 1999 cho kết mơ dịng chảy (Hình 3.5, 3.6, 3.7, 3.8, 3.9, 3.10, 3.11) chạy chƣơng trình tính cho chuỗi số liệu độc lập (Hình 3.12, 3.13) nhƣ sau: Qdb KQ412XI-98 Q (m^3/s) Qtd 1600 1400 1200 1000 800 600 400 200 Ngay 10 15 Ngày 10 11 12 R2 Sai số đỉnh Sai số tổng lƣợng Qdb (m3/s) Qtd (m3/s) 150 155.3 1241.41 609.62 263.58 315.18 219.72 185.81 257.59 92% 12,8% 150 259 1400 896 536 410 335 288 264 33,5% Hình 3.5 Kết mơ trận lũ từ ngày 4/XI đến 12/XI năm 1998 trạm Sơn Giang Qdb Kd21725XI98 Q (m^3/s) Qtd 12000 10000 8000 6000 4000 2000 Ngay 15 17 19 21 23 25 Ngày 17 18 19 20 21 22 23 24 25 R2 Sai số đỉnh Sai số đỉnh Sai số tổng 55 Qdb (m3/s) Qtd (m3/s) 753 753 4465.09 7892.75 1939.01 10335.13 1931.57 927.57 1224.53 74,1% 22,7% 35% 19% 753 600 2450 6100 2940 6720 2450 1400 1080 Hình 3.6 Kết mơ trận lũ từ ngày lƣợng 17/XI đến 25/XI năm 1998 trạm Sơn Giang Q (m^3/s) Qdb Qtd Kd22530XI-98 4500 4000 3500 3000 2500 2000 1500 1000 500 27 28 29 30 R2 Sai số đỉnh Sai số tổng lƣợng Ngay 24 26 28 30 Ngày 25 26 32 Qdb (m3/s) Q(i) (m3/s) 1080 3799.5 1080 1322.1 1164.3 1121.4 1102.9 94,7% 10,7% 3510 1180 874 722 603 16,9% Hình 3.7 Kết mô trận lũ từ ngày 25/XI đến 30/XI năm 1998 trạm Sơn Giang Qdb Q (m^3/s) X2227-99 Qtd 3500 3000 2500 2000 1500 1000 500 Ngay 20 22 24 26 28 Hình 3.8 Kết mơ trận lũ từ ngày 22/X đến 27/X năm 1999 trạm Sơn Giang 56 Ngày 22 23 24 25 26 27 R2 Sai số đỉnh Sai số tổng lƣợng Qdb (m3/s) Qtd (m3/s) 357 1152.91 2871.72 1076.82 971.05 776.64 94,3% 17,1% 8% 357 1450 2380 1070 722 695 Q (m^3/s) Qdb XI18-99 Qtd 9000 8000 7000 6000 5000 4000 3000 2000 1000 Ngay 10 Ngày R2 Sai số đỉnh Sai số đỉnh Sai số tổng lƣợng Hình 3.9 Kết mô trận lũ từ ngày Qdb (m3/s) Qtd (m3/s) 1555.09 3736.7 6216.88 6277.82 8169.62 2526.46 1168.16 1054.68 63% 35% 967 3710 4050 3020 5460 3090 1360 965 33,2% 26,4% 1/XI đến 8/XI năm 1999 trạm Sơn Giang Q (m^3/s) Qdb Qtd K1218XII-99 3000 Ngày 12 13 2500 14 2000 15 1500 16 17 1000 18 500 R2 Ngay 10 12 14 16 18 Sai số đỉnh Sai số đỉnh 20 Qdb (m3/s) Qtd (m3/s) 720 847.23 1044.63 839.4 2433.25 968.63 720 62% 2,4% 51,9% Sai số tổng lƣợng Hình 3.10 Kết mô trận lũ từ ngày 12/XII đến 18/XII năm 1999 trạm Sơn Giang 57 16,9% 720 982 1020 797 1170 908 694 Qtd Qdb Kh1319XI-98 Q (m^3/s) 8000 Ngày 13 7000 25 6000 37 5000 49 4000 61 3000 73 2000 85 1000 97 Gio 109 0 50 100 150 200 121 133 145 R2 Sai số đỉnh Qdb (m3/s) Qtd (m3/s) 264 264 1480 264.05 4210 4711.01 4640 7061.47 2420 3932.72 2090 4042.21 1400 1992.09 1030 842.21 866 505.34 681 385.53 630 335.94 589 304.36 569 293.24 70,1% 34,3% Hình 3.11 Kết mơ trận lũ từ ngày 13/XI đến ngày 19/XI năm 1998 Sai số tổng lƣợng trạm Sơn Giang Qdb Kda913XII-98 Q (m^3/s) Qtd 4500 4000 3500 3000 2500 2000 1500 Ngày 10 11 12 13 R2 Sai số đỉnh 16,3% Qdb (m3/s) Qtd (m3/s) 896 896 1239.07 1900 3882.54 2970 1133.17 1550 1031.31 1000 86,4% 23,5% 1000 500 Ngay 11 13 15 Hình 3.12 Kết tính cho chuỗi số liệu độc Sai số tổng lƣợng lập từ ngày 9/XII đến ngày 13/XII năm 1998 trạm Sơn Giang 58 1,6% Qdb (m3/s) Qtd (m3/s) Ngày 569 569 569 984 13 569.01 2890 19 1471.78 4600 25 5118.29 4800 31 6186.69 7270 37 8449.09 7490 43 7921.83 4820 49 5792.2 4540 55 3445.75 3500 61 2358 2670 67 1545.86 2110 73 1477.09 2670 79 5698.94 6980 85 16730.15 9940 91 13159.11 6190 97 5277.22 3550 103 2754.68 2800 109 1796.67 2320 115 1249.2 1990 121 937.88 1690 127 845.02 1520 133 736.32 1360 139 719.31 1220 145 741.59 1130 R 69% Sai số đỉnh 11,4% Sai số đỉnh 40,6% Sai số tổng lƣợng 6,8% Qdb Qtd Kdh1825-98 Q (m^3/s) 18000 16000 14000 12000 10000 8000 6000 4000 Gio 2000 0 20 40 60 80 100 120 140 160 Hình 3.13 Kết tính cho chuỗi số liệu độc lập từ ngày 19/XI đến 19 ngày 25/XI năm 1998 trạm Sơn Giang Đánh giá sai số: Theo tiêu chuẩn đánh giá sai số tổ chức Khí tƣợng giới (WMO) thơng qua độ hữu hiệu đánh giá qua tiêu R2 Độ hữu hiệu đƣợc xác định nhƣ sau: R2  F02  F 100% F02 N N i 1 i 1  Trong đó: F   Qid  Qit 2 , F02   Qid  Q d  đó: Qiđ lƣu lƣợng thực đo; Qit lƣu lƣợng tính tốn, Qdtb lƣu lƣợng thực đo trung bình thời kỳ tính tốn; N tổng số điểm quan hệ lƣu lƣợng thực đo tính tốn 59 Tiêu chuẩn đánh giá nhƣ sau:  40  65% d ¹ t  R  65%  85% kh ¸   85% tèt  Theo tiêu chuẩn đánh giá sai số tổ chức Khí tƣợng giới với trận lũ chạy mô cho kết mô tả trận lũ thuộc loại tốt với R2>85%, trận lũ đạt loại trận lũ đạt Qua kết tính tốn so sánh dự báo với thực đo thấy thơng số dùng để mô cho trận lũ cho kết đạt áp dụng đƣợc vào lƣu vực sơng Trà Khúc Để có đƣợc thơng số tƣơng đối hồn chỉnh khố luận sử dụng hai trận lũ độc lập để kiểm định thông số Kết chạy cho hai trận lũ độc lập đƣợc thể hình 9.12 hình 9.13 Với trận lũ từ ngày 9/XII đến ngày 13/XII năm 1998 cho kết lƣợng tốt (với sai số tổng lƣợng