3-D microearthquake attenuation tomography at the Northwest Geysers geothermal region, California

Geophysics ◽  
1997 ◽  
Vol 62 (1) ◽  
pp. 149-167 ◽  
Author(s):  
Arturo E. Romero ◽  
Thomas V. McEvilly ◽  
Ernest L. Majer
Keyword(s):  
Geothermal Field ◽  
Text Structure ◽  
P Wave ◽  
Inversion Method ◽  
Reference Spectrum ◽  
Spectral Ratios ◽  
Velocity Inversion ◽  
Velocity Models ◽  
Average Spectrum ◽  
Attenuation Structure

We present the results of the modeling of a 3-D differential attenuation structure [Formula: see text] beneath the Northwest Geysers geothermal field. A set of 480 high‐quality microearthquakes distributed evenly throughout the field and with a minimum of 10 P‐wave arrivals was selected for the study. We constructed spectral ratios by dividing each spectrum with a reference spectrum from each observing station. The reference spectrum was derived from the average spectrum of all events that were recorded at that station to correct for the strong site dependence of the observed spectra. We then estimated the differential attenuation operator from the slopes of the spectral ratios. The velocity models and the raypaths for all events are known from a previous velocity inversion study. The inversion for the differential attenuation structure was carried out using a modified progressive inversion method. The observed [Formula: see text] structure correlates well with mapped geologic units. High [Formula: see text] and lower velocities correlate with Franciscan melange, while lower attenuation and higher velocities correspond to metagraywacke units. High P‐wave [Formula: see text] also underlies the southern region between 2 and 3 km depth where low [Formula: see text] values suggest undersaturation of the reservoir rocks. Most of the steam entries also occur within this region and probably delineate the steam reservoir. These anomalies may be explained by high rock temperatures and the presence of steam and other gases.

The relationship between Vs, Vp, density and depth based on PS-logging data at K-NET and KiK-net sites

2021 ◽  
Author(s):  
Fumiaki Nagashima ◽  
Hiroshi Kawase
Keyword(s):  
Standard Deviation ◽  
Seismic Velocity ◽  
Saturated Soil ◽  
P Wave ◽  
Soil Types ◽  
Depth Range ◽  
Velocity Inversion ◽  
Velocity Models ◽  
Regression Curves ◽  
Ps Logging

Summary P-wave velocity (Vp) is an important parameter for constructing seismic velocity models of the subsurface structures by using microtremors and earthquake ground motions or any other geophysical exploration data. In order to reflect the ground survey information in Japan to the Vp structure, we investigated the relationships among Vs, Vp, and depth by using PS-logging data at all K-NET and KiK-net sites. Vp values are concentrated at around 500 m/s and 1,500 m/s when Vs is lower than 1,000 m/s, where these concentrated areas show two distinctive characteristics of unsaturated and saturated soil, respectively. Many Vp values in the layer shallower than 4 m are around 500 m/s, which suggests the dominance of unsaturated soil, while many Vp values in the layer deeper than 4 m are larger than 1,500 m/s, which suggests the dominance of saturated soil there. We also investigated those relationships for different soil types at K-NET sites. Although each soil type has its own depth range, all soil types show similar relationships among Vs, Vp, and depth. Then, considering the depth profile of Vp, we divided the dataset into two by the depth, which is shallower or deeper than 4 m, and calculated the geometrical mean of Vp and the geometrical standard deviation in every Vs bins of 200 m/s. Finally, we obtained the regression curves for the average and standard deviation of Vp estimated from Vs to get the Vp conversion functions from Vs, which can be applied to a wide Vs range. We also obtained the regression curves for two datasets with Vp lower and higher than 1,200 m/s. These regression curves can be applied when the groundwater level is known. In addition, we obtained the regression curves for density from Vs or Vp. An example of the application for those relationships in the velocity inversion is shown.

Crustal velocity structure beneath the NW Dinarides from 1-D hypocenter-velocity inversion

2021 ◽  
Author(s):  
Gregor Rajh ◽  
Josip Stipčević ◽  
Mladen Živčić ◽  
Marijan Herak ◽  
Andrej Gosar
Keyword(s):  
Velocity Structure ◽  
Velocity Model ◽  
P Wave ◽  
Depth Range ◽  
Velocity Inversion ◽  
Arrival Times ◽  
Simultaneous Inversion ◽  
Velocity Modeling ◽  
Velocity Models ◽  
Station Corrections

<p>The investigated area of the NW Dinarides is bordered by the Adriatic foreland, the Southern Alps, and the Pannonian basin at the NE corner of the Adriatic Sea. Its complex crustal structure is the result of interactions among different tectonic units. Despite numerous seismic studies taking place in this region, there still exists a need for a detailed, smaller scale study focusing mainly on the brittle part of the Earth's crust. Therefore, we decided to investigate the velocity structure of the crust using concepts of local earthquake tomography (LET) and minimum 1-D velocity model. Here, we present the results of the 1-D velocity modeling and the catalogue of the relocated seismicity. A minimum 1-D velocity model is computed by simultaneous inversion for hypocentral and velocity parameters together with seismic station corrections and represents the best fit to the observed arrival times.</p><p>We used 15,579 routinely picked P wave arrival times from 631 well-located earthquakes that occurred in Slovenia and in its immediate surroundings (mainly NW Croatia). Various initial 1-D velocity models, differing in velocity and layering, were used as input for velocity inversion in the VELEST program. We also varied several inversion parameters during the inversion runs. Most of the computed 1-D velocity models converged to a stable solution in the depth range between 0 and 25 km. We evaluated the inversion results using rigorous testing procedures and selected two best performing velocity models. Each of these models will be used independently as the initial model in the simultaneous hypocenter-velocity inversion for a 3-D velocity structure in LET. Based on the results of the 1-D velocity modeling, seismicity distribution, and tectonics, we divided the study area into three parts, redefined the earthquake-station geometry, and performed the inversion for each part separately. This way, we gained a better insight into the shallow velocity structure of each subregion and were able to demonstrate the differences among them.</p><p>Besides general structural implications and a potential to improve the results of LET, the new 1-D velocity models along with station corrections can also be used in fast routine earthquake location and to detect systematic travel time errors in seismological bulletins, as already shown by some studies using similar methods.</p>

Joint inversion of Rayleigh-wave dispersion and P-wave refraction data for laterally varying layered models

Geophysics ◽  
2014 ◽  
Vol 79 (4) ◽  
pp. EN49-EN59 ◽  
Author(s):  
Daniele Boiero ◽  
Laura Valentina Socco
Keyword(s):  
Rayleigh Wave ◽  
Joint Inversion ◽  
A Priori ◽  
P Wave ◽  
Inversion Method ◽  
Model Parameters ◽  
Data Set ◽  
Wave Refraction ◽  
Velocity Models ◽  
Refraction Data

We implemented a joint inversion method to build P- and S-wave velocity models from Rayleigh-wave and P-wave refraction data, specifically designed to deal with laterally varying layered environments. A priori information available over the site and any physical law to link model parameters can be also incorporated. We tested and applied the algorithm behind the method. The results from a field data set revealed advantages with respect to individual surface-wave analysis (SWA) and body wave tomography (BWT). The algorithm imposed internal consistency for all the model parameters relaxing the required a priori assumptions (i.e., Poisson’s ratio level of confidence in SWA) and the inherent limitations of the two methods (i.e., velocity decreases for BWT).

First-arrival traveltime inversion of seismic diving waves observed on undulant surface

2021 ◽  
Vol 225 (2) ◽  
pp. 1020-1031
Author(s):  
Huachen Yang ◽  
Jianzhong Zhang ◽  
Kai Ren ◽  
Changbo Wang
Keyword(s):  
Turning Points ◽  
Analytical Formula ◽  
Velocity Model ◽  
Western China ◽  
Inversion Method ◽  
Mountain Area ◽  
Traveltime Inversion ◽  
Static Correction ◽  
Velocity Models ◽  
First Arrival

SUMMARY A non-iterative first-arrival traveltime inversion method (NFTI) is proposed for building smooth velocity models using seismic diving waves observed on irregular surface. The new ray and traveltime equations of diving waves propagating in smooth media with undulant observation surface are deduced. According to the proposed ray and traveltime equations, an analytical formula for determining the location of the diving-wave turning points is then derived. Taking the influence of rough topography on first-arrival traveltimes into account, the new equations for calculating the velocities at turning points are established. Based on these equations, a method is proposed to construct subsurface velocity models from the observation surface downward to the bottom using the first-arrival traveltimes in common offset gathers. Tests on smooth velocity models with rugged topography verify the validity of the established equations, and the superiority of the proposed NFTI. The limitation of the proposed method is shown by an abruptly-varying velocity model example. Finally, the NFTI is applied to solve the static correction problem of the field seismic data acquired in a mountain area in the western China. The results confirm the effectivity of the proposed NFTI.

