共查询到20条相似文献，搜索用时 47 毫秒
1.
为提高重力及重力梯度数据的正演效率,笔者引入基于MPI (Message Passing Interface)和OpenMP (Open Multi-Processing)的并行计算,通过对比分析不同数据规模的网格数和模型体个数对并行效率和加速比的影响,得出随着正演数据规模的增加,并行效率和加速比均得到提高。同时对比了基于MPI和基于OpenMP的两种并行方式的性能,结果表明,重力及重力梯度数据正演的并行计算中MPI的提速能力优于OpenMP,且在较大规模数据的正演计算中基于MPI的并行效率优于基于OpenMP的并行效率。 相似文献
2.
目前有限元技术的开发及在电磁勘查技术中的应用已日趋成熟.然而,有限元正演模拟主要集中于地面和海洋电磁,航空电磁三维有限元模拟尚没有受到足够重视.以前人工作为基础,利用结构化网格实现了航空电磁系统的三维矢量有限元正演模拟.从二次场双旋度矢量非齐次亥姆霍兹方程出发,应用广义变分原理推导出变分方程,并采用六面体单元剖分,将场置于单元棱边上,对每个单元应用线性插值,最后合成含有稀疏矩阵的线性方程组.针对航空电磁多源性问题,利用MUMPS(multifrontal massively parallel sparse direct solver)直接求解器进行求解,在保证精度的前提下大幅度提高计算速度.利用单个异常体模型检验算法的精度和软件的稳定性,进而通过典型地电模型的模拟验证算法的有效性.对不同地下电性结构正演模拟结果进行对比分析,进一步研究了覆盖层和垂直接触带等典型构造对航空电磁响应的影响特征. 相似文献
3.
《物探化探计算技术》2017,(2)
由于航空瞬变电磁的数据量很大,外加瞬变电磁反演方法耗时较大地电磁更多,为了在后期资料处理时减少不必要的耗时,对正演计算中的频时转换的方法做了一定研究,采用基于拉普拉斯尺度变换性质的Gaver-Stehfest变换做频时转换,适当减少了编写程序代码的行数及程序计算的时间消耗。针对现在的计算机很多都是多核的特点,研究了航空瞬变电磁一维反演的MPI(Message Passing Interface)并行计算,大大地节约了运算时间,分析了并行计算时不同的多核CPU和不同的进程数时的加速比及并行效率。 相似文献
4.
《物探化探计算技术》2017,(1)
半航空时间域电磁法集成了地面和航空电磁法各自的优势,不仅容易实现大发射磁矩,增大勘探深度,同时还可以提高勘探效率。但是由于半航空时间域电磁法在测量过程中,收、发距会发生变化,因此与地面和航空电磁法相比,其数据处理与反演解释更加困难。以长导线电性源半航空电磁法为例,将自适应正则化反演算法应用于半航空时间域电磁数据反演中。首先给出最平缓模型约束条件下的半航空时间域电磁数据自适应正则化反演算法,并引入并行技术实现了该反演算法的并行计算;然后对两种三层地电模型和一种六层地电模型的理论电磁响应加高斯白噪声后的数据进行反演计算;最后对其反演结果和并行计算的效率进行了详细地分析与讨论。实验结果分析表明:自适应正则化反演算法是一种有效的半航空时间域电磁数据反演解释方法,该方法具有很好的收敛性和稳定性;自适应正则化反演算法并行计算能大大减少反演计算时间,提高反演效率。 相似文献
5.
时间域航空电磁法是一种高效的地球物理勘查方法。由于激电效应的存在,晚期道电磁数据经常会出现符号反转现象。此时,传统的不考虑激电效应的成像和反演方法无法恢复地下真实电性结构,因此需要开发含激电效应的电磁正反演模拟技术。本文提出了一种基于非结构有限元的时间域航空电磁激电效应三维正演算法。为描述介质的激电效应,引入Debye模型,利用频时转换将其变成时间域导数形式,进而采用无条件稳定的二阶后推欧拉方法进行时间离散。之后,采用基于非结构四面体网格的矢量有限元法,对麦克斯韦方程进行空间离散,实现含激电效应的时间域航空电磁正演问题求解。通过与均匀半空间模型的一维半解析解进行对比,验证了该算法的准确性。最后,计算了半空间和三维地质体具有不同极化条件时的电磁响应,发现可以根据航空电磁响应曲线的衰减特征,有效识别围岩及异常体的极化特征。此外,激电参数对时间域航空电磁系统的影响特征研究表明,地质体电阻率越高或者充电率越强,极化效应影响越严重。 相似文献
6.
《物探与化探》2015,(5)
采用三角网格剖分的有限元法,研究了2.5维航空瞬变电磁法正演模拟问题。利用时频变换数值方法将时间域电磁场转换到拉氏域,再利用傅里叶变换将三维问题降维变为2.5维问题,然后由有限元法求解得到拉氏域二维电磁场,逆拉氏变换后得到时间域航空瞬变响应。为了回避正演模拟中总感应磁场在场源处的奇异性问题,采用异常场算法,场源响应通过在微分方程中施加背景电磁场实现。由于瞬变电磁信号具有较大的动态范围,而且需要经过两次正、逆拉氏变换和傅里叶变换,每个环节的计算精度和速度要严格控制在较高的水平上,否则积累误差会非常大。模型计算表明均匀大地和层状大地模型解析解与数值解吻合很好。这证明该算法是正确可行的,可作为研究二维复杂地质体的方法手段。 相似文献
7.
复电阻率法二维正演采用二次场算法,结合有限单元法进行离散化,然后引入ColeCole模型,完成数值模拟。通过理论模型试算验证了二维正演程序的正确性。当模型剖分较细时,多发射源、多频率的正演计算效率有待提高。复电阻率法二维数值模拟是按发射源或频率计算,且各发射源之间、各频率之间计算是相互独立且互不影响的,因此通过引入MPI并行计算来提高其计算效率,从而减少其正演所需的时间。采用串行、并行复电阻率法二维正演程序来计算同一理论模型,结果验证了该并行程序的正确性。对并行算法的效率进行分析表明,该并行算法是可靠的、高效的。 相似文献
8.
赵晓博 《物探化探计算技术》2011,33(5):517-521,463
利用Delaunay三角化这种网格非结构化方法。通过编程实现了二维模型的非结构化三角形网格剖分,并编写了中心回线法瞬变电磁2．5维有限元正演程序。与前人计算结果对比,在取得相同计算精度的情况下,与结构化网格相比,非结构化网格所需网格和节点数量大大减少,计算效率更高。通过将非结构化网格法引入到瞬变电磁2．5维正演模拟中,实现了对复杂二维地电模型的有限元数值模拟,提高了现有有限元算法的应用范围。 相似文献
9.
《物探化探计算技术》2017,(1)
地震波场频率域正演是频率域全波形反演的基础。针对反演计算量巨大的问题,利用频率域二维波形正演算法中频率和炮点计算的独立性,开发出粗细粒度结合的MPI并行算法:粗粒度为频率并行,细粒度为单个频率解方程并行。实现方法是将正演频率分组放入不同的MPI通讯域内,每个通讯域内单个频率求解方程过程,采用基于MPI的MUMPS(多波前大规模并行稀疏直接解法器)软件包并行加速。模型测试结果表明:MPI多通讯域并行算法计算结果正确,计算效率显著提高,加速效果稳定。 相似文献
10.
研究了二阶声波方程频域PML边界条件和频域变网格并行计算技术。PML边界是一种较为理想的吸收边界方法,多用在求解时域应力速度方程中,但对于频域声波正演,二阶位移方程更常用。从一阶声波方程PML吸收边界条件导出频域二阶位移方程PML边界条件,模拟算例得到的频率切片、时间切片和地震记录对比都说明该边界条件吸收效果很好。频域单炮正演不同频率间是独立的,据此低频部分采用大网格计算,高频采用小网格,实现变网格步长计算技术,这是较时间域正演的一个优势,在保证模拟质量的同时,减少计算量和内存消耗。 相似文献
11.
为了研究复杂地电模型的航空瞬变电磁法全波形响应特征,需要开发考虑发射波形的三维数值模拟算法。本研究基于非结构四面体网格和位移逆Krylov子空间(Shift-and-Invert Krylov,简称SAI Krylov)方法,采用基于电偶极子离散的场源处理方法模拟场源,在时间域进行计算实现了全波形航空瞬变电磁法矢量有限元三维数值模拟。使用均匀半空间模型在阶跃波、半正弦波、三角波和梯形波激发下的全波形解析解、VTEM实际激发波形的后推欧拉算法计算结果,检验了本研究开发的数值模拟算法的正确性。设计地表起伏异常体模型,计算和分析了航空瞬变电磁响应特征。开发的基于位移逆Krylov子空间的全波形航空瞬变电磁法三维数值模拟算法适合模拟复杂地电模型的响应,具有较高的计算精度。 相似文献
12.
《物探化探计算技术》2019,(5)
有限元求解大地电磁正演问题时,对研究区域的剖分常规做法是采用规则化的网格。但规则化网格在剖分地形、断层、褶皱等复杂模型会产生较大的几何离散误差。针对上述情况,这里采用非结构化的四边形网格对二维地电模型进行网格剖分,并与自适应有限元相结合,由剖分的粗网格出发,利用每个频点下网格单元的后验误差估计值指导网格的局部加密,优化网格质量和数量,从而提高正演模拟的精度。通过一维K型地电模型利用本文方法算出数值解与解析解进行对比分析,验证了自适应有限元法在求解大地电磁二维正演模拟中的有效性。并通过对断层、褶皱模型的正演模拟,分析了其大地电磁正演响应特征。 相似文献
13.
利用有限元方法进行大地电磁正演数值模拟时,由于是在有限网格区域上的数值计算,模拟计算时的网格边界为截断边界,而有限元数值模拟时的大地电磁场边界条件需要在足够远处才能够满足,所以截断边界的存在可能会使大地电磁正演模拟的边界条件无法满足,致使对计算结果和计算精度产生影响。利用有限元二维正演程序,在网格边界处加载一维情况下的大地电磁场,然后固定研究区域的网格剖分,并对一维地电模型和二维地电模型在改变有限元网格边界大小的情况下进行计算。在对一维模型进行模拟计算时,截断边界对边界条件没有影响,边界条件自然满足。而对二维模型进行模拟计算时,截断边界的存在对计算结果有较大影响。利用趋肤深度作为有限元网格边界变化的量度,通过改变网格边界大小,对不同的二维地电模型进行计算比较,总结出适合大地电磁有限元正演模拟的参考网格边界。 相似文献
14.
积分方程法三维模拟井地电法并行算法研究 总被引：1，自引：1，他引：1
在积分方程法模拟三维电磁响应原理的基础上,对井地电法三维正演模拟并行算法做了研究。同时,还利用Fortran 6.5和消息传递接口(MPI)开发了井地电法三维正演模拟并行计算程序。在求解大型线性方程组时,采用的是稳定的双共轭梯度(BICGSTAB)算法,单机正演程序与并行正演程序的计算结果有着很好的一致性,这说明了并行算法的正确性,为进一步开发井地电法三维反演并行程序奠定了基础,也为井地电法三维正、反演程序实用化提供了技术支持。 相似文献
15.
近年来,瑞利波波形反演技术因其避开了常规频散曲线计算,直接进行波场计算和反演不再受水平层状介质理论假设的限制,得到广大学者的高度重视。但瑞利波波形反演过程中需要不断进行波场正演和逆推计算。另外,由于浅地表速度较小,模拟计算时需要较小的网格间距才能避免数值频散,这无疑大大增加了正演模拟的计算量。对于这一问题,通常采用并行化设计来提高正演模拟的计算效率。本文基于消息传递接口（MPI）并行有限差分算法,以区域分解思路将模型区间分解成若干子区域,各区域互相通信,共同完成对模型的正演计算。并详细给出了区域分解、坐标转换、区域通信、波场合并等并行方案中的具体实现方法和实现步骤。通过对弹性模型、Kelvin黏弹性模型和标准线弹性固体（SLS）黏弹性模型不同并行方案的计算结果进行分析,验证了本文并行方案的可行性和有效性。并行计算结果表明,与单处理器计算时间相比,增加处理器数目可以明显减少计算时间,但随着处理器数目的增加,不同处理器之间的通信时间也增大;因此,并行时需要选择合适的处理器数目。对于黏弹性介质模型,SLS黏弹性模型的并行计算效率优于Kelvin黏弹性模型。 相似文献
16.
为了加快大地电磁三维正演的求解速度,本文将一种新型的代数多重网格算法——聚集多重网格（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程序能够提高十几倍的计算速度,尤其适合大规模大地电磁三维正演问题。 相似文献
17.
基于AeroTEM直升机时间域航空电磁响应曲线分析 总被引：1，自引：0，他引：1
在航空电磁实际勘查中,有较多因素影响系统的勘探精度和分辨率.这里以AeroTEM系统为例,基于MAXWELL软件正演模拟算法,对异常体横向距离变化、异常体倾角变化及异常体埋深变化对电磁响应曲线的影响进行分析.正演模拟结果表明:时间域直升机航空电磁系统具有较高的横向分辨率,对不同倾角的异常体都有较好的响应以及具有较大的勘探深度. 相似文献
18.
讨论了亿万单元有限元模型的可扩展并行计算方法。从软件和硬件两个方面提出了前处理、并行计算方法、程序算法、后处理实现等核心问题的解决方案。采用网格加密方法生成一亿单元的有限元模型,利用对偶原始有限元撕裂内联法（FETI-DP）求解系统方程。基于图论理论建立了子区域间的通讯拓扑关系,实现了子区域间点对点通讯,避免速度慢、通信量大的全局通讯。在自主开发程序基础上,增加相应模块,采用面向对象编程技术和MPI消息传递库开发程序。对一个一亿多单元的工程实例运用5 000核并行计算,得到了超线性加速比。计算结果在专用图形工作站上进行后处理,显示和交互操作速度良好。研究在两方面实现了突破：一是将模型规模提高到了一亿多单元;二是同时调用了5 000个计算核来并行运算,并得到了很高效率。高分辨率有限元并行模拟研究成果可为岩土工程中结构特别复杂、计算区域特别大、地质情况复杂等模拟提供很好的技术方法和实现手段。 相似文献
19.
为了更好地模拟地下介质连续变化及开展连续介质的反演,对二维电导率分块线性变化的线源频率域可控源电磁法进行了有限元正演模拟,在剖分单元内同时对电场及电导率参数线性插值,使电导率参数在剖分单元之间保持连续变化。首先,提出有限元正演模拟的边值问题及变分问题,并详细论述了有限元的剖分、插值、单元分析及总体合成的各个步骤;其次,采用稀疏存储及基于不完全LU分解的BICGSTAB算法求解复系数方程组,节省了内存并提高了计算速度;然后,对一个均匀半空间模型进行模拟,计算结果表明,低频及高频的有限元数值解都与解析解吻合,证明了算法的正确性;最后,对水平层状模型及垂直断层模型进行正演计算,视电阻率及相位的等值线图均较好地反映出了异常体,说明文中算法能够对电导率连续变化的线源可控源电磁法进行有效地模拟。 相似文献