基于遙感的環(huán)渤海地區(qū)海岸線變化及驅(qū)動(dòng)力分析＊

吳春生，黃翀，劉高煥，褚琳，劉慶生，趙軍

（1.中國(guó)科學(xué)院地理科學(xué)與資源研究所 資源與環(huán)境信息系統(tǒng)國(guó)家重點(diǎn)實(shí)驗(yàn)室 北京 100101；2.中國(guó)科學(xué)院大學(xué) 北京100049；3.東營(yíng)市人民政府辦公室電子政務(wù)中心 東營(yíng) 257091）

1 引言

海岸線位于大陸與海洋的交匯處，我國(guó)將其定義為“海岸線是海陸分界線，系指多年大潮高潮位時(shí)的海陸界線。”海岸線不單是一條分界曲線，同時(shí)附有一定的生態(tài)環(huán)境信息，隨著自然變化和人工影響，海岸線長(zhǎng)度和位置不斷變化，海岸線兩側(cè)海陸生態(tài)環(huán)境處于非平衡狀態(tài)，對(duì)沿海生物的生存和生活產(chǎn)生影響［1－3］。

自然條件變化和人類活動(dòng)與海岸帶生態(tài)環(huán)境之間長(zhǎng)期相互影響，引起了國(guó)內(nèi)外學(xué)者極大興趣，海岸線作為水陸分界同樣受到關(guān)注，在其類型確定和空間變化方面研究較廣。在國(guó)內(nèi)，馬小峰根據(jù)不同海岸線地貌特點(diǎn)，利用IDL語(yǔ)言實(shí)現(xiàn)了Canny算子、中值濾波器和腐蝕算子，對(duì)大連和營(yíng)口地區(qū)海岸線做了自動(dòng)提取［4］；馮蘭娣利用小波變換的圖像邊緣檢測(cè)方法對(duì)黃河三角洲海岸線進(jìn)行了提取［5］；姜義將歷史資料和遙感影像結(jié)合獲取了渤海灣西岸海岸線近百年來(lái)的變遷［6］；常軍等采用平均高潮線法來(lái)確定黃河口周邊潮灘岸線的位置［7］；同時(shí)在海岸線變化所引起的生態(tài)環(huán)境問(wèn)題方面也做了深入研究［8－12］。國(guó)外對(duì)海岸線的研究也較廣，如利用閾值法提取海岸線位置［13－14］；利用合成孔徑雷達(dá)的后向散射和紋理特征確定水陸交界線［15］；利用多時(shí)相多光譜遙感影像對(duì)海岸線變化進(jìn)行監(jiān)測(cè)［16］；利用激光探測(cè)獲取海岸DEM數(shù)據(jù)來(lái)提取海岸線［17］等。以上提取海岸線位置的技術(shù)為本文開(kāi)展研究提供了方法依據(jù)。

環(huán)渤海地區(qū)港口密集，發(fā)展?jié)摿Υ螅菍?shí)現(xiàn)我國(guó)經(jīng)濟(jì)突破的重點(diǎn)區(qū)域［18］，國(guó)家和地方也頒布了多個(gè)發(fā)展渤海海洋經(jīng)濟(jì)的條例［19］；同時(shí)環(huán)渤海沿岸具有重要的生態(tài)功能，是大量鳥(niǎo)類和魚(yú)類的棲息地，華北和東北地區(qū)南部重要的生態(tài)屏障，為區(qū)域可持續(xù)發(fā)展提供了基礎(chǔ)；由于人類活動(dòng)頻繁，生態(tài)環(huán)境脆弱，該地區(qū)生態(tài)平衡受到嚴(yán)重威脅［20］。因此對(duì)該地區(qū)海岸線變化進(jìn)行探測(cè)，及時(shí)掌握沿海開(kāi)發(fā)狀況，對(duì)合理控制規(guī)劃沿岸土地利用和生態(tài)環(huán)境保護(hù)具有重要的現(xiàn)實(shí)意義。前期對(duì)該地區(qū)集中于渤海灣和黃河三角洲海岸帶整體研究，而對(duì)渤海沿岸海岸線研究較少，對(duì)海岸線也未做具體分類，基于此，本文利用多期遙感影像對(duì)環(huán)渤海岸線進(jìn)行位置提取，并做詳細(xì)分類和分區(qū)統(tǒng)計(jì)各類型海岸線長(zhǎng)度，通過(guò)總結(jié)各類型岸線長(zhǎng)度變化強(qiáng)度和整體岸線分維數(shù)來(lái)展現(xiàn)其時(shí)空變化狀況，結(jié)合社會(huì)經(jīng)濟(jì)資料，從自然和人工因素兩方面對(duì)環(huán)渤海海岸線的變化做驅(qū)動(dòng)力分析。

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

2.1 研究區(qū)

本文以環(huán)渤海海岸作為研究區(qū)，起始端為黃海和渤海分界線兩端，即大連市老鐵山頭（38°43′34.02″N，121°07′58.92″E）和蓬萊市蓬萊角（37°49′56.27″N，120°44′37.18″E）。涉及遼寧省、天津市、河北省和山東省，相關(guān)地市有大連、營(yíng)口、盤錦、錦州、葫蘆島、秦皇島、唐山、天津、滄州、濱州、東營(yíng)、濰坊和煙臺(tái)（圖1）。

圖1 研究范圍

2.2 數(shù)據(jù)源的獲取和預(yù)處理

本文所用數(shù)據(jù)源如下。

（1）主要研究數(shù)據(jù)：研究區(qū)3期TM影像，時(shí)相為2000年、2005年和2010年6－10月，分辨率為30 m×30 m，云量小于10%，每期共7景影像覆蓋，合計(jì)21景（表1）。

表1 TM數(shù)據(jù)列表

續(xù)表

（2）輔助數(shù)據(jù)：包括2010年環(huán)境1號(hào)衛(wèi)星多光譜數(shù)據(jù)、2009－2010年CBERS02b高分辨率全色波段數(shù)據(jù)，以及部分地區(qū)2009年SPOT5和ALOS數(shù)據(jù)，主要是輔助TM影像解譯。

（3）調(diào)查數(shù)據(jù)：環(huán)渤海地區(qū)海岸線野外調(diào)查數(shù)據(jù)和統(tǒng)計(jì)數(shù)據(jù)，為海岸線和圍填海類型解譯精度評(píng)價(jià)提供支持。

對(duì)各遙感數(shù)據(jù)進(jìn)行預(yù)處理，包括幾何精校正、彩色增強(qiáng)、投影變換、圖像融合，拼接與裁剪等；統(tǒng)計(jì)數(shù)據(jù)前期處理包括數(shù)據(jù)空間化和空間統(tǒng)計(jì)分析等。

2.3 海岸線分類與遙感識(shí)別

海岸線并非現(xiàn)實(shí)存在的一條曲線，而是人們對(duì)海陸交界的一種虛擬表示，其類型則根據(jù)海岸物質(zhì)組成或用途進(jìn)行定義。海岸線類型劃分主要依據(jù)前期研究積累，結(jié)合本研究目的和研究區(qū)實(shí)際，綜合海岸線物質(zhì)組成及開(kāi)發(fā)特點(diǎn)，將海岸線分為自然岸線和人工岸線；自然岸線根據(jù)位置和物質(zhì)組成又分為河口岸線、淤泥質(zhì)岸線、砂礫質(zhì)岸線和基巖岸線；人工岸線根據(jù)開(kāi)發(fā)目的分為建設(shè)圍堤、養(yǎng)殖圍堤、農(nóng)田圍堤和鹽田圍堤；對(duì)圍填海主要根據(jù)海域圈圍目的和現(xiàn)狀進(jìn)行分類，共3類，圍海、填海和圍海養(yǎng)殖。

