scholarly journals Fracture compliance estimation using borehole tube waves

Geophysics ◽  
2013 ◽  
Vol 78 (4) ◽  
pp. D249-D260 ◽  
Author(s):  
Sudhish Kumar Bakku ◽  
Michael Fehler ◽  
Daniel Burns
Keyword(s):  
Amplitude Ratio ◽  
Wave Generation ◽  
Transition Frequency ◽  
P Wave ◽  
Fracture Aperture ◽  
Generation Model ◽  
Data Set ◽  
Fracture Properties ◽  
Lateral Extent ◽  
Tube Wave

We tested two models, one for tube-wave generation and the other for tube-wave attenuation at a fracture intersecting a borehole that can be used to estimate fracture compliance, fracture aperture, and lateral extent. In the tube-wave generation model, we consider tube-wave excitation in the borehole when a P-wave is incident on the fracture. The amplitude ratio of the pressure due to the tube wave to that of the incident P-wave is a function of fracture compliance, aperture, and length. Similarly, the attenuation of a tube wave in the borehole as it crosses a fracture intersecting the borehole is also a function of fracture properties. Numerically solving the dispersion relation in the fracture, we study tube-wave generation and the attenuation coefficient as a function of frequency. We observed that measuring amplitude ratios or attenuation near a transition frequency can help constrain the fracture properties. The transition frequency corresponds to the regime in which the viscous skin depth in the fracture is comparable to its aperture. Measurements in the high-frequency limit can place a lower bound on fracture compliance and lateral extent. We evaluated the applicability of the tube-wave generation model to a previously published VSP data set and found that compliance values of the order [Formula: see text]–[Formula: see text] are likely in the field. These observations support scaling of fracture compliance with fracture size.

The 1998 Valhall microseismic data set: An integrated study of relocated sources, seismic multiplets, and S-wave splitting

Geophysics ◽  
2009 ◽  
Vol 74 (5) ◽  
pp. B183-B195 ◽  
Author(s):  
K. De Meersman ◽  
J.-M. Kendall ◽  
M. van der Baan
Keyword(s):  
Seismic Anisotropy ◽  
Amplitude Ratio ◽  
Temporal Variations ◽  
Oil Field ◽  
P Wave ◽  
Wave Form ◽  
S Wave ◽  
Data Set ◽  
Wave Splitting ◽  
Source Mechanisms

We relocate 303 microseismic events recorded in 1998 by sensors in a single borehole in the North Sea Valhall oil field. A semiautomated array analysis method repicks the P- and S-wave arrival times and P-wave polarizations, which are needed to locate these events. The relocated sources are confined predominantly to a [Formula: see text]-thick zone just above the reservoir, and location uncertainties are half those of previous efforts. Multiplet analysis identifies 40 multiplet groups, which include 208 of the 303 events. The largest group contains 24 events, and five groups contain 10 or more events. Within each multiplet group, we further improve arrival-time picking through crosscorrelation, which enhances the relative accuracy of the relocated events and reveals that more than 99% of the seismic activity lies spatially in three distinct clusters. The spatial distribution of events and wave-form similarities reveal two faultlike structures that match well with north-northwest–south-southeast-trending fault planes interpreted from 3D surface seismic data. Most waveform differences between multiplet groups located on these faults can be attributed to S-wave phase content and polarity or P-to-S amplitude ratio. The range in P-to-S amplitude ratios observed on the faults is explained best in terms of varying source mechanisms. We also find a correlation between multiplet groups and temporal variations in seismic anisotropy, as revealed by S-wave splitting analysis. We explain these findings in the context of a cyclic recharge and dissipation of cap-rock stresses in response to production-driven compaction of the underlying oil reservoir. The cyclic nature of this mechanism drives the short-term variations in seismic anisotropy and the reactivation of microseismic source mechanisms over time.

scholarly journals Tube‐wave suppression in single‐well seismic acquisition

Geophysics ◽  
2003 ◽  
Vol 68 (3) ◽  
pp. 863-869 ◽  
Author(s):  
Thomas M. Daley ◽  
Roland Gritto ◽  
Ernest L. Majer ◽  
Phillip West
Keyword(s):  
Wave Amplitude ◽  
Amplitude Ratio ◽  
Fluid Pressure ◽  
Oil Field ◽  
Body Waves ◽  
P Wave ◽  
Reservoir Rock ◽  
Wave Tube ◽  
Single Well ◽  
Tube Wave

