Next Article in Journal
OpenMP Implementation of a Novel Potential-Field-Data Source-Growth-Based Inversion Approach for 3D Salt Imaging in Deepwater Gulf of Mexico
Previous Article in Journal
Pipelined Dynamic Scheduling of Big Data Streams
Previous Article in Special Issue
Structural Damage Detection Based on Real-Time Vibration Signal and Convolutional Neural Network
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stress Wave Propagation through Rock Joints Filled with Viscoelastic Medium Considering Different Water Contents

1
Key Laboratory of Shale Gas and Geoengineering, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
2
Innovation Academy for Earth Science, Chinese Academy of Sciences, Beijing 100029, China
3
College of Earth and Planetary Science, University of Chinese Academy of Sciences, Beijing 100049, China
4
State Key Laboratory of Lithospheric Evolution, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(14), 4797; https://doi.org/10.3390/app10144797
Submission received: 4 June 2020 / Revised: 7 July 2020 / Accepted: 9 July 2020 / Published: 13 July 2020
(This article belongs to the Special Issue Computational Methods in Vibration Problems and Wave Mechanics)

Abstract

:
A rock mass often contains joints filled with a viscoelastic medium of which seismic response is significant to geophysical exploration and seismic engineering design. Using the propagator matrix method, an analytical model was established to characterize the seismic response of viscoelastic filled joints. Stress wave propagation through a single joint highly depended on the water content and thickness of the filling as well as the frequency and incident angle of the incident wave. The increase in the water content enhanced the viscosity (depicted by quality factor) of the filled joint, which could promote equivalent joint stiffness and energy dissipation with double effects on stress wave propagation. There existed multiple reflections when the stress wave propagated through a set of filled joints. The dimensionless joint spacing was the main controlling factor in the seismic response of the multiple filled joints. As it increased, the transmission coefficient first increased, then it decreased instead, and at last it basically kept invariant. The effect of multiple reflections was weakened by increasing the water content, which further influenced the variation of the transmission coefficient. The water content of the joint filling should be paid more attention in practical applications.

1. Introduction

The study on the seismic response of the rock mass is of great significance to geophysical exploration and seismic engineering design. Usually, a rock mass contains joints that are discontinuous interfaces and that significantly influence mechanical behaviors of the rock mass [1,2,3,4,5,6,7]. Due to the weathering effect, the natural joint was often filled with loose and soft materials such as sand with pore space containing normal air with high moisture content which significantly controlled the joint mechanical behaviors [1,2,8,9,10]. It was found that the thickness and water content of the filling had great influences on the seismic response of the filled joint [11,12,13]. It was found that saturated soils behaved viscoelastic in response to the dynamic loads, which can be characterized by the Kelvin viscoelastic model [14,15]. Previous studies also showed that the viscosity had an important effect on the dynamic response of the joint and could also dissipate wave energy [16,17,18,19,20,21,22,23]. Meanwhile, when the density of the filled material could not be neglected compared with the rock density, the initial mass of the filled joint could affect wave propagation [20]. Therefore, it is crucial to investigate the seismic response of the joints filled with a viscoelastic medium.
Many analytical models were proposed to characterize the seismic response of the filled joints. Ref. [24] applied the displacement discontinuity model (DDM) to study wave propagation through a filled rock joint. The model assumed that the stress was continuous, but the displacement was not when the stress wave impinged on the filled joint which was controlled by the fracture compliance [25]. Based on the DDM, [26] proposed a three-phase medium model for the filled rock joint considering the interaction of the three-phase system consisting of gas, liquid and granular solid. Ref. [20] proposed a displacement and stress discontinuity model (DSDM) to study wave propagation through filled joints considering the initial mass effect of the filled material. This model assumes that not only the displacement but also the stress is discontinuous in the front and rear interface of the filled joint. However, the above studies mainly have two limitations: (1) Wave reflection in the filled layer was neglected when it needs to be considered as the ratio of filled thickness to the wavelength is relatively large; (2) There are many parameters in these analytical models which are difficult to determine by laboratory tests, such as the viscosity and stiffness of the filled joint. To overcome these problems, a simple and feasible way is to model the filled joint as a layer of independent continuum solid medium. Properties of the filled joint are controlled by those of the filled medium such as the density, filled thickness (aperture), wave velocity and quality factor (Q) that can be readily determined. Furthermore, multiple reflections in the filled layer can be taken into account. It was found that a joint or an interface between two solids can be represented as a thin-layer interface [27]. Ref. [28] pointed out that a thin-layer interface or a zero-thickness interface for a joint should be essentially the same from the physical point of view. The thin-layer medium model (TLMM) was also proved to be reasonable and feasible to characterize stress wave propagation through filled joints in previous studies [29]. They found that predicted results of stress wave propagation through joints by the purely elastic TLMM are nearly the same as the zero-thickness interface model if the ratio of the wavelength to the filled thickness is large enough. Zhu et al., in 2012 [30] investigated the propagation of the stress waves normally impinging on a single rock joint filled with a viscoelastic medium. By modeling the interface between two solids as a thin viscoelastic layer with stiffness and inertia term, wave propagation through a single filled joint was addressed by [31]. However, thus far, analytical models based on the TLMM were rarely reported to investigate seismic response of viscoelastic filled joints as the water content of the filling varied.
The purpose of this paper is to establish an analytical model on the seismic response of the rock joints filled with a viscoelastic medium, considering the water content, the filled thickness (aperture), viscosity of the filled material (quality factor), the joint spacing (related to joint density) as well as the incident angle and frequency of the incident wave.

2. Method

