南美白对虾品系、密度、放苗时间与养殖效益的统计学分析

李举鹏，史健，胡小弘，周影，王彩蕴，单红云，孙丽敏

(1. 营口市水产科学研究所,辽宁 营口 115003；2. 老边区水产技术推广站,辽宁 营口 115005；3. 营口市水产技术推广站,辽宁 营口 115003)

南美白对虾是世界养殖产量最高的三大优良虾种之一[1]，自1988年引入我国[2]。目前，我国南美白对虾的养殖模式多种多样，不同的养殖模式对虾苗的品系和要求也在不断变化，但随着养殖时间的推移，病害加剧，养殖难度愈来愈大，养殖成功率很低，国内很多学者在养殖池塘水质的理化因子[3]、浮游植物群落结构[4]， 水体微生物群落结构[5]以及对虾饵料[6]等方面，进行了大量的研究，这些研究往往是从试验对比的角度探求对虾发病的原因，该文拟从养殖生产的实际效果调查入手，通过统计分析的方法，以寻求对虾发病的原因， 为此，作者于2017年在营口市高坎镇太平村育苗场进行了不同品系虾苗的标粗淡化，然后将这些淡化虾苗投放到不同的养殖区域进行淡水养殖，并于2018年在浙江宁海华大海昌公司进行了同一品系(SIS)虾苗不同密度及不同季节放苗的养殖生产，通过跟踪调查，回访等方式，将收集到的数据整理，并进行统计学分析，以探求南美白对虾养殖失败的原因.

1 方法

1.1 不同品系虾苗的养殖

2017年4月25日至5月21日，先后从南方育苗场购进科纳湾、SIS和正大三个品系，共11个批次的P4-P5一代对虾苗8 000万尾，经过12～19 d的培育，共培育出体长1.0～1.5 cm的仔虾苗6 250万尾，其中科纳湾虾苗3 830万尾、SIS虾苗1 640万尾、正大虾苗780万尾，然后分别投放到营口市的水源镇、老边区、沟沿和高坎镇的土池进行淡水养殖，每个池塘的面积为0.067～0.670 hm2，放苗密度为4万尾/667 m2左右，养殖总面积约100 hm2，养殖池塘共327个。尽管各批次虾苗刚入池淡化时规格及活力有所不同，但经过十多天的淡化培育，放苗时虾苗的活力以及摄食情况等虾苗质量方面基本相差无几。

1.2 单一品系虾苗的养殖

2018年4月20日开始至2018年8月9日先后24次从海南育苗场购进南美白对虾(SIS)无节幼体5.2亿尾，共培育出体长0.8～1.2 cm的仔虾1.8亿尾，分别投放到浙江、江苏、河南、山东等地的池塘进行养殖。

1.3 数据收集

数据收集主要采取跟踪调查和走访的方式进行，每个养殖池塘作为一个观察单位，观察时间为90 d，分别对三个品系虾苗发病时的养殖天数和养殖效益进行记录，并将收集来的数据进行汇总整理。

1.4 统计分析

统计分析采用SPSS.17.0软件进行，分别对对虾养殖期间发病时的存活天数进行秩和检验、对养殖效益进行χ2检验，不同养殖密度和不同放苗季节对养殖效益的影响采取图表分析，统计学差异水平规定为P<0.05。

2 结果

2.1 秩和检验

Kruskal-Wallis检验是1952年由Kruskal和Wallis二人共同提出的，是二个独立样本检验的一种推广[7]，Kruskal-Wallis检验用于推断计量资料或等级资料的多个独立样本所来自的多个总体分布是否有差别[8]，其基本思想是：将原始资料在不分组的情况下从小到大排顺序，即编秩次，然后分组将秩次相加得到各组的秩和，若组间秩和相差不大，则认为各组间无差异。

其计算公式为：

(1)

式中：H——统计量；

Ri——各组的秩和；

ni——各组对应的例数；

N=∑ni。

2.2 同一品系虾苗在不同区域养殖 虾苗的发病情况

高坎镇和沟沿镇由于养殖面积相对较少，所以合并为一个区域，每个池塘虾苗的存活时间按照对虾发病时的日期与放苗日期之间的天数之差计算。各养殖区以养殖科纳湾虾苗的发病池塘为选取条件；老边区发病池塘较多，从中随机选取30个发病池塘做为样本，各发病池塘虾苗的存活时间如表1。

经SPSS.17.0计算得H=6.250，得P<0.05， 差别有统计学意义，其中，均值为52.95 d，标准差为13.215，即：虾池平均发病时间是53 d，三个区域中虾池最早发病时间是25 d，最晚发病时间是80 d，说明水源镇、老边区、高坎和沟沿镇虾苗的发病时间是有差异的。

2.3 不同品系虾苗在同一地区养殖 虾苗的发病情况

