Nghiên Cứu Mô Hình Delft 3D Để Mô Phỏng Chất Lượng Nước Tại Đầm Phá Tam Giang Cầu Hai

Tài liệu nghiên cứu Luận án triển khai mô hình delft 3d mô phỏng chất lượng nước hệ đầm phá tam giang cầu hai và tính, tổng hợp lý thuyết và thực hành, cung cấp kiến thức chuyên

• Trường học: Trường Đại Học Khoa Học Tự Nhiên
• Chuyên ngành: Môi Trường
• Đề tài: Mô Hình Delft 3D: Phân Tích Chất Lượng Nước Đầm Phá Tam Giang Cầu Hai
• Loại tài liệu: luận án
• Năm xuất bản: 2016
• Địa điểm: Huế

Mục lục chi tiết

• LỜI MỞ ĐẦU
• 1. CHƯƠNG 1: TỔNG QUAN VẤN ĐỀ NGHIÊN CỨU VÀ KHU VỰC NGHIÊN CỨU
• 1.1. Tổng quan về sức tải môi trường
• 1.1.1. Các khái niệm

Tóm tắt

I. Giới thiệu về mô hình Delft 3D

Mô hình Delft 3D là một công cụ mạnh mẽ trong việc mô phỏng các quá trình thủy động lực và chất lượng nước trong các hệ thống thủy vực. Mô hình này cho phép phân tích và dự đoán sự lan truyền của các chất ô nhiễm trong môi trường nước, từ đó hỗ trợ cho việc quản lý và bảo vệ chất lượng nước. Đặc biệt, trong nghiên cứu về đầm phá Tam Giang - Cầu Hai, mô hình này được áp dụng để đánh giá tác động của các nguồn thải đến chất lượng nước. Việc sử dụng mô hình Delft 3D giúp xác định các kịch bản phát triển và ảnh hưởng của chúng đến sức tải môi trường của hệ đầm phá. Theo nghiên cứu, mô hình này đã cho thấy khả năng mô phỏng chính xác các yếu tố như dòng chảy, nồng độ chất ô nhiễm và sự phân bố của chúng trong không gian và thời gian.

1.1. Nguyên lý hoạt động của mô hình

Mô hình Delft 3D hoạt động dựa trên các phương trình thủy động lực học và chất lượng nước. Nó sử dụng các dữ liệu đầu vào như lưu lượng nước, nồng độ chất ô nhiễm và các yếu tố môi trường khác để mô phỏng sự lan truyền và biến đổi của các chất trong nước. Mô hình này có khả năng xử lý các kịch bản khác nhau, từ đó đưa ra các dự báo về chất lượng nước trong tương lai. Việc áp dụng mô hình này trong nghiên cứu đầm phá Tam Giang - Cầu Hai đã cho thấy sự cần thiết trong việc quản lý nguồn thải và bảo vệ môi trường nước, đặc biệt là trong bối cảnh phát triển kinh tế - xã hội tại khu vực này.

II. Phân tích chất lượng nước tại đầm phá Tam Giang Cầu Hai

Chất lượng nước tại đầm phá Tam Giang - Cầu Hai đang chịu ảnh hưởng lớn từ các hoạt động kinh tế - xã hội. Các nguồn thải từ nông nghiệp, nuôi trồng thủy sản và sinh hoạt đã làm gia tăng nồng độ các chất ô nhiễm như BOD5, COD, amoni, phosphat và nitrat. Nghiên cứu cho thấy rằng, nếu không có biện pháp quản lý hiệu quả, chất lượng nước sẽ tiếp tục suy giảm, ảnh hưởng đến hệ sinh thái và sức khỏe cộng đồng. Việc đánh giá chất lượng nước thông qua mô hình Delft 3D đã chỉ ra rằng, nồng độ các chất ô nhiễm vượt quá ngưỡng cho phép, gây ra những tác động tiêu cực đến sinh vật thủy sinh. Điều này nhấn mạnh tầm quan trọng của việc áp dụng các biện pháp quản lý bền vững nhằm bảo vệ chất lượng nước và hệ sinh thái tại khu vực này.

2.1. Các yếu tố ảnh hưởng đến chất lượng nước

Các yếu tố ảnh hưởng đến chất lượng nước tại đầm phá Tam Giang - Cầu Hai bao gồm tải lượng chất thải từ các hoạt động nông nghiệp, nuôi trồng thủy sản và sinh hoạt. Nghiên cứu đã chỉ ra rằng, sự gia tăng nồng độ các chất ô nhiễm như BOD5 và COD là do sự gia tăng hoạt động kinh tế trong khu vực. Các chất dinh dưỡng như nitrat và phosphat cũng gia tăng do việc sử dụng phân bón trong nông nghiệp. Điều này dẫn đến hiện tượng phú dưỡng, làm giảm chất lượng nước và ảnh hưởng đến sự phát triển của các loài thủy sinh. Việc áp dụng mô hình Delft 3D giúp xác định rõ ràng các nguồn ô nhiễm và mức độ ảnh hưởng của chúng đến chất lượng nước, từ đó đưa ra các giải pháp quản lý hiệu quả.

III. Đánh giá sức tải môi trường của đầm phá

Đánh giá sức tải môi trường là một phần quan trọng trong nghiên cứu về đầm phá Tam Giang - Cầu Hai. Sức tải môi trường được xác định dựa trên khả năng tự làm sạch của hệ thống nước và mức độ ô nhiễm hiện tại. Nghiên cứu cho thấy rằng, sức tải môi trường của đầm phá đang bị đe dọa do sự gia tăng chất thải từ các hoạt động kinh tế. Việc xác định ngưỡng sức tải tối đa cho phép là cần thiết để bảo vệ chất lượng nước và duy trì sự cân bằng sinh thái. Mô hình Delft 3D đã cung cấp các kịch bản dự báo về sức tải môi trường, giúp các nhà quản lý đưa ra các quyết định hợp lý trong việc quản lý nguồn thải và bảo vệ môi trường.

3.1. Phương pháp đánh giá sức tải môi trường

Phương pháp đánh giá sức tải môi trường trong nghiên cứu này bao gồm việc sử dụng mô hình Delft 3D để mô phỏng sự lan truyền của các chất ô nhiễm và khả năng tự làm sạch của hệ thống. Các yếu tố như BOD5, COD, nitrat và phosphat được xem xét để xác định sức tải tối đa cho phép. Nghiên cứu đã chỉ ra rằng, sức tải môi trường của đầm phá Tam Giang - Cầu Hai đang ở mức báo động, cần có các biện pháp can thiệp kịp thời để bảo vệ chất lượng nước. Việc áp dụng mô hình này không chỉ giúp đánh giá chính xác sức tải môi trường mà còn cung cấp cơ sở dữ liệu quan trọng cho việc quản lý và phát triển bền vững khu vực.

Nội dung chính

Tổng quan về luận án

Hệ đầm phá Tam Giang - Cầu Hai (tỉnh Thừa Thiên Huế) là hệ thống đầm phá ven bờ lớn nhất khu vực Đông Nam Á với chiều dài 68 km, diện tích mặt nước 216 km², là nơi cư ngụ của gần 1.000 loài động vật và thực vật thủy sinh có giá trị kinh tế và bảo tồn sinh học cao [1]. Nằm tại điểm tiếp nhận nguồn thải hạ lưu từ lưu vực các sông chính như sông Ô Lâu, sông Hương (nhánh Tả Trạch, Hữu Trạch), sông Bồ, sông Truồi, sông Cầu Hai và sông Nong, đầm phá chịu áp lực tích tụ chất ô nhiễm nghiêm trọng từ các hoạt động nhân sinh ven bờ bao gồm đô thị hóa, du lịch, trồng trọt, chăn nuôi và nuôi trồng thủy sản (NTTS). Thực trạng phát triển kinh tế nhanh chóng đã đẩy tải lượng hữu cơ (BOD, COD) và muối dinh dưỡng (N, P) gia tăng liên tục, đe dọa trực tiếp đến cấu trúc sinh thái và chức năng tự làm sạch tự nhiên của thủy vực.

