基于SBM模型的我國海洋漁業(yè)生態(tài)效率的時空演變

韓增林 計雪晴 胡盈 蔡先哲

摘要：隨著“海洋強國”建設(shè)的深入實施，我國海洋經(jīng)濟實力不斷提升，但也出現(xiàn)海洋環(huán)境污染和生態(tài)破壞等制約海洋產(chǎn)業(yè)可持續(xù)發(fā)展的問題，亟須改善海洋生態(tài)環(huán)境和提高海洋產(chǎn)業(yè)生態(tài)效率。文章以海水養(yǎng)殖排污總量作為非期望產(chǎn)出，采用SBM模型、標準差和變異系數(shù)，分析2006—2015年我國海洋漁業(yè)生態(tài)效率的時空演變特征，并對我國沿海地區(qū)的海洋漁業(yè)生態(tài)效率水平進行分類。研究結(jié)果表明：2006—2015年我國海洋漁業(yè)生態(tài)效率整體處于中等水平，且呈波動式下降趨勢;沿海各地區(qū)的海洋漁業(yè)生態(tài)效率差異顯著，且差異呈不斷擴大的趨勢;沿海各地區(qū)的海洋漁業(yè)生態(tài)效率水平可分為高效型、中效型和低效型3個類型，其中廣東省和江蘇省從中效型下降至低效型，主要原因在于其存在生態(tài)恢復(fù)滯后性。

關(guān)鍵詞：海洋漁業(yè);生態(tài)效率;非期望產(chǎn)出;海洋生態(tài)環(huán)境;海洋牧場

中圖分類號：F3074;P74文獻標志碼：A文章編號：1005-9857（2019）12-0003-06

The Spatialtemporal Evolution of Marine Fishery Ecoefficiency Based on SBM Model in China

HAN Zenglin12，JI Xueqing2，HU Ying2，CAI Xianzhe3

（1Marine Economics and Sustainable Development Research Center，Liaoning Normal University，Dalian 116029，China;

2School of Geography，Liaoning Normal University，Dalian 116029，China;

3School of Maritime Economics and Management，Dalian Maritime University，Dalian 116029，China）

Abstract：With the implementation of the strategy of building a “maritime power”，the marine economic strength of China has been constantly improving.But at the same time there have been a series of problems that restrict the sustainable development of the marine primary industry，such as marine environmental pollution and ecological destruction.Therefore，improving the marine ecological environment and the ecological efficiency of marine fisheries are the basic guarantees for achieving sustainable development of the marine economy.Based on the total amount of sewage discharged from mariculture process as an unexpected output，the temporal and spatial evolution characteristics of China′s marine fishery ecological efficiency from 2006 to 2015 were analyzed by using SBM model，standard deviation and coefficient of variation.And classifying the level of ecological efficiency in coastal area.The results showed that during the study period，the average value of marine fishery ecoefficiency in China，which had taken into account the undesired output，was at a moderate level and fluctuating downward，while the regional differences showed a trend of fluctuation and expansion.The marine fishery ecoefficiency level in China was divided into three categories：highefficiency，mediumefficiency and lowefficiency.During the study period，Guangdong and Jiangsu dropped from moderateefficiency to lowefficiency.The main reason lied in the lag of ecological restoration in Guangdong and Jiangsu provinces.

Key words：Marine fisheries，Ecological efficiency，Unexpected output，Marine ecological environment，Marine ranch

0引言

藍色經(jīng)濟作為新時代區(qū)域發(fā)展的動力源泉之一，對提升綜合國力發(fā)揮重要作用[1]。黨的十九大提出“堅持陸海統(tǒng)籌，加快建設(shè)海洋強國”，將海洋強國建設(shè)納入現(xiàn)代化經(jīng)濟體系中，踐行綠色發(fā)展新理念和提高海洋生態(tài)效率是海洋經(jīng)濟可持續(xù)發(fā)展的關(guān)鍵節(jié)點。現(xiàn)代海洋漁業(yè)是世界四大海洋支柱產(chǎn)業(yè)之一，其產(chǎn)值對國家海洋經(jīng)濟實力的提升作用日益顯現(xiàn)。然而隨著沿海地區(qū)經(jīng)濟發(fā)展水平的快速提高以及對海洋自然資源開發(fā)程度的不斷加深，傳統(tǒng)粗放型的海洋經(jīng)濟發(fā)展模式已不符合節(jié)約集約和循環(huán)利用的資源觀，且隨著海洋資源的過量消耗，出現(xiàn)漁業(yè)資源枯竭和海洋生態(tài)環(huán)境惡化等問題[2]，嚴重制約海洋漁業(yè)的可持續(xù)發(fā)展，已成為沿海地區(qū)乃至全世界亟待解決的問題。

