Accurate removal of surface-related multiples remains a challenge in shallow-water cases. One reason is that the success of surface-related multiple estimation (SRME)-related algorithms is sensitive to the quality of the near-offset reconstruction. When it comes to a larger missing gap and a shallower water bottom, the state-of-the-art near-offset gap construction method - the parabolic Radon transform - fails to provide reliable recovery of the shallow reflections due to the limited information from the data and highly curved events at near offsets with strong lateral amplitude variations. Therefore, we have developed a novel workflow that first deploys a deep-learning-based reconstruction of the shallow reflections and then uses the reconstructed data as the input for the subsequent surface multiple removal. In particular, we use a convolutional neural network architecture - U-net that was developed from convolutional autoencoders with extra direct skip connections between different levels of encoders and the corresponding decoders. Instead of using field data directly in network training, the training set is carefully synthesized based on the prior water-layer information of the field data; thus, a fully sampled field data set, which is difficult to obtain, is not needed for training in our workflow. An inversion-based approach - closed-loop SRME - is used for the surface multiple removal, in which the primaries are directly estimated via full-waveform inversion in a data-driven manner. Finally, the effectiveness of our workflow is determined based on 2D North Sea field data in a shallow-water scenario (92.5 m water depth) with a relatively large minimum offset (150 m).@en

The journal is publishing this Expression of Concern to alert readers while we investigate to determine whether further action is required.@en

Strong noise is one of the toughest problems in the controlled-source electromagnetic (CSEM) method, which highly affects the quality of recorded data. The three main types of noise existing in CSEM data are periodic noise, Gaussian white noise, and nonperiodic noise, among which the nonperiodic noise is thought to be the most difficult to remove. We have developed a novel and effective method for removing such nonperiodic noise by formulating an inverse problem that is based on inverse discrete Fourier transform and several time windows in which only Gaussian white noise exists. These critical locations, which we call reconstruction locations, can be found by taking advantage of the continuous wavelet transform (CWT) and the temporal derivative of the scalogram generated by CWT. The coefficients of the nonperiodic noise are first estimated using the new least-squares method, and then they are subtracted from the coefficients of the raw data to produce denoised data. Together with the nonperiodic noise, we also remove Gaussian noise using the proposed method. We validate the methodology using real-world CSEM data.@en

As full waveform inversion (FWI) is a non-unique and typically ill-posed inversion problem, it needs proper regularization to avoid cycle-skipping. To reduce the non-linearity of FWI, Joint Migration Inversion (JMI) is proposed as an alternative, explaining the reflection data with decoupled velocity and reflectivity parameters. However, the velocity update may also suffer from being trapped in local minima. To optimally include geologic information, we propose FWI/JMI with directional total variation as an L1-norm regularization on the velocity. We design the directional total variation operator based on the local dip field, instead of ignoring the local structural direction of the subsurface and only using horizontal- and vertical- gradients in the traditional TV. The local dip field is estimated using plane-wave destruction based on a raw reflectivity model, which is usually calculated from the initial velocity model. With two complex synthetic examples, based on the Marmousi model, we demonstrate that the proposed method is much more effective compared to both FWI/JMI without regularization and FWI/JMI with the conventional TV regularization. In the JMI-based example, we also show that L1 directional TV works better than L2 directional Laplacian smoothing. In addition, by comparing these two examples, it can be seen that the impact of regularization is larger for FWI than for JMI, because in JMI the velocity model only explains the propagation effects and, thereby, makes it less sensitive to the details in the velocity model. @en

Five-dimensional (5D) seismic data reconstruction becomes more appealing in recent years because it takes advantage of five physical dimensions of the seismic data and can reconstruct data with large gap. The low-rank approximation approach is one of the most effective methods for reconstructing 5D dataset. However, the main disadvantage of the low-rank approximation method is its low computational efficiency because of many singular value decompositions (SVD) of the block Hankel/Toeplitz matrix in the frequency domain. In this paper, we develop an SVD-free low-rank approximation method for efficient and effective reconstruction and denoising of the seismic data that contain four spatial dimensions. Our SVD-free rank constraint model is based on an alternating minimization strategy, which updates one variable each time while fixing the other two. For each update, we only need to solve a linear least-squares problem with much less expensive QR factorization. The SVD-based and SVD-free low-rank approximation methods in the singular spectrum analysis (SSA) framework are compared in detail, regarding the reconstruction performance and computational cost. The comparison shows that the SVD-free low-rank approximation method can obtain similar reconstruction performance as the SVD-based method but with a large computational speedup.@en

