集成UAV磁测、Monte Carlo不确定性传播与DBSCAN聚类的高分辨率海岸磁异常检测与环境评估工作流——以葡萄牙Cascais为例

《Results in Engineering》：An Integrated UAV Magnetometry Workflow for Magnetic Anomaly Detection and Environmental Assessment in Cascais (Portugal)

　　

研究人员提出用于葡萄牙Cascais海岸地下调查的高分辨率无人机（UAV）磁测工作流。Monte Carlo（MC）不确定性传播提高了反演解释的可靠性。地质统计学克里金（kriging）能够实现磁场精确重建。DBSCAN可自动识别连贯的磁异常簇。可重复工作流支持环境与基础设施监测。

Miriam Martínez Mu?oz、André Garcia Sá、Machiel Bos和Isabel Cristina Gil García的研究关注浅层磁性源识别与海岸环境地球物理评估。研究背景是：磁测非侵入且对铁磁材料和磁性质差异敏感，但传统地面调查受通达和地形限制，载人航空调查距地表过远，难以分辨浅部小目标。无人机（UAV）磁测通过低空采集、灵活航线和密集采样弥补尺度差距，已用于近地表铁质目标、未爆弹（UXO）相关目标、矿产和考古。但测量质量受平台电磁干扰和悬挂传感器运动/姿态变化影响，且磁解释非唯一，异常定位不能等同于确认埋藏物或管线。Cascais海岸城市发展、交通和海岸工程与异质海岸材料共存，是挑战性测试环境。研究人员因此评估集成工作流，目标包括记录可用UAV磁采集和定位数据，用地统计工具表征空间磁变异，识别空间连贯磁异常候选，并区分异常检测/插值性能与物理源验证。由于无独立地面真值，源解释被作为需后续确认的假设。研究结论强调，保留档案可验证2021年6月9日UAV磁测四个MagDrone R3记录、2,049,937个原始磁样本、200 Hz采样、约5 Hz全球导航卫星系统（GNSS）定位和三轴加速度数据；原始分析报告普通克里金（ordinary kriging）RMSE=4.2 nT、RMSSE=1.02、相关=0.987，以及七个异常簇，但精确阈值、eps、min_samples、缩放和bootstrap设置未保存，不能完全重现。检测特征是磁异常候选，不是确认基础设施。意义在于支持UAV磁测作为复杂海岸环境快速筛查工具，并强调完整参数归档和未来与公用记录、地面磁测、探地雷达（GPR）、电阻率层析成像（ERT）等独立观测配对。

关键技术方法（≤250字）：研究采用DJI Matrice 600 Pro多旋翼，搭载MagDrone R3和双Bartington Mag-03MSES70三轴磁通门磁力仪，非磁悬挂，集成GNSS定位和三轴加速度通道。样本队列来源为2021年6月9日葡萄牙Cascais四个MagDrone R3记录，共2,049,937个原始磁样本，200 Hz磁采样，约5 Hz GNSS更新。处理使用中位绝对偏差（MAD）异常检测，磁补偿、时间校正和去趋势，变差函数与普通克里金，基于密度的带噪声空间聚类（DBSCAN）和轮廓系数，交叉验证及不确定性评估。

