CN106019394B - 海洋大地电磁场非线性共轭梯度三维并行反演方法 - Google Patents
海洋大地电磁场非线性共轭梯度三维并行反演方法 Download PDFInfo
- Publication number
- CN106019394B CN106019394B CN201610268705.0A CN201610268705A CN106019394B CN 106019394 B CN106019394 B CN 106019394B CN 201610268705 A CN201610268705 A CN 201610268705A CN 106019394 B CN106019394 B CN 106019394B
- Authority
- CN
- China
- Prior art keywords
- seafloor
- model
- sea
- data
- impedance
- 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
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
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种海洋大地电磁场非线性共轭梯度三维并行反演方法,将固定电阻率值的海水层加入到初始模型,推导带有海水层模型各节点的大地电磁场场值的插值表达式和计算过程,以实现正演计算;推导带有海水层模型各节点的模型修改量(计算过程中主要包括目标函数、目标函数梯度、查找方向和查找步长等变量参数)的计算表达式和求解过程,以实现反演计算。本发明适用于海洋或其他包含水域的大地电磁场非线性共轭梯度三维反演。
Description
技术领域
本发明涉及大地电磁测深方法的处理技术领域,具体地说是一种海洋大地电磁场非线性共轭梯度三维并行反演方法。
背景技术
大地电磁测深方法广泛应用于深部地质探测和矿产、水、石油、地热等资源勘查。而目前二维反演是其数据处理的主流方法,三维反演方法的研究较少,国外科学家分别提出了快速松弛法、共轭梯度法、非线性共轭梯度法等方法,国内已经实现了快速松弛法和共轭梯度法。
目前,对海洋资源的开发已经成为重要的发展方向。现有的大地电磁三维反演均为陆地资料处理方法,只能反演地表(空气与地面界面)测量的大地电磁场数据,不能处理海底等水底大地测量的数据。
发明内容
有鉴于此,本发明实施例提供一种海洋大地电磁场非线性共轭梯度三维并行反演方法,主要目的是提供一种适用于包含水域数据的大地电磁场反演的方法。
为达到上述目的,本发明主要提供如下技术方案:
一方面,本发明实施例提供了一种海洋大地电磁场非线性共轭梯度三维反演方法,包括如下步骤:
a.在笛卡尔坐标系下沿x、y、z三个坐标轴将模型空间划分成Nx、Ny、Nz个小的长方体网格单元,间距为Δx(i)(i=1,...,Nx)、Δy(j)(j=1,...,Ny)、Δz(k)(k=1,...,Nz),其中模型空间包括空气层、海水层和地下半空间;海底节点(i,j,k)的长方体网格单元的长度、宽度和高度分别为xseafloor-i、yseafloor-j、zseafloor-k,电阻率为ρseafloor-i,j,k,其中电场取在各网格单元边缘的中点,用和表示;磁场取在各网格单元表面的中心,用和表示;
b.通过麦克斯韦方程组和边界条件的离散化得到模型各节点的场值表达式及其计算,实现正演;
c.将正演得到的阻抗与实测阻抗对比,获得阻抗数据偏差;
d.通过阻抗数据偏差计算目标函数,判断目标函数值是否足够小,当目标函数值足够小时结束迭代,否则进入下一步;
e.计算目标函数的梯度,以非线性共轭梯度法作为反演方法获得查找步长和查找方向,从而计算出模型各节点的模型修改量,其中使用阻抗数据误差和当前迭代计算所得的电阻率计算反演方法中的预处理因子;
f.使用模型修改量修正模型,并返回步骤b。
作为优选,其中根据计算机的CPU数量自动将每个频点的数据分配给所有的CPU进行正演以及目标函数梯度的并行计算。
作为优选,节点的场值表达式,Eseafloor-x如下所示:
其中seafloor表示海底及地层网格节点,sea表示海水层网格节点。
作为优选,使用电场值E,通过麦克斯韦方程能够求出各网格节点的磁场值H,便能得到模型的阻抗响应Z,
作为优选,反演目标函数为:
其中,和Zseafloor-n专指海底观测数据和正演数据,εseafloor-n为观测数据误差,第1~N个数据为阻抗各分量的实部,第N+1~2N个数据为阻抗各分量的虚部,W=Wair+Wsea+Wseafloor为预调节矩阵,m=mair+msea+mseafloor为网格剖分模型,air表示空气层部分,sea表示海水层部分,seafloor表示海底及地下半空间部分,
目标函数梯度表示为:
三维阻抗灵敏度张量表示为:
其中,K=Kair+Ksea+Kseafloor为正演系数矩阵,E=Eair+Esea+Eseafloor为场值并分为两种源条件下的场值,
为正演模型有限差分网格至接收点插入两种极化源电磁场向量的线性组合;
得到:
p(i)=-C(i) -1g(i)+β(i)p(i-1) (14)
其中γ为模型电阻率和数据误差相关的变量,从而得到更新的模型:
mseafloor-(i+1)=mseafloor-(i)+α(i)p(i) (18)。
与现有技术相比,本发明的有益效果在于:
本发明的海洋大地电磁场非线性共轭梯度三维反演方法能够将野外采集到的海底(或湖、河等水底)大地电磁场测深数据通过反演计算转换成地下三维电性结构信息,反映更为真实的地下电性结构信息,根据海底测量的数据反映出海底测区地壳浅表、地壳浅层、地壳深层、壳幔边界、上地幔的电阻率信息。该反演方法在保证计算精度的条件下,计算效率很高,适用于目前国内外十分主流的Windows操作系统和Linux、Unix操作系统,并且可以在普通个人电脑上推广使用。该新技术能够为探索海洋成矿区带不同深度尺度上地球物理特征与成矿之间的关系、形成海洋矿产资源立体探测的技术解决方案和深部资源勘查提供依据技术支持,指明海洋找矿方向,重塑深部成矿动力学过程,为建立海洋成矿理论提供更可靠的信息。
附图说明
图1为本发明的海洋大地电磁场非线性共轭梯度三维反演方法的流程框图;
图2为本发明的三维模型的示意图;
图3为图2中的单个长方体网格单元示意图;
图4为本发明方法的应用例1的模型示意图;
图5为本发明方法的应用例2的模型示意图;
图6为本发明方法的应用例2的三维反演模型的二维切片图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步详细描述,但不作为对本发明的限定。
本发明中的海洋大地电磁场是指与陆地大地电磁场相对应的,包括水域数据的大地电磁场,如海底、湖底、河底等的大地电磁场。本发明实施例中未尽之处,均可从现有技术中获得。
图1为本发明的海洋大地电磁场非线性共轭梯度三维反演方法的流程框图。如图1所示,海洋大地电磁场非线性共轭梯度三维反演方法,包括如下步骤:
a.在笛卡尔坐标系下沿x、y、z三个坐标轴将模型空间划分成Nx、Ny、Nz个小的长方体网格单元,间距为Δx(i)(i=1,...,Nx)、Δy(j)(j=1,...,Ny)、Δz(k)(k=1,...,Nz),其中模型空间包括空气层、海水层和地下半空间;海底节点(i,j,k)的长方体网格单元的长度、宽度和高度分别为xseafloor-i、yseafloor-j、zseafloor-k,电阻率为ρseafloor-i,j,k,其中电场取在各网格单元边缘的中点,用和表示;磁场取在各网格单元表面的中心,用和表示;
b.通过麦克斯韦方程组和边界条件的离散化得到模型各节点的场值表达式及其计算,实现正演;
c.将正演得到的阻抗与实测阻抗对比,获得阻抗数据偏差;
d.通过阻抗数据偏差计算目标函数,判断目标函数值是否足够小,当目标函数值足够小时结束迭代,否则进入下一步;
e.计算目标函数的梯度,以非线性共轭梯度法作为反演方法获得查找步长和查找方向,从而计算出模型各节点的模型修改量,其中使用阻抗数据误差和当前迭代计算所得的电阻率计算反演方法中的预处理因子;
f.使用模型修改量修正模型,并返回步骤b。
本发明方法将固定电阻率值(海水的电阻率值是根据实际情况设计的,在正、反演过程中保持不变)的海水层加入到初始模型,推导带有海水层模型各节点的场值插值表达式和计算过程,以实现正演计算;推导带有海水层模型各节点的模型修改量(计算过程中主要包括目标函数、目标函数梯度、查找方向和查找步长等变量参数)的计算表达式和求解过程,以实现反演计算;添加基于频点的并行计算结构,以提高效率,最终达到快速获得准确地下电性结构分布的目的。本发明方法能够将野外采集到的海底(或湖、河等水底)大地电磁测深数据通过反演计算转换成地下三维电性结构信息,反映更为真实的地下深部电性结构信息,根据测量数据能够反映出海底地壳浅表、地壳浅层、地壳深层、壳幔边界、上地幔的电阻率信息。能够为探索海洋成矿区带不同深度尺度上地球物理特征与成矿之间的关系、形成海洋矿产资源立体探测的技术解决方案和海底深部资源勘查提供依据技术支持,指明找矿方向,查明成矿系统的结构,重塑深部成矿动力学过程,为建立海洋地区成矿理论提供更可靠的信息。
本发明实施例中建立的三维模型参见图2和图3。在笛卡尔坐标系下沿x、y、z三个坐标轴将模型空间划分成Nx、Ny、Nz个小的长方体网格单元,间距为Δx(i)(i=1,...,Nx)、Δy(j)(j=1,...,Ny)、Δz(k)(k=1,...,Nz),其中模型空间包括空气层、海水层和地下半空间;海底节点(i,j,k)的长方体网格单元的长度、宽度和高度分别为xseafloor-i、yseafloor-j、zseafloor-k,电阻率为ρseafloor-i,j,k,其中电场取在各网格单元边缘的中点,用和表示;磁场取在各网格单元表面的中心,用和表示。图3为图2中的单个长方体网格单元示意图,图3所示的是编号为(i,j,k)的长方体网格单元。
海底大地电磁场三维正演:
设随时间变化的谐波函数为e-iωt,麦克斯韦方程式组用微分形式表示为:
电场强度E满足矢量方程:
其中,ω是角频率,σ为电导率,μ为磁导率,μ0为真空磁导率。
通过麦克斯韦方程组和边界条件的离散化,可以得到该节点的场值表达式,Eseafloor-x如下所示:
其中seafloor表示海底及地层网格节点,sea表示海水层网格节点;
使用电场值E,通过麦克斯韦方程可以求出各网格节点的磁场值H,便可以得到模型的阻抗响应Z。
海底大地电磁非线性共轭梯度三维反演:
海底三维反演的目标函数表示为:
其中和Zseafloor-n在这里专指海底观测数据和正演数据,εseafloor-n为观测数据误差,第1-N个数据为阻抗各分量的实部,第N+1-2N个数据为阻抗各分量的虚部,W=Wair+Wsea+Wseafloor为预调节矩阵,m=mair+msea+mseafloor为网格剖分模型,air表示空气层部分,sea表示海水层部分,seafloor表示海底及地下半空间部分。
目标函数梯度表示为:
三维阻抗灵敏度张量表示为:
其中,K=Kair+Ksea+Kseafloor为正演系数矩阵,E=Eair+Esea+Eseafloor为场值并分为两种源条件下的场值,
为正演模型有限差分网格至接收点插入两种极化源电磁场向量的线性组合。
得到:
p(i)=-C(i) -1g(i)+β(i)p(i-1) (14)
其中γ为模型电阻率和数据误差相关的变量,从而得到更新的模型:
mseafloor-(i+1)=mseafloor-(i)+α(i)p(i) (18)
重复上述迭代过程直至达到最小,便可得到最优模型。
在OPENMP并行技术基础上,实现了该反演方法的并行运算算法,将不同频点的数据分派给不同的CPU核心分别计算目标函数的梯度,大幅提高了计算效率
本发明实施例中,
第一步,接收反演参数、初始模型参数和数据及误差等参数;
第二步,根据应用计算机的cpu数量自动将每个频点的数据分配给所有的cpu进行正演和目标函数梯度的并行计算;
第三步,并行计算结束,合并数据,计算查找方向和查找步长;
第四步,计算模型改变量并对模型参数进行更新;
第五步,当目标函数足够小时结束迭代,否则令i=i+1,进入第二步。
本发明采用将各频点数据分散于不同线程进行并行运算不会影响计算精度。本发明根据计算机线程数量将各频点以及网格节点数据分散传输于不同的线程,由于该并行反演方法是高效低损耗的方法,对内存和CPU的要求较低,因此本发明方法能够适用于目前主流的PC机上,基本上不受硬件条件制约,能够在多种操作系统上推广使用。
并行的非线性共轭梯度反演算法中的并行计算部分主要是通过类似正演模拟的计算求解不同频点对应的目标函数梯度,因此对于各个计算机线程,分别计算对应频点的查找步长和查找方向,并行计算结束后进行参数规整求和即可。
下面设置不同的海洋地下空间模型,使用交错网格有限差分法计算正演响应并模拟反演使用的实测数据,建立加密的反演初始模型,应用本发明方法进行反演。
模型一
图4为本发明方法的应用例1的模型示意图,如图4所示,(网格为5*5*6,测点9个,频点为0.01s、0.1s、1s、10s、100s)反演结果如下:
数据频率为0.01s、0.1s、1s、10s、100s,9个测点位于网格中心,迭代82次后满足终止条件,RMS为0.056,主频为2.93GHz的8线程计算机运行时间为42s,(a)为y、z方向的电阻率断面,中心600m*1000m范围内电阻率为9Ω·m,其上方网格电阻率为95Ω·m,下方网格电阻率为75Ω·m,左右网格电阻率为91Ω·m;(b)为x、z方向的电阻率断面,中心600m*1000m范围内电阻率为9Ω·m,其上方网格电阻率为99Ω·m,下方网格电阻率为78Ω·m,左右网格电阻率为88Ω·m。可以看出本发明方法的反演收敛速度较快,计算效率和计算精度均有较好表现。
模型二
图5为本发明方法的应用例2的模型示意图,如图5所示,模型包括海水层和地下半空间,海水层电阻率0.3Ω·m,厚度1.8km,地下半空间电阻率100Ω·m,包含2个低阻异常体和2个高阻异常体,异常体电阻率分别为10Ω·m和1000Ω·m,异常体顶深(异常体顶界到海底的距离)1.8km,长10km,宽10km、厚3km、间距4km,数据频率范围0.1-1000s,共40个频点,网格数为30*30*40,400个测点位于网格中心。在本应用例中迭代60次后满足终止条件,RMS为0.2,主频为2.63GHz的4核计算机运行时间为2小时,由结果可知(图6所示),两个低阻异常体和两个高阻异常体均能够较好的分辨出来,异常体中心电阻率分别为9Ω·m和300Ω·m。
以上实施例仅为本发明的示例性实施例,不用于限制本发明,本发明的保护范围由权利要求书限定。本领域技术人员可以在本发明的实质和保护范围内,对本发明做出各种修改或等同替换,这种修改或等同替换也应视为落在本发明的保护范围内。
Claims (2)
1.海洋大地电磁场非线性共轭梯度三维并行反演方法,其特征在于,包括如下步骤:
a.在笛卡尔坐标系下沿x、y、z三个坐标轴将模型空间划分成Nx、Ny、Nz个小的长方体网格单元,间距为Δxi、Δyj、Δzk,其中i=1,...,Nx,j=1,...,Ny,k=1,...,Nz,模型空间包括空气层、海水层和地下半空间;海底节点(i,j,k)的长方体网格单元的长度、宽度和高度分别为xseafloor-i、yseafloor-j、zseafloor-k,电阻率为ρseafloor-i,j,k,其中电场取在各网格单元边缘的中点,用和表示;磁场取在各网格单元表面的中心,用和表示;
b.通过麦克斯韦方程组和边界条件的离散化得到模型各节点的场值表达式及其计算,实现正演;
c.将正演得到的阻抗与实测阻抗对比,获得阻抗数据偏差;
d.通过阻抗数据偏差计算目标函数,判断目标函数值是否足够小,当目标函数值足够小时结束迭代,否则进入下一步;
e.计算目标函数的梯度,以非线性共轭梯度法作为反演方法获得查找步长和查找方向,从而计算出模型各节点的模型修改量,其中使用阻抗数据误差和当前迭代计算所得的电阻率计算反演方法中的预处理因子;
f.使用模型修改量修正模型,并返回步骤b,
其中,
节点的场值表达式,Eseafloor-x如下所示:
其中seafloor表示海底及地层网格节点,sea表示海水层网格节点,ω是角频率,μ0为真空磁导率;
使用电场值E,通过麦克斯韦方程能够求出各网格节点的磁场值H,便能得到模型的阻抗响应Z:
2.根据权利要求1所述的海洋大地电磁场非线性共轭梯度三维并行反演方法,其特征在于,反演目标函数为:
其中,和Zseafloor-n专指海底观测数据和正演数据,εseafloor-n为观测数据误差,第1~N个数据为阻抗各分量的实部,第N+1~2N个数据为阻抗各分量的虚部,W=Wair+Wsea+Wseafloor为预调节矩阵,m=mair+msea+mseafloor为网格剖分模型,air表示空气层部分,sea表示海水层部分,seafloor表示海底及地下半空间部分,目标函数梯度表示为:
三维阻抗灵敏度张量表示为:
其中,K=Kair+Ksea+Kseafloor为正演系数矩阵,E=Eair+Esea+Eseafloor为场值并分为两种源条件下的场值,
为正演模型有限差分网格至接收点插入两种极化源电磁场向量的线性组合;
得到:
p(i)=-C(i) -1g(i)+β(i)p(i-1) (14)
其中γ为模型电阻率和数据误差相关的变量,从而得到更新的模型:
mseafloor-(i+1)=mseafloor-(i)+α(i)p(i) (18)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610268705.0A CN106019394B (zh) | 2016-04-27 | 2016-04-27 | 海洋大地电磁场非线性共轭梯度三维并行反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610268705.0A CN106019394B (zh) | 2016-04-27 | 2016-04-27 | 海洋大地电磁场非线性共轭梯度三维并行反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106019394A CN106019394A (zh) | 2016-10-12 |
CN106019394B true CN106019394B (zh) | 2019-04-05 |
Family
ID=57081662
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610268705.0A Active CN106019394B (zh) | 2016-04-27 | 2016-04-27 | 海洋大地电磁场非线性共轭梯度三维并行反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106019394B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108508490B (zh) * | 2018-03-07 | 2019-12-27 | 吉林大学 | 一种基于解析信号的磁张量梯度数据均衡边界识别方法 |
CN108562896A (zh) * | 2018-04-23 | 2018-09-21 | 广西师范大学 | 一种基于三维正压浅海大陆架模型的深层海流反演方法 |
CN109343347B (zh) * | 2018-10-12 | 2021-10-29 | 哈尔滨工程大学 | 一种海底飞行节点的轨迹跟踪控制方法 |
CN111856596B (zh) * | 2020-08-05 | 2023-03-28 | 中国海洋大学 | 层状介质电阻率各向异性海洋可控源电磁快速反演方法 |
CN111856597B (zh) * | 2020-08-05 | 2023-03-21 | 中国海洋大学 | 拖曳式海洋电磁地层电阻率与接收站位置联合反演方法 |
CN113094915B (zh) * | 2021-04-21 | 2022-05-27 | 中国海洋大学 | 一种海流感生磁场计算方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090150124A1 (en) * | 2007-12-07 | 2009-06-11 | Schlumberger Technology Corporation | Model based workflow for interpreting deep-reading electromagnetic data |
CN102798898A (zh) * | 2012-08-20 | 2012-11-28 | 中国地质科学院矿产资源研究所 | 大地电磁场非线性共轭梯度三维反演方法 |
CN102798897A (zh) * | 2012-08-20 | 2012-11-28 | 中国地质科学院矿产资源研究所 | 坑-井地大地电磁场非线性共轭梯度二维反演方法 |
CN104375195A (zh) * | 2013-08-15 | 2015-02-25 | 中国石油天然气集团公司 | 时频电磁的多源多分量三维联合反演方法 |
-
2016
- 2016-04-27 CN CN201610268705.0A patent/CN106019394B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090150124A1 (en) * | 2007-12-07 | 2009-06-11 | Schlumberger Technology Corporation | Model based workflow for interpreting deep-reading electromagnetic data |
CN102798898A (zh) * | 2012-08-20 | 2012-11-28 | 中国地质科学院矿产资源研究所 | 大地电磁场非线性共轭梯度三维反演方法 |
CN102798897A (zh) * | 2012-08-20 | 2012-11-28 | 中国地质科学院矿产资源研究所 | 坑-井地大地电磁场非线性共轭梯度二维反演方法 |
CN104375195A (zh) * | 2013-08-15 | 2015-02-25 | 中国石油天然气集团公司 | 时频电磁的多源多分量三维联合反演方法 |
Non-Patent Citations (4)
Title |
---|
"一种并行的大地电磁场非线性共轭梯度三维反演方法";张昆 等;《地球物理学报》;20131130;第56卷(第11期);第3922-3931页 * |
"大地电磁场非线性共轭梯度三维反演";张昆 等;《矿床地质》;20121231;第31卷;第891-892页 * |
"海岸效应对近海地区大地电磁测深数据畸变作用研究";张帆 等;《地球物理学报》;20121231;第55卷(第12期);第4023-4035页 * |
"海底大地电磁场非线性共轭梯度二维反演";张昆;《中国优秀硕士学位论文全文数据库-基础科学辑》;20150815(第8期);正文第1-2、20-21页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106019394A (zh) | 2016-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106019394B (zh) | 海洋大地电磁场非线性共轭梯度三维并行反演方法 | |
CN102798898B (zh) | 大地电磁场非线性共轭梯度三维反演方法 | |
Abubakar et al. | Joint inversion approaches for geophysical electromagnetic and elastic full-waveform data | |
Streich | 3D finite-difference frequency-domain modeling of controlled-source electromagnetic data: Direct solution and optimization for high accuracy | |
Grayver et al. | 3D inversion and resolution analysis of land-based CSEM data from the Ketzin CO 2 storage formation | |
CN110058317A (zh) | 航空瞬变电磁数据和航空大地电磁数据联合反演方法 | |
Zhdanov | New advances in regularized inversion of gravity and electromagnetic data | |
CN101194262B (zh) | 在海洋电磁勘测中判断地球垂直电各向异性的方法 | |
CN106980736B (zh) | 一种各向异性介质的海洋可控源电磁法有限元正演方法 | |
CN109633762A (zh) | 基于正弦函数的相关性约束条件联合反演重磁数据的方法 | |
Li et al. | Inversion of controlled‐source electromagnetic data using a model‐based approach | |
CN111638556A (zh) | 基于地空分解策略的大地电磁正演方法及装置、存储介质 | |
CN104123455B (zh) | 大地电磁场非线性共轭梯度三维倾子反演方法 | |
CN105044782B (zh) | 一种海洋地下介质总有机碳含量的获取方法 | |
Li et al. | 3-D marine CSEM forward modeling with general anisotropy using an adaptive finite-element method | |
Yin et al. | Forward modeling of marine DC resistivity method for a layered anisotropic earth | |
Elías et al. | Three-dimensional modelling of controlled source electro-magnetic surveys using non-conforming finite element methods | |
Endo et al. | A multigrid integral equation method for large-scale models with inhomogeneous backgrounds | |
US10656300B2 (en) | Inversion method | |
Jang et al. | Three-dimensional finite-difference modeling of time-domain electromagnetic responses for a large-loop source | |
Bobrov et al. | Physical and mathematical modeling of transient electromagnetic soundings over salt-dome structures | |
Zhang | Parallel goal-oriented adaptive finite element modeling for 3D electromagnetic exploration | |
Zhao et al. | Transient electromagnetic detection of unexploded ordnance buried in underwater sediments | |
Palshin et al. | On-shore amplification of the electric field induced by a coastal sea current | |
CN115563791B (zh) | 基于压缩感知重构的大地电磁数据反演方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |