Tổng quan luận án

Luận án tiến sĩ kỹ thuật của tác giả Thanh Thuy Nguyen, thực hiện dưới sự hướng dẫn của Giáo sư Patrick Willems tại Khoa Khoa học Kỹ thuật (Faculty of Engineering Science), Đại học KU Leuven (Vương quốc Bỉ, bảo vệ năm 2017), tập trung vào đề tài đánh giá chất lượng nước sông bằng phương pháp mô hình hóa dưới điều kiện khí hậu hiện tại và tương lai (Model Based River Water Quality Assessment Under Current and Future Climate Conditions).

Về tính cấp thiết, ô nhiễm nguồn nước sông là vấn đề môi trường nghiêm trọng trên toàn cầu, đe dọa trực tiếp và gián tiếp đến sức khỏe con người. Báo cáo của Liên Hợp Quốc (UN 2000) và WHO/UNICEF (2005) chỉ ra rằng các bệnh liên quan đến nước chiếm tới 50% số bệnh nhân nội trú và gây ra gần 1/5 số ca tử vong ở trẻ dưới 5 tuổi trên thế giới. Biến đổi khí hậu với sự gia tăng nhiệt độ không khí, biến động dòng chảy cực đoan và thay đổi hệ sinh thái dự kiến sẽ làm suy giảm chất lượng nước sông. Để phục vụ công tác quản lý lưu vực tích hợp, việc cung cấp thông tin định lượng chính xác về diễn biến chất lượng nước là điều kiện tiên quyết.

Khoảng trống nghiên cứu được tác giả xác định xuất phát từ những giới hạn thực tế:

  • Dữ liệu quan trắc chất lượng nước thường bị giới hạn về không gian (chỉ đo tại một số vị trí cố định) và thời gian (tần suất đo thưa), đồng thời chịu sai số đo đạc và lấy mẫu.
  • Các bộ mô hình thủy lực và chất lượng nước chi tiết 1D dựa trên phương trình thủy động lực học đầy đủ (như MIKE 11, InfoWorks RS, InfoWorks ICM) đòi hỏi bước thời gian mô phỏng nhỏ để đảm bảo hội tụ và tiêu tốn thời gian tính toán rất lớn. Điều này khiến chúng không khả thi đối với các ứng dụng yêu cầu số lượng kịch bản chạy lớn như phân tích độ bất định (phương pháp Monte Carlo) hoặc đánh giá tác động biến đổi khí hậu dài hạn (với nhiều mô hình khí hậu toàn cầu GCM).
  • Giả định hệ số nhám lòng sông Manning $n$ không đổi theo thời gian trong các mô hình thủy động lực học truyền thống bỏ qua ảnh hưởng của sự phát triển thực vật thủy sinh theo mùa, gây sai lệch lớn đến tính toán mực nước mùa kiệt và vận tốc dòng chảy.
  • Thiếu hụt các nghiên cứu định lượng độc lập về độ bất định cấu trúc mô hình và độ bất định dữ liệu đầu vào (tải lượng ô nhiễm điểm và phân tán) trong mô phỏng chất lượng nước sông.