Khoảng trống nghiên cứu (Research Gap) then chốt trong các công trình khoa học trước đây tại Việt Nam là việc đánh giá sức tải môi trường (STMT) chủ yếu chỉ dựa vào các tiêu chuẩn hành chính tĩnh (như QCVN) hoặc các mô hình sinh thái đơn lẻ cho một vùng nuôi cụ thể (ví dụ mô hình Ecopath tại đầm Sam Chuồn [98] hay phương pháp đồng vị đánh dấu của Võ Duy Sơn và Nguyễn Tác An tại đầm Nha Phu [45]). Các tiếp cận này bỏ qua động lực học trao đổi nước phức tạp qua hai cửa biển Thuận An, Tư Hiền, cơ chế phân tầng dòng chảy hai chiều vào mùa khô, và đặc biệt là chưa xây dựng được ngưỡng chịu tải sinh thái thực sự dựa trên độc học sinh vật thủy sinh và năng lực đồng hóa tự nhiên của thủy vực.

Nhằm giải quyết triệt để khoảng trống này, luận án tập trung giải quyết 3 câu hỏi nghiên cứu và kiểm chứng 3 giả thuyết khoa học:

• Câu hỏi nghiên cứu 1 (RQ1): Tổng tải lượng các chất ô nhiễm hữu cơ và dinh dưỡng (C, N, P) phát sinh từ các nguồn thải lục địa và nội tại đổ vào đầm phá Tam Giang - Cầu Hai là bao nhiêu theo hiện trạng và kịch bản phát triển 2020, 2030?
• Câu hỏi nghiên cứu 2 (RQ2): Quy luật thủy động lực và lan truyền chất ô nhiễm trong hệ đầm phá biến đổi theo không gian và thời gian (mùa khô và mùa mưa) như thế nào dưới tác động của dòng chảy sông và thủy triều?
• Câu hỏi nghiên cứu 3 (RQ3): Giới hạn sức chịu tải tối đa của hệ đầm phá đối với các yếu tố C, N, P là bao nhiêu khi xét theo cả ngưỡng quy chuẩn quốc gia (QCVN 10-MT:2015/BTNMT) và ngưỡng an toàn sinh thái bảo vệ sinh vật thủy sinh?
• Giả thuyết khoa học 1 (H1): Khả năng tiếp nhận và tự làm sạch chất ô nhiễm của đầm phá có sự bất đối xứng sâu sắc giữa mùa mưa và mùa khô do chế độ thủy văn chi phối (lưu lượng sông mùa mưa chiếm 70 - 75% tổng lượng cả năm).
• Giả thuyết khoa học 2 (H2): Các khu vực bán khép kín như đầm Sam Chuồn và đầm Thủy Tú có thời gian lưu nước ($\tau$) lớn, dẫn đến nguy cơ quá tải cục bộ các chất dinh dưỡng (N-NH4+, P-PO43-) ngay cả khi tải lượng tổng thể toàn đầm phá chưa chạm ngưỡng.
• Giả thuyết khoa học 3 (H3): Ngưỡng sức tải xác định dựa trên giới hạn sinh thái (DO $\ge$ 5 mg/L, NH3 < 0,05 mg/L, P-PO43- < 0,045 mg/L) chặt chẽ và phản ánh chính xác rủi ro bùng phát phú dưỡng hơn so với các ngưỡng nồng độ giới hạn cho phép trong quy chuẩn quản lý hiện hành.

Khung lý thuyết của luận án được xây dựng dựa trên sự tích hợp giữa Lý thuyết Cân bằng Vật chất Thủy vực (Vlađimirov, 1982), Nguyên lý Sức tải Môi trường Ven bờ (GESAMP, 1986) và Động học Sinh địa hóa Chuyển hóa Dinh dưỡng (Delft3D-WAQ). Đóng góp đột phá của công trình là lượng hóa được sức tải môi trường chi tiết theo không gian cho 5 thông số cốt lõi (BOD5, COD, N-NH4+, N-NO3-, P-PO43-), cung cấp công cụ mô phỏng số phân giải cao hỗ trợ điều phối không gian phát triển kinh tế xã hội tỉnh Thừa Thiên Huế với độ tin cậy được kiểm chứng qua chuỗi dữ liệu thực địa 24 giờ liên tục qua các mùa.

Literature Review và Positioning

Nghiên cứu về sức chịu tải môi trường nước (Water Environmental Carrying Capacity) đã trải qua ba thập kỷ tiến hóa với ba trường phái học thuật chủ đạo:

Trường phái thứ nhất tập trung vào mô hình sinh địa hóa phục vụ NTTS và sức tải sinh thái đơn loài. Khởi xướng bởi Duarte và cs (2003) [6], Brigolin và cs (2008) [16], Byron và cs (2011) [18], hướng tiếp cận này gắn liền sức tải với mật độ nuôi thả tối đa của các loài thân mềm hai mảnh vỏ hoặc cá mà không gây suy thoái nền đáy. Tại vịnh Goseung (Hàn Quốc), mô hình thủy động lực - phú dưỡng 3 chiều HEM-3D kết hợp mô hình tăng trưởng ngao của Nakata (1983) đã xác định mật độ nuôi tối ưu là 16 cá thể/m³, tương ứng sức tải 1.500 tấn thịt [20, 21]. Tương tự, tại vịnh Hiroshima (Nhật Bản), mô hình hệ sinh thái chỉ ra rằng để nuôi ngao bền vững và duy trì nồng độ Chlorophyll-a ở mức 7 $\mu$g/L, tải lượng phospho tổng số từ sông Ohta phải được kiểm soát nghiêm ngặt ở mức 0,5 tấn/ngày [22]. Tuy nhiên, trường phái này bị hạn chế do chỉ tối ưu hóa quy mô trang trại đơn lẻ, bỏ qua các nguồn thải tổng hợp từ lưu vực lục địa.

Trường phái thứ hai tiếp cận theo hướng Tổng tải lượng xả thải tối đa hàng ngày (TMDL - Total Maximum Daily Load) và Năng lực đồng hóa của thủy vực. Khởi đầu từ định nghĩa nền tảng của GESAMP (1986) [3] với công thức: $$EC_{max} = \frac{(C_{tiêu chuẩn} - C_{hiện tại}) \cdot V}{\tau}$$ trong đó $\tau$ là thời gian lưu nước và $V$ là thể tích thủy vực. Hướng tiếp cận này được phát triển vượt bậc nhờ mô hình số 2D/3D mô phỏng lan truyền chất ô nhiễm. Điển hình là Chang Hee-Lee (2003) ứng dụng hệ thống quản lý TPLMS tại hồ Shihwa (Hàn Quốc) nhằm duy trì nồng độ COD mục tiêu ở mức 2 ppm [12]; Ayeon Lee và cs (2013) sử dụng mô hình HSPF trên lưu vực sông Anyangcheon giúp cắt giảm BOD từ 11,21 mg/L xuống 7,39 mg/L [29]. Tại Trung Quốc, Li KeQiang và cs (2012, 2014) tích hợp mô hình thủy động lực 3D tính toán năng lực môi trường đối với kim loại nặng tại vịnh Jiaozhou và vịnh Jinzhou [25, 26]; Zhao và cs (2012, 2015) thiết lập mô hình cho hồ Fuxian và vịnh Quanzhou, chứng minh nguồn thải sinh hoạt chiếm 45,2 - 46,5% và đề xuất cắt giảm 57 - 68% tải lượng TN, TP, COD [27, 28].

Trường phái thứ ba sử dụng Hệ thống chỉ số tổng hợp và Quy trình phân tích thứ bậc (AHP). Đại diện bởi Yongquan Yin (2010) [38], Meriem Naimi Ait-Aoudia và cs (2016) [14] tại Algiers, và Shuping Wang (2017) [36] tại Thanh Đảo. Mặc dù bao quát được các chiều cạnh kinh tế - xã hội, cách tiếp cận này phụ thuộc nhiều vào đánh giá chuyên gia chủ quan và không phản ánh được cơ chế tương tác thủy lực thời gian thực.

Tiêu chí so sánh	Nghiên cứu quốc tế điển hình (Hồ Shihwa - Chang Hee-Lee, 2003; Vịnh Quanzhou - Zhao et al., 2015)	Nghiên cứu tại Đầm phá Tam Giang - Cầu Hai (Luận án)
Đối tượng thủy vực	Hồ nhân tạo/vịnh mở ven biển, chế độ trao đổi nước tương đối đồng nhất.	Đầm phá ven bờ dạng đới hẹp dài (68 km), thông biển qua 2 cửa (Thuận An, Tư Hiền), có cồn cát chắn bọc ngoài.
Cơ chế thủy văn	Dao động dòng chảy nội lưu ở mức biến thiên trung bình.	Độ bất đối xứng cực đoan: Mùa mưa tiếp nhận 70 - 75% lưu lượng nước ngọt từ 6 hệ thống sông lớn; mùa khô phân tầng 2 chiều sâu sắc.
Công cụ mô hình hóa	WASP5, HSPF, mô hình hộp sinh địa hóa 1D/2D hoặc mô hình hộp 3D.	Tích hợp bộ mô hình số trị Delft3D (Delft3D-Flow kết hợp Delft3D-WAQ) trên lưới cong trực giao bám sát địa hình đáy.
Căn cứ xác lập ngưỡng	Chủ yếu áp dụng tiêu chuẩn chất lượng nước hành chính quốc gia (CNWQS loại I/II hoặc US EPA TMDL).	Kết hợp đồng thời 3 ngưỡng: Ngưỡng quy chuẩn QCVN, Ngưỡng sức tải tối đa, và Ngưỡng độc học sinh thái (DO $\ge$ 5 mg/L, NH3 < 0,05 mg/L, PO43- < 0,045 mg/L).

Đứng trên bình diện học thuật, luận án đã định vị chính xác điểm giao thoa giữa cơ chế thủy động lực học ven bờ và động học chuyển hóa sinh địa hóa, đưa lý thuyết sức tải môi trường vượt ra khỏi giới hạn của các tiêu chuẩn tĩnh để tiến tới mô hình hóa khả năng tự làm sạch thực chất của một hệ sinh thái đầm phá nhiệt đới gió mùa phức tạp.

Đó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 và làm sâu sắc thêm Lý thuyết Sức tải Môi trường của GESAMP (1986) và Furuya (2003) bằng cách tái định nghĩa STMT trong mối tương quan hữu cơ với động thái tự làm sạch tự nhiên (Natural Assimilative Capacity). Luận án khẳng định: "Sức tải môi trường là lượng chất tối đa mà thủy vực có thể chấp nhận sao cho vẫn duy trì được sự sinh trưởng và phát triển của các loài sinh vật". Luận điểm này bác bỏ quan điểm đơn thuần coi sức tải là khoảng trống nồng độ hành chính $(C_{tiêu chuẩn} - C_{hiện tại})$, vốn dễ bị bóp méo bởi các nhóm lợi ích kinh tế khi xây dựng quy chuẩn kỹ thuật.

Khung lý thuyết của đề tài hoàn thiện phương trình cân bằng chất của Vlađimirov (1982) [5]: $$\Delta C = C_{vào} - C_{ra}$$ Trong điều kiện bảo tồn trạng thái bình thường của thủy vực ($0 < C \le NĐTHCP$ và $\Delta C < 0$), luận án đã cụ thể hóa các thành phần của $C_{ra}$ thông qua 3 cơ chế tự làm sạch:

1. Cơ chế vật lý: Pha loãng, xáo trộn thủy triều, khuếch tán dòng chảy và sa lắng bồi tích xuống tầng đáy trầm tích.
2. Cơ chế hóa học: Oxy hóa quang hóa UV, phản ứng oxy hóa khử và trung hòa ion.
3. Cơ chế sinh học ("Tự làm sạch đúng nghĩa"): Khoáng hóa hữu cơ hiếu khí/kỵ khí bởi vi sinh vật, chu trình chuyển hóa nitrat hóa và phản nitrat hóa.

Hệ thống giả thuyết lý thuyết được lượng hóa thành 3 mệnh đề cốt lõi:

• Mệnh đề 1 (Proposition 1): Năng lực đồng hóa hữu cơ tỷ lệ thuận với tốc độ tái cấp khí oxy hòa tan bề mặt và thể tích khối nước trao đổi qua triều, nhưng bị giới hạn bởi nồng độ oxy hòa tan tối thiểu bảo đảm hô hấp sinh vật ($DO \ge 5\text{ mg/L}$ theo Stickney, 2005 [121]).
• Mệnh đề 2 (Proposition 2): Độc tính amoniac trong đầm phá là một hàm phi tuyến phụ thuộc vào trạng thái cân bằng nhiệt dịch ($NH_4^+ \leftrightarrow NH_3 + H^+$). Tại ngưỡng nhiệt độ cao (33°C) và kiềm hóa (pH đạt 9,0), nồng độ amoni tổng số $N\text{-}NH_4^+ = 0,1\text{ mg/L}$ sẽ chuyển hóa giải phóng khí $NH_3$ tự do vượt ngưỡng $0,05\text{ mg/L}$ (Boardman và nnk, 2004 [122]), gây ức chế thần kinh và tỷ lệ tử vong cao ở ấu trùng thủy sản.
• Mệnh đề 3 (Proposition 3): Phospho đóng vai trò là yếu tố hạn chế bậc nhất (limiting nutrient); nồng độ $P\text{-}PO_4^{3-}$ vượt quá $0,045\text{ mg/L}$ (tiêu chuẩn ASEAN [125] và SteffIii, 2003 [124]) là điều kiện động học kích hoạt quá trình siêu phú dưỡng và hiện tượng tảo nở hoa độc hại (Harmful Algal Blooms).

flowchart TD
    subgraph Drivers ["Nguồn thải Lục địa & Nội sinh"]
        D1["Nước thải Dân cư & Du lịch (WHO 1993)"]
        D2["Nước thải Chăn nuôi & Trồng trọt"]
        D3["Nước thải Nuôi trồng Thủy sản"]
        D4["Nước thải Công nghiệp & Rửa trôi đất"]
    end

    subgraph Hydro ["Thủy động lực học Delft3D-Flow"]
        H1["Dòng chảy sông Ô Lâu, Hương, Bồ, Truồi, Cầu Hai"]
        H2["Dao động Thủy triều Vịnh Bắc Bộ qua Cửa Thuận An & Tư Hiền"]
        H3["Trường Vận tốc, Phân tầng Muối & Thời gian lưu nước (τ)"]
    end

    subgraph WQ ["Động học Sinh địa hóa Delft3D-WAQ"]
        W1["Quá trình Vật lý: Khuếch tán, Pha loãng, Lắng đọng"]
        W2["Quá trình Hóa học: Oxy hóa khử, Cân bằng Axit-Bazơ"]
        W3["Quá trình Sinh học: Nitrat hóa, Phản nitrat hóa, Khoáng hóa C-N-P"]
    end

    subgraph Thresholds ["Ngưỡng Đánh giá Sức tải"]
        T1["Ngưỡng QCVN 10-MT:2015/BTNMT"]
        T2["Ngưỡng Tối đa Hành chính"]
        T3["Ngưỡng Độc học Sinh thái: DO ≥ 5 mg/L, NH3 < 0.05 mg/L, PO43- < 0.045 mg/L"]
    end

    Drivers --> Hydro
    Hydro --> WQ
    WQ --> Thresholds
    Thresholds --> Output["Sức tải Môi trường Hữu dụng & Bản đồ Phân vùng Quản trị"]

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

Khung phân tích của luận án tích hợp đa tầng giữa:

1. Lý thuyết Lan truyền Bình lưu - Phân tán (Advection-Dispersion Theory) mô tả quá trình vận chuyển khối chất lỏng.
2. Lý thuyết Động học Sinh địa hóa vùng ven bờ (Coastal Biogeochemical Kinetics) mô phỏng sự tương tác giữa chu trình carbon, nitơ, phospho và oxy hòa tan.
3. Lý thuyết Quản lý Lưu vực Tích hợp (Integrated Watershed Management) kết nối hoạt động phát triển kinh tế xã hội vùng nội lục với sức khỏe vùng nước tiếp nhận ven biển.

