测绘地理信息   2017, Vol. 42 Issue (6): 1-10
0
海洋测绘技术发展现状[PDF全文]
赵建虎1,2, 陆振波1,2, 王爱学1,2    
1. 武汉大学测绘学院,湖北 武汉,430079;
2. 武汉大学海洋研究院,湖北 武汉,430079
摘要: 简要介绍了海洋测绘的定义及内涵,从海洋大地测量、海洋导航定位、水深及水下和海岸带地形测量、海洋遥感、海洋水文测量、海底底质探测、海洋工程测量、电子海图和海洋地理信息系统9个方面介绍了海洋测绘技术的发展现状,认为海洋测绘技术正呈现出立体、高精度、高分辨率、高效的信息获取、处理和应用的态势。
关键词: 海洋测绘     测量技术     电子海图     海洋地理信息     发展现状    
Development Status of Marine Surveying and Mapping Technology
ZHAO Jianhu1,2, LU Zhenbo1,2, WANG Aixue1,2    
1. School of Geodesy and Geomatics, Wuhan University, Wuhan 430079, China;
2. Institute of Marine Science and Technology, Wuhan University, Wuhan 430079, China
Abstract: This paper briefly introduces the definition and connotation of marine surveying and mapping and overviews its development status in nine aspects including marine geodesy, marine navigation and positioning, bathymetry and coastal topographic survey, ocean remote sensing, marine hydrographic surveying, seafloor sediment detection, marine engineering survey, electronic chart and marine geographic information system. A conclusion is reached that marine surveying and mapping technology is showing a situation of three-dimensional, high-precision, high-resolution and high-efficiency in information data acquisition, processing and application.
Key words: marine surveying and mapping     measurement technology     electronic chart     marine geographic information     development status    

海洋测绘是获取和描绘海洋、江河、湖泊等水体和包围水体各对象的基础地理要素及其几何和物理属性信息的理论和技术,是测绘科学与技术的一个重要分支,是一切海洋军事、海洋科学研究及开发和利用活动的基础[1]。海洋测绘主要包括海洋测量、各类测量信息的绘制以及信息的综合管理和利用三大内容。海洋测量是对水体、水底、周围陆地进行测量的理论、技术和方法,按照研究内容和任务可分为海洋大地测量、海洋导航定位、水深测量及水下、海岸带地形测量、海洋遥感、海洋水文测量、海底声学底质测量以及海洋工程测量等内容[2]。海图绘制是综合呈现和表达测量信息的工作,是设计与制作海图的理论、方法和技术的总称。海洋地理信息系统是对海洋的空间信息处理、管理、显示、分析和应用的技术和方法。根据海洋测绘的定义和内涵,结合当前进展,下面介绍各海洋测绘技术的发展现状。

1 海洋大地测量

海洋大地测量以确保海洋测量控制基准为目的,是为海洋测绘建立物理(重力和磁力)和几何(平面和高程)基准体系与维护框架的大地测量技术,是陆地大地测量在海区的扩展。

1.1 海洋物理(重力和磁力)大地测量

海洋重力、磁力测量属于海洋物理大地测量,是建立海洋重力、磁力测量基准,获取海洋重力磁力基础数据、开展海洋重力场模型和地磁场模型构建、应用重力磁力异常变化开展海洋资源和目标探测的基础性海洋测量工作。

海洋重力测量目前呈现立体测量态势,按照测量载体可分为星基、机载、船基和沉箱式海底重力测量。船载重力测量仍是目前主要采用的高精度重力测量方式。近年来,机载重力测量发展迅速,已初步具备实际生产能力。我国的海洋重力仪长期依赖进口,近年在引进、消化和吸收的基础上,国产重力仪研发进展较快,目前已完成实验验证,实现了数据的自动采集和规范处理,性能指标接近国外同类产品[3]。海洋重力测量数据处理进展目前主要体现在四个方面:一是测量数据采集与处理实现全过程自动化与智能化;二是重力仪性能评价实现了技术流程标准化和评价指标的系统化与定量化;三是精细化数据处理方法体系更趋严密,成果精度显著提高;四是构建了多源海洋重力数据融合处理理论[4]。具体进展表现为:研发了基于电子海图的导航、数据采集、处理和成图软件;提出了重力仪零点趋势性漂移、有色观测噪声与随机误差的分离方法,形成了稳定性评估的标准化技术流程;推导了均方根误差、系统偏差、标准差等新的重复测线内符合精度评估公式,为重力仪动态性能评估提出了更精细的评估指标;提出了基于互相关的交叉耦合效应修正方法,对高动态测量数据实施了综合补偿和精处理;基于Tikhonov正则化方法和移去-恢复技术,构建了多源重力数据融合模型,提出了融合多源重力数据的纯解析方法;开展了多源重力数据同化、融合理论和方法研究,给出了中国海域1′×1′垂线偏差模型;联合EGM2008模型确定了1′×1′重力异常模型,在开阔、近海海域分别取得了±3 mGal和±(7~9) mGal的精度。

相对于立体海洋重力测量,海洋磁力测量目前仍以船基拖曳测量方式为主。近年来的技术突破主要表现在如下几个方面:开展了海底地磁日变站布放选址方法研究,解决了远海磁力测量日变改正难题;在拖鱼入水深度计算与控制方面,建立了入水深度与配重、拖缆长度和船速间的关系模型;采用傅立叶谱分析方法,实现了磁平静日变和磁扰改正的分离,解决了强磁扰期日变改正问题;基于多台海洋磁力仪、测深仪和GNSS(global navigation satellite system),构建了一种阵列式海洋磁力测量系统[5];在海岛礁地磁力测量方面,实现了地磁经纬仪、陀螺经纬仪、天文观测和GNSS高精度定位与定向等多系统一体化集成应用,给出了完整的地磁三分量测量技术流程,编制了技术规程,建立了相应的数据处理模型,提高了海洋地磁测量的精度;借助实测磁力数据,开展了几何建模法和基于地磁位的建模法研究,构建了海洋局域地磁场模型,取得了优于10 nT的建模精度,满足了地磁匹配导航的需要[6]

