A study of methods for tomographic velocity estimation in the presence of low‐velocity zones

Geophysics ◽  
1985 ◽  
Vol 50 (6) ◽  
pp. 969-988 ◽  
Author(s):  
Sven Ivansson

This paper deals with the problem of seismic velocity estimation from first‐arrival traveltimes in a two‐dimensional (2-D) cross‐hole geometry where explosions are detonated in one borehole while recordings are made in another borehole and on the surface. Standard tomographic procedures are based on decomposition of the cross‐hole area into a number of cells and a simplifying assumption of straight raypaths. In the presence of significant low‐velocity zones, the resulting images may be contaminated. Different ways of performing tomographic inversion are tested on a number of synthetic examples. Images obtained by direct, unrestricted least‐squares inversion are often seriously distorted. However, methods using more cells and some kind of damping often give more satisfactory results. Because the risk of distorted images is always present in inversion procedures, comparison with synthetic data (forward modeling) is a valuable tool in the interpretation process. With a reasonably good initial solution, improvements can often be achieved by using iterative procedures to take account of ray‐bending affects as proposed in Bois et al.(1971). An alternative way of performing these calculations is described.


Geophysics ◽  
1992 ◽  
Vol 57 (2) ◽  
pp. 334-342 ◽  
Author(s):  
Larry R. Lines ◽  
Kenneth R. Kelly ◽  
John Queen

Layered geological formations with large seismic velocity contrasts can effectively create channel waves in cross‐borehole seismic data. The existence of channel waves for such waveguides can be confirmed by ray tracing, wave equation modeling, and modal analysis. Channel wave arrivals are identified in cross‐borehole data recorded at Conoco’s Newkirk test facility. For these data, where velocity contrasts are about 2 to 1, tomography based on first arrival traveltimes, is limited due to problems with extreme ray bending and seismic shadow zones. However, it may be possible to extract geological information using channel wave information. The seismometer differencing method appears to be a promising approach for detecting waveguide boundaries by use of cross‐borehole data.



2021 ◽  
Author(s):  
Kota Mukumoto ◽  
Takeshi Tsuji

<p>In Japan, seismic velocity structures have been estimated by using first arrival tomography method. Many significant crustal structures such as the coordinate of the subducted Philippine Sea plate has been revealed by seismic tomographic images. In this study, we applied the adjoint tomography including full numerical simulation and finite frequency sensitivity kernels for the area of central Japan. The study area is characterized by the very heterogeneous geologic structures. We used 72 natural earthquakes in this study. Because the dominant phase used in our analysis is the surface wave, only S-wave velocity was inverted. We tried to minimize the time-frequency phase misfit between observed and calculated waveforms with the frequency of 0.033~0.1Hz. Based on the checker bord test, our inversion scenario resolved the upper and lower crust. From the results, we identified more heterogeneous structures compared to those from the first arrival tomography. The estimated S-wave velocity model clearly resolved the low velocity anomalies around the active volcanoes. Furthermore, the velocity boundaries agree with the main tectonic lines in the central Japan.</p>



Geophysics ◽  
2001 ◽  
Vol 66 (2) ◽  
pp. 627-636 ◽  
Author(s):  
Pantelis M. Soupios ◽  
Constantinos B. Papazachos ◽  
Christopher Juhlin ◽  
Gregory N. Tsokas

This paper deals with the problem of nonlinear seismic velocity estimation from first‐arrival traveltimes obtained from crosshole and downhole experiments in three dimensions. A standard tomographic procedure is applied, based on the representation of the crosshole area into a number of cells which have an initial slowness assigned. For the forward modeling, the raypath matrix is computed using the revisited ray bending method, supplemented by an approximate computation of the first Fresnel zone at each point of the ray, hence using physical and not only mathematical rays. Since 3-D ray tracing is incorporated, the inversion technique is nonlinear. Velocity images are obtained by a constrained least‐squares inversion scheme using both “damping” and “smoothing” factors. The appropriate choice of these factors is defined by the use of appropriate criteria such as the L-curve. The tomographic approach is improved by incorporating a priori information about the media to be imaged into our inversion scheme. This improvement in imaging is achieved by projecting a desirable solution onto the null space of the inversion, and including this null‐space contribution with the standard non‐null‐space inversion solution. The efficiency of the inversion scheme is tested through a series of tests with synthetic data. Moreover, application in the area of the Ural Mountains using real data demonstrates that the proposed technique produces more realistic velocity models than those obtained by other standard approaches.



