Luận án tiến sĩ nông nghiệp nghiên cứu mô phỏng quá trình chuyển hóa các chất dinh dưỡng chủ yếu tác động đến môi trường ao nuôi cá tra việt nam

Luận án tiến sĩ nông nghiệp nghiên cứu mô phỏng quá trình chuyển hóa chất dinh dưỡng ảnh hưởng đến môi trường ao nuôi cá tra tại Việt Nam.

• Tác giả: Lê Xuân Thịnh
• Người hướng dẫn: PGS.TS Đặng Xuân Hiển
• Trường học: Trường Đại học Bách Khoa Hà Nội
• Chuyên ngành: Khoa học và Công nghệ môi trường
• Đề tài: Nghiên cứu mô phỏng chuyển hóa chất dinh dưỡng tác động môi trường ao nuôi cá tra Việt Nam
• Loại tài liệu: luận án
• Năm xuất bản: 2022
• Địa điểm: Hà Nội

Mục lục chi tiết

• LỜI CẢM ƠN
• LỜI CAM ĐOAN
• DANH MỤC CÁC KÝ HIỆU, CÁC CHỮ VIẾT TẮT
• DANH MỤC CÁC HÌNH ẢNH
• DANH MỤC CÁC BẢNG
• MỞ ĐẦU
• 1. CHƯƠNG 1: TỔNG QUAN
• 1.1. Tổng quan về cá Tra
• 1.1.1. Đặc điểm hình thái
• 1.1.2. Đặc điểm phân bố
• 1.1.3. Đặc điểm sinh trưởng
• 1.1.4. Dinh dưỡng và thức ăn cá Tra
• 1.1.4.1. Cân bằng dinh dưỡng trong ao nuôi
• 1.1.4.2. Dinh dưỡng và thức ăn cá Tra
• 1.1.5. Quản lý chất lượng nước
• 1.1.6. Phòng trị bệnh
• 1.1.7. Thu hoạch cá Tra
• 1.2. Các vấn đề môi trường trong ao nuôi cá Tra
• 1.2.1. Các thông số thủy lý
• 1.2.1.1. Tổng chất rắn lơ lửng (TSS)
• 1.2.2. Các thông số thủy hóa trong ao nuôi
• 1.2.2.1. Tổng nitơ amôn
• 1.2.2.2. Tổng phốt pho
• 1.2.3. Tính chất của bùn đáy ao nuôi cá Tra
• 1.2.3.1. Nồng độ chất hữu cơ
• 1.2.3.2. Tổng phốt pho (TP)
• 1.2.3.3. Biến động nồng độ Ca và Mg (mg/kg)
• 1.2.3.4. Nồng độ các nguyên tố vi lượng – kim loại nặng
• 1.2.4. Tảo và các chất hữu cơ trong ao
• 1.2.5. Các chu trình sinh học trong ao
• 1.2.5.1. Chu trình nitơ
• 1.2.5.2. Chu trình phốt pho
• 1.2.5.3. Chu trình các bon
• 1.2.6. Các quá trình chuyển hóa trong bùn đáy ao
• 1.3. Enzyme trong thức ăn thủy sản
• 1.3.1. Vai trò và tác dụng của enzyme phytase trong thức ăn thủy sản
• 1.3.2. Các nghiên cứu ứng dụng của enzyme phytase trong nuôi trồng thủy sản
• 1.4. Mô hình và mô hình hóa
• 1.4.1. Mô hình hóa trong nghiên cứu môi trường
• 1.4.2. Các thành phần trong quá trình mô hình hóa môi trường
• 1.4.3. Mô hình hóa chất lượng môi trường nước ao, hồ
• 1.4.3.1. Các tính chất đặc trưng của hồ
• 1.4.3.2. Các quá trình chuyển đổi vật chất trong hồ
• 1.4.4. Các nghiên cứu về mô phỏng mô hình lưới thức ăn
• 1.4.4.1. Giới thiệu chung mô hình lưới thức ăn
• 1.4.4.2. Tính toán, mô phỏng các mối quan hệ trong lưới thức ăn trong ao nuôi cá
• 1.4.5. Một số mô hình được nghiên cứu trên thế giới
• 1.4.5.1. Mô hình PCLake nghiên cứu về hiện tượng phú dưỡng trong hồ nông
• 1.4.5.2. Nghiên cứu về mô hình hóa chu trình photpho gây phú dưỡng hồ Loosdrecht
• 1.4.5.3. Nghiên cứu mô hình động lực dinh dưỡng trong ao nuôi cá bán thâm canh
• 1.4.5.4. Nghiên cứu mô hình phú dưỡng hồ Washington
• 1.4.5.5. Mô phỏng các chất dinh dưỡng và động học của thực vật phù du trong hồ Bắc Mỹ
• 1.4.6. Các nghiên cứu về xây dựng mô hình mô phỏng trong nuôi trồng thủy sản ở Việt Nam
• 2. CHƯƠNG 2: PHƯƠNG PHÁP NGHIÊN CỨU, THÍ NGHIỆM
• 2.1. Phương pháp nghiên cứu
• 2.2. Thí nghiệm, thu mẫu và phân tích các thông số môi trường ao nuôi
• 2.2.1. Bố trí thí nghiệm
• 2.2.2. Phương pháp thu mẫu và phân tích
• 2.2.2.1. Đánh giá hiệu quả sử dụng thức ăn
• 2.2.2.2. Phương pháp thu mẫu và phân tích
• 2.3. Thiết lập mô hình
• 2.3.1. Cơ sở thiết lập các phương trình
• 2.3.1.1. Mô hình các quá trình sinh học trong đầm hồ và đầm hồ nuôi cá Tra
• 2.3.1.2. Mô hình hóa sự phát triển của cá Tra
• 2.3.2. Mô hình khái niệm ao nuôi cá Tra
• 2.3.2.1. Các giả thiết của mô hình
• 2.3.2.2. Mô hình khái niệm
• 2.3.2.3. Ma trận tác động tương hỗ
• 2.3.2.4. Cơ sở thiết lập mối quan hệ toán học giữa các cấu tử và các quá trình sinh học xảy ra trong ao nuôi
• 2.3.3. Xây dựng các phương trình toán học
• 2.3.3.1. Phương trình sự phát triển của cá trong hồ nuôi cá Tra
• 2.3.3.2. Các phương trình tính toán các thông số môi trường
• 2.3.4. Giải hệ phương trình vi phân
• 2.3.5. Giải bằng phương pháp số trên cơ sở ngôn ngữ máy tính
• 2.4. Phương pháp phân tích độ nhạy và hiệu chỉnh mô hình
• 2.4.1. Phân tích độ nhạy
• 3. CHƯƠNG 3: KẾT QUẢ VÀ THẢO LUẬN
• 3.1. Ảnh hưởng của enzyme phytase đến môi trường trong ao nuôi
• 3.1.1. Biến động nhiệt độ, pH, DO giữa đối chứng và thí nghiệm
• 3.1.2. Biến động của BOD và COD
• 3.1.3. Biến động của TSS
• 3.1.4. Biến động của TAN, NH3 và NH4+
• 3.1.5. Biến động N-NO2-, N-NO3-
• 3.1.6. Biến động P-PO43- và TP
• 3.1.7. Biến động TN, TP bùn đáy ao
• 3.1.8. Đánh giá ảnh hưởng của bổ sung enzym phytase lên hiệu quả sử dụng thức ăn và môi trường ao nuôi cá Tra
• 3.1.8.1. Đánh giá tốc độ tăng trưởng, tỉ lệ sống và hiệu quả sử dụng thức ăn của cá Tra ở thí nghiệm
• 3.1.8.2. Hiệu quả của bổ sung enzyme phytase lên môi trường
• 3.2. Số liệu chạy mô hình
• 3.3. Kết quả phân tích độ nhạy và hiệu chỉnh mô hình
• 3.3.1. Phân tích độ nhạy
• 3.4. Mô phỏng quá trình phát triển động thực vật phù du trong ao cá
• 3.5. Mô phỏng quá trình phát triển và chuyển hóa các chất dinh dưỡng chủ yếu trong ao
• 3.5.1. Quá trình phát triển của cá
• 3.5.1.1. Quá trình phát triển của cá tại ao nuôi thực địa
• 3.5.1.2. Mô phỏng quá trình phát triển của cá trong ao nuôi
• 3.5.2. Biến động ôxy hòa tan trong ao nuôi
• 3.5.2.1. Biến động ôxy hòa tan trong ao nuôi thực địa
• 3.5.2.2. Mô phỏng quá trình biến thiên của ôxy trong ao nuôi
• 3.5.3. Tổng nitơ và tổng phốt pho trong cột nước
• 3.5.3.1. Tổng nitơ và tổng phốt pho trong cột nước trong ao nuôi thực địa
• 3.5.3.2. Mô phỏng tổng nitơ và tổng phốt pho trong cột nước trong ao nuôi cá Tra
• 3.5.4. Tổng nitơ và tổng phốt pho trong bùn lắng
• 3.5.4.1. Tổng nitơ và tổng phốt pho trong bùn lắng ao nuôi thực địa
• 3.5.4.2. Mô phỏng Tổng nitơ và tổng phốt pho trong bùn lắng trong ao nuôi thực địa
• KẾT LUẬN VÀ KIẾN NGHỊ
• TÀI LIỆU THAM KHẢO
• DANH MỤC CÁC CÔNG TRÌNH ĐÃ CÔNG BỐ CỦA LUẬN ÁN
• SỐ LIỆU MÔI TRƯỜNG AO NUÔI CÁ TRA Ở VIỆT NAM
• MỘT SỐ HÌNH ẢNH THÍ NGHIỆM AO NUÔI THỰC ĐỊA
• MỘT SỐ HÌNH ẢNH GIAO DIỆN KẾT QUẢ CHẠY MÔ HÌNH
• BẢNG GIÁ TRỊ THAM SỐ VÀ CODE SỐ

