Tổng quan về luận án

Nghiên cứu của nghiên cứu sinh Nguyễn Toàn Thắng (2018) tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội, dưới sự hướng dẫn khoa học của PGS. Lê Văn Thiện và PGS. Trần Hồng Thái, với tiêu đề "Đánh giá chất lượng nước, tính toán khả năng chịu tải của sông Nhuệ - Đáy trong bối cảnh biến đổi khí hậu và phát triển kinh tế xã hội" (Chuyên ngành: Môi trường đất và nước, Mã số: 62440303) là một công trình khoa học tiên phong giải quyết bài toán suy thoái môi trường nước mặt tại vùng kinh tế trọng điểm Bắc Bộ.

Về bối cảnh thực tiễn và tính cấp thiết, lưu vực sông Nhuệ – Đáy nằm ở hữu ngạn sông Hồng với diện tích tự nhiên khoảng 7.665 km² (trong đó diện tích tính toán thủy văn trực tiếp là 7.388 km²), dân số năm 2015 đạt xấp xỉ 12 triệu người cùng mật độ dân số trung bình 1.160 người/km², trải dài trên địa bàn 5 tỉnh, thành phố: Hà Nội, Hà Nam, Hòa Bình, Nam Định và Ninh Bình. Áp lực từ quá trình đô thị hóa và công nghiệp hóa nhanh chóng đã biến hệ thống sông này thành nơi tiếp nhận lượng nước thải khổng lồ. Dữ liệu khảo sát của luận án chỉ ra rằng toàn lưu vực có "khoảng 1.662 nguồn thải là cơ sở sản xuất, kinh doanh; 39 nguồn thải là khu công nghiệp, cụm công nghiệp; 137 cơ sở y tế và 144 làng nghề", phát thải trung bình khoảng 621.000 m³/ngày đêm phần lớn chưa qua xử lý đạt chuẩn.

Khoảng trống nghiên cứu (research gap) cốt lõi được xác định rõ ràng: Các công trình trước đây (như Nguyễn Văn Cư và nnk, 2005; Trần Hồng Thái và nnk, 2009; Nguyễn Thanh Sơn và nnk, 2011; Nguyễn Minh Lâm và nnk, 2012) hoặc chỉ tiếp cận đơn lẻ việc mô phỏng chất lượng nước theo hiện trạng thủy lực, hoặc chỉ đánh giá biến đổi dòng chảy do biến đổi khí hậu mà chưa tích hợp đồng bộ giữa mô hình mưa – dòng chảy, mô hình thủy lực 1 chiều không ổn định, mô hình sinh thái phân hủy sinh hóa với cơ chế tính toán sức chịu tải theo các kịch bản phát triển kinh tế – xã hội và kịch bản phát thải khí nhà kính trung bình RCP4.5.

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: Hiện trạng ô nhiễm không gian – thời gian và động lực diễn biến các thông số $DO$, $BOD_5$, $COD$, $NH_4^+$ trên trục chính sông Nhuệ – Đáy diễn ra theo quy luật nào?
  • RQ2: Dòng chảy kiệt và khả năng đồng hóa tự làm sạch của lưu vực biến thiên ra sao dưới tác động kép của biến đổi khí hậu (kịch bản RCP4.5 đến năm 2020, 2030) và các quy hoạch phát triển kinh tế – xã hội?
  • RQ3: Ngưỡng chịu tải tối đa ($TMDL$) và khả năng tiếp nhận nước thải của từng phân đoạn sông là bao nhiêu để đáp ứng các mục tiêu chất lượng nước theo QCVN 08-MT:2015 và Thông tư 02/2009/TT-BTNMT?
  • H1: Cơ chế tự làm sạch của hệ thống sông Nhuệ – Đáy trong mùa kiệt bị suy giảm nghiêm trọng do sự suy kiệt dòng chảy tự nhiên từ sông Hồng qua cống Liên Mạc và tình trạng thiếu hụt oxy hòa tan cục bộ.
  • H2: Tác động của biến đổi khí hậu theo kịch bản RCP4.5 làm gia tăng tính cực đoan của chế độ thủy văn (lũ tăng, kiệt giảm sâu), làm thu hẹp biên độ tiếp nhận tải lượng ô nhiễm từ 15% đến 25% trong mùa khô.
  • H3: Giải pháp kiểm soát nguồn thải đạt tỷ lệ xử lý 80% tại nguồn là điều kiện tiên quyết để phục hồi khả năng chịu tải dương cho mục đích cấp nước tưới tiêu và sinh hoạt.

Khung lý thuyết của công trình tích hợp lý thuyết cân bằng khối lượng, lý thuyết lan truyền khuếch tán Fickian, hệ phương trình thủy động lực học Saint-Venant 1 chiều, động học phân hủy sinh hóa các hợp chất hữu cơ – dinh dưỡng (mở rộng từ mô hình Streeter-Phelps) và khung đánh giá tải lượng xả thải tối đa theo ngày ($TMDL$ của US EPA). Luận án có ý nghĩa khoa học và thực tiễn to lớn, cung cấp luận cứ định lượng chính xác để các cơ quan quản lý phân bổ hạn ngạch xả thải, quy hoạch lưu vực và bảo vệ nguồn nước bền vững.

Literature Review và Positioning

Lịch sử phát triển của các nghiên cứu mô hình hóa chất lượng nước và đánh giá khả năng tự làm sạch của các thủy vực đã trải qua nhiều giai đoạn tiến hóa quan trọng. Khởi đầu từ năm 1925, Streeter và Phelps đã đặt nền móng lý thuyết đầu tiên với mô hình suy giảm oxy (oxygen sag curve) trên sông Ohio (Hoa Kỳ), mô phỏng sự tương tác giữa quá trình tiêu thụ oxy do phân hủy chất hữu cơ bậc nhất ($BOD$) và quá trình tái cấp khí bề mặt ($DO$). Trong các thập niên 1960–1970, các mô hình sinh thái nhiều chiều bắt đầu tích hợp quá trình lắng đọng bùn cát, quang hợp và hô hấp của thực vật thủy sinh (Chen & Orlob; Thomann & Fitzpatrick). Đến giai đoạn 1980–1990, các mô hình phức hợp như QUAL2E, QUAL2K, AQUATOX và MIKE 11 (DHI) ra đời, kết hợp chặt chẽ giữa thủy lực dòng chảy hở và động thái chuyển hóa dinh dưỡng nitơ, phốt pho (Thomann & Mueller, 1989; Law & Chalup, 1990).

Trong trường phái nghiên cứu tự làm sạch (self-purification) tại châu Âu, các công trình kinh điển từ đầu thế kỷ XX đã phân tích vai trò phân hủy của vi sinh vật và động học dòng chảy. Frank khi nghiên cứu sông Spree (Đức) đã chứng minh mật độ vi khuẩn tăng từ 10.000 khuẩn lạc/cm³ lên hàng triệu khuẩn lạc khi qua Berlin trước khi giảm mạnh nhờ quá trình tự làm sạch; Goldschmidt và Pransmidt tính toán khoảng cách tự làm sạch trên sông Isar (Munich) là 27 km trong thời gian 8 giờ; Kluse và Lossen xác định trên sông Rhine là 25–26 km trong 5,5–6 giờ. Tại Hoa Kỳ, Jordan (1900) đã phát hiện mối tương quan chặt chẽ giữa nhiệt độ môi trường nước và tốc độ tự làm sạch trên sông Illinois; trong khi các nghiên cứu tại sông Teshio (Nhật Bản) của Seki và Ebara (1987) chỉ ra sự phân tầng mức độ tự làm sạch theo các lớp nước trước khi đổ ra biển.

Trong bức tranh học thuật đó, tồn tại hai luồng quan điểm đối thoại sâu sắc:

  1. Trường phái thủy động lực học thuần túy: Cho rằng quá trình tự làm sạch chủ yếu do hiện tượng xáo trộn, pha loãng và truyền tải khuếch tán cơ học của khối nước chi phối, phụ thuộc trực tiếp vào lưu lượng và vận tốc dòng chảy.
  2. Trường phái sinh thái – sinh hóa học: Nhấn mạnh vai trò đồng hóa, hấp thụ và chuyển hóa chất ô nhiễm của quần xã vi sinh vật, thực vật thủy sinh (vi tảo, bèo tây, bèo cái) và động vật đáy, biến đổi chất hữu cơ và kim loại nặng thông qua các chu trình khoáng hóa và lắng đọng trầm tích.

