共查询到20条相似文献,搜索用时 31 毫秒
1.
《物探化探计算技术》2016,(3)
用fortran语言编程实现了有限元三维大地电磁正演,通过层状介质模型、二维棱柱体模型及三维低阻体模型结果的对比,验证了所编写程序的正确性。首先通过加入第一类边界条件,减少了最终求解方程组的维数,同时对系数矩阵的存储采用非零存储技术,大大降低了对计算机内存的需求;最后在串行程序的基础上,基于MPI实现了频点间的并行,并对一个三维模型进行计算,并行后开启4进程时加速比达到了2.328,有效地减少了所需时间。 相似文献
2.
直流电阻率测深中二维与三维反演结果的对比与分析 总被引:1,自引:0,他引:1
直流电阻率测深中,用二维反演程序对三维地质体进行了反演,并与三维反演结果进行对比和分析。首先,对比二维和三维最小二乘反演在正演模拟算法和先验信息的确定方法;然后,对若干比较典型的模型进行反演实例对比。由于二维程序反演仅考虑单个剖面的电阻率信息,无论是在异常位置、形态及电阻率特性上,其反演精度都比较低;三维电阻率反演综合了多个测深剖面的电阻率信息,其结果与实际模型吻合得非常好。 相似文献
3.
4.
余涛朱自强鲁光银曹书锦 《物探化探计算技术》2013,(4):446-451
重力张量是重力位的二阶空间导数,对密度体的变化和细节部份反映更为灵敏。对于复杂的重力密度模型来说,张量的解析公式是很难推导的。为了模拟这种情况,将复杂模型进行有限元剖分,计算每个单元的重力异常对测点张量的影响,最后叠加得到整个复杂模型在测点处的重力张量。采用Delaunay四面体非结构化网格对密度体进行剖分,经分析表明,这种剖分方式具有较高的精度。有限元剖分的方式拟合复杂形体缺点之一是计算量比较大,而并行化能有效解决这个问题。这里基于OpenMP并行模型编写了重力张量并行正演程序,分析了不同情况下并行的执行性能,为大规模复杂模型的重力张量正演提供了一种并行策略和思路。 相似文献
5.
重力张量是重力位的二阶空间导数,对密度体的变化和细节部份反映更为灵敏.对于复杂的重力密度模型来说,张量的解析公式是很难推导的.为了模拟这种情况,将复杂模型进行有限元剖分,计算每个单元的重力异常对测点张量的影响,最后叠加得到整个复杂模型在测点处的重力张量.采用Delaunay四面体非结构化网格对密度体进行剖分,经分析表明,这种剖分方式具有较高的精度.有限元剖分的方式拟合复杂形体缺点之一是计算量比较大,而并行化能有效解决这个问题.这里基于OpenMP并行模型编写了重力张量并行正演程序,分析了不同情况下并行的执行性能,为大规模复杂模型的重力张量正演提供了一种并行策略和思路. 相似文献
6.
《岩土力学》2020,(Z1)
混凝土面板与堆石体之间由于变形不协调导致的面板脱空是面板堆石坝工程中常见的现象,严重时可引起面板的开裂并威胁大坝的整体安全,非常有必要计算分析面板是否会出现脱空以及脱空值的大小。针对以往的计算方法难以在直观性、计算精度和容易实施之间取得平衡的现状,尝试将脱空视为物体内部不连续界面发生的相对位移,用扩展有限元法(XFEM)进行模拟分析。首先将一种复合式XFEM接触本构模型扩展到三维,实现在三维条件下XFEM对不连续界面发生张开和错动的表述,再通过一个滑块模型算例用自编程序和ABAQUS的接触力学方法进行模拟和对比,验证该三维界面接触算法的正确性,最后以天生桥面板堆石坝为算例计算了面板部位的脱空情况,并与实测值和对偶mortar方法的计算结果进行对比,验证该方法计算大型实际工程的能力。 相似文献
7.
由于地质体和矿体的形态非常复杂,使用长方体网格离散建立正演模型时可能和真实情况有很大差别,因此计算结果可靠性差。本文提出一种基于约束Delaunay网格剖分的方法对地质体进行离散并进行重力建模,在模型边界等复杂区域使用网格自适应加密技术,将三维地质体离散为有限个四面体;并详细推导出针对四面体网格的重力正演公式,实现了基于约束Delaunay网格剖分技术的三维重力数值模拟;最后,针对一个合成数据模型,将计算解与解析解对比。结果表明,细化网格的模拟结果比粗糙网格更好,满足数值模拟的精度要求。将该方法应用到金川矿区实际地质体建模中,根据局部需要,建立各处网格密度不均匀的三维模型,并计算该模型的地表重力场,而后对比模拟数据与实测数据,结果表明Delaunay网格建模方法具有很强的适用性,能够模拟复杂的地质体重力异常。 相似文献
8.
矿区疏干地下水三维流与二维流对比研究 总被引:1,自引:0,他引:1
分别从解析解、数值解和电网络三个方面,对三维流和二维流模型求解方法和结果进行了理论分析和对比,讨论了在参数拟合、水量水位预测方面的差别及处理对策。最后给出元宝山露天煤矿疏干排水三维流模型与二维流模型的计算实例及分析,并提出了矿区大降深、大流量疏排地下水属明显的三维流过程,必需用三维流模型进行模拟计算。 相似文献
9.
为了加快大地电磁三维正演的求解速度,本文将一种新型的代数多重网格算法——聚集多重网格(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程序能够提高十几倍的计算速度,尤其适合大规模大地电磁三维正演问题。 相似文献
10.
11.
积分方程法三维模拟井地电法并行算法研究 总被引:1,自引:1,他引:1
在积分方程法模拟三维电磁响应原理的基础上,对井地电法三维正演模拟并行算法做了研究。同时,还利用Fortran 6.5和消息传递接口(MPI)开发了井地电法三维正演模拟并行计算程序。在求解大型线性方程组时,采用的是稳定的双共轭梯度(BICGSTAB)算法,单机正演程序与并行正演程序的计算结果有着很好的一致性,这说明了并行算法的正确性,为进一步开发井地电法三维反演并行程序奠定了基础,也为井地电法三维正、反演程序实用化提供了技术支持。 相似文献
12.
13.
14.
反应运移模拟在实际应用中的主要制约因素是其计算量巨大,对模拟程序并行化可以拓展反应运移模拟的应用范围。在阐述反应运移模拟耦合过程的基础上,介绍了TOUGHREACT基于MPI的并行实现方法。以ECO2N模块为例,通过分析程序的流程图、加速比和并行效率,探讨TOUGHREACT并行程序的计算性能,并通过加密网格的方式对其加速效果进行验证。结果表明:基于MPI的TOUGHREACT并行程序采用METIS区域分解方法和AZTEC并行迭代解法器库,有效地实现了对TOUGHREACT的并行加速;在计算集群上调用2~128个处理器时,可以达到2.00~68.09倍的相对加速,显著提高了程序的计算效率;网格的剖分精度越高,该程序的并行效果就越好,表明该程序对复杂问题的并行加速效果更优。 相似文献
15.
体积分方程法模拟复电阻率三维体电磁响应 总被引:1,自引:0,他引:1
利用体积分方程法计算了均匀半空间中复电阻率(激电)效应和电磁效应同时存在时的三维体响应。在计算中,对于需要计算三重积分的张量格林函数电荷项一次部分,应用一种差分近似的方法求解,这种方法在保证计算精度的同时更加便于计算机实现;采用二次剖分的算法解决了计算张量格林函数时的奇异值问题;计算含有贝塞尔函数的积分项时,利用一种结合连分式展开的高斯求积代替常规的快速汉克尔变换方法。验证了计算结果并分析了三维体复电阻率模型(Cole-Cole)参数对正演结果的影响,为三维体复电阻率及其参数反演提供了正演依据。 相似文献
16.
在成熟应用的缝洞型碳酸盐岩油藏模拟器KarstSim的基础上,针对更加精细和大规模网格的模拟需求,开发了KarstSim的并行模拟器。并行化主要通过网格分区、大型线性方程解法和处理器之间信息通讯3个部分实现,同时借助了Aztec库和METIS软件。基于Buckley-Leverett两相纵向流动理论和塔河油田某区块的模型,对并行版与解析解、单机版计算结果的对比,验证了并行模拟结果的准确性以及计算的效率性。将并行模拟器应用至64万网格非均质的缝洞型油藏,在最高采用128核服务器中模型计算时间最低为7.3h,计算效率提升显著。最后对不同网格规模的均质油藏进行了模拟,网格规模从1万逐渐增加到1千万规模,结果显示对于千万级网格的均质油藏,并行模拟器计算时间为85min,证实了并行模拟器在大规模油藏模拟中的实用性。 相似文献
17.
18.
19.
20.
变步长激发极化法有限元数值模拟 总被引:1,自引:0,他引:1
用改变剖分步长的方法对三维地电模型进行激发极化法正演模拟计算.首先,在给出三维地电模型异常电位的边值问题与变分问题的基础上,通过对三维激电模型进行变步长的六面体剖分,将异常体部分用较短步长的单元进行加密剖分,并在各个单元上进行插值,于是得到其总刚度矩阵;然后用变带宽方式对其进行编号存储,并选用成熟的LDLT方法求解大型稀疏矩阵,运用等效视电阻率方法求取地表视极化率值.通过用中间梯度装置对几例典型的激电模型进行了正演计算,结果表明,该方法能够达到满意的精度和效果,从而为激发极化正演解释提供有益的帮助. 相似文献