趙 昕，雷 亮，彭 楠，丁黎黎

（1.中國(guó)海洋大學(xué) 經(jīng)濟(jì)學(xué)院，山東 青島 266100；2.西南財(cái)經(jīng)大學(xué) 工商管理學(xué)院，成都 611130）

0 引言

從傳統(tǒng)的觀點(diǎn)看，狹義的海洋漁業(yè)指包括海水養(yǎng)殖、海洋捕撈在內(nèi)的生產(chǎn)活動(dòng)，屬于海洋第一產(chǎn)業(yè)范疇，是海洋經(jīng)濟(jì)的重要組成部分，也是我國(guó)四大傳統(tǒng)海洋產(chǎn)業(yè)之一。盡管廣義的海洋漁業(yè)還包括水產(chǎn)育苗、漁業(yè)工業(yè)、休閑漁業(yè)等在內(nèi)的上下游產(chǎn)業(yè)，但是水產(chǎn)品的捕撈、養(yǎng)殖等仍是海洋漁業(yè)基礎(chǔ)的和主導(dǎo)的生產(chǎn)活動(dòng)，因此本文對(duì)于海洋漁業(yè)的討論限于狹義的海洋漁業(yè)生產(chǎn)領(lǐng)域。限定狹義的海洋漁業(yè)范疇不僅有利于具體分析海洋漁業(yè)主體經(jīng)濟(jì)狀況，也避免了因海洋漁業(yè)相關(guān)產(chǎn)業(yè)歸屬的模糊而帶來(lái)的海洋漁業(yè)數(shù)據(jù)統(tǒng)計(jì)口徑不一致的問(wèn)題。就狹義的海洋漁業(yè)而言，我國(guó)海洋漁業(yè)發(fā)展尚存在以下問(wèn)題：近海水域環(huán)境污染嚴(yán)重，一些漁場(chǎng)常遭到赤潮等海洋災(zāi)害的侵襲[1]；漁業(yè)科研投入少，技術(shù)成果轉(zhuǎn)化率低[2,3]；漁船、網(wǎng)具等生產(chǎn)設(shè)備老化[4]。這些因素導(dǎo)致了漁業(yè)生產(chǎn)效率的低下。為了提高海洋漁業(yè)生產(chǎn)效率，多數(shù)學(xué)者都認(rèn)識(shí)到了海洋漁業(yè)轉(zhuǎn)型升級(jí)的重要性和緊迫性[5,6]，這也是本文的出發(fā)點(diǎn)。

為了研究省內(nèi)海洋漁業(yè)發(fā)展的區(qū)域差異，特別是區(qū)縣一級(jí)的微觀海洋漁業(yè)的投入產(chǎn)出在技術(shù)效率、規(guī)模報(bào)酬和要素利用方面的具體差異，本文選取代表性省份——山東省作為研究對(duì)象，通過(guò)分析省內(nèi)不同沿海地區(qū)海洋漁業(yè)的投入產(chǎn)出情況，利用DEA方法測(cè)算海洋漁業(yè)發(fā)展的區(qū)域差異，進(jìn)而為地方政府制定海洋漁業(yè)轉(zhuǎn)型政策提供參考。

1 模型設(shè)定與指標(biāo)選取

1.1 模型設(shè)定

