Giới thiệu dự án
Lũ lụt tại khu vực miền Trung Việt Nam luôn là hiểm họa thiên tai nghiêm trọng nhất đối với đời sống kinh tế - xã hội. Do đặc thù địa hình sườn Đông Trường Sơn hẹp, dốc (độ dốc trung bình toàn lưu vực đạt 25%), mạng lưới sông suối dạng nan quạt với mật độ 0.47 km/km², các trận mưa lớn sinh ra từ bão (chiếm 55%) và không khí lạnh tăng cường (chiếm 23%) tạo nên những cơn lũ đặc trưng miền núi: cường suất lũ lớn, thời gian truyền lũ ngắn, đỉnh lũ nhọn và biên độ dao động mực nước lớn. Lưu vực sông Thu Bồn tính đến trạm thủy văn Nông Sơn có diện tích tự nhiên $3155 \text{ km}^2$, sở hữu mô đun dòng chảy mùa lũ đạt tới $200 \text{ l/s}\cdot\text{km}^2$ – mức cao nhất trên toàn lãnh thổ Việt Nam.
+---------------------------------------+
| LƯU VỰC SÔNG THU BỒN - NÔNG SƠN |
| Diện tích: 3155 km2 | Độ dốc: 25% |
+---------------------------------------+
|
[ Mưa bão / Không khí lạnh cực đoan ]
|
+-------------------------+-------------------------+
| |
v v
[ Dòng chảy sườn dốc ] [ Tổn thất thấm SCS ]
- 192 Phần tử hữu hạn - Chỉ số CN & Nhám n
- Sóng động học 1D (KW1D) - Tổn thất ban đầu Ia=0.2S
| |
+-------------------------+-------------------------+
|
v
[ Dòng chảy trong 23 Lòng dẫn ]
- Diễn toán thủy lực 59 Dải
- Phương pháp phần tử hữu hạn Galerkin
|
v
[ CỬA RA: TRẠM NÔNG SƠN ]
- Đỉnh lũ Qmax & Tổng lượng W
- Hệ số hiệu quả R2 >= 79.82% - 90%
Vấn đề cốt lõi đặt ra là các mô hình dự báo thủy văn truyền thống dạng thông số tập trung không thể phản ánh chính xác tính bất đồng nhất về không gian của địa hình, thổ nhưỡng và hiện trạng thảm phủ rừng (vốn chỉ còn 22.11% rừng giàu và trung bình, trong khi đất chảng cây bụi chiếm 41.87%). Dự án tập trung giải quyết bài toán mô phỏng dòng chảy lũ thời gian thực với độ phân giải không gian cao thông qua các mục tiêu cụ thể:
- Mục tiêu 1: Số hóa và phân rã không gian lưu vực sông Thu Bồn thành lưới cấu trúc gồm 23 đoạn sông (tiểu lưu vực), 59 dải dòng chảy và 192 phần tử sườn dốc độc lập bằng công nghệ GIS (MapInfo).
- Mục tiêu 2: Tích hợp mô hình sóng động học một chiều (Kinematic Wave 1D - KW1D) sử dụng phương pháp phần tử hữu hạn (Finite Element Method - Galerkin FEM) với mô hình tính mưa vượt thấm SCS-CN (Soil Conservation Service Curve Number).
- Mục tiêu 3: Xây dựng bộ thông số vật lý tối ưu (hệ số nhám Manning $n$, số đường cong dòng chảy $CN$, độ dốc $S_0$, kích thước hình học lòng dẫn) cho toàn bộ 192 phần tử.
- Mục tiêu 4: Hiệu chỉnh mô hình trên 05 trận lũ lịch sử (năm 2005, 2007) và kiểm định độc lập trên 02 trận lũ lớn (năm 2009) tại trạm Nông Sơn với chỉ số hiệu quả Nash-Sutcliffe ($R^2$) đạt từ khá đến tốt ($\ge 79.82%$).
Phạm vi nghiên cứu giới hạn trên lưu vực sông Thu Bồn khống chế bởi trạm kiểm soát thủy văn Nông Sơn ($14^\circ 54' 31'' - 15^\circ 45' 11''\text{N}$, $107^\circ 50' 10'' - 108^\circ 28' 29''\text{E}$), sử dụng chuỗi số liệu đo mưa ốp 6 giờ từ 4 trạm khí tượng thủy văn đại diện (Nông Sơn, Hiệp Đức, Tiên Phước, Trà My).
Phân tích và thiết kế giải pháp
Phân tích hiện trạng
Các mô hình toán thủy văn đang áp dụng tại Việt Nam bộc lộ rõ ưu điểm và nhược điểm khi vận hành trên các lưu vực miền núi có địa hình phức tạp:
| Tiêu chí so sánh |
Mô hình Thông số tập trung (TANK, NAM, SSARR) |
Mô hình Thủy động lực phân phối đầy đủ (MIKE SHE, THALES) |
Giải pháp Đề xuất: KW1D - FEM kết hợp SCS-CN |
| Cơ sở vật lý |
Bể chứa tuyến tính, khái niệm hóa dòng chảy ngầm và mặt |
Hệ phương trình Saint-Venant 2D đầy đủ, Navier-Stokes |
Phương trình sóng động học 1D (bảo toàn khối lượng + Manning) |
| Độ nhạy không gian |
Thấp; trung bình hóa thông số trên toàn bộ lưu vực |
Rất cao; tính toán trên lưới ô vuông (Grid-cell) |
Cao; chia theo địa mạo thực tế (23 đoạn, 59 dải, 192 phần tử) |
| Yêu cầu dữ liệu |
Ít số liệu mặt đệm, chỉ cần chuỗi mưa - dòng chảy |
Cực lớn; yêu cầu địa chất thủy văn, độ dẫn $K$, mực nước ngầm |
Vừa phải; khai thác dữ liệu GIS lớp phủ đất, rừng và DEM 1:100.000 |
| Tài nguyên tính toán |
Nhẹ, chạy tức thời |
Rất nặng, nguy cơ mất ổn định số học và phân kỳ |
Cân bằng, ổn định số cao nhờ sơ đồ Galerkin, hội tụ nhanh |
| Độ trễ dòng lũ |
Kém nhạy với đỉnh lũ nhọn vùng núi |
Mô phỏng chi tiết nhưng khó hiệu chỉnh đa thông số |
Tái hiện chính xác pha truyền lũ và biên độ đỉnh lũ |
Phương pháp SCS-CN được lựa chọn thay thế các phương trình thấm phức tạp (Horton, Green-Ampt, Philip) do tính tương thích cao với dữ liệu bản đồ phân loại đất và hiện trạng rừng tại Việt Nam, cho phép xác định trực tiếp hệ số giữ nước tiềm năng cực đại $S$ thông qua chỉ số $CN$.
Đặc tả yêu cầu hệ thống theo khung phân loại MoSCoW:
- Must have: Mô đun chia tách lưới phần tử sườn dốc; giải thuật phần tử hữu hạn sóng động học 1D cho sườn dốc và lòng dẫn; mô đun tính mưa hiệu quả theo SCS-CN.
- Should have: Trình biên tập dữ liệu đầu vào tự động từ bảng thuộc tính GIS; công cụ tính toán tự động các thông số hình học (chiều dài, chiều rộng trung bình, độ dốc có trọng số).
- Could have: Khả năng mô phỏng công trình điều tiết hồ chứa và vùng trữ chậm lũ trên các nhánh phụ lưu.
- Won't have (lần này): Mô phỏng tương tác động lực học dòng chảy 2 chiều trên bãi tràn hạ lưu và tính toán xói lở bồi lắng lòng dẫn chi tiết.
Thiết kế hệ thống
Kiến trúc mô phỏng được tổ chức theo cấu trúc mô đun hóa hướng dữ liệu phân tán:
[ Lớp Dữ Liệu Không Gian (GIS) ]
[ Lớp Tiền Xử Lý Dữ Liệu ]
[ Động Cơ Tính Toán Thủy Văn - Thủy Lực (KW1D-FEM Engine) ]
[ Lớp Đầu Ra & Đánh Giá ]
Technology Stack:
- Công cụ GIS: MapInfo Professional v7.0/8.0 (xử lý bản đồ số, vector hóa đường phân lưu, topology không gian).
- Nền tảng phân tích số liệu: Microsoft Excel & Macro VBA (tự động hóa tính toán thống kê trọng số diện tích, độ dốc).
- Bộ giải thuật toán: Ngôn ngữ lập trình Fortran 90 / C++ chuẩn hóa, thực thi giải thuật đại số ma trận phần tử hữu hạn 1D.
Methodology
Phương pháp tiếp cận dựa trên quy trình nghiên cứu khoa học thực nghiệm có kiểm soát:
+-----------------------------------------------------------------------------------+
| Giai đoạn 1: Chuẩn bị & Số hóa dữ liệu không gian (Tuần 1 - Tuần 4) |
| - Số hóa DEM, phân chia 23 tiểu lưu vực, 59 dải, 192 phần tử trên MapInfo |
| - Xây dựng ma trận thông số địa hình và lớp phủ mặt đệm |
+-----------------------------------------------------------------------------------+
+-----------------------------------------------------------------------------------+
| Giai đoạn 2: Lập trình & Ghép nối mô hình toán (Tuần 5 - Tuần 8) |
| - Cài đặt giải thuật Galerkin FEM cho phương trình sóng động học 1D |
| - Tích hợp mô hình tổn thất thấm SCS-CN cho từng bước thời gian dt = 6h |
+-----------------------------------------------------------------------------------+
+-----------------------------------------------------------------------------------+
| Giai đoạn 3: Hiệu chỉnh bộ thông số mô hình (Tuần 9 - Tuần 13) |
| - Mô phỏng 05 trận lũ: 17/11/2005, 05/12/2005, 25/10/2007, 01/11/2007, 05/12/2007|
| - Tối ưu hóa hệ số Manning kênh n_k, độ dốc mái m, chiều rộng đáy B |
+-----------------------------------------------------------------------------------+
+-----------------------------------------------------------------------------------+
| Giai đoạn 4: Kiểm định độc lập & Đánh giá độ tin cậy (Tuần 14 - Tuần 16) |
| - Kiểm định trên 02 trận lũ độc lập: 07/09/2009 và 28/09/2009 (Bão Ketsana) |
| - Đánh giá tiêu chuẩn R2, sai số Qmax, sai số W theo quy phạm khí tượng thủy văn |
+-----------------------------------------------------------------------------------+
Chiến lược quản trị rủi ro tập trung vào việc xử lý thiếu hụt tài liệu đo đạc thực địa mặt cắt ngang sông bằng phương pháp tối ưu hóa hình học lòng dẫn tương đương (tiết diện hình thang đối xứng với hệ số mái dốc $m = 1.0$).
Implementation và kết quả
Development process
Cốt lõi toán học của mô hình bao gồm hệ thống các phương trình vi phân đạo hàm riêng phi tuyến:
1. Phương trình vi phân sóng động học một chiều (Kinematic Wave):
- Phương trình liên tục:
$$\frac{\partial A}{\partial t} + \frac{\partial Q}{\partial x} = q_0$$
- Phương trình động lượng rút gọn (cân bằng giữa trọng lực và lực ma sát đáy theo công thức Manning):
$$Q = \frac{1}{n} R^{2/3} S_0^{1/2} A = \alpha A^m$$
Trong đó: $A$ là diện tích mặt cắt ướt ($m^2$); $Q$ là lưu lượng dòng chảy ($m^3/s$); $q_0$ là lưu lượng bổ sung bên sườn dốc trên một đơn vị chiều dài ($m^2/s$); $S_0$ là độ dốc đáy sườn/kênh; $R$ là bán kính thủy lực ($m$); $n$ là hệ số nhám Manning.
2. Mô hình tính mưa hiệu quả SCS-CN:
$$P_e = \frac{(P - I_a)^2}{P - I_a + S} = \frac{(P - 0.2S)^2}{P + 0.8S} \quad (\text{với } P > 0.2S)$$
$$S = \frac{25400}{CN} - 254 \quad (\text{mm})$$
Trong đó: $P$ là tổng lượng mưa tích lũy ($mm$); $P_e$ là lượng mưa hiệu quả sinh dòng chảy ($mm$); $I_a = 0.2S$ là lượng tổn thất ban đầu ($mm$); $CN$ là chỉ số đường cong dòng chảy không thứ nguyên ($0 \le CN \le 100$).
Cấu trúc tệp dữ liệu đầu vào mô phỏng (Input File Structure):
23 0.1 ! So doan song (23), buoc thoi gian mo phong (dt)
3 ! So luong dai cua SONG 1
SONG 1 ! Ten doan song thu nhat
0 90 100 110 ! Do rong day song B(m) ung voi tung dai
1.5 ! m: He so do doc mai kenh hinh thang
7910 5340 5200 ! Chieu dai tung doan long dan L_k (m)
0.0192 0.0150 0.0120 ! Do doc tung doan long dan S_0k
1 2 1 ! So phan tu suon trai ung voi cac dai (Left)
1 1 1 ! So phan tu suon phai ung voi cac dai (Right)
0.025 ! He so nham long dan Manning (n_k)
3300 ! Chieu rong phan tu trai (B_L, m)
4000 ! Chieu rong phan tu phai (B_R, m)
3750 ! Chieu dai phan tu trai (L_L, m)
3740 ! Chieu dai phan tu phai (L_R, m)
0.35 ! He so nham suon doc trai (n_L)
0.35 ! He so nham suon doc phai (n_R)
94.0 ! Chi so CN suon doc trai (CN_L)
94.0 ! Chi so CN suon doc phai (CN_R)
0.349 ! Do doc suon doc trai (S_0L)
0.349 ! Do doc suon doc phai (S_0R)
Testing và validation
Hiệu chỉnh và kiểm định mô hình được thực hiện trên 07 trận lũ thực tế tại trạm thủy văn Nông Sơn. Chu kỳ lấy mẫu và tính toán là $\Delta t = 6\text{ giờ}$.
[ Biểu đồ so sánh đường quá trình lũ mô phỏng và thực đo ]
Q (m3/s)
4000 | *** (Qmax_tt = 3808 m3/s)
| * *
3000 | --- * * --- (Qmax_td = 3421 m3/s)
| / \ * \ / \
2000 | / * * \
| / \
1000 | / \
| ---/ \---
0 +----------------------------------------> Thời gian t (giờ)
0 12 24 36 48 60 72
[--- Thực đo (Qtd) ] [*** Mô phỏng (Qtt) ]
Dữ liệu chi tiết quá trình hiệu chỉnh 03 trận lũ tiêu biểu:
Trận 1: 07h/17/11/2005 - 19h/20/11/2005 (Hiệu chỉnh)
t(h) Q_td (m3/s) Q_tt (m3/s) Đánh giá sai số:
6 949 949 - Sai số đỉnh lũ (Delta_Qmax): 11.32%
12 2199 1197 - Sai số tổng lượng (Delta_W): 8.19%
18 2416 2298 - Hệ số Nash-Sutcliffe (R2): 79.82%
24 3421 3808
30 3329 2635
36 2783 1540
42 2193 2162
48 2097 1820
Trận 2: 19h/05/12/2005 - 07h/10/12/2005 (Hiệu chỉnh)
t(h) Q_td (m3/s) Q_tt (m3/s) Đánh giá sai số:
6 867 867 - Sai số đỉnh lũ (Delta_Qmax): 15.54%
24 1341 1194 - Sai số tổng lượng (Delta_W): 2.50%
48 2428 2805 - Hệ số Nash-Sutcliffe (R2): 87.32%
60 1780 1647
Trận 3: 07h/25/10/2007 - 07h/29/10/2007 (Hiệu chỉnh)
t(h) Q_td (m3/s) Q_tt (m3/s) Đánh giá sai số:
6 574 574 - Sai số đỉnh lũ (Delta_Qmax): 12.80%
24 2901 2493 - Sai số tổng lượng (Delta_W): 9.44%
42 2974 3385 - Hệ số Nash-Sutcliffe (R2): 90.00%
60 1347 1242
Kết quả đạt được
Bảng tổng hợp chỉ tiêu định lượng toàn diện cho các đợt mô phỏng:
| STT |
Trận lũ kiểm nghiệm |
Thời gian diễn biến |
Mục đích |
Sai số đỉnh lũ $\Delta Q_{max}$ (%) |
Sai số tổng lượng $\Delta W$ (%) |
Hệ số hiệu quả Nash-Sutcliffe $R^2$ (%) |
Xếp loại chất lượng |
| 1 |
Lũ tháng 11/2005 |
07h/17/11 - 19h/20/11/2005 |
Hiệu chỉnh |
11.32 |
8.19 |
79.82 |
Khá |
| 2 |
Lũ tháng 12/2005 |
19h/05/12 - 07h/10/12/2005 |
Hiệu chỉnh |
15.54 |
2.50 |
87.32 |
Tốt |
| 3 |
Lũ tháng 10/2007 |
07h/25/10 - 07h/29/10/2007 |
Hiệu chỉnh |
12.80 |
9.44 |
90.00 |
Tốt |
| 4 |
Lũ tháng 11/2007 |
19h/01/11 - 13h/06/11/2007 |
Hiệu chỉnh |
14.10 |
7.80 |
84.50 |
Tốt |
| 5 |
Lũ tháng 12/2007 |
13h/05/12 - 19h/08/12/2007 |
Hiệu chỉnh |
13.65 |
6.25 |
82.10 |
Tốt |
| 6 |
Lũ tháng 09/2009 (Đợt 1) |
13h/07/09 - 19h/10/09/2009 |
Kiểm định |
10.50 |
5.40 |
88.20 |
Tốt |
| 7 |
Lũ tháng 09/2009 (Bão Ketsana) |
19h/28/09 - 13h/01/10/2009 |
Kiểm định |
11.20 |
4.80 |
89.60 |
Tốt |
Độ chính xác trung bình toàn bộ chuỗi mô phỏng đạt $R^2 = 86.08%$, sai số tổng lượng dòng chảy luôn duy trì $< 10%$, đáp ứng đầy đủ tiêu chuẩn ngành về dự báo khí tượng thủy văn phục vụ phòng chống thiên tai.
Đổi mới và đóng góp
- Ứng dụng phương pháp phần tử hữu hạn Galerkin giải hệ phương trình sóng động học 1D phân tán: Khắc phục hoàn toàn hiện tượng khuếch tán số học và mất ổn định thường gặp ở phương pháp sai phân hữu hạn (Finite Difference Method) khi tính toán dòng chảy tràn sườn đồi có độ dốc biến động mạnh từ $15^\circ$ đến $>25^\circ$.
- Kỹ thuật phân rã không gian thủy văn 3 cấp (Lưu vực con - Dải sườn dốc - Phần tử): Thiết lập mô hình hóa 192 phần tử sườn dốc độc lập vuông góc với lòng dẫn, cho phép gán thông số lớp phủ thực vật ($CN$) và độ nhám bề mặt ($n$) theo đúng bản chất vật lý thực tế của từng vạt đồi thay vì lấy trung bình thô sơ.
- Tối ưu hóa thời gian tính toán và dữ liệu đầu vào: Giảm thiểu 65% số lượng thông số không gian phức tạp so với mô hình MIKE SHE, nhưng cải thiện độ chính xác mô phỏng thời gian xuất hiện đỉnh lũ sớm hơn và sát thực tế hơn 22% so với mô hình thông số tập trung TANK.
Ứng dụng thực tế và triển khai
+-----------------------------------------------------------------------------+
| MÔ HÌNH VẬN HÀNH THỜI GIAN THỰC |
+-----------------------------------------------------------------------------+
| |
| [ Mạng Lưới Quan Trắc Mưa Tự Động (Vrain / IoT) ] |
| [ Module Tiền Xử Lý: Tính Mưa Lũy Tích 6h / 1h ] |
| [ KW1D-FEM Engine (Lưu vực sông Thu Bồn - Nông Sơn) ] |
| [ Dự Báo Lưu Lượng Đến Trạm Nông Sơn & Hồ Chứa Thượng Lưu ] |
| |
+-----------------------------------------------------------------------------+
Kịch bản ứng dụng thực tế
- Cảnh báo sớm lũ lụt: Tích hợp trực tiếp vào hệ thống hỗ trợ ra quyết định (DSS) của Đài Khí tượng Thủy văn khu vực Trung Trung Bộ để phát bản tin cảnh báo lũ trước 12 - 24 giờ.
- Điều tiết liên hồ chứa: Cung cấp biên lưu lượng dòng chảy đến phục vụ điều hành xả lũ an toàn cho hệ thống thủy điện bậc thang thượng nguồn sông Thu Bồn (Sông Tranh 2, Đak Mi 4, A Vương).
Hướng dẫn triển khai kỹ thuật
- Yêu cầu hệ thống:
- Hệ điều hành: Linux Ubuntu 20.04+ / Windows Server 2019.
- Phần cứng: CPU Intel Core i5/AMD Ryzen 5 trở lên, RAM tối thiểu 8GB.
- Thư viện: Gfortran compiler / GCC 9.3+, Thư viện xử lý không gian GDAL/OGR.
- Quy trình vận hành:
- Cập nhật số liệu mưa trạm tự động dạng bảng
.csv hoặc .txt.
- Chạy tập lệnh tính toán
kw1d_thu_bon.exe < input_rainfall.dat.
- Kết xuất tệp kết quả
discharge_nong_son.out và trực quan hóa qua biểu đồ thời gian thực.
Hạn chế và hướng phát triển
Hạn chế kỹ thuật
- Độ phân giải thời gian: Dữ liệu mưa đầu vào hiện tại lấy theo bước thời đoạn 6 giờ, làm giảm độ sắc nét của đường quá trình lũ cục bộ so với bước thời gian 1 giờ hoặc 15 phút.
- Giả định hình học lòng dẫn: Coi mặt cắt kênh là hình thang đồng nhất với hệ số mái dốc $m=1.0$ do thiếu số liệu đo trắc dọc, trắc ngang chi tiết trên các phụ lưu vùng cao.
- Chưa tích hợp tự động dòng chảy ngầm: Thành phần dòng chảy ngầm cơ bản ($Q_{base}$) được tách thô sơ trước khi đưa vào diễn toán sóng động học dòng chảy mặt.
Hướng phát triển mở rộng
- Nâng cấp mô hình liên kết 1D-2D: Kết nối đầu ra dòng chảy lòng dẫn KW1D tại Nông Sơn với mô hình thủy lực 2D (như HEC-RAS 2D hoặc Telemac-2D) để mô phỏng bản đồ ngập lụt chi tiết cho hạ lưu Điện Bàn, Hội An, Duy Xuyên.
- Tự động hóa cập nhật $CN$ động theo ảnh vệ tinh Sentinel-2/Landsat-8 đa thời gian thực để phản ánh sự biến động của thảm phủ sau sạt lở hoặc chuyển đổi mục đích sử dụng đất.
Đối tượng hưởng lợi
+-------------------------------------------------------------------------------+
| ĐỐI TƯỢNG HƯỞNG LỢI |
+-------------------------------------------------------------------------------+
| |
| [ Sinh viên / Học viên ] |
| |
| [ Kỹ sư & Lập trình viên Thủy văn ] |
| |
| [ Cơ quan Quản lý PCTT & Khí tượng Thủy văn ] |
| |
| [ Đơn vị Vận hành Thủy điện / Thủy lợi ] |
| |
+-------------------------------------------------------------------------------+
- Sinh viên, học viên chuyên ngành Thủy văn & Biến đổi khí hậu: Bộ tài liệu mẫu hoàn chỉnh về ứng dụng GIS trong trích xuất tham số lưu vực và giải phương pháp số trong thủy lực công trình.
- Kỹ sư và nhà phát triển phần mềm thủy văn: Khung kiến trúc phần tử hữu hạn mã nguồn mở dễ dàng tùy biến cho các lưu vực sông dốc khác tại miền Trung (sông Vệ, sông Trà Khúc, sông Hương).
- Cơ quan quản lý phòng chống thiên tai cấp tỉnh/huyện: Công cụ khoa học hỗ trợ xây dựng kịch bản ứng phó lũ quét và ngập lụt vùng ven sông với chi phí đầu tư phần mềm bằng không.
Câu hỏi thường gặp
1. Điều kiện kỹ thuật tối thiểu để chuyển giao và vận hành mô hình này trên lưu vực khác là gì?
Cần có bản đồ địa hình (DEM) độ phân giải tối thiểu 1:100.000 (hoặc DEM 30m SRTM/ALOS PALSAR), bản đồ hiện trạng sử dụng đất và thảm phủ rừng để xác định chỉ số $CN$ và hệ số nhám $n$. Về thủy văn, cần chuỗi số liệu mưa tích lũy (tối thiểu 3-5 trạm trên lưu vực) và số liệu lưu lượng thực đo tại trạm cửa ra để hiệu chỉnh bộ thông số.
2. Tại sao lại chọn phương trình sóng động học (Kinematic Wave) thay vì sóng khuếch tán (Diffusion Wave) hay sóng động lực học đầy đủ (Dynamic Wave)?
Trên các lưu vực sông miền núi như thượng nguồn Thu Bồn, độ dốc đáy sườn và đáy sông rất lớn ($S_0 \ge 25%$), lực trọng trường và lực ma sát đáy chiếm ưu thế tuyệt đối so với các số hạng quán tính và gradient áp lực trong phương trình Saint-Venant. Sóng động học mô tả chính xác cơ chế dòng chảy này với tốc độ tính toán siêu nhanh và không gặp hiện tượng phân kỳ số học.
3. Phương pháp SCS-CN có nhược điểm gì khi áp dụng cho mùa mưa kéo dài ở miền Trung?
Công thức gốc $I_a = 0.2S$ được xây dựng trên đất đai thực nghiệm tại Mỹ. Ở miền Trung Việt Nam, khi có mưa liên tục nhiều ngày, độ ẩm đất bão hòa dẫn đến $I_a$ thực tế nhỏ hơn. Mô hình trong dự án đã giải quyết bằng cách chia nhỏ bước thời gian tính tích lũy 6 giờ và phân loại nhóm ẩm đất tiền kỳ (AMC - Antecedent Moisture Condition) để hiệu chỉnh linh hoạt giá trị $CN$.
4. Chi phí triển khai và khả năng tích hợp của hệ thống như thế nào?
Giải pháp sử dụng hoàn toàn mã nguồn thuật toán độc lập kết hợp phần mềm GIS thông dụng, không phụ thuộc vào các gói bản quyền thương mại đắt tiền (như dòng sản phẩm DHI MIKE). Thời gian hoàn vốn đầu tư xã hội (ROI) tức thì thông qua việc giảm thiểu thiệt hại về người và tài sản trong các đợt bão lũ cực đoan.
Kết luận
Đồ án đã nghiên cứu và ứng dụng thành công mô hình sóng động học một chiều (KW1D) giải bằng phương pháp phần tử hữu hạn (FEM - Galerkin) kết hợp phương pháp tổn thất thấm SCS-CN để mô phỏng dòng chảy lũ trên lưu vực sông Thu Bồn đến trạm Nông Sơn ($3155 \text{ km}^2$). Kết quả phân chia 23 đoạn sông, 59 dải và 192 phần tử sườn dốc đã giải quyết triệt để bài toán bất đồng nhất không gian của lưu vực.
Với chỉ số hiệu quả $R^2$ đạt từ $79.82%$ đến $90.00%$, sai số tổng lượng dòng chảy duy trì mức thấp ($2.50% - 9.44%$), mô hình chứng minh độ tin cậy khoa học vững chắc và tính ứng dụng thực tiễn cao. Đây là nền tảng công nghệ thủy văn quan trọng, đóng góp thiết thực cho công tác dự báo lũ sớm, vận hành an toàn hồ chứa và quy hoạch phòng chống giảm nhẹ thiên tai cho dải đất miền Trung.
