CN110135022B - 一种基于极化介质模型的广义趋肤深度计算方法 - Google Patents
一种基于极化介质模型的广义趋肤深度计算方法 Download PDFInfo
- Publication number
- CN110135022B CN110135022B CN201910348025.3A CN201910348025A CN110135022B CN 110135022 B CN110135022 B CN 110135022B CN 201910348025 A CN201910348025 A CN 201910348025A CN 110135022 B CN110135022 B CN 110135022B
- Authority
- CN
- China
- Prior art keywords
- medium
- skin depth
- conductivity
- underground
- expression
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Remote Sensing (AREA)
- Geology (AREA)
- Geophysics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种基于地下极化介质的频率域广义趋肤深度计算方法,目的在于提高地下极化介质的探测深度和解释精度。本发明主要基于地下极化介质的柯尔‑柯尔分数阶电导率模型,推导频率域广义趋肤深度计算公式,建立了极化率与广义趋肤深度的关系。在时域电磁法测量中,先对实测数据进行噪声滤波、取样叠加等处理,采用广义S变换将时域数据变换到频率域,利用反演法提取地下介质的电导率、极化率、充电率和频散系数,利用波数定义式求解衰减和相移常数,最后在频率域计算地下极化介质的广义趋肤深度。与经典趋肤深度的计算方法相比,本发明提出的极化介质广义趋肤深度,更符合电磁波在实际大地的传播规律,可以更精确的计算可极化金属矿的深度。
Description
技术领域
本发明涉及地球物理勘探领域,尤其涉及一种基于极化介质模型的广义趋肤深度计算方法。
背景技术
趋肤深度是电磁勘探领域中一个非常重要的概念,通常将其用来估算电磁法的勘探深度。在实际电法勘探中,为了解释实测数据,需要计算视电导率与探测深度数值,并表示成视电导率-深度断面的形式,因此,需要计算在实测不同频率或不同时刻情况下的探测深度。
目前有关于电磁探测深度的研究,Szarka、Fischer(1989)提出,采用1.5倍趋肤深度来估计大地电磁的理论探测深度。Yan S(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中,地下极化介质的广义趋肤深度为:
进一步地,
其中,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ω的负分数阶项,进行有理化整理,获得多孔极化介质的电导率正分数阶简化表达式。
σ=(1-η)(iωτ)c
分子分母同时乘以,
则:σ=σ′+σ″(iω)c。
(4)根据波数表达式求解虚部和实部表达式,进行简化整理得出衰减常数及相移常数表达式。
进一步地,根据波数表达式k2=ω2με-iωμσ(ω),其中ω为角频率,σ为电导率,ε为介电常数。简化后的极化介质电导率σ=σ′+σ″(iω)c代入波数表达式中整理虚部和实部,ic=a+ib,k=α+iβ,写为:α2-β2=ω2με+ω1+cμσ″b,2αβ=-(ωμσ′+ω1+cμσ″α)。进行整理化简,得到衰减常数α和相移常数β。表达式为
(5)借鉴均匀半空间模型的趋肤深度推导思想,重新推导了基于分数阶电导率模型的广义趋肤深度公式,根据多孔极化介质的电导率正分数阶简化表达式,获得了地下极化介质的广义趋肤深度频域公式。
进一步地,根据均匀半空模型趋肤深度的定义方式,趋肤深度d等于场量振幅衰减到其表面值的e-1是所经过的距离由此e-αd=e-1,引入简化后的柯尔-柯尔模型σ=σ′+σ″(iω)c,根据波数表达式得到地下极化介质的广义趋肤深度公式为:其中,为趋肤深度,其中Re(α)是代表取α的实部,α为衰减常数。
图3为采用图1所示的本发明一个实施例,计算含极化情况下,传统趋肤深度计算方法与新定义的极化介质广义趋肤深度计算方法对比结果,传统趋肤深度计算方法并不能体现极化率所带来的影响,因此无法对极化介质进行计算。
图4(a)(b)为采用图1所示的本发明一个实施例,计算基于极化介质广义趋肤深度与传统趋肤深度的相对误差和绝对误差,两者最大误差主要集中在低频部分,高频部分误差较小。因此,对于极化介质趋肤深度的计算,更适合在低频下进行。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种基于极化介质模型的广义趋肤深度计算方法,其特征在于,包括如下步骤:
(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μσ″α),进行整理化简,得到衰减常数α和相移常数β,表达式为:
3.按照权利要求1所述的一种基于极化介质模型的广义趋肤深度计算方法,其特征在于:
采用反演法从频域中提取地下介质电导率、极化率、充电率、频散系数包括:求解目标函数最小值,目标函数表示为:
其中,M代表的是原始模型,n为实测的频点个数;pαi为第i个频点的实测响应值,其为视电阻率、相位或者电场与磁场值,Mj,j=1,2,3,…m为模型参数;pci是由预测模型正演计算所得到的第i个频率的响应参数,α=2为范数;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910348025.3A CN110135022B (zh) | 2019-04-28 | 2019-04-28 | 一种基于极化介质模型的广义趋肤深度计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910348025.3A CN110135022B (zh) | 2019-04-28 | 2019-04-28 | 一种基于极化介质模型的广义趋肤深度计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110135022A CN110135022A (zh) | 2019-08-16 |
CN110135022B true CN110135022B (zh) | 2022-10-14 |
Family
ID=67575386
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910348025.3A Active CN110135022B (zh) | 2019-04-28 | 2019-04-28 | 一种基于极化介质模型的广义趋肤深度计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110135022B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114487687B (zh) * | 2022-04-15 | 2022-07-01 | 华北电力大学 | 一种非接触式快速测量金属板透入深度的装置及方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107657137A (zh) * | 2017-11-09 | 2018-02-02 | 吉林大学 | 一种有理函数逼近的分数阶电磁反常扩散三维模拟方法 |
CN108169802A (zh) * | 2018-03-02 | 2018-06-15 | 吉林大学 | 一种粗糙介质模型的时域电磁数据慢扩散成像方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9857316B2 (en) * | 2012-01-13 | 2018-01-02 | University Of Notre Dame Du Lac | Methods and apparatus for electromagnetic signal polarimetry sensing |
-
2019
- 2019-04-28 CN CN201910348025.3A patent/CN110135022B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107657137A (zh) * | 2017-11-09 | 2018-02-02 | 吉林大学 | 一种有理函数逼近的分数阶电磁反常扩散三维模拟方法 |
CN108169802A (zh) * | 2018-03-02 | 2018-06-15 | 吉林大学 | 一种粗糙介质模型的时域电磁数据慢扩散成像方法 |
Non-Patent Citations (5)
Title |
---|
Steven L. Jacques等.Imaging skin pathology with polarized light.《Journal of Biomedical Optics 》.2002, * |
YANJU JI等.Three-Dimensional Transient Electromagnetic Modeling Based on Fictitious Wave Domain Methods.《Pure and Applied Geophysics》.2017, * |
嵇艳鞠等.基于神经网络电性源半航空视电阻率反演研究.《电波科学学报》.2014,(第05期), * |
王志刚等.时频电磁技术的新进展及应用效果.《石油地球物理勘探》.2016, * |
符超等.基于Cole-Cole模型的中间极化水平层大地电磁IP效应研究.《地球物理学进展》.2016,(第02期), * |
Also Published As
Publication number | Publication date |
---|---|
CN110135022A (zh) | 2019-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang | Multichannel matching pursuit for seismic trace decomposition | |
Wu et al. | Removal of multisource noise in airborne electromagnetic data based on deep learning | |
US20050075789A1 (en) | Method and apparatus for inversion processing of well logging data in a selected pattern space | |
CN110133733B (zh) | 一种基于粒子群优化算法的电导-极化率多参数成像方法 | |
CN102832908A (zh) | 基于小波变换与变步长lms自适应滤波的信号降噪方法 | |
CN108680966B (zh) | 海洋可控源电磁勘探噪声降噪效果评估方法 | |
Won | A wide-band electromagnetic exploration method—Some theoretical and experimental results | |
GB2581087A (en) | Determination of dispersive resistive and dielectric properties of earth formations using multi-frequency measurements | |
US20110166842A1 (en) | Layer stripping method | |
CN111273367B (zh) | 一种大地电磁阻抗的估算方法 | |
CN110135022B (zh) | 一种基于极化介质模型的广义趋肤深度计算方法 | |
Yue et al. | Suppression of periodic interference during tunnel seismic predictions via the Hankel-SVD-ICA method | |
CN109188542B (zh) | 一种波区相关性检测的远参考大地电磁阻抗计算方法 | |
Jegen et al. | Using homotopy to invert geophysical data | |
CN108376204B (zh) | 一种基于地下粗糙介质模型的电磁广义趋肤深度计算方法 | |
Moilanen | Modern methods of airborne electromagnetic survey | |
Beiki et al. | Leveling HEM and aeromagnetic data using differential polynomial fitting | |
CN116879964A (zh) | 一种时频电磁频率域数据自约束稳健电阻率反演方法 | |
Qi et al. | Full waveform modeling of transient electromagnetic response based on temporal interpolation and convolution method | |
Lu et al. | Robust estimation of instantaneous phase using a time-frequency adaptive filter | |
Zhu et al. | Suppression of motion noise based on a linear‐homomorphic filtering algorithm in airborne electromagnetic survey | |
Zhao et al. | Fault Diagnosis of Oil Pumping Machine Retarder Based on Sound Texture‐Vibration Entropy Characteristics and Gray Wolf Optimization‐Support Vector Machine | |
CN115016008A (zh) | 基于神经网络的电性源感应-极化共生效应多参数成像方法 | |
CN116430464A (zh) | 提高时频电磁法频率域资料处理质量的方法 | |
Liu et al. | 3-D Full-Waveform Modeling and Analysis of Induced Polarization and Magnetic Viscosity Effect in Time-Domain Electromagnetic Method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |