基于分数阶时间导数的黏弹性衰减VTI介质中平面波传播

目前在地震勘探频带范围内通常假设品质因子Q与频率无关,且呈衰减各向同性.事实上,相比较速度各向异性,介质的衰减各向异性同样不可忽视.本文将衰减各向异性和速度各向异性二者与常Q模型相结合,建立了黏弹性衰减VTI介质模型,并基于分数阶时间导数理论,给出了对应的本构关系和波动方程.利用均匀平面波分析和Poynting定理,推...     

Propagation of guided waves in a fluid layer bounded by two viscoelastic transversely isotropic solids

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.     

Free-field response from inclined body waves in a viscoelastic random medium

This paper deals with the free-field response of the in-plane motion resulting from a combination of inclined incident body waves. The amplification of waves in a viscoelastic layer with stochastic changes in the elastic properties and density is investigated. The method used is that of Karal and Keller and is based on the idea of the fundamental matrix. The third order correlations are neglected. The resulting integro-differential equations for the average displacements are solved by the Laplace transform. Generally, analysis indicates that the stochastic changes in the shear modulus and density enhance the damping in a significant manner. However, increases in the waves' amplification can arise in the case of a small dimensionless frequency and uncorrelated stochastic changes of material parameters.     

黏弹介质中勒夫波频散问题的统一解及其动态特征分析

目前完全弹性介质中面波频散特征的研究已较为完善,多道面波分析技术(MASW)在近地表勘探领域也取得了较好的效果,但黏弹介质中面波的频散特征研究依然较少.本文基于解析函数零点求解技术,给出了完全弹性、常Q黏弹和Kelvin-Voigt黏弹层状介质中勒夫波频散特征方程的统一求解方法.对于每个待计算频率,首先根据传递矩阵理论得到勒夫波复频散函数及其偏导的解析递推式,然后在复相速度平面上利用矩形围道积分和牛顿恒等式将勒夫波频散特征复数方程的求根问题转化为等价的连带多项式求解问题,最后通过求解该连带多项式的零点得到多模式勒夫波频散曲线与衰减系数曲线.总结了地层速度随深度递增和夹低速层条件下勒夫波频散特征根在复相速度平面上的运动规律和差异.证明了频散曲线交叉现象在复相速度平面上表现为:随频率增加,某个模式特征根的移动轨迹跨越了另一个模式特征根所在的圆,并给出了这个圆的解析表达式.研究还表明,常Q黏弹地层中的基阶模式勒夫波衰减程度随频率近似线性增加,而Kelvin-Voigt黏弹地层中的基阶模式勒夫波衰减程度随频率近似指数增加,且所有模式总体衰减程度强于常Q黏弹地层中的情况.     

Magnetoelastic shear wave propagation in pre-stressed anisotropic media under gravity

The present study investigates the propagation of shear wave (horizontally polarized) in two initially stressed heterogeneous anisotropic (magnetoelastic transversely isotropic) layers in the crust overlying a transversely isotropic gravitating semi-infinite medium. Heterogeneities in both the anisotropic layers are caused due to exponential variation (case-I) and linear variation (case-II) in the elastic constants with respect to the space variable pointing positively downwards. The dispersion relations have been established in closed form using Whittaker’s asymptotic expansion and were found to be in the well-agreement to the classical Love wave equations. The substantial effects of magnetoelastic coupling parameters, heterogeneity parameters, horizontal compressive initial stresses, Biot’s gravity parameter, and wave number on the phase velocity of shear waves have been computed and depicted by means of a graph. As a special case, dispersion equations have been deduced when the two layers and half-space are isotropic and homogeneous. The comparative study for both cases of heterogeneity of the layers has been performed and also depicted by means of graphical illustrations.     

黏弹性介质瑞雷波有限差分模拟与特性分析

本文通过数值模拟研究了介质黏弹性对瑞雷波传播的影响.模拟采用结合了交错Adams-Bashforth时间积分法、应力镜像法和多轴完美匹配层的标准交错网格高阶有限差分方案.通过模拟结果和理论结果对比,测试了方法的精度,验证了结果的正确性.在均匀半空间模型中,分别从波场快照、波形曲线及频散能量图三个角度,对黏弹性介质瑞雷波衰减和频散特性进行了详细分析.两层速度递增模型被用于进一步分析瑞雷波在黏弹性层状介质中的特性.结果表明：由于介质的黏弹性,瑞雷波振幅发生衰减,高频成分比低频成分衰减更剧烈,衰减程度随偏移距增大而增强;瑞雷波相速度发生频散,且随频率增大而增大,频散能量的分辨率有所降低;黏弹性波动方程中的参考频率,不会影响瑞雷波振幅衰减和相速度频散的程度,但决定了黏弹性和弹性介质瑞雷波相速度相等的频率位置.本研究有助于人们更好地理解地球介质中瑞雷波的行为,并为瑞雷波勘探的应用和研究提供了科学和有价值的参考.     

