Tổng quan về luận án

Vấn đề suy thoái chất lượng nước và ô nhiễm môi trường biển ven bờ đang là thách thức an ninh sinh thái mang tính toàn cầu, đặc biệt tại các khu vực chịu áp lực phát triển kinh tế đa ngành. Vịnh Cửa Lục (diện tích mặt nước khoảng 18 $\text{km}^2$) thuộc tỉnh Quảng Ninh giữ vai trò là một vụng sinh thái bán kín mang tính chất chuyển tiếp, tiếp nhận toàn bộ mạng lưới thủy văn từ lưu vực đất liền trước khi thông ra vịnh Hạ Long qua eo Cửa Lục hẹp tại vị trí cầu Bãi Cháy. Đây là khu vực chịu tác động tổng hợp phức tạp từ hoạt động khai thác than vùng Hòn Gai – Cẩm Phả, phát triển hệ thống cảng biển nước sâu Cái Lân, các cụm công nghiệp xi măng – nhiệt điện, đô thị hóa nhanh chóng tại thành phố Hạ Long và huyện Hoành Bồ, cùng với hoạt động du lịch biển quy mô lớn.

Khoảng trống nghiên cứu (research gap) cốt lõi mà công trình này giải quyết là sự thiếu hụt phương pháp luận tích hợp giữa mô hình thủy văn lưu vực và mô hình động lực học chất lượng nước 2 chiều để lượng hóa chính xác tải lượng các nguồn phát thải, đánh giá sức chịu tải môi trường (environmental carrying capacity) và khả năng tự làm sạch tự nhiên của một vịnh bán kín có chế độ thủy triều phức tạp. Các nghiên cứu trước đây tại khu vực vịnh Hạ Long và vịnh Cửa Lục (như nghiên cứu của JICA, 1999; Hoàng Việt và cộng sự, 1998; Trần Lưu Khanh, 2006) chủ yếu tiếp cận theo phương pháp quan trắc gián đoạn, thống kê đơn lẻ hoặc phân vùng chỉ số chất lượng nước mà chưa thiết lập được mối quan hệ nhân quả động học giữa quy hoạch không gian kinh tế - xã hội trên lưu vực với sự biến động trường nồng độ các chất ô nhiễm trong không gian và thời gian.

Nghiên cứu được định hướng bởi hệ thống câu hỏi nghiên cứu và giả thuyết khoa học chặt chẽ:

  • Câu hỏi nghiên cứu 1 (RQ1): Tổng tải lượng các chất ô nhiễm chính phát sinh từ lưu vực đổ vào vịnh Cửa Lục ở hiện trạng (2018) và dự báo theo quy hoạch phát triển (2030) là bao nhiêu?
  • Câu hỏi nghiên cứu 2 (RQ2): Quy luật động lực học nào chi phối quá trình lan truyền, khuếch tán và phân bố nồng độ của các chất ô nhiễm hữu cơ, dinh dưỡng và kim loại nặng trong chế độ nhật triều của vịnh?
  • Câu hỏi nghiên cứu 3 (RQ3): Ngưỡng chịu tải tối đa của thủy vực vịnh Cửa Lục đối với từng thông số ô nhiễm là bao nhiêu và những giải pháp công nghệ kỹ thuật - phi công trình nào mang tính khả thi cao nhất để bảo vệ chất lượng nước?
  • Giả thuyết khoa học 1 (H1): Sức chịu tải của vịnh Cửa Lục đối với các chất hữu cơ và dinh dưỡng đã vượt ngưỡng cho phép do sự bất đối xứng giữa tốc độ gia tăng phát thải từ đất liền và khả năng tự làm sạch của chế độ nhật triều.
  • Giả thuyết khoa học 2 (H2): Chế độ thủy động lực theo mùa (mùa mưa và mùa khô) tạo ra các vùng ô nhiễm cục bộ đặc thù, trong đó dòng chảy triều dâng mùa khô đẩy chất ô nhiễm tích tụ sát bờ, còn dòng xả mùa mưa vận chuyển ô nhiễm trực tiếp ra vùng lõi di sản vịnh Hạ Long.

Khung lý thuyết của luận án được xây dựng dựa trên sự liên kết giữa khung phân tích DPSIR (Động lực - Áp lực - Hiện trạng - Tác động - Đáp ứng) và lý thuyết động lực học chất lưu môi trường (Environmental Fluid Dynamics). Luận án xác định và mô phỏng 06 thông số chất lượng nước then chốt: Nhu cầu oxy sinh hóa ($\text{BOD}_5$), Amoni ($\text{NH}_4^+$), Nitrat ($\text{NO}_3^-$), Photphat ($\text{PO}_4^{3-}$), Sắt ($\text{Fe}$) và Mangan ($\text{Mn}$). Đóng góp mang tính đột phá của nghiên cứu là xác lập định lượng mức vượt tải môi trường hiện trạng: $\text{BOD}_5$ vượt 114,41%, $\text{NH}_4^+$ vượt 528,65%, $\text{PO}_4^{3-}$ vượt 170,81%, và $\text{Fe}$ vượt 241,14%, đồng thời dự báo bức tranh tải lượng và diễn biến ô nhiễm đến năm 2030, làm cơ sở khoa học thiết kế các trạm xử lý nước thải công nghệ tiên tiến FRP-MBBR và cụm xử lý nước thải mỏ than chuyên dụng.

Literature Review và Positioning

Nghiên cứu về chất lượng nước mặt và ô nhiễm biển ven bờ đã phát triển qua gần một thế kỷ, khởi nguồn từ mô hình kinh điển một chiều của Streeter và Phelps (1925) về phương trình suy giảm oxy hòa tan trên sông Ohio. Trải qua các giai đoạn hoàn thiện thuật toán số, các công cụ mô phỏng chất lượng nước đã nâng cấp từ mô hình một chiều, trạng thái tĩnh sang các hệ thống mô hình động lực học không gian hai chiều và ba chiều tích hợp như WASP, QUAL2E, CE-QUAL-W2, EFDC và bộ công cụ MIKE của Viện Thủy lực Đan Mạch (DHI). Về mặt thủy văn lưu vực, mô hình SWAT (Soil and Water Assessment Tool) do Jeff Arnold và Srinivasan (1998) phát triển tại Bộ Nông nghiệp Hoa Kỳ (USDA) đã trở thành công cụ tiêu chuẩn toàn cầu để lượng hóa dòng chảy mặt, xói mòn đất và tải lượng chất ô nhiễm từ các nguồn diện nông - lâm nghiệp và nhân sinh.

Trong bức tranh học thuật quốc tế và khu vực, hai luồng quan điểm đối lập đã hình thành xoay quanh cơ chế suy thoái chất lượng nước biển ven bờ:

  • Luồng quan điểm thứ nhất: Nhấn mạnh áp lực từ các hoạt động hàng hải, nạo vét luồng tàu và nuôi trồng thủy sản tại chỗ trên biển là nguồn phát tán ô nhiễm tức thời và trực tiếp nhất (GESAMP, 1990; BFAR, 2004).
  • Luồng quan điểm thứ hai: Cho rằng nguồn ô nhiễm có nguồn gốc lục địa (Land-based sources) thông qua hệ thống lưu vực sông mới đóng vai trò quyết định, chiếm từ 50% đến 80% tổng tải lượng chất ô nhiễm đưa ra đại dương (UNEP, GPA, 2006).

Một cuộc tranh luận khoa học khác liên quan trực tiếp đến vịnh Cửa Lục là mức độ ảnh hưởng của hàm lượng tổng chất rắn lơ lửng (TSS) và kim loại nặng. Báo cáo đánh giá ban đầu năm 1997 cho rằng TSS rửa trôi từ lưu vực chỉ ở mức thấp (1 – 2 mg/l tại Cửa Lục và Bãi Cháy) và không phải nguồn ô nhiễm chính cho vịnh Hạ Long. Trái lại, nghiên cứu của Hoàng Việt và cộng sự (1998) đã chỉ ra sự biến động bất thường với nồng độ TSS tăng vọt lên 45,7 – 98,2 mg/l, vượt xa quy chuẩn môi trường do quá trình bạt đồi khai thác than và san lấp mặt bằng đô thị.

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

  • Nghiên cứu 12 vịnh ven biển tại Philippines: Báo cáo của PEMSEA và Cục Nghề cá Philippines (BFAR, 2004) tại các vịnh như vịnh Manila, vịnh San Miguel, vịnh Ragay đã chỉ ra hiện tượng phú dưỡng hóa nghiêm trọng do nồng độ các chất dinh dưỡng ($\text{NO}_3^-$, $\text{PO}_4^{3-}$) và thuốc bảo vệ thực vật tăng cao từ nông nghiệp, kết hợp quá tải kim loại nặng từ khai thác khoáng sản ven bờ.
  • Nghiên cứu các vùng biển ven bờ châu Âu: Đánh giá của Cơ quan Môi trường Châu Âu (EEA) và GESAMP tại biển Baltic và Địa Trung Hải cho thấy 60% – 80% lượng nước thải ven bờ không qua xử lý đạt chuẩn dẫn đến tích lũy nitơ, photpho và các chất ô nhiễm hữu cơ bền (POPs), đe dọa trực tiếp các rạn san hô và sinh vật đáy.

