1985-2017年環(huán)渤海地區(qū)圍填海演化及驅(qū)動力分析

溫馨燃， 王建國,， 王雨婷， 劉志川， 王明常

(1.吉林大學(xué) 地球科學(xué)學(xué)院, 吉林 長春 130021; 2.吉林大學(xué) 地球探測科學(xué)與技術(shù)學(xué)院, 吉林 長春 130026)

沿海地區(qū)因自然環(huán)境優(yōu)越、海陸資源豐富、信息交流便捷而成為世界上人口密度最大、經(jīng)濟(jì)發(fā)展最快的區(qū)域。沿海國家為緩解人口壓力、開拓發(fā)展空間，紛紛采取圍填海的方式來開發(fā)和利用海洋。圍填?；顒又饕植加跂|南亞沿岸、波斯灣沿岸、歐洲和美洲沿岸。縱觀圍填海歷程，荷蘭、日本、韓國等國均經(jīng)歷了養(yǎng)殖圍?！ㄔO(shè)填?！獓詈９芸氐陌l(fā)展路線。人口和經(jīng)濟(jì)的快速發(fā)展促使用海方式由養(yǎng)殖轉(zhuǎn)向建設(shè)，而生態(tài)環(huán)境的持續(xù)惡化則呼吁著人們科學(xué)用海，對圍填海加以管控從而實(shí)現(xiàn)人海和諧發(fā)展。20世紀(jì)80年代以來，我國圍填海規(guī)模不斷擴(kuò)大，經(jīng)歷了由養(yǎng)殖圍海向建設(shè)填海轉(zhuǎn)變的過程，同時(shí)也加劇了海洋環(huán)境惡化、海洋災(zāi)害頻發(fā)、海洋生物多樣性銳減等一系列問題[1-2]。為協(xié)調(diào)人海矛盾，恢復(fù)沿海生境，國務(wù)院于2018年年初印發(fā)了《關(guān)于加強(qiáng)濱海濕地保護(hù)嚴(yán)格管控圍填海的通知》，嚴(yán)控新增圍填海項(xiàng)目，取消圍填海地方年度計(jì)劃指標(biāo)，除國家重大戰(zhàn)略項(xiàng)目外，全面停止新增圍填海項(xiàng)目審批。渤海做為繼珠三角、長三角后又一圍填海集中區(qū)，該地區(qū)日益嚴(yán)峻的生態(tài)、環(huán)境、資源問題引起了學(xué)界的普遍關(guān)注。基于此，國內(nèi)外學(xué)者對環(huán)渤海進(jìn)行了大量研究：從研究海岸線的定義與提取方法[3-6]，分析海岸線演化特征[7-10]，到建立海岸帶土地利用分類體系[11-12]，探討海岸帶演化特征[13-17]，再到分析、預(yù)測土地利用變化的驅(qū)動機(jī)制和響應(yīng)機(jī)制[1-2,18-20]，由線及面的探討了沿海地區(qū)的演化歷程。不難發(fā)現(xiàn)，現(xiàn)有研究多是圍繞海岸線和海岸帶展開的，缺乏以圍填?；顒訛楹诵牡尼槍π匝芯?，更缺乏在國家政策限制的條件下對圍填海發(fā)展現(xiàn)狀、演化特征、驅(qū)動因素、發(fā)展水平及發(fā)展?jié)摿Φ纳钊胙芯?。為此，本文?985—2017年環(huán)渤海圍填海區(qū)域?yàn)檠芯繉ο?，以圍填海岸線擺動區(qū)為研究區(qū)，采用基于規(guī)則的面向?qū)ο蠓诸惙ㄌ崛詈Ｍ恋乩脭?shù)據(jù)，結(jié)合分形維數(shù)、圍墾強(qiáng)度系數(shù)、開發(fā)利用強(qiáng)度和轉(zhuǎn)移矩陣定量分析渤海圍填海時(shí)空演化特征，采用灰色關(guān)聯(lián)度法定量分析社會經(jīng)濟(jì)因子對圍填?；顒拥尿?qū)動作用，此外，針對國家政策和自然條件等驅(qū)動因素進(jìn)行定性討論，旨在為加強(qiáng)圍填海綜合管治，實(shí)現(xiàn)區(qū)域可持續(xù)發(fā)展提供一定的參考依據(jù)。

