Simultaneous reconstruction of seismic reflections and diffractions using a global hyperbolic Radon dictionary

Geophysics ◽  
2018 ◽  
Vol 83 (6) ◽  
pp. V315-V323 ◽  
Author(s):  
Amr Ibrahim ◽  
Paolo Terenghi ◽  
Mauricio D. Sacchi

We have developed a new transform with basis functions that closely resemble seismic reflections and diffractions. The new transform is an extension of the classic hyperbolic Radon transform and accounts for the apex shifts of the seismic reflection hyperbolas and the asymptote shifts of the seismic diffraction hyperbolas. The adjoint and forward operators of the proposed transform are computed using Stolt operators in the frequency domain to increase the computational efficiency of the transform. This new transform is used, in conjunction with a sparse inversion algorithm, to reconstruct common-shot gathers. Our tests indicate that this new transform is an efficient tool for interpolating coarsely sampled seismic data in cases in which one cannot use small data windows to validate the linear event assumption that is often made by Fourier-based reconstruction methods.

Geophysics ◽  
1988 ◽  
Vol 53 (7) ◽  
pp. 894-902 ◽  
Author(s):  
Ruhi Saatçilar ◽  
Nezihi Canitez

Amplitude‐ and frequency‐modulated wave motion constitute the ground‐roll noise in seismic reflection prospecting. Hence, it is possible to eliminate ground roll by applying one‐dimensional, linear frequency‐modulated matched filters. These filters effectively attenuate the ground‐roll energy without damaging the signal wavelet inside or outside the ground roll’s frequency interval. When the frequency bands of seismic reflections and ground roll overlap, the new filters eliminate the ground roll more effectively than conventional frequency and multichannel filters without affecting the vertical resolution of the seismic data.


Geophysics ◽  
2015 ◽  
Vol 80 (1) ◽  
pp. R31-R41 ◽  
Author(s):  
Andrea Zunino ◽  
Klaus Mosegaard ◽  
Katrine Lange ◽  
Yulia Melnikova ◽  
Thomas Mejer Hansen

Determination of a petroleum reservoir structure and rock bulk properties relies extensively on inference from reflection seismology. However, classic deterministic methods to invert seismic data for reservoir properties suffer from some limitations, among which are the difficulty of handling complex, possibly nonlinear forward models, and the lack of robust uncertainty estimations. To overcome these limitations, we studied a methodology to invert seismic reflection data in the framework of the probabilistic approach to inverse problems, using a Markov chain Monte Carlo (McMC) algorithm with the goal to directly infer the rock facies and porosity of a target reservoir zone. We thus combined a rock-physics model with seismic data in a single inversion algorithm. For large data sets, the McMC method may become computationally impractical, so we relied on multiple-point-based a priori information to quantify geologically plausible models. We tested this methodology on a synthetic reservoir model. The solution of the inverse problem was then represented by a collection of facies and porosity reservoir models, which were samples of the posterior distribution. The final product included probability maps of the reservoir properties in obtained by performing statistical analysis on the collection of solutions.


Geophysics ◽  
2010 ◽  
Vol 75 (6) ◽  
pp. WB103-WB111 ◽  
Author(s):  
Side Jin

Regularizing inadequate and irregularly sampled seismic data is one of important problems in seismic data processing. An improvement to existing methods to solve this problem is proposed by applying a 5D regularization/interpolation scheme with a damped least-norm Fourier inversion. Under the assumption of planar seismic events within small data windows, the spatial spectrum of regularized data for a fixed frequency should be sparse and have minimum damped norm. The inversion scheme consists in finding a set of regularly spaced spatial Fourier coefficients by minimizing its damped norm for each frequency, subject to the condition that the resulting spatial Fourier coefficients also faithfully reconstruct the original data. The damping factors are automatically derived from the amplitude spectra of the regularized low-frequency data. With the guidance of the damping factors and automatic adjustment of wavenumber ranges according to the Nyquist sampling theory, the proposed inversion algorithm naturally yields a one-step solution for both stabilization and antialiasing of the interpolation problem. A distinctive feature of the method is that it uses high-dimensional nonuniform fast Fourier transforms to evaluate expensive discrete Fourier transforms involved in conjugate gradient iterations. This improves the computational efficiency. The results of applying this algorithm to synthetic and field data demonstrate that it performs well when applied to highly irregular data and outperforms lower dimensional interpolation schemes.


1992 ◽  
Vol 129 (5) ◽  
pp. 633-636 ◽  
Author(s):  
N. R. Goulty ◽  
M. Leggett ◽  
T. Douglas ◽  
C. H. Emeleus

AbstractWe have conducted a seismic reflection test over a short profile on the granite of the Skye Tertiary central intrusive complex. From previous gravity modelling work it had been inferred that the granite is approximately 1.5 km thick and overlies basic rocks. The seismic data indicate that the granite is at least 2 km thick at the test location. Reflection events of alternating polarity between depths of 2.1 and 2.4 km suggest that basic and acidic sheets are interlayered at the base of the granitic mass.


Geophysics ◽  
2001 ◽  
Vol 66 (6) ◽  
pp. 1838-1842 ◽  
Author(s):  
C. M. Schmeissner ◽  
K. T. Spikes ◽  
D. W. Steeples

Ultrashallow seismic reflection surveys require dense spatial sampling during data acquisition, which increases their cost. In previous efforts to find ways to reduce these costs, we connected geophones rigidly to pieces of channel iron attached to a farm implement. This method allowed us to plant the geophones in the ground quickly and automatically. The rigidly interconnected geophones used in these earlier studies detected first‐arrival energy along with minor interfering seismic modes, but they did not detect seismic reflections. To examine further the feasibility of developing rigid geophone emplacement systems to detect seismic reflections, we experimented with four pieces of channel iron, each 2.7 m long and 10 cm wide. Each segment was equipped with 18 geophones rigidly attached to the channel iron at 15‐cm intervals, and the spikes attached to all 18 geophones were pushed into the ground simultaneously. The geophones detected both refracted and reflected energy; however, no significant signal distortion or interference attributable to the rigid coupling of the geophones to the channel iron was observed in the data. The interfering seismic modes mentioned from the previous experiments were not detected, nor was any P‐wave propagation noted within the channel iron. These results show promise for automating and reducing the cost of ultrashallow seismic reflection and refraction surveys.


2021 ◽  
Author(s):  
Pimpawee Sittipan ◽  
Pisanu Wongpornchai

Some of the important petroleum reservoirs accumulate beneath the seas and oceans. Marine seismic reflection method is the most efficient method and is widely used in the petroleum industry to map and interpret the potential of petroleum reservoirs. Multiple reflections are a particular problem in marine seismic reflection investigation, as they often obscure the target reflectors in seismic profiles. Multiple reflections can be categorized by considering the shallowest interface on which the bounces take place into two types: internal multiples and surface-related multiples. Besides, the multiples can be categorized on the interfaces where the bounces take place, a difference between long-period and short-period multiples can be considered. The long-period surface-related multiples on 2D marine seismic data of the East Coast of the United States-Southern Atlantic Margin were focused on this research. The seismic profile demonstrates the effectiveness of the results from predictive deconvolution and the combination of surface-related multiple eliminations (SRME) and parabolic Radon filtering. First, predictive deconvolution applied on conventional processing is the method of multiple suppression. The other, SRME is a model-based and data-driven surface-related multiple elimination method which does not need any assumptions. And the last, parabolic Radon filtering is a moveout-based method for residual multiple reflections based on velocity discrimination between primary and multiple reflections, thus velocity model and normal-moveout correction are required for this method. The predictive deconvolution is ineffective for long-period surface-related multiple removals. However, the combination of SRME and parabolic Radon filtering can attenuate almost long-period surface-related multiple reflections and provide a high-quality seismic images of marine seismic data.


2018 ◽  
Vol 123 (12) ◽  
pp. 10,810-10,830
Author(s):  
Michael Dentith ◽  
Huaiyu Yuan ◽  
Ruth Elaine Murdie ◽  
Perla Pina-Varas ◽  
Simon P. Johnson ◽  
...  

Geophysics ◽  
1996 ◽  
Vol 61 (5) ◽  
pp. 1351-1362 ◽  
Author(s):  
B. A. Hardage ◽  
D. L. Carr ◽  
D. E. Lancaster ◽  
J. L. Simmons ◽  
D. S. Hamilton ◽  
...  

