混合粒子群算法ADPSO对月球物理参数反演的改进

计算物理 ›› 2017, Vol. 34 ›› Issue (6): 740-746.

混合粒子群算法ADPSO对月球物理参数反演的改进

钟振^1,2

1. 贵州师范大学物理与电子科学学院, 贵阳 550001;
2. 中科院国家天文台-贵州师范大学天文研究与教育中心, 贵阳 550001

收稿日期:2016-09-23 修回日期:2016-12-18 出版日期:2017-11-25 发布日期:2017-11-25
作者简介:钟振(1982-),男,博士,副教授,从事月球重力场和月球内部结构研究,E-mail:zhenzhongmail@163.com
基金资助:
国家自然科学基金（41404021）、贵州省科学技术基金（黔科合J字[2014]2128）和贵州师范大学博士科研项目资助

Lunar Geophysical Parameter Inversion with Admixed Particle Swarm Optimization ADPSO

ZHONG Zhen^1,2

1. School of Physics and Electronic Science, Guizhou Normal University, Guiyang 550001, China;
2. NAOC-GZNU Center for Astronomy Research and Education, Guizhou Normal University, Guiyang 550001, China

Received:2016-09-23 Revised:2016-12-18 Online:2017-11-25 Published:2017-11-25

摘要/Abstract

摘要： 常规粒子群算法PSO对低敏感参数较难进行最优估计，本文提出联合自适应惯性权重和变异的混合粒子群算法ADPSO.结果表明：最佳拟合效果不随粒子数和变异概率的增大而加强.就四参数估计而言，粒子数取400个，变异概率小于0.03，可实现最优拟合.采用最佳粒子数和较小的变异概率，利用新近月球重力场模型GL0990d和LRO激光测高数据，对月球南半球部分高地的物理参数进行最优估计.发现研究区域的模型导纳谱很好地拟合了观测值，模型重力异常与观测重力异常的残差较小，符合参数估计的要求，验证了算法的有效性.ADPSO适合于月球物理参数估计，为大规模的反演运算提供了参考.

关键词: 月球物理, 参数估计, 混合粒子群, 导纳谱

Abstract: Ordinary particle swarm optimization PSO fails frequently in estimation of low sensitivity selenophysics parameters. With an adaptive inertia weight and a mutation factor, an admixed particle swarm optimization ADPSO is proposed. It is found that misfit is not reduced as increasing number of particles and probability of mutation. For four-parameter inversion, the best-fitting parameters can be estimated as considering about 400 particles and a low probability (<0.03) of mutation. With employment of gravity filed model GL0990d and altimeter data from LRO, we apply the method into selenophysical parameter inversion on southern highland of the moon. It shows a best-fit between modeled admittance spectral and observed values. Small residual of gravity anomaly verifies success of parameter inversion and validity of the method. The approach can be used in selenophysics research and could provide a reference for large-scale estimation of parameters.

Key words: selenophysics, parameter estimation, admixed particle swarm, admittance spectral

中图分类号:

P223

钟振. 混合粒子群算法ADPSO对月球物理参数反演的改进[J]. 计算物理, 2017, 34(6): 740-746.

ZHONG Zhen. Lunar Geophysical Parameter Inversion with Admixed Particle Swarm Optimization ADPSO[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2017, 34(6): 740-746.

参考文献

[1] CROSBY A, MCKENZIE D. Measurements of the elastic thickness under ancient lunar terrain[J]. Icarus, 2005, 173(1):100-107.
[2] 李斐, 柯宝贵, 王文睿, 等. 利用重力地形导纳估计月壳厚度[J]. 地球物理学报, 2009, 52(8):2001-2007.
[3] HUANG Q, WIECZOREK M A. Density and porosity of the lunar crust from gravity and topography[J]. J Geophy Res, 2012, 117(E5):1-9.
[4] ZUBER M T, SMITH D E, WATKINS M M, et al. Gravity field of the moon from the gravity recovery and interior laboratory (GRAIL) mission[J]. Science, 2013, 339(6120):668-671.
[5] LEMOINE F G, GOOSSENS S, SABAKA T J, et al. High-degree gravity models from GRAIL primary mission data[J]. J Geophy Res, 2013, 118(8):1676-1698.
[6] KONOPLIV A S, PARK R S, YUAN D N, et al. The JPL lunar gravity field to spherical harmonic degree 660 from the GRAIL primary mission[J]. J Geophy Res, 2013, 118(7):1415-1434.
[7] KONOPLIV A S, PARK R S, YUAN D N, et al. High resolution lunar gravity fields from the GRAIL primary and extended mission[J]. Geophy Res Lett, 2014, 41(5):1452-1458.
[8] SMITH D E, ZUBER M T, NEUMANN G A, et al. Initial observations from the lunar orbiter laser altimeter (LOLA)[J]. Geophys Res Lett, 2010, 37(18):109-118.
[9] ZHONG Z, LI F, YAN J, et al. 2014. Lunar geophysical parameters inversion based on gravity/topography admittance and particle swarm optimization[J]. Adv in Space Res, 2014, 54(4):770-779.
[10] 钟振, 晏鹏, 叶茂, 等. 基于粒子群PSO算法的月球物理参数反演[J]. 计算物理, 2013,30(5):706-713.
[11] EBERHART R C, KENNEDY J. A new optimizer using particles swarm theory[C]//Proc 6th International Symposium on Micro Machine and Human Science. Piscataway, NJ:IEEE Press, 1995, 39-43.
[12] YANG X M, YUAN J, YUAN J, et al. A modified particle swarm optimizer with dynamic adaptation[J]. Applied Mathematics and Computation, 2007, 189(2):1205-1213.
[13] 温正. 精通MATLAB智能算法[M]. 北京:清华大学出版社, 2015.
[14] XIAO Y, CUI G, PENG F, et al. An improved particle swarm optimization for precocious phenomenon in nonlinear system engineering[J]. Chinese J Comput Phys, 2015, 32(6):693-700.
[15] WIECZOREK M A, NEUMANN G A, NIMMO F, et al. The crust of the moon as seen by GRAIL[J]. Science, 2013, 339(6120):671-675.
[16] 熊德永,钟振,刘高福. 高分辨率卫星重力和激光测高数据的月球岩石圈有效弹性厚度估计的分析及应用[J]. 地球物理学进展,2016,32(2):0622-0628.