VpVs from mode‐converted P-SV reflections

Geophysics ◽  
1989 ◽  
Vol 54 (7) ◽  
pp. 843-852 ◽  
Author(s):  
William P. Iverson ◽  
Bill A. Fahmy ◽  
Scott B. Smithson

P-SV reflections are generated by a compressional‐wave source and result from P waves that are converted to shear (SV) waves upon reflection. Recording both the P and SV components yields compressional and shear data simultaneously. Verifying that the easily detected events really are P-SV reflections is accomplished by noting the good correlation of surface CDP data with vertical seismic profile (VSP) reflections. Stacking velocities from P-SV CDP gathers determine the [Formula: see text] product when source‐to‐receiver offset is less than the depth of the reflector but data from synthetic models show that P-SV reflections are nonhyperbolic for shallow reflections or when source‐to‐receiver offset is too large. Shear velocity [Formula: see text] can be calculated from P-SV reflections by one of two techniques: comparison of stacked section P-P and P-SV reflection times or by using the P-P and P-SV stacking velocities. Unfortunately, most P-SV reflections on a P-SV seismic section do not necessarily originate from exactly the same depth as P-P reflections. When this depth discrepancy occurs, the reflection‐time comparison technique fails. In addition, [Formula: see text] cannot be calculated from P-SV reflections, and we must settle for the [Formula: see text] product from P-SV reflection stacking velocities. When P-SV stacking velocities are input to the familiar Dix equation, the resulting interval velocities yield the [Formula: see text] product.

Geophysics ◽  
1995 ◽  
Vol 60 (4) ◽  
pp. 955-967 ◽  
Author(s):  
Hiroshi Amano

An analytical solution to separate P‐waves and S‐waves in vertical seismic profile (VSP) wavefields is derived using combinations of certain terms of the formal solution for forward VSP modeling. Some practical applications of this method to synthetic seismograms and field data are investigated and evaluated. Little wave distortion is recognized, and the weak wavefield masked by dominant wavetrains can be extracted with this method. The decomposed wavefield is expressed in the frequency‐depth (f-z) domain as a linear combination of up to the third‐order differential of traces, which is approximated by trace differences in the practical separation process. In general, five traces with single‐component data are required in this process, but the same process is implemented with only three traces in the acoustic case. Two‐trace extrapolation is applied to each edge of the data gather to enhance the accuracy of trace difference. Since the formulas are developed in the f-z domain, the influence of anelasticity can be taken into account, and the calculation is carried out fast enough with the benefit of the fast Fourier transform (FFT).


1985 ◽  
Vol 25 (1) ◽  
pp. 254
Author(s):  
T.J.C. Prudence ◽  
J. Flentri

The Kanpa 1A Vertical Seismic Profile (VSP) was conducted for Shell by Schlumberger and incorporated variable time and depth sampling, different source offsets and recording in cased and uncased hole. Processing was performed using Shell proprietary programs, with particular attention to:Editing and resampling of the data setSeparation of upgoing and downgoing waves using FK and median filtersComparison of gain recovery based on modelling the amplitude decay of direct arrivals with an averaging process (AGC)Sensitivity of the final VSP stack to blanking of residual tube wave noiseZero-phase whitening of the VSP traceStacked VSP traces for Kanpa 1A were compared with a zero-phase seismic section and synthetic seismogram at the well. The VSP/seismic match is good and, due to poor synthetic/seismic correlation, was the basis for the final seismic/well tie. Interpretation of deep VSP data enabled the estimation of formation boundaries below the total depth of the well.It is concluded that VSPs can be invaluable in establishing well ties where seismic is poor or when detailed correlation is required (e.g. stratigraphic traps). Reflectors "ahead of the bit" can be interpreted from VSPs based on assumed velocities and VSP/seismic tie, and the predicted thickness and seismic character of the target interval. A consistent field configuration is recommended for acquisition with attention to tube wave suppression and adequate spatial and temporal sampling. Previous processing experience is advantageous if quick and reliable VSP results are required for decisions while drilling.


Geophysics ◽  
2003 ◽  
Vol 68 (3) ◽  
pp. 837-848 ◽  
Author(s):  
Gislain B. Madiba ◽  
George A. McMechan

Simultaneous elastic impedance inversion is performed on the 2D North Viking Graben seismic data set used at the 1994 SEG workshop on amplitude variation with offset and inversion. P‐velocity (Vp), S‐velocity (Vs), density logs, and seismic data are input to the inversion. The inverted P‐impedance and S‐impedance sections are used to generate an approximate compressional‐to‐shear velocity ratio (Vp/Vs) section which, in turn, is used along with water‐filled porosity (Swv) derived from the logs from two wells, to generate fluid estimate sections. This is possible as the reservoir sands have fairly constant total porosity of approximately 28 ± 4%, so the hydrocarbon filled porosity is the total porosity minus the water‐filled porosity. To enhance the separation of lithologies and of fluid content, we map Vp/Vs into Swv using an empirical crossplot‐derived relation. This mapping expands the dynamic range of the low end of the Vp/Vs values. The different lithologies and fluids are generally well separated in the Vp/Vs–Swv domain. Potential hydrocarbon reservoirs (as calibrated by the well data) are identified throughout the seismic section and are consistent with the fluid content estimations obtained from alternative computations. The Vp/Vs–Swv plane still does not produce unique interpretation in many situations. However, the critical distinction, which is between hydrocarbon‐bearing sands and all other geologic/reservoir configurations, is defined. Swv ≤ 0.17 and Vp/Vs ≤ 1.8 are the criteria that delineate potential reservoirs in this area, with decreasing Swv indicating a higher gas/oil ratio, and decreasing Vp/Vs indicating a higher sand/shale ratio. As these criteria are locally calibrated, they appear to be valid locally; they should not be applied to other data sets, which may exhibit significantly different relationships. However, the overall procedure should be generally applicable.


2000 ◽  
Vol 37 (2-3) ◽  
pp. 291-306 ◽  
Author(s):  
J Martignole ◽  
A J Calvert ◽  
R Friedman ◽  
P Reynolds

Results of deep seismic reflection survey along a 375 km long transect of the Grenville Province in western Quebec are combined with a review of geological observations and published isotopic ages. The seismic profile offers a remarkably clear image of the crust-mantle boundary and a good definition of the various crustal blocks. Crust about 44 km thick beneath the Grenville Front zone thins abruptly to ca. 36 km southeastward, perhaps the result of extension on southeast-dipping surfaces extending to the Moho. Other zones of relatively thin crust, although less pronounced, occur where Proterozoic crust overlies Archean crust, and beneath the Morin anorthosite complex. The thickest crust is found at the extreme southeast of the transect, east of the Morin anorthosite. From northwest to southeast, three main crustal subdivisions are (1) deformed Archean rocks with southeast-dipping reflectors in the Grenville Front zone, (2) an Archean parautochthon with northwest-dipping reflectors extending to the lower crust, and (3) an overlying three-layer crust interpreted as accreted Proterozoic terranes. The boundary between (2) and (3) is a major, southeast-dipping, crustal-scale ramp (Baskatong ramp) interpreted to have accommodated strain during and after accretion. U-Pb and Pb-Pb ages on detrital zircons show that metasedimentary rocks of the allochthons (Mont-Laurier, Réservoir Cabonga, and Lac Dumoine terranes) range from Archean to as young as 1.21 Ga. A single zone with 1.4 Ga old Sm-Nd model ages appears to lack Archean components and may be considered as a fragment of juvenile Mesoproterozoic crust pinched in a shear zone (Renzy shear zone) that could be raised to the status of terrane (Renzy terrane). In the allochthons, U-Pb ages of metamorphic zircon and monazite cluster around 1.17 Ga (Mont-Laurier and Réservoir Cabonga terrane) and 1.07 Ga (Renzy and Lac Dumoine terrane) and are interpreted to record late and post-accretion crustal reworking, a common feature of the Grenville orogen. A final high-grade metamorphic event (ca. 1.0 Ga) documented only in the parautochthon and the Grenville Front zone records large-scale, piggyback-style thrusting of allochthonous slabs onto the parautochthon. The age of transcurrent displacement following peak metamorphism affecting both the allochthons and the parautochthon decreases northwestward from 1.07 to 1.00 Ga. Dating thus shows that Grenvillian deformation in western Quebec occurred in pulses over an interval of 180 million years, with a tendency to propagate from the inner part of the orogen toward the Grenville Front. Reworked migmatites from the parautochthon cooled from the ca. 1.0 Ga peak of metamorphism through about 450°C (Ar closure in hornblende) at ca. 0.96 Ga with calculated cooling rates of about 6°C per million years, and unroofing rates of 0.33 km per million years. The cooling-unroofing history of the allochthons is not so straightforward, probably due to tectonic disturbances related to allochthon emplacement. Cooling through 450°C occurred between 1.04 and 1.01 Ga, at least 50 million years earlier than cooling in the parautochthon; this contrast agrees with the northwestward propagation of the orogen.