1.2 海洋控制网建立

海洋控制网建立是布设、施测覆盖海岸、岛礁、水体和海底控制网的工作,是大地控制网的组成部分,是陆地平面坐标框架网向海洋的延伸。海岸、岛-陆、岛礁控制网建设和施测方法相对成熟,目前主要借助GNSS技术来实现。海底控制网建设尚处于起步阶段,受到了各国的高度重视。海底控制网建设研究目前主要聚焦于网址选择、网型优化、绝对平面基准水下传递、控制网联合测量及处理、网点维护等方面。具体的技术进展主要包括声速反演和局域声速场的建立、高精度高效声线跟踪、顾及深度约束的平面绝对基准传递、联合测量方法及顾及水深的联合网平差方法等[7]。海底控制网建设主要借助声学测量技术来实现,涉及的声学设备有换能器、应答器及水听器,目前测量中仍以国外设备为主,国内科研院所和公司,如哈尔滨工业大学和中国科学院声学研究所,研发进展较快,设备已投入了实际应用。

1.3 海洋垂直基准

海洋垂直基准是海洋垂直测量成果的起算参考,包括陆地高程基准、平均海平面和深度基准面。海洋垂直基准通常借助潮位站潮位数据来确定,随着卫星测高、GNSS等相关技术的发展,海洋垂直基准确定采用的数据源和表达方式发生了深刻的变革[8]

GNSS高精度垂直解已广泛应用于海洋测绘中,在GNSS一体化测深、GNSS潮位测量等技术中发挥着重要的作用,但均需解决如何将GNSS实测大地高转换为基于深度基准面的海图高问题。目前,海洋垂直基准转换以平均海平面或似大地水准面为中介,实现大地高到正常高、正常高到海图高的两步转换,也可实现大地高到海图高的一步转换。无论采用何种转换,均需解决三个问题,即垂直参考基准面选取,深度基准面的无缝化和垂直基准间的无缝转换。椭球面或平均海平面常被选做垂直参考基准面。深度基准面定义于验潮站,是离散、跳变和不连续的,其无缝化问题是三个问题中的关键。河口水域潮位站密集,各站上大地高、正常高和海图高关系明确,可实现彼此转换。潮位站包围水域的深度基准面无缝化问题可借助几何法或潮差比法内插获得[9]。据此机理已在多个水域实现了深度基准面无缝化及垂直基准的无缝转换。全球潮汐模型和全球平均海平面模型为任意海域垂直基准面的无缝构建及垂直基准间的无缝转换提供了条件。借助全球潮汐模型可获得深度基准面相对平均海平面的差距,基于全球平均海平面模型可得到基于椭球面的平均海平面值,联合二者可实现椭球面、平均海平面和深度基准面之间的无缝转换[10]。但在近海岸,全球平均海平面模型和全球潮汐模型的精化问题应引起足够重视。

2 海洋导航定位技术

一切海洋活动均无法离开位置服务。海上位置服务目前主要借助GNSS来提供。一般船只导航多采用GNSS单点定位技术;中小比例尺水下地形测量中,导航定位多采用GNSS广域差分或星际差分技术;在高精度测量中,定位主要采用GNSS RTK(real-time kinematics)、PPK(post-processing kinematics)和PPP (precise point positioning)定位技术[11]

水下导航定位多采用水下声学定位系统,如LBL(long baseline)、SBL(short baseline)和USBL(ultra-short baseline)。LBL、SBL和USBL均采用交会定位方法,且常组合使用[12]。此外,为了提高水下导航定位的精度和稳健性,声学定位技术常与惯性导航系统、航位推算系统等组合使用。近年来,为确保水下潜器的隐蔽性、导航的连续性和长时性,惯性导航系统常与海底几何场(地形、地貌)或物理场(重力、磁力)的匹配导航技术组合,形成(无源)自主导航定位系统,为水下潜器导航服务。

3 水深测量及水下和海岸带地形获取

海底、海岸带地形测量是海洋测量的基本内容之一。传统海底地形测量需要通过深度、潮位、吃水、涌浪等测量才能获得,现代水下、海岸带地形测量因相关技术的进步已发生深刻的变革,呈现出立体、高效、高精度测量态势。

3.1 水深测量

水深测量是水下地形测量的重要组成部分,目前主要采用单波束测深、多波束测深系统(multibeam bathymetric system, MBS)和机载激光(LiDAR)全覆盖测深等技术。单波束一次测量只能获取一个测深点,适用于中小比例尺或小区域大比例尺水下地形测量;MBS一次测量可在航行正交扇面内获得上百到几百个测点,实现对海底全覆盖扫测。无论是单波束还是MBS,我国已具有自主研发能力,目前已形成国内外设备联合使用的格局。

3.2 船基海底地形测量

目前的水下地形测量呈现传统与现代测量模式并存的态势。GNSS一体化水深测量技术是现代船基水深测量的代表,集单波束、多波束测深技术、GNSS RTK、PPK、PPP高精度定位技术、POS(position and orientation system)技术和声速在航测量技术等于一体,在航实现多源信息综合采集和融合,最大限度地削弱各项误差影响,提高海底地形测量的精度和效率。GNSS三维解在该技术中起着为测深提供平面和垂直基准的重要作用。联合GNSS三维解、船姿、船体方位及换能器和GNSS天线在船体坐标系下的坐标,实时获得换能器的三维坐标。GNSS一体化测深技术的优势在于无需开展潮位站潮位观测,较彻底地补偿了波浪和声速等因素对测量的影响[13];在航实现海底测点三维坐标确定,因此显著提高了测量精度和效率。船基海底地形测量的无人化和自动化充分体现在无人船海底地形测量技术中。除此之外,无人船还配备自动操控、避碰、无线电等系统[14]。无人船地形测量降低了作业成本,提高了效率,在海况良好的大水域测量、浅滩等危险或困难水域测量中作用明显[15]

