共查询到20条相似文献，搜索用时 421 毫秒
1.
本文对声波与弹性波方程进行有限元法离散,构造有限元法频散关系的一般特征值问题,分析了时间离散格式为中心差分的三角网格有限元法声波与弹性波模拟的频散特性. 比较了三种质量矩阵即分布式质量矩阵、集中质量矩阵和混合质量矩阵对有限元法频散的影响;选取四种典型三角网格,分析了混合质量矩阵有限元（MFEM）频散的方向各向异性;数值频散、方向各向异性随插值阶数的增加逐渐减弱,当空间为三阶插值时,频散主要表现为随采样率的变化而几乎无明显方向各向异性, 其频散幅值也较小. 控制其他影响因素不变的情况下,研究了不同波速比介质中弹性波的数值频散. 最后给出了三角网格MFEM的数值耗散性. 相似文献
2.
波动方程有限差分法正演模拟,对认识地震波传播规律、进行地震属性研究、地震资料地质解释、储层评价等,均具有重要的理论和实际意义.但有限差分法本身固有存在着数值频散问题,数值频散在正演模拟中是一种严重的干扰,会降低波场模拟的精度与分辨率.针对TI介质波场模拟的交错网格有限差分方法,本文从空间网格离散、时间网格离散和算子近似等三个方面对其产生的数值频散进行了分析,并结合其他学者的研究成果给出了TI介质波场模拟中压制数值频散的方法与策略：在已知介质频散关系时,对差分算子可实施算子校正;通过提高差分方程的阶数来提高波场模拟精度;采用流体力学中守恒式方程的通量校正传输方法来压制波场模拟中的数值频散;在实际正演模拟时,采用交错网格高阶有限差分方程,不仅在空间上采用高阶差分,而且在时间上也要采用高阶差分,否则只在单一方向上（空间或时间）提高方程的阶数对压制数值频散也不会取得理想的效果. 相似文献
3.
4.
分数阶微分算子具有描述历史依赖性和全域相关性的特质,本文利用这种特质描述双相介质固体骨架的黏弹性特征.基于Kjartansson常Q理论将含有分数阶时间导数的黏弹固体骨架各向异性本构关系与双相介质理论有机地结合起来,并引入流变学本构关系描述孔隙流体的黏滞性力学行为,提出一种新的基于分数阶时间导数常Q黏弹本构关系的含黏滞流体双相VTI模型.推导了相应的时间域波传播方程,然后对该方程进行了数值模拟.对整数阶导数采用高阶交错网格有限差分算法,对分数阶时间导数采用短时记忆中心差分算法,进行了不同相界、不同品质因子组及双层地质结构情况下该类介质中波场的数值模拟与特征分析.模拟结果表明：将含有分数阶时间导数的常Q黏弹固体骨架各向异性本构关系及孔隙流体的黏滞性本构关系引入双相介质理论是可行的,二者的结合能更好地反映地下介质的黏弹性特征,对于进一步认识波在黏弹各向异性孔隙介质中的传播机理具有重要意义,为反演和重构地下油气储层和结构奠定正演理论基础. 相似文献
5.
含流体孔隙介质中的波能量耗散通常由多种力学机制造成.传统Biot理论中的能量耗散仅仅考虑了固流两相相对运动引起的摩擦耗散,无法准确预测波在孔隙介质中低频段出现的高频散与强衰减现象.为了建立一个能准确预测地震波频段高频散与强衰减现象的动力学模型,我们在Biot理论的基础上引入黏弹性机制,并利用分数阶导数刻画黏弹性本构关系,最终获得了一种新的孔隙介质波传播模型.与传统的Biot模型相比,新模型考虑了含流体孔隙介质中固体骨架的内耗散,对波能量耗散的刻画更为精准.通过数值算例,我们研究了分数阶导数的阶数参数对快P波和S波频散和衰减的影响,并通过来自不同地区且具有不同物理性质的几组流体饱和岩芯实验数据,对比研究了新模型的有效性.结果表明,文章提出的新模型能更准确地预测快P波和S波在低频段出现的高频散和强衰减现象. 相似文献
6.
间断有限元（Discontinuous Galerkin：DG）方法具有低数值频散、网格剖分灵活、能模拟地震波在复杂介质中传播等优点.因此,本文将一种新的DG方法推广到双相和黏弹性等复杂介质的地震波场模拟,发展了求解Biot弹性波方程和D'Alembert介质波动方程的DG方法.首先通过引入辅助变量将Biot双相介质弹性波方程和D'Alembert介质波动方程转化为关于时间-空间的一阶偏微分方程组,然后对该方程组进行DG空间离散,得到半离散化的常微分方程组.最后,对此常微分方程组,应用加权的Runge-Kutta格式进行时间推进计算.数值结果表明,DG方法可以有效地求解Biot双相介质弹性波方程和D'Alembert介质波动方程,并能很好地压制因离散求解波动方程而产生的数值频散,获得清晰的各种地震波震相. 相似文献
7.
为了探明介质的非完全弹性对瑞利型槽波传播的影响以及利用瑞利型槽波品质因子Q_R反演煤层碳化程度、裂隙、煤层厚度等煤层属性,在本文中,我们将复速度引入瑞利型槽波频散方程,计算了煤层瑞利型槽波的品质因子Q_R曲线;采用控制变量法,分析了瑞利型槽波品质因子Q_R随煤层厚度、煤层纵横波品质因子以及围岩纵横波品质因子的变化规律。研究表明:瑞利型槽波品质因子Q_R曲线与群速度曲线变化趋势一致;埃里相(Airy-phase)频率附近瑞利型槽波的品质因子Q_R最小,且埃里相频率随煤层厚度增加而减小;瑞利型槽波品质因子Q_R随煤层横波品质因子Q_(S2)的增大而增大。利用瑞利型槽波品质因子Q_R曲线可以实现对瑞利型槽波的吸收衰减进行补偿、对煤层厚度进行预测以及岩性成像。 相似文献
8.
为提高频率域弹性波动方程数值求解的计算效率,本文引入近似解析离散化(NAD)方法将其进行数值离散并得到大型线性代数方程组.在详细分析了相应系数矩阵的稀疏分块结构与数学性质之后,本文提出采用不精确旋转分块三角预处理子加速Krylov子空间迭代方法来快速求解该线性方程组,并利用数值试验证实这种方法在弹性波场模拟方面的数值效率.通过与另外两种经典数值方法(常规有限差分方法和交错网格有限差分方法)对多种介质模型进行波场模拟、数值频散分析以及与解析解的波形对比,NAD方法显示了其在压制数值频散和提高计算效率方面的优势以及对复杂介质模型弹性波场数值模拟的有效性. 相似文献
9.
《地球物理学进展》2017,(4)
TTI(Tilted Transversely Isotropic)各向异性是对地下岩石中广泛存在的规则发育的裂缝和层理的一种有效的弹性近似,基于TTI介质的地震波数值模拟技术是分析地震波在复杂各向异性介质中的传播机理的有效工具.同时,高精度的数值模拟算法也能为后续的逆时偏移技术提供重要的技术支撑.由于TTI介质中地震波方程的弹性参数众多且变化复杂,常规有限差分技术在解决TTI介质正演模拟问题时往往会产生严重的数值频散现象,降低了数值模拟精度.通量校正传输(FluxCorrected Transport,FCT)技术能够有效地压制由空间离散产生的数值频散.本文将FCT技术用于TTI介质中弹性波方程的交错网格高阶精度差分正演,在数值模拟过程中通过对波场进行漫射和反漫射校正实现了空间网格频散的压制.模型模拟结果表明,与常规有限差分算法相比,本文算法能够有效的压制大网格条件下的数值频散,提高模拟精度. 相似文献
10.
本文通过数值模拟研究了介质黏弹性对瑞雷波传播的影响.模拟采用结合了交错Adams-Bashforth时间积分法、应力镜像法和多轴完美匹配层的标准交错网格高阶有限差分方案.通过模拟结果和理论结果对比,测试了方法的精度,验证了结果的正确性.在均匀半空间模型中,分别从波场快照、波形曲线及频散能量图三个角度,对黏弹性介质瑞雷波衰减和频散特性进行了详细分析.两层速度递增模型被用于进一步分析瑞雷波在黏弹性层状介质中的特性.结果表明：由于介质的黏弹性,瑞雷波振幅发生衰减,高频成分比低频成分衰减更剧烈,衰减程度随偏移距增大而增强;瑞雷波相速度发生频散,且随频率增大而增大,频散能量的分辨率有所降低;黏弹性波动方程中的参考频率,不会影响瑞雷波振幅衰减和相速度频散的程度,但决定了黏弹性和弹性介质瑞雷波相速度相等的频率位置.本研究有助于人们更好地理解地球介质中瑞雷波的行为,并为瑞雷波勘探的应用和研究提供了科学和有价值的参考. 相似文献
11.
目前在地震勘探频带范围内通常假设品质因子Q与频率无关,且呈衰减各向同性.事实上,相比较速度各向异性,介质的衰减各向异性同样不可忽视.本文将衰减各向异性和速度各向异性二者与常Q模型相结合,建立了黏弹性衰减VTI介质模型,并基于分数阶时间导数理论,给出了对应的本构关系和波动方程.利用均匀平面波分析和Poynting定理,推导出准压缩波qP、准剪切波qSV和纯剪切波SH的复速度、相速度、能量速度以及品质因子的解析表达式.对模型的正确性进行了数值验证,并分析了qP,qSV和SH波在介质中的传播特性.数值试验结果表明:本模型能够实现理想的恒定Q行为,表现了品质因子和速度的各向异性特征,显示出黏弹性增强将导致能量速度和相速度的频散曲线变化剧烈;速度和衰减各向异性参数与传播角度之间的耦合效应对qP,qSV和SH波的速度和能量影响明显;qP,qSV和SH波的频散曲线和波前面随着衰减各向异性强度的改变发生显著变化,其中耦合在一起的qP和qSV波变化趋势相同,而SH波与它们呈现相反的变化规律.本研究为从常Q模型角度分析地震波在衰减各向异性黏弹性介质中的传播特征奠定了理论基础. 相似文献
12.
13.
The simplified macro‐equations of porous elastic media are presented based on Hickey's theory upon ignoring effects of thermomechanical coupling and fluctuations of porosity and density induced by passing waves. The macro‐equations with definite physical parameters predict two types of compressional waves (P wave) and two types of shear waves (S wave). The first types of P and S waves, similar to the fast P wave and S wave in Biot's theory, propagate with fast velocity and have relatively weak dispersion and attenuation, while the second types of waves behave as diffusive modes due to their distinct dispersion and strong attenuation. The second S wave resulting from the bulk and shear viscous loss within pore fluid is slower than the second P wave but with strong attenuation at lower frequencies. Based on the simplified porous elastic equations, the effects of petrophysical parameters (permeability, porosity, coupling density and fluid viscosity) on the velocity dispersion and attenuation of P and S waves are studied in brine‐saturated sandstone compared with the results of Biot's theory. The results show that the dispersion and attenuation of P waves in simplified theory are stronger than those of Biot's theory and appear at slightly lower frequencies because of the existence of bulk and shear viscous loss within pore fluid. The properties of the first S wave are almost consistent with the S wave in Biot's theory, while the second S wave not included in Biot's theory even dies off around its source due to its extremely strong attenuation. The permeability and porosity have an obvious impact on the velocity dispersion and attenuation of both P and S waves. Higher permeabilities make the peaks of attenuation shift towards lower frequencies. Higher porosities correspond to higher dispersion and attenuation. Moreover, the inertial coupling between fluid and solid induces weak velocity dispersion and attenuation of both P and S waves at higher frequencies, whereas the fluid viscosity dominates the dispersion and attenuation in a macroscopic porous medium. Besides, the heavy oil sand is used to investigate the influence of high viscous fluid on the dispersion and attenuation of both P and S waves. The dispersion and attenuation in heavy oil sand are stronger than those in brine‐saturated sandstone due to the considerable shear viscosity of heavy oil. Seismic properties are strongly influenced by the fluid viscosity; thus, viscosity should be included in fluid properties to explain solid–fluid combination behaviour properly. 相似文献
14.
深层-超深层油气地震勘探涉及高温介质地震波传播问题,热弹介质参数对地震波传播有重要影响.含弛豫时间修正项的Lord-Shulman双曲型耦合热弹波动方程从理论上预测了热弹性介质中存在快纵波、慢纵波(一种准静态慢纵波,简称热波)和横波的传播,两个纵波为热耗散衰减波而横波不受介质热特性的影响.本文结合平面波频散分析和格林函数法数值模拟,详细研究两个热耗散衰减波的频散和衰减特征,着重分析热导率、热膨胀系数及比热的变化对波速和衰减的影响.研究表明热导率作为主要参数决定了波速与衰减的临界变化,热膨胀系数对波速和衰减的幅度有明显影响,比热则兼顾了前两个热弹系数的影响特征.最后,利用热弹性动力学频率域的二阶格林函数进行波场快照数值模拟,展示热弹性介质中纵波、横波和热波的传播行为. 相似文献
15.
本文首先由Christoffel方程推导出黏弹性EDA介质中均匀、 非均匀P波、 SV波和SH波的相速度表达式, 然后参照极端各向异性介质的相关计算方法, 推导出EDA介质中均匀、 非均匀地震波相衰减系数和群衰减系数的表达式, 并通过数值计算分析了相速度、 相衰减系数、 群衰减系数与裂隙方位的关系. 结果表明: 均匀介质中SH波的相速度和相衰减系数均可指示裂隙的走向; 非均匀介质中SH波相衰减系数随非均匀角的增大而增大, 且其对称轴与介质对称轴的夹角也相应增加; 由于地震波振幅的衰减随岩石物理性质的变化比地震波速度的变化更为灵敏, 而且携带了更多的岩石物理性质信息, 因此可用来探明裂隙走向、 密度及含水特性, 进而应用于预测、 预防地下工程地质灾害事故. 相似文献
16.
Hua-You Chai Kok-Kwang Phoon Chang-Fu Wei Ying-Fa Lu 《Journal of Applied Geophysics》2011,73(1):49-58
In the wave field induced by active sources, the observed phase velocity of surface waves is influenced by both mode incompatibility (i.e. non-planar spread of surface waves is idealized as plane waves) and body waves. Effects of sources are usually investigated based on numerical simulations and physical models. Several methods have been proposed to mitigate the effects. In application, however, these methods may also have difficulties since the energy of the body waves depends on soil stratification and parameters. There are multiple modes of surface waves in layered media, among which the higher modes dominate the wave field for soils with the irregular shear velocity profiles. Considering the mode incompatibility and the higher modes, we derive analytical expressions for the effective phase velocity of the surface waves based on the thin layer stiffness method, and investigate the effects of the body waves on the observed phase velocity through the phase analysis of the vibrations of both the surface waves and the body waves. The results indicate that the effective phase velocity of the surface waves in layered media varies with the frequency and the spread distance, and is underestimated compared to that of the plane surface waves in the spread range less than about one wavelength. The oscillations that appeared in the observed phase velocity are due to the involvement of the body waves. The mode incompatibility can be ignored in the range beyond one wavelength, while the influence range of the body waves is far beyond one wavelength. The body waves have a significant influence on the observed phase velocity of the surface waves in soils with a soft layer trapped between the first and the second layers because of strong reflections. 相似文献
17.
随机弹性介质中地震波散射衰减分析(英文) 总被引：2，自引：0，他引：2
地震波衰减一直是许多学科研究的热点,因为可以反映介质的特性。导致地震波衰减的因素很多,如:传播过程中由于能量扩散导致的几何衰减,固体岩石内部晶粒间相对滑移导致的摩擦衰减,岩石结构不均匀引起的地震波散射衰减。本文主要从统计的观点出发,通过多次数值模拟的方法研究纵波散射在随机弹性介质中所引发的衰减。首先用随机理论建立了二维空间随机弹性介质模型,然后用错格伪谱法的数值方法模拟了波在随机介质中的传播,再通过波场中虚拟检波器的记录,用谱比法估计了弹性波在随机介质中的散射衰减。不同非均匀程度随机弹性介质中的数值结果表明:介质不均匀程度越高,散射衰减越大;在散射体尺寸小于波长的前提下,不同散射体尺寸的计算结果说明:散射体尺寸越大,弹性波衰减越明显。最后提出了一种不均匀孔隙介质中流体流动衰减的方法。通过对随机孔隙介质中地震波的总衰减和散射衰减分别进行了计算,并定量得出了随机孔隙介质中流体流动衰减,结果表明:在实际地震频段下,当介质不均匀尺度101米量级时,散射衰减比流体流动衰减要大,散射衰减是地震波在实际不均匀岩石孔隙介质中衰减的主要原因。 相似文献
18.
弹性波在储层渗流场中的传播与衰减规律是研究波场强化采油动力学机理的重要基础.基于等效流体理论和饱和静态流体弹性波传播Biot理论,建立油水两非混相流体渗流条件下储层多孔介质中弹性波传播的动力学模型,通过算例求解与分析,发现含油水两相渗流储层多孔介质中同时存在着3种纵波P1、P2、P3和1种横波S;受频率和含水饱和度的影响,各波波速和品质因子呈现出不同变化规律,4种体波波速与频率、饱和度正相关,P1、P2波品质因子与饱和度正相关,P3和S波品质因子与饱和度负相关;最后,通过与传统静态弹性波模型结果对比,进一步分析了宏观渗流场对弹性波传播特征的影响规律,为揭示低频人工地震波辅助强化采油技术的动力学机理和工艺参数优化提供了重要理论依据. 相似文献
19.
Jixin Deng Shangxu Wang Gengyang Tang Jianguo Zhao Xiangyang Li 《Studia Geophysica et Geodaetica》2013,57(3):482-506
In sedimentary rocks attenuation/dispersion is dominated by fluid-rock interactions. Wave-induced fluid flow in the pores causes energy loss through several mechanisms, and as a result attenuation is strongly frequency dependent. However, the fluid motion process governing the frequency dependent attenuation and velocity remains unclear. We propose a new approach to obtain the analytical expressions of pore pressure, relative fluxes distribution and frame displacement within the double-layer porous media based on quasi-static poroelastic theory. The dispersion equation for a P-wave propagating in a porous medium permeated by aligned fractures is given by considering fractures as thin and highly compliant layers. The influence of mesoscopic fluid flow on phase velocity dispersion and attenuation is discussed under the condition of varying fracture weakness. In this model conversion of the compression wave energy into Biot slow wave diffusion at the facture surface can result in apparent attenuation and dispersion within the usual seismic frequency band. The magnitude of velocity dispersion and attenuation of P-wave increases with increasing fracture weakness, and the relaxation peak and maximum attenuation shift towards lower frequency. Because of its periodic structure, the fractured porous media can be considered as a phononic crystal with several pass and stop bands in the high frequency band. Therefore, the velocity and attenuation of the P-wave show an oscillatory behavior with increasing frequency when resonance occurs. The evolutions of the pore pressure and the relative fluxes as a function of frequency are presented, giving more physical insight into the behavior of P-wave velocity dispersion and the attenuation of fractured porous medium due to the wave-induced mesoscopic flow. We show that the specific behavior of attenuation as function of frequency is mainly controlled by the energy dissipated per wave cycle in the background layer. 相似文献
20.
The analysis of Stoneley wave propagation in a fracture is essential for the identification and evaluation of fracture parameters from the borehole Stoneley wave. Also, it is important for many geophysics considerations, e.g. for tremor and long-period events observed in volcanoes and geothermal areas. In this paper, we investigate the guided waves propagation in a fluid layer lying between two viscoelastic vertically transversely isotropic media. The viscoelastic mechanism models the attenuation due to the presence of fluid saturation in the rock. A model based on the superposition of three inhomogeneous partial plane waves: one in the fluid and two heterogeneous waves in the solid is developed. The dispersion equation is obtained for this case. A numerical solution is carried out to obtain the guided wave velocity and attenuation coefficient. The results of this investigation show that there is a strong correlation between the velocity dispersion and attenuation of Stoneley wave and the anisotropic parameters of the medium especially in a sandstone (fast) medium. 相似文献