Deep learning spectral enhancement considering features of seismic field data

Geophysics ◽  
2021 ◽  
pp. 1-60
Author(s):  
Yonggyu Choi ◽  
Yeonghwa Jo ◽  
Soon Jee Seol ◽  
Joongmoo Byun ◽  
Young Kim
Keyword(s):  
Machine Learning ◽  
Seismic Data ◽  
Input Data ◽  
Field Data ◽  
A Priori ◽  
Resolution Enhancement ◽  
Training Data ◽  
A Priori Information ◽  
Spectral Enhancement ◽  
Priori Information

The resolution of seismic data dictates the ability to identify individual features or details in a given image, and the temporal (vertical) resolution is a function of the frequency content of a signal. To improve thin-bed resolution, broadening of the frequency spectrum is required; this has been one of the major objectives in seismic data processing. Recently, many researchers have proposed machine learning based resolution enhancement and showed their applicability. However, since the performance of machine learning depends on what the model has learned, output from training data with features different from the target field data may be poor. Thus, we present a machine learning based spectral enhancement technique considering features of seismic field data. We used a convolutional U-Net model, which preserves the temporal connectivity and resolution of the input data, and generated numerous synthetic input traces and their corresponding spectrally broadened traces for training the model. A priori information from field data, such as the estimated source wavelet and reflectivity distribution, was considered when generating the input data for complementing the field features. Using synthetic tests and field post-stack seismic data examples, we showed that the trained model with a priori information outperforms the models trained without a priori information in terms of the accuracy of enhanced signals. In addition, our new spectral enhancing method was verified through the application to the high-cut filtered data and its promising features were presented through the comparison with well log data.

Resolution Enhancement of Composite Spectra with Fractal Noise in Derivative Spectrometry

2000 ◽  
Vol 54 (5) ◽  
pp. 721-730 ◽  
Author(s):  
S. S. Kharintsev ◽  
D. I. Kamalova ◽  
M. Kh. Salakhov
Keyword(s):  
A Priori ◽  
Resolution Enhancement ◽  
Numerical Differentiation ◽  
A Priori Information ◽  
Experimental Ir ◽  
Composite Spectra ◽  
Self Similar ◽  
Priori Information ◽  
Fractal Noise ◽  
Derivative Spectrometry

The problem of improving the resolution of composite spectra with statistically self-similar (fractal) noise is considered within the framework of derivative spectrometry. An algorithm of the numerical differentiation of an arbitrary (including fractional) order of spectra is produced by the statistical regularization method taking into account a priori information on statistical properties of the fractal noise. Fractal noise is analyzed in terms of the statistical Hurst method. The efficiency and expedience of this algorithm are exemplified by treating simulated and experimental IR spectra.

scholarly journals A method of combining coherence-constrained sparse coding and dictionary learning for denoising

Geophysics ◽  
2017 ◽  
Vol 82 (3) ◽  
pp. V137-V148 ◽  
Author(s):  
Pierre Turquais ◽  
Endrias G. Asgedom ◽  
Walter Söllner
Keyword(s):  
Dictionary Learning ◽  
Seismic Data ◽  
Random Noise ◽  
A Priori ◽  
Seismic Signal ◽  
Noise Variance ◽  
A Priori Information ◽  
Seismic Section ◽  
Learned Dictionary ◽  
Priori Information

We have addressed the seismic data denoising problem, in which the noise is random and has an unknown spatiotemporally varying variance. In seismic data processing, random noise is often attenuated using transform-based methods. The success of these methods in denoising depends on the ability of the transform to efficiently describe the signal features in the data. Fixed transforms (e.g., wavelets, curvelets) do not adapt to the data and might fail to efficiently describe complex morphologies in the seismic data. Alternatively, dictionary learning methods adapt to the local morphology of the data and provide state-of-the-art denoising results. However, conventional denoising by dictionary learning requires a priori information on the noise variance, and it encounters difficulties when applied for denoising seismic data in which the noise variance is varying in space or time. We have developed a coherence-constrained dictionary learning (CDL) method for denoising that does not require any a priori information related to the signal or noise. To denoise a given window of a seismic section using CDL, overlapping small 2D patches are extracted and a dictionary of patch-sized signals is trained to learn the elementary features embedded in the seismic signal. For each patch, using the learned dictionary, a sparse optimization problem is solved, and a sparse approximation of the patch is computed to attenuate the random noise. Unlike conventional dictionary learning, the sparsity of the approximation is constrained based on coherence such that it does not need a priori noise variance or signal sparsity information and is still optimal to filter out Gaussian random noise. The denoising performance of the CDL method is validated using synthetic and field data examples, and it is compared with the K-SVD and FX-Decon denoising. We found that CDL gives better denoising results than K-SVD and FX-Decon for removing noise when the variance varies in space or time.

