CN101731172B - 一种长吻鮠养殖动态投饲表的建立方法 - Google Patents

Abstract

本发明公开了一种长吻鮠养殖动态投饲表的建立方法，其步骤是：1.通过生长实验，研究摄食水平、体重和水温与能量收支式中最大摄食率、排粪率、排泄率、体增热、标准代谢和活动代谢的动态关系；2.在生物能量学模型的基础上，建立长吻鮠最适摄食率模型；3.利用计算机编程，输入养殖天数、水温、鱼初始体重、鱼体初始能值、饲料蛋白水平和饲料能值，经过计算，输出日体重生长、日最大摄食量和最适摄食量，形成长吻鮠养殖动态投饲表。这种动态投饲表能够及时地、精准地预测长吻鮠养殖过程中不同环境下的日生长和日摄食量，既不会造成饲料的浪费而增加饲料成本和饲料的次生污染，也不会因饲料投喂不足而影响鱼的生长。

Description

一种长吻鮠养殖动态投饲表的建立方法

技术领域

本发明涉及水产养殖业技术领域，更具体涉及一种长吻鮠养殖动态投饲表的建立方法，尤其适用于长吻鮠养殖。 

背景技术

准确、适宜的饲料投喂量是水产养殖业的关键因子之一，也是养殖技术中最重要的一环，是降低饲料系数的关键因素。投饲量不足，鱼处于半饥饿状态，不能满足鱼类能量和营养需要，生长发育缓慢，甚至使鱼类不能维持体重而减产，严重影响水产养殖效益；投饲量过大，不但饲料利用率低，造成饲料浪费，加大残饵对养殖水体的再次污染，而且病害增多，养殖效益大幅下降。 

目前，利用配合饲料，进行高投入高产出的长吻鮠集约化水产养殖中，水产饲料投喂量的确定基本是依据养殖户的经验，常用方法主要有以下几种： 

1.以鱼类净增重倍数和饵料系数计算水产饲料的年投饲量、月投饲量、日投饲量。饵料系数定义为饵料消耗量与鱼体增重量的比值。(1)年投饲量：根据鱼类净增重倍数和饵料系数进行推算，即鱼种放养量×净增重倍数×饵料系数。如果几种饲料交替使用，则分别以各自的饵料系数计算出使用量，然后相加即为年投饲量；(2)月投饲量：年投饲量×当月饲料分配百分比。一般春季放养鱼种，从3月份开始投喂至11月份结束，期间每月投饵量比例约为：3月1％、4月4％、5月8％、6月15％、7月20％、8月20％、9月20％、10月9％和11月3％；(3)日投饲量：根据月投饲量分上、中、下三旬安排。3～8月，上旬的日投饲量是当月投饲量日平均数的80％，中旬为平均数，下旬为平均数的120％；从9月起，上旬的日投饲量为当月投饲量日平均数的120％，下旬为平均数的80％。 

2.测量鱼类存塘数，乘以不同水温下的日投饲率，计算日投饲量。一般水温15～20℃时，日投饲量为鱼体重的1％～3％；水温20℃以上时，日投饲量为鱼体重的3％～5％。采用此法，须准确求得存塘鱼重量后，才可精确计算出日投饲量。 

3.日投饲量的调整。日投饲量的调整主要根据季节、天气、水温、水质和鱼摄食情况灵活调整(谢信桐，鱼用配合饲料投喂技术，江西饲料，2004年第2期，23-25页)。(1)根据养殖季节调整。一年养殖周期之中的饲料投喂主要依据“早开食，抓中间，带两头”的规律，集中在6～9月份，4月份以前投喂工作尽量提前，10月份以后，应延长投喂，做到上市前停食，维持养殖对象的体重；(2)根据天气情况调整。天气晴朗，水中溶氧量高，鱼群摄食旺盛，应适应多投；反之，天气闷热，连续阴雨，水中溶氧量低，鱼群食欲不振，应少投或不投；(3)根据池塘水温情况调整。鱼类摄食量显著受到水温变化的影响。在适温范围内，水温升高对养殖鱼摄食强度有显著的促进作用；水温降低，鱼代谢水平下降，食欲减退。在高温季节超过适宜温度时，鱼类食欲减退，应减少投饵量；(4)根据池塘水质情况调整。水质清爽，鱼群摄食旺盛，应多投；水质不好，过肥、过浓，鱼群食欲不振，而且残饵容易使池水变坏，应少投；水质很坏，鱼已浮头时，应禁止投喂；(5)根据养殖鱼摄食情况调整。每次投饵量一般以鱼吃到七、八成饱为准，大部分鱼吃饱游走，仅有少量鱼在表层索饵。 

国内现在主要使用的这些投饲量的确定及调整的方法是建立在养殖户长期经验的基础上，只是生产中的经验总结，是一个指导性的概念，人为主观因素影响很大，不能真实的反映养殖鱼类实际需要的饲料量，容易受到环境变化和发育阶段的影响而导致较大的误差；而且由于缺乏科学、合理的投饲依据，不能及时、准确地调整日投饲量，造成饲料开支大，饲料系数高，饲料对养殖水体的次生污染、鱼病暴发的机会以及死亡率的增加等问题。不仅影响到水产品产量和养殖效益，而且带来环境保护和健康养殖的问题。 

国外现在确定养殖鱼类投饲量比较先进的方法主要有两种：利用经验模型和生物能量学模型来预测养殖鱼类的投饲量。 

1.经验模型。经验模型主要指用逐步回归分析方法来模拟鱼类摄食和不同环境因子之间的直接关系而得到的模型。已经在不少鱼类中建立了经验模型，如河鳟、湖鳟、虹鳟、大鳞大马哈鱼和大西洋鳟鲑(Cho，Feeding systems for rainbow trout and other salmonids with reference to current estimates of energy and protein requirements.Aquaculture，1992：100，107-123)。高首鲟养殖的适宜投饲量可以通过经验模型得到，与环境因子的主要关系为：lnR_opt＝-2.88-0.25lnW+0.4T-0.0077T²，式中T为水温(℃)，W为体重(g)，R_opt为适应投饲率(％BW/d)(Cui & Hung，A prototype feeding-growth table for white sturgeon.Journal of Applied Aquaculture，1995：5(4)，25-34)。 

经验模型在模型归纳的数据范围之内预测比较精确。但是这种模型存在着严重的缺陷： 

(1)在模型成立的条件之外，外推结果会导致严重偏差；(2)无法估算代谢废物排泄量；(3)当饲料和养殖环境改变时，投饲量无法及时调整；(4)主要基于体长和鲜重的增加或饲料转化效率，而不是结合鱼体蛋白和能量储存的饲料可利用能及饲料营养物含量，不适用于现代水产养殖的高能量饲料。 

