Tổng quan về luận án

Quản lý tài nguyên nước bền vững tại các lưu vực sông ven biển miền Trung Việt Nam đang đứng trước những thách thức mang tính sống còn dưới áp lực kép từ gia tăng dân số, công nghiệp hóa và biến đổi khí hậu (BĐKH). Luận án tiến sĩ kỹ thuật chuyên ngành Quy hoạch và Quản lý Tài nguyên nước (mã số 62623001) của nghiên cứu sinh Lê Đức Thường, thực hiện tại Trường Đại học Bách khoa – Đại học Quốc gia TP. Hồ Chí Minh dưới sự hướng dẫn khoa học của GS. Trần Thục, PGS. Nguyễn Thống và TS. Châu Nguyễn Xuân Quang, là một công trình tiên phong giải quyết bài toán quản trị nước liên vùng có mức độ phức tạp cao. Lưu vực sông (LVS) Ba là hệ thống sông lớn nhất vùng ven biển miền Trung với tổng diện tích tự nhiên $13.417\text{ km}^2$ (đạt $14.140\text{ km}^2$ nếu tính cả phụ lưu sông Bàn Thạch), trải dài trên địa bàn ba tỉnh Gia Lai, Đắk Lắk và Phú Yên.

                                +--------------------------------------------+
                                |  ÁP LỰC KÉP TRÊN LƯU VỰC SÔNG BA           |
                                |  - Gia tăng dân số, thủy điện, tưới tiêu   |
                                |  - Biến đổi khí hậu, nước biển dâng        |
                                +---------------------+----------------------+
                                                      |
                                                      v
                                +--------------------------------------------+
                                |  CHUỖI MÔ HÌNH TOÁN TÍCH HỢP               |
                                |  NAM -> CROPWAT -> MIKE BASIN -> MIKE 11   |
                                +---------------------+----------------------+
                                                      |
                                                      v
                                +--------------------------------------------+
                                |  ĐÁNH GIÁ TỔN THƯƠNG (VI) & CÂN BẰNG NƯỚC  |
                                |  - 8 nhóm chỉ số áp lực, tài nguyên & QL   |
                                |  - Mô phỏng cân bằng cung - cầu 2020       |
                                +---------------------+----------------------+
                                                      |
                                                      v
                                +--------------------------------------------+
                                |  QUẢN LÝ BỀN VỮNG & QUY HOẠCH PHÂN BỔ      |
                                |  - Khung thể chế Hội đồng Lưu vực sông     |
                                |  - Quy trình vận hành liên hồ chứa         |
                                +--------------------------------------------+

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực trạng các công trình trước đây (tiêu biểu như nghiên cứu của Nguyễn Văn Cư, 2006; Ngô Đình Tuấn, 2006; Huỳnh Thị Lan Hương, 2012) chủ yếu tiếp cận theo góc độ cân bằng nước tĩnh hoặc giải quyết tranh chấp cục bộ đơn ngành, chưa tích hợp đồng bộ các kịch bản BĐKH vào toàn bộ chu trình lập quy hoạch phân bổ tài nguyên nước (TNN) và vận hành liên hồ chứa đa mục tiêu. Nhằm lấp đầy khoảng trống này, luận án xác lập ba câu hỏi nghiên cứu (Research Questions - RQ) và giả thuyết khoa học (Hypotheses - H) tương ứng:

  • RQ1: Động thái dòng chảy và cán cân cung - cầu nguồn nước trên LVS Ba biến động như thế nào dưới các kịch bản BĐKH và phát triển kinh tế - xã hội (KTXH) đến năm 2020? Giả thuyết H1: Tác động tích hợp giữa gia tăng nhu cầu tưới tiêu, chuyển nước thủy điện ngoại lưu vực và BĐKH sẽ làm trầm trọng thêm tình trạng suy kiệt dòng chảy mùa khô, gia tăng nguy cơ xâm nhập mặn (XNM) tại hạ lưu đồng bằng Tuy Hòa.
  • RQ2: Mức độ dễ bị tổn thương của TNN LVS Ba phân bố ra sao theo không gian và thời gian khi xem xét các nhân tố tự nhiên, xã hội và năng lực quản lý? Giả thuyết H2: Tính dễ bị tổn thương TNN không chỉ phụ thuộc vào lượng nước sẵn có mà chịu sự chi phối quyết định bởi hệ số phân bổ bất hợp lý và năng lực thể chế quản trị lưu vực chưa hoàn thiện.
  • RQ3: Khung thể chế và quy trình kỹ thuật nào cho phép tích hợp hiệu quả kịch bản BĐKH vào công tác quản lý tổng hợp tài nguyên nước (QLTHTNN) trên phạm vi toàn lưu vực? Giả thuyết H3: Việc thiết lập quy trình tích hợp BĐKH 5 bước kết hợp với cơ cấu Hội đồng LVS Ba liên tỉnh sẽ tối ưu hóa phân bổ tài nguyên và bảo đảm dòng chảy môi trường.

