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

本发明专利技术涉及一种基于地下极化介质的频率域广义趋肤深度计算方法，目的在于提高地下极化介质的探测深度和解释精度。本发明专利技术主要基于地下极化介质的柯尔‑柯尔分数阶电导率模型，推导频率域广义趋肤深度计算公式，建立了极化率与广义趋肤深度的关系。在时域电磁法测量中，先对实测数据进行噪声滤波、取样叠加等处理，采用广义S变换将时域数据变换到频率域，利用反演法提取地下介质的电导率、极化率、充电率和频散系数，利用波数定义式求解衰减和相移常数，最后在频率域计算地下极化介质的广义趋肤深度。与经典趋肤深度的计算方法相比，本发明专利技术提出的极化介质广义趋肤深度，更符合电磁波在实际大地的传播规律，可以更精确的计算可极化金属矿的深度。

A Generalized Skin Depth Calculation Method Based on Polarized Medium Model

【技术实现步骤摘要】
一种基于极化介质模型的广义趋肤深度计算方法
本专利技术涉及地球物理勘探领域，尤其涉及一种基于极化介质模型的广义趋肤深度计算方法。
技术介绍
趋肤深度是电磁勘探领域中一个非常重要的概念，通常将其用来估算电磁法的勘探深度。在实际电法勘探中，为了解释实测数据，需要计算视电导率与探测深度数值，并表示成视电导率-深度断面的形式，因此，需要计算在实测不同频率或不同时刻情况下的探测深度。目前有关于电磁探测深度的研究，Szarka、Fischer(1989)提出，采用1.5倍趋肤深度来估计大地电磁的理论探测深度。YanS(2009)采用解析法、有限差分时域法和时频法对电磁测深探测深度进行了研究。目的是为现场工作设计提供依据，以达到预期目标。薛国强(2014)以接地长导线源为例研究了可控源电磁场全域的有效趋肤深度。CN108376204A公开了一种基于地下粗糙介质模型的电磁广义趋肤深度计算方法，在趋肤深度定义方式的基础上，根据粗糙介质广义电导率的定义,推导出适用于符合地下粗糙介质电导率的电磁广义趋肤深度公式，根据不同岩石层粗糙度参数，计算出相应的探测深度。但以上探测深度计算方法是基于均匀介质以及粗糙介质理论下的计算公式，而实际地下介质体结构极为复杂，尤其地下多孔极化介质或者多金属矿探测时，传统的趋肤深度定义不再适用，因此，重新定义基于分数阶电导率模型的广义趋肤深度具有重要意义。
技术实现思路
本专利技术的目的在于针对传统均匀半空间模型的趋肤深度计算方法，已经不满足地下极化介质的频散特性等传播规律，根据极化介质的分数阶电导率，以及趋肤深度的定义方式，提供一种基于地下极化介质模型的电磁广义趋肤深度计算方法。本专利技术在传统趋肤深度定义方式的基础上，根据极化介质电导率的定义，推导出适用于符合地下极化介质电导率的电磁广义趋肤深度公式，根据不同地下极化介质的电导率、极化率、充电率、频散系数，计算出探测深度。本专利技术是这样实现的，一种基于地下极化介质的广义趋肤深度频域计算方法，包括如下步骤：(1)先对实测数据进行噪声滤波、取样叠加处理，采用广义S变换将数据进行时频变换，获得频域形式的数据；(2)采用反演法对步骤1中处理后的电磁数据，进行地下介质电导率、极化率、充电率、频散系数的提取；(3)将柯尔-柯尔分数阶电导率模型iω的负分数阶项，进行有理化整理，获得多孔极化介质的电导率正分数阶简化表达式；(4)电导率正分数阶简化表达式代入到波数表达式求解虚部和实部表达式，进行简化整理得出衰减常数及相移常数表达式；(5)对衰减常数取实部求取地下极化介质的广义趋肤深度。进一步地，步骤3中，根据柯尔-柯尔模型，地下极化介质分数阶电导率表示为将地下极化介质分数阶电导率表达式中的iω的负分数阶项，进行有理化整理，转化为正分数阶的简化形式，表示为：σ＝σ′+σ″(iω)c，其中σ∞为高频交流电导率，τ为充电时间常数，η为极化率，c为频散系数，ω为角频率。进一步地，步骤4中，根据波数表达式k2＝ω2με-iωμσ(ω)，其中ω为角频率，σ为电导率，ε为介电常数，步骤3中简化后的极化介质电导率代入波数表达式中整理虚部和实部，令ic＝a+ib，k＝α+iβ，写为：α2-β2＝ω2με+ω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.一种基于地下极化介质的广义趋肤深度频域计算方法，其特征在于，包括如下步骤：(1)先对实测数据进行噪声滤波、取样叠加处理，采用广义S变换将数据进行时频变换，获得频域形式的数据；(2)采用反演法对步骤1中处理后的电磁数据，进行地下介质电导率、极化率、充电率、频散系数的提取；(3)将柯尔‑柯尔分数阶电导率模型iω的负分数阶项，进行有理化整理，获得多孔极化介质的电导率正分数阶简化表达式；(4)电导率正分数阶简化表达式代入到波数表达式求解虚部和实部表达式，进行简化整理得出衰减常数及相移常数表达式；(5)对衰减常数取实部求取地下极化介质的广义趋肤深度。

【技术特征摘要】
1.一种基于地下极化介质的广义趋肤深度频域计算方法，其特征在于，包括如下步骤：(1)先对实测数据进行噪声滤波、取样叠加处理，采用广义S变换将数据进行时频变换，获得频域形式的数据；(2)采用反演法对步骤1中处理后的电磁数据，进行地下介质电导率、极化率、充电率、频散系数的提取；(3)将柯尔-柯尔分数阶电导率模型iω的负分数阶项，进行有理化整理，获得多孔极化介质的电导率正分数阶简化表达式；(4)电导率正分数阶简化表达式代入到波数表达式求解虚部和实部表达式，进行简化整理得出衰减常数及相移常数表达式；(5)对衰减常数取实部求取地下极化介质的广义趋肤深度。2.按照权利要求1所述的一种基于地下极化介质的广义趋肤深度频域计算方法，其特征在于：步骤3中，根据柯尔-柯尔模型，地下极化介质分数阶电导率表示为将地下极化介质分数阶电导率表达式中的iω的负分数阶项，进行有理化整理，转化为正分数阶的简化形式，表示为：σ＝σ′+σ″(iω)c，其中σ∞为高频交流电导率，τ为充电时间常数，η为极化率，c为频散系数，ω为角频率。3.按照权利要求2所述的一种基于地下极化介质的广义趋肤深度频域计算方法，其特征在于：步骤4中，根据波数表达式k2＝ω2με-iωμσ(ω)，其中ω为角频率，σ为电导...

【专利技术属性】
技术研发人员：嵇艳鞠，孟祥东，吴琼，吴燕琪，黎东升，栾卉，
申请(专利权)人：吉林大学，
类型：发明
国别省市：吉林,22