Tổng quan về luận án

Bảo vệ môi trường biển và đới ven bờ là thách thức mang tính chiến lược toàn cầu, đặc biệt tại các vùng vịnh kín chịu áp lực kép từ quá trình đô thị hóa nhanh và phát triển công nghiệp khai khoáng. Luận án tiến sĩ kỹ thuật chuyên ngành Công nghệ môi trường nước và nước thải (Mã số: 9520320) của nghiên cứu sinh Trần Đức Dũng, dưới sự hướng dẫn khoa học của GS.NGND. Trần Hiếu Nhuệ và PGS.TS. Trương Văn Bốn tại Trường Đại học Xây dựng (2020), mang tựa đề: "Nghiên cứu, đánh giá một số chỉ tiêu ô nhiễm bằng phương pháp mô hình hóa, làm cơ sở đề xuất một số giải pháp cải thiện chất lượng nước biển ven bờ Vịnh Cửa Lục, Vịnh Hạ Long, Quảng Ninh". Đây là công trình tiên phong thiết lập phương pháp luận định lượng sức chịu tải môi trường và khả năng tiếp nhận chất ô nhiễm của vùng vịnh bán khép kín kết nối trực tiếp với Di sản Thiên nhiên Thế giới Vịnh Hạ Long.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực trạng các công trình trước đây tại khu vực vịnh Cửa Lục – vịnh Hạ Long (như nghiên cứu của JICA, 1999; Hoàng Doanh Sơn và Vũ Văn Thành, 2000; Trần Lưu Khanh và cộng sự, 2006) chủ yếu dừng lại ở việc kiểm kê nguồn thải đơn lẻ, phân tích thống kê số liệu quan trắc gián đoạn hoặc đánh giá thủy lực độc lập, thiếu sự liên kết động lực học giữa quá trình thủy văn - xói mòn lưu vực với sự lan truyền - chuyển hóa sinh hóa các chất ô nhiễm trong môi trường biển ven bờ. Luận án giải quyết triệt để khoảng trống này bằng cách tích hợp chuỗi mô hình không gian đa quy mô SWAT (Soil and Water Assessment Tool) và MIKE 21 (với các phân hệ HD và ECOLAB).

Hệ thống câu hỏi và giả thuyết nghiên cứu của luận án được xác lập chặt chẽ:

  • Câu hỏi nghiên cứu 1 (RQ1): Tải lượng các chất ô nhiễm chính (BOD5, NH4+, NO3-, PO43-, Fe, Mn) phát sinh từ các phân vùng lưu vực đổ vào vịnh Cửa Lục ở hiện trạng (2018) và theo kịch bản quy hoạch phát triển kinh tế - xã hội đến năm 2030 là bao nhiêu?
  • Câu hỏi nghiên cứu 2 (RQ2): Quy luật thủy động lực học dòng triều và cơ chế lan truyền, pha trộn, tự làm sạch các chỉ tiêu ô nhiễm trong khối nước vịnh Cửa Lục diễn ra như thế nào theo các pha triều và theo mùa (mưa - khô)?
  • Câu hỏi nghiên cứu 3 (RQ3): Ngưỡng sức chịu tải môi trường nước vịnh Cửa Lục đối với từng thông số ô nhiễm là bao nhiêu, những chỉ tiêu nào đã vượt ngưỡng tiếp nhận và giải pháp công nghệ - phi công trình nào khả thi để kiểm soát ô nhiễm?
  • Giả thuyết nghiên cứu 1 (H1): Quá trình bồi lắng trầm tích lơ lửng kết hợp với tải lượng ô nhiễm hữu cơ và dinh dưỡng từ lục địa vượt quá khả năng tự làm sạch tự nhiên của vịnh Cửa Lục, gây suy thoái chất lượng nước vùng đệm trước khi thoát ra vịnh Hạ Long.
  • Giả thuyết nghiên cứu 2 (H2): Chế độ nhật triều và hiệu ứng "nút thắt cổ chai" tại eo Cửa Lục tạo ra sự bất đối xứng dòng chảy, làm tích tụ các chất ô nhiễm dạng vết (kim loại nặng Fe, Mn) và chất dinh dưỡng trong mùa khô cao hơn mùa mưa tại các vùng xoáy quẩn ven bờ.
  • Giả thuyết nghiên cứu 3 (H3): Việc kết hợp đồng bộ các mô-đun xử lý sinh học phân tán tiên tiến (AAO-MBBR) cho nước thải đô thị và quy trình hóa lý - lọc màng (Keo tụ - Lắng lamen - Lọc mangan - Nanofilter) cho nước thải mỏ than sẽ đưa tải lượng xả thải về dưới ngưỡng sức chịu tải cho phép của thủy vực đến năm 2030.

Khung lý thuyết tiếp cận của luận án dựa trên mô hình DPSIR (Driving Forces - Pressures - State - Impact - Response) kết hợp Lý thuyết Sức chịu tải môi trường (Environmental Carrying Capacity Theory) và Động học chuyển hóa sinh hóa nước mặt. Luận án xác định phạm vi nghiên cứu bao gồm toàn bộ diện tích mặt nước vịnh Cửa Lục (khoảng 18 km2), hệ thống 5 lưu vực sông phụ cận (sông Diễn Vọng diện tích 258 km2, sông Man diện tích 94,4 km2, sông Trới diện tích 181 km2) và vùng ven bờ vịnh Hạ Long trong bán kính 1,5 km tính từ cầu Bãi Cháy, bao quát chuỗi thời gian từ năm 2011 đến 2018 và dự báo động thái chất lượng nước đến năm 2030.

Literature Review và Positioning

Lịch sử mô hình hóa chất lượng nước mặt đã phát triển qua ba giai đoạn lớn kể từ công trình kinh điển của Streeter và Phelps (1925) về mô hình suy giảm oxy hòa tan trên sông Ohio. Từ các mô hình giải tích 1 chiều trạng thái tĩnh, khoa học môi trường thế giới đã tiến tới các hệ thống mô hình 2D/3D phi tĩnh phức hợp mô phỏng tương tác sinh thái như WASP7 (Ambrose et al.), QUAL2E (Brown & Barnwell, USEPA), EFDC (Hamrick, 1992), Delft3D-WAQ (Deltares) và hệ phần mềm MIKE 21/3 (Danish Hydraulic Institute - DHI). Kannel và cộng sự (2007) đã chứng minh tính ưu việt của việc tích hợp động lực học dòng chảy với các mô-đun sinh thái trong đánh giá tải lượng tối đa cho phép ngày (TMDL).

Trong bức tranh nghiên cứu ô nhiễm vùng ven bờ, các báo cáo của Tổ chức Nhóm Chuyên gia về Khía cạnh Khoa học của Bảo vệ Môi trường Biển (GESAMP, 1990; 2001) và Chương trình Môi trường Liên Hợp Quốc (UNEP/GPA, 2006) chỉ rõ 60 - 80% tải lượng ô nhiễm biển có nguồn gốc từ lục địa thông qua 6 con đường dẫn nhập chính. Tuy nhiên, các cuộc tranh luận học thuật lớn vẫn tồn tại xoay quanh hai luồng quan điểm:

  1. Trường phái Thủy văn - Trầm tích: Cho rằng quá trình xói mòn cơ học, rửa trôi đất và tải lượng chất rắn lơ lửng (TSS) từ các hoạt động khai khoáng, san lấp ven biển là nhân tố chi phối chính làm biến đổi nền đáy và suy giảm chất lượng nước ven bờ (điển hình như các đánh giá của Nguyễn Đức Cự, 1998 về suy giảm diện tích bãi triều gây mất năng lực tự làm sạch).
  2. Trường phái Sinh hóa - Phú dưỡng: Lại nhấn mạnh rằng ô nhiễm hữu cơ hòa tan (BOD, COD) và mất cân bằng tỷ lệ N:P do nước thải sinh hoạt đô thị và nuôi trồng thủy sản chưa qua xử lý mới là mối đe dọa sinh thái cốt tử gây hiện tượng phú dưỡng, giảm DO và tạo mầm bệnh (tương đồng với nhận định của Báo cáo Hiện trạng Môi trường Quốc gia, 2010 và dự báo gia tăng 14% tải lượng nitơ toàn cầu vào năm 2030 của UNEP).

Luận án của Trần Đức Dũng định vị chính xác điểm giao thoa của hai dòng quan điểm trên: trong một thủy vực bán kín như vịnh Cửa Lục, động lực học dòng triều đóng vai trò điều tiết không gian khuếch tán, trong khi tải lượng trầm tích mỏ than và tải lượng dinh dưỡng đô thị cùng tác động cộng hưởng làm suy giảm sức chịu tải.

So sánh với các nghiên cứu quốc tế:

  • Nghiên cứu vịnh Manila và 11 vịnh ven biển Philippines (PEMSEA & BFAR, 2004): Đánh giá tình trạng phú dưỡng do NO3-, NH4+, PO43- và kim loại nặng từ nông nghiệp, nước thải sinh hoạt bằng phương pháp cân bằng khối đơn giản, nhưng thiếu mô hình thủy động lực học 2D độ phân giải cao để mô phỏng tương tác triều - sông chi tiết theo từng mùa.
  • Nghiên cứu Biển Baltic và Địa Trung Hải theo Khung Chỉ thị Chiến lược Biển Châu Âu (EU Marine Strategy Framework Directive - MSFD): Các nghiên cứu tại đây chỉ ra nồng độ Cadmium, Thủy ngân, Chì vượt ngưỡng sinh thái dù đã kiểm soát nguồn điểm; tuy nhiên các vịnh ven biển nhiệt đới gió mùa có biên độ triều lớn (như vịnh Cửa Lục biên độ triều cực đại đạt 4,7 m) sở hữu cơ chế xáo trộn hoàn toàn khác biệt với các biển kín ít triều tại châu Âu.

Công trình của tác giả đã phát triển vượt bậc so với các nghiên cứu trước đó tại vịnh Cửa Lục (như nghiên cứu năm 1997 của Viện Tài nguyên và Môi trường biển; Hoàng Việt và cộng sự, 1998; Trần Lưu Khanh, 2006) nhờ việc lượng hóa chính xác tỷ lệ phần trăm vượt tải của từng chất ô nhiễm dựa trên Tiêu chuẩn và Quy chuẩn kỹ thuật Quốc gia (QCVN 10-MT:2015/BTNMT và QCVN 08-MT:2015/BTNMT).

Đó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 Sức chịu tải Thủy vực Ven bờ (Coastal Environmental Carrying Capacity) thông qua việc xây dựng mối quan hệ hàm số giữa tải lượng phát thải từ lưu vực sông, động lực học dòng chảy triều không ổn định và tốc độ phân hủy/chuyển hóa nội tại của 06 thông số chỉ thị: nhu cầu oxy hóa sinh hóa (BOD5), amoni (NH4+), nitrat (NO3-), photphat (PO43-), sắt tổng số (Fe) và mangan (Mn).

Mô hình lý thuyết thiết lập phương trình cân bằng khối lượng bảo toàn động học đối với chất ô nhiễm trong thể tích kiểm soát $V$ của vịnh kín:

$$\frac{\partial (hC)}{\partial t} + \frac{\partial (uhC)}{\partial x} + \frac{\partial (vhC)}{\partial y} = \frac{\partial}{\partial x}\left(h D_x \frac{\partial C}{\partial x}\right) + \frac{\partial}{\partial y}\left(h D_y \frac{\partial C}{\partial y}\right) + S + \sum R_k$$

Trong đó $h$ là độ sâu cột nước; $u, v$ là các thành phần vận tốc dòng chảy trung bình theo độ sâu; $D_x, D_y$ là hệ số khuếch tán rối ngang; $S$ là tải lượng nguồn thải biên (từ mô hình SWAT và các cống xả trực tiếp); $\sum R_k$ biểu diễn các quá trình động học chuyển hóa sinh hóa (quá trình nitrate hóa phân hủy NH4+, hấp phụ/giải hấp PO43-, oxy hóa kết tủa Fe2+/Mn2+ và phân hủy hiếu khí BOD5).

Đóng góp đột phá về mặt lý thuyết thể hiện ở việc xác định rõ sự chuyển dịch pha (paradigm shift) từ quản lý nồng độ xả thải cuối đường ống (end-of-pipe concentration management) sang kiểm soát tổng lượng xả thải tối đa dựa trên sức chịu tải môi trường (Total Maximum Allocated Load) có xét đến dao động thủy văn theo mùa.

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

Khung phân tích của luận án tích hợp liên hoàn 3 lý thuyết và công cụ nền tảng:

  1. Lý thuyết Cân bằng Thủy văn Lưu vực (Hydrological Water Balance Framework): Sử dụng mô hình SWAT để phân tích chu trình nước trong pha đất, phân chia lưu vực thành các Đơn vị Phản hồi Thủy văn (Hydrological Response Units - HRUs) dựa trên bản đồ thổ nhưỡng, độ dốc và loại hình sử dụng đất. Mô hình SWAT tính toán chính xác dòng chảy mặt, dòng ngầm và tải lượng rửa trôi phân tán từ các bãi thải khai thác than, đất rừng và khu dân cư thuộc huyện Hoành Bồ và TP. Hạ Long.
  2. Lý thuyết Thủy động lực học Vùng nước Nông 2 chiều (2D Shallow Water Equations): Hệ phương trình Saint-Venant 2 chiều giải bằng phương pháp sai phân hữu hạn/phần tử hữu hạn trên lưới phi cấu trúc của mô hình MIKE 21 HD, mô tả chính xác chế độ nhật triều, dòng triều dâng/rút và vận tốc biến thiên qua eo hẹp Bãi Cháy.
  3. Lý thuyết Động thái Sinh thái Thủy vực (Aquatic Ecological Kinetics): Mô-đun ECOLAB giải quyết đồng thời các phương trình vi phân mô tả tốc độ tiêu thụ oxy, phân rã chất hữu cơ bậc một, chu trình chuyển hóa nitơ - photpho và động thái oxy hóa - lắng đọng của kim loại nặng.

Các điều kiện biên và giới hạn lý thuyết được tường minh hóa: mô hình giả định nước xáo trộn hoàn toàn theo chiều đứng (2D depth-averaged), thích hợp cho vùng vịnh nông có độ sâu trung bình dưới 10 m chịu xáo trộn triều mạnh; các phản ứng sinh hóa được chuẩn hóa ở nhiệt độ tham chiếu 20°C với hệ số điều chỉnh nhiệt độ Arrhenius $\theta$.

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

Thiết kế nghiên cứu

Luận án tuân thủ nghiêm ngặt triết lý nghiên cứu thực chứng duy lý (positivism) thông qua việc kết hợp chặt chẽ giữa quan trắc thực địa, phân tích phòng thí nghiệm đạt chuẩn và kỹ thuật mô hình hóa số trị tiên tiến. Thiết kế nghiên cứu phân tầng đa quy mô (multi-level design):

  • Cấp độ 1 - Lưu vực (Basin Level): Quy mô không gian 533,4 km2 gồm 3 lưu vực sông chính (sông Diễn Vọng, sông Man, sông Trới) và các tiểu lưu vực phụ cận. SWAT được thiết lập với dữ liệu DEM địa hình, bản đồ thổ nhưỡng feralit/phù sa/đất mặn ven biển và chuỗi số liệu khí tượng - lượng mưa 1800 - 2400 mm/năm.
  • Cấp độ 2 - Cửa sông & Vịnh biển (Estuarine & Coastal Bay Level): Miền tính toán thủy động lực học và chất lượng nước bao trùm toàn bộ vịnh Cửa Lục và vùng nước mở vịnh Hạ Long.

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý số liệu được triển khai theo các giao thức khoa học chuẩn mực:

  • Khảo sát - Đo đạc thực địa: Triển khai 02 đợt đo đạc đồng bộ chuyên sâu tại hiện trường: Đợt 1 vào mùa khô (tháng 11/2017) và Đợt 2 vào mùa mưa (tháng 5/2018). Mỗi đợt đo được tiến hành lặp lại tại 02 thời điểm đặc trưng: đỉnh triều (nước lớn) và chân triều (nước ròng).
  • Mạng lưới trạm: Đo lưu lượng và mực nước thực tế liên tục tại 02 mặt cắt kiểm soát quan trọng là cầu Bang (sông Diễn Vọng) và eo Cửa Lục (cầu Bãi Cháy); bố trí trạm đo vận tốc dòng chảy tầng mặt và tầng đáy tại trạm giữa vịnh Cửa Lục (Trạm V).
  • Mạng lưới quan trắc chất lượng nước: Lấy mẫu tại các điểm cửa sông chính, vùng trung tâm vịnh, luồng Cái Lân và khu vực giáp ranh vịnh Hạ Long (theo Bảng 2.13 và các Bảng phụ lục PL 2.1, PL 2.2).
  • Phân tích phòng thí nghiệm: Tuân thủ các tiêu chuẩn TCVN và ISO; thực hiện tại Phòng thử nghiệm - Trung tâm Phân tích FPD (TP. Hạ Long) và Phòng Phân tích Chất lượng Môi trường - Công ty CP Kỹ thuật Môi trường Đại Việt.
  • Kiểm định độ tin cậy và chuẩn hóa mô hình: Sử dụng hệ số hiệu quả mô hình Nash-Sutcliffe (NASH/NSE) và hệ số tương quan $R^2$ để đánh giá mức độ khớp giữa đường mực nước, vận tốc và nồng độ chất ô nhiễm thực đo so với kết quả mô phỏng số.

Data và phân tích

Dữ liệu đầu vào và các kiểm tra độ vững (robustness checks) đạt mức độ chi tiết cao:

  • Hiệu chỉnh mô hình thủy lực: Bộ thông số đáy và độ nhám Manning ($n$) được hiệu chỉnh chi tiết (Bảng 3.22); kết quả so sánh đường mực nước và vận tốc dòng chảy thực đo tháng 11/2017 và tháng 5/2018 tại trạm Bãi Cháy và trạm V cho thấy sự trùng khớp rất cao cả về pha và biên độ triều (Hình 3.3 đến Hình 3.6).
  • Hiệu chỉnh mô hình chất lượng nước: Hiệu chuẩn thành công các thông số trong mô-đun ECOLAB (Bảng 3.23), bao gồm hệ số phân hủy BOD ($k_d = 0,1 - 0,25 \text{ ngày}^{-1}$), tốc độ nitrate hóa ($k_n = 0,05 - 0,15 \text{ ngày}^{-1}$), tốc độ lắng trầm tích lơ lửng, và hệ số oxy hóa sắt/mangan. Kết quả mô phỏng độ mặn, BOD5, Fe, Mn được kiểm chứng tương thích cao với dữ liệu thực nghiệm mùa khô và mùa mưa (Hình 3.8 đến Hình 3.15).

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

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

Kết quả mô phỏng và tính toán tải lượng mang lại 5 phát hiện cốt lõi mang tính định lượng xác thực:

  1. Bức tranh vượt tải nghiêm trọng tại hiện trạng (Năm 2018): Đánh giá sức tải môi trường cho thấy có 04/06 chỉ tiêu chất lượng nước đã vượt ngưỡng chịu tải tự nhiên của vịnh Cửa Lục theo quy chuẩn nước biển ven bờ (QCVN 10-MT:2015/BTNMT):
    • Nhu cầu oxy sinh hóa ($BOD_5$): Đạt mức tải 114,41% (vượt tải 14,41%).
    • Amoni ($NH_4^+$): Vượt tải cực kỳ nghiêm trọng, đạt mức tải lên tới 528,65% (gấp hơn 5,2 lần sức chịu tải của vịnh).
    • Photphat ($PO_4^{3-}$): Đạt mức tải 170,81% (vượt tải 70,81%).
    • Sắt tổng số ($Fe$): Đạt mức tải 241,14% (vượt tải 141,14%).
    • Hai chỉ tiêu còn lại ($NO_3^-$ và $Mn$) vẫn nằm trong giới hạn tiếp nhận cho phép ở thời điểm 2018.
  2. Nhận diện 02 điểm nóng ô nhiễm nghiêm trọng (Critical Hotspots):
    • Khu vực ô nhiễm 1: Khu vực cửa sông Diễn Vọng và vùng tiếp nhận nước thải khai thác than Hòn Gai - Cẩm Phả; đặc trưng bởi hàm lượng chất rắn lơ lửng TSS, Fe và Mn rất cao do nước tháo khô mỏ than và bãi thải đất đá bị rửa trôi (sản lượng nước thải mỏ hàng chục triệu m3/năm).
    • Khu vực ô nhiễm 2: Vùng ven bờ phía Tây và Tây Nam (khu vực đô thị Bãi Cháy, Hà Khánh, Cao Xanh, Cột 5 - Cột 8 và khu vực cảng Cái Lân); chịu tác động tập trung của nước thải sinh hoạt đô thị, dịch vụ du lịch và tàu thuyền cảng biển với nồng độ BOD5, Coliform và Amoni cao đột biến.
  3. Quy luật phân bố thủy động lực học phi đối xứng: Vận tốc dòng triều rút tại eo Cửa Lục cao hơn dòng triều dâng từ 1,5 đến 2 lần; vào kỳ triều cường cao gấp 2,5 - 3 lần kỳ triều kiệt. Vào mùa mưa, lưu lượng lũ từ sông Diễn Vọng (mô đun đỉnh lũ lên đến 51,0 l/s.km2) cộng hưởng với pha triều xuống đẩy nhanh khối nước ô nhiễm ra ngoài vịnh Hạ Long. Ngược lại, vào mùa khô, lượng nước ngọt thượng lưu suy giảm nghiêm trọng khiến nước mặn xâm nhập sâu, vận tốc dòng yếu tạo ra các vùng xoáy tù đọng tại các vũng ven bờ, làm gia tăng thời gian lưu chất ô nhiễm và đẩy nồng độ NH4+, PO43- tích lũy cục bộ lên mức cực đại.
  4. Dự báo bức tranh quá tải môi trường đến năm 2030 (Scenario 2): Nếu không có các công trình can thiệp xử lý nước thải triệt để theo quy hoạch, đến năm 2030, tổng tải lượng phát sinh từ sinh hoạt, du lịch và công nghiệp sẽ tiếp tục tăng mạnh (Bảng 3.20, Bảng 3.21). Đến năm 2030, có 04 thông số dự báo vượt tải nghiêm trọng là $NH_4^+$, $PO_4^{3-}$, $Fe$ và $Mn$. Đặc biệt, Mangan từ trạng thái còn khả năng tiếp nhận năm 2018 sẽ chuyển sang trạng thái quá tải do tích lũy kéo dài từ các nguồn thải công nghiệp mỏ và nhiệt điện.
  5. Cơ chế bồi lắng và xói mòn: Nghiên cứu xác lập mối liên hệ giữa xói mòn các bãi thải than phía Đông lưu vực với sự gia tăng độ dày lớp trầm tích bồi lắng tại lòng vịnh Cửa Lục (Hình 3.22, Hình 3.24), làm giảm dung tích chứa và thu hẹp mặt cắt lưu thông nước của vịnh qua từng năm.

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp mô hình toán tích hợp liên kết lưu vực - vịnh biển kín, hoàn thiện phương pháp tính toán sức tải cho các vùng bờ nhiệt đới gió mùa có chế độ nhật triều biên độ lớn.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình liên kết SWAT - MIKE 21 ECOLAB, có thể chuyển giao và áp dụng trực tiếp cho các hệ thống vũng vịnh ven bờ tương tự tại Việt Nam (như vịnh Cam Ranh, vịnh Đà Nẵng, đầm phá Tam Giang - Cầu Hai).
  • Về giải pháp công trình thực tiễn:
    • Xử lý nước thải sinh hoạt đô thị: Đề xuất mô hình phân tán kết hợp tập trung sử dụng công nghệ vi sinh bùn hoạt tính kết hợp giá thể di động cải tiến AAO-MBBR (Mô-đun FRP-MBBR-100) có hiệu suất khử Nitơ và Photpho vượt trội, phù hợp với các dải đô thị đồi dốc quanh vịnh Cửa Lục.
    • Xử lý nước thải mỏ than hầm lò: Đề xuất chuỗi công nghệ tích hợp Keo tụ - Lắng lamen - Lọc mangan trọng lực - Lọc màng Nanofilter (NF) nhằm loại bỏ triệt để huyền phù than mịn, ion kim loại nặng Fe2+, Mn2+ và sulfate trước khi xả thải ra hệ thống sông Diễn Vọng.
  • Về chính sách và quản lý: Cung cấp luận cứ khoa học để Sở Tài nguyên và Môi trường tỉnh Quảng Ninh và Ban Quản lý Vịnh Hạ Long điều chỉnh quy hoạch phân vùng xả thải, thiết lập hạn ngạch phát thải cho từng lưu vực nhánh, và hoàn thiện trạm quan trắc tự động liên tục tại eo Cửa Lục.

Limitations và Future Research

Luận án thẳng thắn thừa nhận một số hạn chế khoa học do giới hạn nguồn lực và điều kiện kỹ thuật:

  1. Giới hạn chiều không gian mô hình: Việc áp dụng mô hình 2 chiều trung bình theo độ sâu (2D depth-averaged) chưa mô tả được đầy đủ hiện tượng phân tầng mật độ (stratification) và nêm mặn (salt wedge) xảy ra cục bộ tại các hố sâu đáy vịnh vào mùa mưa.
  2. Mật độ chuỗi dữ liệu thực nghiệm: Mới chỉ thực hiện được 02 đợt khảo sát thực tế chuyên sâu (tháng 11/2017 và tháng 5/2018). Mặc dù dữ liệu này đủ độ tin cậy để hiệu chỉnh mô hình nhưng chưa bao quát được toàn diện các biến cố thời tiết cực đoan như bão lớn hay lũ lịch sử.
  3. Phạm vi thông số ô nhiễm: Chưa tích hợp mô phỏng quá trình tích tụ sinh học (bioaccumulation) của kim loại nặng trong chuỗi thức ăn thủy sinh và vi nhựa (microplastics).

Chương trình nghiên cứu tiếp theo (Future Research Agenda) được đề xuất gồm 4 hướng:

  • Hướng 1: Nâng cấp mô hình thủy động lực học lên không gian 3 chiều (MIKE 3D / Delft3D 3D) kết hợp khảo sát chi tiết phân tầng nhiệt - mặn theo mùa.
  • Hướng 2: Ứng dụng công nghệ đồng hóa dữ liệu viễn thám độ phân giải cao (Sentinel-2, Landsat-8, VNREDSat-1) với thuật toán OC2/OC4 để giám sát liên tục trường phân bố Chlorophyll-a và TSS tầng mặt theo thời gian thực.
  • Hướng 3: Đánh giá rủi ro sinh thái và độc học môi trường (Toxicity Identification Evaluation - TIE) đối với hệ sinh thái rừng ngập mặn và rạn san hô vùng lõi Vịnh Hạ Long dưới tác động dài hạn của các kim loại nặng dạng vết.
  • Hướng 4: Thiết lập mô hình kinh tế tuần hoàn tài nguyên nước, nghiên cứu tái sử dụng toàn diện nước thải mỏ than sau xử lý màng NF phục vụ công nghiệp và đô thị tỉnh Quảng Ninh.

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

Luận án tạo ra tác động sâu rộng trên cả phương diện học thuật lẫn thực tiễn quản trị môi trường:

  • Giá trị học thuật: Kết quả luận án đóng góp nguồn dữ liệu định lượng chuẩn xác và mô hình tham chiếu chất lượng cao cho các nghiên cứu hải văn, sinh thái biển ven bờ tại Việt Nam; tạo tiền đề cho các công bố quốc tế về cơ chế trao đổi chất qua eo biển hẹp.
  • Tác động chính sách cấp tỉnh: Cung cấp cơ sở định lượng trực tiếp để UBND tỉnh Quảng Ninh phê duyệt các đề án bảo vệ môi trường, kiểm soát cấp phép xả thải công nghiệp đối với Tập đoàn Công nghiệp Than - Khoáng sản Việt Nam (TKV), đồng thời điều chỉnh Quy hoạch chung TP. Hạ Long mở rộng sau khi sáp nhập huyện Hoành Bồ.
  • Bảo tồn di sản toàn cầu: Đóng góp trực tiếp vào nỗ lực thực thi các khuyến nghị của UNESCO và IUCN trong việc kiểm soát ô nhiễm vùng đệm, bảo đảm tính toàn vẹn của hệ sinh thái Di sản Thiên nhiên Thế giới Vịnh Hạ Long.

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

  • Nghiên cứu sinh & Giới học thuật: Tiếp cận phương pháp luận liên kết mô hình SWAT - MIKE 21 ECOLAB chuẩn mực, khung đánh giá sức tải môi trường chi tiết và bộ số liệu thủy hải văn thực nghiệm quý giá tại vịnh Cửa Lục.
  • Cơ quan Quản lý Nhà nước (Bộ TN&MT, Sở TN&MT Quảng Ninh, BQL Vịnh Hạ Long): Sở hữu công cụ số hóa mô phỏng để dự báo sự cố môi trường, cấp phép xả thải dựa trên hạn ngạch sức tải và xây dựng mạng lưới quan trắc thông minh.
  • Các doanh nghiệp công nghiệp & Khai thác mỏ (TKV, Tổng công ty Đông Bắc): Ứng dụng dây chuyền công nghệ xử lý nước thải mỏ than tối ưu (Keo tụ - Lắng lamen - Lọc mangan - Nanofilter) nhằm tuân thủ tiêu chuẩn xả thải và hướng tới kinh tế tuần hoàn nước.
  • Cộng đồng dân cư & Ngành du lịch địa phương: Hưởng lợi từ chất lượng môi trường nước biển được cải thiện, bảo vệ nguồn lợi thủy sản, giảm thiểu nguy cơ dịch bệnh và nâng cao giá trị dịch vụ du lịch sinh thái bền vững.

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ì và đã mở rộng lý thuyết nào?
Trả lời: Luận án đã mở rộng Lý thuyết Sức chịu tải Môi trường Ven bờ (Coastal Carrying Capacity Theory) bằng cách tích hợp liên hoàn phương trình truyền khối động học 2D với chuỗi phản ứng sinh hóa của 06 chỉ tiêu ô nhiễm đặc trưng trong điều kiện vịnh kín chịu chi phối của nhật triều biên độ lớn. Luận án đã định lượng hóa chính xác tỷ lệ phần trăm vượt tải của từng chất (BOD5: 114,41%, NH4+: 528,65%, PO43-: 170,81%, Fe: 241,14%) thay vì chỉ dừng ở các ước lượng định tính.

