我国从20世纪80年代初就已经开始了近海资源增殖试验和较大规模的生产性种苗放流及海珍品的底播增殖，取得了很好的经济效益。近年来，由于过度捕捞和环境污染，我国近海生态环境受到严重破坏，渔业资源日渐枯竭，优质鱼比例越来越小，渔业资源结构极不合理，广大渔民收入降低。在目前渔业资源普遍衰退的情况下，加强资源增殖放流以恢复渔业资源、提高渔业产量及质量是我国渔业发展的大趋势，资源增殖放流对于我国顺利推行渔民转产转业的措施有重大意义。

但是，我国有关资源增殖放流的基础研究工作明显滞后，没有形成有效的增殖放流研究应用技术，以往开展的规模生产性种苗放流缺乏科学指导，在某种意义上带有很大的盲目性，对回捕率的悬殊变化无法作出科学解释^[1]。

本文提出一套“放流效果统计量评估法”的增殖放流效果评估新方法。该方法已应用于2005年广东省渔业资源增殖放流效果跟踪监测评估和标志放流研究。

1.   材料与方法

依据渔业资源评估原理^[2-8]，结合渔业资源增殖放流的特点，选用适用的渔业资源评估模型，估算未建立生长方程的增殖放流种类的生长参数及其自然死亡系数，并估算增殖放流种类的合理放流数量。提出确认渔获物中来自放流种苗数量的方法。推导计算捕捞死亡系数的公式，以及按时间序列计算放流群体残存量、回捕量、回捕率和回捕效益等的公式。

2.   结果

2.1   利用年龄-体长资料估算增殖放流种类的生长参数

对于未建立生长方程的增殖放流种类，收集其年龄-体长资料，应用GULLAND和HOLT图解法^[2-4]估算出生长参数l_∞和K：

可得：$K=-b, l_{\infty}=-a / b$

因已知不同年龄的体长，用下式估算各年龄的t₀值，然后取各年龄所求t₀的平均值。

式中：l_t为时间t时的体长，△l_t为时间差△t内的体长差，a、b为回归系数，l_∞为渐近体长，K为生长曲线的平均曲率。

2.2   自然死亡系数M的估算

对于没有研究过自然死亡系数的增殖放流种类，用PAULY^[6]的经验公式估算自然死亡系数。PAULY根据175种鱼类的有关数据和资料，用多元线性回归方法建立了自然死亡系数M的估算式：

式中：l_∞为渐近体长(cm)，K为生长曲线的平均曲率，T为放流种类栖息环境的年平均水温。

2.3   合理放流数量的估算

放流数量与渔获量不完全成正比关系，放流量超过一定数量后, 因栖息空间变小和自身代谢产物污染生活环境以及饵料生物密度可能降低等影响了放流生物的生长和生存，产量可能下降。

GULLAND^[7]认为，渔业资源的最适开发率为E=0.5，即捕捞死亡系数F等于自然死亡系数M。若认为以往某一时期(如20世纪60年代)的渔业资源处于没有捕捞过度状态及已调查过资源密度，利用当时的调查资料估计当时的平均大量捕捞年龄，以目前放流时种苗的日龄作为当时的补充日龄，在E=0.5的前提下，用残存率逆算当时补充日龄的资源量，可认为该资源量是目前的合理放流数量。

自然海区中的渔业资源种类由于自然死亡和捕捞死亡，数量不断减少，残存率^[2-4]为：

式中：S_t为t龄时的残存率，t_r为补充时的年龄或开始计算残存率的年龄，M为自然死亡系数，F为捕捞死亡系数。

2.4   确认渔获物中来自放流种苗的数量

渔获物中的标志放流种类，由于带有标志，容易确认。

对于没有标志的大规模增殖放流：研究增殖种类放流前和放流后的资源密度情况^[8]，依据放流种类的年龄生长规律、天然繁殖习性和人工培育种苗、放流种苗等资料，估算时间序列的放流种苗和天然种苗在海区中的同时生长情况，区分两者，或区分两者的比例，从调查渔获物中确认来自放流种苗的数量。

2.4.1   放流种类早期幼鱼的体长与日龄关系

式中：l_D为D日龄时幼鱼的体长，a、b为系数。

2.4.2   放流种类的体长体重关系

式中：l_t、W_t分别为t龄时的体长和体重，a为生长的条件因子，b为幂指数系数。

2.4.3   不同时期的生长情况用von Bertalanffy生长方程估算

式中：W_∞、l_∞分别为渐近体重和渐近体长，K为生长曲线的平均曲率，t为年龄，t₀为理论上体长和体重等于零时的年龄，b为幂指数系数。

2.5   放流后残存尾数的估算

估算增殖放流种类性成熟时的残存尾数，可以根据放流种类的个体繁殖力估算繁殖下一代的总量，从而评估其产生的呈良性循环的生态效益。

若某一增殖放流品种的放流数量为R，此时年龄为t_r，假设放流后到开捕年龄t_C之前放流品种没有受到捕捞，只有自然死亡，则t_C龄时海区中残存尾数为^[2-5]：

