0 引言
瞬变电磁法作为一种常见的电磁测深物探方法,广泛应用于水文勘查、工程勘查以及矿产勘查[1 ,2 ] 。在目前的工程应用中,除“烟圈”反演快速成像法外,一维反演结果仍然是地质解释的主要依据。对于中心回线装置,可以解析推导出一维正演响应的计算公式[3 ] ,其数值计算的核心是,内层积分的汉克尔变换[4 ] 、外层积分的余弦变换[5 ] 均可采用数值滤波算法,随着滤波系数的不断推陈出新,计算量明显减少,大大提高了正演的计算效率和精度。在反演方面,目前常见的反演方法,如阻尼最小二乘法(马奎特法)[6 ] 、Occam法[7 ,8 ,9 ,10 ] 、正则化法[11 ,12 ,13 ] 等,都需要计算正演响应对模型的导数,即雅克比矩阵,因此雅克比矩阵的计算效率和效果就成为反演计算的关键因素之一。计算雅克比矩阵最简便的方法是差分法,该算法程序实现较为简单,缺点是不仅调用正演次数多,计算效率低,而且有可能增加累积误差。因此,本文详细推导了中心回线装置瞬变电磁一维正演全区视电阻率的雅克比矩阵的解析计算公式,并通过Occam反演,对比分析了理论模型和实测数据的反演效果,证明推导正确,值得采用。
1 瞬变电磁正演响应之全区视电阻率
目前,瞬变电磁仪器大多采集的数据为感应电动势,对其归一化后可转换为晚期视电阻率,然而晚期视电阻率是假设在时间趋于无穷大时近似推导所得,因此在响应早期,视电阻率曲线变形严重,呈明显高阻假象,成为了地质解释重要的干扰因素之一,因此很多学者提出了全区视电阻率(也称全期视电阻率)的概念[14 ] 。由于通过磁场强度计算的全区视电阻率更为接近真正的视电阻率定义[15 ] ,可通过感应电动势转换为磁场强度计算得出。因此,本文将该全区视电阻率作为Occam反演的拟合参数。
对于中心回线瞬变电磁的一维正演,可直接求解磁场强度。假设n 层地电模型,第j 层电阻率为ρj 、层厚度为hj ,频率域的磁场强度响应为:
(1) H z ( ω ) = aI ( ω ) ∫ 0 ∞ λ 2 λ + u ( 1 ) J 1 ( λa ) d λ ,
(2) u j = λ 2 - k j 2 , k j 2 = - i ω μ 0 / ρ j , j = 1,2 , 3 , … , n u ( j ) = u j u ( j + 1 ) + u j tanh ( u j h j ) u j + u ( j + 1 ) tanh ( u j h j ) , u ( n ) = u n 。
対于发射波形的选择,大多瞬变电磁仪普遍采用占空比1∶1的方波,如V8电法工作站的TD50方波,在正演模拟中只需考虑该波形的半个周期,可视为阶跃波:
(3) I ( t ) = I 0 t < 0 0 t > 0
式中I 0 为发射电流。对上式进行傅里叶变换,得到频率域中发射电流表达式:
(4) I ( ω ) = ∫ - ∞ ∞ I ( t ) e j ωt d t = - 1 j ω · I 0 。
(5) H z ( t ) = 2 π ∫ 0 ∞ Im H z ( ω ) ω cos ( ωt ) d ω = 2 I 0 π ∫ 0 ∞ Im 1 ω ∫ 0 ∞ λ 2 λ + u ( 1 ) J 1 ( λa ) d λ cos ( ωt ) d ω
令内层积分核函数为func= λ 2 λ + u ( 1 ) ,汉克尔变换用“HT {}”表示,余弦变换用“CT {}”表示,于是式(5)可简写成:
(6) H z ( t ) = 2 I 0 π · CT { HT { Im ( func ) } / ω } 。
对于均匀半空间地电模型,可以推导出Hz 的解析表达式[3 ] :
(7) H z ( t , x ) = I 0 2 a 3 π 1 x e - x 2 + 1 - 3 2 x 2 erf ( x ) ,
(8) Z ( x ) = 3 π 1 x e - x 2 + 1 - 3 2 x 2 erf ( x ) = 4 a π · CT { HT { Im ( func ) } / ω } ,
最后求取Z (x )的反函数x ,得出基于磁场强度的全区视电阻率:
(9) ρ a = μ 0 a 2 4 x 2 t 。
由于函数Z (x )为隐函数,无法推导出其反函数,因此无法通过解析算法求解,但可利用函数在0~ 1之间单调递增的特点,通过二分查找的办法得到数值解。笔者曾采用该方法实现了全区视电阻率的计算以及对应的“烟圈”反演,取得了不错的效果[15 ] 。
可以看到,瞬变电磁正演响应计算可简化为:首先计算核函数,然后再对内层积分(计算频率域的电磁响应)做汉克尔变换,对外层积分(频率域到时间域的转换)做余弦变换(式(6)),得到Hz 和Z (x ),最后采用二分查找法通过Z (x )求出x ,从而得到正演响应的全区视电阻率。
2 雅克比矩阵的解析推导
式(9)对模型参数求导,根据复合函数的链式求导法则,层状介质模型任意层j 的模型电阻率ρj 及厚度hj 的导数为:
(10) d ρ a d ρ j = d ρ a d x · d x d Z · d Z d ρ j d ρ a d h j = d ρ a d x · d x d Z · d Z d h j ,
(11) d ρ a d x = - 2 μ 0 a 2 4 · t · x 3 。
(12) d x d Z = 1 / ( d Z / d x ) = 1 / 3 x 3 erf ( x ) - 2 π e - x 2 2 + 3 x 2 ,
(13) d Z d ρ j = 4 a π · CT HT Im dfunc d ρ j / ω d Z d h j = 4 a π · CT HT Im dfunc d h j / ω 。
将式(11)~ 式(13)带入式(10),化简得到:
(14) d ρ a d ρ j = A · CT 1 ω HT Im dfunc d ρ j d ρ a d h j = A · CT 1 ω HT Im dfunc d h j 。
(15) A = 2 μ 0 a 3 π t · 1 2 π x ( 2 x 2 + 3 ) e - x 2 - 3 · erf ( x ) 。
(16) dfunc d ρ j = dfunc d u j ( 1 ) · d u j ( 1 ) d ρ j dfunc d h j = dfunc d u j ( 1 ) · d u j ( 1 ) d h j
(17) d u j ( 1 ) d ρ j = d u j ( 1 ) d u j ( 2 ) · d u j ( 2 ) d u j ( 3 ) · … · d u j ( j - 1 ) d u j ( j ) · d u j ( j ) d ρ j d u j ( 1 ) d h j = d u j ( 1 ) d u j ( 2 ) · d u j ( 2 ) d u j ( 3 ) · … · d u j ( j - 1 ) d u j ( j ) · d u j ( j ) d h j
(18) func U 1 j = dfunc d u j ( 1 ) = - λ 2 ( λ + u j ( 1 ) ) 2 ,
(19) UU 1 j = d u j ( 1 ) d u j ( 2 ) · d u j ( 2 ) d u j ( 3 ) · … · d u j ( j - 1 ) d u j ( j ) , ( j ≥ 2 ) 。
(20) d u j ( j ) d u j ( j + 1 ) = u j 2 [ 1 - tan h 2 ( u j h j ) ] [ u j + u ( j + 1 ) · tanh ( u j h j ) ] 2 ,
(21) d ρ a d ρ j = A · CT 1 ω HT Im func U 1 j · UU 1 j · d u j ( j ) d ρ j d ρ a d h j = A · CT 1 ω HT Im func U 1 j · UU 1 j · d u j ( j ) d h j ,
(22) d u j ( j ) d h j = u j 2 [ 1 - tan h 2 ( u j h j ) ] · [ u j 2 - [ u ( j + 1 ) ] 2 ] [ u j + u ( j + 1 ) tanh ( u j h j ) ] 2 ,
(23) d u j d ρ j = - i ω μ 0 2 · u j · ρ j 2 ,
(24) d u ( j ) j d ρ j = [ u j 2 + ( u ( j + 1 ) ) 2 ] · tanh ( u j h j ) + 2 u j · u ( j + 1 ) · tan h 2 ( u j h j ) + u j [ u j 2 - ( u ( j + 1 ) ) 2 ] · [ 1 - tan h 2 ( u j h j ) ] · h j · d u j d ρ j / [ u j + u ( j + 1 ) tanh ( u j h j ) ] 2 。
根据正演公式(式(1)、式(2)),进而推出第n 层的模型导数关系:
(25) d u n d ρ n = - i ω μ 0 2 u n ρ n 2 ,
(26) dfun c n d ρ n = dfun c n d u n · d u n d ρ n = ( - λ 2 ) ( λ + u n ) 2 · ( - i ω μ 0 ) 2 2 u n ρ n 2 。
3 Occam反演
Constable等提出Occam反演法[7 ] ,其基本原理是寻找在目标函数相对极小情况下的光滑模型,它具有迭代过程稳定、收敛较快且对初始模型依赖度较小等优点。该方法在反演过程中每迭代一次需要对模型求一次导数,因此由模型导数构成的雅克比矩阵的计算精度和计算效率将决定Occam反演的计算效果和计算效率。Occam反演原理简述如下。
(27) d j = F j [ m ] , j = 1,2 , … , M
(28) χ 2 = ‖ Wd - WF [ m ] ‖ 2 。
根据约束极小化理论,利用拉格朗日乘子形成一个无约束的函数U :
(29) U = ‖ ∂ m ‖ 2 + μ - 1 { ‖ Wd - WF [ m ] ‖ 2 } - χ * 2 ,
(30) ∂ = 0 0 - 1 1 - 1 1 … 0 - 1 1 N × N 。
式中右端第一项是粗糙度,第二项是由拉格朗日乘子加权的拟合差;取梯度,引起U 不变的向量m 遵循
(31) μ - 1 ( WJ ) T WJm - μ - 1 ( WJ ) T Wd + ∂ T ∂ m = 0 ,
(33) J = ∂ F ∂ m = ∂ F 1 ∂ m 1 … ∂ F 1 ∂ m N ∂ F M ∂ m 1 … ∂ F M ∂ m N 。
上式可通过式(21)、(22)、(26)计算得出。
最终,假设第k 次反演迭代已完成,k +1次迭代的模型向量为:
(34) m k + 1 ( μ ) = [ μ ∂ T ∂ + ( W J k ) T W J k ] - 1 ( W J k ) T W d k ,
(35) χ k + 1 ( μ ) = ‖ Wd - WF [ m k + 1 ( μ ) ] ‖ ,
最终找到一个μ 值相对最小,拟合差相对最小,光滑度相对高的模型,即是反演结果。
4 理论模型算例
(36) ∂ F i ( m ) ∂ m j = F i ( m j + Δ m j ) - F i ( m j - Δ m j ) 2 Δ m j
做对比。从理论上说,差分步长越短,结果越接近真实值,但步长过短,会由于正演的精度不够而偏差增大,一般采用的经验差分步长(Δm/m )取0.001。
对典型三层地电模型做一维Occam反演:不失一般性,Occam反演的初始模型选择均匀半空间,电阻率取响应的平均值,模型层采用对数等间距递增,起始厚度5 m,对数递增间距为0.1,共分10层,最大反演深度687.2 m。
H型地电模型:层电阻率分别为100、10、100 Ω·m,层厚度分别为60、30 m,计算结果见图1 、表1 。K型地电模型:层电阻率分别为10、100、10 Ω·m,层厚度分别为30、60 m,计算结果见图2 、表2 。
图1
图1
H型反演结果对比
Fig.1
Comparison of H-type layered Model inversion results
图2
图2
K型反演结果对比
Fig.2
Comparison of K-type layered Model inversion results
从上述算例可以看出,差分算法不仅计算效率低,而且曲线拟合不佳,反演模型与真实模型差距大,而解析算法则表现优秀。差分算法效果不佳的原因很可能与全区视电阻率的计算有关,因为全区视电阻率在计算隐函数Z(x)的反函数时,是通过二分查找后插值所得,本身误差较大,再对其进行差分求导,放大了误差。因为若差分步长较大,差分本身的误差较大;而步长较小,正演响应Hz 变化不大,二分查找的结果非常接近而趋于相等,差分结果趋近于零,总之会放大误差。随着反演不断迭代,误差会进一步累积放大,最终导致曲线拟合困难、反演模型震荡的结果(图1 、图2 )。因此,对于全区视电阻率的反演不建议使用差分算法。
5 实测数据实例
解析算法对理论模型的反演表现良好。为进一步验证解析算法Occam反演对实测数据的有效性,采用对实测剖面数据进行逐点一维反演的方法,即拟二维反演,绘制成电阻率反演断面图。与全区视电阻率拟断面、“烟圈”反演的结果进行效果对比,并通过地质钻探来验证方法的有效性。
5.1 地质概况
原始数据采集自内蒙古自治区凉城县西水塘镇附近,如图3 所示,为内蒙古凉城县某地热勘查项目中的一条剖面(400线),测线位于山前倾斜平原,布置该测线的目的是为划分工作区地层,推断基底埋深,并进一步验证地质推断的构造,为地热找水提供物探依据。工作区地形平坦,空旷且开阔,无工业电流干扰,非常有利于开展地面瞬变电磁工作。
图3
图3
工作区地质图
1—中低山; 2—山前倾斜平原; 3—第四系全新统冲洪积物; 4—侏罗系火山岩; 5—太古宇桑干群片麻岩; 6—太古宙似斑状花岗岩; 7—断裂破碎带; 8—测试瞬变电磁测深测线; 9—地质界线; 10—地貌界线; 11—已有钻孔;12—验证地热钻孔(DR2); 13—瞬变电磁发射线框位置; 14—瞬变电磁接收线圈位置
Fig.3
Geological map of work area
1—middle low mountain; 2—piedmont inclined plain; 3—Quaternary Holocene alluvial proluvial; 4—Jurassic volcanic rock; 5—Archaean Sanggan group gneiss; 6—Archaean porphyry like granite; 7—fracture fracture Fracture Zone; 8—test TEM sounding line; 9—geological boundary; 10—geomorphic boundary; 11—existing drilling;12—verification of geothermal drilling (DR2); 13—TEM emission line frame position; 14—TEM receiving coil position
根据区域地质资料显示[16 ] ,工作区地层主要有太古宇界桑干群,中生界侏罗系及新生界古近系、新近系、第四系(图3 )。区域地质构造属阴山东西向构造带的南缘,主要有NE和NW两个方向,其中NE向构造属高角度压性断层,断层连续性强,是地下热水和温泉形成的主要通道。侵入岩主要分布在蛮汉山区及平原区深部,主要为太古宙代早期侵入岩,岩性以苏长岩和似斑状花岗岩为主。
测线所在工作区全部为第四系覆盖。据区内已有钻孔显示(图3 ,图4 a),岩性从下至上分别为:太古宙似斑状花岗岩(γ),多为灰色中粗粒似斑状花岗岩;古近系(E)为深棕红色泥岩、泥质砂岩、砂质泥岩夹砂砾岩;新近系(N)上部为浅棕红色泥岩、泥质砂岩、砂质泥岩夹砂砾岩,下部为土黄色、灰白色砂砾岩,与古近系地层呈假整合接触;第四系(Q)岩性为浅黄棕色、黄褐色粉土、粉质黏土与砂砾卵石、含卵砾中粗砂、中粗砂等。
图4
图4
工作区钻孔柱状图
Fig.4
Borehole histogram of work area
5.2 Occam反演及地质解释
工作装置的设置:根据工作区内已知钻孔显示(图4 a),工作区地层分层明显,侵入岩、沉积岩具有明显的纵向电性差异,具备了瞬变电磁测深施工的物性条件。测线所在位置基底埋深推断为在200~300 m之间,且地层中包括较厚的一套泥岩地层,方法试验也证实了这套低阻泥岩的存在。 为穿透这套低阻地层探测到基底花岗岩,在充分利用场地宽度的情况下,设置发射装置中心回线为400 m×600 m。
本着从已知到未知的推断原则,设定Occam反演初始模型。初始模型采用深度对数等间隔模型,让反演模型与瞬变电磁响应的采样间隔方式一致,这样一方面有利于曲线拟合,另一方面可以尽可能多分层,提高反演的纵向分辨能力。经过多次实验,以兼顾反演计算效率和突出局部异常为原则,确定初始模型起始深度为20 m,避开瞬变电磁响应早期的盲区,对数间隔1.15,分25层,最小层厚3 m,最大层厚65 m,最大反演深度498 m;初始模型电阻率参考全区视电阻率的均值,定为60 Ω·m。
瞬变电磁测深视电阻率曲线大多表现为KH、KQH型,反演结果总体可分为5个大电性层(图5 )。第一层为中低阻,电阻率在 50 Ω·m左右,反演深度在0~30 m左右,与钻孔揭露的第四系粉土、黏土对应。第二层为中高阻层,电阻率在50~80 Ω·m之间,反演深度在30~80 m之间,对应以砂岩、砂砾岩、粉砂岩等为主的第四系地层,透水性好,为含水层。第三层为低阻层,电阻率缓慢下降,在10~50 Ω·m之间,反演深度80~150 m左右,与第四系地层粉质黏土对应,隔水性较好,为隔水层。第四层为低层,电阻率在10 Ω·m左右,反演深度150~300 m之间,与新近系,古近系的红色泥岩对应,为隔水层。第五层为高阻,电阻率大于100 Ω·m,对应基底太古代似斑状花岗岩。
图5
图5
钻孔投影点820点反演结果
Fig.5
Inversion result of site 820 where is projectionposition of borehole
上述地质推断,经钻孔DR2验证(图4 b),投影到剖面位置为820点,与实际情况基本吻合。该井定位于构造中,钻探深度 308.6 m,成井后自流,自流量为4 483.64 m3 /d,水温 36 ℃,矿化度 1.01 g/L,是内蒙古凉城县地热勘查中取得的突出成果[17 ] 。
5.3 反演结果对比分析
利用本实测数据地层分层明确、纵向电性差异明显的特点,可测试反演算法在纵向上的分辨能力,因此绘制了全区视电阻率拟断面和“烟圈”反演断面图,与Occam反演结果对比分析(图6 ):拟断面和“烟圈”反演总体表现为“高—低—高”三层模型,中间层表现为巨厚的低阻层,对地层的分辨能力有限,而Occam反演则反演出第四系次高阻层,对应含砂岩、粉砂岩较多的地层,该类地层渗透性好,一般为含水层,对找水工作意义重大;对于高阻基底,上述两种的方法的反演结果基本一致,与钻探结果基本吻合。
图6
图6
400线瞬变电磁反演电阻率断面对比
a—拟断面;b—“烟圈”反演结果;c—Occam反演结果
Fig.6
Comparison of TEM inversion sections of line 400
a—pseudo section; b—smoke ring inversion section; c—Occam inversion section
6 结论
采用基于瞬变电磁全区视电阻率的Occam反演,有效避免了晚期视电阻率早期变形严重而造成的假高阻。在理论正演模型算例中显示:雅克比矩阵的解析算法在计算效果和计算效率两方面都优于差分算法;在实测数据算例中可以看到:Occam反演结果与钻孔更为吻合,对地层的纵向分辨能力高于烟圈反演,特别是在第四系低阻地层中次高阻砂岩层的反映,对瞬变电磁找水工作具有较高的实用价值。上述算例证明,解析算法公式推导正确,基于全区视电阻率的Occam反演不仅反演效率高,而且效果好。
此外,本文不仅推导了Occam反演中雅克比矩阵所需的对模型层电阻率的求导公式,还推导出了对模型层厚度的求导公式,可应用于对模型层厚度反演的方法,如马奎特反演等。
致谢
感谢内蒙古国土资源勘查开发院以及内蒙古自治区第一水文地质工程地质勘查院提供的宝贵实测数据及相关地质资料;感谢审稿老师提出的宝贵意见。
参考文献
View Option
[1]
牛之琏 . 时间域电磁法原理 [M]. 长沙 : 中南大学出版社 , 2007 .
[本文引用: 1]
Niu Z L . Theory of time domain electromagnetic method [M]. Changsha : Press of Central South University , 2007 .
[本文引用: 1]
[2]
蒋邦远 . 实用近区磁源瞬变电磁法勘探 [M]. 北京 : 地质出版社 , 1998 .
[本文引用: 1]
Jiang B Y . Practical area near magnetic source transient electromagnetic method prospecting [M]. Beijing : Geology Publishing House , 1998 .
[本文引用: 1]
[3]
米萨克·N·纳比吉安 . 勘查地球物理-电磁法 [M].赵经祥译. 北京 : 地质出版社 , 1992 : 195 -209 .
[本文引用: 3]
Misac N N . Electromagnetic method in applied geophysics [M]. Beijing : Geological Publishing House , 1992 : 195 -209 .
[本文引用: 3]
[4]
Guptasarma D , Singh B . New digital linear filters for Hankel J0 and J1 transforms
[J]. Geophysical Prospecting , 1997 ,45 :745 -762 .
[本文引用: 1]
[5]
王华军 . 正余弦变换的数值滤波算法
[J]. 工程地球物理学报 , 2004 ,1 (4 ):329 -335 .
[本文引用: 1]
Wang H J . Digital filter algorithm of the sine and cosine transform
[J]. Chinese Journal of Engineering Geophysics , 2004 ,1 (4 ):329 -335 .
[本文引用: 1]
[6]
杨云见 , 何展翔 , 王绪本 . 直流电测深法与中心回线瞬变电磁法联合反演
[J]. 物探与化探 , 2008 ,32 (4 ):442 -444 .
[本文引用: 1]
Yang Y J , He Z X , Wang X B . A discussion of the inversion of combined vertical electrical sounding and central-loop TEM data
[J]. Geophysical and Geochemical Exploration , 2008 ,32 (4 ):442 -444 .
[本文引用: 1]
[7]
Constable S C , Parker R L , Catherine G . Constable Occam’ sinversion: apractical algorithm for generating smooth models from electromagnetic sounding data
[J]. Geophysics , 1987 ,52 (3 ):289 -300 .
[本文引用: 2]
[8]
李斌 , 赵虎 , 钟邱平 , 等 . 公路深埋隧道音频大地电磁法勘察中的反演方法简介及应用效果分析
[J]. 西南公路 , 2018 (2 ):6 -12 .
[本文引用: 1]
Li B , Zhao H , Zhong Q P , et al . Study on inversion and application of AMT in deeply-buried tunnel exploration
[J]. Outhwest Highway , 2018 (2 ):6 -12 .
[本文引用: 1]
[9]
冯思臣 , 王绪本 , 阮帅 . 一维大地电磁测深几种反演算法的比较研究
[J]. 石油地球物理勘探 , 2004 ,39 (5 ):594 -599 .
[本文引用: 1]
Feng S C , Wang X B , Ruan S . Comparison among several inversion algorithm of 1D MT
[J]. Oil geophysical prospecting , 2004 ,39 (5 ):594 -599 .
[本文引用: 1]
[10]
郭睿 , 武欣 , 肖攀 , 等 . 直升机航空瞬变电磁Occam反演及软件实现
[J]. 电子测量技术 , 2017 ,40 (9 ):169 -173 .
[本文引用: 1]
Guo R , Wu X , Xiao P , et al . Helicopter transient electromagnetic occam inversion and software development
[J]. Electronic measurement technology , 2017 ,40 (9 ):169 -173 .
[本文引用: 1]
[11]
陈小斌 , 赵国泽 , 汤吉 , 等 . 大地电磁自适应正则化反演算法
[J]. 地球物理学报 , 2005 ,48 (4 ):937 -946 .
[本文引用: 1]
Chen X B , Zhao G Z , Tang J , et al . An adaptive regularized inversion algorithm for magnetotelluric data
[J]. Chinese Journal of Geophysics , 2005 ,48 (4 ):937 -946 .
[本文引用: 1]
[12]
毛立峰 , 王绪本 , 陈斌 . 直升机航空瞬变电磁自适应正则化一维反演方法研究
[J]. 地球物理学进展 , 2011 ,26 (1 ):300 -305 .
[本文引用: 1]
Mao L F , Wang X B , Chen B . Study on an adaptive regularized 1D inversion method of helicopter TEM data
[J]. Progress in Geophysics , 2011 ,26 (1 ):300 -305 .
[本文引用: 1]
[13]
李刚 , 潘和平 , 王智 , 等 . 回线源瞬变电磁法一维反演算法
[J]. 煤田地质与勘探 , 2017 ,45 (5 ):161 -166,172 .
[本文引用: 1]
Li G , Pan H P , Wang Z , et al . One-dimensional inversion for loop source transient electromagnetic method
[J]. Coal Geology & Exploration , 2017 ,45 (5 ):161 -166,172 .
[本文引用: 1]
[14]
白登海 , Maxwell A Meju , 卢健 , 等 . 时间域瞬变电磁法中心回线方式全程视电阻率的数值计算
[J]. 地球物理学报 , 2003 ,46 (5 ):698 -704 .
[本文引用: 1]
Bai D H , Maxwell A Meju , Lu J , et al . Numerical calculation of all-time apparent resistivity for the central loop transient electromagnetic method
[J]. Chinese journal of geophysics , 2003 ,46 (5 ):698 -704 .
[本文引用: 1]
[15]
郭嵩巍 , 王绪本 . 瞬变电磁全区视电阻率数值计算方法研究
[J]. 物探化探计算技术 , 2010 ,32 (5 ):500 -507,455 .
[本文引用: 2]
Guo S W , Wang X B . Study on numerical calculation methods of transient electromagnetic all-time apparent resistivity
[J]. Computing Technique for Geophysical and Geochemical Exploration , 2010 ,32 (5 ):500 -507,455 .
[本文引用: 2]
[16]
内蒙古地质矿产局 . 内蒙古自治区岩石地层 [M]. 武汉 : 中国地质大学出版社 , 1996 .
[本文引用: 1]
Inner Mongolia Bureau of Geology and Mineral Resources . Lithostratigraphic of Inner Mongolia . Wuhan : China University of Geosciences Press , 1996 .
[本文引用: 1]
[17]
任阿龙 . 凉城县岱海地区地热资源预可行性勘查项目取得突出成果
[J]. 西部资源 , 2014 (3 ):32 .
[本文引用: 1]
Ren A L . Outstanding achievements in the pre Feasibility Exploration Project of geothermal resources in Daihai area, Liangcheng County
[J]. Westen Resource , 2014 (3 ):32 .
[本文引用: 1]
1
2007
... 瞬变电磁法作为一种常见的电磁测深物探方法,广泛应用于水文勘查、工程勘查以及矿产勘查[1 ,2 ] .在目前的工程应用中,除“烟圈”反演快速成像法外,一维反演结果仍然是地质解释的主要依据.对于中心回线装置,可以解析推导出一维正演响应的计算公式[3 ] ,其数值计算的核心是,内层积分的汉克尔变换[4 ] 、外层积分的余弦变换[5 ] 均可采用数值滤波算法,随着滤波系数的不断推陈出新,计算量明显减少,大大提高了正演的计算效率和精度.在反演方面,目前常见的反演方法,如阻尼最小二乘法(马奎特法)[6 ] 、Occam法[7 ,8 ,9 ,10 ] 、正则化法[11 ,12 ,13 ] 等,都需要计算正演响应对模型的导数,即雅克比矩阵,因此雅克比矩阵的计算效率和效果就成为反演计算的关键因素之一.计算雅克比矩阵最简便的方法是差分法,该算法程序实现较为简单,缺点是不仅调用正演次数多,计算效率低,而且有可能增加累积误差.因此,本文详细推导了中心回线装置瞬变电磁一维正演全区视电阻率的雅克比矩阵的解析计算公式,并通过Occam反演,对比分析了理论模型和实测数据的反演效果,证明推导正确,值得采用. ...
1
2007
... 瞬变电磁法作为一种常见的电磁测深物探方法,广泛应用于水文勘查、工程勘查以及矿产勘查[1 ,2 ] .在目前的工程应用中,除“烟圈”反演快速成像法外,一维反演结果仍然是地质解释的主要依据.对于中心回线装置,可以解析推导出一维正演响应的计算公式[3 ] ,其数值计算的核心是,内层积分的汉克尔变换[4 ] 、外层积分的余弦变换[5 ] 均可采用数值滤波算法,随着滤波系数的不断推陈出新,计算量明显减少,大大提高了正演的计算效率和精度.在反演方面,目前常见的反演方法,如阻尼最小二乘法(马奎特法)[6 ] 、Occam法[7 ,8 ,9 ,10 ] 、正则化法[11 ,12 ,13 ] 等,都需要计算正演响应对模型的导数,即雅克比矩阵,因此雅克比矩阵的计算效率和效果就成为反演计算的关键因素之一.计算雅克比矩阵最简便的方法是差分法,该算法程序实现较为简单,缺点是不仅调用正演次数多,计算效率低,而且有可能增加累积误差.因此,本文详细推导了中心回线装置瞬变电磁一维正演全区视电阻率的雅克比矩阵的解析计算公式,并通过Occam反演,对比分析了理论模型和实测数据的反演效果,证明推导正确,值得采用. ...
1
1998
... 瞬变电磁法作为一种常见的电磁测深物探方法,广泛应用于水文勘查、工程勘查以及矿产勘查[1 ,2 ] .在目前的工程应用中,除“烟圈”反演快速成像法外,一维反演结果仍然是地质解释的主要依据.对于中心回线装置,可以解析推导出一维正演响应的计算公式[3 ] ,其数值计算的核心是,内层积分的汉克尔变换[4 ] 、外层积分的余弦变换[5 ] 均可采用数值滤波算法,随着滤波系数的不断推陈出新,计算量明显减少,大大提高了正演的计算效率和精度.在反演方面,目前常见的反演方法,如阻尼最小二乘法(马奎特法)[6 ] 、Occam法[7 ,8 ,9 ,10 ] 、正则化法[11 ,12 ,13 ] 等,都需要计算正演响应对模型的导数,即雅克比矩阵,因此雅克比矩阵的计算效率和效果就成为反演计算的关键因素之一.计算雅克比矩阵最简便的方法是差分法,该算法程序实现较为简单,缺点是不仅调用正演次数多,计算效率低,而且有可能增加累积误差.因此,本文详细推导了中心回线装置瞬变电磁一维正演全区视电阻率的雅克比矩阵的解析计算公式,并通过Occam反演,对比分析了理论模型和实测数据的反演效果,证明推导正确,值得采用. ...
1
1998
... 瞬变电磁法作为一种常见的电磁测深物探方法,广泛应用于水文勘查、工程勘查以及矿产勘查[1 ,2 ] .在目前的工程应用中,除“烟圈”反演快速成像法外,一维反演结果仍然是地质解释的主要依据.对于中心回线装置,可以解析推导出一维正演响应的计算公式[3 ] ,其数值计算的核心是,内层积分的汉克尔变换[4 ] 、外层积分的余弦变换[5 ] 均可采用数值滤波算法,随着滤波系数的不断推陈出新,计算量明显减少,大大提高了正演的计算效率和精度.在反演方面,目前常见的反演方法,如阻尼最小二乘法(马奎特法)[6 ] 、Occam法[7 ,8 ,9 ,10 ] 、正则化法[11 ,12 ,13 ] 等,都需要计算正演响应对模型的导数,即雅克比矩阵,因此雅克比矩阵的计算效率和效果就成为反演计算的关键因素之一.计算雅克比矩阵最简便的方法是差分法,该算法程序实现较为简单,缺点是不仅调用正演次数多,计算效率低,而且有可能增加累积误差.因此,本文详细推导了中心回线装置瞬变电磁一维正演全区视电阻率的雅克比矩阵的解析计算公式,并通过Occam反演,对比分析了理论模型和实测数据的反演效果,证明推导正确,值得采用. ...
3
1992
... 瞬变电磁法作为一种常见的电磁测深物探方法,广泛应用于水文勘查、工程勘查以及矿产勘查[1 ,2 ] .在目前的工程应用中,除“烟圈”反演快速成像法外,一维反演结果仍然是地质解释的主要依据.对于中心回线装置,可以解析推导出一维正演响应的计算公式[3 ] ,其数值计算的核心是,内层积分的汉克尔变换[4 ] 、外层积分的余弦变换[5 ] 均可采用数值滤波算法,随着滤波系数的不断推陈出新,计算量明显减少,大大提高了正演的计算效率和精度.在反演方面,目前常见的反演方法,如阻尼最小二乘法(马奎特法)[6 ] 、Occam法[7 ,8 ,9 ,10 ] 、正则化法[11 ,12 ,13 ] 等,都需要计算正演响应对模型的导数,即雅克比矩阵,因此雅克比矩阵的计算效率和效果就成为反演计算的关键因素之一.计算雅克比矩阵最简便的方法是差分法,该算法程序实现较为简单,缺点是不仅调用正演次数多,计算效率低,而且有可能增加累积误差.因此,本文详细推导了中心回线装置瞬变电磁一维正演全区视电阻率的雅克比矩阵的解析计算公式,并通过Occam反演,对比分析了理论模型和实测数据的反演效果,证明推导正确,值得采用. ...
... 其中:a 为发射线框等效半径;u 的递推公式[3 ] 为 ...
... 对于均匀半空间地电模型,可以推导出Hz 的解析表达式[3 ] : ...
3
1992
... 瞬变电磁法作为一种常见的电磁测深物探方法,广泛应用于水文勘查、工程勘查以及矿产勘查[1 ,2 ] .在目前的工程应用中,除“烟圈”反演快速成像法外,一维反演结果仍然是地质解释的主要依据.对于中心回线装置,可以解析推导出一维正演响应的计算公式[3 ] ,其数值计算的核心是,内层积分的汉克尔变换[4 ] 、外层积分的余弦变换[5 ] 均可采用数值滤波算法,随着滤波系数的不断推陈出新,计算量明显减少,大大提高了正演的计算效率和精度.在反演方面,目前常见的反演方法,如阻尼最小二乘法(马奎特法)[6 ] 、Occam法[7 ,8 ,9 ,10 ] 、正则化法[11 ,12 ,13 ] 等,都需要计算正演响应对模型的导数,即雅克比矩阵,因此雅克比矩阵的计算效率和效果就成为反演计算的关键因素之一.计算雅克比矩阵最简便的方法是差分法,该算法程序实现较为简单,缺点是不仅调用正演次数多,计算效率低,而且有可能增加累积误差.因此,本文详细推导了中心回线装置瞬变电磁一维正演全区视电阻率的雅克比矩阵的解析计算公式,并通过Occam反演,对比分析了理论模型和实测数据的反演效果,证明推导正确,值得采用. ...
... 其中:a 为发射线框等效半径;u 的递推公式[3 ] 为 ...
... 对于均匀半空间地电模型,可以推导出Hz 的解析表达式[3 ] : ...
New digital linear filters for Hankel J0 and J1 transforms
1
1997
... 瞬变电磁法作为一种常见的电磁测深物探方法,广泛应用于水文勘查、工程勘查以及矿产勘查[1 ,2 ] .在目前的工程应用中,除“烟圈”反演快速成像法外,一维反演结果仍然是地质解释的主要依据.对于中心回线装置,可以解析推导出一维正演响应的计算公式[3 ] ,其数值计算的核心是,内层积分的汉克尔变换[4 ] 、外层积分的余弦变换[5 ] 均可采用数值滤波算法,随着滤波系数的不断推陈出新,计算量明显减少,大大提高了正演的计算效率和精度.在反演方面,目前常见的反演方法,如阻尼最小二乘法(马奎特法)[6 ] 、Occam法[7 ,8 ,9 ,10 ] 、正则化法[11 ,12 ,13 ] 等,都需要计算正演响应对模型的导数,即雅克比矩阵,因此雅克比矩阵的计算效率和效果就成为反演计算的关键因素之一.计算雅克比矩阵最简便的方法是差分法,该算法程序实现较为简单,缺点是不仅调用正演次数多,计算效率低,而且有可能增加累积误差.因此,本文详细推导了中心回线装置瞬变电磁一维正演全区视电阻率的雅克比矩阵的解析计算公式,并通过Occam反演,对比分析了理论模型和实测数据的反演效果,证明推导正确,值得采用. ...
正余弦变换的数值滤波算法
1
2004
... 瞬变电磁法作为一种常见的电磁测深物探方法,广泛应用于水文勘查、工程勘查以及矿产勘查[1 ,2 ] .在目前的工程应用中,除“烟圈”反演快速成像法外,一维反演结果仍然是地质解释的主要依据.对于中心回线装置,可以解析推导出一维正演响应的计算公式[3 ] ,其数值计算的核心是,内层积分的汉克尔变换[4 ] 、外层积分的余弦变换[5 ] 均可采用数值滤波算法,随着滤波系数的不断推陈出新,计算量明显减少,大大提高了正演的计算效率和精度.在反演方面,目前常见的反演方法,如阻尼最小二乘法(马奎特法)[6 ] 、Occam法[7 ,8 ,9 ,10 ] 、正则化法[11 ,12 ,13 ] 等,都需要计算正演响应对模型的导数,即雅克比矩阵,因此雅克比矩阵的计算效率和效果就成为反演计算的关键因素之一.计算雅克比矩阵最简便的方法是差分法,该算法程序实现较为简单,缺点是不仅调用正演次数多,计算效率低,而且有可能增加累积误差.因此,本文详细推导了中心回线装置瞬变电磁一维正演全区视电阻率的雅克比矩阵的解析计算公式,并通过Occam反演,对比分析了理论模型和实测数据的反演效果,证明推导正确,值得采用. ...
正余弦变换的数值滤波算法
1
2004
... 瞬变电磁法作为一种常见的电磁测深物探方法,广泛应用于水文勘查、工程勘查以及矿产勘查[1 ,2 ] .在目前的工程应用中,除“烟圈”反演快速成像法外,一维反演结果仍然是地质解释的主要依据.对于中心回线装置,可以解析推导出一维正演响应的计算公式[3 ] ,其数值计算的核心是,内层积分的汉克尔变换[4 ] 、外层积分的余弦变换[5 ] 均可采用数值滤波算法,随着滤波系数的不断推陈出新,计算量明显减少,大大提高了正演的计算效率和精度.在反演方面,目前常见的反演方法,如阻尼最小二乘法(马奎特法)[6 ] 、Occam法[7 ,8 ,9 ,10 ] 、正则化法[11 ,12 ,13 ] 等,都需要计算正演响应对模型的导数,即雅克比矩阵,因此雅克比矩阵的计算效率和效果就成为反演计算的关键因素之一.计算雅克比矩阵最简便的方法是差分法,该算法程序实现较为简单,缺点是不仅调用正演次数多,计算效率低,而且有可能增加累积误差.因此,本文详细推导了中心回线装置瞬变电磁一维正演全区视电阻率的雅克比矩阵的解析计算公式,并通过Occam反演,对比分析了理论模型和实测数据的反演效果,证明推导正确,值得采用. ...
直流电测深法与中心回线瞬变电磁法联合反演
1
2008
... 瞬变电磁法作为一种常见的电磁测深物探方法,广泛应用于水文勘查、工程勘查以及矿产勘查[1 ,2 ] .在目前的工程应用中,除“烟圈”反演快速成像法外,一维反演结果仍然是地质解释的主要依据.对于中心回线装置,可以解析推导出一维正演响应的计算公式[3 ] ,其数值计算的核心是,内层积分的汉克尔变换[4 ] 、外层积分的余弦变换[5 ] 均可采用数值滤波算法,随着滤波系数的不断推陈出新,计算量明显减少,大大提高了正演的计算效率和精度.在反演方面,目前常见的反演方法,如阻尼最小二乘法(马奎特法)[6 ] 、Occam法[7 ,8 ,9 ,10 ] 、正则化法[11 ,12 ,13 ] 等,都需要计算正演响应对模型的导数,即雅克比矩阵,因此雅克比矩阵的计算效率和效果就成为反演计算的关键因素之一.计算雅克比矩阵最简便的方法是差分法,该算法程序实现较为简单,缺点是不仅调用正演次数多,计算效率低,而且有可能增加累积误差.因此,本文详细推导了中心回线装置瞬变电磁一维正演全区视电阻率的雅克比矩阵的解析计算公式,并通过Occam反演,对比分析了理论模型和实测数据的反演效果,证明推导正确,值得采用. ...
直流电测深法与中心回线瞬变电磁法联合反演
1
2008
... 瞬变电磁法作为一种常见的电磁测深物探方法,广泛应用于水文勘查、工程勘查以及矿产勘查[1 ,2 ] .在目前的工程应用中,除“烟圈”反演快速成像法外,一维反演结果仍然是地质解释的主要依据.对于中心回线装置,可以解析推导出一维正演响应的计算公式[3 ] ,其数值计算的核心是,内层积分的汉克尔变换[4 ] 、外层积分的余弦变换[5 ] 均可采用数值滤波算法,随着滤波系数的不断推陈出新,计算量明显减少,大大提高了正演的计算效率和精度.在反演方面,目前常见的反演方法,如阻尼最小二乘法(马奎特法)[6 ] 、Occam法[7 ,8 ,9 ,10 ] 、正则化法[11 ,12 ,13 ] 等,都需要计算正演响应对模型的导数,即雅克比矩阵,因此雅克比矩阵的计算效率和效果就成为反演计算的关键因素之一.计算雅克比矩阵最简便的方法是差分法,该算法程序实现较为简单,缺点是不仅调用正演次数多,计算效率低,而且有可能增加累积误差.因此,本文详细推导了中心回线装置瞬变电磁一维正演全区视电阻率的雅克比矩阵的解析计算公式,并通过Occam反演,对比分析了理论模型和实测数据的反演效果,证明推导正确,值得采用. ...
Constable Occam’ sinversion: apractical algorithm for generating smooth models from electromagnetic sounding data
2
1987
... 瞬变电磁法作为一种常见的电磁测深物探方法,广泛应用于水文勘查、工程勘查以及矿产勘查[1 ,2 ] .在目前的工程应用中,除“烟圈”反演快速成像法外,一维反演结果仍然是地质解释的主要依据.对于中心回线装置,可以解析推导出一维正演响应的计算公式[3 ] ,其数值计算的核心是,内层积分的汉克尔变换[4 ] 、外层积分的余弦变换[5 ] 均可采用数值滤波算法,随着滤波系数的不断推陈出新,计算量明显减少,大大提高了正演的计算效率和精度.在反演方面,目前常见的反演方法,如阻尼最小二乘法(马奎特法)[6 ] 、Occam法[7 ,8 ,9 ,10 ] 、正则化法[11 ,12 ,13 ] 等,都需要计算正演响应对模型的导数,即雅克比矩阵,因此雅克比矩阵的计算效率和效果就成为反演计算的关键因素之一.计算雅克比矩阵最简便的方法是差分法,该算法程序实现较为简单,缺点是不仅调用正演次数多,计算效率低,而且有可能增加累积误差.因此,本文详细推导了中心回线装置瞬变电磁一维正演全区视电阻率的雅克比矩阵的解析计算公式,并通过Occam反演,对比分析了理论模型和实测数据的反演效果,证明推导正确,值得采用. ...
... Constable等提出Occam反演法[7 ] ,其基本原理是寻找在目标函数相对极小情况下的光滑模型,它具有迭代过程稳定、收敛较快且对初始模型依赖度较小等优点.该方法在反演过程中每迭代一次需要对模型求一次导数,因此由模型导数构成的雅克比矩阵的计算精度和计算效率将决定Occam反演的计算效果和计算效率.Occam反演原理简述如下. ...
公路深埋隧道音频大地电磁法勘察中的反演方法简介及应用效果分析
1
2018
... 瞬变电磁法作为一种常见的电磁测深物探方法,广泛应用于水文勘查、工程勘查以及矿产勘查[1 ,2 ] .在目前的工程应用中,除“烟圈”反演快速成像法外,一维反演结果仍然是地质解释的主要依据.对于中心回线装置,可以解析推导出一维正演响应的计算公式[3 ] ,其数值计算的核心是,内层积分的汉克尔变换[4 ] 、外层积分的余弦变换[5 ] 均可采用数值滤波算法,随着滤波系数的不断推陈出新,计算量明显减少,大大提高了正演的计算效率和精度.在反演方面,目前常见的反演方法,如阻尼最小二乘法(马奎特法)[6 ] 、Occam法[7 ,8 ,9 ,10 ] 、正则化法[11 ,12 ,13 ] 等,都需要计算正演响应对模型的导数,即雅克比矩阵,因此雅克比矩阵的计算效率和效果就成为反演计算的关键因素之一.计算雅克比矩阵最简便的方法是差分法,该算法程序实现较为简单,缺点是不仅调用正演次数多,计算效率低,而且有可能增加累积误差.因此,本文详细推导了中心回线装置瞬变电磁一维正演全区视电阻率的雅克比矩阵的解析计算公式,并通过Occam反演,对比分析了理论模型和实测数据的反演效果,证明推导正确,值得采用. ...
公路深埋隧道音频大地电磁法勘察中的反演方法简介及应用效果分析
1
2018
... 瞬变电磁法作为一种常见的电磁测深物探方法,广泛应用于水文勘查、工程勘查以及矿产勘查[1 ,2 ] .在目前的工程应用中,除“烟圈”反演快速成像法外,一维反演结果仍然是地质解释的主要依据.对于中心回线装置,可以解析推导出一维正演响应的计算公式[3 ] ,其数值计算的核心是,内层积分的汉克尔变换[4 ] 、外层积分的余弦变换[5 ] 均可采用数值滤波算法,随着滤波系数的不断推陈出新,计算量明显减少,大大提高了正演的计算效率和精度.在反演方面,目前常见的反演方法,如阻尼最小二乘法(马奎特法)[6 ] 、Occam法[7 ,8 ,9 ,10 ] 、正则化法[11 ,12 ,13 ] 等,都需要计算正演响应对模型的导数,即雅克比矩阵,因此雅克比矩阵的计算效率和效果就成为反演计算的关键因素之一.计算雅克比矩阵最简便的方法是差分法,该算法程序实现较为简单,缺点是不仅调用正演次数多,计算效率低,而且有可能增加累积误差.因此,本文详细推导了中心回线装置瞬变电磁一维正演全区视电阻率的雅克比矩阵的解析计算公式,并通过Occam反演,对比分析了理论模型和实测数据的反演效果,证明推导正确,值得采用. ...
一维大地电磁测深几种反演算法的比较研究
1
2004
... 瞬变电磁法作为一种常见的电磁测深物探方法,广泛应用于水文勘查、工程勘查以及矿产勘查[1 ,2 ] .在目前的工程应用中,除“烟圈”反演快速成像法外,一维反演结果仍然是地质解释的主要依据.对于中心回线装置,可以解析推导出一维正演响应的计算公式[3 ] ,其数值计算的核心是,内层积分的汉克尔变换[4 ] 、外层积分的余弦变换[5 ] 均可采用数值滤波算法,随着滤波系数的不断推陈出新,计算量明显减少,大大提高了正演的计算效率和精度.在反演方面,目前常见的反演方法,如阻尼最小二乘法(马奎特法)[6 ] 、Occam法[7 ,8 ,9 ,10 ] 、正则化法[11 ,12 ,13 ] 等,都需要计算正演响应对模型的导数,即雅克比矩阵,因此雅克比矩阵的计算效率和效果就成为反演计算的关键因素之一.计算雅克比矩阵最简便的方法是差分法,该算法程序实现较为简单,缺点是不仅调用正演次数多,计算效率低,而且有可能增加累积误差.因此,本文详细推导了中心回线装置瞬变电磁一维正演全区视电阻率的雅克比矩阵的解析计算公式,并通过Occam反演,对比分析了理论模型和实测数据的反演效果,证明推导正确,值得采用. ...
一维大地电磁测深几种反演算法的比较研究
1
2004
... 瞬变电磁法作为一种常见的电磁测深物探方法,广泛应用于水文勘查、工程勘查以及矿产勘查[1 ,2 ] .在目前的工程应用中,除“烟圈”反演快速成像法外,一维反演结果仍然是地质解释的主要依据.对于中心回线装置,可以解析推导出一维正演响应的计算公式[3 ] ,其数值计算的核心是,内层积分的汉克尔变换[4 ] 、外层积分的余弦变换[5 ] 均可采用数值滤波算法,随着滤波系数的不断推陈出新,计算量明显减少,大大提高了正演的计算效率和精度.在反演方面,目前常见的反演方法,如阻尼最小二乘法(马奎特法)[6 ] 、Occam法[7 ,8 ,9 ,10 ] 、正则化法[11 ,12 ,13 ] 等,都需要计算正演响应对模型的导数,即雅克比矩阵,因此雅克比矩阵的计算效率和效果就成为反演计算的关键因素之一.计算雅克比矩阵最简便的方法是差分法,该算法程序实现较为简单,缺点是不仅调用正演次数多,计算效率低,而且有可能增加累积误差.因此,本文详细推导了中心回线装置瞬变电磁一维正演全区视电阻率的雅克比矩阵的解析计算公式,并通过Occam反演,对比分析了理论模型和实测数据的反演效果,证明推导正确,值得采用. ...
直升机航空瞬变电磁Occam反演及软件实现
1
2017
... 瞬变电磁法作为一种常见的电磁测深物探方法,广泛应用于水文勘查、工程勘查以及矿产勘查[1 ,2 ] .在目前的工程应用中,除“烟圈”反演快速成像法外,一维反演结果仍然是地质解释的主要依据.对于中心回线装置,可以解析推导出一维正演响应的计算公式[3 ] ,其数值计算的核心是,内层积分的汉克尔变换[4 ] 、外层积分的余弦变换[5 ] 均可采用数值滤波算法,随着滤波系数的不断推陈出新,计算量明显减少,大大提高了正演的计算效率和精度.在反演方面,目前常见的反演方法,如阻尼最小二乘法(马奎特法)[6 ] 、Occam法[7 ,8 ,9 ,10 ] 、正则化法[11 ,12 ,13 ] 等,都需要计算正演响应对模型的导数,即雅克比矩阵,因此雅克比矩阵的计算效率和效果就成为反演计算的关键因素之一.计算雅克比矩阵最简便的方法是差分法,该算法程序实现较为简单,缺点是不仅调用正演次数多,计算效率低,而且有可能增加累积误差.因此,本文详细推导了中心回线装置瞬变电磁一维正演全区视电阻率的雅克比矩阵的解析计算公式,并通过Occam反演,对比分析了理论模型和实测数据的反演效果,证明推导正确,值得采用. ...
直升机航空瞬变电磁Occam反演及软件实现
1
2017
... 瞬变电磁法作为一种常见的电磁测深物探方法,广泛应用于水文勘查、工程勘查以及矿产勘查[1 ,2 ] .在目前的工程应用中,除“烟圈”反演快速成像法外,一维反演结果仍然是地质解释的主要依据.对于中心回线装置,可以解析推导出一维正演响应的计算公式[3 ] ,其数值计算的核心是,内层积分的汉克尔变换[4 ] 、外层积分的余弦变换[5 ] 均可采用数值滤波算法,随着滤波系数的不断推陈出新,计算量明显减少,大大提高了正演的计算效率和精度.在反演方面,目前常见的反演方法,如阻尼最小二乘法(马奎特法)[6 ] 、Occam法[7 ,8 ,9 ,10 ] 、正则化法[11 ,12 ,13 ] 等,都需要计算正演响应对模型的导数,即雅克比矩阵,因此雅克比矩阵的计算效率和效果就成为反演计算的关键因素之一.计算雅克比矩阵最简便的方法是差分法,该算法程序实现较为简单,缺点是不仅调用正演次数多,计算效率低,而且有可能增加累积误差.因此,本文详细推导了中心回线装置瞬变电磁一维正演全区视电阻率的雅克比矩阵的解析计算公式,并通过Occam反演,对比分析了理论模型和实测数据的反演效果,证明推导正确,值得采用. ...
大地电磁自适应正则化反演算法
1
2005
... 瞬变电磁法作为一种常见的电磁测深物探方法,广泛应用于水文勘查、工程勘查以及矿产勘查[1 ,2 ] .在目前的工程应用中,除“烟圈”反演快速成像法外,一维反演结果仍然是地质解释的主要依据.对于中心回线装置,可以解析推导出一维正演响应的计算公式[3 ] ,其数值计算的核心是,内层积分的汉克尔变换[4 ] 、外层积分的余弦变换[5 ] 均可采用数值滤波算法,随着滤波系数的不断推陈出新,计算量明显减少,大大提高了正演的计算效率和精度.在反演方面,目前常见的反演方法,如阻尼最小二乘法(马奎特法)[6 ] 、Occam法[7 ,8 ,9 ,10 ] 、正则化法[11 ,12 ,13 ] 等,都需要计算正演响应对模型的导数,即雅克比矩阵,因此雅克比矩阵的计算效率和效果就成为反演计算的关键因素之一.计算雅克比矩阵最简便的方法是差分法,该算法程序实现较为简单,缺点是不仅调用正演次数多,计算效率低,而且有可能增加累积误差.因此,本文详细推导了中心回线装置瞬变电磁一维正演全区视电阻率的雅克比矩阵的解析计算公式,并通过Occam反演,对比分析了理论模型和实测数据的反演效果,证明推导正确,值得采用. ...
大地电磁自适应正则化反演算法
1
2005
... 瞬变电磁法作为一种常见的电磁测深物探方法,广泛应用于水文勘查、工程勘查以及矿产勘查[1 ,2 ] .在目前的工程应用中,除“烟圈”反演快速成像法外,一维反演结果仍然是地质解释的主要依据.对于中心回线装置,可以解析推导出一维正演响应的计算公式[3 ] ,其数值计算的核心是,内层积分的汉克尔变换[4 ] 、外层积分的余弦变换[5 ] 均可采用数值滤波算法,随着滤波系数的不断推陈出新,计算量明显减少,大大提高了正演的计算效率和精度.在反演方面,目前常见的反演方法,如阻尼最小二乘法(马奎特法)[6 ] 、Occam法[7 ,8 ,9 ,10 ] 、正则化法[11 ,12 ,13 ] 等,都需要计算正演响应对模型的导数,即雅克比矩阵,因此雅克比矩阵的计算效率和效果就成为反演计算的关键因素之一.计算雅克比矩阵最简便的方法是差分法,该算法程序实现较为简单,缺点是不仅调用正演次数多,计算效率低,而且有可能增加累积误差.因此,本文详细推导了中心回线装置瞬变电磁一维正演全区视电阻率的雅克比矩阵的解析计算公式,并通过Occam反演,对比分析了理论模型和实测数据的反演效果,证明推导正确,值得采用. ...
直升机航空瞬变电磁自适应正则化一维反演方法研究
1
2011
... 瞬变电磁法作为一种常见的电磁测深物探方法,广泛应用于水文勘查、工程勘查以及矿产勘查[1 ,2 ] .在目前的工程应用中,除“烟圈”反演快速成像法外,一维反演结果仍然是地质解释的主要依据.对于中心回线装置,可以解析推导出一维正演响应的计算公式[3 ] ,其数值计算的核心是,内层积分的汉克尔变换[4 ] 、外层积分的余弦变换[5 ] 均可采用数值滤波算法,随着滤波系数的不断推陈出新,计算量明显减少,大大提高了正演的计算效率和精度.在反演方面,目前常见的反演方法,如阻尼最小二乘法(马奎特法)[6 ] 、Occam法[7 ,8 ,9 ,10 ] 、正则化法[11 ,12 ,13 ] 等,都需要计算正演响应对模型的导数,即雅克比矩阵,因此雅克比矩阵的计算效率和效果就成为反演计算的关键因素之一.计算雅克比矩阵最简便的方法是差分法,该算法程序实现较为简单,缺点是不仅调用正演次数多,计算效率低,而且有可能增加累积误差.因此,本文详细推导了中心回线装置瞬变电磁一维正演全区视电阻率的雅克比矩阵的解析计算公式,并通过Occam反演,对比分析了理论模型和实测数据的反演效果,证明推导正确,值得采用. ...
直升机航空瞬变电磁自适应正则化一维反演方法研究
1
2011
... 瞬变电磁法作为一种常见的电磁测深物探方法,广泛应用于水文勘查、工程勘查以及矿产勘查[1 ,2 ] .在目前的工程应用中,除“烟圈”反演快速成像法外,一维反演结果仍然是地质解释的主要依据.对于中心回线装置,可以解析推导出一维正演响应的计算公式[3 ] ,其数值计算的核心是,内层积分的汉克尔变换[4 ] 、外层积分的余弦变换[5 ] 均可采用数值滤波算法,随着滤波系数的不断推陈出新,计算量明显减少,大大提高了正演的计算效率和精度.在反演方面,目前常见的反演方法,如阻尼最小二乘法(马奎特法)[6 ] 、Occam法[7 ,8 ,9 ,10 ] 、正则化法[11 ,12 ,13 ] 等,都需要计算正演响应对模型的导数,即雅克比矩阵,因此雅克比矩阵的计算效率和效果就成为反演计算的关键因素之一.计算雅克比矩阵最简便的方法是差分法,该算法程序实现较为简单,缺点是不仅调用正演次数多,计算效率低,而且有可能增加累积误差.因此,本文详细推导了中心回线装置瞬变电磁一维正演全区视电阻率的雅克比矩阵的解析计算公式,并通过Occam反演,对比分析了理论模型和实测数据的反演效果,证明推导正确,值得采用. ...
回线源瞬变电磁法一维反演算法
1
2017
... 瞬变电磁法作为一种常见的电磁测深物探方法,广泛应用于水文勘查、工程勘查以及矿产勘查[1 ,2 ] .在目前的工程应用中,除“烟圈”反演快速成像法外,一维反演结果仍然是地质解释的主要依据.对于中心回线装置,可以解析推导出一维正演响应的计算公式[3 ] ,其数值计算的核心是,内层积分的汉克尔变换[4 ] 、外层积分的余弦变换[5 ] 均可采用数值滤波算法,随着滤波系数的不断推陈出新,计算量明显减少,大大提高了正演的计算效率和精度.在反演方面,目前常见的反演方法,如阻尼最小二乘法(马奎特法)[6 ] 、Occam法[7 ,8 ,9 ,10 ] 、正则化法[11 ,12 ,13 ] 等,都需要计算正演响应对模型的导数,即雅克比矩阵,因此雅克比矩阵的计算效率和效果就成为反演计算的关键因素之一.计算雅克比矩阵最简便的方法是差分法,该算法程序实现较为简单,缺点是不仅调用正演次数多,计算效率低,而且有可能增加累积误差.因此,本文详细推导了中心回线装置瞬变电磁一维正演全区视电阻率的雅克比矩阵的解析计算公式,并通过Occam反演,对比分析了理论模型和实测数据的反演效果,证明推导正确,值得采用. ...
回线源瞬变电磁法一维反演算法
1
2017
... 瞬变电磁法作为一种常见的电磁测深物探方法,广泛应用于水文勘查、工程勘查以及矿产勘查[1 ,2 ] .在目前的工程应用中,除“烟圈”反演快速成像法外,一维反演结果仍然是地质解释的主要依据.对于中心回线装置,可以解析推导出一维正演响应的计算公式[3 ] ,其数值计算的核心是,内层积分的汉克尔变换[4 ] 、外层积分的余弦变换[5 ] 均可采用数值滤波算法,随着滤波系数的不断推陈出新,计算量明显减少,大大提高了正演的计算效率和精度.在反演方面,目前常见的反演方法,如阻尼最小二乘法(马奎特法)[6 ] 、Occam法[7 ,8 ,9 ,10 ] 、正则化法[11 ,12 ,13 ] 等,都需要计算正演响应对模型的导数,即雅克比矩阵,因此雅克比矩阵的计算效率和效果就成为反演计算的关键因素之一.计算雅克比矩阵最简便的方法是差分法,该算法程序实现较为简单,缺点是不仅调用正演次数多,计算效率低,而且有可能增加累积误差.因此,本文详细推导了中心回线装置瞬变电磁一维正演全区视电阻率的雅克比矩阵的解析计算公式,并通过Occam反演,对比分析了理论模型和实测数据的反演效果,证明推导正确,值得采用. ...
时间域瞬变电磁法中心回线方式全程视电阻率的数值计算
1
2003
... 目前,瞬变电磁仪器大多采集的数据为感应电动势,对其归一化后可转换为晚期视电阻率,然而晚期视电阻率是假设在时间趋于无穷大时近似推导所得,因此在响应早期,视电阻率曲线变形严重,呈明显高阻假象,成为了地质解释重要的干扰因素之一,因此很多学者提出了全区视电阻率(也称全期视电阻率)的概念[14 ] .由于通过磁场强度计算的全区视电阻率更为接近真正的视电阻率定义[15 ] ,可通过感应电动势转换为磁场强度计算得出.因此,本文将该全区视电阻率作为Occam反演的拟合参数. ...
时间域瞬变电磁法中心回线方式全程视电阻率的数值计算
1
2003
... 目前,瞬变电磁仪器大多采集的数据为感应电动势,对其归一化后可转换为晚期视电阻率,然而晚期视电阻率是假设在时间趋于无穷大时近似推导所得,因此在响应早期,视电阻率曲线变形严重,呈明显高阻假象,成为了地质解释重要的干扰因素之一,因此很多学者提出了全区视电阻率(也称全期视电阻率)的概念[14 ] .由于通过磁场强度计算的全区视电阻率更为接近真正的视电阻率定义[15 ] ,可通过感应电动势转换为磁场强度计算得出.因此,本文将该全区视电阻率作为Occam反演的拟合参数. ...
瞬变电磁全区视电阻率数值计算方法研究
2
2010
... 目前,瞬变电磁仪器大多采集的数据为感应电动势,对其归一化后可转换为晚期视电阻率,然而晚期视电阻率是假设在时间趋于无穷大时近似推导所得,因此在响应早期,视电阻率曲线变形严重,呈明显高阻假象,成为了地质解释重要的干扰因素之一,因此很多学者提出了全区视电阻率(也称全期视电阻率)的概念[14 ] .由于通过磁场强度计算的全区视电阻率更为接近真正的视电阻率定义[15 ] ,可通过感应电动势转换为磁场强度计算得出.因此,本文将该全区视电阻率作为Occam反演的拟合参数. ...
... 由于函数Z (x )为隐函数,无法推导出其反函数,因此无法通过解析算法求解,但可利用函数在0~ 1之间单调递增的特点,通过二分查找的办法得到数值解.笔者曾采用该方法实现了全区视电阻率的计算以及对应的“烟圈”反演,取得了不错的效果[15 ] . ...
瞬变电磁全区视电阻率数值计算方法研究
2
2010
... 目前,瞬变电磁仪器大多采集的数据为感应电动势,对其归一化后可转换为晚期视电阻率,然而晚期视电阻率是假设在时间趋于无穷大时近似推导所得,因此在响应早期,视电阻率曲线变形严重,呈明显高阻假象,成为了地质解释重要的干扰因素之一,因此很多学者提出了全区视电阻率(也称全期视电阻率)的概念[14 ] .由于通过磁场强度计算的全区视电阻率更为接近真正的视电阻率定义[15 ] ,可通过感应电动势转换为磁场强度计算得出.因此,本文将该全区视电阻率作为Occam反演的拟合参数. ...
... 由于函数Z (x )为隐函数,无法推导出其反函数,因此无法通过解析算法求解,但可利用函数在0~ 1之间单调递增的特点,通过二分查找的办法得到数值解.笔者曾采用该方法实现了全区视电阻率的计算以及对应的“烟圈”反演,取得了不错的效果[15 ] . ...
1
1996
... 根据区域地质资料显示[16 ] ,工作区地层主要有太古宇界桑干群,中生界侏罗系及新生界古近系、新近系、第四系(图3 ).区域地质构造属阴山东西向构造带的南缘,主要有NE和NW两个方向,其中NE向构造属高角度压性断层,断层连续性强,是地下热水和温泉形成的主要通道.侵入岩主要分布在蛮汉山区及平原区深部,主要为太古宙代早期侵入岩,岩性以苏长岩和似斑状花岗岩为主. ...
1
1996
... 根据区域地质资料显示[16 ] ,工作区地层主要有太古宇界桑干群,中生界侏罗系及新生界古近系、新近系、第四系(图3 ).区域地质构造属阴山东西向构造带的南缘,主要有NE和NW两个方向,其中NE向构造属高角度压性断层,断层连续性强,是地下热水和温泉形成的主要通道.侵入岩主要分布在蛮汉山区及平原区深部,主要为太古宙代早期侵入岩,岩性以苏长岩和似斑状花岗岩为主. ...
凉城县岱海地区地热资源预可行性勘查项目取得突出成果
1
2014
... 上述地质推断,经钻孔DR2验证(图4 b),投影到剖面位置为820点,与实际情况基本吻合.该井定位于构造中,钻探深度 308.6 m,成井后自流,自流量为4 483.64 m3 /d,水温 36 ℃,矿化度 1.01 g/L,是内蒙古凉城县地热勘查中取得的突出成果[17 ] . ...
凉城县岱海地区地热资源预可行性勘查项目取得突出成果
1
2014
... 上述地质推断,经钻孔DR2验证(图4 b),投影到剖面位置为820点,与实际情况基本吻合.该井定位于构造中,钻探深度 308.6 m,成井后自流,自流量为4 483.64 m3 /d,水温 36 ℃,矿化度 1.01 g/L,是内蒙古凉城县地热勘查中取得的突出成果[17 ] . ...