Microseismic methods are crucial for real-Time monitoring of the hydraulic fracturing dynamic status during the development of unconventional reservoirs. However, unlike the active-source seismic events, the microseismic events usually have low signal-To-noise ratio (SNR), which makes its data processing challenging. To overcome the noise issue of the weak microseismic events, we propose a new workflow for high-resolution microseismic event detection. For the preprocessing, fix-sized segmentation with a length of 2∗wavelength is used to divide the data into segments. Later on, 191 features have been extracted and used as the input data to train the support vector machine (SVM) model. These features include 63 1-D time/spectral-domain features, and 128 2-D texture features, which indicate the continuity, smoothness, and irregularity of the events/noise. The proposed feature extraction maximally exploits the limited information of each segment. Afterward, we use a combination of univariate feature selection and random-forest-based recursive feature elimination for feature selection to avoid overfitting. This feature selection strategy not only finds the best features, but also decides the optimal number of features that are needed for the best accuracy. Regarding the training process, SVM with a Gaussian kernel is used. In addition, a cross-validation (CV) process is implemented for automatic parameter setting. In the end, a group of synthetic and field microseismic data with different levels of complexity show that the proposed workflow is much more robust than the state-of-The-Art short-Term-Average over long-Term-Average ratio (STA/LTA) method and also performs better than the convolutional-neural-networks (CNN), for this case where the amount of training data sets is limited. A demo for the synthetic example is available: https://github.com/shanqu91/ML_event_detection_microseismic.@en

Prestack seismic inversion has emerged as a powerful technique for reconstructing parameters attribute to the subsurface properties and building the geophysical parameter models. However, the inversion algorithms always suffer from spatial blur and low resolution. Total variation (TV) regularization preserves the spatial variation boundary of data by highlighting the sparsity of the first-order difference, which is regarded as an important technical means for image restoration. However, when the data do not change along the spatial grid direction, TV regularization is prone to a staircase effect. In this article, a directional TV (DTV) method is proposed to conduct the prestack amplitude variation with offset/angle (AVO/AVA) inversion. The method consists of three essential steps: estimating the seismic slope attribute from the seismic data, introducing seismic slope attribute to the TV regularization to establish the objective function, and optimizing the objective function by the split-Bregman algorithm. Finally, the conventional and proposed methods are applied to the synthetic and the real seismic data. The comparison of different methods demonstrates that the proposed method is applicable to reveal the detailed subsurface models, alleviate the staircase effect or artifact substantially, and further upgrade the quality of prestack inversion results.@en

Surface-related multiple elimination (SRME) has already been proven as a solid multiple and primary estimation tool for decades due to its data-driven property and strong physics behind. However, surface-related multiple leakage is still commonly seen in the SRME processed results, which might arise from the imperfect sampling and the balancing effect of the adaptive subtraction. Local primary-and-multiple orthogonalization (LPMO) is recently proposed to mitigate the multiple leakage. LPMO framework includes two separate steps: an initial multiple and primary estimation step via conservative SRME and an external multiple leakage extraction step via LPMO. Although decent leakage extraction can be achieved, LPMO requires a large computational cost due to many conjugate-gradient iterations within the shaping regularization based inversion framework. Assuming that the scalar LPMO weight is closely related to its neighboring time-and-space points, a scaled point-by-point division can be used to avoid the iterative inversion of LPMO. Therefore, we propose a fast LPMO (FLPMO) for surface-related multiple estimation. Applications on two different field data sets demonstrate the nearly same multiple leakage extraction performances for both LPMO and FLPMO, while showing that, the FLPMO is much faster than LPMO. @en