测深数据处理一直是海底地形测量中研究的热点问题之一,进展集中体现在:①声速及声线跟踪方法、声速剖面简化。声线跟踪主要有层内常声速、常梯度声线跟踪方法、误差修正法和等效声速剖面法[16]。SVP(sound velocity profile)简化是目前该领域研究的另一个热点之一[17]。SVP简化目前主要有滑动平均法、D-P(Douglas-Peucker)法、MOV(maximum offset of sound velocity)法和面积差法。SVP简化极大地缩短了声线跟踪时间,提高了深海批量多波束测深数据处理效率,但在简化率和可靠性的均衡研究仍需深入。②测深数据滤波。受海洋环境、系统噪声和测量船噪声等影响,测深数据滤波备受关注。常用的方法主要有COP法、Ware、Knight & Wells法、Eag(RDANH)法、趋势面法、抗差估计法、Bayes估计法、中/均值滤波、选权迭代加权平均、CUBE(combined uncertainty bathymetry estimation)滤波等。CUBE滤波因其高效、可靠、抗差、稳健等特点被广泛采用[18],但在碎石区滤波中滤波效果欠佳,需给予改进。③残余误差综合影响削弱。除受自身测量误差影响外,多波束测深还受声速、姿态、安装偏差等的影响。虽对其严格测定和补偿,但残余误差仍给测深带来系统性综合影响,导致“哭脸”或“笑脸”地形。基于地形频谱特征的削弱法认为,该影响仅改变了测深结果对地形变化趋势的反映,未影响对微地形的呈现。据此,利用高精度中央波束测深结果构建边缘波束地形趋势,联合实测微地形,合成边缘波束地形,从机理上削弱了残余误差的综合影响[19]。国内测深数据处理软件研制已打破了国外垄断,并开展了应用。

3.3 潜基/深拖海底地形测量技术

为提高海底地形信息获取的分辨率和精度,满足海洋科学研究和工程应用需要,以AUV(antonomous underwater vehicle)、ROV(remotely operated vehicle)、深拖系统为平台,携载多波束测深系统、压力传感器、USBL等设备的潜基海底地形地貌测量系统已经面世,并在我国一些重点勘测水域和工程中得到了应用,并受到海事、水下考古、海洋调查等部门的高度重视[20]。潜基海底地形地貌测量系统借助USBL、罗经、姿态传感器和压力传感器为平台提供绝对平面和垂直坐标,利用多波束获取海底地形,并将信息通过电缆传输到船载存储和处理单元,综合计算获得海底地形。潜基测量技术适用于深海地形测量。随着我国深海调查活动的深入,其应用必将越来越广泛。

3.4 反演技术 3.4.1 卫星遥感反演水深

卫星遥感反演水深是借助可见光在水中传播和反射后的光谱变化,结合实测水深,构建反演模型,实现大面积水深反演,再结合遥感成像时刻水位反算得到海底地形[21]。目前可用卫星主要有IRS、Ikonos、QuickBird、AVIRIS、Sentinel-2、Landsat、TM、SPOT等。卫星遥感反演水深具有经济、灵活等优点,但反演精度及范围需提高。

反演水深的关键是构建不同波段或组合波段与水深间的反演模型,主要包括波段优选、波段组合及反演模型构建三部分[22]。波段优选是借助主成分分析法或相关法提取与水深强相关波段的工作。波段组合是分析不同显著波段组合对反演水深精度的改善程度,进而确定最优组合波段的工作。反演模型构建是建立显著波段或组合波段与实测水深关系模型的工作,如线性模型、附加幂函数非线性修正的线性模型、基于底反射建立的单、双、三波段模型等[23]。对于不同卫星和不同水域,最优波段选择、最佳波段组合及反演模型均存在差异。

3.4.2 重力反演海底地形

重力异常和海底地形在一定波段内存在高度相关,据此可反演大尺度海底地形。重力反演地形经历了从一维到二维线性滤波的发展,其核心是反演模型构建。

反演模型构建经历了直接建模和修正建模过程,目前多采用修正建模。如利用ETOPO5模型、GMT海岸线数据、卫星测高重力异常和船测水深,一些学者建立了海底地形模型[24]。采用垂直重力梯度异常可以反演得到独立于重力异常的海底地形模型。在不同海底模型假设基础上,如在椭圆形海山模型假设基础上,利用垂直重力梯度异常、采用非线性反演法对全球的海山分布进行了反演,结合高斯海山模型,通过分析地壳密度、岩石圈有效弹性厚度及截断波长对反演的影响,采用垂直重力梯度异常反演得到海底地形。

借助重力地质法(gravity geologic method,GGM),利用大地水准面数据,在频域内采用二维反演技术,以迭代法处理海底地形和大地水准面的高次项问题,削弱岩石圈挠曲强度的误差影响,改善反演精度;采用快速模拟退火法,利用重力垂直梯度可反演海底地形[25]。比较GGM、导纳法、SAS (Smith and Sandwell)法、垂直重力梯度异常法和最小二乘配置法,GGM和SAS法反演精度较高,GGM、垂直重力梯度异常法和最小二乘配置法适宜开展大面积海底地形反演。

3.4.3 声呐图像反演高分辨率海底地形

高分辨率海底地形在沉船打捞、油气勘探和环境监测等海洋工程和科学研究中发挥着重要作用。在浅水水域,高精度和高分辨率海底地形主要借助多波束获得,但在深水水域,其测深分辨率会随波束入射角和水深增大而显著降低。侧扫声呐通过深拖可获得20~100倍于测深分辨率的海底声呐图像,但缺少高程信息。基于侧扫声呐成像机理及光照理论,借助SFS(shape from shading)方法可实现基于声呐图像的海底高分辨率地形反演。SFS方法是基于声波在海床表面遵循的海底反射理论,通过构建回波强度与入射方向、地形梯度等因素之间的关系,对模型求解即可得到海床地形[26]。SFS反演仅能得到相对形状,需借助外部测深数据或侧扫声呐测量中提取出的水深数据约束,才能实现绝对海底地形的恢复[27, 28]

3.5 海岸带一体化地形测量技术

