共查询到20条相似文献，搜索用时 218 毫秒
1.
由于激发极化效应的影响,时间域航空电磁晚期道信号经常会出现变号现象.基于电阻率的传统反演方法无法对变号数据进行正确反演,因此通常在数据处理中予以剔除.为深入了解极化介质的电磁扩散特征,认识航空瞬变电磁负响应的产生机理,本文研究时间域航空电磁系统的电磁扩散特征.我们以均匀极化、非极化半空间及层状介质模型为例,通过直接积分的方法求解频率域电场响应,并由欧姆定律得到电流响应,再经过汉克尔变换得到时间域电流响应.通过研究电流随时间在地下极化介质中的传播特征研究电磁扩散过程;通过对比不同激电参数对电磁扩散的影响,研究极化介质中感应电流与极化电流的扩散规律,从而合理地解释极化介质中负响应的产生机理.基于本文研究和分析结果,可加深对时间域航空电磁法中激电效应的认识. 相似文献
2.
3.
地电体对频率域电磁波激发源的响应为电磁感应和激电效应的综合响应.传统CSAMT法进行数据正反演时认为大地介质电阻率是与频率无关的实数,而实际上因为激电效应,地下可极化体的电阻率是一个与频率相关的复数.为推进二者总体响应研究,并扩展激电法的应用范围,同时提高电磁法勘探的精度,本文基于Dias模型,以复电阻率代替不考虑地电体极化效应的直流电阻率,对CSAMT场源一维层状模型进行了正演模拟,为提取CSAMT信号中所含激电信息提供理论基础.结果表明,考虑激电参数后,视电阻率及相位响应曲线出现明显异常（包括远场、过渡场、近场）;极化前后振幅比值异常峰值、相位差值异常峰值可直观体现激电异常;异常峰值与极化层层厚、埋深以及电阻率变化有连续的对应关系.认为从频率域电磁法信号中提取激电信息有乐观的前景. 相似文献
4.
利用数值模式匹配技术,研究并建立层状各向异性倾斜地层中多分量感应测井响应的快速数值模拟算法.首先将位于井轴上的三个相互正交磁偶极子转化成关于极角θ的三个谐变分量的叠加,以便将三个正交磁偶极子电磁场的正演问题完全简化成三个谐变分量电磁场的轴对称定解问题,并给出电磁场各个谐变分量在井轴上满足的边界条件,保证电磁场在井轴附近仍然可解.然后利用数值模式匹配技术建立电磁场各个谐变分量的正演过程,得到层状各向异性地层中磁流源并矢Green函数的半解析表达式,给出计算层状各向异性倾斜地层中多分量感应测井响应的具体方法,最后通过数值计算结果证明该算法的有效性并考察几种不同情况下多分量感应测井响应特征. 相似文献
5.
利用数值模式匹配技术,研究并建立层状各向异性倾斜地层中多分量感应测井响应的快速数值模拟算法.首先将位于井轴上的三个相互正交磁偶极子转化成关于极角θ的三个谐变分量的叠加,以便将三个正交磁偶极子电磁场的正演问题完全简化成三个谐变分量电磁场的轴对称定解问题,并给出电磁场各个谐变分量在井轴上满足的边界条件,保证电磁场在井轴附近仍然可解.然后利用数值模式匹配技术建立电磁场各个谐变分量的正演过程,得到层状各向异性地层中磁流源并矢Green函数的半解析表达式,给出计算层状各向异性倾斜地层中多分量感应测井响应的具体方法,最后通过数值计算结果证明该算法的有效性并考察几种不同情况下多分量感应测井响应特征. 相似文献
6.
7.
8.
9.
常见弹性波动理论的建立是基于介质均匀这一基本假设,实际介质的非均匀性非常普遍.为研究连续介质中波的传播特征,本文从弹性力学中建立弹性波动方程的三个基本方程出发,考虑连续介质弹性参数的空变特征,建立非均匀介质的弹性波动方程,利用Alkhalifah声学近似思想建立位移表征的纵波波动方程,利用本征值问题求解方法建立标量波频率-波数域传播算子,从而建立描述纵波传播的标量波方程,其中波函数为纵波位移的散度,不同于均匀介质标量波方程的波函数为位移势.随后推导含PML边界波动方程差分格式并建立不同模型数值模拟进行数值试算,与均匀假设标量波方程和变密度方程对比证明本方法的准确性和稳定性. 相似文献
10.
《地球物理学进展》2015,(3)
海洋可控源电磁法(MCSEM)已日益成为海洋油气勘探的重要方法,但在传统MCSEM的数值模拟和反演解释中,通常考虑的是非极化海底介质模型,这与油气藏会产生激电效应的实际情况并不相符.在前人的研究基础上,本文总结了油气藏产生激电效应的两种模式:外生黄铁矿极化模式和油气藏自身极化模式.在典型的海洋4层地电模型基础上,本文设计了包含外生黄铁矿的海洋6层地电模型.基于这两种模型,通过推导MCSEM一维正演公式,引入Cole-Cole模型,本文分别计算了上述两种极化模式下MCSEM的激电响应,分析了各参数对两种模式下激电响应的影响特征,并对两种极化模式共同作用所产生的激电效应特征做了分析.结果表明,外生黄铁矿极化和油气层自身极化都会产生较明显的激电效应,而激电效应的强弱与观测频率、模型参数和激电参数有关.在相同的激电参数下,外生黄铁矿极化所产生激电效应的影响主要集中在小收发距处,而油气层自身极化所产生的激电效应则主要影响较大收发距处的场值.两种极化模式共同作用的效果相当于两种极化模式分别作用时的叠加. 相似文献
11.
12.
三维陆地可控源电磁法有限元快速正演的主要瓶颈在于多频率大型稀疏方程组求解问题.本文引入一种基于模型降阶的Krylov子空间投影算法,推导了有限元刚度矩阵的模型降阶形式,构建了频率域传递函数;采用标准正交向量序列,构建一个远远小于有限元刚度矩阵维度的矩阵,该矩阵与频率无关,通过一次模型降阶即可实现多频点有限元方程快速求解.采用基于电场的变分方程,加入散度校正条件,以消除伪解;引入伪δ函数,消除了源点的奇异性,可适用于复杂背景模型三维有限元数值模拟,并为多源的求解奠定了基础;以层状介质模型解析解为标准,通过和基于Pardiso直接求解器的有限元算法（3DFEM）进行比较,模型降阶法计算时间小于前者的1/10,平均相对误差在1.72%,在满足精度要求下,实现了高效率三维有限元数值求解;分别设计了横向高低阻模型和纵向高低阻模型,分析了从近区到远区电场和卡尼亚视电阻率的变化规律,假极值的表现特征,阴影效应的影响等,从而也验证了该算法的正确性.最后,建立了一个地层陷落柱模型,通过模型降阶有限元正演模拟,发现视电阻率断面图在陷落柱上方出现"凹陷",与模型设计吻合,表明该算法对复杂地层模拟具有同样的适用性. 相似文献
13.
考虑地球介质电导率任意各向异性且随空间位置连续变化的情况,本文实现了直接求解电磁场的可控源音频大地电磁测深（CSAMT）三维有限元数值模拟.首先给出了电导率任意各向异性介质中CSAMT二次电场满足的控制方程及其相应变分问题,然后采用任意六面体单元对研究区域进行剖分,在网格单元中对任意各向异性电导率进行线性插值,解决了实际工作中岩矿石电导率各向异性且连续变化的情况,将变分问题转化为线性代数方程组的求解.电导率各向异性且连续变化一维模型三维有限元数值模拟结果与电导率各向异性且分层均匀渐进模型解析解结果对比验证了方法的有效性;三维地电模型电导率随位置线性变化且各向同性、主轴各向异性、方位各向异性和倾斜各向异性的数值模拟结果表明,电导率各向异性且连续变化对CSAMT视电阻率和相位数据均有明显的影响. 相似文献
14.
全局弱式无单元法是在有限单元法基础上发展起来的一种数值模拟方法,它采用局部支持域内的节点信息来构造形函数实现局部精确逼近,摆脱了单元,仅依赖于节点信息,具有预处理简单、模拟精度高、灵活性强的特点,适用于复杂地电条件下直流电阻率正演模拟.本文采用RPIM构造直流电阻率全局弱式无单元法形函数,利用RPIM形函数推导了直流电阻率全局弱式无单元法方程.然后,编制了直流电阻率全局弱式无单元法正演模拟Fortran程序,利用该程序对典型的地电模型进行了正演模拟,并将正演结果与基于线性插值的FEM正演结果及解析解进行对比,结果表明采用RPIM形函数的全局弱式无单元法用于直流电阻率正演模拟的正确性及有效性,且在同等条件下,全局弱式无单元法模拟精度高于矩形剖分的FEM,更有利于指导电法勘探的数据解译;利用该程序对复杂地电模型进行了正演模拟,结果表明全局弱式无单元法对复杂地电模型模拟效果良好,适应性强,灵活性高,可任意加密节点提高模拟精度. 相似文献
15.
本文实现了2.5维电导率正交各向异性海洋可控源电磁等参有限元数值模拟.利用傅里叶变换导出了电导率正交各向异性2.5维海洋可控源电磁法波数域电磁场耦合方程,采用伽里金加权余量法推导了相应的有限元方程;采用任意四边形单元对研究区域进行剖分,在单元中进行双二次插值,将有限元方程化为线性代数方程组;最后,求解线性方程组并进行反傅里叶变换获得空间域电磁场值.这个方法可以模拟海底起伏地形条件下地下任意形状电导率正交各向异性的复杂模型.与一维模型的数值模拟结果对比表明,电磁场数值解与解析解吻合.二维模型的计算结果与二维自适应非结构有限元模拟结果也吻合.水平海底二维地电模型考察了不同各向异性系数对海洋可控源电磁响应的影响特征.海底起伏地形地电模型的数值结果表明,电导率各向异性对海洋可控源电磁响应影响明显,有可能淹没海底地形和高阻油气藏引起的异常. 相似文献
16.
本文用边界单元法模拟二维地形对大地电磁场的影响.首先用格林公式将二维大地电磁的边值问题转变成积分方程,然后用边界单元法解积分方程,得到地形上的大地电磁场和它的法向导数,由此可计算电阻率.与有限单元法相比,本方法剖分后的地形与实际地形的拟合程度高,向计算机输入原始数据的工作简单,可以在微机上计算地形引起的视电阻率.本方法适合于在野外生产现场进行大地电磁法的地形改正.计算表明,地形对H_x型波的影响比对E_x型波严重得多;随着周期的增长,地形对H_x型波的影响变得稳定. 相似文献
17.
复电阻率法在矿产、油气勘探调查中发挥着重要作用.为了认识复杂构造的复电阻率法电磁场的变化规律,本文基于自适应有限元方法,采用非结构化网格,引入Cole-Cole模型,实现了电偶源2.5D复电阻率法电磁场正演,可以模拟复杂地形和地电结构,正演结果更符合野外实际地质情况.通过将本文的计算结果与半空间模型解析解、层状介质和起伏模型结果进行对比,验证了本文算法的正确性.最后,基于复杂地电模型,通过正演模拟,系统分析了地形、激电参数、复杂构造对复电阻率法电磁场的影响特征. 相似文献
18.
点源三维地电场的求解是一大型数值计算问题．本文用有限差分方法求解,最后形成一个线性方程组Ax＝b,这里A是大型稀疏的带状对称矩阵．解大型稀疏方程组的完全Cholesky分解直接算法,一般要求巨大的机器内存来存储系数矩阵A,而且计算速度极慢．因此引入不完全Cholesky共轭梯度（ICCG）算法及按行索引的稀疏存储模式,充分利用A的稀疏性,使得计算速度大大提高,而内存要求则大大减少,因此ICCG算法是地电三维正演的强有力工具． 相似文献
19.
系统地论述了用有限单元法研究复杂地形条件下三维直流电阻率的正演计算技术.首先给出了三维构造中点源电场的边值问题以及相应的变分问题;然后利用有限单元法求解变分问题,采用四面体单元对研究区域进行剖分,在单元中进行三线性函数插值,将变分方程化为线性代数方程组;最后,考虑到节约计算时间,利用对称超松弛顸条件共轭梯度迭代算法求解大型线性方程组,得到了各节点的电位值,进而计算出地表的视电阻率.通过理论模型的计算检验了算法的可行性之后,给出了几种常见纯地形异常的数值模拟结果和一个组合模型的计算结果,其研究工作为研究三维直流电阻率反演奠定了基础. 相似文献
20.
谱元法已成为区域性乃至大陆性尺度地震波场模拟的重要工具.对于区域或大陆尺度层析成像而言,地球曲率不可忽略,此时模拟地震波传播采用球坐标系更为合适.本文从球坐标系下弹性波动方程弱形式出发,基于球坐标系变分原理给出了球坐标系下求解三维地震波方程的谱元法.另一方面,计算Fréchet敏感核是进行全波形反演的关键,本文借助伴随原理,推导了全波走时层析成像三维Fréchet敏感核表达式.为了验证球坐标系下谱元法的精度,我们将数值模拟结果与normal mode方法得到的解析解在1-D PREM模型下进行了对比.同时,我们将此方法应用到华北克拉通区域,以期获得地球内部结构精确成像.基于3-D全球径向各向异性地幔模型S362ANI和3-D地壳模型Crust1.0,我们建立了华北克拉通初始3-D背景模型,并将数值模拟结果与实际观测台站记录波形资料进行对比分析,利用互相关方法提取走时残差,最后给出了Fréchet敏感核在3-D空间中的分布,这些工作为下一步开展球坐标系下三维大尺度全波形反演奠定了基础. 相似文献