Tổng quan về luận án

An ninh nguồn nước (ANNN - Water Security) trong thế kỷ 21 đã trở thành một thách thức mang tính sống còn đối với sự phát triển bền vững toàn cầu, đặc biệt tại các quốc gia đang phát triển sở hữu hệ thống lưu vực sông xuyên biên giới chịu tác động kép từ biến đổi khí hậu (BĐKH) và áp lực nhân sinh. Luận án tiến sĩ kỹ thuật "Nghiên cứu an ninh nguồn nước cho phát triển bền vững lưu vực sông Mã" của tác giả Nguyễn Thị Mùi (chuyên ngành Môi trường đất và nước, Mã số: 9.03, Đại học Thủy lợi, năm 2019, người hướng dẫn: GS.TS Lê Đình Thành) là công trình nghiên cứu tiên phong tại Việt Nam thiết lập một khung đánh giá định lượng toàn diện về ANNN cấp lưu vự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 nghiên cứu trước đây tại Việt Nam phần lớn mới tiếp cận ANNN ở góc độ khái niệm chung, phân tích khung chính sách đơn lẻ (Phạm Thị Tính, 2015), đánh giá thử nghiệm quy mô đô thị hành chính biệt lập (Nguyen Duc Hai, 2015; Babel và cộng sự, 2017), hoặc chỉ khảo sát nguy cơ cho từng đối tượng công trình như thủy điện bậc thang (Nguyễn Trúc Lê, Nguyễn Mạnh Cường, 2016). Trên lưu vực sông Mã – hệ thống sông lớn thứ 3 Việt Nam với tổng diện tích $28.400\text{ km}^2$ ($17.600\text{ km}^2$ thuộc Việt Nam, chiếm 62%; $10.800\text{ km}^2$ thuộc Lào, chiếm 38%) – các nghiên cứu quy hoạch thủy lợi trước đây (Viện Quy hoạch Thủy lợi, 2011, 2012, 2015; Lương Ngọc Chung, 2014; Nguyễn Thanh Hùng, 2016) chỉ tập trung giải quyết các bài toán thủy văn, cân bằng nước cục bộ, biến động lòng dẫn hoặc phòng chống thiên tai đơn ngành mà hoàn toàn vắng bóng một bộ chỉ số tích hợp đa tiêu chí phản ánh bức tranh ANNN tổng thể.

Luận án tập trung giải quyết 3 câu hỏi nghiên cứu (Research Questions - RQ) và kiểm chứng 3 giả thuyết khoa học (Hypotheses - H):

  • RQ1: Cấu trúc thành phần và các tiêu chí đo lường nào phản ánh chính xác nhất bản chất ANNN trên lưu vực sông nhiệt đới gió mùa tại Việt Nam?
    • H1: ANNN lưu vực sông tại Việt Nam là một hàm đa biến tương tác động giữa 6 nhóm trụ cột: Khả năng sẵn có của nguồn nước, Đảm bảo nước cho dân sinh, Đảm bảo nước cho kinh tế, Môi trường sinh thái, Giảm thiểu rủi ro thiên tai, và Năng lực quản trị - quản lý.
  • RQ2: Mức độ đảm bảo ANNN phân hóa như thế nào theo không gian giữa các tiểu vùng và theo thời gian giữa các kịch bản hiện trạng và tương lai?
    • H2: Mặc dù tổng lượng dòng chảy dồi dào, sự phân bố dị thường theo mùa và áp lực phát triển kinh tế sẽ gây suy giảm nghiêm trọng chỉ số ANNN tổng hợp tại các vùng hạ du ven biển.
  • RQ3: Các giải pháp công trình và phi công trình nào có tính khả thi cao nhất nhằm tối ưu hóa bộ chỉ số ANNN đến năm 2030?
    • H3: Việc điều tiết liên hồ chứa kết hợp kiểm soát ô nhiễm nguồn thải và quy hoạch phân bổ hạn ngạch nước có thể cải thiện đáng kể bậc phân cấp ANNN của toàn lưu vực.

Khung lý thuyết của luận án được xây dựng trên sự giao thoa giữa Lý thuyết hệ thống (Systems Theory), Lý thuyết Quản lý tổng hợp tài nguyên nước (IWRM - Integrated Water Resources Management) và Khung phát triển bền vững Rio 1992 tích hợp 3 trụ cột: Kinh tế, Xã hội và Môi trường.

Phạm vi nghiên cứu bao quát phần lưu vực sông Mã thuộc lãnh thổ Việt Nam trên địa bàn 5 tỉnh: Điện Biên, Sơn La, Hòa Bình, Nghệ An và Thanh Hóa; với chuỗi dữ liệu thủy văn thực đo lịch sử 1961–2011, dữ liệu chất lượng nước quan trắc chi tiết năm 2015 (6 đợt khảo sát) và dự báo kịch bản phát triển kinh tế xã hội, BĐKH đến năm 2030 (KBPT 2030 + BĐKH). Đóng góp mang tính đột phá của luận án là đã lượng hóa mức độ ANNN cho lưu vực có tổng dòng chảy mặt lên tới $18\text{ tỷ m}^3/\text{năm}$, chuyển dịch từ quản lý cảm tính sang công cụ định lượng khoa học có khả năng nhân rộng toàn quốc.

Literature Review và Positioning

Khái niệm ANNN đã trải qua quá trình tiến hóa học thuật sâu sắc trên trường quốc tế. Giai đoạn sơ khởi, các nghiên cứu tập trung vào khía cạnh vật lý của sự khan hiếm nước thông qua Chỉ số Căng thẳng Nguồn nước của Falkenmark và cộng sự (1989), đo lường ngưỡng nước mặt bình quân đầu người ($>1.700\text{ m}^3/\text{người/năm}$ là an toàn, $1.000\text{--}1.700\text{ m}^3$ là căng thẳng, $<500\text{ m}^3$ là cực kỳ khan hiếm). Tiếp đó, Raskin và cộng sự (1997) định lượng chỉ số căng thẳng bằng tỷ lệ phần trăm lượng nước khai thác so với lượng nước sẵn có ($>40%$ là căng thẳng nghiêm trọng).

Tuy nhiên, các mô hình đơn biến trên bộc lộ hạn chế lớn khi bỏ qua yếu tố chất lượng nước, năng lực hạ tầng kỹ thuật và khả năng thích ứng của xã hội. Bước chuyển giao lý thuyết diễn ra khi Tổ chức Hợp tác vì Nước Toàn cầu (GWP, 2000) đưa ra định nghĩa mang tính nhân văn: "ANNN là sự đảm bảo an toàn nguồn nước ở mọi cấp độ từ gia đình đến toàn cầu, có nghĩa là mọi người đều được cung cấp đủ nước sạch với chi phí phải chăng phục vụ các nhu cầu thiết yếu trong cuộc sống, đồng thời đảm bảo môi trường tự nhiên được bảo tồn và phát huy". Sadoff (2007) mở rộng định nghĩa gắn liền với quản trị rủi ro: "ANNN là lượng nước sẵn có, đảm bảo về trữ lượng và chất lượng cho sức khỏe, sinh hoạt, hệ sinh thái và sản xuất đã tính đến khả năng xảy ra rủi ro cho con người, môi trường và nền kinh tế".

Trên bình diện phương pháp luận, tồn tại hai luồng quan điểm đối lập:

  1. Trường phái Không gian Địa vật lý thuần túy (tiêu biểu là Vörösmarty và cộng sự, 2010): Sử dụng khung không gian định lượng 23 chỉ số áp lực môi trường đe dọa ANNN của con người và đa dạng sinh học thủy vực toàn cầu, chú trọng các yếu tố nhiễu loạn lưu vực, ô nhiễm, dòng chảy gián đoạn.
  2. Trường phái Thể chế - Kinh tế Xã hội và Rủi ro đa chiều (tiêu biểu là OECD, 2013; AWDO, 2013; Dunn và cộng sự, 2014): Nhấn mạnh ANNN là việc chấp nhận và kiểm soát 4 rủi ro lớn (thiếu hụt nước, chất lượng không đạt, lũ lụt dư thừa, và suy giảm phục hồi hệ sinh thái), đề cao sự tham gia của các bên liên quan và phân bổ kinh tế.
TIẾN TRÌNH TIẾP CẬN LÝ THUYẾT AN NINH NGUỒN NƯỚC (ANNN)
┌───────────────────────────────┐
│ Falkenmark (1989)             │ ──> Căng thẳng vật lý đơn biến (m³/người/năm)
│ Raskin (1997)                 │ ──> Tỷ lệ khai thác/trữ lượng nước sẵn có (%)
└──────────────┬────────────────┘
               │
┌──────────────▼────────────────┐
│ GWP (2000), Sadoff (2007),    │ ──> Tiếp cận tích hợp IWRM, sinh kế,
│ OECD (2013)                   │     kinh tế - xã hội & quản trị rủi ro
└──────────────┬────────────────┘
               │
┌──────────────▼────────────────┐
│ Vörösmarty (2010), AWDO (2013)│ ──> Mô hình hóa đa tiêu chí cấp lưu vực & quốc gia
│ Babel (2017), Wu (2018)       │     (5 nhóm chiều không gian: Dân sinh, Kinh tế,
└──────────────┬────────────────┘     Môi trường, Thiên tai, Quản trị)
               │
┌──────────────▼────────────────┐
│ LUẬN ÁN NGUYỄN THỊ MÙI (2019) │ ──> Bản địa hóa khung WSI cho lưu vực sông Việt Nam
│ Sông Mã (28.400 km²)          │     Tích hợp MIKE-NAM + WEAP + WQI + Tennant + AHP
└───────────────────────────────┘

Trong bức tranh tổng quan đó, nghiên cứu của Nguyễn Thị Mùi được định vị tại điểm giao thoa giữa các khung quốc tế tiên tiến và đặc thù lưu vực sông nhiệt đới chịu ảnh hưởng xuyên biên giới. So sánh với hai nghiên cứu quốc tế điển hình:

  • Nghiên cứu của S. Xiao và cộng sự (2008) tại lưu vực sông Heihe/Shiyanghe (Tây Bắc Trung Quốc) xây dựng 21 chỉ số chia làm 4 nhóm cấp 1 và 8 nhóm cấp 2 sử dụng phương pháp AHP để giải quyết mâu thuẫn nước vùng khô hạn. Luận án kế thừa cấu trúc phân cấp đa tiêu chí của Xiao nhưng bổ sung chiều kích phân tích xâm nhập mặn cửa sông và điều tiết bậc thang thủy điện vùng cận nhiệt đới.
  • Nghiên cứu của Babel và cộng sự (2017) đánh giá ANNN cho lưu vực sông Banas (Ấn Độ), Chao Phraya (Thái Lan) và sông Hồng (Việt Nam) đã xây dựng 5 nhóm với 8 chỉ số cho cấp lưu vực nhưng mới dừng lại ở việc tính điểm riêng lẻ cho từng nhóm mà chưa lượng hóa thành một chỉ số tổng hợp WSI duy nhất. Luận án đã vượt lên giới hạn này bằng cách thiết lập thuật toán xác định trọng số và tính toán chỉ số ANNN tổng hợp cuối cùng cho từng tiểu vùng không gian.

Đó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 và kiểm chứng lý thuyết Quản lý tổng hợp tài nguyên nước (IWRM) và Lý thuyết Hệ thống Tài nguyên Nước - Xã hội (Socio-Hydrological Systems) trong bối cảnh các lưu vực sông có cấu trúc dòng chảy phụ thuộc nguồn nước ngoại biên. Nghiên cứu đã chứng minh rằng: "Khái niệm an ninh nguồn nước không thể đồng nhất với sự dồi dào về trữ lượng thủy văn tổng thể, mà là trạng thái cân bằng động giữa khả năng cung ứng được kiểm soát, chất lượng nước theo mục đích sử dụng, sức tải sinh thái và năng lực chống chịu thiên tai".

Mô hình lý thuyết của luận án xác lập 3 mệnh đề (Propositions):

  • Mệnh đề 1 (P1): Tính bất đối xứng không - thời gian của dòng chảy tự nhiên là yếu tố quyết định mức độ dễ bị tổn thương về ANNN hơn là tổng lượng dòng chảy trung bình năm ($W_0$).
  • Mệnh đề 2 (P2): Sự can thiệp của hệ thống công trình tích nước thượng nguồn làm thay đổi căn bản động lực học dòng chảy kiệt và cân bằng mặn hạ du, chuyển đổi rủi ro từ dạng tự nhiên sang rủi ro điều tiết kỹ thuật.
  • Mệnh đề 3 (P3): Chỉ số ANNN tổng hợp là một hàm phi tuyến tính của các chỉ số thành phần, trong đó suy thoái chất lượng nước cục bộ (WQI) và rủi ro mặn có thể triệt tiêu hoàn toàn lợi thế về trữ lượng nước dồi dào.

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

Khung phân tích của luận án tích hợp đồng thời 3 nền tảng lý thuyết: Lý thuyết Thủy văn động lực, Kinh tế học tài nguyên nước và Lý thuyết Đánh giá đa tiêu chí (MCDA). Cấu trúc bộ chỉ số ANNN lưu vực sông Việt Nam được thiết lập gồm 6 nhóm chỉ số trụ cột ($WSI_1$ đến $WSI_6$) bao gồm 17 chỉ số thành phần:

  1. Nhóm 1 - Nguồn nước ($WSI_1$): Đo lường tiềm năng tài nguyên nước mặt và nước dưới đất bình quân đầu người, mức độ phụ thuộc nguồn nước bên ngoài ($38%$ diện tích lưu vực nằm trên lãnh thổ Lào).
  2. Nhóm 2 - Nước cho đời sống dân sinh ($WSI_2$): Đo lường tỷ lệ dân số đô thị và nông thôn được cấp nước sạch hợp vệ sinh, mức độ tiêu thụ nước sinh hoạt bình quân ($L/\text{người/ngày}$).
  3. Nhóm 3 - Nước cho các ngành kinh tế ($WSI_3$): Đo lường mức độ đáp ứng nhu cầu nước cho sản xuất nông nghiệp (tưới lúa $303.238\text{ ha}$), sản xuất công nghiệp (các KCN Bỉm Sơn, Lễ Môn, Đình Hương, Lam Sơn, KKT Nghi Sơn) và thủy sản.
  4. Nhóm 4 - Môi trường và chất lượng nước ($WSI_4$): Đo lường tỷ lệ dòng chảy môi trường tối thiểu ($DCMT$ tính theo phương pháp Tennant), chỉ số chất lượng nước mặt ($WQI$), và tỷ lệ thu gom/xử lý nước thải sinh hoạt - công nghiệp.
  5. Nhóm 5 - Rủi ro thiên tai liên quan đến nước ($WSI_5$): Đánh giá mức độ tổn thương do lũ lụt, diện tích hạn hán, và mức độ xâm nhập mặn sâu vào nội đồng.
  6. Nhóm 6 - Quản trị và thể chế tài nguyên nước ($WSI_6$): Đo lường mức độ hoàn thiện quy hoạch, hiệu lực quy trình vận hành liên hồ chứa và hệ thống quan trắc giám sát tự động.

Toàn bộ các chỉ số được chuẩn hóa trên thang đo 5 mức độ đảm bảo: Mức 1 (Rất kém/Kém an ninh nhất), Mức 2 (Kém), Mức 3 (Trung bình), Mức 4 (Tốt), và Mức 5 (Rất tốt/Đảm bảo an ninh cao). Trọng số ($\omega_i$) của các nhóm và chỉ số được xác định thông qua ma trận phân tích thứ bậc AHP kết hợp tham vấn chuyên gia theo kỹ thuật Delphi.

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 bản thể luận hiện thực phản biện (Critical Realism) kết hợp phương pháp luận thực chứng (Positivism), kết hợp chặt chẽ giữa mô hình hóa toán học thủy văn định lượng chính xác cao và phương pháp đánh giá thứ cấp kinh tế - xã hội. Thiết kế nghiên cứu đa tầng (Multi-level design) chia không gian lưu vực sông Mã thành các vùng tính toán cân bằng nước đặc trưng: Vùng thượng lưu sông Mã, Trung lưu sông Mã, Lưu vực sông Chu, Sông Bưởi, Sông Cầu Chày, Sông Hoạt và Vùng đồng bằng hạ du (phân lưu sông Lèn và sông Lạch Trường).

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu tuân thủ tiêu chuẩn học thuật khắt khe thông qua 4 giai đoạn logic:

  1. Thu thập và thẩm định dữ liệu: Sử dụng chuỗi số liệu khí tượng - thủy văn đồng bộ tại 12 trạm thủy văn chính (Xã Là, Cẩm Thủy, Cửa Đạt, Xuân Khánh, Vụ Bản...) và 11 trạm đo mưa đặc trưng (Lạc Sơn, Tuần Giáo, Mường Lát, Hồi Xuân, Bái Thượng, Thanh Hóa, Tĩnh Gia...). Số liệu chất lượng nước được kiểm chuẩn từ 6 đợt quan trắc môi trường năm 2015 tại 9 vị trí trên dòng chính sông Mã, 7 vị trí trên sông Chu, 3 vị trí trên sông Cầu Chày và sông Bưởi.
  2. Mô hình hóa dòng chảy và thủy văn: Sử dụng mô hình mưa - dòng chảy MIKE-NAM để tính toán tái tạo và dự báo chuỗi dòng chảy đến tại các nút tính toán trên toàn mạng lưới sông ứng với tần suất nước kiệt thiết kế $P = 85%$.
  3. Mô hình hóa cân bằng nước toàn diện: Ứng dụng mô hình WEAP (Water Evaluation and Planning System) do Viện Môi trường Stockholm (SEI) phát triển. Cấu trúc mô hình WEAP lưu vực sông Mã được thiết lập với hệ thống các nút cung cấp (nguồn nước đến từ sông chính, sông nhánh, các hồ chứa lớn: Cửa Đạt $1.450\times 10^6\text{ m}^3$, Hủa Na $569,35\times 10^6\text{ m}^3$, Trung Sơn $348,5\times 10^6\text{ m}^3$, Sông Mực $200\times 10^6\text{ m}^3$, Yên Mỹ $84,4\times 10^6\text{ m}^3$) và các nút nhu cầu sử dụng nước (nông nghiệp, công nghiệp, dân sinh, dòng chảy sinh thái EWR).
  4. Mô hình hóa xâm nhập mặn: Ứng dụng phần mềm thủy lực - chất lượng nước 1 chiều MIKE 11 để mô phỏng lan truyền mặn trên hệ thống sông Mã, sông Yên, sông Lèn và sông Lạch Trường.
QUY TRÌNH MÔ HÌNH HÓA VÀ TÍNH TOÁN CHỈ SỐ ANNN
┌──────────────────────────────────────────────────────────┐
│              DỮ LIỆU ĐẦU VÀO ĐA NGUỒN                    │
│  KTTV (1961-2011) │ Chất lượng nước (WQI) │ KTXH, GIS    │
└──────────────┬───────────────────────────────┬───────────┘
               │                               │
┌──────────────▼─────────────┐   ┌─────────────▼───────────┐
│     MÔ HÌNH MIKE-NAM       │   │     MÔ HÌNH MIKE 11     │
│  Mưa - dòng chảy kiệt P=85%│   │ Thủy lực & Xâm nhập mặn │
└──────────────┬─────────────┘   └─────────────┬───────────┘
               │                               │
               └───────────────┬───────────────┘
                               │
                ┌──────────────▼──────────────┐
                │        MÔ HÌNH WEAP         │
                │  Cân bằng nước hệ thống     │
                │  (Nhu cầu vs Cung cấp)      │
                │  Dòng chảy MT (Tennant)     │
                └──────────────┬──────────────┘
                               │
                ┌──────────────▼──────────────┐
                │   CHUẨN HÓA BỘ CHỈ SỐ WSI   │
                │   6 Nhóm - 17 Chỉ số WSI    │
                │   Trọng số AHP & Delphi     │
                └──────────────┬──────────────┘
                               │
                ┌──────────────▼──────────────┐
                │   CHỈ SỐ ANNN TỔNG HỢP      │
                │  Đánh giá phân vùng 2015    │
                │  Dự báo kịch bản 2030+BĐKH  │
                └─────────────────────────────┘

Data và phân tích

Tính toán dòng chảy môi trường ($DCMT$) được luận án áp dụng phương pháp thực nghiệm nổi tiếng của Tennant (1976) – xác định tỷ lệ phần trăm theo lưu lượng dòng chảy năm trung bình nhiều năm ($Q_0$). Cụ thể, trong mùa kiệt quy định ngưỡng $10%\text{--}20% Q_0$ để duy trì sinh thái thủy sinh tối thiểu và $30%\text{--}50% Q_0$ trong mùa lũ.

