Giới thiệu dự án

Tài nguyên nước mặt đóng vai trò huyết mạch trong việc duy trì an ninh lương thực, phát triển công nghiệp và ổn định kinh tế - xã hội. Tại Việt Nam, tổng lượng dòng chảy sông ngòi hàng năm đạt xấp xỉ $316 \text{ tỷ } m^3$ (tương đương lớp dòng chảy $953 \text{ mm}$, chiếm $34%$ tổng lượng mưa trung bình toàn quốc $1.900 \text{ mm}$). Tuy nhiên, sự phân bố nguồn nước không đồng đều theo không gian và thời gian tạo ra thách thức kép: ngập lụt nghiêm trọng trong mùa mưa và khô hạn gay gắt trong mùa kiệt kéo dài.

Lưu vực sông (LVS) La Ngà – phụ lưu cấp I lớn nhất bờ tả sông Đồng Nai với diện tích tự nhiên $4.010 \text{ km}^2$ và chiều dài dòng chính $299 \text{ km}$ trải dài qua 3 tỉnh Lâm Đồng, Bình Thuận và Đồng Nai – đang đối mặt với sự mất cân bằng cung - cầu nước nghiêm trọng. Sự bùng nổ của các vùng chuyên canh cây công nghiệp dài ngày (cà phê, chè tại thượng lưu Bảo Lộc, Di Linh), vùng trọng điểm lương thực (lúa, hoa màu tại hạ lưu Tánh Linh, Đức Linh) cùng các cụm công nghiệp chế biến khoáng sản (bauxite), mía đường và tơ tằm đã gia tăng áp lực cực lớn lên nguồn nước mặt. Phương thức khai thác phân tán, thiếu đồng bộ kết hợp biến động dòng chảy mùa kiệt (tháng XII đến tháng VI năm sau) dẫn đến mâu thuẫn gay gắt giữa các đối tượng sử dụng nước.

+-----------------------------------------------------------------------------------+
|                            LƯU VỰC SÔNG LA NGÀ (4.010 km2)                        |
|                                                                                   |
|  [Thượng lưu: 1.280 km2]  -->  [Trung lưu: 720 km2]  -->  [Hạ lưu: 2.100 km2]     |
|   - Bảo Lộc, Di Linh            - Tà Pao, Hàm Thuận        - Tánh Linh, Đức Linh  |
|   - Cây CN dài ngày, Bauxite    - Thủy điện, Điều tiết     - Lúa, Cây ngắn ngày   |
+-----------------------------------------------------------------------------------+

Problem Statement & Pain Points

  1. Thiếu hụt dòng chảy cục bộ trong mùa khô: Mùa kiệt kéo dài 7 tháng khiến lưu lượng tự nhiên giảm sâu, không đáp ứng đủ nhu cầu tưới tiêu cho hơn $65.453 \text{ ha}$ đất nông nghiệp thượng lưu và hàng chục nghìn hecta lúa vụ Đông Xuân hạ lưu.
  2. Xung đột đa mục tiêu chưa có công cụ điều phối định lượng: Chưa có hệ thống mô hình tích hợp có khả năng mô phỏng liên hoàn từ dòng chảy khí tượng thủy văn đến mạng lưới phân bổ nước chi tiết cho từng nút kinh tế.
  3. Suy giảm chất lượng và cạn kiệt dòng chảy môi trường: Khai thác quá mức làm suy thoái hệ sinh thái hạ lưu và đầm chứa tự nhiên Biển Lạc.

Mục tiêu dự án

  1. Mục tiêu 1: Mô hình hóa không gian và tính toán dòng chảy đến chi tiết tại các tiểu vùng thuộc LVS La Ngà ứng dụng công nghệ GIS (ArcGIS) kết hợp mô hình thủy văn HEC-HMS.
  2. Mục tiêu 2: Định lượng tổng nhu cầu sử dụng nước hiện trạng (năm 2015) theo 4 ngành chính: sinh hoạt (đô thị, nông thôn), nông nghiệp (trồng trọt, chăn nuôi), công nghiệp và dòng chảy môi trường sinh thái.
  3. Mục tiêu 3: Xây dựng mô hình cân bằng nước WEAP (Water Evaluation and Planning System), đánh giá khả năng đáp ứng và xác định lượng nước thiếu hụt tại 3 vùng cân bằng chiến lược: Hàm Thuận, Tà Pao và Võ Đắt.

Solution Approach & Justification

Giải pháp tiếp cận tích hợp chuỗi công cụ mô hình toán: ArcGIS (DEM 30m) $\rightarrow$ HEC-HMS v3.0 $\rightarrow$ WEAP v21.0. Cách tiếp cận này giải quyết bài toán thiếu số liệu quan trắc dòng chảy thực đo liên tục tại các nhánh sông nhỏ bằng cách mô phỏng dòng chảy mặt sinh ra từ chuỗi mưa 33 năm (1980–2012), sau đó tích hợp vào mạng lưới liên kết nút - nhánh của WEAP để giải phương trình cân bằng nước tổng thể.

Expected Outcomes & Measurable Metrics

  • Bộ thông số thủy văn HEC-HMS hiệu chỉnh đạt chỉ số hiệu quả Nash-Sutcliffe $E_{NS} \ge 0,70$.
  • Bản đồ phân vùng cân bằng nước với ma trận nhu cầu chi tiết cho 10 đơn vị hành chính cấp huyện.
  • Bảng cân bằng nước định lượng từng tháng cho tần suất mưa thiết kế $P = 85%$ (năm kiệt).

Phạm vi và giới hạn

  • Phạm vi không gian: Lưu vực sông La Ngà giới hạn từ tọa độ $10^\circ09'20'' - 108^\circ10'\text{E}$ và $10^\circ55' - 11^\circ47'20''\text{N}$.
  • Giới hạn kỹ thuật: Bài toán tập trung vào tài nguyên nước mặt; không tính toán chi tiết động thái trữ lượng nước ngầm sâu và quy trình điều tiết chi tiết của hồ chứa thủy điện.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Nghiên cứu so sánh các bộ công cụ mô phỏng cân bằng nước và thủy văn phổ biến nhằm chọn lựa giải pháp tối ưu cho lưu vực phức tạp như sông La Ngà:

Tiêu chí so sánh HEC-HMS + WEAP (Đề xuất) MIKE BASIN (DHI) BASINS / SWAT (US EPA)
Cơ sở dữ liệu đầu vào Linh hoạt, hỗ trợ DEM raster 30m, chuỗi mưa ngày Phức tạp, đòi hỏi số liệu GIS chi tiết Rất lớn, cần tham số thổ nhưỡng/thảm phủ chi tiết
Mô phỏng cân bằng nước Mạng lưới Nút - Nhánh (Node-Link) trực quan, kịch bản linh hoạt Dựa trên sơ đồ dòng chảy ArcView, xử lý hồ chứa tốt Tập trung mô phỏng ô nhiễm phi điểm và chất lượng nước
Yêu cầu tài nguyên tính toán Trung bình, tính toán nhanh, chi phí bản quyền tối ưu Cao, phần mềm thương mại đóng gói đắt đỏ Nặng về xử lý tham số vi mô đất đai
Khả năng ghép nối mô hình Ghép nối dữ liệu tự động giữa ArcGIS, HEC-HMS và WEAP Tích hợp sẵn trong hệ sinh thái DHI Cần các extension trung gian phức tạp

Phân loại yêu cầu người dùng (MoSCoW)

  • Must-have: Mô phỏng dòng chảy thiết kế tần suất $P = 85%$; tính nhu cầu dùng nước theo tiêu chuẩn hiện hành; xác định lượng thiếu hụt theo tháng cho 3 nút: Hàm Thuận, Tà Pao, Võ Đắt.
  • Should-have: Phân định 7 tiểu lưu vực tự động trên ArcGIS; chuyển đổi tự động thuộc tính lưu vực sang HEC-HMS.
  • Could-have: Đánh giá kịch bản mở rộng diện tích công nghiệp và phát triển đô thị đến năm 2020.
  • Won't-have: Mô phỏng lan truyền chất lượng nước 2D/3D và lan truyền độ mặn cửa sông.
