联合多代测高数据建立中国海域重力异常模型

本文联合T/P数据、T/P新轨道数据、ERS数据、GFO数据、GeosatGM数据和ERS-1/168数据,用测高卫星记录点的位置信息直接计算沿轨大地水准面的方向导数,结合测线轨迹方向的方位角在交叉点处推求垂线偏差,然后利用逆Vening-Meinesz公式计算了中国近海(0o～41oN,105o～132oN)2′×2′格网分辨率的海域重力异常模型。将其与CLS_SHOW99重力异常模型比较,统计结果表示与该模型差异的RMS为8.15mgal,在剔除差值大于20mgal的点(剔除3.3%)以后,RMS为4.72mgal;与某海区船测重力异常比较的RMS为8.91mgal。     

联合TOPEX/Poseidon,ERS2和Geosat卫星测高资料确定中国近海重力异常

处理了 TOPEX/Poseidon(第 9周期至第 2 4 9周期 ) ,ERS2 (第 0周期至第 44周期 )和Geosat/GM(第 1周期至第 2 5周期 )以及 Geosat ERM(第 1周期至第 66周期 )卫星测高资料 ,求解了各自卫星任务的交叉点和垂线偏差 ,利用逆 Vening- Meinesz公式确定了 2 .5′×2 .5′中国近海海洋重力异常 ,并与我国南海船测重力异常作了比较 ,其精度为± 9.3m Gal( 1 Gal=1 cm/s2 )。本文同时导出了严密的 2维平面卷积公式 ,它与 1维严密卷积公式计算结果差值的标准差为± 0 .1 m Gal,而 2维球面公式为± 0 .5 m Gal     

自主海洋测高卫星串飞模式的设计与重力场反演精度分析

根据自主海洋测高卫星发展需求,设计了双星串飞运行模式,该运行模式下2.3 a时间可满足全球海洋区域1'×1'分辨率的地面轨迹覆盖要求。首先,将测高卫星重力场反演分为不考虑轨道运行特点（思路1）和考虑串飞轨道运行特点（思路2）两种思路,利用逆Vening-Meinesz方法开展了正态分布下随机误差传播的仿真计算,获得了两种思路下对应的误差指标。以该误差指标为基础,分别计算了双星串飞模式下两种重力场反演思路对应的精度指标。其中,反演思路2充分利用了串飞模式双星东西方向地面观测值可以进行相对定轨的特点,并考虑到近距离条件下传播误差、地球物理改正误差的系统误差特性,因此反演思路2的垂线偏差精度较反演思路1有了一定的提高,其重力场反演也具有一定的优势。理论计算结果表明,利用思路1的反演方法,2.3 a时间可获得1'×1'重力异常精度为6~10 mGal,4.6 a时间可达到4.2~7.1 mGal;利用思路2的反演方法,2.3 a时间可获得1'×1'重力异常精度为3.9 mGal,4.6 a时间可达到2.8 mGal。     

利用卫星测高数据反演海洋重力异常研究

全面研究了利用卫得测高数据反演海洋重力异常3种主要方法（即Stokes数据解析反解以及逆Vening-Meinesz公式）的技术特点，建立了3种算法的数学模型及其谱计算式，在以1440阶次位模型定义的标准场中完成了3种算法的数值比较和内部检核，通过仿真试验实现了3种算法的可靠性和稳定性检验，最后，本文利用卫得测高实测对南中国海地区的海洋重力异常进行了实际反演，并将反演结果同船测数据进行了比较。     

高度计测距精度对沿轨迹重力异常反演的影响

应用求解沿轨迹重力异常的垂线偏差法以及求解空间分辨率的交叉谱分析法,建立了高度计测距精度与沿轨迹重力异常反演精度以及空间分辨率的关联性模型。首先依据卫星测高原理,给出了沿轨迹重力异常的误差传播公式,然后以此为基础通过推导交叉谱分析中一致性系数与信噪比的数学表达式,建立了高度计测距精度与空间分辨率的解析关系。数值仿真结果表明:雷达高度计测距精度与沿轨迹重力异常反演精度成正比关系,与空间分辨率成幂函数关系,即高度计测距精度提高m倍,沿轨迹重力异常反演精度提高m倍,全球海域平均空间分辨率提高m0.464 4倍。将数值仿真结果与相关文献中对实际测高数据的处理结果进行比较,验证了理论分析及模型的正确性。     

联合卫星重力和卫星测高数据确定稳态海洋动力地形

新一代卫星重力探测任务GRACE大大提高了地球重力场模型中长波分量的精度,使得联合卫星测高平均海面分离更精细稳态海洋动力地形成为可能。本文利用T/P(1994年~2003年)和JASON-1(2003年~2005年)卫星测高数据确定了全球30′×30′平均海面高;基于重力场模型WHU-GM-05,计算得到对应于海面高的GRACE海洋大地水准面格网值;利用“移去-恢复法”和高斯滤波求得全球稳态海面地形。与EGM96、R io05、ECCO和GGM02模型进行比较,检验结果表明GRACE任务有效的改善了海洋大地水准面的精度,使得稳态海洋动力地形能够呈现更多细部。     

卫星测高反演海洋重力异常的精度分析

