Tổng quan về luận án

Nghiên cứu của tác giả Nguyễn Đại Trung (2024) với đề tài "Nghiên cứu đánh giá mức độ căng thẳng nguồn nước cho thành phố Đà Nẵng" thuộc chuyên ngành Kỹ thuật Tài nguyên nước (Mã số: 9580212), được thực hiện dưới sự hướng dẫn khoa học của TS. Nguyễn Anh Đức và GS.TS. Nguyễn Trung Việt tại Trường Đại học Thủy Lợi, đánh dấu một bước tiến quan trọng trong việc định lượng rủi ro an ninh nguồn nước tại các đô thị ven biển Việt Nam. Trong bối cảnh đô thị hóa nhanh và biến đổi khí hậu diễn biến phức tạp, các trung tâm kinh tế biển đối mặt với bài toán nan giải: tổng lượng nước tự nhiên dồi dào nhưng mức độ căng thẳng nguồn nước (Water Stress - WSI) sinh hoạt và sản xuất lại tăng đột biến vào mùa cạn.

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 nghiên cứu tài nguyên nước trước đây tại Việt Nam (Nguyễn Mai Đăng và cs., 2016; Viện Khoa học Thủy lợi Việt Nam, 2013) chủ yếu tiếp cận đơn lẻ theo góc độ thủy văn thuần túy, cân bằng nước vật lý hoặc đánh giá chất lượng nước độc lập mà thiếu vắng một khung phân tích tích hợp đa chiều. Các công cụ truyền thống chưa lượng hóa được mối quan hệ hữu cơ giữa áp lực khai thác thượng nguồn, diễn biến xâm nhập mặn cửa sông, yêu cầu dòng chảy môi trường và năng lực ứng phó thể chế của một đô thị hạ du.

Để giải quyết khoảng trống này, luận án xây dựng 3 câu hỏi nghiên cứu (RQs) và 3 luận điểm bảo vệ cốt lõi:

  1. RQ1: Làm thế nào để xây dựng một khung và bộ chỉ số định lượng toàn diện mức độ căng thẳng nguồn nước đặc thù cho đô thị ven biển hạ lưu sông liên tỉnh?
  2. RQ2: Mức độ căng thẳng nguồn nước tại Đà Nẵng biến động như thế nào giữa hiện trạng năm 2020 và tương lai năm 2030 dưới tác động kép của biến đổi khí hậu (kịch bản RCP4.5) và gia tăng dân số?
  3. RQ3: Những kịch bản can thiệp công trình và phi công trình nào có khả năng tối ưu hóa chỉ số WSI nhằm bảo đảm an ninh cấp nước bền vững?

Nghiên cứu được neo giữ trên nền tảng lý thuyết Quản lý Tổng hợp Tài nguyên Nước (IWRM) và Khung phân tích Tính dễ bị tổn thương sinh thái - xã hội (Socio-Ecological Vulnerability Framework). Về mặt phạm vi, luận án khảo sát toàn diện không gian hạ lưu lưu vực sông Vu Gia - Thu Bồn (diện tích lưu vực tính đến Ái Nghĩa đạt $5.180\text{ km}^2$, tổng lượng dòng chảy năm $W_0 = 8,55 \times 10^9\text{ m}^3$) và lưu vực sông Cu Đê ($472\text{ km}^2$), tập trung vào 8 đơn vị hành chính của thành phố Đà Nẵng (diện tích tự nhiên 128.488 ha, dân số năm 2020 là 1.169,5 nghìn người với 87,3% dân số đô thị). Mốc thời gian đánh giá bao phủ chuỗi số liệu thủy văn lịch sử 1980–2016 và dự báo động thái cân bằng nước giai đoạn 2016–2035.

Đóng góp đột phá của công trình nằm ở việc phát triển thành công bộ chỉ số WSI gồm 4 nhóm và 25 chỉ số thành phần, được lượng hóa qua chuỗi mô hình toán học tích hợp (MIKE NAM, MIKE HYDRO BASIN, MIKE 11 HD/AD). Kết quả giúp làm rõ hiện tượng "nghịch lý khan hiếm trong thừa thãi" tại Đà Nẵng, cung cấp bằng chứng định lượng chính xác để chuyển dịch từ quản lý ứng phó thụ động sang kiểm soát rủi ro chủ động theo Kết luận số 36/KL-TW của Bộ Chính trị.

Literature Review và Positioning