An important imaging challenge is creating reliable seismic images without internal multiple crosstalk, especially in cases with strong overburden reflectivity. Several data-driven methods have been proposed to attenuate the internal multiple crosstalk, for which fully sampled data in the source and receiver side are usually required. To overcome this acquisition constraint, model-driven full-wavefield migration (FWM) can automatically include internal multiples and only needs dense sampling in either the source or receiver side. In addition, FWM can correct for transmission effects at the reflecting interfaces. Although FWM has been shown to work effectively in compensating for transmission effects and suppressing internal multiple crosstalk compared to conventional least-squares primary wavefield migration (PWM), it tends to generate relatively weaker internal multiples during modeling. Therefore, some leaked internal multiple crosstalk can still be observed in the FWM image, which tends to blend in the background and can be misinterpreted as real geology. Thus, we adopted a novel framework using local primary-and-multiple orthogonalization (LPMO) on the FWM image as a postprocessing step for leaked internal multiple crosstalk estimation and attenuation. Due to their opposite correlation with the FWM image, a positive-only LPMO weight can be used to estimate the leaked internal multiple crosstalk, whereas a negative-only LPMO weight indicates the transmission effects that need to be retained. Application to North Sea field data validates the performance of the proposed framework for removing the weak but misleading leaked internal multiple crosstalk in the FWM image. Therefore, with this new framework, FWM can provide a reliable solution to the long-standing issue of imaging primaries and internal multiples automatically, with proper primary restoration.@en

The local signal-to-noise orthogonalization algorithm has been widely used in the community of seismic processing and imaging. It helps orthogonalize the signal-and-noise components in an elegant way so that the noise does not contain the signal leakage in seismic denoising. The traditional local signal-to-noise orthogonalization is based on solving a highly underdetermined, ill-posed inverse problem with local smoothness constraint. Due to the inversion nature, the local orthogonalization method requires a large number of iterations and thus is computationally demanding in large-scale applications. Here, we proposed a much accelerated signal-and-noise orthogonalization method, where we design an efficient way for calculating the orthogonalization weight. When new samples are involved in the calculation, we calculate the orthogonalization weight of the new samples by connecting them with the calculated weights of the previous samples. The orthogonalization weight needs to be smoothed and scaled after all samples have been processed to make the resulted orthogonalization weight smooth across the seismic data and match the amplitude level of the initially suppressed noise. In this way, we avoid iterations when calculating the orthogonalization weight. We apply the proposed method to several synthetic and field data examples, have a benchmark comparison with state-of-the-art algorithms, and demonstrate its much accelerated efficiency compared with the traditional local signal-and-noise orthogonalization.@en

Microseismic methods are crucial for real-time monitoring of the hydraulic fracturing dynamic status. However, unlike the active-source seismic events, the microseismic events usually have very low signal-to-noise ratio (SNR), which makes its processing challenging. To overcome the noise issue of the weak microseismic events, we propose a novel method for microseismic event detection based on the support vector machine (SVM) classification with a Gaussian kernel. For the preprocessing, fix-sized segmentation with a length of 2* wavelength is used to divide the data into segments. 123 features in total, which are used as input data to train the SVM model, have been extracted. These features include 63 1D time/spectral-domain features, and 60 2D texture features. Afterwards, we use a combination of univariate feature selection and random-forest-based recursive feature elimination for feature selection. This feature selection strategy not only finds the best features but also decides the number of features that are needed for the best accuracy. Regarding the essential training process, a C-SVM model, where coefficient C is used to control the tolerance of error item, is considered and a cross-validation (CV) process is implemented for automatic parameter setting. In the end, a group of synthetic and real microseismic data with different levels of complexity show the effectiveness of the proposed method. @en

We propose a new low-rank based noise attenuation method using an efficient algorithm for tracking subspaces from highly corrupted seismic observations. The subspace tracking algorithm requires only basic linear algebraic manipulations. The algorithm is derived by analysing incremental gradient descent on the Grassmannian manifold of subspaces. When the multidimensional seismic data are mapped to a low-rank space, the subspace tracking algorithm can be directly applied to the input low-rank matrix to estimate the useful signals. Since the subspace tracking algorithm is an online algorithm, it is more robust to random noise than traditional truncated singular value decomposition (TSVD) based subspace tracking algorithm. Compared with the state-of-The-Art algorithms, the proposed denoising method can obtain better performance. More specifically, the proposed method outperforms the TSVD-based singular spectrum analysis method in causing less residual noise and also in saving half of the computational cost. Several synthetic and field data examples with different levels of complexities demonstrate the effectiveness and robustness of the presented algorithm in rejecting different types of noise including random noise, spiky noise, blending noise, and coherent noise.@en

