山東省海洋牧場經(jīng)濟(jì)效益的實(shí)證分析

張婧一，宋 欣*，張海英，王俊鵬

（1.山東省淡水漁業(yè)研究院，山東濟(jì)南 250013；2.平度市自然資源局，山東青島 266700）

山東省是海洋大省，毗鄰海域面積約159 500 km2，海洋資源豐度指數(shù)居全國首位［1］。2005年，《山東省漁業(yè)資源修復(fù)行動(dòng)規(guī)劃》開始實(shí)施，山東省沿海城市大范圍開展海洋牧場建設(shè)［2］。2014年，山東省提出“投礁型、田園型、底播型、裝備型、游釣型”五大類海洋牧場協(xié)同發(fā)展的新理念［3］。2017年，《山東省海洋牧場建設(shè)規(guī)劃（2017—2020年）》正式印發(fā)。2018年，《山東省現(xiàn)代化海洋牧場建設(shè)綜合試點(diǎn)方案》獲國家批復(fù)［4－5］。海洋牧場不僅可以修復(fù)海洋生態(tài)環(huán)境，發(fā)揮其生態(tài)效益，解決漁民就業(yè)問題，發(fā)揮其社會(huì)效益，還可以更加高效地發(fā)展海洋漁業(yè)產(chǎn)業(yè)，發(fā)揮其經(jīng)濟(jì)效益［6］。目前，山東省海洋牧場建設(shè)盡管取得了初步成效，但是在發(fā)展中仍存在用海面積過大等諸多問題［7－8］，因此，山東省在發(fā)展海洋牧場的過程中需要不斷優(yōu)化其發(fā)展模式，促進(jìn)海洋牧場由高速發(fā)展轉(zhuǎn)向高質(zhì)量發(fā)展，實(shí)現(xiàn)海域資源的集約使用和可持續(xù)發(fā)展。本文選取40處山東省海洋牧場，利用2017—2019年海洋牧場收入、利潤、海域面積和海洋牧場類型等數(shù)據(jù)，對(duì)不同類型海洋牧場經(jīng)營效益的差異和同類型海洋牧場經(jīng)營效益的差異進(jìn)行分析，以期為山東省海洋牧場發(fā)展提供一定參考依據(jù)。

1 山東省海洋牧場總體情況

截至2020年7月，山東省已建成國家級(jí)海洋牧場44處，省級(jí)海洋牧場105處，居全國首位。這些海洋牧場分布在山東省煙臺(tái)、威海、青島、日照、東營、濰坊和濱州7個(gè)沿海城市，其中省級(jí)海洋牧場數(shù)量在煙臺(tái)和威海最多，青島和日照次之，東營、濰坊和濱州最少。具體分布情況如表1所示。

表1 山東省海洋牧場分布表Tab.1 Distribution of marine ranches in Shandong Province （處）

本文以原山東省海洋與漁業(yè)廳《關(guān)于公布2018年省級(jí)海洋牧場建設(shè)示范項(xiàng)目名單的通知（魯海漁函〔2018〕285號(hào)）》所公示的72處海洋牧場作為研究總體［9］，海洋牧場海域面積、類型來源于此文件。2017—2019年海洋牧場收入、利潤數(shù)據(jù)取自于2020年度山東省海洋軟科學(xué)研究課題“海洋牧場統(tǒng)計(jì)體系建設(shè)研究”的調(diào)研數(shù)據(jù)。本研究所涉及的海洋牧場收入主要包括海洋牧場產(chǎn)出的魚、蝦、蟹、參和貝等海珍品的收入，人工魚礁釣捕的收入和休閑海釣的收入；利潤主要是指收入扣除維持海洋牧場正常運(yùn)營所需的費(fèi)用。

按照海洋牧場功能和建設(shè)方式的不同，海洋牧場分為游釣型海洋牧場、投礁型海洋牧場、底播型海洋牧場、裝備型海洋牧場和田園型海洋牧場5種類型［10］。其中，從功能定位來講，底播型海洋牧場側(cè)重于底棲貝類和底棲海珍品的增殖；游釣型海洋牧場側(cè)重于休閑海釣；田園型海洋牧場側(cè)重于立體、循環(huán)和生態(tài)養(yǎng)殖；投礁型海洋牧場側(cè)重于投放人工魚礁，改善海洋生態(tài)環(huán)境；裝備型海洋牧場側(cè)重于在深遠(yuǎn)海建設(shè)大型智能網(wǎng)箱、深海養(yǎng)殖工船和養(yǎng)殖平臺(tái)?；跀?shù)據(jù)的可獲取性和完整性，選取了40處山東省海洋牧場作為研究對(duì)象。根據(jù)相關(guān)文件及調(diào)研情況，將40處海洋牧場進(jìn)行了類型分類，具體分布情況如表2所示。從表中可以看出，底播型海洋牧場主要分布在濰坊、東營和濱州3市，其他類型主要分布在煙臺(tái)、威海和日照3市。

