陳思露，劉百橋，胡蓓蓓

（天津師范大學(xué)地理與環(huán)境科學(xué)學(xué)院，天津300387）

近年來，構(gòu)建與資源環(huán)境承載力相適應(yīng)的社會經(jīng)濟(jì)發(fā)展模式逐漸成為社會各界的共識，生態(tài)用海理念已深入人心[1-4].落實到海域管理領(lǐng)域，這一共識可以理解為必須根據(jù)海洋資源環(huán)境承載力，合理配置海域資源，統(tǒng)籌協(xié)調(diào)海洋保護(hù)和開發(fā)活動[5-7].為此，我們要深入認(rèn)識海域的自然稟賦，及時、全面地掌握海域開發(fā)現(xiàn)狀，在此基礎(chǔ)上科學(xué)測評海洋資源環(huán)境承載力與開發(fā)負(fù)載之間的協(xié)調(diào)程度.海域使用是指持續(xù)使用特定海域3個月以上的排他性用海活動[8].目前，國內(nèi)學(xué)者對于海域使用強(qiáng)度的概念并無統(tǒng)一的內(nèi)涵界定，評估方法也多種多樣，較多的研究工作基于海岸線位置和形態(tài)等的變化程度衡量海岸開發(fā)強(qiáng)度，一方面利用填海面積與占用岸線長度之比作為圍填海開發(fā)強(qiáng)度的度量指標(biāo)[9-12]，另一方面利用某一時間段內(nèi)海岸線長度變化率、岸線曲折度和分形維數(shù)等反映海岸開發(fā)強(qiáng)度[13-15].還有一些研究工作根據(jù)不同岸線類型和景觀類型所對應(yīng)的人類干擾程度差異評估海岸區(qū)域內(nèi)的人類活動影響強(qiáng)度[16-19]，針對海灣內(nèi)的各岸線類型和景觀類型，分別進(jìn)行指數(shù)賦值，計算海岸人工化強(qiáng)度和海灣人力干擾強(qiáng)度.此外，張君玨等[20]根據(jù)海岸帶景觀的不均勻性和破碎程度反映人類活動影響強(qiáng)度，利用網(wǎng)格化的地表不透水面數(shù)據(jù)，采用泰爾指數(shù)表征沿海土地開發(fā)利用空間差異，并據(jù)此分析海岸帶開發(fā)強(qiáng)度.

現(xiàn)有研究從不同角度反映了海岸帶開發(fā)的強(qiáng)度特征，但與海域資源承載力監(jiān)測預(yù)警和海域使用管控業(yè)務(wù)等的現(xiàn)實需求還存在一定差距，主要體現(xiàn)在以下幾個方面：一是所考慮的用?；顒又饕獮樘詈Ｔ斓?，而未考慮圍海、構(gòu)筑物和開放式等方式的用?；顒?，這些未被考慮的用?；顒釉诂F(xiàn)實中客觀存在，占用較廣的海域面積，對海域環(huán)境均造成一定程度的壓力；二是未區(qū)分用海產(chǎn)業(yè)，不能反映不同產(chǎn)業(yè)用海對海域空間資源利用程度的差異性；三是不能很好地利用海域管理業(yè)務(wù)數(shù)據(jù)，需要以專門的研究數(shù)據(jù)為支撐，研究工作限于局部地區(qū)，難以適應(yīng)全國或海區(qū)范圍內(nèi)常態(tài)化快速評估的需要.因此，本研究以海域使用動態(tài)監(jiān)測業(yè)務(wù)數(shù)據(jù)為信息源，設(shè)計了一套根據(jù)海域使用范圍、類型和方式評估海域使用負(fù)載的空間表征方法，該方法以海域使用負(fù)載的空間密度表征海域使用強(qiáng)度的思路和方法，為大尺度范圍內(nèi)海域使用強(qiáng)度的快速高效評估提供方法參考，為海域資源承載力監(jiān)測預(yù)警、海域使用管控業(yè)務(wù)等工作提供支持.

1 基本思路

