2021年 第10期
海洋开发与管理
39
占懿娟，吴雨凝，李畅，胡业翠
(中国地质大学(北京)土地科学技术学院 北京 100089
)收稿日期:2020-10-17；修订日期:2021-09-08基金项目：国家自然科学基金项目(面上项目，重点项目，重大项目)(41877034)；中国地质大学(北京)大学生创新创业训练计划项目(202011415319).
作者简介：占懿娟，研究方向为土地资源管理
讲真的歌词是什么意思通信作者：胡业翠，教授，博士，研究方向为土地资源评价与规划
摘要：为科学开发利用海洋资源和促进粤港澳大湾区的可持续发展，文章通过目视解译遥感影像提取该区域的海岸线信息，分析海岸线的变迁强度㊁分形维数和自然岸线保有率的变化以及海岸线变迁的驱动
力㊂研究结果表明:1991 2018年粤港澳大湾区海岸线的时空变迁特征明显，总体呈向外扩张的趋势；海岸线长度增长155.76k m ，总体变迁强度为0.43%,
陆域扩张面积达482.68k m 2；海岸线平均分形维数为1.14，呈持续上升的趋势，但变化速率逐渐减缓；自然岸线保有率总体呈下降趋势，总下降幅度为24.05%，各城市呈较大差异但存在规律，可分4个类型分析海岸线变迁的驱动力；人为因素是粤港澳大湾区海岸线变迁的主要驱动因素，促使海岸线变迁的方向和速率产生巨大变化㊂
关键词：粤港澳大湾区；海岸线；自然岸线保有率；人工岸线；分形维数
中图分类号:P 748 文献标志码:A 文章编号:1005-9857(2021)10-0039-06
T h e S p a c e -t i m eC h a n g e o fC o a s t l i n e o f t h eG u a n g d o n g -H o n g K o n g -M a c a oG r e a t e rB a y A
r e aF r o m1991t o 2018Z H A N Y i j u a n ,WU Y u n i n g ,L IC h a n g
,HU Y e c u i (S c h o o l o fL a n dS c i e n c e a n dT e c h n o l o g y ,C h i n aU n i v e r s i t y o fG e o s c i e n c e s (B e i j i n g ),B e i j i n g 1
00089,C h i n a )A b s t r a c t :I no r d e r t od e v e l o p a n du t i l i z em a r i n e r e s o u r c e s s c i e n t i f i c a l l y a
n d p r o m o t e t h e s u s t a i n -a b l ed e v e l o p m e n t o f G u a n g d o n g -H o n g K o n g -M a c a oG r e a t e r B a y A r e a ,t h i s p a p
e r e x t r a c t e d c o a s t -l i n e i n
f o r m a t i o n t h r o u
g hv i s u a l i n t e r p r e t a t i o no f r e m o t e s e n s i n g i m a g e s ,a n d a n a l y z e d t
h e c h a n g
e i n t e n s i t y ,
f r a c t a l d i m e n s i o n ,n a t u r a l s h o r e l i n e p r e s e r v a t i o n r a t e o f c o a s t l i n e a n d t h e d r i v i n
g f o r c e s o f c o a s t l i n ec
h a n g e .T h er e s u l t ss h o w e dt h a tt h es h o r e l
i n eo f G u a n g d o n g -H o n g K o n g -M a c a o G r e a t e r B a y A r e a h a d o b v i o u s s p a t i o -t e m p
o r a l e v o l u t i o n c h a r a c t e r i s t i c s f r o m1991t o 2018.T h e o -v e r a l ls h o r e l i n e h a d b e e n e x p a n d i n g o u t w a r d .T h e l e n g t h o f t h e s h o r e l i n e h a d i n c r e a s e d 155.76k m ,w i t h t h eo v e r a l l c h a n g e i n t e n s i t y o f 0.43%,a n d t h e l a n de x p a n s i o na r e ah a d r e a c h e d 482.68k m 2.T h e a v e r a g e f r a c t a l d i m e n s i o nw a s 1.14,s h o w i n g a
t r e n do f c o n t i n u o u s i n c r e a s e ,b u t t h e c h a n g e r a t e g r a d u a l l y s l o w e dd o w n .T h en a t u r a l s h o r e l i n e p r e s e r v a t i o n r a t e g e n e r a l l y s港台娱乐新闻
h o w e d . All Rights Reserved.
40海洋开发与管理2021年ad o w n w a r d t r e n d,w i t h a t o t a l d e c l i n e o f24.05%.T h e r ew e r e g r e a t d i f f e r e n c e s a m o n g c i t i e s,b u t t h e r ew e r e r e g u l a r i t i e s,w h i c hc o u l db ed i v i d e d i n t o4c a t e g o r i e s t oa n a l y z e t h ed r i v i n g f o r c e so f c o a s t l i n ec h
a n g e.H u m a nf a c t o r sw e r et h e m a i nd r i v i n g f a c t o r so f t h ec o a s t l i n ec h a n g e s i nt h e G u a n g d o n g-H o n g K o n g-M a c a oG r e a t e r B a y A r e a,w h i c h g r e a t l y c h a n g e d t h e d i r e c t i o n a n d r a t e o f t h e c o a s t l i n e c h a n g e s.
K e y w o r d s:G u a n g d o n g-H o n g K o n g-M a c a o G r e a t e r B a y A r e a,C o a s t l i n e,N a t u r a l s h o r e l i n e r e t e n t i o n r a t e,A r t i f i c i a l s h o r e l i n e,F r a c t a l d i m e n s i o n
0引言
湾区经济是当今国际经济版图的亮点，具有高度开放㊁包容集聚㊁创新引领和宜居宜业等特点，是世界一流滨海城市的显著标志㊂粤港澳大湾区是我国在综合实力雄厚的珠江三角洲基础上划定的新发展区域，是我国参与全球竞争的重要空间载体[1],2019年‘粤港澳大湾区发展规划纲要“体现我国对粤港澳大湾区发展的高度重视㊂近年来，对粤港澳大湾区的研究内容主要包括气候㊁生态和环境[2-4]㊂海岸带是海洋与陆地的过渡地区以及海洋经济发展的重要空间载体，其环境变化会影响海岸线的空间形态，因此研究海岸线的变迁不仅对海岸线的开发利用和保护具有指导作用[5]，而且对海岸带环境变化的研究具有重要意义[6]㊂本研究以粤港澳大湾区为研究对象，以其海岸线变迁强度㊁海岸线分形维数和自然岸线保有率等数据为基础，科学分析粤港澳大湾区海岸线的时空变迁及其驱动力，旨在为粤港澳大湾区海岸带的开发利用和保护提供科学依据，促进粤港澳大湾区的可持续发展㊂
1研究方法
1.1研究区域
粤港澳大湾区由广州㊁深圳㊁珠海㊁佛山㊁惠州㊁东莞㊁中山㊁江门和肇庆9个珠三角城市以及我国香港和澳门2个特别行政区组成，大地构造部位大致属于华南准地台(一级单位)，呈断裂发育，故其被山地和丘陵所环绕，中心区为由西江㊁北江㊁东江和潭江等干支流冲积而成的湾内充填式三角洲，沿海地区新构造运动较为活跃[7]㊂由于填海造地项目的建设，珠江口地区大多为人工海岸，珠江口以东和珠江口以西地区以基岩海岸为主，海岸线蜿蜒曲折㊂粤港澳大湾区是我国开放程度最高和经济活力最强的区域之一，在我国经济发展格局中具有举足轻重的地位，本研究的研究区域为拥有海岸线的9个城市(除佛山和肇庆)㊂
1.2数据来源和处理
本研究以4年(1991年㊁2000年㊁2011年和2018年)的L a n d s a t遥感影像㊁相关行政区划图和G o o g l eE a r t h卫星地图为基础数据源㊂在地理空间数据云网站获取遥感影像，因研究区域范围较大，覆盖2个条带号和2个行编号(表1)㊂
表1遥感影像信息
轨道数据集
日期
(年-月-日)
分辨
率/m
条带号:121
行编号:44
L a n d s a t4-5TM1991-10-0930
L a n d s a t7E T M S L C-o n2000-06-0315
L a n d s a t4-5TM2011-01-0130
L a n d s a t8O L I_T I R S2018-03-0915
条带号:121
行编号:45
L a n d s a t4-5TM1991-10-0930
L a n d s a t7E T M S L C-o n2000-11-0115
苏见信的老婆L a n d s a t4-5TM2011-11-0130
L a n d s a t8O L I_T I R S2018-04-1015
条带号:122
行编号:44
L a n d s a t4-5TM1991-09-1430
L a n d s a t7E T M S L C-o n2000-11-0115
L a n d s a t4-5TM2011-07-0130
L a n d s a t8O L I_T I R S2018-02-1215
条带号:122
行编号:45
L a n d s a t4-5TM1991-09-1430
L a n d s a t7E T M S L C-o n2000-10-0915
L a n d s a t4-5TM2011-09-2130
L a n d s a t8O L I_T I R S2018-02-1215
在E N V I5.2软件的支持下，对遥感影像进行融合㊁纠正㊁拼接㊁补云㊁镶嵌㊁裁剪和匀等预处理，得到研究区域4个年份拼接的R G B影像㊂将研
. All Rights Reserved.
第10期占懿娟，等:1991 2018年粤港澳大湾区海岸线的时空变迁41
究区域的海岸线分为人工岸线和自然岸线2个类型，通过解译标志对遥感影像进行目视解译，得到海岸线的相关数据[8]㊂
由于数据条件的限制，参考已有研究经验[9-10]，采用目视解译方法提取瞬时水边线㊂由于无法确保该水边线即海岸线，为减少数据误差，尽量选取10-12月的遥感影像㊂利用A r c g i s10.5软件的线性构造功能，提取瞬时水边线的确切位置，同时确定海岸线的类型，得到各期海岸线长度及其变化的矢量数据㊂利用A r c g i s10.5软件的矢量数据长度计算功能，得到各期和各类型的海岸线总长度㊂
1.3研究方法
1.3.1海岸线变迁强度
海岸线长度的年均变化占比可表示某时段内各地貌岸区的海岸线变迁强度[11-12]，便于更客观地呈现各单元岸区海岸线变迁的时空差异，计算公式为:
L C I i j=L j-L i
L i(j-i)ˑ100%(1)式中:L C I i j表示某单元岸区第i年至第j年的海岸线变迁强度;L i和L j分别表示第i年和第j年的海岸线长度㊂
当L C I i j为负值时表示海岸线缩短，当L C I i j为正值时表示海岸线增长,L C I i j的绝对值与海岸线变迁强度呈正相关㊂
1.3.2海岸线分形维数
海岸线分形维数可量化表示海岸线的空间形态并展现海岸线的弯曲和复杂程度，且海岸线的弯曲和复杂程度与分形维数呈正相关㊂
基于已有研究成果[13-14]，运用网格法和M a t l a b计算研究区域各期海岸线的分形维数㊂运用A r c g i s10.5软件中的要素转栅格工具将海岸线
矢量数据转换为栅格数据，并将数据进行二值化处理，从而生成覆盖研究区域海岸线的全部正方形网格；统计不同年份的网格边长对应的网格数目[15],
对数据进行拟合得到斜率的绝对值即分形维数㊂相关表达式为:
l n N(ε)=K l nε+b(2)式中:N表示网格数目;ε表示网格边长;K表示海岸线分形维数;b表示对应常数㊂
网格的大小会对分形维数的计算产生一定的影响㊂本研究参考基本比例尺地形图数字化过程中的转换公式[16]以及比例尺为1︰90万的地图换算得到网格长度，根据地图精度粗视化最小尺度选择180m，
允许误差为0.2mm㊂由于相关资料和分形理论对各网格的取值间隔没有严格规定，根据本研究的具体情况，以最大值覆盖整条曲线为限，在网格边长序列值间隔较均匀的基础上，选取网格边长分别为180m㊁500m㊁1500m㊁3000m㊁6000m㊁9000m㊁15000m㊁21000m和50000m㊂
1.3.3自然岸线保有率
自然岸线在人类开发利用过程中被转化为人工岸线，自然岸线保有率[17]可反映人类活动对自然岸线的影响程度，计算公式为:
R=M L(3)式中:R表示自然岸线保有率;M表示自然岸线长度;L表示海岸线总长度㊂
2粤港澳大湾区海岸线的时空变迁
2.1海岸线长度和变迁强度
根据海岸线遥感监测和数据矢量化的结果, 1991 2018年粤港澳大湾区海岸线(大陆岸线)不断向海推进，空间位置变化显著，且长度增长155.76k m(表2)㊂
表21991—2018年粤港澳大湾区海岸线长度及其变化
地区
长度/k m长度变化量/k m
1991年2000年2011年2018年1991 2000年2000 2011年2011 2018年
澳门14.1415.6216.7616.761.481.140
东莞19.2919.1920.1320.13-0.100.940
广州38.9347.8448.2145.328.910.37-2.89 . All Rights Reserved.
42海洋开发与管理2021年
续表2
地区
长度/k m长度变化量/k m
1991年2000年2011年2018年1991 2000年2000 2011年2011 2018年
惠州219.25219.98217.78217.750.73-2.20-0.03江门214.12221.82223.79224.717.701.970.92深圳231.
84243.14267.15267.8311.3024.010.68香港360.78371.37369.56371.0010.59-1.811.44中山38.2155.9457.8757.8717.731.930珠海179.62212.98247.22250.5733.3634.243.35合计1316.181407.881468.471471.9491.7060.593.47
1991 2018年粤港澳大湾区的海岸线变迁强度为0.43%，其中1991 2000年较大(0.78%)，后逐渐减小,2000 2011年为0.39%,2011 2018年仅为0.03%㊂各地区的海岸线变迁强度由大到小依次为中山(1.91%)㊁珠海(1.46%)㊁澳门(0.69%)㊁广州(0.61%)㊁深圳(0.57%)㊁江门(0.18%)㊁东莞(0.16%)㊁香港(0.10%)㊁惠州(0.03%)㊂
2.2海岸线变迁引起的陆海格局变化
由于海岸线数据为非闭合的矢量数据，须利用固定的边界计算海岸线变迁造成的海岸带面积变化，进而计算岸滩面积㊂将1991年海岸线作为海岸带的内侧边界，将现有海岸线作为海岸带的外侧边界，所形成的闭合曲线就是岸滩变化区域㊂
总体来看,1991 2018年粤港澳大湾区海岸线变迁引起的陆域扩张面积为482.68k m2，侵蚀面积为16.20k m2，净增面积为466.48k m2㊂其中, 1991 2000年陆域面积变化最大，扩张面积为300.31k m2，侵蚀面积为1.74k m2，净增面积为298.57k m2，净增速为33.17k m2/a;2000 2011年陆域扩张面积为163.38k m2，侵蚀面积为6.45k m2，净增面积为156.93k m2，净增速为14.27k m2/a; 2011 2018
年陆域面积变化最小，扩张面积为17.24k m2，侵蚀面积为8.00k m2，净增面积为9.24k m2，净增速为1.32k m2/a㊂
2.3海岸线分形维数
总体来看，由于自然营力和人类活动的综合影响,1991 2018年粤港澳大湾区的海岸线分形维数平均值为1.14，且呈持续上升的趋势，但变化速率逐渐减缓，其中最大值为2018年的1.15㊂
海岸线分形维数的市级差异较大㊂①东莞㊁广州和澳门总体较小，且在1.00左右波动㊂其中：广州最小，但后期有所上升；东莞前期下降，但2000年后呈持续上升状态；澳门前期缓慢上升，但2011年由1.04降至0.94㊂②中山㊁江门和深圳为1.00~ 1.50㊂其中：中山一直上升；江门前期较小，且在上升后又持续下降；深圳在1.08左右有较小波动，近期表现为下降㊂③香港㊁惠州和珠海较大且大于平均值㊂其中：香港稳定在1.20左右；惠州有较小波动，但相对稳定在1.15左右；珠海由1.11持续上升至1.16㊂
2.4自然岸线保有率胜利床照
总体来看,1991 2018年粤港澳大湾区的自然岸线保有率呈下降趋势，总下降幅度为24.05%㊂其中,1991 2000年下降幅度最大(18.37%)，此后下降幅度骤降至3.23%,2011 2018年下降幅度有小幅回升(3.86%)㊂
漏电保护器跳闸自然岸线保有率的市级差异较大，但在差异中存在规律，主要可分为4个类型㊂①平稳且保持较高水平，自然岸线保有率均高于0.5(香港);②平稳但保持较低水平(东莞㊁广州和中山);③总体呈下降趋势，下降速率前期较高，而后期趋于0(江门和深圳);④总体呈下降趋势，下降速率前期较高，而中后期有小幅度回升(惠州和珠海)㊂值得说明的是，澳门的海岸线较短，且对现有精度的数据进行分析后认为其均为人工岸线，即自然岸线保有率为0，故不作具体讨论㊂
. All Rights Reserved.
第10期占懿娟，等:1991 2018年粤港澳大湾区海岸线的时空变迁43
3粤港澳大湾区海岸线变迁的驱动力
3.1人为因素
海岸线变迁反映土地利用强度变化，同时土地利用强度可在一定程度上用自然岸线保有率表征㊂因此，本研究在自然岸线保有率分类的基础上分析海岸线变迁驱动力的人为因素㊂
(1)平稳且保持较高水平，即香港㊂香港和东莞同处珠江口地区，自然因素的影响机制和影响程度相近，但二者的自然岸线保有率差异巨大㊂究其原因，二者的经济发展阶段不同，并体现在政策导向不同㊂香港较高的自然岸线保有率得益于其长期以来注重环境保护：一方面，避免采用低效和不可持续
的土地利用方式，提高土地利用效率，如积极发展高效的驳船运输[18]；另一方面，严格保护自然岸线，于1990年完成移沙补滩工程，使长期以来受到侵蚀的优良海滨沙滩 浅水湾重新焕发生机，并贯彻可持续发展的利用方针㊂多措并举，保持高水平的自然岸线保有率和优良的生态环境㊂
(2)平稳但保持较低水平，包括东莞㊁广州和中山㊂研究期内3个城市的自然岸线保有率均无显著变化，表明自然岸线和人工岸线保持相对平衡㊂3个城市同属珠江口地区,1995 2005年珠江口地区的围填海面积达208k m2，以渔业养殖为主，其次为建设用地，面积占比分别为75.3%和17.3%[19]㊂然而珠江口滩涂湿地土壤受重金属污染严重，不宜直接作为耕地或用于水产养殖[20]，从客观来看这也是研究期内3个城市的自然岸线保有率保持平稳的原因之一㊂2000 2010年珠江口地区海岸线向海扩张的趋势大幅减缓[21]，这种减缓的趋势一直持续到2017年[22]㊂
(3)总体呈下降趋势，下降速率前期较高，而后期趋于0，包括江门和深圳㊂在研究期前期,2个城市的经济均呈现高速发展的势头,1978 2015年江门以围垦养殖和工程建设为主的海岸线开发利用活动持续增加，截至2015年的海岸线占比分别为25.88%和21.64%[9]；深圳作为我国首个经济特区，在研究期前期处于高速发展时期,1988 1998年工程建设海岸线占比增长近1倍,1998 2008年围垦养殖海岸线占比下降且人工岸线整体呈缩短趋势,
此后人工岸线仍呈增长趋势㊂
abba
(4)总体呈下降趋势，下降速率前期较高，而中后期有小幅度回升，包括惠州和珠海㊂在研究期前期，此类数据的特征与第三类数据大致相似，其人为驱动力也以经济发展为主，但在研究期中后期的情况有所不同㊂例如：惠州的经济发展导致海岸带生态环境被严重破坏,2010年后惠州海岸带的生态安全指数大多明显下降，并于2013年达到低谷, 2013年的生态安全等级处于较危险和预警级别[23];2011年‘惠州市环境保护和生态建设 十二五 规划“发布并实施，但惠州的自然岸线保有率再次大幅下降㊂
3.2自然因素
海岸线变迁易受自然因素的影响，粤港澳大湾区海岸线变迁的自然因素主要包括3个方面㊂①地质构造运动㊂广州㊁东莞和江门处于地壳较稳定地区，地质构造运动稀疏且活动性弱；而香港和澳门地处断裂带，地质构造运动频繁，主要表现为地震和沿海地面沉降，对不同地区的海岸线产生不同程度的影响㊂②海平面上升侵蚀㊂粤港澳大湾区的地势西北高㊁东南低，且海岸带地势较低平，易受温室效应引起的冰川融化和海平面上升的影响；加上台风的影响，海水对海岸带的侵蚀强度加大，导致海岸线变迁㊂③泥沙入海㊂粤港澳大湾区拥有我国第二大河流 珠江，珠江年平均悬移质输沙量为8359万t；由于气候变化㊁龙滩截流和西江流域多级水库联合调水[24]，珠江的入海泥沙量大幅减少，对珠江口地区海岸线的影响尤为显著㊂
4结语
本研究提取1991 2018年粤港澳大湾区的海岸线长度数据，进而计算海岸线变迁强度和分形维数；通过辨识人工岸线和自然岸线，从人为因素和自然因素2个方面分析海岸线变迁的驱动力㊂粤港澳大湾区是世界第四大湾区，发展潜力不可估量，但围垦造陆等人类活动已造成海岸带生态环境出现问题，可能阻碍其发展㊂因此，分析海岸线的时空变迁可对海岸带的开发利用和保护提供科学指导㊂根据研究结果，本研究对粤港澳大湾区的发展提出建议㊂①保证自然岸线保有率总体不
. All Rights Reserved.
发布评论