Reflection/refraction at the interface of an elastic solid and a partially saturated porous solid containing liquid filled bound pores and a connected pore space saturated by two-phase fluid

The phenomena of reflection and refraction of plane waves incident obliquely at a plane interface between uniform elastic solid half-space and porous solid containing liquid filled bound pores and two-phase fluid in connected pores has been analyzed. The amplitude ratios of the reflected and refracted waves to that of the incident wave are calculated as a non-singular system of linear algebraic equations. These amplitude ratios are used further to derive the expressions for the partition of incident energy among the reflected and refracted waves. Partition of incident energy among the reflected and refracted waves is studied for incidence of P and SV waves. The conservation of the energy across the interface is verified. The effect of gas saturation, wave frequency, capillary pressure and bound liquid film on the amplitude ratios and energy partitions are studied in the numerical example.


INTRODUCTION
A poroelastic solid is considered as an elastic matrix with a Newtonian fluid filling its connected pores.Dynamic behavior of fluids saturated porous media has been the center of study due to their importance in modelling the sedimentary materials in the field of acoustics, oil exploration, earthquake engineering, soil dynamics and hydrology.The dynamic equations formulated by Biot (1956aBiot ( , b, 1962) ) are, generally, used to derive the mathematical models for wave propagation studies in a poroelastic solid saturated completely with a single-fluid phase.Pride et al. (2004) extended the Biot's single pore fluid formulation to the porous media saturated by multiple fluids.
Mixture theory seems to be convenient in studying the wave propagation in a porous solid saturated by multiphase fluid.It was Brutsaert (1964), who predicted the existence of third compressional wave due to presence of second fluid in pores.In this, the constituent phases are assumed to exist everywhere, but without any integration between them.Bowen (1976Bowen ( , 1980) ) has given an extensive study of mixture theory.The work of Bedford and Drumheller (1983) deals with mixture of immiscible constituents.Garg and Nayfeh (1986) have studied the propagation of compressional waves in porous media saturated with chemically non-reactive miscible liquid/gas mixture.Santos et al. (1990a, b) derived the governing equations and presented a method to calculate elastic constants for isotropic porous solids saturated by two-phase fluids.Garg and Nayfeh (1986), and Tuncay and Corapcioglu (1997) formulated a comprehensive procedure relevant to wave propagation in porous solids saturated with multiple fluids.A recent mathematical model presented by Lo et al. (2005) is also based on continuum mixture theory.It is general enough to account for changes in capillary pressure and viscous/inertial coupling among the constituents.In the absence of inertial coupling, this model reduces to that of Tuncay and Corapcioglu (1997).Most of studies on wave propagation in porous media prefer to use the elastodynamics of Biot's theories and the boundary conditions of Deresiewicz and Skalak (1963), for example, Pride et al. (1992), Kaynia and Banerjee (1993), Gurevich and Schoenberg (1999) and Denneman et al. (2002).Burridge and Keller (1981) used two space method of homogenization to derive the constitutive equations for poroelasticity from microstructure.
In the contemporary times, acoustic wave propagation in porous media has also got importance in the oil explorations and medical field.During the propagation of a sound wave in such a medium, interactions between these two phases of different nature take place, giving various physical properties that are unusual in classical media.The large contact area between solid and fluid, which is the main characteristic of porous media induces new phenomena of diffusion and transport in the fluid, in relation to micro-geometry of the pore space.Many applications are concerned with understanding the behavior of acoustic waves in such media.In geophysics, we are interested in the propagation of acoustic waves in porous rocks, for information on soil composition and their fluid content.Acoustic characterization of materials is often achieved by measuring the attenuation coefficient and phase velocity in the frequency domain.The most recent theoretical and experimental methods developed by the authors for the acoustic characterization of porous materials are shown in Fellah et al. (2013).
Besides, seismic waves in the earth's crust are influenced by the properties of the strata through which they travel and the reflection/refraction at discontinuities between different layers.Recorded signals of these waves provide information about the internal structures of the Earth which is used further in devising an effective strategy for the exploration of minerals and hydrocarbons.Reflection and refraction of elastic waves at the boundaries of fluid saturated porous materials is a process, which has direct relevance in the studies on geophysical exploration.The problem of reflection and refraction of plane elastic waves striking at the plane interface between an elastic solid and a poroelastic solid saturated by a single fluid/two immiscible fluids have been attempted by many researchers.The latest book by Carcione (2007) is referred to for relevant references and detailed procedures.
In recent years, some studies (Tomar and Arora, 2006;Arora andTomar, 2007, 2010;Yeh et al., 2010) have considered the reflection and refraction of plane harmonic waves at the boundaries of porous media saturated by two immiscible fluids.In this study, pore-fluids were assumed non-viscous Latin American Journal of Solids and Structures 12 (2015) 1870-1900 so as to avoid the involvement of attenuation.For the same reason, the incidence was restricted to pre-critical angles.Unfortunately, this is in contrast to the realistic flow mechanics in crustal rocks where the equilibration of fluid pressure produces a great deal of seismic attenuation (Sams et al., 1997).However, Sharma and Kumar (2011) and Kumar and Saini (2012) ignored all these restrictions.Sharma and Saini (2012) studied the wave propagation in porous solid containing liquid filled bound pores and two-phase fluid in connected pores.Kumar and Sharma (2013) studied the reflection and transmission of attenuated waves at the boundary between two dissimilar poroelastic solids saturated with two immiscible viscous fluids.Kumar and Kumari (2014) studied the reflection of attenuated waves at the surface of fractured porous solids saturated with two immiscible viscous fluids.
The present work generalizes the reflection at the free surface studied by Sharma and Saini (2012) to the reflection-refraction phenomenon at the plane interface between uniform elastic solid and porous solid saturated with two miscible/immiscible fluids containing liquid filled bound pores.A porous medium is considered to be dissipative due to the presence of the viscosity in the pore fluids.Hence, the waves refracted to the dissipative porous medium are identified as inhomogeneous waves with the attenuation always normal to the interface.An energy matrix is calculated, which defines the shares of two waves reflected to elastic solid and four waves refracted to saturated porous solid containing liquid filled bound pores.This matrix enables to identify the interaction energy among the refracted waves, which is required to ensure conservation of energy at the interface.Numerical example is considered to study the nature of dependence of amplitude ratios and energy ratios on angle of incidence of the incident wave.The conservation of the energy across the interface is verified.The effects of gas saturation, wave frequency, capillary pressure and bound liquid film on the amplitude ratios and energy partitions are depicted graphically and discussed.

