1、基于多偶极子磁性模型的基于多偶极子磁性模型的卫星复杂磁场建模方法卫星复杂磁场建模方法孟立飞1，肖琦2，王国强1，陈金刚2，易忠2，吴明雨1，张铁龙1,3（1.哈尔滨工业大学（深圳），深圳518055；2.北京卫星环境工程研究所，北京100094；3.奥地利科学院空间研究所，格拉茨A-8042）摘要：飞行器、舰船和潜艇等典型可移动平台具有复杂的磁场，获取这些磁场信息对于提高磁场探测精度具有重要意义。文章在简要总结典型可移动平台复杂磁场的建模方法基础上，结合卫星自身特点，提出基于多偶极子磁性模型的卫星复杂磁场快速建模方法。其技术途径有两种：一是基于单机磁试验数据建立单机磁性模型，再组合成整星多偶极
2、子磁性模型；二是基于卫星磁试验数据直接建立整星多偶极子磁性模型。采用这两种技术途径完成某型号卫星的磁场建模、计算与测量验证，测试结果表明多偶极子法可很好地适用于磁洁净卫星复杂磁场的快速建模计算，精度优于 0.2nT。关键词：卫星；磁试验；磁场建模；多偶极子法中图分类号：TB115;V474.1文献标志码：A文章编号：1673-1379(2023)04-0356-05DOI:10.12126/see.2023104A method for modeling complex magnetic field of satellite based onmulti-dipole magnetic mode
3、lMENGLifei1,XIAOQi2,WANGGuoqiang1,CHENJingang2,YIZhong2,WUMingyu1,ZHANGTielong1,3(1.HarbinInstituteofTechnology(Shenzhen),Shenzhen518055,China;2.BeijingInstituteofSpacecraftEnvironmentEngineering,Beijing100094,China;3.SpaceResearchInstitute,AustrianAcademyofSciences,Graz8042,Austria)Abstract:Typical
4、movableplatformssuchasaircrafts,ships,andsubmarineshavecomplexmagneticfields.Itissignificanttoobtainthemagneticfieldinformationforimprovingthemagneticfielddetectionaccuracy.Basedonabriefsummaryofmethodsformodelingcomplexmagneticfieldsoftheseplatforms,andtakingintoaccountthecharacteristicsofsatellite
5、s,amethodformodelingthecomplexmagneticfieldsofsatellitebasedonthemulti-dipolemagneticmodelwasproposedinthispaper.Therearetwotechnicalapproaches:Onewastocreateasingle-machinemagneticmodelbasedonthesingle-machinemagnetictestdataandthencombineallthemodelsintoawhole-satellitemulti-dipolemagneticmodel;th
6、eotherwastodirectlyconstructthemulti-dipolemagnetic model based on the test data of satellite magnetic.The magnetic field modeling,calculation andmeasurementverificationofasatellitehavebeencompletedalongthetwotechnicalpathsabove.Theresultsshow that the multi-dipole method is suitable for the rapid m
7、agnetic modeling and calculation of complexmagneticfieldsofthemagneticcleansatellitewithanaccuracyofbetterthan0.2nT.Keywords:satellite;magnetictest;magneticfieldmodeling;multi-dipolemethod收稿日期：2023-04-03；修回日期：2023-07-31基金项目：国家重点研发计划（编号：2022YFA1604600）引用格式：孟立飞,肖琦,王国强,等.基于多偶极子磁性模型的卫星复杂磁场建模方法J.航天器环境工程,
8、2023,40(4):356-360MENG L F,XIAO Q,WANG G Q,et al.A method for modeling complex magnetic field of satellite based on multi-dipole magneticmodelJ.Spacecraft Environment Engineering,2023,40(4):356-360Vol.40,No.4航天器环境工程第40卷第4期356SPACECRAFTENVIRONMENTENGINEERING2023年8月https:/E-mail:Tel:(010)68116407,6811
9、6408,68116544 0 引言引言为提高卫星、飞机等飞行器平台的磁场探测精度，需对源自平台的磁场干扰加以控制。进行控制的前提是通过磁场建模获取平台的磁场信息，而仿真分析是磁场建模的常用手段，更是复杂磁场建模方法的关键。自 20 世纪中期开始，国内外已对飞机、舰船及潜艇等典型平台的磁场开展了大量的理论及实验测试工作，并发展了多种磁场建模方法1-4。然而，关于卫星平台磁场的研究起步较晚，相关建模方法仍有待进一步完善。本文首先简要介绍典型平台的磁场建模方法，分析这些系统的复杂磁场建模方法在卫星上的适用性；然后结合卫星自身特性，研究卫星复杂磁场的建模方法；最后用真实卫星的磁场参数验证所提出的复杂
10、磁场建模方法的精度。1 复杂磁场建模方法复杂磁场建模方法1.1飞机干扰磁场建模方法在航空磁测中，由于磁力仪安装在飞机内部，飞行平台产生的磁场会对目标磁测产生很大干扰，所以要精确测量环境磁场，就需要将飞机的干扰磁场剔除。Tolles 等建立了干扰磁场的数学模型，该模型给出了由飞机机动产生的固定磁场、感应磁场和涡流磁场的数学表达式，称为 Tolles-Lawson方程5；后续大量的研究都是在此理论基础上开展的。飞机干扰磁场的矢量表达式为B=Bre+Bin+Bvo=(axex+ayey+azez)+(bxxBx+bxyBy+bxzBz)ex+(byxBx+byyBy+byzBz)ey+(bzxBx+
11、bzyBy+bzzBz)ez+(cxxdBxdt+cxydBxdt+cxzdBxdt)ex+(cyxdBydt+cyydBydt+cyzdBydt)ey+(czxdBzdt+czydBzdt+czzdBzdt)ez，(1)式中：ex、ey和 ez为方向矢量；Bx、By和 Bz为磁场测量分量；其余参数为各干扰系数。飞机干扰磁场的标量表达式为B=3i=1aiUi+BE3i=13j=1bijUiUj+BE3i=13j=1cijUiUj，(2)Uj式中：Ui和 Uj分别为地磁场和干扰磁场在飞机坐标系中三个坐标轴的方向余弦；为 Uj对时间的导数；ai、bij和 cij为干扰系数；BE为地磁场。通过设计多
12、次复杂的飞行轨迹，记录飞机在不同飞行状态、姿态和位置下磁力仪接收的磁场信息，便可计算出数学模型的干扰系数，从而提高后续飞行的磁场测量精度。该方法的优点是：无需面向飞机建立复杂的磁性模型，而是设计飞行轨迹，通过实际测量的磁场数据计算出干扰系数，从而消除飞机平台对环境磁场测量的影响，最终达到高精度磁场探测的目的。1.2舰船、潜艇磁场建模方法舰船、潜艇的壳体由钢铁材料加工而成，对内部设备有较大的磁屏蔽作用，因此舰船、潜艇的外部磁场主要是壳体的剩磁场和感应磁场，其磁场建模主要面向金属壳体。壳体具有界面连续、材料布局相对均匀、壁薄和磁导率较大等特征，适宜应用数值法进行建模计算。对于潜艇，甚至可将其定义为
13、标准的旋转椭球体，应用解析法进行建模计算，文献 4 就采用了该方法。磁偶极子阵列模型是另一种应用较多的潜艇磁场建模方法6。该模型的基本原理是在壳体上均匀设置 n 个磁偶极子，位置 ri已知，磁矩 Mi未知，下标 i=1,2,n。在 m 个测量点测量潜艇产生的磁场值 Bj，则Bzj=04ni=1(GxjiMxi+GyjiMyi+GzjiMzi)，j=1,2,m。(3)式中：G 矩阵的各个分量可根据偶极子磁场的计算公式推导出来，是偶极子位置 ri和测量点位置 rj的函数。由于磁场测量点的位置均已知，G 矩阵各个分量值均为常量。根据式(3)可知，只需要测量m 个位置的磁场 z 分量，就可以得到 m
14、个线性方程组，从而计算出所有 Mi值，最终得到需要的磁偶极子阵列模型(M,r)。1.3卫星磁场建模的特点相对于飞机而言，由于卫星在既定轨道上运第4期孟立飞等：基于多偶极子磁性模型的卫星复杂磁场建模方法357行，很难实现在空间中自由飞行或随意变换姿态，所以飞机干扰磁场的建模方法不适用于卫星。基于轻量化的要求，卫星的舱体一般使用铝蜂窝或者碳结构等无磁材料制造，所以卫星的磁性主要来源于其内部近百台单机中大量的电子元器件以及离散的软硬磁材料。加之卫星磁源比舰船、潜艇平台的更为复杂，因此将磁偶极子阵列模型应用于卫星复杂磁场的建模需要结合卫星自身特点进行优化。例如，卫星平台上磁源的位置 ri和磁矩 Mi均
15、未知，式(3)由线性方程变化为非线性方程，因此需要测量三分量的圆周磁场数据，从而获得足够的外部磁场信息，通过反演算法得到符合精度要求的多偶极子磁场模型(M,r)，再依据多偶极子磁场模型计算卫星任意 r位置的磁场值：B(r)=04ni=03(rri)(rri)Mi|rri|5Mi|rri|3。(4)2 卫星多偶极子磁场建模方法卫星多偶极子磁场建模方法卫星多偶极子磁场模型有两种建模方式：单机组合法和整星测量数据反演法。2.1单机组合法该方法通过测量卫星单机磁场分布，计算出单机磁偶极矩，获得所有单机的磁偶极矩及其在卫星坐标系下的方位，进而组合得到整星的多偶极子磁场模型7，用Mi(i=1,2,N)和r
16、i(i=1,2,N)表示，其中 N 为单机总数。单机磁偶极矩 Mi的计算可参照航天器磁矩测量的标准方法，如偶极子法、球面作图法或赤道作图法8。单机组合法的优点为：单机的磁场测量非常容易；依据所有单机磁场模型组合成卫星磁场模型，无须整星测量，没有尺寸限制；磁偶极子和单机一一对应，与卫星磁源实际分布符合度高。单机组合法的缺点为：由于卫星复杂的构型、内部电缆网的影响、互感的影响以及单机长期的磁性状态变化等，造成磁偶极子磁场模型与卫星真实结果有一定偏差，尤其当磁洁净卫星伸杆远端磁场很小时会有更大的相对偏差；所有单机需要磁试验。单机组合法常用于卫星磁性预估、设计优化和过程磁性控制，但一般不能作为最终的卫
17、星磁洁净指标验收手段。2.2整星测量数据反演法首先测量卫星圆周磁场分布 Bm（m 个测量点），设置一组多偶极子初始磁场模型，代入式(4)，计算出 m 个测量点位置的磁场计算值 B。定义磁场计算值和测量值的方差函数为C(M,r)=eTe=(BmB)T(BmB)。(5)C(M,r)其中磁偶极子磁场模型参量 M 和 r 的优化值由函数的最小化原理确定，MC=0rC=0。(6)通过对多偶极子磁矩和位置的双重迭代优化，最终得到符合精度要求的磁偶极子磁场模型(M,r)，具体理论国内外都有文献报道9-11。整星测量数据反演法的优点为：1）根据卫星整星测量数据建模计算，可信度较高；2）通过测量卫星近处的磁场分
18、布能建模计算远处的微弱磁场值。整星测量数据反演法的缺点为：1）反演难度大，多偶极子磁矩和位置参数自由多重迭代，存在非线性运算，可能由于测量数据误差等因素造成难收敛甚至不收敛等问题；2）结果不唯一，最终的多偶极子模型与实际磁源不是一一对应的。整星测量数据反演法可用于卫星磁试验数据分析和磁洁净指标验证。2.3建模方法的卫星应用与结果校验以某磁洁净卫星为例（如图 1 所示）：星载磁力仪探头安装在长伸杆的顶端，在该位置上磁场洁净度控制指标非常严格，要求不超过 0.5nT。在卫星研制阶段，对该卫星实施了严格的磁洁净控制措施，并应用多偶极子法进行了磁场建模计算。首先对卫星所有单机进行了磁试验，计算出每台单
19、机的磁偶极矩。然后按照单机组合法，在卫星坐标系下组合成多偶极子磁场模型，再依据磁偶极子模型计算出卫星本体在磁力仪探头位置上产生的磁场值 B1。磁力仪探头安装位置磁力仪探头安装位置图1磁洁净卫星构型Fig.1Configurationofamagneticcleansatellite358航天器环境工程第40卷此外，对卫星进行了整星磁试验，磁场测量探测器布局如图 2 所示，其中 8 号探头布置在星载磁力仪探头安装位置。卫星被放置在可旋转 360的水平无磁转台上，转台每转动 10，测量一组磁场数据。123567849西南r1r2=240 cmr3w=67 cmh1=160 cmh2=230 cmh
20、=70 cmh=70 cmL1L2L3=506 cm图2磁场测量探测器布置示意Fig.2Layoutdiagramofmagneticmeasurementdetectors利用 2 号、4 号、5 号探头的近场圆周分布测量数据，根据反演法构建了卫星的多偶极子磁场模型，偶极子数目为 16 个。如图 3 所示，离散点为3 个探头 9 个磁场分量的实际测量值，曲线为多偶极子磁场模型的计算值，二者符合性很好。根据多磁偶极子磁场模型计算得到卫星本体在磁力仪探头位置上产生的磁场值 B2。角度/()3503002502001501005010.015.05.00.0-5.0-10.0-15.0磁场值/nT
21、 图3磁场测量值与多偶极子磁场模型的计算值Fig.3Measured values of magnetic field and calculatedvaluesofmulti-dipolemagneticmodel依据 8 号探头的实测数据得到卫星本体在该磁力仪探头位置上产生的磁场值 B3。对比了组合法和反演法计算结果，以及直接测量结果，如表 1所示。由于卫星本体磁场在磁力仪位置处的值非常小，所以利用磁力仪直接测量的数据存在 0.1nT 量级的误差。即便如此，从表 1 中仍可发现，反演法与直接测量的最大偏差不超过 0.13nT；组合法与直接测量的最大偏差为 0.32nT。该结果表明：反演法有更
22、好的精度；组合法计算结果虽然有一定的偏差，但是仍然能够有效指导卫星研制过程的磁洁净度控制，确保最终指标受控。表1磁偶极子模型的计算值与磁场实测值对比Table1Comparison between calculated values of multi-dipole magnetic model and measured values ofmagneticfield磁场值x分量/nTy分量/nTz分量/nTB1（组合法）-0.040.070.15B2（反演法）-0.140.130.47B3（直接测量）-0.010.050.47 3 结束语结束语我国未来规划的专业磁测卫星和深空磁场探测任务对卫星磁
23、洁净度将有更高的要求，这也促使卫星磁场建模计算技术不断发展。目前基于反演法的磁洁净卫星磁场建模计算的精度可优于 0.2nT，为进一步提高计算精度，将从以下三个方面开展多偶极子建模方法的深入研究：1）在单机组合法方面，开展单机磁场模型的细化研究，增加外磁场感应和单机互感模型，以进一步提高组合模型的计算精度。2）在整星测量数据反演法方面，继续开展非线性迭代算法的研究，优化测量数据的布点，寻求规避不收敛奇异点的方法，提高运算效率。3）结合组合法和反演法的优点，在单机组合模型的基础上，进一步利用卫星圆周近场测量数据，研究多偶极子磁场模型的优化方法，得到反演运算更快收敛、磁源分布更符合实际、计算精度更高
24、的卫星磁场优化模型。参考文献（References）LINCS,ZHOUJJ,YANGZY.Amethodtosolvetheaircraft magnetic field model basing on geomagneticenvironment simulationJ.Journal of Magnetism andMagneticMaterials,2015,384:314-3191NGUYENTS,GUICHONJM,CHADEBECO,etal.Ships magnetic anomaly computation with integralequationandfastmultip
25、olemethodJ.IEEETransactionsonMagnetics,2011,47(5):1414-14172吕俊伟,于振涛,樊利恒,等.舰船磁场模型适用范围研究J.海军航空工程学院学报,2012,27(5):563-566LV J W,YU Z T,FAN L H,et al.Research on the3第4期孟立飞等：基于多偶极子磁性模型的卫星复杂磁场建模方法359applicability of warship magnetic modelJ.Journal ofNaval Aeronautical and Astronautical University,2012,27(
26、5):563-566赵建扬,林春生,孙玉绘.直升机平台背景磁干扰建模与特性分析J.海军工程大学学报,2016,28(1):37-40ZHAO J Y,LIN C S,SUN Y H.Modeling andcharacterizationofhelicopterplatformbackgroundmagnetic interferenceJ.Journal of Naval University ofEngineering,2016,28(1):37-404TOLLESWE,MINEOLANY.Compensationofaircraftmagneticfields:US2692970R,19
27、54:78-865周峰,刘卫东,林轶群.舰船磁场磁偶极子阵列数学模型研究J.计算机测量与控制,2008,16(10):1499-1501ZHOU F,LIU W D,LIN Y Q.Research on themathematical model of ship magnetic dipole arrayJ.Computer Measurement&Control,2008,16(10):1499-15016TAKAAKIN,SATOSHIS,SHIGERUA.Aclosed-formformulaformagneticdipolelocalizationbymeasurement7of i
28、ts magnetic field and spatial gradientsJ.IEEETransactionsonMagnetics,2006,42(10):3291-3293齐燕文.空间磁环境模拟技术J.航天器环境工程,2005,22(1):19-23QI Y W.Space magnetic environment simulationtechnologyJ.SpacecraftEnvironmentEngineering,2005,22(1):19-238肖琦,耿晓磊,陈金刚,等.低轨磁测卫星干扰磁场标定方法研究J.地球物理学报,2018,61(8):3134-3138XIAO Q,
29、GENG X L,CHEN J G,et al.CalibrationmethodsoftheinterferencemagneticfieldforlowEarthorbit(LEO)magnetic satelliteJ.Chinese Journal ofGeophysics,2018,61(8):3134-31389HANSEN R O,LAURA S.Multiple-source EulerdeconvolutionJ,Geophysics,2002,67(2):525-53510BULENTO.Locationanddepthestimationofpoint-dipoleandlineofdipolesusinganalyticsignalsofthemagneticgradient tensor and magnitude of vector componentsJ,JournalofAppliedGeophysics,2010,70(1):27-3711（编辑：许京媛）一作简介：孟立飞，研究员，主要从事空间磁学及应用技术研究。360航天器环境工程第40卷