Evaluation of 3C microelectromechanical system data on a 2D line: Direct comparison with conventional vertical-component geophone arrays and PS-wave analysis

Geophysics ◽  
2011 ◽  
Vol 76 (3) ◽  
pp. B79-B87 ◽  
Author(s):  
Christian Stotter ◽  
Erika Angerer

A 2D vibroseis line was acquired in the Vienna Basin (Austria) for the purpose of comparing the data of digital multicomponent single sensors based on microelectromechanical system (MEMS) sensors alongside conventional vertical-component geophone arrays. For efficient removal of coherent noise during processing, all source points were recorded in single-sweep mode, i.e., no vertical stacking was performed in the field. On this densely sampled data set, several noise-reduction techniques, such as digital array forming, frequency-wavenumber (f-k) filtering in shot and receiver domains, and polarization filters, proved to be valuable in reducing source-generated noise. The results showed that, with the use of single-sweep recording and polarization filter techniques, it is possible to produce seismic sections for a single-receiver three-component (3C) MEMS line that are comparable to a conventional geophone array line in signal-to-noise ratio. However, the higher number of single geophones and hence the stronger attenuation of random noise in the conventional array resulted in an advantage for the analog geophone data set. The second goal for this survey was to evaluate additional information contained in the horizontal components of the MEMS data. The multicomponent data allowed for the processing of mode-converted shear-wave data, performed for the first time in the Vienna Basin. Azimuthal anisotropy related to horizontal stresses was observed in the Neogene section of the shear-wave data set. A PP-PS event correlation allowed the identification of major shallow horizons. Interpretation of the final sections confirmed that the PS data are useful to distinguish between gas reservoirs and high-porosity water sands, which can cause similar P-wave amplitude variation with offset (AVO) effects.

Geophysics ◽  
1999 ◽  
Vol 64 (4) ◽  
pp. 1126-1138 ◽  
Author(s):  
Leon Thomsen ◽  
Ilya Tsvankin ◽  
Michael C. Mueller

Alford rotation analysis of 2C × 2C shear‐wave data (two source components, two receiver components) for azimuthal anisotropy is valid only when the orientation of that azimuthal anisotropy is invariant with depth. The Winterstein and Meadows method of layer stripping vertical seismic profiling (VSP) data relaxes this restriction for coarse‐layer variation of the orientation of the anisotropy. Here we present a tensor generalization of the conventional convolutional model of scalar wave propagation and use it to derive generalizations of Winterstein and Meadows layer stripping, valid for 2C × 2C data and for the restricted 2C-only case, in the VSP and reflection contexts. In the 2C × 2C VSP application, the result reduces to that of Winterstein and Meadows in the case where both fast and slow shear modes have the same attenuation and dispersion; otherwise, a balancing of mode spectra and amplitudes is required. The 2C × 2C reflection result differs from the 2C × 2C VSP result since two applications of the mode‐balancing and mode‐advance operations are required (since the waves travel up as well as down). Application to a synthetic data set confirms these results. The 2C × 2C reflection algorithm enables the exploration for sweet spots of high fracture intensity ahead of the bit without the restrictive assumption that the anisotropy orientation is depth invariant.


2022 ◽  
Vol 41 (1) ◽  
pp. 47-53
Author(s):  
Zhiwen Deng ◽  
Rui Zhang ◽  
Liang Gou ◽  
Shaohua Zhang ◽  
Yuanyuan Yue ◽  
...  

The formation containing shallow gas clouds poses a major challenge for conventional P-wave seismic surveys in the Sanhu area, Qaidam Basin, west China, as it dramatically attenuates seismic P-waves, resulting in high uncertainty in the subsurface structure and complexity in reservoir characterization. To address this issue, we proposed a workflow of direct shear-wave seismic (S-S) surveys. This is because the shear wave is not significantly affected by the pore fluid. Our workflow includes acquisition, processing, and interpretation in calibration with conventional P-wave seismic data to obtain improved subsurface structure images and reservoir characterization. To procure a good S-wave seismic image, several key techniques were applied: (1) a newly developed S-wave vibrator, one of the most powerful such vibrators in the world, was used to send a strong S-wave into the subsurface; (2) the acquired 9C S-S data sets initially were rotated into SH-SH and SV-SV components and subsequently were rotated into fast and slow S-wave components; and (3) a surface-wave inversion technique was applied to obtain the near-surface shear-wave velocity, used for static correction. As expected, the S-wave data were not affected by the gas clouds. This allowed us to map the subsurface structures with stronger confidence than with the P-wave data. Such S-wave data materialize into similar frequency spectra as P-wave data with a better signal-to-noise ratio. Seismic attributes were also applied to the S-wave data sets. This resulted in clearly visible geologic features that were invisible in the P-wave data.


Geophysics ◽  
2000 ◽  
Vol 65 (5) ◽  
pp. 1446-1454 ◽  
Author(s):  
Side Jin ◽  
G. Cambois ◽  
C. Vuillermoz

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.


