反距离权重插值与克里金插值在冀鲁豫交界区地磁场数据分布中的对比研究
Comparison of Inverse Distance Weighting and Kriging Interpolation in the Distribution of Geomagnetic Field Data in Juncture Area of Hebei, Shandong and Henan
DOI: 10.12677/ag.2025.158110, PDF,    科研立项经费支持
作者: 李 博:河北省地震局保定地震监测中心站,河北 保定
关键词: 地磁场冀鲁豫地区反距离权重插值克里金插值Geomagnetic Field Area of Hebei Shandong and Henan Inverse Distance Weighting Kriging Interpolation
摘要: 为提升冀鲁豫交界地区(35.5˚N~37.5˚N, 113.5˚~116.5˚E)地磁场总强度空间分布的刻画精度,本文基于2025年4月35个地磁测点的总强度日变通化数据,系统对比了克里金插值(Kriging)与反距离权重插值(IDW)的空间建模性能。通过交叉验证法计算残差,结合平均绝对误差(MAE)、均方根误差(RMSE)和确定系数(R2)定量评估精度。结果表明:两种方法均呈现地磁场强度西北高、东南低的趋势,但IDW对局部极值(如孤立峰谷)表现更突出,而Kriging结果更平滑;Kriging的MAE (80.66 nT)、RMSE (94.65 nT)显著低于IDW (MAE = 127.30 nT, RMSE = 159.27 nT),且R2更高(0.90 > 0.71),表明其模型拟合度更优;IDW残差离散性强(−357.68~388.93 nT),存在明显区域偏移;Kriging残差分布集中(−200~200 nT),近似正态分布,验证其更有效捕捉空间相关性。综上,在冀鲁豫地区地磁场建模中,克里金插值精度与可靠性优于反距离权重插值,可为区域磁异常识别和深部构造研究提供更可靠的数据基础。未来需进一步优化插值参数并拓展多方法对比。
Abstract: To improve the precision of spatial distribution characterization of the geomagnetic total intensity in the juncture area of Hebei, Shandong and Henan (35.5˚N~37.5˚N, 113.5˚E~116.5˚E), this study systematically compares the spatial modeling performance of Kriging interpolation and Inverse Distance Weighting (IDW) interpolation based on the total intensity diurnal variation data from 35 geomagnetic measurement points in April 2025. Cross-validation is used to calculate residuals, and the accuracy is quantitatively evaluated using Mean Absolute Error (MAE), Root Mean Square Error (RMSE), and Coefficient of Determination (R2). The results show that both methods exhibit a trend of higher geomagnetic intensity in the northwest and lower in the southeast. However, IDW highlights the local extrema better (e.g., isolated peaks and valleys), while Kriging results are smoother. Kriging has significantly lower MAE (80.66 nT) and RMSE (94.65 nT) compared to IDW (MAE = 127.30 nT, RMSE = 159.27 nT), with a higher R2 (0.90 > 0.71), indicating better model fitting. IDW residuals show strong dispersion (−357.68~388.93 nT) with obvious regional shifts, whereas Kriging residuals are concentrated (−200~200 nT) and approximately normally distributed, validating its better capture of spatial correlation. In conclusion, Kriging interpolation demonstrates higher precision and reliability than IDW for geomagnetic modeling in the juncture area of Hebei, Shandong and Henan, providing a more reliable data foundation for regional magnetic anomaly identification and deep structural studies. Future work should focus on optimizing interpolation parameters and expanding multi-method comparisons.
文章引用:李博. 反距离权重插值与克里金插值在冀鲁豫交界区地磁场数据分布中的对比研究[J]. 地球科学前沿, 2025, 15(8): 1187-1194. https://doi.org/10.12677/ag.2025.158110

参考文献

[1] 张海洋, 李博, 苏树朋, 等. 冀鲁豫交界区岩石圈磁场异常分布与磁化率结构分析[J]. 地震地磁观测与研究, 2019, 40(6): 47-52.
[2] 王红军, 徐浩然, 赵文元. 电磁频谱地图构建中的空间插值方法综述[J]. 电讯技术, 2025, 65(1): 141-151.
[3] 李海涛, 邵泽东. 空间插值分析算法综述[J]. 计算机系统应用, 2019, 28(7): 1-8.
[4] 杨雨亭, 尚松浩, 李超. 土壤水分空间插值的克里金平滑效应修正方法[J]. 水科学进展, 2010, 21(2): 208-213.
[5] 张优, 王娟, 张杰, 等. GIS与地统计学的土壤水分空间插值方法[J]. 四川师范大学学报(自然科学版), 2019, 42(5): 703-710.
[6] 王怀计, 祝俊慧, 胡莹, 等. 顾及特征水深点距离重分配的反距离加权插值算法[J]. 海洋通报, 2022, 41(1): 9-18.
[7] 梁玉柔, 伍红玲, 王伟鹏, 等. PM2.5时空序列缺失数据的反距离权重插值方法补缺研究[J]. 环境与职业医学, 2025, 42(2): 171-178.
[8] 徐志萍, 徐顺强, 姜磊, 等. 冀鲁豫交界区流动地磁场总强度异常特征及地震活动性分析[J]. 山西地震, 2017(3): 28-33, 39.
[9] 王朝景, 李博, 苏树朋, 等. 2018年冀南地区岩石圈磁场变化分析[J]. 高原地震, 2020, 32(3): 9-12, 29.
[10] Su, S., Chen, S. and Zhao, H. (2022) Taylor Polynomial Spatial Reference Field Method for Field Geomagnetic Diurnal Variation Reduction. Izvestiya, Physics of the Solid Earth, 58, 981-991. [Google Scholar] [CrossRef
[11] 陈琳, 任春颖, 王宗明, 等. 基于克里金插值的耕地表层土壤有机质空间预测[J]. 干旱区研究, 2017, 34(4): 798-805.
[12] 冯波, 陈明涛, 岳冬冬, 等. 基于两种插值算法的三维地质建模对比[J]. 吉林大学学报(地球科学版), 2019, 49(4): 1200-1208.
[13] 付艳平, 幸松, 王妍, 等. 基于三种空间插值方法的气温适应性分析——以江西省为例[J]. 绿色科技, 2023, 25(4): 51-56.
[14] 陈鹏, 邓飞, 刘思廷. 三维空间属性插值方法的研究[J]. 电脑知识与技术, 2015, 11(7): 235-239.
[15] 刘光孟, 汪云甲, 王允. 反距离权重插值因子对插值误差影响分析[J]. 中国科技论文在线, 2010, 5(11): 879-884.