南美白對(duì)蝦品系、密度、放苗時(shí)間與養(yǎng)殖效益的統(tǒng)計(jì)學(xué)分析

李舉鵬，史健，胡小弘，周影，王彩蘊(yùn)，單紅云，孫麗敏

(1. 營(yíng)口市水產(chǎn)科學(xué)研究所,遼寧 營(yíng)口 115003；2. 老邊區(qū)水產(chǎn)技術(shù)推廣站,遼寧 營(yíng)口 115005；3. 營(yíng)口市水產(chǎn)技術(shù)推廣站,遼寧 營(yíng)口 115003)

南美白對(duì)蝦是世界養(yǎng)殖產(chǎn)量最高的三大優(yōu)良蝦種之一[1]，自1988年引入我國(guó)[2]。目前，我國(guó)南美白對(duì)蝦的養(yǎng)殖模式多種多樣，不同的養(yǎng)殖模式對(duì)蝦苗的品系和要求也在不斷變化，但隨著養(yǎng)殖時(shí)間的推移，病害加劇，養(yǎng)殖難度愈來(lái)愈大，養(yǎng)殖成功率很低，國(guó)內(nèi)很多學(xué)者在養(yǎng)殖池塘水質(zhì)的理化因子[3]、浮游植物群落結(jié)構(gòu)[4]， 水體微生物群落結(jié)構(gòu)[5]以及對(duì)蝦餌料[6]等方面，進(jìn)行了大量的研究，這些研究往往是從試驗(yàn)對(duì)比的角度探求對(duì)蝦發(fā)病的原因，該文擬從養(yǎng)殖生產(chǎn)的實(shí)際效果調(diào)查入手，通過(guò)統(tǒng)計(jì)分析的方法，以尋求對(duì)蝦發(fā)病的原因， 為此，作者于2017年在營(yíng)口市高坎鎮(zhèn)太平村育苗場(chǎng)進(jìn)行了不同品系蝦苗的標(biāo)粗淡化，然后將這些淡化蝦苗投放到不同的養(yǎng)殖區(qū)域進(jìn)行淡水養(yǎng)殖，并于2018年在浙江寧海華大海昌公司進(jìn)行了同一品系(SIS)蝦苗不同密度及不同季節(jié)放苗的養(yǎng)殖生產(chǎn)，通過(guò)跟蹤調(diào)查，回訪等方式，將收集到的數(shù)據(jù)整理，并進(jìn)行統(tǒng)計(jì)學(xué)分析，以探求南美白對(duì)蝦養(yǎng)殖失敗的原因.

1 方法

1.1 不同品系蝦苗的養(yǎng)殖

2017年4月25日至5月21日，先后從南方育苗場(chǎng)購(gòu)進(jìn)科納灣、SIS和正大三個(gè)品系，共11個(gè)批次的P4-P5一代對(duì)蝦苗8 000萬(wàn)尾，經(jīng)過(guò)12～19 d的培育，共培育出體長(zhǎng)1.0～1.5 cm的仔蝦苗6 250萬(wàn)尾，其中科納灣蝦苗3 830萬(wàn)尾、SIS蝦苗1 640萬(wàn)尾、正大蝦苗780萬(wàn)尾，然后分別投放到營(yíng)口市的水源鎮(zhèn)、老邊區(qū)、溝沿和高坎鎮(zhèn)的土池進(jìn)行淡水養(yǎng)殖，每個(gè)池塘的面積為0.067～0.670 hm2，放苗密度為4萬(wàn)尾/667 m2左右，養(yǎng)殖總面積約100 hm2，養(yǎng)殖池塘共327個(gè)。盡管各批次蝦苗剛?cè)氤氐瘯r(shí)規(guī)格及活力有所不同，但經(jīng)過(guò)十多天的淡化培育，放苗時(shí)蝦苗的活力以及攝食情況等蝦苗質(zhì)量方面基本相差無(wú)幾。

1.2 單一品系蝦苗的養(yǎng)殖