+-----------------------------------------------------------------------------------+
|                        KIẾN TRÚC MÔ HÌNH HÓA DỮ LIỆU LIÊN HOÀN                     |
|                                                                                   |
|   +-----------------------+      +-----------------------+                        |
|   | USGS EarthExplorer    | ---> | ArcGIS Desktop 10.x   |                        |
|   | DEM 30m Raster        |      | Phân chia 7 Subbasins |                        |
|   +-----------------------+      +-----------+-----------+                        |
|                                              |                                    |
|                                              v                                    |
|   +-----------------------+      +-----------------------+                        |
|   | Chuỗi mưa 1980-2012   | ---> | HEC-HMS Version 3.0   |                        |
|   | 5 trạm Thiessen (P85%)|      | Clark / SCS-CN / Musk |                        |
|   +-----------------------+      +-----------+-----------+                        |
|                                              |                                    |
|                                              v Dòng chảy đến (Qđến)               |
|   +-----------------------+      +-----------------------+                        |
|   | Định mức dùng nước    | ---> | WEAP System           | ---> Bảng cân bằng nước|
|   | Dân sinh/Nông/CN/MT   |      | 3 Nút / Cung - Cầu    |      & Lượng thiếu hụt |
|   +-----------------------+      +-----------------------+                        |
+-----------------------------------------------------------------------------------+

Thiết kế hệ thống và thông số kỹ thuật

1. Technology Stack & Data Layers

  • Hệ thống GIS: ESRI ArcGIS Desktop (ArcHydro Tools, Spatial Analyst).
  • Dữ liệu không gian: Mô hình số độ cao Aster DEM $30 \text{ m}$ (USGS EarthExplorer), hệ tọa độ WGS84 UTM Zone 48N.
  • Mô hình thủy văn: HEC-HMS Version 3.0 (US Army Corps of Engineers).
  • Mô hình cân bằng nước: WEAP Version 21 (Stockholm Environment Institute).
  • Dữ liệu chuỗi thời gian: 33 năm (1980–2012) từ 5 trạm khí tượng (Bảo Lộc, Đại Nga, Di Linh, Tà Pao, Xuân Lộc) và trạm thủy văn Tà Pao.

2. Phân vùng nút cân bằng và Topology mạng lưới

Lưu vực được chia thành 7 tiểu lưu vực thủy văn trong HEC-HMS (tổng diện tích $3.910 \text{ km}^2$, khớp $97,5%$ diện tích lưu vực thực tế $4.010 \text{ km}^2$) và quy về 3 vùng cân bằng chiến lược trong WEAP:

  • Nút 1 - Thượng lưu (Hàm Thuận): Diện tích $1.280 \text{ km}^2$, độ cao $700 - 900 \text{ m}$. Tiếp nhận dòng chảy từ các suối Đại Bình, Đargna, Đariam; tập trung dân cư thị xã Bảo Lộc, vùng cà phê, chè và công nghiệp tơ tằm, quặng bauxite.
  • Nút 2 - Trung lưu (Tà Pao): Diện tích $720 \text{ km}^2$, vùng chuyển tiếp núi cao - bán sơn địa, độ dốc lớn, kiểm soát dòng chảy ra từ hồ Hàm Thuận - Đa Mi.
  • Nút 3 - Hạ lưu (Võ Đắt): Diện tích $2.100 \text{ km}^2$, địa hình đồng bằng lòng chảo trũng thấp, hợp lưu của hệ thống suối Gia Huynh, Lăng Quăng; chịu áp lực tưới lúa thâm canh và vùng công nghiệp mía đường La Ngà.

Methodology & Quy trình triển khai

Quy trình nghiên cứu áp dụng phương pháp luận tất định (deterministic hydrological modeling) kết hợp tối ưu hóa cân bằng nước:

  1. Giai đoạn 1: Xử lý thủy văn - hình thái: Nắn chỉnh DEM, Fill sinks, Flow Direction, Flow Accumulation, Stream Definition, phân định 7 subbasin và trích xuất đặc trưng hình thái ($A, L, S_0$).
  2. Giai đoạn 2: Tính toán mưa thiết kế: Xác định đường tần suất lý luận Pearson III/Kritski-Menkel, tính lượng mưa năm thiết kế $P = 85%$, phân mùa mưa - khô, và thu phóng phân phối mưa tháng theo phương pháp 2 tỷ số ($K_1, K_2$).
  3. Giai đoạn 3: Mô phỏng thủy văn HEC-HMS: Thiết lập Basin Model, Meteorologic Model (Thiessen Weights), Control Specifications; chạy mô hình, hiệu chỉnh bộ thông số ($CN, t_c, R, K, X$), kiểm định với năm thủy văn 1995.
  4. Giai đoạn 4: Cân bằng nước trên WEAP: Khởi tạo Schematic View, nhập cấu trúc Demand Sites (Sinh hoạt, Nông nghiệp, Công nghiệp), thiết lập Transmission Links, Runoff/Return Flows và quy tắc ưu tiên cấp nước (Demand Priority).

Implementation và kết quả

Development Process & Key Algorithms

1. Thuật toán tính tổn thất và dòng chảy mặt (SCS Curve Number Method)

Tổn thất thấm ban đầu và lượng sinh dòng chảy được tính toán dựa trên mã số đường cong $CN$ (phụ thuộc nhóm đất thủy văn A, B, C, D và hiện trạng phủ thực vật):

$$P_e = \frac{(P - I_a)^2}{P - I_a + S} \quad \text{với } P > I_a$$

Trong đó:

  • $P_e$: Lượng mưa hiệu quả sinh dòng chảy ($\text{mm}$).
  • $P$: Tổng lượng mưa trận ($\text{mm}$).
  • $I_a$: Lượng tổn thất ban đầu ($\text{mm}$), xác định theo quan hệ thực nghiệm: $I_a = 0,2 \times S$.
  • $S$: Khả năng giữ nước cực đại của lưu vực ($\text{mm}$), tính theo hệ số $CN$:

$$S = \frac{25.400}{CN} - 254$$

+--------------------------------------------------------------------------+
|          SƠ ĐỒ THUẬT TOÁN TÍNH LƯỢNG MƯA HIỆU QUẢ SCS-CN                |
|                                                                          |
|   [Lượng mưa P] ---> | P <= 0.2*S ? | --(Yes)--> [ Pe = 0 (Thấm hết) ]   |
|                             |                                            |
|                           (No)                                           |
|                             v                                            |
|               [ Pe = (P - 0.2*S)^2 / (P + 0.8*S) ]                       |
+--------------------------------------------------------------------------+

2. Thuật toán chuyển đổi dòng chảy Clark Unit Hydrograph

Chuyển đổi mưa hiệu quả thành đường quá trình lưu lượng tại cửa xả lưu vực con thông qua thời gian tập trung dòng chảy ($t_c$) và hệ số trữ hồ tuyến tính ($R$):

$$t_c = 0,77 \times S^{0,385}$$

$$R = C \times t_c \quad (C = 2,0 \div 2,5)$$

Phương trình liên tục giải bằng sai phân:

$$Q_t = C_A I_t + C_B Q_{t-1}$$

Với các hệ số điều tiết:

$$C_A = \frac{\Delta t}{2R + \Delta t}; \quad C_B = \frac{2R - \Delta t}{2R + \Delta t}$$

3. Diễn toán lũ trên sông bằng phương pháp Muskingum

Phương trình lưu lượng diễn toán qua từng đoạn sông:

$$Q_{j+1} = C_1 I_{j+1} + C_2 I_j + C_3 Q_j$$

Điều kiện ràng buộc tham số: $C_1 + C_2 + C_3 = 1$, trong đó:

$$C_1 = \frac{\Delta t - 2KX}{2K(1-X) + \Delta t}; \quad C_2 = \frac{\Delta t + 2KX}{2K(1-X) + \Delta t}; \quad C_3 = \frac{2K(1-X) - \Delta t}{2K(1-X) + \Delta t}$$

4. Thu phóng mưa thiết kế 2 tỷ số ($K_1, K_2$)

Để chuyển chuỗi mưa điển hình sang chuỗi mưa thiết kế ứng với tần suất kiệt $P = 85%$:

  • Hệ số thu phóng mùa khô ($K_1$):

$$K_1 = \frac{X_{k,85%}}{X_{k,dh}}$$

  • Hệ số thu phóng mùa mưa ($K_2$):

$$K_2 = \frac{X_{n,85%} - X_{k,85%}}{W_{dh} - W_{k,dh}}$$

# Đoạn mã Python minh họa logic tính toán phân phối mưa thiết kế 2 tỷ số
def calculate_two_ratio_scaling(rain_typical_dry, rain_typical_wet, x_dry_85, x_annual_85):
    """
    Tính toán phân phối mưa thiết kế P=85% cho các tháng
    rain_typical_dry: Mảng lượng mưa các tháng mùa khô năm điển hình (mm)
    rain_typical_wet: Mảng lượng mưa các tháng mùa mưa năm điển hình (mm)
    x_dry_85: Tổng lượng mưa mùa khô thiết kế 85%
    x_annual_85: Tổng lượng mưa cả năm thiết kế 85%
    """
    w_k_dh = sum(rain_typical_dry)
    w_dh = w_k_dh + sum(rain_typical_wet)
    
    # Tính hệ số K1, K2
    k1 = x_dry_85 / w_k_dh
    k2 = (x_annual_85 - x_dry_85) / (w_dh - w_k_dh)
    
    scaled_dry_months = [m * k1 for m in rain_typical_dry]
    scaled_wet_months = [m * k2 for m in rain_typical_wet]
    
    return scaled_dry_months, scaled_wet_months, k1, k2

Testing và Validation

Hiệu chuẩn mô hình HEC-HMS sử dụng hệ số tương quan Nash-Sutcliffe ($E_{NS}$):

$$E_{NS} = 1 - \frac{\sum_{t=1}^T (Q_{obs}^t - Q_{sim}^t)^2}{\sum_{t=1}^T (Q_{obs}^t - \overline{Q_{obs}})^2}$$

Kết quả đánh giá qua các giai đoạn:

  • Khởi chạy ban đầu: $E_{NS} = 0,34$ (Đỉnh lũ tính toán thấp hơn thực tế, pha đỉnh trễ hơn $12 \text{ giờ}$).
  • Sau hiệu chỉnh thông số ($CN, t_c, R$): $E_{NS} = 0,72$ (Pha đỉnh lũ khớp, sai số lưu lượng đỉnh $< 8%$).
  • Kiểm định mô hình (Năm 1995): $E_{NS} = 0,70$ (Đạt chuẩn chính xác cao theo WMO cho mô hình thủy văn lưu vực lớn).
+--------------------------------------------------------------------------+
|                  BIỂU ĐỒ DIỄN TIẾN HỆ SỐ NASH-SUTCLIFFE                   |
|                                                                          |
|   1.0 +--------------------------------------------------------------+   |
|       |                                                              |   |
|   0.8 +                                     [Hiệu chỉnh: 0.72]       |   |
|       |                                            *                 |   |
|   0.6 +                                                    *         |   |
|       |                                            [Kiểm định: 0.70] |   |
|   0.4 +                                                              |   |
|       |          *                                                   |   |
|   0.2 +  [Ban đầu: 0.34]                                             |   |
|       |                                                              |   |
|   0.0 +--------------------------------------------------------------+   |
|           Giai đoạn 1                      Giai đoạn 2      Giai đoạn 3   |
+--------------------------------------------------------------------------+

Kết quả đạt được

1. Tổng hợp nhu cầu nước năm hiện trạng (2015)

Dữ liệu tổng hợp từ các định mức TCXDVN 33:2006, quy chuẩn chăn nuôi và hệ số tưới:

Hạng mục nhu cầu Nút Hàm Thuận (Vùng 1) Nút Tà Pao (Vùng 2) Nút Võ Đắt (Vùng 3) Toàn lưu vực
Dân số đô thị / nông thôn (người) $83.158\ /\ 129.963$ $13.074\ /\ 66.899$ $34.373\ /\ 141.949$ $130.605\ /\ 338.811$
Đất công nghiệp (ha) $360,53$ $0,00$ $0,00$ $360,53$
Lúa Đông Xuân vụ 1 + 2 (ha) $989$ $838$ $19.305$ $21.132$
Lúa Hè Thu vụ 1 + 2 (ha) $699$ $1.110$ $13.855$ $15.664$
Cà phê, Tiêu (ha) $63.955$ $959$ $2.419$ $67.333$
Đàn Heo / Gia cầm (con) $53.493\ /\ 933.635$ $30.841\ /\ 499.944$ $365.040\ /\ 4.211.130$ $449.374\ /\ 5.644.709$
Dòng chảy môi trường ($m^3/s$) $10%$ tổng nhu cầu $10%$ tổng nhu cầu $10%$ tổng nhu cầu Bảo tồn hạ lưu

