共查询到18条相似文献，搜索用时 343 毫秒
1.
为了更好地模拟地下介质连续变化及开展连续介质的反演,对二维电导率分块线性变化的线源频率域可控源电磁法进行了有限元正演模拟,在剖分单元内同时对电场及电导率参数线性插值,使电导率参数在剖分单元之间保持连续变化。首先,提出有限元正演模拟的边值问题及变分问题,并详细论述了有限元的剖分、插值、单元分析及总体合成的各个步骤;其次,采用稀疏存储及基于不完全LU分解的BICGSTAB算法求解复系数方程组,节省了内存并提高了计算速度;然后,对一个均匀半空间模型进行模拟,计算结果表明,低频及高频的有限元数值解都与解析解吻合,证明了算法的正确性;最后,对水平层状模型及垂直断层模型进行正演计算,视电阻率及相位的等值线图均较好地反映出了异常体,说明文中算法能够对电导率连续变化的线源可控源电磁法进行有效地模拟。 相似文献
2.
岩矿石在形成过程中受湿度、含水率、温度、压力等多种外界因素的影响,这些因素会导致岩矿石的电导率连续变化,于是电导率连续变化地电模型可能更符合实际地质情况。对电导率连续变化二维介质大地电磁场分布规律,进行了有限元正演数值模拟研究。基于自适应有限元模拟技术和非结构网格剖分技术,使得新的算法可以真实地模拟电导率连续变化的复杂二维地电模型,并能够提供高精度的数值解。算例表明,海底近地表垂向电导率发生变化能够对高频大地电磁场响应产生较大影响。 相似文献
3.
利用有限元方法进行大地电磁正演数值模拟时,由于是在有限网格区域上的数值计算,模拟计算时的网格边界为截断边界,而有限元数值模拟时的大地电磁场边界条件需要在足够远处才能够满足,所以截断边界的存在可能会使大地电磁正演模拟的边界条件无法满足,致使对计算结果和计算精度产生影响。利用有限元二维正演程序,在网格边界处加载一维情况下的大地电磁场,然后固定研究区域的网格剖分,并对一维地电模型和二维地电模型在改变有限元网格边界大小的情况下进行计算。在对一维模型进行模拟计算时,截断边界对边界条件没有影响,边界条件自然满足。而对二维模型进行模拟计算时,截断边界的存在对计算结果有较大影响。利用趋肤深度作为有限元网格边界变化的量度,通过改变网格边界大小,对不同的二维地电模型进行计算比较,总结出适合大地电磁有限元正演模拟的参考网格边界。 相似文献
4.
《物探化探计算技术》2019,(5)
有限元求解大地电磁正演问题时,对研究区域的剖分常规做法是采用规则化的网格。但规则化网格在剖分地形、断层、褶皱等复杂模型会产生较大的几何离散误差。针对上述情况,这里采用非结构化的四边形网格对二维地电模型进行网格剖分,并与自适应有限元相结合,由剖分的粗网格出发,利用每个频点下网格单元的后验误差估计值指导网格的局部加密,优化网格质量和数量,从而提高正演模拟的精度。通过一维K型地电模型利用本文方法算出数值解与解析解进行对比分析,验证了自适应有限元法在求解大地电磁二维正演模拟中的有效性。并通过对断层、褶皱模型的正演模拟,分析了其大地电磁正演响应特征。 相似文献
5.
目前有限元技术的开发及在电磁勘查技术中的应用已日趋成熟.然而,有限元正演模拟主要集中于地面和海洋电磁,航空电磁三维有限元模拟尚没有受到足够重视.以前人工作为基础,利用结构化网格实现了航空电磁系统的三维矢量有限元正演模拟.从二次场双旋度矢量非齐次亥姆霍兹方程出发,应用广义变分原理推导出变分方程,并采用六面体单元剖分,将场置于单元棱边上,对每个单元应用线性插值,最后合成含有稀疏矩阵的线性方程组.针对航空电磁多源性问题,利用MUMPS(multifrontal massively parallel sparse direct solver)直接求解器进行求解,在保证精度的前提下大幅度提高计算速度.利用单个异常体模型检验算法的精度和软件的稳定性,进而通过典型地电模型的模拟验证算法的有效性.对不同地下电性结构正演模拟结果进行对比分析,进一步研究了覆盖层和垂直接触带等典型构造对航空电磁响应的影响特征. 相似文献
6.
高频大地电磁法有限元数值模拟 总被引：1，自引：0，他引：1
由于中深度工程地球物理勘探的迫切需要,以美国EH-4电导率成像系统为代表的高频率大地电磁法,在我国地球物理勘探行业的应用越来越广泛。其勘探深度在地下1000m以内,频率范围在10Hz～100kHz。这里研究了高频率大地电磁法有限元高精度,快速数值模拟。首先利用广义变分原理推导出电磁场的有限元方程,编制了双线性和双二次插值有限元程序,并进行了层状模型的验证性计算,然后对典型地堑模型的数值模拟结果进行了分析,得出了一些有意义的结论。 相似文献
7.
为了加快大地电磁三维正演的求解速度,本文将一种新型的代数多重网格算法——聚集多重网格（aggregation-based algebraic multigrid, AGMG）算法引入大地电磁三维正演模拟中。首先从准静态条件下的麦克斯韦方程出发,利用交错网格有限体积法进行离散,并采用第一类Dirichlet边界条件形成大型稀疏复线性方程组;然后阐述AGMG算法的粗化策略和套迭代技术,并实施3种不同的AGMG求解算法：1）传统的V循环AGMG算法;2）AGMG预处理共轭梯度（AGMG-CG）;3）AGMG预处理广义共轭残差法（AGMG-GCR）。最终实现大地电磁法三维正演模拟。对典型地电模型进行正演模拟,并与已有的大地电磁三维正反演程序（ModEM）进行结果对比,以验证本文算法的准确性。另外,不同剖分网格和极化方式正演模拟结果与准残量最小化（QMR）迭代算法的对比表明,AGMG预处理求解算法（AGMG-CG、AGMG-GCR）不仅能够改善算法的稳定性,而且能够快速有效地求解正演问题;其中AGMG-GCR迭代次数更少,求解速度更快,误差衰减曲线更光滑,在144×152×104网格剖分情况下,相对于现有ModEM程序能够提高十几倍的计算速度,尤其适合大规模大地电磁三维正演问题。 相似文献
8.
电阻率/激发极化率数据的二维反演程序 总被引：25，自引：2，他引：23
本程序在目前日本流行的电阻率／激发极化法二维反演程序［４］基础上,将电导率和极化率分块均匀变为分块连续变化,并在目标函数中加入了更多的先验信息,使反演结果的唯一性和真实性大大增加。在有限元正演模拟中,又采用三角单元剖分,使实测数据在反演以前不需要进行地形改正。本文介绍了该方法的基本原理和程序的使用方法,并给出了几个理论算例和野外实例 相似文献
9.
电各向异性在自然界中普遍存在,特别是沉积盆地中的部分岩层经过压实变质作用,表现出很强的各向异性导电性,采用基于各向同性模型对实测资料进行正反演解释必然会造成困难甚至结果的错误,只有基于各向异性理论的正反演解释才更为合理准确。本文通过自适应有限元对大地电磁各向异性进行正演计算,分别模拟水平、垂直和倾斜各向异性介质在不同偏转角和主轴电阻率下的响应结果。结果表明：自适应有限元能够在后验误差的控制下得到合理的网格,使计算结果更加接近解析解;在各向异性介质中,大地电磁TE极化模式的视电阻率和阻抗相位与垂直于层面的电阻率无关;二维电各向异性结构中,大地电磁TM极化模式响应结果总是由主轴上的电阻率在y轴方向上的分量所决定。 相似文献
10.
目前大地电磁测深法的一维正演理论,是基于平面波垂直入射水平层状介质的假设模型。但由于地球是一个球体,因此有必要研究基于球状介质模型的大地电磁正演理论。这里详细推导了基于球体层状介质模型的大地电磁正演公式,计算了若干理论模型。通过同基于水平层状介质模型的大地电磁正演结果对比,验证了正演公式的正确性。同时,指出当探测周期增加到上万秒时,阻抗相位会增大;而当探测周期增加到数十万秒时,视电阻率会减小。 相似文献
11.
巷?孔瞬变电磁法在存在明显电导率各向异性的勘探区会产生较大的解释误差。基于时域有限差分算法,通过引入电导率各向异性张量构建控制方程、将矩形回线源电流密度加入Maxwell方程安培环路定理实现任意电流源的加载,以差分代替微分对控制方程进行离散,实现巷?孔瞬变电磁三维模型正演计算。在与解析解对比验证算法计算精度的基础上,构建全空间模型、层状模型和三维块状模型进行正演并分析电导率各向异性对巷?孔瞬变电磁三分量响应的影响程度与方式。结果表明:垂直轴电导率对巷?孔瞬变电磁三分量响应基本没有影响,水平轴电导率对巷?孔瞬变电磁三分量响应影响较大,其中,?Bx/?t响应主要受y轴电导率影响,?By/?t响应主要受x轴电导率影响;通过三分量响应的形态和幅值的相互关系可以辨别各向异性介质所在方位及主轴电导率方向;当异常体所处全空间介质为电导率各向异性时,异常体产生的异常响应会被全空间介质电导率各向异性产生的异常响应所淹没。电导率轴向各向异性特征在巷?孔瞬变电磁法解释过程中不可忽略,研究成果为巷?孔瞬变电磁法各向异性解释提供指导,也为各向异性反演提供参考。 相似文献
12.
大地电磁正演的精度对后续反演意义重大。决定大地电磁有限元正演精度的因素,主要有插值函数、方程组的求解误差,以及网格的剖分程度。而网格剖分的合适与否,为有限元的求解提供了先决条件。因此,从边界条件和横向网格以及纵向网格三方面,讨论了网格剖分对大地电磁正演精度的影响。经研究表明,对于一维和左右对称的二维介质,左右和下边界都无需放置在无穷远处,并且横向网格对精度影响不大。而由于纵向网格与近地表的最小趋肤深度密切相关,所以地表网格不得大于1/4个趋肤深度。 相似文献
13.
频率域电磁剖面有限差分法2.5维正演数值模拟 总被引：1，自引：0，他引：1
对于实际生产中遇到的2.5维问题,从麦克斯韦方程组出发,把地电参数变化小的走向方向转化到波数域,用一系列波数模拟三维源的特征,并在波数域中,得到2组关于Hx (kx,y,z)和Ex (kx,y,z)的偏微分方程。选取适当的kx值,用有限差分法在y-z平面的网格中求解,再通过反傅里叶变换得到空间域中的电磁场。在验证了算法的正确性之后,对不同埋深的直立异常体、倾斜异常体及断陷模型进行了数值模拟,其结果直观地显示了异常体磁异常的形态,同时研究了程序对低阻异常体的横向分辨率。对层状大地模型,用井间收发方式进行了模拟研究,并取得了较好的结果,这对今后实际勘探应用奠定了基础。 相似文献
14.
潜艇在现代海军中占有举足轻重的地位,所以反潜技术也越来越重要,相对于传统的声波、磁法等反潜技术,电磁法探测是近几年才开始起步的。这里研究了电磁法探测水下二维体的问题,基于天然电磁场的有限单元正演,同时考虑电阻率与磁导率,采用三角剖分,通过模型试算验证了天然电磁场的有限元算法在水下二维体探测中,具有较高的精度和适用性,模拟结果 TM模式优于TE模式,可以为电磁法探测水下潜艇打下坚实基础。 相似文献
15.
High-performance algorithms and software have been designed for fast 2D forward modeling and inversion of high-frequency electromagnetic logs to recover the patterns of electrical conductivity and relative permittivity (dielectric constant) of rocks. Forward modeling implies linear representation, with regard to eddy current, of attenuations and phase shifts of electromagnetic fields propagating in 2D conducting formations. In order to improve the modeling quality, the background model is chosen proceeding from the eddy current pattern in a heterogeneous subsurface. The inverse solution is obtained by pseudoinversion of the sensitivity matrix of relative amplitudes and phases. Joint 2D inversion of phase shifts and attenuations improves the reliability of the inferred spatial conductivity distribution and allows estimating permittivity. The reported numerical experiments on inversion of synthetic and field EM logs prove the possibility of estimating the permittivity of complex terrigenous and carbonate reservoirs with water- and oil-filled porosity. 相似文献
16.
赵晓博 《物探化探计算技术》2011,33(5):517-521,463
利用Delaunay三角化这种网格非结构化方法。通过编程实现了二维模型的非结构化三角形网格剖分,并编写了中心回线法瞬变电磁2．5维有限元正演程序。与前人计算结果对比,在取得相同计算精度的情况下,与结构化网格相比,非结构化网格所需网格和节点数量大大减少,计算效率更高。通过将非结构化网格法引入到瞬变电磁2．5维正演模拟中,实现了对复杂二维地电模型的有限元数值模拟,提高了现有有限元算法的应用范围。 相似文献
17.
18.
为了分析阶梯钻头冲击挤密作用下土体的变形机理,将土体视为DP材料,服从Drucker-Prager屈服准则,通过钻具与土体之间的接触以及施加冲击载荷来模拟钻进过程,建立了潜孔锤冲击挤密钻进的有限元模型。以非线性瞬态动力学有限元分析为手段,利用有限元软件ANSYS9.0对阶梯钻头冲击挤密作用下土体的变形机理进行了模拟分析。结果表明,在冲击载荷作用下,冲击功主要用于竖向压缩土体,水平方向的挤压作用相对较小。阶梯钻头对土体的破坏作用可以看作是每个台阶的微剪切作用,当台阶布置合理时,土体沿等效锥面整体破坏,破坏后的土体在等效锥面的作用下被挤向孔壁四周,挤密效果较好。 相似文献