Recognition and reconstruction of coherent energy with application to deep seismic reflection data

Reflections in deep seismic reflection data tend to be visible on only a limited number of traces in a common midpoint gather. To prevent stack degeneration, any noncoherent reflection energy has to be removed. In this paper, a standard classification technique in remote sensing is presented to enhance data quality. It consists of a recognition technique to detect and extract coherent energy in both common shot gathers and fi- nal stacks. This technique uses the statistics of a picked seismic phase to obtain the likelihood distribution of its presence. Multiplication of this likelihood distribution with the original data results in a “cleaned up” section. Application of the technique to data from a deep seismic reflection experiment enhanced the visibility of all reflectors considerably. Because the recognition technique cannot produce an estimate of “missing” data, it is extended with a reconstruction method. Two methods are proposed: application of semblance weighted local slant stacks after recognition, and direct recognition in the linear tau-p domain. In both cases, the power of the stacking process to increase the signal-to-noise ratio is combined with the direct selection of only specific seismic phases. The joint application of recognition and reconstruction resulted in data images which showed reflectors more clearly than application of a single technique

Signal processing techniques for the enhancement of marine seismic data

This thesis presents several signal processing techniques applied to the enhancement of marine seismic data. Marine seismic exploration provides an image of the Earth's subsurface from reflected seismic waves. Because the recorded signals are contaminated by various sources of noise, minimizing their effects with new attenuation techniques is necessary. A statistical analysis of background noise is conducted using Thomson’s multitaper spectral estimator and Parzen's amplitude density estimator. The results provide a statistical characterization of the noise which we use for the derivation of signal enhancement algorithms. Firstly, we focus on single-azimuth stacking methodologies and propose novel stacking schemes using either enhanced weighted sums or a Kalman filter. It is demonstrated that the enhanced methods yield superior results by their ability to exhibit cleaner and better defined reflected events as well as a larger number of reflections in deep waters. A comparison of the proposed stacking methods with existing ones is also discussed. We then address the problem of random noise attenuation and present an innovative application of sparse code shrinkage and independent component analysis. Sparse code shrinkage is a valuable method when a noise-free realization of the data is generated to provide data-driven shrinkages. Several models of distribution are investigated, but the normal inverse Gaussian density yields the best results. Other acceptable choices of density are discussed as well. Finally, we consider the attenuation of flow-generated nonstationary coherent noise and seismic interference noise. We suggest a multiple-input adaptive noise canceller that utilizes a normalized least mean squares alg orithm with a variable normalized step size derived as a function of instantaneous frequency. This filter attenuates the coherent noise successfully when used either by itself or in combination with a time-frequency median filter, depending on the noise spectrum and repartition along the data. Its application to seismic interference attenuation is also discussed

Unrest at Domuyo Volcano, Argentina, detected by geophysical and geodetic data and morphometric analysis

New volcanic unrest has been detected in the Domuyo Volcanic Center (DVC), to the east of the Andes Southern Volcanic Zone in Argentina. To better understand this activity, we investigated new seismic monitoring data, gravimetric and magnetic campaign data, and interferometric synthetic aperture radar (InSAR) deformation maps, and we derived an image of the magma plumbing system and the likely source of the unrest episode. Seismic events recorded during 2017-2018 nucleate beneath the southwestern flank of the DVC. Ground deformation maps derived from InSAR processing of Sentinel-1 data exhibit an inflation area exceeding 300 km2, from 2014 to at least March 2018, which can be explained by an inflating sill model located 7 km deep. The Bouguer anomaly reveals a negative density contrast of ~35 km wavelength, which is spatially coincident with the InSAR pattern. Our 3D density modeling suggests a body approximately 4-6 km deep with a density contrast of -550 kg/m3. Therefore, the geophysical and geodetic data allow identification of the plumbing system that is subject to inflation at these shallow crustal depths. We compared the presence and dimensions of the inferred doming area to the drainage patterns of the area, which support long-established incremental uplift according to morphometric analysis. Future studies will allow us to investigate further whether the new unrest is hydrothermal or magmatic in origin.Fil: Astort, Ana. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Estudios Andinos "Don Pablo Groeber". Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Estudios Andinos "Don Pablo Groeber"; ArgentinaFil: Walter, Thomas R. German Research Centre for Geosciences; AlemaniaFil: Ruiz, Francisco. Universidad Nacional de San Juan. Facultad de Ciencias Exactas, Físicas y Naturales. Instituto Geofísico Sismológico Volponi; ArgentinaFil: Sagripanti, Lucía. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Estudios Andinos "Don Pablo Groeber". Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Estudios Andinos "Don Pablo Groeber"; ArgentinaFil: Nacif, Andres Antonio. Universidad Nacional de San Juan. Facultad de Ciencias Exactas, Físicas y Naturales. Instituto Geofísico Sismológico Volponi; ArgentinaFil: Acosta, Gemma. Universidad Nacional de San Juan. Facultad de Ciencias Exactas, Físicas y Naturales. Instituto Geofísico Sismológico Volponi; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Folguera Telichevsky, Andres. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Estudios Andinos "Don Pablo Groeber". Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Estudios Andinos "Don Pablo Groeber"; Argentin

APPLICATION OF SEISMIC RADIAL ANISOTROPY FOR NEAR-SURFACE FRACTURES IDENTIFICATION

Fractures significantly control the groundwater flow and solute transport in geological settings of low-permeable rocks. Fractures also affect seismic wave propagation. For instance, they can create a directional dependence of seismic velocity with respect to their orientations, known as seismic anisotropy. Seismic radial anisotropy as used here is the difference between the velocity of a vertically polarized S-wave (SV) and one polarized horizontally (SH). In this thesis, seismic radial anisotropy was used to evaluate its usefulness for correlating with near-surface fractures. The seismic radial anisotropy models were obtained at two sites from dispersion analyses of the Rayleigh waves, with vertical polarization, and Love waves, with horizontal polarization, using the Multichannel Analysis of Surface Waves (MASW) method. The seismic radial anisotropies at these two sites in different geological settings (one metamorphic-igneous bedrock and the other sedimentary), shows a strong correlation of seismic radial anisotropy with near surface fractures, and hence, can be used to characterize near-surface fractures

Hypothesis-based machine learning for deep-water channel systems

2020 Spring.Includes bibliographical references.Machine learning algorithms are readily being incorporated into petroleum industry workflows for use in well-log correlation, prediction of rock properties, and seismic data interpretation. However, there is a clear disconnect between sedimentology and data analytics in these workflows because sedimentologic data is largely qualitative and descriptive. Sedimentology defines stratigraphic architecture and heterogeneity, which can greatly impact reservoir quality and connectivity and thus hydrocarbon recovery. Deep-water channel systems are an example where predicting reservoir architecture is critical to mitigating risk in hydrocarbon exploration. Deep-water reservoirs are characterized by spatial and temporal variations in channel body stacking patterns, which are difficult to predict with the paucity of borehole data and low quality seismic available in these remote locations. These stacking patterns have been shown to be a key variable that controls reservoir connectivity. In this study, the gap between sedimentology and data analytics is bridged using machine learning algorithms to predict stratigraphic architecture and heterogeneity in a deep-water slope channel system. The algorithms classify variables that capture channel stacking patterns (i.e., channel positions: axis, off-axis, and margin) from a database of outcrop statistics sourced from 68 stratigraphic measured sections from outcrops of the Upper Cretaceous Tres Pasos Formation at Laguna Figueroa in the Magallanes Basin, Chile. An initial hypothesis that channel position could be predicted from 1D descriptive sedimentologic data was tested with a series of machine learning algorithms and classification schemes. The results confirmed this hypothesis as complex algorithms (i.e., random forest, XGBoost, and neural networks) achieved accuracies above 80% while less complex algorithms (i.e., decision trees) achieved lower accuracies between 60%-70%. However, certain classes were difficult for the machine learning algorithms to classify, such as the transitional off-axis class. Additionally, an interpretive classification scheme performed better (by around 10%-20% in some cases) than a geometric scheme that was devised to remove interpretation bias. However, outcrop observations reveal that the interpretive classification scheme may be an over-simplified approach and that more heterogeneity likely exists in each class as revealed by the geometric scheme. A refined hypothesis was developed that a hierarchical machine learning approach could lend deeper insight into the heterogeneity within sedimentologic classes that are difficult for an interpreter to discern by observation alone. This hierarchical analysis revealed distinct sub-classes in the margin channel position that highlight variations in margin depositional style. The conceptual impact of these varying margin styles on fluid flow and connectivity is shown

Direct And Evolutionary Approaches For Optimal Receiver Function Inversion