海岸带地形是海底地形与岸上地形的连接部分,也是海洋测量中的一个难点和热点。传统海岸带地形测量多采用全站仪或RTK人工完成,但效率低、测量困难。近年来,机载LiDAR、航空摄影、水下近岸一体化测量等技术得到广泛应用。机载LiDAR可实现水下和干出地形一体化测量,但激光穿透能力和测量精度受海水浑浊度的影响较大。航空摄影测量在无植被的干出部分效果显著,但在芦苇或杂草覆盖的海岸带则受到限制。利用设备的非接触、面扫测特点,出现了集多波束、激光扫描仪、稳定平台、POS等于一体,安装在测量船或气垫船上的一体化测量系统,同步实现水深及岸边地形的测量[29]。一体化测量系统在堤坝、码头等水域有较好的应用,但在一般浅滩地带则存在测量盲区。

4 海洋遥感 4.1 卫星遥感

借助卫星遥感可以对海洋进行实时、全方位的立体监测,长期获得稳定可靠的海洋观测资料[30]。目前可用的卫星遥感数据源有MODIS、MERIS、GOCI、Sentinel-2、Landsat 8、Landsat TM、SPOT、Hyperion、QuickBird、WorldView-2、Pleiades-1等。借助这些数据,可开展海岸带、植被、海洋生态监测及水深反演等方面的工作。

岸线监测算法研究方面,采用边缘提取算子如Canny算子、Sobel算子或利用归一化水体指数、修复归一化水体指数法进行水陆分离,利用轮廓边界跟踪技术提取岸线;利用归一化水体指数确定岸线位置;根据海岸线的多源多时相影像特征和空间特征,与多种非遥感信息资料相结合,采用对照分析法,实现海岸线长度、位置、类型提取[31];根据岸线的时空和分形维数变化,可研究岸线的演变[32]。植被监测算法研究方面,利用中等空间分辨率的单时相、多时相或融合的多光谱影像可研究海岸带植被覆盖度、分类、生物量、碳储量、净初级生产力等。水色监测算法研究方面,通过获取海面上行的离水辐亮度,经大气校正和水色反演,可得到水体中浮游植物的色素浓度、悬浮颗粒浓度等信息[33];水色反演已从单一的叶绿素浓度扩展到透明度、悬浮泥沙浓度、黄色物质浓度等参数;水色反演算法的区域性特征明显,基于大量现场数据同步采集的基础上,不断优化现有算法是未来研究的方向。

4.2 机载遥感测量技术

机载遥感测量主要借助机载可见光相机、可见光摄像机、红外相机、高光谱成像仪、LiDAR、SAR(synthetic aperture radar)、合成孔径雷达等开展海岸带地形测量,岸线、植被、水色等监测,采用技术与卫星遥感近似。机载激光雷达测深系统ALB(airborne LiDAR bathymetry)是一种即可用于地形测量又可用于环境监测的技术,近年来,无论在系统还是应用研究等方面,均取得了长足发展。

ALB系统主要由激光扫描仪、POS、高度计等组成。按照设计测线,飞行过程中快速扫描实现海底地形全覆盖测量。以飞机平台为基准,ALB分别借助红、绿激光获得海面和海底高,测量中同步采集飞机位置、姿态和航向,飞机位置借助GNSS RTK、PPK、PPP获得,姿态和航向借助POS获得。此外,联合外部获得的海水浑浊度及部分船基实测海底地形,构建修正模型,提高ALB测量精度。

ALB系统研发已从实用化迈入商业化,脉冲发射频率、探测能力得到了提高,系统体积、重量和能耗显著减小,机动性和续航时间增强[34]。典型的ALB系统主要有Optech的SHOALS 200/400/1000/3000、CZMIL和Aquarius系列产品,AHAB的HawkEye Ⅱ/Ⅲ和Chiroptera,Fugro的LADS MK3,RIEGL的VQ-820-G和VQ-880-G。我国目前也在研制ALB系统,多处于研发阶段。

目前,ALB数据处理软件多为各公司研制的随机软件,国内尚无相应的软件系统。ALB数据处理研究进展主要表现在如下几个方面:①ALB测深理论[2]。激光测深能力与水体散射系数和衰减系数比值强相关,借助唯像理论可建立激光测量的唯像雷达方程;激光束虽具有一定发射角,但其传输规律仍可用准直光束传输特性来描述,据此可建立准直光束在海水中传输的唯像理论模型;影响水底回波振幅的因素主要有水底反射率和脉冲展宽,基于激光辐射传输模型可对水底回波振幅进行校正。②归位计算。根据GNSS提供的激光扫描仪三维绝对坐标,结合飞机姿态、激光扫描模式及扫描角、往返测量时间,可归算海面点的三维坐标;根据红、绿激光测量时差、海水折射率、波束扫描角,归算绿激光海底圆斑的三维坐标。③波形识别。波形识别是检测激光回波、获取水面和海底波束传播时间,进而计算深度的关键。目前采用的技术主要有:为抑制白天强背景噪声,更精确地提取激光回波信号,对回波信号首先开展高通滤波滤除低频信号,再识别两种高频脉冲;利用回波信号的上升时间及振幅等特征,采用半波峰法识别海表和海底回波信号,进而估算水深;采用窄脉冲、高速探测器、小接收视场、窄带干涉滤光片和正交偏振方式接收,改善浅水海表和海底反射信号叠加;采用双高斯脉冲拟合,从极浅海水回波中分离海表和海底脉冲,实现水深提取。④浑浊度反演。海水浑浊度会引起激光能量衰减,影响激光回波波形;反之,根据激光水体回波特征可估计海水浑浊度。提取ALB原始波形数据后,分析激光水体后向散射波形,估计有效衰减系数,进而反演海水浑浊度。⑤绿激光高度修正及单一绿激光测量。受海表渗透深度影响,绿激光海表测量存在不可靠性。利用红、绿激光测量结果分析绿激光水表渗透深度空间变化,利用统计法对绿激光海表高程进行修正,提高绿激光海表测量的精度[35]。采用逐步回归法建立关于泥沙含量、波束扫描角和传感器高度的绿激光水表渗透深度模型,推导绿激光高度修正模型,据此对绿激光海表和海底高程进行修正,实现基于单一绿激光的高精度海底地形测量[36]。⑥深度偏差修正。几何发散和多次散射使绿激光底回波产生脉冲展宽效应,引起波峰位移,导致测深产生偏差。深度偏差主要与ALB系统测量参数(波束扫描角、传感器高度)和海水水文参数(水深、浑浊度)有关。采用逐步回归法建立关于水深、波束扫描角、传感器高度和海水浑浊度的深度偏差模型,据此对ALB进行深度修正,实现ALB高精度水深测量。

