岩体夹层应力波能量演化及应力响应特征分析

夹层是常见的地质结构,在地震或爆破荷载作用下,夹层对应力波的传播以及对岩体的响应具有重要的影响。以往对应力波在夹层介质中传播的研究多集中于夹层对应力波的隔振或透射性能,而对应力波在夹层中的多次折、反射过程中能量的演化规律缺乏讨论,对夹层介质的应力响应与破坏没有开展较好地分析。通过理论方法对应力波在夹层内部传播过程中的能量系数的变化规律进行了研究,分析了岩体介质波阻抗和应力波入射角对夹层内外介质中累积波动能量系数的影响规律,以及平面型边坡中软弱夹层的应力响应特征和动态稳定性。研究发现,应力波在夹层内部往复反射过程中,夹层内剩余应力波能量随折、反射发生次数呈指数曲线下降,第4次折、反射后产生的应力波能量可以忽略;夹层内外介质中应力波的累积能量系数的差异随着夹层内外介质波阻抗的相对差异的增大而增大。在平面谐波入射下,边坡内部的夹层中的剪应力和抗剪强度呈波动变化;相对P波,SV波入射会产生较高水平的剪应力,对边坡稳定性影响最大。且SV波入射时,边坡的安全系数对夹层的倾角变化更为敏感,随倾角增大而迅速降低。     

不同饱和度土层分界面上剪切波的反射与透射

基于多相孔隙介质弹性理论,给出了非饱和土中不同弹性波的传播方程。根据分界面上的边界条件,建立了各势函数波幅值之间的关系式,讨论了入射剪切波在不同饱和度土层分界面上的反射与透射问题。在无限空间非饱和土体中存在3种压缩波和1种剪切波,因此,当剪切波传播到不同饱和度的非饱和土层分界面上将分别在上、下土层激发产生4种反射波和4种透射波。推导出不同反射波和透射波的振幅比例系数和能量比例系数的理论表达式,并且在此基础上进行数值分析。在数值算例中分别研究了各反射波与透射波的能量比例系数（即能量反射率和能量透射率）受入射频率、入射角度以及上、下土层土体饱和度变化的影响情况。计算结果表明：各能量反射率和能量透射率不仅与入射角和入射频率有关,而且其受上、下土层饱和度变化的影响也同样不能忽视。     

斜入射线弹性节理应力波传播特征的数值模拟

依据弹性波理论,应力波斜入射线弹性节理时会发生波场分解。根据入射、透射及反射各波形的不同到时,运用离散元软件UDEC模拟应力波在含倾斜节理岩体中的传播并计算其透射、反射系数,并分析其波型转换规律。应力波斜入射单节理时,模拟得到的透射、反射系数随节理刚度、入射角度的变化规律,与已有的理论解是吻合的。应力波斜入射一组平行节理时,随着节理间距的增大,其同类波的透射系数Tpp、Tss先增大后减小,最后趋于稳定值;节理条数越多,Tpp、Tss越小。此外,不同条数的节理,透射系数达到最大值的临界节理间距值基本一致,但趋于稳定时的节理间距值随节理条数的增加而逐渐增大。     

粘弹性介质地震波传播特征及反射特征研究

理想弹性介质并不能解决许多复杂的实际问题。在实际介质中,地震波的传播过程存在吸收作用,使地震波能量发生了损耗、频带变窄、相位延迟、频率降低,尤其是高频部份。这种特性可以用传播矢量和衰减矢量来表示,衰减矢量越大,介质的吸收越强。这里通过引入广义平面波、纵横波品质因子以及衰减角,利用应力应变关系和连续性条件,推导了粘弹性介质精确的Zoeppritz方程,并通过求解精确的Zoeppritz方程,分析了粘弹性介质中传播矢量、衰减矢量、纵波反射系数、转换波反射系数与频率、衰减角、入射角、品质因子、速度等物理量之间的关系,详细地了解了粘弹性介质中地震波传播特征及反射特征,为粘弹性介质中储层预测问题的研究提供了理论基础。     

Effect of initial stress on reflection at the free surface of anisotropic elastic medium

The propagation of plane waves is considered in a general anisotropic elastic medium in the presence of initial stress. The Christoffel equations are solved into a polynomial of degree six. The roots of this polynomial represent the vertical slowness values for the six quasi-waves resulting from the presence of a discontinuity in the medium. Three of these six values are identified with the three quasi-waves traveling in the medium but away from its boundary. Reflection at the free plane surface is studied for partition of energy among the three reflected waves. For post-critical incidence, the reflected waves are inhomogeneous (evanescent) waves. Numerical examples are considered to exhibit the effects of initial stress on the phase direction, attenuation and reflection coefficients of the reflected waves. The phase velocities and energy shares of the reflected waves change significantly with initial stress as well as anisotropic symmetry. The presence of initial stress, however, has a negligible effect on the phase directions of reflected waves.     