scholarly journals Seismic Structure of Local Crustal Earthquakes beneath the Zipingpu Reservoir of Longmenshan Fault Zone

2011 ◽  
Vol 2011 ◽  
pp. 1-6 ◽  
Author(s):  
Haiou Li ◽  
Xiwei Xu ◽  
Wentao Ma ◽  
Ronghua Xie ◽  
Jingli Yuan ◽  
...  
Keyword(s):  
Fault Zone ◽  
Three Dimensional ◽  
P Wave ◽  
The Tibetan Plateau ◽  
Seismic Structure ◽  
Infiltration Depth ◽  
Low Velocity ◽  
Longmenshan Fault Zone ◽  
Velocity Models ◽  
Longmenshan Fault

Three-dimensional P wave velocity models under the Zipingpu reservoir in Longmenshan fault zone are obtained with a resolution of 2 km in the horizontal direction and 1 km in depth. We used a total of 8589 P wave arrival times from 1014 local earthquakes recorded by both the Zipingpu reservoir network and temporary stations deployed in the area. The 3-D velocity images at shallow depth show the low-velocity regions have strong correlation with the surface trace of the Zipingpu reservoir. According to the extension of those low-velocity regions, the infiltration depth directly from the Zipingpu reservoir itself is limited to 3.5 km depth, while the infiltration depth downwards along the Beichuan-Yingxiu fault in the study area is about 5.5 km depth. Results show the low-velocity region in the east part of the study area is related to the Proterozoic sedimentary rocks. The Guanxian-Anxian fault is well delineated by obvious velocity contrast and may mark the border between the Tibetan Plateau in the west and the Sichuan basin in the east.

Fracture detection using P-wave and S-wave vertical seismic profiling at The Geysers

Geophysics ◽  
1988 ◽  
Vol 53 (1) ◽  
pp. 76-84 ◽  
Author(s):  
E. L. Majer ◽  
T. V. McEvilly ◽  
F. S. Eastwood ◽  
L. R. Myer
Keyword(s):  
Fractured Rock ◽  
Geothermal Field ◽  
P Wave ◽  
Velocity Variation ◽  
Fracture Detection ◽  
S Wave ◽  
Fracture Geometry ◽  
Seismic Profiling ◽  
Vertical Seismic Profiling ◽  
The Geysers

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.

A method for correcting acoustic finite-difference amplitudes for elastic effects

Geophysics ◽  
2014 ◽  
Vol 79 (4) ◽  
pp. T243-T255 ◽  
Author(s):  
James W. D. Hobro ◽  
Chris H. Chapman ◽  
Johan O. A. Robertsson
Keyword(s):  
Wave Equation ◽  
Finite Difference ◽  
Effective Means ◽  
Cost Effective ◽  
P Wave ◽  
S Waves ◽  
Velocity Models ◽  
Acoustic Simulation ◽  
Wide Range ◽  
Elastic Simulation

We present a new method for correcting the amplitudes of arrivals in an acoustic finite-difference simulation for elastic effects. In this method, we selectively compute an estimate of the error incurred when the acoustic wave equation is used to approximate the behavior of the elastic wave equation. This error estimate is used to generate an effective source field in a second acoustic simulation. The result of this second simulation is then applied as a correction to the original acoustic simulation. The overall cost is approximately twice that of an acoustic simulation but substantially less than the cost of an elastic simulation. Because both simulations are acoustic, no S-waves are generated, so dispersed converted waves are avoided. We tested the characteristics of the method on a simple synthetic model designed to simulate propagation through a strong acoustic impedance contrast representative of sedimentary geology. It corrected amplitudes to high accuracy for reflected arrivals over a wide range of incidence angles. We also evaluated results from simulations on more complex models that demonstrated that the method was applicable in realistic sedimentary models containing a wide range of seismic contrasts. However, its accuracy was reduced for wide-angle reflections from very high impedance contrasts such as a shallow top-salt interface. We examined the influence of modeling at coarse grid resolutions, in which converted S-waves in the equivalent elastic simulation are dispersed. These results provide some validation for the accuracy of the method when applied using finite-difference grids designed for acoustic modeling. The method appears to offer a cost-effective means of modeling elastic amplitudes for P-wave arrivals in a useful range of velocity models. It has several potential applications in imaging and inversion.

scholarly journals Fluid activity detection in geothermal areas using a single seismic station by monitoring horizontal-to-vertical spectral ratios

