Apparent layering in common‐midpoint stacked images of two‐dimensionally heterogeneous targets

Geophysics ◽  
1990 ◽  
Vol 55 (11) ◽  
pp. 1466-1477 ◽  
Author(s):  
Bruce S. Gibson ◽  
Alan R. Levander

Model studies with finite‐difference synthetic data demonstrate a fundamental spatial bias in the appearance of common‐midpoint (CMP) stacked images. The CMP stack of data recorded over a target having 2-D random variations in velocity shows numerous short reflection segments; similar reflection patterns in field data are often interpreted in terms of 1-D fine‐scale layering. The stacked image appears layered because of enhanced lateral continuity attributable to the well‐known dip filter of the stacking process. The stack filter can be characterized using the formulation of Bolondi et al. (1982). Lateral correlation in the target and its seismic image is quantified with a measure based on the spectral coefficient of coherence. Broadband primary reflectivity (defined as the vertical‐incidence, primaries‐only reflection coefficients of the 2-D target) is often taken as an ideal seismic image. The primary reflectivity section of a 2-D random target, however, shows greater apparent lateral correlation than is present in the random structure. This apparent increase in lateral continuity is attributable to the fact that reflectivity measured from the surface depends on the vertical derivative of velocity but depends on horizontal changes in velocity directly. The dip‐filtering effects of stacking cannot be reversed by poststack migration; the synthetic data demonstrate the necessity of migration before stack or equivalent processing (such as dip moveout correction). A field data example illustrates the effects of CMP stack filtering using lateral coherence functions measured on stacked and unstacked sections.

2010 ◽  
Vol 14 (3) ◽  
pp. 545-556 ◽  
Author(s):  
J. Rings ◽  
J. A. Huisman ◽  
H. Vereecken

Abstract. Coupled hydrogeophysical methods infer hydrological and petrophysical parameters directly from geophysical measurements. Widespread methods do not explicitly recognize uncertainty in parameter estimates. Therefore, we apply a sequential Bayesian framework that provides updates of state, parameters and their uncertainty whenever measurements become available. We have coupled a hydrological and an electrical resistivity tomography (ERT) forward code in a particle filtering framework. First, we analyze a synthetic data set of lysimeter infiltration monitored with ERT. In a second step, we apply the approach to field data measured during an infiltration event on a full-scale dike model. For the synthetic data, the water content distribution and the hydraulic conductivity are accurately estimated after a few time steps. For the field data, hydraulic parameters are successfully estimated from water content measurements made with spatial time domain reflectometry and ERT, and the development of their posterior distributions is shown.


Geophysics ◽  
2019 ◽  
Vol 84 (2) ◽  
pp. N29-N40
Author(s):  
Modeste Irakarama ◽  
Paul Cupillard ◽  
Guillaume Caumon ◽  
Paul Sava ◽  
Jonathan Edwards

Structural interpretation of seismic images can be highly subjective, especially in complex geologic settings. A single seismic image will often support multiple geologically valid interpretations. However, it is usually difficult to determine which of those interpretations are more likely than others. We have referred to this problem as structural model appraisal. We have developed the use of misfit functions to rank and appraise multiple interpretations of a given seismic image. Given a set of possible interpretations, we compute synthetic data for each structural interpretation, and then we compare these synthetic data against observed seismic data; this allows us to assign a data-misfit value to each structural interpretation. Our aim is to find data-misfit functions that enable a ranking of interpretations. To do so, we formalize the problem of appraising structural interpretations using seismic data and we derive a set of conditions to be satisfied by the data-misfit function for a successful appraisal. We investigate vertical seismic profiling (VSP) and surface seismic configurations. An application of the proposed method to a realistic synthetic model shows promising results for appraising structural interpretations using VSP data, provided that the target region is well-illuminated. However, we find appraising structural interpretations using surface seismic data to be more challenging, mainly due to the difficulty of computing phase-shift data misfits.


2019 ◽  
Vol 221 (1) ◽  
pp. 87-96
Author(s):  
S Malecki ◽  
R-U Börner ◽  
K Spitzer

SUMMARY We present a procedure for localizing underground positions using a time-domain inductive electromagnetic (EM) method. The position to be localized is associated with an EM receiver placed inside the Earth. An EM field is generated by one or more transmitters located at known positions at the Earth’s surface. We then invert the EM field data for the receiver positions using a trust-region algorithm. For any given time regime and source–receiver geometry, the propagation of the electromagnetic fields is determined by the electrical conductivity distribution within the Earth. We show that it is sufficient to use a simple 1-D model to recover the receiver positions with reasonable accuracy. Generally, we demonstrate the robustness of the presented approach. Using confidence ellipses and confidence intervals we assess the accuracy of the recovered location data. The proposed method has been extensively tested against synthetic data obtained by numerical experiments. Furthermore, we have successfully carried out a location recovery using field data. The field data were recorded within a borehole in Alberta (Canada) at 101.4 m depth. The recovered location of the borehole receiver differs from the actual location by 0.70 m in the horizontal plane and by 0.82 m in depth.