國內(nèi)外對于海洋漁業(yè)經(jīng)濟效率的評價研究較多，如運用CD生產(chǎn)函數(shù)法[3]、DEA模型[4-5]和Malmquist指數(shù)[6]等方法對漁業(yè)經(jīng)濟效率進行測算，然而對海洋漁業(yè)生態(tài)效率的評價研究較少。國外海洋生態(tài)效率研究多集中于海港生態(tài)效率方面[7]，國內(nèi)海洋生態(tài)效率研究則多集中于海洋生態(tài)經(jīng)濟可持續(xù)發(fā)展[8]、海洋產(chǎn)業(yè)[9]和海洋環(huán)境效率[10]等方面，而鮮有對海洋漁業(yè)生態(tài)效率的測算和評價。本研究借鑒現(xiàn)有理論基礎(chǔ)和研究成果，嘗試將生態(tài)效率中的經(jīng)濟發(fā)展和生態(tài)系統(tǒng)聚焦于海洋漁業(yè)經(jīng)濟和海洋漁業(yè)生態(tài)系統(tǒng)，提出“海洋漁業(yè)生態(tài)效率”的概念，體現(xiàn)海洋漁業(yè)經(jīng)濟效益與海洋漁業(yè)資源環(huán)境消耗的定量關(guān)系;運用超效率SBM模型，以海水養(yǎng)殖排污總量作為非期望產(chǎn)出，對2006—2015年我國沿海11?。ㄗ灾螀^(qū)、直轄市）海洋漁業(yè)生態(tài)效率的時空演變特征進行分析，以期為我國提升海洋漁業(yè)生態(tài)效率水平和實現(xiàn)海洋漁業(yè)經(jīng)濟可持續(xù)發(fā)展提供具有可行性的參考。

1研究方法和指標數(shù)據(jù)

11研究方法

在測算決策單元的效率方面，Tone[11]提出的考慮非期望產(chǎn)出的SBM（SlackBased Measure）模型直接在目標函數(shù)中加入松弛變量，既解決投入產(chǎn)出中的松弛性問題，又解決存在非期望產(chǎn)出時的效率值評價問題[12]。因此，本研究選用基于非期望產(chǎn)出的超效率SBM模型，對我國海洋漁業(yè)生態(tài)效率進行測算。參考相關(guān)文獻[13-14]，設(shè)ρ*為海洋漁業(yè)生態(tài)效率值，設(shè)定06<ρ*≤1為高等效率水平、04<ρ*≤06為中等效率水平、0<ρ*≤04為低等效率水平。

12指標選取和數(shù)據(jù)來源

運用SBM模型計算效率的準確性多依賴于在測算過程中選用的投入和產(chǎn)出指標[15-17]。本研究測算我國第一產(chǎn)業(yè)中海洋漁業(yè)的生態(tài)效率值，選取沿海11?。ㄗ灾螀^(qū)、直轄市）的指標數(shù)據(jù)。在海洋漁業(yè)的投入指標中，參考已有研究成果[18-19]，選用海水養(yǎng)殖面積、國家級水產(chǎn)種質(zhì)資源保護區(qū)（僅包括海水資源保護區(qū)）面積和海水魚苗數(shù)量3項指標，同時將海洋捕撈漁船功率等資本和勞動力投入也納入投入指標。

非期望產(chǎn)出是伴隨期望產(chǎn)出而產(chǎn)生的，是對總體目標即生態(tài)環(huán)境無益和不符合期待的產(chǎn)出。根據(jù)《水產(chǎn)養(yǎng)殖業(yè)污染源產(chǎn)排污系數(shù)手冊》[20]，選取各地區(qū)海水養(yǎng)殖產(chǎn)量占比約為80%的主要海產(chǎn)品，依據(jù)海水養(yǎng)殖業(yè)排污系數(shù)測算其養(yǎng)殖過程中氮、磷和化學需氧量3種污染物的排放量，以此作為非期望產(chǎn)出，測算生態(tài)效率。其中，貝類養(yǎng)殖所產(chǎn)生的氮、磷排污系數(shù)為負值，這歸因于貝類特殊的生理機能[21]，故在計算貝類養(yǎng)殖排污總量時僅測算化學需氧量的排放量。