數(shù)據(jù)包絡(luò)分析（DEA）是研究不同生產(chǎn)部門(mén)，又稱(chēng)決策單元（Decision Making Unit，DMU）相對(duì)有效性的非參數(shù)評(píng)價(jià)方法。技術(shù)效率、規(guī)模報(bào)酬、要素利用率是利用DEA方法評(píng)價(jià)地區(qū)海洋漁業(yè)投入產(chǎn)出效率差異的三個(gè)方面。技術(shù)效率（TE）是對(duì)決策單元的資源配置能力、使用效率等多方面能力的綜合評(píng)價(jià)[7]，又可分為純技術(shù)效率（PTE）和規(guī)模效率（SE）。PTE衡量了DMU處在最優(yōu)規(guī)模時(shí)投入要素的生產(chǎn)效率，SE衡量了DMU實(shí)際生產(chǎn)規(guī)模與最優(yōu)規(guī)模的差距，且有關(guān)系式：TE=PTE×SE。TE、PTE、SE的取值均在0和1之間，值越大，表示效率越高，當(dāng)其值等于1時(shí)被分別被稱(chēng)為技術(shù)有效、純技術(shù)有效、規(guī)模有效。利用CCR模型和BCC模型可求解TE、PTE和SE。CCR模型原理如下：設(shè)n個(gè)具有相同任務(wù)、環(huán)境和目標(biāo)的DMU，其輸入向量X=(x1，x2…xm)T,輸出向量為Y=(y1，y2…ys)T。這些DMU構(gòu)成了一個(gè)參考集，在生產(chǎn)可能性集滿(mǎn)足凸性、錐性、無(wú)效性和最小性的假定下，以被評(píng)價(jià)的DMU的效率θ*為目標(biāo)函數(shù)，以參考集中全部DMU的效率系數(shù)為約束條件，構(gòu)建投入方向的CCR模型如下：

CCR、BCC、SBM模型是DEA方法在海洋漁業(yè)投入產(chǎn)出應(yīng)用的基礎(chǔ)，然而它們沒(méi)有剔除環(huán)境因素對(duì)漁業(yè)生產(chǎn)效率的影響，也沒(méi)有把海洋漁業(yè)生產(chǎn)前沿面的隨機(jī)性納入考慮，因此需要對(duì)上述模型中可能存在的管理無(wú)效率進(jìn)行檢驗(yàn)和修正，通常采用三階段DEA方法[12]來(lái)處理：第一階段利用傳統(tǒng)CCR、BCC、SBM模型計(jì)算投入冗余。第二階段利用第一階段求出的各要素投入冗余，構(gòu)建隨機(jī)前沿模型（SFA），估計(jì)管理無(wú)效率和隨機(jī)噪聲，并剔除掉原投入數(shù)據(jù)中的掉環(huán)境因素和隨機(jī)因素，第三階段將修正后的數(shù)據(jù)利用CCR、BCC、SBM模型再次求解。第一、三階段與傳統(tǒng)DEA應(yīng)用無(wú)異，不再贅述。第二階段原理如下，首先構(gòu)建SFA回歸模型：

其中，Yij表示第i個(gè)決策單元中第j種要素投入冗余，Z=(z1，z2…zn)T表示環(huán)境變量向量，β=(β1，β2…βn)T表示回歸參數(shù)向量。Uji～N(0，σu)表示管理無(wú)效率項(xiàng)，一般為截?cái)喾植?。Vji～N(0，σu)表示隨機(jī)噪聲。再利用極大似然法估計(jì)SFA模型中的回歸參數(shù)β、混合誤差σ0、變差率γ。再利用管理無(wú)效率公式[12]計(jì)算Uji和Vji：

其中，η是環(huán)境項(xiàng)和隨機(jī)項(xiàng)的方差比值，εi=Uji+Vji

其中，Xji表示第i個(gè)決策單元中第j種投入，X*ji表示調(diào)整后的新投入，第一個(gè)括號(hào)里面代表了環(huán)境調(diào)整量，第二個(gè)括號(hào)里面代表了隨機(jī)誤差項(xiàng)的調(diào)整量。通過(guò)剔除環(huán)境和隨機(jī)誤差因素，即可使得所有DMU面臨相同的外部環(huán)境。

1.2 指標(biāo)選取

產(chǎn)出方面，選取海洋漁業(yè)增加值作為產(chǎn)出指標(biāo)。該指標(biāo)去除了海洋漁業(yè)總產(chǎn)值中多次重復(fù)計(jì)算的轉(zhuǎn)移價(jià)值，能夠反映海洋漁業(yè)經(jīng)濟(jì)活動(dòng)在一定時(shí)期內(nèi)所創(chuàng)造的價(jià)值。投入方面，選取海洋漁業(yè)從業(yè)人員和海水養(yǎng)殖面積作為“勞動(dòng)”和“土地”投入。在“資本”方面，一些研究將漁業(yè)固定資產(chǎn)投資作為“資本”投入[7]，本文未予采納。一是統(tǒng)計(jì)數(shù)據(jù)沒(méi)有統(tǒng)計(jì)該指標(biāo)，二是投資的滯后效應(yīng)，三是固定投資是資本增量，其對(duì)產(chǎn)出的貢獻(xiàn)相對(duì)于資本存量而言要小得多。漁船作為海洋漁業(yè)最主要的固定資本，直接體現(xiàn)了海洋漁業(yè)，尤其是遠(yuǎn)洋漁業(yè)的作業(yè)能力和裝配水準(zhǔn)。因此，相關(guān)文獻(xiàn)改用漁船數(shù)或漁船總動(dòng)力作為“資本”投入[13]?？紤]到漁船大小不一，性能各異，采用漁船數(shù)量難以反映漁船的捕撈能力，故選取漁船總動(dòng)力作為“資本”投入。在環(huán)境變量方面，目前的一些文獻(xiàn)關(guān)于環(huán)境指標(biāo)的選取標(biāo)準(zhǔn)沒(méi)有明確的規(guī)定。一般選取影響海洋漁業(yè)生產(chǎn)活動(dòng)但不受海洋漁業(yè)生產(chǎn)活動(dòng)控制的重要外生變量。本文從能夠顯著影響到海洋漁業(yè)產(chǎn)出和三種投入要素出發(fā)，充分權(quán)衡數(shù)據(jù)可獲得性和代表性，最終選取人均GDP（影響海洋漁業(yè)產(chǎn)出水平）、人均土地面積（制約海洋漁業(yè)生產(chǎn)空間）、人均固定資產(chǎn)投資額（影響海洋漁業(yè)投入水平）作為環(huán)境變量。本文采用2010年山東省7個(gè)沿海地區(qū)共35個(gè)沿海區(qū)/縣統(tǒng)計(jì)數(shù)據(jù)，數(shù)據(jù)來(lái)源于各地統(tǒng)計(jì)局。表示混合誤差項(xiàng)，φ(.)和φ(.)表示標(biāo)準(zhǔn)正態(tài)分布的分布函數(shù)和密度函數(shù)。最后利用調(diào)整公式[12]對(duì)原投入進(jìn)行同質(zhì)化調(diào)整：

2 山東省海洋漁業(yè)投入產(chǎn)出區(qū)域差異

2.1 準(zhǔn)備階段的數(shù)據(jù)預(yù)處理

第一階段，在求投入冗余之前，需要保證無(wú)異常數(shù)據(jù)影響結(jié)果的可靠性，可先通過(guò)刀切法（Jackknifing method）來(lái)檢驗(yàn)CCR、BCC、SBM模型的穩(wěn)健性：逐一剔除DEA模型中有效DMU，重新測(cè)度參考集中其他DMU的效率值并排序，將DMU新的效率值與原模型比對(duì)，測(cè)算次序值間的相關(guān)系數(shù)。結(jié)果顯示，相關(guān)系數(shù)均很高，表明去掉有效DMU的檢驗(yàn)?zāi)Ｐ蛯?duì)DMU排序影響很小，可認(rèn)為無(wú)異常數(shù)據(jù)。

第二階段，利用第一階段計(jì)算的投入冗余，估計(jì)SFA回歸方程。由表1可知，9個(gè)SFA回歸方程的似然比值均在1%的置信度下通過(guò)檢驗(yàn)，表明半正態(tài)分布假設(shè)適用于本文選用的SFA模型。變差率γ均趨于1，說(shuō)明原DEA模型存在較大的管理無(wú)效率，需要剔除環(huán)境因素和隨機(jī)因素；除海洋漁業(yè)從業(yè)人員投入冗余在BBC模型下的SFA回歸方程未通過(guò)參數(shù)檢驗(yàn)，其余8個(gè)模型全部回歸參數(shù)均在1%和5%置信度下通過(guò)檢驗(yàn)，其中絕大部分回歸參數(shù)在1%置信度下通過(guò)檢驗(yàn)。這表明，本文所選用的環(huán)境變量對(duì)海洋漁業(yè)各投入要素冗余有顯著影響。

