Passive seismic imaging using microearthquakes

Geophysics ◽  
1995 ◽  
Vol 60 (4) ◽  
pp. 1178-1186 ◽  
Author(s):  
M. Reza Daneshvar ◽  
Clarence S. Clay ◽  
Martha K. Savage

We have developed a method of processing seismic signals generated by microearthquakes to image local subsurface structure beneath a recording station. This technique uses the autocorrelation of the vertically traveling earthquake signals to generate pseudoreflection seismograms that can be interpreted for subsurface structure. Processed pseudoreflection data, from microearthquakes recorded in the island of Hawaii, show consistent reflectivity patterns that are interpreted as near‐surface horizontal features. Forward modeling of the pseudoreflection data results in a P‐wave velocity model that shows reasonable agreement with the velocity model derived from a refraction study in the region. Usable signal‐to‐noise ratio is obtained down to 2 s. A shear‐wave velocity model was also generated by applying this technique to horizontal component data.

1998 ◽  
Vol 41 (4) ◽  
Author(s):  
G. Iannaccone ◽  
L. Improta ◽  
P. Capuano ◽  
A. Zollo ◽  
G. Biella ◽  
...  

This paper describes the results of a seismic refraction profile conducted in October 1992 in the Sannio region, Southern Italy, to obtain a detailed P-wave velocity model of the upper crust. The profile, 75 km long, extended parallel to the Apenninic chain in a region frequently damaged in historical time by strong earthquakes. Six shots were fired at five sites and recorded by a number of seismic stations ranging from 41 to 71 with a spacing of 1-2 km along the recording line. We used a two-dimensional raytracing technique to model travel times and amplitudes of first and second arrivals. The obtained P-wave velocity model has a shallow structure with strong lateral variations in the southern portion of the profile. Near surface sediments of the Tertiary age are characterized by seismic velocities in the 3.0-4.1 km/s range. In the northern part of the profile these deposits overlie a layer with a velocity of 4.8 km/s that has been interpreted as a Mesozoic sedimentary succession. A high velocity body, corresponding to the limestones of the Western Carbonate Platform with a velocity of 6 km/s, characterizes the southernmost part of the profile at shallow depths. At a depth of about 4 km the model becomes laterally homogeneous showing a continuous layer with a thickness in the 3-4 km range and a velocity of 6 km/s corresponding to the Meso-Cenozoic limestone succession of the Apulia Carbonate Platform. This platform appears to be layered, as indicated by an increase in seismic velocity from 6 to 6.7 km/s at depths in the 6-8 km range, that has been interpreted as a lithological transition from limestones to Triassic dolomites and anhydrites of the Burano formation. A lower P-wave velocity of about 5.0-5.5 km/s is hypothesized at the bottom of the Apulia Platform at depths ranging from 10 km down to 12.5 km; these low velocities could be related to Permo-Triassic siliciclastic deposits of the Verrucano sequence drilled at the bottom of the Apulia Platform in the Apulia Foreland.


Geophysics ◽  
2004 ◽  
Vol 69 (2) ◽  
pp. 345-351 ◽  
Author(s):  
Geoff J.M. Moret ◽  
William P. Clement ◽  
Michael D. Knoll ◽  
Warren Barrash

P‐wave velocity information obtained from vertical seismic profiles (VSPs) can be useful in imaging subsurface structure, either by directly detecting changes in the subsurface or as an aid to the interpretation of seismic reflection data. In the shallow subsurface, P‐wave velocity can change by nearly an order of magnitude over a short distance, so curved rays are needed to accurately model VSP traveltimes. We used a curved‐ray inversion to estimate the velocity profile and the discrepancy principle to estimate the data noise level and to choose the optimum regularization parameter. The curved‐ray routine performed better than a straight‐ray inversion for synthetic models containing high‐velocity contrasts. The application of the inversion to field data produced a velocity model that agreed well with prior information. These results show that curved‐ray inversion should be used to obtain velocity information from VSPs in the shallow subsurface.


1984 ◽  
Vol 74 (4) ◽  
pp. 1263-1274
Author(s):  
Lawrence H. Jaksha ◽  
David H. Evans

Abstract A velocity model of the crust in northwestern New Mexico has been constructed from an interpretation of direct, refracted, and reflected seismic waves. The model suggests a sedimentary section about 3 km thick with an average P-wave velocity of 3.6 km/sec. The crystalline upper crust is 28 km thick and has a P-wave velocity of 6.1 km/sec. The lower crust below the Conrad discontinuity has an average P-wave velocity of about 7.0 km/sec and a thickness near 17 km. Some evidence suggests that velocity in both the upper and lower crust increases with depth. The P-wave velocity in the uppermost mantle is 7.95 ± 0.15 km/sec. The total crustal thickness near Farmington, New Mexico, is about 48 km (datum = 1.6 km above sea level), and there is evidence for crustal thinning to the southeast.


Geophysics ◽  
2021 ◽  
pp. 1-52
Author(s):  
Yuzhu Liu ◽  
Xinquan Huang ◽  
Jizhong Yang ◽  
Xueyi Liu ◽  
Bin Li ◽  
...  

Thin sand-mud-coal interbedded layers and multiples caused by shallow water pose great challenges to conventional 3D multi-channel seismic techniques used to detect the deeply buried reservoirs in the Qiuyue field. In 2017, a dense ocean-bottom seismometer (OBS) acquisition program acquired a four-component dataset in East China Sea. To delineate the deep reservoir structures in the Qiuyue field, we applied a full-waveform inversion (FWI) workflow to this dense four-component OBS dataset. After preprocessing, including receiver geometry correction, moveout correction, component rotation, and energy transformation from 3D to 2D, a preconditioned first-arrival traveltime tomography based on an improved scattering integral algorithm is applied to construct an initial P-wave velocity model. To eliminate the influence of the wavelet estimation process, a convolutional-wavefield-based objective function for the preprocessed hydrophone component is used during acoustic FWI. By inverting the waveforms associated with early arrivals, a relatively high-resolution underground P-wave velocity model is obtained, with updates at 2.0 km and 4.7 km depth. Initial S-wave velocity and density models are then constructed based on their prior relationships to the P-wave velocity, accompanied by a reciprocal source-independent elastic full-waveform inversion to refine both velocity models. Compared to a traditional workflow, guided by stacking velocity analysis or migration velocity analysis, and using only the pressure component or other single-component, the workflow presented in this study represents a good approach for inverting the four-component OBS dataset to characterize sub-seafloor velocity structures.


2019 ◽  
Vol 24 (1) ◽  
pp. 151-158
Author(s):  
Ionelia Panea

Results are presented for shallow seismic reflection measurements performed southwest of Săcel village in Romania for the purpose of obtaining information about the geological structure in the near subsurface. The P-wave and S-wave velocity distributions were also obtained below the soil surface. The measurements were performed along a nearly linear profile on the top of an elongated hill. Most of the shot gathers were characterized by a good signal-to-noise ratio. A depth-converted migrated section was obtained after the processing of shot gathers, on which an image of sedimentary deposits with various thicknesses, separated by shallow faults until a depth of about 80 m, were observed. The P-wave and S-wave velocity-depth models for two segments were of considerable interest for a geotechnical study proposed for the construction of a windmill park. The two- and three-layered P-wave velocity-depth models were comparable until depths of about 10 m after first-arrival traveltime inversions. The lateral variations in the subsurface geological structure and lithology reflected the variations in the P-wave velocity values from both models. The S-wave velocity-depth models for comparable depth intervals were similar to those from the P-wave velocity-depth models. Reliable S-wave velocity distributions were obtained after inversion of fundamental-mode and higher-mode surface waves.


2019 ◽  
Vol 92 ◽  
pp. 18006
Author(s):  
Yannick Choy Hing Ng ◽  
William Danovan ◽  
Taeseo Ku

Seismic cross-hole tomography has been commonly used in oil and gas exploration and the mining industry for the detection of precious resources. For near-surface geotechnical site investigation, this geophysical method is relatively new and can be used to supplement traditional methods such as the standard penetration test, coring and sampling, thus improving the effectiveness of site characterization. This paper presents a case study which was carried out on a reclaimed land in the Eastern region of Singapore. A seismic cross-hole test was performed by generating both compressional waves and shear waves into the ground. The signals were interpreted by using first-arrival travel time wave tomography and the arrival times were subsequently inverted using Simultaneous Iterative Reconstruction Technique (SIRT). A comparison with the borehole logging data indicated that P-wave velocity model cannot provide sufficient information about the soil layers, especially when the ground water table is near the surface. The S-wave velocity model seemed to agree quite well with the variation in the SPT-N value and could identify to a certain extent the interface between the different soil layers. Finally, P-wave and S-wave velocities are used to compute the Poisson's ratio distribution which gave a good indication of the degree of saturation of the soil.


2019 ◽  
Vol 133 ◽  
pp. 01011
Author(s):  
Jakub Kokowski ◽  
Zbigniew Szreder ◽  
Elżbieta Pilecka

In the study, the determining of the reference velocity of the P-wave in coal seams used in seismic profiling to assess increases and decreases in relative stresses at large depths has been presented. The seismic profiling method proposed by Dubinski in 1989 covers a range of depth up to 970 m. At present, coal seams exploitation in Polish coal mines is conducted at greater depths, even exceeding 1200 m, which creates the necessity for a new reference velocity model. The study presents an empirical mathematical model of the change of the P-wave velocity in coal seams in the geological conditions of the Jastrzebie coal mine. A power model analogous to the Dubinski’s one was elaborated with new constants. The calculations included the results from 35 measurements of seismic profiling carried out in various coal seams of the Jastrzebie mine at depths from 640 to 1200 m. The results obtained cause changes in the result of calculations of seismic anomalies. Future validation of the proposed model with larger data set will be required.


2016 ◽  
Author(s):  
Beatriz Benjumea ◽  
Albert Macau ◽  
Anna Gabàs ◽  
Sara Figueras

Abstract. We combine geophysical well logging and passive seismic measurements to characterize the near surface geology of an area located in Hontomin, Burgos (Spain). This area has some near-surface challenges for a geophysical study. The irregular topography is characterized by limestone outcrops and unconsolidated sediments areas. Additionally, the near surface geology includes an upper layer of pure limestones overlying marly limestones and marls (Upper Cretaceous). These materials lie on top of Low Cretaceous siliciclastic sediments (sandstones, clays, gravels). In any case, decreasing seismic velocity with depth is expected. The geophysical datasets used in this study include sonic and gamma ray logs at two boreholes and passive seismic measurements: 224 H/V stations and 3 arrays. Well logging data defines two significant changes in the P-wave velocity log within the Upper Cretaceous layer and one more at the Upper to Lower Cretaceous contact. This technique has also used for refining the geological interpretation. The passive seismic measurements provide a map of sediment thickness with maximum of around 40 m and shear-wave velocity profiles from the array technique. A comparison between seismic velocity coming from well logging and array measurements defines the resolution limits of the passive seismic techniques and helps for its interpretation. This study shows how these low-cost techniques can provide useful information about near-surface complexity that could be used for designing a geophysical field survey or for seismic processing steps such as statics or imaging.


Sign in / Sign up

Export Citation Format

Share Document