綜合上述指標的選取，建立我國海洋漁業(yè)生態(tài)效率測算指標體系（表1）。

由于SBM模型一般要求決策單元的數(shù)量不少于投入和產(chǎn)出指標數(shù)量的乘積，同時不少于投入和產(chǎn)出指標數(shù)量的3倍，本研究采用熵值法對各項二級指標進行加權(quán)求和，再用熵值法將投入的一級指標融為投入指標，最終匯總為投入、期望產(chǎn)出和非期望產(chǎn)出3項指標數(shù)據(jù)。所需數(shù)據(jù)來自2007—2016年《中國漁業(yè)統(tǒng)計年鑒》《中國海洋統(tǒng)計年鑒》和我國農(nóng)業(yè)農(nóng)村部發(fā)布的相關(guān)文件[22]。

2我國海洋漁業(yè)生態(tài)效率演變

21海洋漁業(yè)生態(tài)效率差異

211時間差異

利用DEA軟件，計算基于非期望產(chǎn)出的2006—2015年我國沿海地區(qū)海洋漁業(yè)生態(tài)效率的平均值（圖1）。

由圖1可以看出，研究期內(nèi)我國海洋漁業(yè)生態(tài)效率的平均值在048～058之間波動，處于中等效率水平，且呈現(xiàn)不穩(wěn)定下降趨勢。其中，2006—2007年呈上升趨勢，但在2008年降至最低值（048），2009—2012年又呈波動上升趨勢，并于2012年達到最高值（058），之后下降至前一階段水平。

研究期內(nèi)我國海洋漁業(yè)生態(tài)效率平均值2008年為最低值，2012年為最高值。2008年我國受金融危機滯后性的影響，在相關(guān)國家政策的支持下通過拉動內(nèi)需恢復(fù)經(jīng)濟增長活力，但隨之而來的資源環(huán)境問題逐漸顯現(xiàn)，工業(yè)廢水超標排放和資源過度開發(fā)早已超過海洋生態(tài)資源環(huán)境承載力，傳統(tǒng)粗放式的經(jīng)濟發(fā)展模式制約海洋漁業(yè)的可持續(xù)發(fā)展;此外，2008年夏季登陸我國大陸的臺風偏多，且登陸期集中、影響范圍廣，帶來的強降雨導致海岸帶資源破壞以及沿海地區(qū)漁業(yè)生產(chǎn)損失嚴重。

2012年山東省和浙江省的海洋漁業(yè)生態(tài)效率值均為1，即達到完全有效狀態(tài)。2012年山東省和浙江省的海洋漁業(yè)產(chǎn)值分別為2011年的17倍和15倍;2012年山東省增設(shè)無棣中國毛蝦等4處國家級海產(chǎn)種質(zhì)資源保護區(qū)，浙江省增設(shè)1處國家級海產(chǎn)種質(zhì)資源保護區(qū)，使?jié)O業(yè)資源枯竭的趨勢得到緩解，同時促進海洋漁業(yè)生態(tài)環(huán)境的改善。

212空間差異

依據(jù)計算出的我國海洋漁業(yè)生態(tài)效率值，利用標準差和變異系數(shù)計算公式，定量分析我國各沿海地區(qū)海洋漁業(yè)生態(tài)效率的絕對差異和相對差異（圖2）。

由圖2可以看出：我國海洋漁業(yè)生態(tài)效率區(qū)域間的絕對差異數(shù)值較小，而相對差異數(shù)值較大，二者變化趨勢相同，且均呈波動增長態(tài)勢，表明區(qū)域發(fā)展不平衡，差異日益加劇;其中，2008年差異最小，2012年差異最大，2012年后差異變化趨勢較穩(wěn)定，但2015年又呈擴大趨勢。

我國海洋漁業(yè)生態(tài)效率的區(qū)域差異不平穩(wěn)，且呈擴大趨勢。究其原因，我國大部分沿海地區(qū)受夏季臺風登陸影響，海洋漁業(yè)受災(zāi)嚴重，其中天津市、福建省和遼寧省因漁業(yè)災(zāi)情等原因遭受的經(jīng)濟損失巨大，環(huán)境破壞嚴重，從而形成與未受災(zāi)地區(qū)的效率差異。

22海洋漁業(yè)生態(tài)效率類型