Mục tiêu nghiên cứu của luận án được đánh số cụ thể gồm ba mục tiêu chính:

  1. Đánh giá và so sánh định lượng các mô hình chất lượng nước sông chi tiết hiện có (MIKE 11, InfoWorks RS, InfoWorks ICM); làm rõ đóng góp của sự khác biệt trong phương trình chuyển hóa lý - hóa - sinh, thông số và vai trò của kết quả thủy động lực học (bao gồm biến thiên nhám theo mùa do thực vật thủy sinh).
  2. Phát triển mô hình chất lượng nước sông dạng khái niệm CORIWAQ dựa trên nền tảng phần mềm InfoWorks RS (CORIWAQ-RS), mở rộng cấu trúc mô hình linh hoạt từ phiên bản CORIWAQ-MIKE11 trước đó, giúp giảm thời gian tính toán trong khi vẫn bảo đảm độ chính xác so với mô hình tham chiếu chi tiết.
  3. Phân tích độ bất định và đánh giá tác động của biến đổi khí hậu đến chất lượng nước sông: xác định các nguồn dữ liệu đầu vào nhạy cảm nhất, định lượng đóng góp của độ bất định đầu vào vào tổng sai số mô hình, và đánh giá tác động của 30 kịch bản biến đổi khí hậu đến dòng chảy cực đoan cùng nồng độ các chất ô nhiễm.

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

  • Đối tượng nghiên cứu: Quá trình truyền triệt và biến đổi lý - hóa - sinh của các thông số chất lượng nước sông (nhiệt độ nước, oxy hòa tan DO, nhu cầu oxy sinh hóa $\text{BOD}_5$, amoni $\text{NH}_4\text{-N}$, nitrat $\text{NO}_3\text{-N}$, nitơ hữu cơ $\text{ON-N}$, nitơ tổng $\text{TN}$, orthophosphate $\text{PO}_4\text{-P}$, phospho dạng hạt $\text{PP-P}$); các phương pháp mô hình hóa thủy văn, thủy động lực học, độ nhám lòng sông và độ bất định mô hình.
  • Phạm vi không gian: Lưu vực sông Molse Neet thuộc lưu vực sông Grote Nete (vùng Flanders, Vương quốc Bỉ) với diện tích khoảng 24,54 $\text{km}^2$, trải dài qua ba đô thị Lommel (tỉnh Limburg), Balen và Mol (tỉnh Antwerp).
  • Phạm vi thời gian: Giai đoạn quan trắc lịch sử phục vụ hiệu chỉnh và kiểm định mô hình là 2000–2010 (dữ liệu chi tiết năm tham chiếu 2006); giai đoạn tương lai đánh giá biến đổi khí hậu là mốc thời gian 2050 (chu kỳ 2050–2060).

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

Luận án đã tổng quan các hướng nghiên cứu lý thuyết và ứng dụng trong mô hình hóa chất lượng nước và đánh giá tác động môi trường:

  • Mô hình hóa chất lượng nước truyền thống và phân loại mô hình: Nền tảng phương trình Streeter và Phelps (1925) về tương quan DO-BOD; hệ thống phân loại mô hình từ đơn giản (như TOMCAT), trung bình (như QUAL2E) đến phức tạp (như Delft 3D) theo Cox (2003); các mô hình mô phỏng dòng chảy ngầm và bổ cập dòng chảy (Peters et al. 2003; Willems 2014).
  • Mô hình hóa mô phỏng rút gọn (Emulation modeling): Cách tiếp cận xấp xỉ bậc thấp từ các mô hình chi tiết dựa trên vật lý nhằm giảm tải thời gian tính toán (Castelletti et al. 2012); các phương pháp lựa chọn biến số dựa trên dữ liệu (Bowden et al. 2008; Hejazi và Cai 2009) hoặc dựa trên cấu trúc mô hình (Cox et al. 2006; Crout et al. 2009).
  • Độ bất định trong mô hình hóa chất lượng nước: Phân loại độ bất định thành cấu trúc mô hình, thông số, dữ liệu đầu vào và sai số đo đạc (Beck 1987; Tung và Yen 2006; Harmel và Smith 2007; Mannina và Viviani 2010). Các nghiên cứu định lượng độ bất định cấu trúc bằng cách so sánh nhiều cấu trúc mô hình khả dĩ (Van der Perk 1997) hoặc phân tách phương sai phần dư (Radwan et al. 2004; Freni và Mannina 2010; Willems 2012).
  • Xử lý dữ liệu đầu vào khuyết thiếu: Phương pháp truyền thống như nội suy trung bình, hồi quy (Little et al. 2006; Kalteh và Hjorth 2009) và phương pháp hiện đại như ước lượng hợp lý cực đại, ước lượng Bayes và đa nội suy (Enders 2010; NIST/SEMATECH 2012). Phương pháp vùng hóa lưu vực không quan trắc (Young 2006; Parajka et al. 2007). Phân tích độ nhạy cục bộ và toàn cục (Saltelli et al. 2008; Norton 2015).
  • Biến thiên nhám lòng sông do thực vật thủy sinh: Ảnh hưởng của thảm thực vật thủy sinh (macrophytes) đến sức cản dòng chảy, vận tốc và khả năng tự làm sạch (Watson 1987; Bakry et al. 1992; Gurnell và Midgley 1994; Schulz et al. 2003; De Doncker et al. 2011; Keupers et al. 2015).
  • Tác động của biến đổi khí hậu đến chất lượng nước: Ứng dụng mô hình SWAT (Neitsch et al. 2011; Bouraoui et al. 2002; Van Liew et al. 2012), INCA-N (Whitehead et al. 1998; Wilby et al. 2006), phân tích chu kỳ nhiệt độ tích lũy (Growing Degree Days - GDD) và thời gian sinh trưởng của cây trồng (Length of Growing Period - LGP) (Miller et al. 2001).