Chỉ số chất lượng nước tổng hợp ($WQI$) được tính toán chặt chẽ theo Quyết định số 46/QĐ-TCMT (2011) của Tổng cục Môi trường với 5 mức phân loại màu sắc quy chuẩn: Xanh nước biển ($91\text{--}100$: cấp nước sinh hoạt tốt), Xanh lá cây ($76\text{--}90$: cấp nước sinh hoạt qua xử lý), Vàng ($51\text{--}75$: tưới tiêu thủy lợi), Da cam ($26\text{--}50$: giao thông thủy), và Đỏ ($0\text{--}25$: ô nhiễm nặng).

Các kịch bản cân bằng nước được xây dựng chi tiết:

  • Kịch bản Hiện trạng (KB HT 2015).
  • Kịch bản Phát triển 2030 (KBPT 2030): Tăng trưởng công nghiệp $24,3%/\text{năm}$, mở rộng diện tích đô thị và thủy sản.
  • Kịch bản Phát triển 2030 có tính đến Biến đổi khí hậu (KBPT 2030 + BĐKH): Giả định nền nhiệt độ tăng, lượng bốc hơi tăng, lượng mưa mùa khô suy giảm từ $10%\text{--}15%$, và nước biển dâng làm gia tăng ranh giới mặn $1\text{‰}$ và $4\text{‰}$.

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

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

  1. Bản chất mâu thuẫn giữa tổng lượng nước dồi dào và sự thiếu hụt nghiêm trọng trong mùa kiệt: Tổng lượng dòng chảy mặt bình quân nhiều năm của lưu vực sông Mã đạt $18\text{ tỷ m}^3$ (tương ứng lưu lượng trung bình $Q_0 = 570\text{ m}^3/\text{s}$, mô đun dòng chảy $M_0 = 20,1\text{ L/s}\cdot\text{km}^2$). Tuy nhiên, dòng chảy phân bố cực kỳ bất đối xứng: mùa lũ (tháng VI đến tháng XI) chiếm $70%\text{--}90%$ tổng lượng nước, trong khi 3 tháng kiệt nhất (tháng II, III, IV) chỉ chiếm vỏn vẹn $8%\text{--}8,5%$ tổng lượng dòng chảy năm (riêng tháng III chỉ đạt $2,5%\text{--}2,7%$, mô đun dòng chảy kiệt tụt xuống $M_0 = 6,0\text{--}6,5\text{ L/s}\cdot\text{km}^2$). Điển hình trong năm hạn lịch sử 2010, lưu lượng kiệt tại trạm Sét Thôn (sông Mã) chỉ còn $60\text{ m}^3/\text{s}$, và tại trạm Phong Mục (sông Lèn) chỉ còn $3\text{ m}^3/\text{s}$, gây thiếu nước tưới trên diện rộng.

  2. Sự tổn thương sâu sắc của các phân lưu hạ du trước nguy cơ xâm nhập mặn và ô nhiễm hóa sinh: Kết quả tính toán $WQI$ năm 2015 bộc lộ sự suy giảm chất lượng nước nghiêm trọng cục bộ:

    • Trên dòng chính sông Mã, vào đợt 4 (mùa mưa lũ cuốn theo hàm lượng chất rắn lơ lửng và chất ô nhiễm tràn bề mặt cực cao), chỉ số $WQI$ tại Cầu Na Sài tụt xuống $21,13$, Phà La Hán đạt $21,11$, Cửa Hà đạt $21,40$ (đều ở mức Màu Đỏ - Ô nhiễm nặng).
    • Vùng hạ lưu cảng Lễ Môn và Cửa Hới thường xuyên bị suy giảm chất lượng nước do tiếp nhận nước thải công nghiệp từ các KCN Đình Hương - Tây Ga, Lễ Môn và đô thị TP. Thanh Hóa ($WQI$ đợt 6 tại Cảng Lễ Môn chỉ đạt $55,74$; $WQI$ đợt 4 tại Cửa Hới chỉ đạt $38,30$).
    • Về độ mặn: Khi lưu lượng xả thượng nguồn suy giảm, ranh giới mặn xâm nhập sâu vào nội đồng. Điển hình năm 2010 mặn xâm nhập qua trạm Giàng với độ mặn cực đại đo được lên tới $6,1\text{‰}$. Các sông Lèn và sông Lạch Trường là các điểm nóng nhạy cảm nhất về ANNN do vừa tiếp nhận nguồn nước thải vừa chịu áp lực triều mặn từ biển.
BIẾN THIÊN CHỈ SỐ CHẤT LƯỢNG NƯỚC (WQI) TẠI CÁC ĐIỂM NÓNG DÒNG CHÍNH SÔNG MÃ (2015)
WQI
100 ┌──────────────────────────────────────────────────────────┐
 90 │ █  █  █        █  █     █  █  █        █        █  █  █  │
 80 │ █  █  █     █  █  █     █  █  █     █  █  █     █  █  █  │
 70 │ █  █  █     █  █  █     █  █  █     █  █  █     █  █  █  │
 60 │ █  █  █  █  █  █  █     █  █  █     █  █  █  █  █  █  █  │
 50 │ █  █  █  █  █  █  █     █  █  █     █  █  █  █  █  █  █  │
 40 │ █  █  █  █  █  █  █     █  █  █     █  █  █  █  █  █  █  │
 30 │ █  █  █  █  █  █  █     █  █  █     █  █  █  █  █  █  █  │
 20 │ █  █  █ [█] █  █  █     █  █  █ [█] █  █  █     █  █  █  │ (Đợt 4: Mưa lũ cuốn trôi,
 10 │ █  █  █  █  █  █  █     █  █  █  █  █  █  █     █  █  █  │  WQI tụt đáy mức Ô nhiễm nặng)
  0 └─┴──┴──┴──┴──┴──┴──┴─────┴──┴──┴──┴──┴──┴──┴─────┴──┴──┴──┴─┘
      Đ1 Đ2 Đ3 Đ4 Đ5 Đ6       Đ1 Đ2 Đ3 Đ4 Đ5 Đ6       Đ1 Đ2 Đ3 Đ4 Đ5 Đ6
       [Cầu Na Sài]             [Phà La Hán]             [Cầu Bản Lát]

  1. Chênh lệch cấu trúc an ninh nước dân sinh giữa thành thị và nông thôn: Tổng dân số lưu vực là $4.074.686\text{ người}$, trong đó dân cư nông thôn chiếm đại đa số với $3.500.374\text{ người}$ ($86%$), thành thị chỉ chiếm $574.312\text{ người}$ ($14%$). Mật độ phân bố phân hóa cực hạn từ $2.407\text{ người/km}^2$ tại TP. Thanh Hóa xuống $31\text{ người/km}^2$ tại huyện vùng cao Sốp Cộp. Tỷ lệ tiếp cận nước sạch đạt chuẩn QCVN và công trình cấp nước tập trung ở vùng nông thôn miền núi còn ở mức thấp, khiến nhóm chỉ số $WSI_2$ tại khu vực miền núi chỉ đạt mức trung bình yếu (Mức 2 - Mức 3).

  2. Dự báo lượng nước thiếu hụt gia tăng mạnh trong Kịch bản 2030 + BĐKH: Mô hình WEAP chỉ ra rằng đến năm 2030, sự gia tăng nhu cầu nước của các ngành kinh tế (công nghiệp, khu kinh tế Nghi Sơn, tưới nông nghiệp công nghệ cao) kết hợp với BĐKH sẽ khiến lượng nước thiếu cục bộ trong mùa kiệt tại các tiểu vùng sông Cầu Chày, hạ lưu sông Lèn và sông Mã gia tăng từ $20%\text{--}35%$ so với năm 2015 nếu không có sự can thiệp điều tiết liên hồ chứa hiệu quả.

Implications đa chiều

  • Về mặt Lý thuyết: Luận án cung cấp minh chứng thực nghiệm khẳng định tính tương thích của khung đánh giá ANNN 6 trụ cột đối với các lưu vực sông có tính chất phức hợp tại khu vực Đông Nam Á, đóng góp vào lý thuyết quản lý thủy văn xuyên biên giới.
  • Về mặt Phương pháp luận: Thiết lập quy trình chuẩn kết hợp chuỗi công cụ mô hình MIKE-NAM + WEAP + MIKE 11 + AHP + GIS, tạo tiền đề phương pháp cho các nghiên cứu đánh giá ANNN trên 8 hệ thống sông lớn còn lại của Việt Nam (sông Hồng, sông Mê Kông, sông Đồng Nai, sông Ba...).
  • Về mặt Ứng dụng Thực tiễn và Chính sách:
    • Cung cấp cơ sở khoa học định lượng để UBND tỉnh Thanh Hóa và Bộ Nông nghiệp & PTNT, Bộ Tài nguyên & Môi trường ban hành quy chế phối hợp quản lý lưu vực sông liên tỉnh (Hòa Bình - Thanh Hóa - Sơn La).
    • Tối ưu hóa quy trình vận hành liên hồ chứa trên lưu vực sông Mã (Hồ Cửa Đạt, Trung Sơn, Bá Thước 1 & 2, Hủa Na) nhằm đảm bảo lưu lượng xả duy trì dòng chảy môi trường ($DCMT$) và đẩy mặn cho các trạm bơm cấp nước sinh hoạt vùng hạ du.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận 4 hạn chế khoa học (Boundary Conditions & Limitations):

  1. Thiếu hụt dữ liệu thủy văn xuyên biên giới: Do $38%$ diện tích lưu vực ($10.800\text{ km}^2$) và $3,9\text{ tỷ m}^3$ dòng chảy sản sinh trên lãnh thổ CHDCND Lào, việc thiếu chuỗi số liệu quan trắc khí tượng thủy văn và quy hoạch phát triển kinh tế thực tế từ phía thượng nguồn Lào buộc mô hình phải sử dụng các giả định tính toán dòng chảy đến biên giới, tạo ra độ bất định nhất định.
  2. Mức độ tích hợp dữ liệu nước dưới đất: Luận án tập trung chủ yếu vào tài nguyên nước mặt; mối quan hệ thủy lực tương hỗ giữa nước mặt và các tầng chứa nước lỗ hổng (dày $10\text{--}35\text{ m}$ tại đồng bằng Thanh Hóa) và nước khe nứt vùng karst chưa được mô phỏng chi tiết bằng mô hình dòng ngầm chuyên sâu (như MODFLOW).
  3. Tính động của trọng số chuyên gia: Phương pháp AHP và Delphi mang tính chất chuyên gia tại thời điểm nghiên cứu (2015–2018), chưa tích hợp cơ chế học máy (Machine Learning) để cập nhật trọng số thích ứng tự động theo biến động kinh tế - xã hội thời gian thực.
  4. Giới hạn chuỗi quan trắc chất lượng nước: Dữ liệu chất lượng nước quan trắc định kỳ 6 đợt/năm tại các trạm chưa phản ánh liên tục 24/7 diễn biến ô nhiễm hóa chất công nghiệp và sự cố tràn dầu/hóa chất đột xuất.

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

  • Mở rộng ứng dụng mô hình kết hợp Thủy văn - Sinh thái - Thể chế (Hydro-Ecological-Institutional Model) cho toàn bộ lưu vực xuyên biên giới Việt - Lào.
  • Xây dựng hệ thống hỗ trợ ra quyết định (DSS - Decision Support System) dựa trên WebGIS và dữ liệu viễn thám vệ tinh (Remote Sensing) giám sát ANNN theo thời gian thực.
  • Nghiên cứu cơ chế chi trả dịch vụ môi trường rừng (PES) và thị trường hạn ngạch trao đổi quyền khai thác nước mặt liên tỉnh trên lưu vực sông Mã.

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

Công trình luận án tiến sĩ của NCS. Nguyễn Thị Mùi tạo ra các tác động sâu rộng trên nhiều bình diện:

  • Tác động Học thuật (Academic Impact): Kết quả nghiên cứu đã được công bố trên các tạp chí khoa học chuyên ngành uy tín trong nước (Tạp chí Khoa học & Công nghệ Thủy lợi, Tạp chí Nông nghiệp & PTNT) và kỷ yếu hội thảo quốc tế. Luận án trở thành tài liệu tham khảo chuẩn mực cho các chương trình đào tạo sau đại học ngành Kỹ thuật tài nguyên nước, Quản lý môi trường tại Đại học Thủy lợi, Đại học Quốc gia Hà Nội và Đại học Hồng Đức.
  • Tác động Thể chế và Chính sách (Policy Influence): Luận án là luận cứ khoa học trực tiếp phục vụ việc lập "Quy hoạch tổng thể tài nguyên nước lưu vực sông Mã đến năm 2030, tầm nhìn 2050" và hỗ trợ công tác thẩm định hạn ngạch cấp phép khai thác tài nguyên nước mặt của Sở Tài nguyên và Môi trường tỉnh Thanh Hóa.
  • Lợi ích Kinh tế - Xã hội và Môi trường: Giúp giảm thiểu rủi ro thiệt hại kinh tế do hạn mặn hàng trăm tỷ đồng mỗi năm cho ngành sản xuất lúa và nuôi trồng thủy sản hạ du; bảo vệ nguồn nước an toàn cho hơn 4 triệu dân sinh sống trên toàn lưu vực.

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

CƠ CẤU ĐỐI TƯỢNG HƯỞNG LỢI TỪ CÔNG TRÌNH NGHIÊN CỨU
┌───────────────────────────────┬────────────────────────────────────────┐
│ NHÓM ĐỐI TƯỢNG               │ GIÁ TRỊ VÀ LỢI ÍCH TIẾP NHẬN           │
├───────────────────────────────┼────────────────────────────────────────┤
│ 1. Nhà nghiên cứu & NCS      │ Khung phương pháp luận WSI chuẩn hóa,   │
│    (Doctoral Researchers)     │ bộ công cụ MIKE-NAM + WEAP + Tennant   │
├───────────────────────────────┼────────────────────────────────────────┤
│ 2. Cơ quan hoạch định QLNN   │ Cơ sở phân bổ hạn ngạch tài nguyên     │
│    (Bộ NN&PTNT, Bộ TN&MT)     │ nước, quy hoạch vận hành liên hồ chứa  │
├───────────────────────────────┼────────────────────────────────────────┤
│ 3. Đơn vị QL khai thác CTTL  │ Lịch trình xả nước đẩy mặn tối ưu,     │
│    (Công ty Thủy nông, Thủy điện) │ giải quyết xung đột lợi ích nước      │
├───────────────────────────────┼────────────────────────────────────────┤
│ 4. Cộng đồng dân cư lưu vực  │ Đảm bảo an ninh cấp nước sạch sinh     │
│    (>4 triệu người)           │ hoạt, giảm thiểu hiểm họa hạn - mặn    │
└───────────────────────────────┴────────────────────────────────────────┘

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?

Luận án đã mở rộng Lý thuyết Khung An ninh Nguồn nước đa chiều (Multi-dimensional Water Security Framework) của Sadoff (2007) và AWDO (2013) bằng việc thiết lập mô hình toán tích hợp 6 nhóm chỉ số ($WSI$) thích ứng riêng cho lưu vực sông nhiệt đới xuyên biên giới tại Việt Nam. Đóng góp độc đáo nhất là việc đưa chỉ số Dòng chảy môi trường theo phương pháp Tennant kết hợp với Chỉ số chất lượng nước tổng hợp WQI vào làm biến số ràng buộc cốt lõi trong mô hình cân bằng nước WEAP, bác bỏ quan điểm truyền thống chỉ đo lường an ninh nước bằng tổng thể tích dòng chảy đơn thuần.

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 S. Xiao và cộng sự (2008) tại Trung Quốc (chỉ dùng AHP tĩnh trên các số liệu hành chính) và nghiên cứu của Babel và cộng sự (2017) tại sông Hồng/sông Banas (chưa tổng hợp thành chỉ số đại diện duy nhất), luận án đã thực hiện một bước nhảy vọt về phương pháp:

  • Kết nối liền mạch chuỗi mô hình động lực: Mô hình toán thủy văn (MIKE-NAM) $\rightarrow$ Mô hình cân bằng nước không gian (WEAP) $\rightarrow$ Mô hình lan truyền mặn (MIKE 11) $\rightarrow$ Chuẩn hóa đa tiêu chí (AHP/Delphi).
  • Cho phép mô phỏng dự báo định lượng chỉ số ANNN dưới tác động đồng thời của kịch bản BĐKH toàn cầu và sự phát triển kinh tế xã hội nội tại đến năm 2030.

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

Phát hiện bất ngờ nhất là hiện tượng "Nghịch lý mùa lũ gây mất an ninh chất lượng nước": Tại các trạm trung lưu dòng chính sông Mã (Cầu Na Sài, Phà La Hán, Cửa Hà), chất lượng nước mặt bị suy thoái nghiêm trọng nhất không phải vào mùa cạn kiệt lưu lượng mà rơi vào đợt quan trắc mùa mưa lũ (Đợt 4), chỉ số $WQI$ tụt giảm xuống mức nguy hại ($21,11\text{--}21,40$). Nguyên nhân do dòng chảy mặt mùa lũ cuốn theo khối lượng khổng lồ chất thải rắn lơ lửng, dư lượng thuốc bảo vệ thực vật từ đất dốc canh tác nông nghiệp và chất thải sinh hoạt chưa qua xử lý đổ dồn vào lòng sông.

4. Luận án có cung cấp quy trình nhân rộng (Replication Protocol) không?

Có. Luận án cung cấp tường minh toàn bộ quy trình 5 bước: (1) Thu thập và chuẩn hóa dữ liệu khí tượng - thủy văn - chất lượng nước - KTXH; (2) Thiết lập và hiệu chỉnh thông số mô hình MIKE-NAM và WEAP; (3) Tính toán phân ngưỡng 17 chỉ số theo thang điểm 1-5; (4) Xác định ma trận trọng số AHP thông qua tham vấn chuyên gia; (5) Tính toán chỉ số tổng hợp $WSI$ và hiển thị bản đồ hóa GIS. Quy trình này có thể áp dụng nguyên vẹn cho bất kỳ lưu vực sông nào tại Việt Nam.

5. Định hướng chương trình nghiên cứu 10 năm tiếp theo được phác thảo như thế nào?

Chương trình 10 năm tập trung vào 3 trụ cột: (i) Hợp tác quốc tế nghiên cứu mạng lưới quan trắc khí tượng thủy văn tự động xuyên biên giới trên phần lưu vực sông Mã thuộc Lào; (ii) Tích hợp mô hình số trị dòng ngầm 3D để đánh giá cạn kiệt tầng chứa nước dưới đất vùng đồng bằng ven biển; (iii) Xây dựng nền tảng Trí tuệ nhân tạo (AI) và Bản sao số (Digital Twin) lưu vực sông Mã phục vụ điều hành xả lũ và chống hạn mặn theo thời gian thực.

Kết luận

Công trình nghiên cứu tiến sĩ của tác giả Nguyễn Thị Mùi đã giải quyết trọn vẹn một bài toán khoa học cấp thiết với 6 đóng góp cốt lõi:

  1. Hệ thống hóa và bản địa hóa thành công cơ sở lý luận, phương pháp luận về an ninh nguồn nước phù hợp với điều kiện thể chế và tự nhiên tại Việt Nam.
  2. Xây dựng hoàn chỉnh Bộ chỉ số ANNN lưu vực sông Việt Nam gồm 6 nhóm, 17 chỉ số thành phần với thang đánh giá 5 bậc phân cấp định lượng rõ ràng.
  3. Ứng dụng thành công chuỗi công cụ mô hình hóa tiên tiến (MIKE-NAM, WEAP, MIKE 11, WQI, Tennant) mô phỏng chính xác cân bằng nước và mức độ ANNN lưu vực sông Mã cho năm hiện trạng 2015 và kịch bản 2030 có xét đến BĐKH.
  4. Phát hiện và lượng hóa các điểm nóng dễ bị tổn thương về ANNN tại các vùng hạ lưu phân lưu sông Lèn, sông Lạch Trường và sông Cầu Chày do áp lực kép từ xâm nhập mặn ($6,1\text{‰}$) và ô nhiễm cục bộ.
  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, trong đó trọng tâm là quy chế vận hành tối ưu hóa liên hồ chứa bậc thang kết hợp kiểm soát xả thải công nghiệp.
  6. Mở ra hướng tiếp cận liên ngành mới trong quản trị lưu vực sông, chuyển dịch căn bản từ tư duy "quản lý khai thác nguồn nước đơn thuần" sang mô hình "quản trị an ninh nguồn nước tích hợp vì sự phát triển bền vững".