White's model for wave propagation in partially saturated rocks: Comparison with poroelastic numerical experiments

Geophysics ◽  
2003 ◽  
Vol 68 (4) ◽  
pp. 1389-1398 ◽  
Author(s):  
José M. Carcione ◽  
Hans B. Helle ◽  
Nam H. Pham

We use a poroelastic modeling algorithm to compute numerical experiments of wave propagation in White's partial saturation model. The results are then compared to the theoretical predictions. The model consists of a homogeneous sandstone saturated with brine and spherical gas pockets. White's theory predicts a relaxation mechanism, due to pressure equilibration, causing attenuation and velocity dispersion of the wavefield. We vary gas saturation either by increasing the radius of the gas pocket or by increasing the density of gas bubbles. Despite that the modeling is two dimensional and interaction between the gas pockets is neglected in White's model, the numerical results show the trends predicted by the theory. In particular, we observe a similar increase in velocity at high frequencies (and low permeabilities). Furthermore, the behavior of the attenuation peaks versus water saturation and frequency is similar to that of White's model. The modeling results show more dissipation and higher velocities than White's model due to multiple scattering and local fluid‐flow effects. The conversion of fast P‐wave energy into dissipating slow waves at the patches is the main mechanism of attenuation. Differential motion between the rock skeleton and the fluids is highly enhanced by the presence of fluid/fluid interfaces and pressure gradients generated through them.

2017 ◽  
Vol 144 ◽  
pp. 94-103 ◽  
Author(s):  
Jie Liu ◽  
Dameng Liu ◽  
Yidong Cai ◽  
Quan Gan ◽  
Yanbin Yao

Geophysics ◽  
2001 ◽  
Vol 66 (6) ◽  
pp. 1838-1842 ◽  
Author(s):  
C. M. Schmeissner ◽  
K. T. Spikes ◽  
D. W. Steeples

Ultrashallow seismic reflection surveys require dense spatial sampling during data acquisition, which increases their cost. In previous efforts to find ways to reduce these costs, we connected geophones rigidly to pieces of channel iron attached to a farm implement. This method allowed us to plant the geophones in the ground quickly and automatically. The rigidly interconnected geophones used in these earlier studies detected first‐arrival energy along with minor interfering seismic modes, but they did not detect seismic reflections. To examine further the feasibility of developing rigid geophone emplacement systems to detect seismic reflections, we experimented with four pieces of channel iron, each 2.7 m long and 10 cm wide. Each segment was equipped with 18 geophones rigidly attached to the channel iron at 15‐cm intervals, and the spikes attached to all 18 geophones were pushed into the ground simultaneously. The geophones detected both refracted and reflected energy; however, no significant signal distortion or interference attributable to the rigid coupling of the geophones to the channel iron was observed in the data. The interfering seismic modes mentioned from the previous experiments were not detected, nor was any P‐wave propagation noted within the channel iron. These results show promise for automating and reducing the cost of ultrashallow seismic reflection and refraction surveys.


Author(s):  
Suresh Dande ◽  
◽  
Robert R. Stewart ◽  
Nikolay Dyaur ◽  
◽  
...  

Laboratory physical models play an important role in understanding rock properties and wave propagation, both theoretically and at the field scale. In some cases, 3D-printing technology can be adopted to construct complex rock models faster, more inexpensively, and with more specific features than previous model-building techniques. In this study, we use 3D-printed rock models to assist in understanding the effects of various fluids (air, water, engine oil, crude oil, and glycerol) on the models’ elastic properties. We first used a 3D-printed, 1-in. cube-shaped layered model. This model was created with a 6% primary porosity and a bulk density of 0.98 g/cc with VTI anisotropy. We next employed a similar cube but with horizontal inclusions embedded in the layered background, which contributed to its total 24% porosity (including primary porosity). For air to liquid saturation, P-velocities increased for all liquids in both models, with the highest increase being with glycerol (57%) and an approximately 45% increase for other fluids in the inclusion model. For the inclusion model (dry and saturated), we observed a greater difference between two orthogonally polarized S-wave velocities (Vs1 and Vs2) than between two P-wave velocities (VP0 and VP90). We attribute this to the S2-wave (polarized normal to both the layering and the plane of horizontal inclusions), which appears more sensitive to horizontal inclusions than the P-wave. For the inclusion model, Thomsen’s P-wave anisotropic parameter (ɛ) decreased from 26% for the air case to 4% for the water-saturated cube and to 1% for glycerol saturation. The small difference between the bulk modulus of the frame and the pore fluid significantly reduces the velocity anisotropy of the medium, making it almost isotropic. We compared our experimental results with theory and found that predictions using Schoenberg’s linear slip theory combined with Gassmann’s anisotropic equation were closer to actual measurements than Hudson’s isotropic calculations. This work provides insights into the usefulness of 3D-printed models to understand elastic rock properties and wave propagation under various fluid saturations.


2014 ◽  
Vol 988 ◽  
pp. 502-507 ◽  
Author(s):  
Shao Bo Chai ◽  
Jian Chun Li ◽  
Hai Bo Li ◽  
Ya Qun Liu

According to the displacement discontinuity method and the conservation of momentum at the wave fronts, analysis for cylindrical P-wave propagation across a linear elastic rock joint is carried out. Considering the energy variation for wave propagation in one medium, the wave propagation equation was derived and expressed in an iterative form. The transmission and reflection coefficients are then obtained from the equation. By verification, the results agree very well with those from the existing results.


2016 ◽  
Vol 40 ◽  
pp. 1660061
Author(s):  
Vadim Baru

The recent measurements of the reactions [Formula: see text] and [Formula: see text] by the ANKE collaboration at COSY are analyzed with the focus on the p-wave pion production amplitudes. These amplitudes are known to provide an important connection between [Formula: see text] and other low-energy few-nucleon reactions. The results of the recent partial wave analysis of the ANKE data are discussed and compared with the theoretical predictions.


2021 ◽  
Author(s):  
Thomas Möller ◽  
Wolfgang Friederich

<p>Modeling waveforms of teleseismic body waves requires the solution of the seismic wave equation in the entire Earth. Since fully-numerical 3D simulations on a global scale with periods of a few seconds are far too computationally expensive, we resort to a hybrid approach in which fully-numerical 3D simulations are performed only within the target region and wave propagation through the rest of the Earth is modeled using methods that are much faster but apply only to spherically symmetric Earth models.</p><p>We present a hybrid method that uses GEMINI to compute wave fields for a spherically symmetric Earth model up to the boundaries of a regional box. The wavefield is injected at the boundaries, where wave propagation is continued using SPECFEM-Cartesian. Inside the box, local heterogeneities in the velocity distribution are allowed, which can cause scattered and reflected waves. To prevent these waves from reflecting off the edges of the box absorbing boundary conditions are specifically applied to these parts of the wavefields. They are identified as the difference between the wavefield calculated with SPECFEM at the edges and the incident wavefield.</p><p>The hybrid method is applied to a target region in and around the Alps as a test case. The region covers an area of 1800 by 1350 km centered at 46.2°N and 10.87°E and includes crust and mantle to a depth of 600 km. We compare seismograms with a period of up to ten seconds calculated with the hybrid method to those calculated using GEMINI only for identical 1D earth models. The comparison of the seismograms shows only very small differences and thus validates the hybrid method. In addition, we demonstrate the potential of the method by calculating seismograms where the 1D velocity model inside the box is replaced by a velocity model generated using P-wave traveltime tomography.</p>


Sign in / Sign up

Export Citation Format

Share Document