Analysis of Stochastic Seismic Wave Interaction with a Slippery Rock Fault

Stochastic seismic wave interaction with a slippery rock fault is studied, based on the principle of conservation of momentum at the wave fronts along the fault. By using the displacement discontinuity method, the wave propagation equations are derived for incident longitudinal-(P-) and shear-(S-) waves, respectively. This is an extension of the study by Li and Ma (2010) for blast-induced wave propagation across a linear rock joint. Stochastic seismic waves are generated from a frequency spectrum and used to analyze the seismic wave interaction with a rock fault having a Coulomb-slip behavior. Parametric studies are carried out to investigate the effect of the intensity and impinging angle of the incident seismic waves on wave propagation across a slippery rock fault. Results show that the transmission of the incident P-wave is almost not affected by the fault, on the contrary, this is not the case for an incident S-wave, due to the occurrence of a relative slip which is related to the impinging angle of the incident S-wave. A quantitative study is presented which is of help in understanding the propagation and attenuation laws of seismic waves in discontinuous rock masses.     

P‐wave transmission across fractures with nonlinear deformational behaviour

Stress wave attenuation across fractured rock masses is a great concern of underground structure safety. When the wave amplitude is large, fractures experience nonlinear deformation during the wave propagation. This paper presents a study on normal transmission of P‐wave across parallel fractures with nonlinear deformational behaviour (static Barton–Bandis model). The results show that the magnitude of transmission coefficient is a function of incident wave amplitude, nondimensional fracture spacing and number of fractures. Two important indices of nondimensional fracture spacing are identified, and they divide the area of nondimensional fracture spacing into three parts (individual fracture area, transition area and small spacing area). In the different areas, the magnitude of transmission coefficient has different trends with nondimensional fracture spacing and number of fractures. In addition, the study reveals that under some circumstances, the magnitude of transmission coefficient increases with increasing number of fractures, and is larger than 1. Copyright © 2006 John Wiley & Sons, Ltd.     

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.     

A Numerical Study on Wave Transmission Across Multiple Intersecting Joint Sets in Rock Masses with UDEC

This paper presents a numerical study on wave transmission across jointed rock masses with UDEC, where multiple intersecting joint sets exist. The capability of UDEC of studying wave transmission across rock joints is validated through comparison with analytical solutions and experimental data. Through parametric studies on wave transmission across jointed rock masses, it is found that joint mechanical and spatial parameters including joint normal and shear stiffnesses, nondimensional joint spacing, joint spacing ratio, joint intersecting angle, incident angle, and number of joint sets together determine the wave transmission. And for P wave incidence, compared with other parameters, joint normal stiffness, nondimensional joint spacing, and joint intersecting angle have more significant effects on wave transmission. The physical reasons lying behind those phenomena are explained in detail. Engineering applications and indications of the modeling results are also mentioned.     

Validation study of the distinct lattice spring model (DLSM) on P-wave propagation across multiple parallel joints

A validation study of the distinct lattice spring model (DLSM) for wave propagation problems is performed. DLSM is a microstructure-based numerical model, which is meshless and has advantages in modelling dynamic problems where stress wave propagation is important. To verify the applicability of DLSM to modelling wave propagation through a discontinuous medium, the virtual wave source (VWS) method is used to obtain analytical solutions for wave propagation across a jointed rock mass. Numerical modelling results of the commercial code UDEC are selected as the reference. The effects of particle size and lattice rotation angle on wave propagation are first studied. Then, the results of wave transmission across a single joint with a different joint stiffness and across multiple parallel joints with different joint spacings are derived with DLSM, UDEC and VWS. These results are in good agreement with each other. Therefore, the capability of DLSM to model P-wave propagation across jointed rock mass is verified, which provides confidence for the further application of DLSM to modelling more complex problems.     

节理面透射模型及其隔振性能研究

振动波入射到节理界面上,与软弱断层一样会产生波场的分解,能够起到隔振作用。根据线性变形节理的透反射规律研究发现,隔振后各点的波场可由透射P波与透射S波叠加而成,其能量分配服从Zoeppritz方程,不连续位移造成了透射波的衰减;法向刚度增加透射纵波系数会增加,而切向刚度与透射纵波无关,法向与切向刚度对透射S波均有影响,刚度增加透射S波振幅逐渐降低,逐步逼近于0;由于透射后波场相互干涉,不同波场延时不同,透射合成波与入射波相比,振幅降低,持续时间增加,波形尾段有畸变;总能量系数与节理面的刚度无关,而随入射角增加不断降低。     

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

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

双重孔隙介质中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波相反。     

岩体中应力波幅值随时空衰减的关系

把岩石抽象为Kelvin-Voigt黏弹性体,应用复数理论建立了应力波在岩石中传播时的振动频率和波动频率的关系,在此基础上,研究了岩石和含有一组平行结构面的岩体中应力波幅值随时间和随空间衰减的关系。研究表明,应力波在岩石的传播过程中,频率是影响应力波幅值衰减的的主要因素;应力波在含有1组平行结构面岩体的传播过程中,由于结构面的存在,影响了应力波能量的传播,其幅值随时空衰减系数主要由结构面间距及其反射系数决定,结构面间距减小和反射系数增加,幅值随传播距离的衰减系数增加和随时间的衰减系数减小,而结构面的黏性系数增加,幅值随传播距离的衰减系数增加,但随时间的衰减系数变化很小。上述结论对岩体结构探测具有重要的指导意义。     

