Luận án TS Thủy văn: Phân bổ hợp lý nguồn nước mặt lưu vực sông Vu Gia - Thu Bồn

Tổng quan về luận án

Nghiên cứu của tác giả Thân Văn Đón (2021) với tiêu đề "Nghiên cứu phân bổ hợp lý nguồn nước mặt cho lưu vực sông Vu Gia - Thu Bồn" (Luận án Tiến sĩ Thủy văn học, Mã số: 9440 224, Viện Khoa học Khí tượng Thủy văn và Biến đổi khí hậu) giải quyết bài toán cốt lõi trong quản trị tài nguyên nước: thiết lập cơ sở khoa học và pháp lý định lượng cho việc chia sẻ, điều tiết nguồn nước mặt giữa các thực thể hành chính (tỉnh Quảng Nam và thành phố Đà Nẵng) và các ngành kinh tế cạnh tranh trực tiếp. Lưu vực sông Vu Gia - Thu Bồn có diện tích tự nhiên 10.350 km² (trong đó 560,5 km² thuộc tỉnh Kon Tum, phần lớn còn lại thuộc Quảng Nam và Đà Nẵng, tọa độ 16°03’ - 14°55’ vĩ độ Bắc và 107°15’ - 108°24’ kinh độ Đông). Đây là hệ thống thủy văn phức hợp sườn Đông dãy Trường Sơn, nơi 95% lượng nước từ vùng thượng nguồn đổ về hạ lưu giữ vai trò quyết định tới an ninh nguồn nước của 1.518.000 dân cư (mật độ trung bình 207 người/km², cá biệt 6 quận nội thành Đà Nẵng lên đến 3.458 người/km²) cùng hệ thống 11 khu công nghiệp tập trung và 59.789 ha đất canh tác nông nghiệp.

                  ┌─────────────────────────────────────────────────────────┐
                  │ TỔNG LƯỢNG NƯỚC MẶT LƯU VỰC: ~24 TỶ M³/NĂM (100%)       │
                  └───────────────────────────┬─────────────────────────────┘
                                              │
                    ┌─────────────────────────┴─────────────────────────┐
                    ▼                                                   ▼
┌───────────────────────────────────────┐   ┌───────────────────────────────────────┐
│     MÙA LŨ (Tháng IX/X - XII)         │   │       MÙA KHÔ (Tháng I - VIII)        │
│  - Chiếm 63% - 68% tổng lượng nước    │   │  - Chiếm 32% - 37% tổng lượng nước    │
│  - Lưu lượng đỉnh:                    │   │  - Tháng IV kiệt nhất: chỉ 2% - 3%    │
│    + Nông Sơn: 5.600 m³/s             │   │  - Lưu lượng kiệt kỷ lục:             │
│    + Thành Mỹ: 3.778 m³/s             │   │    + Thành Mỹ: 11,5 m³/s (1988)       │
└───────────────────────────────────────┘   └───────────────────┬───────────────────┘
                                                                │
                                            ┌───────────────────┴───────────────────┐
                                            ▼                                       ▼
┌───────────────────────────────────────────────────────┐ ┌────────────────────────────────────────┐
│ CĂNG THẲNG THỦY ĐIỆN VÀ CHUYỂN DÒNG LIÊN LƯU VỰC       │ │ HỆ LỤY VÀ XUNG ĐỘT HẠ LƯU              │
│ - Mật độ dày đặc: Đăk Mi 4, A Vương, Sông Bung 4,...  │ │ - Nhà máy nước Cầu Đỏ bị nhiễm mặn     │
│ - Đăk Mi 4 chuyển nước từ Vu Gia sang Thu Bồn         │ │ - Tranh chấp gay gắt: Quảng Nam - Đà Nẵng│
│ - Làm suy kiệt hạ lưu sông Vu Gia trong mùa kiệt      │ │ - Thiếu nước tưới cho 40.332 ha lúa    │
└───────────────────────────────────────────────────────┘ └────────────────────────────────────────┘

Khoảng trống nghiên cứu (Research Gap) xuất phát từ thực trạng các quy hoạch tài nguyên nước trước đây (Lê Thế Trung, 2016; Viện Quy hoạch Thủy lợi, 2010) chủ yếu tiếp cận theo hướng "quản lý cung" (supply-side management) đơn ngành, chưa tích hợp được mối liên kết động giữa cơ chế pháp lý (Luật Tài nguyên nước 2012, Thông tư 04/2020/TT-BTNMT) với công cụ mô hình toán phân bổ tối ưu trong điều kiện hạn hán, thiếu nước cực đoan. Đặc biệt, việc công trình thủy điện Đăk Mi 4 chuyển dòng liên lưu vực từ sông Vu Gia sang sông Thu Bồn đã làm biến đổi cấu trúc thủy lực tự nhiên, làm trầm trọng thêm tình trạng suy giảm dòng chảy kiệt và xâm nhập mặn sâu (ranh mặn 1‰ lấn sâu 15-25 km trên sông Vĩnh Điện, đẩy mặn tại trạm bơm Cầu Đỏ - nguồn cấp nước sinh hoạt chính của Đà Nẵng).

Nghiên cứu đặt ra hai câu hỏi học thuật và hệ thống giả thuyết tương ứng:

• Câu hỏi nghiên cứu 1 (RQ1): Trong điều kiện biến động thủy văn hiện tại và tương lai (năm 2030) dưới tác động của BĐKH, lượng nước mặt có thể phân bổ ($W_{pb}$) của từng tiểu vùng trên lưu vực sông Vu Gia - Thu Bồn là bao nhiêu khi đã khấu trừ lượng nước bắt buộc cho dòng chảy tối thiểu ($W_{tt}$) và nhu cầu thiết yếu ($W_{ty}$)?
  • Giả thuyết 1 (H1): Vào mùa kiệt ở các năm có tần suất dòng chảy $P=85%$ và $P=95%$, nguồn nước tự nhiên không đáp ứng đủ tổng nhu cầu dùng nước của các ngành nếu không có sự tái phân bổ và điều tiết dung tích hữu ích từ hệ thống hồ chứa thượng nguồn.
• Câu hỏi nghiên cứu 2 (RQ2): Thứ tự ưu tiên và tỷ lệ phân bổ lượng nước thiếu hụt giữa các ngành dùng nước (sinh hoạt, môi trường sinh thái, nông nghiệp, công nghiệp, thủy điện) cùng cơ chế chia sẻ lợi ích - trách nhiệm liên tỉnh cần được thiết lập như thế nào để tối đa hóa hiệu quả kinh tế - xã hội và bảo đảm tính bền vững môi trường?
  • Giả thuyết 2 (H2): Việc áp dụng khung ma trận phân bổ đa tiêu chí tích hợp mô hình Mike NAM - WEAP sẽ giảm thiểu tổn thất kinh tế do thiếu hụt nước tới mức tối đa, đồng thời duy trì dòng chảy sinh thái đẩy mặn tốt hơn so với cơ chế phân bổ hành chính tĩnh.

Khung lý thuyết của luận án xây dựng trên nền tảng Lý thuyết Quản lý Tổng hợp Tài nguyên Nước (Integrated Water Resources Management - IWRM), Lý thuyết Tối ưu hóa Kinh tế - Thủy văn (Hydro-Economic Optimization), và Khung lý thuyết Phân bổ Nguồn nước Bền vững (Sustainable Water Allocation Framework) theo định nghĩa chuẩn hóa của Viện Quản lý Nước Quốc tế (IWMI): "Phân bổ nguồn nước (water allocation) là sự chia sẻ nguồn nước giữa các ngành dùng nước với nhau trên lưu vực sông để đảm bảo lợi ích của tất cả người dùng". Đóng góp đột phá của công trình là lượng hóa chính xác lượng nước có thể phân bổ hàng năm thông qua phương trình cấu trúc:

$$W_{pb} = W_{sd} - (W_{tt} + W_{ty})$$

Trong đó: $W_{sd}$ là lượng nước có thể sử dụng; $W_{tt}$ là lượng nước đảm bảo dòng chảy tối thiểu; $W_{ty}$ là lượng nước đáp ứng nhu cầu sinh hoạt thiết yếu.

Phạm vi nghiên cứu bao quát chuỗi số liệu thủy văn 40 năm (1977 - 2016), phân tích hiện trạng 2014 và dự báo nhu cầu đến 2030 cho toàn bộ hệ thống 10.350 km² lưu vực sông Vu Gia - Thu Bồn.

Literature Review và Positioning

Lịch sử nghiên cứu phân bổ nguồn nước trên thế giới đã tiến hóa qua nhiều giai đoạn: từ các quy tắc phân chia nước cục bộ thời cổ đại như Bộ luật Hammurabi (1795-1750 TCN) tại lưu vực Tigris - Euphrates (Dinar et al., 1997; Speed et al., 2013), đến các hiệp ước pháp lý phân chia hạn ngạch cứng giữa các bang hoặc quốc gia. Điển hình như Hiệp ước Sông Colorado năm 1922/2012 giữa Mỹ và Mexico (phân chia cố định 9,25 tỷ m³ cho hạ lưu trong 10 năm và xả tối thiểu 1,85 tỷ m³/năm qua biên giới) bộc lộ nhược điểm nghiêm trọng khi không thể thích ứng với hạn hán kéo dài và hoàn toàn bỏ quên dòng chảy môi trường (Garrick et al., 2012; WWF, 2013). Tương tự, Hiệp định Phân bổ Nước sông Indus năm 1991 (Water Apportionment Accord) giữa các tỉnh của Pakistan (tổng lượng phân bổ 144 tỷ m³, trong đó tỉnh Punjab nhận 55,9 tỷ m³, Sindh nhận 48,8 tỷ m³) cũng vấp phải tranh cãi gay gắt do thiếu cơ chế chia sẻ dòng chảy sinh thái hạ lưu và thiếu tính dự báo động (Briscoe & Qamar, 2006; IWMI, 2007).

┌───────────────────────────────────────────────────────────────────────────────────────────────────┐
│                                 TIẾN TRÌNH TIẾP CẬN PHÂN BỔ NGUỒN NƯỚC                             │
└───────────────────────────────────────────────────────────────────────────────────────────────────┘
                                                  │
 ┌────────────────────────────────────────────────┴────────────────────────────────────────────────┐
 ▼                                                                                                 ▼
TRƯỜNG PHÁI PHÂN CHIA HẠN NGẠCH CỐ ĐỊNH                              TRƯỜNG PHÁI QUẢN TRỊ THÍCH ỨNG & ĐA TIÊU CHÍ
- Quản lý từ trên xuống (Top-down)                                   - Tích hợp thủy văn - kinh tế - sinh thái (IWRM)
- Phân bổ tĩnh theo địa giới hành chính                              - Tối ưu hóa phân bổ theo $W_{pb} = W_{sd} - (W_{tt} + W_{ty})$
- Thiếu điều chỉnh theo năm khô hạn                                  - Ràng buộc dòng chảy tối thiểu & đẩy mặn
- Tiêu biểu: Colorado Compact (1922), Indus Accord (1991)            - Tiêu biểu: Murray-Darling (2008), Incomati (2002)
                                                  │
                                                  ▼
                         ┌─────────────────────────────────────────────────┐
                         │ ĐỊNH VỊ CỦA CÔNG TRÌNH NGHIÊN CỨU               │
                         │ Áp dụng mô hình động Mike NAM - WEAP            │
                         │ Tích hợp đa tiêu chí ưu tiên và cơ chế chia sẻ  │
                         │ lợi ích liên tỉnh cho Vu Gia - Thu Bồn          │
                         └─────────────────────────────────────────────────┘

Trong bức tranh học thuật quốc tế, có hai trường phái đối lập chính trong phân bổ nguồn nước:

1. Trường phái Quản lý phân bổ tập trung theo hạn ngạch pháp lý/hành chính (Administrative/Command-and-Control Approach): Đại diện bởi các nghiên cứu ban đầu về sông Colorado và sông Nile (Allan, 2002; Wolf, 1998), ưu tiên tính ổn định pháp lý nhưng kém linh hoạt khi xảy ra biến động nguồn nước đến.
2. Trường phái Kinh tế - Sinh thái thị trường và Quản trị thích ứng (Hydro-economic & Adaptive Allocation School): Đại diện bởi trường phái CSIRO (2008) tại lưu vực sông Murray-Darling (Úc) với mô hình phân bổ 18 tiểu vùng kết hợp công cụ Mike Basin - GMS, cắt giảm 20% lượng nước dùng để phục hồi hệ sinh thái và chính phủ chi 3,1 tỷ AUD mua lại quyền khai thác nước; hay mô hình phân chia theo mức độ ưu tiên trên lưu vực sông Incomati (Nam Phi - Mozambique - Swaziland) giải quyết công bằng xã hội và hạn ngạch sinh thái (van der Zaag, 2005; Le Quesne et al., 2010); và mô hình động thái LERMA 6 module cho 17 tiểu lưu vực sông Lerma-Chapala tại Mexico (de Silva et al., 2004).

Tại Việt Nam, các nghiên cứu trước đây đã từng bước tiếp cận bài toán cân bằng nước. Nghiên cứu của Tô Trung Nghĩa và Lê Mạnh Hùng (2008), Ringler & Huy (2004) áp dụng mô hình tối ưu hóa GAMS cho lưu vực sông Hồng và sông Đồng Nai. Dương Văn Khánh (2011) xây dựng bộ tiêu chí quy hoạch tài nguyên nước cho lưu vực sông Ba; Nguyễn Chí Công (2014) xác định tỷ lệ chia sẻ nước mùa kiệt cho sông Đồng Nai ở các tần suất $P=75%$ và $P=85%$; Nguyễn Ngọc Hà (2016) kết hợp WEAP, WRAM và LINGO để tối ưu hóa giá trị kinh tế nước trên sông Vệ. Riêng tại lưu vực sông Vu Gia - Thu Bồn, các công trình của Viện Quy hoạch Thủy lợi (2002, 2010), Viện Khoa học Thủy lợi Việt Nam (2015), Trần Thị Diễn (2020) đã khảo sát mâu thuẫn vận hành hồ chứa thủy điện lớn (A Vương, Sông Bung 4, Đăk Mi 4, Sông Tranh 2) và tính toán dòng chảy tối thiểu.

Tuy nhiên, phần lớn các công trình trong nước chỉ dừng lại ở bài toán "cân bằng cung - cầu thụ động" hoặc tính toán thủy lực đơn lẻ mà chưa thiết lập được:

• Phương pháp luận phân bổ theo lượng nước có thể phân bổ thực tế ($W_{pb}$).
• Phân bổ theo chuỗi thời gian nhiều kịch bản tần suất dòng chảy kiệt ($P=85%$, $P=95%$).
• Khung chia sẻ lợi ích, trách nhiệm mang tính ràng buộc thể chế đa ngành, liên tỉnh. Luận án của Thân Văn Đón định vị chuẩn xác tại điểm giao thoa này, đưa quản trị lưu vực Vu Gia - Thu Bồn tiệm cận trình độ quản lý lưu vực tiên tiến của thế giới.

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

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

Công trình mở rộng và cụ thể hóa ba cấu trúc lý thuyết nền tảng trong khoa học thủy văn và quản lý tài nguyên nước:

┌──────────────────────────────────────────────────────────────────────────────────────────────────┐
│                           HỆ THỐNG ĐÓNG GÓP LÝ THUYẾT CỦA CÔNG TRÌNH                             │
└──────────────────────────────────────────────────────────────────────────────────────────────────┘
                                                 │
      ┌──────────────────────────────────────────┼──────────────────────────────────────────┐
      ▼                                          ▼                                          ▼
┌───────────────────────────┐      ┌───────────────────────────┐      ┌───────────────────────────┐
│ LÝ THUYẾT QUẢN LÝ TỔNG HỢP │      │ LÝ THUYẾT DÒNG CHẢY SINH  │      │ LÝ THUYẾT THỂ CHẾ VÀ      │
│ TÀI NGUYÊN NƯỚC (IWRM)    │      │ THÁI & SỨC CHỊU TẢI       │      │ CHIA SẺ LỢI ÍCH           │
│ (Biswas, 2004; GWP)       │      │ (Poff et al., 1997)       │      │ (Sadoff & Grey, 2002)     │
│ - Chuyển dịch từ quản trị │      │ - Định lượng $W_{tt}$     │      │ - Giải quyết xung đột     │
│   cung sang quản trị cầu  │      │   ràng buộc cứng          │      │   chuyển dòng Đăk Mi 4    │
│ - Tích hợp thể chế - kỹ   │      │ - Nội hàm hóa ngăn mặn    │      │ - Cơ chế đền bù và phối   │
│   thuật - kinh tế lượng   │      │   vào phương trình $W_{pb}$│      │   hợp vận hành liên tỉnh  │
└───────────────────────────┘      └───────────────────────────┘      └───────────────────────────┘

1. Phát triển Lý thuyết Quản lý Tổng hợp Tài nguyên Nước (IWRM - Global Water Partnership, 2000; Biswas, 2004): Luận án chuyển đổi quan niệm từ "phân phối nguồn nước có sẵn" sang "phân bổ theo dung lượng thích ứng của hệ thống tự nhiên". Lý thuyết hệ thống tài nguyên nước phức hợp được chuẩn hóa bằng định nghĩa rõ ràng: "Hệ thống tài nguyên nước là một hệ thống phức tạp bao gồm nguồn nước, các công trình khai thác tài nguyên nước, các yêu cầu về nước cùng với mối quan hệ tương tác giữa chúng và sự tác động của môi trường lên nó".
2. Làm sâu sắc Lý thuyết Dòng chảy Sinh thái và Sức chịu tải Thủy văn (Environmental Flows Theory - Poff et al., 1997; Richter et al., 2003): Công trình xác lập nguyên tắc dòng chảy tối thiểu ($W_{tt}$) và nhu cầu sinh hoạt cơ bản ($W_{ty}$) không phải là đối tượng tham gia cạnh tranh cắt giảm, mà là "ngưỡng sinh thái bất biến" phải khấu trừ tiên quyết trước khi xác định lượng nước phân bổ cho phát triển sản xuất.
3. Mở rộng Lý thuyết Thể chế và Chia sẻ Lợi ích (Benefit-Sharing Theory - Sadoff & Grey, 2002): Chứng minh rằng việc giải quyết mâu thuẫn thủy điện - hạ du không thể chỉ giải quyết bằng ranh giới hành chính mà phải thông qua cơ chế đền bù kinh tế đa ngành và chia sẻ quyền kiểm soát vận hành hồ chứa theo thời gian thực.

Mô hình lý thuyết đề xuất bao gồm 3 mệnh đề khoa học (Propositions):

• Mệnh đề 1 ($P_1$): Lượng nước mặt có thể phân bổ ($W_{pb}$) của một lưu vực có tính chất chuyển dòng biến đổi phi tuyến tính theo chu kỳ xả của các hồ chứa bậc thang thủy điện.
• Mệnh đề 2 ($P_2$): Khi xảy ra hạn hán với tần suất $P \ge 85%$, thứ bậc ưu tiên phân bổ bắt buộc phải chuyển đổi từ tối đa hóa hiệu quả kinh tế (thủy điện, công nghiệp) sang bảo vệ an ninh nước sinh hoạt và duy trì dòng chảy đẩy mặn hạ du.
• Mệnh đề 3 ($P_3$): Tổn thất phúc lợi xã hội đạt mức cực tiểu khi và chỉ khi tỷ lệ cắt giảm nguồn nước giữa các tiểu vùng được xác lập thông qua ma trận chi phí cơ hội thay vì cắt giảm cơ học đồng loạt.

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

Khung phân tích của luận án kết hợp đồng bộ 3 trụ cột phương pháp luận:

• Tích hợp liên mô hình Mike NAM - WEAP: Sử dụng Mike NAM để mô phỏng dòng chảy mặt từ lượng mưa phân bố theo không gian (phương pháp đa giác Thiessen), sau đó kết nối toàn bộ dữ liệu dòng chảy đến làm đầu vào biên cho mô hình đánh giá và lập kế hoạch nước WEAP.
• Tiếp cận cấu trúc chuỗi khối phân vùng kép: Phân chia lưu vực Vu Gia - Thu Bồn thành các tiểu vùng tính toán dựa trên ranh giới lưu vực tự nhiên và các nút khống chế thủy văn (trạm Thành Mỹ, Ái Nghĩa, Nông Sơn, Giao Thủy).
• Ranh giới áp dụng (Boundary Conditions): Khung phân tích áp dụng nghiêm ngặt cho điều kiện thủy văn mùa cạn (tháng I đến tháng VIII), đặc biệt tập trung vào các tháng khô hạn gay gắt nhất (tháng IV, V, VI và tháng IX chuyển tiếp) với chuỗi số liệu được đồng nhất hóa giai đoạn 1977 - 2016.

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

┌──────────────────────────────────────────────────────────────────────────────────────────────────┐
│                         QUY TRÌNH PHƯƠNG PHÁP NGHIÊN CỨU TOÀN DIỆN                               │
└──────────────────────────────────────────────────────────────────────────────────────────────────┘
                                                 │
 ┌───────────────────────────────────────────────┼───────────────────────────────────────────────┐
 ▼                                               ▼                                               ▼
THU THẬP & CHUẨN HÓA DỮ LIỆU            MÔ PHỎNG DÒNG CHẢY MẶT (MIKE NAM)       ĐÁNH GIÁ & PHÂN BỔ NGUỒN NƯỚC (WEAP)
- Khí tượng - thủy văn: 40 năm          - Phân bố mưa: Đa giác Thiessen         - Mạng lưới sơ đồ nút lưu vực
  (1977-2016) tại 10+ trạm              - Hiệu chỉnh: Trạm Nông Sơn, Thành Mỹ   - Cài đặt thứ bậc ưu tiên (1 -> 5)
- Không gian: GIS bản đồ 1/50.000       - Kiểm định: Chỉ số NSE > 0.85          - Tính toán cân bằng theo các kịch bản
- Khảo sát thực địa: 15-30/7/2017       - Xuất dòng chảy đến cho các tiểu vùng    (P=85%, P=95% cho năm 2014 & 2030)
                                                 │
                                                 ▼
┌──────────────────────────────────────────────────────────────────────────────────────────────────┐
│ LỰA CHỌN PHƯƠNG ÁN TỐI ƯU BẰNG MA TRẬN ĐA TIÊU CHÍ (MCDA) VÀ THIẾT LẬP CƠ CHẾ CHIA SẺ LIÊN TỈNH │
└──────────────────────────────────────────────────────────────────────────────────────────────────┘

Thiết kế nghiên cứu

Luận án tuân thủ thế giới quan thực chứng kết hợp hiện thực phê phán (Positivism & Critical Realism), sử dụng phương pháp hỗn hợp (Mixed Methods) lấy mô hình hóa toán học làm trung tâm, kết hợp điều tra khảo sát xã hội học và thể chế kinh tế.

Thiết kế nghiên cứu đa tầng (Multi-level design) được cấu trúc chặt chẽ:

• Cấp độ 1 (Vĩ mô - Toàn lưu vực): Phân tích tương quan tổng tài nguyên nước mặt trung bình nhiều năm đạt 24 tỷ m³/năm (tương đương 24 km³), chiếm 2,6% tổng lượng dòng chảy toàn quốc và 6,8% lượng dòng chảy nội địa sinh ra trên lãnh thổ Việt Nam.
• Cấp độ 2 (Trung gian - Tiểu vùng thủy văn): Phân chia thành các tiểu vùng tính toán dựa trên các phân lưu Quảng Huế (dẫn nước Vu Gia sang Thu Bồn) và sông Vĩnh Điện (dẫn nước Thu Bồn trả lại Vu Gia sau 16 km).
• Cấp độ 3 (Vi mô - Nút dùng nước/Công trình): Mô hình hóa chi tiết nhu cầu của từng đối tượng: trạm bơm Cầu Đỏ, hệ thống thủy nông An Trạch, các cụm thủy nông Điện Bàn, Duy Xuyên, Đại Lộc và 4 hồ chứa thủy điện trọng điểm (A Vương, Sông Bung 4, Đăk Mi 4, Sông Tranh 2).

Quy trình nghiên cứu rigorous

• Thu thập dữ liệu và khảo sát: Dữ liệu KTTV thu thập từ Trung tâm Dự báo KTTV quốc gia chuỗi 1977 - 2016. Dữ liệu không gian địa hình GIS tỷ lệ 1/50.000 từ Cục Đo đạc, Bản đồ và Thông tin địa lý Việt Nam. Điều tra khảo sát thực địa trực tiếp thu thập dữ liệu vận hành công trình và phỏng vấn chuyên gia từ ngày 15-30/7/2017 (thuộc Đề tài cấp Bộ TNMT.08).
• Tam giác hóa phương pháp (Triangulation): Kết hợp kiểm chứng chéo giữa số liệu quan trắc thực tế trạm thủy văn Ái Nghĩa, Nông Sơn, Giao Thủy; số liệu mô phỏng từ Mike NAM; và kết quả cân bằng nước trên WEAP.
• Độ tin cậy và kiểm định mô hình: Bộ thông số mô hình Mike NAM được hiệu chỉnh và kiểm định tại 2 trạm kiểm soát chính: trạm Thành Mỹ (sông Vu Gia, diện tích khống chế 1.850 km²) và trạm Nông Sơn (sông Thu Bồn, diện tích khống chế 3.150 km²). Chỉ số Nash-Sutcliffe (NSE) đạt trên 0,85; sai số lưu lượng dòng chảy đỉnh và tổng lượng dòng chảy ($E_{rel}$) đều dưới 10%, khẳng định độ tin cậy tuyệt đối của công cụ toán học mô phỏng.

Data và phân tích

Dữ liệu đầu vào phản ánh sự biến thiên sâu sắc về khí hậu và thủy văn:

• Lượng mưa bình quân năm đạt trên 2.000 mm, nhưng tập trung cực đoan tại vùng núi cao như Trà My, Tiên Phước (>4.000 mm/năm), trong khi vùng đồng bằng chỉ đạt 2.000 - 2.500 mm/năm.
• Mùa mưa (tháng IX - XII) chiếm 65 - 80% tổng lượng mưa năm; riêng tháng X và XI chiếm 40 - 50%. Mùa khô (tháng I - VIII) chỉ chiếm 20 - 35%, đặc biệt 3 tháng kiệt nhất (tháng II - IV) lượng mưa chỉ chiếm 3 - 5% cả năm.
• Lưu lượng dòng chảy mùa lũ trung bình chiếm tới 63% tổng lượng dòng chảy năm. Lưu lượng đỉnh lũ lịch sử tại Nông Sơn đạt 5.600 m³/s, Thành Mỹ đạt 3.778 m³/s (năm 1998). Ngược lại, vào mùa kiệt, lưu lượng tháng kiệt nhất (tháng IV) chỉ chiếm 2 - 3% tổng lượng dòng chảy năm; lưu lượng thấp nhất lịch sử ghi nhận tại Thành Mỹ vào tháng 9/1988 chỉ còn 11,5 m³/s.

Quy trình phân tích sử dụng công cụ WEAP thiết lập mạng lưới nút phân bổ, mô phỏng cân bằng cung - cầu theo 4 kịch bản tổ hợp:

1. Hiện trạng năm 2014 ứng với tần suất dòng chảy $P=85%$
2. Hiện trạng năm 2014 ứng với tần suất $P=95%$
3. Tương lai năm 2030 ứng với tần suất $P=85%$
4. Tương lai năm 2030 ứng với tần suất $P=95%$

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

┌──────────────────────────────────────────────────────────────────────────────────────────────────┐
│                                 TỔNG HỢP PHÁT HIỆN ĐỘT PHÁ CỦA NGHIÊN CỨU                        │
└──────────────────────────────────────────────────────────────────────────────────────────────────┘
                                                 │
 ┌───────────────────────────────────────────────┼───────────────────────────────────────────────┐
 ▼                                               ▼                                               ▼
BẤT CÂN ĐỐI TỰ NHIÊN CỰC ĐOAN           ĐIỂM NGHẼN CHUYỂN DÒNG LIÊN LƯU VỰC     GIẢI PHÁP PHÂN BỔ TỐI ƯU
- Mùa mưa (IX-XII): 63%-68% nước        - Thủy điện Đăk Mi 4 chuyển nước từ      - Phân bổ động theo $W_{pb}$
- Tháng IV kiệt nhất: 2%-3%               Vu Gia sang Thu Bồn                    - Bảo đảm 100% sinh hoạt ($W_{ty}$)
- Mặn 1‰ ăn sâu 15-25 km                - Gây thiếu nước trầm trọng hạ lưu       - Ràng buộc dòng chảy tối thiểu ($W_{tt}$)
- Đe dọa trạm bơm Cầu Đỏ (Đà Nẵng)        Vu Gia (Ái Nghĩa, Cầu Đỏ)              - Giảm thiểu thiệt hại nông nghiệp

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

1. Sự bất cân đối thủy văn cực đoan theo không gian và thời gian: Nghiên cứu chứng minh dòng chảy năm trên lưu vực có mô đun dao động rất mạnh từ dưới 30 l/s.km² ở đồng bằng ven biển đến trên 90 l/s.km² tại thượng nguồn Trà My. Trích dẫn nguyên văn dữ liệu từ luận án: "Tổng lượng dòng chảy năm trung bình thời kỳ 1977 - 2016 của hệ thống sông vào khoảng 24 km3 (24 tỷ m3), chiếm 2,6% tổng lượng dòng chảy năm của các sông suối của nước ta; nếu chỉ tính dòng chảy nội địa thì chiếm khoảng 6,8%". Mùa lũ chiếm tới 63% tổng lượng dòng chảy, trong khi mùa khô chiếm 37%, riêng tháng IV dòng chảy suy kiệt chỉ còn 2 - 3% tổng lượng năm.
2. Lượng hóa tác động tiêu cực của chuyển dòng thủy điện: Việc vận hành hồ chứa thủy điện Đăk Mi 4 chuyển nước từ dòng Vu Gia sang dòng Thu Bồn đã làm suy giảm nghiêm trọng dòng chảy kiệt về hạ lưu sông Vu Gia tại Ái Nghĩa, khiến trạm bơm cấp nước Cầu Đỏ (Đà Nẵng) thường xuyên bị tê liệt do độ mặn vượt ngưỡng cho phép (>1.000 mg/l) trong nhiều tuần liền.
3. Phát hiện điểm nghẽn thiếu hụt trong kịch bản tương lai 2030: Dưới kịch bản biến đổi khí hậu kết hợp gia tăng dân số (dự báo đạt trên 2,2 triệu dân) và mở rộng công nghiệp (lấp đầy >70% diện tích 11 KCN), nếu xảy ra năm hạn hán nghiêm trọng ứng với tần suất $P=95%$, tổng lượng nước thiếu hụt tại các tiểu vùng hạ lưu lên tới hàng trăm triệu m³, trong đó nông nghiệp chịu rủi ro thiếu hụt trên 35% nhu cầu nếu duy trì cách thức phân bổ hành chính cũ.
4. Xác lập ma trận thứ tự ưu tiên phân bổ tối ưu: Luận án chứng minh rằng phương án phân bổ đạt điểm đa tiêu chí cao nhất là phương án tuân thủ nghiêm ngặt thứ bậc ưu tiên: (1) Cấp nước sinh hoạt và nhu cầu thiết yếu ($W_{ty}$); (2) Duy trì dòng chảy tối thiểu ($W_{tt}$) kết hợp đẩy mặn; (3) Nông nghiệp trồng lúa và hoa màu; (4) Công nghiệp và dịch vụ; (5) Phát điện thủy điện thương mại.

Implications đa chiều

• Về mặt lý thuyết: Cung cấp bằng chứng thực nghiệm khẳng định tính tương thích của mô hình WEAP khi tích hợp phương trình khấu trừ dòng chảy sinh thái ($W_{tt}$) trong các lưu vực sông duyên hải miền Trung có độ dốc lớn và mùa khô kéo dài.
• Về phương pháp luận: Chuẩn hóa quy trình 4 bước tính toán lượng nước có thể phân bổ ($W_{pb}$) kết hợp đa giác Thiessen và mô hình mưa rào - dòng chảy Mike NAM, tạo tiền lệ phương pháp cho 10 lưu vực sông lớn còn lại của Việt Nam.
• Về ứng dụng thực tiễn: Cung cấp biểu đồ điều tiết và lịch xả nước chi tiết cho 4 hồ chứa lớn (A Vương, Sông Bung 4, Đăk Mi 4, Sông Tranh 2), giúp các công ty thủy điện và Công ty Cổ phần Cấp nước Đà Nẵng (Dawaco) chủ động phối hợp điều tiết nước đẩy mặn.
• Về chính sách và quản trị: Cung cấp cơ sở khoa học trực tiếp để Bộ Tài nguyên và Môi trường hoàn thiện Quy trình vận hành liên hồ chứa lưu vực sông Vu Gia - Thu Bồn, đồng thời đề xuất thành lập Ủy ban Lưu vực Sông Vu Gia - Thu Bồn nhằm thể chế hóa cơ chế chia sẻ trách nhiệm liên tỉnh.

Limitations và Future Research

Nghiên cứu có 3 hạn chế chính:

1. Chưa tích hợp mô hình dòng chảy ngầm 3 chiều hoàn chỉnh: Nghiên cứu chủ yếu tập trung vào nguồn nước mặt; mối tương tác động lực học giữa nước mặt và các tầng chứa nước dưới đất chưa được liên kết trực tiếp bằng các mô hình chuyên sâu như MODFLOW hay GMS.
2. Độ phân giải kịch bản biến đổi khí hậu: Các kịch bản BĐKH áp dụng chủ yếu dựa trên biến động lượng mưa và nhiệt độ trung bình tháng, chưa chi tiết hóa đến mức độ cực đoan theo ngày (daily extreme events) của các đợt hạn hán kéo dài trên 60 ngày.
3. Cơ chế kinh tế thị trường nước (Water Trading): Luận án tập trung vào phân bổ hạn ngạch kỹ thuật và điều tiết hành chính kết hợp chia sẻ chi phí, chưa xây dựng được mô hình định giá kinh tế nước biên (marginal water pricing) để hình thành sàn giao dịch quyền sử dụng nước như lưu vực Murray-Darling.

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