When the stress wave impinges on a viscoelastic filled joint, only one portion of the stress wave can pass through the filled joint while the other portion will be reflected. Meanwhile, wave conversion happens for the case of the oblique incidence at the joint [18], as shown in Figure 1. There are multiple reflections at joint interfaces [32,33]. The overall transmitted wave, after propagating through a set of filled joints, is the superposed wave consisting of the first-arriving transmitted wave and the later-arriving transmitted waves induced by multiple reflections at each joint. In previous studies, three types of analytical methods were mainly proposed to solve the multiple-reflection problem, i.e., the method of characteristics (MC) [33,34]; the virtual wave source method (VWSM) [35,36]; the recursive method (RM) [21,22,23,30]. MC was mainly applied to solve the 1-D problem of stress wave propagation through joints. For the VWSM, it is necessary to clarify the analytical solution of the wave propagation through a single joint. The propagator matrix method (PMM) is one of the RMs, that was widely applied to solve wave propagation problems of stratified medium [37,38]. The algorithm needs not to identify the analytical solution of wave propagation through a single interface or joint. Ref. [39] first introduced the PMM to study the propagation of the stress wave through unfilled joints satisfying purely elastic deformational behavior. Ref. [21] extended the PMM from the purely elastic case to the viscoelastic case and studied the seismic response of the rock mass containing viscoelastic unfilled joints. Combining with the TLMM, the PMM is also used to derive formulas of wave propagation through rock mass containing joints filled with viscoelastic medium in this paper.
As shown in Figure 1, a 2-D plane is considered and both the local coordinates (xi, zi) and the global coordinates (X, Z) are applied. The rock mass is divided by the multiple parallel filled joints with identical joint spacing hr as shown in Figure 1a. The filled joints are numbered from 1 to n and the rock layers from 1 to n + 1, respectively, with Z increasing towards the lower layers. When a P wave impinges on the upper boundary, wave conversions, reflection and transmission occur. θ denotes the incident and the reflection angle of the incident P wave and φ represents the reflection angle of the converted SV wave. An element surrounded by the dash frame (see Figure 1a) is chosen to derive the stress and displacement relations between two adjacent rock layers. The sandwiched element consists of two layers rocks and one filled layer with a thickness of hf, which are numbered as i, f and i + 1 from the upper to the low layer, respectively.
Each rock layer is assumed to be a linearly elastic, homogeneous and isotropic medium with the density of ρr, P- and S-wave velocity of Vpr and Vsr. Deformation of each filled medium is characterized by the Kelvin viscoelastic model, which has the density of ρf, P- and S-wave velocity of Vpp and Vsp, and P- and S-wave quality factor of Qp and Qs, where Vpp and Vsp denote the purely elastic wave velocity when the viscosity is not considered. The interface between the rock and the filled medium is assumed to be welded, in which both stresses and displacements are continuous. Two local coordinates xf-of-zf and xi+1-oi+1-zi+1 are applied as shown in Figure 1b. The symbols “+” and “−“ represent the upper boundary and the lower boundary of a layer (e.g., in the (i + 1) th layer, “+” represents the upper boundary (zi+1 = 0), whereas “−“ represents the low boundary (zi+1 = hr).
For the linearly elastic rock layer, the 2-D wave equation can be written as:
ρ r 2 u x t 2 = σ x x + τ x z z
ρ r 2 u z t 2 = τ x z x + σ z z
σ z = λ r u x x + ( λ r + 2 μ r ) u z z
τ x z = μ r ( u z x + u x z )
where ui denotes the particle displacement and σij denotes the stress, λr and μr are Lamé constants for rock and λr = ρr (Vpr)2-2ρr (Vsr)2, μr = ρr (Vsr)2.
For the Kelvin viscoelastic filled layer, the 2-D wave equation can be written as:
ρ f 2 u x t 2 = σ x x + τ x z z
ρ f 2 u z t 2 = τ x z x + σ z z
σ z t = λ p 2 u x t x + ( λ p + 2 μ p ) 2 u z t z + λ f 3 u x x t 2 + ( λ f + 2 μ f ) 3 u z z t 2
τ x z t = μ p 2 u z t x + μ p 2 u x t z + μ f 3 u z x t 2 + μ f 3 u x z t 2
where λp and μp are the purely elastic Lamé constants of the filled layer when the viscosity is not considered, λp = ρf (Vpp)2-2ρf (Vsp)2, μp = ρf (Vsp)2 and λfʹ and μfʹ denote the viscoelastic modulus for the filled medium. Equations (1)–(8) are transformed into frequency domains by Fourier transform.
For the rock layer,
ρ r ω 2 u ˜ x = ( σ ˜ x x + τ ˜ x z z ) ,
ρ r ω 2 u ˜ z = ( τ ˜ x z x + σ ˜ z z )
σ ˜ z = λ r u ˜ x x + ( λ r + 2 μ r ) u ˜ z z
τ ˜ x z = μ r ( u ˜ x z + u ˜ z x ) ,
where, variables with the symbol “~” over them denotes the Fourier transform.
For the filled layer,
ρ f ω 2 u ˜ x = ( σ ˜ x x + τ ˜ x z z )
ρ f ω 2 u ˜ z = ( τ ˜ x z x + σ ˜ z z )
σ ˜ z = ( λ p + i ω λ f ) u ˜ x x + [ ( λ p + 2 μ p ) + i ω ( λ f + 2 μ f ) ] u ˜ z z
τ ˜ x z = ( μ p + i ω μ f ) ( u ˜ x z + u ˜ z x )
where ω is angular frequency and i = 1 .
The quality factors related to energy attenuation can be defined as [40]:
Q p = λ p + 2 μ p λ f + 2 μ f 1 ω Q s = μ p μ f 1 ω
When Qp→∞ and Qs→∞, Equations (13)–(16) represent the pure elastic case without viscosity energy loss. When Qp→0 and Qs→0, the filled medium has a very large viscosity. Based on Equations (15)–(17), we can obtain corresponding complex equivalent Lamé constants and an equivalent wave velocity of the filled medium as follows:
λ f = λ p + i ω λ f = λ p [ 1 + i Q s V p p 2 2 Q p V s p 2 Q p Q s ( V p p 2 V s p 2 ) ] μ f = μ p + i ω μ f = μ p ( 1 + i 1 Q s )
V p f = λ f + 2 μ f ρ f V s f = μ f ρ f
Substituting the Equation (18) into Equations (15) and (16), we can obtain the following equations:
σ ˜ z = λ f u ˜ x x + ( λ f + 2 μ f ) u ˜ z z
τ ˜ x z = μ f ( u ˜ x z + u ˜ z x )
It can be seen that Equations (20) and (21) have the same form as those of purely elastic wave equations. Based on the Helmholtz decomposition theory [41], the 2-D dynamic equation can be expressed as two displacement potential functions. Assume that the φ ˜ i (the displacement potential of the P wave) and ψ ˜ 2 i (the second component of the displacement potential of the S wave) are the potential functions of the ith layer in the frequency domain which can be expressed by following two Equations [21,39].
φ ˜ i ( x , z , ω ) = U p i e i ( k x i x + k z p i z ) + D p i e i ( k x i x k z p i z )
ψ ˜ 2 i ( x , z , ω ) = U s i e i ( k x i x + k z s i z ) + D s i e i ( k x i x k z s i z )
where kxi, kzpi are the P-wave numbers of the ith layer in the x and z directions, and kzsi is the S-wave number of the ith layer in the z direction. As shown in Figure 1b, the sandwiched element consists of two types of materials in turn. When the wave propagates in the rock layer, the wave numbers kxi, kzpi and kzsi can be expressed as kxr = kpr sin θ, kzpr = kpr cos θ, k z s r = ( ω / V s r ) 2 k x r 2 , where kpr = ω/Vpr is the wave number of the P wave. Similarly, when the wave propagates in the filled layer, the wave numbers kxi, kzpi and kzsi can be expressed as kxf = kpf sin θ, kzpf = kpf cos θ and k z s f = ( ω / V s f ) 2 k x f 2 , where kpf = ω/Vpf is the wave number of the P-wave propagation through the filled medium.
The particle displacement can be written as:
u ˜ x i = φ ˜ i x ψ ˜ 2 i z
u ˜ z i = φ ˜ i z + ψ ˜ 2 i x
Substituting Equations (22) and (23) into Equations (24) and (25), we can form:
u ˜ x i = [ ( i k x i e i k z p i z ) U p i + ( i k x i e i k z p i z ) D p i + ( i k z s i e i k z s i z ) U s i + ( i k z s i e i k z s i z ) ] D s i ] e i k x i x
u ˜ z i = [ ( i k z p i e i k z p i z ) U p i + ( i k z p i e i k z p i z ) D p i + ( i k x i e i k z s i z ) U s i + ( i k x i e i k z s i z ) ] D s i ] e i k x i x
Based on the constitutive Equations (20) and (21) and Equations (26) and (27), we can also form:
σ ˜ z i = [ ( λ i ω 2 V p i 2 2 μ i k z p i 2 ) e i k z p i z U p i + ( λ i ω 2 V p i 2 2 μ i k z p i 2 ) e i k z p i z D p i + 2 μ i k x i k z s i e i k z s i z U s i + 2 μ i ( k x i k z s i ) e i k z s i z D s i ] e i k x i x
τ ˜ x z = [ 2 μ i k x i k z p i e i k z p i z U p i 2 μ i k x i k z p i e i k z p i z D p i μ i ( k x i 2 k z s i 2 ) e i k z s i z U s i μ i ( k x i 2 k z s i 2 ) e i k z s i z D s i ] e i k x i x
Equations (26)–(29) can be written as the matrix form:
[ u ˜ x u ˜ z σ ˜ z τ ˜ x z ] i = [ i k x i e i k z p i z i k x i e i k z p i z i k z s i e i k z s i z i k z s i e i k z s i z i k z p i e i k z p i z i k z p i e i k z p i z i k x i e i k z s i z i k x i e i k z s i z ( λ i ω 2 V p i 2 2 μ i k z p i 2 ) e i k z p i z ( λ i ω 2 V p i 2 2 μ i k z p i 2 ) e i k z p i z 2 μ i k x i k z s i e i k z s i z 2 μ i k x i k z s i e i k z s i z 2 μ i k x i k z p i e i k z p i z 2 μ i k x i k z p i e i k z p i z μ i ( k x i 2 k z s i 2 ) e i k z s i z μ i ( k x i 2 k z s i 2 ) e i k z s i z ] [ U p D p U s D s ] i e i k x i x = P i Q i [ U p D p U s D s ] i e i k x i x
For the rock layer,
P i = P r = [ i k x r i k x r i k z s r i k z s r i k z p r i k z p r i k x r i k x r ( λ r ω 2 V p r 2 2 μ r k z p r 2 ) ( λ r ω 2 V p r 2 2 μ r k z p r 2 ) 2 μ r k x r k z s r 2 μ r k x r k z s r 2 μ r k x r k z p r 2 μ r k x r k z p r μ r ( k x r 2 k z s r 2 ) μ r ( k x r 2 k z s r 2 ) ]
and
Q i = Q r = [ e i k z p r z 0 0 0 0 e i k z p r z 0 0 0 0 e i k z s r z 0 0 0 0 e i k z s r z ]
the subscripts ‘i’ denote the layer number.
For the filled layer,
P i = P f = [ i k x f i k x f i k z s f i k z s f i k z p f i k z p f i k x f i k x f ( λ f ω 2 V p f 2 2 μ f k z p f 2 ) ( λ f ω 2 V p f 2 2 μ f k z p f 2 ) 2 μ f k x f k z s f 2 μ f k x f k z s f 2 μ f k x f k z p f 2 μ f k x f k z p f μ f ( k x f 2 k z s f 2 ) μ f ( k x f 2 k z s f 2 ) ]
and
Q i = Q f = [ e i k z p f z 0 0 0 0 e i k z p f z 0 0 0 0 e i k z s f z 0 0 0 0 e i k z s f z ]
Equation (30) establishes the relation between stress and displacement and the potential amplitudes in each layer. Note that the matrix Qi is related to z coordinate. As shown in Figure 1b, two local coordinate systems (xf-of-zf and xi + 1-oi+1-zi+1) are adopted, respectively. For the upper boundary of a layer, z = 0 and Qi can be expressed as Qi+. For the low boundary of a layer, z = hi and Qi can be expressed as Qi. Therefore, we can form:
Q r + = Q f + = [ 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 ] Q r = [ e i k z p r h r 0 0 0 0 e i k z p r h r 0 0 0 0 e i k z s r h r 0 0 0 0 e i k z s r h r ]
and
Q f = [ e i k z p f h f 0 0 0 0 e i k z p f h f 0 0 0 0 e i k z s f h f 0 0 0 0 e i k z s f h f ]
The relation of the potential amplitudes between two adjacent rock layers can be derived from the sandwiched element as shown in Figure 1b. As mentioned before, the stress and the displacement at the interface between the rock layer and the filled layer are continuous. In addition, based on Equation (30), we can easily form:
[ u ˜ x u ˜ z σ ˜ z τ ˜ x z ] i = [ u ˜ x u ˜ z σ ˜ z τ ˜ x z ] f + = P f Q f + [ U p D p U s D s ] f e i k x f x
[ u ˜ x u ˜ z σ ˜ z τ ˜ x z ] f = [ u ˜ x u ˜ z σ ˜ z τ ˜ x z ] i + 1 + = P f Q f [ U p D p U s D s ] f e i k x f x
[ u ˜ x u ˜ z σ ˜ z τ ˜ x z ] i + 1 = P r Q r [ U p D p U s D s ] i + 1 e i k x r x
From (31)–(33), we can form:
[ u ˜ x u ˜ z σ ˜ z τ ˜ x z ] i + 1 = R i [ u ˜ x u ˜ z σ ˜ z τ ˜ x z ] i
where, R i = P r Q r P r 1 P f Q f P f 1 .
Note that Ri can be regarded as the stiffness matrix of the filled joint which is related to the thickness, density, P- and S- wave velocity, and viscosity (quality factor) of the filled medium. It has the similar form compared with that of the unfilled joints [21,39], which additionally considers the effect of the wave impedance of the filling.
The relation of Equation (34) can be extended from the first layer to the nth layer by chain rule:
[ u ˜ x u ˜ z σ ˜ z τ ˜ x z ] n = R n 1 R 2 R 1 [ u ˜ x u ˜ z σ ˜ z τ ˜ x z ] 1
When the wave impinges on the lower boundary of the first layer, the propagation distance is 0, therefore:
[ u ˜ x u ˜ z σ ˜ z τ ˜ x z ] 1 = P r [ U p D p U s D s ] 1 e i k x r x
Meanwhile, after the wave propagating through n filled joints, the wave propagates into the (n + 1)th layers in fact. The displacement and stress vector at the upper boundary of the (n + 1) rock layer should be expressed.
[ u ˜ x u ˜ z σ ˜ z τ ˜ x z ] n = [ u ˜ x u ˜ z σ ˜ z τ ˜ x z ] f + = P f Q f + [ U p D p U s D s ] f e i k x f x
[ u ˜ x u ˜ z σ ˜ z τ ˜ x z ] n + 1 + = [ u ˜ x u ˜ z σ ˜ z τ ˜ x z ] f = P f Q f [ U p D p U s D s ] f e i k x f x
From Equations (37) and (38), we can form:
[ u ˜ x u ˜ z σ ˜ z τ ˜ x z ] n = P f Q f 1 P f 1 [ u ˜ x u ˜ z σ ˜ z τ ˜ x z ] n + 1 + = P f Q f 1 P f 1 P r Q r + [ U p D p U s D s ] n + 1 e i k x r x
Substituting Equations (36) and (39) into (35), we can form:
[ U p D p U s D s ] n + 1 = P r 1 P f Q f P f 1 R n 1 R 2 R 1 P r [ U p D p U s D s ] 1
For the (n + 1)th layer, there only exists down-going P- and S- waves, therefore Up = 0 and Us = 0. Because only the P-wave impinges on the first layer, Ds is equal to 0. Based on Equation (40), the reflected and transmitted waves can be calculated. We can obtain reflection–transmission coefficients that are defined as ratios of reflected and transmitted wave amplitudes to the incident wave amplitude. Specifically, the reflection–transmission coefficients can be written as Rp→p = max(abs(Up1)), Rp→s = max(abs(Us1)), Tp→p = max(abs(Dpn+1)) and Tp→s = max(abs(Dsn+1)). For these symbols, the first subscript letter represents the incident wave and the second subscript letter stands for the transmitted or reflected wave. Details can be found in [21]. During calculation, the wave has to be first transferred into frequency domain signals by the fast Fourier transform (FFT) to calculate the reflected and transmitted waves. Finally, these waves in the frequency domain are transferred into signals in the time domain by the inverse fast Fourier transform (IFFT). In the following section, the experimental incident wave is first used as the input incident wave for model validation and then one cycle of sinusoidal wave is used as the input incident wave for parametric studies. According to [21], the start time to oscillate needed to have enough delay rather than from 0 to eliminate the aliasing effect as much as possible in the processing of FFT. Meanwhile, some amplitude compensations are also needed.