Single‐well seismic imaging is significantly hampered by the presence of borehole tube waves. A tube‐wave suppressor has been tested using single‐well seismic equipment at the Lost Hills (California) oil field. The suppressor uses a gas‐filled bladder kept slightly above borehole fluid pressure. Field tests show a measurable reduction in tube‐wave energy as compared to body waves propagating in the surrounding reservoir rock. When using a high‐frequency (500–4000 Hz) piezoelectric source, the P‐wave–tube‐wave amplitude ratio was increased by 33 dB. When using a lower frequency (50–350 Hz) orbital vibrator source, the S‐wave–tube‐wave amplitude ratio was increased by 21 dB while the P‐wave–tube‐wave amplitude ratio was increased by 23 dB. These reductions in tube‐wave amplitudes significantly improve single‐well data quality.

P-wave backazimuth anomalies observed by a small-aperture seismic array at Pinyon Flat, southern California: Implications for structure and source location

1996 ◽  
Vol 86 (2) ◽  
pp. 470-476 ◽  
Author(s):  
Cheng-Horng Lin ◽  
S. W. Roecker
Keyword(s):  
Southern California ◽  
Seismic Array ◽  
P Wave ◽  
Small Aperture ◽  
Data Set ◽  
Single Interface ◽  
Seismic Networks ◽  
Dense Seismic Array ◽  
Beamforming Technique

Abstract Seismograms of earthquakes and explosions recorded at local, regional, and teleseismic distances by a small-aperture, dense seismic array located on Pinyon Flat, in southern California, reveal large (±15°) backazimuth anomalies. We investigate the causes and implications of these anomalies by first comparing the effectiveness of estimating backazimuth with an array using three different techniques: the broadband frequency-wavenumber (BBFK) technique, the polarization technique, and the beamforming technique. While each technique provided nearly the same direction as a most likely estimate, the beamforming estimate was associated with the smallest uncertainties. Backazimuth anomalies were then calculated for the entire data set by comparing the results from beamforming with backazimuths derived from earthquake locations reported by the Anza and Caltech seismic networks and the Preliminary Determination of Epicenters (PDE) Bulletin. These backazimuth anomalies have a simple sinelike dependence on azimuth, with the largest anomalies observed from the southeast and northwest directions. Such a trend may be explained as the effect of one or more interfaces dipping to the northeast beneath the array. A best-fit model of a single interface has a dip and strike of 20° and 315°, respectively, and a velocity contrast of 0.82 km/sec. Application of corrections computed from this simple model to ray directions significantly improves locations at all distances and directions, suggesting that this is an upper crustal feature. We confirm that knowledge of local structure can be very important for earthquake location by an array but also show that corrections computed from simple models may not only be adequate but superior to those determined by raytracing through smoothed laterally varying models.

Source dynamics of the 1979 Imperial Valley earthquake from near-source observations (of ground acceleration and velocity)

1982 ◽  
Vol 72 (6A) ◽  
pp. 1957-1968
Author(s):  
Mansour Niazi
Keyword(s):  
Particle Acceleration ◽  
Particle Velocity ◽  
High Frequency ◽  
Observation Point ◽  
P Wave ◽  
Imperial Valley ◽  
Ground Acceleration ◽  
Data Set ◽  
Ground Shaking ◽  
Fault Surface

abstract Two sets of observations obtained during the 15 October 1979 Imperial Valley earthquake, MS 6.9, are presented. The data suggest different dynamic characteristics of the source when viewed in different frequency bands. The first data set consists of the observed residuals of the horizontal peak ground accelerations and particle velocity from predicted values within 50 km of the fault surface. The residuals are calculated from a nonlinear regression analysis of the data (Campbell, 1981) to the following empirical relationships, PGA = A 1 ( R + C 1 ) − d 1 , PGV = A 2 ( R + C 2 ) − d 2 in which R is the closest distance to the plane of rupture. The so-calculated residuals are correlated with a positive scalar factor signifying the focusing potential at each observation point. The focusing potential is determined on the basis of the geometrical relation of the station relative to the rupture front on the fault plane. The second data set consists of the acceleration directions derived from the windowed-time histories of the horizontal ground acceleration across the El Centro Differential Array (ECDA). The horizontal peak velocity residuals and the low-pass particle acceleration directions across ECDA require the fault rupture to propagate northwestward. The horizontal peak ground acceleration residuals and the high-frequency particle acceleration directions, however, are either inconclusive or suggest an opposite direction for rupture propagation. The inconsistency can best be explained to have resulted from the incoherence of the high-frequency radiation which contributes most effectively to the registration of PGA. A test for the sensitivity of the correlation procedure to the souce location is conducted by ascribing the observed strong ground shaking to a single asperity located 12 km northwest of the hypocenter. The resulting inconsistency between the peak acceleration and velocity observations in relation to the focusing potential is accentuated. The particle velocity of Delta Station, Mexico, in either case appears abnormally high and disagrees with other observations near the southeastern end of the fault trace. From the observation of a nearly continuous counterclockwise rotation of the plane of P-wave particle motion at ECDA, the average rupture velocity during the first several seconds of source activation is estimated to be 2.0 to 3.0 km/sec. A 3 km upper bound estimate of barrier dimensions is tentatively made on the basis of the observed quasiperiodic variation of the polarization angles.

Shear‐wave velocity and density estimation from PS-wave AVO analysis: Application to an OBS dataset from the North Sea

Geophysics ◽  
2000 ◽  
Vol 65 (5) ◽  
pp. 1446-1454 ◽  
Author(s):  
Side Jin ◽  
G. Cambois ◽  
C. Vuillermoz
Keyword(s):  
Wave Velocity ◽  
North Sea ◽  
Random Noise ◽  
P Wave ◽  
S Wave ◽  
Data Set ◽  
Avo Analysis ◽  
The North ◽  
Avo Inversion ◽  
The North Sea

S-wave velocity and density information is crucial for hydrocarbon detection, because they help in the discrimination of pore filling fluids. Unfortunately, these two parameters cannot be accurately resolved from conventional P-wave marine data. Recent developments in ocean‐bottom seismic (OBS) technology make it possible to acquire high quality S-wave data in marine environments. The use of (S)-waves for amplitude variation with offset (AVO) analysis can give better estimates of S-wave velocity and density contrasts. Like P-wave AVO, S-wave AVO is sensitive to various types of noise. We investigate numerically and analytically the sensitivity of AVO inversion to random noise and errors in angles of incidence. Synthetic examples show that random noise and angle errors can strongly bias the parameter estimation. The use of singular value decomposition offers a simple stabilization scheme to solve for the elastic parameters. The AVO inversion is applied to an OBS data set from the North Sea. Special prestack processing techniques are required for the success of S-wave AVO inversion. The derived S-wave velocity and density contrasts help in detecting the fluid contacts and delineating the extent of the reservoir sand.

scholarly journals Crustal velocity structure of the Apennines (Italy) from P-wave travel time tomography

1996 ◽  
Vol 39 (6) ◽  
Author(s):  
C. Chiarabba ◽  
A. Amato
Keyword(s):  
Velocity Structure ◽  
P Wave ◽  
Depth Interval ◽  
Low Velocity ◽  
Data Set ◽  
Arrival Times ◽  
Velocity Anomaly ◽  
Minimum Number ◽  
Velocity Images ◽  
Lower Crustal

In this paper we provide P-wave velocity images of the crust underneath the Apennines (Italy), focusing on the lower crustal structure and the Moho topography. We inverted P-wave arrival times of earthquakes which occurred from 1986 to 1993 within the Apenninic area. To overcome inversion instabilities due to noisy data (we used bulletin data) we decided to resolve a minimum number of velocity parameters, inverting for only two layers in the crust and one in the uppermost mantle underneath the Moho. A partial inversion of only 55% of the overall dataset yields velocity images similar to those obtained with the whole data set, indicating that the depicted tomograms are stable and fairly insensitive to the number of data used. We find a low-velocity anomaly in the lower crust extending underneath the whole Apenninic belt. This feature is segmented by a relative high-velocity zone in correspondence with the Ortona-Roccamonfina line, that separates the northern from the southern Apenninic arcs. The Moho has a variable depth in the study area, and is deeper (more than 37 km) in the Adriatic side of the Northern Apennines with respect to the Tyrrhenian side, where it is found in the depth interval 22-34 km.