Khoảng trống luận án giải quyết: Luận án lấp đầy khoảng trống về công cụ mô hình hóa trung gian có khả năng kế thừa độ chính xác của các mô hình thủy động lực học chi tiết (InfoWorks RS và MIKE 11) nhưng đạt tốc độ tính toán nhanh để phục vụ phân tích độ bất định Monte Carlo và chạy 30 kịch bản biến đổi khí hậu; đồng thời cung cấp phương pháp luận định lượng ảnh hưởng của nhám lòng sông biến thiên theo mùa và phân tích đóng góp của từng nguồn thải đến sai số dự báo chất lượng nước.


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

Luận án tích hợp các khung lý thuyết và phương pháp nghiên cứu định lượng sau:

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

  • Lý thuyết thủy động lực học và vận chuyển chất: Phương trình Saint-Venant 1D mô tả dòng chảy không ổn định; phương trình bình lưu - khuếch tán (Advection-Diffusion - AD) mô tả sự lan truyền chất ô nhiễm.
  • Động học biến đổi lý - hóa - sinh: Hệ phương trình chuyển hóa chất ô nhiễm mô tả các quá trình: tái sục khí bề mặt (reaeration), phân hủy chất hữu cơ ($\text{BOD}$), nitrat hóa amoni thành nitrit/nitrat ($\text{NH}_4 \rightarrow \text{NO}_2 \rightarrow \text{NO}_3$), khử nitrat (denitrification), quang hợp và hô hấp của thực vật/tảo, lắng đọng và hấp thụ dinh dưỡng.
  • Khái niệm hóa dòng chảy và chất lượng nước (CORIWAQ): Biểu diễn mạng lưới sông bằng chuỗi các bể chứa dung tích biến đổi (cascades of reservoirs), gộp (lump) các quá trình bình lưu - khuếch tán và chuyển hóa sinh hóa theo từng đoạn sông tương ứng với từng bể chứa.

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

  1. Mô hình hóa thủy văn và tải lượng dinh dưỡng:
    • Mô hình thủy văn PDM (Probability Distribution Model): sử dụng phân bố Pareto mô tả biến thiên không gian của dung tích trữ ẩm đất, định tuyến dòng chảy nhanh qua 2 bể chứa tuyến tính và dòng chảy cơ bản qua bể chứa phi tuyến bậc 3.
    • Mô hình VHM (Versalgemeend conceptueel Hydrologisch Model): mô hình thủy văn khái niệm hóa tổng quát hỗ trợ phân tách các thành phần dòng chảy.
    • Mô hình SENTWA (System for the Evaluation of Nutrient Transport to Water): ước tính tải lượng nitơ từ các hoạt động nông nghiệp vào mạng lưới sông.
  2. Mô hình hóa chất lượng nước chi tiết và khái niệm hóa:
    • Thiết lập mô hình chi tiết trên MIKE 11 (DHI), InfoWorks RS và InfoWorks ICM (Innovyze).
    • Phát triển thuật toán CORIWAQ-RS trong môi trường MATLAB, hiệu chỉnh thông số bể chứa dựa trên kết quả trích xuất từ InfoWorks RS.
  3. Phương pháp phân tích độ nhạy và độ bất định:
    • Phân tích độ nhạy cục bộ (Local Sensitivity Analysis) và phân tích độ nhạy đồ thị để đánh giá ảnh hưởng của các thông số và cấu trúc mô hình.
    • Phương pháp hồi quy ngẫu nhiên (Stochastic regression): thêm chuỗi sai số ngẫu nhiên theo phân bố xác định $\mathcal{N}(\mu, \sigma)$ vào phương trình hồi quy để mô hình hóa độ bất định dữ liệu đầu vào.
    • Kỹ thuật mô phỏng Monte Carlo kết hợp phương pháp phân tách phương sai (Variance decomposition) để tính toán chỉ số nhạy $S_i$ và lượng hóa đóng góp của từng nguồn thải vào phương sai sai số đầu ra.
  4. Phương pháp đánh giá tác động biến đổi khí hậu:
    • Thu nhỏ thống kê (statistical downscaling) lượng mưa, nhiệt độ không khí và bốc thoát hơi nước tiềm năng ($\text{ET}_p$) từ 30 kịch bản chạy mô hình khí hậu toàn cầu (GCM) thuộc 4 kịch bản phát thải khí nhà kính.
    • Công cụ WETSPRO (Water Engineering Time Series PROcessing tool) và phương pháp trích xuất đỉnh vượt ngưỡng (Peak-Over-Threshold - POT) để phân tích lưu lượng cực đoan và nồng độ chất ô nhiễm cực đoan theo chu kỳ lặp lại (return periods).

