2. 电极排列与装置系数核对
装置系数为什么必须逐点核对
视电阻率不是直接测量量,而是由实测电位差 ΔV、供电电流 I 和装置系数 K 换算得到:ρa = K · ΔV / I。K 只取决于四个电极的相对位置,与地下介质无关。因此一旦 K 用错,整条剖面的 ρa 会被同一个比例因子整体放大或缩小,表现为"系统性偏差"——曲线形态不变、异常位置不变,但幅值整体偏高或偏低。这正是装置系数错误最典型的指纹,也是与坏点、接地不良等随机误差区分的关键。
两类常用装置的 K 公式(四极 AMNB,A、B 为供电电极,M、N 为测量电极):
- 温纳装置(Wenner α):AM = MN = NB = a,K = 2πa。
- 偶极-偶极装置(Dipole-Dipole):供电偶极 AB = a,测量偶极 MN = a,两偶极间距为 n·a,K = π·n·a·(n+1)(n+2)。
逐点核对实例
设电极距 a = 5 m,测线电极编号 1~10,坐标 xi = 5·(i−1) m。取温纳排列 A=1、M=2、N=3、B=4(即 AM=MN=NB=5 m),以及偶极-偶极排列 A=1、B=2、M=3、N=4(n=1)和 A=1、B=2、M=4、N=5(n=2)。用通用式逐点计算并与公式值比对:
| 装置 | 电极 (A,M,N,B) | AM / AN / BM / BN (m) | 通用式计算 K (m) | 公式值 K (m) | 判定 |
|---|---|---|---|---|---|
| 温纳 α | 1,2,3,4 | 5 / 10 / 5 / 10 | 2π/(0.2−0.1−0.2+0.1)=31.42 | 2πa = 31.42 | 一致 |
| 偶极-偶极 n=1 | 1,2,3,4 | 5 / 10 / 5 / 10 | 31.42 | π·1·5·2·3 = 94.25 | 不一致 |
| 偶极-偶极 n=2 | 1,2,4,5 | 5 / 15 / 10 / 20 | 2π/(0.2−0.0667−0.1+0.05)=37.70 | π·2·5·3·4 = 376.99 | 不一致 |
上表第二、三行暴露了常见错误来源:把偶极-偶极的电极坐标按温纳的等间距方式填写(A、B、M、N 顺序错位),导致通用式算出的 K 与装置公式值相差 3 倍、10 倍。正确做法是:偶极-偶极中 A、B 相邻、M、N 相邻,两偶极之间隔 n·a。以 n=1 为例,A=1、B=2、M=3、N=4 时 AM=10、AN=15、BM=5、BN=10,通用式得 K = 2π/(0.1−0.0667−0.2+0.1) = 94.25 m,与 π·n·a·(n+1)(n+2) 完全吻合。
系统性偏差的识别与排查流程
拿到一条剖面后,若怀疑 K 有误,按以下顺序排查,可快速定位问题:
- 看偏差形态:整条剖面 ρa 与邻线或历史数据相比整体偏高/偏低,但异常位置、曲线形态一致 → 高度怀疑 K 错误。
- 抽查 K 值:从数据文件中取 3~5 个不同 n 值的测点,用通用式复算 K,与文件中的 K 列比对。若比值恒定(如恒为 3、10),即为装置系数错误。
- 核对电极坐标:检查数据头中 A、B、M、N 的坐标是否与野外实际排列一致,重点看偶极-偶极的 n·a 间距是否被误填为 a。
- 修正并重算:用正确 K 重算 ρa,对比修正前后剖面,确认异常幅值回到合理区间。
核对记录表(可直接填写)
| 测点号 | 装置 | n | A/B/M/N 坐标 (m) | 文件 K (m) | 复算 K (m) | 比值 | 结论 |
|---|---|---|---|---|---|---|---|
| W-01 | 温纳 α | — | 0/15/5/10 | 31.42 | 31.42 | 1.00 | 正常 |
| D-01 | 偶极-偶极 | 1 | 0/5/10/15 | 31.42 | 94.25 | 0.33 | K 偏小 3 倍 |
| D-02 | 偶极-偶极 | 2 | 0/5/15/20 | 37.70 | 376.99 | 0.10 | K 偏小 10 倍 |
表中比值一列是排查的核心判据:比值恒为 1 说明 K 正确;比值恒定但不为 1,说明整条剖面存在系统性比例误差,需按正确装置公式重算全部 K 值后再反演。完成修正后,建议用修正前后剖面的异常幅值对比作为验证——若高阻异常带的 ρa 峰值回到与地质背景相符的量级,即可确认修正有效。