共查询到20条相似文献，搜索用时 158 毫秒
1.
烃类储集层是一种复合多相介质，在固体颗粒的空隙中含有气体或液体. 研究弹性波在该类地层中的传播规律对于油气勘探开发，特别对于全波列声波测井有重要意义. 为了提高孔隙弹性介质数值模拟的计算效率，本文采用改进显式交错网格有限差分算法取代常用的空间域四阶和时间域二阶的速度 - 应力有限差分算法，算法的空间域为八阶、时间域为二阶. 虽然计算的时间步长略小于空间域四阶的情形，但高阶有限差分算法可以选择较粗糙的网格，因此补偿了计算的低效；同时高阶交错网格有限差分算法的空间频散性比低阶算法小. 利用该算法计算了一个两层模型的波场，同时还模拟了等效弹性和孔隙弹性模型中波的传播. 结果表明慢波及其影响明显，尽管慢波衰减很快，但被某一界面反射后，转换形成的P波和S波仍以正常的方式传播，且比慢波衰减小. 相似文献
2.
在数值模拟中,隐式有限差分具有较高的精度和稳定性.然而,传统隐式有限差分算法大多由于需要求解大型矩阵方程而存在计算效率偏低的局限性.本文针对一阶速度-应力弹性波方程,构建了一种优化隐式交错网格有限差分格式,然后将改进格式由时间-空间域转换为时间-波数域,利用二范数原理建立目标函数,再利用模拟退火法求取优化系数.通过对均匀模型以及复杂介质模型进行一阶速度-应力弹性波方程数值模拟所得单炮记录、波场快照分析表明:这种优化隐式交错网格差分算法与传统的几种显式和隐式交错网格有限差分算法相比不但降低了计算量,而且能有效的压制网格频散,使弹性波数值模拟的精度得到有效的提高. 相似文献
3.
地震波场数值模拟在地球物理勘探和地震学中具有重要的支撑作用.本文将组合型紧致差分格式用于声波和弹性波方程的数值模拟中.根据泰勒级数展开和声波方程,建立了位移场时间四阶离散格式,并将组合型紧致差分格式用于位移场空间导数的求取,然后对该差分格式进行了精度分析、误差分析、频散分析和稳定性分析.理论研究结果表明：①该差分格式为时间四阶、空间六阶精度,与常规七点六阶中心差分和五点六阶紧致差分相比,具有更小的截断误差和更高的模拟精度;②每个波长仅需要5.6个采样点,且满足稳定性条件的库郎数为0.792,可以使用粗网格和较大时间步长进行计算.所以该方法具有占用内存少、计算效率高和低数值频散等优势.最后,本文进行了二维各向同性完全弹性介质的声波和弹性波方程的数值模拟,实验结果表明本文提出的方法具有更高的计算精度,能够大幅度的节约计算量和内存需求,对于三维大尺度模型问题具有更好的适应性. 相似文献
4.
时间域的波场延拓方法在本质上都可以归结为对一个空间-波数域算子的近似.本文基于一阶波数-空间混合域象征,提出一种新的方法求解解耦的二阶位移弹性波方程.该方法采用交错网格,连续使用两次一阶前向和后向拟微分算子,推导得到了解耦的二阶位移弹性波方程的波场延拓算子.由于该混合域象征在伪谱算子的基础上增加了一个依赖于速度模型的补偿项,可以补偿由于采用二阶中心差分计算时间微分项带来的误差,有效地减少模拟结果的数值频散,提高模拟精度.然而,在非均匀介质中,直接计算该二阶的波场延拓算子,每一个时间步上需要做N次快速傅里叶逆变换,其中N是总的网格点数.为了减少计算量,提出了交错网格低秩分解方法;针对常规有限差分数值频散问题,本文将交错网格低秩方法与有限差分法结合,提出了交错网格低秩有限差分法.数值结果表明,交错网格低秩方法和交错网格低秩有限差分法具有较高的精度,对于复杂介质的地震波数值模拟和偏移成像具有重要的价值. 相似文献
5.
《地球物理学进展》2015,(5)
交错网格有限差分方法已经被广泛应用到数值模拟和地震波传播的研究中.传统交错网格有限差分方法中,一阶空间导数的高阶差分系数是通过Taylor级数展开求取的,这种表示空间导数的方法会导致数值频散的产生.本文针对时间二阶空间十阶交错网格有限差分算法,采用最小二乘法通过改变积分区间求取一系列一阶空间导数的差分系数,分析该差分系数和传统方法求取的差分系数的频散关系.选取效果最佳的最小二乘法进行数值模拟,并与传统方法相比较.数值频散分析和弹性波场模拟分析表明:介质弹性参数和离散参数相同的情况下,采用最佳积分区间的最小二乘法更能有效地压制数值频散,比Taylor级数展开法具有更高的数值模拟精度. 相似文献
6.
本文对董良国、马在田等人发表在<地球物理学报>上的一阶弹性波方程交错网格高阶差分解法一文进行了差分格式的改进.发现原文对一阶空间导数采用2N阶差分精度,而对高阶导数却采用简单的中心差分,即对同一方程中的各阶导数的处理存在有的采用高阶格式、有的采用低阶格式的不一致性问题,本文的改进点正是消除了这种不一致性,建立了时间和空间均为四阶精度的差分格式.该方法具有编程简便、易于计算机实现和精度较高的特点. 相似文献
7.
地震波场数值模拟是理解地震波在地下介质中的传播特点，帮助解释观测数据的有效手段，而提高计算精度和运算效率是所有波场数值模拟方法研究所追求的目标.有限差分技术是求解波动方程计算效率最高、应用最为广泛的方法之一.但传统的有限差分技术计算过程中的数值频散问题影响了该技术的计算精度与计算效率.本文通过交错网格高阶有限差分技术与通量校正传输方法（Flux|corrected transport method，FCT）相结合, 对横向各向同性介质（Transverse isotropic medium，TI）一阶速度|应力弹性波动方程组进行了数值求解研究.波场快照数值模拟结果表明，本文研究的数值模拟方法与波动方程二阶有限差分方法、交错网格四阶有限差分方法相比，在压制网格数值频散方面有明显的优势，计算精度提高，而且可以利用较大的空间步长，提高计算效率. 相似文献
8.
叠前逆时偏移是目前成像精度最高的地震偏移方法之一,其实现过程中的一个重要步骤是数值求解全波方程,所以快速有效求解全波方程的数值算法对逆时偏移至关重要. 四阶近似解析辛可分Runge-Kutta (NSPRK) 方法是近年发展的一种具有高效率、高精度的数值求解波动方程的保辛差分方法, 能在粗网格条件下有效压制数值频散, 从而提高计算效率, 节省计算机内存需求量. 本文利用四阶NSPRK方法构造的基本思想,发展了具有六阶空间精度的NSPRK方法,并对新的六阶NSPRK方法进行了详细的稳定性和数值频散分析,以及计算效率比较和波场模拟. 同时将该方法用于声波叠前逆时偏移中, 得到一种时间上保辛、空间具有六阶精度、低数值频散、可应用大步长进行波场延拓并能长时计算的叠前逆时偏移方法,对Sigsbee2B模型进行了偏移成像, 并和四阶NSPRK方法、传统的六阶差分方法、四阶Lax-Wendroff correction (LWC) 方法进行了对比. 数值结果表明, 基于六阶NSPRK方法的叠前逆时偏移能得到更好的成像结果, 是一种优于四阶NSPRK方法、传统的六阶差分方法、四阶LWC叠前逆时偏移的方法, 尤其是在粗网格情况下具有更明显的优越性. 相似文献
9.
给出了任意倾斜各向异性介质中二维三分量一阶应力速度弹性波方程交错网格任意偶数阶精度有限差分格式及其稳定性条件,并推导出了二维任意倾斜各向异性介质完全匹配吸收层法边界条件公式和相应的交错网格任意偶数阶精度差分格式. 数值模拟结果表明,该方法模拟精度高,计算效率高,边界吸收效果好. 各向异性介质中弹性波波前面形态复杂, 且qP波波速不总是比qS波波速快. qS波波前面和同相轴的三分叉现象普遍, 且其同相轴一般不是双曲线型. 当TI介质倾斜时,3个分量上均能够观测到横波分裂现象, 而且各波形的同相轴变得不对称. 相似文献
10.
提高计算精度和运算效率是所有波场正演方法所追求的目标,本文通过将速度 （应力）对时间的奇数阶高阶寻数转化为应力（速度）对空间的导数,运用时间和空间差分精度 均可达任意阶的高阶差分法,通过交错网格技术,对一阶速度－应力弹性波方程进行了数值求 解．波场快照以及实际模型的正演结果表明,这种求解一阶弹性波方程的高阶差分解法,和 常规的差分法相比网格频散显著减小,精度明显提高,而且可以取较大的空间步长,提高计算 效率。 相似文献
11.
A mesh grading approach based on investigated lump method has been presented for simulating wave propagation in high velocity-contrast media. Unstructured fine grids are used to discretize the low wave-velocity medium in order to ensure the accuracy of numerical computation, and unstructured coarse grids are used for the high wave-velocity medium in order to substantially reduce the computational cost. On the interface, one coarse grid can match the fine grids of arbitrary odd number. The key feature of the proposed method is the constructions of investigated lumps on the interfaces of media. The transition zone, which is commonly used in the discontinuous grid scheme based on the staggered-grid finite-difference method, will not be used any more. Moreover, the computational instability that the discontinuous grid schemes frequently encountered does not arise in the proposed method. The comparisons with the analytical solutions and the application in studying the effects of sedimentary basin demonstrated that the mesh grading approach is a valid, accurate, convenient and flexible algorithm in simulating wave propagations in high velocity-contrast media with irregular interfaces. 相似文献
12.
传统有限差分方法在处理起伏地表时存在一些困难,而坐标变换法可将起伏地表映射为水平地表以克服此缺点.但同时,地下构造被变换得更加复杂,导致了波传播和成像的不准确.本文提出了一种分层的坐标变换方法,并应用到了弹性波逆时偏移中,此方法既可以克服起伏地表的影响,又可以不破坏地下构造.波场正向延拓、逆时延拓和分离是在辅助坐标系下完成的,而成像是在笛卡尔坐标系下完成的.通过对简单起伏模型和中原起伏模型的试算证明了本文提出方法的准确性.同时,对两种极端起伏地层高程不准确的情况进行测试可以看出:分层坐标变换逆时偏移方法的成像效果远好于传统坐标变换方法. 相似文献
13.
14.
15.
Numerical simulation in coupled elastic and poroelastic media is important in oil and gas exploration. However, the interface between elastic and poroelastic media is a challenge to handle. In order to deal with the coupled model, the first-order velocity–stress wave equations are used to unify the elastic and poroelastic wave equations. In addition, an arbitrary high-order discontinuous Galerkin method is used to simulate the wave propagation in coupled elastic–poroelastic media, which achieves same order accuracy in time and space domain simultaneously. The interfaces between the two media are explicitly tackled by the Godunov numerical flux. The proposed forms of numerical flux can be used efficiently and conveniently to simulate the wave propagation at the interfaces of the coupled model and handle the absorbing boundary conditions properly. Numerical results on coupled elastic–poroelastic media with straight and curved interfaces are compared with those from a software that is based on finite element method and the interfaces are handled by boundary conditions, demonstrating the feasibility of the proposed scheme in dealing with coupled elastic–poroelastic media. In addition, the proposed method is used to simulate a more complex coupled model. The numerical results show that the proposed method is feasible to simulate the wave propagation in such a media and is easy to implement. 相似文献
16.
Yongming Lu Qiancheng Liu Jianfeng Zhang Kai Yang Hui Sun 《Geophysical Prospecting》2019,67(5):1296-1311
With the progress in computational power and seismic acquisition, elastic reverse time migration is becoming increasingly feasible and helpful in characterizing the physical properties of subsurface structures. To achieve high-resolution seismic imaging using elastic reverse time migration, it is necessary to separate the compressional (P-wave) and shear (S-wave) waves for both isotropic and anisotropic media. In elastic isotropic media, the conventional method for wave-mode separation is to use the divergence and curl operators. However, in anisotropic media, the polarization direction of P waves is not exactly parallel to the direction of wave propagation. Also, the polarization direction of S-waves is not totally perpendicular to the direction of wave propagation. For this reason, the conventional divergence and curl operators show poor performance in anisotropic media. Moreover, conventional methods only perform well in the space domain of regular grids, and they are not suitable for elastic numerical simulation algorithms based on non-regular grids. Besides, these methods distort the original wavefield by taking spatial derivatives. In this case, a new anisotropic wave-mode separation scheme is developed using Poynting vectors. This scheme can be performed in the angle domain by constructing the relationship between group and polarization angles of different wave modes. Also, it is performed pointwise, independent of adjacent space points, suitable for parallel computing. Moreover, there is no need to correct the changes in phase and amplitude caused by the derivative operators. By using this scheme, the anisotropic elastic reverse time migration is more efficiently performed on the unstructured mesh. The effectiveness of our scheme is verified by several numerical examples. 相似文献
17.
交错网格高阶差分方法是一种在保持效率的前提下提高弹性波模拟精度的有效方法.本文将可变空间网格与变化的时间步长技术引入到交错网格高阶差分弹性波模拟中,提出一种空间网格可任意奇数倍变化与时间步长任意变化的交错网格高阶差分弹性波模拟方法.一系列数值试验表明,该方法能够在保证模拟精度的同时,通过有效降低空间与时间维度上的过采样来显著提高弹性波模拟的效率.同时,该方法还能够精细刻画含孔缝洞介质以及横向变化剧烈介质的局部细微结构,减小弹性波模拟误差,提高介质细微结构处的弹性波传播模拟精度. 相似文献
18.
Perfectly matched layer-absorbing boundary condition for finite-element time-domain modeling of elastic wave equations 总被引：2，自引：0，他引：2
The perfectly matched layer (PML) is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second-order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite-element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media. 相似文献
19.
When treating the forward full waveform case, a fast and accurate algorithm for modelling seismic wave propagation in anisotropic inhomogeneous media is of considerable value in current exploration seismology. Synthetic seismograms were computed for P-SV wave propagation in transversely isotropic media. Among the various techniques available for seismic modelling, the finite-difference method possesses both the power and flexibility to model wave propagation accurately in anisotropic inhomogeneous media bounded by irregular interfaces. We have developed a fast high-order vectorized finite-difference algorithm adapted for the vector supercomputer. The algorithm is based on the fourth-order accurate MacCormack-type splitting scheme. Solving the equivalent first-order hyperbolic system of equations, instead of the second-order wave equation, avoids computation of the spatial derivatives of the medium's anisotropic elastic parameters. Examples indicate that anisotropy plays an important role in modelling the kinematic and the dynamic properties of the wave propagation and should be taken into account when necessary. 相似文献
20.
A hybrid method combining finite element and 4th-order finite difference techniques is developed to model SH and P-SV seismic wave propagation in a 2D elastic medium with irregular surface topography. Both the classic staggered grid finite difference scheme and the partially staggered grid scheme are tested. The accuracy of the hybrid method is studied by comparison with a semi-analytical and another numerical method. Subsequently, to study the amplification, numerical simulations of seismic wave propagation in a series of hills are carried out and compared with the single-hill case. Depending on the position of the source in relation to the topography, the ratio between the heights and lengths of the hills or the ratio between the lengths of the hills and the wavelength, the presence of several hills as opposed to a single one can increase the amplification effect due to topography. This study highlights the fact that, when evaluating topographic site effects, surrounding topography must be taken into account in addition to local topography. 相似文献