共查询到19条相似文献，搜索用时 125 毫秒
1.
殷文 《吉林大学学报(地球科学版)》2008,38(1):144-151
在弹性波频率空间域有限差分数值模拟方面,差分网格及边界条件是影响弹性波模拟成功与否的关键,为了压制数值模拟中的网格频散,采用25点有限差分算子,建立了有限差分矩阵方程,且借鉴匹配层衰减边界条件思想,设计了弹性波频率空间域有限差分数值模拟算法。由于采用高阶有限差分法来提高差分格式的精度,将会导致计算量显著增加,为此,对频率空间域有限差分弹性波数值模拟方法,采用流水线技术与分治策略进行了并行算法研究,提高了计算效率,使得在合理的计算时间内更精确地模拟弹性波在弹性介质中的传播过程。 相似文献
2.
3.
为高效求解扩展型Boussinesq水波方程,建立了基于有限差分和有限体积方法的混合数值格式。将一维控制方程写为守恒形式,方程中通量部分采用有限体积方法求解,剩余部分采用有限差分方法求解。其中,有限体积方法采用Godunov类高分辨率格式,并结合HLL(Harten-Lax and van Leer)式黎曼问题近似解求界面数值通量,黎曼问题界面左右变量通过高精度状态插值方法(MUSCL)构筑。有限差分方法则采用具有二阶精度的中心差分公式进行。采用具有TVD(Total Variation Diminishing)性质的三阶龙格-库塔多步积分法进行时间积分。对数值模式进行了验证,数值结果同解析解或实验数据吻合良好。 相似文献
4.
交错网格有限差分算法以其高效、精确、实用等优点在地震波数值模拟中得到广泛应用。目前交错网格有限差分的精度已达到时间4阶、空间2N阶;然而在求空间三次导数时,差分格式实际上并未达到所谓的2N阶精度,而是采用了低阶的差分格式及差分系数,这样有利于提高大尺度空间正演时的计算效率;但从计算精度的角度考虑,有必要推导出准确的满足2N阶精度的交错网格有限差分格式及差分系数,以得到更高精度的正演结果。笔者利用Taylor公式展开首次推导出了可导函数任意次导数的任意偶数阶精度的差分近似式及相应的差分系数,从而完善了常规高精度交错网格有限差分算法。采用新推导的交错网格有限差分格式得到的正演波形与解析解进行了对比,证明了新推导的差分格式的正确性,并与常规差分格式的正演波形进行了比较,结果显示,新推导出的交错网格有限差分格式模拟结果稳定性好,精度更高。 相似文献
5.
常用的求解地下水渗流模型有限差分法，三角形单元Δikj内水文地质参数相同，以相同参数的三角形单元进行参数分区，这种方法，对某些情形可以达到解决地下水渗流场模拟和预报问题，但对非均质各向异性有一定的局限性。本文阐述了以三角形单元的棱边控制面积作为参数分区的最小单元的参数分区方法，建立了单元非等参有限差分方程，给出了实际应用例子。该方法可更准确刻画非均质各向异性问题，同时兼容以往的差分方程，可退化成一般有限差分格式。 相似文献
6.
我们已推出一种二维有限差分时间域（ＦＤＴＤ）偏移算法,有限差分方法的优点在于它能很地应用于各种电导率背景模型,比较ＦＤＴＤ算法和积分方程（ＩＥ）算法的电磁偏移计算结果表明,有限差分法结果比ＩＥ结果能更好地反映电导异常,对于含有榀 的数据也可提供稳定的图像。Ｍｉｎｄｃｏ的地球物理学家们,跨跃日本的Ａｗａｊｉ岛的Ｎｏｊｉｍａ断裂带,开展了ＴＤＥＭ法工作,并利用有限差分偏移方法测定了与断裂有关的导电破碎 相似文献
7.
通过对声波方程采用二阶有限差分与四阶有限差分近似,进行数值模拟及波场分析,边界条件利用吸收边界和增加衰减带混合使用的方法。计算结果表明,四阶有限差分的精度明显高于二阶有限差分,并且边界吸收效果好。对于单炮记录难于模拟复杂模型的情况下,实现了单次叠加剖面和爆炸反射面的零炮检距剖面,并进行了偏移处理。地震剖面更加清晰地反映了复杂构造的信息,更好地实现了对复杂地质模型的正演模拟。 相似文献
8.
基于球坐标系下有限差分的地磁测深三维正演 总被引：2，自引：0，他引：2
为了计算全球尺度电磁感应的响应,本文介绍地磁测深频率域三维正演。正演算法采用球坐标系下的交错网格有限差分方法,从Maxwell方程的积分形式出发,采用PARDISO对离散后的方程组求解,避免了迭代求解的散度校正。为了验证本文结果的正确性和精度,与前人的有限元和有限差分方法进行了对比,一维层状模型的三维交错网格有限差分数值结果和解析解相对误差小于5%,双半球模型的计算结果与前人的计算结果完全吻合。三维"棋盘模型"计算表明磁场分量对异常体的大小和位置具有很好的分辨能力。 相似文献
9.
非饱和土热湿耦合传输问题的数值解——理论及一维问题的解 总被引：1，自引：0，他引：1
本文介绍了求解非饱和土壤中热量和水分耦合传输问题的一种数值方法——积分有限差分方法(IFDM)。基于菲利普-迪弗瑞斯(Philip-De Vries)多孔介质热湿耦合流动模型,采用积分有限差分方法,编制了求解热湿传输问题的计算程序 HM1,该程序可用于求解各类工程中遇到的多孔介质一维传热传湿问题。文中还给出了计算结果与实测结果相比较的实例。 相似文献
10.
11.
贮量集中有限元方法在渗流计算中的优点 总被引：2，自引：0，他引：2
普通有限元方法和有限差分方法是渗流计算中常用的数值方法,而贮量集中有限元方法则很少被实际应用。为了使渗流计算者全面了解贮量集中有限元方法,并能够较多地应用这种方法,本文将渗流计算中的贮量集中有限元方法、普通有限元方法和有限差分方法进行了对比,具体比较和分析了有限元方程组的系数矩阵特征、计算工作量,通过具体的算例分析了它们的计算精度,由此指出:贮量集中有限元方法的计算精度与有限差分法基本相同,而远好于普通有限元法,它的计算工作量小于其他二种方法 相似文献
12.
13.
14.
15.
碎石化沥青加铺路面结构中不同层位材料性质差异显著,难以由同一种力学模型加以描述。针对这一问题,将沥青加铺层及碎石化混凝土视为连续及非连续颗粒介质,分别采用有限差分法(DEM)与颗粒离散单元法(PCC)模拟沥青加铺层及碎石结构层的力学行为,并采用耦合的方式,通过两域对应边界的相容条件与平衡条件的控制,实现了连续域与离散域的平滑过渡。作为应用,结合现场情况,将破碎混凝土简化成碎粒层、碎石表层和碎石底层,通过模拟计算,分析了破碎粒径不同组合条件下沥青加铺层对荷载的响应。结果表明,碎粒层粒径越大,沥青层的拉应力、应变越大;碎石表层粒径越小,即多锤头破碎机(MHB)破碎的越碎,原水泥混凝土面板的强度下降的越多,整体性就越差,沥青层底的拉应力、应变也就越大,因此要将破碎状态控制在一定程度为宜。 相似文献
16.
为了较好地模拟给定模型在随机源激发下产生的稳定波场信息,为微动信号的检验和数据特征分析等提供理论依据,采用交错网格有限差分方法模拟微动信号,并根据瑞雷波和勒夫波的产生机理,通过测区内震源的分布情况估测微动信号中面波的组成成分。首先建立交错网格有限差分法模拟微动信号的数值模型,然后给出估测微动信号中面波组成成分的算式,最后采用典型方法数值模拟微动信号,应用空间自相关方法提取模拟信号中面波的频散曲线进行分析对比。结果表明,交错网格有限差分方法模拟得到的微动信号更符合实测结果。 相似文献
17.
介绍了有限差分方法的原理,叙述了砂岩型铀矿体上氡气场的分布,并建立了氡气场分布的理论方程,通过有限差分算法求解氡气场分布方程,为砂岩型铀矿体上氡气场的正演研究提供了数值计算方法。针对砂岩型铀矿体上氡气场低值区对应矿体而高值异常场对应矿体两端的特点,以矿体端点高值异常峰值点为中心进行二维等步长剖分,得出差分方程,为问题的进一步研究提供一种算法思路。 相似文献
18.
19.
A number of phenomena and processes in geosciences can be summarized by second order partial differential equations. The major numerical methods for their solution include the classical finite difference method and the finite element method newly developed in the last two or three decades. Since 1977 the author has proved that for the Laplace and Poisson equations, these two methods are identical and are different only in the process of formulation. For transient problems, such as heat conduction in the earth and the groundwater and oil-gas unsteady flow in porous media, there are some differences in resulting linear algebraic euqations. In general, two methods give similar results, but when the time step is decreased to some extent, the resulting algebraic equation will be consistent with the anti-heat conduction equation rather than the original heat conduction equation. This is the reason why unrealistic potentials are produced by the finite element method. Such a problem can be overcome by using the 相似文献