Poroelastic Theory Research Articles

Emergent vegetation flow with varying vertical porosity

The vertical distribution of an emergent vegetation in wetlands and marshes is generally non-uniform due to the variations of the stem thickness and the leaf density. This study focuses on the effect of the varying vertical porosity on the hydrodynamic characteristics of the emergent vegetation flow. The new governing equation for the emergent vegetation flow with varying vertical porosity is established by applying the poroelastic media flow theory. According to the movement mechanism of the water flow in the porous media and the experimental data, the best fitted expression of the permeability is established for the emergent vegetation flow. The velocity distribution is obtained with the finite analytic method. Then, a flume experiment is performed using truncated cones to simulate the actual vegetation with varying vertical porosity. The calculated velocity distribution is compared with the measured data, which shows that the theoretical results are in good agreement with the experimental data. Finally, the influence of the vertical porosity distribution on the vegetation flow characteristics is analyzed graphically. The study can provide a useful reference and technical support for further study of the flow with a complex shape vegetation.

Biochemomechanical modeling of vascular collapse in growing tumors

The collapse of blood vessels are widely observed in solid tumors, but the mechanisms underpinning this abnormal behavior remain unclear. In this paper, we investigate the stability of blood vessels embedded in a growing solid tumor by using a chemomechanical poroelastic theory. Linear stability analysis is first made to give the critical condition of vascular buckling. Thereby, pattern diagrams are provided to predict the collapsed shape of a blood vessel from the mechanical and geometric parameters of the system. Due to the mechano-chemo-biological coupling effect in tumors which involve both fluid transport and tissue growth, the buckling of blood vessels exhibits distinctly different features from such previously reported tubular tissues as airways and esophagi. We show that interstitial fluid pressure tends to drive blood vessels to buckle with a lower buckling mode, while perivascular cell proliferation favors a higher mode. Furthermore, a nonlinear biochemomechanical finite element method is formulated to track the morphological evolution of blood vessels during postbuckling. It is found that the blood vessels with the second mode of buckling may readily evolve into a crack-like shape. This method can reproduce the essential features of vascular collapse observed in in vivo tumors and provide clues for vascular normalization and stress alleviation in cancer treatments.

Mechanistic analysis of coal permeability evolution data under stress-controlled conditions

At present, two types of experiments under stress-controlled conditions were normally conducted to measure coal permeability: constant confining pressure (CCP) tests and constant effective stress (CES) ones. The original rationale of this situation was to assume that the impacts of effective stresses and gas sorption-induced matrix swelling/shrinking on coal permeability could be separated and investigated individually. In this study, we collected coal permeability data measured under both conditions with a purpose to see if this original rationale was appropriate. This goal was achieved through collection of experimental permeability data under the CCP conditions; collection of experimental permeability data under the CES conditions; and comparison of those experimental data with solutions of the poroelastic theory. For CCP tests, the permeability ratios change from reductions (less than 1.0) to enhancements (greater than 1). These changes are bounded by an upper envelope and a lower one. The upper envelope is corresponding to the solution of free-swelling while the lower one zero-swelling. For CES tests, the permeability ratios also change within an upper envelope and a lower one. The upper envelope is equal to 1.0 corresponding to the solution of free-swelling while the lower one zero-swelling. Through these comparisons, we found that permeability data for both types of tests are confined within the poroelastic solutions for two extreme boundary conditions: free-swelling and zero-swelling. These findings suggest that permeability ratios for both constant confining tests and constant effective stress tests are primarily determined by the matrix-fracture interactions, including sorption-induced swelling/shrinking, through transient effective stresses in matrixes and fractures.

Quantifying Compressible Groundwater Storage by Combining Cross‐Hole Seismic Surveys and Head Response to Atmospheric Tides

Groundwater specific storage varies by orders of magnitude, is difficult to quantify, and prone to significant uncertainty. Estimating specific storage using aquifer testing is hampered by the nonuniqueness in the inversion of head data and the assumptions of the underlying conceptual model. We revisit confined poroelastic theory and reveal that the uniaxial specific storage can be calculated mainly from undrained poroelastic properties, namely, uniaxial bulk modulus, loading efficiency, and the Biot‐Willis coefficient. In addition, literature estimates of the solid grain compressibility enables quantification of subsurface poroelastic parameters using field techniques such as cross‐hole seismic surveys and loading efficiency from the groundwater responses to atmospheric tides. We quantify and compare specific storage depth profiles for two field sites, one with deep aeolian sands and another with smectitic clays. Our new results require bulk density and agree well when compared to previous approaches that rely on porosity estimates. While water in clays responds to stress, detailed sediment characterization from a core illustrates that the majority of water is adsorbed onto minerals leaving only a small fraction free to drain. This, in conjunction with a thorough analysis using our new method, demonstrates that specific storage has a physical upper limit of m−1. Consequently, if larger values are derived using aquifer hydraulic testing, then the conceptual model that has been used needs reappraisal. Our method can be used to improve confined groundwater storage estimates and refine the conceptual models used to interpret hydraulic aquifer tests.

Three-Dimensional Numerical Investigation of Coupled Flow-Stress-Damage Failure Process in Heterogeneous Poroelastic Rocks

The failure mechanism of heterogeneous rocks (geological materials), especially under hydraulic conditions, is important in geological engineering. The coupled mechanism of flow-stress-damage should be determined for the stability of rock mass engineering under triaxial stress states. Based on poroelasticity and damage theory, a three-dimensional coupled model of the flow-stress-damage failure process is studied, focusing mainly on the coupled characteristics of permeability evolution and damage in nonhomogeneous rocks. The influences of numerous mesoscale mechanical and hydraulic properties, including homogeneity, residual strength coefficient, loading rates, and strength criteria, on the macro mechanical response are analyzed. Results reveal that the stress sensitive factor and damage coefficient are key variables for controlling the progress of permeability evolution, and these can reflect the hydraulic properties under pre-peak and post-peak separately. Moreover, several experiments are conducted to evaluate the method in terms of permeability evolution and failure process and to verify the proposed two-stage permeability evolution model. This model can be used to illustrate the failure mechanics under hydraulic conditions and match different rock types. The relation of permeability with strain can also help confirm appropriate rock mass hydraulic parameters, thereby enhancing our understanding of the coupled failure mechanism in rock mass engineering.

Poromechanical Properties of a Sandstone Under Different Stress States

Poromechanical properties of a sandstone from an underground gas storage site are investigated with the use of a neutral gas to control the pore pressure. The poroelastic theory of anisotropic medium is used to evaluate elastic properties and coupling coefficients under different stress states. In-situ CT observations have been made under uniaxial compression test and they are used to underline some effects due to the occurrence of damage and microcraks. Under hydrostatic loading, pores and micro-cracks are gradually compressed. The initial state of the sandstone is “slightly” transversely isotropic but in a first approach the Biot’s tensor reduces to a scalar. This coefficient decreases with the increase in confining pressure that can be attributed to the closure of micro-cracks. The effect of damage and cracking is investigated with two series of conventional triaxial tests conducted at different confining pressures, either with an increase in pore pressure (series 1) or with a decrease in pore pressure (series 2), to evaluate the coupling coefficients. The results, derived from these two test series, are consistent as regards the effects of cracking on the material behavior. There is an obvious damage due to the triaxial loading which induces a parallel decrease in the ratio E3/(1 − ν3) and in the modulus H3. Axial compaction results in the continual increase in the H1 modulus observed for every sample whatever the measurement technique was.

A multi-layered poroelastic slab model under cyclic loading for a single osteon

BackgroundAn osteon consists of a multi-layered bone matrix and interstitial fluid flow in the lacunar–canalicular system. Loading-induced interstitial fluid flow in the lacunar–canalicular system is critical for osteocyte mechanotransduction and bone remodelling.MethodsTo investigate the effects of the lamellar structure and heterogeneous material properties of the osteon on the distributions of interstitial fluid flow and seepage velocity, an osteon is idealized as a hollow two-dimensional poroelastic multi-layered slab model subjected to cyclic loading. Based on poroelastic theory, the analytical solutions of interstitial fluid pressure and seepage velocity in lacunar–canalicular pores were obtained.ResultsThe results show that strain magnitude has a greater influence on interstitial fluid pressure than loading frequency. Interestingly, the heterogeneous distribution of permeability produces remarkable variations in interstitial fluid pressure and seepage velocity in the cross-section of cortical bone. In addition, interstitial fluid flow stimuli to osteocytes are mostly controlled by the value of permeability at the surface of the osteon rather than at the inner wall of the osteon.ConclusionInterstitial fluid flow induced by cycling loading stimuli to an osteocyte housed in a lacunar–canalicular pore is not only correlated with strain amplitude and loading frequency, but also closely correlated with the spatial gradient distribution of permeability. This model can help us better understand the fluid flow stimuli to osteocytes during bone remodelling.

