Luận án tiến sĩ kỹ thuật về mô hình toán mô phỏng phú dưỡng hồ Cự Chính - Hà Nội của Tạ Đăng Thuần

Phú dưỡng là hiện tượng tích tụ quá mức chất dinh dưỡng (nitơ, phốt pho) trong môi trường nước. Điều này kích thích sự phát triển bùng nổ của thực vật thủy sinh, đặc biệt là tảo. Quá trình này làm suy giảm chất lượng nước nghiêm trọng, ảnh hưởng đến đa dạng sinh học. Nguyên nhân chính bao gồm xả thải công nghiệp, nông nghiệp và sinh hoạt chưa qua xử lý vào các thủy vực. Nước thải chứa hàm lượng dinh dưỡng cao, gây mất cân bằng hệ sinh thái. Hiện tượng phú dưỡng dẫn đến giảm oxy hòa tan, cá chết hàng loạt và phát sinh mùi hôi. Các hồ tĩnh nông đặc biệt dễ bị phú dưỡng do khả năng trao đổi nước hạn chế. Mức độ phú dưỡng cần được đánh giá định kỳ. Việc hiểu rõ cơ chế hình thành phú dưỡng là cần thiết để xây dựng giải pháp bền vững. Các chỉ số như nồng độ chlorophyll-a, tổng phốt pho, tổng nitơ được sử dụng để định lượng mức độ phú dưỡng. Tình trạng này gây ra thách thức lớn trong quản lý môi trường nước đô thị.

Nghiên cứu mô hình phú dưỡng trên thế giới và Việt Nam có lịch sử phát triển lâu dài. Các mô hình toán học là công cụ hiệu quả để dự báo và quản lý hiện tượng này. Trên thế giới, nhiều mô hình phức tạp đã được phát triển, tích hợp các quá trình vật lý, hóa học và sinh học. Các mô hình này giúp dự đoán sự biến đổi chất lượng nước theo thời gian và không gian. Ở Việt Nam, các nghiên cứu tập trung vào việc áp dụng và điều chỉnh các mô hình hiện có cho điều kiện cụ thể của các hồ. Nhiều công trình khoa học đã sử dụng mô hình để đánh giá mức độ phú dưỡng. Tuy nhiên, việc xây dựng một mô hình chuyên biệt cho hồ tĩnh nông, có tính đến các đặc thù của địa phương, vẫn còn nhiều thách thức. Các mô hình cần được hiệu chỉnh và kiểm định kỹ lưỡng với dữ liệu thực nghiệm. Sự phát triển mô hình cần một cách tiếp cận đa ngành.

Hồ Cự Chính là một điển hình cho hồ tĩnh nông chịu ảnh hưởng phú dưỡng tại Hà Nội. Khu vực này đối mặt với áp lực đô thị hóa mạnh mẽ. Nước thải sinh hoạt từ các khu dân cư xung quanh đổ trực tiếp vào hồ, mang theo lượng lớn chất dinh dưỡng. Điều này gây ô nhiễm nghiêm trọng, làm suy giảm chất lượng nước. Hồ Cự Chính thường xuyên xuất hiện hiện tượng tảo nở hoa, khiến nước có màu xanh đậm và bốc mùi khó chịu. Mức độ phú dưỡng tại hồ Cự Chính đạt ngưỡng cao, ảnh hưởng đến cảnh quan, môi trường sống và sức khỏe cộng đồng. Việc nghiên cứu hồ Cự Chính cung cấp dữ liệu thực tế quý giá cho việc phát triển và kiểm định mô hình. Phân tích thực trạng hồ giúp xác định các yếu tố đầu vào quan trọng cho mô hình toán. Các giải pháp được phát triển cho hồ Cự Chính có tiềm năng áp dụng cho các hồ tĩnh nông tương tự trên cả nước, góp phần vào quản lý môi trường bền vững.

Mô hình toán mô phỏng phú dưỡng dựa trên cơ sở lý thuyết về động học chất dinh dưỡng và tăng trưởng sinh vật. Các yếu tố chính ảnh hưởng đến quá trình phú dưỡng bao gồm nhiệt độ nước, ánh sáng mặt trời, nồng độ nitơ và phốt pho. Tốc độ tăng trưởng của tảo phụ thuộc vào sự có mặt của các chất dinh dưỡng giới hạn. Nồng độ oxy hòa tan cũng là một chỉ thị quan trọng của chất lượng nước. Các giả thiết trong mô hình đơn giản hóa các quá trình phức tạp trong hồ. Ví dụ, hồ được coi là hệ thống khuấy trộn hoàn toàn hoặc chia thành các ngăn riêng biệt. Các quá trình sinh học như quang hợp, hô hấp, lắng đọng và phân hủy chất hữu cơ đều được xem xét. Hiểu rõ các yếu tố này giúp xây dựng một mô hình chính xác và hiệu quả. Mô hình cần phản ánh chân thực các phản ứng của hệ sinh thái hồ.

Mô hình sử dụng các phương trình vi phân để mô tả sự thay đổi nồng độ các chất trong hồ theo thời gian. Các phương trình này tích hợp động học của tảo, nitơ vô cơ, phốt pho vô cơ và oxy hòa tan. Phương trình cho tảo thường bao gồm các thành phần tăng trưởng (quang hợp), suy giảm (hô hấp, lắng, tử vong do tự nhiên). Phương trình dinh dưỡng tính toán lượng chất dinh dưỡng tiêu thụ bởi tảo và lượng chất dinh dưỡng tái tạo từ phân hủy. Mô hình cải tiến có thể bổ sung các yếu tố như dòng chảy, trao đổi với trầm tích. Việc lựa chọn các hàm giới hạn (Limiting functions) cho tăng trưởng của tảo là rất quan trọng. Các hàm này thường là hàm Michaelis-Menten hoặc các biến thể của chúng. Mục tiêu là phát triển một hệ phương trình mô tả đầy đủ các quá trình sinh hóa trong hồ.

Quá trình phát triển mô hình bao gồm nhiều bước tuần tự. Bước đầu tiên là xác định mục tiêu và phạm vi của mô hình. Tiếp theo, xây dựng sơ đồ khái niệm và xác định các biến trạng thái, các tham số. Sau đó, lựa chọn các phương trình toán học phù hợp để mô tả các quá trình. Giai đoạn lập trình mô hình thường sử dụng các phần mềm chuyên dụng như Matlab. Dữ liệu thực nghiệm là cần thiết để hiệu chỉnh và kiểm định mô hình. Hiệu chỉnh mô hình điều chỉnh các tham số để kết quả mô phỏng phù hợp với dữ liệu thực tế. Kiểm định mô hình đánh giá khả năng dự báo của mô hình với dữ liệu độc lập. Cuối cùng, mô hình được sử dụng để chạy các kịch bản khác nhau. Quá trình này đảm bảo tính tin cậy và khả năng ứng dụng của mô hình.

Các phương trình vi phân trong mô hình phú dưỡng thường không có lời giải giải tích trực tiếp. Do đó, cần sử dụng các phương pháp số để tìm lời giải xấp xỉ. Phương pháp Runge-Kutta là một trong những phương pháp phổ biến và hiệu quả, đặc biệt cho việc giải hệ phương trình vi phân thường. Phương pháp này cung cấp xấp xỉ chính xác cho lời giải, giảm thiểu sai số tính toán. Việc lựa chọn bước thời gian phù hợp là rất quan trọng để đảm bảo độ ổn định và chính xác của lời giải số. Các công cụ phần mềm như Matlab tích hợp sẵn các hàm giải phương trình vi phân bằng phương pháp Runge-Kutta, giúp đơn giản hóa quá trình tính toán. Áp dụng phương pháp số giúp mô hình có thể xử lý các hệ thống phức tạp với nhiều biến phụ thuộc và tương tác. Kết quả giải phương trình cho biết sự thay đổi nồng độ các chất trong hồ theo thời gian một cách chi tiết.

Để hiệu chỉnh và kiểm định mô hình, cần có dữ liệu thực nghiệm đáng tin cậy. Quy trình lấy mẫu, đo đạc chất lượng nước tuân thủ các tiêu chuẩn kỹ thuật. Các thông số cần đo bao gồm nhiệt độ nước, pH, oxy hòa tan, độ đục. Đặc biệt, nồng độ chlorophyll-a, tổng nitơ, tổng phốt pho là các chỉ số quan trọng cho đánh giá phú dưỡng. Dữ liệu biến ngoại sinh như cường độ bức xạ mặt trời, lượng mưa cũng cần được thu thập. Việc lấy mẫu được thực hiện định kỳ trong suốt thời gian nghiên cứu tại các điểm khác nhau trong hồ. Phân tích mẫu nước được tiến hành tại các phòng thí nghiệm chuyên sâu, đảm bảo độ chính xác cao. Dữ liệu thực nghiệm này cung cấp cơ sở đối chiếu vững chắc cho kết quả mô phỏng. Độ chính xác của dữ liệu thực nghiệm ảnh hưởng trực tiếp đến độ tin cậy của mô hình.

Phần mềm Matlab đóng vai trò quan trọng trong việc phát triển và vận hành mô hình. Matlab cung cấp môi trường lập trình linh hoạt và các công cụ toán học mạnh mẽ. Mô hình được lập trình dưới dạng các script và hàm trong Matlab, cho phép tùy chỉnh cao. Giai đoạn hiệu chỉnh mô hình (calibration) sử dụng dữ liệu thực nghiệm để tối ưu các tham số của mô hình. Các thuật toán tối ưu hóa có thể được áp dụng để tìm ra bộ tham số phù hợp nhất, giảm thiểu sai số giữa mô phỏng và thực tế. Sau khi hiệu chỉnh, mô hình được kiểm định (validation) bằng một bộ dữ liệu độc lập, chưa từng được sử dụng. Việc kiểm định đánh giá khả năng dự báo của mô hình trên các điều kiện chưa từng thấy. Kết quả kiểm định được đánh giá bằng các chỉ số thống kê như sai số trung bình bình phương (RMSE) hoặc hệ số tương quan. Giao diện công cụ GA trong Matlab cũng hỗ trợ quá trình này.

Kết quả phân tích chất lượng nước tại hồ Cự Chính cho thấy tình trạng phú dưỡng nghiêm trọng. Nồng độ các chất dinh dưỡng (nitơ, phốt pho) vượt ngưỡng cho phép nhiều lần. Chỉ số chlorophyll-a cao, cho thấy sự phát triển mạnh mẽ của tảo. Hàm lượng oxy hòa tan thường giảm mạnh vào ban đêm, gây áp lực lên các loài thủy sinh. Quần xã thực vật nổi tại hồ thường xuyên bị chi phối bởi các loài tảo lam, dẫn đến hiện tượng tảo nở hoa. Tình trạng này phản ánh rõ nét tác động của hoạt động nhân sinh và đô thị hóa. Dữ liệu thực tế cung cấp bức tranh tổng thể về hiện trạng môi trường hồ. Mức độ phú dưỡng được đánh giá dựa trên các tiêu chuẩn và chỉ số phú dưỡng quốc tế, khẳng định sự cấp thiết của các biện pháp can thiệp. Các mối liên hệ giữa các thông số chất lượng nước được phân tích kỹ lưỡng.

Phân tích độ nhạy xác định mức độ ảnh hưởng của từng tham số mô hình đến kết quả mô phỏng. Các tham số có độ nhạy cao cần được xác định chính xác hơn, tập trung vào chúng trong quá trình hiệu chỉnh. Kết quả phân tích độ nhạy giúp hiểu rõ hơn về sự tương tác giữa các yếu tố trong mô hình. Sau đó, quá trình hiệu chỉnh mô hình được thực hiện một cách tỉ mỉ. Các tham số của mô hình được điều chỉnh để kết quả mô phỏng phù hợp nhất với dữ liệu thực tế đã thu thập. Các giá trị tham số được tối ưu hóa trong một khoảng giá trị hợp lý, dựa trên các nghiên cứu trước đây. Biểu đồ so sánh kết quả mô phỏng và đo đạc thực tế được sử dụng để đánh giá sự phù hợp của mô hình. Quá trình hiệu chỉnh đảm bảo mô hình phản ánh đúng động thái của hệ thống hồ, tăng cường độ tin cậy. Độ chính xác của mô hình sau hiệu chỉnh đạt yêu cầu về thống kê, sẵn sàng cho các phân tích tiếp theo.

Mô hình được sử dụng để mô phỏng sự phát triển của tảo và chất lượng nước dưới các kịch bản khác nhau. Các kịch bản này bao gồm việc giảm tải lượng chất dinh dưỡng đầu vào, thay đổi chế độ thủy văn của hồ, hoặc ảnh hưởng của biến đổi khí hậu. Ví dụ, mô phỏng ảnh hưởng của việc giảm nồng độ phốt pho đầu vào từ 20% đến 50%. Kết quả mô phỏng cho thấy tác động của từng kịch bản lên mức độ phú dưỡng và nồng độ chlorophyll-a. Mô hình có khả năng dự báo xu hướng thay đổi chất lượng nước trong tương lai, cung cấp thông tin quý giá. Điều này hỗ trợ các nhà quản lý đưa ra quyết định dựa trên bằng chứng khoa học, lựa chọn các biện pháp can thiệp hiệu quả nhất. Các kịch bản giúp đánh giá hiệu quả của các biện pháp xử lý và phòng ngừa trước khi triển khai thực tế. Sự phát triển của tảo được mô phỏng chi tiết dưới các hàm giới hạn khác nhau, cho thấy phản ứng của hệ sinh thái.

Mô hình toán mô phỏng phú dưỡng này có những đóng góp mới đáng kể. Mô hình được phát triển chuyên biệt cho hồ tĩnh nông, điều chỉnh phù hợp với điều kiện khí hậu và đặc điểm thủy văn của Việt Nam. Cấu trúc mô hình tích hợp sâu rộng các quá trình vật lý, hóa học, sinh học ảnh hưởng đến phú dưỡng một cách toàn diện. Mô hình cho phép dự báo chính xác hơn về sự biến động của chất lượng nước và nồng độ tảo trong hồ. Khả năng mô phỏng các kịch bản khác nhau giúp cung cấp cơ sở khoa học vững chắc cho các quyết định quản lý. Việc sử dụng công cụ Matlab giúp mô hình dễ dàng triển khai, điều chỉnh và mở rộng. Các tham số mô hình được hiệu chỉnh dựa trên dữ liệu thực tế của hồ Cự Chính, đảm bảo tính ứng dụng cao. Mô hình cung cấp một công cụ mạnh mẽ để đánh giá, kiểm soát và dự báo hiện tượng phú dưỡng, hỗ trợ công tác bảo vệ môi trường.

Mặc dù mô hình đạt được những kết quả khả quan và có nhiều đóng góp, vẫn còn một số hạn chế cần được xem xét. Mô hình có thể chưa bao gồm tất cả các yếu tố phức tạp của hệ sinh thái hồ, ví dụ như tác động của trầm tích đến việc giải phóng dinh dưỡng hoặc ảnh hưởng của chuỗi thức ăn bậc cao hơn. Việc thu thập dữ liệu thực nghiệm đầy đủ và liên tục với tần suất cao vẫn là một thách thức lớn về thời gian và nguồn lực. Định hướng nghiên cứu tiếp theo có thể mở rộng mô hình đa chiều (2D hoặc 3D) để mô tả chi tiết hơn sự phân bố chất lượng nước. Phát triển mô hình tích hợp khả năng dự báo biến đổi khí hậu và các yếu tố thời tiết cực đoan. Nghiên cứu sâu hơn về tác động của các loài tảo cụ thể và sự tương tác giữa chúng. Tích hợp mô hình kinh tế-xã hội để đánh giá toàn diện hiệu quả và chi phí của các giải pháp quản lý. Hợp tác với các cơ quan quản lý để triển khai ứng dụng mô hình rộng rãi hơn trong thực tế.

Dựa trên kết quả mô hình và đánh giá phú dưỡng, nhiều kiến nghị được đưa ra nhằm cải thiện chất lượng nước hồ. Cần kiểm soát chặt chẽ nguồn thải chất dinh dưỡng (nitơ, phốt pho) từ các hoạt động sinh hoạt, nông nghiệp và công nghiệp vào hồ. Xây dựng và nâng cấp hệ thống thu gom, xử lý nước thải sinh hoạt cho khu vực xung quanh hồ Cự Chính. Áp dụng các biện pháp sinh học hoặc vật lý để kiểm soát sự phát triển của tảo, như thu hoạch tảo hoặc sử dụng các loài thủy sinh ăn tảo. Nạo vét bùn đáy định kỳ có thể giảm tải lượng dinh dưỡng nội sinh tích tụ trong trầm tích. Nâng cao nhận thức cộng đồng về bảo vệ môi trường hồ, khuyến khích các hành vi thân thiện với môi trường. Thiết lập hệ thống giám sát chất lượng nước tự động để theo dõi liên tục các thông số quan trọng. Các giải pháp cần được triển khai đồng bộ và bền vững, có sự phối hợp giữa các bên liên quan. Việc sử dụng mô hình dự báo sẽ hỗ trợ tối ưu hóa các chiến lược quản lý hồ, đảm bảo hiệu quả lâu dài. Bảo vệ môi trường hồ tĩnh nông là một nhiệm vụ cấp thiết cho sự phát triển bền vững.