BASIC EQUATIONS
The composite porous medium consists of four constituents, i.e., solid grains, bound liquid film, poreliquid, pore-gas, which are identified with indices ' ' s , ' '  , ' ' l , ' ' g respectively.Out of the total porosity   f of the medium, a fraction α is occupied by bound liquid film and the remaining part is the connected porosity  .Then, the volume fractions of the constituents are defined as where  is the fraction of gas saturation in connected pore-space.These volume fractions are scaling functions which are used to relate partial and intrinsic values of any characteristic of the medium.For example, the product s  defines the contribution of solid grains in the aggregate density    of multiphase mixture.
In applying continuum models to treat multiphase media, it is assumed that the local variables can be replaced by mixture variables averaged over a region, which is quite large in comparison to grain-size but very small when compared to sample-size.It is further assumed that for each of the phases, i) partial stress tensors are symmetric, ii) external body forces are absent, and iii) deformations are infinitesimal.The gas is assumed to be soluble in liquid but no mass exchange is allowed between the solid matrix and the twin-phase pore-fluid.Inertial coupling between the mixture constituents is Latin American Journal of Solids and Structures 12 (2015) 1870-1900 excluded.Following Garg and Nayfeh (1986), the equations of motion for the low-frequency vibrations of constituent particles in isotropic porous solid are given by where the superscript ' ' d is used to denote drained porous solid frame.The particles of elastic skeleton and bound liquid have the same displacement and pressure.Hence, the drained porous matrix is considered a single continuum which behaves viscoelastic to wave propagation (Edelman, 1997).
'  s are used to define stresses and '  s are intrinsic densities.i u , i v and i w denote the components of displacements of the drained solid, liquid and gas particles, respectively.The indices (other than s ,  , l , g ) can take values 1, 2 and 3. A repetition of these indices implies summation.Dot over a variable implies partial derivative with time and comma before an index implies partial space differentiation.
Darcy's law relates viscous dissipation to the motion of gas and liquid particles relative to the pore-walls.The assumption of Poiseuille flow, necessary for this law, breaks down if the frequency exceeds a certain value.The present work is specifically restricted to low frequency such that viscousfluid dissipation does not depend on frequency.Following Garg and Nayfeh (1986), dissipation coefficients for liquid ( l d ) and gas ( g d ) are defined as follows: where k  and k  define the viscosity and the relative permeability of fluid phase k.  denotes the intrinsic permeability of the porous medium.Garg and Nayfeh (1986) employed the concepts from the theory of interacting continua (Bedford and Drumheller, 1983) to formulate separate constitutive models for different constituents of porous medium.In terms of intrinsic stress tensors and densities, the constitutive relations for porous matrix, liquid and gas are defined as follows: ) 3 where ij  is Kronecker symbol.The elastic constants '  s derived from the elastic moduli of the constituents, are given in the Appendix.
In Kelvin-Voigt model of linear viscoelasticity, the elastic solid element and viscous fluid element are assumed to be parallel and thus subjected to same strain.Then, an effective elastic modulus of Latin American Journal of Solids and Structures 12 (2015) 1870-1900 the composite is obtained as the sum of partial values of the modulus for different constituents (Wong and Bollampally, 1999).Following Edelman (1997), the time-dependent rigidity modulus p  of viscoelastic porous frame relates to the rigidity ( s  ) of solid grains as follows: where   is the dynamic shear viscosity and Re  is acoustic Reynolds number for bound liquid film.Edelman (1997) has defined this number as 2 1 Re    .Non-dimensional parameter  is expressed in terms of fluid viscosity, bulk modulus, density and medium permeability (Nikolaevskiy, 1990;Maksimov et al., 1994).

