[发明专利]一种大回线源瞬变电磁全域视电阻率确定方法及处理终端

说明书

本发明涉及一种大回线装置瞬变电磁全域视电阻率确定方法及处理终端，所述方法包括如下步骤：步骤1：将大回线分割成若干小段，计算各电偶极子的理论二次磁场，将各个理论二次磁场累加，得到整个大回线的理论二次磁场；步骤2：预设参考电阻率下，将理论二次磁场转换为时间上的离散序列，以离散序列的理论二次磁场作为节点、时间作为函数值，构造三次样条差值函数；将待求大回线的二次磁场Bc代入三次样条差值函数，得到二次磁场Bc对应的时间tx；步骤3：获得二次磁场Bc对应的采样时间tc，该大回线的视电阻率根据tc、tx和预设参考电阻率按公式计算得到。本发明的视电阻率是全域视电阻率，能够准确反映地下介质真实电阻率。

技术领域

本发明涉及确定视电阻率技术领域，具体是一种大回线源瞬变电磁全域视电阻率确定方法及处理终端。

背景技术

视电阻率是电磁勘探方法中的重要组成部分，其对勘探具有重要的意义。目前，在瞬变电磁领域定义视电阻率有两类方法，其一是基于磁场对时间的偏导数(也即是感应电动势)，其二是直接基于磁场。无论哪种方法，在现有文献中，关于瞬变电磁视电阻率的计算，主要集中于中心点位置的晚期或全时期定义，对于非中心点的视电阻率计算涉及极少。相比利用中心点晚期近似公式计算非中心点的视电阻率，直接计算非中心点视电阻率的方法更精确。在定义视电阻率的两类方法中，基于感应电动势定义的视电阻率有其固有缺陷：感应电动势和电阻率之间不是单调函数，基于此定义的视电阻率可能会出现双解，并且可能存在假极值(Raiche，1983；赵越等，2015)。目前对二次磁场定义的大回线非中心视电阻率的研究太少，因此，需要不区分接收点是否在大回线中心点、不区分接收数据是在二次磁场早期还是晚期的视电阻率确定方法，也即是需要全域视电阻率确定方法。

发明内容

针对现有技术的不足，本发明的目的之一是提供一种大回线源瞬变电磁全域视电阻率确定方法，其能够解决全域内确定视电阻率的问题；

本发明的目的之二是提供一种处理终端，其能够解决全域内确定视电阻率的问题。

实现本发明的目的之一的技术方案为：一种大回线源瞬变电磁全域视电阻率确定方法，包括如下步骤：

步骤1：将大回线分割成N个小段，使每段线的长度小于收发距的若干倍，从而将每个小段视为一个电偶极子，均匀半空间情况下，第i个电偶极子的理论二次磁场记为i＝1,2,…,N，

整个大回线产生的理论二次磁场B_z为公式①：

理论二次磁场为关于自变量为电阻率ρ和时间t的函数，从而公式①可变换为公式②：

B_z＝f(ρt)------②

式中，ρ表示电阻率，t表示电流关断后的时间；

步骤2：根据公式②，预设参考电阻率ρ₀，参考电阻率ρ₀下的接收点的理论二次磁场B_z转换为时间t[j]上的离散序列的理论二次磁场B_z[j]，B_z[j]表示第j个f(ρt)函数值，t[j]表示第j个时间，j＝1,2,3,…，t[j]和B_z[j]一一对应；

以理论二次磁场B_z[j]作为节点，以时间t[j]作为函数值，构造三次样条差值函数，

获得待求视电阻率的大回线的二次磁场B_c，二次磁场B_c为差值点代入三次样条差值函数，计算出二次磁场B_c对应的时间t_x；