Tomographic imaging of velocity and Q, with application to crosswell seismic data from the Gypsy Pilot Site, Oklahoma

Geophysics ◽  
1997 ◽  
Vol 62 (6) ◽  
pp. 1804-1811 ◽  
Author(s):  
Qingbo Liao ◽  
George A. McMechan
Keyword(s):  
Velocity Ratio ◽  
Synthetic Data ◽  
P Wave ◽  
Reconstruction Technique ◽  
Shift Method ◽  
New Information ◽  
P Wave Velocity ◽  
Iterative Reconstruction Technique ◽  
Interwell Connectivity

The centroid frequency shift method is implemented, tested with synthetic data, and applied to field data from three contiguous crosswell seismic experiments at the Gypsy Pilot in northern Oklahoma. The similtaneous iterative reconstruction technique is used for tomographic estimations of both P‐wave velocity and Q. No amplitude corrections or spreading loss corrections are needed for the Q estimation. The estimated in‐situ velocity and Q distributions correlate well with log data and local lithology. The Q/velocity ratio appears to correlate with the sand/shale ratio (ranging from an average of ∼15 s/km for the sand‐dominated lithologies to an average of ∼8.5 s/km for the shale‐dominated ones), with the result that new information is provided on interwell connectivity.

scholarly journals Simulation of the effect of stress-induced anisotropy on borehole compressional wave propagation

Geophysics ◽  
2014 ◽  
Vol 79 (4) ◽  
pp. D205-D216 ◽  
Author(s):  
Xinding Fang ◽  
Michael C. Fehler ◽  
Arthur Cheng
Keyword(s):  
Wave Propagation ◽  
Stress Concentration ◽  
Elastic Properties ◽  
Wave Velocity ◽  
Wave Amplitude ◽  
P Wave ◽  
Azimuthal Variation ◽  
P Wave Velocity ◽  
Velocity Region

Formation elastic properties near a borehole may be altered from their original state due to the stress concentration around the borehole. This can lead to an incorrect estimation of formation elastic properties measured from sonic logs. Previous work has focused on estimating the elastic properties of the formation surrounding a borehole under anisotropic stress loading. We studied the effect of borehole stress concentration on sonic logging in a moderately consolidated Berea sandstone using a two-step approach. First, we used an iterative approach, which combines a rock-physics model and a finite-element method, to calculate the stress-dependent elastic properties of the rock around a borehole subjected to an anisotropic stress loading. Second, we used the anisotropic elastic model obtained from the first step and a finite-difference method to simulate the acoustic response of the borehole. Although we neglected the effects of rock failure and stress-induced crack opening, our modeling results provided important insights into the characteristics of borehole P-wave propagation when anisotropic in situ stresses are present. Our simulation results were consistent with the published laboratory measurements, which indicate that azimuthal variation of the P-wave velocity around a borehole subjected to uniaxial loading is not a simple cosine function. However, on field scale, the azimuthal variation in P-wave velocity might not be apparent at conventional logging frequencies. We found that the low-velocity region along the wellbore acts as an acoustic focusing zone that substantially enhances the P-wave amplitude, whereas the high-velocity region caused by the stress concentration near the borehole results in a significantly reduced P-wave amplitude. This results in strong azimuthal variation of P-wave amplitude, which may be used to infer the in situ stress state.

The diminishing residual matrices method for surface-consistent statics — A review

Geophysics ◽  
2017 ◽  
Vol 82 (4) ◽  
pp. V257-V274
Author(s):  
Necati Gülünay
Keyword(s):  
Synthetic Data ◽  
Reconstruction Technique ◽  
Data Set ◽  
Arrival Times ◽  
Long Wavelength ◽  
Iterative Reconstruction Technique ◽  
Land Survey ◽  
First Arrival ◽  
Waveform Distortions ◽  
The Individual

