论文部分内容阅读
海洋重力测量技术和航空重力测量技术是获取重力场信息的有效手段。当前,在海空重力测量仪器精度并无明显改善的情况下,对测量数据进行精细化处理研究具有十分重要的理论和实践意义。本文在前人研究基础上,重点研究了重力异常或重力扰动信号提取中FIR(Finite Impulse Response)、IIR(infinite impulse response)和小波阈值滤波的设计方法,以及重力扰动归算的空间改正法、球谐系数法和测线网交叉点快速搜索方法。论文的主要研究内容及成果概括如下:1.基于牛顿第二运动定律,详细推导了海空矢量和标量重力测量的数学模型,分别给出了重力测量系统在当地水平坐标系和载体坐标系下的观测方程;逐一推导了海空重力测量中的各项改正,并进行了数值计算与分析。2.研究了FIR、IIR和小波阈值滤波等海空重力测量数据精细滤波技术,并利用实测数据对各种滤波方法进行了验证及对比分析。主要成果有:(1)基于最佳一致逼近法和频率采样法设计了两种FIR低通滤波器。结果表明:在相同设计指标下,采用最佳一致逼近法设计的滤波器比采用频率采样法设计的滤波器具有更好的低通滤波效果。将设计的低通滤波器应用于海洋实际重力测量数据分析,重力异常精度与窗函数法滤波器结果吻合较好。(2)详细对比分析了不同小波、阈值确定及其施加方式对滤波效果的影响。结果表明:db6、db7、db8、db9、db10、sym6、sym7、sym8、sym9、sym10、coif3、coif4和coif5小波较适用于海洋重力测量应用,小波分解层次可取至8层或9层,采用史坦无偏风险阈值滤波效果较好;小波分解至第8、9层时,阈值滤波结果与经截止频率为0.005Hz、0.0033 Hz Butterworth低通滤波器的滤波结果吻合较好,但小波阈值滤波结果比Butterworth低通滤波器的滤波结果更加平滑,两者滤波结果差值的RMSE在0.25mGal以内。小波阈值滤波较Butterworth低通滤波器更容易分解出噪声成分,可以更有效地消除重力畸变。3.研究了海空重力测量数据精密归算方法。研究探讨了利用地球重力场模型将航线高度处重力扰动观测值归算至平均飞行高度面的方法,以及海面重力异常观测值归算至大地水准面的方法。针对目前航空重力测量空间改正方法的不完备性问题,提出了重力扰动一步归算法-球谐系数法,用Somigliana严密公式对其正确性进行了验证,并以球谐系数法为参考,对比分析了不同空间改正公式对重力扰动归算的影响。结果表明,美国大地测量局(NGS)改进的三阶公式精度最高,H&M(Heiskanen and Moritz,1967)二阶公式次之,我国学者广为采用的二阶公式①精度较低。因此,在1mGal或更好精度的航空重力测量,或者高海拔测量,建议采用球谐系数法或NGS改进的三阶公式进行归算,以提高航空重力测量成果精度。4.研究了海空重力测量测线网内符合精度评估技术。针对目前航空重力测量数据处理中常规交叉点搜索方法存在搜索速度慢、正确率低、适用范围窄的不足,研究分析了测线交叉点快速搜索方法,提出一种以主测线点为搜索中心,采用一定大小的搜索半径对主副测线点进行搜索,将搜索到的主副测线上相邻两点组成线段序列,运用行列式法、投影法、面积法三种判断准则精确求取交叉点及不符值的方法,可称为滑动窗口求解法。通过实际算例,从搜索时间和正确率两方面对滑动窗口求解算法进行了讨论分析。结果表明,本文所提方法简单有效,能明显提高测线交叉点搜索效率,减少工作量。