Prestack Depth Migration Research Articles

Three methods of visco‐acoustic migration based on the De Wolf approximation and comparison of their migration images

AbstractThe viscosity of a medium affects the amplitude attenuation and velocity dispersion of seismic waves. Therefore, it is necessary to consider these factors during migration. First, to eliminate the viscous effect of a medium, we combine the Futterman model with the integral equation of the De Wolf approximation to construct a compensation operator of the De Wolf approximation for a visco‐acoustic medium. Next, we use the visco‐acoustic screen approximation method to realize the continuation operator then establish a prestack depth migration algorithm. Finally, an error analysis, impulse response test and model test are performed. The results show that three different generalized visco‐acoustic screen methods (phase screen method, generalized screen method and extended local Born Fourier method) can satisfactorily compensate for the attenuation of deep interface amplitude. Among these methods, the visco‐acoustic extended local Born Fourier method has the highest accuracy and the best compensation effect.

Rejuvenating exploration opportunities in Egypt's Mediterranean region using basin-scale 3D seismic reprocessing and subsurface evaluations

The Eastern Mediterranean region has proven to be a key province for natural gas resources over the last two decades. Gas discoveries in the Nile Delta Basin, Levantine Basin, and over the Eratosthenes platform have proven that an efficient working petroleum system extends to the deepwater areas and into the distal parts of the Nile Delta and Levantine basins. While the suprasalt biogenic post-Messinian play is adequately imaged in various areas in the Eastern Mediterranean region, seismic imaging of the deeper subsalt pre-Messinian plays is challenged by legacy seismic acquisition and processing techniques. Moreover, the complexity of the Messinian salt structures and the heterogeneity of the lithology of the Messinian layer across the Mediterranean basins represent a key subsurface challenge. In phase 1 of the study presented in this paper, we reprocessed 12 legacy 3D seismic data sets over 12,000 km2 and constructed a consistent basin-scale earth model from the shelf areas offshore Sinai to the deepwater areas around the Zohr carbonate platform. The data processing and imaging sequences incorporate full-waveform inversion and tomography in the earth model building process followed by a basin-to-lead evaluation to delineate exploration opportunities in phase 1. The evaluations resulted in delineating Cretaceous and Oligo-Miocene untested structures. Furthermore, potential deepwater Lower Miocene sand fairways were outlined, supported by improved regional seismic imaging and consistent seismic attributes over multiple legacy surveys. The authors of this paper expect to update the results of their work after completing the earth model building and the prestack depth migration workflows of an additional 13 legacy surveys over 35,000 km2 in phase 2 of the study.

A three-dimensional immersed boundary method for accurate simulation of acoustic wavefields with complex surface topography

Abstract Irregular topography of the free surface significantly affects seismic wavefield modelling, especially when employing finite-difference methods on rectangular grids. These methods represent the free surface as discrete points, resulting in a boundary that resembles a ‘staircase’. This approximation inaccurately represents surface topography, introducing errors in surface reflection traveltimes and generating artificial diffractions in wavefield simulation. We introduce a stable three-dimensional immersed boundary method (3DIBM) employing Cartesian coordinates to address these challenges. The 3DIBM enables the simulation of acoustic waves in media with complex topography through standard finite difference, extending the two-dimensional immersed boundary approach to compute spatial coordinates for ghost and mirror points in a three-dimensional space. Wavefield values at these points are obtained by three-dimensional spatial iterative symmetric interpolation, specifically through the Kaiser-windowed sinc method. By implicitly implementing the free surface boundary condition in three dimensions, this method effectively reduces artificial diffractions and enhances the accuracy of reflection traveltime. The effectiveness and accuracy of 3DIBM are validated through numerical tests and pre-stack depth migration imaging with simulated data, demonstrating its superiority as a modelling engine for migration imaging and waveform inversion in three-dimensional land seismic analysis.

Quantitative Elimination of Seismic Pseudofaults and Fine Analysis of True Faults Underlying Igneous Rocks of No-Well Areas: A Case Study of Shuntuoguole Uplift in Tarim Basin

Abstract After multistage tectonic movement and evolution, large superimposed oil and gas basins generally developed many igneous rocks in the early rifting stages. The lithology and lithofacies of igneous rocks are complex, which is easy to lead to the distortion of the underlying migration velocity field and thus the response of seismic pseudofaults. Also, because of the obvious shielding and absorption effect of igneous rocks on seismic waves, the waveform quality of underlying strata is poor and the seismic response characteristics of faults are fuzzy. Currently, relevant studies have shown that the influence of igneous rock can be eliminated by the prestack depth migration with an accurate igneous rock velocity model. However, improving the accuracy of the velocity model needs to be corrected by well-logging data, resulting in poor applicability of the existing velocity modeling technology underlying igneous rocks without well, which is an obvious technical bottleneck. In this paper, the secondary strike-slip fault in Shuntuoguole low uplift of Tarim Basin, which has great oil and gas exploration potential but a very low degree of drilling, is selected as the research object. Aiming at difficult fault detection underlying igneous rocks caused by lack of drilling, the accuracy of fault seismic identification is improved by “interpretative fault preprocessing” and “fault sensitive attribute optimization.” In addition, through the “extreme hypothesis method” to maximize the complex migration velocity and simulate the underlying target layer distortion maximization, we realize the quantitative elimination of seismic pseudofaults. The practical application shows that this technology can determine the true and fake underlying faults quantitatively without establishing an accurate igneous rock velocity model. It is crucial not only for exploring oil and gas in the Tarim Basin’s secondary strike-slip faults but also for offering a method and technical guide for identifying faults in other basins affected by igneous rocks.

Seismic net-to-gross estimation for a geologic model update: A case study from a turbidite lobe reservoir in the deepwater of the Niger Delta

Geologic model updates are routinely performed in mature fields to obtain improved descriptions of facies distributions and reservoir properties such as volume of shale, net to gross (NTG), and porosity, for better understanding of the static and dynamic behaviors of reservoirs for effective well placement, improved production, and monitoring. A simple but integrated seismic NTG estimation approach, using detuned seismic amplitudes is used in guiding NTG modeling during the geologic model update of a thin turbidite lobe reservoir in a mature oil field in the deep offshore Niger Delta. The objective is to address NTG overestimation and gross rock volume (GRV) uncertainty in a previous model arising from seismic amplitude tuning effects. A seismic NTG approach is chosen relative to sophisticated deterministic or stochastic inversion techniques to avoid tuning effects, which usually bias NTG estimates in thin turbidite reservoirs. The primary data set is a 2019 reprocessed prestack depth-migrated (PSDM) seismic data vintage, which had better resolution, higher signal-to-noise ratio, and more appropriate angle-stack apertures for amplitude variation with angle fidelity, than the older 2011 PSDM seismic data that are used in the previous geologic model. The methodology involved rock-physics analysis, seismic data quality checks (QCs), tuned area determination, detuning of composite seismic amplitudes of the top and base reservoir, and their direct calibration to NTG at wells. Good correlations are obtained between the detuned composite seismic amplitudes and NTG at wells. The seismic NTG map shows good calibrations at wells and provides a robust trend for net sand modeling in the oil pool and aquifer. Static model QCs and dynamic simulations prove that the seismic NTG attribute addressed the GRV uncertainty in the earlier model, thus giving confidence for using the updated model for the planning and geosteering of infill wells, sand completions, reservoir monitoring, and production.

A depth-variant wavelet extraction method based on the orthogonal matching pursuit for direct inversion

The results of pre-stack depth migration reflect subsurface structures directly. However, it is an urgent geophysical problem that how to deal with this non-stationarity exhibited by depth-domain seismic data and achieve high precision reservoir prediction. To address this issue, we develop a depth-variant wavelet extraction method based on the orthogonal matching pursuit and attain direct inversion of the depth-domain acoustic impedance using the depth-variant wavelets. First, the parameters of source wavelets estimated by complex seismic trace analysis are utilized to construct the overcomplete dictionary. Secondly, the overcomplete dictionary is used to decompose the depth-domain seismic traces by orthogonal matching pursuit. Thirdly, we extract the depth-variant Ricker wavelets from the decomposed optimal atoms. Test of examples demonstrates that the proposed depth-variant wavelet method has high accuracy. The synthetic seismogram obtained by the proposed method has high Pearson correlation coefficients of 99%, which verifies that our method is feasible and accurate. Next, a depth-domain impedance trend constraint is introduced into the objective function of the conventional basis pursuit inversion to enhance the lateral continuity. Finally, the objective function with impedance trend constraint is solved by basis pursuit for depth-domain acoustic impedance. We test the direct depth-domain acoustic impedance inversion method on the synthetic and field data, which indicates that our method has high accuracy and robustness. The mean relative error of its inversion result is only 0.9% for the noise-free example and about 5% for the strong noisy example. Therefore, our method can be effectively applied to the seismic-well calibration and seismic impedance inversion in the depth domain. And it has powerful potential in reservoir characterization, fluid prediction and attribute extraction in the depth domain. Additionally, the energy distribution features of the Ricker wavelet are derived and analyzed, which can guide the application of the Ricker wavelet in seismic signal analysis.