实际上，放流后从开捕年龄t_C开始，既有自然死亡也有捕捞死亡，总死亡系数Z为自然死亡系数M和捕捞死亡系数F之和。在自然死亡和捕捞死亡的共同作用下，到t龄时海区中残存尾数为：

2.6   捕捞死亡系数F的估算

由巴拉诺夫渔获量方程^[2-4]，放流后从开捕年龄t_C到年龄t时的总回捕尾数为：

则开捕以后任一时间间隔(t₁~t₂)的回捕尾数可以用下式估算：

其中：N_tC为放流后开捕年龄t_C时的残存尾数，M为自然死亡系数，F为捕捞死亡系数。

利用标志放流回捕资料、放流前本底调查和放流后跟踪调查资料(市场调查及渔民座谈、渔业资源调查)、渔捞日志登记资料及渔业生产统计资料，均可以使用上式迭代运算得出对应资料的F值。

2.7   放流后不同时期回捕尾数和回捕率的估算

回捕率是回捕尾数与放流苗数的百分比。增殖放流后，在自然死亡和捕捞死亡的共同作用下尾数不断减少，回捕率也在不断变化。

由巴拉诺夫渔获量方程^[2-4]，放流后从开捕年龄t_C到年龄t时的总回捕尾数为：

放流后从开捕年龄t_C到年龄t时的总回捕率为：

从开捕年龄t_C以后，第D天1天内的回捕尾数为：

其中：N_tC为放流后开捕年龄t_C时的残存尾数，M为自然死亡系数，F为捕捞死亡系数，t_D为放流后第D天时的年龄，t_(D+1)为放流后第D+1天时的年龄。

2.8   放流后不同时期回捕生物量和回捕效益估算

增殖放流后，在自然死亡和捕捞死亡的共同作用下尾数不断减少，又由于个体重量的生长使总生物量增加，两者共同作用的结果，使增殖放流种类的生物量随着年龄的增加而不断变化。

将巴拉诺夫渔获量方程^[2-4]和von Bertalanffy体重生长方程^[2-4]相结合推导，从开捕年龄t_C以后，第D天1天内的回捕生物量为：

从开捕年龄t_C以后，第D天1天内的回捕效益为：

从开捕年龄t_C以后，第D天以前的总回捕生物量：

从开捕年龄t_C以后，第D天以前的总回捕效益：

其中：N_tC为放流后开捕年龄t_C时的残存尾数，t_r为放流时种苗的年龄，M为自然死亡系数，F为捕捞死亡系数，W_∞为渐近体重，K为生长曲线的平均曲率，t₀为理论上体长和体重等于零时的年龄，b为幂指数系数，t_D为放流后第D天时的年龄，t_(D+1)为放流后第D+1天时的年龄，V为某一增殖放流种类单位生物量的市场平均价格。

2.9   各种类型渔船的捕捞努力量标准化

对于渔捞日志登记资料，要对各种类型渔船的捕捞努力量标准化，再根据渔捞日志登记船的增殖放流种类回捕资料，估算所有以增殖放流种类为捕捞对象的捕捞力量的总回捕量。

效能比是某任一渔船与所选定的标准船在相同渔场、相同资源密度、捕捞时间相同的条件下其CPUE(单位网次产量或单位时间产量)的比值^[2-4]。设有A、B、C、D等几种类型的渔船在同一海区作业，其各类渔船的捕捞能力各不相同，而假设每一类中的所有渔船捕捞能力相同，根据渔业统计资料选取A类船为标准船，则其余各类渔船的效能比为：

r_A=(CPUE)_A/(CPUE)_A

r_B=(CPUE)_B/(CPUE)_A

r_C=(CPUE)_C/(CPUE)_A

r_D=(CPUE)_D/(CPUE)_A

从而可以根据标准船的回捕情况估算其它类型渔船的回捕情况。

3.   讨论

目前，国内外对标志放流进行了许多研究，但国外学者们往往仅从理论上探讨，而从实际试验所得的数据和资料加以应用和分析的报道极少；国内有关根据标志放流的重捕资料进行群体生物统计量计算的至今为止仅见于费鸿年和张诗全^[2]的报道，其计算结果与生产实际也存在较大差距。对于大规模增殖放流，尚未见有群体生物统计量计算的报道。

本文提出的“放流效果统计量评估法”，是以已有或估算的放流种类的生物学等数据为输入参数来估算海区的最适增殖放流数量；利用提出的方法从渔获物中确认放流种苗的回捕数量后，使用推导公式估算标志放流回捕资料、放流前本底调查和放流后跟踪调查资料(市场调查及渔民座谈、渔业资源调查)、渔捞日志登记及渔业生产统计资料等各种调查资料对应的捕捞死亡系数，定量估算出时间序列的放流群体的残存量、回捕量、回捕率和回捕效益等，从而评估生产性或研究性放流后不同时期的增殖放流效果。