共查询到20条相似文献，搜索用时 62 毫秒
1.
本文以基于改进BISQ模型的二维双相各向同性介质一阶速度-应力方程为基础，推导出了曲线坐标系下对应的方程，然后采用低频散、低耗散的同位网格MacCormack有限差分法来离散方程，并采用紧致的单边MacCormack差分格式结合牵引力镜像法来施加自由地表边界条件，实现了地震波场数值模拟.曲线网格有限差分法采用贴体网格来描述自由表面，地表的网格线紧贴地形，避免了台阶近似造成的数值散射.数值模拟结果表明，在双相介质起伏自由地表和分界面处，各类波型复杂的反射透射规律可以清晰展现，曲线网格有限差分法可以精确地解决地震波在含起伏地表的双相各向同性介质中的传播问题. 相似文献
2.
基于应力、速度混合变量弹性波方程及任意四边形网格差分算子,给出了交错计算应力及速度的非规则网格弹性波应力一速度差分法该方法融合了有限元法能适应复杂形状边界及差分法无需计算刚度阵的特点,具有较高的计算精度,所需计算机存储空间较少,计算效率也很高．基于积分平衡方程引入了任意形状自由表面的边界条件,且通过局部滤波改善了自由表面边界条件的稳定性,使得该方法可应用于考虑地表形状影响的地震波数值模拟 相似文献
3.
为了便于研究双相介质固流相混合弹性波场中纵横波波场的传播规律,提出了基于交错网格的Biot双相各向同性介质弹性波动方程高精度波场分离正演数值模拟方法.采用高阶交错网格有限差分法来构建一阶双曲型双相各向同性介质弹性波动方程正演算子实现波场正演,并在每一步递推过程中,分别计算出同相和流相分量相应的散度场(纯纵波场)和旋度场... 相似文献
4.
The acoustic wave velocity varies with fluid saturation and pore-fluid distribution. We use a P-wave source and the staggered grid finite-difference method, with second-order accuracy in time and eighth-order accuracy in space, to simulate the acoustic wave field in a fractured medium that is saturated with a two-phase pore fluid (gas & water). Further, we analyze the variation of acoustic wave velocity with saturation for different pore-fluid distribution modes. The numerical simulation method is simple and yields accurate results. 相似文献
5.
In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equations. In this study, we compare two kinds of such wave equations: the first-order (velocity–stress) and the second-order (displacement–stress) separate elastic wave equations, with the first-order (velocity–stress) and the second-order (displacement–stress) full (or mixed) elastic wave equations using a high-order staggered grid finite-difference method. Comparisons are given of wavefield snapshots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corresponding first-order or second-order full elastic wave equations. These mixed equations are computationally slightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-component processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements. 相似文献
6.
地震波传播过程中,质点的振动不仅包括三个独立的平移部分,还包括三个独立的旋转部分.本文基于一阶速度-应力弹性波方程,采用分裂完全匹配层(SPML)的吸收边界条件,推导了时间导数二阶精度和空间导数高阶精度的交错网格有限差分格式的弹性波速度与应力各分量计算公式,模拟了各向同性介质中均匀模型和层状模型下的六分量波场,并对二维各向同性层状模型下的三个分量地震记录做高分辨率线性拉东变换得到各自的频散能谱.数值模拟分析结果表明:(1)旋转分量的能量要比平动分量弱的多;(2)在平动分量上,面波能量强,频率低,反射P波能量较强,反射S波能量稍弱;在旋转分量上,反射P波能量很弱,S波能量强;(3)与平动分量相比,旋转分量的频散能谱效果更好,能看到基阶和完整的高阶面波,即旋转分量能反映更多的地下介质信息. 相似文献
7.
三维起伏地表条件下的地震波走时计算技术是研究三维起伏地表地区很多地震数据处理技术的基础性工具.为了获得适应于任意三维起伏地表且计算精度高的走时算法,提出三维不等距迎风差分法.该方法采用不等距网格剖分三维起伏地表模型,通过在迎风差分格式中引入不等距差分格式、Huygens原理及Fermat原理来建立地表附近的局部走时计算公式,并通过在窄带技术中设定新的网格节点类型来获得三维起伏地表条件下算法的整体实现步骤.精度及算例分析表明:三维不等距迎风差分法具有很高的计算精度且能够适应于任意三维起伏地表模型. 相似文献
8.
We have pursued two-dimensional (2D) finite-difference (FD) modelling of seismic scattering from free-surface topography. Exact free-surface boundary conditions for the particle velocities have been derived for arbitrary 2D topographies. The boundary conditions are combined with a velocity–stress formulation of the full viscoelastic wave equations. A curved grid represents the physical medium and its upper boundary represents the free-surface topography. The wave equations are numerically discretized by an eighth-order FD method on a staggered grid in space, and a leap-frog technique and the Crank–Nicholson method in time.
In order to demonstrate the capabilities of the surface topography modelling technique, we simulate incident point sources with a sinusoidal topography in seismic media of increasing complexities. We present results using parameters typical of exploration surveys with topography and heterogeneous media. Topography on homogeneous media is shown to generate significant scattering. We show additional effects of layering in the medium, with and without randomization, using a von Kármán realization of apparent anisotropy. Synthetic snapshots and seismograms indicate that prominent surface topography can cause back-scattering, wave conversions and complex wave patterns which are usually discussed in terms of inter-crust heterogeneities. 相似文献
In order to demonstrate the capabilities of the surface topography modelling technique, we simulate incident point sources with a sinusoidal topography in seismic media of increasing complexities. We present results using parameters typical of exploration surveys with topography and heterogeneous media. Topography on homogeneous media is shown to generate significant scattering. We show additional effects of layering in the medium, with and without randomization, using a von Kármán realization of apparent anisotropy. Synthetic snapshots and seismograms indicate that prominent surface topography can cause back-scattering, wave conversions and complex wave patterns which are usually discussed in terms of inter-crust heterogeneities. 相似文献
9.
We present a discrete modelling scheme which solves the elastic wave equation on a grid with vertically varying grid spacings. Spatial derivatives are computed by finite-difference operators on a staggered grid. The time integration is performed by the rapid expansion method. The use of variable grid spacings adds flexibility and improves the efficiency since different spatial sampling intervals can be used in regions with different material properties. In the case of large velocity contrasts, the use of a non-uniform grid avoids spatial oversampling in regions with high velocities. The modelling scheme allows accurate modelling up to a spatial sampling rate of approximately 2.5 gridpoints per shortest wavelength. However, due to the staggering of the material parameters, a smoothing of the material parameters has to be applied at internal interfaces aligned with the numerical grid to avoid amplitude errors and timing inaccuracies. The best results are obtained by smoothing based on slowness averaging. To reduce errors in the implementation of the free-surface boundary condition introduced by the staggering of the stress components, we reduce the grid spacing in the vertical direction in the vicinity of the free surface to approximately 10 gridpoints per shortest wavelength. Using this technique we obtain accurate results for surface waves in transversely isotropic media. 相似文献
10.
利用传统有限差分方法对基于Biot理论的双相介质波动方程进行数值求解时,由于慢纵波的存在,数值频散效应较为明显,影响模拟精度.相对于声学近似方程及普通弹性波方程,Biot双相介质波动方程在同等数值求解算法和精度要求条件下,其地震波场正演模拟需要更多的计算时间.本文针对Biot一阶速度-应力方程组发展了一种变阶数优化有限差分数值模拟方法,旨在同时提高其正演模拟的精度和效率.首先结合交错网格差分格式推导Biot方程的数值频散关系式.然后基于Remez迭代算法求取一阶空间偏导数的优化差分系数,并用于Biot方程的交错网格有限差分数值模拟.在此基础上把三类波的平均频散误差参数限制在给定的频散误差阈值和频率范围内,此时优化有限差分算子的长度就能自适应非均匀双相介质模型中的不同速度区间.数值频散曲线分析表明:基于Remez迭代算法的优化有限差分方法相较传统泰勒级数展开方法在大波数范围对频散误差的压制效果更明显;可变阶数的优化有限差分方法能取得与固定阶数优化有限差分方法相近的模拟精度.在均匀介质和河道模型的数值模拟实验中将本文变阶数优化有限差分算法与传统泰勒展开算法、最小二乘优化算法进行比较,进一步证明其在复杂地下介质中的有效性和适用性. 相似文献
11.
如何有效压制数值频散是有限差分正演模拟研究中的关键问题之一.近年来,许多学者对二阶声波方程的差分算子开展了大量的优化工作,在压制频散方面取得不错的效果.一阶压强-速度方程广泛用于研究地震波在地下变密度模型中传播规律,目前针对一阶方程的优化工作大多只是在空间差分算子上展开.本文在前人研究的基础上,推导出一阶声波方程中压强场与偏振速度场之间的解析关系,据此在传统交错网格基础上给出一种高精度的显式时间递推格式,该递推格式将时间差分与空间差分算子结合在一起,并采用共轭梯度法得到精确时间递推匹配系数,实现时空差分算子的同时优化.在编程实现算法的基础上,通过频散分析与三个典型模型测试表明:本文方法能够较为有效地压制时间频散与空间频散,提高数值计算精度;同时对复杂模型也有很好适用性. 相似文献
12.
本文提出了一种新的基于平均导数优化方法(average-derivative optimal method,简称ADM)的二维VTI介质qP波波动方程频率空间域二阶9点格式,这种新算法将二维VTI介质qP波波动方程中中心空间导数项的差分近似表示为正交方向上3个网格点的加权平均形式.通过最小二乘优化方法求取空间导数项和加速度项的加权优化系数从而使数值频散达到极小化,每个波长所需要的网格点数在1%的误差范围内仅为3.57个网格点数,而VTI介质常规9点差分格式在相同的误差范围内则需要约12个网格点数,新方法的计算精度明显提高.复杂BP2007 2D VTI海洋标准模型数值模拟结果也验证了本文VTI介质9点ADM算法的有效性和准确性. 相似文献
13.
在数值模拟中,隐式有限差分具有较高的精度和稳定性.然而,传统隐式有限差分算法大多由于需要求解大型矩阵方程而存在计算效率偏低的局限性.本文针对一阶速度-应力弹性波方程,构建了一种优化隐式交错网格有限差分格式,然后将改进格式由时间-空间域转换为时间-波数域,利用二范数原理建立目标函数,再利用模拟退火法求取优化系数.通过对均匀模型以及复杂介质模型进行一阶速度-应力弹性波方程数值模拟所得单炮记录、波场快照分析表明:这种优化隐式交错网格差分算法与传统的几种显式和隐式交错网格有限差分算法相比不但降低了计算量,而且能有效的压制网格频散,使弹性波数值模拟的精度得到有效的提高. 相似文献
14.
探讨了井巷隧道工程反射波法超前地质预报成像问题.从一阶速度－应力弹性波方程出发,推导了二维各向同性介质情况下弹性波逆时传播的高阶差分格式,实现了弹性波在数值空间中的逆时延拓.构建反射波法隧道超前探测中的断层、软弱夹层等介质模型,以反射波探测的正演记录作为初始条件,并从程函方程出发,采用逆时差分格式求取介质模型网格空间中各点的直达波旅行时作为弹性波逆时偏移的成像条件,实现多波多分量资料的逆时偏移.偏移结果表明,逆时偏移能够使隧道壁接收到的波场准确归位,提高隧道反射波超前探测的资料处理的精度. 相似文献
15.
波动方程数值模拟是深入研究地震波传播规律的有效方法.有限差分法因其方法简单、精度高而得到了广泛的应用.但其缺点是不能准确模拟具有复杂几何形态的物性界面.因而当遇到起伏地表或复杂构造时，求解精度低.为了准确模拟起伏地形、复杂构造和复杂介质条件下的地震波场，本文采用有限元法模拟二维声波方程.用三角形单元模拟地形和速度界面；把单元内的场和波速均看作单元上的线性函数，以适应复杂介质压制边角散射；采用吸收边界条件去除来自截断边界上的反射；采用集中质量矩阵和集中阻尼矩阵使得显式时间递推无需对矩阵求逆，提高了计算效率.对模型的计算表明该方法正确有效. 相似文献
16.
从各向同性介质中波场数值模拟的褶积微分算子法出发,推导出了各向异性双相介质中波场传播数值计算的褶积新算法.将常见的二阶微分Biot波动方程用等效的一阶速度—应力双曲方程表示,其中未知的波场向量包括固相和流体的速度分量和应力分量,由此对方程的时间项使用交错网格差分方法计算,而对空间项则采用褶积微分算法进行求解.对各向异性双相介质在单层介质模型和双层介质模型中的波场特征进行了研究.研究的结果显示,在两层介质分界面上当地震波产生反射时能观测到两类纵波和横波,并且在衰减系数大的介质里慢纵波很难见到. 相似文献
17.
各向异性介质纯P波方程完全不受横波的干扰,在一定程度上可以减缓由于介质各向异性引起的数值不稳定,本文推导了具有垂直对称轴的横向各向同性(VTI)介质纯P波一阶速度-应力方程.由于纯P波方程存在一个分数形式的伪微分算子,无法直接采用有限差分法求解.针对该问题,本文采用伪谱法和高阶有限差分法联合求解波动方程,重点分析了混合法求解纯P波一阶速度-应力方程的稳定性问题,并给出了混合法求解纯P波方程的稳定性条件.数值模拟结果表明纯P波方程伪谱法和高阶有限差分混合法能够进行复杂介质的正演模拟,在强变速度、变密度的地球介质中仍然具有较好的稳定性. 相似文献
18.
海洋勘探环境可以抽象为下伏固体与上覆流体相互耦合的介质，本文针对流-固边界耦合介质提出了一种高效、稳定的多参数（速度和密度）全波形反演方法.本文采用弹性波一阶位移-应力方程作为过渡层耦合声波压力方程与弹性波位移方程来模拟耦合环境，相比于传统的交错网格建模方法或者构建连续性条件，本文提出的方法在正演精度和稳定性上凸显出很大优势，极大降低了计算内存.反演策略对多参数全波形反演至关重要，由于不同参数之间的相互耦合使得密度在多参数全波形反演中较难获得，因此本文将非均匀流-固边界耦合介质多参数全波形反演分为两个步骤完成：第一步利用变密度声波方程结合推导出的密度梯度算子进行纵波速度和密度的双参数反演；第二步根据链式法则求取横波速度的梯度，结合第一步的反演结果使用流-固边界耦合方程反演横波速度.最后通过与声波动方程数值模拟结果对比证明正演算法的准确性；上覆流体的Marmousi-2模型的数值试验测试说明反演方法的有效性和适应性. 相似文献
19.
Seismic wavefield modeling in media with fluid-filled fractures and surface topography 总被引：3，自引：0，他引：3
We present a finite difference (FD) method for the simulation of seismic wave fields in fractured medium with an irregular (non-flat) free surface which is beneficial for interpreting exploration data acquired in mountainous regions. Fractures are introduced through the Coates-Schoenberg approach into the FD scheme which leads to local anisotropic properties of the media where fractures are embedded. To implement surface topography, we take advantage of the boundary-conforming grid and map a rectangular grid onto a curved one. We use a stable and explicit second-order accurate finite difference scheme to discretize the elastic wave equations (in a curvilinear coordinate system) in a 2D heterogeneous transversely isotropic medium with a horizontal axis of symmetry (HTI). Efficiency tests performed by different numerical experiments clearly illustrate the influence of an irregular free surface on seismic wave propagation in fractured media which may be significant to mountain seismic exploration. The tests also illustrate that the scattered waves induced by the tips of the fracture are re-scattered by the features of the free surface topography. The scattered waves provoked by the topography are re-scattered by the fractures, especially Rayleigh wave scattering whose amplitudes are much larger than others and making it very difficult to identify effective information from the fractures. 相似文献