王兆蘭，王利，2，趙彪

(1. 遼寧師范大學城市與環(huán)境學院 大連 116029； 2. 遼寧師范大學海洋經(jīng)濟可持續(xù)發(fā)展研究中心 大連 116029)

?

近年來大連市圍填海區(qū)域動態(tài)演化研究*

王兆蘭1，王利1，2，趙彪1

(1. 遼寧師范大學城市與環(huán)境學院 大連 116029； 2. 遼寧師范大學海洋經(jīng)濟可持續(xù)發(fā)展研究中心 大連 116029)

基于遙感和GIS技術，以2000年、2005年、2010年LandsatTM數(shù)據(jù)和2014年Landsat OLI數(shù)據(jù)共4期的遙感影像為數(shù)據(jù)源，完成了對大連市海岸線變遷及圍填海區(qū)域演變信息的提取，分析了大連市圍填海區(qū)域的動態(tài)演化特征，并探討了圍填海區(qū)域變化的驅動因素。研究表明：2000—2014年，大連市海岸線長度以及圍填海面積逐年增加，2005年以后進入快速增長時期，但增速呈由快變慢的趨勢；圍填?；顒又饕性谕叻康旰褪休爡^(qū)岸段；人工岸線增加，岸線年均長度變化強度由大到小依次為普蘭店岸段、莊河岸段、瓦房店岸段、市轄區(qū)岸段；從驅動因素上看，圍填海面積變化是主要是社會經(jīng)濟因素作用的結果，與養(yǎng)殖業(yè)、工業(yè)以及港口物流業(yè)的發(fā)展密切相關。

圍填海區(qū)域；填海造地；海岸線；動態(tài)演化；大連市

1 引言

大連市是三面環(huán)海的半島城市，是東北地區(qū)對外開放的窗口和最大的港口城市。隨著經(jīng)濟社會的快速發(fā)展，大連市人多地少的矛盾日益突出，面對日益嚴重的“土地赤字”問題，圍填海成為緩解用地緊張的重要途徑。圍填海是指通過人類在海岸線外進行圍海和填?；顒?，使指定海域失去海洋屬性，從而對其進行有效利用的方式[1]；大規(guī)模的圍填海工程使得海岸帶發(fā)生巨大變化，直接改變潮間帶灘涂資源量及海岸帶環(huán)境，引起海岸帶多種資源與生態(tài)過程的改變，影響沿海人民的生存與發(fā)展[2-3]。因此，及時準確地監(jiān)測沿海地區(qū)圍填海變化情況，對于研究圍填?；顒訉Ｑ蟓h(huán)境造成的影響，推進沿海地區(qū)可持續(xù)發(fā)展具有重要意義。

近年來國內(nèi)外學者對圍填海區(qū)域變化及其驅動因素進行了深入的研究，國外部分學者分別對英國東部海岸、土耳其東部地區(qū)海岸帶、阿聯(lián)酋阿布扎比海岸帶、尼羅河三角洲、韓國 Gomso 灣淤泥質海岸、印度卡納塔克海岸、埃及尼羅河三角洲東北部海岸進行了研究[4-8]。國內(nèi)學者的研究主要集中在對渤海灣西北岸、遼寧省、廈門市、深圳灣、膠州灣等地區(qū)的圍填海情況進行了監(jiān)測研究[9-13]。通過對相關文獻梳理發(fā)現(xiàn)，現(xiàn)有研究多是從省區(qū)的層面研究海岸帶的變化，對微觀海域面積的研究較少，缺乏對大連市近年來的圍填海面積變化的研究；解譯方法較為單一，對用地類型分辨較為模糊。因此，本研究以大連市為研究對象，結合遙感影像和Google Earth實時影像，采用半自動-交互式目視判讀方法進行解譯，對2000-2014年大連市圍填海動態(tài)變化進行監(jiān)測研究，以期為大連市海洋資源合理利用、環(huán)境保護、城市規(guī)劃等科學決策提供相關依據(jù)。

2 數(shù)據(jù)來源與研究方法

2.1 研究區(qū)與數(shù)據(jù)源

