首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
标量CSAMT只适合一维及测量方向与构造方向垂直的二维情况,对于复杂的三维地电结构,CSAMT需采用张量测量。本文试图采用矢量有限元法实现三维张量CSAMT的正演模拟。为了验证算法的正确性,本文在层状介质中计算了三维CSAMT远区的电场,磁场及阻抗张量,并且与层状介质中的理论解进行了比较,接着还模拟了均匀半空间中含有三维异常体的模型,并且分析了四个阻抗张量、视电阻率及阻抗相位的响应特征。得出如下结论:采用矢量有限元法来模拟三维张量CSAMT,其电磁场及阻抗张量的实虚部计算精度都比较高,并且该方法本身满足电场法向不连续,不用进行散度校正。  相似文献   

2.
三维CSAMT法非结构化网格有限元数值模拟   总被引:2,自引:1,他引:1       下载免费PDF全文
考虑到可控源音频大地电磁法(CSAMT)电偶极发射源与地下介质的三维结构特点,本文采用非结构化网格剖分技术,开展了三维CSAMT方法有限元数值模拟研究,将三维电磁场的背景场和异常场分别求解,避免了电偶极发射源的奇异性问题,并减小了计算区域.推导了三维异常电场遵循的有限元方程,加入散度条件进行约束以消除电场伪解;对非结构化网格单元采用高斯加权平均算法,得到了精度较高的异常磁场.针对层状介质模型,与积分方程法对比,验证了有限元算法的正确性;计算分析了典型三维地质模型的电磁响应,异常体反映明显.结果表明本文算法正确、可靠,适用于三维地质模型的CSAMT方法正反演研究.  相似文献   

3.
在地质情况复杂化以及地质勘探精细化的要求下,有必要开展张量CSAMT正反演研究.考虑到各向异性对正反演结果影响较大,本文开展基于主轴各向异性介质的张量CSAMT正演研究.根据含源的电场控制方程,引入主轴各向异性介质的张量电阻率,将电场总场分解为一次场和二次场,一次场由Key的算法得到,二次场利用有限差分法计算,即实现了主轴各向异性介质的张量CSAMT三维正演.通过与二维有限元算法结果比较来验证算法的正确性;通过设计几组主轴各向异性三维模型,分析卡尼亚电阻率和相位的响应特征.结果表明卡尼亚电阻率ρxy和ρyx分别对异常体X方向和Y方向的电阻率变化敏感,说明张量CSAMT方法对识别各向异性具有优势.本研究为开展任意各向异性介质的张量CSAMT正反演研究提供基础.  相似文献   

4.
局部畸变问题曾经困扰大地电磁资料反演解释几十年,大地电磁三维数值模拟技术的发展为剖析局部畸变特点和得到可靠的反演成像结果提供了技术基础。本文采用三维数值模拟成像方法对典型三维局部畸变模型进行模拟分析。三维数值模拟结果显示:电场分量垂直电性分界面的极化模式视电阻率曲线(对应二维情况下TM模式)在穿越低阻异常体界面时,曲线会先上移后下移,而在穿越高阻异常体界面时,曲线会先下移后上移,这与电性分界面处积累面电荷产生的二次电场有关。三维模型中XY模式、YX模式视电阻率和相位在三维异常体附近的水平变化是呈现近似垂向对称的,该现象与电场垂直跨越电性界面时视电阻率的变化规律是吻合的,当测线分别沿X方向和Y方向展布时,三维情况下的XY和YX模式分别对应二维情况下的TM模式。低阻小异常体对区域构造响应的畸变影响比高阻小异常体要严重。低阻小异常体对二维区域响应的两种极化模式视电阻率和相位都有非常明显的畸变影响,相比较而言对TE模式的畸变要大于TM模式,因此我们在做二维反演解释时,可优先考虑拟合TM模式数据。位于小异常体中心上方测点的三维畸变响应虽然与对应真实二维区域响应的差异比较大,但可以等效于某种二维模型响应,这种由局部畸变造成的假二维响应在实际野外数据的解释中是需要注意的。   相似文献   

