Tổng quan luận án

Tính cấp thiết của đề tài

Tại các vùng nước tĩnh nông (ao hồ tự nhiên, hồ nông nghiệp có độ sâu trung bình dưới 5 m), do mức độ trao đổi nước hạn chế cùng với ảnh hưởng từ gió bề mặt, các chất dinh dưỡng và chất ô nhiễm từ lưu vực bị cuốn trôi theo dòng chảy tràn tích tụ dần trong thủy vực, dẫn đến suy giảm chất lượng nước nghiêm trọng. Vấn đề môi trường phổ biến nhất tại các thủy vực này là hiện tượng phú dưỡng. Phú dưỡng gây ra sự bùng phát quá mức của rong, tảo và thực vật thủy sinh, gia tăng vi khuẩn lam và tảo độc. Khi sinh khối này chết đi, quá trình phân hủy làm suy giảm nồng độ oxy hòa tan (DO), tạo ra các hợp chất độc hại như nitrite ($\text{NO}_2\text{-N}$), nitrate ($\text{NO}_3\text{-N}$), đồng thời giải phóng các khí gây mùi độc hại trong điều kiện yếm khí như hydro sulfua ($\text{H}_2\text{S}$) và mercaptane ($\text{CH}_3\text{SH}$). Xác thực vật thủy sinh lắng đọng làm lòng hồ ngày càng nông dần, suy giảm dung tích trữ nước.

Để kiểm soát và dự báo hiện tượng này, việc ứng dụng mô hình toán mô phỏng quá trình phú dưỡng mang lại nhiều ưu thế: thời gian tính toán nhanh, chi phí thấp, cho phép mô phỏng các kịch bản quản lý chất lượng nước và khắc phục khó khăn của các thí nghiệm hiện trường chịu nhiều yếu tố nhiễu. Tuy nhiên, phần lớn các mô hình phú dưỡng trên thế giới được xây dựng cho các hồ sâu, diện tích lớn tại vùng khí hậu ôn đới, chưa thực sự phù hợp với các hồ nhỏ, nông ở vùng nhiệt đới và cận nhiệt đới. Tại Hà Nội, các hồ nội đô đóng vai trò điều hòa vi khí hậu, trữ nước mưa và tạo cảnh quan, phần lớn có kích thước nhỏ và nông. Nhiều hồ dù đã được tách và ngăn chặn hoàn toàn nguồn nước thải sinh hoạt trực tiếp vẫn xảy ra hiện tượng phú dưỡng và tảo nở hoa do nguồn ô nhiễm phân tán (lắng đọng khí quyển, nước mưa chảy tràn bề mặt) và giải phóng dinh dưỡng nội tại từ trầm tích. Do đó, đề 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" được thực hiện nhằm giải quyết các yêu cầu thực tiễn và khoa học nêu trên.

Mục tiêu nghiên cứu

Luận án xác định hai mục tiêu nghiên cứu cụ thể:

  1. Phát triển công cụ mô hình toán mô phỏng quá trình phú dưỡng ở các vùng nước tĩnh nông.
  2. Ứng dụng mô hình phú dưỡng phát triển vào một hồ tự nhiên nông đang bị ảnh hưởng bởi phú dưỡng ở nội thành Hà Nội (hồ Cự Chính).

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Quá trình phú dưỡng và các thông số sinh khối thực vật nổi (tảo lục, vi khuẩn lam, tảo silic), động vật phù du và các thông số chất lượng nước bao gồm các dạng chất dinh dưỡng của nitơ ($\text{NH}_4\text{-N}$, $\text{NO}_2\text{-N}$, $\text{NO}_3\text{-N}$, nitơ hữu cơ), phốt pho (phốt pho vô cơ hòa tan - DIP, phốt pho hữu cơ - OP), carbon (carbon hữu cơ dạng hạt - POC, carbon hữu cơ hòa tan - DOC) và nồng độ oxy hòa tan (DO) ở vùng nước tĩnh nông.
  • Phạm vi nghiên cứu:
    • Không gian: Hồ Cự Chính, nằm giáp ranh giữa hai phường Thượng Đình và Nhân Chính, quận Thanh Xuân, thành phố Hà Nội. Hồ có diện tích mặt nước khoảng 4.000 $\text{m}^2$, chu vi khoảng 250 m, chiều dài lớn nhất khoảng 75 m, chiều rộng lớn nhất khoảng 65 m, độ sâu trung bình từ 1,5–1,7 m, dung tích nước trung bình khoảng 4.600 $\text{m}^3$, khuôn viên bê tông xung quanh hồ rộng khoảng 1.500 $\text{m}^2$.
    • Thời gian và dữ liệu khí tượng: Chuỗi số liệu khí tượng quan trắc tại trạm Láng (Hà Nội, cách khu vực nghiên cứu khoảng 2 km) từ năm 2002 đến năm 2018; các đợt khảo sát, lấy mẫu phân tích chất lượng nước và thực vật nổi thực địa tại hồ Cự Chính phục vụ giai đoạn hiệu chỉnh và kiểm định mô hình.

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu mô hình phú dưỡng trên thế giới và tại Việt Nam

