High-frequency seismic response during permeability reduction due to biopolymer clogging in unconsolidated porous media

Geophysics ◽  
2013 ◽  
Vol 78 (6) ◽  
pp. EN117-EN127 ◽  
Author(s):  
Tae-Hyuk Kwon ◽  
Jonathan B. Ajo-Franklin

The accumulation of biopolymers in porous media, produced by stimulating either indigenous bacteria or artificially introduced microbes, readily blocks pore throats and can effectively reduce bulk permeability. Such a microbial clogging treatment can be used for selective plugging of permeable zones in reservoirs and is considered a potentially promising approach to enhance sweep efficiency for microbial enhanced oil recovery (MEOR). Monitoring in situ microbial growth, biopolymer formation, and permeability reduction in the reservoir is critical for successful application of this MEOR approach. We examined the feasibility of using seismic signatures (P-wave velocity and attenuation) for monitoring the in situ accumulation of insoluble biopolymers in unconsolidated sediments. Column experiments, which involved stimulating the sucrose metabolism of Leuconostoc mesenteroides and production of the biopolymer dextran, were performed while monitoring changes in permeability and seismic response using the ultrasonic pulse transmission method. We observed that L. mesenteroides produced a viscous biopolymer in sucrose-rich media. Accumulated dextran, occupying 4%–6% pore volume after [Formula: see text] days of growth, reduced permeability more than one order of magnitude. A negligible change in P-wave velocity was observed, indicating no or minimal change in compressive stiffness of the unconsolidated sediment during biopolymer formation. The amplitude of the P-wave signals decreased [Formula: see text] after [Formula: see text] days of biopolymer production; spectral ratio analysis in the 0.4–0.8-MHz band showed an approximate 30%–50% increase in P-wave attenuation ([Formula: see text]) due to biopolymer production. A flow-induced loss mechanism related to the combined grain/biopolymer structure appeared to be the most plausible mechanism for causing the observed increase in P-wave attenuation in the ultrasonic frequency range. Because permeability reduction is also closely linked to biopolymer volume, P-wave attenuation in the ultrasonic frequency range appears to be an effective indicator for monitoring in situ biopolymer accumulation and permeability reduction and could provide a useful proxy for regions with altered transport properties.

Geophysics ◽  
2014 ◽  
Vol 79 (4) ◽  
pp. D205-D216 ◽  
Author(s):  
Xinding Fang ◽  
Michael C. Fehler ◽  
Arthur Cheng

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.


Geophysics ◽  
1989 ◽  
Vol 54 (1) ◽  
pp. 76-81 ◽  
Author(s):  
D. Goldberg ◽  
B. Zinszner

We computed compressional‐wave velocity [Formula: see text] and attenuation [Formula: see text] from sonic log waveforms recorded in a cored, 30 m thick, dolostone reservoir; using cores from the same reservoir, laboratory measurements of [Formula: see text] and [Formula: see text] were also obtained. We used a resonant bar technique to measure extensional and shear‐wave velocities and attenuations in the laboratory, so that the same frequency range as used in sonic logging (5–25 kHz) was studied. Having the same frequency range avoids frequency‐dependent differences between the laboratory and in‐situ measurements. Compressional‐wave attenuations at 0 MPa confining pressure, calculated on 30 samples, gave average [Formula: see text] values of 17. Experimental and geometrical errors were estimated to be about 5 percent. Measurements at elevated effective pressures up to 30 MPa on selected dolostone samples in a homogeneous interval showed mean [Formula: see text] and [Formula: see text] to be approximately equal and consistently greater than 125. At effective stress of 20 MPa and at room temperature, the mean [Formula: see text] over the dolostone interval was 87, a minimum estimate for the approximate in‐situ conditions. We computed compressional‐wave attenuation from sonic log waveforms in the 12.5–25 kHz frequency band using the slope of the spectral ratio of waveforms recorded 0.914 m and 1.524 m from the source. Average [Formula: see text] over the interval was 13.5, and the mean error between this value and the 95 percent confidence interval of the slope was 15.9 percent. The laboratory measurements of [Formula: see text] under elevated pressure conditions were more than five times greater than the mean in‐situ values. This comparison shows that additional extrinsic losses in the log‐derived measurement of [Formula: see text], such as scattering from fine layers and vugs or mode conversion to shear energy dissipating radially from the borehole, dominate the apparent attenuation.


Geophysics ◽  
1986 ◽  
Vol 51 (12) ◽  
pp. 2193-2199 ◽  
Author(s):  
G. W. Purnell

The velocity and attenuation of a wave transmitted through a two‐phase material are functions of the material’s composition. In physical model experiments, I used suspensions of grains in a silicone rubber matrix to reduce or avoid uncertainties about framework elastic constants, porosity, and permeability that result from using fluid‐saturated grain frameworks. I varied the composition to produce materials that are useful in physical seismic modeling. In the tested suspensions, ultrasonic P-wave velocity, velocity dispersion, and attenuation all increase with grain concentration and frequency. I compared seven published mathematical models for wave propagation in two‐phase media. One given by Mehta most closely agrees with the P-wave velocities I observed. The agreement is sufficiently close to merit use of Mehta’s model in the design of physical model materials. The observed P-wave attenuation generally increases approximately linearly with frequency. This approximate linearity leads to reliable constant-Q estimates, ranging from 187 to 16 for grain concentrations from 0 to 0.49. I conclude that relative motion between the grains and the rubber matrix contributes most of the observed attenuation at lower concentrations, whereas scattering losses become much more important at higher concentrations and frequencies.


Geophysics ◽  
2012 ◽  
Vol 77 (3) ◽  
pp. B125-B134 ◽  
Author(s):  
Xiujuan Wang ◽  
Myung Lee ◽  
Shiguo Wu ◽  
Shengxiong Yang

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.


2021 ◽  
Vol 74 (4) ◽  
pp. 521-528
Author(s):  
André Cezar Zingano ◽  
Paulo Salvadoretti ◽  
Rafael Ubirajara Rocha ◽  
João Felipe Coimbra Leite Costa

Geology ◽  
2012 ◽  
Vol 40 (7) ◽  
pp. 651-654 ◽  
Author(s):  
Shengsi Sun ◽  
Shaocheng Ji ◽  
Qian Wang ◽  
Matthew Salisbury ◽  
Hartmut Kern

Geophysics ◽  
2015 ◽  
Vol 80 (1) ◽  
pp. D51-D63 ◽  
Author(s):  
Zizhen Wang ◽  
Ruihe Wang ◽  
Ralf J. Weger ◽  
Tianyang Li ◽  
Feifei Wang

The relationship between P-wave velocity and porosity in carbonate rocks shows a high degree of variability due to the complexity of the pore structure. This variability introduces high uncertainties to seismic inversion, amplitude variation with offset analysis, porosity estimation, and pore-pressure prediction based on velocity data. Elastic wave propagation in porous media is numerically modeled on the pore scale to investigate the effects of pore structure on P-wave velocities in carbonate rocks. We built 2D models of porous media using pore structure information and the similarity principle. Then, we simulated normal incidence wave propagation using finite element analysis. Finally, the velocity was determined from received modeled signals by means of crosscorrelation. The repeatability and accuracy of this modeling process was verified carefully. Based on the modeling results, a simple formulation of Sun’s frame flexibility factor ([Formula: see text]), aspect ratio (AR, the ratio of the major axis to the minor axis), and pore density was developed. The numerical simulation results indicated that the P-wave velocity increases as a power function as the AR increases. Pores with small AR ([Formula: see text]) or large [Formula: see text] created softening effects that decrease P-wave velocity significantly. The P-wave velocity of carbonate rocks was dispersive; it depends on the ratio of the wavelength to pore size ([Formula: see text]). Such scale-dependent dispersion was more evident for carbonate rocks with higher porosity, lower AR, and/or lower P-wave impedance of pore fluids. The P-wave velocity of carbonate rocks with complicated pore geometries (low AR, high [Formula: see text], small [Formula: see text]) was much lower than that of rocks with simple pore geometries (high AR, small [Formula: see text], large [Formula: see text]) at low and high [Formula: see text]. The pore-scale modeling of elastic wave properties of porous rocks may explain the poor velocity-porosity correlation in carbonate rocks.