Geophysics ◽  
2016 ◽  
Vol 81 (6) ◽  
pp. P57-P70 ◽  
Author(s):  
Shaun Strong ◽  
Steve Hearn

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.


Geophysics ◽  
2016 ◽  
Vol 81 (4) ◽  
pp. ID59-ID71 ◽  
Author(s):  
Kyle Basler-Reeder ◽  
John Louie ◽  
Satish Pullammanappallil ◽  
Graham Kent

Joint seismic and gravity analyses of the San Emidio geothermal field in the northwest Basin and Range province of Nevada demonstrate that joint optimization changes interpretation outcomes. The prior 0.3–0.5 km deep basin interpretation gives way to a deeper than 1.3 km basin model. Kirchoff prestack depth migrations reveal that joint optimization ameliorates shallow velocity artifacts, flattening antiformal reflectors that could have been interpreted as folds. Furthermore, joint optimization provides a clearer picture of the rangefront fault by increasing the depth of constrained velocities, which improves reflector coherency at depth. This technique provides new insight when applied to existing data sets and could replace the existing strategy of forward modeling to match gravity data. We have achieved stable joint optimization through simulated annealing, a global optimization algorithm that does not require an accurate initial model. Balancing the combined seismic-gravity objective function is accomplished by a new approach based on analysis of Pareto charts. Gravity modeling uses an efficient convolution model, and the basis of seismic modeling is the highly efficient Vidale eikonal equation traveltime generation technique. Synthetic tests found that joint optimization improves velocity model accuracy and provides velocity control below the deepest headwave raypath. Restricted offset-range migration analysis provides insights into precritical and gradient reflections in the data set.


Geophysics ◽  
1999 ◽  
Vol 64 (4) ◽  
pp. 1312-1328 ◽  
Author(s):  
Heloise B. Lynn ◽  
Wallace E. Beckham ◽  
K. Michele Simon ◽  
C. Richard Bates ◽  
M. Layman ◽  
...  

Reflection P- and S-wave data were used in an investigation to determine the relative merits and strengths of these two data sets to characterize a naturally fractured gas reservoir in the Tertiary Upper Green River formation. The objective is to evaluate the viability of P-wave seismic to detect the presence of gas‐filled fractures, estimate fracture density and orientation, and compare the results with estimates obtained from the S-wave data. The P-wave response to vertical fractures must be evaluated at different source‐receiver azimuths (travelpaths) relative to fracture strike. Two perpendicular lines of multicomponent reflection data were acquired approximately parallel and normal to the dominant strike of Upper Green River fractures as obtained from outcrop, core analysis, and borehole image logs. The P-wave amplitude response is extracted from prestack amplitude variation with offset (AVO) analysis, which is compared to isotropic‐model AVO responses of gas sand versus brine sand in the Upper Green River. A nine‐component vertical seismic profile (VSP) was also obtained for calibration of S-wave reflections with P-wave reflections, and support of reflection S-wave results. The direction of the fast (S1) shear‐wave component from the reflection data and the VSP coincides with the northwest orientation of Upper Green River fractures, and the direction of maximum horizontal in‐situ stress as determined from borehole ellipticity logs. Significant differences were observed in the P-wave AVO gradient measured parallel and perpendicular to the orientation of Upper Green River fractures. Positive AVO gradients were associated with gas‐producing fractured intervals for propagation normal to fractures. AVO gradients measured normal to fractures at known waterwet zones were near zero or negative. A proportional relationship was observed between the azimuthal variation of the P-wave AVO gradient as measured at the tops of fractured intervals, and the fractional difference between the vertical traveltimes of split S-waves (the “S-wave anisotropy”) of the intervals.


Geophysics ◽  
1999 ◽  
Vol 64 (4) ◽  
pp. 1172-1180 ◽  
Author(s):  
W. Scott Leaney ◽  
Colin M. Sayers ◽  
Douglas E. Miller

Multioffset vertical seismic profile (VSP) experiments, commonly referred to as walkaways, enable anisotropy to be measured reliably in the field. The results can be fed into modeling programs to study the impact of anisotropy on velocity analysis, migration, and amplitude versus offset (AVO). Properly designed multioffset VSPs can also provide the target AVO response measured under optimum conditions, since the wavelet is recorded just above the reflectors of interest with minimal reflection point dispersal. In this paper, the multioffset VSP technique is extended to include multioffset azimuths, and a multiazimuthal multiple VSP data set acquired over a carbonate reservoir is analyzed for P-wave anisotropy and AVO. Direct arrival times down to the overlying shale and reflection times and amplitudes from the carbonate are analyzed. Data analysis involves a three‐term fit to account for nonhyperbolic moveout, dip, and azimuthal anisotropy. Results indicate that the overlying shale is transversely isotropic with a vertical axis of symmetry (VTI), while the carbonate shows 4–5% azimuthal anisotropy in traveltimes. The fast direction is consistent with the maximum horizontal stress orientation determined from break‐out logs and is also consistent with the strike of major faults. AVO analysis of the reflection from the top of the carbonate layer shows a critical angle reduction in the fast direction and maximum gradient in the slow direction. This agrees with modeling and indicates a greater amplitude sensitivity in the slow direction—the direction perpendicular to fracture strike. In principle, 3-D surveys should have wide azimuthal coverage to characterize fractured reservoirs. If this is not possible, it is important to have azimuthal line coverage in the minimum horizontal stress direction to optimize the use of AVO for fractured reservoir characterization. This direction can be obtained from multiazimuthal walkaways using the azimuthal P-wave analysis techniques presented.


Geophysics ◽  
1998 ◽  
Vol 63 (2) ◽  
pp. 692-706 ◽  
Author(s):  
Subhashis Mallick ◽  
Kenneth L. Craft ◽  
Laurent J. Meister ◽  
Ronald E. Chambers

In an azimuthally anisotropic medium, the principal directions of azimuthal anisotropy are the directions along which the quasi-P- and the quasi-S-waves propagate as pure P and S modes. When azimuthal anisotropy is induced by oriented vertical fractures imposed on an azimuthally isotropic background, two of these principal directions correspond to the directions parallel and perpendicular to the fractures. S-waves propagating through an azimuthally anisotropic medium are sensitive to the direction of their propagation with respect to the principal directions. As a result, primary or mode‐converted multicomponent S-wave data are used to obtain the principal directions. Apart from high acquisition cost, processing and interpretation of multicomponent data require a technology that the seismic industry has not fully developed. Anisotropy detection from conventional P-wave data, on the other hand, has been limited to a few qualitative studies of the amplitude variation with offset (AVO) for different azimuthal directions. To quantify the azimuthal AVO, we studied the amplitude variation with azimuth for P-wave data at fixed offsets. Our results show that such amplitude variation with azimuth is periodic in 2θ, θ being the orientation of the shooting direction with respect to one of the principal directions. For fracture‐induced anisotropy, this principal direction corresponds to the direction parallel or perpendicular to the fractures. We use this periodic azimuthal dependence of P-wave reflection amplitudes to identify two distinct cases of anisotropy detection. The first case is an exactly determined one, where we have observations from three azimuthal lines for every common‐midpoint (CMP) location. We derive equations to compute the orientation of the principal directions for such a case. The second case is an overdetermined one where we have observations from more than three azimuthal lines. Orientation of the principal direction from such an overdetermined case can be obtained from a least‐squares fit to the reflection amplitudes over all the azimuthal directions or by solving many exactly determined problems. In addition to the orientation angle, a qualitative measure of the degree of azimuthal anisotropy can also be obtained from either of the above two cases. When azimuthal anisotropy is induced by oriented vertical fractures, this qualitative measure of anisotropy is proportional to fracture density. Using synthetic seismograms, we demonstrate the robustness of our method in evaluating the principal directions from conventional P-wave seismic data. We also apply our technique to real P-wave data, collected over a wide source‐to‐receiver azimuth distribution. Computations using our method gave an orientation of the principal direction consistent with the general fracture orientation in the area as inferred from other geological and geophysical evidence.


Geophysics ◽  
2000 ◽  
Vol 65 (5) ◽  
pp. 1613-1621 ◽  
Author(s):  
Richard L. Gibson ◽  
Stephen Theophanis ◽  
M. Nafi Toksöz

Fractured reservoirs are an important target for exploration and production geophysics, and the azimuthal anisotropy often associated with these reservoirs can strongly influence seismic wave propagation. We created a physical model of a fractured reservoir to simulate some of these propagation effects. The reservoir is represented by a phenolite disk that is thin with respect to the elastic wavelengths in the experiment, creating model dimensions that are representative of realistic reservoirs. Phenolite is strongly anisotropic with orthorhombic symmetry, which suggests that azimuthal amplitude versus offset (AVO) effects should be obvious in data. We acquired both SH- and P-wave data in common‐offset gathers with a near offset and a far offset and found that although the SH-wave data show clear azimuthal variations in AVO, the P-wave signals show no apparent changes with azimuth. We then applied numerical modeling to analyze the data. Because ray methods cannot model diffractions from the disk edge, we first used a ray‐Born technique to simulate variations in waveforms associated with such scattering. The synthetic seismograms reproduced variations in the SH-wave waveforms accurately, though the amplitude contrast between acquisition azimuths was overestimated. Assuming a laterally homogeneous model, we then applied ray methods to simulate tuning effects in SH- and P-wave data and confirmed that in spite of the large contrasts in elastic properties, the tuning of the P-wave reflections from the thin disk changed so there was negligible contrast in AVO with azimuth. Models of field scale reservoirs showed that the same effects could be expected for field applications.


Sign in / Sign up

Export Citation Format

Share Document