各類型岸線的影像特征不同，本研究依據(jù)TM2、3、4和5波段組合進(jìn)行提取。基巖岸線質(zhì)地為巖石，波浪與巖石碰撞使岸線在影像中顏色較淺，形狀多凹凸和彎曲，自然條件下多期影像無(wú)太大變化，岸線位置定于白色邊界向海一側(cè)［4，21］；河口岸線位于河道與海水相接處，根據(jù)水體顏色變化及河道寬度的改變來(lái)確定其位置，一般將河流入海口突然變寬處河口兩岸連線作為河口岸線［22］。淤泥質(zhì)岸線位于河口周邊，由河流攜帶泥沙淤積形成，岸線提取較難，本文采用常軍的平均高潮線法來(lái)確定位置［7］；砂礫質(zhì)岸線在影像上比較平滑，具有一定弧度，其紋理特征是判斷的主要依據(jù)，亮度較大［4］。人工岸線的人工信息明顯；養(yǎng)殖圍堤位于淤泥質(zhì)或沙礫質(zhì)海岸上，基巖岸線外圍也有分布，養(yǎng)殖池呈多排長(zhǎng)方形并列狀，顏色與海水相近，養(yǎng)殖圍堤平直，一般無(wú)弧度；鹽田圍堤位于淤泥質(zhì)海岸上，一般為方形且面積較大，蒸發(fā)池顏色較深，與養(yǎng)殖池主要從形狀上區(qū)分，結(jié)晶池顏色較亮，多為白色，鹽田圍堤也較平直，無(wú)弧度；農(nóng)田圍堤內(nèi)部地塊排列密集，與海域分界明顯；建設(shè)圍堤主要是港口和碼頭向海一側(cè)岸線，一般伸向海域，建設(shè)特征明顯，有船舶停靠或進(jìn)出，無(wú)固定形狀。圍海和填海根據(jù)圈圍海域是否已經(jīng)填充進(jìn)行區(qū)分，未填充為圍海，已開(kāi)始填充則為填海，而圍海和圍海養(yǎng)殖以是否具有養(yǎng)殖池作為區(qū)分依據(jù)。

首先對(duì)2010年TM影像進(jìn)行解譯，劃分各岸線類型，并參考高分影像對(duì)分類結(jié)果進(jìn)行修改。由于野外調(diào)查數(shù)據(jù)時(shí)間為2012年10月，故先將其與google earth中2010年和2012年岸線對(duì)比，提取未發(fā)生改變部分并對(duì)變化部分進(jìn)行修改使其代表2010年野外數(shù)據(jù)；然后進(jìn)行分類精度檢驗(yàn)，得出海岸線總體分類精度為89.3%，滿足本文要求。以2010年岸線類型為基準(zhǔn)，參照2005年和2000年TM影像對(duì)產(chǎn)生變化的岸線做修改，獲取2005年和2000年海岸線，避免3期影像獨(dú)立解譯產(chǎn)生的人為錯(cuò)誤。

2.4 分維數(shù)計(jì)算方法

海岸線時(shí)空變化不正體現(xiàn)在長(zhǎng)度和空間位置的變動(dòng)，其形狀結(jié)構(gòu)變化也是一種表現(xiàn)方式，作為海岸線分形理論中常用的定量指數(shù)，分維數(shù)的大小可表達(dá)海岸線形狀復(fù)雜程度，值越大，海岸線形狀越復(fù)雜［2］。本文分維數(shù)計(jì)算采用網(wǎng)格法，用不同尺度網(wǎng)格覆蓋海岸線，不同尺度對(duì)應(yīng)不同網(wǎng)格總數(shù)，取網(wǎng)格數(shù)的對(duì)數(shù)與網(wǎng)格尺度的對(duì)數(shù)的斜率的絕對(duì)值作為海岸線的分維數(shù)［23－24］。網(wǎng)格尺度的選取參照文獻(xiàn)［25］。分維數(shù)計(jì)算公式如下

式中：A為待定常數(shù)截距；ε為網(wǎng)格邊長(zhǎng)；N（ε）為邊長(zhǎng)為ε的網(wǎng)格數(shù)；－D為斜率；D即為分維數(shù)。

3 結(jié)果與分析

3.1 環(huán)渤海地區(qū)海岸線空間分布特征

