海底地形探测和模型研制现状及精度分析

海底地形影响着海洋自然灾害、资源开发以及军事国防.海底地形的探测手段由过去传统的船载声呐技术发展到依靠重力数据反演水深.国外从1980年代就开始利用船载和卫星测高数据开发全球水深模型,截至目前发布的水深模型空间分辨率达到15″.本文总结了声呐、激光测深雷达、卫星测高重力数据和遥感影像反演水深的研究现状,并分析每种测深手...     

基于重力梯度的模拟退火法反演中国南海海底地形

模拟退火法是一种全局寻优算法,在地球物理反演中得到广泛应用.本文应用快速模拟退火法对中国南海海域的重力垂直梯度进行反演来获得中国南海高分辨率海底地形,首先推导出密度界面与重力垂直梯度的对应关系,然后利用理论模型试验证明快速模拟退火法能成功地根据重力垂直梯度来计算密度界面的变化,最后将其应用于由Geosat,ERS-1/2,T/P,Jason-1,EnviSat-1卫星测高数据计算得到的中国南海海域重力垂直梯度异常,反演得到的中国南海海底地形与LDEO船测海深相接近,两者之间差的最大值为0.24 km.     

顾及海底地形坡度的南海海底地形分区GGM反演

利用重力地质法(Gravity-Geologic Method, GGM)反演海底地形时,海水与海底洋壳的密度差异常数是影响反演精度的一项关键参数.由于海底地形复杂程度的影响,不同区域密度差异常数存在差异.针对于此,本文以南海局部海域(113°E—119°E, 12°N—19°N)为实验区,依据区域内坡度及地形分布,将研究区域划分为6个子区域,分别求取每个子区域的最佳密度差异常数,以改善海底地形反演的精度.基于HY-2A测高数据获得的重力异常数据,反演了研究区域格网间隔1′×1′的SGGM海底地形模型.结果表明,将SGGM模型与检核点实测水深进行比较,其差值的标准差为101.46 m,相对精度为2.82%,优于GEBCO＿2021模型、V19.1模型对比实测水深差值的标准差(131.50 m、129.81 m)与相对精度(3.64%、3.59%).结合各子区域不同的坡度及船测点分布,统计各子区域相对精度为2%～4%,并分析得出反演较差结果集中分布在地形突变区域,而坡度平缓区域的标准差可达到42.20 m.此外,同未进行区域划分而采用一个密度差异常数反演的结果相比,反演精度提高了9.68...     

计算测高卫星的正常点及其平均海面高

利用测高卫星的冻结轨道特征计算了正常点轨迹,利用三次B样条函数和Stacking技术内插计算了T／P、ERS-1／2和Geosat测高数据在正常点的平均海面高,改善了测高数据的精度和分辨率．最后,利用正常点上的平均海面高反演了多种分辨率的海域重力异常,其中所反演的海域重力异常的标准差为mGal级．     

测高数据反演海底地形研究进展与展望

针对卫星测高技术发展研究现状,重点突出海洋重力场建模和海底地形反演技术,分析两者相关关系,介绍了国际上公开发布的两款海洋重力模型,阐述了国内相关技术发展现状;对测高数据反演海底地形的重力地质方法、S&S法、导纳函数法、垂直重力梯度异常的解析算法研究进展情况进行了梳理,给出了算法基本原理和数理基础.     

应用多种来源重力异常编制中国海陆及邻区空间重力异常图及重力场解读

本文通过分析陆地实测空间重力异常数据、海洋船载测量空间重力异常数据、卫星测高重力异常,布格重力异常数据、EGM2008地球重力模型数据等多种来源数据的性质和精度,并对相关数据进行对比,研究了编制1:500万中国海陆空间重力异常图的数据使用方案和技术方法.在地形较为平坦、实测数据分布均匀的陆区,使用实测数据,在地形复杂,实测数据稀少以及没有实测数据的陆区或岛屿,利用布格重力异常反推空间异常的方法合成平均空间重力数据,西藏地区的数据对比实验证明合成平均空间重力异常数据是一种有效的数据补充.利用三观测列方差分解法在南海地区对船载测量空间重力数据和美国SS系列及丹麦DNSC08GRA卫星重力数据进行了方差分解计算,结果表明不同来源的卫星测高重力数据具有很大的一致性,数据精度较以往有了很大的提高.海区空间重力数据使用原则是在船载重力测量数据校准下,全面使用卫星测高重力数据进行编图.海陆过渡区的异常处理应以EGM2008地球重力模型重力场为基准参考场,实现海陆异常平缓过渡,无缝连接.对中国海陆空间重力异常场进行了小波变换处理,对空间重力异常场进行了解读,勾画出三横四竖的一级重力梯级带及其所围限的8个一级重力异常区,并划分了二级重力异常区和梯级带,为块体构造学体系中大地构造格架的建立提供了地球物理证据.     