4.3 声呐遥感 4.3.1 侧扫声呐测量

侧扫声呐系统SSS(side-scan sonar)是常用的条带式海底成像设备,借助拖鱼上左、右舷换能器阵列发射的宽扫幅波束,并在走航过程中对海底进行线扫描,进而形成可反映水体、海底目标分布和地貌特征的条带图像。目前SSS系统向多频段、多脉冲、多波束、深拖及同时具备测深及成像功能方向发展,如EdgeTech系列、Klein系列、Triton C3D系列以及适用于深拖的Kongsberg Maritime DeepTow系列。我国的SSS研制起步较晚,但发展迅速,目前已实现商业化,并在生产中取得了较好的应用效果。SSS图像主要用于寻找目标,近年来在图像精处理及目标探测和识别方面取得了长足发展,具体表现在条带图像均衡化、海底线综合追踪、多条带声呐图像均衡化、基于地理编码和共视特征匹配的多条带图像拼接、大区域海底图像生成及基于精细化处理后声呐图像的目标自动检测和识别等方面[37]

4.3.2 多波束成像技术

MBS具有测深功能,还具有接收海底和水体回波,形成海底回波图像和水体图像的能力。MBS底回波图像有三类,即由平均波束强度、波束序列片段(Snippet)强度和伪侧扫回波强度形成的声呐图像。Snippet数据的垂航分辨率远高于平均波束强度,尤其在大入射角区。考虑高分辨率图像特点,Snippet图像方面的研究相对较多,主要聚焦于Snippet图像的AR改正、图像均衡化、Snippet图像与SSS图像的融合等方面[38]。MBS水体回波可形成水柱图像,进而实现水体中目标的探测。围绕水柱图像的研究主要聚焦于水柱图像消噪、基于噪声特点的水体中目标自动探测和识别几个方面。

4.3.3 合成孔径声呐

合成孔径声呐(synthetic aperture sonar,SAS)是一种利用合成孔径技术的侧扫式主动成像声呐。与SSS相比,SAS图像具有更高的径向分辨率,且与距离无关。装备有高、低频换能器的合成孔径声呐可同时获得高、低频声呐图像,可以清晰地呈现海底地貌,以及海床下一定深度的目标,可全面地反映管线的分布。SAS的作业模式与SSS相同。目前SAS正在向小型化发展,测量技术和性能不断完善。我国的SAS技术基本与国外处于同级,自主产品已经投入实际应用。SAS数据处理研究进展主要表现在条带图像的快速生成、处理、拼接、大区域地貌图像生成、高低频图像融合及基于SAS图像的目标探测和识别方面。

4.3.4 二维声呐成像

二维声呐成像近年来在工程中应用越来越广泛。通常采用基于船基悬臂或座底工作模式,近场(小于100 m)对目标开展快速线扫描成像。二维声呐成像系统目前主要采用国外设备,如Kongsberg MS1000(机械)、BlueViewM900(实时)等。与MBS水柱数据处理类似,相关研究主要聚焦于高质量影像获取和基于图像的目标检测。因二维声呐成像只能获取平面图像,无法反映目标三维形态,研究目标三维形态参数挖掘算法是未来的一个发展方向。

4.4 光学近景摄影技术

二维声呐成像虽满足了对近场目标和环境信息的获取需要,但存在分辨率偏低、对细微特征难以捕捉等不足。二维声呐成像系统常与水下光学成像系统(相机)配套使用,即利用二维声呐成像系统快速发现目标,再利用光学成像系统接近目标,获取分辨率和清晰度更高的目标图像。水下光学成像质量受水质影响较大,只能在清澈的海水环境下采用。

5 海洋水文测量

水文测量是海洋测量的组成部分,包括的要素主要有水位、流速、温度、盐度、水色、透明度、含沙量、浑浊度等。诸要素中,水位和流速与其他海洋测量相关度较高,发展较快。

水位获取研究方面,近年来开展了基于GNSS RTK/PPK/PPP的在航、锚定潮位测量和提取方法研究,实现了潮位的高精度获取[39];借助GNSS潮位,通过基准关系研究,实现了基于GNSS潮位的平均海平面和深度基准面传递;开展了基于潮位站潮位、潮汐模型的潮位预报方法研究,余水位推算方法研究,实现了水位的多样化、高精度获取。

流速测量是获取流速流向,构建流场,反映流速分布和变化规律的工作。流速测量按照测量对象的不同分为固定水层流速测量和垂线流速测量两类。前者主要借助流速流向仪测量,后者主要借助ADCP(acoustic doppler current profiler)来获取。ADCP因测量效率和精度高、可获取垂线流速等特点在生产中广泛应用。围绕流速的研究目前主要聚焦于基于外部传感器的ADCP流速测量、基于有限个座底ADCP的宽断面流量实时估计、流速数据的精处理(如流速数据质量控制、盲区流速推算、等距等时流速内插等)、基于径向基函数的潮流分离和局域流场构建等方面[40]

6 声学底质探测

底质探测是获取海床表面及浅表层沉积物类型、分布等信息的技术,是海洋动力学研究、海洋矿产资源开发和利用、舰船锚泊、水下潜器座底等海洋科学、经济、军事的基础数据,是海洋测量的内容之一。海底底质通常借助采样器取样或钻孔取芯,通过实验室分析获得,存在效率低、成本高等缺陷。声学底质测量借助声波回波特征与底质的相关性实现底质探测,具有探测底质效率和分辨率高的特点,是传统底质取样探测的一种很好的补充方法。声学底质探测研究近年来发展迅速,集中体现在底质声学测量和声学底质分类两个方面。

