Effects of probiotics on the dynamic of bacteria in marine shrimp pond
-
摘要:
2001年5~7月，在深圳东海岸水产公司南澳半封闭式斑节对虾Penaeus monodom精养基地，进行了有益微生物的应用对虾塘总异养细菌和弧菌数量影响的调查研究。结果表明，对虾养成过程中，施用有益微生物的实验组虾塘水体中总异养细菌和弧菌数量明显较未施用的对照组虾塘低，但各虾塘水体中总异养细菌和弧菌数量的变化特征基本相同，即在养殖前期，虾塘水体中总异养细菌和弧菌数量均较稳定，并处于相对较低的水平; 而在养殖中、后期，总异养细菌和弧菌数量均急剧升高，尤以对照组虾塘的升幅最大，其弧菌数量甚至超过100×102CFU · mL-1的对虾发病之弧菌数量临界值。实验组与对照组虾塘中表层沉积物中，总异养细菌和弧菌数量差异不大，总异养细菌和弧菌数量变化特征也基本相同，但实验组与对照组之间有所差异，其中实验组表现为双峰型的变化特征，峰值出现在养殖前期或中期和养殖将结束时; 对照组则呈单峰型变化，峰值均出现在养殖中期。有益微生物的应用对虾塘水体中总异养细菌和弧菌数量的抑制效果较对沉积环境中的效果明显。
Abstract:In the Shenzhen East Sea Aquaculture Company, the effects of probiotics on the dynamic of heterotrophic bacteria and vibrio in semi-enclosed intensive marine shrimp aquaculture system were studied during the adult aquaculture period of Penaeus monodom from May to July, 2001. The results show that heterotrophic bacteria and vibrio density in water of experimental ponds, where probiot-ics were used, are obviously higher than that in control pond where no probiotics were used, and their dynamic characteristics are al-most the same. During the prophase of shrimp culture, the heterotrophic bacteria and vibrio density in water kept in a relatively low lev-el, while at the middle and final phase, both increased rapidly, especially in the control pond, the vibrio density even exceeded the threshold of 100×102 CFU · mL-1, that was considered to cause disease problem easily. There are no obvious differences for hetero-trophic bacteria and vibrio density in sediment between experimental ponds and control ponds, and their dynamic variation characteristics are almost similar, but there are some differences between experimental ponds and control pond. In experimental ponds, there were two peaks of individual density appeared in the prophase or middle phase and near the harvest time, while in control pond, the peak of individual density appeared at middle phase. It is obvious that the application of probiotics can depress the increase of heterotrophic bacteria, especially vibrio in shrimp pond water are more significant than that in sediment.
-
Keywords:
- marine shrimp pond /
- probiotics /
- dynamic of bacteria
-
实际种群分析(virtual population analysis，VPA)又称有效种群分析和现实种群分析，是渔业资源评估的经典方法之一。VPA被广泛应用来评估渔业的历史数据，求解种群资源量和捕捞死亡率[1-3]。在FRY等[2]前人研究工作的基础上，GULLAND[4]提出并应用了这个方法。
然而，鱼类年龄数据的不易获得[5]，对年龄结构VPA的应用造成了困难。由于更易获得体重数据，采用体重结构的VPA(WVPA)对于实际渔业来说更具有优越性。尤其是对于远洋渔业，像大西洋鳕(Gadus macrocephalus)，青鳕(Theragra chalcogramma)等，一般都是在渔船上按照重量进行分类，然后进行冷藏，对于上岸统计来说，只需将每箱的重量乘以箱数即可得出相应的渔获量，简单易行。伴随着捕捞努力量的增大，我国近海的渔业资源呈现越来越严重的低龄化问题[6]。由此导致年龄结构的VPA模型难以获得比较长的世代年龄序列，使得估算过程困难。而采用WVPA模型，将考察资源数量变动的时间单位由“年”改变为以鱼的个体生长间隔△W所经历的时间，将资源数量分为多个“伪世代”，使得计算步骤简单。
本文通过Monte Carlo模拟数据估算捕捞死亡系数，研究在不同的白色噪音下WVPA的可靠性及稳定性，进而回溯渔业开发史以及资源量状况[7]。并将WVPA的结果与年龄结构的VPA对比，初步探究WVPA在渔业资源评估中的可用性。
1. 方法和数据
1.1 方法-年龄结构的VPA(age-based virtual population analysis)
GULLAND[8]根据Baranov的渔获量方程提出了VPA的公式如下：
$$ N_{i, j}=\frac{F_{i, j}+M}{F_{i, j}\left(1-e^{-\left(F_{i, j}+M\right)}\right)} C_{i, j} $$ (1) $$ \frac{N_{i, j}}{N_{i+1, j+1}}=\frac{F_{i, j}\left(1-e^{-\left(F_{i, j}+M\right)}\right)}{\left(F_{i, j}+M\right) \times e^{-\left(F_{i, j}+M\right)}} $$ (2) 模型当中Ni，j是i年j龄的资源量，Ni+1，j+1是下一年该世代的资源量，Ci，j是i年该世代的渔获量，M是自然死亡系数，Fi，j是i年j龄的捕捞死亡系数，N、C的单位都是尾数。
VPA的计算流程如下：通过假设最大年龄的捕捞死亡系数F，加上与之对应的渔获量Ci，j，进而估算出Ni，j，然后将Ni，j及Ci-1，j-1带入公式(2)求算出Fi-1，j-1。以此类推，以达到重现渔业历史的目的。
1.2 模拟数据
Monte Carlo模拟仿真又称随机模拟方法(random simulation)，有时也称统计实验(statistical testing)，是人们对随机事件分析的一种数学方法，它使用不同的随机数反映随机过程的涨落现象。Monte Carlo模拟仿真常被应用于检验结果的有效性和在一定程度上分析具体的模型和渔业数据[9]。在本文的Monte Carlo模拟过程中我们采用了相同的操作模型和评估模型，其简单的执行过程如下：(1)根据原始数据求出模型的各个参数真值；(2)以此模型作为操作模型，并根据参数真值计算出模拟数据，在其上加入白色噪音[其不同水平用变异系数(coefficient of variation，CV)表示]后，即得到“伪观测数据”；(3)用模型分析伪观测数据，计算出种群的生物学参数估计值，并与参数真值比较得到相对误差。
在Monte Carlo模拟过程中我们设定了原始资源量(N1)、自然死亡系数和捕捞死亡系数。假定它们的初始值分别为：N1等于1 000，F在0.1~1.5之间；在此我们构造了2个具有不同生长特点和自然死亡率的模拟种群[10]，一个具有8龄的寿命，即：tM=8，M=0.32；另一个寿命相对短，约为3龄，tM=3，M=0.86。
在模拟过程中，资源量(N)和渔获量(C)分别用下式计算[2]：
$$ N_{t+1}=N_t e^{-(F+M)} $$ (3) $$ C_t=\frac{F}{F+M} N_t\left(1-e^{-(F+M)}\right) $$ (4) 本文设定5个不同水平的白色噪音[11]加在资源量上(CVN=1%，5%，10%，20%和30%)，而将资源量噪音水平的一半加在渔获量上(CVC=0.5%，2.5%，5%，10%和15%)。正态分布的随机变量通过Box-Mueller方法[12]产生。选择2个单一分布的随机数U1和U2并设定：
$$ Z_1=\sqrt{-2 \log \left(U_1\right)} \cos \left(2 \pi U_2\right) $$ (5a) $$ Z_2=\sqrt{-2 \log \left(U_1\right)} \sin \left(2 \pi U_2\right) $$ (5b) Z1和Z2便是平均值等于0，方差为1的正态随机变量。
本研究重复模拟1 000次，以获得稳定的解。在求算出不同年份的捕捞死亡系数F后，本文计算了参数真值和估计值的相对估计误差REE：
$$ \mathrm{REE}=\frac{\sum\limits_i^n\left|F_{\text {true }}-F_{\text {esti }}\right|}{n F_{\text {true }}} \times 100 \% $$ (6) n是模拟重复的次数，这里n=1 000。REE愈小表明估计值愈精确。
1.3 实际渔业数据
本文采用赤鲷Pagrus pagrus为实验对象。赤鲷是一种底栖性鱼类，在大西洋西部多分布在美国的纽约到阿根廷，在大西洋东部，从英国诸岛到塞内加尔一直到非洲最西端的佛得角都有它的分布[13]。赤鲷的寿命较长，PAJUELO报告最高为14龄，自然死亡系数M采用PAULY[14]的公式计算得M=0.32。
本文采用VAUGHAN等[15]在美国东南部沿海统计的1978~1997年的1~8龄渔获量数据，将最高年龄的捕捞死亡系数设为F= 1.2。通过对各龄的捕捞死亡系数的估算，进而估算资源量，达到对资源历史再现的目的，为渔业管理提供参考数据。
1.4 体重结构的VPA (weight-based virtual population analysis，WVPA)
WVPA是将第i年的鱼按体重分为j级，则Ci，j表示第i年j重量级的渔获量，它的计算公式如下，
$$ C_{i, j}=\frac{\delta_{i, j} V_i}{\gamma_j} $$ (7) δi，j：第i年的第j重量级的比例。
Vi：第i年的总产量。
γj：第j重量级的中值。
之后的过程与VPA相似，由此计算得到捕捞死亡系数和资源量。
2. 结果与讨论
2.1 应用WVPA对2个模拟种群不同白色噪音下的F的估算
表 1总结了对2个模拟种群应用WVPA计算捕捞死亡系数的结果。结果表明，即使在较高的白色噪音水平下(30%)，WVPA仍然可以较准确的估算捕捞死亡系数[相对估算误差REE=3.81%(长寿命种群)，和2.75%(短寿命种群)]。这可能指示这种方法更适用于寿命短的鱼种。
表 1 应用WVPA对2个模拟种群在不同白色噪音(CV)下捕捞死亡系数F的平均估计值和相对估计误差Table 1. Estimated average F values and relative estimated errors for two simulated populations for different white noises(CV)a F CVN=1% CVN=5% CVN=10% CVN=20% CVN=30% 0.1 0.0995 (0.49) 0.0995 (0.48) 0.0994 (0.55) 0.0995 (0.47) 0.0996 (0.40) 0.6 0.5994 (0.09) 0.5992 (0.13) 0.6028 (0.47) 0.6034 (0.53) 0.6082 (1.36) 1.0 0.1995 (0.05) 1.0100 (1.00) 1.0011 (0.11) 1.0148 (1.48) 1.0337 (3.37) 1.5 1.5038 (0.25) 1.5047 (0.31) 1.5087 (0.54) 1.5551 (3.67) 1.6518 (10.12) 平均REE/%
average0.22 0.48 0.42 1.543.81 b F CVN=1% CVN=5% CVN=10% CVN=20% CVN=30% 0.1 0.0995 (0.50) 0.0995 (0.47) 0.0997 (0.35) 0.0997 (0.35) 0.0990 (0.99) 0.6 0.5995 (0.09) 0.5996 (0.06) 0.5997 (0.04) 0.6035 (0.59) 0.6075 (1.25) 1.0 0.9993 (0.01) 0.9995 (0.05) 1.004 (0.41) 1.0115 (1.15) 1.0290 (2.91) 1.5 1.4992 (0.05) 1.4908 (0.61) 1.5205 (1.37) 1.5424 (2.68) 1.5879 (5.86) 平均REE/%
average0.16 0.30 0.54 1.19 2.75 注：a: 长寿命种群tM=8, M = 0.32; b: 短寿命种群tM=3, M = 0.86
Note: a: long-lived population tM=8, M=0.32; b: short-lived population tM=3, M =0.862.2 应用VPA与WVPA估算赤鲷的捕捞死亡系数和资源量
图 1显示的是采用体重级别对上岸鱼进行分组得到的渔获量组成。图 2是应用WVPA估算的资源量的组成状况。资源量在1980年前后出现高峰，渔获量也在这之后达到最高，但随着1982年的捕捞死亡系数的增大，使得资源量在此之后一蹶不振。在将1992~1997年的资源量由尾数转变为重量后，WVPA的计算结果与2002年6月的赤鲷资源状况报告的数据基本相符[16]。图 3为应用VPA与WVPA估算赤鲷的捕捞死亡系数，在某些年份(1985、1986、1993~1997)，F的差异较大(REE>10%)，其他年份的捕捞死亡系数差别不大。
![]()
图 1 1978~1997年赤鲷按照体重构成的渔获数量Figure 1. Catch-at-weight Ci, , j of red porgy in 1978~1997![]()
图 2 采用体重结构的VPA(WVPA)估算1978~1997年赤鲷的资源量Figure 2. Estimated population abundance of red porgy over 300 g using WVPA![]()
图 3 采用VPA与WVPA 2种方法对赤鲷1978~1997年的捕捞死亡系数的评估结果Figure 3. Estimated fishing mortality of red porgy in 1978~1997 using VPA and WVPA图 4显示的为采用年龄结构的VPA(>3龄)与体重结构的VPA(＞300 g)估算的赤鲷资源量状况的比较结果。2种方法估算的资源量均在1980年达到最大值。从1978~1997年采用VPA较采用WVPA估算的资源量偏高的原因可能是渔获量数据中大于3龄的比大于300 g的数量大的缘故。在1984~1997年估测的资源量中采用WVPA较VPA获得的结果变化趋势较缓和，这或许是由于一个世代中包含几个重量级，或者一个重量级包含几个年龄世代的缘故。这是因为在实际的渔业中，鱼类的补充并不是由年龄，而是由体长或体重决定的[17]。因此，采用体重结构的VPA，即WVPA，在渔业中更具有实际的意义。
![]()
图 4 采用年龄结构的VPA(>3龄)与体重结构的VPA(＞300 g)估算的赤鲷资源量Figure 4. Comparison of red porgy abundance estimated using WVPA (over 300 g) and VPA (over age 3)此外，还有一些问题，例如WVPA与捕获比率及鱼类生长曲线的相关，对上岸量进行重量级别分组需要采用多大的间隔等等，还有待于日后的进一步研究。
2.3 误差来源分析
(1) 本研究假设模型没有过程误差，使用该方法时渔获量数据含有误差是结果出现误差的主要原因。渔业中一些渔获物没有被统计到的情况时常发生，例如私人消费渔获物以及海上抛弃[18-19]。仅凭渔民捕捞日志的资料有时不能得到可靠的渔获数据，这可能低估了产量，这样捕捞死亡系数的估计就会偏高。
(2) 在评估过程中，将各龄的自然死亡率确定为同一值的做法，也会对结果造成一定的影响。
(3) 一般来说，在采用VPA来估算资源量的过程中，我们获得的捕捞数据是根据年龄测量得到的数据，样本的数量受到很大的限制。而采用WVPA样本的数量几乎可以等同于上岸鱼的数量。因此在减小误差方面，WVPA会更有优势。
3. 总结
本文演示了由体重结构的VPA估算鱼类捕捞死亡率和资源量的方法。WVPA的计算精度近似于年龄结构的VPA，但它对数据的要求要低一些。所以当实际渔业中的年龄数据不存在或不够精确时，可以采用本文的方法。
-
![]()
图 2 虾塘水体中总异养菌数量变化
Figure 2. Density variation of heterotrophic bacteria in the water column
![]()
图 4 虾塘表层沉积物中总异养菌数量变化
Figure 4. Density variation of heterotrophic bacteria in the sediment
表 1 对虾养殖各阶段虾塘水体总异养细菌平均数量
Table 1 Average density of heterotrophic bacteria in the water column at different culture stages
×104CFU · mL-1 养殖阶段
aquaculture phase虾塘号 pond no. 3 4 5 前期 prophase 7.2 22.2 10.2 中期 middle phase 10.0 11.7 31.7 后期 final phase 46.3 64.7 71.0 
下载: 导出CSV
表 2 对虾养殖各阶段虾塘水体总弧菌平均数量
Table 2 Average numbers of vibrio in the water column at different culture stages
×102CFU · mL-1 养殖阶段
aquaculture phase虾塘号 pond no. 3 4 5 前期 prophase 1.8 2.3 5.0 中期 middle phase 2.9 5.8 8.8 后期 final phase 39.0 24.0 97.7
下载: 导出CSV
表 3 对虾养殖各阶段虾塘表层沉积物中总异养细菌平均数量
Table 3 Average numbers of heterotrophic bacteria in the sediment at different culture stages
×106CFU · g-1 养殖阶段
aquaculture phase虾塘号 pond no. 3 4 5 前期 prophase 6.4 7.3 4.6 中期 middle phase 4.9 16.0 10.0 后期 final phase 7.0 9.7 2.3
下载: 导出CSV
表 4 对虾养殖各阶段虾塘表层沉积物中弧菌平均数量
Table 4 Average numbers of vibrio in the sediment at different culture stages
×103CFU · g-1 养殖阶段
aquaculture phase虾塘号 pond no. 3 4 5 前期 prophase 48.6 16.0 20.0 中期 middle phase 4.5 36.0 36.0 后期 final phase 38.7 90.0 8.3
下载: 导出CSV
-
[1] 江育林. 对虾养殖、环境、病害[J]. 中国水产, 1994(5): 35-36. http://ir.ihb.ac.cn/handle/152342/5044 [2] PAEZ-OSUNA F, GUERRERO-GALVAN S R, RUIZ-FERNAN-DEZ A C. The environmental impact of shrimp aquaculture and the coastal pollution in Mexico[J]. Mar Poll Bull, 1998, 36(1): 35-36. doi: 10.1016/S0025-326X(98)90035-2
[3] PAEZ-OSUNA F, GUERRERO-GALVAN S R, RUIZ-FERNAN-DEZ A C. Discharge of nutrients from shrimp farming to coastal wa-ters of the gulf of California[J]. Mar Poll Bull, 1999, 38(7): 585-592. doi: 10.1016/S0025-326X(98)00116-7
[4] 董双林, 潘克厚. 海水养殖对沿岸生态环境影响研究进展[J]. 青岛海洋大学学报, 2000, 30(4): 575-582. doi: 10.3969/j.issn.1672-5174.2000.04.030 [5] 吴伟. 应用复合微生物制剂控制养殖水质因子初探[J]. 湛江海洋大学学报, 1997, 17(1): 16-20. https://kns.cnki.net/kcms2/article/abstract?v=NK8hpUzgeRWFVBtcOFEXhEOpl3NDN5n2c48c7mvGsb-K3d2FYgPuXkgXBJvtBEZyGlhS9zq8ZO9fulUBLNZhay4Ty8iFSllisXU_mTk506k-xprXYQksh4syWQk2vyVu52CHdRiYXf8zfHGHvYadTB0yJHtJVGtxOmdFj0XFCBCfmxhC_Xes_aH5JCUwuvLK&uniplatform=NZKPT&language=CHS [6] 李卓佳, 张庆, 陈康德. 池塘养殖中应用微生物防治斑节对虾疾病的初步研究[M]//何其锐. 南海资源开发研究. 广州: 广东经济出版社, 1998: 703-709. [7] 李卓佳, 张庆, 陈康德, 等. 应用微生物健康养殖斑节对虾的研究[J]. 中山大学学报: 自然科学版, 2000, 39(增刊): 229-232. doi: 10.3321/j.issn:0529-6579.2000.z1.049 [8] 刘卫东, 苏浩, 邓立康. 微生物在水产养殖中的应用[J]. 水产科学, 2001, 20(2): 28-31. doi: 10.3969/j.issn.1003-1111.2001.02.014 [9] 吴琴瑟. 微生态制剂对虾类健康养殖和环境保护的作用[M]//张本. 虾类养殖研究. 北京: 海洋出版社, 2002: 42-45. [10] 叶乐, 杨莺莺, 吴开畅, 等. 益生菌在凡纳滨对虾育苗中的应用[J]. 南方水产, 2006, 2(4): 13-18. doi: 10.3969/j.issn.2095-0780.2006.04.003 [11] 张汉华, 李卓佳, 郭志勋, 等. 有益微生物对海水养虾池浮游物生态特征的影响[J]. 南方水产, 2005, 1(2): . 7-14. doi: 10.3969/j.issn.2095-0780.2005.02.002 [12] 李卓佳, 林亮, 杨莺莺, 等. 芽孢杆菌制剂对凡纳滨对虾Litopenaeus nannamei肠道微生物群落的影响[J]. 南方水产, 2005, 1(3): 54-59. doi: 10.3969/j.issn.2095-0780.2005.03.009 [13] VISSCHER P T, DUERR E O. Water quality and microbial dy-namics in shrimp ponds receiving bagasse-based feed[J]. J World Aquac Soc, 1991, 22(1): 56-76. doi: 10.1111/j.1749-7345.1991.tb00718.x
[14] 高尚德, 陈旭仁, 吴以平. 中国对虾养成期间虾塘水体和底质中细菌含量的变化[J]. 水产学报, 1994, 18(2): 138 -142. doi: 10.1007/BF02007173 [15] 郭平, 许美美. 对虾养殖池水域环境细菌的动态变化[J]. 海洋与湖沼, 1994, 25(6): 625-630. https://kns.cnki.net/kcms2/article/abstract?v=NK8hpUzgeRWJ5yZg-01fr0BCtXSHyJbeG2oWfIzQGFb2_Jegrr5UT4WAT08N6AOZgSSdI2TcqSKpxdVqIMGzQGGO84F340bLT-6Jqss3pnxZ2zcZ6vx8tKKCiX96TXiVgUnDEAKa1dxbVhfw-DwgQbVSL9BFx-Uaj9YDpHqIKwY-PT4Y8VDr6MU9mZEKqEJ4&uniplatform=NZKPT&language=CHS [16] 刘国才, 李德尚, 徐怀恕, 等. 海水养虾塘细菌数量动态及细菌生产力的研究[J]. 应用与环境微生物学报, 1997, 3(4): 340-344. https://kns.cnki.net/kcms2/article/abstract?v=NK8hpUzgeRWzhzCnzWQEsSqcDXOC98CrugbzF1wwBa-J_-qHjV13aiLrnPSSv4osnYlL-u9o1WrExI5uOj9XVibqh_pUsqB4immzoLsOTeH_VRl2mQOMQ0OzIKIt3Bc0dgeERAQGkaFkOABZ7eyYSxUcvWtxuFoXEiDr-rNVqnPhnppcC2Pv3v7NDzqPPXfW&uniplatform=NZKPT&language=CHS [17] 李秋芬, 曲克明, 陈碧鹃, 等. 老化虾塘生态系中几类主要细菌的季节变化特征[J]. 海洋水产研究, 2002, 23(2): 12-18. https://kns.cnki.net/kcms2/article/abstract?v=NK8hpUzgeRXjR1ib4osPezaBmKOmuMyPEVJQjGvGZ7HIAwdTOAlk0-VcG6tbh5Y3RQgEkHalDZYcEEicbHcuuuoDH4XGh5Xurw_oq3loh7yPBPSNjbylO6YZVQQZ8V6akAoAtdSvMR7kCiezapQ0dH95N04fTwyC_mup3rz6h7ZAF6FXhLucTvLmB036HUSM&uniplatform=NZKPT&language=CHS [18] 李卓佳, 郭志勋, 张汉华, 等. 斑节对虾养殖池塘藻-菌关系初探[J]. 中国水产科学, 2003, 10(3): 262-264. doi: 10.3321/j.issn:1005-8737.2003.03.016 [19] 黄洪辉, 林钦, 李卓佳, 等. 对虾养成过程中水质与异养细菌变化特征及其关系研究[M]//张本. 虾类养殖研究. 北京: 海洋出版社, 2002: 339-343. [20] 海洋监测规范编辑委员会. 海洋监测规范[S]. 北京: 海洋出版社, 1991. [21] 于占国, 林凤翱, 贺杰. 异养细菌与虾病关系的研究[J]. 海洋学报, 1995, 17(3): 85-91. https://d.wanfangdata.com.cn/periodical/Ch9QZXJpb2RpY2FsQ0hJTmV3UzIwMjQxMTA1MTcxMzA0Eg5RSzE5OTUwMDQzMTAzMhoIZ2U5NTF6YnA%3D [22] 田黎, 李光友. 海洋生境芽胞杆菌(Bacillus spp. )的培养条件及产生的胞外抗菌蛋白[J]. 海洋学报, 2001, 23(4): 87-92. doi: 10.3321/j.issn:0253-4193.2001.04.011 [23] 俞勇, 李会荣, 李筠, 等. 益生菌制剂在水产养殖中的应用[J]. 中国水产科学, 2001, 8(2): 92-96. doi: 10.3321/j.issn:1005-8737.2001.02.024
计量
- 文章访问数: 5117
- HTML全文浏览量: 132
- PDF下载量: 3782
粤公网安备 44010502001741号
