于謹(jǐn)凱，李思達(dá)

(中國(guó)海洋大學(xué)經(jīng)濟(jì)學(xué)院，山東青島 266100)

海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度的時(shí)空分異研究
——以沿海11省市為例

于謹(jǐn)凱，李思達(dá)

(中國(guó)海洋大學(xué)經(jīng)濟(jì)學(xué)院，山東青島 266100)

以中國(guó)沿海11省市2003－2013年面板數(shù)據(jù)為樣本，在引入海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)耦合協(xié)調(diào)綜合發(fā)展評(píng)價(jià)指標(biāo)體系的基礎(chǔ)上，利用海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境耦合協(xié)調(diào)度的時(shí)空分異模型，從時(shí)間和空間兩個(gè)維度對(duì)海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)耦合協(xié)調(diào)度進(jìn)行了研究。結(jié)果表明:從時(shí)序上看，我國(guó)沿海省際海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度總體上升顯著且有趨同的趨勢(shì)；從局部空間看，江蘇－浙江－上海區(qū)域協(xié)調(diào)度差異始終存在，上海市協(xié)調(diào)度的穩(wěn)步提升使得山東－江蘇－浙江－上海區(qū)域整體協(xié)調(diào)度提升至較高水平。據(jù)此，提出了加強(qiáng)海洋漁業(yè)布局規(guī)劃，充分發(fā)揮資源環(huán)境系統(tǒng)的自我凈化能力等對(duì)策，以期實(shí)現(xiàn)我國(guó)省際海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)及資源環(huán)境系統(tǒng)高度協(xié)調(diào)及省際差距的進(jìn)一步縮小。

海洋漁業(yè)；經(jīng)濟(jì)系統(tǒng)；協(xié)調(diào)度；時(shí)空分異

0 引 言

作為我國(guó)傳統(tǒng)漁業(yè)資源產(chǎn)業(yè)，海洋漁業(yè)近年來發(fā)展迅速。根據(jù)《中國(guó)海洋統(tǒng)計(jì)年鑒》的數(shù)據(jù)，2012年海洋漁業(yè)增加值達(dá)3 560.5億元，超過了2000年海洋漁業(yè)的總產(chǎn)值3 240.2億元?！度珖?guó)漁業(yè)發(fā)展第十三個(gè)五年規(guī)劃(2016－2020年)》[1]提出，“十三五”漁業(yè)發(fā)展要牢固樹立創(chuàng)新、協(xié)調(diào)、綠色、開放、共享發(fā)展理念，以提質(zhì)增效、減量增收、綠色發(fā)展、富裕漁民為目標(biāo)，堅(jiān)持生態(tài)優(yōu)先、創(chuàng)新驅(qū)動(dòng)、“走出去”戰(zhàn)略、以人為本、依法治漁，大力推進(jìn)漁業(yè)供給側(cè)結(jié)構(gòu)性改革，加快轉(zhuǎn)變漁業(yè)發(fā)展方式，加快實(shí)現(xiàn)漁業(yè)現(xiàn)代化。基于此，海洋漁業(yè)是否與其所處的生態(tài)環(huán)境相適應(yīng)、相協(xié)調(diào)理應(yīng)得到重點(diǎn)關(guān)注。為明晰該海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)與資源環(huán)境系統(tǒng)間的協(xié)調(diào)度在時(shí)空維度上的分布情況，進(jìn)而研究如何促進(jìn)海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)與資源環(huán)境系統(tǒng)相互協(xié)調(diào)，本文在求解該協(xié)調(diào)度的基礎(chǔ)上，進(jìn)一步利用空間相關(guān)分析等方法對(duì)海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度的時(shí)空分異進(jìn)行研究，并提出相應(yīng)的對(duì)策建議。

1 文獻(xiàn)綜述

1.1 關(guān)于經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度的研究

