Tổng quan về luận án

Biến đổi khí hậu toàn cầu và sự gia tăng của các hiện tượng thời tiết cực đoan đang đặt các dải đồng bằng ven biển miền Trung Việt Nam trước những thách thức sống còn về an ninh tài nguyên nước và quản lý rủi ro thiên tai. Luận án tiến sĩ chuyên ngành Thủy văn học và Kỹ thuật Tài nguyên nước với đề tài "Nghiên cứu tính toán thiệt hại do lũ vùng hạ lưu sông Vu Gia - Thu Bồn" giải quyết trực tiếp khoảng trống học thuật (research gap) cốt lõi: sự thiếu hụt một khung mô hình toán tích hợp liên kết động giữa cơ chế thủy văn - thủy lực dòng chảy tràn một chiều với hệ thống thông tin địa lý (GIS) và hàm lượng giá thiệt hại kinh tế thực chứng tại một lưu vực có địa hình chia cắt phức tạp, độ dốc lớn ($25%$) và mạng lưới phân lưu chuyển dòng phi tuyến tính.

       +-------------------------------------------------------------+
       |                  MÔ HÌNH THỦY VĂN MIKE NAM                  |
       |  Mô phỏng mưa - dòng chảy 4 bể trữ (Nielsen & Hansen, 1973) |
       +------------------------------+------------------------------+
                                      | Biên lưu lượng nhập lưu Q(t)
                                      v
       +-------------------------------------------------------------+
       |                  MÔ HÌNH THỦY LỰC MIKE 11 HD                |
       |     Phương trình Saint-Venant 1D (Saint-Venant, 1871)       |
       |   Mạng lưới sông: Phân lưu Quảng Huế <-> Sông Vĩnh Điện     |
       +------------------------------+------------------------------+
                                      | Chuỗi mực nước H(t) & Diện ngập
                                      v
       +-------------------------------------------------------------+
       |                  KHUNG TÍCH HỢP MIKE 11 GIS                 |
       |    Chồng xếp bản đồ số độ cao DEM & Không gian hóa ngập lũ   |
       +------------------------------+------------------------------+
                                      | Độ sâu ngập & Bản đồ ngập p%
                                      v
       +-------------------------------------------------------------+
       |                 LƯỢNG GIÁ THIỆT HẠI KINH TẾ                 |
       |    Đường cong Hmax - Thiệt hại quy đổi (Khung ECLAC, 2007)  |
       |      Áp dụng cho Huyện Điện Bàn (Nhà cửa, Lúa, Hạ tầng)     |
       +-------------------------------------------------------------+

Hệ thống câu hỏi nghiên cứu và giả thuyết khoa học được xác lập chặt chẽ:

  • Câu hỏi nghiên cứu 1 (RQ1): Cơ chế thủy động lực của dòng lũ biến đổi như thế nào khi truyền qua mạng lưới sông phân lưu phức tạp Vu Gia - Thu Bồn dưới các kịch bản tần suất lũ thiết kế $p = 1%$, $2%$, $5%$, và $10%$?
  • Câu hỏi nghiên cứu 2 (RQ2): Mức độ tương tác giữa dòng chảy lũ thượng nguồn và khả năng thoát lũ vùng cửa sông tạo ra diện tích và độ sâu ngập lụt thực tế ra sao tại hạ du?
  • Câu hỏi nghiên cứu 3 (RQ3): Mối quan hệ định lượng giữa các đặc trưng ngập lụt (mực nước đỉnh lũ $H_{\max}$, diện tích ngập) và tổn thất kinh tế đa ngành (nhà cửa, nông nghiệp, hạ tầng giao thông) được mô hình hóa như thế nào?
  • Giả thuyết khoa học 1 (H1): Việc ghép nối chuỗi mô hình MIKE NAM - MIKE 11 HD với công cụ phân tích không gian MIKE 11 GIS kiểm soát được sai số truyền lũ trên mạng sông có phân lưu nhân tạo/tự nhiên (sông Quảng Huế và sông Vĩnh Điện), cho phép tái hiện chính xác đường quá trình ngập lụt lịch sử với chỉ số hiệu quả Nash-Sutcliffe ($NSE$) và hệ số tương quan đạt độ tin cậy cao.
  • Giả thuyết khoa học 2 (H2): Thiệt hại kinh tế do lũ tại vùng hạ lưu không tuyến tính đơn thuần với lượng mưa mà phụ thuộc phi tuyến vào mực nước đỉnh lũ $H_{\max}$ thông qua đường cong tần suất - thiệt hại thực chứng được quy đổi theo giá trị hiện tại ròng.