Tóm tắt

I. Mở đầu

Nghiên cứu về chuyển hóa chất dinh dưỡng trong ao nuôi cá Tra (Pangasianodon hypophthalmus) tại Việt Nam đã chỉ ra rằng loài cá này có khả năng thích nghi tốt với điều kiện môi trường khắc nghiệt. Tuy nhiên, việc nuôi cá Tra cũng gây ra nhiều vấn đề về tác động môi trường. Theo thống kê, mỗi tấn cá Tra thành phẩm thải ra khoảng 4.023m3 nước, trong đó một phần lớn là chất thải từ thức ăn dư thừa. Điều này dẫn đến ô nhiễm nguồn nước và ảnh hưởng đến sức khỏe của cá. Nghiên cứu này nhằm mục đích mô phỏng và phân tích các quá trình chuyển hóa chất dinh dưỡng và tác động sinh thái trong ao nuôi cá Tra, từ đó đưa ra các giải pháp quản lý bền vững.

II. Đặc điểm sinh học và dinh dưỡng của cá Tra

Cá Tra có đặc điểm sinh học nổi bật với khả năng tăng trưởng nhanh và phổ dinh dưỡng rộng. Nghiên cứu cho thấy, khi cho cá ăn 100% lượng thức ăn công nghiệp, chỉ khoảng 37,5% lượng dinh dưỡng được tích lũy trong cơ thể cá. Phần còn lại chủ yếu là chất thải, dẫn đến sự tích lũy chất dinh dưỡng trong bùn đáy ao. Dinh dưỡng và thức ăn cá Tra là yếu tố quan trọng trong việc quản lý chất lượng nước. Việc quản lý hiệu quả thức ăn không chỉ giúp tăng trưởng cá mà còn giảm thiểu ô nhiễm môi trường. Các nghiên cứu đã chỉ ra rằng, việc bổ sung enzyme phytase vào thức ăn có thể cải thiện hiệu quả sử dụng thức ăn và giảm thiểu tác động tiêu cực đến môi trường.