Vị trí học thuật của luận án được xác lập rõ nét: Khắc phục hạn chế của các nghiên cứu trước đây tại Việt Nam vốn chỉ áp dụng đơn lẻ mô hình thủy lực ven biển hoặc quan trắc hóa học tĩnh, luận án tiên phong thiết lập chuỗi mô hình ghép đôi đồng bộ (Coupled modeling chain) giữa mô hình thủy văn lưu vực SWAT và mô hình thủy động lực - sinh thái chất lượng nước 2 chiều MIKE 21 HD/Ecolab. Cách tiếp cận này giúp bao quát từ nguồn phát sinh trên lưu vực $533,4\text{ km}^2$ của 3 hệ thống sông chính (sông Diễn Vọng, sông Trới, sông Man) đến sự biến đổi sinh hóa trong không gian $18\text{ km}^2$ của vùng nước vịnh Cửa Lục.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đóng góp quan trọng vào lý thuyết sức chịu tải môi trường nước (Water Environmental Carrying Capacity Theory) và lý thuyết lan truyền đối lưu - khuếch tán chất ô nhiễm (Advection-Diffusion-Reaction Theory) trong các thủy vực bán kín chịu chế độ nhật triều biến động lớn. Bằng việc lượng hóa chính xác cơ chế tương tác giữa dòng chảy sông nước ngọt từ lưu vực đồi núi dốc và khối nước triều mặn từ biển, luận án đã mở rộng mô hình cân bằng khối lượng kinh điển sang trạng thái biên động lực học phi tuyến tính.

Mô hình lý thuyết của luận án thiết lập các mệnh đề khoa học:

  • Mệnh đề 1 ($P_1$): Sức chịu tải môi trường của một vịnh bán kín tỷ lệ nghịch với thời gian lưu nước (residence time) và phụ thuộc phi đối xứng vào chu kỳ triều cường – triều kiệt theo mùa.
  • Mệnh đề 2 ($P_2$): Sự kết hợp giữa quá trình khoáng hóa chất hữu cơ và lắng đọng phù sa lơ lửng tại vùng ranh giới mặn - ngọt tạo ra các điểm dị thường nồng độ (concentration anomalies) kim loại nặng Fe, Mn độc lập với khoảng cách đến nguồn thải điểm.
       +-----------------------------------------------------------+
       |               KHUNG PHÂN TÍCH TỔNG THỂ (DPSIR)            |
       +-----------------------------------------------------------+
                                     |
       +-----------------------------v-----------------------------+
       | ĐỘNG LỰC (DRIVING FORCES):                                |
       | - Khai thác than Hòn Gai - Cẩm Phả                        |
       | - Đô thị hóa TP. Hạ Long, Huyện Hoành Bồ                  |
       | - Du lịch Bãi Cháy, Cảng Cái Lân, Nuôi trồng thủy sản     |
       +-----------------------------+-----------------------------+
                                     |
       +-----------------------------v-----------------------------+
       | ÁP LỰC (PRESSURE) - TẢI LƯỢNG Ô NHIỄM:                    |
       | - Mô hình lưu vực SWAT: Mô phỏng dòng chảy mặt & xói mòn  |
       | - Kiểm kê phát thải: Sinh hoạt, Công nghiệp mỏ, Chăn nuôi |
       +-----------------------------+-----------------------------+
                                     |
       +-----------------------------v-----------------------------+
       | HIỆN TRẠNG (STATE) & TÁC ĐỘNG (IMPACT):                   |
       | - Mô hình MIKE 21 HD: Trường vận tốc dòng chảy & mực nước |
       | - Module MIKE 21 Ecolab: Mô phỏng BOD5, NH4+, NO3-,       |
       |   PO43-, Fe, Mn (Đánh giá vượt tải môi trường)            |
       +-----------------------------+-----------------------------+
                                     |
       +-----------------------------v-----------------------------+
       | ĐÁP ỨNG (RESPONSE) - GIẢI PHÁP CẢI THIỆN:                 |
       | - Công trình: Module FRP-MBBR-100 & Cụm xử lý mỏ than NF  |
       | - Phi công trình: Phân vùng sức tải, kiểm toán nguồn thải |
       +-----------------------------------------------------------+

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 nền tảng lý thuyết và công cụ phân tích:

  1. Khung quản lý môi trường DPSIR: Đóng vai trò là khung phân tích logic định hướng toàn bộ chuỗi đánh giá nguyên nhân - hệ quả từ phát triển kinh tế đến suy thoái thủy vực.
  2. Mô hình thủy văn bán phân bố SWAT: Chia tách lưu vực thành các Đơn vị phản hồi thủy văn (Hydrologic Response Units - HRU) dựa trên sự đồng nhất về loại đất, độ dốc và thảm phủ thực vật, tính toán chính xác dòng chảy kiệt và đỉnh lũ từ lưu vực sông Diễn Vọng (diện tích lưu vực $258\text{ km}^2$, độ dốc 17,8%), sông Trới ($181\text{ km}^2$, độ dốc 14,8%), và sông Man ($94,4\text{ km}^2$, độ dốc 21,7%).
  3. Hệ thống phương trình vi phân đạo hàm riêng thủy động lực và chất lượng nước 2 chiều nông (2D Shallow Water Equations & Ecolab): Mô tả các quá trình đối lưu, khuếch tán rối, lắng đọng, quang hợp, nitrat hóa, phản nitrat hóa và phân hủy sinh hóa của các cấu tử ô nhiễm.