1980 ◽  
Vol 20 (05) ◽  
pp. 317-326 ◽  
Author(s):  
E.A. Koerperich

Abstract Acoustic waveforms from long- and short-spacedsonic logs were investigated to determine ifshort-spaced tools give accurate measurements of shear wave velocity. Compressional wave interference canaffect shear velocities from both tools adversely.However, the short-spaced tool was useful over awider range of conditions. Introduction The areas where shear velocity data can be appliedtheoretically or empirically are diverse. Most of theseinvolve use of the dynamic elastic rock constants, which can be computed from shear (S) velocity[along with compressional (P) velocity and bulkdensity, which are obtained readily from existingwireline logging devices]. Some of these applicationareas are (1) seismic amplitude calibration andinterpretation, (2) sand control,(3) formationfracturing, reservoir material balance and subsidencestudies(through relationships between rock andpore-volume changes with stress),(4) lithologyand porosity, 14 and (5) geopressure prediction. While rich in possible application areas, shearvelocity is difficult to measure automatically withconventional acoustic devices and detection schemes.Except in limited lithology-logging conditions, manual examination of waveforms commonly isrequired to extract shear velocity. Even then there has been considerable uncertainty in shear arrivals onshort-spaced tools due to P-wave interference. Insofter rocks, conventional tools simply do nottransmit distinct shear arrivals. Current axial transmitter-receiver (T-R) toolsare designed primarily for detection of P waves.Downhole amplifiers adjusted to accentuate the firstP-wave arrival normally saturate through the shearand late P regions of the waveform. When downholegain is reduced to eliminate amplifier saturation, initial shear arrivals generally are superimposed onlate P arrivals. This interference makes automaticdetection difficult and leads to a concern about theconsistency and dependability of this arrival fordetermining shear velocity. The interference effect iscompounded in that the initial shear energycommonly is not extremely high relative to P-waveenergy. Rather shear amplitudes are generally lowinitially and increase with succeeding arrivals. Theshear breaking point, therefore, almost always isobscured by P-wave interference. In somelithologies, such as low-porosity carbonates, an earlyshear arrival (probably the second or third shearhalf-cycle)sometimes has relatively high amplitudecompared with superimposed P arrivals. This"high-amplitude" event is commonly used to determineshear velocity. SPEJ P. 317^


2010 ◽  
Vol 2 (2) ◽  
pp. 307-329 ◽  
Author(s):  
C. Juhlin ◽  
B. Lund

Abstract. Reflection seismic data were acquired along a ca. 22 km long profile over the end-glacial Burträsk Fault with a nominal receiver and source spacing of 20 m. A steeply dipping reflection can be correlated to the Burträsk Fault, indicating that the fault dips at about 55° to the southeast near the surface. The reflection from the fault is rather poorly imaged, probably due to a jump in the fault and the crookedness of the seismic profile in the vicinity of the fault. A more pronounced steeply dipping reflection is observed about 4 km southeast of the Burträsk Fault. Based on its correlation with a topographic low at the surface this reflection is interpreted to originate from a fracture zone. There are no signs of large displacements along this fault as the glacial ice receded, but it may be active today. Other reflections on the processed seismic section may originate from changes in lithological variations in the supra-crustal rocks or from intrusions of more mafic rock. Constraints on the fault geometry provided by the reflection seismic data will help determine what stresses were required to activate the fault when the major rupture along it occurred.


Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. C107-C123
Author(s):  
Ivan Pšenčík ◽  
Bohuslav Růžek ◽  
Petr Jílek

