西北太平洋柔魚漁業(yè)生物經(jīng)濟(jì)模型及其管理策略

劉金立，陳新軍 ，李 綱，李思亮

(1．上海海洋大學(xué)海洋科學(xué)學(xué)院，上海 201306;2．上海海洋大學(xué)圖書館，上海 201306;

3．上海海洋大學(xué)大洋漁業(yè)可持續(xù)開發(fā)省部共建教育部重點(diǎn)實(shí)驗(yàn)室，上海 201306;

4．國家遠(yuǎn)洋漁業(yè)工程技術(shù)研究中心，上海 201306;5．風(fēng)脈(武漢)可再生能源技術(shù)有限責(zé)任公司，武漢 430073)

柔魚(Ommastrephes bartramii)是重要的經(jīng)濟(jì)頭足類之一，廣泛分布在西北太平洋海域［1］。開發(fā)柔魚資源的主要國家和地區(qū)有中國大陸、日本及中國臺(tái)灣省。中國大陸于1993年開始開發(fā)該資源，之后捕撈規(guī)模和作業(yè)海域不斷擴(kuò)大，年產(chǎn)量穩(wěn)定在6—10萬t間［2］;1993年公海大型流刺網(wǎng)作業(yè)禁止后，日本以大型專業(yè)魷釣船為主進(jìn)行捕撈，年均產(chǎn)量為3—4萬t;中國臺(tái)灣省在2000年前的魷釣?zāi)昃a(chǎn)量約為2萬t，之后年產(chǎn)量均小于1萬t［1］，主要以秋刀魚漁業(yè)兼捕船隊(duì)為主。近幾年來西北太平洋柔魚的產(chǎn)量日趨下降，該資源可能已處于充分開發(fā)的狀態(tài)［3］。國內(nèi)外學(xué)者對(duì)西北太平洋柔魚進(jìn)行了較深入的研究，主要包括生物學(xué)特性、漁場(chǎng)形成及其開發(fā)狀況、資源量變動(dòng)以及漁業(yè)資源評(píng)估等方面［2-10］，但不同捕撈船隊(duì)、捕撈成本、漁獲價(jià)格等對(duì)柔魚資源開發(fā)的影響則沒有相應(yīng)的研究。本文利用多船隊(duì)的漁業(yè)生物經(jīng)濟(jì)模型［11-14］，考慮具有不同捕撈能力的中國大陸船隊(duì)、日本船隊(duì)及中國臺(tái)灣船隊(duì)同時(shí)捕撈西北太平洋柔魚時(shí)，在技術(shù)上的相互影響，并結(jié)合柔魚漁業(yè)的相關(guān)生物學(xué)和經(jīng)濟(jì)學(xué)參數(shù)，模擬不同管理方案下50a內(nèi)柔魚資源開發(fā)利用過程中，其資源量、捕撈努力量、漁獲產(chǎn)量及漁業(yè)利潤的動(dòng)態(tài)變化，為科學(xué)制定西北太平洋柔魚漁業(yè)資源管理策略提供科學(xué)依據(jù)。

1 材料與方法

1．1 材料來源