A multidisciplinary team, composed of stratigraphers, petrophysicists, reservoir engineers, and geophysicists, studied a portion of Boonsville gas field in the Fort Worth Basin of North‐Central Texas to determine how modern geophysical, geological, and engineering techniques could be combined to understand the mechanisms by which fluvio‐deltaic depositional processes create reservoir compartmentalization in a low‐ to moderate‐accommodation basin. An extensive database involving well logs, cores, production, and pressure data from 200‐plus wells, [Formula: see text] [Formula: see text] of 3-D seismic data, vertical seismic profiles (VSPs), and checkshots was assembled to support this investigation. The reservoir system we studied was the Bend Conglomerate, a productive series of gas reservoirs composed of Middle Pennsylvanian fluvio‐deltaic clastics 900 to 1300 ft (275 to 400 m) thick in our project area. We were particularly interested in this reservoir system because evidence suggested that many of the sequences in this stratigraphic interval were deposited in low‐accommodation conditions (that is, in an environment where there was limited vertical space available for sediment accumulation), and our objective was to investigate how fluvio‐deltaic reservoirs were compartmentalized by low‐accommodation depositional processes. Using an extensive well log database (200 plus wells) and a core‐calibrated calculation of rock facies derived from these logs, we divided the Bend Conglomerate interval into ten genetic sequences, with each sequence being approximately 100 ft (30 m) thick. We then used local VSP and checkshot control to transform log‐measured depths of each sequence boundary to seismic two‐way time coordinates and identified narrow seismic data windows encompassing each sequence across the [Formula: see text] [Formula: see text] 3-D seismic grid. A series of seismic attributes was calculated in these carefully defined data windows to determine which attributes were reliable indicators of the presence of productive reservoir facies and which attributes could, therefore, reveal distinct reservoir compartments and potentially show where infield wells should be drilled to reach previously uncontacted gas reservoirs. Our best success was the seismic attribute correlations we found in the Upper and Lower Caddo sequences, at the top of the Bend Conglomerate. These sequences were deposited in a low‐accommodation setting, relative to other Boonsville sequences, and we found that reflection amplitude and instantaneous frequency, respectively, were reliable indicators of the areal distribution of reservoir facies in these low‐accommodation sequences.


2021 ◽  
Author(s):  
Piotr Krzywiec ◽  
Łukasz Słonka ◽  
Quang Nguyen ◽  
Michał Malinowski ◽  
Mateusz Kufrasa ◽  
...  

<p>In 2016, approximately 850 km of high-resolution multichannel seismic reflection data of the BALTEC survey have been acquired offshore Poland within the transition zone between the East European Craton and the Paleozoic Platform. Data processing, focused on removal of multiples, strongly overprinting geological information at shallower intervals, included SRME, TAU-P domain deconvolution, high resolution parabolic Radon demultiple and SWDM (Shallow Water De-Multiple). Entire dataset was Kirchhoff pre-stack time migrated. Additionally, legacy shallow high-resolution multichannel seismic reflection data acquired in this zone in 1997 was also used. All this data provided new information on various aspects of the Phanerozoic evolution of this area, including Late Cretaceous to Cenozoic tectonics and sedimentation. This phase of geological evolution could be until now hardly resolved by analysis of industry seismic data as, due to limited shallow seismic imaging and very strong overprint of multiples, essentially no information could have been retrieved from this data for first 200-300 m. Western part of the BALTEC dataset is located above the offshore segment of the Mid-Polish Swell (MPS) – large anticlinorium formed due to inversion of the axial part of the Polish Basin. BALTEC seismic data proved that Late Cretaceous inversion of the Koszalin – Chojnice fault zone located along the NE border of the MPS was thick-skinned in nature and was associated with substantial syn-inversion sedimentation. Subtle thickness variations and progressive unconformities imaged by BALTEC seismic data within the Upper Cretaceous succession in vicinity of the Kamień-Adler and the Trzebiatów fault zones located within the MPS documented complex interplay of Late Cretaceous basin inversion, erosion and re-deposition. Precambrian basement of the Eastern, cratonic part of the study area is overlain by Cambro-Silurian sedimentary cover. It is dissected by a system of steep, mostly reverse faults rooted in most cases in the deep basement. This fault system has been regarded so far as having been formed mostly in Paleozoic times, due to the Caledonian orogeny. As a consequence, Upper Cretaceous succession, locally present in this area, has been vaguely defined as a post-tectonic cover, locally onlapping uplifted Paleozoic blocks. New seismic data, because of its reliable imaging of the shallowest substratum, confirmed that at least some of these deeply-rooted faults were active as a reverse faults in latest Cretaceous – earliest Paleogene. Consequently, it can be unequivocally proved that large offshore blocks of Silurian and older rocks presently located directly beneath the Cenozoic veneer must have been at least partly covered by the Upper Cretaceous succession; then, they were uplifted during the widespread inversion that affected most of Europe. Ensuing regional erosion might have at least partly provided sediments that formed Upper Cretaceous progradational wedges recently imaged within the onshore Baltic Basin by high-end PolandSPAN regional seismic data. New seismic data imaged also Paleogene and younger post-inversion cover. All these results prove that Late Cretaceous tectonics substantially affected large areas located much farther towards the East than previously assumed.</p><p>This study was funded by the Polish National Science Centre (NCN) grant no UMO-2017/27/B/ST10/02316.</p>


Sign in / Sign up

Export Citation Format

Share Document