根據(jù)我國沿海地區(qū)海洋漁業(yè)生態(tài)效率值以及“十一五”和“十二五”2個時期的平均值（表2），參考已有研究成果[13-14]，將沿海11省（自治區(qū)、直轄市）的效率水平劃分為高效型、中效型和低效型3個類型（表3）。

由表2和表3可以看出：①“十一五”時期，浙江、山東和海南的海洋漁業(yè)生態(tài)效率平均值均大于06，屬于高效型;福建、廣西、遼寧、廣東、江蘇和河北的海洋漁業(yè)生態(tài)效率平均值均大于04但小于06，屬于中效型;上海和天津的海洋漁業(yè)生態(tài)效率平均值均小于04，屬于低效型。②“十二五”時期，浙江、山東和海南的海洋漁業(yè)生態(tài)效率平均值均大于06，屬于高效型;福建、遼寧、廣西和河北的海洋漁業(yè)生態(tài)效率平均值均大于04但小于06，屬于中效型;天津、江蘇、廣東和上海的海洋漁業(yè)生態(tài)效率平均值均小于04，屬于低效型。③對比“十一五”和“十二五”2個時期，廣東和江蘇從中效型下降至低效型。

（1）高效型地區(qū)包括浙江、山東和海南。浙江是我國傳統(tǒng)漁業(yè)大省，在當?shù)貎?yōu)勢資源的基礎(chǔ)上，“十一五”期間深入推進生態(tài)省建設(shè)，扎實開展“811”環(huán)境污染整治行動，認真實施資源節(jié)約與環(huán)境保護行動計劃，且于2007年建立多個國家級水產(chǎn)種質(zhì)資源保護區(qū);山東也是我國傳統(tǒng)漁業(yè)大省，自2005年起實施《山東省漁業(yè)資源修復(fù)規(guī)劃》，在全省沿海大范圍開展海洋牧場和人工魚礁建設(shè)，并于2014年印發(fā)并實行《關(guān)于推進海上糧倉建設(shè)的實施意見》，獲得良好的經(jīng)濟效益和生態(tài)效益;海南具有得天獨厚的地理位置，“十二五”期間實施“綠化寶島”大行動，積極推進省域“多規(guī)合一”工作，劃定海南島生態(tài)保護紅線，以期改善海洋生態(tài)環(huán)境，雖然局部環(huán)境污染依然存在，但整體處于高效水平。

（2）中效型地區(qū)包括福建、廣西、遼寧和河北。福建海洋經(jīng)濟基礎(chǔ)雄厚，提倡實施資源開發(fā)與環(huán)境保護并行的理念，但污染物排放強度處于較高水平，海洋生態(tài)環(huán)境的修復(fù)稍滯后于經(jīng)濟的總體發(fā)展;廣西海洋經(jīng)濟以第一產(chǎn)業(yè)為主，且產(chǎn)業(yè)結(jié)構(gòu)以資源型為主，排污強度處于較高水平，“十二五”時期確立“生態(tài)立區(qū)，綠色崛起”戰(zhàn)略，積極推進生態(tài)文明建設(shè)，以期改善海洋生態(tài)環(huán)境惡化趨勢;遼寧具有渤海灣優(yōu)越的地理區(qū)位，但重工業(yè)帶動經(jīng)濟高速發(fā)展也對資源環(huán)境產(chǎn)生巨大壓力，在調(diào)整工業(yè)生產(chǎn)方式和推進污染整治工作的同時，“十二五”后期積極推進海洋牧場示范區(qū)建設(shè);河北處于渤海灣西部，長期的工業(yè)生產(chǎn)活動加之資源稟賦不足，制約海洋生態(tài)環(huán)境的修復(fù)和保護，但其積極推進生態(tài)補償機制，使受損海洋資源和海洋生態(tài)功能得到初步恢復(fù)。