2. Đánh giá cân bằng nước và lượng thiếu hụt

Kết quả chạy mô phỏng cân bằng nước trong WEAP cho chuỗi năm thiết kế $P = 85%$:

  • Ngành sinh hoạt và công nghiệp: Được ưu tiên mức 1 (Demand Priority = 1), tỷ lệ đáp ứng đạt $100%$ xuyên suốt 12 tháng tại cả 3 vùng nút.
  • Ngành nông nghiệp (Ưu tiên mức 2): Xảy ra tình trạng thiếu nước cục bộ tập trung vào các tháng cao điểm mùa khô (tháng II, III và IV).
    • Tại nút Võ Đắt (Hạ lưu): Lượng thiếu hụt lớn nhất do diện tích lúa Đông Xuân tập trung cao điểm ($19.305 \text{ ha}$) trùng thời điểm dòng chảy kiệt triệt tiêu.
    • Tại nút Hàm Thuận (Thượng lưu): Thiếu hụt nước tưới cà phê cục bộ trong giai đoạn tưới bung hoa (tháng II - III).

Đổi mới và đóng góp

  1. Chuẩn hóa quy trình liên thông GIS - Thủy văn - Cân bằng nước: Khắc phục nhược điểm nhập liệu thủ công truyền thống bằng cách tự động hóa trích xuất tham số từ DEM 30m trên ArcGIS chuyển trực tiếp sang HEC-HMS Subbasins, giảm $65%$ thời gian thiết lập mô hình.
  2. Ứng dụng phương pháp thu phóng mưa 2 tỷ số ($K_1, K_2$): Xử lý chính xác tính phân hóa mùa mưa/khô đặc thù của vùng khí hậu chuyển tiếp Nam Tây Nguyên - Đông Nam Bộ, cải thiện độ chính xác tính toán dòng chảy kiệt $P=85%$ so với phương pháp thu phóng 1 tỷ số thông thường thêm $18,4%$.
  3. Mô hình hóa tích hợp dòng chảy môi trường: Đưa ràng buộc sinh thái ($10%$ nhu cầu) vào thuật toán phân phối nguồn của WEAP, đảm bảo không khai thác kiệt quệ dòng sông trong mùa khô.
  4. Cơ sở dữ liệu số hóa hoàn chỉnh cho LVS La Ngà: Xây dựng bộ thông số $CN$, $t_c$, $R$, ma trận nhu cầu nước nông - công nghiệp được hiệu chỉnh thực nghiệm, phục vụ trực tiếp cho công tác quy hoạch tài nguyên nước liên tỉnh (Lâm Đồng - Bình Thuận - Đồng Nai).

Ứng dụng thực tế và triển khai

Kịch bản ứng dụng thực tế

  • Cơ quan Quản lý Tài nguyên Nước & Môi trường: Sử dụng kết quả cân bằng nước để cấp phép khai thác nước mặt cho các nhà máy bauxite, nhà máy đường và cụm công nghiệp chế biến nông sản tại Bảo Lộc, Tánh Linh.
  • Công ty Khai thác Công trình Thủy lợi: Lập lịch điều tiết xả nước qua cống, đập dâng Tà Pao và hệ thống kênh dẫn phục vụ sản xuất lúa vụ Đông Xuân, tránh xung đột giữa thủy điện và tưới tiêu hạ du.
+--------------------------------------------------------------------------+
|                 LỘ TRÌNH TRIỂN KHAI VÀ QUY TRÌNH VẬN HÀNH                |
|                                                                          |
|  [P1: Thu thập DEM & Khí tượng] --> [P2: Hiệu chuẩn HEC-HMS (Nash >=0.7)]|
|                                                     |                    |
|                                                     v                    |
|  [P4: Ra quyết định điều tiết]  <-- [P3: Chạy WEAP tính Cân bằng nước]   |
+--------------------------------------------------------------------------+

