Viscoelasticity of precompacted unconsolidated sand with viscous cement

Geophysics ◽  
2006 ◽  
Vol 71 (2) ◽  
pp. T31-T40 ◽  
Author(s):  
Klaus C. Leurer ◽  
Jack Dvorkin
Keyword(s):  
Viscous Fluid ◽  
Contact Stiffness ◽  
Frequency Dispersion ◽  
P Wave ◽  
Seismic Velocities ◽  
Effective Pressure ◽  
S Wave ◽  
Normal Stiffness ◽  
Cemented Sand ◽  
Wave Velocities

The elastic properties of sand strongly depend on the grains’ contact stiffness, which can be increased significantly by solid matter and, depending on frequency, viscous fluid acting as contact cement. To calculate seismic velocities in precompacted fluid-cemented sand, we examine how a small amount of viscous fluid at the grain contacts influences their normal and tangential stiffnesses as a function of effective pressure. Using the Hertz-Mindlin approach and considering oscillatory loading in addition to precompaction of a combination of two elastic spheres, we extend the dry-contact elastic theory by a viscoelastic formulation. Here, we describe the radial flow of the fluid cement induced by the oscillations of the grains’ surfaces around the direct contact, a process that leads to a complex normal stiffness and stiffness/frequency dispersion. In the resulting combined model, the low-frequency real part of the complex normal stiffness identical is to the original Hertz-Mindlin expression. The magnitude of the dispersion is governed by the amount of viscous cement; magnitude decreases as effective pressure increases. The frequency of the maximum imaginary part of the normal stiffness is determined mainly by cement viscosity and contact geometry. The tangential contact stiffness virtually is not influenced by the viscous fluid. Comparison of predicted results with data from pulse transmission experiments (500 kHz) on glass beads with two different fluids shows an excellent fit in P-wave velocities [Formula: see text], whereas S-wave velocities [Formula: see text] are systematically overestimated by the model. The experimental results confirm, however, the predicted change with effective pressure in the [Formula: see text] ratio for both examined cases as well as reflect the predicted increase in [Formula: see text] and [Formula: see text], respectively, between the two cases. This implies that our viscoelastic formulation represents a reasonable way to describe the role of viscous cement in sand.

Seismic Velocities and Attenuation from Borehole Measurements near the Parkfield Prediction Zone, Central California

1989 ◽  
Vol 5 (3) ◽  
pp. 513-537 ◽  
Author(s):  
James F. Gibbs ◽  
Edward F. Roth
Keyword(s):  
P Wave ◽  
Depth Interval ◽  
Seismic Velocities ◽  
S Wave ◽  
Wave Source ◽  
Central California ◽  
Unconsolidated Sands ◽  
Wave Velocities ◽  
Average Shear ◽  
Parkfield Earthquake

Shear (S)- and compressional (P)- wave velocities were measured to a depth of 195 m in a borehole near the San Andreas fault where a recurrence of a moderate Parkfield earthquake is predicted. S-wave velocities determined from orthogonal directions of the S-wave source show velocity differences of approximately 20 percent. An average shear-wave Q of 4 was determined in relatively unconsolidated sands and gravels of the Paso Robles Formation in the depth interval 57.5-102.5 m.

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 LOCAL LOW PRESSURE AREAS IN ANTICLINE STRUCTURES

2015 ◽  
Vol 33 (2) ◽  
Author(s):  
Boris P. Sibiryakov ◽  
Lourenildo W.B. Leite ◽  
Egor P. Sibiryakov ◽  
Wildney W.S. Vieira
Keyword(s):  
Oil And Gas ◽  
Target Surface ◽  
Sedimentary Basins ◽  
P Wave ◽  
Necessary Condition ◽  
Seismic Velocities ◽  
S Wave ◽  
Wave Velocities ◽  
Horizontal Boundary ◽  
Constitutive Parameters