Geophysics ◽  
2017 ◽  
Vol 82 (5) ◽  
pp. W31-W45 ◽  
Author(s):  
Necati Gülünay

The old technology [Formula: see text]-[Formula: see text] deconvolution stands for [Formula: see text]-[Formula: see text] domain prediction filtering. Early versions of it are known to create signal leakage during their application. There have been recent papers in geophysical publications comparing [Formula: see text]-[Formula: see text] deconvolution results with the new technologies being proposed. These comparisons will be most effective if the best existing [Formula: see text]-[Formula: see text] deconvolution algorithms are used. This paper describes common [Formula: see text]-[Formula: see text] deconvolution algorithms and studies signal leakage occurring during their application on simple models, which will hopefully provide a benchmark for the readers in choosing [Formula: see text]-[Formula: see text] algorithms for comparison. The [Formula: see text]-[Formula: see text] deconvolution algorithms can be classified by their use of data which lead to transient or transient-free matrices and hence windowed or nonwindowed autocorrelations, respectively. They can also be classified by the direction they are predicting: forward design and apply; forward design and apply followed by backward design and apply; forward design and apply followed by application of a conjugated forward filter in the backward direction; and simultaneously forward and backward design and apply, which is known as noncausal filter design. All of the algorithm types mentioned above are tested, and the results of their analysis are provided in this paper on noise free and noisy synthetic data sets: a single dipping event, a single dipping event with a simple amplitude variation with offset, and three dipping events. Finally, the results of applying the selected algorithms on field data are provided.


Geophysics ◽  
2020 ◽  
Vol 85 (5) ◽  
pp. J85-J98
Author(s):  
Shuang Liu ◽  
Xiangyun Hu ◽  
Dalian Zhang ◽  
Bangshun Wei ◽  
Meixia Geng ◽  
...  

Natural remanent magnetization acts as a record of the previous orientations of the earth’s magnetic field, and it is an important feature when studying geologic phenomena. The so-called IDQ curve is used to describe the relationship between the inclination ( I) and declination ( D) of remanent magnetization and the Köenigsberger ratio ( Q). Here, we construct the IDQ curve using data on ground and airborne magnetic anomalies. The curve is devised using modified approaches for estimating the total magnetization direction, e.g., identifying the maximal position of minimal reduced-to-the-pole fields or identifying correlations between total and vertical reduced-to-the-pole field gradients. The method is tested using synthetic data, and the results indicate that the IDQ curve can provide valuable information on the remanent magnetization direction based on available data on the Köenigsberger ratio. Then, the method is used to interpret field data from the Yeshan region in eastern China, where ground anomalies have been produced by igneous rocks, including diorite and basalt, which occur along with magnetite and hematite ore bodies. The IDQ curves for 24 subanomalies are constructed, and these curves indicate two main distribution clusters of remanent magnetization directions corresponding to different structural units of magma intrusion and help identify the lithologies of the magnetic sources in areas covered by Quaternary sediments. The estimated remanent magnetization directions for Cenozoic basalt are consistent with measurements made in paleomagnetism studies. The synthetic and field data indicate that the IDQ curve can be used to efficiently estimate the remanent magnetization direction from a magnetic anomaly, which could help with our understanding of geologic processes in an area.


Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. KS127-KS138 ◽  
Author(s):  
Yujin Liu ◽  
Yue Ma ◽  
Yi Luo

Locating microseismic source positions using seismic energy emitted from hydraulic fracturing is essential for choosing optimal fracking parameters and maximizing the fracturing effects in hydrocarbon exploitation. Interferometric crosscorrelation migration (ICCM) and zero-lag autocorrelation of time-reversal imaging (ATRI) are two important passive seismic source locating approaches that are proposed independently and seem to be substantially different. We have proven that these two methods are theoretically identical and produce very similar images. Moreover, we have developed cross-coherence that uses normalization by the spectral amplitude of each of the traces, rather than crosscorrelation or deconvolution, to improve the ICCM and ATRI methods. The adopted method enhances the spatial resolution of the source images and is particularly effective in the presence of highly variable and strong additive random noise. Synthetic and field data tests verify the equivalence of the conventional ICCM and ATRI and the equivalence of their improved versions. Compared with crosscorrelation- and deconvolution-based source locating methods, our approach shows a high-resolution property and antinoise capability in numerical tests using synthetic data with single and multiple sources, as well as field data.


Geophysics ◽  
2013 ◽  
Vol 78 (4) ◽  
pp. T79-T87 ◽  
Author(s):  
A. Oelke ◽  
D. Alexandrov ◽  
I. Abakumov ◽  
S. Glubokovskikh ◽  
R. Shigapov ◽  
...  

