Seismic imaging using microearthquakes induced by hydraulic fracturing

Geophysics ◽  
1994 ◽  
Vol 59 (1) ◽  
pp. 102-112 ◽  
Author(s):  
Lisa V. Block ◽  
C. H. Cheng ◽  
Michael C. Fehler ◽  
W. Scott Phillips
Keyword(s):  
Hydraulic Fracturing ◽  
Wave Velocity ◽  
Seismic Imaging ◽  
P Wave ◽  
Basement Rock ◽  
S Wave ◽  
Fractured Zone ◽  
Wave Velocities ◽  
The Absolute ◽  
Wave Arrival

Seismic imaging using microearthquakes induced by hydraulic fracturing produces a three-dimensional (3-D), S-wave velocity model of the fractured zone, improves the calculated locations of the microearthquakes, and may lead to better estimates of fractureplane orientations, fracture density, and water flow paths. Such information is important for predicting the amount of heat energy that may be extracted from geothermal reservoir. A fractured zone was created at the Los Alamos Hot Dry Rock Reservoir in north-central New Mexico within otherwise impermeable basement rock by injecting [Formula: see text] of water into a borehole under high pressure at a depth of 3.5 km. Induced microearthquakes were observed using four borehole seismometers. The P-wave and S-wave arrival times have been inverted to find the 3-D velocity structures and the microearthquake locations and origin times. The inversion was implemented using the separation of parameters technique, and constraints were incorporated to require smooth velocity structures and to restrict the velocities within the fractured region to be less than or equal to the velocities of the unfractured basement rock. The rms amval time residuals decrease by 11–15 percent during the joint hypocenter-velocity inversion. The average change in the microearthquake locations is 20–27 m, depending on the smoothing parameter used. Tests with synthetic data imply that the absolute locations may improve by as much as 35 percent, while the relative locations may improve by 40 percent. The general S-wave velocity patterns are reliable, but the absolute velocity values are not uniquely determined. However, studies of inversions using various degrees of smoothing suggest that the S-wave velocities decrease by at least 13 percent in the most intensely fractured regions of the reservoir. The P-wave velocities are poorly constrained because the P-wave traveltime perturbations caused by the fluid-filled fractures are small compared to the amval time noise level. The significant difference in the relative signal-to-noise levels of the P-wave and S-wave arrival time data, coupled with the limited ray coverage, can produce a bias in the computed [Formula: see text] ratios, and corresponding systematic rotation of the microearthquake cluster. These adverse effects were greatly reduced by applying a [Formula: see text] lower bound based on the [Formula: see text] ratio of the unfractured basement rock.

Wavefield simulation and data-acquisition-scheme analysis for LWD acoustic tools in very slow formations

Geophysics ◽  
2011 ◽  
Vol 76 (3) ◽  
pp. E59-E68 ◽  
Author(s):  
Hua Wang ◽  
Guo Tao
Keyword(s):  
Wave Velocity ◽  
Cutoff Frequency ◽  
P Wave ◽  
Flexural Wave ◽  
S Wave ◽  
Low Frequencies ◽  
Wave Velocities ◽  
Dipole System ◽  
P Wave Velocity ◽  
Wavefield Simulation

Propagating wavefields from monopole, dipole, and quadrupole acoustic logging-while-drilling (LWD) tools in very slow formations have been studied using the discrete wavenumber integration method. These studies examine the responses of monopole and dipole systems at different source frequencies in a very slow surrounding formation, and the responses of a quadrupole system operating at a low source frequency in a slow formation with different S-wave velocities. Analyses are conducted of coherence-velocity/slowness relationships (semblance spectra) in the time domain and of the dispersion characteristics of these waveform signals from acoustic LWD array receivers. These analyses demonstrate that, if the acoustic LWD tool is centralized properly and is operating at low frequencies (below 3 kHz), a monopole system can measure P-wave velocity by means of a “leaky” P-wave for very slow formations. Also, for very slow formations a dipole system can measure the P-wave velocity via a leaky P-wave and can measure the S-wave velocity from a formation flexural wave. With a quadrupole system, however, the lower frequency limit (cutoff frequency) of the drill-collar interference wave would decrease to 5 kHz and might no longer be neglected if the surrounding formation becomes a very slow formation, with S-wave velocities at approximately 500 m/s.