III. Các vấn đề môi trường trong ao nuôi cá Tra

Môi trường ao nuôi cá Tra đang phải đối mặt với nhiều thách thức, bao gồm ô nhiễm do chất thải và thức ăn dư thừa. Các thông số như tổng chất rắn lơ lửng (TSS), tổng nitơ (TN) và tổng phốt pho (TP) đều có xu hướng gia tăng. Nghiên cứu cho thấy, nồng độ nitơ và phốt pho trong nước và bùn đáy ao nuôi cá Tra cao hơn mức cho phép, gây ra hiện tượng phú dưỡng. Điều này không chỉ ảnh hưởng đến chất lượng nước mà còn làm giảm năng suất và sức khỏe của cá. Việc quản lý chất lượng nước và giảm thiểu ô nhiễm là rất cần thiết để đảm bảo tính bền vững của nghề nuôi cá Tra.

IV. Mô phỏng và mô hình hóa trong nghiên cứu môi trường

Mô phỏng các quá trình chuyển hóa chất dinh dưỡng trong ao nuôi cá Tra là một công cụ quan trọng để hiểu rõ hơn về các tác động môi trường. Các mô hình hóa chất lượng nước ao, hồ giúp dự đoán các biến động trong nồng độ chất dinh dưỡng và ô nhiễm. Nghiên cứu đã áp dụng mô hình PCLake để phân tích hiện tượng phú dưỡng trong ao nuôi. Kết quả cho thấy, việc mô phỏng không chỉ giúp đánh giá tình trạng môi trường mà còn cung cấp thông tin hữu ích cho việc quản lý và cải thiện chất lượng nước trong ao nuôi cá Tra.

V. Kết luận và kiến nghị

Nghiên cứu về chuyển hóa chất dinh dưỡng và tác động môi trường trong ao nuôi cá Tra đã chỉ ra rằng, việc quản lý chất lượng nước và thức ăn là rất quan trọng để đảm bảo tính bền vững của nghề nuôi cá. Cần có các giải pháp hợp lý để giảm thiểu ô nhiễm và cải thiện hiệu quả sử dụng thức ăn. Đề xuất các nghiên cứu tiếp theo nên tập trung vào việc phát triển các mô hình quản lý bền vững, nhằm bảo vệ môi trường và nâng cao năng suất nuôi cá Tra.

Nội dung chính

Tổng quan luận án

Luận án tiến sĩ kỹ thuật với đề tài "Nghiên cứu mô phỏng quá trình chuyển hóa các chất dinh dưỡng chủ yếu tác động đến môi trường ao nuôi cá Tra Việt Nam" do nghiên cứu sinh Lê Xuân Thịnh thực hiện dưới sự hướng dẫn khoa học của PGS.TS. Đặng Xuân Hiển và PGS.TS. Trần Văn Nhân, tại Viện Khoa học và Công nghệ Môi trường – Trường Đại học Bách Khoa Hà Nội, hoàn thành vào tháng 12 năm 2022. Nghiên cứu nhận được sự hỗ trợ kinh phí thực địa từ Liên minh Châu Âu (EU) thông qua Chương trình Switch-Asia trong khuôn khổ dự án "Xây dựng chuỗi cung ứng cá tra bền vững ở Việt Nam – SUPA", cùng sự phối hợp phân tích từ Khoa Thủy sản – Trường Đại học Cần Thơ.

Tính cấp thiết và khoảng trống nghiên cứu

Cá Tra (Pangasianodon hypophthalmus) là đối tượng nuôi trồng thủy sản nước ngọt chủ lực tại Đồng bằng sông Cửu Long (ĐBSCL), tập trung tại các tỉnh Cần Thơ, An Giang, Đồng Tháp, Tiền Giang, Bến Tre... Năm 2019, diện tích nuôi đạt 6.600 ha, sản lượng đạt 1,42 triệu tấn và kim ngạch xuất khẩu đạt 1,61 tỷ USD sang 131 thị trường. Tuy nhiên, phương thức nuôi thâm canh với mật độ cao (48–70 con/m²) và năng suất lớn (250–300 tấn/ha/vụ) làm phát sinh lượng chất thải khổng lồ.

Để sản xuất 1 tấn cá thương phẩm cần 4.023 m³ nước và trung bình 1,56 kg thức ăn cho mỗi kg cá. Khi sử dụng 100% thức ăn công nghiệp, cá chỉ tích lũy được khoảng 37,5% vật chất dinh dưỡng vào cơ thịt (nitơ tích lũy 32,6–42,6%; phốt pho tích lũy 29,3–29,8%), trong khi 67,4% vật chất khô bị thải ra môi trường (45,63% tích tụ ở bùn đáy, 5,04% hòa tan trong nước và 16,74% mất đi do bay hơi, rò rỉ). Một ao nuôi 300 tấn/ha/vụ thải ra khoảng 2.677 tấn bùn ướt (tương đương 937 tấn bùn khô), khiến nước ao thường xuyên vượt ngưỡng cho phép của QCVN 08:2008/BTNMT (nước mặt) và QCVN 40:2011/BTNMT (nước thải công nghiệp).

Khoảng trống nghiên cứu được xác định: Tại Việt Nam, các nghiên cứu trước đây chủ yếu dừng lại ở việc quan trắc hiện trạng chất lượng nước hoặc đánh giá thực nghiệm đơn lẻ. Rất thiếu các công trình nghiên cứu phát triển mô hình số có khả năng tính toán, mô phỏng đồng thời động học của các chất dinh dưỡng chính (C, N, P, O₂), sự phát triển của chuỗi sinh thái (thực vật phù du, động vật phù du, cá Tra) và tích lũy bùn đáy để phục vụ quản lý, giảm thiểu ô nhiễm ao nuôi.

Mục tiêu nghiên cứu

Luận án xác định các mục tiêu cụ thể gồm:

1. Đánh giá các đặc trưng chất lượng môi trường thủy lý, thủy hóa và bùn đáy của ao nuôi cá Tra thâm canh tại ĐBSCL.
2. Thiết lập bộ số liệu thực nghiệm đánh giá ảnh hưởng của việc bổ sung enzyme phytase vào thức ăn đối với sinh trưởng của cá và khả năng giảm phát thải phốt pho ra môi trường ao nuôi.
3. Xây dựng mô hình toán và chương trình số mô phỏng quá trình chuyển hóa các chất dinh dưỡng chủ yếu (C, N, P, O₂), các quá trình sinh học và tương tác môi trường trong ao nuôi cá Tra ở ĐBSCL.

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

• Đối tượng nghiên cứu: Các chất dinh dưỡng và thông số chất lượng nước (DO, BOD, COD, TSS, TAN, NH₃, NH₄⁺, N-NO₂⁻, N-NO₃⁻, P-PO₄³⁻, TN, TP) trong cột nước và bùn đáy; các quá trình sinh học và quần thể sinh vật (tảo lục, tảo lam, tảo cát, copepods, cladocerans, cá Tra).
• Phạm vi không gian: Các ao nuôi cá Tra thâm canh nước ngọt tại khu vực Đồng bằng sông Cửu Long (Cần Thơ, An Giang, Đồng Tháp).
• Phạm vi thời gian: Theo dõi chu kỳ nuôi cá Tra thương phẩm kéo dài từ 6 đến 8 tháng.

---

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

Các hướng nghiên cứu trước đây

Tác giả đã tổng hợp và điểm luận các nhóm công trình nghiên cứu trong và ngoài nước:

• Đặc tính sinh học, dinh dưỡng và thức ăn cá Tra: Các nghiên cứu về nhu cầu đạm, lipid và cân bằng năng lượng của cá Tra ở các giai đoạn phát triển (Trần Thị Thanh Hiền, 2004, 2010; Trương Hà Phương, 2004; Huỳnh Kim Khánh, 2005; dự án ACIAR, 2007). Nhu cầu protein tối ưu dao động 25–35% tùy kích cỡ cá; thức ăn công nghiệp chiếm 76–80,8% tổng chi phí nuôi.
• Cân bằng vật chất và phát thải ô nhiễm ao nuôi: Nghiên cứu của Nguyen Phu Hoa (2010), Pham Minh Duc (2010), Phan Thi Phuong Trang (2010), Boyd (1990, 1998) và Tucker & Hargreaves (2004) xác định lượng đạm và lân thất thoát vào môi trường, biến động tích lũy hữu cơ, TN, TP trong bùn đáy ao cá Tra vượt trội so với ao nuôi cá nheo Mỹ hay tôm sú.
• Ứng dụng enzyme phytase trong nuôi trồng thủy sản: Các nghiên cứu trên cá da trơn, cá trê, cá rô phi, cá chép và cá Tra của Jackson et al. (1996), Van Weerd et al. (1999), Debnath et al. (2005), Phùng Thị Hồng Gấm (2017), Nguyễn Thị Thùy (2019) chứng minh enzyme phytase giúp thủy phân phytate thực vật, giải phóng phốt pho và khoáng chất, nâng cao hệ số sử dụng protein (PER) và giảm hệ số chuyển đổi thức ăn (FCR).
• Mô hình hóa sinh thái ao hồ và chuỗi thức ăn: Các mô hình quốc tế như PCLake nghiên cứu hiện tượng phú dưỡng hồ nông (Janse, 1997, 2005); mô hình chu trình phốt pho hồ Loosdrecht; mô hình hồ Washington (Arhonditsis & Brett, 2005); mô hình cân bằng năng lượng lưới thức ăn hồ Ontario bằng phương pháp logic tích hợp (Efraim Halfon, 1996); mô hình cấu trúc dinh dưỡng đầm lầy bằng phần mềm STELLA (Douglas, 2003); phần mềm Ecopath mô phỏng cân bằng sinh khối lưới thức ăn (Polovina, 1984; Christensen & Pauly, 1992; áp dụng tại hồ Kelavarapalli năm 2009).

Khoảng trống luận án lựa chọn giải quyết

Phần lớn các mô hình sinh thái thủy vực trên thế giới được xây dựng cho các hồ tự nhiên vùng ôn đới hoặc đầm lầy bán thâm canh, không phản ánh được đặc thù của ao nuôi cá Tra thâm canh vùng nhiệt đới ĐBSCL (mật độ nuôi cực cao, độ sâu ao 3–4 m, chế độ cho ăn nhân tạo cường độ lớn, chu kỳ thay nước 30–40%/ngày và khả năng thở khí trời của cá). Luận án giải quyết khoảng trống này bằng cách thiết lập một mô hình toán học tích hợp đầy đủ các tương tác động học giữa cá Tra, thực vật phù du, động vật phù du, oxy hòa tan và hai chu trình dinh dưỡng Nitơ - Phốt pho trong cả cột nước lẫn bùn lắng.

---

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

Khung lý thuyết và mô hình khái niệm

Luận án áp dụng lý thuyết mô hình hóa môi trường sinh thái (theo phân loại của Bùi Tá Long và Jørgensen), bao gồm 5 thành phần: biến trạng thái, hàm điều khiển (biến ngoại sinh), các phương trình toán học mô tả quá trình, tham số và hằng số.

                                                     [Mảnh vụn hữu cơ / Bùn đáy ao]

Khung phân tích dựa trên nguyên lý bảo toàn vật chất (cân bằng khối lượng) cho các chu trình sinh địa hóa trong ao nước ngọt:

• Chu trình Nitơ: Các quá trình đồng hóa của tảo, bài tiết amôn (TAN/NH₄⁺), nitrit hóa tạo NO₂⁻, nitrat hóa tạo NO₃⁻, khử nitrat, thủy phân đạm hữu cơ hòa tan (DON) và phân hủy lắng đọng bùn (TN).
• Chu trình Phốt pho: Quá trình hấp thụ orthophosphate (P-PO₄³⁻), giải phóng phốt pho hữu cơ hòa tan (DOP), khoáng hóa, hấp phụ - giải hấp ở lớp bùn đáy và tích lũy TP.
• Cân bằng Oxy: Nguồn cung cấp từ quang hợp của tảo và trao đổi khí quyển; tiêu thụ oxy do hô hấp của cá, hô hấp của vi sinh vật, thực vật phù du, động vật phù du, quá trình oxy hóa sinh hóa chất hữu cơ (BOD/COD) và nhu cầu oxy của bùn đáy (SOD).

