Tổng quan về luận án
Nghiên cứu của NCS. Nguyễn Đức Phong với đề tài "Nghiên cứu tác động của các loại nguồn xả thải đến chất lượng nước vùng Bán đảo Cà Mau" (Chuyên ngành: Môi trường đất và nước, Mã số: 9 44 03 03, Viện Khoa học Thủy lợi Việt Nam, 2023, người hướng dẫn: GS.TS Tăng Đức Thắng và GS.TS Nguyễn Vũ Việt) đánh dấu bước tiến tiên phong trong việc lượng hóa tương tác đa nguồn ô nhiễm trong hệ thống sông kênh chịu tác động phức tạp của chế độ bán nhật triều Biển Đông và nhật triều Biển Tây. Trong bối cảnh biến đổi khí hậu, sự chuyển dịch cơ cấu nông nghiệp từ ngọt sang lợ - mặn (từ trồng lúa sang nuôi tôm công nghiệp sau giai đoạn 1999–2000) cùng tốc độ đô thị hóa nhanh đã tạo ra xung đột gay gắt về tài nguyên nước.
Khoảng trống nghiên cứu (Research gap) cốt lõi được nhận diện: Trước đây, các nghiên cứu tại Đồng bằng sông Cửu Long (ĐBSCL) chủ yếu tiếp cận theo hướng mô phỏng thủy lực hoặc chất lượng nước đơn lẻ mà chưa tích hợp được công cụ phân tách nguồn gốc đóng góp của từng loại hình xả thải (sinh hoạt, công nghiệp, trồng trọt, nuôi trồng thủy sản) cũng như chưa lượng hóa được cơ chế tương tác giữa các nguồn nước biên (nguồn ngọt sông Hậu, biên mặn Biển Đông, biên mặn Biển Tây) với các nguồn nước thải nội vùng. Luận án đặt ra 3 câu hỏi nghiên cứu và hệ giả thuyết khoa học:
- RQ1: Tải lượng và động thái phát thải của các nguồn thải điểm và nguồn thải diện phân tán $\ge 50\text{ m}^3/\text{ngày đêm}$ trên địa bàn Bán đảo Cà Mau (BĐCM) biến động như thế nào qua các mốc thời gian 2016, 2020 và viễn cảnh 2030?
- RQ2: Quy luật lan truyền không gian - thời gian của các thành phần nguồn nước từ biên và nguồn ô nhiễm nội vùng diễn ra như thế nào dưới tác động của thủy triều hai biển và hạ tầng kiểm soát nước?
- RQ3: Mức độ đóng góp tương đối của từng nhóm nguồn thải đến sự suy giảm các chỉ tiêu oxy hòa tan ($\text{DO}$), nhu cầu oxy sinh hóa ($\text{BOD}_5$), amoni ($\text{NH}_4^+$) và chỉ số chất lượng nước ($\text{WQI}$) tại từng tiểu vùng thủy lợi là bao nhiêu?
Khung lý thuyết của nghiên cứu được xây dựng trên nền tảng Lý thuyết lan truyền các thành phần nguồn nước kết hợp cùng hệ phương trình thủy động lực học Saint-Venant 1D, phương trình tải - khuếch tán Advection-Dispersion ($\text{AD}$), động học sinh hóa trong mô hình $\text{Ecolab}$, kỹ thuật trung bình mô hình Bayes ($\text{Bayesian Model Averaging - BMA}$) và các thuật toán học máy nâng cao. Phạm vi không gian nghiên cứu bao phủ toàn bộ BĐCM với diện tích tự nhiên khoảng 1,6 triệu ha, thuộc 6 tỉnh/thành phố (Cần Thơ, Hậu Giang, Sóc Trăng, Bạc Liêu, Cà Mau, một phần Kiên Giang), dân số 6.680.174 người (mật độ 408 người/$\text{km}^2$), lấy năm 2016 làm năm điển hình cho chu kỳ hạn - mặn lịch sử khi dòng chảy thượng nguồn sông Mê Công suy giảm tới 30% so với trung bình nhiều năm.
Literature Review và Positioning
Tổng quan nghiên cứu quốc tế chỉ ra rằng việc kiểm soát ô nhiễm đa nguồn trên lưu vực hạ lưu sông và vùng cửa biển là thách thức toàn cầu:
- Tại Trung Quốc, Chen et al. và các báo cáo môi trường chỉ ra việc dư thừa 40% phân bón hóa học đã xả thải trực tiếp hàng năm 10 triệu tấn dinh dưỡng chưa xử lý, làm 41,7% trong tổng số 140.000 km sông bị ô nhiễm phú dưỡng nghiêm trọng.
- Tại Ấn Độ, nghiên cứu lưu vực sông Hằng cho thấy 80% rác và nước thải đô thị đổ thẳng ra nguồn tiếp nhận, làm suy giảm 4% GDP quốc gia do tổn thất môi trường, tước đi khả năng cấp nước sinh hoạt và tưới tiêu trong mùa kiệt.
- Nghiên cứu của Ngân hàng Phát triển Châu Á (ADB, 2011) trên sông Citarum (Indonesia) – nơi cấp 80% nước cho 14 triệu dân Jakarta – chỉ ra áp lực từ 2.000 nhà máy và dư lượng thuốc bảo vệ thực vật (BVTV) khiến dòng sông trở thành "bãi rác di động".
- Nghiên cứu của Evans et al. (2012) trên sông Marilao (Philippines) khẳng định các kim loại nặng và chất hữu cơ từ chế xuất công nghiệp đã hủy hoại nghiêm trọng hệ sinh thái ven bờ vịnh Manila.
Tại Việt Nam, các mô hình toán truyền thống như DELTA (Nguyễn Tất Đắc), KOD-WQPS (Nguyễn Ân Niên), VRSAP (Nguyễn Như Khuê), HYDROGIS (Nguyễn Hữu Nhân), F28 (Lê Song Giang) đã đóng góp lớn cho bài toán dòng chảy và xâm nhập mặn. Chu Thái Hoành (2001, 2005) đã ứng dụng GIS và mô hình thủy lực phân tích xung đột giữa nông nghiệp và nuôi trồng thủy sản (NTTS). Timothy (2007) và Đoàn Quang Trí (2018) đã chỉ ra hiện tượng suy thoái chất lượng nước tại các cửa sông ven biển BĐCM do các thông số $\text{DO}, \text{BOD}_5, \text{NH}_4^+$ vượt quy chuẩn kỹ thuật quốc gia QCVN 08-MT:2015/BTNMT nhiều lần.
Về mặt học thuật, tồn tại cuộc tranh luận lớn giữa hai trường phái: (1) Trường phái sinh thái thủy văn truyền thống coi dòng tiếp nhận là hệ đệm có khả năng tự làm sạch tự nhiên và tập trung kiểm soát nồng độ cuối nguồn; (2) Trường phái quản lý tích hợp tài nguyên nước đòi hỏi phải bóc tách tỷ lệ đóng góp của từng thành phần dòng chảy ($P_i%$) để áp đặt trách nhiệm xả thải theo tải lượng. Luận án định vị nghiên cứu vào trường phái thứ hai, lấp đầy khoảng trống khi các nghiên cứu trước đây tại BĐCM chưa tách bạch được tác động lan truyền của từng nguồn thải cụ thể dưới chế độ trao đổi nước phức tạp của hệ thống kênh mương đan xen dày đặc nối liền hai chế độ triều đối lập.
Đóng góp lý thuyết và khung phân tích
Đóng góp cho lý thuyết
Luận án kế thừa và phát triển mở rộng Lý thuyết lan truyền các thành phần nguồn nước (Water Source Component Transport Theory) do GS.TS. Nguyễn Ân Niên và GS.TS. Tăng Đức Thắng khởi xướng (1999, 2002). Luận án đã mở rộng khung lý thuyết này từ bài toán lan truyền khối nước tự nhiên thuần túy sang bài toán lan truyền chất ô nhiễm gắn với nguồn thải nhân sinh theo phương trình cân bằng nồng độ thể tích:
$$\frac{\partial (A \cdot P_i)}{\partial t} + \frac{\partial (Q \cdot P_i)}{\partial x} - \frac{\partial}{\partial x}\left(A \cdot D \cdot \frac{\partial P_i}{\partial x}\right) = q_i$$
Trong đó: $P_i$ là tỷ lệ thành phần nguồn nước thứ $i$ ($0 \le P_i \le 100%$, $\sum P_i = 100%$); $A$ là diện tích mặt cắt ngang lòng dẫn; $Q$ là lưu lượng toàn dòng; $D$ là hệ số khuếch tán; $q_i$ là lưu lượng dòng gia nhập bên của nguồn nước thành phần $i$. Luận án chứng minh một bước chuyển biến nhận thức (paradigm shift): Chất lượng nước tại một mặt cắt trong mạng lưới kênh triều không chỉ phụ thuộc vào nồng độ xả thải tại chỗ mà bị chi phối bởi cấu trúc tỷ lệ phần trăm các khối nước thành phần hội tụ tại điểm đó.
+-----------------------------------------------------------------------------------+
| KHUNG PHÂN TÍCH TÍCH HỢP CỦA LUẬN ÁN |
+-----------------------------------------------------------------------------------+
|
+---------------------------------+---------------------------------+
| |
v v
+-------------------------------+ +---------------------------------+
| 1. ĐIỀU KIỆN BIÊN THỦY LỰC | | 2. INVENTORY NGUỒN XẢ THẢI |
| - Biên sông Hậu (Nước ngọt) | | - Nước thải sinh hoạt đô thị |
| - Biển Đông (Bán nhật triều) | | - Nước thải công nghiệp chế biến|
| - Biển Tây (Nhật triều) | | - Nông nghiệp (Lúa, hoa màu) |
+-------------------------------+ | - Nuôi trồng thủy sản (Tôm) |
| +---------------------------------+
+---------------------------------+---------------------------------+
|
v
+---------------------------------------------------+
| 3. BỘ CÔNG CỤ TÍCH HỢP MÔ PHỎNG |
| - MIKE 11-HD: Mô hình dòng chảy thủy lực Saint- |
| Venant 1D trên toàn mạng kênh BĐCM |
| - MIKE 11-AD & TPN: Lan truyền thành phần nguồn |
| nước và tải lượng các chất ô nhiễm |
| - MIKE 11-Ecolab: Mô phỏng động thái biến đổi |
| hóa sinh của DO, BOD5, NH4+, N, P |
+---------------------------------------------------+
|
+---------------------------------+---------------------------------+
| |
v v
+-------------------------------+ +---------------------------------+
| 4. TỐI ƯU HÓA BIẾN (BMA) | | 5. HỌC MÁY DỰ BÁO CHỈ SỐ WQI |
| - Bayesian Model Averaging | | - Cây quyết định (Decision Tree)|
| - Tiêu chuẩn BIC & Bayes | | - Thuật toán tăng cường |
| Factor (BF) khử biến dư | | (Gradient Boosting/Ensemble) |
+-------------------------------+ +---------------------------------+
|
+---------------------------------+
|
v
+-----------------------------------+
| 6. ĐÁNH GIÁ TÁC ĐỘNG & BẢN ĐỒ HÓA |
| - Không gian phân bố tải lượng |
| - Bản đồ nguy cơ suy thoái nước |
| - Giải pháp công trình/phi công |
| trình kiểm soát theo tiểu vùng |
+-----------------------------------+
Khung phân tích độc đáo
Khung phân tích của luận án tích hợp liên hoàn 3 lý thuyết nền tảng:
- Lý thuyết Thủy lực Mạng sông kênh hở: Mô tả trường vận tốc $v(x,t)$ và mực nước $Z(x,t)$ theo hệ phương trình vi phân phi tuyến Saint-Venant 1 chiều.
- Lý thuyết Động học Môi trường Nước (Ecolab Kinetics): Mô phỏng chu trình chuyển hóa các hợp chất hữu cơ và dinh dưỡng ($\text{BOD} - \text{DO} - \text{Nitrogen cycle}$).
- Lý thuyết Xác suất Thống kê Bayes và Học máy: Tích hợp Bayesian Model Averaging ($\text{BMA}$) thông qua chỉ số Schwarz Bayesian Information Criterion ($\text{BIC}$) để loại bỏ hiện tượng đa cộng tuyến và thừa biến trong chuỗi số liệu quan trắc, từ đó nạp vào mô hình Cây quyết định (Decision Tree) và thuật toán Tăng cường (Boosting) để dự báo chỉ số $\text{WQI}$.
Điều kiện biên lý thuyết xác định rõ: Khung mô hình áp dụng cho hệ thống kênh dẫn hở vùng đồng bằng triều có độ dốc đáy nhỏ, không xét hiện tượng phân tầng mật độ thẳng đứng nghiêm trọng trong lòng kênh và giả định chất ô nhiễm hòa tan hoàn toàn trong pha nước.
Phương pháp nghiên cứu tiên tiến
Thiết kế nghiên cứu
Nghiên cứu được thiết kế theo trường phái Thực chứng luận kết hợp Hiện thực phản biện (Critical Realism), ứng dụng phương pháp hỗn hợp quy mô lớn (Multi-level Mixed-Methods Design). Thiết kế nghiên cứu phân cấp theo 3 cấp độ:
- Cấp độ 1 (Vĩ mô): Toàn vùng BĐCM (1,6 triệu ha) chịu chi phối bởi biên Sông Hậu, Biển Đông và Biển Tây.
- Cấp độ 2 (Trung mô): Phân chia 7 tiểu vùng sinh thái - thủy lợi và 120 ô thủy lợi khép kín/hở.
- Cấp độ 3 (Vi mô): Các điểm xả thải cụ thể có lưu lượng $Q_{nt} \ge 50\text{ m}^3/\text{ngày đêm}$ dọc theo các trục kênh chính (Quản Lộ - Phụng Hiệp, Cà Mau - Bạc Liêu, Xà No, sông Gành Hào, sông Ông Đốc).
+-----------------------------------------------------------------------------------+
| SƠ ĐỒ PHÂN CẤP THIẾT KẾ NGHIÊN CỨU VÙNG BÁN ĐẢO CÀ MAU |
+-----------------------------------------------------------------------------------+
[ CẤP ĐỘ 1: VĨ MÔ - TOÀN BỘ BÁN ĐẢO CÀ MAU (1,6 TRIỆU HA) ]
[ CẤP ĐỘ 2: TRUNG MÔ - 7 TIỂU VÙNG THỦY LỢI & 120 Ô THỦY LỢI ]
[ CẤP ĐỘ 3: VI MÔ - NGUỒN XẢ THẢI ĐIỂM & DIỆN (Q_nt >= 50 m3/ngày đêm) ]
Quy trình nghiên cứu rigorous
Quy trình thu thập và xử lý số liệu tuân thủ nghiêm ngặt các quy chuẩn kỹ thuật quốc gia:
- Quan trắc và phân tích: Lấy mẫu theo chuỗi thời gian mùa khô (tháng 4) và mùa mưa (tháng 10/2016) trên toàn mạng lưới trạm quan trắc của các Sở Tài nguyên và Môi trường Cần Thơ, Hậu Giang, Sóc Trăng, Bạc Liêu, Cà Mau, Kiên Giang. Các chỉ tiêu đo đạc gồm $\text{pH}, \text{DO}, \text{TSS}, \text{BOD}_5, \text{COD}, \text{NH}_4^+, \text{NO}_2^-, \text{NO}_3^-, \text{PO}_4^{3-}, \text{Coliform}$ theo đúng phương pháp chuẩn Standard Methods for the Examination of Water and Wastewater (APHA).
- Tam giác đạc dữ liệu (Data Triangulation): Kết hợp đồng thời 3 nguồn dữ liệu: (1) Số liệu đo đạc thực tế tại hiện trường; (2) Số liệu thống kê niên giám phát triển kinh tế - xã hội, định mức sử dụng phân bón, thuốc BVTV; (3) Kết quả viễn thám GIS giải đoán diện tích mặt nước nuôi tôm và đất canh tác lúa.
- Độ tin cậy của mô hình: Bộ thông số mô hình Mike11-HD được hiệu chỉnh và kiểm định với hệ số Nash-Sutcliffe ($NSE$) đạt từ $0{,}82$ đến $0{,}91$, sai số biên đỉnh triều $< 5\text{ cm}$, pha triều sai lệch $< 30\text{ phút}$.
Data và phân tích
Số liệu kiểm kê nguồn thải phát sinh cho thấy áp lực khổng lồ lên hệ thống thủy vực:
- Nước thải sinh hoạt: Năm 2016 phát sinh hàng trăm nghìn $\text{m}^3/\text{ngày đêm}$, chỉ $10%$ lượng nước thải đô thị được xử lý sơ bộ qua trạm tập trung, $100%$ các khu dân cư nông thôn và ven đô xả thẳng trực tiếp vào kênh rạch.
- Nước thải công nghiệp: 49 nguồn lớn (34 KCN và 15 cụm công nghiệp), tỷ lệ xử lý đạt chuẩn tại các KCN chỉ đạt khoảng $60%$. 5 nhà máy đường có lưu lượng trung bình $340\text{ m}^3/\text{ngày đêm/nhà máy}$ xả tập trung trong vụ ép mía (tháng 9 đến tháng 3 năm sau).
- Nước thải nông nghiệp và hóa chất: Suất sử dụng phân bón vô cơ tại BĐCM cao hơn $35{,}3%$ so với mức trung bình toàn quốc; suất sử dụng thuốc BVTV hóa học cao hơn $71{,}9%$ (cá biệt có tỉnh cao gấp 3 lần), trung bình nông dân phun thuốc trừ sâu 5 lần/vụ lúa. Ước tính có khoảng 70.000 kg và hơn 40.000 lít thuốc BVTV cùng 70.000 kg vỏ bao bì hóa chất thải bỏ ra môi trường mỗi năm; lượng dư thừa $\text{NPK}$ không được cây trồng hấp thụ đạt từ $50 - 70%$, bị rửa trôi trực tiếp ra sông kênh.
- Nước thải nuôi trồng thủy sản: Các đợt xả nước ao tôm và bùn đáy định kỳ (5–10 kỳ/vụ nuôi 4–5 tháng) mang theo hàm lượng hữu cơ, amoni, mầm bệnh cực lớn xả tự do vào cùng tuyến kênh cấp nước do thiếu hệ thống kênh cấp - thoát riêng biệt.
Xử lý thống kê và học máy: Sử dụng mô hình BMA nhằm phân tích không gian xác suất của các mô hình hồi quy để chọn lựa tập hợp biến độc lập có ảnh hưởng trọng yếu đến $\text{WQI}$. Kết quả lựa chọn các thông số cốt lõi gồm $\text{DO}, \text{BOD}_5, \text{NH}_4^+, \text{TSS}, \text{Coliform}$. Bốn mô hình học máy (Decision Tree, Random Forest, Gradient Boosting, Support Vector Machine) được thử nghiệm đối chứng trên tập dữ liệu train/test ($80/20$), kết quả thuật toán Cây quyết định tăng cường cho độ chính xác dự báo $\text{WQI}$ cao nhất với $R^2 > 0{,}88$ và $RMSE < 4{,}5$.
Phát hiện đột phá và implications
Những phát hiện then chốt
Luận án mang lại 5 phát hiện cốt lõi có giá trị khoa học và thực tiễn sâu sắc:
+------------------------------------------------------------------------------------+
| CÁC PHÁT HIỆN ĐỘT PHÁ CỦA LUẬN ÁN VÙNG BÁN ĐẢO CÀ MAU |
+------------------------------------------------------------------------------------+
1. NGUỒN PHÂN TÁN CHIẾM TỶ TRỌNG ÁP ĐẢO VÀ NGUY CƠ NHIỄM BẨN CHÉO
2. HIỆN TƯỢNG "NƯỚC ĐỨNG" VÀ Ô NHIỄM DO XUNG ĐỘT TRIỀU HAI BIỂN
3. SUY THOÁI CHẤT LƯỢNG NƯỚC NGHIÊM TRỌNG TẠI CÁC ĐÔ THỊ VÀ VÙNG CHUYỂN ĐỔI
4. DỊCH CHUYỂN BẤT LỢI CỦA KHỐI NƯỚC NGỌT BIÊN SÔNG HẬU TRONG MÙA KHÔ
5. ĐỘC TÍNH VÀ THIỆT HẠI KINH TẾ QUY MÔ LỚN TRONG THỦY SẢN
- Sự thống trị của tải lượng ô nhiễm phân tán và cơ chế tự gây nhiễm (Self-contamination): Trái với quan niệm phổ biến cho rằng ô nhiễm nước chủ yếu do công nghiệp, kết quả tính toán định lượng chỉ ra nguồn thải phân tán từ nông nghiệp và NTTS chiếm tỷ trọng áp đảo về lưu lượng và tải lượng hữu cơ. Do hạ tầng thủy lợi thiết kế phục vụ trồng lúa trước đây, các hộ nuôi tôm sử dụng chung kênh cấp và kênh tiêu, dẫn đến hiện tượng nước thải xả từ ao này ngay lập tức bị hút vào làm nước cấp cho ao lân cận, lan truyền dịch bệnh trên diện rộng (gây thiệt hại 12.410 ha tôm năm 2018).
- Cơ chế lan truyền tù đọng do vùng giao thoa triều hai biển (Tidal Stagnation Zone): Dưới tác động của bán nhật triều Biển Đông (biên độ lớn 3,0–4,0 m) truyền qua các cửa sông Mỹ Thanh, Gành Hào và nhật triều Biển Tây (biên độ nhỏ 0,8–1,2 m) truyền qua sông Cái Lớn, sông Ông Đốc, vùng trung tâm BĐCM (khu vực Ngã Năm, Phước Long, Hồng Dân) hình thành các "vùng nước đứng". Tại đây, vận tốc dòng chảy triệt tiêu, tỷ lệ nguồn nước thải NTTS và sinh hoạt lưu giữ lại trong hệ thống sau 7 đến 30 ngày xả thải vẫn duy trì ở mức cao ($> 20 - 40%$), không thể tự thau rửa ra biển.
- Mức độ suy thoái chất lượng nước vượt ngưỡng chịu tải sinh thái: Toàn bộ hệ thống sông kênh chính tiếp nhận nước thải đô thị và công nghiệp (kênh Xà No, sông Gành Hào, sông Đốc, kênh Bạc Liêu – Cà Mau, Quản Lộ – Phụng Hiệp) có các chỉ tiêu $\text{BOD}_5, \text{COD}, \text{NH}_4^+, \text{Coliform}$ vượt QCVN 08-MT:2015 từ 2 đến 12 lần. Đánh giá theo chỉ số $\text{WQI}$ năm 2020 cho thấy: Sóc Trăng có 10/19 vị trí quan trắc $\text{WQI} \le 50$; Bạc Liêu có 6/8 vị trí $\text{WQI} \le 70$ (trong đó cống Cái Cùng và cống Đầu Bằng $\text{WQI} \le 25$ - ô nhiễm rất nặng); Cà Mau có 10/20 vị trí $\text{WQI} \le 50$, chỉ duy nhất 1 vị trí cửa sông Cửa Lớn đạt mức khá.
- Quy luật lan truyền thành phần nước ngọt Sông Hậu suy giảm trong mùa khô: Khi xả thải diễn ra trong mùa kiệt, tỷ lệ thành phần nước ngọt Sông Hậu ($P_{\text{Hậu}}$) tại các vùng cuối hệ thống ngọt hóa (như vùng Long Phú - Tiếp Nhật, Nam Mang Thít) giảm xuống dưới $15%$, thay vào đó là thành phần nước biển mặn xâm nhập ($P_{\text{biển}} > 80%$). Việc đóng cống ngăn mặn kéo dài nhiều tháng biến hệ thống kênh nội đồng thành các "túi chứa nước thải khép kín", làm chỉ số $\text{DO}$ giảm sâu xuống $< 2\text{ mg/l}$, bùng phát hiện tượng phú dưỡng và ô nhiễm hữu cơ nội tại.
- Đột biến ô nhiễm theo mùa vụ công nghiệp và nông nghiệp: Các đợt xả thải của 5 nhà máy mía đường (lưu lượng $340\text{ m}^3/\text{ngày}$) vào giai đoạn tháng 9 đến tháng 3 năm sau tạo nên những "vệt đen ô nhiễm" kéo dài hàng chục km (điển hình sự cố sông Cái Nước Đục tại thị xã Long Mỹ tháng 5/2019 làm tê liệt toàn bộ nhà máy nước sinh hoạt).
Implications đa chiều
- Về lý thuyết: Làm giàu thêm lý thuyết thủy lực - chất lượng nước vùng đồng bằng châu thổ sông Mê Công thông qua việc tích hợp định lượng thành phần nguồn nước vào chu trình động học các chất ô nhiễm.
- Về phương pháp luận: Thiết lập quy trình phối hợp hiệu quả giữa mô hình hóa số trị 1D phân tán ($\text{MIKE 11}$) với phương pháp chọn lọc biến Bayes ($\text{BMA}$) và thuật toán học máy, mở ra công cụ chi phí thấp nhưng độ chính xác cao trong quan trắc và dự báo chất lượng nước.
- Về thực tiễn kỹ thuật: Cung cấp cơ sở khoa học để tái cơ cấu quy trình vận hành đóng/mở hệ thống cống liên tỉnh (Quản Lộ - Phụng Hiệp, cống Cái Lớn - Cái Bé) theo thời gian thực: vận hành mở cống tạo dòng chảy một chiều để "thau chua, rửa bẩn", đẩy khối nước ô nhiễm ra phía biển Tây hoặc biển Đông vào các chân triều thấp.
- Về chính sách: Đề xuất lộ trình bắt buộc tách rời mạng lưới thủy lợi cấp - thoát nước cho các vùng nuôi tôm công nghiệp thâm canh tại Sóc Trăng, Bạc Liêu, Cà Mau và Kiên Giang; ban hành quy chuẩn kỹ thuật địa phương siết chặt tiêu chuẩn xả thải tại chỗ cho các cơ sở chế biến thủy sản ngoài khu công nghiệp.
Limitations và Future Research
Luận án thẳng thắn thừa nhận các giới hạn nghiên cứu:
- Giới hạn số liệu và đối tượng xả thải: Chỉ tập trung mô hình hóa chi tiết các nguồn thải điểm có lưu lượng $Q_{nt} \ge 50\text{ m}^3/\text{ngày đêm}$. Hàng nghìn điểm xả thải phân tán quy mô nhỏ dưới ngưỡng (nông hộ, tàu thuyền giao thông thủy) được gộp chung vào lưu lượng bên ($q$) mà chưa phân tách chi tiết tọa độ.
- Bỏ qua động học dòng chảy tràn nước mưa: Chưa xét đến quá trình rửa trôi ô nhiễm mặt đất trong các trận mưa rào đầu mùa mưa – thời điểm bùng phát axit hóa từ đất phèn tiềm tàng ($\text{Acid Sulfate Soils}$) và cuốn trôi lượng lớn dư lượng thuốc BVTV.
- Giới hạn không gian mô hình 1D: Hệ thống Mike11 1 chiều giả định chất ô nhiễm hòa trộn đều trên mặt cắt ngang, chưa phản ánh được hiện tượng phân tầng mặn - ngọt hoặc phân tầng nhiệt độ/oxy hòa tan tại các khúc sông sâu và khu vực cửa biển ngập triều.
Chương trình nghiên cứu tiếp theo (Future Research Agenda):
- Mở rộng phát triển mô hình ghép nối 1D-2D-3D ($\text{MIKE 11 - MIKE 21/3D}$) tại các khu vực cửa sông quan trọng (Gành Hào, Mỹ Thanh, Cửa Lớn, Cái Lớn) để mô phỏng chính xác quá trình lan truyền chất ô nhiễm trong không gian đa chiều.
- Tích hợp module sinh thái học đất phèn để lượng hóa đồng thời quá trình giải phóng ion nhôm ($\text{Al}^{3+}$), sắt ($\text{Fe}^{2+}/\text{Fe}^{3+}$) và sunfat ($\text{SO}_4^{2-}$) tương tác với các nguồn xả thải hữu cơ.
- Ứng dụng công nghệ cảm biến viễn thám kết hợp mạng lưới Internet vạn vật ($\text{IoT}$) để cập nhật liên tục dữ liệu đầu vào thời gian thực cho mô hình dự báo học máy.
- Đánh giá rủi ro sức khỏe sinh thái và sức khỏe con người ($\text{Ecological and Human Health Risk Assessment}$) theo chuỗi thức ăn từ nguồn nước mặt ô nhiễm sang sản phẩm tôm - lúa xuất khẩu.
Tác động và ảnh hưởng
Nghiên cứu tạo ra tác động sâu rộng trên nhiều phương diện:
- Ảnh hưởng học thuật: Luận án là tài liệu tham khảo khoa học chuẩn mực cho các cơ sở nghiên cứu chuyên ngành Thủy văn, Thủy lực, Môi trường đất và nước; cung cấp bộ tham số mô hình hóa đáng tin cậy đã được hiệu chỉnh cho điều kiện thủy văn khắc nghiệt của vùng ĐBSCL.
- Chuyển đổi ngành công nghiệp và nông nghiệp: Cung cấp luận cứ khoa học để các tỉnh BĐCM tái cấu trúc quy hoạch vùng nuôi tôm công nghiệp, chuyển từ hình thức nuôi xả thải trực tiếp sang mô hình kinh tế tuần hoàn nước, ứng dụng công nghệ biofloc và ao lắng xử lý tuần hoàn khép kín.
- Hoạch định chính sách quản lý nhà nước: Đóng góp trực tiếp vào cơ sở dữ liệu của Bộ Nông nghiệp & PTNT và Bộ Tài nguyên & Môi trường trong việc xây dựng Quy hoạch tài nguyên nước lưu vực sông Cửu Long giai đoạn 2021–2030, tầm nhìn 2050; định hình quy trình vận hành liên hồ chứa và hệ sinh thái công trình thủy lợi Cái Lớn - Cái Bé.
- Lợi ích xã hội và hội nhập quốc tế: Giúp bảo vệ an ninh nguồn nước cấp cho hơn 6,6 triệu dân BĐCM, giảm thiểu nguy cơ bùng phát dịch bệnh thủy sản trị giá hàng trăm triệu USD mỗi năm, đồng thời đóng góp bài học thực tiễn giá trị cho Ủy hội sông Mê Công quốc tế ($\text{MRC}$) trong việc quản lý chất lượng nước xuyên biên giới.
Đối tượng hưởng lợi
- Nghiên cứu sinh và học viên cao học: Tiếp cận phương pháp luận bóc tách thành phần nguồn nước tích hợp mô hình số trị và thuật toán học máy giải quyết bài toán ô nhiễm sông ngòi vùng đồng bằng triều phức tạp.
- Các nhà khoa học chuyên ngành Tài nguyên nước & Môi trường: Kế thừa bộ cơ sở dữ liệu quan trắc, phân tích và tham số hiệu chỉnh mô hình thủy lực - chất lượng nước chuẩn hóa trên toàn BĐCM.
- Doanh nghiệp R&D ngành Thủy sản và Chế biến thực phẩm: Nắm bắt quy luật lan truyền và vùng phân bố ô nhiễm để lựa chọn vị trí lấy nước mặt an toàn, tối ưu hóa công nghệ xử lý nước cấp và nước thải.
- Cơ quan quản lý nhà nước (Sở NN&PTNT, Sở TN&MT các tỉnh ĐBSCL): Sở hữu công cụ dự báo và cảnh báo sớm nguy cơ ô nhiễm nguồn nước, xây dựng kế hoạch đóng mở cống ngăn mặn/tiêu ô nhiễm tối ưu theo từng mùa vụ sản xuất.
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ì?
Đóng góp lý thuyết độc đáo nhất là việc mở rộng và chuẩn hóa Lý thuyết lan truyền các thành phần nguồn nước (GS.TS Nguyễn Ân Niên và GS.TS Tăng Đức Thắng) để giải bài toán tích hợp giữa động lực học dòng chảy triều hai biển với động thái phân rã - chuyển hóa của các chất ô nhiễm từ 4 nhóm nguồn xả thải riêng biệt (sinh hoạt, công nghiệp, trồng trọt, NTTS). Luận án đã biến lý thuyết từ công cụ phân tích nguồn nước tự nhiên thành công cụ quản lý định lượng tải lượng môi trường cho từng đối tượng gây ô nhiễm.
2. Phương pháp luận của luận án có điểm gì đổi mới so với các nghiên cứu trước đây?
So với các nghiên cứu trước đây (như mô hình VRSAP của Chu Thái Hoành chỉ tập trung vào xâm nhập mặn và chua phèn, hay mô hình Mike11 truyền thống của Đoàn Quang Trí chỉ mô phỏng chất lượng nước tổng hợp), luận án đã:
- Thiết lập quy trình bóc tách từng nguồn thải thành phần $P_i%$ trong môi trường mô phỏng phân tán Mike11-AD.
- Tích hợp phương pháp trung bình mô hình Bayes ($\text{BMA}$) qua tiêu chuẩn $\text{BIC}$ để lọc sạch biến nhiễu trước khi đưa vào các thuật toán học máy (Cây quyết định, Boosting), giúp tăng tốc độ dự báo $\text{WQI}$ với độ chính xác cao ($R^2 > 0{,}88$).
3. Phát hiện bất ngờ nhất từ số liệu thực nghiệm là gì?
Phát hiện bất ngờ nhất là nguồn thải sinh hoạt đô thị và công nghiệp chế biến thủy sản không phải là nguyên nhân chính gây suy thoái diện rộng trên toàn bán đảo, mà chính hoạt động nuôi tôm thâm canh và canh tác lúa với suất dùng thuốc BVTV cao gấp 1,7 lần toàn quốc mới là tác nhân chi phối làm suy sụp chỉ số WQI. Đặc biệt, tại vùng tiếp giáp hai chế độ triều, các chất ô nhiễm không bị đẩy ra biển mà dao động lắc lư tại chỗ hơn 30 ngày, biến các kênh cấp nước thành nguồn phát tán mầm bệnh ngược lại cho các ao nuôi.
4. Luận án có cung cấp quy trình tái lập (Replication Protocol) không?
Có. Luận án mô tả chi tiết toàn bộ các bước thiết lập mạng lưới tính toán Mike11 với 120 ô thủy lợi, phương pháp xác định điều kiện biên thủy lực (mực nước triều Biển Đông tại trạm Gành Hào, Biển Tây tại trạm Rạch Giá, lưu lượng sông Hậu tại trạm Cần Thơ), hệ số nhám Manning $n$, bảng tham số phân rã động học trong Ecolab, quy trình trích xuất dữ liệu huấn luyện cho BMA và thuật toán học máy dự báo WQI trong các phụ lục chi tiết.
5. Định hướng nghiên cứu 10 năm tới được vạch ra như thế nào?
Định hướng tập trung vào: (1) Xây dựng hệ thống bản sao số (Digital Twin) mạng lưới sông kênh BĐCM tích hợp AI và viễn thám để điều hành cống tự động theo thời gian thực; (2) Nghiên cứu cơ chế lan truyền các chất ô nhiễm mới nổi (Microplastics, kháng sinh trong NTTS); (3) Đánh giá tác động tương hỗ giữa sụt lún đất do khai thác nước ngầm với diễn biến lan truyền ô nhiễm nước mặt.
Kết luận
- Lượng hóa toàn diện phát thải: Đã kiểm kê và mô hình hóa chi tiết 49 nguồn xả thải công nghiệp lớn, 5 nhà máy đường, hệ thống nước thải sinh hoạt đô thị cùng hàng triệu $\text{m}^3$ nước thải/bùn thải từ nông nghiệp và nuôi tôm công nghiệp trên phạm vi 1,6 triệu ha BĐCM.
- Làm sáng tỏ cơ chế thủy động lực ô nhiễm: Chứng minh quy luật tích tụ ô nhiễm hữu cơ và dinh dưỡng tại vùng giao thoa triều hai biển Biển Đông – Biển Tây, nơi chất ô nhiễm có thời gian lưu giữ kéo dài từ 7 đến trên 30 ngày.
- Đột phá công cụ mô hình hóa: Ứng dụng xuất sắc Lý thuyết lan truyền các thành phần nguồn nước tích hợp Mike11-Ecolab và phương pháp học máy Bayes, bóc tách chính xác tỷ lệ đóng góp của từng nguồn nước biên và nguồn thải nội đồng.
- Cảnh báo hiện trạng suy thoái nước nghiêm trọng: Chỉ ra các thông số $\text{BOD}_5, \text{COD}, \text{NH}_4^+, \text{Coliform}$ vượt QCVN từ 2 đến 12 lần, phần lớn sông kênh nội đồng có $\text{WQI} \le 50$ không đáp ứng tiêu chuẩn cấp nước sinh hoạt và an toàn nuôi trồng thủy sản.
- Đề xuất giải pháp kiểm soát nguồn thải thực thi: Đưa ra cơ sở khoa học để tái cơ cấu hệ thống công trình thủy lợi, vận hành hệ thống cống một chiều đẩy ô nhiễm và quy hoạch tách biệt triệt để kênh cấp - kênh thoát cho vùng nuôi tôm Bán đảo Cà Mau.