Định vị nghiên cứu của luận án nằm ở điểm giao thoa giữa hai trường phái này trên một hệ thống thủy văn nhiệt đới gió mùa phức tạp chịu tác động mạnh của chế độ thủy triều và điều tiết nhân tạo. So với hai nghiên cứu quốc tế điển hình:

  • Mô hình EPDRiv1 trên lưu vực sông Coosa và sông Chattahoochee (Hoa Kỳ): Ứng dụng xác định $TMDL$ trong điều kiện mạng lưới dòng chảy tương đối ổn định và hệ thống dữ liệu quan trắc tự động dày đặc.
  • Nghiên cứu sông Selenga đổ vào hồ Baikal (Nga) của Makushkin (1998): Tập trung vào cơ chế chuyển hóa sinh hóa của vi sinh vật trong điều kiện khí hậu ôn đới lạnh.

Luận án của Nguyễn Toàn Thắng đã tạo bước đột phá khi giải quyết bài toán thiếu hụt chuỗi dữ liệu lưu lượng thực đo kéo dài ở các trạm nội lưu bằng cách kết nối chuỗi mô hình: sử dụng mô hình mưa – dòng chảy NAM để khôi phục biên thủy văn, lồng ghép vào mô hình thủy lực – chất lượng nước MIKE 11 HD-AD-Ecolab trên toàn bộ mạng lưới 8 con sông liên kết, và áp dụng công thức $TMDL$ kết hợp quy chuẩn Thông tư 02/2009/TT-BTNMT trong bối cảnh kịch bản biến đổi khí hậu toàn cầu.

Đó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à làm sâu sắc thêm lý thuyết sức chịu tải môi trường (environmental capacity) của GESAMP (1986) và khung lý thuyết $TMDL$ của Cơ quan Bảo vệ Môi trường Hoa Kỳ (US EPA) khi ứng dụng vào hệ thống sông ngòi đồng bằng sông Hồng. Tác giả đã mở rộng công thức tổng quát: $$\text{TMDL} = \sum \text{WLAs} + \sum \text{LAs} + \text{MOS}$$ (Trong đó $\text{WLAs}$ là tải lượng nguồn điểm, $\text{LAs}$ là tải lượng nguồn diện, $\text{MOS}$ là hệ số dự phòng an toàn).

Công trình đóng góp 3 mệnh đề lý thuyết (Theoretical Propositions - P) cốt lõi:

  • P1: Trong hệ thống sông phân lưu chịu ảnh hưởng của các công trình ngăn đập và suy giảm dòng chảy thượng nguồn, khả năng tự làm sạch vật lý (pha loãng – khuếch tán) suy giảm nhanh hơn tốc độ phân hủy sinh hóa, dẫn đến hiện tượng "tê liệt sinh thái" tại các đoạn sông tiếp nhận nước thải đô thị đậm đặc.
  • P2: Khả năng hấp thụ sinh học của hệ thực vật thủy sinh ($M_{VITAO}$, $M_{TVTS}$) và động vật thủy sinh ($M_{DVTS}$) tuân theo hàm toán học phụ thuộc trực tiếp vào sinh khối loài và nồng độ chất ô nhiễm, đóng vai trò như một hệ số đệm quan trọng trong phương trình cân bằng tải lượng: $$M_{TVTS} = \sum_{i=1}^{K} K_i \times f_i$$ (với $K_i$ là khối lượng loài thứ $i$, $f_i$ là khả năng hấp thụ của loài).
  • P3: Dưới tác động của biến đổi khí hậu theo kịch bản RCP4.5, sự phân bố lại lượng mưa và nhiệt độ nước làm phi tuyến hóa tốc độ phản ứng nitrat hóa ($NH_4^+ \to NO_2^- \to NO_3^-$) và hằng số tái cấp khí, làm dịch chuyển không gian của vùng lõm oxy ($DO\text{ sag}$) về phía hạ lưu.

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