Phương pháp tiếp cận giải quyết bài toán biên phức tạp của đầm phá nhiệt đới: Chiều dài lớn (68 km), lòng đầm nông (trung bình 1,5 - 2,5 m), chế độ triều yếu (nhật triều không đều với biên độ nhỏ) nhưng chịu lũ đột biến từ vùng núi Bạch Mã với vũ lượng trên 3.000 mm/năm. Luận án đã thiết lập điều kiện biên ràng buộc chặt chẽ, cho phép tách biệt rõ ràng giữa khả năng tự làm sạch biểu kiến và tự làm sạch thực chất.

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ủ thế giới quan thực chứng (Positivism) kết hợp phương pháp luận mô hình hóa số trị (Numerical Modeling) và kiểm chứng thực nghiệm hiện trường nghiêm ngặt. Hệ thống nghiên cứu được thiết kế đa cấp độ (Multi-level design):

• Cấp độ vi mô: Động học phân hủy sinh hóa các hợp chất C, N, P trong cột nước và lớp ranh giới bùn - nước.
• Cấp độ trung mô: Thủy động lực học các phân vùng đầm phá đặc thù (phá Tam Giang, đầm Thủy Tú, đầm Cầu Hai, đầm Sam Chuồn).
• Cấp độ vĩ mô: Toàn bộ lưu vực đầm phá Tam Giang - Cầu Hai diện tích mặt nước 216 km² cùng dải đệm đất liền thuộc 5 huyện/thị xã ven đầm phá (Phong Điền, Quảng Điền, Hương Trà, Phú Vang, Phú Lộc).

Quy trình nghiên cứu rigorous

Quy trình thu thập và phân tích dữ liệu tuân thủ nghiêm ngặt các quy chuẩn kỹ thuật quốc gia và chuẩn mực quốc tế:

• Khung pháp lý kỹ thuật: Thông tư số 34/2010/TT-BTNMT quy định kỹ thuật điều tra, khảo sát hải văn, hóa học và môi trường vùng ven bờ và hải đảo; Thông tư 76/2017/TT-BTNMT về đánh giá khả năng tiếp nhận nước thải và sức chịu tải nguồn nước mặt.
• Chiến dịch khảo sát thực địa: Tiến hành 2 đợt đo đạc quy mô lớn đại diện cho 2 mùa tương phản:
  • Mùa mưa: Đo liên tục 24 giờ từ 17h00 ngày 25/11/2011 đến 17h00 ngày 26/11/2011.
  • Mùa khô: Đo liên tục 24 giờ từ 16h00 ngày 19/5/2012 đến 16h00 ngày 20/5/2012.
  • Tần suất quan trắc: 4 giờ/ốp (tổng cộng 6 ốp quan trắc hoàn chỉnh trong 24 giờ) tại 3 trạm thủy văn - chất lượng nước then chốt: Trạm phá Tam Giang (đại diện vùng phân lưu phía Bắc), Trạm đầm Sam (đại diện vùng nuôi trồng tập trung và tiếp nhận nước sông Hương), Trạm đầm Cầu Hai (đại diện vùng nước lợ phía Nam giao thoa cửa Tư Hiền).
• Quy chuẩn thu và bảo quản mẫu: Thu tầng mặt (cách mặt nước 0,5 m) và tầng đáy (cách đáy 0,5 m) bằng thiết bị lấy mẫu chuyên dụng Bathomet theo TCVN 5998:1995 (ISO 5667-9:1992). Bảo quản mẫu lạnh tức thì theo TCVN 5993:1995 (ISO 5667-3:1985).

+---------------------------------------------------------------------------------------------------+
|                                  QUY TRÌNH THỰC NGHIỆM VÀ MÔ HÌNH HÓA                            |
+---------------------------------------------------------------------------------------------------+
|  [ĐIỀU TRA NGUỒN THẢI VEN BỜ]       [KHẢO SÁT HIỆN TRƯỜNG 24H LIÊN TỤC]    [ĐO ĐẠC KHÍ TƯỢNG - THỦY VĂN]  |
|  - Dân cư (WHO: 26.24 kg COD/ng)    - 3 Trạm liên tục: Tam Giang, Sam,      - 6 Lưu vực sông đổ vào      |
|  - Chăn nuôi (Trâu bò, Lợn, Cầm)      Cầu Hai (4h/ốp, tầng mặt & đáy)       - Trạm khí tượng Bạch Mã     |
|  - NTTS (Tôm sú, Cá lồng)           - Thiết bị lấy mẫu Bathomet              - Dữ liệu triều Thuận An,    |
|  - Công nghiệp & Du lịch            - Phân tích HACH DR/2800, Winkler         Tư Hiền                     |
+----------------------------------+----------------------------------+-----------------------------+
                                   |                                  |
                                   v                                  v
+---------------------------------------------------------------------------------------------------+
|                                   THIẾT LẬP MÔ HÌNH SỐ TRỊ DELFT-3D                               |
|  +-------------------------------------------------+--------------------------------------------+ |
|  | Module Thủy động lực (Delft3D-Flow)             | Module Chất lượng nước (Delft3D-WAQ)       | |
|  | - Lưới cong trực giao bám sát đường bờ          | - 21 biến số chất lượng nước               | |
|  | - Hiệu chuẩn mực nước, lưu lượng, độ mặn        | - Động học phân hủy C, N, P, DO            | |
|  +-------------------------------------------------+--------------------------------------------+ |
+---------------------------------------------------------------------------------------------------+
                                   |
                                   v
+---------------------------------------------------------------------------------------------------+
|                             TÍNH TOÁN SỨC TẢI VÀ MÔ PHỎNG KỊCH BẢN                                |
|  - Kịch bản hiện trạng (2011 - 2012) & Kiểm chứng 2016                                             |
|  - Kịch bản phát triển kinh tế - xã hội 2020 & 2030                                               |
|  - Đánh giá theo 3 ngưỡng: QCVN 10-MT, Ngưỡng tối đa, Ngưỡng độc thái sinh học                   |
+---------------------------------------------------------------------------------------------------+

Data và phân tích

Thiết bị và phương pháp phân tích phòng thí nghiệm:

• Nhiệt độ nước và không khí: Nhiệt kế thủy ngân chuyên dụng (độ chính xác $\pm 0,1^\circ\text{C}$).
• Độ muối ($S‰$): Khúc xạ kế cầm tay (Hand Refractometer, độ chính xác $\pm 1‰$).
• pH: Máy đo pH điện cực thủy tinh chuẩn hóa định kỳ (độ chính xác $\pm 0,01\text{ đơn vị}$).
• Oxy hòa tan (DO): Đo máy đo oxy quang học kết hợp chuẩn độ Winkler đối chứng (độ chính xác $\pm 0,01\text{ mg/L}$).
• Nhu cầu oxy sinh hóa ($BOD_5$): Phương pháp ủ tối trực tiếp ở $20^\circ\text{C}$ trong 5 ngày.
• Nhu cầu oxy hóa học ($COD$): Oxy hóa bằng Kali Permanganat ($KMnO_4$) trong môi trường kiềm nóng.
• Muối dinh dưỡng vô cơ hòa tan ($P\text{-}PO_4^{3-}, N\text{-}NO_2^-, N\text{-}NH_4^+ + NH_3$): Phân tích so màu trên máy quang phổ kế tử ngoại khả kiến DR/2800 HACH (USA).

Hệ thống tính toán tải lượng ô nhiễm nguồn phát sinh:

• Nước thải sinh hoạt dân cư và du lịch: Áp dụng định mức phát thải WHO (1993) [126, 127, 128]: $COD = 26,24\text{ kg/người/năm}$, $BOD_5 = 16,4\text{ kg/người/năm}$, $N\text{-Tổng} = 2,2\text{ kg/người/năm}$, $P\text{-Tổng} = 0,4\text{ kg/người/năm}$, $N\text{-NH}_4^+ = 1,20\text{ kg/người/năm}$, $P\text{-PO}_4^{3-} = 0,12\text{ kg/người/năm}$.
• Nước thải chăn nuôi: Tính trên quy mô tổng đàn gia súc, gia cầm với hệ số WHO (1993): Trâu bò phát sinh $262,4\text{ kg COD/con/năm}$ và $164\text{ kg BOD/con/năm}$; Lợn phát sinh $52,64\text{ kg COD/con/năm}$ và $32,9\text{ kg BOD/con/năm}$; Gia cầm phát sinh $2,57\text{ kg COD/con/năm}$.
• Nước thải nuôi trồng thủy sản: Dựa trên nghiên cứu của San Diego - McGlone (2000) và tài liệu chuyên ngành [130, 133, 134]: Nuôi thâm canh tôm sú phát thải $28,4\text{ kg COD/tấn sản phẩm/năm}$, $8,1\text{ kg BOD/tấn/năm}$, $5,2\text{ kg N-T/tấn/năm}$, $4,7\text{ kg P-T/tấn/năm}$; Nuôi cá lồng phát thải $15,9\text{ kg COD/tấn/năm}$, $4,5\text{ kg BOD/tấn/năm}$, $2,9\text{ kg N-T/tấn/năm}$, $2,6\text{ kg P-T/tấn/năm}$.
• Nước thải công nghiệp: Tính toán theo hệ số phát thải sản phẩm ngành (bia $150\text{ mg COD/L}$, bột mì $1.500\text{ mg COD/L}$, chế biến hải sản đông lạnh $1.950\text{ mg COD/L}$ với định mức $95\text{ m}^3\text{ nước thải/tấn sản phẩm}$) và ngành dệt sợi tổng hợp [132].

Mô hình toán số trị: Ứng dụng gói phần mềm bản quyền Delft3D (Viện Thủy lực Delft, Hà Lan) gồm module thủy động lực Delft3D-Flow thiết lập trên lưới cong trực giao và module chất lượng nước Delft3D-WAQ mô phỏng 21 biến trạng thái. Mô hình được hiệu chỉnh và thẩm định thông qua sai số tuyệt đối trung bình (MAE) và hệ số hiệu quả Nash-Sutcliffe (NSE) đạt chuẩn tin cậy cao trước khi thực hiện các kịch bản mô phỏng tương lai.

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

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

Thứ nhất, quy luật bất đối xứng thủy văn và phân hóa không gian của sức tải môi trường: Chế độ dòng chảy sông đổ vào đầm phá có sự phân hóa cực đoan. Lưu lượng mùa mưa (tháng 9 - 12) chiếm tới 70 - 75% tổng lượng dòng chảy năm (sông Tả Trạch đạt 1,58 tỷ m³/năm, sông Hữu Trạch 1,18 tỷ m³/năm, sông Bồ 1,42 tỷ m³/năm, sông Ô Lâu 0,545 tỷ m³/năm). Vào mùa lũ, thời gian lưu nước ($\tau$) tại phá Tam Giang giảm xuống mức tối thiểu (dưới 3 - 5 ngày), năng lực rửa trôi và đẩy chất ô nhiễm ra biển qua cửa Thuận An đạt đỉnh. Ngược lại, vào mùa khô (tháng 5), lưu lượng sông suy kiệt sâu sắc (sông Ô Lâu chỉ đạt 7,27 m³/s), hiện tượng nêm mặn lấn sâu gây phân tầng mật độ hai chiều tại cửa sông Hương và khu vực đầm Thủy Tú, kéo dài thời gian lưu nước lên trên 25 - 30 ngày, khiến sức tải tiếp nhận ô nhiễm hữu cơ và dinh dưỡng sụt giảm từ 60% đến 80% so với mùa mưa.

Thứ hai, phát hiện sự gia tăng đáng báo động của muối dinh dưỡng và áp lực phú dưỡng: Kết quả khảo sát và tính toán chỉ ra hàm lượng Nitrat ($N\text{-}NO_3^-$) trong nước đầm phá dao động từ 88,5 $\mu$g/L đến 204,2 $\mu$g/L, vượt giới hạn tiêu chuẩn chất lượng nước biển ven bờ ASEAN từ 1,4 đến 3,4 lần. Hàm lượng Silicat ($Si_2O_3^{2-}$) dao động từ 718 đến 4.400 $\mu$g/L, tăng từ 1,0 - 2,1 lần trong mùa khô và từ 2,2 - 3,4 lần trong mùa mưa giai đoạn 1993 - 2009. Hàm lượng Phosphat ($P\text{-}PO_4^{3-}$) dao động từ 3,6 đến 44,1 $\mu$g/L, tiệm cận sát ngưỡng bùng phát tảo nở hoa ($45\text{ }\mu\text{g/L} = 0,045\text{ mg/L}$).

Thứ ba, sự xung đột giữa chuẩn nồng độ hành chính và độc học amoniac phi ion hóa: Mặc dù hàm lượng amoni tổng số ($N\text{-}NH_4^+ + NH_3$) trong nước dao động từ 17,1 đến 112,5 $\mu$g/L (hầu hết nằm trong giới hạn cho phép theo QCVN 10-MT:2015/BTNMT), các phân tích độc học nhiệt động chỉ ra rằng vào các đợt nắng nóng mùa khô (nhiệt độ nước chạm $33^\circ\text{C}$, pH kiềm hóa đạt 8,5 - 9,0), tỷ lệ phân ly tạo khí $NH_3$ tự do vượt ngưỡng độc tính $0,05\text{ mg/L}$ tại các vũng nuôi tôm thâm canh đầm Sam Chuồn. Đây là phát hiện phản trực giác giải thích nguyên nhân các đợt tôm cá chết hàng loạt vào mùa hè dù các chỉ tiêu quan trắc đơn lẻ đều đạt quy chuẩn QCVN.

Thứ tư, dự báo quá tải sức tải theo các kịch bản 2020 và 2030: Mô phỏng các kịch bản phát triển kinh tế xã hội cho thấy nếu không có các biện pháp can thiệp thu gom và xử lý nước thải tập trung, đầm Thủy Tú và đầm Sam Chuồn sẽ rơi vào tình trạng vượt ngưỡng sức tải hữu dụng đối với $COD$ và $P\text{-}PO_4^{3-}$ lần lượt là 135% và 160% vào năm 2020, và tăng lên trên 210% vào năm 2030, dẫn tới nguy cơ sụp đổ cục bộ dịch vụ hệ sinh thái đầm phá.

graph LR
    subgraph Spatial_Dynamics ["Phân hóa Thủy văn & Không gian"]
        A["Mùa Mưa (T9-T12)<br>Dòng chảy sông chiếm 70-75%<br>Thời gian lưu τ < 5 ngày"] -->|Rửa trôi mạnh| B["Sức tải hữu dụng TĂNG CAO<br>Ô nhiễm trôi nhanh ra biển"]
        C["Mùa Khô (T5)<br>Dòng chảy sông kiệt quệ<br>Thời gian lưu τ > 25-30 ngày"] -->|Nêm mặn & Phân tầng 2 chiều| D["Sức tải hữu dụng GIẢM 60-80%<br>Tích tụ hữu cơ & Dinh dưỡng"]
    end

    subgraph Ecological_Risks ["Điểm nghẽn Sinh thái Cục bộ"]
        E["Nhiệt độ cao (33°C) + pH kiềm (8.5-9.0)"] -->|Dịch chuyển cân bằng hóa học| F["Khí độc NH3 > 0.05 mg/L<br>(Độc học thủy sản)"]
        G["P-PO43- tiệm cận 0.045 mg/L"] -->|Phospho là yếu tố hạn chế| H["Nguy cơ Tảo nở hoa & Phú dưỡng"]
    end

Implications đa chiều