2018年4月20日開始至2018年8月9日先后24次從海南育苗場(chǎng)購(gòu)進(jìn)南美白對(duì)蝦(SIS)無(wú)節(jié)幼體5.2億尾，共培育出體長(zhǎng)0.8～1.2 cm的仔蝦1.8億尾，分別投放到浙江、江蘇、河南、山東等地的池塘進(jìn)行養(yǎng)殖。

1.3 數(shù)據(jù)收集

數(shù)據(jù)收集主要采取跟蹤調(diào)查和走訪的方式進(jìn)行，每個(gè)養(yǎng)殖池塘作為一個(gè)觀察單位，觀察時(shí)間為90 d，分別對(duì)三個(gè)品系蝦苗發(fā)病時(shí)的養(yǎng)殖天數(shù)和養(yǎng)殖效益進(jìn)行記錄，并將收集來(lái)的數(shù)據(jù)進(jìn)行匯總整理。

1.4 統(tǒng)計(jì)分析

統(tǒng)計(jì)分析采用SPSS.17.0軟件進(jìn)行，分別對(duì)對(duì)蝦養(yǎng)殖期間發(fā)病時(shí)的存活天數(shù)進(jìn)行秩和檢驗(yàn)、對(duì)養(yǎng)殖效益進(jìn)行χ2檢驗(yàn)，不同養(yǎng)殖密度和不同放苗季節(jié)對(duì)養(yǎng)殖效益的影響采取圖表分析，統(tǒng)計(jì)學(xué)差異水平規(guī)定為P<0.05。

2 結(jié)果

2.1 秩和檢驗(yàn)

Kruskal-Wallis檢驗(yàn)是1952年由Kruskal和Wallis二人共同提出的，是二個(gè)獨(dú)立樣本檢驗(yàn)的一種推廣[7]，Kruskal-Wallis檢驗(yàn)用于推斷計(jì)量資料或等級(jí)資料的多個(gè)獨(dú)立樣本所來(lái)自的多個(gè)總體分布是否有差別[8]，其基本思想是：將原始資料在不分組的情況下從小到大排順序，即編秩次，然后分組將秩次相加得到各組的秩和，若組間秩和相差不大，則認(rèn)為各組間無(wú)差異。

其計(jì)算公式為：

(1)

式中：H——統(tǒng)計(jì)量；

Ri——各組的秩和；

ni——各組對(duì)應(yīng)的例數(shù)；

N=∑ni。

2.2 同一品系蝦苗在不同區(qū)域養(yǎng)殖 蝦苗的發(fā)病情況

高坎鎮(zhèn)和溝沿鎮(zhèn)由于養(yǎng)殖面積相對(duì)較少，所以合并為一個(gè)區(qū)域，每個(gè)池塘蝦苗的存活時(shí)間按照對(duì)蝦發(fā)病時(shí)的日期與放苗日期之間的天數(shù)之差計(jì)算。各養(yǎng)殖區(qū)以養(yǎng)殖科納灣蝦苗的發(fā)病池塘為選取條件；老邊區(qū)發(fā)病池塘較多，從中隨機(jī)選取30個(gè)發(fā)病池塘做為樣本，各發(fā)病池塘蝦苗的存活時(shí)間如表1。

經(jīng)SPSS.17.0計(jì)算得H=6.250，得P<0.05， 差別有統(tǒng)計(jì)學(xué)意義，其中，均值為52.95 d，標(biāo)準(zhǔn)差為13.215，即：蝦池平均發(fā)病時(shí)間是53 d，三個(gè)區(qū)域中蝦池最早發(fā)病時(shí)間是25 d，最晚發(fā)病時(shí)間是80 d，說(shuō)明水源鎮(zhèn)、老邊區(qū)、高坎和溝沿鎮(zhèn)蝦苗的發(fā)病時(shí)間是有差異的。

2.3 不同品系蝦苗在同一地區(qū)養(yǎng)殖 蝦苗的發(fā)病情況

