CN113887106A - 一种基于Chikazumi模型的感应-磁化效应三维数值模拟方法 - Google Patents
一种基于Chikazumi模型的感应-磁化效应三维数值模拟方法 Download PDFInfo
- Publication number
- CN113887106A CN113887106A CN202111182187.8A CN202111182187A CN113887106A CN 113887106 A CN113887106 A CN 113887106A CN 202111182187 A CN202111182187 A CN 202111182187A CN 113887106 A CN113887106 A CN 113887106A
- Authority
- CN
- China
- Prior art keywords
- magnetic field
- induction
- magnetic
- control equation
- numerical simulation
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 230000000694 effects Effects 0.000 title claims abstract description 29
- 238000004088 simulation Methods 0.000 title claims abstract description 27
- 230000005291 magnetic effect Effects 0.000 claims abstract description 54
- 230000004044 response Effects 0.000 claims abstract description 22
- 239000011159 matrix material Substances 0.000 claims abstract description 20
- 230000035699 permeability Effects 0.000 claims abstract description 20
- 238000004364 calculation method Methods 0.000 claims abstract description 15
- 230000008569 process Effects 0.000 claims abstract description 8
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 7
- 238000001914 filtration Methods 0.000 claims abstract description 7
- 238000013507 mapping Methods 0.000 claims abstract description 7
- 239000006185 dispersion Substances 0.000 claims abstract description 4
- 230000005684 electric field Effects 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 3
- 230000014509 gene expression Effects 0.000 claims description 3
- 238000011160 research Methods 0.000 abstract description 2
- 230000007547 defect Effects 0.000 abstract 1
- 230000005415 magnetization Effects 0.000 description 7
- 238000001514 detection method Methods 0.000 description 5
- 229910052500 inorganic mineral Inorganic materials 0.000 description 3
- 239000011707 mineral Substances 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 230000001052 transient effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000006698 induction Effects 0.000 description 2
- 239000006249 magnetic particle Substances 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000004020 conductor Substances 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000005294 ferromagnetic effect Effects 0.000 description 1
- 229910052595 hematite Inorganic materials 0.000 description 1
- 239000011019 hematite Substances 0.000 description 1
- LIKBJVNGSGBSGK-UHFFFAOYSA-N iron(3+);oxygen(2-) Chemical compound [O-2].[O-2].[O-2].[Fe+3].[Fe+3] LIKBJVNGSGBSGK-UHFFFAOYSA-N 0.000 description 1
- SZVJSHCCFOBDDC-UHFFFAOYSA-N iron(II,III) oxide Inorganic materials O=[Fe]O[Fe]O[Fe]=O SZVJSHCCFOBDDC-UHFFFAOYSA-N 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Operations Research (AREA)
- Evolutionary Computation (AREA)
- Algebra (AREA)
- Geometry (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
本发明涉及一种基于Chikazumi模型的感应‑磁化效应三维数值模拟方法,针对磁化率模型在时域差分离散困难,无法直接进行时域三维数值模拟的问题,将Chikazumi复磁化率模型引入频域Maxwell方程组,提出了复磁导率的矩阵运算形式,优化控制方程的计算方式,根据差分离散后磁场在网格上的位置关系,重新建立了磁导率和磁场响应的科学映射矩阵,定义了磁异常体特别的边界设置方式,实现了利用频域有限差分方法对感应‑磁化效应的三维求解。并利用正余弦变换数值滤波算法,将频域结果变换到时域,最终实现了感应‑磁化效应时域三维数值模拟。本发明目的在于,可以克服目前研究方法仅能进行感应‑磁化效应的一维数值模拟,实现对感应‑磁化效应衰减过程的三维数值模拟。
Description
技术领域
本发明涉及一种地球物理勘探领域的基于Chikazumi模型的感应-磁化效应三维数值模拟方法。
背景技术
航空瞬变电磁法(Airbrone Electromagnetic Method,AEM)具有快速、高效、探测范围广等优点,已被广泛应用于矿产资源探勘查,环境监测,基础地质调查等多个领域(Lei,2006;Zhang,2005)。瞬变电磁法在探测过程中,发射电流的on-time阶段会让地下的磁性粒子改变方向,获得磁化;在电流关断后,磁性粒子会失去磁化强度,接收系统会观测到近似-1次幂率的衰减信号,这一过程被称为感应-磁化效应。国内外实验结果表明,感应-磁化效应不仅存在于多种的岩石及矿物中(Buselli,1982;Neumann,2006)。在0.1-15um半径范围内的磁铁矿和赤铁矿颗粒也可产生感应-磁化效应(Bournas,2017;Dabas andSkinner,2002)。随着机载发射机偶极矩的增加以及系统噪声水平的数量级的降低,电磁传播过程中的感应-磁化效应已可以被正常检测到(Kwan et al.,2015)。目前科学家们在多地均检测到感应-磁化效应的存在(Macnae,2017;Zhang,2020)。而传统的探测理论,往往会忽略感应-磁化效应的影响,引入较大的误差,无法对探测信息做出正确的解释。因此,对感应-磁化效应进行有效观测和准确解释尤为重要。
Chikazumi根据观测到的基本衰减规律,推导出了被广泛应用的Chikazumi磁化率模型(Chikazumi,1978)。Lee利用Chikazumi模型进行理论计算,验证了基于Chikazumi磁化率模型模拟出的电磁响应呈-1次幂率衰减,贴合实际观测到的情况(Lee,1984)。Kamnev等通过观测磁性岩石样品,确定了响应晚期呈-1.07~-1.13幂率衰减(Kamnev,2015)。Kozhevnikov等针对地面长导线源激发的地下均匀磁化介质与地面回线源激发的层状磁化介质的数值模拟进行了研究,分析了零频磁化率等性质对于感应-磁化效应的影响(Antonov,2017;Kozhevnikov,2008,2011,2018)。Cowan等推导了地面大回线源下感应-磁化效应垂直和径向分量的解析式,研究有关感应响应与磁化响应的交叉时间的估计方法(Cowan,2017)。Zhang等研究了长导线源下的感应-磁化效应,进行了层状模型下的一维以及二维数值模拟。并研究了岩、矿石样品的磁化率,验证了磁化效应对于铁磁性矿物探测的有效性(Zhang,2020)。然而,由于基于磁化率随时间变化的磁化率模型直接进行时域三维数值模拟存在困难,目前相关正演研究多数集中于一维数值模拟方面。
目前,常用的瞬变电磁三维数值模拟方法,主要有有限差分法(Wang,1993)、有限元法(Assous,1993)以及无网格法(Yang,2014)等。其中,有限差分方法是数值模拟最早采用的方法,与其他方法相比,虽然在计算复杂目标体会产生较大误差,但是较为直观,且具有较强的实用性。由于Chikazumi磁化率模型在时域的非线性,无法直接在时域进行运算。若采用频率域有限差分方法,由于一般的频率域有限差分方法使用真空磁导率,所以无需考虑磁导率的位置对应问题,在把磁导率变成复数模型直接进行运算时,需要将磁化率展开成矩阵进行除法运算,会导致交错网格位置对应错误以及计算困难,得到错误结果。因此,如何实现感应-磁化效应的三维数值模拟是待本领域技术人员迫切解决的一个技术问题。
发明内容
本发明所要解决的技术问题在于提供一种基于Chikazumi模型的感应-磁化效应三维数值模拟方法,针对磁化率模型在时域差分离散困难,无法直接进行时域三维数值模拟的问题,将Chikazumi复磁化率模型引入频域Maxwell方程组,磁场响含有复磁导率与复频率乘积的形式,传统频域有限差分方法无法求解复磁导率Maxwell方程组,提出了复磁导率的矩阵运算形式,优化控制方程的计算方式,根据差分离散后磁场在网格上的位置关系,重新建立了磁导率和磁场响应的科学映射矩阵,定义了磁异常体特别的边界设置方式,实现了利用频域有限差分方法对感应-磁化效应的三维求解。并利用正余弦变换数值滤波算法,将频域结果变换到时域,最终实现了感应-磁化效应时域三维数值模拟。本发明是这样实现的,一种基于Chikazumi模型的感应-磁化效应三维数值模拟方法:
包括如下步骤:
1)、采用频率域加源Maxwell旋度方程作为电、磁场分量Ex、Ey、Ez、Hx、Hy、Hz的控制方程;
2)、引入Chikazumi磁化率模型表征磁化介质磁导率变化过程,建立磁导率分量μx(ω)、μy(ω)、μz(ω)与磁场响应Hx、Hy、Hz的映射矩阵,并带入控制方程;
3)、采用非均匀三维Yee氏网格对计算区域进行剖分,基于频域有限差分法,推导控制方程在整个计算区域的差分格式,并将其改写为离散项;
4)、加载Dirichlet边界条件,将磁导率分量与磁场分量进行元素相乘后,保证磁场控制方程的优先运算,将其带入电场控制方程中;
5)、将控制方程整理成Ax=b的形式,通过计算x=A\b得到磁场各分量频率域响应结果;
6)、利用正余弦变换数值滤波算法,将磁场响应变换到时域,并对计算结果进行显示;
其中步骤2)中,Chikazumi磁化率模型表达式为:
式(1)中,ω表示角频率,i为虚部,χ0表示零频磁化率,τ1、τ2为时间常数τ的最小值和最大值;
对式(1)的复磁化率模型带入控制方程中,可以得出:
其中Ex、Ey、Ez为电场强度E在x、y、z三个方向上的分量,Hx、Hy、Hz为磁场强度H在x、y、z三个方向上的分量,Jx、Jy、Jz为导电电流密度J在x、y、z三个方向上的电流源密度,σx、σy、σz代表电导率σ在x,y,z三个方向上的分量;μx(ω)、μy(ω)、μz(ω)代表磁导率μ(ω)在x,y,z三个方向上的分量;
将差分矩阵进行离散,并加载Dirichlet边界条件,将磁场控制方程带入电场控制方程中,可以得到步骤4)中的运算矩阵为:
其中是电场强度E的差分矩阵,代表电场分量Ex、Ey、Ez在不同方向上的微分形式;是磁场H的差分矩阵,代表磁场分量Hx、Hy、Hz在不同方向上的微分形式;对于回线源J,只需要对Jx、Jy进行赋值,Jz分量总为零;将μx(ω)、μy(ω)、μz(ω)根据差分离散后的磁场分量建立映射矩阵,并带入控制方程进行运算;对于磁异常体,网格平面上的磁导率可以用两个相邻立方体的磁导率做面平均进行替代;
其中步骤5)中,将式(4)整理成Ax=b的形式,可得到:
其中,A11、A22、A33分别为:
将计算结果代入步骤6)中的正余弦变换数值滤波算法,可得到感应-磁化响应的时域解。
本发明与现有技术相比,有益效果在于:可以克服基于Chikazumi磁化率模型感应-磁化效应三维数值模拟困难的问题,实现基于频域有限差分方法的感应-磁化效应衰减过程的三维数值模拟。
附图说明
图1是一种基于Chikazumi模型的感应-磁化效应三维数值模拟方法示意图;
图2是磁化半空间模型下三维数值模拟模拟结果与解析解响应及误差对比;
图3复杂磁化模型下,三维数值模拟感应电动势水平切片图;
图4复杂磁化模型下,三维数值模拟感应电动势垂直切片图;
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
实施例
参见图1,一种基于Chikazumi模型的感应-磁化效应三维数值模拟方法,包括:
1)、采用频域Maxwell旋度方程作为磁场分量的控制方程,准静态条件下、结合非均匀Yee式网格和Dirichlet边界条件,构建网格差分格式,并且网格数目为38×38×32,其中x、y方向上的网格数目均为38个,z方向上网格数32个,相邻网格的比例系数为2,最小和最大网格步长分别为3m和2560m。
2)在整个计算区域内设置电导率、零频磁化率等参数,算例中背景电导率设置为0.01西门子/米,分别在距离地面10m和40m的深度设置一小一大两个磁异常体,分别位于发射线圈的两侧,异常体1边长为12m×12m×10m,异常体2边长为36m×36m×40m,设置异常体零频磁化率χ0=0.01,τ1=10-6,τ1=106。在计算过程中,网格边缘的电导率可以用相邻的四个立方体的电导率做体平均进行替代,网格平面上的磁导率可以用两个相邻立方体的磁导率做面平均进行替代。
3)、根据磁场响应差分离散后网格上的位置关系,建立磁导率和磁场响应的科学映射矩阵。
4)、在10-1Hz-106Hz范围内等间隔选取70个频点进行计算,利用控制方程计算磁感应强度B的频域解。
5)、判断所有频点是否计算完成,若未完成,重复步骤3),若完成则输出响应矩阵。
6)、选取时间范围为10-5-10-2s,将响应矩阵带入正余弦变换滤波算法中,逐点进行频时变换。
7)、判断是否完成全部频点的频时变换,若未完成,重复步骤6),若完成则输出时域响应结果。
8)、根据所选时刻的磁场值,对响应的感应电动势进行结果显示。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种基于Chikazumi模型的感应-磁化效应三维数值模拟方法其特征在于,包括如下步骤:
1)、采用频率域加源Maxwell旋度方程作为电、磁场分量Ex、Ey、Ez、Hx、Hy、Hz的控制方程;
2)、引入Chikazumi磁化率模型表征磁化介质磁导率变化过程,建立磁导率分量μx(ω)、μy(ω)、μz(ω)与磁场响应Hx、Hy、Hz的映射矩阵,并带入控制方程;
3)、采用非均匀三维Yee氏网格对计算区域进行剖分,基于频域有限差分法,推导控制方程在整个计算区域的差分格式,并将其改写为离散项;
4)、加载Dirichlet边界条件,将磁导率分量与磁场分量进行元素相乘后,保证磁场控制方程的优先运算,将其带入电场控制方程中;
5)、将控制方程整理成Ax=b的形式,通过计算x=A\b得到磁场各分量频率域响应结果;
6)、利用正余弦变换数值滤波算法,将磁场响应变换到时域,并对计算结果进行显示;
其中步骤2)中,Chikazumi磁化率模型表达式为:
式(1)中,ω表示角频率,i为虚部,χ0表示零频磁化率,τ1、τ2为时间常数τ的最小值和最大值;
对式(1)的复磁化率模型带入控制方程中,可以得出:
其中Ex、Ey、Ez为电场强度E在x、y、z三个方向上的分量,Hx、Hy、Hz为磁场强度H在x、y、z三个方向上的分量,Jx、Jy、Jz为导电电流密度J在x、y、z三个方向上的电流源密度,σx、σy、σz代表电导率σ在x,y,z三个方向上的分量;μx(ω)、μy(ω)、μz(ω)代表磁导率μ(ω)在x,y,z三个方向上的分量;
将差分矩阵进行离散,并加载Dirichlet边界条件,将磁场控制方程带入电场控制方程中,可以得到步骤4)中的运算矩阵为:
其中是电场强度E的差分矩阵,代表电场分量Ex、Ey、Ez在不同方向上的微分形式;是磁场H的差分矩阵,代表磁场分量Hx、Hy、Hz在不同方向上的微分形式;对于回线源J,只需要对Jx、Jy进行赋值,Jz分量总为零;将μx(ω)、μy(ω)、μz(ω)根据差分离散后的磁场分量建立映射矩阵,并带入控制方程进行运算;对于磁异常体,网格平面上的磁导率可以用两个相邻立方体的磁导率做面平均进行替代;
其中步骤5)中,将式(4)整理成Ax=b的形式,可得到:
其中,A11、A22、A33分别为:
将计算结果代入步骤6)中的正余弦变换数值滤波算法,可得到感应-磁化响应的时域解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111182187.8A CN113887106B (zh) | 2021-10-11 | 2021-10-11 | 一种基于Chikazumi模型的感应-磁化效应三维数值模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111182187.8A CN113887106B (zh) | 2021-10-11 | 2021-10-11 | 一种基于Chikazumi模型的感应-磁化效应三维数值模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113887106A true CN113887106A (zh) | 2022-01-04 |
CN113887106B CN113887106B (zh) | 2024-04-12 |
Family
ID=79006000
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111182187.8A Active CN113887106B (zh) | 2021-10-11 | 2021-10-11 | 一种基于Chikazumi模型的感应-磁化效应三维数值模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113887106B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114970289A (zh) * | 2022-07-25 | 2022-08-30 | 中南大学 | 三维大地电磁各向异性正演数值模拟方法、设备及介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070061123A1 (en) * | 2005-09-09 | 2007-03-15 | Fujitsu Limited | Electromagnetic field simulator and electromagnetic field simulation program storage medium |
CN105893678A (zh) * | 2016-04-01 | 2016-08-24 | 吉林大学 | 一种时域有限差分的三维感应-极化双场数值模拟方法 |
CN112285788A (zh) * | 2020-10-29 | 2021-01-29 | 吉林大学 | 一种基于电磁波动方程的cpml吸收边界条件加载方法 |
WO2021068527A1 (zh) * | 2019-10-12 | 2021-04-15 | 中南大学 | 一种足迹引导的高效航空电磁法数值模拟方法 |
-
2021
- 2021-10-11 CN CN202111182187.8A patent/CN113887106B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070061123A1 (en) * | 2005-09-09 | 2007-03-15 | Fujitsu Limited | Electromagnetic field simulator and electromagnetic field simulation program storage medium |
CN105893678A (zh) * | 2016-04-01 | 2016-08-24 | 吉林大学 | 一种时域有限差分的三维感应-极化双场数值模拟方法 |
WO2021068527A1 (zh) * | 2019-10-12 | 2021-04-15 | 中南大学 | 一种足迹引导的高效航空电磁法数值模拟方法 |
CN112285788A (zh) * | 2020-10-29 | 2021-01-29 | 吉林大学 | 一种基于电磁波动方程的cpml吸收边界条件加载方法 |
Non-Patent Citations (3)
Title |
---|
严波;李予国;韩波;梅晓勇;贺少波;许思维;: "任意方位电偶源的MCSEM电磁场三维正演", 石油地球物理勘探, no. 04, 15 August 2017 (2017-08-15) * |
关珊珊;殷长春;嵇艳鞠;黎东升;耿毅男;: "三维随机断裂带的航空时域电磁响应数值模拟", 地球物理学报, no. 11, 15 November 2018 (2018-11-15) * |
马炳镇;: "起伏地形下地面瞬变电磁法三维正演数值模拟研究", 物探与化探, no. 04, 15 August 2018 (2018-08-15) * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114970289A (zh) * | 2022-07-25 | 2022-08-30 | 中南大学 | 三维大地电磁各向异性正演数值模拟方法、设备及介质 |
CN114970289B (zh) * | 2022-07-25 | 2022-10-25 | 中南大学 | 三维大地电磁各向异性正演数值模拟方法、设备及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN113887106B (zh) | 2024-04-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101194184A (zh) | 利用多个频率空间上解释电磁数据的方法 | |
Özyıldırım et al. | Two-dimensional inversion of magnetotelluric/radiomagnetotelluric data by using unstructured mesh | |
CN113887106B (zh) | 一种基于Chikazumi模型的感应-磁化效应三维数值模拟方法 | |
Rahmani et al. | Sensitivity of dipole magnetic tomography to magnetic nanoparticle injectates | |
Heagy et al. | Investigating the potential of using conductive or permeable proppant particles for hydraulic fracture characterization | |
Zlobinskiy et al. | Applying TM-polarization geoelectric exploration for study of low-contrast three-dimensional targets | |
Grigorev et al. | On the applicability of electromagnetic monitoring of hydraulic fracturing | |
Smith et al. | The moments of the impulse response: A new paradigm for the interpretation of transient electromagnetic data | |
CN110119586B (zh) | 轴向电导率各向异性瞬变电磁三分量三维fdtd正演方法 | |
Huang et al. | Detecting metal objects in magnetic environments using a broadband electromagnetic method | |
Ji et al. | Calculation and application of full-wave airborne transient electromagnetic data in electromagnetic detection | |
Lewis et al. | The effect of host rock on transient electromagnetic fields | |
Ruo et al. | The analysis of CSAMT responses of dyke embedded below conductive overburden | |
CN104612661B (zh) | 一种随钻电磁波测井装置和方法 | |
Smith | Induced-polarization effects in airborne electromagnetic data: Estimating chargeability from shape reversals | |
Li et al. | Land-based TEM data processing: from turn-off ramp to full waveform | |
Dang et al. | Long-Distance Crosswell EM Logging of Copper Ore Using Borehole-Surface Current Injection in Slim Holes | |
Zhang et al. | Three-dimensional transient electromagnetic response of multiple grounded-wire sources | |
Zhao et al. | Transient electromagnetic detection of unexploded ordnance buried in underwater sediments | |
Marchant et al. | Modelling IP effects in airborne time domain electromagnetics | |
CN115935746A (zh) | 一种时域感应-磁化-极化效应三维数值模拟方法 | |
Yi et al. | Three-Dimensional Inversion of ATEM Responses for Grounded Source | |
CN113779853B (zh) | 一种时域电磁感应-磁化效应分数阶三维数值模拟方法 | |
Shubitidze et al. | The Crank-Nicolson FDTD Method in Cylindrical Coordinates and Its Application to Underwater UXO Detection and Classification Problems | |
Liu* et al. | 3D forward modeling and analysis for loop-source transient electromagnetic data based on finite-volume method in an arbitrarily anisotropy earth |
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 |