共查询到20条相似文献，搜索用时 171 毫秒
1.
波动方程方法是解决地震正反演问题的基本工具之一.无单元法作为一种新兴的偏微分方程数值计算方法,已经在材料力学、热传导等领域取得了显著的成功.由于抛弃了单元的概念及采用滑动最小二乘的拟合方法,使得无单元法具有前处理简单、精度高、独立变量解高次连续等优点.本文首先介绍无单元法求解波动方程的原理,指出影响其精度的主要因素.在算例的基础上详细讨论了无单元法用于实际波动问题的效果,并进一步尝试利用无单元法进行地震波数值模拟和反演成像的研究.模型计算的结果表明,无单元法能够较好的处理地震模拟和成像问题,精度和稳定性是令人满意的. 相似文献
2.
无单元Galerkin法采用滑动最小二乘法拟合场函数,只需节点无需单元,具有前处理简单、精度高、解高次连续等优点,被用于求解探地雷达(GPR)正问题.本文从Maxwell方程出发,推导了GPR正演需满足的波动方程;详细介绍了滑动最小二乘法形函数的构造方法.针对EFGM不满足插值条件导致强加边界条件的处理变复杂的特性,采用罚因子法对强加边界条件进行了处理;同时为了消除EFGM进行GPR正演模拟时来自截断边界处的超强反射,采用透射边界条件把GPR波在截断边界处的反射波透射出去,进而压制了来自截断边界处的反射波.然后,编制了EFGM的GPR正演模拟Matlab程序,应用该程序对典型GPR地电模型进行了正演模拟,并把该正演剖面图与基于线性插值FEM正演剖面图进行了对比,结果表明了EFGM用于GPR正演计算的正确性及有效性,并且在相同节点数条件下,EFGM比矩形剖分的FEM的精度要高,更有利于指导雷达剖面的数据解译. 相似文献
3.
有限差分方法是波场数值模拟的一个重要方法,交错网格差分格式比规则网格差分格式稳定性更好,但方法本身都存在因网格化而形成的数值频散效应,这会降低波场模拟的精度与分辨率.为了缓解有限差分算子的数值频散效应,精确求解空间偏导数,本文把求解波动方程的线性化方法推广到用于求解弹性波方程交错网格有限差分系数;同时应用最大最小准则作为模拟退火(SA)优化算法求解差分系数的数值频散误差判定标准来求解有限差分系数.通过上述两种方法,分别利用均匀各向同性介质和复杂构造模型进行了数值正演模拟和数值频散分析,并与传统泰勒展开算法、最小二乘算法进行比较,验证了线性化方法和模拟退火方法都能有效压制数值频散,并比较了各个算法的特点. 相似文献
4.
5.
《地球物理学进展》2015,(5)
交错网格有限差分方法已经被广泛应用到数值模拟和地震波传播的研究中.传统交错网格有限差分方法中,一阶空间导数的高阶差分系数是通过Taylor级数展开求取的,这种表示空间导数的方法会导致数值频散的产生.本文针对时间二阶空间十阶交错网格有限差分算法,采用最小二乘法通过改变积分区间求取一系列一阶空间导数的差分系数,分析该差分系数和传统方法求取的差分系数的频散关系.选取效果最佳的最小二乘法进行数值模拟,并与传统方法相比较.数值频散分析和弹性波场模拟分析表明:介质弹性参数和离散参数相同的情况下,采用最佳积分区间的最小二乘法更能有效地压制数值频散,比Taylor级数展开法具有更高的数值模拟精度. 相似文献
6.
利用传统有限差分方法对基于Biot理论的双相介质波动方程进行数值求解时,由于慢纵波的存在,数值频散效应较为明显,影响模拟精度.相对于声学近似方程及普通弹性波方程,Biot双相介质波动方程在同等数值求解算法和精度要求条件下,其地震波场正演模拟需要更多的计算时间.本文针对Biot一阶速度-应力方程组发展了一种变阶数优化有限差分数值模拟方法,旨在同时提高其正演模拟的精度和效率.首先结合交错网格差分格式推导Biot方程的数值频散关系式.然后基于Remez迭代算法求取一阶空间偏导数的优化差分系数,并用于Biot方程的交错网格有限差分数值模拟.在此基础上把三类波的平均频散误差参数限制在给定的频散误差阈值和频率范围内,此时优化有限差分算子的长度就能自适应非均匀双相介质模型中的不同速度区间.数值频散曲线分析表明:基于Remez迭代算法的优化有限差分方法相较传统泰勒级数展开方法在大波数范围对频散误差的压制效果更明显;可变阶数的优化有限差分方法能取得与固定阶数优化有限差分方法相近的模拟精度.在均匀介质和河道模型的数值模拟实验中将本文变阶数优化有限差分算法与传统泰勒展开算法、最小二乘优化算法进行比较,进一步证明其在复杂地下介质中的有效性和适用性. 相似文献
7.
8.
电阻率法有限单元正演模拟中,采用第三类边界条件时为保证精度仍要求较大范围的计算域.无单元法为地球物理领域的新兴正演模拟方法,其计算效率低,但其中采用的移动最小二乘(MLS)形函数相比于有限单元法形函数具有良好的连续性,模拟精度高.本文将MLS形函数应用于电阻率法有限单元2.5维正演的第三类边界条件处理,提出电阻率法有限单元-移动最小二乘(FEM-MLS)耦合正演方法.通过不同正演方法的模型算例模拟结果对比,验证了本文算法的有效性,并讨论了各个参数选择对模拟结果的影响.本文数值模拟结果表明采用第三类边界条件时,在同等计算精度前提下,FEM-MLS耦合法相比于有限单元法可进一步缩小计算域并提高了计算效率,相比于采用较大计算域满足边界条件的有限单元法计算效率提高了约一倍,相比于采用相同小范围计算域的有限单元法平均精度提高了约一倍. 相似文献
9.
一般偏移算法是用反演算子通过解析方法求解.最小二乘偏移方法采用另一种思路,即采用数值方法,通过解一个线性离散反问题来索求解.这样我们试着寻找一个模型匹配地震数据并能表现出其某些特点来得到偏移图像.最小二乘法能减少偏移赝像,得到更精确的偏移效果.Kirchhoff算子在最小二乘偏移方法中应用较广,但需要较多的迭代次数,而且具有Kirchhoff偏移的缺点.本文把最小二乘方法运用到基于波长延拓的波动方程偏移方法中,为提高最小二乘偏移的效率,可采用效率较高的正传播算子和反传播算子.我们利用效率较高,能适应剧烈横向变速的傅立叶有限差分正传播和反传播算子来做叠后最小二乘偏移.数值实例表明,通过少数的共轭梯度法迭代,就能得到与真实模型差别很微小的偏移效果.对于傅式变换我们采用了数值软件FFTW,其变换速度比常规FFT算法一般要快六倍以上,进一步提高了效率.本文算法很容易在并行机上实现,这些特点在处理大型数据时大有裨益. 相似文献
10.
11.
THE APPLICATION OF PML BOUNDARY CODITIONS IN THE SIMULATION OF SEISMIC WAVE BY THE HIGH-ORDER STAGGERED-GRID FINITE DIFFERENCES AT THE TRANSVERSELY ISOTROPIC MEDIA 
下载免费PDF全文
下载免费PDF全文 In the realm of the numerical simulation, finite difference method and finite element method are more intuitive and effective than other simulation methods. In the process of simulating seismic wave propagation, the finite differences method is widely used because of its high computational efficiency and the advantage of the algorithm is more efficient. With the demand of precision, more and more researchers have proposed more effective methods of finite differences, such as the high-order staggered-grid finite differences method, which can restore the actual process of wave propagation on the premise of ensuring accuracy and improving the efficiency of operation. In the past numerical simulation of seismic wave field, different models of isotropic medium are mostly used, but it is difficult to reflect the true layer situation. With the research demand of natural seismology and seismic exploration, the research on anisotropic media is more and more extensive. Transversely isotropic(TI)media can well simulate the seismic wave propagation in the formation medium, such as gas-bearing sandstone, mudstone, shale et al., the character of TI media is reflected by introducing the Thomsen parameters to reflect its weak anisotropy of vertical direction by using Thomson parameter. Therefore, studying the process of seismic wave propagation in TI media can restore the true information of the formation to the greatest extent, and provide a more reliable simulation basis for the numerical simulation of seismic wave propagation. In the geodynamic simulation and the numerical simulation of the seismic wave field, under the limited influence of the calculation area, if no boundary conditions are added, a strong artificial boundary reflection will be generated, which greatly reduces the validity of the simulation. In order to minimize the influence of model boundaries on the reflection of seismic waves, it is often necessary to introduce absorbing boundary conditions. At present, there are three types of absorption boundary conditions: one-way wave absorption boundary, attenuation absorption boundary, and perfectly matched layer(PML)absorption boundary. In terms of numerical simulation of seismic waves, the boundary absorption effect of PML is stronger than the first two, which is currently the most commonly used method, and it also represents the cutting-edge development direction of absorption boundary technology. The perfectly matched layer absorbing boundary is effectively applied to eliminating the reflective waves from model boundaries, but for transversely isotropic medium, the effect of the absorbing is not very well. For this reason, the elastic dynamic wave equations in transversely isotropic media are derived, and we describe a second-order accurate time, tenth-order accurate space, formulation of the Madariaga-Virieux staggered-grid finite difference methods with the perfectly matched layer(PML)are given. In addition, we have established vertical transversely isotropic(VTI)media and arbitrary inclined tilted transversely isotropic(TTI)media models, using a uniform half-space velocity model and a two-layer velocity model, respectively. By combining the actual geoscience background, we set the corresponding parameters and simulation conditions in order to make our model more research-oriented. When setting model parameters, different PML thickness, incident angle, source frequency and velocity layer models were transformed to verify the inhibition of boundary reflection effect by PML absorption boundary layer. The implementations of this simulation show that the formula is correct and for the transversely isotropic(TI)media of any angular symmetry axis, when the thickness of the PML layer reaches a certain value, the seismic wave reflection effect generated by the artificial boundary can be well suppressed, and the absorption effect of PML is not subject to changes in incident angle and wave frequency. Therefore, the results of our study indicate that our research method can be used to simulate the propagation process of seismic waves in the transversely isotropic(TI)media without being affected by the reflected waves at the model boundary to restore the actual formation information and more valuable geological research. 相似文献
12.
Introduction There are various factors affecting the propagation of earthquake waves. Usually these fac- tors include the earthquake source, the path and the site effects. Simulations of the earthquake wave propagation and the strong ground motion are an adequate way to exhibit these factors (Virieux, 1986; Levander, 1988; Frankel and Vidale, 1992; Frankel, 1993; Graves, 1993, 1996a; Olsen et al, 1995; Pitarka et al, 1996; Pitarka, 1999). The most popular numerical method in use is the gri… 相似文献
13.
14.
引入计算空气声学领域的选择性滤波同位网格有限差分算法（SFFD法）用于二维地震波数值模拟．SFFD法使用经过优化的11点DRP同位网格差分格式,对空间一阶导数进行离散近似,同时采用选择性滤波方法来消除同位网格差分所产生的格点高频振荡,它既提高了数值模拟的精度, 又保证了求解过程的稳定性．数值实验结果表明,SFFD法能够达到O(Delta;x8, Delta;t4)阶交错网格算法同样的精度,同时该方法还具有很强的适应性,能够应用于存在着强泊松比差异的介质模型中,完整地模拟地震波传播过程中各类型的波场,并且对复杂非均匀介质的适应能力也很好．此外,由于避免了交错网格算法在曲线坐标系和一般各向异性介质的数值模拟时所需进行的复杂的插值运算, SFFD法在这些问题上也有着很好的应用前景． 相似文献
15.
针对目前广泛用于解决工程地质问题的浅层地震勘探特点,通过对浅层小尺度介质进行地震波数值模拟,揭示地震波在不同浅层介质模型中的波场特征和传播规律.为浅层地震勘探野外数据采集前合理选择观测系统提供指导性的方案,并能对地震勘探成果解释的合理性给予客观评价.在分析各种数值模拟方法特点的基础上,选择具有精度高、易模拟复杂结构的有... 相似文献
16.
17.
针对粘弹性裂隙型单斜介质,本文应用高阶交错网格差分法,对其传播快照以及多方位地面三分量记录进行模拟,结果证明了裂隙填充物的性质引起的各向异性与粘弹性性质对单斜介质波场有明显的影响,总结了不同观测方位地震记录呈现出的规律性,这为进一步了解波场在介质中的传播情况,提出更合理的理论模型,提供一种尝试. 相似文献
18.
间断Galerkin有限元法（DG-FEM）作为一种有效的高阶有限元法受到了国内外学者的广泛关注.本文基于任意高阶间断Galerkin有限元法对弹性波方程进行空间离散,并将离散后所得的非齐次线性常微分方程系统齐次化,最后结合针对齐次问题的强稳定性保持龙格库塔（SSP Runge-Kutta）算法,将DG-FEM推广至时间任意高阶精度.另外,借鉴近最佳匹配层（NPML）的思想,基于复频移（CFS）拉伸坐标变换推导了一种新的PML吸收边界条件（简称为CFS-NPML）,该CFS-NPML能够与DG-FEM算法很好地结合,形成有效的起伏地表地震波传播数值模拟技术.数值试验结果表明,DG-FEM具有高阶精度,可以适应任意复杂起伏地表和复杂构造情况下的弹性波传播数值模拟.同时,CFS-NPML对包括面波等震相的人为边界反射都具有良好的吸收效果. 相似文献
19.