基于重力地质法的南中国海海底地形反演

根据重力地质法（GGM）,利用南中国海海域内63179个船测控制点水深将测高自由空间重力异常划分为长波参考场和短波残差场,并反演出了该海域112°E-119°E,12°N-20°N范围的1’×1’海底地形模型,该过程中使用的海水和海底洋壳密度差异常数1.32 g·cm^-3通过实测水深估计得到.利用反演得到的GGM模型对剩余的10529个检核点船测水深插值计算后与实测水深进行比较,其较差结果的均值为-1.64 m,标准差为76.95 m,相对精度为4.06%.此外,根据船测点数量、分布和海底地形的不同,选择了三个海域进行统计,结果表明：在船测控制点分布均匀的海域,GGM模型精度优于ETOPO1模型,在控制点过于分散的海域其精度会有所下降,但好于船测水深的直接格网化结果.为进一步探究检核点的较差结果中出现较大数值的成因,本文对精度较差的点位进行了单独分析,选择了两条船测航迹剖面进行了研究,并分析了检核点的水深较差、相对精度与水深和重力异常的关系,结果表明：GGM模型精度受水深和重力异常的相关性影响较小,受海底地形复杂程度影响较大,地形坡度变化平缓海域的预测精度明显高于海山地区.最后,综合GGM模型和ETOPO1模型优势,利用所有船测水深作为控制,生成了综合的海底地形模型.     

单层/双层地壳物理参数对重力数据反演海底地形影响分析

针对利用海面重力数据构建海底地形模型质量受诸多地壳地球物理参数影响现状,以重力异常反演海底地形的频率域算法模型为基础,讨论了不同地球物理参数对单层/双层地壳仿真结构和不同均衡补偿模式下导纳函数各频段影响程度,推导给出了不同地球物理参数及其组合对海底地形反演结果影响的定量数学表达.仿真试验结果表明,单层地壳的密度变化导致单层导纳函数和双层导纳函数明显分离;双层地壳仿真结构中海山载荷密度对导纳函数影响显著且几乎为全频段影响;因此,建议应用解析法开展数值试验分析时,重点考虑海山载荷密度(双层地壳)或者地壳密度(单层地壳)变化.就当前海底地形反演精度而言,运用相关解析算法建立海底地形数值模型时,应重点关注有效弹性厚度、平均海深、海山和海水密度差异参数取值准确性,特别是顾及地壳均衡补偿恢复海底地形,岩石圈有效弹性厚度(通过有效弹性厚度可获取挠曲刚度)取值准确程度将较大程度影响最终海底地形模型效能高低.     

我国海域大地水准面的计算及其与大陆大地水准面拼接的研究和实施

研究和实施了由卫星测高数据计算垂线偏差，用莫洛 金斯基（Molodensky）公式反演 大地水准面高，由此求得我国海域大地水准面高. 为了检核，将测高垂线偏差利用逆维宁迈 纳斯（Vening Meinesz）公式反演重力异常，与海上船测重力值进行了外部检核；同时还用 司托克斯（Stokes）公式，将上述反演的重力异常计算大地水准面高，与莫洛金斯基公式直 接解得的相应结果进行比较作为内部检核. 在积分计算中充分应用了FFT的严格公式.由重力和GPS水准数据确定的陆地大地水准面，和主要由卫星测高数据确定的海洋大地水准 面，二者之间一般都存在以系统误差为主的拼接差，本文分析了产生这一现象的主要原因, 并结合我国在陆海大地水准面拼接区重力资料稀疏的实际，提出了新的拼接技术，最后将拟 合参数校正中国全部海域的重 力大地水准面，以最大限度地削弱拼接点和制约测高海洋大地水准面可能存在的系统误差.     

沿海雷达卫星测高波形重定多子波参数方法和重力异常恢复

雷达卫星测高数据质量是测高在大地测量学、地球物理学和海洋学等领域应用研究的关键。在沿海海域,测高波形受到地形和地球物理环境等污染,测高数据精度低于开阔海域的测高数据。提出了一种新的测高波形重定方法,即多子波参数重定方法（MSPR）,用于改善沿海海域的测高数据质量。由测高数据,利用最小二乘配置法解算沿海重力异常,给出了对应的方差-协方差矩阵。利用MSPR对台湾岛周边海域的Geosat／GM波形数据进行重定,与台湾岛大地水准面进行比较,MSPR重定结果精度要优于常用的β-5参数方法重定结果和GDR数据。利用重定后的测高数据,由最小二乘法计算了台湾岛近海海域的重力异常。与船测重力数据进行比较,由重定数据得到的结果精度优于GDR结果。     

中国海及邻域重力场特征及其构造解释

利用中国海及邻域 2’× 2’由Seaset,Geosat,ERS 1及TOPEX/POSEIDON等卫星测高资料解算的自由空气重力异常 ,结合ETOPO5全球海底地形资料计算了中国海及邻域布格重力异常 ,并反演了得到了中国海及邻域地壳厚度 .经对重力异常特征进行分析 ,讨论了重力异常与大地构造及其活动的相关性 .     

基于泊松小波径向基函数融合多源数据的局部重力场建模

融合多源数据的高精度、高分辨率的局部重力场建模是物理大地测量学的前沿和热点问题.本文研究了基于径向基函数融合多源数据的局部重力场建模方法,利用Monte-Carlo方差分量估计实现了不同类型的观测数据的合理定权,引入了最小标准差法确定基函数的适宜网络,分析了地形因素对于基函数网络确定及局部重力场建模精度的影响.以泊松小波基函数为构造基函数,结合残差地形模型,融合实测的陆地重力异常、船载重力异常及航空重力扰动数据构建了局部区域陆海统一的似大地水准面模型.研究结果表明:引入残差地形模型平滑了地形质量引入的高频扰动信号,简化了基函数的网络设计;并提高了重力似大地水准面的精度,平原地区其精度提高了4mm,地形起伏较大的山区其精度提高了约5cm.总体而言,基于"三步法"构建的局部重力似大地水准面在荷兰、比利时及德国相关区域,其精度分别达到1.12cm、2.80cm以及2.92cm.     

联合多代卫星测高和多源重力数据的局部大地水准面精化方法

本文研究了基于泊松小波径向基函数融合多代卫星测高及多源重力数据精化大地水准面模型的方法.分别以沿轨垂线偏差和大地水准面高高差作为卫星测高观测量,研究了使用不同类型测高数据对于大地水准面建模精度的影响.针对全球潮汐模型在浅水区域及部分开阔海域精度较低的问题,引入局部潮汐模型研究了不同潮汐模型对于大地水准面的影响.数值分析表明:相比于使用沿轨垂线偏差作为测高观测量,基于沿轨大地水准面高高差解算得到的大地水准面模型的精度更高,特别是在海域区域,其精度提高了2.3cm.由于使用沿轨大地水准面高高差作为测高观测量削弱了潮汐模型长波误差的影响,采用不同潮汐模型对大地水准面解算的影响较小.总体而言,船载重力及测高观测数据在海洋重力场的确定中呈现互补性关系,联合两类重力场观测量可以提高局部重力场的建模精度.     

基于轨道扰动引力谱分析的方法确定低低观测重力卫星反演重力场的空间分辨率

基于低低卫-卫跟踪重力卫星的轨道特性,从垂直和水平两个方向计算了重力卫星高空扰动引力,并根据其谱特性及星载加速度的测量噪声水平分析了重力卫星能反演重力场的阶数.利用EGM96重力场模型分别计算了400 km、450 km和500 km 轨道高度处重力卫星受到的扰动引力谱及扰动引力谱的平均量级,分析其垂直特性表明:在三个轨道高度处能分别能反演150、140和130阶的重力场模型.利用两颗同轨重力卫星相距220 km的特性,计算了400 km、450 km和500 km 轨道高度处纬度相差2°的两颗卫星纬向扰动引力差,即扰动引力水平分量,分析其谱特性,表明:重力卫星能反演至117阶的地球重力场模型.     

基于Lp范数正则化的V字型密度界面重力反演

V字型密度界面是一类常见的密度界面,如海沟、半地堑以及俯冲带之下的莫霍面,利用重力数据刻画此类密度界面形态对于区域构造研究、油气勘探以及物理海洋学等都具有重要意义.本文首先建立了L_p-范数形式的模型约束函数,并利用正则化原理将其与重力数据误差函数和已知深度约束函数结合形成V字型密度界面反演的目标函数,推导了目标函数的梯度表达式,并以非线性共轭梯度法为核心给出了反演流程.二维简单模型试算结果表明p=5时该方法能准确地刻画V字型密度界面起伏特征,且亦能准确地应用于二维复杂密度界面和三维界面的反演.最后将反演方法应用于挑战者深渊及邻区的实际资料处理之中,利用研究区海底地形数据和沉积层厚度数据对自由空间重力异常逐层剥离而得到莫霍面引起的重力异常,用本文方法对此重力异常进行反演,结果呈现了板块俯冲作用引起的V字型莫霍面起伏特征.     

基于卫星测高的海域大地水准面

利用测高数据的一次差分计算海域垂线偏差，有效降低了动力海面地形和系统残差对垂线偏差的影响；然后根据扰动场元间的协方差函数是具有各态历经性的平稳随机函数这一特征，提出了利用垂线偏差精确逼近海域大地水准面的协方差函数. 而海域大地水准面的精确确定，为从测高数据中精确分离动力海面地形提供了条件. 本文还利用Topex/Poseidon、ERS 1/2测高数据计算了全球海域大地水准面和动力海面地形，证明了本文所述方法是科学合理的.     

川滇地区陆地流动重力测网场源分辨能力评估

陆地重力测量是研究地球内部物质运移和质量变化的重要手段,场源分辨能力是评估流动重力测网监测效能和由场及源应用研究的关键指标。本文以构造运动和地震活动活跃的川滇地区为例,在设定场源异常体参数下,通过Tesseroid网格模型建立不同分辨率的重力异常扰动模型检测板,并基于川滇地区陆地流动重力测网测点分布,对比反距离加权法、 Krigng插值法和最小二乘配置法对模型异常检测场的恢复效果。在此基础上,基于不同分辨率的模型检测板和模型抗噪分析结果,研究川滇地区陆地流动重力测网的空间分辨率和场源反演分辨能力,进而评估该地区重力测网的场源分辨能力,并对该地区2014至2017年测网实测数据开展应用研究。结果表明,最小二乘配置法对川滇地区重力场恢复效果最佳,川滇地区陆地流动重力测网大部分地区的空间分辨率优于0.75°×0.75°,部分地区可达0.25°×0.25°;在Ⅰ级活动地块边界带为0.50°×0.50°,在Ⅱ级活动地块边界带优于0.75°×0.75°;进一步根据活动地块边界和重力测点分布,对川滇地区分区开展陆地流动测网场源分辨能力评估。     

由多卫星测高大地水准面推断西太平洋边缘海构造动力格局

由于卫星观测技术的发展，现在已能利用测高大地水准面简捷可靠地研究海底构造与动力问题. 根据Geosat T2/ERM、Topex/Poseidon 和ERS 1/2测高数据给定的0°N ～45°N、100°E～150°E范围内4′×4′大地水准面，采用全球地形/均衡的重力位效应改正，计算布格、格莱尼和均衡大地水准面. 由格莱尼大地水准面反演Moho面埋深，再从均衡大地水准面起伏推算小尺度地幔流应力场. 结果表明，菲律宾海和南海显示了与磁条带、扩张脊对应或斜交的高频成分大地水准面起伏条带. 各边缘海盆的Moho面埋深有往南变浅的趋势，与菲律宾海各海盆的Moho面埋深大致相当，说明琉球—台湾—菲律宾岛弧两侧的构造动力强度基本相近. 从各种构造特征和大、中、小尺度的地幔流应力场的驱动机制，可以证明岛弧内侧的边缘海具有不同于大洋、大陆的独特构造动力格局和特征.     

EGM2008地球重力模型数据在中国大陆地区的精度分析

本文介绍了5′×5 ′的EGM2008地球重力模型及其在全球的精度评价.按照地形变化规律,将中国大陆大致分为7个区域,在10 km网度上,将EGM2008地球重力模型数据与中国地面实测空间重力网格数据进行了对比.由于数据源的问题,中国大陆的模型数据精度普遍低于北美和欧洲.二种数据在地形平坦的东部地区差别较小,向西随着地形复杂程度的增加,二种数据之间的标准差从小于10 mGal增大到50多mGal.畸变点分析表明精度极低的网格点均分布在地形起伏大的地区.总体而言,5′×5′的EGM2008地球重力模型数据在中国大陆将近80％的面积上的精度可达10 mGal之内,可用于小比例尺重力编图和构造研究.在地形起伏较大的中国西部以青藏高原为例进一步比较了EGM2008重力模型和重力测点数据,结果表明在重力点分布稀疏不均匀的地区,平面网格数据难以准确表达重力场信息.由于缺少地面重力数据控制,EGM2008重力模型数据在中国西部精度较低,但模型数据依然在很大程度上提高了空间重力异常信息的丰富程度.将中国区域重力调查成果数据应用于地球模型的构建是一项有意义的工作.     

用GRACE卫星跟踪数据反演地球重力场

利用141天GRACE卫星观测资料,包括K波段、星载加速度和卫星轨道数据,反演了80阶地球重力场模型IGGGRACE01S,该模型在半波长为500km的空间分辨率上,确定大地水准面的精度约为0012m,中长波(<80阶)精度优于重力卫星发射以前研制的重力场模型. 与EIGEN_GRACE02S、EIGEN_CHAMP03S和EGM96模型的位系数相比,该模型系数最接近于EIGEN_GRACE02S,与另两个模型差异较大. 比较几种模型确定的全球重力异常和大地水准面起伏,结果发现IGGGRACE01S与EIGEN_GRACE02S模型的计算结果比较接近,与EGM96模型结果差异较大,差别较大地区主要在南极等地区. 对于中国大陆,比较IGGGRACE01S模型（前72阶）计算的重力异常和NIMA重力异常数据（25°×25°网格）,两者之间的标准偏差为48mGal.