• Về mặt lý thuyết: Đặt nền móng cho việc chuẩn hóa phương pháp luận tính toán sức tải môi trường cho các đầm phá ven biển tại các quốc gia nhiệt đới gió mùa, chuyển đổi tư duy quản lý từ "kiểm soát nồng độ cuối đường ống" sang "kiểm soát tổng tải lượng dựa trên năng lực đồng hóa sinh thái".
• Về mặt phương pháp luận: Cung cấp quy trình tích hợp giữa điều tra phát thải diện rộng (WHO emission factors) với hệ công cụ mô hình số trị Delft3D có khả năng chuyển giao và áp dụng trực tiếp cho 11 hệ đầm phá còn lại tại miền Trung Việt Nam (từ Quảng Bình đến Ninh Thuận).
• Về mặt chính sách quản lý: Cung cấp luận cứ khoa học chuẩn xác phục vụ việc triển khai Thông tư 76/2017/TT-BTNMT tại tỉnh Thừa Thiên Huế; định lượng hạn ngạch cấp phép xả thải (TMAL) cho từng tiểu lưu vực sông Hương, sông Bồ và sông Ô Lâu.
• Về mặt ứng dụng thực tiễn: Khuyến nghị tái cấu trúc quy hoạch nuôi trồng thủy sản đầm Sam Chuồn và Cầu Hai, chuyển đổi từ nuôi tôm thâm canh mật độ cao sang mô hình đồng quản lý nuôi sinh thái đa loài kết hợp bảo tồn thảm cỏ biển và rừng ngập mặn nhằm nâng cao năng lực lọc sinh học tự nhiên.

Limitations và Future Research

Luận án thừa nhận 4 giới hạn nghiên cứu xuất phát từ điều kiện nguồn lực và kỹ thuật:

1. Dữ liệu chuỗi thời gian thực địa: Mặc dù các đợt khảo sát 24 giờ liên tục phản ánh được biến thiên ngày - đêm của các thông số, việc thiếu mạng lưới trạm quan trắc tự động liên tục (Real-time telemetry sensors) nhiều năm khiến việc mô phỏng các biến cố thời tiết cực đoan (như bão lịch sử, lũ quét đột biến) còn phụ thuộc vào dữ liệu nội suy.
2. Thông lượng trao đổi dinh dưỡng nền đáy (Benthic Flux): Quá trình giải phóng nitơ và phospho từ lớp trầm tích đáy vào cột nước trong điều kiện thiếu khí đáy mới chỉ được tham số hóa thông qua các hệ số động học tổng quan trong module Delft3D-WAQ chứ chưa được đo đạc trực tiếp bằng các buồng ủ đáy (Benthic chambers) tại hiện trường.
3. Phạm vi thông số ô nhiễm: Nghiên cứu mới tập trung chủ yếu vào chu trình C, N, P; các chất ô nhiễm vi lượng như kim loại nặng (Cu, Pb, Zn, Cd, As, Hg) và hóa chất bảo vệ thực vật cơ clo (Lindan, Aldrin, Endrin, DDT) dù đã được khảo sát sơ bộ nhưng chưa được đưa vào mô hình lan truyền sức tải động học.
4. Biến đổi khí hậu và nước biển dâng: Các kịch bản mô phỏng 2020 và 2030 chưa tích hợp đầy đủ biến số nước biển dâng và hiện tượng thay đổi cửa bồi lấp/mở mới tại cửa biển Thuận An và Tư Hiền.

Chương trình nghiên cứu tiếp nối (Future Research Agenda) cần tập trung vào 4 hướng chiến lược:

• Xây dựng mô hình ghép nối thủy động lực - trầm tích - sinh địa hóa 3D phân giải siêu cao tích hợp thông lượng chất đáy (Benthic nutrient fluxes) đo đạc thực tế.
• Đánh giá sức tải môi trường đối với ô nhiễm vi nhựa (Microplastics) và dư lượng kháng sinh từ hoạt động nuôi tôm công nghiệp.
• Ứng dụng công nghệ học máy (Machine Learning / Deep Neural Networks) kết hợp với Delft3D để dự báo cảnh báo sớm nguy cơ suy giảm oxy hòa tan và bùng phát tảo độc thời gian thực.
• Nghiên cứu tác động kinh tế sinh thái (Ecological-Economic modeling) nhằm định giá dịch vụ hấp thụ carbon và đồng hóa chất ô nhiễm của hệ sinh thái đầm phá Tam Giang - Cầu Hai.

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

• Ảnh hưởng học thuật (Academic Impact): Luận án xác lập tài liệu tham khảo chuẩn mực đầu tiên tại Việt Nam về phương pháp luận tính toán sức tải đầm phá quy mô lớn, tạo tiền đề cho hàng loạt công trình nghiên cứu sinh tiến sĩ và các đề tài cấp Nhà nước về hải dương học và môi trường biển.
• Chuyển đổi ngành Thủy sản (Industry Transformation): Kết quả nghiên cứu trực tiếp cung cấp cơ sở dữ liệu giúp ngành nông nghiệp Thừa Thiên Huế kiểm soát mật độ lồng bè cá và diện tích đầm nuôi tôm sú, thúc đẩy chuyển dịch sang mô hình nuôi trồng thủy sản tuần hoàn và an toàn sinh học.
• Hoạch định chính sách (Policy Influence): Là cơ sở khoa học then chốt hỗ trợ Sở Tài nguyên và Môi trường tỉnh Thừa Thiên Huế xây dựng "Quy hoạch phân vùng quản lý tổng hợp vùng bờ" và "Kế hoạch hành động bảo vệ môi trường đầm phá Tam Giang - Cầu Hai tầm nhìn 2030".
• Lợi ích xã hội và môi trường toàn cầu: Bảo tồn trực tiếp hệ sinh thái đất ngập nước có tầm quan trọng quốc tế, bảo vệ sinh kế bền vững cho hơn 300.000 cư dân thuộc 33 xã/thị trấn sinh sống ven đầm phá, đóng góp vào mục tiêu phát triển bền vững của Liên Hợp Quốc (UN SDG 14: Life Below Water).

Đố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 mẫu mực về tích hợp mô hình phát thải kinh tế xã hội với mô hình số trị thủy động lực chất lượng nước Delft3D; tiếp nhận các bộ tham số động học sinh địa hóa đã được hiệu chỉnh cho vùng nước ven bờ nhiệt đới.
• Cơ quan Quản lý Nhà nước (Bộ TN&MT, Chi cục BVMT Thừa Thiên Huế): Sở hữu công cụ định lượng chính xác để thực hiện cấp phép xả thải theo hạn ngạch TMDL, kiểm soát ô nhiễm liên huyện và phân bổ kinh phí xử lý môi trường trúng đích.
• Doanh nghiệp và Cộng đồng Nuôi trồng Thủy sản: Nhận biết rõ các rủi ro độc học môi trường theo mùa, từ đó tối ưu hóa lịch thời vụ thả giống, giảm thiểu thiệt hại dịch bệnh do thiếu oxy hòa tan hoặc ngộ độc amoniac.
• Các Tổ chức Quốc tế và Ban Quản lý Khu Bảo tồn: Sử dụng dữ liệu sức tải để thiết kế các dự án phục hồi rừng ngập mặn, rạn san hô, thảm cỏ biển và thúc đẩy các sáng kiến tài chính xanh / tín chỉ đa dạng sinh học.

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: Đóng góp độc đáo nhất là việc mở rộng Lý thuyết Sức tải Môi trường của GESAMP (1986) và Lý thuyết Cân bằng Ô nhiễm của Vlađimirov (1982) từ dạng công thức hộp tĩnh sang mô hình động học không gian 2D/3D gắn liền với ngưỡng độc học sinh thái ($DO \ge 5\text{ mg/L}, NH_3 < 0,05\text{ mg/L}, P\text{-}PO_4^{3-} < 0,045\text{ mg/L}$). Luận án chứng minh rằng sức tải không phải là một hằng số cố định theo thể tích thủy vực mà là một hàm biến thiên động theo chu kỳ trao đổi triều và tốc độ chuyển hóa sinh địa hóa tự nhiên.

