scholarly journals Higher-Resolution Determination of Zero-Offset Common-Reflection-Surface Stack Parameters

2011 ◽  
Vol 2011 ◽  
pp. 1-10 ◽  
Author(s):  
Endrias G. Asgedom ◽  
Leiv J. Gelius ◽  
Martin Tygel

We developed a higher resolution method for the estimation of the three travel-time parameters that are used in the 2D zero-offset, Common-Reflection-Surface stack method. The underlying principle in this method is to replace the coherency measure performed using semblance with that of MUSIC (multiple signal classification) pseudospectrum that utilizes theeigenstructureof the data covariance matrix. The performance of the two parameter estimation techniques (i.e., semblance and MUSIC) was investigated using both synthetic seismic diffraction and reflection data corrupted with white Gaussian noise, as well as a multioffset ground penetrating radar (GPR) field data set. The estimated parameters employing MUSIC were shown to be superior of those from semblance.

Geophysics ◽  
1993 ◽  
Vol 58 (2) ◽  
pp. 240-256 ◽  
Author(s):  
Xiang‐Yang Li ◽  
Stuart Crampin

Most published techniques for analyzing shear‐wave splitting tend to be computing intensive, and make assumptions, such as the orthogonality of the two split shear waves, which are not necessarily correct. We present a fast linear‐transform technique for analyzing shear‐wave splitting in four‐component (two sources/ two receivers) seismic data, which is flexible and widely applicable. We transform the four‐component data by simple linear transforms so that the complicated shear‐wave motion is linearized in a wide variety of circumstances. This allows various attributes to be measured, including the polarizations of faster split shear waves and the time delays between faster and slower split shear waves, as well as allowing the time series of the faster and slower split shear waves to be separated deterministically. In addition, with minimal assumptions, the geophone orientations can be estimated for zero‐offset verticle seismic profiles (VSPs), and the polarizations of the slower split shear waves can be measured for offset VSPs. The time series of the split shear‐waves can be separated before stack for reflection surveys. The technique has been successfully applied to a number of field VSPs and reflection data sets. Applications to a zero‐offset VSP, an offset VSP, and a reflection data set will be presented to illustrate the technique.


2019 ◽  
Vol 219 (3) ◽  
pp. 1876-1884 ◽  
Author(s):  
M Glöckner ◽  
J Walda ◽  
S Dell ◽  
D Gajewski ◽  
J Karstens ◽  
...  

SUMMARY Standard seismic acquisition and processing require appropriate source–receiver offsets. P-cable technology represents the opposite, namely, very short source–receiver offsets at the price of increased spatial and lateral resolution with a high-frequency source. To use this advantage, a processing flow excluding offset information is required. This aim can be achieved with a processing tuned to diffractions because point diffractions scatter the same information in the offset and midpoint direction. Usually, diffractions are small amplitude events and a careful diffraction separation is required as a first step. We suggest the strategy to use a multiparameter stacking operator, for example, common-reflection surface, and stack along the midpoint direction. The obtained kinematic wave-front attributes are used to calculate time-migration velocities. A diffractivity map serves as a filter to refine the velocities. This strategy is applied to a 3-D P-cable data set to obtain a time-migrated image.


Geophysics ◽  
2004 ◽  
Vol 69 (1) ◽  
pp. 265-274 ◽  
Author(s):  
Eric Duveneck

Kinematic information for constructing velocity models can be extracted in a robust way from seismic prestack data with the common‐reflection‐surface (CRS) stack. This data‐driven process results, in addition to a simulated zero‐offset section, in a number of wavefront attributes—wavefront curvatures and normal ray emergence angles—associated with each simulated zero‐offset sample. A tomographic inversion method is presented that uses this kinematic information to determine smooth, laterally heterogeneous, isotropic subsurface velocity models for depth imaging. The input for the inversion consists of wavefront attributes picked at a number of locations in the simulated zero‐offset section. The smooth velocity model is described by B‐splines. An optimum model is found iteratively by minimizing the misfit between the picked data and the corresponding modeled values. The required forward‐modeled quantities are obtained during each iteration by dynamic ray tracing along normal rays pertaining to the input data points. Fréchet derivatives for the tomographic matrix are calculated by ray perturbation theory. The inversion procedure is demonstrated on a 2D synthetic prestack data set.


Geophysics ◽  
1997 ◽  
Vol 62 (3) ◽  
pp. 884-894 ◽  
Author(s):  
Weijian Mao ◽  
Graham W. Stuart

