Seismic wave velocity investigation at The Geysers‐Clear Lake geothermal field, California

Geophysics ◽  
1982 ◽  
Vol 47 (5) ◽  
pp. 819-824 ◽  
Author(s):  
Harsh K. Gupta ◽  
Ronald W. Ward ◽  
Tzeu‐Lie Lin

Analysis of P‐ and S‐waves from shallow microearthquakes in the vicinity of The Geysers geothermal area, California, recorded by a dense, telemetered seismic array operated by the U.S. Geological Survey (USGS) shows that these phases are easily recognized and traced on record sections to distances of 80 km. Regional average velocities for the upper crust are estimated to be [Formula: see text] and [Formula: see text] for P‐ and S‐waves, respectively. Poisson’s ratio is estimated at 23 locations using Wadati diagrams and is found to vary from 0.13 to 0.32. In general, the Poisson’s ratio is found to be lower at the locations close to the steam production zones at The Geysers and Clear Lake volcanic field to the northeast. The low Poisson ratio corresponds to a decrease in P‐wave velocity in areas of high heat flow. The decrease may be caused by fracturing of the rock and saturation with gas or steam.

Geofluids ◽  
2021 ◽  
Vol 2021 ◽  
pp. 1-22
Author(s):  
Yunjiang Sun ◽  
Jianping Zuo ◽  
Yue Shi ◽  
Zhengdai Li ◽  
Changning Mi ◽  
...  

Ultrasonic wave velocity is effective to evaluate anisotropy property and predict rock failure. This paper investigates the correlation between dynamic ultrasonic and mechanical properties of sandstones with different buried depths subjected to uniaxial compression tests. The circumferential anisotropy and axial wave velocity of sandstone are obtained by means of ultrasonic wave velocity measurements. The mechanical properties, including Young’s modulus and uniaxial compressive strength, are positively correlated with the axial P wave velocity. The average angles between the sandstone failure plane and the minimum and maximum wave directions are 35.8° and 63.3°, respectively. The axial P wave velocity almost keeps constant, and the axial S wave velocity has a decreasing trend before the failure of rock specimen. In most rock samples under uniaxial compression, shear failure occurs in the middle and splitting appears near both sides. Additionally, the dynamic Young’s modulus and dynamic Poisson’s ratio during loading are obtained, and the negative values of the Poisson’s ratio occur at the initial compression stage. Distortion and rotation of micro/mesorock structures may be responsible for the negative Poisson’s ratio.


2020 ◽  
Vol 222 (1) ◽  
pp. 582-594
Author(s):  
Thomas Forbriger ◽  
Lingli Gao ◽  
Peter Malischewsky ◽  
Matthias Ohrnberger ◽  
Yudi Pan

SUMMARY Other than commonly assumed in seismology, the phase velocity of Rayleigh waves is not necessarily a single-valued function of frequency. In fact, a single Rayleigh mode can exist with three different values of phase velocity at one frequency. We demonstrate this for the first higher mode on a realistic shallow seismic structure of a homogeneous layer of unconsolidated sediments on top of a half-space of solid rock (LOH). In the case of LOH a significant contrast to the half-space is required to produce the phenomenon. In a simpler structure of a homogeneous layer with fixed (rigid) bottom (LFB) the phenomenon exists for values of Poisson’s ratio between 0.19 and 0.5 and is most pronounced for P-wave velocity being three times S-wave velocity (Poisson’s ratio of 0.4375). A pavement-like structure (PAV) of two layers on top of a half-space produces the multivaluedness for the fundamental mode. Programs for the computation of synthetic dispersion curves are prone to trouble in such cases. Many of them use mode-follower algorithms which loose track of the dispersion curve and miss the multivalued section. We show results for well established programs. Their inability to properly handle these cases might be one reason why the phenomenon of multivaluedness went unnoticed in seismological Rayleigh wave research for so long. For the very same reason methods of dispersion analysis must fail if they imply wave number kl(ω) for the lth Rayleigh mode to be a single-valued function of frequency ω. This applies in particular to deconvolution methods like phase-matched filters. We demonstrate that a slant-stack analysis fails in the multivalued section, while a Fourier–Bessel transformation captures the complete Rayleigh-wave signal. Waves of finite bandwidth in the multivalued section propagate with positive group-velocity and negative phase-velocity. Their eigenfunctions appear conventional and contain no conspicuous feature.


Author(s):  
M J Funnell ◽  
A H Robinson ◽  
R W Hobbs ◽  
C Peirce

