基于贝叶斯证据的全球地磁场模型主磁场阶数分析.pdf
《基于贝叶斯证据的全球地磁场模型主磁场阶数分析.pdf》由会员分享,可在线阅读,更多相关《基于贝叶斯证据的全球地磁场模型主磁场阶数分析.pdf(9页珍藏版)》请在咨信网上搜索。
1、Chinese).Chinese Journal ofpaceScience,22023,43(4):600-608.D0I:10.11.04.2022-0009MA Sen,MA Jiahui,TONG Jizhou,LI Yunlong.Analysis of Global Geomagnetic Main Field Model Order Based on BayesianEvidence(in0254-6124/20243(4)-0600-09Chin.J.Space2Sci.研究论文|ResearchArticles空间科学学报基于贝叶斯证据的全球地磁场模型主磁场阶数分析马森1,2
2、马嘉卉1,21,31,3佟继周李云龙1(中国科学院国家空间科学中心北京100190)2(中国科学院大学北京100049)3(国家空间科学数据中心北京101407)摘要全球主磁场模型描述了主磁场的时空分布特征。模型的主磁场阶数是构建主磁场模型的关键问题之一。使用贝叶斯推理分析全球主磁场模型,依据贝叶斯证据比较模型阶数,为主磁场的阶数分析提供一种统计依据。利用Swarm卫星磁测数据,估计不同阶数的主磁场模型的证据。结果表明,在1 2 0 的模型阶数中,阶数N=12具有全局最佳证据。参照Jefrreys scale的阅值区间,数据对阶数N=12的偏好显著优于其他阶数。实验表明,主磁场阶数的证据推理可
3、用于研究主磁场贡献,结果与1 4阶球谐函数的功率谱分析相匹配。关键词同地磁场,主磁场模型,贝叶斯证据,模型比较,地磁功率谱中图分类号P353.1Analysis of Global Geomagnetic Main Field ModelOrder Based on Bayesian EvidenceMA Sen1,2MA Jiahuil.2TONG Jizhou1,3LI Yunlong1,31(National Space Science Center,Chinese Academy of Sciences,Beijing 100190)2(University of Chinese Ac
4、ademy of Sciences,Beijing 100049)3(National Space Science Data Center,Beijing 101407)AbstractThe global main magnetic field model describes the space-time distribution characteristicsof the main magnetic field.The order of the main magnetic field in the model is one of the key issues tobuild the mai
5、n magnetic field model.This paper used Bayesian inference to analyze the global geomag-netic main field model and compares model orders based on Bayesian evidence.It provided a statisticalperspective for the main field order selection.Using magnetic observations from Swarm satellites,the ev-2022-04-
6、20收到原稿,2 0 2 2-0 7-2 4收到修定稿E-mail:The Author(s)2023.This is an open access article under the CC-BY 4.0 License(https:/creativecommons.org/licenses/by/4.0/)马601森等:基于贝叶斯证据的全球地磁场模型主磁场阶数分析idence for different orders of the main field model was estimated.The results show that order N=12has the global bes
7、t evidence in model from 1 to 20.Referring to the threshold interval given by Jefrreysscale,the data preference for order N=12 is significantly better than other orders.The experimentshows that the evidence reasoning of the main magnetic field order can be used to study the main mag-netic field cont
8、ribution,and the results match the power spectrum analysis of spherical harmonic of or-der N=14.Key wordsGeomagnetic field,Main field model,Bayesian evidence,Model comparison,Geomagnetic power spectrum0引言地磁场是地球重要的固有物理场。地磁场模型可用于描述已知观测、预估磁场强度,揭示地球磁场内部机制等研究。构建全球地磁场模型的高斯球谐分析法(Spherical Harmonic Analysis
9、,SH A)使用阶数累加的球谐函数来表示地球内源场的两种成分主磁场和地壳场。其中低阶项主磁场成分占优,高阶项地壳场成分占优。低阶项中主磁场对应的最大阶数称为主磁场阶数。选择合理的主磁场阶数可以更准确地描述主磁场的贡献,是构建地球主磁场模型的核心点之一。主磁场阶数选择的理论依据是地磁场的空间功率谱。功率谱的拐点是不同物理过程的指示器,意味着前后两段功率谱对应着不同的物理起源。这是一种行之有效的选择模型中主磁场阶数的方法,被广泛应用于地磁场模型的构建。随着卫星磁测的观测精度与时空均匀覆盖能力的不断提高以及模型反演能力的不断改进,国际上相关的研究机构推出了一系列全球地磁场模型,表1 给出了部分基于S
10、HA方法的全球地磁场模型,并列出了其描述的磁场类型及主磁场对应的阶数。全球地磁场模型的构建需要综合考虑磁测数据精度、预估误差和主磁场长期变化等因素,模型中主磁场阶数的使用情况存在差异,可以考虑引人统计方法对主磁场阶数进行分析。本文将主磁场阶数不同的全球主磁场模型视为相互竞争的候选模型。选择合理的主磁场阶数,即从候选中选择最能真实反应观测数据的模型。由此主磁场阶数的选择转换为一个模型比较问题,可利用模型比较的方法进行统计分析和评估。目前,常用的模型比较准则是赤池信息准则(A k a i k e In f o r m a t i o n Cr i t e r i o n,A IC)和贝叶斯信息准则
11、(Bayesian Information Criterion,BIC),两者都是贝叶斯推理中模型贝叶斯证据的近似值7.8。贝叶斯推理是模型比较的前沿方法,已经在物理学问题,特别在天体物理学和宇宙学中得到了广泛的应用。Beltran等9 利用贝叶斯推理比较了包含等曲率扰动模态与纯绝热条件(不包括曲率扰动)两种情况下的宇宙学模型证据,并发现贝叶斯模型的选择倾向于纯绝热模型。Trottal0 介绍了Savage-Dickey密度比来确定两个嵌套模型的贝叶斯证据,进而依据其进行模型选择。Arregui 等分别基于密度尺度高度三个参数和磁通管理论反演出两个日冕环内多模谐扭结振荡模型,然后使用贝叶斯推理
12、来比较两个模型的可信度。贝叶斯推理也可用于分析天文学中的模型参量和心理学的假设检验可靠性1 2-1 5。当当前,贝叶斯推表1基于SHA方法的全球地磁场模型及其主磁场对应阶数Table 1Geomagnetic field model and itscorresponding main field order basedon SHA method模型名称磁场类型主磁场阶数发布时间IGRF-13(2MM:1132020CGGMI13MM:1152021COV-OBS.x214MM:1142020CHAOS-715M+L+EM:1202020EMM2017M+LM:1122017POMME-10M+
13、L+EM:1152016CM510M+LM:1202015WMM2020M+LM:1202015注M代表主磁场,L代表地壳场,E代表外源场。6022023,43(4)Chin.J.SpaceSci.空间科学学报理的方法在地磁场模型研究中的应用较为少见,Swarm卫星的高精度地磁观测量给数据分析和处理带来了机遇和挑战。与基于地磁空间功率谱的传统阶数选择方法相比,贝叶斯推理只依赖于观测数据和模型结构,能够客观地比较数据对不同阶数全球主磁场模型的偏好程度,给主模型阶数选择提供一种统计依据。1贝叶斯推理贝叶斯推理可以量化模型与参数的不确定性,度量数据对模型的偏好程度,其原理源自贝叶斯定理的深人 1 6
14、,1 7 O假定提出一种合理的模型,即可通过拟合模型参数来描述已知的观测数据,模型参数的合理性由参数后验概率给出。这个概率通过贝叶斯公式展开可得P(D|w,H,)P(w|H,)P(w|D,H,)=(1)P(D|H,)其中,H;代表一个一般化的模型(H表示一个模型的理论推测或假设),w为待拟合的参数向量,D表示数据集。式(1)左侧是模型H;参数向量的后验概率。给定w和H,等式右侧P(Dlw,H,)为似然函数,代表模型产生数据的概率。模型参数的先验P(w|H,)定义了模型参数较为合理的取值范围。式(1)中分母P(D|H,)为归一化常数,称为边缘似然或贝叶斯证据。式(1)给出了参数后验概率分布与似然
15、函数、参数先验和贝叶斯证据之间的计算关系。在模型拟合的过程中,证据P(D|H)常常被忽略,因为其与参数w的选择无关。给定数据集D,数据对模型的偏好程度为P(H;ID),又称模型后验,可表示为P(D|H,)P(H,)P(H,ID)=(2)P(D)归一项P(D)与模型无关,即P(H;|D)P(D|H,)P(H,).(3)因此,P(H;|D)正比于证据P(D|H,)与模型先验P(H,)的乘积。P(H)为模型假设是正确的先验概率,其表示与数据无关的模型偏好,即推导模型所用理论的合理性。这里假设不认为存在一个或某些特定的方法适合于某一类特定问题,即模型与数据的信息结构无关,那么证据就是决定因素,P(DH
16、,)可以等效为数据对模型的偏好。证据可通过对式(1)进行积分得到,即Pevidence=P(D|H,)=P(D|w,H,)P(w|H,)dw.(4)式(4)证据的积分需构建模型的似然函数与参数先验。在实际应用中,证据积分的解析值不可得,可采用基于蒙特卡洛的采样算法来进行近似估计,如模拟退火、嵌套采样等。需要注意的是,证据采样算法的计算代价通常是相当昂贵的,模型复杂度、数据相关性和数据集大小都将显著地影响证据估计的时间与收敛情况。2数据集与地磁球谐模型2.1磁测数据集磁测数据集为构建CHAOS系列模型所使用的地磁数据集(以下称为CHAOS数据集)。CHAOS模型数据集中近2 0 年的卫星磁测数据
17、来源于Orsted,CHAMP,SAC-C,Cryo-Sat-2和 Swarm。具体地,Swarm数据使用MAGX_LR_1B1Hz校准数据,基线0 50 5/0 50 6,并对 Swarm的Alpha,Bravo,Char-lie三颗卫星进行了1 min的初采样。为避免地磁与太阳活动对观测数据的干扰,CHAOS数据集应用了理论完备的筛选准则来获得高质量的磁静期数据。数据集还包括了由CHAOS-6x9模型模拟的非主磁场贡献修正量,以便于主磁场模型的研究。2.2全球地磁场球谐模型在电磁场的准静态近似下,矢量磁场B可以表示为标量位势的负梯度B=VV=Vint+Vext。对于地磁内源场,在地心地固(
18、ECEF)球坐标系中,采用球谐展开可得到地磁标量位势的解1 8,即Nn+1gm(t)cosm入+n=1m=0hm(t)sinm|Pm(cos0).(5)其中,代表地球表面的参考半径,一般取=6371.2km;(r,0,入)表示球坐标下的地心距、余纬度和经度;Pm(c o s 0)为Schmidt半规范化关联Le-gendre函数;gm,hm为待拟合的球谐系数;N为球谐展开的最大阶数;m和n为球谐函数的阶和次。式(5)可单独或综合表示主磁场与地壳场。在构603马森等:基于贝叶斯证据的全球地磁场模型主磁场阶数分析建主磁场模型时,N即为主磁场阶数。CHAOS数据集中的模拟值可以用来修正观测数据中的地
19、壳场与外源场Vext成分。地磁场的矢量强度(X,Y,Z)可以表示为标量位势的负梯度,地磁标量强度F的参数模型可由矢量强度三分量得到F=VX?+Y2+Z2=oV22V6rsinga入r3主磁场模型证据估计式(4)表明,主磁场球谐模型的证据估计是一个高维的积分运算。首先需要推导似然函数P(Dw,H,)的表达式,然后确定参数先验P(w|H,)。结合二者可以得到证据积分的解析式,最后使用采样算法近似估计积分的真实值。3.1地磁标量强度的似然函数本文主要使用CHAOS数据集中的地磁场标量强度F数据及其观测噪声误差。为推导模型的似然函数,需构建地磁场标量强度F的噪声模型,有Fots=Fmodel+Eero
20、r(7)其中,Fobs表示实际观测量;Fmodel表示模型预测值,其参数表达可由式(5)和式(6)得到;Eeror为模型的预测误差,理论的预测误差源于观测时的噪声。在仅关心主磁场部分时,非主磁场成分同样被视为噪声error由3部分噪声组成,即Eerror=eobs+Ccrust+Cext.(8)式中,eobs为真实观测噪声,其分布为均值为零且独立分度的高斯噪声;erust为地壳场成分的修正;eext为外源场成分的修正。参照修正值的分布情况(见图1)可将其建模为高斯分布,且相互独立(见图2)。因此有Cos N(0,Oobs),Cestd N(uerust,Ostda),Cext N(uext,C
21、ext).(9)联合分布为Eeror N(uerust+Uext,obs?+Cerust?2+Oext2(10)因此,似然统计量x可表示为NobsFFmodel(uerust+uext)12x=Fobs,i(11)22=1式中,i表示第i个观测似然函数一般采用对数形式,即-lgLr(Fobsl(r,0,),A,w)=Nobs1NXD12(2元)+g121Nobsgg;+Nobslgoerust+Nobslg ext:(12)工其中,Nobs为数据量的总数。3.2球谐系数的先验当地磁场模型的球谐函数形式确定后,应给定球谐系数的先验分布,系数先验会影响模型对可观测量的预测。最常用的参数先验是零信息
22、先验(Non-informativePrior),先验在一个区间上为常数,故又80N70N60N50550N40N50030N49520N10N490048510S20S48030SS40S47550S47060S70S80S465图1观测数据的空间与高度分布Fig.1Spatial and height distribution of observation data6042023,43(4)Chin.J.SpaceSci.空间科学学报0.0001645(a)(b)0.0001440350.00012300.00010250.00008200.00006150.00004100.000025
23、002500035000450005500000.050.100.15F/nTNoisea/nT(c)0.16(d)0.300.140.250.120.200.100.150.080.060.100.040.050.0200-6-4-202460510152025Fstd/nTFext/nT图2数据直方分布与密度曲线。(a)地磁标量强度F,(b)观测噪声方差,(c)地壳场贡献模拟修正量,(d)外源场贡献模拟修正量Fig.2Data histogram distribution and density curve.(a)Geomagnetic scalar intensity F,(b)obse
24、rved noisevariance,(c)simulation correction of crustal field,(d)simulation correction of external field称为平坦先验,表示对模型的系数没有预先的偏好。平坦先验等价于均匀分布,即P(alA,R II:(13)其对数形式为lgP(w|A,R)=C.(14)式(1 3)和式(1 4)中,w为参数向量,A代表模型结构(A g r i c u l t u r e),R 代表先验的结构形式,d;表示宽度,C为一常数。参数先验决定了证据积分运算的参数空间大小。参数均匀分布的中心点与宽度d,可参考CHAOS和
25、IGRF等系列模型对应阶次的球谐系数值,可有效减少参数空间过大或先验区间偏离最佳后验参数值太远而导致的不必要计算。限制参数gi的先验分布全部位于负区间,可以避免因F开方使系数出现双根。3.3模型证据地磁场球谐模型系数的后验分布由如下贝叶斯定理给出Lr P(wN,R)P(w|N,D,R)=(15)p(DIN,R)其中,球谐函数的阶数N代替了模型结构A,R 代表先验的形式,Lr为似然函数。代人似然与平坦先验,对数证据可由积分得到lgP(D|N)=J Cr Cdw.(16)忽略常数C,即在平坦先验下,似然与参数后验等价。可得如下对数证据的积分解析形式NobsNobsIgP(DIN)=JZNobs1g
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 贝叶斯 证据 全球 地磁场 模型 磁场 分析
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。