Nguồn tư liệu và dữ liệu thực nghiệm

  • Dữ liệu địa hình và mạng lưới: Bản đồ độ cao số, bình đồ các mặt cắt ngang sông Molse Neet cách nhau trung bình 50 m, thông số các công trình thủy công (cống, đập tràn, cầu).
  • Dữ liệu trầm tích đáy: Mẫu trầm tích lấy vào tháng 8/2015 xác định cấp phối hạt từ 63 đến 2000 $\mu\text{m}$, đường kính trung vị $D_{50} = 200\ \mu\text{m}$, bề dày lớp bùn đáy từ 0 cm đến 85 cm.
  • Dữ liệu khí tượng: Lượng mưa, nhiệt độ không khí, $\text{ET}_p$ giai đoạn 2000–2010 cung cấp bởi Viện Khí tượng Hoàng gia Bỉ (KMI/RMI).
  • Dữ liệu quan trắc chất lượng nước: Dữ liệu quan trắc từ Cơ quan Môi trường Flanders (VMM) tại trạm 333100 (thượng lưu) và trạm 333000 (hạ lưu, cách biên thượng lưu 6 km) với tần suất đo khoảng 1 lần/tháng.
  • Dữ liệu nguồn thải: Dữ liệu tải lượng công nghiệp năm cơ sở 2006 từ mô hình PEGASE và VMM; dữ liệu dân số phát thải sinh hoạt không qua xử lý; số liệu canh tác và chăn nuôi nông nghiệp.

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

Chapter 1: Introduction

Chương 1 thiết lập bối cảnh nghiên cứu về ô nhiễm nguồn nước sông, tác động tiềm tàng của biến đổi khí hậu và vai trò của mô hình hóa toán học trong quản lý lưu vực tích hợp. Tác giả tổng thuật các phương pháp mô hình hóa hiện có, nêu rõ các hạn chế về thời gian tính toán của mô hình chi tiết, ảnh hưởng của độ nhám do thực vật thủy sinh, sự thiếu hụt dữ liệu đầu vào và các loại độ bất định. Chương này xác định rõ 3 mục tiêu cụ thể và cấu trúc tổng thể 9 chương của luận án.

Chapter 2: Molse Neet catchment

Chương 2 mô tả chi tiết các đặc trưng địa lý, khí hậu, thủy văn và nguồn thải ô nhiễm của lưu vực sông Molse Neet (24,54 $\text{km}^2$, độ dốc trung bình 0,3‰, cao độ từ 82 m TAW giảm dần về 26 m TAW). Sông Molse Neet hợp lưu với sông Schoppelijke Neet tại km 11,5. Lưu vực có lớp đất cát dày (~0,3 m), mực nước ngầm nông (thay đổi từ 0,4–1,0 m dưới mặt đất vào mùa hè và chỉ cách mặt đất 0,15 m vào các mùa khác). Thảm thực vật thủy sinh phát triển mạnh từ mùa xuân đến mùa thu, phổ biến gồm các loài: Lemna minor, Bur-Reed, rong mái chèo (pondweed), Sterrenkroos.

Chương tổng hợp số liệu về ba nhóm nguồn thải chính tác động lên chất lượng nước sông Molse Neet:

Nhóm nguồn thải Tên nguồn / Mã nguồn Đặc điểm xả thải / Quy mô
Công nghiệp Ajinomoto Omnichem Nhà máy dược phẩm, xả thải trực tiếp, quan trắc bởi VMM
Philips EMGO Nhà máy sản xuất bóng đèn/thủy tinh, xả thải trực tiếp
MORA Food Nhà máy chế biến thực phẩm, tính toán theo thuế/đô thị
Sinh hoạt chưa xử lý UE-114363 Nguồn điểm: 112 dân số tương đương (IE)
UE-114372 Nguồn điểm: 125 IE
UE-114373 Nguồn điểm: 60 IE
UD-501 (VHA501) Nguồn phân tán: 3.510 IE
Nông nghiệp Trồng trọt & chăn nuôi Đồng cỏ cắt cỏ (37,9%), ngô (22,9%), bãi chăn thả (15,9%), khoai tây (1,0%)