Điều kiện biên (Boundary Conditions) được xác lập khắt khe: Biên ngoài phía biển tại eo Cửa Lục (Cầu Bãi Cháy) tiếp nhận dao động triều thực đo; biên sông phía thượng lưu tiếp nhận chuỗi lưu lượng và nồng độ chất ô nhiễm từ mô phỏng SWAT và dữ liệu quan trắc mùa khô (tháng 11/2017) và mùa mưa (tháng 5/2018).

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ủ triết lý khoa học thực chứng (Positivism) kết hợp mô hình hóa số tất định (Deterministic numerical modeling). Thiết kế nghiên cứu đa tầng (Multi-level design) bao gồm:

  • Tầng 1: Thu thập, chuẩn hóa dữ liệu viễn thám, GIS, địa hình đáy biển, bản đồ thổ nhưỡng và chuỗi số liệu khí tượng - thủy văn nhiều năm.
  • Tầng 2: Khảo sát thực địa đồng bộ thủy hải văn và chất lượng nước theo chu kỳ triều.
  • Tầng 3: Xây dựng lưới tính toán số trị, hiệu chuẩn (calibration) và kiểm định (validation) mô hình thủy lực và chất lượng nước.
  • Tầng 4: Xây dựng các kịch bản mô phỏng hiện trạng (2018) và dự báo phát triển (2030) để tính toán khả năng tiếp nhận tải lượng ô nhiễm.

Quy trình nghiên cứu rigorous

Quy trình thu thập và phân tích dữ liệu được tiến hành với độ chuẩn xác học thuật cao:

  • Quan trắc hiện trường: Thực hiện 02 đợt đo đạc chi tiết đại diện cho 2 mùa tương phản: Đợt 1 vào tháng 11/2017 (mùa khô) và Đợt 2 vào tháng 5/2018 (mùa mưa). Mỗi đợt tiến hành đo đạc 02 lần tương ứng với hai pha triều đặc trưng: đỉnh triều (nước lớn) và chân triều (nước ròng).
  • Mạng lưới quan trắc: Thiết lập 2 mặt cắt thủy hải văn chính tại cầu Bang (sông Diễn Vọng) và cầu Bãi Cháy (eo Cửa Lục) để đo lưu lượng, mực nước và vận tốc dòng chảy tầng mặt - tầng đáy; cùng hệ thống các điểm lấy mẫu chất lượng nước phân bố đều từ thượng lưu các cửa sông đến tâm vịnh Cửa Lục và vùng đệm vịnh Hạ Long (bán kính 1,5 km tính từ cầu Bãi Cháy).
  • Phân tích phòng thí nghiệm: Toàn bộ mẫu nước biển và nước sông được bảo quản và phân tích theo Tiêu chuẩn Quốc gia (TCVN) và Tiêu chuẩn Quốc tế (Standard Methods for the Examination of Water and Wastewater) tại Trung tâm Phân tích FPD Quảng Ninh và Phòng Phân tích Môi trường Công ty CP Kỹ thuật Môi trường Đại Việt.
  • Tam giác đạc (Triangulation): Dữ liệu được đối chiếu chéo giữa chuỗi quan trắc thực đo, dữ liệu quan trắc định kỳ của Sở Tài nguyên và Môi trường tỉnh Quảng Ninh giai đoạn 2011–2018, và ảnh vệ tinh viễn thám độ phân giải trung bình/cao (Landsat-8, VNREDSat-1).

Data và phân tích

