Empirical studies of some of the seismic phenomena of Hawaii*

1938 ◽  
Vol 28 (4) ◽  
pp. 313-337
Author(s):  
Austin E. Jones
Keyword(s):  
P Wave ◽  
Secondary Peak ◽  
S Wave ◽  
S Waves ◽  
Mauna Loa ◽  
Earthquake Records ◽  
Type K ◽  
Sea Bottom ◽  
The Difference ◽  
The Waves

Summary and Conclusions A comparison was made of all the periods of local earthquakes entered in the record books, and this showed that the P wave of 0.3-sec. period occurred a maximum of 156 times, and a secondary peak for the period of 0.5 sec. occurred 89 times. The S wave of 0.5-sec. period had a maximum of 129 occurrences, and a secondary peak for 0.8-sec. period had 100. This suggested that in any earthquake the ratio of the period of the S to the P wave was inversely as their velocities, or as the square root of three. The maxima just given appear to hold for such waves from all depths of origin. It had been noted previously that large amplitudes and periods occur together. The upper limits of the amplitudes of the P and S waves of local shocks were found to vary with the cube of the periods. Different results were found for the variation of epicentral shocks in California and Japan. The difference may be caused by the difference in physical characteristics of the underlying crustal rock. While these studies in Hawaii were made on shocks of intensities I to IV, Rossi-Forel, they show promise of giving information about the waves to be expected in destructive earthquakes. The sectorial lines may be raised by new data, but in each region should approach some unknown lines as a limit. Formulas were used to correct the observed waves to those of standard displacement and consequent period. These periods were plotted with respect to distance and depth, with no reliable result. A tendency was shown for the period of P waves to increase with distance more rapidly than the period of S waves, whereas observations of more distant earthquakes would suggest the opposite. Study of the ratios of the amplitudes of the P to the S wave (AP/SS) showed no distance effect. The formulas from the previous amplitude-period study suggest that this ratio should not vary with the local distance. For Hawaii the ratio averages about 15 per cent. About 60 per cent of the foci are less than 5 km. deep, 70 per cent less than 10 km. deep. Very few appear to have originated at 60 or more km. depth. The decline in numbers of earthquakes with depth is a rapidly decreasing exponential function. Most of the deep earthquakes are under Mauna Loa and the Kilauea southeast rift zone. A large number of the located shallow foci are in and near the Kilauea crater. Possibly this is an increase that should be expected near any active volcanic crater, but it may be due to the close network of stations about Kilauea crater. The magnitude of the shock is not a function of the location either areally or in depth; that is, large earthquakes may be expected in any part of the island and near-by sea bottom and at all depths to at least 60 km. A method of classifying the earthquake records is based on the number of P or S waves shown on the seismogram, which indicate the key number from one to seven. A map of Hawaii was constructed showing the areas in which the different types of shock had originated. The first type, K-1, occurs either central to Mauna Loa or within 50 to 60 km. radius of the seismograph. Type K-2 is not recorded from northwest Hawaii. Type K-3 does not occur close to the instruments. Types K-4 to K-7 are noted to occur at somewhat greater distances, and to date have been observed only from small outlying areas. Earthquake records of simple character are generally near the area of deep-focus shocks and near the seismographs, so that the waves come in at a steep angle. Earthquakes under Kilauea crater are generally simple. As the foci become more distant and shallow they also become more complicated in type. These criteria should help in designating phases and consequent locations, but they are not final, and may be of no help beyond 100 km. The number of phases in some of the records of outlying earthquakes suggest a complexity of structure in the island mass and the near-by sea bottom. The locations near and on the extension of rifts and in pronounced lines and zones suggest a larger and more numerous system of rifts than has previously been mapped. The resulting pattern of rifts about Mauna Loa is roughly an asterisk. The main accent is on the visible active rifts to the southwest and the east-northeast of Mokuaweoweo. These rifts have apparently controlled most of the island's seismicity in the immediate past.

Azimuthal AVO for tilted TI medium

Geophysics ◽  
2017 ◽  
Vol 82 (1) ◽  
pp. C21-C33 ◽  
Author(s):  
Hongwei Wang ◽  
Suping Peng ◽  
Wenfeng Du
Keyword(s):  
Isotropic Medium ◽  
Symmetry Axis ◽  
Transversely Isotropic ◽  
P Wave ◽  
S Wave ◽  
Transverse Waves ◽  
Propagation Angle ◽  
S Waves ◽  
The Difference ◽  
Dip Angle