根據(jù)解譯結(jié)果，統(tǒng)計(jì)各類型岸線3期長(zhǎng)度和空間分布（表2）。2010年岸線總長(zhǎng)度為3 049.7 km，自然岸線和人工岸線分別為682.3 km和2 367.4 km；基巖岸線和沙礫質(zhì)岸線為主要的自然岸線，長(zhǎng)度分別為241.0 km和234.4 km，基巖岸線集中分布于遼東半島沿岸，砂礫質(zhì)岸線在萊州灣東部、秦皇島和葫蘆島沿岸較多；其次是淤泥質(zhì)岸線，主要位于黃河口周邊，長(zhǎng)度為179.4 km；河口岸線占比較少，不到1%；建設(shè)圍堤和養(yǎng)殖圍堤分別為1 145.9 km和1 108.4 km，農(nóng)田圍堤和鹽田圍堤長(zhǎng)度較小。2005年海岸線總長(zhǎng)度2 666.8 km，自然岸線和人工岸線比例各占29.0%和71.0%；自然岸線中砂礫質(zhì)岸線最多，長(zhǎng)度為269.9 km，其次為基巖岸線，長(zhǎng)度為243.5 km，淤泥質(zhì)岸線較2010年長(zhǎng)，為234.0 km，河口岸線依然最少；建設(shè)圍堤明顯少于2010年，為635.6 km，養(yǎng)殖圍堤變化不大。2000年岸線總長(zhǎng)度為2 455.8 km，自然岸線和人工岸線分別為921.2 km和1 534.6 km；砂礫質(zhì)岸線和淤泥質(zhì)岸線較其他類型長(zhǎng)，分別為348.9 km和287.0 km，基巖岸線長(zhǎng)度為249.5 km，養(yǎng)殖圍堤最長(zhǎng)為1 111.3 km。

不同時(shí)期各省區(qū)海岸線總長(zhǎng)度如表3所示，為表達(dá)方便，將河北省和天津市海岸線合并統(tǒng)計(jì)。遼寧省海岸線長(zhǎng)度和比例始終大于其他地區(qū)，2010年各省區(qū)海岸線均達(dá)到最長(zhǎng)，遼寧省為1 404.5 km，津冀地區(qū)為792.7 km，山東省為852.5 km，分別占當(dāng)年海岸線總長(zhǎng)度的46.0%、 26.0%和28.0%。

表2 2000年、2005年和2010年不同海岸線類型長(zhǎng)度和比例

表3 2000年、2005年和2010年不同省區(qū)海岸線長(zhǎng)度和比例

3.2 近10年環(huán)渤海海岸線變化特征

3.2.1 長(zhǎng)度變化

2000年、2005年和2010年岸線逐漸增加，自然岸線減少，人工岸線增長(zhǎng)明顯。為比較清晰地對(duì)比各時(shí)段范圍內(nèi)變化差異，利用主要海岸線的年均變化百分比表示其變化強(qiáng)度（圖2）。

圖2 2000年、2005年和2010年間各類型海岸線變化強(qiáng)度

2005－2010年岸線變化強(qiáng)度高于2000－2005年，分別為2.9%和1.7%。10年間岸線共增長(zhǎng)593.9 km，增加率為24.2%，前5年和后5年分別增加211.0 km和382.9 km。人工岸線前后5年變化強(qiáng)度分別為4.7%和5.0%，10年間增加832.8 km，增加率達(dá)54.3%，前后5年分別增加357.6 km和475.2 km。人工岸線高強(qiáng)度變化依賴于建設(shè)圍堤增加，2000－2010年建設(shè)圍堤變化強(qiáng)度為20.6%，長(zhǎng)度增加771.3 km，前后5年分別增加261.1 km和510.2 km，遠(yuǎn)超其他岸線。自然岸線前后5年變化強(qiáng)度為－3.2%和－2.4%，10 年間 減 少 238.9 km，損 失 率 為25.9%，長(zhǎng)度變化最大的為砂礫質(zhì)岸線，減少114.5 km，其次為淤泥質(zhì)岸線107.6 km，其他自然岸線減少量較小。

各省區(qū)岸線長(zhǎng)度均呈逐年增長(zhǎng)趨勢(shì)，遼寧省前后5年變化強(qiáng)度分別為0.8%和3.2%，增長(zhǎng)集中于2005－2010年，增加長(zhǎng)度194.9 km；津冀地區(qū)各時(shí)間段內(nèi)海岸線變化強(qiáng)度大于其他地區(qū)，前后5年分別為4.6%和3.4%，增加集中于2000－2005年，增加長(zhǎng)度為125.8 km；山東省海岸線變化強(qiáng)度則相對(duì)較為平穩(wěn)，10年間共增加113.3 km。