Acoustic measurements in unconsolidated sands at low effective pressure and overpressure detection

Geophysics ◽  
2002 ◽  
Vol 67 (2) ◽  
pp. 405-412 ◽  
Author(s):  
Manika Prasad
Keyword(s):  
Wave Velocity ◽  
P Wave ◽  
Acoustic Measurements ◽  
S Wave ◽  
Deepwater Drilling ◽  
Unconsolidated Sands ◽  
Wave Velocities ◽  
P Wave Velocity ◽  
Shallow Water Flows ◽  
Low Pressures

Shallow water flows and over‐pressured zones are a major hazard in deepwater drilling projects. Their detection prior to drilling would save millions of dollars in lost drilling costs. I have investigated the sensitivity of seismic methods for this purpose. Using P‐wave information alone can be ambiguous, because a drop in P‐wave velocity (Vp) can be caused both by overpressure and by presence of gas. The ratio of P‐wave velocity to S‐wave velocity (Vp/Vs), which increases with overpressure and decreases with gas saturation, can help differentiate between the two cases. Since P‐wave velocity in a suspension is slightly below that of the suspending fluid and Vs=0, Vp/Vs and Poisson's ratio must increase exponentially as a load‐bearing sediment approaches a state of suspension. On the other hand, presence of gas will also decrease Vp but Vs will remain unaffected and Vp/Vs will decrease. Analyses of ultrasonic P‐ and S‐wave velocities in sands show that the Vp/Vs ratio, especially at low effective pressures, decreases rapidly with pressure. At very low pressures, Vp/Vs values can be as large as 100 and higher. Above pressures greater than 2 MPa, it plateaus and does not change much with pressure. There is significant change in signal amplitudes and frequency of shear waves below 1 MPa. The current ultrasonic data shows that Vp/Vs values can be invaluable indicators of low differential pressures.

EXPERIMENTAL EVALUATION OF ULTRASONIC VELOCITIES AND ANISOTROPY IN THE TUSCALOOSA MARINE SHALE FORMATION

2020 ◽  
pp. 1-62 ◽  
Author(s):  
Jamal Ahmadov ◽  
Mehdi Mokhtari
Keyword(s):  
Wave Velocity ◽  
Mineralogical Composition ◽  
Organic Content ◽  
P Wave ◽  
S Wave ◽  
Velocity Anisotropy ◽  
Wave Velocities ◽  
Marine Shale ◽  
P Wave Velocity ◽  
Degree Of Anisotropy

Tuscaloosa Marine Shale (TMS) formation is a clay- and organic-rich emerging shale play with a considerable amount of hydrocarbon resources. Despite the substantial potential, there have been only a few wells drilled and produced in the formation over the recent years. The analyzed TMS samples contain an average of 50 wt% total clay, 27 wt% quartz and 14 wt% calcite and the mineralogy varies considerably over the small intervals. The high amount of clay leads to pronounced anisotropy and the frequent changes in mineralogy result in the heterogeneity of the formation. We studied the compressional (VP) and shear-wave (VS) velocities to evaluate the degree of anisotropy and heterogeneity, which impact hydraulic fracture growth, borehole instabilities, and subsurface imaging. The ultrasonic measurements of P- and S-wave velocities from five TMS wells are the best fit to the linear relationship with R2 = 0.84 in the least-squares criteria. We observed that TMS S-wave velocities are relatively lower when compared to the established velocity relationships. Most of the velocity data in bedding-normal direction lie outside constant VP/VS lines of 1.6–1.8, a region typical of most organic-rich shale plays. For all of the studied TMS samples, the S-wave velocity anisotropy exhibits higher values than P-wave velocity anisotropy. In the samples in which the composition is dominated by either calcite or quartz minerals, mineralogy controls the velocities and VP/VS ratios to a great extent. Additionally, the organic content and maturity account for the velocity behavior in the samples in which the mineralogical composition fails to do so. The results provide further insights into TMS Formation evaluation and contribute to a better understanding of the heterogeneity and anisotropy of the play.