Nghiên cứu được triển khai trên toàn bộ lưu vực sông Vu Gia - Thu Bồn với diện tích tự nhiên $10.350\text{ km}^2$ (chiếm $87,5%$ diện tích tỉnh Quảng Nam và thành phố Đà Nẵng), tập trung lượng giá chi tiết tại huyện Điện Bàn – trọng điểm ngập lụt hạ du. Chuỗi số liệu khí tượng thủy văn được khai thác liên tục trong giai đoạn 1976–2010 từ 24 trạm quan trắc (trong đó có 2 trạm đo lưu lượng Thành Mỹ, Nông Sơn và 7 trạm đo mực nước), kết hợp kiểm định qua các trận lũ lịch sử 1998, 1999, 2002, 2004 và 2009. Nghiên cứu mang tính đột phá khi định lượng hóa mức độ rủi ro: năm 2009, bão lũ đã cướp đi sinh mạng của 52 người, làm 220 người bị thương và cuốn trôi $3.500$ tỷ đồng, tương đương gần như toàn bộ tổng sản phẩm trên địa bàn (GDP 6 tháng đầu năm 2009 của Quảng Nam đạt $4.140$ tỷ đồng).

Literature Review và Positioning

Các nghiên cứu thủy văn - thủy lực và lượng giá tổn thất lũ lụt trên thế giới đã trải qua quá trình tiến hóa mạnh mẽ. Trường phái mô phỏng thủy văn truyền thống dựa trên các mô hình tất định thông số tập trung như mô hình Tank (Sugawara, 1956), mô hình giáng thủy - dòng chảy NAM (Nielsen & Hansen, 1973; Abbott & Refsgaard, 1996; Refsgaard & Knudsen, 1997) hay mô hình bán phân bố HEC-HMS của Trung tâm Kỹ thuật Thủy văn Quân đội Hoa Kỳ (USACE).

Về mặt thủy lực, các nghiên cứu quốc tế đã chứng minh tính hiệu quả của phương trình Saint-Venant 1D trong dự báo thời gian thực và quản lý lưu vực lớn. Điển hình, Daniel Jilles & Matthew Moore (2010) khi nghiên cứu hệ thống dự báo lũ tại Hà Lan, Bỉ và Anh Quốc đã khẳng định mô hình thủy lực 1D (MIKE 11, HEC-RAS) là công cụ tối ưu nhất về độ ổn định số học và tài nguyên tính toán cho công tác cảnh báo khẩn cấp quy mô quốc gia. Song song đó, Knebla cùng các cộng sự (2005) đã tích hợp HEC-HMS, HEC-RAS với dữ liệu radar NEXRAD và ArcHydro GIS cho lưu vực sông San Antonio ($10.000\text{ km}^2$, Texas, Hoa Kỳ), mở ra hướng tiếp cận "Map to Map". Tại Bỉ, Vanderkimpen đã ứng dụng mô hình kết nối 1D-2D (MIKE FLOOD) để xác lập bản đồ ngập lụt ven biển phục vụ sơ tán dân cư. Viện Nghiên cứu Công chính Nhật Bản (PWRI, 2007) phát triển hệ thống IFAS phân tích dòng chảy phân bố tích hợp mưa vệ tinh trên nền tảng DEM.

Công trình nghiên cứu Phạm vi & Phương pháp Ưu điểm cốt lõi Hạn chế / Khoảng trống học thuật
Knebla et al. (2005) Sông San Antonio, Texas ($10.000\text{ km}^2$); HEC-HMS + HEC-RAS + NEXRAD Radar + ArcHydro GIS. Tự động hóa chuỗi xử lý dữ liệu không gian "Map to Map"; mô phỏng nhanh lưu vực rộng. Chưa xây dựng hàm thiệt hại kinh tế đa ngành chi tiết theo từng cấp độ sâu ngập thực tế.
Vanderkimpen (2010) Đồng bằng ven biển Bỉ; Mô hình MIKE FLOOD (khớp nối thủy lực 1D-2D). Mô phỏng chi tiết trường vận tốc 2D vùng bãi tràn ven biển phục vụ di tản khẩn cấp. Chi phí tính toán cao; đòi hỏi số liệu địa hình DEM độ phân giải siêu cao, khó áp dụng cho lưu vực nghèo số liệu.
Phạm Ngọc (2000) Đồng bằng sông Cửu Long; Mô hình VRSAP kết hợp điều tra thiệt hại sau lũ năm 2000. Xây dựng đường cong tần suất - thiệt hại cho 4 loại hình sử dụng đất (dân cư, thương nghiệp, công nghiệp, hạ tầng). Mô hình VRSAP mã nguồn đóng khó tùy biến liên kết GIS động; địa hình ĐBSCL bằng phẳng, không có tính chất lũ quét như miền Trung.
Nguyễn Văn Tuấn et al. Hạ lưu sông Ba; Ứng dụng mô hình HEC-FDA lượng hóa giá trị thiệt hại kinh tế. Kết hợp kết quả thủy lực và phân tích kinh tế rủi ro đa kịch bản. Chưa tích hợp mô hình phân lưu phức tạp hai chiều giữa hai dòng sông chính như Vu Gia - Thu Bồn.
Luận án này Toàn bộ lưu vực Vu Gia - Thu Bồn ($10.350\text{ km}^2$); MIKE NAM + MIKE 11 HD + MIKE 11 GIS + Khung ECLAC (2007). Giải quyết trọn vẹn bài toán truyền lũ phân lưu kép (Quảng Huế - Vĩnh Điện) và hàm thiệt hại quy đổi $H_{\max}$ thực chứng. Cung cấp công cụ mô phỏng ngập lụt và lượng giá kinh tế trực tiếp phục vụ quy hoạch phòng chống thiên tai cấp tỉnh.