scholarly journals Upscaling of elastic properties in carbonates: A modeling approach based on a multiscale geophysical data set

2020 ◽  
Author(s):  
Jerome Fortin ◽  
Cedric Bailly ◽  
Mathilde Adelinet ◽  
Youri Hamon
Keyword(s):  
Elastic Properties ◽  
Wave Velocity ◽  
Representative Elementary Volume ◽  
P Wave ◽  
Elementary Volume ◽  
Data Set ◽  
Wave Velocities ◽  
P Wave Velocity ◽  
Ultrasonic Measurements ◽  
Seismic Scale

<p>Linking ultrasonic measurements made on samples, with sonic logs and seismic subsurface data, is a key challenge for the understanding of carbonate reservoirs. To deal with this problem, we investigate the elastic properties of dry lacustrine carbonates. At one study site, we perform a seismic refraction survey (100 Hz), as well as sonic (54 kHz) and ultrasonic (250 kHz) measurements directly on outcrop and ultrasonic measurements on samples (500 kHz). By comparing the median of each data set, we show that the P wave velocity decreases from laboratory to seismic scale. Nevertheless, the median of the sonic measurements acquired on outcrop surfaces seems to fit with the seismic data, meaning that sonic acquisition may be representative of seismic scale. To explain the variations due to upscaling, we relate the concept of representative elementary volume with the wavelength of each scale of study. Indeed, with upscaling, the wavelength varies from millimetric to pluri-metric. This change of scale allows us to conclude that the behavior of P wave velocity is due to different geological features (matrix porosity, cracks, and fractures) related to the different wavelengths used. Based on effective medium theory, we quantify the pore aspect ratio at sample scale and the crack/fracture density at outcrop and seismic scales using a multiscale representative elementary volume concept. Results show that the matrix porosity that controls the ultrasonic P wave velocities is progressively lost with upscaling, implying that crack and fracture porosity impacts sonic and seismic P wave velocities, a result of paramount importance for seismic interpretation based on deterministic approaches.</p><p>Bailly, C., Fortin, J., Adelinet, M., & Hamon, Y. (2019). Upscaling of elastic properties in carbonates: A modeling approach based on a multiscale geophysical data set. Journal of Geophysical Research: Solid Earth, 124. https://doi.org/10.1029/2019JB018391</p>

Survey design for coal-scale 3D-PS seismic reflection

Geophysics ◽  
2016 ◽  
Vol 81 (6) ◽  
pp. P57-P70 ◽  
Author(s):  
Shaun Strong ◽  
Steve Hearn
Keyword(s):  
Survey Design ◽  
High Amplitude ◽  
Azimuthal Anisotropy ◽  
Mine Planning ◽  
P Wave ◽  
Depth Range ◽  
Low Velocity ◽  
Converted Wave ◽  
Data Set ◽  
Economic Considerations

Survey design for converted-wave (PS) reflection is more complicated than for standard P-wave surveys, due to raypath asymmetry and increased possibility of phase distortion. Coal-scale PS surveys (depth [Formula: see text]) require particular consideration, partly due to the particular physical properties of the target (low density and low velocity). Finite-difference modeling provides a pragmatic evaluation of the likely distortion due to inclusion of postcritical reflections. If the offset range is carefully chosen, then it may be possible to incorporate high-amplitude postcritical reflections without seriously degrading the resolution in the stack. Offsets of up to three times target depth may in some cases be usable, with appropriate quality control at the data-processing stage. This means that the PS survey design may need to handle raypaths that are highly asymmetrical and that are very sensitive to assumed velocities. A 3D-PS design was used for a particular coal survey with the target in the depth range of 85–140 m. The objectives were acceptable fold balance between bins and relatively smooth distribution of offset and azimuth within bins. These parameters are relatively robust for the P-wave design, but much more sensitive for the case of PS. Reduction of the source density is more acceptable than reduction of the receiver density, particularly in terms of the offset-azimuth distribution. This is a fortuitous observation in that it improves the economics of a dynamite source, which is desirable for high-resolution coal-mine planning. The final-survey design necessarily allows for logistical and economic considerations, which implies some technical compromise. However, good fold, offset, and azimuth distributions are achieved across the survey area, yielding a data set suitable for meaningful analysis of P and S azimuthal anisotropy.

Inversion of differences in frequency components of azimuthal seismic data for indicators of oil-bearing fractured reservoirs based on an attenuative cracked model

Geophysics ◽  
2020 ◽  
Vol 85 (3) ◽  
pp. R163-R175
Author(s):  
Huaizhen Chen ◽  
Junxiao Li ◽  
Kristopher A. Innanen
Keyword(s):  
Reflection Coefficient ◽  
Stiffness Matrix ◽  
Fractured Reservoirs ◽  
P Wave ◽  
Parameter Estimates ◽  
Fluid Viscosity ◽  
Mathematical Form ◽  
Data Set ◽  
Frequency Dependent ◽  
Frequency Components

Based on a model of attenuative cracked rock, we have derived a simplified and frequency-dependent stiffness matrix associated with (1) a rock volume containing aligned and partially saturated cracks and (2) a new indicator of oil-bearing fractured reservoirs, which is related to pressure relaxation in cracked rocks and influenced by fluid viscosity and saturation. Starting from the mathematical form of a perturbation in this stiffness matrix across a reflecting interface separating two attenuative cracked media, we set up a linearized P-wave to P-wave reflection coefficient as an azimuthally and frequency-dependent function of dry rock elastic properties, dry fracture weaknesses, and the new indicator. By varying this reflection coefficient with azimuthal angle, we derive a further expression referred to as the quasidifference in elastic impedance, or [Formula: see text], which is primarily affected by the dry fracture weaknesses and the new indicator. An inversion approach is established to use differences in frequency components of seismic amplitudes to estimate these weaknesses and the indicator based on the derived [Formula: see text]. In synthetic inversion tests, we determine that the approach produces interpretable parameter estimates in the presence of data with a moderate signal-to-noise ratio (S/N). Testing on a real data set suggests that reliable fracture weakness and indicator are generated by the approach; fractured and oil-bearing reservoirs are identified through a combination of the dry fracture weakness and the new indicator.

scholarly journals Potential shallow aquifers characterization through an integrated geophysical method: multivariate approach by means of k-means algorithms

2017 ◽  
Vol 6 (2) ◽  
Author(s):  
Stefano Bernardinetti ◽  
Stefano Maraio ◽  
Pier Paolo Gennaro Bruno ◽  
Valentina Cicala ◽  
Serena Minucci ◽  
...  
Keyword(s):  
Groundwater Resources ◽  
Geophysical Methods ◽  
Optimal Number ◽  
P Wave ◽  
Quaternary Period ◽  
Data Set ◽  
Alluvial Deposits ◽  
Wave Refraction ◽  
Bivariate Data ◽  
Shallow Basin

The need to obtain a detailed hydrogeological characterization of the subsurface and its interpretation for the groundwater resources management, often requires to apply several and complementary geophysical methods. The goal of the approach in this paper is to provide a unique model of the aquifer by synthesizing and optimizing the information provided by several geophysical methods. This approach greatly reduces the degree of uncertainty and subjectivity of the interpretation by exploiting the different physical and mechanic characteristics of the aquifer. The studied area, into the municipality of Laterina (Arezzo, Italy), is a shallow basin filled by lacustrine and alluvial deposits (Pleistocene and Olocene epochs, Quaternary period), with alternated silt, sand with variable content of gravel and clay where the bottom is represented by arenaceous-pelitic rocks (Mt. Cervarola Unit, Tuscan Domain, Miocene epoch). This shallow basin constitutes the unconfined superficial aquifer to be exploited in the nearly future. To improve the geological model obtained from a detailed geological survey we performed electrical resistivity and P wave refraction tomographies along the same line in order to obtain different, independent and integrable data sets. For the seismic data also the reflected events have been processed, a remarkable contribution to draw the geologic setting. Through the k-means algorithm, we perform a cluster analysis for the bivariate data set to individuate relationships between the two sets of variables. This algorithm allows to individuate clusters with the aim of minimizing the dissimilarity within each cluster and maximizing it among different clusters of the bivariate data set. The optimal number of clusters “K”, corresponding to the individuated geophysical facies, depends to the multivariate data set distribution and in this work is estimated with the Silhouettes. The result is an integrated tomography that shows a finite number of homogeneous geophysical facies, which therefore permits to distinguish and interpret the porous aquifer in a quantitative and objective way.

Sign in / Sign up

Export Citation Format

Share Document