經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)的定量研究主要集中于相關(guān)區(qū)域與環(huán)境系統(tǒng)，或是產(chǎn)業(yè)經(jīng)濟(jì)系統(tǒng)與環(huán)境系統(tǒng)的互動(dòng)影響，研究的形式主要包括可持續(xù)發(fā)展評(píng)價(jià)、綜合發(fā)展評(píng)價(jià)、生態(tài)或環(huán)境效率評(píng)價(jià)、測(cè)算經(jīng)濟(jì)系統(tǒng)與環(huán)境系統(tǒng)的協(xié)調(diào)度等類型。從區(qū)域經(jīng)濟(jì)與環(huán)境協(xié)調(diào)的角度來看，國(guó)內(nèi)研究主要著眼于對(duì)區(qū)域內(nèi)部經(jīng)濟(jì)整體和環(huán)境協(xié)調(diào)現(xiàn)狀進(jìn)行測(cè)度，分析其時(shí)空差異性[2]，并對(duì)未來的協(xié)調(diào)度進(jìn)行預(yù)測(cè)[3]；一些學(xué)者進(jìn)一步引入能源系統(tǒng)，測(cè)算區(qū)經(jīng)濟(jì)、環(huán)境、能源形成的3E系統(tǒng)協(xié)調(diào)度，利用空間計(jì)量方法分析區(qū)間協(xié)調(diào)度空間交互效應(yīng)及分布特征[4]。此外，一些國(guó)內(nèi)研究者對(duì)產(chǎn)業(yè)經(jīng)濟(jì)系統(tǒng)與環(huán)境是否協(xié)調(diào)進(jìn)行測(cè)度、預(yù)測(cè)并分析其時(shí)序和空間分布特征[5]。具體到海洋經(jīng)濟(jì)及相關(guān)產(chǎn)業(yè)系統(tǒng)，國(guó)內(nèi)學(xué)者多致力于海洋環(huán)境評(píng)價(jià)指標(biāo)體系的構(gòu)建，如:從海洋可持續(xù)發(fā)展生產(chǎn)力、科技和生態(tài)環(huán)境三個(gè)方面對(duì)山東省2000－2008年海洋經(jīng)濟(jì)可持續(xù)發(fā)展能力進(jìn)行測(cè)度[6]；以壓力－狀態(tài)－響應(yīng)(P－S－R)為主線將海水養(yǎng)殖生態(tài)健康評(píng)價(jià)系統(tǒng)下的25個(gè)指標(biāo)分別歸類于系統(tǒng)壓力、系統(tǒng)狀態(tài)、系統(tǒng)響應(yīng)三個(gè)子系統(tǒng)[7]；建立我國(guó)海水養(yǎng)殖的生態(tài)預(yù)警評(píng)價(jià)體系對(duì)我國(guó)2001－2009年海水養(yǎng)殖業(yè)可持續(xù)發(fā)展程度進(jìn)行測(cè)度，并得出我國(guó)海水養(yǎng)殖業(yè)生態(tài)系統(tǒng)與經(jīng)濟(jì)社會(huì)發(fā)展已出現(xiàn)了一定程度上的不協(xié)調(diào)[8]；利用狀態(tài)空間法構(gòu)建海洋承載力評(píng)價(jià)體系，并對(duì)我國(guó)2003－2009年海域承載力進(jìn)行測(cè)算[9]。在指標(biāo)體系建立的基礎(chǔ)上，國(guó)內(nèi)學(xué)者從不同角度對(duì)海洋經(jīng)濟(jì)與環(huán)境間的協(xié)調(diào)進(jìn)行研究，如:驗(yàn)證海洋生態(tài)系統(tǒng)在其經(jīng)濟(jì)、環(huán)境、社會(huì)三個(gè)子系統(tǒng)的交互脅迫作用下呈“S型發(fā)展機(jī)制”[10]；對(duì)我國(guó)海水養(yǎng)殖與海洋生態(tài)環(huán)境的協(xié)調(diào)度進(jìn)行求證，得出二者并非總是協(xié)調(diào)的結(jié)論[11]；對(duì)海洋產(chǎn)業(yè)空間集聚及其協(xié)調(diào)發(fā)展進(jìn)行測(cè)度，得出我國(guó)海洋產(chǎn)業(yè)已基本形成集聚區(qū)，且產(chǎn)業(yè)系統(tǒng)與環(huán)境資源系統(tǒng)基本協(xié)調(diào)[12]；從時(shí)間和空間兩個(gè)角度測(cè)度了我國(guó)海洋經(jīng)濟(jì)與環(huán)境的脫鉤效應(yīng)。結(jié)果顯示，海洋經(jīng)濟(jì)增長(zhǎng)與海洋環(huán)境壓力系統(tǒng)間尚未達(dá)到強(qiáng)脫鉤狀態(tài)，海洋環(huán)境系統(tǒng)仍然面臨壓力[13]。

1.2 關(guān)于時(shí)空分異的研究現(xiàn)狀