研究结果：
3.1 研究区与采集证据：研究区位于葡萄牙西部Cascais沿海，坐标约38.67803–38.69592° N和9.46047–9.41876° W。四个记录含2,049,937个原始磁样，200 Hz；GNSS约5 Hz；GNSS高度约10.0–66.9 m，不等于离地高度（AGL）。典型运动地速约9.8–10.0 m s^?1。加速度标准差约0.075–0.095 g（X）、0.116–0.170 g（Y）、0.335–0.399 g（Z），支持运动/振荡是相关噪声机制。档案可验证采样频率、坐标范围等，但不足以重建全部后续处理决策，因此“异常高空间分辨率”被缓和。
3.2 处理后磁场特征：补偿、时间校正、去趋势和不确定性评估后，总磁场强度（TMI）约43,686–44,124 nT，总范围约438 nT。原始分析报告残余磁场标准差约28.4 nT，但不等同仪器不确定性；组合不确定性1.5±0.5 nT是误差预算，二者概念不同。
3.3 变差函数与空间结构：实验变差函数显示清晰空间相关，球状模型最满意。块金约18 nT²，偏基台约1,247 nT²，实际变程约42.3 m。方向变差函数仅有限各向异性，故采用全向球状模型。
3.4 克里金插值与三维建模：普通克里金生成连续磁场。交叉验证ME=?0.3 nT，RMSE=4.2 nT，RMSSE=1.02，相关=0.987，表明原分析插值内部一致，但不验证物理源。三维表面显示多个局部峰。
3.5 自动异常检测：原始工作流报告4,712个异常节点、23个候选组、7个保留簇，平均轮廓系数0.73。因处理输入表和DBSCAN参数未保留，不能独立重现，也不作为bootstrap确认的高置信选择。研究区分测量磁变异、空间/插值和聚类结果、物理源验证三级证据，仅有前两级。
3.6 磁异常簇解释：结合几何、振幅、簇持续性、地统计和环境背景，解释为可能源类而非确定识别。A1紧凑高幅，可能为局部磁性源或埋藏金属碎片；A2空间足迹较广，可能反映源尺寸、深度、几何或邻近多源；A3紧凑浅部，可能为人为铁磁源；A4和A5拉长，可能与线性源如公用走廊、管道、电缆、工程特征或拉长地质/人为反差相容。无独立野外观察确认材料或功能。
3.7 讨论：
3.7.1 与先前UAV磁测研究比较：先前研究覆盖矿产、近地表和考古目标；Cascais增加海岸/城市案例，但比较需区分磁异常检测与确认目标识别。主要方法贡献是集成采集、地统计插值、不确定性考虑和密度聚类；单个组件不新，价值在组合及明确可重现性边界。
3.7.2 提出方法学的优势：UAV磁采集可与空间建模和自动异常分组结合。修订缓和“完全可重现框架”，原始采集文档较好，但若干处理参数未归档。可转移性是工作流层面，贡献是明确链接采集证据、插值诊断、异常分组和不确定性/验证限制。
3.7.3 研究局限：无独立地面真值，无公用地图、GPR、地面磁测或挖掘等确认源身份；仅依赖磁测；磁解释非唯一；平台电磁干扰和运动诱导效应均可能贡献。结果与低空UAV磁测可分辨局部磁异常的文献一致，新颖性限于集成案例工作流和证据层级分离。
3.7.4 未来研究方向：应集成UAV磁测与GPR、ERT、电磁感应等互补地球物理，探索机器学习异常分类和重复UAV调查的时间监测，以增强UAV地球物理在环境管理、基础设施监测和海岸可持续评估中的作用。

讨论部分总结：研究人员认为，UAV磁测适合复杂海岸环境快速筛查，但成功采集不等于物理源确认；插值统计、聚类稳定性与源验证必须分开报告；完整参数归档对可重复性关键。研究结论翻译：保留的Cascais档案验证了2021年6月9日UAV磁测，包括四个MagDrone R3记录、2,049,937个200 Hz原始磁样本、约5 Hz GNSS定位和三轴加速度数据。档案坐标定义调查海岸段，GNSS和加速度通道提供采集几何及传感器/平台动力学直接证据。原始地统计分析报告普通克里金RMSE=4.2 nT、RMSSE=1.02、观测与预测相关=0.987，这些仅表征插值性能。原始DBSCAN工作流还报告七个保留磁异常簇，但精确阈值、eps、min_samples、缩放、bootstrap设置和处理输入表未保存，因此不能描述为从保留档案完全可重现。最重要的是，无独立地面真值数据集。检测特征因此是磁异常候选，提出的源类为假设，不是确认埋藏基础设施。早期版本中的数值源深度和磁矩同样需要原始源模型和方程，才能视为可重现估计。研究支持UAV磁测作为复杂海岸环境空间磁变异快速筛查工具，同时显示完整参数归档的重要性。未来工作应将UAV磁测与公用记录、地面磁测、GPR、ERT或其他独立观测配对，并归档端到端重现所需所有处理脚本和参数。