3.2.2 形狀變化

對(duì)3期岸線按照不同尺度進(jìn)行網(wǎng)格分割，獲取各期岸線網(wǎng)格總數(shù)，并分省統(tǒng)計(jì)，得到2000年、2005年和2010年整體海岸線的分維數(shù)分別為1.078 3、1.086和1.103 6，分維數(shù)逐年增加，海岸線形狀越來(lái)越復(fù)雜，局部地區(qū)海岸線的彎曲度增加，分維數(shù)大小的變化也反映出對(duì)原有岸線空間和長(zhǎng)度上的改變程度（表4）。

表4 2000年、2005年和2010年環(huán)渤海地區(qū)海岸線分維數(shù)

各地區(qū)分維數(shù)與整體岸線趨勢(shì)相同，呈增長(zhǎng)態(tài)勢(shì)，岸線長(zhǎng)度和形狀復(fù)雜度越來(lái)越大；山東省岸線分維數(shù)變化較小，2010年萊州灣南部岸線向海凸出，增加了局部岸線復(fù)雜度；津冀地區(qū)海岸線變化最為劇烈，海岸線向海域延伸頻度增大，空間變廣，極大增加了海岸線彎曲迂回程度。遼寧省海岸線分維數(shù)大于山東省和津冀地區(qū)，與其自身海岸形狀有關(guān)，且遼寧海岸多基巖岸線分布，養(yǎng)殖圍堤也較多，增加了整體岸線復(fù)雜度。

3.3 海岸線變化驅(qū)動(dòng)力分析

本文從自然和人工因素對(duì)海岸線變化進(jìn)行驅(qū)動(dòng)力分析，自然因素?zé)o規(guī)律，具有偶然性，主要影響自然岸線；人工因素主要是經(jīng)濟(jì)建設(shè)活動(dòng)，可能會(huì)完全改變?cè)凶匀换蛉斯ぐ毒€的類型和空間位置。

3.3.1 自然蝕淤

淤積和蝕退是導(dǎo)致海岸線變化的主要自然因素；自然淤積是河流入海攜帶的泥沙在入海口周邊沉積形成淤泥質(zhì)岸線或沙礫質(zhì)岸線，同時(shí)改變自然岸線的位置和長(zhǎng)度；自然蝕退受海潮影響較重，遭侵蝕岸線主要為淤泥質(zhì)、砂礫質(zhì)和河口岸線，部分人工岸線也會(huì)受影響。據(jù)前人研究成果和相關(guān)部門的統(tǒng)計(jì)分析，黃渤海每年多達(dá)60次不同強(qiáng)度的風(fēng)暴潮對(duì)環(huán)渤海岸線及海岸生態(tài)環(huán)境帶來(lái)嚴(yán)重威脅，尤其是山東沿海。對(duì)比各地區(qū)岸線受自然因素影響程度，山東省最重，經(jīng)現(xiàn)場(chǎng)調(diào)查和遙感分析可知，各時(shí)間段海岸蝕退大于淤積（表5和圖3）。

表5 山東省渤海海岸蝕淤狀況 km2

圖3 山東省渤海海岸蝕淤分布