Tổng quan y văn quốc tế cho thấy sự tiến hóa của các chỉ số đo lường tài nguyên nước trải qua bốn thập kỷ. Khởi đầu từ chỉ số Falkenmark (1989) lấy ngưỡng lượng nước sẵn có bình quân đầu người ($>1.700\text{ m}^3/\text{người/năm}$ là không căng thẳng, $1.000 - 1.700\text{ m}^3$ là căng thẳng, và $<500\text{ m}^3$ là khan hiếm tuyệt đối), các nghiên cứu ban đầu chỉ phản ánh khía cạnh thủy văn vật lý tĩnh. Vörösmarty và cs. (2000) mở rộng đánh giá quy mô toàn cầu khi kết hợp tổng lượng dòng chảy mặt với nhu cầu tiêu thụ, chỉ ra sức ép gia tăng dân số đến năm 2025. Tuy nhiên, các cách tiếp cận này bộc lộ mâu thuẫn học thuật lớn: xem con người là trung tâm khai thác mà xem nhẹ chức năng bảo tồn hệ sinh thái, đồng thời bỏ qua sự biến động theo mùa và ranh giới mặn - ngọt ở vùng cửa sông ven biển (White, 2012; Damkjaer & Taylor, 2017).

Trong bối cảnh đô thị, hai luồng quan điểm đối lập định hình các mô hình đánh giá hiện đại:

  • Luồng quan điểm thứ nhất (Kỹ thuật - Cân bằng cung cầu vật lý): Đại diện bởi các nghiên cứu của OECD (2003) và Smakhtin và cs. (2004), tập trung vào tỷ lệ khai thác trên tổng lượng nước sẵn có ($WSI = Q_{\text{khai thác}} / Q_{\text{tổng}}$). Khi tỷ lệ vượt $40%$, hệ thống rơi vào trạng thái căng thẳng nghiêm trọng. Nhược điểm của luồng quan điểm này là không tính đến năng lực trữ nước, hiệu quả phân phối hay sự gián đoạn chất lượng do nhiễm bẩn hoặc nhiễm mặn.
  • Luồng quan điểm thứ hai (Đa chiều sinh thái - thể chế): Đại diện bởi các khung đánh giá an ninh nguồn nước đô thị (UWSI). Tiêu biểu như nghiên cứu của Firdaus Ali (2012) tại Jakarta (Indonesia) thiết lập khung 3 thành phần với 8 chỉ số, nhưng tồn tại hạn chế là bỏ qua động thái xâm nhập mặn và xem nhẹ trọng số của từng chỉ số thành phần. Tương tự, nghiên cứu của Al-Nozaily và cs. (2022) tại thành phố Ibb (Yemen) phát triển 4 khía cạnh với 24 chỉ số nhưng gán trọng số đồng đều, không phản ánh được tính nhạy cảm đặc thù của vùng địa hình ven biển có nguồn nước phụ thuộc vào vận hành liên hồ chứa bậc thang thượng nguồn.
TIẾN HÓA CỦA CÁC KHUNG ĐÁNH GIÁ NGUỒN NƯỚC
Falkenmark (1989) ──> Vörösmarty (2000) ──> Firdaus Ali (2012) / Ibb (2022) ──> WSI Đà Nẵng (2024)
(Chỉ số vật lý 1D)    (Thủy văn - Dân số)      (UWSI đa chiều, không trọng số)   (Tích hợp 4 nhóm 25 chỉ số + AHP/Delphi/MIKE)

Tại Việt Nam, các nghiên cứu tại TP. Hồ Chí Minh (2016) chỉ áp dụng các công thức WSI đơn lẻ nên dẫn đến kết luận mâu thuẫn: không căng thẳng nếu tính tổng năm toàn lưu vực, nhưng lại khan hiếm gay gắt vào mùa khô (tháng 2–4 ứng với tần suất $P=95%$). Nghiên cứu tại Hà Nội của Nguyễn Mai Đăng và cs. (2016) tiếp cận theo an ninh nguồn nước gồm 7 khía cạnh nhưng thiếu vắng cấu phần thủy lực xâm nhập mặn.

Luận án của Nguyễn Đại Trung định vị chính xác vào điểm giao thoa còn khuyết: Xây dựng khung WSI chuẩn hóa cho đô thị duyên hải có nguồn nước phụ thuộc vào hạ du sông liên tỉnh (hệ thống sông Vu Gia - Hàn). Nghiên cứu so sánh trực tiếp và khắc phục triệt để các khiếm khuyết của mô hình Jakarta (Firdaus Ali, 2012) và mô hình Ibb (Al-Nozaily và cs., 2022) bằng cách bổ sung chỉ số lan truyền mặn thực đo kết hợp mô phỏng thủy lực 1D/AD, đồng thời tích hợp dòng chảy môi trường theo phương pháp Tennant (1976) và gán trọng số khoa học thông qua ma trận so sánh cặp AHP kiểm định qua thuật toán Delphi-KAMET.