Tomographic imaging of South Rewa basin, Central India: Implications of Gondwana rifting and Late Cretaceous volcanism

We have imaged thick (3.0–5.0 km) sub-trappean Gondwana sediments deposited in a typical rift-environment of the south Rewa basin in Central India masked by the Deccan basalts. The basalts intrude through the basin forming several dykes, sills, and other intrusions along with large-scale undulations of the basement showing complex geology and tectonic settings due to Late Cretaceous volcanism (∼65 Ma). To delineate complex subsurface geological structures and understand the tectonic settings, we use long-offset seismic reflection data along the N-S trending 127 km Kuvari-Shahdol profile of the basin. The conventional stack and pre-stack time migration techniques fail to image the complex subsurface structures dominated by steeply dipping faults. The steeply dipping faults with small-scale structures are influenced by strong lateral and vertical velocity variations that cause problems in imaging the basin. Hence, a robust tomographic inversion with pre-stack depth migration technique is adopted to image fine-scale subtle subsurface geological structures with smooth velocity variations. The tomographic velocity model and pre-stack depth migration seismic image could decipher basaltic trap thickness, numerous dyke intrusions, alternate horsts and grabens with deposition of thick hydrocarbon-bearing sub-trappean Gondwana rocks along with the basement configuration. The deep basinal faults with highly distorted basement structures represent the intense tectonic activity of this Gondwana rift-basin. The deep-seated hydrocarbon reservoir is discovered with an excellent trapping mechanism forming both pre-rift and syn-rift settings of the rift-basin. The seismic image is further constrained and corroborated by the inversion of residual Bouguer gravity anomaly data to obtain corresponding density model derived from tomographic velocity model. The large basement faults, horsts and grabens, several intrusives like dyke, sill, laccolith, and thick sub-trappean Gondwana formations with Deccan trap cover indicate complex geology and tectonic settings of south Rewa basin forming the Late Cretaceous Volcanic Province of India.

1300-m-High gas plume from pockmarks in the north Nansha waters, South China sea

The northeastern region of Nansha Waters, located in the southern margin of the South China Sea, is abundant in gas resources. However, it remains unclear if the central and northern areas of Nansha Waters possess gas resources and if there are any associated leaks. The Jiuzhang Basin is located in the northern Nansha Waters. Pre-stack depth migration, analysis of seismic attributes and multibeam data processing were conducted to image the seabed topography and backscatter intensity, as well as the seismic structure of the water and sedimentary layers in the Jiuzhang Basin. The result reveals a presence of large plumes, approximately 1500 m in width and up to 1300 m in height, from a pockmark. The plume is characterized by a strong seismic amplitude featuring chaotic reflections, a high instantaneous frequency (100–150 Hz), and a strong instantaneous amplitude. These seismic features commonly appear in hydrocarbon-bearing seawater, indicating that the plume is composed of gas. The relationship between the volume fraction of bubbles in the gas plumes and the amplitude of multi-channel seismic waves suggests that the volume fraction of bubbles within the gas plume can reach up to 5.0 × 10−3. Meanwhile, a new pockmark with a length of 248 m, a width of 90 m, and a depth of 20.4 m, and an abrupt 22 m drop in water depth of an existing pockmark are observed in the gas plume area. These observations, combined with the gas plume revealed by the multi-channel seismic collected in 2013 and multibeam collected in 2020, suggest a very strong and recently active gas leaks in the pockmark area of the northern Jiuzhang Basin. Furthermore, in the pockmark area, there are gas-bearing formations with weak reflections, and mud diapirs charactered by low seismic velocity, low amplitude chaotic reflections, and pulled up reflections on both sides, and faults extend to the seabed. These observations suggest that gas plume originate from deep strata of the northern Jiuzhang Basin. Large gas plumes and mud diapirs are often indicative of the existence of deep gas reservoirs. Therefore, the presence of such active gas plumes and mud diapirs in the northern region of the Jiuzhang Basin signifies the presence of abundant gas resources.

A depth-variant seismic wavelet extraction method for basis pursuit inversion with an impedance trend constraint

The seismic images produced by prestack depth migration indicate more accurate subsurface structures than time images, resulting in a growing need for depth-domain inversion. However, due to the strong nonstationarity exhibited by depth-domain seismic data, time-domain inversion methods based on the convolutional model cannot be directly applied in the depth domain. To address this issue, we develop a method for extracting a depth-variant seismic wavelet, which is then combined with a nonstationary convolutional model to enable direct inversion of the depth-domain acoustic impedance (AI). First, we extend the Morlet wavelet to the depth domain and develop an orthogonal matching pursuit spectral decomposition method using the depth-domain Morlet wavelet. We then investigate the waveforms and wavenumber spectra similarities between the depth-domain Morlet wavelet and depth-domain Ricker wavelet and extract depth-variant Ricker wavelets from the depth-wavenumber spectrum. We add a depth-domain impedance trend constraint to the conventional basis pursuit inversion to enhance the lateral continuity of the inversion results. Then, we attain direct inversion of the depth-domain AI. Tests of synthetic and field data demonstrate that our method achieves high-accuracy inversion results while maintaining high computational efficiency, highlighting our approach’s effectiveness and strong reservoir characterization potential.

A Wavelet Extraction Method of Attenuation Media for Direct Acoustic Impedance Inversion in Depth Domain

The seismic image produced by pre-stack depth migration is more accurate and has clearer geological significance than the time image. However, the waveform of the depth-domain seismic image is affected not only by depth-dependent velocity variation but also by media attenuation, resulting in strong spectral variation of depth-domain seismic data. Therefore, depth-domain seismic inversion is still challenging. We propose a wavelet extraction method of attenuation media based on the generalized seismic wavelet, to address this issue. Then, the estimated depth-domain wavelets were applied to the direct acoustic impedance inversion. First, we investigated the effect of attenuation media on depth-domain source wavelets and derived an analytical formula for the depth-domain wavelets of attenuation media. Next, the time-domain generalized seismic wavelet was extended to the depth domain, which was utilized to study the feasibility of using the generalized seismic wavelet to characterize the seismic wavelet of the depth-domain attenuation media. Based on the orthogonal matching pursuit, we propose a method to extract the depth-domain generalized seismic wavelet directly from depth-domain seismic data. Finally, we applied this method to the depth-domain direct acoustic impedance inversion of a 3D field data example. Tests on the synthetic and 3D field datasets show that the proposed method can correctly extract the depth-domain seismic wavelet of attenuation media and attain direct inversion of the depth-domain acoustic impedance with high accuracy. Therefore, our method is effective and has robust potential in reservoir characterization, fluid prediction, and attribute extraction in the depth domain.

Imaging the shallow velocity structure of the slow‐spreading ridge of the South China Sea with downward continued multichannel seismic data

AbstractHigh‐resolution shallow oceanic crust velocity models provide crucial information on the tectonothermal history of the oceanic crust. The ocean bottom seismometers record wide‐angle seismic reflection and refraction data to image deeper structures compared with streamer data set. However, most ocean bottom seismometers experiments produce low‐resolution velocity models with limited shallow crustal structure due to sparse ocean bottom seismometers spacing. Multichannel seismic data recorded by towed streamers provide complementary seismic images of the oceanic crust but yield little information on subseafloor velocity because most subseafloor refractions are masked by seafloor reflections. Therefore, it is difficult to obtain fine‐scale velocity structure of shallow upper oceanic crust with both ocean bottom seismometers and multichannel seismic data. Downward continuation technique redatumed the shots and receivers to the seafloor to collapse the seafloor reflections to the zero offset and extract refractions as first arrivals from nearly zero offset, enabling dense ray coverage at the shallow crust. We applied the downward continuation and traveltime tomography methods to two synthetic models, Marmousi and SEAM Phase I Salt models, to demonstrate the performance of the strategy in the situations of flat seafloor and rough seafloor topography. We conducted the first‐arrival traveltime tomography on downward continued towed‐streamer multichannel seismic data across a slow‐spreading ridge of the South China Sea, providing unprecedented details of shallow velocity structure in the sediments. The low velocity sediments revealed by traveltime tomography match well with the prestack depth migration profile.

Seismic Reprocessing Overcomes Uncertainty in Offshore Egyptian Field

_ This article, written by JPT Technology Editor Chris Carpenter, contains highlights of paper SPE 214066, “Seismic Reprocessing Leads to New Breakthroughs—A Successful Case in the ASH Field, AG Basin in Egypt,” by Mohamed Elokr and Ahmed Lotfy, Kuwait Energy Egypt, and Wei Xing, United Energy Group, et al. The paper has not been peer reviewed. _ The ASH oil field is in the eastern portion of the AG Basin in Egypt. The Lower Cretaceous Alam El Bueib is the main oil-producing formation. Because of the high heterogeneity of the Abu Roash succession, in addition to the influence of thick limestone of the Upper Cretaceous and the influence of multiple complex faults, the quality of seismic data is very poor, requiring seismic reprocessing. Two key techniques were used to achieve the goal: New vertical seismic profile (VSP) well data were acquired to adjust the velocity model, and common reflection angle migration (CRAM) prestack depth migration (PSDM) was used for reprocessing. Introduction The ASH field produces oil from the Aptian Lower Cretaceous reservoir. The ASH structure is interpreted as an elongated east-northeast/west-southwest anticline with intensive breaching on the downthrown side of northwest/southeast to west-northwest/east-southeast-oriented normal faults with obvious signature of synsedimentary activities. The entire area was subject to a series of tectonic events affecting the Western Desert, which led to a complex fault regime clearly recognized in the ASH structure by a series of horizontal displacement across the northwest/southeast faults. Strike/slip to oblique-slip faults also are observed from thicknesses of the cretaceous formations, creating significant difficulty in fault definition. The 3D seismic data acquired in 2007 suffered from losing amplitude of the Lower Cretaceous and deep Jurassic reflectors because of the complex structure. The strong carbonate markers disappear north of the field, where high, dense fault intersection is present closing from the main bounding fault. Loss of seismic amplitude could have severe effects on of the future economics of the field. The seismic survey was reprocessed two times without significant enhancement of the fault imaging. PSDM was run without velocity control, resulting in an uncertain velocity model with inaccurate fault imaging and a mismatching between actual and prognosed depths of the drilled wells of greater than 50 m. This contradiction in reservoir mapping meant that the northern portion of the field remained unexplored. Acquiring VSP in the recently drilled well encouraged reprocessing of the 3D seismic survey using well velocity as a control point in adjusting the seismic velocity model, and application of CRAM technology to improve fault imaging and appraise the uncertain area of the ASH field. CRAM Methodology CRAM was developed as a type of beam migration. Asymptotic ray tracing is performed, assuming that one-way diffraction rays from the subsurface image point to the surface sources and receivers. The takeoff angles from the image point are measured around a given local normal to a background reflection surface for each source ray and receiver ray. A system of source and receiver ray pairs is formed, and the recorded seismic data can be mapped into a local angle domain based on the reflection angles at the subsurface image points.​​​​

Lateral shift from turbidite- to contourite-dominated continental slope, a case study from southeast Brazil slope

New bathymetric and subsurface data from the Santos Basin, offshore southeast Brazil, reveal sharp lateral morphological variations facing the margin-parallel (along-margin) flow of the Antarctic Intermediate Water. The data, derived from 3D Pre-Stack Depth Migrated (PSDM) seismic reflection volumes, show on a ~100 km long Southern segment slope gullies, mass transport deposits, and scars. These morphologies sharply shift laterally to a ~50 km long Northern segment, typified by a smooth relief built of plastered contourite drifts with occasional elongated depressions. We examine this sharp lateral shift in the morphology of bottom current indicators in view of the hydrodynamical regime, sediment sources, and their distribution pattern. Results indicate that the sharp lateral morphological shift is attributed to (1) a 30° change in slope orientation, and (2) a northward increase in the across-margin velocity component of the Intermediate Water Boundary Current (IWBC). Subsurface data suggest this flow pattern has persisted since the mid-Miocene. In this case study, the lateral morphological shifts under contemporary hydrodynamic conditions emphasize the importance of understanding the coupling between ocean bottom currents and the seafloor in any marine environment.

Prestack Depth Migration supported by Madagascar software - Summary of various methods through application to Sigsbee2a model -

Prestack Depth Migration supported by Madagascar software - Summary of various methods through application to Sigsbee2a model -

Improving prestack time migration by introducing a new velocity-related parameter: Parameter picking and 3D real data application

Prestack time migration (PSTM), a commonly used tool for seismic imaging, has been widely applied in 3D seismic data processing. However, the conventional PSTM algorithms use only one effective velocity parameter (i.e., root-mean-square [rms] velocity) for each imaging point, which may not be accurate when stronger lateral variations occur in seismic velocities. In this paper, we introduce a new parameter called the velocity variation factor that considers velocity variations in inhomogeneous media to improve PSTM. This new parameter, together with the rms velocity, describes the propagation Green’s function at an imaging point with two effective parameters rather than one effective parameter as in conventional PSTMs. This provides a more accurate traveltime calculation for the wave propagating through media with moderate lateral velocity variation. Unlike the conventional bending-ray PSTM, the additional effective parameter is fully independent of the rms velocities. We estimate the two effective parameters at each imaging point by flattening the neighboring image gathers with a global optimization algorithm. The objective function is built at each imaging point using a selective crosscorrelation-based time shift, which can quantitatively describe the slight bending of events in the local migrated gathers regardless of the quality of the gathers. We estimate the two effective parameters using the very fast simulated annealing algorithm and multiscale approach, thus avoiding the local minimum caused by the noises in the migrated gathers. We apply our two-parameter PSTM to a real 3D land data set to demonstrate its industrial applicability. A comparison of the new imaging result with the conventional prestack depth migration also is presented.

IMPROVED SEISMIC IMAGING THROUGH PETROPHYSICAL APPROACH USING THOMSEN'S PARAMETER ESTIMATION IN ANISOTROPIC PRESTACK DEPTH MIGRATION (PSDM)

Research has been carried out using the Prestack Depth Migration (PSDM) anisotropy method for 2D reflection seismic data. This research aims to estimate anisotropic parameter values with a petrophysical approach to the Prestack Depth Migration method to obtain better imaging results. The type of anisotropic medium used is Vertical Transverse Isotropy (VTI) because that can explain the effect of anisotropy in a simple form on the sediment layer. In theory, 2 parameters are needed to describe this parameter, namely ε and δ. δ is an anisotropic parameter that describes the velocity variation towards a nearly vertical direction and is the control depth of seismic. In contrast, ε describes the velocity variation towards the near horizontal direction. The stages of data processing are divided into two; the first is PSDM isotropy, which flattens gather at near offset (angle mute <30 degrees), and PSDM anisotropy to flatten-gather at far offset (mute angle> 30 degrees). The results of this study showed that the relation of ε and δ in shale formations was (δ = 0.4958ε - 0.0152) and sand formation (δ = 0.9082ε - 0.0203) and the range of anisotropic parameters δ and ε were -0.02 up to 0.13 and 0 up to 0.3. From all ranges of anisotropic parameter values, this study belongs to the "weak anisotropy" category Thomsen (0 to 0.5). The value of anisotropy parameters obtained can give the results of seismic images more accurately and clearly than PSDM isotropy. It can correct errors in the depth of seismic up to 51 m.

The lithosphere–asthenosphere boundary structure at 11–21 Ma from wide-angle seismic data in the equatorial Atlantic Ocean

SUMMARYThe lithosphere–asthenosphere boundary (LAB) separates the rigid lithospheric plate above with the ductile and convective asthenosphere below and plays a fundamental role in plate tectonic processes. The LAB has been imaged using passive geophysical methods, but these methods only provide low-resolution images. Recently, seismic reflection imaging method has provided high-resolution images of the LAB, but imaging of the LAB at younger ages has been difficult. Here, we present the image of the LAB using wide-angle seismic reflection data covering 11–21 Ma old lithosphere in the equatorial Atlantic Ocean. Using ocean bottom seismometers (OBSs), we have observed wide-angle reflections between 150 and 400 km offsets along with crustal and mantle refraction arrivals. We first performed traveltime tomography to obtain the velocity in the crust and upper mantle. The Pn arrivals provide the information about P-wave velocity down to 4 km below the Moho. The disappearance of Pn arrivals beyond 130 km offset suggests that vertical P-wave velocity gradient is negligible or negative below this depth. We extended these velocities down to 90 km depth and then applied two imaging techniques to wide-angle reflection data, namely traveltime mapping of picked reflection arrivals and pre-stack depth migration of full wavefield data. We find that these reflections originate between 34 and 67 km depth, possibly from the LAB system. We have carried out extensive modelling to show that these reflections are real and not artefacts of imaging. Comparison of our results with coincident passive seismological and magnetotelluric results suggests that wide-angle imaging technique can be successfully used to study the lithosphere and the LAB system. We find that the LAB gradually deepens with age, but becomes very deep at 17–19 Ma, which we interpret to be due to the anomalous geology along this part of the profile.

ISTRICI – Tools for facilitating seismic depth imaging and velocity analysis with seismic unix

ISTRICI was developed to facilitate the use of Seismic Unix, a free software designed to process seismic data. Seismic Unix is a powerful tool for performing pre-stack depth migration and using the residual parameters to update the seismic velocity model based on the pre-stack depth migration results. What is missing are codes and scripts that allow practical, interactive, and easy application of the Seismic Unix algorithms. Therefore, ISTRICI was developed to fill this gap. ISTRICI consists of three workflows (INTER, CIG and TRAD) for interactively performing residual velocity analysis and pre-stack depth migration, switching from one workflow to the next depending on the characteristics and type of seismic data. To demonstrate the flexibility and efficiency of our package, we describe the application of this tool on marine seismic data.

Wide‐spectrum reconstruction of a velocity model based on the wave equation reflection inversion method and its application

AbstractMost seismic data used for conventional exploration lack far‐offset signals and effective low‐frequency components. This makes it difficult for conventional full waveform inversion methods to recover the low‐wavenumber components of mid‐ and deep‐layer models. To resolve the bottleneck of ray‐theoretical reflection traveltime tomography, wave equations, such as reflection inversion methods, are receiving increasing attention. We reconstructed a wide‐spectrum velocity model based on the multiscale wave‐equation reflection inversion method for model decomposition. First, using the dynamic image warping method to obtain reflection traveltime residuals, the wave‐equation reflection traveltime inversion was used to recover the low‐wavenumber components of the background model, which also solved the cycle‐skipping problem. Second, the medium‐wavenumber component of the model is then supplemented by the reflection waveform inversion, while the reflectivity model is obtained relying on the least‐squares reverse time migration method. Based on this method, the background and perturbation models are updated alternately and iteratively. At the same time, the stratum structural tensor information was extracted using the perturbed model image to construct a preconditioned operator for the stratum structural constraint, suppress unreasonably high‐wavenumber components in the generalized gradient and improve the geological consistency of the inversion results. Testing of the model using the Sigsbee2b model and seismic field data from the East China Sea showed that, compared with the conventional reflection traveltime tomography method based on prestack depth migration, the wave‐equation reflection inversion strategy with well‐matched information from traveltimes to waveforms significantly improved the accuracy of the middle and deep velocity modelling and enhanced the imaging quality of the whole body.

Evaluation of the Possibilities Validation of Interval Velocity Models Using Non-Seismic Data and Its Impact on Geological Interpretation of PreSDM Results

The paper presents the problem of generation and validation of Velocity Interval Depth (VID) models with the application of non-seismic geophysical and geological data. The study area is a part of the Carpathian Foredeep located close to its contact with the Carpathian Overthrust. In this area of complicated geological structure, hydrocarbon deposits have been successfully explored for decades with seismic methods and drilling. The research applied the Simultaneous Joint Inversion (SJI) of independent geophysical data, which is a modern methodology of geophysical data processing, that is still under development. Such an attempt was necessary due to the lack of a sufficiently dense grid of wells in the study area, in which seismic velocities would be correctly recorded. Such data would be then applied for the generation of relevant VID models, which in turn, could be used to perform the Prestack Depth Migration (PreSDM) procedures. The application of procedures taking advantage of independent geophysical and geological data enabled researchers to control the generation process of the spatial VID model in the areas without wells. The analyses aimed to verify the correctness of VID model evaluation and its influence on the quality of seismic imaging in the area of the Carpathian Overthrust. Precisely, the influence was tested of such non-standard generation procedure of seismic velocity fields, not only on the PreSDM results but also on the geological interpretation of both the Rączyna and the Jodłówka gas deposits. The latter aspect of the presented results seems to be crucial to the effectiveness of petroleum exploration in the transition zone between the Carpathian Orogen and the Carpathian Foredeep.