2000－2010年山東省海岸蝕退159.5 km2，淤積42.0 km2，海岸的蝕退和淤積造成海岸線在空間上向陸地后撤或向海域突進(jìn)，最長(zhǎng)距離分別為6.8 km和7.3 km，自然岸線增加113.27 km，變化的主要是淤泥質(zhì)岸線；2000－2005年間渤海灣發(fā)生強(qiáng)風(fēng)暴潮的次數(shù)增多，強(qiáng)度變大，5年間海岸蝕退117.8 km2，導(dǎo)致自然岸線后撤，而總長(zhǎng)度卻增加40.98 km；除此之外，北部刁口河河口自1976年黃河改道后，斷絕了泥沙和河水補(bǔ)給，造成海水倒灌，海岸侵蝕嚴(yán)重，現(xiàn)黃河唇形入海口南部區(qū)域遭侵蝕原因相同，而北部由于黃河繼續(xù)挾帶泥沙如海，河口出現(xiàn)淤積；2005－2010年，海洋災(zāi)害減弱，黃河口和萊州灣南部河流入海泥沙沉積增加，被侵蝕海岸出現(xiàn)淤積，淤積面積增至57.7 km2，海岸線呈向海推進(jìn)；2005年以來(lái)，國(guó)務(wù)院和黃委會(huì)商討將刁口河作為黃河備用流路，于2009年7月借助黃河調(diào)水調(diào)沙工程進(jìn)行注水試驗(yàn)，保持入海水量和泥沙量，有效阻止海水的進(jìn)一步侵蝕。渤海沿岸其他河流如海河和遼河入海處以及淤泥質(zhì)岸線分布區(qū)均存在蝕退和淤積狀況，但面積和范圍都較黃河三角洲地區(qū)小，且近些年其流路并未發(fā)生變化，所受自然因素影響明顯小于人工因素，蝕退和淤積量較少，故未做詳細(xì)分析。

3.3.2 圍填海影響

隨著人口增多，土地資源優(yōu)勢(shì)已不明顯，經(jīng)濟(jì)發(fā)展壓力及國(guó)家和地區(qū)政策的實(shí)施迫使沿海地區(qū)加大近海開(kāi)發(fā)，海岸帶人類活動(dòng)頻繁，尤其工業(yè)建設(shè)帶來(lái)的近海格局變化導(dǎo)致海岸線類型和區(qū)位發(fā)生改變，圍填海作為人類對(duì)海岸開(kāi)發(fā)的外在表現(xiàn)，能夠從定量角度展示岸線變化程度。

通過(guò)解譯，獲取研究區(qū)圍填海面積和空間分布信息（圖4），2000－2010年圍填海面積達(dá)1 558.7 km2，圍海、填海和圍海養(yǎng)殖分別為410.0 km2、656.3 km2和486.2 km2；圍填海空間上多分布于平直海岸，港口和碼頭建設(shè)用圍海和填海向海域延伸較長(zhǎng)；填海多出現(xiàn)在原有養(yǎng)殖圍堤和建設(shè)圍堤向海一側(cè)，對(duì)自然岸線影響不明顯，而圍海和圍海養(yǎng)殖在原自然和人工岸線上均有分布，尤其是對(duì)自然岸線中的砂質(zhì)岸線和淤泥質(zhì)岸線的占用，使兩者在10年間連續(xù)減少，且減少長(zhǎng)度較大，同時(shí)自然岸線大量縮短，人工岸線明顯增加，而養(yǎng)殖圍堤長(zhǎng)度有增有減，整體長(zhǎng)度未發(fā)生太大變化，建設(shè)圍堤增加771.3 km。從表1中各類型岸線解譯標(biāo)志可知，建設(shè)圍堤形狀復(fù)雜度明顯大于其他岸線，建設(shè)圍堤的大量增加導(dǎo)致岸線分維數(shù)逐步變大。

從區(qū)位上看，填海主要分布在津冀地區(qū)，圍海和圍海養(yǎng)殖則在遼寧和山東海岸較多，這與國(guó)家和地區(qū)發(fā)展規(guī)劃有關(guān)，津冀地區(qū)尤其是天津海岸港口經(jīng)濟(jì)是區(qū)域發(fā)展的重點(diǎn)，10年間新建港口有曹妃甸港、漢沽港、天津港、天津大港和黃驊港，建設(shè)圍堤增長(zhǎng)341.2 km，增加了近2.5倍。而遼寧和山東沿海則以海水養(yǎng)殖、海洋旅游以及鹽化工為主，并未把港口建設(shè)作為重點(diǎn)建設(shè)項(xiàng)目，位于研究區(qū)內(nèi)的港口規(guī)模較小；2005－2010年，因?yàn)H坊濱海生態(tài)旅游度假區(qū)建設(shè)開(kāi)工，圈圍大面積海域，造成此時(shí)間段內(nèi)山東岸段圍海面積陡增以及建設(shè)圍堤突然變長(zhǎng)，但沿海其他區(qū)域并未發(fā)生太大變化。后5年圍填海總面積比前5年多703.0 km2，填海增加372.2 km2，圍海增加347.9 km2，圍海養(yǎng)殖減少17.2 km2；圍海和填海面積增加反映海岸開(kāi)發(fā)活動(dòng)激烈程度加大，建設(shè)圍堤長(zhǎng)度后5年增加量超過(guò)前5年，分別為510.2 km和261.1 km。填海增加主要來(lái)自津冀地區(qū)港口建設(shè)，圍海在遼東灣北部和萊州灣南部增加明顯，而圍海養(yǎng)殖在整個(gè)研究區(qū)內(nèi)小幅減少。

圖4 環(huán)渤海地區(qū)圍填海面積統(tǒng)計(jì)

4 結(jié)論

本文利用遙感方法提取了不同類型海岸線和圍填海分布，從自然因素和人工因素兩方面對(duì)海岸線變化做了驅(qū)動(dòng)力分析，主要結(jié)論有：

環(huán)渤海地區(qū)海岸線逐年增長(zhǎng)，2000年、2005年和2010年總長(zhǎng)度分別為2 455.8 km、2 666.8 km和3 049.7 km，自然岸線大量減少，人工岸線顯著增加；建設(shè)圍堤變化強(qiáng)度最高；整體海岸線分維數(shù)逐年變大，海岸線形狀越來(lái)越復(fù)雜，特別是津冀地區(qū)復(fù)雜度最高。

海潮侵蝕和河流入海泥沙淤積對(duì)自然岸線空間分布影響較大，現(xiàn)黃河入海口和故道入海處淤積和侵蝕對(duì)比較為明顯，海岸線突進(jìn)和后撤最大距離分別為7.3 km和6.8 km，致使自然岸線增加。受人口增加和土地資源緊缺的壓力，各地區(qū)對(duì)海洋開(kāi)發(fā)加快，使圍填海10年間總面積變化較大，后5年開(kāi)發(fā)建設(shè)較前5年更為激烈。2000－2005年間以圍海養(yǎng)殖為主，面積為238.1 km2，填海和圍海面積較小；2005－2010年間各地區(qū)加大開(kāi)發(fā)力度，以填海類型為主，面積為501.7 km2，港口、碼頭和城鎮(zhèn)建設(shè)等已具有一定規(guī)模，圍海養(yǎng)殖面積較前5年小，圍海面積較前5年大，即未來(lái)建設(shè)活動(dòng)還將進(jìn)一步增多。

津冀地區(qū)港口和碼頭建設(shè)最頻繁，岸線和圍填海變化最大，港口面積成倍增加；遼寧岸線最長(zhǎng)，根據(jù)當(dāng)?shù)睾Ｑ蠼?jīng)濟(jì)發(fā)展規(guī)劃，人工岸線以養(yǎng)殖圍堤為主，并逐年增長(zhǎng)，部分岸段存在港口建設(shè)，但規(guī)模較小；山東海岸受岸線性質(zhì)和穩(wěn)定性因素影響，開(kāi)發(fā)力度較其他兩地區(qū)小，岸線變化主要由自然因素導(dǎo)致，用于保護(hù)油氣資源開(kāi)發(fā)而建設(shè)的堤壩較多，致使建設(shè)圍堤有所增長(zhǎng)，圍填海類型主要為圍海養(yǎng)殖。

［1］高義，王輝，蘇奮振，等 .中國(guó)大陸海岸線近30a的時(shí)空變化分析［J］.海洋學(xué)報(bào)，2013，35（6）：31－42.

［2］徐進(jìn)勇，張?jiān)鱿椋w曉麗，等.2000－2012年中國(guó)北方海岸線時(shí)空變化分析［J］.地理學(xué)報(bào)，2013，68（5）：651－660.

［3］姚曉靜，高義，杜云艷，等 .基于遙感技術(shù)的近30a海南島海岸線時(shí)空變化［J］.自然資源學(xué)報(bào)，2013，28（1）：114－125.

［4］馬小峰，趙冬至，邢小罡，等 .海岸線衛(wèi)星遙感提取方法研究［J］.海洋環(huán)境科學(xué)，2007，26（2）：185－189.

