劉金立，陳新軍 ，李 綱

（上海海洋大學(xué)1.海洋科學(xué)學(xué)院；2.圖書館；3.國家遠(yuǎn)洋漁業(yè)工程技術(shù)研究中心；4.漁業(yè)可持續(xù)開發(fā)省部共建教育部重點(diǎn)實(shí)驗(yàn)室，上海 201306）

柔魚（Ommastrephesbartramii）為西北太平洋海域重要的經(jīng)濟(jì)頭足類之一［1］。我國于1993年對該資源進(jìn)行開發(fā)，之后捕撈規(guī)模和作業(yè)海域不斷擴(kuò)大，產(chǎn)量穩(wěn)定在6～10萬t間［2］，是我國遠(yuǎn)洋漁業(yè)重要的捕撈對象。國內(nèi)外學(xué)者在柔魚漁業(yè)生物學(xué)，包括種群結(jié)構(gòu)、洄游、年齡與生長等，漁場形成與變動、資源量變動以及漁業(yè)資源評估等方面進(jìn)行了深入的研究［2－9］，陳新軍等［10－11］利用 Delury衰減模型、基于貝葉斯統(tǒng)計(jì)方法的Schaefer模型對西北太平洋柔魚資源量進(jìn)行評估，并對其漁業(yè)管理策略做了風(fēng)險(xiǎn)分析；曹杰等［12］利用世代分析法估算了不同自然死亡系數(shù)下柔魚冬春生西部群體的初始資源量，以及以該群體為對象的漁業(yè)管理參考點(diǎn)。但在柔魚漁業(yè)資源生物經(jīng)濟(jì)學(xué)方面則沒有研究。目前，國內(nèi)利用生物經(jīng)濟(jì)學(xué)模型來研究漁業(yè)資源優(yōu)化配置及管理的文獻(xiàn)較少，主要集中在東黃海鮐魚方面［13－15］，且傳統(tǒng)的生物經(jīng)濟(jì)模型只考慮貼現(xiàn)率α＝0和α＝∞ 兩種特殊情況，未考慮貼現(xiàn)率變化對漁業(yè)資源優(yōu)化配置的影響［16］。漁業(yè)資源開發(fā)是1個系統(tǒng)工程，不僅涉及資源數(shù)量本身，而且還包括投入成本、經(jīng)濟(jì)效益以及社會就業(yè)、生態(tài)影響等。因此，本研究依據(jù)漁業(yè)資源經(jīng)濟(jì)學(xué)的理論和方法［16－17］，利用我國魷釣船在西北太平洋海域捕撈的柔魚產(chǎn)量、捕撈努力量和經(jīng)濟(jì)成本與漁獲數(shù)據(jù)，探討不同貼現(xiàn)率情況下柔魚漁業(yè)資源生物經(jīng)濟(jì)學(xué)的優(yōu)化配置，為西北太平洋柔魚漁業(yè)資源管理策略提供科學(xué)依據(jù)。

1 材料與方法

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

漁獲量數(shù)據(jù)來源于1996—2011年中國大陸魷釣船漁業(yè)生產(chǎn)統(tǒng)計(jì)數(shù)據(jù)，包括作業(yè)天數(shù)、日漁獲量、作業(yè)船次等，作業(yè)海域?yàn)?8°N～46°N，150°E～165°E，由中國遠(yuǎn)洋漁業(yè)協(xié)會上海海洋大學(xué)魷釣技術(shù)組提供。作業(yè)成本數(shù)據(jù)和柔魚價(jià)格數(shù)據(jù)以舟山市寧泰遠(yuǎn)洋漁業(yè)有限公司為基準(zhǔn)，近年來每天單船作業(yè)成本約為0.6萬元，近5年來柔魚平均價(jià)格約為1.0萬元/t。

1.2 分析方法

1.2.1 Gordon-Schaefer生物經(jīng)濟(jì)模型 一般的Gordon-Schaefer生物經(jīng)濟(jì)模型如下：

