Anisotropic velocity tomography: A case study in a near‐surface rock mass

Geophysics ◽  
1993 ◽  
Vol 58 (12) ◽  
pp. 1748-1763 ◽  
Author(s):  
R. G. Pratt ◽  
W. J. McGaughey ◽  
C. H. Chapman
Keyword(s):  
Rock Mass ◽  
Fault Zone ◽  
Seismic Velocity ◽  
Anisotropic Media ◽  
Rock Properties ◽  
Initial Assessment ◽  
Borehole Data ◽  
Near Surface ◽  
Crown Pillar ◽  
Velocity Tomography

Cross‐borehole data were acquired in the surface crown pillar of a massive sulfide ore mine. The data consist of five, two‐dimensional (2-D), cross‐borehole panels, each with approximately 900 source‐receiver pairs. The panels were located within the crown pillar at either side of and within a major subvertical fault zone that intersects the orebody. An initial analysis of the data indicates that the bedrock containing the orebody is seismically anisotropic. A rigorous analysis of the traveltimes using anisotropic velocity tomography confirms the initial assessment that anisotropy exists within the crown pillar rock mass. Anisotropic velocity tomography is the generalization of tomographic methods to anisotropic media. As in any geophysical problem, the data are insufficient to completely resolve the distributions of the rock properties at all scale lengths; we use external constraints on the roughness of the final solution to ensure an algebraically well‐posed problem. Plots of the data residuals (the “traveltime surfaces”) are an essential tool in determining an optimal level of constraint. Of equal importance are plots of the relationship between the solution roughness and the rms level of the residuals. The final results of anisotropic velocity tomography are a set of images (tomograms) of the velocity and selected anisotropy parameters for the five panels. Our images do not contain the distortions typically exhibited when using isotropic tomography in anisotropic media. The velocity tomograms clearly show the geometry of the overburden contact at the top of the bedrock. The anisotropy tomograms show a decrease in anisotropy with depth on two of the panels. They also show a decrease in anisotropy with proximity to the fault zone. These features of the seismic velocity anisotropy are consistent with observations of fracture orientation and distribution. The results of the crosshole data interpretation contribute to the overall site investigation and provide a reliable interrogation of the bulk properties of the rock mass.

scholarly journals Characterization of an Intraplate Seismogenic Zone Using Geophysical and Borehole Data: The Vila Franca de Xira Fault, Portugal

2020 ◽  
Vol 91 (4) ◽  
pp. 2287-2297
Author(s):  
João Carvalho ◽  
Daniela Alves ◽  
João Cabral ◽  
Ranajit Ghose ◽  
José Borges ◽  
...  
Keyword(s):  
Fault Zone ◽  
P Wave ◽  
Seismogenic Zone ◽  
S Wave ◽  
Borehole Data ◽  
Near Surface ◽  
Lower Tagus Valley ◽  
Seismogenic Source ◽  
Neogene Sediments ◽  
Tagus Valley

Abstract The Vila Franca de Xira (VFX) fault is a regional fault zone located about 25 km northeast of Lisbon, affecting Neogene sediments. Recent shear-wave seismic studies show that this complex fault zone is buried beneath Holocene sediments and is deforming the alluvial cover, in agreement with a previous work that proposes the fault as the source of the 1531 Lower Tagus Valley earthquake. In this work, we corroborate these results using S-wave, P-wave, geoelectric, ground-penetrating radar and borehole data, confirming that the sediments deformed by several fault branches are of Upper Pleistocene to Holocene. Accumulated fault vertical offsets of about 3 m are estimated from the integrated interpretation of geophysical and borehole data, including 2D elastic seismic modeling, with an estimated resolution of about 0.5 m. The deformations affecting the Tagus alluvial sediments probably resulted from surface or near-surface rupture of the VFX fault during M∼7 earthquakes, reinforcing the fault as the seismogenic source of regional historical events, as in 1531, and highlighting the need for preparedness for the next event.

Time variance of seismic velocity from multiple explosive sources southeast of Hollister, California

