Estimation of Dry-Rock Elastic Moduli Based on the Simulation of Mud-Filtrate Invasion Effects on Borehole Acoustic Logs

2009 ◽  
Vol 12 (06) ◽  
pp. 898-911 ◽  
Author(s):  
Tobiloluwa B. Odumosu ◽  
Carlos Torres-Verdín ◽  
Jesús M. Salazar ◽  
Jun Ma ◽  
Benjamin Voss ◽  
...  

Summary Reliable estimates of dry-rock elastic properties are critical to the accurate interpretation of the seismic response of hydrocarbon reservoirs. We describe a new method for estimating elastic moduli of rocks in-situ based on the simulation of mud-filtrate invasion effects on resistivity and acoustic logs. Simulations of mud-filtrate invasion account for the dynamic process of fluid displacement and mixing between mud-filtrate and hydrocarbons. The calculated spatial distributions of electrical resistivity are matched against resistivity logs by adjusting the underlying petrophysical properties. We then perform Biot-Gassmann fluid substitution on the 2D spatial distributions of fluid saturation with initial estimates of dry-bulk (kdry) modulus and shear rigidity (µdry) and a constraint of Poisson's ratio (?d) typical of the formation. This process generates 2D spatial distributions of compressional and shear-wave velocities and density. Subsequently, sonic waveforms are simulated to calculate shear-wave slowness. Initial estimates of the dry-bulk modulus are progressively adjusted using a modified Gregory-Pickett (1963) solution of Biot's (1956) equation to estimate a shear rigidity that converges to the well-log value of shear-wave slowness. The constraint on dynamic Poisson's ratio is then removed and a refined estimate of the dry-bulk modulus is obtained by both simulating the acoustic log (monopole) and matching the log-derived compressional-wave slowness. This technique leads to reliable estimates of dry-bulk moduli and shear rigidity that compare well to laboratory core measurements. Resulting dry-rock elastic properties can be used to calculate seismic compressional-wave and shear-wave velocities devoid of mud-filtrate invasion effects for further seismic-driven reservoir-characterization studies.

Geophysics ◽  
1976 ◽  
Vol 41 (5) ◽  
pp. 895-921 ◽  
Author(s):  
A. R. Gregory

The influence of saturation by water, oil, gas, and mixtures of these fluids on the densities, velocities, reflection coefficients, and elastic moduli of consolidated sedimentary rocks was determined in the laboratory by ultrasonic wave propagation methods. Twenty rock samples varying in age from Pliocene to early Devonian and in porosity from 4 to 41 percent were tested under uniform pressures equivalent to subsurface depths of 0 to 18,690 ft. Fluid saturation effects on compressional‐wave velocity are much larger in low‐porosity than in high‐porosity rocks; this correlation is strengthened by elevated pressures but is absent at atmospheric pressure. At a frequency of 1 MHz, the shear‐wave velocities do not always decrease when liquid pore saturants are added to rocks as theorized by Biot; agreement with theory is dependent upon pressure, porosity, fluid‐mineral chemical interactions, and the presence of microcracks in the cementing material. Experimental results support the belief that lower compressional‐wave velocities and higher reflection coefficients are obtained in sedimentary rocks that contain gas. Replacing pore liquids with gas markedly reduces the elastic moduli of rocks, and the effect is enhanced by decreasing pressure. As a rule, the moduli decrease as the porosity increases; Poisson’s ratio is an exception to the rule. Liquid and gas saturation in consolidated rocks can also be distinguished by the ratio of compressional and shear wave velocities [Formula: see text]. The potential diagnostic value of elastic moduli in seismic exploration may stimulate interest in the use of shear‐wave reflection methods in the field.


Geophysics ◽  
1985 ◽  
Vol 50 (4) ◽  
pp. 571-581 ◽  
Author(s):  
J. P. Castagna ◽  
M. L. Batzle ◽  
R. L. Eastwood

New velocity data in addition to literature data derived from sonic log, seismic, and laboratory measurements are analyzed for clastic silicate rocks. These data demonstrate simple systematic relationships between compressional and shear wave velocities. For water‐saturated clastic silicate rocks, shear wave velocity is approximately linearly related to compressional wave velocity and the compressional‐to‐shear velocity ratio decreases with increasing compressional velocity. Laboratory data for dry sandstones indicate a nearly constant compressional‐to‐shear velocity ratio with rigidity approximately equal to bulk modulus. Ideal models for regular packings of spheres and cracked solids exhibit behavior similar to the observed water‐saturated and dry trends. For dry rigidity equal to dry bulk modulus, Gassmann’s equations predict velocities in close agreement with data from the water‐saturated rock.