2.生物能量学模型。早在二十世纪八十年代，北美就已经发展出应用生物能量学方法确定养殖鱼类的饲料投喂量，即根据鱼类的总能量需求评估鱼类的日摄食，进而评估鱼类的日投饲量(郭学武和唐启升，鱼类摄食量的研究方法，海洋水产研究，2004年25卷第1期，68-78页)。生物能量学是研究能量在生物体内转换的学科。生物能量学模型是根据生物能量学原理建立的预测鱼类生长和摄食的模型，这类模型的基本方程为能量收支式：G＝C-F-U-SDA-R_s-R_a。式中，G为贮存在鱼体内的能量(生长能)，C为从食物中获取的能量(摄食能)，F为从粪便中损失的能量(排粪能)，U为从排泄物中损失的能量(排泄能)，SDA为与食物在体内转换、利用有关的能量消耗(体增热)，R_s为在饥饿、静止状态下的能量消耗(标准代谢)，R_a为与游泳等活动有关的能量消耗(活动代谢)。研究能量收支式各组分与影响摄食和生长的主要因子(摄食水平、体重和水温)间的关系，及各组分之间的相互关系，是鱼类生物能量学模型的主要内容(崔奕波，鱼类生物能量学的理论与方法，水生生物学报，1989年13卷第4期，369-383页)。 

以美国鱼类学家Kitchel为代表的威斯康星大学的研究者对鱼类生物能量学及其模型的发展做出了杰出贡献(Kitchell等，Application of a bioenergetics model to yellow perch and walleye.Journal of the Fisheries Research Board of Canada，1977：34，1922-1935)，其威斯康星模型(Wisconsin Model)已计算机程序化，在北美的渔业研究中得到了广泛应用(Ney，Bioenergetics modeling today：Growing pains on the cutting edge.Transactions of the American Fisheries Society，1993：122，736-748.)。威斯康星模型的一项重要贡献就是，鱼类在任何体重和温度条件下的日摄食量(C)皆可以用其最大日摄食量(C_max)来表达，即：C＝C_max×P×rc。式中，C_max＝a×bW是特定体重(W)的鱼在最适温度下的日摄食量，a和b是回归常数；P是均衡常数，取值范围0到1，用以调整日摄食量以符合体重生长的观测曲线；rc是温度标量，变化范围也是从0到1。在威斯康星模型的基础上，加拿大的安大略省自然资源部(Ontario Ministry of Natural Resources)开发出了基于生物能量学模型的计算机模拟系统Fish-PrFEQ Program。这套系统可以预测不同条件下养殖鱼类的生长和摄食，目前已经应用到鲑鳟和尖吻鲈等鱼类中。 

国外已经广泛应用在渔业中的利用生物能量学模型建立养殖鱼类投喂表的方法，存在 以下缺陷：(1)很多子模型是估算的或来自不同实验室的数据，而非通过直接的实验数据所得，导致模型的预测结果误差较大；(2)很多模型忽略鱼体能值随发育阶段而变化；(3)这些生物能量学模型的应用基本针对鲑鳟鱼类，不能应用在其他鱼类中。 

本发明的目的是在生物能量学模型的基础上，提供一种长吻鮠养殖动态投饲表，准确地预测不同养殖环境下，长吻鮠养殖过程中的日摄食量，形成长吻鮠养殖的日投饲表，显著提高饲料效率，降低饲料成本和渔业污染。 

发明内容

本发明的目的是在于提供了一种长吻鮠养殖动态投饲表的建立方法，方法易行，操作简便，长吻鮠的生物能量学模型严格建立在实验数据的基础上，由此形成的长吻鮠动态投饲表能够及时地、精准地预测养殖过程中不同水温下，不同生长阶段长吻鮠的日生长和日摄食量。 

为了实现上述的目的，本发明采用以下技术措施： 

通过生长实验，研究能量收支式各项与摄食水平、水温和体重之间的回归关系，建立长吻鮠的生物能量学模型和最适摄食量模型，预测养殖过程中的日生长和日摄食量，进而形成长吻鮠养殖动态投饲表。 

一种长吻鮠养殖动态投饲表的建立方法，其步骤如下： 

1.生长实验的条件 

1.1 摄食水平对长吻鮠摄食、生长和能量收支的影响 

经过8周的生长实验，体重为6.6克的长吻鮠幼鱼在饥饿及初始体重的0.8％/天、1.6％/天、2.4％/天和饱食摄食水平下的摄食、生长和能量收支。 

1.1.1 特定生长率与摄食水平 

长吻鮠的生长用特定生长率表示：特定生长率(SGR，％/天)＝100×(lnW_t-lnW_o)/t，式中W_t(克)和W_o(克)分别为终末和初始体重，t(天)为实验持续时间。特定生长率SGR与摄食水平RL(％体重/天)之间的关系为： 

SGR＝-0.66+2.65×ln(RL+1)，n＝15，R²＝0.9998。 

1.1.2饲料转化效率与摄食水平 

饲料利用用饲料转化效率表示，饲料转化效率(FCE，％)＝100×(W_t-W_o)/FI，式中W_t(克)和W_o(克)分别为终末和初始体重，摄食量FI(克)为长吻鮠摄食饲料的重量。饲料转化效率FCE与摄食水平RL(％体重/天)之间的关系为： 

FCE＝97.1+69.3×RL-28.9×RL²，n＝15，R²＝0.675； 

1.1.3 能量收支与摄食水平 

饱食水平下，长吻鮠的能量收支式为：100C＝12.9×F+5.5×U+40.1×R+41.5×G，式中C为摄食能(千焦/尾/天)，F为粪便能(千焦/尾/天)，U为排泄能(千焦/尾/天)，R为代谢能(千焦/尾/天)和G为生长能(千焦/尾/天)。 

1.1.4 排粪能与摄食能 

长吻鮠的排粪能F(千焦/尾/天)与摄食能C(千焦/尾/天)之间的关系为： 

F＝0.1261×C+0.0215，n＝15，R²＝0.9406。 

1.1.5排泄能与摄食能 

长吻鮠的排泄能U(千焦/尾/天)与摄食能C(千焦/尾/天)之间的关系为： 

U＝0.0704×C-0.0199，n＝15，R²＝0.9909。 

1.1.6 标准代谢和活动代谢与摄食水平 

长吻鮠的代谢能分为体增热SDA，标准代谢R_s和活动代谢R_a。长吻鮠标准代谢和活动代谢与摄食能的比值(R_a+R_s)/C与摄食水平RL(％体重/天)之间的关系为： 

(R_a+R_s)/C＝0.0558×RL²-0.2844×RL+0.8855，n＝15，R²＝0.7074； 

由此可得，最适摄食水平时长吻鮠的代谢能R_opt(千焦/尾/天)和最大摄食水平时的代谢能R(千焦/尾/天)的关系为： 

R_opt＝(SDA+R_a+R_s)×61％。 

1.2 体重和水温对长吻鮠的生长和摄食的影响 

17、36、67和303克四个体重组的长吻鮠，分别在20℃、24℃、28℃和32℃四组水温条件下进行4个月的生长实验。 

1.2.1 鱼体能值与体重和水温 

长吻鮠能值E_t(千焦/尾)和体重W(克)及水温T(℃)之间的关系为： 

