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 tĩnh nông đô thị đang là thách thức sinh thái môi trường nghiêm trọng trên toàn cầu, đặc biệt tại các vùng khí hậu nhiệt đới gió mùa có tốc độ đô thị hóa nhanh như Hà Nội. Luận án tiến sĩ kỹ thuật chuyên ngành Môi trường Đất và Nước (Mã số: 9.03) của NCS. Tạ Đăng Thuần dưới sự hướng dẫn khoa học của PGS.TS. Bùi Quốc Lập tại Trường Đại học Thủy Lợi (2019) 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" đã giải quyết một khoảng trống nghiên cứu (research gap) quan trọng: sự thiếu hụt các công cụ toán học mô phỏng hệ sinh thái thủy vực nông nhiệt đới có tích hợp nguồn ô nhiễm phân tán từ lắng đọng khí quyển và dòng chảy tràn bề mặt đô thị.

Các nghiên cứu trước đây tại Việt Nam và trên thế giới chủ yếu tập trung vào hồ sâu vùng ôn đới hoặc áp dụng các mô hình cân bằng khối lượng đơn biến (Vollenweider, 1968; Lorenzen, 1976). Khi ứng dụng vào hệ thống hồ nội đô Hà Nội – nơi nhiều hồ dù đã được kè bờ và tách hoàn toàn hệ thống nước thải sinh hoạt điểm (point-source pollution) như hồ Hoàn Kiếm, Thiền Quang, Thủ Lệ, Thành Công, Hai Bà Trưng và Cự Chính – hiện tượng nở hoa nước (algal bloom) do vi khuẩn lam (Cyanobacteria) vẫn diễn biến phức tạp. Nguyên nhân xuất phát từ sự tích tụ dinh dưỡng phân tán (non-point sources) và quá trình tái giải phóng nội tại từ trầm tích đáy.

Luận án thiết lập hệ thống câu hỏi nghiên cứu và giả thuyết khoa học cụ thể:

  1. Câu hỏi nghiên cứu 1 (RQ1): Động học chuyển hóa các dạng hợp chất nitơ, phốt pho, carbon hữu cơ và oxy hòa tan trong vùng nước tĩnh nông nhiệt đới chịu sự chi phối như thế nào bởi các biến ngoại sinh khí tượng (nhiệt độ, bức xạ mặt trời, lượng mưa) và nguồn tải lượng phân tán?
  2. Câu hỏi nghiên cứu 2 (RQ2): Làm thế nào để mô hình hóa sự cạnh tranh sinh thái giữa 3 nhóm thực vật nổi chủ đạo (tảo lục, vi khuẩn lam, tảo silic) và tương tác bậc dinh dưỡng với động vật phù du trong điều kiện hồ nông pha trộn hoàn toàn?
  3. Giả thuyết khoa học 1 (H1): Việc bổ sung trực tiếp thông lượng dinh dưỡng từ khí quyển ($L_{atm}$) và nước mưa chảy tràn đô thị ($L_{runoff}$) vào hệ phương trình vi phân thường (ODE) động thái dinh dưỡng sẽ cải thiện đáng kể độ chính xác mô phỏng so với mô hình sinh thái kín.
  4. Giả thuyết khoa học 2 (H2): Trong điều kiện nhiệt độ cao ($>28^\circ\text{C}$) và tỷ lệ $\text{TN}/\text{TP} < 7,2$, vi khuẩn lam chiếm ưu thế áp đảo so với tảo lục và tảo silic nhờ cơ chế thích nghi quang học và động học hấp thụ dinh dưỡng cạnh tranh.

Khung lý thuyết của luận án được xây dựng dựa trên sự tích hợp của thuyết cân bằng hóa học sinh thái Redfield (Redfield stoichiometry), động học enzyme Michaelis–Menten mở rộng cho nhiều cơ chất giới hạn và mô hình tăng trưởng tảo Monod. Phạm vi thực nghiệm tập trung tại hồ Cự Chính (quận Thanh Xuân, Hà Nội) với diện tích mặt nước $4.000,\text{m}^2$, chu vi $250,\text{m}$, độ sâu trung bình $1,5 - 1,7,\text{m}$, dung tích trữ nước $4.600,\text{m}^3$, được theo dõi liên tục trong chu kỳ mùa khô và mùa mưa qua chuỗi quan trắc 2017–2018 với 53 tham số sinh thái học được nhận dạng và tối ưu hóa.

Literature Review và Positioning

Tổng quan y văn quốc tế ghi nhận sự phát triển của mô hình phú dưỡng qua ba giai đoạn tiến hóa chính:

  1. Trường phái mô hình thực nghiệm và cân bằng vật chất đơn giản: Khởi xướng bởi Vollenweider (1968, 1976), tiếp nối bởi Dillon & Rigler (1974), Larsen & Mercier (1976), Chapra (1975) và Lorenzen (1976). Các mô hình này giả định hồ ở trạng thái pha trộn hoàn toàn, ổn định dài hạn và lấy tổng phốt pho ($\text{TP}$) làm yếu tố giới hạn duy nhất ($1\text{D}/1\text{-box model}$). Mặc dù hữu ích cho các hồ lớn vùng ôn đới Bắc Mỹ và Châu Âu, cách tiếp cận này thất bại khi áp dụng cho các hồ nông nhiệt đới biến động nhanh do bỏ qua tương tác nitơ, carbon, sự phân tầng nhiệt vi mô và giải phóng dinh dưỡng nội tại từ bùn đáy (Chapra & Canale, 1991; Hakanson et al., 2007).

  2. Trường phái mô hình thủy động lực học - chất lượng nước phức tạp: Điển hình là EFDC (Hamrick, 1992), CE-QUAL-W2 (Cole & Wells, 2000) và CAEDYM (Hipsey et al., 2006). Các mô hình 3D này giải hệ phương trình Navier–Stokes kết hợp vận chuyển khuếch tán. Tuy nhiên, rào cản lớn khi triển khai tại các nước đang phát triển là yêu cầu khối lượng dữ liệu địa hình, khí tượng biên cực lớn, chi phí bản quyền cao và hiện tượng "quá tham số hóa" (over-parameterization) dẫn đến độ bất định cao khi áp dụng cho các hồ đô thị quy mô nhỏ ($<5,\text{ha}$).

  3. Trường phái mô hình sinh thái dinh dưỡng - chuỗi thức ăn (Ecological dynamic models): Nổi bật là mô hình SALMO (Recknagel et al., 2008; Benndorf, 1979) mô phỏng tương tác dinh dưỡng - thực vật nổi - động vật phù du; mô hình hồ nông Taihu của Zhang et al. (2008); mô hình sinh thái hồ nông nghiệp của Do Thuy Nguyen et al. (2010).

   [GIAI ĐOẠN 1: MÔ HÌNH THỰC NGHIỆM ĐƠN LỚP]
   Vollenweider (1968), Larsen (1976), Lorenzen (1976)
   - Trọng tâm: Cân bằng TP, giả định hồ ôn đới sâu
   - Hạn chế: Bỏ qua N, C, chuỗi thức ăn & nguồn phân tán đô thị
   [GIAI ĐOẠN 2: MÔ HÌNH 3D PHỨC TẠP]
   EFDC (Hamrick, 1992), CAEDYM (Hipsey, 2006)
   - Trọng tâm: Thủy động lực học Navier-Stokes 3D
   - Hạn chế: Quá tải tham số, đòi hỏi biên dữ liệu cực lớn
   [GIAI ĐOẠN 3: ĐỘNG HỌC SINH THÁI HỆ NÔNG]
   SALMO (Recknagel, 2008), Do Thuy Nguyen et al. (2010)
   - Trọng tâm: Đa nhóm tảo + ĐVPD + Động học dinh dưỡng
   - Hạn chế: Giả định hệ kín, thiếu nguồn bổ cập khí quyển/chảy tràn
   [VỊ TRÍ ĐỊA VỊ KHOA HỌC CỦA LUẬN ÁN TẠ ĐĂNG THUẦN (2019)]
   Mô hình sinh thái động học vùng nước tĩnh nông nhiệt đới
   - Cải tiến: Tích hợp 14 biến trạng thái + Nguồn mưa/chảy tràn + Thuật giải GA tối ưu 53 tham số

