一种基于极化介质模型的广义趋肤深度计算方法与流程

本发明涉及地球物理勘探领域，尤其涉及一种基于极化介质模型的广义趋肤深度计算方法。

背景技术：

趋肤深度是电磁勘探领域中一个非常重要的概念，通常将其用来估算电磁法的勘探深度。在实际电法勘探中，为了解释实测数据，需要计算视电导率与探测深度数值，并表示成视电导率-深度断面的形式，因此，需要计算在实测不同频率或不同时刻情况下的探测深度。

目前有关于电磁探测深度的研究，szarka、fischer(1989)提出，采用1.5倍趋肤深度来估计大地电磁的理论探测深度。yans(2009)采用解析法、有限差分时域法和时频法对电磁测深探测深度进行了研究。目的是为现场工作设计提供依据，以达到预期目标。薛国强(2014)以接地长导线源为例研究了可控源电磁场全域的有效趋肤深度。

cn108376204a公开了一种基于地下粗糙介质模型的电磁广义趋肤深度计算方法，在趋肤深度定义方式的基础上，根据粗糙介质广义电导率的定义,推导出适用于符合地下粗糙介质电导率的电磁广义趋肤深度公式，根据不同岩石层粗糙度参数，计算出相应的探测深度。但以上探测深度计算方法是基于均匀介质以及粗糙介质理论下的计算公式，而实际地下介质体结构极为复杂，尤其地下多孔极化介质或者多金属矿探测时，传统的趋肤深度定义不再适用，因此，重新定义基于分数阶电导率模型的广义趋肤深度具有重要意义。

技术实现要素：

本发明的目的在于针对传统均匀半空间模型的趋肤深度计算方法，已经不满足地下极化介质的频散特性等传播规律，根据极化介质的分数阶电导率，以及趋肤深度的定义方式，提供一种基于地下极化介质模型的电磁广义趋肤深度计算方法。

本发明在传统趋肤深度定义方式的基础上，根据极化介质电导率的定义，推导出适用于符合地下极化介质电导率的电磁广义趋肤深度公式，根据不同地下极化介质的电导率、极化率、充电率、频散系数，计算出探测深度。

本发明是这样实现的，一种基于地下极化介质的广义趋肤深度频域计算方法，包括如下步骤：

(1)先对实测数据进行噪声滤波、取样叠加处理，采用广义s变换将数据进行时频变换，获得频域形式的数据；

(2)采用反演法对步骤1中处理后的电磁数据，进行地下介质电导率、极化率、充电率、频散系数的提取；

(3)将柯尔-柯尔分数阶电导率模型iω的负分数阶项，进行有理化整理，获得多孔极化介质的电导率正分数阶简化表达式；

(4)电导率正分数阶简化表达式代入到波数表达式求解虚部和实部表达式，进行简化整理得出衰减常数及相移常数表达式；

(5)对衰减常数取实部求取地下极化介质的广义趋肤深度。

进一步地，步骤3中，根据柯尔-柯尔模型，地下极化介质分数阶电导率表示为将地下极化介质分数阶电导率表达式中的iω的负分数阶项，进行有理化整理，转化为正分数阶的简化形式，表示为：σ＝σ′+σ″(iω)^c，其中σ∞为高频交流电导率，τ为充电时间常数，η为极化率，c为频散系数，ω为角频率。

进一步地，步骤4中，根据波数表达式k²＝ω²με-iωμσ(ω)，其中ω为角频率，σ为电导率，ε为介电常数，步骤3中简化后的极化介质电导率代入波数表达式中整理虚部和实部，令i^c＝a+ib，k＝α+iβ，写为：α²-β²＝ω²με+ω^1+cμσ″b，2αβ＝-(ωμσ′+ω^1+cμσ″α)，进行整理化简，得到衰减常数α和相移常数β，表达式为：

进一步地，步骤5中，地下极化介质的广义趋肤深度为：

其中为趋肤深度，re(α)是代表取α的实部，α为衰减常数。

进一步地，

采用反演法从频域中提取地下介质电导率、极化率、充电率、频散系数包括：求解目标函数最小值，目标函数表示为：

其中，m代表的是原始模型，n为实测的频点个数；pαi为第i个频点的实测响应值，其为视电阻率、相位或者电场与磁场值，mj(j＝1，2，3，…m)为模型参数；pci是由预测模型正演计算所得到的第i个频率的响应参数，α＝2为范数；

将上述目标函数趋于最小，对于各个频点i(i＝1，2，3，…n)的目标函数均满足下面一组线性方程：对线性方程进行求解可以得到反演m的修改量δm，从而获得新模型，通过对线性方程的求解提取出地下介质电导率、极率化率、充电率等多参数。

有益效果：本发明与现有技术相比，在地下多孔极化介质或多金属矿结构进行电磁探测时，可以建立地下极化率与广义趋肤深度之间的关系，能够准确计算地下极化介质的趋肤深度参数。本方法为我国开展感应电磁探测寻找矿产资源提供准确的深度定位，更有利于提高感应电磁勘探方法的解释精度。

附图说明

图1是基于地下极化介质模型的电磁广义趋肤深度计算方法流程图；

图2是电磁实测数据预处理流程图；

图3是本发明一个实施例同极化率不同频率下，传统趋肤深度与极化介质广义趋肤深度对比结果图；

图4是本发明一个实施例,同极化率不同频率下，传统趋肤深度与极化介质趋肤深度相对误差(a)和绝对误差(b)。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图和具体实施方式，对本发明进行进一步详细说明。本发明的核心构思之一在于，发明提供了极化介质广义电导率的探测深度计算方法，能够对电磁探测数据进行准确解释，从而获得地质目标体的准确位置信息。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

实施例

结合图1所示，一种基于地下极化介质模型的电磁广义趋肤深度计算方法，包括：

(1)针对地下极化介质建立柯尔-柯尔分数阶电导率模型，根据测区实测数据获取电磁响应参数；

如图1所示，先对电磁实测数据进行处理，包括噪声滤波、频时转换等预处理；

预处理首先对原始数据进行去噪处理，由于电磁探测现场环境的复杂性，数据会受到很多因素干扰从而存在噪声，主要干扰因素包括工业干扰、无线发射塔，天气等。为了提高数据质量，要对噪声特点分析，并对实测数据进行噪声滤波。其次，由于极化电导率表达式是以频率形式给出，因此需要将数据转换到频域形式。

进一步，如图2所示的实测数据预处理流程图。先对原始数据进行噪声滤波，去除噪声干扰，再对数据进行频时转换，将数据转换到频域形式，检查原始数据是否合乎质量要求。经过预处理和质量统计合格的数据，可以进行极化率、电导率等多参数的提取。

(2)采用反演法对预处理后的电磁数据，进行地下介质电导率、极化率、充电率、频散系数的提取。

采用反演法从频域中提取地下介质电导率、极化率、充电率、频散系数，反演法可归结于求解目标函数最小值，目标函数表示为：其中m代表的是原始模型，n为实测的频点个数；pαi为第i个频点的实测响应值，其可以为视电阻率、相位或者电场与磁场值。mj(j＝1，2，3，…m)为模型参数；pci是由预测模型正演计算所得到的第i个频率的响应参数，α＝2为范数。将上述目标函数趋于最小，对于各个频点i(i＝1，2，3，…n)的目标函数均满足下面一组线性方程：对线性方程进行求解可以得到反演m的修改量δm，从而获得新模型，通过对线性方程的求解提取出地下介质电导率、极率化率、充电率等多参数。

(3)将柯尔-柯尔分数阶电导率模型的iω的负分数阶项，进行有理化整理，获得多孔极化介质的电导率正分数阶简化表达式。

进一步地，根据多孔极化介质广义电导率公式其中σ∞为高频交流电导率，τ为充电时间常数，η为极化率，c为频散系数，ω为角频率。将电导率表达式进行整理化简，将括号打开后，分子分母同时乘以σ＝(1-η)(iωτ)^c，

令

则：σ＝σ′+σ″(iω)^c。

(4)根据波数表达式求解虚部和实部表达式，进行简化整理得出衰减常数及相移常数表达式。

进一步地，根据波数表达式k²＝ω²με-iωμσ(ω)，其中ω为角频率，σ为电导率，ε为介电常数。简化后的极化介质电导率σ＝σ′+σ″(iω)^c代入波数表达式中整理虚部和实部，i^c＝a+ib，k＝α+iβ，写为：α²-β²＝ω²με+ω^1+cμσ″b，2αβ＝-(ωμσ′+ω^1+cμσ″α)。进行整理化简，得到衰减常数α和相移常数β。表达式为

(5)借鉴均匀半空间模型的趋肤深度推导思想，重新推导了基于分数阶电导率模型的广义趋肤深度公式，根据多孔极化介质的电导率正分数阶简化表达式，获得了地下极化介质的广义趋肤深度频域公式。

进一步地，根据均匀半空模型趋肤深度的定义方式，趋肤深度d等于场量振幅衰减到其表面值的e^-1是所经过的距离由此e^-αd＝e^-1，引入简化后的柯尔-柯尔模型σ＝σ′+σ″(iω)^c，根据波数表达式得到地下极化介质的广义趋肤深度公式为：其中，为趋肤深度，其中re(α)是代表取α的实部，α为衰减常数。

图3为采用图1所示的本发明一个实施例，计算含极化情况下，传统趋肤深度计算方法与新定义的极化介质广义趋肤深度计算方法对比结果，传统趋肤深度计算方法并不能体现极化率所带来的影响，因此无法对极化介质进行计算。

图4(a)(b)为采用图1所示的本发明一个实施例，计算基于极化介质广义趋肤深度与传统趋肤深度的相对误差和绝对误差，两者最大误差主要集中在低频部分，高频部分误差较小。因此，对于极化介质趋肤深度的计算，更适合在低频下进行。

以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。