lnE_t＝1.65091+1.02946×lnW-0.0014×T×lnW 

n＝60，R²＝0.9940，P＜0.01 

1.2.2 最大摄食量与体重和水温 

长吻鮠最大摄食率FI_max(克/尾/天)与体重W(克)和水温T(℃)之间的关系为： 

lnFI_max＝-14.34+1.26×lnW+0.72×T-0.012×T²-0.018×T×lnW 

n＝48，R²＝0.9225，P＜0.05 

对于体重为100克的长吻鮠，获得最大摄食量的水温是在26.5℃；而对于体重为400 克的长吻鮠，获得最大摄食量的水温则是在25.5℃。 

2.模型的建立 

在上述生长实验结果的基础上，建立长吻鮠的生物能量学模型，主要包括以下子模型：鱼体能值子模型、最大摄食能子模型、排粪子模型、排泄子模型、体增热子模型、标准代谢和活动代谢子模和日生长子模型；在生物能量学模型的基础上建立最适摄食量子模型。 

2.1 生物能量学模型 

2.1.1 鱼体能值子模型 

不同水温下长吻鮠的能值E_t随体重的增加而显著上升(P＜0.01)。长吻鮠能值E_t(千焦/尾)与体重W(克)和水温T(℃)之间的关系可以用下式方程进行回归模拟： 

lnE_t＝1.65091+1.02946×lnW-0.0014×T×lnW 

n＝60，R²＝0.9940，P＜0.01 

2.1.2 最大摄食能子模型 

长吻鮠的最大摄食能在不同水温下随体重的增加而增加。长吻鮠最大摄食量FI_max(克/尾/天)与体重W(克)、水温T(℃)之间的关系可以用以下方程模拟： 

lnFI_max＝-14.34+1.26×lnW+0.72×T-0.012×T²-0.018×T×lnW 

n＝48，R²＝0.9225，P＜0.05 

由于长吻鮠的摄食量F1(克/尾/天)与摄食水平RL(％体重/天)和体重W(克)之间存在下式关系：FI＝W×RL/100； 

由此，长吻鮠的最大摄食水平RL_max(％体重/天)与体重W(克)和水温T(℃)之间的关系： 

RL_max＝100×Exp(-14.34+1.26×lnW+0.72×T-0.012×T²-0.018×T×lnW)/W 

而最大摄食能C_max(千焦/尾/天)与最大摄食量FI_max(克/尾/天)和饲料能值E_d(千焦/克)的关系为：C_max＝FI_max×E_d； 

因此，长吻鮠的最大摄食能C_max与体重W(克)、水温T(℃)和饲料能值E_d之间的关系为： 

C_max＝(Exp(-14.34+1.26×lnW+0.72×T-0.012×T²-0.018×T×lnW))×E_d

2.13 排粪子模型 

在鱼类摄食不受限制的情况下，鱼类的排粪能占摄食能的比例恒定。从生长实验结果可得，长吻鮠的排粪能F(千焦/尾/天)和摄食能C(千焦/尾/天)之间的回归关系为： 

F＝0.1261×C+0.0215 

n＝15，R²＝0.9406，P＜0.05 

长吻鮠排粪能没有受到摄食水平的显著影响，最适摄食水平下长吻鮠的排粪能F_opt(千焦/尾/天)和摄食能C_opt(千焦/尾/天)之间的回归关系为： 

F_opt＝0.1261×C_opt+0.0215 

2.1.4排泄子模型 

在鱼类摄食不受限制的情况下，鱼类的排泄能占摄食能的比例恒定。从生长实验结果可得，长吻鮠的排泄能U(千焦/尾/天)和摄食能C(千焦/尾/天)之间的回归关系为： 

U＝0.0704×C-0.0199 

n＝15，R²＝0.9909，P＜0.05 

长吻鮠排泄能没有受到摄食水平的显著影响，最适摄食水平下长吻鮠的排泄能U_opt(千焦/尾/天)和摄食能C_opt(千焦/尾/天)之间的回归关系为： 

U_opt＝0.0704×C_opt-0.0199。 

2.1.5 体增热子模型 

体增热SDA和摄食率能C(千焦/尾/天)之间的关系为： 

SDA＝SDA_C×C/100 

SDA_C＝9.03+0.0502×P_d-0.0541W 

SDA＝(9.03+0.0502×P_d-0.0541W)×C/100 

式中，SDA_C为体增热系数；P_d为饲料的蛋白含量(％)；W为长吻鮠体重(g)。 

2.1.6 标准代谢和活动代谢子模型 

鱼类的代谢能分为标准代谢(R_s)、活动代谢(R_a)和体增热(SDA)。从生长实验结果得出，除去体增热SDA后，长吻鮠的标准代谢和活动代谢(R_s+R_a)占总摄食能C的比例与摄食水平RL(％体重/天)之间的回归关系为： 

(R_a+R_s)/C＝0.0558×RL²-0.2844×RL+0.8855 

n＝15，R²＝0.7074，P＜0.05 

2.1.7 日生长子模型 

长吻鮠摄食不受限制时的生长能G(千焦/尾/天)由能量收支式计算得到： 

G＝C-F-U-SDA-(R_a+R_s) 

长吻鮠的第t天的能值E_t(千焦/尾)是第t天的日生长能G_t和第(t-1)天的能值E_t-1之和：E_t＝E_t-1+G_t。 

长吻鮠的第t天的体重W_t(g)由下式方程得到： 

lnW_t＝(lnE_t-1.65091)/(1.02946-0.0014×T) 

n＝60，R²＝0.9940，P＜0.01。 

2.2最适摄食量模型 

最适摄食量定义为：不影响最大生长的情况下，长吻鮠摄食的最少的饲料量。长吻鮠的最适摄食能C_opt(千焦/尾/天)由能量收支式方程右边四个组成部分(G、R_opt、F_opt和U_opt)相加得到： 

C_opt＝G+R_opt+F_opt+U_opt； 

由上述可得，式中G、R_opt、F_opt和U_opt： 

G＝C-F-U-SDA-(R_a+R_s)； 

R_opt＝(SDA+R_a+R_s)×61％； 

F_opt＝0.1261×C_opt+0.0215； 

U_opt＝0.0704×C_opt-0.0199； 

由此得到的最适摄食能C_opt(千焦/尾/天)与长吻鮠最适摄食量FI_opt(克/尾/天)之间的关系为：C_opt＝FI_opt×E_d；式中E_d为饲料能量含量(千焦/克)。 

3.长吻鮠养殖动态投饲表 

