于謹凱，李思達

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

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

于謹凱，李思達

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

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

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

0 引 言

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

1 文獻綜述

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

其中，D為耦合協(xié)調(diào)度，T為海洋漁業(yè)經(jīng)濟與海洋資源環(huán)境協(xié)調(diào)發(fā)展的綜合評價指數(shù)；a、b為待定系數(shù)，考慮到兩子系統(tǒng)的綜合發(fā)展具有相同的重要性，取a＝b＝0.5；F1、F2分別為海洋漁業(yè)經(jīng)濟子系統(tǒng)與海洋資源系統(tǒng)的綜合評價值。借鑒前人的相關(guān)研究，本文對海洋漁業(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)濟系統(tǒng)環(huán)境耦合協(xié)調(diào)度的時空分異測度

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

2.3.1 變異系數(shù)

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

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

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

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

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

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

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

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

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

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

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

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

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

各主成分評價綜合評價公式具體如下所示:

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

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

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

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

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

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

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

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

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

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

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

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

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

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

圖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)濟系統(tǒng)環(huán)境適應(yīng)耦合協(xié)調(diào)度的空間分異特征

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

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

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

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

4 結(jié)論和建議

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

[20]張榮天，焦華富.泛長三角城市發(fā)展效率時空格局演化與驅(qū)動機制[J].經(jīng)濟地理，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é)任編輯 時明芝)

2017－03－01

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

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