ABSTRACT. In order to localize low pressures zones in sedimentary basins for oil and gas exploration, it is necessary to know P and S wave velocities for medium. Strictly speaking, we need to know the rock densities for all layers, and in addition there are many correlation tables between seismic velocities and densities; besides, density is a parameter admitted to change slowly with depth to the top of the target interface. P wave velocities are considered a conventional asset, and S wave velocities can be obtained from special field survey, in particular from converted P-S waves registered by VSP technology, and by petrophysical measurements. The theory in this paper deals with stress prediction in the subsurface, and takes in consideration the constitutive parameters (density and Lame’s), and the geometry of the reservoir target surface. The model does not separate the different contributions (porosity, fluids) to the rock velocities controlled by the constitutive parameters. It is not a necessary condition that an anticline be a potential structure for oil and gas accumulation. This role can be played by horizontal structures if there is a positive γ = VS VPratio discontinuity, or a negative discontinuity of the Poisson, σ, ratio across the horizontal boundary. These conditions are responsible for producing a pressure discontinuity, such that beneath the boundary there will be a sufficiently lower pressure zone than above the boundary. In this case, the lower horizontal boundary is saidto be an attractor surface for fluids of the any kind; in the opposite case, this boundary does not have fluid attractor properties.Keywords: seismic structured media, porous media, anticline structures, pressure prediction.RESUMO. Com o objetivo de localizar zonas de baixa pressão em bacias sedimentares voltadas à exploração de óleo e gás, é necessário conhecer as velocidades das ondas P e S para o meio. Mais especificamente, precisamos conhecer a densidade das rochas em todas as camadas, e aditivamente existem várias tabelas de correlação entre velocidade e densidade das rochas; além disso, é um parâmetro que varia lentamente com a profundidade até o topo da interface-alvo. A velocidade das ondas P é considerada uma informação convencional, e a velocidade das ondas S pode ser obtida por levantamentos especiais de campo, em particular a partir da conversão P-S registrada por tecnologia VSP, e por medidas petrofísicas. A teoria deste trabalho trata da predição de tensão na subsuperfície, e leva em consideração os parâmetros constitutivos (densidade e de Lamé), e a topografia em superfície do reservatório-alvo. O modelo não separa as diferentes contribuições (porosidade e fluidos) para estabelecer as velocidades nas rochas controladas pelos parâmetros constitutivos. Não é uma condição necessária que uma superfície anticlinal seja uma estrutura potencial para o acúmulo de óleo e gás. Esta condição pode ser representada por uma superfície horizontal, se existir uma discontinuidade na razão γ = VS VP, ou uma descontinuidade negativa na razão de Poisson, σ, através da superfície. Estas condições são responsáveis por produzir uma descontinuidade de pressão, de tal forma que abaixo da interface existirá uma zona de pressão mais baixa do que há acima da mesma. Neste caso, a parte inferior da superfície é considerada como um atrativo de fluidos de qualquer tipo; e no caso oposto, esta superfície não ´e dotada de propriedades de atração de fluidos.Palavras-chave: meios sísmicos estruturados, meios porosos, meios fraturados.

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.

Seismic velocities in transversely isotropic media, II

Geophysics ◽  
1980 ◽  
Vol 45 (1) ◽  
pp. 3-17 ◽  
Author(s):  
Franklyn K. Levin
Keyword(s):  
Wave Velocity ◽  
Transversely Isotropic ◽  
P Wave ◽  
Seismic Velocities ◽  
Sh Wave ◽  
S Wave ◽  
Wave Velocities ◽  
P Wave Velocity ◽  
Surface Spread ◽  
Isotropic Solids

P‐wave, SV‐wave, and SH‐wave velocities are computed for transversely isotropic solids formed from two isotropic solids. The combinations are shale‐sandstone and shale‐limestone solids of an earlier paper (Levin, 1979), but one velocity of the nonshale component is allowed to vary over the range of Poisson’s ratios σ = 0 to σ = 0.45, i.e., from a rigid solid to a near‐liquid. When the S‐wave velocity of either the sandstone or limestone is varied, the ratio of horizontal P‐wave velocity to vertical P‐wave velocity goes through a maximum as σ increases and subsequently falls to values less than unity as σ approaches 0.5. The P‐wave velocity that would be found with a short surface spread also goes through a maximum and, at σ = 0.5, is less than the P‐wave velocity of either isotropic component. SV‐wave velocities found for data from a short spread are unreasonably large; SH‐wave velocities decrease monotonically as σ increases, but the ratio of horizontal SH‐wave velocity to vertical SH‐wave velocity goes through a minimum of unity.

