Squirt Flow Research Articles

Effects of fluid saturation and viscosity on seismic dispersion characteristics in Berea sandstone

Despite additional availability of laboratory data from water-saturated sandstone at seismic frequencies, measurements of rock samples saturated with high viscous fluids, particularly at partial saturation, are still rare. To quantify the effects of fluid viscosity and saturation levels on seismic dispersion and attenuation characteristics, we conducted two comparative forced-oscillation measurements in partially saturated sandstone with varying fluid viscosity (e.g., water, glycerin) at seismic frequencies (2–400 Hz). The results demonstrate that fluid viscosity and saturation levels substantially influence the dispersion and attenuation characteristics at the measured frequencies. Significant dispersion and attenuation are observed in the presence of a relatively small amount of gas (approximately 6%–8%) for glycerin and water saturation cases but vary in their magnitudes and characteristic frequencies. Specifically, the maximum extensional attenuation (approximately 0.024) occurs at approximately 200 Hz for water-saturated rock at 94% saturation, whereas at approximately 30 Hz with a peak of 0.032 for glycerin-saturated rock at 92% saturation. Based on theoretical modeling analysis, we suggest that mesoscopic fluid flow might be a dominant mechanism accounting for the observed attenuation in partial water or glycerin saturation, while the microscopic (squirt) flow mechanism possibly dominates the fully saturated cases.

Experimental investigation of pore-filling substitution effect on frequency-dependent elastic moduli of Berea sandstone

SUMMARY Based on both forced oscillation and ultrasonic pulse transmission methods, we investigated solid pore infill influences on rock elastic moduli in a broad frequency range $[ {1 - 3000,\,\,{{10}^6}} ]$ Hz for different differential pressures. For a Berea sandstone sample, filled sequentially by solid (${22\,\,^{\rm{o}}}{\rm{C}}$), quasi-solid (${26\,\,^{\rm{o}}}{\rm{C}}$) and liquid (${34\,\,^{\rm{o}}}{\rm{C}}$) octadecane, a frequency-dependence was found for the Poisson's ratio, Young's modulus and bulk modulus, nevertheless, these elastic parameters were strongly suppressed by increasing pressures. Experimental measurements showed that shear wave velocity and modulus of solid-octadecane-filled samples are significantly larger than those of the dry and liquid-octadecane-filled ones, implying the potential stiffening effects related to solid infill in compliant pores. A three porosity structure model, which describes the solid stiffening effects related to equant, compliant and the intermediate pores with aspect ratios larger than those of compliant pores but much less than those of stiff pores, was used to compare against the experimentally measured elastic properties for octadecane pore infill, together with several other fluid/solid substitution theories. The agreement between experimental measurements and theoretical predictions is reasonably good for the sandstone tested, providing that the three porosity model can be applied for pressure- and frequency-dependent elastic moduli estimations for a viscoelastic pore-infill-saturated sandstone. Evaluating the combined squirt flow mechanism responsible for the observed moduli dispersion and attenuation is of great importance to reduce potential errors in seismic AVO inversion and 4-D seismic monitoring of gas-hydrate or bitumen-saturated reservoir, especially for reservoir rocks with complex microstructures and heterogeneous pore types.

Shear-wave attenuation anisotropy: a new constraint on mantle melt near the Main Ethiopian Rift

The behaviour of fluids in preferentially aligned fractures plays an important role in a range of dynamic processes within the Earth. In the near-surface, understanding systems of fluid-filled fractures is crucial for applications such as geothermal energy production, monitoring CO2 storage sites, and exploration for metalliferous sub-volcanic brines. Mantle melting is a key geodynamic process, exerting control over its composition and dynamic processes. Upper mantle melting weakens the lithosphere, facilitating rifting and other surface expressions of tectonic processes.Aligned fluid-filled fractures are an efficient mechanism for seismic velocity anisotropy, requiring very low volume fractions, but such rock physics models also predict significant shear-wave attenuation anisotropy. In comparison, the attenuation anisotropy expected for crystal preferred orietation mechanisms is negligible or would only operate outside of the seismic frequency band.Here we demonstrate a new method for measuring shear-wave attenuation anisotropy, apply it to synthetic examples, and make the first measurements of SKS attenuation anisotropy using data recorded at the station FURI, in Ethiopia. At FURI we measure attenuation anisotropy where the fast shear-wave has been more attenuated than the slow shear-wave. This can be explained by the presence of aligned fluids, most probably melts, in the upper mantle using a poroelastic squirt flow model. Modelling of this result suggests that a 1% melt fraction, hosted in aligned fractures dipping ca. 40° that strike perpendicular to the Main Ethiopian Rift, is required to explain the observed attenuation anisotropy. This agrees with previous SKS shear-wave splitting analysis which suggested a 1% melt fraction beneath FURI. The interpreted fracture strike and dip, however, disagrees with previous work in the region which interprets sub-vertical melt inclusions aligned parallel to the Main Ethiopian Rift which only produce attenuation anisotropy where the slow shear-wave is more attenuated. These results show that attenuation anisotropy could be a useful tool for detecting mantle melt, and may offer strong constraints on the extent and orientation of melt inclusions which cannot be achieved from seismic velocity anisotropy alone.

Simulating squirt flow in realistic rock models using graphical processing units (GPUs)

SUMMARY Understanding the underlying mechanisms of seismic attenuation and moduli dispersion in fluid-saturated cracked porous rocks is of great importance for the development of non-invasive methods to characterize the subsurface. Wave-induced fluid flow at the pore scale, so-called squirt flow, is responsible for seismic attenuation and moduli dispersion at sonic and ultra-sonic frequencies and may be relevant at seismic frequencies. The squirt flow associated attenuation is usually quantified using analytical models. However, numerical experiments suggest that the squirt flow related dissipation is sensitive to fine details of the pore geometry, which can only be modelled numerically. Most of the existing numerical studies explore this phenomenon using simplified models, and there is a lack of numerical studies that model the physics in realistic pore geometries with sufficient numerical resolution. As a result, the impact of wave-induced fluid flow on the effective static and time-dependent mechanical characteristics in realistic settings is still poorly understood. I address these issues by developing a numerical method to model the effective mechanical properties of a hydromechanically coupled system at the pore scale suitable for graphical processing units. A numerical evaluation of attenuation and modulus dispersion due to squirt flow in models based on 3-D microtomography images of cracked Carrara marble is presented. It is shown that the local hydraulic conductivity can be quantitatively estimated from the numerically evaluated effective properties. The accuracy of the numerical results is carefully analysed. This study improves the understanding of the underlying mechanisms of attenuation and moduli dispersion in fluid-saturated cracked rocks. The new method can be applied to model squirt flow for entire laboratory samples in the centimetre scale which was not possible a decade ago.

Frequency-dependent elastic properties of fracture-induced VTI rocks in a fluid-saturated porous and microcracked background

Fractures are widely distributed underground. The stiffness matrix of fractured rocks has been extensively investigated in a fluid-saturated porous background medium. However, the existing stiffness models only incorporate the attenuation mechanism of wave-induced fluid flow (WIFF). For macroscopic fractures, the elastic scattering (ES) of fractures cannot be ignored. To alleviate this issue, a frequency-dependent stiffness matrix model is developed, including the mesoscopic WIFF between fractures and background (FB-WIFF), the microscopic squirt flow, and the macroscopic ES from the fractures. By combining the far-field scattered wavefields of normal incident P and SV waves with the linear slip theory, the dynamic full-stiffness matrices for fracture-induced effective vertically transversely isotropic rocks in a fluid-saturated porous and microcracked background are derived. Then, the P, SV, and SH-wave velocities and attenuation can be obtained through the Kelvin-Christoffel equation. The results indicate that the FB-WIFF mechanism significantly affects the velocities and attenuation of the P and SV waves but has almost no effect on the SH wave, whereas the squirt flow and ES mechanisms affect the velocities and attenuation of all the P, SV, and SH waves. For validation, the model is compared with existing models and previous experimental ultrasonic data.

The effects of pore structure on wave dispersion and attenuation due to squirt flow: A dynamic stress-strain simulation on a simple digital pore-crack model

Squirt flow, a phenomenon typically observed in porous cracked rocks, occurs due to the contrasting compressibility between the pores and cracks, leading to the pore pressure diffusion and dissipation of wave energy. Understanding the influence of pore structure on wave dispersion and attenuation signatures due to squirt flow is essential for interpreting seismic and sonic logging data in various fields of earth and energy sciences, such as hydrocarbon exploration, geothermal energy exploitation, and CO2 sequestration. In this study, we construct a simple digital pore-crack model with varying pore structures and use a dynamic stress-strain simulation approach to characterize wave dispersion and attenuation signatures due to squirt flow. Numerical simulation suggests that, in addition to the commonly considered parameters such as porosity, crack density, and crack aspect ratio, additional pore structure parameters, such as pore size, pore aspect ratio, crack orientation, crack length, and crack width, significantly affect the wave dispersion and attenuation signatures. Increasing the pore size leads to the pronounced enhancement of attenuation magnitude and a decrease in characteristic frequency. We demonstrate that variations in crack length have a pronounced impact on attenuation magnitude, whereas crack width is decisive in controlling the characteristic frequency. Furthermore, it is found that the saturation paths (the gas filling the pore first or gas filling the crack first) considerably affect the velocity-saturation and attenuation-saturation relationship, suggesting that the coupling effects of pore structure and fluid distribution complicate the fluid pressure diffusion and wave attenuation behaviors. The presented results offer insights for deciphering pore structure parameters using attenuation- or dissipation-related seismic attributes.

Laboratory experiments and theoretical study of pressure and fluid influences on acoustic response in tight rocks with pore microstructure

AbstractWave‐induced fluid flow is considered to be a major source of seismic attenuation and dispersion in porous rocks. From the physical description of partially saturated reservoirs, numerous analytical solutions based on upscaling homogenization theories have been employed to calculate equivalent frequency‐dependent poroelastic media. Nevertheless, dispersion and attenuation predictions are often not reasonably consistent with laboratory and field measurements in a broad frequency range, particularly due to influences of biphasic fluids and their distribution, presence of heterogeneities on various length scales, and pore microstructure. We investigate the role of pore microstructure on pressure and fluid saturation dependence of elastic velocities in tight sandstones. Previous work points out that differentiating the impacts of heterogeneities at various scales on dispersion within seismic exploration and sonic frequencies can be very difficult. In practice, this is because fluid‐related dispersion mechanisms are impossible to be independent. Thus, it is important for a theoretical and more quantitative analysis of the relative contribution of interrelated energy dissipation processes through a better understanding of combined influences due to the presence of microscopic and mesoscopic heterogeneities. Based on microscopic squirt flow and mesoscopic flow in a partially saturated medium, we develop a poroelastic model that allows evaluating the overall frequency‐dependent dispersion via considering a random distribution of fluid heterogeneities as well as the broadly distributed aspect ratio of compliant pores. Experimental validation of the model is accomplished via a comprehensive comparison of predictions with measurements of partially saturated velocities versus pressure and fluid for sandstones with specific pore microstructures.

Experimental Study of Seismic Wave Attenuation in Carbonate Rocks

Summary Seismic wave attenuation has a great potential for studying saturated and fractured media, due to its high sensitivity to the physical properties of geological media. However, accurately estimating this parameter can be challenging due to its sensitivity to signal noise, particularly in heterogeneous media such as carbonate rocks. This explains the paucity of attenuation studies carried out in carbonate rocks compared with sandstones, and the ambiguity around its mechanisms and its relationship with petrophysical properties. To investigate further, we conducted an experimental study of ultrasonic waveform signals (0.5–3 MHz) reordered under dry and fully saturation conditions in 13 samples covering a wide range of petrophysical values and subjected them to differential pressure reaching reservoir pressure. The resulting increase in attenuation magnitudes and their variation with pressure due to brine saturation were more pronounced than in velocity magnitudes, confirming the higher sensitivity of attenuation to fluid content. However, understanding the relationship between attenuation and petrophysical properties required a careful examination of the results and more elucidation about attenuation mechanisms. We suggested that multiple attenuation mechanisms coexist, including scattering, cracks slipping, solid frictional relative motion, and global and squirt flow. This explains the frequency dependence of attenuation, with higher magnitudes at sonic frequencies, where the squirt flow mechanism may be dominant. In contrast to sandstone, the magnitude of compressional to shear attenuation ratio (QP−1/QS−1) was found to be greater than unity in both dry and brine fully saturated carbonate samples at ultrasonic frequencies. This result may be due to the complex porosity structure of carbonate rocks, which makes it not appropriate to the sandstone rock physics models.

Stress and frequency dependence of wave velocities in saturated rocks based on acoustoelasticity with squirt-flow dissipation

SUMMARY We perform seismic and ultrasonic measurements in carbonate and shaley sandstone samples as a function of differential pressure. The velocities show a strong frequency and pressure dependence. The dispersion disappears with increasing pressure and the squirt flow in turn inhibits the pressure dependence. To model these effects, we combine the Gurevich's squirt-flow model with the Mori–Tanaka scheme and the David Zimmerman model, and extend it with third-order elastic constants, to obtain a frequency-dependent acoustoelasticity model. Comparisons between measurements from this study and literature and modelling results show that the P-wave velocity increases non-linearly first and then nearly linearly, dominated by crack closure and acoustoelasticity, respectively. The pressure dependence of wave velocities is reduced by liquid substitution and further by the squirt-flow mechanism. The effects of fluid properties and crack closure on P-wave velocity decrease with differential pressure. The results will feed a new model and help better understanding the wave propagation in pre-stressed rocks at different scales.

PP-Wave Reflection Coefficient Equation for HTI Media Incorporating Squirt Flow Effect and Frequency-Dependent Azimuthal AVO Inversion for Anisotropic Fluid Indicator

PP-Wave Reflection Coefficient Equation for HTI Media Incorporating Squirt Flow Effect and Frequency-Dependent Azimuthal AVO Inversion for Anisotropic Fluid Indicator

A simple and accurate model for attenuation and dispersion caused by squirt flow in isotropic porous rocks

Seismic waves propagating in fluid-saturated porous rocks exhibit attenuation and velocity dispersion in a broad range of frequencies. At sonic and ultrasonic frequencies, the attenuation is predominantly caused by fluid flow in cracks and grain contacts, so-called squirt flow. This physical mechanism for attenuation also may be relevant at seismic frequencies. We develop a simple and accurate analytical model for attenuation and dispersion caused by squirt flow in isotropic porous rocks. The input material properties for a specific rock model can be directly measured in a laboratory or calculated using analytical and numerical approaches. The results from our squirt flow model are compared with inherently accurate 3D numerical solutions for the same pore geometries. The analytical and numerical results are in good agreement. Furthermore, we observe that our analytical model is more accurate than the currently available analytical solution for squirt flow in isotropic porous rocks. MATLAB routines to reproduce the presented results are made available.

P- and SV-wave dispersion and attenuation in saturated microcracked porous rock with aligned penny-shaped fractures

P- and SV-wave dispersion and attenuation have been extensively investigated in saturated poroelastic media with aligned fractures. However, there are few existing models that incorporate the multiple wave attenuation mechanisms from the microscopic scale to the macroscopic scale. Hence, in this work, we developed a unified model to incorporate the wave attenuation mechanisms at different scales, which includes the microscopic squirt flow between the microcracks and pores, the mesoscopic wave-induced fluid flow between fractures and background (FB-WIFF), and the macroscopic Biot's global flow and elastic scattering (ES) from the fractures. Using Tang's modified Biot's theory and the mixed-boundary conditions, we derived the exact frequency-dependent solutions of the scattering problem for a single penny-shaped fracture with oblique incident P- and SV- waves. We then developed theoretical models for a set of aligned fractures and randomly oriented fractures using the Foldy approximation. The results indicated that microcrack squirt flow considerably influences the dispersion and attenuation of P- and SV-wave velocities. The coupling effects of microcrack squirt flow with the FB-WIFF and ES of fractures cause much higher velocity dispersion and attenuation for P waves than for SV waves. Randomly oriented fractures substantially reduce the attenuation caused by the FB-WIFF and ES, particularly for the ES attenuation of SV waves. Through a comparison with existing models in the limiting cases and previous experimental measurements, we validated our model.

Laboratory acousto-mechanical study into moisture-induced changes of elastic properties in intact granite

The water adsorption into pore spaces in brittle rocks affects wave velocity and transmitted amplitude of elastic waves. Experimental and theoretical studies have been performed to characterize moisture-induced elastodynamic variations due to macroporous effects; however, little attention has been paid to the manner in which wetting of nanopores affects elastic wave transmission. In this work, we extend our understanding of moisture-induced elastic changes in a microcracked nanopore-dominated medium where 80 % of the surface area exhibits pore diameters (also include microcrack widths) below 10 nm. We study acousto-mechanical response resulting from a gradual wetting on a free-standing intact Herrnholz granite specimen over 98 h using time-lapse ultrasonic and digital imaging techniques. Linkages between ultrasonic attributes and adsorption-induced stress/strain are established during the approach of the wetting front. We find that Gassmann theory, previously validated in channel-like nanoporous media, does not work properly in predicting the P-wave velocity increase of microcracked nanopore-dominated media at ultrasonic frequency. However, squirt flow – a theory recognized to characterize wave velocity increase and attenuation in microcracked macropore-dominated media at the pore scale – also accounts for the observed increase of P-wave velocity in microcracked nanopore-dominated media. The transmitted amplitude changes in direct P waves are explained and predicted by the elastic wave propagation within P-wave first Fresnel zone and reflection/refraction on the wetting front.

Experimental Validation and Calibration of the Galvin Model with Artificial Tight Sandstones with Controlled Fractures

The study of fractures in the subsurface is very important in unconventional reservoirs since they are the main conduits for hydrocarbon flow. For this reason, a variety of equivalent medium theories have been proposed for the estimation of fracture and fluid properties within reservoir rocks. Recently, the Galvin model has been put forward to model the frequency-dependent elastic moduli in fractured porous rocks and has been widely used to research seismic wave propagation in fractured rocks. We experimentally investigated the feasibility of applying the Galvin model in fractured tight stones. For this proposal, three artificial fractured tight sandstone samples with the same background porosity (11.7% ± 1.2%) but different fracture densities of 0.00, 0.0312, and 0.0624 were manufactured. The fracture thickness was 0.06 mm and the fracture diameter was 3 mm in all the fractured samples. Ultrasonic P- and S-wave velocities were measured at 0.5 MHz in a laboratory setting in dry and water-saturated conditions in directions at 0°, 45°, and 90° to the fracture normal. The results were compared with theoretical predictions of the Galvin model. The comparison showed that model predictions significantly underestimated P- and S- wave velocities as well as P-wave anisotropy in water-saturated conditions, but overestimated P-wave anisotropy in dry conditions. By analyzing the differences between the measured results and theoretical predictions, we modified the Galvin model by adding the squirt flow mechanism to it and used the Thomsen model to obtain the elastic moduli in high- and low-frequency limits. The modified model predictions showed good fits with the measured results. To the best of our knowledge, this is the first study to validate and calibrate the frequency-dependent equivalent medium theories in tight fractured rocks experimentally.

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

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, we consider 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 our 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 our method and model. We link 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.

Macroscopic mechanical properties of fluid-saturated sandstone at variable temperatures

Rocks can be viewed as composites of solid minerals with pores or cracks filled with softer material such as pore fluids, kerogen, bitumen, and other organic matter. The mechanical properties of highly viscous soft phases are highly sensitive to ambient temperatures and lead to temperature-dependent static and dynamic observations with the composite rock. However, the constituents operate by forming some effective (averaged) mechanical properties of the composite rock, and yet these averaged properties are still little known. To reveal such macroscopic temperature-dependent mechanical properties and measure their values in rock samples, a double-porosity model of porous rock with nonlinear viscosity is developed. The model is based on rigorous continuum mechanics with physically meaningful, real-valued, and time- and frequency-independent material properties and elegantly unifies the existing frequency-dependent microscopic squirt flow and mesoscopic wave-induced fluid flow models. The approach is used to accurately model the broad attenuation peaks and Young’s modulus dispersion observed in previously published laboratory experiments with glycerol-saturated Berea sandstone and invert for its mechanical properties. The observations are explained as mainly due to the temperature-dependent elastic coupling caused by non-Newtonian fluid within microcracks. Several hitherto unknown mechanical properties of the rock are constrained quantitatively: the average porosity of the microcracks, the effective high-pressure bulk modulus of the drained frame, the internal stiffness defect within the rock frame, the solid viscosities associated with bulk and shear deformations, and an exponent of nonlinearity for viscosity. These parameters constitute a Biot-consistent mechanical model of the rock, which can be used to simulate its behavior in arbitrary experimental environments. The rigorous first-principle model can be used in many applications: detailed and physically accurate interpretations of laboratory experiments, numerical wavefield simulations and seismic data inversion, reservoir characterization, geothermal exploration, thermal-enhanced oil recovery, and exploration for deep oil and gas resources in high-temperature environments.

Introducing realistic distributions of gas bubble size and inclusion aspect ratio to model the coexistence of gas hydrate and free gas

The coexisting free gas within the hydrate stability zone is recognized by laboratory and field studies. The conventional resistivity method fails to obtain hydrate saturation in the sediment with gas hydrate and free gas because both of them are electrical insulators. Hydrate saturation can be obtained from a velocity log based on a certain rock physics model. However, few studies of rock physics have been conducted on gas and hydrate coexistence. In this study, a more realistic distribution of gas bubble size and inclusion aspect ratio was introduced into an effective medium model to elucidate how such coexistence affects velocity and attenuation at a broad-band frequency. The modeling results obtained indicate that such coexistence decreases P-wave velocity at low frequencies, while increases P- and S-wave velocities at high frequencies. The velocity decrease of P-wave at low frequencies is attributed to the overwhelming effect of gas bubble vibration, while it seems not to affect P-wave velocity at high frequencies significantly. Moreover, coexisting gas hydrate and free gas may strengthen P-wave attenuation at the broad-band frequency range. It is found that the damping of oscillating gas bubbles dominates P-wave attenuation at low frequencies, while squirt flow in the porous hydrate is more likely to contribute to P-wave attenuation at high frequencies. Hydrate morphology is also a significant factor affecting P-wave velocity and attenuation at high frequencies. This research reveals that rock physics modeling can provide insight into how to characterize the coexisting gas hydrate and free gas as referring to the acoustic velocity and attenuation at seismic, sonic, and ultrasonic frequencies.

Measurements of Wave‐Induced Attenuation in Saturated Metapelite and the Band‐Limitation of Low‐Frequency Earthquakes

AbstractThe most common explanation for the depletion of high frequency waves that defines low‐frequency earthquakes (LFEs) and very low‐frequency earthquakes (VLFEs) is that fault rupture and slip are slower than typical earthquakes. However, it is difficult to rule out the possibility that the high frequency waves are produced during slip, but attenuated near the LFE source. One reason this hypothesis has been poorly tested is that there are no measurements of attenuation on the relevant rocks. We present the results of forced oscillation experiments that measure the frequency‐dependent attenuation of a chlorite‐rich metapelitic schist, a lithology found along the subduction plate boundary where LFEs and VLFEs have been documented. Experiments were run on dry and water‐saturated schist at effective pressures of 2–10 MPa and at frequencies of 2 × 10−5–30 Hz. We find that pore fluids and low effective pressure result in the attenuation of high frequencies. The frequency‐dependent attenuation is consistent with the concomitant operation of two wave‐induced fluid flow mechanisms, squirt flow, and patchy saturation. When the effects of these mechanisms are extrapolated to geologic conditions using rock physics models, our results predict that attenuation is capable of completely diminishing the frequencies depleted in LFEs and VLFEs. Therefore, LFEs and VLFEs may not necessarily record slow fault slip, but possibly the presence of high fluid pressure.

Effects of pore fluids on the hardening of shear modulus in tight dolomite: experimental observations and theoretical modelling

SUMMERY Understanding the influence of pore fluids on acoustic characteristics of rocks containing pores of low aspect ratio, that is, cracks is crucial for both laboratory and field observations of seismic exploration. Unlike the petrophysical aspects of conventional rock species, those of tight dolomites are relatively more complex and poorly understood. To address this problem, we investigated the ultrasonic velocity of dry and saturated tight dolomite in the effective pressure range of 2–60 MPa and analysed the variation in shear modulus upon fluid saturation. Our results showed that while undergoing the change in form from the dry state to being saturated with water, the rock may weaken or harden, and the shear modulus may correspondingly either decrease or increase. Shear modulus hardening was more apparent in dolomites, with microcracks mostly dispersed between and within particles of a high proportion of soft pores with a wide distribution range of aspect ratio. This study used the Walsh–Hudson model to satisfactorily simulate the hardening of shear modulus in fluid-saturated tight dolomites. The suppression of squirt flow that is responsible for the increased shear modulus of fluid-saturated tight dolomites tested at sufficiently high frequencies. As the change in shear modulus affects the prediction of longitudinal and shear wave velocities, using Gassmann's theory to fluid substitution in unconventional reservoirs such as tight dolomite may yield inaccurate results. Therefore, it may be necessary to consider variation in shear modulus while selecting the appropriate simulation model.

Stress Relaxing Simulation on Digital Rock: Characterize Attenuation Due To Wave‐Induced Fluid Flow and Scattering

AbstractDigital rock physics (DRP) is essential to understand the wave elastic and attenuation characteristics of complex porous rocks, offering insights into interpreting seismic and sonic logging data. Most of the existing DRP technology can only simulate wave‐induced fluid flow (WIFF) or scattering effects separately in a relatively narrow or single‐frequency band. To incorporate both the effects of WIFF and elastic scattering simultaneously, we develop a stress relaxing simulation (SRS) technique on digital rocks to characterize the P wave dispersion and attenuation signatures in a broad frequency band. We propose to use the combined methodology of curve resampling and wavelet‐packet denoising to suppress the high‐frequency noise of the stress relaxation curve. For a digital rock containing pore‐crack elements, we demonstrate that SRS is capable of characterizing the attenuation due to squirt flow and scattering as well as their interactions. The simulation results of wave dispersion and attenuation using SRS are overall consistent with other well‐developed DRP technology of transmitted wave simulation (high‐frequency band), forced oscillating simulation (intermediate‐frequency band), and static simulation (low‐frequency band), suggesting the reliability and robustness of the developed methodology. We also apply the SRS on a realistic 2D digital rock of porous sandstone. It is found that the attenuation due to elastic scattering of the microscopic heterogeneities is significantly higher than that caused by squirt flow, which might be caused by the lack of pore connectivity impeding the fluid exchange at the microscopic scale due to the limitation of 2D digital porous rocks.