基于主轴各向异性介质的三维张量CSAMT正演模拟

在地质情况复杂化以及地质勘探精细化的要求下,有必要开展张量CSAMT正反演研究.考虑到各向异性对正反演结果影响较大,本文开展基于主轴各向异性介质的张量CSAMT正演研究.根据含源的电场控制方程,引入主轴各向异性介质的张量电阻率,将电场总场分解为一次场和二次场,一次场由Key的算法得到,二次场利用有限差分法计算,即实现了...     

各向异性弹性波场的有限元数值模拟

研究了各向异性弹性波有限元正演系统的精度和效率问题，提出了一种透射加衰减的组合人工边界方案（吸收边办条件），它对各向异性波具有较好的吸收效果，并且有较好的稳定性能，均匀TI介质中的模拟获得了非常清晰的波场快照，其波场特征与理论分析能够准确吻合，各向异性介质模型的地表地震记录表明，各向异性波炮集记录在波的类型、同相轴形态、能量分布和相位等方面与各向同性波都有很大差别。     

三维三分量CSAMT法有限元正演模拟研究初探

首先从麦克斯韦方程出发,用伽里金方法推导了三维三分量CSAMT法的有限元方程.在研究过程中,认识到加入散度条件的必要性,在公式中强加了散度条件,提高了解的完备性.其次将成功应用于二维线源频率域电磁法有限元模拟中的两种技术推广到三维中,一是边界条件统一采用一阶吸收边界,使线源产生的电磁波在边界上按波的传播规律被吸收,以降低平面波假设造成的影响;二是总体系数矩阵的存储,用两个二维数组分别记录总体系数矩阵的非零元素及其在总体结点编号中所处的位置,使总体系数矩阵的存储量达到最小的同时,物理意义明确,迭代求解时迅速简便.最后用均匀半空间模型进行了验证.     

基于非结构网格有限元的三维井地电阻率法任意各向异性正演模拟

地层介质的电各向异性增加了井地电阻率法响应的复杂性, 开展基于电各向异性介质模型的井地电阻率法响应规律研究对于正确解释各向异性显著地区的观测数据特征至关重要.针对垂直线源井地电阻率法的任意各向异性响应模拟问题, 本文提出了一种基于非结构网格有限元三维正演算法, 通过引入3×3的对称正定张量来表征任意各向异性的电导率, 采用非结构四面体网格有限元方法来离散电位的边值问题, 通过将垂直线源等效为一系列点源问题, 进而实现了任意各向异性介质中井地电阻率法的高效数值计算.通过与三个地电模型的解的对比, 验证了本文数值解算法的精度和有效性.针对线源远离和垂直穿过异常体的两类模型, 分别考察了当围岩或异常体为各向异性介质时的井地视电阻率响应特征.结果表明, 对于各向异性地层, 围岩和异常体的主轴电阻率值和旋转角均会对井地视电阻率的幅值及分布产生显著的影响.研究结果对于提高井地电阻率法的认识和资料解释水平具有重要的理论和实际意义.

     

任意各向异性介质中三维可控源音频大地电磁正演模拟

在一些地层层理发育的地区,地下介质存在显著的电各向异性,此时基于各向同性模型解释含各向异性效应的可控源音频大地电磁（CSAMT）测深观测数据会导致错误的结果.本文通过引入3×3的对称正定张量表征电导率各向异性,采用非结构四面体网格和矢量有限元方法离散电场满足的矢量Helmholtz方程,并将电磁场源等效为系列电偶极子,实现任意各向异性介质中CSAMT高效数值模拟.本文首先通过层状各向异性模型检验三维有限元算法的精度和有效性,进一步建立三维地电模型研究异常体各向异性和围岩各向异性对CSAMT响应的影响,最后使用视电阻率极性图来识别各向异性电导率主轴方向.数值模拟结果表明,各向异性电导率对CSAMT视电阻率幅值及分布规律都有很大影响,视电阻率极性图能够很好地识别各向异性主轴方向.

     

电导率任意各向异性海洋可控源电磁三维矢量有限元数值模拟

现有海洋可控源电磁三维数值模拟方法大多基于电导率各向同性介质理论,不能模拟海底地层电导率各向异性的实际情况.本文给出了电导率各向异性三维介质中电性源海洋可控源电磁二次电场的边值问题以及相应的变分问题,采用长方体单元对研究区域剖分,将场分量定义在剖分单元的边上,利用矢量有限单元法求解变分问题,实现了电导率任意各向异性海洋可控源电磁三维矢量有限元数值模拟.这个新的正演方法可以计算电导率任意各向异性三维地电模型的海洋可控源电磁响应,基于二次场矢量有限元法直接求解电磁场,避免了传统有限元方法可能遇到的伪解问题和难于处理电场法向分量不连续的问题,提高了数值模拟计算精度.一维电导率各向异性模型电磁场数值解与解析解吻合得相当好,无论在源附近还是远离源处相对误差均不超过1%.电导率各向异性二维模型的计算结果与已有文献采用的非结构有限元模拟结果十分吻合.三维地电模型数值模拟结果显示,电导率各向异性张量电导率主轴分量和欧拉角对不同装置海洋可控源电磁响应均有着明显的影响.

     

CSAMT三维交错采样有限差分数值模拟

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

某煤矿典型CSAMT法视电阻率曲线的一维模拟

