Euclid and the art of wavelet estimation, Part II: Robust algorithm and field‐data examples

Geophysics ◽  
1997 ◽  
Vol 62 (6) ◽  
pp. 1939-1946 ◽  
Author(s):  
Eike Rietsch
Keyword(s):  
Field Data ◽  
Synthetic Data ◽  
Noise Resistance ◽  
Zero Eigenvalue ◽  
Seismic Line ◽  
Robust Algorithm ◽  
Seismic Wavelet ◽  
Time Zones ◽  
The Common ◽  
Common Signal

In this second part of a two‐part work, a more robust algorithm is derived and used for the estimation of the seismic wavelet as the common signal of two or more seismic traces. It is based on the properties of the eigenvectors with zero eigenvalue of a matrix derived in the first part, whose elements are the samples of the autocorrelation functions and crosscorrelation functions of these seismic traces for a number of lags. The noise resistance of this algorithm is illustrated by means of a synthetic‐data example and then demonstrated on field data. In one field‐data example, the so‐called Euclid wavelet is compared with one derived deterministically by means of an impedance log. The other example relates three quite different Euclid wavelets determined in three different time zones on a seismic line to one another by showing that their differences can be explained by absorption.

scholarly journals Coupled hydrogeophysical parameter estimation using a sequential Bayesian approach

2010 ◽  
Vol 14 (3) ◽  
pp. 545-556 ◽  
Author(s):  
J. Rings ◽  
J. A. Huisman ◽  
H. Vereecken
Keyword(s):  
Water Content ◽  
Particle Filtering ◽  
Field Data ◽  
Synthetic Data ◽  
Time Domain Reflectometry ◽  
Second Step ◽  
Parameter Estimates ◽  
Posterior Distributions ◽  
Data Set ◽  
Resistivity Tomography

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.

scholarly journals Fractional integration of seismic wavelets in anelastic media to recover multiscale properties of impedance discontinuities

Geophysics ◽  
2018 ◽  
Vol 83 (2) ◽  
pp. V61-V71 ◽  
Author(s):  
Stephan Ker ◽  
Yves Le Gonidec
Keyword(s):  
Wavelet Transform ◽  
Fractional Integration ◽  
Synthetic Data ◽  
Seismic Source ◽  
Seismic Attributes ◽  
Data Set ◽  
Shaping Filter ◽  
Gaussian Derivative ◽  
Intrinsic Attenuation ◽  
Seismic Wavelet

Multiscale seismic attributes based on wavelet transform properties have recently been introduced and successfully applied to identify the geometry of a complex seismic reflector in an elastic medium. We extend this quantitative approach to anelastic media where intrinsic attenuation modifies the seismic attributes and thus requires a specific processing to retrieve them properly. The method assumes an attenuation linearly dependent with the seismic wave frequency and a seismic source wavelet approximated with a Gaussian derivative function (GDF). We highlight a quasi-conservation of the Gaussian character of the wavelet during its propagation. We found that this shape can be accurately modeled by a GDF characterized by a fractional integration and a frequency shift of the seismic source, and we establish the relationship between these wavelet parameters and [Formula: see text]. Based on this seismic wavelet modeling, we design a time-varying shaping filter that enables making constant the shape of the wavelet allowing retrieval of the wavelet transform properties. Introduced with a homogeneous step-like reflector, the method is first applied on a thin-bed reflector and then on a more realistic synthetic data set based on an in situ acoustic impedance sequence and a high-resolution seismic source. The results clearly highlight the efficiency of the method in accurately restoring the multiscale seismic attributes of complex seismic reflectors in anelastic media by the use of broadband seismic sources.

scholarly journals Zika mosquito vectors: the jury is still out

F1000Research ◽  
2016 ◽  
Vol 5 ◽  
pp. 2546 ◽  
Author(s):  
Walter S. Leal
Keyword(s):  
Aedes Aegypti ◽  
Yellow Fever ◽  
Culex Quinquefasciatus ◽  
Zika Virus ◽  
Field Data ◽  
Vector Competence ◽  
Sexual Transmission ◽  
Mosquito Vectors ◽  
Yellow Fever Mosquito ◽  
The Common

After a 40-year hiatus, the International Congress of Entomology (ICE 2016) convened in Orlando, Florida (September 25-30, 2016). One of the symposia at ICE 2016, the Zika Symposium, covered multiple aspects of the Zika epidemic, including epidemiology, sexual transmission, genetic tools for reducing transmission, and particularly vector competence. While there was a consensus among participants that the yellow fever mosquito, Aedes aegypti, is a vector of the Zika virus, there is growing evidence indicating that the range of mosquito vectors might be wider than anticipated. In particular, three independent groups from Canada, China, and Brazil presented and discussed laboratory and field data strongly suggesting that the southern house mosquito, Culex quinquefasciatus, also known as the common mosquito, is highly likely to be a vector in certain environments.

A new method for absolute underground positioning based on transient electromagnetics

2019 ◽  
Vol 221 (1) ◽  
pp. 87-96
Author(s):  
S Malecki ◽  
R-U Börner ◽  
K Spitzer
Keyword(s):  
Field Data ◽  
Synthetic Data ◽  
Trust Region ◽  
Trust Region Algorithm ◽  
Location Data ◽  
The Earth ◽  
Time Regime ◽  
Em Method ◽  
Em Field ◽  
Confidence Ellipses

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.

Signal leakage in f-x deconvolution algorithms

Geophysics ◽  
2017 ◽  
Vol 82 (5) ◽  
pp. W31-W45 ◽  
Author(s):  
Necati Gülünay
Keyword(s):  
Field Data ◽  
New Technologies ◽  
Filter Design ◽  
Synthetic Data ◽  
Data Sets ◽  
Amplitude Variation ◽  
Amplitude Variation With Offset ◽  
Backward Design ◽  
Use Of Data ◽  
Domain Prediction

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.

Near-surface scattering phenomena and implications for tunnel detection

2015 ◽  
Vol 3 (1) ◽  
pp. SF43-SF54 ◽  
Author(s):  
Shelby L. Peterie ◽  
Richard D. Miller
Keyword(s):  
Synthetic Data ◽  
P Wave ◽  
S Wave ◽  
Seismic Line ◽  
Data Set ◽  
Near Surface ◽  
Tunnel Detection ◽  
Impedance Contrast ◽  
Diffraction Imaging ◽  
Imaging Algorithms

Tunnel locations are accurately interpreted from diffraction sections of focused mode converted P- to S-wave diffractions from a perpendicular tunnel and P-wave diffractions from a nonperpendicular (oblique) tunnel. Near-surface tunnels are ideal candidates for diffraction imaging due to their small size relative to the seismic wavelength and large acoustic impedance contrast at the tunnel interface. Diffraction imaging algorithms generally assume that the velocities of the primary wave and the diffracted wave are approximately equal, and that the diffraction apex is recorded directly above the scatterpoint. Scattering phenomena from shallow tunnels with kinematic properties that violate these assumptions were observed in one field data set and one synthetic data set. We developed the traveltime equations for mode-converted and oblique diffractions and demonstrated a diffraction imaging algorithm designed for the roll-along style of acquisition. Potential processing and interpretation pitfalls specific to these diffraction types were identified. Based on our observations, recommendations were made to recognize and image mode-converted and oblique diffractions and accurately interpret tunnel depth, horizontal location, and azimuth with respect to the seismic line.

The IDQ curve: A tool for evaluating the direction of remanent magnetization from magnetic anomalies

Geophysics ◽  
2020 ◽  
Vol 85 (5) ◽  
pp. J85-J98
Author(s):  
Shuang Liu ◽  
Xiangyun Hu ◽  
Dalian Zhang ◽  
Bangshun Wei ◽  
Meixia Geng ◽  
...  
Keyword(s):  
Remanent Magnetization ◽  
Field Data ◽  
Eastern China ◽  
Synthetic Data ◽  
Magnetic Anomalies ◽  
Natural Remanent Magnetization ◽  
Magma Intrusion ◽  
Magnetization Direction ◽  
Ore Bodies ◽  
Using Data

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.

Source location with cross-coherence migration

Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. KS127-KS138 ◽  
Author(s):  
Yujin Liu ◽  
Yue Ma ◽  
Yi Luo
Keyword(s):  
Field Data ◽  
Random Noise ◽  
Synthetic Data ◽  
Seismic Energy ◽  
Seismic Source ◽  
Spectral Amplitude ◽  
Multiple Sources ◽  
Passive Seismic ◽  
Similar Images ◽  
Resolution Property

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.

scholarly journals A new method for interpolating linear features in aeromagnetic data

Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. JM15-JM24 ◽  
Author(s):  
Tomas Naprstek ◽  
Richard S. Smith
Keyword(s):  
Field Data ◽  
Interpolation Method ◽  
Source Code ◽  
Synthetic Data ◽  
Aeromagnetic Data ◽  
Derivative Expansion ◽  
Linear Features ◽  
Interpolation Methods ◽  
Minimum Curvature ◽  
Linear Nature

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.

Sign in / Sign up

Export Citation Format

Share Document