式中：π為利潤；TR為總收入；TC為總成本；p為價(jià)格；Y為產(chǎn)量；f為捕撈努力量；q為可捕系數(shù)；K為環(huán)境負(fù)載容量；r為種群內(nèi)稟增長率；c為單船作業(yè)成本［16］。根據(jù)Ichii等［18］和曹杰［19］的研究結(jié)果，柔魚種群內(nèi)稟增長率r為1.19；根據(jù)陳新軍等［11］研究結(jié)果，環(huán)境負(fù)載容量K為57.2萬t，可捕系數(shù)q為0.5×10－5。

由此模型可以推出最大可持續(xù)產(chǎn)量（Maximum Sustainable Yield，MSY）、最大經(jīng)濟(jì)產(chǎn)量（Maximum E-conomic Yield，MEY）和生物經(jīng)濟(jì)平衡點(diǎn)（Bio-economic Equilibrium，BE）及其所對應(yīng)的捕撈努力量fMSY、fMEY和fBE（見表1）。

表1 MSY、MEY和BE及對應(yīng)的捕撈努力量和產(chǎn)量Table 1 The arithmetic of fishing effort and yield corresponding to MSY，MEY and BE

1.2.2 最優(yōu)種群數(shù)量模型 以Schaefer剩余產(chǎn)量模型為基礎(chǔ)構(gòu)建最優(yōu)種群數(shù)量模型，可求得最優(yōu)種群資源量水平［20］。最優(yōu)種群數(shù)量（B＊）模型如下：

式中：B＊為最優(yōu)種群資源數(shù)量；Y＊為最優(yōu)持續(xù)產(chǎn)量；α為貼現(xiàn)率。

1.2.3 確定貼現(xiàn)率的近似方法 在漁業(yè)資源開發(fā)時，通常要權(quán)衡長遠(yuǎn)利益和當(dāng)前利益的關(guān)系，這種關(guān)系與貼現(xiàn)率有關(guān)［21］。葉昌臣和王文波［22］提出了確定貼現(xiàn)率的近似方法：

式中：D為當(dāng)前利益占的比重，D∈［0，1］；BD＊為根據(jù)D值變化的種群數(shù)量。

1.2.4 貼現(xiàn)率供給曲線對資源配置的影響 若漁業(yè)資源開發(fā)以MEY為追求目標(biāo)，則其貼現(xiàn)供給曲線就是最優(yōu)控制漁業(yè)的供給曲線［16］，假定需求是無限彈性，即漁獲價(jià)格為常數(shù)，以Schaefer模型為基礎(chǔ)，根據(jù)公式（3）和（4）可推出最優(yōu)持續(xù)產(chǎn)量Y＊是p的函數(shù)，繪制成曲線就是最優(yōu)控制漁業(yè)的供給曲線，由于Y＊也受貼現(xiàn)率α的影響，因此也稱為貼現(xiàn)供給曲線。結(jié)合西北太平洋柔魚漁業(yè)經(jīng)濟(jì)數(shù)據(jù)，可得到柔魚魚價(jià)與最優(yōu)種群數(shù)量B＊、最優(yōu)持續(xù)產(chǎn)量Y＊的關(guān)系。

2 結(jié)果

2.1 MSY、MEY和BE點(diǎn)的產(chǎn)量及其對應(yīng)捕撈努力量

根據(jù)Gordon-Schaefer生物模型估算，MSY約為17.02萬t，對應(yīng)的捕撈努力量約為11.9萬船次；MEY和BE所對應(yīng)的產(chǎn)量分別為16.29萬t和11.28萬t，所對應(yīng)的捕撈努力量分別為9.4萬船次和18.81萬船次。

2.2 不同貼現(xiàn)率情況下的柔魚資源量及產(chǎn)量估算

在貼現(xiàn)率α為0、0.05、0.1等11種情況下模擬了柔魚資源開發(fā)狀況（見圖1）。由圖1可知，當(dāng)貼現(xiàn)率為α＝0時，柔魚資源的最優(yōu)種群數(shù)量（B＊）為34.60萬t，對應(yīng)的最優(yōu)持續(xù)產(chǎn)量（Y＊）為16.29萬t，即為YMEY，捕撈努力量為9.40萬船次；當(dāng)α＝∞時，B＊＝BBE＝17.3萬t，對應(yīng)的產(chǎn)量YBE＝11.28萬t，其捕撈努力量為18.81萬船次；當(dāng)α＝0.43時，B＊＝BMSY＝28.60萬t，對應(yīng)的產(chǎn)量YMSY＝17.02萬t，其捕撈努力量為11.90萬船次。當(dāng)貼現(xiàn)率α為0.3～0.5時，其柔魚資源量相對穩(wěn)定，且達(dá)到最優(yōu)種群數(shù)量水平（見圖1）。

圖1 不同貼現(xiàn)率情況下柔魚資源的B＊值、Y＊值及其捕撈努力量Fig.1 The optimum populations（B＊），optimum sustainable yield（Y＊）and fishing effort corresponding to the different discounting rates

2.3 貼現(xiàn)率的確定

當(dāng)D＝0時，即不考慮當(dāng)前利益，以長遠(yuǎn)利益為主，其模擬的結(jié)果是柔魚漁業(yè)的利潤最高，產(chǎn)量中等，捕撈努力量最少，為9.40萬船次。當(dāng)D＝1，即不考慮長遠(yuǎn)利益，以當(dāng)前利益為主，其模擬結(jié)果捕撈努力量最高，為18.81萬船次，即該漁業(yè)就業(yè)機(jī)會最多，但產(chǎn)量最低，漁業(yè)利潤等于0。當(dāng)D＝0.279時，貼現(xiàn)率α＝0.46時，即當(dāng)前利益占長遠(yuǎn)利益的27.9%，此時，漁業(yè)產(chǎn)量達(dá)到了最大可持續(xù)產(chǎn)量，因此可選擇MSY為柔魚漁業(yè)的管理目標(biāo)（見圖2，3）。

圖2 不同D值情況下柔魚資源的B＊值、Y＊值及其捕撈努力量Fig.2 The optimum populations（B＊），optimum sustainable yield（Y＊）and fishing effort corresponding to the different Dvalues

圖3 不同D值下柔魚資源開發(fā)過程中的貼現(xiàn)率、作業(yè)成本及利潤Fig.3 The discount rate，cost and profit under the different Dvalues for Ommastrephes bartramii fishery

2.4 貼現(xiàn)率供給曲線對柔魚資源配置的影響

分別取貼現(xiàn)率α為0、0.3、0.5及∞繪制成柔魚漁業(yè)的貼現(xiàn)供給曲線（見圖4），α＝0時，為受控漁業(yè)；α＝∞時，為開放式漁業(yè)。由圖4可知，在一個開放式漁業(yè)中，柔魚漁獲價(jià)格小于2 097元/t時，其柔魚持續(xù)產(chǎn)量為0t；當(dāng)柔魚價(jià)格從2 097元/t逐漸增加時，其持續(xù)產(chǎn)量也隨之增加，當(dāng)漁獲價(jià)格增加到4 194元/t時（最大持續(xù)產(chǎn)量的價(jià)格，PMSY），其持續(xù)產(chǎn)量達(dá)到了最大可持續(xù)產(chǎn)量MSY。以后漁獲價(jià)物格繼續(xù)增加，其持續(xù)產(chǎn)量反而減少，這時產(chǎn)生了生物學(xué)捕撈過度（見圖4）。由圖4可知，在一個受控漁業(yè)中，其曲線不彎曲，而是漸近于MSY。當(dāng)α＝0.3時，漁獲產(chǎn)量達(dá)到MSY時的柔魚價(jià)格約為12 500元/t；當(dāng)α＝0.5時，漁獲產(chǎn)量達(dá)到 MSY時的柔魚價(jià)格約為9 195元/t。從圖4可知，曲線的彎曲程度與貼現(xiàn)率有關(guān)，貼現(xiàn)率越大曲線的彎曲程度也就越大。此外，貼現(xiàn)率對柔魚價(jià)格影響較大，貼現(xiàn)率越高，在產(chǎn)量達(dá)到MSY處的魚價(jià)越低（見表2）。

圖4 西北太平洋柔魚漁業(yè)貼現(xiàn)率供給曲線Fig.4 The discounted supply curve of the Ommastrephes bartramii fishery in the Northwestern Pacific Ocean

表2 西北太平洋柔魚魚價(jià)與B＊、Y＊的關(guān)系（α＝0.3和0.5）Table 2 The relationship among unit price and B＊，Y＊for Ommastrephes bartramii in the Northwestern Pacific Ocean（α＝0.3and 0.5）

3 討論

3.1 柔魚資源開發(fā)現(xiàn)狀分析

我國于1993年開始對西北太平洋柔魚資源進(jìn)行開發(fā)利用［2］，是目前捕撈柔魚最主要的國家和地區(qū)之一，1996—2008年產(chǎn)量基本上維持在8～12萬t，但在2009年，柔魚產(chǎn)量急劇下降，達(dá)到歷史最低水平。本文結(jié)合貼現(xiàn)率的生物經(jīng)濟(jì)模型，通過對比1996—2011年16年間的西北太平洋柔魚的產(chǎn)量數(shù)據(jù)，1999年漁獲最高，約為13.2萬t，捕撈努力量約為6.46萬船次，2009年漁獲產(chǎn)量為歷史最低水平約為3.68萬t，捕撈努力量約為2.63萬船次。本研究估算得到西北太平洋柔魚漁業(yè)的最大可持續(xù)產(chǎn)量MSY為17.02萬t，對應(yīng)的捕撈努力量為11.9萬船次。結(jié)合我國的生產(chǎn)實(shí)際，可知當(dāng)前西北太平洋柔魚資源狀況仍處在良好狀態(tài)，并未遭受過度捕撈，這一結(jié) 果與陳 新軍等［10－11］和曹 杰等［12］結(jié)果相似。目前我國在該海域投入捕撈努力量遠(yuǎn)低于漁獲量達(dá)到MSY時的捕撈努力量，因此尚可在該海域適當(dāng)投入較多的捕撈努力量。

3.2 貼現(xiàn)率和價(jià)格變動分析

貼現(xiàn)率對漁業(yè)資源開發(fā)將會產(chǎn)生很大的影響，但其影響程度也取決于不同的種類［17］，內(nèi)稟增長率大的種群，貼現(xiàn)率對其資源開發(fā)利用的影響不敏感，反之，則敏感性大［16］。因此，結(jié)合柔魚資源本身的生物經(jīng)濟(jì)學(xué)特性，合理確定貼現(xiàn)率可確保漁業(yè)資源的可持續(xù)利用。本研究通過對貼現(xiàn)率對西北太平洋柔魚資源量最優(yōu)種群數(shù)量和最優(yōu)持續(xù)產(chǎn)量的影響進(jìn)行分析，得知西北太平洋柔魚的最大可持續(xù)產(chǎn)量MSY為17.02萬t，最大經(jīng)濟(jì)產(chǎn)量 MEY為16.29萬t，當(dāng)貼現(xiàn)率α＝0.43時，最優(yōu)持續(xù)產(chǎn)量Y＊等于最大可持續(xù)產(chǎn)量MSY。這說明貼現(xiàn)率對西北太平洋柔魚資源種群數(shù)量水平的影響不明顯，但貼現(xiàn)率對柔魚價(jià)格影響較大。研究表明α在0.3～0.5間，同時D＝0.2～0.3，價(jià)格在9 195～12 500元/t間，柔魚資源種群資源量相對穩(wěn)定，其能達(dá)到最優(yōu)種群數(shù)量水平，可保證柔魚漁獲量始終處在MSY之下，從而實(shí)現(xiàn)柔魚資源的可持續(xù)開發(fā)和利用。根據(jù)調(diào)研的經(jīng)濟(jì)數(shù)據(jù)顯示，2008—2012年西北太平洋柔魚價(jià)格在9 250～12 500元/t間，其價(jià)格接近貼現(xiàn)率α為0.3～0.5之間的PMSY的價(jià)格。但本研究中未考慮海洋環(huán)境等因素對柔魚資源年間資源補(bǔ)充量的影響，已有研究表明［23－25］，厄爾尼諾、產(chǎn)卵場表溫等因素對柔魚資源補(bǔ)充量會產(chǎn)生一定的影響，從而導(dǎo)致柔魚漁業(yè)漁獲產(chǎn)量及其價(jià)格的劇烈波動。

