大連海域黃鮟鱇的生長特征與資源合理利用

尹增強(qiáng)　許傳才　陳勇

摘要[目的]制定大連海域黃鮟鱇開捕規(guī)格和總允許捕撈量。[方法]通過分析黃鮟鱇的體長—體重關(guān)系式和體長生長方程、體重生長方程、瞬時(shí)死亡率，計(jì)算體重生長的拐點(diǎn)年齡和臨界年齡，估算大連海域黃鮟鱇的資源量。[結(jié)果]體長生長方程為Lt=576.4×[1-e-0.28（t+0.44）]，體重生長方程為Wt=5 344.3×[1-e-0.28（t+0.44）]3.005 2，黃鮟鱇體重生長的拐點(diǎn)年齡為4.2齡，臨界年齡為4.7齡，大連海域黃鮟鱇的資源生物量為8 216.7 t。[結(jié)論]黃鮟鱇開捕體長419.2～439.7 mm，開捕體重2 052.2～2 369.5 g，總允許捕撈量1 182 t。

關(guān)鍵詞黃鮟鱇；生長方程；臨界年齡；資源生物量；大連海域

中圖分類號(hào)S932文獻(xiàn)標(biāo)識(shí)碼

A文章編號(hào)0517-6611（2015）28-125-03

Growth and Rational Utilization of Lophius litulon in Dalian Sea

YIN Zengqiang， XU Chuancai， CHEN Yong （Marine Science and Environment College， Dalian Fisheries University， Dalian， Liaoning 116023）

Abstract [Objective] To formulate Lophius litulon catchable size and total allowable catch in Dalian sea. [Methods] Based on analyzing the lengthweight formula， body length and weight growth equation， and instantaneous mortality rate， to calculate weight growth inflexion age， the critical age and the biomass of Lophius litulon in Dalian sea. [Results] The growth equation of length and weight are respectively Lt=576.4×[1-e-0.28（t+0.44）] and Wt=5 344.3×[1-e-0.28（t+0.44）]3.005 2， and the weight growth inflection age and critical age of Lophius litulon are 4.2 age and 4.7 age. The biomass is about 8 216.7 t. [Conclusion] In terms of Lophius litulon stock in Dalian sea， the catchable length is about 419.2-439.7 mm， and catchable weight is about 2 052.2-2 369.5 g. Total allowable catch is about 1 182 t.

Key words Lophius litulon； VBGF； Critical age； Biomass； Dalian sea

黃鮟鱇（Lophius litulon）隸屬鮟鱇目（Lophiiformes）、鮟鱇科（Lophiidae）黃鮟鱇屬F（Lophius），俗稱海蛤蟆，為近海底層魚類，主要分布于西北太平洋和印度洋，在我國渤海、黃海、東海和南海均有產(chǎn)出[1]。黃鮟鱇肉質(zhì)鮮美，干制品甚受歡迎，除具有較高的食用價(jià)值外，還具有較高的藥用價(jià)值，在日韓等國外市場(chǎng)非常走俏，市場(chǎng)潛力巨大，是重要海產(chǎn)經(jīng)濟(jì)種類[2]。目前，有關(guān)黃鮟鱇資源研究主要集中于黃海南部和東海資源群體的研究[2-4]，針對(duì)大連海域黃鮟鱇研究的文章尚未報(bào)道。根據(jù)大連海域黃鮟鱇的調(diào)查資料，筆者對(duì)黃鮟鱇生長特征和資源管理進(jìn)行研究，希望對(duì)大連海域黃鮟鱇資源的合理利用有所裨益。

1材料與方法

1.1數(shù)據(jù)來源

2010年10月、2013年10月和2014年10月分別在大連大黑石碼頭、黑石礁碼頭和孫家溝水產(chǎn)品市場(chǎng)隨機(jī)取樣，將樣品帶到實(shí)驗(yàn)室進(jìn)行生物學(xué)測(cè)定（測(cè)定項(xiàng)目包括體長、體重和性別鑒定等），3次共采樣53條黃鮟鱇。