大連市東臨黃海，西鄰渤海，南與山東半島隔海相望，北與東北內(nèi)陸相連，處于121°44′E—121°49′E、39°01′N—39°04′N。海岸線范圍為瓦房店市浮渡河口至莊河市南尖鎮(zhèn)，海岸類型為基巖海岸、沙礫質海岸、岬灣沙礫海岸、岬灣淤泥海岸。近年來，隨著大連海洋經(jīng)濟的快速持續(xù)發(fā)展，沿海地區(qū)的土地資源緊缺，大連岸段海岸線的變化基本上源于圍填海工程。本研究采用2000年、2005年、2010年TM影像(30 m)、2014年OLI影像(30 m)共4個時期覆蓋全區(qū)的16景Landsat影像，以及大連政區(qū)圖、地形圖、海圖、大連市城市總體規(guī)劃(2010-2020年)、《大連市人民政府辦公廳關于印發(fā)大連海洋功能區(qū)劃通知》[大政辦(2006)165號]和歷年統(tǒng)計年鑒為數(shù)據(jù)來源。由于圍填海區(qū)域的岸線一般為人工岸線，成像季節(jié)對研究影響不大，因此選取的16景影像成像時間均為9-10月，且均為無云的高質量影像。

2.2 數(shù)據(jù)處理

綜合研究區(qū)的海岸類型和各種海岸線提取方法的利弊，本研究選取人工目視解譯法，遙感影像目視解譯法提取海岸線的誤差可控制在30 m內(nèi)，面積量算誤差為6 km2[21]，因此采用此法勾出各時相海岸線可行。本研究對圍填海區(qū)域的提取主要是在GIS技術的支撐下，并參照其他輔助數(shù)據(jù)采用半自動-交互式目視判讀解譯。圍填海區(qū)域的提取具體實施步驟為：在RS影像處理技術的支持下對16景影像分別做快速大氣校正，以便對地物光譜特征的識別。TM影像選取5、4、3波段，OLI數(shù)據(jù)選取6、5、4波段進行彩色合成，這樣的組合有利于陸地/水體的目視判讀。對已得到的全部影像分幅做幾何校正和地形校正，之后對每一時相的四景影像進行鑲嵌，形成一幅覆蓋大連全域的影像文件，并且坐標統(tǒng)一為WGS84 坐標系統(tǒng) UTM-51N投影。最后以2000年鑲嵌完整的影像作為基準影像對其他時相影像做配準，對配準后的各個時相的影像，在參考大連行政區(qū)劃圖以及海圖的基礎上進行裁剪，確定研究區(qū)域界限。本研究將鹽田和養(yǎng)殖池視為不透水構筑物，以鹽田和養(yǎng)殖池圍堤向海測的平均大潮高潮時水陸分界的痕跡線作為人工岸線[20]，且不考慮防波堤和島嶼構成的岸線。利用處理后30 m分辨率的影像進行海岸線的目視解譯，并勾勒出各個時相大連海岸線。在GIS和Excel等軟件的支持下，完成圍填海面積變化的定量分析、空間分析以及表達。

3 圍填海動態(tài)變化結果分析

依據(jù)大連行政區(qū)劃，將大連海岸分為4個岸段，即瓦房店岸段(包括長興島臨港工業(yè)區(qū))、普蘭店岸段、市轄區(qū)岸段(包括中山區(qū)、西崗區(qū)、沙河口區(qū)、甘井子區(qū)、旅順口區(qū)、金州新區(qū))、莊河岸段。根據(jù)遙感影像提取結果，分別統(tǒng)計了各個岸段以及整體岸段的海岸線長度和圍填海面積的變化量(表1)。

表1 2000-2014年大連市海岸線長度和圍填海面積變化量

3.1 大連市海岸線長度的時空變化特征

