Experimental and texture-derived P-wave anisotropy of principal rocks from the TRANSALP traverse: An aid for the interpretation of seismic field data

2006 ◽  
Vol 414 (1-4) ◽  
pp. 97-116 ◽  
Author(s):  
Klaus Ullemeyer ◽  
Siegfried Siegesmund ◽  
Patrick N.J. Rasolofosaon ◽  
Jan H. Behrmann
Keyword(s):  
Field Data ◽  
P Wave

scholarly journals Finite-difference modelling to evaluate seismic P-wave and shear-wave field data

Solid Earth ◽  
2015 ◽  
Vol 6 (1) ◽  
pp. 33-47 ◽  
Author(s):  
T. Burschil ◽  
T. Beilecke ◽  
C. M. Krawczyk
Keyword(s):  
Wave Field ◽  
Shear Wave ◽  
Field Data ◽  
P Wave ◽  
Sh Wave ◽  
Surface Shear ◽  
Near Surface ◽  
Reflection Seismic ◽  
Surface Shear Wave ◽  
Seismic Measurements

Abstract. High-resolution reflection seismic methods are an established non-destructive tool for engineering tasks. In the near surface, shear-wave reflection seismic measurements usually offer a higher spatial resolution in the same effective signal frequency spectrum than P-wave data, but data quality varies more strongly. To discuss the causes of these differences, we investigated a P-wave and a SH-wave seismic reflection profile measured at the same location on the island of Föhr, Germany and applied seismic reflection processing to the field data as well as finite-difference modelling of the seismic wave field. The simulations calculated were adapted to the acquisition field geometry, comprising 2 m receiver distance (1 m for SH wave) and 4 m shot distance along the 1.5 km long P-wave and 800 m long SH-wave profiles. A Ricker wavelet and the use of absorbing frames were first-order model parameters. The petrophysical parameters to populate the structural models down to 400 m depth were taken from borehole data, VSP (vertical seismic profile) measurements and cross-plot relations. The simulation of the P-wave wave-field was based on interpretation of the P-wave depth section that included a priori information from boreholes and airborne electromagnetics. Velocities for 14 layers in the model were derived from the analysis of five nearby VSPs (vP =1600–2300 m s-1). Synthetic shot data were compared with the field data and seismic sections were created. Major features like direct wave and reflections are imaged. We reproduce the mayor reflectors in the depth section of the field data, e.g. a prominent till layer and several deep reflectors. The SH-wave model was adapted accordingly but only led to minor correlation with the field data and produced a higher signal-to-noise ratio. Therefore, we suggest to consider for future simulations additional features like intrinsic damping, thin layering, or a near-surface weathering layer. These may lead to a better understanding of key parameters determining the data quality of near-surface shear-wave seismic measurements.

scholarly journals Investigation of ULF magnetic pulsations, air conductivity changes, and infra red signatures associated with the 30 October Alum Rock M5.4 earthquake

2009 ◽  
Vol 9 (2) ◽  
pp. 585-603 ◽  
Author(s):  
T. Bleier ◽  
C. Dunson ◽  
M. Maniscalco ◽  
N. Bryant ◽  
R. Bambery ◽  
...  
Keyword(s):  
Field Data ◽  
Low Frequency ◽  
Field Tests ◽  
P Wave ◽  
Electromagnetic Signal ◽  
San Jose ◽  
Data Set ◽  
Bay Area ◽  
Conductivity Sensor ◽  
Magnetic Pulsations

