简介:标量CSAMT只适合一维及测量方向与构造方向垂直的二维情况,对于复杂的三维地电结构,CSAMT需采用张量测量。本文试图采用矢量有限元法实现三维张量CSAMT的正演模拟。为了验证算法的正确性,本文在层状介质中计算了三维CSAMT远区的电场,磁场及阻抗张量,并且与层状介质中的理论解进行了比较,接着还模拟了均匀半空间中含有三维异常体的模型,并且分析了四个阻抗张量、视电阻率及阻抗相位的响应特征。得出如下结论:采用矢量有限元法来模拟三维张量CSAMT,其电磁场及阻抗张量的实虚部计算精度都比较高,并且该方法本身满足电场法向不连续,不用进行散度校正。
简介:传统上,有限差分的差分系数一般可以通过泰勒级数展开法或优化方法来极小化频散误差得到。基于泰勒级数展开的差分法在有限的波数范围内精度较高,但在这个范围之外会产生较强的数值频散;基于最小二乘的优化有限差分法能在更大的波数范围内达到较高的精度,并可以在较小的计算需求内获得全局最优解。本文将基于最小二乘的优化有限差分法从二维正演模拟推广到三维,形成了计算效率高、高精度范围宽、适合并行计算的三维声波优化有限差分方法。频散分析及正演模拟表明本文发展的有限差分方法可以很好地压制数值频散。最后,将本文发展的有限差分方法应用到三维逆时偏移的震源波场延拓和检波点波场延拓中,并结合有效边界存储策略与checkpointing技术在GPU集群上实现三维逆时偏移以提高计算效率、减少存储量。三维逆时偏移试算结果表明本文三维优化有限差分方法与传统的有限差分法相比可以获得更高精度的偏移成像结果。
简介:标量CSAMT只适合一维及测量方向与构造方向垂直的二维情况,对于复杂的三维地电结构,CSAMT需采用张量测量。本文试图采用矢量有限元法实现三维张量CSAMT的正演模拟。为了验证算法的正确性,本文在层状介质中计算了三维CSAMT远区的电场,磁场及阻抗张量,并且与层状介质中的理论解进行了比较,接着还模拟了均匀半空间中含有三维异常体的模型,并且分析了四个阻抗张量、视电阻率及阻抗相位的响应特征。得出如下结论:采用矢量有限元法来模拟三维张量CSAMT,其电磁场及阻抗张量的实虚部计算精度都比较高,并且该方法本身满足电场法向不连续,不用进行散度校正。
简介:无网格法形函数构造不依赖预定义的单元,具有计算精度高、处理复杂模型便利等优点。本文介绍了无单元Galerkin法(EFGM)、点插值法(PIM)与径向基点插值法(RPIM)三种全域弱式无网格法的近似原理及特点;以二维泊松方程为例研究了支持域无量纲尺寸、场节点与背景网格设置对无网格法计算精度的影响。将RPIM与EFGM应用于频率域线源二维正演,给出了RPIM形状参数的推荐值;分析了均匀介质模型大地电磁(MT)二维正演无网格法边界条件直接加载与罚函数法加载的精度差异,结合PIM与RPIM边界条件加载便利及EFGM计算复杂模型精度高的优势,提出了EFG—PIM及EFG.RPIM耦合算法,数值计算结果验证了耦合算法的有效性。研究发现:无网格法及其耦合方法适用于电磁法数值模拟;支持域无量纲尺寸取1.0时无网格法精度与效率高,场节点与背景网格重合时计算效果佳;泊松方程求解PIM及RPIM精度较EFGM低,计算均匀介质MT响应精度较EFGM高;RPIM改善了PIM计算涉及的奇异性问题,对应支持域无量纲尺寸选择空间大。
简介:地震波场数值模拟中不可避免地会出现边界反射,一般采用吸收边界条件以压制人工边界反射。目前常用的分裂式完全匹配层(PML)边界条件需要在边界处进行特殊处理,尤其是在三维情况下需要将变量分裂为三个分量,增加了数值模拟的计算时间和内存占用量。与分裂式PML吸收边界条件相比,混合吸收边界条件(HABC)具有易于实现、计算量小和吸收效果好等优点,可以提高三维波动方程数值模拟的计算效率。本文将基于一阶Higdon单程波方程的混合吸收边界条件从二维计算域发展到三维,提出了适用于三维弹性波数值模拟的混合吸收边界条件。均匀模型以及复杂模型的三维数值模拟结果表明,混合吸收边界条件与传统的完全匹配层边界条件相比,具有效率高、吸收效果好的优势。
简介:随着系统仿真技术的快速发展,复杂电磁环境(CEME)已然引起人们的广泛关注,CEME仿真对于雷达预警、体系对抗、移动通信、GPS定位等都有重要作用。首先对于二维抛物方程(2DPE)的推导过程及离散过程进行介绍,并对2D电磁传播的衰减过程进行了仿真。其次,基于2DPE法,研究了并行度高的准三维抛物方程(quasi-3DPE)仿真方法。最后,通过仿真案例,验证了准三维法仿真的三维效果和并行效率,总结比较了3DPE和quasi-3DPE法的计算效率和精度,并借助商业电磁仿真软件WirelessInsie验证了quasi-3DPE法的传播损耗仿真精度。
简介:摘要 目的 比较CT二维法与三维法测量股骨颈前倾角与解剖法的一致性,探讨临床评估前倾角的方法。方法 成人股骨干骨标本30例,分别应用解剖法、CT二维法、CT三维法测量股骨颈前倾角,运用配对t检验比较两种方法测量股骨颈前倾角与解剖法相比的准确性,运用组内相关系数(ICC)评估两种方法测量股骨颈前倾角与解剖法的一致性。结果 解剖法、CT二维法、CT三维法测量的股骨前倾角分别为10.47±3.79°、12.82±4.50°、9.80±4.16°。CT二维法、CT三维法测得的股骨前倾角与解剖法测量的结果比较无显著差异(P<0.05)。 CT二维法与解剖法的ICC为0.872(95%CI, 0.770~0.932, P<0.05);CT三维法与解剖法的ICC为0.943 (95%CI, 0.855~0.977, P<0.05)。结论 CT二维法、CT三维法均能比较准确地测量股骨前倾角,CT三维法测量准确性相对更高。