Effects of fluids and dual-pore systems on pressure-dependent velocities and attenuations in carbonates

Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. N35-N47 ◽  
Author(s):  
Remy Agersborg ◽  
Tor Arne Johansen ◽  
Morten Jakobsen ◽  
Jeremy Sothcott ◽  
Angus Best
Keyword(s):  
Fluid Flow ◽  
Confining Pressure ◽  
Substitution Effects ◽  
Fluid Substitution ◽  
Pore System ◽  
S Wave ◽  
Low Velocity ◽  
Wave Velocities ◽  
Three Samples ◽  
Local Fluid Flow

The effects of fluid substitution on P- and S-wave velocities in carbonates of complex texture are still not understood fully. The often-used Gassmann equation gives ambiguous results when compared with ultrasonic velocity data. We present theoretical modeling of velocity and attenuation measurements obtained at a frequency of [Formula: see text] for six carbonate samples composed of calcite and saturated with air, brine, and kerosene. Although porosities (2%–14%) and permeabilities [Formula: see text] are relatively low, velocity variations are large. Differences between the highest and lowest P- and S-wave velocities are about 18% and 27% for brine-saturated samples at 60 and [Formula: see text] effective pressure, respectively. S-wave velocities are measured for two orthogonal polarizations; for four of six samples, anisotropy is revealed. TheGassmann model underpredicts fluid-substitution effects by [Formula: see text] for three samples and by as much as 5% for the rest of the six samples. Moreover, when dried, they also show decreasing attenuation with increasing confining pressure. To model this behavior, we examine a pore model made of two pore systems: one constitutes the main and drainable porosity, and the other is made of undrained cracklike pores that can be associated with grain-to-grain contacts. In addition, these dried rock samples are modeled to contain a fluid-filled-pore system of grain-to-grain contacts, potentially causing local fluid flow and attenuation. For the theoretical model, we use an inclusion model based on the [Formula: see text]-matrix approach, which also considers effects of pore texture and geometry, and pore fluid, global- and local-fluid flow. By using a dual-pore system, we establish a realistic physical model consistently describing the measured data.

scholarly journals How Nanopores Influence Dry-Frame VP Pressure Sensitivity

2021 ◽  
Vol 9 ◽  
Author(s):  
Rohit Raj ◽  
Priyank Jaiswal ◽  
Yulun Wang ◽  
G. Michael Grammer ◽  
Ralf J. Weger
Keyword(s):  
Fluid Model ◽  
Confining Pressure ◽  
P Wave ◽  
Pore Shape ◽  
Pressure Sensitivity ◽  
S Wave ◽  
Shape Distribution ◽  
Transit Times ◽  
Three Samples ◽  
The Individual

This paper investigates how nanopore size distribution influences dry-frame P-wave velocity (VP) pressure sensitivity. The study uses a set of twenty-three samples belonging to a single vertical core from the Mississippian-age Meramec formation of the mid-continent US. Individual samples had their facies interpreted, composition estimated, He-gas porosity (ΦHe) determined, and P-wave and S-wave transit times systematically measured for dry core-plugs in a 5–40 MPa loading and unloading cycle. Data from the unloading cycle were linearized in the log scale, and the slope of the best fitting line was considered as a representative of the dry-frame VP pressure sensitivity. A series of photomicrographs from each sample were analyzed using image processing methods to obtain the shape and size of the individual pores, which were mostly in the nanopore (10−6–10–9 m) scale. At the outset, the pore-shape distribution plots were used to identify and discard samples with excessive cracks and complex pores. When the remaining samples were compared, it was found that within the same facies and pore-shape distribution subgroups VP pressure sensitivity increased as the dominant pore-size became smaller. This was largely independent of ΦHe and composition. The paper postulates that at the nanopore scale in the Meramec formation, pores are mostly isolated, and an increase in the confining pressure increased the bulk moduli of the fluids in the isolated pores, which in turn increased the VP pressure sensitivity. The study proposes incorporating this effect quantitatively through a dual-fluid model where the part of the fluid in unconnected pores is considered compressible while the remaining is considered incompressible. Results start to explain the universal observation of why the presence of microporosity quintessentially enhances VP pressure sensitivity.

Prediction of seismic-wave velocities in rock at various confining pressures based on unconfined data

Geophysics ◽  
2014 ◽  
Vol 79 (4) ◽  
pp. D235-D242 ◽  
Author(s):  
Ali Reza Najibi ◽  
Mohammad Reza Asef
Keyword(s):  
Rock Specimen ◽  
Confining Pressure ◽  
The United States ◽  
Wellbore Stability ◽  
Pressure Curve ◽  
S Wave ◽  
Confining Stress ◽  
Rock Density ◽  
Wave Velocities ◽  
Velocity Pressure

Laboratory measurement of P- and S-wave velocities ([Formula: see text] and [Formula: see text], respectively) under confining pressure indicates that with an increase in confining pressure, [Formula: see text] and [Formula: see text] will increase. The trend is exponential at low pressures, transitioning to linear above a critical pressure. However, the trend of the velocity-pressure curve for each rock specimen may be determined knowing the coefficients of this curve. We first studied how the coefficients of the velocity-pressure curve were expected to be functions of elastic moduli. Then, four empirical equations were used to estimate four coefficients of the velocity-pressure curve, using the rock density and [Formula: see text] and [Formula: see text] at atmospheric pressure (unconfined conditions). This analysis was carried out based on laboratory experiments on 285 rock specimens of different lithology from around the world, namely the United States, China, Germany, Iran, and deep-sea-drilling projects. For each rock specimen, [Formula: see text] and [Formula: see text] were measured at different confining stress levels, rendering more than 4000 data points. The accuracy of the estimated wave velocities was on the order of 2%–3% of the measured values on average. This methodology is especially valuable for prediction and analysis of the rock behavior at deep well conditions. This is applicable for predicting geophysical properties of the earth’s crust at depth, geomechanical study of hydrocarbon and geothermal reservoirs, wellbore stability analysis, and in situ stress determination.

scholarly journals The influence of pore system change during CO2 storage on 4D seismic interpretation

2019 ◽  
Vol 74 ◽  
pp. 81
Author(s):  
Lin Li ◽  
Jinfeng Ma
Keyword(s):  
Co2 Storage ◽  
Seismic Inversion ◽  
Pore System ◽  
S Wave ◽  
Wave Velocities ◽  
Rock Testing ◽  
Amplitude Versus Offset ◽  
Weyburn Field ◽  
4D Seismic ◽  
Amplitude Versus

A 4D seismic forward model constitutes the foundation of 4D seismic inversion. Here, in combination with the Gassmann equation, the Digby model is improved to calculate the S-wave velocity, and the resulting equation is verified using rock testing results. Then, considering the influences of changes in the pore pressure, CO2 saturation and porosity on the P- and S- wave velocities, rock testing results from a CO2 injection area in the Weyburn field are used to calculate the P- and S-wave velocities of the reservoir. These P- and S-wave velocities are found to overlap under different pressure conditions with or without considering porosity variations. Therefore, two-layer models and well models are developed to simulate synthetic seismograms; the models considering porosity variations may provide greater seismic responses and different Amplitude Versus Offset (AVO) trends in the synthetic seismogram profiles than those without considering porosity variations. Thus, porosity variations must be considered when establishing 4D seismic forward models.

scholarly journals An experimental study to investigate the physical and dynamic elastic properties of Eagle Ford shale rock samples

2021 ◽  
Author(s):  
Faisal Altawati ◽  
Hossein Emadi ◽  
Rayan Khalil
Keyword(s):  
Elastic Properties ◽  
Injection Pressure ◽  
Confining Pressure ◽  
Xrd Analysis ◽  
Low Permeability ◽  
S Wave ◽  
Eagle Ford ◽  
Core Samples ◽  
Wave Velocities ◽  
Dynamic Elastic Properties

AbstractUnconventional resources, such as Eagle Ford formation, are commonly classified for their ultra-low permeability, where pore sizes are in nano-scale and pore-conductivity is low, causing several challenges in evaluating unconventional-rock properties. Several experimental parameters (e.g., diffusion time of gas, gas injection pressure, method of permeability measurement, and confining pressure cycling) must be considered when evaluating the ultra-low permeability rock's physical and dynamic elastic properties measurements, where erroneous evaluations could be avoided. Characterizing ultra-low permeability samples' physical and elastic properties helps researchers obtain more reliable information leading to successful evaluations. In this study, 24 Eagle Ford core samples' physical and dynamic elastic properties were evaluated. Utilizing longer diffusion time and higher helium injection pressure, applying complex transient method, and cycling confining pressure were considered for porosity, permeability, and velocities measurements. Computerized tomography (CT) scan, porosity, permeability, and ultrasonic wave velocities were conducted on the core samples. Additionally, X-ray Diffraction (XRD) analysis was conducted to determine the mineralogical compositions. Porosity was measured at 2.07 MPa injection pressure for 24 h, and the permeability was measured using a complex transient method. P- and S-wave velocities were measured at two cycles of five confining pressures (up to 68.95 MPa). The XRD analysis results showed that the tested core samples had an average of 81.44% and 11.68% calcite and quartz, respectively, with a minor amount of clay minerals. The high content of calcite and quartz in shale yields higher velocities, higher Young's modulus, and lower Poisson's ratio, which enhances the brittleness that is an important parameter for well stimulation design (e.g., hydraulic fracturing). The results of porosity and permeability showed that porosity and permeability vary between 5.3–9.79% and 0.006–12 µD, respectively. The Permeability–porosity relation of samples shows a very weak correlation. P- and S-wave velocities results display a range of velocity up to 6206 m/s and 3285 m/s at 68.95 MPa confining pressure, respectively. Additionally, S-wave velocity is approximately 55% of P-wave velocity. A correlation between both velocities is established at each confining pressure, indicating a strong correlation. Results illustrated that applying two cycles of confining pressure impacts both velocities and dynamic elastic moduli. Ramping up the confining pressure increases both velocities owing to compaction of the samples and, in turn, increases dynamic Young's modulus and Poisson's ratio while decreasing bulk compressibility. Moreover, the results demonstrated that the above-mentioned parameters' values (after decreasing the confining pressure to 13.79 MPa) differ from the initial values due to the hysteresis loop, where the loop is slightly opened, indicating that the alteration is non-elastic. The findings of this study provide detailed information about the rock physical and dynamic elastic properties of one of the largest unconventional resources in the U.S.A, the Eagle Ford formation, where direct measurements may not be cost-effective or feasible.

Impact of change in pore-fill material on P-wave velocity

Geophysics ◽  
2014 ◽  
Vol 79 (6) ◽  
pp. D399-D407 ◽  
Author(s):  
Nishank Saxena ◽  
Gary Mavko
Keyword(s):  
Shear Modulus ◽  
Wave Velocity ◽  
P Wave ◽  
Exact Equation ◽  
Fluid Substitution ◽  
Seismic Velocities ◽  
S Wave ◽  
Wave Velocities ◽  
Fill Material ◽  
P Wave Velocity

The problem of predicting the change in seismic velocities (P-wave and S-wave) upon the change in pore-fill material properties is commonly known as substitution. For isotropic rocks, P- and S-wave velocities are fundamentally linked to the effective P-wave and shear moduli. The change in the S-wave velocity or shear modulus upon fluid substitution can be predicted with Gassmann’s equations starting only with the initial S-wave velocity. However, predicting changes in P-wave velocity or the P-wave modulus requires knowledge of the initial P- and S-wave velocities. We initiated a rigorous derivation of the P-wave modulus for fluid and solid substitution in monomineralic isotropic rocks for cases in which an estimate of the S-wave velocity or shear modulus is not available. For the general case of solid substitution, the exact equation for the P-wave modulus depends on parameters that are usually unknown. However, for fluid substitution, fewer parameters are required. As Poisson’s ratio increases for the mineral in the rock frame, the dependence of exact substitution on these unknown parameters decreases. As a result, in the absence of shear velocity, P-wave modulus fluid substitution can, for example, be performed with higher confidence for rocks with a calcite or dolomite frame than it can for rocks with quartz frame. We evaluated a recipe for applying the new P-wave modulus fluid substitution. This improves on existing work and is recommended for practice.

scholarly journals Subsonic near-surface P-velocity and low S-velocity observations using propagator inversion

Geophysics ◽  
2005 ◽  
Vol 70 (4) ◽  
pp. R15-R23 ◽  
Author(s):  
Robbert van Vossen ◽  
Andrew Curtis ◽  
Jeannot Trampert
Keyword(s):  
Shear Velocity ◽  
Confining Pressure ◽  
Compressional Wave ◽  
P Wave ◽  
Soil Conditions ◽  
Detailed Knowledge ◽  
S Wave ◽  
Near Surface ◽  
Unconsolidated Sands ◽  
Wave Velocities

Detailed knowledge of near-surface P- and S-wave velocities is important for processing and interpreting multicomponent land seismic data because (1) the entire wavefield passes through and is influenced by the near-surface soil conditions, (2) both source repeatability and receiver coupling also depend on these conditions, and (3) near-surface P- and S-wave velocities are required for wavefield decomposition and demultiple methods. However, it is often difficult to measure these velocities with conventional techniques because sensitivity to shallow-wave velocities is low and because of the presence of sharp velocity contrasts or gradients close to the earth's free surface. We demonstrate that these near-surface P- and S-wave velocities can be obtained using a propagator inversion. This approach requires data recorded by at least one multicomponent geophone at the surface and an additional multicomponent geophone at depth. The propagator between them then contains all information on the medium parameters governing wave propagation between the geophones at the surface and at depth. Hence, inverting the propagator gives local estimates for these parameters. This technique has been applied to data acquired in Zeist, the Netherlands. The near-surface sediments at this site are unconsolidated sands with a thin vegetation soil on top, and the sediments considered are located above the groundwater table. A buried geophone was positioned 1.05 m beneath receivers on the surface. Propagator inversion yielded low near-surface velocities, namely, 270 ± 15 m/s for the compressional-wave velocity, which is well below the sound velocity in air, and 150 ± 9 m/s for the shear velocity. Existing methods designed for imaging deeper structures cannot resolve these shallow material properties. Furthermore, velocities usually increase rapidly with depth close to the earth's surface because of increasing confining pressure. We suspect that for this reason, subsonic near-surface P-wave velocities are not commonly observed.

Bulk modulus dispersion and attenuation in sandstones

Geophysics ◽  
2015 ◽  
Vol 80 (2) ◽  
pp. D111-D127 ◽  
Author(s):  
L. Pimienta ◽  
J. Fortin ◽  
Y. Guéguen
Keyword(s):  
Fluid Flow ◽  
Bulk Modulus ◽  
Low Frequency ◽  
Confining Pressure ◽  
S Wave ◽  
Standard Samples ◽  
Stress Strain ◽  
Fontainebleau Sandstone ◽  
Attenuation Characteristic ◽  
Confining Pressures

We report experimental data on the frequency dependence of bulk elastic modulus in porous sandstones. A new methodology was developed to investigate the dispersion/attenuation phenomena on a rock’s bulk modulus [Formula: see text] for varying confining pressures in the range of 1–50 MPa and fluids of varying viscosities (i.e., air, glycerin, and water). This methodology combined (1) ultrasonic (i.e., [Formula: see text]) P- and S-wave velocity measurements, leading to the high-frequency (HF) [Formula: see text], (2) stress-strain measurements from forced periodic oscillations of confining pressure at low-frequency (LF) ranges (i.e., [Formula: see text]), leading to [Formula: see text] and [Formula: see text], and (3) pore-pressure measurement to document the induced fluid-flow in the LF range (i.e., [Formula: see text]). The stress-strain method was first checked using three standard samples: glass, gypsum, and Plexiglas samples. Over the frequency and pressure range of the apparatus [Formula: see text] was stable and accurate and the lowest measurable LF attenuation was [Formula: see text]. The methodology was applied to investigate Fontainebleau sandstone samples of 7% and 9% porosity. The [Formula: see text] and [Formula: see text] exhibited correlated variations, which also correlated with an experimental evidence of frequency-dependent fluid-flow out of the sample. Attenuation peaks as high as [Formula: see text] and [Formula: see text] are measured. The attenuation/dispersion measured under glycerin saturation was compared to Biot-Gassmann predictions. The overall behavior of one sample was consistent with a dispersion/attenuation characteristic of the drained/undrained transition. On the reverse, the other sample exhibited exotic behaviors as the measurements were underestimated by the drained/undrained transition and indicated a direct transition from drained to unrelaxed domain. These different behaviors were consistent with the values of the critical frequencies expected for the drained/undrained (i.e., [Formula: see text]) and relaxed/unrelaxed (i.e., [Formula: see text]) transitions.

3D-printed rock models: Elastic properties and the effects of penny-shaped inclusions with fluid substitution

Geophysics ◽  
2016 ◽  
Vol 81 (6) ◽  
pp. D669-D677 ◽  
Author(s):  
Long Huang ◽  
Robert R. Stewart ◽  
Nikolay Dyaur ◽  
Jose Baez-Franceschi
Keyword(s):  
3D Printing ◽  
Elastic Properties ◽  
Wave Velocity ◽  
Bedding Plane ◽  
P Wave ◽  
Fluid Substitution ◽  
S Wave ◽  
Wave Velocities ◽  
P Wave Velocity ◽  
3D Printed

3D printing techniques (additive manufacturing) using different materials and structures provide opportunities to understand porous or fractured materials and fluid effects on their elastic properties. We used a 3D printer (Stratasys Dimension SST 768) to print one “solid” cube model and another with penny-shaped inclusions. The 3D printing process builds materials, layer by layer, producing a slight “bedding” plane, somewhat similar to a sedimentary process. We used ultrasonic transducers (500 kHz) to measure the P- and S-wave velocities. The input printing material was thermoplastic with a density of [Formula: see text], P-wave velocity of [Formula: see text], and S-wave velocity of [Formula: see text]. The solid cube had a porosity of approximately 6% and a density of [Formula: see text]. Its P-wave velocity was [Formula: see text] in the bedding direction and [Formula: see text] normal to bedding. We observed S-wave splitting with fast and slow velocities of 879 and [Formula: see text], respectively. Quality factors for P- and S-waves were estimated using the spectral-ratio method with [Formula: see text] ranging from 15 to 17 and [Formula: see text] from 24 to 27. By introducing penny-shaped inclusions along the bedding direction in a 3D printed cube, we created a more porous volume with density of [Formula: see text] and porosity of 24%. The inclusions significantly decreased the P-wave velocity to 1706 and [Formula: see text] parallel and normal to the bedding plane. The fast and slow S-wave velocities also decreased to 812 and [Formula: see text]. A fluid substitution experiment, performed with water, increased (20%–46%) P-wave velocities and decreased (9%–10%) S-wave velocities. Theoretical predictions using Schoenberg’s linear-slip theory and Hudson’s penny-shaped theory were calculated, and we found that both theories matched the measurements closely (within 5%). The 3D printed material has interesting and definable properties and is an exciting new material for understanding wave propagation, rock properties, and fluid effects.

Dual attenuation peaks revealing mesoscopic and microscopic fluid flow in partially oil-saturated Fontainebleau sandstones

2020 ◽  
Vol 224 (3) ◽  
pp. 1670-1683
Author(s):  
Liming Zhao ◽  
Genyang Tang ◽  
Chao Sun ◽  
Jianguo Zhao ◽  
Shangxu Wang
Keyword(s):  
Fluid Flow ◽  
Reservoir Characterization ◽  
Seismic Analysis ◽  
Confining Pressure ◽  
Time Lapse ◽  
Seismic Attenuation ◽  
Fluid Distribution ◽  
Oil Viscosity ◽  
Measurement Results ◽  
And Shifts

SUMMARY We conducted stress–strain oscillation experiments on dry and partially oil-saturated Fontainebleau sandstone samples over the 1–2000 Hz band at different confining pressures to investigate the wave-induced fluid flow (WIFF) at mesoscopic and microscopic scales and their interaction. Three tested rock samples have similar porosity between 6 and 7 per cent and were partially saturated to different degrees with different oils. The measurement results exhibit a single or two attenuation peaks that are affected by the saturation degree, oil viscosity and confining pressure. One peak, exhibited by all samples, shifts to lower frequencies with increasing pressure, and is mainly attributed to grain contact- or microcrack-related squirt flow based on modelling of its characteristics and comparison with other experiment results for sandstones. The other peak is present at smaller frequencies and shifts to higher frequencies as the confining pressure increases, showing an opposite pressure dependence. This contrast is interpreted as the result of fluid flow patterns at different scales. We developed a dual-scale fluid flow model by incorporating the squirt flow effect into the patchy saturation model, which accounts for the interaction of WIFFs at microscopic and mesoscopic scales. This model provides a reasonable interpretation of the measurement results. Our broad-frequency-band measurements give physical evidence of WIFFs co-existing at two different scales, and combining with modelling results, it suggests that the WIFF mechanisms, related to pore microstructure and fluid distribution, interplay with each other and jointly control seismic attenuation and dispersion at reservoir conditions. These observations and modelling results are useful for quantitative seismic interpretation and reservoir characterization, specifically they have potential applications in time-lapse seismic analysis, fluid prediction and reservoir monitoring.

Sign in / Sign up

Export Citation Format

Share Document