1977 ◽  
Vol 67 (5) ◽  
pp. 1339-1354 ◽  
Author(s):  
L. G. Peake ◽  
J. H. Healy ◽  
J. C. Roller
Keyword(s):  
Fault Zone ◽  
San Andreas Fault ◽  
Seismic Velocity ◽  
Temporal Distribution ◽  
Rock Properties ◽  
Arrival Times ◽  
Time Variance ◽  
San Andreas ◽  
Shot Point ◽  
Explosive Charges

abstract As part of an intensive effort to monitor the section of the San Andreas fault south of Hollister, California, for any seismic velocity time variance associated with earthquakes, 17 explosive charges were fired between August 1974 and May 1975, at a point 35 km east of the fault zone and recorded by a special network of 51 seismometers within and west of the San Andreas fault zone. The shot point was selected to minimize differences in rock properties for successive shot locations, and all the data were recorded on a single tape recorder dedicated to the experiment. By cross-correlating the records from successive shots, it was possible to determine differential travel times to an accuracy of 4 msec or better using a 250-kg explosive charge recorded in the distance range between 30 to 50 km. No definitive changes in seismic velocity were detected during the time covered by this experiment. The spatial distribution of stations and the temporal distribution of shots were designed to reveal anomalies that might be associated with a magnitude 4 or greater earthquake. Nine small earthquakes (M = 2.9 to 3.3) did occur during the experiment, but anomalies associated with these earthquakes could have escaped detection. This work demonstrates the feasibility of monitoring large regions with a high density of instruments to detect systematic changes in arrival times from artifical sources with an accuracy approaching 1 msec and a cost commensurate with the task of earthquake prediction.

scholarly journals Model-Based Probabilistic Inversion Using Magnetic Data: A Case Study on the Kevitsa Deposit

Geosciences ◽  
2021 ◽  
Vol 11 (4) ◽  
pp. 150
Author(s):  
Nilgün Güdük ◽  
Miguel de la Varga ◽  
Janne Kaukolinna ◽  
Florian Wellmann
Keyword(s):  
Structural Parameters ◽  
Rock Properties ◽  
Magnetic Data ◽  
Geological Model ◽  
Borehole Data ◽  
Likelihood Functions ◽  
Different Types ◽  
Geophysical Information ◽  
Different Sources

Structural geological models are widely used to represent relevant geological interfaces and property distributions in the subsurface. Considering the inherent uncertainty of these models, the non-uniqueness of geophysical inverse problems, and the growing availability of data, there is a need for methods that integrate different types of data consistently and consider the uncertainties quantitatively. Probabilistic inference provides a suitable tool for this purpose. Using a Bayesian framework, geological modeling can be considered as an integral part of the inversion and thereby naturally constrain geophysical inversion procedures. This integration prevents geologically unrealistic results and provides the opportunity to include geological and geophysical information in the inversion. This information can be from different sources and is added to the framework through likelihood functions. We applied this methodology to the structurally complex Kevitsa deposit in Finland. We started with an interpretation-based 3D geological model and defined the uncertainties in our geological model through probability density functions. Airborne magnetic data and geological interpretations of borehole data were used to define geophysical and geological likelihoods, respectively. The geophysical data were linked to the uncertain structural parameters through the rock properties. The result of the inverse problem was an ensemble of realized models. These structural models and their uncertainties are visualized using information entropy, which allows for quantitative analysis. Our results show that with our methodology, we can use well-defined likelihood functions to add meaningful information to our initial model without requiring a computationally-heavy full grid inversion, discrepancies between model and data are spotted more easily, and the complementary strength of different types of data can be integrated into one framework.

scholarly journals Near-surface real-time seismic imaging using parsimonious interferometry

2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Sherif M. Hanafy ◽  
Hussein Hoteit ◽  
Jing Li ◽  
Gerard T. Schuster
Keyword(s):  
Fluid Flow ◽  
Real Time ◽  
Temporal Resolution ◽  
Seismic Imaging ◽  
Time Lapse ◽  
Rock Properties ◽  
Recording Time ◽  
Near Surface ◽  
Arrival Times ◽  
Order Of Magnitude