根据上述模型的公式，采用计算机编程建立了长吻鮠的最适摄食量模型。计算时采用的步长为1天，流程为：(1)程式开始；(2)输入模拟次数NRUN和模拟天数NDAY(天)；(3)输入水温T(℃)、初始体重W_o(g)、鱼体初始能值E_o(千焦/尾)、饲料蛋白含量P_d(％)和饲料能值E_d(千焦/克)；(4)计算日最大摄食能C(千焦/尾/天)和最大摄食量FI_max(克/尾/天)；(5)计算日排粪能F(千焦/尾/天)和排泄能U(千焦/尾/天)；(6)计算日体增热SDA；(7)计算日标准代谢和活动代谢(R_s+R_a)；(8)根据公式G＝C-F-U-SDA-(R_a+R_s)，计算日能量生长G(千焦/尾)；(9)根据第t天的能值E_t(千焦/尾)、第t天的日生长能G_t和第(t-1)天的能值E_t-1的关系：E_t＝E_t-1+G_t和公式lnW_t＝(lnE_t-1.65091)/(1.02946-0.0014×T)计算日体重W_t(克)；(10)计算日最适代谢能R_opt(千焦/尾/天)；(11)计算日最适摄食能C_opt(千焦/尾/天)；(12)计算日最适摄食量FI_opt(克/尾/天)；(13)打印日最大摄食量FI_max、日最适摄食量FI_opt和日体重W_t；(14)打印整个养殖周期的体重、天数、累积最大摄食量、累积最适摄食量和累积能量收支；(15)程式结束。请见：表4长吻鮠养殖的动态投饲表，这种动态投饲表能够动态地、精准地预测长吻鮠养殖过程中不同环境下的日摄食量和日生长，不会造成饲料的浪费而增加饲料成本及饲料对养殖水体的次生污染或因饲料投喂不足而 影响鱼的生长。 

本发明与现有技术相比，具有以下优点：方法易行，操作简便，长吻鮠的生物能量学模型严格建立在实验数据的基础上，由此形成的长吻鮠动态投饲表，能够及时地、精准地预测养殖过程中不同水温下、不同生长阶段长吻鮠的日生长和日摄食量，从而提高饲料效率和养殖效益，降低饲料成本和养殖次生污染，对渔业管理和养殖环境的保护具有重要意义。 

附图说明

图1是特定生长率SGR与摄食水平RL之间的关系。 

图2是饲料转化效率FCE与摄食水平RL之间的关系。 

图3是长吻鮠的排粪能F与摄食能C之间的关系。 

图4是长吻鮠的排泄能U与摄食能C之间的关系。 

图5是长吻鮠的标准代谢和活动代谢与摄食能的比值(R_a+R_s)/C与摄食水平RL的关系。 

图6是不同水温下长吻鮠的能值E_t(千焦/尾)与体重W(克)之间的关系。 

其中：A：20℃长吻鮠的能值E_t与体重的关系；B：24℃长吻鮠的能值E_t与体重的关系；C：28℃长吻鮠的能值E_t与体重的关系；D：32℃长吻鮠的能值E_t与体重的关系。 

图7是不同水温下长吻鮠最大摄食量FI_max(克/尾/天)与体重W(克)之间的关系。 

其中：A：20℃长吻鮠最大摄食量FI_max与体重的关系；B：24℃长吻鮠最大摄食量FI_max与体重的关系；C：28℃长吻鮠最大摄食量FI_max与体重的关系；D：32℃长吻鮠最大摄食量FI_max与体重的关系。 

图8是长吻鮠最适摄食量模型的计算机程序流程图。 

图中1.程式开始，2.输入模拟次数(NRUN)和天数(NDAY)，3.输入水温T、初始体重W_o、初始能值E_o、饲料蛋白P_d和饲料能值E_d，4.计算最大摄食能C和日最大摄食量RL_max，5.计算排粪能F和排泄能U，6.计算体增热SDA，7.计算标准代谢和活动代谢(R_s+R_a)，8.计算日能量生长G，9.计算日体重生长W_t，10.计算最适总代谢R_opt，11.计算最适摄食能C_opt，12.计算日最适摄食量RL_opt，13.打印日最大摄食量RL_max、日最适摄食量RL_opt和终末体重W_t，14.打印天数、累积摄食率和累积能量收支，15.程式结束。 

图9长吻鮠养殖104天后模型预测的生长和实测值的关系。 

具体实施方式

下面结合附图和实施例对本发明进一步说明： 

一种长吻鮠养殖动态投饲表的建立方法，其步骤如下： 

1.生长实验的条件 

1.1摄食水平对长吻鮠摄食、生长和能量收支的影响 

研究不同摄食水平对长吻鮠摄食、饲料利用和能量收支的影响。实验共设计5组摄食水平，分别为：饥饿、初始体重的0.8、1.6、2.4％/天和饱食投喂组，每个摄食水平有三个重复。实验在15个鱼缸(水体积120升)构成的循环水养殖系统中进行，水中溶氧高于7.5毫克/升，水温维持在29±2℃左右，pH值在7.7-8.0，氨氮浓度低于0.1毫克/升。 

实验开始前，长吻鮠饥饿一天。每个实验缸中随机转入12尾初始体重6.6g的幼鱼。实验期间，长吻鮠每天用实验饲料按照设计的摄食水平投喂两次(9:00，15:00)。每两周称鱼一次，以调整各组摄食率。实验饲料配方见表1，饲料中含有45％的粗蛋白和7.5％的粗脂肪。饲料颗粒直径大小为2mm左右，在4℃下保存。实验持续8周。 

表1 实验饲料配方 

所述的鱼粉为美国海鲜公司的白鱼粉，购自武汉高龙饲料有限公司； 

所述的鱼油为秘鲁鯷鱼油，购自武汉高龙饲料有限公司； 

所述的矿物盐预混物为(克/千克饲料)：氯化钠NaCl，1；硫酸镁MgSO₄.7H₂O，15；磷酸二氢钾KH₂PO₄，32；磷酸二氢钙Ca(H₂PO₄)₂，20；柠檬酸铁FeC₆H₅O₇.5H₂O，2.5；乳酸钙C₆H₁₀CaO₆.5H₂O，3.5；硫酸锌ZnSO₄.7H₂O，0.353；硫酸锰MnSO₄.4H₂O，0.162；硫酸铜CuSO₄.5H₂O，0.031；氯化钴CoCl.6H₂O，0.001；碘酸钾KIO₃，0.003；微晶纤维素，0.45； 

所述的维生素预混物为(克/千克饲料)：维生素A，110；维生素B₁，20；维生素B₂，20；维生素B₆，20；维生素B₁₂，0.020；维生素C，100；维生素D，20；维生素E，50；维生素K，10；叶酸，5；泛酸钙，50；肌醇，100；烟酸，100；生物素，0.1；面粉，1645.2； 

所述的纤维素为微晶纤维素，购自山东聊城阿华制药有限公司。 

1.1.1特定生长率与摄食水平 

长吻鮠的生长用特定生长率表示：特定生长率(SGR，％/天)＝100×(lnW_t-lnW_o)/t，式中W_t(克)和W_o(克)分别为终末和初始体重，t(天)为实验持续时间。特定生长率SGR与摄食水平RL(％体重/天)之间的关系如图1所示，是减速增长的曲线模式： 

SGR＝-0.66+2.65×ln(RL+1)，n＝15，R²＝0.9998。 