Khung phân tích của luận án là một kiến trúc mô hình hóa 4 tầng tích hợp:

  1. Tầng mô hình mưa – dòng chảy (NAM): Mô phỏng liên tục hàm lượng ẩm trong 5 bể chứa thẳng đứng (bể mặt $U$, bể rễ cây $L$, bể nước ngầm tầng trên và tầng dưới) nhằm khôi phục chuỗi dữ liệu dòng chảy biên tại các trạm thủy văn Ba Thá (sông Đáy) và Hưng Thi (sông Bôi).
  2. Tầng mô hình thủy động lực 1 chiều (MIKE 11 HD): Giải hệ phương trình Saint-Venant mô tả dòng chảy không ổn định dọc theo mạng lưới sông phức hợp.
  3. Tầng mô hình truyền tải – khuếch tán và sinh thái (MIKE 11 AD & Ecolab): Vận hành phương trình truyền tải vi phân theo định luật Fick, kết hợp 6 cấp độ sinh thái mô phỏng 13 biến trạng thái chất lượng nước ($DO$, $BOD$, $COD$, $NH_4^+$, $NO_3^-$, $Coliform$, nhiệt độ).
  4. Tầng mô hình đánh giá khả năng tiếp nhận và sức chịu tải: Tính toán tải lượng tối đa cho phép ($L_{td}$) và khả năng tiếp nhận nước thải ($L_{tn}$) theo từng đoạn phân vùng chức năng nguồn nước (Cột A2, B1, B2 theo QCVN 08-MT:2015).

Điều kiện biên và phạm vi không gian được xác lập trên toàn bộ mạng lưới gồm 8 sông: sông Hồng (131 mặt cắt từ Sơn Tây đến Ba Lạt), sông Nhuệ (43 mặt cắt từ cống Liên Mạc đến Phủ Lý), sông Đáy (96 mặt cắt từ Ba Thá đến Như Tân), sông Hoàng Long (14 mặt cắt), sông Đào (10 mặt cắt), sông Ninh Cơ (25 mặt cắt), sông Châu Giang (11 mặt cắt) và sông Luộc (5 mặt cắt).

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

Thiết kế nghiên cứu

Luận án lựa chọn thế giới quan bản thể luận thực chứng (Positivism) kết hợp phương pháp nghiên cứu định lượng đa cấp độ (multi-level quantitative design). Quy trình nghiên cứu kết hợp chặt chẽ giữa quan trắc thực địa, phân tích hóa nghiệm phòng thí nghiệm tiêu chuẩn, xử lý dữ liệu không gian bằng hệ thông tin địa lý (GIS) và mô phỏng số trị động lực học.

Quy mô không gian bao quát diện tích 7.388 km², phân chia hệ thống sông Nhuệ – Đáy thành các đoạn sông chức năng đồng nhất dựa trên lưu lượng dòng chảy, vị trí các điểm nhập lưu và mục đích sử dụng nước theo quy hoạch phát triển kinh tế – xã hội của 5 tỉnh, thành phố.

Quy trình nghiên cứu rigorous

Chiến lược thu thập dữ liệu và quan trắc thực nghiệm được tổ chức đồng bộ trong các năm 2014 và 2015 qua hai chiến dịch lấy mẫu quy mô lớn vào mùa khô (tháng 5/2014 và tháng 5/2015) trên toàn tuyến:

  • Tuyến sông Nhuệ: Các vị trí trọng yếu gồm Cống Liên Mạc, Cổ Nhuế, Cầu Hà Đông, Phúc La, Thanh Liệt, Cự Đà, Đò Quan, Cống Thần, Đồ Kiều, Cống Nhật Tựu, Cầu Phủ Lý.
  • Tuyến sông Đáy: Các vị trí gồm Ba Thá, Tuế Tiêu, Vân Đình, Cầu Quế, Hồng Phú, Ngoại Độ, Trung Hiếu Hạ.

Phương pháp phân tích trong phòng thí nghiệm tuân thủ nghiêm ngặt các Tiêu chuẩn Quốc gia (TCVN) và thiết bị chuẩn hóa:

  • Nhiệt độ: Đo trực tiếp bằng máy đo đa chỉ tiêu điện tử HQ40d.
  • Oxy hòa tan ($DO$): Phương pháp chuẩn độ Winkler theo TCVN 7325:2004.
  • Nhu cầu oxy sinh hóa ($BOD_5$ ở 20°C): Phương pháp ủ 5 ngày theo TCVN 6001-1:2008.
  • Amoni ($NH_4^+$): Phương pháp trắc quang phổ kế theo TCVN 6491:1999.

