Reverse-time demigration using the extended-imaging condition

Geophysics ◽  
2014 ◽  
Vol 79 (3) ◽  
pp. WA97-WA105 ◽  
Author(s):  
Wiktor Waldemar Weibull ◽  
Børge Arntsen

The forward and inverse process of seismic migration and demigration or remodeling has many useful applications in seismic data processing. We evaluated a method to reobtain the seismic reflection data after migration, by inverting the common image point gathers produced by reverse-time migration (RTM) with an extended-imaging condition. This provided a transformation of the results of seismic data processing in the image domain back to the data domain. To be able to reconstruct the data with high fidelity, we set up demigration as a least-squares inverse problem and we solved it iteratively using a steepest-descent method. Because we used an extended-imaging condition, the method is not dependent on an accurate estimate of the migration-velocity field, and it is able to accurately reconstruct both primaries and multiples. At the same time, because the method is based on RTM, it can accurately handle seismic reflection data acquired over complex geologic media. Numerical results showed the feasibility of the method and highlighted some of its applications on 2D synthetic and field data sets.

Geophysics ◽  
1967 ◽  
Vol 32 (2) ◽  
pp. 207-224 ◽  
Author(s):  
John D. Marr ◽  
Edward F. Zagst

The more recent developments in common‐depth‐point techniques to attenuate multiple reflections have resulted in an exploration capability comparable to the development of the seismic reflection method. The combination of new concepts in digital seismic data processing with CDP techniques is creating unforeseen exploration horizons with vastly improved seismic data. Major improvements in multiple reflection and reverberation attenuation are now attainable with appropriate CDP geometry and special CDP stacking procedures. Further major improvements are clearly evident in the very near future with the use of multichannel digital filtering‐stacking techniques and the application of deconvolution as the first step in seismic data processing. CDP techniques are briefly reviewed and evaluated with real and experimental data. Synthetic data are used to illustrate that all seismic reflection data should be deconvolved as the first processing step.


Geophysics ◽  
2019 ◽  
Vol 84 (2) ◽  
pp. KS29-KS38
Author(s):  
Carlos Almagro Vidal ◽  
Joost van der Neut ◽  
Kees Wapenaar

Time-lapse changes in the subsurface can be analyzed by comparing seismic reflection data from two different states, one serving as the base survey and the second as the monitor survey. Conventionally, reflection data are acquired by placing active seismic sources at the acquisition surface. Alternatively, these data can be acquired from passive sources in the subsurface, using seismic interferometry (SI). Unfortunately, the reflection responses as retrieved by SI inherit an imprint of the passive-source distribution; therefore, monitoring with SI requires high passive-source repeatability, which is very often not achievable in practice. We have developed an alternative by using active seismic data for the base survey and a single passive source (e.g., a seismic tremor produced by induced seismicity) for the monitor survey. By constraining the source-radiation pattern of the (active) base survey according to the characteristics of the (passive) monitor survey, we succeed in extracting the time-lapse response in the image domain. We tested our method with numerically modeled 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 ◽  
...  

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>


Geophysics ◽  
1995 ◽  
Vol 60 (2) ◽  
pp. 341-353 ◽  
Author(s):  
Xiao‐Gui Miao ◽  
Wooil M. Moon ◽  
B. Milkereit

A multioffset, three‐component vertical seismic profiling (VSP) experiment was carried out in the Sudbury Basin, Ontario, as a part of the LITHOPROBE Sudbury Transect. The main objectives were determination of the shallow velocity structure in the middle of the Sudbury Basin, development of an effective VSP data processing flow, correlation of the VSP survey results with the surface seismic reflection data, and demonstration of the usefulness of the VSP method in a crystalline rock environment. The VSP data processing steps included rotation of the horizontal component data, traveltime inversion for velocity analysis, Radon transform for wavefield separation, and preliminary analysis of shear‐wave data. After wavefield separation, the flattened upgoing wavefields for both P‐waves and S‐waves display consistent reflection events from three depth levels. The VSP-CDP transformed section and corridor stacked section correlate well with the high‐resolution surface reflection data. In addition to obtaining realistic velocity models for both P‐ and S‐waves through least‐square inversion and synthetic seismic modeling for the Chelmsford area, the VSP experiment provided an independent estimation for the reflector dip using three component hodogram analysis, which indicates that the dip of the contact between the Chelmsford and Onwatin formations, at an approximate depth of 380 m in the Chelmsford borehole, is approximately 10.5° southeast. This study demonstrates that multioffset, three‐component VSP experiments can provide important constraints and auxiliary information for shallow crustal seismic studies in crystalline terrain. Thus, the VSP technique bridges the gap between the surface seismic‐reflection technique and well‐log surveys.


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 ◽  
1998 ◽  
Vol 63 (4) ◽  
pp. 1332-1338 ◽  
Author(s):  
Gregory S. Baker ◽  
Don W. Steeples ◽  
Matt Drake

A 300-m near‐surface seismic reflection profile was collected in southeastern Kansas to locate a fault(s) associated with a recognized stratigraphic offset on either side of a region of unexposed bedrock. A substantial increase in the S/N ratio of the final stacked section was achieved by muting all data arriving in time after the airwave. Methods of applying traditional seismic data processing techniques to near‐surface data (200 ms of data or less) often differ notably from hydrocarbon exploration‐scale processing (3–4 s of data or more). The example of noise cone muting used is contrary to normal exploration‐scale seismic data processing philosophy, which is to include all data containing signal. The noise cone mute applied to the data removed more than one‐third of the total data volume, some of which contains signal. In this case, however, the severe muting resulted in a higher S/N ratio in the final stacked section, even though some signal could be identified within the muted data. This example supports the suggestion that nontraditional techniques sometimes need to be considered when processing near‐surface seismic data.


Geophysics ◽  
2012 ◽  
Vol 77 (5) ◽  
pp. WC163-WC171 ◽  
Author(s):  
Musa S. D. Manzi ◽  
Mark A. S. Gibson ◽  
Kim A. A. Hein ◽  
Nick King ◽  
Raymond J. Durrheim

As expensive as 3D seismic reflection surveys are, their high cost is justified by improved imaging of certain ore horizons in some of the Witwatersrand basin gold mines. The merged historical 3D seismic reflection data acquired for Kloof and South Deep mines forms an integral part of their Ventersdorp Contact Reef mine planning and development programme. The recent advances in 3D seismic technology have motivated the reprocessing and reinterpretation of the old data sets using the latest algorithms, therefore significantly increasing the signal-to-noise ratio of the data. In particular, the prestack time migration technique has provided better stratigraphic and structural imaging in complex faulted areas, such as the Witwatersrand basin, relative to older poststack migration methods. Interpretation tools such as seismic attributes have been used to identify a number of subtle geologic structures that have direct impact on ore resource evaluation. Other improvements include more accurate mapping of the depths, dip, and strike of the key seismic horizons and auriferous reefs, yielding a better understanding of the interrelationship between fault activity and reef distribution, and the relative chronology of tectonic events. The 3D seismic data, when integrated with underground mapping and borehole data, provide better imaging and modeling of critical major fault systems and zones of reef loss. Many faults resolve as multifault segments that bound unmined blocks leading to the discovery and delineation of resources in faulted areas of the mines.


Geophysics ◽  
1993 ◽  
Vol 58 (3) ◽  
pp. 419-428 ◽  
Author(s):  
Arthur E. Barnes

Fourier power spectra are often usefully characterized by average measures. In reflection seismology, the important average measures are center frequency, spectral bandwidth, and dominant frequency. These quantities have definitions familiar from probability theory: center frequency is the spectral mean, spectral bandwidth is the standard deviation about that mean, and dominant frequency is the square root of the second moment, which serves as an estimate of the zero‐crossing frequency. These measures suggest counterparts defined with instantaneous power spectra in place of Fourier power spectra, so that they are instantaneous in time though they represent averages in frequency. Intuitively reasonable requirements yield specific forms for these instantaneous quantities that can be computed with familiar complex seismic trace attributes. Instantaneous center frequency is just instantaneous frequency. Instantaneous bandwidth is the absolute value of the derivative of the instantaneous amplitude divided by the instantaneous amplitude. Instantaneous dominant frequency is the square root of the sum of the squares of the instantaneous frequency and instantaneous bandwidth. Instantaneous bandwidth and dominant frequency find employment as additional complex seismic trace attributes in the detailed study of seismic data. Instantaneous bandwidth is observed to be nearly always less than instantaneous frequency; the points where it is larger may mark the onset of distinct wavelets. These attributes, together with instantaneous frequency, are perhaps, of greater use in revealing the time‐varying spectral properties of seismic data. They can help in the search for low frequency shadows or in the analysis of frequency change due to effects of data processing. Instantaneous bandwidth and dominant frequency complement instantaneous frequency and should find wide application in the analysis of seismic reflection data.


Sign in / Sign up

Export Citation Format

Share Document