了解煤矿开采区的地层地质情况,指明断层、破碎带、陷落柱等构造及其含水性,并分析其对开采过程中可能造成的安全威胁是煤矿安全生产的重要保证．可控源音频大地电磁测深(CSAMT)法是浅层地球物理勘探的一种有力的手段,可以用来探明从地表到地下1500米范围的地质情况．在开采区完成三维地震勘探,对构造有了较系统的控制后,利用电法CSAMT手段重点来分析构造及地层的含水及富水特性,是一种科学的工作流程．本文采用一维正演方法对该矿区典型的CSAMT视电阻率曲线进行模拟并与实际曲线进行了对比．通过大量的模型拟合,得到了较好的拟合效果,这对反演解释提供了有用信息。     

大地电磁二维对称各向异性介质的有限元数值模拟

假定垂直轴为二维对称各向异性介质主轴之一，构造走向与另一主轴方向成任一夹角，用伽勒金（Ｇａｌｅｒｋｉｎ）有限元法和矩形网格，优化地合成总体刚度矩阵，使二维各向异性介质的基本方程形成有限元代数方程组，求出各节点场值，并利用ＭＯＭ法求出辅助场，进而求出张量阻抗等响应函数．同时还对有关文献的计算模型进行了数值模拟检验     

二维有限元大地电磁正演模拟在Matlab上的实现

Matlab高级语言开发平台对矩阵的计算简便快速,并且绘图输出文件方便。以这些特点为契机,在Matlab平台上,讨论利用有限单元法实现二维大地电磁正演的问题,推导出相应的公式,计算地表或海底视电阻率及相位,并与解析解以及主流方法结果对比,证实了本文方法的正确性。     

有限元法2.5维CSAMT数值模拟

目前国内大多数CSAMT资料处理都用一维反演方法,即假定地球介质是水平层状地电结构,这对于横向变化较大的介质结构的反演则是不真实的．为了合理的模拟可控信号源的三维特征应该使用三维有限元方法．受计算机内存和计算速度的限制,三维方法的实用化受到约束．在许多情况下,地电结构沿走向变化很小,只沿倾向发生变化．这种地电结构是二维的,而人工源是三维的,因此CSAMT资料的观测可用2．5维有限元方法进行数值模拟．本文从麦克斯韦方程组出发,建立了2．5维有限元CSAMT数值模拟方法,其核心是把地电参数变化小的走向方向转化成波数域,用一系列波数模拟三维源的特征．用一个横向均匀的三层地电结构模型展示了2．5维数值模拟的特征,并和一维模拟结果进行了比较,证实了有限元方法2．5维CSAMT数值模拟的可靠性．在此基础上对一个地电结构已知的实际模型进行了2．5维数值模拟并和该剖面的野外实测剖面数据进行了比较,进一步有效地说明了本文介绍的2．5维CSAMT数值模拟方法是仿真和可靠的,为在此基础上的2．5维反演打下了良好的基础．     

双轴各向异性介质中回线源瞬变电磁三维拟态有限体积正演算法

采用模拟离散的有限体积法实现了双轴各向异性地层回线源瞬变电磁三维正演.首先引入内积定义,采用自然边界条件,将瞬变电磁法的控制方程转化为弱形式表示.将计算区域划分为一系列的控制体积单元,采用交错网格对控制方程进行模拟有限体积空间离散,包括旋度算子离散和空间内积离散.基于斯托克斯定理的旋度积分定义公式实现旋度算子离散.中点平均实现电导率双轴各向异性的空间内积离散,从而得到离散化的控制方程.时间步迭代采用无条件稳定的欧拉后向差分格式.并通过均匀全空间中稳定电流回线源的磁场解析表达式得到回线源初始时刻的电磁场分布.为了同时保证计算精度和效率,本文采用分段等间隔的时间步迭代,利用直接法求解器PARDISO实现其快速求解.最后通过对比层状模型和各向异性半空间模型的正演计算结果,验证了本文算法的计算精度和计算效率;计算三维双轴各向异性模型的正演响应可知,水平方向电导率变化对电磁响应产生显著影响,而垂直方向的电导率变化对电磁响应几乎没有影响.产生这一现象的主要原因是回线源产生的感应电流主要是水平方向的,因此响应主要受到水平方向电导率的影响,垂直方向的电导率影响很小.

     

三维CSAMT法非结构化网格有限元数值模拟

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

     

频率域线源大地电磁法有限元正演模拟

本文介绍了频率域线源大地电磁法有限元正演模拟的研究结果. 在外边界上统一应用适合于人工源的一阶吸收边界条件来形成边值问题,可减小基于平面波假设造成的人为截断边界的影响. 程序编辑中设计了两个二维数组分别存储总体系数矩阵的非零元素和在总体结点编号中的位置,使内存占用量减少,且物理意义明确,方便用高斯-赛德尔等迭代法解有限元方程时调用. 采用视δ函数模拟线源,提高了解方程组的稳定性. 最后通过对1个简单模型和1个复杂模型的模拟,证明所用的方法对异常体能够有明显的反映,说明了该方法的可靠性和有效性.     

用矢量有限元方法模拟随钻测井仪在倾斜各向异性地层中的电磁响应

用矢量有限元法模拟和分析了电各向异性倾斜沉积岩层的随钻测井电磁响应,为很好地模拟钻铤、井眼和激励线圈等,采用了基于圆柱坐标下的六面体剖分;同时采用了高阶叠层矢量基函数来描述每个剖分单元内的场,有利于减少剖分单元数和未知量.三维有限元程序采用数值模式匹配(当处于简单地层时)以及时域有限差分法计算的结果进行了验证,显示吻合非常好.最后仿真了不同各向异性系数和不同倾角的各向异性地层的随钻电阻率测井响应,有助于对各向异性及地层倾角的校正,以及对储层作出正确评价.     

复杂介质有限元法2.5维可控源音频大地电磁法数值模拟

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

起伏地形下CSAMT二维正反演研究与应用

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