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 về luận án

Nghề nuôi cá Tra (Pangasianodon hypophthalmus) tại vùng Đồng bằng sông Cửu Long (ĐBSCL) đóng vai trò xương sống trong cơ cấu kinh tế thủy sản Việt Nam, đóng góp hơn 1,61 tỷ USD kim ngạch xuất khẩu tới 131 thị trường quốc tế và giải quyết việc làm cho hàng triệu lao động. Tuy nhiên, sự phát triển nóng với mật độ thâm canh cực cao (năng suất trung bình 250 – 300 tấn/ha/vụ) đã tạo áp lực sinh thái nghiêm trọng lên các thủy vực tiếp nhận. Luận án tiến sĩ kỹ thuật môi trường "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" của tác giả Lê Xuân Thịnh (Trường Đại học Bách Khoa Hà Nội, hướng dẫn bởi PGS.TS Đặng Xuân Hiển và PGS.TS Trần Văn Nhân) là công trình tiên phong giải quyết bài toán cân bằng sinh thái - kỹ thuật thông qua công cụ mô hình hóa số học kết hợp giải pháp sinh học phân tử trong thức ăn.

Nghiên cứu xuất phát từ khoảng trống tri thức (research gap) rõ rệt: hầu hết các nghiên cứu trước đây (như nghiên cứu của Boyd, Tucker & Hargreaves trên cá nheo Mỹ hay Janse trên mô hình PCLake) chủ yếu tập trung vào hồ tự nhiên hoặc ao nuôi nước nông vùng ôn đới, thiếu vắng các mô hình cơ chế (mechanistic models) mô phỏng động học tương hỗ giữa cột nước sâu (3 – 4 m) và trầm tích đáy trong điều kiện nhiệt đới gió mùa có sinh khối nuôi siêu thâm canh. Các câu hỏi và giả thuyết nghiên cứu được xác lập chặt chẽ:

• Câu hỏi nghiên cứu 1 (RQ1): Động lực biến thiên không - thời gian của các thông số thủy lý, thủy hóa (DO, TAN, $\text{NO}_2^-$, $\text{NO}_3^-$, $\text{PO}_4^{3-}$, TN, TP) và sự phát triển của các quần xã thủy sinh (tảo lục, tảo lam, tảo cát, giáp xác râu ngành, chân chèo) trong ao nuôi diễn ra theo quy luật nào?
• Câu hỏi nghiên cứu 2 (RQ2): Việc bổ sung enzyme phytase vào khẩu phần thức ăn công nghiệp ảnh hưởng như thế nào đến khả năng tiêu hóa, hiệu quả sử dụng thức ăn và tỷ lệ phát thải phốt pho ra môi trường ao?
• Câu hỏi nghiên cứu 3 (RQ3): Làm thế nào để thiết lập một mô hình toán học giải tích - số trị mô phỏng chính xác chu trình chuyển hóa C-N-P-O và sinh khối cá nhằm phục vụ tối ưu hóa vận hành?
• Giả thuyết nghiên cứu (H1, H2, H3): Bổ sung phytase với liều lượng tối ưu 0,01% trọng lượng khô (tương đương 1500 UI/kg) sẽ tăng tỷ lệ tích lũy phốt pho (PHR), giảm phát thải phốt pho trên 40%; hệ phương trình vi phân phi tuyến bậc một tích hợp động học sinh hóa có khả năng tái hiện diễn biến chất lượng nước với sai số bình phương trung bình gốc (RMSE) chấp nhận được.

Khung lý thuyết của luận án được định hình từ lý thuyết cân bằng vật chất sinh thái (Ecological Stoichiometry), động học sinh hóa vi sinh vật Monod/Michaelis-Menten và lý thuyết dòng năng lượng trong mạng lưới thức ăn (Trophic Food Web Dynamics). Về quy mô, nghiên cứu bao phủ các ao nuôi thực địa tại các vùng trọng điểm Cần Thơ, An Giang, Đồng Tháp trong chu kỳ nuôi 6 – 8 tháng kết hợp các thí nghiệm đối chứng nghiêm ngặt trong hệ thống bể composite.

Literature Review và Positioning

Tổng quan tài liệu quốc tế và trong nước cho thấy sự phân hóa rõ rệt giữa hai luồng quan điểm trong quản lý chất lượng nước ao nuôi thâm canh. Luồng quan điểm truyền thống (Boyd, 1990; Gross et al., 2000) nhấn mạnh vào việc kiểm soát môi trường thông qua các giải pháp ngoại vi như tăng cường tốc độ trao đổi nước (water exchange rate từ 30 – 50%/ngày) và sử dụng hóa chất xử lý đáy ao. Ngược lại, luồng quan điểm sinh thái học hiện đại (Hargreaves, 1998; Janse, 2005; Verdegem, 2013) khẳng định rằng việc thay nước chỉ đơn thuần chuyển dịch tải lượng ô nhiễm ra sông ngòi tự nhiên, đồng thời đề xuất giải pháp cốt lõi phải xuất phát từ việc kiểm soát dòng dinh dưỡng nội tại (nutrient flux) và tối ưu hóa hệ số chuyển hóa thức ăn thông qua dinh dưỡng học chính xác.

Tranh luận học thuật lớn thứ hai liên quan đến hiệu lực sinh học của enzyme phytase trong đường tiêu hóa của động vật thủy sản. Trong khi một số nghiên cứu trên cá hồi vân (Oncorhynchus mykiss) của Hughes & Soares (1998) hay cá da trơn Mỹ của Robinson et al. (2002) cho rằng phytase không mang lại hiệu quả rõ rệt khi khẩu phần có tỷ lệ đạm thực vật cao do thời gian lưu giữ thức ăn trong ruột ngắn và pH đường tiêu hóa biến động, các công trình của Baruah et al. (2007) và Trần Thị Thanh Hiền et al. (2010) lại chỉ ra rằng ở các loài cá da trơn nhiệt đới, phytase thủy phân hiệu quả liên kết myo-inositol hexakisphosphate, giải phóng phốt pho vô cơ và khoáng vi lượng.

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

1. Mô hình lưới thức ăn hồ Ontario của Halfon (1996): Sử dụng hệ phương trình cân bằng dòng năng lượng sinh thái $Tiêu\ thụ = Hô\ hấp + SDA + Sản\ xuất + Bài\ tiết$ cho thủy vực hở lớn nhưng chưa tích hợp cơ chế tương tác đa tầng giữa bùn lắng đáy ao và cột nước mặt.
2. Mô hình phú dưỡng hóa PCLake (Janse et al., 2005; 2010): Mô hình hóa chi tiết chu trình photpho và hiện tượng chuyển pha sinh thái trong các hồ nông châu Âu, tuy nhiên hệ số tải lượng dinh dưỡng đầu vào thấp hơn hàng chục lần so với mật độ nuôi cá Tra thâm canh tại ĐBSCL.

Luận án định vị chính xác khoảng trống học thuật này bằng cách kết hợp cơ chế dinh dưỡng vi mô (bổ sung phytase) với mô hình động học sinh thái vĩ mô (mô hình hóa toán học 15 biến trạng thái) dành riêng cho ao nuôi nhiệt đới siêu thâm canh.

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

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

Luận án mở rộng lý thuyết Động thái dinh dưỡng sinh thái (Ecological Nutrient Dynamics) và lý thuyết Cân bằng vật chất (Mass Balance Theory) khi áp dụng vào hệ sinh thái nhân tạo có áp lực trao đổi chất cực đoan. Công trình chứng minh rằng trong ao nuôi cá Tra, bùn đáy đóng vai trò là bể chứa (sink) lẫn nguồn phát tán ngược (source) phốt pho và nitơ tùy thuộc vào gradient nồng độ oxy hòa tan (DO) và hoạt động đảo trộn tầng nước của cá.

Hệ thống luận điểm/đề xuất lý thuyết được lượng hóa:

• Proposition 1: Động học phân hủy phytate bởi enzyme ngoại sinh trong đường ruột cá Tra tỷ lệ thuận với khả năng ức chế quá trình hình thành phức hợp phytate-khoáng không tan, trực tiếp làm giảm tải lượng phốt pho bài tiết qua phân.
• Proposition 2: Tốc độ khoáng hóa hữu cơ ở lớp bùn đáy chịu sự chi phối phi tuyến của tỷ lệ C:N:P trong thức ăn thừa và sinh khối tảo tàn, trong đó tốc độ giải phóng orthophosphate ($\text{PO}_4^{3-}$) vào nước bị khống chế bởi khả năng hấp phụ trao đổi cation của bùn sét phù sa ĐBSCL.

graph TD
    A["Thức ăn bổ sung 0,01% Phytase"] --> B["Cá Tra (Pangasianodon hypophthalmus)"]
    B -->|"Tích lũy sinh khối (DM: 32,6%; N: 42,6%; P: 29,3%)"| C["Thu hoạch thương phẩm"]
    B -->|"Bài tiết hòa tan (TAN, NH3, PO4 3-)"| D["Cột nước (Water Column)"]
    B -->|"Phân thải & Thức ăn thừa"| E["Bùn đáy ao (Sediment)"]
    D <-->|"Quang hợp / Hô hấp / Tảo tàn"| F["Thực vật phù du (Tảo lục, lam, cát)"]
    F <-->|"Chuỗi thức ăn"| G["Động vật phù du (Copepods, Cladocerans)"]
    E <-->|"Khuếch tán ngược / Hấp phụ / Phản nitrat hóa"| D

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

Khung phân tích của luận án tích hợp ba trụ cột lý thuyết: Lý thuyết cân bằng hóa học lượng sinh thái, Động học enzyme vi sinh vật và Thủy sinh học hệ thống. Ma trận tác động tương hỗ được thiết lập giữa 15 biến trạng thái (Cá, TVPD, DVPD, DO, BOD, COD, TSS, TAN, $\text{NH}_3$, $\text{NH}_4^+$, $\text{NO}_2^-$, $\text{NO}_3^-$, $\text{PO}_4^{3-}$, TN, TP) và các tham số môi trường.

Điều kiện biên (boundary conditions) được xác định rõ: Ao nuôi nước ngọt nhiệt đới có độ sâu 3 – 4 m, nhiệt độ nước dao động 27,5 – 32 °C, độ mặn < 10‰, pH 6,3 – 9,2, chế độ thay nước định kỳ từ 30 – 40%/ngày sau 2 tháng nuôi đầu tiên.

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

Thiết kế nghiên cứu

Luận án theo đuổi thế giới quan thực chứng (Positivism) với thiết kế thực nghiệm định lượng kết hợp mô hình hóa giải tích số trị đa cấp độ (Multi-level mechanistic design). Khung nghiên cứu bao gồm: Cấp độ vi mô (thí nghiệm dinh dưỡng bể composite kiểm soát các biến sinh hóa), Cấp độ trung mô (khảo sát ao nuôi thực địa đa điểm tại Cần Thơ, An Giang, Đồng Tháp) và Cấp độ vĩ mô (xây dựng phần mềm mô phỏng bằng ngôn ngữ MATLAB).

Cỡ mẫu và tiêu chí lựa chọn:

• Thí nghiệm bể: Nghiệm thức đối chứng (CT) và nghiệm thức thí nghiệm (ET) bổ sung phytase từ Bacillus subtilis Ba 58 ở các mức 500, 1000, 1500 UI/kg thức ăn (mỗi nghiệm thức lặp lại 3 lần, n = 30 cá/bể).
• Khảo sát thực địa: Đánh giá tại các ao nuôi công nghiệp mật độ 48 – 70 con/m², độ sâu trung bình 4 m, diện tích $1.000\ \text{m}^2 - 5.000\ \text{m}^2$ qua các chu kỳ nuôi 6 – 8 tháng.

flowchart LR
    subgraph ThuThap [Thu thập dữ liệu]
        A1[Thí nghiệm bổ sung Phytase] --> A3[Phân tích mẫu nước & bùn đáy]
        A2[Khảo sát ao thực địa ĐBSCL] --> A3
    end
    subgraph MoHinh [Xây dựng mô hình]
        A3 --> B1[Thiết lập ma trận tương hỗ 15 biến]
        B1 --> B2[Hệ phương trình vi phân phi tuyến ODE]
        B2 --> B3[Giải thuật Runge-Kutta bậc 4 trên MATLAB]
    end
    subgraph KiemChung [Hiệu chỉnh & Đánh giá]
        B3 --> C1[Phân tích độ nhạy tham số]
        C1 --> C2[Hiệu chỉnh mô hình qua RMSE/MSE]
        C2 --> C3[Dự báo & Đề xuất quản lý vận hành]
    end

Quy trình nghiên cứu rigorous

Quy trình thu mẫu và phân tích tuân thủ nghiêm ngặt Tiêu chuẩn Việt Nam (TCVN) và Standard Methods for the Examination of Water and Wastewater (APHA):

• Thông số thủy lý: Đo tầng mặt và tầng đáy bằng máy đo đa chỉ tiêu (Horiba U-50, Hanna HI98194) cho nhiệt độ, pH, DO, EC, TDS, độ đục.
• Thông số thủy hóa: Xác định BOD theo phương pháp ủ 5 ngày ở 20 °C; COD theo phương pháp chuẩn độ Bicromat; TSS bằng phương pháp sấy ở 105 °C; TAN, $\text{NO}_2^-$, $\text{NO}_3^-$, $\text{PO}_4^{3-}$ bằng phương pháp so màu quang phổ kế UV-Vis; TN và TP theo phương pháp phân hủy mẫu bằng persulfate và nung Kjeldahl.
• Trầm tích bùn đáy: Thu bằng gàu cuốc Ekman, phân tích hàm lượng chất hữu cơ (mất khi nung LOI ở 550 °C), TN, TP, Ca, Mg và kim loại nặng (Cu, Zn, Fe, Pb) bằng quang phổ hấp thụ nguyên tử (AAS).
• Tam giác đạc (Triangulation): Kết hợp kiểm chứng chéo giữa dữ liệu hóa sinh thực nghiệm, quan trắc chu kỳ ngày - đêm tại ao thực tế và mô phỏng số học.

Data và phân tích

Hệ phương trình vi phân thường phi tuyến (Ordinary Differential Equations - ODEs) mô tả tốc độ biến thiên sinh khối và nồng độ chất dinh dưỡng được xây dựng dựa trên nguyên lý bảo toàn khối lượng:

$$\frac{dC_i}{dt} = \sum R_{sinh_trưởng} - \sum R_{tử_vong} - \sum R_{bài_tiết} \pm \sum R_{chuyển_hóa} \pm F_{trao_đổi}$$

