Band-limited scattered wavefield reconstruction beneath complex overburdens using the Marchenko method

Author(s):  
David Vargas ◽  
Ivan Vasconcelos ◽  
Matteo Ravasi

<p>Structural imaging beneath complex overburdens, such as sub-salt or sub-basalt, typically characterized by high-impedance contrasts represents a major challenge for state-of-the-art seismic methods. Reconstructing complex geological structures in the vicinity of and below salt bodies is challenging not only due to uneven, single-sided illumination of the target area but also because of the imperfect removal of surface and internal multiples from the recorded data, as required by traditional migration algorithms. In such tectonic setups, most of the downgoing seismic wavefield is reflected toward the surface when interacting with the overburden's top layer. Similarly, the sub-salt upcoming energy is backscattered at the salt's base. Consequently, the actual energy illuminating the sub-salt reflectors, recorded at the surface, is around the noise level. In diapiric trap systems, conventional seismic extrapolation techniques do not guarantee sufficient quality to reduce exploration and production risks; likewise, seismic-based reservoir characterization and monitoring are also compromised. In this regard, accurate wavefield extrapolation techniques based on the Marchenko method may open up new ways to exploit seismic data.</p><p>The Marchenko redatuming technique retrieves reliable full-wavefield information in the presence of geologic intrusions, which can be subsequently used to produce artefact-free images by naturally including all orders of multiples present in seismic reflection data. To achieve such a goal, the method relies on the estimation of focusing operators allowing the synthesis of virtual surveys at a given depth level. Still, current Marchenko implementations do not fully incorporate available subsurface models with sharp contrasts, due to the requirements regarding the initialization of the focusing functions. Most importantly, in complex media, even a fairly accurate estimation of a direct wave as a proxy for the required initial focusing functions may not be enough to guarantee sufficiently accurate wavefield reconstruction.</p><p>In this talk, we will discuss a scattering-based Marchenko redatuming framework which improves the redatuming of seismic surface data in highly complex media when compared to other Marchenko-based schemes. This extended version is designed to accommodate for band-limited, multi-component, and possibly unevenly sampled seismic data, which contain both free-surface and internal multiples, whilst requiring minimum pre-processing steps. The performance of our scattering Marchenko method will be evaluated using a comprehensive set of numerical tests on a complex 2D subsalt model.</p>

2018 ◽  
Vol 26 (2) ◽  
pp. 229-241 ◽  
Author(s):  
Dehua Wang ◽  
Jinghuai Gao ◽  
Hongan Zhou

AbstractAcoustic impedance (AI) inversion is a desirable tool to extract rock-physical properties from recorded seismic data. It plays an important role in seismic interpretation and reservoir characterization. When one of recursive inversion schemes is employed to obtain the AI, the spatial coherency of the estimated reflectivity section may be damaged through the trace-by-trace processing. Meanwhile, the results are sensitive to noise in the data or inaccuracies in the generated reflectivity function. To overcome the above disadvantages, in this paper, we propose a data-driven inversion scheme to directly invert the AI from seismic reflection data. We first explain in principle that the anisotropic total variation (ATV) regularization is more suitable for inverting the impedance with sharp interfaces than the total variation (TV) regularization, and then establish the nonlinear objective function of the AI model by using anisotropic total variation (ATV) regularization. Next, we solve the nonlinear impedance inversion problem via the alternating split Bregman iterative algorithm. Finally, we illustrate the performance of the proposed method and its robustness to noise with synthetic and real seismic data examples by comparing with the conventional methods.


2021 ◽  
pp. 1-29
Author(s):  
Papia Nandi ◽  
Patrick Fulton ◽  
James Dale