3. Results

3.1. Experimental Calibration

In this section, the analytical model is calibrated according to the experimental data from [12]. In their study, a layer of quartz sand with a fixed thickness of 2 mm but varied water content of 0%, 2%, 5%, 10%, and 15% (saturated) was sandwiched between two norite bars to simulate filled joints with different water contents. The stress wave was induced by a split Hopkinson rock bar (SHRB) system. The quartz sand had a porosity of around 0.4 and material density of 2620 kg∙m−3 (solid density). The density of fillings with different water contents can be calculated and listed in Table 1. The norite had a P-wave velocity of 6000 m∙s−1 and density of 2900 kg∙m−3 [12,42]. The S-wave velocity of the norite was taken as 2/3 of its P-wave velocity (4000 m∙s−1). The experimental incident wave was adopted as the incident wave to predict transmitted waves after propagating through filled joints with different water contents. The start time of the wave signal to oscillate has a delay of around 0.0012 s compared with that in [12]. Table 1 shows the calibrated parameters in cases of different water contents. Except the density of the filled medium, other parameters all decrease with increasing water content. According to Section 2, the decrease in the quality factor indicates the increase in the viscosity of the filled joint.
Figure 2 shows the experimental incident wave, transmitted waves through a single filled joint with different water contents and predicted transmitted waves. The amplitude of the experimental transmitted wave obviously decreases with increasing water content. Predicted transmitted waves are basically consistent with experimental results except that their durations are a bit larger than those in the experiment. Figure 3 shows a comparison of the transmission coefficient between the experimental and predicted results. We can see that the experimental transmission coefficient decreases with increasing water content, which agrees well with the experimental result. From Figure 2 and Figure 3, it can be inferred that our analytical model basically has the capability to characterize the seismic response of the filled joint considering different water contents.