Đó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 (Integrated Water Resources Management - IWRM) và Lý thuyết Khả năng Phục hồi Đô thị (Urban Resilience Theory) thông qua 4 luận điểm lý thuyết đột phá:

  1. Tái cấu trúc ranh giới khái niệm Căng thẳng nguồn nước (Water Stress): Luận án chứng minh rằng trạng thái căng thẳng nguồn nước đô thị không đồng nhất với "khan hiếm nước vật lý" (Water Scarcity). Bằng chứng thực nghiệm tại Đà Nẵng cho thấy: mặc dù tổng lượng dòng chảy mặt lưu vực sông Vu Gia - Thu Bồn đạt tới $8,55 \times 10^9\text{ m}^3/\text{năm}$, thành phố vẫn chịu căng thẳng nguồn nước nghiêm trọng vào mùa kiệt do suy giảm chất lượng nước (nhiễm mặn) và xung đột điều tiết thủy điện thượng lưu.
  2. Lý thuyết hóa mối quan hệ bất đối xứng Cung - Cầu qua thời gian và không gian: Thiết lập mệnh đề lý thuyết cho thấy sự tập trung hóa công trình thu nước tại hạ du (như NMN Cầu Đỏ chiếm 89% công suất toàn đô thị) làm gia tăng độ nhạy cảm tổn thương của hệ sinh thái đô thị trước biến động lưu lượng xả kiệt từ các hồ chứa thủy điện bậc thang thượng nguồn (A Vương, Đak Mi 4, Sông Bung 4, Sông Tranh 2).
  3. Định lượng hóa yêu cầu sinh thái hạ du trong cấu trúc WSI: Tích hợp dòng chảy môi trường tối thiểu (Tennant, 1976) như một biến số ràng buộc bắt buộc, chấm dứt quan điểm xem dòng chảy hạ du ra biển là lượng nước dư thừa.
  4. Chuyển dịch hình thái đánh giá (Paradigm Shift): Chuyển đổi mô hình đánh giá từ tĩnh (static, dựa trên thống kê bình quân năm) sang mô hình động (dynamic, dựa trên chuỗi thời gian phân giải cao và mô phỏng thủy lực kịch bản).

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

Khung phân tích WSI được cấu trúc hóa theo 4 nhóm trụ cột ($WSI_1 \to WSI_4$) bao hàm 25 chỉ số thành phần độc lập nhưng tương tác phi tuyến:

$$\text{WSI}{\text{Tổng hợp}} = \sum{i=1}^{4} \mathbf{W}i \cdot \text{WSI}i = \sum{i=1}^{4} \mathbf{W}i \left( \sum{j=1}^{n_i} w{ij} \cdot I_{ij} \right)$$

Trong đó: $\mathbf{W}i$ là trọng số nhóm cấp 1; $w{ij}$ là trọng số chỉ số thành phần cấp 2; $I_{ij}$ là giá trị chuẩn hóa của chỉ số $j$ thuộc nhóm $i$ trên thang đo $[1, 5]$.

                    KHUNG ĐÁNH GIÁ CĂNG THẲNG NGUỒN NƯỚC ĐÔ THỊ (WSI)
                                            │
    ┌───────────────────────┬───────────────┴───────────────┬───────────────────────┐
    ▼                       ▼                               ▼                       ▼
┌──────────────────┐    ┌──────────────────┐    ┌──────────────────┐    ┌──────────────────┐
│  Nhóm WSI_1      │    │  Nhóm WSI_2      │    │  Nhóm WSI_3      │    │  Nhóm WSI_4      │
│  Nguồn nước &    │    │  Hệ sinh thái &  │    │  Cấp nước SH từ  │    │  Năng lực ứng    │
│  KTSDN           │    │  Môi trường      │    │  CTCNTT          │    │  phó thể chế     │
│  (Trọng số W1)   │    │  (Trọng số W2)   │    │  (Trọng số W3)   │    │  (Trọng số W4)   │
└────────┬─────────┘    └────────┬─────────┘    └────────┬─────────┘    └────────┬─────────┘
         │                       │                       │                       │
         ├─ Lượng nước sẵn có    ├─ Dòng chảy MT         ├─ Tỷ lệ cấp nước SH    ├─ Dung tích trữ
         ├─ Biến động Cv         ├─ Độ che phủ rừng      ├─ Thất thoát nước      ├─ Thể chế chính sách
         ├─ Tỷ lệ khai thác      └─ Suy thoái HST        ├─ Tần suất nhiễm mặn   ├─ Quan trắc cảnh báo
         └─ Chất lượng VN_WQI                            └─ Ổn định nguồn cấp    └─ Kế hoạch ứng phó

Các điều kiện biên của khung phân tích:

  • Về tài nguyên: Chỉ tính toán cho nguồn nước mặt; không tính nước ngầm do Đà Nẵng đã dừng cấp phép khai thác từ ngày 01/01/2015 và tiến hành trám lấp giếng khoan theo Quyết định số 7517/QĐ-UBND nhằm chống sụt lún và sụt giảm tầng chứa nước Holocen ($qh$) và Pleistocen ($qp$).
  • Về thủy lực: Biên động học phụ thuộc hoàn toàn vào chế độ bán nhật triều không đều tại vịnh Đà Nẵng (mực nước triều cao nhất $+143\text{ cm}$, thấp nhất $-40\text{ cm}$) và 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.

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

Thiết kế nghiên cứu