Abstract. Several electromagnetic signal types were observed prior to and immediately after 30 October 2007 (Local Time) M5.4 earthquake at Alum Rock, Ca with an epicenter ~15 km NE of San Jose Ca. The area where this event occurred had been monitored since November 2005 by a QuakeFinder magnetometer site, unit 609, 2 km from the epicenter. This instrument is one of 53 stations of the QuakeFinder (QF) California Magnetometer Network-CalMagNet. This station included an ultra low frequency (ULF) 3-axis induction magnetometer, a simple air conductivity sensor to measure relative airborne ion concentrations, and a geophone to identify the arrival of the P-wave from an earthquake. Similar in frequency content to the increased ULF activity reported two weeks prior to the Loma Prieta M7.0 quake in 1989 (Fraser-Smith, 1990, 1991), the QF station detected activity in the 0.01–12 Hz bands, but it consisted of an increasing number of short duration (1 to 30 s duration) pulsations. The pulsations peaked around 13 days prior to the event. The amplitudes of the pulses were strong, (3–20 nT), compared to the average ambient noise at the site, (10–250 pT), which included a component arising from the Bay Area Rapid Transit (BART) operations. The QF station also detected different pulse shapes, e.g. negative or positive only polarity, with some pulses including a combination of positive and negative. Typical pulse counts over the previous year ranged from 0–15 per day, while the count rose to 176 (east-west channel) on 17 October, 13 days prior to the quake. The air conductivity sensor saturated for over 14 h during the night and morning prior to the quake, which occurred at 20:29 LT. Anomalous IR signatures were also observed in the general area, within 50 km of the epicenter, during the 2 weeks prior to the quake. These three simultaneous EM phenomena were compared with data collected over a 1–2-year period at the site. The data was also compared against accounts of air ionization reported to be associated with radon emission from the ground (Ouzounov, 2007), and a series of laboratory rock stressing experiments (Freund, 2006, 2007a, b, c) to determine if field data was consistent either of these accounts. We could not find a data set with pre-earthquake radon measurements taken near the Alum Rock epicenter to compare against our field data. However, based on the Alum Rock data set example and another data set at Parkfield, the field tests are at least consistent with the lab experiments showing currents, magnetic field disturbances, air conductivity changes, and IR signatures. This is encouraging, but more instrumented earthquake examples are needed to prove a repeating pattern for these types of pre-earthquake EM signatures. For more information on QuakeFinder please view http://www.quakefinder.com.

Integration of surface-based tomographic models for zonation and multimodel guided extrapolation of sparsely known petrophysical parameters

Geophysics ◽  
2013 ◽  
Vol 78 (4) ◽  
pp. EN43-EN53 ◽  
Author(s):  
Barbara Hachmöller ◽  
Hendrik Paasche
Keyword(s):  
Cluster Analysis ◽  
Field Data ◽  
Spectral Cluster ◽  
P Wave ◽  
Fuzzy Membership ◽  
Geophysical Model ◽  
Natural Gamma ◽  
Refraction Seismic ◽  
Spatially Varying ◽  
Membership Information

We integrate the information of multiple tomographic models acquired from the earth’s surface by modifying a statistical approach recently developed for the integration of cross-borehole tomographic models. In doing so, we introduce spectral cluster analysis as the new core of the model integration procedure to capture the spatial heterogeneity present in all considered tomographic models and describe this heterogeneity in a fuzzy sense. Because spectral cluster algorithms analyze model structure locally, they are considered relatively robust with regard to systematically and spatially varying imaging capabilities typical for geophysical tomographic surveys conducted on the earth’s surface. Using a synthetic aquifer example, a fuzzy spectral cluster algorithm can be used to integrate the information provided by 2D tomographic refraction seismic and DC resistivity surveys. The integrated information in the fuzzy membership domain is then used to derive an integrated zonal geophysical model outlining the major structural units present in both input models. We also explain how the fuzzy membership information can be used to identify optimal locations for sparse logging of additional target parameters, i.e., porosity information in our synthetic example. We demonstrate how this sparse porosity information can be extrapolated based on all tomographic input models. The resultant 2D porosity model matches the original porosity distribution reasonably well within the spatial resolution limits of the underlying tomographic models. Consecutively, we apply this approach to a field data base acquired over a former river channel. Sparse information about natural gamma radiation is available and extrapolated on the basis of the fuzzy membership information obtained by spectral cluster analysis of 2D P-wave velocity and electrical resistivity models. This field data shows that the presented parameter extrapolation procedure is robust, even if the locations of target parameter acquisition have not been optimized with regard to the fuzzy membership information.

Feasibility of estimating vertical transverse isotropy from microseismic data recorded by surface monitoring arrays

Geophysics ◽  
2011 ◽  
Vol 76 (6) ◽  
pp. WC117-WC126 ◽  
Author(s):  
Davide Gei ◽  
Leo Eisner ◽  
Peter Suhadolc
Keyword(s):  
Vertical Velocity ◽  
Field Data ◽  
P Wave ◽  
Time Inversion ◽  
Source Depth ◽  
Effective Anisotropy ◽  
Arrival Times ◽  
Microseismic Data ◽  
Surface Monitoring ◽  
Wave Arrival

Microseismic data recorded by surface monitoring arrays can be used to estimate the effective anisotropy of the overburden and reservoir. In this study we used the inversion of picked P-wave arrival times to estimate the Thomsen parameter [Formula: see text] and the anellipticity coefficient [Formula: see text]. This inversion employs an analytic equation of P-wave traveltimes as a function of offset in homogeneous, transversely isotropic media with a vertical axis of symmetry. We considered a star-like distribution of receivers and, for this geometry, we analyzed the sensitivity of the inversion method to picking noise and to uncertainties in the P-wave vertical velocity and source depth. Long offsets, as well as a high number of receivers per line, improve the estimation of [Formula: see text] and [Formula: see text] from noisy arrival times. However, if we do not use the correct value of the P-wave vertical velocity or source depth, the long-offset may increase the inaccuracy in the estimation of the anisotropic parameters. Such inaccuracy cannot be detected from time residuals. We also applied this inversion to field data acquired during the hydraulic fracturing of a gas shale reservoir and compared the results with the anisotropic parameters estimated from synthetic arrival times computed for an isotropic layered medium. The effective anisotropy from the inversion of the field data cannot be explained by layering only and is partially due to the intrinsic anisotropy of the reservoir and/or overburden. This study emphasizes the importance of using accurate values of the vertical velocity and source depth in the P-wave arrival time inversion for estimating anisotropic parameters from passive seismic data.

scholarly journals Seismic interferometry using multidimensional deconvolution and crosscorrelation for crosswell seismic reflection data without borehole sources

Geophysics ◽  
2011 ◽  
Vol 76 (1) ◽  
pp. SA19-SA34 ◽  
Author(s):  
Shohei Minato ◽  
Toshifumi Matsuoka ◽  
Takeshi Tsuji ◽  
Deyan Draganov ◽  
Jürg Hunziker ◽  
...  
Keyword(s):  
Seismic Reflection ◽  
Field Data ◽  
Time Lapse ◽  
P Wave ◽  
Source Distribution ◽  
Seismic Interferometry ◽  
Imaging Method ◽  
Seismic Reflection Data ◽  
Reflection Data ◽  
Receiver Arrays

Crosswell reflection method is a high-resolution seismic imaging method that uses recordings between boreholes. The need for downhole sources is a restrictive factor in its application, for example, to time-lapse surveys. An alternative is to use surface sources in combination with seismic interferometry. Seismic interferometry (SI) could retrieve the reflection response at one of the boreholes as if from a source inside the other borehole. We investigate the applicability of SI for the retrieval of the reflection response between two boreholes using numerically modeled field data. We compare two SI approaches — crosscorrelation (CC) and multidimensional deconvolution (MDD). SI by MDD is less sensitive to underillumination from the source distribution, but requires inversion of the recordings at one of the receiver arrays from all the available sources. We find that the inversion problem is ill-posed, and propose to stabilize it using singular-value decomposition. The results show that the reflections from deep boundaries are retrieved very well using both the CC and MDD methods. Furthermore, the MDD results exhibit more realistic amplitudes than those from the CC method for downgoing reflections from shallow boundaries. We find that the results retrieved from the application of both methods to field data agree well with crosswell seismic-reflection data using borehole sources and with the logged P-wave velocity.

scholarly journals Exploitation of Field Drilling Data with an Innovative Drilling Simulator: Highly Effective Simulation of Rotating and Sliding Mode

2021 ◽  
Author(s):  
Alexis Koulidis ◽  
Vassilios Kelessidis ◽  
Shehab Ahmed
Keyword(s):  
Field Data ◽  
Sliding Mode ◽  
P Wave ◽  
Drilling Process ◽  
Continuous Change ◽  
Rate Of Penetration ◽  
Drilling Performance ◽  
Drilling Parameters ◽  
Rock Drillability ◽  
Drilling Data

