Abstract
The region around the reservoir sometimes encounters the seismic response of ground. The phenomena of Love waves propagation also caused by seismic activity which can further significantly influenced by an impulsive line source present in the rock structures of ground. The current study elucidates the propagation behavior of Love waves induced by an impulsive line source within the anisotropic layered heterogeneous double porous rock structures commonly found in proximity to the reservoir areas. The expression of the dispersion relation for the Love waves propagation due to an impulsive line source has been derived by using the Fourier transformation and Green’s function technique. The above-mentioned rock composite structure is consisted of two different portions: upper transversely isotropic double porous (TIDP) layer medium and lower heterogeneous isotropic double porous half-space (IDP) medium. In the aforesaid half-space medium, both types of heterogeneities viz. linear and quadratic have been considered. The impacts of heterogeneity parameters associated with the lower heterogeneous IDP half-space, porosity parameter (corresponds to both upper TIDP layer and lower IDP half-space) and anisotropic parameter of TIDP layer on Love waves phase velocity have also been studied. Some special outlines have also been depicted graphically.
Similar content being viewed by others
Introduction
The study of seismic waves propagation caused by an earthquake phenomena are crucial to examine the interior of Earth as well as to analyze the seismic activity at the different margins of the Earth. The Earth consists of various kind of layers which are subjected to seismic waves due to the occurrence of earthquake. The propagation characteristic of seismic waves in layered structures is of primary importance due to its applications in geophysical prospective, civil engineering, mechanical engineering, etc. Ewing et al1. elaborated the different type of Earth structure according to the dispersive property of elastic layered media. Many investigations2,3have been done over the propagation behaviour of elastic waves in layered anisotropic media. A number of literature4,5provide the explanation on characterization of elastic waves in elastic layered media. Love waves are one of the devastating seismic surface waves produced by the earthquake. It may have an enormous effect on engineering designs, reservoirs and other man made structures6,7. In order to minimize the catastrophic effects of Love waves on human life and engineering structures, analysis of its propagation characteristics on different types of earth structures is crucial. Chattopadhyay et al8. analysed the propagation behaviour of Love waves in an isotropic elastic media lying over a semi-infinite elastic media and concluded a uniform decrease in the Love waves phase velocity with increasing non-dimensional wave number. Palermo et al9. investigated the control of Love wave dispersion through resonator parameter variation. Saha et al10. incorporated a numerical modelling of SH-wave propagation in initially-stressed multilayered composite structures.
In view of the problems with elastic stability for anisotropic media, the investigation of wave propagation in fluid saturated porous layer has been found to be very crucial. A porous medium is a type of material medium made up of the solid skeleton including with pore space that is often filled with fluid, is regarded as the poroelastic medium11. As the layer of the Earth is made up of various kind of rocks such as sandstone, lime stones, shale, etc., exhibiting the poroelastic nature. These poroelastic rocks are subjected to seismic waves during the occurrence of earthquake phenomena. Due to this, the investigation of the seismic waves propagation in poroelastic media is a matter of great interest for geoscientists and engineers. In view of above-mentioned reasons, many researchers made a lot of efforts to study the influence of properties of the poroelastic rock structures on the seismic wave’s propagation. The basic concepts regarding the waves propagation in poroelastic rock medium have been given by Biot12. Konczak13analysed the nature of Love wave propagation in the anisotropic porous layer. Also investigations have been done over the propagation of SH- wave in a poroelastic layer and on the dynamic interfacial crack effect in poroelastic strip14,15.
Generally the Earth medium are not homogeneous. Some kind of heterogeneity are found due to the variation of elastic properties in vertical depth. The heterogeneity can occur due to the phenomena of formation, growth, and coalescence of microcracks inside the rock medium16. The literature gives substantial evidence that the rock structures in the earth medium contains several sorts of heterogeneity, which can be expressed mathematically by different functions17,18. Based upon this idea, a lot of investigations have performed on the impact of various forms of heterogeneity viz. linear, quadratic, exponential or any other non-linear function of depths on the SH-wave’s propagation in an isotropic layered structure19,20. It has been examined by many studies that the Love wave’s propagation is strongly influenced by the heterogeneity of the structures. Kumari et al21. and Kundu et al22. investigated the effect of heterogeneity parameter on the Love wave’s propagation in a layered isotropic medium and heterogeneous micropolar media respectively. Gupta et al23. analysed the influence of propagation behaviour of Love wave on heterogeneity anisotropic porous media. Saha et al24. analyzed the impact of inhomogeneity on SH-type wave propagation in an initially stressed composite structure. Also, Sattari et al25. analyzed the body waves propagation through heterogeneous and discontinuous media. The problem having non-homogeneous boundary conditions may be homogenized with the use of some variable changes and will be able to provide solutions. Such type of non-homogeneous boundary value problem can be easily solved by using Green’s function methodology. Covert26was the researcher, who evolve a fresh method for finding green function for developed bodies. Also, Saha et al27. and Mario et al28. analyzed the behaviour of SH-waves propagation subjected to the curved boundary and imperfect boundary conditions respectively. Further, Venkatesan29 studied the Love-type wave propagation due to an impulsive point source by using the Green’s function technique.
Earthquakes are triggered by body forces, modeled as space-time dependent impulsive line sources, typically represented by the Dirac delta function. The resulting elastic displacement in layered structures is determined using the Green’s function technique, a proven method for solving elastodynamic problems involving impulsive sources. Several investigations30,31have been carried out over the propagation behaviour of Love waves utilizing Green’s function approach to handle impulsive point source problem in a heterogeneous medium, viscoelastic FGM half-space media respectively. Gupta et al32. analysed the effect of point source on propagation of SH-wave through an orthotropic crustal layer. Further, Gupta et al33. incorporated a comparative analysis due to the effect of point source on generation of SH wave. Hoop and Hijden34investigated the acoustic wave motion which depends on time-space are proliferated due to the impulsive line source. Also, Chattopadhyay et al35. investigated the propagation behaviour of SH-wave generated by impulsive line source in elastic media by considering Green’s function approach. Recently, Mahanty et al36. investigated the shear wave propagation in a heterogeneous poroelastic layered structure subjected to an impulsive line source. In many porous rocks, interconnected pores allow fluid to flow freely. However, fluid can sometimes become trapped, and prolonged entrapment can lead to mineral dissolution, creating patches within the rock. This process can form unique rock types where spherical inclusions are embedded within a porous host rock. As a result, two types of porosity emerge, storage porosity and transport porosity37. These rock formations, known as double porous structures, are commonly found near reservoir regions38. The formation of a large reservoir can cause the occurrence of an earthquake in the surrounding area39,40. Therefore, it is vital to investigate the behaviour of Love wave’s propagation in these double porous rock structures. Many investigations41,42have been carried out over the propagation behaviour of Love waves in various kind of double porosity structure. Recently, Gupta et al43. analyzed the propagation characteristic of Love wave in heterogeneous porosity media having discontinuity.
A critical review of the available literature reveals a gap in knowledge regarding the impact of impulsive line sources on Love wave propagation in heterogeneous double-porosity composite rock structures. Notably, no prior studies have investigated the propagation behavior of Love waves in a layered composite double-porosity structure consisting of a transversely isotropic double-porosity (TIDP) layer overlying a heterogeneous isotropic double-porosity (IDP) half-space influenced by an impulsive line source. It motivates the authors for performing the present investigation, which aims to formulate a problem that more closely reflects realistic scenarios by incorporating an impulsive line source in a double-porosity composite structure. This study assumes a non-linear function of vertically depths and stiffness to represent the material’s heterogeneity. To achieve an efficient analytical solution, the dispersion equation is derived in closed form by employing Green’s function technique in conjunction with Fourier transforms. The investigation is motivated with following research questions which aroused after carrying out a detail literature survey.
Key questions
-
What are the influence of heterogeneity parameter (linear as well as non linear quadratic heterogeneity) of the IDP half-space on the phase velocity of Love waves in the double porous rock structure?
-
How the porosity parameter (corresponds to both upper TIDP layer and lower IDP half-space) can effect the phase velocity of the Love waves in the double porous rock structure?
-
What are the effect of anisotropic parameter of TIDP layer on the phase velocity of Love waves in the double porous rock structure?
Mathematical model of dynamical problem
In the present study, we assume the Love wave’s propagation caused by an impulsive line source (S) in a double porous rock structure comprises an upper transversely isotropic double porous layer overlying the lower heterogeneous isotropic double porous half-space medium. Let H is the thickness of layer measured from the origin O. Here, the half-space medium is not homogeneous because of the consideration of varying stiffness and density. The heterogeneity in half-space can be expressed as the non-linear function of vertically downward depth. The Cartesian co-ordinate system containing (x, z) axes has been selected from this perspective that the waves propagation is considered in x-direction while z-axis is considered to be placed in the positive direction of vertical depth. The location of a line source S (impulsive in nature) is considered at a distance H below from the origin. It also represents the interface (\(z=H\)) between the upper layer and lower half-space medium as depicted in the Figure 1. The particle displacements in x, y, z direction are \(u_j\), \(v_j\) and \(w_j\) respectively, as for the layer j=1 and for the half-space medium j=2 respectively.
In view of the properties of Love waves in xz-plane causing the particle displacement in y-axis only, the displacement components for Love wave’s propagation can be represented as
where \(u_k^{(l)}, v_k^{(l)}, w_k^{(l)}\) and \(U_k^{(l)}, V_k^{(l)}, W_k^{(l)}\) are defined in Nomenclature respectively.
Dynamics of the upper transversely isotropic homogeneous double porous (TIDP) layer influenced by the line source S:
The constitutive equations governing the TIDP layer can be expressed as44
Considering equations (1) and (2), the non-vanishing equations of motion describing the Love wave’s propagation within a fluid-saturated TIDP layer can be stated as41
Now, introducing the function associated with the distribution of force density \(F_1(r,t)\) caused by an impulsive line source at the interface between TIDP layer and IDP half-space acting inside the medium. The existing equation of motion describing Love wave propagation in the TIDP layer caused by an impulsive line source derived from equations (1- 3) may be implemented as
where, r represents the distance between origin and coordinate points where the force is applied and
In view of harmonic nature of Love wave’s propagation, the following transformation can be taken
Using equation (6), equation (4) yields
We can now express the force density distribution function, \(F_1(r)\), caused by the impulsive line source disturbance by using the Dirac delta function.
Using equation (8) in equation (7), the only existing equation of motion impacted by the line source for the TIDP layer can be represented as
Considering the following Fourier transform of \(\bar{v}_1{(x,z)}\) as
and the related inverse Fourier transform as
where \(\Theta\) denotes the transformation parameter.
Using equation (10) in equation (9), we get
where \(4\pi F_1{(z)}\) = \(2\delta {(z-H)}/N_1\) and \(\alpha _1^2\)=\(\Theta ^2-{(\omega /\beta _1)^2}\) such that
denotes the velocity of Love waves corresponding to the TIDP layer and \(\gamma ^2\)=\({\frac{L_1}{N_1}}\) denotes the anisotropic parameter of TIDP layer.
Dynamics of isotropic heterogeneous double porous (IDP) half-space
The constitutive equations governing the IDP half-space is being expressed by45,
Considering equations (1) and (14), the non-vanishing equations of motion describing the Love wave’s propagation for IDP half-space can be represented as follows41
The presence of material heterogeneity can significantly impact the behaviour of the medium’s elastic constants and density. On account of, considering the heterogeneity in the form of linear as well as the quadratic function of zin vertical depths and stiffness as19
where \(R ={\epsilon _1(z-H) + \epsilon _2(z-H)^2}\).
Here, \(\epsilon _1\) and \(\epsilon _2\) denote the linear and quadratic heterogeneity parameters respectively. \(\epsilon _1\) and \(\epsilon _1\) are very small such that \(\epsilon _1<<1\) and \(\epsilon _1<<1\). \(N_{02}\) denotes the stiffness for the isotropic double porous half-space.
On account of equations (1), (14), (15) and (16), the only existing equation of motion for the Love waves propagation in the IDP half-space may be stated as
Now, using the following transformation (due to the harmonic nature of waves), we have
Using equation (18), equation (17) yields
On using the following Fourier transform for \(\bar{v}_2{(x,z)}\)
and the related inverse Fourier transform can be given as
Now, using the equation (20) in equation (19), we get
where \(\alpha _2^2\) = \(\Theta ^2-{(\omega /\beta _2)^2}\) in which
is represented as the velocity of the love waves corresponding the half-space medium, where
and
Boundary conditions
On considering the schematic diagram 1 of present problem and equations (2, 5, 10, 14, 16), the following possible boundary conditions can be provided
-
The upper transversely isotropic double porous layer has traction free upper surface, i.e.
$$\begin{aligned} \frac{dV_1}{dz}=0\,\,\,\ \text {at}\,\,\,\ z=0 \end{aligned}$$(25) -
At the common interface between the double porous layer and the half-space, there are continuous components for both stress and displacement. i.e, at \(z=H\),
$$\begin{aligned} L_1\frac{dV_1}{dz}&= N_{02}\frac{dV_2}{dz}\,\,\,\ \end{aligned}$$(26)$$\begin{aligned} V_1&=V_2\,\,\,\,\,\,\, \end{aligned}$$(27)
Solution of problem
We will now use the Green’s function technique for solving equations (12) and (22) with the aid of the necessary boundary conditions (25-27).
Let us introduce the Green’s function \(G_1{(z,z_0)}\) for any arbitrary \(z_0\) for the upper TIDP layer holding the boundary condition \(\frac{dG_1}{dz}\)=0 at \(z=0\) and at \(z=H\). Therefore, \(G_1{(z,z_0)}\) satisfies the equation
where \(z_0\) locates any where in the TIDP layer.
By multiplying equation (12) and equation (28) by \(G_1{(z,z_0)}\) and \(V_1(z_0)\) respectively, subtracting the resultant equations, and integrating within the range z=0 to \(z=H\), we get
Replacing \(z_0\) with z and utilizing Green’s function symmetry \(G_1(H,z) = G_1(z,H)\), equation (29) yields
Similarly, for IDP half-space, let \(G_2{(z,z_0)}\) be the Green’s function with fulfilling the following equations:
where \(z_0\) locates any where in the IDP half-space medium.
The similar procedure as discussed for the case of upper TIDP layer can be used for the lower IDP half-space. The equations (22) and (31) are multiplied by \(G_2{(z,z_0)}\) and \(V_2\) respectively; subtracting the resultant equations, and then integrating within the range z=H to \(z=\infty\), and using Green’s function property, we get
where \(T = 4\pi F_2(z_0)\).
On using equations (30), (32) and boundary conditions (26-27), we get
where \(A = \frac{1}{\gamma ^2}G_1{(H,H)}+\frac{L_1}{N_{02}} G_2{(H,H)}\).
Using the equations (33) in equation (30) and putting the value of \(4\pi F_2{(z_0)}\) from equation (24), we get
where
Furthermore, using the equations (33) and boundary condition(26), equation (32) leads to the following form
Now, putting the value of T in the above equation, we get
where
Now, \({V_2{(z)}}\) in equation (36) can be simplified by using the procedure of approximation. The obtained expression can be further used in the determination of \(V_1(z)\). In view of this, the first order approximation for \(V_2(z)\) can be derived by avoiding the higher powers of both parameters \(\epsilon _1\) and \(\epsilon _2\) present in equation (36) as
Equation (37) mathematically describes the displacement \(V_2(z)\) for any arbitrary position in the considered IDP half-space medium.
On account of equations (24), (34) and (37), the following expressions of \(V_1{(z)}\) can be obtained
In order to find the final value of \(V_1{(z)}\), it is needed to calculate the value of \(G_1\) and \(G_2\), which is accomplished by solving the following differential equation
Let \(S_1=e^{\alpha _1z}\) and \(S_2=e^{-\alpha _1z}\) represent the two distinct solutions to the differential equation (39) which vanish at \(z=\infty\) and \(z=-\infty\), respectively.
Therefore, the solution of equation (39) for an infinite medium can be taken as \({(e^{-\alpha _1|z-z_0|}/2\alpha _1)}\).
Consequently, \(S_1\) and \(S_2\)are combined to represent Green’s function for upper homogeneous TIDP layer as46
where P and Q are undetermined coefficients which can be found by following conditions
On using equation (41) into equation (40), we get
Using the similar approach, the value of \(G_2{(z,z_0)}\) can be calculated from equation (31), we get
Dispersion relation
The substitution of \(G_1(z,z_0)\) and \(G_2(z,z_0)\) from equations (42-43) into the equation (38) provides
Since, we have \(\epsilon _1<<1\) and \(\epsilon _2<<1\), therefore, the higher powers of \(\epsilon _1\) and \(\epsilon _2\) can be ignored. In view of this, the above equation (44) can be further expressed as
By applying the inverse Fourier transform mentioned in equation (11) into equation (45), we may further determine the displacement component at any position in TIDP layer as
The integral value of equation (46) is completely influenced by the integrand’s poles, which are situated at the equation’s roots
On replacing \(\Theta\) by k and \(\alpha _1\) by \(i \alpha _1^{\prime}\), equation (47) yields
where
Equation (48) characterizes the dispersion relation for Love waves propagating in a homogeneous TIDP layer overlying a heterogeneous IDP half-space.
Particular cases
Case (1)
When the TIDP layer becomes the single poroelastic layer and IDP half-space becomes the single poro-elastic half space i.e (\(\frac{(\rho _{02}^{2})^{(l)}}{\rho _{22}^{(l)}}\rightarrow 0\)) in the considered structure, the expression for dispersion relation in (48) reduces to
where
and
where (49) express the dispersion relation for love waves which propagates in a layered heterogeneous single poroelastic structure.
Case (2)
When the TIDP layer and IDP half-space in the considered layered double porous structure becomes the isotropic elastic layer i.e \(L_1\) = \(N_1\), (\(\frac{(\rho _{0k}^{2})^{(1)}}{\rho _{kk}^{(1)}}\rightarrow 0\)) and the isotropic elastic half-space (\(\frac{(\rho _{0k}^{2})^{(2)}}{\rho _{kk}^{(2)}}\rightarrow 0\)) respectively, then the dispersion relation (48) becomes
where
and
Equation (52) provides the dispersion relation for the Love wave’s propagation caused by the line source in the layered rock structure made up of isotropic elastic layer overlying the isotropic elastic half-space. The result obtained in (52) matches with the result (dispersion relation) obtained by Kumar et al19. for the case when the uppermost layer is removed \((H\,\rightarrow \,0\), H is the thickness of the uppermost layer) and the structure is comprised of only isotropic elastic layer overlying isotropic half-space.
Case (3)
When the considered layered double porous structure becomes the isotropic layered elastic structure i.e \(N_1=L_1\), (\(\frac{(\rho _{0k}^{2})^{(l)}}{\rho _{kk}^{(l)}}\rightarrow 0\)) and the isotropic elastic half-space is free from linear heterogeneity i.e \(\epsilon _1=0\), then the expression for dispersion relation (48) becomes
Equation (55) provides the dispersion relation for the Love wave’s propagation caused by the line source in the layered rock structure made up of transversely isotropic double porous layer lying over the isotropic double porous half-space having quadratic heterogeneity only. The result obtained in equation (55) matches with the result (dispersion relation) obtained by Chattopadhyay et al35. after ignoring viscosity.
Case (4)
When the considered layered double porous structure becomes the isotropic layered elastic structure and the isotropic elastic half-space i.e , \(N_1=L_1\), (\(\frac{(\rho _{0k}^{2})^{(l)}}{\rho _{kk}^{(l)}}\rightarrow 0\)) is free from quadratic heterogeneity i.e \(\epsilon _2=0\), then the expression for dispersion relation in (48) becomes
Equation (56) provides the dispersion relation for the Love wave’s propagation caused by the line source in the layered rock structure made up of isotropic elastic layer overlying the isotropic elastic half-space having linear heterogeneity only. The result obtained in equation (56) matches comparatively with the result (dispersion relation) obtained by Chattopadhyay et al47..
Case (5)
When the considered layered double porous structure becomes the isotropic layered elastic structure and the isotropic elastic homogeneous half-space i.e \(\epsilon _1=0\) and \(\epsilon _2=0\), (\(\frac{(\rho _{0k}^{2})^{(l)}}{\rho _{kk}^{(l)}}\rightarrow 0\)), then the expression for dispersion relation in (48) becomes
where
The expression obtained in equation (57) represents the propagation of Love wave in layered isotropic elastic structure made up of isotropic elastic upper layer overlying the isotropic elastic half-space. The expression obtained in equation (57) matches exactly with the standard Love wave’s equation1.
Numerical results and discussion
For the purpose of numerical simulations and graphical depiction, the data of the following physio mechanical properties have been considered:
For Upper Transversely Isotropic Double-Porous Layer44,45
For Lower Isotropic Double-Porous Half-Space44,45
Note: Aggregate density \(d^{(l)}\) = \(\rho _{00}^{(l)}\)+\(2\left( \rho _{01}^{(l)}+\rho _{02}^{(l)}\right)\) + \(\rho _{11}^{(l)}+ \rho _{22}^{(l)}\), as for the layer \(l=1\) and for the half-space \(l=2\) respectively.
Figure 2 represents the effect of linear heterogeneity parameter \((\nu _1=\frac{\epsilon _1 H}{N_{02}})\) of a IDP half-space on the phase velocity \((\frac{c_{ph}}{\beta _1})\) of Love waves in the layered double porous structures containing TIDP layer lying over IDP half-space. It is noted from Fig. 2 that as the linear heterogeneity parameter \((\nu _1)\) increases, the Love wave’s phase velocity decreases in the aforementioned structure.
Figure 3 shows the influence of quadratic heterogeneity parameter \((\nu _2=\frac{\epsilon _2 H^2}{N_{02}})\) of a IDP half space on the phase velocity \((\frac{c_{ph}}{\beta _1})\) of Love wave. It can be seen from the Fig. 3 that as the quadratic heterogeneity parameter \((\nu _2)\) of IDP half-space rises, the Love wave’s phase velocity in the structure under consideration falls.
It is observed from both figures that the heterogeneity parameters (both linear and quadratic) have adverse impact on the phase velocity of Love waves in the considered double porous structures. It may be due to the reason that as the heterogeneity (associated with the heterogeneity parameters) increases, the size of the grains of the rock materials varies. It may lead to the gapping at some positions which further may results into the decrease of phase velocity of Love waves.
Figure 4 demonstrate the impact of porosity parameter of TIDP layer \((d_l=\frac{\rho ^(1)}{d^(1)})\) on the phase velocity \((\frac{c_{ph}}{\beta _1})\) of Love wave. It’s been noticed from Fig. 4 that as the porosity parameter of the TIDP layer increases (i.e., porosity decreases), the Love wave’s phase velocity increases significantly in the studied double porous rock structure. Now, since most of the Love wave’s energy is centred on the finite layer, so, in the examined rock formation, the phase velocity of the Love waves increases as the upper layer becomes more compressed (i.e. porosity decreases).
Figure 5 explains the effectiveness of porosity parameter of a IDP half-space \((d_h=\frac{\rho ^(2)}{d^(2)})\) on the phase velocity \((\frac{c_{ph}}{\beta _1})\) of Love waves. It has been noticed from Fig. 5 that as the porosity parameter of the IDP half-space increases (i.e., porosity decreases), a significant reduction has occurred in the Love wave’s phase velocity. Probably because, due to the heterogeneity in the half-space, a feeble compaction has been found comparatively in the lower half-space. Also, as the energy of Love waves centred on the finite layer, only a minor portion of it propagates in the lower half-space which could be the reason behind the drop in phase velocity of Love waves.
Figure 6 exhibits the impact of anisotropic parameter (\(\gamma\)) of TIDP layer on the phase velocity \((\frac{c_{ph}}{\beta _1})\) of Love waves. It has been revealed From the Fig. 6 that as the anisotropic parameter (\(\gamma\)) increases, the Love wave’s phase velocity decreases considerably in the regarded double-porous rock structure, it happens because as the ratio of anisotropic parameter (\(L_1\)/\(N_1\)) associated with the TIDP layer increases, the energy of Love waves gets dispersed in many direction, causing the Love wave’s phase velocity reduces.
Figure 7 exhibits a comparative study of the effect to the phase velocity \((\frac{c_{ph}}{\beta _1})\) of Love wave in different kind of scenario. The figure demonstrates that the phase velocity is maximum for case (5) i.e for the case of isotropic layered elastic structure and the isotropic elastic homogeneous half-space and minimum for case (3) i.e for the case of isotropic layered elastic structure and the isotropic elastic half-space is free from linear heterogeneity.
Conclusion
This study investigates Love wave propagation in a heterogeneous double-porosity structure influenced by an impulsive line source. The considered structure comprises a homogeneous transversely isotropic double porous (TIDP) layer overlying a heterogeneous isotropic double porous (IDP) half-space with depth-dependent stiffness. Green’s function technique yields the Love wave dispersion relation, and the influence of heterogeneity, porosity, and anisotropy on phase velocity is analyzed. Numerical simulations and graphical representations are included. The single-porosity case is derived as a special instance. Validation is achieved through comparisons with existing literature. Key findings are presented subsequently.
-
1.
The increase in the linear heterogeneity parameter of the IDP half-space causes reduction in phase velocity of Love wave’s in the double porous rock structure.
-
2.
The IDP half-space’s non-linear quadratic heterogeneity parameter adversely affects Love wave phase velocity in the double porous rock structure.
-
3.
The augmentation in the porosity parameter (i.e decrease in the porosity) of the TIDP layer causes the upsurge in the phase velocity of Love waves in the studied double porous structure. On the contrary, the increment in the porosity parameter of the IDP half space significantly resist the Love wave’s phase velocity.
-
4.
The phase velocity of Love wave’s propagation decreases with the enhancement of anisotropic parameter (associated with the extent of anisotropy) of the TIDP layer in the double porous rock structure.
Engineering applications
The regions near the reservoir sometimes encounter the seismic activity. The impulsive line source in the rock materials can significantly affect the behavior of seismic wave propagation. The results obtained from the current study can have significant possible applications to examine the dynamic strength of the rock structure that can be found near the reservoirs. From the current analysis, the medium parameters that are responsible for the rise and fall of phase velocity can be examined. It can be useful to select the appropriate rocks for the construction of more safe engineering structure in the earthquake prone area nearby pools, pond, and reservoirs regions.
It has been examined from the obtained results that higher the porosity parameter (i.e. lower the porosity) of upper TIDP layer causes the increment in the phase velocity of Love waves, while an increase in the anisotropy of the upper TIDP layer is associated with a decrease in the phase velocity of Love waves. On the other side, if the porosity is less and the heterogeneity is more in the lower half-space medium then the Love wave’s velocity reduces in the considered double porous structure. In view of above, the double porous rock blocks for the construction purpose have been selected in such a way that its upper portion is of higher porosity and higher anisotropy while the lower bigger portion is of lower porosity and higher heterogeneity. Due to this, the velocity of Love waves (i.e the damaging effect of Love waves) can be reduced. In this way, the engineering structures influenced by Love waves near the reservoir region can be constructed with the more efficient and safe design.
Data availability
All data relevant to this study has been provided with in the manuscript.
Abbreviations
- \(\tau _{ij}^{(l)}\) :
-
Incremental stresses induced in solid phase of TIDP layer (l1) and IDP half-space (l2).
- \(N_1, L_1\) :
-
Modulus of rigidity in the direction of the symmetric axis and in the upper layer transversely isotropic plane, respectively.
- \(\phi _l^k\) :
-
Host medium (k1) porosity and inclusion spherical shape (k2) porosity for TIDP layer (l1) and IDP half-space (l2).
- \(p_l^{f_1}\) and \(p_l^{f_2}\) :
-
Fluid pressure corresponds to the host medium (\(f_{1}\)) and fluid pressure is related with the spherical inclusions (\(f_{2}\)) in TIDP layer (l1) and IDP half-space (l2), respectively.
- \(\lambda _1, \mu _1, I_1, S_1, S_2, S_3, S_4, \zeta _1, k_1\) and \(k_2\) :
-
Elastic constants for the TIDP layer.
- \(\rho _{00}^{(l)}\) :
-
Inertial coefficient of TIDP layer’s (l1) and IDP half-space (l2) for solid skeleton.
- \(\rho _{0k}^{(l)}\) :
-
Inertial coefficient associated with the coupling between the host medium’s solid and liquid phases (k1) and for the inclusions spherical shapes (k2) in the TIDP layer (l1) and IDP half-space (l2), respectively.
- \(\rho _{kk}^{(l)}\) :
-
Inertial coefficient associated with the fluid contained in the host medium (k1) and the inertial coefficient related with the fluid contained in the spherical inclusions (k2) corresponds to TIDP layer’s (l1) and the IDP half-space’s (l2), respectively.
- \(u_k^{(l)}, v_k^{(l)}, w_k^{(l)}\) :
-
Solid phase particle displacement for the host medium (l1) as well as for the spherical inclusion (l2) for both the TIDP layer (k1) and the IDP half-space (k2) respectively.
- \(U_k^{(l)}, V_k^{(l)}, W_k^{(l)}\) :
-
Fluid phase particle displacement for the host medium (l1) as well as for the spherical inclusion (l2) for both the TIDP layer (k1) and the IDP half-space (k2) respectively.
- \(\lambda _2, N_2\) :
-
Modulus of rigidity for the IDP half-space.
- \(S_5, S_6, k_3\) and \(k_4\) :
-
IDP half space elastic constants.
References
Ewing, W. M., Jardetzky, W. S., Press, F. & Beiser, A. Elastic waves in layered media. Physics Today 10(12), 27–28 (1957).
Anderson, D. L. Elastic wave propagation in layered anisotropic media. Journal of Geophysical Research 66(9), 2953–2963 (1961).
Alam, P., Kundu, S. & Gupta, S. Effect of magneto-elasticity, hydrostatic stress and gravity on Rayleigh waves in a hydrostatic stressed magneto-elastic crystalline medium over a gravitating half-space with sliding contact. Mechanics Research Communications 89, 11–17 (2018).
Ostachowicz, W., Kudela, P., Krawczuk, M. & Zak, A. Introduction to the theory of elastic waves, Guided Waves in Structures for SHM, Chichester, West Sussex, United Kingdom, John Wiley & Sons Ltd, 1–46 (2012).
Gupta, S. et al. Dynamic analysis of wave propagation and buckling phenomena in carbon nanotubes (cnts). Wave Motion 104, 102730 (2021).
Eringen, A. & Suhubi, E. Elastodynamics Vol. ii (Academic), New York, 1975).
Achenbach, J. Wave propagation in elastic solids 1st edn. (Elsevier, New-York, 2012).
Chattopadhyay, A., Singh, P., Kumar, P. & Singh, A. K. Study of Love-type wave propagation in an isotropic tri layers elastic medium overlying a semi-infinite elastic medium structure. Waves in Random and Complex Media 28(4), 643–669 (2018).
Palermo, A. & Marzani, A. Control of Love waves by resonant metasurfaces. Scientific reports 8(1), 7234 (2018).
Saha, S., Chattopadhyay, A. & Singh, A. K. Numerical modelling of SH-wave propagation in initially-stressed multilayered composite structures: A casewise study. Engineering Computations 36(1), 271–306 (2019).
Kumar, P., Singh, A. K. & Chattopadhyay, A. Influence of an impulsive source on shear wave propagation in a mounted porous layer over a foundation with dry sandy elastic stratum and functionally graded substrate under initial stress. Soil Dynamics and Earthquake Engineering 142, 106536 (2021).
Biot, M. A. Theory of propagation of elastic waves in a fluid-saturated porous solid. ii. higher frequency range, The Journal of the acoustical Society of america 28 (2) 179–191 (1956).
Kończak, Z. The propagation of Love waves in a fluid-saturated porous anisotropic layer. Acta Mechanica 79(3), 155–168 (1989).
Son, M. S. & Kang, Y. J. Propagation of shear waves in a poroelastic layer constrained between two elastic layers. Applied Mathematical Modelling 36(8), 3685–3695 (2012).
Negi, A., Singh, A. & Yadav, R. Analysis on dynamic interfacial crack impacted by SH-wave in bi-material poroelastic strip. Composite Structures 233, 111639 (2020).
Sato, H., Fehler, M. C. & Maeda, T. Seismic Wave Propagation and Scattering in the Heterogeneous Earth: Second Edition, Springer, (2012).
Meissner, R. The little book of planet Earth (Springer, New-York, 2002).
Mahanty, M., Chattopadhyay, A., Kumar, P. & Singh, A. K. Effect of initial stress, heterogeneity and anisotropy on the propagation of seismic surface waves. Mechanics of Advanced Materials and Structures 27(3), 177–188 (2020).
Kumar, P., Mahanty, M., Chattopadhyay, A. & Kumar Singh, A. Green’s function technique to study the influence of heterogeneity on horizontally polarised shear-wave propagation due to a line source in composite layered structure, Journal of Vibration and Control 26 (9-10) 701–712 (2020).
Arif, M. & Alam, P. SH-wave interaction with MTI strip coated with Newtonian viscous liquid laid over self-weighted heterogeneous half-space: a comparative study, Mechanics of Solids 1–27 (2024).
Kumari, N., Anand Sahu, S., Chattopadhyay, A. & Kumar Singh, A. Influence of heterogeneity on the propagation behavior of Love-type waves in a layered isotropic structure, International Journal of Geomechanics 16 (2) 04015062 (2016).
Kundu, S., Kumari, A., Pandit, D. K. & Gupta, S. Love wave propagation in heterogeneous micropolar media. Mechanics Research Communications 83, 6–11 (2017).
Gupta, S., Pramanik, A. & Ahmed, M. Impact of pre-stress, inhomogeneity and porosity on the propagation of Love wave. Acta Geophysica 66, 855–866 (2018).
Saha, S., Chattopadhyay, A. & Singh, A. Impact of inhomogeneity on SH-type wave propagation in an initially stressed composite structure. Acta Geophysica 66, 1–19 (2018).
Sattari, A. S., Rizvi, Z. H., Aji, H. D. & Wuttke, F. Study of wave propagation in discontinuous and heterogeneous media with the dynamic lattice method. Scientific reports 12(1), 6343 (2022).
Covert, E. E. A note on an approximate calculation of green’s functions for built-up bodies. Journal of Mathematics and Physics 37(1–4), 58–65 (1958).
Saha, S., Singh, A. K. & Chattopadhyay, A. Impact of curved boundary on the propagation characteristics of Rayleigh-type wave and SH-wave in a prestressed monoclinic media. Mechanics of Advanced Materials and Structures 28(12), 1274–1287 (2021).
Mario A, J. S. & Alam, P. Impacts on SH-waves regulating through a FGPM plate clamped between a temperature dependent plate and a microstructural coupled stressed substrate subjected to the perfect and imperfect boundary conditions, Journal of Vibration Engineering & Technologies 1–15 (2024).
Venkatesan, P. & Alam, P. A multi-layered model of poroelastic, HSTI, and inhomogeneous media to study the Love-type wave propagation due to an impulsive point source: a Green?s function approach. Acta Mechanica 235(1), 409–428 (2024).
Kundu, S., Gupta, S., Vaishnav, P. K. & Manna, S. Propagation of Love waves in a heterogeneous medium over an inhomogeneous half-space under the effect of point source. Journal of Vibration and Control 22(5), 1380–1391 (2016).
Kundu, S., Kumhar, R., Maity, M. & Gupta, S. Influence of point source on love-type waves in anisotropic layer overlying viscoelastic FGM half-space: Green?s function approach. International Journal of Geomechanics 20(1), 04019141 (2020).
Gupta, S. et al. Inquisitive analysis of the point source effect on propagation of SH wave through an orthotropic crustal layer. Journal of Solid Mechanics 10(4), 831–844 (2018).
Gupta, S., Pramanik, S., Smita, & Verma, A. K. A comparative analysis due to the effect of point source on generation of SH wave, Journal of Earth System Science 128, 1–18 (2019).
de Hoop, A. T. & Van der Hijden, J. H. Generation of acoustic waves by an impulsive line source in a fluid/solid configuration with a plane boundary. The Journal of the Acoustical Society of America 74(1), 333–342 (1983).
Chattopadhyay, A., Gupta, S., Kumari, P. & Sharma, V. K. Effect of point source and heterogeneity on the propagation of SH-waves in a viscoelastic layer over a viscoelastic half space. Acta Geophysica 60, 119–139 (2012).
Mahanty, M., Kumar, P., Singh, A. K. & Chattopadhyay, A. Green?s function analysis of shear wave propagation in heterogeneous poroelastic sandwiched layer influenced by an impulsive source. Wave Motion 107, 102821 (2021).
Olny, X. & Boutin, C. Acoustic wave propagation in double porosity media. The Journal of the Acoustical Society of America 114(1), 73–89 (2003).
Choquet, C. Derivation of the double porosity model of a compressible miscible displacement in naturally fractured reservoirs. Applicable Analysis 83(5), 477–499 (2004).
Simpson, D. Triggered earthquakes. Annual Review of Earth and Planetary Sciences 14(1), 21–42 (1986).
Gupta, H. K. A review of recent studies of triggered earthquakes by artificial water reservoirs with special emphasis on earthquakes in Koyna. India, Earth-Science Reviews 58(3–4), 279–310 (2002).
Pal, M. K., Negi, A. & Singh, A. K. Propagation of Love-type wave in an imperfectly bonded double-porous composite rock structure impacted by liquid loading. International Journal of Geomechanics 22(12), 06022033 (2022).
Rajak, B. P., Kundu, S., Gupta, S. & Kumar, D. Love wave propagation characteristics in a fluid-saturated cracked double porous layered structure. Mechanics of Advanced Materials and Structures 31(11), 2349–2361 (2024).
Gupta, S., Dutta, R. & Das, S. Love-type wave propagation in an inhomogeneous cracked porous medium loaded by heterogeneous viscous liquid layer. Journal of Vibration Engineering & Technologies 9, 433–448 (2021).
Sharma, M. Effect of local fluid flow on the propagation of elastic waves in a transversely isotropic double-porosity medium. Geophysical Journal International 200(3), 1423–1435 (2015).
Ba, J., Carcione, J. & Nie, J. Biot-Rayleigh theory of wave propagation in double-porosity media, Journal of Geophysical Research: Solid Earth 116 (B6) (2011).
Chattopadhyay, A., Chakraborty, M. & Kushwaha, V. On the dispersion equation of Love waves in a porous layer. Acta mechanica 58(3), 125–136 (1986).
Chattopadhyay, A., Pal, A. & Chakraborty, M. SH waves due to a point source in an inhomogeneous medium. International journal of non-linear mechanics 19(1), 53–60 (1984).
Acknowledgements
The authors express their sincere thanks to VIT-AP University, Andhra Pradesh to provide the research facility to carried out this research work.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of Interest
The authors confirm that there is no potential conflict of interest with this study.
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
About this article
Cite this article
Mishra, A., Negi, A. Investigation on the transmission of Love waves due to an impulsive line source in a heterogeneous double porous rock structure. Sci Rep 14, 27499 (2024). https://doi.org/10.1038/s41598-024-77203-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-024-77203-1
Keywords
This article is cited by
-
Love Wave Propagation at a Multi-fluid Porous Layer and Initially Stressed Elastic Half-Space with Triangular Irregularity
International Journal of Applied and Computational Mathematics (2025)