3.2. Parametric Studies

In this section, reflection–transmission coefficients when a P wave obliquely impinges on a single filled joint and set parallel filled joints with identical spacing were investigated, respectively, considering the incident angle θ and the frequency f of the incident wave, as well as the filled thickness hf and the quality factor of the filled material Qp and Qs. One of these parameters was studied while other parameters were fixed. Parameters of the rock material are taken as those of the norite introduced in Section 3.1. Parameters of the filled medium for different water contents are taken according to Table 1.

3.2.1. Case of One Single Filled Joint

When a P wave obliquely impinges on a single filled joint, there are reflected P and SV waves as well as transmitted P and SV waves as shown in Figure 4.
Figure 5 shows reflection–transmission coefficients versus the filled thickness hf considering different water contents of the filled medium. It displays that Rp→p first increases abruptly with the increase in hf, then its increase tends to be gentle while Tp→p decreases basically in a negative exponential form. For the conversion waves, Rp→sv first increases sharply with hf increasing then it basically increases in a linear manner while Tp→sv basically keeps invariant. In general, Rp→sv varies more violently with the increase in hf compared with Tp→sv. As the water content increases, Tp→p decreases evenly, and corresponding curves become divergent; Rp→p and Rp→sv have an increasing trend on the whole; Tp→sv curves for different water contents basically overlap with each other.
Figure 6 shows the reflection –transmission coefficients varying with wave frequency f. We can see that Tp→p decreases basically in a negative exponential form with increasing f. Rp→p and Rp→sv first increase abruptly with increasing f, then they only have slight changes. Compared with the other three types of coefficients, Tp→sv is relatively insensitive to f. Cases of different water contents present similar variation of reflection–transmission coefficients versus f except for the magnitude difference. As the water content increases, Tp→p decreases evenly and corresponding curves become divergent; Both Rp→p and Rp→sv generally have an increasing trend of which the variation is smaller than that of Tp→p; Tp→sv curves for different water contents basically overlap with each other which indicates that it is independent of the water content of the filled medium.
Figure 7 shows the reflection–transmission coefficients versus the incident angle θ when a P wave obliquely impinges on a single filled joint considering the different water contents. It shows that Tp→p does not change much with increasing θ until θ approaches the critical angle 90°, and it suddenly decreases to 0. When θ is smaller than around 40°, Rp→p first decreases with increasing θ. As θ is larger than around 40°, it gradually increases to a peak and then it deceases instead. When θ is larger than around 83°, Rp→p increases rapidly. Tp→sv changes gently with the increase in θ. Rp→sv first increases to the maximum value at θ of approximately 40° and then it decreases with increasing θ instead. Generally, Tp→sv is much smaller than Tp→p. As the water content increases, Tp→p decreases evenly and related curves become divergent when θ is smaller than 83°; When the water content increases from 0% to 5%, both Rp→p and Rp→sv have an obvious increase while they have very limited changes as the water content increases from 5% to 10%; Tp→sv is basically independent on the water content of the filled medium.
Figure 8 shows the reflection –transmission coefficients versus the quality factor when a P wave impinges on a single filled joint with a thickness of 40 mm. For convenience, a variable Q is introduced to express Qp and Qs as Qp = Q and Qs = 2/3Q, respectively. To only investigate the influence of the quality factor of the filled medium, the density, thickness and wave velocity of the filled medium were fixed as constant values (Vpp = 185 m∙s−1, Vsp = 123 m∙s−1 and ρf = 1807 kg∙m−3). It shows that Tp→p decreases abruptly with Q increasing, after achieving a minimum value, it increases with an increment of Q instead and eventually it basically keeps invariant. Tp→sv does not change much with increasing Q. Both Rp→p and Rp→sv increase rapidly with Q increasing and then they increase gently, and they basically keep invariant as Q is larger than 20. Generally, Q has a more significant effect on Tp→p, Rp→p and Rp→sv than Tp→sv.
Due to the viscosity of the filled medium, the wave energy dissipated after propagating through the filled joint, which can be defined as follows [21,40]:
E l o s s = 1 T p p 2 R p p 2 T p s v 2 R p s v 2
Figure 9 shows the energy loss versus the quality factor Q corresponding to Figure 8. It displays that the energy loss Eloss first increases abruptly with the increase in Q. After achieving a peak, Eloss decreases rapidly with an increment of Q. When Q is larger than 25, Eloss changes gently with Q increasing.

