小波变换在船载海洋重力测量数据处理与分析中的应用_刘备.pdf

第 卷第期 年月系统工程与电子技术 文章编号：（）网址：收稿日期：；修回日期：；网络优先出版日期：。网络优先出版地址：基金项目：国家自然科学基金（，）；海军工程大学自主立项项目（）；广西空间信息与测绘重点实验室开放基金（）资助课题通讯作者引用格式：刘备，边少锋，纪兵，等小波变换在船载海洋重力测量数据处理与分析中的应用系统工程与电子技术，（）：，（）：小波变换在船载海洋重力测量数据处理与分析中的应用刘备，边少锋，纪兵，贤鹏飞（海军工程大学电气工程学院，湖北 武汉 ；广西空间信息与测绘重点实验室，广西 桂林 ）摘要：船载重力测量是获取海洋地球重力场的重要方法之一，为提高船载重力测量数据的精度，采用 小波和 小波两种小波变换方法，分别对添加了两种随机噪声的某海域重力测量数据进行降噪处理与分析。结果显示，经 小波处理的数据绝对精度分别提高了 和 ，相对精度分别提高了 和；采用 小波变换后，绝对精度分别提高了 和 ，相对精度分别提高了 和。两种小波变换方法处理后的数据均能够很好地保留数据基本变化趋势、为重力勘探判别物体提供判别依据。关键词：船载重力测量；多贝西小波；小波；误差比率中图分类号：文献标志码：，（.，；.，）：，：；（）；引言海洋重力测量作为地球重力场研究的一个重要组成部分，为海洋资源勘探、军事应用、基础海洋研究等领域提供了重要的数据支撑 。目前，海洋重力测量的主要方法包括卫星重力测量、航空重力测量和船载重力测量，船载海洋重力测量是获取高精度以及高分辨率海洋重力场信息的重要方法之一 。常规的重力测量包括静态重力测量和动态重力测量。但由于海洋环境的特殊性，静态重力测量无法满足当前海洋开发和未来海战的需求。动态重力测量以第期刘备等：小波变换在船载海洋重力测量数据处理与分析中的应用 其能够快速获取高精度、高分辨率的优势成为获取海洋重力场信息的主要手段之一 。动态重力测量获取的高精度、高分辨率海洋重力场信息能够为海洋开发、海底地形反演、未来海战等提供一些关键信息 。卫星重力测量作为动态重力测量的主要手段之一，其空间分辨率目前可以达到 ，测量精度可达到 ，若需要更高精度的海洋重力场测量成果，则需要通过船载重力测量方式获取重力场数据对卫星重力场数据进行补充 。船载重力测量由于其测量环境的特殊性，其数据处理成为制约测量成果精度的关键问题之一。对于高精度海洋重力测量数据处理与分析，文献 将 重力仪搭载在不同测量船只上进行测量任务，并对测量结果进行处理分析，搭载在两艘测量船上的测量精度分别为 、，并利用测量结果给出了波罗的海地区局部大地水准面构建的初步结果。等 对北大西洋地区的重力测量数据进行了处理分析，结果为交叉点处的标准差为 ，经平差处理后提高到 。蔡劭琨等 将经验模态分解用于动态重力测量数据处理当中，将重 力 异 常 的 内 符 合 精 度 从 提 高 到 了 。董庆亮等 对丹麦技术大学重力数据用于检查海洋重力测量成果的可行性进行了研究，对海洋重力测量数据的检核方法起到了补充作用。吴燕雄等以激光干涉绝对重力仪为基础设计了一套船载绝对重力仪测量系统，并给出了该系统工作的误差修正方法、动态限制条 件 和 修 正 精 度，经 修 正 后 该 系 统 测 量 精 度 优 于 。黄谟涛等 对海空重力测量误差进行深入分析，提出了一种适用于补偿各类海空重力仪动态效应剩余影响的通用模型，经验证，该模型将海洋重力测量的内符合精度从 提高到 。对于船载海洋重力测量来说，外界环境会产生大量的噪声。如何减弱噪声、消除噪声对测量成果的影响，提高测量精度，具有非常重要的意义 。目前常用的降噪方法是采用有限冲激响应（，）低通滤波器，可以降低高频噪声，但滤波器最优参数的普适性不高 。孙鹤泉等 将最大重叠离散小波变换（，）用于海洋重力观测航行数据的降噪中，并对该方法的可靠程度进行了验证。基于小波变换具有适用性高的优点，本文将小波变换方法引入到船载重力场数据处理分析中，并利用实测数据验证小波变换方法用于船载重力场数据降噪中的精度和可靠性。船载重力测量本文用于验证实验的船载海洋重力场数据是将 型船载重力仪搭载到广海局“海洋四号”测量船上在某海域开展重力实验测量任务所测得的重力场数据。型船载重力仪（见图）是一款针对船载平台的动态重力仪，系统由重力仪本体、不间断电源（，）、全球导航卫星系统（，）接收机、显控设备实时显控软件等组成。其中，重力传感器采用高精度石英挠性加速度计，姿态测量传感器采用高精度光学陀螺。重力仪在线实时完成数据采集与处理，通过可视化软件实时显示重力异常、位置、速度、姿态等信息。图 重力仪 该次重力实验测量任务开始前，测量人员在珠江口外海进行了仪器精度测试，完成了重复线和交叉点测量，在确保仪器精度满足要求之后，测量船再去测区开展作业，测线图如图所示。图测线图 测线包括东西测线、南北测线、和 斜线。测线间距。其中、测线各有两条重复线。数据处理方法及评价指标 小波变换理论小波变换的出现弥补了傅里叶变换中的不足。小波是一种长度有限、平均值为的波形，主要特点有：在时域具有紧支性或近似紧支性；直流分量为。小波变换的数学公式可以表示如下：系统工程与电子技术第 卷 （，）（），（）（）（）式中：为尺度因子；反映位移。在不同尺度下小波的持续时间随式（）右边值的加大而增宽，幅值则与式（）右边值反比减少，但波的形状保持不变。常用的小波基函数有 小波、多贝西（，）小波、小波、墨西哥小波（，）和 双 正 交 样 条 小 波（，）等。本文选择了 小波和 小波两种基函数用于数据的变换处理，并对变换结果进行对比分析。误差比率误差比率是指先将相邻重力异常值作差，然后再将两个相邻的差值相比得到的数值，具体计算如下：（）误差比率可以显示重力异常差值的变化趋势及规律。本文选取误差比率作为指标旨在验证信噪比较低的数据经过小波降噪后是否可以保留原始数据的变化规律及趋势。数据处理 漂移检查重力仪固有的一大缺点就是存在零点漂移问题。但是，只要漂移变化的幅度不大，并且存在一定的规律性，漂移就可以通过一定的方法被检测出来并加以改正。通过对数据处理分析得出该航次重力仪漂移量，如表所示，仪器的日漂移为 ，月漂移为 ，该船次重力测量总漂移量为 。表漂移统计 前校时间后校时间漂移 总漂移日漂移月漂移 ：测线内符合精度对重复测线计算其内符合精度，计算公式如下：（）（）式中：为重复测线测量值均方根误差；为重复测线的第一次测量的第个点的值；为重复测线的第二次测量第个 点 的 值。具 体 各 测 线 的 内 符 合 精 度 如 表所示。表各测线内符合精度表 测线均方根误差 续表 测线均方根误差 通过各测线的标准差统计结果可以看出 重力仪的测量精度在海况良好的情况下可以满足科考和其他海洋勘探任务的测量精度要求。数据降噪本文进行数据降噪处理是对所有测线总体数据进行整体处理，小波变换方法选取了 和 两种小波基函数。小波选择的是 小波基，分解层数为层；小波分解层数同样为层。降噪处理结果如图图 以及表、表所示。图中，数据中添加的倍随机噪声服从均匀分布。图中添加的正态分布随机噪声属于一般正态分布。图中，为倍随机噪声与原始数据之差；为经 小波变换后的倍随机噪声与原始数据之差；为正态分布随机噪声与原始数据之差；为经 小波变换后的正态分布随机噪声与原始数据之差。图中，为倍随机噪声与原始数据之差；为经 小波 变 换 后 的倍 随 机 噪 声 与 原 始 数 据 之 差；为正态分布随机噪声与原始数据之差；为经 小波变换后的正态分布随机噪声与原始数据之差。图原始数据与随机误差 图正态分布随机误差 第期刘备等：小波变换在船载海洋重力测量数据处理与分析中的应用 图 小波变换 图 小波变换 图 小波变换后误差对比图 图 小波变换后误差对比图 图 小波变换后误差比率图 图 小波变换后误差比率图 表标准差统计表 误差类型未变换 倍随机误差 正态随机误差 表相对精度统计表 误差类型 倍随机误差 正态随机误差 系统工程与电子技术第 卷表误差比率统计表 误差类型未变换 倍随机误差 正态随机误差 对比原始数据、噪声数据以及经过降噪处理的数据可以发现，加入倍随机噪声和正态分布随机噪声的数据经过 小波变换处理后数据绝对精度分别提高了 和 ；采用 小波变换后，分别提高了 和 。通过图、图以及表可知，经过小波变换处理后的数据相对精度至少提高，随着噪声的增大，最高可提高。对图、图 以及表中的误差比率变化情况进行分析，可以看出经过降噪后的重力异常数据能够很好地保持数据的基本变化趋势不变，并且经过小波变换后的数据视图更加简洁，能够较为轻松地判别测量过程中是否存在异常体导致的测量结果突变，为利用重力异常进行探测提供判别依据。结论本文研究了小波变换方法在船载重力测量数据降噪处理与分析中的应用，将 小波和 小波两种小波基函数引入到重力场数据降噪中，并利用 重力仪实测数据验证降噪方法的精度与可靠性，结果表明小波变换方法可以很好地提高数据精度，加入倍随机误差和正态分布随机误差的数据经过 小波变换处理后数据绝对精度分别提高了 和 ，相对精度分别提高了 和；采 用 小 波 变 换 后，分 别 提 高 了 和 ，相对精度分别提高了 和。对误差比率变化情况进行分析可以推断，经过小波变换方法降噪后的数据可以保留重力异常变化趋势，并且更加简洁，能够为重力探测判别是否存在异常体提供判别依据。参考文献黄谟涛，翟国君，管铮，等海洋重力场确定及其应用北京：测绘出版社，：，：，（）：陆飚，钟波，罗志才波罗的海地区高精度海洋重力测量数据处理与分析武汉大学学报（信息科学版），（）：，（）：，（）：陈上及海洋数据处理分析方法及其应用北京：海洋出版社，：，孙鹤洋，金涛勇，赵闯融合卫星测高和船测重力数据的非均匀密度点质量方法测绘地理信息，（）：，（）：刘志鑫，孟小红，王俊，等海域重力场解算中垂线偏差方法的并行化改进地球物理学进展，（）：，（）：吴燕雄，滕云田，吴琼，等船载绝对重力仪测量系统的误差修正模型及 不 确 定 度 分 析 武 汉 大 学 学 报（信 息 科 学 版），（）：，（）：，（）：宁津生，黄谟涛，欧阳永忠，等海空重力测量技术进展海洋测绘，（）：，（）：赵德军，孙中苗，赵东明基于因子克里金分析方法的航空重力梯度降噪中国惯性技术学报，（）：，（）：，（），（）：，：，（）：钟波基于 卫星重力测量技术确定地球重力场的研究武汉：武汉大学，：，第期刘备等：小波变换在船载海洋重力测量数据处理与分析中的应用 瞿庆亮，常晓涛，朱广彬，等基于地面重力的卫星重力梯度检校方法地球物理学报，（）：，（）：，（）：蔡劭琨，吴美平，张开东，等经验模态分解在动态重力测量数据处理中的应用海洋测绘，（）：，（）：董庆亮，陈洁，潘乐，等基于 重力数据检查海洋重力测量成果的质量海洋测绘，（）：，（）：黄谟涛，刘敏，孙岚，等海洋重力测量动态环境效应分析与补偿海洋测绘，（）：，（）：欧阳永忠海空重力测量数据处理关键技术研究武汉：武汉大学，：，蔡劭琨航空重力矢量测量及误差分离方法研究长沙：国防科技大学，：孙中苗，夏哲仁 低通差分器的设计及其在航空重力测量中的应用地球物理学报，（）：，（）：严飞，李达，李中，等一种用于重力梯度动态测量的载体环境引力梯度补偿方法中国惯性技术学报，（）：，（）：黄谟涛，陈欣，邓凯亮，等补偿海空重力测量动态效应剩余影响的通用模 型测绘学报，（）：，（）：孙鹤泉，金绍华，张宇基于 变换的海洋重力观测航行数据滤波方法海洋通报，（）：，（）：作者简介刘备（），博士 研究 生，主 要 研 究 方 向 海 洋 重 力 场 探 测 与导航。边少锋（），教授，博士，主要研究方向大地测量学。纪兵（），副教授，博士，主要研究方向海洋物理场导航与探测。贤鹏飞（），硕士研究生，主要研究方向重力梯度探测导航。