朱姣， 殷长春， 任秀艳， 刘云鹤， 惠哲剑， 谷宇

吉林大学地球探测科学与技术学院， 长春 130026

0 引言

直流电阻率法是电法勘探领域一个重要的分支已在环境和工程领域获得广泛应用.该方法利用接地电极向地下供电、建立地下稳定电流场，通过观测测量电极间的电位差以达到推断和解释地下目标体的目的.在野外勘探中我们常常采用剖面或者测深装置进行高密度电阻率测量.为此，需要在地面布设数十至数百个测量电极，通过多芯电缆连接到观测设备，从而实现多种电极装置的分布式测量(Loke and Barker, 1996).由于环境工程领域勘探目标大多为三维地质体，采用这种阵列测量方法可以准确地重构出三维地质体的空间分布，因此高密度电阻率法可获得丰富的地质信息(Loke et al., 2013).

随着高密度电阻率法在环境、工程、地下水资源等勘探领域的广泛应用，三维反演和解释技术研究受到广泛关注(Udphuay et al., 2011; Kneisel et al., 2014; Neyamadpour, 2019; Nthaba et al., 2020).同时，随着电阻率法向三维精细化探测方向发展，能够提取更多的细节信息，为数据解释打下良好的基础.然而，对精细结构的刻画将导致庞大的三维数据体，对这些三维数据体进行反演解释需要高效和高精度的三维正演模拟技术(Günther et al., 2006; Rücker et al., 2006).常用的电阻率三维模拟算法包括积分方程法(IE)、有限差分法(FD)和有限元法(FE)(Wang and Signorelli, 2004; Yoon et al., 2016; Ren et al., 2018a).积分方程法仅对异常体进行剖分，计算效率高，但不适合复杂异常体的数值计算(Sheng et al., 2006; Ueda and Zhdanov, 2006).有限差分算法简单、求解速度快，但由于要求计算域被剖分成规则单元，对……