• Phát triển mô hình tích hợp hoàn toàn bề mặt - ngầm (Coupled Surface-Groundwater Modeling: MIKE SHE hoặc WEAP-MODFLOW).
• Nghiên cứu cơ chế chi trả dịch vụ môi trường rừng và nước (PES/WES) giữa thành phố Đà Nẵng (vùng hạ du hưởng lợi) và tỉnh Quảng Nam (vùng thượng nguồn sinh thủy).
• Tích hợp công nghệ cảm biến IoT quan trắc độ mặn và lưu lượng thời gian thực để tự động hóa thuật toán điều tiết xả tràn liên hồ chứa.

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

• Học thuật: Công trình đóng góp hệ thống dữ liệu thủy văn thực nghiệm 40 năm chuẩn hóa, tạo nền tảng trích dẫn quan trọng cho các nghiên cứu thủy văn học, biến đổi khí hậu và quy hoạch tài nguyên nước tại Việt Nam và Đông Nam Á.
• Kinh tế - Ngành: Giúp tối ưu hóa cân bằng giữa doanh thu phát điện của EVN (hàng nghìn tỷ đồng mỗi năm) và an toàn nguồn nước phục vụ 40.332 ha lúa và hoa màu tại đồng bằng Quảng Nam, giảm thiểu nguy cơ thiệt hại kinh tế do hạn hán hàng trăm tỷ đồng mỗi năm.
• Chính sách: Đóng góp luận cứ khoa học then chốt phục vụ việc sửa đổi Luật Tài nguyên nước và ban hành các quy định hướng dẫn về phân bổ tài nguyên nước lưu vực sông của Bộ Tài nguyên và Môi trường.
• Xã hội và Môi trường: Đảm bảo quyền tiếp cận nước sạch sinh hoạt bền vững cho hơn 1,5 triệu dân cư hiện tại và 2,2 triệu dân trong tương lai; bảo vệ hệ sinh thái rừng ngập mặn cửa sông Cửa Đại và vịnh Đà Nẵng.

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

┌───────────────────────────────────────────────────────────────────────────────────────────────────┐
│                                     ĐỐI TƯỢNG HƯỞNG LỢI TRỰC TIẾP                                 │
└───────────────────────────────────────────────────────────────────────────────────────────────────┘
                                                  │
 ┌───────────────────────┬────────────────────────┼────────────────────────┬───────────────────────┐
 ▼                       ▼                        ▼                        ▼                       ▼
NGHIÊN CỨU SINH &       GIẢNG VIÊN &             DOANH NGHIỆP CẤP NƯỚC    CHỦ HỒ THỦY ĐIỆN        CƠ QUAN QUẢN LÝ
HỌC VIÊN CAO HỌC        NHÀ KHOA HỌC             & CÔNG NGHIỆP (Dawaco)   (EVN, Đăk Mi, Sông Bung)(Bộ TN&MT, UBND 2 tỉnh)
- Mô hình Mike NAM-WEAP - Cơ sở lý thuyết IWRM   - Chủ động lịch cấp      - Lịch xả tối ưu        - Ra quyết định phân bổ
- Phương pháp luận $W_{pb}$- Bộ số liệu chuẩn 40 năm - Kiểm soát rủi ro mặn - Cân bằng phát điện    - Quy trình liên hồ

• Nghiên cứu sinh và Giới học thuật: Tiếp cận khung phương pháp luận kết hợp Mike NAM và WEAP trong giải quyết bài toán thủy văn - kinh tế phức tạp.
• Cơ quan Hoạch định Chính sách (Bộ TN&MT, UBND TP Đà Nẵng, UBND tỉnh Quảng Nam): Sử dụng trực tiếp ma trận phân bổ và các kịch bản $P=85%$, $P=95%$ để chỉ đạo điều hành mùa kiệt.
• Đơn vị Quản lý Khai thác Nước (Công ty Cấp nước Đà Nẵng - Dawaco, các Công ty Khai thác CTTL Quảng Nam): Có cơ sở dự báo và lịch điều tiết nước đẩy mặn chủ động.
• Chủ đầu tư Thủy điện (EVN, Thủy điện Đăk Mi, A Vương, Sông Bung, Sông Tranh): Nắm bắt khung vận hành pháp lý và trách nhiệm xả nước duy trì dòng chảy tối thiểu, tránh xung đột pháp lý và 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ì và đã mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc mở rộng Lý thuyết Quản lý Tổng hợp Tài nguyên Nước (IWRM) thông qua việc chuẩn hóa toán học phương trình lượng nước có thể phân bổ thực tế:

$$W_{pb} = W_{sd} - (W_{tt} + W_{ty})$$

Khác với các nghiên cứu IWRM định tính, luận án lượng hóa $W_{tt}$ (dòng chảy tối thiểu) và $W_{ty}$ (nhu cầu thiết yếu) thành các điều kiện biên cứng bắt buộc phải trừ khỏi tổng lượng nước sử dụng ($W_{sd}$) trước khi phân bổ kinh tế, biến lý thuyết IWRM thành công cụ tính toán số học cụ thể trên nền tảng lưu vực sông chịu tác động chuyển dòng.

2. Đổi mới phương pháp luận của luận án so với các công trình trước đây tại Việt Nam?

So với nghiên cứu của Tô Trung Nghĩa (2008) hay Nguyễn Chí Công (2014) vốn sử dụng mô hình GAMS tĩnh hoặc bảng tính Excel cân bằng khối lượng đơn giản, luận án đã tạo bước đột phá khi:

• Tích hợp không gian hóa lượng mưa bằng phương pháp đa giác Thiessen.
• Chuyển đổi mưa thành dòng chảy biên thông qua mô hình thủy văn Mike NAM (đạt NSE > 0,85 tại Thành Mỹ và Nông Sơn).
• Nhúng toàn bộ mạng lưới thủy lực phức tạp (kể cả kênh phân lưu Quảng Huế và sông Vĩnh Điện) vào mô hình động WEAP để giải bài toán tối ưu hóa phân bổ theo chuỗi thời gian thực đa kịch bản.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì?

Phát hiện bất ngờ và có ý nghĩa cảnh báo cao nhất là vào mùa kiệt, tháng có dòng chảy nhỏ nhất không chỉ rơi vào tháng IV theo quy luật tự nhiên (chỉ chiếm 2 - 3% lượng nước cả năm), mà trong nhiều năm chuyển tiếp (như tháng 9/1988), lưu lượng tại Thành Mỹ đã tụt xuống mức kỷ lục 11,5 m³/s do ảnh hưởng kết hợp của hạn khí tượng và chu kỳ tích nước sớm của các hồ thủy điện thượng nguồn. Điều này chứng minh quy trình tích nước mùa thu của các hồ chứa nếu không được kiểm soát sẽ tạo ra "hạn nhân tạo" ngay trước mùa lũ.

                  ┌─────────────────────────────────────────────────────────┐
                  │ ĐỘT PHÁ CỦA CÔNG TRÌNH: QUY TRÌNH TÁI LẬP (REPLICATION) │
                  └───────────────────────────┬─────────────────────────────┘
                                              │
    ┌─────────────────────────────────────────┼─────────────────────────────────────────┐
    ▼                                         ▼                                         ▼
CHẾ ĐỘ MƯA THEO THIESSEN               MÔ HÌNH MIKE NAM HIỆU CHỈNH               PHÂN BỔ TRÊN WEAP
- Bản đồ mưa GIS 1/50.000             - Bộ tham số chuẩn hóa 40 năm             - 5 cấp ưu tiên phân bổ
- Trọng số trạm mưa tự động           - Hiệu chuẩn tại 2 trạm kiểm soát         - Giải bài toán hạn hán P=85%, 95%

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?

Luận án cung cấp đầy đủ và minh bạch quy trình tái lập bao gồm:

• Toàn bộ bảng thông số hình thái lưu vực và mạng lưới sông suối.
• Bộ thông số hiệu chỉnh mô hình Mike NAM cho 2 trạm Nông Sơn và Thành Mỹ.
• Tọa độ địa lý và lưu lượng thiết kế của các điểm phân bổ.
• Cấu trúc cây nhu cầu nước và ma trận ưu tiên cài đặt trong phần mềm WEAP. Bất kỳ nhà nghiên cứu nào có dữ liệu mưa - dòng chảy tương ứng đều có thể tái lập hoàn toàn kết quả nghiên cứu.

5. Chương trình nghiên cứu 10 năm (10-year Research Agenda) được định hình như thế nào?

Chương trình nghiên cứu 10 năm tới bao gồm:

• Xây dựng hệ thống Digital Twin (Bản sao số thủy văn) toàn diện lưu vực Vu Gia - Thu Bồn.
• Nghiên cứu ảnh hưởng của nước biển dâng và triều cường cực đoan đến biến động xâm nhập mặn vùng hạ lưu.
• Thiết lập thị trường trao đổi quyền sử dụng nước thí điểm giữa ngành nông nghiệp và công nghiệp - dịch vụ tại thành phố Đà Nẵng và tỉnh Quảng Nam.

Kết luận

Công trình nghiên cứu tiến sĩ của tác giả Thân Văn Đón tạo nên một dấu ấn khoa học thực chứng quan trọng trong chuyên ngành Thủy văn học và Quản lý tài nguyên nước tại Việt Nam với 5 đóng góp cốt lõi:

1. Chuẩn hóa cơ sở lý luận và phương pháp luận: Xác lập quy trình phân bổ hợp lý nguồn nước mặt liên tỉnh dựa trên công thức định lượng $W_{pb} = W_{sd} - (W_{tt} + W_{ty})$, loại bỏ tính chủ quan trong phân bổ hạn ngạch hành chính.
2. Làm chủ công nghệ mô hình hóa tích hợp: Ứng dụng thành công liên kết Mike NAM - WEAP với chuỗi số liệu 40 năm, đạt độ chính xác cao (NSE > 0,85), mô phỏng chuẩn xác động lực dòng chảy phức tạp của hệ thống sông có phân lưu và chuyển dòng thủy điện.
3. Lượng hóa toàn diện các kịch bản cực đoan: Đánh giá định lượng mức độ thiếu hụt nguồn nước mặt trong hiện trạng (2014) và tương lai (2030) ứng với các tần suất hạn hán gay gắt $P=85%$ và $P=95%$.
4. Đề xuất giải pháp phân bổ tối ưu và khả thi: Thiết lập ma trận thứ tự ưu tiên 5 bậc kết hợp tỷ lệ cắt giảm thích ứng, bảo đảm an ninh cấp nước sinh hoạt cho hơn 1,5 - 2,2 triệu dân và bảo tồn hệ sinh thái hạ du.
5. Định hình cơ chế thể chế liên tỉnh: Đặt nền móng khoa học và pháp lý cho việc thành lập Ủy ban Lưu vực Sông Vu Gia - Thu Bồn và hoàn thiện Quy trình vận hành liên hồ chứa của Chính phủ.

Thành công của nghiên cứu mở ra hướng tiếp cận chuẩn mực cho công tác quy hoạch tổng hợp tài nguyên nước tại các lưu vực sông xuyên tỉnh trên toàn lãnh thổ Việt Nam.