Subsurface Shear Wave Velocity Characterization Using P-Wave Seismograms in Central and Eastern North America

2016 ◽  
Vol 32 (1) ◽  
pp. 143-169 ◽  
Author(s):  
Byungmin Kim ◽  
Youssef M. A. Hashash ◽  
Ellen M. Rathje ◽  
Jonathan P. Stewart ◽  
Sidao Ni ◽  
...  
Keyword(s):  
North America ◽  
Wave Velocity ◽  
P Wave ◽  
Motion Prediction ◽  
Eastern North America ◽  
Ground Motion Prediction ◽  
Profile Data ◽  
S Wave ◽  
Wave Velocities ◽  
Linear Relationships

The time-averaged shear ( S) wave velocity in the upper 30 meters of sediment ( V S30) is a widely used site parameter for ground motion prediction. When unavailable from measurements, as is often the case at accelerograph stations in Central and Eastern North America (CENA), V S30 is typically estimated from proxies. We propose an alternative for CENA based on a theoretical relationship between S-wave velocity and the ratio of radial to vertical components of the compressional ( P)-wave–dominated portion of the velocity time series. This method is applied to 31 CENA accelerograph sites having measured S-wave velocity profiles. Time-averaged S-wave velocities to depth z ( V SZ) from the proposed method agree well with those from measurements. We develop linear relationships between V SZ and V S30 using CENA S-wave velocity profile data. Values of V S30 established from the proposed method (including depth extrapolation) have lower dispersion relative to data ( σln V = 0.43) than do estimates from available CENA proxies.

Determining P- and S-wave velocities and Q-values from single ultrasound transmission measurements performed on cylindrical rock samples: it’s possible, when…

2020 ◽  
Author(s):  
Marc S. Boxberg ◽  
Mandy Duda ◽  
Katrin Löer ◽  
Wolfgang Friederich ◽  
Jörg Renner
Keyword(s):  
Wave Velocity ◽  
P Wave ◽  
S Wave ◽  
Standard Material ◽  
Rock Samples ◽  
Intrinsic Attenuation ◽  
Finite Element Software ◽  
Wave Velocities ◽  
P Wave Velocity

<p>Determining elastic wave velocities and intrinsic attenuation of cylindrical rock samples by transmission of ultrasound signals appears to be a simple experimental task, which is performed routinely in a range of geoscientific and engineering applications requiring characterization of rocks in field and laboratory. P- and S-wave velocities are generally determined from first arrivals of signals excited by specifically designed transducers. A couple of methods exist for determining the intrinsic attenuation, most of them relying either on a comparison between the sample under investigation and a standard material or on investigating the same material for various geometries.</p><p>Of the three properties of interest, P-wave velocity is certainly the least challenging one to determine, but dispersion phenomena lead to complications with the consistent identification of frequency-dependent first breaks. The determination of S-wave velocities is even more hampered by converted waves interfering with the S-wave arrival. Attenuation estimates are generally subject to higher uncertainties than velocity measurements due to the high sensitivity of amplitudes to experimental procedures. The achievable accuracy of determining S-wave velocity and intrinsic attenuation using standard procedures thus appears to be limited.</p><p>We pursue the determination of velocity and attenuation of rock samples based on full waveform modeling and inversion. Assuming the rock sample to be homogeneous - an assumption also underlying standard analyses - we quantify P-wave velocity, S-wave velocity and intrinsic P- and S-wave attenuation from matching a single ultrasound trace with a synthetic one numerically modelled using the spectral finite-element software packages SPECFEM2D and SPECFEM3D. We find that enough information on both velocities is contained in the recognizable reflected and converted phases even when nominal P-wave sensors are used. Attenuation characteristics are also inherently contained in the relative amplitudes of these phases due to their different travel paths. We present recommendations for and results from laboratory measurements on cylindrical samples of aluminum and rocks with different geometries that we also compare with various standard analysis methods. The effort put into processing for our approach is particularly justified when accurate values and/or small variations, for example in response to changing P-T-conditions, are of interest or when the amount of sample material is limited.</p>

