深圳市1988-2018年海岸線時空演變分析

柏葉輝，李向新，鐘舒怡

?

深圳市1988-2018年海岸線時空演變分析

柏葉輝，李向新，鐘舒怡

（昆明理工大學 國土資源工程學院，云南 昆明 650000）

以深圳為研究區域并對其進行岸灘分區，利用1988年、1998年、2008年和2018年四個時相的TM和OLI衛星影像數據，利用面向對象的方法提取各期海岸線，系統分析海岸線時空演化特征和岸線開發利用格局變化。結果表明：研究時段內，深圳市海岸線長度先增長后縮短，海岸線的整體變化強度為0.39%，西部地區海岸線變化強度要明顯大于東部地區。大量的自然岸線轉變為平直的人工岸線，人工岸線在逐年增加，自然岸線在不斷縮減，深圳灣岸段現已基本都已變為人工岸線。與此同時，海岸線的開發強度持續增長，且表現出不斷的上升趨勢，由1988年的0.32變為2018年的0.53，深圳灣和珠江口岸段開發強度明顯高于大鵬灣和大亞灣岸段。

海岸線；時空演變；遙感；大數據分析

0 引言

海岸帶地區是海洋系統與陸地系統之間相互作用形成的獨特生態過渡帶[1]，其生態環境極其敏感、脆弱，易受到自然變化與人類活動的影響[2,3]。海岸線不僅標識了沿海地區的水陸分界線，而且蘊含著豐富的環境信息，其變化直接影響潮間帶灘涂資源量及海岸帶環境，將引起海岸帶多重資源與生態環境過程的改變，影響沿海人民的生存發展[4]。在自然因素和人類活動的影響下，海岸線處于不斷的變化中[5]，這種變化既是海岸對各種動力作用的響應，也是海岸環境演變的直接體現[6]。因此，快速準確地監測海岸線變化，從而為決策部門提供科學的、有效的信息[7,8]，對進行海岸資源開發利用和海洋環境災害預防與評估等都具有十分重要的意義[9,10]。

伴隨著大數據時代的來臨，包括大數據分析在內的大數據理念和方法應用于諸多方面[11-14]，遙感大數據擁有時空、光譜、分辨率等數據特征和瀏覽、顯示、存取等應用特點[15]，對于多時相數據分析具有明顯的優勢。遙感技術擁有強大的獲取數據的能力，在海岸線調查中具有明顯的優勢[5]，通過對遙感影像的解譯成為海岸線提取的新手段[16,17]。面向對象法是一種興新的遙感圖像解譯方法，該方法通過對影像的分割，使同質像元組成大小不同的對象，從而實現較高層次的遙感圖像分類和目標地物提 取[18]。本文以深圳市海岸帶為研究區域，并進行岸灘分區，采用LandsatTM/OLI影像為數據源，采用面向對象的方法，建立分類規則集提取海岸線，系統的分析了海岸線動態演變特征，包括海岸線長度和海岸類型的變化，以及對各海岸開發利用強度進行評價并探討其原因，以期為深圳市海洋的可持續發展提供參考。

1 材料與方法

1.1 研究區概況

深圳市（22°27′―22°52′N，113°46′―114°37′E）是中國南部濱海城市，位于廣東省東南部，西臨珠江口，東瀕大亞灣和大鵬灣，南接香港新界，北靠東莞和惠州兩市。地貌類型多樣，地勢東南高、西北低，海岸類型豐富，其海岸帶被九龍半島分割為東、西兩部分，東部大鵬灣、大亞灣岸線曲折多為基巖質和沙礫質，西部深圳灣和珠江口主要為淤泥質和人類圍墾修筑岸線。海岸地貌類型的多樣化，使得海岸帶開發呈現多元化的特征。本文參考了各國開展海岸帶調查時常用的海岸線向陸延伸10 KM的空間范圍作為海岸帶調查范圍，將深圳市行政界線海域一側界限為基準向內陸一側10 KM緩沖區作為研究區（圖1）。由于研究內容的需要，本文將根據地貌的差異對研究區域進行岸灘分區，分為深圳灣、大鵬灣、大亞灣、珠江口四個研究單元。

1.2 研究數據來源