5.
在实际复杂的三维地电条件下,倾子作为大地电磁法的重要解释参数之一,其正演模拟和响应分析对提高大地电磁法的探测精度具有重要作用.本文在简要阐述大地电磁三维正演基本理论的基础上,利用交错网格有限差分法开展了大地电磁三维倾子正演模拟研究,首先正演模拟单个低、高阻两种地电模型的倾子响应,结果表明倾子资料能够较好反映地下异常体的空间分布;在此基础上,设计了含有两个低(高)阻异常体和低阻垂直断层的复杂组合模型并计算了其三维倾子响应,结果表明,倾子具有对复杂的地电结构良好的空间分辨率,且倾子响应既保留了模型中每单个异常体的倾子响应形态,又表现出了异常体之间相互作用的整体响应情况,进一步加深了对倾子的响应特征和规律、以及倾子资料对异常体边界的识别能力的认识,为倾子资料从定性解释向定量解释发展提供参考和依据.  相似文献   

6.
可控源音频大地电磁法(CSAMT)目前实际应用多以单个源采集标量数据为主,这种工作方式限制了CSAMT探测复杂地质构造的能力.随着CSAMT张量观测系统的推出和精细电法勘探的需求增多,CSAMT张量采集逐步受到重视.为了深化对CSAMT三维张量阻抗响应特征的认识,本文基于电场含源控制方程,采用交错网格有限差分法实现了CSAMT的三维正演算法;与二维CSAMT正演结果对比,验证了三维正演的正确性;设计了代表性三维地电模型,计算了标量和张量阻抗,对比和分析了标量和张量阻抗的响应特征.  相似文献   

7.
应用CSAMT方法进行标量测量时,一般要求测线沿垂直构造走向的方向布测.在实际测量的时候由于地表条件复杂经常出现无法判断地下构造走向的情况,例如,在长白山地区进行CSAMT测量,由于表面覆盖有玄武岩无法判断地下构造走向,对得到的标量数据进行一维、二维反演处理时无法判断哪种反演方式的结果更准确.鉴于此,建立三维地质模型,应用CSAMT三维正演程序进行数值模拟.研究测线与构造走向平行,呈45°及垂直三种情况下的一维、二维反演,对结果做比较分析.实验结果表明,测线平行构造走向时,一维、二维反演结果近似;测线与构造走向呈45°或垂直时,二维反演结果比一维反演结果对异常体的形状,大小,埋深及电性参数反演的更准确.  相似文献   

8.
电磁场数值模拟的背景场/异常场算法是三维正演的有效策略之一,优点为采用解析法计算电磁场背景场代替场源项、克服了场源奇异性,缺点为不适用于发射源布置于起伏地表或背景模型复杂的情形.总场算法是直接对电磁场总场开展数值模拟,其难点是有效加载场源、保证近区与过渡区数值解精度.本文以水平电偶源形式分段加载接地长导线源,并以电场总场Helmholtz方程为矢量有限元法控制方程,实现了基于非结构化四面体网格剖分的接地长导线源频率域电磁法三维正演.通过与均匀全空间中水平电偶源产生的电场解析解对比,验证了本文算法的正确性,并分析了四面体外接圆半径与其最短棱边的最大比值和四面体二面角最小值对数值解精度的影响规律.通过与块状高导体地电模型的积分方程法、有限体积法和基于磁矢量势Helmholtz方程的有限元法数值解对比,进一步验证了本文算法正确性,同时说明了非结构化四面体网格能够更加精细地剖分电性异常体,利于获得精确数值解.  相似文献   

9.
起伏地形下CSAMT二维正反演研究与应用   总被引:13,自引:5,他引:8       下载免费PDF全文
雷达 《地球物理学报》2010,53(4):982-993
CSAMT在山区金属矿勘查中,采用各种滤波和相位积分之类的处理方法,校正因地形起伏和局部电性不均匀引起的静态效应,往往难保奏效,开发消除静态效应的新方法是提高CSAMT资料处理与解释水平和方法应用效果的重要研究课题. 本文以如何消除地形影响为重点,对起伏地形下CSAMT二维大地三维源地电模型,采用加权余弦数值积分法,进行波数域电磁场二维有限单元法正演. 为模拟复杂地形地电模型,选取交叉对称网格三角形剖分法,实现了在国内常用赤道电偶极装置的CSAMT二维正演计算;在二维正演的基础上,开发了基于奥克姆反演法的CSAMT二维反演技术,研制出一套起伏地形下CSAMT二维正反演处理与解释方法技术系统. 通过理论模型试算和实测数据处理证实,本系统能有效地削减起伏地形影响. 在找矿应用中,该系统反演的电阻率断面,极大地消除了起伏地形影响和静态效应,突显出清晰的控矿构造和矿体的异常,取得了重要成效.  相似文献   

10.
从电偶源三维地电断面可控源电磁法的二次电场边值问题及其变分问题出发,采用任意六面体单元对研究区域进行剖分,并且在单元分析中同时对电导率及二次电场进行三线性插值,实现电导率分块连续变化情况下,基于二次场的可控源电磁三维有限元数值模拟.这个新的可控源电磁三维正演方法可以模拟实际勘探中地下任意形状及电性参数连续变化的复杂模型.理论模型的计算结果表明,均匀大地计算的视电阻率误差和相位误差分别为0.002%和0.0005°.分层连续变化模型的有限元计算结果表明,其与对应的分层均匀模型解析结果有明显差异.三维异常体组合模型以及倾斜异常体等复杂模型的有限元计算结果也有效地反映了异常形态.  相似文献   

11.
复杂介质有限元法2.5维可控源音频大地电磁法数值模拟   总被引:38,自引:13,他引:25       下载免费PDF全文
利用有限元2.5维可控源音频大地电磁法(简称CSAMT)数值模拟方法,对100Ωm均匀半空间介质中有限长度的电偶极源产生的电场、磁场及视电阻率、相位特征进行了数值模拟,研究了场的空间变化规律. 在一个象限中,场的特征存在双叶现象,当收发距大于4个趋肤深度时,电阻率较接近介质真实的电阻率,这些结果为观测系统和收发距的选择提供了依据. 波数域场的特征表明,低波数对源的贡献占较大的比例. 有限元法2.5维CSAMT数值模拟的优势在于能较准确地获得复杂介质结构的波场特征. 本文结合直立异常体、倾斜异常体及断陷模型对CSAMT电阻率、相位剖面特征及频率曲线特征的可靠性进行了研究. 数值模拟结果直观地给出了异常体的剖面异常形态. 通过对比研究异常体的剖面异常形态和半空间场的特征进一步说明本文方法和软件在模拟复杂介质结构场特征时是可靠的. 这为认识观测数据,指导反演解释提供了较好的依据.  相似文献   

12.
三维电磁模拟技术及其在频率测深法中应用   总被引:3,自引:0,他引:3       下载免费PDF全文
利用现有的并矢格林函数和积分方程理论对轴向频率测深三维问题进行正演模拟,并讨论其精度.三维良导体或高阻体的存在对频率测深电场分量E定义的视电阻率响应(振幅、相位)分布特征的影响各有其特殊的规律性.选择测量电极中心作为频率测深三维问题正演计算结果的记录点,既使视电阻率响应在异常体正上方获得最大异常,同时又使正演计算的拟断面图的异常分布与异常体之间有很好的对应关系.  相似文献   

13.
三维电磁模拟技术及其在频率测深法中应用   总被引:4,自引:0,他引:4       下载免费PDF全文
利用现有的并矢格林函数和积分方程理论对轴向频率测深三维问题进行正演模拟,并讨论其精度.三维良导体或高阻体的存在对频率测深电场分量Ex定义的视电阻率响应(振幅、相位)分布特征的影响各有其特殊的规律性.选择测量电极中心作为频率测深三维问题正演计算结果的记录点,既使视电阻率响应在异常体正上方获得最大异常,同时又使正演计算的拟断面图的异常分布与异常体之间有很好的对应关系.  相似文献   

14.
海洋可控源电磁三维非结构矢量有限元数值模拟   总被引:14,自引:12,他引:2       下载免费PDF全文
杨军  刘颖  吴小平 《地球物理学报》2015,58(8):2827-2838
本文实现了海洋可控源电磁三维矢量有限元数值模拟.由于采用非结构四面体单元进行三维网格剖分,该方法可以模拟复杂电性异常体和海底地形.一维模型的数值模拟结果表明,电场实、虚部均与解析解吻合得相当好,计算误差基本小于1%.二维模型的计算结果与已有的二维自适应非结构有限元模拟结果吻合很好.带地形模型的数值模拟结果显示,海底地形对电场影响很大,有可能掩盖海底油气藏产生的异常.  相似文献   

15.
CSAMT三维交错采样有限差分数值模拟   总被引:5,自引:1,他引:4  
系统分析大地电磁三维交错采样有限差分算法的基础上,根据可控源音频大地电磁场特征,采用将总场分解为一次场和二次场计算,一次场利用快速汉克尔变换,二次场利用数值模拟的思路.从CSAMT满足的麦克斯韦方程组积分形式出发,利用交错采样有限差分算法推导了电场和磁场的离散关系式,提出了简洁的边界条件和合理的剖分方案,所实现的CSA...  相似文献   

16.
地空电磁法已经成为深部资源勘探的重要地球物理方法,但对频率域地空系统的三维多源电磁响应特征研究较少.本文设计了多种激励源组合方式,采用非结构化有限元数值模拟方法,对三维地电模型的空中垂直磁场的响应特征进行了研究.首先推导了基于电场的双旋度公式及其变分形式,加入罚项以减少伪解的影响.接着把有限元稀疏矩阵方程转换为频率的函数,采用Krylov子空间投影方法,通过模型降阶算法降低稀疏矩阵的阶数,实现多频点的快速计算.建立了三维低阻体模型、高阻体模型以及两个相邻低阻体模型,分别采用单源、双源、三源和四源激励模式,从垂直磁场的总场、二次场响应和全域视电阻率等方面进行分析比较.结果表明:多源地空电磁法不仅可以增加总场的强度,而且可以改变异常体的二次电磁响应分布规律.各电偶源延长线呈正三角形分布的三源和矩形分布的四源激励模式在增强信号强度以及削弱异常体的边界效应方面具有一定的优势,是一种优化的多源激励方式.  相似文献   

17.
本文研究并建立了一种模拟各向异性地层中三维电性异常体电磁响应的积分方程算法.首先讨论了并矢Green函数及其相关积分的计算,将水平层状各向异性地层中的电场并矢Green函数分解成含有奇异项的直达波与非奇异的来自各个层界面的反射和透射波两个部分,再应用等效体积单元和表面积分技术对积分方程的奇异核进行离散化处理以便提高离散方程的精度.然后为了节省计算机内存以及计算时间,引入基于Krylov子空间的迭代算法求解积分方程的离散化矩阵方程.最后通过与现有文献中的结果作对比从而检验了所述算法的有效性,并结合具体算例考察分析了地层的各向异性对三维电性异常体电磁响应的影响特征和规律.  相似文献   

18.
可控源音频大地电磁测深法(CSAMT)基本上是基于电性各向同性理论的假设,而实际地层普遍表现为电性各向异性现象,因此研究各向异性地层条件下CSAMT方法的电磁响应具有重要的理论意义和应用价值.本文从有源Maxwell方程出发,推导出了二维各向异性地层条件下无限长线源频率域的电磁场响应,得到了一组平行于线源方向的电场分量E_x和磁场分量H_x的偏微分方程,并采用有限元法进行求解.模型的计算结果表明:不论是对均匀半空间还是一维层状介质中的各向异性二维体模型的计算,均能从结果中明显地看出各向同性背景场和各向异性异常体的存在;在计算中通过改变各向异性系数或各向异性旋转角,其结果在视电阻率和相位曲线上均存在显著区别.计算结果表明了在各向异性地层中,仍然采用各向同性的假设进行CSAMT资料的处理与解释,可能会带来比较严重的误差,这种情况下必须采用基于各向异性的地质模型对CSAMT资料进行处理和解释.研究结果对于提高CSAMT方法的勘探效果和应用水平具有重要的指导意义.  相似文献   

19.
利用三维(3D)全空间积分方程准线性解作为数值模拟手段,对包含电离层、大气层和地球介质层的“地-电离层”典型异常目标体多层介质模型进行了数值模拟,得到了偶极源长度100 km、电流200 A、收发距离远达1600 km的合理的异常电性目标体的电阻率-频率响应结果.通过对不同埋深异常目标体的电阻率-频率特征的讨论,认为当考虑电离层和大气层的“地-电离层”大尺度深层横向不均匀复杂介质模拟时,电磁场对深部目标体仍有很好的异常响应,但当异常体电阻率及同一电阻率的岩石埋深不同时,其电阻率-频率的响应特征有很大不同,这正反映了不同岩石物性参数的电阻率-频率响应特征.文中结果表明从长偶极大功率源(WEM)激励的电磁场观测资料区分岩石的电阻率参数是可能的,这为利用WEM观测资料建立电性参数和岩石/地层参数经验关系奠定了基础.  相似文献   

20.
地-井瞬变电磁法是在地面发射,井(钻孔)中接收的装置形式,能够利用已有钻孔使接收探头深入地下更加接近矿体,获得更加可靠的目标体信息.文中针对复杂地形地区矿区深部找矿面临的实际问题,提出了电性源地-井瞬变电磁法的思路,并以FDTD正演模拟来研究均匀半空间瞬变场的空间分布以及一维、三维模型钻孔中的三分量瞬变响应特征,以期获得理论模型三分量响应特征规律,为实际勘查提供参考.理论模型计算结果表明,电性源地-井TEM方法三分量瞬变响应曲线对于电性界面、异常体具有良好的反应.  相似文献   

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

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