傳統(tǒng)的時(shí)空差異的研究主要是對(duì)加權(quán)的變異系數(shù)、錫爾系數(shù)、基尼系數(shù)計(jì)算并分析[14]；傳統(tǒng)的馬爾科夫鏈法常用于時(shí)序問題的分析、預(yù)測(cè)，一些學(xué)者對(duì)其修正，運(yùn)用空間馬爾科夫轉(zhuǎn)移分析法來對(duì)囊括時(shí)空因素的問題加以研究，如利用馬爾科夫鏈分析法描述中國(guó)外商直接投資區(qū)位分布的空間演變過程[15]；ESDA，即探索性空間數(shù)據(jù)分析，其核心為空間自相關(guān)分析。近年來，國(guó)內(nèi)外學(xué)者多使用ESDA來對(duì)數(shù)據(jù)空間分布上的關(guān)聯(lián)性、異質(zhì)性、影響因素、驅(qū)動(dòng)機(jī)制進(jìn)行研究，如對(duì)長(zhǎng)三角城市群R＆D投入相關(guān)數(shù)據(jù)分析，得出該區(qū)域R＆D投入的空間分異及演化特征[16]。國(guó)內(nèi)學(xué)者主要利用ESDA對(duì)效率評(píng)價(jià)、綜合發(fā)展評(píng)價(jià)等指標(biāo)的區(qū)差異進(jìn)行分析。實(shí)踐中，Moran’s I和熱點(diǎn)分析常被用來描述空間全局和局部自相關(guān)狀態(tài)和程度，如利用Moran’s I和熱點(diǎn)分析對(duì)中國(guó)級(jí)以上城市四化協(xié)調(diào)度的空間自相關(guān)加以分析，并運(yùn)用重心移動(dòng)來反映四化協(xié)調(diào)度的結(jié)構(gòu)變化[17]；利用空間自相關(guān)分析得出浙江海洋漁業(yè)的空間集聚趨勢(shì)較為顯著，中心－外圍的空間布局模式已經(jīng)形成[18]；利用空間自相關(guān)模型和重心遷移模型分析長(zhǎng)三角土地集約利用和經(jīng)濟(jì)發(fā)展的空間差異特征以及空間相關(guān)性[19]；利用R/S分型法對(duì)泛長(zhǎng)江城市發(fā)展效率時(shí)序特征的持續(xù)性進(jìn)行進(jìn)一步研究，并得出結(jié)論:泛長(zhǎng)江城市發(fā)展效率在未來一段時(shí)間內(nèi)將較穩(wěn)定沿襲以往的趨勢(shì)[20]。

綜上所述，以往協(xié)調(diào)度的研究，或僅從時(shí)序角度研究，或主要從空間角度研究，而時(shí)空角度同時(shí)進(jìn)行的研究較少；具體到海洋相關(guān)產(chǎn)業(yè)，前人研究多從海洋產(chǎn)業(yè)整體視角出發(fā)，對(duì)海洋產(chǎn)業(yè)與環(huán)境的協(xié)調(diào)關(guān)系探究，較少對(duì)單個(gè)產(chǎn)業(yè)集中研究；空間尺度的選擇上，研究多以一國(guó)整體或局部區(qū)域?yàn)橹黧w進(jìn)行探究，基于一國(guó)省際尺度的研究較少。所以，本文在借鑒前人成果的基礎(chǔ)之上，使用全局時(shí)序主成分分析法與耦合協(xié)調(diào)度模型對(duì)中國(guó)省際海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度進(jìn)行測(cè)度，并利用變異系數(shù)和ESDA等方法對(duì)其時(shí)序和空間特征進(jìn)行分析。

2 海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境耦合協(xié)調(diào)度的時(shí)空分異模型構(gòu)建

2.1 海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境耦合協(xié)調(diào)度時(shí)空分異的指標(biāo)體系構(gòu)建

海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)耦合協(xié)調(diào)度時(shí)空分異模型的構(gòu)建需要兼顧海洋漁業(yè)經(jīng)濟(jì)子系統(tǒng)和海洋資源環(huán)境子系統(tǒng)兩大系統(tǒng)，并體現(xiàn)出各子系統(tǒng)內(nèi)指標(biāo)間的互動(dòng)關(guān)系。因此，本文在考慮指標(biāo)的科學(xué)性及數(shù)據(jù)可得性的基礎(chǔ)上，借鑒前人經(jīng)驗(yàn)[8]，分別從經(jīng)濟(jì)規(guī)模、經(jīng)濟(jì)水平、經(jīng)濟(jì)結(jié)構(gòu)三個(gè)方面選取代表性指標(biāo)構(gòu)建海洋漁業(yè)經(jīng)濟(jì)子系統(tǒng)的綜合發(fā)展評(píng)價(jià)體系，并從海洋環(huán)境承載、海洋災(zāi)害、海洋環(huán)境治理與保護(hù)三個(gè)方面建立起海洋資源環(huán)境子系統(tǒng)的綜合發(fā)展評(píng)價(jià)體系。

表1 海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)耦合協(xié)調(diào)綜合發(fā)展評(píng)價(jià)指標(biāo)體系構(gòu)建

2.2 海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境耦合協(xié)調(diào)度的測(cè)度

耦合度的概念來源于物理學(xué)，是度量系統(tǒng)結(jié)構(gòu)中各要素間相互聯(lián)系緊密程度的一種方法。兩個(gè)或兩個(gè)以上系統(tǒng)在運(yùn)行過程中，由于系統(tǒng)的各要素間相互關(guān)聯(lián)及作用，導(dǎo)致系統(tǒng)彼此影響，最終實(shí)現(xiàn)系統(tǒng)間各要素緊密結(jié)合、相互依存。

海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)與其所在環(huán)境的相互作用擬合了耦合原理。就海洋漁業(yè)而言，其生產(chǎn)過程勢(shì)必對(duì)海洋資源有所消耗，并對(duì)海洋環(huán)境造成一定的影響；就海洋資源環(huán)境系統(tǒng)而言，海洋資源是否充裕，海洋環(huán)境態(tài)勢(shì)是否良好將對(duì)海洋漁業(yè)經(jīng)濟(jì)活動(dòng)產(chǎn)生不同的影響。可見海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)及海洋資源環(huán)境系統(tǒng)在運(yùn)行及融合過程中，各要素密切相關(guān)、相互影響，擬合耦合原理的機(jī)制。

因此，為了反映出海洋漁業(yè)經(jīng)濟(jì)與海洋環(huán)境兩個(gè)系統(tǒng)的協(xié)調(diào)互動(dòng)關(guān)系，本文采用物理學(xué)的容量耦合度概念構(gòu)建系統(tǒng)間耦合度函數(shù)[12]，具體見公式(1)，并利用全局主成分分析法對(duì)指標(biāo)進(jìn)行處理，并據(jù)此獲得兩個(gè)子系統(tǒng)的綜合評(píng)價(jià)值。

其中，S1、S2分別對(duì)應(yīng)海洋漁業(yè)經(jīng)濟(jì)子系統(tǒng)、海洋環(huán)境子系統(tǒng)的綜合評(píng)價(jià)值的歸一化值，C表示兩個(gè)子系統(tǒng)耦合度。C∈[0，1]，且C越接近1表明子共振耦合越好，子系統(tǒng)之間的結(jié)構(gòu)越有序。

耦合度雖然反映出了兩系統(tǒng)的協(xié)調(diào)程度，但無法區(qū)分系統(tǒng)處于協(xié)調(diào)水平的高低。因此，需要進(jìn)一步引入耦合協(xié)調(diào)模型來界定協(xié)調(diào)度水平的高低，對(duì)系統(tǒng)協(xié)調(diào)狀態(tài)做出更準(zhǔn)確的描述，具體見公式(2)[12]。

其中，D為耦合協(xié)調(diào)度，T為海洋漁業(yè)經(jīng)濟(jì)與海洋資源環(huán)境協(xié)調(diào)發(fā)展的綜合評(píng)價(jià)指數(shù)；a、b為待定系數(shù)，考慮到兩子系統(tǒng)的綜合發(fā)展具有相同的重要性，取a＝b＝0.5；F1、F2分別為海洋漁業(yè)經(jīng)濟(jì)子系統(tǒng)與海洋資源系統(tǒng)的綜合評(píng)價(jià)值。借鑒前人的相關(guān)研究，本文對(duì)海洋漁業(yè)環(huán)境適應(yīng)協(xié)調(diào)度所處狀態(tài)的界定見表2[3]。

表2 海洋漁業(yè)環(huán)境適應(yīng)協(xié)調(diào)度界定表

2.3 海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境耦合協(xié)調(diào)度的時(shí)空分異測(cè)度

本文通過變異系數(shù)來對(duì)海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境耦合協(xié)調(diào)度的時(shí)序分異進(jìn)行測(cè)度，并利用局部空間自相關(guān)檢驗(yàn)來對(duì)海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境耦合協(xié)調(diào)度的局部空間分異進(jìn)行測(cè)度。

2.3.1 變異系數(shù)

為明晰省際協(xié)調(diào)度的差異在時(shí)序上的變化，本文引入了變異系數(shù)，具體見公式(3)。

2.3.2 局部空間自相關(guān)檢驗(yàn)

局部自相關(guān)具體探求協(xié)調(diào)度在何處集聚，并采用局部Moran’s I測(cè)度。局部Moran’s I可以檢測(cè)某理單元與鄰近理單元的屬性相似度，本文借此來描述相鄰省份耦合協(xié)調(diào)度的關(guān)聯(lián)度，并用LISA圖對(duì)協(xié)調(diào)度的集聚進(jìn)行描述。局部Moran’s I計(jì)算式如下:

其中，xi、xj分別代表區(qū)域i、j的耦合協(xié)調(diào)度觀測(cè)值；n為區(qū)總數(shù)；wij為空間權(quán)值矩陣。Ii∈[－1，1]，Ii＞0表示空間正相關(guān)，Ii＞0表示空間負(fù)相關(guān)，Ii＝0表示空間不相關(guān)。

3 海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)耦合協(xié)調(diào)度的時(shí)空分異分析

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

本文以2003－2013年數(shù)據(jù)為例，對(duì)中國(guó)沿海11省、市海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度的時(shí)空分異進(jìn)行研究。本文指標(biāo)的數(shù)據(jù)來源于2004－2014年《中國(guó)漁業(yè)統(tǒng)計(jì)年鑒》(整理獲得海水養(yǎng)殖產(chǎn)值、海洋捕撈產(chǎn)值、海洋機(jī)動(dòng)漁船年末擁有量、海洋漁業(yè)從業(yè)人員以及漁業(yè)從業(yè)人員、漁民人均收入、水產(chǎn)技術(shù)推廣機(jī)構(gòu)人員經(jīng)費(fèi)、海水養(yǎng)殖總產(chǎn)量及其中藻類產(chǎn)量、海水養(yǎng)殖深水網(wǎng)箱養(yǎng)殖產(chǎn)量、海水養(yǎng)殖工廠化養(yǎng)殖產(chǎn)量等數(shù)據(jù))和2004－2014年《中國(guó)海洋統(tǒng)計(jì)年鑒》(整理獲得工業(yè)廢水直接入海量、海洋災(zāi)害導(dǎo)致經(jīng)濟(jì)損失、當(dāng)年廢水治理竣工項(xiàng)目數(shù)、海洋保護(hù)區(qū)面積等的相關(guān)數(shù)據(jù))。