Phương pháp nghiên cứu thực nghiệm và thu thập dữ liệu

• Khảo sát thực địa: Bố trí theo dõi, thu mẫu nước và bùn đáy tại các ao nuôi cá Tra thương phẩm ở ĐBSCL theo định kỳ các tháng trong vụ nuôi.
• Thử nghiệm bổ sung enzyme phytase: Thiết lập nghiệm thức đối chứng (CT - Control Treatment: cho ăn thức ăn công nghiệp tiêu chuẩn) và nghiệm thức thí nghiệm (ET - Experiment Treatment: bổ sung enzyme phytase nguồn gốc Bacillus subtilis Ba 58 với tỷ lệ 0,01% khối lượng khô, tương đương 1500 UI/kg thức ăn).
• Phương pháp phân tích phòng thí nghiệm: Áp dụng Tiêu chuẩn Việt Nam (TCVN) và Standard Methods for the Examination of Water and Wastewater để xác định: pH, DO, nhiệt độ, EC, TSS, BOD₅, COD, TAN, NH₃, NH₄⁺, N-NO₂⁻, N-NO₃⁻, P-PO₄³⁻, TP, TN trong nước và hàm lượng hữu cơ, khoáng, kim loại nặng, TN, TP trong bùn đáy. Đánh giá các chỉ tiêu sinh trưởng: Tốc độ tăng trưởng ngày (DWG), Tỷ lệ sống, Hệ số chuyển đổi thức ăn (FCR), Hệ số sử dụng protein (PER), Tỷ lệ tích lũy phốt pho (PHR).

Phương pháp mô hình hóa toán học và thuật toán số

• Hệ phương trình vi phân: Xây dựng hệ phương trình vi phân thường phi tuyến mô tả tốc độ biến thiên theo thời gian $\frac{dC_i}{dt}$ của từng biến trạng thái (sinh khối cá, mật độ tảo, động vật phù du, nồng độ DO, NH₄⁺, NO₃⁻, PO₄³⁻, TN, TP nước và TN, TP bùn đáy).
• Công cụ giải số: Lập trình giải hệ phương trình vi phân bằng phương pháp Runge-Kutta bậc 4 trên phần mềm MATLAB.
• Phương pháp phân tích độ nhạy và hiệu chỉnh: Tiến hành phân tích độ nhạy cục bộ để sàng lọc các tham số có ảnh hưởng lớn nhất đến đầu ra mô hình; sử dụng thuật toán tối ưu hóa giảm thiểu sai số bình phương trung bình (MSE) và sai số căn bậc hai trung bình (RMSE) giữa giá trị mô phỏng và giá trị đo thực tế để hiệu chỉnh bộ tham số.

---

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

Chương 1: Tổng quan

Chương này tổng kết toàn diện các nghiên cứu về sinh học cá Tra (Pangasianodon hypophthalmus), kỹ thuật nuôi thâm canh và hiện trạng phát sinh chất thải tại ĐBSCL. Tác giả phân tích cơ chế dinh dưỡng, nhu cầu đạm (25–35%), lipid (5–11%) và vai trò của enzyme phytase trong việc giải phóng phytate-phosphorus. Đồng thời, chương 1 điểm luận chi tiết các cơ sở lý thuyết về mô hình hóa chất lượng nước ao hồ, phân tầng nhiệt, chu trình C-N-P-O₂, các mô hình hồ quốc tế (PCLake, STELLA, Ecopath) và khẳng định sự cần thiết phải xây dựng mô hình số chuyên biệt cho ao cá Tra Việt Nam.

Chương 2: Phương pháp nghiên cứu, thí nghiệm

Chương 2 mô tả chi tiết phương pháp tiếp cận kết hợp giữa thực nghiệm và mô hình hóa toán học:

• Quy trình bố trí thí nghiệm ao nuôi đối chứng và thí nghiệm bổ sung enzyme phytase (1500 UI/kg thức ăn).
• Danh mục phương pháp lấy mẫu, bảo quản và quy chuẩn phân tích các thông số thủy lý, thủy hóa, bùn đáy.
• Thiết lập mô hình khái niệm, ma trận tác động tương hỗ giữa các cấu tử sinh học và dinh dưỡng trong ao.
• Xây dựng hệ phương trình vi phân mô tả sinh trưởng cá, động lực học thực vật phù du, động vật phù du, cân bằng oxy, biến thiên các dạng đạm và lân trong cột nước và bùn đáy.
• Quy trình code hóa thuật toán Runge-Kutta bậc 4 trong môi trường MATLAB, cùng phương pháp phân tích độ nhạy tham số và kỹ thuật hiệu chỉnh số liệu.

Chương 3: Kết quả và thảo luận

Chương 3 trình bày các kết quả thực nghiệm và kết quả mô phỏng số, bao gồm:

1. Đánh giá ảnh hưởng của enzyme phytase: Việc bổ sung 0,01% phytase (1500 UI/kg) giúp cải thiện đáng kể các chỉ số sinh trưởng và chất lượng môi trường ao nuôi:
   • Trọng lượng cá thu hoạch đạt 41,10 ± 1,31 g/con (so với đối chứng 31,81 ± 1,14 g/con).
   • Tốc độ tăng trưởng đặc biệt đạt 1,68 ± 0,03 %/ngày (đối chứng đạt 1,39 ± 0,04 %/ngày).
   • FCR giảm từ 1,66 ± 0,09 (đối chứng) xuống 1,29 ± 0,06 (nghiệm thức phytase).
   • PER tăng từ 1,58 ± 0,08 lên 2,01 ± 0,09.
   • Tỷ lệ tích lũy phốt pho (PHR) tăng từ 47,83% lên 66,46%, giúp giảm phát thải phốt pho ra nước từ 21,82 ± 0,96 g P/kg cá tăng trọng xuống còn 10,39 ± 0,74 g P/kg cá (giảm 47,61% tải lượng P phát thải, $P < 0,05$).