The diminishing residual matrices (DRM) method can be used to surface-consistently decompose individual trace statics into source and receiver components. The statics to be decomposed may either be first-arrival times after the application of linear moveout associated with a consistent refractor as used in refraction statics or residual statics obtained by crosscorrelating individual traces with corresponding model traces (known as pilot traces) at the same common-midpoint (CMP) location. The DRM method is an iterative process like the well-known Gauss-Seidel (GS) method, but it uses only source and receiver terms. The DRM method differs from the GS method in that half of the average common shot and receiver terms are subtracted simultaneously from the observations at each iteration. DRM makes the under-constrained statics problem a constrained one by implicitly adding a new constraint, the equality of the contribution of shots and receivers to the solution. The average of the shot statics and the average of the receiver statics are equal in the DRM solution. The solution has the smallest difference between shot and receiver statics profiles when the number of shots and the number of receivers in the data are equal. In this case, it is also the smallest norm solution. The DRM method can be derived from the well-known simultaneous iterative reconstruction technique. Simple numerical tests as well as results obtained with a synthetic data set containing only the field statics verify that the DRM solution is the same as the linear inverse theory solution. Both algorithms can solve for the long-wavelength component of the statics if the individual picks contain them. Yet DRM method is much faster. Application of the method to the normal moveout-corrected CMP gathers on a 3D land survey for residual statics calculation found that pick-decompose-apply-stack stages of the DRM method need to be iterated. These iterations are needed because of time and waveform distortions of the pilot traces due to the individual trace statics. The distortions lessen at every external DRM iteration.

scholarly journals Determination of rock densities in the Carpathian-Pannonian Basin lithosphere: based on the CELEBRATION 2000 experiment

2016 ◽  
Vol 46 (4) ◽  
pp. 269-287 ◽  
Author(s):  
Barbora Šimonová ◽  
Miroslav Bielik
Keyword(s):  
Lower Crust ◽  
Pannonian Basin ◽  
Average Density ◽  
East European Craton ◽  
P Wave ◽  
Seismic Velocities ◽  
Upper Crust ◽  
East European ◽  
P Wave Velocity

Abstract The international seismic project CELEBRATION 2000 brought very good information about the P-wave velocity distribution in the Carpathian-Pannonian Basin litosphere. In this paper seismic data were used for transformations of in situ P-wave velocities to in situ densities along all profiles running across the Western Carpathians and the Pannonian Basin: CEL01, CEL04, CEL05, CEL06, CEL09, CEL11 and CEL12. The calculation of rock densities in the crust and lower lithosphere was done by the transformation of seismic velocities to densities using the formulae of Sobolev-Babeyko, Christensen-Mooney and in the lower lithosphere also by Lachenbruch-Morgan’s formula. The density of the upper crust changes significantly in the vertical and horizontal directions, while the interval ranges of the calculated lower crust densities narrow down prominently. The lower lithosphere is the most homogeneous - the intervals of the calculated densities for this layer are already very narrow. The average density of the upper crust (ρ̅ = 2.60 g · cm−3) is the lowest in the Carpathian Foredeep region. On the contrary, the highest density of this layer (ρ̅ = 2.77 g · cm−3) is located in the Bohemian Massif. The average densities ρ̅ of the lower crust vary between 2.90 and 2.98 g · cm−3. The Palaeozoic Platform and the East European Craton have the highest density (ρ̅ = 2.98 g · cm−3 and ρ̅ = 2.97 g · cm−3, respectively). The lower crust density is the lowest (ρ̅ = 2.90 g · cm−3) in the Pannonian Basin. The range of calculated average densities ρ̅ for the lower lithosphere is changed in the interval from 3.35 to 3.40 g · cm−3. The heaviest lower lithosphere can be observed in the East European Craton (ρ̅ = 3.40 g · cm−3). The lower lithosphere of the Transdanubian Range and the Palaeozoic Platform is characterized by the lowest density ρ̅ = 3.35 g · cm−3.

Identification of gas hydrate dissociation from wireline-log data in the Shenhu area, South China Sea

Geophysics ◽  
2012 ◽  
Vol 77 (3) ◽  
pp. B125-B134 ◽  
Author(s):  
Xiujuan Wang ◽  
Myung Lee ◽  
Shiguo Wu ◽  
Shengxiong Yang
Keyword(s):  
Wave Velocity ◽  
Gas Hydrate ◽  
Seismic Data ◽  
Pore Space ◽  
P Wave ◽  
Well Log ◽  
Hydrate Dissociation ◽  
Free Gas ◽  
P Wave Velocity