1.2方法

（1）分析黃鮟鱇的平均體長、平均體重、體長組成和體重組成。

（2）求解體長與體重關(guān)系式：黃鮟鱇體長和體重的關(guān)系呈冪函數(shù)相關(guān)[5]，其表達(dá)式為W=a·Lb（式中，W為體重，L為體長，a、b為參數(shù)）。

（3）求解體長生長方程和體重生長方程：黃鮟鱇生長方程模型選用Von Bertalanffy生長方程。對(duì)最大漸進(jìn)體長L∞和生長曲率K的估計(jì)采用FISATⅡ中的體長頻率分析法ELEFAN I（Electronic length frequency analysis I）估算；理論初始年齡t0則采用Pauly經(jīng)驗(yàn)公式[6]計(jì)算，ln（-t0）=-0.392 2-0.275 2lnL∞- 1.038lnK；最大漸進(jìn)體重W∞由公式W∞=a×L∞b求得。

（4） 求解黃鮟鱇的資源量：采用體長股分析法求解，主要公式[7]為：Nt=Ct·（M+Ft）/Ft；Ci=Ni+Δt·（Fi/（M+Fi）·（e（M+Fi）Δti-1（式中，Δti=（ti+1-ti）；其中，ti=t0-（1/K）·ln（1-（Li/L∞））；Ni=Ni+Δt·eM+Fi。

式中，Nt為黃鮟鱇最大體長時(shí)的資源量，Ct為黃鮟鱇最大體長時(shí)的漁獲量，M為自然死亡系數(shù)（由Pauly[8]經(jīng)驗(yàn)公式：lnM=-0.006 6-0.279lnL∞+0.654 3lnK+0.463 4lnT，T為棲息海域年均水溫，該試驗(yàn)取20 ℃），F(xiàn)t為黃鮟鱇最大體長的捕撈死亡系數(shù)，Ni、Ni+△t為黃鮟鱇i齡和i+△t齡時(shí)的資源量，Ci為黃鮟鱇i齡時(shí)的漁獲量，F(xiàn)i為黃鮟鱇i齡的捕撈死亡系數(shù)。

2結(jié)果與分析

2.1黃鮟鱇的平均體長、平均體重、體長組成和體重組成

大連海域黃鮟鱇群體體長范圍為90～491 mm，均值在220.5 mm左右，優(yōu)勢(shì)體長組為151～200 mm，占調(diào)查黃鮟鱇總數(shù)的28.3%（圖1）； 大連海域黃鮟鱇群體體重范圍為20～3 097 g，均值在540.3 g左右。體重分布主要集中在0～200 g，占群體總數(shù)的52.8%以上，優(yōu)勢(shì)體重組為101～200 g（圖2）。

2.2體長與體重的關(guān)系式

在取樣的53尾黃鮟鱇中，能夠鑒別性別的尾數(shù)為49尾。協(xié)方差分析（ANCOVA）表明，雌性和雄性黃鮟鱇的體長與體重對(duì)應(yīng)關(guān)系無明顯差異（P>0.05），因此，體長與體重的關(guān)系按雌、雄合并擬合（圖3），其表達(dá)式為 W=2.7×10-5L3.005 2（R2=0.993 3 ，n=53）。

2.3體長與體重生長

2.3.1體長生長方程與體重生長方程。

根據(jù)的大連海域黃鮟鱇的體長數(shù)據(jù)，求解大連海域黃鮟鱇的體長生長方程和體重生長方程，按10 mm間隔統(tǒng)計(jì)出體長頻率數(shù)據(jù)。根據(jù)體長頻率數(shù)據(jù)，應(yīng)用ELEFAN I方法，求出Von Bertalanffy 方程的生長參數(shù)。生長曲率K 值為0.28a-1，最大漸進(jìn)體長L∞為576.4 mm，理論初始年齡t0為-0.44a，由體長與體重的關(guān)系式可以求得W∞為5 344.3 g。因此，大連海域黃鮟鱇的Von Bertalanffy生長方程為：體長生長方程為Lt=576.4×[1-e-0.28（t+0.44）]，體重生長方程為Wt=5 344.3×[1-e-0.28（t+0.44）]3.005 2。