Chỉ số theo dõi	Nghiệm thức đối chứng (CT)	Nghiệm thức bổ sung Phytase (ET - 1500 UI/kg)	Mức độ cải thiện / Biến động
Trọng lượng cá trung bình (g/con)	31,81 ± 1,14	41,10 ± 1,31	Tăng 29,2% ($P < 0,05$)
Tốc độ tăng trưởng (%/ngày)	1,39 ± 0,04	1,68 ± 0,03	Tăng 20,8% ($P < 0,05$)
Hệ số chuyển đổi thức ăn (FCR)	1,66 ± 0,09	1,29 ± 0,06	Giảm 22,3% ($P < 0,05$)
Hiệu quả sử dụng Protein (PER)	1,58 ± 0,08	2,01 ± 0,09	Tăng 27,2% ($P < 0,05$)
Tỷ lệ tích lũy Phốt pho (PHR, %)	47,83 ± 2,31	66,46 ± 2,38	Tăng 38,9% ($P < 0,05$)
Lượng P thải (g P/kg cá tăng trọng)	21,82 ± 0,96	10,39 ± 0,74	Giảm 47,61% ($P < 0,05$)

2. Hiện trạng môi trường ao nuôi thực địa:

   • Nhiệt độ nước ổn định 27,5–32 °C; độ sâu trung bình 4 m hạn chế phân tầng nhiệt.
   • Độ đục tăng dần và vượt 100 NTU từ giữa đến cuối vụ; TSS dao động 38–86 mg/L, cực đại đạt 274 mg/L.
   • DO biến động ngày-đêm sâu sắc: tầng mặt đạt 15,5 mg/L (200% bão hòa) vào xế chiều do tảo quang hợp mạnh và giảm xuống dưới 2,0 mg/L (thậm chí 0,44–0,58 mg/L) vào sáng sớm.
   • Nồng độ dinh dưỡng tăng vọt ở cuối vụ: TAN đạt 4,0–8,0 mg/L; NO₃⁻ dao động 0,12–19,5 mg/L; TN nước từ 4–5 mg/L đầu vụ tăng lên xấp xỉ 21–60 mg/L cuối vụ; TP nước từ 0,05 mg/L tăng lên 3,78–14,4 mg/L; P-PO₄³⁻ đạt 0,02–2,28 mg/L. Bùn đáy tích lũy hữu cơ 10,5–11,7%, TN chiếm 0,5% và TP chiếm 0,22%.

3. Kết quả mô phỏng và hiệu chỉnh mô hình số:

   • Quá trình phân tích độ nhạy đã xác định nhóm các tham số nhạy nhất (tốc độ quang hợp cực đại của tảo, hệ số hô hấp của cá, hệ số khoáng hóa hữu cơ, tốc độ trầm tích của phốt pho SEDP...) và nhóm tham số ít nhạy nhất.
   • Kết quả chạy mô hình trên MATLAB mô phỏng chính xác sự biến thiên của sinh khối cá, động thái bùng nổ và suy tàn của tảo lục/tảo lam ở 28 °C và 30 °C, sự suy giảm DO vào ban đêm, cũng như đường tích lũy TN, TP trong cột nước và bùn lắng qua 180 ngày nuôi, khớp với số liệu đo thực tế tại ao nuôi ĐBSCL.

---

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

Đóng góp mới về mặt khoa học và lý luận

1. Luận án đã làm rõ cơ chế biến động thủy lý, thủy hóa và động học tích lũy chất dinh dưỡng đặc trưng trong hệ sinh thái ao nuôi cá Tra thâm canh nước ngọt ở điều kiện nhiệt đới Việt Nam.
2. Xây dựng thành công mô hình toán học tích hợp mô tả các quá trình sinh thái - sinh hóa phức tạp trong ao cá Tra, liên kết đồng thời giữa chuỗi dinh dưỡng (thực vật phù du - động vật phù du - cá Tra) với các chu trình chuyển hóa C, N, P và cân bằng oxy hòa tan.
3. Ứng dụng thuật toán số Runge-Kutta bậc 4 trên ngôn ngữ MATLAB cùng quy trình phân tích độ nhạy tham số, cung cấp công cụ mô phỏng định lượng có độ tin cậy khoa học để dự báo chất lượng nước ao nuôi.

Đóng góp mới về mặt thực tiễn

1. Xác định tỷ lệ bổ sung enzyme phytase tối ưu: Đưa ra khuyến cáo khoa học về tỷ lệ bổ sung enzyme phytase vào thức ăn cá Tra là 0,01% theo trọng lượng khô (tương đương 1500 UI/kg thức ăn). Kết quả lượng hóa thực nghiệm chứng minh giải pháp này giúp giảm tới 47,61% lượng phốt pho phát thải ra môi trường nước, giảm FCR xuống 1,29 và tăng hiệu quả sử dụng protein lên 2,01.
2. Bộ công cụ hỗ trợ quản lý vận hành: Mô hình số cho phép người nuôi và nhà quản lý tính toán dự báo trước diễn biến suy giảm oxy, tích lũy amoni và tích tụ bùn lắng, từ đó tối ưu hóa chế độ cho ăn, điều chỉnh lịch thay nước (30–40%/ngày), thời điểm sục khí và chu kỳ nạo vét bùn đáy, giúp cắt giảm chi phí sản xuất và giảm tải lượng ô nhiễm xả ra các lưu vực sông tại ĐBSCL.

---

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

Hạn chế của luận án

