基于正演的拖曳式航磁全张量梯度数据软补偿方法
【技术领域】
[0001] 本发明属于航磁软补偿领域，尤其涉及一种基于正演的拖曳式航磁全张量梯度数 据软补偿方法。
【背景技术】
[0002] 航磁软补偿主要目的在于消除航磁测量中飞机由于自身磁场（恒定场、感应场和 涡流场）以及飞行姿态变化等对磁测值的影响。
[0003] 传统航磁全张量梯度测量方法主要采用机载的方式，将飞机自身磁干扰场以及飞 机姿态变化等对测量值影响用线性方程组表示，各项干扰的系数称为补偿系数，通过学习 飞行获得补偿系数进而消除其影响。该种测量方式下由于磁梯度测量装置与机身距离很 近，磁干扰场组成及其变化都十分复杂。随着航空全张量磁梯度测量装置的研宄发展，航空 磁测方式也发展到采用直升机拖拽式测量，该种测量方法优点是磁梯度测量装置离机身较 远，直升机磁干扰减少，但是由于在飞行过程中直升机与全张量磁梯度测量装置相对位置 都时刻变化，传统机载式磁补偿的方法并不适用。另外以往拖曳式测量中的直升机补偿参 数多通过在地面围绕机身周围所测磁场值求取，而实际测量过程中磁梯度测量装置位于机 身下方，这在补偿参数求取以及补偿方面都存在较大误差。
【发明内容】
[0004] 为了解决上述问题，本发明提供一种基于正演的拖曳式航磁全张量梯度数据软补 偿方法。针对航磁全张量梯度拖曳式测量方式设计了一种基于正演的拖曳式航磁全张量梯 度数据软补偿方法，相对于现有技术而言，实现了对航磁梯度测量过程中直升机和吊舱姿 态、相对位置以及直升机磁场均动态变化情况下的航磁全张量梯度数据软补偿，并且减少 相对误差。
[0005] 为了实现以上目的，本发明采用的技术方案是：
[0006] -种基于正演的拖曳式航磁全张量梯度数据软补偿方法，所述方法用于消除航 磁测量中飞机由于自身磁场以及飞行姿态变化对磁测值的影响，其特征在于，包括以下步 骤：
[0007] 步骤1)建立航磁全张量梯度测量正演补偿模型；
[0008] 步骤2)根据步骤1)中建立的模型设计学习飞行方案进行学习飞行；
[0009] 步骤3)根据上述步骤2)的学习飞行数据，确定补偿参数；
[0010] 步骤4)利用上述步骤3)中确定的补偿参数对航磁全张量梯度数据进行软补偿。
[0011] 优先地，步骤4)所述根据确定的补偿参数对航磁全张量梯度数据进行软补偿包 括一下步骤：
[0012] 步骤4. 1)根据确定的补偿参数，分析出在直升机姿态、自身磁场变化影响下当地 水平坐标系下吊舱位置三分量及梯度张量补偿值；
[0013] 步骤4. 2)根据吊舱提供的姿态信息，将航磁三分量及梯度张量实际测量值从吊 舱坐标系改正到当地水平坐标系下，得到修正值；
[0014] 步骤4. 3)根据上述步骤4. 2)得出的航磁三分量及梯度张量实际测量值从吊舱坐 标系改正到当地水平坐标系下的修正值减去上述步骤4. 1)分析得出的当地水平坐标系下 吊舱位置三分量及梯度张量补偿值，确定补偿结果。
[0015] 优先地，步骤2)所述设计学习飞行方案，该学习飞行方案进行三次学习飞行设 计：
[0016] 学习飞行一，选定平静磁场区氏一平稳飞行；
[0017] 学习飞行二，另选定平静磁场区H2二平稳飞行；
[0018] 学习飞行三，选定变化磁场区H3适当采取机动动作飞行。
[0019] 优先地，所述步骤3)补偿参数包括恒定磁矩、感磁参数和涡流参数。
[0020] 优先地，所述步骤根据补偿参数及直升机姿态信息分析出在直升机姿态、自身磁 场变化影响下第k个测点当地水平坐标系下吊舱位置梯度张量补偿值：
【主权项】
1. 一种基于正演的拖曳式航磁全张量梯度数据软补偿方法，所述方法用于消除航磁测 量中飞机由于自身磁场以及飞行姿态变化对磁测值的影响，其特征在于，所述方法包括以 下步骤： 步骤1)建立航磁全张量梯度测量正演补偿模型； 步骤2)根据步骤1)中建立的模型设计学习飞行方案进行学习飞行； 步骤3)根据上述步骤2)的学习飞行数据，确定补偿参数； 步骤4)利用上述步骤3)中确定的补偿参数对航磁全张量梯度数据进行软补偿。
2. 根据权利要求1所述的软补偿方法，其特征在于，步骤4)所述根据确定的补偿参数 对航磁全张量梯度数据进行软补偿包括一下步骤： 步骤4. 1)根据确定的补偿参数，分析出在直升机姿态、自身磁场变化影响下当地水平 坐标系下吊舱位置三分量及梯度张量补偿值； 步骤4. 2)根据吊舱提供的姿态信息，将航磁三分量及梯度张量实际测量值从吊舱坐 标系改正到当地水平坐标系下，得到修正值； 步骤4. 3)根据上述步骤4. 2)得出的航磁三分量及梯度张量实际测量值从吊舱坐标系 改正到当地水平坐标系下的修正值减去上述步骤4. 1)分析得出的当地水平坐标系下吊舱 位置三分量及梯度张量补偿值，确定补偿结果。
3. 根据权利要求2所述的软补偿方法，其特征在于，步骤2)所述设计学习飞行方案，该 学习飞行方案进行按以下三次顺序学习飞行设计： 学习飞行设计一，选定平静磁场区氏一平稳飞行； 学习飞行设计二，另选定平静磁场区H2二平稳飞行； 学习飞行设计三，选定变化磁场区H3适当采取机动动作飞行。
4. 根据权利要求3所述的软补偿方法，其特征在于，所述步骤3)补偿参数包括恒定磁 矩、感磁参数和涡流参数。
5. 根据权利要求4所述的软补偿方法，其特征在于，所述步骤4. 1)根据补偿参数及直 升机姿态信息分析出在直升机姿态、自身磁场变化影响下第k个测点当地水平坐标系下吊 舱位置梯度张量补偿值：
所述Bxx,Bxy,Bxz,Byx,Byy,Byz,Bzx,Bzy,Bzz为当地水平坐标系下第k个测点全张 量梯度，U〇为真空磁导率；X，y，z分别为第k个测点以直升机位置为原地当地水平坐标 系下吊舱位置；Hx，Hy，扎分别为第k个测点在直升机坐标系下x、y、z三方向上外磁场值； 分别为第k个测点在直升机坐标系下x、y、z三方向上外磁场变化率；r1 at dr at 为第k个测点对应的直升机姿态矩阵T逆矩阵，（r1)1为T逆矩阵的转置矩阵，a，e，Y 分别为该测点对应的直升机偏向角、俯仰角和横滚角；M"，M_M^eix，eiy，eiy，evx，evy，evz分 别为第k个测点直升机在直升机坐标系下x、y、z三方向上恒定磁矩、感磁参数、涡流参数； Mx,My, 1\^分别为第k个测点直升机总磁矩M在直升机坐标系下X、y、Z三方向的分量。
6.根据权利要求5所述的软补偿方法，其特征在于，所述步骤4. 2)将航磁三分量及梯 度张量实际测量值从吊舱坐标系改正到当地水平坐标系下，通过吊舱姿态变换矩阵T。进行 坐标系改成，所述T。由一时刻测量的吊舱偏航角a。、俯仰角和横滚角Y。组成：
所述忑1为矩阵T。逆矩阵，T。吊舱姿态变换矩阵；Hax。&、^ ^义^为吊舱位置经过吊 舱姿态改正后的三分量值 ，BxX-Car、BXy-Car、BXZ-Car、ByX-car、Byy-car、ByZ-car、Bzx -car、Bzy-car、吊 ?置吊idi、改正后的梯度张里值；Hax_measm:e、Hay_measm:e、Za_measm:e、Bxx_measm:e、BXy_measure、 Bxz-measure、ByX-measure、Byy-measure、ByZ-measure、Bzx-measure、BZy-measure、Bzz-measure 为吊 头际/则里值，根据 得出的上述值对各点进行分别补偿。
7.根据权利要求6所述的软补偿方法，其特征在于，所述补偿参数确定：
其中Mi、M2分别为学习飞行一、学习飞行二直升机坐标系下直升机总磁矩，113为学习飞 行三中直升机坐标系下的直升机涡流参数；hx1，hx2，hx3均为列向量，分别为学习飞行一、二 和三中各测点在当地水平坐标系下对应的地磁x分量，77^, 77^分别为Hxl，hx2的平均值； SHmi与SHax2分别为n行3列矩阵，n为测点总数；H^^与H"_。"2分别为学习飞行一和学习 飞行二的实测数据进行吊舱姿态改正后的值；SHax3Sn行3列矩阵，n为测点总数；Hax。"3为 实测数据进行吊舱姿态改正后的值，Hay为直升机恒定磁场及姿态变化对各测点三分量中 的X方向分量Hax造成的干扰值，Haxi为直升机感应磁场及姿态变化对各测点三分量中的X 方向分量Hax造成的干扰值，Hx3为各测点地磁场X方向分量，均为n行1列矩阵； 通过上述表达式可得出：
所述eix,eiy,eiy，evx,evy,evz，分别为直升机在直升机坐标系下x、y、z三方向 恒定磁矩、感磁参数、涡流参数;Mxl、Myl、MZJPMx2、My2、Mz2分别为直升机总磁矩Mi、M2在直升 机坐标系下X、y、z三方向上的分量；77^,77^',77^和T^',7^',77；'分别为稳定磁场一 和稳定磁场二中直升机坐标系下的地磁三分量在X,y,z方向上的平均值。
8.根据权利要求7所述的软补偿方法，其特征在于，所述步骤1)中补偿模型涉及四 种坐标系，分别是直升机坐标系、吊舱坐标系、直升机处为原点的当地水平坐标系，以及GPS 系统所采用的坐标系；所述四种坐标系均为三维坐标系。
【专利摘要】本发明属于航磁软补偿领域，尤其涉及一种基于正演的拖曳式航磁全张量梯度数据软补偿方法。针对航磁全张量梯度拖曳式测量方式设计了一种基于正演的拖曳式航磁全张量梯度数据软补偿方法，相对于现有技术而言，实现了对航磁梯度测量过程中直升机和吊舱姿态、相对位置以及直升机磁场均动态变化情况下的航磁全张量梯度数据软补偿，并且减少相对误差。
【IPC分类】G06F19-00
【公开号】CN104809352
【申请号】CN201510236807
【发明人】修春晓, 孟小红, 郭良辉, 陈召曦, 张盛, 王俊, 张兴东, 姚长利, 李淑玲, 郑元满, 刘国峰, 邹韵, 周俊杰, 马亚伟, 刘默
【申请人】中国地质大学（北京）
【公开日】2015年7月29日
【申请日】2015年5月11日