（3）低效型地區(qū)包括上海和天津，江蘇和廣東由中效型降為低效型。作為直轄市，上海和天津一直是我國重要的沿海樞紐港和綜合運輸體系中心，第二和第三產(chǎn)業(yè)發(fā)展突出，而對漁業(yè)等第一產(chǎn)業(yè)的依賴性較低;同時，由于大型船舶運輸頻繁、工業(yè)廢水排放和自然災(zāi)害等原因，相關(guān)海域受到不同程度的污染，導致海洋漁業(yè)生態(tài)效率處于低效水平。廣東依靠捕撈業(yè)和養(yǎng)殖業(yè)快速發(fā)展海洋漁業(yè)經(jīng)濟，同時海洋生態(tài)環(huán)境遭到相當程度的破壞，隨著“河長制”的全面實行和治污行動的嚴格落實，海洋漁業(yè)資源環(huán)境有所修復(fù)和改善，但前期污染較嚴重，且漁業(yè)經(jīng)濟發(fā)展的速度快于資源環(huán)境治理的速度，受生態(tài)修復(fù)滯后性的影響，其生態(tài)效率水平呈下降趨勢。江蘇沿海漁場面積廣大，但漁業(yè)養(yǎng)殖水域污染嚴重，海洋生態(tài)平衡遭到破壞，盡管已積極開展?jié)O業(yè)資源增殖放流與人工魚礁建設(shè)，并加強陸源污染物排海監(jiān)管，但其近海海域長期處于亞健康狀態(tài)，恢復(fù)周期較長，導致生態(tài)效率有所下降。

3結(jié)語

本研究結(jié)合地區(qū)發(fā)展實情，對我國海洋漁業(yè)生態(tài)效率的時空演變進行分析，得到主要結(jié)論：①在考慮非期望產(chǎn)出的背景下，我國海洋漁業(yè)生態(tài)效率整體處于中等水平且波動式下降，結(jié)合效率值的標準差和變異系數(shù)，海洋漁業(yè)生態(tài)效率的區(qū)域差異呈現(xiàn)波動增大趨勢;②沿海地區(qū)海洋漁業(yè)生態(tài)效率可分為高效型、中效型和低效型3個類型，其中浙江、山東和海南屬于高效型，福建、廣西、遼寧和河北屬于中效型，上海和天津?qū)儆诘托?，江蘇和廣東由中效型降為低效型。

基于上述結(jié)論，建議：高效型地區(qū)應(yīng)與周邊地區(qū)形成聯(lián)動機制，產(chǎn)生溢出和帶動作用，同時提高漁業(yè)技術(shù)水平，加強漁業(yè)科技創(chuàng)新，進一步推廣和轉(zhuǎn)變漁業(yè)養(yǎng)殖方式，積極發(fā)展種養(yǎng)結(jié)合稻田養(yǎng)殖和海洋牧場立體養(yǎng)殖等健康養(yǎng)殖模式;中效型地區(qū)應(yīng)進一步實施生態(tài)補償機制，同時由政府出資引導和鼓勵養(yǎng)殖節(jié)水減排改造，重點投入支持循環(huán)用水等環(huán)保設(shè)施的升級改造，對養(yǎng)殖水質(zhì)進行監(jiān)測，制定廢水排放標準并嚴格落實，調(diào)整并相應(yīng)減少近海過于集中的網(wǎng)箱養(yǎng)殖，根據(jù)養(yǎng)殖品種和養(yǎng)殖容量科學布局，加強養(yǎng)殖技術(shù)創(chuàng)新;低效型地區(qū)應(yīng)學習先進地區(qū)的管理經(jīng)驗和治理模式，根據(jù)漁業(yè)資源稟賦和生態(tài)環(huán)境現(xiàn)狀，在治理原有污染海域的同時，科學有序地利用漁業(yè)資源，加大治理環(huán)節(jié)的財政投入，保護灘涂生態(tài)環(huán)境，實施近海養(yǎng)殖水域生態(tài)環(huán)境修復(fù)工程，向深遠海拓展養(yǎng)殖空間，形成水域和灘涂資源綜合利用與保護的新格局。

為促進我國海洋漁業(yè)經(jīng)濟的可持續(xù)發(fā)展和提升海洋漁業(yè)生態(tài)效率，沿海地區(qū)應(yīng)在實施海洋資源環(huán)境保護的同時，綜合考量地區(qū)發(fā)展實情，培育新型職業(yè)漁民，增強漁民的創(chuàng)業(yè)能力和就業(yè)技能，并大力開展區(qū)域協(xié)作和經(jīng)驗交流，促進形成互補機制，加強漁業(yè)產(chǎn)業(yè)分工與高新技術(shù)等的良性關(guān)系，從而實現(xiàn)海洋漁業(yè)資源最佳生態(tài)效益與經(jīng)濟效益的良性循環(huán)，進一步推進我國海洋強國建設(shè)。

參考文獻

[1]王澤宇，郭萌雨，韓增林.基于集對分析的海洋綜合實力評價研究[J].資源科學，2014，36（2）：351-360.