Amplitude-preserving data processing is an important and challenging topic in many scientific fields. The amplitude-variation details in seismic data are especially important because the amplitude variation is directly related with the subsurface wave impedance and fluid characteristics. We propose a novel seismic noise attenuation approach that is based on local plane-wave assumption of seismic events and the amplitude preserving capability of the orthogonal polynomial transform (OPT). The OPT is a way for representing spatially correlative seismic data as a superposition of polynomial basis functions, by which the random noise is distinguished from the useful energy by the high orthogonal polynomial coefficients. The seismic energy is the most correlative along the structural direction and thus the OPT is optimally performed in a flattened gather. We introduce in detail the flattening operator for creating the flattened dimension, where the OPT can be applied subsequently. The flattening operator is created by deriving a plane-wave trace continuation relation following the plane-wave equation. We demonstrate that both plane-wave trace continuation and OPT can well preserve the strong amplitude variation existing in seismic data. In order to obtain a robust slope estimation performance in the presence of noise, a robust slope estimation approach is introduced to substitute the traditional method. A group of synthetic, pre-stack and post-stack field seismic data are used to demonstrate the potential of the proposed framework in realistic applications.@en

Accurate multiple removal remains an important step in seismic data processing sequences. Most multiple removal methods, such as surface-related multiple elimination (SRME), consist of a multiple prediction step and an adaptive subtraction step. Due to imperfect circumstances (e.g., coarse data sampling) or built-in assumptions (e.g., 2D method versus 3D data), multiple leakage is commonly observed in the results. More aggressive adaptive multiple subtraction can reduce the leakage problem, for example, by using small local windows and a long filter length, but at the risk of severely damaging the primaries due to overfitting. In contrast, conservative adaptive subtraction with large or global windows and a short filter length can preserve most primary energy while tending to have more multiple leakage because of underfitting. Assuming that the primaries and multiples do not correlate locally in the time-space domain, our solution to this problem is to extract the leaked multiples from the initially estimated primaries using local primary-and-multiple orthogonalization (LPMO) rather than restoring the damaged primaries. Our framework consists of two steps: an initial primary estimation step and a multiple leakage extraction step. The initial step corresponds to conservative SRME (or an equivalent method) that produces the initially estimated primary and multiple models. The second step is based on LPMO to retrieve the leaked multiples from the estimated primaries via a time- and space-varying weight function that is estimated from the local correlation of predicted multiples and residual multiples in the estimated primaries with the help of shaping regularization. In this way, we can obtain a better primary model that has much less leaked multiple energy and less primary damage at the same time. We find good performance of our framework via two synthetic data examples and one field data example.@en

Microseismic methods are crucial for real-timemonitoring of the hydraulic fracturing dynamic status during the development of unconventional reservoirs. However, unlike the active-source seismic events, the microseismic events usually have low signal-to-noise ratio (SNR), which makes its data processing challenging. To overcome the noise issue of the weak microseismic events, we propose a new workflow for high-resolution microseismic event detection. For the preprocessing, fix-sized segmentation with a length of 2.wavelength is used to divide the data into segments. Later on, 191 features have been extracted and used as the input data to train the support vector machine (SVM) model. These features include 63 1-D time/spectraldomain features, and 128 2-D texture features, which indicate the continuity, smoothness, and irregularity of the events/noise. The proposed feature extractionmaximally exploits the limited information of each segment. Afterward, we use a combination of univariate feature selection and random-forest-based recursive feature elimination for feature selection to avoid overfitting. This feature selection strategy not only finds the best features, but also decides the optimal number of features that are needed for the best accuracy. Regarding the training process, SVM with a Gaussian kernel is used. In addition, a cross-validation (CV) process is implemented for automatic parameter setting. In the end, a group of synthetic and field microseismic data with different levels of complexity show that the proposed workflow is much more robust than the state-of-the-art short-term-average over long-term-average ratio (STA/LTA) method and also performs better than the convolutional-neural-networks (CNN), for this case where the amount of training data sets is limited. A demo for the synthetic example is available: https://github.com/shanqu91/ML event detection microseismic.@en