2. Điểm cải tiến phương pháp luận vượt bậc so với các nghiên cứu trước đây là gì? Trả lời: So với nghiên cứu của Võ Duy Sơn và Nguyễn Tác An (2001) chỉ dùng đồng vị vết đo trao đổi nước đơn giản, hay nghiên cứu IMOLA (2005) chỉ dùng mô hình Ecopath tĩnh cho đầm Sam Chuồn, luận án đã liên kết toàn diện chuỗi dữ liệu phát thải lục địa đa nguồn (sinh hoạt, chăn nuôi, công nghiệp, NTTS, rửa trôi) với bộ mô hình Delft3D-Flow và Delft3D-WAQ trên lưới cong trực giao, giải quyết đồng thời bài toán thủy lực phân tầng 2 chiều và 21 phản ứng sinh địa hóa trong cột nước.

3. Phát hiện bất ngờ và phản trực giác nhất từ dữ liệu thực nghiệm là gì? Trả lời: Phát hiện bất ngờ nhất là nồng độ amoni tổng số ($N\text{-}NH_4^+$) luôn nằm trong giới hạn cho phép của QCVN 10-MT:2015/BTNMT, nhưng vào mùa khô nắng nóng, sự kết hợp giữa nhiệt độ cao ($33^\circ\text{C}$) và pH kiềm hóa ($8,5 - 9,0$) đã chuyển dịch cân bằng ion hóa, sinh ra hàm lượng khí $NH_3$ tự do vượt ngưỡng độc học sinh thái $0,05\text{ mg/L}$, tạo thành "bẫy độc học vô hình" đối với sinh vật thủy sinh mà hệ thống quan trắc truyền thống bỏ sót.

4. Luận án có cung cấp đầy đủ giao thức tái lập (Replication Protocol) không? Trả lời: Có. Luận án mô tả chi tiết toàn bộ thông số thiết lập lưới tính, điều kiện biên thủy văn các sông, hằng số phân hủy động học ($BOD, COD, N, P$), các hệ số phát thải đơn vị chuẩn hóa theo WHO và San Diego - McGlone, cùng giao thức quan trắc hiện trường 4h/lần theo TCVN 5998:1995 và chuẩn phân tích HACH DR/2800, cho phép các nhà khoa học khác tái lập hoàn toàn quy trình trên các hệ đầm phá khác.

5. Chương trình nghị sự nghiên cứu 10 năm tới được phác thảo như thế nào? Trả lời: Lộ trình 10 năm tập trung vào: (1) Thiết lập hệ thống cảm biến thông minh IoT quan trắc tự động online kết hợp mô hình Delft3D đồng hóa dữ liệu thời gian thực; (2) Tích hợp mô hình biến đổi khí hậu và tác động của việc vận hành các hồ chứa thủy điện thượng nguồn sông Hương (Tả Trạch, Hương Điền, Bình Điền); (3) Mở rộng tính toán sức tải sang các chất ô nhiễm mới nổi (vi nhựa, dược phẩm, kim loại nặng dạng keo).

Kết luận

1. Luận án đã hệ thống hóa và làm sáng tỏ cơ sở lý luận về sức tải môi trường thủy vực ven bờ, chuyển dịch mô thức từ đánh giá định tính sang lượng hóa số trị phân giải cao dựa trên động học tự làm sạch tự nhiên.
2. Thiết lập thành công bộ công cụ mô hình hóa Delft3D mô phỏng chính xác chế độ thủy động lực và lan truyền chất ô nhiễm cho toàn bộ 216 km² hệ đầm phá Tam Giang - Cầu Hai trong cả mùa mưa và mùa khô.
3. Định lượng toàn diện tải lượng ô nhiễm phát sinh từ 5 nguồn thải lục địa chính và dự báo xu thế biến động đến năm 2020 và 2030 theo các kịch bản phát triển kinh tế xã hội.
4. Lần đầu tiên xác định hệ thống ngưỡng sức tải môi trường đa cấp độ cho các thông số $BOD_5, COD, N\text{-}NH_4^+, N\text{-}NO_3^-$ và $P\text{-}PO_4^{3-}$, đặc biệt nhấn mạnh ngưỡng an toàn độc thái học sinh vật thủy sinh.
5. Cung cấp luận cứ khoa học trực tiếp và công cụ phân vùng quản lý phục vụ đắc lực cho công tác bảo vệ môi trường, kiểm soát cấp phép xả thải của tỉnh Thừa Thiên Huế và mở ra hướng đi mẫu mực cho công tác quản trị vùng bờ bền vững tại Việt Nam.

Trích đoạn nội dung tài liệu

Tính cấp thiết của đề tài Hệ đầm phá Tam Giang - Cầu Hai (thuộc tỉnh Thừa Thiên Huế) là đầm phá lớn nhất ở ven bờ Việt Nam với gần 1000 loài động vật, thực vật thủy sinh có giá trị kinh tế [1]. Các hoạt động kinh tế - xã hội vùng đầm phá đang diễn ra hết sức sôi động bao gồm nông nghiệp, nghề cá và khai thác biển, giao thông - cảng, du lịch – dịch vụ v. Hệ đầm phá là nơi tiếp nhận các nguồn thải ven bờ không những của các huyện giáp ranh mà còn cả của các khu vực miền núi. Khả năng suy thoái chất lượng môi trường, cạn kiệt nguồn giống sẽ xảy ra nếu không có những biện pháp quản lý hệ thống đầm phá.

Mỗi một hệ thống tự nhiên có một khả năng chịu tải nhất định. Vượt quá ngưỡng đó, hệ thống sẽ bị thay đổi kéo theo sự thay đổi chức năng của hệ thống. Trong khi đó, các hoạt động phát triển kinh tế - xã hội ven bờ đã dẫn đến tải lượng hữu cơ và dinh dưỡng đưa vào hệ đầm phá không ngừng gia tăng mà không có biện pháp bảo vệ hoặc cảnh báo. Trước sức ép phát triển kinh tế của khu vực, nghiên cứu sinh đã chọn đề tài “Nghiên cứu đánh giá sức tải một số yếu tố môi trường (C, N, P) khu vực đầm phá Tam Giang - Cầu Hai (tỉnh Thừa Thiên Huế)” làm luận án nghiên cứu của mình.

Do các nguồn thải đưa vào hệ đầm phá Tam Giang - Cầu Hai (TG – CH) chủ yếu là các chất thải sinh hoạt từ dân cư, khách du lịch, chăn nuôi và nuôi trồng thủy sản (NTTS) trong vùng nên các yếu tố môi trường được chọn để đánh giá sức chịu tải là các chất hữu cơ (BOD5, COD), các chất dinh dưỡng (amoni, phosphat, nitrat, nitrit, T-N, P-P). Đây là các hợp chất có mặt thường xuyên trong nước thải sinh hoạt, chăn nuôi và NTTS. Sức tải môi trường là một hướng tiếp cận có tính phòng ngừa trong đó khả năng chịu tải môi trường của mỗi hệ thống tự nhiên là hữu hạn. Khi các chất ô nhiễm đi vào các thủy vực ven bờ, chúng sẽ tham gia vào các quá trình tự làm sạch tự nhiên bao gồm hóa học, lý học và sinh học [2].

Nếu lượng chất ô nhiễm đi vào thủy vực lớn hơn khả năng tự làm sạch tự nhiên của nó, thủy vực sẽ bị thay đổi về cấu trúc và chức năng, thậm chí mất khả năng tự phục hồi. Nghiên cứu sức tải môi trường có một số hướng tiếp cận và được áp dụng trong một số lĩnh vực như trong NTTS, đánh giá sức 1 tải du lịch, quản lý nguồn thải, quản lý hệ sinh thái. Hướng tiếp cận của luận án tập trung vào nghiên cứu, quản lý nguồn thải, góp phần bảo vệ chất lượng môi trường nước và hệ sinh thái. Hiểu và đánh giá đúng sức chịu tải môi trường có ý nghĩa quan trọng trong việc đưa ra các chính sách phát triển kinh tế - xã hội và bảo vệ môi trường.

Mục tiêu của luận án - Đánh giá được sức chịu tải của một số yếu tố môi trường (C, N, P) trong hệ đầm phá TG - CH làm cơ sở cho quản lý, phát triển bền vững hệ đầm phá. Nội dung nghiên cứu - Phân tích, đánh giá và dự báo tải lượng ô nhiễm từ các nguồn đưa vào hệ đầm phá. - Mô phỏng lan truyền các chất ô nhiễm trong hệ đầm phá TG - CH theo các kịch bản cơ sở (năm 2011 – 2012) và kịch bản 2020, 2030, kịch bản đột xuất. - Nghiên cứu, tính toán sức tải hệ đầm phá TG - CH đối với các chất hữu cơ và chất dinh dưỡng theo các ngưỡng của quy chuẩn Việt Nam, ngưỡng sức tải tối đa và ngưỡng gây bất lợi đối với sinh vật thủy sinh.

Đối tượng và phạm vi nghiên cứu - Đối tượng nghiên cứu: các hợp chất C, N, P trong nước hệ đầm phá TG - CH; các nguồn thải đưa vào hệ thống đầm phá. - Phạm vi không gian: không gian nghiên cứu là hệ đầm phá TG – CH và các vùng xung quanh đưa các chất ô nhiễm vào đầm phá. - Phạm vi thời gian: mùa mưa (tháng 11) và mùa khô (tháng 5) của các năm 2011, 2012, 2016; dự báo cho các năm 2020, 2030 5. Phương pháp nghiên cứu - Nghiên cứu tổng quan, hồi cứu tài liệu: Thu thập các tài liệu đã công bố trên thế giới và trong nước về sức tải môi trường; các tài liệu về đặc điểm tự nhiên, kinh tế- xã hội, chất lượng môi trường hệ đầm phá TG – CH.

- Điều tra, khảo sát và thực nghiệm ngoài hiện trường: đo đạc dòng chảy, thu và phân tích mẫu chất lượng nước tại hệ đầm phá TG – CH. - Mô hình hóa: mô phỏng chế độ thủy động lực và sự lan truyền chất ô nhiễm trong hệ đầm phá sử dụng phần mềm Delft – 3D. Ý nghĩa khoa học và thực tiễn của đề tài 2 - Hệ thống hóa được phương pháp đánh giá sức chịu tải môi trường cho một thủy vực ven bờ Việt Nam. - Góp phần xác định sức tải môi trường các yếu tố BOD5, COD, N-NH4+, N- NO3- và P-PO43- cho hệ đầm phá Tam Giang – Cầu Hai, có thể sử dụng làm nguồn tham khảo cho công tác quản lý môi trường.

Những kết quả khoa học đạt được và đóng góp mới của luận án - Đã xác định nguồn thải và ước tính lượng thải các chất ô nhiễm C, N, P từ các hoạt động kinh tế - xã hội đưa vào hệ đầm phá TG-CH. - Đã hiệu chỉnh mô hình chất lượng nước hệ đầm phá TG-CH và mô phỏng chất lượng nước theo các kịch bản phát triển đến năm 2020 và 2030. - Đã tính được sức tải môi trường cho hệ đầm phá TG-CH theo các ngưỡng của QCVN và theo khả năng tự làm sạch (đồng hóa) của đầm phá. TỔNG QUAN VẤN ĐỀ NGHIÊN CỨU VÀ KHU VỰC NGHIÊN CỨU 1.

Tổng quan về sức tải môi trường 1. Các khái niệm Sức tải môi trường (STMT) là một hướng nghiên cứu của khoa học môi trường, mới được phát triển trong hơn ba thập kỷ gần đây. Bắt đầu từ cuối những năm 1980, Nhóm các chuyên gia về khoa học ô nhiễm biển (GESAMP) [3] đã tìm hiểu về khái niệm, cách tiếp cận liên quan đến STMT. Sau này, với việc phát triển của các công cụ mô hình hóa, hướng nghiên cứu này đã được áp dụng tại một số lĩnh vực như trong nuôi trồng thủy sản, du lịch, kiểm soát nguồn thải v.

Đặc biệt, trong khoảng 15 năm trở lại đây, dưới sức ép của các quy hoạch phát triển và sự suy thoái của tài nguyên và môi trường, hướng nghiên cứu STMT được phát triển và áp dụng khá rộng rãi. Vấn đề đặt ra là cần làm rõ các khái niệm, hiểu đúng bản chất của STMT và có phương pháp tính toán đúng đắn để có thể áp dụng vào trong thực tiễn quản lý nguồn thải và bảo vệ môi trường. Ngoài khái niệm về STMT, các khái niệm về khả năng tự làm sạch và cân bằng chất ô nhiễm trong thủy vực cũng có mối liên quan chặt chẽ với nhau, bổ sung cho nhau. Khả năng tự làm sạch của thủy vực: Các chất ô nhiễm khi đi vào các thủy vực tự nhiên như sông, hồ và biển sẽ bị biến mất dần dần theo thời gian.

Việc loại bỏ các chất ô nhiễm từ các thủy vực mà không có bất kỳ can thiệp nào của con người được gọi là quá trình tự làm sạch, hoặc sự làm sạch tự nhiên [4]. Cơ chế của quá trình tự làm sạch tự nhiên có thể chia làm 3 nhóm: vật lý, hóa học và sinh học [2]. Quá trình vật lý đóng góp vào việc loại bỏ chất ô nhiễm từ thủy vực tự nhiên bao gồm pha loãng/trộn lẫn bằng các dòng chảy, khuyếch tán chất ô nhiễm trong nước và kết tủa/lắng đọng chất ô nhiễm xuống trầm tích đáy. Sự bay hơi của các chất ô nhiễm dễ bay hơi vào bầu khí quyển cũng làm giảm chất ô nhiễm trong nước.

Quá trình hóa học liên quan đến việc loại bỏ các chất ô nhiễm khỏi vực nước là sự oxi hóa bởi các chất ô xi hóa như ozon và oxy, sự oxi hóa bởi tia tử ngoại; quá trình khử và sự trung hòa. 4 Quá trình sinh học bao gồm sự thoái hóa/chuyển đổi các chất ô nhiễm hữu cơ bởi vi khuẩn dưới các điều kiện hiếu khí hoặc yếm khí, cũng như sự tham gia của các vi khuẩn trong quá trình nitrat hóa và phản nitrat hóa tương ứng. Quá trình sinh học giữ một vai trò quan trọng trong cơ chế làm sạch nước ở các thủy vực tự nhiên. Việc loại bỏ chất ô nhiễm khỏi thủy vực bằng phương pháp sinh học thường được gọi là “tự làm sạch đúng nghĩa” và sự tự làm sạch tổng số “lý, hóa, sinh” được gọi là “tự làm sạch biểu kiến”.

Cân bằng chất ô nhiễm trong thủy vực: Nghiên cứu về cân bằng các chất ô nhiễm cung cấp cơ sở để dự báo động thái ô nhiễm của các vùng nước biển và đại dương theo những giá trị cho trước về phát thải các chất ô nhiễm, xây dựng những khuyến cáo về chế độ phát thải tối ưu và xác định các mức chịu tải cho phép tới hạn hay dung lượng dung hòa của biển với các chất ô nhiễm. Vlađimirov trong cuốn “Bảo vệ Môi trường” [5], về định tính, cân bằng các chất ô nhiễm đối với biển nói chung có thể biểu diễn bằng sơ đồ sau: ∆C = C vào – C ra Ở đây ∆C là gia lượng nồng độ chất ô nhiễm sau thời gian t, C vào là phát thải chất ô nhiễm vào thủy vực (bao gồm phát thải từ bờ, phát thải do sông mang ra, phát thải trực tiếp từ biển như tàu, giàn khoan v.v; lượng chất tới trong quá trình trao đổi nước, lượng chất xâm nhập từ khí quyển vào nước) và C ra là lượng chất ô nhiễm ra khỏi thủy vực trong thời gian t (bao gồm lượng chất đi do trao đổi nước, lượng chất đi từ biển vào không khí, lượng chất bị phân hủy hóa học, sinh hóa và lượng chất lắng đọng vào bùn đáy). Sự bảo tồn trạng thái bình thường của môi trường biển đòi hỏi phải thỏa mãn những điều kiện hạn chế: 0< C ≤ NĐTHCP và ∆C < 0 trong đó NĐTHCP là nồng độ tới hạn cho phép của các chất ô nhiễm.