[发明专利]一种面向分层大地的视在电阻率测算方法、系统及存储介质

权利要求书

1.一种面向分层大地的视在电阻率测算方法，其特征在于，包括以下步骤：

S1、根据温纳四极法，初始化设置极距a，对电流极输入测试电流I，获取对应极距a的测试电压ΔV，计算视在电阻率测量值ρ_am；

S2、引入欧拉变换法获得两个电压极间的地表电位差ΔV′(a)，根据地表电位差ΔV′(a)计算分层大地对应极距a的视在电阻率计算值ρ_as；其中，采用欧拉变换法对地表电位差ΔV′(a)进行计算时，对极距a对应的参数以预设的间隔及预设的数量进行采样；

S3、计算测量值ρ_am和计算值ρ_as间的相对误差ε，根据相对误差ε是否超过预设阈值确认是否更新极距a：若超过预设阈值，则确认更新极距a并转入步骤S1；若不超过预设阈值，则不更新极距a，结束计算并输出视在电阻率计算值ρ_as。

2.根据权利要求1所述的一种面向分层大地的视在电阻率测算方法，其特征在于，所述步骤S1中，测试电流I为0.1mA-10mA。

3.根据权利要求1所述的一种面向分层大地的视在电阻率测算方法，其特征在于，所述步骤S1中，所述极距a为1m以上。

4.根据权利要求1所述的一种面向分层大地的视在电阻率测算方法，其特征在于，所述步骤S2中，所述欧拉变换法包括以下步骤：

S2.1、对应极距a构建关于地表电位差ΔV′(a)的模型：

其中，λ为积分变量，表示电磁场传播系数；ρ₁表示水平分层大地第1层大地的电阻率；J₀为第一类零阶贝塞尔函数；R₁为表达地表特征阻抗的函数，其形式如下：

其中，R_i表示水平分层大地第i层的特征阻抗；ρ_i表示第i层大地的电阻率；h_i表示第i层大地的分层厚度；m表示水平分层大地最深层所在层数；

S2.2、对式(1)进行汉克尔变换，所述汉克尔变换结果表达式为：

其中，J₀为第一类零阶贝塞尔函数；

S2.3、引入欧拉变换法，所述式(3.1)和式(3.2)变换结果为：

其中,λ＝e^-y；a＝e^x；x,y∈(-∞,∞)；

S2.4、按间隔Δ对式(4)中e^x-y进行采样，所述式(4)输出为离散信号如下：

其中，H_n为滤波系数；Δ_n为采样点位置；此时积分核函数f＝R₁；

S2.5、采取90个采样点形式，对式(5)进行截断求和，更新地表电位差ΔV′(a)模型，其表达式为：

其中，Δ_k为预设的位置系数，H_k为预设的权重系数。

5.根据权利要求4所述的一种面向分层大地的视在电阻率测算方法，其特征在于，所述步骤S2.4中，所述采样间隔Δ小于0.5倍积分核函数f的截止频率。

6.根据权利要求4所述的一种面向分层大地的视在电阻率测算方法，其特征在于，所述步骤S2.5中预设的Δ_k和预设的H_k如下表所示：

7.根据权利要求1所述的一种面向分层大地的视在电阻率测算方法，其特征在于，所述步骤S3中，所述相对误差ε表达式为：