As rising ocean temperatures can destabilize gas hydrate, identifying and characterizing large shallow hydrate bodies is increasingly important in order to understand their hazard potential. In the southwestern Gulf of Mexico, reanalysis of 3D seismic reflection data reveals evidence for the presence of six potentially large gas hydrate bodies located at shallow depths below the seafloor. We originally interpreted these bodies as salt, as they share common visual characteristics on seismic data with shallow allochthonous salt bodies, including high-impedance boundaries and homogenous interiors with very little acoustic reflectivity. However, when seismic images are constructed using acoustic velocities associated with salt, the resulting images were of poor quality containing excessive moveout in common reflection point (CRP) offset image gathers. Further investigation reveals that using lower-valued acoustic velocities results in higher quality images with little or no moveout. We believe that these lower acoustic values are representative of gas hydrate and not of salt. Directly underneath these bodies lies a zone of poor reflectivity, which is both typical and expected under hydrate. Observations of gas in a nearby well, other indicators of hydrate in the vicinity, and regional geologic context, all support the interpretation that these large bodies are composed of hydrate. The total equivalent volume of gas within these bodies is estimated to potentially be as large as 1.5 gigatons or 10.5 TCF, considering uncertainty for estimates of porosity and saturation, comparable to the entire proven natural gas reserves of Trinidad and Tobago in 2019.


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 ◽  
2000 ◽  
Vol 65 (2) ◽  
pp. 377-389 ◽  
Author(s):  
Paul J. Hatchell

Transmission distortions are observed on prestack seismic data at two locations in the Gulf of Mexico. These distortions produce anomalous amplitude versus offset (AVO) signatures. The locations of the distortion zones are determined using acquisition geometry and ray tracing. No obvious reflection events, such as shallow gas zones, are observed at the predicted locations of the distortion zones. Instead, the distortion zones correlate with buried faults and unconformities. It is postulated that the distortions are produced by velocity changes across buried faults and unconformities. The distortions result from an interference pattern resulting from seismic waves arriving from different sides of the faults. A simple model is developed to explain many of the characteristics of the distortion pattern.


2018 ◽  
Vol 6 (1) ◽  
pp. 122
Author(s):  
Okoli Austin ◽  
Onyekuru Samuel I. ◽  
Okechukwu Agbasi ◽  
Zaidoon Taha Abdulrazzaq

Considering the heterogeneity of the reservoir sands in the Niger Delta basin which are primary causes of low hydrocarbon recovery efficiency, poor sweep, early breakthrough and pockets of bypassed oil there arises a need for in-depth quantitative interpretation and more analysis to be done on seismic data to achieve a reliable reservoir characterization to improve recovery, plan future development wells within field and achieve deeper prospecting for depths not penetrated by the wells and areas far away from well locations. An effective tool towards de-risking prospects is seismic inversion which transforms a seismic reflection data to a quantitative rock-property description of a reservoir. The choice of model-based inversion in this study was due to well control, again considering the heterogeneity of the sands in the field. X-26, X-30, and X-32 were used to generate an initial impedance log which is used to update the estimated reflectivity from which we would obtain our inverted volumes. Acoustic impedance volumes were generated and observations made were consistent with depth trends established for the Niger Delta basin, inverted slices of Poisson impedances validated the expected responses considering the effect of compaction. This justifies the use of inversion method in further characterizing the plays identified in the region.


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.


2016 ◽  
Vol 4 (4) ◽  
pp. T507-T519 ◽  
Author(s):  
Yousf Abushalah ◽  
Laura Serpa

The Mamuniyat petroleum reservoir in southwestern Libya is comprised of clean sandstones and intercalated shale and sand facies that are characterized by spatial porosity variations. Seismic reflection data from the field exhibit relatively low vertical seismic resolution, side lobes of reflection wavelets, reflection interference, and low acoustic impedance contrast between the reservoir and the units underneath the reservoir, which make mapping those facies a difficult task. In the absence of broadband seismic data, optimizing frequency bands of bandlimited data can be used to suppress pseudoreflectors resulting from side-lobe effects and help to separate the clean sandstone facies of the reservoir. We have optimized the data based on our investigation of seismic frequency bands and used instantaneous frequency analysis to reveal the reflection discontinuity that is mainly associated with the reservoir boundary of the sandstone facies of the clean Mamuniyat reservoir. We also preformed rock-physics diagnostic modeling and inverted the seismic data using spectral-based colored inversion into relative acoustic impedance. The inverted impedance matches the up-scaled impedance from the well data and the inversion of relative acoustic impedance confirms the conclusion that was drawn from the instantaneous frequency results. The interpretation of facies distributions based on the instantaneous frequency was supported by the inversion results and the rock-physics models.


Geophysics ◽  
2016 ◽  
Vol 81 (3) ◽  
pp. R57-R74 ◽  
Author(s):  
Santi Kumar Ghosh ◽  
Animesh Mandal

Because seismic reflection data are band limited, acoustic impedance profiles derived from them are nonunique. The conventional inversion methods counter the nonuniqueness either by stabilizing the answer with respect to an initial model or by imposing mathematical constraints such as sparsity of the reflection coefficients. By making a nominal assumption of an earth model locally consisting of a stack of homogeneous and horizontal layers, we have formulated a set of linear equations in which the reflection coefficients are the unknowns and the recursively integrated seismic trace constitute the data. Drawing only on first principles, the Zoeppritz equation in this case, the approach makes a frontal assault on the problem of reconstructing reflection coefficients from band-limited data. The local layer-cake assumption and the strategy of seeking a singular value decomposition solution of the linear equations counter the nonuniqueness, provided that the objective is to reconstruct a smooth version of the impedance profile that includes only its crude structures. Tests on synthetic data generated from elementary models and from measured logs of acoustic impedance demonstrated the efficacy of the method, even when a significant amount of noise was added to the data. The emergence of consistent estimates of impedance, approximating the original impedance, from synthetic data generated for several frequency bands has inspired our confidence in the method. The other attractive outputs of the method are as follows: (1) an accurate estimate of the impedance mean, (2) an accurate reconstruction of the direct-current (DC) frequency of the reflectivity, and (3) an acceptable reconstruction of the broad outline of the original impedance profile. These outputs can serve as constraints for either more refined inversions or geologic interpretations. Beginning from the restriction of band-limited data, we have devised a method that neither requires a starting input model nor imposes mathematical constraints on the earth reflectivity and still yielded significant and relevant geologic information.


1991 ◽  
Vol 28 (4) ◽  
pp. 542-556 ◽  
Author(s):  
A. J. Calvert ◽  
R. M. Clowes

Multichannel deep seismic reflection data from the subduction zone of western Canada delineate the wedge of accereted sediments and the principal terranes (Crescent, Pacific Rim, and Wrangellia) that form the convergent margin. The top of the igneous oceanic crust is defined by subhorizontal reflections extending at least 100 km landward of the deformation front. Upon incorporation into the accretionary wedge, the clearly defined stratigraphy of the incoming oceanic sedimentary section is destroyed over a distance of about 10 km. Initially, an unreflective zone, which correlates well with maximum fluid expulsion, is formed. Farther landward, a predominantly landward-dipping reflectivity exists. A number of reflections are thrust faults, which appear to merge at depth with the subhorizontal reflections, but most have another origin. These reflections may be related to the movement of fluids generated by the compaction of sediments or possibly by the dehydration of the subducting plate. They are strongest in a region of depressed seismic velocities beneath the continental slope, where an analysis of reflection amplitude with offset implies that a high Poisson's ratio exists; this is consistent with the presence of elevated pore pressures. Thus, pore pressure variations associated with the migration of fluids may be the cause of much of the reflectivity within the accreted wedge, although the precipitation of minerals from rising fluids could also be important. Evidence from the seismic data also indicates that fluids from the accretionary prism are being expelled into the sediments of the overlying Tofino basin. A number of anomalously strong reflections and disruption of the horizontally stratified sediments within the lower levels of the basin probably represent fluids that migrated upward from the accreted wedge and were trapped against impermeable barriers created through the deposition of sediments on the continental slope and in the basin.


Sign in / Sign up

Export Citation Format

Share Document