1.1.2  饲料转化效率与摄食水平 

饲料利用用饲料转化效率表示，饲料转化效率(FCE，％)＝100×(W_t-W_o)/FI，式中W_t(克)和W_o(克)分别为终末和初始体重，摄食量FI(克)为长吻鮠摄食饲料的重量。如图2所示，饲料转化效率FCE与摄食水平RL(％体重/天)之间的关系呈二次曲线： 

FCE＝97.1+69.3×RL-28.9×RL²，n＝15，R²＝0.675； 

在饲料转化效率最大时，即长吻鮠摄食最少的饲料量，得到最佳的生长，定义为最适摄食水平。长吻鮠的最适摄食水平为初始体重的1.2％体重/天。 

1.1.3  能量收支与摄食水平 

摄食水平对长吻鮠能量收支有显著的影响(表2)。摄食水平为0.8％体重/天的长吻鮠排泄能占摄食能(U/C)的比例显著低于其他摄食水平组(P＜0.05)。摄食水平为2.4％体重/天的长吻鮠代谢能占摄食能(R/C)的比例显著低于其他摄食水平组。食物能中分配给生长能的部分随摄食水平而上升。 

饱食水平下，长吻鮠的能量收支式为：100C＝12.9×F+5.5×U+40.1×R+41.5×G，式中C为摄食能(千焦/尾/天)，F为粪便能(千焦/尾/天)，U为排泄能(千焦/尾/天)，R为代谢能(千焦/尾/天)和G为生长能(千焦/尾/天)。 

表2  不同摄食水平下长吻鮠的能量收支式 

*表中各数值(平均数±标准误)后的不同字母表示差异显著(P＜0.05)。 

1.1.4 排粪能与摄食能 

随着摄食水平的上升，粪便能占摄入饲料总能(F/C)的比例没有显著差异，大约在13-14％。如图3所示，长吻鮠的排粪能F(千焦/尾/天)与摄食能C(千焦/尾/天)之间的关系为直线回归的关系： 

F＝0.1261×C+0.0215，n＝15，R²＝0.9406。 

1.1.5排泄能与摄食能 

随着摄食水平的上升，排泄能占摄入饲料总能(U/C)的比例没有显著差异，大约在4-5％。如图4所示，长吻鮠的排泄能U(千焦/尾/天)与摄食能C(千焦/尾/天)之间的关系为直线回归的关系： 

U＝0.0704×C-0.0199，n＝15，R²＝0.9909。 

1.1.6 标准代谢和活动代谢与摄食水平 

长吻鮠的代谢能分为体增热SDA，标准代谢R_s和活动代谢R_a。如图5所示，长吻鮠标准代谢和活动代谢与摄食能的比值(R_a+R_s)/C与摄食水平RL(％体重/天)之间呈显著的二次曲线关系： 

(R_a+R_s)/C＝0.0558×RL²-0.2844×RL+0.8855，n＝15，R²＝0.7074； 

由此可得，摄食水平对长吻鮠的代谢能有显著的影响，在最适摄食水平时长吻鮠的代谢能R_opt(千焦/尾/天)和最大摄食水平时的代谢能R(千焦/尾/天)的关系为： 

R_opt＝(SDA+R_a+R_s)×61％。 

1.2 体重和水温对长吻鮠的生长和摄食的影响 

研究长吻鮠最大摄食率、鱼体能值和体重及水温之间的关系。实验采用17、36、67和303克四个体重组，每个体重组有三个重复，分四个实验分别在20℃、24℃、28℃和32℃四组水温条件下进行。每组水温实验进行1个月，共持续4个月。实验在室内循环水养殖系统的12个养殖缸(直径1.5米，容积1200升)中进行。实验开始前，各实验缸水温通过加热棒和温控仪来加热和调节，以2-3℃/天升到实验温度并维持水温的恒定。 

实验开始前，长吻鮠饥饿一天。为了保证各处理组的生物量相似，每组实验缸中随机转入不同尾数的长吻鮠，四个体重组分别转入200尾(17克)、100尾(36克)、50尾(67 克)和10尾(303克)的长吻鮠。实验期间，每天用实验饲料饱食投喂两次(9:00，15:00)。实验饲料配方见表3，实验饲料根据实验鱼大小制成大中小三种颗粒，最小的颗粒直径大小为2mm左右，保存在4℃下直至实验完成。 

表3 实验饲料配方 

1.2.1 鱼体能值与体重和水温 

长吻鮠能值E_t(千焦/尾)和体重及水温之间的关系如图6所示。不同水温下长吻鮠的能值E_t随体重的增加而显著上升(P＜0.01)。多元回归分析显示体重W(克)和水温T(℃)与长吻鮠能值E_t之间的关系可以用下式方程进行回归模拟： 

lnE_t＝1.65091+1.02946×lnW-0.0014×T×lnW 

n＝60，R²＝0.9940，P＜0.01 

通过逐步回归分析表明，体重和水温对长吻鮠鱼体能值E_t的影响存在显著的交互作用(P＜0.01)。 

1.2.2 最大摄食量与体重和水温 

不同水温下，长吻鮠最大摄食量随体重和水温的变化如图7所示。长吻鮠的最大摄食量在不同水温下随体重的增加而增加。多元回归分析显示体重W(克)、水温T(℃)和长吻鮠最大摄食量FI_max(克/尾/天)之间的关系可以用以下方程模拟： 

lnFI_max＝-14.34+1.26×lnW+0.72×T-0.012×T²-0.018×T×lnW 

n＝48，R²＝0.9225，P＜0.05 

通过逐步回归分析表明，体重和水温在对长吻鮠最大摄食量的影响存在显著的交互作用(P＜0.05)。对于体重为100克的长吻鮠，获得最大摄食量的水温是在26.5℃；而对于体 重为400克的长吻鮠，获得最大摄食量的水温则是在25.5℃。 

2.模型的建立 

在上述生长实验结果的基础上，建立长吻鮠的生物能量学模型，主要包括以下子模型：鱼体能值子模型、最大摄食能子模型、排粪子模型、排泄子模型、体增热子模型、标准代谢和活动代谢子模和日生长子模型；在生物能量学模型的基础上建立最适摄食量子模型。 

2.1 生物能量学模型 

2.1.1 鱼体能值子模型 

不同水温下长吻鮠的能值E_t随体重的增加而显著上升(P＜0.01)。多元回归分析显示长吻鮠能值E_t(千焦/尾)与体重W(克)和水温T(℃)之间的关系可以用下式方程进行回归模拟： 

lnE_t＝1.65091+1.02946×lnW-0.0014×T×lnW 

n＝60，R²＝0.9940，P＜0.01 

通过逐步回归分析表明，体重和水温对长吻鮠鱼体能值E_t的影响存在显著的交互作用(P＜0.01)。 

2.1.2 最大摄食能子模型 

长吻鮠的最大摄食量在不同水温下随体重的增加而增加。多元回归分析表明长吻鮠最大摄食量FI_max(克/尾/天)与体重W(克)、水温T(℃)之间的关系可以用以下方程模拟： 