Nghiên cứu phát triển mô hình phú dưỡng hồ bắt đầu từ thập niên 1970 và phát triển mạnh mẽ qua các giai đoạn:

  • Mô hình thực nghiệm: Điển hình là mô hình của Vollenweider (1968, 1975, 1976), thiết lập mối tương quan đồ thị giữa tải lượng phốt pho và trạng thái dinh dưỡng dựa trên số liệu của hơn 200 hồ ôn đới tại châu Âu và Bắc Mỹ. Mô hình giả thiết hồ pha trộn hoàn toàn, ở trạng thái ổn định và mặc định phốt pho là yếu tố giới hạn. Các mô hình dạng này dễ sử dụng nhưng cho sai số lớn khi áp dụng ngoài phạm vi hiệu chỉnh.
  • Mô hình cân bằng vật chất: Tiếp cận theo nguyên lý cân bằng khối lượng các chất dinh dưỡng, tiêu biểu là các nghiên cứu của Imboden (1974), O’Melia (1972), Larsen (1973), Lorenzen (1976), Chapra (1977, 1997), Ruley và Rusch (2002). Các mô hình này phát triển từ việc xem lắng đọng là tổn thất một chiều tới việc mô tả hệ thống 2 lớp (cột nước và trầm tích đáy) có xét đến quá trình khuếch tán ngược chất dinh dưỡng từ bùn đáy lên cột nước (Chapra và Canale, 1991; Lorenzen, 1976).
  • Mô hình sinh thái và động thái chuỗi thức ăn: Mô tả động thái tương tác sinh - hóa phức tạp giữa các dạng chất dinh dưỡng (hòa tan/dạng hạt, vô cơ/hữu cơ, bền/không bền) và các bậc dinh dưỡng sinh học. Tiêu biểu gồm:
    • Mô hình Thomann (1975), Di Toro (1971), Jørgensen (1976), Cleaner (Scavia, 1979), Salmo (Benndorf, 1979; Recknagel et al., 2008; Cao et al., 2008).
    • Mô hình CAEDYM phân chia chi tiết các dạng nitơ ($\text{NH}_4$, $\text{NO}_2$, $\text{NO}_3$), phốt pho, các nhóm tảo và động vật phù du ăn mồi.
    • Mô hình Reynolds (1993, 2001) ứng dụng động học Michaelis–Menten mô tả phản ứng của các nhóm tảo với ánh sáng, nhiệt độ.
    • Mô hình kết hợp thủy động lực học và sinh thái như CE-QUAL-W2 (Cole và Wells, 2002), Fragoso et al. (2009), James (2005, 2006) tại hồ Okeechobee, Cerco và Cole (1993) tại vịnh Chesapeake, Zhang et al. (2008) tại hồ Taihu.
    • Nghiên cứu của Do Thuy Nguyen et al. (2010) phát triển mô hình cho hồ nông nghiệp tĩnh nông mô phỏng 14 biến trạng thái, 2 biến ngoại sinh (nhiệt độ, bức xạ mặt trời), nhưng giả thiết không có nguồn dinh dưỡng bổ sung từ bên ngoài.