Luận án tuân thủ lập trường triết học Thực chứng định lượng (Positivism) kết hợp Hiện thực phản biện (Critical Realism), kết hợp đa phương pháp (Mixed-methods) với thiết kế đa tầng bậc (Multi-level design):

  • Cấp độ 1 - Lưu vực liên tỉnh: Phân tích tương tác hoàn lưu khí quyển, lượng mưa, dòng chảy lũ - kiệt trên diện tích $10.350\text{ km}^2$ của hệ thống sông Vu Gia - Thu Bồn.
  • Cấp độ 2 - Đô thị hạ du: Đánh giá mạng lưới thủy lực phân lưu tại đập dâng An Trạch, sông Yên, sông Cầu Đỏ, sông Cẩm Lệ, sông Vĩnh Điện và sông Cu Đê.
  • Cấp độ 3 - Điểm tiêu thụ: Khảo sát 04 hệ thống cấp nước đô thị tập trung (Cầu Đỏ - Sân Bay công suất $267.000\text{ m}^3/\text{ngày}$, Hải Vân $3.600\text{ m}^3/\text{ngày}$, Sơn Trà $6.900\text{ m}^3/\text{ngày}$, Hòa Trung $10.000\text{ m}^3/\text{ngày}$) và các trạm cấp nước nông thôn Hòa Vang.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu được thực hiện qua các bước chuẩn hóa nghiêm ngặt:

                            QUY TRÌNH NGHIÊN CỨU RIGOROUS
                                          │
┌─────────────────────────────────────────┴─────────────────────────────────────────┐
│ Bước 1: Khảo sát thực địa & Thu thập dữ liệu KTTV (Chuỗi số liệu 1980 - 2020)     │
└─────────────────────────────────────────┬─────────────────────────────────────────┘
                                          │
┌─────────────────────────────────────────┴─────────────────────────────────────────┐
│ Bước 2: Sàng lọc chỉ số qua phương pháp Delphi & Quy tắc KAMET (Hội đồng chuyên gia)│
└─────────────────────────────────────────┬─────────────────────────────────────────┘
                                          │
┌─────────────────────────────────────────┴─────────────────────────────────────────┐
│ Bước 3: Xác lập ma trận so sánh cặp AHP & Kiểm định tỷ số nhất quán (CR < 0.1)   │
└─────────────────────────────────────────┬─────────────────────────────────────────┘
                                          │
┌─────────────────────────────────────────┴─────────────────────────────────────────┐
│ Bước 4: Thiết lập chuỗi mô hình toán học (MIKE NAM ──> HYDRO BASIN ──> MIKE 11)   │
└─────────────────────────────────────────┬─────────────────────────────────────────┘
                                          │
┌─────────────────────────────────────────┴─────────────────────────────────────────┐
│ Bước 5: Hiệu chuẩn & Kiểm định mô hình (Nash-Sutcliffe NSE > 0.80, % sai số < 10%) │
└─────────────────────────────────────────┬─────────────────────────────────────────┘
                                          │
┌─────────────────────────────────────────┴─────────────────────────────────────────┐
│ Bước 6: Mô phỏng kịch bản hiện trạng (2020) & Kịch bản BĐKH RCP4.5 (2030)         │
└───────────────────────────────────────────────────────────────────────────────────┘
  1. Sàng lọc chỉ số: Sử dụng kỹ thuật Delphi 2 vòng tham vấn 25 chuyên gia đầu ngành từ các viện nghiên cứu, trường đại học và cơ quan quản lý (Đại học Thủy lợi, Viện Quy hoạch Thủy lợi, Sở TN&MT Đà Nẵng), xử lý dữ liệu qua quy tắc sàng lọc KAMET để chọn lọc 25 chỉ số có độ tin cậy và tính khả thi dữ liệu cao nhất từ 38 chỉ số đề xuất ban đầu.
  2. Xác định trọng số: Ứng dụng quy trình Phân tích Thứ bậc (AHP - Saaty). Ma trận so sánh cặp được kiểm tra độ nhất quán nghiêm ngặt, đảm bảo Tỷ số nhất quán ($CR - Consistency\ Ratio$) đạt $CR < 0,10$.
  3. Đảm bảo độ tin cậy (Reliability & Validity): Triệt tiêu sai số qua phương pháp tam giác đạc dữ liệu (Data Triangulation), đối chiếu số liệu quan trắc trạm thủy văn quốc gia (Nông Sơn, Thạnh Mỹ, Ái Nghĩa) với dữ liệu quan trắc tự động của Công ty Cổ phần Cấp nước Đà Nẵng (DAWACO).

Data và phân tích