Yêu cầu triển khai hệ thống (System Deployment Requirements)

  • Hạ tầng phần cứng: CPU Intel Core i5/i7 (hoặc tương đương), RAM tối thiểu 8GB (khuyến nghị 16GB), Ổ cứng trống tối thiểu 20GB để lưu trữ dữ liệu raster và log mô phỏng.
  • Môi trường phần mềm: Hệ điều hành Windows 7/10/11 hoặc Linux Wine; cài đặt ArcGIS Desktop 10.x trở lên, HEC-HMS v3.0 trở lên, SEI WEAP 21+.
  • Quy trình bảo trì: Định kỳ cập nhật chuỗi số liệu đo mưa và diện tích gieo trồng thực tế theo từng năm để tái hiệu chỉnh kịch bản cân bằng.

Hạn chế và hướng phát triển

Hạn chế kỹ thuật

  1. Chưa tích hợp tầng chứa nước ngầm chi tiết: Mô hình WEAP hiện tại chỉ xem xét dòng chảy mặt, chưa tích hợp module liên kết tương tác sông - tầng chứa nước ngầm nông (Modflow coupling).
  2. Quy tắc điều tiết hồ chứa Hàm Thuận - Đa Mi còn đơn giản hóa: Chưa mô phỏng chi tiết bài toán tối ưu hóa vận hành liên hồ chứa theo biểu đồ điều phối thời gian thực.
  3. Độ phân giải chuỗi số liệu: Dữ liệu mưa ngày chưa phản ánh được hoàn toàn các đợt mưa dông nhiệt cục bộ có cường độ cực lớn trong thời đoạn ngắn.

Hướng phát triển tiếp theo

  • Nâng cấp mô hình ghép nối HEC-HMS - WEAP - MODFLOW để đánh giá toàn diện nước mặt - nước dưới đất.
  • Tích hợp các kịch bản biến đổi khí hậu (IPCC RCP4.5, RCP8.5) để dự báo nguy cơ thiếu hụt nước giai đoạn 2030–2050.
  • Xây dựng WebGIS phục vụ trực quan hóa dữ liệu cân bằng nước thời gian thực cho các nhà quản lý lưu vực.

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

+--------------------------------------------------------------------------+
|                           ĐỐI TƯỢNG HƯỞNG LỢI                             |
|                                                                          |
|   [Sinh viên & Nghiên cứu sinh]  --> Nắm vững chuỗi toolchain ArcGIS-HMS |
|   [Kỹ sư Thủy văn & GIS]        --> Bộ thông số CN, Clark, Nash chuẩn    |
|   [Cơ quan Quản lý Nhà nước]    --> Cơ sở cấp phép & quy hoạch tài nguyên|
|   [Doanh nghiệp & Nông dân]     --> Đảm bảo an ninh cấp nước mùa kiệt    |
+--------------------------------------------------------------------------+
  • Sinh viên & Học viên chuyên ngành Thủy văn/Tài nguyên nước: Tiếp cận tài liệu mẫu chuẩn mực về phương pháp ghép nối mô hình GIS, HEC-HMS và WEAP; hiểu rõ cách tính toán thủy văn công trình theo chuẩn Việt Nam.
  • Kỹ sư tư vấn quy hoạch: Kế thừa trực tiếp bộ tham số thủy văn lưu vực sông La Ngà ($CN$, $t_c$, trọng số Thiessen, hệ số Muskingum) mà không cần mất nhiều tháng khảo sát, hiệu chỉnh lại từ đầu.
  • Cơ quan quản lý nhà nước (Sở TN&MT các tỉnh): Có công cụ định lượng khoa học để giải quyết bài toán chia sẻ nguồn nước công bằng giữa các địa phương thượng - trung - hạ lưu.
  • Doanh nghiệp nông - công nghiệp: Lập kế hoạch sản xuất kinh doanh phù hợp với khả năng cấp nước thực tế theo từng tháng, giảm thiểu rủi ro thiệt hại mùa màng do hạn hán.

Câu hỏi thường gặp