老边区养殖面积相对较大，发病比例也最高，在此区域内，每个品系的虾苗随机抽取15个发病池塘，每个池塘虾苗的存活时间也是按照对虾发病时的日期与放苗日期之间的天数之差计算，各品系虾苗存活时间如表2。

表2 同一区域不同品系虾苗养殖存活天数

经SPSS.17.0计算得H=1.527，P>0.05，差别没有统计学意义，其中，均值为59.69 d，标准差为12.989，说明三个品系的虾苗在老边区养殖平均的发病时间是57 d，最早发病时间是29 d，最晚是80 d，但科纳湾、SIS和正大3个品系的虾苗发病时间总体分布是没有差异的。

2.4 χ2检验

χ2检验是英国统计学家Pearson于1900年提出的一种用途很广的计数资料的假设检验方法[9]。属于非参数检验，主要是比较两个及两个以上样本率(构成比)以及两个分类变量的关联性分析。根本思想在于比较理论频数和实际频数的吻合程度或者拟合优度，其计算公式为：

(2)

式中：A——实际频数；

T——理论频数。

2.5 不同区域、不同品系对虾养殖效益情况

对虾养殖的最终目的是为了获得经济效益，将不同区域，不同品系虾苗的养殖效益进行汇总分析，更有利于揭示它们之间的关系，该次共调查了327个养殖池塘，其中，老边区发病严重，调查了144个养虾池塘。科纳湾虾苗在各区域养殖的效益情况见表3，不同品系虾苗在同一区域养殖的效益情况见表4.

表3 科纳湾虾苗在不同区域的养殖效益情况(以池塘个数计)

经SPSS.17.0计算得χ2=31.792，P<0.05， 差别有统计学意义，从统计结果来看，高坎和沟沿镇养殖效益最好，老边区养殖效果最差，水源镇介于二者之间，3个区域养殖效益差异特别明显。

表4 不同品系虾苗在老边区的养殖效益情况(以池塘个数计)

经SPSS.17.0计算得χ2=1.892,P>0.05， 差别没有统计学意义，统计结果说明3个品系的虾苗在同一区域养殖，就总体而言，其养殖效益差异不明显。

2.6 不同养殖密度下池塘的发病情况

2018年，由于放苗面积较广，路途较远，有些调查内容难以获得，因此只对各批次的虾苗是否发病及养殖密度进行了简单的调查，养殖密度分为4个梯度，即：1～2万尾/667 m2、2～3万尾/667 m2、3～4万尾/667 m2、4～6万尾/667 m2。池塘虾苗发病比率分别是：38%，41%，47%，68%，如图1所示。

图1 不同养殖密度下发病池塘所占比率

结果显示，随着养殖密度的增大，虾池的发病比例也在随之增加，尤其高温时期放苗，当养殖密度达到4万尾以上/667 m2时，土池养殖的发病比率非常高。

2.7 不同放苗时间，池塘的发病情况

今年从4月20日开始到8月9日结束，先后24批次，共培育出0.8～1.2 cm的仔虾1.8亿尾，分别投放到浙江、江苏、、河南、山东等地的池塘进行养殖，放苗结束后，随时收集记录各批次虾苗的发病情况，根据各批次虾苗的发病情况，计算其发病比率，绘制成图，各批次虾苗发病比率如图2。

图2 不同批次虾苗发病情况

不同季节放苗，各批次虾苗的发病比率是不一样的，最高为93%，最低为18%，4月的幼体苗，发病池塘比例在78%～83%，5月和6月的幼体苗，发病池塘比例较低，发病池塘比例维持在18%～3%，7月以后的幼体苗，发病虾池比例为64%～93%，而且，高温期放苗往往伴随着虾苗的生长缓慢和参差不齐。

3 分析与讨论

3.1 同一品系虾苗在不同区域养殖虾苗的发病情况

同一品系的虾苗投放到不同的养殖区域进行养殖，发病时间是存在差异的，从表1的统计结果可以看出，3个品系的虾苗在不同养殖区域的发病时间具有明显的差异，这与养殖的实际情况是吻合的。秦志浩等[10]在对不同品系南美白对虾室内工厂化养殖的研究中认为：我国南美白对虾的养殖模式多种多样，不同养殖模式对虾苗的品系和要求也在不断变化，在工厂化养殖模式下，正大品系虾苗最适合养殖，产量较高，个体较大；但杨文慧等[11]曾对天津地区南美白对虾产业的调查分析指出：南美白对虾养殖过程中，确实存在地区化差异和年份差异，使得南美白对虾在不同区域养殖效果表现不一。该次调查的结果与杨文慧等人的观点是一致的，从营口地区的养虾实际情况来看，对虾养殖的发病时间确实具有区域性的特点，即：在某一时间段同一区域的绝大部分养殖池塘同时或相继发病，这很可能与养殖区域的水质以及气候状况有关。