Tên mô hình Số phân đoạn/lớp Số chiều Số biến trạng thái Chất dinh dưỡng mô phỏng Hiệu chỉnh (C) / Kiểm định (V)
Vollenweider 1 lớp 1 1 P (N) C + V
Imboden 2 lớp 1 2 P C + V
O’Melia 1 lớp 1 2 P C
Larsen 1 lớp 1 3 P C
Lorenzen 1 lớp 1 2 P C + V
Thomann 2 lớp 1 8 P, N, C C + V
Di Toro 1 lớp 1 7 P, N C + V
Jørgensen 1–2 chiều 1–2 17 P, N, C C + V
Cleaner Vài lớp - 40 P, N, C, Si C
Salmo 2 lớp 1 17 P, N C + V
  • Nghiên cứu tại Việt Nam:
    • Nguyễn Thanh Hùng (2010) xây dựng mô hình tính toán quá trình phú dưỡng cho hệ thống hồ chứa á nhiệt đới (hồ Tây, hồ Đại Lải, hồ Tuyên Quang) giải hệ 2 phương trình vi phân thường cân bằng tổng phốt pho (TP) bằng phương pháp Runge–Kutta bậc 4. Tuy nhiên mô hình chỉ áp dụng cho hồ chứa sâu, dung tích lớn và chỉ xét biến TP.
    • Nguyễn Xuân Tình et al. (2012, 2013, 2015) ứng dụng mô hình 3D mã nguồn mở EFDC mô phỏng thủy động lực học và chất lượng nước tại Hồ Tây, nhưng mô hình đòi hỏi dữ liệu biên phức tạp, chưa phân tách các quá trình động học chi tiết của từng nhóm tảo và thiếu bộ thông số hiệu chỉnh tối ưu cho hồ nông.
    • Một số nghiên cứu khác ứng dụng mô hình thực nghiệm Vollenweider và Jørgensen tại hồ Bảy Mẫu và Hồ Tây để tính toán tải lượng P và chỉ số Carlson TSI.

Khoảng trống nghiên cứu luận án lựa chọn giải quyết

  1. Các mô hình sẵn có trên thế giới chủ yếu áp dụng cho hồ sâu, diện tích lớn ở vùng khí hậu ôn đới, hoặc là các mô hình thương mại dạng "hộp đen" đắt tiền, yêu cầu dữ liệu đầu vào phức tạp.
  2. Tại Việt Nam, chưa có mô hình sinh thái toán học chuyên biệt mô tả đầy đủ bản chất các quá trình động học vật lý, hóa học, sinh học của hiện tượng phú dưỡng cho các thủy vực nước tĩnh nhỏ và nông ở đô thị nhiệt đới.
  3. Các nghiên cứu trước đây thường bỏ qua hoặc giả thiết triệt tiêu nguồn bổ cập dinh dưỡng phân tán từ khí quyển (lắng đọng khô, lắng đọng ướt) và dòng chảy tràn bề mặt đệm đô thị, trong khi đây là nguồn cung cấp nitơ và phốt pho chủ yếu duy trì trạng thái siêu phú dưỡng tại các hồ nội đô đã kiểm soát cống xả thải sinh hoạt.

Cơ sở lý thuyết và phương pháp nghiên cứu

Cơ sở lý thuyết và khung phân tích

  • Cơ chế quang hợp và nhu cầu dinh dưỡng: Sinh khối thực vật nổi được tổng hợp qua quá trình quang hợp sử dụng năng lượng mặt trời, $\text{CO}_2$, nước và các muối vô cơ hòa tan. Tùy thuộc vào nguồn nitơ vô cơ, phản ứng được mô tả theo hai phương trình:
    • Khi sử dụng amoni ($\text{NH}_4^+$): $$106\text{CO}2 + 16\text{NH}4^+ + \text{HPO}4^{2-} + 106\text{H}2\text{O} \rightarrow \text{C}{106}\text{H}{263}\text{O}{110}\text{N}{16}\text{P} + 106\text{O}_2 + 15\text{H}^+$$
    • Khi sử dụng nitrate ($\text{NO}_3^-$): $$106\text{CO}2 + 16\text{NO}3^- + \text{HPO}4^{2-} + 122\text{H}2\text{O} + 17\text{H}^+ \rightarrow \text{C}{106}\text{H}{263}\text{O}{110}\text{N}{16}\text{P} + 138\text{O}_2$$ Tỷ lệ nguyên tử chuẩn Redfield $\text{C}:\text{N}:\text{P}$ được lý tưởng hóa theo khối lượng khô là $40 : 7{,}2 : 1$.
  • Yếu tố giới hạn dinh dưỡng: Tỷ lệ $\text{TN}/\text{TP}$ quyết định chất dinh dưỡng giới hạn sự phát triển của tảo trong nước ngọt:
    • $\text{TN}/\text{TP} < 4{,}5$: Nitơ là yếu tố giới hạn.
    • $4{,}5 \le \text{TN}/\text{TP} \le 6$: Cả nitơ và phốt pho cùng là yếu tố giới hạn.
    • $\text{TN}/\text{TP} > 6$ (hoặc $> 7{,}2$ theo tỷ lệ khối lượng): Phốt pho là yếu tố giới hạn.
  • Phương pháp đánh giá mức độ phú dưỡng:
    • Phương pháp Hakanson et al. (2007): Dựa trên các ngưỡng giá trị nồng độ TP, TN và Chlorophyll-a (Chl.a).
    • Phương pháp chỉ số trạng thái dinh dưỡng Carlson (TSI) và bổ sung của Kratzer - Brezonik (1981): $$\text{TSI}(\text{TP}) = 14{,}42 \times \ln(\text{TP}) + 4{,}15 \quad (\text{TP tính theo } \mu\text{g/L})$$ $$\text{TSI}(\text{Chl.a}) = 9{,}81 \times \ln(\text{Chl.a}) + 30{,}6 \quad (\text{Chl.a tính theo } \mu\text{g/L})$$ $$\text{TSI}(\text{TN}) = 14{,}43 \times \ln(\text{TN}) + 54{,}45 \quad (\text{TN tính theo mg/L})$$ Phân loại: $\text{TSI} < 40$ (nghèo dinh dưỡng); $40 \le \text{TSI} \le 50$ (trung dưỡng); $50 < \text{TSI} \le 70$ (phú dưỡng); $\text{TSI} > 70$ (siêu phú dưỡng).