Wireline logs were acquired in eight wells during China’s first gas hydrate drilling expedition (GMGS-1) in April–June of 2007. Well logs obtained from site SH3 indicated gas hydrate was present in the depth range of 195–206 m below seafloor with a maximum pore-space gas hydrate saturation, calculated from pore water freshening, of about 26%. Assuming gas hydrate is uniformly distributed in the sediments, resistivity calculations using Archie’s equation yielded hydrate-saturation trends similar to those from chloride concentrations. However, the measured compressional (P-wave) velocities decreased sharply at the depth between 194 and 199 mbsf, dropping as low as [Formula: see text], indicating the presence of free gas in the pore space, possibly caused by the dissociation of gas hydrate during drilling. Because surface seismic data acquired prior to drilling were not influenced by the in situ gas hydrate dissociation, surface seismic data could be used to identify the cause of the low P-wave velocity observed in the well log. To determine whether the low well-log P-wave velocity was caused by in situ free gas or by gas hydrate dissociation, synthetic seismograms were generated using the measured well-log P-wave velocity along with velocities calculated assuming both gas hydrate and free gas in the pore space. Comparing the surface seismic data with various synthetic seismograms suggested that low P-wave velocities were likely caused by the dissociation of in situ gas hydrate during drilling.

scholarly journals Estimating uniaxial compressive strength, density and porosity of rocks from the p-wave velocity measurements in-situ and in the laboratory

2021 ◽  
Vol 74 (4) ◽  
pp. 521-528
Author(s):  
André Cezar Zingano ◽  
Paulo Salvadoretti ◽  
Rafael Ubirajara Rocha ◽  
João Felipe Coimbra Leite Costa
Keyword(s):  
Compressive Strength ◽  
Uniaxial Compressive Strength ◽  
Wave Velocity ◽  
P Wave ◽  
Velocity Measurements ◽  
P Wave Velocity

P-wave velocity differences between surface-derived and core samples from the Sulu ultrahigh-pressure terrane: Implications for in situ velocities at great depths

Geology ◽  
2012 ◽  
Vol 40 (7) ◽  
pp. 651-654 ◽  
Author(s):  
Shengsi Sun ◽  
Shaocheng Ji ◽  
Qian Wang ◽  
Matthew Salisbury ◽  
Hartmut Kern
Keyword(s):  
Wave Velocity ◽  
P Wave ◽  
Ultrahigh Pressure ◽  
Core Samples ◽  
P Wave Velocity

High‐resolution crosswell imaging of a west Texas carbonate reservoir: Part 5—Core analysis

Geophysics ◽  
1995 ◽  
Vol 60 (3) ◽  
pp. 712-726 ◽  
Author(s):  
Richard C. Nolen‐Hoeksema ◽  
Zhijing Wang ◽  
Jerry M. Harris ◽  
Robert T. Langan
Keyword(s):  
Wave Velocity ◽  
Field Experiments ◽  
Velocity Ratio ◽  
Carbonate Reservoir ◽  
P Wave ◽  
Core Analysis ◽  
S Wave ◽  
Traveltime Tomography ◽  
West Texas ◽  
P Wave Velocity

We conducted a core analysis program to provide supporting data to a series of crosswell field experiments being carried out in McElroy Field by Stanford University’s Seismic Tomography Project. The objective of these experiments is to demonstrate the use of crosswell seismic profiling for reservoir characterization and for monitoring [Formula: see text] flooding. For these west Texas carbonates, we estimate that [Formula: see text] saturation causes P‐wave velocity to change by −1.9% (pooled average, range = −6.3 to +0.1%), S‐wave velocity by +0.6% (range = 0 to 2.7%), and the P‐to‐S velocity ratio by −2.4% (range = −6.4 to −0.3%). When we compare these results to the precisions we can expect from traveltime tomography (about ±1% for P‐ and S‐wave velocity and about ±2% for the P‐to‐S velocity ratio), we conclude that time‐lapse traveltime tomography is sensitive enough to resolve changes in the P‐wave velocity, S‐wave velocity, and P‐to‐S velocity ratio that result from [Formula: see text] saturation. We concentrated here on the potential for [Formula: see text] saturation to affect seismic velocities. The potential for [Formula: see text] saturation to affect other seismic properties, not discussed here, may prove to be more significant (e.g., P‐wave and S‐wave impedance).

scholarly journals Analyzing the early structural build-up of accelerated cement pastes

2021 ◽  
Vol 54 (2) ◽  
Author(s):  
Tobias Dorn ◽  
Tamino Hirsch ◽  
Dietmar Stephan
Keyword(s):  
P Wave ◽  
Major Influence ◽  
X Ray Diffraction ◽  
Cement Pastes ◽  
Ettringite Formation ◽  
Needle Penetration ◽  
P Wave Velocity ◽  
Early Reaction ◽  
Main Factor

