黄诚,吴能友,龙军桥,汤伟,辛卓,罗生龙,胡旋,颜历,邢景峰,陈龙,张云锺,郭泽俊,李央
1. 中国地质调查局海口海洋地质调查中心,海口 571127
2. 中国海洋大学,青岛 266100
3. 中国地质调查局青岛海洋地质研究所,青岛 266237
4. 东华理工大学,南昌 330013
传统的海岛地质调查方法如登岛人工调查、船载作业等,受地形、地貌复杂因素和开发程度影响,面临较大困难[1]。20 世纪90 年代以来,卫星遥感技术逐步在海岛地质调查工作中获得了较好的应用效果,如利用卫星遥感影像进行大比例尺海岛(礁)测绘[2-3]、利用高分辨率卫星遥感技术对海岛动态变化进行监测[4]、应用遥感影像对人类活动特征进行识别[5]等。卫星遥感技术的发展为海岛综合调查带来了极大便利,具有空间分辨率高、定位精确和地物信息丰富等优势[6],但也存在诸如历史数据质量良莠不齐(受云层和潮位影响)、现势性和分辨率不能满足高精度和高时效调查需要、只能获取正射影像等不足[7-8]。无人机低空倾斜摄影是近年来发展起来的一项新的测量技术,取得了很好的实际应用效果,展示了广泛的应用前景[9-10]。
与卫星遥感技术相比,无人机低空倾斜摄影技术在海岛综合地质调查方面优势明显[8,11-13]:①低空飞行,可以有效规避云层和潮汐影响,获取海岛高质量光学影像;②搭载多个摄像机,获取不同角度高分辨率影像,实现三维实景建模,增强信息提取的效率;③可以针对具体目标进行反复拍摄,提高重点调查精度;④操作简单,低成本,设计好航线后,可以实现自动采集相关数据;⑤可实施高风险区调查任务。对于人员无法到达的危险区域,应用无人机技术有独特优势。
本文采用低空无人机倾斜摄影技术,对三亚市蜈支洲岛开展全方位倾斜摄影测量工作,并应用Smart3D 将获取的不同角度数据进行分析处理,建立海岛三维实景模型。基于该模型,结合海岛实测资料,提取有效基础地质、环境地质、旅游地质及灾害地质信息,为海岛空间规划、开发利用和生态保护修复提供决策依据。
蜈支洲岛位于海南岛三亚东部海棠湾,是国家5A 级景区,又称牛奇洲(图1)。蜈支洲岛西北角距离海南本岛约2.5 km,北侧与椰子洲岛遥相呼应,南临亚龙湾,全岛呈不规则“爱心状”,面积约1.01 km2,东西长1 500 m,南北宽1 100 m[14-15]。蜈支洲岛中心地理坐标位置为18°18′43″N、109°45′44″E,海拔高约79.3 m,岸线总长度约6.91 km。蜈支洲岛北部发育第四系烟墩组(Qh3y)灰白色亚砂土、含砾黏土、灰黄色砾质砂及海滩岩;东南、西南两侧出露燕山早期花岗岩,主要岩性为中粗粒含斑黑云母花岗岩;中部出露燕山晚期花岗岩,主要岩性为细粒黑云母花岗岩。
图1 研究区所在地理位置Fig.1 Geographical location of the study area
蜈支洲岛地处热带,属热带海洋气候,受南部季风影响,雨量充沛,岛上拥有85 科2 700 多种原生植物,包括有“地球植物老寿星”之称龙血树等珍稀植物[16]。蜈支洲岛海洋牧场是我国第一个热带休闲旅游型海洋牧场国家级示范区,近岛珊瑚礁区域是休闲渔业开发的核心区[17]。与此同时,三亚蜈支洲岛海洋牧场靠近琼东沿岸上升流核心区,其近岸珊瑚礁位于我国珊瑚礁过渡区,上升流可以减弱区域内海水温度上升对珊瑚胁迫影响,使得岛周珊瑚礁保持了较高的生物多样性[18]。
无人机采用飞马D20 机型,搭载1.2 亿像素的D-OP3000 模块相机镜头(其中1 个为正射镜头,其余为倾斜镜头)。D20 及相机镜头参数如表1。
表1 飞马D20 无人机和D-OP3000 相机硬件参数Table 1 Hardware parameters of FEIMA D20 UAV and DOP3000 cameras
本次工作无人机倾斜摄影时间为2020 年8 月30—31 日,历时2 天,飞行期间天气晴朗,风速3 级,飞行条件良好。无人机倾斜摄影测区面积约为1.8 km2,在踏勘的基础上进行了航线规划设计,航线相关参数按照《GB/T 39612-2020 低空数字航摄与数据处理规范》执行[19]。采用国家2000 坐标系统和2000 国家高程基准,飞行高度190 m,航向216°,无人机飞行姿态俯仰角6~8°,航向重叠率70%,旁向重叠率60%,影像平均地面分辨率优于0.03 m。在蜈支洲岛中部较为空旷、通视条件较好的电动车站设置无人机起降点,自北向南依次等间距部署21 条航线(图2),为保证全面覆盖设计范围,航线做了必要外扩,本次共采集照片5 535 张。
图2 无人机倾斜摄影航线规划图Fig.2 Route planning of the UAV tilt photographing
无人机倾斜摄影工作过程中采用飞马D20 自带的RTK 模块进行实时厘米级定位,并根据地形起伏进行变高飞行,以保证影像一致性。同时,按照航摄要求,提前在蜈支洲岛遮挡较少的空旷区域或有标志性建筑的区域布置像片控制点,本次工作布设像片控制点10 个(表2),以保证影像数据质量。
表2 像控点信息表Table 2 Information at the image control points
同一地面点会在多张影像上成像,这些点的集合称为连接点,连接点生成在摄影测量中非常重要,为捆绑调整提供必要的输入信息量[20]。数据处理前,对所有采集到的5 535 张照片进行标定,标定成功的照片数量为3 160 张,获取的连接点为24.8 万个。未标定成功的照片中大部分为处于航线的外延部分。应用Smart3D 软件对标定成功的3 160 张照片进行处理,经过差分解算、处理分区和空三计算,得到蜈支洲岛及其周边部分海域的数字正射影像图(DOM)和实景三维模型(图3)。
图3 蜈支洲岛实景三维模型图Fig.3 3D real scene model of the Wuzhizhou Island
空三计算是通过无人机倾斜摄影外业测量的地面像控点位置信息,在室内无人机影像处理过程中辅助确定全部影像的位置信息元素的一种测量方法[21]。空三计算的关键取决于航测过程中获取POS(机载定位定向系统)数据的精度,POS 数据包括该组镜头航测瞬间所处位置的高程、经纬度、航向角、俯仰角及翻滚角信息[22]。本次研究中无人机倾斜摄影采用前、后、左、右和正射5 个镜头同时拍摄,5 个镜头的POS 数据同时参与计算,显著提高了计算精度,提升了构建实景三维模型的质量和精度。空三计算完成之后,Smart3D 软件可自动完成实景三维模型构建。
实景三维模型不仅可以多角度、立体化实现研究区野外场景重现,利于对研究区现状进行宏观和微观相结合的观察、理解,而且可借助无人机倾斜摄影测量成果对关注的地质点、线、面、体等重点地质信息进行识别、量测和分析[21]。具体包括地质体空间坐标信息、地质接触界面、节理及断裂产状、现状岸线侵蚀淤积状况、地质灾害及隐患区空间分布范围等。本次工作基于高精度无人机倾斜摄影数据构建的实景三维模型,结合数字正射影像图(DOM),对蜈支洲岛进行综合地质解译。
本次倾斜摄影测量的第5 组镜头为正射影像镜头,利用Smart3D 软件对正射数据经差分解算、控制点平差和空三计算等步骤形成正射影像。正射影像图可视同为近地遥感数据,兼具地图的几何特征和影像特征,以高精度、直观易读和缩放方便等优点,广泛地应用于地图信息处理、自然资源保护、灾害防护治理等领域[23]。通过对正射影像图进行解译,可以比较直观研究蜈支洲岛海岛线性构造以及与断裂密切相关的球状风化等地质现象。
通过对数字正射影像进行解译,可以清晰识别出NW、NE 和近SN 向3 组线性断裂及部分环形断裂(图4)。这些线性断裂多沿着节理发育,规模不一,其中NW 向断裂占据主导地位。NW 向断裂走向301°~333°,倾角近直立;NE 向断裂走向为40°~48°,倾角78°~86°;近SN 向断裂走向342°~353°,倾角76°~82°。沿NW 向和近SN 向断裂发育了多处基性岩脉,其宽度为0.5~2.4 m,表现为岩墙群形态,岩性为辉绿岩(图5)。环形断裂可能是蜈支洲岛上后期侵位的花岗岩体所致。
图4 蜈支洲岛线性断裂分布图基于DOM 数据。Fig.4 Distribution of linear fracture in Wuzhizhou Island Based on DOM.
图5 蜈支洲岛东南部发育的基性岩脉基于实景三维模型。Fig.5 Basic dykes developed in southeast of the Wuzhizhou Island Based on real-scene 3D Model.
蜈支洲岛是一个花岗岩型海岛,除北部被少量第四系沉积物覆盖以外,其余部分均为花岗岩。花岗岩的球状风化现象在自然界比较常见,一般认为,其成因为岩体早期受三维裂隙控制形成四方块体,由于棱角突出,易受风化,故棱角逐渐缩减,最终趋向球形[24]。除了普通的球状风化,花岗岩的球状风化还可表现为如洋葱般层层剥落,可能是花岗岩在循环失水-吸水过程中,反复不均匀收缩、膨胀促使微观岩体结构扰动和拉伸,最终形成洋葱状球状风化现象[25]。蜈支洲岛岩滩部分花岗岩普遍发育球状风化,且普通球状风化和洋葱状球状风化均有发育,特别是在海浪能够作用到的范围易形成层层剥落的球状风化。应用实景三维模型,在蜈支洲岛上能够比较快捷地发现球状风化发育的各个阶段(图6)。球状风化不仅对工程地质条件影响较大,而且是加速基岩岸线侵蚀的一个重要潜在因素,对蜈支洲岛来说是一个不良地质现象。
图6 蜈支洲岛岛岸球状风化4 个不同发育阶段Fig.6 Different development stages of spherical weathering alone the coastal Wuzhizhou Island
基于DOM 数据遥感解译,结合实地调查,对海岛岸线类型及空间分布和海岛土地利用现状图等信息进行识别。应用Arcgis 等软件还可统计各开区块面积、各类岸线长度信息等。蜈支洲岛岸线类型主要为基岩岸线、砂质岸线和人工岸线(图7)。其中基岩岸线除海岛西北角以外几乎遍布全岛,总长度约5.38 km,基岩岩性为细粒黑云母花岗岩和中粗粒含斑黑云母花岗岩。砂质岸线分布于海岛北侧及西北角,总长度约1.08 km,砂质洁白,蜈支洲岛重要的旅游资源。人工岸线主要为夏、冬两季码头修筑时所形成的地基及少部分人工构筑的防波提,分布于海岛西北角。
图7 蜈支洲岛岸线类型及空间分布(基于DOM 数据)Fig.7 Coastline types and spatial distribution of the Wuzhizhou Island
蜈支洲岛的岸滩主要为岩滩和沙滩,几乎不含砾石滩,蜈支洲岛的沙滩是世界上极少数不含或仅含少量砾石和其他岩屑的海滩[14](图8),这与位于三亚市宁远河入海口处养生园岛岛滩明显不同(图9)。反映了其不同的沉积物源供应特点。
图8 蜈支洲岛北部沙滩Fig.8 Sandy beach of the northern Wuzhizhou Island
图9 崖州湾含砾石沙滩Fig.9 Sandy beach with pebbles and cobbles in the Yazhou Bay
蜈支洲岛是海南乃至全国著名的旅游海岛,海岛的开发利用主要围绕旅游和生态保护进行。按照GB/ T21010-2017 土地利用现状分类方案[26],将蜈支洲岛土地利用状况划分为乔木林地、灌木林地、建设用地等10 类。其中乔木林地分布面积最大,为637 630 m2;其次为灌木林地,统计面积为120 478.16 m2;建设用地面积为62 294.2 m2。对各用地类型空间分布及面积进行了统计(图10),以乔木林地、灌木林地和草地总面积,计算蜈支洲岛的植被覆盖率约为77.8%。
图10 蜈支洲岛土地开发利用现状图Fig.10 Current situation of land development and utilization in the Wuzhizhou Island
海岛地质灾害是指对海岛地质、生态环境和地貌景观造成破坏,或对海岛人民生命财产造成直接损失的,分布于岛陆、岛岸、环岛近岸海域的地质现象和地质作用[27]。随着海岛开发力度和强度的提升,潜在的地质灾害也日益突出[28]。常见的海岛地质灾害主要包括滑坡、崩塌、海岸侵蚀、湿地退化、地面沉降、风沙灾害、沙滩泥化等,其中以滑坡、海岸侵蚀、海水入侵和滨海湿地退化最为典型[29-30]。按照分布位置,海岛地质灾害又可分为岛陆、岛岸和近岸海底地质灾害。总体来看,岛岸地质灾害最多,也是受灾害影响最严重区域。通常,海岛岛岸地区植被覆盖度低,通视情况较好,应用遥感技术能够取得较好的应用效果。
本次工作基于实景三维模型,获取了蜈支洲岛岸崩塌、断裂活动、海岸侵蚀等地质灾害信息,共发现崩塌风险隐患点5 处(图11-①⑤⑥⑦⑯),断裂活动风险隐患点3 处(图11-①②⑰);海岸侵蚀隐患点7 处(图11-①③④⑫⑬⑭⑮)。
图11 蜈支洲岛地质灾害隐患点分布图Fig.11 Distribution of geological hazards in the Wuzhizhou Island
3.3.1 蜈支洲岛崩塌及断层活动地质灾害
岛岸的崩塌主要发生在基岩裸露,且节理密集发育的区域,受热带气候强烈阳光照射,在昼夜温差和热带季风的双重作用下,节理密集发育区的基岩逐渐发生崩解(图12)。
图12 蜈支洲岛东侧崩塌隐患点Fig.12 A hidden risk point of bluff collapse on the eastern side of the Wuzhizhou Island
断裂活动造成岩石错动或塌陷,在海浪的反复侵蚀和冲刷下,断裂活动的危害性进一步加剧,不仅损坏道路,也会对海岛上的建筑产生破坏。蜈支洲岛最南端转角处,长约8 m 的道路已经出现裂痕,并产生约20 cm 的错断,影响了观光车及游客安全(图13)。
图13 蜈支洲岛南部崩塌隐患点Fig.13 A hidden risk point of bluff collapse in the southern Wuzhizhou Island
3.3.2 蜈支洲岛海岸侵蚀地质灾害
蜈支洲岛岸滩整体保护较好,但在登岛实际调查过程中也发现几处岩滩存在侵蚀及部分沙滩存在淤积的状况。受海岛管理和通行条件影响,系统的实地调查难以实现,但基于无人机三维实景模型则可在室内实现系统分析研究。详细的研究表明蜈支洲岛存在7 处较为明显的岩滩侵蚀点和4 处沙滩淤积点(图11)。
岩滩侵蚀部位主要发生在节理或断裂密集发育处,长期受到海浪和潮流冲刷影响,逐步向岛陆后退,与此同时,海岛基岩在断裂、海风、海浪、太阳辐射的共同作用下逐渐崩解,被海流带入海洋(图14)。还有部分侵蚀发生在人工建筑物周边。这些侵蚀作用对海岛开发利用带来很大安全隐患。长此以往,在无其他沉积物源补充的情况下,海岛岩滩部分逐步缩小。
分布在海岛北侧的沙滩常在地形变化较大的转折处发生淤积,淤积如果发生在河道,通常造成河道和水路改道或淤堵,是一种不良地质现象[31],但这种淤积对于海岛的沙滩养护来说有其积极的一面。通常海岸的侵蚀淤积是平衡的,一个地段的侵蚀可能造成另一地段的淤积。在⑧号地质点可以清晰观察到海底岩石(附着有海草)的定向排列显示了海流的运动方向。基岩尖嘴凸出部分改变了海流的运动方向,从而导致海洋沉积物在湾口淤积,与此同时,海洋环流的作用会显著增强岩石尖嘴部分的侵蚀作用(图15)。
图15 蜈支洲岛北部砂质岸线淤积段Fig.15 A siltation section of sandy coastline in the northern Wuzhizhou Island
蜈支洲岛虽然是一个沿岸海岛,距离大陆岸线最近的距离也有2.5 km,但从本质上而言,蜈支洲岛原是大陆的一部分,地貌构架基本上是大陆向海的延伸,因海平面上升或地面沉降与大陆分离。按成因分类,蜈支洲岛属于大陆岛,其岩性组成与相邻陆地完全相同[32-33]。因此,蜈支洲岛的地质构造特征与海南岛的地质演化过程息息相关。
我国海岛大多受NE、NNE 向断裂带、坳陷和隆起的控制,其走向分布和运动方式与中生代燕山运动时形成的构造保持一致[34]。空间分布上,蜈支洲岛与其SW 方向的野薯岛,NE 方向的分界洲岛、南洲仔岛、大洲岛等海南岛近岸海岛在空间上呈NE-SW 向展布,受NE 向断裂控制作用明显。根据1∶25 万乐东县幅、陵水县幅区域地质调查报告,蜈支洲岛上的花岗岩大致形成于燕山期。有学者认为,燕山早期,由于区域岩石圈裂解拉张,在三亚地体形成了与陆内拉张有关的侵入岩,但在燕山早期,蜈支洲岛所在的三亚地体与琼中、琼北微地块地质构造环境不同。至燕山晚期,三亚地体才与琼中块体沿九所-陵水断裂发生拼贴对接[35]。
蜈支洲岛上NE、NW 和近SN 向三组断裂具备如下特征:①断裂断面平直光滑;②断裂产状稳定,沿走向延伸较远;③空间上呈等距离排列,且没有明显位移;④充填有基性岩脉的断裂面出现转弯、分叉,表现为张裂特征(图16);⑤NW 向和近SN 向断裂裂开的程度比NE 向断裂要大。蜈支洲岛上发育的3 组线性断裂表现出共轭“X”型剪节理系特征,很可能是受燕山晚期三亚地体与琼中块体发生拼贴时产生的SN 向压应力作用而形成。
图16 蜈支洲岛东南侧近SN 向张性断裂素描图Fig.16 Sketch of a tensional fracture in near S-N direction in the southeastern side of Wuzhizhou Island
调查中还发现,基性辉绿岩脉大多发育在NW向和近SN 向断裂中,局部地段密集出现(图5)。海南岛广泛发育基性岩墙群,如琼海市文市乡、三亚市南山村以及西瑁洲岛等地。研究表明三亚市南山村的辉绿岩基性岩墙群形成于90 Ma 左右,指示了海南岛南部在该时期经历了强烈的区域性拉张作用[36]。表明蜈支洲岛上的基性岩墙群很可能也是在此地质背景下形成。
一般近岸海岛面积较小,岸线蜿蜒曲折,传统遥感手段受干扰条件较多,主要用于活动区海岛整体影像,不易精细识别不同海岛岸线所涵盖的复杂地物类型[37]。无人机低空遥感技术则可以快速获取小尺度范围内高分辨率遥感影像,是海岛岸线精细化探测的有效技术手段。应用无人机高光谱遥感和机载激光雷达数据(LiDAR)开展自动或半自动海岛岸线已经取得了较好应用[37-38]。本次工作尝试应用无人机正射影像结合实景三维模型对蜈支洲岛岸线进行识别。
蜈支洲岛岸线类型主要包括基岩岸线、砂质岸线和人工岸线。其中人工岸线主要为港口码头和防波提,岸线信息提取相对容易。对于砂质岸线和基岩岸线受潮汐影响,若采用海水与岛滩(沙滩或岩滩)交接线作为岸线,则在不同时间获取的影像岸线位置完全不同,对于地形较为平缓的岛滩,差异性非常大。需要获知遥感影像拍摄的具体时间(精确到小时),再引入潮汐数据并结合DEM 数据进行校正,流程较为繁琐。因此,直接采用海水交接线来确定岸线位置比较困难,且效率较低。海岛经常受海水冲刷影响,长期以来在海岛周围形成了一条稳定的海岛痕迹线。以海岛痕迹线作为岸线界定标准时,在同等条件下,受潮汐影响较小,能够提高岸线识别和提取效率[39]。对于砂质岸线,海岛痕迹线在干滩与潮滩交界处;对于基岩岸线,海岛痕迹线在基岩的明暗交界处。这两类海岛痕迹线在DOM 数据和三维实景影像上可以清晰识别出来(图14),但由于分辨率的问题,在卫星遥感影像上难以识别潮滩与干滩的交界线(图17)。与卫星遥感影像相比,海岛痕迹线的识别具有显著优势。
图17 砂质岸段海岛痕迹线的识别Fig.17 Identification of island trace line in sandy coastline
造成海岸线变化的主要因素包括气象、海浪作用、海平面上升等自然变化以及人工堤坝、围垦、采砂等人为活动[40]。利用蜈支洲岛过往遥感影像数据,经过处理、对比分析可以比较清晰获取海岛岸线变化情况。总体而言,蜈支洲岛基岩岸线变化不大;北部砂质岸线有淤积的倾向,但受限于物源不足,形态略有变化,总体稳定;变化较大的是蜈支洲岛西北角的砂质岸线(图18)。
图18 蜈支洲岛2005 年以来岸线变化Fig.18 The coastline changes of the Wuzhizhou Island since 2005
结合实地调查发现,蜈支洲岛西北角变化较大的原因是人为开发行为与海浪侵蚀共同作用的结果。西北角是人工浴场所在地,蜈支洲岛人工浴场经历了多次搬迁,人为将海砂进行了搬运。
蜈支洲岛基岩岸线从2005 年到现在总体形态变化不大,显示了基岩岸线的稳定性。但在海浪、断裂和气象的长期共同作用下,基岩岸线形态表现出规律性特征。东南侧基岩岸线向海突出或凹陷的部分,总体方位表现为NW-SE 向;而西南侧基岩岸线向海突出或凹陷部分,总体方位为近SN 向(图18),这与蜈支洲岛的线性断裂系统表现出来的特征基本一致,表明断裂构造是控制蜈支洲岛基岩岸线形态的主控因素。亚热带气候的昼夜温差对基岩的风化作用(球状风化)以及海浪和潮汐的作用加剧了这一过程(图6)。虽然海蚀作用时刻在发生,但要观察到基岩岸线的显著变化,在人为不干预及海平面较稳定的情况下,需要更长的时间尺度。
(1)蜈支洲岛主要发育NW、NE 和近SN 向3 组线性断裂,断裂面较陡直近于直立。断裂发育的早期表现为共轭“X”型剪切性质。燕山晚期(约90 Ma),海南岛南部在经历了强烈的区域性拉张作用,导致基性岩浆上侵,并形成辉绿岩墙群。
(2)蜈支洲岛的主要土地利用类型包括:乔木林地、灌木林地、裸岩石砾地(基岩海岸、自然景观)、建设用地(旅游设施、酒店、人文景观等)。以乔木林地、灌木林地和草地总面积,计算蜈支洲岛的植被覆盖率约为77.8%。蜈支洲岛的整体土地利用率较好。但一些潜在旅游地质景观,如海蚀柱、海蚀崖、基性岩墙群等还有待于进一步开发。
(3)蜈支洲岛岛陆地质灾害主要包括岛岸崩塌、断裂活动、岸线侵蚀3 种类型。本次工作应用无人机倾斜摄影成果有效识别出岛岸崩塌隐患点5 处,这些崩塌隐患点均被分布在旅游路线上,风险性极大,建议下步及时安装防护网或人工清除。识别出明显的断裂活动点3 处(包括2 处辉绿岩脉点),断裂活动对旅游道路已造成明显的破坏,下步随着海浪的侵蚀还会进一步加剧,应尽早做好防范。识别出岛滩侵蚀点7 处,其中蜈支洲岛西北角的砂质岸线侵蚀点(图11-12)对海岛的生态地质和旅游地质带来较大影响,应尽早人工干预进行养护;对于基岩岸线的侵蚀,也需及时采取必要的养护手段,防止危及海岛旅游路线。
(4)无人机倾斜摄影测量在海岛综合地质调查中具有非常明显的优势,包括分辨率高、多维度信息提取、重点突出、低成本、安全可靠等。在海岛岸线(海岛痕迹线)识别、海岛地质灾害信息提取方面有独特优势。