3.2.2. Case of Multiple Parallel Filled Joints

When the stress wave propagates though a rock mass containing a set of joints, multiple reflections occur at each joint [33,34,43]. In this section, involved filled joints are assumed to have the same properties such as filled thickness, density, wave velocity, quality factor of the filled medium, etc. The rock mass is assumed to contain a set of joints with the identical joint spacing. A cycle of sinusoidal P wave with unit amplitude and frequency of 50 Hz was used to be the incident wave. The start time of the wave signal to oscillate has a delay of around 0.1 s. According to [21], the actual central frequency of the incident wave is 42.5 Hz for one cycle of sinusoidal wave with an input frequency of 50 Hz. It is reasonable to use the central frequency of 42.5 Hz to represent the frequency of the incident wave. The ratio between the joint spacing and wavelength have an important effect on the transmission coefficient because of the multiple reflections at filled joints [21,33,34,35]. Therefore, it is convenient to study the dimensionless joint spacing which is defined as ξ = hr/χ = (hrf)/Vpr, where χ is the P-wave wavelength.
Figure 10 shows transmitted waveforms (P-wave component) after a P wave obliquely propagates through four parallel filled joints with an identical filled thickness hf of 40 mm and joint spacing hr of 50 m, considering different water contents. It displays that each transmitted wave is made up of several cycles of waveforms arriving at different times induced by multiple reflections at joints. We can clearly see the wave superposition presenting multiple peaks and valleys in turn. As the water content increases, the transmitted wave decreases in the amplitude and the component of the later-arriving wave gradually fades away. When the water content is 15%, it is difficult to observe the later-arriving transmitted wave (after 0.18 s).
Figure 11 shows that variation of transmission coefficients Tp→p and Tp→sv versus ξ after the P wave obliquely propagates through four parallel joints considering different water contents. It shows that there exists two threshold lines ξ1 and ξ2 dividing each curve into three parts, which is similar to results in published papers [21,33,34,35,44]. Tp→p first decreases with increasing ξ (ξ < ξ1) and then it decreases instead (ξ1 < ξ < ξ2) and at last it basically keeps invariant (ξ > ξ2). With the water content increasing, Tp→p decreases and the position for the first threshold point moves to the left while the second one moves to the right. The increasing part of the curve (ξ < ξ1) shortens while the decreasing part (ξ1 < ξ< ξ2) widens. Tp→sv first increases with increasing ξ, then it decreases after achieving a peak in general. When the water content is low, there exist fluctuations on the curve. As the water content becomes higher, these fluctuations gradually fade away and the curve becomes more and more flat. From Figure 10 and Figure 11, it indicates that the effects of multiple reflections at joints on the transmitted wave and transmission coefficient were significantly reduced when the water content increased.

4. Discussion

