共查询到18条相似文献，搜索用时 361 毫秒
1.
《物探与化探》2015,(5)
采用三角网格剖分的有限元法,研究了2.5维航空瞬变电磁法正演模拟问题。利用时频变换数值方法将时间域电磁场转换到拉氏域,再利用傅里叶变换将三维问题降维变为2.5维问题,然后由有限元法求解得到拉氏域二维电磁场,逆拉氏变换后得到时间域航空瞬变响应。为了回避正演模拟中总感应磁场在场源处的奇异性问题,采用异常场算法,场源响应通过在微分方程中施加背景电磁场实现。由于瞬变电磁信号具有较大的动态范围,而且需要经过两次正、逆拉氏变换和傅里叶变换,每个环节的计算精度和速度要严格控制在较高的水平上,否则积累误差会非常大。模型计算表明均匀大地和层状大地模型解析解与数值解吻合很好。这证明该算法是正确可行的,可作为研究二维复杂地质体的方法手段。 相似文献
2.
利用G-S逆拉氏变换法计算瞬变测深正演问题 总被引：3，自引：0，他引：3
求解时间域电磁法问题,通常有两种途径。第一,直接求解时间域问题;第二,先解相应的频率域问题,然后变换为时间域问题。本文采用第二种方法,利用Gaver-stohfest逆拉氏变换法将频率测深正演问题转化为对应的瞬变测深正演问题。实际计算表明,该方法是将频率域问题转化为时间域问题的有效方法。 相似文献
3.
TEM正演响应计算的几种频-时域转换方法对比 总被引：1，自引：0，他引：1
在瞬变电磁的一维正演响应模拟中,常用的方法是先在频率域中求解,之后将结果转换到时间域。但该方法在晚期的计算精度通常不高,因此,使用5种频-时域转换方法(正、余弦变换的数值滤波算法,G-S逆拉普拉斯算法、正、余弦变换的折线逼近法)进行了计算,与解析解对比,得出余弦变换的数值滤波算法在晚期计算中精度最高的结论,并对这五种转换方法产生误差问题的原因进行了讨论和分析。本研究有利于瞬变电磁一维正演响应的高精度计算,使其在多维计算中得到更好的应用。 相似文献
4.
国外瞬变电磁法中正反演问题研究现状 总被引：2，自引：0，他引：2
国外在瞬变电磁法正反演问题研究中,运用有限差分法、有限单元法及积分方程法直接在时间域中求解以及运用富氏变换法、衰减谱法和数字滤波法从频率域变换等数值计算方法,并对其在瞬变电磁资料解释中的效果作了评述. 相似文献
5.
《物探化探计算技术》2021,(3)
地空瞬变电磁法(Ground-airborne transient electromagnetic method)是一种新兴的物探方法,其探测范围大,工作效率高,在地形复杂的区域有着传统物探方法无法比拟的优势。这里通过对地空瞬变电磁的一维正演问题的研究,采用Hankel变换、Gauss积分和频-时转换实现了地空瞬变电磁时间域响应正演计算,通过解析解和数值解的对比验证了正演方法的准确性和稳定性,并研究了其响应特征。利用最小构造模型约束下的自适应正则化反演方法,进行了地空瞬变电磁法反演模拟研究,证实了该方法对低阻层的反映比较真实,对高阻层的反映略有误差的特点。最后用该方法对某滑坡实测数据进行了反演解释,反演结果与工区地质情况符合较好。 相似文献
6.
7.
张秋瑞 《华北地质矿产杂志》1998,13(3):276-284
本文以垂直磁偶极子发射,接地导线测量这样一种发一收装置为例,讨论了求解时间减电磁场问题的计算方法,具体思路是：用线性滤波法,先求解频率域问题,然后且Ｇａｖｅｒ－Ｓｔｅｈｆｅｓｔ逆拉氏变换法将频率测深正演问题转化为相对应瞬变测深（时间域）正演问题,得到了垂直磁极子场中,水平层状介质表面瞬变电磁场分量的近似表达式,由此计算了视电阻率曲线,并从理论上探讨了瞬变测深视电阻率曲线的垂向分辨能力,本文试图对中 相似文献
8.
电偶源瞬变电磁测深研究（二）——瞬变电磁场的求解方法 总被引：4，自引：1，他引：3
主要指出了瞬变电磁场求解的几种方法,重点阐述频率域的方法,也就是利用已有的频率域的解,通过傅里叶变换或拉普拉斯变换的逆变换求得时间域的瞬态解。 相似文献
9.
根据不同时域电磁场的频率特性,从频率域到时间域的转换一般有二种方法,即线性滤波算法(余弦变换)和快速傅氏变换(FFT)。利用高密度采样的线性滤波算法计算余弦变换,对层状介质的甚早期瞬变电磁进行模拟,并与均匀半空间的闭合解析式结果进行对比,证明了该方法的有效性。另外,从层状介质的高频电磁场响应出发,利用FFT算法,有效地模拟了不同频率和地电参数的GPR射线,以及层状介质的剖面图。甚早期瞬变电磁和低频GPR的频率,都跨越了传统电磁法中所忽略的中间频率。这里的主要目的就是通过对比选择有效的模拟方法,利用模拟结果来分析该频段时域电磁场的特征。 相似文献
10.
11.
为提高电磁波场数值计算精度,对时间域有限差分法与时间域伪谱法进行了对比研究。时间域有限差分法是一种目前流行的电磁场时域数值计算方法,已被广泛应用于求解与时间有关的偏导数方程。对于大规模数值计算,时间域有限差分法需要较多的内存空间。时间域伪谱法基于时间域有限差分法,该方法使用快速傅里叶度换来计算麦克斯韦方程中的空间导数。由于傅里叶变换的准确性,时间域伪谱法使计算精度提高,数值计算时所需的格子数大大减少,这极大地节省了计算机内存空间,适合于大规模正反演问题的数值计算。 相似文献
12.
在考虑横观各向同性含液饱和多孔介质固体骨架和流体可压缩性以及固体骨架的黏弹性特征下,基于横观各向同性含液饱和多孔介质u-w形式的三维动力控制方程,以固相位移u、液相相对位移w为基本未知量,综合运用Laplace变换、双重Fourier变换等方法,在直角坐标系下通过引入中间变量,将六元2阶动力控制方程组化为两组各含4个未知变量的常微分方程组,给出了直角坐标系下横观各向同性含液饱和多孔介质三维黏弹性动力反应的积分形式一般解;作为理论推导的验证,通过引入初始条件和边界条件,对横观各向同性含液饱和多孔介质半空间黏弹性瞬态反应问题进行了求解。解答的退化验证表明,所推导的理论解是正确的。 相似文献
13.
斜阶跃波层状介质航空瞬变电磁响应数值计算 总被引：1，自引：0，他引：1
目前,航空瞬变电磁方法数据解释主要采用阶跃波均匀半空间模型计算视电导率值,而实际航空电磁系统发射波形的下降沿多为斜阶跃,导致解释时计算的视电导率值存在较大偏差。为此,笔者研究了航空瞬变电磁系统发射电流为斜阶跃时的电磁响应数值计算,将发射波形进行拉氏变换,利用G-S逆拉氏变换与241点汉克尔变换相结合的方法,实现斜阶跃波关断后的层状大地模型航空瞬变电磁响应计算;并对均匀半空间和层状大地模型下,不同关断时间和不同飞行高度对电磁响应的影响进行分析。得出结论:不同关断时间,关断后取样延时2 ms时,均匀半空间电磁响应的平均偏差为27.78%,三层模型的平均偏差为32.16%;当飞行高度从20 m增加到60 m时,均匀半空间和三层模型的感应电动势分别减小了43.6%和83.2%。 相似文献
14.
Transient dynamic response of multilayered saturated media subjected to impulsive loadings 
下载免费PDF全文
下载免费PDF全文 This paper presents a stable and efficient method for calculating the transient solution of layered saturated media subjected to impulsive loadings by means of the analytical layer element method. Starting with the field equations based on Biot's linear theory for porous, fluid‐saturated media, and the seepage continuity equation, an analytical layer element for a single layer is established by applying Laplace‐Hankel integral transform. The global stiffness matrix in the transform domain for a layered saturated half‐space subjected to a transient circular patch loading is obtained by assembling the layer elements of each layer. The displacements in the time domain are derived by Laplace‐Hankel inverse transform of the global stiffness matrix. Numerical examples are conducted to verify the accuracy of the method and to demonstrate the influences of type of transient loading, buried depth of loading, permeability, and stratification of materials on the transient response of the multilayered saturated poroelastic media. 相似文献
15.
Unsaturated soils are considered as porous continua, composed of porous skeleton with its pores filled by water and air. The governing partial differential equations (PDE) are derived based on the mechanics for isothermal and infinitesimal evolution of unsaturated porous media in terms of skeleton displacement vector, liquid, and gas scalar pressures. Meanwhile, isotropic linear elastic behavior and liquid retention curve are presented in terms of net stress and capillary pressure as constitutive relations. Later, an explicit 3D Laplace transform domain fundamental solution is obtained for governing PDE and then closed‐form analytical transient 3D fundamental solution is presented by means of analytical inverse Laplace transform technique. Finally, a numerical example is presented to validate the assumptions used to derive the analytical solution by comparing them with the numerically inverted ones. The transient fundamental solutions represent important features of the elastic wave propagation theory in the unsaturated soils. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
16.
准饱和黄土中暗穴水动力扩展响应的研究 总被引：1，自引：0，他引：1
为分析强降水过程中暗穴扩展机理,在野外调查资料的基础上,建立了准饱和黄土暗穴的扩展计算模型,概化为弹性准饱和土层中无限长圆柱形孔洞表面受水压力的动力响应问题。通过引入势函数,得到了Lap lace变换域中的应力、位移及孔隙水压力的解析表达式,并利用数值逆变换方法求得时域解,分析了降雨历时及饱和度等因子对暗穴动力响应的影响。结果表明,水动力是引起黄土暗穴扩展的主要因素,水压力作用历时对位移及应力有较大的影响,同时饱和度的细微变化对径向位移有显著影响。 相似文献
17.
针对探地雷达子波提取问题,提出了一种基于分数阶傅里叶变换的子波提取算法。假设探地雷达地下介质模型为线性系统并是时变的,并且子波在传播过程中不变,把子波和地层信息描述为一个时间的积分形式。分数阶傅里叶变换是酉变换,具有保范性,把子波的提取问题转化为在分数阶傅里叶变换域上的最优滤波问题。对分数阶傅里叶变换的阶数p在(－2,2]范围以1/300为步长进行了搜索,在分数阶傅里叶变换域上找到最优解,最后经过分数阶傅里叶变换反变换得到子波估计的时域信号。仿真实验首先验证了分数阶傅里叶变换最优滤波算法的正确性,然后利用探地雷达实验数据提取了子波,验证了该算法在提取子波中的有效性。 相似文献