• Mô hình toán học được thiết lập dựa trên giả thiết thủy vực ao nuôi khuấy trộn hoàn toàn (0 chiều/1 chiều đơn giản), chưa xét đến sự phân tầng vi mô cục bộ hoặc chuyển động dòng chảy phức tạp do cánh quạt nước/máy bơm phun tạo ra trong không gian 2D/3D.
• Các kịch bản mô phỏng chạy với nhiệt độ gán cố định (28 °C, 30 °C), chưa tích hợp chuỗi số liệu khí tượng biến thiên liên tục theo từng giờ trong ngày hoặc các đợt mưa lớn bất thường.
• Nghiên cứu chưa mô hình hóa chi tiết động lực học phát sinh và chuyển hóa của khí độc H₂S ở tầng đáy ao nuôi sâu.

Hướng nghiên cứu tiếp theo

• Nâng cấp mô hình toán lên không gian 2 chiều (2D) hoặc 3 chiều (3D) kết hợp động lực học chất lưu (CFD) để mô phỏng chính xác vùng lưu thông nước và lắng cặn trong ao.
• Tích hợp thêm module mô phỏng chu trình lưu huỳnh (khí độc H₂S) và tương tác với các biện pháp xử lý sinh học (chế phẩm vi sinh probiotic).
• Kết nối mô hình ao nuôi đơn lẻ với mô hình thủy lực - chất lượng nước lưu vực sông Tiền và sông Hậu để đánh giá sức tải môi trường tổng thể phục vụ quy hoạch vùng nuôi bền vững.

---

Giá trị tham khảo

Luận án là tài liệu chuyên khảo có giá trị ứng dụng và tham khảo cho nhiều nhóm đối tượng:

• Nghiên cứu sinh, học viên cao học và giảng viên: Chuyên ngành Kỹ thuật Môi trường, Khoa học Môi trường, Nuôi trồng Thủy sản và Sinh thái học thủy vực; có thể tham khảo trực tiếp cấu trúc hệ phương trình vi phân sinh thái, ma trận tác động tương hỗ và phương pháp lập trình giải số bằng MATLAB tại Chương 2 và phần phụ lục mã nguồn.
• Các nhà máy sản xuất thức ăn thủy sản và doanh nghiệp nuôi cá Tra: Tham khảo quy trình và tỷ lệ phối trộn enzyme phytase (0,01% khối lượng khô) để tối ưu hóa công thức thức ăn, hạ giá thành sản xuất và đạt các chứng nhận nuôi trồng thủy sản bền vững (VietGAP, GlobalGAP, ASC).
• Cơ quan quản lý nhà nước (Sở Tài nguyên và Môi trường, Chi cục Thủy sản các tỉnh ĐBSCL): Sử dụng các hệ số phát thải thực nghiệm và mô hình dự báo của luận án để tính toán tải lượng ô nhiễm N, P, xây dựng quy chuẩn kỹ thuật xả thải và quy hoạch vùng nuôi an toàn sinh thái.

---

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

1. Tỷ lệ bổ sung enzyme phytase được luận án xác định tối ưu là bao nhiêu và mang lại lợi ích cụ thể gì?

Tỷ lệ bổ sung enzyme phytase tối ưu được xác định là 0,01% theo trọng lượng khô (tương đương 1500 UI/kg thức ăn). Nghiệm thức này giúp:

• Tăng trọng lượng trung bình cá từ 31,81 g lên 41,10 g sau 2 tháng ương.
• Giảm hệ số chuyển đổi thức ăn (FCR) từ 1,66 xuống 1,29.
• Tăng hiệu quả sử dụng protein (PER) từ 1,58 lên 2,01.
• Tăng tỷ lệ tích lũy phốt pho (PHR) từ 47,83% lên 66,46%.
• Giảm lượng phốt pho thải ra môi trường nước 47,61% (từ 21,82 g P/kg cá tăng trọng xuống còn 10,39 g P/kg cá tăng trọng).

2. Nguyên nhân chính khiến chất lượng nước ao nuôi cá Tra thâm canh suy giảm nghiêm trọng ở cuối vụ nuôi là gì?

Vào cuối vụ nuôi, sinh khối cá trong ao đạt rất cao (năng suất 250–300 tấn/ha) kéo theo lượng thức ăn công nghiệp đưa vào ao mỗi ngày rất lớn. Cá chỉ hấp thu khoảng 37,5% vật chất dinh dưỡng, lượng còn lại tích lũy dưới đáy ao và hòa tan vào nước. Điều này dẫn đến sự gia tăng bùng nổ của tổng đạm amôn (TAN lên tới 4,0–8,0 mg/L), nitrat (NO₃⁻ lên tới 19,5 mg/L), tổng đạm (TN lên tới 21–60 mg/L) và tổng lân (TP lên tới 3,78–14,4 mg/L). Sự dư thừa dinh dưỡng kích thích tảo nở hoa với mật độ 1,6–10,5 triệu cá thể/L, làm độ đục vượt 100 NTU, gây dao động DO cực đoan (tăng 200% bão hòa ban ngày nhưng tụt sâu dưới 2 mg/L vào ban đêm) và làm gia tăng phân hủy yếm khí ở bùn đáy.

3. Mô hình toán của luận án được xây dựng và giải bằng công cụ nào?

Mô hình toán học của luận án bao gồm hệ phương trình vi phân thường phi tuyến biểu diễn các định luật bảo toàn khối lượng và tốc độ chuyển hóa sinh học giữa các thành phần (cá Tra, tảo, động vật phù du, DO, N-NH₄⁺, N-NO₃⁻, P-PO₄³⁻, TN, TP nước và bùn). Hệ phương trình này được giải số bằng thuật toán Runge-Kutta bậc 4 và được viết mã nguồn trên nền tảng ngôn ngữ lập trình MATLAB.

4. Quá trình phân tích độ nhạy trong mô hình nhằm mục đích gì?

Phân tích độ nhạy cục bộ nhằm đánh giá mức độ ảnh hưởng của từng tham số sinh hóa (tốc độ quang hợp, hệ số hô hấp, hằng số phân hủy, tốc độ lắng trầm tích...) đối với sự thay đổi của các biến trạng thái đầu ra (nồng độ DO, N, P, sinh khối sinh vật). Qua đó, tác giả sàng lọc được các tham số có độ nhạy lớn nhất để tập trung hiệu chỉnh bằng phương pháp tối ưu hóa sai số (MSE, RMSE) so với số liệu đo thực tế, đảm bảo mô hình đạt độ chính xác cần thiết khi vận hành mô phỏng.

---

Kết luận

Luận án của tác giả Lê Xuân Thịnh đã giải quyết trọn vẹn cả hai khía cạnh thực nghiệm môi trường và mô hình hóa toán học trong ao nuôi cá Tra thâm canh tại Đồng bằng sông Cửu Long. Nghiên cứu đã chứng minh hiệu quả của việc bổ sung 0,01% enzyme phytase giúp giảm 47,61% lượng phát thải phốt pho và nâng cao hiệu quả sử dụng thức ăn của cá. Đồng thời, mô hình số xây dựng trên nền tảng MATLAB bằng thuật toán Runge-Kutta bậc 4 đã mô phỏng thành công động học của các chất dinh dưỡng N, P, cân bằng DO và chuỗi sinh thái ao nuôi. Đây là cơ sở khoa học quan trọng phục vụ công tác quản lý chất thải, tối ưu hóa quy trình nuôi và thúc đẩy sự phát triển bền vững của ngành công nghiệp cá Tra Việt Nam.

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

Chương 1.1 Tổng quan về cá Tra 1.1 Phân loại Cá Tra là một trong số 11 loài thuộc họ cá tra (Pangasiidae) đã được xác định ở ĐBSCL. Hệ thống phân loại cá tra như sau: Giới: Animalia Nghành: Chordata Lớp: Actinopterygii Bộ: Siluriformes Họ: Pangasiidae Giống: Pangasius Loài: P. Cá Tra (Pangasianodon hypophthalmus) Tên loài Pangasianodon hypophthalmus được Rainboth sử dụng lần đầu vào năm 1996 [6] để chỉ định cho loài cá Tra và sau đó được nhiều tác giả khác sử dụng phổ biến đến nay. Hiện nay đã có 11 loài thuộc họ cá Tra được tìm thấy ở Việt Nam, trong đó có 5 loài là đối tượng nuôi quan trọng trong ao và bè.

Cá sống chủ yếu ở nước ngọt, chịu được nước lợ nhẹ (độ muối dưới 10%) chịu đựng được phèn có pH >4. Cá Tra có cơ quan hô hấp phụ nên có thể sống được ở những ao hồ chặt hẹp và phát triển bình thường trong điều kiện ao tù nước đọng, nhiều chất hữu cơ, thiếu oxy, nên nuôi được mật độ rất cao. Trong vòng đời của cá, giai đoạn cá bột hết noãn hoàng thì thích ăn mồi tươi sống, ăn các loài động vật phù du, cá bột còn ăn thịt lẫn nhau trong bể ương nuôi [7].2 Đặc điểm hình thái Theo [8] cá tra có những đặc điểm hình thái như sau: đầu rộng, dẹp bằng, mõm ngắn, nhìn từ trên xuống chót mõm tròn. Miệng trước rộng ngang, không co duỗi được có dạng hình vòng cung và nằm trên mặt phẳng ngang.

Răng nhỏ mịn, răng vòm miệng chia thành 4 đám nhỏ, mỏng, nằm trên đường vòng cung, đôi khi bị che lấp bởi nếp da vòm miệng. Lỗ mũi sau gần lỗ mũi trước hơn mắt và nằm trên đường thẳng kẻ từ lỗ mũi trước đến cạnh trên của mắt. Có hai đôi râu, râu mép kéo dài chưa chạm đến gốc vi ngực, râu cằm ngắn hơn. Thân thon dài, phần sau dẹp bên.

Đường bên hoàn toàn và Trang | 16 phân nhánh, bắt đầu từ mép trên của lỗ mang đến điểm giữa của gốc vi đuôi. Mặt sau của vi lưng, vi ngực có răng cưa hướng xuống gốc vi. Vi bụng kéo dài chưa chạm đến khởi điểm của gốc vi hậu môn.3 Đặc điểm phân bố Cá Tra phân bố ở lưu vực sông Mê kông, có mặt ở cả 4 nước Lào, Việt Nam, Cămpuchia và Thái lan. Ở Thái Lan còn gặp cá tra ở lưu vực sông Mekong và Chao Phraya.

Theo đặc tính thích ứng môi trường, cá Tra không phân bố thuộc các tỉnh phía Bắc, cá Tra thuộc đối tượng cá nuôi vùng nhiệt đới. Ở nước ta những năm trước đây khi chưa có cá sinh sản nhân tạo, cá bột và cá giống Tra được vớt trên sông Tiền và sông Hậu. Cá trưởng thành chỉ thấy trong ao nuôi, rất ít gặp trong tự nhiên địa phận Việt nam, do cá có tập tính di cư ngược dòng sông Mekong để sinh sống và tìm nơi sinh sản tự nhiên. Khảo sát chu kỳ di cư của cá tra ở địa phận Campuchia cho thấy cá ngược dòng từ tháng 10 đến tháng 5 và di cư về hạ lưu từ tháng 5 đến tháng 9 hàng năm [8].4 Đặc điểm sinh trưởng Cá tra có tốc độ tăng trưởng tương đối nhanh, còn nhỏ cá tăng nhanh về chiều dài.

Cá ương trong ao sau 2 tháng đã đạt chiều dài 10-12cm (14-15gam). Từ khoảng 2,5 kg trở đi, mức tăng trọng lượng nhanh hơn so với mức tăng chiều dài cơ thể. Cỡ cá trên 10 tuổi trong tự nhiên (ở Campuchia) tăng trọng rất ít. Cá tra trong tự nhiên có thể sống trên 20 năm.

Trong tự nhiên có cá đạt 18 kg hoặc có mẫu cá dài tới 1,8m. Trong ao nuôi vỗ, cá bố mẹ cho đẻ đạt tới 25 kg ở cá 10 năm tuổi.