老邊區(qū)養(yǎng)殖面積相對(duì)較大，發(fā)病比例也最高，在此區(qū)域內(nèi)，每個(gè)品系的蝦苗隨機(jī)抽取15個(gè)發(fā)病池塘，每個(gè)池塘蝦苗的存活時(shí)間也是按照對(duì)蝦發(fā)病時(shí)的日期與放苗日期之間的天數(shù)之差計(jì)算，各品系蝦苗存活時(shí)間如表2。

表2 同一區(qū)域不同品系蝦苗養(yǎng)殖存活天數(shù)

經(jīng)SPSS.17.0計(jì)算得H=1.527，P>0.05，差別沒(méi)有統(tǒng)計(jì)學(xué)意義，其中，均值為59.69 d，標(biāo)準(zhǔn)差為12.989，說(shuō)明三個(gè)品系的蝦苗在老邊區(qū)養(yǎng)殖平均的發(fā)病時(shí)間是57 d，最早發(fā)病時(shí)間是29 d，最晚是80 d，但科納灣、SIS和正大3個(gè)品系的蝦苗發(fā)病時(shí)間總體分布是沒(méi)有差異的。

2.4 χ2檢驗(yàn)

χ2檢驗(yàn)是英國(guó)統(tǒng)計(jì)學(xué)家Pearson于1900年提出的一種用途很廣的計(jì)數(shù)資料的假設(shè)檢驗(yàn)方法[9]。屬于非參數(shù)檢驗(yàn)，主要是比較兩個(gè)及兩個(gè)以上樣本率(構(gòu)成比)以及兩個(gè)分類變量的關(guān)聯(lián)性分析。根本思想在于比較理論頻數(shù)和實(shí)際頻數(shù)的吻合程度或者擬合優(yōu)度，其計(jì)算公式為：

(2)

式中：A——實(shí)際頻數(shù)；

T——理論頻數(shù)。

2.5 不同區(qū)域、不同品系對(duì)蝦養(yǎng)殖效益情況