lnFI_max＝-14.34+1.26×lnW+0.72×T-0.012×T²-0.018×T×lnW 

n＝48，R²＝0.9225，P＜0.05 

由于长吻鮠的摄食量FI(克/尾/天)与摄食水平RL(％体重/天)和体重W(克)之间存在下式关系：FI＝W×RL/100； 

由此，长吻鮠的最大摄食水平RL_max与体重W(克)和水温T(℃)之间的关系为： 

RL_max＝100×Exp(-14.34+1.26×lnW+0.72×T-0.012×T²-0.018×T×lnW)/W。 

而最大摄食能C_max(千焦/尾/天)与最大摄食量FI_max(克/尾/天)和饲料能值E_d(千焦/克)的关系为：C_max＝FI_max×E_d； 

因此，长吻鮠的最大摄食能C_max与体重W(克)、水温T(℃)和饲料能值E_d之间的关系为： 

C_max＝(Exp(-14.34+1.26×lnW+0.72×T-0.012×T²-0.018×T×lnW))×E_d。 

2.1.3排粪子模型 

在鱼类摄食不受限制的情况下，鱼类的排粪能占摄食能的比例恒定。从生长实验结果 可得，长吻鮠的排粪能F(千焦/尾/天)和摄食能C(千焦/尾/天)之间的回归关系为： 

F＝0.1261×C+0.0215 

n＝15，R²＝0.9406，P＜0.05 

从表2可以看出，不同摄食水平下长吻鮠排粪能在总摄食能中所占的比例大约在13-14％，没有受到摄食水平的显著影响。因此，最适摄食水平下长吻鮠的排粪能F_opt(千焦/尾/天)和最适摄食能C_opt(千焦/尾/天)之间的回归关系为： 

F_opt＝0.1261×C_opt+0.0215 

2.1.4 排泄子模型 

在鱼类摄食不受限制的情况下，鱼类的排泄能占摄食能的比例恒定。从生长实验结果可得，长吻鮠的排泄能U(千焦/尾/天)和摄食能C(千焦/尾/天)之间的回归关系为： 

U＝0.0704×C-0.0199 

n＝15，R²＝0.9909，P＜0.05 

从表2可以看出，不同摄食水平下长吻鮠排泄能在总摄食能中所占的比例大约在4-5％，没有受到摄食水平的显著影响。因此，最适摄食水平下长吻鮠的排泄能U_opt(千焦/尾/天)和最适摄食能C_opt(千焦/尾/天)之间的回归关系为： 

U_opt＝0.0704×C_opt-0.0199。 

2.1.5  体增热子模型 

有学者引入了体增热系数(SDA_C)来描述体增热SDA和摄食能C之间的关系(Ross等The effects of dietary protein content，lipid content and ration level on oxygen consumption and specific dynamic action in Oreochromis niloticus.Comparative Biochemistry and Physiology，1992：103A，573-578)。体增热SDA和摄食率能C(千焦/尾/天)之间的关系为： 

SDA＝SDA_C×C/100 

SDA_C＝9.03+0.0502×P_d-0.0541W 

SDA＝(9.03+0.0502×P_d-0.0541W)×C/100 

式中，SDA_C为体增热系数；P_d为饲料的蛋白含量(％)；W为长吻鮠体重(g)。 

2.1.6 标准代谢和活动代谢子模型 

鱼类的代谢能分为标准代谢(R_s)、活动代谢(R_a)和体增热(SDA)。从生长实验结果得出，除去体增热SDA后，长吻鮠的标准代谢和活动代谢(R_s+R_a)占总摄食能C的比例与摄食水平RL(％体重/天)之间的回归关系为： 

(R_a+R_s)/C＝0.0558×RL²-0.2844×RL+0.8855 

n＝15，R²＝0.7074，P＜0.05 

2.1.7 日生长子模型 

长吻鮠摄食不受限制时的生长能G(千焦/尾/天)由能量收支式计算得到： 

G＝C-F-U-SDA-(R_a+R_s) 

长吻鮠的第t天的能值E_t(千焦/尾)是第t天的日生长能G_t和第(t-1)天的能值E_t-1之和：E_t＝E_t-1+G_t。 

长吻鮠的第t天的体重W_t(g)由下式方程得到： 

lnW_t＝(lnE_t-1.65091)/(1.02946-0.0014×T) 

n＝60，R²＝0.9940，P＜0.01。 

2.2 最适摄食量模型 

最适摄食量定义为：不影响最大生长的情况下，长吻鮠摄食的最少的饲料量。长吻鮠的最适摄食能C_opt(千焦/尾/天)由能量收支式方程右边四个组成部分(G、R_opt、F_opt和U_opt)相加得到： 

C_opt＝G+R_opt+F_opt+U_opt； 

由上述可得，式中G、R_opt、F_opt和U_opt： 

G＝C-F-U-SDA-(R_a+R_s)； 

R_opt＝(SDA+R_a+R_s)×61％； 

F_opt＝0.1261×C_opt+0.0215； 

U_opt＝0.0704×C_opt-0.0199； 

由此得到的最适摄食能C_opt(千焦/尾/天)与长吻鮠最适摄食量FI_opt(克/尾/天)之间的关系为：C_opt＝FI_opt×E_d；式中E_d为饲料能量含量(千焦/克)。 

3.长吻鮠养殖动态投饲表 

根据上述模型的公式，采用计算机编程建立了长吻鮠的最适摄食率模型(图8)。模型需要输入的初始值为：养殖天数、水温、鱼初始体重、鱼体初始能值、饲料蛋白水平和饲料能值，经过计算，输出长吻鮠养殖的最终体重、最大摄食率、最适摄食率和能量收支。计算时采用的步长为1天，基本流程为：1.程式开始，2.输入模拟次数NRUN、模拟天数NDAY(天)，3.输入水温T(℃)、初始体重W_o(g)、鱼体初始能值E_o(千焦/尾)、饲料蛋白水平P_d(％)和饲料能值E_d(千焦/克)，4.计算最大摄食能C(千焦/尾/天)和日最大摄食量FI_max(克/尾/天)，5.计算排粪能F(千焦/尾/天)和排泄能U(千焦/尾/天)，6.计算体增热SDA，7.计算标准代谢和活动代谢(R_s+R_a)，8.根据公式G＝C-F-U-SDA-(R_a+R_s)，计算日能 量生长G(千焦/尾)，9.根据第t天的能值E_t(千焦/尾)、第t天的日生长能G_t和第(t-1)天的能值E_t-1的关系：E_t＝E_t-1+G_t和公式lnW_t＝(lnE_t-1.65091)/(1.02946-0.0014×T)计算日体重W_t(克)，10.计算最适总代谢R_opt(千焦/尾/天)，11.计算最适摄食能C_opt(千焦/尾/天)，12.计算日最适摄食量RL_opt(克/尾)，13.打印日最大摄食量RL_max、日最适摄食量RL_opt和终末体重W_t，14.打印天数、累积摄食率和累积能量收支，15.程式结束。 

4.验证实验 

长吻鮠动态投饲表的建立就是用于指导长吻鮠养殖。在实验数据基础上建立的模型对生产实践具有多大的指导意义，需要验证实验结果的证实。为此，设计了水泥池养殖的验证实验，探讨了三种投饲策略(T_A组：最大投饲水平组，T_B组：生物能量学模型预测的最适投饲水平组，T_C组：最适投饲水平80％的处理组)对长吻鮠生长、摄食和氮磷排放的影响。每种策略作为一个处理，每个处理三个重复。 

实验用鱼长吻鮠购自湖北石首市国家级长吻鮠良种场。验证实验在9个水泥池(7×4×0.6m，容量16.8m³)中进行。每周用曝气过的自来水交换30％的池水。实验开始前，长吻鮠饥饿1天，随即选择平均体重50g左右的长吻鮠投放到9个水泥池中(65尾/池)。如表4所示，实验期间每天按照长吻鮠动态投饲表，用商品饲料饱食投喂两次(06:00和18:00)。专用商品饲料购自湖北省武汉市高龙饲料公司，饲料粗蛋白含量为42.8％，粗脂肪含量为4％，饲料能值为18.4千焦/克。 

实验期间，水温为21.5℃～29.9℃，氨氮浓度低于0.5毫克/升，溶氧大于7.8毫克/升，pH平均为8.1，自然光照。实验持续104天。实验结束时，长吻鮠饥饿1天全池称重。 

表4 长吻鮠养殖的动态投饲表 

* FI_max为最适摄食量模型预测的日最大摄食量，FI_opt为模型预测的日最适摄食量；养殖周期为104天，初始体重为50克，水温为21.5℃～29.9℃，鱼体初始能值为251.8千焦/尾，商品饲料蛋白水平为42.8％，饲料能值为18.4千焦/克，饲料干物质含量93.8％。 

如表5所示，经过104天的养殖，最大摄食水平下长吻鮠T_A组的终末体重(193.6±13.0克)和最适摄食水平下T_B组的体重(175.0±6.5克)没有显著差异，这和最适摄食量模型预测的终末体重(178.1克)基本一致(图9)；T_A组的摄食率显著高于T_B组，而T_A组的饲料转化效率低于T_B组的饲料转化效率。 

表5 投喂策略对长吻鮠生长和饲料利用的影响 

*表中各数值(平均数±标准误)后的不同字母表示差异显著(P＜0.05)；T_A组：最大投饲水平组，T_B组：生物能量学模型预测的最适投饲水平组，T_C组：最适投饲水平80％的处理组。 

验证实验结果表明，最适摄食量模型能够精准的预测不同环境下的长吻鮠日体重生长和日摄食量，由此形成的长吻鮠动态投饲表在长吻鮠养殖中具有广泛的应用前景。 

Claims (1)

1.一种长吻鮠养殖动态投饲表的建立方法，其步骤如下：

A、生长实验的条件

1.摄食水平对长吻鮠摄食、生长和能量收支，经过8周的生长，体重为6.6克的长吻鮠幼鱼在饥饿及初始体重的0.8％/天、1.6％/天、2.4％/天和饱食摄食水平下的摄食、生长和能量收支；

1.1 生长率与摄食水平，长吻鮠的生长用生长率表示：生长率SGR＝100×(lnW_t-lnW_o)/t，式中W_t和W_o分别为终末和初始体重，t为实验持续时间，SGR的单位为％/天，W_t和W_o单位为克，t的单位为天；生长率SGR与摄食水平RL之间的关系为：

SGR＝-0.66+2.65×ln(RL+1)，n＝15，R²＝0.9998；

SGR的单位为％/天，RL的单位为％体重/天；

1.2 饲料转化效率与摄食水平，饲料利用用饲料转化效率表示，饲料转化效率＝100×(W_t-W_o)/FI，式中W_t和W_o分别为终末和初始体重，摄食量FI为长吻鮠摄食饲料的重量，W_t和W_o单位为克，FI的单位为克/尾/天；饲料转化效率FCE与摄食水平RL之间的关系为：

FCE＝97.1+69.3×RL-28.9×RL²，n＝15，R²＝0.675；

FCE的单位为％，RL的单位为％体重/天；

能量收支与摄食水平，饱食水平下，长吻鮠的能量收支式为：100 C＝12.9×F+5.5×U+40.1×R+41.5×G，式中C为摄食能，F为排粪能，U为排泄能，R为代谢能，G为生长能；

C、F、U、R和G的单位为千焦/尾/天；

1.3排粪能与摄食能，长吻鮠的排粪能F与摄食能C之间的关系为：

F＝0.1261×C+0.0215，n＝15，R²＝0.9406；

F的单位为千焦/尾/天，C的单位为千焦/尾/天；

1.4排泄能与摄食能，长吻鮠的排泄能U与摄食能C之间的关系为：

U＝0.0704×C-0.0199，n＝15，R²＝0.9909；

U的单位为千焦/尾/天，C的单位为千焦/尾/天；

1.5 标准代谢和活动代谢与摄食水平，长吻鮠的代谢能分为体增热SDA，标准代谢R_s和活动代谢R_a，长吻鮠标准代谢和活动代谢与摄食能C的比值与摄食水平RL之间的关系为：

$\frac{\mathrm{Ra}+\mathrm{Rs}}{C}=0.0558\×{\mathrm{RL}}^2-0.2844\×\mathrm{RL}+0.8855,$ n＝15，R²＝0.7074；

C的单位为千焦/尾/天，RL的单位为％体重/天；

由此，在最适摄食水平时长吻鮠的代谢能R_opt和最大摄食水平时的代谢能R的关系为：

R_opt＝(SDA+R_a+R_s)×61％；

R_opt、R和SDA的单位为千焦/尾/天；

2.体重和水温对长吻鮠的生长和摄食的影响：

17、36、67和303克四个体重组的长吻鮠，分别在20℃、24℃、28℃和32℃四组水温条件下进行4个月的生长实验：

2.1鱼体能值与体重和水温，长吻鮠能值E_t和体重W及水温T之间的关系为：

lnE_t＝1.65091+1.02946×lnW-0.0014×T×lnW

n＝60，R²＝0.9940，P＜0.01；

E_t的单位为千焦/尾，W的单位为克，T的单位为℃；

2.2最大摄食量与体重和水温，长吻鮠最大摄食量FI_max与体重W和水温T之间的关系为：

lnFI_max＝-14.34+1.26×lnW+0.72×T-0.012×T²-0.018×T×lnW

n＝48，R²＝0.9225，P＜0.05；

FI_max的单位为克/尾/天，W的单位为克，T的单位为℃；

B、模型的建立

1.生物能量学模型

1.1鱼体能值子模型，不同水温下长吻鮠的能值E_t随体重的增加上升，长吻鮠能值E_t与体重W和水温T之间的关系用下式方程进行回归模拟：