黏弹性节理岩体的位移不连续模型研究

选用三参数标准线性固体作为岩石本构,提出了一种考虑岩体黏弹性的位移不连续模型;根据一维黏弹性波的特征线法,推导了节理处质点速度、应力和应变递推公式。首先,基于分离式霍普金森压杆（SHPB）对砂层进行试验,得到其应力-应变关系,并换算出砂层节理的法向刚度;接着,通过一维强间断黏弹性波的波速公式、高频波衰减系数以及任一频率下的衰减系数,确定出数值算法中的三参数。最后,基于自制的摆锤装置,探讨了一维应力波在节理岩体中的传播规律,试验中以两根长1 000 mm、直径为68.50 m的岩杆作为入射和透射杆,以3 mm砂层模拟节理。试验和数值结果吻合度良好,进一步验证了该方法的可靠性。     

Transmission of Elastic P-waves across Single Fractures with a Nonlinear Normal Deformational Behavior

Summary This paper presents a theoretical study on normally incident elastic P-wave transmission across single dry fractures with a nonlinear normal deformational behavior. The effects of nonlinear fracture normal behavior on P-wave transmission are examined without the mixture of fracture shear behavior. The linear displacement discontinuity model for wave propagation across fractures is extended to a nonlinear model – the hyperbolic elastic model (BB model). Numeric solutions of magnitudes of transmission (|T _non|) and reflection (|R _non|) coefficients, for normally incident P-wave transmission across the nonlinear deformable fractures, are obtained and related to the closure behavior of fractures. Parametric studies are conducted to acquire an insight into the effects of the nonlinear fracture normal deformation on P-wave transmission, in terms of initial normal stiffness and the ratio of current maximum closure to maximum allowable closure of the fractures, as well as the incident wave amplitude and frequency. Comparisons between the linear and nonlinear models are presented. It is shown that, |T _lin| and |R _lin| for the linear model are special solutions of |T _non| and |R _non| for the nonlinear model, when the incident wave amplitude is so low that the current maximum closure of fracture incurred during the wave transmission is much smaller, relative to the maximum allowable closure. In addition, the nonlinear fracture behavior gives rise to a phenomenon of higher harmonics during the wave transmission across the fracture. The higher harmonics contribute to the increase of |T _non| from |T _lin|.     

P波在弹性介质与饱和冻土介质分界面上的透反射问题研究

基于弹性波在冻结饱和多孔介质与单相弹性介质中的传播理论,研究了平面P波在饱和冻土介质与单相弹性介质分界面上的透反射问题。利用Helmholtz矢量分解定理,根据分界面上的边界条件,获得了平面P波从单相弹性介质入射到饱和冻土介质分界面上透反射振幅比的理论表达式。通过数值计算,分析了在不同入射频率、胶结参数、孔隙率、饱和度和接触参数下,弹性波的透反射振幅比随入射角变化的关系。研究结果表明：P波从单相弹性介质垂直入射到饱和冻土介质分界面上时只有反射P波和3种透射P波产生,当掠入射时只产生反射而没有透射现象发生。入射频率、胶结参数、孔隙率、饱和度以及接触参数等参数对反射波和透射波的振幅比影响显著。     

非均匀介质地震AVO响应模拟及分析

地震AVO （amplitude versus offset）技术是一项利用振幅信息研究岩性、检测油气的地震勘探技术。常规方法基于Zoeppritz方程计算模型界面处的反射系数,而实际地下非均匀介质中地震反射特征不仅与入射角度、物性差异有关,还与入射波频率、地层厚度、薄互层结构等因素有关。为此,应用传播矩阵理论充分考虑与这些因素有关的调谐干涉等传播效应,针对实际数据计算非均匀地下介质的高精度合成地震记录,对比Zoeppritz方程、Shuey二项近似方程、Shuey三项近似方程以及传播矩阵算法的模拟效果。研究发现：在小角度入射时Zoeppritz方程、Shuey二项近似方程和Shuey三项近似方程的反射振幅和波形基本一致,大角度入射时Zoeppritz方程与Shuey三项近似方程接近;Zoeppritz算法的模拟结果在小角度入射和浅层情况下与传播矩阵算法差别较小,而在大角度入射和深层情况下与传播矩阵差别较大,说明层间多次波的调谐干涉等传播效应不可忽略。     

非饱和土地基中复合多层波阻板对S波的隔离效应

基于弹性波在非饱和多孔介质与单相弹性介质中的传播理论,考虑在非饱和土地基中设置一定厚度的复合多层波阻板（复合多层波阻板以3层为例）,利用Helmholtz矢量分解定理,推导了非饱和土地基中S波通过复合多层波阻板的透射、反射振幅比的解析解。通过数值算例,分析了层间波阻板剪切模量和密度等物理力学参数对非饱和土地基中S波通过复合多层波阻板时传播特性的影响规律。结果表明：复合多层波阻板中层间波阻板材料的剪切模量和密度对透反射系数影响显著。复合多层波阻板是一种有效的隔振屏障,严格控制层间波阻板的剪切模量和密度可以获得最佳隔振效果,这为复合多层波阻板在地基振动控制领域中的应用提供理论指导。     

Effect of pore connectivity on reflection amplitudes of an inhomogeneous wave in a composite porous solid saturated by two immiscible fluids

Present paper aims to study the phenomenon of reflection and transmission when an inhomogeneous wave strikes some discontinuity in a composite porous medium saturated by two immiscible viscous fluids. The incident wave splits into six reflected and six transmitted waves at the interface. All reflected and transmitted waves are inhomogeneous in nature with different directions of propagation vector and attenuation vector. A dimensionless parameter \(\varsigma \in [0, 1]\) is introduced to represent the extent of connection among the pores at the interface. Expression of Umov–Poynting vector is derived to obtain energy flux vector. Continuity of energy flux vector at the interface gives the required boundary conditions for the system. Connecting parameter \(\varsigma \) is also employed in boundary conditions to model the partial connection of pores at the interstices of two media. For numerical discussion we consider a porous medium composed of sandstone and ice, saturated with oil and water. The effect of parameter \(\varsigma \) and angle of incidence is determined numerically on the amplitude and the energy ratios of reflected and transmitted waves.     

位移不连续模型在黏弹性节理刚度计算的应用

考虑节理质量对应力波传播的影响,运用波的位移位函数推导了谐波在厚黏弹性节理的透、反射系数计算公式,采用波形相关系数描述子波穿过黏弹性节理的波形变化,讨论具有一定厚度的黏弹性节理简化位移不连续模型的适用条件。设厚黏弹性节理模型和位移不连续模型的透射波波形相关系数为0.9时对应的节理厚度为临界厚度,岩体与节理的阻抗比对临界厚度的影响很小;临界厚度随子波中心频率增大呈负指数减小,入射角越大,临界厚度随中心频率减小得越慢。试验数据分析表明：当节理厚度为0.03 m时,采用位移不连续模型和厚黏弹性模型计算得到的节理力学参数非常接近,随节理厚度和子波中心频率增加,运用位移不连续模型的计算结果偏差越大,试验结果与理论分析是一致的。     

P波作用下饱和土中输流管道动力响应

根据Biot波动理论,采用复变函数与多级坐标法,求解了P波作用下饱和土体中地下输流管道的波动散射方程,分析了管道中流体介质性质、入射波角度及管道埋深等对地下输流管道周边动应力集中系数及孔压集中系数分布的影响。计算结果表明：在低频弹性波入射时,管道周边动应力集中系数及孔压集中系数分布相对均匀,而随入射频率的增加,其分布变得复杂化,但其峰值有所减小;中低频波入射作用下,管道内为水、石油介质时,动应力集中系数和孔压集中系数较空气介质时小,而在高频波作用时情况相反;对于管道内流体介质为水时,随着入射角度的变化,应力集中分布也沿一定角度的方向发生转动,入射波自管道下方垂直入射时,管道周边动应力集中系数峰值相对较大;随埋深的增加,动应力集中系数和孔压集中系数均呈震荡减小的趋势。     

Multi-Modal Monitoring of Slip Along Frictional Discontinuities

Seismic wave transmission and digital image correlation (DIC) are employed to study slip processes along frictional discontinuities. A series of biaxial compression experiments are performed on gypsum specimens with non-homogeneous contact surfaces. The specimens are composed of two blocks with perfectly mated contact surfaces with a smooth surface with low frictional strength on the upper half and a rough surface with high frictional strength on the lower half. Compressional, P, and shear, S, wave pulses were transmitted through the discontinuity while digital images of the specimen surface were acquired during the test. A distinct peak in the amplitude of transmitted wave occurs prior to the peak shear strength and is considered a “precursor” to the failure. Precursors indicate that slip initiates from the smooth surface and extends to the rough surface as the shear load is increased. From the DIC data, slip is identified as a jump in the displacement field along the fracture that initiates from the smooth surface and propagates to the rough surface. Precursors are associated with an increase in the rate of slip across the discontinuity and are a measure of the reduction in the fracture shear stiffness.     

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.     

Analytical study of ground motion caused by seismic wave propagation across faulted rock masses

Studying seismic wave propagation across rock masses and the induced ground motion is an important topic, which receives considerable attention in design and construction of underground cavern/tunnel constructions and mining activities. The current study investigates wave propagation across a rock mass with one fault and the induced ground motion using a recursive approach. The rocks beside the fault are assumed as viscoelastic media with seismic quality factors, Q_p and Q_s. Two kinds of interactions between stress waves and a discontinuity and between stress waves and a free surface are analyzed, respectively. As the result of the wave superposition, the mathematical expressions for induced ground vibration are deduced. The proposed approach is then compared with the existing analysis for special cases. Finally, parametric studies are carried out, which includes the influences of fault stiffness, incident angle, and frequency of incident waves on the peak particle velocities of the ground motions.