横向各向同性地层井孔伪瑞利波、弯曲波、螺旋波的低频极限速度

在多极源声测井中, 低频弯曲波或螺旋波被广泛应用于测量地层的横波速度. 前人的研究已证明在各向同性地层中井孔伪瑞利波、弯曲波及螺旋波的低频极限速度都等于地层横波速度. 大量的数值计算结果似乎表明此结论在横向各向同性（TI）地层情况下也能成立, 但缺乏理论证明. 本文在井孔平行于TI弹性地层对称轴的模型下, 考察了井内声源激发的波在流-固边界上的反射和透射情况, 阐述了非泄漏模式导波产生的必要条件并讨论了其速度上限值. 我们发现在各向异性参数满足一定条件的TI地层中,导波速度的低频极限值小于沿井孔方向上的横波速度. 通过对井孔导波的频散分析以及对时域全波列的数值模拟, 论证了在这类地层中进行多极源声测井时,不可能根据弯曲波或螺旋波的到达时间准确获取地层的横波速度值.     

Vibration isolation by trenches in continuously nonhomogeneous soil by the BEM

Vibration isolation of structures from ground-transmitted waves by open trenches in isotropic, linearly elastic or viscoelastic soil with a shear modulus varying continuously with depth is numerically studied. Both an exponential and a linear shear modulus variation with depth are used in this work. Waves produced by the harmonic motion of a rigid surface machine foundation are considered. The problem is solved by the frequency domain boundary element method employing the Green's function of Kausel-Peek-Hull for a thin layered half-space. Thus only the trench perimeter and the soil-foundation interface need essentially to be discretized. The proposed methodology is first tested for accuracy by solving two Rayleigh wave propagation problems in nonhomogeneous soil with known analytical solutions and/or for which experimental results are available. Then the method is applied to vibration isolation problems and the effect of the inhomogeneity on the wave screening effectiveness of trenches is studied.     

层状黏弹性介质中Rayleigh波频散曲线“交叉”现象分析

Rayleigh波勘探方法在探测近地表横波速度、动力学特征等环境与工程地球物理领域获得了广泛应用.这种方法以弹性层状介质理论为基础,然而实际介质具有黏弹性,研究面波在层状黏弹性介质中的传播特征,将为近地表面波勘探提供有益帮助.在某些弹性层状介质模型中,例如存在低速夹层和强波阻抗差异地层模型,Rayleigh波相邻两条频散曲线彼此会非常靠近,产生看似彼此"交叉"的现象,即"osculation"现象,但对于黏弹性介质中的这种现象并没有进行相关的研究.本文利用Muller法计算层状黏弹性介质Rayleigh波频散方程,基于层状介质模型中Rayleigh波频散和衰减曲线连续的性质,结合本征位移曲线特征,分析二层黏弹性介质模型中Rayleigh波频散曲线"交叉"现象以及"交叉"点附近的波动特性.结果表明:与弹性介质相比,黏弹性介质中Rayleigh波的波动特性存在明显差异,随着介质对地震波的损耗越来越强,将导致Rayleigh波频散曲线发生"交叉"现象.     

Pseudospectral modeling and dispersion analysis of Rayleigh waves in viscoelastic media

Multichannel Analysis of Surface Waves (MASW) is one of the most widely used techniques in environmental and engineering geophysics to determine shear-wave velocities and dynamic properties, which is based on the elastic layered system theory. Wave propagation in the Earth, however, has been recognized as viscoelastic and the propagation of Rayleigh waves presents substantial differences in viscoelastic media as compared with elastic media. Therefore, it is necessary to carry out numerical simulation and dispersion analysis of Rayleigh waves in viscoelastic media to better understand Rayleigh-wave behaviors in the real world. We apply a pseudospectral method to the calculation of the spatial derivatives using a Chebyshev difference operator in the vertical direction and a Fourier difference operator in the horizontal direction based on the velocity–stress elastodynamic equations and relations of linear viscoelastic solids. This approach stretches the spatial discrete grid to have a minimum grid size near the free surface so that high accuracy and resolution are achieved at the free surface, which allows an effective incorporation of the free surface boundary conditions since the Chebyshev method is nonperiodic. We first use an elastic homogeneous half-space model to demonstrate the accuracy of the pseudospectral method comparing with the analytical solution, and verify the correctness of the numerical modeling results for a viscoelastic half-space comparing the phase velocities of Rayleigh wave between the theoretical values and the dispersive image generated by high-resolution linear Radon transform. We then simulate three types of two-layer models to analyze dispersive-energy characteristics for near-surface applications. Results demonstrate that the phase velocity of Rayleigh waves in viscoelastic media is relatively higher than in elastic media and the fundamental mode increases by 10–16% when the frequency is above 10 Hz due to the velocity dispersion of P and S waves.     

A wave propagation model with the Biot and the fractional viscoelastic mechanisms

Energy loss in porous media containing fluids is typically caused by a variety of dynamic mechanisms.In the Biot theory,energy loss only includes the frictional dissipation between the solid phase and the fluid phase,resulting in underestimation of the dispersion and attenuation of the waves in the low frequency range.To develop a dynamic model that can predict the high dispersion and strong attenuation of waves at the seismic band,we introduce viscoelasticity into the Biot model and use fractional derivatives to describe the viscoelastic mechanism,and finally propose a new wave propagation model.Unlike the Biot model,the proposed model includes the intrinsic dissipation of the solid frame.We investigate the effects of the fractional order parameters on the dispersion and attenuation of the P-and S-waves using several numerical experiments.Furthermore,we use several groups of experimental data from different fluid-saturated rocks to testify the validity of the new model.The results demonstrate that the new model provides more accurate predictions of high dispersion and strong attenuation of different waves in the low frequency range.     

Dynamic stiffness and kinematic response of single piles in inhomogeneous soil

The effect of soil inhomogeneity on dynamic stiffness and kinematic response of single flexural elastic piles to vertically-propagating seismic SH waves is explored. A generalized parabolic function is employed to describe the variable shear wave propagation velocity in the inhomogeneous stratum. A layered soil with piece-wise homogeneous properties is introduced to approximate the continuous inhomogeneity in the realm of a Beam-on-Dynamic-Winkler-Foundation model. The problem is treated numerically by means of a layer transfer-matrix (Haskell–Thompson) formulation, and validated using available theoretical solutions and finite-element analyses. The role of salient model parameters such as pile-head fixity conditions, pile-to-soil stiffness ratio, surface-to-base shear wave velocity ratio and rate of inhomogeneity is elucidated. A new normalization scheme for inertial and kinematic response of such systems is presented based on an average Winkler wavenumber. With reference to long piles in moderately inhomogeneous soils, results indicate that: (a) kinematic pile response is essentially governed by a single dimensionless frequency parameter accounting for pile-to-soil stiffness ratio, pile slenderness and soil inhomogeneity and (b) definition of a characteristic pile wavelength allows an approximate estimation of pile elastodynamic response for preliminary design or analysis. Issues related to domain discretization and Winkler moduli are discussed.     

Propagation of shear waves in viscoelastic medium at irregular boundaries

The aim of the paper is to study the shear wave propagation in a viscoelastic layer over a semi-infinite viscoelastic half space due to irregularity in the viscoelastic layer. It is of great interest to study the propaga-tion of shear waves in the assumed medium having a non planar boundary due to its similarity to most of the real situations. The perturbation method is applied to find the displacement field. The effect of complex wave number on dissipation factor is analysed. Finally, as an application, the result obtained has been used to get the reflected field in viscoelastic layer when the shear wave is incident on an irregular boundary in the shape of parabolic irregularity as well as triangular notch. It is observed that the amplitude of this reflected wave decreases with increasing length of the notch, and increases with increasing depth of the irregularity.     

Propagation of seismic waves in a viscoelastic medium and the effect of topography

In this paper we have simulated the propagation of seismic waves in viscoelastic medium and calculated the effect of irregular surface topography by using finite element method. Several types of elastic and viscoelastic medium models, such as block structure, vidges and/or valley with and without a soft layer underneath have been studied. The distribution of maximum amplitude of displacement and acceleration on surface has been calculated in the case of vertically incident SH and P waves from the basement. The possible application in earthquake engineering is also discussed. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 442–449, 1991.     

Construction of constant-Q viscoelastic model with three parameters