Trong nước, các nghiên cứu của Đinh Đức Trường & Lê Hà Thanh (2013) tại Thừa Thiên Huế đã tiếp cận khung lượng giá môi trường và kinh tế sau thiên tai của Ủy ban Kinh tế Châu Mỹ Latinh và Caribe (ECLAC, 2007; UNDAC, 2006). Tuy nhiên, việc gắn kết hữu cơ giữa mô hình thủy lực thủy văn với đường cong thiệt hại thực tế trên lưu vực sông Vu Gia - Thu Bồn vẫn tồn tại các tranh biện học thuật:

  • Quan điểm thứ nhất: Cho rằng chỉ cần sử dụng mô hình thủy lực 1D truyền thống kết hợp điều tra xã hội học hồi cố là đủ để quy hoạch đê điều.
  • Quan điểm đối lập: Khẳng định tính chất lũ miền Trung "lên nhanh, rút chậm" do lòng sông dốc ($2‰$ từ Thành Mỹ đến Ái Nghĩa, nhưng giảm đột ngột xuống $0,08‰$ từ Ái Nghĩa đến Câu Lâu và $0,04‰$ từ Câu Lâu ra biển) cùng tác động nghẽn dòng tại các cửa biển đòi hỏi phải mô hình hóa không gian ngập lụt trên nền GIS và tính toán tổn thất quy đổi tài chính theo từng cấp tần suất.

Luận án định vị chính xác vào giao điểm này: vừa giải quyết bài toán thủy lực mạng lưới sông phân lưu phức tạp (Vu Gia chuyển dòng sang Thu Bồn qua sông Quảng Huế, sau đó $16\text{ km}$ hạ lưu, sông Vĩnh Điện lại dẫn nước ngược từ Thu Bồn về Vu Gia), vừa xây dựng bộ hàm quan hệ thiệt hại kinh tế $H_{\max}$ - Thiệt hại thực chứng cho từng đối tượng dễ bị tổn thương.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng và kiểm chứng lý thuyết thủy động lực học Saint-Venant (1871) và lý thuyết cân bằng ẩm đa tầng trong chu kỳ thủy văn của Nielsen & Hansen (1973) trong điều kiện khí hậu nhiệt đới gió mùa khắc nghiệt sườn Đông dãy Trường Sơn:

$$\frac{\partial Q}{\partial x} + \frac{\partial A}{\partial t} = q$$

$$\frac{\partial Q}{\partial t} + \frac{\partial}{\partial x}\left(\frac{\alpha Q^2}{A}\right) + gA\frac{\partial h}{\partial x} + gA S_f = 0$$

Trong đó:

  • $Q$: Lưu lượng dòng chảy mặt ($\text{m}^3/\text{s}$)
  • $A$: Diện tích mặt cắt ướt lòng dẫn ($\text{m}^2$)
  • $q$: Lưu lượng gia nhập dọc bờ ($\text{m}^2/\text{s}$)
  • $h$: Mực nước dòng sông ($\text{m}$)
  • $S_f$: Độ dốc ma sát thủy lực, được xác định qua hệ số nhám Manning $n$: $S_f = \frac{n^2 Q |Q|}{A^2 R^{4/3}}$
  • $\alpha$: Hệ số hiệu chỉnh động lượng Boussinesq; $g$: Gia tốc trọng trường ($9,81\text{ m/s}^2$)
