Tổng quan về luận án

Nghiên cứu "Nghiên cứu xây dựng phương pháp xác định dòng chảy tối thiểu – Áp dụng cho hệ thống sông Vu Gia – Thu Bồn" của tác giả Nguyễn Thị Kim Dung (Viện Khoa học Thủy lợi Việt Nam, 2018) thuộc chuyên ngành Kỹ thuật Tài nguyên nước (Mã số: 62 58 02 12), dưới sự hướng dẫn khoa học của PGS.TS. Nguyễn Vũ Việt và PGS.TS. Nguyễn Quang Trung, là công trình tiên phong giải quyết bài toán xung đột tài nguyên nước phức tạp tại lưu vực duyên hải miền Trung Việt Nam.

                  ┌─────────────────────────────────────────────────────────┐
                  │          KHUNG XÁC ĐỊNH DÒNG CHẢY TỐI THIỂU (Q_TT)       │
                  └────────────────────────────┬────────────────────────────┘
                                               │
               ┌───────────────────────────────┼───────────────────────────────┐
               ▼                               ▼                               ▼
┌─────────────────────────────┐ ┌─────────────────────────────┐ ┌─────────────────────────────┐
│  Dòng chảy duy trì sông     │ │  Dòng chảy sinh thái        │ │  Dòng chảy khai thác sử dụng│
│           (Q_dts)           │ │           (Q_st)            │ │          (Q_ktsd)           │
├─────────────────────────────┤ ├─────────────────────────────┤ ├─────────────────────────────┤
│• Phân tích thủy văn         │ │• Phương pháp chu vi ướt     │ │• Mô hình MIKE BASIN         │
│• Đường cong duy trì FDC     │ │  (Wetted Perimeter)         │ │  (Cân bằng nước hệ thống)   │
│• Chỉ số Q90, Q95, 7Q10      │ │• Mô hình MIKE 11 Ecolab     │ │• Mô hình MIKE 11 HD & AD    │
│• Duy trì hình thái lòng dẫn │ │  (DO, BOD5, NH4+, NO3-)     │ │  (Mực nước & đẩy mặn)       │
└──────────────┬──────────────┘ └──────────────┬──────────────┘ └──────────────┬──────────────┘
               │                               │                               │
               └───────────────────────────────┼───────────────────────────────┘
                                               │
                                               ▼
                  ┌─────────────────────────────────────────────────────────┐
                  │    TỔ HỢP ĐIỀU HỢP THEO ĐIỂM KIỂM SOÁT (ĐKS) VÀ THỜI ĐOẠN │
                  │            Q_TT = Max(Q_dts, Q_st, Q_ktsd)              │
                  └─────────────────────────────────────────────────────────┘

Về bối cảnh khoa học, tài nguyên nước tại lưu vực sông Vu Gia – Thu Bồn đóng vai trò quyết định đối với an ninh nguồn nước, kinh tế nông nghiệp và hệ sinh thái vùng hạ du thuộc tỉnh Quảng Nam và thành phố Đà Nẵng. Mặc dù tổng lượng dòng chảy trung bình năm đạt khoảng 20,4 tỷ $m^3$, nhưng sự phân bố dòng chảy theo mùa diễn biến bất đối xứng cực kỳ nghiêm trọng: lượng dòng chảy trong 8 tháng mùa khô chỉ chiếm từ 30% đến 35% tổng lượng cả năm. Nghiên cứu khẳng định: "Việc chuyển nước của thủy điện Đắk Mi 4 từ sông Vu Gia sang sông Thu Bồn khiến dòng chảy mùa kiệt trên sông Vu Gia suy giảm mạnh, mực nước hạ thấp, mặn xâm nhập cao ảnh hưởng lớn đến hoạt động của các nhà máy cấp nước sinh hoạt cho TP. Đà Nẵng, các trạm bơm lấy nước tưới dọc sông."

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tiễn: các quy định pháp lý như Luật Tài nguyên nước 2012, Nghị định 112/2008/NĐ-CP, Nghị định 120/2008/NĐ-CP và Thông tư 64/2017/TT-BTNMT đã quy định bắt buộc phải duy trì dòng chảy tối thiểu (DCTT), đồng thời Quyết định số 1537/QĐ-TTg ngày 07/09/2015 của Thủ tướng Chính phủ đã ban hành 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 theo thời đoạn 10 ngày. Tuy nhiên, quy trình vận hành này chủ yếu dựa trên mực nước khống chế tại hai trạm Ái Nghĩa và Giao Thủy mà chưa có cơ sở khoa học định lượng DCTT chi tiết cho từng đoạn sông đặc thù, chưa tích hợp đồng thời giữa yếu tố duy trì hình thái dòng sông, bảo tồn sinh thái thủy sinh và đáp ứng nhu cầu khai thác sử dụng nước vùng hạ du.