對(duì)蝦養(yǎng)殖的最終目的是為了獲得經(jīng)濟(jì)效益，將不同區(qū)域，不同品系蝦苗的養(yǎng)殖效益進(jìn)行匯總分析，更有利于揭示它們之間的關(guān)系，該次共調(diào)查了327個(gè)養(yǎng)殖池塘，其中，老邊區(qū)發(fā)病嚴(yán)重，調(diào)查了144個(gè)養(yǎng)蝦池塘?？萍{灣蝦苗在各區(qū)域養(yǎng)殖的效益情況見表3，不同品系蝦苗在同一區(qū)域養(yǎng)殖的效益情況見表4.

表3 科納灣蝦苗在不同區(qū)域的養(yǎng)殖效益情況(以池塘個(gè)數(shù)計(jì))

經(jīng)SPSS.17.0計(jì)算得χ2=31.792，P<0.05， 差別有統(tǒng)計(jì)學(xué)意義，從統(tǒng)計(jì)結(jié)果來(lái)看，高坎和溝沿鎮(zhèn)養(yǎng)殖效益最好，老邊區(qū)養(yǎng)殖效果最差，水源鎮(zhèn)介于二者之間，3個(gè)區(qū)域養(yǎng)殖效益差異特別明顯。

表4 不同品系蝦苗在老邊區(qū)的養(yǎng)殖效益情況(以池塘個(gè)數(shù)計(jì))

經(jīng)SPSS.17.0計(jì)算得χ2=1.892,P>0.05， 差別沒(méi)有統(tǒng)計(jì)學(xué)意義，統(tǒng)計(jì)結(jié)果說(shuō)明3個(gè)品系的蝦苗在同一區(qū)域養(yǎng)殖，就總體而言，其養(yǎng)殖效益差異不明顯。

2.6 不同養(yǎng)殖密度下池塘的發(fā)病情況

2018年，由于放苗面積較廣，路途較遠(yuǎn)，有些調(diào)查內(nèi)容難以獲得，因此只對(duì)各批次的蝦苗是否發(fā)病及養(yǎng)殖密度進(jìn)行了簡(jiǎn)單的調(diào)查，養(yǎng)殖密度分為4個(gè)梯度，即：1～2萬(wàn)尾/667 m2、2～3萬(wàn)尾/667 m2、3～4萬(wàn)尾/667 m2、4～6萬(wàn)尾/667 m2。池塘蝦苗發(fā)病比率分別是：38%，41%，47%，68%，如圖1所示。

圖1 不同養(yǎng)殖密度下發(fā)病池塘所占比率

結(jié)果顯示，隨著養(yǎng)殖密度的增大，蝦池的發(fā)病比例也在隨之增加，尤其高溫時(shí)期放苗，當(dāng)養(yǎng)殖密度達(dá)到4萬(wàn)尾以上/667 m2時(shí)，土池養(yǎng)殖的發(fā)病比率非常高。

2.7 不同放苗時(shí)間，池塘的發(fā)病情況

今年從4月20日開始到8月9日結(jié)束，先后24批次，共培育出0.8～1.2 cm的仔蝦1.8億尾，分別投放到浙江、江蘇、、河南、山東等地的池塘進(jìn)行養(yǎng)殖，放苗結(jié)束后，隨時(shí)收集記錄各批次蝦苗的發(fā)病情況，根據(jù)各批次蝦苗的發(fā)病情況，計(jì)算其發(fā)病比率，繪制成圖，各批次蝦苗發(fā)病比率如圖2。

圖2 不同批次蝦苗發(fā)病情況

不同季節(jié)放苗，各批次蝦苗的發(fā)病比率是不一樣的，最高為93%，最低為18%，4月的幼體苗，發(fā)病池塘比例在78%～83%，5月和6月的幼體苗，發(fā)病池塘比例較低，發(fā)病池塘比例維持在18%～3%，7月以后的幼體苗，發(fā)病蝦池比例為64%～93%，而且，高溫期放苗往往伴隨著蝦苗的生長(zhǎng)緩慢和參差不齊。

3 分析與討論

3.1 同一品系蝦苗在不同區(qū)域養(yǎng)殖蝦苗的發(fā)病情況

同一品系的蝦苗投放到不同的養(yǎng)殖區(qū)域進(jìn)行養(yǎng)殖，發(fā)病時(shí)間是存在差異的，從表1的統(tǒng)計(jì)結(jié)果可以看出，3個(gè)品系的蝦苗在不同養(yǎng)殖區(qū)域的發(fā)病時(shí)間具有明顯的差異，這與養(yǎng)殖的實(shí)際情況是吻合的。秦志浩等[10]在對(duì)不同品系南美白對(duì)蝦室內(nèi)工廠化養(yǎng)殖的研究中認(rèn)為：我國(guó)南美白對(duì)蝦的養(yǎng)殖模式多種多樣，不同養(yǎng)殖模式對(duì)蝦苗的品系和要求也在不斷變化，在工廠化養(yǎng)殖模式下，正大品系蝦苗最適合養(yǎng)殖，產(chǎn)量較高，個(gè)體較大；但楊文慧等[11]曾對(duì)天津地區(qū)南美白對(duì)蝦產(chǎn)業(yè)的調(diào)查分析指出：南美白對(duì)蝦養(yǎng)殖過(guò)程中，確實(shí)存在地區(qū)化差異和年份差異，使得南美白對(duì)蝦在不同區(qū)域養(yǎng)殖效果表現(xiàn)不一。該次調(diào)查的結(jié)果與楊文慧等人的觀點(diǎn)是一致的，從營(yíng)口地區(qū)的養(yǎng)蝦實(shí)際情況來(lái)看，對(duì)蝦養(yǎng)殖的發(fā)病時(shí)間確實(shí)具有區(qū)域性的特點(diǎn)，即：在某一時(shí)間段同一區(qū)域的絕大部分養(yǎng)殖池塘同時(shí)或相繼發(fā)病，這很可能與養(yǎng)殖區(qū)域的水質(zhì)以及氣候狀況有關(guān)。

3.2 不同品系蝦苗在同一地區(qū)養(yǎng)殖蝦苗的發(fā)病情況

很多養(yǎng)殖戶認(rèn)為，對(duì)蝦的發(fā)病與蝦苗的品系有關(guān)，對(duì)蝦養(yǎng)殖期間的發(fā)病往往都是由于蝦苗品系不好造成的，但從表2的統(tǒng)計(jì)結(jié)果來(lái)看，3個(gè)品牌蝦苗在同一區(qū)域養(yǎng)殖，其發(fā)病時(shí)間是沒(méi)有統(tǒng)計(jì)學(xué)差異的。張嘉晨等[12]在進(jìn)行高密度養(yǎng)殖凡納濱對(duì)蝦(Litopenaeusvannamei)新品系選育時(shí)的研究認(rèn)為：以不同國(guó)家和地區(qū)的8個(gè)凡納濱對(duì)蝦養(yǎng)殖群體作為親本，構(gòu)建8個(gè)子代群體，美國(guó)凡納濱對(duì)蝦改良公司(SIS)♀×SIS♂生長(zhǎng)優(yōu)勢(shì)明顯，生長(zhǎng)比速較群體均值提高10.4%～17.9%;海大選育(HD)♀×粵海選育(YH)♂成活優(yōu)勢(shì)明顯，成活率較群體均值提高9.3%～13.7%;YH♀×中興1號(hào)(ZX)♂和YH♀×科納灣種蝦(K)♂生長(zhǎng)比速、成活率均高于群體均值，表現(xiàn)出生長(zhǎng)快成活率高雙重優(yōu)良性狀，該次調(diào)查的結(jié)論與上述觀點(diǎn)不一致，是否與各地的養(yǎng)殖環(huán)境不同有關(guān)，還有待于進(jìn)一步的研究。

3.3 不同區(qū)域、不同品系蝦苗養(yǎng)殖效益分析

從表3、表4的分析結(jié)果看出，同一品系蝦苗在不同區(qū)域養(yǎng)殖，蝦苗的發(fā)病時(shí)間存在差異，導(dǎo)致養(yǎng)殖效益也存在差異；不同品系蝦苗在同一區(qū)域養(yǎng)殖，養(yǎng)殖效益沒(méi)有顯著性差異。陳曉明等[13]在進(jìn)行微生態(tài)制劑防控南美白對(duì)蝦病害試驗(yàn)發(fā)現(xiàn)：經(jīng)濟(jì)效益的高低與放苗的品種有直接關(guān)系，即：一代苗與土苗對(duì)比，養(yǎng)殖效益差異特別明顯；但不同品系的一代苗之間(SIS和桂海)，養(yǎng)殖效益差異不大，這一結(jié)論與該次調(diào)查結(jié)果基本是吻合的，基于這一結(jié)論，建議在對(duì)蝦養(yǎng)殖的實(shí)際生產(chǎn)中著重注意水環(huán)境的改善以及養(yǎng)殖技術(shù)的提高。

3.4 不同養(yǎng)殖密度下池塘的發(fā)病情況分析

由圖1可以看出，隨著養(yǎng)殖密度的增加，發(fā)病池塘所占的比例越來(lái)越大，在養(yǎng)殖過(guò)程中，特別是高溫季節(jié)，隨著個(gè)體的增大，尤其對(duì)蝦體長(zhǎng)達(dá)到6～7 cm以后，養(yǎng)殖水體中的殘餌、糞便大量增加，水環(huán)境嚴(yán)重惡化，加速了對(duì)蝦發(fā)病的幾率。李玉全等[16]認(rèn)為：凡納濱對(duì)蝦(體長(zhǎng)<7.6 cm或體質(zhì)量<6.1 g)養(yǎng)殖密度為150～900尾/m3時(shí)，造成對(duì)蝦生長(zhǎng)和非特異性免疫因子差異的主要原因是密度脅迫，水質(zhì)因子的作用是次要的。張國(guó)新[15]曾就不同養(yǎng)殖密度對(duì)南美白對(duì)蝦生長(zhǎng)的影響做過(guò)試驗(yàn)，試驗(yàn)設(shè)置養(yǎng)殖密度分別為75尾/m2、150尾/m2、300尾/m2、450尾/m2、900尾/m25種梯度，結(jié)果表明：不同密度對(duì)南美白對(duì)蝦的成活率有顯著影響，各不同試驗(yàn)組的成活率分別為93.7%、94.7%、84.7%、77%、69%，隨密度增加，成活率呈下降趨勢(shì)，各組間成活率最大相差25.7%[14]。張呂平等報(bào)道了在我國(guó)熱帶地區(qū)集約化養(yǎng)殖凡納濱對(duì)蝦(L.vannamei)4種密度下(70、100、120、150尾/m)養(yǎng)殖成活率分別為80.44%、75.16%、51.74%和44.43%。上述觀點(diǎn)在實(shí)際養(yǎng)殖生產(chǎn)中也得到了印證，大面積低密度的大圍池塘養(yǎng)殖成功率明顯高于小面積的精養(yǎng)池塘。

3.5 不同放苗時(shí)間，池塘的發(fā)病情況分析

由圖2可以看出，同一品系的蝦苗，5、6月放養(yǎng)的各批苗，發(fā)病比率比較低，4月和7月以后放養(yǎng)的蝦苗，發(fā)病比例都比較高，主要是由于前期氣溫較低，溫差變化較大，7月進(jìn)入高溫期以后，由于受氣候和天氣狀況的影響，水中的理化因子及藻類密度、結(jié)構(gòu)變化劇烈，導(dǎo)致水質(zhì)不穩(wěn)定，從而造成蝦病的發(fā)生。彭榮南等[17]在臺(tái)風(fēng)引起養(yǎng)殖水質(zhì)突變的研究中發(fā)現(xiàn)：風(fēng)切變引起水環(huán)流擾亂，有毒物質(zhì)擴(kuò)散，大風(fēng)引起對(duì)蝦過(guò)度應(yīng)激反應(yīng);同時(shí)，伴隨陰雨寡照使得水體光合作用減少(DO減少)，DO達(dá)到最低，對(duì)蝦病害明顯；沈明明[18]在對(duì)蝦工廠化養(yǎng)殖水體微藻群落演替的主要因子及其與疾病發(fā)生的關(guān)系中指出：微藻的密度和生物量顯著影響對(duì)蝦發(fā)病程度，過(guò)高的微藻密度尤其是養(yǎng)殖后期藍(lán)藻優(yōu)勢(shì)種的密度，與養(yǎng)殖后期對(duì)蝦持續(xù)發(fā)病關(guān)系密切；過(guò)低的微藻密度和生物量則會(huì)造成微藻對(duì)蝦池生態(tài)系統(tǒng)調(diào)控能力低下，導(dǎo)致對(duì)蝦生長(zhǎng)緩慢，個(gè)頭參差不齊。

對(duì)蝦病害的發(fā)生及養(yǎng)殖效益，受多種因素的影響，其中，養(yǎng)殖環(huán)境的污染導(dǎo)致病害的發(fā)生，國(guó)內(nèi)外很多學(xué)者、研究人員進(jìn)行了大量的研究[19-21]，從該次調(diào)查的結(jié)果來(lái)看，蝦苗品系不是影響?zhàn)B殖效益的根本因素，養(yǎng)殖密度、放苗時(shí)間、養(yǎng)殖水環(huán)境以及天氣等因素都會(huì)導(dǎo)致對(duì)蝦病害的發(fā)生，從而影響?zhàn)B殖效益。