lnE_t＝1.65091+1.02946×lnW-0.0014×T×lnW

n＝60，R²＝0.9940，P＜0.01；

E_t的单位为千焦/尾，W的单位为克，T的单位为℃；

1.2最大摄食能子模型，长吻鮠的最大摄食量在不同水温下随体重的增加，长吻鮠最大摄食量FI_max与体重W、水温T之间的关系用以下方程模拟：

lnFI_max＝-14.34+1.26×lnW+0.72×T-0.012×T²-0.018×T×lnW

n＝48，R²＝0.9225，P＜0.05；

FI_max的单位为克/尾/天，W的单位为克，T的单位为℃；

长吻鮠的摄食量FI(克/尾/天)与摄食水平RL(％体重/天)和每尾鱼体重W(克)之间存在下式关系：FI＝W×RL/100；由此，长吻鮠的最大摄食水平RL_max与体重W、水温T之间的关系为：

RL_max＝100×Exp(-14.34+1.26×lnW+0.72×T-0.012×T²-0.018×T×lnW)/W；

RL_max的单位为％体重/天，W的单位为克，T的单位为℃；

最大摄食能C_max(千焦/尾/天)与最大摄食量FI_max(克/尾/天)和饲料能值E_d(千焦/克)的关系为：C_max＝FI_max×E_d；由此，长吻鮠的最大摄食能C_max与体重W、水温T和饲料能值E_d之间的关系为：

C_max＝(Exp(-14.34+1.26×lnW+0.72×T-0.012×T²-0.018×T×lnW))×E_d；

C_max的单位为千焦/尾/天，W的单位为克，T的单位为℃，E_d的单位为千焦/克；

1.3排粪子模型，长吻鮠的排粪能F和摄食能C之间的回归关系为：

F＝0.1261×C+0.0215

n＝15，R²＝0.9406，P＜0.05；

F的单位为千焦/尾/天，C的单位为千焦/尾/天；

长吻鮠排粪能不受摄食水平的影响，最适摄食水平下长吻鮠的排粪能F_opt和摄食能C_opt之间的关系为：

F_opt＝0.1261×C_opt+0.0215；

F_opt的单位为千焦/尾/天，C_opt的单位为千焦/尾/天；

1.4 排泄子模型，长吻鮠的排泄能U和摄食能C之间的回归关系为：

U＝0.0704×C-0.0199

n＝15，R²＝0.9909，P＜0.05；

U的单位为千焦/尾/天，C的单位为千焦/尾/天；

长吻鮠排泄能不受摄食水平的影响，最适摄食水平下长吻鮠的排泄能U_opt和摄食能C_opt之间的关系为：

U_opt＝0.0704×C_opt-0.0199；

U_opt的单位为千焦/尾/天，C_opt的单位为千焦/尾/天；

1.5 体增热子模型，体增热SDA和摄食能C之间的关系为：

SDA＝SDA_C×C/100

SDA_C＝9.03+0.0502×P_d-0.0541W

SDA＝(9.03+0.0502×P_d-0.0541W)×C/100；

式中，SDA_C为体增热系数；P_d为饲料的蛋白含量％；W为长吻鮠体重g；SDA和C的单位为千焦/尾/天，W的单位为克，P_d的单位为％；

1.6标准代谢和活动代谢子模型，鱼类的代谢能分为标准代谢R_s、活动代谢R_a和体增热SDA，长吻鮠的标准代谢和活动代谢占总摄食能C的比例与摄食水平RL之间的关系为：

$\frac{\mathrm{Ra}+\mathrm{Rs}}{C}=0.0558\×{\mathrm{RL}}^2-0.2844\×\mathrm{RL}+0.8855$

n＝15，R²＝0.7074，P＜0.05；

C的单位为千焦/尾/天，RL的单位为％体重/天；

1.7日生长子模型，能量收支式计算得到长吻鮠的生长能G：

G＝C-F-U-SDA-(R_a+R_s)；

G、C、F、U和SDA的单位为千焦/尾/天；

长吻鮠的第t天的能值E_t是第t天的日生长能G_t和第t-1天的能值E_t-1之和：E_t＝E_t-1+G_t，

长吻鮠的第t天的体重W_t由下式方程得到：

lnW_t＝(lnE_t-1.65091)/(1.02946-0.0014×T)

n＝60，R²＝0.9940，P＜0.01；

W_t的单位为克，E_t的单位为千焦/尾，T的单位为℃；

2.最适摄食量模型

最适摄食量FI_opt定义为：不影响生长的情况下，长吻鮠摄食的最少的饲料量，能量收支式右边四个组成部分G、R_opt、F_opt和U_opt相加得到长吻鮠的最适摄食能C_opt：

C_opt＝G+R_opt+F_opt+U_opt；

式中G、R_opt、F_opt和U_opt：

G＝C-F-U-SDA-(R_a+R_s)；

R_opt＝(SDA+R_a+R_s)×61％；

F_opt＝0.1261×C_opt+0.0215；

U_opt＝0.0704×C_opt-0.0199；

G、C、F、U、SDA、R_opt、F_opt、U_opt和C_opt的单位为千焦/尾/天；

长吻鮠最适摄食能C_opt与最适摄食量FI_opt之间的关系为：C_opt＝FI_opt×E_d；式中E_d为饲料能量含量；

C_opt的单位为千焦/尾/天，FI_opt的单位为克/尾/天，E_d的单位为千焦/克；

C、长吻鮠养殖动态投饲表

根据上述模型的公式，采用计算机编程建立了长吻鮠的最适摄食量模型，计算时采用的步长为1天，流程为：(1)程式开始；(2)输入模拟次数NRUN和模拟天数NDAY(天)；(3)输入水温T(℃)、初始体重W_o(g)、鱼体初始能值E_o(千焦/尾)、饲料蛋白含量P_d(％)和饲料能值E_d(千焦/克)；(4)计算日最大摄食能C(千焦/尾/天)和最大摄食量FI_max(克/尾/天)；(5)计算日排粪能F(千焦/尾/天)和排泄能U(千焦/尾/天)；(6)计算日体增热SDA；(7)计算日标准代谢和活动代谢(R_s+R_a)；(8)根据公式G＝C-F-U-SDA-(R_a+R_s)，计算日能量生长G(千焦/尾)；(9)根据第t天的能值E_t(千焦/尾)、第t天的日生长能G_t和第(t-1)天的能值E_t-1的关系：E_t＝E_t-1+G_t和公式lnW_t＝(lnE_t-1.65091)/(1.02946-0.0014×T)计算日体重W_t(克)；(10)计算日最适代谢能R_opt(千焦/尾/天)；(11)计算日最适摄食能C_opt(千焦/尾/天)；(12)计算日最适摄食量FI_opt(克/尾/天)；(13)打印日最大摄食量FI_max、日最适摄食量FI_opt和日体重W_t；(14)打印整个养殖周期的体重、天数、累积最大摄食量、累积最适摄食量和累积能量收支；(15)程式结束。

Images