本研究旨在為海域資源承載力監(jiān)測預(yù)警和海域使用管控等工作提供可供業(yè)務(wù)運行的海域使用強(qiáng)度評估方法，方法的便捷性和結(jié)果的有效性是主要目標(biāo).海域使用強(qiáng)度空間表征方法立足于3個基本要求：①評估指標(biāo)在內(nèi)涵上應(yīng)體現(xiàn)人們對區(qū)域內(nèi)海域使用強(qiáng)度高低的一般認(rèn)識.由于用海活動對海洋資源環(huán)境具有一定的壓力，可視為海域所承受的一種負(fù)擔(dān)，本研究稱之為海域使用負(fù)載.對某一海域而言，用海項目布局越緊湊，對空間資源的利用越充分，對海域自然屬性的改變程度越大，意味著其承載的海域使用負(fù)載就越大，從承載對象的角度看，可理解為海域使用強(qiáng)度越高.②評估方法應(yīng)直接利用海域動態(tài)監(jiān)測業(yè)務(wù)數(shù)據(jù)，在較短時間內(nèi)實現(xiàn)全海域的使用強(qiáng)度評估；③評估結(jié)果應(yīng)有效反映海域使用強(qiáng)度的空間分異特征.

本研究以海域使用動態(tài)監(jiān)測業(yè)務(wù)數(shù)據(jù)為數(shù)據(jù)源，根據(jù)海域使用范圍、類型和方式評估海域使用負(fù)載，以海域使用負(fù)載的空間密度表征海域使用強(qiáng)度，并歸納其空間分布特征.主要思路為：

（1）根據(jù)某點的海域使用類型和用海方式表征該點的海域使用負(fù)載.海域使用管理包括海域使用類型的管理和用海方式的管理.不同海域使用類型在利用海域資源要素、發(fā)揮資源價值的綜合程度和充分程度方面存在差異，即海域使用類型一定程度上對應(yīng)著用海活動對海域資源的利用程度，反映其對海域資源的索取壓力；用海方式對應(yīng)著用?；顒訉Ｓ蜃匀粚傩缘挠绊懗潭龋从称鋵Ｓ蛸Y源的損耗壓力.因此，可以按照一定的規(guī)則，將各海域使用類型和用海方式分別折合為海域資源利用程度系數(shù)和海域自然屬性影響程度指標(biāo).針對具體點位，可根據(jù)該點實際的海域使用類型和用海方式，將對應(yīng)系數(shù)與指標(biāo)的乘積作為該點海域使用負(fù)載的指標(biāo)值.

（2）用某點一定鄰域內(nèi)海域使用負(fù)載指標(biāo)值的空間密度值表征該點的海域使用強(qiáng)度.如前所述，海域使用強(qiáng)度是對海域使用的區(qū)域聚集程度、資源利用程度和自然屬性影響程度的綜合反映，是衡量一定范圍內(nèi)海域使用空間分布態(tài)勢的指標(biāo).事實上，海洋開發(fā)活動既在用海范圍內(nèi)對海域空間構(gòu)成直接的壓力，也對鄰近一定范圍的海域產(chǎn)生輻射影響，某一點的海域使用強(qiáng)度既來自該點海域使用的直接作用，也來自一定鄰域內(nèi)其他點位的海域使用影響，對其評估需要考慮鄰域內(nèi)海域使用的總體情況.本研究根據(jù)前述對海域使用強(qiáng)度的描述，選用地理信息核密度分析方法，用計算所得海域使用負(fù)載的核密度值表征海域使用強(qiáng)度.

（3）通過空間聚類分析揭示海域使用強(qiáng)度分布的區(qū)域分異特征.為了更好地服務(wù)管理業(yè)務(wù)，及時提供區(qū)域分類和指示信息，尚需進(jìn)一步歸納海域使用強(qiáng)度空間分布特征.可按一定規(guī)則，對海域使用強(qiáng)度評估結(jié)果進(jìn)行空間聚類，提取海域使用高強(qiáng)度區(qū)域.

2 方法介紹

2.1 海域柵格化處理

