面向电离层反演的参数平滑迭代层析方法  

Iterative tomography method for ionosphere inversion via parameters smoothing

作　　者：余接情[1] 王颖 朱砚宇 黄英麒 吴立新[2] YU JieQing;WANG Ying;ZHU YanYu;HUANG YingQi;WU LiXin(School of Environment Science and Spatial Informatics,China University of Mining and Technology,Xuzhou Jiangsu221116,China;School of Geosciences and Info-physics,Central South University,Changsha410012,China)

机构地区：[1]中国矿业大学环境与测绘学院,江苏徐州221116 [2]中南大学地球科学与信息物理学院,长沙410012

出　　处：《地球物理学报》2023年第6期2263-2274,共12页Chinese Journal of Geophysics

基　　金：国家重点研发计划课题(2018YFB0505304);国家自然科学基金(41771416);国家重点研发课题(2018YFC1503505);中国博士后基金(2017M611949)共同资助。

摘　　要：层析是电离层三维结构反演的重要技术.增加约束可缓解层析模型的病态问题.然而,迭代层析中,大量无射线穿过的网格易对施加的约束产生与真实意图相违背的约束结果(称负面约束).针对该问题,提出了一种参数平滑的迭代层析方法(iterative tomography method via parameters smoothing, ITPS).每次迭代中,先用Chapman函数及最小二乘法拟合及改正每一垂直剖面,获取关于Chapman函数参数的二维图像,然后用移动窗口法平滑各参数图像并以此改正各网格的电子密度.实验表明:ITPS方法可一定程度上减少负面约束并抑制层像的扰动.相对于MART和CMART算法,ITPS方法在垂直剖面、F2层临界频率(f_(o)F_(2))、F_(2)层峰值高度(h_(m)F_(2))、斜电子总含量(Slant Total Electron Content, STEC)及h_(m)F_(2)以上电子密度等方面均具有更佳的表现.相对于精度较高的CMART算法,ITPS方法在f_(o)F_(2)与h_(m)F_(2)的平均优化率分别为7.49%及6.60%,STEC的平均优化率为5.19%,h_(m)F_(2)以上电子密度的平均优化率为11.41%.Tomography is an important technology for three-dimensional ionosphere inversion.The well-known problem of ill-posed in ionosphere tomography is usually addressed by introducing additional constraints into tomographic models.However,in the method of iterative tomography,the existence of the massive unpassed grids facilitates the method to produce an unexpected result for the added constraints(named negative constraint hereafter).This paper is to develop a new ionosphere tomography method,named iterative tomography method via parameters smoothing(ITPS),to lower down the negative constraint.In ITPS method,each vertical profile of the ionosphere is fitted and corrected by a Chapman function via a least square fitting method.After that,an image for each parameter of the Chapman is generated.The image is smoothed by a moving window method and then used to correct the electron density of each voxel.Real experiments showed that ITPS method reduces the negative effect of the constraints and the fluctuations of the tomographic image.When compared with the MART and CMART algorithms,ITPS shows a better performance on the vertical profiles,the peak frequency of F_(2)layer(f_(o)F_(2)),the peak height of F_(2)layer(h_(m)F_(2)),slant total electron content(STEC)and the electron density above h_(m)F_(2).The average improvement of the ITPS over the CMART algorithm,the better one,are 7.49%,6.60%,5.19%and 11.41%,respectively,on the aspects of f_(o)F_(2),h_(m)F_(2),STEC and the electron density above h_(m)F_(2).

关 键 词：电离层电子密度 电离层层析 迭代层析 Chapman函数 水平光滑 MART算法 

分 类 号：P353[天文地球—空间物理学]