Saturated porous solid
Propagation of harmonic plane waves is considered in a partially saturated porous solid containing liquid filled bound pores and a connected pore space saturated by two-phase fluid, identified as gas and liquid.Following Sharma and Saini (2012), in a partially saturated porous solid containing liquid filled bound pores and a connected pore space saturated by two-phase fluid, three dilatational waves and one shear wave propagate.For convenience in discussion, the three longitudinal waves with are named as I P , II P , III P waves, respectively.The lone transverse wave is identified as S wave.For the displacements   , , j j j u v w , the system of Christoffel equations in Sharma and Saini ( 2012) is solved to define the complex velocities ( j V , 1, 2, 3, 4 j  ) of four attenuated waves in the medium.Corresponding to each wave, the polarizations ( S , L , G ) for material particles are calculated to define the displacements of material particles as given in Sharma and Saini (2012).

Elastic solid
Equation of motion for isotropic elastic solid is given by The isotropic stress-strain relation in the elastic medium is given as where λ and  are the lame's constants.i U is the displacement of particles in the elastic solid.In terms of the displacement components, the equations of motion are expressed as follows: Latin American Journal of Solids and Structures 12 (2015) 1870-1900 To seek the harmonic solution of system of equations ( 8), for the propagation of plane waves, the displacement components are written as follows: where  is angular frequency and   corresponds to the polarization for the motion of solid particles in the elastic solid medium.Substituting ( 9) in (8) yields a system of three equations, given by Now, after re-adjusting the above equation, we get which are the Christoffel's equations for the propagation of harmonic plane waves in the elastic solid medium.The coefficients used in the above the above relation are In terms of velocity V , the slowness is defined as The dual (complex) vector  N represents the directions of propagation and attenuation of a wave in the porous medium.In terms of  N and V , the Christoffel equations ( 11) are expressed as The non-trivial solution for Christoffel equations is ensured by vanishing the determinant ( This condition translates into two equations as follows: The first one (i.e., 1 0 The root of this square equation define the velocity ( 1 V ) in the elastic medium.In this case the polarization vector   1 2 3 , , , S S S    corresponding to equation ( 12), is calculated to be parallel to 

N
and hence the wave identified with velocity 1 V is longitudinal wave.Another equation (i.e., 2 0   ) yields Latin American Journal of Solids and Structures 12 (2015) 1870-1900 which implies a wave with velocity 2 V    .The corresponding polarization vector   is represented through a singular matrix   T    I N N .So, the polarization vector may be parallel to a column (or, row) vector of this symmetric matrix.This defines the direction of polarisation in a plane, which is normal to the propagation vector  N .This implies that the wave with velocity 2 V is a transverse wave.The polarisation vector  S defines the polarisation of solid particles in the elastic medium.

FORMULATION OF THE PROBLEM
Consider an elastic solid and a porous solid containing liquid filled bound pores and a connected pore space saturated by two-phase fluid having a common boundary.In the cartesian co-ordinate system   , , x x x , let the plane 3 0, x  define this common boundary which is separating the two different media (say, 1 M and 2 M ).A wave (P or SV) travels through the elastic medium 1 M (i.e., 3 0 x  ) with velocity 0 V and incident at the interface making an angle 0  to the 3 x -axis pointing in to this medium.For two dimensional motion in the 1 3 x x  plane, unit vector   0 0 sin , 0, cos   represents the phase direction of the incident wave.The incident angle may vary from 0 to 2  .Such an incidence results in two waves reflected back in to elastic medium 1 M and four waves refracted to the continuing porous medium 2 M .The primed quantities separate the medium 1 M from 2 M .
The displacement in the elastic medium 1 M are expressed as where values 1-2 of the index l represent the P and SV waves, respectively.The l f  are relative excitation factors for two reflected waves.We have     1 l l j j n n    and from Snell's law, ( ) Similarly, for waves refracted in medium 2 M , the displacements are expressed as where  L AS and G  BS are as given in Sharma and Saini (2012).
The l f are relative excitation factors for refracted waves.We have n n  , and from Snell's law,     Latin American Journal of Solids and Structures 12 (2015) 1870-1900

BOUNDARY CONDITIONS
We assume that two half spaces separated by a plane interface 3 0, x  are in perfect contact.Therefore, the boundary conditions are continuity of stress components and displacement components along the interface plus one more condition which restrict the flow of two fluids of porous solid in to uniform elastic solid, i.e, at 3 0.
where superposed dot represent the temporal derivative.
The above boundary conditions are satisfied through a system of six linear inhomogeneous equations in 1 f , 2 f , 3 f , 4 f , 5 f and 6 f .This system of equations is given by where the coefficients ij a , ( 1, 2, 3, 4, 5, 6 i  ) are given as follows: Latin American Journal of Solids and Structures 12 (2015) 1870-1900

ENERGY RATIOS
We now consider the partitioning of energy between different reflected and refracted waves at the surface element of unit area.Following Achenbach (1973), the rate at which the energy is transferred per unit area of the surface is given by the scalar product of surface traction and particle velocity, denoted by * P .The time average of * P over a period, denoted by * P , represents the average energy transmission per unit surface area per unit time.Thus, on the surface with normal along 3 xdirection, the average energy intensities of the waves in the uniform elastic solid medium are defined by with the help of expressions for two arbitrary complex functions f and g , we obtain the energy ratios giving the rate of average energy transmission corresponding to each of the reflected and refracted waves to that of the incident wave.These energy ratios ), for the reflected P and SV waves, respectively, are defined as follows: where Latin American Journal of Solids and Structures 12 (2015) 1870-1900 For a saturated porous solid half space with normal along the 3 x -direction, the average energy flux is represented through the components ij P given by To explain the distribution of incident energy at the free surface of a dissipative porous medium, a matrix defined with its element given by, where bar over an entity implies its complex conjugate.Elements of the matrix P in ( 27) are defined as follows: The matrices used in the above expressions are defined as where the superscript '( )' j on matrices A and B means the matrices are evaluated for slowness vector p of the corresponding refracted wave represented with a value of   1, 2, 3, 4 j  .An energy matrix (27) calculates the distribution of the energy among four waves travelling into the dissipative porous medium saturated by two miscible/immiscible fluids.The diagonal entries of the energy matrix ij E represent the energy share of the four refracted waves in the medium.The terms 11 E , 22 E , 33 E , 44 E are identified as the refraction (energy) coefficients for I P , II P , III P and SV waves, respectively.The sum of all non-diagonal entries of this energy matrix gives the share of the interaction energy among all the refracted waves in the medium.This part of energy is given by , yields the conservation of the incident energy across the interface through relation 1 Latin American Journal of Solids and Structures 12 (2015) 1870-1900

NUMERICAL EXAMPLE
A reservoir rock (sandstone) saturated with a liquid and gas is chosen for the numerical model of porous medium (Garg and Nayfeh, 1986).The solid grains of the rock with bulk modulus s K  36 GPa , rigidity modulus s   9 GPa, and density s   2650 kg/m 3 form a porous frame of porosity The connected pore space is filled with the bubbles of gas of bulk modulus g K  0.0037 GPa and density g   100 kg/m 3 mixed in a liquid of bulk modulus l K  2.3 GPa and density l   .The elastic medium ( 1M ) has the parameters,   2650 kg/m 3 ,   16.9 GPa, and   28.3 GPa as the values chosen for the density and elastic constants for granite.

Reflection and refraction coefficients
We find that the energy ratios of reflected and refracted waves due to either an incident P wave or SV wave depend upon the angle of incidence (i.e., 0  ).The energy ratios ) and the energy matrices ij E ( , 1, 2, 3, 4 i j  ) defined in the previous section are calculated for a given value of the incident angle 0  varying from 0  to 90  .In the present case, the incidence is considered only for Pand SVwaves.For each incidence, the energy partitions are computed and the conservation of energy is ensured at the interface of an elastic solid half space and the porous solid half space containing liquid filled bound pores and a connected pore space saturated by two-phase fluid.In the present discussion, 1 E and 2 E denote the reflection coefficients of P and SV waves, respectively, whereas 11 E , 22 E , 33 E , 44 E denote the refraction coefficients of the I P , II P , III P , and SV waves, respectively.The amplitude ratios for the reflected and refracted waves are donated as the energy ratios except these are denoted by notation ( represent the amplitude ratios of refracted I P , II P , III P , -SV , reflected P and SV-waves, respectively.Henceforth in the discussions the notations for the energy and amplitude ratios will be used for sake of convenience.Fig. 2 shows the variation of amplitude ratios with gas saturation for incident P -wave. 1 Z , 2 Z and 11 Z increases with increase in the gas saturation.44 Z is almost ineffective to the variation in the gas saturation.22 Z and 33 Z are very small yet they show variation with  as depicted in Fig. 2.
Fig. 3 shows the variation of amplitude ratios with change in the frequency for incident P -wave.

2
Z and 11 Z gets strengthen with increase of frequency.Unlike 2 Z , the amplitude ratio 44 Z decreases in magnitude with increase in the frequency.1 Z remains ineffective to the change in the frequency.As 22 Z and 33 Z are negligibly small yet their amplitude ratios get strengthen with increase in frequency.
Fig. 4 shows the variation of the amplitude ratios with change in the capillary pressure for incident P-wave.It is evident from the figure that capillary pressure has a little effect upon the variation of amplitude ratios.Whatever effect, if seen, is at the small scale, which is clearly visible in the amplitude Latin American Journal of Solids and Structures 12 (2015) 1870-1900  Figure ,   0.2, 0.5, 0.8, and different angles of incidence of P-wave.
ratios 22 Z and 33 Z .Where, 22 Z gets weaken with increase in the capillary pressure and 33 Z shows an anomalous behaviour with increase of capillary pressure.
, and different values of  .
1 Z , first decreases until 58 o , thereafter, it gets little strengthen with increase in the fraction of bound liquid film.Thus, 1 Z shows a mixed behavior with variation in the fraction of bound liquid film.Again, 22 Z and 33 Z are very small, but yet their amplitude ratios get weaken with increase in the fraction of bound liquid film.Structures 12 (2015) 1870-1900  Figure

Latin American Journal of Solids and
Z and 44 Z ) with angle of incidence of P-wave for , and different values of cap K .
Fig. 6 shows the variation of amplitude ratios with gas saturation for incident SV-wave. 1 Z has a mixed behaviour with increase in the gas saturation, i.e., it gets strengthen until the angle of incidence Latin American Journal of Solids and Structures 12 (2015) 1870-1900  Figure Z and 44 Z ) with angle of incidence of P-wave for   0.4, , and different values of  .
42 o and thereafter, it starts weakening with increase in the gas saturation. 2 Z has no variation till the angle of incidence 28 o , after that it increases in strength with increase in the gas saturation till 63 o , beyond this some weakening is observed with increase in the gas saturation.11 Z decreases in strength with increase in the gas saturation σ until the angle of incidence 42 o , after that increase is observed with  .Similarly, 44 Z increases little with increase in the gas saturation until 48 o , afterwards, it decreases a little with increase in the gas saturation.Structures 12 (2015) 1870-1900 Figure 6: The same as Fig. 2, but variation for incidence of SV-wave.Fig. 7 shows the variation of the amplitude ratios with change in the frequency  for incident SVwave. 1 Z strengthen with increase in the frequency until 42  , after this it remains ineffective to the change in frequency.Therefore, mixed behavior is shown with variation in the frequency.44 Z gain in strength with increase in  . 2 Z increases with increase in  , change is most prominent after 48 o .While 11 Z shows almost no variation with  .22 Z and 33 Z are very small, but it strengthen with increase in the frequency.Structures 12 (2015) 1870-1900 .Fig. 8 shows the variation of amplitude ratios with change of capillary pressure for incident SV-wave.

Latin American Journal of Solids and
Like Fig. 4, 1 Z , 2 Z , 11 Z and 44 Z show almost no variation with increase in the capillary pressure as the variation is observed at the small scale only, which is visible in 22 Z and 33 Z .
Fig. 9 shows the variation in the amplitude ratios with fraction of bound liquid film for incidence of SV-wave.44 Z decreases with increase in  . 1 Z shows a mixed behaviour, i.e., it gets weaken with increase in the fraction of bound liquid film until the angle of incidence 42 o , thereafter, it increases with  until 55 o , henceforth, no variation till grazing incidence.Therefore, mixed behavior is observed.Z increases with increase in  till 42 o , after that, it decreases with increase in the fraction of bound liquid film.22 Z and 33 Z are again very small and shows a weakening with increase in the fraction of bound liquid film.
Fig. 12 shows the effect of the capillary pressure on the energy partitions for the incident P-wave.
The variation is observed at the small scale. 1 E and 44 E do not show any change with capillary pressure.2 E shows an increase in the energy with increase in the capillary pressure.Unlike 2 E , 11 E shows a decrease in energy with increase in the capillary pressure.Fig. 13 shows the effect of bound liquid film on the energy partitions for the incident P-wave. 1 E shows a little variation that too after 65  , i.e., dips with increase in the fraction of the bound liquid film.E increases in magnitude with increase in the frequency till the critical incidence 39 o . 2 E is insensitive to the change of frequency.11 E decreases too with increase in the frequency.But, 44 E shows an increase with increase in the frequency.However, 22 E , 33 E are small but show an increase in the magnitude with increase in the frequency.RR E shows a decrease in the magnitude with increase in the frequency till 39  and beyond which the pattern is reversed, i.e., increases with increase in the frequency.decreases with increase in the fraction of the bound liquid film.22 E and 33 E show a decrease with increase in the value of  , however, the variation is at small scale.RR E increases till 39 o and beyond which it decreases with increase in the value of the  .

CONCLUDING REMARKS
The work presented here study the reflection/refraction at the interface of an elastic solid and a partially saturated porous solid containing liquid filled bound pores and a connected pore space saturated by two-phase fluid and the study is for low frequency regime.The porous medium is dissipative due to presence of viscous fluids in the connected pores.The four attenuated waves in the porous medium are identified with complex velocities.The variable gas share in pores enables to represent the pore saturation from all liquid to all gas.Besides, the method used in this paper is not Latin American Journal of Solids and Structures 12 (2015) 1870-1900 based on elastic Lame's potentials, but looks directly for the solution of the elastodynamic equation in term of displacement vectors.The study through this method has advantage over the traditional potential method.Firstly, it could be generalized to anisotropic media for the given model.Secondly, it could be used for the inhomogeneous media, but lame's potential method could not be used for both the cases.Some main observations from the numerical example may be important and hence are explained as follows: 1) The amplitude ratios 1 Z and 2 Z show variation with change of the saturation of gas, both for the incident P and SV-waves.2) For incidence of P-wave, 2 Z , 11 Z and 44 Z show change with change of frequency.While for incidence of SV-wave, 1 Z , 2 Z , and 44 Z show a significant change with change of frequency.3) All the amplitude ratios are insensitive to the change of the capillary pressure.4) All amplitude ratios are sensitive to change of fraction of bound liquid film.5) The variation with saturation of gas is observed for 22 E and 44 E for incidence of P-wave and for incident SV-wave, the variation is prominent for 1 E , 11 E and 44 E waves only.
6) The variation of energy with frequency for incidence P-wave is most prominent for reflected SV and refracted SV-wave.For SV incidence, the significant variation is observed for reflected P, refracted I P and refracted SV .
7) The capillary pressure has an insignificant effect on the energy ratios variation.Any effect, if there, is very small and could be observed for small energy ratios 22 E and 33 E .8) Variation with bound liquid film is dominant only in case of 2 E , 11 E and 44 E waves for incident P-wave.For incidence of SV-wave, 11 E and 44 E alone show variation with change of fraction of bound liquid film.9) The energy ratios 22 E and 33 E are very small.10) Angle of incidence render a significant effect on the energy partitions across the surface.11) The sum of all the energies i.e. reflected as well as refracted at the interface is unity.This shows that there is no dissipation of energy at the interface.Hence, energy is conserved.

Appendix:
The elastic constants used in (4) are expressed in terms of the elastic constants of the constituents, as follows:   , where h  is Henry's constant (Garg and Nayfeh, 1986) to represent the mixing of two pore-fluids.For immiscible pore-fluids, i.e., 0 h   , ik  are reduced to ik  , which are expressed as follows: Latin American Journal of Solids and Structures 12 (2015) 1870-1900 , where cap K is equivalent bulk modulus for macroscopic capillary pressure (Garg and Nayfeh, 1986) and d K is the bulk modulus for drained porous solid.j K denotes the bulk modulus of phase ( , , , ) j s l g   .
1 c K  represents the effective compressibility of mixture of pore-fluids and c K reduces to l K when 0   .The present work is specifically restricted to the propagation low frequency harmonic waves so as to follow Poiseuille flow.As a result the capillary pressure in pores is assumed to be independent of frequency and hence a constant cap K .

Figure 1 :
Figure 1: Shows the schematic diagram of the incidence, reflection and refraction of waves.
m .Both the pore-fluids are viscous and the values chosen for dissipation coefficients are l d  1 MPa-s/m 2 and g d  0.04 MPa-s/m 2 .The same liquid with viscosity    10 -12 GPa-s and Reynolds number Re   100 is assumed in the bound pores.The value of ratio cap l

Fig. 5 Figure 3 :
Fig. 5 shows the variation of the amplitude ratios with the change in the fraction of bound liquid film for incident P-wave. 2 Z and 11 Z show a weakening in the amplitude ratio with increase in the fraction of bound liquid film.44 Z increases with increase in the fraction of bound liquid film.While Latin American Journal of Solids and Structures 12 (2015) 1870-1900

Figure 7 :
Figure 7:The same as Fig.3, but variation for incidence of SV-wave.

Figure 8 :Z
Figure 8: The same as Fig. 4, but variation for incidence of SV-wave.

Figure 9 :
Figure 9:The same as Fig.5, but variation for incidence of SV-wave.

Fig. 11
Fig. 11 depicts the effect of the frequency on the energy partitions for incident P-wave.One could easily draw inference from the figure that 1 E and 11 E are insensitive to the change in the frequencyand if any little is seen that is seen after 78  .The share of reflected SV-wave (i.e., 2

Fig. 15
Fig. 15 depicts the effect of the frequency on the energy partitions for incident SV-wave. 1

Figure 14 :Figure 15 :
Figure 14: The same as Fig. 10, but variation for incidence of SV-wave.

Figure 16 :Figure 17 :
Figure 16:The same as Fig.12, but variation for incidence of SV-wave.

Latin American Journal of Solids and Structures 12 (2015) 1870-1900
2 E decreases with increase in the value of  .44 E increases gradually with increase in the E increases with increase in the saturation of gas. 2 E is insensitive to the change of gas saturation.11 E decreases with increase in the gas saturation and the decrease is quite significant.44 E increases with increase in the gas saturation.22 E , 33 E and RR E are very small and their variation is as shown in the figure.

Latin American Journal of Solids and Structures 12 (2015) 1870-1900
. 16 shows the effect of capillary pressure on the energy partitions for incidence of SV-wave. 1 Fig