We have analyzed the angle-dependent reflectivity of microseismic wavefields at a hydraulic fracture, which we modeled as an ideal thin fluid layer embedded in an elastic, isotropic solid rock. We derived full analytical solutions for the reflections of an incident P-wave, the P-P and P-S reflection coefficients, as well as for an incident S-wave, and the S-S and S-P reflection coefficients. The rather complex analytical solutions were then approximated and we found that these zero-thickness limit approximations are in good agreement with the linear slip model, representing a fracture at slip contact. We compared the analytical solutions for the P-P reflections with synthetic data that were derived using finite-difference modeling and found that the modeling confirmed our theoretical results. For typical parameters of microseismic monitoring by hydraulic fracturing, e.g., a layer thickness of [Formula: see text] and frequencies of [Formula: see text], the reflection coefficients depend on the Poisson’s ratio. Furthermore, the reflection coefficients of an incident S-wave are remarkably high. Theoretical results suggested that it is feasible to image hydraulic fractures using microseismic events as a source and to solve the inverse problem, that is, to interpret reflection coefficients extracted from microseismic data in terms of reservoir properties.


Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. JM15-JM24 ◽  
Author(s):  
Tomas Naprstek ◽  
Richard S. Smith

When aeromagnetic data are interpolated to make a gridded image, thin linear features can result in “boudinage” or “string of beads” artifacts if the anomalies are at acute angles to the traverse lines. These artifacts are due to the undersampling of these types of features across the flight lines, making it difficult for most interpolation methods to effectively maintain the linear nature of the features without user guidance. The magnetic responses of dikes and dike swarms are typical examples of the type of geologic feature that can cause these artifacts; thus, these features are often difficult to interpret. Many interpretation methods use various enhancements of the gridded data, such as horizontal or vertical derivatives, and these artifacts are often exacerbated by the processing. Therefore, interpolation methods that are free of these artifacts are necessary for advanced interpretation and analysis of thin, linear features. We have developed a new interpolation method that iteratively enhances linear trends across flight lines, ensuring that linear features are evident on the interpolated grid. Using a Taylor derivative expansion and structure tensors allows the method to continually analyze and interpolate data along anisotropic trends, while honoring the original flight line data. We applied this method to synthetic data and field data, which both show improvement over standard bidirectional gridding, minimum curvature, and kriging methods for interpolating thin, linear features at acute angles to the flight lines. These improved results are also apparent in the vertical derivative enhancement of field data. The source code for this method has been made publicly available.


Geophysics ◽  
2009 ◽  
Vol 74 (6) ◽  
pp. WCA199-WCA209 ◽  
Author(s):  
Guojian Shan ◽  
Robert Clapp ◽  
Biondo Biondi

We have extended isotropic plane-wave migration in tilted coordinates to 3D anisotropic media and applied it on a Gulf of Mexico data set. Recorded surface data are transformed to plane-wave data by slant-stack processing in inline and crossline directions. The source plane wave and its corresponding slant-stacked data are extrapolated into the subsurface within a tilted coordinate system whose direction depends on the propagation direction of the plane wave. Images are generated by crosscorrelating these two wavefields. The shot sampling is sparse in the crossline direction, and the source generated by slant stacking is not really a plane-wave source but a phase-encoded source. We have discovered that phase-encoded source migration in tilted coordinates can image steep reflectors, using 2D synthetic data set examples. The field data example shows that 3D plane-wave migration in tilted coordinates can image steeply dipping salt flanks and faults, even though the one-way wave-equation operator is used for wavefield extrapolation.


Geophysics ◽  
1984 ◽  
Vol 49 (4) ◽  
pp. 379-397 ◽  
Author(s):  
Bruce Gibson ◽  
Ken Larner

Predictive deconvolution is commonly applied to seismic data generated with a Vibroseisr® source. Unfortunately, when this process invokes a minimum‐phase assumption, the phase of the resulting trace will not be correct. Nonetheless, spiking deconvolution is an attractive process because it restores attenuated higher frequencies, thus increasing resolution. For detailed stratigraphic analyses, however, it is desirable that the phase of the data be treated properly as well. The most common solution is to apply a phase‐shifting filter that corrects for errors attributable to a zero‐phase source. The phase correction is given by the minimum‐phase spectrum of the correlated Vibroseis wavelet. Because no minimum‐phase spectrum truly exists for this bandlimited wavelet, white noise is added to its amplitude spectrum in order to design the phase‐correction filter. Different levels of white noise, however, produce markedly different results when field data sections are filtered. A simple argument suggests that the amount of white noise used should match that added in designing the (minimum‐phase) spiking deconvolution operator. This choice, however, also produces inconsistent results; field data again show that the phase treatment is sensitive to the amount of added white noise. Synthetic data tests show that the standard phase‐correction procedure breaks down when earth attenuation is severe. Deterministically reducing the earth‐filter effects before deconvolution improved the resulting phase treatment for the synthetic data. After application of the inverse attenuation filter to the field data, however, phase differences again remain for different levels of added white noise. These inconsistencies are attributable to the phase action of spiking deconvolution. This action is dependent upon the shape of the signal spectrum as well as the spectral shape and level of contaminating noise. Thus, in practice the proper treatment of phase in data-dependent processing requires extensive knowledge of the spectral characteristics of both signal and noise. With such knowledge, one could apply deterministic techniques that either eliminate the need for statistical deconvolution or condition the data so as to satisfy better the statistical model assumed in data‐dependent processing.


Sign in / Sign up

Export Citation Format

Share Document