中國大陸船隊(duì)在西北太平洋的柔魚漁業(yè)生產(chǎn)統(tǒng)計(jì)數(shù)據(jù)由中國遠(yuǎn)洋漁業(yè)協(xié)會(huì)上海海洋大學(xué)魷釣技術(shù)組提供;作業(yè)成本和柔魚價(jià)格以舟山市寧泰遠(yuǎn)洋漁業(yè)有限公司等調(diào)查統(tǒng)計(jì)數(shù)據(jù)為基準(zhǔn)，近5年來每天單船作業(yè)成本約為0．6萬元，近5年來柔魚平均價(jià)格約為1．0萬元/t。中國臺(tái)灣省船隊(duì)的柔魚漁業(yè)統(tǒng)計(jì)數(shù)據(jù)來源于中國臺(tái)灣省對(duì)外漁業(yè)合作發(fā)展協(xié)會(huì)(http://www．ofdc．org．tw/index．htm)，每天單船作業(yè)成本為0．476萬元［15-16］;日本大型專業(yè)魷釣船船隊(duì)漁業(yè)統(tǒng)計(jì)數(shù)據(jù)來源于日本國際漁業(yè)資源現(xiàn)狀網(wǎng)站(http://kokushi．job．a(chǎn)ffrc．go．jp/index．html)，并假設(shè)單船作業(yè)成本為中國臺(tái)灣船隊(duì)的2倍。根據(jù)Ichii等［17］、曹杰［18］的研究結(jié)果，柔魚種群內(nèi)稟增長率 r設(shè)為1．19，根據(jù)陳新軍等［3］研究結(jié)果，環(huán)境負(fù)載容量K設(shè)為55.1萬t，1995年的西北太平洋柔魚的資源量為25萬t，初始資源量B0假定為1996年的資源量B1996，結(jié)合Schaefer資源量動(dòng)態(tài)模型，推算1996年資源量B1996為26．63萬t。根據(jù)中國臺(tái)灣魷釣船隊(duì)、日本魷釣船隊(duì)的年產(chǎn)量，結(jié)合中國魷釣船隊(duì)歷年的CPUE，來折算中國臺(tái)灣省和日本魷釣船隊(duì)的捕撈努力量，并考慮到中國臺(tái)灣船隊(duì)和日本船隊(duì)的作業(yè)特點(diǎn)，假設(shè)其可捕系數(shù)分別為中國大陸船隊(duì)的1/2和3/2(表1)。

表1 多船隊(duì)作業(yè)生物經(jīng)濟(jì)模型的相關(guān)參數(shù)設(shè)置Table 1 Parameters set for themultifleet bio-econom ic model

1．2 分析方法

1．2．1 多船隊(duì)作業(yè)的捕撈努力量動(dòng)態(tài)變化模型

在漁業(yè)資源開發(fā)過程中，各船隊(duì)投入到漁業(yè)生產(chǎn)中的捕撈努力量由其漁業(yè)的經(jīng)濟(jì)利益來決定的［11］，其捕撈努力量動(dòng)態(tài)變化模型如下:

式中，π(t)表示利潤函數(shù)，p表示漁獲價(jià)格，q表示可捕系數(shù)，c為單船作業(yè)成本，B(t)表示資源量變動(dòng)函數(shù)，f(t)表示捕撈努力量變動(dòng)函數(shù)，Y(t)表示漁獲產(chǎn)量函數(shù)，當(dāng)利潤為正值，作業(yè)船隊(duì)將會(huì)進(jìn)入該漁業(yè)，當(dāng)利潤為負(fù)值，作業(yè)船隊(duì)將會(huì)退出該漁業(yè);φ為作業(yè)船隊(duì)進(jìn)入或退出漁業(yè)的動(dòng)態(tài)參數(shù)，表示作業(yè)船隊(duì)長期的動(dòng)態(tài)變化，與捕撈努力量變動(dòng)所引起的成本變化密切相關(guān)，反映作業(yè)船隊(duì)利潤的變動(dòng)情況［19］。根據(jù) Anderson 等［12］和 Seijo 等［19］對(duì)海洋無脊椎動(dòng)物漁業(yè)在工業(yè)化捕撈過程中，作業(yè)船隊(duì)進(jìn)入或退出漁業(yè)的動(dòng)態(tài)參數(shù)設(shè)置，假定 φ 為0．5×10－5，本文不探討其影響程度。

1．2．2 多船隊(duì)作業(yè)的資源量動(dòng)態(tài)變化模型

漁業(yè)中，具有不同捕撈能力的3個(gè)作業(yè)船隊(duì)捕撈同一種群時(shí)，其漁業(yè)資源動(dòng)態(tài)模型表示如下:

式中，Y1(t)、Y2(t)、Y3(t)、f1(t)、f2(t)、f3(t)分別表示中國大陸船隊(duì)、中國臺(tái)灣船隊(duì)及日本船隊(duì)的產(chǎn)量和捕撈努力量;q1、q2、q3分別表示中國大陸船隊(duì)、中國臺(tái)灣船隊(duì)及日本船隊(duì)的可捕系數(shù)。由此模型可模擬出西北太平洋柔魚資源量及作業(yè)船隊(duì)的漁獲產(chǎn)量隨時(shí)間的動(dòng)態(tài)變化。

1．2．3 各船隊(duì)的漁業(yè)利潤隨捕撈努力量的動(dòng)態(tài)變化

利潤方程如下:

式中，π1(t)、π2(t)、π3(t)分別表示中國大陸船隊(duì)、中國臺(tái)灣船隊(duì)及日本船隊(duì)的利潤;c1、c2、c3分別表示中國大陸船隊(duì)、中國臺(tái)灣船隊(duì)及日本船隊(duì)的單船作業(yè)成本;p為柔魚價(jià)格。

1．2．4 模擬不同管理方案下作業(yè)船隊(duì)的動(dòng)態(tài)變化以基準(zhǔn)管理方案為標(biāo)準(zhǔn)，改變不同的生物和經(jīng)濟(jì)參數(shù)，模擬不同管理方案下各船隊(duì)50a內(nèi)的作業(yè)狀況及其漁業(yè)資源的動(dòng)態(tài)變化(表2)。

表2 不同管理方案下的參數(shù)設(shè)置Table 2 Parameters set corresponding to differentmanagement programs

2 結(jié)果

2．1 不同管理方案下各船隊(duì)捕撈努力量及資源量的動(dòng)態(tài)變化

根據(jù)捕撈努力量動(dòng)態(tài)變化模型及多船隊(duì)生物經(jīng)濟(jì)模型，模擬了不同管理方案下50a內(nèi)(1997—2047年)中國大陸船隊(duì)、中國臺(tái)灣省船隊(duì)及日本船隊(duì)的捕撈努力量及柔魚資源量隨時(shí)間的動(dòng)態(tài)變化情況(圖 1，圖2)。

對(duì)各船隊(duì)50a內(nèi)捕撈努力量的動(dòng)態(tài)變化進(jìn)行分析，中國大陸船隊(duì)投入的捕撈努力量占居絕對(duì)優(yōu)勢(shì)，其次是日本船隊(duì)，中國臺(tái)灣船隊(duì)最少。在管理方案1下，中國大陸船隊(duì)在1997—2003年間捕撈努力量逐年增加，達(dá)到最大值(6．16萬船次)后，開始逐年減少，當(dāng)減少至4．82萬船次后，又開始緩慢增加捕撈努力量，到2037年后捕撈努力量穩(wěn)定在約5．45萬船次;日本船隊(duì)在1997—2003年間捕撈努力量逐年緩慢增加，當(dāng)達(dá)到最大值(2．31萬船次)后，開始逐年減少，當(dāng)減少至1．52萬船次后略有回升，隨后捕撈努力量維持在約1．7—1．8萬船次間;中國臺(tái)灣船隊(duì)在1997—2002年間捕撈努力量略有增加，當(dāng)?shù)竭_(dá)0．5萬船次后，逐年緩慢減少，在2047年時(shí)捕撈努力量約為0．3萬船次(圖1)。其它9種方案的捕撈努力量和資源的動(dòng)態(tài)變化見圖1。

對(duì)不同管理方案下50a內(nèi)柔魚資源量的動(dòng)態(tài)變化進(jìn)行分析可知，在各方案下，柔魚漁業(yè)資源量的變化趨勢(shì)大致相同，即在1997—2007年內(nèi)資源量急劇下降，在2007—2017年間下降到資源量的最低水平，隨后資源量處于波動(dòng)和恢復(fù)階段，最后穩(wěn)定在一定的資源量水平(圖2)。其中，方案2在2007—2017年間資源量接近崩潰，之后隨著捕撈努力量的動(dòng)態(tài)調(diào)整資源量開始恢復(fù)，2042—2047年資源量維持在1．7—2．0萬 t間;方案 3、方案 5及方案 7在2037—2047年資源量維持在3萬t左右;方案1、方案6、方案8、方案9及方案10在2037—2047年資源量維持在4—5萬t間;方案4在2032—2047年資源量維持在較高的水平，約為6萬t(圖2)。

2．2 不同管理方案下各船隊(duì)漁獲產(chǎn)量及其利潤的動(dòng)態(tài)變化

從50a內(nèi)各船隊(duì)漁獲產(chǎn)量的動(dòng)態(tài)變化來看，各方案下，中國大陸船隊(duì)的漁獲產(chǎn)量最大，日本船隊(duì)次之，中國臺(tái)灣船隊(duì)最少(圖3)。不同方案下各船隊(duì)漁獲產(chǎn)量的變化趨勢(shì)大致相同，即在1997—2007年間各船隊(duì)漁獲產(chǎn)量呈急劇下降趨勢(shì)，在2007—2012年產(chǎn)量持續(xù)緩慢下降并達(dá)到最小值，之后隨著各船隊(duì)捕撈努力量的動(dòng)態(tài)調(diào)整漁獲產(chǎn)量出現(xiàn)波動(dòng)，最后趨于穩(wěn)定。中國大陸船隊(duì)在方案2下，在1997—2002年產(chǎn)量急劇下降，在2007—2022年產(chǎn)量下降到1萬t以下，之后產(chǎn)量逐漸增加，當(dāng)增加到7．4萬t后，又有逐年下降趨勢(shì)，至2047年時(shí)其產(chǎn)量約為2萬t;其它方案下，中國大陸船隊(duì)50a內(nèi)的產(chǎn)量變動(dòng)趨勢(shì)基本相似，其中，到2047年時(shí)，方案6約為4．2萬t，方案1、方案3、方案8及方案10的產(chǎn)量約為3萬t，方案4和方案9的產(chǎn)量約為2．5萬t，方案5和方案7的產(chǎn)量約為1．5萬t(圖3)。日本船隊(duì)在方案4下，到2047年時(shí)其產(chǎn)量約為3．2萬t，在方案1、方案5及方案7—方案10下，其產(chǎn)量在1．2—1．6萬 t間，方案2、方案3及方案5下，其產(chǎn)量在1萬t以下，接近于退出該漁業(yè)的水平。中國臺(tái)灣船隊(duì)，在方案9下，至2047年是產(chǎn)量為約為0．35萬t外，其余方案均為接近或直接退出該漁業(yè)的水平(圖3)。

圖1 不同管理方案下各作業(yè)船隊(duì)捕撈努力量隨時(shí)間的動(dòng)態(tài)變化Fig．1 The variation of fishing efforts for China mainland，Japan and Chinese Taiwan squid－jigging fleets over 50 years under different management program s

圖2 不同管理方案下柔魚資源量隨時(shí)間的動(dòng)態(tài)變化Fig．2 The biomass variation of Ommastrephes Bartramii over 50 years under differentmanagement programs

圖3 不同管理方案下各作業(yè)船隊(duì)柔魚產(chǎn)量隨時(shí)間的動(dòng)態(tài)變化Fig．3 The variation of catches for China mainland，Japan and Chinese Taiwan squid-jigging fleets over 50 years under different management program s

圖4 不同管理方案下各作業(yè)船隊(duì)捕撈柔魚的漁獲利潤隨時(shí)間的動(dòng)態(tài)變化Fig．4 The variation of fishing profits for China mainland，Japan and Chinese Taiwan squid-jigging fleets over 50 years under different management program s

從50a內(nèi)各船隊(duì)利潤的動(dòng)態(tài)變化來看，各方案下，1997—2002年各船隊(duì)均處于盈利狀態(tài)，且利潤逐年急劇遞減，其中，中國大陸船隊(duì)獲得的利潤最大，日本船隊(duì)次之，中國臺(tái)灣船隊(duì)最小(圖4)。在方案1下，中國大陸船隊(duì)和日本船隊(duì)在2003—2022年內(nèi)長期處于虧損狀態(tài)，在之后的25a內(nèi)長期在盈利與虧損之間波動(dòng)，而中國臺(tái)灣船隊(duì)則在2002年之后長期處于虧損狀態(tài)(圖4)。在方案2下，中國大陸船隊(duì)在2002—2026年間長期處于虧損狀態(tài)，且在2007—2012年內(nèi)嚴(yán)重虧損，隨后長期在盈利和虧損之間波動(dòng);日本船隊(duì)和中國臺(tái)灣船隊(duì)則長期在盈利和虧損之間波動(dòng)。在方案3下，中國大陸船隊(duì)從2005年之后長期在盈利和虧損之間波動(dòng);日本船隊(duì)在2003年之后長期在盈利和虧損之間波動(dòng)，而中國臺(tái)灣船隊(duì)則從2003年開始長期處于虧損狀態(tài)。在方案4下，中國大陸船隊(duì)和中國臺(tái)灣船隊(duì)在2002年之后長期處于虧損狀態(tài)，且中國大陸船隊(duì)在2007—2012年間嚴(yán)重虧損，而日本船隊(duì)則在2005—2017年間處于虧損狀態(tài)，在隨后的30年內(nèi)長期處于盈利狀態(tài)。在方案5下，中國大陸船隊(duì)和日本船隊(duì)在2002年之后長期處在盈利和虧損之間波動(dòng)，而中國臺(tái)灣船隊(duì)在2002年之后長期處在虧損狀態(tài)。在方案6下，中國大陸船隊(duì)在2003年之后長期處在盈利和虧損之間波動(dòng)，日本船隊(duì)和中國臺(tái)灣船隊(duì)則在2003年之后長期處于虧損狀態(tài)。在方案7下，中國大陸船隊(duì)和日本船隊(duì)在2003年之后長期處在盈利和虧損之間波動(dòng)，到2027年后為盈利階段，日本船隊(duì)獲得的利潤略高于中國大陸船隊(duì)，而中國臺(tái)灣船隊(duì)在2002年之后長期處于虧損狀態(tài)。在方案8下，中國大陸船隊(duì)和日本船隊(duì)在2003年之后長期處在盈利和虧損之間波動(dòng)，而中國臺(tái)灣船隊(duì)在2000年之后就長期處于虧損狀態(tài)。在方案9及方案10下，中國大陸船隊(duì)、中國臺(tái)灣船隊(duì)及日本船隊(duì)均在2003年之后長期處在盈利和虧損之間波動(dòng)(圖4)。

2．3 不同管理方案下柔魚漁業(yè)短期及長期累計(jì)效益分析

不同管理方案下累計(jì)前10a(1997—2007年)的產(chǎn)量及利潤，中國大陸船隊(duì)在方案6下的累計(jì)產(chǎn)量最高，約為 99．7 萬 t，方案 5 最低，約為 56．3 萬 t，方案6獲得的累計(jì)利潤最大，約為58．3億元，方案4最小，約為20．2億元;中國臺(tái)灣船隊(duì)在方案9下的累計(jì)產(chǎn)量最高，約為7．2萬 t，方案2最低，約為2．0萬t，在方案9下獲得的累計(jì)利潤最大，約為4．1億元，方案8最小，約為－1．4億元;日本船隊(duì)在方案4下的累計(jì)產(chǎn)量最高，約為43．1萬t，方案2最低，約為20．7萬t，在方案7下獲得的累計(jì)利潤最大，約為25．0億元，方案2最小，約為3．3億元。累計(jì)3個(gè)作業(yè)船隊(duì)前10a(1997—2007年)的總產(chǎn)量及總利潤，在方案6下累計(jì)總產(chǎn)量及利潤最高，分別為134．4萬t和70．8億元，方案2最低，分別為88．5萬 t和31．6億元(表3，表4)。對(duì)前30a(1997—2027年)和前50a(1997—2047年)的累計(jì)產(chǎn)量及利潤詳見表3和表4。

表3 不同管理方案下各船隊(duì)短期、中期及長期的累計(jì)產(chǎn)量/萬tTable 3 The cumulative production of each squid jigging fleet in short，medium and long-term under differentmanagement programs

表4 不同管理方案下各船隊(duì)短期、中期及長期的累計(jì)利潤/億元Table 4 The cumulative profits of each squid jigging fleet in short，medium and long－term under differentmanagement program s

3 討論與分析

3．1 柔魚資源開發(fā)現(xiàn)狀

中國大陸、日本及中國臺(tái)灣省是西北太平洋柔魚主要的生產(chǎn)國家和地區(qū)。1993年全面禁止流刺網(wǎng)作業(yè)后，主要作業(yè)方式以魷釣為主，年產(chǎn)量基本上在10—20萬 t間［1］。其中中國大陸船隊(duì) 2000—2008年間年產(chǎn)量基本上穩(wěn)定在8—10萬t間，但是2009年之后年柔魚產(chǎn)量急劇下降，約為5萬t;從捕撈努力量來看，1996年約為4．44萬船次，之后逐年增加，至1999年達(dá)到最大值，約為6．46萬船次，之后捕撈努力量呈減少趨勢(shì)，2007年之后明顯減少，且CPUE也有明顯的下降。中國臺(tái)灣船隊(duì)1995—1999年柔魚年產(chǎn)量約為2萬t，其中1998年達(dá)到最大值(3．48萬t)，之后產(chǎn)量急劇下降，至2006年后年產(chǎn)量在1萬t以下;從捕撈努力量來看，1996年為4478船次，之后逐年增加，至1998達(dá)到最大值，約為1．08萬船次，之后捕撈努力量大幅減少。日本為大型專業(yè)魷釣作業(yè)船隊(duì)，1993年以后其年產(chǎn)量不足1萬t，至1998年產(chǎn)量增加到5．36萬t，之后其柔魚產(chǎn)量明顯下降，年均產(chǎn)量不足2萬t;從捕撈努力量來看，1996年約為1．44萬船次，1998年達(dá)到最大值，約1．5萬船次，之后捕撈努力量呈下降趨勢(shì)。近年來，由于中國在西北太平洋的柔魚產(chǎn)量急劇下降，且CPUE出現(xiàn)下降趨勢(shì)，結(jié)合歷年中國大陸魷釣船隊(duì)在西北太平洋海域生產(chǎn)的漁獲產(chǎn)量及其變化情況［3］，可以初步判斷該資源可能已處于充分利用的狀態(tài)，但尚未遭受過度捕撈，與陳新軍等［20］和曹杰等［21］的研究結(jié)果相一致。因此，制定柔魚漁業(yè)管理措施時(shí)，需要綜合考慮生態(tài)效益、經(jīng)濟(jì)效益和社會(huì)效益等諸多方面的因素，以確保西北太平洋柔魚資源的可持續(xù)利用。

3．2 不同管理方案的模擬分析

通過對(duì)不同管理方案50a內(nèi)各船隊(duì)作業(yè)情況的動(dòng)態(tài)模擬，從捕撈努力量、產(chǎn)量及利潤等方面來看，在方案1(基準(zhǔn)方案)中，各船隊(duì)的作業(yè)情況幾乎與實(shí)際漁業(yè)相符，中國大陸船隊(duì)捕撈努力量占居絕對(duì)優(yōu)勢(shì)，約為日本船隊(duì)的3倍，日本船隊(duì)居中，而中國臺(tái)灣船隊(duì)最少，且長期呈遞減趨勢(shì)，直至最后退出漁業(yè);方案4中，各船隊(duì)能夠形成較好的良性競(jìng)爭(zhēng)，日本船隊(duì)的捕撈努力量最終接近并且可能超過中國大陸船隊(duì)，中國臺(tái)灣船隊(duì)也能保持較好的捕撈努力量的投入，且至2047年時(shí)該方案柔魚資源量最高，約為6萬t。方案9中，中國大陸船隊(duì)仍然占據(jù)主導(dǎo)地位，日本船隊(duì)次之，中國臺(tái)灣船隊(duì)最少，但捕撈努力量逐年呈遞增趨勢(shì)，且能獲得較好的產(chǎn)量及利潤。一般而言，對(duì)于多個(gè)不同國家或地區(qū)的船隊(duì)同時(shí)捕撈同一種漁業(yè)資源時(shí)，由于其捕撈能力及作業(yè)成本不同(捕撈能力越大其作業(yè)成本越高)，在漁業(yè)資源的最初利用階段，捕撈能力大的船隊(duì)占據(jù)優(yōu)勢(shì)，在公開入漁的情況下，資源量將逐步減少，并導(dǎo)致捕撈能力大的船隊(duì)退出漁業(yè)，最終被作業(yè)成本低的船隊(duì)所替代。因此，在柔魚資源優(yōu)化配置過程中，需綜合考慮日本、中國大陸及中國臺(tái)灣省等各作業(yè)船隊(duì)的實(shí)際情況，在確保各船隊(duì)能獲得相應(yīng)的生產(chǎn)效益的同時(shí)，也要確保柔魚資源免遭過度開發(fā)。

就長期效益而言，通過對(duì)50a的模擬結(jié)果可知，方案1(基準(zhǔn)方案)下的累計(jì)總產(chǎn)量及利潤均中等偏上，分別約為296．4萬t和33．1億元，其中，中國大陸船隊(duì)和日本船隊(duì)能夠獲得較好的效益，而中國臺(tái)灣船隊(duì)在中長期呈現(xiàn)虧損狀態(tài)，最終可能會(huì)退出漁業(yè);方案4的累計(jì)總產(chǎn)量最高，達(dá)到349．6萬t，但其總利潤較低，約為14．9億元，該方案下，日本船隊(duì)獲得了較好的利益，但犧牲了中國大陸船隊(duì)和中國臺(tái)灣船隊(duì)的中長期利益，中國大陸船隊(duì)和中國臺(tái)灣船隊(duì)在中長期均呈現(xiàn)虧損狀態(tài);方案6的累計(jì)總產(chǎn)量不是最高，約為329．6萬t，其總利潤卻達(dá)到了最大值，約為50．2億元，該方案下，中國大陸船隊(duì)獲得了的經(jīng)濟(jì)利益最高，但犧牲了日本船隊(duì)和中國臺(tái)灣船隊(duì)的長期利益;方案9的累計(jì)總產(chǎn)量和總利潤均為中等，分別約為289萬t和32．4億元，方案10的累計(jì)總產(chǎn)量和總利潤均為中等偏上，分別為295．2萬t和35．4億元，這兩種方案下均兼顧了各個(gè)船隊(duì)的利益，但方案10下中國臺(tái)灣船隊(duì)獲得的中長期效益較低。

3．3 柔魚漁業(yè)資源管理措施的探討

通過模擬10種方案下各船隊(duì)50a內(nèi)柔魚漁業(yè)作業(yè)狀況及其漁業(yè)資源的動(dòng)態(tài)變化，根據(jù)模擬結(jié)果，綜合考慮不同國家及地區(qū)柔魚漁業(yè)的生態(tài)效益、經(jīng)濟(jì)效益及社會(huì)效益等管理目標(biāo)，并結(jié)合模擬的短期和長期效益等，來科學(xué)制定柔魚漁業(yè)的管理措施。

在漁業(yè)生產(chǎn)的實(shí)際過程中，各國船隊(duì)或漁業(yè)企業(yè)均以追求漁業(yè)利潤的最大化為目標(biāo)。在柔魚漁業(yè)管理過程中，也不能以犧牲某一船隊(duì)的利益來換取另一船隊(duì)的利益，需要協(xié)調(diào)不同國家或地區(qū)的作業(yè)船隊(duì)之間相互利益，使之達(dá)到一定的平衡，因此，建議以方案9(提高中國臺(tái)灣船隊(duì)的捕撈系數(shù))作為柔魚漁業(yè)的參考管理目標(biāo)，可以兼顧各作業(yè)船隊(duì)的經(jīng)濟(jì)利益和社會(huì)效益。

從柔魚資源優(yōu)化配置的角度來看，既要使各國船隊(duì)獲得較好的漁業(yè)利潤，又能使柔魚資源免遭過度捕撈，保持較好的資源量水平，在模擬50a內(nèi)資源量的動(dòng)態(tài)變化過程中，至25a后方案4的資源量水平保持最好，約保持在6萬t以上，因此，建議以方案4作為參考管理目標(biāo)，并在此基礎(chǔ)上，對(duì)方案4進(jìn)行適當(dāng)優(yōu)化，即提高中國大陸船隊(duì)作業(yè)成本，同時(shí)適當(dāng)降低日本船隊(duì)的可捕系數(shù)，可以確保柔魚資源量保持在更高的水平。同時(shí)，建議組建西北太平洋柔魚漁業(yè)國際管理組織，協(xié)調(diào)管理各國船隊(duì)的漁業(yè)活動(dòng)，共同開發(fā)和管理柔魚資源，確保該資源的可持續(xù)利用。

本文應(yīng)用的多船隊(duì)動(dòng)態(tài)生物經(jīng)濟(jì)模型主要是針對(duì)單一種群的漁業(yè)資源開發(fā)的動(dòng)態(tài)模擬，由于柔魚是短生命周期種類，其資源量極易受到海洋環(huán)境波動(dòng)的影響［4，22-25］，本研究中沒有考慮環(huán)境變化對(duì)柔魚資源開發(fā)的影響。漁業(yè)資源開發(fā)是一個(gè)綜合系統(tǒng)，在漁業(yè)資源生物經(jīng)濟(jì)模型的研究過程中需要考慮的因子也越來越多，多魚種、多船隊(duì)、海洋環(huán)境因子以及社會(huì)效益、價(jià)格和成本變動(dòng)等因素都被考慮到了漁業(yè)資源優(yōu)化配置的模型中，從而逐漸發(fā)展成為一個(gè)復(fù)雜的生物經(jīng)濟(jì)社會(huì)綜合模型。因此，需要充分掌握柔魚資源生物學(xué)、生態(tài)學(xué)、經(jīng)濟(jì)學(xué)及社會(huì)統(tǒng)計(jì)等數(shù)據(jù)，利用計(jì)算機(jī)模擬開展西北太平洋柔魚資源生物經(jīng)濟(jì)學(xué)模型及其資源優(yōu)化配置的研究。

［1］ Wang Y G，Chen X J．The Oceanic Economic Squid Resources and Fisheries in the World． Beijing:Ocean Press， 2005:124-159．

［2］ Chen X J，Chen Y，Tian S Q，Liu B L，Qian W G．An assessment of the west winter-spring cohort of neon flying squid(Ommastrephes bartramii) in the Northwest Pacific Ocean．Fisheries Research，2008，92(2/3):221-230．

［3］ Chen X J，Cao J，Liu B L，Lu H J，Tian S Q，Ma J．Stock assessment and management of Ommartrephesbartramii by using a Bayesian Schaefer model in the Northwestern Pacific Ocean．Journal of Fisheries of China，2011，35(10):1572-1581．

［4］ Yatsu A，Mori J．Early growth of the autumn cohort of neon flying squid，Ommastrephes bartramii，in the North Pacific Ocean．Fisheries Research，2000，45(2):189-194．

［5］ Chen X J，Xu L X，Tian SQ．Spatial and temporal analysis of Ommastrephe bartrami resources and its fishing ground in North Pacific Ocean．Journal of Fisheries of China，2003，27(4):334-342．

［6］ WangW Y，Shao QQ，Xue Y C，Zhang TY．On the relationship between the resources of Ommastrephes Bartrami and marine environment in the Northwest Pacific Ocean based on GIS．Geo-Information Science，2003，5(1):39-44．

［7］ Ichii T，Mahapatra K，Sakai M，Inagake D，Okada Y．Differing body size between the autumn and the winter-spring cohorts of neon flying squid (Ommastrephes bartramii) related to the oceanographic regime in the North Pacific: A hypothesis．Fisheries Oceanography，2004，13(5):295-309．

［8］ Tang F H，Wu Y M，F(xiàn)an W．Population structure and reproductive biology of mature Neon Flying Squid in the North Pacific Ocean．Periodical of Ocean University of China，2011，41(7/8):72-78．

［9］ Ma J，Chen X J，Liu B L，Tian SQ，Li S L，Cao J．Review of fisheries biology of neon flying squid(Ommastrephesbartramii)in the North Pacific Ocean．Journal of Shanghai Ocean University，2011，20(4):563-570．

［10］ Fan W．A Study on Application of Satellite Remote Sensing in Marine Fishing-Ground Analysis and Fishing Condition Forecasting-A Case of Ommastrephes bartrami Fisheries in Northwest Pacific Ocean［D］．Shanghai:East China Normal University，2004．

［11］ Chen X J． Fisheries Resource Economics． Beijing:China Agriculture Press，2004:62-105．

［12］ Anderson LG，Seijo JC．Bioeconomics of FisheriesManagement．Ames，Iowa:A John Wiley ＆ Sons，Ltd．，2010:189-229．

［13］ Seijo J C，Defeo O，Salas S．Fisheries Bioeconomics Theory，Modelling and Management．Rome: Food and Agriculture Organization of the United Nations Rome，1998．

［14］ Zhang G W，Chen X J，Li S L，Li G．A multi-fleet and bioeconomic simulationmodel and its application to Chub Mackerel in the East China Sea and the Yellow Sea．Resources Science，2010，32(8):1627-1633．

［15］ Yang L C．Economic Analysis on the Production of Deep-sea Squid Fishery in Taiwan［D］．Keelung:National Taiwan Ocean University，2005．

［16］ Lin P C．A Study on the Adjustmentof Deep-sea Squid Fishery in Taiwan［D］．Keelung:National Taiwan Ocean University，2007．

［17］ Ichii T，Mahapatra K，Okamura H，Okadac Y．Stock assessment of the autumn cohort of neon flying squid (Ommastrephes bartramii)in the North Pacific based on past large-scale high seas driftnet fishery data．Fisheries Research，2006，78(2/3):286-297．

［18］ Cao J．Stock Assessment and Risk Analysis of Management Strategies for neon flying squid(Ommastrephes bartramii)in the Northwest Pacific ocean［D］．Shanghai: Shanghai Ocean University，2010．

［19］ Seijo J C，Defeo O．Dynamics of resource and fishermen behaviour in coastal invertebrate fisheries//Antona M，Catanzano J，Sutinen J，eds．Proceedings of the Sixth Conference of the International Institute of Fisheries Economics and Trade．Paris，F(xiàn)rancia，1994，(1):209-222．

［20］ Chen X J，Qian W G，Liu B L，Tian S Q．An assessment of western winter-spring stock of neon flying squid，Ommastrephes bartramii，in the Northwest Pacific Ocean using the depletion model．Transactions of Oceanology and Limnology，2008，(2):130-140．

［21］ Cao J，Chen X J，Tian S Q，Liu B L．Assessment of western winter-spring stock of Neon Flying Squid (Ommastrephes bartramii) in the Northwest Pacific Ocean based on cohort analysis．Periodical of Ocean University of China，2010，40(3):37-42．

［22］ Rodhouse P G．Managing and forecasting squid fisheries in variable environments．Fisheries Research，2001，54(1):3-8．

［23］ Shao Q Q，Ma W W，Chen Z Q，You Z M，Wang W Y．Relationship between kuroshiomeander pattern and Ommastrephes bartrami CPUE in Northwest Pacific Ocean．Oceanologia et Limnologia Sinica，2005，36(2):111-122．

［24］ Chen X J，Liu B L，Tian SQ，Qian W G，Li G．Forecasting the fishing ground of Ommastrephes bartramii with SST-based habitat suitability modelling in Northwestern Pacific．Oceanologia et Limnologia Sinica，2009，40(6):707-713．

［25］ Cao J，Chen X J，Liu B L，Tian SQ，Tian W G．Review on the relationship between stock recruitment of squid and oceanographic environment．Journal of Shanghai Ocean University，2010，19(2):232-239．

參考文獻(xiàn):

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

［3］ 陳新軍，曹杰，劉必林，陸化杰，田思泉，馬金．基于貝葉斯Schaefer模型的西北太平洋柔魚資源評(píng)估與管理．水產(chǎn)學(xué)報(bào)，2011，35(10):1572-1581．

［5］ 陳新軍，許柳雄，田思泉．北太平洋柔魚資源與漁場(chǎng)的時(shí)空分析．水產(chǎn)學(xué)報(bào)，2003，27(4):334-342．

［6］ 王文宇，邵全琴，薛允傳，仉天宇．西北太平洋柔魚資源與海洋環(huán)境的GIS空間分析．地球信息科學(xué)，2003，5(1):39-44．

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

［9］ 馬金，陳新軍，劉必林，田思泉，李思亮，曹杰．北太平洋柔魚漁業(yè)生物學(xué)研究進(jìn)展．上海海洋大學(xué)學(xué)報(bào)，2011，20(4):563-570．

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

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

［14］ 張廣文，陳新軍，李思亮，李綱．基于多船隊(duì)作業(yè)的東、黃海鮐魚生物經(jīng)濟(jì)模型及管理策略．資源科學(xué)，2010，32(8):1627-1633．

［15］ 楊立中．臺(tái)灣遠(yuǎn)洋魷釣漁業(yè)生產(chǎn)之經(jīng)濟(jì)分析［D］．基隆:臺(tái)灣國立海洋大學(xué)，2005．

［16］ 林報(bào)儉．臺(tái)灣遠(yuǎn)洋魷釣漁業(yè)產(chǎn)業(yè)調(diào)適之研究［D］．基隆:臺(tái)灣國立海洋大學(xué)，2007．

［18］ 曹杰．西北太平洋柔魚資源評(píng)估與管理［D］．上海:上海海洋大學(xué)，2010．(未鏈接到本條文獻(xiàn)信息)

［20］ 陳新軍，錢衛(wèi)國，劉必林，田思泉．利用衰減模型評(píng)估柔魚西部冬春生群體資源．海洋湖沼通報(bào)，2008，(2):130-140．

［21］ 曹杰，陳新軍，田思泉，劉必林．基于世代分析法的西北太平洋柔魚冬春生西部群體資源評(píng)估．中國海洋大學(xué)學(xué)報(bào):自然科學(xué)版，2010，40(3):37-42．

［23］ 邵全琴，馬巍巍，陳卓奇，游智敏，王文宇．西北太平洋黑潮路徑變化與柔魚CPUE的關(guān)系研究．海洋與湖沼，2005，36(2):111-122．

［24］ 陳新軍，劉必林，田思泉，錢衛(wèi)國，李綱．利用基于表溫因子的棲息地模型預(yù)測(cè)西北太平洋柔魚(Ommastrephes bartramii)漁場(chǎng)．海洋與湖沼，2009，40(6):707-713．

［25］ 曹杰，陳新軍，劉必林，田思泉，錢衛(wèi)國．魷魚類資源量變化與海洋環(huán)境關(guān)系的研究進(jìn)展．上海海洋大學(xué)學(xué)報(bào)，2010，19(2):232-239．