对虾是一种具有较好经济效益的水产养殖种类[1]。南美白对虾(Litopenaeus vannamei)因具有抗病性强、生长快和肉质鲜美等优点,已成为对虾主要养殖品种之一[2]。提升养殖产量、降低肉饲比是对虾养殖的重要研究内容。目前,在对虾养殖中,对虾主要依赖饲料托盘观察摄食[3],由于虾塘水体环境浑浊复杂,通过人工很难准确掌握对虾摄食情况,需要使用其他技术手段加以观测。
近年来,被动声学技术逐渐应用到对虾养殖的饲料投喂与生长管理中。Silva等[4]研究表明,南美白对虾摄食时大颚闭合碰撞摩擦会产生声音,摄食前后声音信号特征不同[5-6],因此,可通过对虾摄食声音信号特征判断其摄食及行为状态。被动声学技术作为一种对生物无干扰监测方法,受虾塘水体环境影响较少,可应用于包括对虾在内的水中生物的发声研究[7-9]。目前,基于被动声学技术的自动声学投饵器已应用于实际生产,有效改善了饲料管理能力,提高了对虾的养殖效益[10-12]。研究对虾摄食发声信号特征,判断对虾摄食情况与行为状态[9,13],了解对虾的摄食发声行为规律,有助于构建科学智能的饲养管理体系,提升对虾养殖产量[14-16]。
已有研究表明,对虾的数量密度与体质量[17]及饲料的质量与规格[18-20]等均会影响对虾的摄食发声信号。此外,不同体长的对虾因大颚咬合力不同也会影响其摄食发声信号[21]。Soares等[22]通过对2种体长的南美白对虾摄食发声信号研究发现,大体长对虾摄食产生较高的能量;Silva等[4]进行了3种体长的南美白对虾摄食发声试验,结果发现,大体长对虾摄食速度更快,具有更稳定的摄食次数。上述研究结果均表明,对虾体长会影响对虾摄食发声信号。目前,已有的不同体长对虾摄食发声研究仅是单一的大、小体长对虾摄食发声信号对比,而不同体长对虾摄食发声信号是否存在规律,其生长信息与摄食发声信号是否存在关联尚不明确。在养殖实践中,对虾通过蜕壳实现生长[23-24],且存在对虾吃壳这一常见现象,但目前尚无对虾摄食虾壳的相关研究。若能够区分对虾摄食颗粒饲料与虾壳发声信号的不同,就可以捕捉对虾处于蜕壳后的生长时期信息,并有利于水下实时监测对虾的摄食生长状况、饲料投喂情况及蜕壳后的发育状态。
本研究中,以5种不同体长的南美白对虾为研究对象,采集南美白对虾摄食颗粒饲料与虾壳的声音信息,通过声学信号及单因素方差分析,确定5种体长的南美白对虾摄食颗粒饲料与虾壳时发声信号的变化规律及二者差异特征,以期为精细化对虾养殖及智能投饵机等的应用开发提供数据基础。
试验用虾选取养殖大池(直径3.4 m、高1.2 m)中培育的从幼虾至成虾5个生长阶段(体长2~3、4~5、6~7、8~9、9~10 cm)的南美白对虾,养殖水体温度、盐度、溶解氧和饲料无机盐均维持在适宜范围内,即24 h持续通氧,温度为28 ℃。试验前,将对虾放置在与养殖大池水体环境一致的水箱(44 cm×28 cm×30 cm)中暂养,禁食20 h,确保试验时对虾有足够的摄食欲望。试验时,随机挑选一组对虾(3尾)放入声学试验玻璃水箱(44 cm×28 cm×30 cm),保持与之前一样的水环境,分别投喂粒径为1.2 mm的颗粒饲料及虾壳,一组体长对虾试验结束后,更换另一组体长对虾进行试验。试验玻璃水箱外表面覆盖黑色薄膜,同时在地面垫一层泡沫以减少环境噪声干扰,每个水箱中均配有一个加热棒和增氧管(图1)。
图1 试验装置
Fig.1 Equipment of experiment
试验所用装置为大疆水下摄像机和声音采集装置。声音采集装置由带有储存卡的SM4型声音记录仪和1个HTI-96-Min水听器组成(图2)。水听器可接收带宽(2 Hz~30 kHz)声音,可以捕捉到对虾摄食时的声音,SM4型声音记录仪用于储存对虾摄食声音。南美白对虾从幼虾至成虾的最高摄食发声频率为40~48 kHz,因此,本研究中采样频率设置为96 kHz,储存格式为.wav。
图2 声音采集装置
Fig.2 Sound acquisition device
1.2.1 数据采集 采集数据时关闭增氧机,开启与水听器连接的SM4型声音记录仪,将水下摄像机置于相应位置并开启录制。等待约5 min对虾适应环境后,取15颗粒径为1.2 mm的颗粒饲料(或虾壳)润湿后投入水听器下方,记录对虾摄食时间,摄食后8~10 min时结束信号采集,之后进入下一组体长对虾摄食试验。通过对虾摄食饲料与虾壳的音频与视频信息对比,确定对虾摄食信号与行为活动的关系。
1.2.2 发声信号特征分析 使用 Matlab R2019b软件分析对虾摄食颗粒饲料和虾壳的音频信息,挑选有效音频数据对其时间序列信号采用傅里叶变换,分析信号中不同频率成分及信号特征[5]。离散时间傅里叶变换计算公式[25]为
(1)
式中:n为变换长度;ω为角频率;X为时域信号x的傅里叶变换(即频谱)。本文中,信号处理的采样频率为96 kHz,增益为20 dB,进行时频分析时按信号波形总时间的1/50分时窗,并利用Matlab中的Spectrogram 函数进行时频图的绘制。最后确定代表信号特征的各类声学参数。
使用Matlab进行函数定义,将对应对虾摄食声学参数嵌套入单因素方差分析数学模型Analysis of Variance中[26],显著性水平设为0.05,即
(2)
式中:X为测量值;μ为真值;ε为误差。
由于使用单因素方差分析时,数据需符合正态分布条件,因此,首先对所得到的5种体长对虾摄食颗粒饲料与虾壳发声的最低频率、最高频率、峰值频率及单个脉冲信号时长等4项声学参数值的统计数据进行Lilliefors[27]检验,即
(3)
式中:D*为检验统计量;F^(x)为平均值为0、标准差为1的正态分布累积分布函数;G(x)为经验分布函数值。
将不同体长对虾摄食的音频数据经过傅里叶变换后得到其摄食的时频、波形和频谱图。对各组9~10 cm的对虾观察采集到的信号发现,投饵前发声信号非常少,投饵后发声信号数量显著增加(图3(a)),可以确定为对虾摄食发声信号。选取其中一个波形信号放大观察(图3(b)),对虾单个摄食发声信号时长约22 ms,结合时频图和频谱图(图3(c)、(d)),确定对虾摄食发声信号的主要频率为5~13 kHz,峰值频率约7 kHz。同时还观察到,对虾摄食期间会产生较高的能量(图3(c)),波形和频率均有一定起伏变化。采用同样的方式,对其余4种体长对虾摄食声学数据进行分析,得到所有体长对虾的摄食声学参数。
图3 体长9~10 cm对虾的声学信号特征
Fig.3 Acoustic signals characteristics of the shrimp with body length of 9-10 cm
根据Lilliefors检验结果(表1),获取的所有体长对虾的4项声学参数均满足正态分布,因此,可以对试验数据进行单因素方差分析。
表1 Lilliefors检验
Tab.1 Lilliefors test
摄食类型feeding type对虾体长/cmbodylength最低频率minimum frequency最高频率maximum frequency峰值频率peak frequency单个脉冲信号时长duration of a single pulse signalh 测试结果hypothes test resultP值P valueh 测试结果hypothes test resultP值P valueh 测试结果hypothes test resultP值P valueh 测试结果hypothes test resultP值P value颗粒饲料pellet2~300.440 500.258 400.079 700.500 04~500.500 000.321 400.500 000.500 06~700.500 000.440 500.500 000.409 38~900.500 000.409 300.500 000.210 89~1000.500 000.441 700.500 000.500 0虾壳shrimp shell2~300.500 000.500 000.159 000.469 44~500.500 000.172 200.274 600.072 66~700.274 600.458 400.469 400.500 08~900.500 000.404 000.310 200.500 09~1000.500 000.500 000.500 000.159 0
注:h=0意味着数据可能符合正态分布(0.05显著性水平);如果不用MCTol,P值范围为0.001至0.50。
Note:h=0 means we don’t reject the null hypothesis at 0.05 level;Without MCTol,P-value ranges from 0.001 to 0.50.
对5种体长(2~3、4~5、6~7、8~9、9~10 cm)对虾摄食同一粒径颗粒饲料与虾壳的4项声学参数数据进行方差分析显示,对虾体长均极显著影响对虾摄食颗粒饲料与虾壳的4项声学参数(P<0.01)(表2和表3)。
表2 5种体长对虾摄食颗粒饲料声学参数的单因素方差分析
Tab.2 Single factor variance analysis of acoustic parameters in the shrimp with 5 body length fed the pellets
声学参数类型acoustic parameter type方差来源source of variance偏差平方和sum of square自由度free degree方差varianceF值F value显著性significance因素A357.36489.3435.4520.000最低频率minimum frequency误差e50.40202.52总和407.7624因素A2 906.204726.5641.7560.000最高频率maximum frequency误差e348.002017.40总和3 254.2024因素A1 230.164307.5461.0200.000峰值频率peak frequency误差e100.80 205.04总和1 330.96 24因素A353.84488.4644.2300.000单个脉冲信号时长duration of a single pulse signal误差e40.00202.00总和393.8424
表3 5种体长对虾摄食虾壳声学参数的单因素方差分析
Tab.3 Single factor variance analysis of acoustic parameters in the shrimp with 5 body length fed the shrimp shell
声学参数类型acoustic parameter type方差来源source of variance偏差平方和sum of square自由度free degree方差varianceF值F value显著性significance因素A105.04426.2638.6180.000最低频率minimum frequency误差e13.60200.68总和118.6424因素A391.84497.968.1630.000最高频率maximum frequency误差e240.002012.00总和631.8424因素A216.16454.0425.4910.000峰值频率peak frequency误差e42.40202.12总和258.5624因素A182.56445.6446.5710.000单个脉冲信号时长duration of a single pulse signal误差e19.60200.98总和202.1624
2.3.1 峰值频率变化趋势 从图4可见:对虾摄食颗粒饲料的峰值频率随对虾体长的增加而降低(拟合方程为1/y=0.181 1-0.140 65/x,R=0.999 6);对虾体长为2~3 cm时,发声峰值频率最高(约25 kHz),明显高于其他体长;对虾体长为9~10 cm时,发声峰值频率最低(约7 kHz);对虾体长从2~3 cm变化至4~5 cm时,发声峰值频率降低较显著,对虾体长从4~5 cm变化至9~10 cm时,发声峰值频率变化较小,总体呈缓慢降低趋势。
图4 峰值频率与对虾体长变化的关系
Fig.4 Relationship between peak frequency and body length of the shrimp
对虾摄食虾壳的发声峰值频率随体长增加无明显变化规律(拟合方程为y=8.856 4 e0.077 682x,R=0.432 2);对虾体长为2~3、8~9、9~10 cm时,发声峰值频率相对较高(分别约13.4、14.8、14.4 kHz);对虾体长为4~5、6~7 cm时,发声峰值频率相对较低(约9 kHz)(图4)。
2.3.2 发声频率范围变化趋势 从图5可见:对虾摄食颗粒饲料的频率范围随对虾体长的增加而缩小;发声最大频率随体长的增加显著降低(拟合方程为y=11.539 2 e1.361 8/x,R=0.99);发声最小频率也随体长的增加而降低(拟合方程为y=3.101 6 e1.412 6/x,R=0.963 7),体长从4~5 cm变化至9~10 cm时,发声最小频率无明显变化;对虾体长为2~3 cm时,发声频率范围最大(14~44 kHz),对虾体长为9~10 cm时,发声频率范围最小(5~15 kHz)。
图5 对虾摄食颗粒饲料频率范围与对虾体长的关系
Fig.5 Relationship between frequency range of feeding pellets and body length of the shrimp
从图6可见:对虾摄食虾壳的频率范围随对虾体长的增加也无明显变化规律;发声最大频率在对虾体长为6~7 cm时最低(约15 kHz)(拟合方程为y=17.671 8 e0.051 959x,R=0.413 8);发声最小频率随对虾体长的增加而增加(拟合方程为y=4.189 2 e0.150 73x,R=0.800 5);对虾体长为6~7 cm时,发声频率范围最小(5~15 kHz),对虾体长为4~5、8~9 cm时,频率范围最大(5~24 kHz)。
图6 对虾摄食虾壳频率范围与对虾体长的关系
Fig.6 Relationship between frequency range of feeding shrimp shell and body length of the shrimp
从图7可见:对虾摄食颗粒饲料的单个脉冲信号时长随体长的增加而增加(拟合方程为y=6.311 1 e0.200 81x,R=0.927 1);对虾体长为9~10 cm时,单个脉冲信号时长最长(约18 ms),明显长于其他体长;对虾体长为2~3 cm时,单个脉冲信号时长最短(约7 ms);对虾体长从4~5变化至8~9 cm时,单个脉冲信号时长无明显变化。
图7 单个脉冲信号时长与对虾体长变化的关系
Fig.7 Relationship between duration of a single pulse signal and body length of the shrimp
对虾摄食虾壳的单个脉冲信号时长随体长的增加而增加(拟合方程为y=16.471 3 e-1.066 5/x,R=0.990 1),对虾体长为2~3 cm时,单个脉冲信号时长最短(约6 ms);对虾体长为9~10 cm时,单个脉冲信号时长最长(约13 ms);对虾体长从6~7 cm变化至9~10 cm时,单个脉冲信号时长随体长增加略有增加。
不同体长对虾摄食颗粒饲料与虾壳的声学参数差异明显。峰值频率方面,对虾体长为2~3 cm时,摄食颗粒饲料的峰值频率较摄食虾壳时明显增加,而对虾体长为8~9、9~10 cm时,摄食颗粒饲料的峰值频率较摄食虾壳时明显降低,对虾体长为4~5、6~7 cm时,摄食颗粒饲料与摄食虾壳的峰值频率差异不大(图4)。在频率范围方面,摄食颗粒饲料的频率范围随对虾体长的增加呈明显缩小的趋势,摄食虾壳的频率范围随对虾体长的增加无明显变化规律(图5、图6)。单个脉冲信号方面,对虾摄食颗粒饲料与虾壳的信号时长均随体长增加而增加,区别在于小体长对虾(2~3、4~5 cm)与大体长对虾(9~10 cm)摄食颗粒饲料的单个脉冲信号时长明显长于摄食虾壳的信号时长(图7)。
对虾摄食的咀嚼机制导致了发声信号特征的不同,不同体长对虾大颚咀嚼的咬合力不同,从而影响其摄食过程中产生信号的声学参数[21]。小体长对虾大颚发育不成熟导致咬合力弱,需要增加单位时间咬合次数以促进颗粒饲料的物理分解,故出现摄食发声频率高而单次摄食发声时间短的特点;随着对虾生长,大颚发育逐渐成熟,咬合力逐渐增强,单位时间内颗粒饲料的物理分解速率提高,因此,咬合次数不断减少,发声频率降低,而单次咬合时间则增加。本研究中,5种不同体长对虾摄食颗粒饲料发声信号的峰值频率和频率范围均随对虾体长的增加而降低,单个脉冲信号时长随对虾体长的增加而增加。Soares等[22]研究表明,大、小2种体长的对虾在平均声能上存在显著差异。Silva等[4]研究表明,3种体长(9~10、11~12、13~14 cm)对虾摄食发声信号的4项声学参数特征无显著性差异。由于其试验对虾体长均较大,未包含对虾从小到大的快速生长时期,不足以概括出对虾摄食声学信号特征随体长变化的规律。Reis等[28]综述了南美白对虾摄食发声信号在养殖中的应用情况,其中归纳了多篇论文中不同体质量(0.04~0.5、10.2、20.5、31.4 g)的对虾摄食声学信号特征,表明不同体质量对虾摄食颗粒饲料发声信号的峰值频率、最低频率和最高频率有差别,即体质量大的对虾发声频率小,体质量小的对虾发声频率大,但并未考察其变化规律。另外,国外研究中所用颗粒饲料与中国的颗粒饲料类型不同,因此,根据中国常用饲料类型及南美白对虾的生长阶段,确定发声信号的变化规律更具有实际应用意义。
此外,不同种类对虾摄食发声信号也有所不同。Berk[29]发现,白滨对虾(Litopenaeus setiferus)摄食发声信号时长约40 ms,是南美白对虾摄食时长(18 ms)的2.7倍;Daniel等[7]发现,斑节对虾(Penaeus monodon)摄食发声频率达到50 kHz,均高于现有相关南美白对虾摄食发声研究结果。不同种类对虾摄食时发声信号存在明显的不同特征,主要在于对虾种类不同、摄食发声机制存在差异。中国对虾养殖种类较多,在某些地区还会采用混养方式,其他种类对虾的摄食发声信号是否存在规律,能否反映对虾生长状态有待进一步研究。
本研究中,5种不同体长对虾摄食虾壳发声信号的峰值频率、频率范围与对虾体长相关,但随着体长增加并无明显变化规律,单个脉冲信号时长则随着体长的增加而增加,4项声学参数特征较对虾摄食颗粒饲料时明显不同。曹正良等[20]研究表明,对虾(体长9~10 cm)摄食沙蚕与不同粒径颗粒饲料时的声学特征存在明显差异,摄食不同粒径颗粒饲料的频率为5~45 kHz,峰值频率为10~17 kHz;摄食沙蚕的频率为5~20 kHz,峰值频率为7 kHz。在虾塘对虾养殖过程中,为了给对虾提供更充足的营养,往往在育苗或特殊情况下补充沙蚕等生物饵料,而对虾生长过程也会出现摄食虾壳的行为。本研究中,对虾(9~10 cm)摄食虾壳的频率范围为10~24 kHz,峰值频率约14 kHz。可以发现,南美白对虾摄食颗粒饲料、沙蚕和虾壳时的声学信号特征存在显著性差异。差异原因可能是三者质地结构不同,沙蚕、虾壳较颗粒饲料坚硬且不易被水浸泡变软,会对南美白对虾摄食沙蚕、虾壳时步足的抓取及大颚的咬合产生一定影响。对虾通过蜕壳实现生长,对虾摄食虾壳发声蕴含着对虾蜕壳后的生长状态,通过研究对虾摄食虾壳发声信号特征可以了解对虾蜕壳生长及异常摄食情况,根据对虾的生长状态可及时更改饲料投喂策略。
曹正良等[20]研究表明,虾塘中检测出摄食信号大约在投饵后的2 min甚至更长时间,而实验室大约在投饵后的10 s,差异原因可能与对虾识别饲料并到达饲料处所需要的时间有关。本研究中,对虾被放置在空间较小的水箱中,因此,饥饿状态的对虾确定饲料较快,所得到的对虾摄食发声信号频率与虾塘具有较高的相似性,虾塘中对虾摄食发声信号的峰值频率为5~15 kHz[5],与本研究中体长为2~3 cm对虾摄食颗粒饲料的峰值频率不一致,可能与虾塘中对虾体长较大且大小分布不均、存在其他水生生物摄食发声等有关。而本研究中单个脉冲信号时长与虾塘中的结果存在显著差异,可能是受到试验环境的影响,虾塘声音环境复杂,受多种其他水生生物摄食饲料的声音影响,同时试验玻璃水箱存在声音反射等因素。因此,在后续虾塘摄食发声试验中,要更多考虑对虾摄食虾塘残余饲料和其他饲料(鱼、虾等),以及虾塘中其他水生生物及群虾的摄食发声。
1)对虾摄食颗粒饲料时,随着对虾体长的增加,摄食发声信号的峰值频率降低,而频率范围和单个脉冲信号时长均增加。这些变化呈现一定的线性变化规律。
2)对虾摄食虾壳时,随着对虾体长的增加,摄食发声信号的峰值频率呈先降低后上升,频率范围呈先增加后降低再增加再降低的波动变化,单个脉冲信号时长则增加。与摄食颗粒饲料相比,这些摄食发声信号无明显的线性变化规律。
3)体长2~3 cm的对虾摄食颗粒饲料时峰值频率(25 kHz)显著高于摄食虾壳时的峰值频率(13 kHz);体长8~9 cm的对虾摄食颗粒饲料的峰值频率(8 kHz)较摄食虾壳(16 kHz)时明显降低;体长9~10 cm的对虾摄食颗粒饲料的单个脉冲信号时长(18 ms)明显长于摄食虾壳的信号时长(13 ms)。说明对虾摄食颗粒饲料与虾壳的发声信号特征明显不同。
[1] BONDAD-REANTASO M G,SUBASINGHE R P,JOSUPEIT H,et al.The role of crustacean fisheries and aquaculture in global food security:past,present and future[J].Journal of Invertebrate Pathology,2012,110(2):158-165.
[2] 农业农村部渔业渔政管理局编制.中国渔业统计年鉴2017[M].北京:中国农业出版社,2017.
Ministry of Agriculture and Rural Affairs of the People’s Republic of CHINA.China fishery statistical yearbook[M].Beijing:China Agriculture Press,2017.(in Chinese)
[3] CASILLAS-HERNNDEZ R,MAGALLN-BARAJAS F,PORTILLO-CLARCK G,et al.Nutrient mass balances in semi-intensive shrimp ponds from Sonora,Mexico using two feeding strategies:trays and mechanical dispersal[J].Aquaculture,2006,258(1/2/3/4):289-298.
[4] SILVA J F,HAMILTON S,ROCHA J V,et al.Acoustic characterization of feeding activity of Litopenaeus vannamei in captivity[J].Aquaculture,2019,501:76-81.
[5] 张丽珍,陆天辰,杨加庆,等.南美白对虾进食前后发声信号特点分析[J].声学技术,2020,39(4):413-418.
ZHANG L Z,LU T C,YANG J Q,et al.Analysis on sound characteristics of Penaeus vannamei before and after feeding[J].Technical Acoustics,2020,39(4):413-418.(in Chinese)
[6] BERK I M,EVANS W E,BENSON R H,et al.The use of passive sonar to detect sound production and calculate population densities of penaeid shrimp in the Gulf of Mexico[J].The Journal of the Acoustical Society of America,1996,99(4):2533-2574.
[7] SMITH D V,TABRETT S.The use of passive acoustics to measure feed consumption by Penaeus monodon (giant tiger prawn) in cultured systems[J].Aquacultural Engineering,2013,57:38-47.
[8] LAMMERS M O,MUNGER L M.From shrimp to whales:biological applications of passive acoustic monitoring on a remote Pacific coral reef[M]//Modern Acoustics and Signal Processing.New York:Springer,2016:61-81.
[9] SMITH D V,SHAHRIAR M S.A context aware sound classifier applied to prawn feed monitoring and energy disaggregation[J].Knowledge-Based Systems,2013,52:21-31.
[10] 曹正良,张丽珍,胡庆松.南美白对虾养殖虾塘中水下声景的应用研究初探[C]//中国声学学会水声学分会2019年学术会议论文集.南京:中国声学学会水声学分会,2019:346-348.
CAO Z L,ZHANG L Z,HU Q S.Preliminary study on the application of underwater soundscape in the farming shrimp pond of Penaeus vannamei[C]//Proceedings of the 2019 Academic Conference of the Underwater Acoustics Branch of the Acoustical Society of China.Nanjing:Academic Conference of the Underwater Acoustics Branch of the Acoustical Society of China,2019:346-348.(in Chinese)
[11] JESCOVITCH L N,ULLMAN C,RHODES M,et al.Effects of different feed management treatments on water quality for Pacific white shrimp Litopenaeus vannamei[J].Aquaculture Research,2018,49(1):526-531.
[12] ULLMAN C,RHODES M,HANSON T,et al.Effects of four different feeding techniques on the pond culture of Pacific white shrimp,Litopenaeus vannamei[J].Journal of the World Aquaculture Society,2019,50(1):54-64.
[13] WEI M C,CHEN K Y,LIN Y T,et al.Recognition of behavior state of Penaeus vannamei based on passive acoustic technology[J].Frontiers in Marine Science,2022,9:973284.
[14] DE LIMA P P,PONTES C S,ARRUDA M D F.Activity pattern of the marine shrimp Litopenaeus vannamei (Boone 1931) in laboratory as a function of different feeding frequencies[J].Aquaculture Research,2009,41(1):53-60.
[15] SILVA P F,DE SOUZA MEDEIROS M,ALVES SILVA H P,et al.A study of feeding in the shrimp Farfantepenaeus subtilis indicates the value of species level behavioral data for optimizing culture management[J].Marine and Freshwater Behaviour and Physiology,2012,45(2):121-134.
[16] BARDERA G,USMAN N,OWEN M,et al.The importance of behaviour in improving the production of shrimp in aquaculture[J].Reviews in Aquaculture,2019,11(4):1104-1132.
[17] SANTIAG O,HAMILTO N.The loud crowd:interactions between stocking density and acoustic feeding activity of different size classes of Litopenaeus vannamei[J].Aquaculture,2023,563:738904.
[18] PEIXOTO S,SOARES R,ALLEN DAVIS D.An acoustic based approach to evaluate the effect of different diet lengths on feeding behavior of Litopenaeus vannamei[J].Aquacultural Engineering,2020,91:102114.
[19] SILVI O,PEIXOT O.Acoustic activity of Litopenaeus vannamei fed pelleted and extruded diets[J].Aquaculture,2020,525:735307.
[20] 曹正良,沈梦庭,李钊丞,等.摄食不同粒径颗粒饲料的凡纳滨对虾发声信号特征[J].南方水产科学,2022,18(6):26-34.
CAO Z L,SHEN M T,LI Z C,et al.Characteristics of feeding acoustic signals of Litopenaeus vannamei fed with pellets of different sizes[J].South China Fisheries Science,2022,18(6):26-34.(in Chinese)
[21] BOHNENSTIEHL D R,LILLIS A,EGGLESTON D B.The curious acoustic behavior of estuarine snapping shrimp:temporal patterns of snapping shrimp sound in sub-tidal oyster reef habitat[J].PLoS One,2016,11(1):e0143691.
[22] SOARES R,PEIXOTO S,GALKANDA-ARACHCHIGE H S C,et al.Growth performance and acoustic feeding behavior of two size classes of Litopenaeus vannamei fed pelleted and extruded diets[J].Aquaculture International,2021,29(1):399-415.
[23] 申玉春,陈作洲,刘丽,等.盐度和营养对凡纳滨对虾蜕壳和生长的影响[J].水产学报,2012,36(2):290-299.
SHEN Y C,CHEN Z Z,LIU L,et al.The effects of salinity and nutrition on molt and growth of Litopenaeus vannamei[J].Journal of Fisheries of China,2012,36(2):290-299.(in Chinese)
[24] DALL W,HILL B J,ROTHLISBERG P C,et al.The biology of the Penaeidae[M]//Advances in Marine Biology.New York:Academic Press,1990:1-489.
[25] MERCHANT N D,FRISTRUP K M,JOHNSON M P,et al.Measuring acoustic habitats[J].Methods in Ecology and Evolution,2015,6(3):257-265.
[26] 王岩,隋思涟.试验设计与MATLAB数据分析[M].北京:清华大学出版社,2012.
WANG Y,SUI S L.Experimental design and MATLAB data analysis[M].Beijing:Tsinghua University Press,2012.(in Chinese)
[27] LILLIEFORS H W.On the Kolmogorov-Smirnov test for the exponential distribution with mean unknown[J].Journal of the American Statistical Association,1969,64(325):387-389.
[28] REIS J,PEIXOTO S,SOARES R,Passive acoustic monitoring as a tool to assess feed response and growth of shrimp in ponds and research systems[J].Aquaculture,2022,546:737326.
[29] BERK I M.Sound production by white shrimp (Penaeus setiferus),analysis of another crustacean-like sound from the Gulf of Mexico,and applications for passive sonar in the shrimping industry[J].Journal of Shellfish Research,1998,17(5):1497-1500.