一种卫星重力梯度数据的ARMA最优滤波模型构建方法

由于卫星重力梯度观测的有色噪声特性和海量观测特征，在利用直接法进行重力场模型的最小二乘求解时，观测值的协方差阵为超大型的非对角阵，这给数值求解带来了极大困难.本文提出了一种基于先验误差功率谱密度的最优ARMA滤波模型构建方法，结合法方程的分块求解策略，可实现对卫星重力梯度观测值的高效滤波处理.数值仿真结果表明，利用最优ARMA滤波器进行时域滤波后，法方程的态性得到了明显改善，重力梯度观测值中的有色噪声得到了有效的"白化"处理，大地水准面精度得到了显著提升.     

利用GOCE模拟观测反演重力场的Torus法

在介绍Torus方法反演地球重力场模型的基本原理和方法的基础上,基于圆环面上均匀分布的卫星引力梯度模拟观测值解算了200阶次的地球重力场模型,在无误差情况下,Torus方法解算模型的阶误差RMS小于10-16,验证了该方法的严密性。利用61dGOCE卫星轨道上无误差的模拟引力梯度观测值解算了200阶次的地球重力场模型,分析了格网化误差、极空白对解算精度的影响,迭代3次后,在不考虑低次系数情况下,模型的大地水准面阶误差和累积误差均较小,最大值仅为0.022mm和0.099mm。在沿轨卫星引力梯度模拟数据中加入5mE/Hz1/2的白噪声,基于Torus方法和空域最小二乘法解算了200阶次的地球重力场模型,Torus方法的精度略低于空域最小二乘法的精度,在不考虑低次项的情况下,两种方法解算模型的大地水准面阶误差最大值分别为1.58cm和1.45cm,累积误差最大值分别为6.37cm和5.55cm。但由于采用了二维快速傅里叶技术和块对角最小二乘法,极大地提高了计算效率。本文数值结果说明Torus方法是一种独立有效的方法,可用于GOCE任务海量卫星引力梯度观测值反演重力场的快速解算。     

联合EGM2008模型重力异常和GOCE观测数据构建超高阶地球重力场模型SGG-UGM-1

本文研究了联合卫星观测数据和重力异常数据确定超高阶重力场模型的理论方法,并使用EGM2008模型重力异常和GOCE（gravity field and ocean circulation explorer）观测数据构建了重力场模型SGG-UGM-1。重点研究了由球面格网重力异常快速构建超高阶重力场模型的块对角最小二乘方法,将OpenMP技术引入到块对角最小二乘中以提高计算效率,并基于模拟数据验证了方法及算法和软件模块的正确性。采用本文制定的联合解算策略,利用GOCE重力卫星观测数据构建的220阶次法方程和EGM2008模型重力异常构建的2159阶次块对角法方程,联合求解了2159阶次的重力场模型SGG-UGM-1。将SGG-UGM-1与EGM2008、EIGEN-6C2、EIGEN-6C4等超高阶模型在频谱域内进行了比较分析,结果表明SGG-UGM-1相对参考模型的系数误差较小,且在220阶次内的系数精度相比EGM2008模型有了提高。采用中国与美国的GPS/水准数据和毛乌素测区的航空重力观测数据对这些模型进行了外符合精度的检验。检核结果表明,在中国区域,SGG-UGM-1模型大地水准面的精度在EIGEN-6C2和EIGEN-6C4两个模型之间,优于GOSG-EGM模型和EGM2008模型,与美国区域几个模型的精度相当。利用毛乌素测区的航空重力数据对几个模型进行了检核,结果表明SGG-UGM-1模型计算的重力扰动精度与EGM2008、EIGEN-6C4模型相当,优于GOSG-EGM模型和EIGEN-6C2模型。     

GRACE与SLR J_2项时间序列的小波分析

利用小波功率谱分析和离散小波变换,对CNES/GRGS、CSR、GFZ等研究机构发布的GRACE时变重力位模型的J2项时间序列进行了周期分析和特征提取,并与利用SLR跟踪数据获得的J2项时间序列的分析结果进行了比较,结果显示:SLR跟踪数据获得的J2项相较GRACEJ2项具有更明显的季节性变化特征.GRACE 时变重力位模型的J2项中混淆进了161天的周期信号,是由海洋潮汐模型中存在的太阳半日分潮S2项的误差引起的.     