2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Kyosuke Okamoto ◽  
Hiroshi Asanuma ◽  
Hiro Nimiya
Keyword(s):  
Carbon Capture ◽  
Seismic Station ◽  
Carbon Capture And Storage ◽  
Fluid Flows ◽  
Temporal Variations ◽  
Geothermal Field ◽  
Spectral Ratios ◽  
Subsurface Structures ◽  
Geothermal Fields ◽  
Single Station

AbstractSubsurface structure survey based on horizontal-to-vertical (H/V) spectral ratios is widely conducted. The major merit of this survey is its convenience to obtain a stable result using a single station. Spatial variations of H/V spectral ratios are well-known phenomena, and it has been used to estimate the spatial fluctuation in subsurface structures. It is reasonable to anticipate temporal variations in H/V spectral ratios, especially in areas like geothermal fields, carbon capture and storage fields, etc., where rich fluid flows are expected, although there are few reports about the temporal changes. In Okuaizu Geothermal Field (OGF), Japan, dense seismic monitoring was deployed in 2015, and continuous monitoring has been consistent. We observed the H/V spectral ratios in OGF and found their repeated temporary drops. These drops seemed to be derived from local fluid activities according to a numerical calculation. Based on this finding, we examined a coherency between the H/V spectral ratios and fluid activities in OGF and found a significance. In conclusion, monitoring H/V spectral ratios can enable us to grasp fluid activities that sometimes could lead to a relatively large seismic event.

scholarly journals The offset-midpoint traveltime pyramid in 3D transversely isotropic media with a horizontal symmetry axis

Geophysics ◽  
2015 ◽  
Vol 80 (1) ◽  
pp. T51-T62 ◽  
Author(s):  
Qi Hao ◽  
Alexey Stovas ◽  
Tariq Alkhalifah
Keyword(s):  
Symmetry Axis ◽  
Transversely Isotropic ◽  
P Wave ◽  
Analytic Representation ◽  
Velocity Inversion ◽  
Transversely Isotropic Media ◽  
Horizontal Symmetry ◽  
Time Domains ◽  
Isotropic Media ◽  
Hti Media

Analytic representation of the offset-midpoint traveltime equation for anisotropy is very important for prestack Kirchhoff migration and velocity inversion in anisotropic media. For transversely isotropic media with a vertical symmetry axis, the offset-midpoint traveltime resembles the shape of a Cheops’ pyramid. This is also valid for homogeneous 3D transversely isotropic media with a horizontal symmetry axis (HTI). We extended the offset-midpoint traveltime pyramid to the case of homogeneous 3D HTI. Under the assumption of weak anellipticity of HTI media, we derived an analytic representation of the P-wave traveltime equation and used Shanks transformation to improve the accuracy of horizontal and vertical slownesses. The traveltime pyramid was derived in the depth and time domains. Numerical examples confirmed the accuracy of the proposed approximation for the traveltime function in 3D HTI media.

scholarly journals PRELIMINARY GULLY HAZARD ASSESSMENT USING 2D-ERT AND 2D-SRT GEOPHYSICAL SURVEYS: A CASE STUDY IN SOUTHEASTERN NIGERIA

2021 ◽  
pp. 49-65
Author(s):  
GN Egwuonwu ◽  
EI Okoyeh ◽  
DC Agarana ◽  
EG Nwaka ◽  
OB Nwosu ◽  
...  
Keyword(s):  
Seismic Refraction ◽  
P Wave ◽  
Lateritic Soil ◽  
Risk Level ◽  
Clayey Sand ◽  
Southeastern Nigeria ◽  
Velocity Models ◽  
Refraction Tomography ◽  
Geophysical Surveys ◽  
Oriented Parallel

Two-dimensional Electrical Resistivity Tomography (2DERT) and Seismic Refraction Tomography (2DSRT) were concurrently applied in assessment of a gully site with the view of assessing its stability and risk level. Eight profile lines oriented parallel and perpendicular to the boundary of the gully were surveyed. As a result, apparent resistivity model tomograms in the range of 1-9,000 and p-wave velocity models in the range of 300-700 were obtained from the two techniques respectively. Interpretation of the models obtained show predominance of unconsolidated clay, shale intercalates, clayey sand, sandy clay and weathered lateritic soil at shallow depths. Low amplitude undulating refracting layers, landslide slip subsurface and lose horizons were also delineated at shallow depths. The predominance of weak, clayey and unconsolidated lithology at the gully site suggests evidence of unstable gravitational equilibrium which imply environmental hazard. The plausible deductions made from the two

Sign in / Sign up

Export Citation Format

Share Document