指標(biāo)體系中，海洋漁業(yè)經(jīng)濟(jì)總產(chǎn)值以海洋捕撈與海水養(yǎng)殖總產(chǎn)值的和代替；海洋漁業(yè)勞動(dòng)力在區(qū)漁業(yè)勞動(dòng)力中占比為海洋漁業(yè)從業(yè)人員同漁業(yè)從業(yè)人員之比；海洋災(zāi)害導(dǎo)致經(jīng)濟(jì)損失以風(fēng)暴潮導(dǎo)致的經(jīng)濟(jì)損失來計(jì)；海洋漁業(yè)受災(zāi)面積由風(fēng)暴潮導(dǎo)致漁業(yè)受災(zāi)面積代替。個(gè)別數(shù)據(jù)在統(tǒng)計(jì)年鑒中有缺漏，本文以插值法求得，計(jì)算式為:xt＝(xt＋1＋xt－1)/2，xt為t年度的數(shù)據(jù)。

3.2 海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度的時(shí)序分異分析

3.2.1 變異系數(shù)的測(cè)算

為利用CV指數(shù)對(duì)海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度的時(shí)序分異測(cè)度，首先要計(jì)算海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)耦合度，進(jìn)而得到海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度。本文通過全局時(shí)序主成分分析法計(jì)算子系統(tǒng)綜合發(fā)展評(píng)價(jià)值，并基于此測(cè)度海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)耦合度。

在進(jìn)行主成分分析之前，本文首先對(duì)數(shù)據(jù)進(jìn)行KMO檢驗(yàn)和Bartlett球形檢驗(yàn)。檢驗(yàn)結(jié)果顯示，海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)與海洋環(huán)境資源系統(tǒng)的Bartlett球形檢驗(yàn)的顯著值均小于0.01；KMO檢驗(yàn)分別為0.723與0.636，皆大于0.60。以上結(jié)果說明，全局主成分分析法適于對(duì)本文的數(shù)據(jù)分析。繼而，用全局主成分法對(duì)兩子系統(tǒng)下的二級(jí)指標(biāo)分析并提取各子系統(tǒng)的主成分。根據(jù)累計(jì)貢獻(xiàn)率大于80%的原則，本文提取海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)的前三個(gè)主成分及海洋資源環(huán)境系統(tǒng)的前四個(gè)主成分。兩系統(tǒng)分別對(duì)應(yīng)的主成分、累計(jì)貢獻(xiàn)率、特征根如表3所示。

表3 主成分分析結(jié)果

各主成分評(píng)價(jià)綜合評(píng)價(jià)公式具體如下所示:

據(jù)上述公式，將計(jì)算得出的各主成分值根據(jù)貢獻(xiàn)率在所提取主成分的累計(jì)貢獻(xiàn)率中的占比確定權(quán)重，并得到兩系統(tǒng)綜合評(píng)價(jià)式(12)、(13)。

中國(guó)沿海11省市海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)耦合度見表4。由表4可知，我國(guó)省際海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)耦合度均有所提升。其中天津、河北、上海、廣西提升幅度較大，其海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)同資源環(huán)境系統(tǒng)在磨合中趨于協(xié)調(diào)。

表4 海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)耦合度

本文根據(jù)耦合度計(jì)算出海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度值，具體各省份協(xié)調(diào)度實(shí)證結(jié)果見表5。

表5 海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度

在表5所得結(jié)果的基礎(chǔ)上，2003－2013年中國(guó)海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度年度變異系數(shù)結(jié)果見表6。表6顯示，2003－2013年，中國(guó)海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度年度變異系數(shù)總體水平較低:除2003－2005年三年外，其余年份的變異系數(shù)值均低于變異系數(shù)均值0.175；其中2003年，變異系數(shù)值最高達(dá)0.29，2012年最低，僅為0.12。從變化趨勢(shì)上看，協(xié)調(diào)度的年度變異系數(shù)值整體趨于下降:除2005－2006年顯著降低外，變異系數(shù)值在年度間波動(dòng)較小。

表6 2003－2013中國(guó)海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度年度變異系數(shù)表

3.2.2 海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度的時(shí)序分異特征

基于表5、表6和圖1、圖2，我國(guó)沿海11省海洋漁業(yè)環(huán)境適應(yīng)協(xié)調(diào)度的時(shí)序特征可歸納為:2003－2013年，除福建省外，我國(guó)沿海省份海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度總體上都有顯著提升。我國(guó)海洋漁業(yè)協(xié)調(diào)度年度均值由0.463上升至0.656，上升了0.194，約躍遷了一個(gè)協(xié)調(diào)等級(jí)，總體從勉強(qiáng)協(xié)調(diào)發(fā)展至總體中度協(xié)調(diào)。此外，各省份海洋漁業(yè)與海洋資源環(huán)境協(xié)調(diào)發(fā)展的差異仍然存在。至2013年末，山東省先于其他10省達(dá)到優(yōu)質(zhì)協(xié)調(diào)狀態(tài)。遼寧、上海、浙江、江蘇、廣東、廣西、海南7省處于中度協(xié)調(diào)狀態(tài)，天津、河北、福建三省處于勉強(qiáng)協(xié)調(diào)狀態(tài)。從趨勢(shì)上看，各省份海洋漁業(yè)環(huán)境適應(yīng)協(xié)調(diào)度在總體上升的背景下趨同的趨勢(shì)漸為明顯，因?yàn)檠睾?1省海洋漁業(yè)環(huán)境適應(yīng)協(xié)調(diào)度的變異系數(shù)在2003－2013年間有了較大幅度的下降，自0.291降至0.171，區(qū)域間協(xié)調(diào)度的空間差異在逐年減小。

圖1 海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度

圖2 各省份海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度

3.4 海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)耦合協(xié)調(diào)度的空間分異分析

3.4.1 海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)耦合協(xié)調(diào)度的空間分異測(cè)度

通過計(jì)算各省市局部Moran’s I統(tǒng)計(jì)值，并借助OpenGeoda軟件建立空間矩陣得到2003年、2009年、2013年各省份海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度的LISA圖，并據(jù)此對(duì)海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)耦合協(xié)調(diào)度的局部空間分異進(jìn)行測(cè)度。

圖3中，High－High表示高－高，代表協(xié)調(diào)度高的區(qū)聚集在一起；Low－Low表示低－低，代表協(xié)調(diào)度低的區(qū)聚集在一起；High－Low表示高－低，代表本區(qū)海洋漁業(yè)環(huán)境適應(yīng)協(xié)調(diào)度高，但是周邊區(qū)環(huán)境適應(yīng)協(xié)調(diào)度低；Low－High表示低－高，代表本區(qū)海洋漁業(yè)環(huán)境適應(yīng)協(xié)調(diào)度低，但是周邊區(qū)環(huán)境適應(yīng)協(xié)調(diào)度高；Not Significant則代表該區(qū)協(xié)調(diào)度集聚關(guān)系不顯著。

圖3 2003年(左)、2009年(中)、2013年(右)沿海11省海洋漁業(yè)環(huán)境適應(yīng)協(xié)調(diào)度LISA集聚圖

3.4.2 海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)耦合協(xié)調(diào)度的空間分異特征

根據(jù)圖3，中國(guó)省際海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度所呈現(xiàn)的局部空間分異特征可歸結(jié)為:

(1)2003年，華北平原中的沿海一帶(即山東省、河北省、江蘇省)海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度差異顯著。如圖3所示，山東省為高－低區(qū)，其中，山東省海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度較高，而臨近的河北省和江蘇省海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度則相對(duì)較低。

(2)上海市海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)耦合度、協(xié)調(diào)度穩(wěn)步提升。特別，在2009年各省份協(xié)調(diào)度普遍下行的情況下，上海市海洋漁業(yè)環(huán)境適應(yīng)協(xié)調(diào)度仍具明顯的上升勢(shì)頭。2013年江蘇省成為“高－高區(qū)”，意味著以江蘇省為中心的山東省、江蘇省、浙江省、上海市四省市的海洋漁業(yè)環(huán)境適應(yīng)協(xié)調(diào)度整體處于高水平狀態(tài)，這一轉(zhuǎn)變也折射出上海市海洋漁業(yè)環(huán)境適應(yīng)協(xié)調(diào)度的提升勢(shì)頭。

(3)2003－2013年，長(zhǎng)江三角洲經(jīng)濟(jì)區(qū)中沿海一帶的江、浙、滬區(qū)域海洋漁業(yè)經(jīng)濟(jì)生態(tài)化發(fā)展始終存在差距。由圖3可知，2003年、2009年、2013年上海市皆為“低－高區(qū)”，即在該區(qū)域內(nèi)部，上海市海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度雖穩(wěn)步提升但與江蘇、浙江相比仍處于較低水平。

4 結(jié)論和建議