scholarly journals Inverse Problems and Artificial Intelligence

2021 ◽  
pp. 33-43
Author(s):  
С.И. Кабанихин
Keyword(s):  
Artificial Intelligence ◽  
Machine Learning ◽  
A Priori ◽  
Stability And Convergence ◽  
Computational Algorithms ◽  
Conditional Stability ◽  
A Priori Information ◽  
Sensitivity Data ◽  
Priori Information ◽  
The Relationship

В данной работе приведен анализ взаимосвязей теории обратных и некорректных задач и математических аспектов искусственного интеллекта. Показано, что при анализе вычислительных алгоритмов, которые условно можно отнести к вычислительному искусственному интеллекту (машинное обучение, природоподобные алгоритмы, методы анализа и обработки данных), возможно, а подчас и необходимо, использовать результаты и подходы, развитые в теории и численных методах решения обратных и некорректных задач, такие как регуляризация, условная устойчивость и сходимость, использование априорной информации, идентифицируемость, чувствительность, усвоение данных. This paper analyzes the relationship between the theory of inverse and incorrect problems and the mathematical aspects of artificial intelligence. It is shown that computational algorithms that can be categorized as computational artificial intelligence (machine learning, nature-like algorithms, data analysis and processing) can or should be analyzed with the approaches developed for the theory and numerical methods for solving inverse and incorrect problems. They are regularization, conditional stability and convergence, the use of a priori information, identifiability, sensitivity, data assimilation.

On data-independent multicomponent interpolators and the use of priors for optimal reconstruction and 3D up/down separation of pressure wavefields

Geophysics ◽  
2010 ◽  
Vol 75 (6) ◽  
pp. WB39-WB51 ◽  
Author(s):  
Kemal Özdemir ◽  
Ali Özbek ◽  
Dirk-Jan van Manen ◽  
Massimiliano Vassallo
Keyword(s):  
Particle Motion ◽  
Seismic Data ◽  
A Priori ◽  
Sparse Sampling ◽  
Effective Bandwidth ◽  
A Priori Information ◽  
Sampling Theorems ◽  
Multichannel Sampling ◽  
Priori Information ◽  
Seismic Wavefield

In marine acquisition, the interference between the upgoing and downgoing wavefields introduces a receiver ghost which reduces the effective bandwidth of the seismic wavefield. A two-component streamer provides means for removing the receiver ghost by measuring pressure and vertical particle velocity. However, due to nonuniform and relatively sparse sampling in the crossline direction, the seismic data are usually severely aliased in the crossline direction and the deghosting may not be feasible in a true 3D sense. A true multicomponent streamer measures all components of the particle motion wavefield in addition to the pressure wavefield. This enables solving the 3D deghosting and crossline reconstruction problems simultaneously, without making assumptions on the wavefield or the subsurface. We havedeveloped two data-independent algorithms suited for multicomponent acquisition. The first algorithm reconstructs the total pressure wavefield in the crossline direction by using the pressure and the crossline component of particle motion simultaneously. The second algorithm reconstructs the upgoing pressure wavefield by using the pressure, the crossline, and the vertical components of particle motion simultaneously. Both algorithms are optimal in the minimum-mean-squares-error sense and are ideally suited for a small number of irregularly spaced samples, as is common in towed marine acquisition. We find that by using the spectrum of the wavefield as a priori information, these algorithms have the potential to overcome higher-order aliasing than what is predicted by multichannel sampling theorems. Such a priori information can be extracted from an unaliased portion of the seismic data in novel and robust manners.

A constrained global inversion method using an overparameterized scheme: Application to poststack seismic data

Geophysics ◽  
2001 ◽  
Vol 66 (2) ◽  
pp. 613-626 ◽  
Author(s):  
Xin‐Quan Ma
Keyword(s):  
Simulated Annealing ◽  
Seismic Data ◽  
A Priori ◽  
Constant Thickness ◽  
Inversion Method ◽  
Model Parameters ◽  
A Priori Information ◽  
Initial Model ◽  
Seismic Waveform ◽  
Priori Information

A global optimization algorithm using simulated annealing has advantages over local optimization approaches in that it can escape from being trapped in local minima and it does not require a good initial model and function derivatives to find a global minimum. It is therefore more attractive and suitable for seismic waveform inversion. I adopt an improved version of a simulated annealing algorithm to invert simultaneously for acoustic impedance and layer interfaces from poststack seismic data. The earth’s subsurface is overparameterized by a series of microlayers with constant thickness in two‐way traveltime. The algorithm is constrained using the low‐frequency impedance trend and has been made computationally more efficient using this a priori information as an initial model. A search bound of each parameter, derived directly from the a priori information, reduces the nonuniqueness problem. Application of this technique to synthetic and field data examples helps one recover the true model parameters and reveals good continuity of estimated impedance across a seismic section. This approach has the capability of revealing the high‐resolution detail needed for reservoir characterization when a reliable migrated image is available with good well ties.