The popularly used viscoelastic models have some shortcomings in describing relationship between quality factor (Q) and frequency, which is not consistent with the observation data. Based on the theory of viscoelasticity, a new approach to construct constant-Q viscoelastic model in given frequency band with three parameters is developed. The designed model describes the frequency-independence feature of quality factor very well, and the effect of viscoelasticity on seismic wave field can be studied relatively accurate in theory with this model. Furthermore, the number of required parameters in this model has been reduced fewer than that of other constant-Q models, this can simplify the solution of the viscoelastic problems to some extent. At last, the accuracy and application range have been analyzed through numerical tests. The effect of viscoelasticity on wave propagation has been briefly illus-trated through the change of frequency spectra and waveform in several different viscoelastic models.     

Surface waves in a micropolar elastic solid containing voids

The present paper investigates the effect of voids on the propagation of surface waves in a homogeneous micropolar elastic solid medium which contains a distribution of vacuous pores (voids). The general theory for surface wave propagation in micropolar elastic media containing voids has been presented. Particular cases of surface waves (Rayleigh’s, Love’s and Stoneley’s) in micropolar media which contain vacuous pores have been deduced from the above general theory. Discussions have been made in each case to highlight the effect of voids and micropolar character of the material medium separately. Their joint effect has also been studied in details. Modulation of Rayleigh wave velocity has been studied numerically. It is observed that Love waves are not affected by the presence of voids.     

Construction of constant-Q viscoelastic mode with three parameters

The popularly used viscoelastic models have some shortcomings in describing relationship between quality factor （Q） and frequency, which is not consistent with the observation data. Based on the theory of viscoelasticity, a new approach to construct constant-Q viscoelastic model in given frequency band with three parameters is developed. The designed model describes the frequency-independence feature of quality factor very well, and the effect of viscoelasticity on seismic wave field can be studied relatively accurate in theory with this model. Furthermore, the number of required parameters in this model has been reduced fewer than that of other constant-Q models, this can simplify the solution of the viscoelastic problems to some extent. At last, the accuracy and application range have been analyzed through numerical tests. The effect of viscoelasticity on wave propagation has been briefly illustrated through the change of frequency spectra and waveform in several different viscoelastic models.     

Love waves in a sedimentary layer lying over an anisotropic inhomogeneous half space

Summary The effects of anisotropy and inhomogeneity on the propagation of Love waves in a sedimentary layer, overlying the inhomogeneous and transversely isotropic half space, are studied in this paper. The results of numerical analysis show an appreciable variation of phase- and group-velocity of Love waves in low frequency region compared to high frequency region due to the presence of transverse isotropy and inhomogeneity in the half space. The higher values for phase velocity are found for the increasing values of anisotropy factor as well as for the greater power of density variation. However, the presence of higher anisotropy factor and inhomogeneity in the half space reduce group velocity considerably in the lower frequency region.     

适用于更广渗透率范围含泥质砂岩的广义粘弹性BISQ模型

To simultaneously take into account the Biot-flow mechanism, the squirt-flow mechanism, and the frame-viscoelasticity mechanism, a generalized viscoelastic BISQ （Biot/squirt） model is developed for wave propagation in clay-bearing sandstones based on Dvorkin＇s elastic BISQ model. The present model is extended to a wide range of permeabilities （k 〉 0.05 mD） by introducing a dimensionless correction factor for viscoelastic parameters, defined as a function of the permeability and the clay content. We describe the frame＇s stress-strain relationship of the clay-bearing sandstones by the differential constitutive equations of generalized viscoelasticity and then derive the viscoelastic-wave dynamic equations. With the assumption of a plane-wave solution, we finally yield the phase velocities and the attenuation coefficients by solving the dynamic wave equations in the frequency and wave number domain. The comparison of numerical results and experimental data shows that the generalized viscoelastic BISQ model is applicable for modeling the wave propagation in most of the sandstones mainly bearing kaolinite clay.     

含泥质低孔渗各向异性黏弹性介质中的波频散和衰减研究

本文定义了各向异性黏弹性参数修正因子,并将其引入到黏弹性模型中以体现泥质含量对黏弹性机制的影响,同时将波传播过程中孔隙介质骨架黏弹性力学机制与两种孔隙流体流动力学机制（Biot流动和喷射流动机制）有机地统一起来处理,从而给出了描述含泥质低孔渗孔隙各向异性介质中波传播规律的黏弹性Biot/squirt (BISQ)模型.数值计算结果表明,入射波的方位角、各向异性渗透率以及泥质含量等对含流体复杂孔隙介质中波频散和衰减的影响具有显著的方位各向异性特征,在低频范围内（地震波勘探频率）黏弹性力学机制对波传播能量的衰减起主导作用.