Search

Geometric heterogeneities in tight reservoir rocks saturated with a fluid mixture may exhibit different scale distribution characteristics. Conventional models of rock physics based on poroelasticity, which usually consider single‐scale pore structure and fluid patches, are inadequate for describing elastic wave responses. A major challenge is to establish the relationship between the wave response at different spatial scales and frequencies. To address this problem, three sets of observational data over a wide frequency range were obtained from a tight oil reservoir in the Ordos Basin, China. Ultrasonic measurements were made on eight sandstone samples at partial oil‐water saturation at 0.55 MHz. Data from six borehole measurements and seismic profiles were acquired and analyzed at about 10 kHz and 30 Hz, respectively. Analysis of the cast thin sections shows that dissolution pores and microcracks generally develop, with fractal dimensions of the pores ranging from 2.45 to 2.67 for the samples with porosities between 5.1% and 10.2%. Compressional wave velocity and attenuation were estimated from the observed data. The results show that the velocity dispersion from seismic to ultrasonic frequencies is 10.02%, mostly occurring between sonic and ultrasonic frequencies. The attenuation is stronger at higher oil saturation. The relationships between velocity, attenuation, and wavelength were established and can be used for further forward modeling and seismic interpretation studies. A partial saturation model has been derived based on effective differential medium theory and a double double‐porosity model, assuming that the medium contains fractal cracks and fluid patches. The effects of scale and saturation on wave responses are prevalent. Modeling results consistent with observed data show that the radii of cracks and fluid patches range from 0.1 μm to 2.8 mm, affecting ultrasonic, acoustic, and seismic attenuation. The multiscale data and proposed model quantify the relationship between fracture and fluid distributions and attenuation and could be useful for upscaling to the reservoir scale. The study helps improve the understanding of seismic wave propagation in partially saturated rocks, which has potential applications in seismic exploration, hydrocarbon production in reservoirs, and CO2 sequestration in aquifers. Plain Language Summary: The physical properties of the rock and fluid can be inferred from the measured elastic wave responses and energy dissipation characteristics. However, the effects of heterogeneities of different sizes and at different frequencies can hinder studies to quantify wave responses in a partially saturated porous medium, which are usually based on laboratory measurements. A major problem is the difference between observed frequencies and scales: megahertz in the laboratory, 10 of kilohertz in the borehole scale, and hertz in the seismic exploration scale. In this work, the frequency‐ and saturation‐dependent compressional velocity and attenuation are investigated using three geophysical data sets from the same tight reservoirs. A strong velocity dispersion over the measured frequency range is observed. The stronger attenuation at partial saturation may be caused by the multiscale heterogeneities of the pore structure and fluid patch distribution. A fractal poroelasticity model is developed by gradually inserting inclusions of different sizes with compliant pores and liquid patches into a homogeneous host skeleton. The wave responses are significantly affected by scale distribution and saturation. The proposed model, verified by the measured data, can be useful in interpreting the anelasticity of tight heterogeneous reservoirs in a broadband range. Key Points: We investigate how fluid saturation affects wave attenuation and dispersion over a wide frequency range in tight sandstone layersA partial saturation model describes the fractal properties of the layersThe size of the cracks and fluid patches are determined by applying the fractal model to the multiscale data [ABSTRACT FROM AUTHOR]

Wave propagation in anelastic rocks is a relevant scientific topic in basic research with applications in exploration geophysics. The classical Biot theory laid the foundation for wave propagation in porous media composed of a solid frame and a saturating fluid, whose constitutive relations are linear. However, reservoir rocks may have a high-viscosity fluid, which exhibits a non-Newtonian (nN) behaviour. We develop a poroelasticity theory, where the fluid stress-strain relation is described with a Kelvin–Voigt mechanical model, thus incorporating viscoelasticity. First, we obtain the differential equations from first principles by defining the strain and kinetic energies and the dissipation function. We perform a plane-wave analysis to obtain the wave velocity and attenuation. The validity of the theory is demonstrated with three examples, namely, considering a porous solid saturated with a nN pore fluid, a nN fluid containing solid inclusions and a pure nN fluid. The analysis shows that the fluid may cause a negative velocity dispersion of the fast P (S)-wave velocities, that is velocity decreases with frequency. In acoustics, velocity increases with frequency (anomalous dispersion in optics). Furthermore, the fluid viscoelasticity has not a relevant effect on the wave responses observed in conventional field and laboratory tests. A comparison with previous theories supports the validity of the theory, which is useful to analyse wave propagation in a porous medium saturated with a fluid of high viscosity. [ABSTRACT FROM AUTHOR]

Seismic wave propagation is affected by wave‐induced local fluid flow between stiff pores and multiscale fractures. To investigate this phenomenon, forced vibration (1–100 Hz) and ultrasonic (106 Hz) measurements are performed on two tight sandstones with complex pore geometry in dry and water‐saturated scenarios. Porosity, permeability, and ultrasonic velocities were also measured at different differential pressures. The results indicate that the nonlinear behavior of these properties is strongly influenced by the presence of cracks, and the correlations between the permeability/ultrasonic velocities and porosity are different. A wave propagation model is then developed in which penny‐shaped cracks are inclusions introduced stepwise into a porous medium to describe the wave anelasticity in a wide frequency band at different pressures. As a result, the model provides good agreement with the measured P‐wave velocity dispersion, and the pore aspect ratio spectrum and crack radii are determined. We then compare the estimated crack radii and pore size distributions from nuclear magnetic resonance spectroscopy. Published data of a tight sandstone and a low porosity sandstone in the frequency range (2–200, 106) and (1–3,000) Hz at different differential pressures are also analyzed to validate the model. The aspect ratios, volume fractions, and radii of pores/cracks are used to describe the measured permeability. The present work can provide new insights into the geophysical properties of reservoir rocks with complex pore geometry. Plain Language Summary: Multiscale cracks exert an important influence on the dispersion and attenuation of broadband acoustic waves. To better understand the underlying mechanism, we measured frequency‐dependent elastic moduli and pressure‐dependent porosity, permeability, and ultrasonic wave velocities of two tight sandstones. These data show that velocity dispersion is observable in water‐saturated rock and that permeability and ultrasonic velocities are differentially correlated with porosity. A proposed wave propagation model, in which penny‐shaped cracks are inclusions introduced stepwise into a porous medium, successfully describes the experimental data and other published data on tight sandstone and low‐porosity sandstone, and aspect ratios, volume fractions, and radii of pores/cracks are determined. Key Points: We propose a model to investigate the effect of multiscale cracks on seismic waves in tight sandstonesThe pore aspect ratio and crack‐radius spectra are estimated at different differential pressuresThe permeability is obtained according to the pore geometry [ABSTRACT FROM AUTHOR]

The wave-induced local fluid flow mechanism is relevant to the complex heterogeneity of pore structures in rocks. The analysis of the local fluid flow mechanism is useful for accurately describing the wave propagation characteristics in reservoir rocks. In the exploration and production of hydrocarbon reservoirs, the real stratum may be partially saturated with a multi-phase fluid mixture in general. Therefore, it is of great significance to investigate the wave velocity dispersion and attenuation features in relation to pore structures and fluids. In this work, the characteristics of fabric microstructures are obtained on the basis of pressure dependency of dry rock moduli using the effective medium theory. A novel anelasticity theoretical model for the wave propagation in a partially-saturated medium is presented by combining the extended Gurevich squirt-flow model and White patchy-saturation theory. Numerical simulations are used to analyze wave propagation characteristics that depend on water saturation, external patchy diameter, and viscosity. We consider a tight sandstone from the Qingyang area of the Ordos Basin in west China and perform ultrasonic measurements under partial saturation states and different confining pressures, where the basic properties of the rock are obtained at the full gas saturation. The comparison of experimental data and theoretical modeling results shows a fairly good agreement, indicating that the new theory is effective. [ABSTRACT FROM AUTHOR]

The exploration and exploitation of shale oil is an important aspect in the oil industry. Seismic properties and well‐log data are essential to establish wave‐propagation models. Specifically, the description of wave dispersion and attenuation under complex geological conditions needs proper lithological and petrophysical information. This complex physical mechanism has to be considered if a traditional modeling approach is adopted. In this sense, machine learning (ML) techniques provide new possibilities for this purpose. We compare two deep‐neural‐network (DNN)‐based wave propagation models. In the first (pure data‐driven), a DNN is trained to connect seismic attributes, such as wave velocities, to multivariate functions of rock‐physics properties. By training DNNs with different initial parameters, the uncertainty of the proposed method can be quantified. The second method assumes the form of the wave equations. Then, the elastic constants of the constitutive relations are predicted by DNNs. The resulting dynamical equations describe the dispersion and attenuation and wavefield simulations can be performed to obtain more information. On the basis of a test, the two kinds of wave‐propagation models yield acceptable estimations of the seismic properties, with the second approach showing a broader application because the DNN is trained without S wave data. The methodologies illustrate that the new wave‐propagation model based on ML has high precision and can be general in terms of rheological description. Plain Language Summary: We develop a model to estimate the underground conditions, based on two machine‐learning techniques by using measurement data as input. The first uses pure data‐driven surrogate models, while the second is based on the classical wave propagation model in fluid‐saturated porous medium proposed by Biot, where deep neural networks predict the elastic coefficients so that the theoretical seismic properties match the actual measurements. The two models are tested on the basis of data from five wells located in a shale‐oil reservoir area. Compared with the pure data‐driven model, the wave‐equation method may be used to predict more reservoir properties. Key Points: Machine‐learning methods are proposed to drive wave propagation based on seismic dataOptimal prediction results of seismic attributes are obtained by feeding rock‐physics properties into the trained modelsThe proposed method has the potential to obtain additional physical quantities and shows a broad application prospect [ABSTRACT FROM AUTHOR]

SUMMARY: We consider surface-wave propagations at an interface separating a fluid layer and a double-porosity medium embedded with cracks. The theory is based on a generalization of the Biot-Rayleigh model from spherical cavities to penny-shaped cracks randomly embedded into a host medium, where mesoscopic local fluid flow (LFF) plays an important role. We derive closed-form dispersion equations of surface waves, based on potentials and suitable boundary conditions (BCs), to obtain the phase velocity and attenuation by using numerical iterations. Two special cases are considered by letting the thickness of the fluid (water) layer to be zero and infinity. We obtain pseudo-Rayleigh and pseudo-Stoneley waves for zero and infinite thickness and high-order surface modes for finite nonzero thickness. Numerical examples confirm that the LFF affects the propagation at low frequencies, causing strong attenuation, whereas the impact of BCs is mainly observed at high frequencies, due to the propagation of slow wave modes. The crack density mainly affects the level of attenuation, whereas the aspect ratio the location of the relaxation peak. The fundamental mode undergoes a significant velocity dispersion, whose location moves to low frequencies as the thickness increases. In all cases, there also exist two slower surface modes that resemble the two slow body waves, only present for sealed BCs. [ABSTRACT FROM AUTHOR]

Mesoscopic P-wave attenuation in layered, partially saturated thermo-poroelastic media is analyzed by combining the theories of Biot poroelasticity and Lord–Shulman thermoelasticity (BLS). The attenuation is quantified by estimating the quality factor Q. The mesoscopic attenuation effect, commonly referred to as wave-induced fluid flow (WIFF), is the process that converts fast compressional and shear waves into slow diffusive Biot waves at mesoscopic heterogeneities larger than the pore scale, but much smaller than the dominant wavelengths. This effect was first modeled in White's isothermal theory by quantifying the seismic response of a periodic sequence of planar porous layers that are alternately saturated with gas or water. This work presents a numerical extension of White's theory for the non-isothermal case in this type of sequence. For this purpose, an initial-boundary-value problem (IBVP) for the BLS wave propagation equations is solved using the finite element method, where the particle velocity field is recorded at uniformly distributed receivers. The quality factor is estimated using spectral-ratio and frequency-shift methods. The Q-estimates show that thermal effects influence the attenuation of the P-wave and the velocity dispersion compared to the isothermal case. [ABSTRACT FROM AUTHOR]

We simulate variations in the P- and S-wave velocities due to CO2 injection in basalt, caused by water replacement and mineralization in the form of carbonate cements. Carbonization can be monitored with seismic data on the basis of changes in velocity. We consider basalts composed of (Fe, Ca and Mg) minerals mainly, and the pore space is assumed to be initially saturated with water. The mineralization is based on a model that estimates the volume of precipitated cement and resulting porosity loss, and the velocities as a function of injection time are calculated with a patchy cement model. In the example, there are two stages: 1 week of CO2 injection and 29 weeks of frame cementation. After 1 week, water has been replaced by CO2, where density and P-wave velocity decrease, and cementation implies an increase in bulk density, P-wave velocity, whereas porosity and surface area are reduced. To our knowledge there are no models relating the effects of carbon mineralization on seismic properties. The present model can be applied to other rocks, where CO2 injection induces mineralization and storage. [ABSTRACT FROM AUTHOR]

We study the wave anelasticity (attenuation and velocity dispersion) of a periodic set of three flat porous layers saturated by two immiscible fluids. The fluids are very dissimilar in properties, namely gas, oil, and water, and, at most, three layers are required to study the problem from a general point of view. The sequence behaves as viscoelastic and transversely isotropic (VTI) at wavelengths much longer than the spatial period. Wave propagation causes fluid flow and slow P modes, inducing anelasticity. The fluids are characterized by capillary forces and relative permeabilities, which allow for the existence of two slow modes and the presence of dissipation, respectively. The methodology to study the physics is based on a finite-element uspcaling approach to compute the complex and frequency-dependent stiffnesses of the effective VTI medium. The results of the experiments indicate that there is higher dissipation and anisotropy compared to the widely used model based on an effective fluid that ignores the effects of surface tension (capillarity) and viscous flow interference between the two fluid phases. [ABSTRACT FROM AUTHOR]