Chapter 3: Existing models and methods

Chương 3 trình bày cơ sở toán học và cấu trúc của các mô hình được sử dụng trong luận án:

  • Mô hình thủy văn PDM và VHM để chuyển đổi lượng mưa và bốc thoát hơi nước thành dòng chảy mặt và dòng chảy ngầm cung cấp cho mô hình chất lượng nước.
  • Các phương trình bình lưu - khuếch tán và phương trình động học chuyển hóa sinh hóa trong ba gói phần mềm MIKE 11, InfoWorks RS và InfoWorks ICM.
  • Phương pháp ước tính dữ liệu đầu vào bị thiếu cho các nguồn thải công nghiệp và sinh hoạt; mô hình SENTWA mô phỏng tải lượng đạm từ nông nghiệp.
  • Giới thiệu bộ công cụ biến đổi khí hậu dựa trên phương pháp nhiễu động (perturbation tool) và phương pháp trích xuất đỉnh vượt ngưỡng POT (Peak-Over-Threshold) bằng công cụ WETSPRO.

Chapter 4: Development of a COnceptual River WAter Quality model (CORIWAQ) with flexible model structure

Chương 4 trình bày thuật toán xây dựng mô hình khái niệm CORIWAQ-RS mở rộng từ nền tảng InfoWorks RS, kế thừa công trình ban đầu của Keupers (2014) trên MIKE 11. Sông được chia thành chuỗi các bể chứa nối tiếp; các quá trình thủy lực và chuyển hóa sinh hóa được gộp theo từng đoạn bể chứa.

  • Kết quả hiệu chỉnh và kiểm định trên lưu vực sông Molse Neet cho thấy CORIWAQ-RS và CORIWAQ-MIKE11 mô phỏng các thông số nhiệt độ, DO, $\text{NH}_4\text{-N}$, $\text{NO}_3\text{-N}$, $\text{BOD}_5$, $\text{PO}_4\text{-P}$ và $\text{PP-P}$ với độ chính xác tương đương các mô hình chi tiết gốc (chỉ số hiệu quả Nash-Sutcliffe NSE duy trì ở mức cao qua các năm khô, trung bình và ẩm).
  • Thời gian tính toán của mô hình khái niệm CORIWAQ nhanh gấp 10 lần so với các mô hình thủy động lực học đầy đủ ban đầu.

Chapter 5: Seasonal variation of river bed roughness impacts on water level

Chương 5 nghiên cứu ảnh hưởng của sự biến thiên theo mùa của độ nhám lòng sông (do thực vật thủy sinh phát triển) đến mực nước bằng hai phần mềm MIKE 11 và InfoWorks RS:

  • Hệ số Manning $n$ được hiệu chỉnh theo từng khoảng thời gian 5 ngày dựa trên mực nước quan trắc tại trạm Meerhout. Giá trị Manning $n$ tăng cao rõ rệt vào mùa hè khi thảm thực vật thủy sinh phát triển đạt đỉnh.
  • Nếu sử dụng hệ số Manning $n$ cố định (bỏ qua thực vật thủy sinh), mô hình đánh giá thấp (underelevation) đáng kể mực nước trong mùa hè.
  • Trong các giai đoạn dòng chảy đỉnh (lũ lớn), ảnh hưởng của nhám thực vật trở nên không đáng kể do thảm thực vật bị uốn rạp hoặc ngập sâu. Tuy nhiên, do nồng độ chất ô nhiễm cao nhất thường xuất hiện vào thời kỳ mực nước thấp (mùa kiệt), việc mô hình hóa biến thiên nhám theo mùa đóng vai trò quyết định đến độ chính xác của mô hình chất lượng nước.
  • Kết quả so sánh giữa InfoWorks RS và MIKE 11 cho thấy cả hai phần mềm đều mô phỏng tốt xu thế biến thiên này, trong đó InfoWorks RS cho kết quả bám sát hơn một chút nhưng mức độ chênh lệch giữa hai mô hình là tương đối nhỏ.

Chapter 6: Influence of model structure uncertainty on river water quality assessment

Chương 6 tiến hành phân tích độ bất định cấu trúc mô hình chất lượng nước thông qua so sánh ba gói phần mềm MIKE 11, InfoWorks RS và InfoWorks ICM:

  • Tác giả thực hiện sàng lọc các quá trình hóa sinh và phân tích độ nhạy cục bộ (Local sensitivity analysis) để tìm ra các quá trình then chốt gây sai khác kết quả.
  • Kết quả chỉ ra rằng mô đun thủy động lực học, tốc độ truyền nhiệt bề mặt và công thức tái sục khí tự nhiên (reaeration) là những nguyên nhân chính dẫn đến sự khác biệt lớn giữa các mô hình đối với nhiệt độ nước và nồng độ DO.
  • Việc bỏ qua quá trình vận chuyển bùn cát, thực vật phù du (phytoplankton) và vi khuẩn dẫn đến ước tính nồng độ chất hữu cơ cao hơn và nồng độ oxy hòa tan thấp hơn.
  • Cả ba cấu trúc mô hình đạt được sự đồng thuận cao khi xác định các nguồn ô nhiễm chính chi phối nồng độ chất hữu cơ và nitrat, nhưng lại có sự bất đồng đáng kể trong việc xác định các yếu tố chính chi phối nồng độ DO.

Chapter 7: CORIWAQ input uncertainty analysis

Chương 7 ứng dụng mô hình khái niệm CORIWAQ-RS để phân tích định lượng độ bất định dữ liệu đầu vào:

  • Xử lý số liệu khuyết tật và mô hình hóa ngẫu nhiên chuỗi sai số đầu vào của các nguồn thải công nghiệp (nhà máy Ajinomoto, Philips), tải lượng sinh hoạt và nông nghiệp bằng kỹ thuật hồi quy ngẫu nhiên (stochastic regression).
  • Áp dụng mô phỏng Monte Carlo để lan truyền độ bất định qua mô hình CORIWAQ và tính toán các chỉ số nhạy $S_i$ cho các thông số chất lượng nước tại trạm quan trắc 333000 theo các mùa trong năm.
  • Áp dụng phương pháp phân tách phương sai để định lượng mức đóng góp của từng nguồn xả vào tổng phương sai phần dư của mô hình, giúp xác định các nguồn thải trọng yếu cần ưu tiên thu thập thêm dữ liệu quan trắc để nâng cao độ tin cậy của mô hình với chi phí tối ưu.

Chapter 8: Climate change impacts on river water quality

Chương 8 phân tích tác động của biến đổi khí hậu đến chất lượng nước sông Molse Neet bằng cách kết hợp mô hình CORIWAQ-RS với kết quả thu nhỏ thống kê từ 30 kịch bản GCM (giai đoạn 2000–2010 so với 2050):

  • Tín hiệu biến đổi khí hậu được lan truyền qua các mô hình thủy văn, quan hệ hồi quy và phương pháp ngẫu nhiên để xác định sự thay đổi dòng chảy tràn, tải lượng nitơ từ lưu vực, nhiệt độ nước và nồng độ DO bão hòa.
  • Nhiệt độ tăng làm rút ngắn chu kỳ sinh trưởng của cây trồng (LGP) dựa trên tổng nhiệt hữu hiệu (GDD) và nhiệt độ cơ sở ($Temp_{base}$), dẫn đến thời gian đất không được che phủ kéo dài hơn và làm tăng tổn thất nitơ từ phân bón nông nghiệp vào mùa đông.
  • Kết quả chỉ ra biến đổi khí hậu có tác động tiêu cực mạnh mẽ đến nồng độ oxy hòa tan (DO giảm sâu trong mùa kiệt). Mức độ tác động của biến đổi khí hậu có biên độ bất định rất lớn, đặc biệt là đối với các biến cố có chu kỳ lặp lại cao (các hiện tượng cực đoan).

Chapter 9: Conclusions and perspectives

Chương 9 tổng kết các phát hiện chính, khẳng định tính hiệu quả của mô hình khái niệm CORIWAQ-RS, lượng hóa các nguồn bất định và tác động của biến đổi khí hậu; đồng thời đề xuất các định hướng phát triển trong tương lai cho nền tảng CORIWAQ (mở rộng khả năng tích hợp mô hình chất lượng nước đầy đủ hơn, nâng cao độ phân giải dữ liệu và mở rộng cho các lưu vực khác).


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