Permeability mapping of gelatin methacryloyl hydrogels

We report the development of an efficient, customized spherical indentation-based testing method to systematically estimate the hydraulic permeability of gelatin methacryloyl (GelMA) hydrogels fabricated in a wide range of mass concentrations and photocrosslinking conditions. Numerical simulations and Biot’s theory of poroelasticity were implemented to calibrate our experimental data. We correlated elastic moduli and permeability coefficients with different GelMA concentrations and crosslinking densities. Our model could also predict drug release rates from the GelMA hydrogels and diffusion of biomolecules into the three-dimensional GelMA hydrogels. The results potentially provide a design map for choosing desired GelMA-based hydrogels for use in drug delivery, tissue engineering, and regenerative medicine, which may be further expanded to predicting the permeability behaviors of various other hydrogel types. Statement of SignificanceGelMA hydrogels have attracted increasing attention in recent years as matrices for cell cultures and biomolecule delivery. This inexpensive polymer is derived from gelatin functionalized with methacryloyl groups that can be crosslinked by photochemical reactions. Here we report the development of an efficient, customized testing method to systematically estimate the hydraulic permeability of GelMA hydrogels. Hydraulic permeability indicates the resistance of GelMA hydrogels to the movement of saturated fluid. We used the model to measure the elastic moduli and permeability coefficients, providing a permeability map for various GelMA hydrogel formulations.

Analysis of radial vibrations in thick walled hollow poroelastic cylinder in the framework of Biot’s extension theory

Purpose In this paper, wave propagation in a poroelastic thick-walled hollow cylinder is investigated in the framework of Biot’s extension theory. Biot’s theory of poroelasticity is valid for isotropic porous solids saturated with non-viscous fluid. The bulk and shear viscosities are not considered in the classical Biot’s theory. Biot’s extension theory takes all these into an account. Biot’s extension theory is applied here to investigate the radial vibrations in thick-walled hollow poroelastic cylinder. The paper aims to discuss these issues. Design/methodology/approach By considering the stress-free boundaries, the frequency equation is obtained in the presence of dissipation. Limiting case when the ratio between thickness and inner radius is very small is investigated numerically. In the limiting case, the asymptotic expansions of Bessel functions are employed so that frequency equation is separated into two parts which gives attenuation coefficient and phase velocity. If the shear viscosity is neglected, then the problem reduces to that of the classical Biot’s theory. Findings For the numerical purpose, the solids Berea sandstone and bone are used. The results are presented graphically. Originality/value Radial vibrations of thick-walled hollow poroelastic cylinder are investigated in the framework of Biot’s extension theory. Due to the mathematical complexity, limiting case is considered. The complex valued frequency equation is discussed numerically which gives the attenuation coefficient and phase velocity. If shear viscosity is neglected, then the problem reduces to that of the classical Biot’s theory. The comparison has been made between the current results and that of classical results.

Poroelastic modeling of cutting rock in pressurized condition

Cutting rock in pressurized conditions is essentially different as compared to the ambient conditions. Under pressurized condition, cutter is not only fragmenting the rock matrix, but also driving pore fluid to flow forward. Due to the coupling of stresses and pore pressure, the pore pressure change caused by cutter will affect effective stresses in the rock and eventually affect mechanical specific energy (MSE). The change in stresses and pore pressure caused by cutter movement are the poroelastic effects in cutting process.To study pressurized cutting process, we developed a cutting model based on the theory of linear poroelasticity. The model can predict coupled stresses and pore pressure in the rock and give a better understanding of the poroelastic effect in the cutting process. To obtain the desired solution, a Fourier Transform is employed. After the analytic solution in Fourier space is obtained, a Discrete Fourier Transform is employed to numerically invert the solution back to the real space. The Mohr-Coulomb failure criterion is also introduced to determine failure of the rock ahead of the cutter. Eventually, the model can predict cutting forces and MSE in pressurized cutting process. The model is also compared to the published experimental data and a good consistency between model and experiment results was found.The influence of hydrostatic pressure and the existence of cuttings on MSE are studied and poroelastic effect during cutting process is also discussed. The results in this paper can be applied to bit-rock interaction in well drilling.

Modeling of Biot’s coefficient for a clay-bearing sandstone reservoir

The exploration of unconventional reservoirs has been vigorously developed with the application of hydraulic fracturing technology in horizontal wells. The success of fracturing depends on the accuracy of the effective stress, which could be considerably influenced by Biot’s coefficient. Biot’s coefficient is mainly influenced by the clay content, porosity, and fluid saturation. Previous studies regarding the calculation of Biot’s coefficient considered only the porosity or porosity structure and are irrational for sandstones with a high clay content and low permeability. A method was conducted to estimate Biot’s coefficient from wave velocities based on the Biot and Squirt (BISQ) poroelasticity theory. It is suggested that the clay content surrounding the skeleton, rather than that in the skeleton, affects the Biot’s coefficient significantly. A clay content coefficient and fluid saturation are considered in the BISQ model. Experiments with samples from four oil fields show that the Biot’s coefficient calculated using our method matches well with those obtained from in situ experiments. It is suggested that the new method for determining Biot’s coefficient is much more suitable for sandstones with low permeability, low porosity, and high clay content, compared to the methods put forward by other studies. In addition, it is proposed to calculate Biot’s coefficient by using seismic or acoustic logging data.

Estimation of fluid indicator and dry fracture compliances using azimuthal seismic reflection data in a gas-saturated fractured reservoir

Amplitude variation with offset and azimuth (AVOA) inversion is a powerful tool used to predict the fluid and fracture information by utilizing the observable seismic reflection amplitudes. Understanding the effect of in situ fluid and fracture characteristics on seismic response is significant in a fractured reservoir. Our goal is to demonstrate an AVOA inversion approach to utilize the observable azimuthal seismic reflection data to predict the fluid and fracture properties in a gas-saturated fractured reservoir. Starting from Gassmann’s anisotropic poroelasticity theory and the linear-slip model, we first derive a linearized PP-wave reflection coefficient in terms of fluid, rigidity, density, and dry fracture compliance parameters based on the scattering function and perturbations in dry stiffness matrix of a horizontal transversely isotropic (HTI) rock containing a system of aligned vertical fractures. The derived equation builds a linearized relationship between the characteristics of fluid and fracture and the seismic response of azimuthal seismic reflection amplitudes. Incorporating the convolution model and Bayesian inversion, we propose a feasible Bayesian AVOA inversion method to estimate the fluid and fracture parameters directly. Cauchy and Gaussian probability distribution are used for the a priori information of model parameters and the likelihood function, respectively, and the nonlinear iteratively reweighted least squares (IRLS) strategy is finally utilized to solve the maximum a posteriori solutions of model parameters. The synthetic examples containing a moderate noise demonstrate the feasibility of inversion method, and the real data illustrates the inversion stabilities of fluid and fracture parameters in a gas-saturated fractured reservoir.

A phase-field modeling approach of fracture propagation in poroelastic media

This paper proposes a phase field model for fracture in poroelastic media. The porous medium is modeled based on the classical Biot poroelasticity theory and the fracture behavior is controlled by the phase field model. Moreover, the fracture propagation is driven by the elastic energy where the phase field is used as an interpolation function to transit fluid property from the intact medium to the fully broken one. We use a segregated (staggered) scheme and implement our approach in Comsol Multiphysics. The proposed model is verified by a single-phase solid subjected to tension and a 2D specimen subjected to an increasing internal pressure. We also compare our results with analytical solutions. Finally, we show 2D and 3D examples of internal fluid injection to illustrate the capability of the proposed approach.

Velocity and attenuation of elastic wave in a developed layer with the initial inner percolation in the pores

An improved model based on Biot poroelastic theory is presented incorporating the initial seepage flow in the matrix. The proposed model quantifies wave propagation in the oil field development process as VSP, transient well tests, and seismic production technology. Porosity variation, fluid–solid compressibility, and pseudo-threshold pressure gradient in low permeability reservoir are simultaneously considered, resulting in a modified form of continuity and motion equations. Integrating the equations and decomposition for the fluid–solid displacements helps to derive the analytical wave vectors of the fast P, slow P, and S waves in porous media. The sensitivity of affecting factors, such as petrophysical parameter, fluid property, and vibration frequency, on the wave velocity and attenuation is subsequently evaluated. Furthermore, the Biot model can be considered a special situation when the initial flow rate in this model tends to zero. The increase in initial percolating rate, expressed by the ratio of initial flow rate and solid scalar potential, causes an obvious decrease in the fast P wave velocity and an increase in its quality factor. Propagation of the slow P wave and the S wave is scarcely influenced. Low vibration frequency and low permeability also contribute to a large difference between the wave propagation parameters of the improved and Biot models. Cognition of the proportion in the inertia and viscosity effects is helpful in analyzing the complicated multiple mechanisms in wave propagation in the actual developed layer.

Dynamic characteristics of a large‐diameter pile in saturated soil and its application

SummaryThis study theoretically investigates the dynamic response of an end‐bearing pile embedded in saturated soil considering the transverse inertial effect of the pile. The saturated soil surrounding the pile is described by Biot poroelastic theory, and the pile is represented by a Rayleigh‐Love rod because both the vertical and radial displacements at the soil‐pile interface are considered. The potential function decomposition method and variable separation method are introduced to solve the governing equations of the soil, in which the vertical and radial displacement components are coupled. The governing equation of the pile is solved using the continuity conditions at the pile‐soil interface. Next, the velocity admittance in the frequency domain and the velocity response in the time domain at the pile top are presented based on the Laplace transform and inverse Fourier transform, respectively. Subsequently, the reduced solution is compared with a 1‐dimensional model solution to verify the validity, and the influences of the slenderness ratio of the pile on the transverse inertial effect of the pile are analyzed. Moreover, Poisson ratio, the slenderness ratio of the pile, and the pile‐soil modulus ratio are studied. Finally, the theoretical and measured curves in the engineering project are compared, and the results demonstrate the good application prospects of the solution presented in this article.

Boundary-element modelling of dynamics in external poroviscoelastic problems

A problem of a spherical cavity in porous media is considered. Porous media are assumed to be isotropic poroelastic or isotropic poroviscoelastic. The poroviscoelastic formulation is treated as a combination of Biot’s theory of poroelasticity and elastic-viscoelastic correspondence principle. Such viscoelastic models as Kelvin–Voigt, Standard linear solid, and a model with weakly singular kernel are considered. Boundary field study is employed with the help of the boundary element method. The direct approach is applied. The numerical scheme is based on the collocation method, regularized boundary integral equation, and Radau stepped scheme.

The state of stress in SW Iran and implications for hydraulic fracturing of a naturally fractured carbonate reservoir

The first part of this study provides an interpretation of the present-day stress state and mechanical parameters for the giant Bangestan reservoir in SW Iran. A total of 28 drilling induced tensile fractures (DITF) and 172 breakouts were mapped in seven wells in the Bangestan reservoir using image logs. The average maximum principal horizontal stress (SH) direction for each well is given a quality ranking according to the World Stress Map (WSM) criteria. The mean orientation of SH finds to be N44°(±9.95°). A correlation is seen between this data and stress orientations from earthquake focal mechanism solution. This suggests that stresses are related to tectonic forces generated by Arabia–Eurasia collision. Estimates of vertical principal stress (Sv) gradients vary between 23.3 and 24.7 kPa/m (1.03–1.09 psi/ft), based on empirical methods and density log data, respectively. Minimum horizontal stress (Sh) gradient estimates vary from 15.2 to 17.4 kPa/m (0.67–0.77 psi/ft) using extended leak-off test (XLOT) and step-rate test results, respectively. Uncertainties remain in estimating SH from acid fracturing in the deep naturally fractured carbonate reservoirs, even though numerous papers have been written on the subject. Herein, an approach called Deformation/Diffusion/Thermal (DDT) is suggested to estimate SH in the Bangestan Reservoir. The method combines poroelastic theories for the circular underground cavities, bilinear theory for fluid flow through the preexisting induced fracture and thermal stress. Using this method, SH gradient is calculated to be 18.9 kPa/m (0.83 psi/ft). For comparison, horizontal stress magnitudes and mechanical parameters were computed from sonic wave velocity using poroelastic-tectonic stress equations over the full reservoir layer thickness. All data have been calibrated with field tests and lab results. The measured SH was in good agreement with the previous analytical solution. The study suggests that the contemporary state of stress in the reservoir is characterized by a normal faulting (NF) stress regime. Next, current and prospective hydraulic fracturing operations in the field are discussed. Based on the NF stress regime, a vertical plane fracture is expected and the wells drilled in the direction of SH are likely to intersect the most natural fractures.