在中國，對海域使用現(xiàn)狀信息實施“一張圖”管理[21].在海域使用動態(tài)監(jiān)測業(yè)務(wù)系統(tǒng)中，所有已知項目用海的空間分布和屬性信息儲存在一個多邊形矢量數(shù)據(jù)圖層中，其中每個項目的用海范圍記錄為一組多邊形，每個多邊形又有相應(yīng)的海域使用類型和用海方式屬性.海域使用強(qiáng)度評估需要利用地理信息系統(tǒng)的柵格圖層運算工具.評估前需要按照一定的柵格大小，如50 m×50 m，對評估海域進(jìn)行柵格化處理，并根據(jù)使用現(xiàn)狀矢量數(shù)據(jù)圖層，生成反映海域使用現(xiàn)狀信息的柵格圖層，其中主體屬于海域使用范圍內(nèi)的柵格會有相應(yīng)的海域使用類型和用海方式屬性，其他柵格則無此類屬性.

為滿足海域使用強(qiáng)度評估工作需要，運算過程中還需建立若干個柵格數(shù)據(jù)圖層，分別記載海域自然屬性影響程度指標(biāo)、海域資源利用程度系數(shù)、海域使用負(fù)載指標(biāo)和海域使用強(qiáng)度指標(biāo).

2.2 海域使用負(fù)載評估

2.2.1 海域自然屬性影響程度指標(biāo)賦值

根據(jù)各柵格對應(yīng)的用海方式，對記載海域自然屬性影響程度指標(biāo)的柵格圖層數(shù)據(jù)進(jìn)行賦值.設(shè)柵格（i，j）上的用海方式為αij，則海域自然屬性影響程度指標(biāo) Eij賦值為 e（αij）；如柵格（i，j）上無用?；顒樱瑒tEij賦值為0.本研究按照《海域使用分類》（HY/T123-2009）取填海造地、構(gòu)筑物、圍海和開放式4種用海方式，其他用海方式根據(jù)對海域自然屬性的影響程度歸入上述4種方式.出于實驗?zāi)康?，?jīng)咨詢專家，e（αij）取值為

2.2.2 海域資源利用程度系數(shù)賦值

根據(jù)各柵格對應(yīng)的海域使用類型，對記載海域資源利用程度系數(shù)的柵格圖層數(shù)據(jù)進(jìn)行賦值.設(shè)柵格（i，j）上的海域使用類型為βij，則海域資源利用程度系數(shù)Wij賦值為 w（βij）；如柵格（i，j）上無用?；顒?，則 Wij賦值為0.鑒于《海域使用分類》（HY/T123-2009）中海域使用類型較多，本研究根據(jù)空間資源利用程度相近原則，將海域使用類型歸為漁鹽業(yè)用海、工業(yè)城鎮(zhèn)用海、交通運輸用海、旅游娛樂用海和其他用海共5大類.根據(jù)專家打分結(jié)果，w（βij）取值為

2.2.3 海域使用負(fù)載指標(biāo)值計算

通過圖層運算，對記載海域使用負(fù)載指標(biāo)的柵格圖層數(shù)據(jù)進(jìn)行賦值，柵格（i，j）的海域使用負(fù)載指標(biāo)取值為

2.3 海域使用強(qiáng)度評估

將柵格圖層數(shù)據(jù)轉(zhuǎn)換為以柵格中心為點坐標(biāo)的點層數(shù)據(jù)，進(jìn)行核密度分析，計算海域使用負(fù)載指標(biāo)Pij的核密度值，以此反映海域使用強(qiáng)度的相對水平.

設(shè)柵格（i，j）的中心位置為xij，則柵格的核密度值為

式（4）中：k為核函數(shù)，本研究取Silverman[22]中的二次核函數(shù)；Pt為第t個觀測位置上的海域使用負(fù)載指標(biāo)值；dt為第t個觀測位置距xij位置的距離；h為衰減距離閾值，距離xij越近的點的海域使用負(fù)載對xij處的核密度值貢獻(xiàn)越大，當(dāng)某點與xij的路徑距離達(dá)到或超過h時，則該點對xij處的核密度值貢獻(xiàn)衰減為0；T為與xij的距離小于或等于h且具有海域使用負(fù)載指標(biāo)值的柵格數(shù).

目前，衰減距離閾值h的取值暫無理想的參考標(biāo)準(zhǔn)，仍需參照人們感性上衡量海域開發(fā)利用強(qiáng)度時著眼的鄰域空間大小，一般可按用?；顒訉Ｓ蚩臻g資源的最大影響距離取值.現(xiàn)實中，海域使用的布局比土地利用寬松很多.海域使用中的“相鄰”和“集聚”概念不能機(jī)械理解為海域使用一定要緊鄰、密布.對于海域使用間的空間關(guān)系，通常人們心目中有一個大體的空間尺度標(biāo)準(zhǔn)L0，當(dāng)海域使用間的距離大于L0時，認(rèn)為它們分離，小于L0時，認(rèn)為它們相鄰，且隨著海域使用間的距離不斷縮小，人們會認(rèn)為它們具有越來越明顯的鄰近、集聚特征.本研究定義的某點海域使用強(qiáng)度是人們對該點位周圍L0尺度范圍內(nèi)海域使用劇烈程度的一種感受.因此，h的取值應(yīng)與上述理解中的尺度L0相當(dāng).《海域使用論證技術(shù)導(dǎo)則》（國海發(fā)[2010]22號）要求論證范圍應(yīng)覆蓋項目用?？赡苡绊懙降娜繀^(qū)域，并規(guī)定一般情況下一級論證的論證范圍以項目用海外緣線為起點向外擴(kuò)展15 km進(jìn)行劃定，該標(biāo)準(zhǔn)相當(dāng)于給出了尺度L0=15 km，鑒于此，本研究也取h=15 km.

2.4 海域使用強(qiáng)度空間特征歸納

借鑒空間自相關(guān)原理[23]，利用自然斷裂法，根據(jù)組內(nèi)差異最小、組間差異最大的準(zhǔn)則，對核密度柵格圖層數(shù)據(jù)進(jìn)行空間聚類分析，可以歸納出海域使用強(qiáng)度的相對高值區(qū)域，并將這些區(qū)域理解為海域使用負(fù)載的集聚區(qū)域.

3 試驗與結(jié)果

3.1 試驗區(qū)域與數(shù)據(jù)情況

本研究選擇渤海海域為試驗區(qū)，嘗試開展海域使用強(qiáng)度計算并分析其空間分異特征.試驗區(qū)涉及遼寧、河北、天津和山東4省（市）的13市.試驗區(qū)海域向陸一側(cè)的邊界為各?。ㄊ校┕嫉墓芾戆毒€.

試驗采用國家海域使用動態(tài)監(jiān)視監(jiān)測系統(tǒng)中的海域使用管理業(yè)務(wù)數(shù)據(jù)，包括2015年底的海域使用確權(quán)登記數(shù)據(jù)和海域使用本底庫數(shù)據(jù)，數(shù)據(jù)覆蓋試驗區(qū)內(nèi)所有已知的海域開發(fā)活動.

3.2 海域使用負(fù)載計算結(jié)果

試驗區(qū)內(nèi)總的海域使用負(fù)載如圖1所示，其中漁鹽業(yè)用海、工業(yè)城鎮(zhèn)用海、交通運輸用海和旅游娛樂用海4種海域使用類型的負(fù)載如圖2所示.

圖1 渤海區(qū)總海域使用負(fù)載的示意圖Fig.1 Schematic diagram for load of total sea area use in Bohai Sea

圖2 渤海區(qū)4種類型的海域使用負(fù)載示意圖Fig.2 Schematic diagram for load of four types of sea area use in Bohai Sea

