CN113933905A - 一种圆锥型场源瞬变电磁反演方法 - Google Patents
一种圆锥型场源瞬变电磁反演方法 Download PDFInfo
- Publication number
- CN113933905A CN113933905A CN202111157516.3A CN202111157516A CN113933905A CN 113933905 A CN113933905 A CN 113933905A CN 202111157516 A CN202111157516 A CN 202111157516A CN 113933905 A CN113933905 A CN 113933905A
- Authority
- CN
- China
- Prior art keywords
- inversion
- model
- field source
- data
- transient electromagnetic
- 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 32
- 230000001052 transient effect Effects 0.000 title claims abstract description 28
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 14
- 238000004364 calculation method Methods 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 23
- 230000004044 response Effects 0.000 claims description 15
- 230000005674 electromagnetic induction Effects 0.000 claims description 7
- 230000006698 induction Effects 0.000 claims description 7
- 230000005672 electromagnetic field Effects 0.000 claims description 6
- 230000035699 permeability Effects 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 230000005284 excitation Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000011218 segmentation Effects 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 6
- 238000001514 detection method Methods 0.000 abstract description 5
- 238000003672 processing method Methods 0.000 abstract description 3
- 239000000779 smoke Substances 0.000 description 4
- 238000005457 optimization Methods 0.000 description 3
- 241000219000 Populus Species 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 238000012897 Levenberg–Marquardt algorithm Methods 0.000 description 1
- 241001415846 Procellariidae Species 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
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/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
- G01V3/26—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with magnetic or electric fields produced or modified either by the surrounding earth formation or by the detecting device
- G01V3/28—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with magnetic or electric fields produced or modified either by the surrounding earth formation or by the detecting device using induction coils
-
- 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
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- 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)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Remote Sensing (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Geophysics (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Geology (AREA)
- Computational Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Geometry (AREA)
- Computing Systems (AREA)
- Electromagnetism (AREA)
- Evolutionary Computation (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种圆锥型场源瞬变电磁反演方法,属于瞬变电磁探测技术领域。本发明基于OCCAM反演算法,提出了一种圆锥型场源瞬变电磁反演方法。OCCAM法是一种带平滑约束的最小二乘反演,受初始模型影响小,运算稳定收敛,是一种有效的数据反演处理方法。本发明的反演方法不需要依赖于初始模型。对初始模型没有特别要求,且计算效率高。本发明采用OCCAM法得到的反演电阻率更接近模型电阻率。本发明反演采用简化后的拉格朗日乘子算法使得反演仍稳定、速度更快、结果更真实。本发明方法使得反演深度更接近层界面,尤其对目标层底界面的分辨能力获得加强,表现出较好的分层效果。
Description
技术领域
本发明涉及一种圆锥型场源瞬变电磁反演方法,具体涉及一种基于OCCAM反演算法的圆锥型场源瞬变电磁数据的反演方法,属于瞬变电磁探测技术领域。
背景技术
瞬变电磁法(TEM)逐渐被应用于矿山巷道、交通隧道和地下空间探测,受到有限地下空间的限制,只能采用小尺寸发射-接收装置,如多匝小回线装置。杨海燕等发明了一种圆锥型场源装置,如图1所示,与多匝小回线装置相比,圆锥型装置具有互感低、关断时间短和“盲区”范围小的优点。但是,用于处理圆锥型场源探测资料的方法较少。
目前,用于圆锥型场源瞬变电磁的资料处理方法,主要有如下几种:
(1)杨海燕,李锋平,岳建华,刘旭华,赵海娇.基于“烟圈”理论的圆锥型场源瞬变电磁优化反演.中国矿业大学学报,2016,45(06):1230-1237.该论文提出了将“烟圈”反演结果转换到频率域作为初始模型,采用最小二乘正则化方法对反演结果进行改进的反演方案。
(2)杨海燕,李锋平,Chen Shen-En,岳建华,郭福生,陈晓,张华.An inversion oftransient electromagnetic data from a conical source.Applied Geophysics,2018,15(Z1):545-555+570.该论文提出了根据“烟圈”反演结果设定初始模型的阻尼最小二乘反演方法。对瞬变电磁正演所获得的视电阻率进行“烟圈”反演,以“烟圈”反演结果设定的初始模型在频率域内进行了阻尼最小二乘反演。同时,提出了以“烟圈”反演曲线的极值点和拐点为基点,建立阻尼最小二乘反演初始模型的方法。
(3)黄刚,杨海燕,李锋平,彭峰,焦俊俊,刘建鹏.圆锥型场源瞬变电磁数据AWPSO算法优化反演.地球物理学进展:1-13[2021-09-08].http://kns.cnki.net/kcms/detail/11.2982.p.20201109.1151.072.html.该论文提出了采用粒子群优化算法(PSO)和神经网络算法(BP)对圆锥型场源探测数据进行反演的方案,改进了一种基于神经网络算法Sigmoid函数的自适应加权粒子群优化(AWPSO)算法。
但上述反演方法的分层能力比较弱,为了使反演深度更接近层界面,加强对目标层底界面的分辨能力,获得较好的分层效果,需要提出一种新的圆锥型场源瞬变电磁反演方法。
发明内容
本发明所要解决的技术问题是:提供一种圆锥型场源瞬变电磁反演方法,该方法使得反演深度更接近层界面,尤其对目标层底界面的分辨能力获得加强,表现出较好的分层效果。
本发明为解决上述技术问题采用以下技术方案:
一种圆锥型场源瞬变电磁反演方法,该方法包括如下步骤:
步骤2,对于第k次迭代,计算当前模型mk对应的正演数据F[mk]、雅克比矩阵Jk以及数据拟合差X2;
步骤3,搜索使数据拟合差X2最小时的模型mk+1,同时计算模型mk+1对应的粗糙度矩阵Rk+1以及步长ak;
步骤4,判断步骤3得到的最小数据拟合差是否小于反演所要求达到的拟合差若是,则进入步骤5,否则,判断步骤3得到的最小数据拟合差是否小于第k-1次迭代得到的最小数据拟合差;若是,则进入步骤6,否则使步长ak减半并返回步骤3;
步骤5,确定当数据拟合差X2等于反演所要求达到的拟合差时的拉格朗日乘子及对应的模型m′k+1,根据模型m′k+1计算粗糙度矩阵R′k+1,判断R′k+1是否大于Rk+1,若是,则使步长ak减半并返回步骤3,否则进入步骤6;
步骤6,保存步长ak和模型mk+1;
作为本发明的一种优选方案,步骤1所述根据初始模型m1构造粗糙度矩阵为:
N为地层模型分段得到的总层数。
作为本发明的一种优选方案,所述步骤2的具体过程如下:
对于第k次迭代,向圆锥型场源通入电流,对于圆锥型场源的第i匝圆形发射线圈,其中心处垂直方向磁场的频率域响应为:
式中,H(ω)表示频率域响应,ri为第i匝圆形发射线圈的半径,I(ω)为发射电流的傅式谱,λ为积分变量,hi为第i匝圆形发射线圈离地高度,u0为真空中的磁导率,Z0为表层介质的输入波阻抗,Z(1)为第一层的输入阻抗,J1为第一类零阶贝塞尔函数;Z1=-iωμ0/u1,k1为第一层介质的波数,σ1为第一层介质的电导率,ω为角频率,μ0为真空中磁导率;
对频率域响应公式进行140点汉克尔变换得到:
设定激励电流I(t)为单位阶跃电流:
将140点汉克尔变换得到的频率域,通过傅里叶变换求得时间域中响应:
式中,t为时间,f(t)为时间域电磁场分布,F(ω)为频率域的电磁场分布,j为虚数单位;
将频率域响应公式代入上式即得时间域的磁场响应:
根据法拉第电磁感应定律,得到二次场的感应电位为vi,将圆锥型场源产生的瞬变电磁感应电位视为n个线圈感应电位的叠加,则圆锥型场源下层状模型的瞬变电磁感应电位V(t)为:
式中,n为圆锥型场源缠绕的圆形发射线圈匝数,Δti是第i匝线圈与第1匝线圈距离产生的延时,对于第1匝线圈,Δt1=0;采用全区视电阻率算法将感应电位V(t)转换成反映电性变化的的视电阻率,即当前模型对应的正演数据;
雅克比矩阵Jk为:
数据拟合差X2为:
X2=‖Wd-WF[mk]‖2
其中,W=diag{1/σ1,1/σ2,…,1/σM},σj为第j个数据的标准差,j=1,…,M,M为场源装置测得的数据个数,d为场源装置测得的数据集合。
作为本发明的一种优选方案,步骤3所述步长ak的计算公式为:
作为本发明的一种优选方案,所述模型mk+1定义为:
其中,为第k+1次迭代的拉格朗日乘子,为N×N的矩阵,N为地层模型分段得到的总层数,W=diag{1/σ1,1/σ2,…,1/σM},σj为第j个数据的标准差,j=1,…,M,M为场源装置测得的数据个数,Jk为雅克比矩阵,dk=d-F[mk]+Jkmk,d为场源装置测得的数据集合,F[mk]表示当前模型mk对应的正演数据。
本发明采用以上技术方案与现有技术相比,具有以下技术效果:
1、本发明的反演方法不需要依赖于初始模型,对初始模型没有特别要求,且计算效率高,较好地反映了模型的地质结构,对模型具有平滑效果。
2、本发明采用OCCAM法得到的反演电阻率更接近各层介质的真实电阻率,对电阻率的改进效果良好。
3、本发明反演采用简化后的拉格朗日乘子算法使得反演仍稳定、速度更快、结果更真实。
4、本发明在一定步长下逐次递减拉格朗日乘子,使反演仍稳定,速度更快,结果更真实。
附图说明
图1是一种圆锥型场源装置的结构图。
图2是本发明一种圆锥型场源瞬变电磁反演方法的流程图。
图3是D型模型反演曲线。
图4是G型模型反演曲线。
图5是H型模型反演曲线。
图6是K型模型反演曲线。
图7是KH型模型反演曲线。
具体实施方式
下面详细描述本发明的实施方式,所述实施方式的示例在附图中示出。下面通过参考附图描述的实施方式是示例性的,仅用于解释本发明,而不能解释为对本发明的限制。
本发明基于OCCAM反演算法,提出了一种圆锥型场源瞬变电磁反演方法,具体流程如图2所示。OCCAM法是一种带平滑约束的最小二乘反演,受初始模型影响小,运算稳定收敛,是一种有效的数据反演处理方法。
首先,定义粗糙度,代表了模型的光滑程度。粗糙度可表示为模型相对某一坐标的一阶或二阶导数的平方的积分,如对z方向,则:
式中,z为深度,m为模型相应深度下所对应的电性参数,一般是地层的电阻率,或者电阻率的对数。模型的最底部为均匀半空间,将地层模型连续分段处理得到:
m(z)=mi,zi-1<z<zi+1,i=1,2,…,N (2)
式中,z0=0,zi-1/zi取为小于1的常数,mi是电阻率或电导率值。
将粗糙度的微分形式离散化,得到:
粗糙度用矩阵表示为:
根据加权最小二乘法原理,定义模型数据与实测数据的拟合差X为:
式中,M为观测数据个数,dj是第j个观测数据,σj为第j个数据的标准差,Fj[m]为模型m对应的正演数据。
正演部分如下:
对于半径ri,通入电流I的第i匝圆形发射线圈,其中心处垂直方向磁场的频率域响应为:
式中,λ为积分变量,I(w)为发射电流的傅式谱,hi为第i匝线圈离地高度,Z(1)为第一层的输入阻抗,J1为第一类零阶贝塞尔函数。对公式(7)进行汉克尔变换得到:
式中,I为回线源中通入的电流,m为抽样点的位置,在140点J1汉克尔变换的滤波器系数中,a=-7.91001919000e+00,s=8.79671439570e-02。然后进行频率域向时间域的转换。假设激励电流为单位阶跃电流:
公式(8)通过傅里叶变换求得时域中响应:
将公式(7)代入上式即可得到时间域的磁场响应:
式中Im为虚部,Re为实部,为角频率,t为时间。式(11)、(12)含有特殊函数,内层贝塞尔函数积分采用140点快速汉克尔变换计算。根据法拉第电磁感应定律,可得到二次场的感应电位为vi。圆锥形场源产生的瞬变电磁感应电位可视为多个线圈感应电位的叠加,因此,得到该场源下层状模型的瞬变电磁感应电位为:
式中,Δti是是各匝线圈与第1匝线圈距离不同而产生的延时,对于第1匝线圈i=1,Δt1=0。进而,采用全区视电阻率算法将感应电位V(t)转换成反映电性变化的的视电阻率ρa(t)。
这样求解的数学问题就转化为:对于给定的数据集d及其误差,寻求当X2达到可接受的值时R1或R2尽可能小的模型m。
数据拟合差写为:
X2=‖Wd-WF[m]‖2 (14)
式中,W=diag{1/σ1,1/σ2,…,1/σM}。
引入拉格朗日乘子μ-1,构造无约束的目标函数:
式中,J是M×N阶的雅克比矩阵:
在迭代过程开始前,给定一个初始模型m1,如果F在m1可微,则对于足够小的向量Δ有:
F[m1+Δ]=F[m1]+J1Δ (18)
式中,Δ=m2-m1,将式(18)代入式(15),则在m2会出现如下的线性问题:
上式右端第二项圆括号中的表达式是一种数据向量,记作d1。
进而,定义模型m2,使U取得最小值,有
选取μ以达到期望的拟合差,在迭代中m2由系列m1,m2,…中后面的模型依次取代,系列中每一个前面的向量都是下个向量的初始近似值。
当第k次迭代完成后,定义向量:
用一系列的μ值计算模型mk+1(μ)的真正拟合差:
Xk+1(μ)=‖Wd-WF[mk+1(μ)]‖ (22)
为防止拟合问题发散,定义拉格朗日乘子函数的新模型集,即:
G(μ)=(1-a)mk+amk+1(μ) (23)
其中,a为步长,其将被连续地平分直到获得一个合适的模型。
层状模型更能检验反演算法的有效性,尤其能反映出算法对层界面的分辨能力。分别建立典型的2层(D型和G型)、3层(H型、K型)以及4层(KH型)介质模型来验证。图3、图4、图5、图6和图7分别为层状模型的反演结果。图中左上部分是预设的初始模型及各项参数,左下为模型的一维正演结果,横坐标为时间,纵坐标为电阻率;右边为反演结果,横坐标为视电阻率,纵坐标为深度。图中反演曲线的形态反映出对应介质模型的特征,且对层界面的分辨较好。
以上实施例仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明保护范围之内。
Claims (5)
1.一种圆锥型场源瞬变电磁反演方法,其特征在于,该方法包括如下步骤:
步骤2,对于第k次迭代,计算当前模型mk对应的正演数据F[mk]、雅克比矩阵Jk以及数据拟合差X2;
步骤3,搜索使数据拟合差X2最小时的模型mk+1,同时计算模型mk+1对应的粗糙度矩阵Rk+1以及步长ak;
步骤4,判断步骤3得到的最小数据拟合差是否小于反演所要求达到的拟合差若是,则进入步骤5,否则,判断步骤3得到的最小数据拟合差是否小于第k-1次迭代得到的最小数据拟合差;若是,则进入步骤6,否则使步长ak减半并返回步骤3;
步骤5,确定当数据拟合差X2等于反演所要求达到的拟合差时的拉格朗日乘子及对应的模型m′k+1,根据模型m′k+1计算粗糙度矩阵R′k+1,判断R′k+1是否大于Rk+1,若是,则使步长ak减半并返回步骤3,否则进入步骤6;
步骤6,保存步长ak和模型mk+1;
3.根据权利要求1所述圆锥型场源瞬变电磁反演方法,其特征在于,所述步骤2的具体过程如下:
对于第k次迭代,向圆锥型场源通入电流,对于圆锥型场源的第i匝圆形发射线圈,其中心处垂直方向磁场的频率域响应为:
式中,H(ω)表示频率域响应,ri为第i匝圆形发射线圈的半径,I(ω)为发射电流的傅式谱,λ为积分变量,hi为第i匝圆形发射线圈离地高度,u0为真空中的磁导率,Z0为表层介质的输入波阻抗,Z(1)为第一层的输入阻抗,J1为第一类零阶贝塞尔函数;Z1=-iωμ0/u1,k1为第一层介质的波数,σ1为第一层介质的电导率,ω为角频率,μ0为真空中磁导率;
对频率域响应公式进行140点汉克尔变换得到:
设定激励电流I(t)为单位阶跃电流:
将140点汉克尔变换得到的频率域,通过傅里叶变换求得时间域中响应:
式中,t为时间,f(t)为时间域电磁场分布,F(ω)为频率域的电磁场分布,j为虚数单位;
将频率域响应公式代入上式即得时间域的磁场响应:
根据法拉第电磁感应定律,得到二次场的感应电位为vi,将圆锥型场源产生的瞬变电磁感应电位视为n个线圈感应电位的叠加,则圆锥型场源下层状模型的瞬变电磁感应电位V(t)为:
式中,n为圆锥型场源缠绕的圆形发射线圈匝数,Δti是第i匝线圈与第1匝线圈距离产生的延时,对于第1匝线圈,Δt1=0;采用全区视电阻率算法将感应电位V(t)转换成反映电性变化的的视电阻率,即当前模型对应的正演数据;
雅克比矩阵Jk为:
数据拟合差X2为:
X2=||Wd-WF[mk]||2
其中,W=diag{1/σ1,1/σ2,...,1/σM},σj为第j个数据的标准差,j=1,...,M,M为场源装置测得的数据个数,d为场源装置测得的数据集合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111157516.3A CN113933905B (zh) | 2021-09-30 | 2021-09-30 | 一种圆锥型场源瞬变电磁反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111157516.3A CN113933905B (zh) | 2021-09-30 | 2021-09-30 | 一种圆锥型场源瞬变电磁反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113933905A true CN113933905A (zh) | 2022-01-14 |
CN113933905B CN113933905B (zh) | 2023-09-19 |
Family
ID=79277761
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111157516.3A Active CN113933905B (zh) | 2021-09-30 | 2021-09-30 | 一种圆锥型场源瞬变电磁反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113933905B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110927458A (zh) * | 2019-11-11 | 2020-03-27 | 中国电子科技集团公司第十一研究所 | 多载流子体系的测试及拟合方法 |
CN114675337A (zh) * | 2022-03-22 | 2022-06-28 | 扬州蓝德森科技有限公司 | 一种基于多匝线圈和瞬变电磁法的地下测深方法 |
CN115508898A (zh) * | 2022-10-11 | 2022-12-23 | 四川中水成勘院工程物探检测有限公司 | G-s变换的接地长导线源瞬变电磁快速正反演方法及系统 |
CN115983105A (zh) * | 2022-12-12 | 2023-04-18 | 成都理工大学 | 基于深度学习加权决策的Occam反演拉格朗日乘子优化方法 |
CN116879964A (zh) * | 2023-08-14 | 2023-10-13 | 成都理工大学 | 一种时频电磁频率域数据自约束稳健电阻率反演方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005085909A1 (en) * | 2004-02-24 | 2005-09-15 | Kjt Enterprises, Inc. | Combined surface and wellbore electromagnetic measurement system and method for determining formation fluid properties |
US20100235100A1 (en) * | 2009-03-16 | 2010-09-16 | Bruce Alan Hobbs | Method for determining resistivity anisotropy from earth electromagnetic responses |
CN104914473A (zh) * | 2015-05-27 | 2015-09-16 | 中国石油天然气集团公司 | 一种瞬变电磁电阻率的反演方法与装置 |
CN106501867A (zh) * | 2016-10-19 | 2017-03-15 | 中国科学院电子学研究所 | 一种基于横向平滑约束的瞬变电磁反演方法 |
US20170075030A1 (en) * | 2015-09-15 | 2017-03-16 | Brent D. Wheelock | Accelerated Occam Inversion Using Model Remapping and Jacobian Matrix Decomposition |
CN108680966A (zh) * | 2018-03-21 | 2018-10-19 | 中国石油大学(华东) | 海洋可控源电磁勘探噪声降噪效果评估方法 |
CN209167552U (zh) * | 2019-01-08 | 2019-07-26 | 东华理工大学 | 一种瞬变电磁法圆锥形场源装置 |
CN110764154A (zh) * | 2019-11-08 | 2020-02-07 | 吉林大学 | 基于改进粒子群算法的时间域航空电磁一维反演方法 |
AU2020102911A4 (en) * | 2020-10-21 | 2020-12-17 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Ground-borehole electrical source short offset transient electromagnetic detection method |
CN112949134A (zh) * | 2021-03-09 | 2021-06-11 | 吉林大学 | 基于非结构有限元方法的地-井瞬变电磁反演方法 |
-
2021
- 2021-09-30 CN CN202111157516.3A patent/CN113933905B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005085909A1 (en) * | 2004-02-24 | 2005-09-15 | Kjt Enterprises, Inc. | Combined surface and wellbore electromagnetic measurement system and method for determining formation fluid properties |
US20100235100A1 (en) * | 2009-03-16 | 2010-09-16 | Bruce Alan Hobbs | Method for determining resistivity anisotropy from earth electromagnetic responses |
CN104914473A (zh) * | 2015-05-27 | 2015-09-16 | 中国石油天然气集团公司 | 一种瞬变电磁电阻率的反演方法与装置 |
US20170075030A1 (en) * | 2015-09-15 | 2017-03-16 | Brent D. Wheelock | Accelerated Occam Inversion Using Model Remapping and Jacobian Matrix Decomposition |
CN106501867A (zh) * | 2016-10-19 | 2017-03-15 | 中国科学院电子学研究所 | 一种基于横向平滑约束的瞬变电磁反演方法 |
CN108680966A (zh) * | 2018-03-21 | 2018-10-19 | 中国石油大学(华东) | 海洋可控源电磁勘探噪声降噪效果评估方法 |
CN209167552U (zh) * | 2019-01-08 | 2019-07-26 | 东华理工大学 | 一种瞬变电磁法圆锥形场源装置 |
CN110764154A (zh) * | 2019-11-08 | 2020-02-07 | 吉林大学 | 基于改进粒子群算法的时间域航空电磁一维反演方法 |
AU2020102911A4 (en) * | 2020-10-21 | 2020-12-17 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Ground-borehole electrical source short offset transient electromagnetic detection method |
CN112949134A (zh) * | 2021-03-09 | 2021-06-11 | 吉林大学 | 基于非结构有限元方法的地-井瞬变电磁反演方法 |
Non-Patent Citations (12)
Title |
---|
HAIYAN YANG,等: "Transient Electromagnetic Response With a Ramp Current Excitation Using Conical Source", 《IEEE》 * |
HAIYAN YANG,等: "Transient Electromagnetic Response With a Ramp Current Excitation Using Conical Source", 《IEEE》, 6 May 2019 (2019-05-06), pages 63829 - 63836 * |
R. A. S. GEHRMANN等: "Marine Mineral Exploration With Controlled Source Electromagnetics at the TAG Hydrothermal Field, 26°N Mid‐Atlantic Ridge", 《GEOPHYSICAL RESEARCH LETTERS》 * |
R. A. S. GEHRMANN等: "Marine Mineral Exploration With Controlled Source Electromagnetics at the TAG Hydrothermal Field, 26°N Mid‐Atlantic Ridge", 《GEOPHYSICAL RESEARCH LETTERS》, 5 June 2019 (2019-06-05), pages 5808 - 5816 * |
李刚等: "回线源瞬变电磁法一维反演算法", 《煤田地质与勘探》 * |
李刚等: "回线源瞬变电磁法一维反演算法", 《煤田地质与勘探》, no. 05, 31 December 2017 (2017-12-31), pages 165 - 176 * |
李峰平: "瞬变电磁法圆锥型场源响应特征理论研究", 《中国优秀博硕士学位论文全文数据库(硕士)》 * |
李峰平: "瞬变电磁法圆锥型场源响应特征理论研究", 《中国优秀博硕士学位论文全文数据库(硕士)》, 15 December 2017 (2017-12-15), pages 1 - 64 * |
杨海燕,等: "基于"烟圈"理论的圆锥型场源瞬变电磁优化反演", 《中国矿业大学学报》 * |
杨海燕,等: "基于"烟圈"理论的圆锥型场源瞬变电磁优化反演", 《中国矿业大学学报》, 30 November 2016 (2016-11-30), pages 1230 - 1237 * |
汪凌: "圆锥型场源瞬变电磁Occam反演技术研究", 《中国优秀博硕士学位论文全文数据库(硕士)》 * |
汪凌: "圆锥型场源瞬变电磁Occam反演技术研究", 《中国优秀博硕士学位论文全文数据库(硕士)》, 31 December 2022 (2022-12-31), pages 1 - 44 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110927458A (zh) * | 2019-11-11 | 2020-03-27 | 中国电子科技集团公司第十一研究所 | 多载流子体系的测试及拟合方法 |
CN114675337A (zh) * | 2022-03-22 | 2022-06-28 | 扬州蓝德森科技有限公司 | 一种基于多匝线圈和瞬变电磁法的地下测深方法 |
CN115508898A (zh) * | 2022-10-11 | 2022-12-23 | 四川中水成勘院工程物探检测有限公司 | G-s变换的接地长导线源瞬变电磁快速正反演方法及系统 |
CN115508898B (zh) * | 2022-10-11 | 2023-11-10 | 四川中水成勘院工程物探检测有限公司 | G-s变换的接地长导线源瞬变电磁快速正反演方法及系统 |
CN115983105A (zh) * | 2022-12-12 | 2023-04-18 | 成都理工大学 | 基于深度学习加权决策的Occam反演拉格朗日乘子优化方法 |
CN115983105B (zh) * | 2022-12-12 | 2023-10-03 | 成都理工大学 | 基于深度学习加权决策的Occam反演拉格朗日乘子优化方法 |
CN116879964A (zh) * | 2023-08-14 | 2023-10-13 | 成都理工大学 | 一种时频电磁频率域数据自约束稳健电阻率反演方法 |
CN116879964B (zh) * | 2023-08-14 | 2024-04-26 | 成都理工大学 | 一种时频电磁频率域数据自约束稳健电阻率反演方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113933905B (zh) | 2023-09-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113933905A (zh) | 一种圆锥型场源瞬变电磁反演方法 | |
CN105893678B (zh) | 一种时域有限差分的三维感应-极化双场数值模拟方法 | |
Fan et al. | Multidomain pseudospectral time-domain simulations of scattering by objects buried in lossy media | |
Wang et al. | Optimizing GIS partial discharge pattern recognition in the ubiquitous power internet of things context: A MixNet deep learning model | |
CN110058315B (zh) | 一种三维各向异性射频大地电磁自适应有限元正演方法 | |
Chiu et al. | Comparative study of some population-based optimization algorithms on inverse scattering of a two-dimensional perfectly conducting cylinder in dielectric slab medium | |
CN112949134B (zh) | 基于非结构有限元方法的地-井瞬变电磁反演方法 | |
CN110133733A (zh) | 一种基于粒子群优化算法的电导-极化率多参数成像方法 | |
Guo et al. | A three‐dimensional direct lightning strike model for lightning protection of the substation | |
CN104766090A (zh) | 一种基于bemd和sofm的探地雷达数据可视化方法 | |
Chiu et al. | Comparison of particle swarm optimization and asynchronous particle swarm optimization for inverse scattering of a two-dimensional perfectly conducting cylinder | |
CN106154259A (zh) | 一种随机集理论下的多传感器自适应管控方法 | |
CN106446470A (zh) | 一种高效并行的非均匀介质频域有限差分方法 | |
CN104849573A (zh) | 海洋环境下舰船平台强电磁脉冲环境仿真预测方法 | |
CN103235888B (zh) | 一种精确计算双各向异性介质球电磁散射的方法 | |
Li et al. | Estimation of frequency domain soil parameters of horizontally multilayered earth by using Cole–Cole model based on the parallel genetic algorithm | |
CN104915326A (zh) | 基于等效原理的区域分解阶数步进时域积分方法 | |
CN116879964B (zh) | 一种时频电磁频率域数据自约束稳健电阻率反演方法 | |
CN109710995A (zh) | 一种针对随机排布线缆的串扰噪声预测方法 | |
Xue et al. | The multilevel fast physical optics method for calculating high frequency scattered fields | |
CN115238566A (zh) | 一种基于频域电磁响应的人工空洞识别方法 | |
Wang et al. | GIS partial discharge pattern recognition via a novel capsule deep graph convolutional network | |
Wang et al. | Calculation of 3D transient Eddy current by the face‐smoothed finite element–boundary element coupling method | |
CN105760343A (zh) | 分析复杂多金属目标电磁散射特性的时域积分方法 | |
Xiong | 2.5 D forward for the transient electromagnetic response of a block linear resistivity distribution |
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 |