表2 40處山東省海洋牧場分布情況表Tab.2 Distribution of 40 marine ranches in Shandong Province （處）

2 實(shí)證分析

2.1 收入、利潤的描述性統(tǒng)計(jì)

近年來，山東在發(fā)展思路、發(fā)展模式和發(fā)展方法上大膽嘗試，積極作為，海洋牧場發(fā)展迅速。2017—2019年，40處海洋牧場收入平均值、利潤平均值逐年增加（表3）。2018年收入平均值環(huán)比增長率為7.23%，2019年收入平均值環(huán)比增長率為6.18%，收入平均值增速放緩。2018年利潤平均值環(huán)比增長率為42.83%，2019年利潤平均值環(huán)比增長率為17.57%，利潤平均值增速放緩。同時(shí)期利潤平均值的環(huán)比增長率大于收入平均值的環(huán)比增長率，說明海洋牧場平均收入利潤率逐年提高［11］。2017—2019年收入平均值的標(biāo)準(zhǔn)差逐年降低，說明收入平均值反映數(shù)據(jù)集中趨勢的代表性顯著；利潤平均值的標(biāo)準(zhǔn)差逐年升高，說明利潤平均值反映數(shù)據(jù)集中趨勢的代表性不顯著。

表3 收入、利潤的描述性分析Tab.3 Descriptive analysis of income and profit

偏度是對(duì)分布偏斜方向和程度的判斷與測定，同時(shí)在一定程度上也反映出均值作為集中趨勢指標(biāo)的代表性和數(shù)據(jù)分布的離散程度［12］。但要準(zhǔn)確地測定分布的偏斜程度并進(jìn)行比較分析，通常需要計(jì)算分布的偏態(tài)系數(shù)，本文采用的偏態(tài)系數(shù)的計(jì)算方法如下：

式（1）中，α3為偏態(tài)系數(shù)，σ3為標(biāo)準(zhǔn)差的三次方，Xi為第i個(gè)樣本數(shù)據(jù)，k表示樣本容量，為樣本平均值，F(xiàn)i為第i個(gè)樣本數(shù)據(jù)出現(xiàn)的次數(shù)，N為樣本個(gè)數(shù)。

當(dāng)α3大于0，表示數(shù)據(jù)分布呈現(xiàn)右偏態(tài)，位于均值左邊的數(shù)據(jù)比位于右邊的多，左邊的離散度比右邊弱；當(dāng)α3小于0，表示數(shù)據(jù)分布呈現(xiàn)左偏態(tài)，位于均值左邊的數(shù)據(jù)比位于右邊的少，左邊的離散度比右邊強(qiáng)；當(dāng)α3等于0，數(shù)據(jù)分布符合正態(tài)分布。2017—2019年海洋牧場收入平均值、利潤平均值均呈現(xiàn)右偏態(tài)，說明左邊的離散度比右邊弱。

峰度用于衡量數(shù)據(jù)分布相對(duì)于正態(tài)分布而言是陡峭還是平緩。測定峰度狀態(tài)的統(tǒng)計(jì)指標(biāo)為峰度系數(shù)，它是離差四次方的平均數(shù)再除以標(biāo)準(zhǔn)差的四次方，即：

式（2）中，α4為峰度系數(shù)，σ4為標(biāo)準(zhǔn)差的四次方，Xi為第i個(gè)樣本數(shù)據(jù)，k表示樣本容量，為樣本平均值，F(xiàn)i為第i個(gè)樣本數(shù)據(jù)出現(xiàn)的次數(shù)，N為樣本個(gè)數(shù)。

α4等于3時(shí)，數(shù)據(jù)分布符合正態(tài)分布；α4大于3時(shí)，峰的形狀比較尖，比正態(tài)分布峰要陡峭；α4小于3時(shí)，峰的形狀比較緩，比正態(tài)分布峰要平緩。2017—2019年40處海洋牧場收入平均值、利潤平均值α4均大于3，說明40處海洋牧場收入平均值、利潤平均值的數(shù)據(jù)分布相對(duì)正態(tài)分布而言比較陡峭。

2.2 收入、利潤和海域面積偏相關(guān)分析

很多情況下，兩個(gè)變量的數(shù)值會(huì)同時(shí)受其他變量的影響，這時(shí)候就需要把其他變量控制住，計(jì)算出控制其他變量影響后的相關(guān)系數(shù)（稱為偏相關(guān)系數(shù)），這就是偏相關(guān)分析［13］。偏相關(guān)系數(shù)是真正反映兩個(gè)變量相關(guān)關(guān)系的統(tǒng)計(jì)量。P值是針對(duì)原假設(shè)H0：假設(shè)兩變量無線性相關(guān)而言的。本文將顯著性水平設(shè)置為0.01，所以只需要將P值和0.01進(jìn)行比較：如果P值小于0.01，就拒絕原假設(shè)H0，說明兩變量有線性相關(guān)的關(guān)系；如果大于0.01，則一般認(rèn)為無線性相關(guān)關(guān)系。本文利用SPSS24.0軟件中的偏相關(guān)分析功能，分析了2017—2019年山東省40處海洋牧場收入、利潤受海域面積的影響程度（表4）。

表4 2017—2019年收入、利潤和海域面積偏相關(guān)分析結(jié)果Tab.4 Partial correlation analysis results of income，profit and sea area during 2017—2019

2017年不控制海域面積時(shí)，收入與利潤的相關(guān)系數(shù)是0.558，相關(guān)系數(shù)檢驗(yàn)的概率P值是0.000，小于顯著性水平0.01；控制海域面積時(shí)，收入與利潤的偏相關(guān)系數(shù)是0.559，偏相關(guān)系數(shù)檢驗(yàn)的概率P值是0.000，小于顯著性水平0.01。偏相關(guān)系數(shù)0.559大于相關(guān)系數(shù)0.558，所以剔除海域面積的影響后，2017年收入與2017年利潤的相關(guān)關(guān)系稍弱。

2018年不控制海域面積時(shí)，收入與利潤的相關(guān)系數(shù)是0.496，相關(guān)系數(shù)檢驗(yàn)的概率P值是0.001，小于顯著性水平0.01；控制海域面積時(shí)，收入與利潤的偏相關(guān)系數(shù)是0.496，偏相關(guān)系數(shù)檢驗(yàn)的概率P值是0.001，小于顯著性水平0.01。偏相關(guān)系數(shù)等于相關(guān)系數(shù)，可見，海域面積對(duì)2018年收入與2018年利潤的相關(guān)關(guān)系沒有影響。

2019年不控制海域面積時(shí)，收入與利潤的相關(guān)系數(shù)是0.583，相關(guān)系數(shù)檢驗(yàn)的概率P值是0.000，小于顯著性水平0.01；控制海域面積時(shí)，收入與利潤的偏相關(guān)系數(shù)是0.582，偏相關(guān)系數(shù)檢驗(yàn)的概率P值是0.000，小于顯著性水平0.01。相關(guān)系數(shù)0.583大于偏相關(guān)系數(shù)0.582，所以剔除海域面積的影響后，2019年收入與2019年利潤相關(guān)性更強(qiáng)［14］。

從上面分析可以看出，2017年收入與利潤的相關(guān)性受海域面積影響更大，2018年收入與利潤的相關(guān)性不受海域面積影響，2019年收入與利潤的相關(guān)性受海域面積影響減弱?？梢?，山東省海洋牧場經(jīng)濟(jì)效益可能會(huì)受到海洋牧場海域面積的影響，但海域面積并不是影響海洋牧場經(jīng)濟(jì)效益的必然要素。

2.3 收入、利潤和類型分析

底播型、田園型海洋牧場多建設(shè)于淺海地區(qū)，用于淺海養(yǎng)殖，基本設(shè)施配備相對(duì)簡單；投礁型、裝備型和游釣型海洋牧場需投入人工魚礁，前期投入較大；裝備型海洋牧場建設(shè)于深海地區(qū)，投入較大。如表5所示，游釣型海洋牧場平均收入、平均利潤和平均收入利潤率逐年向好發(fā)展；投礁型海洋牧場2018年的平均利潤和平均收入利潤率比2017年和2019年高；底播型海洋牧場平均利潤和平均收入利潤率逐年提高；田園型海洋牧場平均收入和平均利潤逐年提高，平均收入利潤率比較穩(wěn)定；裝備型海洋牧場由2017—2018年連續(xù)兩年的虧損，2019年開始盈利，發(fā)展良好。按年度來看，底播型海洋牧場平均收入和平均利潤占同年40處省級(jí)海洋牧場平均收入和平均利潤總值的比例最高，裝備型海洋牧場平均收入和平均利潤占同年40處海洋牧場平均收入和平均利潤的比例最低。由此可見，不同類型海洋牧場經(jīng)營效益存在較大差異［15－18］。

表5 40處海洋牧場不同類型收入、利潤和利潤率情況表Tab.5 Income，profit and profit margin of different types of 40 marine ranches

以40處海洋牧場2017—2019年平均收入、平均利潤為基礎(chǔ)，利用SPSS 24.0軟件對(duì)40處海洋牧場進(jìn)行聚類分析，得到各個(gè)海洋牧場的分類表［19］。聚類分析是將分類對(duì)象置于一個(gè)多維空間中，按照它們空間關(guān)系的親疏程度進(jìn)行分類。通俗地講，聚類分析就是根據(jù)事物彼此不同的屬性進(jìn)行辨認(rèn)，將具有相似屬性的事物聚為一類，使得同一類的事物具有高度的相似性。聚類分析主要有二階分析、K中心聚類和層次聚類。本文采用層次聚類方法，層次聚類法（hierarchical clustering method）基本思想是：先將各個(gè)樣本各自看成一類，然后規(guī)定樣本之間的距離和類與類之間的距離；然后選擇距離最小的一對(duì)并成一個(gè)新類，計(jì)算新類和其他類的距離；再將距離最小的兩類合并，這樣每次減少一類，直至所有的樣本都成為一類為止。40處海洋牧場具體分類情況見圖1。

圖1直觀反映了40處海洋牧場逐步合并的過程。如果把40處海洋牧場分成兩類，第27處底播型海洋牧場和第30處田園型海洋牧場為一類，其余海洋牧場為一類；如果把40處海洋牧場分成三類，第27處底播型海洋牧場為一類，第30處田園型海洋牧場為一類，其余海洋牧場為一類；如果把40處海洋牧場分成四類，第23處投礁型海洋牧場和第26處底播型海洋牧場為一類，第27處底播型海洋牧場為一類，第30處田園型海洋牧場為一類，其他海洋牧場為一類；如果分成五類或者更多類，分類情況可以從圖中看到。由此可見，同種類型海洋牧場經(jīng)營情況也存在一定差異［20］。

圖1 聚類分析樹狀圖Fig.1 Tree diagram of cluster analysis

3 討論

近年來，為應(yīng)對(duì)日益嚴(yán)重的海洋生態(tài)荒漠化問題，促進(jìn)傳統(tǒng)漁業(yè)轉(zhuǎn)型升級(jí)，建設(shè)海洋牧場正成為我國海洋發(fā)展的重要選擇，山東省作為海洋大省，發(fā)展海洋牧場具有得天獨(dú)厚的自然條件。通過開展海洋牧場經(jīng)濟(jì)效益的實(shí)證分析，可以為山東省海洋牧場發(fā)展提供一定的參考。

第一，海洋牧場海域面積對(duì)經(jīng)濟(jì)效益的影響。海域面積可能會(huì)對(duì)海洋牧場經(jīng)濟(jì)效益產(chǎn)生影響，但并不是必然的要素。一方面部分用海單位為享受海洋牧場項(xiàng)目補(bǔ)貼，盲目申請(qǐng)大面積用海，審批下來后卻對(duì)海域面積開發(fā)利用不徹底，海域資源利用率不高，造成海洋牧場粗放式發(fā)展；另一方面由于海洋牧場用海面積較大，迅速發(fā)展的海洋牧場勢必會(huì)占用大量海域空間，沒有科學(xué)依據(jù)的建設(shè)海洋牧場會(huì)擠占其他行業(yè)用海，容易造成海域資源浪費(fèi)［21］。因此，需要建立科學(xué)合理的海洋牧場用海面積標(biāo)準(zhǔn)，提高海洋牧場海域資源利用率，實(shí)現(xiàn)海域資源的集約使用和可持續(xù)發(fā)展，提升海洋牧場經(jīng)濟(jì)效益，推動(dòng)海洋牧場現(xiàn)代化發(fā)展。

第二，海洋牧場類型對(duì)經(jīng)濟(jì)效益的影響。不同類型海洋牧場有不同的功能定位和發(fā)展方向，經(jīng)營情況不平衡，經(jīng)濟(jì)效益有差距；相同類型海洋牧場有不同的管理模式和資源稟賦，經(jīng)營情況也不平衡，經(jīng)濟(jì)效益有差異。因此，海洋牧場建設(shè)應(yīng)科學(xué)規(guī)劃，充分發(fā)揮自然條件、科研技術(shù)資源等優(yōu)勢，因地制宜，合理發(fā)展不同類型海洋牧場，從而提升海洋牧場經(jīng)濟(jì)效益［22］，實(shí)現(xiàn)山東省沿海城市海洋牧場均衡發(fā)展。