Propagation of Love wave in fiber‐reinforced medium lying over an initially stressed orthotropic half‐space

In this paper, the propagation of Love waves in a fiber‐reinforced layered medium lying over an elastic orthotropic half‐space under initial stress has been investigated. We have obtained the velocity equations for Love wave in this media. It is observed that propagation of Love wave is influenced by reinforced parameters and initial stressed parameter. The velocity of Love wave has been computed for three different cases. Our computed equation of Love wave coincides with the standard equation of Love wave for the case of homogeneous layer and homogeneous half‐space (AEH Love, 1911). To study the effect of reinforced and initial stressed parameters, we have computed the numerical values for phase velocity and plotted in several figures. It is observed that the phase velocity decreases with the increases of reinforced parameters and initial stressed parameter. Using MATLAB software, GUI has been developed to generalize the effect of various parameters. Copyright © 2014 John Wiley & Sons, Ltd.     

Analysis of Blast Wave Interaction with a Rock Joint

The interaction between rock joints and blast waves is crucial in rock engineering when rock mass is suffered from artificial or accidental explosions, bursts or weapon attacks. Based on the conservation of momentum at the wave fronts and the displacement discontinuity method, quantitative analysis for the interaction between obliquely incident P- or S-blast wave and a linear elastic rock joint is carried out in the present study, so as to deduce a wave propagation equation. For some special cases, such as normal or tangential incidence, rigid or weak joint, the analytical solution of the stress wave interaction with a rock joint is obtained by simplifying the wave propagation equation. By verification, it is found that the transmission and reflection coefficients from the wave propagation equation agree very well with the existing results. Parametric studies are then conducted to evaluate the effects of the joint stiffness and incident waves on wave transmission and reflection. The wave propagation equation derived in the present study can be straightforwardly extended for different incident waveforms and nonlinear rock joints to calculate the transmitted and reflected waves without mathematical methods such as the Fourier and inverse Fourier transforms.     

Propagation Characteristics of Blasting Stress Waves in Layered and Jointed Rock Caverns

The existence of joint fissures makes explosive actions between rock masses more complex. Therefore, it is of great significance to carry out experiments studying blasting stress waves propagating in jointed rock masses. Based on the Froude Similarity principle, the geological mechanical models of intact rock masses and jointed rocks have been proposed. A blasting vibration experiment was carried out and demonstrated that the propagation of the blasting stress waves and changing structures have an important relationship. A numerical simulation of the blasting stress wave propagation law in a layered jointed rock mass was carried out. This study found that with an increase in the joint angle, the peak velocity of blasting stress wave, transmission coefficient and reflection coefficient all gradually increased, while the attenuation coefficient gradually decreased. With an increase in joint spacing, the attenuation rate of the blasting stress waves increased.     

The study of reflection/transmission phenomena in a corrugated interface between two magnetoelastic transversely isotropic media

In this paper, an integrated approach has been carried out with an intent to study the reflection and transmission phenomena of plane SH-type wave on a corrugated interface separating two magnetoelastic transversely isotropic half-space. Closed form expressions for reflection and transmission coefficients have been derived for both plane and corrugated surface. Rayleigh’s method of approximation have been incorporated to deduce equations for the first- and second-order approximation of corrugation. Analytical solutions for both the half-spaces have been worked out. All possible cases have been apprehended specially for anisotropic and an isotropic medium. The effects of magnetoelastic coupling parameter, angle at which the wave crosses the magnetic field, corrugation amplitude, frequency factor and wave number have been explained by collaborating with graphical analysis.     

Numerical study on S-wave transmission across a rough,filled discontinuity

This paper presents a numerical simulation of S-wave propagation across a rough, filled discontinuity using the universal distinct element code (UDEC). The ability of UDEC to simulate a stress wave across a smooth and planar discontinuity filled with an elastic material is validated through comparisons with analytical solutions. Next, the effect of the plastic deformation of the fill on the wave propagation is investigated. The model is extended to further study S-wave propagation across a filled discontinuity with rough interfaces, which is described using a sawtooth. The transmission coefficient defined by the energy is used to measure the wave attenuation. Finally, a parametric study is conducted to investigate the influences of the filled thickness, asperity angle, and incident amplitude on the transmission waves and transmission coefficients. The asperity angle and filled thickness together determine the transmitted waveform and transmission coefficient. The transmitted wave may be cut off when the incident wave amplitude exceeds a threshold value. The transmission coefficient decreases with a different trend with the incident wave amplitude increasing when the asperity angle varies. Compared with planar discontinuity, a filled discontinuity with rough interfaces is more sensitive to the amplitude of the incident wave. The causes of these phenomena are analyzed in detail. In addition, the deformation of the fill material is strongly related to the wave attenuation.     

Effect of irregularity on torsional surface waves in an initially stressed anisotropic porous layer sandwiched between homogeneous and non-homogeneous half-space

The present paper is concerned with the propagation of torsional surface waves in an initially stressed anisotropic porous layer sandwiched between homogeneous and non-homogeneous half-space. We assume the quadratic inhomogeneity in rigidity and density in the lower half-space and irregularity is taken in the form of rectangle at the interface separating the layer from the lower half-space. The dispersion equation for torsional waves has been obtained in a closed form. Velocity equation is also obtained in the absence of irregularity. The study reveals that the presence of irregularity, initial stress, porosity, inhomogeneity and anisotropy factor in the dispersion equation approves the significant effect of these parameters in the propagation of torsional waves in porous medium. It has also been observed that for a uniform media, the velocity equation reduces to the classical result of Love wave.     

Reflection and refraction of SH‐waves at a corrugated interface between two monoclinic elastic half‐spaces

A problem of reflection and transmission of shear wave incident upon a corrugated interface between two monoclinic solid half‐spaces have been investigated. Rayleigh's method of approximation is used to investigate the reflection and transmission coefficients for first and second approximation of the corrugation. For a special interface, closed‐form expressions of these coefficients for the first order approximation of the corrugation are obtained. It is found that these coefficients are functions of amplitude of corrugation, elastic parameters of the media, frequency of the incident wave and angle of incidence. The numerical computations are performed for a specific model and the results obtained are presented graphically. It is found that the reflection and transmission coefficients are strongly influenced by the corrugation and the elastic properties of the media. Results of some earlier workers in this field have been reduced as particular case from the present formulation. Copyright © 2004 John Wiley & Sons, Ltd.     

双重孔隙介质中P1波在自由边界上的反射

对双重孔隙介质波动中的体波及P1波在半空间自由透水和不透水边界上的反射问题进行了详细的理论求解。通过数值计算,对双重孔隙介质中4种体波的弥散和衰减特性,以及不同边界上反射波振幅反射系数随P1波入射角和频率的变化进行了分析。结果表明：（1）P1、P2、P3和S波均具有弥散和衰减特性,P1波传播最快而衰减最慢,P3波传播最慢而衰减最快;（2）随着入射角的增大,反射P2、P3和S波的振幅反射系数均先增大后减小,反射P1波则先减小后增大,且边界透水条件对反射P2和P3波的振幅反射系数影响显著,对P1和S波的影响较小,但透水和不透水边界下振幅反射系数随入射角的变化趋势一致;（3）反射P2、P3波的振幅反射系数均随频率的增大而增大,且不透水边界下反射P2波的振幅反射系数大于透水边界下的反射系数,而反射P3波则相反;（4）随着频率增大,透水边界下反射P1波的振幅反射系数先增大后减小,不透水边界下先减小后增大再减小,反射S波的振幅反射系数随频率变化情况则与P1波相反。     

基于离散颗粒模型的岩体结构面动-静态刚度系数对比的数值模拟研究

岩体工程计算分析中结构面刚度系数是至关重要的力学参数,计算分析的精度和可靠程度与这个参数密不可分,结构面刚度系数取值仍然是一个难点.岩体中应力波传播至结构面处将会发生反射和透射现象,利用应力波透射系数可反演结构面动态刚度系数.本文从细观力学角度运用颗粒离散元方法,开发分段线性接触模型及应力波吸收边界模型,开展宏观岩体中...     

REFLECTION AND TRANSMISSION OF SHEAR WAVES IN MONOCLINIC MEDIA

This paper deals with the propagation, reflection and transmission of shear waves in monoclinic media. The dispersion equation for a monoclinic layer overlying a monoclinic half-space has been obtained and curves are plotted. The amplitude ratios for both the reflected and transmitted waves due to reflection of shear waves at the interface of two monoclinic half-spaces have also been computed and the numerical results are presented graphically. The results are compared with the isotropic case. It has been observed that, in monoclinic media, the amplitude ratios for reflected and transmitted wave increases approximately by 25 and 50 per cent respectively, in comparison to the isotropic case. © 1997 John Wiley & Sons, Ltd.     

含煤地层转换波叠后横波波阻抗反演

叠后地震反演方法适用于纵波和横波地震资料的反演,但对转换波地震资料却不适用。为了实现利用转换波剖面得到横波波阻抗信息,通过建立转换波反射系数和横波反射系数之间的关系,直接从转换波地震数据中获取横波反射系数,把转换波叠后波阻抗反演问题转化为横波波阻抗反演问题。基于测井约束反演原理,利用横波反射系数和转换波解释层位建立初始横波波阻抗模型,优选反演处理参数,实现转换波叠后横波波阻抗反演。含煤地层转换波横波波阻抗反演结果,为进一步研究岩性参数比属性奠定了基础。和横波勘探相比,转换波勘探施工方便,工作成本相对低廉,因此利用转换波叠后地震反演获得横波波阻抗,是一种经济的办法。但由于在计算横波反射系数时采用了近似公式,所以反演结果精度会受到一定影响。     

Magnetoelastic shear waves in an infinite self-reinforced plate

The propagation of magnetoelastic shear waves in an infinite self-reinforced plate is studied in this paper. The frequency spectrum for SH-modes in the plate is drawn and it is shown that if the dimensionless wave number is real, the branches are hyperbolas while for imaginary wave numbers, they are ellipses. It is also observed that for higher modes, the dimensionless group velocity approaches 9·9754 as the wave number increases. The group velocity is normalized with respect to the phase velocity of shear waves. For a homogeneous and isotropic medium, on the other hand, the dimensionless group velocity approaches unity. So it is approximately 10 times larger in case of magnetoelastic shear waves in an infinite self-reinforced plate. It is also noted that in a reinforced medium, if both the values of the coupling parameter and the angle at which the wave crosses the primary magnetic field are increased, the values of the dimensionless frequency and group velocity are decreased.     

Analysis of Wave Propagation Through a Filled Rock Joint

An analytical and experimental study on a longitudinal wave (P-wave) transmission normally across a filled rock joint is presented in this paper. The dynamic property of the filling material for the artificial rock joints is derived from a series of modified split Hopkinson pressure bar (SHPB) tests. The incident and transmitted waves in granitic pressure bars are calculated by wave separations of the strain gauge readings. The incident wave is approximated by a series of half-sinusoidal waves, and an analytical model on wave propagation across a filled rock joint is then deduced. The derived wave transmission coefficients across the filled joint agree very well with those from the test results. Both the analytical and test results show that the wave transmission coefficients are influenced by the mechanical properties and the input energy of the incident waves. Analytical parametric studies with respect to pre-compaction of the filling material, the frequency and amplitude of the incident wave have also been conducted.     

Elastic waves at a corrugated interface between two dissimilar fibre‐reinforced elastic half‐spaces

The reflection and transmission phenomena of elastic waves incident at a corrugated interface between two dissimilar fibre‐reinforced elastic half‐spaces have been analysed. Using Rayleigh method of approximation, the expressions of the reflection and transmission coefficients are obtained in closed form for the plane interface as well as for the first order approximation of the periodic interface ζ = d cos px. All these reflection and transmission coefficients of regular and irregular waves are found to be the functions of angle of incidence and elastic parameters of the media. Moreover, the coefficients of irregularly reflected and transmitted waves are found to be proportional to the amplitude of the corrugated interface and are functions of the frequency of the incident wave. Numerical computations have been performed for a specific model to compute these coefficients and results obtained are shown graphically. The results of Singh and Singh (Sadhana 2004; 29 :249–257) and Ben‐Menahem and Singh (Seismic Waves and Sources. Springer: New York) have been derived from our analysis as particular cases. Copyright © 2006 John Wiley & Sons, Ltd.     

Propagation of Rayleigh waves in fluid‐saturated non‐homogeneous soils with the graded solid skeleton distribution

The propagation characteristic of Rayleigh waves in a fluid‐saturated non‐homogeneous poroelastic half‐plane is addressed. Based on Biot's theory for fluid‐saturated media, which takes the inertia, fluid viscosity, mechanical coupling, compressibility of solid grains, and fluid into account, the dispersion equations of Rayleigh waves in fluid‐saturated non‐homogeneous soils/rocks are established. By considering the shear modulus of solid skeleton variation with depth exponentially, a small parameter, which reflects the relative change of shear modulus, is introduced. The asymptotic solution of the dispersion equation expressing the relationship between the phase velocity and wave number is obtained by using the perturbation method. In order to analyze the effects of non‐homogeneity on the propagation characteristic of Rayleigh waves, the variation of the phase velocity with the wave number is presented graphically and discussed through numerical examples. Copyright © 2016 John Wiley & Sons, Ltd.