3.3 漁業(yè)管理策略探討

由于我國在西北太平洋柔魚生產(chǎn)的漁獲產(chǎn)量約占該種類總產(chǎn)量的80%以上［2］，經(jīng)修正后在該海域的最高年產(chǎn)量為1999年的16.5萬t［11］，仍低于本研究得出在該海域的最大可持續(xù)產(chǎn)量MSY，可以采用MSY作為西北太平洋柔魚漁業(yè)的管理目標(biāo)，目前在該海域投入的年均捕撈努力量約為4萬船次，投入捕撈努力量最高的年份（1999年）也僅有6.46萬船次，遠(yuǎn)小于fMSY（11.9萬船次），因此，漁業(yè)管理部門可適當(dāng)?shù)卦黾硬稉婆α浚杂行ч_發(fā)柔魚資源。

同樣，柔魚漁業(yè)資源可持續(xù)開發(fā)也可以采用經(jīng)濟(jì)手段來調(diào)節(jié)。但由于柔魚價(jià)格受市場和資源量的影響，人們對其調(diào)節(jié)的能力較低。漁業(yè)管理及其有關(guān)部門可通過貼現(xiàn)率來進(jìn)行調(diào)節(jié)，例如，選擇適當(dāng)?shù)馁N現(xiàn)率（α＝0.3～0.5）可對柔魚的捕撈產(chǎn)量進(jìn)行調(diào)控，綜合考慮當(dāng)前利益與長遠(yuǎn)利益的關(guān)系（D＝0.2～0.3），確保柔魚資源保持在最優(yōu)種群數(shù)量水平，從而達(dá)到可持續(xù)利用西北太平洋柔魚資源的目的。

本研究應(yīng)用的基于貼現(xiàn)率的Gordon-Schaefer模型是單一種群的生物經(jīng)濟(jì)模型，沒有考慮種群間的競爭關(guān)系、捕食與被捕食關(guān)系，以及環(huán)境因子對西北太平洋柔魚資源的影響。后續(xù)的研究中應(yīng)系統(tǒng)開展基于生態(tài)和環(huán)境因素的西北太平洋柔魚生物經(jīng)濟(jì)學(xué)模型及其資源優(yōu)化配置的研究，并充分考慮r、K及q值的不確定性，但這是今后一個長期的研究過程。

［1］王堯耕，陳新軍.世界大洋性經(jīng)濟(jì)柔魚類資源及其漁業(yè)［M］.北京：海洋出版社，2005：152-153.

［2］Chen X J，Chen Y，Tian S Q，et al.An assessment of the west winter-spring cohort of neon flying squid （Ommastrephesbartramii）in the Northwest Pacific Ocean［J］.Fisheries Research，2008，92（2-3）：221-230.

［3］馬金，陳新軍，劉必林，等.北太平洋柔魚漁業(yè)生物學(xué)研究進(jìn)展［J］.上海海洋大學(xué)學(xué)報(bào)，2011，20（4）：563-570.

［4］李建華，陳新軍，劉必林，等.夏秋季西北太平洋柔魚漁業(yè)生物學(xué)的初步研究［J］.上海海洋大學(xué)學(xué)報(bào)，2011，20（6）：890-894.

［5］唐峰華，伍玉梅，樊偉.北太平洋柔魚生殖群體結(jié)構(gòu)特征及繁殖生物學(xué)［J］.中國海洋大學(xué)學(xué)報(bào)：自然科學(xué)版，2011，41（增2）：72-78.

