4. 视电阻率曲线的数据质量检查与预处理
野外采集到的视电阻率曲线,在进入反演或定性解释之前,必须先经过一轮系统的质量检查与预处理。原始曲线中混杂着电极接地不良引起的跳点、地形起伏造成的假异常、随机噪声导致的毛刺,以及个别测点因耦合失效产生的畸变段。如果直接把这些数据送入反演,轻则分层界面模糊,重则出现虚假的低阻或高阻异常体,误导地质推断。本节围绕“如何判断曲线是否可信”和“如何把不可信的部分处理到可接受”两条主线,讲清楚平滑、剔点、地形校正和接地电阻识别四类实操方法,并给出一套可执行的整理规范。
一、先看整体,再抠细节:质量检查的顺序
拿到一条测线的视电阻率曲线,第一步不是急着平滑,而是先做整体浏览。把同一测线各测点的ρs曲线按测点号排列,观察三件事:曲线首尾是否收敛、相邻测点之间形态是否连续、是否存在整段平移或突跳。相邻测点在同一极距上的ρs值通常变化平缓,若某测点整体抬高或压低一个数量级,多半是供电电流记录错误或电极接地电阻突变,而不是真实地质变化。第二步才是逐点检查单条曲线的局部形态,识别毛刺、台阶和尖峰。这个“先整体后局部”的顺序能避免把区域性系统误差误判为随机噪声而错误平滑掉。
二、曲线平滑:压住毛刺,不抹掉真实异常
平滑的目的是压制高频随机噪声,但必须守住一条底线——不能把有意义的小幅异常抹平。工程上常用滑动平均和中值滤波两种方式。滑动平均对高斯型噪声效果好,但对孤立尖峰会把尖峰“摊开”,反而污染邻近点;中值滤波对脉冲型干扰更稳健,适合先用来压制跳点。实际操作中,建议窗口取3点或5点,且只对确认属于随机噪声的区段施加平滑,对曲线拐点、极值附近应保留原始值或减小窗口。判断平滑是否过度,可以对比平滑前后曲线的极值位置和幅值:若某个局部高阻峰被削低了20%以上,就要怀疑平滑窗口过大。平滑后的曲线应另存为派生数据,原始曲线必须保留,以便追溯。
三、畸变点剔除:用统计准则加人工复核
畸变点通常表现为单点或连续两三点相对邻点出现数倍以上的突跳。可先用相邻比值法初筛:计算某点ρs与其前后各两点均值的比值,若偏离超过设定阈值(例如大于2.5或小于0.4),标记为可疑点。再用3σ准则在局部窗口内复核,减少误判。需要强调的是,统计准则只能“标记”,不能自动“删除”。野外常见的陡倾角薄层、断裂带附近本来就会出现真实的大梯度变化,若一律按统计阈值剔除,会把真实构造信息一起删掉。正确做法是把可疑点标在曲线图上,结合地质草图、测点地形和电极接地记录逐点复核,确认属于观测错误的才剔除,剔除后在该极距位置用相邻点线性插值补值,并在处理记录中注明“剔除原因+补值方式”。
四、地形校正:把起伏地表的几何影响分离出去
视电阻率是在假定地表水平的前提下定义的。当测点位于山坡、沟谷或起伏较大的基岩裸露区时,电流场的实际分布与水平半空间假设不符,会在曲线上叠加地形假异常,典型表现是山脊处出现假高阻、沟谷处出现假低阻。地形校正的基本思路是:先根据测点高程和电极排列的实际位置,建立二维或三维地形模型,用有限元或有限差分正演计算纯地形影响下的视电阻率响应,再用实测ρs除以地形响应系数,得到校正后的视电阻率。实操中,若测线高差在极距的10%以内,可只做一维近似修正;高差更大或地形突变处,必须做二维校正。校正后要检查曲线形态是否更连续、假异常是否减弱,同时警惕校正过度引入新的畸变。
五、电极接地电阻影响识别:从供电回路找线索
电极接地电阻过大是野外最常见的质量隐患。当供电电极接地不良时,供电电流下降,测量电位差随之偏小,计算出的ρs会系统性偏低;若接地电阻不稳定,还会造成同一测点重复观测值离散。识别方法很直接:检查野外记录中的供电电流I和供电电压V,若某测点I明显低于同测线正常值(例如下降超过30%),或V异常升高,就应怀疑接地电阻问题。同时观察该点ρs是否整体偏低、重复观测是否超限。处理上,优先通过浇水、增加电极并联数量、换用低阻位置重新接地来重测;确实无法重测的测点,应在数据集中标记为“低置信度”,在反演时降权或剔除,而不能简单用平滑掩盖。
六、可信曲线数据集的整理规范
经过上述处理,每条曲线都应形成“原始数据—处理记录—派生数据”三件套。建议的数据集字段包括:测点号、极距、原始ρs、供电电流、接地电阻状态、是否剔除、剔除原因、平滑方式与窗口、地形校正系数、最终采用ρs、置信度等级。置信度可分三级:A级为原始质量良好或经重测确认;B级为经平滑或插值处理但形态可信;C级为存在接地或畸变疑点、仅作参考。整理时坚持三条原则:原始数据只读不改;每一步处理可追溯;同一测线内处理参数尽量统一,避免因处理方式不一致制造人为的横向假异常。只有把质量检查做在前面,后续的曲线类型划分和反演解释才有可靠基础。