We have developed a practical concept of compressional wave (P-wave) traveltime inversion in weakly to moderately anisotropic media of arbitrary symmetry and orientation. The concept provides sufficient freedom to explain and reproduce observed anisotropic seismic signatures to a high degree of accuracy. The key to this concept is the proposed P-wave anisotropy parameterization (A-parameters) that, together with the use of the weak-anisotropy approximation, leads to a significantly simplified theory. Here, as an example, we use a simple and transparent formula relating P-wave traveltimes to 15 P-wave A-parameters describing anisotropy of arbitrary symmetry. The formula is used in the inversion scheme, which does not require any a priori information about anisotropy symmetry and its orientation, and it is applicable to weak and moderate anisotropy. As the first step, we test applicability of the proposed scheme on a blind inversion of synthetic P-wave traveltimes generated in vertical seismic profile experiments in homogeneous models. Three models of varying anisotropy are used: tilted orthorhombic and triclinic models of moderate anisotropy (approximately 10%) and an orthorhombic model of strong anisotropy (>25%) with a horizontal plane of symmetry. In all cases, the inversion yields the complete set of 15 P-wave A-parameters, which make reconstruction of corresponding phase-velocity surfaces possible with high accuracy. The inversion scheme is robust with respect to noise and the source distribution pattern. Its quality depends on the angular illumination of the medium; we determine how the absence of nearly horizontal propagation directions affects inversion accuracy. The results of the inversion are applicable, for example, in migration or as a starting model for inversion methods, such as full-waveform inversion, if a model refinement is desired. A similar procedure could be designed for the inversion of S-wave traveltimes in anisotropic media of arbitrary symmetry.


Geophysics ◽  
1985 ◽  
Vol 50 (4) ◽  
pp. 530-538 ◽  
Author(s):  
P. M. Carrion ◽  
S. Hassanzadeh

Conventional velocity analysis of seismic data is based on normal moveout of common‐depth‐point (CDP) traveltime curves. Analysis is done in a hyperbolic framework and, therefore, is limited to using the small‐angle reflections only (muted data). Hence, it can estimate the interval velocities of compressional waves only, since mode conversion is negligible when small‐angle arrivals are concerned. We propose a new method which can estimate the interval velocities of compressional and mode‐converted waves separately. The method is based on slant stacking or plane‐wave decomposition (PWD) of the observed data (seismogram), which transforms the data from the conventional T-X domain into the intercept time‐ray parameter domain. Since PWD places most of the compressional energy into the precritical region of the slant‐stacked seismogram, the compressional‐wave interval velocities can be estimated using the “best ellipse” approximation on the assumption that the elliptic array velocity (stacking velocity) is approximately equal to the root‐mean‐square (rms) velocity. Similarly, shear‐wave interval velocities can be estimated by inverting the traveltime curves in the region of the PWD seismogram, where compressional waves decay exponentially (postcritical region). The method is illustrated by examples using synthetic and real data.


1983 ◽  
Vol 73 (3) ◽  
pp. 765-780
Author(s):  
Jia-Ju Lee ◽  
Charles A. Langston

abstract A three-dimensional ray method is used to compute three components of ground motion for complex structures involving curved boundaries. The method of principal curvature is developed to compute geometrical spreading of rays. This method, commonly used in electromagnetic wave propagation problems, employs phase matching at model interfaces and analysis of the wave front surface metric as the ray propagates throughout the model. It is an elegant way to examine the characteristics of three-dimensional caustics. Results computed for a two-dimensional canonical basin model with a plane SH-wave source are compared and are found to be in good agreement with those previously obtained by other independent numerical methods. Relaxing the restriction that the incident wave be perpendicular to the basin symmetry axis gives rise to large amplitude vertical and radial motions for incident SH waves and large tangential motions for incident P waves. As in the two-dimensional case, seismic energy is geometrically focused in the central region of the basin but strong later arrivals from the curved boundaries are not well developed in the three-dimensional case. The method is of direct use in analyzing three-dimensional crustal structure from off-azimuth P to S and S to P conversions.


Sign in / Sign up

Export Citation Format

Share Document