2. Điểm đổi mới về mặt phương pháp luận so với các nghiên cứu trước đây là gì?
Trả lời: So với các nghiên cứu trước đây (như nghiên cứu của JICA 1999 chỉ dùng mô hình đơn lẻ; Trần Hữu Tuyên 2005 dùng HSCTM-2D chỉ mô phỏng độ mặn), luận án đã thiết lập quy trình tích hợp đa quy mô hai chiều: mô hình SWAT tính toán tải lượng ô nhiễm từ toàn bộ 533,4 km2 lưu vực đất liền làm điều kiện biên đầu vào liên tục cho mô hình thủy động lực - chất lượng nước MIKE 21 HD/ECOLAB trên mạng lưới ô lưới phi cấu trúc chi tiết hóa toàn bộ vịnh Cửa Lục và vùng nước mở vịnh Hạ Long.

3. Phát hiện nào gây bất ngờ nhất từ kết quả mô hình hóa và thực nghiệm?
Trả lời: Phát hiện bất ngờ nhất là sự vượt tải cực kỳ nghiêm trọng của ion Amoni ($NH_4^+$) lên tới 528,65% sức chịu tải của vịnh ở hiện trạng năm 2018, cùng với sự đảo chiều của thông số Mangan ($Mn$) từ mức an toàn năm 2018 sang trạng thái vượt tải nghiêm trọng trong kịch bản dự báo năm 2030 nếu nguồn nước thải mỏ và công nghiệp tiếp tục gia tăng.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?
Trả lời: Có. Luận án cung cấp đầy đủ các bảng thông số hiệu chỉnh mô hình (Bảng 3.22 cho thủy lực, Bảng 3.23 cho ECOLAB), hệ số nhám Manning, hệ số phân hủy động học, bản đồ phân vùng thổ nhưỡng/sử dụng đất, vị trí tọa độ các trạm đo và toàn bộ chuỗi số liệu quan trắc thực địa tại các phụ lục (PL 2.1 đến PL 5.24), cho phép các nhà khoa học khác tái lập hoàn toàn mô hình.

5. Chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?
Trả lời: Luận án vạch ra lộ trình 10 năm tập trung vào: (1) Nâng cấp hệ thống mô phỏng 3D thời gian thực tích hợp trí tuệ nhân tạo (AI-driven hydrodynamic modeling); (2) Tự động hóa giám sát chất lượng nước qua vệ tinh viễn thám; (3) Mở rộng phạm vi đánh giá rủi ro sinh thái độc học kim loại nặng đến toàn bộ quần thể Di sản Thiên nhiên Thế giới Vịnh Hạ Long và Vịnh Bái Tử Long.

Kết luận

  1. Luận án đã thiết lập thành công hệ thống phương pháp luận mô hình hóa tích hợp đa quy mô (SWAT liên kết MIKE 21 HD/ECOLAB), giải quyết trọn vẹn bài toán đánh giá lan truyền ô nhiễm từ lưu vực lục địa ra vịnh biển kín.
  2. Xác lập bức tranh toàn cảnh hiện trạng chất lượng nước năm 2018 với 04/06 chỉ tiêu vượt ngưỡng chịu tải tự nhiên: $BOD_5$ (114,41%), $NH_4^+$ (528,65%), $PO_4^{3-}$ (170,81%) và $Fe$ (241,14%).
  3. Nhận diện và khoanh vùng chính xác 02 vùng ô nhiễm nghiêm trọng tại cửa sông Diễn Vọng (nguồn khoáng sản mỏ than) và dải ven bờ đô thị - cảng Bãi Cháy - Cái Lân - Cột 5 - Cột 8.
  4. Dự báo chính xác kịch bản chất lượng nước đến năm 2030, cảnh báo sự vượt tải bổ sung của thông số Mangan ($Mn$) và đề xuất ngưỡng giới hạn tải lượng tiếp nhận cho từng tiểu lưu vực.
  5. Đề xuất hệ thống giải pháp công nghệ kỹ thuật khả thi, gồm công nghệ vi sinh AAO-MBBR cho nước thải sinh hoạt đô thị và dây chuyền Keo tụ - Lắng lamen - Lọc mangan - Nanofilter cho nước thải mỏ than hầm lò.
  6. Thiết lập cơ sở khoa học và pháp lý vững chắc phục vụ công tác quy hoạch bảo vệ môi trường, kiểm soát sức chịu tải thủy vực, đóng góp thiết thực vào sự phát triển bền vững và bảo tồn Di sản Thiên nhiên Thế giới Vịnh Hạ Long.