Tổng quan về luận án

Hiện tượng phú dưỡng (eutrophication) tại các thủy vực nước ngọt nội đô đang trở thành thách thức môi trường cấp bách trên phạm vi toàn cầu, đặc biệt tại các quốc gia thuộc vùng khí hậu nhiệt đới gió mùa có tốc độ đô thị hóa nhanh chóng như Việt Nam. Luận án tiến sĩ kỹ thuật chuyên ngành Môi trường Đất và Nước với đề tài "Nghiên cứu phát triển mô hình toán mô phỏng quá trình phú dưỡng ở các vùng nước tĩnh nông, ứng dụng cho hồ Cự Chính - Hà Nội" do nghiên cứu sinh Tạ Đăng Thuần thực hiện dưới sự hướng dẫn khoa học của PGS.TS. Bùi Quốc Lập (Trường Đại học Thủy Lợi, 2019) là công trình nghiên cứu tiên phong trong việc giải quyết bài toán mô phỏng động học sinh địa hóa ở các thủy vực nước tĩnh quy mô nhỏ.

Về bối cảnh khoa học, phần lớn các mô hình phú dưỡng kinh điển trên thế giới được phát triển tập trung cho các hồ tự nhiên có diện tích lưu vực và độ sâu lớn tại các vùng khí hậu ôn đới (Bắc Mỹ, Tây Âu, Bắc Âu). Nền tảng lý thuyết sinh thái học thủy vực khẳng định: "Ở các vùng nước tĩnh thì hồ nông với độ sâu trung bình nhỏ hơn 5 mét là nơi thường xuyên xảy ra hiện tượng phú dưỡng" [1]. Tuy nhiên, các hồ nội thành Hà Nội có các đặc tính thủy động lực học và sinh thái học rất dị biệt: diện tích nhỏ ($< 1\text{ ha}$ đến vài chục hecta), nông ($h_0 \approx 1,5 - 2,5\text{ m}$), chịu tác động phân tầng nhiệt ngày-đêm rõ rệt, không có dòng chảy liên tục và tiếp nhận lượng lớn chất ô nhiễm phân tán. Khoảng trống nghiên cứu (research gap) cốt lõi tồn tại ở ba điểm:

  1. Thiếu các mô hình động thái sinh thái (ecological dynamic models) mô tả chi tiết tương tác đồng thời giữa nhiều chu trình dưỡng chất (N, P, C, O) và sự cạnh tranh sinh thái giữa các nhóm vi tảo nhiệt đới.
  2. Sự phụ thuộc vào các phần mềm thương mại ngoại nhập ("mô hình hộp đen" đòi hỏi số lượng dữ liệu biên khổng lồ, chi phí bản quyền cao, không linh hoạt trong việc tùy biến phương trình động học).
  3. Bỏ qua nguồn dinh dưỡng ngoại sinh từ lắng đọng khí quyển (khuếch tán khô, nước mưa ướt) và dòng chảy tràn bề mặt đô thị trong các phương trình vi phân động học.

Nhằm giải quyết triệt để các khoảng trống trên, luận án thiết lập 3 câu hỏi nghiên cứu (Research Questions - RQ) và 3 giả thuyết khoa học (Hypotheses - H) tương ứng:

  • RQ1: Cấu trúc hệ phương trình vi phân mô tả quá trình biến đổi sinh học - hóa học nào phản ánh chính xác nhất động học phú dưỡng trong điều kiện hồ nước tĩnh nông nhiệt đới?
    • H1: Một hệ phương trình vi phân tích hợp 14 biến trạng thái với các hàm giới hạn phi tuyến tính đa yếu tố (dinh dưỡng, nhiệt độ, bức xạ) sẽ mô tả được sự bùng phát của từng nhóm tảo chuyên biệt.
  • RQ2: Nguồn dinh dưỡng bổ sung từ khí quyển và dòng chảy tràn đóng vai trò định lượng như thế nào đối với động thái phú dưỡng hồ nội đô khi không có nguồn thải điểm trực tiếp?
    • H2: Lắng đọng khí quyển và nước mưa chảy tràn là tác nhân kích hoạt chính làm thay đổi tỷ lệ $TN/TP$, thúc đẩy hiện tượng nở hoa của vi khuẩn lam vào mùa hè.
  • RQ3: Phương pháp tối ưu hóa tham số nào khắc phục được tính bất định (equifinality) của hệ thống đa tham số sinh thái phức tạp trong điều kiện chuỗi dữ liệu thực nghiệm hạn chế?
    • H3: Ứng dụng thuật giải di truyền (Genetic Algorithm - GA) kết hợp phân tích độ nhạy toàn cục sẽ tối ưu hóa tự động bộ tham số động học, nâng cao độ chính xác dự báo so với phương pháp thử sai truyền thống.

Khung lý thuyết của nghiên cứu được xây dựng trên sự tích hợp giữa Định luật tối thiểu Liebig (Liebig's Law of the Minimum), Động học hấp thu enzyme Michaelis-Menten, Mô hình ức chế ánh sáng Steele, Động học phụ thuộc nhiệt độ Arrhenius/Gaussian và Mô hình dinh dưỡng chuỗi thức ăn Lotka-Volterra. Phạm vi nghiên cứu tập trung vào hồ Cự Chính (quận Thanh Xuân, Hà Nội) với diện tích mặt nước $A = 4.000\text{ m}^2$, độ sâu trung bình $h_0 = 1,65\text{ m}$, chuỗi số liệu quan trắc hiệu chỉnh kéo dài 140 ngày trong năm 2017 và kiểm định độc lập trong năm 2018.

Literature Review và Positioning

Lịch sử phát triển mô hình phú dưỡng hồ trải qua ba dòng nghiên cứu chính (research streams):

  1. Dòng mô hình thực nghiệm và cân bằng khối lượng đơn giản: Khởi xướng bởi Vollenweider (1968, 1976), Dillon & Rigler (1974), Larsen & Mercier (1976), Lorenzen (1979), Chapra & Canale (1991). Các mô hình này giả định hồ là một bể khuấy trộn hoàn chỉnh ở trạng thái dừng, mặc định phốt pho tổng số (TP) là nhân tố hạn chế duy nhất và sử dụng đồ thị tải lượng phốt pho để dự báo phản ứng phục hồi dài hạn của hồ. Nhược điểm chí tử là không thể dự báo động thái ngắn hạn (vài ngày đến vài tuần) và bỏ qua tương tác sinh hóa chi tiết.
  2. Dòng mô hình sinh thái động học đa biến: Phát triển bởi Di Toro et al. (1971), Thomann et al. (1975), Jørgensen (1976), Cerco & Cole (1993 - vịnh Chesapeake), James (2005 - hồ Okeechobee), Cao et al. (2008) và Recknagel et al. (2008 - mô hình SALMO), Zhang et al. (2008 - hồ Erie). Dòng mô hình này phân rã chu trình dinh dưỡng thành các cấu phần hòa tan/hạt, hữu cơ/vô cơ, đồng thời tích hợp các nhóm tảo cạnh tranh và động vật phù du (zooplankton).
  3. Dòng chỉ số trạng thái dinh dưỡng (TSI): Carlson (1977) thiết lập chỉ số $TSI(TP)$ và $TSI(Chla)$, sau đó Kratzer & Brezonik (1981) bổ sung $TSI(TN)$ cho các thủy vực bị hạn chế bởi nitơ, và Håkanson et al. (2007) chuẩn hóa bảng phân loại phú dưỡng theo $TN, TP, Chla$.
                                 TIẾN TRÌNH TIẾP CẬN LÝ THUYẾT VÀ ĐỊNH VỊ NGHIÊN CỨU
                                 
  [1. Mô hình Cân bằng Khối lượng]     [2. Chỉ số Trạng thái Dinh dưỡng]     [3. Mô hình Thủy động lực 3D]
  - Vollenweider (1968), Lorenzen (1979) - Carlson (1977), Håkanson (2007)     - EFDC, CAEDYM, CE-QUAL-W2
  - Đơn biến (TP), trạng thái dừng       - Phân loại tĩnh qua TN, TP, Chla     - Phức tạp, đòi hỏi dữ liệu lớn
                         [Khoảng trống nghiên cứu cho Hồ Nông Nội Đô Nhiệt Đới]
                         - Thiếu mô hình động học mở tích hợp ngoại sinh (mưa, chảy tràn)
                         - Thiếu cơ chế phân rã 3 nhóm vi tảo cạnh tranh nhiệt đới
                         - Thiếu công cụ tự tối ưu tham số không gian nghiệm lớn (GA)
                                    [MÔ HÌNH PHÁT TRIỂN CỦA LUẬN ÁN]
                         - Hệ 14 PTVP thường (ODE-45), 88 tham số động học
                         - Tích hợp hàm ngoại sinh mưa, lắng đọng khí quyển, bức xạ, nhiệt độ
                         - Tối ưu hóa 15 tham số nhạy nhất bằng Genetic Algorithm (Matlab)

Trong bức tranh học thuật, hai cuộc tranh luận lý thuyết lớn đang diễn ra:

  • Tranh luận 1: Đơn hạn chế Phốt pho (P-limitation) đối lập với Đồng hạn chế Nitơ - Phốt pho (N & P co-limitation): Trường phái Vollenweider - Schindler cho rằng kiểm soát duy nhất phốt pho là đủ để chặn đứng phú dưỡng. Ngược lại, Paerl et al. (2011) và Lewis & Wurtsbaugh (2008) chứng minh rằng ở các vùng nước tĩnh nông nhiệt đới, tỷ lệ $TN/TP$ biến động mạnh, nitơ thường xuyên trở thành yếu tố hạn chế sinh trưởng trong mùa khô và chuyển sang đồng hạn chế hoặc hạn chế phốt pho trong mùa mưa.
  • Tranh luận 2: Mô hình phân giải không gian 3D phức tạp đối lập với Mô hình sinh thái tập trung (0D/1D lumped dynamic models): Việc ứng dụng các mô hình 3D như EFDC (được US EPA tài trợ) vào hồ Tây (Nguyễn Xuân Tình et al.) hay hồ Erie (Zhang et al., 2008) tuy giải quyết được trường vận tốc và phân tầng không gian nhưng đòi hỏi khối lượng dữ liệu địa hình biên và công suất tính toán rất lớn, dễ rơi vào bế tắc khi áp dụng cho các hồ đô thị vừa và nhỏ.

Tại Việt Nam, nghiên cứu của Nguyễn Thanh Hồng (2010) bước đầu phát triển mô hình cân bằng khối lượng TP cho hồ chứa sâu (hồ Tây, hồ Đại Lải, hồ Tuyên Quang) bằng phương pháp giải số Runge-Kutta bậc 4, nhưng hoàn toàn bỏ qua chu trình nitơ, carbon, nồng độ oxy hòa tan (DO) và sự phân chia nhóm vi tảo.

So sánh với hai nghiên cứu quốc tế điển hình:

  • Mô hình hồ Erie của Zhang et al. (2008): Kết hợp thủy động lực học 2D/3D với chuỗi thức ăn bao gồm động vật đáy (trai vằn Dreissena polymorpha). Mô hình này tập trung vào tác động lọc của sinh vật đáy ở hồ lớn ôn đới, trong khi luận án của Tạ Đăng Thuần tinh chỉnh cấu trúc hướng vào động thái hồ nông nhiệt đới, nơi quá trình phân hủy kỵ khí lớp bùn đáy và hấp phụ/giải phóng dưỡng chất bề mặt trầm tích chiếm ưu thế.
  • Mô hình hồ Okeechobee của James (2005): Phân chia cột nước thành 3 tầng và trầm tích 2 lớp (hiếu khí/yếm khí). Mô hình của luận án kế thừa phương trình trao đổi dinh dưỡng bề mặt đáy nhưng cải tiến vượt bậc bằng cách đưa trực tiếp thông lượng dưỡng chất từ nước mưa chảy tràn đô thị vào hệ phương trình vi phân thường (ODEs).

Đó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à đóng góp sâu sắc cho lý thuyết mô hình hóa sinh thái học thủy vực (Aquatic Ecological Modeling Theory) thông qua các điểm đột phá:

  1. Mở rộng phương trình động học dưỡng chất vô cơ: Bổ sung thành phần tải lượng ngoại sinh phân tán tức thời từ nước mưa và dòng chảy tràn bề mặt đô thị trực tiếp vào phương trình biến thiên của Carbon hữu cơ dạng hạt ($POC$), Phốt phát vô cơ hòa tan ($DIP$), Amoni ($NH_4\text{-}N$), Nitrit ($NO_2\text{-}N$) và Nitrat ($NO_3\text{-}N$).
  2. Cụ thể hóa lý thuyết cạnh tranh sinh thái của vi tảo: Phân chia thực vật nổi ($TVN$) thành 3 nhóm taxon có đặc tính sinh lý riêng biệt: Tảo lục (Chlorophyta - $y_1$), Vi khuẩn lam (Cyanobacteria - $y_2$) và Tảo silic (Bacillariophyta/Diatom - $y_3$). Tốc độ sinh trưởng của từng nhóm được mô tả bằng tích của ba hàm giới hạn phi tuyến độc lập: $$\mu_k = \alpha_{1,k} \cdot f_1(N, P)_k \cdot f_2(I)_k \cdot f_3(T)_k$$

Trong đó:

  • Hàm giới hạn dinh dưỡng $f_1(N,P)k = \min\left(\frac{C{DIN}}{K_{HN,k} + C_{DIN}}, \frac{C_{DIP}}{K_{HP,k} + C_{DIP}}\right)$ tuân thủ nghiêm ngặt Định luật Liebig.
  • Hàm giới hạn ánh sáng $f_2(I)k$ tích hợp độ suy giảm ánh sáng qua cột nước $h_0$ theo công thức thực nghiệm cải tiến của hệ số tắt quang học tích hợp nồng độ diệp lục: $$K{ess} = 0,825 + 0,043\sum(Chla_C_k \cdot y_k) + 0,299\left[\sum(Chla_C_k \cdot y_k)\right]^{2/3}$$
  • Hàm giới hạn nhiệt độ $f_3(T)k$ sử dụng dạng hàm mũ kép đối xứng quanh nhiệt độ tối ưu $T{opt,k}$: $$f_3(T)k = \begin{cases} \exp\left[-K{1,k}(T - T_{k})^2\right] & \text{khi } T < T_{k} \ \exp\left[-K_{2,k}(T - T_{k})^2\right] & \text{khi } T \ge T_{k} \end{cases}$$

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

Mô hình toán học phát triển gồm hệ 14 phương trình vi phân thường (Ordinary Differential Equations - ODEs) phi tuyến liên kết chặt chẽ:

  • Nhóm sinh học (4 biến): Sinh khối Tảo lục ($y_1$), Vi khuẩn lam ($y_2$), Tảo silic ($y_3$), Động vật phù du ($y_4$) [đơn vị: $\text{mg C/L}$].
  • Nhóm carbon hữu cơ (3 biến): Tổng carbon hữu cơ ($TOC - y_5$), Carbon hữu cơ dạng hạt ($POC - y_9$), Carbon hữu cơ hòa tan ($DOC - y_{10}$) [$\text{mg C/L}$].
  • Nhóm dưỡng chất và khí hòa tan (7 biến): Tổng phốt pho ($TP - y_6$), Tổng nitơ ($TN - y_7$), Oxy hòa tan ($DO - y_8$), Phốt phát vô cơ hòa tan ($DIP/PO_4\text{-}P - y_{11}$), Amoni ($NH_4\text{-}N - y_{12}$), Nitrit ($NO_2\text{-}N - y_{13}$), Nitrat ($NO_3\text{-}N - y_{14}$) [$\text{mg/L}$].
                                  SƠ ĐỒ KHUNG ĐỘNG HỌC MÔ HÌNH PHÚ DƯỠNG (14 BIẾN)
                                  
    [Ngoại sinh: Mưa (X), Nhiệt độ (T), Bức xạ (I), Lắng đọng khí quyển, Chảy tràn đô thị]

Điểm độc đáo trong cấu trúc phương trình thể hiện ở việc đưa trực tiếp xung lượng mưa $X$ ($\text{mm}$) và tỷ lệ diện tích lưu vực thu gom $A_{catch} / A_{lake}$ vào các số hạng nguồn: $$\frac{dy_{11}}{dt} = \sum_{k=1}^3 APC_k(B_{3,k} - B_{1,k}) - APC_{ZP} B_4 + APC_{POC} B_{10} + APC_{DOC} B_{13} + B_{14} + \left[0,15 \cdot \frac{X}{1.000} \cdot 1.500\right]$$ $$\frac{dy_{12}}{dt} = -(1 - r_1)\sum_{k=1}^3 ANC_k B_{1,k} + \sum_{k=1}^3 ANC_k B_{3,k} - ANC_{ZP} B_4 + ANC_{POC} B_{10} + ANC_{DOC} B_{13} + B_{15} - B_{16} + \left[0,75 \cdot \frac{X}{1.000} \cdot 1.500\right]$$

Điều kiện biên và phạm vi áp dụng của mô hình được xác định tường minh: áp dụng cho các vùng nước tĩnh nông ngọt ($h_0 < 5\text{ m}$), khuấy trộn đều, không có phân tầng nhiệt vĩnh cửu, không có dòng chảy xuất nhập liên tục dạng kênh hở, và chịu tác động áp đảo của các nguồn ô nhiễm phân tán đô thị.

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

Thiết kế nghiên cứu

Nghiên cứu đứng trên lập trường nhận thức luận thực chứng - cơ chế (Positivist-Mechanistic Paradigm), kết hợp chặt chẽ giữa điều tra thực địa quy chuẩn, phân tích hóa lý nghiệm sinh trong phòng thí nghiệm và mô hình hóa toán học số trị (Numerical Mathematical Modeling).

Thiết kế nghiên cứu đa tầng bao gồm:

  • Tầng động lực khí quyển/ngoại cảnh: Đo đạc bức xạ mặt trời ($I_0$), nhiệt độ không khí, lượng mưa ngày ($X$) tại Trạm Khí tượng Láng (Hà Nội).
  • Tầng động lực cột nước: Quan trắc các thông số lý - hóa - sinh tại hiện trường và phòng thí nghiệm.
  • Tầng động học trầm tích: Mô tả thông lượng giải phóng phốt pho ($B_{14}$), amoni ($B_{15}$) và nhu cầu oxy đáy ($SOD - B_{19}$) phụ thuộc nhiệt độ: $$B_{14} = \frac{\alpha_{14} \cdot \exp(\beta_3 \cdot T)}{h_0}; \quad B_{19} = \frac{\alpha_{19} \cdot \exp(\beta_5 \cdot T)}{h_0} \cdot 10^{-3}$$

Quy trình nghiên cứu rigorous

Quy trình thu thập và phân tích mẫu tuân thủ nghiêm ngặt các Tiêu chuẩn Việt Nam (TCVN) và Standard Methods for the Examination of Water and Wastewater (APHA):

  • Thu mẫu thực địa: Lấy mẫu tích hợp tầng mặt và tầng đáy tại hồ Cự Chính trong 9 đợt quan trắc năm 2017 (các ngày thứ 1, 15, 28, 47, 78, 96, 112, 127, 140) đại diện cho các chu kỳ thời tiết chuyển mùa.
  • Phân tích trong phòng thí nghiệm:
    • $DO, pH, Nhiệt độ, Độ dẫn điện (EC)$ đo trực tiếp tại hiện trường bằng máy đo đa chỉ tiêu cầm tay HACH HQ40D.
    • $BOD_5$ đo bằng phương pháp chai đo ủ 5 ngày ở $20^\circ\text{C}$.
    • $TN, TP$ xác định bằng phương pháp phá mẫu Persulfate kết hợp so màu quang phổ UV-Vis.
    • $NH_4\text{-}N, NO_2\text{-}N, NO_3\text{-}N, PO_4\text{-}P$ phân tích bằng phương pháp trắc quang theo TCVN.
    • Định lượng tế bào thực vật nổi: Cố định mẫu bằng dung dịch Lugol 1%, sử dụng buồng đếm Sedgewick-Rafter dưới kính hiển vi quang học Olympus BX53 ở độ phóng đại $200\text{X} - 400\text{X}$, quy đổi thể tích sinh khối tế bào sang đơn vị Carbon ($\text{mg C/L}$) theo công thức hình học chuẩn của Hillebrand et al. (1999).
  • Phân tích chất lượng nước mưa: Dữ liệu nồng độ dinh dưỡng trong nước mưa tại trạm Cự Khê và nội thành Hà Nội được chuẩn hóa: nồng độ trung bình $NO_3\text{-}N$ đạt $0,96\text{ mg/L}$ ($0,1 - 8,6\text{ mg/L}$), $NH_4\text{-}N$ đạt $0,33\text{ mg/L}$ ($0,03 - 0,86\text{ mg/L}$), $NO_2\text{-}N$ đạt $0,23\text{ mg/L}$ ($0 - 1,398\text{ mg/L}$), và $TN$ nước mưa trung bình $0,9\text{ mg/L}$, $TP$ trung bình $0,015\text{ mg/L}$ theo công bố của Minju Lee et al. (2017) [17] và FISRWG [16].

Data và phân tích

Hệ phương trình vi phân phi tuyến bậc 1 được giải số học bằng thuật toán Runge-Kutta bậc 4-5 tự thích ứng bước thời gian thông qua hàm solver ode45 trên nền tảng phần mềm MATLAB R2018a.

Tổng số tham số của mô hình là 88 tham số. Trong đó, 35 tham số là hằng số sinh hóa đã được chuẩn hóa quốc tế (tỷ lệ phân tử $APC, ANC, ADC$, hệ số nhiệt động $\theta, \beta$). Đối với 53 tham số biến thiên còn lại, tác giả thực hiện phân tích độ nhạy (Sensitivity Analysis - SA) toàn cục. Kết quả sàng lọc được 15 tham số có độ nhạy chi phối cao nhất đối với các biến trạng thái chính (như bảng tổng hợp sau):

Ký hiệu tham số Diễn giải tham số sinh địa hóa Đơn vị Khoảng tìm kiếm Giá trị tối ưu GA
$\alpha_{1,1}$ Tốc độ sinh trưởng cực đại của Tảo lục ở $20^\circ\text{C}$ $\text{ngày}^{-1}$ $1,5 - 3,5$ $2,400$
$I_{k,1}$ Cường độ ánh sáng bão hòa của Tảo lục $\text{W/m}^2$ $150 - 250$ $185,5$
$T_{k,1}$ Nhiệt độ nước tối ưu cho Tảo lục sinh trưởng $^\circ\text{C}$ $26,0 - 32,0$ $28,5$
$K_{HN,1}$ Hằng số nửa bão hòa hấp thụ Nitơ của Tảo lục $\text{mg/L}$ $0,01 - 0,09$ $0,045$
$K_{HP,1}$ Hằng số nửa bão hòa hấp thụ Phốt pho của Tảo lục $\text{mg/L}$ $0,001 - 0,06$ $0,018$
$\alpha_{1,2}$ Tốc độ sinh trưởng cực đại của Vi khuẩn lam $\text{ngày}^{-1}$ $1,8 - 4,0$ $3,200$
$I_{k,2}$ Cường độ ánh sáng bão hòa của Vi khuẩn lam $\text{W/m}^2$ $80 - 140$ $98,2$
$T_{k,2}$ Nhiệt độ nước tối ưu cho Vi khuẩn lam $^\circ\text{C}$ $30,0 - 36,0$ $33,5$
$K_{HN,2}$ Hằng số nửa bão hòa hấp thụ Nitơ của Vi khuẩn lam $\text{mg/L}$ $0,01 - 0,08$ $0,035$
$K_{HP,2}$ Hằng số nửa bão hòa hấp thụ Phốt pho của Vi khuẩn lam $\text{mg/L}$ $0,001 - 0,04$ $0,012$
$\alpha_{1,3}$ Tốc độ sinh trưởng cực đại của Tảo silic $\text{ngày}^{-1}$ $1,0 - 2,8$ $1,750$
$T_{k,3}$ Nhiệt độ nước tối ưu cho Tảo silic $^\circ\text{C}$ $18,0 - 24,0$ $21,0$
$\alpha_{19}$ Tốc độ tiêu thụ oxy đáy ($SOD$) ở $20^\circ\text{C}$ $\text{g/m}^2\cdot\text{ngày}$ $0,5 - 4,0$ $2,250$
$K_a$ Hệ số tái sục khí khuếch tán oxy bề mặt $\text{ngày}^{-1}$ $0,1 - 1,5$ $0,469$
$\alpha_{4}$ Tốc độ tiêu thụ tảo của Động vật phù du $\text{L/mg C}\cdot\text{ngày}$ $0,1 - 0,8$ $0,310$

Quy trình tối ưu hóa bộ tham số được thực thi tự động thông qua công cụ Genetic Algorithm (GA) trong MATLAB Global Optimization Toolbox. Hàm mục tiêu (Fitness function) được thiết lập là giá trị trung bình của sai số chuẩn quân phương ($RMSE$) trên toàn bộ 14 biến trạng thái: $$\text{Fitness} = \frac{1}{14} \sum_{i=1}^{14} RMSE_i = \frac{1}{14} \sum_{i=1}^{14} \sqrt{\frac{1}{M} \sum_{m=1}^M \left[y_{ob}(i, m) - y_{sim}(i, m)\right]^2}$$

Thuật toán GA được vận hành qua 38 thế hệ lặp độc lập nhằm triệt tiêu bẫy cực trị địa phương (local minima). Bốn chỉ số thống kê tiêu chuẩn quốc tế được áp dụng nghiêm ngặt để đánh giá độ tin cậy mô hình: Sai số quân phương chuẩn hóa ($RSR$), Chỉ số hiệu quả Nash-Sutcliffe ($NSE$), Hệ số phần trăm sai số ($PBIAS$) và Hệ số tương quan Pearson/Spearman ($r$).

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

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

  1. Thực trạng siêu phú dưỡng (Hypereutrophic) nghiêm trọng tại hồ Cự Chính: Kết quả quan trắc chỉ ra nồng độ các chất dinh dưỡng luôn ở mức báo động. Tổng nitơ ($TN$) dao động từ $1,578\text{ mg/L}$ đến $5,400\text{ mg/L}$, Tổng phốt pho ($TP$) dao động từ $0,141\text{ mg/L}$ đến $0,567\text{ mg/L}$, $Chla$ trung bình đạt $58,3 - 74,3\text{ }\mu\text{g/L}$. Theo thang phân loại Håkanson et al. (2007) ($TP > 0,09\text{ mg/L}$ và $Chla > 25\text{ }\mu\text{g/L}$) và chỉ số Carlson ($TSI(TP) = 75,5 - 95,6 > 70$), hồ Cự Chính liên tục nằm trong trạng thái siêu phú dưỡng.
  2. Quy luật biến dịch chất dinh dưỡng giới hạn sinh trưởng theo mùa: Phân tích tỷ lệ phân tử $TN/TP$ làm sáng tỏ sự chuyển pha sinh thái độc đáo. Vào mùa khô và đầu mùa hè, giá trị $TN/TP < 4,5$, Nitơ đóng vai trò là yếu tố giới hạn sinh trưởng của vi tảo. Khi các trận mưa lớn xuất hiện, do nồng độ nitơ trong nước mưa và nước chảy tràn rất cao ($TN/TP$ nước mưa đạt $\approx 60$), tỷ lệ trong hồ tăng vọt lên $TN/TP > 6,0$, chuyển dịch hệ sinh thái sang trạng thái hạn chế bởi Phốt pho.
  3. Cơ chế nở hoa và ưu thế sinh thái của Vi khuẩn lam: Quần xã thực vật nổi ghi nhận 4 nhóm chính gồm Tảo lục (Chlorophyta), Vi khuẩn lam (Cyanobacteria), Tảo silic (Bacillariophyta) và Tảo mắt (Euglenophyta). Trong các tháng nhiệt độ cao ($T_{nước} > 30^\circ\text{C}$, bức xạ $I_0 > 200\text{ W/m}^2$), Vi khuẩn lam (đặc biệt là hai chi độc tố Microcystis và Anabaena) chiếm ưu thế tuyệt đối với mật độ đạt đỉnh $1,269\text{ mg C/L}$ (chiếm $> 65%$ tổng sinh khối thực vật nổi). Mô hình toán chứng minh điều này xuất phát từ hằng số bão hòa ánh sáng thấp ($I_{k,2} = 98,2\text{ W/m}^2$) và nhiệt độ tối ưu cao ($T_{k,2} = 33,5^\circ\text{C}$).
  4. Tác động chi phối của tải lượng lắng đọng khí quyển và chảy tràn: Dữ liệu mô phỏng chứng minh rằng các đợt mưa rào mùa hạ không làm loãng hồ mà ngược lại, bổ sung tức thời một lượng lớn $NO_3\text{-}N$ và $DIP$, tạo cú hích dưỡng chất (nutrient pulse) kích hoạt tảo bùng phát nở hoa ngay sau khi mưa kết thúc từ 2 đến 5 ngày.
  5. Độ chính xác vượt trội của mô hình hiệu chỉnh bằng GA: Quá trình kiểm định mô hình với bộ dữ liệu độc lập năm 2018 cho thấy sự trùng khớp rất cao giữa đường mô phỏng và số liệu thực đo: chỉ số $NSE$ của Tảo lục đạt $0,78$, Vi khuẩn lam đạt $0,84$, Tảo silic đạt $0,72$, $DO$ đạt $0,81$, $PBIAS$ của tất cả các biến trạng thái đều dao động trong biên độ xuất sắc $\pm 15%$, $RSR < 0,50$.
                              SO SÁNH CÁC KỊCH BẢN QUẢN LÝ VÀ PHỤC HỒI HỒ CỰ CHÍNH
                              
  Chỉ tiêu Mô phỏng                Kịch bản Nền        KB1: Tách dòng chảy tràn   KB2: Nạo vét bùn đáy    KB3: Tích hợp Đa giải pháp
  Sinh khối Vi khuẩn lam (mg C/L)      1,24 - 1,27               0,85 (-32%)              0,68 (-45%)             0,32 (-74%)
  Oxy hòa tan DO tầng đáy (mg/L)       1,50 - 3,20               3,80 (+19%)              5,10 (+59%)             6,40 (+100%)
  Nồng độ DIP cột nước (mg/L)          0,15 - 0,57               0,11 (-31%)              0,06 (-62%)             0,03 (-80%)
  Chỉ số Trạng thái Dinh dưỡng (TSI)   75,5 (Siêu phú dưỡng)     64,2 (Phú dưỡng)         56,8 (Phú dưỡng vừa)    46,5 (Trung dưỡng)
  Khả năng xuất hiện Tảo nở hoa        Rất cao (thường xuyên)    Trung bình               Thấp                    Được kiểm soát triệt để

Implications đa chiều

  • Về mặt học thuật và lý thuyết: Công trình đã phát triển hoàn chỉnh một mô hình sinh thái thủy vực mở (Open Ecological Model) giải quyết thành công bài toán tương tác phi tuyến đa dưỡng chất cho vùng nước tĩnh nông nhiệt đới, bổ sung một nhánh lý thuyết quan trọng cho khoa học limnology khu vực Đông Nam Á.
  • Về mặt phương pháp luận: Thiết lập quy trình chuẩn hóa kết hợp phương pháp mô phỏng số giải phương trình vi phân thường (ODEs) với thuật giải tối ưu hóa thông minh nhân tạo (Genetic Algorithm), cung cấp giải pháp vượt qua rào cản thiếu hụt chuỗi dữ liệu lớn trong quan trắc môi trường.
  • Về ứng dụng thực tiễn và chính sách: Kết quả mô phỏng các kịch bản kiểm soát chất lượng nước cung cấp cơ sở định lượng trực tiếp cho Sở Tài nguyên và Môi trường Hà Nội cùng Công ty TNHH MTV Thoát nước Hà Nội:
    1. Giải pháp kỹ thuật công trình: Xây dựng hệ thống rãnh thu gom tách riêng nước mưa đợt đầu quanh hồ để ngăn chặn tải lượng dinh dưỡng bùng phát.
    2. Kiểm soát nội tại: Nạo vét định kỳ lớp trầm tích giàu hữu cơ trước tháng 4 hàng năm nhằm triệt tiêu quá trình giải phóng phốt pho nội tại ($B_{14}$) khi nhiệt độ nước bước vào chu kỳ tăng cao mùa hè.
    3. Hỗ trợ sinh thái: Lắp đặt hệ thống bè thủy sinh nhân tạo kết hợp sục khí cưỡng bức cục bộ tầng đáy nhằm duy trì $DO > 4\text{ mg/L}$, ngăn chặn quá trình khử kỵ khí giải phóng khí độc $H_2S$ và photphin.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn chỉ ra 4 hạn chế khoa học xuất phát từ các giả định đơn giản hóa:

  1. Giả thiết phân bố không gian 0 chiều (0D - Lumped model): Mô hình giả định nồng độ các chất và sinh khối tảo khuấy trộn đồng nhất theo toàn bộ thể tích hồ. Dù phù hợp với hồ Cự Chính ($A = 4.000\text{ m}^2$), giả định này sẽ bộc lộ sai số nếu áp dụng cho các hồ có dạng hình học phức tạp, eo ngách hoặc có diện tích lớn hàng trăm hecta như hồ Tây.
  2. Đơn giản hóa cấu trúc quần xã bậc cao: Thành phần động vật phù du ($y_4$) mới được mô tả như một biến gộp duy nhất, chưa phân tách thành các nhóm ăn lọc kích thước lớn (Cladocera như Daphnia) và nhóm ăn lọc nhỏ (Rotifera, Copepoda), do đó chưa mô tả trọn vẹn hiện tượng kiểm soát từ trên xuống (top-down control).
  3. Mô phỏng động thái trầm tích đáy dạng thông số tĩnh phụ thuộc nhiệt độ: Các quá trình trao đổi phốt pho và tiêu thụ oxy đáy được tham số hóa thông qua các hàm mũ nhiệt độ bán thực nghiệm ($B_{14}, B_{19}$) thay vì mô hình hóa tường minh các lớp phản ứng sinh địa hóa (Diagenetic sediment model) trong cột bùn.
  4. Thiếu vắng module động thái độc tố tảo: Chưa tích hợp phương trình biến thiên nồng độ độc tố vi khuẩn lam (Microcystin-LR, Saxitoxin) hòa tan phát tán khi tế bào tảo bị phân hủy.

Định hướng nghiên cứu tiếp theo trong 5-10 năm tới bao gồm:

  • Mở rộng liên kết mô hình sinh thái 14 biến với module thủy động lực học 2 chiều nông (2D Saint-Venant/Shallow Water Equations) để mô phỏng sự lan truyền không gian của vệt tảo nở hoa dưới tác dụng của trường gió đô thị.
  • Tích hợp động thái độc tố sinh học và chu trình biến đổi silic vô cơ hòa tan ($DSi$) nhằm nâng cao độ chính xác dự báo động thái cạnh tranh giữa tảo silic và vi khuẩn lam.
  • Ứng dụng kỹ thuật mạng thần kinh nhân tạo tích hợp vật lý (Physics-Informed Neural Networks - PINNs) kết hợp thuật giải di truyền để tự động hóa quá trình đồng hóa dữ liệu quan trắc viễn thám vệ tinh (Sentinel-2, Landsat-8) vào mô hình dự báo thời gian thực.

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

Công trình luận án mang lại những giá trị tác động sâu rộng trên nhiều bình diện:

  • Tác động học thuật (Academic Impact): Cung cấp một bộ mã nguồn tính toán mở bằng MATLAB với đầy đủ cấu trúc giải thuật, hệ số động học đã được hiệu chỉnh và kiểm định chi tiết. Đây là tài liệu tham khảo học thuật cốt lõi cho các chương trình đào tạo tiến sĩ, thạc sĩ chuyên ngành Kỹ thuật Môi trường, Thủy văn và Quản lý Tài nguyên Nước.
  • Chuyển đổi công nghệ và dịch vụ kỹ thuật môi trường: Tạo ra công cụ phần mềm nội địa có độ tin cậy tương đương các module sinh thái của MIKE 11/21 hay WASP nhưng với chi phí triển khai bằng 0, giao diện tùy biến linh hoạt, mở ra hướng ứng dụng rộng rãi cho các viện nghiên cứu, công ty tư vấn môi trường tại Việt Nam.
  • Hiệu quả kinh tế - xã hội: Giúp các cơ quan quản lý đô thị tối ưu hóa ngân sách xử lý ô nhiễm hồ. Thay vì áp dụng tràn lan các biện pháp hóa chất đắt đỏ (như rải hóa chất diệt tảo sunfat đồng $CuSO_4$ gây độc hệ sinh thái), mô hình cho phép dự báo chính xác thời điểm bùng phát để can thiệp sớm bằng các giải pháp sinh thái chi phí thấp, bảo vệ cảnh quan và sức khỏe cộng đồng dân cư nội đô.

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

  • Nghiên cứu sinh và Giới học thuật (Doctoral & Senior Researchers): Tiếp cận phương pháp luận mô hình hóa sinh thái toán học hoàn chỉnh, kế thừa bộ 15 tham số nhạy đã được chuẩn hóa cho điều kiện nhiệt đới, mở rộng nghiên cứu sang các hệ thống hồ chứa và lưu vực sông.
  • Kỹ sư R&D Môi trường và Doanh nghiệp Thoát nước: Sử dụng bộ mã nguồn Matlab làm lõi thuật toán để xây dựng hệ thống phần mềm cảnh báo sớm hiện tượng tảo nở hoa và suy giảm oxy hòa tan trên nền tảng WebGIS hoặc ứng dụng di động.
  • Cơ quan Hoạch định Chính sách và Quản lý Nhà nước: Sở Tài nguyên và Môi trường, Chi cục Bảo vệ Môi trường tại các đô thị loại I có cơ sở khoa học xác đáng để ban hành quy chuẩn kỹ thuật địa phương về kiểm soát nước mưa chảy tràn đô thị và tiêu chuẩn vận hành nạo vét hồ chứa.

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?
Trả lời: Đóng góp độc đáo nhất là việc mở rộng Hệ phương trình động học sinh địa hóa kinh điển (mô hình dòng Di Toro - Jørgensen) bằng cách tích hợp trực tiếp xung lượng ô nhiễm phân tán khí quyển và nước mưa chảy tràn đô thị vào hệ phương trình vi phân biến thiên tức thời của 14 biến trạng thái, phá vỡ giả định hệ kín của các mô hình nước tĩnh truyền thống.

2. Điểm cải tiến phương pháp luận cốt lõi so với các nghiên cứu trước đây là gì?
Trả lời: Khắc phục triệt để phương pháp thử sai thủ công truyền thống bằng việc thiết lập khung tích hợp: Phân tích độ nhạy toàn cục để rút gọn không gian 53 tham số xuống 15 tham số then chốt, sau đó áp dụng Thuật giải di truyền (Genetic Algorithm) với 38 thế hệ lặp để tự động tìm kiếm nghiệm tối ưu toàn cục, nâng chỉ số $NSE$ từ mức không đạt ($< 0,4$) lên mức xuất sắc ($> 0,75 - 0,84$).

3. Phát hiện bất ngờ nhất từ chuỗi số liệu thực nghiệm là gì?
Trả lời: Trái ngược với quan niệm phổ biến cho rằng các trận mưa rào mùa hạ sẽ làm loãng và cải thiện chất lượng nước hồ, mô hình và dữ liệu chứng minh nước mưa tại Hà Nội mang hàm lượng $NO_3\text{-}N$ và $TN$ rất cao, làm đảo chiều tỷ lệ $TN/TP$ từ $< 4,5$ lên $> 6,0$, trực tiếp kích hoạt hiện tượng bùng nổ sinh khối của loài vi khuẩn lam độc hại Microcystis aeruginosa chỉ sau 48-120 giờ sau mưa.

4. Luận án có cung cấp đầy đủ giao thức tái lập (Replication Protocol) không?
Trả lời: Toàn bộ giao thức nghiên cứu được cung cấp chi tiết tuyệt đối trong các phụ lục: từ số liệu quan trắc khí tượng ngày, kết quả phân tích mẫu nước mưa, bảng tọa độ lấy mẫu, ma trận quy đổi sinh khối tảo sang đơn vị Carbon, đến toàn bộ đoạn mã nguồn lập trình MATLAB (he_pt_vi_phan.m và thuật toán giải ode45) trong Phụ lục 11, cho phép bất kỳ nhà nghiên cứu nào cũng có thể tái lập và chạy lại mô hình.

5. Lộ trình nghiên cứu 10 năm được vạch ra như thế nào?
Trả lời: Lộ trình 10 năm định hướng phát triển mô hình qua 3 giai đoạn: (1) 2020-2023: Mở rộng tích hợp module trầm tích 2 lớp phản ứng sinh hóa; (2) 2024-2026: Nâng cấp liên kết không gian 2 chiều (2D hydrodynamic-water quality coupling); (3) 2027-2030: Ứng dụng trí tuệ nhân tạo (AI/PINNs) kết hợp dữ liệu lớn từ trạm quan trắc tự động và ảnh vệ tinh để vận hành hệ thống dự báo phú dưỡng trực tuyến thời gian thực cho mạng lưới hồ toàn quốc.

Kết luận

  1. Phát triển thành công mô hình toán sinh thái chuyên biệt: Xây dựng hoàn chỉnh hệ 14 phương trình vi phân thường mô phỏng đồng thời 22 quá trình động học hóa sinh, tương tác giữa 3 nhóm vi tảo và chu trình dinh dưỡng trong hồ nước tĩnh nông nhiệt đới.
  2. Đột phá trong kỹ thuật tối ưu hóa tham số: Ứng dụng xuất sắc Thuật giải di truyền (GA) để hiệu chỉnh tự động 15 tham số nhạy nhất trong tổng số 88 tham số mô hình, đạt độ chính xác kiểm định vượt bậc ($NSE = 0,72 - 0,84$; $PBIAS < \pm 15%$).
  3. Làm sáng tỏ cơ chế sinh thái học phú dưỡng hồ nội đô: Định lượng được vai trò kích hoạt nở hoa vi khuẩn lam của nguồn dinh dưỡng từ lắng đọng khí quyển và nước mưa chảy tràn đô thị, chỉ ra sự biến dịch tính chất hạn chế dinh dưỡng ($TN/TP$) theo mùa.
  4. Cung cấp công cụ dự báo và đánh giá kịch bản thực tiễn: Ứng dụng thành công tại hồ Cự Chính (Hà Nội), chứng minh kịch bản nạo vét bùn đáy kết hợp kiểm soát nước mưa đợt đầu giúp giảm $74%$ sinh khối tảo độc và đưa hồ từ trạng thái siêu phú dưỡng ($TSI = 75,5$) về trung dưỡng bền vững ($TSI = 46,5$).
  5. Đóng góp nền tảng công nghệ mở: Chuyển giao toàn bộ thuật toán và mã nguồn tính toán, tạo tiền đề vững chắc cho việc tự chủ công nghệ mô hình hóa môi trường nước tại Việt Nam, hội nhập toàn diện với các tiêu chuẩn nghiên cứu thủy vực học quốc tế.