Nghiên cứu ứng dụng bộ công cụ mô hình toán thủy lực - thủy văn chuyên sâu của Viện Thủy lực Đan Mạch (DHI):

  • Mô hình MIKE NAM: Mô phỏng quá trình mưa - dòng chảy trên 14 tiểu lưu vực sông Vu Gia - Thu Bồn và 3 tiểu lưu vực sông Cu Đê. Kết quả hiệu chỉnh giai đoạn 1981–1996 và kiểm định giai đoạn 1997–2014 tại các trạm Thạnh Mỹ, Nông Sơn, Thượng Nhật cho thấy hệ số hiệu quả Nash-Sutcliffe ($NSE$) đạt từ $0,81$ đến $0,92$, sai số lượng dòng chảy tổng cộng ($PBIAS$) nằm trong ngưỡng $\pm 6,5%$, khẳng định độ tin cậy rất cao.
  • Mô hình MIKE HYDRO BASIN: Mô phỏng cân bằng nước toàn lưu vực, tích hợp dung tích hữu ích và quy trình vận hành xả kiệt của hệ thống thủy điện thượng nguồn (A Vương: $266,5 \times 10^6\text{ m}^3$, Đak Mi 4: $158 \times 10^6\text{ m}^3$, Sông Bung 4: $234 \times 10^6\text{ m}^3$, Sông Tranh 2: $521 \times 10^6\text{ m}^3$) cùng nhu cầu cấp nước chi tiết cho hạ du.
  • Mô hình MIKE 11 (HD & AD): Mô phỏng thủy động lực học dòng chảy và lan truyền mặn trên toàn mạng lưới sông hạ du. Mô hình phân giải chi tiết sự dịch chuyển của đường ranh mặn $1,0‰$ và $0,25‰$ (ngưỡng quy chuẩn QCVN 08-MT:2015/BTNMT). Kết quả kiểm định mực nước tại Ái Nghĩa và Cẩm Lệ đạt hệ số $NSE > 0,85$; sai số mô phỏng độ mặn cực đại tại cửa thu Cầu Đỏ nhỏ hơn $8,2%$.

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

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

  • Phát hiện 1 - Hiện tượng suy giảm nguồn nước cục bộ cực đoan: Số liệu thực nghiệm chỉ ra sự bất đối xứng nghiêm trọng giữa dòng chảy mùa lũ và mùa cạn. Dòng chảy 3 tháng mùa kiệt (tháng 3 đến tháng 5) chỉ chiếm $11,2 - 14,8%$ tổng lượng dòng chảy năm, đẩy lưu lượng tại trạm Ái Nghĩa có thời điểm xuống dưới $45\text{ m}^3/\text{s}$.
  • Phát hiện 2 - Xâm nhập mặn vượt ngưỡng lịch sử tại cửa thu nước chính: Luận án trích xuất minh chứng thực địa xác thực: "Sáng ngày 29/3/2021, mực nước sông Yên tại trạm bơm phòng mặn An Trạch hạ xuống mức +1,44m (rất thấp so với thiết kế là +2,00m) gây khó khăn cho việc vận hành bình thường của các máy bơm, độ mặn trên sông Cẩm Lệ tại cửa thu nước của nhà máy nước Cầu Đỏ tăng cao lên đến 4.846mg/l, gấp gần 20 lần so với ngưỡng an toàn theo quy định (250 mg/l)". Độ mặn cao nhất từng ghi nhận tại đây đạt tới $7.700\text{ mg/l}$, làm tê liệt khả năng cấp nước tự chảy trong nhiều ngày liên tiếp.
  • Phát hiện 3 - Căng thẳng nguồn nước diễn biến phân hóa mạnh theo không gian: Điểm đánh giá nhóm $WSI_1$ (Nguồn nước & Khai thác) và $WSI_3$ (Cung cấp nước sinh hoạt) tại khu vực trung tâm và hạ lưu sông Cầu Đỏ, sông Vĩnh Điện luôn ở mức Căng thẳng Cao (mức 4) đến Rất cao (mức 5), trong khi vùng thượng lưu Cu Đê ở mức Trung bình (mức 2–3).
  • Phát hiện 4 - Xu thế gia tăng căng thẳng nguồn nước đến năm 2030: Dưới kịch bản RCP4.5 kết hợp mức tăng nhu cầu nước đô thị bình quân $12%/\text{năm}$ (tương ứng mức tăng dân số $5,9%/\text{năm}$), nếu không có giải pháp can thiệp, chỉ số tổng hợp $WSI_{\text{Tổng hợp}}$ của toàn thành phố Đà Nẵng sẽ tăng từ $2,84$ (năm 2020 - mức Căng thẳng Trung bình) lên $3,62$ (năm 2030 - tiệm cận ngưỡng Căng thẳng Nghiêm trọng).
BIẾN ĐỘNG CHỈ SỐ WSI TỔNG HỢP QUA CÁC GIAI ĐOẠN VÀ KỊCH BẢN
  5 ┬───────────────────────────────────────────────────────────────
    │                                                    [Mức 4: Căng thẳng cao]
  4 ┼──────────────────────────────────── 3.62 ─────────────
    │                                   (2030 Không can thiệp)
  3 ┼────────── 2.84 ─────────────────────────────────── 2.45 ──────
    │    (2020 Hiện trạng)                           (2030 Có giải pháp)
  2 ┼───────────────────────────────────────────────────────────────
    │                                                    [Mức 2: Ít căng thẳng]
  1 ┴───────────────────────────────────────────────────────────────

Implications đa chiều

  • Hàm ý lý thuyết: Khẳng định tính ưu việt của mô hình WSI 4 nhóm trong việc giải thích cơ chế xung đột nước hạ du sông liên tỉnh; bổ sung công cụ định lượng vào lý thuyết phân bổ nước công bằng.
  • Hàm ý phương pháp luận: Cung cấp quy trình tích hợp thành công giữa chuỗi mô hình toán học thủy văn - thủy lực DHI với kỹ thuật ra quyết định đa tiêu chí (MCDA - AHP/Delphi), tạo khung mẫu chuẩn cho các nghiên cứu tiếp theo.
  • Hàm ý thực tiễn và chính sách: Luận án đề xuất gói giải pháp tích hợp mang tính đột phá:
    1. Công trình: Xây dựng đập dâng ngăn mặn trên sông Cầu Đỏ / hạ lưu sông Hàn; xây dựng hồ chứa nước Sông Bắc trên lưu vực sông Cu Đê (dung tích $W_{\text{hiệu ích}} \approx 35 \times 10^6\text{ m}^3$); đưa Nhà máy nước Hòa Liên ($120.000\text{ m}^3/\text{ngày}$) vào vận hành chính thức để đa dạng hóa nguồn cấp, giảm tải sự phụ thuộc $89%$ vào cửa thu Cầu Đỏ.
    2. Phi công trình: Điều chỉnh Quy trình vận hành liên hồ chứa lưu vực Vu Gia - Thu Bồn (Quy trình 1865/QĐ-TTg) theo hướng ràng buộc lưu lượng xả kiệt liên tục về đập dâng An Trạch tối thiểu duy trì mực nước $>+2,0\text{ m}$; nâng cao tỷ lệ che phủ rừng đầu nguồn huyện Hòa Vang từ $61,2%$ lên $65,5%$; kéo giảm tỷ lệ thất thoát nước sạch đô thị từ $18,5%$ (năm 2020) xuống dưới $12%$ vào năm 2030. Kết quả mô phỏng chứng minh: khi áp dụng đồng bộ các giải pháp này, chỉ số $WSI$ năm 2030 giảm mạnh từ $3,62$ xuống còn $2,45$ (mức an toàn).

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án thẳng thắn chỉ ra các giới hạn nghiên cứu:

  • Giới hạn về phạm vi tài nguyên: Nghiên cứu chủ động loại trừ tài nguyên nước dưới đất do các quy định pháp lý hành chính đặc thù của địa phương, do đó mô hình chưa phản ánh được động thái tương tác thủy lực giữa nước mặt và các tầng chứa nước ven biển ($qh, qp$).
  • Giới hạn dữ liệu xuyên biên giới và quy hoạch thủy điện: Chưa lượng hóa chi tiết tác động của việc chuyển dòng từ sông Đak Mi sang sông Thu Bồn của thủy điện Đak Mi 4 trong các kịch bản thời tiết cực đoan kéo dài hơn 6 tháng.
  • Điều kiện biên về kinh tế - xã hội: Chuỗi số liệu năm 2020 chịu sự nhiễu loạn cục bộ do đại dịch Covid-19 làm sụt giảm GRDP thành phố $7,99%$, dẫn đến nhu cầu nước của ngành dịch vụ du lịch tạm thời thấp hơn quy luật phát triển dài hạn.

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

  1. Mở rộng tích hợp mô hình dòng ngầm 3D (như MODFLOW/FEFLOW) để đánh giá quá trình xâm thực mặn vào tầng ngậm nước ngầm dưới áp lực biến đổi khí hậu và nước biển dâng.
  2. Ứng dụng trí tuệ nhân tạo (AI/Machine Learning) kết hợp mạng thần kinh nhân tạo (ANN/LSTM) để dự báo độ mặn thời gian thực (Real-time salinity forecasting) tại cửa thu Cầu Đỏ trước 48–72 giờ.
  3. Nghiên cứu cơ chế chi trả dịch vụ môi trường rừng (PES) và cơ chế bồi hoàn kinh tế liên tỉnh giữa Quảng Nam và Đà Nẵng trong quản trị lưu vực sông Vu Gia - Thu Bồn.

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

  • Tác động học thuật: Luận án thiết lập chuẩn mực phương pháp luận mới trong đánh giá căng thẳng nguồn nước tại Việt Nam; cung cấp bộ dữ liệu hiệu chuẩn mô hình DHI chuẩn xác cho lưu vực Vu Gia - Thu Bồn, mở ra tiềm năng trích dẫn cao trong các tạp chí chuyên ngành quốc tế (ISI/Scopus) thuộc nhóm Q1/Q2 về Kỹ thuật Tài nguyên nước và Thủy văn Đô thị.
  • Chuyển đổi công nghiệp & Cấp nước đô thị: Cung cấp cơ sở tính toán kinh tế - kỹ thuật trực tiếp cho Công ty Cổ phần Cấp nước Đà Nẵng (DAWACO) trong việc tối ưu hóa lịch bơm phòng mặn An Trạch, giảm chi phí năng lượng bơm trung gian hàng tỷ đồng mỗi năm.
  • Tác động chính sách công: Là căn cứ khoa học trực tiếp phục vụ UBND Thành phố Đà Nẵng và Bộ Tài nguyên & Môi trường trong việc rà soát, điều chỉnh Quy trình vận hành liên hồ chứa lưu vực sông Vu Gia - Thu Bồn, hiện thực hóa mục tiêu an ninh nguồn nước tại Kết luận số 36/KL-TW của Bộ Chính trị.
  • Lợi ích xã hội: Đảm bảo nguồn cấp nước sinh hoạt an toàn, liên tục cho hơn 1,2 triệu người dân đô thị và hàng chục ngàn doanh nghiệp du lịch - công nghiệp trọng điểm, xóa bỏ triệt để rủi ro mất nước sinh hoạt diện rộng trong các tháng cao điểm mùa khô.

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

                              CÁC NHÓM ĐỐI TƯỢNG HƯỞNG LỢI
                                            │
    ┌───────────────────────┬───────────────┴───────────────┬───────────────────────┐
    ▼                       ▼                               ▼                       ▼