Tính giá trị nội tại (internal validity) và độ tin cậy (reliability) của chuỗi số liệu được bảo đảm thông qua việc hiệu chuẩn chéo giữa số liệu quan trắc định kỳ của Tổng cục Môi trường giai đoạn 2011–2015 và các đợt khảo sát thực tế chuyên sâu của đề tài.

Data và phân tích

Dữ liệu khí tượng thủy văn được tổng hợp từ chuỗi đo đạc lịch sử giai đoạn 1972–2015 tại các trạm Sơn Tây, Hà Nội, Hưng Yên, Ba Thá, Hưng Thi, Phủ Lý, Ninh Bình, Như Tân, Ba Lạt. Phương pháp đa giác Thiessen được áp dụng để tính toán trọng số mưa bình quân lưu vực thượng nguồn sông Hoàng Long và thượng nguồn sông Đáy.

Quy trình hiệu chỉnh và kiểm nghiệm mô hình toán được tiến hành với độ chuẩn xác cao:

  • Mô hình NAM: Hiệu chỉnh bộ thông số ($U_{max}, L_{max}, CQOF, CQIF, CK_{1,2}, CKBF$) tại trạm Ba Thá (giai đoạn 1972–1973) và Hưng Thi (1973–1975); kiểm nghiệm tại Ba Thá (1976–1978) và Hưng Thi (1976–1977). Đánh giá độ trùng khớp giữa đường quá trình thực đo và tính toán bằng chỉ số hiệu quả Nash-Sutcliffe ($R^2$ đạt từ 0,78 đến 0,86), chứng minh độ tin cậy vượt trội trong việc khôi phục dòng chảy biên 2010–2015.
  • Mô hình thủy lực MIKE 11 HD: Hiệu chỉnh chuỗi mực nước mùa kiệt các năm 2010–2013 và kiểm nghiệm độc lập cho giai đoạn mùa kiệt 2013–2015 tại 4 trạm trọng điểm: Hà Nội, Hưng Yên, Phủ Lý, Ninh Bình. Kết quả sai số mực nước đỉnh và chân triều không vượt quá 5–10 cm.
  • Mô hình chất lượng nước MIKE 11 AD-Ecolab: Hiệu chỉnh với bộ thông số đo đạc thực tế tháng 5/2014 và kiểm định độc lập với chuỗi số liệu tháng 5/2015 cho các thông số $DO$, nhiệt độ, $BOD_5$ và $NH_4^+$, cho thấy sự tương quan chặt chẽ giữa đường tính toán mô phỏng và số liệu thực nghiệm dọc theo trắc dọc sông.

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

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

Kết quả nghiên cứu của luận án cung cấp bức tranh định lượng chính xác và mang tính cảnh báo cao về hiện trạng ô nhiễm và sức chịu tải của lưu vực:

  1. Tình trạng suy thoái oxy và ô nhiễm hữu cơ nghiêm trọng: Trên sông Nhuệ, đoạn từ hạ lưu cống Hà Đông đến Phủ Lý, nồng độ $DO$ thường xuyên dao động ở mức cực thấp ($0,5 - 2,0$ mg/L), nồng độ $BOD_5$ và $COD$ vượt quy chuẩn QCVN 08-MT:2015 cột A2 từ 3 đến 8 lần. Nồng độ $NH_4^+$ có những thời điểm vượt tiêu chuẩn cho phép hàng chục lần, biến sông Nhuệ thành mương thoát nước thải hở mang tính chất dòng chảy kỵ khí vào mùa kiệt.
  2. Sự thâm hụt sức chịu tải ($L_{tn} < 0$): Kết quả tính toán tải lượng ô nhiễm tối đa cho thấy trên hầu hết các phân đoạn của sông Nhuệ (đặc biệt là đoạn từ Liên Mạc qua Hà Đông, Cự Đà đến Cống Thần) và sông Đáy (đoạn tiếp nhận nước sông Nhuệ tại Phủ Lý), khả năng tiếp nhận nước thải đối với các chỉ tiêu $BOD_5$, $COD$ và $NH_4^+$ theo mục đích sử dụng A2 (cấp nước sinh hoạt) và B1 (tưới tiêu) đều mang giá trị âm sâu. Điều này có nghĩa là nguồn nước sông không còn khả năng tự làm sạch để tiếp nhận thêm bất kỳ tải lượng ô nhiễm nào từ các nguồn xả mới.
  3. Tác động khuếch đại tiêu cực của biến đổi khí hậu: Dưới kịch bản biến đổi khí hậu RCP4.5 vào các năm 2020 và 2030, lượng mưa mùa khô trên lưu vực có xu hướng suy giảm trong khi nhiệt độ trung bình tăng, dẫn đến dòng chảy kiệt bình quân tại trạm Ba Thá giảm từ 10% đến 18%. Hiện tượng này làm suy giảm nghiêm trọng lưu lượng pha loãng tự nhiên, khiến sức chịu tải của sông Đáy và sông Nhuệ giảm thêm 15–20% so với điều kiện thủy văn hiện tại.
  4. Tương quan kịch bản xử lý nước thải 80% vs 0%: Nếu đến năm 2020 và 2030, toàn bộ các nguồn thải lớn (đô thị, KCN, làng nghề) được thu gom và xử lý đạt 80% trước khi xả ra sông, chất lượng nước sông Nhuệ – Đáy sẽ có bước cải thiện vượt bậc: nồng độ $DO$ tăng lên trên 4 mg/L, chỉ tiêu $BOD_5$ và $NH_4^+$ giảm từ 50% đến 70%, từng bước khôi phục khả năng chịu tải dương cho mục đích tưới tiêu (B1). Ngược lại, với kịch bản 0% xử lý và phát triển KTXH tự nhiên, toàn bộ hệ thống hạ lưu sẽ rơi vào trạng thái "chết sinh học".

Implications đa chiều

  • Đóng góp lý thuyết: Công trình cung cấp minh chứng thực nghiệm khẳng định giới hạn chịu tải phi tuyến của hệ thống thủy vực nhiệt đới, làm phong phú thêm lý thuyết quản lý tổng hợp tài nguyên nước (IWRM) trong điều kiện bất định của khí hậu.
  • Đổi mới phương pháp luận: Thiết lập một quy trình chuẩn hóa kết nối mô hình thủy văn NAM – mô hình 1D MIKE 11 – đánh giá $TMDL$, có khả năng chuyển giao và áp dụng cho các lưu vực sông liên tỉnh khác tại Việt Nam như lưu vực sông Cầu, sông Mã hay sông Vu Gia – Thu Bồn.
  • Ứng dụng thực tiễn và quản lý: Cung cấp cơ sở dữ liệu số hóa và bảng hạn ngạch xả thải ($WLA$) chi tiết cho từng đoạn sông, giúp Sở Tài nguyên và Môi trường các tỉnh Hà Nội, Hà Nam, Nam Định, Ninh Bình, Hòa Bình kiểm soát việc cấp phép xả thải vào nguồn nước.
  • Khuyến nghị chính sách: Đề xuất cấp bách việc hoàn thiện thể chế của Ủy ban Quản lý Lưu vực sông Nhuệ – Đáy, thiết lập cơ chế điều tiết tiếp nước cưỡng bức từ sông Hồng qua cống Liên Mạc và đập Đáy trong mùa kiệt nhằm duy trì dòng chảy môi trường tối thiểu.

Limitations và Future Research

Tác giả đã chỉ ra một số hạn chế nội tại của đề tài:

  • Do điều kiện mạng lưới quan trắc, mô hình MIKE 11 là mô hình 1 chiều, giả định mặt cắt ngang đồng nhất, chưa mô tả chi tiết được hoàn toàn các dòng xoáy cục bộ và hiện tượng phân tầng chất lượng nước tại các khu vực hợp lưu phức tạp (như ngã ba Phủ Lý, Gián Khẩu).
  • Tải lượng từ các nguồn thải phân tán (nguồn diện từ nông nghiệp trồng trọt, chăn nuôi gia hộ) được ước tính thông qua hệ số phát thải trung bình và diện tích canh tác, chưa đo đạc liên tục theo thời gian thực tại các cống xả nhỏ.
  • Động thái giải phóng chất dinh dưỡng từ lớp bùn đáy trầm tích tích tụ lâu năm chưa được mô hình hóa như một module động học độc lập mà mới dừng lại ở các tham số hằng số tốc độ trong Ecolab.

Chương trình nghiên cứu tiếp theo (Future Research Agenda) cần tập trung vào 4 hướng:

  1. Phát triển mô hình thủy lực – chất lượng nước 2D/3D (như Delft3D hoặc MIKE 21) cho các phân đoạn sông nhạy cảm và vùng cửa sông chịu triều mặn.
  2. Tích hợp mạng lưới cảm biến tự động IoT kết hợp viễn thám vệ tinh quan trắc chất lượng nước theo thời gian thực để cập nhật biên đầu vào liên tục cho mô hình.
  3. Nghiên cứu sâu hơn về cơ chế hấp phụ – giải hấp các vi chất ô nhiễm mới nổi (dược phẩm tồn dư, hóa chất bảo vệ thực vật, hạt vi nhựa) trong trầm tích sông Nhuệ – Đáy.
  4. Đánh giá kinh tế học môi trường toàn diện nhằm xác định chi phí – lợi ích của việc thực thi hạn ngạch xả thải $TMDL$ trên quy mô liên tỉnh.

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

  • Tác động học thuật: Luận án là tài liệu tham khảo chuẩn mực trong đào tạo sau đại học chuyên ngành Khoa học Môi trường và Kỹ thuật Môi trường. Phương pháp luận kết nối chuỗi mô hình của luận án đã được công bố trên các tạp chí chuyên ngành uy tín trong nước và là tiền đề cho nhiều công trình nghiên cứu sâu hơn.
  • Tác động chính sách: Kết quả tính toán sức chịu tải của luận án đóng góp trực tiếp làm luận cứ khoa học phục vụ việc rà soát "Đề án tổng thể bảo vệ môi trường lưu vực sông Nhuệ – sông Đáy đến năm 2020", hỗ trợ Bộ Tài nguyên và Môi trường trong việc xây dựng các văn bản hướng dẫn thi hành Luật Bảo vệ Môi trường 2014 và Luật Tài nguyên nước 2012.
  • Tác động kinh tế – xã hội: Định hướng các dự án đầu tư nhà máy xử lý nước thải tập trung tại Hà Nội và Hà Nam, bảo vệ trực tiếp nguồn nước tưới cho hàng trăm nghìn hecta đất nông nghiệp vùng hạ du và nguồn nước thô phục vụ đời sống dân sinh của gần 12 triệu người dân.

Đố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 mô hình hóa thủy văn – thủy lực – chất lượng nước đồng bộ, cấu trúc bài toán $TMDL$ và quy trình hiệu chuẩn thông số mô hình chuyên sâu.
  • Các nhà khoa học và giảng viên: Cung cấp hệ thống dữ liệu thực nghiệm phong phú về trắc dọc hình học, chế độ thủy văn và diễn biến chất lượng nước sông Nhuệ – Đáy phục vụ giảng dạy và nghiên cứu phát triển mô hình sinh thái.
  • Cơ quan quản lý nhà nước (Bộ TN&MT, Chi cục BVMT lưu vực sông Nhuệ – Đáy, Sở TN&MT các tỉnh): Sử dụng các bảng phân vùng tải lượng tiếp nhận tối đa làm công cụ kỹ thuật phục vụ thẩm định báo cáo ĐTM và cấp giấy phép xả thải vào nguồn nước.
  • Các kỹ sư quy hoạch đô thị và hạ tầng môi trường: Sử dụng dự báo nhu cầu dùng nước và tải lượng ô nhiễm để thiết kế công suất các hệ thống thu gom và xử lý nước thải đô thị, KCN đến năm 2030.

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 khung lý thuyết tải lượng tối đa theo ngày ($TMDL$ của US EPA) và lý thuyết tự làm sạch truyền thống sang điều kiện thủy văn biến đổi phi tuyến do biến đổi khí hậu (kịch bản RCP4.5) trên hệ thống sông đồng bằng chịu ảnh hưởng của chế độ nhân tạo hóa và dòng chảy kiệt cực đoan. Công trình đã lượng hóa thành công mối quan hệ giữa quá trình động học phân hủy sinh hóa với tải lượng ô nhiễm giới hạn trên từng đoạn sông chức năng.