2.3.2體長與體重的瞬時(shí)生長速度。

由體長生長方程Lt=576.4×[1-e-0.28（t+0.44）]，得出體長生長速度方程dLt/dt=161.4×e-0.28（t+0.44）；由體重生長方程Wt=5 344.3×[1-e-0.28（t+0.44）]3.005 2，得出體重生長速度方程dWt/dt=4 496.99× e-0.28（t+0.44）×[1-e-0.18（t+0.44）]2。根據(jù)體長和體重生長速度方程，分別繪制體長生長速度和體重生長速度對(duì)年齡的曲線圖。從圖4可以看出，體長生長速度隨年齡的增大而逐漸減小，體重生長速度隨年齡的增大先逐漸增大，達(dá)到393.23 g/a后開始逐漸減小。體重生長速度由快變慢的年齡即為體重生長拐點(diǎn)年齡，可由d2Wt/dt2=0求得，即t=4.2齡。將體重拐點(diǎn)年齡帶入生長方程，進(jìn)而求出體重生長拐點(diǎn)年齡時(shí)的體長和體重分別為419.2 mm和2 052.2 g。

2.4臨界年齡

臨界年齡是指一個(gè)世代的漁業(yè)資源群體在沒有捕撈的情況下生物量達(dá)到最大時(shí)的年齡，或重量的相對(duì)生長速度等于瞬時(shí)自然死亡率時(shí)的年齡[7]。根據(jù)Pauly[8]經(jīng)驗(yàn)公式：lnM=-0.006 6-0.279lnL∞+0.654 3lnK+0.463 4lnT，其中T為棲息海域年均水溫，該試驗(yàn)取20 ℃，求得M≈0.290 9 a-1；再根據(jù)Von Bertalanffy體重生長方程和瞬時(shí)自然死亡率M值，由dWt/（Wt×dt）≈ dB/（B×dt）=M，計(jì)算可得黃鮟鱇的臨界年齡為4.7齡，代入生長方程，得出相應(yīng)臨界體長為439.8 mm，相應(yīng)的臨界體重為2 369.5 g。

2.5黃鮟鱇的資源量的估算

2.5.1瞬時(shí)捕撈死亡率的估算。

從圖5可以看出，經(jīng)線性回歸，擬合的直線方程為ln（N/△t）=-0.991t+4.201，相關(guān)系數(shù)R為-0.822 4，方程的斜率為-0.991，因此瞬時(shí)總死亡率Z的估計(jì)量為0.991；已知M≈0.290 9a-1，則瞬時(shí)捕撈死亡率F=Z-M=0.700 1a-1

2.5.2黃鮟鱇的資源量估算。

根據(jù)黃鮟鱇的體長頻率數(shù)據(jù)，設(shè)黃鮟鱇最大體長的瞬時(shí)捕撈死亡率Ft為0.700 1；由近年來大連市黃鮟鱇產(chǎn)量的統(tǒng)計(jì)數(shù)據(jù)，設(shè)采樣量為黃鮟鱇年產(chǎn)量的1/10 000，運(yùn)用體長股分析法通過FISATⅡ軟件可估算大連海域黃鮟鱇的資源數(shù)量約為1 308.4萬尾，黃鮟鱇的資源量約為8 126.7 t。

3討論

魚類體長和體重是衡量魚類生長的重要指標(biāo)。大連海域黃鮟鱇體長為90～491 mm，平均體長為220.5 mm，體重為20～3 097 g，平均體重為540.3 g。大連海域黃鮟鱇生物學(xué)資料甚少，難以縱向比較研究其生長變化情況。與東海北部和黃海南部的黃鮟鱇群體相比，魚體規(guī)格整體偏小[4]。這可能與取樣時(shí)間、調(diào)查海域等不同有關(guān)。