表1 三種DEA模型第二階段SFA方程回歸結(jié)果

但為穩(wěn)妥起見(jiàn)，本文仍不使用BCC模型下的SFA回歸方程來(lái)調(diào)整原始數(shù)據(jù)，而選用經(jīng)CCR-SFA回歸方程調(diào)整后的數(shù)據(jù)測(cè)算技術(shù)效率和規(guī)模報(bào)酬，選用經(jīng)SBM-SFA回歸方程調(diào)整后的數(shù)據(jù)測(cè)算要素利用率。第三階段輸出結(jié)果如表2所示。

2.2 技術(shù)效率的區(qū)域差異

在技術(shù)效率方面，TE、SE區(qū)域差異性分布較為明顯，且兩者呈現(xiàn)協(xié)同分布。TE區(qū)域差異性分布主要表現(xiàn)在：青島即墨、膠南兩市TE值分別高于均值62.64%和112.71%，而黃島、嶗山兩區(qū)僅為均值的23.89%和52.98%。煙臺(tái)市11區(qū)縣中，有兩市TE值為1，兩市略高于均值，其余7地均低于均值，其中最低的芝罘區(qū)僅為均值的35.86%。威海有三地TE值高于均值，兩地TE值低于均值。東營(yíng)市河口區(qū)TE值為1，而東營(yíng)區(qū)、利津縣、廣饒縣TE值低于均值。濰坊、日照、濱州轄內(nèi)絕大部分地區(qū)TE值都比較低。SE區(qū)域差異性分布主要表現(xiàn)在：青島市即墨、膠南兩市SE值分別高于均值78.62%和94.67%，而黃島、嶗山兩區(qū)僅為均值的22.16%和49.99%。煙臺(tái)有6地SE值高于均值，其中長(zhǎng)島和蓬萊SE值為1，另外5地SE值低于均值，最低的萊山區(qū)僅為均值的0.8%。東營(yíng)市河口區(qū)SE值為1，而東營(yíng)區(qū)、利津縣、廣饒縣SE值低于均值。威海經(jīng)發(fā)區(qū)SE值僅為均值的30.5%，其余四地SE值高于均值，最高的乳山市SE值高于均值58.20%。濰坊、日照、濱州轄內(nèi)絕大部分地區(qū)SE值都比較低，且大部分低于均值。PTE分布較為集中，其均值為0.727，顯著高于TE、SE均值，方差為0.044，顯著小于TE、PTE方差。

表2 技術(shù)效率、規(guī)模報(bào)酬、要素利用率輸出結(jié)果

綜合TE、PTE、SE的分布來(lái)看，可以將山東省沿海地區(qū)海洋漁業(yè)的發(fā)展模型分為三類(lèi)。一是兩級(jí)分化類(lèi)：這些地區(qū)轄內(nèi)區(qū)縣部分技術(shù)有效，處在最優(yōu)生產(chǎn)規(guī)模上，部分技術(shù)落后、管理欠缺，包括青島、煙臺(tái)、威海。二是一枝獨(dú)秀類(lèi)：這些地區(qū)轄內(nèi)少數(shù)區(qū)縣發(fā)展較好，多數(shù)區(qū)縣發(fā)展落后，代表地區(qū)是東營(yíng)。三是發(fā)展落后類(lèi)：這些地區(qū)轄內(nèi)區(qū)縣都處于技術(shù)無(wú)效和規(guī)模無(wú)效狀況，包括濰坊、日照、濱州。為了縮小海洋漁業(yè)技術(shù)效率在區(qū)域間的差異，促進(jìn)區(qū)域海洋漁業(yè)協(xié)調(diào)發(fā)展，應(yīng)根據(jù)各地的實(shí)際情況提出改善意見(jiàn)。針對(duì)兩級(jí)分化類(lèi)地區(qū)，要提高發(fā)展滯后地區(qū)勞動(dòng)力的待遇，避免人才用腳投票，同時(shí)鼓勵(lì)勞動(dòng)力流動(dòng)，建立一對(duì)一幫扶政策，加強(qiáng)分化地區(qū)的交流與合作，促進(jìn)周?chē)貐^(qū)漁業(yè)技術(shù)普及和管理經(jīng)驗(yàn)的分享，縮小地區(qū)間技術(shù)、管理方面的差距。針對(duì)一枝獨(dú)秀類(lèi)地區(qū)，這些地區(qū)少數(shù)發(fā)展較好的區(qū)縣可以看作是區(qū)域增長(zhǎng)極，應(yīng)增強(qiáng)增長(zhǎng)極的擴(kuò)散效應(yīng)。具體而言，在周?chē)貐^(qū)適量擴(kuò)大海水養(yǎng)殖面積、放寬海水捕撈管理政策，促進(jìn)生產(chǎn)要素的擴(kuò)散。針對(duì)發(fā)展落后類(lèi)，應(yīng)集中有限的資金、人才、政策等資源，重點(diǎn)培育具有地方特色和優(yōu)勢(shì)的海洋漁業(yè)發(fā)展方向和較優(yōu)秀的相關(guān)企業(yè)，另外可以在海洋漁業(yè)的下游產(chǎn)業(yè)上尋找突破口，通過(guò)政策優(yōu)惠吸引投資建立水產(chǎn)品加工廠、漁文化產(chǎn)業(yè)園，擴(kuò)大下游市場(chǎng)容量，增進(jìn)產(chǎn)業(yè)聯(lián)動(dòng)，提高海洋漁業(yè)抗風(fēng)險(xiǎn)能力。

2.3 規(guī)模報(bào)酬與要素利用率

在規(guī)模報(bào)酬方面，從RTS統(tǒng)計(jì)分布看，處在規(guī)模報(bào)酬不變的地區(qū)有3個(gè)，占總數(shù)8.57%，其中煙臺(tái)市占2個(gè)，分別為長(zhǎng)島縣和蓬萊市，另一個(gè)是東營(yíng)河口區(qū)。處在規(guī)模報(bào)酬遞減的地區(qū)共2個(gè)，占總數(shù)5.71%，分別為青島膠南市、煙臺(tái)龍口市、威海榮成市。處在規(guī)模報(bào)酬遞增的地區(qū)共30個(gè)，占總數(shù)85.71%。這表明山東省超過(guò)五分之四的地區(qū)海洋漁業(yè)生產(chǎn)規(guī)模仍存在增加的空間。在要素利用率方面，漁船總動(dòng)力利用率（UTP）均值為44.85%，UTP值在均值以下的共計(jì)22個(gè)，占總數(shù)62.86%。UTP值在80%（高利用率）以上共7個(gè)，僅占總數(shù)五分之一。海洋漁業(yè)從業(yè)人員利用率（UFW）均值68.83%，UFW值在均值以下的共計(jì)17個(gè)，占總數(shù)約一半，UFW值在80%以上共12個(gè)，占比34.29%。海水養(yǎng)殖面積利用率（UAA）均值60.30%，UAA值在均值以下的共計(jì)19個(gè)，占比54.30%。而UAA在80%以上共12個(gè)，占比34.29%。由此可見(jiàn)，三種投入要素利用率均不高，其中，漁船總動(dòng)力平均利用率最低，這與漁船設(shè)備老化、環(huán)渤海地區(qū)漁業(yè)資源枯竭有關(guān)。其次是海水養(yǎng)殖面積利用率，長(zhǎng)期以來(lái)，山東省海洋漁業(yè)養(yǎng)殖布局不合理，高密度養(yǎng)殖和魚(yú)藥濫用現(xiàn)象突出，水環(huán)境污染嚴(yán)重。海洋漁業(yè)從業(yè)人員平均利用率較另外兩種投入高些，這與近年來(lái)漁業(yè)人才素質(zhì)不斷提高有關(guān)。