Khung lý thuyết của luận án được định vị trên nền tảng Quản lý Tổng hợp Tài nguyên Nước theo Định chế Đối tác Nước Toàn cầu (GWP, 2000), 4 Nguyên tắc Dublin (1992), Chương trình Nghị sự 21 (Agenda 21) của Liên Hợp Quốc, kết hợp cùng Khung đánh giá tổn thương khí hậu của Ủy ban Liên chính phủ về Biến đổi Khí hậu (IPCC, 2007). Phạm vi nghiên cứu bao quát toàn bộ diện tích $13.417\text{ km}^2$, tập trung vào nguồn nước mặt và nước mưa trên chuỗi số liệu quan trắc khí tượng thủy văn (KTTV) 1977–2011, đánh giá tiềm năng nước dưới đất, và dự phóng kịch bản cân bằng nước đến năm 2020 với mật độ dân số khu vực dao động từ $51\text{ người/km}^2$ (Tây Nguyên) đến $476\text{ người/km}^2$ (đồng bằng Tuy Hòa).

Literature Review và Positioning

Tổng quan y văn quốc tế và trong nước về quản trị lưu vực sông phản ánh ba dòng lý thuyết và phương pháp tiếp cận chủ đạo:

Thứ nhất, Trường phái Quản trị Nguồn nước Tổng hợp (IWRM Theory) xuất phát từ Hội nghị Quốc tế về Nước và Môi trường tại Dublin (1992) và Chương 18 - Agenda 21, xác định: "Quản lý tổng hợp tài nguyên nước dựa trên nhận thức là một bộ phận nội tại của hệ sinh thái, một nguồn tài nguyên thiên nhiên và một loại hàng hóa kinh tế và xã hội, mà số lượng và chất lượng quyết định bản chất của việc sử dụng." Tổ chức GWP (2000) tiếp tục chuẩn hóa định nghĩa: "Quản lý tổng hợp tài nguyên nước là một quá trình đẩy mạnh, phối hợp phát triển và quản lý nguồn nước, đất đai và tài nguyên liên quan, để tối đa hoá lợi ích kinh tế và phúc lợi xã hội một cách công bằng mà không phương hại đến tính bền vững của các hệ sinh thái thiết yếu." Dòng lý thuyết này nhấn mạnh sự phối hợp đa ngành, phân cấp quản lý theo ranh giới thủy văn thay vì ranh giới hành chính.

Thứ hai, Trường phái Mô hình hóa Thủy văn - Hệ thống Kinh tế Nước (Hydro-Economic Systems Modeling) được phát triển bởi các viện nghiên cứu thủy lực hàng đầu thế giới (DHI, FAO), ứng dụng các công cụ mô phỏng động lực dòng chảy (mô hình NAM, MIKE 11) kết hợp với các mô hình cân bằng hệ thống (MIKE BASIN, WEAP). Tại Việt Nam, hướng tiếp cận này được thể hiện qua các công trình của Huỳnh Thị Lan Hương (2012) về hệ thống hỗ trợ kỹ thuật giải quyết tranh chấp TNN LVS Ba, và Vũ Thành Tâm cùng cộng sự (2012) khi sử dụng mô hình WEAP.

Thứ ba, Trường phái Đánh giá Tổn thương và Thích ứng Biến đổi Khí hậu (Climate Vulnerability & Adaptation Framework) dựa trên nền tảng của IPCC (2001, 2007) và Sullivan (2002) với Chỉ số Nghèo đói Nguồn nước (Water Poverty Index - WPI).

+----------------------------------------------------------------------------------------------------+
|                                    BẢNG TỔNG HỢP TIẾP CẬN Y VĂN                                    |
+--------------------------+------------------------------+------------------------------------------+
| Trường phái / Tác giả    | Trọng tâm nghiên cứu         | Hạn chế / Điểm nghẽn lý thuyết           |
+--------------------------+------------------------------+------------------------------------------+
| IWRM (GWP, 2000;         | Nguyên tắc quản trị tích     | Mang tính nguyên tắc chung, thiếu công cụ|
| Dublin, 1992)            | hợp, đa ngành, đa bên        | định lượng thích ứng BĐKH ở quy mô lưu vực|
+--------------------------+------------------------------+------------------------------------------+
| Thủy văn - Kinh tế       | Mô phỏng cân bằng nước hệ    | Thường giả định điều kiện biên tĩnh,     |
| (DHI; Vũ Thành Tâm, 2012)| thống bằng công cụ mô hình   | chưa lượng hóa tính dễ bị tổn thương     |
+--------------------------+------------------------------+------------------------------------------+
| Tổn thương & Thích ứng   | Đánh giá chỉ số nhạy cảm và  | Tách rời khỏi các công cụ điều tiết      |
| (IPCC, 2007; IMHEN, 2010)| thích ứng khí hậu vĩ mô      | thủy lực và quy hoạch phân bổ liên hồ    |
+--------------------------+------------------------------+------------------------------------------+
| Luận án Lê Đức Thường    | Hợp nhất mô hình toán chuỗi  | Giải quyết triệt để bài toán tích hợp    |
| (2015)                   | với chỉ số tổn thương VI     | BĐKH trong QLTHTNN liên tỉnh             |
+--------------------------+------------------------------+------------------------------------------+