本文在構(gòu)建多指標(biāo)評(píng)價(jià)體系的基礎(chǔ)上，建立海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度時(shí)空分異模型進(jìn)行實(shí)證研究，得出海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度時(shí)空分異特征:(1)時(shí)序上，2003－2013年，中國(guó)沿海11省份海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度總體提升顯著(福建省除外)，但省際海洋漁業(yè)與海洋資源環(huán)境協(xié)調(diào)發(fā)展的差異仍然存在，且省際海洋漁業(yè)環(huán)境適應(yīng)協(xié)調(diào)度在總體上升的背景下趨同趨勢(shì)漸為明顯。(2)從局部區(qū)域上看，2003年，山東－河北－江蘇區(qū)域內(nèi)部海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度差異顯著；2003－2013上海市海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)耦合度、協(xié)調(diào)度穩(wěn)步提升，使得山東－江蘇－浙江－上海局部區(qū)域整體的海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)耦合度、協(xié)調(diào)度水平提升至較高水平；2003－2013年，長(zhǎng)江三角洲經(jīng)濟(jì)區(qū)中沿海一帶的江、浙、滬區(qū)域海洋漁業(yè)經(jīng)濟(jì)生態(tài)化發(fā)展差距始終存在。

綜上，我國(guó)海洋漁業(yè)經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度總體提升顯著，但局部區(qū)域差異仍舊存在。未來應(yīng)當(dāng)從以下幾個(gè)方面著力實(shí)現(xiàn)海洋經(jīng)濟(jì)系統(tǒng)環(huán)境適應(yīng)協(xié)調(diào)度的提升，并縮小其在區(qū)域間的差距:(1)加強(qiáng)海洋漁業(yè)布局規(guī)劃，優(yōu)化海洋漁業(yè)結(jié)構(gòu)。大力發(fā)展海洋生態(tài)漁業(yè)，助研漁業(yè)高新技術(shù)，確定合理養(yǎng)殖容量，保障生產(chǎn)污染不超過海域承載力，以期推動(dòng)海洋漁業(yè)轉(zhuǎn)型發(fā)展；(2)充分發(fā)揮資源環(huán)境系統(tǒng)的自我凈化能力。積極研發(fā)污水治理設(shè)備與方法，推進(jìn)養(yǎng)殖結(jié)構(gòu)優(yōu)化，加大藻類貝類養(yǎng)殖比例；(3)嚴(yán)控陸源污染入海量，助推排污許可證制度實(shí)行，加強(qiáng)污染管理法制建設(shè)；(4)加強(qiáng)海洋災(zāi)害應(yīng)急響應(yīng)體系，減少突發(fā)性污染事故發(fā)生，保障海洋功能實(shí)現(xiàn)。

[1]農(nóng)業(yè)部.全國(guó)漁業(yè)發(fā)展第十三個(gè)五年規(guī)劃[N].中國(guó)漁業(yè)報(bào)，2017－01－09(A01).

[2]郭付友，李誠(chéng)固，陳才，等.2003年以來東北區(qū)人口城鎮(zhèn)化與土城鎮(zhèn)化時(shí)空耦合特征[J].經(jīng)濟(jì)地理，2015，35(9):49－56.

[3]宋馬林，徐龍，王舒鴻，等.小樣本數(shù)據(jù)下安徽省經(jīng)濟(jì)－環(huán)境系統(tǒng)協(xié)調(diào)發(fā)展的實(shí)證分析[J].地理研究，2013，32(8):1402－1410.

[4]蘇靜，胡宗義，唐李偉.我國(guó)能源－經(jīng)濟(jì)－環(huán)境(3E)系統(tǒng)協(xié)調(diào)度的理空間分布與動(dòng)態(tài)演進(jìn)[J].經(jīng)濟(jì)地理，2013，33(9):19－24，30.

[5]張廣海，劉真真，王新越.中國(guó)沿海區(qū)域旅游化與生態(tài)環(huán)境耦合度分析及預(yù)測(cè)[J].生態(tài)環(huán)境學(xué)報(bào)，2013，22(5):792－800.

[6]邵桂蘭，韓菲，李晨.基于主成分分析的海洋經(jīng)濟(jì)可持續(xù)發(fā)展能力測(cè)算——以山東省2000－2008年數(shù)據(jù)為例[J].中國(guó)海洋大學(xué)學(xué)報(bào)(社會(huì)科學(xué)版)，2011(6):18－22.

[7]蒲新明，傅明珠，王宗靈，等.海水養(yǎng)殖生態(tài)系統(tǒng)健康綜合評(píng)價(jià):方法與模式[J].生態(tài)學(xué)報(bào)，2012，32(19):6210－6222.

[8]李京梅，郭斌.我國(guó)海水養(yǎng)殖的生態(tài)預(yù)警評(píng)價(jià)指標(biāo)體系與方法[J].海洋環(huán)境科學(xué)，2012，31(3):448－452.

[9]許冬蘭，李玉強(qiáng).基于狀態(tài)空間法的海洋生態(tài)環(huán)境承載力評(píng)價(jià)[J].統(tǒng)計(jì)與決策，2013(18):58－60.

[10]高樂華，高強(qiáng).海洋生態(tài)經(jīng)濟(jì)系統(tǒng)交互脅迫關(guān)系驗(yàn)證及其協(xié)調(diào)度測(cè)算[J].資源科學(xué)，2012，34(1):173－184.