1 研究區(qū)概況

研究區(qū)北起遼東半島的老鐵山角，南至山東半島北岸的蓬萊邊界，沿線自北向南分布有大連、營口、盤錦、錦州、葫蘆島、秦皇島、唐山、天津、滄州、濱州、東營、濰坊、煙臺13個(gè)地級市，分別隸屬于遼寧、河北、天津和山東4個(gè)省及直轄市[17]。環(huán)渤海地區(qū)人口承載力強(qiáng)[21]，全國海洋經(jīng)濟(jì)產(chǎn)值占比高，是中國北部沿海的黃金海岸，在中國對外開放的沿海發(fā)展戰(zhàn)略中占有重要地位。

2 數(shù)據(jù)與方法

2.1 數(shù)據(jù)來源與預(yù)處理

研究選取1985，1990，1995，2000，2005，2010年(Landsat 4-5TMC1Level-1)和2015，2017年(Landsat 8OLI/TIRSC1Level-1)共8個(gè)時(shí)相(7個(gè)時(shí)期)，56景Landsat系列衛(wèi)星遙感影像為數(shù)據(jù)源。數(shù)據(jù)下載于美國地質(zhì)調(diào)查局網(wǎng)站(USGS)(http:∥glovis.usgs.gov/)，并利用ENVI軟件進(jìn)行影像的預(yù)處理。

2.2 圍填海信息提取

2.2.1 圍填海岸線提取 研究中將海岸線分為自然型岸線和圍填型岸線并建立對應(yīng)的解譯標(biāo)志。結(jié)合岸線類型，分別提取植被線[22]、水陸分界線[23]、圍海邊界線作為圍填海岸線。當(dāng)同時(shí)存在自然型岸線和圍填型岸線時(shí)，優(yōu)先提取圍填型岸線。以1985年為基準(zhǔn)，遙感解譯各時(shí)相的圍填海岸線[8]；再將各時(shí)相圍填海岸線進(jìn)行疊加[24]，確定圍填研究區(qū)。

2.2.2 圍填海分類提取 結(jié)合地物間的光譜、紋理、空間屬性差異，采用基于規(guī)則的面向?qū)ο蠓诸惙?，將圍填研究區(qū)分為植被用地、自然海域、未利用地、養(yǎng)殖水域、鹽田用地、建設(shè)用地6類。其中，植被用地是指草地、林地等有植被覆蓋的地塊；自然海域是指未經(jīng)圍填利用的自然海面；未利用地是指沿海灘涂以及有明顯圈圍邊界但尚未開發(fā)的地塊；養(yǎng)殖水域是指用于養(yǎng)殖的圈圍水面；鹽田用地是指用于鹽業(yè)生產(chǎn)的地塊；建設(shè)用地是指用于港口、交通、住宅、產(chǎn)業(yè)園區(qū)、公共基礎(chǔ)設(shè)施等建設(shè)項(xiàng)目的地塊。結(jié)合Google Earth衛(wèi)星影像對分類結(jié)果進(jìn)行精度驗(yàn)證，分類結(jié)果均符合要求。

2.3 圍填海指數(shù)計(jì)算

2.3.1 圍填海岸線分形維數(shù) 采用網(wǎng)格法[7]計(jì)算圍填海岸線的分形維數(shù)(fractal dimension，D)。分別取邊長為r(r=30，60，90，120，150，180，210，240，270，300 m)的網(wǎng)格覆蓋圍填海岸線，覆蓋所需連續(xù)網(wǎng)格數(shù)為N(r)，利用公式(2)采取最小二乘法對分形維數(shù)進(jìn)行擬合分析。

根據(jù)分形理論：

N(r)∝r-D

(1)

對公式(1)左右同時(shí)取對數(shù)則有：

lnN(r)=-Dlnr+C

(2)

式中：C為待定常數(shù)；D為圍填海岸線的分形維數(shù)。分形維數(shù)與岸線的復(fù)雜曲折程度成正比，即D值越大，岸線越曲折。