Thông số Nghèo dinh dưỡng Dinh dưỡng trung bình Phú dưỡng Siêu phú dưỡng
Tổng phốt pho - TP (mg/L) $< 0{,}008$ $0{,}008 - 0{,}025$ $0{,}025 - 0{,}06$ $> 0{,}06$
Tổng nitơ - TN (mg/L) $< 0{,}06$ $0{,}06 - 0{,}18$ $0{,}18 - 0{,}43$ $> 0{,}43$
Chlorophyll-a - Chl.a ($\mu\text{g/L}$) $< 2$ $2 - 6$ $6 - 20$ $> 20$

Phương pháp nghiên cứu và phát triển mô hình toán

  • Cấu trúc hệ phương trình vi phân mô hình: Mô hình toán gồm hệ phương trình vi phân thường (PTVPT - ODE) mô tả động thái tương tác của các biến trạng thái: sinh khối 3 nhóm thực vật nổi (tảo lục, vi khuẩn lam, tảo silic), sinh khối động vật phù du (ĐVPD), các hợp chất nitơ ($\text{NH}_4\text{-N}$, $\text{NO}_2\text{-N}$, $\text{NO}_3\text{-N}$, nitơ hữu cơ), phốt pho (phốt pho vô cơ hòa tan - DIP, phốt pho hữu cơ - OP), carbon hữu cơ (POC, DOC) và oxy hòa tan (DO).
  • Điểm cải tiến trong phương trình động học: Bổ sung các số hạng đại diện cho nồng độ và tải lượng chất dinh dưỡng gia nhập từ khí quyển (lắng đọng khô, ướt qua nước mưa) và dòng chảy tràn bề mặt đệm đô thị trực tiếp vào phương trình cân bằng nồng độ các dạng nitơ và phốt pho trong hồ tĩnh nông.
  • Phương pháp giải số: Sử dụng phương pháp Runge-Kutta kết hợp bậc 4 và bậc 5 (ODE45) để tích phân hệ PTVPT theo bước thời gian ngày. Toàn bộ thuật toán được lập trình bằng phần mềm Matlab.
  • Phương pháp phân tích độ nhạy và tối ưu hóa tham số:
    • Thực hiện phân tích độ nhạy thống kê nhằm xác định 15 tham số có mức độ ảnh hưởng lớn nhất trong tổng số 53 tham số mô hình.
    • Sử dụng thuật toán di truyền (Genetic Algorithm - GA) trong môi trường Matlab để tự động tìm kiếm bộ giá trị tối ưu toàn cục cho 53 tham số, triệt tiêu nhược điểm rơi vào điểm cực tiểu cục bộ của phương pháp thử sai truyền thống.
  • Tiêu chí đánh giá chất lượng mô hình: Sử dụng 4 chỉ số thống kê gồm: Sai số quân phương (RMSE), Hệ số hiệu quả mô hình Nash-Sutcliffe (NSE), Tỷ lệ sai số quân phương trên độ lệch chuẩn thực đo (RSR), và Phần trăm sai số (PBIAS).

Nội dung chính theo từng chương

Chương 1: Tổng quan vấn đề nghiên cứu

Chương 1 hệ thống hóa cơ sở lý luận về hiện tượng phú dưỡng, nguyên nhân, cơ chế và phương pháp đánh giá (Hakanson, Carlson). Tác giả phân tích lịch sử phát triển các thế hệ mô hình phú dưỡng trên thế giới (thực nghiệm, cân bằng vật chất, sinh thái) và tổng kết các nghiên cứu mô hình hóa chất lượng nước hồ tại Việt Nam.

Chương này cũng trình bày kết quả khảo sát thực tế hiện trạng chất lượng nước của hệ thống hồ nội thành Hà Nội. Dữ liệu quan trắc cho thấy nhiều hồ dù đã được kè bờ, cải tạo tách nước thải sinh hoạt như hồ Hoàn Kiếm, hồ Thủ Lệ, hồ Thành Công, hồ Hai Bà Trưng, hồ Thiền Quang, hồ Cần, hồ Nghĩa Tân và hồ Cự Chính vẫn ở trạng thái phú dưỡng đến siêu phú dưỡng. Điển hình, hồ Hoàn Kiếm có $\text{TSI}$ thuộc mức siêu phú dưỡng, vi khuẩn lam thuộc chi Microcystis chiếm tới 90% tổng số thực vật nổi, sản sinh độc tố gan microcystins từ 2,1–46 $\mu\text{g/L}$; hồ Thành Công có $\text{TN} = 2{,}2\text{ mg/L}$, $\text{TP} = 0{,}4\text{ mg/L}$, $\text{TSI} = 65{,}2$; hồ Thiền Quang có $\text{TN} = 2{,}47\text{ mg/L}$, $\text{TP} = 0{,}41\text{ mg/L}$, $\text{TSI} = 64{,}1$.

Về điều kiện tự nhiên khu vực nghiên cứu hồ Cự Chính (quận Thanh Xuân), chuỗi số liệu khí tượng nhiều năm (2002–2018) tại trạm Láng được thống kê chi tiết:

  • Nhiệt độ trung bình năm đạt $24{,}5^\circ\text{C}$; tháng 6 cao nhất ($30{,}1^\circ\text{C}$), tháng 1 thấp nhất ($16{,}5^\circ\text{C}$); từ tháng 5 đến tháng 9 nhiệt độ đều trên $28^\circ\text{C}$, rất thuận lợi cho tảo lục và vi khuẩn lam phát triển.
  • Tổng số giờ nắng trung bình năm là 1.292 giờ; tháng 6 cao nhất (156 giờ), tháng 3 thấp nhất (42 giờ).
  • Độ ẩm tương đối trung bình năm đạt $78%$.
  • Lượng mưa trung bình năm đạt 1.652 mm; mùa mưa tập trung từ tháng 5 đến tháng 9 (chiếm $80%$ lượng mưa cả năm), tháng 8 có lượng mưa cao nhất (348 mm), tạo nguồn rửa trôi dinh dưỡng lớn vào hồ.
Thông số khí tượng (Trạm Láng, 2002–2018) Th.1 Th.2 Th.3 Th.4 Th.5 Th.6 Th.7 Th.8 Th.9 Th.10 Th.11 Th.12 Trung bình / Tổng
Nhiệt độ trung bình ($^\circ\text{C}$) 16,5 18,7 20,9 24,9 28,3 30,1 29,8 29,0 28,2 26,4 22,9 18,7 24,5
Số giờ nắng trung bình (giờ) 59 53 42 81 147 156 152 139 130 122 118 93 1.292
Độ ẩm tương đối trung bình ($%$) 77 81 82 82 78 76 78 81 79 74 74 73 78
Lượng mưa trung bình (mm) 25 19 37 72 205 228 308 348 229 100 60 21 1.652

Chương 2: Phát triển mô hình toán mô phỏng quá trình phú dưỡng trong hồ tĩnh nông

Chương 2 trình bày chi tiết quy trình 8 bước phát triển mô hình toán. Tác giả xây dựng các phương trình vi phân mô tả các quá trình: quang hợp, hô hấp, bài tiết, tử vong và lắng đọng của 3 nhóm thực vật nổi (tảo lục, vi khuẩn lam, tảo silic) và động vật phù du; các quá trình khoáng hóa, thủy phân, nitrat hóa, khử nitrat, tiêu thụ và tái sinh oxy hòa tan; trao đổi chất dinh dưỡng giữa trầm tích và cột nước.

Mô hình được cải tiến bằng cách tích hợp nồng độ các chất dinh dưỡng trong nước mưa rơi trực tiếp trên mặt hồ và trong dòng chảy tràn từ khuôn viên bê tông 1.500 $\text{m}^2$ xung quanh hồ. Chương trình mô phỏng được lập trình hoàn chỉnh trên Matlab, sử dụng thuật toán Runge-Kutta bậc 4-5 giải hệ ODE và công cụ GA Tool để tối ưu hóa tham số. Các bước thu thập số liệu thực nghiệm, xác định điều kiện ban đầu, điều kiện biên ngoại sinh (nhiệt độ nước, cường độ bức xạ mặt trời, lượng mưa ngày, chất lượng nước mưa) và thiết lập kịch bản mô phỏng tại hồ Cự Chính được chuẩn hóa chi tiết.

+-----------------------------------------------------------------------------------+
|               QUY TRÌNH 8 BƯỚC PHÁT TRIỂN MÔ HÌNH PHÚ DƯỠNG LUẬN ÁN               |
+-----------------------------------------------------------------------------------+
  Bước 1: Xác định bài toán, mục tiêu, đối tượng và phạm vi nghiên cứu
  Bước 2: Kế thừa và cải tiến hệ phương trình vi phân mô phỏng quá trình phú dưỡng
  Bước 3: Giải hệ phương trình ODE bằng phương pháp Runge-Kutta 4-5 trên Matlab
  Bước 4: Quan trắc, đo đạc dữ liệu thực nghiệm và thu thập chuỗi số liệu khí tượng
  Bước 5: Phân tích độ nhạy thống kê (Xác định 15 tham số nhạy nhất)
  Bước 6: Hiệu chỉnh và tối ưu hóa 53 tham số mô hình bằng Thuật giải di truyền (GA)
  Bước 7: Kiểm định mô hình với tập dữ liệu độc lập (Đánh giá RMSE, NSE, RSR, PBIAS)
  Bước 8: Mô phỏng các kịch bản quản lý, kiểm soát phú dưỡng hồ

Chương 3: Kết quả nghiên cứu và thảo luận

Chương 3 trình bày các kết quả thực nghiệm và mô phỏng chính:

  1. Chất lượng nước và đặc trưng sinh học hồ Cự Chính: Kết quả phân tích chất lượng nước thực tế xác nhận hồ Cự Chính ở trạng thái siêu phú dưỡng theo thang Hakanson và Carlson TSI. Thành phần thực vật nổi bị chi phối áp đảo bởi vi khuẩn lam và tảo lục. Nồng độ $\text{DIN}$, $\text{DIP}$, nhiệt độ nước và bức xạ mặt trời có tương quan thuận rõ rệt với mật độ tế bào thực vật nổi.
  2. Phân tích độ nhạy và tối ưu hóa tham số: Phân tích độ nhạy đã trích xuất được 15 tham số có mức độ ảnh hưởng lớn nhất đối với các biến trạng thái. Quá trình tối ưu hóa bằng GA qua 38 lần chạy lặp đã tìm ra bộ 53 tham số tối ưu, giúp giá trị hàm mục tiêu đạt cực tiểu ổn định.
  3. Hiệu chỉnh và kiểm định mô hình: Kết quả so sánh giữa giá trị mô phỏng và số liệu thực đo cho thấy sự phù hợp cao đối với các biến trạng thái: sinh khối vi khuẩn lam, tảo lục, tảo silic, động vật phù du, oxy hòa tan và nồng độ DIP. Các chỉ số sai số (RMSE, NSE, RSR, PBIAS) trong cả hai giai đoạn hiệu chỉnh và kiểm định đều đạt tiêu chí thống kê từ tốt đến rất tốt.
  4. Mô phỏng các kịch bản kiểm soát: Luận án tính toán mô phỏng các kịch bản can thiệp thực tế nhằm đánh giá mức độ suy giảm sinh khối thực vật nổi, phục hồi oxy hòa tan và kiểm soát nồng độ DIP trung bình trong hồ Cự Chính, cung cấp cơ sở định lượng cho công tác quản lý môi trường hồ.

Kết quả và những đóng góp mới

Những đóng góp mới về lý luận và khoa học

  • Xây dựng và phát triển thành công mô hình toán sinh thái mô phỏng quá trình phú dưỡng phù hợp cho các vùng nước tĩnh nông ở điều kiện khí hậu nhiệt đới.
  • Cải tiến hệ phương trình động học chất dinh dưỡng thông qua việc bổ sung thành phần nồng độ và tải lượng chất dinh dưỡng gia nhập từ khí quyển (lắng đọng khô, ướt) và lượng nước mưa chảy tràn trên lưu vực bề mặt đệm đô thị.

Những đóng góp mới về thực tiễn

  • Ứng dụng thành công mô hình đã phát triển cho hồ Cự Chính (quận Thanh Xuân, Hà Nội), xác lập được bộ 53 tham số tối ưu (trong đó có 15 tham số nhạy nhất) được hiệu chỉnh tự động bằng thuật giải di truyền GA.
  • Cung cấp bộ công cụ mô phỏng hỗ trợ các nhà quản lý môi trường đô thị dự báo diễn biến chất lượng nước, sự bùng phát của tảo và đánh giá hiệu quả các giải pháp kiểm soát phú dưỡng hồ nông nội thành.

Kiến nghị và giải pháp đề xuất

  • Áp dụng các giải pháp kỹ thuật kiểm soát và thu gom triệt để nước mưa chảy tràn quanh khuôn viên hồ trước khi chảy vào hồ.
  • Thực hiện định kỳ các biện pháp nạo vét, xử lý lớp bùn đáy nhằm hạn chế sự khuếch tán ngược của phốt pho và nitơ từ trầm tích trong các tháng mùa nóng.
  • Tăng cường quan trắc liên tục các biến ngoại sinh (nhiệt độ nước, bức xạ mặt trời) để kịp thời đưa ra cảnh báo nguy cơ nở hoa của vi khuẩn lam và tảo độc.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của nghiên cứu

  • Mô hình phát triển dựa trên giả thiết hồ ở chế độ pha trộn hoàn toàn 1 lớp (0D/1D đơn giản hóa), chưa tích hợp đầy đủ mô hình thủy động lực học không gian 2 chiều (2D) hoặc 3 chiều (3D) đối với các hồ có hình thái phức tạp hoặc có dòng chảy lưu thông.
  • Số liệu kiểm định thực nghiệm mới chỉ thực hiện trên đối tượng hồ Cự Chính trong một chuỗi thời gian nhất định; chưa kiểm chứng trên nhiều loại hình hồ tĩnh nông khác nhau.
  • Mô hình chưa phân tách sâu các nhóm động vật phù du kích thước khác nhau cũng như động thái chi tiết của quần xã vi sinh vật đáy.

Hướng nghiên cứu tiếp theo

  • Tiếp tục mở rộng và tích hợp mô hình sinh thái với các mô hình thủy động lực học 2D/3D để áp dụng cho các hồ có diện tích lớn hơn, dòng chảy phức tạp.
  • Ứng dụng và kiểm định bộ tham số mô hình cho nhiều hồ nông đô thị khác tại Hà Nội và các tỉnh thành nhiệt đới.
  • Nghiên cứu bổ sung tác động dài hạn của biến đổi khí hậu và các hiện tượng thời tiết cực đoan đến quá trình phú dưỡng hồ đô thị.

Giá trị tham khảo

  • Đối với nghiên cứu sinh, học viên cao học và giảng viên: Là tài liệu tham khảo chuyên sâu trong các ngành Môi trường Đất và Nước, Kỹ thuật Môi trường, Quản lý Tài nguyên Nước và Tin học Môi trường. Phần đáng tham khảo nhất là quy trình 8 bước thiết lập hệ phương trình vi phân mô phỏng sinh thái hồ và kỹ thuật tối ưu hóa tham số mô hình bằng thuật giải di truyền (GA) trên Matlab.
  • Đối với cơ quan quản lý và hoạch định chính sách: Cung cấp cơ sở khoa học cho Sở Tài nguyên và Môi trường Hà Nội, Công ty Thoát nước Hà Nội và các đơn vị quản lý công viên cây xanh - hồ nước trong việc đánh giá tải lượng ô nhiễm phân tán, dự báo sự nở hoa của vi khuẩn lam và lựa chọn phương án xử lý phú dưỡng hồ đô thị.

Câu hỏi thường gặp

1. Tại sao nhiều hồ nội đô Hà Nội đã tách hoàn toàn nước thải sinh hoạt nhưng vẫn duy trì trạng thái siêu phú dưỡng?

Theo kết quả khảo sát và phân tích của luận án tại các hồ như Hoàn Kiếm, Thành Công, Thiền Quang và Cự Chính, nguyên nhân là do hồ nông vẫn liên tục nhận các nguồn chất dinh dưỡng phân tán từ bên ngoài gồm lắng đọng khí quyển (nước mưa có nồng độ TN, TP cao), dòng chảy tràn rửa trôi bề mặt đệm đô thị bị bê tông hóa, cùng với nguồn dinh dưỡng nội tại giải phóng từ quá trình phân hủy xác thực vật thủy sinh và sự khuếch tán chất dinh dưỡng từ trầm tích bùn đáy dưới điều kiện nhiệt độ cao.

2. Điểm cải tiến kỹ thuật chính trong các phương trình toán học của mô hình luận án là gì?

Điểm cải tiến cốt lõi là việc đưa thêm các số hạng nồng độ chất dinh dưỡng từ khí quyển (lắng đọng ướt trong nước mưa) và dòng chảy tràn bề mặt đệm đô thị trực tiếp vào phương trình vi phân mô phỏng cân bằng động thái các dạng nitơ ($\text{NH}_4\text{-N}$, $\text{NO}_2\text{-N}$, $\text{NO}_3\text{-N}$, $\text{ON}$) và phốt pho ($\text{DIP}$, $\text{OP}$), khắc phục hạn chế bỏ qua nguồn ngoại lai phân tán của các mô hình hồ tĩnh trước đây.

3. Tác giả đã giải hệ phương trình vi phân và tối ưu hóa bộ tham số mô hình bằng phương pháp nào?

Hệ phương trình vi phân thường (ODE) được tích phân số bằng phương pháp Runge-Kutta bậc 4-5 (ODE45) trên phần mềm Matlab. Bộ 53 tham số của mô hình được tự động hiệu chỉnh và tối ưu hóa toàn cục thông qua Thuật giải di truyền (Genetic Algorithm - GA) với công cụ GA Tool trong Matlab, kết hợp phân tích độ nhạy thống kê để nhận diện 15 tham số có mức độ ảnh hưởng lớn nhất.

4. Các chỉ số nào được sử dụng để đánh giá trạng thái phú dưỡng và độ chính xác của mô hình?

  • Đánh giá trạng thái phú dưỡng: Sử dụng phương pháp phân loại của Hakanson et al. (2007) và Chỉ số trạng thái dinh dưỡng Carlson ($\text{TSI}$ đối với TP, Chl.a) kết hợp công thức Kratzer - Brezonik ($\text{TSI}$ đối với TN).
  • Đánh giá sai số mô hình: Sử dụng 4 chỉ số thống kê gồm Sai số quân phương (RMSE), Hệ số hiệu quả mô hình Nash-Sutcliffe (NSE), Tỷ lệ sai số quân phương trên độ lệch chuẩn thực đo (RSR) và Phần trăm sai số (PBIAS).

5. Hồ Cự Chính có những đặc điểm hình thái nào phù hợp với các giả thiết của mô hình phát triển?

Hồ Cự Chính có diện tích nhỏ (khoảng 4.000 $\text{m}^2$), độ sâu trung bình nông (1,5–1,7 m), dung tích khoảng 4.600 $\text{m}^3$, có bờ kè chắc chắn không có dòng ngầm lớn, không có cống xả nước thải sinh hoạt trực tiếp, nguồn nước bổ cập chủ yếu là nước mưa rơi trên mặt hồ và chảy tràn từ khuôn viên bê tông 1.500 $\text{m}^2$. Các đặc điểm này hoàn toàn phù hợp với giả thiết vùng nước tĩnh nông, pha trộn đều và chịu tác động chính bởi nguồn dinh dưỡng phân tán.

Kết luận

Luận án của tác giả Tạ Đăng Thuần đã phát triển thành công mô hình toán sinh thái mô phỏng quá trình phú dưỡng cho các thủy vực nước tĩnh nông nhiệt đới, tích hợp nguồn bổ cấp dinh dưỡng từ nước mưa và dòng chảy tràn đô thị. Bằng việc kết hợp phương pháp số Runge-Kutta bậc 4-5 và thuật giải di truyền GA trên nền tảng Matlab, mô hình đã giải quyết hiệu quả bài toán hiệu chỉnh 53 tham số phức tạp. Kết quả ứng dụng tại hồ Cự Chính (Hà Nội) cho thấy độ tin cậy cao qua các chỉ số kiểm định thống kê và phản ánh chính xác tương tác giữa các nhóm thực vật nổi với điều kiện dinh dưỡng và khí tượng. Công trình đóng góp cơ sở khoa học và công cụ thực tiễn quan trọng phục vụ công tác kiểm soát, quản lý và phục hồi chất lượng nước các hồ đô thị.