6.1 底质声学测量

底质声学测量是借助声学换能器测量来自海床表面或海底浅表层底质层界的回波强度的工作。近年来,研究测量主要聚焦于声呐测量和回波强度数据处理两个方面。

海床表面底质声学信息可借助单波束回声测深仪、MBS和SSS来获取,浅表层声学信息主要借助浅地层剖面仪或单道地震来获取。为从以上设备接收的回波强度信息中提取出底质特征,需对回波强度进行质量控制、各项补偿和修正的处理[41]。对于海床表面回波,数据处理进展包括:①质量控制。常采用统计方法如滑动平均来消除回波强度序列中的异常观测值。②声速改正及传播损失改正。根据声速和声波的往返程传播时间,通过声线跟踪,计算波束传播到海底的实际声程,进而用于声波在海水中传播损失的计算和补偿。③波束模式改正。根据波束在海底反射模式随入射角的变化而变化这一机理,小入射角(≤25°)时采用镜反射改正模型或线性改正模型,大入射角(>25°)时采用兰伯特改正模型。④海底地形坡度改正。根据测深数据提供的地形坡度计算和补偿地形因素对声强的影响量。对于浅表层回波强度,数据处理进展包括:①数据预处理(振幅的希尔伯特变换、滑动滤波、强度到灰度的转换等);②层位综合拾取(灰度突变法、层位追踪法和层位生长法);③数据处理软件的研制。

6.2 声学底质分类

借助海底底质的声学回波强度特征参数或统计特征参数进行底质划分。目前主要采用的方法有声学参数反演法和声波回波强度统计特征分类法。

声学参数反演法基于底质对声波信号相干分量贡献的差异,通过反演海底沉积物声阻抗、声吸收系数等声学参数,结合不同沉积物密度、声速、孔隙率和颗粒度等物理参数,构建经验模型,实现不同声学特征参数到实际底质类型的映射模型构建,进而实现底质声学分类[42]

声波回波强度统计特征分类法利用不同底质回波强度或振幅的统计特征参量,借助一定的聚类分析方法,通过构建分类器实现不同底质类型的划分[43]。统计特征参数主要有:反映声阻抗变化和界面粗糙度的平均值、标准差和高阶矩等;评估回波强度分布的分位数和直方图;描述回波强度功率谱的特征量;用于纹理分析的灰度共生矩阵;对于回波强度和深度变化的分布和结构非常敏感的分形维;回波强度的角度依赖性特征等特征参数。采用的聚类分析法主要有模板匹配分类法、判别函数分类法、神经网络分类法和聚类分析法。

按照是否具备先验底质样本,声学底质分类法又分为监督分类法和非监督分类法。监督分类法通过寻找先验底质样本与对应位置回波强度间的关系实现底质分类,采用的分类方法有模板匹配法、判别函数法和神经网络分类法;非监督分类无需先验底质样本,根据预分类底质的回波强度间的相似性关系实现聚类,常采用自组织神经网络分类法[44]、聚类分析法。

7 海洋工程测量

海洋工程测量是海洋工程建设中勘查设计、施工建造和运行管理过程的测量技术,是海洋测量的组成部分。测量内容几乎包括所有海洋测绘内容,既有单一属性要素测量的特点,又有围绕工程服务的特殊性。实际测量中可根据工程需要组合测量,也可为满足测量要求对现有的测量方法改进,并通过多源测量信息融合,最终实现测量对象信息的全方位获取。

近年来,随着海洋工程活动的增加和作业难度的增大,海洋工程测量技术受到施工单位的高度重视,技术研发进展迅速。如港珠澳大桥工程海底隧道段中的管节精确安放是整个工程的一个难点,在浅水区需要开展RTK定位、全站仪自动测量、水下声学定位、管节姿态和方位测量以及为管节拼接研制的拉线系统开展的测距和定向测量,根据不同阶段测量精度和各测量方法特点,对多源测量信息融合,综合实现管节的沉放和对接;在深水区,则组合研制的水下声学定位系统、管节姿态和方位测量系统、拉线系统,综合实现管节的安放和对接。水下坝体检测中,常组合二维声呐、水下相机、声学定位和潜器定姿等系统,通过测量和多源信息融合,实现坝体裂缝检测和定位[45]。码头建设中,常综合开展水下地形测量、水位观测、底质探测、水质调查等测量内容。桥墩冲刷检测中,常开展水下地形、声呐成像、潮位测量等内容,还需根据桥墩冲刷形态对测深设备进行改进,实现冲刷形状的复原。海上钻机平台安装和监测中,需开展海底控制网建设、水下定位、姿态测量、声学和光学成像等测量,实现安装和平台运行过程中的各项参数监控。

8 电子海图绘制

海图是所有测量要素的综合承载体,目前纸质海图虽仍在沿用,但电子海图更为普及,其进展主要体现在如下几个方面。

符号及标记方面,基于IHO相关标准和规范,采用文本描述法,设计了“所见即所得”的海图符号编辑器;提出了一种基于字符颜色扩展的海图水深注记表示方法;研究了海岛礁符号分类,实现了海岛礁符号的科学编码,确立了海岛礁符号设计方法,并开发了符号库系统。

在海图主体标准研究方面,分析了S-100的框架与机制,比较了S-100与S-57的不同,展望了基于S-100的S-101电子海图生产规范的发展情况。对水深表面产品规范标准S-102进行了研究,提出了改进建议,探讨了S-99 1.1.0版相比S-99 1.0.0的主要变化[46]

生产工艺方面,研究了海图配准、电子海图中数字接边、数字化海图制图中点状要素注记自动配置和航海图书生产流程中的色彩管理方案等问题,提出了一种英版航海通告信息自动搜集与处理技术,进一步完善了电子海图生产工艺。

应用方面,提出了电子海图云服务概念,设计了云环境下的海图集合论数据模型,提出了海图集合的云存储策略,建立了云环境下的空间索引模型,提出了全球电子海图的云可视化服务方案,研究了云计算环境下电子海图网络服务的部署方法。以自主知识产权电子海图控件为显示核心,建立电子海图功能服务,并按照网络地图服务标准发布,实现了浏览器/服务器模式的电子海图的发布。利用MapServer平台的开源、开放、跨平台特性以及支持S-57电子海图数据的功能,深入研究基于MapFile的海图数据访问、制图表达等关键问题,更好地促进Web电子海图的应用。研究了电子海图在AUV区域搜索任务中的应用,开发了电子海图遥感溢油识别显示应用平台,设计了基于电子海图基础平台的海洋调查方案辅助生成系统,深化了电子海图的专业化应用。

在电子海图导航方面,开展了中国海区e-航海原型系统技术架构研究,对中国海区e-航海建设进行了全面论述,提出了以e-航海系统为关键环节的“智慧港口”概念,提出了以服务天津港复式航道通航安全为核心的天津港e-航海试点工程建设的总体设想[47]

9 海洋地理信息

深入研究了海洋地理信息系统(marine geographic information system, MGIS)理论构成体系中的时空数据模型、时空场特征分析、信息可视化和信息服务等技术,通过Multipatch格式扩充CDC(changed data capture)数据,实现了从二维CDC格式数字海图和海洋测量数据快速构建三维空间的方法;研究了数字海洋系统中电子海图数据融合可视化问题,提出了温跃层数据的自动提取和三维表达的理论与实现方法,实现了可视化海洋环境空间数据的动态演示,形象地表达了海洋环境空间分布。

在数字MGIS建设方面,沿用S-57标准数据结构的部分特性,采用面向对象的思想,设计了满足ENC_SDE要求的系统电子海图空间数据库的空间数据模型,支持了电子海图空间数据的统一管理;优化了两级空间索引算法,设计了数据库存储文件的空间数据组织结构;提出了港口航行信息数据集成的组织方法,构建了港口信息数据模型;提出了海洋测绘产品的标准化、海洋测绘质量管理体系的标准化和海洋测绘生产体系的标准化等构想。

在应用方面,从数据特征和用户需求出发,研发了集成数据管理与查询、数据处理与分析和数据可视化功能于一体的南海海洋信息集成服务系统;提出了“虚拟港湾”的概念,并以天津海岸带“虚拟港湾”仿真平台建设为原型,详细说明了“虚拟港湾”仿真平台建设的技术原理和技术路线[48];积极推进了“数字海洋”建设,实现了数据采集、全景图像生成技术、三维全景实景建库等关键技术,研发了数据库服务、三维全景实景显示漫游和渔政地图等子系统。研制了海洋多源异构数据转换系统,设计了可实现海洋数据解译与再存储的统一数据存储结构,搭建了海洋水文环境要素可视化系统,基于面向数字海洋应用的虚拟海洋三维可视化仿真引擎——i4Ocean,模拟了海上溢油现象[49]

10 结束语

本文系统地介绍了海洋测绘的定义及内涵,总结了海洋大地测量、海洋导航定位、水深及水下、海岸带地形测量、海洋遥感、海洋水文测量、海底底质探测、海洋工程测量、电子海图和海洋地理信息系统等海洋测绘任务的定义及各个领域的发展现状,认为当前的海洋测绘技术正呈现出立体、高精度、高分辨率、高效的信息获取、处理和应用的发展态势。

参考文献
[1] 翟国君, 黄谟涛, 欧阳永忠, 等. 关于海道测量与海洋测量的定义问题[J]. 海洋测绘, 2012, 32(3): 65–72
[2] 王克平, 欧阳永忠, 蒋红燕, 等. 2014年《海洋测绘》论文综述[J]. 海洋测绘, 2015, 35(6): 78–82
[3] 宁津生, 黄谟涛, 欧阳永忠, 等. 海空重力测量技术进展[J]. 海洋测绘, 2014, 34(3): 67–72, 76
[4] 邓凯亮, 暴景阳, 黄谟涛. 航空重力数据向下延拓的Tikhonov正则化法仿真研究[J]. 武汉大学学报·信息科学版, 2010, 35(12): 1 414–1 417
[5] 宗发保, 邓瑞辉, 任来平. 一种阵列式海洋磁力测量系统[J]. 海洋测绘, 2015, 35(1): 45–47
[6] 秦清亮. 海洋磁力测量技术设计关键技术研究[D]. 北京: 中国地质大学, 2015
[7] 赵建虎, 邹亚靖, 吴永亭, 等. 深度约束的海底控制网点坐标确定方法[J]. 哈尔滨工业大学学报, 2016, 48(10): 137–141 DOI: 10.11918/j.issn.0367-6234.2016.10.020
[8] 暴景阳, 翟国君, 许军. 海洋垂直基准及转换的技术途径分析[J]. 武汉大学学报·信息科学版, 2016, 41(1): 52–57
[9] 柯灏. 海洋无缝垂直基准构建理论和方法研究[D]. 武汉: 武汉大学, 2012
[10] 柯灝. 利用调和常数内插的局域无缝深度基准面构建方法[J]. 武汉大学学报·信息科学版, 2014, 39(5): 616–620
[11] 张全德, 范京生. 我国卫星导航定位技术应用及发展[J]. 导航定位学报, 2016, 4(3): 82–88
[12] 阳凡林. 海洋导航定位技术及其应用与展望[J]. 海洋测绘, 2006, 26(1): 71–74
[13] 何广源, 吴迪军, 李剑坤. GPS无验潮多波束水下地形测量技术的分析与应用[J]. 地理空间信息, 2013, 11(2): 155–156
[14] 黄国良, 徐恒, 熊波, 等. 内河无人航道测量船系统设计[J]. 水运工程, 2016, (1): 162–168
[15] 金久才, 张杰, 马毅, 等. 一种无人船水深测量系统及试验[J]. 海洋测绘, 2013, 33(2): 53–56
[16] 孙强. 多波束测深系统声线跟踪方法对比分析[J]. 海洋测绘, 2015, 35(2): 48–51
[17] 赵荻能. 声速剖面精简运算的改进D-P算法及其评估[J]. 测绘学报, 2014, 43(7): 681–689
[18] 王德刚, 叶银灿. CUBE算法及其在多波束数据处理中的应用[J]. 海洋学研究, 2008, 26(2): 82–88
[19] 赵建虎, 张红梅, 严俊, 等. 削弱残余误差对多波束测深综合影响的方法研究[J]. 武汉大学学报·信息科学版, 2013, 38(10): 1 184–1 187
[20] 马小川, 栾振东, 张鑫, 等. 基于ROV的近海底地形测量及其在马努斯盆地热液区的应用[J]. 海洋学报, 2017, 39(3): 76–84
[21] 滕惠忠, 辛宪会, 李军, 等. 卫星遥感水深反演技术的发展与应用[C]. 高分辨率对地观测学术年会, 北京, 2013
[22] 黄文骞, 苏奋振, 杨晓梅, 等. 多光谱遥感水深反演及其水下碍航物探测技术[J]. 海洋测绘, 2015, 35(3): 16–19
[23] 曹斌, 邱振戈, 曹彬才. 四种遥感浅海水深反演算法的比较[J]. 测绘科学技术学报, 2016, 33(4): 388–393
[24] 胡敏章, 李建成, 金涛勇, 等. 联合多源数据确定中国海及周边海底地形模型[J]. 武汉大学学报·信息科学版, 2015, 40(9): 1 266–1 273
[25] 李倩倩, 鲍李峰. 高精度测高重力场反演南海海底地形[J]. 海洋测绘, 2016, 36(2): 1–5
[26] Coiras E, Petillot Y, Lane D M. Multiresolution 3-D Reconstruction from Side-Scan Sonar Images[J]. IEEE Transactions on Image Processing, 2007, 16(2): 382–390 DOI: 10.1109/TIP.2006.888337
[27] 王爱学, 赵建虎, 尚晓东, 等. 单波束水深约束的侧扫声呐图像三维微地形反演[J]. 哈尔滨工程大学学报, 2017, 38(5): 739–745
[28] 赵建虎, 尚晓东, 张红梅. 水深数据约束下的声呐图像海底地形恢复方法[J]. 中国矿业大学学报, 2017, 42(2): 443–448
[29] 元建胜. 海岸带地形一体化测图生产技术研究[J]. 海洋测绘, 2011, 31(4): 75–78
[30] 路晓庆. 我国卫星海洋遥感监测[J]. 海洋预报, 2008, 25(4): 85–89 DOI: 10.11737/j.issn.1003-0239.2008.04.013
[31] 涂晔昕. 使用多源遥感影像监测深圳市海岸线变迁[J]. 海洋开发与管理, 2016, 33(10): 83–88 DOI: 10.3969/j.issn.1005-9857.2016.10.018
[32] 徐进勇. 2000-2012年中国北方海岸线时空变化分析[J]. 地理学报, 2013, 68(5): 651–660 DOI: 10.11821/xb201305007
[33] 李清泉. 海岸带地理环境遥感监测综述[J]. 遥感学报, 2016, 20(5): 1 216–1 229
[34] 秦海明, 王成, 习晓环, 等. 机载激光雷达测深技术与应用研究进展[J]. 遥感技术与应用, 2016, 31(4): 617–624
[35] Mandlburger G, Pfennigbauer M, Pfeifer N. Analyzing Near Water Surface Penetration in Laser Bathymetry—A Case Study at the River Pielach[J]. ISPRS Annals of Photogrammetry, Remote Sensing and Spatial Information Sciences, 2013, 1(2): 175–180
[36] Zhao J Z, Zhao X L, Zhang H M, et al. Shallow Water Measurements Using a Single Green Laser Corrected by Building a Near Water Surface Penetration Model[J]. Remote Sensing, 2017, 9(5): 426–443 DOI: 10.3390/rs9050426
[37] 赵建虎, 王晓, 张红梅, 等. 侧扫声呐图像分割的中性集合与量子粒子群算法[J]. 测绘学报, 2016, 45(8): 935–942 DOI: 10.11947/j.AGCS.2016.20150555
[38] 严俊, 张红梅, 赵建虎, 等. 多波束声呐后向散射数据角度响应模型的改进算法[J]. 测绘学报, 2016, 45(11): 1 301–1 307
[39] 秦海波. 基于GNSS技术的近海岸潮位提取关键技术研究[D]. 南昌: 东华理工大学, 2015
[40] 黄家勇, 张红梅, 陈志高, 等. 流场构建中径向基函数节点选取方法研究[J]. 海洋通报, 2016, 35(5): 516–522 DOI: 10.11840/j.issn.1001-6392.2016.05.005
[41] 何林帮, 冯杰. 基于回波强度的海底底质分类系统关键技术[J]. 大地测量与地球动力学, 2015, 35(1): 140–144
[42] 郭永刚. 海底声参数反演研究与应用[D]. 青岛: 中国海洋大学, 2004
[43] 唐秋华. 济州岛南部海域海底声呐图像分析与声学底质分类[J]. 海洋学报, 2014, 36(7): 133–141
[44] 章皖秋. 地物波谱的自组织神经网络分类[J]. 遥感信息, 2004, (4): 26–28
[45] 罗意, 扬凯, 吴坤泽. 多波束水下三维成像系统在工程检测中的应用[J]. 中国西部科技, 2015, 14(10): 80–82 DOI: 10.3969/j.issn.1671-6396.2015.10.021
[46] 王志云. S101标准与SS57标准电子海图产品初步比较[J]. 航海, 2015, (2): 51–54
[47] 张安民. e-航海中的动态信息服务若干关键技术研究[D]. 武汉: 武汉大学, 2013
[48] 吴文周, 苏奋振, 左秀玲, 等. 基于Web服务的海洋GIS功能设计与原型系统实现[J]. 海洋环境科学, 2015, 34(4): 601–605
[49] 李久松. 基于i4Ocean平台的海上溢油动态可视化及应急决策支持系统的研究与设计[D]. 青岛: 中国海洋大学, 2013