针对卫星测高技术中由于大地水准面取值受各项误差影响导致精度较低的问题,该文联合多源多代卫星测高数据,基于逆Vening-Meinesz公式确定海洋重力异常,进一步对海洋重力异常进行内部和外部检核。结果表明,卫星测高反演的海洋重力异常与EGM2008比较的精度为±7.116mgal;与船测重力异常比较的精度为7.417mgal,这与国际上对测高重力异常与船测重力异常比较精度一致。     

联合多种测高数据确定中国海及其邻域1．5′×1．5′重力异常

联合多种测高资料和Geosat/GM波形重构数据,基于EIGEN_CG01C重力场模型,采用沿轨迹加权最小二乘方法和逆Vening-Meinesz公式,确定了中国海及其邻域1.5′×1.5′重力异常。将计算结果与最新船测资料进行了比较,标准差为3.37×10-5m.s-2。     

西太平洋海域卫星测高重力垂直梯度分布

从地球重力场基本理论出发,推导了重力垂直梯度计算的数值积分公式,并利用由多颗测高卫星联合反演的2'×2'海洋重力异常资料,解算了西太平洋海域重力垂直梯度.将计算所得的重力异常垂直梯度与现有资料进行比较,其差值均方根为±10.08E,表明两者精度相当.此外,还对差值大小和空间分布进行了深入分析.     

卫星测高沿轨海面高的精细结构及其应用技术初探

首先介绍了频域和多尺度域中的两种沿轨测高海面高精细结构，进而提出了利用沿轨测高海面高精细结构计算平均海面高与海平面距平的新方法，最后进一步讨论了沿轨平均海面高和海平面距平的频谱和多尺度分解，以及它们随时间的变化特性。     

利用重力异常反演马里亚纳海沟海底地形

为了尽可能利用海面重力数据所含地形相关信息, 提升依据重力数据反演海底地形结果质量, 选择南中国海为试验海区, 详细研讨了频率域不同阶次海底地形正演恢复海面重力异常和重力异常垂直梯度在不同海底地形环境下的影响。数值分析试验结果表明, 一次项海底地形正演结果对重力信息贡献起主要作用, 随着海底地形阶次的增加, 相应阶次海底地形正演海面重力信息的幅度不断减弱; 海底地形起伏剧烈且起伏幅度较大海区, 高阶次海面重力异常和重力异常垂直梯度变化明显, 而地形平坦海区正演的海面重力信息变化微弱; 重力异常垂直梯度相比重力异常对海底地形高频部分更加敏感, 在试验海区依据海面重力数据反演海底地形过程中, 若输入数据源为重力异常, 建议顾及至二次项海底地形, 若输入数据源为重力异常垂直梯度, 建议顾及至三次项或者四次项海底地形, 从而尽可能利用海面重力数据所包含信息。     

CHAMP重力卫星结果在中国海及邻近海域的初步分析

利用世界上第一个采用高低卫_卫跟踪技术的CHAMP重力卫星计划导出的全新的高精度全球长波重力场模型EIGEN_1S结果 ,根据二维高斯滤波原理 ,基于相同空域尺度对卫星重力结果和由卫星测高解算的几种版本的海洋重力异常数据进行了长波部分的分析研究。结果表明 ,在中国海及其邻近海域卫星重力数据与卫星测高解算的海洋重力异常数据之间存在明显的偏差。     

利用卫星测高资料反演中国近海海洋重力

本文利用Topex/Poseidon卫星测高资料,从快速Hartley变换(FHT)基本概念入手,给出了Hotine公式在平面近似、球面近似、Molodenskii近似下,反演中国近海海洋重力的数学模型。另对FHT处理中所需的坐标转换以及边缘效应等问题进行了讨论。同时,为了改善长波特性的重力场信息,引入了M阶次的OSU91A参考重力场对上述Molodenskii模型进行了改化。     

多代卫星测高数据联合平差及重力场反演

为了解决多代卫星测高数据之间的不协调问题,基于误差验后补偿理论,提出了将传统的交叉点平差整体解法简化为两步处理法,即首先使用条件平差法对交叉点观测方程进行平差计算,然后沿测高轨迹进行海面高滤波和推估。使用新方法可大大简化联合平差的计算过程,且有利于提高计算结果的稳定性和可靠性,对规则和不规则的区域网和全球网交叉点平差问题都同样适用。实际算例证明了新方法的有效性。     

Determination of Geoid over South China Sea and Philippine Sea from Multi-satellite Altimetry Data

A new computational procedure for derivation of marine geoid on a 2.5′×2.5′grid in a non-tidal system over the South China Sea and the Philippine Sea from multi-satellite altimeter sea surface heights is discussed. Single-and dual-satellite crossovers were performed, and components of deflections of the vertical were determined at the crossover positions using Sand-well's computational theory, and gridded onto a 2.5′×2.5′resolution grid by employing the Shepard's interpolation procedure. 2.5′×2.5′grid of EGM96-derived components of deflections of the vertical and geoid heights were then used as reference global geopotential model quantities in a remove-restore procedure to implement the Molodensky-like formula via 1D-FFT technique to predict the geoid heights over the South China Sea and the Philippine Sea from the gridded altimeter-derived components of deflec-tions of the vertical. Statistical comparisons between the altimeter-and the EGM96- derived geoid heights showed that there was a root-mean-square agreement of ±0.35 m between them in a region of less tectonically active geological structures. However, over areas of tectonically active structures such as the Philippine trench, differences of about -19.9 m were obtained.