With the incident P-wave, we derive approximate formulas for amplitudes and polarizations of waves reflected from and transmitted through a planar, horizontal boundary between an overlying isotropic medium and an underlying tilted transversely isotropic (TTI) medium assuming that the directions of the phase and group velocities are consistent. Provided that the velocities in the isotropic medium are equal to the velocities along the symmetry axis direction, we derive the relational expression between the propagation angle in the TTI medium and the propagation angle in the hypothetical isotropic medium, under the condition that the horizontal slowness is the same, and then we update the approximate formula of the polarization in the TTI medium. Provided that the slow and fast transverse waves (qS and SH) are generated simultaneously in the anisotropic interface, we linearize for a six-order Zoeppritz equation, derive the azimuthal formula of longitudinal and S-waves, and determine their detailed expressions within the symmetry axis plane. According to the derived azimuthal AVO formula, we establish medium models, compare the derived AVO with the precision, and obtain the following conclusions: (1) The dip angle for the symmetry axis with respect to the vertical may have a sufficiently large impact on AVO, and the vertical longitudinal wave can generate an S-wave. (2) For the derived AVO formula, within the symmetry axis plane, the fitting effect of the approximate and exact formulas is good; however, within the other incident planes, taking the azimuth angle 45° as an example, the approximation is suitable for the large impedance contrast if the anisotropic parameters are set properly. (3) The error between the approximation and precision is mainly caused by the difference between the reflected and transmitted angles, the velocities’ derivation with respect to azimuth, and the division of approximation into isotropic and anisotropic parts.

Polarization and coherence of 5 to 30 Hz seismic wave fields at a hard-rock site and their relevance to velocity heterogeneities in the crust

1990 ◽  
Vol 80 (2) ◽  
pp. 430-449 ◽  
Author(s):  
William Menke ◽  
Arthur L. Lerner-Lam ◽  
Bruce Dubendorff ◽  
Javier Pacheco
Keyword(s):  
Hard Rock ◽  
Scale Effects ◽  
Spatial Coherence ◽  
Compressional Wave ◽  
P Wave ◽  
Transverse Motion ◽  
S Wave ◽  
Chemical Explosions ◽  
Aperture Arrays ◽  
The Waves

Abstract Except for its very onset, the P wave of earthquakes and chemical explosions observed at two narrow-aperture arrays on hard-rock sites in the Adirondack Mountains have a nearly random polarization. The amount of energy on the vertical, radial, and transverse components is about equal over the frequency range 5 to 30 Hz, for the entire seismogram. The spatial coherence of the seismograms is approximately exp(−cfΔx), where c is in the range 0.4 to 0.7 km−1Hz−1, f is frequency and Δx is the distance between array elements. Vertical, radial, and transverse components were quite coherent over the aperture of the array, indicating that the transverse motion of the compressional wave is a property of relatively large (106 m3) volumes of rock, and not just an anomaly caused by a malfunctioning instrument, poor instrument-rock coupling, or out-crop-scale effects. The spatial coherence is approximately independent of component, epicentral azimuth and range, and whether P- or S-wave coda is being considered, at least for propagation distances between 5 and 170 km. These results imply a strongly and three-dimensionally heterogeneous crust, with near-receiver scattering in the uppermost crust controlling the coherence properties of the waves.

Numerical simulation of azimuthal acoustic logging in a borehole penetrating a rock formation boundary

Geophysics ◽  
2016 ◽  
Vol 81 (3) ◽  
pp. D283-D291 ◽  
Author(s):  
Peng Liu ◽  
Wenxiao Qiao ◽  
Xiaohua Che ◽  
Xiaodong Ju ◽  
Junqiang Lu ◽  
...  
Keyword(s):  
Domain Model ◽  
P Wave ◽  
S Wave ◽  
Angle Variation ◽  
P Waves ◽  
Acoustic Logging ◽  
S Waves ◽  
Rock Formation ◽  
Logging Tool ◽  
Difference Time