Simultaneous source technology can accelerate data acquisition and improve subsurface illumination. But those advantages are compromised due to dense interference. To address the intense interference in simultaneous source data, we have investigated a method based on a deep neural network. The designed architecture consists of convolutional and deconvolutional networks. The convolutional network can learn the local features of the training data set, and the deconvolutional network constructs the output using the extracted features to match the ground truth. Because the main computational cost results from the optimization of the network parameters, the trained network can separate simultaneous source data efficiently. Besides, with the given dithering code, we embed the trained network into an iterative framework that can further improve the deblending. A numerical test on synthetic data demonstrates that the iterative framework with the trained network can obtain comparable performance with high efficiency compared to the conventional method. Next, we test our method with two different trained networks (one is from a synthetic data set, and the other is from a field data set) on field data. The test results confirm the performance of our method. @en

Amplitude-preserving data processing is an important and challenging topic in many scientific fields. The amplitude-variation details in seismic data are especially important because the amplitude variation is directly related with the subsurface wave impedance and fluid characteristics. We propose a novel seismic noise attenuation approach that is based on local plane-wave assumption of seismic events and the amplitude preserving capability of the orthogonal polynomial transform (OPT). The OPT is a way for representing spatially correlative seismic data as a superposition of polynomial basis functions, by which the random noise is distinguished from the useful energy by the high orthogonal polynomial coefficients. The seismic energy is the most correlative along the structural direction and thus the OPT is optimally performed in a flattened gather. We introduce in detail the flattening operator for creating the flattened dimension, where the OPT can be applied subsequently. The flattening operator is created by deriving a plane-wave trace continuation relation following the plane-wave equation. We demonstrate that both plane-wave trace continuation and OPT can well preserve the strong amplitude variation existing in seismic data. In order to obtain a robust slope estimation performance in the presence of noise, a robust slope estimation approach is introduced to substitute the traditional method. A group of synthetic, pre-stack and post-stack field seismic data are used to demonstrate the potential of the proposed framework in realistic applications.@en

We propose a new low-rank based noise attenuation method using an efficient algorithm for tracking subspaces from highly corrupted seismic observations. The subspace tracking algorithm requires only basic linear algebraic manipulations. The algorithm is derived by analysing incremental gradient descent on the Grassmannian manifold of subspaces. When the multidimensional seismic data are mapped to a low-rank space, the subspace tracking algorithm can be directly applied to the input low-rank matrix to estimate the useful signals. Since the subspace tracking algorithm is an online algorithm, it is more robust to random noise than traditional truncated singular value decomposition (TSVD) based subspace tracking algorithm. Compared with the state-of-the-art algorithms, the proposed denoising method can obtain better performance. More specifically, the proposed method outperforms the TSVDbased singular spectrum analysis method in causing less residual noise and also in saving half of the computational cost. Several synthetic and field data examples with different levels of complexities demonstrate the effectiveness and robustness of the presented algorithm in rejecting different types of noise including random noise, spiky noise, blending noise, and coherent noise.@en

Surface-related multiple elimination (SRME) is a solid and effective approach for primary estimation. However, due to the imperfections in data and method multiple energy leakage is commonly seen in the results of SRME-predicted primaries. Assuming that the primaries and multiples do not correlate locally in the time-space domain, we are able to extract the leaked multiples from the initially estimated primaries using multi-domain local primary-and-multiple orthogonalization. The proposed framework consists of two steps: an initial primary/multiple estimation step and a multiple-leakage extraction step. The initial step corresponds to SRME, which produces the initial estimated primary and multiple models. The second step is based on multi-domain local primary-and-multiple orthogonalization to retrieve the leaked multiples. Multi-domain indicates that we first extract the leaked multiples in shot domain, and then the residual can be further extracted in common-offset domain. Thus, we can obtain a better primary output which has much less leaked multiple energy. We demonstrate a good performance of our proposed framework on both synthetic and field data, where it repairs the leakage of standard global adaptive subtraction.@en