Tại Việt Nam, nghiên cứu của Nguyễn Thanh Hùng (2010) chỉ giải bài toán cân bằng phốt pho 2 lớp cho các hồ chứa lớn, sâu (Đại Lải, Tuyên Quang) phục vụ cấp nước. Nghiên cứu của Nguyễn Xuân Tình và cs. (2013, 2015) ứng dụng EFDC cho Hồ Tây nhưng chưa mô tả đầy đủ các quá trình trao đổi chất, động học đa nhóm tảo và thiếu bộ tham số thực nghiệm tối ưu. Đặc biệt, mô hình 14 biến trạng thái của Do Thuy Nguyen et al. (2010) cho hồ nông nghiệp giả thiết không có nguồn chất dinh dưỡng từ bên ngoài bổ cập, dẫn đến sai số lớn khi mô phỏng $\text{NH}_4\text{-N}$, $\text{NO}_3\text{-N}$ và $\text{PO}_4\text{-P}$.

Luận án của Tạ Đăng Thuần định vị tại điểm giao thoa giữa độ phức tạp sinh thái học và tính khả thi tính toán, phát triển một mô hình sinh thái mở 14 biến trạng thái, hiệu chỉnh trực tiếp cho điều kiện khí hậu nhiệt đới gió mùa của Hà Nội.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đóng góp quan trọng vào lý thuyết động thái chất lượng nước và sinh thái học thủy vực tĩnh thông qua việc mở rộng các lý thuyết nền tảng:

  • Mở rộng lý thuyết tỷ lệ Redfield trong môi trường phú dưỡng cực hạn: Tỷ lệ nguyên tử kinh điển $106\text{C}:16\text{N}:1\text{P}$ tương ứng với tỷ lệ khối lượng $\text{C}:\text{N}:\text{P} = 40:7,2:1$ được luận án kiểm chứng thực nghiệm. Kết quả chỉ ra rằng khi $\text{TN}/\text{TP} < 4,5$, nitơ đóng vai trò là yếu tố giới hạn tuyệt đối; trong khoảng $4,5 \le \text{TN}/\text{TP} \le 6,0$, cả nitơ và phốt pho đồng giới hạn; và khi $\text{TN}/\text{TP} > 6,0$, phốt pho trở thành chất dinh dưỡng giới hạn sinh trưởng.
  • Thiết lập phương trình quang hợp phân hóa theo cơ chất nitơ vô cơ: Luận án vận dụng hai cơ chế quang hợp dựa trên nguồn hấp thụ amoni ($\text{NH}_4^+$) và nitrat ($\text{NO}_3^-$):

$$106\text{CO}2 + 16\text{NH}4^+ + \text{HPO}4^{2-} + 106\text{H}2\text{O} \xrightarrow{\text{Năng lượng, NTVL}} \text{C}{106}\text{H}{181}\text{O}{45}\text{N}{16}\text{P} + 106\text{O}_2 + 15\text{H}^+$$

$$106\text{CO}2 + 16\text{NO}3^- + \text{HPO}4^{2-} + 122\text{H}2\text{O} + 17\text{H}^+ \xrightarrow{\text{Năng lượng, NTVL}} \text{C}{106}\text{H}{181}\text{O}{45}\text{N}{16}\text{P} + 138\text{O}_2$$

Chứng minh rằng quá trình đồng hóa $\text{NO}_3^-$ giải phóng lượng oxy hòa tan cao hơn $30,18%$ so với đồng hóa $\text{NH}_4^+$, ảnh hưởng trực tiếp đến động thái DO ban ngày trong tầng nước mặt.

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

Khung phân tích của mô hình bao gồm hệ thống 14 biến trạng thái phi tuyến liên kết chặt chẽ:

  • 02 biến carbon hữu cơ: Carbon hữu cơ dạng hạt ($\text{POC}$), carbon hữu cơ hòa tan ($\text{DOC}$).
  • 05 biến dinh dưỡng vô cơ và dưỡng khí: Oxy hòa tan ($\text{DO}$), amoni ($\text{NH}_4\text{-N}$), nitrite ($\text{NO}_2\text{-N}$), nitrate ($\text{NO}_3\text{-N}$), phosphate vô cơ hòa tan ($\text{PO}_4\text{-P}$ hay $\text{DIP}$).
  • 05 biến sinh học quần xã: Sinh khối tảo lục ($\text{A}_1$ - Chlorophyta), sinh khối vi khuẩn lam ($\text{A}_2$ - Cyanobacteria), sinh khối tảo silic ($\text{A}_3$ - Bacillariophyta), sinh khối tảo lông hai roi ($\text{A}_4$), sinh khối động vật phù du ($\text{CZP}$).
  • 02 biến tổng hợp: Tổng nitơ ($\text{TN}$) và Tổng phốt pho ($\text{TP}$).

Điểm đột phá cấu trúc là phương trình vi phân bổ sung các số hạng nguồn phân tán ngoại sinh:

$$\frac{dC_i}{dt} = \sum (\text{Phản ứng sinh hóa nội tại}) + \frac{L_{atm,i} \cdot A_{lake}}{V} + \frac{L_{runoff,i} \cdot A_{catchment} \cdot P \cdot \varphi}{V} - K_{sed,i} C_i$$

Trong đó: $C_i$ là nồng độ biến trạng thái $i$; $L_{atm,i}$ là suất tải lắng đọng khí quyển; $L_{runoff,i}$ là nồng độ trong nước mưa chảy tràn; $A_{lake}, A_{catchment}$ lần lượt là diện tích mặt hồ ($4.000,\text{m}^2$) và lưu vực thu nước ($1.500,\text{m}^2$ sân bê tông); $P$ là lượng mưa ngày ($\text{mm/ngày}$); $\varphi$ là hệ số dòng chảy mặt ($\varphi = 0,9$ đối với bề mặt bê tông); $V$ là dung tích hồ ($4.600,\text{m}^3$).

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

Thiết kế nghiên cứu

Nghiên cứu áp dụng lập trường nhận thức luận thực chứng (Positivism / Post-positivism) kết hợp mô hình hóa tất định (Deterministic ecological modeling) và kỹ thuật tối ưu hóa ngẫu nhiên (Stochastic optimization).

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

  1. Quan trắc thực địa và phân tích hóa nghiệm chất lượng nước định kỳ theo mùa.
  2. Định lượng sinh thái phân loại học thực vật nổi và vi khuẩn lam độc hại.
  3. Rời rạc hóa hệ phương trình vi phân thường phi tuyến bằng thuật toán Runge–Kutta bậc 4–5 thích ứng bước thời gian ($\text{ODE45}$).
  4. Tối ưu hóa toàn cục bộ tham số mô hình thông qua Thuật giải Di truyền (Genetic Algorithm - GA).

Quy trình nghiên cứu rigorous

Quy trình thu mẫu và phân tích 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):

  • Thiết bị đo nhanh hiện trường: Nồng độ DO, nhiệt độ nước, pH, độ dẫn điện (EC) được đo trực tiếp bằng máy đo đa chỉ tiêu HACH HQ40D.
  • Phân tích hóa nghiệm phòng thí nghiệm:
    • $\text{TN}$: Vô cơ hóa mẫu bằng Persulfate, đo quang phổ hấp thụ phân tử UV-Vis.
    • $\text{TP}$: Vô cơ hóa bằng Kali persulfate trong nồi hấp tiệt trùng ($121^\circ\text{C}, 30,\text{phút}$), tạo phức màu xanh Molybdate, đo quang ở bước sóng $880,\text{nm}$.
    • $\text{NH}_4\text{-N}$: Phương pháp so màu Phenate.
    • $\text{NO}_3\text{-N}$: Phương pháp trắc quang dùng thanh khử Cadmium.
    • $\text{PO}_4\text{-P}$: Phương pháp Axit Ascorbic.
    • $\text{Chl.a}$: Lọc mẫu qua màng lọc sợi thủy tinh GF/C ($0,45,\mu\text{m}$), chiết xuất sắc tố bằng cồn Acetone $90%$, đo mật độ quang ở bước sóng $664,\text{nm}, 647,\text{nm}, 630,\text{nm}$.
  • Phân tích sinh học: Định danh chi loài thực vật nổi dưới kính hiển vi quang học Olympus BX53 ở độ phóng đại $400\times - 1000\times$, đếm mật độ tế bào bằng buồng đếm Sedgewick-Rafter.

Data và phân tích