本文研究采用的主要數據源包括30m空間分辨率的1988、1998、2008年的Landsat TM衛星影像和2018年的Landsat OLI衛星影像，共四期12景；深圳市地形圖及相關資料。針對四期影像數據，利用ENVI軟件對其進行了預處理，包括幾何校正、影像融合以及影像增強等，分別求得各自的緩沖區覆蓋范圍。

圖1 研究區范圍

1.3 海岸線提取

鑒于Landsat衛星影像過境為當地時間上午10點左右，研究區潮位較高，衛星影像可用于提取海岸線。本研究針對不同類型的岸段分別提取，在連接各研究單元內的提取結果。參照馬小峰[19]，孫偉福[20]，孫美仙[7]等人提出的海岸線提取方法，對于基巖和人工海岸，根據其對光譜的高反射特性，選擇灰度值區分海陸，其邊界作為岸線；沙質海岸提取灘脊痕跡線上限；對于淤泥質岸線，利用NDVI區分濕生植被與灘涂，靠陸一側分界線作為岸線。所有岸線都是利用eCognition 軟件根據不同地物在影像上表現出不同的光譜，形狀和紋理特征來提取出來的，同時結合Googel Earth高清影像對所有提取結果進行人工目視解譯并修正。最終分出四期海岸線。

1.4 岸線長度變化強度

為了更客觀的的對比各研究單元不同時段海岸線長度變化的時空差異，再此采用某一時段內某一研究單元海岸線長度的年均變化百分比來表示其海岸線的變化強度[10]。

1.5 岸線開發利用強度

岸線開發利用強度是定量表征不同的海岸類型對海岸帶資源的影響強弱的量，能夠更好的了解其空間分布規律，挖掘海岸線資源的開發利用潛力，有助于海岸帶資源的可持續開發利用[21]。具體公式如下：

表1 各類型岸線的資源環境影響因子

Tab.1 Resources and environment impact factor of all kinds of the shoreline

2 研究結果與討論

2.1 岸線時空演化特征

2.1.1 岸線長度演變特征

根據遙感監測提取的海岸線成果顯示（圖2），1988年至2018年間，深圳市海岸線發生較大的變化，且不斷的向海推進，深圳市西部地區海岸線變化在范圍和程度上都大于東部地區，特別是西部的珠江口岸區和深圳灣岸區變化最為明顯，東部的大鵬灣和大亞灣僅部分岸段發生明顯變化。

圖2 1988-2018年深圳市海岸線動態變化

不同岸段的長度及變化也有顯著地差別。在這30年間，深圳市海岸線由213.70 km變為了239.02 km，增長了25.32 km。其中，珠江口岸區增長了18.1km，深圳灣岸區增長了0.99 km，大鵬灣岸區增長了7.74 km，大亞灣岸區減少了1.51 km。深圳市1998年較1988年增長了17.78 km，平均速度1.79 km/a，此階段岸線增長速度最快；2008年較1998年增長了14.18 km，平均速度1.42 km/a；2018年較2008年減少了6.64 km，平均速度0.66 km/a。

1988至2018年深圳市四個岸段海岸線長度一直保持大亞灣岸段>大鵬灣岸段>珠江口岸段>深圳灣岸段的分布格局。其主要原因是各岸段的地理位置差異所導致的海岸線長度有差異，因此，可以采用1988-1998年、1998-2008年、2008-2018年、1988-2018年三個時段和大亞灣岸段、大鵬灣岸段、珠江口岸段、深圳灣岸段四個研究單元海岸線長度的年均變化百分比來表示其海岸線的變化強度（圖3）。

圖3 1988-2018年深圳市各岸灘海岸線變化強度

1988-2018年研究區整體海岸線的變化強度為0.39%，從時間序列上看，1988-1998年是研究去海岸線長度變化強度最大的時期，海岸線變化強度為0.83%；其次為1998-2008年，海岸線變化強度為0.61%；2008年后海岸線進入快速減少時期，2008-2018年海岸線開始縮減，其變化強度為-0.27%。從研究單元看，1988-2018年珠江口岸段海岸線的長度變化最為劇烈，海岸線變化強度為1.36%；其次是大鵬灣，為0.43%；再次是深圳灣和大亞灣，海岸線變化強度分別為0.12%和-0.06%。

對于西部地區，1998-2008年是珠江口岸段海岸線變化最為激烈的一個時期，變化強度為2.01%；其次為1988-1998年，海岸線變化強度為1.53%；2008-2018年，海岸線變化強度最小，為0.17%。 1988-1998年是深圳灣岸段海岸線變化最為明顯的一個時期，變化強度為0.65%，其次為1998-2008年，海岸線變化強度為-0.36%，2008-2018年，海岸線變化強度最小，為0.09%。總體而言，西部地區海岸線變化強度較大，主要原因是珠江口區域，有大量的泥沙不斷淤積使灘涂向海推進[22]，岸線平直，岸線平緩，是灘涂圍墾的主要區域。深圳灣是一個半封閉型的淺水海灣，主要為淤泥質岸線，是填海造陸的主要區域。

對于東部地區，1988-1998年是大鵬灣岸段海岸線變化最強的一個時期，變化強度為1.39%；其次為2008-2018年，變化強度為-0.44%；變化最小為1998-2008年，變化強度為0.39%。2008-2018年是大亞灣岸線海岸線變化最為明顯的一年，變化強度為-0.57%；其次為1998-2008年，變化強度為0.29%；變化最小的為1988-1998年，變化強度為0.11%。總體而言，東部地區海岸線變化強度沒有西部大，主要原因大鵬灣和大亞灣岸段主要為基巖岸線，只有部分區域建設碼頭和網箱養殖，岸線小幅延伸，海岸線由于人為建設裁彎取直。

2.1.2 各類型岸線變化特征

根據遙感影像特點，將1988、1998、2008、2018年各時期的海岸線分為自然岸線和人工岸線，其中自然岸線分為生物岸線、基巖岸線、沙質岸線、淤泥質岸線，人工岸線包括建設圍堤岸線、養殖圍墾岸線、農業圍墾岸線。由于人類對岸線的圍墾開發，近30年來，深圳市海岸帶的自然岸線呈現出不斷縮短的趨勢，而人工岸線不斷的增加。1988年，深圳市人工岸線長度為83.98 km，占總岸線的37%，經過30年的城市發展，截止到2018年，人工岸線長度已達135.6 km，占總岸線的57%，比例上升了20%。與之相反的是，深圳市自然岸線在不斷的減少，且減少速度也逐漸變緩，1988-1998階段速度為-2.62 km/a，為30年中減少最快的階段；1998-2008階段速度為-1.24 km/a，2008-2018階段速度為-0.32 km/a。

具體各類型岸線的變化見圖4。從圖看出，自然岸線中，基巖岸線呈現出逐年下降的趨勢，更多的基巖海岸被開發用來建設城市住宅用地和城市工業用地等，30年間，基巖海岸共減少了19.39 km，減少速度達到了0.65 km/a。沙質海岸岸線的長度在先減少后增加后，最終2018年與1988年基本一致，這是由于隨著旅游業的不斷發展，岸線受到了更多的保護和利用，減少沙粒被海浪沖蝕的現象。淤泥質海岸呈現出逐年減少的趨勢，這主要是由于更多的淤泥質海岸被人為圍墾開發，轉化為建設圍堤岸線、養殖基地或農田。

圖4 1988-2018年深圳市各類型岸線長度變化

人工岸線中，農業圍墾岸線在逐年減少，由1988年的33.27 km，變為2018年的8.81 km。而建設圍堤岸線卻在不斷增加，到2008年達到最大值，且2018年和2008年基本持平。說明隨著工業、旅游業及海上貿易的發展，人類對土地的需求量不斷的增加，越來越多的農業用地被用來發展工業。養殖圍墾岸線從1998年后基本維持不變。

2.2 岸線開發強度評價

依據岸線開發利用強度公式計算出深圳市海岸線開發強度（圖5），在過去的30年里，深圳市海岸線的開發強度持續增長，且表現出不斷的上升趨勢，由1988年的0.32提高為2018年的0.53，這表明人類對岸線資源的開發利用越來越迫切。各岸段的岸線開發利用強度也存在顯著的差異，具體來說，深圳灣岸段開發強度一直很高，珠江口岸段開發強度緊隨其后，但最近幾年超越了深圳灣岸段，大鵬灣和大亞灣岸段岸線開發強度都較低，但隨著時間的變化，也呈現出不斷增大的趨勢。