［5］馮蘭娣，孫效功，胥可輝 .利用海岸帶遙感圖像提取岸線的小波變換方法［J］.青島海洋大學(xué)學(xué)報(bào)，2002，32（5）：777－781.

［6］姜義，李建芬，康慧，等 .渤海灣西岸近百年來(lái)海岸線變遷遙感分析［J］.國(guó)土資源遙感，2003（58）：54－58.

［7］常軍，劉高煥，劉慶生 .黃河三角洲海岸線遙感動(dòng)態(tài)監(jiān)測(cè)［J］.地球信息科學(xué)，2004，6（1）：94－98.

［8］雷坤，孟偉，鄭丙輝，等 .渤海灣海岸帶生境退化診斷方法［J］.環(huán)境科學(xué)研究，2009，22（12）：1361－1365.

［9］聶紅濤，陶建華 .渤海灣海岸帶開(kāi)發(fā)對(duì)近海水環(huán)境影響分析［J］.海洋工程，2008，26（3）：44－50.

［10］秦延文，鄭丙輝，李小寶，等 .渤海灣海岸帶開(kāi)發(fā)對(duì)近岸沉積物重金屬的影響［J］.環(huán)境科學(xué)，2012，33（7）：2359－2367.

［11］索安寧，張明慧，于永海，等 .曹妃甸圍填海工程的環(huán)境影響回顧性評(píng)價(jià)［J］.中國(guó)環(huán)境監(jiān)測(cè)，2012，28（2）：105－110.

［12］趙鑫，孫群，魏皓 .圍填海工程對(duì)渤海灣風(fēng)浪場(chǎng)的影響［J］.海洋科學(xué)，2013，37（1）：7－15.

［13］CHEN L C ，RAN J Y.Detection of shoreline changes for tide-land areas using multi-temporal satellite images［J］.Remote Sensing，1998，19（17）：3383－3397.

［14］PTATT W K.Digital image processing［M］.New York：John Wiley ＆Sons Inc，2001：551－581.

［15］PALMER P L，PETROU M.Locating boundaries of textured regions［J］.IEEE Transactions on Geosicience and Remote Sensing，1997，35（5）：1367－1371.

［16］EKERCIN S.Coastline change assessment at the aegean sea coasts in turkey using multitemporal landsat imagery［J］.Journal of Coastal Research，2009，23（3）：691－698.

［17］WHITE S A，WANG Y.Utilizing dems derived from LIDAR data to analyze morphologic change in the Carolina coastline［J］.Remote Sensing of Environment，2003，85（1）：39－47.

［18］田東娜，馮筱婧，江海旭 .環(huán)渤海地區(qū)海洋旅游發(fā)展研究［J］.海洋開(kāi)發(fā)與管理，2014，31（2）：111.

［19］黃盛 .環(huán)渤海地區(qū)海洋產(chǎn)業(yè)結(jié)構(gòu)調(diào)整優(yōu)化研究［D］.青島：中國(guó)海洋大學(xué)，2013.

［20］索安寧，于永海，苗麗娟 .渤海海域生態(tài)系統(tǒng)功能服務(wù)價(jià)值評(píng)價(jià)［J］.海洋經(jīng)濟(jì)，2011，1（4）：42－47.

［21］孫偉富，馬毅，張杰，等 .不同類型海岸線遙感解譯標(biāo)志建立與提取方法研究［J］.測(cè)繪通報(bào)，2011（3）：41－44.

［22］馬小峰，趙冬至，張豐收，等 .海岸線衛(wèi)星遙感提取方法研究進(jìn)展［J］.遙感技術(shù)與應(yīng)用，2007，22（4）：575－580.

［23］Grassberger P.On efficient box counting algorithms［J］.International Journal of Modern Physics C，1983，4（3）：515－523.

［24］韓雪培，傅小毛，湯景燕，等 .圖上曲線長(zhǎng)度量算的分維糾正法［J］.華東師范大學(xué)學(xué)報(bào)，2006（6）：34－40.

［25］高義，蘇奮振，周成虎，等 .基于分形的總過(guò)大陸海岸線尺度效應(yīng)研究［J］.地理學(xué)報(bào)，2011，66（3）：331－339.