Trong đó, các phương trình thành phần bao gồm: phương trình tăng trưởng cá theo mô hình năng lượng sinh học von Bertalanffy mở rộng; phương trình quang hợp và tử vong của thực vật phù du theo hàm ánh sáng Steele và động học dinh dưỡng Michaelis-Menten; phương trình tiêu thụ oxy tổng (SOD); phương trình động học chu trình Nitơ ($\text{NH}_4^+ \rightarrow \text{NO}_2^- \rightarrow \text{NO}3^-$) và chu trình Phốt pho ($\text{P}{hữu_cơ} \rightarrow \text{PO}_4^{3-}$).

Hệ phương trình được lập trình giải số bằng thuật toán Runge-Kutta bậc 4 (RK4) trong môi trường MATLAB. Quá trình phân tích độ nhạy (Sensitivity Analysis) được thực hiện để xác định các tham số có tác động chi phối lớn nhất đến hàm mục tiêu, tiếp nối bởi bước hiệu chỉnh mô hình (Calibration) nhằm tối thiểu hóa 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 đạc thực tế.

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

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

[!IMPORTANT] Dữ liệu cân bằng vật chất thực nghiệm chỉ ra rằng: "sản xuất 1 tấn cá Tra thành phẩm là 4.023m3 nước... với ao nuôi đạt năng suất 300 tấn/ha/vụ thì mỗi vụ nuôi sẽ thải ra môi trường khoảng 2.677 tấn bùn ướt (tương đương 937 tấn bùn khô)... khi cho cá ăn 100% lượng thức ăn công nghiệp thì lượng vật chất dinh dưỡng mà cá tích lũy được trong cơ thịt chỉ chiếm khoảng 37,5%".

Công trình đã khám phá các quy luật then chốt sau:

1. Hiệu quả đột phá của bổ sung Phytase: Ở nghiệm thức bổ sung 1500 UI enzyme phytase/kg thức ăn (tỷ lệ 0,01% trọng lượng khô), tốc độ tăng trưởng đặc biệt của cá đạt $1,68 \pm 0,03\ %/\text{ngày}$ (vượt trội so với đối chứng $1,39 \pm 0,04\ %/\text{ngày}$, $p < 0,05$). Trọng lượng trung bình sau 2 tháng đạt $41,10 \pm 1,31\ \text{g/con}$ (đối chứng chỉ đạt $31,81 \pm 1,14\ \text{g/con}$). Hệ số FCR giảm mạnh xuống mức $1,29 \pm 0,06$ so với $1,66 \pm 0,09$ ở nhóm đối chứng; hiệu quả sử dụng protein (PER) tăng lên $2,01 \pm 0,09$ so với $1,58 \pm 0,08$.

2. Lượng hóa mức độ giảm phát thải Phốt pho: Khẩu phần bổ sung phytase 1500 UI/kg giúp nâng tỷ lệ tích lũy phốt pho (PHR) từ $47,83 \pm 2,31%$ lên $66,46 \pm 2,38%$. Lượng phốt pho thải ra môi trường nước giảm từ $21,82 \pm 0,96\ \text{g P/kg}$ tăng trọng xuống còn $10,39 \pm 0,74\ \text{g P/kg}$ tăng trọng, tương đương mức cắt giảm phát thải thực tế lên đến 47,61% ($p < 0,05$).

3. Hiện tượng tích lũy dinh dưỡng cực đoan tại đáy ao: Trong ao nuôi thương phẩm không bổ sung enzyme, lượng chất thải dạng rắn chiếm đến 45,6% vật chất khô, 50,4% tổng nitơ và 65,57% tổng phốt pho đưa vào qua thức ăn. Hàm lượng hữu cơ trong bùn đáy tích lũy từ 10,5 – 11,7%, TN chiếm 0,5% và TP chiếm 0,22%.

4. Động thái suy giảm chất lượng nước cuối vụ: Vào giai đoạn cuối chu kỳ nuôi (tháng thứ 6 – 8), nồng độ TAN đạt đỉnh 6,83 – 8,0 mg/L (vượt gấp 5 – 10 lần ao nuôi thủy sản thông thường), $\text{NO}_3^-$ đạt 18,0 – 19,5 mg/L (cao gấp 6 – 7 lần ao nuôi cá nheo Mỹ), chlorophyll-a vượt ngưỡng $600\ \mu\text{g/L}$ gây hiện tượng tảo nở hoa dày đặc và sụp đổ oxy hòa tan xuống dưới 2 mg/L vào rạng sáng.

