Luận án tiến sĩ 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ĩ nghiên cứu 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, phát triển phương pháp mới, đánh giá hiệu quả ứng dụng

• 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 ảnh hưởng đến môi trường ao nuôi cá tra Việt Nam
• Loại tài liệu: luận án
• Đị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
• Mục đích nghiên cứu
• Đối tượng và phạm vi nghiên cứu
• Đối tượng nghiên cứu
• Phạm vi nghiên cứu
• Nội dung nghiên cứu
• Các đóng góp của luận án
• Ý nghĩa khoa học và thực tiễn
• 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.2. Đá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.2.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.2.2. Hiệu quả của bổ sung enzyme phytase lên môi trường
• 3.3. Số liệu chạy mô hình
• 3.4. Kết quả phân tích độ nhạy và hiệu chỉnh mô hình
• 3.4.1. Phân tích độ nhạy
• 3.5. Mô phỏng quá trình phát triển động thực vật phù du trong ao cá
• 3.6. 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.6.1. Quá trình phát triển của cá
• 3.6.1.1. Quá trình phát triển của cá tại ao nuôi thực địa
• 3.6.1.2. Mô phỏng quá trình phát triển của cá trong ao nuôi
• 3.6.2. Biến động ôxy hòa tan trong ao nuôi thực địa
• 3.6.3. Mô phỏng quá trình biến thiên của ôxy trong ao nuôi
• 3.6.4. Tổng nitơ và tổng phốt pho trong cột nước
• 3.6.4.1. Tổng nitơ và tổng phốt pho trong cột nước trong ao nuôi thực địa
• 3.6.4.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.6.5. Tổng nitơ và tổng phốt pho trong bùn lắng
• 3.6.5.1. Tổng nitơ và tổng phốt pho trong bùn lắng ao nuôi thực địa
• 3.6.5.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. Nghiên cứu mô phỏng

Nghiên cứu mô phỏng là phương pháp chính được sử dụng để phân tích chuyển hóa chất dinh dưỡng trong môi trường ao nuôi cá tra Việt Nam. Các mô hình mô phỏng được xây dựng dựa trên các phương trình toán học, mô tả các quá trình sinh học và hóa học xảy ra trong ao nuôi. Các mô hình này giúp dự đoán sự biến động của các chất dinh dưỡng như nitơ, phốt pho, và oxy hòa tan, từ đó đánh giá ảnh hưởng dinh dưỡng đến chất lượng nước và sức khỏe cá. Các kết quả mô phỏng cung cấp cơ sở khoa học để tối ưu hóa nuôi trồng và quản lý chất thải hiệu quả.

1.1. Mô hình hóa chất lượng nước

Mô hình hóa chất lượng nước là công cụ quan trọng trong nghiên cứu, giúp mô phỏng các quá trình chuyển hóa vật chất trong ao nuôi. Các mô hình như PCLake và các mô hình động lực dinh dưỡng được sử dụng để phân tích sự biến động của các thông số môi trường như DO, BOD, COD, TSS, và các chất dinh dưỡng chính. Kết quả mô phỏng cho thấy sự tương quan giữa lượng thức ăn, chất thải, và chất lượng nước, từ đó đề xuất các biện pháp quản lý ao nuôi hiệu quả.

1.2. Mô phỏng chu trình dinh dưỡng

Mô phỏng chu trình dinh dưỡng tập trung vào các quá trình chuyển hóa nitơ, phốt pho, và carbon trong ao nuôi. Các mô hình này giúp hiểu rõ cách thức các chất dinh dưỡng được tích lũy, chuyển hóa, và thải ra môi trường. Kết quả cho thấy, chỉ khoảng 37.5% chất dinh dưỡng từ thức ăn được cá tích lũy, phần còn lại thải ra môi trường, gây ô nhiễm. Điều này nhấn mạnh sự cần thiết của việc tối ưu hóa dinh dưỡng cá tra và quản lý chất thải.

II. Chuyển hóa chất dinh dưỡng

Chuyển hóa chất dinh dưỡng là quá trình trung tâm trong nuôi trồng thủy sản, đặc biệt là trong ao nuôi cá tra. Các chất dinh dưỡng như nitơ, phốt pho, và carbon được chuyển hóa thông qua các chu trình sinh học và hóa học, ảnh hưởng trực tiếp đến chất lượng nước và sức khỏe cá. Nghiên cứu chỉ ra rằng, việc sử dụng enzyme phytase trong thức ăn giúp cải thiện hiệu quả sử dụng phốt pho, giảm thiểu chất thải và tác động môi trường. Điều này có ý nghĩa quan trọng trong việc tối ưu hóa nuôi trồng và quản lý ao nuôi bền vững.

2.1. Chu trình nitơ và phốt pho

Chu trình nitơ và phốt pho là hai quá trình chính trong chuyển hóa chất dinh dưỡng. Nitơ được chuyển hóa qua các dạng như amoni, nitrit, và nitrat, trong khi phốt pho được tích lũy trong bùn đáy ao. Nghiên cứu cho thấy, việc bổ sung enzyme phytase giúp giảm lượng phốt pho thải ra môi trường, cải thiện chất lượng nước và giảm tác động môi trường. Điều này đóng góp vào việc quản lý chất thải hiệu quả hơn.

2.2. Ảnh hưởng của enzyme phytase

Enzyme phytase đóng vai trò quan trọng trong việc cải thiện hiệu quả sử dụng phốt pho trong thức ăn. Nghiên cứu chỉ ra rằng, việc bổ sung enzyme này giúp giảm lượng phốt pho thải ra môi trường, đồng thời cải thiện tốc độ tăng trưởng và tỷ lệ sống của cá. Điều này không chỉ nâng cao hiệu quả kinh tế mà còn giảm thiểu tác động môi trường, góp phần vào sự phát triển bền vững của nuôi trồng thủy sản.

III. Môi trường ao nuôi cá tra Việt Nam

Môi trường ao nuôi cá tra Việt Nam đang đối mặt với nhiều thách thức do tác động môi trường từ chất thải nuôi trồng. Nghiên cứu chỉ ra rằng, lượng chất thải từ ao nuôi cá tra gây ô nhiễm nguồn nước, đặc biệt là các chất dinh dưỡng như nitơ và phốt pho. Việc quản lý chất lượng nước và quản lý chất thải hiệu quả là yếu tố then chốt để đảm bảo sự bền vững của nghề nuôi cá tra. Các giải pháp như sử dụng enzyme phytase và tối ưu hóa dinh dưỡng được đề xuất để giảm thiểu tác động môi trường và nâng cao hiệu quả nuôi trồng.

3.1. Chất lượng nước ao nuôi

Chất lượng nước ao nuôi là yếu tố quyết định đến sức khỏe và năng suất của cá tra. Nghiên cứu cho thấy, các thông số như DO, BOD, COD, và TSS có sự biến động lớn trong quá trình nuôi, ảnh hưởng trực tiếp đến môi trường ao nuôi. Việc quản lý chất lượng nước thông qua các biện pháp như sử dụng enzyme phytase và tối ưu hóa dinh dưỡng giúp cải thiện chất lượng nước, giảm thiểu tác động môi trường.

3.2. Quản lý chất thải

Quản lý chất thải là vấn đề cấp bách trong nuôi trồng thủy sản, đặc biệt là trong ao nuôi cá tra. Nghiên cứu chỉ ra rằng, lượng chất thải từ ao nuôi gây ô nhiễm nguồn nước, đặc biệt là các chất dinh dưỡng như nitơ và phốt pho. Các giải pháp như sử dụng enzyme phytase và tối ưu hóa dinh dưỡng được đề xuất để giảm thiểu tác động môi trường và đảm bảo sự bền vững của nghề nuôi cá tra.

Nội dung chính

Tổng quan về luận án

Ngành nuôi cá Tra (Pangasianodon hypophthalmus) thâm canh tại Đồng bằng sông Cửu Long (ĐBSCL) đạt quy mô sản xuất vượt 1,42 triệu tấn trên diện tích 6.600 ha, đóng góp hơn 1,61 tỷ USD kim ngạch xuất khẩu nhưng đồng thời tạo áp lực sinh thái nghiêm trọng lên lưu vực sông Mekong. Thực trạng thâm canh hóa với mật độ thả cực cao (48–70 con/m²) và năng suất đạt 250–300 tấn/ha/vụ kéo theo sự tích lũy khổng lồ của các hợp chất chứa Nitơ (N), Phốt pho (P) và các chất hữu cơ lơ lửng.

[Thức ăn công nghiệp (100% N, P, DM)]

Khoảng trống nghiên cứu (Research Gap) cốt lõi nằm ở việc thiếu vắng các mô hình động thái sinh thái - hóa học có khả năng lượng hóa chi tiết cơ chế tương tác đa chiều giữa các chu trình sinh sinh học (C-N-P-DO), chuỗi thức ăn thủy sinh (thực vật phù du, động vật phù du) và sự tích tụ trầm tích đáy dưới tác động của các can thiệp enzyme dinh dưỡng sinh học. Các nghiên cứu trước đây (Boyd, 1998; Gross et al., 2000; Nguyễn Thanh Phương et al., 2010) chủ yếu dừng lại ở các quan trắc đơn lẻ hoặc cân bằng khối lượng tĩnh.

Luận án xác lập 03 câu hỏi nghiên cứu và 03 giả thuyết khoa học tương ứng:

• RQ1: Động thái biến thiên thủy lý - thủy hóa (DO, TAN, $\text{NO}_2^-$, $\text{NO}_3^-$, $\text{PO}_4^{3-}$, TSS) trong cột nước và bùn đáy ao cá Tra diễn ra theo quy luật nào dưới mật độ thâm canh cao? $\rightarrow$ H1: Sự suy giảm oxy hòa tan và tích lũy TAN/TP vượt ngưỡng chịu tải sinh thái của thủy vực do tải lượng thức ăn dư thừa và phân giải yếm khí ở trầm tích.
• RQ2: Can thiệp bổ sung enzyme phytase ngoại sinh vào khẩu phần thức ăn tác động thế nào đến hệ số chuyển hóa thức ăn (FCR) và bài tiết phốt pho? $\rightarrow$ H2: Enzyme phytase thủy phân hiệu quả liên kết myo-inositol hexakisphosphate, nâng cao độ khả dụng sinh học của P và giảm phát thải phốt pho ra nước ao trên 40%.
• RQ3: Có thể mô phỏng toán học chính xác các chu trình biến đổi chất dinh dưỡng và tăng trưởng sinh khối cá Tra bằng hệ phương trình vi phân phi tuyến giải số không? $\rightarrow$ H3: Mô hình tích hợp cân bằng khối lượng và động học Monod giải bằng thuật toán Runge-Kutta bậc 4 phản ánh đúng diễn biến thực địa với sai số RMSE tối thiểu.

Khung lý thuyết (Theoretical Framework) của luận án kết hợp Lý thuyết Cân bằng Dinh dưỡng Sinh học (Bioenergetics Theory), Động học Biến đổi Dưỡng chất (Nutrient Transformation Kinetics) và Mô hình Cân bằng Khối lượng Hồ sinh thái (Mass Balance Ecological Modeling). Luận án được thực hiện trên hệ thống ao nuôi thực địa tại các tỉnh trọng điểm (Cần Thơ, An Giang, Đồng Tháp) kết hợp thực nghiệm bể composite kiểm soát nghiêm ngặt, tạo ra bước đột phá trong quản trị môi trường thủy sản bền vững theo Quyết định 339/QĐ-TTg của Thủ tướng Chính phủ.

---

Literature Review và Positioning

Nghiên cứu về chuyển hóa chất dinh dưỡng và chất lượng nước ao nuôi thủy sản trên thế giới chia thành ba dòng nghiên cứu chính:

                  TIẾN TRÌNH NGHIÊN CỨU MÔ HÌNH HÓA THỦY SẢN
[1. Cân bằng Tĩnh]          [2. Động lực Dinh dưỡng]     [3. Lưới thức ăn Đa bậc]
Tucker & Hargreaves (2004)    Janse (2005) - PCLake        Polovina (1984) - Ecopath
Gross & Boyd (2000)           Chapra (1997) - Water Qual   Halfon (1996) - Ontario
                        [LUẬN ÁN LÊ XUÂN THỊNH]
        Tích hợp Sinh hóa C-N-P + Phytase Ngoại sinh + ODEs RK4 MATLAB

Dòng nghiên cứu 1: Mô hình cân bằng tải lượng khối lượng tĩnh (Mass-balance Box Models) tiêu biểu với các công trình của Boyd (1998), Gross et al. (2000) và Tucker & Hargreaves (2004) trên cá nheo Mỹ (Ictalurus punctatus). Các tác giả tập trung vào tính toán lượng N và P bài tiết cố định trên một đơn vị sinh khối sản xuất. Hạn chế của trường phái này là xem ao nuôi như một "hộp đen" đồng nhất, bỏ qua động học biến thiên phi tuyến theo thời gian thực và sự phân tầng chất lượng nước.

Dòng nghiên cứu 2: Mô hình hóa hiện tượng phú dưỡng và chu trình dinh dưỡng động học (Dynamic Nutrient Cycling Models) nổi bật với mô hình PCLake (Janse, 2005) áp dụng cho các hồ nông châu Âu và nghiên cứu phú dưỡng hồ Washington của Arhonditsis & Brett (2005). Các mô hình này giải quyết xuất sắc động học của C, N, P, oxy hòa tan và phiêu sinh vật, nhưng chưa được tham số hóa cho các thủy vực nhân tạo có mật độ nuôi siêu thâm canh và tốc độ trao đổi chất cực cao như ao cá Tra.

Dòng nghiên cứu 3: Mô hình cấu trúc mạng lưới thức ăn (Trophic Food Web Models) với phần mềm Ecopath của Polovina (1984), được cải tiến bởi Christensen & Pauly (1992), Halfon (1996) tại hồ Ontario (Canada), và Douglas (2003) trên hệ đầm lầy Olentangy River Wetlands (ORWs - Mỹ) bằng phần mềm STELLA. Các nghiên cứu này mô tả dòng năng lượng thông qua các bậc dinh dưỡng nhưng thiếu tính năng mô phỏng chi tiết các quá trình thủy hóa vô cơ tức thời.

                              BẢN ĐỒ ĐỊNH VỊ HỌC THUẬT
           Thấp             Mức độ chuyên biệt hóa ao Thâm canh Nhiệt đới  Cao

Tranh luận học thuật cốt lõi diễn ra giữa hai quan điểm đối nghịch:

1. Nhóm nghiên cứu sinh thái hồ cổ điển (Scheffer, 1998; Janse, 2005) cho rằng nồng độ phốt pho hòa tan ($\text{PO}_4^{3-}$) trong nước được kiểm soát chủ yếu bởi quá trình giải phóng/hấp phụ từ phức hợp Fe-P và Ca-P ở lớp bùn lắng đáy thông qua thế oxy hóa khử (ORP).
2. Nhóm công nghệ nuôi trồng thủy sản (Tucker, 2000; Piedrahita, 2003) lại lập luận rằng trong các hệ thống nuôi thâm canh cao độ, dòng phospho ngoại sinh từ thức ăn công nghiệp và hoạt động bài tiết phân cá mới là động lực chi phối áp đảo toàn bộ cân bằng dinh dưỡng.

Luận án đã định vị chính xác điểm giao thoa này: Thiết lập mô hình động học phi tuyến tích hợp chu trình dưỡng chất (N, P, C, DO) với sự can thiệp dinh dưỡng chính xác bằng enzyme phytase, giải quyết triệt để sự thiếu hụt mô hình tính toán chuyên sâu cho loài cá da trơn vùng nhiệt đới.

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

• So với hệ thống ao nuôi cá nheo Mỹ vùng Tây Bắc Mississippi (Tucker & Hargreaves, 2004) với hàm lượng TSS dao động 80–135 mg/L và TAN 0,4–2,7 mg/L, ao nuôi cá Tra ĐBSCL có tải lượng ô nhiễm vượt trội với TSS đạt 38–274 mg/L và TAN đạt 4,0–8,0 mg/L.
• So với mô hình năng lượng hồ Kelavarapalli (Ấn Độ) của Saravana Bhavan et al. (2009) vốn sử dụng phương pháp cân bằng tĩnh Ecopath, mô hình của luận án áp dụng hệ phương trình vi phân đạo hàm động học trên nền tảng MATLAB, cho phép dự báo biến thiên theo từng bước thời gian $\Delta t$ với độ nhạy tham số cao.

---

Đó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 học Sinh thái Hồ nông (Shallow Lake Ecological Theory) của Scheffer và Lý thuyết Năng lượng Sinh học Thủy sản (Fish Bioenergetics Theory) của Kitchell et al. Luận án bác bỏ giả định môi trường tĩnh trong các hệ thống nuôi thâm canh bằng việc xác lập hệ thống tương tác động học đa tầng:

               KHUNG LÝ THUYẾT VÀ TƯƠNG TÁC ĐỘNG HỌC ĐA TẦNG
                      
                        [BỨC XẠ MẶT TRỜI & NHIỆT ĐỘ]
 (Tảo lục, lam, cát)         (Copepoda, Cladocera)           (P. hypophthalmus)
 (Khoáng hóa kỵ khí: C, TN, TP bùn)

Các mệnh đề lý thuyết cốt lõi được phát triển:

• Mệnh đề P1: Động học giải phóng phốt pho nội sinh phụ thuộc phi tuyến vào hoạt tính phytase đường ruột và tỷ lệ bổ sung phytase ngoại sinh trong thức ăn, làm thay đổi tỷ lệ $P_{\text{hòa tan}} / P_{\text{trầm tích}}$ trong mô hình cân bằng khối lượng.
• Mệnh đề P2: Sự biến thiên oxy hòa tan (DO) trong ao thâm canh nhiệt đới chịu sự chi phối áp đảo của quá trình hô hấp thực vật phù du ban đêm và nhu cầu oxy bùn đáy (SOD), vượt qua tốc độ khuếch tán tái tạo oxy từ bề mặt khí quyển.

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

Khung phân tích của luận án tích hợp đồng thời 03 lý thuyết nền tảng:

1. Lý thuyết Động học Enzyme Michaelis-Menten (thủy phân phức hợp phytate-protein bởi phytase).
2. Lý thuyết Động học Giới hạn Tăng trưởng Monod (sự hấp thụ chất dinh dưỡng N, P của tảo lục, tảo lam và tảo cát).
3. Lý thuyết Động lực học Quần thể Lotka-Volterra mở rộng (tương tác dinh dưỡng giữa thực vật phù du, động vật phù du và cá Tra).

                            MA TRẬN TÁC ĐỘNG TƯƠNG HỖ
Ghi chú: (+) Tác động kích thích/Tăng; (-) Tác động ức chế/Tiêu thụ; (0) Không tác động trực tiếp; (--) Biến tự tương quan.

Điều kiện biên (Boundary Conditions): Mô hình xác lập cho các thủy vực nước ngọt nhiệt đới nông ($\text{độ sâu } h = 3,5 \text{--} 4,5\text{ m}$), biên độ nhiệt độ nước $27,5 \text{--} 32^\circ\text{C}$, tỷ lệ thay nước định kỳ $30 \text{--} 40%/\text{ngày}$ và chế độ cho ăn sử dụng thức ăn công nghiệp viên nổi ($22 \text{--} 28%\text{ crude protein}$).

---

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái Thực chứng luận hậu kỳ (Post-positivism), kết hợp phương pháp quan sát - đo đạc thực nghiệm định lượng cao với mô hình hóa số trị lý thuyết.

                           SƠ ĐỒ KHỐI THIẾT KẾ NGHIÊN CỨU
     [THỰC NGHIỆM ĐỐI CHỨNG]                       [MÔ HÌNH HÓA TOÁN HỌC]
   Ao thực địa & Bể Composite                     Thiết lập ODEs C-N-P-DO-Biomass
                            [PHÂN TÍCH ĐỘ NHẠY & HIỆU CHỈNH]
                             MSE, RMSE đối soát Data Thực
                            [MÔ HÌNH HOÀN THIỆN & ỨNG DỤNG]

Thiết kế đa tầng bao gồm:

• Cấp độ 1 (Micro-mesocosm): Thí nghiệm kiểm soát trên bể composite tại Đại học Cần Thơ nhằm đánh giá tác động sinh học của các mức enzyme phytase ($0$, $500$, $1000$, $1500$, $2000\text{ UI/kg}$ thức ăn) lên cá Tra giai đoạn giống và thương phẩm.
• Cấp độ 2 (Macro-field): Khảo sát, quan trắc dọc (longitudinal field monitoring) tại các ao nuôi thực địa đại diện ở Cần Thơ, An Giang, Đồng Tháp trong chu kỳ nuôi 6–8 tháng.

Quy trình nghiên cứu rigorous

Chiến lược thu mẫu và phân tích hóa nghiệm tuân thủ Tiêu chuẩn Chuẩn hóa Quốc tế (APHA, 2012; TCVN):

• Thu mẫu nước: Thu định kỳ 2 tuần/lần tại 3 tầng nước (tầng mặt $-0,5\text{ m}$, tầng giữa $-2,0\text{ m}$, tầng đáy cách bùn $0,5\text{ m}$) bằng thiết bị lấy mẫu nước chuyên dụng Niskin.
• Thu mẫu bùn đáy: Sử dụng gàu cuốc Ekman lấy mẫu trầm tích ở độ sâu $0 \text{--} 15\text{ cm}$ để phân tích tổng đạm bùn (TN bùn), tổng phốt pho bùn (TP bùn) và hàm lượng chất hữu cơ (OM).
• Quy trình phân tích:
  • Nhiệt độ, pH, DO, Độ dẫn điện (EC) đo trực tiếp bằng máy đo đa chỉ tiêu cầm tay TOA-DKK WQC-24D.
  • $\text{BOD}_5$ xác định bằng phương pháp ủ 5 ngày ở $20^\circ\text{C}$; $\text{COD}$ xác định bằng phương pháp chuẩn độ Bicromat.
  • $\text{TAN}$ phân tích bằng phương pháp so màu Phenate; $\text{NO}_2^-$ bằng phương pháp Diazot hóa; $\text{NO}_3^-$ bằng phương pháp khử Cadmium.
  • $\text{PO}_4^{3-}$ và $\text{TP}$ phân tích bằng phương pháp so màu Ascorbic Acid sau khi công phá bằng Persulfate; Chlorophyll-a chiết xuất bằng Acetone $90%$ và đo quang phổ.

Tam giác đạc (Triangulation): Dữ liệu được đối soát chéo giữa phân tích hóa nghiệm phòng thí nghiệm, số liệu đo tự động tại chỗ và kết quả mô phỏng số học từ phần mềm.

Data và phân tích

Hệ phương trình vi phân mô tả sự biến thiên các cấu tử được xây dựng dựa trên định luật bảo toàn khối lượng:

$$\frac{dC_i}{dt} = \sum \text{Nguồn sinh} - \sum \text{Tiêu hao} \pm \text{Vận chuyển biên (Trao đổi nước)}$$

Ví dụ phương trình vi phân biến thiên Oxy hòa tan ($\text{DO}$): $$\frac{d[\text{DO}]}{dt} = K_L ([\text{DO}]{\text{sat}} - [\text{DO}]) + P{\text{max}} \cdot f(I) \cdot [\text{TVPD}] - R_{\text{TVPD}} [\text{TVPD}] - R_{\text{DVPD}} [\text{DVPD}] - R_{\text{Fish}} [\text{Fish}] - \text{NOD} - \text{SOD}$$

Phương trình biến thiên Tổng Đạm Amôn ($\text{TAN}$): $$\frac{d[\text{TAN}]}{dt} = E_{\text{Fish}} + D_{\text{Sludge_N}} + \text{Hyd}{\text{DON}} - U{\text{TVPD}} - K_{\text{nit}} [\text{TAN}] \cdot \theta^{(T-20)} - \frac{Q_{\text{out}}}{V}[\text{TAN}]$$

Phương trình biến thiên Sinh khối cá Tra ($W$): $$\frac{dW}{dt} = \left( \alpha \cdot I_{\text{feed}} \cdot f(T) \cdot f(\text{DO}) - R_{\text{maint}} - \text{SDA} \right) W$$

Hệ phương trình vi phân phi tuyến gồm 12 biến trạng thái liên đới được lập trình và giải số bằng thuật toán Runge-Kutta bậc 4 (RK4) trên môi trường MATLAB R2018b:

$$y_{n+1} = y_n + \frac{h}{6}(k_1 + 2k_2 + 2k_3 + k_4)$$

Phân tích độ nhạy cục bộ (Local Sensitivity Analysis) được thực hiện để sàng lọc các tham số nhạy cảm nhất. Đánh giá độ chính xác hiệu chuẩn bằng Sai số bình phương trung bình (MSE) và Sai số căn bậc hai trung bình (RMSE):

$$\text{RMSE} = \sqrt{\frac{1}{N} \sum_{t=1}^N \left( Y_{\text{obs}}(t) - Y_{\text{sim}}(t) \right)^2}$$

---

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

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

                         KẾT QUẢ ĐỘT PHÁ CỦA ENZYME PHYTASE
Ghi chú: (*) Khác biệt có ý nghĩa thống kê ở mức P < 0,05.

1. Hiệu năng giải phóng phospho sinh học vượt bậc: Bổ sung enzyme phytase (từ Bacillus subtilis Ba 58) ở mức tối ưu $1500\text{ UI/kg}$ thức ăn (tương đương $0,01%$ trọng lượng khô) giúp nâng cao tỷ lệ tích lũy phospho (PHR) từ $47,83 \pm 2,31%$ lên $66,46 \pm 2,38%$ ($P < 0,05$). Trích dẫn nguyên văn dữ liệu luận án: "Thức ăn có bổ sung phytase (1500 U/kg thức ăn) giúp giảm thải ra ngoài môi trường nước 10,39 ± 0,74 g P/kg WG (47,61%) so với không bổ sung phytase trong thức ăn tăng lên 21,82 ± 0,96 g P/kg WG (100%)".
2. Cải thiện triệt để chỉ số FCR và PER: FCR giảm mạnh từ $1,66 \pm 0,09$ xuống $1,29 \pm 0,06$, đồng thời hiệu quả sử dụng protein (PER) tăng từ $1,58 \pm 0,08$ lên $2,01 \pm 0,09$ ($P < 0,05$), chứng minh phytase không chỉ giải phóng liên kết phytate-P mà còn phá vỡ cấu trúc phức hợp phytate-protein, giải phóng axit amin tự do cho cá hấp thụ.
3. Lượng hóa sự mất cân đối chuyển hóa vật chất: Luận án chỉ rõ thực trạng: "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%, trong đó vật chất khô: nitơ, phốt pho tích lũy trong cá chiếm lần lượt là 32,6%; 42,7%; 29,8%. Trong khi vật chất khô thải ra môi trường chiếm đến 67,4% (5% trong nước, 45,63% trong bùn đáy và 16,74% mất đi do bay hơi hoặc thẩm thấu)".
4. Quy luật biến thiên oxy hòa tan và ô nhiễm cực đoan: Hàm lượng DO tầng mặt dao động dữ dội từ $< 2\text{ mg/L}$ ($25%$ bão hòa) vào lúc sáng sớm (6:00) vọt lên $> 15,5\text{ mg/L}$ ($200%$ bão hòa) vào xế chiều (15:00). Vào cuối vụ nuôi, tổng đạm amôn (TAN) đạt ngưỡng kỷ lục $8,0\text{ mg/L}$, Nitrat ($\text{NO}_3^-$) đạt $19,5\text{ mg/L}$ và Tổng Nitơ (TN) xấp xỉ $60\text{ mg/L}$, vượt xa Quy chuẩn QCVN 08:2008/BTNMT và mức B của QCVN 40:2011/BTNMT.
5. Độ tin cậy của mô hình số mô phỏng: Kết quả chạy mô phỏng trên MATLAB khớp với số liệu đo thực địa với hệ số tương quan $R^2 > 0,88$ và chỉ số sai số RMSE $< 12%$ đối với toàn bộ các biến trạng thái DO, TAN, $\text{NO}_3^-$, $\text{PO}_4^{3-}$, TN nước, TP nước và sinh khối bùn đáy.

Implications đa chiều

• Về mặt Lý thuyết: Cung cấp bộ tham số động học chuẩn hóa (Standardized Kinetic Parameters) cho hệ sinh thái ao nuôi cá da trơn nhiệt đới, làm phong phú thêm cơ sở dữ liệu sinh thái học quốc tế.
• Về mặt Phương pháp: Thiết lập quy trình tích hợp giữa mô hình số trị MATLAB và các công cụ thực nghiệm mesocosm, tạo khuôn mẫu cho nghiên cứu mô phỏng chất lượng nước trong nuôi trồng thủy sản.
• Về mặt Thực tiễn: Đưa ra giải pháp trực tiếp cắt giảm $15 \text{--} 20%$ chi phí thức ăn thông qua việc hạ thấp FCR, đồng thời giảm chi phí điện năng sục khí nhờ dự báo chính xác pha suy giảm oxy hòa tan.
• Về mặt Chính sách: Cung cấp cơ sở khoa học để Bộ Nông nghiệp và Môi trường hoàn thiện các quy chuẩn xả thải nước thải thủy sản (thay thế/bổ sung QCVN 02-20:2014/BNNPTNT và QCVN 40:2011/BTNMT).

---

Limitations và Future Research

1. Giới hạn không gian và phân tầng thủy lực: Mô hình toán giả định nước ao ở trạng thái hòa trộn hoàn toàn nhờ sự xáo trộn của cá (Completely Stirred Tank Reactor - CSTR), chưa xét đến vi phân không gian 3 chiều (3D Hydrodynamics) trong các ao có diện tích lớn $> 1\text{ ha}$ hoặc độ sâu $> 5\text{ m}$.
2. Biến động hệ vi sinh vật bản địa: Nghiên cứu chưa tích hợp động học quần xã vi khuẩn chuyển hóa N (nhóm Nitrosomonas, Nitrobacter) ở mức độ giải trình tự gen mà chỉ lượng hóa qua các hằng số tốc độ phản ứng bậc một ($K_{\text{nit}}$).
3. Điều kiện biên về khí tượng cực đoan: Bộ tham số mô hình được hiệu chuẩn trong điều kiện khí hậu nhiệt đới gió mùa điển hình của ĐBSCL, cần tái hiệu chỉnh khi áp dụng cho các vùng khí hậu có mùa đông lạnh như Đồng bằng Sông Hồng.

Chương trình nghiên cứu tiếp nối (Future Agenda):

• Tích hợp module Thủy động lực học 3D (sử dụng nền tảng Delft3D hoặc EFDC kết hợp MATLAB).
• Ứng dụng kỹ thuật phân tích eDNA và metagenomics để mô hình hóa chính xác vai trò của hệ vi sinh vật đáy trong chu trình phản nitrat hóa.
• Xây dựng hệ thống Hỗ trợ Ra Quyết định thời gian thực (Real-time Decision Support System - DSS) kết nối mạng lưới cảm biến IoT đo online với mô hình số để tự động điều khiển van thay nước và quạt nước.

---

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

                               CHUỖI TÁC ĐỘNG ĐA LĨNH VỰC

• Tác động Học thuật: Đóng góp bài toán mô hình hóa sinh thái nhiệt đới có chất lượng cao, tiềm năng trích dẫn lớn trong các tạp chí chuyên ngành hàng đầu như Aquacultural Engineering, Ecological Modelling, Aquaculture.
• Chuyển đổi Ngành công nghiệp: Với quy mô sản xuất 1,42 triệu tấn cá Tra/năm, việc áp dụng công thức bổ sung phytase $0,01%$ giúp cắt giảm hàng chục ngàn tấn phốt pho thải ra hệ thống sông rạch ĐBSCL, tiết kiệm ước tính hơn 1.200 tỷ đồng chi phí thức ăn toàn ngành nhờ giảm FCR từ 1,66 xuống 1,29.
• Hỗ trợ Chứng nhận Sinh thái: Giúp các vùng nuôi đáp ứng nghiêm ngặt các tiêu chuẩn quốc tế (ASC, GlobalGAP, BAP), tháo gỡ các rào cản kỹ thuật khi xuất khẩu sang các thị trường khó tính như EU và Hoa Kỳ.

---

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

• Nghiên cứu sinh / Học viên Cao học: Kế thừa cấu trúc mã nguồn mô hình vi phân MATLAB, phương pháp giải số RK4 và bộ hằng số động học để mở rộng cho các đối tượng thủy sản khác (tôm thẻ chân trắng, cá rô phi, cá chẽm).
• Các Nhà Khoa học Thủy sản & Môi trường: Tiếp cận phương pháp luận liên ngành kết hợp giữa hóa học môi trường, enzyme học và mô hình toán tin.
• Doanh nghiệp Sản xuất Thức ăn Thủy sản (Feed Mills): Ứng dụng quy chuẩn phối trộn enzyme phytase ($1500\text{ UI/kg}$) nhằm giảm tỷ lệ bột cá đắt tiền và giảm bổ sung dicanxi photphat (DCP) vô cơ.
• Chủ Trang trại Nuôi cá Tra: Áp dụng lịch trình thay nước, quạt nước tối ưu và kiểm soát bùn đáy theo dự báo của mô hình nhằm hạ giá thành nuôi, hạn chế dịch bệnh vi khuẩn (Aeromonas, Edwardsiella ictaluri).
• Cơ quan Quản lý Môi trường (Sở TN&MT, Chi cục Thủy sản): Sử dụng mô hình để đánh giá sức tải môi trường (carrying capacity) của các phân vùng nuôi tập trung dọc sông Tiền và sông Hậu.

---

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?

Đóng góp lý thuyết độc đáo nhất là việc tích hợp thành công động học thủy phân phytate của enzyme phytase ngoại sinh vào hệ phương trình vi phân động thái phốt pho trong mô hình sinh thái ao nuôi thâm canh. Công trình này mở rộng Lý thuyết Động thái Dinh dưỡng Hồ nông của Janse (mô hình PCLake) bằng cách bổ sung cơ chế kiểm soát dinh dưỡng nội sinh từ đường tiêu hóa của cá, phá vỡ thế bế tắc của các mô hình sinh thái truyền thống vốn chỉ xem sinh khối cá là biến số tiêu thụ thụ động.

2. Đột phá phương pháp luận so với các nghiên cứu quốc tế trước đây?

So với nghiên cứu cân bằng năng lượng hồ Ontario của Halfon (1996) sử dụng ma trận năng lượng tĩnh và nghiên cứu đầm lầy của Douglas (2003) trên phần mềm STELLA với thuật toán Euler sai số lớn, luận án đã:

1. Xây dựng hệ phương trình vi phân phi tuyến liên đới bậc cao biểu diễn đầy đủ 12 biến trạng thái sinh hóa.
2. Ứng dụng thuật toán số trị Runge-Kutta bậc 4 (RK4) lập trình trên MATLAB, đảm bảo độ hội tụ cao và triệt tiêu sai số tích lũy qua chu kỳ mô phỏng dài 180–240 ngày.
3. Thực hiện quy trình đối chuẩn kép (dual-calibration) đồng thời trên hệ thống bể mesocosm kiểm soát và hệ thống ao thực địa quy mô lớn.

3. Phát hiện bất ngờ nhất có số liệu thực nghiệm chứng minh?

Phát hiện bất ngờ nhất là sự mất cân đối cực đoan trong phân bổ dinh dưỡng bùn đáy và hiệu quả nghịch đảo của việc thay nước ở cuối vụ. Dù người nuôi tăng tỷ lệ thay nước lên đến $30 \text{--} 50%/\text{ngày}$, hàm lượng TAN và TN trong nước vẫn tiếp tục tăng vọt lên mức đỉnh điểm ($8,0\text{ mg/L}$ TAN và $60\text{ mg/L}$ TN). Dữ liệu phân tích cho thấy bùn đáy tích lũy tới $50,4%$ lượng Nitơ và $65,57%$ lượng Phốt pho xả thải. Vận động bơi lội của đàn cá mật độ cao ở tầng đáy liên tục khuấy đảo lớp trầm tích, tạo ra dòng hồi lưu dinh dưỡng nội sinh hòa tan ngược vào cột nước lớn hơn nhiều so với lượng chất bẩn được loại bỏ qua đường thoát nước.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?

Luận án cung cấp đầy đủ và chi tiết toàn bộ quy trình tái lập khoa học:

• Bảng tham số động học đầy đủ với giá trị cực đại, giá trị cực tiểu và giá trị hiệu chuẩn cuối cùng ($K_{\text{half}}$, tốc độ tăng trưởng cực đại $\mu_{\text{max}}$, hệ số hô hấp $R$, hằng số phân hủy $K_d$).
• Cấu trúc ma trận tác động tương hỗ giữa các cấu tử sinh thái.
• Toàn bộ mã nguồn thuật toán giải số Runge-Kutta bậc 4 lập trình trên ngôn ngữ MATLAB được công bố tại phần phụ lục.
• Quy trình chuẩn hóa pha trộn enzyme phytase vào thức ăn và quy trình quan trắc lấy mẫu phân tích 3 tầng nước theo tiêu chuẩn APHA.

5. Định hướng chương trình nghiên cứu 10 năm tới (10-Year Research Agenda)?

1. Giai đoạn 2026–2028: Phát triển mô hình tích hợp Thủy động lực học 3D - Sinh thái - Trầm tích (Coupled 3D Hydrodynamic-Ecological-Sediment Model) cho các hồ nuôi quy mô công nghiệp.
2. Giai đoạn 2028–2031: Tích hợp module Trí tuệ Nhân tạo / Học máy (AI/Machine Learning) nhằm tự động đồng hóa dữ liệu viễn thám và trạm IoT quan trắc nước online vào mô hình số.
3. Giai đoạn 2031–2035: Mở rộng mô hình sinh thái chuỗi giá trị tích hợp: Tính toán dấu chân Carbon (Carbon Footprint) và dấu chân Nước (Water Footprint) của toàn bộ chuỗi ngành hàng cá Tra Việt Nam phục vụ mục tiêu Net-Zero 2050.

---

Kết luận

1. Lượng hóa toàn diện cân bằng dưỡng chất: Xác lập chính xác tỷ lệ phân bố vật chất khô (DM), Nitơ (N) và Phốt pho (P) trong ao cá Tra thâm canh, chứng minh thực nghiệm rằng $67,4%$ DM, $57,3%$ N và $70,2%$ P từ thức ăn bị thải trực tiếp ra môi trường nước và bùn đáy.
2. Tối ưu hóa can thiệp sinh học Phytase: Chứng minh tỷ lệ bổ sung enzyme phytase $1500\text{ UI/kg}$ thức ăn ($0,01%$ trọng lượng khô) mang lại hiệu quả vượt trội: giảm FCR từ $1,66$ xuống $1,29$, nâng PER từ $1,58$ lên $2,01$, tăng tỷ lệ tích lũy phospho từ $47,83%$ lên $66,46%$, và giảm $47,61%$ lượng P bài tiết ra môi trường nước ($P < 0,05$).
3. Đột phá mô hình hóa số trị: Xây dựng thành công mô hình toán học giải số trên MATLAB bằng thuật toán Runge-Kutta bậc 4 mô phỏng chính xác động học C-N-P-DO, biến thái quần xã phiêu sinh vật và tăng trưởng sinh khối cá Tra với độ tin cậy cao ($R^2 > 0,88$, $\text{RMSE} < 12%$).
4. Định hình bước chuyển Paradigm: Chuyển đổi tư duy quản lý ao nuôi từ "xử lý ô nhiễm thụ động cuối đường ống" sang "kiểm soát chủ động từ dinh dưỡng đầu vào kết hợp dự báo động thái sinh thái theo thời gian thực".
5. Mở ra các hướng nghiên cứu liên ngành mới: Thiết lập nền móng cho 03 hướng nghiên cứu tiếp nối: Hệ thống tự động hóa điều tiết môi trường ao nuôi dựa trên IoT-DSS, Mô hình hóa phát thải khí nhà kính ($\text{CH}_4$, $\text{N}_2\text{O}$) từ bùn đáy ao thủy sản thâm canh, và Công nghệ sinh học enzyme thế hệ mới trong kinh tế tuần hoàn thủy sản.

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ỗ Trang | 16 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à 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. Trang | 17 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 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 (CO 2, 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 NH 3, 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 Trang | 18 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.