首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到17条相似文献,搜索用时 171 毫秒
1.
数值模拟是研究探地雷达电磁波在多层介质中传播规律的有效手段,传输线算法是一种基于Huygens原理的时域分析电磁场方法。这里阐述了TLM方法的基本原理及基本关系式,用TLM方法对多层介质进行了正演模拟,将模拟结果与FDTD方法合成结果进行了对比,结果表明了TLM用于GPR正演计算的正确性及有效性。应用该算法对已有模型进行了定性分析,验证了随着介电常数和电导率增大,雷达电磁波衰减加快的结论。  相似文献   

2.
大量观测资料表明,地下介质的物性普遍存在非各向同性现象,只有建立非各向同性模型才能真实地反应地下介质情况。在分析前人研究成果的基础上,从实际地质情况出发,以层状介质的电学性质为基础引入非各向同性系数到对称非各向同性体的研究当中,推导出对称非各向同性体二维大地电磁场的变分方程,实现了模拟对称非各向同性体的有限单元法。将计算结果与Pek的有限差分法的计算结果进行对比验证了程序的正确性。通过固定非各向同性系数在不同倾角下的模型,说明倾角对非各向同性效应的影响;通过固定倾角在不同非各向同性系数下的模型,说明非各向同性系数对非各向同性效应的影响;通过不同倾角和不同非各向同性系数模型的正演拟断面图,更加形象地展示两者与非各向同性效应之间的关系。  相似文献   

3.
空气介质和电离层是电导率不同的分层媒质,当研究由空气层中电流或磁矩激发的电磁场时,地面接收到的电磁波不仅受到空气介质、地下介质的影响,而且还受到电离层的影响。首先,推导出了电磁场一般边值问题,然后,将电离层、空气及地壳视为均匀、线性、各向同性的三层媒质,考虑地壳、空气、电离层的相互耦合,利用分离变量法推导求得了位于地壳表面水平谐变电偶极子的电磁响应表达式,通过这些表达式能够方便地求出电离层影响下均匀半空间空气层中的电磁场场值,同时,为后续的数值计算奠定了理论基础。  相似文献   

4.
在地震传播理论中,地震波场的正演数值模拟一直是研究的热点。而在正演数值模拟方法的研究中,计算精度和计算效率是评价此方法的有效性及优越性的二个关键问题。伪谱法在计算精度与计算效率方面优越性十分明显,常用的有Fourier变换法和Hartley变换法;虽然Hartley变换法在求取导数时较Fourier变换法复杂,但由于Fourier变换法的计算同时涉及复数的实部与虚部,在计算速度和占用内存方面不如Hartley变换法。这里详细地阐述了利用Hartley变换求解正交各向异性介质波动方程的数值方法,并模拟了多种正交各向异性介质模型,对比分析了地震波在各向同性介质与正交各向异性介质中的传播差异。结论表明,该方法能正确、高效、直观地反应地震波在正交各向异性介质中的传播规律。  相似文献   

5.
由于弹性波逆时偏移更符合实际情况,而且转换波的成像结果有更高的分辨率,因此弹性波逆时偏移越来越受到重视。然而,弹性波逆时偏移需要多波多分量数据,为了减少成像结果中的串扰假象,在逆时偏移过程中进行P波和S波分解就变得非常必要。结合基于向量的激发振幅成像条件,我们把基于解耦传播的波场分解方法应用到弹性波逆时偏移中,并对比了其在各向同性介质和横向各向同性(VTI)介质中的应用效果。结果说明,该方法可以在各向同性介质中完全分解P波和S波,并保留波场的向量信息。虽然在VTI介质中有较小的分解残余,但是该分解残余不会在逆时偏移结果中产生明显的串扰;因此,这种波场分解方法可以应用于各向同性介质和VTI介质弹性波逆时偏移。该方法是在时间空间域实现的,可以在波场传播过程中直接对P波和S波进行分离,应用方便,计算效率高。与利用Helmholtz分解的弹性波逆时偏移相比,该方法避免了在PS波成像结果中的极性反转问题。在复杂Hess VTI模型的逆时偏移结果中,高速岩体和断层的成像清晰,甚至是两个低速薄夹层也能较好成像;这说明该方法对复杂介质具有较好的适应性。PS波成像结果中的各向异性体成像清晰,说明各向异性介质弹性波逆时偏移可以对传统逆时偏移不能很好成像的构造进行成像。  相似文献   