2.3.2 圍墾強(qiáng)度系數(shù) 本研究中將吳文挺[25]等人提出的圍墾強(qiáng)度系數(shù)(coastal reclamation index， CRI)中的海岸線替換為圍填海岸線。通過計(jì)算，定量分析不同時(shí)期、不同區(qū)劃內(nèi)的岸線圍墾強(qiáng)度。

計(jì)算公式為：

(3)

式中：CRI為圍墾強(qiáng)度系數(shù)；St-t0為t至t0期間的圍填海面積；t為目標(biāo)年份；t0為起始年份；lt為目標(biāo)年份的圍填海岸線長度；lt0為起始年份的圍填海岸線長度。

2.3.3 開發(fā)強(qiáng)度指數(shù) 將各時(shí)期年均圍填開發(fā)面積與圍填研究區(qū)面積的比值定義為開發(fā)強(qiáng)度指數(shù)(development intensity index， DII)，反映人類干預(yù)下各個(gè)時(shí)期及不同用地類型間的圍填開發(fā)強(qiáng)度。

計(jì)算公式為：

(4)

(5)

(6)

DII=DIIA+DIIS+DIIC

(7)

式中：DIIA為養(yǎng)殖水域開發(fā)強(qiáng)度指數(shù)； DIIS為鹽田用地開發(fā)強(qiáng)度指數(shù)； DIIC為建設(shè)用地開發(fā)強(qiáng)度指數(shù)； DII為總開發(fā)強(qiáng)度指數(shù)；SAt-t0為t至t0期間的養(yǎng)殖水域面積；SSt-t0為t至t0期間的鹽田用地面積；SCt-t0為t至t0期間的建設(shè)用地面積；t為目標(biāo)年份，t0為起始年份；SR為圍填研究區(qū)面積。

2.3.4 灰色關(guān)聯(lián)度分析 采用灰色關(guān)聯(lián)度法[26]計(jì)算社會經(jīng)濟(jì)因子與圍填海利用類型間的關(guān)聯(lián)度，以量化社會經(jīng)濟(jì)因素的驅(qū)動作用。采用均值化歸一法實(shí)現(xiàn)社會經(jīng)濟(jì)因子的無量綱化；基于Mathlab軟件計(jì)算灰色關(guān)聯(lián)度系數(shù)及關(guān)聯(lián)系數(shù)均值，形成關(guān)聯(lián)序。

原理及計(jì)算公式為：

均值化歸一：

(8)

關(guān)聯(lián)度系數(shù)：

(9)

式中：ρ為分辨系數(shù)，取值為0.5。

3 結(jié)果與分析

3.1 圍填海岸線演化

1985—2017年，環(huán)渤海圍填海岸線共增長1 734 km，增長幅度為62%。各時(shí)期的年均增長速率呈降低—升高—再降低的趨勢。1990—1995年，圍填海岸線增長77 km，同1990年相比年均增長速率僅為0.5%；2005—2010年，圍填海岸線增長400 km，同2005年相比年均增長速率為2.3%，2010—2015年，圍填海岸線增長606 km，同2010年相比年均增長速率達(dá)3.2%，10 a間岸線增長量占總增長量的58%；其余各時(shí)期的年均增長速率在1.0%～1.5%之間。從空間分布上看，山東省岸線增長最多，為616 km，其次是河北省增長471 km，遼寧省增長396 km，天津市最少，增長251 km。但從增長幅度上看，同1985年比較，天津市增幅高達(dá)196.9%，即2017年天津市岸線長度約為1985年的3倍，其次是河北省增幅134.9%，山東省增幅70%，與環(huán)渤海整體基本持平，遼寧省增幅最低，為27.5%。研究期內(nèi)，圍填海岸線分形維數(shù)整體呈上升趨勢，從1985年的1.171 0增長到2017年的1.225 5，在2005年出現(xiàn)了極低值1.060 8。但實(shí)際上，2005年岸線長度并沒有減少，反而增長了192.6 km，說明分形維數(shù)的下降并不是由海岸線縮短導(dǎo)致的。對比2000年與2005年的岸線發(fā)現(xiàn)，此期間內(nèi)，大規(guī)模自然型岸線以“截彎取直”的方式轉(zhuǎn)為圍填型岸線，岸線趨向平直化，說明圍填海利用是導(dǎo)致分形維數(shù)大幅度下降的原因。這一現(xiàn)象在濱州至濰坊一帶表現(xiàn)的較為顯著。