+-----------------------------------------------------------------------------------------+
|                  ĐẶC TRƯNG HÌNH THÁI VÀ THỦY VĂN LƯU VỰC VỊNH CỬA LỤC                   |
+---------------+-------------------+----------------+-------------+----------------------+
| Sông chính    | Diện tích LV (km²)| Chiều dài (km) | Độ dốc (%)  | Lưu lượng Qp 50% (m³/s) |
+---------------+-------------------+----------------+-------------+----------------------+
| Diễn Vọng     | 258,0             | 32,0           | 17,8        | 8,68                 |
| Trới          | 181,0             | 31,0           | 14,8        | 4,69                 |
| Man           | 94,4              | 22,5           | 21,7        | 5,04                 |
+---------------+-------------------+----------------+-------------+----------------------+
  • Mô hình SWAT: Thiết lập dữ liệu đầu vào gồm mô hình số độ cao (DEM), bản đồ sử dụng đất, bản đồ đất feralit, đất dốc tụ và đất mặn ven bờ kết hợp dữ liệu khí tượng mưa ngày (1800 – 2400 mm/năm) từ các trạm Bãi Cháy, Hoành Bồ.
  • Mô hình MIKE 21 HD & Ecolab: Miền tính toán bao trùm toàn bộ vịnh Cửa Lục và vùng nước phụ cận vịnh Hạ Long với hệ thống lưới phi cấu trúc (Flexible Mesh) tối ưu hóa độ phân giải tại các luồng lạch hẹp và cửa sông.
  • Hiệu chuẩn và kiểm định: Đánh giá độ tin cậy của mô hình bằng hệ số hiệu quả mô hình Nash-Sutcliffe (NASH), sai số tuyệt đối trung bình (MAE) và hệ số tương quan ($R^2$). Kết quả so sánh giữa đường mực nước, vận tốc dòng chảy, độ mặn và nồng độ $\text{BOD}_5$, $\text{Fe}$, $\text{Mn}$ tính toán với số liệu thực đo tại trạm Bãi Cháy cho thấy độ trùng khớp cao, phản ánh chính xác quy luật thủy động lực học của khu vực.

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

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

  1. Sự vượt tải nghiêm trọng của 4/6 chỉ tiêu chất lượng nước (Hiện trạng 2018): Kết quả tính toán sức chịu tải môi trường nước vịnh Cửa Lục cho thấy khả năng tiếp nhận của thủy vực đã bị phá vỡ nghiêm trọng ở 4 thông số:

    • Nhu cầu oxy sinh hóa ($\text{BOD}_5$): Vượt tải 114,41%
    • Amoni ($\text{NH}_4^+$): Vượt tải đặc biệt nghiêm trọng 528,65%
    • Photphat ($\text{PO}_4^{3-}$): Vượt tải 170,81%
    • Tổng Sắt ($\text{Fe}$): Vượt tải 241,14%
    • Hai chỉ tiêu còn nằm trong ngưỡng chịu tải là $\text{NO}_3^-$ và $\text{Mn}$.

  2. Nhận diện 02 vùng ô nhiễm môi trường nghiêm trọng cục bộ:

    • Khu vực ô nhiễm nghiêm trọng 1: Vùng nước phía Đông và Đông Bắc vịnh Cửa Lục, nơi tiếp nhận lưu lượng nước thải mỏ than qua sông Diễn Vọng và các suối thoát nước khu vực Hòn Gai – Cẩm Phả. Khu vực này đặc trưng bởi nồng độ tổng sắt ($\text{Fe}$), mangan ($\text{Mn}$) và chất rắn lơ lửng rất cao.
    • Khu vực ô nhiễm nghiêm trọng 2: Vùng ven bờ đô thị Bãi Cháy, Cái Lân và khu vực cửa sông Trới, nơi tập trung nước thải sinh hoạt đô thị chưa qua xử lý triệt để, phát sinh nồng độ $\text{BOD}_5$ và $\text{NH}_4^+$ vượt quy chuẩn QCVN 10-MT:2015/BTNMT nhiều lần.

  3. Quy luật phân bố nồng độ ô nhiễm theo cơ chế mùa và pha triều:

    • Mùa khô: Dòng chảy sông kiệt (lưu lượng sông Diễn Vọng chỉ đạt 1,11 – 2,69 $\text{m}^3/\text{s}$), năng lượng thủy triều từ vịnh Hạ Long chiếm ưu thế tuyệt đối. Dòng triều dâng ép các khối nước ô nhiễm vào sát bờ và tích tụ trong các vụng kín phía trong vịnh, làm gia tăng thời gian lưu chất ô nhiễm.
    • Mùa mưa: Lưu lượng sông tăng đột biến (sông Diễn Vọng đạt 16,7 – 26,5 $\text{m}^3/\text{s}$, sông Trới đạt 8,3 – 12,9 $\text{m}^3/\text{s}$). Dòng xả lũ kết hợp với pha triều rút tạo ra động năng cực đại hướng ra eo Bãi Cháy, vận chuyển một khối lượng khổng lồ bùn cát lơ lửng, chất hữu cơ và kim loại nặng thoát trực tiếp ra vịnh Hạ Long, đe dọa vùng lõi di sản.

  4. Dự báo nguy cơ suy thoái chất lượng nước đế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, 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 gia tăng. Khả năng tự làm sạch của vịnh Cửa Lục bị thu hẹp đáng kể, dẫn đến 4 thông số $\text{NH}_4^+$, $\text{PO}_4^{3-}$, $\text{Fe}$ và $\text{Mn}$ đều rơi vào tình trạng mất kiểm soát về sức chịu tải.