［6］Yatsu A，Tanaka H，Mori J.Population structure of the neon flying squid，Ommastrephesbartramii，in the North Pacific Ocean［C］.//Okutani T.Contributed Papers to International Symposium on Large Pelagic Squids.Tokyo：Japan Marine Fishery Resources Research Center，Tokyo，1998：31-48.

［7］樊偉.衛(wèi)星遙感漁場漁情分析應(yīng)用研究－以西北太平洋柔魚漁業(yè)為例［D］.上海：華東師范大學(xué)，2004.

［8］陳新軍，田思泉.西北太平洋柔魚資源豐度時空分布的GAM模型分析［J］.集美大學(xué)學(xué)報(bào)：自然科學(xué)版，2006，11（4）：295-300.

［9］范江濤，陳新軍，曹杰，等.西北太平洋柔魚漁場變化與黑潮的關(guān)系［J］.上海海洋大學(xué)學(xué)報(bào)，2010，19（3）：378-384.

［10］陳新軍，錢衛(wèi)國，劉必林，等.利用衰減模型評估柔魚西部冬春生群體資源［J］.海洋湖沼通報(bào)，2008（2）：130-140.

［11］陳新軍，曹杰，劉必林，等.基于貝葉斯Schaefer模型的西北太平洋柔魚資源評估與管理［J］.水產(chǎn)學(xué)報(bào)，2011，35（10）：1572-1581

［12］曹杰，陳新軍，田思泉，等.基于世代分析法的西北太平洋柔魚冬春生西部群體資源評估［J］.中國海洋大學(xué)學(xué)報(bào)：自然科學(xué)版，2010，40（3）：37-42.

［13］張廣文，陳新軍，李綱.東黃海鮐魚生物經(jīng)濟(jì)模型及管理策略探討［J］.上海海洋大學(xué)學(xué)報(bào)，2009，4：447-452.

［14］王雅麗，陳新軍，李綱.基于貼現(xiàn)率的東黃海鮐魚動態(tài)生物經(jīng)濟(jì)模型分析［J］.資源科學(xué)，2011，11：2157-2161.

［15］王從軍，陳新軍，李綱.東、黃海鮐魚生物經(jīng)濟(jì)社會綜合模型的優(yōu)化配置研究［J］.上海海洋大學(xué)學(xué)報(bào)，2013（4）：623-628.

［16］陳新軍.漁業(yè)資源經(jīng)濟(jì)學(xué)［M］.北京：中國農(nóng)業(yè)出版社，2004：62-105.

［17］Clark，C W.Bioeconomic Modeling and Fisheries Management［M］.New York：Tohn Wiley，1985.

［18］Ichii T，Mahapatra K，Okamura H，et al.Stock assessment of the autumn cohort of neon flying squid（Ommastrephesbartramii）in the North Pacific based on past large-scale high seas driftnet fishery data［J］.Fisheries Research，2006，78（2）：286-297.

［19］曹杰.西北太平洋柔魚資源評估與管理［D］.上海：上海海洋大學(xué)，2010.

［20］Clark，C W.Mathematical Bioeconomics：The Optimal Manage-ment of Renewable Resources［M］.New York：John Wiley＆Sons，1990.

［21］馬中.環(huán)境與資源經(jīng)濟(jì)學(xué)［M］.北京：高等教育出版社，2004.

［22］葉昌臣，王文波.漁業(yè)管理中確定貼現(xiàn)率的一種近似方法［J］.水產(chǎn)科學(xué)，1990，9（1）：23-27.

［23］邵全琴，馬巍巍，陳卓奇，等.西北太平洋黑潮路徑變化與柔魚CPUE的關(guān)系研究［J］.海洋與湖沼，2005，36（2）：111-122.

［24］陳新軍，田思泉.利用GAM模型分析表溫和時空因子對西北太平洋海域柔魚資源狀況的影響［J］.海洋湖沼通報(bào)，2007（2）：104-113.

［25］田思泉，陳新軍，馮波，等.西北太平洋柔魚資源豐度與棲息環(huán)境的關(guān)系及其時空分布［J］.上海海洋大學(xué)學(xué)報(bào)，2009，18（5）：586-592.