Thông số môi trường	Ao cá Tra ĐBSCL (Cuối vụ)	Ao cá nheo Mỹ (I. punctatus)	Ngưỡng an toàn / QCVN
Tổng đạm amôn (TAN)	$4,0 - 8,0\ \text{mg/L}$	$0,4 - 2,7\ \text{mg/L}$	$< 1,0\ \text{mg/L}$ (QCVN 08)
Tổng Nitơ (TN)	$21,0 - 60,0\ \text{mg/L}$	$4,3 - 10,6\ \text{mg/L}$	$< 20\ \text{mg/L}$ (Cột A, QCVN 40)
Tổng Phốt pho (TP)	$1,6 - 3,78\ \text{mg/L}$	$0,2 - 0,8\ \text{mg/L}$	$< 0,2\ \text{mg/L}$ (Khuyến nghị ao)
Hàm lượng Chlorophyll-a	$200 - 600+\ \mu\text{g/L}$	$39 - 128\ \mu\text{g/L}$	$50 - 200\ \mu\text{g/L}$ (Mức sinh thái)
Lượng bùn khô phát sinh	$937\ \text{tấn/ha/vụ}$	$50 - 120\ \text{tấn/ha/vụ}$	Cần thu gom định kỳ

5. Độ chính xác và tính tương thích của mô hình số: Kết quả mô phỏng trên MATLAB sau khi hiệu chỉnh tham số cho thấy quỹ đạo tính toán của DO, TAN, $\text{NO}_3^-$, $\text{PO}_4^{3-}$, TN, TP và sinh khối cá bám sát diễn biến thực tế với hệ số tương quan $R^2 > 0,85$, chứng minh tính khả thi của việc sử dụng mô hình để dự báo chất lượng nước ao.

Implications đa chiều

• Ý nghĩa học thuật: Bổ sung cơ sở dữ liệu số hóa và bộ tham số động học sinh thái học nhiệt đới cho ngành kỹ thuật môi trường nước mặt và mô hình hóa sinh thái thế giới.
• Đổi mới phương pháp luận: Thiết lập quy trình chuẩn kết hợp giữa công nghệ enzyme phân tử thức ăn và giải pháp mô phỏng số học bậc cao trong kiểm soát tải lượng nguồn thải nuôi trồng thủy sản.
• Ứng dụng thực tiễn: Cho phép các trang trại cá Tra chủ động thiết lập lịch trình cho ăn, tối ưu hóa thời điểm xi-phông hút bùn đáy, điều chỉnh lưu lượng bơm quạt nước cung cấp oxy và xác định chính xác chu kỳ thay nước (30 – 40%) để tránh sốc môi trường và giảm chi phí năng lượng.
• Chính sách môi trường: Cung cấp cơ sở khoa học định lượng phục vụ việc rà soát, hoàn thiện Quy chuẩn kỹ thuật quốc gia về nước thải nuôi trồng thủy sản (QCVN 02-20:2014/BNNPTNT và QCVN 40:2011/BTNMT), hỗ trợ thực thi Quyết định số 339/QĐ-TTg của Thủ tướng Chính phủ về Chiến lược phát triển Thủy sản Việt Nam bền vững.

Limitations và Future Research

Nghiên cứu thể hiện tính trung thực học thuật cao khi chỉ rõ các giới hạn nội tại:

1. Giả định xáo trộn hoàn toàn (Completely Mixed Tank Assumption): Mô hình số học hiện tại xem ao nuôi như một bể phản ứng đồng nhất 1D theo phương thẳng đứng, chưa tích hợp hoàn chỉnh mô hình thủy động lực học 3D (3D hydrodynamics) để mô phỏng chi tiết các góc chết dòng chảy do vị trí đặt máy bơm tạo dòng gây ra.
2. Quy mô quần thể sinh vật phù du: Các nhóm tảo và động vật phù du được gộp thành các nhóm chức năng lớn (tảo lục, tảo lam, tảo cát; copepods, cladocerans) mà chưa mô hình hóa chi tiết động thái cạnh tranh tới từng loài đặc hữu.
3. Biến số khí hậu cực đoan: Các thông số môi trường đầu vào của mô hình chưa bao gồm kịch bản xâm nhập mặn sâu và biến đổi nhiệt độ bất thường do biến đổi khí hậu tại hạ lưu sông Mekong.

Chương trình nghiên cứu tương lai cần hướng tới:

• Phát triển mô hình tích hợp 3 chiều Thủy động lực học - Hóa sinh - Sinh thái (Coupled 3D Hydrodynamic-Biogeochemical Models) sử dụng nền tảng mã nguồn mở (như Delft3D hoặc EFDC).
• Nghiên cứu cơ chế tác động tương tác giữa enzyme phytase và các phức hợp enzyme ngoại sinh khác (protease, xylanase) lên hệ vi sinh vật đường ruột cá Tra thông qua kỹ thuật giải trình tự gen Metagenomics.
• Ứng dụng trí tuệ nhân tạo (AI) và mạng thần kinh nhân tạo (ANN) để tự động hóa quá trình đồng hóa dữ liệu quan trắc từ các cảm biến IoT thời gian thực vào mô hình số nhằm cảnh báo sớm hiện tượng thiếu hụt oxy và bùng phát khí độc trong ao nuôi.

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

Luận án mở ra tiềm năng trích dẫn học thuật cao trên các tạp chí quốc tế chuyên ngành kỹ thuật môi trường và thủy sản (như Aquaculture, Ecological Modelling, Water Research). Về phương diện kinh tế - kỹ thuật, việc áp dụng bổ sung enzyme phytase (0,01%) và vận hành ao nuôi theo chỉ dẫn của mô hình số ước tính giúp toàn ngành cá Tra ĐBSCL:

• Tiết kiệm 15 – 20% chi phí thức ăn nhờ giảm hệ số FCR từ 1,56 – 1,66 xuống 1,29 – 1,35.
• Cắt giảm trực tiếp hàng chục ngàn tấn phốt pho thải ra hệ thống sông Tiền và sông Hậu mỗi năm.
• Giảm chi phí điện năng vận hành máy bơm và hóa chất xử lý đáy ao từ 20 – 30%.
• Thúc đẩy ngành nuôi cá Tra Việt Nam đạt chứng nhận quốc tế bền vững (ASC, GlobalGAP, BAP), nâng cao năng lực cạnh tranh xuất khẩu tại các thị trường khó tính như EU, Mỹ và Nhật Bản.

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

• Nghiên cứu sinh và Giảng viên đại học: Tiếp cận bộ khung toán học hoàn chỉnh, phương pháp giải số RK4 trên MATLAB và bộ cơ sở dữ liệu thực nghiệm chi tiết về động học dinh dưỡng nhiệt đới.
• Nhà máy sản xuất thức ăn thủy sản: Cơ sở khoa học vững chắc để chuẩn hóa công thức phối trộn thức ăn công nghiệp bổ sung enzyme phytase tối ưu hàm lượng đạm/năng lượng và giảm thiểu phối trộn phốt pho vô cơ đắt đỏ.
• Chủ trang trại và Kỹ sư vận hành ao nuôi: Công cụ mô phỏng dự báo hỗ trợ ra quyết định hàng ngày về lượng thức ăn, thời gian sục khí, lịch xả thải bùn đáy và thay nước.
• Nhà hoạch định chính sách và Quản lý môi trường: Căn cứ định lượng để quy hoạch phân vùng nuôi, cấp phép xả thải và kiểm soát tải lượng môi trường lưu vực sông ĐBSCL.

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

Đóng góp độc đáo nhất là việc mở rộng Lý thuyết cân bằng hóa học lượng sinh thái (Ecological Stoichiometry) vào hệ sinh thái thủy vực nhân tạo siêu thâm canh nhiệt đới. Luận án đã tích hợp đồng thời cơ chế thủy phân phytate ở cấp độ phân tử bên trong ruột cá với quá trình biến đổi pha và chuyển hóa C-N-P-O đa cấp trong cột nước và lớp bùn đáy dày đặc chất hữu cơ.

2. Phương pháp nghiên cứu đổi mới như thế nào so với các nghiên cứu trước?

So với mô hình PCLake của Janse et al. (chủ yếu cho hồ tự nhiên châu Âu) hay nghiên cứu cân bằng dinh dưỡng tĩnh của Tucker & Hargreaves (chỉ đánh giá nồng độ độc lập trên cá nheo), luận án đã:

• Thiết lập hệ thống ma trận tác động tương hỗ 15 biến trạng thái tương thích với điều kiện ao nước sâu (3 – 4 m) vùng nhiệt đới.
• Ứng dụng giải thuật số Runge-Kutta bậc 4 trên ngôn ngữ máy tính MATLAB kết hợp phân tích độ nhạy cục bộ để giải hệ phương trình vi phân phi tuyến có độ phức tạp cao, cho phép dự báo liên tục biến thiên chất lượng nước theo thời gian thự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à lượng phốt pho tích lũy trong cá ở ao nuôi đối chứng chỉ đạt 29,3 – 29,8%, trong khi lượng phốt pho lắng đọng trong bùn đáy chiếm tới 65,57% tổng lượng phốt pho đưa vào từ thức ăn. Tuy nhiên, chỉ cần bổ sung một tỷ lệ cực nhỏ enzyme phytase (0,01% trọng lượng khô hay 1500 UI/kg thức ăn), lượng phốt pho bài tiết vào nước đã giảm tới 47,61% mà không làm thay đổi tiêu cực đến sinh lý của cá.

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

Có. Luận án mô tả chi tiết quy trình bố trí thí nghiệm, công thức thức ăn (22% và 28% protein), phương pháp thu mẫu nước và trầm tích chuẩn APHA/TCVN, mã nguồn code số MATLAB và bảng giá trị tham số sinh hóa (tốc độ tăng trưởng cực đại, hệ số nửa bão hòa, tốc độ phân hủy hữu cơ, hệ số hô hấp) cho phép các nhà nghiên cứu độc lập tái lập hoàn toàn mô hình.

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

Lộ trình 10 năm tập trung vào 3 mũi nhọn:

• Nâng cấp mô hình 1D lên không gian 3D tương tác thủy động lực học thời gian thực.
• Xây dựng hệ thống bản sao số (Digital Twin) tích hợp mạng lưới cảm biến IoT trên các vùng nuôi tập trung ĐBSCL.
• Khảo sát sâu cơ chế tương tác đa enzyme và công nghệ sinh học xử lý bùn đáy ao nuôi cá Tra thành phân bón hữu cơ cao cấp phục vụ kinh tế tuần hoàn.

Kết luận

Luận án tiến sĩ của tác giả Lê Xuân Thịnh đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 5 đóng góp cốt lõi:

1. Xác lập bộ cơ sở dữ liệu toàn diện về hiện trạng thủy lý, thủy hóa, sinh học và động lực tích lũy bùn đáy trong ao nuôi cá Tra thâm canh tại ĐBSCL.
2. Xác định tỷ lệ bổ sung enzyme phytase tối ưu (0,01% trọng lượng khô / 1500 UI/kg thức ăn), chứng minh hiệu quả giảm 47,61% lượng phốt pho phát thải ra môi trường, giảm FCR xuống 1,29 và nâng cao hiệu quả sử dụng protein (PER = 2,01).
3. Phát triển thành công mô hình toán học giải tích - số trị mô phỏng chu trình chuyển hóa C-N-P-O và sinh khối thủy sinh vật chuyên biệt cho ao nuôi cá Tra Việt Nam.
4. Lập trình và hiệu chỉnh hoàn chỉnh bộ công cụ tính toán trên nền tảng MATLAB với độ tin cậy và hệ số tương quan cao, hỗ trợ đắc lực cho công tác quản lý vận hành ao nuôi.
5. Đóng góp luận cứ khoa học quan trọng cho việc xây dựng tiêu chuẩn kỹ thuật môi trường và chính sách phát triển bền vững ngành hàng cá Tra quốc gia trong kỷ nguyên thích ứng biến đổi khí hậu.

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. Nuôi trong ao một năm cá đạt 1-1,5kg/con ở năm đầu tiên, những năm về sau cá tăng trọng nhanh hơn, có khi đạt tới 5-6kg/năm tùy thuộc môi trường sống và sự cung cấp thức ăn cũng như loại thức ăn có hàm lượng đạm nhiều hay ít. Độ béo fulton của cá tăng dần theo trọng lượng và nhanh nhất ở những năm đầu, cá đực thường có độ béo cao hơn cá cái và độ béo thường giảm đi vào mùa sinh sản [13].

Đối với cá Tra ở giai đoạn bột lên hương tốc độ tăng trưởng rất nhanh, sau 20-25 ngày ương thì cá đạt kích cỡ khoảng 2,7-3,0 cm theo quy định kích cỡ cá hương (Tiêu chuẩn Việt Nam 28 TCVN 211:2004).5 Dinh dưỡng và thức ăn cá Tra 1.1 Cân bằng dinh dưỡng trong ao nuôi Sự phân bố nguồn vật chất và cân bằng vật chất dinh dưỡng trong ương nuôi cá Tra đã được công bố qua một số công trình nghiên cứu. Nghiên cứu [9] trên cá giống có chiều cao thân từ 0,5-1,5 cm được nuôi trong hệ thống tuần hoàn và cá được cho ăn theo nhu cầu với thức ăn viên công nghiệp loại viên nổi 40% chất đạm. Kết quả cho thấy có 13% N và 8,45% vật chất khô (DM) từ thức ăn không được tiêu hóa và bị thải qua phân cá. Trong phần thức ăn được tiêu hóa có 35,3% N và 34,1% DM được tích lũy trong cơ thể cá thông qua quá trình sinh trưởng; 42,5% N và 29,1% DM được sử dụng cho các quá trình trao đổi chất và chúng được thải ra môi trường nước dưới dạng chất hòa tan (CO2, NH3, H2O,.) qua mang hay thận.

Vì vậy, sản xuất 1 kg cá giống thì cần cung cấp cho cá 765 g DM và 48,7 g N; trong đó cá tích lũy được 252 g DM và 18,5 g N; 514 g DM và 30,2 g N được thải ra môi trường dưới dạng chất vô cơ hòa tan và chất hữu cơ lơ lửng. Một nghiên cứu [4] về cân bằng vật chất dinh dưỡng ở 3 ao nuôi cá Tra có diện tích 1000 m2, mật độ thả trung bình là 48 con/m2, thời gian nuôi là 8 tháng cho thấy năng suất nuôi trung bình 28,8 tấn/1000 m2 (hay 288 tấn/ha), thức ăn được sử dụng là thức ăn viên công nghiệp nổi, hệ số chuyển hóa thức ăn (FCR) là 1,56. Kết quả xác định sự phân Trang | 17 bố vật chất trong ao nuôi cho thấy cá tích lũy 32,6% DM từ nguồn thức ăn, lượng vật chất thải ra môi trường ở dạng hòa tan là 5,04%, dạng chất rắn là 45,6%, và vật chất khô bị mất đi do bốc hơi và rò rỉ là 16,7%. Tương tự cá tích lũy 42,6% N từ thức ăn; 5,44% N thải ra môi trường dạng hòa tan, 50,4% N dạng rắn (bùn đáy ao và vật chất lơ lửng trong nước) và 1,5% N bốc hơi và rò rỉ.

Tỉ lệ tích lũy phốt pho (P) từ nguồn thức ăn của cá tương đối thấp, 29,3%, thấp hơn so với tỉ lệ tích lũy N và DM; và lượng P thải dưới dạng hòa tan là 1,8% và dạng rắn là 65,57%. Nghiên cứu [4] cho thấy sản xuất 1 kg cá Tra, trung bình cần phải cung cấp cho cá 1,56 kg thức ăn viên công nghiệp, trong đó tổng DM khoảng 1.420 g, cá sử dụng thức ăn và đồng hóa khoảng 460 g (tích lũy trong cơ thể), phần còn lại khoảng 960 g vật chất khô là chất thải thải ra môi trường dưới dạng thức ăn dư thừa tích tụ ở đáy ao và chất thải của cá thải ra qua quá trình hô hấp (CO2, NH3) và bài tiết (phân và muối dinh dưỡng). Về N, để sản xuất 1 kg cá Tra thì thức ăn cần cung cấp 43,8 g N, cá sử dụng khoảng 18,3 g cho cơ tể và phần còn lại khoảng 25,5 g sẽ tích tụ trong bùn đáy ao hoặc hòa tan trong nước. Các chất đạm hòa tan bao gồm NH3, NH4+, NO2- và NO3- có từ quá trình bài tiết của cá và quá trình phân hủy bùn hữu cơ.

Các chất đạm tích lũy trong bùn là các chất đạm hữu cơ có từ thức ăn thừa và phân thải của cá. Tương tự, sản xuất ra 1 kg cá Tra, trong thức ăn cần cung cấp khoảng 18 g P, cá hấp thụ và tích lũy trong cơ thể khoảng 5,23 g, các chất thải của cá và thức ăn dư thừa chứa khoảng 12,8 g P, phần lớn P chứa trong chất thải rắn bị tích tụ trong bùn, một phần nhỏ hòa tan vào trong nước dưới dạng orthophosphate. Nghiên cứu [10] đã ước tính lượng chất thải N và P trên các ao nuôi cá Tra sử dụng hai loại thức ăn viên công nghiệp (commercial pellet) và thức ăn tự chế (farm-made feed), kết quả cho thấy lượng chất thải N gần tương đương nhau 46 kg/tấn cá (46 g/kg cá) và 46,8 kg/tấn cá (46,8 g/kg cá) cho hai loại thức ăn viên công nghiệp và thức ăn tự chế. Tương tự, lượng P chất thải lần lượt là 14,4 kg/tấn cá và 18,4 kg/tấn cá khi dùng thức ăn viên công nghiệp và thức ăn tự chế.

Như vậy, nếu năm 2014 cả nước sản xuất 1,116 triệu tấn cá Tra [11], ước tính lượng chất thải N và P theo kết quả của [4] thì tổng lượng chất thải thải ra môi trường 28. Tương tự, nếu ước tính dựa theo kết quả nghiên cứu của [10] thì lượng chất thải thải ra môi trường là 51. Như vậy, hàm lượng dinh dưỡng tích tụ trong ao nuôi cá Tra thâm canh cao hơn nhiều lần so với ao nuôi thâm canh của các đối tượng nuôi thủy sản khác. Chất lượng nước ao nuôi cá Tra thâm canh hầu hết là vượt quá giới hạn về chất lượng nước dùng cho nuôi trồng thủy sản.

Nhiều yếu tố vượt giới hạn của Quy chuẩn kỹ thuật quốc gia về tiêu chuẩn nước mặt (QCVN 08:2008/BTNMT). Yếu tố TSS và TAN, TN và TP thường xuyên vượt mức A và giai đoạn cuối vụ nuôi hầu như vượt mức B của Quy chuẩn kỹ thuật quốc gia về nước thải công nghiệp (QCVN 40:2011/BTNMT). Hàm lượng vật chất dinh dưỡng tích lũy trong bùn đáy cũng cao, cao gấp nhiều lần so với các thủy vực nuôi trồng thủy sản khác. Vì vậy, lượng chất thải thải ra môi trường cũng vô cùng lớn, vấn đề này cần được quan tâm trong thời gian tới.

Đảm bảo cho sự phát triển bền vững của nghề nuôi cá Tra và bảo vệ môi trường cần nghiên cứu đánh giá tải lượng của các thủy vực ở các vùng nuôi cá Tra thâm canh, quy hoạch vùng nuôi hợp lý và xây dựng các quy trình xử lý chất thải trước khi thải ra môi trường xung quanh.2 Dinh dưỡng và thức ăn cá Tra Trong nuôi thủy sản thức ăn chiếm tỉ lệ cao trong tổng chi phí chung từ 50-80%. Trước đây cá Tra chủ yếu được nuôi bằng thức ăn tự chế nhưng từ năm 1995 trở lại đây hầu hết cá Tra được nuôi bằng thức ăn viên.