┌──────────────────┐    ┌──────────────────┐    ┌──────────────────┐    ┌──────────────────┐
│ Nghiên cứu sinh  │    │ Giảng viên &     │    │ Doanh nghiệp &   │    │ Nhà hoạch định   │
│ & Học viên SĐH   │    │ Nhà khoa học     │    │ Đơn vị vận hành  │    │ Chính sách       │
└────────┬─────────┘    └────────┬─────────┘    └────────┬─────────┘    └────────┬─────────┘
         │                       │                       │                       │
         ├─ Khung WSI 25 chỉ số  ├─ Khung lý thuyết IWRM ├─ Tối ưu lịch bơm      ├─ Điều chỉnh quy trình
         └─ Quy trình tích hợp   └─ Dữ liệu kiểm định    └─ Kế hoạch đầu tư      └─ Phân bổ nguồn nước
  • Nghiên cứu sinh và Học viên cao học: Tiếp cận khung phương pháp luận kết hợp chuẩn mực giữa AHP, Delphi-KAMET và chuỗi mô hình thủy lực phức hợp; khai thác các hướng nghiên cứu mở rộng về an ninh nguồn nước đô thị.
  • Giảng viên và Nhà khoa học: Khai thác cơ sở dữ liệu thủy văn, ranh mặn và mô hình toán đã được kiểm định chính xác cao ($NSE > 0,85$) cho lưu vực Vu Gia - Thu Bồn để phát triển các đề tài khoa học cấp quốc gia và hợp tác quốc tế.
  • Doanh nghiệp cấp nước và Khu công nghiệp: DAWACO và các ban quản lý KCN (Hòa Khánh, Liên Chiểu) sở hữu công cụ lượng hóa rủi ro nguồn nước, chủ động lập phương án dự trữ và điều tiết cấp nước an toàn.
  • Nhà hoạch định chính sách cấp tỉnh và Bộ, Ngành: Sở TN&MT, Sở Xây dựng Đà Nẵng, Cục Quản lý Tài nguyên nước sử dụng kết quả để xây dựng Quy hoạch tài nguyên nước thời kỳ 2021–2030, tầm nhìn 2050 và kế hoạch thích ứng biến đổi khí hậu.

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à nó mở rộng lý thuyết nào? Trả lời: Luận án đã 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 xây dựng cấu trúc giải thích cơ chế phát sinh căng thẳng nguồn nước tại đô thị cửa sông có nguồn nước liên tỉnh. Khác với các mô hình truyền thống xem WSI là tỷ số cân bằng thể tích nước ($Q_{\text{cung}} / Q_{\text{cầu}}$), luận án chứng minh trạng thái căng thẳng tại Đà Nẵng là hàm số phụ thuộc phi tuyến vào 4 chiều kích: biến động dòng chảy kiệt ($C_v$), xâm nhập mặn động học cửa sông, yêu cầu dòng chảy môi trường (Tennant) và năng lực điều tiết trữ nước hạ du.

2. Đổi mới phương pháp luận của luận án thể hiện thế nào khi so sánh với các nghiên cứu quốc tế trước đây? Trả lời: So với nghiên cứu của Firdaus Ali (2012) tại Jakarta (chỉ có 8 chỉ số, không gán trọng số và bỏ qua xâm nhập mặn) và nghiên cứu của Al-Nozaily và cs. (2022) tại thành phố Ibb (24 chỉ số với trọng số cào bằng), luận án của Nguyễn Đại Trung có 2 bước đột phá phương pháp:

  • Thiết lập cấu trúc thứ bậc 25 chỉ số được chọn lọc qua thuật toán Delphi-KAMET và gán trọng số khách quan bằng ma trận so sánh cặp AHP ($CR < 0,1$).
  • Tích hợp dữ liệu mô phỏng liên hoàn từ chuỗi 3 mô hình toán học DHI (MIKE NAM $\to$ MIKE HYDRO BASIN $\to$ MIKE 11 HD/AD) để cấp đầu vào số liệu thủy lực - mặn trực tiếp cho bộ chỉ số WSI, loại bỏ tính chủ quan của các phương pháp chuyên gia truyền thống.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì? Trả lời: Phát hiện bất ngờ nhất là sự suy giảm lưu lượng mùa kiệt tự nhiên không phải nguyên nhân duy nhất gây thiếu nước ngọt, mà yếu tố quyết định là sự sụt giảm mực nước thủy lực tại đập dâng An Trạch (xuống $+1,44\text{ m}$ ngày 29/3/2021 so với mực nước thiết kế $+2,00\text{ m}$). Sự tụt mực nước này khiến trạm bơm phòng mặn không thể vận hành hết công suất, trong khi độ mặn tại cửa thu Cầu Đỏ vọt lên $4.846\text{ mg/l}$ (gấp 20 lần quy chuẩn). Điều này chỉ ra rằng: giải quyết bài toán căng thẳng nước Đà Nẵng cốt lõi là giải bài toán duy trì cột nước thủy lực phòng mặn tại An Trạch và đập ngăn mặn hạ lưu, chứ không đơn thuần là tăng dung tích khai thác.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) rõ ràng không? Trả lời: Có. Toàn bộ thông số thiết lập mô hình (bước thời gian $\Delta t$, hệ số nhám Manning $n$, thông số dòng chảy mặt/ngầm trong MIKE NAM, ma trận trọng số AHP, bảng phân ngưỡng 5 cấp độ của 25 chỉ số và biểu mẫu Delphi-KAMET) đều được văn bản hóa chi tiết trong Chương 2 và Chương 3, cho phép các nhà nghiên cứu khác áp dụng chính xác quy trình này cho các đô thị ven biển khác như Quy Nhơn, Tuy Hòa, Phan Thiết hay Cần Thơ.

5. Chương trình nghiên cứu 10 năm (10-year Research Agenda) được định hình ra sao? Trả lời: Lộ trình nghiên cứu giai đoạn 2024–2034 tập trung vào 3 trụ cột:

  • Giai đoạn 2024–2026: Xây dựng hệ thống bản sao số (Digital Twin) mạng lưới sông Vu Gia - Thu Bồn và tự động hóa cảnh báo sớm xâm nhập mặn dựa trên AI.
  • Giai đoạn 2027–2030: Đánh giá tác động tương hỗ giữa nước mặt và động lực học nước ngầm ven biển dưới kịch bản nước biển dâng cực đoan.
  • Giai đoạn 2031–2034: Thiết lập mô hình kinh tế tài nguyên nước lượng hóa chi phí - lợi ích trong liên kết vùng Quảng Nam - Đà Nẵng về chia sẻ rủi ro khan hiếm nước.

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Đại Trung đã hoàn thành xuất sắc các mục tiêu nghiên cứu thông qua các đóng góp cụ thể:

  1. Chuẩn hóa khung đánh giá WSI: Phát triển thành công bộ chỉ số gồm 4 nhóm và 25 chỉ số thành phần đặc thù cho đô thị ven biển hạ lưu sông liên tỉnh.
  2. Lượng hóa chính xác thực trạng và xu thế: Xác định chỉ số $WSI$ hiện trạng năm 2020 đạt $2,84$ (mức Trung bình) và dự báo năm 2030 đạt $3,62$ (mức Cao) dưới tác động của BĐKH RCP4.5 và đô thị hóa.
  3. Đột phá phương pháp luận mô hình: Tích hợp thành công chuỗi mô hình thủy văn - cân bằng nước - thủy lực lan truyền mặn MIKE (NAM, HYDRO BASIN, 11 HD/AD) với phương pháp chuyên gia Delphi-KAMET và phân tích thứ bậc AHP ($NSE > 0,85$).
  4. Đề xuất gói giải pháp công trình - phi công trình tối ưu: Chứng minh khả năng kéo giảm chỉ số $WSI$ năm 2030 từ $3,62$ xuống $2,45$ thông qua việc xây dựng đập ngăn mặn Cầu Đỏ, hồ chứa Sông Bắc, vận hành NMN Hòa Liên và tối ưu quy trình liên hồ chứa.
  5. Mở ra 3 luồng nghiên cứu mới: Quản trị an ninh nước hạ lưu sông liên tỉnh; Tích hợp mô hình lan truyền mặn thời gian thực vào quy hoạch đô thị; Kinh tế học chia sẻ tài nguyên nước lưu vực.
  6. Hiện thực hóa chính sách: Đóng góp cơ sở khoa học then chốt phục vụ triển khai Luật Tài nguyên Nước 2023 và Kết luận số 36/KL-TW của Bộ Chính trị về bảo đảm an ninh nguồn nước quốc gia.