Luận án thiết lập hệ thống câu hỏi nghiên cứu và giả thuyết khoa học:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để phân rã và định lượng cấu trúc DCTT thành các thành phần chức năng độc lập nhưng tương hỗ trên một lưu vực chịu tác động chuyển dòng nhân tạo?
  • Câu hỏi nghiên cứu 2 (RQ2): Mối quan hệ thủy lực giữa chu vi ướt, chất lượng nước và không gian sống của thủy sinh vật thay đổi như thế nào tại các điểm kiểm soát trọng yếu trong mùa kiệt?
  • Câu hỏi nghiên cứu 3 (RQ3): Lưu lượng xả tối thiểu từ các công trình thượng lưu cần đạt ngưỡng bao nhiêu để kiểm soát xâm nhập mặn và đảm bảo mực nước lấy nước cho các trạm bơm hạ du trong các thời kỳ cấp nước gia tăng và thông thường?
  • Giả thuyết nghiên cứu 1 (H1): DCTT không phải là một giá trị lưu lượng tĩnh đơn lẻ mà là một hàm tổ hợp phi tuyến tính $Q_{TT} = \max(Q_{dts}, Q_{st}, Q_{ktsd})$, biến thiên theo không gian (từng đoạn sông) và thời gian (các thời đoạn mùa kiệt).
  • Giả thuyết nghiên cứu 2 (H2): Việc áp dụng phương pháp chu vi ướt kết hợp mô phỏng chất lượng nước (Ecolab) sẽ xác định được ngưỡng dòng chảy sinh thái tới hạn bảo vệ các loài cá bản địa quý hiếm mà các phương pháp thủy văn truyền thống bỏ qua.

Phạm vi nghiên cứu bao quát toàn bộ diện tích lưu vực sông Vu Gia – Thu Bồn ($10.350\text{ km}^2$), tập trung vào hệ thống dòng chính, 756 hồ chứa và đập dâng, mạng lưới bậc thang thủy điện (A Vương, Sông Bung 4, Đắk Mi 4, Sông Tranh 2) và chuỗi số liệu thủy văn - khí tượng thực đo giai đoạn 1977–2016. Đóng góp đột phá của công trình là xây dựng thành công bộ công cụ tích hợp phương pháp thủy văn, thủy lực chu vi ướt và mô hình toán (MIKE BASIN, MIKE 11 HD/AD/Ecolab), cung cấp luận cứ chuẩn xác để tối ưu hóa Quy trình vận hành liên hồ chứa và bảo đảm an ninh nước vùng duyên hải.

Literature Review và Positioning

Tổng quan y văn quốc tế và trong nước về dòng chảy môi trường (Environmental Flows - E-Flows) và DCTT ghi nhận bốn trường phái tiếp cận chủ đạo:

Trường phái tiếp cận Đại diện tiêu biểu & Năm Nguyên lý cốt lõi Ưu điểm Hạn chế chính
Thủy văn (Hydrological) Tennant (1976), Richter et al. (1997) [RVA] Dựa trên phân tích thống kê chuỗi lưu lượng lịch sử (FDC, $Q_{90}$, $Q_{95}$, 7Q10) Nhanh, chi phí thấp, áp dụng diện rộng Bỏ qua đặc trưng hình thái lòng dẫn và động học sinh thái
Thủy lực (Hydraulic) Gippel & Stewardson (1998) Dựa trên quan hệ giữa chu vi ướt ($P_w$), bán kính thủy lực và lưu lượng ($Q$) Phản ánh hình thái mặt cắt và diện tích đáy ướt sinh thái Chưa tính toán định lượng nhu cầu chất lượng nước và sử dụng nước
Mô phỏng nơi cư trú (Habitat) Bovee (1982) [IFIM/PHABSIM] Mô phỏng vi môi trường sống (vận tốc, độ sâu, giá thể) theo loài chỉ thị Độ chính xác sinh học cao Đòi hỏi dữ liệu sinh thái chi tiết, chi phí quan trắc rất lớn
Toàn diện (Holistic) King & Louw (1998) [BBM], Brown & King (2000) [DRIFT] Tích hợp đa ngành: thủy văn, thủy lực, địa mạo, sinh thái, kinh tế - xã hội Toàn diện, bao quát mọi thành phần dòng chảy Cực kỳ phức tạp, khó áp dụng tại các nước đang phát triển thiếu số liệu

