When a seismic wave travels through a fluid-saturated porous reservoir containing aligned fractures, it induces oscillatory fluid flow between the fractures and the embedding background medium. Although there are numerous theoretical models for quantifying the associated seismic attenuation and velocity dispersion, they rely on certain assumptions, such as infinitesimal fracture thickness and dilute fracture concentration, which rarely hold in real reservoirs. The objective of this work is to overcome some of these limitations and, therefore, improve the applicability of the available theoretical models. To do so, we extend existing models to the finite fracture thickness case for P-waves propagating perpendicular to the fracture plane using the so-called branching function approach. We consider three types of fractures, namely, periodically and randomly spaced planar fractures, as well as penny-shaped cracks. The extended unified model is then tested by comparing with corresponding numerical simulations based on Biot’s theory of poroelasticity. We consider two cases of 2D rock samples with aligned elliptical fractures, one with low fracture density and the other with high fracture density. The results indicate that the influence of the finite fracture thickness on seismic dispersion and attenuation is small at low frequencies when the fluid pressure has enough time to equilibrate between the fractures and background medium. However, this effect is significant at high frequencies when there is not sufficient time for the fluid pressure equilibration. In addition, the theoretical predictions of the penny-shaped crack model are found to match the numerical simulation results very well, even under relatively high fracture density. Analyses of stress distributions suggest that the small discrepancies found between theoretical predictions and numerical simulations are probably due to fracture interactions. In a companion paper, we will extend the analysis for considering the full stiffness matrix and anisotropic properties of such rocks.

The detection and characterisation of domains of intersecting fractures are important goals in several disciplines of current interest, including exploration and production of unconventional reservoirs, nuclear waste storage, CO2 sequestration, and groundwater hydrology, among others. The objective of this study is to propose a theoretical framework for quantifying the effects of fracture intersections on the frequency‐dependent elastic properties of fluid‐saturated porous and fractured rocks. Three characteristic frequency regimes for fluid pressure communication are identified. In the low‐frequency limit, fractures are in full pressure communication with the embedding porous matrix and with other fractures. Conversely, in the high‐frequency limit, fractures are hydraulically isolated from the matrix and from other fractures. At intermediate frequencies, fractures are hydraulically isolated from the matrix porosity but can be in hydraulic communication with each other, depending on whether fracture sets are intersecting. For each frequency regime, the effective stiffness coefficients are derived using the linear‐slip theory and anisotropic Gassmann equations. Explicit mathematical expressions for the two characteristic frequencies that separate the three frequency regimes are also determined. Theoretical predictions are then applied to two synthetic 2D samples, each containing two orthogonal fracture sets: one with and another without intersections. The resulting stiffness coefficients, Thomsen‐style anisotropy parameters, and the transition frequencies show good agreement with corresponding numerical simulations. The theoretical results are applicable not only to 2D but also to 3D fracture systems and are amenable to being employed in inversion schemes designed to characterise fracture systems.

Seismic waves may exhibit significant dispersion and attenuation in reservoir rocks due to pore-scale fluid flow. Fluid flow at the microscopic scale is referred to as squirt flow and occurs in very compliant pores, such as grain contacts or microcracks, that are connected to other stiffer pores. We have performed 3D numerical simulations of squirt flow using a finite-element approach. Our 3D numerical models consist of a pore space embedded into a solid grain material. The pore space is represented by a flat cylinder (a compliant crack) whose edge is connected with a torus (a stiff pore). Grains are described as a linear isotropic elastic material, whereas the fluid phase is described by the quasistatic linearized compressible Navier-Stokes momentum equation. We obtain the frequency-dependent effective stiffness of a porous medium and calculate dispersion and attenuation due to fluid flow from a compliant crack to a stiff pore. We compare our numerical results against a published analytical solution for squirt flow and analyze the effects of its assumptions. Previous interpretation of the squirt flow phenomenon based mainly on analytical solutions is verified, and some new physical effects are identified. The numerical and analytical solutions agree only for the simplest model in which the edge of the crack is subjected to zero fluid pressure boundary condition while the stiff pore is absent. For the more realistic model that includes the stiff pore, significant discrepancies are observed. We identify two important aspects that need improvement in the analytical solution: the calculation of the frame stiffness moduli and the frequency dependence of attenuation and dispersion at intermediate frequencies.

Numerous theoretical models have been proposed for computing seismic wave dispersion and attenuation in rocks with aligned fractures due to wave-induced fluid flow between the fractures and the embedding background. However, all these models rely on certain assumptions, for example, infinitesimal fracture thickness or dilute fracture concentration, which rarely hold in real reservoirs and, thus, limit their applicability. To alleviate this issue, theoretical models for periodically or randomly spaced planar fractures and penny-shaped cracks were recently extended by the authors to the case of finite fracture thickness for P-waves propagating perpendicular to the fracture plane. Theoretical predictions under low and relatively high fracture density were then assessed by comparing with corresponding numerical simulations. However, the case of arbitrary incidence angles as well as the behaviors of S-waves remained unexplored. In this work, we therefore extended the prediction results to the full stiffness matrix through two theoretical approaches. The first approach uses an interpolation between the low- and high-frequency limits using a relaxation function obtained from the normal-incidence solution. The second approach is based on the linear slip theory with a frequency-dependent fracture compliance. Both derivations rely on the fact that all the stiffness coefficients are controlled by the same relaxation function. With the full stiffness matrix, anisotropic seismic properties can then be studied. P- and S-wave velocities and attenuations at different frequencies and incidence angles and also corresponding anisotropy parameters are calculated for one synthetic 2D rock sample. The results indicate that the predictions provided by the two theoretical approaches are in good agreement with each other and also indicate a good agreement with the corresponding numerical simulations. The extended theoretical models presented in this work are easy to apply and computationally much cheaper than numerical simulations and, hence, can be used in the seismic characterization of fractured reservoirs.

Estimating the impact of solid pore fill on effective elastic properties of rocks is important for a number of applications such as seismic monitoring of production of heavy oil or gas hydrates. We develop a simple model relating effective seismic properties of a rock saturated with a liquid, solid, or viscoelastic pore fill, which is assumed to be much softer than the constituent minerals. A key feature of the model is division of porosity into stiff matrix pores and compliant crack‐like pores because the presence of a solid material in thin voids stiffens the rock to a much greater extent than its presence in stiff pores. We approximate a typical compliant pore as a plane circular interlayer surrounded by empty pores. The effect of saturation of the stiff pores is then taken into account using generalized Gassmann's equations. The proposed model provides a good fit to measurements of the shear stiffness and loss factor of the Uvalde heavy‐oil rock at different temperatures and frequencies. When the pore fill is solid, the predictions of the scheme are close to the predictions of the solid squirt model recently proposed by Saxena and Mavko. At the same time, the present scheme also gives a continuous transition to the classic Gassmann's equations for a liquid pore fill at low frequencies and the squirt theory at high frequencies.

No abstract

Vertical seismic profile (VSP) is one of the technologies for monitoring hydrocarbon production and CO2 geosequestration. However, quantitative interpretation of time‐lapse VSP is challenging due to its irregular distribution of source‐receiver offsets. One way to overcome this challenge is to use full waveform inversion (FWI), which does not require regular offsets. We present a workflow of elastic FWI applied to offset vertical seismic profile data for the purpose of identification and estimation of time‐lapse changes introduced by injection of 15,000 t of CO2‐rich gas mixture at 1.5 km depth. Application of this workflow to both synthetic and field data shows that elastic FWI is able to detect and quantify the time‐lapse anomaly in P wave velocity with the magnitude of 100–150 m/s.

No abstract