3.2 岸線圍墾強(qiáng)度演化

1985—2017年，岸線平均圍墾強(qiáng)度為3.6 hm2/(a·km)(表1)，總體上呈現(xiàn)出強(qiáng)—弱—強(qiáng)—弱的變化趨勢；其中，1985—1990年和2005—2010年兩個(gè)時(shí)期圍墾強(qiáng)度高于平均值，分別為6.9，6.2 hm2/(a·km)。三省一市的圍墾強(qiáng)度變化趨勢與渤海的變化趨勢基本一致，第一個(gè)圍墾峰值均出現(xiàn)在1985—1990年；第二個(gè)圍墾峰值在時(shí)間上的表現(xiàn)略有差異，河北省、天津市的圍墾峰值出現(xiàn)在2005—2010年，山東省圍墾峰值略早于其他省市，出現(xiàn)在1995—2005年，而遼寧省圍墾峰值則略晚于其他省市，出現(xiàn)在2010—2015年(表1)。

從空間上看，第一個(gè)圍填高峰期(1985—1990年)內(nèi)，圍填?；顒又饕性诓澈５臇|北方向和西南方向，沿線13個(gè)市中有7個(gè)市的圍墾強(qiáng)度大于區(qū)域平均值，東營、盤錦、濰坊的圍墾強(qiáng)度約為平均值的3倍。第二個(gè)圍填高峰期(2005—2015年)內(nèi)，前期(2005—2010年)圍填?；顒又饕l(fā)生在渤海西部，以唐山、天津、濰坊為主要代表城市；后期(2010—2015年)圍填?；顒影l(fā)生在遼寧省、山東省各市，錦州、盤錦、濰坊為主要代表城市，而在津冀一帶，僅滄州圍墾強(qiáng)度略高于區(qū)域平均值(表1)。

表1 研究區(qū)圍墾強(qiáng)度系數(shù)變化 hm2/(a·km)

3.3 圍填海利用演化

1985—2017年，環(huán)渤海圍填海總面積4.25×105hm2，各時(shí)期內(nèi)圍填面積依次為9.07×104hm2，4.21×104hm2，4.87×104hm2，6.01×104hm2，1.05×105hm2，7.03×104hm2，8 493 hm2，面積變化趨勢與圍墾強(qiáng)度變化趨勢一致。區(qū)域內(nèi)自然海域、未利用地、植被用地分別減少了1.68×105hm2，1.42×105hm2，4.93×104hm2，養(yǎng)殖水域、建設(shè)用地、鹽田用地分別增加了1.82×105hm2，1.29×105hm2，4.80×104hm2，區(qū)域開發(fā)強(qiáng)度達(dá)84.6%，剩余可開發(fā)面積約6.55×104hm2。從開發(fā)強(qiáng)度指數(shù)來看(圖1)，1985—1990年、2000—2005年是養(yǎng)殖圍海的高峰期；1990—1995年、2005—2010年是鹽田圍海的高峰期；2005—2010年是建設(shè)用地圍填海的高峰期。由此可以看出，大規(guī)模的圍海養(yǎng)殖是第一個(gè)圍填高峰期(1985—1990年)產(chǎn)生的主要原因；而城市擴(kuò)張、港口建設(shè)等則是第二個(gè)圍填高峰期(2005—2015年)產(chǎn)生的主要原因。第二個(gè)圍填高峰期的開發(fā)強(qiáng)度遠(yuǎn)大于第一個(gè)圍填高峰期，即建設(shè)用地的開發(fā)規(guī)模與強(qiáng)度遠(yuǎn)高于養(yǎng)殖及鹽田用地。

圖1 研究區(qū)不同時(shí)期開發(fā)強(qiáng)度指數(shù)變化

通過ArcGIS進(jìn)行疊加分析獲得全時(shí)期(1985—2017年)以及兩個(gè)圍填高峰期(1985—1990年、2005—2015年)的土地利用轉(zhuǎn)移矩陣(表2)。從整體上看，1985—2017年，區(qū)域內(nèi)的養(yǎng)殖水域、建設(shè)用地和鹽田用地表現(xiàn)為完全轉(zhuǎn)入，占總面積的比例依次為42.9%，30.4%，11.3%。與之相對的，植被用地、自然海域和未利用地的轉(zhuǎn)出遠(yuǎn)大于轉(zhuǎn)入，其中植被用地占比下降11.6%，轉(zhuǎn)出方式主要為養(yǎng)殖水域和鹽田用地；自然海域占比由41.3%下降至1.8%，轉(zhuǎn)為建設(shè)用地和養(yǎng)殖水域的比例高達(dá)87.6%；未利用地占比下降33.5%，養(yǎng)殖和曬鹽是其主要的利用轉(zhuǎn)移方式。從空間分布上看，養(yǎng)殖水域主要分布在大連、錦州、唐山、濱州、東營各市；建設(shè)用地廣泛分布于營口、盤錦、錦州、唐山、天津、滄州、濰坊、煙臺各市；鹽田用地集中分布在山東省的濱州、東營、濰坊各市。第一個(gè)圍填高峰期內(nèi)呈現(xiàn)出陸域向水域轉(zhuǎn)出的趨勢，植被用地和未利用地占比下降16.9%，而養(yǎng)殖水域和自然海域占比增加16.8%。但地類轉(zhuǎn)換并不是簡單的單向轉(zhuǎn)換，而是復(fù)雜的交互轉(zhuǎn)換，植被用地轉(zhuǎn)出為未利用地、養(yǎng)殖水域和自然海域的比例依次為36.8%，17.2%，14.2%；未利用地轉(zhuǎn)出為自然海域、養(yǎng)殖水域、植被用地的比例依次為21.9%，13.3%，10.2%。該時(shí)期內(nèi)建設(shè)用地和鹽田用地的開發(fā)占比較小。

第二個(gè)圍填高峰期內(nèi)自然海域和未利用地大幅減少，建設(shè)用地、養(yǎng)殖水域、鹽田用地大幅增加。建設(shè)用地的轉(zhuǎn)入途徑主要為自然海域，占比64.7%，其次為未利用地和養(yǎng)殖水域，占比分別為13.8%，10.0%；養(yǎng)殖水域轉(zhuǎn)入的主要途徑為未利用地和自然海域，占比為25.3%和18.4%；鹽田用地的主要來源為未利用地和養(yǎng)殖用地，總占比達(dá)54.7%。截止2015年，圍填研究區(qū)內(nèi)剩余可開發(fā)面積為7.27×104hm2，占比為17.1%。對比各階段土地利用變化矩陣發(fā)現(xiàn)，植被用地的轉(zhuǎn)出方式更傾向于養(yǎng)殖水域和未利用地，自然海域則更傾向于轉(zhuǎn)換為建設(shè)用地和養(yǎng)殖水域，未利用地是養(yǎng)殖池、鹽田的主要轉(zhuǎn)入方式。

表2 研究區(qū)各時(shí)期土地利用轉(zhuǎn)換矩陣 %

3.4 圍填海驅(qū)動力分析

3.4.1 自然基礎(chǔ) 適宜的自然環(huán)境是環(huán)渤海大規(guī)模開展圍填?；顒拥幕?。環(huán)渤海沿岸分布有遼河、凌河、灤河、黃河等入海河流，以及遼東灣、渤海灣、萊州灣等海灣，入海河流帶來了大量泥沙，海灣內(nèi)水體平靜，易于泥沙沉積。在河流與海灣的共同作用下，渤海沿岸灘涂發(fā)育，營養(yǎng)物質(zhì)豐富，為錦州、唐山、濱州、東營等地區(qū)海產(chǎn)品養(yǎng)殖業(yè)的發(fā)展提供了基礎(chǔ)。海灣以其風(fēng)浪擾動小、近岸海水淺等特點(diǎn)成為城市建設(shè)的首選之地。此外，海平面的持續(xù)上升使山東沿岸海水入侵范圍不斷擴(kuò)大、土壤鹽漬化程度升高，為東營、濰坊、濱州等地的鹽業(yè)發(fā)展創(chuàng)造了條件。

3.4.2 政策驅(qū)動 21世紀(jì)以來，國家對海洋經(jīng)濟(jì)的發(fā)展日益重視，相繼批復(fù)了《遼寧沿海地區(qū)發(fā)展規(guī)劃》《天津市濱海新區(qū)行政體制改革方案》《河北沿海地區(qū)發(fā)展規(guī)劃》《黃河三角洲高效生態(tài)經(jīng)濟(jì)區(qū)發(fā)展規(guī)劃》《山東半島藍(lán)色經(jīng)濟(jì)區(qū)發(fā)展規(guī)劃》等重大決策和戰(zhàn)略部署。在此基礎(chǔ)上通過了如大連長興島臨港工業(yè)區(qū)、天津南港工業(yè)區(qū)、曹妃甸經(jīng)濟(jì)區(qū)、濰坊濱海生態(tài)旅游度假區(qū)等多項(xiàng)建設(shè)用海規(guī)劃。截止至2012年10月環(huán)渤海地區(qū)共批復(fù)建設(shè)用海規(guī)劃項(xiàng)目27個(gè)，批準(zhǔn)填海面積7.00×104hm2[21]。這表明國家政策的頒布與實(shí)施是推動環(huán)渤海圍填?；顒影l(fā)展的重要?jiǎng)恿Α?/p>

3.4.3 社會經(jīng)濟(jì)驅(qū)動 社會經(jīng)濟(jì)因素是圍填海活動的主要驅(qū)動力。本文選取2000，2005，2010，2015和2017年5期GDP，第一產(chǎn)業(yè)GDP，第二產(chǎn)業(yè)GDP，第三產(chǎn)業(yè)GDP，人口，固定資產(chǎn)投資，房地產(chǎn)開發(fā)投資，旅游人數(shù)，旅游收入，沿海港口貨物吞吐量，水產(chǎn)品總產(chǎn)量，漁業(yè)產(chǎn)值等12個(gè)驅(qū)動因子，計(jì)算了它們與6種圍填海利用類型的灰色關(guān)聯(lián)度系數(shù)，結(jié)果詳見表3。

表3 三省一市圍填海利用類型與驅(qū)動因子的關(guān)聯(lián)系數(shù)

注：X1為GDP(億元)；X2為第一產(chǎn)業(yè)GDP(億元)；X3為第二產(chǎn)業(yè)GDP(億元)；X4為第三產(chǎn)業(yè)GDP(億元)；X5為人口(萬人)；X6為固定資產(chǎn)投資(億元)；X7為房地產(chǎn)開發(fā)投資(億元)；X8為旅游人數(shù)(萬人)；X9為旅游收入(億元)；X10為沿海港口貨物吞吐量(萬t)；X11為水產(chǎn)品總產(chǎn)量(萬t)；X12為漁業(yè)產(chǎn)值(億元)。

在環(huán)渤海的三省一市中，植被用地、自然海域、未利用地都與水產(chǎn)品產(chǎn)量、區(qū)域人口、漁業(yè)產(chǎn)值間表現(xiàn)出較高的關(guān)聯(lián)度，即水產(chǎn)品產(chǎn)量的增加、區(qū)域人口的增長、漁業(yè)產(chǎn)值的提高會導(dǎo)致植被用地、自然海域和未利用地的減少。各省市的養(yǎng)殖用地驅(qū)動因子間存在差異，但差異較小，各省市養(yǎng)殖用地與水產(chǎn)品產(chǎn)量、第一產(chǎn)業(yè)GDP的關(guān)聯(lián)度都很高，此外，河北、天津地區(qū)的養(yǎng)殖用地與區(qū)域人口關(guān)聯(lián)偏高，遼寧地區(qū)的養(yǎng)殖用地與GDP關(guān)聯(lián)偏高，而山東的養(yǎng)殖用地與第二產(chǎn)業(yè)GDP關(guān)聯(lián)偏高，與植被、海域、未利用地不同，養(yǎng)殖用地的驅(qū)動因子起正向推動作用，即水產(chǎn)品產(chǎn)量、第一產(chǎn)業(yè)GDP的增加會擴(kuò)大養(yǎng)殖用地的規(guī)模。各省市建設(shè)用地驅(qū)動因子間差異較大，遼寧省建設(shè)用地受旅游收入、沿海港口貨物吞吐量和第三產(chǎn)業(yè)GDP影響較大，河北省建設(shè)用地受沿海港口貨物吞吐量、旅游人數(shù)和固定資產(chǎn)投資影響較大，天津建設(shè)用地與固定資產(chǎn)投資、第二產(chǎn)業(yè)GDP和GDP關(guān)聯(lián)較高，山東省建設(shè)用地與旅游收入、固定資產(chǎn)投資及旅游人數(shù)關(guān)聯(lián)較高。總的來看，旅游人數(shù)、固定資產(chǎn)投資與沿海港口貨物吞吐量是影響環(huán)渤海建設(shè)用地增加的主要驅(qū)動因子。各省市鹽田用地間差異最大，沒有表現(xiàn)出明顯的一致性，考慮到鹽田用地在各省市間的分布情況，僅將山東省情況加以說明。山東省鹽田用地受房地產(chǎn)開發(fā)投資、第一產(chǎn)業(yè)GDP及第三產(chǎn)業(yè)GDP影響較大。

