共查询到19条相似文献,搜索用时 156 毫秒
1.
由于受发射线圈电感和电容的影响,瞬变电磁发射电流类似斜阶跃关断,采用阶跃关断"单烟圈"理论进行视电阻率解释,在早、中期产生失真。根据斜阶跃关断二次场感应电压表达式,提出采用"双烟圈"理论进行解释:将完全关断后二次感应电压等效为大小相等、方向相反、先后发射的两个"单烟圈"二次磁感应强度之和。分别采用前向递推和后向递推两种方法,计算出每个烟圈二次磁感应强度,再利用"单烟圈"二次磁感应强度的单调性,计算斜阶跃关断全期视电阻率,消除了失真。最后通过数值模型,验证了方法的正确性,为全程瞬变电磁探测及解释提供理论支持。 相似文献
2.
为了研究复杂地电模型的航空瞬变电磁法全波形响应特征,需要开发考虑发射波形的三维数值模拟算法。本研究基于非结构四面体网格和位移逆Krylov子空间(Shift-and-Invert Krylov,简称SAI Krylov)方法,采用基于电偶极子离散的场源处理方法模拟场源,在时间域进行计算实现了全波形航空瞬变电磁法矢量有限元三维数值模拟。使用均匀半空间模型在阶跃波、半正弦波、三角波和梯形波激发下的全波形解析解、VTEM实际激发波形的后推欧拉算法计算结果,检验了本研究开发的数值模拟算法的正确性。设计地表起伏异常体模型,计算和分析了航空瞬变电磁响应特征。开发的基于位移逆Krylov子空间的全波形航空瞬变电磁法三维数值模拟算法适合模拟复杂地电模型的响应,具有较高的计算精度。 相似文献
3.
首先给出瞬变电磁(TEM)场卷积计算式,并将关断电流变化率函数看成冲击函数,垂直阶跃脉冲激发的场分量看成输入信号,卷积结果就是任意形状关断电流的TEM响应。把电流激发脉冲看成是许多垂直阶跃函数的叠加,通过叠加这些阶跃函数的响应求得任意电流脉冲激发的瞬变电磁响应。当已知TEM场的实测值和关断电流形状,就可按编程的方法计算全期视电阻率。文中给出了2个斜坡关断电流激发的电偶极源TEM垂直磁场算例,一个为均匀半空间的响应及其全期视电阻率,另一个为三层K型地电断面的响应及其全期视电阻率。 相似文献
4.
地?井瞬变电磁法响应规律复杂,现有解释方法以定性分析和半定量解释应用最为广泛,不能直接获取大地电阻率参数。针对这一问题,提出一种基于瞬变冲激时刻的快速定量解释方法。首先给出均匀半空间地?井瞬变电磁响应的表达式,分析地?井瞬变电磁响应的冲激时刻特征。结果表明,接收点深度越大、大地电导率越高,则瞬变冲激时刻越晚。结合已有的研究成果,推导冲激时刻与大地电导率和深度的函数关系,依据反函数理论进行大地视电阻率定义。以获取真实大地电阻率为目标,研究基于地下电磁场扩散速度的改进大地电阻率恢复算法。采用所提出的算法,根据实际常用工作方式,分别设计均匀半空间、二层模型和三层模型进行模拟计算。模型算例和实测数据试算结果表明:基于冲激时刻的视电阻率定义方法能够较好地反映大地电阻率的变化趋势,但具有较强的体积效应;基于电磁场扩散速度的改进算法能够有效地削弱体积效应的影响,更加准确地恢复大地电阻率值和反映电性界面。该算法无需进行复杂模型的迭代正演计算,具有较高的计算效率,能够定量恢复大地电阻率值,适用于地?井瞬变电磁法的快速初步定量解释。但在实际资料解释应用中,还需考虑视电阻率的“overshoot”和“undershoot”现象,避免造成错误解释。 相似文献
5.
《物探化探计算技术》2019,(5)
由于在均匀半空间模型的垂直感应电动势与模型电导率之间存在"二值性",这对进行视电导率求取带来了一定的困难。以中心回线式直升机航空瞬变电磁法为例研究电导率深度成像(CDI),利用分段二分法改进视电导率搜索从而找到一个更符合理论的视电导率值,然后根据扩散深度求得各层的视深度,最后对成像深度进行微校正后建立CDI成像结果。用层状模型进行验算结果表明,该方法对低阻薄层,特别是浅部成像有很好的效果,该方法简单易行且成像效果较为理想。 相似文献
6.
7.
多数地层和裂隙积水区的电导率具有各向异性特性,地层的各向异性对瞬变电磁观测会产生较大的影响,尤其是对水平分量特征更加明显。为研究轴向各向异性地层的瞬变电磁三分量响应特征,基于时域有限差分算法,引入电导率张量构建控制方程,实现电导率轴向各向异性三维瞬变电磁三分量正演。通过与各向同性半空间模型和各向异性半空间模型的一维解析结果对比,验证了算法的准确性;建立各向异性半空间模型、层状模型、含水体模型并计算其回线源瞬变电磁三分量磁场响应。结果表明:水平方向电导率各向异性对三分量磁场响应均产生较大影响,其中x轴各向异性对?By/?t分量响应的影响大于?Bx/?t,y轴各向异性对?Bx/?t分量响应的影响大于?By/?t,z轴各向异性对三分量磁场响应几乎没有影响;浅部地层的电导率各向异性对三分量磁场响应的影响占主导成分;当采集点在x方向距异常体中心小于y方向时,x轴各向异性对三分量响应影响程度大于y轴各向异性,反之亦然。基于本文研究成果,可为各向异性地层条件下瞬变电磁三分量处理解释提供一些有价值的理论借鉴。 相似文献
8.
9.
10.
11.
介绍了层状大地上水平共平面共心线框在阶跃函数或矩形脉冲激发时使用Gaver-Stehfest算法计算互阻抗的方法以及斜坡函数断电效应的修正、给出了重叠回线装置时全时域视电阻率的计算方法和TEM资料的反演算法。 相似文献
12.
13.
时间域航空电磁系统探测深度与采样时间、发射磁矩、大地电导率、仪器背景噪声和灵敏度等诸多因素相关。本文基于时间域航空电磁一维正演对时间域航空电磁系统探测深度进行研究。正演算法从麦克斯韦方程出发,结合准静态近似条件,得到一维层状介质上空中心回线频率域电磁场响应的垂直分量,并采用汉克尔积分进行计算,进而通过时-频变换获得时间域电磁场响应。探测深度研究参考仪器背景噪声水平,通过设定最小可识别信号阈值(本文设定为系统背景噪声三倍)获取最大可探测分离时间,进而利用该分离时间结合平均电导率求取最大探测深度。最后本文以VTEM系统为例,对不同层状介质模型以及不同发射磁矩、背景噪声、发射脉宽和飞行高度下的时间域航空电磁系统最大探测深度进行分析和讨论。 相似文献
14.
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… 相似文献
15.
三维地形大地电磁场的边界元模拟方法 总被引:1,自引:0,他引:1
提出了一种用边界元法计算大地电磁场三维地形影响的数值模拟方法.首先用矢量积分理论和电磁场边界条件, 将上半空间(空气)和下半空间(地下介质)两个区域电磁场边值问题变为仅对地形界面的两个矢量面积分方程, 其中一个计算磁场, 称磁场方程; 另一个计算电场, 称电场方程.然后将对地形界面的积分剖分为一系列的三角单元积分.在三角单元积分中, 假设单元中电磁场为水平均匀大地空间电磁场与地形影响的迭加, 并假设地形影响为常项, 这样既保证了计算精度又使得计算方法简便.通过分解和计算, 每一个矢量面积分方程分解为对应3个坐标方向的3个常量线性方程, 这些线性方程组成了对角占优的线性方程组, 可用SSOR方法求解.文中给出了2个三维地形上大地电磁视电阻率曲线的计算结果. 相似文献
16.
射频大地电磁法(Radio-magnetotelluric, RMT)是浅地表电磁勘探的重要手段之一。由于勘探频段为10~300 kHz,其电磁场传播受地下介质介电常数的影响较大。传统准静态条件下的电磁响应严重制约RMT正演模拟精度,并进一步影响反演成像分辨率。针对这一问题,提出了一种基于Julia并行加速的全电流RMT电磁响应数值模拟方法,利用Julia的分布式计算将各频点的计算发送到不同进程进行求解,从而达到加速计算的目的,同时在计算中考虑位移电流的影响,提升正演的模拟精度。通过计算几种典型高阻/高介电模型的RMT响应,分析并总结了位移电流对射频段电磁场视电阻率及相位响应的影响规律。数值模拟结果表明:当浅部存在高阻覆盖层时,基于准静态假设条件计算的RMT视电阻率和相位响应偏高,且频率越高、覆盖层电阻率越大,响应偏差越大;对于煤炭采空区模型,RMT法能有效反映异常体位置,但忽略位移电流会在采空区及其附近引起较大计算误差;起伏地形的算例表明地形会覆盖地下异常体的RMT数值响应,尤其是地形拐角处;2种不同规模的并行算例对比证明了并行算法的高效性,且随着求解问题规模增大,并行算法效率也随之... 相似文献
17.
以SHAKE2000为代表的频域等效线性化是目前主流的土层地震动计算方法。由于该方法在软土场地计算结果严重不合理,其改进方法成为目前研究热点,主要采用频率相关法,但一直没有取得实效。对频率相关法中土体动剪应变和震动速度成恒定比例关系的基本假定进行研究,从波动方程推导出了均匀全空间单向行波、单一均匀半无限场地、成层场地等3种典型场地中土体剪应变与震动速度关系的精确解答,通过数值试验研究了这一假定的合理性与偏差程度。结果表明:动剪应变和震动速度成恒定比例关系这一假定只有在无边界均匀介质单向行波情形下才成立,对于实际成层土场地,动剪应变与震动速度关系强烈依赖于波的频率和观察点的位置,如果在土层地震反应分析中忽略反射波而采用单向行波恒定比例假定将会使计算结果产生显著偏差;动剪应变与震动速度成恒定比例关系的假定导致的偏差达4个量级,即使对于单一均匀半无限空间模型偏差也十分显著;对实际土层地震动的求解,该假定在理论上定性错误,定量上偏差不可接受,必须放弃而另寻其他途径。 相似文献
18.
为了更好地模拟地下介质连续变化及开展连续介质的反演,对二维电导率分块线性变化的线源频率域可控源电磁法进行了有限元正演模拟,在剖分单元内同时对电场及电导率参数线性插值,使电导率参数在剖分单元之间保持连续变化。首先,提出有限元正演模拟的边值问题及变分问题,并详细论述了有限元的剖分、插值、单元分析及总体合成的各个步骤;其次,采用稀疏存储及基于不完全LU分解的BICGSTAB算法求解复系数方程组,节省了内存并提高了计算速度;然后,对一个均匀半空间模型进行模拟,计算结果表明,低频及高频的有限元数值解都与解析解吻合,证明了算法的正确性;最后,对水平层状模型及垂直断层模型进行正演计算,视电阻率及相位的等值线图均较好地反映出了异常体,说明文中算法能够对电导率连续变化的线源可控源电磁法进行有效地模拟。 相似文献
19.
巷?孔瞬变电磁法在存在明显电导率各向异性的勘探区会产生较大的解释误差。基于时域有限差分算法,通过引入电导率各向异性张量构建控制方程、将矩形回线源电流密度加入Maxwell方程安培环路定理实现任意电流源的加载,以差分代替微分对控制方程进行离散,实现巷?孔瞬变电磁三维模型正演计算。在与解析解对比验证算法计算精度的基础上,构建全空间模型、层状模型和三维块状模型进行正演并分析电导率各向异性对巷?孔瞬变电磁三分量响应的影响程度与方式。结果表明:垂直轴电导率对巷?孔瞬变电磁三分量响应基本没有影响,水平轴电导率对巷?孔瞬变电磁三分量响应影响较大,其中,?Bx/?t响应主要受y轴电导率影响,?By/?t响应主要受x轴电导率影响;通过三分量响应的形态和幅值的相互关系可以辨别各向异性介质所在方位及主轴电导率方向;当异常体所处全空间介质为电导率各向异性时,异常体产生的异常响应会被全空间介质电导率各向异性产生的异常响应所淹没。电导率轴向各向异性特征在巷?孔瞬变电磁法解释过程中不可忽略,研究成果为巷?孔瞬变电磁法各向异性解释提供指导,也为各向异性反演提供参考。 相似文献