共查询到20条相似文献,搜索用时 31 毫秒
1.
二维地形不平条件下均匀外电场的有限差分模拟 总被引:3,自引:0,他引:3
本文提出了一种求解二维不平地形条件下中梯装置电阻率法正演问题的有限差分算法。引入了数学上的某些新进展:采用三角形和矩形相结合的方式对求解区作剖分,以积分插值法导出有限差分格式,用多层网格修正迭代法求解线性方程组,从而得到了比较满意的模拟精度和计算速度。 相似文献
2.
基于一阶速度-应力波动方程,采用高阶交错网格有限差分数值模拟方法,对弹性及黏弹性TTI介质进行正演数值模拟。模拟时采用完全匹配层吸收边界条件(PML)消除边界反射。同时设计了层状介质模型、断层模型,通过模型的正演计算,得到了不同时刻的地震波波场快照及合成地震记录,分析其波场运动学及动力学特征。模拟结果表明,交错网格有限差分法可以很好地完成对复杂介质的波场模拟,具有较高的精度和可靠性。 相似文献
3.
交错网格有限差分算法以其高效、精确、实用等优点在地震波数值模拟中得到广泛应用。目前交错网格有限差分的精度已达到时间4阶、空间2N阶;然而在求空间三次导数时,差分格式实际上并未达到所谓的2N阶精度,而是采用了低阶的差分格式及差分系数,这样有利于提高大尺度空间正演时的计算效率;但从计算精度的角度考虑,有必要推导出准确的满足2N阶精度的交错网格有限差分格式及差分系数,以得到更高精度的正演结果。笔者利用Taylor公式展开首次推导出了可导函数任意次导数的任意偶数阶精度的差分近似式及相应的差分系数,从而完善了常规高精度交错网格有限差分算法。采用新推导的交错网格有限差分格式得到的正演波形与解析解进行了对比,证明了新推导的差分格式的正确性,并与常规差分格式的正演波形进行了比较,结果显示,新推导出的交错网格有限差分格式模拟结果稳定性好,精度更高。 相似文献
4.
贴体网格有限差分正演模拟算法不仅能够精确模拟任意起伏地形下的波场特征,且计算效率较高,是一种很有应用前景的处理西部复杂地表问题的方法;然而,目前求解波动方程时常用的同位网格和标准交错网格,在处理贴体网格起伏地表正演模拟时存在诸多问题。为此,将全交错网格引入到曲线坐标系下,避免了标准交错网格的插值误差和同位网格中奇偶失联引起的高频振荡现象,提高了模拟精度,减小了算法实现的复杂度。在自由边界条件实施时,采用牵引力镜像法计算速度分量,速度自由边界条件配合紧致交错差分格式更新应力分量,得到了较好的效果。随后,重点研究了贴体全交错网格正演模拟算法的影响因素,考虑了网格正交性、网格间距和网格拼接等的影响,并取得了如下认识:算法对网格的正交性没有过分要求;网格间距的突变会引起虚假反射的产生;不同类型的网格拼接对模拟结果不会造成明显的影响。 相似文献
5.
介绍了有限差分方法的原理,叙述了砂岩型铀矿体上氡气场的分布,并建立了氡气场分布的理论方程,通过有限差分算法求解氡气场分布方程,为砂岩型铀矿体上氡气场的正演研究提供了数值计算方法。针对砂岩型铀矿体上氡气场低值区对应矿体而高值异常场对应矿体两端的特点,以矿体端点高值异常峰值点为中心进行二维等步长剖分,得出差分方程,为问题的进一步研究提供一种算法思路。 相似文献
6.
为了加快大地电磁三维正演的求解速度,本文将一种新型的代数多重网格算法——聚集多重网格(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程序能够提高十几倍的计算速度,尤其适合大规模大地电磁三维正演问题。 相似文献
7.
《物探化探计算技术》2019,(5)
有限元求解大地电磁正演问题时,对研究区域的剖分常规做法是采用规则化的网格。但规则化网格在剖分地形、断层、褶皱等复杂模型会产生较大的几何离散误差。针对上述情况,这里采用非结构化的四边形网格对二维地电模型进行网格剖分,并与自适应有限元相结合,由剖分的粗网格出发,利用每个频点下网格单元的后验误差估计值指导网格的局部加密,优化网格质量和数量,从而提高正演模拟的精度。通过一维K型地电模型利用本文方法算出数值解与解析解进行对比分析,验证了自适应有限元法在求解大地电磁二维正演模拟中的有效性。并通过对断层、褶皱模型的正演模拟,分析了其大地电磁正演响应特征。 相似文献
8.
9.
常规变加密网格有限差分波动方程数值模拟方法采用水平分层加密网格,该网格剖分策略在适应地形起伏和近地表速度结构变化特征方面效果较差。针对该问题,提出一种起伏多重变加密网格有限差分波动方程数值模拟方法。该方法根据地形起伏和近地表低速层到高速层的速度分布特征进行网格加密;采用不同网格中的变系数差分格式离散声波方程,在保证波场模拟精度的同时兼顾计算效率;同时,为了进一步保障地表附近波场模拟的精度,地表附近最细网格中的差分格式不做降阶处理,针对位于地表以上的虚像点的波场值,提出一种融入自由地表边界条件的法向虚像外推法。算例分析验证该算法对速度模型不同区域进行的网格多次加密显著提高了计算效率,以黄土塬实际模型为例,耗时为常规1 m×1 m网格耗时的43.3%,并可达到和细网格基本一致的模拟精度,模拟误差控制在10-12范围内,同时表现出很好的近地表散射压制和边界吸收效果,且算法能稳定地适应实际复杂地表介质。 相似文献
10.
《物探化探计算技术》2015,(3)
为了解决网格方法面临的剖分困难,将基于径向基函数的配点型无单元法用于重磁二维正演。主要研究无单元法的基本原理,包括配置-中心节点的设置、形状参数的选取、边界条件的处理、微分方程的离散、矩阵方程的求解等一系列具体问题。通过重磁二维正演结果表明:在相同的剖分尺寸下,无单元法的计算精度高于有限元、有限差分正演;相比于网格法,无单元法还具有节点设置灵活、不用网格剖分、程序编写容易等优点。 相似文献
11.
12.
《物探与化探》2017,(2)
为优化二维各向同性介质中弹性波频率域正演时阻抗矩阵的结构,减小正演所需内存,提高正演效率,在25点差分格式的基础上进行适当的简化,得到了二维弹性波频率域17点差分格式。该格式重新计算了弹性波中偏微分项和加速项的差分算子,减少了计算过程中的网格节点需求,构造了优化阻抗矩阵后的频率域正演矩阵方程;推导了纵波和横波相速度的频散公式,给出了不同泊松比条件下的频散曲线,得到了相速度误差控制范围±1%时每一横波波长内网格数需求。通过对比频散曲线和数值模拟时得到的波场快照及检波点处U、V分量,验证了17点差分格式与25点差分格式相比,具有稍严格的网格间距需求、相当的计算精度、略少的计算时间和更小的阻抗矩阵带宽等特点。 相似文献
13.
赵晓博 《物探化探计算技术》2011,33(5):517-521,463
利用Delaunay三角化这种网格非结构化方法。通过编程实现了二维模型的非结构化三角形网格剖分,并编写了中心回线法瞬变电磁2.5维有限元正演程序。与前人计算结果对比,在取得相同计算精度的情况下,与结构化网格相比,非结构化网格所需网格和节点数量大大减少,计算效率更高。通过将非结构化网格法引入到瞬变电磁2.5维正演模拟中,实现了对复杂二维地电模型的有限元数值模拟,提高了现有有限元算法的应用范围。 相似文献
14.
15.
为了改进计算区域离散化问题,本文利用自适应非结构化网格有限单元法求解二维地电结构下大地电磁场满足的加权余量表达式。在有限元求解电磁场的过程中,网格剖分越精细、计算精度越高,计算量也会越大。此外,结构化网格难以适应任意地形以及复杂地质构造。而自适应非结构化网格在电性变化剧烈的区域会自动加密,在电性缓变的区域则生成粗疏的网格,从而优化网格质量与数量。因此,文中引入COMSOL Multiphysics软件,以实现若干地电模型的构建及非结构化自由四边形单元网格化。将网格数据信息导入本文算法,计算大地电磁场响应,并与解析解及数值解对比。结果表明,基于非结构化网格的正演模拟精度高、适应性强,为计算区域网格化提供了新的方法。 相似文献
16.
17.
18.
MT二维正演的混合剖分算法及其在线源CSAMT正演计算中的应用 总被引:1,自引:1,他引:1
在矩形有限元和三角形有限元方法大地电磁测深(MT)二维正演计算的基础上,设计了矩形单元粗剖分、程序自动进行三角形单元细剖分的混合剖分算法程序,使得MT二维模型正演的计算速度大为提高。该方法同时被应用于线电流源变频电磁测深(CSAMT)TE极化方式的二维正演计算。文章给出了在不同模型上计算的MT和线源CSAMT的响应,分析了计算CSAMT响应时,计算点离源的距离对视电阻率和相位值计算的影响,比较了MT和CSAMT网格剖分的特点,提出了计算CSAMT响应的改进的边界条件设计 相似文献
19.
由于地质体和矿体的形态非常复杂,使用长方体网格离散建立正演模型时可能和真实情况有很大差别,因此计算结果可靠性差。本文提出一种基于约束Delaunay网格剖分的方法对地质体进行离散并进行重力建模,在模型边界等复杂区域使用网格自适应加密技术,将三维地质体离散为有限个四面体;并详细推导出针对四面体网格的重力正演公式,实现了基于约束Delaunay网格剖分技术的三维重力数值模拟;最后,针对一个合成数据模型,将计算解与解析解对比。结果表明,细化网格的模拟结果比粗糙网格更好,满足数值模拟的精度要求。将该方法应用到金川矿区实际地质体建模中,根据局部需要,建立各处网格密度不均匀的三维模型,并计算该模型的地表重力场,而后对比模拟数据与实测数据,结果表明Delaunay网格建模方法具有很强的适用性,能够模拟复杂的地质体重力异常。 相似文献
20.
基于球坐标系下有限差分的地磁测深三维正演 总被引:2,自引:0,他引:2
为了计算全球尺度电磁感应的响应,本文介绍地磁测深频率域三维正演。正演算法采用球坐标系下的交错网格有限差分方法,从Maxwell方程的积分形式出发,采用PARDISO对离散后的方程组求解,避免了迭代求解的散度校正。为了验证本文结果的正确性和精度,与前人的有限元和有限差分方法进行了对比,一维层状模型的三维交错网格有限差分数值结果和解析解相对误差小于5%,双半球模型的计算结果与前人的计算结果完全吻合。三维"棋盘模型"计算表明磁场分量对异常体的大小和位置具有很好的分辨能力。 相似文献