Tổng quan về luận án

Nghiên cứu về động thái oxi hòa tan (Dissolved Oxygen – DO) đóng vai trò nền tảng trong hóa học môi trường nước và quản lý chất lượng thủy vực. Luận án tiến sĩ chuyên ngành Hóa lý thuyết và Hóa lý (Mã số: 62.19) của nghiên cứu sinh Lê Minh Thành với đề tài "Mô hình hóa sự thay đổi nồng độ oxi trong môi trường nước dưới tác động của lớp bùn đáy" (thực hiện tại Viện Hóa học – Viện Hàn lâm Khoa học và Công nghệ Việt Nam dưới sự hướng dẫn của GS. Lê Quốc Hùng và TS. Phạm Hồng Phong, 2015) là công trình tiên phong tại Việt Nam tích hợp cơ chế hóa lý, quá trình khuếch tán phân tử và phương pháp mô phỏng số giải hệ phương trình đạo hàm riêng nhằm làm sáng tỏ quy luật phân bố DO trong không gian và thời gian.

Khoảng trống nghiên cứu (research gap) trọng tâm mà luận án giải quyết xuất phát từ hạn chế của các mô hình chất lượng nước truyền thống: phần lớn các nghiên cứu trước đây chỉ áp dụng phương trình thực nghiệm một chiều dạng hộp đen hoặc mô hình kinh nghiệm của Streeter – Phelps (1925) vốn gộp chung các động thái phức tạp ở đáy thành một thông số hằng số tiêu thụ đáy cố định. Thiếu vắng các mô hình động học vi mô nhiều chiều mô tả chi tiết tương tác khuếch tán – phản ứng tại ranh giới pha bùn – nước và cơ chế phân rã chất hữu cơ bậc cao. Để giải quyết khoảng trống này, luận án đặt ra ba câu hỏi nghiên cứu cốt lõi:

  1. Quy luật động học nào mô tả chính xác nhất tốc độ tiêu thụ DO bởi chất hữu cơ và nhu cầu oxi bùn đáy (Sediment Oxygen Demand – SOD) trong điều kiện ranh giới pha bùn – nước?
  2. Bùn đáy tác động như thế nào đến sự phân tầng nồng độ DO theo độ sâu và thời gian trong các cột nước tĩnh thông qua cơ chế khuếch tán Fickian?
  3. Làm thế nào để thiết lập một mô hình toán học giải tích – số hoàn chỉnh trên nền tảng phương pháp phần tử hữu hạn (Finite Element Method – FEM) có khả năng dự báo chuẩn xác phân bố DO tại các thủy vực tự nhiên?

Hệ thống giả thuyết nghiên cứu được kiểm định bao gồm:

  • Giả thuyết $H_1$: Sự suy giảm nồng độ DO tại lớp nước sát đáy chịu sự chi phối quyết định bởi quá trình tiêu thụ hóa lý và sinh hóa của lớp bùn đáy với tốc độ phản ứng tỷ lệ phi tuyến với gradient nồng độ tại lớp biên khuếch tán.
  • Giả thuyết $H_2$: Hệ phương trình lan truyền – khuếch tán kết hợp điều kiện biên hỗn hợp (Dirichlet, Neumann, Robin) trên phần mềm Comsol Multiphysics phản ánh sát thực nghiệm diễn biến nồng độ DO hơn các mô hình giải tích 1D cổ điển.
  • Giả thuyết $H_3$: Trong điều kiện nước tĩnh, tốc độ khuếch tán phân tử của oxi trong nước ($D_1 \approx 5 \times 10^{-5}\text{ cm}^2/\text{s}$) là yếu tố hạn chế vận tốc (rate-limiting step), khiến quá trình tái cấp khí bề mặt không thể bù đắp lượng DO tiêu hao ở tầng đáy nếu không có xáo trộn thủy lực.

Khung lý thuyết của luận án được xây dựng dựa trên Định luật khuếch tán Fick II, Định luật hòa tan chất khí Henry, động học phản ứng oxy hóa – sinh hóa bậc tổng quát và lý thuyết truyền khối qua màng biên. Phạm vi thực nghiệm của công trình bao gồm hệ mô hình vật lý cột nước kiểm soát trong phòng thí nghiệm (chiều cao cột nước 0,5 m; 1,0 m; 2,0 m) và khảo sát thực địa đa tầng tại 4 hồ tự nhiên: Hồ Bảy Mẫu, Hồ Gươm, Hồ Tây và Hồ Thủy điện Hòa Bình.

Literature Review và Positioning

Lịch sử nghiên cứu chất lượng nước ghi nhận mô hình kinh điển Streeter – Phelps (1925) đặt nền móng cho việc mô tả đường cong suy giảm DO (oxygen sag curve) trên sông dựa trên sự cân bằng giữa tốc độ tái sục khí từ khí quyển ($k_2$) và tốc độ tiêu thụ oxi bởi nhu cầu sinh hóa BOD ($k_1$). Tuy nhiên, trường phái này bộc lộ nhược điểm lớn khi coi hệ thống là đồng nhất 1 chiều dọc dòng chảy và giả định động học phản ứng phân hủy hữu cơ luôn tuân theo bậc 1 ($-\frac{dC_{BOD}}{dt} = k_1 C_{BOD}$).

Các tranh luận khoa học nảy sinh mạnh mẽ giữa hai trường phái:

  • Trường phái động học cổ điển (Fair, 1939; Thomann & Mueller, 1987) duy trì quan điểm tuyến tính hóa để đơn giản hóa việc giải phương trình vi phân.
  • Ngược lại, trường phái động học phi tuyến (Fillos & Molof, 1972; Borsuk et al., 2001; Di Toro, 2001) chứng minh rằng quá trình tiêu thụ oxi tại lớp bùn đáy phụ thuộc chặt chẽ vào nồng độ DO cục bộ theo hàm bán bão hòa Michaelis – Menten ($\frac{k_3 C_{DO}}{k_{DOhalf} + C_{DO}}$), đồng thời phản ứng oxy hóa các hợp chất hữu cơ phức tạp trong trầm tích diễn ra theo động học bậc hai hoặc bậc tổng quát $m$ ($-\frac{dC_{DO}}{dt} = k_1 C_{BOD}^m$).

Khi so sánh với các công trình quốc tế điển hình như mô hình tích hợp QUAL2K của Cơ quan Bảo vệ Môi trường Hoa Kỳ (Chapra et al., 2008) hay mô hình thông lượng trầm tích của Di Toro (2001), các công trình này thường đòi hỏi bộ thông số thủy động lực học khổng lồ và áp dụng cho các thủy vực chuyển động quy mô lớn. Luận án của Lê Minh Thành đã định vị chính xác điểm nghẽn học thuật: mô hình hóa cơ chế vi mô của tương tác ranh giới pha bùn – nước trong môi trường nước ngọt tĩnh lặng, nơi quá trình khuếch tán phân tử chiếm ưu thế tuyệt đối. Nghiên cứu tạo ra bước tiến quan trọng khi chuyển dịch từ mô hình hóa thủy lực vĩ mô sang mô hình hóa hóa lý vi mô chi tiết, lượng hóa chính xác dòng thông lượng oxi đi vào trầm tích đáy dựa trên các đặc tính cấu trúc vật lý của bùn.

Đó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 lý thuyết cân bằng oxi thủy vực thông qua 3 đóng góp lý thuyết căn bản:

  1. Tổng quát hóa phương trình động học tiêu thụ DO: Phát triển phương trình động học tiêu thụ oxi bởi chất hữu cơ từ dạng bậc một kinh điển sang dạng phương trình vi phân đạo hàm riêng liên kết bậc tổng quát, tích hợp đồng thời quá trình lắng đọng ($k_{set}$) và xáo trộn tái phát tán hữu cơ từ bùn vào nước ($k_{rese}$): $$\frac{\partial C_{BOD}}{\partial t} = D_2 \nabla^2 C_{BOD} - k_1 C_{BOD} - k_{set} C_{BOD} + k_{rese} C_{BODsed}$$
  2. Xác lập mô hình tương tác ranh giới pha bùn – nước đa chiều: Luận án chỉ rõ lớp bùn đáy không chỉ là nguồn tiêu thụ (sink) đơn thuần mà là một môi trường xốp liên tục có hệ số khuếch tán hiệu dụng riêng biệt ($D_1^*$). Nhu cầu oxi bùn đáy (SOD) được mô tả chuẩn xác dưới dạng thông lượng qua biên phụ thuộc vào độ xốp ($\varepsilon = 30-70%$), mật độ hạt ($\rho_s$) và tỉ khối bùn khô ($\rho_b$).
  3. Chuyển đổi hệ hình mô hình hóa (Paradigm shift): Chuyển dịch từ cách tiếp cận thông số gộp (lumped-parameter models) sang mô hình trường liên tục phân bố không gian (spatially-distributed continuous field models), chứng minh rằng gradient DO sát mặt bùn (trong phạm vi $0-5\text{ cm}$) dốc hơn nhiều so với các dự báo trước đây của lý thuyết màng tĩnh mỏng.

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