We have developed a new 3D acoustic logging tool (3DAC). To examine the azimuthal resolution of 3DAC, we have evaluated a 3D finite-difference time-domain model to simulate a case in which the borehole penetrated a rock formation boundary when the tool worked at the azimuthal-transmitting-azimuthal-receiving mode. The results indicated that there were two types of P-waves with different slowness in waveforms: the P-wave of the harder rock (P1) and the P-wave of the softer rock (P2). The P1-wave can be observed in each azimuthal receiver, but the P2-wave appears only in the azimuthal receivers toward the softer rock. When these two types of rock are both fast formations, two types of S-waves also exist, and they have better azimuthal sensitivity compared with P-waves. The S-wave of the harder rock (S1) appears only in receivers toward the harder rock, and the S-wave of the softer rock (S2) appears only in receivers toward the softer rock. A model was simulated in which the boundary between shale and sand penetrated the borehole but not the borehole axis. The P-wave of shale and the S-wave of sand are azimuthally sensitive to the azimuth angle variation of two formations. In addition, waveforms obtained from 3DAC working at the monopole-transmitting-azimuthal-receiving mode indicate that the corresponding P-waves and S-waves are azimuthally sensitive, too. Finally, we have developed a field example of 3DAC to support our simulation results: The azimuthal variation of the P-wave slowness was observed and can thus be used to reflect the azimuthal heterogeneity of formations.

Green's function of the Lord–Shulman thermo-poroelasticity theory

2020 ◽  
Vol 221 (3) ◽  
pp. 1765-1776 ◽  
Author(s):  
Jia Wei ◽  
Li-Yun Fu ◽  
Zhi-Wei Wang ◽  
Jing Ba ◽  
José M Carcione
Keyword(s):  
Plane Wave ◽  
Green’S Function ◽  
Green's Function ◽  
Heterogeneous Media ◽  
Numerical Algorithms ◽  
Heat Diffusion ◽  
P Wave ◽  
Wave Analysis ◽  
S Wave ◽  
S Waves

SUMMARY The Lord–Shulman thermoelasticity theory combined with Biot equations of poroelasticity, describes wave dissipation due to fluid and heat flow. This theory avoids an unphysical behaviour of the thermoelastic waves present in the classical theory based on a parabolic heat equation, that is infinite velocity. A plane-wave analysis predicts four propagation modes: the classical P and S waves and two slow waves, namely, the Biot and thermal modes. We obtain the frequency-domain Green's function in homogeneous media as the displacements-temperature solution of the thermo-poroelasticity equations. The numerical examples validate the presence of the wave modes predicted by the plane-wave analysis. The S wave is not affected by heat diffusion, whereas the P wave shows an anelastic behaviour, and the slow modes present a diffusive behaviour depending on the viscosity, frequency and thermoelasticity properties. In heterogeneous media, the P wave undergoes mesoscopic attenuation through energy conversion to the slow modes. The Green's function is useful to study the physics in thermoelastic media and test numerical algorithms.

Scalar reverse‐time depth migration of prestack elastic seismic data

Geophysics ◽  
2001 ◽  
Vol 66 (5) ◽  
pp. 1519-1527 ◽  
Author(s):  
Robert Sun ◽  
George A. McMechan
Keyword(s):  
Seismic Waves ◽  
Velocity Model ◽  
P Wave ◽  
Common Source ◽  
Reverse Time ◽  
S Wave ◽  
Wave Source ◽  
Reflected Waves ◽  
S Waves ◽  
Elastic Data

Reflected P‐to‐P and P‐to‐S converted seismic waves in a two‐component elastic common‐source gather generated with a P‐wave source in a two‐dimensional model can be imaged by two independent scalar reverse‐time depth migrations. The inputs to migration are pure P‐ and S‐waves that are extracted by divergence and curl calculations during (shallow) extrapolation of the elastic data recorded at the earth’s surface. For both P‐to‐P and P‐to‐S converted reflected waves, the imaging time at each point is the P‐wave traveltime from the source to that point. The extracted P‐wave is reverse‐time extrapolated and imaged with a P‐velocity model, using a finite difference solution of the scalar wave equation. The extracted S‐wave is reverse‐time extrapolated and imaged similarly, but with an S‐velocity model. Converted S‐wave data requires a polarity correction prior to migration to ensure constructive interference between data from adjacent sources. Synthetic examples show that the algorithm gives satisfactory results for laterally inhomogeneous models.

Three-parameter prestack seismic inversion based on L1-2 minimization

Geophysics ◽  
2019 ◽  
Vol 84 (5) ◽  
pp. R753-R766 ◽  
Author(s):  
Lingqian Wang ◽  
Hui Zhou ◽  
Yufeng Wang ◽  
Bo Yu ◽  
Yuanpeng Zhang ◽  
...  
Keyword(s):  
Wave Velocity ◽  
Seismic Inversion ◽  
Thin Layers ◽  
P Wave ◽  
S Wave ◽  
Reservoir Prediction ◽  
Linear Inversion ◽  
Prestack Inversion ◽  
The Difference ◽  
Predictive Filtering

Prestack inversion has become a common approach in reservoir prediction. At present, the critical issue in the application of seismic inversion is the estimation of elastic parameters in the thin layers and weak reflectors. To improve the resolution and the accuracy of the inversion results, we introduced the difference of [Formula: see text] and [Formula: see text] norms as a nearly unbiased approximation of the sparsity of a vector, denoted as the [Formula: see text] norm, to the prestack inversion. The nonconvex penalty function of the [Formula: see text] norm can be decomposed into two convex subproblems via the difference of convex algorithm, and each subproblem can be solved efficiently by the alternating direction method of multipliers. Compared with the [Formula: see text] norm regularization, the [Formula: see text] minimization can reconstruct reflectivities more accurately. In addition, the [Formula: see text]-[Formula: see text] predictive filtering was introduced to guarantee the lateral continuity of the location and the amplitude of the reflectivity series. The generalized linear inversion and [Formula: see text]-[Formula: see text] predictive filtering are combined for stable elastic impedance inversion results, and three parameters of P-wave velocity, S-wave velocity, and density can be inverted with the Bayesian linearized amplitude variation with offset inversion. The inversion results of synthetic and real seismic data demonstrate that the proposed method can effectively improve the resolution and accuracy of the inversion results.

scholarly journals A new acoustic assumption for orthorhombic media

2020 ◽  
Vol 223 (2) ◽  
pp. 1118-1129
Author(s):  
Mohammad Mahdi Abedi ◽  
Alexey Stovas
Keyword(s):  
Wave Propagation ◽  
Conventional Approach ◽  
P Wave ◽  
Model Parameters ◽  
Governing Equations ◽  
S Wave ◽  
S Waves ◽  
Exploration Seismology ◽  
And Inversion ◽  
Approximate Equations

SUMMARY In exploration seismology, the acquisition, processing and inversion of P-wave data is a routine. However, in orthorhombic anisotropic media, the governing equations that describe the P-wave propagation are coupled with two S waves that are considered as redundant noise. The main approach to free the P-wave signal from the S-wave noise is the acoustic assumption on the wave propagation. The conventional acoustic assumption for orthorhombic media zeros out the S-wave velocities along three orthogonal axes, but leaves significant S-wave artefacts in all other directions. The new acoustic assumption that we propose mitigates the S-wave artefacts by zeroing out their velocities along the three orthogonal symmetry planes of orthorhombic media. Similar to the conventional approach, our method reduces the number of required model parameters from nine to six. As numerical experiments on multiple orthorhombic models show, the accuracy of the new acoustic assumption also compares well to the conventional approach. On the other hand, while the conventional acoustic assumption simplifies the governing equations, the new acoustic assumption further complicates them—an issue that emphasizes the necessity of simple approximate equations. Accordingly, we also propose simpler rational approximate phase-velocity and eikonal equations for the new acoustic orthorhombic media. We show a simple ray tracing example and find out that the proposed approximate equations are still highly accurate.

Comparison of P‐ and S‐wave velocities and Q’S from VSP and sonic log data

Geophysics ◽  
1994 ◽  
Vol 59 (10) ◽  
pp. 1512-1529 ◽  
Author(s):  
Gopa S. De ◽  
Donald F. Winterstein ◽  
Mark A. Meadows
Keyword(s):  
Velocity Dispersion ◽  
P Wave ◽  
S Wave ◽  
S Waves ◽  
Wave Velocities ◽  
Sonic Log ◽  
Transit Times ◽  
Wave Slope ◽  
Northern Alberta ◽  
Q Values

We compared P‐ and S‐wave velocities and quality factors (Q’S) from vertical seismic profiling (VSP) and sonic log measurements in five wells, three from the southwest San Joaquin Basin of California, one from near Laredo, Texas, and one from northern Alberta. Our purpose was to investigate the bias between sonic log and VSP velocities and to examine to what degree this bias might be a consequence of dispersion. VSPs and sonic logs were recorded in the same well in every case. Subsurface formations were predominantly clastic. The bias found was that VSP transit times were greater than sonic log times, consistent with normal dispersion. For the San Joaquin wells, differences in S‐wave transit times averaged 1–2 percent, while differences in P‐wave transit times averaged 6–7 percent. For the Alberta well, the situation was reversed, with differences in S‐wave transit times being about 6 percent, while those for P‐waves were 2.5 percent. For the Texas well, the differences averaged about 4 percent for both P‐ and S‐waves. Drift‐curve slopes for S‐waves tended to be low where the P‐wave slopes were high and vice versa. S‐wave drift‐curve slopes in the shallow California wells were 5–10 μs/ft (16–33 μs/m) and the P‐wave slopes were 15–30 μs/ft (49–98 μs/m). The S‐wave slope in sandstones in the northern Alberta well was up to 50 μs/ft (164 μs/m), while the P‐wave slope was about 5 μs/ft (16 μs/m). In the northern Alberta well the slopes for both P‐ and S‐waves flattened in the carbonate. In the Texas well, both P‐ and S‐wave drifts were comparable. We calculated (Q’s) from a velocity dispersion formula and from spectral ratios. When the two Q’s agreed, we concluded that velocity dispersion resulted solely from absorption. These Q estimation methods were reliable only for Q values smaller than 20. We found that, even with data of generally outstanding quality, Q values determined by standard methods can have large uncertainties, and negative Q’s may be common.

Velocity anisotropy in shale determined from crosshole seismic and vertical seismic profile data

Geophysics ◽  
1990 ◽  
Vol 55 (4) ◽  
pp. 470-479 ◽  
Author(s):  
D. F. Winterstein ◽  
B. N. P. Paulsson
Keyword(s):  
Seismic Profile ◽  
P Wave ◽  
Isotropic Model ◽  
S Wave ◽  
Vertical Seismic Profile ◽  
Near Surface ◽  
S Waves ◽  
Velocity Anisotropy ◽  
Wave Velocities ◽  
Late Tertiary

Crosshole and vertical seismic profile (VST) data made possible accurate characterization of the elastic properties, including noticeable velocity anisotropy, of a near‐surface late Tertiary shale formation. Shear‐wave splitting was obvious in both crosshole and VSP data. In crosshole data, two orthologonally polarrized shear (S) waves arrived 19 ms in the uppermost 246 ft (75 m). Vertically traveling S waves of the VSP separated about 10 ms in the uppermost 300 ft (90 m) but remained at nearly constant separation below that level. A transversely isotropic model, which incorporates a rapid increase in S-wave velocities with depth but slow increase in P-wave velocities, closely fits the data over most of the measured interval. Elastic constants of the transvesely isotropic model show spherical P- and [Formula: see text]wave velocity surfaces but an ellipsoidal [Formula: see text]wave surface with a ratio of major to minor axes of 1.15. The magnitude of this S-wave anisotropy is consistent with and lends credence to S-wave anisotropy magnitudes deduced less directly from data of many sedimentary basins.

scholarly journals Travel-Time Inversion Method of Converted Shear Waves Using RayInvr Algorithm

2021 ◽  
Vol 11 (8) ◽  
pp. 3571
Author(s):  
Genggeng Wen ◽  
Kuiyuan Wan ◽  
Shaohong Xia ◽  
Huilong Xu ◽  
Chaoyan Fan ◽  
...  
Keyword(s):  
Travel Time ◽  
Wave Velocity ◽  
P Wave ◽  
Inversion Method ◽  
Time Inversion ◽  
Ocean Bottom Seismometer ◽  
S Wave ◽  
S Waves ◽  
Inversion Model ◽  
Travel Time Inversion

The detailed studies of converted S-waves recorded on the Ocean Bottom Seismometer (OBS) can provide evidence for constraining lithology and geophysical properties. However, the research of converted S-waves remains a weakness, especially the S-waves’ inversion. In this study, we applied a travel-time inversion method of converted S-waves to obtain the crustal S-wave velocity along the profile NS5. The velocities of the crust are determined by the following four aspects: (1) modelling the P-wave velocity, (2) constrained sediments Vp/Vs ratios and S-wave velocity using PPS phases, (3) the correction of PSS phases’ travel-time, and (4) appropriate parameters and initial model are selected for inversion. Our results show that the vs. and Vp/Vs of the crust are 3.0–4.4 km/s and 1.71–1.80, respectively. The inversion model has a similar trend in velocity and Vp/Vs ratios with the forward model, due to a small difference with ∆Vs of 0.1 km/s and ∆Vp/Vs of 0.03 between two models. In addition, the high-resolution inversion model has revealed many details of the crustal structures, including magma conduits, which further supports our method as feasible.

Sign in / Sign up

Export Citation Format

Share Document