2000—2014年間大連海岸線長度呈持續(xù)增加的趨勢，從2000年的1 217.5 km增加到2014年的1 439.6 km，共增加222.1 km，年均增加15.6 km。其中，2005—2010年間是大連海岸線增長最為劇烈的時期，共增加了127.9 km，占總增加長度的58%；其次為2000—2005年間，增加61.3 km，占總增加長度的28%；2010—2014年海岸線長度增加最小，占總增加長度的14%。由表1可以看出，4個岸段不同時期的海岸線長度均呈持續(xù)增加趨勢；瓦房店岸段的岸線長度增加主要集中在2000—2010年的10年間， 2010年以來，此岸段的岸線長度增量較小，僅占14年間的7%；市轄區(qū)岸段2000—2005年岸線長度增量最大，接下來的9年變化較小年均且變化程度相當；莊河岸段的岸線長度增加主要發(fā)生在2005—2010年間；普蘭店岸段的岸線長度增量發(fā)生在2005—2010年和2010—2014年的兩個時段內(nèi)且各占一半，而2000—2005年間普蘭店岸段的海岸線長度無變化。

根據(jù)大連市鄉(xiāng)鎮(zhèn)街道的行政界線，將大連海岸線更加細致地分為各個鄉(xiāng)鎮(zhèn)街道岸段(將大連金州灣國際機場所形成的岸線合并到金州城區(qū)岸段內(nèi))。根據(jù)式(1)[22]來計算2000—2014年各個岸段的海岸線變化強度，并制作海岸線變化強度空間分布圖(圖1)，以更加客觀的分析各岸段的海岸線的變化強度以及空間分布特征。

(1)

式中：LCIij為第i年至第j年海岸線長度變化強度 (length change intensity)；Li為第i年海岸線的長度；Lj為第j年海岸線長度。

從空間上看，2000—2014年間，大連市海岸線長度年均變化強度為1.30%。各岸段的海岸線年均長度變化強度由大到小依次呈現(xiàn)普蘭店、莊河、瓦房店、市轄區(qū)的空間分布特征，變化強度分別為2.49%、2.41%、1.65%、0.55%。

3.2 大連市圍填海區(qū)域的時空變化特征

從時間方面看(表1和圖2)，2000—2014年大連圍填海面積共增加了387.07 km2，年均增長27.64 km2，總體呈先升后降的趨勢。2005—2010年為圍填海面積迅速增加時期，約為前后兩個時期增加面積的總和，占14年間總增加面積的47%；這一時期，圍填?；顒又饕性谕叻康旰颓f河岸段，分別占總增量的43%、34%。其次為2010—2014年間，占增加面積的34%，這一時期圍填海面積增加主要集中在瓦房店岸段，占總增量的62%，其他岸段均有較少分布。2000—2005年圍填海增加面積最小，占14年間總增加面積的19%，圍填?；顒又饕性谕叻康旰褪休爡^(qū)岸段，分別占總增加量的50%、41%。

圖2 2000—2014年大連各岸段圍填海面積變化情況

從空間方面看，各岸段圍填海特征主要有：瓦房店岸段3個時期圍填海分布較集中，主要分布在兩個岸段，即李官鎮(zhèn)至永寧鎮(zhèn)岸段(包括李官鎮(zhèn)、土城鄉(xiāng)、永寧鎮(zhèn)海域)和長興島至復州灣鎮(zhèn)岸段(包括長興島、西中島、交流島鄉(xiāng)、謝屯鎮(zhèn)、復州灣鎮(zhèn)海域)。市轄區(qū)圍填海主要分布在金州區(qū)的三十里堡街道、金州城區(qū)岸段、大連開發(fā)區(qū)岸段，大連市內(nèi)的甘井子區(qū)大連灣街道兩側、凌水街道岸段和中山區(qū)北部海域也有圍填海工程分布。莊河岸段圍填海分布較均勻，相對來說大鄭鎮(zhèn)岸段、莊河城區(qū)西部岸段、鞍子山鄉(xiāng)岸段、栗子房鎮(zhèn)岸段圍填海面積較大一些。普蘭店各岸段圍填海分布也是較均勻，相對來說皮口鎮(zhèn)岸段分布較多一些。

4 2000年以來大連市圍填海面積變化原因分析

圖3 2000-2014年大連市圍填海利用類型

