共查询到19条相似文献,搜索用时 250 毫秒
1.
Rayleigh面波的频散特性可以用来研究地表浅层结构. 本文使用时域有限差分法来模拟复杂黏弹性介质中的Rayleigh面波,研究了Q值对面波频散特性的影响.文中采用旋转交错网格有限差分,以非分裂卷积形式的完全匹配层为吸收边界,推出了求解二阶位移-应力各向同性黏弹性波动方程的数值方法.为了检验数值解的精度,首先将简单模型的正演结果与解析解对比,验证了方法的正确性;然后模拟了横向缓变层状介质和含有洞穴的介质中的面波,对弹性和黏弹性介质中的面波的频散特性进行对比分析.模拟结果表明浅层Q值对面波的频散特性有显著的影响;强吸收情况下,高阶面波的能量相对低阶面波能量显著增强. 相似文献
2.
3.
地震波在非弹性介质中的衰减效应常用品质因子Q度量.相对准确的Q模型对提高强衰减介质中地震波成像的质量至关重要.本文提出了黏声介质反射波形反演(QRWI)方法来重建地下宏观Q模型.在缺乏大偏移距和低频地震数据时,该方法以黏声波方程为波场传播引擎,利用反射波核函数对模型中深部的敏感性去提取背景Q值.当速度高、低波数成分均已知时,基于波形拟合的QRWI可以获得较高分辨率的反演结果.由于地下介质速度的高波数扰动很难准确估计,本文通过引入峰值频移目标函数,极大地降低了QRWI对速度高波数成分的依赖.理论合成数据实验结果表明,本文方法反演得到的宏观Q模型可以满足衰减补偿逆时偏移成像的要求. 相似文献
4.
本文基于Kjartansson常Q模型理论,推导了常Q衰减介质中黏声波和黏弹性波的速度-应力方程,并采用基于二项式窗函数的优化交错网格有限差分方法进行了数值模拟,同时引入不分裂的复频移卷积完全匹配层(CPML)吸收边界条件,以消除边界反射.使用基于自适应时间步长记忆方法的中心差分近似时间分数阶导数,与常用的短时记忆方法相比,提高了波动方程的离散化精度和计算效率.通过对比均匀模型下声波的数值解与解析解,验证了算法的精确性,并进一步分析了不同品质因子下地震波的频散及衰减特征.对BP盐丘模型的数值模拟结果可以较好地反映本文数值方法对复杂介质的适应性及频散压制效果. 相似文献
5.
三维声波方程相比二维声波方程能够更好的模拟三维空间的地震波传播,模拟标量近似下的弹性波在三维复杂介质的传播过程.基于非规则网格的正演模拟方法的格子法可以处理很好的刻画起伏地表、速度间断面等复杂构造,但是这类方法需要大量的几何描述来描述网格.本文提出了三维六面体双重网格的格子法来模拟声波方程,一方面该方法继承了格子法能够灵活处理自由表面和速度间断面的特性.另一方面,该方法通过双重网格的实现极大的减少了几何描述文件的大小,可以最大的实现GPU加速,实现粗粒度并行,在节省了几何描述空间的同时达到了很高的加速比. 相似文献
6.
瑞雷面波经常被用来反演地表浅层横波速度,受到越来越广泛的关注。对瑞雷波的研究一般都基于完全弹性介质,而实际地层更接近黏弹性介质,对黏弹性介质中的瑞雷面波进行模拟更具实际意义。本文采用广义标准线性体模型来描述黏弹性介质,并采用交错网格有限差分法对考虑水平自由表面的黏弹介质进行正演模拟,再与弹性介质中的结果进行对比分析。首先采用非线性最优化算法根据期望常数品质因子直接求取松弛时间来拟合常Q模型,并给出广义标准线性固体的具体算例,实施自由表面条件时采用声学-弹性边界近似法,通过剪切模量不变来考虑自由表面上、下横向应力保持连续的条件。对于非自由表面,采用非分裂的多轴卷积完全匹配层来吸收波场。然后对几种典型的数值模型进行正演模拟计算,数值解与解析解的对比验证了本文方法的准确性与有效性,正演结果的对比表明波场尤其是面波频散会受黏弹性影响,因此有必要在面波勘探中考虑黏弹性因素。 相似文献
7.
吸收衰减是地震波在实际地球介质中传播的固有特征.在实际应用中,通常假设表征吸收衰减特征的品质因子Q在地震频带范围内不随频率变化.高阶广义流变模型能够在时间域内精确的表征品质因子Q不随频率变化的特征,为黏弹性介质波动方程精细模拟奠定了基础.基于广义标准线性体模型理论,采用最小二乘拟合方法对Q值不随频率变化特征进行拟合,分析了不同阶次广义标准线性体模型对黏弹性介质Q值特征的拟合程度,在权衡计算精度和三维计算量的基础上,确定了五阶广义标准线性体模型并建立了相应的三维黏弹性波的速度-应力方程,结合CFS-PML边界条件开展了高精度三维黏弹性波正演模拟.通过均匀介质正演模拟,验证算法的正确性,明确了地震波的传播时的吸收衰减特征,对三维盐丘模型进行数值模拟,表明了五阶广义标准线性体可以精确的模拟黏弹性介质地震波吸收衰减特征. 相似文献
8.
《地球物理学进展》2016,(2)
正演模拟是瑞雷面波研究的一个主要方面,目前对于瑞雷面波的正演研究大多是基于均匀层状弹性介质条件下的瑞雷面波频散曲线方程,这只适用于层状模型,无法模拟全波场,而且不考虑实际介质的黏弹效应.本文采用交错网格高阶有限差分法对黏弹性介质中的瑞雷面波进行了高精度全波场模拟,并对频散特征进行了提取与分析.其中采用属于非线性最优化的Levenberg-Marquarat方法直接计算松弛时间来拟合常Q模型,并将应力镜像法与紧致差分格式相结合来准确实施自由表面条件,在其余边界处以非分裂的多轴卷积完全匹配层为吸收边界.然后利用相移法从地震记录中提取频散剖面并对几种典型模型的面波频散特征进行了对比分析.结果表明黏弹性对面波的频散特性有显著影响,面波勘探中有必要考虑黏弹性因素. 相似文献
9.
为更好地研究复杂构造和速度分布条件下地震资料的精确偏移.在波场模拟方面,本文采用了高阶交错网格有限差分[1]推导了二维各向同性介质声波方程的数值模拟公式.并给出了声波方程完全匹配层吸收边界[2-5](PML).在偏移方面,本文介绍了一种能适应介质速度横向变化的非稳态相移算子(NSPS)[6]及其叠后深度偏移方法.并对简单的四层模型和Marmousi模型进行试算.偏移结果表明,该方法可以有效地处理复杂条件下的各向同性介质声波偏移问题. 相似文献
10.
黏声波方程常被用于描述地下介质的黏弹性及波的传播现象,频域有限差分(finite difference frequency domain, FDFD)方法是黏声波和黏弹性波波场模拟的常用工具.目前FDFD黏声波模拟常用的二阶五点方法和优化九点方法在一个波长内的网格点数小于4时误差较大.通过令FDFD系数随一个波长内的网格点数自适应从而提高FDFD方法的精度,本文针对黏声波波场模拟发展了一种适用于不同空间采样间隔之比的通用格式自适应系数FDFD方法.同时,为了验证自适应系数FDFD方法对一般黏声波模型的有效性,本文针对三个典型的黏声波模型,分别采用解析解和基于高阶FDFD的参考解验证了所提出方法的有效性.本方法的FDFD格式通过在传统的二阶FDFD格式的基础上引入相关校正项得到,其中校正项按网格点与中心点的距离进行分类选取,同时校正项对应的自适应FDFD系数不仅和空间采样间隔之比相关,还和一个波长内的采样点数相关.所需的自适应FDFD系数可通过声波方程的数值频散关系和查找表高效给出.数值频散分析表明,在空间采样间隔相等或不等的情况下,以相速度误差不超过1%为标准,通用格式自适应系数FDF... 相似文献
11.
分数阶微分算子具有描述历史依赖性和全域相关性的特质,本文利用这种特质描述双相介质固体骨架的黏弹性特征.基于Kjartansson常Q理论将含有分数阶时间导数的黏弹固体骨架各向异性本构关系与双相介质理论有机地结合起来,并引入流变学本构关系描述孔隙流体的黏滞性力学行为,提出一种新的基于分数阶时间导数常Q黏弹本构关系的含黏滞流体双相VTI模型.推导了相应的时间域波传播方程,然后对该方程进行了数值模拟.对整数阶导数采用高阶交错网格有限差分算法,对分数阶时间导数采用短时记忆中心差分算法,进行了不同相界、不同品质因子组及双层地质结构情况下该类介质中波场的数值模拟与特征分析.模拟结果表明:将含有分数阶时间导数的常Q黏弹固体骨架各向异性本构关系及孔隙流体的黏滞性本构关系引入双相介质理论是可行的,二者的结合能更好地反映地下介质的黏弹性特征,对于进一步认识波在黏弹各向异性孔隙介质中的传播机理具有重要意义,为反演和重构地下油气储层和结构奠定正演理论基础. 相似文献
12.
Seismic wavefield modeling is important for improving seismic data processing and interpretation. Calculations of wavefield propagation are sometimes not stable when forward modeling of seismic wave uses large time steps for long times. Based on the Hamiltonian expression of the acoustic wave equation, we propose a structure-preserving method for seismic wavefield modeling by applying the symplectic finite-difference method on time grids and the Fourier finite-difference method on space grids to solve the acoustic wave equation. The proposed method is called the symplectic Fourier finite-difference (symplectic FFD) method, and offers high computational accuracy and improves the computational stability. Using acoustic approximation, we extend the method to anisotropic media. We discuss the calculations in the symplectic FFD method for seismic wavefield modeling of isotropic and anisotropic media, and use the BP salt model and BP TTI model to test the proposed method. The numerical examples suggest that the proposed method can be used in seismic modeling of strongly variable velocities, offering high computational accuracy and low numerical dispersion. The symplectic FFD method overcomes the residual qSV wave of seismic modeling in anisotropic media and maintains the stability of the wavefield propagation for large time steps. 相似文献
13.
14.
15.
16.
全波形反演是地震资料处理中速度建模的有力工具,相比层析成像等速度建模方法它能够得到速度场的更高频成分.本文给出了基于声波方程格子法正演的时间域全波形反演方法,该方法用非规则、非结构化的三角网格来离散计算区域及模型参数,能实现网格粒度与反演分辨率在空间上的自动匹配,内存需求少,计算效率高;采用L-BFGS优化方法,以分频段变网格的方式实施多尺度反演.以二维Overthrust模型进行了速度反演数值测试,显示了该方法的高效性和潜力. 相似文献
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.
本文定义了各向异性黏弹性参数修正因子,并将其引入到黏弹性模型中以体现泥质含量对黏弹性机制的影响,同时将波传播过程中孔隙介质骨架黏弹性力学机制与两种孔隙流体流动力学机制(Biot流动和喷射流动机制)有机地统一起来处理,从而给出了描述含泥质低孔渗孔隙各向异性介质中波传播规律的黏弹性Biot/squirt (BISQ)模型.数值计算结果表明,入射波的方位角、各向异性渗透率以及泥质含量等对含流体复杂孔隙介质中波频散和衰减的影响具有显著的方位各向异性特征,在低频范围内(地震波勘探频率)黏弹性力学机制对波传播能量的衰减起主导作用. 相似文献