4 討論與結(jié)論

(1) 1985—2017年，環(huán)渤海圍填?？偯娣e4.25×105hm2，年均開發(fā)強(qiáng)度指數(shù)為2.6%，圍填海岸線增長1 734 km，增長幅度為62%。各時(shí)期的年均增長速率呈降低—升高—再降低的趨勢，南部開發(fā)強(qiáng)度遠(yuǎn)高于北部，西部高于東部。

(2) 圍墾強(qiáng)度隨時(shí)間呈現(xiàn)出強(qiáng)—弱—強(qiáng)—弱的發(fā)展趨勢，兩次峰值分別出現(xiàn)在1985—1990年和2005—2015年。第一個(gè)圍填高峰以圍海養(yǎng)殖為主要利用方式，植被、未利用地大量轉(zhuǎn)為養(yǎng)殖水域；第二個(gè)圍填高峰以城市和港口建設(shè)為主，主要由自然海域轉(zhuǎn)入。

(3) 大規(guī)模的圍填?；顒邮亲匀粭l件、國家政策和社會經(jīng)濟(jì)共同作用的結(jié)果。社會經(jīng)濟(jì)發(fā)展是圍填?；顒拥闹饕?qū)動力，其中：水產(chǎn)品產(chǎn)量、區(qū)域人口、漁業(yè)產(chǎn)值是主要驅(qū)動因子。

海水養(yǎng)殖和沿海建設(shè)是近年來圍填海利用開發(fā)的主要方式，截止2017年，研究區(qū)內(nèi)42.9%的面積被用來養(yǎng)殖，30.4%的面積被用于港口和城市建設(shè)，侵占了大量濱海濕地資源，對濱海生態(tài)環(huán)境造成了一系列負(fù)面影響。解決圍填海帶來的問題勢必要從這兩方面入手：一方面對養(yǎng)殖水域進(jìn)行合理規(guī)劃，提高單位海面的產(chǎn)能，調(diào)控漁業(yè)價(jià)格，控制區(qū)域人口增速等都有利于還“養(yǎng)”為“?！钡耐七M(jìn)；一方面促進(jìn)環(huán)渤海城市群一體化，均衡區(qū)域開發(fā)強(qiáng)度，優(yōu)化建設(shè)用地布局，將沿海旅游向省市內(nèi)陸旅游引流等對恢復(fù)沿海自然環(huán)境、提高區(qū)域競爭力、實(shí)現(xiàn)可持續(xù)發(fā)展有重要意義。

本研究還存在一些不足： ①受遙感圖像精度限制，未對建設(shè)用地進(jìn)行細(xì)分，導(dǎo)致計(jì)算與社會經(jīng)濟(jì)因子的灰色關(guān)聯(lián)系數(shù)時(shí)過于籠統(tǒng)； ②2000年前后的統(tǒng)計(jì)數(shù)據(jù)差異較大，20世紀(jì)的統(tǒng)計(jì)信息不夠全面，且獲取難度大，對灰色關(guān)聯(lián)系數(shù)有一定影響； ③填海用地需要有沉降期，像金州灣機(jī)場這樣用途明確的，可直接劃分為建設(shè)用地，但用途不明的可能會被誤分到未利用地或其他地類中去，對分類結(jié)果有一定的影響。