Abstract Drilling challenging wells requires a combination of drilling analytics and comprehensive simulation to prevent poor drilling performance and avoid drilling issues for the upcoming drilling campaign. This work focuses on the capabilities of a drilling simulator that can simulate the directional drilling process with the use of actual field data for the training of students and professionals. This paper presents the results of simulating both rotating and sliding modes and successfully matching the rate of penetration and the trajectory of an S-type well. Monitored drilling data from the well were used to simulate the drilling process. These included weight on bit, revolutions per minute, flow rate, bit type, inclination and drilling fluid properties. The well was an S-type well with maximum inclination of 16 degrees. There were continuous variations from rotating to sliding mode, and the challenge was approached by classifying drilling data into intervals of 20 feet to obtain an appropriate resolution and efficient simulation. The simulator requires formation strength, pore and fracture pressures, and details of well lithology, thus simulating the actual drilling environment. The uniaxial compressive strength of the rock layer is calculated from p–wave velocity data from an offset field. Rock drillability is finally estimated as a function of the rock properties of the drilled layer, bit type and the values of the drilling parameters. It is then converted to rate of penetration and matched to actual data. Changes in the drilling parameters were followed as per the field data. The simulator reproduces the drilling process in real-time and allows the driller to make instantaneous changes to all drilling parameters. The simulator provides the rate of penetration, torque, standpipe pressure, and trajectory as output. This enables the user to have on-the-fly interference with the drilling process and allows him/her to modify any of the important drilling parameters. Thus, the user can determine the effect of such changes on the effectiveness of drilling, which can lead to effective drilling optimization. Certain intervals were investigated independently to give a more detailed analysis of the simulation outcome. Additional drilling data such as hook load and standpipe pressure were analyzed to determine and evaluate the drilling performance of a particular interval and to consider them in the optimization process. The resulting rate of penetration and well trajectory simulation results show an excellent match with field data. The simulation illustrates the continuous change between rotating and sliding mode as well as the accurate synchronous matching of the rate of penetration and trajectory. The results prove that the simulator is an excellent tool for students and professionals to simulate the drilling process prior to actual drilling of the next inclined well.

Seismic characterization of submarine gas-hydrate deposits in the Western Black Sea by acoustic full-waveform inversion of ocean-bottom seismic data

Geophysics ◽  
2019 ◽  
Vol 84 (5) ◽  
pp. B311-B324 ◽  
Author(s):  
Laura Gassner ◽  
Tobias Gerach ◽  
Thomas Hertweck ◽  
Thomas Bohlen
Keyword(s):  
Wave Velocity ◽  
Black Sea ◽  
Gas Hydrate ◽  
Seismic Data ◽  
Field Data ◽  
Waveform Inversion ◽  
P Wave ◽  
Full Waveform Inversion ◽  
Ocean Bottom ◽  
Full Waveform

Evidence for gas-hydrate occurrence in the Western Black Sea is found from seismic measurements revealing bottom-simulating reflectors (BSRs) of varying distinctness. From an ocean-bottom seismic data set, low-resolution traveltime-tomography models of P-wave velocity [Formula: see text] are constructed. They serve as input for acoustic full-waveform inversion (FWI), which we apply to derive high-resolution parameter models aiding the interpretation of the seismic data for potential hydrate and gas deposits. Synthetic tests indicate the applicability of the FWI approach to robustly reconstruct [Formula: see text] models with a typical hydrate and gas signature. Models of S-wave velocity [Formula: see text] containing a hydrate signature can only be reconstructed when the parameter distribution of [Formula: see text] is already well-known. When we add noise to the modeled data to simulate field-data conditions, it prevents the reconstruction of [Formula: see text] completely, justifying the application of an acoustic approach. We invert for [Formula: see text] models from field data of two parallel profiles of 14 km length with a distance of 1 km. Results indicate a characteristic velocity trend for hydrate and gas occurrence at BSR depth in the first of the analyzed profiles. We find no indications for gas accumulations below the BSR on the second profile and only weak indications for hydrate. These differences in the [Formula: see text] signature are consistent with the reflectivity behavior of the migrated seismic streamer data of both profiles in which a zone of high-reflectivity amplitudes is coincident with the potential gas zone derived from the FWI result. Calculating saturation estimates for the potential hydrate and gas zones yields values of up to 30% and 1.2%, respectively.

Application of waveform tomography to marine seismic reflection data from the Queen Charlotte Basin of western Canada

Geophysics ◽  
2011 ◽  
Vol 76 (2) ◽  
pp. B55-B70 ◽  
Author(s):  
E. M. Takam Takougang ◽  
A. J. Calvert
Keyword(s):  
Seismic Reflection ◽  
Field Data ◽  
Velocity Model ◽  
P Wave ◽  
Low Velocity ◽  
Seismic Line ◽  
Seismic Reflection Data ◽  
Reflection Data ◽  
Sonic Log ◽  
Waveform Tomography

To obtain a higher resolution quantitative P-wave velocity model, 2D waveform tomography was applied to seismic reflection data from the Queen Charlotte sedimentary basin off the west coast of Canada. The forward modeling and inversion were implemented in the frequency domain using the visco-acoustic wave equation. Field data preconditioning consisted of f-k filtering, 2D amplitude scaling, shot-to-shot amplitude balancing, and time windowing. The field data were inverted between 7 and 13.66 Hz, with attenuation introduced for frequencies ≥ 10.5 Hz to improve the final velocity model; two different approaches to sampling the frequencies were evaluated. The limited maximum offset of the marine data (3770 m) and the relatively high starting frequency (7 Hz) were the main challenges encountered during the inversion. An inversion strategy that successively recovered shallow-to-deep structures was designed to mitigate these issues. The inclusion of later arrivals in the waveform tomography resulted in a velocity model that extends to a depth of approximately 1200 m, twice the maximum depth of ray coverage in the ray-based tomography. Overall, there is a good agreement between the velocity model and a sonic log from a well on the seismic line, as well as between modeled shot gathers and field data. Anomalous zones of low velocity in the model correspond to previously identified faults or their upward continuation into the shallow Pliocene section where they are not readily identifiable in the conventional migration.

Tunnel signature prediction for a cross‐borehole seismic survey

Geophysics ◽  
1995 ◽  
Vol 60 (1) ◽  
pp. 76-86 ◽  
Author(s):  
Richard D. Rechtien ◽  
Roy J. Greenfield ◽  
Robert F. Ballard
Keyword(s):  
Cross Section ◽  
Field Data ◽  
Rectangular Cross Section ◽  
Circular Cross Section ◽  
Data Interpretation ◽  
P Wave ◽  
Seismic Survey ◽  
S Waves ◽  
Borehole Seismic ◽  
Tunnel Size

Seismic location of tunnels or voids with a cross‐borehole survey is examined with field data and theory. The field data were taken at a site with a 2.2-m high by 2.7-m wide, roughly rectangular cross‐section tunnel, using a newly developed 1 to 5 kHz system employing a P‐wave sparker source. The synthetic records were obtained using a 2.5-D boundary‐valued solution for an explosive point source near a cylindrical void, and the solution was evaluated with the method of steepest descent. The synthetic waveforms compared well to the field data; both showed a maximum reduction of amplitude in the tunnel shadow of 8 dB and a maximum first arrival delay of 0.1 ms. Additional theoretical modeling was used to examine the variations of the received signals with tunnel size and frequency and showed amplitude reduction increased with frequency and tunnel size. Calculations showed that S‐waves scattered from the tunnel are more than 20 dB smaller than the primary P‐wave on hydrophones and more than 12 dB smaller on particle velocity sensors and so could be difficult to see in field data. The close comparison of synthetic waveforms to the field data indicate that the cylindrical model can be used to model data for roughly square cross‐section tunnels or voids, as well as for circular cross‐section tunnels, and thus is useful for data interpretation and survey planning.

Sign in / Sign up

Export Citation Format

Share Document