Đóng góp về mặt lý luận và phương pháp

  • Mở rộng thành công nền tảng mô hình khái niệm CORIWAQ dựa trên cấu trúc quá trình và kết quả của phần mềm InfoWorks RS (CORIWAQ-RS), thiết lập khung mô hình chất lượng nước linh hoạt cho phép chuyển đổi giữa các cấu trúc quá trình khác nhau.
  • Chứng minh rằng mô hình khái niệm hóa dạng chuỗi bể chứa kết hợp gộp các quá trình bình lưu - khuếch tán và chuyển hóa hóa sinh đạt độ chính xác tương đương các mô hình thủy lực chi tiết 1D (MIKE 11, InfoWorks RS), nhưng giảm thời gian tính toán đi 10 lần.
  • Xây dựng quy trình tích hợp kết hợp hồi quy ngẫu nhiên, phân tích độ nhạy toàn cục và phân tách phương sai để định lượng mức độ đóng góp riêng biệt của từng nguồn ô nhiễm đầu vào vào tổng sai số dự báo chất lượng nước sông.
  • Phát triển khung phương pháp hoàn chỉnh để lan truyền các tín hiệu biến đổi khí hậu từ mô hình GCM sang mô hình thủy văn lưu vực, tải lượng dinh dưỡng nông nghiệp (LGP, GDD), nhiệt độ nước và chất lượng nước sông.

Đóng góp về mặt thực tiễn và ứng dụng

  • Chứng minh định lượng việc bỏ qua biến thiên nhám theo mùa do thực vật thủy sinh sẽ dẫn đến đánh giá thấp mực nước mùa hè và sai lệch vận tốc dòng chảy, ảnh hưởng trực tiếp đến tính toán chất lượng nước trong mùa kiệt.
  • Cung cấp bức tranh định lượng về tác động của biến đổi khí hậu đến sông Molse Neet: nồng độ oxy hòa tan (DO) suy giảm nghiêm trọng vào mùa hè ở mốc năm 2050, trong khi tải lượng nitơ có xu hướng gia tăng vào mùa đông.
  • Đưa ra giải pháp hỗ trợ cơ quan quản lý nước (như Cơ quan Môi trường Flanders - VMM) xác định chính xác các nguồn thải trọng điểm cần đầu tư kinh phí quan trắc bổ sung, tối ưu hóa hiệu quả kinh tế trong quản lý chất lượng nước lưu vực.

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

Hạn chế của nghiên cứu (theo văn bản)

  • Mô hình khái niệm CORIWAQ phụ thuộc vào dữ liệu đầu vào thủy động lực học trích xuất từ các mô hình chi tiết ban đầu; các thông số hình học và thủy lực được gộp theo từng đoạn sông.
  • Chuỗi dữ liệu quan trắc chất lượng nước lịch sử tại lưu vực nghiên cứu còn thưa (tần suất khoảng 1 tháng/lần), một số chỉ tiêu như $\text{BOD}_5$ và $\text{NH}_4\text{-N}$ rất khan hiếm, trạm thượng lưu 333100 bị gián đoạn quan trắc trong giai đoạn 2007–2010.
  • Các quá trình liên quan đến vận chuyển bùn cát chi tiết, động học vi khuẩn và thực vật phù du chưa được mô hình hóa toàn diện trong một số cấu trúc mô hình so sánh.

Định hướng nghiên cứu tiếp theo

  • Nâng cấp và hoàn thiện mô hình CORIWAQ: Tích hợp thêm các phương trình động học mô tả chi tiết hơn quá trình bùn đáy, thực vật phù du và trao đổi chất dinh dưỡng nâng cao; tối ưu hóa thuật toán hiệu chuẩn thông số tự động.
  • Cải thiện chất lượng dữ liệu và kết quả dự báo: Áp dụng khung mô hình cho các lưu vực có chuỗi số liệu quan trắc liên tục và độ phân giải cao hơn; nghiên cứu sâu hơn về tính bất định của các kịch bản khí hậu cực đoan với chu kỳ lặp lại lớn.

Giá trị tham khảo

Luận án là tài liệu tham khảo hữu ích cho các nhóm đối tượng học thuật và ứng dụng chuyên ngành:

  • Nghiên cứu sinh và nhà nghiên cứu: Tham khảo phương pháp xây dựng mô hình khái niệm hóa chất lượng nước (CORIWAQ), kỹ thuật phân tích độ nhạy, mô phỏng ngẫu nhiên Monte Carlo và quy trình phân tách phương sai để xác định độ bất định mô hình.
  • Giảng viên và chuyên gia mô hình hóa: Sử dụng làm tài liệu giảng dạy chuyên đề về mô hình hóa thủy lực sông 1D (so sánh MIKE 11, InfoWorks RS, InfoWorks ICM), tương tác giữa thủy lực - thực vật thủy sinh và phương pháp lan truyền kịch bản biến đổi khí hậu vào chất lượng nước lưu vực.
  • Nhà quản lý tài nguyên nước và môi trường (Cơ quan môi trường, Ban quản lý lưu vực): Ứng dụng khung đánh giá để xác định nguồn ô nhiễm trọng điểm, xây dựng kế hoạch quan trắc tối ưu chi phí và hoạch định các chiến lược thích ứng với biến đổi khí hậu trong quản lý chất lượng nước sông.

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