Analysis of Dynamic Fracture Compliance Based on Poroelastic Theory. Part II: Results of Numerical and Experimental Tests

In Part I, a dynamic fracture compliance model (DFCM) was derived based on the poroelastic theory. The normal compliance of fractures is frequency-dependent and closely associated with the connectivity of porous media. In this paper, we first compare the DFCM with previous fractured media theories in the literature in a full frequency range. Furthermore, experimental tests are performed on synthetic rock specimens, and the DFCM is compared with the experimental data in the ultrasonic frequency band. Synthetic rock specimens saturated with water have more realistic mineral compositions and pore structures relative to previous works in comparison with natural reservoir rocks. The fracture/pore geometrical and physical parameters can be controlled to replicate approximately those of natural rocks. P- and S-wave anisotropy characteristics with different fracture and pore properties are calculated and numerical results are compared with experimental data. Although the measurement frequency is relatively high, the results of DFCM are appropriate for explaining the experimental data. The characteristic frequency of fluid pressure equilibration calculated based on the specimen parameters is not substantially less than the measurement frequency. In the dynamic fracture model, the wave-induced fluid flow behavior is an important factor for the fracture–wave interaction process, which differs from the models at the high-frequency limits, for instance, Hudson’s un-relaxed model.

Impact of gas adsorption on apparent permeability of shale fracture and shale gas recovery rate

Shale gas reservoirs can experience adsorption-induced strain due to the in-situ shale gas and the CO2 invasion after anhydrous fracturing. The adsorption-induced strain as well as effective stress and flow regime affect apparent permeability of fractures in shale reservoirs, which makes the impact of adsorption-induced on apparent permeability is difficult to be analyzed and normally ignored in models. In this study, based on poroelastic theory and characteristics of adsorption-induced strain of shale, the analytical apparent permeability model for fractured shale was built to analyze experimental data obtained from various boundary conditions, meanwhile, a experimental method was proposed to investigate the impact of gas adsorption on apparent permeability. With this method, the impact of gas adsorption on organic-rich shale was identified. Finally, the impact of gas adsorption on recovery rate of shale gas reservoirs was analyzed through combining the apparent permeability model and multiscale model of gas flow through fractured media. The analytical model illustrates that the controlling mechanisms of apparent permeability evolutions for different boundary conditions and gases are significantly different. Based on this feature, using the experimental data obtained from non-adsorption phase gas and adsorption-phase gas under conditions of constant confining pressure and constant pore pressure, the impact of internal strain due to adsorption on apparent permeability was identified and cannot be ignored. In additional, compared with the internal strain due to adsorption, the adsorption-induced strain of the whole shale plays a more important role in evolutions of apparent permeability and recovery rate of shale gas reservoirs.

Sound Transmission Loss Enhancement in an Inorganic‐Organic Laminated Wall Panel Using Multifunctional Low‐Density Nanoporous Polyurea Aerogels: Experiment and Modeling

Recently, the authors have reported an exceptional normal incidence sound transmission loss characteristic for a class of low density, highly porous, and mechanically strong polyurea aerogels. Herein, a laminated composite comprising the organic low‐density aerogels bonded with an inorganic compound (e.g., gypsum materials) is considered to investigate the constrained damping effects of the aerogels on the airborne sound insulation behavior of the composite using the standard chamber‐based diffuse sound field measurements. Huge improvement in the sound transmission loss is obtained due to the use of aerogel without a significant increase in the overall weight and thickness of the composite panel (e.g., more than 10 dB increase by reaching 40 dB sound transmission loss at 2 kHz after the implementation of only two 5 mm‐thick aerogel layers at bulk densities 0.15 and 0.25 g cm−3). This uncommon behavior breaks the empirical “Mass Law” nature of the most conventional acoustic materials. In addition, an exact analytical time‐harmonic plane‐strain solution for the diffused wave propagation through the multilayered structure is provided using theories of linear elasticity and Biot's dynamic poroelasticity. The theoretical results are well supported by the experiments which can be utilized for the design of the future light‐weight multifunctional composite structures.