Traveltime tomographic inversion with simultaneous static corrections — Well worth the effort

Geophysics ◽  
2009 ◽  
Vol 74 (6) ◽  
pp. WCB25-WCB33 ◽  
Author(s):  
Ari Tryggvason ◽  
Cedric Schmelzbach ◽  
Christopher Juhlin
Keyword(s):  
Seismic Reflection ◽  
Real Data ◽  
Seismic Velocities ◽  
Low Velocity ◽  
Data Set ◽  
Detailed Model ◽  
Seismic Reflection Data ◽  
Near Surface ◽  
First Arrival ◽  
Static Corrections

We have developed a first-arrival traveltime inversion scheme that jointly solves for seismic velocities and source and receiver static-time terms. The static-time terms are included to compensate for varying time delays introduced by the near-surface low-velocity layer that is too thin to be resolved by tomography. Results on a real data set consisting of picked first-arrival times from a seismic-reflection 2D/3D experiment in a crystalline environment show that the tomography static-time terms are very similar in values and distribution to refraction-static corrections computed using standard refraction-statics software. When applied to 3D seismic-reflection data, tomography static-time terms produce similar or more coherent seismic-reflection images compared to the images using corrections from standard refraction-static software. Furthermore, the method provides a much more detailed model of the near-surface bedrock velocity than standard software when the static-time terms are included in the inversion. Low-velocity zones in this model correlate with other geologic and geophysical data, suggesting that our method results in a reliable model. In addition to generally being required in seismic-reflection imaging, static corrections are also necessary in traveltime tomography to obtain high-fidelity velocity images of the subsurface.

Tomographic determination of velocity and depth in laterally varying media

Geophysics ◽  
1985 ◽  
Vol 50 (6) ◽  
pp. 903-923 ◽  
Author(s):  
T. N. Bishop ◽  
K. P. Bube ◽  
R. T. Cutler ◽  
R. T. Langan ◽  
P. L. Love ◽  
...  
Keyword(s):  
Seismic Reflection ◽  
Seismic Velocity ◽  
Real Data ◽  
Seismic Reflection Data ◽  
Near Surface ◽  
Reflection Data ◽  
Tomographic Method ◽  
Recorded Data ◽  
The Difference

Estimation of reflector depth and seismic velocity from seismic reflection data can be formulated as a general inverse problem. The method used to solve this problem is similar to tomographic techniques in medical diagnosis and we refer to it as seismic reflection tomography. Seismic tomography is formulated as an iterative Gauss‐Newton algorithm that produces a velocity‐depth model which minimizes the difference between traveltimes generated by tracing rays through the model and traveltimes measured from the data. The input to the process consists of traveltimes measured from selected events on unstacked seismic data and a first‐guess velocity‐depth model. Usually this first‐guess model has velocities which are laterally constant and is usually based on nearby well information and/or an analysis of the stacked section. The final model generated by the tomographic method yields traveltimes from ray tracing which differ from the measured values in recorded data by approximately 5 ms root‐mean‐square. The indeterminancy of the inversion and the associated nonuniqueness of the output model are both analyzed theoretically and tested numerically. It is found that certain aspects of the velocity field are poorly determined or undetermined. This technique is applied to an example using real data where the presence of permafrost causes a near‐surface lateral change in velocity. The permafrost is successfully imaged in the model output from tomography. In addition, depth estimates at the intersection of two lines differ by a significantly smaller amount than the corresponding estimates derived from conventional processing.

An example of extreme near-surface variability in shallow seismic reflection data

2016 ◽  
Vol 4 (3) ◽  
pp. SH1-SH9
Author(s):  
Steven D. Sloan ◽  
J. Tyler Schwenk ◽  
Robert H. Stevens
Keyword(s):  
Seismic Reflection ◽  
Field Data ◽  
Low Velocity ◽  
Seismic Reflection Data ◽  
Near Surface ◽  
Shallow Subsurface ◽  
Reflection Data ◽  
Shallow Seismic ◽  
Low Velocity Layer ◽  
Velocity Layer

Variability of material properties in the shallow subsurface presents challenges for near-surface geophysical methods and exploration-scale applications. As the depth of investigation decreases, denser sampling is required, especially of the near offsets, to accurately characterize the shallow subsurface. We have developed a field data example using high-resolution shallow seismic reflection data to demonstrate how quickly near-surface properties can change over short distances and the effects on field data and processed sections. The addition of a relatively thin, 20 cm thick, low-velocity layer can lead to masked reflections and an inability to map shallow reflectors. Short receiver intervals, on the order of 10 cm, were necessary to identify the cause of the diminished data quality and would have gone unknown using larger, more conventional station spacing. Combined analysis of first arrivals, surface waves, and reflections aided in determining the effects and extent of a low-velocity layer that inhibited the identification and constructive stacking of the reflection from a shallow water table using normal-moveout-based processing methods. Our results also highlight the benefits of using unprocessed gathers to pragmatically guide processing and interpretation of seismic data.

Linearized tomographic inversion of first‐arrival times

Geophysics ◽  
1992 ◽  
Vol 57 (11) ◽  
pp. 1482-1492 ◽  
Author(s):  
James L. Simmons ◽  
Milo M. Backus
Keyword(s):  
Reference Model ◽  
Depth Range ◽  
Inversion Algorithm ◽  
Data Set ◽  
Tomographic Inversion ◽  
Seismic Reflection Data ◽  
Near Surface ◽  
Geometric Patterns ◽  
First Arrival ◽  
Thick Layers

A linearized tomographic‐inversion algorithm estimates the near‐surface slowness anomalies present in a conventional, shallow‐marine seismic reflection data set. First‐arrival time residuals are the data to be inverted. The anomalies are treated as perturbations relative to a known, laterally‐invariant reference velocity model. Below the sea floor the reference model varies smoothly with depth; consequently the first arrivals are considered to be diving waves. In the offset‐midpoint domain the geometric patterns of traveltime perturbations produced by the anomalies resemble hyperbolas. Based on simple ray theory, these geometric patterns are predictable and can be used to relate the unknown model to the data. The assumption of a laterally‐invariant reference model permits an efficient solution in the offset‐wavenumber domain which is obtained in a single step using conventional least squares. The tomographic image shows the vertical‐traveltime perturbations associated with the anomalies as a function of midpoint at a number of depths. As implemented, the inverse problem is inherently stable. The first arrivals sample the subsurface to a maximum depth of roughly 500 m (≈ one‐fifth of the spread length). The model is parameterized to consist of fifteen 20-m thick layers spanning a depth range of 80–380 m. One‐way vertical‐traveltime delays as large as 10 ms are estimated. Assuming that these time delays are distributed over the entire 20-m thick layers, velocities much slower than water velocity are implied for the anomalies. Maps of the tomographic images show the spatial location and orientation of the anomalies throughout the prospect for the upper 400 m. Each line is processed independently, and the results are corroborated to a high degree at the line intersections.

scholarly journals Using Seismic Attributes in seismotectonic research: an application to the Norcia's Mw = 6.5 earthquake (30th October 2016) in Central Italy

2019 ◽  
Author(s):  
Maurizio Ercoli ◽  
Emanuele Forte ◽  
Massimiliano Porreca ◽  
Ramon Carbonell ◽  
Cristina Pauselli ◽  
...  
Keyword(s):  
Seismic Reflection ◽  
Central Italy ◽  
Seismic Attributes ◽  
Seismic Attribute ◽  
Data Set ◽  
Epicentral Zone ◽  
Seismic Reflection Data ◽  
Near Surface ◽  
Reflection Data ◽  
Attribute Analysis

Abstract. In seismotectonic studies, seismic reflection data are a powerful tool to unravel the complex deep architecture of active faults. Such tectonic structures are usually mapped at surface through traditional geological surveying whilst seismic reflection data may help to trace their continuation from the near-surface down to hypocentral depth. In this study, we propose the application of the seismic attributes technique, commonly used in seismic reflection exploration by oil industry, to seismotectonic research for the first time. The study area is a geologically complex region of Central Italy, recently struck by a long-lasting seismic sequence including a Mw 6.5 main-shock. A seismic reflection data-set consisting of three vintage seismic profiles, currently the only available across the epicentral zone, constitutes a singular opportunity to attempt a seismic attribute analysis. This analysis resulted in peculiar seismic signatures which generally correlate with the exposed surface geologic features, and also confirming the presence of other debated structures. These results are critical, because provide information also on the relatively deep structural setting, mapping a prominent, high amplitude regional reflector that marks the top basement, interpreted as important rheological boundary. Complex patterns of high-angle discontinuities crossing the reflectors have been also identified. These dipping fabrics are interpreted as the expression of fault zones, belonging to the active normal fault systems responsible for the seismicity of the region. This work demonstrates that seismic attribute analysis, even if used on low-quality vintage 2D data, may contribute to improve the subsurface geological interpretation of areas characterized by high seismic potential.

