Rapid map and inversion of P-SV waves

Geophysics ◽  
1991 ◽  
Vol 56 (6) ◽  
pp. 859-862 ◽  
Author(s):  
Robert R. Stewart

Multicomponent seismic recordings are currently being analyzed in an attempt to improve conventional P‐wave sections and to find and use rock properties associated with shear waves (e.g. Dohr, 1985; Danbom and Dominico, 1986). Mode‐converted (P-SV) waves hold a special interest for several reasons: They are generated by conventional P‐wave sources and have only a one‐way travel path as a shear wave through the typically low velocity and attenuative near surface. For a given frequency, they will have a shorter wavelength than the original P wave, and thus offer higher spatial resolution; this has been observed in several vertical seismic profiling (VSP) cases (e.g., Geis et al., 1990). However, for surface seismic data, converted waves are often found to be of lower frequency than P-P waves (e.g., Eaton et al., 1991).

Geophysics ◽  
1982 ◽  
Vol 47 (7) ◽  
pp. 977-1000 ◽  
Author(s):  
C. C. Lash

A vertical wave test employing the vertical seismic profiling (VSP) technique in southern Mississippi confirmed suspicions that apparent multiple reflections might include converted waves as well as multiply reflected compressional waves. Both compressional (P) and shear (S) waves generated near the source were observed to travel to great depths, and P‐to‐S conversions were apparent in deep zones as well as shallow. P‐wave reflections were observed in agreement with predictions from synthetic records based on the sonic log. Up‐traveling P‐waves were reflected a short distance below the surface, at the base of the low‐velocity layer, and were followed as down‐traveling P‐waves to 200 ft depth by means of a vertical spread. Below 2000 ft and following the first P wave train, the predominate energy appeared to be down‐traveling P‐waves which could not be traced back to the reflection of up‐traveling P‐waves. The continuity of wavelets indicated instead that the strong down‐traveling S‐waves generated near the source produced P‐waves by S‐to‐P conversion somewhere in the zone between 800 and 1400 ft. The interference on the recordings made with an individual seismometer, or a small group of seismometers, using dynamite shots as the source was generally of a low‐frequency nature, so that the signal‐to‐noise (S/N) ratio was improved by the use of a high passband filter. The interference was greatly reduced without the need for a filter on recordings in which the source was a distributed charge of 100 ft length. The distributed charge produced much less shear‐wave energy in the P reflection band, demonstrating that the interference encountered when using a concentrated charge source was the consequence of the generation of S‐waves near the source. The distributed charges were previously chosen as a means for effectively eliminating secondary (ghost) reflections, an unwanted form of multiple reflections.


2021 ◽  
Vol 12 (1) ◽  
Author(s):  
Gaochun Wang ◽  
Hans Thybo ◽  
Irina M. Artemieva

AbstractAll models of the magmatic and plate tectonic processes that create continental crust predict the presence of a mafic lower crust. Earlier proposed crustal doubling in Tibet and the Himalayas by underthrusting of the Indian plate requires the presence of a mafic layer with high seismic P-wave velocity (Vp > 7.0 km/s) above the Moho. Our new seismic data demonstrates that some of the thickest crust on Earth in the middle Lhasa Terrane has exceptionally low velocity (Vp < 6.7 km/s) throughout the whole 80 km thick crust. Observed deep crustal earthquakes throughout the crustal column and thick lithosphere from seismic tomography imply low temperature crust. Therefore, the whole crust must consist of felsic rocks as any mafic layer would have high velocity unless the temperature of the crust were high. Our results form basis for alternative models for the formation of extremely thick juvenile crust with predominantly felsic composition in continental collision zones.


Geophysics ◽  
1988 ◽  
Vol 53 (1) ◽  
pp. 76-84 ◽  
Author(s):  
E. L. Majer ◽  
T. V. McEvilly ◽  
F. S. Eastwood ◽  
L. R. Myer

In a pilot vertical seismic profiling study, P-wave and cross‐polarized S-wave vibrators were used to investigate the potential utility of shear‐wave anisotropy measurements in characterizing a fractured rock mass. The caprock at The Geysers geothermal field was found to exhibit about an 11 percent velocity variation between SH-waves and SV-waves generated by rotating the S-wave vibrator orientation to two orthogonal polarizations for each survey level in the well. The effect is generally consistent with the equivalent anisotropy expected from the known fracture geometry.


2022 ◽  
Vol 41 (1) ◽  
pp. 40-46
Author(s):  
Öz Yilmaz ◽  
Kai Gao ◽  
Milos Delic ◽  
Jianghai Xia ◽  
Lianjie Huang ◽  
...  

We evaluate the performance of traveltime tomography and full-wave inversion (FWI) for near-surface modeling using the data from a shallow seismic field experiment. Eight boreholes up to 20-m depth have been drilled along the seismic line traverse to verify the accuracy of the P-wave velocity-depth model estimated by seismic inversion. The velocity-depth model of the soil column estimated by traveltime tomography is in good agreement with the borehole data. We used the traveltime tomography model as an initial model and performed FWI. Full-wave acoustic and elastic inversions, however, have failed to converge to a velocity-depth model that desirably should be a high-resolution version of the model estimated by traveltime tomography. Moreover, there are significant discrepancies between the estimated models and the borehole data. It is understandable why full-wave acoustic inversion would fail — land seismic data inherently are elastic wavefields. The question is: Why does full-wave elastic inversion also fail? The strategy to prevent full-wave elastic inversion of vertical-component geophone data trapped in a local minimum that results in a physically implausible near-surface model may be cascaded inversion. Specifically, we perform traveltime tomography to estimate a P-wave velocity-depth model for the near-surface and Rayleigh-wave inversion to estimate an S-wave velocity-depth model for the near-surface, then use the resulting pairs of models as the initial models for the subsequent full-wave elastic inversion. Nonetheless, as demonstrated by the field data example here, the elastic-wave inversion yields a near-surface solution that still is not in agreement with the borehole data. Here, we investigate the limitations of FWI applied to land seismic data for near-surface modeling.


Geophysics ◽  
1997 ◽  
Vol 62 (5) ◽  
pp. 1365-1368
Author(s):  
M. Boulfoul ◽  
Doyle R. Watts

The petroleum exploration industry uses S‐wave vertical seismic profiling (VSP) to determine S‐wave velocities from downgoing direct arrivals, and S‐wave reflectivities from upgoing waves. Seismic models for quantitative calibration of amplitude variation with offset (AVO) data require S‐wave velocity profiles (Castagna et al., 1993). Vertical summations (Hardage, 1983) of the upgoing waves produce S‐wave composite traces and enable interpretation of S‐wave seismic profile sections. In the simplest application of amplitude anomalies, the coincidence of high amplitude P‐wave reflectivity and low amplitude S‐wave reflectivity is potentially a direct indicator of the presence of natural gas.


Geophysics ◽  
2016 ◽  
Vol 81 (5) ◽  
pp. T221-T231 ◽  
Author(s):  
Christine E. Krohn ◽  
Thomas J. Murray

The top 6 m of the near surface has a surprisingly large effect on the behavior of P- and S-waves. For unconsolidated sediments, the P-wave velocity gradient and attenuation can be quite large. Computer modeling should include these properties to accurately reproduce seismic effects of the near surface. We have used reverse VSP data and computer simulations to demonstrate the following effects for upgoing P-waves. Near the surface, we have observed a large time delay, indicating low velocity ([Formula: see text]), and considerable pulse broadening, indicating high attenuation ([Formula: see text]). Consequently, shallowly buried geophones have greater high-frequency bandwidth compared with surface geophones. In addition, there is a large velocity gradient in the shallow near surface (factor of 10 in 5 m), resulting in the rotation of P-waves to the vertical with progressively smaller amplitudes recorded on horizontal phones. Finally, we have found little indication of a reflection or ghost from the surface, although downgoing reflections have been observed from interfaces within the near surface. In comparison, the following have been observed for upgoing S-waves: There is a small increase in the time delay or pulse broadening near the surface, indicating a smaller velocity gradient and less change in attenuation. In addition, the surface reflection coefficient is nearly one with a prominent surface ghost.


Author(s):  
Haohao Zhang ◽  
Jun Lu ◽  
Benchi Chen ◽  
Xuejun Ma ◽  
Zhidong Cai

Abstract The considerable depth and complicated structure of the Tahe Oilfield in the Tuofutai area of China make it very difficult to delineate its Ordovician carbonate fracture-cavity reservoir. The resolution of conventional ground seismic data is inadequate to satisfy current exploitation requirements. To further improve the understanding of the carbonate fracture-cavity reservoir of the Tahe Oilfield and to provide predictions of reservoir properties that are more accurate, a walkaround 3D-3C vertical seismic profiling (VSP) survey was conducted. First, we preprocessed raw VSP data and developed a method of joint PP- and PSV-wave prestack time migration. In contrast to ground seismic imaging profiles, VSP imaging profiles have a higher resolution and wider spectrum range that provide more detailed strata information. Then, using the joint PP- and PSV-wave prestack inversion method, we obtained the PP- and PSV-wave impedance and Poisson's ratio parameters of the Ordovician carbonate reservoir. Compared with the P-wave impedance of the ground seismic inversion, we found the VSP inversion results had higher accuracy, which enabled clearer identification of the internal characteristics of the carbonate reservoir. Finally, coupled with the Poisson's ratio attribute, we predicted the distribution of favorable reservoirs and interwell connectivity. The prediction results were verified using both logging and production data. The findings of this study demonstrate the applicability of the proposed technical method for the exploration of deep carbonate fracture-cavity reservoirs.


