引言
高密度电法是一种常用的地球物理勘探方法，广泛应用于地下水资源探测、矿产资源勘探和工程地质调查等领域。反演算法是高密度电法数据处理的核心，直接影响反演结果的精度和可靠性。最小二乘法和佐迪反演法是两种常用的反演算法，各有优缺点。本文针对高阻模型，对比研究了这两种方法的反演效果，旨在为实际应用提供参考。
01 基本原理
1.1 最小二乘法原理

最小二乘法是经典的数据处理方法，核心是通过最小化观测与模型计算数据的误差平方和确定模型最优参数。在地球物理反演中，用于推断地下介质物性参数分布。

在反演中，设观测数据向量d（维度N）、模型参数向量m（维度M）、模型响应向量F(m)，误差向量 。最小二乘反演目标是使e的二范数最小，目标函数。

对于线性反演问题，（G为灵敏度矩阵），目标函数，求解正规方程可得模型参数解

实际中多为非线性问题，采用迭代线性化。将F(m)在初始模型处泰勒展开，保留一阶项，为雅可比矩阵）。代入目标函数求导并令导数为零，得到迭代线性方程组，求解得模型参数更新量，不断迭代至目标函数收敛，得到反演结果。

1.2 佐迪反演法的基本原理

佐迪反演法是基于电磁测深曲线的迭代拟合法，通过调整初始模型，使实测与模型拟断面差异最小，获取地下地质结构电阻率分布模型。

视电阻率数据是反演关键，它受地下介质多种因素影响。其数学模型基于正演模拟，常采用有限差分法、有限元法等数值方法。
设地下地质模型参数为m，正演模拟函数F(m)为理论视电阻率值，实际观测数据为d。佐迪反演目标是找最优m使两者差异最小，常用目标函数为均方误差。

不断调整m使目标函数减小至收敛得反演结果。因地下结构复杂和噪声干扰，反演多解，常结合地质构造、岩石物性等先验信息约束反演，如已知高阻岩层，可在初始模型设置或约束高阻层参数，减少不确定性。
02 高阻模型建 立与分析
 2.1 模型装置与参数
高密度电阻率法装置的类型有许多：温纳装置（Wenner Alpha）、偶极装置（Wenner Beta）、斯龙贝格装置（Schlumber）等装置。这些装置都可以进行数据采集，但这些方法各有其优缺点和相应的限制，因此在实际的野外工作中，应根据地面电气条件、地形条件及区域地质条件等具体情况来解决这些问题。并选择合适方法装置进行工作 。本次模型建立只针对的温纳装置(图2.1)进行模拟和分析。

装置系数K=2^πna，A M=M N=N B=n a（n为隔离系数，a为电极极距）。测量时，A M=M N=N B是电极之间的距离，A、B、M和N同时向右移动并得到第一段线。然后，AM、MN、NB增加一个电极间距并逐点移动到另一个点，以获得另一个剖面线，从而连续地扫描，得到梯形断面。

在高密度电阻率法探测中，采集参数直接影响到数据的质量，而数据是整个探测过程的根本，因此采集参数选择的是非常重要的。具体的参数包括：设备装置类型，极距，所用电极个数，隔离系数（包括最小隔离系数和最大隔离系数），电源供电方式。其中，根据所选择的采集设备的差异，供电模式的选择是合理的模式。选择的原则是，在采集过程中电流和电压值不应过大或过小。通常在高电阻率地区采用恒压模式，低电阻率地区采用恒流模式。极距、使用电极总数和隔离系数共同反映的是采集范围的大小。水平距离等于极距*使用电极的总数，但是我们应该考虑采集深度的水平的变化。采集深度＝极距×（最大隔离系数－最小隔离系数）。
 2.2 正演模型建立与分

采用res2dmod软件建立背景模型大小为：120×40m（长×宽），电阻率300Ω.m，三个高阻模型的电阻率为10000Ω.m：2×2m、4×4m、8×8m的矩形，模型位于测线中间附近位置，中心深度10m。

正演观测装置采用温纳装置60道电极、2m极距、最大隔离系数19层（确保高阻模型在其有效探测范围内）。

通过对比正演模拟结果（图2.2-2.4）可以看出：高阻模型2×2m边长与观测2m极距相等，温纳装置正演无法显示其真实的位置与大小，产生层状的假高阻异常；高阻模型4×4、8×8m，温纳装置都能显示出高阻异常位置与大小，观测异常边界相对于原始模型范围较大些，异常中心深度5m左右相当于理论模型一半的深度。

2.3 最小二乘法反演
采用Res2dinv反演软件对于上述三个高阻模型进行最小二乘法反演。通过对反演模拟结果（图2.5-2.7）可以看出：高阻模型2×2m，温纳装置反演无法显示其真实的位置与大小，产生圆形闭合的假高阻异常与实际结果不相符；高阻模型4×4、8×8m，温纳装置都能显示出高阻异常位置与大小，观测异常边界（圆形）相对于原始模型范围（矩形）相差不多，异常中心深度8m左右，相当于理论模型0.8倍深度。

2.4 佐迪反演法  
采用数据转换软件把模型正演数据写成骄佳软件能够读取的格式，然后采用佐迪反演法进行计算分析。

通过对反演模拟结果（图2.9-2.11）可以看出：高阻模型2×2m，温纳装置反演无法显示其真实的位置与大小，产生假高阻异常与实际结果不相符；高阻模型4×4、8×8m，温纳装置都能显示出高阻异常的位置与大小，观测异常边界相对于原始模型范围相差较大，整体规模相当于原始高阻模型的1.5-2倍左右，异常中心深度10m左右，相当于理论模型深度最为接近。

2.5 小结  

1.高阻模型边界小于等于极距大小的时候，无论是正演还是反演在有效探测范围内都无法对于异常进行有效探测而且都会产生假异常。

2.高阻模型边界大于等2倍极距大小的时候，正反演都能产生有效的异常相应，正演：深度大致相当于与实际深度的一半左右。两种反演方法都能高阻异常位置与大小：最小二乘法观测异常边界（圆形）相对于原始模型范围（矩形）相差不多，相当于理论模型0.8倍深度左右；佐迪反演法观测异常边界相对于原始模型范围相差较大，整体规模相当于原始高阻模型的1.5-2倍左右，异常中心深度10m左右，相当于理论模型深度最为接近。

03 实例应用

本次试验目的是验证高阻采空区正演模拟结果的准确性。采用重庆顶峰地质勘探仪器有限公司生产的EDGMD级联式高密度电法测量系统，布设电极60道，极距3m，并且根据测线以及测道数选取最大隔离系数19层，供电断电采用0.5s，最大供电电压600V，单次测量剖面数据点为570个。

已知采空区位置为测线长度102m、深度12m，采用最小二乘法进行反演，圈定异常位置中心点在测线长度102m、深度4m处；佐迪反演法圈定异常位置的中心点在测线长度102m、深度12m处。本次实际测量数据的处理结果与正演模型的反演结果相同，实验结果准确。

04 结论

①.最小二乘法反演对水平位置和高阻异常形态的刻画较为准确，在不人为调节深度系数的情况下，目标体的深度过度压缩，造成测量结果不准确。

②.佐迪法反演对高阻目标体的水平位置及深度的刻画都较为准确，但该方法会对异常范围进行放大，从而影响对目标体的准确判断。

③.用最小二乘法和佐迪法对高密度数据进行反演，将两种方法的反演结果互相佐证，得出的最终结果会更为准确。