利用GRACE与Sentinel-1反演地下水变化与地表沉降研究

针对地表沉降与地下水变化之间的关系研究,助力资源开发与灾害防治之间突出矛盾的解决。该文利用GRACE/GRACE-FO卫星观测数据反演京津冀地区的地下水变化情况,利用Sentinel-1卫星观测结果,通过小基线集干涉测量技术(SBAS-InSAR)计算该区域的地面沉降速率,在地下水变化分析的基础上研究地下水变化与降雨的相关关系,并进而得到地表沉降与地下水变化的相关关系。结果表明,2003—2020年京津冀地区地下水储量呈减少趋势,变化速率达到-11.4 mm/a,降雨量年际变化较大,对地下水储量亏损存在影响;京津冀地区中,北京市地表形变相较于地下水变化滞后时间约为3个月,地下水变化与地表形变存在较高的相关性,二者相关系数达到0.758。     

东秦岭铅锌银金钼多金属成矿带成矿规律及找矿标志

东秦岭位于华北板块与扬子板块之间的拼合带——秦岭造山带东段。该区是一典型的地球化学急变带与地球物理梯度交叉区，壳幔富含Pb，Zn，Au，Ag，W，Mo元素，为一元古宙一古生代裂陷槽，熊耳群、宽坪群、二郎坪群、耀岭河组火山岩系中的火山喷发Cu，Pb，Zn，Ag，Au矿(化)层；秦岭群、官道口群、栾川群、陶湾群滨海、浅海陆缘碎屑碳酸盐建造中强烈的火山喷(气)发活动形成的硅质条带(团块)硅质岩以及含Ag，Pb，Zn，Au矿化层状矽卡岩带，多期次的构造岩浆活动，使区内Pb，Zn，Ag，Au，W，Mo叠加富集体成矿，从而使东秦岭地区Pb，Zn，Ag，Au，W，Mo多金属矿床在区域分布、成矿空间、时间上呈现出一定的规律性和明显的找矿标志。     

多源卫星重力场模型精度与频谱分析

针对新一代卫星重力探测技术对地球重力场的频谱贡献问题,该文提出了一种基于GPS/水准数据获取多源卫星重力场模型频谱变化特征的方法。采用GPS/水准外符合检验,有效分析评估了多源卫星重力场模型在中国东、西部地区的精度水平。研究结果表明,以CHAMP、GRACE和GOCE卫星为代表的高-低卫星跟踪卫星、低-低卫星跟踪卫星和卫星重力梯度技术,对地球重力场的频谱贡献分别集中在600km以上的长波和中长波、300km以上的中波、200~350km之间的中短波部分。     

基于地面重力的卫星重力梯度检校方法

卫星重力梯度仪在轨检校是提高梯度模式重力卫星观测质量的关键.本文面向中国未来梯度模式重力卫星规划任务,研究提出一种基于地面重力的卫星在轨检校方法,该方法顾及卫星设计指标,从地面先验重力数据的精度、空间尺度以及卫星重力梯度仪的观测噪声等卫星检校要素开展分析研究,成功实现地面数据与卫星观测数据的比对检校.研究结果表明,在12°×12°的地面检校区域内,利用精度优于1 mGal的地面重力观测数据即可实现噪声低于10 mE的重力梯度仪高精度在轨检校.     

高分七号卫星双波束激光测高仪在轨几何检校与试验验证

高分七号卫星装载了国内首台双波束对地观测激光测高仪,用于辅助两线阵光学相机实现1:1万立体测图.由于卫星发射时的振动及入轨后空间环境变化,激光测高仪在轨指向等参数相对于实验室标定值发生较大偏移,直接影响激光测高精度.针对高分七号双波束激光测高仪,本文提出了由粗到精的"两步法"在轨几何检校方案.首先,构建单波束激光在轨几何检校模型,基于特殊波形分析来估算激光足印落点位置信息,实现单波束激光指向粗检校;然后,以单波束激光在轨几何检校模型为基础,构建双波束激光在轨几何检校模型,在大气延迟改正、潮汐改正等基础上,采用地面探测器阵列捕获光斑,实现双波束激光指向测距联合检校.最后,采用平静湖面上的激光测高数据和地面控制数据,分别验证激光检校后相对和绝对高程测量精度.试验结果表明:高分七号卫星激光高程测量相对精度优于0.06 m(1σ),平坦地区的绝对测高精度达到0.10m(1σ).