共查询到20条相似文献，搜索用时 31 毫秒
1.
以双异性点电流源电流场二维地电断面为例，着重讨论了电阻率展析成像反演中有限单元法数值模拟边界条件的选择、网格剖分以及模拟精度等问题，同时也讨论了成像反演中最优化及演方法及拟合参数的评价等。文中提出的新的网格剖分方法和选择的边界条件可有效地提高计算速度和数值模拟的精度。 相似文献
2.
非均匀介质中交错网格高阶有限差分数值模拟 总被引：5，自引：0，他引：5
地震波场的数值模拟一直是地球物理学的一个重要的研究领域,而在数值正演模拟方法的研究中,计算精度和计算效率是评价该方法有效性及优越性的二个关键问题。这里从一阶速度—应力弹性波动方程出发,着重介绍如何构造离散化模型的网格,如何求解空间导数,如何选取边界条件等内容,从而更有效地提高数值计算的精度与计算效率。文中构造了不同类型的介质模型,并在交错网格中,利用高阶有限差分模拟非均匀介质的波场传播。模拟结果表明,该方法实现简单,具有很好地稳定性和较高的精度,能够直观、高效地反映出介质中波场的传播规律。 相似文献
3.
4.
5.
利用有限元方法进行大地电磁正演数值模拟时,由于是在有限网格区域上的数值计算,模拟计算时的网格边界为截断边界,而有限元数值模拟时的大地电磁场边界条件需要在足够远处才能够满足,所以截断边界的存在可能会使大地电磁正演模拟的边界条件无法满足,致使对计算结果和计算精度产生影响。利用有限元二维正演程序,在网格边界处加载一维情况下的大地电磁场,然后固定研究区域的网格剖分,并对一维地电模型和二维地电模型在改变有限元网格边界大小的情况下进行计算。在对一维模型进行模拟计算时,截断边界对边界条件没有影响,边界条件自然满足。而对二维模型进行模拟计算时,截断边界的存在对计算结果有较大影响。利用趋肤深度作为有限元网格边界变化的量度,通过改变网格边界大小,对不同的二维地电模型进行计算比较,总结出适合大地电磁有限元正演模拟的参考网格边界。 相似文献
6.
殷文 《吉林大学学报(地球科学版)》2008,38(1):144-151
在弹性波频率空间域有限差分数值模拟方面,差分网格及边界条件是影响弹性波模拟成功与否的关键,为了压制数值模拟中的网格频散,采用25点有限差分算子,建立了有限差分矩阵方程,且借鉴匹配层衰减边界条件思想,设计了弹性波频率空间域有限差分数值模拟算法。由于采用高阶有限差分法来提高差分格式的精度,将会导致计算量显著增加,为此,对频率空间域有限差分弹性波数值模拟方法,采用流水线技术与分治策略进行了并行算法研究,提高了计算效率,使得在合理的计算时间内更精确地模拟弹性波在弹性介质中的传播过程。 相似文献
7.
针对地下水有限元数值模拟中区域三角网格剖分复杂难以处理的情况,提出适合其特点的Delaunay三角网格自动剖分方法,并对含有多个参数分区的含水层进行网格剖分,最后利用遗传算法和有限元程序相耦合来反演含水层水文地质参数。结果表明此方法可大大简化地下水数值模拟的前处理工作,并能提高有限元网格剖分的有效性和准确性,从而得到令人满意的数值模拟结果。 相似文献
8.
大地电磁正演的精度对后续反演意义重大。决定大地电磁有限元正演精度的因素,主要有插值函数、方程组的求解误差,以及网格的剖分程度。而网格剖分的合适与否,为有限元的求解提供了先决条件。因此,从边界条件和横向网格以及纵向网格三方面,讨论了网格剖分对大地电磁正演精度的影响。经研究表明,对于一维和左右对称的二维介质,左右和下边界都无需放置在无穷远处,并且横向网格对精度影响不大。而由于纵向网格与近地表的最小趋肤深度密切相关,所以地表网格不得大于1/4个趋肤深度。 相似文献
9.
10.
基于一阶速度-应力波动方程,采用高阶交错网格有限差分数值模拟方法,对弹性及黏弹性TTI介质进行正演数值模拟。模拟时采用完全匹配层吸收边界条件(PML)消除边界反射。同时设计了层状介质模型、断层模型,通过模型的正演计算,得到了不同时刻的地震波波场快照及合成地震记录,分析其波场运动学及动力学特征。模拟结果表明,交错网格有限差分法可以很好地完成对复杂介质的波场模拟,具有较高的精度和可靠性。 相似文献
11.
常规变加密网格有限差分波动方程数值模拟方法采用水平分层加密网格,该网格剖分策略在适应地形起伏和近地表速度结构变化特征方面效果较差。针对该问题,提出一种起伏多重变加密网格有限差分波动方程数值模拟方法。该方法根据地形起伏和近地表低速层到高速层的速度分布特征进行网格加密;采用不同网格中的变系数差分格式离散声波方程,在保证波场模拟精度的同时兼顾计算效率;同时,为了进一步保障地表附近波场模拟的精度,地表附近最细网格中的差分格式不做降阶处理,针对位于地表以上的虚像点的波场值,提出一种融入自由地表边界条件的法向虚像外推法。算例分析验证该算法对速度模型不同区域进行的网格多次加密显著提高了计算效率,以黄土塬实际模型为例,耗时为常规1 m×1 m网格耗时的43.3%,并可达到和细网格基本一致的模拟精度,模拟误差控制在10-12范围内,同时表现出很好的近地表散射压制和边界吸收效果,且算法能稳定地适应实际复杂地表介质。 相似文献
12.
13.
数值模拟是研究煤层气藏工程的一种常规方法,建立模型时人们常忽略垂向网格精细程度对模拟的影响,无法准确反映出垂向上流体流动规律、气水分异现象以及压降漏斗展布等,对模拟结果有很大影响。为了研究垂向网格划分精度对煤层气藏数值模拟过程的影响,采用不同模拟器对垂向网格的精细程度进行模拟计算,运用渗透率等效方法、局部网格加密方法和模拟器自带压裂方法模拟煤层压裂缝,总共模拟了3种方法15套方案。结果表明,垂向网格划分精度对煤层气生产影响较大,当网格步长达到1.5m时计算结果较为精确,可以满足模拟需求。网格数量及网格步长的合理划分,能够更好地呈现煤层中气水分异现象,有助于分析气水流动状态,便于历史拟合和产量预测。 相似文献
14.
地下水数值模拟中人为边界的处理方法研究 总被引：10，自引：1，他引：9
边界条件预测是地下水数值模拟中的关键问题。本文针对人为边界受自然和人为因素影响,边界流量变化复杂难于预测的问题,讨论了人为边界条件的处理方式和实现方法,并将其应用到具有人为边界的二维潜水含水层系统。预测结果表明,受区内人工开采影响,边界流量呈逐年上升趋势,部分边界由流出边界转化为流入边界,结果符合地下水系统变化趋势,提高了预测精度。 相似文献
15.
沈金松 《物探化探计算技术》2002,24(1):46-53
自Pelton等发表第一篇直接电测深数据的二维反演文章以来，电测深数据的二维反演解释已广泛应用于资源勘探和环境监测，并取得一定效果。然而，由于该方法中加入了人为确定的平滑约束，其反演解释结果具有很大不确定性，使其在实际生产中普及应用受到限制。因此，如何在现有方法的基础上进一步改进完善成为该方法普及应用的关键。为了使反演解释结果既忠于实测数据，又合理地反映地下介质的变化，作者在本文中讨论了基于最大熵谱原理的剩余信息量最小准则实现电测深数据的二维自动反演。作者在文中还介绍了正演计算和偏导数求取中，提高精度和计算效率的技巧。用三个数值模式算例对本文的自动反演方法进行了验证。数值模拟结果表明，基于剩余信息量电小准则的二维自动反演给出的解释结果基本不受人为因素影响，能反映地下介质电阻率的详细变化。 相似文献
16.
对于模拟示踪试验的溶质运移,传统方法采用细网格离散小尺度表皮层,存在网格数量多和计算时间长等问题。由于表皮层外缘的参数变化剧烈,即使采用细网格,仍然存在显著误差。因此发展基于完整井的新示踪试验溶质运移模型,提出新的瞬态Robin边界条件反映表皮层吸附/解吸溶质的影响,实现表皮层的无网格。应用拉普拉斯变换导出模型的解析解,通过有限元法建立基于非完整井的数值解。结果显示新瞬态Robin边界条件精准反映负表皮层的影响。对于求解反演问题,表皮层宽度的估计范围由0.45～0.54 m缩至0.47～0.48 m,纵向弥散度由0.6～10.0 m缩至6.4～7.7 m,参数估计的可靠度显著提升。表皮层的无网格节省97%的数值解计算时间。 相似文献
17.
针对混凝土质量检测的特点,采用弯曲射线进行超声层析成像。首先建立波速分布的网格模型,并引入最短路径法进行射线追踪,然后采用联合迭代重建技术（SIRT）反演成像。为使算法适应弯曲射线,对SIRT的计算步骤进行调整,增加外层迭代,以便对修正后的路径进行下一轮反演;为提高反演精度,根据射线传播距离、幅度衰减、频率突变等信息建立和调整信号可信度系数,并依据该系数在求取网格慢度修正值时进行加权平均;为缩短计算时间,动态调整网格节点个数,并进行双向最短路径追踪。经过数值模拟及对混凝土基桩的实测,证明该方法在计算准确性和精度上有明显改进,计算速度也可以满足工程应用的需要。 相似文献
18.
19.
潘海滨 《物探化探计算技术》2009,31(4):369-373
交错网格高阶差分解法是地震波场模拟的一种有效方法。较大的时间延拓步长可以提高效率,但可能引起算法不稳定,使模拟无法进行。因此,对于空间网格大小和时间步长的选择要适中。虽然震源子波的主频对模拟精度的影响不大,但适当提高震源子波的主频,可以提高地震波对于薄层的分辨率。FCT方法能够有效地压制数值频散,将FCT方法与交错网格有限差分法相结合,可以提高波场模拟的精度和运算效率。 相似文献
20.
贴体网格有限差分正演模拟算法不仅能够精确模拟任意起伏地形下的波场特征,且计算效率较高,是一种很有应用前景的处理西部复杂地表问题的方法;然而,目前求解波动方程时常用的同位网格和标准交错网格,在处理贴体网格起伏地表正演模拟时存在诸多问题。为此,将全交错网格引入到曲线坐标系下,避免了标准交错网格的插值误差和同位网格中奇偶失联引起的高频振荡现象,提高了模拟精度,减小了算法实现的复杂度。在自由边界条件实施时,采用牵引力镜像法计算速度分量,速度自由边界条件配合紧致交错差分格式更新应力分量,得到了较好的效果。随后,重点研究了贴体全交错网格正演模拟算法的影响因素,考虑了网格正交性、网格间距和网格拼接等的影响,并取得了如下认识：算法对网格的正交性没有过分要求;网格间距的突变会引起虚假反射的产生;不同类型的网格拼接对模拟结果不会造成明显的影响。 相似文献