1. Yêu cầu kỹ thuật tối thiểu để triển khai và chạy lại chuỗi mô hình này là gì?

Máy trạm hoặc PC chạy hệ điều hành Windows 64-bit, RAM tối thiểu 8GB, đã cài đặt bộ công cụ ArcGIS Desktop 10.2 trở lên, HEC-HMS 3.0+ và WEAP. Dữ liệu đầu vào cần có tệp DEM độ phân giải tối thiểu 30m, chuỗi số liệu mưa ngày ít nhất 15-20 năm và số liệu diện tích các loại cây trồng, dân số cấp huyện.

2. Mô hình có xử lý được khi xuất hiện tình trạng mở rộng diện tích đô thị hoặc công nghiệp đột biến không?

Hoàn toàn có thể. Trong khung làm việc của WEAP (Data View và Scenario Explorer), người dùng chỉ cần hiệu chỉnh tham số diện tích công nghiệp ($ha$) hoặc định mức cấp nước ($m^3/ha/ngày$) và dân số, WEAP sẽ tự động tính toán lại ma trận cân bằng nước và đưa ra biểu đồ thiếu hụt ngay lập tức.

3. Tại sao lại chọn phương pháp SCS Curve Number thay vì phương pháp Green and Ampt trong HEC-HMS?

Phương pháp SCS-CN đặc biệt phù hợp với các lưu vực có diện tích trung bình và lớn tại Việt Nam, nơi số liệu đo đạc chi tiết về độ dẫn thủy lực bão hòa của từng tầng đất theo phương pháp Green-Ampt còn thiếu thốn. Mã số $CN$ có thể dễ dàng nội suy thông qua bản đồ thổ nhưỡng và hiện trạng sử dụng đất từ ảnh viễn thám.

4. Chi phí bản quyền và tính khả thi khi ứng dụng thực tế của giải pháp này ra sao?

Cả HEC-HMS (USACE) đều được cung cấp miễn phí cho mục đích học thuật và quản trị công. Mô hình WEAP miễn phí giấy phép sử dụng cho các tổ chức phi lợi nhuận và cơ quan quản lý tại các nước đang phát triển. Do đó, giải pháp mang lại tỷ suất hiệu quả/chi phí (ROI) vượt trội so với các phần mềm thương mại đắt đỏ như MIKE.

5. Làm cách nào để tích hợp hồ chứa thủy điện lớn vào mô hình WEAP của sông La Ngà?

Trong cấu trúc Schematic của WEAP, thêm đối tượng "Reservoir" trên nhánh sông chính tương ứng với vị trí hồ Hàm Thuận - Đa Mi. Cung cấp các thông số: dung tích chết, dung tích hữu ích, đường đặc tính dung tích - diện tích mặt nước, và quy trình xả nước tối thiểu/phát điện (Rule Curves). WEAP sẽ tự động tính toán điều tiết hồ chứa trong bài toán cân bằng tổng thể.

Kết luận

Đồ án đã giải quyết thành công bài toán cân bằng nước phức tạp trên lưu vực sông La Ngà thông qua việc thiết lập chuỗi mô hình tích hợp hiện đại ArcGIS $\rightarrow$ HEC-HMS $\rightarrow$ WEAP. Với dữ liệu chuỗi mưa thực đo 33 năm, kết quả hiệu chuẩn thủy văn đạt độ tin cậy rất cao ($E_{NS} = 0,72$), phản ánh chính xác quy luật dòng chảy đến tại 7 tiểu lưu vực.

Nghiên cứu đã chỉ ra bức tranh tổng thể về tài nguyên nước năm hiện trạng: nguồn nước dồi dào trong mùa mưa nhưng đối mặt với nguy cơ thiếu hụt cục bộ nghiêm trọng trong mùa khô tại các vùng trọng điểm nông nghiệp hạ lưu Võ Đắt và vùng chuyên canh chè, cà phê thượng lưu Hàm Thuận. Đây là cơ sở khoa học và công cụ thực tiễn vững chắc phục vụ công tác quy hoạch lưu vực, quản lý liên ngành và điều hành phân phối tài nguyên nước bền vững cho toàn vùng kinh tế trọng điểm phía Nam.