Tranh luận học thuật cốt lõi diễn ra giữa quan điểm "bảo tồn dòng chảy nguyên sinh" (dựa thuần túy vào tỷ lệ dòng chảy tự nhiên theo Tennant hay chỉ số 7Q10 của Cục Bảo vệ Môi trường Hoa Kỳ - US EPA) và quan điểm "quản lý thích ứng dòng chảy phục vụ phát triển bền vững" (Poff et al., 2010 với khung ELOHA). Tại Việt Nam, các nghiên cứu trước đây của Viện Quy hoạch Thủy lợi hay Viện Khí tượng Thủy văn chủ yếu áp dụng phương pháp thủy văn thuần túy ($Q_{90%}$ hoặc dòng chảy tháng kiệt nhất với tần suất $P = 90% \div 95%$), dẫn đến việc đánh giá thấp nhu cầu thực tế của hạ du hoặc không phản ánh được tác động của các công trình chuyển dòng liên lưu vực.

Nghiên cứu của Nguyễn Thị Kim Dung định vị chính xác khoảng trống học thuật bằng cách so sánh đối chiếu với các tiêu chuẩn quốc tế:

  • So sánh với Khung Tiêu chuẩn Châu Âu (WFD 2000/60/EC) & Scotland (SEPA): Cơ quan Bảo vệ Môi trường Scotland quy định phân hạng hệ sinh thái dựa trên tỷ lệ phần trăm khai thác từ $Q_{95}$ ($15% \div 35%$). Tuy nhiên, chuẩn này chỉ phù hợp với các sông có chế độ dòng chảy điều hòa, không phản ánh được tính chất lũ tiểu mãn và kiệt sâu đặc thù của miền Trung Việt Nam.
  • So sánh với Khung tiếp cận BBM (Nam Phi): Phương pháp Building Block Methodology (BBM) xây dựng các khối dòng chảy đáy và dòng chảy lũ sinh thái. Tác giả đã kế thừa tư duy phân khối chức năng của BBM nhưng cải tiến căn bản: cấu trúc lại thành 3 thành phần đo đếm và kiểm soát được trực tiếp tại hiện trường thông qua hệ thống mô hình toán thủy lực - chất lượng nước phân tán.

Sự tiến bộ học thuật của luận án nằm ở việc chuyển đổi từ cách tiếp cận đơn mục tiêu sang khung tổ hợp đa chỉ tiêu, giải quyết bài toán dòng chảy kiệt trong điều kiện nhân vi (bị chi phối bởi 4 nhà máy thủy điện lớn và 756 công trình lấy nước).

Đó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 trị dòng chảy môi trường thông qua việc thiết lập mô hình khái niệm ba thành phần chức năng của DCTT:

$$Q_{TT}(t, s) = \max \left{ Q_{dts}(t, s),, Q_{st}(t, s),, Q_{ktsd}(t, s) \right}$$

Trong đó, mỗi thành phần đóng vai trò như một biên giới hạn lý thuyết (theoretical boundary condition) đối với trạng thái thủy động lực học của dòng sông:

  1. Dòng chảy duy trì sông ($Q_{dts}$): Đại lượng đảm bảo sự tồn tại liên tục về mặt hình thái và thủy văn của lòng dẫn, ngăn chặn hiện tượng đứt dòng, duy trì khả năng tự làm sạch cơ bản, xác định qua đường cong duy trì lưu lượng (FDCA) ở các mức xác suất $P = 90% \div 95%$.
  2. Dòng chảy duy trì sinh thái ($Q_{st}$): Đại lượng đảm bảo không gian sống tối thiểu cho quần xã thủy sinh vật (động vật đáy, thực vật nổi, cá bản địa), được xác định thông qua điểm uốn của đường quan hệ Thủy lực – Chu vi ướt kết hợp tiêu chuẩn chất lượng nước (DO $\ge 5\text{ mg/l}$, $\text{BOD}_5 \le 6\text{ mg/l}$, $\text{NH}_4^+ \le 0,5\text{ mg/l}$, $\text{NO}_3^- \le 5\text{ mg/l}$).
  3. Dòng chảy khai thác sử dụng ($Q_{ktsd}$): Đại lượng động lực học phức hợp nhằm bảo đảm cao trình mực nước lấy nước tự chảy/bơm tưới, đồng thời duy trì nêm thủy lực để ngăn mặn (nồng độ mặn $S \le 0,8 \div 1,0\text{ g/l}$) tại các vị trí lấy nước sinh hoạt (nhà máy nước Cầu Đỏ) và nông nghiệp (trạm bơm Xuyên Đông, Tứ Câu).
   Chu vi ướt (P_w)
       ^
       │                 / Tuyến quan hệ P_w ~ H (hoặc Q)
       │                /
       │               /
       │   Điểm uốn   *  <-- [ĐIỂM TỚI HẠN SINH THÁI]
       │  (Breakpoint)/       Độ dốc cực đại dP_w/dQ
       │             /        Đảm bảo diện tích đáy ướt tối thiểu
       │            /
       │           /
       │          /
       └─────────┴──────────────────────────────> Mực nước (H) / Lưu lượng (Q)

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

Khung phân tích của luận án tích hợp liên hoàn ba lý thuyết nền tảng: Lý thuyết Động lực học dòng sông liên tục (Vannote et al., 1980 - River Continuum Concept), Lý thuyết Cân bằng Thủy văn - Sinh thái (Ecohydrology), và Lý thuyết Phân phối Tài nguyên nước Tối ưu.

Quy trình giải bài toán tổ hợp được thực hiện qua các bước tường minh:

  • Bước 1: Thiết lập mạng lưới các Điểm kiểm soát (ĐKS) đại diện tại các vị trí nhạy cảm: ngã ba nhập lưu/phân lưu (Bung - Cái, Vu Gia - Quảng Huế), sau các thủy điện bậc thang, và thượng lưu các trạm cấp nước lớn.
  • Bước 2: Định lượng độc lập $Q_{dts}$ theo chuỗi số liệu thủy văn lịch sử tự nhiên 32 năm (1977–2008).
  • Bước 3: Sử dụng mô hình thủy lực 1D tính toán quan hệ hình thái học $P_w - H - Q$ để tìm điểm uốn tới hạn xác định $Q_{st}$, sau đó hiệu chỉnh bằng mô hình chất lượng nước MIKE 11 Ecolab.
  • Bước 4: Áp dụng mô hình MIKE BASIN cân bằng nước tổng thể lưu vực kết hợp MIKE 11 AD mô phỏng lan truyền mặn để xác định $Q_{ktsd}$ theo 2 kịch bản: cấp nước thông thường và cấp nước gia tăng.
  • Bước 5: Tổng hợp ma trận quyết định $Q_{TT}$ cho từng ĐKS theo từng thời đoạn 10 ngày trong mùa cạn.

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

Thiết kế nghiên cứu

Nghiên cứu vận hành trên hệ giác quan thực chứng hiện thực phê phán (critical realism), kết hợp chặt chẽ giữa quan trắc thực địa quy mô lớn, phân tích thống kê chuỗi thời gian dài và mô phỏng toán học đa chiều (Multi-level numerical modeling). Thiết kế nghiên cứu phân cấp hệ thống gồm 3 tầng bậc: tầng lưu vực vĩ mô (MIKE BASIN), tầng mạng lưới sông trung mô (MIKE 11 HD/AD), và tầng mặt cắt vi mô (Phương pháp chu vi ướt và Ecolab).

┌─────────────────────────────────────────────────────────────────────────────┐
│                           THIẾT KẾ NGHIÊN CỨU ĐA TẦNG                        │
├─────────────────────────────────────────────────────────────────────────────┤
│  TẦNG LƯU VỰC VĨ MÔ        Mô hình Cân bằng nước MIKE BASIN                 │
│  (Macro-Basin Level)       • Phân bổ 756 hồ chứa, 4 thủy điện lớn            │
│                            • Nhu cầu 40.000 ha tưới & cấp nước đô thị       │
├─────────────────────────────────────────────────────────────────────────────┤
│  TẦNG MẠNG SÔNG TRUNG MÔ   Mô hình Thủy lực & Lan truyền mặn MIKE 11 HD/AD   │
│  (Meso-Network Level)      • Mạng lưới sông phức hợp Vu Gia - Thu Bồn       │
│                            • Động học thủy triều & ranh mặn 0,8 - 1,0 g/l   │
├─────────────────────────────────────────────────────────────────────────────┤
│  TẦNG MẶT CẮT VI MÔ        Phương pháp Chu vi ướt & MIKE 11 Ecolab          │
│  (Micro-Section Level)     • Đường quan hệ P_w ~ H ~ Q tìm điểm uốn         │
│                            • Nồng độ DO, BOD5, NH4+, NO3- thủy sinh         │
└─────────────────────────────────────────────────────────────────────────────┘

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu tuân thủ nghiêm ngặt các tiêu chuẩn kiểm soát chất lượng khoa học:

  • Thu thập và xử lý số liệu: Khai thác chuỗi dữ liệu khí tượng thủy văn đồng bộ 1977–2014 tại 11 trạm khí tượng và 7 trạm thủy văn chính (Thành Mỹ, Nông Sơn, Ái Nghĩa, Giao Thủy, Câu Lâu, Hội Khách, Cẩm Lệ).
  • Mẫu sinh thái thủy sinh: Khảo sát toàn diện về thành phần loài sinh vật nổi, động vật đáy và 42 loài cá có giá trị kinh tế/bảo tồn trên lưu vực sông Vu Gia – Thu Bồn (như cá bạt, cá chình, cá mè hoa).
  • Hiệu chuẩn và kiểm định mô hình thủy lực:
    • Giai đoạn hiệu chỉnh: Sử dụng chuỗi số liệu thực đo năm 2008 và mùa cạn 2013.
    • Giai đoạn kiểm định: Chuỗi số liệu độc lập mùa cạn năm 2009 và 2014.
    • Chỉ số độ tin cậy: Hệ số tương quan Nash-Sutcliffe ($NSE$) tại các trạm thủy văn chính đều đạt mức $NSE > 0,85 \div 0,92$, sai số lưu lượng đỉnh kiệt và mực nước cực tiểu $< 5%$, chứng minh độ hội tụ và độ chính xác xuất sắc của mô hình toán.

Data và phân tích

Phân tích số liệu áp dụng tổ hợp phần mềm chuyên dụng tiêu chuẩn quốc tế của Viện Thủy động học Đan Mạch (DHI):

  1. MIKE BASIN: Phân chia lưu vực thành các tiểu lưu vực, số hóa 756 hồ đập và mạng lưới sông suối bậc 1 đến bậc 5, tích hợp chuỗi mưa phục vụ cân bằng cung - cầu theo từng 10 ngày.
  2. MIKE 11 HD (Hydrodynamic Module): Giải hệ phương trình Saint-Venant đầy đủ cho dòng chảy không ổn định trên toàn bộ mạng sông phân lưu - nhập lưu phức tạp của Vu Gia – Thu Bồn.
  3. MIKE 11 AD (Advection-Dispersion Module): Giải phương trình truyền dẫn - khuếch tán 1D để mô phỏng ranh giới mặn $1,0\text{ g/l}$ và $0,8\text{ g/l}$ lấn sâu vào hạ lưu sông Vu Gia, Thu Bồn và Vĩnh Điện.
  4. MIKE 11 Ecolab: Mô phỏng 6 trạng thái sinh hóa môi trường nước, kiểm tra trực tiếp khả năng tự làm sạch của sông ứng với các cấp lưu lượng xả kiệt.

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

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

Kết quả mô phỏng và tính toán trong luận án mang lại 5 phát hiện cốt lõi có giá trị khoa học và thực tiễn sâu sắc:

Điểm kiểm soát (ĐKS) Vị trí địa lý & Đặc trưng $Q_{dts}$ ($m^3/s$) $Q_{st}$ ($m^3/s$) $Q_{ktsd}$ ($m^3/s$) [Cấp gia tăng] $Q_{TT}$ Đề xuất ($m^3/s$) Yếu tố chi phối chính
Thành Mỹ Thượng lưu Vu Gia $29,60$ $32,50$ $15,20$ $32,50$ Chu vi ướt sinh thái ($Q_{st}$)
Nông Sơn Thượng lưu Thu Bồn $48,20$ $46,50$ $42,10$ $48,20$ Thủy văn duy trì sông ($Q_{dts}$)
Ái Nghĩa Hạ lưu Vu Gia $42,10$ $38,60$ $65,40$ $65,40$ Đẩy mặn NMN Cầu Đỏ ($Q_{ktsd}$)
Giao Thủy Hạ lưu Thu Bồn $65,80$ $58,20$ $87,50$ $87,50$ Trạm bơm tưới Xuyên Đông ($Q_{ktsd}$)
  1. Sự chi phối của Dòng chảy khai thác sử dụng tại hạ du: Tại các ĐKS hạ lưu (Ái Nghĩa và Giao Thủy), $Q_{ktsd}$ vượt trội hoàn toàn so với $Q_{dts}$ và $Q_{st}$. Cụ thể tại Ái Nghĩa, lưu lượng tối thiểu đòi hỏi phải đạt $65,40\text{ m}^3\text{/s}$ trong thời kỳ cấp nước gia tăng (vụ Hè Thu) nhằm đẩy mặn dưới $0,8\text{ g/l}$ tại cửa thu nhà máy nước Cầu Đỏ, cao hơn $55%$ so với nhu cầu dòng chảy duy trì sông thuần túy ($42,10\text{ m}^3\text{/s}$).
  2. Ngưỡng sinh thái tại thượng nguồn: Tại ĐKS Thành Mỹ (hạ lưu ngã ba Bung - Cái), phương pháp chu vi ướt xác định điểm uốn tương ứng với $Q_{st} = 32,50\text{ m}^3\text{/s}$, cao hơn giá trị dòng chảy kiệt tự nhiên trung bình nhiều năm ($29,60\text{ m}^3\text{/s}$). Điều này chứng minh rằng việc khai thác thủy điện thượng lưu nếu chỉ trả về lưu lượng thủy văn tối thiểu sẽ phá hủy sinh cảnh đáy ướt của các loài thủy sinh.
  3. Hiện tượng nghịch lý phân lưu tại ngã ba Quảng Huế: Tác giả phát hiện luồng phân lưu tự nhiên qua sông Quảng Huế bị biến động mạnh sau khi thủy điện Đắk Mi 4 chuyển dòng sang Thu Bồn, làm suy giảm đến $72%$ lưu lượng tự nhiên chảy về Vu Gia trong mùa cạn nếu không có cơ chế bù dòng cưỡng bức.
  4. Động học xâm nhập mặn phi tuyến: Quan hệ giữa lưu lượng xả thượng lưu và chiều dài xâm nhập mặn thể hiện dạng hàm mũ suy giảm; khi lưu lượng tại Ái Nghĩa hạ xuống dưới $45\text{ m}^3\text{/s}$, ranh mặn $1\text{ g/l}$ xâm nhập sâu vượt qua trạm bơm Cầu Đỏ hơn $8,5\text{ km}$, làm tê liệt hoàn toàn hệ thống cấp nước sinh hoạt thành phố Đà Nẵng.

Implications đa chiều

  • Về mặt lý thuyết: Luận án chứng minh thuyết phục rằng việc áp dụng một tiêu chuẩn tỷ lệ % dòng chảy cố định (như phương pháp Tennant) là hoàn toàn sai lầm đối với các lưu vực có hệ số biến động dòng chảy mùa kiệt sâu và có sự can thiệp công trình liên lưu vực.
  • Về mặt phương pháp luận: Thiết lập quy trình chuẩn hóa kết nối mô hình thủy lực 1D và phương pháp chu vi ướt, cho phép chuyển đổi từ dữ liệu mặt cắt hình thái sang giá trị lưu lượng sinh thái có thể kiểm chứng tại hiện trường.
  • Về mặt quản lý thực tiễn: Kết quả tính toán cung cấp cơ sở kỹ thuật trực tiếp để điều chỉnh Quy trình vận hành liên hồ chứa 1537/QĐ-TTg, chuyển đổi phương thức giám sát từ mực nước thụ động sang kiểm soát lưu lượng xả chủ động theo thời đoạn 10 ngày.
  • Về chính sách: Cung cấp luận cứ khoa học để Bộ Tài nguyên và Môi trường, Bộ Nông nghiệp và PTNT cùng UBND TP. Đà Nẵng và tỉnh Quảng Nam thiết lập cơ chế điều phối liên tỉnh trong chia sẻ nguồn nước mùa cạn.

Limitations và Future Research

Nghiên cứu thể hiện tính trung thực khoa học thông qua việc làm rõ các giới hạn:

  1. Dữ liệu sinh thái chi tiết: Chưa tích hợp đầy đủ mô hình đường cong tương thích sinh thái (Habitat Suitability Curves - HSC) cho từng giai đoạn sinh trưởng riêng biệt (trứng, ấu trùng, cá thể trưởng thành) của từng loài cá đặc hữu do thiếu chuỗi quan trắc sinh thái dài hạn.
  2. Động lực học biến hình lòng dẫn: Mô hình thủy lực MIKE 11 giả định mặt cắt hình thái cố định trong chu kỳ tính toán, chưa phản ánh được hiện tượng xói lở đáy sông và biến động bãi bồi do nạn khai thác cát trái phép tại hạ lưu Vu Gia - Thu Bồn.
  3. Tác động tích hợp của BĐKH và Nước biển dâng: Luận án mới tính toán cho các kịch bản khí hậu hiện tại và kịch bản tĩnh, chưa xem xét kịch bản liên hoàn triều cường kết hợp nước dâng do bão trong giai đoạn chuyển mùa.