+-------------------------------------------------------------------------+
|                  CẤU TRÚC 4 TẦNG TRỮ ẨM MÔ HÌNH MIKE NAM                |
+-------------------------------------------------------------------------+
| [ Lượng mưa P(t) ]                                                      |
|        |                                                                |
|        v                                                                |
| +------------------------------------+                                  |
| | 1. Bể trữ bề mặt (Surface Storage) | ---> [ Dòng chảy mặt QOF ]       |
| |    Dung tích cực đại: Umax         |                                  |
| +-----------------+------------------+                                  |
|                   | Thấm (Infiltration)                                 |
|                   v                                                     |
| +------------------------------------+                                  |
| | 2. Bể trữ tầng rễ (Root Zone)      | ---> [ Bốc thoát hơi ẩm Ea ]     |
| |    Dung tích cực đại: Lmax         |                                  |
| +-----------------+------------------+                                  |
|                   | Thấm sâu (Percolation)                              |
|                   v                                                     |
| +------------------------------------+                                  |
| | 3. Bể trữ ngầm nông (Upper Base)   | ---> [ Dòng chảy sát mặt QIF ]   |
| +-----------------+------------------+                                  |
|                   | Chuyển tiếp                                         |
|                   v                                                     |
| +------------------------------------+                                  |
| | 4. Bể trữ ngầm sâu (Lower Base)    | ---> [ Dòng chảy kiệt QBF ]      |
| +------------------------------------+                                  |
+-------------------------------------------------------------------------+

Về mặt lý thuyết thủy văn, công trình chứng minh cấu trúc 4 tầng trữ ẩm liên tục: bể mặt ($U_{\max}$), bể rễ cây ($L_{\max}$), bể ngầm trên và bể ngầm dưới của mô hình NAM phản ánh xuất sắc trạng thái bão hòa độ ẩm đất của miền Trung trước các đợt mưa lớn tháng 10–11, giải thích nguyên nhân các trận lũ "kép" và "đa đỉnh" có biên độ lũ lên tới $5 - 14\text{ m}$ và cường suất lũ đạt đỉnh $100 - 140\text{ cm/h}$.

Khung phân tích độc đáo

Khung phân tích của luận án là sự hợp nhất của 3 trụ cột lý thuyết:

  1. Lý thuyết lan truyền sóng lũ thủy lực 1 chiều trên mạng sông hở có bãi ngập và ô chứa điều tiết.
  2. Lý thuyết phân tích không gian địa lý trên nền tảng Raster/TIN của GIS để chuyển đổi kết quả mực nước 1D thành trường ngập 2D tương hỗ.
  3. Khung lý thuyết lượng giá thiệt hại kinh tế sau thiên tai của ECLAC (2007), lượng hóa mối quan hệ trực tiếp giữa độ sâu ngập $H$, thời gian ngâm nước $T$, và tỷ lệ suy giảm giá trị tài sản đối với công trình xây dựng, sản xuất nông nghiệp (cây lúa) và mạng lưới hạ tầng kỹ thuật.

Điều kiện biên và phạm vi ứng dụng (boundary conditions) được xác định rõ ràng: mô hình vận hành tối ưu cho các lưu vực sông dốc khu vực nhiệt đới gió mùa có mạng lưới phân lưu và cửa sông chịu ảnh hưởng trực tiếp của chế độ bán nhật triều/nhật triều biển Đông.

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Nghiên cứu đứng trên lập trường bản thể luận hiện thực phản biện (critical realism) và phương pháp luận thực chứng (positivism), kết hợp chặt chẽ giữa mô hình hóa toán học tất định, công nghệ không gian địa lý và phương pháp phân tích thống kê kinh tế lượng.

+-------------------------------------------------------------------------+
|                  THIẾT KẾ NGHIÊN CỨU VÀ DÒNG DỮ LIỆU                    |
+-------------------------------------------------------------------------+
| [ DỮ LIỆU ĐẦU VÀO ]                                                     |
| - Khí tượng: 13 trạm mưa ngày, 3 trạm khí tượng (1976-2010)            |
| - Địa hình: DEM, mặt cắt ngang sông, trắc dọc, ô chứa điều tiết         |
| - Thủy văn: Q(t) tại Nông Sơn, Thành Mỹ; H(t) tại Sơn Trà, Ái Nghĩa     |
|                               |                                         |
|                               v                                         |
| [ HIỆU CHỈNH & KIỂM ĐỊNH MÔ HÌNH ]                                      |
| - Thủy lực MIKE 11: Nhám Manning n = 0,025 - 0,045                      |
| - Bộ chuỗi lũ kiểm định: 1998, 1999 (đại hồng thủy), 2002, 2004         |
|                               |                                         |
|                               v                                         |
| [ KỊCH BẢN THỦY LỰC THEO TẦN SUẤT THIẾT KẾ ]                            |
| - Mô phỏng lũ thiết kế: p = 1% (1999), p = 2%, p = 5%, p = 10%          |
| - Chuyển giao mực nước sang MIKE 11 GIS & ArcGIS 9.1                    |
|                               |                                         |
|                               v                                         |
| [ XUẤT BẢN ĐỒ NGẬP & ĐÁNH GIÁ THIỆT HẠI KINH TẾ ]                       |
| - Bản đồ diện tích, cấp độ sâu ngập huyện Điện Bàn                     |
| - Hàm thiệt hại quy đổi năm 2014 (Nhà cửa, Đất lúa, Giao thông)         |
+-------------------------------------------------------------------------+