6.
基于磁导率的频率域航空电磁法双频反演方法   总被引:1,自引:0,他引:1  
以往的频率域航空电磁法数据反演方法中,通常在模型中以真空磁导率代替介质磁导率,忽略了磁性层影响。但在高磁性区域,这样做将会导致计算电阻率值偏高。针对这一问题,首先根据层状介质模型研究了考虑磁导率和电阻率因素条件下的电磁场正演计算方法,并且分析了不同磁性介质对频率域电磁响应的影响规律。在此基础上,提出了磁导率和电阻率的双频反演方法,电阻率是通过2个频率的虚分量之比计算得到的,计算结果更加可靠,提高了磁性区域视电阻率填图的精度。通过对模型计算和实测数据进行反演计算,表明这种方法具有较高的反演精度,为强磁性区域航空电磁法视电阻率填图提了供一种行之有效的反演方法。  相似文献   

7.
双相介质中纵波方程的高阶有限差分解法   总被引:6,自引:0,他引:6  
张会星  何兵寿  宁书年 《物探与化探》2004,28(4):307-309,313
从双相介质中的纵波方程出发,导出了求解双相各向同性介质中纵波方程的高阶差分格式,给出了吸收边界条件和稳定性条件,在此基础上实现了双相各向同性介质中纵波方程的高阶有限差分法正演模拟,数值模拟结果表明,这种算法能在少量增加计算量的前提下大大提高精度,算法可同时应用于叠前和叠后的数值模拟。  相似文献   

8.
在大地电磁资料处理和解释中,大地电磁的各向异性正演一直是国内、外研究的前沿课题。这里首先从麦克斯韦方程组出发,得出各向异性介质二维大地电磁场的边值问题以及等价的变分问题,进而对计算区域进行完全非结构化三角形剖分,在单元内采用线性插值,将变分方程转化成线性方程组,求解出有限单元法数值解,获得各向异性介质下的电磁场值。讨论了电性主轴与笛卡尔坐标轴之间的夹角以及地形的起伏变化对视电阻率和相位值的影响,结果表明:起伏地形条件下,TE和TM两种极化模式的视电阻率和相位值都会受到影响,且TM模式受地形影响较大;由于介质的各向异性,视电阻率和相位断面图出现明显不同于各向同性的形态,导致横向分界面模糊。  相似文献   

9.
电提取法基础理论研究   总被引:1,自引:1,他引:1  
李金铭  卢军 《物探与化探》1992,16(4):275-285
本文从电解质溶液的基本电化学性质出发,根据传导类电法的电场分布理论,以均匀各向同性介质为例,对点源电流场中离子迁移一定距离所需时间和电极(液态、固态元素接收器)提取的物质质量与哪些因素有关的问题,作了较详细的讨论。并在此基础上,结合时量曲线的模拟实验结果给出了离子淌度、迁移距离和迁移浓度的计算方法。  相似文献   

10.
介绍了联合运用积分方程法与有限元法(简称混合法)来计算三维非均匀介质中电磁场分布的理论方法,并进行了数值模拟实验。混合法的原理是引入一个包围非均匀目标体的虚构边界,在边界内部的场用有限元法模拟,在边界外部(包括边界)的场用积分方程表达,二者在边界上通过场的连续性耦合起来。数值实验结果表明,混合法既能显著地减小网格规模,又能灵活地模拟复杂的介质情况,且计算精度较高。  相似文献   

11.
回线源瞬变电磁测深快速反演方法   总被引:4,自引:3,他引:4  
回线源瞬变电磁测深法(简称瞬变电磁法或TEM法)是以不接地回线为场源,在回线内、外测量电磁场的一种电磁测深方法,已经广泛用于固体矿产、工程地质和水文地质调查工作之中[1,2]。目前,回线源瞬变电磁法主要采用中心回线和重叠回线两种测量装置,测量回线内垂直磁场产生的感应电动势。“烟圈”反演方法是瞬变电磁反演的主要计算方法,“烟圈”反演利用地表瞬变电磁响应与某时刻电流源镜像等效原理,计算勘探深度和对应电阻率。这里给出了回线源瞬变电磁测深一维“烟圈”反演解释方法的基本原理、计算方法和应用实例。结果表明,一维“烟圈”反演解释结果纵向分辨率较高,能反映出地下电性变化情况,可以用于瞬变电磁测深资料解释工作。  相似文献   

12.
By using the numerical method to model the ter-rain effect on the magnetotelluric field,few resultshave been obtained. The finite element method(FEM) was used by Chouteau and Bouchard (1988)and Wannamaker et al .(1986) ,andthe boundary el-ement method (BEM) was used by Xu and Zhou(1997) and Xu (1995) to model 2Dtopographyinflu-ences on magnetotelluric surveys . The BEM methodwas also used to model the 3Dtopographic effect onmagnetotelluric deep sounding (Xu et al .,1997 ;Xu,1995) .In t…  相似文献   

13.
三维地形大地电磁场的边界元模拟方法   总被引:1,自引:0,他引:1  
提出了一种用边界元法计算大地电磁场三维地形影响的数值模拟方法.首先用矢量积分理论和电磁场边界条件, 将上半空间(空气)和下半空间(地下介质)两个区域电磁场边值问题变为仅对地形界面的两个矢量面积分方程, 其中一个计算磁场, 称磁场方程; 另一个计算电场, 称电场方程.然后将对地形界面的积分剖分为一系列的三角单元积分.在三角单元积分中, 假设单元中电磁场为水平均匀大地空间电磁场与地形影响的迭加, 并假设地形影响为常项, 这样既保证了计算精度又使得计算方法简便.通过分解和计算, 每一个矢量面积分方程分解为对应3个坐标方向的3个常量线性方程, 这些线性方程组成了对角占优的线性方程组, 可用SSOR方法求解.文中给出了2个三维地形上大地电磁视电阻率曲线的计算结果.   相似文献   

14.
利用时频分析研究瞬变测深中的瞬变信号,已有多种方法给出了定性的结果。但更深入的获取瞬变测深中的地下信息,还需给出定量描述。本文利用S变换和广义S变换研究了瞬变测深信号,并与其他结果做了比较,给出了S变换下瞬时频率(IF)估计的方法。利用这种方法研究了瞬变测深,同时讨论了瞬变测深信号在S域的去噪问题。理论模型与实例均给出了理想的结果。   相似文献   

15.
多层介质点电流源磁场响应计算   总被引:1,自引:0,他引:1  
介绍了点电流源层状介质模型磁电阻率法响应的模拟方法。由电流磁场的旋度公式出发,推导了多层介质中点电流源磁场的表达式,其为一个关于空间径向坐标和贝塞尔函数空间波数的二重积分。利用贝塞尔函数的性质可将该二重积分化为标准的贝塞尔积分,然后采用直接数值积分的方法计算该积分。在此过程中通过连分式展开加快积分求和的收敛速度。利用上述方法对海底天然气水合物模型点电流源产生的磁场进行了计算。模型计算结果表明,海底天然气水合物的存在能引起磁场较为明显的相对变化。  相似文献   

16.
时间谱电阻率法的二维正演算法   总被引:5,自引:0,他引:5       下载免费PDF全文
孟永良  罗延钟 《地球科学》2000,25(6):656-662
时间谱电阻率(TSR)法的二维正演, 是计算可极化二维地电构造上三维电流偶极源的电场瞬变响应, 属所谓2.5维时间域电磁场数值模拟问题, 是目前国际上未妥善解决的计算地球物理疑难问题.针对现有算法的局限性, 建立了新的算法, 其特点是: (1)采用二维有限单元算法, 在矩形网格中增加两对角线形成三角网格剖分, 同时用高斯消元法消除矩形网格中心结点的待求未知量.这样, 既可较准确地模拟任意二维复杂地电断面, 又可节省计算量.(2)采用直接计算二次场的新算法, 只需计算电场和磁场沿地电构造走向两个分量的一次场, 因而, 不但计算精度较高, 而且不显著增加计算量.(3)采用G -S变换法作逆拉氏变换, 并利用拉氏变换延迟定理在倍增的时间间隔中插值, 从而实现对密集采样时间瞬变过程的快速计算.(4)能对可极化和导电大地(即同时包括IP和EM效应).   相似文献   

17.
时域瞬变场电磁场有限差分法   总被引:1,自引:2,他引:1  
徐凯军  李桐林 《世界地质》2004,23(3):301-305
模拟电磁场在地下介质中的传播规律是理解地下介质响应的重要手段,直接应用有限差分法在时间域对二维、三维瞬变场进行数值分析是一种有效方法。应用反映电磁场基本规律的麦克斯韦方程组,导出时域场量的齐次扩散方程,对所研究的空间区域作差分离散,源作为初始条件加入,采用合适的边界条件,利用差分方程进行计算,从而得到场量的数值结果,展现瞬变场在地下随时间扩散的全过程,对电磁场的瞬态响应有更加直观的了解。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号