Summary The seismic velocity of the oceanic crust is a function of its physical properties that include its lithology, degree of alteration, and porosity. Variations in these properties are particularly significant in young crust, but also occur with age as it evolves through hydrothermal circulation and is progressively covered with sediment. While such variation may be investigated through P-wave velocity alone, joint analysis with S-wave velocity allows the determination of Poisson's ratio, which provides a more robust insight into the nature of change in these properties. Here we describe the independent modelling of P- and S-wave seismic datasets, acquired along an ∼330 km-long profile traversing new to ∼8 Myr-old oceanic crust formed at the intermediate-spreading Costa Rica Rift (CRR). Despite S-wave data coverage being almost four-times lower than that of the P-wave dataset, both velocity models demonstrate correlations in local variability and a long-wavelength increase in velocity with distance, and thus age, from the ridge axis of up to 0.8 and 0.6 km s−1, respectively. Using the Vp and Vs models to calculate Poisson's ratio (σ), it reveals a typical structure for young oceanic crust, with generally high values in the uppermost crust that decrease to a minimum of 0.24 by 1.0–1.5 km sub-basement, before increasing again throughout the lower crust. The observed upper crustal decrease in σ most likely results from sealing of fractures, which is supported by observations of a significant decrease in porosity with depth (from ∼15 to < 2 per cent) through the dyke sequence in Ocean Drilling Program borehole 504B. High Poisson's ratio (>0.31) is observed throughout the crust of the north flank of the CRR axis and, whilst this falls within the ‘serpentinite’ classification of lithological proxies, morphological evidence of pervasive surface magmatism and limited tectonism suggests, instead, that the cause is porosity in the form of pervasive fracturing and, thus, that this is the dominant control on seismic velocity in the newly formed CRR crust. South of the CRR, the values of Poisson's ratio are representative of more typical oceanic crust, and decrease with increasing distance from the spreading centre, most likely as a result of mineralisation and increased fracture infill. This is supported by borehole observations and modelled 3-D seismic anisotropy. Crustal segments formed during periods of particularly low half-spreading rate (<35 mm yr−1) demonstrate high Poisson's ratio relative to the background, indicating the likely retention of increased porosity and fracturing associated with the greater degrees of tectonism at the time of their formation. Across the south flank of the CRR, we find that the average Poisson's ratio in the upper 1 km of the crust decreases with age by ∼0.0084 Myr−1 prior to the thermal sealing of the crust, suggesting that, to at least ∼7 Myr, advective hydrothermal processes dominate early CRR-generated oceanic crustal evolution, consistent with heat flow measurements.


Geophysics ◽  
1995 ◽  
Vol 60 (3) ◽  
pp. 651-659 ◽  
Author(s):  
Mark E. Mathisen ◽  
Paul Cunningham ◽  
Jesse Shaw ◽  
Anthony A. Vasiliou ◽  
J. H. Justice ◽  
...  

S‐wave, P‐wave, and Poisson’s ratio tomograms have been used to interpret the 3-D distribution of rock and fluid properties during an early phase of a California heavy oil sand steamflood. Four lines of good quality crosswell seismic data, with source to receiver offsets ranging from 287 to 551 ft (87 to 168 m), were acquired in a radial pattern around a high temperature cemented receiver cable in four days. Processing, first‐arrival picking, and good quality tomographic reconstructions were completed despite offset‐related variations in data quality between the long and short lines. Interpretation was based on correlations with reservoir models, log, core, temperature, and steam injection data. S‐wave tomograms define the 3-D distribution of the “high flow” turbidite channel facies, the “moderate‐low flow” levee facies, porosity, and structural dip. The S‐wave tomograms also define an area with anomalously low S‐wave velocity, which correlates with low shear log velocities and suggests that pressure‐related dilation and compaction may be imageable. P‐wave tomograms define the same reservoir lithology and structure as the S‐wave tomograms and the 3-D distribution of low compressional velocity zones formed by previous steam‐heat injection and the formation of gas. The low P‐wave velocity zones, which are laterally continuous in the “high flow” channel facies near the top of most zones, indicate that the steam‐heat‐gas distribution is controlled by stratification. The stratigraphic control of gas‐bearing zones inferred from P‐wave tomograms is confirmed by Poisson’s ratio tomograms which display low Poisson’s ratios indicative of gas (<0.35) in the same zones as the low P‐wave velocities. The interpretation results indicate that radial survey tomograms can be tied at a central well and used to develop an integrated 3-D geoscience‐engineering reservoir model despite offset‐related variations in data quality. The laterally continuous, stratification‐controlled, low P‐wave velocity zones, which extend up‐dip, suggest that significant amounts of steam‐heat are not heating the surrounding reservoir volume but are flowing updip along “high flow” channels.


Geophysics ◽  
2021 ◽  
pp. 1-56
Author(s):  
Farbod Khosro Anjom ◽  
Thomas Jules Browaeys ◽  
Laura Valentina Socco

Exploration seismic surveys in hard-to-access areas such as foothills and forests are extremely challenging. The Multiphysics Exploration Technologies Integrated System (METIS) research project was initiated to design an exploration system, facilitating the acquisition in these areas by delivering the receivers from the sky using unmanned aerial vehicles. Air dropping of the sensors in vegetated areas results in an irregular geometry for the acquisition. This irregularity can limit the application of conventional surface wave methods. We have developed a surface wave workflow for estimating the S-wave velocity ( VS) and P-wave velocity ( VP) models and that supports the irregular geometry of the deployed sources and receivers. The method consists of a multimodal surface wave tomography (SWT) technique to compute the VS model and a data transform method (the wavelength/depth [W/D] method) to determine the Poisson’s ratio and VP model. We applied the method to the METIS’s first pilot records, which were acquired in the forest of Papua New Guinea. Application of SWT to the data resulted in the first 90 m of the VS model. The W/D method provided the Poisson’s ratio averaged over the area and the VP model between 10 and 70 m from the surface. The impact of the acquisition scale and layout on the resolution of the estimated model and the advantages of including the higher modes of surface waves in the tomographic inversion are assessed in detail. The presence of shots from diverse site locations significantly improves the resolution of the obtained model. Including the higher modes enhances the data coverage and increases the investigation depth.


Sign in / Sign up

Export Citation Format

Share Document