It was found that not only the displacement but also the stress at the front and rear interfaces of the filled joint were discontinuous when loaded by the stress wave [20,45]. This attributes to the inertial effect of the mass of the filled medium. Both the stress and the displacement at the interface between the rock and filled material are assumed to be continuous in the current study. It can be shown that the sandwiched element surrounded by two interfaces as shown in Figure 1b can characterize the stress and displacement discontinuities at the filled joint according to Equation (34). Therefore, our model can simulate the discontinuous response of stress and the displacement at the filled joint, which has similar physics to DSDM by [20]. Additionally, our model can automatically consider the wave reflection inside the filled layer while DSDM loses the capability. However, the influence of the wave reflection is prominent and cannot be neglected when the ratio of the filled thickness to the wavelength is large enough. As described in Section 3.1, our analytical model was suitable to characterize the seismic response of the filled joint considering different water contents. Properties of the filled joints can be simply described by several common parameters such as the filled thickness, the density, the wave velocity, and the quality factors of the filled medium. These parameters have definite physical meanings and can be obtained relatively more readily compared with those of DSDM [20], i.e., the viscous stiffness.
The water content has a significant effect on the dynamic mechanical characteristics of the filled joint. As shown in Table 1, the increase in the water content can increase the density but reduce the wave velocity and quality factor of the filled medium. As a result, the transmission of the stress wave was significantly influenced. Especially, the variation of the quality factor stands for the variation of the viscosity. As derived in Equations (7) and (8), the viscoelastic wave equations have complex forms composed of a real part and imaginary part. The increase in the imaginary part has double effects, not only can it increase the equivalent dynamic modulus to let more wave energy pass through the filled joint (positive effect), but it can also cause energy loss (negative effect) [21]. For Tp→p, the negative effect first plays the domain role and the positive effect plays the domain role in turn when the quality factor increases. For Rp→p, Rp→sv and Tp→sv, the positive effect always plays the domain role. When the quality factor is large, the viscous effect of the filled medium is not prominent. In this case, the reflection–transmission coefficients are independent of the quality factor. We can now understand the variation of the reflection–transmission coefficients versus the quality factor as shown in Figure 8 as well as the energy loss with quality factor in Figure 9. The transmitted waves after propagating through multiple parallel joints are the superposed waves consisting of the first-arriving wave and the later-arriving waves caused by multiple reflections. When the water content increases, the wave reflection increases while the wave transmission decreases (see Figure 5, Figure 6 and Figure 7). According to this, the multiple reflections at the joint should increase with increasing water content. However, the multiple reflection decreases with increasing water content instead (see Figure 10 and Figure 11). It can be inferred that this phenomenon may result from the increase in the viscosity of the filled medium with the water content increasing, which significantly dissipated the energy of waves by multiple reflections at the filled joints. Thus, the change of the quality factor (viscosity) of the filled medium with increasing the water content has a significant influence on the seismic response of the filled joints which should be paid more attention in practical applications.

5. Conclusions

We established an analytical model to characterize stress wave propagation through rock joints filled with a viscoelastic medium considering different water contents. Some important conclusions were reached. The seismic response of the single filled joint highly depended on the water content, thickness of the filling as well as the frequency and incident angle of the incident wave, etc. The increase in the water content enhanced the viscosity (depicted by quality factor) of the filled joint, which could promote equivalent joint stiffness and energy dissipation with double effects on stress wave propagation. When the stress wave propagated through a set of filled joints, there existed multiple reflections at the joints. The dimensionless joint spacing was the main controlling factor in the seismic response of the multiple parallel filled joints, which exhibited three-stage variation. The transmission coefficient first increased with the dimensionless joint spacing, then it decreased instead, and at last it basically kept invariant. The increase in the water content could weaken the effect of multiple reflections. It mainly resulted from the increase in the viscosity of the filled joint which significantly dissipated the energy of waves by multiple reflections at the filled joints. The water content of the joint filling should be paid more attention in practical applications.

Author Contributions

Conceptualization, S.Q. and X.H.; methodology, X.H., S.Q., Y.L. and B.Z.; validation, H.X., B.Z., L.X., N.L.; formal analysis, S.Q., H.X., Y.L., B.Z.; investigation, B.Z., Y.L., L.X., N.L.; writing—original draft preparation, X.H.; writing—review and editing, S.Q., X.H.; supervision, S.Q.; project administration, S.Q.; funding acquisition, S.Q., X.H., L.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