Chương trình nghiên cứu tương lai (Future Research Agenda) định hướng:

  • Mở rộng ứng dụng mô hình 2D/3D (MIKE 21/3 hoặc Delft3D) tại các vùng ngã ba phân lưu phức tạp và vùng cửa sông Cửa Đại, Cửa Hàn.
  • Ứng dụng kỹ thuật viễn thám và dữ liệu ảnh vệ tinh radar để theo dõi độ ẩm đất ướt và diễn biến bãi cạn sinh thái theo thời gian thực.
  • Tích hợp mô hình kinh tế thủy văn (Hydro-economic modeling) nhằm tối ưu hóa bài toán phân bổ nước giữa giá trị kinh tế phát điện và giá trị sinh thái - nông nghiệp hạ du.

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

  • Ảnh hưởng học thuật: Luận án thiết lập một chuẩn mực mới trong nghiên cứu DCTT tại Việt Nam, đóng góp phương pháp luận có thể chuyển giao cho 9 lưu vực sông lớn còn lại của quốc gia (như sông Hồng, sông Ba, sông Đồng Nai).
  • Tác động ngành thủy điện và nông nghiệp: Cung cấp khung định lượng chính xác để tính toán tổn thất sản lượng điện năng khi các nhà máy thủy điện (A Vương, Đắk Mi 4, Sông Bung 4) phải xả dòng chảy môi trường, giúp các cơ quan quản lý năng lượng và tưới tiêu xây dựng phương án kinh tế tối ưu.
  • Tác động chính sách quốc gia: Trở thành tài liệu tham khảo kỹ thuật then chốt phục vụ việc rà soát, sửa đổi Luật Tài nguyên nước và các Thông tư hướng dẫn quy định DCTT trên sông suối hạ lưu hồ chứa.
  • Lợi ích xã hội và môi trường: Bảo vệ an ninh nước sinh hoạt cho hơn 1,2 triệu dân thành phố Đà Nẵng và bảo đảm nước tưới cho gần 40.000 ha đất canh tác lúa, cây hoa màu tại tỉnh Quảng Nam trong các đợt hạn hán lịch sử.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận phương pháp luận liên ngành kết hợp giữa Thủy văn thống kê, Thủy lực hình thái và Mô hình hóa môi trường nước.
  • Cơ quan quản lý Nhà nước (BTNMT, BNNPTNT, Cục Quản lý Tài nguyên nước): Sở hữu công cụ khoa học để thẩm định quy trình vận hành hồ chứa và cấp phép khai thác tài nguyên nước.
  • Các Ban Quản lý Lưu vực sông & UBND Tỉnh/Thành phố (Quảng Nam, Đà Nẵng): Có căn cứ pháp lý - kỹ thuật minh bạch để giải quyết xung đột lợi ích sử dụng nước liên tỉnh.
  • Doanh nghiệp Thủy điện & Thủy lợi: Chủ động lập kế hoạch phát điện và điều tiết hồ chứa tuân thủ quy chuẩn bảo vệ môi trường, giảm thiểu rủi ro pháp lý và tranh chấp cộng đồng.

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

1. Đóng góp lý thuyết độc đáo nhất của công trình là gì và mở rộng lý thuyết nào? Công trình đã mở rộng Lý thuyết Động lực học Dòng chảy Môi trường (Environmental Flow Dynamics) trong điều kiện nhân vi đặc thù (lưu vực bị chuyển dòng liên lưu vực). Đóng góp độc đáo nhất là xác lập công thức tổ hợp động lực $Q_{TT} = \max(Q_{dts}, Q_{st}, Q_{ktsd})$ theo từng thời đoạn 10 ngày, chứng minh rằng tại lưu vực bị chi phối mạnh mẽ bởi con người, thành phần $Q_{ktsd}$ (chi phối bởi nêm mặn và cao trình bơm) có thể thay thế thành phần thủy văn tự nhiên để trở thành điều kiện biên khống chế toàn hệ thống.

2. Đột phá về phương pháp luận so với các nghiên cứu tiền nhiệm tại Việt Nam và quốc tế? So với các nghiên cứu của Tennant (1976) thuần túy thủy văn hay Richter (1997) thuần túy thống kê, luận án đã đột phá bằng cách tích hợp trực tiếp phương pháp Thủy lực Chu vi ướt (Wetted Perimeter) với mô hình số trị chất lượng nước MIKE 11 Ecolab. Thay vì chọn điểm uốn hình học đơn thuần, nghiên cứu kiểm chứng điểm uốn đó dưới áp lực của các thông số hóa sinh ($DO, BOD_5, NH_4^+, NO_3^-$), tạo ra quy trình xác định dòng chảy sinh thái có độ tin cậy vượt trội.

3. Phát hiện bất ngờ nhất có số liệu thực nghiệm chứng minh là gì? Phát hiện bất ngờ nhất là tại ĐKS Thành Mỹ: mặc dù nằm ở khu vực miền núi thượng lưu Vu Gia, nhu cầu dòng chảy sinh thái ($Q_{st} = 32,50\text{ m}^3\text{/s}$) lại cao hơn dòng chảy kiệt tự nhiên trung bình nhiều năm ($29,60\text{ m}^3\text{/s}$). Điều này chứng minh rằng chế độ dòng chảy kiệt lịch sử từng khiến các loài thủy sinh đáy bị suy thoái cục bộ, và việc duy trì một dòng chảy môi trường phục hồi đòi hỏi lượng xả lớn hơn mức tối thiểu tự nhiên.

4. Luận án có cung cấp quy trình nhân bản (Replication Protocol) không? Có. Luận án cung cấp một quy trình 5 bước tường minh, kèm theo hệ thống tham số thủy lực (độ nhám Manning $n$, hệ số phân tán $D$), sơ đồ mạng lưới MIKE BASIN và MIKE 11 chi tiết, cho phép các kỹ sư tài nguyên nước nhân bản hoàn toàn phương pháp này cho bất kỳ lưu vực sông nào tại Việt Nam có đặc điểm tương đồng.

5. Định hướng chương trình nghiên cứu 10 năm tiếp theo từ kết quả luận án? Xây dựng hệ thống Trợ giúp Ra quyết định Thời gian thực (Real-time Decision Support System - DSS) tích hợp Trí tuệ Nhân tạo (AI/Machine Learning) với mô hình thủy lực 1D-2D liên hợp, kết nối mạng lưới trạm đo mưa tự động và cảm biến độ mặn IoT để tự động hóa hoàn toàn quy trình điều hành xả DCTT của liên hồ chứa Vu Gia – Thu Bồn.

Kết luận

Công trình nghiên cứu của Tiến sĩ Nguyễn Thị Kim Dung đã giải quyết trọn vẹn một bài toán khoa học và thực tiễn cấp bách, tạo lập các giá trị cốt lõi:

  1. Xây dựng thành công cơ sở khoa học và phương pháp luận toàn diện để xác định DCTT trên lưu vực sông thông qua tổ hợp 3 thành phần: duy trì sông ($Q_{dts}$), duy trì sinh thái ($Q_{st}$) và khai thác sử dụng ($Q_{ktsd}$).
  2. Định lượng chính xác bộ giá trị DCTT tại các điểm kiểm soát yết hầu của hệ thống sông Vu Gia – Thu Bồn (Thành Mỹ: $32,50\text{ m}^3\text{/s}$; Nông Sơn: $48,20\text{ m}^3\text{/s}$; Ái Nghĩa: $65,40\text{ m}^3\text{/s}$; Giao Thủy: $87,50\text{ m}^3\text{/s}$ trong thời kỳ cấp nước gia tăng).
  3. Làm sáng tỏ cơ chế thủy lực và động học xâm nhập mặn trong điều kiện chuyển nước nhân tạo của thủy điện Đắk Mi 4, cung cấp cơ sở kỹ thuật vững chắc để bảo vệ nguồn nước sinh hoạt cho thành phố Đà Nẵng.
  4. Chuẩn hóa quy trình tích hợp công cụ mô hình toán (MIKE BASIN, MIKE 11 HD/AD/Ecolab và phương pháp Chu vi ướt), đạt độ tin cậy khoa học cao với hệ số Nash-Sutcliffe $NSE > 0,85$.
  5. Mở ra hướng tiếp cận mới trong quản trị lưu vực sông tích hợp (IWRM), tạo bước chuyển dịch quan trọng từ tư duy quản lý tĩnh sang quản lý động thái thích ứng với biến đổi khí hậu và phát triển thủy điện bậc thang.