2016 ◽  
Vol 3 (3) ◽  
pp. 99-105
Author(s):  
Erick López Reyes ◽  
Paola Peralta Mendoza

 Se presenta una breve síntesis de la historia prehispánica y colonial de la región santaelenense, desde la primera llegada de grupos humanos a la zona hasta los primeros días de la República del Ecuador. Dando cuerpo narrativo a una serie de datos dispersos, producto de las investigaciones arqueológicas e histórico documentales, efectuadas en la región desde inicios del siglo pasado hasta la actualidad. En base a ellos se plantea entonces una necesaria y pertinente reflexión respecto a temas como la continuidad étnica, etnicidades, cultura e identidad de los actuales grupos humanos, descendientes de los anteriores, que hoy la habitan, planteando además la importancia que esto posee dentro de los procesos en curso, y aquellos potenciales de acometer, en pro del adelanto y desarrollo socioeconómico local a través del nuevo turismo o turismo de innovación.  Abstract This paper presents a brief summary of the pre-Hispanic and colonial history of the Santa Elena region, from the first arrival of human groups in the area until the early days of the new Republic of Ecuador. By giving narrative body to a series of data dispersed product of archaeological research, and historical documentaries, made in the region since early last century actually. Based on them then arises a necessary and relevant reflection on issues such as ethnic continuity, ethnicities, culture and identity of the current human groups, descendants of the previous, now inhabit, also raising the importance this has in ongoing processes, and those potential to undertake, for the advancement and local socio-economic development through New tourism or innovation tourism.



Geophysics ◽  
2021 ◽  
pp. 1-35
Author(s):  
M. Javad Khoshnavaz

Building an accurate velocity model plays a vital role in routine seismic imaging workflows. Normal-moveout-based seismic velocity analysis is a popular method to make the velocity models. However, traditional velocity analysis methodologies are not generally capable of handling amplitude variations across moveout curves, specifically polarity reversals caused by amplitude-versus-offset anomalies. I present a normal-moveout-based velocity analysis approach that circumvents this shortcoming by modifying the conventional semblance function to include polarity and amplitude correction terms computed using correlation coefficients of seismic traces in the velocity analysis scanning window with a reference trace. Thus, the proposed workflow is suitable for any class of amplitude-versus-offset effects. The approach is demonstrated to four synthetic data examples of different conditions and a field data consisting a common-midpoint gather. Lateral resolution enhancement using the proposed workflow is evaluated by comparison between the results from the workflow and the results obtained by the application of conventional semblance and three semblance-based velocity analysis algorithms developed to circumvent the challenges associated with amplitude variations across moveout curves, caused by seismic attenuation and class II amplitude-versus-offset anomalies. According to the obtained results, the proposed workflow is superior to all the presented workflows in handling such anomalies.



Geophysics ◽  
2019 ◽  
Vol 85 (1) ◽  
pp. U21-U29
Author(s):  
Gabriel Fabien-Ouellet ◽  
Rahul Sarkar

Applying deep learning to 3D velocity model building remains a challenge due to the sheer volume of data required to train large-scale artificial neural networks. Moreover, little is known about what types of network architectures are appropriate for such a complex task. To ease the development of a deep-learning approach for seismic velocity estimation, we have evaluated a simplified surrogate problem — the estimation of the root-mean-square (rms) and interval velocity in time from common-midpoint gathers — for 1D layered velocity models. We have developed a deep neural network, whose design was inspired by the information flow found in semblance analysis. The network replaces semblance estimation by a representation built with a deep convolutional neural network, and then it performs velocity estimation automatically with recurrent neural networks. The network is trained with synthetic data to identify primary reflection events, rms velocity, and interval velocity. For a synthetic test set containing 1D layered models, we find that rms and interval velocity are accurately estimated, with an error of less than [Formula: see text] for the rms velocity. We apply the neural network to a real 2D marine survey and obtain accurate rms velocity predictions leading to a coherent stacked section, in addition to an estimation of the interval velocity that reproduces the main structures in the stacked section. Our results provide strong evidence that neural networks can estimate velocity from seismic data and that good performance can be achieved on real data even if the training is based on synthetics. The findings for the 1D problem suggest that deep convolutional encoders and recurrent neural networks are promising components of more complex networks that can perform 2D and 3D velocity model building.



