2. 日变改正的物理基础

日变改正是磁法勘探中把“随时间变化的磁场”从观测数据中剥离出去的关键环节。要正确实施改正，必须先理解三件事：日变从何而来、日变站与测点之间是什么关系、以及用什么数学模型把这种关系表达出来。本节围绕这三点展开，并给出可直接用于野外设计的参数判断与精度控制方法。

一、地磁场日变的物理成因

地磁场的日变化主要来自电离层中电流体系的昼夜运动。太阳紫外辐射使高层大气电离，形成导电层；在潮汐风与地球磁场共同作用下，电离层中产生电流涡旋，其强度随地方时、纬度和太阳活动水平变化，从而在地面产生量级约几纳特到数十纳特的规则日变。此外还存在太阳静日变化 Sq 与扰动日变（磁暴、亚暴）两类成分：前者周期稳定、可预测性强；后者幅度大、形态不规则，是野外作业中最需要警惕的干扰源。

从频率特征看，日变属于长周期信号，主周期为 24 小时及其谐波，典型振幅在中低纬度约 10～40 nT，磁静日更小，磁扰日可达上百纳特。这一特征决定了日变改正本质上是一个“低频趋势扣除”问题，也决定了改正效果高度依赖于日变站与测点之间日变形态的一致性。

二、日变站与测点磁场的关系

设日变站位置为 B，测点为 P。在某一时刻 t，测点观测值可分解为：

T_P(t) = T_P0 + ΔT_P(t) + δ_P(t)

其中 T_P0 为测点区域正常场与局部异常之和（即我们真正想保留的部分），ΔT_P(t) 为测点处的日变分量，δ_P(t) 为观测噪声与其它干扰。日变站同步记录 T_B(t) = T_B0 + ΔT_B(t)。日变改正的核心假设是：同一时刻，测点与日变站的日变分量近似相等或存在稳定的比例关系，即 ΔT_P(t) ≈ k·ΔT_B(t)。

这一假设的成立程度取决于两站间的距离、纬度差与地下电性结构差异。工程实践中通常要求日变站与测点距离控制在 50～100 km 以内，且尽量处于同一纬度带、同一构造单元；当测区跨度大或跨越不同纬度时，应布设多个日变站分区改正，而不是用一个站“一刀切”。

三、日变改正的数学模型

最常用的改正模型是线性比例模型：ΔT_P(t) = k·ΔT_B(t) + b，其中 k 为比例系数，b 为常数偏置。当测点与日变站距离较近、日变形态一致时，取 k=1、b=0 即可，即直接以日变站的相对变化量作为测点日变改正量。当两站日变幅度存在系统性差异时，可用重叠观测时段的最小二乘拟合求 k 与 b：

# 以重叠观测段拟合比例系数 k（示例数据，非实测）
# t: 时刻; dB: 日变站相对变化; dP: 测点相对变化
import numpy as np
k, b = np.polyfit(dB, dP, 1)      # 一次多项式拟合
# 若 |k-1| < 0.1 且残差标准差 < 2 nT，可采用 k=1 简化模型
dP_corr = dP - (k * dB + b)       # 改正后测点相对变化

实际作业中还需注意两点：其一，日变站记录本身包含仪器漂移，应先用台站已知基线或重复观测点做漂移校正，再参与改正；其二，测点读数时刻与日变记录时刻必须严格对齐，通常要求同步误差小于采样间隔的一半，否则在日变梯度较大的清晨与傍晚会引入明显误差。

四、改正精度的主要影响因素与判断

日变改正精度并非固定值，而是由多个因素共同决定。工程上可用“改正后同一测点重复观测的均方根差”作为直接检验指标：在测区内选取若干检查点，在不同时刻重复观测，改正后若各次结果的标准差小于设计精度（如 ±2 nT），说明日变改正有效；若超限，应排查站距、磁扰、时间同步或仪器漂移问题。

综上，日变改正的物理基础可归纳为：以电离层电流体系产生的低频日变为对象，以日变站与测点日变分量近似一致为假设，以线性比例模型为数学表达，以站距、磁扰、时间同步与仪器漂移为主要精度控制点。理解这一链条，才能在野外布站、数据预处理与质量检验中做出合理决策。