What can machine learning do for seismic data processing? An interpolation application

Geophysics ◽  
2017 ◽  
Vol 82 (3) ◽  
pp. V163-V177 ◽  
Author(s):  
Yongna Jia ◽  
Jianwei Ma
Keyword(s):  
Machine Learning ◽  
Seismic Data ◽  
Input Data ◽  
Interpolation Method ◽  
Window Size ◽  
Tight Frame ◽  
Training Data ◽  
Low Rank ◽  
Support Vector ◽  
Data Sets

Machine learning (ML) systems can automatically mine data sets for hidden features or relationships. Recently, ML methods have become increasingly used within many scientific fields. We have evaluated common applications of ML, and then we developed a novel method based on the classic ML method of support vector regression (SVR) for reconstructing seismic data from under-sampled or missing traces. First, the SVR method mines a continuous regression hyperplane from training data that indicates the hidden relationship between input data with missing traces and output completed data, and then it interpolates missing seismic traces for other input data by using the learned hyperplane. The key idea of our new ML method is significantly different from that of many previous interpolation methods. Our method depends on the characteristics of the training data, rather than the assumptions of linear events, sparsity, or low rank. Therefore, it can break out the previous assumptions or constraints and show universality to different data sets. In addition, our method dramatically reduces the manual workload; for example, it allows users to avoid selecting the window size parameters, as is required for methods based on the assumption of linear events. The ML method facilitates intelligent interpolation between data sets with similar geomorphological structures, which can significantly reduce costs in engineering applications. Furthermore, we combine a sparse transform called the data-driven tight frame (so-called compressed learning) with the SVR method to improve the training performance, in which the training is implemented in a sparse coefficient domain rather than in the data domain. Numerical experiments show the competitive performance of our method in comparison with the traditional [Formula: see text]-[Formula: see text] interpolation method.

2-D Bayesian deconvolution

Geophysics ◽  
1991 ◽  
Vol 56 (12) ◽  
pp. 2008-2018 ◽  
Author(s):  
Marc Lavielle
Keyword(s):  
Objective Function ◽  
Seismic Data ◽  
Layer Structure ◽  
A Priori ◽  
A Priori Information ◽  
The Earth ◽  
Markov Random ◽  
Posteriori Distribution ◽  
A Posteriori Distribution ◽  
Priori Information

Inverse problems can be solved in different ways. One way is to define natural criteria of good recovery and build an objective function to be minimized. If, instead, we prefer a Bayesian approach, inversion can be formulated as an estimation problem where a priori information is introduced and the a posteriori distribution of the unobserved variables is maximized. When this distribution is a Gibbs distribution, these two methods are equivalent. Furthermore, global optimization of the objective function can be performed with a Monte Carlo technique, in spite of the presence of numerous local minima. Application to multitrace deconvolution is proposed. In traditional 1-D deconvolution, a set of uni‐dimensional processes models the seismic data, while a Markov random field is used for 2-D deconvolution. In fact, the introduction of a neighborhood system permits one to model the layer structure that exists in the earth and to obtain solutions that present lateral coherency. Moreover, optimization of an appropriated objective function by simulated annealing allows one to control the fit with the input data as well as the spatial distribution of the reflectors. Extension to 3-D deconvolution is straightforward.

Structurally tailored 3D anisotropic controlled-source electromagnetic resistivity inversion with cross-gradient criterion and simultaneous model calibration

Geophysics ◽  
2019 ◽  
Vol 84 (6) ◽  
pp. E387-E402 ◽  
Author(s):  
Max A. Meju ◽  
Randall L. Mackie ◽  
Federico Miorelli ◽  
Ahmad Shahir Saleh ◽  
Roger V. Miller
Keyword(s):  
Field Data ◽  
A Priori ◽  
Calibration Data ◽  
A Priori Information ◽  
Data Set ◽  
Low Contrast ◽  
Controlled Source ◽  
Low Resistivity ◽  
Band Limited ◽  
Priori Information