Trong y văn tồn tại hai luồng quan điểm tranh biện sâu sắc: một bên ủng hộ cơ chế thị trường hóa và phân quyền tài chính tự chủ hoàn toàn như mô hình 6 Cục Lưu vực (Agences de l'Eau) của Pháp (ban hành từ Luật Nước 1964, sửa đổi 1992), sử dụng triệt để công cụ thuế tài nguyên và phí xả thải; bên kia là quan điểm duy trì vai trò điều phối tập quyền vĩ mô của nhà nước có tham vấn cộng đồng như Hội đồng Lưu vực Lerma - Chapala tại Mexico (thành lập năm 1993, điều phối bởi Ủy ban Nước Quốc gia CNA) hoặc Ủy ban Lưu vực sông Pungwe giữa Mozambique và Zimbabwe (2002).

Công trình của Lê Đức Thường định vị chính xác tại điểm giao thoa giữa mô hình hóa kỹ thuật phức tạp và thể chế quản trị thích ứng. Nghiên cứu nâng tầm lĩnh vực bằng cách chuyển đổi từ quy hoạch thủy lợi - thủy điện đơn mục tiêu truyền thống sang khung QLTHTNN đa mục tiêu có xét đến BĐKH, vượt trội hơn các công trình nghiên cứu trước đây tại miền Trung Việt Nam khi giải quyết đồng bộ cả ba mắt xích: tính toán thủy văn - đánh giá tổn thương - kiến tạo thể chế.

Đó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 Lý thuyết Quản lý Tổng hợp Tài nguyên Nước (IWRM) trong bối cảnh các lưu vực nhiệt đới gió mùa có độ tương phản địa hình và khí hậu sâu sắc. Nghiên cứu chỉ ra rằng việc áp dụng cơ học 4 Nguyên tắc Dublin tại các quốc gia đang phát triển sẽ gặp thất bại nếu không lượng hóa được tham số tổn thương do khí hậu và cấu trúc phân bổ nguồn nước xuyên địa giới hành chính.

Cụ thể, luận án đóng góp vào lý thuyết quản trị tài nguyên thông qua việc xây dựng Mô hình Chỉ số Tổn thương Tài nguyên Nước ($VI$ - Vulnerability Index) bao gồm 8 phân hệ chỉ thị thành phần:

  1. Hệ số khan hiếm nguồn nước ($VI_1$)
  2. Hệ số sức ép nguồn nước ($VI_2$)
  3. Hệ số tiếp nhận nguồn nước sạch ($VI_3$)
  4. Hệ số ô nhiễm nguồn nước ($VI_4$)
  5. Tổng sản phẩm quốc nội bình quân ($VI_5$ - GDP/người/năm)
  6. Hệ số sử dụng hiệu quả nguồn nước ($VI_6$)
  7. Hệ số khả năng tiếp nhận vệ sinh môi trường ($VI_7$)
  8. Hệ số đánh giá năng lực quản lý lưu vực ($VI_8$)
             +------------------------------------------------------------+
             |            CHỈ SỐ DỄ BỊ TỔN THƯƠNG TNN (VI)                |
             +-----------------------------+------------------------------+
                                           |
                   +-----------------------+-----------------------+
                   |                                               |
                   v                                               v
    +------------------------------+               +------------------------------+
    |  ÁP LỰC VÀ HIỆN TRẠNG TNN    |               |  NĂNG LỰC THÍCH ỨNG & KTXH   |
    |  - Khan hiếm nguồn nước (VI1)|               |  - GDP bình quân (VI5)       |
    |  - Sức ép khai thác (VI2)    |               |  - Hiệu quả dùng nước (VI6)  |
    |  - Tiếp cận nước sạch (VI3)  |               |  - Vệ sinh môi trường (VI7)  |
    |  - Ô nhiễm nguồn nước (VI4)  |               |  - Năng lực quản lý (VI8)    |
    +------------------------------+               +------------------------------+

Nghiên cứu xác lập các mệnh đề lý thuyết cốt lõi:

  • Mệnh đề 1: Tác động của BĐKH đến an ninh nguồn nước tại lưu vực sông dốc ven biển có tính chất phi tuyến, trong đó sự gia tăng nhiệt độ và bốc hơi tiềm năng ($ET_0$) tại trung lưu đóng vai trò kích hoạt khủng hoảng thiếu nước mùa cạn mạnh hơn mức suy giảm lượng mưa cục bộ.
  • Mệnh đề 2: Năng lực thích ứng và cơ chế chia sẻ dòng chảy giữa thượng lưu (thủy điện, nông nghiệp cao nguyên) và hạ lưu (đô thị, nông nghiệp đồng bằng) quyết định tính bền vững của lưu vực hơn là tổng lượng dòng chảy năm.

Khám phá này tạo nên bước chuyển dịch mô hình (paradigm shift) từ "Tư duy khai thác công trình cục bộ" (Supply-side Engineering Path) sang "Quản trị thích ứng hệ thống mềm tích hợp" (Adaptive Soft-path Governance).

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

Khung phân tích của luận án tích hợp ba trụ cột lý thuyết: Khung IWRM của GWP, Khung đánh giá tổn thương thích ứng của IPCC, và Lý thuyết cân bằng hệ thống thủy văn - kinh tế. Cách tiếp cận phân tích mới được cụ thể hóa bằng Quy trình 5 bước tích hợp vấn đề BĐKH vào QLTHTNN LVS Ba:

  • Bước 1: Đánh giá tổng thể hiện trạng TNN, hiện trạng KTXH và hiện trạng khai thác sử dụng nước.
  • Bước 2: Đánh giá tác động của BĐKH đến chế độ dòng chảy, ngập lụt, nhu cầu dùng nước và xâm nhập mặn.
  • Bước 3: Tính toán cân bằng nước hệ thống lưu vực theo các kịch bản hiện trạng và tương lai có xét đến BĐKH.
  • Bước 4: Đánh giá tính dễ bị tổn thương của TNN thông qua bộ chỉ số $VI$.
  • Bước 5: Xây dựng giải pháp thích ứng, lồng ghép vào quy hoạch phân bổ nguồn nước và quy trình vận hành liên hồ chứa.

Điều kiện biên (boundary conditions) của khung phân tích được xác lập trong giới hạn địa lý lưu vực sông Ba thuộc 3 tỉnh Gia Lai, Đắk Lắk, Phú Yên; chuỗi số liệu thủy văn chuẩn hóa giai đoạn 1977–2011; kịch bản BĐKH phát thải trung bình và cao của Bộ Tài nguyên và Môi trường.

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

Thiết kế nghiên cứu

Luận án áp dụng thế giới quan thực chứng hậu hiện đại (Post-positivism) kết hợp với chủ nghĩa hiện thực phản biện (Critical Realism) và tư duy lý thuyết hệ thống (Systems Theory). Thiết kế nghiên cứu hỗn hợp (Mixed-methods design) được triển khai thông qua việc kết hợp các mô hình toán học thủy văn - thủy lực - cân bằng nước phân bố không gian với phương pháp điều tra thực địa xã hội học và tham vấn chuyên gia (Expert elicitation).

+----------------------------------------------------------------------------------------------------+
|                                    CHUỖI CÔNG CỤ MÔ HÌNH HÓA                                       |
+-------------------+--------------------+------------------------------------+----------------------+
| Mô hình toán      | Cơ quan phát triển | Phạm vi ứng dụng trong luận án     | Chỉ số kiểm định     |
+-------------------+--------------------+------------------------------------+----------------------+
| NAM (DHI)         | Viện Thủy lực      | Mô phỏng mưa - dòng chảy, khôi     | Nash-Sutcliffe       |
|                   | Đan Mạch           | phục số liệu dòng chảy toàn lưu vực| $R^2 > 0,75 - 0,85$  |
+-------------------+--------------------+------------------------------------+----------------------+
| CROPWAT           | Tổ chức Nông Lương | Tính toán nhu cầu tưới cây trồng   | Hiệu chỉnh theo chuỗi|
|                   | LHQ (FAO)          | và biến động nhu cầu nước do BĐKH  | khí tượng 1977-2011  |
+-------------------+--------------------+------------------------------------+----------------------+
| MIKE BASIN (DHI)  | Viện Thủy lực      | Mô phỏng mạng lưới cân bằng nước   | Cân bằng chuỗi thời  |
|                   | Đan Mạch           | hệ thống cho 3 giai đoạn           | đoạn 10 ngày / tháng |
+-------------------+--------------------+------------------------------------+----------------------+
| MIKE 11 (DHI)     | Viện Thủy lực      | Thủy lực 1D lan truyền lũ và xâm   | Kiểm định trạm       |
|                   | Đan Mạch           | nhập mặn (XNM) tại hạ lưu Đà Rằng  | Phú Lâm (1988, 1993) |
+-------------------+--------------------+------------------------------------+----------------------+

Thiết kế nghiên cứu đa cấp độ (Multi-level design) chia lưu vực sông Ba thành 3 khu vực chính: Thượng lưu (Gia Lai), Trung lưu (Đắk Lắk - Gia Lai) và Hạ lưu (Phú Yên) với các tiểu lưu vực trọng điểm: Sông chính Ba, Ia Yun, Krông Hnăng, Sông Hinh.

Quy trình nghiên cứu rigorous

Chiến lược thu thập và xử lý dữ liệu được thiết kế chặt chẽ:

  • Thu thập dữ liệu khí tượng: 10 trạm khí tượng (5 trạm nội vùng: An Khê, Ayun Pa, Sơn Hòa, M'Đrăk, Tuy Hòa và 5 trạm lân cận: Quy Nhơn, Hoài Nhơn, Pleiku, Buôn Hồ, Buôn Ma Thuột). Chuỗi số liệu bức xạ, số giờ nắng ($2200–2500\text{ giờ/năm}$), nhiệt độ, lượng bốc hơi tiềm năng và lượng mưa giai đoạn 1977–2011.
  • Thu thập dữ liệu thủy văn: 7 trạm thủy văn quan trắc dòng chảy và mực nước (An Khê, Tuy Hòa/Phú Lâm, Ayun Pa, Pơ Mơ Rê, Củng Sơn, Sông Hinh, Krông Hnăng).
  • Phương pháp tam giác đạc (Triangulation): Đối chiếu chéo giữa chuỗi số liệu quan trắc thực đo lịch sử, kết quả mô phỏng số trị từ chuỗi mô hình toán học (NAM, CROPWAT, MIKE BASIN, MIKE 11), và dữ liệu điều tra thực địa hiện trạng kinh tế - xã hội, công trình thủy lợi - thủy điện.
  • Độ tin cậy và kiểm định mô hình: Các mô hình được hiệu chỉnh và kiểm định độc lập với các trận lũ lịch sử (tháng X/1988, X/1993, XI/2005) và các chuỗi mùa kiệt tại trạm Củng Sơn và Phú Lâm.

Data và phân tích

Đặc trưng dữ liệu thể hiện các thông số thủy văn - khí hậu chi tiết:

  • Lượng mưa bình quân năm phân hóa mạnh: từ $1.297\text{ mm}$ (trạm Ayun Pa) lên đến $2.083\text{ mm}$ (Đăk Đoa), $2.217\text{ mm}$ (Pleiku), $2.400\text{ mm}$ (Sông Hinh).
  • Lượng bốc hơi tiềm năng trung bình năm lớn nhất đạt $1.625\text{ mm/năm}$ tại trạm trung lưu Ayun Pa, trong khi thượng nguồn An Khê chỉ $1.420\text{ mm/năm}$ và M'Đrăk đạt $1.434\text{ mm/năm}$.
  • Nhiệt độ trung bình năm dao động $21,9–26,4^\circ\text{C}$; nhiệt độ cao nhất tuyệt đối đạt $41,7^\circ\text{C}$ (ghi nhận tại Sơn Hòa, tháng IV/1983) và thấp nhất tuyệt đối đạt $9^\circ\text{C}$ (An Khê, tháng I/1984).
  • Về thủy văn: Tổng lượng dòng chảy năm của lưu vực đạt xấp xỉ 10 tỷ $\text{m}^3$, tương ứng module dòng chảy bình quân toàn lưu vực $M_0 = 22,8\text{ l/s}\cdot\text{km}^2$. Tuy nhiên, có sự phân hóa cực đoan giữa các phụ lưu: lưu vực sông Hinh đạt $M_0 = 53,4–58,94\text{ l/s}\cdot\text{km}^2$, sông Krông Hnăng đạt $M_0 = 21,7\text{ l/s}\cdot\text{km}^2$, trong khi sông Ia Yun chỉ đạt $M_0 = 18,9\text{ l/s}\cdot\text{km}^2$.

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

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

  1. Sự bất cân đối không - thời gian cực đoan của tài nguyên nước: Mặc dù tổng lượng dòng chảy hàng năm xấp xỉ 10 tỷ $\text{m}^3$, nhưng từ 70% đến 80% lượng dòng chảy tập trung trong mùa lũ ngắn ngủi (tháng IX–XII ở hạ lưu, tháng V–X ở thượng lưu), trong khi mùa cạn kéo dài 6-8 tháng chỉ chiếm 20-30% tổng lượng nước. Tại khu vực lòng chảo Ayun Pa, lượng mưa chỉ $1.297\text{ mm}$ nhưng bốc hơi tiềm năng lên tới $1.625\text{ mm/năm}$, tạo nên túi hạn hán khốc liệt ngay tại trung lưu.
  2. Tác động kép của BĐKH và quy hoạch thủy điện liên vùng: Kết quả mô phỏng cân bằng nước MIKE BASIN và thủy lực MIKE 11 chỉ ra rằng dưới kịch bản BĐKH và phương án dùng nước phát triển 2020, dòng chảy mùa cạn tại trạm Củng Sơn suy giảm từ 15% đến 22%. Đồng thời, việc chuyển nước của cụm thủy điện An Khê - Ka Nak sang lưu vực sông Kôn làm suy kiệt dòng chảy hạ lưu sông Ba, đẩy ranh mặn 4‰ tại cửa sông Đà Rằng lấn sâu vào nội đồng thêm 2,5 đến 4,2 km, đe dọa trực tiếp các trạm bơm cấp nước sinh hoạt và tưới tiêu của đồng bằng Tuy Hòa.
  3. Nhu cầu nước nông nghiệp chiếm tỷ trọng chi phối tuyệt đối: Nhu cầu sử dụng nước tưới nông nghiệp chiếm trên 80% tổng nhu cầu nước trên toàn lưu vực. Dưới tác động của BĐKH làm gia tăng nhiệt độ và bốc hơi, nhu cầu nước tưới cho cây trồng (lúa, mía, cà phê, cao su) tính toán từ mô hình CROPWAT tăng thêm từ 7,5% đến 12,3% vào các tháng mùa khô (tháng II đến tháng V).
  4. Phân hóa mức độ tổn thương nguồn nước qua chỉ số $VI$: Kết quả lượng hóa chỉ số tổn thương cho thấy các tiểu lưu vực Ayun Pa và hạ lưu Tuy Hòa có chỉ số $VI$ thuộc nhóm cao nhất toàn lưu vực. Tình trạng dễ tổn thương này không đơn thuần do thiếu hụt tự nhiên mà xuất phát từ mật độ dân số cao ($476\text{ người/km}^2$ ở Tuy Hòa), tỷ lệ dùng nước nông nghiệp lớn và sự thiếu vắng cơ chế vận hành liên hồ chứa chia sẻ dòng chảy tối thiểu.
                  BIẾN THIÊN DÒNG CHẢY MÙA KIỆT VÀ XÂM NHẬP MẶN
      +-------------------------------------------------------------------+
      | Tình trạng tự nhiên / Lịch sử                                    |
      | Module dòng chảy: M0 = 22.8 l/s.km2                               |
      | Ranh mặn 4‰ duy trì ngoài cửa Đà Rằng                             |
      +---------------------------------+---------------------------------+
                                        |  Tác động BĐKH & Thủy điện 2020
                                        v
      +-------------------------------------------------------------------+
      | Kịch bản 2020 (BĐKH + An Khê - Ka Nak chuyển dòng)                |
      | Dòng chảy mùa cạn Củng Sơn suy giảm: -15% đến -22%                |
      | Ranh mặn 4‰ lấn sâu vào nội đồng: +2.5 km đến +4.2 km             |
      +-------------------------------------------------------------------+

Implications đa chiều

  • Về mặt lý thuyết: Bổ sung vào hệ thống lý luận IWRM bộ chỉ số $VI$ 8 thành phần và quy trình tích hợp BĐKH chuẩn hóa dành riêng cho lưu vực sông dốc ven biển nhiệt đới.
  • Về mặt phương pháp luận: Thiết lập chuỗi công cụ mô hình hóa nối tiếp (NAM $\rightarrow$ CROPWAT $\rightarrow$ MIKE BASIN $\rightarrow$ MIKE 11) có khả năng chuyển giao và áp dụng rộng rãi cho các lưu vực sông miền Trung Việt Nam có đặc điểm tương đồng như LVS Vu Gia - Thu Bồn, LVS Trà Khúc, LVS Kôn.
  • Về mặt thực tiễn và chính sách:
    1. Kiến nghị tái cấu trúc Quy trình vận hành liên hồ chứa trên LVS Ba (gồm các hồ Sông Ba Hạ, Sông Hinh, Krông Hnăng, Ayun Hạ, An Khê - Ka Nak), bắt buộc duy trì dòng chảy môi trường tối thiểu sau đập trong mùa kiệt nhằm ngăn mặn tại cửa Đà Rằng.
    2. Thành lập Ủy ban/Hội đồng Lưu vực sông Ba với sự tham gia bình đẳng của đại diện 3 tỉnh Gia Lai, Đắk Lắk, Phú Yên cùng các bên sử dụng nước lớn (Tập đoàn Điện lực Việt Nam - EVN, các công ty khai thác công trình thủy lợi, hiệp hội người dùng nước).
    3. Áp dụng các công cụ kinh tế như chi trả dịch vụ môi trường rừng (PES), thu phí cấp quyền khai thác nước mặt và phí xả thải để tạo quỹ tài chính tự chủ bảo vệ nguồn nước lưu vực.

Limitations và Future Research

Tác giả luận án thẳng thắn chỉ ra 4 giới hạn nghiên cứu chính:

  1. Dữ liệu nước dưới đất: Mạng lưới quan trắc địa chất thủy văn trên lưu vực còn rất mỏng và gián đoạn, do đó nghiên cứu mới dừng lại ở mức đánh giá tiềm năng và trữ lượng khai thác tĩnh, chưa thể mô phỏng động lực tương tác 3 chiều giữa nước mặt và nước dưới đất.
  2. Độ phân giải kịch bản BĐKH: Các kịch bản BĐKH kế thừa từ Bộ TN&MT (2010) có độ phân giải không gian ở mức khu vực, chưa phản ánh đầy đủ các tiểu vùng vi khí hậu phức tạp giữa sườn đông và tây dãy Trường Sơn.
  3. Động lực thay đổi lớp phủ đất: Nghiên cứu giả định lớp phủ rừng và thảm thực vật biến động theo quy hoạch tĩnh, chưa mô phỏng động lực chuyển đổi đất nông - lâm nghiệp theo thời gian thực dưới áp lực di dân và biến động thị trường nông sản.
  4. Lượng hóa chi phí - lợi ích kinh tế chi tiết: Chưa xây dựng được hàm tối ưu hóa phân bổ kinh tế biên (Hydro-economic optimization) chi tiết cho từng loại cây trồng và từng tổ máy thủy điện.

Chương trình nghiên cứu 10 năm tiếp theo (Future Research Agenda) cần tập trung vào:

  • Tích hợp mô hình dòng ngầm 3D MODFLOW với mô hình bề mặt SWAT/MIKE SHE để đánh giá toàn diện chu trình thủy văn.
  • Lượng hóa kinh tế giá trị dịch vụ hệ sinh thái nước và xây dựng thị trường hạn ngạch nước nội lưu vực.
  • Ứng dụng trí tuệ nhân tạo (Machine Learning) và mạng lưới cảm biến IoT trong vận hành thời gian thực (real-time operation) hệ thống liên hồ chứa thích ứng với các hiện tượng thời tiết cực đoan.

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

Nghiên cứu tạo ra tác động học thuật sâu rộng với dự phóng trích dẫn cao trong các tài liệu chuyên khảo về quản lý lưu vực sông tại Việt Nam và Đông Nam Á. Về phương diện chuyển đổi ngành, các phát hiện của luận án đã cung cấp luận cứ kỹ thuật trực tiếp cho Tập đoàn Điện lực Việt Nam (EVN) và Bộ Công Thương trong việc tái lập lịch xả nước của các nhà máy thủy điện trên sông Ba, đặc biệt là việc điều tiết xả dòng chảy môi trường từ hồ An Khê trả lại dòng chính sông Ba để cứu hạn cho hạ du Gia Lai và Phú Yên.

+----------------------------------------------------------------------------------------------------+
|                                BẢNG ĐỐI TƯỢNG HƯỞNG LỢI VÀ GIÁ TRỊ                                 |
+--------------------------+-------------------------------------------------------------------------+
| Nhóm đối tượng           | Giá trị khoa học và thực tiễn tiếp nhận                                 |
+--------------------------+-------------------------------------------------------------------------+
| Nghiên cứu sinh &        | Kho tư liệu tính toán KTTV 1977-2011; phương pháp luận chuỗi mô hình    |
| Giới học thuật           | NAM - CROPWAT - MIKE BASIN - MIKE 11 và bộ chỉ số tổn thương VI        |
+--------------------------+-------------------------------------------------------------------------+
| Cơ quan hoạch định       | Luận cứ khoa học ban hành Quy hoạch tài nguyên nước LVS Ba; quy trình   |
| chính sách (Bộ TN&MT,    | vận hành liên hồ chứa đa mục tiêu; đề án thành lập Hội đồng Lưu vực     |
| UBND 3 tỉnh)             |                                                                         |
+--------------------------+-------------------------------------------------------------------------+
| Đơn vị R&D & Doanh       | Tối ưu hóa điều tiết phát điện gắn với trách nhiệm cấp nước hạ du; giải |
| nghiệp (EVN, Thủy nông)  | pháp kiểm soát mặn cho các công trình lấy nước vùng cửa sông            |
+--------------------------+-------------------------------------------------------------------------+
| Cộng đồng dân cư         | Giảm thiểu rủi ro thiếu nước ngọt sinh hoạt và tưới tiêu nông nghiệp;   |
| (1,89 triệu người)       | bảo đảm an ninh lương thực và sinh kế bền vững trước BĐKH               |
+--------------------------+-------------------------------------------------------------------------+

Ở cấp độ chính sách, công trình đóng góp trực tiếp vào quá trình hiện thực hóa Luật Tài nguyên nước 2012 (có hiệu lực 2013), Nghị định 120/2008/NĐ-CP về quản lý lưu vực sông và Nghị định 112/2008/NĐ-CP về quản lý khai thác tổng hợp hồ chứa thủy điện, thủy lợi. Về mặt xã hội, kết quả phân bổ nước tối ưu bảo vệ nguồn sinh kế cho hơn 1,89 triệu người dân (trong đó 44,4% là đồng bào các dân tộc thiểu số Gia Rai, Ba Na, Xê Đăng) sinh sống trên toàn lưu vực.

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

  • Nghiên cứu sinh và Giới nghiên cứu Thủy văn: Tiếp cận hệ thống cơ sở dữ liệu KTTV 1977–2011 chi tiết, bộ thông số mô hình NAM hiệu chỉnh cho các lưu vực sông dốc miền Trung, và phương pháp luận thiết lập chỉ số $VI$.
  • Các nhà Khoa học Cấp cao: Tham khảo mô hình khung phân tích tích hợp BĐKH đa tầng, làm tiền đề mở rộng nghiên cứu sang các lưu vực sông quốc tế lân cận.
  • Khối Doanh nghiệp & Kỹ sư Thủy điện - Thủy lợi: Ứng dụng quy trình tính toán cân bằng nước và vận hành liên hồ chứa để nâng cao hiệu suất phát điện đồng thời giảm thiểu xung đột pháp lý với hạ du.
  • Các Nhà Hoạch định Chính sách & Lãnh đạo Địa phương: Sở TN&MT, Sở NN&PTNT các tỉnh Gia Lai, Đắk Lắk, Phú Yên sở hữu bộ công cụ khoa học chính xác để đàm phán phân bổ chỉ ngạch khai thác nguồn nước và cấp phép xả thải.

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ì? Đó là việc mở rộng Lý thuyết Quản lý Tổng hợp Tài nguyên Nước (IWRM) bằng việc xây dựng Khung định lượng hóa Chỉ số Dễ bị Tổn thương Tài nguyên Nước ($VI$) tích hợp 8 nhóm chỉ thị đa ngành (từ tự nhiên, thủy văn đến kinh tế, xã hội và thể chế quản lý), khắc phục triệt để tính định tính chung chung của các nguyên tắc IWRM truyền thống khi áp dụng tại các lưu vực sông dốc có áp lực chuyển nước ngoại lưu vực.

2. Điểm cải tiến phương pháp luận so với các nghiên cứu trước đây là gì? So với nghiên cứu của Nguyễn Văn Cư (2006) và Ngô Đình Tuấn (2006) vốn chỉ tiếp cận cân bằng nước tĩnh đơn ngành, luận án đã thiết lập một chuỗi mô hình số trị liên hoàn (NAM $\rightarrow$ CROPWAT $\rightarrow$ MIKE BASIN $\rightarrow$ MIKE 11), kết hợp kiểm định thủy lực 1D lan truyền lũ và xâm nhập mặn với kịch bản BĐKH động giai đoạn đến 2020.

3. Phát hiện bất ngờ nhất từ số liệu tính toán thực nghiệm là gì? Module dòng chảy năm của sông Hinh ($M_0 = 53,4–58,94\text{ l/s}\cdot\text{km}^2$) gấp gần 3 lần module bình quân lưu vực sông Ba ($22,8\text{ l/s}\cdot\text{km}^2$) và gấp hơn 3 lần sông Ia Yun ($18,9\text{ l/s}\cdot\text{km}^2$). Đồng thời, trạm trung lưu Ayun Pa có lượng bốc hơi tiềm năng lên tới $1.625\text{ mm/năm}$ vượt xa lượng mưa tại chỗ ($1.297\text{ mm/năm}$), chứng minh hiện tượng "nghịch lý thủy văn" khi vùng trung lưu lại là tâm điểm thiếu nước và có chỉ số tổn thương $VI$ cao nhất chứ không phải vùng thượng nguồn.

            SO SÁNH MODULE DÒNG CHẢY (M0) GIỮA CÁC TIỂU LƯU VỰC
   +-----------------------------------------------------------------------+
   | Tiểu lưu vực Sông Hinh   | [================================] 58.94   |
   | Toàn lưu vực Sông Ba     | [============] 22.80                       |
   | Phụ lưu Krông Hnăng      | [===========] 21.70                        |
   | Phụ lưu Ia Yun           | [==========] 18.90                         |
   +-----------------------------------------------------------------------+
     Đơn vị: l/s.km2

4. Luận án có cung cấp đầy đủ giao thức tái lập nghiên cứu (Replication Protocol) không? Có. Toàn bộ bảng thông số hiệu chuẩn mô hình thủy văn NAM (hệ số trữ ẩm bề mặt $U_{max}$, trữ ẩm tầng rễ $L_{max}$, hệ số thấm qua tầng ngầm $CQOF$, thông số dòng chảy ngầm $CKBF$), sơ đồ mạng lưới nút cân bằng MIKE BASIN cho 3 giai đoạn (1980–1999, 2000–2011, 2020) và hệ thống phương trình tính toán 8 chỉ số thành phần $VI$ đều được chuẩn hóa chi tiết trong luận án.

5. Luận án đã phác thảo chương trình nghiên cứu 10 năm tiếp theo như thế nào? Chương trình nghiên cứu tiếp theo bao gồm việc tích hợp các mô hình thủy văn phân bố hoàn toàn (Fully-distributed hydrological models), xây dựng mô hình kinh tế thủy văn tối ưu hóa phi tuyến đa mục tiêu, và phát triển hệ sinh thái quản trị dữ liệu lớn (Big Data) thời gian thực phục vụ vận hành hồ chứa thích ứng khí hậu.

Kết luận

Công trình nghiên cứu tiến sĩ của Lê Đức Thường tạo nên dấu ấn học thuật sâu sắc với 5 đóng góp cốt lõi:

  1. Xác lập cơ sở khoa học toàn diện và có hệ thống nhất từ trước đến nay về hiện trạng, tiềm năng và quy luật biến động của tài nguyên nước LVS Ba trên chuỗi số liệu quan trắc dài hạn 1977–2011.
  2. Thiết lập chuỗi công cụ mô hình hóa số trị tích hợp (NAM, CROPWAT, MIKE BASIN, MIKE 11) cho phép dự báo định lượng dòng chảy, ngập lụt, nhu cầu tưới và xâm nhập mặn dưới tác động của BĐKH đến năm 2020.
  3. Phát triển Khung đánh giá và lượng hóa Chỉ số Dễ bị Tổn thương Tài nguyên Nước ($VI$) đa tiêu chí, định vị chính xác các điểm nóng khủng hoảng nguồn nước theo không gian lưu vực.
  4. Xây dựng Quy trình 5 bước tích hợp vấn đề BĐKH vào QLTHTNN lưu vực sông, cung cấp công cụ chuyển giao phương pháp luận cho công tác quy hoạch tài nguyên nước tại Việt Nam.
  5. Đề xuất hệ thống giải pháp công trình và phi công trình mang tính khả thi cao, trọng tâm là khung vận hành liên hồ chứa bảo đảm dòng chảy môi trường và mô hình thể chế Hội đồng Lưu vực sông Ba liên tỉnh.
+----------------------------------------------------------------------------------------------------+
|                                    TỔNG KẾT ĐÓNG GÓP LUẬN ÁN                                       |
|  1. Hệ thống hóa cơ sở dữ liệu KTTV 1977-2011 toàn diện lưu vực sông Ba ($13.417\text{ km}^2$)       |
|  2. Chuỗi mô hình hóa liên hoàn: NAM -> CROPWAT -> MIKE BASIN -> MIKE 11                            |
|  3. Bộ chỉ số tổn thương VI 8 nhóm tham số định lượng hóa rủi ro an ninh nguồn nước                |
|  4. Quy trình chuẩn 5 bước tích hợp BĐKH vào quy hoạch quản lý tổng hợp nước lưu vực              |
|  5. Luận cứ thể chế thành lập Hội đồng LVS Ba và vận hành liên hồ chứa bảo đảm dòng chảy sinh thái |
+----------------------------------------------------------------------------------------------------+

Công trình thúc đẩy bước tiến căn bản trong mô hình lý thuyết quản trị nước tại Việt Nam, mở ra ba hướng nghiên cứu chuyên sâu về mô hình hóa liên kết nước mặt - nước ngầm, thị trường hóa quyền sử dụng nước và quản trị lưu vực sông thông minh thích ứng với biến đổi toàn cầu.