Abstract:The propagation of seismic waves in fractured rocks is greatly influenced by the fracture system and fluid content. Seismic anisotropy varies with frequency and fluid saturation. Previous theories on frequency-dependent anisotropy (FDA) are mostly limited to the assumption of single-phase fluid, whereas almost all reservoirs are usually partially or fully saturated with one or more fluids. The inelasticity, heterogeneity, and anisotropy exhibited by real earth media seriously challenge the traditional theory of uniform, fully elastic media. Therefore, based on the Chapman elastic theory, this paper considers the viscoelastic model, the relative mobility of saturated fluids, and the coupled effects of commonly independently considered squirt flow and patch effects on the frequency-dependent anisotropic seismic wave propagation in viscoelastic media. By calculating the frequency-dependent anisotropic elastic coefficients of two-phase immiscible fluid-saturated fractured rocks, a new viscoelastic Chapman extension model for FDA in partially saturated rocks, which includes the unified seismic wave propagation effects of squirt flow and patch effects, is constructed. The effect of relative permeability is significant, and compared to the case of complete saturation, the effective fluid mobility in partially saturated rocks may be lower, which may lead to nonmonotonic changes in elastic modulus and water saturation. Numerical experiments are conducted on the proposed viscoelastic Chapman extension model to discuss the coupling effects of squirt flow and patch effects on the FDA of viscoelastic media in cracked rocks with partially saturated reservoirs, under the conditions of no fractures and the existence of fractures. The simulation results of the P-wave modulus, P-wave velocity, and anisotropy parameters confirm the validity of the proposed method and model. This paper links the anisotropy of viscoelastic media with fluid flow parameters in fractures, which is conducive to improving the understanding of frequency-dependent seismic anisotropy in partially saturated rocks with fractures and the combination of seismology and reservoir engineering.

Forward Modeling of Seismic Wave Fields in Vertical Fractured Strata

Finite Difference Modelling for Vertical Seismic Profile in Elastic Media

2D Frequency-Domain Finite-Difference Acoustic Wave Modeling Using Optimized Perfectly Matched Layers

Borehole and crosshole seismic wave field forward modeling

Waveform Imaging Based on Linear Forward Representations for Scalar Wave Seismic Data

3-D Simulation of Wave Propagation in the Fractured Porous Media for Subsurface Sensing: A Rotated Staggered Finite-Difference Algorithm Based on Equivalent Medium Theory

Modeling seismic wave propagation within complex structures

Forward Modeling of Two-phase Anisotropic Media

Modeling Of Rayleigh-Wave Propagation In Three-Dimensional Media With Topography Using Staggered-Grid Finite Difference Scheme

Simulate the elastic wavefields in media with an irregular surface topography based on staggered grid finite difference

Research on the forward-modeling method of viscoelastic Chapman extension in partially saturated rocks

Numerical simulation of cross-hole seismic wave field

Seismic wave-field modeling of 3D irregularities in layered half-space under P-, S-, and Rayleigh waves with arbitrary incident azimuth

A finite-element algorithm with a perfectly matched layer boundary condition for seismic modelling in a diffusive-viscous medium

Seismic wave modeling of fluid-saturated fractured porous rock: including fluid pressure diffusion effects of discretely distributed large-scale fractures

Perfectly matched layers simulation in 2-D VTI media

Spectral Element Modeling of Elastic Wave Propagation in an Anisotropic Background with Discrete Anisotropic Fractures

Theoretical modelling of seismic dispersion, attenuation and frequency-dependent anisotropy in a fluid-saturated porous rock with intersecting fractures

3-D Forward Modeling of Transient EM Field in Rough Media Using Implicit Time-Domain Finite-Element Method.

Modeling 1-D elastic P-waves in a fractured rock with hyperbolic jump conditions

A benchmark study for quasi-static numerical upscaling of seismic wave attenuation and dispersion in fractured poroelastic rocks