Geophysics ◽  
1999 ◽  
Vol 64 (2) ◽  
pp. 323-325 ◽  
Author(s):  
Gregory S. Baker ◽  
Don W. Steeples ◽  
Chris Schmeissner

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.


2020 ◽  
Vol 222 (1) ◽  
pp. 560-571
Author(s):  
Lingli Gao ◽  
Yudi Pan ◽  
Thomas Bohlen

SUMMARY 2-D full-waveform inversion (FWI) of shallow-seismic wavefields has recently become a novel way to reconstruct S-wave velocity models of the shallow subsurface with high vertical and lateral resolution. In most applications, seismic wave attenuation is ignored or considered as a passive modelling parameter only. In this study, we explore the feasibility and performance of multiparameter viscoelastic 2-D FWI in which seismic velocities and attenuation of P and S waves, respectively, and mass density are inverted simultaneously. Synthetic reconstruction experiments reveal that multiple crosstalks between all viscoelastic material parameters may occur. The reconstruction of S-wave velocity is always robust and of high quality. The parameters P-wave velocity and density exhibit weaker sensitivity and can be reconstructed more reliably by multiparameter viscoelastic FWI. Anomalies in S-wave attenuation can be recovered but with limited resolution. In a field-data application, a small-scale refilled trench is nicely delineated as a low P- and S-wave velocity anomaly. The reconstruction of P-wave velocity is improved by the simultaneous inversion of attenuation. The reconstructed S-wave attenuation reveals higher attenuation in the shallow weathering zone and weaker attenuation below. The variations in the reconstructed P- and S-wave velocity models are consistent with the reflectivity observed in a ground penetrating radar (GPR) profile.


2019 ◽  
Vol 221 (1) ◽  
pp. 60-76 ◽  
Author(s):  
Yongjia Song ◽  
Hengshan Hu ◽  
Bo Han

SUMMARY Understanding the wave propagation in fluid-saturated cracked rocks is important for detecting and characterizing cracked reservoirs and fault zones with applications in geomechanics, hydrogeology, exploration geophysics and reservoir engineering. In sedimentary rocks, microscopic-scale pores are usually filled with fluid. One logical means of modelling the essential features of such rocks is to use poroelasticity theory. But previous models of wave propagation in cracked porous medium are either restricted to low frequencies at which effects of the elastic scattering (scattering into fast-P and S waves via mode conversion at the crack faces) are negligible or to the case that the crack-filling fluid is assumed to be incompressible. To overcome these restrictions, we consider the effects of crack fluid compressibility by extending spring condition into poroelasticity and derive exact solutions of the scattering problem of an incident P wave by a circular crack containing compressible fluid in a porous medium. Based on the solutions, we develop two different effective medium models to estimate frequency-dependent effective velocity and attenuation in a fluid-saturated porous rock with a set of aligned cracks. The mixed-boundary value problem reveals that both the wave-induced fluid flow (WIFF) and elastic wave scattering can cause important velocity dispersion and attenuation. The diffusion-type WIFF dominates the velocity change and attenuation for the low frequency range, while the elastic scattering dominates them for the relatively higher frequency range. The dependences of the P-wave velocity on the crack fluid compressibility are different at different frequencies. For the WIFF-dominated frequency range and Rayleigh-scattering frequency range, the P-wave velocity decreases with the crack fluid compressibility. In contrast, for the Mie scattering frequency range, the opposite occurs (the P-wave velocity increases with the crack fluid compressibility).


Sign in / Sign up

Export Citation Format

Share Document