AbstractExtrusion-based additive manufacturing imposes high requirements on the material stability right after the extrusion. Therefore, a thorough understanding of the chemical reactions that determine the early reduction in processability is necessary. Accelerators are especially considered here, which have a major influence on the early reaction. This study contributes to these issues by analyzing the influence of 0.1 wt% TEA (triethanolamine) and 2.0 wt% Ca(NO3)2 on the hydration of two CEM I 52.5 R. The hydration was analyzed by isothermal heat flow calorimetry and in-situ X-ray diffraction. Vicat needle penetration, a penetrometer of own design, and ultrasonic P-wave velocity development were used to monitor the early change in workability. The obtained results indicate that ettringite formation is the main factor influencing workability during the first 60 min of hydration. Afterwards, the influence of ettringite is exceeded by the formation of C–S–H. Ca(NO3)2 was shown to enhance C–S–H formation and had no significant effect on the workability during the first 60 min of hydration while rapidly decreasing workability during the induction period. TEA was shown to increase brownmillerite dissolution and ettringite formation and, consequently, promote the workability loss during the initial hydration period. The time of initial and final setting determined by the Vicat needle test depended mainly on the formation of C–S–H phases.

scholarly journals In‐situ, high‐frequency P-wave velocity measurements within 1 m of the earth’s surface

Geophysics ◽  
1999 ◽  
Vol 64 (2) ◽  
pp. 323-325 ◽  
Author(s):  
Gregory S. Baker ◽  
Don W. Steeples ◽  
Chris Schmeissner
Keyword(s):  
Wave Velocity ◽  
Speed Of Sound ◽  
P Wave ◽  
Sound Waves ◽  
Low Velocity ◽  
Near Surface ◽  
Wave Velocities ◽  
P Wave Velocity ◽  
P Wave Velocities

Seismic P-wave velocities in near‐surface materials can be much slower than the speed of sound waves in air (normally 335 m/s or 1100 ft/s). Difficulties often arise when measuring these low‐velocity P-waves because of interference by the air wave and the air‐coupled waves near the seismic source, at least when gathering data with the more commonly used shallow P-wave sources. Additional problems in separating the direct and refracted arrivals within ∼2 m of the source arise from source‐generated nonlinear displacement, even when small energy sources such as sledgehammers, small‐caliber rifles, and seismic blasting caps are used. Using an automotive spark plug as an energy source allowed us to measure seismic P-wave velocities accurately, in situ, from a few decimeters to a few meters from the shotpoint. We were able to observe three distinct P-wave velocities at our test site: ∼130m/s, 180m/s, and 300m/s. Even the third layer, which would normally constitute the first detected layer in a shallow‐seismic‐refraction survey, had a P-wave velocity lower than the speed of sound in air.

On: “Reflection and transmission of plane compressional waves by R. D. Tooley, T. W. Spencer, and H. F. Sagori, (GEOPHYSICS, 30, 552, April, 1965).

Geophysics ◽  
1984 ◽  
Vol 49 (12) ◽  
pp. 2195-2195 ◽  
Author(s):  
L. R. Denham ◽  
R. A. R. Palmeira
Keyword(s):  
Wave Velocity ◽  
Seismic Wave ◽  
Velocity Ratio ◽  
The Other ◽  
P Wave ◽  
Angle Of Incidence ◽  
Reflection And Transmission ◽  
Compressional Waves ◽  
P Wave Velocity ◽  
Simplified Formula

One of the most widely reproduced figures for the partition of energy of a seismic wave at an interface between two media is from Tooley et al.’s 1965 paper. Recently we referred to this paper and noticed some anomalies in the curves shown there in Figure 11. We looked at the values of the curves at zero angle of incidence and noted that while the values of a P-wave velocity ratio of 2, 3, and 4 agree with the values given by the simplified formula {(V2−V1)/(V2+V1)}**2, the values shown for the P-wave velocity ratios less than 1 do not agree. For a ratio of 0.25, the coefficient at zero angle should be 0.36; it is shown as about 0.72. A ratio of 0.5 should give a coefficient of 0.11; 0.52 is shown. There are similar discrepancies for all the other ratios less than 1.

Sign in / Sign up

Export Citation Format

Share Document