Estimating shear‐wave velocities from P-wave and converted‐wave data

Geophysics ◽  
1999 ◽  
Vol 64 (2) ◽  
pp. 504-507 ◽  
Author(s):  
Franklyn K. Levin
Keyword(s):  
Wave Velocity ◽  
Seismic Data ◽  
P Wave ◽  
S Wave ◽  
Converted Wave ◽  
Wave Data ◽  
Wave Velocities ◽  
Shear Wave Velocities ◽  
P Wave Velocity ◽  
Growing Body

Tessmer and Behle (1988) show that S-wave velocity can be estimated from surface seismic data if both normal P-wave data and converted‐wave data (P-SV) are available. The relation of Tessmer and Behle is [Formula: see text] (1) where [Formula: see text] is the S-wave velocity, [Formula: see text] is the P-wave velocity, and [Formula: see text] is the converted‐wave velocity. The growing body of converted‐wave data suggest a brief examination of the validity of equation (1) for velocities that vary with depth.

Seismological studies in a rock body at Chalk River, Ontario and their relation to fractures

1982 ◽  
Vol 19 (8) ◽  
pp. 1535-1547 ◽  
Author(s):  
C. Wright
Keyword(s):  
Wave Velocity ◽  
Test Site ◽  
P Wave ◽  
Seismic Velocities ◽  
S Wave ◽  
P Waves ◽  
Rock Body ◽  
Near Surface ◽  
Wave Velocities ◽  
P Wave Velocity

Seismological experiments have been undertaken at a test site near Chalk River, Ontario that consists of crystalline rocks covered by glacial sediments. Near-surface P and S wave velocity and amplitude variations have been measured along profiles less than 2 km in length. The P and S wave velocities were generally in the range 4.5–5.6 and 2.9–3.2 km/s, respectively. These results are consistent with propagation through fractured gneiss and monzonite, which form the bulk of the rock body. The P wave velocity falls below 5.0 km/s in a region where there is a major fault and in an area of high electrical conductivity; such velocity minima are therefore associated with fracture systems. For some paths, the P and 5 wave velocities were in the ranges 6.2–6.6 and 3.7–4.1 km/s, respectively, showing the presence of thin sheets of gabbro. Temporal changes in P travel times of up to 1.4% over a 12 h period were observed where the sediment cover was thickest. The cause may be changes in the water table. The absence of polarized SH arrivals from specially designed shear wave sources indicates the inhomogeneity of the test site. A Q value of 243 ± 53 for P waves was derived over one relatively homogeneous profile of about 600 m length. P wave velocity minima measured between depths of 25 and 250 m in a borehole correlate well with the distribution of fractures inferred from optical examination of borehole cores, laboratory measurements of seismic velocities, and tube wave studies.

scholarly journals Elastic reflection-based waveform inversion with a nonlinear approach

Geophysics ◽  
2017 ◽  
Vol 82 (6) ◽  
pp. R309-R321 ◽  
Author(s):  
Qiang Guo ◽  
Tariq Alkhalifah
Keyword(s):  
Wave Velocity ◽  
Equivalent Stress ◽  
Waveform Inversion ◽  
P Wave ◽  
S Wave ◽  
Velocity Models ◽  
Wave Velocities ◽  
Full Waveform ◽  
Highly Nonlinear ◽  
Elastic Reflection

