共查询到19条相似文献，搜索用时 140 毫秒
1.
斜阶跃波层状介质航空瞬变电磁响应数值计算 总被引：1，自引：0，他引：1
目前,航空瞬变电磁方法数据解释主要采用阶跃波均匀半空间模型计算视电导率值,而实际航空电磁系统发射波形的下降沿多为斜阶跃,导致解释时计算的视电导率值存在较大偏差。为此,笔者研究了航空瞬变电磁系统发射电流为斜阶跃时的电磁响应数值计算,将发射波形进行拉氏变换,利用G-S逆拉氏变换与241点汉克尔变换相结合的方法,实现斜阶跃波关断后的层状大地模型航空瞬变电磁响应计算;并对均匀半空间和层状大地模型下,不同关断时间和不同飞行高度对电磁响应的影响进行分析。得出结论:不同关断时间,关断后取样延时2 ms时,均匀半空间电磁响应的平均偏差为27.78%,三层模型的平均偏差为32.16%;当飞行高度从20 m增加到60 m时,均匀半空间和三层模型的感应电动势分别减小了43.6%和83.2%。 相似文献
2.
《物探与化探》2015,(5)
采用三角网格剖分的有限元法,研究了2.5维航空瞬变电磁法正演模拟问题。利用时频变换数值方法将时间域电磁场转换到拉氏域,再利用傅里叶变换将三维问题降维变为2.5维问题,然后由有限元法求解得到拉氏域二维电磁场,逆拉氏变换后得到时间域航空瞬变响应。为了回避正演模拟中总感应磁场在场源处的奇异性问题,采用异常场算法,场源响应通过在微分方程中施加背景电磁场实现。由于瞬变电磁信号具有较大的动态范围,而且需要经过两次正、逆拉氏变换和傅里叶变换,每个环节的计算精度和速度要严格控制在较高的水平上,否则积累误差会非常大。模型计算表明均匀大地和层状大地模型解析解与数值解吻合很好。这证明该算法是正确可行的,可作为研究二维复杂地质体的方法手段。 相似文献
3.
《物探与化探》2015,(6)
在地-井瞬变电磁正演模拟中,前人均以总场为研究对象,分析响应特征。此次选择地-井瞬变电磁异常场为研究对象,对单异常体和多异常体以及在"定源异井"和"动源定井"两种观测方式下典型地电模型进行地-井瞬变电磁法三维正演模拟。首先推导了二次磁场的扩散方程,并用有限差分法进行数值计算。通过分析以地-井瞬变电磁异常场为研究对象,认为异常场能反映目标体空间几何位置参数。研究结果表明:在不同钻孔中接收的异常响应随时间推移而衰减减弱,但异常变化特征保持不变;在不同方位布置场源得到的异常响应差异明显,根据其测量结果可判定异常体空间位置信息。最后总结影响地-井瞬变电磁异常的主要因素有:目标体的埋藏深度、目标体与钻孔之间的耦合关系以及场源的位置关系。研究结果为相关资料解释提供理论参考。 相似文献
4.
5.
《物探与化探》2021,(3)
Z轴倾子电磁法(ZTEM)是一种新型频率域航空天然场源电磁法,该方法采用倾子作为研究参数,可用于起伏地形下的大规模地质勘探。本文在对TEM三维有限差分正演和数据空间OCCAM反演算法研究的基础上,考虑地形起伏的影响,研究了带地形的频率域三维ZTEM正反演算法。首先,对起伏地表下ZTEM正演算法的准确性进行了验证,计算和分析了纯地形的ZTEM三维异常响应特征。其次,通过对山峰和山谷地形下低阻棱柱体合成算例的反演分析,表明带地形的ZTEM数据空间OCCAM反演算法能够获得比较接近真实导电性结构的地电模型,尤其对地下目标体的横向边界具有理想的约束效果。最后,与不带地形的ZTEM反演结果对比,验证了所开发的带地形ZTEM倾子资料反演算法的有效性。 相似文献
6.
积分方程法三维模拟井地电法并行算法研究 总被引：1，自引：1，他引：1
在积分方程法模拟三维电磁响应原理的基础上,对井地电法三维正演模拟并行算法做了研究。同时,还利用Fortran 6.5和消息传递接口(MPI)开发了井地电法三维正演模拟并行计算程序。在求解大型线性方程组时,采用的是稳定的双共轭梯度(BICGSTAB)算法,单机正演程序与并行正演程序的计算结果有着很好的一致性,这说明了并行算法的正确性,为进一步开发井地电法三维反演并行程序奠定了基础,也为井地电法三维正、反演程序实用化提供了技术支持。 相似文献
7.
《物探与化探》2020,(2)
时间域有限元算法已被广泛应用于航空电磁三维正演模拟当中,然而由于航空电磁测区面积大,且采样密集,造成正演计算量巨大,传统的串行算法已经无法满足计算效率要求,为此,开展了并行加速算法研究以解决计算效率不足的问题。基于航空电磁系统的影响范围有限,采用局部网格技术将计算任务划分成多个子网格,即每个发射源一套网格,各网格的正演计算相互独立,不存在数据依赖性,具有很好的可并行性;利用MPI技术对多个子网格正演任务进行分配,在各个进程上进行并行计算;针对每个正演子网格,在进行时间域有限元算法正演模拟过程中,采用OpenMP技术对单元矩阵进行并行计算。典型地电模型的数值模拟结果表明本文开发的MPI+OpenMP并行正演算法可以有效提高正演速度,最高加速比可达10倍。 相似文献
8.
瞬变电磁法是探测煤田采空区的一种有效方法。瞬变电磁拟地震处理解释方法研究日渐成为电磁勘探领域的热点并且取得了突破性进展,所以可以把瞬变电磁拟地震处理解释运用到煤田采空区的探测中来提高勘探精度。文中在瞬变电磁三维时域有限差分数值模拟两层含水采空区模型的基础上,运用预条件正则化共轭梯度算法实现了瞬变电磁扩散场到虚拟波场的稳定转化,并对虚拟波场进行了脉冲压缩处理;最后通过克希霍夫积分解,对虚拟波场进行延拓成像,实现了两层含水采空区模型的三维解释。 相似文献
9.
《物探化探计算技术》2016,(3)
对三维场源二维地电模型的正演计算称为2.5维数值模拟。这里从麦克斯韦方程组出发,1分别求解电磁场的一次场和二次场,将三维场源降为二维;2利用傅氏变换将空间域方程转化为波数域,应用有限元求解波数域电磁场方程,引入无限元解决无穷远边界的收敛问题;3根据电磁场值实虚部的曲线特征,按对数等间隔选取21个波数。编写代码计算均匀半空间与解析解结果对比,电阻率的均方相对误差均小于0.5%,证明其有效性。计算了三种地电模型的电磁场响应,对单一高、低阻体和高低阻组合体的模拟效果真实,异常中心位置基本吻合。结果证明,这里的CSAMT2.5-D正演算法可以模拟较为复杂的地电模型,并取得良好的效果。 相似文献
10.
郑建波李美艳 《物探化探计算技术》2023,(4):478-483
航空瞬变电磁法以其快速高效的优势已获得广泛应用,然而航空瞬变电磁采样密集,数据量巨大。为了实现航空瞬变电磁观测数据的快速解释,开展航空瞬变电磁快速成像算法研究。利用反函数定理建立成像迭代格式,并在成像过程中考虑发射源高度和观测时间的影响,进而实现航空瞬变电磁全时域全空域快速成像。这里首先给出了全时域全空域成像理论和实现方法,然后对航空瞬变电磁装置进行了介绍,最后建立了典型采空区模型,利用开发的航空瞬变电磁全时域全空域视电阻率成像方法进行数据成像,成像结果与真实模型基本一致,验证了方法的有效性。 相似文献
11.
依据电磁场理论，构建地电模型，利用有限差分对矿井三维地电异常体进行数值模拟，得到异常体在电阻率不同、至巷道距离不同、及测点不同情况下的层测深曲线，并对其分布规律进行分析研究。通过模拟矿井三维地质异常体地电异常响应特征，并与井下实测层测深曲线对比和分析，论证了矿井三维地电体数值模拟的正确性。 相似文献
12.
时间谱电阻率(TSR)法的二维正演, 是计算可极化二维地电构造上三维电流偶极源的电场瞬变响应, 属所谓2.5维时间域电磁场数值模拟问题, 是目前国际上未妥善解决的计算地球物理疑难问题.针对现有算法的局限性, 建立了新的算法, 其特点是: (1)采用二维有限单元算法, 在矩形网格中增加两对角线形成三角网格剖分, 同时用高斯消元法消除矩形网格中心结点的待求未知量.这样, 既可较准确地模拟任意二维复杂地电断面, 又可节省计算量.(2)采用直接计算二次场的新算法, 只需计算电场和磁场沿地电构造走向两个分量的一次场, 因而, 不但计算精度较高, 而且不显著增加计算量.(3)采用G -S变换法作逆拉氏变换, 并利用拉氏变换延迟定理在倍增的时间间隔中插值, 从而实现对密集采样时间瞬变过程的快速计算.(4)能对可极化和导电大地(即同时包括IP和EM效应). 相似文献
13.
在三维电阻率的正演计算中往往涉及到快速、准确求解大型线性方程纽Ax=b的问题。通过采用有限差分法来构造出求解点电源三维地电场的大型稀疏对称线性方程组。并引入Lanczos迭代技术，构造出三对角阵方程组，然后采用正交分解法进行求解，它是Krylov子空间方法中的一种。与传统迭代算法相比，它占用内存少，收敛速度快且稳定。针对大型稀疏矩阵及MATLAB语言的特点，采用简单记录矩阵的非零元素值及其所在行、列值的方法存储大型稀疏矩阵，可大大节省机器内存，提高运算速度。理论分析和计算实例显示，此算法是地电三维正演计算的有效方法，为下一步的反演计算打好基础。 相似文献
14.
为了研究瞬变电磁场分布变化规律，在理论上推导了不同阻尼状态下的瞬变电磁暂态响应方程，并采用Simulink对瞬变电磁场暂态过程进行模拟仿真，设计开发了激励电流、发射回线、地质异常体、接收回线等效电路模块和互感关系模块；系统研究了电感、分布电容和电阻对发射和接收回线等效电路的影响。结果表明:发射回线的阻尼状态决定了瞬变电磁场激励波形特点，接收回线的阻尼状态则直接影响瞬变电磁场的衰减特征和观测质量。电路仿真结果对于抑制非地质因素的影响和指导瞬变电磁观测系统设计具有重要意义。 相似文献
15.
采用易于模拟复杂地形起伏和倾斜界面的非结构三角单元剖分网格,并利用对偶加权后验误差估计指导网格自动细化过程,实现了2.5维直流电阻率法自适应有限元数值模拟。在实例模型分析中,分别计算了层状模型和垂直岩脉模型的直流电阻率响应,并与其解析解进行了比较。对比结果表明,该算法所得数值解精度很高,解的相对误差小于0.5%。最后,计算了起伏地形2.5维地电模型视电阻率异常,并利用比较法进行了地形改正。地形改正结果与水平地形时的结果对比表明,比较法可以较好地消除地形影响,突出局部地质体异常。 相似文献
16.
首先给出瞬变电磁（TEM）场卷积计算式,并将关断电流变化率函数看成冲击函数,垂直阶跃脉冲激发的场分量看成输入信号,卷积结果就是任意形状关断电流的TEM响应。把电流激发脉冲看成是许多垂直阶跃函数的叠加,通过叠加这些阶跃函数的响应求得任意电流脉冲激发的瞬变电磁响应。当已知TEM场的实测值和关断电流形状,就可按编程的方法计算全期视电阻率。文中给出了2个斜坡关断电流激发的电偶极源TEM垂直磁场算例,一个为均匀半空间的响应及其全期视电阻率,另一个为三层K型地电断面的响应及其全期视电阻率。 相似文献
17.
本研究实现了考虑发射电流波形的瞬变电磁一维正反演算法。正演采用矩形回线源,通过Duhamel积分计算梯形波的感应电压;反演采用Occam反演方案,得到一个光滑的层状反演结果模型。理论模型合成算例的计算结果表明：上升沿时间、持续时间和下降沿时间不同程度地影响了垂直磁感应强度对时间的偏导数(dB/dt)响应以及电阻率反演结果。理论模型的正演结果显示,上升沿时间、持续时间和下降沿时间分别主要影响晚期、中晚期和早期的瞬变电磁响应,并且取值越大影响越大。理论模型合成数据的反演算例显示,不同上升沿时间对反演结果几乎没有影响,而下降沿时间和持续时间分别影响浅部和深部的反演结果。 相似文献
18.
19.
地–空瞬变电磁法在煤炭采空区勘探等领域受到了越来越多的关注, 有必要研究发射电流波形参数对地–空瞬变电磁三分量响应特征影响, 为地–空瞬变电磁数据处理与解释提供理论依据。以梯形波为例, 首先, 研究不同发射电流波形的频谱分布情况; 然后, 基于三维时域有限差分正演研究发射波形的上升沿时间、脉宽和关断时间对地–空瞬变电磁三分量磁场响应的影响。结果表明:上升沿时间对三分量二次场响应基本不产生影响; 关断时间对三分量二次场响应的影响主要集中在0.2 ms之前, 且关断时间越长对纯异常响应影响越大; 脉宽对三分量二次场响应的影响主要集中在0.1 ms之后, 且脉宽越短对纯异常响应影响越大。三维采空区模型结果表明:关断时间、脉宽对三分量异常场和背景场响应的影响特征基本一致; 通过三分量纯异常场响应多测道图和时间道图可以判断异常体的分布范围和深度。研究成果可为地–空瞬变电磁激励源波形的参数选取提供有价值的理论借鉴。 相似文献