Geologic interpretation of 3D anisotropic resistivity models from conventional marine controlled-source electromagnetic (CSEM) data inversion faces difficulties in low-resistivity contrast sediments and structurally complex environments that typify the new frontiers for hydrocarbon exploration. Currently, the typically reconstructed horizontal resistivity [Formula: see text] and vertical resistivity [Formula: see text] models often have conflicting depth structures that are difficult to explain in terms of subsurface geology, and the resulting resistivities may not be close to the true formation resistivities required for estimating reservoir parameters. We have investigated the concept that an objective geologically oriented or structurally tailored inversion can be achieved by requiring that the cross-product of the gradient of horizontal resistivity and the gradient of the vertical resistivity is equal to zero at significant geologic boundaries. We incorporate this boundary-shape criterion in our 3D inverse problem formulations, implemented within nonlinear model-space and conjugate-gradient contexts, for cases in which a priori calibration data from wells and/or seismically derived subsurface boundaries are available and for cases in which these are lacking. The resulting fit-for-purpose solutions serve to better analyze the peculiarity of a given data set. We applied these algorithms to synthetic and field CSEM data sets representing a fold-thrust environment with low-resistivity and low-contrast sediments. The resulting [Formula: see text] and [Formula: see text] models from cross-gradient joint inversion of synthetic data of appropriate frequency bandwidth without a priori information are structurally similar and consistent with the test models, whereas those from the inversions of band-limited field data are consistent with the available seismic and resistivity well-log data. This particular approach will thus be useful for lithologic correlation in frontier regions with limited a priori information using broadband CSEM data. For these band-limited field data, we found that the anisotropic bulk resistivities of the low-contrast sediments are better determined by incorporating a priori calibration data from triaxial resistivity logs and seismic horizons.

Improvements to elastic full-waveform inversion using cross-gradient constraints

Geophysics ◽  
2018 ◽  
Vol 83 (2) ◽  
pp. R105-R115 ◽  
Author(s):  
Edgar Manukyan ◽  
Hansruedi Maurer ◽  
André Nuber
Keyword(s):  
Field Data ◽  
A Priori ◽  
Waveform Inversion ◽  
Structural Similarity ◽  
Model Parameters ◽  
Full Waveform Inversion ◽  
A Priori Information ◽  
Individual Parameter ◽  
Full Waveform ◽  
Priori Information

Seismic full-waveform inversion (FWI) is potentially a powerful method for obtaining high-resolution subsurface images, but the results are often distorted by nonlinear effects and parameter trade-offs. Such distortions can be particularly severe in the case of multiparameter FWI, such as elastic FWI, in which inversion is performed simultaneously for P- and S-wave velocities and density. The problem can be alleviated by adding constraints in the form of plausible a priori information. A usually well-justified constraint includes the structural similarity of different model parameters; i.e., an anomalous body likely exhibits variations in all elastic properties, although their magnitudes may be different. To consider such types of a priori information, we have developed a structurally constrained elastic FWI, which is based on minimization of the cross products of gradients of different model parameters. Our synthetic 2D experiments show that structurally constrained FWI can significantly improve model reconstruction. It is also demonstrated that our approach still leads to improved results, even when the structural similarity between the individual parameter types is not exactly met. Inversions of field data show that in comparison to conventional FWI, structurally constrained FWI is able to match the field data equally well while requiring less structural complexity of the subsurface.

3-D HYBRID IMAGING BASED ON GRAVITY MIGRATION AND REGULARIZED FOCUSING INVERSION TO PREDICT THE POYANG BASIN INTERFACE

Geophysics ◽  
2021 ◽  
pp. 1-46
Author(s):  
Zhengwei Xu ◽  
Rui Wang ◽  
Wei Xiong ◽  
Jian Wang ◽  
Dian Wang
Keyword(s):  
High Resolution ◽  
Density Distribution ◽  
Seismic Data ◽  
Sedimentary Basin ◽  
A Priori ◽  
Sedimentary Basins ◽  
Density Contrast ◽  
A Priori Information ◽  
Spatial Direction ◽  
Priori Information

Describing and understanding the basement relief of sedimentary basins is vital for oil and gas exploration. The traditional method to map an interface in each spatial direction is based on three-dimensional (3D) modeling of gravity Bouguer anomalies with variable lateral and vertical density contrasts using a priori information derived from other types of geoscience datasets as constraints (e.g., well and/or seismic data). However, in the pre-exploration stage, vertical gravity, gz, which is sometimes the only available geophysical data, are typically used to recover smooth density contrast distributions under a generic set of constraints. Apparently, the use of the gz component is not sufficient to produce geologically reasonable interpretations with high resolution. To address this, we developed a novel process of hybrid inversion, combining gravity migration and inversion using the same gz dataset, to distinguish the complicated interface between basement and sedimentary basin rocks from a full-space inverted density distribution volume. First, a 3D-migrated model delineating the basic sedimentary basin structure was derived using a focusing gravity iterative migration method, where a priori information is not necessary. Subsequently, under the framework of the regularized focusing conjugate inversion algorithm, a high-resolution density contrast model was inverted for the delineation of the basement boundary by integrating the 3D-migrated density model as a priori information. We examined the method using one synthetic example and a field data case, of which a transformed resolution density matrix was developed from logarithmic space to qualitatively evaluate the practical resolutions. The high resolution of density distribution of Cretaceous basement with clear interface was achieved and verified by limited seismic data and strata markers in limited wells.

Sign in / Sign up

Export Citation Format

Share Document