Full-waveform inversion (FWI) is a highly nonlinear problem due to the complex reflectivity of the earth, and this nonlinearity only increases under the more expensive elastic assumption. In elastic media, we need a good initial P-wave velocity and even better initial S-wave velocity models with accurate representation of the low model wavenumbers for FWI to converge. However, inverting for the low-wavenumber components of P- and S-wave velocities using reflection waveform inversion (RWI) with an objective to fit the reflection shape, rather than produce reflections, may mitigate the limitations of FWI. Because FWI, performing as a migration operator, is preferred of the high-wavenumber updates along reflectors. We have developed an elastic RWI that inverts for the low-wavenumber and perturbation components of the P- and S-wave velocities. To generate the full elastic reflection wavefields, we derive an equivalent stress source made up by the inverted model perturbations and incident wavefields. We update the perturbation and propagation parts of the velocity models in a nested fashion. Applications on the synthetic isotropic models and field data indicate that our method can efficiently update the low- and high-wavenumber parts of the models.

A strategy for nonlinear elastic inversion of seismic reflection data

Geophysics ◽  
1986 ◽  
Vol 51 (10) ◽  
pp. 1893-1903 ◽  
Author(s):  
Albert Tarantola
Keyword(s):  
Wave Velocity ◽  
Seismic Reflection ◽  
Wave Impedance ◽  
P Wave ◽  
S Wave ◽  
Seismic Reflection Data ◽  
S Waves ◽  
Wave Velocities ◽  
Reflection Data ◽  
P Wave Velocity

The problem of interpretation of seismic reflection data can be posed with sufficient generality using the concepts of inverse theory. In its roughest formulation, the inverse problem consists of obtaining the Earth model for which the predicted data best fit the observed data. If an adequate forward model is used, this best model will give the best images of the Earth’s interior. Three parameters are needed for describing a perfectly elastic, isotropic, Earth: the density ρ(x) and the Lamé parameters λ(x) and μ(x), or the density ρ(x) and the P-wave and S-wave velocities α(x) and β(x). The choice of parameters is not neutral, in the sense that although theoretically equivalent, if they are not adequately chosen the numerical algorithms in the inversion can be inefficient. In the long (spatial) wavelengths of the model, adequate parameters are the P-wave and S-wave velocities, while in the short (spatial) wavelengths, P-wave impedance, S-wave impedance, and density are adequate. The problem of inversion of waveforms is highly nonlinear for the long wavelengths of the velocities, while it is reasonably linear for the short wavelengths of the impedances and density. Furthermore, this parameterization defines a highly hierarchical problem: the long wavelengths of the P-wave velocity and short wavelengths of the P-wave impedance are much more important parameters than their counterparts for S-waves (in terms of interpreting observed amplitudes), and the latter are much more important than the density. This suggests solving the general inverse problem (which must involve all the parameters) by first optimizing for the P-wave velocity and impedance, then optimizing for the S-wave velocity and impedance, and finally optimizing for density. The first part of the problem of obtaining the long wavelengths of the P-wave velocity and the short wavelengths of the P-wave impedance is similar to the problem solved by present industrial practice (for accurate data interpretation through velocity analysis and “prestack migration”). In fact, the method proposed here produces (as a byproduct) a generalization to the elastic case of the equations of “prestack acoustic migration.” Once an adequate model of the long wavelengths of the P-wave velocity and of the short wavelengths of the P-wave impedance has been obtained, the data residuals should essentially contain information on S-waves (essentially P-S and S-P converted waves). Once the corresponding model of S-wave velocity (long wavelengths) and S-wave impedance (short wavelengths) has been obtained, and if the remaining residuals still contain information, an optimization for density should be performed (the short wavelengths of impedances do not give independent information on density and velocity independently). Because the problem is nonlinear, the whole process should be iterated to convergence; however, the information from each parameter should be independent enough for an interesting first solution.

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.

Sign in / Sign up

Export Citation Format

Share Document