This work was financially supported by the Second Tibetan Plateau Scientific Expedition and Research Program (STEP) under Grant No. 2019QZKK0904, the National Key Research and Development Plan of China under Grant No. 2019YFC1509701, National Natural Science Foundation of China under Grant Nos. 41807273, 41825018, 41941018, 41672307 and 41902289, Key Deployment Program of the Chinese Academy of Sciences under Grant No. KFZD-SW-422, Key Research Program of the Institute of Geology and Geophysics, CAS under Grant No. IGGCAS-201903402, Science Foundation of Key Laboratory of Shale Gas and Geoengineering, Institute of Geology and Geophysics, and Chinese Academy of Sciences under Grant No. KLSG201702.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Barton, N. A review of the shear strength of filled discontinuities in rock. Nor. Geotech. Inst. Publ. 1974, 105, 1–38. [Google Scholar]
  2. Jaeger, J.C.; Cook, N.G.; Zimmerman, R. Fundamentals of Rock Mechanics; John Wiley & Sons: Hoboken, NJ, USA, 2009. [Google Scholar]
  3. Cook, N.G.W. Natural Joints in Rock: Mechanical, Hydraulic and Seismic Behaviour and Properties Under Normal Stress. Int. J. Rock Mech. Min. Sci. Geomech. Abstr. 1992, 29, 198–223. [Google Scholar] [CrossRef]
  4. Goodman, R.E. Methods of Geological Engineering in Discontinuous Rocks; West Pubublishing Company: Eagan, MN, USA, 1976. [Google Scholar]
  5. Sun, G.Z. Rock Mass Structure Mechanics; Science Press: Beijing, China, 1988. (In Chinese) [Google Scholar]
  6. Guo, S.; Qi, S. Numerical study on progressive failure of hard rock samples with an unfilled undulate joint. Eng. Geol. 2015, 193, 173–182. [Google Scholar] [CrossRef]
  7. Zheng, B.; Qi, S. A new index to describe joint roughness coefficient (jrc) under cyclic shear. Eng. Geol. 2016, 212, 72–85. [Google Scholar] [CrossRef]
  8. Indraratna, B.; Oliveira, D.A.F.; Brown, E.T. A shear-displacement criterion for soil-infilled rock discontinuities. Géotechnique 2010, 60, 623–633. [Google Scholar] [CrossRef]
  9. Papaliangas, T.; Hencher, S.R.; Lumsden, A.C.; Manolopoulou, S. The effect of frictional fill thickness on the shear strength of rock discontinuities. Int. J. Rock Mech. Min. Sci. Geomech. Abstr. 1993, 30, 81–91. [Google Scholar] [CrossRef]
  10. Pereira, J.P. Rolling friction and shear behaviour of rock discontinuities filled with sand. Int. J. Rock Mech. Min. Sci. 1997, 34, 244.e1–244.e17. [Google Scholar] [CrossRef]
  11. Li, J.C.; Ma, G.W. Experimental study of stress wave propagation across a filled rock joint. Int. J. Rock Mech. Min. Sci. 2009, 46, 471–478. [Google Scholar] [CrossRef]
  12. Wu, W.; Zhao, J. Effect of water content on p-wave attenuation across a rock fracture filled with granular materials. Rock Mech. Rock Eng. 2015, 48, 867–871. [Google Scholar] [CrossRef]
  13. Huang, X.; Qi, S.; Yao, W.; Xia, K. Effect of Filling Humidity on the Propagation of High-Amplitude Stress Waves through an Artificial Joint. Geotech. Test. J. 2019, 42, 20170192. [Google Scholar] [CrossRef]
  14. Verruijt, A. An Introduction to Soil Dynamics; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2009; Volume 24. [Google Scholar]
  15. Das, B.; Luo, Z. Principles of Soil Dynamics; Nelson Education: Toronto, ON, Canada, 2016. [Google Scholar]
  16. Yoshioka, N.; Kikuchi, M. Visco-elastic response of joints to transmission waves. Geophys. Res. Lett. 1993, 20, 1143–1146. [Google Scholar] [CrossRef]
  17. Pyrak-Nolte, L.J.; Myer, L.R.; Cook, N.G.W. Transmission of seismic waves across single natural fractures. J. Geophys. Res. Space Phys. 1990, 95, 8617–8638. [Google Scholar] [CrossRef]
  18. Myer, L.; Pyrak-Nolte, L.; Cook, N. Effects of single fractures on seismic wave propagation. Rock Jt. 1990, 15, 467–473. [Google Scholar]
  19. Carcione, J.M.; Picotti, S. Reflection and transmission coefficients of a fracture in transversely isotropic media. Stud. Geophys. et Geod. 2011, 56, 307–322. [Google Scholar] [CrossRef]
  20. Zhu, J.; Perino, A.; Zhao, G.F.; Barla, G.; Li, J.C.; Ma, G.W.; Zhao, J. Seismic response of a single and a set of filled joints of viscoelastic deformational behaviour. Geophys. J. Int. 2011, 186, 1315–1330. [Google Scholar] [CrossRef] [Green Version]
  21. Huang, X.; Qi, S.; Liu, Y.; Zhan, Z. Stress wave propagation through viscous-elastic jointed rock masses using propagator matrix method (PMM). Geophys. J. Int. 2015, 200, 452–470. [Google Scholar] [CrossRef] [Green Version]
  22. Perino, A. Wave Propagation through Discontinuous Media in Rock Engineering. Ph.D. Thesis, Politecnico di Torino, Torino, Italy, 2011. [Google Scholar]
  23. Perino, A.; Orta, R.; Barla, G. Wave Propagation in Discontinuous Media by the Scattering Matrix Method. Rock Mech. Rock Eng. 2012, 45, 901–918. [Google Scholar] [CrossRef]
  24. Li, J.; Ma, G.; Huang, X. Analysis of Wave Propagation Through a Filled Rock Joint. Rock Mech. Rock Eng. 2009, 43, 789–798. [Google Scholar] [CrossRef]
  25. Sayers, C.; Taleghani, A.D.; Adachi, J. The effect of mineralization on the ratio of normal to tangential compliance of fractures. Geophys. Prospect. 2009, 57, 439–446. [Google Scholar] [CrossRef]
  26. Ma, G.W.; Li, J.C.; Zhao, J. Three-phase medium model for filled rock joint and interaction with stress waves. Int. J. Numer. Anal. Methods Geéomeéch. 2010, 35, 97–110. [Google Scholar] [CrossRef]
  27. Desai, C.S.; Zaman, M.M.; Lightner, J.G.; Siriwardane, H.J. Thin-layer element for interfaces and joints. Int. J. Numer. Anal. Methods Geéomeéch. 1984, 8, 19–43. [Google Scholar] [CrossRef]
  28. Sharma, K.G.; Desai, C.S. Analysis and Implementation of Thin? Layer Element for Interfaces and Joints. J. Eng. Mech. 1992, 118, 2442–2462. [Google Scholar] [CrossRef]
  29. Li, J.; Wu, W.; Li, H.; Zhu, J.; Zhao, J. A thin-layer interface model for wave propagation through filled rock joints. J. Appl. Geophys. 2013, 91, 31–38. [Google Scholar] [CrossRef] [Green Version]
  30. Zhu, J.B.; Zhao, X.B.; Wu, W.; Zhao, J. Wave propagation across rock joints filled with viscoelastic medium using modified recursive method. J. Appl. Geophys. 2012, 86, 82–87. [Google Scholar] [CrossRef] [Green Version]
  31. Rokhlin, S.; Wang, Y. Analysis of boundary conditions for elastic wave interaction with an interface between two solids. J. Acoust. Soc. Am. 1991, 89, 503–515. [Google Scholar] [CrossRef]
  32. Hopkins, D.; Myer, L.; Cook, N. Seismic wave attenuation across parallel fractures as a function of fracture stiffness and spacing. EOS Trans. AGU 1988, 69, 1427–1436. [Google Scholar]
  33. Cai, J.G.; Zhao, J. Effects of multiple parallel fractures on apparent attenuation of stress waves in rock masses. Int. J. Rock Mech. Min. Sci. 2000, 37, 661–682. [Google Scholar] [CrossRef]
  34. Zhao, J.; Zhao, X.; Cai, J.G. A further study of P-wave attenuation across parallel fractures with linear deformational behaviour. Int. J. Rock Mech. Min. Sci. 2006, 43, 776–788. [Google Scholar] [CrossRef]
  35. Li, J.; Ma, G.; Zhao, J. An equivalent viscoelastic model for rock mass with parallel joints. J. Geophys. Res. Space Phys. 2010, 115, 115. [Google Scholar] [CrossRef] [Green Version]
  36. Zhu, J.; Zhao, X.B.; Li, J.; Zhao, G.; Zhao, J. Normally incident wave propagation across a joint set with the virtual wave source method. J. Appl. Geophys. 2011, 73, 283–288. [Google Scholar] [CrossRef]
  37. Kennett, B. Seismic Wave Propagation in Stratified Media; ANU Press: Canberra, Australia, 2009. [Google Scholar]
  38. Aki, K.; Richards, P.G. Quantitative Seismology; University Science Books: Mill Valley, CA, USA, 2002; Volume 1. [Google Scholar]
  39. Zhao, X.B.; Zhu, J.; Zhao, J.; Cai, J.G. Study of wave attenuation across parallel fractures using propagator matrix method. Int. J. Numer. Anal. Methods Geéomeéch. 2011, 36, 1264–1279. [Google Scholar] [CrossRef]
  40. Carcione, J.M. Wave Fields in Real Media: Wave Propagation in Anisotropic, Anelastic, Porous and Electromagnetic Media; Elsevier: Amsterdam, The Netherlands, 2007; Volume 38. [Google Scholar]
  41. Bedford, A.; Drumheller, D. Elastic Wave Propagation; Wiley & Sons: Hoboken, NJ, USA, 1994. [Google Scholar]
  42. Wu, W.; Zhu, J.B.; Zhao, J. Dynamic response of a rock fracture filled with viscoelastic materials. Eng. Geol. 2013, 160, 1–7. [Google Scholar] [CrossRef] [Green Version]
  43. Huang, X.; Qi, S.; Guo, S.; Dong, W. Experimental Study of Ultrasonic Waves Propagating Through a Rock Mass with a Single Joint and Multiple Parallel Joints. Rock Mech. Rock Eng. 2013, 47, 549–559. [Google Scholar] [CrossRef]
  44. Zhao, X.B.; Zhao, J.; Hefny, A.M.; Cai, J.G. Normal Transmission of S-Wave Across Parallel Fractures with Coulomb Slip Behavior. J. Eng. Mech. 2006, 132, 641–650. [Google Scholar] [CrossRef]
  45. Wu, W.; Zhu, J.B.; Zhao, J. A further study on seismic response of a set of parallel rock fractures filled with viscoelastic materials. Geophys. J. Int. 2012, 192, 671–675. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (a) A P wave impinges on a rock mass containing joints filled with a viscoelastic medium, (b) the schematic view of the sandwiched element.