AbstractResults are presented for real-time seismic imaging of subsurface fluid flow by parsimonious refraction and surface-wave interferometry. Each subsurface velocity image inverted from time-lapse seismic data only requires several minutes of recording time, which is less than the time-scale of the fluid-induced changes in the rock properties. In this sense this is real-time imaging. The images are P-velocity tomograms inverted from the first-arrival times and the S-velocity tomograms inverted from dispersion curves. Compared to conventional seismic imaging, parsimonious interferometry reduces the recording time and increases the temporal resolution of time-lapse seismic images by more than an order-of-magnitude. In our seismic experiment, we recorded 90 sparse data sets over 4.5 h while injecting 12-tons of water into a sand dune. Results show that the percolation of water is mostly along layered boundaries down to a depth of a few meters, which is consistent with our 3D computational fluid flow simulations and laboratory experiments. The significance of parsimonious interferometry is that it provides more than an order-of-magnitude increase of temporal resolution in time-lapse seismic imaging. We believe that real-time seismic imaging will have important applications for non-destructive characterization in environmental, biomedical, and subsurface imaging.

scholarly journals Seismic velocity structure at the southern termination of the 2016 Kumamoto Earthquake rupture, Japan

2020 ◽  
Vol 72 (1) ◽  
Author(s):  
Yasuhira Aoyagi ◽  
Haruo Kimura ◽  
Kazuo Mizoguchi
Keyword(s):  
Fault Zone ◽  
Velocity Structure ◽  
Seismic Velocity ◽  
Earthquake Rupture ◽  
2016 Kumamoto Earthquake ◽  
Seismic Velocity Structure ◽  
Kumamoto Earthquake ◽  
Seismogenic Layer ◽  
The 2016 Kumamoto Earthquake ◽  
Tectonic Line

Abstract The earthquake rupture termination mechanism and size of the ruptured area are crucial parameters for earthquake magnitude estimations and seismic hazard assessments. The 2016 Mw 7.0 Kumamoto Earthquake, central Kyushu, Japan, ruptured a 34-km-long area along previously recognized active faults, eastern part of the Futagawa fault zone and northernmost part of the Hinagu fault zone. Many researchers have suggested that a magma chamber under Aso Volcano terminated the eastward rupture. However, the termination mechanism of the southward rupture has remained unclear. Here, we conduct a local seismic tomographic inversion using a dense temporary seismic network to detail the seismic velocity structure around the southern termination of the rupture. The compressional-wave velocity (Vp) results and compressional- to shear-wave velocity (Vp/Vs) structure indicate several E–W- and ENE–WSW-trending zonal anomalies in the upper to middle crust. These zonal anomalies may reflect regional geological structures that follow the same trends as the Oita–Kumamoto Tectonic Line and Usuki–Yatsushiro Tectonic Line. While the 2016 Kumamoto Earthquake rupture mainly propagated through a low-Vp/Vs area (1.62–1.74) along the Hinagu fault zone, the southern termination of the earthquake at the focal depth of the mainshock is adjacent to a 3-km-diameter high-Vp/Vs body. There is a rapid 5-km step in the depth of the seismogenic layer across the E–W-trending velocity boundary between the low- and high-Vp/Vs areas that corresponds well with the Rokkoku Tectonic Line; this geological boundary is the likely cause of the dislocation of the seismogenic layer because it is intruded by serpentinite veins. A possible factor in the southern rupture termination of the 2016 Kumamoto Earthquake is the existence of a high-Vp/Vs body in the direction of southern rupture propagation. The provided details of this inhomogeneous barrier, which are inferred from the seismic velocity structures, may improve future seismic hazard assessments for a complex fault system composed of multiple segments.

scholarly journals Physics-Based Relationship for Pore Pressure and Vertical Stress Monitoring Using Seismic Velocity Variations

2021 ◽  
Vol 13 (14) ◽  
pp. 2684
Author(s):  
Eldert Fokker ◽  
Elmer Ruigrok ◽  
Rhys Hawkins ◽  
Jeannot Trampert
Keyword(s):  
Pore Pressure ◽  
Seismic Velocity ◽  
Pressure Head ◽  
Groundwater Table ◽  
Surface Velocity ◽  
Seismic Velocities ◽  
Stress Monitoring ◽  
Near Surface ◽  
Velocity Variations ◽  
Velocity Changes