Hệ thống mẫu và mạng lưới tính toán bao gồm:

  • Toàn bộ mạng lưới sông chính Vu Gia - Thu Bồn, sông Vĩnh Điện, sông Quảng Huế, sông Cẩm Lệ, sông Hàn, sông Duy Vinh, sông Thu Cồn.
  • 13 trạm đo mưa đại biểu (Thành Mỹ, Nông Sơn, Tam Kỳ, Ái Nghĩa, Câu Lâu, Đà Nẵng, Cẩm Lệ, Hội Khách, Hội An, Khâm Đức, Tiên Phước, Trà My, Trao).
  • 2 trạm đo lưu lượng thượng lưu ($F_{\text{khống chế Thành Mỹ}} = 1.850\text{ km}^2$, $F_{\text{khống chế Nông Sơn}} = 3.155\text{ km}^2$).
  • Biên dưới là dao động mực nước triều thực đo tại trạm hải văn Sơn Trà.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu tuân thủ nghiêm ngặt chuẩn mực khoa học thủy văn quốc tế:

  1. Thu thập và chuẩn hóa dữ liệu: Đồng nhất hóa chuỗi quan trắc khí tượng thủy văn 35 năm (1976–2010). Sử dụng phương pháp nội suy không gian và đường đẳng diện tích để phân bổ lượng mưa từ tâm mưa Trà My ($4.000\text{ mm}$) về đồng bằng ($2.000 - 2.200\text{ mm}$).
  2. Thiết lập sơ đồ thủy lực và ô chứa: Xây dựng hệ thống các mặt cắt ngang lòng sông, đê bao, cống điều tiết và các ô chứa lũ vùng đồng bằng hạ lưu trong mô hình MIKE 11 để mô phỏng chính xác hiện tượng tràn bờ vào ruộng đồng khi mực nước vượt bờ tại Ái Nghĩa và Giao Thủy.
  3. Hiệu chỉnh và kiểm định (Calibration & Validation): Hiệu chỉnh bộ thông số mô hình với trận lũ năm 1998, 1999; kiểm định độc lập với trận lũ các năm 2002 và 2004. Đánh giá sự tương thích giữa đường quá trình mực nước thực đo và tính toán tại các trạm Ái Nghĩa, Giao Thủy.
+-------------------------------------------------------------------------+
|           BẢNG THÔNG SỐ VÀ KẾT QUẢ ĐO ĐẠC THỦY VĂN LƯU VỰC              |
+---------------------+-------------------+-------------------------------+
| Đại lượng quan trắc | Trạm Thành Mỹ     | Trạm Nông Sơn                 |
+---------------------+-------------------+-------------------------------+
| Diện tích lưu vực   | 1.850 km²         | 3.155 km²                     |
| Qmax lịch sử        | 7.000 m³/s (1998) | 10.600 m³/s (2007)            |
| Vận tốc lũ max      | 3,77 m/s          | 3,74 m/s                      |
| Biên độ lũ max      | 15,2 m            | 12,0 m                        |
| Thời gian truyền lũ | Đến Ái Nghĩa:     | Đến Giao Thủy:                |
|                     | 5 - 11 giờ        | 3 - 7 giờ                     |
+---------------------+-------------------+-------------------------------+

Data và phân tích

Phân tích số liệu áp dụng hệ thống công cụ tiên tiến:

  • MIKE 11 HD & MIKE NAM: Giải hệ phương trình phi tuyến Saint-Venant bằng sơ đồ sai phân hữu hạn ẩn Abbott-Ionescu 6 điểm, đảm bảo độ hội tụ và ổn định số học cao.
  • MIKE 11 GIS & ArcGIS 9.1: Xử lý mô hình số độ cao DEM, nội suy mặt nước ngập lũ 3D từ các mặt cắt ngang và trừ lưới địa hình để trích xuất bản đồ độ sâu ngập cho 4 kịch bản tần suất $p = 1%$, $2%$, $5%$, $10%$.
  • Lượng giá kinh tế: Sử dụng phương pháp điều tra thống kê thiệt hại giai đoạn 1992–2010, đưa tất cả các giá trị thiệt hại lịch sử về mặt bằng giá cố định năm 2014 thông qua chỉ số giảm phát GDP và CPI, từ đó thiết lập đường tương quan phi tuyến giữa $H_{\max}$ tại Nông Sơn, Ái Nghĩa với tổng thiệt hại kinh tế.

Phát hiện đột phá và implications

Những phát hiện then chốt

  1. Đặc trưng động lực học lũ núi cực đoan: Lưu vực sông Vu Gia - Thu Bồn có cường suất lũ đặc biệt cao ($20 - 50\text{ cm/h}$, lớn nhất lên tới $100 - 140\text{ cm/h}$ tại thượng lưu). Vận tốc dòng chảy lũ lớn nhất tại trạm Thành Mỹ đạt $3,77\text{ m/s}$ (biên độ lũ cực đại $15,2\text{ m}$), tại Nông Sơn đạt $3,74\text{ m/s}$ (biên độ lũ $12,0\text{ m}$). Lưu lượng đỉnh lũ thực đo đạt $7.000\text{ m}^3/\text{s}$ tại Thành Mỹ (20/11/1998) và $10.600\text{ m}^3/\text{s}$ tại Nông Sơn (12/11/2007).
  2. Quy luật lan truyền lũ phi đối xứng: Thời gian truyền lũ từ Thành Mỹ về Ái Nghĩa ($40,5\text{ km}$) dao động từ $5 - 11\text{ giờ}$; từ Nông Sơn về Giao Thủy ($26\text{ km}$) từ $3 - 7\text{ giờ}$; từ Giao Thủy về Câu Lâu ($23\text{ km}$) trung bình $7,4\text{ giờ}$. Độ dốc lòng sông hạ lưu suy giảm đột ngột ($0,08‰$ và $0,04‰$) kết hợp tác động của thủy triều làm lũ "lên nhanh nhưng rút cực kỳ chậm", gây úng ngập kéo dài nhiều ngày tại vùng đồng bằng hạ du.
  3. Phân hóa diện tích ngập lụt theo tần suất thiết kế: Bản đồ ngập lụt xây dựng từ MIKE 11 GIS chỉ ra sự mở rộng phi tuyến của diện ngập trên toàn lưu vực và tại huyện Điện Bàn tương ứng với 4 tần suất thiết kế $p = 1%$, $2%$, $5%$, và $10%$.
Tần suất lũ thiết kế ($p%$) Chu kỳ lặp lại (Năm) Mực nước đỉnh lũ $H_{\max}$ trạm Ái Nghĩa (m) Diện tích ngập lụt huyện Điện Bàn (ha) Tổng thiệt hại kinh tế ước tính huyện Điện Bàn (Tỷ đồng)
$1%$ 100 năm 10,77 14.850 485,6
$2%$ 50 năm 10,12 13.210 362,4
$5%$ 20 năm 9,35 11.450 248,1
$10%$ 10 năm 8,60 9.120 145,8
  1. Hàm quan hệ định lượng giữa $H_{\max}$ và thiệt hại kinh tế quy đổi: Tổng hợp chuỗi số liệu thiệt hại giai đoạn 1992–2010 (quy đổi về giá năm 2014) cho thấy mối tương quan chặt chẽ giữa mực nước đỉnh lũ trạm Ái Nghĩa, Nông Sơn với tổn thất thực tế. Trận lũ lịch sử năm 1999 gây ngập hàng nghìn ngôi nhà cấp 4 và bán kiên cố, làm tê liệt $100%$ các tuyến giao thông liên xã và phá hủy hàng nghìn hecta lúa mùa tại Điện Bàn.
  2. Cơ chế nghẽn dòng và chuyển lưu Quảng Huế - Vĩnh Điện: Nghiên cứu phát hiện sông Quảng Huế đóng vai trò xung yếu trong việc san sẻ lưu lượng từ Vu Gia sang Thu Bồn trong các pha lũ sớm, nhưng khi mực nước sông Thu Bồn dâng cao vượt ngưỡng kiểm soát, dòng chảy qua sông Vĩnh Điện bị nghẽn và chảy ngược, làm trầm trọng thêm mức độ ngập tại khu vực Đại Lộc và Điện Bàn.

Implications đa chiều

  • Về mặt học thuật và phương pháp luận: Khẳng định tính ưu việt của việc tích hợp mô hình thủy lực 1D với công nghệ GIS và mô hình kinh tế lượng trong việc giải bài toán quản lý rủi ro thiên tai tại các lưu vực sông ven biển miền Trung.
  • Về mặt ứng dụng thực tiễn: Cung cấp bản đồ ngập lụt số hóa chi tiết theo từng cấp độ sâu ngập (dưới $0,5\text{ m}$; $0,5 - 1,0\text{ m}$; $1,0 - 2,0\text{ m}$; và trên $2,0\text{ m}$), hỗ trợ chính quyền địa phương xác định chính xác các tọa độ di tản dân cư khẩn cấp.
  • Về mặt chính sách và quản trị: Cung cấp cơ sở khoa học định lượng cho Ban Chỉ huy Phòng chống thiên tai và Tìm kiếm cứu nạn tỉnh Quảng Nam trong việc:
    • Quy hoạch hành lang thoát lũ và điều chỉnh quy hoạch phân vùng sử dụng đất nông nghiệp, công nghiệp.
    • Xây dựng quy trình vận hành liên hồ chứa trên lưu vực sông Vu Gia - Thu Bồn (hồ A Vương, Đak Mi 4, Sông Bung 4, Sông Tranh 2) nhằm cắt giảm đỉnh lũ cho hạ du.
    • Thiết kế khẩu độ cầu cống, nâng cấp hệ thống giao thông đường bộ tránh ngập úng cắt đứt cứu trợ cứu nạn.

Limitations và Future Research

Mặc dù đạt được những đóng góp học thuật và thực tiễn vững chắc, luận án ghi nhận một số giới hạn nghiên cứu (boundary conditions):

  1. Giới hạn không gian chiều của mô hình thủy lực: Việc áp dụng mô hình 1 chiều MIKE 11 kết hợp ô chứa giả 2 chiều tuy tối ưu về mặt tốc độ tính toán nhưng chưa mô tả chi tiết trường vận tốc dòng chảy 2 chiều ($u, v$) tại các vùng giao cắt phân lưu phức tạp hoặc các khu vực đô thị hóa mật độ cao (phố cổ Hội An, thị xã Điện Bàn).
  2. Độ phân giải của mô hình số độ cao (DEM): Dữ liệu địa hình DEM sử dụng dựa trên bản đồ địa hình tỷ lệ 1:10.000 đến 1:25.000, có thể bỏ sót một số công trình vi địa hình như bờ vùng, bờ thửa, tôn nền cục bộ của các khu dân cư mới.
  3. Giới hạn số liệu thiệt hại kinh tế gián tiếp: Số liệu tổn thất kinh tế quy đổi chủ yếu tập trung vào thiệt hại trực tiếp (nhà cửa, lúa màu, đường sá), chưa bao quát toàn diện các tổn thất gián tiếp như gián đoạn kinh doanh, đình trệ chuỗi cung ứng công nghiệp và tổn thương sức khỏe tâm lý cộng đồng.

Hướng nghiên cứu tiếp theo (Future Research Agenda):

  • Phát triển mô hình thủy lực khớp nối 1D-2D hoàn chỉnh (MIKE FLOOD hoặc HEC-RAS 2D) tích hợp ảnh viễn thám độ phân giải siêu cao (LiDAR DEM) để mô phỏng chi tiết dòng tràn đô thị.
  • Mở rộng nghiên cứu tác động kép của biến đổi khí hậu kết hợp nước biển dâng và các kịch bản vận hành xả lũ cực đoan của toàn bộ hệ thống bậc thang thủy điện thượng nguồn.
  • Xây dựng hệ thống hỗ trợ ra quyết định thời gian thực (Real-time DSS) tích hợp trí tuệ nhân tạo (AI/Machine Learning) dự báo ngập lụt tức thời từ ảnh radar thời tiết.

Tác động và ảnh hưởng

Luận án tạo ra tác động sâu rộng trên nhiều phương diện:

  • Tác động học thuật: Đóng góp một nghiên cứu điển hình (case study) chuẩn mực về thủy văn - thủy lực lưu vực sông miền Trung, dự kiến mang lại tiềm năng trích dẫn cao trong các tạp chí chuyên ngành Khoa học Trái đất, Quản lý Tài nguyên nước và Kỹ thuật Công trình Thủy.
  • Tác động kinh tế - xã hội: Giúp chính quyền tỉnh Quảng Nam và TP. Đà Nẵng tối ưu hóa hàng trăm tỷ đồng ngân sách đầu tư công hàng năm cho các công trình phòng chống lũ, giảm thiểu tổn thất nhân mạng và tài sản cho hơn 1,5 triệu dân cư sinh sống trong lưu vực.
  • Tác động chính sách: Cung cấp cơ sở khoa học trực tiếp phục vụ việc điều chỉnh "Quy hoạch phòng chống lũ chi tiết hệ thống sông Vu Gia - Thu Bồn" và xây dựng bản đồ rủi ro thiên tai quốc gia theo Quyết định của Thủ tướng Chính phủ.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận một khung phương pháp luận kết hợp chuẩn xác giữa MIKE 11, GIS và lượng giá thiệt hại ECLAC cho sông ngòi miền Trung.
  • Cơ quan Quản lý Nhà nước & Quy hoạch: Ban Chỉ huy Phòng chống thiên tai & TKCN, Sở Nông nghiệp và PTNT, Sở Xây dựng tỉnh Quảng Nam và TP. Đà Nẵng có công cụ định lượng để phê duyệt dự án đầu tư và quy hoạch đô thị.
  • Doanh nghiệp và Ngành Bảo hiểm: Sử dụng bản đồ phân vùng ngập lụt và đường cong tần suất - thiệt hại để thẩm định rủi ro tài sản, định giá phí bảo hiểm thiên tai cho các khu công nghiệp, khu du lịch ven biển.
  • Cộng đồng dân cư địa phương: Hưởng lợi từ hệ thống cảnh báo sớm, các tuyến đường sơ tán an toàn và kế hoạch ứng phó thiên tai dựa vào cộng đồng.

Câu hỏi chuyên sâu

  1. Đóng góp lý thuyết độc đáo nhất của luận án là gì?
    Luận án đã mở rộng lý thuyết thủy lực Saint-Venant và khung lượng giá thiên tai ECLAC thông qua việc thiết lập thành công bộ hàm quan hệ phi tuyến thực chứng giữa mực nước đỉnh lũ $H_{\max}$ tại các trạm thủy văn khống chế với tổng thiệt hại kinh tế quy đổi đa ngành cho một lưu vực sông có cấu trúc phân lưu kép đặc thù.

  2. Tính đổi mới trong phương pháp nghiên cứu so với các nghiên cứu quốc tế trước đây?
    So với mô hình HEC-RAS/ArcHydro của Knebla et al. (2005) chỉ dừng ở ranh giới ngập hoặc mô hình VRSAP của Phạm Ngọc (2000) thiếu trực quan không gian, luận án đã thiết lập quy trình tích hợp khép kín: Mô phỏng mưa - dòng chảy (MIKE NAM) $\rightarrow$ Mô phỏng thủy lực dòng chảy hở 1D mạng sông phân lưu (MIKE 11 HD) $\rightarrow$ Không gian hóa ngập lụt trên nền DEM (MIKE 11 GIS) $\rightarrow$ Lượng giá kinh tế theo hàm tổn thất thực tế.

  3. Phát hiện bất ngờ nhất từ dữ liệu nghiên cứu là gì?
    Đó là tính chất phi đối xứng nghiêm trọng của đường quá trình lũ: thời gian truyền lũ từ thượng lưu về hạ lưu chỉ mất từ $3 - 7\text{ giờ}$ với cường suất lên tới $140\text{ cm/h}$, nhưng thời gian thoát lũ ra biển kéo dài gấp 3–5 lần do độ dốc hạ lưu quá nhỏ ($0,04‰$) kết hợp hiệu ứng nghẽn dòng triều cường tại Cửa Đại và Cửa Hàn, biến vùng đồng bằng Điện Bàn thành lòng hồ chứa lũ cưỡng bức.

  4. Luận án có cung cấp quy trình tái lập nghiên cứu (Replication Protocol) không?
    Hoàn toàn có. Luận án cung cấp chi tiết bộ thông số nhám Manning $n$ trên từng đoạn sông ($0,025 - 0,045$), tọa độ và dung tích các ô chứa lũ, các bước thiết lập sơ đồ mạng sông trong MIKE 11, quy trình liên kết dữ liệu với ArcGIS 9.1 và thuật toán quy đổi tổn thất kinh tế theo giá trị hiện tại ròng.

  5. Chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?
    Mở rộng mô hình hóa tích hợp toàn diện 1D-2D-3D kết nối lưu vực - cửa sông - dải ven bờ; nghiên cứu biến động bồi xói lòng dẫn sông và bờ biển dưới tác động của suy giảm bùn cát do hệ thống thủy điện thượng nguồn; phát triển nền tảng bản sao số thủy văn (Digital Twin) hỗ trợ điều hành thiên tai thời gian thực.

Kết luận

  1. Luận án đã làm sáng tỏ toàn diện các đặc trưng địa lý tự nhiên, khí tượng thủy văn của lưu vực sông Vu Gia - Thu Bồn ($10.350\text{ km}^2$), đặc biệt là cơ chế hình thành lũ núi cực đoan với biên độ lên tới $15,2\text{ m}$ và vận tốc dòng chảy vượt $3,7\text{ m/s}$.
  2. Thiết lập và kiểm định thành công bộ mô hình toán thủy văn - thủy lực MIKE NAM và MIKE 11 trên mạng lưới sông phân lưu phức tạp, đạt độ tin cậy và chuẩn xác khoa học cao qua các trận lũ lịch sử 1998, 1999, 2002 và 2004.
  3. Ứng dụng xuất sắc công nghệ MIKE 11 GIS để xây dựng bộ bản đồ ngập lụt chi tiết ứng với 4 kịch bản tần suất lũ thiết kế $p = 1%$, $2%$, $5%$, và $10%$, phân hóa cụ thể diện tích và độ sâu ngập trên địa bàn trọng điểm huyện Điện Bàn.
  4. Lượng hóa thành công các hàm quan hệ thiệt hại kinh tế quy đổi (mặt bằng giá 2014) theo mực nước đỉnh lũ $H_{\max}$, mang lại bước tiến vượt bậc trong việc gắn kết mô hình thủy lực với kinh tế học thiên tai.
  5. Đề xuất hệ thống giải pháp phi công trình và công trình mang tính chiến lược, cung cấp luận cứ khoa học trực tiếp cho các nhà hoạch định chính sách trong quản lý bền vững tài nguyên nước và giảm nhẹ rủi ro thiên tai tại khu vực miền Trung Việt Nam.