Khung phân tích của luận án tích hợp đồng thời 3 hệ lý thuyết: Lý thuyết lan truyền – khuếch tán Fick, Động học hóa học dị thể và Động học vận chuyển màng biên. Kiến trúc mô hình hóa được xây dựng theo chuỗi 4 mô hình tiến triển:

[Mô hình 1: Khảo sát lý thuyết lan truyền DO - BOD]
                     │
                     ▼
[Mô hình 2: Tiêu thụ DO do ảnh hưởng chủ đạo bởi Bùn đáy]
                     │
                     ▼
[Mô hình 3: Động thái vi mô tại Ranh giới pha Bùn - Nước]
                     │
                     ▼
[Mô hình 4: Mô hình tổng hợp đa yếu tố theo Độ sâu Thủy vực]

Điều kiện biên của hệ thống được toán học hóa nghiêm ngặt trên miền không gian $\Omega$ có biên $\partial\Omega$:

  • Biên mặt nước ($z = 0$): Thiết lập điều kiện biên Neumann mô tả dòng trao đổi khí với khí quyển tuân theo định luật Henry: $$-D_1 \frac{\partial C_{DO}}{\partial z}\Big|{z=0} = k_2 (C{DO}^{sat} - C_{DO})$$
  • Biên ranh giới bùn – nước ($z = -h$): Thiết lập điều kiện thông lượng tiêu thụ DO do nhu cầu SOD: $$-D_1 \frac{\partial C_{DO}}{\partial z}\Big|_{z=-h} = -k_3 \theta_s^{(T-20)}$$
  • Biên thành bên: Áp dụng điều kiện biên Neumann cách ly (no-flux boundary condition, $\mathbf{n} \cdot \nabla C = 0$).

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

Thiết kế nghiên cứu

Nghiên cứu vận hành theo hệ hình thực chứng (positivism) kết hợp chặt chẽ giữa suy diễn toán học lý thuyết và kiểm chứng thực nghiệm định lượng. Thiết kế nghiên cứu đa tầng bậc bao gồm:

  • Cấp độ vi mô: Đo đạc nồng độ DO ở cự ly $1,0\text{ cm}$ sát bề mặt tiếp xúc bùn – nước.
  • Cấp độ trung mô (meso-scale): Ba hệ thống cột nước thực nghiệm bằng vật liệu cách nhiệt có chiều cao $0,5\text{ m}$; $1,0\text{ m}$; và $2,0\text{ m}$, đường kính tiêu chuẩn, đặt trong điều kiện nhiệt độ phòng kiểm soát nghiêm ngặt ($20 \pm 1^\circ\text{C}$).
  • Cấp độ vĩ mô: Đo đạc trắc diện thẳng đứng nồng độ DO tại 4 hồ tự nhiên có đặc tính phú dưỡng và thủy văn khác biệt.

Quy trình nghiên cứu rigorous

Chiến lược thu mẫu bùn và nước được thực hiện theo tiêu chuẩn phân loại trầm tích Wentworth (kích thước hạt bùn $3,9 - 62,5\text{ }\mu\text{m}$). Quy trình phân tích hóa lý xác định: hàm lượng chất hữu cơ (OM: $4,13%$), tổng nitơ ($0,23%$), tổng photpho ($0,29%$), sắt tổng số ($8579,0\text{ mg/kg}$), pH ($6,64$) và độ dẫn điện EC ($0,31\text{ mS/cm}$).

Hệ thống đo đạc thực nghiệm sử dụng thiết bị đo đa kênh oxi hòa tan (Multi-Channel Dissolved Oxygen Meter – MCDM) tích hợp các cảm biến điện hóa màng thẩm thấu chọn lọc do Phòng Ứng dụng Tin học trong Hóa học (Viện Hóa học) tự chế tạo, kết nối liên tục với máy tính qua phần mềm điều khiển chuyên dụng MultiDO. Thiết bị cho phép thu nhận đồng thời tín hiệu nồng độ DO từ nhiều độ sâu với tần suất đo liên tục theo thời gian thực ($1\text{ lần/phút}$), loại bỏ hoàn toàn các sai số do thao tác khuấy trộn mẫu cơ học.

Độ tin cậy và giá trị nội tại của mô hình được bảo đảm qua quy trình hiệu chỉnh (calibration) và kiểm chứng (validation) độc lập trên các tập dữ liệu thực nghiệm khác nhau, sử dụng hệ số xác định ($R^2$), tổng bình phương phần dư (RSS) và tỉ số phân tích độ nhạy tham số (Sensitivity Ratio – SR).

Data và phân tích

Việc giải hệ phương trình vi phân đạo hàm riêng phi tuyến nhiều chiều được thực hiện trên phần mềm Comsol Multiphysics bằng phương pháp phần tử hữu hạn (FEM). Miền tính toán $\Omega$ được chia lưới phi cấu trúc (unstructured triangular/tetrahedral mesh) với mật độ lưới cực mịn tại các vùng ranh giới pha ($z = 0$ và $z = -h$) nhằm hạn chế sai số số trị.

Thuật toán giải số sử dụng phương pháp tích phân thời gian ẩn (implicit time-stepping) kết hợp bộ giải phi tuyến Newton – Raphson. Phân tích độ nhạy SR đối với các biến đầu vào cơ sở ($X_1 \rightarrow X_2$) và đầu ra ($Y_1 \rightarrow Y_2$) được định nghĩa: $$SR = \frac{(Y_2 - Y_1) / Y_1}{(X_2 - X_1) / X_1}$$ Kết quả chỉ ra rằng mô hình nhạy nhất với hệ số tiêu thụ DO của bùn đáy ($k_3$) và hệ số khuếch tán của oxi trong nước ($D_1$), khẳng định tính chính xác của cấu trúc toán học đã đề xuất.

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

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

  1. Sự hình thành vùng yếm khí cục bộ sát đáy: Dữ liệu thực nghiệm trích xuất trực tiếp từ văn bản cho thấy "Tốc độ khuếch tán của oxi trong nước chỉ bằng khoảng 1/10000 tốc độ khuếch tán của oxi trong khí quyển, chính vì vậy môi trường nước là môi trường có tốc độ khuếch tán oxi thấp". Do đó, tại lớp nước sát mặt bùn ($1,0\text{ cm}$), nồng độ DO sụt giảm nhanh chóng từ mức bão hòa ($8,2\text{ mg/L}$ ở $25^\circ\text{C}$) xuống dưới ngưỡng $2,0\text{ mg/L}$ (hypoxia) và tiến tới $0\text{ mg/L}$ (anoxia) chỉ sau $24-48\text{ giờ}$ trong các cột nước tĩnh không có nguồn cấp khí cưỡng bức.
  2. Ưu thế tuyệt đối của SOD so với BOD trong cột nước nông: Tốc độ tiêu thụ oxi của lớp bùn đáy dao động trong khoảng trích dẫn: "Tốc độ tiêu thụ oxi của bùn đáy nằm trong khoảng 50-500 mg O2/(m2 giờ)". Luận án chứng minh rằng đối với các cột nước có chiều sâu dưới $2,0\text{ m}$, hơn $70-85%$ lượng oxi hòa tan bị suy giảm bắt nguồn từ nhu cầu oxy hóa sinh – hóa của trầm tích đáy chứ không phải do quá trình phân hủy các chất hữu cơ lơ lửng trong cột nước (BOD).
  3. Hiện tượng suy giảm phân tầng phi tuyến theo chiều sâu: Dữ liệu mô phỏng và đo đạc tại các hồ Bảy Mẫu, Hồ Gươm, Hồ Tây cho thấy nồng độ DO biến thiên phi tuyến rõ rệt: nồng độ duy trì cao ở lớp mặt ($0-0,5\text{ m}$) do quá trình tái cấp khí bề mặt và quang hợp, nhưng giảm dốc đột ngột khi tiến gần lớp bùn đáy. Tương quan giữa mô phỏng FEM và đo đạc thực địa đạt $R^2 > 0,92$, chứng minh tính vững vững chắc của mô hình.
  4. Tác động xáo trộn và giải phóng chất khử: Khi có sự xáo trộn cơ học hoặc đối lưu nhiệt làm phát tán chất hữu cơ từ bùn vào nước ($k_{rese}$ tăng), tốc độ tiêu hao DO trong toàn bộ thể tích nước tăng vọt gấp $3-5\text{ lần}$ do sự oxy hóa tức thời các ion khử ($Fe^{2+}, Mn^{2+}, HS^-$) giải phóng từ bùn đáy: $$Fe^{2+} + \frac{1}{4}O_2 + \frac{5}{2}H_2O \rightarrow Fe(OH)_3\downarrow + 2H^+$$

Implications đa chiều

  • Đóng góp học thuật: Làm phong phú lý thuyết hóa lý môi trường nước, cung cấp công cụ toán học chính xác thay thế cho các công thức kinh nghiệm thô sơ khi tính toán thông lượng DO.
  • Đổi mới phương pháp luận: Thiết lập quy trình chuẩn mực kết hợp giữa phần cứng đo đạc đa kênh liên tục (MCDM) và phần mềm mô phỏng phần tử hữu hạn (Comsol Multiphysics).
  • Ứng dụng thực tiễn trong Nuôi trồng Thủy sản: Cung cấp cơ sở định lượng khoa học cho việc định kỳ nạo vét bùn đáy và tính toán công suất hệ thống quạt nước/sục khí tầng đáy trong các ao nuôi cá tra, tôm thâm canh, ngăn ngừa hiện tượng cá chết hàng loạt do ngạt khí vào ban đêm.
  • Khuyến nghị chính sách môi trường: Cung cấp cơ sở khoa học cho Bộ Tài nguyên và Môi trường cùng các cơ quan quản lý đô thị trong việc đánh giá tải lượng ô nhiễm hữu cơ, lập kế hoạch cải tạo và xử lý bùn đáy ô nhiễm tại các hồ nội đô bị phú dưỡng.

Limitations và Future Research

Luận án thừa nhận một cách khách quan các giới hạn nghiên cứu:

  • Giả định thủy vực tĩnh: Mô hình chưa tích hợp trường vận tốc dòng chảy 3 chiều phức tạp (3D hydrodynamic Navier – Stokes equations) và lực cắt do gió mạnh gây xáo trộn quy mô lớn.
  • Đơn giản hóa chu trình quang hợp – hô hấp: Trong các mô hình thành phần 2 và 3, tác động quang hợp của tảo được kiểm soát bằng cách che tối thực nghiệm để cô lập ảnh hưởng của bùn đáy, do đó chu trình ngày đêm sinh học theo hàm lượng diệp lục Chlorophyll-a ($C_{Chl-a}$) chưa được ghép nối toàn diện trong bài toán động lực học 3D.
  • Bỏ qua chuỗi phản ứng khử kị khí sâu: Chưa mô hình hóa chi tiết động thái vận chuyển các khí metan ($CH_4$) và hydro sunfua ($H_2S$) sinh ra từ tầng kị khí sâu trong trầm tích thoát lên cột nước.

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

  1. Ghép nối mô hình hóa lý Comsol với mô hình thủy động lực học 3D đa pha (CFD) có xét đến ảnh hưởng của sóng và dòng chảy đối lưu nhiệt.
  2. Tích hợp module sinh học phong phú mô tả quá trình nở hoa của vi tảo và hiện tượng phú dưỡng cực đoan (như trích dẫn văn bản: "Hiện tượng phú dưỡng có thể làm DO vượt quá giá trị DO bão hòa từ 1,5 đến 2 lần do sự quang hợp của rong tảo diễn ra mạnh vào ban ngày, trong khi đó quá trình hô hấp của chúng vào ban đêm có thể tiêu thụ hoàn toàn lượng oxi hòa tan").
  3. Phát triển thuật toán học máy (Machine Learning / PINNs – Physics-Informed Neural Networks) để đồng hóa dữ liệu cảm biến thời gian thực, phục vụ cảnh báo sớm ô nhiễm nguồn nước.

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

  • Tác động học thuật: Đặt nền móng cho các nghiên cứu mô phỏng chất lượng nước hiện đại tại Việt Nam, mở ra hướng nghiên cứu liên ngành giữa Hóa lý tính toán, Toán ứng dụng và Khoa học Môi trường.
  • Chuyển đổi ngành công nghiệp thủy sản: Ứng dụng mô hình giúp tối ưu hóa chu kỳ sục khí cơ học, tiết kiệm $20-30%$ chi phí năng lượng điện tiêu thụ trong các trang trại nuôi trồng thủy sản quy mô công nghiệp mà vẫn duy trì nồng độ DO trên ngưỡng an toàn sinh học ($>5,0\text{ mg/L}$).
  • Lợi ích xã hội và môi trường: Giúp các nhà quản lý dự báo chính xác thời điểm xảy ra sự cố suy giảm DO trong các hồ chứa cấp nước sinh hoạt, bảo vệ hệ sinh thái thủy sinh và sức khỏe cộng đồng.

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

  • Nghiên cứu sinh và Giảng viên ngành Hóa lý, Môi trường: Tiếp cận khung phương pháp giải số hệ PDE ứng dụng trong động học môi trường và kỹ thuật phân tích độ nhạy.
  • Kỹ sư R&D và Doanh nghiệp Thủy sản: Ứng dụng mô hình tính toán lượng bùn đáy tích lũy để thiết lập quy trình vận hành và xử lý đáy ao tối ưu.
  • Cơ quan quản lý nhà nước và Đô thị: Ứng dụng phần mềm dự báo để thẩm định các dự án nạo vét hồ đô thị và đánh giá tác động môi trường các nguồn thải.

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ì? Trả lời: Đó là việc xây dựng thành công hệ phương trình động học mô phỏng phân bố nồng độ DO theo không gian 3 chiều và thời gian, trong đó lượng hóa chính xác nhu cầu oxi bùn đáy (SOD) như một điều kiện biên thông lượng phi tuyến phụ thuộc cấu trúc vi mô của bùn và nhiệt độ ($\theta_s^{(T-20)}$), vượt qua hạn chế của mô hình hộp Streeter – Phelps cổ điển.
  2. Tính đổi mới trong phương pháp nghiên cứu khi so sánh với các công trình trước đây? Trả lời: Khác với các nghiên cứu của O'Connor (1958) hay Thomann (1987) chỉ giải tích phân 1 chiều gần đúng, luận án đã sử dụng phương pháp phần tử hữu hạn (FEM) trên Comsol Multiphysics kết hợp với hệ thiết bị đo đa kênh MCDM tự chế tạo thu nhận dữ liệu đa tầng liên tục, tạo ra sự liên kết khép kín giữa lý thuyết và thực nghiệm.
  3. Phát hiện nào gây bất ngờ và có ý nghĩa thực tiễn cao nhất? Trả lời: Mức độ suy giảm DO ở lớp biên sát mặt bùn ($1\text{ cm}$) diễn ra cực kỳ khốc liệt và nhanh chóng, tạo thành "vùng chết" (dead zone) yếm khí ngay cả khi lớp nước mặt vẫn ở trạng thái bão hòa DO, nguyên nhân do trở lực khuếch tán phân tử của nước quá lớn ($D_1 \approx 5 \times 10^{-5}\text{ cm}^2/\text{s}$).
  4. Luận án có cung cấp quy trình lặp lại thực nghiệm (replication protocol) không? Trả lời: Có. Luận án mô tả chi tiết sơ đồ thiết kế cột nước vật lý, cấu tạo cảm biến điện hóa, các bước thiết lập module PDE trong Comsol Multiphysics từ khâu chia lưới, nhập điều kiện ban đầu, khai báo điều kiện biên Neumann/Dirichlet đến thuật toán giải số.
  5. Chương trình nghiên cứu 10 năm được định hình như thế nào? Trả lời: Mở rộng mô hình sang thủy vực chuyển động 3 chiều, tích hợp các chu trình sinh địa hóa của Nitơ, Photpho, Lưu huỳnh và phát triển hệ thống cảnh báo tự động chất lượng nước dựa trên trí tuệ nhân tạo.

Kết luận

  1. Đã thiết lập thành công hệ mô hình toán học giải tích – số hoàn chỉnh mô tả sự biến thiên nồng độ oxi hòa tan (DO) trong môi trường nước ngọt dưới tác động chủ đạo của lớp bùn đáy và các thông số hóa lý liên quan.
  2. Đã giải quyết triệt để bài toán biên lan truyền – khuếch tán đa chiều bằng phương pháp phần tử hữu hạn (FEM) trên phần mềm Comsol Multiphysics với độ chính xác cao.
  3. Đã chế tạo và ứng dụng thành công thiết bị đo đa kênh DO (MCDM) và phần mềm MultiDO để thu thập chuỗi số liệu thực nghiệm liên tục, kiểm chứng mô hình đạt hệ số tương quan $R^2 > 0,90$.
  4. Đã làm sáng tỏ vai trò chi phối của nhu cầu oxi bùn đáy (SOD) và trở lực khuếch tán phân tử trong việc hình thành phân tầng DO theo chiều sâu.
  5. Đã kiểm chứng tính ứng dụng thực tiễn của mô hình trên các hồ tự nhiên tại khu vực Hà Nội và Hòa Bình, mở ra công cụ hữu hiệu trong dự báo và quản lý chất lượng môi trường nước.