綜合規(guī)模報(bào)酬和要素利用率可以看出，山東省大多數(shù)地區(qū)海洋漁業(yè)發(fā)展存在規(guī)模報(bào)酬遞增和要素利用不足的現(xiàn)象，表明其海洋漁業(yè)產(chǎn)出水平處在長(zhǎng)期平均總成本曲線(xiàn)的遞減階段，仍存在較大發(fā)展空間，應(yīng)提高要素利用率使落在遞減階段的區(qū)縣前移至規(guī)模報(bào)酬不變的位置。針對(duì)海水養(yǎng)殖面積利用率，要控制養(yǎng)殖規(guī)模，優(yōu)化養(yǎng)殖結(jié)構(gòu)。具體而言，在養(yǎng)殖規(guī)模上，要確定海水養(yǎng)殖容量，對(duì)海水養(yǎng)殖總量和密度進(jìn)行控制，減少魚(yú)藥和飼料的濫用，防止養(yǎng)殖環(huán)境的惡化，提高水產(chǎn)品的質(zhì)量。在養(yǎng)殖結(jié)構(gòu)方面，可構(gòu)建立體多功能的海洋養(yǎng)殖區(qū)域，發(fā)展生態(tài)養(yǎng)殖，釋放海水養(yǎng)殖潛力；同時(shí)對(duì)單一的海水養(yǎng)殖管理模式進(jìn)行革新，推動(dòng)政府、行業(yè)、漁民協(xié)同共治，避免出現(xiàn)“公地悲劇”的情況。在勞動(dòng)力方面，我國(guó)漁業(yè)從業(yè)人員的平均受教育年限僅為7.8年，知識(shí)水平較低，可定期進(jìn)行漁業(yè)技術(shù)培訓(xùn)、舉辦漁業(yè)知識(shí)交流會(huì)，提升從業(yè)人員的操作技能和專(zhuān)業(yè)素養(yǎng)，實(shí)現(xiàn)知識(shí)的正外部性效應(yīng)。在資本設(shè)備方面，目前較多作業(yè)漁船接近或已達(dá)到報(bào)廢年限，可通過(guò)政府給予漁民補(bǔ)貼推動(dòng)老舊漁船的更新?lián)Q代；另一方面，應(yīng)提高漁船的裝備水平，完善漁港物資供應(yīng)船舶維修等經(jīng)營(yíng)性服務(wù)設(shè)施。

3 結(jié)論

本文探討了DEA方法在漁業(yè)投入產(chǎn)出中的應(yīng)用，利用三階段CCR、BCC、SBM模型，彌補(bǔ)了傳統(tǒng)DEA方法在漁業(yè)中應(yīng)用的不足。從技術(shù)效率、規(guī)模報(bào)酬、要素利用三個(gè)角度，實(shí)證測(cè)算了2010年山東省沿海地區(qū)海洋漁業(yè)的區(qū)域差異。綜合TE、PTE、SE分布情況，山東省沿海地區(qū)海洋漁業(yè)發(fā)展方式可分為三類(lèi)：一是兩極分化類(lèi)，代表地區(qū)為青島、煙臺(tái)、威海；二是一枝獨(dú)秀類(lèi)，代表地區(qū)為東營(yíng)；三是發(fā)展落后類(lèi)，代表地區(qū)為濰坊、日照、濱州。綜合產(chǎn)業(yè)規(guī)模報(bào)酬和要素利用率來(lái)看，山東省大多數(shù)地區(qū)海洋漁業(yè)發(fā)展同時(shí)存在資源利用不足和規(guī)模報(bào)酬遞增的問(wèn)題，其海洋漁業(yè)產(chǎn)出水平處在長(zhǎng)期平均總成本曲線(xiàn)的遞減階段，仍存在很大的發(fā)展空間，應(yīng)通過(guò)改善養(yǎng)殖模式，加強(qiáng)漁業(yè)從業(yè)人員專(zhuān)業(yè)素質(zhì)培養(yǎng)、提升漁業(yè)機(jī)械化裝備水平等方式，來(lái)提高要素利用率使產(chǎn)出前移至規(guī)模報(bào)酬不變的水平。