Search

We present a numerical method to simulate seismic wave propagation in coupled fluid and porous media. We developed a numerical finite element–based algorithm to approximate solutions to viscoacoustic and Biot's equations, considering the open pore conditions at the interfaces between both media. The algorithm architecture allows to simulate arbitrary distributions of viscoacoustic and poroelastic regions, facilitating the modelling of heterogeneous systems involving complex geometries. The algorithm includes a double parallelization scheme whose efficiency in terms of computing time and memory requirements was tested for different core distributions and mesh sizes. We validate our proposal by performing a comparison between its results and those obtained with a well‐known freely available code. We test its capabilities by studying two different scenarios with geophysical interest: a lake with an irregular bottom and a fractured porous medium. [ABSTRACT FROM AUTHOR]

ABSTRACT We present a generalized effective poroelastic model for periodically layered media in the mesoscopic scale range, which accounts for both Biot's global and interlayer wave-induced fluid flow, as well as for the anisotropy associated with the layering. Correspondingly, it correctly predicts the existence of the fast and slow P-waves as well as quasi and pure S-waves. The proposed analytical model is validated through comparisons of the P-wave and S-wave phase velocity dispersion and attenuation characteristics with those inferred from a one-dimensional numerical solution of Biot's poroelastic equations of motion. We also compare our model with the classical mesoscopic model of White for a range of scenarios. The results demonstrate that accounting for both wave-induced fluid flow mechanisms is essential when Biot's global flow prevails at frequencies that are comparable or smaller with respect to those governing interlayer flow. This is likely to be the case in media of high permeability, such as, for example, unconsolidated sediments, clean sandstones, karstic carbonates, or fractured rocks. Conversely, when interlayer flow occurs at smaller frequencies with respect to Biot's global flow, the predictions of this model are in agreement with White's model, which is based on quasi-static poroelasticity. [ABSTRACT FROM AUTHOR]

The presence of mesoscopic heterogeneities in fluid-saturated porous rocks can produce measurable seismoelectric signals due to wave-induced fluid flow between regions of differing compressibility. The dependence of these signals on the petrophysical and structural characteristics of the probed rock mass remains largely unexplored. In this work, we derive an analytical solution to describe the seismoelectric response of a rock sample, containing a horizontal layer at its centre, that is subjected to an oscillatory compressibility test. We then adapt this general solution to compute the seismoelectric signature of a particular case related to a sample that is permeated by a horizontal fracture located at its centre. Analyses of the general and particular solutions are performed to study the impact of different petrophysical and structural parameters on the seismoelectric response. We find that the amplitude of the seismoelectric signal is directly proportional to the applied stress, to the Skempton coefficient contrast between the host rock and the layer, and to a weighted average of the effective excess charge of the two materials. Our results also demonstrate that the frequency at which the maximum electrical potential amplitude prevails does not depend on the applied stress or the Skempton coefficient contrast. In presence of strong permeability variations, this frequency is rather controlled by the permeability and thickness of the less permeable material. The results of this study thus indicate that seismoelectric measurements can potentially be used to estimate key mechanical and hydraulic rock properties of mesoscopic heterogeneities, such as compressibility, permeability and fracture compliance. [ABSTRACT FROM AUTHOR]

Oscillatory fluid movements in heterogeneous porous rocks induced by seismic waves cause dissipation of wave field energy. The resulting seismic signature depends not only on the rock compressibility distribution, but also on a statistically averaged permeability. This so-called equivalent seismic permeability does not, however, coincide with the respective equivalent flow permeability. While this issue has been analyzed for one-dimensional (1D) media, the corresponding two-dimensional (2D) and three-dimensional (3D) cases remain unexplored. In this work, this topic is analyzed for 2D random medium realizations having strong permeability fluctuations. With this objective, oscillatory compressibility simulations based on the quasi-static poroelasticity equations are performed. Numerical analysis shows that strong permeability fluctuations diminish the magnitude of attenuation and velocity dispersion due to fluid flow, while the frequency range where these effects are significant gets broader. By comparing the acoustic responses obtained using different permeability averages, it is also shown that at very low frequencies the equivalent seismic permeability is similar to the equivalent flow permeability, while for very high frequencies this parameter approaches the arithmetic average of the permeability field. These seemingly generic findings have potentially important implications with regard to the estimation of equivalent flow permeability from seismic data. [ABSTRACT FROM AUTHOR]

The influence of hydraulic conductivity heterogeneity on tide-induced head fluctuations is presented for a theoretical coastal confined aquifer. The conceptual model assumes that the hydraulic conductivity increases linearly with the distance from the coastline. This type of heterogeneity has been observed in many alluvial coastal aquifers. An exact analytical solution that predicts induced head fluctuations is obtained in terms of a Hankel function. The exact solution can be approximated by a simple mathematical expression, valid for small rates of increase of hydraulic conductivity. Both exact and approximate solutions show significant differences from the classical solution obtained for a homogeneous aquifer. Near the coastline the amplitude of the induced head fluctuation is damped but it is enhanced as the distance to the coast increases. The time-lag between sea tide and induced head fluctuation in the aquifer is not linear; it behaves as a square-root type function leading to a faster transmission of the tidal fluctuation. Hypothetical examples show that the influence of hydraulic conductivity heterogeneity can be significant and should be considered for a correct description of the groundwater response. [ABSTRACT FROM AUTHOR]

A patchy saturated two‐layered porous rock, with each layer filled with a different fluid, is examined. We assume that this system is probed by an elastic wave having a wavelength much larger than the layers' thicknesses. We also assume that the diffusion length is smaller than the thickness of an individual layer, which implies hydraulic disconnection between layers. In the context of Gassmann's fluid substitution theory, we analytically derived an exact expression for the effective fluid bulk modulus assuming that both layers have the same porosity, dry frame, and mineral matrix properties. In addition we derived an approximate solution that works well at relatively high porosities. Both solutions are expressed as a weighted average of the arithmetic and harmonic averages of individual bulk moduli of the pore fluid. These weights are explicitly given as functions of the porosity, the fractional thicknesses of both layers, and the elastic moduli of the constituents. For the approximate solution, one does not require explicit knowledge of the shear modulus of the rock. The comparison with laboratory data showed that, in the case where a porous, isotropic rock is filled with water and gas, the approximate solution can be used to model the measured data for high values of water saturation. Key Points: We derive exact and approximate analytic solutions to model the effective fluid bulk modulus of patchy saturated porous rocksWe compare and validate the solutions with stochastic simulation and laboratory data [ABSTRACT FROM AUTHOR]

In glacial studies, properties such as glacier thickness and the basement permeability and porosity are key to understand the hydrological and mechanical behaviour of the system. The seismoelectric method could potentially be used to determine key properties of glacial environments. Here we analytically model the generation of seismic and seismoelectric signals by means of a shear horizontal seismic wave source on top of a glacier overlying a porous basement. Considering a one‐dimensional setting, we compute the seismic waves and the electrokinetically induced electric field. We then analyse the sensitivity of the seismic and electromagnetic data to relevant model parameters, namely depth of the glacier bottom, porosity, permeability, shear modulus and saturating water salinity of the glacier basement. Moreover, we study the possibility of inferring these key parameters from a set of very low noise synthetic data, adopting a Bayesian framework to pay particular attention to the uncertainty of the model parameters mentioned above. We tackle the resolution of the probabilistic inverse problem with two strategies: (1) we compute the marginal posterior distributions of each model parameter solving multidimensional integrals numerically and (2) we use a Markov chain Monte Carlo algorithm to retrieve a collection of model parameters that follows the posterior probability density function of the model parameters, given the synthetic data set. Both methodologies are able to obtain the marginal distributions of the parameters and estimate their mean and standard deviation. The Markov chain Monte Carlo algorithm performs better in terms of numerical stability and number of iterations needed to characterize the distributions. The inversion of seismic data alone is not able to constrain the values of porosity and permeability further than the prior distribution. In turn, the inversion of the electric data alone, and the joint inversion of seismic and electric data are useful to constrain these parameters as well as other glacial system properties. Furthermore, the joint inversion reduces the uncertainty of the model parameters estimates and provides more accurate results. [ABSTRACT FROM AUTHOR]

Most of the analytical solutions to describe tide-induced head fluctuations assume that the coastal aquifer has a constant thickness. These solutions have been applied in many practical problems ignoring possible changes in aquifer thickness, which may lead to wrong estimates of the hydraulic parameters. In this study, a new analytical solution to describe tide-induced head fluctuations in a wedge-shaped coastal aquifer is presented. The proposed model assumes that the aquifer thickness decreases with the distance from the coastline. A closed-form analytical solution is obtained by solving a boundary-value problem with both a separation of variables method and a change of variables method. The analytical solution indicates that wedging significantly enhances the amplitude of the induced heads in the aquifer. However, the effect on time lag is almost negligible, particularly near the coast. The slope factor, which quantifies the degree of heterogeneity of the aquifer, is obtained and analyzed for a number of hypothetical scenarios. The slope factor provides a simple criterion to detect a possible wedging of the coastal aquifer. [ABSTRACT FROM AUTHOR]