共查询到20条相似文献，搜索用时 31 毫秒
1.
利用有限元方法进行大地电磁正演数值模拟时,由于是在有限网格区域上的数值计算,模拟计算时的网格边界为截断边界,而有限元数值模拟时的大地电磁场边界条件需要在足够远处才能够满足,所以截断边界的存在可能会使大地电磁正演模拟的边界条件无法满足,致使对计算结果和计算精度产生影响。利用有限元二维正演程序,在网格边界处加载一维情况下的大地电磁场,然后固定研究区域的网格剖分,并对一维地电模型和二维地电模型在改变有限元网格边界大小的情况下进行计算。在对一维模型进行模拟计算时,截断边界对边界条件没有影响,边界条件自然满足。而对二维模型进行模拟计算时,截断边界的存在对计算结果有较大影响。利用趋肤深度作为有限元网格边界变化的量度,通过改变网格边界大小,对不同的二维地电模型进行计算比较,总结出适合大地电磁有限元正演模拟的参考网格边界。 相似文献
2.
3.
为提高电磁波场数值计算精度,对时间域有限差分法与时间域伪谱法进行了对比研究。时间域有限差分法是一种目前流行的电磁场时域数值计算方法,已被广泛应用于求解与时间有关的偏导数方程。对于大规模数值计算,时间域有限差分法需要较多的内存空间。时间域伪谱法基于时间域有限差分法,该方法使用快速傅里叶度换来计算麦克斯韦方程中的空间导数。由于傅里叶变换的准确性,时间域伪谱法使计算精度提高,数值计算时所需的格子数大大减少,这极大地节省了计算机内存空间,适合于大规模正反演问题的数值计算。 相似文献
4.
基于一阶速度-应力波动方程,采用高阶交错网格有限差分数值模拟方法,对弹性及黏弹性TTI介质进行正演数值模拟。模拟时采用完全匹配层吸收边界条件(PML)消除边界反射。同时设计了层状介质模型、断层模型,通过模型的正演计算,得到了不同时刻的地震波波场快照及合成地震记录,分析其波场运动学及动力学特征。模拟结果表明,交错网格有限差分法可以很好地完成对复杂介质的波场模拟,具有较高的精度和可靠性。 相似文献
5.
三维磁异常傅里叶正演,能够计算整个场源区域内与磁化率分布网格一致的三维磁异常。通过三维傅里叶变换推导了长方体三维磁谱表达式,当源体埋藏较深或者位于反演区域边缘时,标准FFT正演异常场由于强制周期化,边界震荡等原因,使得正演结果发生较大的畸变,为了减少标准FFT算法引起的误差,引入了3DGauss-FFT技术用于三维磁异常频率域正演。通过简单的模型正演验证,从计算时间、计算精度以及内存需求上与空间域算法及标准3DFFT算法进行比较,结果表明,3DGauss-FFT磁异常正演相比于标准FFT算法在计算精度上有很大提升,显著降低了标准3DFFT由于自身缺陷引起的误差,且在运行时间上,4点3DGauss-FFT磁异常正演算法相对于空间域算法降低了三个数量级,但内存需求有所增加,证明了3DGauss-FFT正演方法在磁异常正演方面的高效性以及准确性。 相似文献
6.
大地电磁法正演中多重网格法求解的广义傅里叶谱分析 总被引：1，自引：0，他引：1
为了准确预测和分析多重网格法用于大地电磁法正演计算的收敛效果,对多重网格法的收敛性进行了广义傅里叶谱分析。通常情况下,系数矩阵的傅里叶谱是复数,为了直观地判断、提取收敛信息,将谱转换到实数域。在实数域内,特征向量谱定量解释了二重网格法收敛慢的原因（最粗网格用Gauss-Seidel法求解）。传统的局部傅里叶谱分析没有考虑边界条件和模型参数的变化,在分析二重网格法求解收敛性时得出的渐进收敛估计与数值解偏差大。针对这一问题提出广义傅里叶谱分析,其结果与数值解接近（比如2-V（0,1）的广义傅里叶谱分析为0.706,真实值为0.710）。对五重网格法求解的高重广义傅里叶谱分析结果表明,随广义傅里叶谱分析分析重数的增加, 所求渐近收敛估计趋于收敛的数值结果。基于此得出高重广义傅里叶谱分析的经验公式。求得低重广义傅里叶谱分析,通过矫正近似得到高重广义傅里叶谱分析,其结果在本文2个例子有效。 相似文献
7.
Parker公式在起伏地层重力场正演方面以其简洁、快速等优点被广泛应用于地形改正以及界面反演等,但传统的Parker公式正演方法存在一定的缺陷,在起伏地层模型重力场的正演计算中,由于正演式中e的指数项数值稳定性较差,导致正演结果精度不高。在前人研究的基础上,对Parker公式进行了改进,通过增加上、下界面平均值的方法,提高正演的数值稳定性,改进后的正演精度得到大幅提升。在此基础上,提出了一种曲面观测算法,将观测面由传统的水平面推广到任意起伏面,实现了上下界面起伏地层在起伏观测面上的快速高精度重力正演计算和带地形的密度界面起伏模型的快速高精度重力正演模拟。为了保障数值精度,采用了Gauss-FFT算法,模型实验均取得很好效果。 相似文献
8.
9.
基于球坐标系下有限差分的地磁测深三维正演 总被引：2，自引：0，他引：2
为了计算全球尺度电磁感应的响应,本文介绍地磁测深频率域三维正演。正演算法采用球坐标系下的交错网格有限差分方法,从Maxwell方程的积分形式出发,采用PARDISO对离散后的方程组求解,避免了迭代求解的散度校正。为了验证本文结果的正确性和精度,与前人的有限元和有限差分方法进行了对比,一维层状模型的三维交错网格有限差分数值结果和解析解相对误差小于5%,双半球模型的计算结果与前人的计算结果完全吻合。三维"棋盘模型"计算表明磁场分量对异常体的大小和位置具有很好的分辨能力。 相似文献
10.
通过对声波方程采用二阶有限差分与四阶有限差分近似,进行数值模拟及波场分析,边界条件利用吸收边界和增加衰减带混合使用的方法。计算结果表明,四阶有限差分的精度明显高于二阶有限差分,并且边界吸收效果好。对于单炮记录难于模拟复杂模型的情况下,实现了单次叠加剖面和爆炸反射面的零炮检距剖面,并进行了偏移处理。地震剖面更加清晰地反映了复杂构造的信息,更好地实现了对复杂地质模型的正演模拟。 相似文献
11.
正演数值模拟是反演成像的基础。为了实现磁法勘探精细化反演成像与定量解释,本文利用把一个大问题分解为多个小问题的思路,提出一种基于泊松方程的高效、高精度空间波数混合域二度体磁异常数值模拟方法。该方法利用傅里叶变换把二度体磁位偏微分方程转换为一维常微分方程,并采用基于二次插值的一维有限单元法求解该方程,进而通过反傅里叶变换得到空间域磁异常。在模型算例中,分别设计截面为矩形的常磁化率和变磁化率二度体模型,针对本文算法的计算精度和计算效率进行了验证。模型算例结果表明：该算法计算精度高,相对误差绝对值均小于1%;计算速度快,网格节点剖分2 501×2 501的模型模拟时间为4.18 s;适用于任意复杂地形模型。 相似文献
12.
三维电阻率法对反演的精度和速度的要求越来越高,而正演是反演的基础,因此直流电阻率三维正演计算的速度和精度是三维电阻率反演实用化的关键。这里利用对称超松弛预条件共轭梯度法(SSOR-PCG),求解有限差分法离散生成的大型稀疏线性方程组,预条件矩阵的选择大大降低了系数矩阵的条件数,结合矩阵的一维非零元素压缩存储模式,使得正演计算速度得以提高,而内存占用量明显减小。在直流电阻率三维正演中采用异常场法,提高了电源点附近的解的精度。利用编制的有限差分正演程序,对两层模型、垂直接触带模型和低阻异常体模型进行了数值模拟,计算结果表明该算法是可行的,且可以明显提高正演计算的速度和精度。 相似文献
13.
交错网格有限差分算法以其高效、精确、实用等优点在地震波数值模拟中得到广泛应用。目前交错网格有限差分的精度已达到时间4阶、空间2N阶;然而在求空间三次导数时,差分格式实际上并未达到所谓的2N阶精度,而是采用了低阶的差分格式及差分系数,这样有利于提高大尺度空间正演时的计算效率;但从计算精度的角度考虑,有必要推导出准确的满足2N阶精度的交错网格有限差分格式及差分系数,以得到更高精度的正演结果。笔者利用Taylor公式展开首次推导出了可导函数任意次导数的任意偶数阶精度的差分近似式及相应的差分系数,从而完善了常规高精度交错网格有限差分算法。采用新推导的交错网格有限差分格式得到的正演波形与解析解进行了对比,证明了新推导的差分格式的正确性,并与常规差分格式的正演波形进行了比较,结果显示,新推导出的交错网格有限差分格式模拟结果稳定性好,精度更高。 相似文献
14.
15.
采用解析层元法对存在地下点热源的岩土工程问题进行解答。首先从热弹性力学三维问题的基本控制方程出发,利用拉普拉斯-傅里叶积分变换推导出其在变换域内单层介质及下卧半空间的解析层元;然后结合有限单元法原理组装得到总刚度矩阵,结合边界条件,得到其在变换域内的解答,最后应用拉普拉斯-傅里叶积分逆变换技术,得到物理域内的解。编制了相应的计算程序,算例分析表明：该结果与已有文献吻合较好,该方法求解层状半空间的热-力耦合响应问题具有较好的适用性和较高的精度;层状岩土介质体系中,热扩散系数对温度及地表隆起的变化过程影响显著,但对其初始值和最终稳定值影响不明显;分层特性对岩土介质温度分布及地表位移变化过程均有显著影响。 相似文献
16.
点源二维电法正演的有限单元法 总被引：3，自引：3，他引：3
点源二维任意地电断面电阻率法的正演计算是电法勘探中的一个重要问题。对于这个问题,我们已经实现了有限差分方法的计算,获得了较好的效果。有限单元法是解决电法正演计算又一重要途径,在1971年Coggon首先将有限单元法引进到电法正演的模拟计算中来,以后Rijo等人又发展了这一方法,有限单元法进行电阻率法正演计算的优点,是比较容易拟合实际复杂的地电断电和起伏不平的地形。在实现点源二维地电断面视电阻率异常正演计算的有限单元法中,我们将混合边界条件引入求解,并采朋了若干优化措施,使有限单元法计算视电阻率异常的精度和速度均有进一步提高,获得较好的计算效果。在文献中我们简单地介绍了这个方法,下面我们较详细 相似文献
17.
介绍了国内、外张量测量技术及数据处理技术的发展情况。把离散余弦变换与有限单元方法结合起来,进行重力张量的正演计算。首先,利用有限单元法求解重力场的一阶导数,然后利用离散余弦变换求解重力张量分量。这样既保留了有限单元法对复杂模型的正演优势,又避免了利用有限元法直接求解重力张量时,因多次插值引起的数据量过大的问题。模型试验证明,用基于DCT的有限元法进行重力张量的正演,不但具有很高的计算精度,而且相对于利用有限单元法直接正演,其计算速度得到了显著的提高。 相似文献
18.
19.
复电阻率法二维正演采用二次场算法,结合有限单元法进行离散化,然后引入ColeCole模型,完成数值模拟。通过理论模型试算验证了二维正演程序的正确性。当模型剖分较细时,多发射源、多频率的正演计算效率有待提高。复电阻率法二维数值模拟是按发射源或频率计算,且各发射源之间、各频率之间计算是相互独立且互不影响的,因此通过引入MPI并行计算来提高其计算效率,从而减少其正演所需的时间。采用串行、并行复电阻率法二维正演程序来计算同一理论模型,结果验证了该并行程序的正确性。对并行算法的效率进行分析表明,该并行算法是可靠的、高效的。 相似文献
20.
井间地震正演模拟技术是研究地震波在井间传播规律的重要手段之一,可以帮助认识井间地震的复杂波场.从二维井间地震波传播波动方程出发,结合初始、边界条件,推导出了交错网格任意偶阶精度差分格式,阐述了非均匀差分网格的实现方法.在此基础上,编制了跨孔波场交错网格有限差分正演程序,并应用该程序对半空间粘弹性跨孔模型进行了正演计算,得到了震源在不同位置时的波场快照及地震记录图.这对于跨孔模型波传播规律的研究,具有重要意义. 相似文献