+-----------------------------------------------------------------------------------------+
|                  ĐÁNH GIÁ SỨC TẢI VÀ KHẢ NĂNG TIẾP NHẬN TẠI VỊNH CỬA LỤC                |
+--------------------+-----------------------+---------------------+----------------------+
| Thông số ô nhiễm   | Trạng thái tải (2018) | Mức độ vượt tải (%) | Dự báo nguy cơ 2030  |
+--------------------+-----------------------+---------------------+----------------------+
| BOD₅               | Vượt ngưỡng           | +114,41%            | Nguy cơ cao          |
| NH₄⁺               | Vượt nghiêm trọng     | +528,65%            | Rất nghiêm trọng     |
| NO₃⁻               | Còn sức tải           | Trong ngưỡng        | Cận ngưỡng           |
| PO₄³⁻              | Vượt ngưỡng           | +170,81%            | Vượt ngưỡng          |
| Fe                 | Vượt nghiêm trọng     | +241,14%            | Rất nghiêm trọng     |
| Mn                 | Còn sức tải           | Trong ngưỡng        | Vượt ngưỡng          |
+--------------------+-----------------------+---------------------+----------------------+

Implications đa chiều

  • Về mặt lý luận: Khẳng định vai trò của mô hình ghép đôi SWAT - MIKE 21 như một công cụ chuẩn mực để phân tích chính xác mối quan hệ giữa các hoạt động nhân sinh trên lưu vực và phản ứng động lực học của thủy vực ven biển.
  • Về mặt phương pháp: Cung cấp quy trình định lượng sức chịu tải môi trường chi tiết cho các vùng vịnh kín ven bờ có chế độ nhật triều tại Việt Nam.
  • Về mặt thực tiễn và kỹ thuật:
    • Xử lý nước thải sinh hoạt đô thị và du lịch: Đề xuất mô hình kết hợp giữa các nhà máy xử lý nước thải tập trung quy mô lớn và các trạm xử lý phân tán tại nguồn sử dụng công nghệ sinh học tích hợp giá thể vi sinh di động FRP-MBBR-100 (công nghệ AAO-MBBR) cho các khu đô thị mới ven vịnh và các cơ sở dịch vụ du lịch.
    • Xử lý nước thải hầm lò ngành than: Thiết kế hệ thống công nghệ chuyên dụng bao gồm các công đoạn: Keo tụ tạo bông $\rightarrow$ Lắng lamen tốc độ cao $\rightarrow$ Cột lọc khử sắt, mangan $\rightarrow$ Hệ thống lọc màng Nanofilter (NF) để loại bỏ triệt để ion kim loại nặng và cặn than lơ lửng trước khi xả ra lưu vực sông Diễn Vọng.
  • Về mặt chính sách: Cung cấp cơ sở khoa học để UBND tỉnh Quảng Ninh điều chỉnh Quy hoạch bảo vệ môi trường, kiểm soát chặt chẽ quy hoạch lấn biển đô thị hóa, và thiết lập hạn ngạch phát thải (Total Maximum Daily Loads - TMDL) cho từng phân lưu vực.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận các giới hạn nghiên cứu:

  1. Giới hạn về chiều không gian mô hình: Luận án áp dụng mô hình thủy động lực học 2 chiều dòng nông tích hợp theo chiều sâu (2D depth-averaged). Mô hình này giả định khối nước đồng nhất theo phương thẳng đứng, do đó chưa mô tả trọn vẹn hiện tượng phân tầng mật độ (stratification) mặn - nhiệt vào thời kỳ giao thời giữa mùa mưa và mùa khô tại các luồng lạch sâu.
  2. Giới hạn về chuỗi số liệu quan trắc chuyên sâu: Do điều kiện kinh phí và kỹ thuật, các đợt khảo sát đo đạc thực tế chỉ tập trung vào 02 đợt đo chính (tháng 11/2017 và tháng 5/2018), chưa thể quan trắc liên tục 24/7 trong nhiều năm liền để ghi nhận các hiện tượng thời tiết cực đoan như bão lớn hay triều dâng dị thường.
  3. Phạm vi chất ô nhiễm: Chưa tích hợp mô phỏng sự lan truyền của ô nhiễm hạt vi nhựa (microplastics) và các hợp chất hữu cơ khó phân hủy (POPs) từ hoạt động tàu thuyền.