3.2 不同品系虾苗在同一地区养殖虾苗的发病情况

很多养殖户认为，对虾的发病与虾苗的品系有关，对虾养殖期间的发病往往都是由于虾苗品系不好造成的，但从表2的统计结果来看，3个品牌虾苗在同一区域养殖，其发病时间是没有统计学差异的。张嘉晨等[12]在进行高密度养殖凡纳滨对虾(Litopenaeusvannamei)新品系选育时的研究认为：以不同国家和地区的8个凡纳滨对虾养殖群体作为亲本，构建8个子代群体，美国凡纳滨对虾改良公司(SIS)♀×SIS♂生长优势明显，生长比速较群体均值提高10.4%～17.9%;海大选育(HD)♀×粤海选育(YH)♂成活优势明显，成活率较群体均值提高9.3%～13.7%;YH♀×中兴1号(ZX)♂和YH♀×科纳湾种虾(K)♂生长比速、成活率均高于群体均值，表现出生长快成活率高双重优良性状，该次调查的结论与上述观点不一致，是否与各地的养殖环境不同有关，还有待于进一步的研究。

3.3 不同区域、不同品系虾苗养殖效益分析

从表3、表4的分析结果看出，同一品系虾苗在不同区域养殖，虾苗的发病时间存在差异，导致养殖效益也存在差异；不同品系虾苗在同一区域养殖，养殖效益没有显著性差异。陈晓明等[13]在进行微生态制剂防控南美白对虾病害试验发现：经济效益的高低与放苗的品种有直接关系，即：一代苗与土苗对比，养殖效益差异特别明显；但不同品系的一代苗之间(SIS和桂海)，养殖效益差异不大，这一结论与该次调查结果基本是吻合的，基于这一结论，建议在对虾养殖的实际生产中着重注意水环境的改善以及养殖技术的提高。

3.4 不同养殖密度下池塘的发病情况分析

由图1可以看出，随着养殖密度的增加，发病池塘所占的比例越来越大，在养殖过程中，特别是高温季节，随着个体的增大，尤其对虾体长达到6～7 cm以后，养殖水体中的残饵、粪便大量增加，水环境严重恶化，加速了对虾发病的几率。李玉全等[16]认为：凡纳滨对虾(体长<7.6 cm或体质量<6.1 g)养殖密度为150～900尾/m3时，造成对虾生长和非特异性免疫因子差异的主要原因是密度胁迫，水质因子的作用是次要的。张国新[15]曾就不同养殖密度对南美白对虾生长的影响做过试验，试验设置养殖密度分别为75尾/m2、150尾/m2、300尾/m2、450尾/m2、900尾/m25种梯度，结果表明：不同密度对南美白对虾的成活率有显著影响，各不同试验组的成活率分别为93.7%、94.7%、84.7%、77%、69%，随密度增加，成活率呈下降趋势，各组间成活率最大相差25.7%[14]。张吕平等报道了在我国热带地区集约化养殖凡纳滨对虾(L.vannamei)4种密度下(70、100、120、150尾/m)养殖成活率分别为80.44%、75.16%、51.74%和44.43%。上述观点在实际养殖生产中也得到了印证，大面积低密度的大围池塘养殖成功率明显高于小面积的精养池塘。

3.5 不同放苗时间，池塘的发病情况分析

由图2可以看出，同一品系的虾苗，5、6月放养的各批苗，发病比率比较低，4月和7月以后放养的虾苗，发病比例都比较高，主要是由于前期气温较低，温差变化较大，7月进入高温期以后，由于受气候和天气状况的影响，水中的理化因子及藻类密度、结构变化剧烈，导致水质不稳定，从而造成虾病的发生。彭荣南等[17]在台风引起养殖水质突变的研究中发现：风切变引起水环流扰乱，有毒物质扩散，大风引起对虾过度应激反应;同时，伴随阴雨寡照使得水体光合作用减少(DO减少)，DO达到最低，对虾病害明显；沈明明[18]在对虾工厂化养殖水体微藻群落演替的主要因子及其与疾病发生的关系中指出：微藻的密度和生物量显著影响对虾发病程度，过高的微藻密度尤其是养殖后期蓝藻优势种的密度，与养殖后期对虾持续发病关系密切；过低的微藻密度和生物量则会造成微藻对虾池生态系统调控能力低下，导致对虾生长缓慢，个头参差不齐。

对虾病害的发生及养殖效益，受多种因素的影响，其中，养殖环境的污染导致病害的发生，国内外很多学者、研究人员进行了大量的研究[19-21]，从该次调查的结果来看，虾苗品系不是影响养殖效益的根本因素，养殖密度、放苗时间、养殖水环境以及天气等因素都会导致对虾病害的发生，从而影响养殖效益。