scholarly journals Time-lapse imaging using 3D ultra-high-frequency marine seismic reflection data

Geophysics ◽  
2020 ◽  
Vol 85 (2) ◽  
pp. P13-P25
Author(s):  
Michael J. Faggetter ◽  
Mark E. Vardy ◽  
Justin K. Dix ◽  
Jonathan M. Bull ◽  
Timothy J. Henstock
Keyword(s):  
High Frequency ◽  
Seismic Reflection ◽  
High Tide ◽  
Time Lapse ◽  
Data Set ◽  
Seismic Reflection Data ◽  
Near Surface ◽  
Ultra High Frequency ◽  
Reflection Data ◽  
Marine Seismic

Time-lapse (4D) seismic imaging is now widely used as a tool to map and interpret changes in deep reservoirs as well as investigate dynamic, shallow hydrological processes in the near surface. However, there are very few examples of time-lapse analysis using ultra-high-frequency (UHF; kHz range) marine seismic reflection data. Exacting requirements for navigation can be prohibitive for acquiring coherent, true-3D volumes. Variable environmental noise can also lead to poor amplitude repeatability and make it difficult to identify differences that are related to real physical changes. Overcoming these challenges opens up a range of potential applications for monitoring the subsurface at decimetric resolution, including geohazards, geologic structures, as well as the bed-level and subsurface response to anthropogenic activities. Navigation postprocessing was incorporated to improve the acquisition and processing workflow for the 3D Chirp subbottom profiler and provide stable, centimeter-level absolute positioning, resulting in well-matched 3D data and mitigating 4D noise for data stacked into [Formula: see text] common-midpoint bins. Within an example 4D data set acquired on the south coast of the UK, interpretable differences are recorded within a shallow gas blanket. Reflections from the top and bottom of a gas pocket are imaged at low tide, whereas at high tide only the upper reflection is imaged. This case study demonstrates the viability of time-lapse UHF 3D seismic reflection for quantitative mapping of decimeter-scale changes within the shallow marine subsurface.

Shallow 3-D seismic reflection surveying: Data acquisition and preliminary processing strategies

Geophysics ◽  
1998 ◽  
Vol 63 (4) ◽  
pp. 1434-1450 ◽  
Author(s):  
Frank Büker ◽  
Alan G. Green ◽  
Heinrich Horstmeyer
Keyword(s):  
Data Acquisition ◽  
Seismic Reflection ◽  
Full Range ◽  
Depth Range ◽  
Data Set ◽  
Processing Scheme ◽  
Seismic Reflection Data ◽  
Near Surface ◽  
Glacial Sediments ◽  
Reflection Data

A comprehensive strategy of 3-D seismic reflection data acquisition and processing has been used in a study of glacial sediments deposited within a Swiss mountain valley. Seismic data generated by a downhole shotgun source were recorded with single 30-Hz geophones distributed at 3 m × 3 m intervals across a 357 m × 432 m area. For most common‐midpoint (CMP) bins, traces covering a full range of azimuths and source‐receiver distances of ∼2 to ∼125 m were recorded. A common processing scheme was applied to the entire data set and to various subsets designed to simulate data volumes collected with lower density source and receiver patterns. Comparisons of seismic sections extracted from the processed 3-D subsets demonstrated that high‐fold (>40) and densely spaced (CMP bin sizes ⩽ 3 m × 3 m) data with relatively large numbers (>6) of traces recorded at short (<20 m) source‐receiver offsets were essential for obtaining clear images of the shallowest (<100 ms) reflecting horizons. Reflections rich in frequencies >100 Hz at traveltimes of ∼20 to ∼170 ms provided a vertical resolution of 3 to 6 m over a depth range of ∼15 to ∼150 m. The shallowest prominent reflection at 20 to 35 ms (∼15 to 27 m depth) originated from the boundary between a near‐surface sequence of clays/silts and an underlying unit of heterogeneous sands/gravels.

Aquifer heterogeneity from SH‐wave seismic impedance inversion

Geophysics ◽  
2002 ◽  
Vol 67 (5) ◽  
pp. 1548-1557 ◽  
Author(s):  
Kevin D. Jarvis ◽  
Rosemary J. Knight
Keyword(s):  
Shallow Aquifer ◽  
Velocity Variation ◽  
Bayesian Inversion ◽  
Cone Penetrometer ◽  
Seismic Velocities ◽  
Sh Wave ◽  
Data Set ◽  
Seismic Reflection Data ◽  
Near Surface ◽  
The Impact

We collected SH‐wave seismic reflection data over a shallow aquifer in southwestern British Columbia to investigate the use of such data in hydrogeologic applications. We used this data set in developing a methodology that uses cone penetrometer data as an integral part of the inversion and interpretation of the seismic data. A Bayesian inversion technique converts the seismic amplitude variations to velocity variations, honoring the probabilities of the priors and adhering to a geologically reasonable sparseness criterion. Velocity measurements acquired with the cone penetrometer provide velocity profiles and vertical seismic profiling (VSP) data, all of which are valuable in properly constraining the Bayesian inversion. The differentiation of lithologies (in this data set, sand and clay) is accomplished by first using a normalization procedure to remove the impact of effective stress, which dominates the velocity variation in the upper 10 to 20 m. The final transformation of seismic velocities to void ratio for the sand‐dominated regions is made using laboratory‐derived measurements; it provides an image of the heterogeneity of the near‐surface aquifer.

ESTIMATION AND CORRECTION OF NEAR‐SURFACE TIME ANOMALIES

Geophysics ◽  
1974 ◽  
Vol 39 (4) ◽  
pp. 441-463 ◽  
Author(s):  
M. Turhan Taner ◽  
F. Koehler ◽  
K. A. Alhilali
Keyword(s):  
Real Data ◽  
Least Square ◽  
Seismic Reflection Data ◽  
Near Surface ◽  
Reflection Data ◽  
Common Depth Point ◽  
Surface Time ◽  
Before And After ◽  
Nmo Correction ◽  
Static Corrections

The problem of computing static corrections for CDP seismic reflection data is discussed. A new approach is presented and it is related to various existing approaches. The approach consists of using crosscorrelation computations to find time shifts which appear to align the traces of each common‐depth‐point. These shifts are expressed in terms of surface corrections, one for each source and receiver position; a residual NMO correction for each common‐depth‐point; and a fixed correction for each common‐depth‐point. These simultaneous equations form an overdetermined set which can be solved for the unknown static and NMO corrections. The least‐square‐error solution to these equations has an important indeterminancy which is discussed. Methods for its resolution are proposed. Application of the technique to real data is illustrated by several examples. Validity of the corrections is demonstrated by velocity analyses before and after correction of the traces.

Prestack time migration of nonplanar data: Improving topography prestack time migration with dip-angle domain stationary-phase filtering and effective velocity inversion

Geophysics ◽  
2017 ◽  
Vol 82 (3) ◽  
pp. S235-S246 ◽  
Author(s):  
Jincheng Xu ◽  
Jianfeng Zhang
Keyword(s):  
Stationary Phase ◽  
Real Data ◽  
Data Set ◽  
Near Surface ◽  
Velocity Inversion ◽  
Effective Velocity ◽  
Time Migration ◽  
Static Corrections ◽  
Prestack Time Migration ◽  
Dip Angle

We have developed a modified prestack time migration (PSTM) approach that can directly image nonplanar data by using two effective velocity parameters above and below a datum. The proposed extension improves the so-called topography PSTM by introducing a dip-angle domain stationary-phase migration (or filtering) and combining effective velocity inversion with the residual static corrections. The stationary-phase migration to constrain the imaging aperture within Fresnel zones significantly improves the signal-to-noise ratio (S/N) of the image gathers, especially in the presence of steeply dipping structures. This helps to extract an accurate residual moveout from the common shot and receiver image gathers, and the surface-consistent residual statics hidden in these image gathers can be simultaneously obtained from an inversion process. As a result, the final migrated images show higher S/N and are better focused than the conventional topography PSTM. The proposed technique can handle rugged topography, especially in the presence of high near-surface velocities, without the need for prior elevation static corrections. The SEG foothills overthrust model and a real data set acquired on a piedmont zone are used to validate the modified topography PSTM. Synthetic and field data examples are obtained with good results.

Sign in / Sign up

Export Citation Format

Share Document