[11]王秀娟，胡求光.中國(guó)海水養(yǎng)殖與海洋生態(tài)環(huán)境協(xié)調(diào)度分析[J].中國(guó)農(nóng)村經(jīng)濟(jì)，2013(11):86－96.

[12]高源，韓增林，楊俊，等.中國(guó)海洋產(chǎn)業(yè)空間集聚及其協(xié)調(diào)發(fā)展研究[J].地理科學(xué)，2015，35(8):946－951.

[13]陳琦，李京梅.我國(guó)海洋經(jīng)濟(jì)增長(zhǎng)與海洋環(huán)境壓力的脫鉤關(guān)系研究[J].海洋環(huán)境科學(xué)，2015，34(6):827－833.

[14]王洪橋，袁家冬，孟祥君.東北三省旅游經(jīng)濟(jì)差異的時(shí)空特征分析[J].地理科學(xué)，2014，34(2):163－169.

[15]肖剛.中國(guó)外商直接投資區(qū)位分布的時(shí)空格局演變[J].當(dāng)代財(cái)經(jīng)，2015(10):97－107.

[16]鄒琳，曾剛，曹賢忠.基于ESDA的長(zhǎng)三角城市群研發(fā)投入空間分異特征及時(shí)空演化[J].經(jīng)濟(jì)地理，2015，35(3):73－79.

[17]胡艷興，潘竟虎，陳蜒，等.基于ESDA和GWR的中國(guó)地級(jí)及以上城市四化協(xié)調(diào)發(fā)展時(shí)空分異格局[J].經(jīng)濟(jì)地理，2015，35 (5):45－54.

[18]王俊元.浙江海洋漁業(yè)空間布局演變及影響因素研究[J].科技與管理，2016，18(2):1－6.

[19]李昌峰，武清華，張落成.土地集約利用與經(jīng)濟(jì)發(fā)展的空間差異研究——以長(zhǎng)江三角洲區(qū)為例[J].經(jīng)濟(jì)地理，2011，31(2): 294－299.

[20]張榮天，焦華富.泛長(zhǎng)三角城市發(fā)展效率時(shí)空格局演化與驅(qū)動(dòng)機(jī)制[J].經(jīng)濟(jì)地理，2014，34(5):48－54.

Spatial－temporal Differentiation of Marine Fishery Economic System Environmental Adaptability——A Case Study of Eleven Coastal Provinces and Cities

YU Jinkai，LI Sida
(School of Economics，Ocean University of China，Qingdao 266100，China)

With the 2003－2013 panel data from eleven China coastal provinces and cities as sample and based on the evaluation index system for the coupling coordination development of marine fishery economic system environmental adaptation，this paper studies marine fishery economic system environmental adaptation coupling coordination degree from the two dimensions of time and space by adopting the corresponding spatial－temporal differentiation model.The results show that from the view of time，the environmental adaptation coordination degree of China inter－provincial marine fishery economic system increases significantly with a trend of convergence while from the view of local space，there exist regional coordination degree differences between Jiangsu，Zhejiang and Shanghai，and the steady improvement of Shanghai coordination degree drives the overall regional coordination degree among Shandong，Jiangsu，Zhejiang and Shanghai to a higher level.Therefore，this paper puts forward such countermeasures as strengthening marine fishery planning and giving a full play of the self purification ability role of resource environment system in order to achieve a high－level coordination between China inter－provincial marine fishery economic system and resource environment system and further narrow the gap between provinces.

marine fishery；economic system；coordination degree；spatial－temporal differentiation

F124.5

:A

:2095－929X(2017)03－0092－09

(責(zé)任編輯 時(shí)明芝)

2017－03－01

國(guó)家自然科學(xué)基金面上項(xiàng)目“海域承載力視角下海洋漁業(yè)空間布局優(yōu)化的模型及應(yīng)用”(71273247)；國(guó)家自然科學(xué)基金項(xiàng)目“基于空間多標(biāo)準(zhǔn)分析的我國(guó)海水養(yǎng)殖空間布局優(yōu)化決策:方法及應(yīng)用”(7167030345)；教育部人文社會(huì)科學(xué)研究項(xiàng)目“我國(guó)海洋產(chǎn)業(yè)空間布局優(yōu)化研究”(15JJDZONGHE024)。

于謹(jǐn)凱，男，山東青島人，經(jīng)濟(jì)學(xué)博士，中國(guó)海洋大學(xué)經(jīng)濟(jì)學(xué)院教授，研究方向:產(chǎn)業(yè)經(jīng)濟(jì)學(xué)、國(guó)際經(jīng)濟(jì)學(xué)；李思達(dá)，男，山西長(zhǎng)治人，中國(guó)海洋大學(xué)經(jīng)濟(jì)學(xué)院碩士生，研究方向:資源、環(huán)境、貿(mào)易戰(zhàn)略與政策，Email:lisida1994＠126.com。