Previous studies examining the relationship between the groundwater table and seismic velocities have been guided by empirical relationships only. Here, we develop a physics-based model relating fluctuations in groundwater table and pore pressure with seismic velocity variations through changes in effective stress. This model justifies the use of seismic velocity variations for monitoring of the pore pressure. Using a subset of the Groningen seismic network, near-surface velocity changes are estimated over a four-year period, using passive image interferometry. The same velocity changes are predicted by applying the newly derived theory to pressure-head recordings. It is demonstrated that the theory provides a close match of the observed seismic velocity changes.

scholarly journals Rock Location and Property Analysis of Lunar Regolith at Chang’E-4 Landing Site Based on Local Correlation and Semblance Analysis

2020 ◽  
Vol 13 (1) ◽  
pp. 48
Author(s):  
Hanjie Song ◽  
Chao Li ◽  
Jinhai Zhang ◽  
Xing Wu ◽  
Yang Liu ◽  
...  
Keyword(s):  
Electromagnetic Waves ◽  
Lunar Regolith ◽  
Landing Site ◽  
Rock Properties ◽  
The Moon ◽  
Local Correlation ◽  
Near Surface ◽  
Property Analysis ◽  
Stratigraphic Division ◽  
Formation And Evolution

The Lunar Penetrating Radar (LPR) onboard the Yutu-2 rover from China’s Chang’E-4 (CE-4) mission is used to probe the subsurface structure and the near-surface stratigraphic structure of the lunar regolith on the farside of the Moon. Structural analysis of regolith could provide abundant information on the formation and evolution of the Moon, in which the rock location and property analysis are the key procedures during the interpretation of LPR data. The subsurface velocity of electromagnetic waves is a vital parameter for stratigraphic division, rock location estimates, and calculating the rock properties in the interpretation of LPR data. In this paper, we propose a procedure that combines the regolith rock extraction technique based on local correlation between the two sets of LPR high-frequency channel data and the common offset semblance analysis to determine the velocity from LPR diffraction hyperbola. We consider the heterogeneity of the regolith and derive the relative permittivity distribution based on the rock extraction and semblance analysis. The numerical simulation results show that the procedure is able to obtain the high-precision position and properties of the rock. Furthermore, we apply this procedure to CE-4 LPR data and obtain preferable estimations of the rock locations and the properties of the lunar subsurface regolith.

Rock mass structure analysis based on seismic velocity and attenuation images

2000 ◽  
Vol 45 (13) ◽  
pp. 1211-1216 ◽  
Author(s):  
Xu Chang ◽  
Yike Liu ◽  
Hui Wang ◽  
Xing Gao
Keyword(s):  
Rock Mass ◽  
Structure Analysis ◽  
Seismic Velocity ◽  
Rock Mass Structure

Channel waves in cross‐borehole data

Geophysics ◽  
1992 ◽  
Vol 57 (2) ◽  
pp. 334-342 ◽  
Author(s):  
Larry R. Lines ◽  
Kenneth R. Kelly ◽  
John Queen
Keyword(s):  
Seismic Velocity ◽  
Test Facility ◽  
Equation Modeling ◽  
Borehole Data ◽  
Geological Information ◽  
Channel Wave ◽  
First Arrival ◽  
Borehole Seismic ◽  
Seismic Velocity Contrasts ◽  
Ray Bending

Layered geological formations with large seismic velocity contrasts can effectively create channel waves in cross‐borehole seismic data. The existence of channel waves for such waveguides can be confirmed by ray tracing, wave equation modeling, and modal analysis. Channel wave arrivals are identified in cross‐borehole data recorded at Conoco’s Newkirk test facility. For these data, where velocity contrasts are about 2 to 1, tomography based on first arrival traveltimes, is limited due to problems with extreme ray bending and seismic shadow zones. However, it may be possible to extract geological information using channel wave information. The seismometer differencing method appears to be a promising approach for detecting waveguide boundaries by use of cross‐borehole data.

Sign in / Sign up

Export Citation Format

Share Document