A multiphase tomographic algorithm is presented that allows 2-D and 3-D slowness (inverse of velocity) and variable reflector depth to be reconstructed simultaneously from both transmission and reflection traveltimes. We analyze the ambiguity in the determination of velocity and depth in transmission and reflection data and realize that depth perturbation is more sensitive to reflection traveltime anomalies than slowness perturbation, whereas the reverse is true of transmission traveltime anomalies. Because of the constraints on velocity and depth provided by the different wave types, this algorithm reduces the ambiguity substantially between velocity and depth prevalent in reflection tomography and also avoids the undetermined problem in transmission tomography. The linearized inversion was undertaken iteratively by decoupling velocity parameters from reflector depths. A rapid 2-D and 3-D ray‐tracing algorithm is used to compute transmission and reflection traveltimes and partial derivatives with respect to slowness and reflector depth. Both depth and velocity are parameterized in terms of cubic B‐spline functions. Synthetic examples indicate the improvement in tomographic results when both transmission and reflection times are included. The method has been applied to a reverse vertical seismic profile (VSP) data set recorded on the British coal measures along a crossed‐linear array. Traveltimes were picked automatically by the simultaneous determination of time delays and stacking weights using a waveform matching technique. The tomographic inversion of the observed reverse VSP images two fault‐zones of lower velocity than the surrounding media. The location of the faults was confirmed by near‐by reflection lines. The technique can be applied to offset VSPs or reverse VSPs and coincident VSP and surface reflection data.


Geophysics ◽  
1999 ◽  
Vol 64 (1) ◽  
pp. 208-221 ◽  
Author(s):  
Tamas Nemeth ◽  
Chengjun Wu ◽  
Gerard T. Schuster

A least‐squares migration algorithm is presented that reduces the migration artifacts (i.e., recording footprint noise) arising from incomplete data. Instead of migrating data with the adjoint of the forward modeling operator, the normal equations are inverted by using a preconditioned linear conjugate gradient scheme that employs regularization. The modeling operator is constructed from an asymptotic acoustic integral equation, and its adjoint is the Kirchhoff migration operator. We tested the performance of the least‐squares migration on synthetic and field data in the cases of limited recording aperture, coarse sampling, and acquisition gaps in the data. Numerical results show that the least‐squares migrated sections are typically more focused than are the corresponding Kirchhoff migrated sections and their reflectivity frequency distributions are closer to those of the true model frequency distribution. Regularization helps attenuate migration artifacts and provides a sharper, better frequency distribution of estimated reflectivity. The least‐squares migrated sections can be used to predict the missing data traces and interpolate and extrapolate them according to the governing modeling equations. Several field data examples are presented. A ground‐penetrating radar data example demonstrates the suppression of the recording footprint noise due to a limited aperture, a large gap, and an undersampled receiver line. In addition, better fault resolution was achieved after applying least‐squares migration to a poststack marine data set. And a reverse vertical seismic profiling example shows that the recording footprint noise due to a coarse receiver interval can be suppressed by least‐squares migration.


Geophysics ◽  
2021 ◽  
pp. 1-49
Author(s):  
Zhiwei Xu ◽  
James Irving ◽  
Yu Liu ◽  
Zhu Peimin ◽  
Klaus Holliger

We present a stochastic inversion procedure for common-offset ground-penetrating radar (GPR) reflection measurements. Stochastic realizations of subsurface properties that offer an acceptable fit to GPR data are generated via simulated annealing optimization. The realizations are conditioned to borehole porosity measurements available along the GPR profile, or equivalent measurements of another petrophysical property that can be related to the dielectric permittivity, as well as to geostatistical parameters derived from the borehole logs and the processed GPR image. Validation of our inversion procedure is performed on a pertinent synthetic data set and indicates that the proposed method is capable of reliably recovering strongly heterogeneous porosity structures associated with surficial alluvial aquifers. This finding is largely corroborated through application of the methodology to field measurements from the Boise Hydrogeophysical Research Site near Boise, Idaho, USA.


Geophysics ◽  
2001 ◽  
Vol 66 (2) ◽  
pp. 540-550 ◽  
Author(s):  
Albane C. Saintenoy ◽  
Albert Tarantola

The three electromagnetic properties appearing in Maxwell’s equations are dielectric permittivity, electrical conductivity and magnetic permeability. The study of point diffractors in a homogeneous, isotropic, linear medium suggests the use of logarithms to describe the variations of electromagnetic properties in the earth. A small anomaly in electrical properties (permittivity and conductivity) responds to an incident electromagnetic field as an electric dipole, whereas a small anomaly in the magnetic property responds as a magnetic dipole. Neither property variation can be neglected without justification. Considering radiation patterns of the different diffracting points, diagnostic interpretation of electric and magnetic variations is theoretically feasible but is not an easy task using ground‐penetrating radar. However, using an effective electromagnetic impedance and an effective electromagnetic velocity to describe a medium, the radiation patterns of a small anomaly behave completely differently with source‐receiver offset. Zero‐offset reflection data give a direct image of impedance variations while large‐offset reflection data contain information on velocity variations.


Geophysics ◽  
2001 ◽  
Vol 66 (1) ◽  
pp. 97-109 ◽  
Author(s):  
Rainer Jäger ◽  
Jürgen Mann ◽  
German Höcht ◽  
Peter Hubral

The common‐reflection‐surface stack provides a zero‐offset simulation from seismic multicoverage reflection data. Whereas conventional reflection imaging methods (e.g. the NMO/dip moveout/stack or prestack migration) require a sufficiently accurate macrovelocity model to yield appropriate results, the common‐reflection‐surface (CRS) stack does not depend on a macrovelocity model. We apply the CRS stack to a 2-D synthetic seismic multicoverage dataset. We show that it not only provides a high‐quality simulated zero‐offset section but also three important kinematic wavefield attribute sections, which can be used to derive the 2-D macrovelocity model. We compare the multicoverage‐data‐derived attributes with the model‐derived attributes computed by forward modeling. We thus confirm the validity of the theory and of the data‐derived attributes. For 2-D acquisition, the CRS stack leads to a stacking surface depending on three search parameters. The optimum stacking surface needs to be determined for each point of the simulated zero‐offset section. For a given primary reflection, these are the emergence angle α of the zero‐offset ray, as well as two radii of wavefront curvatures [Formula: see text] and [Formula: see text]. They all are associated with two hypothetical waves: the so‐called normal wave and the normal‐incidence‐point wave. We also address the problem of determining an optimal parameter triplet (α, [Formula: see text], [Formula: see text]) in order to construct the sample value (i.e., the CRS stack value) for each point in the desired simulated zero‐offset section. This optimal triplet is expected to determine for each point the best stacking surface that can be fitted to the multicoverage primary reflection events. To make the CRS stack attractive in terms of computational costs, a suitable strategy is described to determine the optimal parameter triplets for all points of the simulated zero‐offset section. For the implementation of the CRS stack, we make use of the hyperbolic second‐order Taylor expansion of the stacking surface. This representation is not only suitable to handle irregular multicoverage acquisition geometries but also enables us to introduce simple and efficient search strategies for the parameter triple. In specific subsets of the multicoverage data (e.g., in the common‐midpoint gathers or the zero‐offset section), the chosen representation only depends on one or two independent parameters, respectively.


Geophysics ◽  
2017 ◽  
Vol 82 (2) ◽  
pp. R65-R73 ◽  
Author(s):  
Alexander Bauer ◽  
Benjamin Schwarz ◽  
Dirk Gajewski

Wavefront tomography is known to be an efficient and stable approach for velocity inversion that does not require accurate starting models and does not interact directly with the prestack data. Instead, the original data are transformed to physically meaningful wavefront attribute fields. These can be automatically estimated using local-coherence analysis by means of the common-reflection-surface (CRS) stack, which has been shown to be a powerful tool for data analysis and enhancement. In addition, the zero-offset wavefront attributes acquired during the CRS stack can be used for sophisticated subsequent processes such as wavefield characterization and separation. Whereas in previous works, wavefront tomography has been applied mainly to reflection data, resulting in smooth velocity models suitable for migration of targets with moderately complex overburden, we have emphasized using the diffracted contributions in the data for velocity inversion. By means of simple synthetic examples, we reveal the potential of diffractions for velocity inversion. On industrial field data, we suggest a joint inversion based on reflected and diffracted contributions of the measured wavefield, which confirms the general finding that diffraction-based wavefront tomography can help to increase the resolution of the velocity models. Concluding our work, we compare the quality of a reverse time migrated result using the estimated velocity model with the result based on the inversion of reflections, which reveals an improved imaging potential for a complex salt geometry.


Author(s):  
Peter M. Fischer ◽  
Patrik Klingborg ◽  
Fanny Kärfve ◽  
Fredrika Kärfve ◽  
C. Hagberg ◽  
...  

Determination of the complete occupational sequence of the site, including investigation of pre-12th century levels which were thoroughly studied by P. Åström since the 1970s, is the main task of the planned project. During the course of the expedition (NSCE11) in spring 2010 a ground-penetrating radar survey (GPR) was carried out at Dromolaxia Vizatzia/Hala Sultan Tekke in Area 6, leading to the discovery of a large Late Cypriote complex. The compound is bordered to the north by a substantial wall, against which nine rooms (so far) could be exposed. Two occupational phases have been verified but there are indications of a third. The suggested functions of the various structures of the most recent phase are: living, working, storage and administration spaces. The rich find contexts point to the production of textiles and metal objects, and the locally produced pottery is generally of a high quality. There are also many imports, mainly from the Mycenaean sphere of culture. The locally produced vessels from Phase 2 include the “Creature krater” which is a masterpiece of a high artistic standard. Another piece of elevated artistry is the piece of a “Warrior vase”.


Sign in / Sign up

Export Citation Format

Share Document