Định hướng nghiên cứu tiếp theo (Future Research Agenda):

  • Nâng cấp khung mô hình lên 3 chiều (MIKE 3 HD/Ecolab) để nắm bắt đầy đủ cấu trúc dòng chảy phân tầng thẳng đứng và quá trình trao đổi bùn cát đáy tại eo Bãi Cháy.
  • Tích hợp công nghệ viễn thám độ phân giải cao thời gian thực (Sentinel-2, VNREDSat-1) với các thuật toán học máy (Machine Learning) để hiệu chỉnh trực tuyến các tham số mô hình chất lượng nước.
  • Nghiên cứu cơ chế tích lũy sinh học và lan truyền kim loại nặng trong chuỗi thức ăn thủy sản tại vùng đệm vịnh Hạ Long.

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

  • Tác động học thuật (Academic Impact): Công trình là tài liệu tham khảo mẫu mực về kỹ thuật mô hình hóa môi trường nước ven biển, đóng góp trực tiếp vào các chương trình đào tạo sau đại học ngành Kỹ thuật Môi trường và Động lực học Thủy vực.
  • Chuyển đổi ngành công nghiệp (Industry Transformation): Định hướng Tập đoàn Công nghiệp Than - Khoáng sản Việt Nam (TKV) và Tổng công ty Đông Bắc đầu tư nâng cấp đồng bộ các trạm xử lý nước thải mỏ hầm lò và lộ thiên theo quy chuẩn xả thải loại A, chấm dứt việc xả nước thải có hàm lượng sắt và cặn than cao ra sông Diễn Vọng.
  • Định hình chính sách quản lý nhà nước (Policy Influence): Là luận cứ khoa học để Sở Tài nguyên và Môi trường tỉnh Quảng Ninh ban hành các quy định nghiêm ngặt về đánh giá tác động môi trường các dự án đô thị lấn biển, đồng thời hoàn thiện hệ thống mạng lưới quan trắc môi trường tự động liên tục xung quanh vịnh Cửa Lục và vịnh Hạ Long.
  • Bảo tồn di sản và lợi ích xã hội (Societal & Conservation Benefits): Bảo vệ trực tiếp hệ sinh thái biển, rạn san hô và rừng ngập mặn vịnh Hạ Long – Di sản Thiên nhiên Thế giới, nâng cao giá trị thương hiệu du lịch quốc tế và bảo đảm sinh kế bền vững cho cộng đồng ngư dân nuôi trồng thủy sản ven bờ.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận phương pháp luận ghép nối mô hình SWAT và MIKE 21 Ecolab cùng cơ sở dữ liệu thủy văn - chất lượng nước chi tiết của khu vực ven biển Đông Bắc Việt Nam.
  • Cơ quan Quản lý Nhà nước (Bộ TN&MT, UBND tỉnh Quảng Ninh, Ban Quản lý Vịnh Hạ Long): Sở hữu công cụ dự báo và bản đồ sức chịu tải môi trường để thẩm định quy hoạch đô thị, phân vùng xả thải và ra quyết định ứng phó sự cố ô nhiễm môi trường biển.
  • Khối Doanh nghiệp Công nghiệp và Khai khoáng: Nắm bắt các thông số thiết kế công nghệ xử lý nước thải tiên tiến (AAO-MBBR, lắng lọc màng Nanofilter) để tối ưu hóa chi phí đầu tư bảo vệ môi trường.
  • Cộng đồng dân cư và Khách du lịch: Hưởng lợi từ chất lượng môi trường nước biển được cải thiện, giảm thiểu rủi ro sức khỏe cộng đồng và bảo tồn cảnh quan thiên nhiên.

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ì?

Đó là việc tích hợp thành công lý thuyết sức chịu tải môi trường nước với lý thuyết động lực học lan truyền vật chất trong vịnh bán kín có chế độ nhật triều. Luận án đã mở rộng mô hình Streeter-Phelps truyền thống bằng cách kết hợp khung phân tích DPSIR với hệ phương trình đạo hàm riêng mô tả quá trình vận chuyển, biến đổi hóa - sinh 2 chiều nông của 06 cấu tử ô nhiễm, giải quyết bài toán định lượng khả năng tiếp nhận chất ô nhiễm của thủy vực chịu áp lực đa nguồn thải từ lưu vực đồi núi dốc.

2. Sự đổi mới về mặt phương pháp luận so với các nghiên cứu trước đây thể hiện như nào?

