共查询到16条相似文献，搜索用时 171 毫秒
1.
电阻率法有限单元正演模拟中,采用第三类边界条件时为保证精度仍要求较大范围的计算域.无单元法为地球物理领域的新兴正演模拟方法,其计算效率低,但其中采用的移动最小二乘(MLS)形函数相比于有限单元法形函数具有良好的连续性,模拟精度高.本文将MLS形函数应用于电阻率法有限单元2.5维正演的第三类边界条件处理,提出电阻率法有限单元-移动最小二乘(FEM-MLS)耦合正演方法.通过不同正演方法的模型算例模拟结果对比,验证了本文算法的有效性,并讨论了各个参数选择对模拟结果的影响.本文数值模拟结果表明采用第三类边界条件时,在同等计算精度前提下,FEM-MLS耦合法相比于有限单元法可进一步缩小计算域并提高了计算效率,相比于采用较大计算域满足边界条件的有限单元法计算效率提高了约一倍,相比于采用相同小范围计算域的有限单元法平均精度提高了约一倍. 相似文献
2.
全局弱式无单元法是在有限单元法基础上发展起来的一种数值模拟方法,它采用局部支持域内的节点信息来构造形函数实现局部精确逼近,摆脱了单元,仅依赖于节点信息,具有预处理简单、模拟精度高、灵活性强的特点,适用于复杂地电条件下直流电阻率正演模拟.本文采用RPIM构造直流电阻率全局弱式无单元法形函数,利用RPIM形函数推导了直流电阻率全局弱式无单元法方程.然后,编制了直流电阻率全局弱式无单元法正演模拟Fortran程序,利用该程序对典型的地电模型进行了正演模拟,并将正演结果与基于线性插值的FEM正演结果及解析解进行对比,结果表明采用RPIM形函数的全局弱式无单元法用于直流电阻率正演模拟的正确性及有效性,且在同等条件下,全局弱式无单元法模拟精度高于矩形剖分的FEM,更有利于指导电法勘探的数据解译;利用该程序对复杂地电模型进行了正演模拟,结果表明全局弱式无单元法对复杂地电模型模拟效果良好,适应性强,灵活性高,可任意加密节点提高模拟精度. 相似文献
3.
无单元Galerkin法采用滑动最小二乘法拟合场函数,只需节点无需单元,具有前处理简单、精度高、解高次连续等优点,被用于求解探地雷达(GPR)正问题.本文从Maxwell方程出发,推导了GPR正演需满足的波动方程;详细介绍了滑动最小二乘法形函数的构造方法.针对EFGM不满足插值条件导致强加边界条件的处理变复杂的特性,采用罚因子法对强加边界条件进行了处理;同时为了消除EFGM进行GPR正演模拟时来自截断边界处的超强反射,采用透射边界条件把GPR波在截断边界处的反射波透射出去,进而压制了来自截断边界处的反射波.然后,编制了EFGM的GPR正演模拟Matlab程序,应用该程序对典型GPR地电模型进行了正演模拟,并把该正演剖面图与基于线性插值FEM正演剖面图进行了对比,结果表明了EFGM用于GPR正演计算的正确性及有效性,并且在相同节点数条件下,EFGM比矩形剖分的FEM的精度要高,更有利于指导雷达剖面的数据解译. 相似文献
4.
直流电阻率数值模拟对于实测电阻率数据反演、解释地下真实电阻率结构具有重要意义.为解决复杂地电模型剖分以及处理截断边界时面临的困难,本文提出一种新的三维自然单元-无限元耦合法.该算法在自然单元区域灵活布设自然节点以填充复杂模型,并引入无限单元替代截断边界条件以促使刚度矩阵与电源位置无关.首先推导了稳定电流场基本方程;然后介绍了二、三维自然单元法及其形函数的构建过程;再引入三维多向映射无限单元并实现自然单元-无限元的有效耦合;最后通过4个数值模型,结合传统有限单元法、有限元-无限元耦合法以及COMSOL Multiphysics软件的计算结果,验证了算法的正确性及有效性.数值结果表明,本文所提算法满足精度要求,能有效减小计算区域,无需二次更新刚度矩阵,对复杂模型适应性强. 相似文献
5.
系统地论述了用有限单元法研究复杂地形条件下三维直流电阻率的正演计算技术.首先给出了三维构造中点源电场的边值问题以及相应的变分问题;然后利用有限单元法求解变分问题,采用四面体单元对研究区域进行剖分,在单元中进行三线性函数插值,将变分方程化为线性代数方程组;最后,考虑到节约计算时间,利用对称超松弛顸条件共轭梯度迭代算法求解大型线性方程组,得到了各节点的电位值,进而计算出地表的视电阻率.通过理论模型的计算检验了算法的可行性之后,给出了几种常见纯地形异常的数值模拟结果和一个组合模型的计算结果,其研究工作为研究三维直流电阻率反演奠定了基础. 相似文献
6.
频率域可控源电磁是在大地电磁测深的基础上发展起来的一种人工源电磁测深法,其二维电磁响应的计算须采用数值模拟方法.本文以Matlab为程序编译工具,采用双二次插值的有限单元法,推导出相应的计算公式.为了模拟无穷远边界及满足计算机的内存需求,在保证计算精度的情况下设计了非均匀网格剖分.在程序编制中,只存储有限元系数矩阵的非零元素,大大减少了正演计算的时间.针对频率域可控源电磁法中卡尼亚电阻率在过渡区和近区畸变的问题,给出了全区视电阻率的迭代公式,并对典型的一维层状模型以及简单二维模型进行了计算.过渡区和近区数据经过校正后,可以正确反映出模型的地电特征,证明了线源下近区勘探的可能性. 相似文献
7.
8.
对常见的隐伏地质构造体进行大地电磁场正演模拟,并研究和总结了构造体对大地电磁场的响应特征.由麦克斯韦方程组导出大地电磁场的边界条件,应用有限单元法求解变分问题并得出刚度矩阵,根据线性方程组计算出相应的场值,再由改进的视电阻率定义式求出视电阻率.模拟结果显示:改进的视电阻率定义式相对于传统的Cagniard定义式能有效地提高视电阻率曲线的收敛速度和逼近程度;隐伏地质构造体对大地电磁场的响应具有一定的规律性,通过分析两种极化模式下的正演结果,可以有效地识别地质构造体的存在和分布情况.通过实例分析,得出隐伏地质构造正演模拟与实际观测资料的互相验证关系,为实践应用中的反演工作提供了重要的指导意义. 相似文献
9.
伴随着复电阻率法的广泛应用,发展精确和快速的正演和反演算法成为复电阻率法研究的重点.本文采用基于三角单元剖分的有限单元法进行了复电阻率二维数值模拟研究.为了提高计算速度,对无穷远边界进行了近似处理.整个正演计算过程分为两步,首先采用有限单元法计算四个不同频率的视复电阻率数据,然后对前一步得到的视复电阻率数据采用递推算法计算视Cole-Cole参数.采用这种正演算法与一维正演的结果进行了对比,验证了本文方法的正确性.设计了两个二维极化模型,数值模拟结果表明视复电阻率和Cole-Cole视参数等值线断面图对于异常目标体都有比较明显的反映. 相似文献
10.
本文前一部分叙述了点源二维电阻率法数值模拟正演计算的有限单元法,文中采用混合边界条件、用LL~T分解解线性方程组等优化措施,使二维有限单元法电阻率法正演计算的速度和精度都比目前国外流行的L.Rijo方法和程序有了进一步的提高。文中介绍了有关的方法并引入了若干计算结果。本文后一部分叙述了三维电阻率法数值模拟的积分方程近似解法,由于对K.Dieter等人提出的方法作了一些近似处理,并用迭代法求解积分方程表示式,实算结果说明在保证计算精度的条件下,提高了计算速度,从而提高了三维电场正演计算的有效性和实用性。所提出的方法容易推广到激发极化法的正演计算中。 相似文献
11.
Precision of meshfree methods and application to forward modeling of two-dimensional electromagnetic sources 总被引：1，自引：0，他引：1
Meshfree method offers high accuracy and computational capability and constructs the shape function without relying on predefined elements. We comparatively analyze the global weak form meshfree methods, such as element-free Galerkin method (EFGM), the point interpolation method (PIM), and the radial point interpolation method (RPIM). Taking two dimensional Poisson equation as an example, we discuss the support-domain dimensionless size, the field nodes, and background element settings with respect to their effect on calculation accuracy of the meshfree method. RPIM and EFGM are applied to controlled-source two-dimensional electromagnetic modeling with fixed shape parameters. The accuracy of boundary conditions imposed directly and by a penalty function are discussed in the case of forward modeling of two-dimensional magnetotellurics in a homogeneous medium model. The coupling algorithm of EFG–PIM and EFG–RPIM are generated by integrating the PIM or RPIM and EFGM. The results of the numerical modeling suggest the following. First, the proposed meshfree method and corresponding coupled methods are well-suited for electromagnetic numerical modeling. The accuracy of the algorithm is the highest when the support-domain dimensionless size is 1.0 and the distribution of field nodes is consistent with the nodes of background elements. Second, the accuracy of PIM and RPIM are lower than that of EFGM for the Poisson equation but higher than EFGM for the homogeneous medium MT response. Third, RPIM overcomes the matrix inversion problem of PIM and has a wider selection of support-domain dimensionless sizes as compared to RPIM. 相似文献
12.
本文探讨了将商业有限元软件COMSOL Multiphysics应用于地震电信号选择性三维数值模拟方面可能面临的一些问题.在检验算法可靠性的同时,通过选取不同类型的边界条件、不同大小的计算区域、不同的介质电性参数以及地表介质电性非均匀性和地下特定导电通道等不同电性结构模型等,考察了影响地震电信号选择性数值模拟结果的一些可能因素.结果表明,在开展地震电信号选择性数值模拟研究时需对模型参数的选取进行必要的分析,而地表介质电性非均匀性以及地下特定导电通道等可以较好地解释地震电信号选择性现象,其贡献的大小与模型具体的电性结构和对比度等参数有关. 相似文献
13.
Our aim is to introduce the Coupled Finite-Infinite Element Method (CFIEM) as a new alternative approach to the Earth’s gravity
field modelling. We show that if the computational domain is large enough in radial direction, one can obtain the qualitatively
and quantitatively comparable solution to the solution by the Finite Element Method (FEM). We study the influence of the size
of the computational domain on the final CFIEM solution as well as the successive refinement of the discretization and its
convergence to the exact solution. As an input data we use the synthetic boundary conditions computed from a Synthetic Earth
Gravity Model (SEGM) and we test the CFIEM solution by the data generated directly from SEGM and the solution by the FEM. 相似文献
14.
15.
在三维频率域电磁法的正演模拟方法中,有限元方法具有计算精度高、适应性强的优点,近年来来得到了越来越多的关注.在正演过程中,主要的计算量集中在求解由偏微分方程组离散得到的线性方程组上,因此求解线性方程组关系着正演计算速度以及模拟精度.由于由有限元方法离散得到的复系数线性方程组条件数非常大,使用常规的迭代法和预条件很难收敛.目前大多数的研究工作采用直接解法,需要大量的计算机内存,限制了可求解问题的规模.本文研究了线性方程组的迭代解法,通过将复系数线性方程组转化为其实对称形式,构造分块对角预条件.在应用预条件的过程中,需要求解两个较小的实数方程,通过辅助空间解法求解.本文的算法适用于可控源电磁法和大地电磁法,对一系列的数值算例的模拟结果证明了迭代算法的效率,结果表明迭代算法可以在小于20次迭代内收敛,同时迭代次数与模型电阻率、问题规模和频率无关. 相似文献
16.
本文从总电位法和异常电位法两个方面,分别就直流电法有限元中所采用的三类边界条件对计算结果的影响进行了分析,为采用合适的边界条件及选取最佳边界距离来提高正演精度提供了基础.结果表明,混合边界条件精度最高,可大大缩小求解区域而不影响计算精度,其次是Dirichlet边界条件精度较高,但测点越靠边界误差会越大,必须取足够大的边界区域,齐次边界条件的误差比较大,但如果采用非二极装置,通过电位差计算得到的视电阻率,由于无穷远边界对电位差的影响基本消除,视电阻率计算误差与混合边界条件下的接近,在反演中,为了节省计算时间,经常使用齐次边界条件进行有限元正演. 相似文献