[2]狄乾斌，韓增林.遼寧省海洋經(jīng)濟可持續(xù)發(fā)展的演進特征及其系統(tǒng)耦合模式[J].經(jīng)濟地理，2009，29（5）：799-805.

[3]馮浩，車斌.我國海洋漁業(yè)經(jīng)濟增長方式轉(zhuǎn)變及其影響因素研究[J].中國漁業(yè)經(jīng)濟，2017，35（6）：74-79.

[4]孫康，季建文，李麗丹，等.基于非期望產(chǎn)出的中國海洋漁業(yè)經(jīng)濟效率評價與時空分異[J].資源科學，2017，39（11）：2040-2051.

[5]梁盼盼，俞立平.中國漁業(yè)經(jīng)濟投入產(chǎn)出績效分析：基于1999-2010年面板數(shù)據(jù)的實證[J].科技與管理，2014，16（2）：21-26，43.

[6]于淑華，于會娟.中國沿海地區(qū)漁業(yè)產(chǎn)業(yè)效率實證研究：基于DEA的Malmquist指數(shù)分析[J].中國漁業(yè)經(jīng)濟，2012，30（3）：140-146.

[7]TICHAVSKA M，TOVAR B.Environmental cost and ecoefficiency from vessel emissions in Las Palmas Port[J].Transportation Research Part E，2015，83（3）：126-140.

[8]韓增林，胡偉，鐘敬秋，等.基于能值分析的中國海洋生態(tài)經(jīng)濟可持續(xù)發(fā)展評價[J].生態(tài)學報，2017，37（8）：2563-2574.

[9]狄乾斌，梁倩穎.中國海洋生態(tài)效率時空分異及其與海洋產(chǎn)業(yè)結(jié)構(gòu)響應(yīng)關(guān)系識別[J].地理科學，2018，38（10）：1606-1615.

[10]鄭金花.中國沿海十一省市海洋環(huán)境效率測度及其時空特征分析[D].大連：遼寧師范大學，2018.

[11]TONE K.A slacksbased measure of efficiency in data envelopment analysis[J].European Journal of Operational Research，2001，130（3）：498-509.

[12]TONE K.Dealing with undesirable outputs in DEA：a slacksbased measure（SBM）approach[R].GRIPS Research Report Series，2003.

[13]盧麗文，宋德勇，李小帆.長江經(jīng)濟帶城市發(fā)展綠色效率研究[J].中國人口·資源與環(huán)境，2016，26（6）：35-42.

[14]鄭德鳳，郝帥，孫才志，等.中國大陸生態(tài)效率時空演化分析及其趨勢預(yù)測[J].地理研究，2018，37（5）：1034-1046.

[15]張立新，朱道林，杜挺，等.基于DEA模型的城市建設(shè)用地利用效率時空格局演變及驅(qū)動因素[J].資源科學，2017，39（3）：418-429.

[16]彭紅松，章錦河，韓婭，等.旅游地生態(tài)效率測度的SBMDEA模型及實證分析[J].生態(tài)學報，2017，37（2）：628-638.

[17]潘丹，應(yīng)瑞瑤.中國農(nóng)業(yè)生態(tài)效率評價方法與實證：基于非期望產(chǎn)出的SBM模型分析[J].生態(tài)學報，2013，33（12）：3837-3845.

[18]官文江，田思泉，朱江峰，等.漁業(yè)資源評估模型的研究現(xiàn)狀與展望[J].中國水產(chǎn)科學，2013，20（5）：1112-1120.

[19]陳作志，林昭進，邱永松.基于AHP的南海海域漁業(yè)資源可持續(xù)利用評價[J].自然資源學報，2010，25（2）：249-257.

[20]中國水產(chǎn)科學研究院.第一次全國污染源普查水產(chǎn)養(yǎng)殖業(yè)污染源產(chǎn)排污系數(shù)手冊[R].北京：中國水產(chǎn)科學研究院，2010.

[21]陳一波，宋國寶，趙文星，等.中國海水養(yǎng)殖污染負荷估算[J].海洋環(huán)境科學，2016，35（1）：1-6，12.

[22]中華人民共和國農(nóng)業(yè)農(nóng)村部.農(nóng)業(yè)部公布國家級水產(chǎn)種質(zhì)資源保護區(qū)名單[EBOL].（2012-11-29）[2019-05-01].http：www.moa.gov.cnxwbmdt201211t20121129_3076480.htm.