Geophysics ◽  
1995 ◽  
Vol 60 (6) ◽  
pp. 1627-1633 ◽  
Author(s):  
Bart W. Tichelaar ◽  
Klaas W. van Luik

Borehole sonic waveforms are commonly acquired to produce logs of subsurface compressional and shear wave velocities. To this purpose, modern borehole sonic tools are usually equipped with various types of acoustic sources, i.e., monopole and dipole sources. While the dipole source has been specifically developed for measuring shear wave velocities, we found that the dipole source has an advantage over the monopole source when determining compressional wave velocities in a very slow formation consisting of unconsolidated sands with a porosity of about 35% and a shear wave velocity of about 465 m/s. In this formation, the recorded compressional refracted waves suffer from interference with another wavefield component identified as a leaky P‐wave, which hampers the determination of compressional wave velocities in the sands. For the dipole source, separation of the compressional refracted wave from the recorded waveforms is accomplished through bandpass filtering since the wavefield components appear as two distinctly separate contributions to the frequency spectrum: a compressional refracted wave centered at a frequency of 6.5 kHz and a leaky P‐wave centered at 1.3 kHz. For the monopole source, the frequency spectra of the various waveform components have considerable overlap. It is therefore not obvious what passband to choose to separate the compressional refracted wave from the monopole waveforms. The compressional wave velocity obtained for the sands from the dipole compressional refracted wave is about 2150 m/s. Phase velocities obtained for the dispersive leaky P‐wave excited by the dipole source range from 1800 m/s at 1.0 kHz to 1630 m/s at 1.6 kHz. It appears that the dipole source has an advantage over the monopole source for the data recorded in this very slow formation when separating the compressional refracted wave from the recorded waveforms to determine formation compressional wave velocities.


Sensors ◽  
2020 ◽  
Vol 20 (21) ◽  
pp. 6291
Author(s):  
Jong-Sub Lee ◽  
Geunwoo Park ◽  
Yong-Hoon Byun ◽  
Changho Lee

A modified oedometer cell for measuring the applied stresses and elastic waves at the top and bottom of the specimen is developed to evaluate the effect of the side friction on the stress dependence of the elastic wave velocities. In the modified cell, two load cells are installed at the top and bottom plates, respectively. To generate and detect the compressional and shear waves, a pair of piezo disk elements and a pair of bender elements are mounted at both the top and bottom plates. Experimental results show that the stresses measured at the bottom are smaller than those measured at the top during the loading and vice versa during unloading, regardless of the densities and heights of the specimens. Under nearly saturated conditions, the compressional wave velocities remain almost constant for the entire stress state. With plotting stresses measured at top, the shear wave velocities during unloading are greater than those during loading, whereas with plotting stresses measured at bottom, the shear wave velocities during unloading are smaller than those during loading owing to the side friction. The vertical effective stress may be simply determined from the average values of the stresses measured at the top and bottom of the specimens.


Geophysics ◽  
2007 ◽  
Vol 72 (1) ◽  
pp. E1-E13 ◽  
Author(s):  
Michael A. Zimmer ◽  
Manika Prasad ◽  
Gary Mavko ◽  
Amos Nur

Knowledge of the pressure dependences of seismic velocities in unconsolidated sands is necessary for the remote prediction of effective pressures and for the projection of velocities to unsampled locations within shallow sand layers. We have measured the compressional- and shear-wave velocities and bulk, shear, and P-wave moduli at pressures from [Formula: see text] in a series of unconsolidated granular samples including dry and water-saturated natural sands and dry synthetic sand and glass-bead samples. The shear-wave velocities in these samples demonstrate an average pressure dependence approximately proportional to the fourth root of the effective pressure [Formula: see text], as commonly observed at lower pressures. For the compressional-wave velocities, theexponent in the pressure dependence of individual dry samples is consistently less than the exponent for the shear-wave velocity of the same sample, averaging 0.23 for the dry sands and 0.20 for the glass-bead samples. These pressure dependences are generally consistent over the entire pressure range measured. A comparison of the empirical results to theoretical predictions based on Hertz-Mindlin effective-medium models demonstrates that the theoretical models vastly overpredict the shear moduli of the dry granular frame unless the contacts are assumed to have no tangential stiffness. The models also predict a lower pressure exponent for the moduli and velocities [Formula: see text] than is generally observed in the data. We attribute this discrepancy in part to the inability of the models to account for decreases in the amount of slip or grain rotation occurring at grain-to-grain contacts with increasing pressure.


Geophysics ◽  
2001 ◽  
Vol 66 (3) ◽  
pp. 763-771 ◽  
Author(s):  
Myung W. Lee ◽  
Timothy S. Collett

Downhole‐measured compressional- and shear‐wave velocities acquired in the Mallik 2L-38 gas hydrate research well, northwestern Canada, reveal that the dominant effect of gas hydrate on the elastic properties of gas hydrate‐bearing sediments is as a pore‐filling constituent. As opposed to high elastic velocities predicted from a cementation theory, whereby a small amount of gas hydrate in the pore space significantly increases the elastic velocities, the velocity increase from gas hydrate saturation in the sediment pore space is small. Both the effective medium theory and a weighted equation predict a slight increase of velocities from gas hydrate concentration, similar to the field‐observed velocities; however, the weighted equation more accurately describes the compressional- and shear‐wave velocities of gas hydrate‐bearing sediments. A decrease of Poisson’s ratio with an increase in the gas hydrate concentration is similar to a decrease of Poisson’s ratio with a decrease in the sediment porosity. Poisson’s ratios greater than 0.33 for gas hydrate‐bearing sediments imply the unconsolidated nature of gas hydrate‐bearing sediments at this well site. The seismic characteristics of gas hydrate‐bearing sediments at this site can be used to compare and evaluate other gas hydrate‐bearing sediments in the Arctic.


1979 ◽  
Vol 24 (90) ◽  
pp. 313-319
Author(s):  
Joseph F. Kirchner ◽  
Charles R. Bentley

AbstractSeismic short-refraction studies were carried out at five stations on the Ross Ice Shelf during the 1976–77 summer season as part of the comprehensive Ross Ice Shelf Geophysical and Glaciological Survey. Measurements of the velocities of compressional waves were made at each location. Compressional wave velocities were measured along more than one azimuth at three sites, and shear wave velocities (both components) at two. Travel-time curves were fitted to an exponential expression by means of a non-linear least-squares regression technique. The errors in the apparent velocities are estimated to be about ±50 m s–1 at short distances, diminishing to about ±10 m s–1 near the ends of the profiles. Compressional-wave velocities show only slight variations with azimuth and only over certain depth intervals, showing that constant-velocity surfaces are essentially horizontal. Shear-wave velocities, however, exhibit large variations according to azimuth and polarization, indicating that transverse isotropy is violated at least in the upper 30–40 m of the ice shelf. It is believed that the anisotropy is caused by structural details in the firn perhaps modified by preferred crystal orientation and that it may arise at least partly from anisotropic stresses in the ice shelf.


1974 ◽  
Vol 64 (2) ◽  
pp. 375-385 ◽  
Author(s):  
Richard E. Warrick

abstract Measurements were made of the shear- and compressional-wave velocities necessary to determine the seismic response of a site on the edge of the San Francisco Bay, California. The shear wave is the most difficult to isolate because of the complications of earlier arriving energy. Downhole and surface measurements yielded shear-wave velocities of 90 to 130 m/sec and compressional-wave velocities of 1,400 m/sec for soft mud. The underlying older sediments had interval shear-wave velocities of 270 to 380 m/sec, increasing with depth. The corresponding P velocities for the older sediments were 1,600 to 1,840 m/sec. The horizontal traction technique of shear-wave generation proved to be the most successful of the methods tested.


Geophysics ◽  
1991 ◽  
Vol 56 (12) ◽  
pp. 2129-2138 ◽  
Author(s):  
M. A. Payne

In an effort to understand better the amplitude variation with offset for reflections from an oil sand and the sensitivity of the AVO response to shear‐wave velocity variations, I studied synthetic and field gathers collected from an onshore field in the Gulf of Mexico basin. A wave‐equation‐based modeling program generated the synthetic seismic gathers using both measured and estimated shear‐wave velocities. The measured shear‐wave velocities came from a quadrupole sonic tool. The estimated shear‐wave velocities were obtained by applying published empirical and theoretical equations which relate shear‐wave velocities to measured compressional‐wave velocities. I carefully processed the recorded seismic data with a controlled‐amplitude processing stream. Comparison of the synthetic gathers with the processed field data leads to the conclusion that the model containing the measured shear‐wave velocities matches the field data much better than the model containing the estimated shear‐wave velocities. Therefore, existing equations which relate shear‐wave velocities to compressional‐wave velocities yield estimates which are not sufficiently accurate for making quantitative comparisons of synthetic and field gathers. Even small errors in the shear‐wave velocities can have a large impact on the output. Such errors can lead to an incomplete and perhaps inaccurate understanding of the amplitude‐versus‐offset response. This situation can be remedied by collecting shear‐wave data for use in amplitude‐versus‐offset modeling, and for building databases to generate better shear‐wave velocity estimator equations.


Sign in / Sign up

Export Citation Format

Share Document