圖1和圖2中的海域使用負(fù)載數(shù)據(jù)在一定程度上反映了海域使用強(qiáng)度空間分布的宏觀趨勢.由圖1可以看出，海域使用負(fù)載高值區(qū)主要出現(xiàn)在天津、曹妃甸和盤錦港等沿岸海域，此外，長興島、鲅魚圈和黃驊港-埕口等沿岸海域也有相對較高的負(fù)載值.結(jié)合圖2可知，圖1中突出的負(fù)載值主要源自漁鹽業(yè)用海、工業(yè)城鎮(zhèn)用海以及交通運輸用海的貢獻(xiàn)，特別是集聚在臨港工業(yè)綜合配套區(qū)內(nèi)的工業(yè)城鎮(zhèn)用海以及交通運輸用海.漁鹽業(yè)用海在渤海沿岸的分布相對廣泛，主要布局在河海交匯處和海島附近海域，連片開發(fā)特征明顯，但負(fù)載值不高，較高值集中在山東濱州-東營以及濰坊沿岸海域.工業(yè)城鎮(zhèn)用海呈多點集中分布，主要布局在天津、唐山、盤錦和大連等經(jīng)濟(jì)發(fā)達(dá)的地區(qū)沿岸，海域使用負(fù)載值較高，其中尤以天津臨港和南港工業(yè)區(qū)、曹妃甸、葫蘆島港和盤錦港等沿岸海域最為突出.交通運輸用海主要集中在天津港、曹妃甸、鲅魚圈和黃驊港等沿岸海域，以天津港和黃驊港沿岸海域使用負(fù)載最高.旅游娛樂用?？傮w負(fù)載較小，高值區(qū)主要集中在濰坊濱海生態(tài)旅游休閑娛樂區(qū)，該區(qū)域有較大規(guī)模的填海造地.

但海域使用負(fù)載數(shù)據(jù)僅反映所在點位直接的用?；顒?由圖1和圖2可以看出，海域使用負(fù)載數(shù)據(jù)在空間上并非連續(xù)分布，其數(shù)值分布與海域使用范圍直接對應(yīng)，位于海域使用范圍內(nèi)的點位上，對應(yīng)用海類型和用海方式的海域使用負(fù)載值，此外，海域使用負(fù)載值為0.因此，海域使用負(fù)載數(shù)據(jù)并不能綜合反映一定范圍內(nèi)海域使用空間分布的態(tài)勢，不具備區(qū)域代表性.

3.3 海域使用強(qiáng)度分析結(jié)果

為了反映一定范圍內(nèi)海域使用負(fù)載的綜合作用，本研究進(jìn)一步對海域使用負(fù)載數(shù)據(jù)層進(jìn)行核密度分析，并對所得核密度數(shù)據(jù)層進(jìn)行歸一化處理，得到海域使用強(qiáng)度數(shù)據(jù)層.為了自動識別海域使用強(qiáng)度的相對高值區(qū)域，為海域使用管控業(yè)務(wù)提供區(qū)域分類和指示信息，本研究在此基礎(chǔ)上利用自然斷裂法，按5級分類方案對該層數(shù)據(jù)進(jìn)行空間聚類分析，結(jié)果如圖3所示.由圖3可以看出，通過核密度分析所得海域使用強(qiáng)度數(shù)據(jù)在空間分布上保持與海域使用負(fù)載數(shù)據(jù)相同的趨勢，但在具體點位上，海域使用強(qiáng)度數(shù)據(jù)已與海域使用范圍無直接對應(yīng)關(guān)系，每一點的數(shù)值都綜合了一定范圍內(nèi)海域使用的負(fù)載作用，部分點位即使不在海域使用范圍內(nèi)，仍有海域使用強(qiáng)度值，因此整個數(shù)據(jù)層在空間上呈現(xiàn)連續(xù)分布的特征.

圖3 渤海區(qū)海域使用強(qiáng)度示意圖Fig.3 Schematic diagram for intensity of sea area use in Bohai Sea

按照本研究的設(shè)定可知，圖3中多個海域使用強(qiáng)度較高的聚類區(qū)域內(nèi)存在分布范圍較高、量值較高的海域使用負(fù)載，因此，這些區(qū)域可以理解為海域使用負(fù)載的集聚區(qū).結(jié)合上文中海域使用負(fù)載分析結(jié)論可知，渤海區(qū)海域使用負(fù)載集聚區(qū)主要呈現(xiàn)“高值核心式”聚集和“較高值平攤式”聚集2種布局方式，前者以集中分布的工業(yè)城鎮(zhèn)和交通運輸類填海造地用海為主要成因，后者以連片大規(guī)模漁鹽業(yè)圍海和分散的工業(yè)城鎮(zhèn)、交通運輸類填海造地為主要成因，具體分布如表1所示.

表1 渤海區(qū)主要的海域使用負(fù)載聚集區(qū)Tab.1 Main aggregated area of sea use in Bohai Sea

本研究所得海域使用強(qiáng)度分布特征與高志強(qiáng)等[24]和徐進(jìn)勇等[25]有關(guān)渤海區(qū)海岸線變化相對集中區(qū)和強(qiáng)度劇烈區(qū)的判斷相似，說明本研究海域使用強(qiáng)度評價方法的有效性.

4 結(jié)論與討論

4.1 結(jié)論

本研究針對海域使用管控業(yè)務(wù)需求，從海域使用的區(qū)域聚集程度、資源利用程度及對海域自然屬性的影響程度等方面概括了海域使用強(qiáng)度概念的內(nèi)涵，設(shè)計了一套用于表征海域使用強(qiáng)度空間分布特征的方法，并以2015年渤海區(qū)為例，利用海域使用動態(tài)監(jiān)測數(shù)據(jù)開展了海域使用強(qiáng)度評估，研究結(jié)果表明：

（1）利用本研究方法所得結(jié)論與已有認(rèn)知和研究成果基本一致，可以反映渤海區(qū)海域使用強(qiáng)度的空間分布特征，且評估方法簡單、有效.

（2）渤海區(qū)海域使用強(qiáng)度空間分異明顯，有多個集聚區(qū)域，呈現(xiàn)“高值核心式”聚集和“較高值平攤式”聚集2種分布態(tài)勢.“高值核心式”聚集主要出現(xiàn)在盤錦港（遼濱經(jīng)濟(jì)區(qū)）和天津-曹妃甸沿岸海域；“較高值平攤式”聚集主要出現(xiàn)在長興島-普蘭店、老母豬礁、鲅魚圈、葫蘆島港-大凌河、興城臨海產(chǎn)業(yè)區(qū)-長山寺、灤河口、大清河口、黃驊港-埕口、東營港、黃河口、濰坊、小石島和龍口港沿岸海域.

4.2 討論

（1）本研究旨在為海域資源承載力監(jiān)測預(yù)警、海域使用管控等工作提供可用于業(yè)務(wù)運行的海域使用強(qiáng)度評估方法，方法的便捷性和結(jié)果的有效性是主要目標(biāo).從研究結(jié)果來看，本研究的評估思路和方法為大尺度范圍內(nèi)海域使用強(qiáng)度的快速高效評估提供了可行途徑.

（2）目前，本研究針對于海域資源利用程度系數(shù)和海域自然屬性影響程度的取值主要來源于專家意見，評價結(jié)果有一定的主觀性成分，相關(guān)研究有待深入.

（3）根據(jù)他人研究成果可知[26]，核函數(shù)的選擇及衰減距離閾值的取值對海域使用負(fù)載核密度計算值存在一定影響，相關(guān)研究有待深入.從本文研究結(jié)果看，它們不會改變海域使用負(fù)載空間密度的空間分布趨勢.由于本文最終采用歸一化的海域使用負(fù)載空間密度表征海域使用強(qiáng)度，因此兩者對海域使用強(qiáng)度評估結(jié)果的影響有限.

（4）海域柵格化時，柵格大小的選擇會對海域使用強(qiáng)度評估的空間精細(xì)程度和計算效率產(chǎn)生影響，較小的柵格尺度導(dǎo)致較高的空間精細(xì)程度和較低的計算效率，但不會改變海域使用負(fù)載空間密度的空間分布趨勢.本研究認(rèn)為現(xiàn)實中較小用海工程的海域使用范圍基本在數(shù)十米尺度，按數(shù)十米尺度進(jìn)行海域柵格化，評估的空間精細(xì)程度和計算效率均較適宜.