Receiver functions are time series obtained by deconvolving vertical component seismograms from radial component seismograms. Receiver functions represent the impulse response of the earth structure beneath a seismic station. Generally, receiver functions consist of a number of seismic phases related to discontinuities in the crust and upper mantle. The relative arrival times of these phases are correlated with the locations of discontinuities as well as the media of seismic wave propagation. The Moho (Mohorovicic discontinuity) is a major interface or discontinuity that separates the crust and the mantle. In this research, automatic techniques to determine the depth of the Moho from the earth’s surface (the crustal thickness H) and the ratio of crustal seismic P-wave velocity (Vp) to S-wave velocity (Vs) (ï«= Vp/Vs) were developed. In this dissertation, an optimization problem of inverting receiver functions has been developed to determine crustal parameters and the three associated weights using evolutionary and direct optimization techniques

Seismic Ray Impedance Inversion

This thesis investigates a prestack seismic inversion scheme implemented in the ray parameter domain. Conventionally, most prestack seismic inversion methods are performed in the incidence angle domain. However, inversion using the concept of ray impedance, as it honours ray path variation following the elastic parameter variation according to Snell’s law, shows the capacity to discriminate different lithologies if compared to conventional elastic impedance inversion. The procedure starts with data transformation into the ray-parameter domain and then implements the ray impedance inversion along constant ray-parameter profiles. With different constant-ray-parameter profiles, mixed-phase wavelets are initially estimated based on the high-order statistics of the data and further refined after a proper well-to-seismic tie. With the estimated wavelets ready, a Cauchy inversion method is used to invert for seismic reflectivity sequences, aiming at recovering seismic reflectivity sequences for blocky impedance inversion. The impedance inversion from reflectivity sequences adopts a standard generalised linear inversion scheme, whose results are utilised to identify rock properties and facilitate quantitative interpretation. It has also been demonstrated that we can further invert elastic parameters from ray impedance values, without eliminating an extra density term or introducing a Gardner’s relation to absorb this term. Ray impedance inversion is extended to P-S converted waves by introducing the definition of converted-wave ray impedance. This quantity shows some advantages in connecting prestack converted wave data with well logs, if compared with the shearwave elastic impedance derived from the Aki and Richards approximation to the Zoeppritz equations. An analysis of P-P and P-S wave data under the framework of ray impedance is conducted through a real multicomponent dataset, which can reduce the uncertainty in lithology identification.Inversion is the key method in generating those examples throughout the entire thesis as we believe it can render robust solutions to geophysical problems. Apart from the reflectivity sequence, ray impedance and elastic parameter inversion mentioned above, inversion methods are also adopted in transforming the prestack data from the offset domain to the ray-parameter domain, mixed-phase wavelet estimation, as well as the registration of P-P and P-S waves for the joint analysis. The ray impedance inversion methods are successfully applied to different types of datasets. In each individual step to achieving the ray impedance inversion, advantages, disadvantages as well as limitations of the algorithms adopted are detailed. As a conclusion, the ray impedance related analyses demonstrated in this thesis are highly competent compared with the classical elastic impedance methods and the author would like to recommend it for a wider application

Slab Tearing Underneath the Bransfield Strait, Antarctica

We conduct a P-wave receiver function analysis of the Bransfield Strait (West Antarctica) to determine the lithospheric structure of this back-arc basin, thanks to 31 temporary and permanent stations. Our main finding is a 15 km tear of the Phoenix slab, coinciding with the location of the 2020–2021 Orca earthquake swarm's epicenters. Teleseismic wave modeling reveals that the two major earthquakes occurred at the base of the crust, suggesting that the swarm could have been triggered by active underplating driven by mantle flow through the slab tear. There is evidence for such an underplating layer at least under Deception Island and for a widespread low velocity zone in the mantle wedge probably undergoing partial melting. We found average crustal thickness (30.5 ± 1.0 km) and Vp/Vs (1.81 ± 0.04) values close to average extended continental crust, although results in the South Shetland Islands are significantly more heterogeneous than in the Antarctic Peninsula.Spanish national projects PID2019-109608GB-100/ SRA/10.13039/501100011033CMT2016-77315-R, the Andalusian regional project A-RNM-421-UGR18FPI Grant PRE2020-092556 (funded by MCIN/AEI/10.13039/501100011033 and the European Social Fund