scholarly journals Wave-equation-based travel-time seismic tomography – Part 2: Application to the 1992 Landers earthquake (<i>M</i><sub>w</sub> 7.3) area

Solid Earth ◽  
2014 ◽  
Vol 5 (2) ◽  
pp. 1169-1188 ◽  
Author(s):  
P. Tong ◽  
D. Zhao ◽  
D. Yang ◽  
X. Yang ◽  
J. Chen ◽  
...  

Abstract. High-resolution 3-D P and S wave crustal velocity and Poisson's ratio models of the 1992 Landers earthquake (Mw 7.3) area are determined iteratively by a wave-equation-based travel-time seismic tomography (WETST) technique. The details of data selection, synthetic arrival-time determination, and trade-off analysis of damping and smoothing parameters are presented to show the performance of this new tomographic inversion method. A total of 78 523 P wave and 46 999 S wave high-quality arrival-time data from 2041 local earthquakes recorded by 275 stations during the period of 1992–2013 are used to obtain the final tomographic models, which cost around 10 000 CPU hours. Checkerboard resolution tests are conducted to verify the reliability of inversion results for the chosen seismic data and the wave-equation-based travel-time seismic tomography method. Significant structural heterogeneities are revealed in the crust of the 1992 Landers earthquake area which may be closely related to the local seismic activities. Strong variations of velocity and Poisson's ratio exist in the source regions of the Landers and three other nearby strong earthquakes. Most seismicity occurs in areas with high-velocity and low Poisson's ratio, which may be associated with the seismogenic layer. Pronounced low-velocity anomalies revealed in the lower crust along the Elsinore, the San Jacinto, and the San Andreas faults may reflect the existence of fluids in the lower crust. The recovery of these strong heterogeneous structures is facilitated by the use of full wave equation solvers and WETST and verifies their ability in generating high-resolution tomographic models.

2014 ◽  
Vol 6 (2) ◽  
pp. 2567-2613 ◽  
Author(s):  
P. Tong ◽  
D. Zhao ◽  
D. Yang ◽  
X. Yang ◽  
J. Chen ◽  
...  

Abstract. High-resolution 3-D P and S wave crustal velocity and Poisson's ratio models of the 1992 Landers earthquake (Mw 7.3) area are determined iteratively by a wave-equation based traveltime seismic tomography (WETST) technique as developed in the first paper. The details of data selection, synthetic arrival-time determination, and trade-off analysis of damping and smoothing parameters are presented to show the performance of this new tomographic inversion method. A total of 78 523 P wave and 46 999 S wave high-quality arrival-time data from 2041 local earthquakes recorded by 275 stations during the period of 1992–2013 is used to obtain the final tomographic models which costs around 10 000 CPU h. Checkerboard resolution tests are conducted to verify the reliability of inversion results for the chosen seismic data and the wave-equation based traveltime seismic tomography method. Significant structural heterogeneities are revealed in the crust of the 1992 Lander earthquake area which may be closely related to the local seismic activities. Strong variations of velocity and Poisson's ratio exist in the source regions of the Landers and three other strong earthquakes in this area. Most seismicity occurs in areas with high-velocity and low Poisson's ratio, which may be associated with the seismogenic layer. Pronounced low-velocity anomalies revealed in the lower crust along the Elsinore, the San Jacinto and the San Andreas faults may reflect the existence of fluids in the lower crust. The recovery of these strong heterogeneous structures are facilitated by the use of full wave equation solvers and WETST and verifies their ability in generating high-resolution tomographic models.


Geophysics ◽  
2018 ◽  
Vol 83 (4) ◽  
pp. T209-T234 ◽  
Author(s):  
Jing-Bo Chen ◽  
Jian Cao

Because of its high computational cost, we needed to develop an efficient numerical scheme for the frequency-domain 3D elastic wave equation. In addition, the numerical scheme should be applicable to media with a liquid-solid interface. To address these two issues, we have developed a new average-derivative optimal 27-point scheme with arbitrary directional grid intervals and a corresponding numerical dispersion analysis for the frequency-domain 3D elastic wave equation. The novelty of this scheme is that its optimal coefficients depend on the ratio of the directional grid intervals and Poisson’s ratio. In this way, this scheme can be applied to media with a liquid-solid interface and a computational grid with arbitrary directional grid intervals. For media with a variable Poisson’s ratio, we have developed an effective and stable interpolation method for optimization coefficients. Compared with the classic 19-point scheme, this new scheme reduces the required number of grid points per wavelength for equal and unequal directional grid intervals. The reduction of the number of grid points increases as the Poisson’s ratio becomes larger. In particular, the numerical S-wave phase velocity of this new scheme becomes zero, whereas the classic 19-point scheme produces a spurious numerical S-wave phase velocity, as Poisson’s ratio reaches 0.5. We have performed numerical examples to develop the theoretical analysis.


2005 ◽  
Vol 18 (3) ◽  
pp. 345-353 ◽  
Author(s):  
Zhi Liu ◽  
Xian-kang Zhang ◽  
Fu-yun Wang ◽  
Yong-hong Duan ◽  
Xiao-ling Lai

1992 ◽  
Vol 97 (B13) ◽  
pp. 19993 ◽  
Author(s):  
D. J. White ◽  
B. Milkereit ◽  
M. H. Salisbury ◽  
J. A. Percival

Geophysics ◽  
2000 ◽  
Vol 65 (2) ◽  
pp. 559-564 ◽  
Author(s):  
Ran Bachrach ◽  
Jack Dvorkin ◽  
Amos M. Nur

We determined P- and S-wave velocity depth profiles in shallow, unconsolidated beach sand by analyzing three‐component surface seismic data. P- and S-wave velocity profiles were calculated from traveltime measurements of vertical and tangential component seismograms, respectively. The results reveal two discrepancies between theory and data. Whereas both velocities were found to be proportional to the pressure raised to the power of 1/6, as predicted by the Hertz‐Mindlin contact theory, the actual values of the velocities are less than half of those calculated from this theory. We attribute this discrepancy to the angularity of the sand grains. Assuming that the average radii of curvature at the grain contacts are smaller than the average radii of the grains, we modify the Hertz‐Mindlin theory accordingly. We found that the ratio of the contact radius to the grain radius is about 0.086. The second disparity is between the observed Poisson’s ratio of 0.15 and the theoretical value (0.008 for random pack of quartz spheres). This discrepancy can be reconciled by assuming slip at the grain contacts. Because slip decreases the shearing between grains, Poisson’s ratio increases.


1969 ◽  
Vol 59 (3) ◽  
pp. 1189-1200
Author(s):  
John R. McGinley ◽  
Don L. Anderson

abstract The unified magnitude, the ratio of the amiplitudes of S to P waves, and travel-time residuals were compiled from published data for the five Seismological Observatories, TFO, UBO, BMO, WMO and CBO. Using one of the stations as a reference, a relative measure of the above quantities was calculated for each of the other stations for each of a number of earthquakes. The stations in the Basin and Range Province are consistent with a markedly higher attentuation of P waves and a high attenuation of S relative to P when compared to the other stations. This latter observation indicates a high Poisson's ratio in the mantle under the Basin and Range. The delay times to these stations are also consistent with the high Poisson's ratio and with a low-velocity upper mantle. The ratio of the amplitudes of long-period S waves to short-period P waves varies by a factor of 4 among these stations. BMO, in eastern Oregon, has a high S/P amplitude ratio compared to other stations and a travel-time residual that is comparable to the observatories in the mid-continent. This may be another example of a seismic “window” into the upper mantle that is generated by underthrusting of the oceanic lithosphere.


Energies ◽  
2021 ◽  
Vol 14 (16) ◽  
pp. 5007
Author(s):  
Stian Rørheim ◽  
Mohammad Hossain Bhuiyan ◽  
Andreas Bauer ◽  
Pierre Rolf Cerasi

Carbon capture and storage (CCS) by geological sequestration comprises a permeable formation (reservoir) for CO2 storage topped by an impermeable formation (caprock). Time-lapse (4D) seismic is used to map CO2 movement in the subsurface: CO2 migration into the caprock might change its properties and thus impact its integrity. Simultaneous forced-oscillation and pulse-transmission measurements are combined to quantify Young’s modulus and Poisson’s ratio as well as P- and S-wave velocity changes in the absence and in the presence of CO2 at constant seismic and ultrasonic frequencies. This combination is the laboratory proxy to 4D seismic because rock properties are monitored over time. It also improves the understanding of frequency-dependent (dispersive) properties needed for comparing in-situ and laboratory measurements. To verify our method, Draupne Shale is monitored during three consecutive fluid exposure phases. This shale appears to be resilient to CO2 exposure as its integrity is neither compromised by notable Young’s modulus and Poisson’s ratio nor P- and S-wave velocity changes. No significant changes in Young’s modulus and Poisson’s ratio seismic dispersion are observed. This absence of notable changes in rock properties is attributed to Draupne being a calcite-poor shale resilient to acidic CO2-bearing brine that may be a suitable candidate for CCS.


Sign in / Sign up

Export Citation Format

Share Document