Figure 1. (a) A P wave impinges on a rock mass containing joints filled with a viscoelastic medium, (b) the schematic view of the sandwiched element.
Applsci 10 04797 g001
Figure 2. Variation of waveforms with the water content by experiments (reproduced with permission from [12], Springer Wien, 2015) and predicted results by our analytical model.
Figure 2. Variation of waveforms with the water content by experiments (reproduced with permission from [12], Springer Wien, 2015) and predicted results by our analytical model.
Applsci 10 04797 g002
Figure 3. Variation of the transmission coefficient versus the water content by experiments (reproduced with permission from [12], Springer Wien, 2015) and predicted by our analytical model.
Figure 3. Variation of the transmission coefficient versus the water content by experiments (reproduced with permission from [12], Springer Wien, 2015) and predicted by our analytical model.
Applsci 10 04797 g003
Figure 4. Reflected and transmitted waves when a P wave obliquely impinges on a single filled joint.
Figure 4. Reflected and transmitted waves when a P wave obliquely impinges on a single filled joint.
Applsci 10 04797 g004
Figure 5. Reflection–transmission coefficients of a P wave impinging on a single filled joint versus the filled thickness hf considering different water contents. The incident P wave has the frequency of 50 Hz and the incident angle of 20°.
Figure 5. Reflection–transmission coefficients of a P wave impinging on a single filled joint versus the filled thickness hf considering different water contents. The incident P wave has the frequency of 50 Hz and the incident angle of 20°.
Applsci 10 04797 g005
Figure 6. Variation of reflection–transmission coefficients versus the wave frequency f when a P wave impinges a filled joint with a thickness of 20 mm. The incident wave has the incident angle of 20° and the frequency of 50 Hz.
Figure 6. Variation of reflection–transmission coefficients versus the wave frequency f when a P wave impinges a filled joint with a thickness of 20 mm. The incident wave has the incident angle of 20° and the frequency of 50 Hz.
Applsci 10 04797 g006
Figure 7. Reflection–transmission coefficients versus the incident angle θ when the P wave impinges on a single filled joint with a thickness of 10 mm. The incident wave has the frequency of 50 Hz.
Figure 7. Reflection–transmission coefficients versus the incident angle θ when the P wave impinges on a single filled joint with a thickness of 10 mm. The incident wave has the frequency of 50 Hz.
Applsci 10 04797 g007
Figure 8. Reflection–transmission coefficients of a P wave impinging on a single filled joint versus the quality factor Q. The incident P wave has the frequency of 50 Hz and an incident angle of 20°.
Figure 8. Reflection–transmission coefficients of a P wave impinging on a single filled joint versus the quality factor Q. The incident P wave has the frequency of 50 Hz and an incident angle of 20°.
Applsci 10 04797 g008
Figure 9. The energy loss Eloss versus the quality factor Q corresponding to Figure 8.
Figure 9. The energy loss Eloss versus the quality factor Q corresponding to Figure 8.
Applsci 10 04797 g009
Figure 10. The incident P wave and transmitted waves (P-wave component) after propagating through four parallel filled joints with an identical spacing of 50 m considering different water contents. The incident P wave has the frequency of 50 Hz and the incident angle of 20°.
Figure 10. The incident P wave and transmitted waves (P-wave component) after propagating through four parallel filled joints with an identical spacing of 50 m considering different water contents. The incident P wave has the frequency of 50 Hz and the incident angle of 20°.
Applsci 10 04797 g010
Figure 11. Variation of the transmission coefficient (a) Tp→p and (b) Tp→sv versus dimensionless joint spacing ξ after the P wave propagates through four parallel joints considering different water contents. The incident P wave has a central frequency of 42.5 Hz and the incident angle of 20°.
Figure 11. Variation of the transmission coefficient (a) Tp→p and (b) Tp→sv versus dimensionless joint spacing ξ after the P wave propagates through four parallel joints considering different water contents. The incident P wave has a central frequency of 42.5 Hz and the incident angle of 20°.
Applsci 10 04797 g011
Table 1. Calibrated parameters of filled joints with different water contents.
Table 1. Calibrated parameters of filled joints with different water contents.
Parameters of the Filled MediumWater Content
0%2%5%10%15%
ρf (kg∙m−3)15721603165017291807
Vpp (m∙s−1)296270240220185
Vsp (m∙s−1)197180160146123
Qp3001203061.5
Qs200802041.0

Share and Cite

MDPI and ACS Style

Huang, X.; Qi, S.; Zheng, B.; Liu, Y.; Xue, L.; Liang, N. Stress Wave Propagation through Rock Joints Filled with Viscoelastic Medium Considering Different Water Contents. Appl. Sci. 2020, 10, 4797. https://doi.org/10.3390/app10144797

AMA Style

Huang X, Qi S, Zheng B, Liu Y, Xue L, Liang N. Stress Wave Propagation through Rock Joints Filled with Viscoelastic Medium Considering Different Water Contents. Applied Sciences. 2020; 10(14):4797. https://doi.org/10.3390/app10144797

Chicago/Turabian Style

Huang, Xiaolin, Shengwen Qi, Bowen Zheng, Youshan Liu, Lei Xue, and Ning Liang. 2020. "Stress Wave Propagation through Rock Joints Filled with Viscoelastic Medium Considering Different Water Contents" Applied Sciences 10, no. 14: 4797. https://doi.org/10.3390/app10144797

APA Style

Huang, X., Qi, S., Zheng, B., Liu, Y., Xue, L., & Liang, N. (2020). Stress Wave Propagation through Rock Joints Filled with Viscoelastic Medium Considering Different Water Contents. Applied Sciences, 10(14), 4797. https://doi.org/10.3390/app10144797

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop