# Propagation of Pore Pressure and Stress in Saturated Porous Media Based on a Darcy-Brinkman Formulation

Authors: Duoxing Yang; Lianzhong Zhang
Journal: Geofluids
 (2021)

---

## Abstract

Propagation of pore pressure and stress in water-saturated elastic porous media is theoretically investigated when considering the Darcy-Brinkman law. The wave mode, phase velocity, phase lag, damping factor, and characteristic frequency are found from the updated mathematic model. The Brinkman term describes the fluid viscous shear effects and importantly contributes to the dispersion relation and wave damping. The coincidence of the properties of Biot waves of the first and second kinds occurs at a characteristic frequency, which is remarkably influenced by the Brinkman term. A key finding is that, compared to the Darcy-Brinkman law, Darcy’s law overestimates the phase velocity, damping, and phase lag of the first wave, while underestimates the phase velocity, damping, and phase difference of the second wave. The introduction of the Darcy-Brinkman law yields an improved description of the damping of the compressional wave modes in saturated porous media.


---

## Body

## 1. Introduction

Propagation of acoustic waves through saturated porous media is a phenomenon of great fundamental importance and increasing technological relevance [1], which is significant to obtain better surface acoustic wave device’s performance. Compressional waves in saturated porous media have found geophysical interests to hydraulic fracturing, oil exploration, volcanic eruptions, and magneto-electro-elastic behavior of piezoelectric nanostructures [2]. Different investigations are available in literature focusing on the mechanical behaviors of piezoelectric nanostructures [3, 4] on fields of linear and nonlinear, longitudinal and transverse, free and forced wave propagation in microscale. Initial and couple stress influence on transmission of the Love-type wave [5], and Rayleigh waves [6] in material layers with imperfect interface [7], and surface wave propagation between viscous liquid and initially stressed piezoelectric half-space [8] are investigated by using the spatial varying WKB (Wentzel-Kramers-Brillouin) method [9, 10]. Solutions are obtained by a perturbation technique in closed-form for normal, shear stresses, dielectric displacement, and electric potential [11], and it is found that the complex phase velocity with positive imaginary part increases with time [12].As is well known, the first and second kind of compressional waves exist in saturated porous media, which are predicted by the Biot theory [13, 14]. The wave mode, dispersion relation [15, 16], and transfer function [17] are proposed. Generally, the first kind of compressional wave, which is known as Biot’s wave of the first kind or the first wave, is a weakly decaying high-velocity mode, while the second kind of compressional wave, which is known as Biot’s wave of the second kind or the second wave, is a strongly damping slow-velocity mode [18–21]. The nonlinear waves are characterized with diffusive, oscillatory, and shock types in the two-phase fluid saturated porous media [22–24]. Strong viscous coupling causes the coalescence of the two wave fronts into a single front, and the saturated porous medium behaves a single continuum exhibiting internal dissipation [16]. The coincidence of the properties of Biot waves of the first and second kinds occurs at a characteristic frequency, at which the properties of the first and second waves coincide [25]. Under conditions of two-phase flows, the resonance frequency of the gas is the governing factor for the transition of the wave modes [26]. The strong decay of compressional waves is attributed to the viscous coupling effect [27, 28] and the Darcy viscous resistance [29–31]. Elastic waves in non-Newtonian (Maxwell) fluid-saturated porous media are deeply investigated [32, 33], and it is found that the introduction of a Maxwell fluid into the Biot theory results in an overall increase of the attenuation of three Biot waves in the intermediate frequency domain except in the deeply non-Newtonian regime.Based on the Biot theory, the interaction forces due to the relative motion between fluid and solid phases [14, 17, 34] are formulated by the linear Darcy’s law, which relates macroscopic pressure gradients to the flow velocity vector in the fluid phase, describes the solid viscous drag effects [35], and hence, neglects the internal friction arising from viscous shear stress [36, 37]. For the flow of fluids in dual-porosity media (e.g., macropores and micropores) or miniature-sized porous structures, viscous shear stress is significant [1, 37, 38]. Rigorous usage of spatial averaging methods [39] results in the Darcy-Brinkman equation, which is known as Darcy’s law with Brinkman correction or the Darcy-Brinkman law [35]. The Brinkman term in the Darcy-Brinkman law [36, 37, 39] describes the fluid viscous shear effects that oppose the flow through the porous structures [36]. Besides the drag considered in Darcy’s law, the Darcy-Brinkman law accounts for internal friction [37, 38] which is described by the Stokes equation at the fluid-porous medium boundary, and has considerable influences on the momentum transfer between the solid and the fluid in porous media [40].Biot [14] pointed that the assumption of Darcy’s law breaks down if the frequency exceeds a certain value, known as Biot’s critical frequency above which the compressional wave attenuation is controlled by the viscous shear stress and the inertial effects. However, Darcy’s law does not include the internal friction and the inertia effects. Recently, a macroscopic model for describing compressional wave propagation in saturated porous media was proposed, and the Forchheimer term (e.g., inertial effects) was added to Euler’s equation without the Brinkman term included [41]. A Darcy-Brinkman-Biot formulation [37] was developed, which approximates the Navier-Stokes equations in fluid-filled macropores and resembles the equations for poroelasticity in microporous regions. The Darcy-Brinkman-Biot formulation neglects inertial effects (e.g., Forchheimer term) in a momentum transfer term (drag) between the fluid and the solid phases, given that the Reynolds number within the microporous region is small. There arises a problem in which it is required to identify the effects of internal friction on the compressional wave propagation in saturated porous media. That is why the emphasis here is placed on investigation of the wave propagation and damping to the Darcy-Brinkman law.In this paper, with the help of the Darcy-Brinkman law [35, 37, 39], the momentum conservation laws [14, 17] in their linearized harmonic form, which describes the one-dimensional propagation of compressional waves in saturated porous media, are extended from the Biot theory. The dispersion relation of the wave modes is reformulated. A characteristic frequency for the coincidence of the properties of Biot waves of the first and second kinds is found. We analyze the conditions under which the properties of the two compressional waves coincide. The transfer functions for pore pressures and stresses are constructed. Comparisons between the Darcy-Brinkman law and Darcy’s law in the compressional wave motion are carried out. The introduction of the Darcy-Brinkman law may yield a better understanding of the compressional waves in saturated porous media.


## 2. Mathematical Formulation of the Problem

The governing equations of the Biot theory [14, 17] and the Darcy-Brinkman-Biot formulation [37] in saturated porous media were established through use of the conservation laws. Sufficient detail of the mass and momentum conservation laws in time domain was provided by Biot [14] and by Carrillo and Bourg [37]. Recently, based on the Biot theory, the equations for compressional waves in frequency domain was developed by van der Grinten et al. [17]. Here, primary attention is paid to the introduction of the Brinkman term to the Biot theory in frequency domain. Suppose that caret (^) represents the complex amplitude, where j2=−1. p^ and σ^ are the pore pressure and effective stress.u^ and w^ represent the actual velocities of the fluid and solid phases. φ and κ are the porosity and frequency-independent permeability. ρl and ρp denote the intrinsic densities of the fluid and solid phases, respectively. ω is the angular frequency. Kl and Kp are, respectively, the bulk modulus of the fluid and the constrained modulus of the solid. The added mass coefficient α∞ is related to ρa by ρa=α∞−1φρl. The critical Biot frequency [14] is given by ωc=ηφ/κρl.It is assumed that the intrinsic density of the solid is constant, and the solid elastic deformation is fully described in terms of the porosity [14, 17]. Considering a ejωt dependence for all relevant variables, the equations governing the one-dimensional propagation of compressional waves in the linearized harmonic form [17] are the following.
(1)jωφρlφρl^+ρlφ^+∂w^∂x=0,(2)−jωφ^+1−φ∂u^∂x=0,(3)jωφρlw^=−φ∂p^∂x+jωρa+φ2ηκu^−w^,(4)jω1−φρpu^=−∂σ^∂x−1−φ∂p^∂x−jωρa+φ2ηκu^−w^.Equations (1) and (2) are the mass conservation laws for the fluid and the solid, respectively. Equations (3) and (4) represent the conservation of momentum of the fluid and the solid, respectively, as stated in the Biot theory [14]. The term, φ2η/κu^−w^ in Equations (3) and (4) is the fluid-solid momentum transfer term [14, 37] and describes the frequency-dependent interaction forces between fluid and solid and is expressed as a function of the relative velocities of the two phases [37]. The high frequency added mass, jωρa describes the purely imaginary interaction force jωρau^−w^in the limit for high frequencies [17]. The Brinkman term [39] is defined byφη∂2u^−w^/∂x2, which recovers the viscous drag effects and uses the fluid viscosity coefficient [36–38]. According to the momentum conservation laws of the Biot theory [14], we introduce the Brinkman term to the momentum conservation of the fluid and the solid, as stated in the reference [37]. From Equations (3) and (4), it follows that
(5)jωφρlw^=−φ∂p^∂x+φη∂2u^−w^∂x2+jωρa+φ2ηκu^−w^,(6)jω1−φρpu^=−∂σ^∂x−1−φ∂p^∂x−φη∂2u^−w^∂x2−jωρa+φ2ηκu^−w^.Equations (1), (2), (5), and (6) is the Darcy-Brinkman-Biot formulation in the frequency domain under the rigorous assumptions present in the Biot theory [14, 17]. Inserting a harmonic x-dependence e−jωx/c for all dependent variables in equations (1), (2), (5), and (6), and following the procedure [17], the constitutive equations of pore pressure and stress are given by
(7)p^=Klw^c+1−φφu^c,(8)σ^=Kpu^c.Substituting Equations (7) and (8) into Equations (5) and (6), it yields
(9)w^ρ22−φKlc2−jφ2ηkω+jφηωc2=u^−ρ12+1−φKlc2−jφ2ηkω+jφηωc2,(10)w^−ρ12+1−φKlc2−jφ2ηkω−jφηωc2=u^ρ11−Kpc2−1−φ2Klφc2−jφ2ηkω−jφηωc2.By virtue of Equations (9) and (10), the dispersion relation takes the following form
(11)χc4+βc2+γ=0,with
(12)χ=ρ11ρ22−ρ122−jφ2ηκωρ11+ρ22+2ρ12,β=−ρ22Kp−1−φ2Klρ22φ−φKlρ11+21−φKlρ12+jφηKl+φKpκω+jφηωρ11−ρ22,γ=φKlKp+jφηω−Kp+2φ−1φKl.The termsjφηωρ11−ρ22 and jφηω−Kp+2φ−1/φKl in the coefficients of β and γ implies the contribution due to the Brinkman term. c is the complex phase velocity of compressional waves.Use has been made of the notation(13)ρ11=1−φρp−ρ12,ρ22=φρl−ρ12,ρ12=−ρa.By virtue of(14)λ1=ρ11ρ22−ρ122,λ2=−φ2η/κρ11+ρ22+2ρ12,λ3=−ρ22Kp−1−φ2Klρ22/φ−φKlρ11+21−φKlρ12,λ4=φηKl+φKp/κ,λ5=φηρ11−ρ22,λ6=φKlKp,λ7=φη−Kp+2φ−1/φKl.Equation (15) describes the two independent solutions of Equation (11), representing the first and second waves, which are denoted by the subscripts 1 and 2.
(15)c1,22=−λ3+jλ4ω−1+jλ5ω±λ3+jλ4ω−1+jλ5ω2−4λ1+jλ2ω−1λ6+jλ7ω2λ1+jλ2ω−1.LetΦω=λ3+jλ4ω−1+jλ5ω2−4λ1+jλ2ω−1λ6+jλ7ω. We define
(16)ReΦω=λ32−λ4ω−1+λ5ω2−4λ1λ6+4λ2λ7,ImΦω=2λ3λ4−2λ2λ6ω−1+2λ3λ5−2λ1λ7ω.At a nonzero frequency, if the conditionReΦω=0is satisfied, this implies that
(17)ω=ω∗=λ32−2λ4λ5−4λ1λ6+4λ2λ7+λ32−2λ4λ5−4λ1λ6+4λ2λ72−4λ42λ522λ52.BecauseImΦω=0, it has λ1λ4λ7=λ2λ5λ6. Therefore, if λ1λ4λ7=λ2λ5λ6, then at a characteristic frequency ω=ω∗, the properties of the first and second waves are identical [25], and the biquadratic Equation (11) has a single solution, namely, c1=c2=2/2−λ3−jλ4ω−1−jλ5ω/λ1+jλ2ω−1, which corresponds to the classical wave propagating with the mixture wave speed [15]. That is, in the saturated porous medium, a single compressional wave can propagate [25]. Note that the characteristic frequency of compressional waves represents an intrinsic property of saturated porous media which rigorously obey the criterion λ1λ4λ7=λ2λ5λ6 and only depends on the properties of the fluid and solid phases.Supplemented with the prescribed boundary conditions ofp^ and σ^ at x=0, p^0 and 0, respectively, through use of Equations (7) and (8) and by superposition, the transfer functions [17] for the pore pressure and the stress waves travelling over a distance x are reformulated in the form
(18)Ηpx,ω=p^x,ωp^0ω=1β1−β2β1+1−φφe−jωx/c1−β2+1−φφe−jωx/c2,(19)Ησx,ω=σ^x,ωp^0ω=KpKl1β1−β2e−jωx/c1−e−jωx/c2,with
(20)βi=w^iu^i=−ρ12+1−φKl/ci2−jφ2η/kω+jφηω/ci2ρ22−φKl/ci2−jφ2η/kω+jφηω/ci2.The termjφηω/ci2 denotes the contribution due to the Brinkman term. Equations (11), (18), and (19) constitute the fundamental equations of compressional wave fields of frequency domain in saturated porous media with the Darcy-Brinkman law included.


## 3. Results and Discussion

### 3.1. Model Evaluation and Validation

Experimental benchmark results [15, 17] from a shock tube technique are used to validate the proposed fundamental equations of compressional wave fields of frequency domain with the Darcy-Brinkman law included. In the experimental setup, a vertical shock tube is utilized to trigger a step pressure loading at the top of a porous cylinder located in the test section of the tube. By measuring simultaneously pore pressures and strains at observation points, the in-phase and out-of-phase wave modes and the wave speeds are observed. The following mechanical properties [17] of a sandstone saturated with water are utilized: φ=0.3, ρl=990kg/m3, ρp=2650kg/m3, η=1.0×10−3kgm−1s−1, Kl=2.2GPa, Kp=4.5GPa, κ=5.6×10−11m2, and α∞=1.45. Since the constrained modulus of the solid and the added mass are certainly the important effects, it has been verified from the experiments [15, 17] that the effective constrained modulus of the solid has a value of 4.5 GPa, and the effective added mass parameter is 1.45. Quantitative information on compressional wave speeds is compared with theoretical predictions, as shown in Table 1.Table 1
Comparison of calculated and observed wave velocities.
ModeObserved [15]Observed [17]Darcy’s lawDarcy-Brinkman formulationFirst wave (km/s)2.27±0.062.19±0.092.782.36Second wave (km/s)0.78±0.020.79±0.020.670.78The calculated wave speeds from the Darcy-Brinkman law are well agreement with the observed wave speeds. The introduction of the Brinkman term clearly improves the description of the damping for the second compressional wave [17]. The primary emphasis here is that the fluid density and porosity are rigorously dependent of pore pressure and effective stress by the constitutive Equations (18), that the intrinsic density is a constant, and that the compression or expansion of the solid can be fully described in terms of the porosity.Bouzidi and Schmitt [42] measured the speed of the compressional waves in a water-loaded, porous sintered glass bead plate with a novel pair of ultrasonic transducers consisting of a large transmitter and a near-point receiver. Updated observations of the two compressional wave modes are obtained. The following mechanical properties [42] are utilized: φ=0.391, ρl=995kg/m3, ρp=2445kg/m3, η=1.0×10−3kgm−1s−1, Kl=2.209GPa, Kp=4.826GPa, κ=1.99×10−11m2 and α∞=1.44. The velocities are calculated from Equation (11) given in Table 2.Table 2
Comparison of calculated and observed wave velocities at 780 kHz.
ModeObserved [42]Gassmann’s relationsDarcy’s lawDarcy-Brinkman formulationFirst wave (km/s)2700±55266725682362Second wave (km/s)1010±20—837911Compared to the experimental results, Darcy’s law and the Darcy-Brinkman formulation underestimate the velocities of the compressional waves. One key factor included in the dispersion relation (Equation (11)) is the effective constrained modulus of the solid Kp that is not available in the reference [42]. In the calculation, the frame bulk modulus is used. As stated by Bouzidi and Schmitt [42], other factors such as heterogeneity of the porosity might partially lead to this discrepancy. The key fact is that the fluid density and porosity are rigorously related to pore pressure and effective stress by the constitutive equations and that the deformation of the solid can be fully described in terms of the porosity. In the experiments, the emitted ultrasonic waves cannot cause the change of porosity.Albert [16] simulated wave propagation in a water-saturated porous material. Here, we demonstrate that in the limit where the Brinkman term vanishes, the result for Darcy’s law is recovered. Formation and saturated parameters used are those as reported by Albert [16]. Figure 1 shows that the velocity and attenuation dispersion curves obtained using our model show significant dispersion of the same order as that predicted by the original Biot theory [16].Figure 1
Velocity (a) and attenuation (b) predicted with Darcy’s law for the compressional waves in water-filled sand. The velocity is normalized by the fluid velocity 1517 m/s as a function of the normalized frequency,Ω=ω/ωc.
(a)(b)


### 3.2. Phase Velocity and Damping of the First and Second Kind of Waves

The following mechanical properties [17] of a sandstone saturated with water are used: φ=0.3, ρl=990kg/m3, ρp=2650kg/m3, η=1.0×10−3kgm−1s−1, Kl=2.2GPa, Kp=4.5GPa, κ=5.6×10−11m2, and α∞=1.45. Figures 1–3 illustrate the normalized frequency-dependent (Ω=ω/ωc) dispersion profiles of the phase velocity, damping, and phase lag. As expected, the first wave (see Figure 2) is a slightly decaying high-velocity mode, and the relative motion of fluid and solid phases is the in-phase [26, 32, 33], exhibiting a pressure-relaxation phenomenon [30], whereas the second wave (see Figure 3) is a strongly damping slow-velocity mode and the movement of fluid and solid phases is the out-of-phase [32, 33], showing the feature of diffusion waves [28], in which the pore pressure predominates over the effective stress [26]. As seen in Figure 4, under the condition of λ1λ4λ7=λ2λ5λ6, the properties of the first and second waves coincide. The two wave fronts emerge into a single front [16], and the saturated porous medium behaves a single continuum exhibiting internal dissipation. Under a certain relationship λ1λ4λ7=λ2λ5λ6 between the fluid and solid parameters in which only exists a single wave, the phase velocity is underestimated, while the damping factor and phase lag are overestimated by Darcy’s law. Compared to the Darcy-Brinkman law, Darcy’s law overestimates the phase velocity, damping factor, and phase lag of the first wave, while underestimates the phase velocity, damping factor, and phase difference of the second wave. The Darcy-Brinkman law does not lead to the change of the in-phase motion into out-of-phase motion in the Frenkel-Biot wave modes [26], as can be seen that the second wave is of opposite polarization compared to the first compression wave [32, 33], which is consistent with the Gassmann equation [43]. The effects of the internal friction arising from viscous shear stress [8, 36, 37] on propagation of compressional waves are of significant importance.Figure 2
Phase velocity (a), damping factor (b), and phase difference (c) of the first wave versus normalized frequency (Ω=ω/ωc). The damping is the imaginary part of −ω/ci,i=1,2. The phase velocity is the real part of the complex wave velocity. The red line denotes the Darcy-Brinkman law and the black line Darcy’s law.
(a)(b)(c)Figure 3
Phase velocity (a), damping factor (b), and phase difference (c) of the second wave versus normalized frequency (Ω=ω/ωc). The red line presents the Darcy-Brinkman law and the black line Darcy’s law.
(a)(b)(c)Figure 4
Phase velocity (a), damping factor (b), and phase difference (c) of the wave versus normalized frequency (Ω=ω/ωc) under the condition of λ1λ4λ7=λ2λ5λ6. The red line presents the Darcy-Brinkman law and the black line Darcy’s law.
(a)(b)(c)


### 3.3. Wave Fields of Pore Pressure and Stress

Propagation of pore pressure and stress waves is analytically investigated in a one-dimensional semi-infinite domain subject to the step loading,∆pof pore pressure at t=0, x=0. Figure 4 shows a simulation of the frequency dependence of the amplitude of the pore pressure and the stresses over a wave travelling distance. From the comparison of Darcy’s law with the Darcy-Brinkman law, it is found that the Brinkman term is negligible in the near field, while the viscous shear stress significantly influences the wave motion in the far field (see Figure 5(a)). Figure 5(b) shows the interference pattern of the wave fronts. A key result is that during the wave motion, the amplitude of the pore pressure generally decreases, compared to the shear-induced pore pressure changes [44]. The magnitude of the stress firstly increases and then decreases back to its initial value. The amplitude-curve crossovers occur as the interplay between greater phase lags due to increasing frequency [45]. The two-wave modes [43, 46], dissipative and diffusive patterns [47], are demonstrated clearly.Figure 5
Normalized amplitude of pore pressure and stress wave versus travelling distance, (a) comparison of Darcy’s law with the Darcy-Brinkman law for an angular frequency of 1000 Hz, and (b) based on the Darcy-Brinkman effect with angular frequency as a parameter forη=1.0×10−3kgm−1s−1. The solid line denotes the transfer function of the pore pressure, and the dash-dot line represents the transfer function of the stress. The wave travelling distance is normalized by the viscous penetration length δ=2μ/ωρl at frequency of 2×10−6Hz.
(a)(b)Based on the Darcy-Brinkman law, Figure6(a) shows the simulation of the pore pressure and stress waves for various viscosities, η. The viscous dissipation of energy [35, 37] considerably influences the propagation of the pore pressure and stress waves. The pore pressure and stress waves exhibit dissipative and dispersive behaviors [24, 46–48]. The amplitude-curve crossovers are due to the spatial superposition of the confined waves acting like interference patterns [45, 49]. As can be seen in Figure 6(b), the amplitudes of pore pressure and stress are insensitive to the value of the added-mass coefficient. The regime where the pore pressure (or stress) curves coincide implies that the pore pressure (or stress) is in that case completely independent of α∞ and corresponds to the situationu^=w^ as shown from the field Equations (5) and (6). The pressure-coherent stress follows a single gamma distribution, compared to the flow-induced stress in saturated porous media [50].Figure 6
Normalized amplitude of pore pressure and stress versus wave travelling distance based on Darcy-Brinkman law withη as a parameter (a) and with added-mass coefficient α∞ as a parameter (b) at an angular frequency of 1000 Hz. The solid line denotes the transfer function of the pore pressure, and the dash-dot line represents the transfer function of the stress. The wave travelling distance is normalized by the viscous penetration length δ=2μ/ωρl at frequency of 2×10−6Hz.
(a)(b)


### 3.4. Discussions of the Darcy-Brinkman-Biot Formulation

It is noted that the assumption of Darcy’s law is not valid if the frequency exceeds a Biot’s critical frequency [14] above which the compressional wave attenuation is strongly controlled by the viscous shear stress, and the Darcy-Brinkman law should be considered to address the effects of non-Darcy flow on compressional wave motions. Moreover, for compressional elastic wave propagation in Newtonian fluid-saturated porous media, as pointed by Carrillo and Bourg [37], the Darcy term is not adapted to describe systems with both macropores and micropores, which contain two characteristic length scales. Fluid flow in two-scale porous media is governed by a general coupled system of equations, which approximates the Navier-Stokes equations in fluid-filled macropores and resembles the equations for poroelasticity in microporous regions.Now, it has to be realized that in the theoretical analysis, we have assumed that the nonlinear inertial terms in the momentum equations of the fluid and solid can be neglected for the frequency below the critical Biot frequency. This assumption is also reasonable to a small Reynolds number (Re<1) within the microporous region [37]. In the limit of high frequencies, the frequency-dependent interaction forces between solid and fluid can be described in terms of the high-frequency added mass jωρa. However, the Darcy-Brinkman-Biot formulation [37] does not include the Forchheimer equation.Future work involves the addition of Forchheimer term [51, 52], solid dissipation, and viscoelasticity and relaxation effects to the momentum conservation laws of compressional waves in saturated porous media. Of particular interest is that previous mathematical approaches [52, 53] of the Laplace transform, Fourier transform, and matrix transfer would be extended by solving the one-dimensional Darcy-Brinkman-Biot formulation of frequency domain.


## 3.1. Model Evaluation and Validation

Experimental benchmark results [15, 17] from a shock tube technique are used to validate the proposed fundamental equations of compressional wave fields of frequency domain with the Darcy-Brinkman law included. In the experimental setup, a vertical shock tube is utilized to trigger a step pressure loading at the top of a porous cylinder located in the test section of the tube. By measuring simultaneously pore pressures and strains at observation points, the in-phase and out-of-phase wave modes and the wave speeds are observed. The following mechanical properties [17] of a sandstone saturated with water are utilized: φ=0.3, ρl=990kg/m3, ρp=2650kg/m3, η=1.0×10−3kgm−1s−1, Kl=2.2GPa, Kp=4.5GPa, κ=5.6×10−11m2, and α∞=1.45. Since the constrained modulus of the solid and the added mass are certainly the important effects, it has been verified from the experiments [15, 17] that the effective constrained modulus of the solid has a value of 4.5 GPa, and the effective added mass parameter is 1.45. Quantitative information on compressional wave speeds is compared with theoretical predictions, as shown in Table 1.Table 1
Comparison of calculated and observed wave velocities.
ModeObserved [15]Observed [17]Darcy’s lawDarcy-Brinkman formulationFirst wave (km/s)2.27±0.062.19±0.092.782.36Second wave (km/s)0.78±0.020.79±0.020.670.78The calculated wave speeds from the Darcy-Brinkman law are well agreement with the observed wave speeds. The introduction of the Brinkman term clearly improves the description of the damping for the second compressional wave [17]. The primary emphasis here is that the fluid density and porosity are rigorously dependent of pore pressure and effective stress by the constitutive Equations (18), that the intrinsic density is a constant, and that the compression or expansion of the solid can be fully described in terms of the porosity.Bouzidi and Schmitt [42] measured the speed of the compressional waves in a water-loaded, porous sintered glass bead plate with a novel pair of ultrasonic transducers consisting of a large transmitter and a near-point receiver. Updated observations of the two compressional wave modes are obtained. The following mechanical properties [42] are utilized: φ=0.391, ρl=995kg/m3, ρp=2445kg/m3, η=1.0×10−3kgm−1s−1, Kl=2.209GPa, Kp=4.826GPa, κ=1.99×10−11m2 and α∞=1.44. The velocities are calculated from Equation (11) given in Table 2.Table 2
Comparison of calculated and observed wave velocities at 780 kHz.
ModeObserved [42]Gassmann’s relationsDarcy’s lawDarcy-Brinkman formulationFirst wave (km/s)2700±55266725682362Second wave (km/s)1010±20—837911Compared to the experimental results, Darcy’s law and the Darcy-Brinkman formulation underestimate the velocities of the compressional waves. One key factor included in the dispersion relation (Equation (11)) is the effective constrained modulus of the solid Kp that is not available in the reference [42]. In the calculation, the frame bulk modulus is used. As stated by Bouzidi and Schmitt [42], other factors such as heterogeneity of the porosity might partially lead to this discrepancy. The key fact is that the fluid density and porosity are rigorously related to pore pressure and effective stress by the constitutive equations and that the deformation of the solid can be fully described in terms of the porosity. In the experiments, the emitted ultrasonic waves cannot cause the change of porosity.Albert [16] simulated wave propagation in a water-saturated porous material. Here, we demonstrate that in the limit where the Brinkman term vanishes, the result for Darcy’s law is recovered. Formation and saturated parameters used are those as reported by Albert [16]. Figure 1 shows that the velocity and attenuation dispersion curves obtained using our model show significant dispersion of the same order as that predicted by the original Biot theory [16].Figure 1
Velocity (a) and attenuation (b) predicted with Darcy’s law for the compressional waves in water-filled sand. The velocity is normalized by the fluid velocity 1517 m/s as a function of the normalized frequency,Ω=ω/ωc.
(a)(b)


## 3.2. Phase Velocity and Damping of the First and Second Kind of Waves

The following mechanical properties [17] of a sandstone saturated with water are used: φ=0.3, ρl=990kg/m3, ρp=2650kg/m3, η=1.0×10−3kgm−1s−1, Kl=2.2GPa, Kp=4.5GPa, κ=5.6×10−11m2, and α∞=1.45. Figures 1–3 illustrate the normalized frequency-dependent (Ω=ω/ωc) dispersion profiles of the phase velocity, damping, and phase lag. As expected, the first wave (see Figure 2) is a slightly decaying high-velocity mode, and the relative motion of fluid and solid phases is the in-phase [26, 32, 33], exhibiting a pressure-relaxation phenomenon [30], whereas the second wave (see Figure 3) is a strongly damping slow-velocity mode and the movement of fluid and solid phases is the out-of-phase [32, 33], showing the feature of diffusion waves [28], in which the pore pressure predominates over the effective stress [26]. As seen in Figure 4, under the condition of λ1λ4λ7=λ2λ5λ6, the properties of the first and second waves coincide. The two wave fronts emerge into a single front [16], and the saturated porous medium behaves a single continuum exhibiting internal dissipation. Under a certain relationship λ1λ4λ7=λ2λ5λ6 between the fluid and solid parameters in which only exists a single wave, the phase velocity is underestimated, while the damping factor and phase lag are overestimated by Darcy’s law. Compared to the Darcy-Brinkman law, Darcy’s law overestimates the phase velocity, damping factor, and phase lag of the first wave, while underestimates the phase velocity, damping factor, and phase difference of the second wave. The Darcy-Brinkman law does not lead to the change of the in-phase motion into out-of-phase motion in the Frenkel-Biot wave modes [26], as can be seen that the second wave is of opposite polarization compared to the first compression wave [32, 33], which is consistent with the Gassmann equation [43]. The effects of the internal friction arising from viscous shear stress [8, 36, 37] on propagation of compressional waves are of significant importance.Figure 2
Phase velocity (a), damping factor (b), and phase difference (c) of the first wave versus normalized frequency (Ω=ω/ωc). The damping is the imaginary part of −ω/ci,i=1,2. The phase velocity is the real part of the complex wave velocity. The red line denotes the Darcy-Brinkman law and the black line Darcy’s law.
(a)(b)(c)Figure 3
Phase velocity (a), damping factor (b), and phase difference (c) of the second wave versus normalized frequency (Ω=ω/ωc). The red line presents the Darcy-Brinkman law and the black line Darcy’s law.
(a)(b)(c)Figure 4
Phase velocity (a), damping factor (b), and phase difference (c) of the wave versus normalized frequency (Ω=ω/ωc) under the condition of λ1λ4λ7=λ2λ5λ6. The red line presents the Darcy-Brinkman law and the black line Darcy’s law.
(a)(b)(c)


## 3.3. Wave Fields of Pore Pressure and Stress

Propagation of pore pressure and stress waves is analytically investigated in a one-dimensional semi-infinite domain subject to the step loading,∆pof pore pressure at t=0, x=0. Figure 4 shows a simulation of the frequency dependence of the amplitude of the pore pressure and the stresses over a wave travelling distance. From the comparison of Darcy’s law with the Darcy-Brinkman law, it is found that the Brinkman term is negligible in the near field, while the viscous shear stress significantly influences the wave motion in the far field (see Figure 5(a)). Figure 5(b) shows the interference pattern of the wave fronts. A key result is that during the wave motion, the amplitude of the pore pressure generally decreases, compared to the shear-induced pore pressure changes [44]. The magnitude of the stress firstly increases and then decreases back to its initial value. The amplitude-curve crossovers occur as the interplay between greater phase lags due to increasing frequency [45]. The two-wave modes [43, 46], dissipative and diffusive patterns [47], are demonstrated clearly.Figure 5
Normalized amplitude of pore pressure and stress wave versus travelling distance, (a) comparison of Darcy’s law with the Darcy-Brinkman law for an angular frequency of 1000 Hz, and (b) based on the Darcy-Brinkman effect with angular frequency as a parameter forη=1.0×10−3kgm−1s−1. The solid line denotes the transfer function of the pore pressure, and the dash-dot line represents the transfer function of the stress. The wave travelling distance is normalized by the viscous penetration length δ=2μ/ωρl at frequency of 2×10−6Hz.
(a)(b)Based on the Darcy-Brinkman law, Figure6(a) shows the simulation of the pore pressure and stress waves for various viscosities, η. The viscous dissipation of energy [35, 37] considerably influences the propagation of the pore pressure and stress waves. The pore pressure and stress waves exhibit dissipative and dispersive behaviors [24, 46–48]. The amplitude-curve crossovers are due to the spatial superposition of the confined waves acting like interference patterns [45, 49]. As can be seen in Figure 6(b), the amplitudes of pore pressure and stress are insensitive to the value of the added-mass coefficient. The regime where the pore pressure (or stress) curves coincide implies that the pore pressure (or stress) is in that case completely independent of α∞ and corresponds to the situationu^=w^ as shown from the field Equations (5) and (6). The pressure-coherent stress follows a single gamma distribution, compared to the flow-induced stress in saturated porous media [50].Figure 6
Normalized amplitude of pore pressure and stress versus wave travelling distance based on Darcy-Brinkman law withη as a parameter (a) and with added-mass coefficient α∞ as a parameter (b) at an angular frequency of 1000 Hz. The solid line denotes the transfer function of the pore pressure, and the dash-dot line represents the transfer function of the stress. The wave travelling distance is normalized by the viscous penetration length δ=2μ/ωρl at frequency of 2×10−6Hz.
(a)(b)


## 3.4. Discussions of the Darcy-Brinkman-Biot Formulation

It is noted that the assumption of Darcy’s law is not valid if the frequency exceeds a Biot’s critical frequency [14] above which the compressional wave attenuation is strongly controlled by the viscous shear stress, and the Darcy-Brinkman law should be considered to address the effects of non-Darcy flow on compressional wave motions. Moreover, for compressional elastic wave propagation in Newtonian fluid-saturated porous media, as pointed by Carrillo and Bourg [37], the Darcy term is not adapted to describe systems with both macropores and micropores, which contain two characteristic length scales. Fluid flow in two-scale porous media is governed by a general coupled system of equations, which approximates the Navier-Stokes equations in fluid-filled macropores and resembles the equations for poroelasticity in microporous regions.Now, it has to be realized that in the theoretical analysis, we have assumed that the nonlinear inertial terms in the momentum equations of the fluid and solid can be neglected for the frequency below the critical Biot frequency. This assumption is also reasonable to a small Reynolds number (Re<1) within the microporous region [37]. In the limit of high frequencies, the frequency-dependent interaction forces between solid and fluid can be described in terms of the high-frequency added mass jωρa. However, the Darcy-Brinkman-Biot formulation [37] does not include the Forchheimer equation.Future work involves the addition of Forchheimer term [51, 52], solid dissipation, and viscoelasticity and relaxation effects to the momentum conservation laws of compressional waves in saturated porous media. Of particular interest is that previous mathematical approaches [52, 53] of the Laplace transform, Fourier transform, and matrix transfer would be extended by solving the one-dimensional Darcy-Brinkman-Biot formulation of frequency domain.


## 4. Conclusions

As a conclusion, we show that the propagation of the pore pressure and stress waves in water-saturated elastic porous media with the Darcy-Brinkman law involved. The viscous shear stress significantly influences the phase velocity, damping, and phase lag of the first and second waves. The Darcy-Brinkman law does not lead to the change of the in-phase motion into out-of-phase motion in the Frenkel-Biot wave modes. The second wave is of opposite polarization compared to the first wave. The pore pressure-coherent stress follows a single gamma distribution.Compared to the Darcy-Brinkman law, Darcy’s law overestimates the wave velocity, damping, and phase lag of the first wave, while underestimates the wave velocity, damping, and phase difference of the second wave. The average error overestimated by the Darcy law is 317.01 m/s for the phase velocity, 0.17 m-1 for the damping, and 3.05 (degree) for the phase lag with respect to the first wave. The average error underestimated due to the Darcy law is 79.15 m/s for the phase velocity, 0.36 m-1 for the damping and 3.05 (degree) for the phase lag with respect to the second wave. A characteristic frequency is found to be a decisive condition for the coincidence of the properties of Biot waves of the first and second kinds. The transfer functions of the pore pressure and stress waves are reformulated. The Darcy-Brinkman law describes the dissipative and dispersive behavior of compressional wave propagation in a water-saturated elastic porous medium. The introduction of the Brinkman term yields a better understanding of the Frenkel-Biot waves in saturated porous media. For compressional elastic wave propagation in Newtonian fluid-saturated porous media, when the error of Darcy law is large enough, non-Darcy effects should be considered using the Darcy-Brinkman law.



---

*Source: 1301044-2021-12-27.xml*