Số liệu khí tượng thu thập liên tục từ Trạm Khí tượng Láng (cách hồ Cự Chính $2,\text{km}$):

  • Nhiệt độ trung bình năm: $24,5^\circ\text{C}$ (tháng 6 đạt cực đại $30,1^\circ\text{C}$; tháng 1 thấp nhất $16,5^\circ\text{C}$).
  • Số giờ nắng trung bình năm: $1.292,\text{giờ}$ (tháng 6 đạt $156,\text{giờ}$, tháng 3 thấp nhất $42,\text{giờ}$).
  • Độ ẩm không khí trung bình: $78%$.
  • Lượng mưa trung bình năm: $1.652,\text{mm}$ (mùa mưa từ tháng 5 đến tháng 9 chiếm $>80%$ tổng lượng mưa năm; tháng 8 đạt đỉnh $348,\text{mm}$).

Để giải quyết bài toán tối ưu hóa 53 tham số động học, tác giả xây dựng hàm mục tiêu (Fitness Function) dựa trên tổng sai số bình phương chuẩn hóa giữa giá trị mô phỏng ($Y_i^{sim}$) và thực đo ($Y_i^{obs}$):

$$\text{Fitness} = \sum_{j=1}^{M} w_j \left[ \frac{1}{N_j} \sum_{k=1}^{N_j} \left( \frac{Y_{j,k}^{obs} - Y_{j,k}^{sim}}{\sigma_j^{obs}} \right)^2 \right]$$

Trong đó $w_j$ là trọng số của biến trạng thái $j$; $\sigma_j^{obs}$ là độ lệch chuẩn của chuỗi quan trắc thực tế.

Quy trình tối ưu hóa được lập trình trên nền tảng MATLAB R2017a:

  • Phân tích độ nhạy (Sensitivity Analysis) sàng lọc được 15 tham số nhạy nhất trong tổng số 53 tham số.
  • Thuật giải Di truyền (GA) thực hiện 38 lần chạy lặp độc lập với kích thước quần thể (Population size) = 100, xác suất lai ghép (Crossover rate) = 0,8, xác suất đột biến (Mutation rate) = 0,05.
  • Bộ 4 chỉ số thống kê chuẩn quốc tế được sử dụng để đánh giá chất lượng mô hình theo hướng dẫn của Moriasi et al. (2007):

$$\text{RMSE} = \sqrt{\frac{1}{n}\sum_{i=1}^n (O_i - P_i)^2}$$

$$\text{NSE} = 1 - \frac{\sum_{i=1}^n (O_i - P_i)^2}{\sum_{i=1}^n (O_i - \bar{O})^2}$$

$$\text{RSR} = \frac{\text{RMSE}}{\text{STDEV}{obs}} = \frac{\sqrt{\sum{i=1}^n (O_i - P_i)^2}}{\sqrt{\sum_{i=1}^n (O_i - \bar{O})^2}}$$

$$\text{PBIAS} = \left[ \frac{\sum_{i=1}^n (O_i - P_i)}{\sum_{i=1}^n O_i} \right] \times 100%$$

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

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

  1. Xác lập hiện trạng siêu phú dưỡng và động thái dinh dưỡng giới hạn: Theo phân loại Hakanson (2007), hồ Cự Chính nằm hoàn toàn trong ngưỡng siêu phú dưỡng (Hypereutrophic) với nồng độ $\text{TN} > 0,43,\text{mg/L}$, $\text{TP} > 0,06,\text{mg/L}$ và $\text{Chl.a} > 20,\mu\text{g/L}$. Chỉ số Carlson TSI dao động từ $68,5$ đến $82,4$ (trung bình $\text{TSI} > 70$). Tỷ lệ $\text{TN}/\text{TP}$ trung bình đạt $5,2 \pm 1,1$, chứng minh hồ nằm trong giai đoạn đồng giới hạn hoặc chuyển tiếp sang giới hạn nitơ vào mùa hè nóng ẩm.

  2. Cơ cấu quần xã thực vật nổi và cảnh báo bùng nổ độc tố Microcystins: Phát hiện 4 nhóm thực vật nổi chính với 32 chi loài, trong đó vi khuẩn lam (Cyanobacteria) và tảo lục (Chlorophyta) chiếm ưu thế tuyệt đối ($>85%$ tổng mật độ tế bào). Chi Microcystis (thuộc VKL) chiếm tỷ lệ sinh khối cao nhất trong mùa hè khi nhiệt độ nước $>29^\circ\text{C}$ và bức xạ mặt trời đạt đỉnh ($>18,\text{MJ/m}^2/\text{ngày}$). Dẫn liệu so sánh từ hồ Hoàn Kiếm cho thấy mật độ Microcystis đạt $386,94 \times 10^6,\text{tế bào/mL}$ sinh ra độc tố gan Microcystins dao động từ $2,1 - 46,\mu\text{g/L}$, vượt xa ngưỡng an toàn $1,0,\mu\text{g/L}$ của Tổ chức Y tế Thế giới (WHO).

  3. Độ chính xác vượt trội của mô hình sau hiệu chỉnh tối ưu hóa bằng GA: Quá trình tối ưu hóa 53 tham số bằng GA đã đưa các chỉ số thống kê kiểm định về mức xuất sắc theo tiêu chuẩn Moriasi:

  • Hệ số hiệu quả Nash–Sutcliffe ($\text{NSE}$) của các biến sinh khối tảo lục, vi khuẩn lam, $\text{DO}$ và $\text{PO}4\text{-P}$ đều đạt $> 0,75$ (trong đó $\text{NSE}{\text{Cyanobacteria}} = 0,82$).
  • Tỷ số sai số $\text{RSR} < 0,50$.
  • Phần trăm độ lệch $\text{PBIAS} < \pm 15%$.
BẢNG SO SÁNH HIỆU QUẢ CỦA MÔ HÌNH TRƯỚC VÀ SAU KHI ĐIỀU CHỈNH 53 THAM SỐ QUA THUẬT GIẢI DI TRUYỀN (GA)
  1. Định lượng đóng góp của tải lượng nước mưa chảy tràn đô thị: Kết quả mô phỏng chỉ ra rằng trong các trận mưa rào đầu mùa, lượng dinh dưỡng cuốn trôi từ $1.500,\text{m}^2$ lưu vực bê tông xung quanh đóng góp tới $32,4%$ tổng tải lượng $\text{DIP}$ và $28,6%$ $\text{DIN}$ gia nhập hồ, kích hoạt quá trình bùng nổ sinh khối vi khuẩn lam chỉ sau $3 - 5,\text{ngày}$ hửng nắng.

Implications đa chiều

  • Về lý thuyết: Khẳng định tính tất yếu của việc đưa động học lắng đọng phân tán vào các mô hình chất lượng nước hồ nội đô; mở rộng hiểu biết về cơ chế bù trừ dinh dưỡng giữa vi khuẩn lam và tảo lục dưới biến động nhiệt độ vi khí hậu.
  • Về phương pháp luận: Cung cấp quy trình hoàn chỉnh kết hợp giải thuật vi phân số trị Runge–Kutta với trí tuệ nhân tạo GA trong nhận dạng tham số sinh thái học thủy vực, loại bỏ tính chủ quan của phương pháp thử sai truyền thống.
  • Về thực tiễn quản lý: Mô phỏng 4 kịch bản kỹ thuật can thiệp cho hồ Cự Chính:
    • Kịch bản 1 (Nạo vét bùn đáy định kỳ $20,\text{cm}$): Giảm tốc độ giải phóng $\text{PO}_4\text{-P}$ từ trầm tích, hạ sinh khối tảo trung bình $22,8%$.
    • Kịch bản 2 (Thu gom $100%$ nước mưa chảy tràn sân hồ): Cắt giảm đỉnh bùng phát tảo sau mưa $31,5%$.
    • Kịch bản 3 (Sục khí cưỡng bức tăng cường DO tầng đáy $>4,\text{mg/L}$): Ức chế giải phóng phốt pho yếm khí và khử khí độc $\text{H}_2\text{S}, \text{CH}_3\text{SH}$.
    • Kịch bản 4 (Tích hợp nạo vét + thu gom dòng tràn + thả bè thực vật thủy sinh Thủy trúc): Giảm sinh khối vi khuẩn lam $58,4%$, đưa trạng thái hồ từ siêu phú dưỡng về mức phú dưỡng vừa (TSI $< 55$).

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 các giới hạn khoa học:

  1. Giả định pha trộn hoàn toàn 0D/1D: Mô hình giả thiết hồ đồng nhất về mặt không gian, chưa phản ánh được sự phân tầng nhiệt vi mô cục bộ (Micro-stratification) vào những ngày mùa hè lặng gió trong các góc hồ khuất.
  2. Chưa mô hình hóa chi tiết độc tố sinh học: Luận án mới dừng lại ở việc mô phỏng sinh khối vi khuẩn lam ($A_2$) mà chưa thiết lập phương trình riêng biệt cho quá trình sinh tổng hợp và phân hủy của các biến thể độc tố Microcystin-LR/RR trong môi trường nước.
  3. Độ phân giải chuỗi số liệu thực nghiệm: Tần suất lấy mẫu phân tích dinh dưỡng định kỳ 2 tuần/lần; chưa có hệ thống cảm biến đo đạc tự động liên tục (Real-time IoT sensors) để bắt trọn các xung dinh dưỡng tức thời trong các trận mưa cực đoan.

Định hướng nghiên cứu mở rộng 5–10 năm tới bao gồm:

  • Tích hợp module thủy động lực học 2D nông (Shallow Water Equations) với mô hình sinh thái để mô phỏng sự phân bố dị hướng của thảm tảo dưới tác động của trường gió đô thị.
  • Ứng dụng kỹ thuật viễn thám độ phân giải cao (Sentinel-2, PlanetScope) kết hợp giải thuật học máy (Machine Learning/Deep Learning) để đồng hóa dữ liệu (Data Assimilation) vào mô hình toán.
  • Mở rộng bài toán biến đổi khí hậu: Đánh giá tác động của các đợt sóng nhiệt (Heatwaves) cực đoan đến tần suất nở hoa của vi khuẩn lam độc hại tại hệ thống hồ cảnh quan đô thị miền Bắc Việt Nam.

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

Luận án tạo ra những tác động khoa học và xã hội rõ nét:

  • Giá trị học thuật: Đóng góp công cụ mô hình sinh thái mã nguồn mở viết trên nền tảng MATLAB, có khả năng tùy biến cao cho các viện nghiên cứu và trường đại học chuyên ngành Thủy lợi, Môi trường, Tài nguyên Nước. Dự báo tạo lập nguồn trích dẫn học thuật vững chắc trong lĩnh vực mô hình hóa sinh thái hồ nhiệt đới.
  • Chính sách và quản lý đô thị: Cung cấp cơ sở khoa học định lượng giúp Sở Tài nguyên và Môi trường Hà Nội, Công ty TNHH MTV Thoát nước Hà Nội xây dựng quy trình vận hành và cải tạo hệ thống hồ nội đô (kè hồ sinh thái, thiết kế hố thu nước mưa tách dòng, nạo vét bùn đáy định kỳ).
  • Lợi ích xã hội và sức khỏe cộng đồng: Kiểm soát hiện tượng bùng phát vi khuẩn lam độc, loại bỏ nguy cơ phát tán khí độc $\text{H}_2\text{S}$ và độc tố gan gây ô nhiễm không khí và đe dọa sức khỏe cư dân ven hồ.

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

  • Nghiên cứu sinh và học viên cao học: Tiếp cận khung phương pháp luận kết hợp mô hình vi phân sinh thái và thuật giải di truyền GA; kế thừa bộ 53 tham số sinh thái học làm giá trị tham chiếu chuẩn cho các nghiên cứu hồ nhiệt đới.
  • Nhà khoa học và giảng viên chuyên ngành: Tài liệu giảng dạy chuyên sâu về Mô hình hóa chất lượng nước, Động thái sinh thái học thủy vực nông và Quản lý phú dưỡng.
  • Cơ quan quản lý môi trường đô thị: Khung công cụ ra quyết định (Decision Support Framework) đánh giá hiệu quả chi phí trước khi triển khai các dự án cải tạo hồ tốn kém hàng tỷ đồng.
  • Doanh nghiệp công nghệ xử lý môi trường: Định lượng chính xác công suất thiết kế hệ thống sục khí tầng đáy, bè bè sinh học thực vật thủy sinh phù hợp với tải lượng ô nhiễm nội tại của từng hồ.

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ì? Đó là việc mở rộng hệ phương trình động thái sinh thái học vùng nước tĩnh nông nhiệt đới thông qua việc tích hợp đồng thời chu trình 14 biến trạng thái dinh dưỡng (C, N, P, DO, 3 nhóm tảo, ĐVPD) với nguồn bổ cập xung phân tán từ khí quyển và dòng chảy tràn đô thị bê tông hóa.

  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 Nguyễn Thanh Hùng (2010) chỉ giải bài toán cân bằng TP cho hồ sâu và Do Thuy Nguyen et al. (2010) giả định hệ sinh thái nông nghiệp kín, luận án đã tích hợp thuật toán Runge–Kutta bậc 4–5 với thuật giải di truyền GA tự động tối ưu hóa 53 tham số trên MATLAB, giải quyết triệt để bài toán cực tiểu hóa hàm mục tiêu sai số.

  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ì? Mặc dù hồ Cự Chính không có bất kỳ nguồn nước thải sinh hoạt điểm nào đổ trực tiếp vào lòng hồ, mức độ phú dưỡng vẫn duy trì ở ngưỡng siêu phú dưỡng ($\text{TSI} > 70$) quanh năm. Nguyên nhân chính được chứng minh là do "ký ức dinh dưỡng" từ lớp bùn đáy dày giải phóng liên tục $\text{PO}_4\text{-P}$ kết hợp với các xung nitơ từ nước mưa chảy tràn qua bề mặt sân bê tông.

  4. Nghiên cứu có cung cấp giao thức tái lặp (Replication Protocol) hoàn chỉnh không? Có. Luận án cung cấp chi tiết toàn bộ hệ thống 14 phương trình vi phân thường, bảng giá trị và khoảng biến thiên của 53 tham số sinh thái, sơ đồ khối thuật toán tối ưu hóa GA trong MATLAB và quy trình phân tích mẫu chuẩn hóa TCVN/APHA, cho phép tái lập quy trình trên bất kỳ hồ đô thị nào.

  5. Chương trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào? Chuyển đổi mô hình từ 0D/1D sang 2D thủy động lực học - sinh thái tích hợp, phát triển hệ thống giám sát cảnh báo sớm tảo nở hoa theo thời gian thực dựa trên mạng lưới cảm biến IoT kết hợp trí tuệ nhân tạo và ảnh viễn thám vệ tinh phân giải siêu cao.

Kết luận

  1. Xây dựng thành công mô hình toán sinh thái mở 14 biến trạng thái mô phỏng toàn diện quá trình phú dưỡng trong các vùng nước tĩnh nông nhiệt đới, có tính toán đến tải lượng lắng đọng khí quyển và dòng chảy tràn đô thị.
  2. Ứng dụng xuất sắc thuật giải di truyền GA trên nền tảng MATLAB R2017a để tối ưu hóa toàn cục 53 tham số động học, nâng cao chỉ số hiệu quả Nash–Sutcliffe ($\text{NSE}$) đạt mức $>0,78 - 0,86$, vượt trội so với các phương pháp thử sai truyền thống.
  3. Làm sáng tỏ bức tranh sinh thái hồ Cự Chính: trạng thái siêu phú dưỡng ($\text{TSI} = 68,5 - 82,4$), cơ chế đồng giới hạn dinh dưỡng ($\text{TN}/\text{TP} = 5,2$) và nguy cơ nở hoa vi khuẩn lam Microcystis sản sinh độc tố gan Microcystins vào mùa nóng.
  4. Đề xuất và mô phỏng định lượng hệ giải pháp kỹ thuật tích hợp (nạo vét bùn đáy định kỳ, kiểm soát dòng tràn bề mặt, sục khí cưỡng bức và bè thủy trúc), chứng minh khả năng cắt giảm $58,4%$ sinh khối tảo độc, khôi phục chất lượng môi trường nước hồ đô thị bền vững.
  5. Đặt nền móng khoa học vững chắc cho trường phái mô hình hóa sinh thái thủy vực tại Việt Nam, mở ra hướng ứng dụng công nghệ số và trí tuệ nhân tạo trong bảo tồn tài nguyên nước đô thị trước thách thức của biến đổi khí hậu toàn cầu.