1. Mô hình khái niệm CORIWAQ-RS có ưu điểm gì vượt trội so với các mô hình thủy lực chi tiết như InfoWorks RS hay MIKE 11?
CORIWAQ-RS đạt độ chính xác mô phỏng tương đương với mô hình chi tiết gốc nhưng có tốc độ tính toán nhanh hơn gấp 10 lần. Điều này giúp mô hình thực hiện được các tác vụ đòi hỏi hàng nghìn lần chạy mô phỏng như phân tích độ bất định ngẫu nhiên Monte Carlo và đánh giá 30 kịch bản biến đổi khí hậu dài hạn mà mô hình chi tiết không thể đáp ứng được về mặt thời gian.

2. Sự phát triển của thực vật thủy sinh theo mùa ảnh hưởng như thế nào đến chế độ thủy lực sông Molse Neet?
Thực vật thủy sinh phát triển mạnh nhất vào giữa mùa hè, làm gia tăng hệ số nhám lòng sông Manning $n$. Nếu giả định nhám không đổi, mô hình sẽ tính toán thấp hơn thực tế đối với mực nước mùa hè. Trong các đợt lũ đỉnh, ảnh hưởng của thảm thực vật trở nên không đáng kể. Tuy nhiên, do nồng độ các chất ô nhiễm đạt mức nghiêm trọng nhất vào mùa kiệt (mực nước thấp), việc tính toán đúng biến thiên nhám theo mùa là yếu tố bắt buộc để mô phỏng chính xác vận tốc và chất lượng nước.

3. Ba gói phần mềm MIKE 11, InfoWorks RS và InfoWorks ICM có sự thống nhất và khác biệt gì trong mô phỏng chất lượng nước?
Cả ba mô hình đều đạt sự đồng thuận cao khi xác định các nguồn ô nhiễm chính giải thích nồng độ chất hữu cơ ($\text{BOD}$) và nitrat ($\text{NO}_3\text{-N}$). Tuy nhiên, các mô hình có sự bất đồng đáng kể về các yếu tố chi phối nồng độ oxy hòa tan (DO) do sự khác biệt trong công thức tính toán tái sục khí bề mặt, tốc độ truyền nhiệt và cách thức mô tả các quá trình bùn đáy, thực vật phù du.

4. Biến đổi khí hậu tác động như thế nào đến tải lượng nitơ nông nghiệp và oxy hòa tan (DO) trong lưu vực?
Nhiệt độ không khí tăng làm tăng tích ôn hữu hiệu (GDD) và rút ngắn thời gian sinh trưởng của cây trồng (LGP), khiến thời gian đất trống không canh tác kéo dài hơn, dẫn đến rửa trôi và thất thoát nitơ từ đất nông nghiệp vào sông nhiều hơn vào mùa đông. Vào mùa hè, nhiệt độ nước tăng làm giảm nồng độ oxy bão hòa (DOS) và thúc đẩy tốc độ phản ứng phân hủy sinh hóa, dẫn đến tình trạng suy giảm nghiêm trọng nồng độ DO trong nước sông.


Kết luận

Luận án của tác giả Thanh Thuy Nguyen đã giải quyết thành công bài toán cân bằng giữa độ chính xác vật lý và hiệu quả tính toán trong mô hình hóa chất lượng nước sông thông qua việc phát triển nền tảng mô hình khái niệm CORIWAQ-RS. Công trình đã lượng hóa rõ nét ảnh hưởng của biến thiên nhám lòng sông theo mùa do thực vật thủy sinh, đánh giá toàn diện độ bất định cấu trúc giữa các phần mềm thương mại (MIKE 11, InfoWorks RS, InfoWorks ICM) và định lượng đóng góp của độ bất định dữ liệu đầu vào. Kết quả nghiên cứu cung cấp cơ sở khoa học tin cậy và công cụ tính toán hiệu quả phục vụ công tác quy hoạch, quản lý chất lượng nước lưu vực sông và đánh giá rủi ro môi trường dưới tác động của biến đổi khí hậu.