Geology ◽  
2019 ◽  
Vol 47 (8) ◽  
pp. 695-699 ◽  
Author(s):  
Victor E. Camp

AbstractBimodal volcanism and rhyolite migration along the High Lava Plains in central Oregon (United States) lie above a broader feature defined by low seismic velocity in the upper mantle that emanates from the Yellowstone hotspot (northwest United States) and extends westward across the northern Basin and Range. It was emplaced by a westward current, driven in part by rapid buoyancy-driven flow across the east-west cratonic boundary of North America. Geothermometry studies and geochemical considerations suggest that the low-velocity feature may be composed of moderately hot, low-density mantle derived from the Yellowstone plume but diluted by thermomechanical erosion and entrainment of colder mantle lithosphere. Finger-like conduits of plume-modified mantle beneath Quaternary eruption sites delineate flow-line channels that have developed across the broader mantle structure since 2 Ma. These channels have allowed low-density mantle to accumulate against the Cascades arc, thus providing a heated mantle source for mafic magmatism in the Newberry (Oregon) and Medicine Lake (California) volcanic fields.



2019 ◽  
Vol 219 (2) ◽  
pp. 1447-1462 ◽  
Author(s):  
Alexandre P Plourde ◽  
Michael G Bostock

SUMMARY We introduce a new relative moment tensor (MT) inversion method for clusters of nearby earthquakes. The method extends previous work by introducing constraints from S-waves that do not require modal decomposition and by employing principal component analysis to produce robust estimates of excitation. At each receiver, P and S waves from each event are independently aligned and decomposed into principal components. P-wave constraints on MTs are obtained from a ratio of coefficients corresponding to the first principal component, equivalent to a relative amplitude. For S waves we produce constraints on MTs involving three events, where one event is described as a linear combination of the other two, and coefficients are derived from the first two principal components. Nonlinear optimization is applied to efficiently find best-fitting tensile-earthquake and double-couple solutions for relative MT systems. Using synthetic data, we demonstrate the effectiveness of the P and S constraints both individually and in combination. We then apply the relative MT inversion to a set of 16 earthquakes from southern Alaska, at ∼125 km depth within the subducted Yakutat terrane. Most events are compatible with a stress tensor dominated by downdip tension, however, we observe several pairs of earthquakes with nearly antiparallel slip implying that the stress regime is heterogeneous and/or faults are extremely weak. The location of these events near the abrupt downdip termination of seismicity and the low-velocity zone suggest that they are caused by weakening via grain-size and volume reduction associated with eclogitization of the lower crustal gabbro layer.



Geophysics ◽  
1981 ◽  
Vol 46 (8) ◽  
pp. 1088-1099 ◽  
Author(s):  
Robert B. Rice ◽  
Samuel J. Allen ◽  
O. James Gant ◽  
Robert N. Hodgson ◽  
Don E. Larson ◽  
...  

Advances in exploration geophysics have continued apace during the last six years. We have entered a new era of exploration maturity which will be characterized by the extension of our technologies to their ultimate limits of precision. In gravity and magnetics, new inertial navigation systems permit the very rapid helicopter‐supported land acquisition of precise surface gravity data which is cost‐effective in regions of severe topography. Considerable effort is being expended to obtain airborne gravity data via helicopter which is of exploration quality. Significant progress has also been made in processing and interpreting potential field data. The goal of deriving the maximum amount of accurate subsurface information from seismic data has led to much more densely sampled and precise 2- and 3-D land data acquisition techniques. Land surveying accuracy has been greatly improved. The number of individually recorded detector channels has been increased dramatically (up to 1024) in order to approximate much more accurately a point‐source, point‐detector system. Much more powerful compressional‐wave vibrators can now maintain full force while sweeping up or down from 5 Hz to over 200 Hz. In marine surveying, new streamer cables and shipboard instrumentation permit the recording and limited processing of 96 to 480 channels. Improvements have also been made in marine sources and arrays. The most important developments in seismic data processing—wave‐equation based imaging and inversion methods—may be the forerunners of a totally new processing methodology. Wave‐equation methods have been formulated for migration before and after stack, multiples suppression, datum and replacement statics, velocity estimation, and seismic inversion. Inversion techniques which provide detailed acoustic‐impedance or velocity estimates have found widespread commercial application. Wavelet processing has greatly expanded our stratigraphic analysis capabilities. Much more sophisticated 1-, 2-, and 3-D modeling techniques are being used effectively to guide data acquisition and processing, as direct interpretation aids, and to teach basic interpretation concepts. Some systems can now handle vertical and lateral velocity changes, inelastic attenuation, curved reflection horizons, transitional boundaries, time‐variant waveforms, ghosting, multiples, and array‐response effects. Improved seismic display formats and the extensive use of color have been valuable in data processing, modeling, and interpretation. Stratigraphic interpretation has evolved into three major categories: (1) macrostratigraphy, where regional and basinal depositional patterns are analyzed to describe the broad geologic depositional environment; (2) qualitative stratigraphy, where specific rock units and their properties are analyzed qualitatively to delineate lithology, porosity, structural setting, and areal extent and shape; and (3) quantitative stratigraphy, where anomalies are mapped at a specific facies level to define net porosity‐feet distribution, gas‐fluid contacts, and probable pore fill. In essence, what began as direct hydrocarbon‐indicator technology applicable primarily to Upper Tertiary clastics has now matured to utility in virtually every geologic province. Considerable effort has been expended on the direct generation and recording of shear waves in an attempt to obtain more information about stratigraphy, porosity, and oil and gas saturation. Seismic service companies now offer shear‐wave prospecting using vibrator, horizontal‐impact, or explosive sources. Well logging has seen the acceleration of computerization. Wellsite tape recorders and minicomputers with relatively simple interpretation algorithms are routinely available. More sophisticated computerized interpretation methods are offered as a service at data processing centers.



Geophysics ◽  
1992 ◽  
Vol 57 (8) ◽  
pp. 1034-1047 ◽  
Author(s):  
Biondo Biondi

Imaging seismic data requires detailed knowledge of the propagation velocity of compressional waves in the subsurface. In conventional seismic processing, the interval velocity model is usually derived from stacking velocities. Stacking velocities are determined by measuring the coherency of the reflections along hyperbolic moveout trajectories in offset. This conventional method becomes inaccurate in geologically complex areas because the conversion of stacking velocities to interval velocities assumes a horizontally stratified medium and mild lateral variations in velocity. The tomographic velocity estimation proposed in this paper can be applied when there are dipping reflectors and strong lateral variations. The method is based on the measurements of moveouts by beam stacks. A beam stack measures local coherency of reflections along hyperbolic trajectories. Because it is a local operator, the beam stack can provide information on nonhyperbolic moveouts in the data. This information is more reliable than traveltimes of reflections picked directly from the data because many seismic traces are used for computing beam stacks. To estimate interval velocity, I iteratively search for the velocity model that best predicts the events in beam‐stacked data. My estimation method does not require a preliminary picking of the data because it directly maximizes the beam‐stack’s energy at the traveltimes and surface locations predicted by ray tracing. The advantage of this formulation is that detection of the events in the beam‐stacked data can be guided by the imposition of smoothness constraints on the velocity model. The optimization problem of maximizing beam‐stack energy is solved by a gradient algorithm. To compute the derivatives of the objective function with respect to the velocity model, I derive a linear operator that relates perturbations in velocity to the observed changes in the beam‐stack kinematics. The method has been successfully applied to a marine survey for estimating a low‐velocity anomaly. The estimated velocity function correctly predicts the nonhyperbolic moveouts in the data caused by the velocity anomaly.



Sign in / Sign up

Export Citation Format

Share Document