圖5 深圳市1988-2018年各岸段海岸線開發強度指數

總的來說，深圳市西部海岸的珠江口岸段和深圳灣岸段的開發利用強度要明顯高于東部海岸的大鵬灣岸段和大亞灣岸段。其中岸線開發利用強度最大的為這主要因為深圳市東西部海岸線的自然條件存在明顯分異，西部海岸線幾乎全部為珠江河口地區的淤泥質海岸線，而東部海岸線大部分為基巖海岸和沙質海岸[23]。西部地區造地成本低，填海條件優越，現已全部被改造為生產和生活岸線，而東部地區由于岸線底質條件大部分都不適合填海，開發的較慢，現生態岸線居多，只有部分由于城市化和工業化被改造為生產和生活岸線。

3 結論

本文利用RS和GIS技術獲取了1988年-2018年30年間深圳市“三灣一口”的海岸線分布情況；計算了岸線長度變化強度和岸線開發利用強度，分析了海岸線時空長度、類型變化特征以及海岸線開發利用強度。主要結論如下：

（1）1988-2008年研究區海岸線長度持續增加，2008-2018年研究區海岸長度線有所輕微減少。1988至2018年深圳市四個岸段海岸線長度一直保持大亞灣岸段>大鵬灣岸段>珠江口岸段>深圳灣岸段的分布格局。1988-2018年研究區整體海岸線的變化強度為0.39%，西部地區海岸線變化強度要明顯大于東部地區海岸線變化強度，主要是人為開發原因造成的。

（2）由于人類對岸線的圍墾開發，近30年來，深圳市海岸帶的自然岸線在逐年縮減，而人工岸線卻在逐年增加，大部分自然岸線轉化為人工岸線，人工岸線中的各類別也會相互轉換，但主要表現為農業圍墾岸線轉化為建設圍堤岸線。自然岸線中，基巖岸線、淤泥質岸線呈現出逐年下降的趨勢；人工岸線中，農業圍墾岸線在逐年減少，而建設圍堤岸線卻在不斷增加。但是不同的四個岸段也各自呈現出不同的變化特點。

（3）1988-2018年，深圳市海岸線的開發強度持續增長，且表現出不斷的上升趨勢，由1988年的0.32提高為2018年的0.53。深圳市西部海岸的珠江口岸段和深圳灣岸段的開發利用強度要明顯高于東部海岸的大鵬灣岸段和大亞灣岸段，其原因與不同類型海岸的地貌環境差異密不可分。

[1] 王棟, 買合木提, 玉永雄, 等. 我國海岸帶生態現狀研究進展[J]. 河北漁業, 2007(09): 10-13.

[2] 陳夢熊. 沿海地區地質環境特征與地質環境系統--兼論“人地系統”[J]. 中國地質災害與防治學報, 1998(S1): 84-90.

[3] 宋玉環張壽全王思敬. 中國沿海地質環境與區域持續發展的若干問題探討[J]. 工程地質學報, 1996(03): 24-29.

[4] 蔡則健吳曙亮. 江蘇海岸線演變趨勢遙感分析[J]. 國土資源遙感, 2002(03): 19-23.

[5] 劉寶銀, 蘇振奮. 中國海岸帶與海島遙感調查[M]. 海洋出版社, 2005.

[6] 黃鵠, 胡自寧, 陳新庚, 等. 基于遙感和GIS相結合的廣西海岸線時空變化特征分析[J]. 熱帶海洋學報, 2006(01): 66-70.

[7] 孫美仙張偉. 福建省海岸線遙感調查方法及其應用研究[J]. 臺灣海峽, 2004(02): 213-218.

[8] 李學杰. 應用遙感方法分析珠江口伶仃洋的海岸線變遷及其環境效應[J]. 地質通報, 2007(02): 215-222.

[9] 于永海苗豐民張永華尤玉明李四海. 遼—冀與津—冀臨界區海岸演變衛星遙感分析[J]. 國土資源遙感, 2003(03): 30-32.

[10] 徐進勇, 張增祥, 趙曉麗, 等. 2000-2012年中國北方海岸線時空變化分析[J]. 地理學報, 2013, 68(05): 651-660.

[11] 馬凱航, 高永明, 吳止鍰, 等. 大數據時代數據管理技術研究綜述[J]. 軟件, 2015, 36(10): 46-49.

[12] 賀寧, 湯明偉, 賀萌. 大數據課程在高職學院中的建設和發展[J]. 軟件, 2015, 36(10): 50-52.

[13] 諶迅. 大數據資產管理系統的設計與實現[J]. 軟件, 2016, 37(02): 50-53.

[14] 李向陽. 大數據時代國防科技信息工作機遇和挑戰[J]. 軟件, 2018, 39(01): 101-104.

[15] 錢育蓉, 于炯, 英昌甜, 等. 云計算環境下新疆遙感應用數據中心的挑戰與機遇[J]. 軟件, 2015, 36(04): 58-61.

[16] Moigne J L, Tilton J C. Refining image segmentation by integration of edge and region data[J]. Geoscience & Remote Sensing IEEE Transactions on, 1992, 33(3): 605-615.

[17] Dellepiane S, Laurentiis R D, Giordano F. Coastline extraction from SAR images and a method for the evaluation of the coastline precision[J]. Pattern Recognition Letters, 2004, 25(13): 1461-1470.

[18] 陳云浩, 馮通, 史培軍, 等. 基于面向對象和規則的遙感影像分類研究[J]. 武漢大學學報(信息科學版), 2006(04): 316-320.

[19] 馬小峰, 趙冬至, 邢小罡, 等. 海岸線衛星遙感提取方法研究[J]. 海洋環境科學, 2007(02): 185-189.

[20] 孫偉富, 馬毅, 張杰, 等. 不同類型海岸線遙感解譯標志建立和提取方法研究[J]. 測繪通報, 2011(03): 41-44.

[21] 劉百橋, 孟偉慶, 趙建華, 等. 中國大陸1990—2013年海岸線資源開發利用特征變化[J]. 自然資源學報, 2015, 30(12): 2033-2044.

[22] 朱大奎. 深圳海岸海洋環境與空間規劃研究[J]. 第四紀研究, 2005(01): 45-53.

[23] 彭建王仰麟劉松吳健生李衛鋒. 海岸帶土地持續利用景觀生態評價[J]. 地理學報, 2003(03): 363-371.

Dynamic Change of the Coastline in Shenzhen during 1988-2018

BAI Ye-hui, LI Xiang-xin, ZHONG Shu-yi

(Kunming University of Science and Technology, Institute of land and Resources Engineering, kunming, Yunnan, 650000)

Taking Shenzhen, a region of rapid urbanization, as the study area and dividing the shoreline into different zones, using TM and OLI satellite image data of 1988, 1998, 2008 and 2018, the shoreline of each period was extracted by object-oriented method, and the temporal and spatial evolution characteristics of the shoreline and the pattern changes of shoreline development and utilization were analyzed systematically. The results show that the coastline length of Shenzhen increases first and then decreases. The overall change intensity of the coastline is 0.39%. The change intensity of the coastline in the western region is obviously greater than that in the eastern region. A large number of natural shorelines have been transformed into flat artificial shorelines. Artificial shorelines have been increasing year by year. Natural shorelines have been shrinking. At the same time, the development intensity of the coastline continues to grow, and shows a rising trend, from 0.32 in 1988 to 0.53 in 2018. The development intensity of Shenzhen Bay and Pearl River estuary is obviously higher than that of Dapeng Bay and Daya Bay.

Coastline; Dynamic change; Remote sensing; Big data analysis

TP79

A

10.3969/j.issn.1003-6970.2018.10.026

柏葉輝(1994-)，碩士研究生，主要研究方向為海岸帶生態環境變化；李向新(1963-)，副教授，主要研究方向為3S技術集成及工程應用、GIS水文模型；鐘舒怡(1993-)，碩士研究生，主要研究方向為城市擴張。

柏葉輝，李向新，鐘舒怡. 深圳市1988-2018年海岸線時空演變分析[J]. 軟件，2018，39（10）：133-138