2. Đổi mới phương pháp luận của luận án khi so sánh với ít nhất 2 nghiên cứu đi trước?
So với nghiên cứu của Nguyễn Văn Cư (2005) chỉ mô phỏng chất lượng nước tĩnh và nghiên cứu của Nguyễn Minh Lâm (2012) trên sông Vàm Cỏ Đông chỉ xét kịch bản KTXH thông thường, luận án của Nguyễn Toàn Thắng đã giải quyết triệt để bài toán thiếu số liệu thủy văn bằng mô hình NAM, tích hợp mô hình 1 chiều MIKE 11 HD-AD-Ecolab với 335 mặt cắt ngang toàn lưu vực, và lồng ghép đồng thời kịch bản biến đổi khí hậu RCP4.5 với kịch bản xử lý nước thải 80% tại nguồn.

3. Phát hiện bất ngờ và phản trực giác nhất từ dữ liệu thực nghiệm là gì?
Phát hiện bất ngờ nhất là tại nhiều đoạn trên sông Nhuệ, nồng độ $DO$ tiệm cận mức 0 mg/L kéo dài hàng chục km ngay cả trong điều kiện dòng chảy mùa mưa nếu không có nguồn nước pha loãng bổ cập từ sông Hồng. Điều này chứng minh rằng tải lượng hữu cơ tích tụ trong lòng dẫn và nguồn xả tại chỗ đã vượt qua điểm tới hạn (tipping point), làm vô hiệu hóa hoàn toàn cơ chế tái cấp khí tự nhiên của mặt nước.

4. Luận án có cung cấp quy trình tái lập (replication protocol) minh bạch không?
Luận án cung cấp chi tiết toàn bộ quy trình tái lập: từ danh mục 131 mặt cắt sông Hồng, 43 mặt cắt sông Nhuệ, 96 mặt cắt sông Đáy; bảng thông số hiệu chỉnh mô hình NAM ($U_{max}, L_{max}, CQOF, CQIF$); quy trình lấy mẫu và phân tích theo TCVN 7325:2004, TCVN 6001-1:2008, TCVN 6491:1999; cùng các bước thiết lập điều kiện biên và module Ecolab trong MIKE 11.

5. Luận án đã phác thảo chương trình nghiên cứu 10 năm như thế nào?
Chương trình nghiên cứu dài hạn được định hướng vào việc xây dựng hệ thống trợ giúp ra quyết định (DSS) theo thời gian thực cho toàn bộ lưu vực sông Hồng – Thái Bình, nâng cấp mô hình chất lượng nước 2D/3D tích hợp cơ chế trao đổi bùn đáy – nước mặt, và đánh giá tác động sinh thái học độc chất ($ecotoxicology$) đối với chuỗi thức ăn thủy vực.

Kết luận

Luận án tiến sĩ của Nguyễn Toàn Thắng đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 5 đóng góp cốt lõi:

  1. Hoàn thiện phương pháp luận đánh giá toàn diện khả năng chịu tải và tiếp nhận nước thải của lưu vực sông liên tỉnh trong điều kiện kết hợp giữa phát triển kinh tế – xã hội và biến đổi khí hậu.
  2. Xây dựng thành công chuỗi mô hình tích hợp NAM – MIKE 11 (HD, AD, Ecolab) với độ chính xác cao (chỉ số Nash-Sutcliffe đạt từ 0,78 đến 0,86), khôi phục thành công chuỗi dữ liệu thủy văn và mô phỏng chi tiết 13 thông số chất lượng nước.
  3. Lượng hóa chính xác hiện trạng suy thoái chất lượng nước và chỉ ra tình trạng thâm hụt sức chịu tải ($L_{tn} < 0$) trầm trọng trên hầu hết trục chính sông Nhuệ và hạ lưu sông Đáy.
  4. Dự báo định lượng biến động sức chịu tải đến năm 2020 và 2030 theo kịch bản RCP4.5, chứng minh tính cấp thiết của ngưỡng xử lý tối thiểu 80% nước thải tại nguồn để cứu vãn hệ sinh thái lưu vực.
  5. Đề xuất hệ thống giải pháp công trình và phi công trình đồng bộ, tạo cơ sở khoa học vững chắc phục vụ công tác quy hoạch, quản lý tài nguyên nước và bảo vệ môi trường bền vững cho vùng thủ đô và hạ lưu sông Hồng.