Shallow velocity structure and poisson's ratio at the Tarzana, California, strong-motion accelerometer site

1996 ◽  
Vol 86 (6) ◽  
pp. 1704-1713 ◽  
Author(s):  
R. D. Catchings ◽  
W. H. K. Lee
Keyword(s):  
Urban Areas ◽  
Velocity Structure ◽  
Strong Motion ◽  
P Wave ◽  
S Wave ◽  
Near Surface ◽  
Velocity Models ◽  
Wave Velocities ◽  
Poisson's Ratios ◽  
Poisson’S Ratios

Abstract The 17 January 1994, Northridge, California, earthquake produced strong ground shaking at the Cedar Hills Nursery (referred to here as the Tarzana site) within the city of Tarzana, California, approximately 6 km from the epicenter of the mainshock. Although the Tarzana site is on a hill and is a rock site, accelerations of approximately 1.78 g horizontally and 1.2 g vertically at the Tarzana site are among the highest ever instrumentally recorded for an earthquake. To investigate possible site effects at the Tarzana site, we used explosive-source seismic refraction data to determine the shallow (<70 m) P-and S-wave velocity structure. Our seismic velocity models for the Tarzana site indicate that the local velocity structure may have contributed significantly to the observed shaking. P-wave velocities range from 0.9 to 1.65 km/sec, and S-wave velocities range from 0.20 and 0.6 km/sec for the upper 70 m. We also found evidence for a local S-wave low-velocity zone (LVZ) beneath the top of the hill. The LVZ underlies a CDMG strong-motion recording site at depths between 25 and 60 m below ground surface (BGS). Our velocity model is consistent with the near-surface (<30 m) P- and S-wave velocities and Poisson's ratios measured in a nearby (<30 m) borehole. High Poisson's ratios (0.477 to 0.494) and S-wave attenuation within the LVZ suggest that the LVZ may be composed of highly saturated shales of the Modelo Formation. Because the lateral dimensions of the LVZ approximately correspond to the areas of strongest shaking, we suggest that the highly saturated zone may have contributed to localized strong shaking. Rock sites are generally considered to be ideal locations for site response in urban areas; however, localized, highly saturated rock sites may be a hazard in urban areas that requires further investigation.

scholarly journals A method of geo-acoustic parameter inversion in shallow sea by the Bayesian theory and the acoustic pressure field

2019 ◽  
Vol 283 ◽  
pp. 06003
Author(s):  
Guangxue Zheng ◽  
Hanhao Zhu ◽  
Jun Zhu
Keyword(s):  
Sound Pressure ◽  
Pressure Field ◽  
Wave Attenuation ◽  
Acoustic Parameter ◽  
Bayesian Theory ◽  
P Wave ◽  
Shallow Sea ◽  
S Wave ◽  
Wave Velocities ◽  
Parameter Inversion

A method of geo-acoustic parameter inversion based on the Bayesian theory is proposed for the acquisition of acoustic parameters in shallow sea with the elastic seabed. Firstly, the theoretical prediction value of the sound pressure field is calculated by the fast field method (FFM). According to the Bayesian theory, we establish the misfit function between the measured sound pressure field and the theoretical pressure field. It is under the assumption of Gaussian data errors which are in line with the likelihood function. Finally, the posterior probability density (PPD) of parameters is given as the result of inversion. Our research is conducted in the light of Metropolis sample rules. Apart from numerical simulations, a scaled model experiment has been taken in the laboratory tank. The results of numerical simulations and tank experiments show that sound pressure field calculated by the result of inversion is consistent with the measured sound pressure field. Besides, s-wave velocities, p-wave velocities and seafloor density have fewer uncertainties and are more sensitive to complex sound pressure than s-wave attenuation and p-wave attenuation. The received signals calculated by inversion results are keeping with received signals in the experiment which verify the effectiveness of this method.

EXTRACTING SUB-SURFACE INFORMATION FROM LONG OFFSET P-WAVE DATA—A CHEAPER ALTERNATIVE TO MULTICOMPONENT OBC FOR EXPLORATION?

2002 ◽  
Vol 42 (1) ◽  
pp. 627
Author(s):  
R.G. Williams ◽  
G. Roberts ◽  
K. Hawkins
Keyword(s):  
Seismic Energy ◽  
Higher Order ◽  
P Wave ◽  
S Wave ◽  
Additional Information ◽  
Wave Data ◽  
Wave Velocities ◽  
Avo Inversion ◽  
Long Offset ◽  
Multicomponent Data

Seismic energy that has been mode converted from pwave to s-wave in the sub-surface may be recorded by multi-component surveys to obtain information about the elastic properties of the earth. Since the energy converted to s-wave is missing from the p-wave an alternative to recording OBC multi-component data is to examine p-wave data for the missing energy. Since pwave velocities are generally faster than s-wave velocities, then for a given reflection point the converted s-wave signal reaches the surface at a shorter offset than the equivalent p-wave information. Thus, it is necessary to record longer offsets for p-wave data than for multicomponent data in order to measure the same information.A non-linear, wide-angle (including post critical) AVO inversion has been developed that allows relative changes in p-wave velocities, s-wave velocities and density to be extracted from long offset p-wave data. To extract amplitudes at long offsets for this inversion it is necessary to image the data correctly, including correcting for higher order moveout and possibly anisotropy if it is present.The higher order moveout may itself be inverted to yield additional information about the anisotropy of the sub-surface.

EFFECT OF ALIGNMENT ON THE QUALITY OF BENDER ELEMENT PROCEDURE

2015 ◽  
Vol 76 (2) ◽  
Author(s):  
Badee Alshameri ◽  
Ismail Bakar ◽  
Aziman Madun ◽  
Edy Tonnizam Mohamad
Keyword(s):  
Vertical Axis ◽  
Poor Quality ◽  
P Wave ◽  
Secondary Wave ◽  
Horizontal Distance ◽  
Bender Element ◽  
S Wave ◽  
Testing Procedures ◽  
Wave Velocities

One of the main geophysical tools (seismic tools) in the laboratory is the bender element. This tool can be used to measure some dynamic soil properties (e.g. shear and Young’s modulus). However, even if it relatively simple to use the bender element, inconsistent testing procedures can cause poor quality in the bender element data. One of the bender element procedure that always neglected is the alignment (different positions of bender element receiver to the transmitter in the vertical axis). The alignment effect was evaluated via changing the horizontal distance between transmitter and receiver starting from 0 to 110 mm for two sizes of the sample's thickness (i.e. 63.17 mm and 91.51 mm). Five methods were applied to calculate the travel times. Those methods were as the following: visually, first-peak, maximum-peak, CCexcel and CCGDS. In general, the experiments indicated uncertain results for both of the P-wave (primary wave) and S-wave (secondary wave) velocities at zone of Dr:D above 0.5:1 (where Dr is the horizontal distance of the receiver from the vertical axis and D is the thickness of the sample). On the other hand, both the visual and first-peak methods show the wave velocities results are higher than obtained from other methods. However, the ratio between the amplitude of transmitter signals to receiver amplitude signal was taken to calculate the damping-slope of the P-wave and S-wave. Thus the results from damping slope show steeply slope when the ratio of  Dr:D is above 0.5:1 compare with gentle slope below ratio 0.5:1 at the sample with thickness equal to 91.51 mm, while there is no variation at a slope in sample with thickness equal to 63.17 mm.

Sign in / Sign up

Export Citation Format

Share Document