So với các nghiên cứu trước đây tại vịnh Hạ Long vốn chỉ dựa trên quan trắc điểm gián đoạn (JICA, 1999) hoặc mô phỏng thủy lực đơn lẻ (Hoàng Doanh Sơn, 2000), luận án đã thiết lập quy trình mô hình hóa hai tầng liên hoàn: Sử dụng mô hình SWAT để mô phỏng chính xác dòng chảy mặt và tải lượng ô nhiễm từ toàn bộ lưu vực $533,4\text{ km}^2$, sau đó chuyển tiếp làm điều kiện biên động học cho mô hình MIKE 21 HD/Ecolab để mô phỏng chi tiết trường nồng độ và sức tải của vịnh Cửa Lục.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì?

Phát hiện bất ngờ nhất là nồng độ Amoni ($\text{NH}_4^+$) tại vịnh Cửa Lục đã vượt ngưỡng chịu tải lên đến 528,65%, mức vượt tải cao nhất trong tất cả các chỉ tiêu nghiên cứu. Nguyên nhân cốt lõi không chỉ đến từ nước thải sinh hoạt đô thị chưa xử lý mà còn do sự phân hủy yếm khí của thảm thực vật rừng ngập mặn bị vùi lấp trong các dự án san lấp lấn biển, kết hợp với hiện tượng tích tụ chất dinh dưỡng trong bùn đáy tại các vụng trũng kín gió.

4. Luận án có cung cấp quy trình nhân rộng (Replication Protocol) không?

Có. Luận án cung cấp chi tiết toàn bộ quy trình thiết lập mô hình bao gồm: Phương pháp số hóa địa hình đáy biển, phân chia lưới phi cấu trúc, thông số hóa 22 tham số động học trong module MIKE 21 Ecolab (hệ số phân hủy $\text{BOD}_5$, tốc độ nitrat hóa, hệ số lắng cặn), và các bước hiệu chuẩn theo chỉ số NASH. Quy trình này hoàn toàn có thể áp dụng cho các vùng vịnh bán kín khác tại Việt Nam như vịnh Cam Ranh, vịnh Nha Trang, vịnh Quy Nhơn hay vịnh Tiên Yên.

5. Lộ trình nghiên cứu 10 năm tiếp theo được vạch ra ra sao?

Lộ trình nghiên cứu hướng tới: (1) Xây dựng hệ thống Digital Twin (Bản sao số) thời gian thực cho toàn bộ lưu vực vịnh Cửa Lục – vịnh Hạ Long; (2) Mô hình hóa động lực học vận chuyển hạt vi nhựa và kim loại nặng dạng vết trong trầm tích biển; (3) Đánh giá tác động cộng hưởng của kịch bản biến đổi khí hậu - nước biển dâng và triều cường cực đoan đến sức chịu tải của vịnh.

Kết luận

Luận án tiến sĩ của tác giả Trần Đức Dũng đã giải quyết xuất sắc một bài toán khoa học và thực tiễn cấp bách về môi trường nước biển ven bờ thông qua 5 kết luận cốt lõi:

  1. Xác lập hệ phương pháp luận tích hợp: Tiên phong kết hợp thành công mô hình lưu vực SWAT và mô hình thủy động lực - sinh thái MIKE 21 Ecolab trong việc đánh giá tải lượng và chất lượng nước vịnh bán kín.
  2. Bức tranh hiện trạng định lượng: Làm sáng tỏ hiện trạng vượt tải nghiêm trọng của 4/6 chỉ tiêu đánh giá ($\text{BOD}_5$ vượt 114,41%, $\text{NH}_4^+$ vượt 528,65%, $\text{PO}_4^{3-}$ vượt 170,81%, $\text{Fe}$ vượt 241,14%) tại vịnh Cửa Lục năm 2018.
  3. Phát hiện quy luật động lực học theo mùa: Chứng minh cơ chế lan truyền ô nhiễm phụ thuộc chặt chẽ vào tương tác giữa dòng chảy lũ sông mùa mưa và áp lực triều dâng mùa khô, xác định 02 vùng ô nhiễm nghiêm trọng cục bộ phía Đông và ven Bãi Cháy.
  4. Dự báo sức tải đến năm 2030: Cảnh báo nguy cơ quá tải môi trường diện rộng đối với $\text{NH}_4^+$, $\text{PO}_4^{3-}$, $\text{Fe}$ và $\text{Mn}$ theo các kịch bản phát triển kinh tế - xã hội nếu thiếu các giải pháp công nghệ can thiệp.
  5. Hệ thống giải pháp khả thi: Đề xuất đồng bộ các giải pháp phi công trình (chính sách, hạn ngạch phát thải TMDL, hoàn thiện trạm quan trắc tự động) và giải pháp công trình mang tính đột phá kỹ thuật (công nghệ FRP-MBBR-100 xử lý nước thải sinh hoạt phân tán và hệ thống keo tụ - lắng lamen - lọc Nanofilter xử lý nước thải mỏ than), góp phần bảo vệ bền vững di sản thiên nhiên thế giới vịnh Hạ Long.