Geophysics ◽  
1952 ◽  
Vol 17 (2) ◽  
pp. 218-228 ◽  
Author(s):  
H. M. Thralls ◽  
R. W. Mossman

The arbitrary application of any set type of near‐surface corrections to seismic data can lead to erroneous results. The determination of the type of correction to be used must be based, in part, on the type of formations present in the near‐surface. Case studies are offered to illustrate conditions arising in areas of youthful and mature topography. Specifically, they deal with a complex low velocity layer problem in a river valley, a pre‐glacial topography in the Illinois Basin, a problem arising in a mature topography in Kansas, and a youthful topography in central Wyoming. In such cases, the use of a “floating” elevation reference plane is advocated for the “Correction Zone” lying immediately below the surface.


Geophysics ◽  
2002 ◽  
Vol 67 (6) ◽  
pp. 1877-1885 ◽  
Author(s):  
Xin‐Quan Ma

A new prestack inversion algorithm has been developed to simultaneously estimate acoustic and shear impedances from P‐wave reflection seismic data. The algorithm uses a global optimization procedure in the form of simulated annealing. The goal of optimization is to find a global minimum of the objective function, which includes the misfit between synthetic and observed prestack seismic data. During the iterative inversion process, the acoustic and shear impedance models are randomly perturbed, and the synthetic seismic data are calculated and compared with the observed seismic data. To increase stability, constraints have been built into the inversion algorithm, using the low‐frequency impedance and background Vs/Vp models. The inversion method has been successfully applied to synthetic and field data examples to produce acoustic and shear impedances comparable to log data of similar bandwidth. The estimated acoustic and shear impedances can be combined to derive other elastic parameters, which may be used for identifying of lithology and fluid content of reservoirs.


Author(s):  
Glenn Thackray ◽  
Mark Zellman ◽  
Jason Altekruse ◽  
Bruno Protti ◽  
Harrison Colandera

Preliminary results from seismic data collected at two sites on the Teton fault reveal shallow sub-surface fault structure and a basis for evaluating the post-glacial faulting record in greater detail. These new data include high-resolution shallow 2D seismic refraction and Interferometric Multi-Channel Analysis of Surface Waves (IMASW) (O’Connell and Turner 2010) depth-averaged shear wave velocity (Vs). The Teton fault, a down-to-the east normal fault, is expressed as a distinct topographic escarpment along the base of the eastern front of the Teton Range in Wyoming. The average fault scarp height cut into deglacial surfaces in several similar valleys and an assumed 14,000 yr BP deglaciation indicates an average postglacial offset rate of 0.82 m/ka (Thackray and Staley, in review). Because the fault is located almost entirely within Grand Teton National Park (GTNP), and in terrain that is remote and difficult to access, very few subsurface studies have been used to evaluate the fault. As a result, many uncertainties exist in the present characterization of along-strike slip rate, down-dip geometry, and rupture history, among other parameters. Additionally, questions remain about the fault dip at depth. Shallow seismic data were collected at two locations on the Teton fault scarp to (1) use a non-destructive, highly portable and cost-effective data collection system to image and characterize the Teton fault, (2) use the data to estimate vertical offsets of faulted bedrock and sediment, and (3) estimate fault dip in the shallow subsurface. Vs data were also collected at three GTNP facility structures to provide measured 30 m depth-averaged Vs (Vs30) for each site. Seismic data were collected using highly portable equipment packed into each site on foot. The system utilizes a sensor line 92 m long that includes 24 geophones (channels) at 4 m intervals. At both the Taggart Lake and String Lake sites, P-wave refraction data were collected spanning the fault scarp and perpendicular to local fault strike, as well as IMASW Vs seismic lines positioned on the hanging wall to provide Vs vs. Depth profiles crossing and perpendicular to the refraction survey lines. The Taggart Lake and String Lake 2D P-wave refraction profile and IMASW Vs plots reveal buried velocity structure that is vertically offset by the Teton fault. At Taggart Lake, we interpret the velocity horizon to be the top of dense glacial sediment (possibly compacted till), which is overlain by younger, slower, sediments. This surface is offset ~13 m (down-to-the-east) across the Teton fault. The vertical offset is in agreement with the measured height of the corresponding topographic scarp (~12 - 15 m). Geomorphic analysis of EarthScope (2008) LiDAR reveals small terraces, slope inflections and an abandoned channel on the footwall side of the scarp. At String Lake, the shallow buried velocity structure is inferred as unconsolidated alluvium (till, colluvium, alluvium); this relatively low velocity zone (


Sign in / Sign up

Export Citation Format

Share Document