黃鮟鱇資源近年來引起漁業(yè)管理部門和漁民的關(guān)注。科學(xué)有效管理黃鮟鱇資源尤為重要。合理的開捕規(guī)格和總允許漁獲量是實(shí)現(xiàn)其可持續(xù)利用的基礎(chǔ)。①開捕規(guī)格的制定。從生物學(xué)角度來看，不但要保證魚類個(gè)體的生長潛力，而且也要顧及魚類世代生物量的變化。從魚類個(gè)體生長潛能考慮，漁業(yè)利用開捕年齡應(yīng)在拐點(diǎn)年齡之后[9]；從魚類世代生物量變化的角度考慮，臨界年齡才是最佳開捕年齡，此時(shí)魚類世代的瞬間生物量最大，但大連海域黃鮟鱇可根據(jù)漁業(yè)管理部門的實(shí)際情況安排捕撈時(shí)間，因此為減少因自然死亡而造成的漁業(yè)損失，實(shí)際開捕年齡應(yīng)定在臨界年齡之前。根據(jù)計(jì)算所得的拐點(diǎn)年齡和臨界年齡，確定黃鮟鱇開捕年齡為4.2～4.7齡比較合適。為便于操作，實(shí)踐中可按相應(yīng)的體長（或體重）來控制，即黃鮟鱇開捕體長（或體重）控制在419.2～439.7 mm（或2 052.2～2 369.5 g）為宜。②捕撈產(chǎn)量的確定，需要根據(jù)現(xiàn)有資源量，綜合考慮海域情況后確定。計(jì)算的黃鮟鱇資源生物量為8 126.7 t。假設(shè)現(xiàn)在資源生物量為獲取最大持續(xù)產(chǎn)量MSY時(shí)生物量（即1/2B0），即根據(jù)FMSY=M，MSY=FMSYBMSY=0.5 MK=0.5 MB0，因此由MSY=0.5 MB0可得MSY=1 182 t，即年最大可持續(xù)捕撈量，可作為漁業(yè)管理部門漁業(yè)資源管理的參考依據(jù)。

參考文獻(xiàn)

[1] 蘇錦祥，李春生.中國動(dòng)物志：硬骨魚綱，鲀形目，海蛾魚目，喉盤魚目，鮟鱇目[M].北京：科學(xué)出版社，2002：350.

[2] 徐開達(dá)，賀舟挺，李鵬飛，等.東海北部、黃海南部黃鮟鱇的年齡和生長[J].漁業(yè)科學(xué)進(jìn)展，2010，31（6）：10-14.

[3] 張學(xué)健，沈偉，李建生，等.黃海南部和東海中北部黃鮟鱇數(shù)量分布及其與環(huán)境的關(guān)系[J].海洋漁業(yè)，2010，32（2）：59-65.

[4] 徐開達(dá)，賀舟挺，朱文斌，等.東海北部和黃海南部黃鮟鱇的數(shù)量分布及其群體結(jié)構(gòu)特征[J].大連海洋大學(xué)學(xué)報(bào)，2010，25（5）：465-470.

[5] 費(fèi)鴻年，張?jiān)娙?水產(chǎn)資源學(xué)[M].北京：中國科學(xué)技術(shù)出版社，1990.

[6] PAULY D.ELEFANI： Users instructional program listings[M].Manila：ICLARM，1980.

[7] PAULY D.Fish population dynamics in tropical waters： A manual for use with programmable calculators[J].ICLARM Stud Rev，1984（8）：325.

[8] PAULY D.On the interrelationships between natural mortality，growth parameters and mean environmental temperature in 175 fish stocks [J].J Cons Int Exp Mer，1980，39（2）：175-192.

[9] 尹增強(qiáng)，章守宇，汪振華，等.浙江嵊泗人工魚礁區(qū)小黃魚生長特征與資源合理利用的初步研究[J].中國生態(tài)農(nóng)業(yè)學(xué)報(bào)，2010（3）：588-594.