Using P- to S- wave conversions from controlled sources to determine the shear-wave velocity structure along Hikurangi Margin Forearc, New Zealand

Author(s):  
Pasan Herath ◽  
Tim Stern ◽  
Martha Savage ◽  
Dan Bassett ◽  
Stuart Henrys ◽  
...  

<p>The Hikurangi subduction margin offshore of the east coast of New Zealand displays along-strike variations in subduction-thrust slip behavior. Geodetic observations show that the subduction-thrust of the southern segment of the margin is locked on the 30-100 year scale and the northern segment displays periodic slow-slip on the 1-2 year scale. It is hypothesised that spatial variations in pore-pressure may play a role in this contrasting phenomenon. Higher pore-pressures would result in lower effective stresses, which promote slow-slip of the subduction-thrust. In addition, the presence of a sedimentary wedge with very low shear wave-speeds in the northern Hikurangi margin has been proposed to fit the ultra-long duration of ground motions observed following the 2016 Kaikoura earthquake. Compressional (P-) wave velocities (V<sub>p</sub>) of the subsurface provide useful information about the lithological composition. Combined with shear (S-) wave velocities (V<sub>s</sub>), the V<sub>p</sub>/V<sub>s</sub> ratio which is directly related to Poisson’s ratio can be obtained. This is a diagnostic property of a rock’s consolidation and porosity. Typical V<sub>p</sub>/V<sub>s</sub> ratio of consolidated and crystalline rocks range from 1.6 to 1.9 and that of unconsolidated sediments can range from 2.0 to 4.0.</p><p>We use the controlled sources of R/V Marcus G Langseth recorded by a profile of 49 multi-component ocean bottom seismometers (OBS) along the Hikurangi margin forearc for the Seismogenesis at Hikurangi Integrated Research Experiment (SHIRE) to derive the V<sub>s</sub> structure and estimate the V<sub>p</sub>/V<sub>s </sub>ratio. The orientations of the horizontal components of each OBS are found by a hodogram analysis and by an eigenvalue-decomposition of the covariance matrix. Using the orientations, the horizontal components of each OBS are rotated into radial and transverse components. P to S converted phases are identified on the radial and transverse components considering their linear moveout, polarisation angle, and ellipticity. We confirm incoming S-waves to OBSs by comparing them with their hydrophone components. We identify both PPS (up-going P-wave after reflection or refraction converts to an S-wave at an interface) and PSS (down-going P-wave from the controlled source converts to an S-wave at an interface) type conversions. The identified conversion interfaces are the sediment-basement interface and the top of the subducting crust. The travel-time delay of a PPS type conversion relative to its P-wave arrival is indicative of V<sub>s</sub> above the converting interface. The linear-moveout of PSS type conversions are indicative of V<sub>s</sub> along the raypath after the conversion. Preliminary results from the southern Hikurangi margin suggest V<sub>p</sub>/V<sub>s</sub> ratios of ~1.70 for the basement rocks above the subducting crust and ~1.90 for the sediments overlying the basement rocks. These values indicate that the basement rocks are consolidated and less porous than the overlying sediments.</p><p>We expect to estimate the V<sub>p</sub>/V<sub>s</sub> ratios in the northern Hikurangi margin to assess the role played by pore-pressure in the along-strike variation in subduction-thrust slip behavior. We also expect to ascertain the presence and estimate the thickness of the low-velocity sediment wedge in the northern Hikurangi margin.</p>

Geophysics ◽  
2006 ◽  
Vol 71 (6) ◽  
pp. F161-F164 ◽  
Author(s):  
Myung W. Lee

Prediction of shear-wave velocity plays an important role in seismic modeling, amplitude analysis with offset, and other exploration applications. This paper presents a method for predicting S-wave velocity from the P-wave velocity on the basis of the moduli of dry rock. Elastic velocities of water-saturated sediments at low frequencies can be predicted from the moduli of dry rock by using Gassmann’s equation; hence, if the moduli of dry rock can be estimated from P-wave velocities, then S-wave velocities easily can be predicted from the moduli. Dry rock bulk modulus can be related to the shear modulus through a compaction constant. The numerical results indicate that the predicted S-wave velocities for consolidated and unconsolidated sediments agree well with measured velocities if differential pressure is greater than approximately [Formula: see text]. An advantage of this method is that there are no adjustable parameters to be chosen, such as the pore-aspect ratios required in some other methods. The predicted S-wave velocity depends only on the measured P-wave velocity and porosity.


2003 ◽  
Author(s):  
Dan Ebrom ◽  
Phil Heppard ◽  
Mike Mueller ◽  
Leon Thomsen

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.


1996 ◽  
Vol 86 (6) ◽  
pp. 1704-1713 ◽  
Author(s):  
R. D. Catchings ◽  
W. H. K. Lee

Abstract The 17 January 1994, Northridge, California, earthquake produced strong ground shaking at the Cedar Hills Nursery (referred to here as the Tarzana site) within the city of Tarzana, California, approximately 6 km from the epicenter of the mainshock. Although the Tarzana site is on a hill and is a rock site, accelerations of approximately 1.78 g horizontally and 1.2 g vertically at the Tarzana site are among the highest ever instrumentally recorded for an earthquake. To investigate possible site effects at the Tarzana site, we used explosive-source seismic refraction data to determine the shallow (<70 m) P-and S-wave velocity structure. Our seismic velocity models for the Tarzana site indicate that the local velocity structure may have contributed significantly to the observed shaking. P-wave velocities range from 0.9 to 1.65 km/sec, and S-wave velocities range from 0.20 and 0.6 km/sec for the upper 70 m. We also found evidence for a local S-wave low-velocity zone (LVZ) beneath the top of the hill. The LVZ underlies a CDMG strong-motion recording site at depths between 25 and 60 m below ground surface (BGS). Our velocity model is consistent with the near-surface (<30 m) P- and S-wave velocities and Poisson's ratios measured in a nearby (<30 m) borehole. High Poisson's ratios (0.477 to 0.494) and S-wave attenuation within the LVZ suggest that the LVZ may be composed of highly saturated shales of the Modelo Formation. Because the lateral dimensions of the LVZ approximately correspond to the areas of strongest shaking, we suggest that the highly saturated zone may have contributed to localized strong shaking. Rock sites are generally considered to be ideal locations for site response in urban areas; however, localized, highly saturated rock sites may be a hazard in urban areas that requires further investigation.


1988 ◽  
Vol 11 ◽  
pp. 198 ◽  
Author(s):  
S. Anandakrishnan

Detailed seismic short-refraction profiling was conducted on Ice Stream Β (UpB) during the 1983–84 austral summer. A new high-resolution data logger, developed at the University of Wisconsin, recorded both compressional- and shear-wave arrivals. We report here on P-wave and S-wave profiles recorded along a line parallel to the axis of the ice stream. Source-receiver separations up to 720 m yielded seismic velocity-depth curves to below the firn-ice transition zone (slightly greater than 30 m at UpB). For the compressional-wave profile, geophones were separated by 2.5 m, which yielded a velocity-depth curve with a granularity of ∼1 m. The corresponding density-depth curve agrees well with direct density measurements obtained from a core extracted nearby (Alley and Bentley 1988, this volume). Discontinuities in the velocity gradient do not appear at the “critical densities” as they did at Byrd Station, Antarctica, and elsewhere (Kohnen and Bentley 1973 , Robertson and Bentley 1975). Two shear-wave profiles were recorded, both with geophone spacings of 5 m, one with longitudinal polarization (SV) and the other with transverse polarization (SH). There is a marked difference in velocity between the SH and SV waves, particularly in the shallow firn. We suggest that a strong vertical shape-and-bonding fabric in the shallow firn, as observed in cores collected at UpB, would account for this disparity.


2019 ◽  
Vol 283 ◽  
pp. 06003
Author(s):  
Guangxue Zheng ◽  
Hanhao Zhu ◽  
Jun Zhu

A method of geo-acoustic parameter inversion based on the Bayesian theory is proposed for the acquisition of acoustic parameters in shallow sea with the elastic seabed. Firstly, the theoretical prediction value of the sound pressure field is calculated by the fast field method (FFM). According to the Bayesian theory, we establish the misfit function between the measured sound pressure field and the theoretical pressure field. It is under the assumption of Gaussian data errors which are in line with the likelihood function. Finally, the posterior probability density (PPD) of parameters is given as the result of inversion. Our research is conducted in the light of Metropolis sample rules. Apart from numerical simulations, a scaled model experiment has been taken in the laboratory tank. The results of numerical simulations and tank experiments show that sound pressure field calculated by the result of inversion is consistent with the measured sound pressure field. Besides, s-wave velocities, p-wave velocities and seafloor density have fewer uncertainties and are more sensitive to complex sound pressure than s-wave attenuation and p-wave attenuation. The received signals calculated by inversion results are keeping with received signals in the experiment which verify the effectiveness of this method.


2002 ◽  
Vol 42 (1) ◽  
pp. 627
Author(s):  
R.G. Williams ◽  
G. Roberts ◽  
K. Hawkins

Seismic energy that has been mode converted from pwave to s-wave in the sub-surface may be recorded by multi-component surveys to obtain information about the elastic properties of the earth. Since the energy converted to s-wave is missing from the p-wave an alternative to recording OBC multi-component data is to examine p-wave data for the missing energy. Since pwave velocities are generally faster than s-wave velocities, then for a given reflection point the converted s-wave signal reaches the surface at a shorter offset than the equivalent p-wave information. Thus, it is necessary to record longer offsets for p-wave data than for multicomponent data in order to measure the same information.A non-linear, wide-angle (including post critical) AVO inversion has been developed that allows relative changes in p-wave velocities, s-wave velocities and density to be extracted from long offset p-wave data. To extract amplitudes at long offsets for this inversion it is necessary to image the data correctly, including correcting for higher order moveout and possibly anisotropy if it is present.The higher order moveout may itself be inverted to yield additional information about the anisotropy of the sub-surface.


2015 ◽  
Vol 76 (2) ◽  
Author(s):  
Badee Alshameri ◽  
Ismail Bakar ◽  
Aziman Madun ◽  
Edy Tonnizam Mohamad

One of the main geophysical tools (seismic tools) in the laboratory is the bender element. This tool can be used to measure some dynamic soil properties (e.g. shear and Young’s modulus). However, even if it relatively simple to use the bender element, inconsistent testing procedures can cause poor quality in the bender element data. One of the bender element procedure that always neglected is the alignment (different positions of bender element receiver to the transmitter in the vertical axis). The alignment effect was evaluated via changing the horizontal distance between transmitter and receiver starting from 0 to 110 mm for two sizes of the sample's thickness (i.e. 63.17 mm and 91.51 mm). Five methods were applied to calculate the travel times. Those methods were as the following: visually, first-peak, maximum-peak, CCexcel and CCGDS. In general, the experiments indicated uncertain results for both of the P-wave (primary wave) and S-wave (secondary wave) velocities at zone of Dr:D above 0.5:1 (where Dr is the horizontal distance of the receiver from the vertical axis and D is the thickness of the sample). On the other hand, both the visual and first-peak methods show the wave velocities results are higher than obtained from other methods. However, the ratio between the amplitude of transmitter signals to receiver amplitude signal was taken to calculate the damping-slope of the P-wave and S-wave. Thus the results from damping slope show steeply slope when the ratio of  Dr:D is above 0.5:1 compare with gentle slope below ratio 0.5:1 at the sample with thickness equal to 91.51 mm, while there is no variation at a slope in sample with thickness equal to 63.17 mm.


Geophysics ◽  
2020 ◽  
Vol 85 (5) ◽  
pp. MR297-MR308
Author(s):  
Roberta Ruggieri ◽  
Fabio Trippetta

Variations in reservoir seismic properties can be correlated to changes in saturated-fluid properties. Thus, the determination of variation in petrophysical properties of carbonate-bearing rocks is of interest to the oil exploration industry because unconventional oils, such as bitumen (HHC), are emerging as an alternative hydrocarbon reserve. We have investigated the temperature effects on laboratory seismic wave velocities of HHC-bearing carbonate rocks belonging to the Bolognano Formation (Majella Mountain, central Italy), which can be defined as a natural laboratory to study carbonate reservoir properties. We conduct an initial characterization in terms of porosity and density for the carbonate-bearing samples and then density and viscosity measurements for the residual HHC, extracted by HCl dissolution of the hosting rock. Acoustic wave velocities are recorded from ambient temperature to 90°C. Our acoustic velocity data point out an inverse relationship with temperature, and compressional (P) and shear (S) wave velocities show a distinct trend with increasing temperature depending on the amount of HHC content. Indeed, samples with the highest HHC content show a larger gradient of velocity changes in the temperature range of approximately 50°C–60°C, suggesting that the bitumen can be in a fluid state. Conversely, below approximately 50°C, the velocity gradient is lower because, at this temperature, bitumen can change its phase in a solid state. We also propose a theoretical model to predict the P-wave velocity change at different initial porosities for HHC-saturated samples suggesting that the velocity change mainly is related to the absolute volume of HHC.


Sign in / Sign up

Export Citation Format

Share Document