參照《海域使用分類體系》，將圍填海的利用類型分為5類，即養(yǎng)圈用海、工業(yè)用地、港口用地、城市生活及服務業(yè)用地、其他用地(包括旅游用地，房地產(chǎn)開發(fā)用地，用途不明用地)。根據(jù)解譯結果(圖3)，2000—2014年間，養(yǎng)圈海是最大的利用類型，占總圍填海面積的67.36%。由于工業(yè)用地、港口用地、城市生活及服務業(yè)用地以及其他用地屬于填海造地區(qū)域，固填海造地占14年間總圍填海面積的32.64%。養(yǎng)殖用海主要分布在瓦房店岸段(西楊鄉(xiāng)、仙浴鎮(zhèn)灣、三臺子-金州葫蘆套屯鹽田養(yǎng)殖功能區(qū))，市轄區(qū)岸段(石河街道-大魏家街道)，普蘭店岸段(楊樹房鎮(zhèn))，莊河岸段(碧流河口-莊河達拉腰子養(yǎng)殖綜合功能區(qū)，青堆子灣養(yǎng)殖功能區(qū))；填海造地中的工業(yè)用地所占比例較大且分布廣泛，主要是臨港工業(yè)區(qū)拓展以及國有廢棄鹽田上工業(yè)區(qū)的建設等；其次為港口用地，主要有旅順新港擴建工程，大窯灣港、莊河港等港口建設等；再次為其他用地，其主要是用來房地產(chǎn)開發(fā)和旅游業(yè)建設，例如金渤海岸的房地產(chǎn)用地，金石國際運動中心的旅游用海等；所占比例最小的是城市生活及服務業(yè)用地，主要有凌水灣高新園區(qū)總部經(jīng)濟基地、大連金州灣國際機場、小窯灣航運綜合服務區(qū)、東港城市公共服務區(qū)的建設等；14年間填海造地工程大量分布在市轄區(qū)岸段和瓦房店岸段。

基于GIS數(shù)據(jù)分析工具，提取2000—2014年變化的區(qū)域，將變化區(qū)域與2014年影像疊加，并結合Google Earth、《大連海洋功能區(qū)劃》以及實地調查解譯圍填海區(qū)域的利用類型以及空間分布情況(圖4)。

圖4 2000—2014年大連市主要填海區(qū)域及具體利用類型

近年來，大連著力發(fā)展海洋經(jīng)濟，沿海經(jīng)濟帶的經(jīng)濟水平大大加強，在海洋與漁業(yè)強市戰(zhàn)略推動下，2013年大連市實現(xiàn)海洋經(jīng)濟總產(chǎn)值2 652億元[23]4年實現(xiàn)2 990億元的海洋經(jīng)濟總產(chǎn)值。在經(jīng)濟因素和城市發(fā)展的驅動下，海洋漁業(yè)，沿海工業(yè)不斷發(fā)展，城市建設用地規(guī)模不斷擴張。海水養(yǎng)殖產(chǎn)量從2000年的118.6萬t增加到2013年的141.56萬t，由此反映出此期間養(yǎng)殖池的修建，養(yǎng)殖業(yè)的大力發(fā)展，勢必導致海域面積的減少，圍海利用的面積大量增加。例如，在瓦房店、莊河岸段進行了大規(guī)模的海參養(yǎng)殖活動，使得養(yǎng)殖池面積大大增加。大連市GDP從2000年的1 110.8億元增加到2013年的7 650.8億元，港口吞吐量從2000年的9 699萬t增加到2013年的4.1億t，人口從2000年的551.5萬人增加到2013年的591.4萬人。臨港工業(yè)的建設、港口物流業(yè)的發(fā)展以及人口的迅速增加都將導致沿海地區(qū)土地資源的稀缺，使得填海造地活動大規(guī)模進行，進而導致填海面積的不斷增加。

由此可見，近年來大連圍填海的動態(tài)演變主要由經(jīng)濟因素和建設用地需求的驅動，其中養(yǎng)殖業(yè)的發(fā)展、工業(yè)的建設、港口物流業(yè)的發(fā)展是使海域范圍向海擴張的主要動因。

5 結論

本研究利用2000—2014年LandsatTM/OLI影像數(shù)據(jù)，參考其他輔助資料，基于遙感和GIS技術，利用目視解譯法，對圍填海區(qū)域的用地類型進行了解譯，獲取了2000—2014年大連圍填海時空變化情況，統(tǒng)計了各個類型的面積分布情況，并分析了圍填海動態(tài)演變的主要驅動因素。主要結論有：海岸線變化方面，近年來大連海岸線長度總體呈增長趨勢，共增加了222.1 km；2005—2010年間是大連海岸線長度增長相當于前后兩個時期增加量的總和，由大到小依次呈現(xiàn)出2005—2010年、2000—2005年、2010—2014年的時間分布特征；各岸段海岸線長度年均變化強度由大到小依次呈現(xiàn)普蘭店、莊河、瓦房店、市轄區(qū)的空間分布特征。在圍填海面積變化方面，近年來大連圍填海面積共增加了387.07 km2，圍填?；顒忧皟蓚€時期呈快速增長的趨勢，后一時期有所減慢；2005—2010年為圍填海面積迅速增加時期。圍填海活動主要分布在瓦房店和市轄區(qū)岸段；用地類型方面，大連圍填海的主要用地類型為養(yǎng)殖用海、工業(yè)用地和港口用地；大連圍填海的動態(tài)演變主要由經(jīng)濟因素和社會因素驅動；養(yǎng)殖業(yè)的興起，工業(yè)的發(fā)展，港口的建設導致圍填海工程大量進行是驅動海域面積變化的直接因素。本研究僅對大連近年來海岸線以及圍填海的動態(tài)演變情況做了初步的分析，如何結合海洋功能區(qū)劃、主體功能區(qū)劃分結果，生態(tài)保護區(qū)劃分等結果進行適宜性評價，明確圍填海區(qū)域用地類型的合理性，提出具體的改進建議還有待于進一步深入研究。

[1] 高志強，劉向陽，寧吉才，等.基于遙感的近 30年中國海岸線和圍填海面積變化及成因分析[J].農(nóng)業(yè)工程學報，2014，12(30):141.

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

[3] 陶明剛.Landsat-TM遙感影像岸線變遷解譯研究：以九龍江河口地區(qū)為例[J].水文地質工程地質，2006(1):107-110.

[4] DONOGHUE D N，THOMAS D C，ZONG Y.Mapping and Monitoring the intertidal zone of the east coast of England using remote-sensing techniques and a coastal monitoring GIS[J].Marine Technology Society Journal，1994，28(2):19-29.

[5] SESLI F A，KARSLI F，COLKESEN I.et a1.Monitoring the changing position of coastlines using aerial and satellite image data：An example from the eastem coast of Trabzon，Turkey[J].Environmental Monitoring and Assessment，2009，15(3):391-403.

[6] YAGOUB M M，KOLAN G R. Monitoring coastal zone land use and land cover changes of Abu Dhabi using remote sensing[J].Pho-tonirvachak-Journal of the Indian Society of Remote Sensing，2006，34(1):57-68.

[7] KEVIN W，HESHAM M.Monitoring changing position of coastlines using Thematic Mapper imagery， an example from the Nile Delta[J]. Geomorphology，1999，29：93-105.

[8] RYU J H，WON J S，MIN K D.Waterline extraction from Landsat TM data in a tidal flat：a case study in Gomso Bay， Korea[J]. Remote Sensing of Environment，2002，83(3):442-456.

[9] 朱高儒，許學工.渤海灣西北岸1974—2010年逐年填海造陸進程分析[J].地理科學，2009，32(8):1006-1012.

[10] 柯麗娜，王權明.基于RS的遼寧省海岸線1990-2005年動態(tài)變化及驅動力分析.海洋開發(fā)與管理[J].2012， 29(7):54-56.

[11] 王琳，徐涵秋，李勝.廈門島及其鄰域海岸線變化的遙感動態(tài)監(jiān)測[J].遙感技術與應用，2005，20(4):404-409.

[12] 于海波，莫多聞，吳健生.深圳填海造地動態(tài)變化及其驅動因素分析[J].地理科學進展，2009，28(4):584-590.

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

教育部人文社會科學重點研究基地2009年度重大項目“基于GIS技術的遼寧沿海經(jīng)濟帶產(chǎn)業(yè)布局調整與物流環(huán)境優(yōu)化協(xié)調研究”(09JJD790058).

K902；P74

A

1005-9857(2015)08-0043-05